We present our third and final generation joint P and S global adjoint tomography (GLAD) model, GLAD-M35, and quantify its uncertainty based on a low-rank approximation of the inverse Hessian. Starting from our second-generation model, GLAD-M25, we added 680 new earthquakes to the database for a total of 2160 events. New P-wave categories are included to compensate for the imbalance between P- and S-wave measurements, and we enhanced the window selection algorithm to include more major-arc phases, providing better constraints on the structure of the deep mantle and more than doubling the number of measurement windows to 40 million. Two stages of a Broyden–Fletcher–Goldfarb–Shanno (BFGS) quasi-Newton inversion were performed, each comprising five iterations. With this BFGS update history, we determine the model’s standard deviation and resolution length through randomized singular value decomposition.

We present a suite of major element stable isotope (δ13C, δ18O, δ44/40Ca, δ26Mg), and selected trace element (Sr/Ca and Mg/Ca) data from Pleistocene sediments from the Great Barrier Reef (IODP Expedition 325), as well as Holocene surface sediments from the Bahamas (Triple Goose Creek, Andros Island) to identify geochemical fingerprints associated with early marine and meteoric diagenesis. Sediments from both sites exhibit co-variation in δ13C and δ18O values, depletion in trace elements, and distinct geochemical trends in δ26Mg and δ44/40Ca values that reflect differences between diagenetic alteration in marine and meteoric fluids. While marine diagenesis results in lower Sr/Ca ratios, higher δ44/40Ca values, and little effect on bulk sediment δ26Mg values, meteoric diagenesis leads to lower Sr/Ca ratios, lower δ44/40Ca values, and lower δ26Mg values. Using a numerical model of diagenesis, we show how diagenetic alteration by meteoric fluids must occur after an initial period of diagenetic alteration by marine fluids, a two-stepped diagenetic history that complicates the interpretation of geochemical data in meteorically altered marine carbonate sediments. Finally, we discuss how paired metal isotopes may serve as a robust indicator of meteoric alteration in ancient shallow-water marine carbonate sediments. © 2024 Elsevier Ltd

Radiogenic heat production is fundamental to the energy budget of planets. Roughly half of the heat that Earth loses through its surface today comes from the three long-lived, heat-producing elements (potassium, thorium, and uranium). These three elements have long been believed to be highly lithophile and thus concentrate in the mantle of rocky planets. However, our study shows that they all become siderophile under the pressure and temperature conditions relevant to the core formation of large rocky planets dubbed super-Earths. Mantle convection in super-Earths is then primarily driven by heating from the core rather than by a mix of internal heating and cooling from above as in Earth. Partitioning these sources of radiogenic heat into the core remarkably increases the core-mantle boundary (CMB) temperature and the total heat flow across the CMB in super-Earths. Consequently, super-Earths are likely to host long-lived volcanism and strong magnetic dynamos. Entrainment of heat-producing elements in super-Earths’ cores produces intense, long-lasting volcanism and strong magnetic fields.

Crosstalk-free source-encoded elastic full-waveform inversion (FWI) using time-domain solvers demonstrates skill and efficiency at conducting seismic inversions involving multiple sources and receivers with limited computational resources. A drawback of common formulations of the procedure is that, by sweeping through the frequency domain randomly at a rate of one or a few sparsely sampled frequencies per shot, it is difficult to simultaneously incorporate time-selective data windows, as necessary for the targeting of arrivals or wave packets during the various stages of the inversion. Here, we solve this problem by using the Laplace transform of the data. Using complex-valued frequencies allows for damping the records with flexible decay rates and temporal offsets that target specific traveltimes. We present the theory of crosstalk-free source-encoded FWI in the Laplace domain, develop the details of its implementation, and illustrate the procedure with numerical examples relevant to exploration-scale scenarios. © 2024 Society of Exploration Geophysicists. All rights reserved.

We present a computational technique to model hydroacoustic waveforms from teleseismic earthquakes recorded by mid-column MERMAID floats deployed in the Pacific, taking into consideration bathymetric effects that modify seismo-Acoustic conversions at the ocean bottom and acoustic wave propagation in the ocean layer, including reverberations. Our approach couples axisymmetric spectral-element simulations performed for moment-Tensor earthquakes in a 1-D solid Earth to a 2-D Cartesian fluid solid coupled spectral-element simulation that captures the conversion from displacement to acoustic pressure at an ocean-bottom interface with accurate bathymetry. We applied our w orkflo w to 1129 seismograms for 682 earthquakes from 16 MERMAID s (short for Mobile Earthquake Recording in Marine Areas by Independent Divers) owned by Princeton University that were deployed in the Southern Pacific as part of the South Pacific Plume Imaging and Modeling (SPPIM) project. We compare the modelled synthetic waveforms to the observed records in indi viduall y selected frequency bands aimed at reducing local noise levels while maximizing earthquake-generated signal content. The modelled waveforms match the observations very well, with a median correlation coefficient of 0.72, and some as high as 0.95. We compare our correlation-based traveltime measurements to measurements made on the same data set determined by automated arri v al-Time picking and ra y-traced tra veltime predictions, with the aim of opening up the use of MERMAID records for global seismic tomography via full-waveform inversion. © The Author(s) 2024.