Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model communities”

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 523 records · Page 29

Assessing Clouds Using Satellite Observations Through Three Generations of Global Atmosphere Models

Abstract Clouds are parameterized in climate models using quantities on the model grid‐scale to approximate the cloud cover and impact on radiation. Because of the complexity of processes involved with clouds, these parameterizations are one of the key challenges in climate modeling. Differences in parameterizations of clouds are among the main contributors to the spread in climate sensitivity across models. In this work, the clouds in three generations of an atmosphere model lineage are evaluated against satellite observations. Satellite simulators are used within the model to provide an appropriate comparison with individual satellite products. In some respects, especially the top‐of‐atmosphere cloud radiative effect, the models show generational improvements. The most recent generation, represented by two distinct branches of development, exhibits some regional regressions in the cloud representation; in particular the southern ocean shows a positive bias in cloud cover. The two branches of model development show how choices during model development, both structural and parametric, lead to different cloud climatologies. Several evaluation strategies are used to quantify the spatial errors in terms of the large‐scale circulation and the cloud structure. The Earth mover's distance is proposed as a useful error metric for the passive satellite data products that provide cloud‐top pressure‐optical depth histograms. The cloud errors identified here may contribute to the high climate sensitivity in the Community Earth System Model, version 2 and in the Energy Exascale Earth System Model, version 1.

54 ENVIRONMENTAL SCIENCES↗

Modeling, Load Profile Validation, and Assessment of Solar-Rooftop Energy Potential for Low-and-Moderate-Income Communities in the Caribbean

This document presents the modeling of load profile consumption for Low-and-Moderate-Income (LMI) communities in the Caribbean Islands, as well as an assessment of the solar-rooftop energy potential. In this work, real data, together with synthetic and electricity bill data, were collected to validate and improve the load profile models. The solar-rooftop energy potential was obtained through a National Renewable Energy Laboratory (NREL) software called the PVWatts calculator, and mathematical analysis. The analysis of rooftop solar energy potential was conducted to enable the minimum size of solar power systems to fit the energy demand in the community. The results obtained allow estimation of the capacity of the energy system for each house or an entire community.

14 SOLAR ENERGY↗

Persistent ocean anomalies as a response to Northern Hemisphere heating induced by biomass burning variability

Biomass burning aerosol (BBA) emissions in the Coupled Model Intercomparison Project Phase 6 (CMIP6) historical forcing fields have enhanced temporal variability during the years 1997–2014 compared to earlier periods. Recent studies document that the corresponding inhomogeneous shortwave forcing over this period can cause changes in clouds, permafrost, and soil moisture, which contribute to a net terrestrial Northern Hemisphere warming relative to earlier periods. Here, we investigate the ocean response to the hemispherically asymmetric warming, using a 100-member ensemble of the Community Earth System Model version 2 Large Ensemble forced by two different BBA emissions (CMIP6 default and temporally smoothed over 1990–2020). Differences between the two subensemble means show that ocean temperature anomalies occur during periods of high BBA variability and subsequently persist over multiple decades. In the North Atlantic, surface warming is efficiently compensated for by decreased northward oceanic heat transport due to a slowdown of the Atlantic Meridional Overturning Circulation. In the North Pacific, surface warming is compensated for by an anomalous cross-equatorial cell (CEC) that reduces northward oceanic heat transport. The heat that converges in the South Pacific through the anomalous CEC is shunted into the subsurface and contributes to formation of long-lasting ocean temperature anomalies. The anomalous CEC is maintained through latitude-dependent contributions from narrow western boundary currents and basin-wide near-surface Ekman transport. Furthermore, these results indicate that interannual variability in forcing fields may significantly change the background climate state over long timescales, presenting a potential uncertainty in CMIP6-class climate projections forced without interannual variability.

54 ENVIRONMENTAL SCIENCES↗

The effects of bias, drift, and trends in calculating anomalies for evaluating skill of seasonal-to-decadal initialized climate predictions

In initialized seasonal to decadal (S2D) predictions, model hindcasts rapidly drift away from the initial observed state and converge toward a preferred state characterized by systematic error, or bias. Bias and drift are among the greatest challenges facing initialized prediction today. Differences in trends between initial states and drifted states, combined with bias and drift, introduce complexities in calculating anomalies to assess skill of initialized predictions. We examine several methods of calculating anomalies using the Decadal Prediction Large Ensemble (DPLE) using the Community Earth System Model (CESM) initialized hindcasts and focus on Pacific and Atlantic SSTs to illustrate issues with anomaly calculations. Three methods of computing anomalies, one as differences from a long term model climatology, another as bias-adjusted differences from the previous 15 year average from observations, and a third as differences from the previous 15 year average from the model, are contrasted and each is shown to have limitations. For the first, trends in bias and drift introduce higher skill estimates earlier and later in the hindcast period due to the trends that contribute to skill. For the second, higher skill can be introduced in situations where low frequency variability in the observations is large compared to the hindcasts on timescales greater than 15 years, while lower skill can result if the predicted signal is small and the bias-correction itself produces a transition of SST anomalies to the opposite sign of those that are observed. The third method has somewhat lower skill compared to each of the others, but has less difficulties with not only the long term trends in the model climatology, but also with the unrealistic situational skill from using observations as a reference. However, the first 15 years of the hindcast period cannot be evaluated due to having to wait to accumulate the previous 15 year model climatology before the method can be applied. Here, the IPO transition in the 2014–2016 time frame from negative to positive (predicted by Meehl et al. in in Nat Commun, 10.1038/NCOMMS11718, 2016) did indeed verify using all three methods, though each provides somewhat different skill values as a result of the respective limitations. There is no clear best method, as all are roughly comparable, and each has its own set of limitations and caveats. However, all three methods show generally higher overall skill in the AMO region compared to the IPO region.

54 ENVIRONMENTAL SCIENCES↗

Aerosol-Cloud Interactions and the Intersection of Climate Forcing and Feedback

The goal of this project has been to understand at the process level how the interaction of aerosols and cloud microphysics may alter the cloud response to climate change from the scale of turbulence to the global scale. Atmospheric aerosols (suspended particulates) have a unique role to play in the climate system. Anthropogenic aerosols may be an air pollution hazard and degrade human health. Their presence also results in the direct scattering or absorption of sunlight, with small net impacts on climate. But aerosols are also a key part of cloud formation, as the locations on which cloud droplets and ice crystals form. Understanding these Aerosol Cloud Interactions (ACI) is a critical uncertainty for climate models, including the Energy Exascale Earth System Model (E3SM), and affects how well models like E3SM can simulate the past and future. This work aims to understand cloud microphysics as well as the interaction of aerosol forcing and cloud feedback through improving our understanding of aerosol-cloud interactions (ACI) by looking at their uncertainty in critical environments at the process level. It will also contribute to the development of E3SM. Specific goals include: first, high level analysis of Perturbed Parameter Ensemble (PPE) simulations already available from the Community Earth System Model (CESM) and E3SM. Second, detailed analysis of recent changes to temperature and the role of aerosols.

54 ENVIRONMENTAL SCIENCES↗

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗

Dry Deposition Methods Based on Turbulence Kinetic Energy: 1. Evaluation of Various Resistances and Sensitivity Studies Using a Single‐Point Model

Different functions are used to account for turbulence strength in the atmospheric boundary layer for different stability regimes. These functions are one of the sources for differences among different atmospheric models' predictions and associated biases. Also, turbulence strength is underrepresented in some of the resistance formulations. Here, to address these issues with dry deposition, firstly we take advantage of three‐dimensional (3‐D) turbulence information in estimating resistances by proposing and validating a 3‐D turbulence velocity scale that is relevant for different stability regimes of boundary layer. Secondly, we hypothesize and validate that friction velocity measured by 3‐D sonic anemometer can be effectively replaced by the new turbulence velocity scale multiplied by the von Karman constant. Finally, we (a) present a set of resistance formulations for ozone (O 3 ) based on the 3‐D turbulence velocity scale; (b) intercompare estimations of such resistances with those obtained using existing formulations; and, (c) evaluate simulated O 3 fluxes using a single‐point dry deposition model against long‐term observations of O 3 fluxes at the Harvard Forest (MA) site. Results indicate that the new resistance formulations work very well in simulating surface latent heat and O 3 fluxes when compared to respective existing formulations and measurements at a decadal time scale. Findings from this research may help to improve the capability of dry deposition schemes for better estimation of dry deposition fluxes and create opportunities for the development of a community dry deposition model for use in regional/global air quality models.

54 ENVIRONMENTAL SCIENCES↗

Climate change increases risk of extreme rainfall following wildfire in the western United States

Post-wildfire extreme rainfall events can have destructive impacts in the western United States. Using two climate model large ensembles, we assess the future risk of extreme fire weather events being followed by extreme rainfall in this region. By mid-21st century, in a high warming scenario (RCP8.5), we report large increases in the number of extreme fire weather events followed within 1 year by at least one extreme rainfall event. By 2100, the frequency of these compound events increases by 100% in California and 700% in the Pacific Northwest in the Community Earth System Model v1 Large Ensemble. We further project that more than 90% of extreme fire weather events in California, Colorado, and the Pacific Northwest will be followed by at least three spatially colocated extreme rainfall events within five years. Our results point to a future with substantially increased post-fire hydrologic risks across much of the western United States.

54 ENVIRONMENTAL SCIENCES↗

Biodiversity and Stoichiometric Plasticity Increase Pico‐Phytoplankton Contributions to Marine Net Primary Productivity and the Biological Pump

Earth System Models generally predict increasing upper ocean stratification from 21st century warming, which will cause a decrease in the vertical nutrient flux forcing declines in marine net primary productivity (NPP) and carbon export. Recent advances in quantifying marine ecosystem carbon to nutrient stoichiometry have identified large latitudinal and biome variability, with low-latitude oligotrophic systems harboring pico-sized phytoplankton exhibiting large phosphorus to carbon cellular plasticity. The climate forced changes in nutrient flux stoichiometry and phytoplankton community composition are thus likely to alter the ocean's biogeochemical response and feedback with the carbon-climate system. We have added three pico-phytoplankton functional types within the Biogeochemical Elemental Cycling component of the Community Earth System Model while incorporating variable cellular phosphorus to carbon stoichiometry for all represented phytoplankton types. The model simulates Prochlorococcus and Synechococcus populations that dominate the productivity and sinking carbon export of the tropical and subtropical ocean, and pico-eukaryote populations that contribute significantly to productivity and export within the subtropical to mid-latitude transition zone, with the western subtropical regions of each basin supporting the most P-poor stoichiometries. Subtropical gyre recirculation regions along the poleward flanks of surface western boundary currents are identified as regional hotspots of enhanced carbon export exhibiting C-rich/P-poor stoichiometry, preferentially inhabited by pico-eukaryotes and diatoms. Collectively, pico-phytoplankton contribute ∼58% of global NPP and ∼46% of global particulate organic carbon export below 100 m through direct and ecosystem processing pathways. Biodiversity and cellular nutrient plasticity in marine pico-phytoplankton combine to increase their contributions to ocean productivity and the biological carbon pump.

Environmental Sciences & Ecology↗

From soil to sequence: filling the critical gap in genome-resolved metagenomics is essential to the future of soil microbial ecology

Abstract Soil microbiomes are heterogeneous, complex microbial communities. Metagenomic analysis is generating vast amounts of data, creating immense challenges in sequence assembly and analysis. Although advances in technology have resulted in the ability to easily collect large amounts of sequence data, soil samples containing thousands of unique taxa are often poorly characterized. These challenges reduce the usefulness of genome-resolved metagenomic (GRM) analysis seen in other fields of microbiology, such as the creation of high quality metagenomic assembled genomes and the adoption of genome scale modeling approaches. The absence of these resources restricts the scale of future research, limiting hypothesis generation and the predictive modeling of microbial communities. Creating publicly available databases of soil MAGs, similar to databases produced for other microbiomes, has the potential to transform scientific insights about soil microbiomes without requiring the computational resources and domain expertise for assembly and binning.

59 BASIC BIOLOGICAL SCIENCES↗

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗

Influences of North Pacific Ocean Domain Extent on the Western U.S. Winter Hydroclimatology in Variable-Resolution CESM

Variable-resolution global climate models (VRGCMs) are a dynamical downscaling method that can reach spatiotemporal scales needed for regional climate assessments. Over the years, several users of VRGCMs have assumed where the location and extent of the refinement domain should be based on knowledge of the prevailing storm tracks and resolution dependence of important regional climate processes (e.g., atmospheric rivers [ARs] and orographic uplift), but the effect of high-resolution domain size and extent on the simulation of downstream hydroclimatic phenomena has not been systematically evaluated. Here, we use variable resolution in the Community Earth System Model (VR-CESM) to perform such a test. To do this, three VR-CESM grids were generated that span the entire, two thirds, and one third of the North Pacific and evaluated for a 30-year climatology using Atmospheric Model Intercomparison Project protocols. Simulations are compared with reanalysis products offshore (fifth-generation of the European Centre for Medium-Range Weather Forecasts atmospheric reanalysis [ERA5]) and onshore (Livneh, 2015, https://doi.org/10.1038/sdata.2015.42, and Parameter-elevation Regressions on Independent Slopes Model [PRISM]) of the western United States. The westward expansion of refinement domain influenced integrated vapor transport (IVT), which was generally high biased but minimally impacted AR characteristics. Due to slight differences in landfalling AR counts in the western United States, California winter precipitation generally improved with westward expansion of the refinement domain. Western U.S. mountain snowpack and surface temperatures were insensitive to refinement domain size and were more influenced by changes in topographic resolution and/or land surface model version. Here, given minimal dependence of simulated western U.S. hydroclimate on refinement domain size over the North Pacific, we advise future VR-CESM studies to focus grid resolution on better resolving land surface heterogeneity.

54 ENVIRONMENTAL SCIENCES↗

A unified thermionic and thermionic-field emission (TE–TFE) model for ideal Schottky reverse-bias leakage current

We present a unified thermionic emission (TE) and thermionic-field emission (TFE) model for the ideal reverse-bias leakage current in Schottky junctions. The unified TE–TFE analytical model advances upon previous analytical TFE models by Murphy–Good and Padovani–Stratton, which are the two most widely adopted models by the community, in two major aspects: (i) the applicability of the TFE expression therein is extended to near-zero surface electric fields by an error-function correction, allowing for the calculation of the total current by a nontrivial sum of TE (over-the-barrier current) and TFE (below-the-barrier current) contributions; therefore, an accurate description of the TE-to-TFE transition region is captured analytically for the first time; (ii) image-force lowering is considered with much-simpler correction terms. Comparisons with the reference numerical model show that the unified TE–TFE model has excellent accuracy, as well as a 10 000× reduction in computation time. The unified model is further tested against experimental data from Schottky barrier diodes based on Si, 4H-SiC, GaN, and β-Ga 2 O 3 , revealing accurate extractions of barrier heights and correct descriptions of the ideal reverse leakage characteristics. With the extended applicable range, improved accuracy, and computational efficiency, the unified TE–TFE model is highly valuable for the design and analysis of devices based on Schottky junctions, as well as for potential integration in technology computer-aided design (TCAD) tools.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Decadal Thermal Variability of the Upper Southern Ocean: Zonal Asymmetry

Abstract As the major sink of anthropogenic heat, the Southern Ocean has shown quasi-symmetric, deep-reaching warming since the mid-twentieth century. In comparison, the shorter-term heat storage pattern of the Southern Ocean is more complex and has notable impacts on regional climate and marine ecosystems. By analyzing observational datasets and climate model simulations, this study reveals that the Southern Ocean exhibits prominent decadal (>8 years) variability extending to ∼700-m depth and is characterized by out-of-phase changes in the Pacific and Atlantic–Indian Ocean sectors. Changes in the Pacific sector are larger in magnitude than those in the Atlantic–Indian Ocean sectors and dominate the total heat storage of the Southern Ocean on decadal time scales. Instead of heat uptake through surface heat fluxes, these asymmetric variations arise primarily from wind-driven heat redistribution. Pacemaker and preindustrial simulations of the Community Earth System Model version 1 (CESM1) suggest that these variations in Southern Ocean winds arise primarily from natural variability of the tropical Pacific, as represented by the interdecadal Pacific oscillation (IPO). Through atmospheric teleconnection, the positive phase of the IPO gives rise to higher-than-normal sea level pressure and anticyclonic wind anomalies in the 50°–70°S band of the Pacific sector. These winds lead to warming of 0–700 m by driving the convergence of warm water. The opposite processes, involving cyclonic winds and upper-layer divergence, occur in the Atlantic–Indian Ocean sector. These findings aid our understanding of the time-varying heat storage of the Southern Ocean and provide useful implications on initialized decadal climate prediction.

Meteorology & Atmospheric Sciences↗

Ensemble Spread Behavior in Coupled Climate Models: Insights From the Energy Exascale Earth System Model Version 1 Large Ensemble

AbstractAssessing uncertainty in future climate projections requires understanding both internal climate variability and external forcing. For this reason, single‐model initial condition large ensembles (SMILEs) run with Earth System Models (ESMs) have recently become popular. Here we present a new 20‐member SMILE with the Energy Exascale Earth System Model version 1 (E3SMv1‐LE), which uses a “macro” initialization strategy choosing coupled atmosphere/ocean states based on inter‐basin contrasts in ocean heat content (OHC). The E3SMv1‐LE simulates tropical climate variability well, albeit with a muted warming trend over the twentieth century due to overly strong aerosol forcing. The E3SMv1‐LE's initial climate spread is comparable to other (larger) SMILEs, suggesting that maximizing inter‐basin ocean heat contrasts may be an efficient method of generating ensemble spread. We also compare different ensemble spread across multiple SMILEs, using surface air temperature and OHC. The Community Earth system Model version 1, the only ensemble which utilizes a “micro” initialization approach perturbing only atmospheric initial conditions, yields lower spread in the first ∼30 years. The E3SMv1‐LE exhibits a relatively large spread, with some evidence for anthropogenic forcing influencing spread in the late twentieth century. However, systematic effects of differing “macro” initialization strategies are difficult to detect, possibly resulting from differing model physics or responses to external forcing. Notably, the method of standardizing results affects ensemble spread: control simulations for most models have either large background trends or multi‐centennial variability in OHC. This spurious disequlibrium behavior is a substantial roadblock to understanding both internal climate variability and its response to forcing.

Stevenson, Samantha↗

Diagnosing and Predicting Mixed-Culture Fermentations with Unicellular and Guild-Based Metabolic Models

Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

59 BASIC BIOLOGICAL SCIENCES↗