Pub Date : 2024-06-03DOI: 10.1016/j.epsl.2024.118806
Guang-Yi Wei, Feifei Zhang
Early marine and meteoric diagenesis is commonly proposed to alter the primary geochemical signals of shallow-marine carbonates, hindering their use as a reliable archive for ancient seawater chemistry. However, there are still uncertainties on the extent and controlling factors of early diagenetic alteration for shallow-marine carbonates. In this study, we compile and compare radiogenic Sr, Li, and U isotope data of carbonate island reefs and platforms from the South China Sea and the Bahamas to better constrain the early diagenetic effects on geochemical signals in carbonates. Multiple diagenetic zones are identified in carbonate drillcores from the South China Sea, and temporally correlated with similar zones in the Bahamas. The radiogenic Sr, Li and U isotope records from shallow-marine carbonates in these two geographically distant regions are affected by multiple factors, including compositions of diagenetic fluids, local depositional environments, fluid advection/diffusion rates, and carbonate recrystallization rates. These factors ultimately determined the diagenetic extents and regimes (fluid-buffered or sediment-buffered) of shallow-marine carbonates, which vary across different isotopic systems. A full understanding of early diagenetic history and its effects on different isotopic proxies will validate the use of shallow-marine carbonate as an archive for ancient seawater information, and contribute to tracking the temporal changes in paleo-climate and -environment.
人们普遍认为,早期海洋和陨石成岩作用改变了浅海碳酸盐的初级地球化学信号,阻碍了它们作为古海水化学可靠档案的使用。然而,浅海碳酸盐早期成岩蚀变的程度和控制因素仍存在不确定性。在本研究中,我们汇编并比较了南海和巴哈马群岛碳酸盐岛礁和平台的放射性 Sr、Li 和 U 同位素数据,以更好地确定早期成岩作用对碳酸盐地球化学信号的影响。在中国南海的碳酸盐岩钻屑中发现了多个成岩带,并与巴哈马群岛的类似成岩带在时间上相互关联。这两个地理位置遥远地区浅海碳酸盐的放射性 Sr、Li 和 U 同位素记录受到多种因素的影响,包括成岩流体成分、当地沉积环境、流体平流/扩散速率以及碳酸盐再结晶速率。这些因素最终决定了浅海碳酸盐岩的成岩范围和机制(流体缓冲或沉积物缓冲),不同同位素系统的成岩范围和机制各不相同。充分了解早期成岩历史及其对不同同位素代用指标的影响,将验证将浅海碳酸盐作为古海水信息档案的有效性,并有助于追踪古气候和古环境的时间变化。
{"title":"Pristine or altered, what can early diagenesis tell us in shallow-water carbonates?","authors":"Guang-Yi Wei, Feifei Zhang","doi":"10.1016/j.epsl.2024.118806","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118806","url":null,"abstract":"<div><p>Early marine and meteoric diagenesis is commonly proposed to alter the primary geochemical signals of shallow-marine carbonates, hindering their use as a reliable archive for ancient seawater chemistry. However, there are still uncertainties on the extent and controlling factors of early diagenetic alteration for shallow-marine carbonates. In this study, we compile and compare radiogenic Sr, Li, and U isotope data of carbonate island reefs and platforms from the South China Sea and the Bahamas to better constrain the early diagenetic effects on geochemical signals in carbonates. Multiple diagenetic zones are identified in carbonate drillcores from the South China Sea, and temporally correlated with similar zones in the Bahamas. The radiogenic Sr, Li and U isotope records from shallow-marine carbonates in these two geographically distant regions are affected by multiple factors, including compositions of diagenetic fluids, local depositional environments, fluid advection/diffusion rates, and carbonate recrystallization rates. These factors ultimately determined the diagenetic extents and regimes (fluid-buffered or sediment-buffered) of shallow-marine carbonates, which vary across different isotopic systems. A full understanding of early diagenetic history and its effects on different isotopic proxies will validate the use of shallow-marine carbonate as an archive for ancient seawater information, and contribute to tracking the temporal changes in paleo-climate and -environment.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141241990","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-06-01DOI: 10.1016/j.epsl.2024.118799
Jiangyang Zhang , Fan Zhang , Jian Lin , Xiang Gao , Chen Cai , Zhiyuan Zhou
A primary pathway for surface water to enter the Earth's interior is through faulting of oceanic plates. Mantle hydration of subducting plate is considered to play an important role on the water flux of subduction zone. Previous studies have found that the mantle hydration was related to plate bending and faulting, however the controlling mechanism of plate bending on mantle hydration at subduction zone is still highly unclear. In this study, we use the latest obtained reduction in uppermost mantle P-wave velocity (Vp) profiles beneath the subduction zone outer rise region, which is thought to be the result of mantle serpentinization, to examine the degrees of mantle serpentinization at several subduction zones and analyzed their relationship to plate age, bending curvature, and the sediment thickness near the trench. Results of analysis revealed that the average degree of the estimated mantle serpentinization increases with both the plate age and bending curvature, while it is hampered by sedimentation. The thick sediment can almost entirely prevent water from entering the mantle. Importantly, we found a good linear correlation between the reduction of uppermost mantle Vp and the theoretical brittle extensional strain determined by plate age and curvature. Based on this, a mechanism is proposed to illustrate the governing mechanism of the combined effect of the plate age and bending curvature in controlling the mantle serpentinization degree for global subduction zones.
{"title":"Mantle serpentinization of subducting plate are controlled by combined effect of plate age and bending curvature","authors":"Jiangyang Zhang , Fan Zhang , Jian Lin , Xiang Gao , Chen Cai , Zhiyuan Zhou","doi":"10.1016/j.epsl.2024.118799","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118799","url":null,"abstract":"<div><p>A primary pathway for surface water to enter the Earth's interior is through faulting of oceanic plates. Mantle hydration of subducting plate is considered to play an important role on the water flux of subduction zone. Previous studies have found that the mantle hydration was related to plate bending and faulting, however the controlling mechanism of plate bending on mantle hydration at subduction zone is still highly unclear. In this study, we use the latest obtained reduction in uppermost mantle <em>P-</em>wave velocity (<em>V</em><sub>p</sub>) profiles beneath the subduction zone outer rise region, which is thought to be the result of mantle serpentinization, to examine the degrees of mantle serpentinization at several subduction zones and analyzed their relationship to plate age, bending curvature, and the sediment thickness near the trench. Results of analysis revealed that the average degree of the estimated mantle serpentinization increases with both the plate age and bending curvature, while it is hampered by sedimentation. The thick sediment can almost entirely prevent water from entering the mantle. Importantly, we found a good linear correlation between the reduction of uppermost mantle <em>V</em><sub>p</sub> and the theoretical brittle extensional strain determined by plate age and curvature. Based on this, a mechanism is proposed to illustrate the governing mechanism of the combined effect of the plate age and bending curvature in controlling the mantle serpentinization degree for global subduction zones.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141240628","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-06-01DOI: 10.1016/j.epsl.2024.118767
Subhajit Ghosh , Holger Stünitz , Hugues Raimbourg , Jacques Précigout , Ida Di Carlo , Renée Heilbronner , Laurette Piani
When H2O is present along grain boundaries, the deformation processes responsible for plasticity in silicate mineral aggregates can deviate from what may be conventionally expected. Although a necessary component of understanding crustal deformation processes, there is no theoretical framework that incorporates grain boundary processes into polycrystalline quartz rheology. To address this issue, we carried out high-pressure and high-temperature deformation experiments on fine-grained quartz aggregates. Our study illustrates that grain boundary migration (GBM) through dissolution-precipitation (in the presence of an aqueous fluid) and grain boundary sliding (GBS) may act as accommodation mechanisms to prevent hardening from dislocation glide. GBM and GBS can relax incompatibilities resulting from an inadequate number of independent slip systems, plastic anisotropy between neighbouring grains, and non-planar grain boundaries together with grain boundary junctions. As demonstrated earlier in the literature, GBM may act as a recrystallization mechanism counteracting hardening, but also is a potential mechanism that allow H2O to enter in the quartz crystal (hydrolization) at the experimental time-scale. The above serial processes occur over a range of more than two orders of magnitude in grain size (∼3 to 200 μm) and explain a grain-size-insensitive stress exponent (n = 2) and low activation energy (Q = 110 kJ/mol). In the absence of a switch to grain size sensitive deformation mechanisms induced by grain size reduction, our results imply that only a modest weakening (∼5 times the strength of the protolith) is needed (or possible) to localize shear zones in the Earth's crust.
{"title":"Importance of grain boundary processes for plasticity in the quartz-dominated crust: Implications for flow laws","authors":"Subhajit Ghosh , Holger Stünitz , Hugues Raimbourg , Jacques Précigout , Ida Di Carlo , Renée Heilbronner , Laurette Piani","doi":"10.1016/j.epsl.2024.118767","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118767","url":null,"abstract":"<div><p>When H<sub>2</sub>O is present along grain boundaries, the deformation processes responsible for plasticity in silicate mineral aggregates can deviate from what may be conventionally expected. Although a necessary component of understanding crustal deformation processes, there is no theoretical framework that incorporates grain boundary processes into polycrystalline quartz rheology. To address this issue, we carried out high-pressure and high-temperature deformation experiments on fine-grained quartz aggregates. Our study illustrates that grain boundary migration (GBM) through dissolution-precipitation (in the presence of an aqueous fluid) and grain boundary sliding (GBS) may act as accommodation mechanisms to prevent hardening from dislocation glide. GBM and GBS can relax incompatibilities resulting from an inadequate number of independent slip systems, plastic anisotropy between neighbouring grains, and non-planar grain boundaries together with grain boundary junctions. As demonstrated earlier in the literature, GBM may act as a recrystallization mechanism counteracting hardening, but also is a potential mechanism that allow H<sub>2</sub>O to enter in the quartz crystal (hydrolization) at the experimental time-scale. The above serial processes occur over a range of more than two orders of magnitude in grain size (∼3 to 200 μm) and explain a grain-size-insensitive stress exponent (<em>n</em> = 2) and low activation energy (<em>Q</em> = 110 kJ/mol). In the absence of a switch to grain size sensitive deformation mechanisms induced by grain size reduction, our results imply that only a modest weakening (∼5 times the strength of the protolith) is needed (or possible) to localize shear zones in the Earth's crust.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244459","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-05-31DOI: 10.1016/j.epsl.2024.118780
Charles D. Beard , Adrian A. Finch , Anouk M. Borst , Kathryn M. Goodenough , William Hutchison , Ian L. Millar , Tom Andersen , Helen M. Williams , Owen M. Weller
Alkaline-silicate complexes host some of the world's largest resources of rare-earth elements and high-field-strength elements (REE & HFSE) and represent the most fractionated magmatic systems on our planet. Geochemical evidence indicates that they are mantle melts, but while various studies highlight a role for lithospheric mantle, we do not know the precise origin of their contained REE and HFSE, and whether enrichment of the mantle source for these magmas can be attributed to specific geodynamic processes or events.
We present new Nd-Hf isotope measurements (/ & / ) made by LA-MC-ICP-MS, as well as a compilation of existing isotopic data for a suite of alkaline igneous rocks from the Gardar Province, a Mesoproterozoic continental rift in southern Greenland. Neodymium and hafnium isotopes are unaffected by crystal fractionation and can directly fingerprint the source of REE and HFSE. The dataset covers both phases of Gardar magmatism (1325–1261 and 1184–1140 Ma) and incorporates mafic dyke swarms and km-scale intrusive complexes, including Ilimmaasaq (Ilímaussaq) and Motzfeldt, which host some of the world's largest REE and HFSE deposits. The majority of Gardar complexes have a narrow range of positive median initial εNd (0 to +3.3) and εHf values (+0.2 to +6.0). Only two granite intrusions and the Eriksfjord basaltic lavas have crustally contaminated Nd-Hf isotope compositions, with the vast majority of Gardar igneous rocks preserving the isotope signature of their mantle source. Considering the diversity of rock types in the Gardar Province, initial εNd -εHf compositions are remarkably homogeneous, indicating a derivation of the Gardar's REE and HFSE from a laterally-extensive mantle melt source.
Several Gardar systems have low initial εHf for a given εNd (εHf to -9.7), a distinctive signature as few geological processes decouple the Nd and Hf isotope systems. The decoupled Nd-Hf isotope signatures are consistent with contributions from isotopically-matured phlogopite-bearing metasomatic veins (commonly known as PIC: phlogopite-ilmenite-clinopyroxene) in the lithospheric mantle. The metasomatising fluids that formed these source rocks were introduced via Palaeoproterozoic subduction, but the Gardar isotopic signatures indicate that REE and HFSE enrichment of these metasomes was not derived from subducted sediment; instead it is likely that metals were scavenged from the mantle wedge overlying the ancient subduction zone. The Gardar Nd-Hf isotope evolution trends overlap with a global compilation of kimberlites through time and allow us to tie the origin of the PIC metasomes to the regional geodynamic history of South Greenland. We identify PIC metasomes as a key metal s
{"title":"A phlogopite-bearing lithospheric mantle source for Europe's largest REE-HFSE belt: Gardar Rift, SW Greenland","authors":"Charles D. Beard , Adrian A. Finch , Anouk M. Borst , Kathryn M. Goodenough , William Hutchison , Ian L. Millar , Tom Andersen , Helen M. Williams , Owen M. Weller","doi":"10.1016/j.epsl.2024.118780","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118780","url":null,"abstract":"<div><p>Alkaline-silicate complexes host some of the world's largest resources of rare-earth elements and high-field-strength elements (REE & HFSE) and represent the most fractionated magmatic systems on our planet. Geochemical evidence indicates that they are mantle melts, but while various studies highlight a role for lithospheric mantle, we do not know the precise origin of their contained REE and HFSE, and whether enrichment of the mantle source for these magmas can be attributed to specific geodynamic processes or events.</p><p>We present new Nd-Hf isotope measurements (<figure><img></figure>/<figure><img></figure> & <figure><img></figure>/<figure><img></figure> ) made by LA-MC-ICP-MS, as well as a compilation of existing isotopic data for a suite of alkaline igneous rocks from the Gardar Province, a Mesoproterozoic continental rift in southern Greenland. Neodymium and hafnium isotopes are unaffected by crystal fractionation and can directly fingerprint the source of REE and HFSE. The dataset covers both phases of Gardar magmatism (1325–1261 and 1184–1140 Ma) and incorporates mafic dyke swarms and km-scale intrusive complexes, including Ilimmaasaq (Ilímaussaq) and Motzfeldt, which host some of the world's largest REE and HFSE deposits. The majority of Gardar complexes have a narrow range of positive median initial <em>ε</em><sub>Nd</sub> (0 to +3.3) and <em>ε</em><sub>Hf</sub> values (+0.2 to +6.0). Only two granite intrusions and the Eriksfjord basaltic lavas have crustally contaminated Nd-Hf isotope compositions, with the vast majority of Gardar igneous rocks preserving the isotope signature of their mantle source. Considering the diversity of rock types in the Gardar Province, initial <em>ε</em><sub>Nd</sub> -<em>ε</em><sub>Hf</sub> compositions are remarkably homogeneous, indicating a derivation of the Gardar's REE and HFSE from a laterally-extensive mantle melt source.</p><p>Several Gardar systems have low initial <em>ε</em><sub>Hf</sub> for a given <em>ε</em><sub>Nd</sub> (<figure><img></figure><em>ε</em><sub>Hf</sub> to -9.7), a distinctive signature as few geological processes decouple the Nd and Hf isotope systems. The decoupled Nd-Hf isotope signatures are consistent with contributions from isotopically-matured phlogopite-bearing metasomatic veins (commonly known as PIC: phlogopite-ilmenite-clinopyroxene) in the lithospheric mantle. The metasomatising fluids that formed these source rocks were introduced via Palaeoproterozoic subduction, but the Gardar isotopic signatures indicate that REE and HFSE enrichment of these metasomes was not derived from subducted sediment; instead it is likely that metals were scavenged from the mantle wedge overlying the ancient subduction zone. The Gardar Nd-Hf isotope evolution trends overlap with a global compilation of kimberlites through time and allow us to tie the origin of the PIC metasomes to the regional geodynamic history of South Greenland. We identify PIC metasomes as a key metal s","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24002139/pdfft?md5=ff5d376be64564e07f6cecbb0a301879&pid=1-s2.0-S0012821X24002139-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244462","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-05-31DOI: 10.1016/j.epsl.2024.118795
Lori N Willhite, Valerie A Finlayson, Richard J Walker
Highly siderophile element abundances and 182W/184W and 187Os/188Os were determined for a suite of Mauna Kea lavas from the Hawaiian Scientific Drilling Project phase 2 drill core. The new analyses, combined with previous measurements, compose the largest database for µ182W (the parts-per-million deviation of 182W/184W from a terrestrial standard) for a single volcano (n= 16). Although most lavas analyzed are characterized by negative µ182W values, lavas with values similar to the modern bulk silicate Earth are found throughout the entire stratigraphic column. This suggests that components with normal µ182W are collocated with components that host µ182W deficits in the plume. Negative µ182W values are associated with elevated 3He/4He, as well as elevated Ti and Nb. These correlations may link µ182W anomalies to ancient deep mantle crystal-liquid fractionation processes. Consistent with previously measured 3He/4He (R/RA) in the drill core, the magnitude of negative µ182W values was greatest when Mauna Kea was close to the plume axis then generally decreased over the ∼400 kyr captured by the stratigraphic section. The component with anomalous µ182W was either concentrated near the plume axis, or was more effectively sampled by melting near the plume axis where the temperature excess was greatest, suggesting it was less fusible than the dominant plume components. The process leading to the generation of a mantle component with a negative µ182W anomaly could either be related to some form of core-mantle isotopic equilibration, or early-Earth fractionation within the silicate Earth. At present each possibility remains viable.
{"title":"Evolution of tungsten isotope systematics in the Mauna Kea volcano provides new constraints on anomalous µ182W and high 3He/4He in the mantle","authors":"Lori N Willhite, Valerie A Finlayson, Richard J Walker","doi":"10.1016/j.epsl.2024.118795","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118795","url":null,"abstract":"<div><p>Highly siderophile element abundances and <sup>182</sup>W/<sup>184</sup>W and <sup>187</sup>Os/<sup>188</sup>Os were determined for a suite of Mauna Kea lavas from the Hawaiian Scientific Drilling Project phase 2 drill core. The new analyses, combined with previous measurements, compose the largest database for µ<sup>182</sup>W (the parts-per-million deviation of <sup>182</sup>W/<sup>184</sup>W from a terrestrial standard) for a single volcano (<em>n</em> <em>=</em> 16). Although most lavas analyzed are characterized by negative µ<sup>182</sup>W values, lavas with values similar to the modern bulk silicate Earth are found throughout the entire stratigraphic column. This suggests that components with normal µ<sup>182</sup>W are collocated with components that host µ<sup>182</sup>W deficits in the plume. Negative µ<sup>182</sup>W values are associated with elevated <sup>3</sup>He/<sup>4</sup>He, as well as elevated Ti and Nb. These correlations may link µ<sup>182</sup>W anomalies to ancient deep mantle crystal-liquid fractionation processes. Consistent with previously measured <sup>3</sup>He/<sup>4</sup>He (R/R<sub>A</sub>) in the drill core, the magnitude of negative µ<sup>182</sup>W values was greatest when Mauna Kea was close to the plume axis then generally decreased over the ∼400 kyr captured by the stratigraphic section. The component with anomalous µ<sup>182</sup>W was either concentrated near the plume axis, or was more effectively sampled by melting near the plume axis where the temperature excess was greatest, suggesting it was less fusible than the dominant plume components. The process leading to the generation of a mantle component with a negative µ<sup>182</sup>W anomaly could either be related to some form of core-mantle isotopic equilibration, or early-Earth fractionation within the silicate Earth. At present each possibility remains viable.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244463","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-05-31DOI: 10.1016/j.epsl.2024.118801
Peter Köhler , Luke C. Skinner , Florian Adolphi
Carbon cycle models used to calculate the marine reservoir age of the non-polar surface ocean (called Marine20) out of IntCal20, the compilation of atmospheric C, have so far neglected a key aspect of the millennial-scale variability connected with the thermal bipolar seesaw: changes in the strength of the Atlantic meridional overturning circulation (AMOC) related to Dansgaard/Oeschger and Heinrich events. Here we implement such AMOC changes in the carbon cycle box model BICYCLE-SE to investigate how model performance over the last 55 kyr is affected, in particular with respect to available 14C and CO2 data. Constraints from deep ocean 14C data suggest that the AMOC in the model during Heinrich stadial 1 needs to be highly reduced or even completely shutdown. Ocean circulation and sea ice coverage combined are the processes that almost completely explain the simulated changes in deep ocean 14C age, and these are also responsible for a glacial drawdown of ∼60 ppm of atmospheric CO2. We find that the implementation of abrupt reductions in AMOC during Greenland stadials in the model setup that was previously used for the calculation of Marine20 leads to differences of less than ±100 14C yrs. The representation of AMOC changes therefore appears to be of minor importance for deriving non-polar mean ocean radiocarbon calibration products such as Marine20, where atmospheric carbon cycle variables are forced by reconstructions. However, simulated atmospheric CO2 exhibits minima during AMOC reductions in Heinrich stadials, in disagreement with ice core data. This mismatch supports previous suggestions that millennial-scale changes in CO2 were probably not driven directly by the AMOC, but rather by biological and physical processes in the Southern Ocean and by contributions from variable land carbon storage.
{"title":"Simulated radiocarbon cycle revisited by considering the bipolar seesaw and benthic 14C data","authors":"Peter Köhler , Luke C. Skinner , Florian Adolphi","doi":"10.1016/j.epsl.2024.118801","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118801","url":null,"abstract":"<div><p>Carbon cycle models used to calculate the marine reservoir age of the non-polar surface ocean (called Marine20) out of IntCal20, the compilation of atmospheric <span><math><msup><mrow><mi>Δ</mi></mrow><mrow><mn>14</mn></mrow></msup></math></span>C, have so far neglected a key aspect of the millennial-scale variability connected with the thermal bipolar seesaw: changes in the strength of the Atlantic meridional overturning circulation (AMOC) related to Dansgaard/Oeschger and Heinrich events. Here we implement such AMOC changes in the carbon cycle box model BICYCLE-SE to investigate how model performance over the last 55 kyr is affected, in particular with respect to available <sup>14</sup>C and CO<sub>2</sub> data. Constraints from deep ocean <sup>14</sup>C data suggest that the AMOC in the model during Heinrich stadial 1 needs to be highly reduced or even completely shutdown. Ocean circulation and sea ice coverage combined are the processes that almost completely explain the simulated changes in deep ocean <sup>14</sup>C age, and these are also responsible for a glacial drawdown of ∼60 ppm of atmospheric CO<sub>2</sub>. We find that the implementation of abrupt reductions in AMOC during Greenland stadials in the model setup that was previously used for the calculation of Marine20 leads to differences of less than ±100 <sup>14</sup>C yrs. The representation of AMOC changes therefore appears to be of minor importance for deriving non-polar mean ocean radiocarbon calibration products such as Marine20, where atmospheric carbon cycle variables are forced by reconstructions. However, simulated atmospheric CO<sub>2</sub> exhibits minima during AMOC reductions in Heinrich stadials, in disagreement with ice core data. This mismatch supports previous suggestions that millennial-scale changes in CO<sub>2</sub> were probably not driven directly by the AMOC, but rather by biological and physical processes in the Southern Ocean and by contributions from variable land carbon storage.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24002346/pdfft?md5=02e7f6decab49d3433b42b2126685fb7&pid=1-s2.0-S0012821X24002346-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244461","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-05-30DOI: 10.1016/j.epsl.2024.118798
J.A. Parera-Portell , F.d.L. Mancilla , J. Morales , J. Díaz
The mantle transition zone (MTZ) beneath Iberia and NW Maghreb has been precisely mapped using more than 56000 high-quality P-wave receiver functions calculated from the data collected by permanent seismic networks and multiple temporary deployments in the region. Three-dimensional depth migration using both regional and global tomographic models has allowed us to obtain robust and continuous measurements of the MTZ thickness and the depth of the 410 and 660 discontinuities. We found the MTZ thickened by as much as in the Mediterranean coast due to the effect of the cold Gibraltar-Alboran and Alpine-Tethys slabs. Coinciding with expected water-enriched MTZ areas near the subducted slabs there is evidence for partial melting atop the 410 in at least three low-velocity layers (LVL). Partial melting is also likely in a LVL in the uppermost lower mantle under the Alpine-Tethys slab, while we attribute other intra-MTZ LVL to increased mineralogical heterogeneity. We link a thinning of the MTZ at the rear of the Gibraltar-Alboran slab to mantle upwelling, and a band of depressed 410 along its southern boundary as an area of hot toroidal flow. A discontinuous region of depressed 410 following the Atlas Mountains also supports mantle upwelling beneath this range. Areas with LVL atop a depressed 410 discontinuity correlate well with active intraplate volcanism, suggesting a possible MTZ source. We also found that deep-focus seismicity occurs where the 660 discontinuity starts to deepen at the westernmost edge of the Gibraltar-Alboran slab.
{"title":"High-resolution mapping of the mantle transition zone and its interaction with subducted slabs in the Ibero-Maghrebian region","authors":"J.A. Parera-Portell , F.d.L. Mancilla , J. Morales , J. Díaz","doi":"10.1016/j.epsl.2024.118798","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118798","url":null,"abstract":"<div><p>The mantle transition zone (MTZ) beneath Iberia and NW Maghreb has been precisely mapped using more than 56000 high-quality P-wave receiver functions calculated from the data collected by permanent seismic networks and multiple temporary deployments in the region. Three-dimensional depth migration using both regional and global tomographic models has allowed us to obtain robust and continuous measurements of the MTZ thickness and the depth of the 410 and 660 discontinuities. We found the MTZ thickened by as much as <figure><img></figure> in the Mediterranean coast due to the effect of the cold Gibraltar-Alboran and Alpine-Tethys slabs. Coinciding with expected water-enriched MTZ areas near the subducted slabs there is evidence for partial melting atop the 410 in at least three low-velocity layers (LVL). Partial melting is also likely in a LVL in the uppermost lower mantle under the Alpine-Tethys slab, while we attribute other intra-MTZ LVL to increased mineralogical heterogeneity. We link a thinning of the MTZ at the rear of the Gibraltar-Alboran slab to mantle upwelling, and a band of depressed 410 along its southern boundary as an area of hot toroidal flow. A discontinuous region of depressed 410 following the Atlas Mountains also supports mantle upwelling beneath this range. Areas with LVL atop a depressed 410 discontinuity correlate well with active intraplate volcanism, suggesting a possible MTZ source. We also found that deep-focus seismicity occurs where the 660 discontinuity starts to deepen at the westernmost edge of the Gibraltar-Alboran slab.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24002310/pdfft?md5=dc10956e0279ba2a777a561180ea2c41&pid=1-s2.0-S0012821X24002310-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244464","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-05-29DOI: 10.1016/j.epsl.2024.118794
Samuel Angiboust , Paraskevi Io Ioannidi , Iskander Muldashev
A paradox exists between the great number of intermediate-depth earthquakes occurring along active subduction interfaces worldwide and the extreme scarcity of paleo-seismic events recorded in exhumed metasediments from ancient subducted slabs. Recrystallization associated with exhumation-related overprinting generally contributes to the nearly-complete erasing of markers of unstable slip events in metamorphic rocks. We herein focus on a sample from an ancient deep thrust from a Cretaceous High-Pressure paleo-accretionary complex in Chilean Patagonia. A representative, moderately foliated micaschist exhibits broken garnet crystals that host a dense network of healed micro-fractures. While garnet fragments appear thoroughly disaggregated along the main foliation, the rock matrix that completely recrystallized has lost the record of brittle deformation. We employ a 2D visco-elasto-plastic numerical modelling approach in order to investigate the mechanical conditions that enable the fracturing of isolated garnet grains in a relatively weak matrix. The rupture of these stiff grains is achieved in our models at strain rates faster than 10−10 /s to 10−12 /s for elevated pore fluid pressures (80 to 99 % of the lithostatic value, respectively). Since high pore fluid pressures prevail in deep subduction interface settings, it is suggested that the rupture of these garnet crystals occurred through cataclastic deformation via (transient) slip rate acceleration, perhaps as a consequence of localized slip associated with slow to conventional earthquakes. Upon slip rate deceleration, viscous disaggregation of the broken garnet clasts occurred along with the erasing of the matrix cataclastic fabric.
{"title":"Garnet fracturing reveals ancient unstable slip events hosted in plate interface metasediments","authors":"Samuel Angiboust , Paraskevi Io Ioannidi , Iskander Muldashev","doi":"10.1016/j.epsl.2024.118794","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118794","url":null,"abstract":"<div><p>A paradox exists between the great number of intermediate-depth earthquakes occurring along active subduction interfaces worldwide and the extreme scarcity of paleo-seismic events recorded in exhumed metasediments from ancient subducted slabs. Recrystallization associated with exhumation-related overprinting generally contributes to the nearly-complete erasing of markers of unstable slip events in metamorphic rocks. We herein focus on a sample from an ancient deep thrust from a Cretaceous High-Pressure paleo-accretionary complex in Chilean Patagonia. A representative, moderately foliated micaschist exhibits broken garnet crystals that host a dense network of healed micro-fractures. While garnet fragments appear thoroughly disaggregated along the main foliation, the rock matrix that completely recrystallized has lost the record of brittle deformation. We employ a 2D visco-elasto-plastic numerical modelling approach in order to investigate the mechanical conditions that enable the fracturing of isolated garnet grains in a relatively weak matrix. The rupture of these stiff grains is achieved in our models at strain rates faster than 10<sup>−10</sup> /s to 10<sup>−12</sup> /s for elevated pore fluid pressures (80 to 99 % of the lithostatic value, respectively). Since high pore fluid pressures prevail in deep subduction interface settings, it is suggested that the rupture of these garnet crystals occurred through cataclastic deformation via (transient) slip rate acceleration, perhaps as a consequence of localized slip associated with slow to conventional earthquakes. Upon slip rate deceleration, viscous disaggregation of the broken garnet clasts occurred along with the erasing of the matrix cataclastic fabric.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24002279/pdfft?md5=53ba3d5d4823cfd61ea847c40560ef1d&pid=1-s2.0-S0012821X24002279-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244458","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-05-29DOI: 10.1016/j.epsl.2024.118779
Birger Rasmussen , Jian-wei Zi , Andrey Bekker
Major environmental changes during the Great Oxidation Episode are recorded in 2.5–2.0 Ga sedimentary successions across the world. In North America, the Huronian Supergroup is the most complete succession, preserving three glacial intervals with possible equivalents in the Lake Superior region. The third and youngest glacial unit (Gowganda Formation, Cobalt Group) of the Huronian Supergroup has been correlated with glacial deposits of the Chocolay Group in the Lake Superior region. Here we present the results of in situ SHRIMP U-Pb geochronology of zircon in thin tuff layers from the basal Wewe Slate, Chocolay Group, from drill-core in the Marquette Range area, upper Michigan. Zircon U-Pb data from two intervals ∼16 m apart yielded weighted mean 207Pb/206Pb ages of 2174 ± 9 Ma and 2172 ± 6 Ma, which define the depositional ages of the Wewe Slate and the conformably underlying Kona Dolomite. In the Huronian Supergroup, SHRIMP U-Pb dating of zircon in a tuff bed from the Gordon Lake Formation, Cobalt Group, yielded a weighted mean 207Pb/206Pb age of 2318 ± 8 Ma, interpreted to be the age of deposition. Re-examination of previously published U-Pb zircon data from the Gordon Lake Formation (Rasmussen et al., 2013) yields a slightly older weighted mean 207Pb/206Pb age of 2310 ± 5 Ma. Our results show that the Wewe Slate, and probably the Kona Dolomite, are ∼150 million years younger than the ∼2.32 Ga Gordon Lake Formation, with which they were previously correlated. Our results suggest that the glacial deposits in the Lake Superior region represent a fourth Paleoproterozoic glaciation in North America, whose age range corresponds to that of the fourth and final early Paleoproterozoic glaciation in southern Africa. We propose that the Chocolay Group was deposited after the Huronian Supergroup and is broadly coeval with the rift-passive margin succession (Cycle 1 of the Kaniapiskau Supergroup) in the New Québec Orogen, Canada, whose deposition has been linked to rifting and breakup along the eastern margin of the Superior craton at ∼2.22 Ga.
{"title":"New U-Pb zircon tuff ages and revised stratigraphic correlations in the Superior craton during the Great Oxidation Episode","authors":"Birger Rasmussen , Jian-wei Zi , Andrey Bekker","doi":"10.1016/j.epsl.2024.118779","DOIUrl":"10.1016/j.epsl.2024.118779","url":null,"abstract":"<div><p>Major environmental changes during the Great Oxidation Episode are recorded in 2.5–2.0 Ga sedimentary successions across the world. In North America, the Huronian Supergroup is the most complete succession, preserving three glacial intervals with possible equivalents in the Lake Superior region. The third and youngest glacial unit (Gowganda Formation, Cobalt Group) of the Huronian Supergroup has been correlated with glacial deposits of the Chocolay Group in the Lake Superior region. Here we present the results of in situ SHRIMP U-Pb geochronology of zircon in thin tuff layers from the basal Wewe Slate, Chocolay Group, from drill-core in the Marquette Range area, upper Michigan. Zircon U-Pb data from two intervals ∼16 m apart yielded weighted mean <sup>207</sup>Pb/<sup>206</sup>Pb ages of 2174 ± 9 Ma and 2172 ± 6 Ma, which define the depositional ages of the Wewe Slate and the conformably underlying Kona Dolomite. In the Huronian Supergroup, SHRIMP U-Pb dating of zircon in a tuff bed from the Gordon Lake Formation, Cobalt Group, yielded a weighted mean <sup>207</sup>Pb/<sup>206</sup>Pb age of 2318 ± 8 Ma, interpreted to be the age of deposition. <em>Re</em>-examination of previously published U-Pb zircon data from the Gordon Lake Formation (Rasmussen et al., 2013) yields a slightly older weighted mean <sup>207</sup>Pb/<sup>206</sup>Pb age of 2310 ± 5 Ma. Our results show that the Wewe Slate, and probably the Kona Dolomite, are ∼150 million years younger than the ∼2.32 Ga Gordon Lake Formation, with which they were previously correlated. Our results suggest that the glacial deposits in the Lake Superior region represent a fourth Paleoproterozoic glaciation in North America, whose age range corresponds to that of the fourth and final early Paleoproterozoic glaciation in southern Africa. We propose that the Chocolay Group was deposited after the Huronian Supergroup and is broadly coeval with the rift-passive margin succession (Cycle 1 of the Kaniapiskau Supergroup) in the New Québec Orogen, Canada, whose deposition has been linked to rifting and breakup along the eastern margin of the Superior craton at ∼2.22 Ga.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0012821X24002127/pdfft?md5=a0f986f587dfcf98cb9caa7bcf1ad2e7&pid=1-s2.0-S0012821X24002127-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141172653","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-05-29DOI: 10.1016/j.epsl.2024.118793
Yujian Wang, Jingao Liu
Geochemically distinguishing the products of melt depletion from refertilization and constraining the timing of such mantle-melt interactions in the subcontinental lithosphere (SCLM) remain outstanding research issues. Here we utilize alphaMELTS thermodynamic modeling of both partial melting and refertilization to revisit the origins of the Lherz mantle rocks from the Pyrenean orogenic mantle massifs. Thermodynamic modeling reveals subtle but critical differences between refertilization (i.e., elevated TiO2/Al2O3 and higher HREE in both whole rocks and clinopyroxenes) and partial melting processes. Harzburgites from the Lherz massif comprise predominantly pristine residues of partial melting and subordinately refertilized products with low melt/rock ratios. The Lherz lherzolites that show close spatial associations with olivine websterite layers represent secondary rocks, derived from refertilization involving the pristine, refractory harzburgites and upwelling N-MORB-like melts. Lherzolites with no intimate spatial association to olivine websterites are open to an additional origin, i.e., via stationary cooling at the base of lithosphere after moderate adiabatic upwelling of asthenospheric mantle. The Re-depletion model ages (TRD) of refractory harzburgites yield a systematic peak melting age of 2.0 Ga. We have developed a novel solution for constraining the relatively ancient age of refertilization (∼ 1.5 – 2.0 Ga) through the approach of an adapted percolation model to assess the behavior of Re-Os and Lu-Hf isotopic systems. Together with a comprehensive dataset of global on- and off-cratonic SCLM, this study has successfully distinguished silicate-melt induced refertilization from partial melting on elemental level, demonstrated the contribution of different melting mechanisms to the distinctive SCLM compositional evolution in the history, and highlighted how profound a role that refertilization has played on the variations of geochemical buoyancy and mechanical robustness and eventually on the stability and longevity of the ancient SCLM.
{"title":"Thermodynamic modelling and chronometric dating of melting and refertilization in the subcontinental lithospheric mantle","authors":"Yujian Wang, Jingao Liu","doi":"10.1016/j.epsl.2024.118793","DOIUrl":"https://doi.org/10.1016/j.epsl.2024.118793","url":null,"abstract":"<div><p>Geochemically distinguishing the products of melt depletion from refertilization and constraining the timing of such mantle-melt interactions in the subcontinental lithosphere (SCLM) remain outstanding research issues. Here we utilize alphaMELTS thermodynamic modeling of both partial melting and refertilization to revisit the origins of the Lherz mantle rocks from the Pyrenean orogenic mantle massifs. Thermodynamic modeling reveals subtle but critical differences between refertilization (i.e., elevated TiO<sub>2</sub>/Al<sub>2</sub>O<sub>3</sub> and higher HREE in both whole rocks and clinopyroxenes) and partial melting processes. Harzburgites from the Lherz massif comprise predominantly pristine residues of partial melting and subordinately refertilized products with low melt/rock ratios. The Lherz lherzolites that show close spatial associations with olivine websterite layers represent secondary rocks, derived from refertilization involving the pristine, refractory harzburgites and upwelling N-MORB-like melts. Lherzolites with no intimate spatial association to olivine websterites are open to an additional origin, i.e., via stationary cooling at the base of lithosphere after moderate adiabatic upwelling of asthenospheric mantle. The Re-depletion model ages (T<sub>RD</sub>) of refractory harzburgites yield a systematic peak melting age of 2.0 Ga. We have developed a novel solution for constraining the relatively ancient age of refertilization (∼ 1.5 – 2.0 Ga) through the approach of an adapted percolation model to assess the behavior of Re-Os and Lu-Hf isotopic systems. Together with a comprehensive dataset of global on- and off-cratonic SCLM, this study has successfully distinguished silicate-melt induced refertilization from partial melting on elemental level, demonstrated the contribution of different melting mechanisms to the distinctive SCLM compositional evolution in the history, and highlighted how profound a role that refertilization has played on the variations of geochemical buoyancy and mechanical robustness and eventually on the stability and longevity of the ancient SCLM.</p></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":null,"pages":null},"PeriodicalIF":5.3,"publicationDate":"2024-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141244460","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}