Geochemistry Geophysics Geosystems最新文献

The effect of mantle plumes is secondary to that of subducting slabs for modern plate tectonics when considering plate driving forces. However, the impact of plumes on tectonics and planetary surface evolution may nonetheless have been significant. We use numerical mantle convection models in a 3-D spherical chunk geometry with damage rheology to study some of the dynamics of plume-slab interactions. Substantiating our earlier 2-D results, we observe a range of interaction scenarios, and that the plume-driven subduction terminations we had identified earlier persist in more realistic convective flow. We analyze the dynamics of plume affected subduction, including in terms of their geometry, frequency, and the overall effect of plumes on surface dynamics as a function of the fraction of internal to bottom heating. Some versions of such plume-slab interplay may be relevant for geologic events, for example, for the inferred ∼183 Ma Karoo large igneous province formation and associated slab disruption. More recent examples may include the impingement of the Afar plume underneath Africa leading to disruption of the Hellenic slab, and the current complex structure imaged for the subduction of the Nazca plate under South America. Our results imply that plumes may play a significant role not just in kick-starting plate tectonics, but also in major modifications of slab-driven plate motions, including for the present-day mantle.

The late-stage union of Pangea was associated with the convergence of Siberia with Laurussia, but the exact timing remains unclear. The orogenic duration of the Kazakhstan block can provide geochronological constraints as it connects Siberia, Baltica, and Tarim. Zircon petrochronology offers a reliable approach for ascertaining the lifespan of an ancient orogen. In this study, we explore three phases of magmatism recorded in detrital zircons from late Paleozoic and Mesozoic sandstone-siltstones in Western Tianshan, that is, 400–470 Ma, 320–380 Ma, and 280–320 Ma. Based on their age-propagated Hf isotopes, melt SiO₂ contents, and crustal thicknesses, our findings suggest that the southern limb of Kazakhstan underwent the early Paleozoic amalgamation of microcontinents with arcs, the late Paleozoic maturation of an Andean-like continental arc, and the late Carboniferous collision of Kazakhstan with the Junggar oceanic basin and the Tarim craton. Such characteristics manifest the long-term orogenic progression of Kazakhstan. Combining published timelines and paleolatitudes of major orogens and blocks, we propose that the Kazakhstan block welded northeastern Pangea along with the cessation of these orogenic activities around it. Consequently, by docking Kazakhstan with surrounding cratons, the fundamental configuration of Pangea could have been established in the late Carboniferous.

Most of the current helium (He) reserves originate from fortuitous discoveries, mainly made during oil and gas exploration in sedimentary basins. As helium generation depends on U and Th α-decay, old geological provinces gather key ingredients for high He accumulation. However, numerous He-rich springs have also been documented in much younger rocks, such as Variscan granites (320–250 Ma). These latter discoveries question the current exploration guidelines and require revisiting some of the longstanding paradigms. Here, is investigated He migration along a major fault rooted in the Corso-Sardinian batholith (France). Two thermal springs, Caldanelle and Guagno-Les-Bains, show significant outgassing activities of crustal sourced He with concentrations up to 1.45 vol% and flow rates of 110 m³ STP ⁴He/year. Besides He, the gas phase is dominated by N₂ (≈98 vol%) and minor CH₄. Based on a survey employing multidisciplinary methodologies, it is revealed that (a) Variscan rocks represent efficient ⁴He source rocks, (b) the main source of He comes from the underlying Eo-Variscan basement, (c) A deeply rooted fault and dense fractures networks drain the He, (d) the helium loss is limited, (e) faults and fractures may act as partial traps, and finally (f) the presence of an efficient trap could promote a He-rich reservoir with high flux but low reserves. In that sense, young post-orogenic granites represent promising helium plays. The geological context in which Caldanelle and Guagno-Les-Bains are embedded is ubiquitous in European Variscan batholiths. This case study is therefore intended to serve as a guide for helium exploration and to provide insights into helium behavior within a Variscan geological context.

There are various hotspots in the Atlantic Ocean, which are underlain by mantle plumes that likely cross the mantle and originate at the core-mantle boundary. We use teleseismic core-diffracted shear waves to look for an Ultra-Low Velocity Zone (ULVZ) at the potential base of central Atlantic mantle plumes. Our data set shows delayed postcursory phases after the core-diffracted shear waves. The observed patterns are consistent in frequency dependence, delay time, and scatter pattern with those caused by mega-ULVZs previously modeled elsewhere. Synthetic modeling of a cylindrical structure on the core-mantle boundary below St. Helena provides a good fit to the data. The preferred model is 600 km across and 20 km high, centered at approximately 15° South, 15° West, and with a 30% S-wave velocity reduction. Significant uncertainties and trade-offs do remain to these parameters, but a large ULVZ is needed to explain the data. The location is west of St. Helena and south of Ascension. Helium and neon isotopic systematics observed in samples from this region could point to a less-outgassed mantle component mixed in with the dominant signature of recycled material. These observations could be explained by a contribution from the Large Low Shear Velocity Province (LLSVP). Tungsten isotopic measurements would be needed to understand whether a contribution from the mega-ULVZ is also required at St. Helena or Ascension.

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Computing the velocity field is an expensive process for mantle convection codes. This has implications for particle methods used to model the advection of quantities such as temperature or composition. A common choice for the numerical treatment of particle trajectories is classical fourth-order explicit Runge–Kutta (ERK4) integration, which involves a velocity computation at each of its four stages. To reduce the cost per time step, it is possible to evaluate the velocity for a subset of the four time integration stages. We explore two such alternative schemes, in which velocities are only computed for: (a) stage 1 on odd-numbered time steps and stages 2–4 for even-numbered time steps, and (b) stage 1 for all time steps. A theoretical analysis of stability and accuracy is presented for all schemes. It was found that the alternative schemes are first-order accurate with stability regions different from that of ERK4. The efficiency and accuracy of the alternate schemes were compared against ERK4 in four test problems covering isothermal, thermal, and thermochemical flows. Exact solutions were used as reference solutions when available. In agreement with theory, the alternate schemes were observed to be first-order accurate for all test problems. Accordingly, they may be used to efficiently compute solutions to within modest error tolerances. For small error tolerances, however, ERK4 was the most efficient.

While basaltic volcanism is dominant during rifting and continental breakup, felsic magmatism may be a significant component of some rift margins. During International Ocean Discovery Program (IODP) Expedition 396 on the continental margin of Norway, a graphite-garnet-cordierite bearing dacitic unit (the Mimir dacite) was recovered in two holes within early Eocene sediments on Mimir High (Site U1570), a marginal high on the Vøring Transform Margin. Here, we present a comprehensive textural, petrological, and geochemical study of the Mimir dacite in order to assess its origin and discuss the geodynamic implications. The major mineral phases (garnet, cordierite, quartz, plagioclase, alkali feldspar) are hosted in a fresh rhyolitic, vesicular, glassy matrix that is locally mingled with sediments. The major element chemistry of garnet and cordierite, the presence of zircon inclusions with inherited cores, and thermobarometric calculations all support an upper crustal metapelitic origin. While most magma-rich margin models favor crustal anatexis in the lower crust, thermobarometric calculations performed here show that the Mimir dacite was produced at upper-crustal depths (<5 kbar, 18 km depth) and high temperature (750–800°C) with up to 3 wt% water content. In situ U-Pb analyses on zircon inclusions give a magmatic crystallization age of 54.6 ± 1.1 Ma, consistent with emplacement that post-dates the Paleocene-Eocene Thermal Maximum. Our results suggest that the opening of the Northeast Atlantic was associated with a phase of low-pressure, high-temperature crustal anatexis preceding the main phase of magmatism.

Significant and readily accessible orogenic gold deposits have been previously recognized, exploited, and progressively depleted. Innovative approaches are required to discover new and deeply buried deposits. Recently, trace element variations in apatite have been used to distinguish fertile and barren environments as reliable mineral exploration tools. In this study, machine learning models using Random Forest and Deep Neutral Network are utilized to assess the fertility of quartz veins and altered zones in the orogenic gold systems. The two models have been trained using trace element data of apatite, and the performance of both models yield good classification accuracy (∼90% on average) with low false positive rates. Feature importance analysis (Gini decrease and hidden layer weights) suggest that Pb, As, U, Sr, Eu, Mn, and Fe are the important parameters. Arsenic, U, Eu, Mn, and Fe are redox-sensitive elements, with their concentrations responding to changes in fluid redox conditions. Strontium primarily originates from the breakdown of plagioclase, which is more likely to occur under oxidizing fluid conditions. Therefore, we infer that the main controlling factor of the models is the redox conditions. The two distinct models consistently highlight the most significant contribution of Pb to this differentiation, even though Pb is not a redox-sensitive element and can only substitute for Ca²⁺ in apatite as Pb²⁺. We infer that the high contribution of Pb may be attributed to the potential transportation of Au in the form of a Pb-(Bi)-Au melt, and the Pb content in apatite is influenced by the Pb content in the melt, fluid oxygen, and sulfur fugacity. We also propose a novel discriminant plot using Linear Discriminant Analysis (LDA) to assess the mineralization potential in quartz veins and alteration zones based on apatite trace element data. The machine learning and LDA results suggest that apatite trace elements could be used effectively in the future orogenic gold deposit exploration.

The southeastern United States Coastal Plain ecosystem contains baldcypress (Taxodium distichum) swamps and grass-dominated marshes. These ecosystems also occurred on the exposed continental shelf during lower sea levels but are rarely preserved due to the mechanically erosive nature of transgression and regression. Two presently marine sites on the northeastern Gulf of Mexico's continental shelf contain well-preserved woody terrestrial sediments that were the subject of previous studies. This study continues the investigation using geochemical (δ¹³C, δ¹⁵N, δ³⁴S) and palynological characteristics of these formerly terrestrial sediments to determine if swamps and/or marshes existed at the time of deposition. The first site is located ∼20 km southeast of Horn Island, Mississippi (MS) and the core has terrestrial sediments radiocarbon dated to 11,066–10,228 (2σ) calibrated years BP (early Holocene). The second site is the “Alabama Underwater Forest” located ∼13 km south of Gulf Shores, Alabama (AL) and the cores have terrestrial sediments optically stimulated luminescence dated to 63 ka (±10 ka, 2σ) to 72 ka (±16 ka, 2σ) (late Pleistocene). Geochemical results for the MS sediments indicate a swamp-to-freshwater marsh transitional series, whereas the AL sediments indicate a swamp-to-saltwater marsh transitional series, both supported by palynological results. Further exploration of the geochemical results using linear discriminant analysis, trained with published geochemical data, supports the swamp and marsh interpretations. We conclude that the near-pristine preservation of these woody deposits is not solely due to physical mechanisms, such as rapid burial, but is also coupled with anoxia- and euxinia-driven biogeochemical reactions promoting wood and woody debris preservation in swamp and marsh environments.

Various types of Global Navigation Satellite System (GNSS) data are used for a wide range of applications. When modeled correctly, millimeter precision daily GNSS position time-series yield velocities and other derived products that can be used for investigations of lithospheric processes and properties. In this review paper, we describe the specific types of GNSS data and data products that are valuable for studies of the lithosphere, such as coseismic offsets, post-seismic decay in time-series, seasonal signals, secular velocities, and strain rates, and how those data are derived. We also discuss the applications of several types of GNSS data and data products. We provide open access resources for precision GNSS daily position time-series, quality secular velocity solutions, and daily GNSS RINEX files for researchers interested in processing their own data.