Pub Date : 2024-11-16DOI: 10.1016/j.epsl.2024.119095
Chris E. Mohn , Razvan Caracas , Clinton P. Conrad
Proton self diffusion coefficients for bridgmanite at lower mantle conditions are calculated from ab initio molecular dynamics simulations. We find that the proton self diffusion coefficient, is nearly constant ∼ 10−8 m2 s−1 along the lower mantle geotherm but increases by nearly one order of magnitude from ∼10−10 m2 s−1 to ∼ 10−9 m2 s−1 along a cold slab geotherm to about 1800 km depth. These rates imply that the proton diffusion length scale is less than 10 km in lower mantle peridotite in the 150-200 million years timescale for slab material to sink through the lower mantle. Cold wet slabs probably lose less than one percent of their total water content to the ambient mantle on their journey through the lower mantle, indicating that recycled water is far from homogeneously distributed since slab delivery is highly heterogeneous. We estimate that 0.1 to 0.3 ocean masses (<100 ppm wt%) of recycled water may be currently stored in slab remnant materials within the lower mantle. This water is likely not entrained by plumes but is instead captured by background mantle flow before returning to the mid-ocean ridges. By contrast, deep-rooted mantle plumes may entrain materials containing primordial-like water from the lowermost mantle or the core, and preserve these anomalies in fairly small-scale heterogeneities. Over the age of the Earth, the proton diffusion length scale is a few tens of km, which places constraints on the size of possible primordial water reservoirs isolated from convective mixing, and indicates little flux of water across the core-mantle boundary.
{"title":"Lower mantle water distribution from ab initio proton diffusivity in bridgmanite","authors":"Chris E. Mohn , Razvan Caracas , Clinton P. Conrad","doi":"10.1016/j.epsl.2024.119095","DOIUrl":"10.1016/j.epsl.2024.119095","url":null,"abstract":"<div><div>Proton self diffusion coefficients for bridgmanite at lower mantle conditions are calculated from <em>ab initio</em> molecular dynamics simulations. We find that the proton self diffusion coefficient, <span><math><msup><mrow><mi>D</mi></mrow><mrow><mtext>self</mtext></mrow></msup></math></span> is nearly constant ∼ 10<sup>−8</sup> m<sup>2</sup> <!-->s<sup>−1</sup> along the lower mantle geotherm but increases by nearly one order of magnitude from ∼10<sup>−10</sup> m<sup>2</sup> <!-->s<sup>−1</sup> to ∼ 10<sup>−9</sup> m<sup>2</sup> <!-->s<sup>−1</sup> along a cold slab geotherm to about 1800 km depth. These rates imply that the proton diffusion length scale is less than 10 km in lower mantle peridotite in the 150-200 million years timescale for slab material to sink through the lower mantle. Cold wet slabs probably lose less than one percent of their total water content to the ambient mantle on their journey through the lower mantle, indicating that recycled water is far from homogeneously distributed since slab delivery is highly heterogeneous. We estimate that 0.1 to 0.3 ocean masses (<100 ppm wt%) of recycled water may be currently stored in slab remnant materials within the lower mantle. This water is likely not entrained by plumes but is instead captured by background mantle flow before returning to the mid-ocean ridges. By contrast, deep-rooted mantle plumes may entrain materials containing primordial-like water from the lowermost mantle or the core, and preserve these anomalies in fairly small-scale heterogeneities. Over the age of the Earth, the proton diffusion length scale is a few tens of km, which places constraints on the size of possible primordial water reservoirs isolated from convective mixing, and indicates little flux of water across the core-mantle boundary.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"649 ","pages":"Article 119095"},"PeriodicalIF":4.8,"publicationDate":"2024-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142658541","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-11-15DOI: 10.1016/j.epsl.2024.119112
Francis Nimmo , Thorsten Kleine , Alessandro Morbidelli , David Nesvorny
The nucleosynthetic isotope signatures of meteorites and the bulk silicate Earth (BSE) indicate that Earth consists of a mixture of “carbonaceous” (CC) and “non-carbonaceous” (NC) materials. We show that the fraction of CC material recorded in the isotopic composition of the BSE varies for different elements, and depends on the element's tendency to partition into metal and its volatility. The observed behavior indicates that the majority of material accreted to the Earth was NC-dominated, but that CC-dominated material enriched in moderately volatile elements by a factor of ∼10 was delivered during the last ∼2–10% of Earth's accretion. The late delivery of CC material to Earth contrasts with dynamical evidence for the early implantation of CC objects into the inner solar system during the growth and migration of the giant planets. This, together with the NC-dominated nature of both Earth's late veneer and bulk Mars, suggests that material scattered inwards had the bulk of its mass concentrated in a few, large CC embryos rather than in smaller planetesimals. We propose that Earth accreted a few of these CC embryos while Mars did not, and that at least one of the CC embryos impacted Earth relatively late (when accretion was 90–98% complete). This scenario is consistent with the subsequent Moon-forming impact of a large NC body, as long as this impact did not re-homogenize the entire Earth's mantle.
陨石和大块硅酸盐地球(BSE)的核合成同位素特征表明,地球由 "碳质"(CC)和 "非碳质"(NC)物质混合组成。我们的研究表明,BSE 的同位素组成中记录的 CC 物质的比例因不同元素而异,并取决于该元素分化为金属的倾向及其挥发性。观察到的行为表明,地球吸积的大部分物质以NC为主,但以CC为主的物质在地球吸积的最后∼2-10%期间富含中度挥发性元素,富集系数为∼10。CC物质很晚才被输送到地球,这与巨行星生长和迁移过程中CC物体很早就被植入内太阳系的动力学证据形成了鲜明对比。这一点,再加上地球晚期的贴面和火星主体都以NC为主的性质,表明向内散布的物质的大部分质量都集中在少数大型CC胚胎中,而不是较小的行星碎片中。我们提出,地球吸积了其中几个 CC 胚胎,而火星没有,至少有一个 CC 胚胎在相对较晚的时候(吸积完成 90-98% 时)撞击了地球。只要这种撞击没有使整个地球的地幔重新同质化,那么这种情况就与随后一个大型NC天体撞击月球形成的情况是一致的。
{"title":"Mechanisms and timing of carbonaceous chondrite delivery to the Earth","authors":"Francis Nimmo , Thorsten Kleine , Alessandro Morbidelli , David Nesvorny","doi":"10.1016/j.epsl.2024.119112","DOIUrl":"10.1016/j.epsl.2024.119112","url":null,"abstract":"<div><div>The nucleosynthetic isotope signatures of meteorites and the bulk silicate Earth (BSE) indicate that Earth consists of a mixture of “carbonaceous” (CC) and “non-carbonaceous” (NC) materials. We show that the fraction of CC material recorded in the isotopic composition of the BSE varies for different elements, and depends on the element's tendency to partition into metal and its volatility. The observed behavior indicates that the majority of material accreted to the Earth was NC-dominated, but that CC-dominated material enriched in moderately volatile elements by a factor of ∼10 was delivered during the last ∼2–10% of Earth's accretion. The late delivery of CC material to Earth contrasts with dynamical evidence for the early implantation of CC objects into the inner solar system during the growth and migration of the giant planets. This, together with the NC-dominated nature of both Earth's late veneer and bulk Mars, suggests that material scattered inwards had the bulk of its mass concentrated in a few, large CC embryos rather than in smaller planetesimals. We propose that Earth accreted a few of these CC embryos while Mars did not, and that at least one of the CC embryos impacted Earth relatively late (when accretion was 90–98% complete). This scenario is consistent with the subsequent Moon-forming impact of a large NC body, as long as this impact did not re-homogenize the entire Earth's mantle.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119112"},"PeriodicalIF":4.8,"publicationDate":"2024-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660316","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-11-10DOI: 10.1016/j.epsl.2024.119110
Simon J. Lloyd , Andrew J. Biggin , Henry Halls , Steve Denyszyn
The Paleoproterozoic era is the longest in Earth's history, with significant changes hypothesised to have occurred in the deep Earth's physical and chemical conditions at this time. It has been suggested that the paleomagnetic field became weaker at this time (∼2.4 Ga) and remained weak for the next billion years. Paleomagnetism is intrinsically linked to, and is able to inform on, ancient deep Earth processes; a weak dipole strength sustained over this time period may have implications for both core and mantle evolution.
We test this hypothesis here in a two-fold approach: (1) A paleointensity study on the widespread ca. 1.6 Ga diabase/dolerite Melville Bugt dyke swarm. The swarm extends along the west coast of Greenland for more than 1000 km and intruded over ∼13 million years, capturing polarity reversals of Earth's magnetic field. (2) A detailed statistical analysis on the long-term trend in average dipole moment from an improved paleointensity dataset (PINT.org) that has recently undergone a major update.
Five of the Greenland dykes produce paleointensity results ranging from 1.4 μT to 5.1 μT (virtual dipole moment range 0.3–1.2 × 1022 Am2) during the mid-point of this extended period of ‘dipole low’. Our statistical study robustly confirms that this one-billion-year period was indeed associated with an anomalously weak dipole moment (2.7 × 1022 Am2) relative to 500-million-year intervals before and after, which were almost twice as strong. Sampling of more geographically diverse rocks from this time is needed to yield a clear picture of the long-term time evolution of the dipole moment.
{"title":"Weak paleointensities from 1.6 Ga Greenland dykes: Further evidence for a billion-year period of paleomagnetic dipole low during the Paleoproterozoic","authors":"Simon J. Lloyd , Andrew J. Biggin , Henry Halls , Steve Denyszyn","doi":"10.1016/j.epsl.2024.119110","DOIUrl":"10.1016/j.epsl.2024.119110","url":null,"abstract":"<div><div>The Paleoproterozoic era is the longest in Earth's history, with significant changes hypothesised to have occurred in the deep Earth's physical and chemical conditions at this time. It has been suggested that the paleomagnetic field became weaker at this time (∼2.4 Ga) and remained weak for the next billion years. Paleomagnetism is intrinsically linked to, and is able to inform on, ancient deep Earth processes; a weak dipole strength sustained over this time period may have implications for both core and mantle evolution.</div><div>We test this hypothesis here in a two-fold approach: (1) A paleointensity study on the widespread ca. 1.6 Ga diabase/dolerite Melville Bugt dyke swarm. The swarm extends along the west coast of Greenland for more than 1000 km and intruded over ∼13 million years, capturing polarity reversals of Earth's magnetic field. (2) A detailed statistical analysis on the long-term trend in average dipole moment from an improved paleointensity dataset (PINT.org) that has recently undergone a major update.</div><div>Five of the Greenland dykes produce paleointensity results ranging from 1.4 μT to 5.1 μT (virtual dipole moment range 0.3–1.2 × 10<sup>22</sup> Am<sup>2</sup>) during the mid-point of this extended period of ‘dipole low’. Our statistical study robustly confirms that this one-billion-year period was indeed associated with an anomalously weak dipole moment (2.7 × 10<sup>22</sup> Am<sup>2</sup>) relative to 500-million-year intervals before and after, which were almost twice as strong. Sampling of more geographically diverse rocks from this time is needed to yield a clear picture of the long-term time evolution of the dipole moment.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119110"},"PeriodicalIF":4.8,"publicationDate":"2024-11-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660315","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-11-09DOI: 10.1016/j.epsl.2024.119105
Jaime Araya Vargas , Jorge Sanhueza , Daniel Díaz , María José Segovia , Francisco Pastén-Araya , Katarzyna Ślęzak
The Patagonian slab window (PSW) is a region of the Southamerican subduction zone where the absence of subducted slabs is interpreted, due to the subduction of the Chile Mid-Ocean Ridge at the Chile Triple Junction. Here we report the results of a long-period magnetotelluric (MT) study conducted in two 300 km-long trench-parallel transects crossing the northern boundary of the PSW in the proximal backarc. We modeled the MT data using 3-D inversion, obtaining an electrical resistivity model of the continental crust and upper mantle up to a depth of ∼150 km. Our model shows a heterogeneous resistivity structure in the uppermost mantle, dominated by resistivities >300 Ωm below the array of sites even within the PSW, and some low-resistivity zones (LRZs, <10 Ωm) mainly at the edge of the array. Using petrophysical models, we estimated the mantle temperature, water content, melt fraction, and viscosity based on obtained resistivity values and a preexistent model of P-wave velocity (Vp) at 50 km and 100 km depth. These estimates suggest that the uppermost mantle within the PSW region is heterogeneous and dominated by high-viscosity blocks, compatible with the continental mantle lithosphere or even subducted slabs. Based on relatively hot and low-viscosity zones estimated in the periphery of LRZs, we interpret the presence of asthenospheric mantle in areas where LRZs coincide with relatively low Vp. According to this interpretation, asthenospheric upwelling in the study area at depths ≤150 km would be localized rather than ubiquitous over the interpreted extent for the PSW. Such localized asthenosphere upwelling processes in the past could explain the scattered distribution of Neogene basaltic lavas in the southern Patagonia backarc. The continental crust exhibits LRZs in the upper and lower crust. Remarkably, ensembles of LRZs at different crustal depths within the presumable area of the PSW were found below the General Carrera Lake, and towards the North Patagonian Icefield, likely indicating the presence of hidden intraplate magmatic and/or hydrothermal systems.
{"title":"Magnetotelluric image of the Patagonian slab window: Constraints on upper mantle physical properties and sources of intraplate magmatism","authors":"Jaime Araya Vargas , Jorge Sanhueza , Daniel Díaz , María José Segovia , Francisco Pastén-Araya , Katarzyna Ślęzak","doi":"10.1016/j.epsl.2024.119105","DOIUrl":"10.1016/j.epsl.2024.119105","url":null,"abstract":"<div><div>The Patagonian slab window (PSW) is a region of the Southamerican subduction zone where the absence of subducted slabs is interpreted, due to the subduction of the Chile Mid-Ocean Ridge at the Chile Triple Junction. Here we report the results of a long-period magnetotelluric (MT) study conducted in two 300 km-long trench-parallel transects crossing the northern boundary of the PSW in the proximal backarc. We modeled the MT data using 3-D inversion, obtaining an electrical resistivity model of the continental crust and upper mantle up to a depth of ∼150 km. Our model shows a heterogeneous resistivity structure in the uppermost mantle, dominated by resistivities >300 Ωm below the array of sites even within the PSW, and some low-resistivity zones (LRZs, <10 Ωm) mainly at the edge of the array. Using petrophysical models, we estimated the mantle temperature, water content, melt fraction, and viscosity based on obtained resistivity values and a preexistent model of P-wave velocity (<em>Vp</em>) at 50 km and 100 km depth. These estimates suggest that the uppermost mantle within the PSW region is heterogeneous and dominated by high-viscosity blocks, compatible with the continental mantle lithosphere or even subducted slabs. Based on relatively hot and low-viscosity zones estimated in the periphery of LRZs, we interpret the presence of asthenospheric mantle in areas where LRZs coincide with relatively low <em>Vp</em>. According to this interpretation, asthenospheric upwelling in the study area at depths ≤150 km would be localized rather than ubiquitous over the interpreted extent for the PSW. Such localized asthenosphere upwelling processes in the past could explain the scattered distribution of Neogene basaltic lavas in the southern Patagonia backarc. The continental crust exhibits LRZs in the upper and lower crust. Remarkably, ensembles of LRZs at different crustal depths within the presumable area of the PSW were found below the General Carrera Lake, and towards the North Patagonian Icefield, likely indicating the presence of hidden intraplate magmatic and/or hydrothermal systems.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119105"},"PeriodicalIF":4.8,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660314","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-11-08DOI: 10.1016/j.epsl.2024.119107
Chengli Liu , Yefei Bai , Thorne Lay , Ping He , Yangmao Wen , Xiaoran Wei , Neng Xiong , Xiong Xiong
A damaging MW 7.5 earthquake struck the western coast of Japan along the Noto Peninsula on January 1, 2024. The initiation of large shallow earthquakes along the Noto Peninsula, particularly above deeper long-duration patchy seismic swarms, presents an unusual seismic phenomenon that warrants in-depth investigation of their interactions. The 2024 earthquake nucleated with an initial low average rupture velocity of 0.5–1.0 km s-1 near the up-dip end of a long-lasting seismic swarm that commenced in November 2020. Analysis of dense seismic, geodetic, and tsunami observations provides good resolution of large shallow slip in the crust below the peninsula and extending offshore to the northeast, revealing a heterogeneous slip distribution characterized by bilateral two-stage rupture expansion during the faulting. Up to 8 m of slip occurred in several patches along ∼150 km of the southeastward-dipping thrust fault, which extends to near the seafloor along the northwest side of the peninsula. Up to 5 m of uplift occurred along the peninsula's northwestern coast. Up-dip fluid migration appears to have weakened the shallow fault prior to failure and influenced the initial slow rupture expansion, highlighting the need to monitor the evolution of worldwide swarms.
{"title":"Shallow crustal rupture in a major MW 7.5 earthquake above a deep crustal seismic swarm along the Noto Peninsula in western Japan","authors":"Chengli Liu , Yefei Bai , Thorne Lay , Ping He , Yangmao Wen , Xiaoran Wei , Neng Xiong , Xiong Xiong","doi":"10.1016/j.epsl.2024.119107","DOIUrl":"10.1016/j.epsl.2024.119107","url":null,"abstract":"<div><div>A damaging <em>M<sub>W</sub></em> 7.5 earthquake struck the western coast of Japan along the Noto Peninsula on January 1, 2024. The initiation of large shallow earthquakes along the Noto Peninsula, particularly above deeper long-duration patchy seismic swarms, presents an unusual seismic phenomenon that warrants in-depth investigation of their interactions. The 2024 earthquake nucleated with an initial low average rupture velocity of 0.5–1.0 km s<sup>-1</sup> near the up-dip end of a long-lasting seismic swarm that commenced in November 2020. Analysis of dense seismic, geodetic, and tsunami observations provides good resolution of large shallow slip in the crust below the peninsula and extending offshore to the northeast, revealing a heterogeneous slip distribution characterized by bilateral two-stage rupture expansion during the faulting. Up to 8 m of slip occurred in several patches along ∼150 km of the southeastward-dipping thrust fault, which extends to near the seafloor along the northwest side of the peninsula. Up to 5 m of uplift occurred along the peninsula's northwestern coast. Up-dip fluid migration appears to have weakened the shallow fault prior to failure and influenced the initial slow rupture expansion, highlighting the need to monitor the evolution of worldwide swarms.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119107"},"PeriodicalIF":4.8,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660312","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-11-08DOI: 10.1016/j.epsl.2024.119109
Line Colin , Chloé Michaut , Stéphane Labrosse , Bernard Bourdon
The energy of the giant impact was large enough to generate an initially fully molten Moon. During the solidification of this lunar magma ocean (LMO), an anorthosite crust formed by flotation of light anorthite crystals. Lunar anorthosites show crystallization ages as young as 4.360 Gyr, suggesting a long-lived LMO or a rather young Moon. Existing models for LMO solidification are for a specific phase diagram based on one compositional model. However, the LMO solidification timescale depends on the lunar bulk composition and on the appearance of anorthite in the crystallization sequence.
Here, we propose a physically robust 1D model for LMO evolution based on a simple anorthite/olivine-pyroxene eutectic phase diagram. Cumulates first settle at the ocean base for about a thousand years. This first stage results in an unstable thermal profile for the cumulates that can lead to their overturn. In the second stage, simultaneous crystallization of anorthite and cumulates leads to the formation of a buoyant lid that considerably slows down LMO cooling.
We explore the impact of an initially hydrated composition, which reduces the stability of plagioclase, of the eutectic position and of the crust thermal conductivity. We show that cumulates overturn may reduce or extend the LMO solidification time depending on its duration. The total LMO solidification timescale ranges between 45 and 250 Myr. Given the most reliable age of 4.360 Gyr for FAN sample 60025, which derives from more than 99% of crystallization, we estimate an age of 4400 to 4560 Myr for the Moon.
{"title":"Thermal evolution of the lunar magma ocean","authors":"Line Colin , Chloé Michaut , Stéphane Labrosse , Bernard Bourdon","doi":"10.1016/j.epsl.2024.119109","DOIUrl":"10.1016/j.epsl.2024.119109","url":null,"abstract":"<div><div>The energy of the giant impact was large enough to generate an initially fully molten Moon. During the solidification of this lunar magma ocean (LMO), an anorthosite crust formed by flotation of light anorthite crystals. Lunar anorthosites show crystallization ages as young as 4.360 Gyr, suggesting a long-lived LMO or a rather young Moon. Existing models for LMO solidification are for a specific phase diagram based on one compositional model. However, the LMO solidification timescale depends on the lunar bulk composition and on the appearance of anorthite in the crystallization sequence.</div><div>Here, we propose a physically robust 1D model for LMO evolution based on a simple anorthite/olivine-pyroxene eutectic phase diagram. Cumulates first settle at the ocean base for about a thousand years. This first stage results in an unstable thermal profile for the cumulates that can lead to their overturn. In the second stage, simultaneous crystallization of anorthite and cumulates leads to the formation of a buoyant lid that considerably slows down LMO cooling.</div><div>We explore the impact of an initially hydrated composition, which reduces the stability of plagioclase, of the eutectic position and of the crust thermal conductivity. We show that cumulates overturn may reduce or extend the LMO solidification time depending on its duration. The total LMO solidification timescale ranges between 45 and 250 Myr. Given the most reliable age of 4.360 Gyr for FAN sample 60025, which derives from more than 99% of crystallization, we estimate an age of 4400 to 4560 Myr for the Moon.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119109"},"PeriodicalIF":4.8,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660311","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-11-08DOI: 10.1016/j.epsl.2024.119099
O. Bernard, F. Costa
The pre- and syn-eruptive magma decompression rate is recognized as a key parameter modulating eruption dynamics, with explosive eruptions being generally associated with much larger decompression rates than effusive ones. Magma decompression rates cannot be directly measured and thus are typically inferred from petrological, geochemical, numerical modelling, and seismic data. Most studies use petrological information of volatile element diffusive equilibration in glass and crystals to infer a single value for the magma ascent rate for a given eruption, even though numerical volcano conduit simulations show that changes of velocity are expected during magma ascent. Here we integrate magma ascent conduit models with diffusion chronometry of volatiles in melt embayments and phenocrysts to obtain a more comprehensive understanding of magma ascent rates. We find that incorporating a more realistic boundary condition that depends on the magma ascent path with variable velocities gives time estimates that can be up to a factor of 7 longer than from the standard assumption of constant magma ascent rate. Therefore, previous magma ascent rates from diffusion chronometry of volatiles in crystals and melts with a fixed boundary condition may be significantly overestimated. Overall, we show that coupling of magma ascent models with diffusion chronometry can provide more robust inferences of magma ascent and thus improve the understanding of the role of this parameter into the explosive and effusive eruption controls.
{"title":"Coupling magma ascent models with volatile diffusion chronometry","authors":"O. Bernard, F. Costa","doi":"10.1016/j.epsl.2024.119099","DOIUrl":"10.1016/j.epsl.2024.119099","url":null,"abstract":"<div><div>The pre- and syn-eruptive magma decompression rate is recognized as a key parameter modulating eruption dynamics, with explosive eruptions being generally associated with much larger decompression rates than effusive ones. Magma decompression rates cannot be directly measured and thus are typically inferred from petrological, geochemical, numerical modelling, and seismic data. Most studies use petrological information of volatile element diffusive equilibration in glass and crystals to infer a single value for the magma ascent rate for a given eruption, even though numerical volcano conduit simulations show that changes of velocity are expected during magma ascent. Here we integrate magma ascent conduit models with diffusion chronometry of volatiles in melt embayments and phenocrysts to obtain a more comprehensive understanding of magma ascent rates. We find that incorporating a more realistic boundary condition that depends on the magma ascent path with variable velocities gives time estimates that can be up to a factor of 7 longer than from the standard assumption of constant magma ascent rate. Therefore, previous magma ascent rates from diffusion chronometry of volatiles in crystals and melts with a fixed boundary condition may be significantly overestimated. Overall, we show that coupling of magma ascent models with diffusion chronometry can provide more robust inferences of magma ascent and thus improve the understanding of the role of this parameter into the explosive and effusive eruption controls.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119099"},"PeriodicalIF":4.8,"publicationDate":"2024-11-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660313","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-11-07DOI: 10.1016/j.epsl.2024.119102
Yongliang Bai , Kunpeng Cheng , Dongdong Dong , Xinjian Song , Sanzhong Li , Shiguo Wu , Zhenjie Wang
Mantle plume and mid-ocean ridge (MOR) interactions have a profound impact on geological processes, such as large igneous province eruptions, mid-ocean ridge jumps and plume flow towards the MOR. The Walvis Ridge and Rio Grande Rise provide an excellent opportunity to study these interactions, as they have a record of during- and post-interaction stages. We have estimated the hotspot melt production rate along the Walvis Ridge - Rio Grande Rise. Our findings are as follows: (1) The hotspot melt production rate during plume-MOR interaction is seven times higher than that after the interaction, primarily because a thinner lithospheric coverage over the plume allows for greater decompression melting and facilitates the ascent of melt through the lithospheric mantle to accrete to the crust. (2) The waning of the Tristan-Gough plume is indicated by a decreasing melt production rate and an increase in inter-volcano spacing. (3) The linear ridges on the western side of the Guyot Province likely result from mantle plume flow beneath the lithosphere from the hotspot to the Mid-Atlantic Ridge. (4) Fracture zones enhance hotspot melt production rates.
{"title":"Mantle plume – mid ocean ridge interactions revealed from hotspot melt production rate along the Walvis Ridge-Rio Grande Rise","authors":"Yongliang Bai , Kunpeng Cheng , Dongdong Dong , Xinjian Song , Sanzhong Li , Shiguo Wu , Zhenjie Wang","doi":"10.1016/j.epsl.2024.119102","DOIUrl":"10.1016/j.epsl.2024.119102","url":null,"abstract":"<div><div>Mantle plume and mid-ocean ridge (MOR) interactions have a profound impact on geological processes, such as large igneous province eruptions, mid-ocean ridge jumps and plume flow towards the MOR. The Walvis Ridge and Rio Grande Rise provide an excellent opportunity to study these interactions, as they have a record of during- and post-interaction stages. We have estimated the hotspot melt production rate along the Walvis Ridge - Rio Grande Rise. Our findings are as follows: (1) The hotspot melt production rate during plume-MOR interaction is seven times higher than that after the interaction, primarily because a thinner lithospheric coverage over the plume allows for greater decompression melting and facilitates the ascent of melt through the lithospheric mantle to accrete to the crust. (2) The waning of the Tristan-Gough plume is indicated by a decreasing melt production rate and an increase in inter-volcano spacing. (3) The linear ridges on the western side of the Guyot Province likely result from mantle plume flow beneath the lithosphere from the hotspot to the Mid-Atlantic Ridge. (4) Fracture zones enhance hotspot melt production rates.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119102"},"PeriodicalIF":4.8,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660310","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-11-07DOI: 10.1016/j.epsl.2024.119108
C. Zuccari , F. Mazzarini , E. Tavarnelli , G. Viola , L. Aldega , V. Moretto , R. Xie , G. Musumeci
We report the first absolute deformation ages for the Mykonos Detachment that juxtaposes, in the context of the current Aegean rifting, Miocene siliciclastic deposits in the hanging wall against metabasites and synkinematic granites in the footwall. We identified and characterised 16 brittle structural facies (BSFs) within the detachment fault architecture through fieldwork, optical microscopy, and X-ray diffraction. BSFs dating by K-Ar of synkinematic mixed layer illite-smectite shows that the preserved BSFs formed during repeated slip events, thus constraining protracted faulting between 13.5 and 6.5 Ma. Dating, structural and mineralogical characterisation allowed for the time-constrained evaluation of the activation/de-activation of the involved deformation mechanisms and of the processes, including (i) fault zone nucleation, (ii) deformation partitioning, (iii) cataclasis, gouge formation and (iv) the final deformation localisation, that govern the evolution of brittle detachments, shaping the local current upper crustal structure. Our results provide new insights into the understanding of wide active and fossil rift systems.
我们首次报告了米克诺斯支脉的绝对变形年龄,在当前爱琴海断裂的背景下,该支脉将悬壁的中新世硅质沉积物与底壁的偏闪长岩和同生花岗岩并列在一起。我们通过实地考察、光学显微镜和 X 射线衍射,确定并描述了剥离断层结构中的 16 个脆性结构面(BSF)。通过同步混合层伊利石-闪长岩的 K-Ar 测定,BSFs 的年代显示,保存下来的 BSFs 是在反复滑动过程中形成的,从而确定了 13.5 至 6.5 Ma 之间的漫长断层。通过年代测定、结构和矿物学特征描述,可以对所涉及的变形机制的激活/去激活过程以及包括(i) 断层区成核、(ii) 变形分区、(iii) 角砾岩、沟槽形成和(iv) 最终变形定位在内的各种过程进行有时间限制的评估,这些过程控制着脆性脱离的演化,塑造了当地当前的上地壳结构。我们的研究结果为了解广阔的活动断裂和化石断裂系统提供了新的视角。
{"title":"How brittle detachments form and evolve through space and time","authors":"C. Zuccari , F. Mazzarini , E. Tavarnelli , G. Viola , L. Aldega , V. Moretto , R. Xie , G. Musumeci","doi":"10.1016/j.epsl.2024.119108","DOIUrl":"10.1016/j.epsl.2024.119108","url":null,"abstract":"<div><div>We report the first absolute deformation ages for the Mykonos Detachment that juxtaposes, in the context of the current Aegean rifting, Miocene siliciclastic deposits in the hanging wall against metabasites and synkinematic granites in the footwall. We identified and characterised 16 brittle structural facies (BSFs) within the detachment fault architecture through fieldwork, optical microscopy, and X-ray diffraction. BSFs dating by K-Ar of synkinematic mixed layer illite-smectite shows that the preserved BSFs formed during repeated slip events, thus constraining protracted faulting between 13.5 and 6.5 Ma. Dating, structural and mineralogical characterisation allowed for the time-constrained evaluation of the activation/de-activation of the involved deformation mechanisms and of the processes, including (i) fault zone nucleation, (ii) deformation partitioning, (iii) cataclasis, gouge formation and (iv) the final deformation localisation, that govern the evolution of brittle detachments, shaping the local current upper crustal structure. Our results provide new insights into the understanding of wide active and fossil rift systems.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119108"},"PeriodicalIF":4.8,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660309","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-11-06DOI: 10.1016/j.epsl.2024.119094
Yunzhe Chen , Peng Wu , Tiago Angelo , Christopher Spencer , Long Li
Subducting sediments contain much more abundant nitrogen (N) than the other slab components. The fate of sedimentary N along subduction directly determines the quantity of surficial N that can be released back to the atmosphere/crust or further delivered to the deep mantle, which has a strong impact on the long-term evolution of the N budgets and isotopic signatures of Earth's major reservoirs. Previous studies on the fate of subducted sedimentary N have been intensively focused on assessing the extent of metamorphic N devolatilization from metasedimentary rocks varying from ultrahigh-pressure eclogites facies (subducted to as deep as 90 km) in cold subduction zones to medium-pressure amphibolite facies in extremely hot subduction zones. However, the fate of sedimentary N after these metamorphic stages, particularly during sediment melting and subsequent magmatic crystallization, has not been examined so far. Here, we present the N data of the Cretaceous peraluminous granitoids in the Oman-United Arab Emirates area, which were formed by melting of subducted pelagic sediments followed by intrusion and crystallization in the lithospheric mantle before the entire lithospheric section was obducted to form the Oman-United Arab Emirates ophiolites. The 27 studied samples show surprisingly low N contents (4.0 - 13.7 ppm; mean: 6.9 ± 2.3 ppm; 1σ) and a narrow δ15N range (mostly in the range of –1.4‰ to +2.2‰; mean: +0.4 ± 1.5‰; 1σ). Data modeling indicates that the N in these granitoids contains little mantle contribution and was predominantly inherited from their source rocks. Compared with potential source rocks, these granitoids display orders of magnitude lower N contents but comparable δ15N values, indicating that the significant N loss during the genesis of these granitoids was associated with little N isotope fractionation. This suggests that magmatic N2 degassing (with minor N isotope fractionation) rather than metamorphic N devolatilization (with much larger N isotope fractionation) is the dominant mechanism driving the N loss from subducted sediments.
俯冲沉积物中含有比其他板块成分丰富得多的氮。沉积氮在俯冲过程中的去向直接决定了表层氮能够释放回大气/地壳或进一步输送到地幔深处的数量,这对地球主要储层的氮预算和同位素特征的长期演化有很大影响。以往对俯冲沉积 N 的归宿的研究主要集中在评估变质 N 从变质岩中脱落的程度,这些变质岩从寒冷俯冲带的超高压夕卡岩面(俯冲深度达 90 千米)到极热俯冲带的中压闪长岩面都有。然而,迄今为止还没有研究过沉积 N 在这些变质阶段之后的去向,特别是在沉积物熔融和随后的岩浆结晶过程中的去向。在此,我们介绍了阿曼-阿联酋地区白垩纪围岩花岗岩的氮数据,这些花岗岩是在整个岩石圈断面被俯冲形成阿曼-阿联酋蛇绿岩之前,由俯冲沉积物熔融后在岩石圈地幔中侵入和结晶形成的。所研究的 27 个样本显示出令人惊讶的低 N 含量(4.0 - 13.7 ppm;平均值:6.9 ± 2.3 ppm;1σ)和狭窄的 δ15N 范围(大部分在-1.4‰至 +2.2‰之间;平均值:+0.4 ± 1.5‰;1σ)。数据建模表明,这些花岗岩中的氮含量很少来自地幔,而主要来自其源岩。与潜在的源岩相比,这些花岗岩中的氮含量低了几个数量级,但δ15N值却相当,这表明这些花岗岩在成因过程中的大量氮损失与很少的氮同位素分馏有关。这表明,岩浆N2脱气(N同位素分馏小)而非变质N脱溶(N同位素分馏大得多)是驱动N从俯冲沉积物中流失的主要机制。
{"title":"Magmatic N2 degassing dominates nitrogen loss during sediment subduction and granitoid genesis","authors":"Yunzhe Chen , Peng Wu , Tiago Angelo , Christopher Spencer , Long Li","doi":"10.1016/j.epsl.2024.119094","DOIUrl":"10.1016/j.epsl.2024.119094","url":null,"abstract":"<div><div>Subducting sediments contain much more abundant nitrogen (N) than the other slab components. The fate of sedimentary N along subduction directly determines the quantity of surficial N that can be released back to the atmosphere/crust or further delivered to the deep mantle, which has a strong impact on the long-term evolution of the N budgets and isotopic signatures of Earth's major reservoirs. Previous studies on the fate of subducted sedimentary N have been intensively focused on assessing the extent of metamorphic N devolatilization from metasedimentary rocks varying from ultrahigh-pressure eclogites facies (subducted to as deep as 90 km) in cold subduction zones to medium-pressure amphibolite facies in extremely hot subduction zones. However, the fate of sedimentary N after these metamorphic stages, particularly during sediment melting and subsequent magmatic crystallization, has not been examined so far. Here, we present the N data of the Cretaceous peraluminous granitoids in the Oman-United Arab Emirates area, which were formed by melting of subducted pelagic sediments followed by intrusion and crystallization in the lithospheric mantle before the entire lithospheric section was obducted to form the Oman-United Arab Emirates ophiolites. The 27 studied samples show surprisingly low N contents (4.0 - 13.7 ppm; mean: 6.9 ± 2.3 ppm; 1σ) and a narrow δ<sup>15</sup>N range (mostly in the range of –1.4‰ to +2.2‰; mean: +0.4 ± 1.5‰; 1σ). Data modeling indicates that the N in these granitoids contains little mantle contribution and was predominantly inherited from their source rocks. Compared with potential source rocks, these granitoids display orders of magnitude lower N contents but comparable δ<sup>15</sup>N values, indicating that the significant N loss during the genesis of these granitoids was associated with little N isotope fractionation. This suggests that magmatic N<sub>2</sub> degassing (with minor N isotope fractionation) rather than metamorphic N devolatilization (with much larger N isotope fractionation) is the dominant mechanism driving the N loss from subducted sediments.</div></div>","PeriodicalId":11481,"journal":{"name":"Earth and Planetary Science Letters","volume":"648 ","pages":"Article 119094"},"PeriodicalIF":4.8,"publicationDate":"2024-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142593868","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}