Department of Earth and Planetary Science

The northern lowlands of early Mars could have contained a significant quantity of liquid water. However, the ocean hypothesis remains controversial due to the lack of conclusive evidence from the Martian subsurface. We use data from the Zhurong Rover Penetrating Radar on the southern Utopia Planitia to identify subsurface dipping reflectors indicative of an ancient prograding shoreline. The reflectors dip unidirectionally with inclinations in the range 6° to 20° and are imaged to a thickness of 10 to 35 m along an uninterrupted 1.3 km northward shoreline-perpendicular traverse. The consistent dip inclinations, absence of dissection by fluvial channels along the extended traverse, and low permittivity of the sediments are consistent with terrestrial coastal deposits-and discount fluvial, aeolian, or magmatic origins favored elsewhere on Mars. The structure, thickness, and length of the section support voluminous supply of onshore sediments into a large body of water, rather than a merely localized and short-lived melt event. Our findings not only provide support for the existence of an ancient Martian ocean in the northern plains but also offer crucial insights into the evolution of the ancient Martian environment.

The triple oxygen isotope composition of sulphate minerals has been used to constrain the evolution of Earths surface environment (e.g., pO2, pCO2 and gross primary productivity) throughout the Proterozoic Eon. This approach presumes the incorporation of atmospheric O2 atoms into riverine sulphate via the oxidative weathering of pyrite. However, this is not borne out in recent geological or modern sulphate records, where an atmospheric signal is imperceptible and where terrestrial pyrite weathering occurs predominantly in bedrock fractures that are physically more removed from atmospheric O2. To better define the transition from a Proterozoic to a modern-like weathering regime, here we present new measurements from twelve marine evaporite basins spanning the Phanerozoic. These data display a step-like transition in the triple oxygen isotope composition of evaporite sulphate during the mid-Paleozoic (420 to 387.7 million years ago). We propose that the evolution of early root systems deepened the locus of pyrite oxidation and reduced the incorporation of O2 into sulphate. Further, the early Devonian proliferation of land plants increased terrestrial organic carbon burial, releasing free oxygen that fueled increased redox recycling of soil-bound iron and resulted in the final rise in pO2 to modern-like levels.

Applying full-waveform methods to image small-scale structures of geophysical interest buried within the Earth requires the computation of the seismic wavefield over large distances compared to the target wavelengths. This represents a considerable computational cost when using state-of-the-art numerical integration of the equations of motion in three-dimensional earth models. Box Tomography is a hybrid method that breaks up the wavefield computation into three parts, only one of which needs to be iterated for each model update, significantly saving computational time. To deploy this method in remote regions containing a fluid-solid boundary, one needs to construct artificial sources that confine the seismic wavefield within a small region that straddles this boundary. The difficulty arises from the need to combine the solid-fluid coupling with a hybrid numerical simulation in this region. Here, we report a reconciliation of different displacement potential expressions used for solving the acoustic wave equation and propose a unified framework for hybrid simulations. This represents a significant step towards applying Box Tomography in arbitrary regions inside the Earth, achieving a thousand-fold computational cost reduction compared to standard approaches without compromising accuracy. We also present examples of benchmarks of the hybrid simulations in the case of target regions at the ocean floor and the core-mantle boundary.

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.

This study investigates the factors influencing runoff response to atmospheric rivers (ARs) over the U.S. West Coast. We focused on runoff time series variations impacted by AR characteristics (e.g., category and frequency) and land preconditions during Northern Hemisphere cool seasons in the period of 1940–2023. Results show that high-category ARs significantly increase local runoff with higher hourly precipitation rates leading to a greater incremental rate and peak runoff. Extreme runoff increases greatly with the AR category with an increase rate up to 12.5 times stronger than non-extreme runoff. Besides the AR category, land preconditions such as soil moisture and snowpack also play crucial roles in modulating runoff response. We found that runoff induced by weak-category ARs is more sensitive to land preconditions than high-category ARs, with high peak runoff occurring when soil is nearly saturated. Additionally, more than 50% of high-peak-runoff events in snow-covered grid cells are associated with rain-on-snow events particularly for the events associated with weaker ARs. Regression analysis reveals that AR precipitation and land preconditions jointly influence runoff, emphasizing the importance of including soil moisture and snowpack levels in AR impact assessments. The study also highlights the intensified runoff response to back-to-back ARs with short intervals, which may become more frequent with climate warming, posing increased flood risks via facilitating wet soil conditions. Our findings have significant implications for AR risk predictions and the development of prediction models for AR-induced runoff.

Rubisco is the primary CO2-fixing enzyme of the biosphere1, yet it has slow kinetics2. The roles of evolution and chemical mechanism in constraining its biochemical function remain debated3,4. Engineering efforts aimed at adjusting the biochemical parameters of rubisco have largely failed5, although recent results indicate that the functional potential of rubisco has a wider scope than previously known6. Here we developed a massively parallel assay, using an engineered Escherichia coli7 in which enzyme activity is coupled to growth, to systematically map the sequence-function landscape of rubisco. Composite assay of more than 99% of single-amino acid mutants versus CO2 concentration enabled inference of enzyme velocity and apparent CO2 affinity parameters for thousands of substitutions. This approach identified many highly conserved positions that tolerate mutation and rare mutations that improve CO2 affinity. These data indicate that non-trivial biochemical changes are readily accessible and that the functional distance between rubiscos from diverse organisms can be traversed, laying the groundwork for further enzyme engineering efforts.

The Red Planet is a magnetic planet. The Martian crust contains strong magnetization from a core dynamo that likely was active during the Noachian period when the surface may have been habitable. The evolution of the dynamo may have played a central role in the evolution of the early atmosphere and the planets transition to the current cold and dry state. However, the nature and history of the dynamo and crustal magnetization are poorly understood given the lack of well-preserved, oriented, ancient samples with geologic context available for laboratory study. Here, we describe how magnetic measurements of returned samples could transform our understanding of six key unknowns about Mars planetary evolution and habitability. Such measurements could i) determine the history of the Martian dynamo fields intensity; ii) determine the history of the Martian dynamo fields direction; iii) test the hypothesis that Mars experienced plate tectonics or true polar wander; iv) constrain the thermal and aqueous alteration history of the samples; v) identify sources of Martian crustal magnetization and vi) characterize sedimentary and magmatic processes on Mars. We discuss how these goals can be achieved using future laboratory analyses of samples acquired by the Perseverance rover.

Specialized or secondary metabolites are small molecules of biological origin, often showing potent biological activities with applications in agriculture, engineering and medicine. Usually, the biosynthesis of these natural products is governed by sets of co-regulated and physically clustered genes known as biosynthetic gene clusters (BGCs). To share information about BGCs in a standardized and machine-readable way, the Minimum Information about a Biosynthetic Gene cluster (MIBiG) data standard and repository was initiated in 2015. Since its conception, MIBiG has been regularly updated to expand data coverage and remain up to date with innovations in natural product research. Here, we describe MIBiG version 4.0, an extensive update to the data repository and the underlying data standard. In a massive community annotation effort, 267 contributors performed 8304 edits, creating 557 new entries and modifying 590 existing entries, resulting in a new total of 3059 curated entries in MIBiG. Particular attention was paid to ensuring high data quality, with automated data validation using a newly developed custom submission portal prototype, paired with a novel peer-reviewing model. MIBiG 4.0 also takes steps towards a rolling release model and a broader involvement of the scientific community. MIBiG 4.0 is accessible online at https://mibig.secondarymetabolites.org/.

Abstract: Many geoscience departments in the United States (US) are working to recruit faculty from underrepresented groups. However, there is little information about how hiring practices are perceived by candidates. Here we address this gap by interviewing 19 geoscientists who identify as an underrepresented race, ethnicity, or gender who recently declined a tenure-track faculty job offer in the US about their faculty job searches, with an emphasis on their decisions to accept or decline an offer. We find that many participants experienced hiring practices inconsistent with existing recommendations to increase faculty diversity, and some participants were subject to uncivilized, even potentially discriminatory, practices. Therefore, we leverage our results to provide actionable recommendations for improving faculty recruitment efforts. We highlight that departments may doubly benefit from improving their culture: in addition to benefiting current members, it may also help with recruitment. Overall, our findings emphasize the need for continued evaluation of faculty hiring practices.

Abstract. Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (re) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed −0.34 ± 0.23, −0.22 ± 0.13, and −0.09 ± 0.11 W m−2, respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS re–LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.