Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discover”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Methane Plume Emissions Associated With Puget Sound Faults in the Cascadia Forearc

Methane gas plumes have been discovered to issue from the seafloor in the Puget Sound estuary. These gas emission sites are co-located over traces of three major fault zones that fracture the entire forearc crust of the Cascadia Subduction Zone. Multibeam and single-beam sonar data from cruises conducted in 2011, 2018, 2019, 2020 and 2021 identified the acoustic signature of over 330 bubble plumes. Dissolved gas from the plumes combines to elevate seawater methane concentrations of the entire Puget Sound estuary. Fluid samples from adjacent terrestrial hot springs and deep-water wells surrounding the estuary contain a helium-3 isotope signature, indicating a deep fluid source located near the underlying Cascadia Subduction Zone. However, Puget Sound seawater emission sites lack chemical isotope signatures and elevated thermal anomalies expected from association with a deep plate-interface reservoir. The existence of vigorous marine methane plumes located in areas of thin sediment cover and above the traces of deeply-penetrating forearc fault zones – but still presenting no thermal or chemical anomalies that are found in many other similar forearc environments and remain an unresolved paradox.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Overview of the Morphology and Chemistry of Diagenetic Features in the Clay-Rich Glen Torridon Unit of Gale Crater, Mars

The clay-rich Glen Torridon region of Gale crater, Mars, was explored between sols 2300 and 3007. Here, we analyzed the diagenetic features observed by Curiosity, including veins, cements, nodules, and nodular bedrock, using the ChemCam, Mastcam, and Mars Hand Lens Imager instruments. We discovered many diagenetic features in Glen Torridon, including dark-toned iron- and manganese-rich veins, magnesium- and fluorine-rich linear features, Ca-sulfate cemented bedrock, manganese-rich nodules, and iron-rich strata. We have characterized the chemistry and morphology of these features, which are most widespread in the higher stratigraphic members in Glen Torridon, and exhibit a wide range of chemistries. These discoveries are strong evidence for multiple generations of fluids from multiple chemical endmembers that likely underwent redox reactions to form some of these features. In a few cases, we may be able to use mineralogy and chemistry to constrain formation conditions of the diagenetic features. For example, the dark-toned veins likely formed in warmer, highly alkaline, and highly reducing conditions, while manganese-rich nodules likely formed in oxidizing and circumneutral conditions. We also hypothesize that an initial enrichment of soluble elements, including fluorine, occurred during hydrothermal alteration early in Gale crater history to account for elemental enrichment in nodules and veins. The presence of redox-active elements, including Fe and Mn, and elements required for life, including P and S, in these fluids is strong evidence for habitability of Gale crater groundwater. Hydrothermal alteration also has interesting implications for prebiotic chemistry during the earliest stages of the crater’s evolution and early Mars.

58 GEOSCIENCES↗
DOE OSTI2021

PreMevE Update: Forecasting Ultra-Relativistic Electrons Inside Earth's Outer Radiation Belt

Energetic electrons inside Earth's Van Allen belts pose a major radiation threat to spaceborne electronics that often play vital roles in modern society. Ultra-relativistic electrons with energies greater than or equal to two megaelectron-volt (MeV) are of particular interest, and thus forecasting these ≥2 MeV electrons has a significant meaning to all space sectors. Here, we update the latest development of the predictive model for MeV electrons in the outer radiation belt. The new version, called PREdictive MEV Electron (PreMevE)-2E, forecasts ultra-relativistic electron flux distributions across the outer belt, with no need for in situ measurements of the trapped MeV electron population except at the geosynchronous orbit (GEO). Model inputs include precipitating electrons observed in low-Earthorbits by NOAA satellites, upstream solar wind speeds and densities from solar wind monitors, as well as ultra-relativistic electrons measured by one Los Alamos GEO satellite. We evaluated 32 supervised machine learning models that fall into four different classes of linear and neural network architectures, and successfully tested ensemble forecasting by using groups of top-performing models. All models are individually trained, validated, and tested by in situ electron data from NASA's Van Allen Probes mission. It is shown that the final ensemble model outperforms individual models at most L-shells, and this PreMevE-2E model can provide 25-h (~1-day) and 50-h (~2-day) forecasts with high mean performance efficiency and correlation values. Our results also suggest that this new model is dominated by nonlinear components at L-shells <~4 for ultra-relativistic electrons, different from the dominance of linear components for 1 MeV electrons as previously discovered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effect of Nucleation Heterogeneity on Mineral Precipitation in Confined Environments

Abstract The formation of new mineral phases in confined environments, especially in porous media, is crucial for various geological processes like mineralization and diagenesis. The nucleation and precipitation of minerals are initiated at the microscale through fluid‐rock interaction, where dissolution of primary phases leads to supersaturated conditions and nucleation and growth of secondary ones. Previous research has focused primarily on either precipitation or nucleation, without fully exploring their combined impact. Our study introduces a computational framework that integrates classical nucleation theory with the micro‐continuum method. We validated our model by comparing with experiments, and discovered that different surface nucleation rate changes the mode of precipitation from a preferential to uniform precipitate textures. Furthermore, our study uncovered that the conventional deterministic precipitation method tends to underestimate the permeability of the porous matrix. In contrast, the new framework significantly improves model accuracy by incorporating preferential precipitation and heterogeneous nucleation.

58 GEOSCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2024

Reduced‐Order Modeling for Linearized Representations of Microphysical Process Rates

Abstract Representing cloud microphysical processes in large scale atmospheric models is challenging because many processes depend on the details of the droplet size distribution (DSD, the spectrum of droplets with different sizes in a cloud). While full or partial statistical moments of droplet size distributions are the typical variables used in bulk models, prognostic moments are limited in their ability to represent microphysical processes across the range of conditions experienced in the atmosphere. Microphysical parameterizations employing prognostic moments are known to suffer from structural uncertainty in their representations of inherently higher dimensional cloud processes, which limit model fidelity and lead to forecasting errors. Here we investigate how data‐driven reduced‐order modeling can be used to learn predictors for microphysical process rates in bulk microphysics schemes in an unsupervised manner from higher dimensional bin distributions. Using simulations characteristic of marine stratiform clouds, we simultaneously learn lower dimensional representations of droplet size distributions and predict the evolution of the microphysical state of the system. Droplet collision‐coalescence, the main process for generating warm rain, is estimated to have an intrinsic dimension of three. This intrinsic dimension provides a lower limit on the number of degrees of freedom needed to accurately represent collision‐coalescence in models. We demonstrate how deep learning based reduced‐order modeling can be used to discover intrinsic coordinates describing the microphysical state of the system, where process rates such as collision‐coalescence are globally linearized. These implicitly learned representations of the DSD retain more information about the DSD than typical moment‐based representations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impacts of Permeability Heterogeneity and Background Flow on Supercritical CO 2 Dissolution in the Deep Subsurface

Motivated by CO 2 capture and sequestration (CCS) design considerations, we consider the coupled effects of permeability heterogeneity and background flow on the dissolution of a supercritical CO 2 lens into an underlying deep, confined aquifer. We present the results of a large-scale Monte Carlo simulation study examining the interaction of background flow rate and three parameters describing multi-Gaussian log-permeability fields: mean, variance, and correlation length. Hundreds of high-resolution simulations were performed using the PFLOTRAN finite volume software to model CO 2 dissolution in a kilometer-scale aquifer over 1,000 years. Predictive dimensionless scaling relationships relating CO 2 dissolution rate to heterogeneity statistics, Rayleigh (Ra) and Péclet (Pe) numbers were developed for both gravitationally dominated free convection to background flow-dominated forced convection regimes. An empirical criterion, Pe = Ra 3/4 , was discovered for regime transition. All simulations converged quickly to a quasi-steady, approximately linear dissolution rate. However, this rate displayed profound variability between permeability field realizations sharing the same heterogeneity statistics, even under mild permeability heterogeneity. In general, increased heterogeneity was associated with a lower mean and higher variance of dissolution rate, undesirable from a CCS design perspective. The relationship between dissolution rate and background flow was found to be complex and nonlinear. Dimensionless scaling relationships were uncovered for a number of special cases. Results call into question the validity of the Boussinesq approximation in the context of modest-to-high background flow rates and the general applicability of numerical simulations without background flow.

42 ENGINEERING↗
DOE OSTI2024

DASEventNet: AI‐Based Microseismic Detection on Distributed Acoustic Sensing Data From the Utah FORGE Well 16A (78)‐32 Hydraulic Stimulation

Abstract Distributed acoustic sensing (DAS) has emerged as a promising seismic technology for monitoring microearthquakes (MEQs) with high spatial resolution. Efficient algorithms are needed for processing large DAS data volumes. This study introduces a deep learning (DL) model based on a Residual Convolutional Neural Network (ResNet) for detecting MEQs using DAS data, named as DASEventNet. The test data were collected from the Utah FORGE 16A (78)‐32 hydraulic stimulation experiments conducted in April 2022. The DASEventNet model achieves a remarkable accuracy of 100% when discriminating MEQs from noise in the raw test set of 260 examples. Surprisingly, the model identified weak MEQ signatures that have been manually categorized as noise. The decision‐making process with the model is decoded by the classic activation map, which illuminates learning features of the DASEventNet model. These features provide clear illustrations of weak MEQs and varied noise types. Finally, we apply the trained model to the entire period (∼7 days) of continuous DAS recordings and find that it discovers >5,700 new MEQs, previously unregistered in the public Silixa DAS catalog. The DASEventNet model significantly outperforms the traditional seismic method Short‐Term Average/Long‐Term Average (STA/LTA), which detected only 1,307 MEQs. The DASEventNet detection threshold is M w −1.80 compared to the minimum magnitude of M w −1.14 detected by STA/LTA. The spatiotemporal distribution of the newly identified MEQs defines an extensive stimulation zone and more accurately characterizes fracture geometry. Our results highlight the potential of DL for long‐term, real‐time microseismic monitoring that can improve enhanced geothermal systems and other activities that include subsurface hydraulic fracturing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

A Climatology and Life‐Cycle Characteristics of Atmospheric Fronts and Their Associated Precipitation

Abstract Atmospheric fronts are one of the main sources of mid‐latitude variability. We employ a novel method for identifying and tracking fronts and frontal precipitation. Thermal and dynamical variables are used to identify fronts as areal objects in space, which are tracked in time using the open‐source TempestExtremes software package. Precipitation objects are co‐located to identify frontal precipitation. The method is subjected to validation and sensitivity tests using manually curated data from the National Weather Service. Climatologies of fronts and frontal precipitation are computed from reanalysis and observations; fronts are present upwards of 14% of the time in the storm tracks, and represent the majority (up to 90%) of total and extreme precipitation. Novel aspects of the method are showcased through the lifetime characteristics of fronts across North America. Three sets of warm and cold fronts were discovered, and their duration, distance‐traveled, and translation velocity are examined. Plain Language Summary Mid‐latitude low‐pressure systems and weather fronts are important for our day‐to‐day experience of weather events, particularly in the mid‐latitudes. This work makes use of standardized atmospheric data and creates a method of automatically tracking these important atmospheric features and their precipitation to quantify their relative role in global precipitation. Weather fronts are persistent in the mid‐latitudes and are associated with the majority of precipitation–particularly the most intense precipitation. Trajectories of fronts over North America are categorized to create a set of archetypal fronts that occur in that region. The differences between these types of fronts are characterized. Key Points An automated, efficient, and skillful frontal detection algorithm is developed and validated Fronts contribute a larger fraction of extreme precipitation than all precipitation in mid‐latitude storm tracks Fronts across North America have substantial variation in characteristics depending on their origin location

extratropical cyclone↗
DOE OSTI2025

Survey and Modeling of Windblown Ejecta Deposits on Venus

Venus' thick atmosphere rotates in the same direction as the solid body, but ∼60 times faster. This atmospheric superrotation has produced dozens of windblown ejecta deposits (“parabolas”) on the surface of Venus. The formation and modification of parabolas is an interplay between impacts, aeolian modification, and atmospheric dynamics. We conducted a survey to explore the nature of these sedimentary surface features. First, we observe trends in parabolas' morphology that shed light on how they are deposited and gradated. Changes in the size and radar albedo of parabolas are likely linked to the height and density (respectively) of ejecta plumes at time of formation. Next, we discovered that parabolas show orientations inconsistent with present atmospheric dynamics. This discrepancy may record a change in these dynamics or geologically recent true polar wander at a rate of ∼1° Myr −1 , which is similar to that observed on Earth over the past century. These results highlight how overlapping observations at different radar wavelengths provide important insights into the history and character of geologic processes on Venus. Overall, atmospheric superrotation has probably persisted for at least the age of Venus' surface.

Venus↗
DOE OSTI2025

Causal Directions Matter: How Environmental Factors Drive Convective Cloud Detrainment Heights

This study investigates how environmental factors influence the level of maximum detrainment (LMD) in deep convective clouds. Through a novel application of the Linear Non‐Gaussian Acyclic Model (LiNGAM), we discover causal structures between environmental variables and LMD, observed at six tropical sites operated by the Atmospheric Radiation Measurement (ARM) user facility. LiNGAM effectively identifies causal directions among variables of interest, revealing robust relationships such as those among the lifting condensation level (LCL), level of free convection (LFC), and convective inhibition (CIN), aligning with prior knowledge. Relative humidity is shown to directly influence LMD; however, this relationship exhibits strong nonlinearity and becomes difficult to detect when the contrast between oceanic and continental environments is excluded from the analysis. This study highlights the importance of establishing causal relationships before performing statistical inference.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ubiquity and Causes of Soil Water Preferential Flow Across 17 Ecoregions

Abstract Preferential flow (PF) in soil causes the rapid transport of water, nutrients, and contaminants into the subsurface, influencing groundwater recharge and streamflow. Data scarcity has hindered the quantification of PF occurrence and the identification of its drivers across diverse ecoregions. We address this gap by analyzing high‐frequency, multi‐depth soil moisture data across 17 ecoregions in the USA, using ∼1,500 sensors at 40 sites. We discovered that PF is widespread, with sites experiencing PF in up to 60% of rainfall events ≥2 mm. Multiple approaches consistently show that PF is more likely to occur with increased peak rainfall intensity, finer textured material, low soil moisture variability, humid climate, and higher net primary productivity. This suggests that PF patterns could shift with projected climate changes, increasing uncertainty in predictions of groundwater recharge, water quality, and streamflow generation. Plain Language Summary Water can bypass part of the soil's matrix through a process called preferential flow (PF). This quick transport of water through the soil brings with it nutrients and contaminants and eventually makes it to groundwater and streams. To ensure ample amounts of good quality groundwater and surface water we need to understand when and where PF occurs. We inferred when PF occurred at 40 different sites across 17 ecoregions in the USA using soil moisture and rainfall data. We found that PF happened at all sites and in up to 60% of rainfall events ≥2 mm. Preferential flow was most likely at sites with high rainfall intensities, high clay content in soils, low variability in soil moisture, and high vegetation productivity. As rainfall intensities are predicted to increase due to climate change and vegetation becomes more productive, PF occurrence becomes more important for predicting groundwater recharge, water quality, and streamflow generation. Key Points Preferential flow (PF) is ubiquitous across the USA and occurs in up to 60% of all rainfall events ≥2 mm Rainfall intensity, soil texture, and antecedent soil moisture emerge as critical in generating PF across diverse ecoregions Two different PF detection approaches show similar relationships between key drivers and occurrence of PF

Li, Bonan↗
DOE OSTI2025

Quantifying the Effect of Pore‐Size Dependent Wettability on Relative Permeabilities Using Lattice Boltzmann Simulation

Relative permeability is a crucial two-phase property in porous media that can be significantly impacted by wettability conditions. While traditional research has predominantly examined homogeneous wettability, this work explores the less studied pore-size dependent (PSD) wettability, featured by a pore-size dependent wettability distribution. Leveraging high-fidelity Lattice Boltzmann simulations on CT-scanned porous samples, we demonstrate how PSD wettability would impact relative permeability at the pore scale. Our findings reveal that the deviation of relative permeability from the homogeneous wettabilitity induced by PSD wettability can be 5%–20%. The deviation of relative permeability curves increases as the spanning range of the contact angle increases. We also find that this impact is less pronounced as the capillary number increases. By adopting a pore-size-dependent contact angle relationship, our approach provides a more accurate and nuanced understanding of how PSD wettability would impact two-phase flow. These findings discovered at the pore scale may also provide valuable insights on relative permeability at the core to reservoir scales.

25 ENERGY STORAGE↗
DOE OSTI2025

The Influence of African Easterly Waves on Atlantic Tropical Cyclone Tracks and Landfall in Large Ensembles

African easterly waves (AEWs) are an important precursor or “seed” for Atlantic tropical cyclones (TCs), with 60%–80% of major hurricanes observed to originate from AEWs. However, climate model simulations indicate that AEWs are not necessary to maintain annual Atlantic TC frequency. Furthermore, small ensembles suggest that AEWs may impact the spatial distribution and landfall of Atlantic TCs. Here, we investigated the influence of AEWs on the spatial distribution of Atlantic TC tracks and landfall using 50-member ensembles of TC-permitting regional model simulations for five hurricane seasons characterized by different levels of TC activity. The control simulations are seasonal hindcasts in which AEWs were prescribed through the eastern lateral boundary condition using reanalysis. In the experiments, we suppressed AEWs by applying a 2–10 day filter to the eastern lateral boundary condition. In response to AEW suppression, we discovered statistically significant increases in Atlantic TC frequency (10%–26%) and landfall (16%–44%), a westward shift in TC genesis location and tracks with increased landfall over the Caribbean Islands, southwestern Gulf Coast, and southeastern US coast, and increases in mid-tropospheric relative humidity in the main development region. In addition, we evaluated TC genesis mechanisms in the absence of AEWs and found evidence that the intertropical convergence zone (ITCZ) intensified and extended northward, resulting in ITCZ wave-breaking that shed vortices which served as TC seeds. By uncovering the connections between TC seed types and the likelihood of TC landfall, this research can provide Atlantic coastal and island communities with useful information to prepare for TC impacts.

African easterly waves↗
DOE OSTI2025

Learning Physically Interpretable Atmospheric Models From Data With WSINDy

The multiscale and turbulent nature of Earth's atmosphere has historically rendered accurate weather modeling a hard problem. Recently, there has been an explosion of interest surrounding data-driven approaches to weather modeling, which in many cases show improved forecasting accuracy and computational efficiency when compared to traditional methods. However, many of the current data-driven approaches employ highly parameterized neural networks, often resulting in uninterpretable models and limited gains in scientific understanding. In this work, we address the interpretability problem by explicitly discovering partial differential equations governing atmospheric phenomena, identifying symbolic mathematical models with direct physical interpretations. The purpose of this paper is to demonstrate that, in particular, the weak-form sparse identification of nonlinear dynamics (WSINDy) algorithm can learn effective atmospheric models from both simulated and assimilated data. Our approach adapts the standard WSINDy algorithm to work with high-dimensional fluid data of arbitrary spatial dimension.

58 GEOSCIENCES↗
DOE OSTI2021

CRL4A DTL degrades DNA-PKcs to modulate NHEJ repair and induce genomic instability and subsequent malignant transformation

Genomic instability induced by DNA damage and improper DNA damage repair is one of the main causes of malignant transformation and tumorigenesis. DNA double strand breaks (DSBs) are the most detrimental form of DNA damage, and nonhomologous end-joining (NHEJ) mechanisms play dominant and priority roles in initiating DSB repair. A well-studied oncogene, the ubiquitin ligase Cullin 4A (CUL4A), is reported to be recruited to DSB sites in genomic DNA, but whether it regulates NHEJ mechanisms of DSB repair is unclear. Here, we discovered that the CUL4A-DTL ligase complex targeted the DNA-PKcs protein in the NHEJ repair pathway for nuclear degradation. Overexpression of either CUL4A or DTL reduced NHEJ repair efficiency and subsequently increased the accumulation of DSBs. Moreover, we demonstrated that overexpression of either CUL4A or DTL in normal cells led to genomic instability and malignant proliferation. Consistent with the in vitro findings, in human precancerous lesions, CUL4A expression gradually increased with increasing malignant tendency and was negatively correlated with DNA-PKcs and positively correlated with γ-H2AX expression. Collectively, this study provided strong evidence that the CUL4A-DTL axis increases genomic instability and enhances the subsequent malignant transformation of normal cells by inhibiting NHEJ repair. These results also suggested that CUL4A may be a prognostic marker of precancerous lesions and a potential therapeutic target in cancer.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Experimental evolution reveals nitrate tolerance mechanisms in Desulfovibrio vulgaris

Elevated nitrate in the environment inhibits sulfate reduction by important microorganisms of sulfate-reducing bacteria (SRB). Several SRB may respire nitrate to survive under elevated nitrate, but how SRB that lack nitrate reductase survive to elevated nitrate remains elusive. To understand nitrate adaptation mechanisms, we evolved 12 populations of a model SRB (i.e., Desulfovibrio vulgaris Hildenborough, DvH) under elevated NaNO 3 for 1000 generations, analyzed growth and acquired mutations, and linked their genotypes with phenotypes. Nitrate-evolved (EN) populations significantly (p < 0.05) increased nitrate tolerance, and whole-genome resequencing identified 119 new mutations in 44 genes of 12 EN populations, among which six functional gene groups were discovered with high mutation frequencies at the population level. In this work, we observed a high frequency of nonsense or frameshift mutations in nitrosative stress response genes (NSR: DVU2543, DVU2547, and DVU2548), nitrogen regulatory protein C family genes (NRC: DVU2394-2396, DVU2402, and DVU2405), and nitrate cluster (DVU0246-0249 and DVU0251). Mutagenesis analysis confirmed that loss-of-functions of NRC and NSR increased nitrate tolerance. Also, functional gene groups involved in fatty acid synthesis, iron regulation, and two-component system (LytR/LytS) known to be responsive to multiple stresses, had a high frequency of missense mutations. Additionally, mutations in those gene groups could increase nitrate tolerance through regulating energy metabolism, barring entry of nitrate into cells, altering cell membrane characteristics, or conferring growth advantages at the stationary phase. This study advances our understanding of nitrate tolerance mechanisms and has important implications for linking genotypes with phenotypes in DvH.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Abstract Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR− and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

59 BASIC BIOLOGICAL SCIENCES↗