Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “yield estimation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Enhanced detection of paramagnetic fluorine-19 magnetic resonance imaging agents using zero echo time sequence and compressed sensing

Fluorine-19 ( 19 F) magnetic resonance imaging (MRI) is an emerging technique offering specific detection of labeled cells in vivo. Lengthy acquisition times and modest signal-to-noise ratio (SNR) makes three-dimensional spin-density–weighted 19 F imaging challenging. Recent advances in tracer paramagnetic metallo-perfluorocarbon (MPFC) nanoemulsion probes have shown multifold SNR improvements due to an accelerated 19 F T 1 relaxation rate and a commensurate gain in imaging speed and averages. However, 19 F T 2 -reduction and increased linewidth limit the amount of metal additive in MPFC probes, thus constraining the ultimate SNR. To overcome these barriers, we describe a compressed sampling (CS) scheme, implemented using a “zero” echo time (ZTE) sequence, with data reconstructed via a sparsity-promoting algorithm. Our CS-ZTE scheme acquires k-space data using an undersampled spherical radial pattern and signal averaging. Image reconstruction employs off-the-shelf sparse solvers to solve a joint total variation and -norm regularized least square problem. To evaluate CS-ZTE, we performed simulations and acquired 19 F MRI data at 11.7 T in phantoms and mice receiving MPFC-labeled dendritic cells. For MPFC-labeled cells in vivo, we show SNR gains of ~6.3 × with 8-fold undersampling. We show that this enhancement is due to three mechanisms including undersampling and commensurate increase in signal averaging in a fixed scan time, denoising attributes from the CS algorithm, and paramagnetic reduction of T 1 . Importantly, 19 F image intensity analyses yield accurate estimates of absolute quantification of 19 F spins. Overall, the CS-ZTE method using MPFC probes achieves ultrafast imaging, a substantial boost in detection sensitivity, accurate 19 F spin quantification, and minimal image artifacts.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Validation of material models for puncture of 7075-T651 aluminum plate

Plate puncture simulations are challenging computational tasks that require advanced material models including high strain rate and thermal-mechanical effects on both deformation and failure, plus finite element techniques capable of representing large deformations and material failure. The focus of this work is on the material issues, which require large sets of experiments, flexible material models and challenging calibration procedures. In this study, we consider the puncture of 12.7 mm thick, 7075-T651 aluminum alloy plates by a cylindrical punch with a hemispherical nose and diameter of 12.7 mm. The plasticity and ductile failure models were isotropic with calibration data obtained from uniaxial tension tests at different temperatures and strain rates plus quasi-static notched tension tests and shear-dominated tests described here. Sixteen puncture experiments were conducted to identify the threshold penetration energy, mode of puncture and punch acceleration during impact, The punch was mounted on a 139 kg mass and dropped on the plates with different impact speeds. Since the mass was the same in all tests, the quantity of interest was the impact speed. The axis and velocity of the punch were perpendicular to the plate surface. The mean threshold punch speed was 3.05 m/s, and the mode of failure was plugging by thermal-mechanical shear banding accompanied by scabbing fragments. Application of the material models in simulations of the tests yielded accurate estimates of the threshold puncture speed and of the mode of failure. Time histories of the punch acceleration compared well between simulation and test. Remarkably, the success of the simulations occurred in spite of even the smallest element used being larger than the width of the shear bands.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A quantitative risk assessment framework for fault reactivation in underground hydrogen storage: Coupled simulation and deep learning approach

Underground hydrogen storage (UHS) is emerging as a critical solution for large-scale energy storage. However, like all subsurface fluid injection activities, UHS poses the risk of injection-induced fault reactivation. Accurate risk assessment is essential to ensuring the safety and efficiency of UHS operations. This study presents the development of deep-learning surrogate models for fault reactivation prediction in UHS, trained on a comprehensive database of fully coupled fluid flow-geomechanics simulations. Our findings reveal that analytical models often yield unreliable estimates, with errors up to 54% in the allowable injection pressure, potentially leading to a 40% reduction in UHS operational capacity. The developed surrogate models were incorporated into a quantitative risk assessment (QRA) framework, enabling probabilistic evaluation of fault reactivation risk while accounting for uncertainties in the input variables. Site-specific features, such as horizontal stress gradients, fault’s dip and strike angles, and operational parameters like bottom-hole injection pressure and well-fault distance, were identified as the primary drivers of fault reactivation across various stress regimes. Whereas other hydraulic, geological, and poroelastic reservoir properties were found to have a secondary impact. Notably, we observed that the risk of fault reactivation for a critically oriented fault with a static friction coefficient greater than 0.55 remains below 10% in a normal faulting stress regime. However, the risk significantly increases as the stress regime transitions from normal to strike-slip and ultimately to reverse faulting conditions. These findings underscore the importance of rigorous site characterization and comprehensive QRA evaluations to optimize UHS performance and minimize geomechanical risks.

25 ENERGY STORAGE↗
DOE OSTI2024

Informing solar blind radioluminescence imaging through a calibrated spectrum

While direct radiation detection methods offer great insight into the origin of ionizing particles and photons, their use to locate contaminated areas or concealed radioactive sources can lead to undue exposure of personnel and equipment to ionizing radiation or the potential for contamination. These same sources induce ultraviolet (UV) optical photon fluorescence in air – a process referred to as radioluminescence – that may be imaged from low dose regions over larger attenuation lengths than ionizing radiation. However, most optical detection methods are limited to low lighting conditions to image the more abundant ultraviolet-A (UV-A) photons. To extend this capability to room light or daytime conditions, the solar blind region (Ultraviolet-C (UV-C), <280 nm) can be tapped. Though the emission yield of UV-C photons is roughly two orders of magnitude lower than that in the UV-A regime, the UV-C offers dramatic improvements in signal-to-noise ratios under bright lighting conditions due to decreased background interferences. The yield of specific UV-C lines, if present in the literature at all, varies widely, which has a large impact in modeling and analyzing standoff UV-C measurement scenarios. Thus, we have captured improved radioluminescence spectra over 250–400 nm and identified observed emission peaks in ambient air. Many of the UV-C photons produced by ionizing radiation excitation result from high-energy, molecular nitrogen Gaydon-Herman transitions which have had limited study to date for this application. Relating these findings to published UV-A yields, we estimate emissions between 0.15–0.19 photons/MeV over 250–280 nm. Additionally, we also use a commercial corona-discharge imaging camera to demonstrate outdoor UV-C radiation mapping of alpha and gamma emitters from 50 and 75 m standoffs, respectively. The imaged “counts” are compared to optically modelled values and show the same trend over distance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Atomic masses with machine learning for the astrophysical r process

The astrophysical r process plays a vital role in the production of heavy elements. Modeling of the r process is sensitive to masses and further requires knowledge of masses beyond current experimental reach. Therefore, simulations of the r process offer a unique test bed for predicting mass extrapolations. We take a Machine-Learning (ML) approach to model the masses across the entire chart of nuclides. For the first time, we simulate r-process nucleosynthesis with a physics-based ML mass model. We compare simulated abundances to solar data in order to evaluate the model's performance far from stability. The resulting r-process abundances up to thorium and uranium qualitatively match those of the observed solar system abundance pattern, with the characteristic peaks well positioned. We propagate the mass uncertainties obtained from the ML model to r-process abundance yields to estimate an uncertainty band associated with our approach. The size of the uncertainty band is approximately one order of magnitude which aligns with the uncertainty reported using alternative techniques.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

pvlib iotools—Open-source Python functions for seamless access to solar irradiance data

Access to accurate solar resource data is critical for numerous applications, including estimating the yield of solar energy systems, developing radiation models, and validating irradiance datasets. However, lack of standardization in data formats and access interfaces across providers constitutes a major barrier to entry for new users. pvlib python’s iotools subpackage aims to solve this issue by providing standardized Python functions for reading local files and retrieving data from external providers. All functions follow a uniform pattern and return convenient data outputs, allowing users to seamlessly switch between data providers and explore alternative datasets. The pvlib package is community-developed on GitHub: https://github.com/pvlib/pvlib-python. As of pvlib python version 0.9.5, the iotools subpackage supports 12 different datasets, including ground measurement, reanalysis, and satellite-derived irradiance data. The supported ground measurement networks include the Baseline Surface Radiation Network (BSRN), NREL MIDC, SRML, SOLRAD, SURFRAD, and the US Climate Reference Network (CRN). Additionally, satellite-derived and reanalysis irradiance data from the following sources are supported: PVGIS (SARAH & ERA5), NSRDB PSM3, and CAMS Radiation Service (including McClear clear-sky irradiance).

14 SOLAR ENERGY↗
DOE OSTI2023

Bullying among children with heart conditions, National Survey of Children’s Health, 2018–2020

Abstract Children with chronic illnesses report being bullied by peers, yet little is known about bullying among children with heart conditions. Using 2018–2020 National Survey of Children’s Health data, the prevalence and frequency of being bullied in the past year (never; annually or monthly; weekly or daily) were compared between children aged 6–17 years with and without heart conditions. Among children with heart conditions, associations between demographic and health characteristics and being bullied, and prevalence of diagnosed anxiety or depression by bullying status were examined. Differences were assessed with chi-square tests and multivariable logistic regression using predicted marginals to produce adjusted prevalence ratios and 95% confidence intervals. Weights yielded national estimates. Of 69,428 children, 2.2% had heart conditions. Children with heart conditions, compared to those without, were more likely to be bullied (56.3% and 43.3% respectively; adjusted prevalence ratio [95% confidence interval] = 1.3 [1.2, 1.4]) and bullied more frequently (weekly or daily = 11.2% and 5.3%; p < 0.001). Among children with heart conditions, characteristics associated with greater odds of weekly or daily bullying included ages 9–11 years compared to 15–17 years (3.4 [2.0, 5.7]), other genetic or inherited condition (1.7 [1.0, 3.0]), ever overweight (1.7 [1.0, 2.8]), and a functional limitation (4.8 [2.7, 8.5]). Children with heart conditions who were bullied, compared to never, more commonly had anxiety (40.1%, 25.9%, and 12.8%, respectively) and depression (18.0%, 9.3%, and 4.7%; p < 0.01 for both). Findings highlight the social and psychological needs of children with heart conditions.

Cardiovascular System & Cardiology↗
DOE OSTI2025

Czochralski Growth and Characterization of a Compositionally Complex Rare Earth Aluminum Garnet Scintillator: (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce

Compositionally complex oxides have garnered increasing interest for their enhanced phase stability and tunable functional properties, yet their development as bulk single crystal scintillators remains limited. Herein, we report the Czochralski growth and characterization of (Gd 1/4 Y 1/4 Tb 1/4 Lu 1/4 ) 3 Al 5 O 12 :Ce (GYTLAG), a compositionally complex garnet incorporating four dodecahedrally coordinated principal rare earth elements. The garnet phase was confirmed by powder and single crystal X-ray diffraction, and macroscopic defects are described. X-ray absorption near-edge structure measurements confirm the 3+ oxidation state of all rare earths and support their occupation of the same crystallographic site; white line intensity variations correlate with the anticipated segregation behavior. Elemental segregation is quantified by SEM/EDS and ICP-OES, and a linear trend was established between the segregation coefficient and the difference between each rare earth’s ionic radius (r) and the average ionic radius (AIR) of the dodecahedral site. This trend offers a predictive framework for compositional control in future REAG crystals grown by the Czochralski method. Photoluminescence and radioluminescence measurements reveal both Ce 3+ and Tb 3+ emission. Scintillation pulses exhibit four-component decay with dominant ~230 µs and ~1.2 ms components, and the light yield is estimated to be up to 43,000 ph/MeV under 137Cs γ-ray excitation. GYTLAG also demonstrates a strong radioluminescence efficiency and 50% lower afterglow at 20 ms compared to a LuAG:Ce reference, underscoring its promise for scintillator applications.

Compositionally complex oxide, high entropy oxide,↗
DOE OSTI2025

Process Feasibility Analysis of Waste Biomass Valorization to Biochar and Bio-Oil via Slow and Fast Pyrolysis

The United States has abundant biomass and waste feedstock to support the nation's energy addition and affordability targets. Pyrolysis, a thermochemical conversion process, decomposes lignocellulosic feedstocks into liquid, solid, and gaseous fuels that can contribute to the domestic production of biofuels, biopower, and bioproducts. Growing private sector interest in this technology is a key motivation for this comprehensive techno-economic process modeling analysis of a respective biorefinery that includes feedstock preprocessing, slow and fast pyrolysis, and product separation to bio-oil, biochar, and syngas hydrocarbons. Results show that biochar from slow pyrolysis could achieve minimum selling prices (MSPs) of $\$$188-$\$$260/t, competitive with reported market values, while bio-oil from fast pyrolysis is estimated to yield MSPs of $\$$6.49-$\$$9.68/GGE, approximately twice conventional fuel benchmarks. Sensitivity analysis identifies feedstock cost, product yield, and scale as primary cost drivers, while scenarios involving biochar carbon credits and high value applications may substantially improve economics. Overall, these results suggest that continued innovation in feedstock logistics, process integration, and market development will be critical to achieving economically viable and scalable bioproducts.

09 BIOMASS FUELS↗
DOE OSTI2020

The Ground State Electronic Energy of Benzene

We report here on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground-state energy of the benzene molecule in a standard correlation-consistent basis set of double-ζ quality. As a broad international endeavor, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around -863 mE H . However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE H ), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Toward a Unified Kinetic Model of Nitrogenase Catalysis

The microbial enzyme nitrogenase catalyzes the MgATP-dependent reduction of N 2 to 2NH 3 , a transformation central to the global nitrogen cycle. While the canonical Thorneley−Lowe (TL) kinetic model has long served as a mechanistic framework, it does not incorporate several recent insights. Here, we present an updated kinetic model for Monitrogenase that incorporates these new findings. A significant insight is that electron transfer (ET) from the reduced Fe protein to the FeMo-cofactor is gated by MgATP-dependent conformational transitions and can be described as a probabilistic event that is dependent on the ligand bound to the active-site metallocofactor. The updated kinetic model quantitatively reproduces steady-state product formation rates across a broad range of experimental conditions, yielding revised estimates for key rate constants. It is demonstrated that under N 2 turnover, the probability of productive ET to the active site decreases by ∼60%, resulting in a significant fraction of Fe protein cycles that are unproductive for electron delivery. This mechanistic feature explains the observed rate limitation in N 2 reduction and implies a revised minimum energetic cost of approximately 25 MgATP per N 2 reduced. Integrating these new features into the revised kinetic model provides a more complete and usable foundation for understanding nitrogenase catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Episodic Magma Hammers for the 15 January 2022 Cataclysmic Eruption of Hunga Tonga‐Hunga Ha'apai

Understanding the forces and magma system dynamics on timescales of seconds to minutes remains challenging. In the January 2022 phreatoplinian Hunga Tonga-Hunga Ha'apai eruption, four remarkably similar seismic subevents within a 5-min interval occurred during the intensifying early eruptive phase. The subevents are similar in waveforms and durations (~25 s each). Each subevent begins with an unusual negative P-wave polarity which is inferred, using full-wave seismic modeling, to be caused by an upward single-force mechanism at the volcano created by a magma hammer likely in response to magma flow blockage/constriction during the early part of the eruption as discharge rapidly increased over orders of magnitude with concomitant conduit geometry evolution and instability. Our proposed episodic magma hammer model is consistent with thermodynamic and phase properties of the magmatic mixture, and yields an estimate of conduit mass flow in agreement with vent discharge rates derived from satellite imagery of plume heights.

58 GEOSCIENCES↗
DOE OSTI2023

Impacts of a Cascadia Subduction Zone Earthquake on Water Levels and Wetlands of the Lower Columbia River and Estuary

Abstract Subsidence after a subduction zone earthquake can cause major changes in estuarine bathymetry. Here, we quantify the impacts of earthquake‐induced subsidence on hydrodynamics and habitat distributions in a major system, the lower Columbia River Estuary, using a hydrodynamic and habitat model. Model results indicate that coseismic subsidence increases tidal range, with the smallest changes at the coast and a maximum increase of ∼10% in a region of topographic convergence. All modeled scenarios reduce intertidal habitat by 24%–25% and shifts ∼93% of estuarine wetlands to lower‐elevation habitat bands. Incorporating dynamic effects of tidal change from subsidence yields higher estimates of remaining habitat by multiples of 0–3.7, dependent on the habitat type. The persistent tidal change and chronic habitat disturbance after an earthquake poses strong challenges for estuarine management and wetland restoration planning, particularly when coupled with future sea‐level rise effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Assessing the Effect of a Deep‐Rooted Grass on Belowground Carbon Storage in Cultivated Land: Insights From a Multi‐Site US Study

Agriculture depletes soil organic carbon (SOC), partly due to the exclusion of deep-rooted perennials. Reintroducing deep-rooted perennials to cultivated land may help to mitigate SOC loss. We quantified the effect of deep roots on SOC by comparing 8 to 30 year-old stands of switchgrass ( Panicum virgatum L .) with paired annual row crop fields at 12 sites across the central and eastern USA. We hypothesized that switchgrass would store more root C and SOC than neighboring shallow-rooted annual crops, and that these effects would extend deeper than 30 cm. We also evaluated whether switchgrass stimulates decomposition of SOC at depth using radiocarbon ( 14 C). Finally, we explored whether the effect of switchgrass on SOC is moderated by soil chemical and physical properties. While the effect of switchgrass on SOC in the surface 100 cm was positive at most sites, the average effect was not statistically significant (difference in SOC = 0.6 kg C m −2 [95% CI −0.8 to +1.9 kg C m −2 ]). By contrast, we found that root C was consistently more abundant under switchgrass, yielding an estimated additional 0.6 kg C m −2 in the surface 100 cm of soil [95% CI +0.5 to +0.7 kg C m −2 ]. 14 C measurements suggested that root C inputs were adding to existing SOC without stimulating decomposition. The effect of switchgrass on belowground C was not strongly related to any of the soil properties that we evaluated. Our observations show that root C can contribute substantially to belowground C stocks when deep-rooted perennials replace shallow-rooted crops.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Coupled Climate Simulations With E3SM‐MMF

Simulations of the recent historical period from 1950 to 2014 are conducted with E3SM‐MMF, which uses an embedded 2D cloud resolving model that runs efficiently on GPUs in place of traditional parameterizations for cloud and turbulence. Analysis of the climate and variability reveal several aspects where E3SM‐MMF produces smaller biases compared to E3SMv2, including better agreement with the observed evolution of global mean surface temperature, although the representation of ENSO is too weak and fast. Three idealized abrupt CO 2 experiments were also conducted to assess climate sensitivity and feedbacks. These yield three estimates of effective climate sensitivity (4.38, 5.21, and 6.06 K), with a corresponding spread in the shortwave cloud feedbacks. These estimates are on the higher end of sensitivity estimates from CMIP ensembles, and the spread indicates substantial state‐dependent feedbacks. These results demonstrate how multiscale modeling framework (MMF) models can be used for climate relevant experiments and projections by leveraging modern GPU enabled computational platforms. The unique qualities of E3SM‐MMF shown in previous literature are largely still present, but various instances of reduced biases suggest that MMF models have utility in improving future projections.

E3SM↗
DOE OSTI2023

The learnability of Pauli noise

Recently, several quantum benchmarking algorithms have been developed to characterize noisy quantum gates on today’s quantum devices. A fundamental issue in benchmarking is that not everything about quantum noise is learnable due to the existence of gauge freedom, leaving open the question what information is learnable and what is not, which is unclear even for a single CNOT gate. Here we give a precise characterization of the learnability of Pauli noise channels attached to Clifford gates using graph theoretical tools. Our results reveal the optimality of cycle benchmarking in the sense that it can extract all learnable information about Pauli noise. We experimentally demonstrate noise characterization of IBM’s CNOT gate up to 2 unlearnable degrees of freedom, for which we obtain bounds using physical constraints. In addition, we show that an attempt to extract unlearnable information by ignoring state preparation noise yields unphysical estimates, which is used to lower bound the state preparation noise.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

59 BASIC BIOLOGICAL SCIENCES↗