Department of Earth and Planetary Science

We introduce five novel types of Monte Carlo (MC) moves that brings the number of moves of ensemble MC calculations from three to eight. So far such calculations have relied on affine invariant stretch moves that were originally introduced by Christen (2007) [8], walk moves by Goodman and Weare (2010) [16] and quadratic moves by Militzer (2023) [31,32]. Ensemble MC methods have been very popular because they harness information about the fitness landscape from a population of walkers rather than relying on expert knowledge. Here we modified the affine method and employed a simplex of points to set the stretch direction. We adopt the simplex concept to quadratic moves. We also generalize quadratic moves to arbitrary order. Finally, we introduce directed moves that employ the values of the probability density while all other types of moves rely solely on the location of the walkers. We apply all algorithms to the Rosenbrock density in 2 and 20 dimensions and to the ring potential in 12 and 24 dimensions. We evaluate their efficiency by comparing error bars, autocorrelation time, travel time, and the level of cohesion that measures whether any walkers were left behind. Our code is open source.

BACKGROUND: Underground research laboratories (URLs) provide a window on the deep biosphere and enable investigation of potential microbial impacts on nuclear waste, CO2 and H2 stored in the subsurface. We carried out the first multi-year study of groundwater microbiomes sampled from defined intervals between 140 and 400 m below the surface of the Horonobe and Mizunami URLs, Japan. RESULTS: We reconstructed draft genomes for > 90% of all organisms detected over a four year period. The Horonobe and Mizunami microbiomes are dissimilar, likely because the Mizunami URL is hosted in granitic rock and the Horonobe URL in sedimentary rock. Despite this, hydrogen metabolism, rubisco-based CO2 fixation, reduction of nitrogen compounds and sulfate reduction are well represented functions in microbiomes from both URLs, although methane metabolism is more prevalent at the organic- and CO2-rich Horonobe URL. High fluid flow zones and proximity to subsurface tunnels select for candidate phyla radiation bacteria in the Mizunami URL. We detected near-identical genotypes for approximately one third of all genomically defined organisms at multiple depths within the Horonobe URL. This cannot be explained by inactivity, as in situ growth was detected for some bacteria, albeit at slow rates. Given the current low hydraulic conductivity and groundwater compositional heterogeneity, ongoing inter-site strain dispersal seems unlikely. Alternatively, the Horonobe URL microbiome homogeneity may be explained by higher groundwater mobility during the last glacial period. Genotypically-defined species closely related to those detected in the URLs were identified in three other subsurface environments in the USA. Thus, dispersal rates between widely separated underground sites may be fast enough relative to mutation rates to have precluded substantial divergence in species composition. Species overlaps between subsurface locations on different continents constrain expectations regarding the scale of global subsurface biodiversity. CONCLUSIONS: Our analyses reveal microbiome stability in the sedimentary rocks and surprising microbial community compositional and genotypic overlap over sites separated by hundreds of meters of rock, potentially explained by dispersal via slow groundwater flow or during a prior hydrological regime. Overall, microbiome and geochemical stability over the study period has important implications for underground storage applications.

A method is presented for high-precision chemical detection that integrates quantum sensing with droplet microfluidics. Using nanodiamonds (ND) with fluorescent nitrogen-vacancy (NV) centers as quantum sensors, rapidly flowing microdroplets containing analyte molecules are analyzed. A noise-suppressed mode of optically detected magnetic resonance is enabled by pairing controllable flow with microwave control of NV electronic spins, to detect analyte-induced signals of a few hundredths of a percent of the ND fluorescence. Using this method, paramagnetic ions in droplets are detected with low limit-of-detection using small analyte volumes, with exceptional measurement stability over >10³ s. In addition, these droplets are used as microconfinement chambers by co-encapsulating ND quantum sensors with various analytes such as single cells, suggesting wide-ranging applications including single-cell metabolomics and real-time intracellular measurements from bioreactors. Important advances are enabled by this work, including portable chemical testing devices, amplification-free chemical assays, and chemical imaging tools for probing reactions within microenvironments.

Cyanobacteria are highly abundant in the marine photic zone and primary drivers of the conversion of inorganic carbon into biomass. To date, all studied cyanobacterial lineages encode carbon fixation machinery relying upon form I Rubiscos within a CO2-concentrating carboxysome. Here, we report that the uncultivated anoxic marine zone (AMZ) IB lineage of Prochlorococcus from pelagic oxygen-deficient zones (ODZs) harbors both form I and form II Rubiscos, the latter of which are typically noncarboxysomal and possess biochemical properties tuned toward low-oxygen environments. We demonstrate that these cyanobacterial form II enzymes are functional in vitro and were likely acquired from proteobacteria. Metagenomic analysis reveals that AMZ IB are essentially restricted to ODZs in the Eastern Pacific, suggesting that form II acquisition may confer an advantage under low-O2 conditions. AMZ IB populations express both forms of Rubisco in situ, with the highest form II expression at depths where oxygen and light are low, possibly as a mechanism to increase the efficiency of photoautotrophy under energy limitation. Our findings expand the diversity of carbon fixation configurations in the microbial world and may have implications for carbon sequestration in natural and engineered systems.

High-resolution full waveform seismic tomography of the Earths mantle beneath the south and central Atlantic Ocean brings into focus a series of asthenospheric low shear velocity channels, or fingers on both sides of the southern and central mid-Atlantic ridge (MAR), elongated in the direction of absolute plate motion with a spacing of [Formula: see text]1,800 to 2,000 km, and associated with bands of shallower residual seafloor depth anomalies that suggest channeled flow over thousands of kilometers. Each of the three most clearly resolved fingers on the African side of the MAR corresponds to a separate group of whole mantle plumes rooted in distinct patches at the core-mantle boundary, feeding hotspots, and volcanic lines with distinct isotopic signatures. Plumes of a given group appear to merge at the top of the lower mantle before separating again, suggesting interaction of deep mantle flow with a more vigorous mesoscale circulation in the upper mantle. The corresponding hotspots are generally offset from the location of the deep mantle plume roots. The distinct isotopic signatures of these hotspot groups are also detected in the mid-ocean ridge basalts at the location where the fingers meet the ridge. Meanwhile, at least some of the variability within each plume group could originate in the upper mantle and extended transition zone where plumes in a given group appear to merge and pond. This study also adds to mounting evidence that the African large low shear velocity province is not a uniform, unbroken pile of dense material rising high above the core-mantle boundary, but rather a collection of mantle plumes rooted in patches of distinct composition.

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

Mafic intrusions, lava flows, and felsic plutons in southwestern Laurentia have been hypothesized to be associated with the emplacement of a late Mesoproterozoic (Stenian Period) large igneous province. Improved geochronologic data resolve distinct episodes of mafic magmatism in the region. The ca. 1,098 Ma main pulse of southwestern Laurentia large igneous province (SWLLIP) magmatism is recorded by mafic intrusions across southeastern California to central Arizona. A younger episode of volcanism resulted in eruptions that formed the ca. 1,082 Ma Cardenas Basalt, which is the uppermost unit of the Unkar Group in the Grand Canyon. With the updated geochronological constraints, we develop new paleomagnetic data from mafic sills in the SWLLIP. Overlapping poles between the Death Valley sills and rocks of similar age in the Midcontinent Rift are inconsistent with large-scale Cenozoic vertical axis rotations in Death Valley. We also develop a new paleomagnetic pole from the ca. 1,082 Ma Cardenas Basalt (pole longitude = 183.9°E, pole latitude = 15.9°N, (Formula presented.) = 7.4°, N = 18). The new paleomagnetic data are consistent with the pole path developed from time-equivalent rocks of the Midcontinent Rift, supporting interpretations that changing pole positions are the result of rapid equatorward motion. These data add to the record of Laurentia's rapid motion from ca. 1,110 to 1,080 Ma that culminated in collisional Grenvillian orogenesis and the assembly of Rodinia.

UNLABELLED: As the closest living relatives of animals, choanoflagellates offer insights into the ancestry of animal cell physiology. Here, we report the isolation and characterization of a colonial choanoflagellate from Mono Lake, California. The choanoflagellate forms large spherical colonies that are an order of magnitude larger than those formed by the closely related choanoflagellate Salpingoeca rosetta. In cultures maintained in the laboratory, the lumen of the spherical colony is filled with a branched network of extracellular matrix and colonized by bacteria, including diverse Gammaproteobacteria and Alphaproteobacteria. We propose to erect Barroeca monosierra gen. nov., sp. nov. Hake, Burkhardt, Richter, and King to accommodate this extremophile choanoflagellate. The physical association between bacteria and B. monosierra in culture presents a new experimental model for investigating interactions among bacteria and eukaryotes. Future work will investigate the nature of these interactions in wild populations and the mechanisms underpinning the colonization of B. monosierra spheres by bacteria. IMPORTANCE: The diversity of organisms that live in the extreme environment of Mono Lake (California, USA) is limited. We sought to investigate whether the closest living relatives of animals, the choanoflagellates, exist in Mono Lake, a hypersaline, alkaline, arsenic-rich environment. We repeatedly isolated members of a new species of choanoflagellate, which we have named Barroeca monosierra. Characterization of B. monosierra revealed that it forms large spherical colonies containing diverse co-isolated bacteria, providing an opportunity to investigate mechanisms underlying physical associations between eukaryotes and bacteria.

Large volumes of liquid water transiently existed on the surface of Mars more than 3 billion years ago. Much of this water is hypothesized to have been sequestered in the subsurface or lost to space. We use rock physics models and Bayesian inversion to identify combinations of lithology, liquid water saturation, porosity, and pore shape consistent with the constrained mid-crust (∼11.5 to 20 km depths) seismic velocities and gravity near the InSight lander. A mid-crust composed of fractured igneous rocks saturated with liquid water best explains the existing data. Our results have implications for understanding Mars' water cycle, determining the fates of past surface water, searching for past or extant life, and assessing in situ resource utilization for future missions.

Present and future climatic trends are expected to markedly alter water fluxes and stores in the hydrologic cycle. In addition, water demand continues to grow due to increased human use and a growing population. Sustainably managing water resources requires a thorough understanding of water storage and flow in natural, agricultural, and urban ecosystems. Measurements of stable isotopes of water (hydrogen and oxygen) in the water cycle (atmosphere, soils, plants, surface water, and groundwater) can provide information on the transport pathways, sourcing, dynamics, ages, and storage pools of water that is difficult to obtain with other techniques. However, the potential of these techniques for practical questions has not been fully exploited yet. Here, we outline the benefits and limitations of potential applications of stable isotope methods useful to water managers, farmers, and other stakeholders. We also describe several case studies demonstrating how stable isotopes of water can support water management decision-making. Finally, we propose a workflow that guides users through a sequence of decisions required to apply stable isotope methods to examples of water management issues. We call for ongoing dialogue and a stronger connection between water management stakeholders and water stable isotope practitioners to identify the most pressing issues and develop best-practice guidelines to apply these techniques.