首页期刊导航|Earth and Planetary Science Letters
期刊信息/Journal information
Earth and Planetary Science Letters
North-Holland
Earth and Planetary Science Letters

North-Holland

0012-821X

Earth and Planetary Science Letters/Journal Earth and Planetary Science LettersSCIEIAHCIISTP
正式出版
收录年代

    In situ oxidation of sulfide minerals supports widespread sulfate reducing bacteria in the deep subsurface of the Witwatersrand Basin (South Africa): Insights from multiple sulfur and oxygen isotopes

    Li, LongWei, SiwenLollar, Barbara SherwoodWing, Boswell...
    12页
    查看更多>>摘要:Dissolved sulfate is a crucial electron acceptor for the subsurface biosphere, particularly for the living microbial ecosystems in the long-isolated (on the order of millions to billions of years) deep subsurface fracture waters in Precambrian cratons, e.g., in the Witwatersrand Basin of the Kaapvaal Craton, South Africa. Aiming to understand the role of sulfate in the sustainability of the subsurface habitat and the spatial extent of the terrestrial subsurface biosphere, we carried out a basin-scale examination on the sources and producing mechanisms of dissolved sulfur (sulfate and sulfide) in the subsurface fracture waters in the Witwatersrand Basin using multiple sulfur isotopes (S-32,S-33,S-34,S-36) and oxygen isotopes (O-16,O-18). Dissolved sulfates in14 fracture water samples collected from 900 to 3,413 meters below the surface at 6 sites show isotopic ranges from -5.3 parts per thousand to +19.4 parts per thousand for delta S-34, -0.40 parts per thousand to +0.50 parts per thousand for Delta S-33, and -1.1 parts per thousand to +10.9 parts per thousand for delta O-18. These isotopic signatures aredistinct to those of the sulfate minerals (e.g., the carbonate-associated sulfate in local Transvaal Supergroup dolomite sequences: delta S-34 =+31.4 parts per thousand to +39.2 parts per thousand; Delta S-33 =-0.01 parts per thousand to +0.16 parts per thousand), but identical to those of the sulfide minerals in the host rocks. This indicates that the dissolved sulfate in the fracture waters were dominantly generated by in situ oxidation of sulfide minerals at basin scale, although mixing of a small amount of surface sulfate at some sites cannot be completely ruled out. The dissolved sulfates inthe less deep and less saline fracture waters are in oxygen isotope disequilibrium with their host waters, a surprising result given that the water residence times are orders of magnitude longer than the time required for oxygen isotope exchange to reach equilibrium. This implies that vigorous in situ sulfate production has occurred after the fracture waters were isolated. In contrast, the dissolved sulfate in the deeper, more saline waters are in apparent oxygen isotope equilibration or close to equilibration with their host waters. This might be attributed to a combined effect of faster oxygen isotope exchange between sulfate and water at higher temperatures, larger extent of sulfate reduction, and/or less efficient sulfate production. The dissolved sulfide in the fracture waters has similar Delta S-33 values but is 3.0 parts per thousand to 26.4 parts per thousand lower in delta S-34 than coexisting sulfate, suggesting that the dissolved sulfide is mostly generated from bacterial sulfate reduction, which is consistent with the widespread existence of sulfate-reducing bacteria down to > 3.4 km below surface in the Witwatersrand Basin. Overall, these novel isotopic results demonstrate that geological processes can provide a steady long-term sulfate source in deep fracture fluids by in situ oxidation of sulfide minerals in their host rocks, and thereby a mechanism to sustain the terrestrial subsurface biosphere, even in deep high-temperature, long-isolated water systems. Thus, sulfate as a terminal electron acceptor is not the limiting factor for the spatial expansion of terrestrial subsurface biosphere. (C) 2021 Elsevier B.V. All rights reserved.

    Sedimentary and crustal structure of the US Gulf Coast revealed by Rayleigh wave and teleseismic P coda data with implications for continent rifting

    Miao, WenpeiNiu, FenglinLi, GuoliangLevander, Alan...
    13页
    查看更多>>摘要:We have developed an S-wave model of the south-central US focusing on the Gulf Coast sedimentary basin and its crust to understand continental rifting and regional tectonics. The model was derived by a joint inversion of Rayleigh wave phase velocities, Z/Hratios and P-coda data. The surface- and body-wave measurements were made, respectively, from ambient noise and teleseismic earthquakes recorded by 215 USArray stations in a rectangular area of 100.-87. west and 28.-37. north. We employed a cross-convolution function and H-.analysis to better constrain sedimentary and Moho structure. We find that the southern edge of the Ouachita fold-and-thrust belt (OFTB) appears as a boundary in measured phase velocities, Z/Hratios, sediment basement depths, Moho depths, average crustal V-s and V-p/V-s ratios. The model shows southeastward thickening of the sedimentary basin, accompanied by thinning of the crystalline crust. The Moho gradient suggests that early rifting between North America and the Yucatan block commenced in a SE direction and involved most of the Pangea crust south of the OFTB boundary (i.e., Gondwana crust). We believe that a high velocity feature in the lowermost crust and upper mantle parallel to the southeast Texas coast was emplaced as a mafic body and is the source of the Houston magnetic anomaly. The seismic structures of the crust and uppermost mantle observed beneath the Mississippi Embayment and the Mississippi Valley Graben are consistent with plume induced Cretaceous uplift of the Mississippi Embayment as North America passed over the Bermuda hotspot. (C) 2021 Elsevier B.V. All rights reserved.

    What comes first: The fault or the ductile shear zone?

    Scholz, Christopher H.Choi, Eunseo
    8页
    查看更多>>摘要:Crustal scale fault zones extend below the brittle-ductile transition as ductile shear zones. Here we address the question of which regime, brittle or ductile, initiates and controls the overall system of shear localization. Observations of crustal scale but low displacement conjugate strike-slip faults show that they are typically nearly orthogonal, as expected from plastic shear criteria. Sub-crustal scale conjugate strike-slip faults, however, have acute dihedral angles in accordance with the Coulomb fracture criterion as Anderson's theory predicts. We modeled the crustal scale system with strain weakening rheologies that follow the Coulomb and von Mises criteria respectively, within the brittle and ductile regimes. We find that when the strain weakening rate in the ductile regime exceeds a critical value the entire system shears in the von Mises mode with orthogonal conjugate shears forming at all depths, in accordance with the observations. There is, as observed, no deflection of fault orientation at the brittle-ductile transition. Anderson's theory of faulting thus breaks down for crustal scale faults. This difference in behavior was already evident in data presented in Anderson's 1951 book, but its significance was not understood at that time. (C) 2021 The Author(s). Published by Elsevier B.V.

    Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure: Implications for stability and seismic signature of Fe-rich melts in the mantle

    Xu, ManJing, ZhichengYu, TonyAlp, E. Ercan...
    12页
    查看更多>>摘要:Iron-rich silicate melts play an important role in the magmatic history of the Earth and the Moon. However, their elastic properties at high pressures, especially the sound velocities, are poorly understood. Here we determined the ultrasonic sound velocities of two silicate melts along the hedenbergite (Hd, CaFeSi2O6) - diopside (Di, CaMgSi2O6) join at high pressure and temperature conditions up to 6GPa and 2329K, using the high-pressure ultrasonic technique combined with synchrotron radiation in a multianvil apparatus. Our results show that Fe can significantly reduce the sound velocity while increasing the density of silicate melts. Comparing the melts of Di, Hd, and a mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd(50)Di(50)), we find that although densities of the Hd-Di melts can be well-described by a linear mixing law empirically at high pressures, sound velocities do not vary with composition linearly. Assuming that the low-velocity regions in Earth's upper mantle are mainly due to the presence of partial melt, we applied our results to study the gravitational stability and seismic signature of Fe-rich silicate melts in the mantle, combined with melt geometry and compaction models. For the lowvelocity zone (LVZ) in mantle asthenosphere, although the degree of seismic velocity reduction can be explained by the presence of a small amount of partial melt distributed in film/band geometry along grain boundaries, silicate melts formed in this depth range are too light to be gravitationally stable, but may be completely entrained in the convecting mantle due to the small melt fraction and low meltmatrix separation velocity. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-rich melts (with FeO>similar to 10 wt%) distributed in textural equilibrium with the ambient mantle is a plausible explanation. (C) 2021 Elsevier B.V. All rights reserved.

    Solid-phase transfer into the forearc mantle wedge: Rutile and zircon xenocrysts fingerprint subducting sources

    Aitchison, Jonathan C.Cluzel, DominiqueIreland, Trevor R.Zhou, Renjie...
    10页
    查看更多>>摘要:Voluminous flux of a hydrous component from subducting oceanic lithosphere into the forearc mantle and its contribution to arc magmas have long been recognised but the transfer of solid-phase materials has received less consideration. Although xenocrystic zircons are known from some arc magmas and ophiolitic chromitites their origins remain enigmatic. How and when such materials are transferred into the overlying wedge of mantle lithosphere and the length of their residency therein are unclear. Using zircons and rutiles recovered from ophiolitic rocks in New Caledonia, we demonstrate unambiguous association with subducting sediments and hence evidence for inter-plate transfer. The `fingerprints' of age spectra for both minerals can be matched to lithospheric slab components subducted during Eocene time. Their occurrence in forearc harzburgite suggests that they were incorporated into their mantle wedge host at subduction depths beyond the onset of serpentinisation at ca. 300 degrees C but before the ca. 600 degrees C destabilisation temperature of the rutile U-Pb system. Relocation from slab to forearc mantle is therefore likely to have occurred at ca. 50-60 km depth within the subduction channel. Fluid-assisted conveyance of high field strength elements within solid-state accessory mineral phases into suprasubduction mantle wedges may represent a little-recognised but significant phenomenon that is both predictable and pervasive at convergent margins globally. (C) 2021 The Author(s). Published by Elsevier B.V.

    The role of preexisting upper plate strike-slip faults during long-lived (ca. 30 Myr) oblique flat slab subduction, southern Alaska

    Lease, R. O.Roeske, S. M.Benowitz, J. A.O'Sullivan, P. B....
    12页
    查看更多>>摘要:Upper plates of subduction zones commonly respond to flat slab subduction by structural reactivation, magmatic arc disruption, and foreland basin inversion. However, the role of active strike-slip faults in focusing convergent deformation and magmatism in response to oblique flat slab subduction remains less clear. Here, we present new detrital apatite fission-track (dAFT) ages from 12 modern catchments in the eastern Alaska Range, Alaska, USA, to reveal how the dextral Denali fault system has facilitated bedrock exhumation and topographic growth during ca. 30 Ma-to-present oblique flat slab subduction of the Yakutat oceanic plateau. Additionally, a 940 ka (Ar-40/Ar-39 whole rock) basalt flow is spatially associated with Cenozoic structures, locally reset AFT ages and provides the first evidence for Quaternary volcanism along the southern flank of the eastern Alaska Range. We integrate our new data with other thermochronologic, geochronologic, and regional geologic datasets to show that (1) most high topography regions in southern Alaska have undergone rapid bedrock cooling and exhumation since ca. 30 Ma; (2) elevated terrain and young cooling are spatially associated with long-lived active strike-slip fault systems; (3) topographic growth associated with strike-slip fault deformation led to local inversion of basin systems and drainage reorganization; (4) the onset of oblique oceanic plateau subduction is coeval with a southward shift in arc magmatism from one region of active strike-slip faulting to another above the northeastern edge of the flat slab; and (5) Quaternary volcanism marks the revival of magmatism in the eastern Alaska Range above the geophysically imaged northeastern edge of the flat slab. Our analysis of the post-30 Ma geologic evolution of southern Alaska demonstrates that strike-slip fault systems that were active at the time of slab flattening evolved into transpression zones that focused bedrock cooling, rock exhumation, and topographic growth. (C) 2021 The Authors. Published by Elsevier B.V.

    Linking surface deformation to thermal and mechanical magma chamber processes

    Townsend, Meredith
    14页
    查看更多>>摘要:Surface deformation at active volcanoes reflects a wide variety of magmatic and crustal processes, and the typical kinematic models employed to interpret these geodetic data are highly nonunique. This work presents a new model for surface deformation linked to magma chamber evolution in response to magma injection, eruptions, cooling and crystallization, volatile exsolution, and viscoelastic deformation of the crust. The model is applied to investigate surface displacements during 1) magma chamber crystallization and second boiling without magma injection; 2) a single episode of magma injection to the chamber; and 3) steady magma injection over multiple eruption cycles. Results indicate that the magnitude of surface velocities associated with crystallization and second boiling should not exceed similar to 1 mm/yr unless the magma cooling rate is substantially elevated above normal conductive cooling rates. The surface response to an episode of magma injection depends on magma compressibility and how the timescale of injection compares to the timescale for viscous relaxation. When the injection time is short compared to the relaxation time, and when the magma is equally or less compressible compared to the wall rocks, post-injection deflation can occur without any volume contraction of the chamber. For more compressible magmas, or longer-duration injections, inflation continues over a protracted time following injection. Over multiple eruption cycles, uplift accumulates when more magma is added to the chamber than lost by eruptions; because both the eruption frequency and uplift rate are sensitive to the time-averaged magma supply and chamber volume, geologic data from the eruptive record can be integrated with geodetic data to place tighter constraints on the size and state of subvolcanic magma systems. (C) 2021 Elsevier B.V. All rights reserved.

    A speciation model linking the fate of carbon and hydrogen during core - magma ocean equilibration

    Gaillard, FabriceMalavergne, ValerieBouhifd, Mohamed AliRogerie, Gregory...
    10页
    查看更多>>摘要:The core - mantle differentiation in the magma ocean constitutes a major planetary event that involved two elements that are essential to life: carbon (C) and hydrogen (H). These two elements are conventionally classified as volatiles (i.e. atmophile), but they can also evolve into being siderophile and lithophile at the extreme conditions found in planetary magma oceans. We report here a model for H and C species dissolved in silicate melts in equilibrium with iron-rich alloys under variable pressure, temperature and redox conditions. This speciation model is able to reconcile and reproduce a large body of experimental data on metal-silicate partitioning for H and C at carbon-saturation and in C-undersaturated systems. At low pressure, we conclude that the prevailing species in a silicate magma ocean are CO2, CO, H2O, H-2, whereas CH4 appears to dominate at high pressure. These speciation changes explain recent experimental observations that (i) C evolves from being strongly siderophile at low pressure to moderately siderophile at high pressure, and (ii) H is not siderophile at low pressure but becomes increasingly so as pressure rises. Moreover, it shows that H becomes increasingly siderophile as the total H content of the silicate melt and C-activity are lowered. Despite it offers a promising reconciliation of a large set of experimental and molecular dynamics observations, this model still suffers from large uncertainties when extrapolated to high pressure. In particular, endmember and mixing properties in both the silicate melt and the molten metal must be independently deciphered. The enhanced CH4 stability in the silicate melt at high pressure couples the fate of C and H in deep magma oceans. In such cases, the solubility of C in the basal ocean depends on the H-content and is higher than C-solubility at low pressure. This implies an increase in C activity as the ascending convective cells of the magma ocean, which may cause C saturation as graphite or diamond unless Fe-metal droplets, having a great C-solubility, are present in the shallow magma ocean. Under certain conditions, enhanced siderophile behaviour for H can lead via a runaway process to the desiccation of the magma ocean. (C) 2021 Elsevier B.V. All rights reserved.

    Earth's inner core rotation, 1971 to 1974, illuminated by the inner core scattered waves

    Wang, WeiVidale, John E.
    8页
    查看更多>>摘要:The solid inner core is at the center of the Earth, gravitationally held within the liquid outer core. It is one of the most dynamic parts of Earth's interior. Since the initial claim of inner core differential rotation relative to the mantle, its existence and rate have been challenged for over two decades. Here, we re-examine the seismic records of two megaton nuclear tests in Novaya Zemlya, Russia, three years apart, from the Large Aperture Seismic Array in Montana, USA. Using an improved static time correction from an antipodal earthquake, we refine the resolution of the beamforming of PKiKP and its inner-core scattered coda. Then, we measure the slight time shifts (tenths of seconds) between the inner-corescattered waves from the two events with moving-time-window cross-correlation. Applying a novel back-projection method, we locate the inner-core regions that scatter the energy within the PKiKP coda based on its slowness and the lapse time. We then measure the inner core rotation, first assuming alignment with Earth's rotation axis, then finding the best-fitting differential rotation axis with a grid search. The rotation rates here are robust and consistent across the many scattered arrivals throughout the inner core scattering wavetrain. Our results indicate 0.10 degrees/year inner core super-rotation rate from 1971 to 1974 aligned with Earth's rotation axis, or 0.125 degrees/year with the rotation axis tilting about 8 degrees from the Earth's rotation axis, which yields a marginally better fit to the observed time shifts. (C) 2021 Elsevier B.V. All rights reserved.

    Modeling of zircon nucleation and growth rates using crystal size distributions in a cooling magmatic intrusion

    Sorokin, M. A.Melnik, O. E.Bindeman, I. N.
    11页
    查看更多>>摘要:With a rapid increase in the use of zircons for geochronology and in situisotopic and chemical analyses it is important to quantify growth rates and nucleation rates of zircons in natural magmatic systems. Here we present a mathematical model of the nucleation and growth of a population of zircon crystals in a cooling magmatic intrusion and apply nucleation laws to derive zircon crystal size distributions (CSD) and particle number density (n per cm(3) of melt) from first principles. The model is based on a numerical solution of the one-dimensional diffusion equation that accounts for the nucleation and growth of individual crystals and the dependence of the equilibrium concentration of Zr on temperature. As experimental studies of zircon nucleation at natural concentration levels of Zr are nearly impossible, we rely on measured CSDs of age-invariant (within error of the method <+/- 5ky) zircon populations in rapidly quenched eruptive products of known crystallization duration to calibrate our model. The CSDs of zircon crystals are calculated at different rates of cooling of the magma chamber and are compared to measurements. The combination of the zircon growth software with measurements of zircon CSDs and particle number density permits direct estimation of the zircon nucleation rates and their evolution with time in cooling magmatic systems. We observe a delayed zircon nucleation of 500-2000 years after the melt becomes zircon saturated. Zr supersaturations that drive grown in simulations range from 1 to similar to 30ppmand vary non-monotonically as crystallization proceeds. Initially, the nucleation rates increase exponentially from 10-3 to similar to 1crystals/cm(3)/yrbut then decrease by an order of magnitude as the distance between zircons decreases. Growth rates stay in 10-15 cm/srange for the most duration of zircon growth. The ratio of growth to nucleation rates in our model also varies insignificantly and is usedto estimate crystallization time based on the CSD slope, which decreases with the increased duration of cooling. Flatter slopes correspond to longer crystallization durations as is also noted for major phases. Due to non-monotonic variations in growth and nucleation rates with time, the observed deficiency of small zircons resulting in concave up CSD shapes can be reproduced from first principles assuming monotonic cooling of the magmatic system. Variations in the mode of cooling from linear to monotonic during conductive heat loss of a cooling intrusion, and minor coeval precipitation (< 10vol%) of major phases has little impact on zircon CSD and its evolution. This study has implication to individual zircons IDTIMS dates and demonstrates that most zircons crystallize much later than magma becomes saturated with Zr, despite sometimes having nominally long tail crystallization rate in the inherited cores. For natural CSDs from eight different eruptions randomly selected zircons predate eruption age only by similar to 300-1000 years. (C) 2021 Elsevier B.V. All rights reserved.