Pub Date : 2024-09-21DOI: 10.1016/j.epsl.2024.119010
Samuel Ebert , Kazuhide Nagashima , Alexander N. Krot , Shigeru Wakita , Jean-Alix Barrat , Addi Bischoff
Calcium-aluminum-rich inclusions (CAIs) commonly observed in chondritic meteorites are the oldest dated solids formed in the Solar System. Short-lived isotope chronologies (26Al-26Mg, 182Hf-182W) suggest a ∼2 Ma gap between the formation of CAIs and the accretion of the final chondrite parent bodies. One thin section, 3.27 cm2 in size, of an ordinary chondrite NWA 3358 (H3.1) studied contains 52 refractory inclusions (CAIs and amoeboid olivine aggregates (AOAs)) comprising 0.14 % of its area, which is the highest abundance of refractory inclusions among non-carbonaceous chondrites containing on average ∼0.009 area % of CAIs and AOAs. In combination with a low chondrule/matrix ratio of ∼1.5, this makes NWA 3358 a unique ordinary chondrite. The aqueously-formed fayalites (Fa>99) in NWA 3358 have the inferred initial 53Mn/55Mn ratio of (5.56 ± 0.44) × 10−6 which is the highest measured value for secondary minerals in chondrites and corresponds to the formation time of ∼1.0–1.5 Ma after CAIs. Based on the 53Mn-53Cr chronology of fayalite formation and the thermal modeling, we infer that the first-generation of an H chondrite parent body, ∼6–12 km in diameter, accreted within 1.0 Ma after formation of CAIs, filling the gap of ∼2 Ma between CAIs and the earliest chondrite parent bodies. This early accretion provides a possible mechanism of CAIs/AOAs storage in the inner solar nebula and could explain the high amount of refractory inclusions in NWA 3358. A later destruction of these first-generation bodies may also explain the presence of CAIs and chondrules of different ages within later formed chondrite parent bodies.
{"title":"Early generation of a refractory inclusions-enriched H-chondritic parent body: A safe harbor for Ca, Al-rich inclusions","authors":"Samuel Ebert , Kazuhide Nagashima , Alexander N. Krot , Shigeru Wakita , Jean-Alix Barrat , Addi Bischoff","doi":"10.1016/j.epsl.2024.119010","DOIUrl":"10.1016/j.epsl.2024.119010","url":null,"abstract":"<div><p>Calcium-aluminum-rich inclusions (CAIs) commonly observed in chondritic meteorites are the oldest dated solids formed in the Solar System. Short-lived isotope chronologies (<sup>26</sup>Al-<sup>26</sup>Mg, <sup>182</sup>Hf-<sup>182</sup>W) suggest a ∼2 Ma gap between the formation of CAIs and the accretion of the final chondrite parent bodies. One thin section, 3.27 cm<sup>2</sup> in size, of an ordinary chondrite NWA 3358 (H3.1) studied contains 52 refractory inclusions (CAIs and amoeboid olivine aggregates (AOAs)) comprising 0.14 % of its area, which is the highest abundance of refractory inclusions among non-carbonaceous chondrites containing on average ∼0.009 area % of CAIs and AOAs. In combination with a low chondrule/matrix ratio of ∼1.5, this makes NWA 3358 a unique ordinary chondrite. The aqueously-formed fayalites (Fa<sub>>99</sub>) in NWA 3358 have the inferred initial <sup>53</sup>Mn/<sup>55</sup>Mn ratio of (5.56 ± 0.44) × 10<sup>−6</sup> which is the highest measured value for secondary minerals in chondrites and corresponds to the formation time of ∼1.0–1.5 Ma after CAIs. Based on the <sup>53</sup>Mn-<sup>53</sup>Cr chronology of fayalite formation and the thermal modeling, we infer that the first-generation of an H chondrite parent body, ∼6–12 km in diameter, accreted within 1.0 Ma after formation of CAIs, filling the gap of ∼2 Ma between CAIs and the earliest chondrite parent bodies. This early accretion provides a possible mechanism of CAIs/AOAs storage in the inner solar nebula and could explain the high amount of refractory inclusions in NWA 3358. A later destruction of these first-generation bodies may also explain the presence of CAIs and chondrules of different ages within later formed chondrite parent bodies.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119010"},"PeriodicalIF":4.8,"publicationDate":"2024-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004424/pdfft?md5=8b09b36f84f5eb27d51cdd486c15c50e&pid=1-s2.0-S0012821X24004424-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142272575","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-21DOI: 10.1016/j.epsl.2024.119016
Heng-Ci Tian , Fang-Zhen Teng , Xin-Yang Chen , Ilya N. Bindeman , Jeffrey G. Ryan
Arc lavas display significant chemical and isotopic heterogeneity mainly due to recycled materials from subducting slabs. However, the extent to which different types of subducted sediments and oceanic crust contribute to the petrogenesis of arc magmas, as well as the roles of the mantle wedge and overlying crust, remain debated. Potassium (K) isotopes have the potential to provide new insights into the processes and sources of arc magmatism because sediments and altered oceanic crust are highly enriched in K and have distinct δ41K values compared with the mid-ocean ridge basalts and upper mantle (-0.42 ± 0.08‰, 2SD). Here we report K isotopic compositions of 32 well-characterized arc lavas from the circum-Pacific margins. We find low δ41K values (-0.86‰ ∼ -0.38‰) in the Setouchi arc samples, which we interpret as the result of incorporation of isotopically light sediments into the subarc mantle. The Kurile and Panama arc lavas have high δ41K values (-0.36‰ ∼ 0.02‰) and their δ41K values correlate positively with the Ba/Th ratios, indicating ∼0.5–2% fluid additions from dehydrated altered crust. Adakites have variable but overall heavy K isotope compositions (-0.44‰ ∼ -0.01‰) , which correlate with the diagnostic trace-element indicators of Sr/Y, La/Yb and K/U. These results indicate limited K isotope fractionation during metamorphic dehydration of hot slab subduction, and/or fluid metasomatism before the slab melting. Our study suggests that fluids and melts originating from subducting sediments and slabs exerted distinct influences on the origins of different types of arc magmas, demonstrating the significance of K isotopes in studying the petrogenesis of arc magmas.
{"title":"Tracing island arc petrogenesis using potassium isotopes","authors":"Heng-Ci Tian , Fang-Zhen Teng , Xin-Yang Chen , Ilya N. Bindeman , Jeffrey G. Ryan","doi":"10.1016/j.epsl.2024.119016","DOIUrl":"10.1016/j.epsl.2024.119016","url":null,"abstract":"<div><div>Arc lavas display significant chemical and isotopic heterogeneity mainly due to recycled materials from subducting slabs. However, the extent to which different types of subducted sediments and oceanic crust contribute to the petrogenesis of arc magmas, as well as the roles of the mantle wedge and overlying crust, remain debated. Potassium (K) isotopes have the potential to provide new insights into the processes and sources of arc magmatism because sediments and altered oceanic crust are highly enriched in K and have distinct δ<sup>41</sup>K values compared with the mid-ocean ridge basalts and upper mantle (-0.42 ± 0.08‰, 2SD). Here we report K isotopic compositions of 32 well-characterized arc lavas from the circum-Pacific margins. We find low δ<sup>41</sup>K values (-0.86‰ ∼ -0.38‰) in the Setouchi arc samples, which we interpret as the result of incorporation of isotopically light sediments into the subarc mantle. The Kurile and Panama arc lavas have high δ<sup>41</sup>K values (-0.36‰ ∼ 0.02‰) and their δ<sup>41</sup>K values correlate positively with the Ba/Th ratios, indicating ∼0.5–2% fluid additions from dehydrated altered crust. Adakites have variable but overall heavy K isotope compositions (-0.44‰ ∼ -0.01‰) , which correlate with the diagnostic trace-element indicators of Sr/Y, La/Yb and K/U. These results indicate limited K isotope fractionation during metamorphic dehydration of hot slab subduction, and/or fluid metasomatism before the slab melting. Our study suggests that fluids and melts originating from subducting sediments and slabs exerted distinct influences on the origins of different types of arc magmas, demonstrating the significance of K isotopes in studying the petrogenesis of arc magmas.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119016"},"PeriodicalIF":4.8,"publicationDate":"2024-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004485/pdfft?md5=da4db8d8b9d80e12ff368d07083c6390&pid=1-s2.0-S0012821X24004485-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142312112","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-20DOI: 10.1016/j.epsl.2024.118973
Loraine Gourbet , Sean F. Gallen , Vincent Famin , Laurent Michon , Miangaly Olivia Ramanitra , Eric Gayer
Climate's role in governing landscape evolution has been intensely studied for several decades, but few studies clearly document climate-landscape interactions in natural landscapes. This study aims to improve understanding of climate-landscape linkages using hotspot volcanic islands in the tropics as natural laboratories. Relatively uniform lithology, strong precipitation and climate gradients, and known initial topographic conditions on Réunion and Mauritius islands (Réunion hotspot) and Kaua'i (Hawaii hotspot) enable us to explore the impact of climate on erosion rates and geomorphic process. We reconstruct paleo-topography of drainage basins based on preserved remnants of relict topography from past volcanic events that repaved the landscapes that are differenced from the modern-day topography to determine eroded volumes. Existing geochronology of the volcanic flows allows us to constrain the timing of repaving (a proxy for the initiation of erosion) and basin average erosion rates. The initial and final conditions and the duration of erosion are used to calibrate a simple stream power model for bedrock river incision for each basin using a Bayesian inversion. We compare the erosion rate and calibrated stream power parameters to precipitation and climate data for each drainage basin on each island to explore potential relationships. Results show that basin average erosion rates for basins eroding < ∼1 mm/yr show a positive relationship with mean annual precipitation (MAP) and a negative relationship with the duration of erosion. Importantly, MAP and erosion duration are correlated, so we infer that the negative correlation between erosion rate and duration is coincidental. The stream power slope exponent and erodibility coefficient only exhibit significant correlations with climate parameters for Réunion Island, particularly mean annual cyclonic precipitation. Our results demonstrate that both mean annual precipitation and extreme events control long-term landscape evolution on volcanic islands.
{"title":"Climate impact on river incision on hotspot volcanoes: insights from paleotopographic reconstructions and numerical modelling","authors":"Loraine Gourbet , Sean F. Gallen , Vincent Famin , Laurent Michon , Miangaly Olivia Ramanitra , Eric Gayer","doi":"10.1016/j.epsl.2024.118973","DOIUrl":"10.1016/j.epsl.2024.118973","url":null,"abstract":"<div><p>Climate's role in governing landscape evolution has been intensely studied for several decades, but few studies clearly document climate-landscape interactions in natural landscapes. This study aims to improve understanding of climate-landscape linkages using hotspot volcanic islands in the tropics as natural laboratories. Relatively uniform lithology, strong precipitation and climate gradients, and known initial topographic conditions on Réunion and Mauritius islands (Réunion hotspot) and Kaua'i (Hawaii hotspot) enable us to explore the impact of climate on erosion rates and geomorphic process. We reconstruct paleo-topography of drainage basins based on preserved remnants of relict topography from past volcanic events that repaved the landscapes that are differenced from the modern-day topography to determine eroded volumes. Existing geochronology of the volcanic flows allows us to constrain the timing of repaving (a proxy for the initiation of erosion) and basin average erosion rates. The initial and final conditions and the duration of erosion are used to calibrate a simple stream power model for bedrock river incision for each basin using a Bayesian inversion. We compare the erosion rate and calibrated stream power parameters to precipitation and climate data for each drainage basin on each island to explore potential relationships. Results show that basin average erosion rates for basins eroding < ∼1 mm/yr show a positive relationship with mean annual precipitation (MAP) and a negative relationship with the duration of erosion. Importantly, MAP and erosion duration are correlated, so we infer that the negative correlation between erosion rate and duration is coincidental. The stream power slope exponent and erodibility coefficient only exhibit significant correlations with climate parameters for Réunion Island, particularly mean annual cyclonic precipitation. Our results demonstrate that both mean annual precipitation and extreme events control long-term landscape evolution on volcanic islands.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 118973"},"PeriodicalIF":4.8,"publicationDate":"2024-09-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004060/pdfft?md5=ab2e278bef72957cd077b05c8a12cf0e&pid=1-s2.0-S0012821X24004060-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142272956","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-19DOI: 10.1016/j.epsl.2024.119005
P. Baccheschi , J.M. Confal , S. Pondrelli
The region between central Europe and the centre of the Mediterranean is characterised by complex tectonics and kinematics. Here, the interaction between thickened crust, subducting lithosphere and surrounding asthenosphere produces strong and pervasive anisotropy in the upper mantle. Shear wave splitting measurements, the most adopted method to image seismic anisotropy so far, when interpreted in a ray-based framework result in little or no depth resolution, hampering a correct image of the anisotropy distribution with depth. In this study, we aim to better constrain the depth-dependent seismic anisotropy beneath Italy and surrounding regions, by isolating for the first time the source region of anisotropy at different depths. To do that, we perform an anisotropy tomography, adopting the splitting intensity inversion method. It is entirely based on the finite-frequency effect in the splitting of SKS waves. We first computed the splitting intensity using SKS waves recorded at all available permanent and temporary stations over the region, obtaining a huge dataset of measurements used as an input for the tomographic inversion. The large-scale 3D model of seismic anisotropy obtained with the inversion shows a clear change of anisotropy properties in terms of fast polarisation direction and intensity for different depths, thus improving the characterization of the main sources of anisotropy in the mantle as a function of depth. Shallower layers (70–100 km depth) are characterised by a complex and variable oriented pattern of anisotropy fast direction and intensity, which becomes progressively more organised with depth (100–300 km). This pattern suggests a strong control exerted by the geometry and motion of the different slab segments and the large-scale asthenospheric flow generated by subduction and roll-back processes. The strength of anisotropy increases with depth, with high values affecting the bulge of the Alps and Apennines chains and the southern Tyrrhenian subduction system. On the contrary, weaker anisotropy characterises the transition zone from the Apennines to Alpine domains beneath the Po plain, and both the Adriatic and European domains.
欧洲中部和地中海中部之间的地区具有复杂的构造和运动学特征。在这里,增厚的地壳、俯冲岩石圈和周围的岩浆层之间的相互作用在上地幔中产生了强烈而普遍的各向异性。剪切波分裂测量是迄今为止最常用的地震各向异性成像方法,但在基于射线的框架下解释时,深度分辨率很低或根本没有,从而阻碍了对各向异性随深度分布的正确成像。在本研究中,我们旨在通过首次分离不同深度的各向异性源区,更好地约束意大利及周边地区地下随深度变化的地震各向异性。为此,我们采用分裂强度反演法进行了各向异性层析成像。该方法完全基于 SKS 波分裂中的有限频率效应。我们首先利用该地区所有可用的永久和临时台站记录的 SKS 波计算分裂强度,获得了一个庞大的测量数据集,作为层析反演的输入。通过反演获得的地震各向异性大尺度三维模型显示,不同深度的各向异性在快速极化方向和强度方面发生了明显变化,从而改善了地幔各向异性主要来源随深度变化的特征。较浅地层(70-100 千米深度)的特征是各向异性快速方向和强度的复杂多变的定向模式,这种模式随着深度(100-300 千米)的增加而逐渐变得更有组织。这种模式表明,不同板块的几何形状和运动以及俯冲和滚回过程产生的大尺度星震层流具有很强的控制力。各向异性的强度随深度的增加而增加,高值影响到阿尔卑斯山脉和亚平宁山脉的隆起以及南第勒尼安海俯冲系统。相反,在波河平原下从亚平宁山脉向阿尔卑斯山脉的过渡区以及亚得里亚海和欧洲海域,各向异性较弱。
{"title":"Splitting intensity tomography to image depth-dependent seismic anisotropy patterns beneath the Italian Peninsula and surrounding regions","authors":"P. Baccheschi , J.M. Confal , S. Pondrelli","doi":"10.1016/j.epsl.2024.119005","DOIUrl":"10.1016/j.epsl.2024.119005","url":null,"abstract":"<div><p>The region between central Europe and the centre of the Mediterranean is characterised by complex tectonics and kinematics. Here, the interaction between thickened crust, subducting lithosphere and surrounding asthenosphere produces strong and pervasive anisotropy in the upper mantle. Shear wave splitting measurements, the most adopted method to image seismic anisotropy so far, when interpreted in a ray-based framework result in little or no depth resolution, hampering a correct image of the anisotropy distribution with depth. In this study, we aim to better constrain the depth-dependent seismic anisotropy beneath Italy and surrounding regions, by isolating for the first time the source region of anisotropy at different depths. To do that, we perform an anisotropy tomography, adopting the splitting intensity inversion method. It is entirely based on the finite-frequency effect in the splitting of SKS waves. We first computed the splitting intensity using SKS waves recorded at all available permanent and temporary stations over the region, obtaining a huge dataset of measurements used as an input for the tomographic inversion. The large-scale 3D model of seismic anisotropy obtained with the inversion shows a clear change of anisotropy properties in terms of fast polarisation direction and intensity for different depths, thus improving the characterization of the main sources of anisotropy in the mantle as a function of depth. Shallower layers (70–100 km depth) are characterised by a complex and variable oriented pattern of anisotropy fast direction and intensity, which becomes progressively more organised with depth (100–300 km). This pattern suggests a strong control exerted by the geometry and motion of the different slab segments and the large-scale asthenospheric flow generated by subduction and roll-back processes. The strength of anisotropy increases with depth, with high values affecting the bulge of the Alps and Apennines chains and the southern Tyrrhenian subduction system. On the contrary, weaker anisotropy characterises the transition zone from the Apennines to Alpine domains beneath the Po plain, and both the Adriatic and European domains.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119005"},"PeriodicalIF":4.8,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004370/pdfft?md5=c4bfcdd6b9836484625a7993e7c67092&pid=1-s2.0-S0012821X24004370-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142242826","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-19DOI: 10.1016/j.epsl.2024.119004
Ping Wang , Huiying Wang , Gang Hu , Tao Liu , Jintang Qin , Yizhou Wang , Bo Xu , Jinyu Zhang , Xiaonan Shi , Jinhua Liu
The immense outburst floods that have occurred on the Tibetan Plateau during the Late Quaternary are closely linked to tectonic and climatic factors. These floods likely induced very rapid, short-term geomorphic impacts on the evolution of mountain drainage systems and patterns of sedimentary movement. In this study, four glacially-dammed lake outburst flood events that occurred along the middle reaches of the Yarlung Tsangpo River since the Middle Pleistocene were reconstructed by combining comprehensive geomorphic, stratigraphic and geochronologic investigations. The most recent outburst flood sequence occurred at ∼10.5 ka, with a peak discharge of ∼5 × 105 m3/s. The tilted uplift of the Cona Normal Fault has resulted in localized topographic lift and the formation of river knickpoints, contributing to the development and stabilization of glacial dams. River damming and outburst events have also been influenced by glacial-interglacial climate fluctuations since the Middle Pleistocene. The focused erosion and extensive mobilization of sediment by these low-frequency, high-energy floods have resulted in a repeated pattern of material transport and deposition from the Tibetan Plateau interior to its exterior. Coupled with glacial activity, the primary factor impacting the sustained stability of knickpoints on the Yarlung Tsangpo River along the Tibetan plateau's southern margin has been differential rock uplift, which results in a distinct geomorphic pattern characterized by knickpoints, glacial dams and alternating wide valleys and deep gorges.
{"title":"Catastrophic outburst floods along the middle Yarlung Tsangpo River: Responses to coupled fault and glacial activity on the Southern Tibetan Plateau","authors":"Ping Wang , Huiying Wang , Gang Hu , Tao Liu , Jintang Qin , Yizhou Wang , Bo Xu , Jinyu Zhang , Xiaonan Shi , Jinhua Liu","doi":"10.1016/j.epsl.2024.119004","DOIUrl":"10.1016/j.epsl.2024.119004","url":null,"abstract":"<div><p>The immense outburst floods that have occurred on the Tibetan Plateau during the Late Quaternary are closely linked to tectonic and climatic factors. These floods likely induced very rapid, short-term geomorphic impacts on the evolution of mountain drainage systems and patterns of sedimentary movement. In this study, four glacially-dammed lake outburst flood events that occurred along the middle reaches of the Yarlung Tsangpo River since the Middle Pleistocene were reconstructed by combining comprehensive geomorphic, stratigraphic and geochronologic investigations. The most recent outburst flood sequence occurred at ∼10.5 ka, with a peak discharge of ∼5 × 10<sup>5</sup> m<sup>3</sup>/s. The tilted uplift of the Cona Normal Fault has resulted in localized topographic lift and the formation of river knickpoints, contributing to the development and stabilization of glacial dams. River damming and outburst events have also been influenced by glacial-interglacial climate fluctuations since the Middle Pleistocene. The focused erosion and extensive mobilization of sediment by these low-frequency, high-energy floods have resulted in a repeated pattern of material transport and deposition from the Tibetan Plateau interior to its exterior. Coupled with glacial activity, the primary factor impacting the sustained stability of knickpoints on the Yarlung Tsangpo River along the Tibetan plateau's southern margin has been differential rock uplift, which results in a distinct geomorphic pattern characterized by knickpoints, glacial dams and alternating wide valleys and deep gorges.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119004"},"PeriodicalIF":4.8,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142242825","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-19DOI: 10.1016/j.epsl.2024.118980
Arianna Olivelli , Maxence Paul , Hui Xu , Katharina Kreissig , Barry J. Coles , Rebekah E.T. Moore , Luke Bridgestock , Micha Rijkenberg , Rob Middag , Maeve C. Lohan , Dominik J. Weiss , Mark Rehkämper , Tina van de Flierdt
Stable lead (Pb) isotopes have been regarded as tracers of ocean circulation, both in the present time and geological past. Here we present a new dataset of seawater Pb concentrations and isotope compositions for ten depth profiles from the South Atlantic Ocean (GEOTRACES cruises GA02 and GA10). By comparing Pb isotope data collected on the two cruises, and by modelling the distribution of Pb with an extended optimum multiparameter analysis, we find evidence of vertical transport of anthropogenic Pb pollution due to reversible scavenging. Surface to depth transfer of polluted Pb is aided by high suspended particulate matter loads at the Brazil – Malvinas Confluence and along ∼40°S in the South Atlantic. Overall, our findings caution the use of Pb isotope ratios as ventilation tracers in the South Atlantic and emphasize the importance of particle-seawater interaction for biogeochemical cycles.
{"title":"Vertical transport of anthropogenic lead by reversible scavenging in the South Atlantic Ocean","authors":"Arianna Olivelli , Maxence Paul , Hui Xu , Katharina Kreissig , Barry J. Coles , Rebekah E.T. Moore , Luke Bridgestock , Micha Rijkenberg , Rob Middag , Maeve C. Lohan , Dominik J. Weiss , Mark Rehkämper , Tina van de Flierdt","doi":"10.1016/j.epsl.2024.118980","DOIUrl":"10.1016/j.epsl.2024.118980","url":null,"abstract":"<div><p>Stable lead (Pb) isotopes have been regarded as tracers of ocean circulation, both in the present time and geological past. Here we present a new dataset of seawater Pb concentrations and isotope compositions for ten depth profiles from the South Atlantic Ocean (GEOTRACES cruises GA02 and GA10). By comparing Pb isotope data collected on the two cruises, and by modelling the distribution of Pb with an extended optimum multiparameter analysis, we find evidence of vertical transport of anthropogenic Pb pollution due to reversible scavenging. Surface to depth transfer of polluted Pb is aided by high suspended particulate matter loads at the Brazil – Malvinas Confluence and along ∼40°S in the South Atlantic. Overall, our findings caution the use of Pb isotope ratios as ventilation tracers in the South Atlantic and emphasize the importance of particle-seawater interaction for biogeochemical cycles.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 118980"},"PeriodicalIF":4.8,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004126/pdfft?md5=b2bd24798be54fe416b91b92712f41be&pid=1-s2.0-S0012821X24004126-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142247878","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-19DOI: 10.1016/j.epsl.2024.119012
Ruirui Meng , Zuozhen Han , Xiao Gao , Yanyang Zhao , Chao Han , Yu Han , Renchao Yang , Shuangjian Li , Fang Liu , Maurice E. Tucker , Yanjing Chen
Marine carbonate rocks are the major reservoir of carbon through CaCO3 and CaMg(CO3)2 precipitation extracting CO2/HCO3- from the atmosphere/oceans; hence dolomite is one of the major sinks in the global carbon cycle. Most ancient dolomite has been considered as mainly precipitated under Earth surface conditions. Previous studies have demonstrated that microbes can mediate dolomite formation. However, the “microbial dolomite” model is not sufficient to explain the dolomite that shows no or little evidence of a microbial origin. Although attempted for decades, the synthesis of inorganically-produced dolomite at normal temperatures (<50 °C) has been relatively unsuccessful. Hence, this "dolomite enigma" has been one principal research focus this century. Here we demonstrate that NH3 catalyzes proto-dolomite precipitation inorganically with higher Mg/Ca molar ratios and CO32- activity at normal temperatures of 30 °C and 40 °C. NH3 dissolution increases the alkalinity of the solution and transformes into NH4+ ions, which prefer to bond with H2O on Mg[(H2O)6]2+ rather than free H2O, thus releasing Mg2+ to facilitate proto-dolomite nucleation. Furthermore, the low dielectric constant and low dipole moment allow NH4+ absorbed on crystal surfaces to lower the energy barrier of Mg[(H2O)6]2+ dehydration, promoting proto-dolomite nucleis growth. The system for proto-dolomite precipitation in our experiments closely simulates the natural aqueous environment. This study brings new insights to understanding the mechanisms of dolomite precipitation in natural waters.
{"title":"Dissolved ammonia catalyzes proto-dolomite precipitation at Earth surface temperature","authors":"Ruirui Meng , Zuozhen Han , Xiao Gao , Yanyang Zhao , Chao Han , Yu Han , Renchao Yang , Shuangjian Li , Fang Liu , Maurice E. Tucker , Yanjing Chen","doi":"10.1016/j.epsl.2024.119012","DOIUrl":"10.1016/j.epsl.2024.119012","url":null,"abstract":"<div><p>Marine carbonate rocks are the major reservoir of carbon through CaCO<sub>3</sub> and CaMg(CO<sub>3</sub>)<sub>2</sub> precipitation extracting CO<sub>2</sub>/HCO<sub>3</sub><sup>-</sup> from the atmosphere/oceans; hence dolomite is one of the major sinks in the global carbon cycle. Most ancient dolomite has been considered as mainly precipitated under Earth surface conditions. Previous studies have demonstrated that microbes can mediate dolomite formation. However, the “microbial dolomite” model is not sufficient to explain the dolomite that shows no or little evidence of a microbial origin. Although attempted for decades, the synthesis of inorganically-produced dolomite at normal temperatures (<50 °C) has been relatively unsuccessful. Hence, this \"dolomite enigma\" has been one principal research focus this century. Here we demonstrate that NH<sub>3</sub> catalyzes proto-dolomite precipitation inorganically with higher Mg/Ca molar ratios and CO<sub>3</sub><sup>2-</sup> activity at normal temperatures of 30 °C and 40 °C. NH<sub>3</sub> dissolution increases the alkalinity of the solution and transformes into NH<sub>4</sub><sup>+</sup> ions, which prefer to bond with H<sub>2</sub>O on Mg[(H<sub>2</sub>O)<sub>6</sub>]<sup>2+</sup> rather than free H<sub>2</sub>O, thus releasing Mg<sup>2+</sup> to facilitate proto-dolomite nucleation. Furthermore, the low dielectric constant and low dipole moment allow NH<sub>4</sub><sup>+</sup> absorbed on crystal surfaces to lower the energy barrier of Mg[(H<sub>2</sub>O)<sub>6</sub>]<sup>2+</sup> dehydration, promoting proto-dolomite nucleis growth. The system for proto-dolomite precipitation in our experiments closely simulates the natural aqueous environment. This study brings new insights to understanding the mechanisms of dolomite precipitation in natural waters.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119012"},"PeriodicalIF":4.8,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24004448/pdfft?md5=996a34f1f457951db14bb28351bd22e6&pid=1-s2.0-S0012821X24004448-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142247877","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-18DOI: 10.1016/j.epsl.2024.119008
Hairuo Fu, Stein B. Jacobsen
The canonical Moon-forming giant-impact models allow for substantial chemical differences between the bulk silicate Moon and Earth due to incomplete mixing of the impactor and the proto-Earth. In comparison, the emerging high-energy giant-impact (Synestia) model requires the refractory element compositions of the Earth and Moon to be nearly identical, owing to extensive chemical homogenization of the Moon-forming disk in a vigorously mixed silicate fluid. These distinct chemical predictions make the lunar refractory element composition crucial for testing Moon-formation hypotheses, yet it remains highly controversial and necessitates new approaches to resolve. In this study, we develop a novel method using the composition of pristine lunar anorthosite samples to constrain the Moon's refractory lithophile element compositions. We obtained a very close match of refractory major and trace element compositions for the lunar magma ocean model, suggesting indistinguishable refractory element abundances between the bulk silicate Moon and Earth. This striking refractory element similarity is difficult to reconcile with the relatively poor mixing conditions of the canonical giant-impact models. The compatibility of this result with disk equilibration models other than the Synestia has yet to be quantitatively verified. Our results further constrain that the formation of the Earth-Moon system requires a thoroughly-mixed protolunar disk of chemical and isotopic homogenization with an initially fully-molten Moon, as enabled by emerging models like the Synestia.
{"title":"Earth-Moon refractory element similarity constrains a thoroughly-mixed Moon-forming disk","authors":"Hairuo Fu, Stein B. Jacobsen","doi":"10.1016/j.epsl.2024.119008","DOIUrl":"10.1016/j.epsl.2024.119008","url":null,"abstract":"<div><p>The canonical Moon-forming giant-impact models allow for substantial chemical differences between the bulk silicate Moon and Earth due to incomplete mixing of the impactor and the proto-Earth. In comparison, the emerging high-energy giant-impact (Synestia) model requires the refractory element compositions of the Earth and Moon to be nearly identical, owing to extensive chemical homogenization of the Moon-forming disk in a vigorously mixed silicate fluid. These distinct chemical predictions make the lunar refractory element composition crucial for testing Moon-formation hypotheses, yet it remains highly controversial and necessitates new approaches to resolve. In this study, we develop a novel method using the composition of pristine lunar anorthosite samples to constrain the Moon's refractory lithophile element compositions. We obtained a very close match of refractory major and trace element compositions for the lunar magma ocean model, suggesting indistinguishable refractory element abundances between the bulk silicate Moon and Earth. This striking refractory element similarity is difficult to reconcile with the relatively poor mixing conditions of the canonical giant-impact models. The compatibility of this result with disk equilibration models other than the Synestia has yet to be quantitatively verified. Our results further constrain that the formation of the Earth-Moon system requires a thoroughly-mixed protolunar disk of chemical and isotopic homogenization with an initially fully-molten Moon, as enabled by emerging models like the Synestia.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 119008"},"PeriodicalIF":4.8,"publicationDate":"2024-09-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142243045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
<div><p>Rubidium (Rb) is a strongly incompatible and highly fluid-mobile element and Rb isotopes have the potential to track subducted material recycling to the mantle source of arc magmas. However, the behavior of Rb isotopes during slab subduction and associated processes remains unconstrained. This study for the first time presents Rb isotope data for the altered oceanic crust (AOC) from IODP Site U1365, subducted sediments, and lavas from the Southern Mariana arc. The δ<sup>87</sup>Rb of the AOC varies between –0.16 ± 0.03 ‰ to 0.07 ± 0.02 ‰ with a weighted average of –0.02 ± 0.06 ‰, which is higher than that of fresh mid-ocean ridge basalts (MORB; –0.12 ± 0.08 ‰). Such Rb isotope fractionation may be attributed to the preferential loss of <sup>85</sup>Rb into seawater during the dissolution of primary phases and the preferential structural incorporation of <sup>87</sup>Rb into secondary phases and preferential <sup>87</sup>Rb adsorption by clays. The δ<sup>87</sup>Rb of the sediments varies from –0.20 ± 0.07 ‰ to –0.03 ± 0.02 ‰ with an average of –0.11 ± 0.12 ‰, identical to that of the upper continental crust (UCC; –0.14 ± 0.08 ‰). The correlations observed between δ<sup>87</sup>Rb and sediment depth (and loss on ignition; LOI) suggest that hydrodynamic sorting of sediments which enriches the upper stratigraphic intervals in clays, exerts a first-order control on the Rb isotope fractionation in subducted sediments. The formation of authigenic clays, which act as a sink of <sup>85</sup>Rb from seawater, could partially account for the elevated δ<sup>87</sup>Rb of seawater (0.14 ± 0.12 ‰) relative to the Bulk Silicate Earth (BSE; –0.12 ± 0.06 ‰). For the Mariana lavas, the δ<sup>87</sup>Rb decreases from fore-arc lavas (–0.03 ± 0.04 ‰ to 0.09 ± 0.02 ‰) to frontal arc lavas (–0.12 ± 0.04 ‰ to –0.01 ± 0.03 ‰), which are all higher than the average δ<sup>87</sup>Rb of fresh MORB (–0.12 ± 0.08 ‰). The elevated δ<sup>87</sup>Rb of the lavas reflects the addition of AOC-derived aqueous fluids with a high δ<sup>87</sup>Rb (∼0.08 ‰) to the magma source, which is higher than the average δ<sup>87</sup>Rb (–0.02 ± 0.18 ‰) of the Site-1365 AOC. The variation of δ<sup>87</sup>Rb with <sup>143</sup>Nd/<sup>144</sup>Nd<sub>i</sub> values of the lavas indicates that AOC-derived aqueous fluids are increasingly mixed with a component with a low δ<sup>87</sup>Rb that most likely represents sediment melts. The estimated δ<sup>87</sup>Rb of sediment melts is about –0.27 ‰, lower than the average δ<sup>87</sup>Rb (–0.11 ± 0.12 ‰) of the Mariana sediments. Thus, our study suggests that Rb isotopes can be fractionated during AOC dehydration and sediment partial melting. Overall, the across-arc Rb isotope variation reflects slab dehydration and melting at different subduction depths. Our findings demonstrate that Rb isotopes are a powerful novel tool for probing deep recycling of subducted slabs. Furthermore, they can distinguish contributions from AOC-derived fluids
{"title":"Rubidium isotopes reveal dehydration and melting of the subducting slab beneath the Mariana arc","authors":"Ding-Sheng Jiang , Xiao-Tong Peng , Wen-Jing Xu , Xia Hu , Saskia Erdmann , Xi-Sheng Xu , Guo-Liang Zhang , Chong-Jin Pang , Hao-Chen Duan , Fang Huang","doi":"10.1016/j.epsl.2024.118978","DOIUrl":"10.1016/j.epsl.2024.118978","url":null,"abstract":"<div><p>Rubidium (Rb) is a strongly incompatible and highly fluid-mobile element and Rb isotopes have the potential to track subducted material recycling to the mantle source of arc magmas. However, the behavior of Rb isotopes during slab subduction and associated processes remains unconstrained. This study for the first time presents Rb isotope data for the altered oceanic crust (AOC) from IODP Site U1365, subducted sediments, and lavas from the Southern Mariana arc. The δ<sup>87</sup>Rb of the AOC varies between –0.16 ± 0.03 ‰ to 0.07 ± 0.02 ‰ with a weighted average of –0.02 ± 0.06 ‰, which is higher than that of fresh mid-ocean ridge basalts (MORB; –0.12 ± 0.08 ‰). Such Rb isotope fractionation may be attributed to the preferential loss of <sup>85</sup>Rb into seawater during the dissolution of primary phases and the preferential structural incorporation of <sup>87</sup>Rb into secondary phases and preferential <sup>87</sup>Rb adsorption by clays. The δ<sup>87</sup>Rb of the sediments varies from –0.20 ± 0.07 ‰ to –0.03 ± 0.02 ‰ with an average of –0.11 ± 0.12 ‰, identical to that of the upper continental crust (UCC; –0.14 ± 0.08 ‰). The correlations observed between δ<sup>87</sup>Rb and sediment depth (and loss on ignition; LOI) suggest that hydrodynamic sorting of sediments which enriches the upper stratigraphic intervals in clays, exerts a first-order control on the Rb isotope fractionation in subducted sediments. The formation of authigenic clays, which act as a sink of <sup>85</sup>Rb from seawater, could partially account for the elevated δ<sup>87</sup>Rb of seawater (0.14 ± 0.12 ‰) relative to the Bulk Silicate Earth (BSE; –0.12 ± 0.06 ‰). For the Mariana lavas, the δ<sup>87</sup>Rb decreases from fore-arc lavas (–0.03 ± 0.04 ‰ to 0.09 ± 0.02 ‰) to frontal arc lavas (–0.12 ± 0.04 ‰ to –0.01 ± 0.03 ‰), which are all higher than the average δ<sup>87</sup>Rb of fresh MORB (–0.12 ± 0.08 ‰). The elevated δ<sup>87</sup>Rb of the lavas reflects the addition of AOC-derived aqueous fluids with a high δ<sup>87</sup>Rb (∼0.08 ‰) to the magma source, which is higher than the average δ<sup>87</sup>Rb (–0.02 ± 0.18 ‰) of the Site-1365 AOC. The variation of δ<sup>87</sup>Rb with <sup>143</sup>Nd/<sup>144</sup>Nd<sub>i</sub> values of the lavas indicates that AOC-derived aqueous fluids are increasingly mixed with a component with a low δ<sup>87</sup>Rb that most likely represents sediment melts. The estimated δ<sup>87</sup>Rb of sediment melts is about –0.27 ‰, lower than the average δ<sup>87</sup>Rb (–0.11 ± 0.12 ‰) of the Mariana sediments. Thus, our study suggests that Rb isotopes can be fractionated during AOC dehydration and sediment partial melting. Overall, the across-arc Rb isotope variation reflects slab dehydration and melting at different subduction depths. Our findings demonstrate that Rb isotopes are a powerful novel tool for probing deep recycling of subducted slabs. Furthermore, they can distinguish contributions from AOC-derived fluids","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 118978"},"PeriodicalIF":4.8,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142243886","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-09-16DOI: 10.1016/j.epsl.2024.118989
Yishen Zhang , Bernard Charlier , Stephanie B. Krein , Timothy L. Grove , Olivier Namur , Francois Holtz
The latest stages of the lunar magma ocean (LMO) crystallization led to the formation of ilmenite-bearing cumulates and urKREEP, residual melts enriched in K, rare earth elements (REEs), P, and other incompatible elements. Those highly evolved lithologies had major impacts on the petrogenesis of lunar volcanic rocks and the compositional diversity of post-LMO magmatism resulting from mantle remelting. Here, we present new experimental results constraining the composition of the very last liquids produced during LMO crystallization. To test the potential role of silicate liquid immiscibility in the formation of urKREEP, synthetic samples representative of residual melts of bulk Moon compositions were placed in double platinum-graphite capsules at 1020–980 °C and 0.08–0.10 GPa in an internally-heated pressure vessel. The produced silicate liquids are multiply saturated with plagioclase, augite, silica phases, and ilmenite (± fayalitic olivine ± pigeonite). Our experiments show that the liquid line of descent reaches a two-liquid field at 1000 °C and >97% crystallization for a range of whole-Moon compositions. Under these conditions, a small proportion of silica-rich melt (70.0–71.4 wt.% SiO2, 6.4–7.3 wt.% FeO, 5.4–6.1 wt.% K2O, 0.2–0.3 wt.% P2O5) coexists within an abundant Fe-rich melt (42.6–44.1 wt.% SiO2, 27.6–28.8 wt.% FeO, 0.9–1.0 wt.% K2O, 2.8–3.2 wt.% P2O5) with sharp two-liquid interfaces. Our experimental results also constrain the relative onset of ilmenite crystallization compared to the development of immiscibility and indicate that an ilmenite-bearing layer formed in the lunar interior before immiscibility was attained. Using a self-consistent physicochemical LMO model, we constrain the thickness and depth of the ilmenite-bearing layer during LMO differentiation. The immiscible K-Si-rich and P-Fe-rich melts together also produced an immiscible urKREEP layer ∼2–6 km thick and ∼30–50 km deep depending on the trapped liquid fraction in the cumulate column (≤10%) and the thickness of the buoyant anorthosite crust (30–50 km). We provide constraints on the relationship between the compositions of immiscible urKREEP melts and those of KREEPy rocks. By modeling the mixing of KREEP-poor basalt and the immiscible melt pairs, we reproduce the K and P enrichments and apparent decoupling of K from P in KREEPy rocks. Our results highlight that processes such as the assimilation of evolved heterogeneous mantle lithologies may be involved in hybridization during post-LMO magmatism. The immiscible K-Si-rich lithology may also have contributed to lunar silicic magmatism.
{"title":"The very late-stage crystallization of the lunar magma ocean and the composition of immiscible urKREEP","authors":"Yishen Zhang , Bernard Charlier , Stephanie B. Krein , Timothy L. Grove , Olivier Namur , Francois Holtz","doi":"10.1016/j.epsl.2024.118989","DOIUrl":"10.1016/j.epsl.2024.118989","url":null,"abstract":"<div><p>The latest stages of the lunar magma ocean (LMO) crystallization led to the formation of ilmenite-bearing cumulates and urKREEP, residual melts enriched in K, rare earth elements (REEs), P, and other incompatible elements. Those highly evolved lithologies had major impacts on the petrogenesis of lunar volcanic rocks and the compositional diversity of post-LMO magmatism resulting from mantle remelting. Here, we present new experimental results constraining the composition of the very last liquids produced during LMO crystallization. To test the potential role of silicate liquid immiscibility in the formation of urKREEP, synthetic samples representative of residual melts of bulk Moon compositions were placed in double platinum-graphite capsules at 1020–980 °C and 0.08–0.10 GPa in an internally-heated pressure vessel. The produced silicate liquids are multiply saturated with plagioclase, augite, silica phases, and ilmenite (± fayalitic olivine ± pigeonite). Our experiments show that the liquid line of descent reaches a two-liquid field at 1000 °C and >97% crystallization for a range of whole-Moon compositions. Under these conditions, a small proportion of silica-rich melt (70.0–71.4 wt.% SiO<sub>2</sub>, 6.4–7.3 wt.% FeO, 5.4–6.1 wt.% K<sub>2</sub>O, 0.2–0.3 wt.% P<sub>2</sub>O<sub>5</sub>) coexists within an abundant Fe-rich melt (42.6–44.1 wt.% SiO<sub>2</sub>, 27.6–28.8 wt.% FeO, 0.9–1.0 wt.% K<sub>2</sub>O, 2.8–3.2 wt.% P<sub>2</sub>O<sub>5</sub>) with sharp two-liquid interfaces. Our experimental results also constrain the relative onset of ilmenite crystallization compared to the development of immiscibility and indicate that an ilmenite-bearing layer formed in the lunar interior before immiscibility was attained. Using a self-consistent physicochemical LMO model, we constrain the thickness and depth of the ilmenite-bearing layer during LMO differentiation. The immiscible K-Si-rich and P-Fe-rich melts together also produced an immiscible urKREEP layer ∼2–6 km thick and ∼30–50 km deep depending on the trapped liquid fraction in the cumulate column (≤10%) and the thickness of the buoyant anorthosite crust (30–50 km). We provide constraints on the relationship between the compositions of immiscible urKREEP melts and those of KREEPy rocks. By modeling the mixing of KREEP-poor basalt and the immiscible melt pairs, we reproduce the K and P enrichments and apparent decoupling of K from P in KREEPy rocks. Our results highlight that processes such as the assimilation of evolved heterogeneous mantle lithologies may be involved in hybridization during post-LMO magmatism. The immiscible K-Si-rich lithology may also have contributed to lunar silicic magmatism.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"646 ","pages":"Article 118989"},"PeriodicalIF":4.8,"publicationDate":"2024-09-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142233680","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}