Geochemistry Geophysics Geosystems最新文献

Paleosecular variation analysis is a primary tool for characterizing ancient geomagnetic behavior and its evolution through time. This study presents a new high-quality directional data set, paleosecular variation of the Paleogene (PSVP), with and without correction for serial correlation, compiled from 1,667 sites from 45 different localities from the Paleogene and late Cretaceous (84–23 Ma). The data set is used to study the variability, structure, and latitude dependence of the geomagnetic field during that period by varying selection criteria and PSV models. Modeled values for the equatorial virtual geomagnetic pole (VGP) dispersion have over-lapping uncertainty intervals within their uncertainty bounds between 8.3° and 18.6° for the past 250 Ma. We investigate the suitability of two descriptive models of PSV, Model G-style quadratic fits and covariant Giant Gaussian Process models, and find that both styles of model fail to satisfactorily reproduce the latitude dependent morphology of PSV, but suggest that estimates of the equatorial VGP dispersion may still robustly characterize aspects of Earth's long-term field morphology. During this time where the PSV behavior has not changed substantially, the reversal frequency has varied widely. The lack of a clear relationship between PSV behavior and reversal frequency is not trivially explained in the context of published findings regarding numerical geodynamo simulations.

The proliferation of microbial carbon fixation is a key control on the evolution of the biosphere and global carbon cycle. Most records of these metabolisms in ancient rocks come from organic matter or fossils, which are not always preserved. Here, we present a potential proxy for microbial carbon fixation (autotrophy) based on the isotopic composition of carbonate minerals. Autotrophs influence carbonate chemistry in the cellular microenvironment by decreasing CO₂ concentration and increasing the carbonate saturation state. This can induce rapid precipitation of carbonate minerals that are out of isotopic equilibrium with their environment. Recent work has identified disequilibrated dual clumped isotope compositions (∆₄₇ and ∆₄₈) in the skeletal fossils of marine calcifying organisms. Here we test whether the same is true of non-skeletal carbonate fabrics associated with microbial autotrophs in modern and Eocene lakes. We found that microbial carbonate formed via autotrophic metabolism recorded lower ∆₄₇ and higher ∆₄₈ values (−∆₄₇/+∆₄₈) than predicted for thermodynamic equilibrium mineral formation. Our findings are supported by models of dual clumped isotope kinetics in the DIC system, and disequilibrium in the oxygen isotope system. We hypothesize that the inverse trajectory away from the equilibrium line (+∆₄₇/−∆₄₈) should be recorded by carbonates formed in association with alkalinizing heterotrophs, such as sulfate reducers. If so, carbonate dual clumped isotopes could be a powerful tool to identify the proliferation and rate of heterotrophic and autotrophic metabolisms in the carbonate rock record on Earth and (perhaps) other planets.

Anomalously overturned thrust faults, lineaments and segmentation causing cross-cutting basement structures characterize the tectonic setting of Sikkim Himalaya. However, its seismotectonics is poorly constrained along with the speculated northward extension of the Dhubri-Chungthang Fault Zone (DCFZ) causing segmentation. Here, we utilize the precise location of newly acquired local earthquake data and fault plane solutions using full-waveform moment tensor inversion to better constrain seismically active zones. Transtensional shearing along the Main Himalayan Thrust in central Sikkim is possibly incited by fluid-rich upper-crust. Cessation of the mapped 20 km wide mid-crustal seismogenic zone of DCFZ at Chungthang and its northward discontinuation into the Higher Himalayas is confirmed by the striking variation in focal mechanisms. Earthquakes along imbricated segments in the lower-crust originate possibly in response to crustal shortening. Extensional shearing along the Moho triggers seismicity to the northwest of Sikkim. Such complex tectonic dynamics instigating persistent seismicity indicates high potential for future great earthquakes in Sikkim Himalaya.

The 2021 La Palma eruption (Tajogaite) was unprecedented in magnitude, duration, and degree of monitoring compared to historical volcanism on La Palma. Here, we provide data on melt inclusions in samples from the beginning and end of the eruption to compare the utility of both melt and fluid inclusions as recorders of magma storage. We also investigated compositional heterogeneities within the magmatic plumbing system. We found two populations of olivine crystals: a low Mg# (78–82) population present at the beginning and end of eruption, recording the maximum volatile contents (2.5 wt % H₂O, 1,800 ppm F, 700 ppm Cl, 3,800 ppm S) and a higher Mg# (83–86) population sampled toward the end of the eruption, with lower volatile contents. Despite their host composition, melt inclusions share the same maximum range of CO₂ concentrations (1.2–1.4 wt %), indicating olivine growth and inclusion capture at similar depths. Overall, both melt and fluid inclusions record similar pressures (450–850 MPa, ∼15–30 km), and when hosted in the same olivine crystal pressures are indistinguishable within error. At these mantle pressures, CO₂ is expected to be an exsolved phase explaining the similar range of CO₂ between the two samples, but other volatile species (F, Cl, S) behave incompatibly, and thus, the increase between the two olivine populations can be explained by fractional crystallization prior to eruption. Finally, based on our new data, we provide estimates on the total volatile emission of the eruption.

The shapes of mantle plumes are sensitive to mantle viscosity, density structure, and flow patterns. Increasingly, global tomographic models reveal broad plume conduits in the lower mantle and highly tilting conduits in the mid and upper mantle. Previous studies mostly relied on 2D slices to analyze plume shapes, but fully investigating the complexity of 3D plume structures requires more effective visualization methods. Here, we use immersive headset-based virtual reality (VR) to visualize the full-waveform global tomographic models SEMUCB-WM1 and GLAD-M25. We develop criteria for the identification of plume conduits based on the relationship between the plume excess temperature and the V_S anomaly (δV_S). We trace 20 major plume conduits, measure the offsets of the conduits in azimuth and distance with respect to the hotspots, calculate the tilt angle, and evaluate the δV_S along all traced conduits. We compare our traced conduits with the conduits predicted by global mantle convection models and vertical conduits. The wavespeed variations along conduits traced from each tomographic model are slower than modeled or vertical conduits, regardless of which tomographic model they are evaluated in. The shapes of traced conduits tend to differ greatly from modeled conduits. Plume ponding and the emergence of secondary plumes, which could result from a combination of compositional variations, phase transitions, small-scale convection, and variations in viscosity, can contribute to the complex observed plume shapes. The variation of δV_S along the traced conduits and complex plume shapes suggest a thermochemical origin of many plumes.

Understanding the behavior of mercury (Hg) in organic-rich sediments as they undergo thermal maturation is important, for example, because enrichment of Hg in sedimentary deposits has become a widely used proxy for volcanism from Large Igneous Provinces (LIPs). In this study, we evaluate the effects of such processes on sedimentary Hg concentrations by investigating a common stratigraphic interval in three drill cores with different levels of thermal maturity (immature, mature and post-mature) in Toarcian sediments (Posidonienschiefer Formation) from the Lower Saxony Basin, Germany. We present Hg concentrations, bulk organic geochemistry, and total sulfur data. Mercury concentrations in the mature and post-mature sediments are increased >2-fold relative to the immature material, which is greater than any potential differences in original Hg concentrations in the studied successions prior to burial. Organic-carbon and host-rock mass loss during thermal maturation may have concentrated Hg in the mature sediments to some extent, provided Hg is considered effectively immobile. The increased Hg, TOC-normalized Hg, and TS-normalized Hg are most likely linked to the “closed system” behavior of Hg in sedimentary basins and the relatively low temperatures (70–260°C) during maturation that resulted in limited Hg mobility. More speculatively, a certain degree of redistribution of Hg within the mature sediments is suggested by its enrichment in distinct stratigraphic levels. Regardless of the exact mechanisms at play, the elevated Hg concentrations in mature sediments amplify both Hg/TOC and Hg/TS, implying that thermal effects must be considered when using normalized Hg as a proxy for far-field volcanic activity.

(Mg,Fe)O ferropericlase-magnesiowüstite has been proposed to host the majority of Earth's sodium, but the mechanism and capacity for incorporating the alkali cation remain unclear. In this work, experiments in the laser-heated diamond anvil cell and first-principles calculations determine the solubility of sodium and favorability of sodium incorporation in iron-rich magnesiowüstite relative to (Mg,Fe)SiO₃ bridgmanite. Reaction of Mg/(Mg + Fe) (Mg#) 55 and 28 olivine with NaCl at 33–128 GPa and 1600–3000 K produces iron-rich magnesiowüstite containing several percent sodium, while iron-rich bridgmanite contains little to no detectable sodium. In sodium-saturated magnesiowüstite, sodium number [Na/(Na + Mg + Fe)] is 2–5 atomic percent at pressures below 60 GPa and drastically increases to 10–20 atomic percent at deep lower mantle pressures. For these two compositions, there is no significant dependence of the results on Mg#. Our calculations not only show consistent results with experiments but further indicate that such an increase in solubility and partitioning of Na into magnesiowüstite is driven by the spin transition in iron. These results provide fundamental constraints on the crystal chemistry of sodium at lower-mantle conditions. If the sodium capacity of (Mg,Fe)O is not strongly dependent on Mg#, (Mg,Fe)O in the lower mantle may have the capacity to store the entire sodium budget of the Earth.

We analyzed 49,592 teleseismic receiver functions (RFs) recorded by 278 CEArray stations to image the mantle transition zone (MTZ) beneath the South China Block to understand the origins of deep velocity anomalies and their potential links to subduction and intraplate volcanism. We employed a fast-marching method and a high-resolution 3-D velocity model (FWEA18) derived from full waveform inversion in computing P-to-S conversion times to better image the 410- and 660-km discontinuities. Our results indicate that the common-conversion-point stacking of RFs using 3-D conversion times yielded better migration images of the two discontinuities. The images revealed a slightly depressed 410-km with a few small uplifted patches, and showed that the 660-km beneath the western Yangtze Craton is depressed by 10–25 km, which is likely caused by the stagnant Paleo-Pacific slab. The 660-km beneath the southern Cathaysia Block has a 5–15 km high plateau with a topographic low at its central part. The lateral dimension of the topographic low is ∼150 km and is located beneath the central Pearl River Mount Basin near Hong Kong. We speculate that the topographic low occurs within the Hainan plume with a temperature excess of ∼300–400 K and is caused by the garnet phase transition. The displaced deep plume enters the MTZ and spreads nearly horizontally at the base. The plume evolves into two channels with a minor one toward the northeast and a major one toward the southwest, which keep moving upward to the 410-km. The southwest channel is likely the source that feeds the Hainan volcanoes.

The seismicity of the Philippines and Taiwan provides insight into the tectonics and seismic hazards of a region characterized by subduction and collision. We summarize the seismotectonics of the Philippines and Taiwan by documenting the distribution of hypocenters for earthquakes of magnitude M ≥ 4.6 and focal mechanisms for earthquakes of magnitude M ≥ 5.0 over ∼21 years. We quantify seismicity rates (earthquake frequency) and compare seismicity distributions with proposed tectonic and faulting models. 6,187 earthquakes of magnitude M ≥ 4.6 occurred between 1 January 2000 and 31 March 2021, 79% at depths <70 km and 70% having magnitudes M < 5.0. Approximately 88 earthquakes of magnitude M ≥ 5.0 occur per year, with 12 events of magnitude M ≥ 7.0 occurring since January 2000. Seismic activity decreases exponentially between 50 and 210 km depth at a rate ∼10% faster than the global average. Intermediate and deep earthquakes at depths >70 km trace the Wadati-Benioff zones of subducting slabs, most of which are only seismically active to depths of ∼250 km. The distribution of earthquakes at depths >70 km is likely influenced by the subduction of young lithosphere, slab tearing, and phase boundary interactions between depths of 410 and 660 km. Shallow earthquakes at depths ≤70 km are generated by megathrust, crustal, and intraslab faulting. Crustal thrust and strike-slip faulting are the most abundant and prevalent sources of damaging earthquakes. The Philippines and Taiwan are subject to high seismic risk, similar to nearby Indonesia.

Fluctuations in the geomagnetic field occur over a broad range of timescales. Short-period fluctuations are called secular variation, whereas excursions and reversals are viewed as anomalous transient events. An open question is whether distinct mechanisms are required to account for these different forms of variability. Clues are sought in trends b of the axial dipole moment from six time-dependent geomagnetic field models. Variability in b has a well-defined dependence on the time interval (or window) for the trend. The variance of b reveals a simple relationship to trends during excursions and reversals. This connection hints at a link between reversals, excursions and secular variation. Stochastic models exhibit a similar behavior in response to random fluctuations in dipole generation. We find that excursions, reversals and secular variation can be distinguished on the basis of trend durations rather than differences in the underlying physical process. While this analysis does not rule out distinct physical mechanisms, the paleomagnetic observations suggest that such distinctions are not required.