SEARCH · Engineering Papers

Results for “Energy Exascale Earth System Model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Disentangling the chemistry and transport impacts of the quasi-biennial oscillation on stratospheric ozone

The quasi-biennial oscillation (QBO) in tropical winds perturbs stratospheric ozone throughout much of the atmosphere via changes in transport of ozone and other trace gases, as well as via temperature changes, both of which alter ozone chemistry. Attributing these causes of QBO–ozone variability may provide insights into model-to-model differences that contribute to ozone simulation. Here we develop a novel metric of steady-state ozone (SSO) to separate these effects: SSO calculates the local steady-state response of ozone due to the changes in temperature, chemical species, and overhead ozone column; the response due to circulation change is presumed when SSO shows no response. It is applied to the nudged Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2) with interactive ozone chemistry to demonstrate its validity. The E3SMv2 simulations nudged to reanalysis data produced reasonable wind and ozone patterns, especially in the tropics. Consistent with previous studies, we find clear demarcations with pressure. Ozone perturbations in the upper stratosphere (<6 hPa) are predicted by temperature changes; those between 6 and 20 hPa are predicted by NO y changes, and those in the lower stratosphere show no temperature or NO y response and are presumably driven by circulation changes. These results are important for diagnosing model-to-model discrepancy in QBO–ozone response and enhancing the reliability of ozone projections.

Xie, Jinbo [Lawrence Livermore National Laboratory↗

DOE OSTI2024

Droplet collection efficiencies inferred from satellite retrievals constrain effective radiative forcing of aerosol–cloud interactions

Abstract. Process-oriented observational constraints for the anthropogenic effective radiative forcing due to aerosol–cloud interactions (ERFaci) are highly desirable because the uncertainty associated with ERFaci poses a significant challenge to climate prediction. The contoured frequency by optical depth diagram (CFODD) analysis supports the evaluation of model representation of cloud liquid-to-rain conversion processes because the slope of a CFODD, generated from joint MODerate Resolution Imaging Spectroradiometer (MODIS)-CloudSat cloud retrievals, provides an estimate of cloud droplet collection efficiency in single-layer warm liquid clouds. Here, we present an updated CFODD analysis as an observational constraint on the ERFaci due to warm rain processes and apply it to the U.S. Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). A series of sensitivity experiments shows that E3SMv2 droplet collection efficiencies and ERFaci are highly sensitive to autoconversion, i.e., the rate of mass transfer from cloud liquid to rain, yielding a strong correlation between the CFODD slope and the shortwave component of ERFaci (ERFaciSW; Pearson's R=-0.91). E3SMv2's CFODD slope (0.20 ± 0.04) is in agreement with observations (0.20 ± 0.03). The strong sensitivity of ERFaciSW to the CFODD slope provides a useful constraint on highly uncertain warm rain processes, whereby ERFaciSW, constrained by MODIS-CloudSat, is estimated by calculating the intercept of the linear association between the ERFaciSW and the CFODD slopes, using the MODIS-CloudSat CFODD slope as a reference.

Beall, Charlotte M. (ORCID:0000000311370835)↗

DOE OSTI2025

Evaluating ecosystem water use efficiency and recovery dynamics during flash droughts: insights from observations and model simulations

Flash droughts (FD), rapidly emerging in a warming future, disrupt ecosystems, agriculture, and water security. Ecosystem water use efficiency (WUE), the ratio of gross primary production (GPP) to actual evapotranspiration (AET), balances carbon assimilation and water loss. FD rapidly disrupts this balance, making WUE critical for assessing plant stress and recovery. Here, this study investigates the dynamics of landscape-scale WUE, and the components of GPP and AET under FD utilizing both observed data from the Missouri Ozark AmeriFlux site (US-MOz) and version 2 of the U.S. Department of Energy’s Earth, Energy, Exascale System Model (E3SM) Land Model (ELMv2). Observations and simulations reveal GPP as dominant for WUE during earlier FD events (2005, 2007, 2012), shifting to AET in recent events (2014, 2018). This agreement indicates that the ELM can capture the shifting dynamics of GPP and AET in regulating WUE under FD conditions. However, the ELM systematically underestimates both GPP and AET and does so in a manner that does not preserve their ratio. As a result, WUE is also underestimated, suggesting that GPP is more strongly underestimated than AET. Furthermore, the ELM also underestimates the speed of GPP recovery, producing an artificially prolonged GPP recovery time following FD events. Observed environmental drivers such as vapor pressure deficit (VPD), soil moisture (SM), and predawn leaf water potential (PLWP) effectively predict WUE, but ELM primarily highlights SM, underestimating VPD’s role. This study demonstrates that relying solely on soil moisture fails to capture the rapid hydraulic recovery observed in PLWP, underscoring the necessity of integrating plant hydraulics into land surface models to improve flash drought predictability.

Evapotranspiration↗

DOE OSTI2026

ELM‐MOSART‐DOC: A Large‐Scale Riverine Dissolved Organic Carbon Model and Its Application Over the United States

Riverine dissolved organic carbon (DOC), primarily sourced from soil organic carbon (SOC), plays a crucial role in regional and global carbon cycles. However, the complexities of the underlying mechanisms and limited observations present significant challenges for predictive understanding of DOC at regional or larger scales. Recently, we developed a machine learning‐based (ML) map of DOC transformation rates, bridging the gap between SOC and DOC leaching flux and simplifying terrestrial DOC representation. Building on this advancement, we introduce ELM‐MOSART‐DOC, a DOC module integrated into the riverine component of the Energy Exascale Earth System Model (E3SM)—the Model for Scale Adaptive River Transport (MOSART). ELM‐MOSART‐DOC simulates DOC transport and transformation across both headwater streams and river networks, including those managed. Model validation demonstrates the ability of ELM‐MOSART‐DOC to accurately capture long‐term average DOC concentrations, with Kling‐Gupta Efficiency (KGE) scores of 0.58 and 0.76 at large and local stations, respectively. We further assess the impact of reservoirs through different simulation schemes, revealing that reservoirs significantly alter DOC fluxes by regulating streamflow patterns and promoting DOC mineralization. Model simulations indicate that reservoirs reduce total DOC flux from the Mississippi River into the ocean by 7.5%, with the long‐term average annual export decreasing from 3.34 to 3.14 teragrams (Tg) per year. ELM‐MOSART‐DOC integrates process‐based modeling with ML parameterization to enhance the predictive understanding of riverine biogeochemical processes. This approach reduces uncertainties in modeling regional and global carbon cycle ESMs and provides new insights into carbon cycling and its implications for global environmental change.

Li, Lingbo [Univ. of Houston, TX (United States); ↗

DOE OSTI2023

Peatland hydrology across scale: a probabilistic framework for confronting variability, heterogeneity, and uncertainty (Final Report)

Our project aimed to guide the representations of peatland hydrology and water-carbon feedbacks within Earth System Models, by improving process understanding and developing parsimonious models that elucidate the effects of hydroclimatic variability and spatial heterogeneity. Our research questions centered on the three key drivers of peatland hydrology: (1) seasonal and interannual hydroclimatic fluctuations, (2) spatial heterogeneity of peatland microtopography, and (3) hydrological connectivity across landscape units within a peatland watershed. Our analyses and model testing used existing long-term datasets within the Marcell Experimental Forest (MEF), with new data collected whenever necessary. To date, we have demonstrated the strong influence of water table elevations on the temperature sensitivity of CH 4 emissions, using a newly developed eddy covariance dataset spanning eleven years at MEF. Specifically, higher water tables dampen the increase in CH 4 emissions in the spring as well as their decrease in the fall, resulting in hysteresis. These results imply that any hydroclimatological changes in peatlands that shift seasonal water availability from winter to summer will increase annual CH 4 emissions, even if soil temperature remains unchanged. To further investigate water table elevation change within the peatland microtopography, we installed automated water table gauges across four bog-forest transects (the “lagg”) in two peatland watersheds at MEF, which are hotspots of intense biogeochemical activity. These measurements will give us new information about the extent of lagg expansion and contraction during snowmelt and high intensity rainfall events, as well as the directionality of water and nutrient flow into the bog from the surrounding forests. Finally, we investigated the hydrological connectivity across the peatland watershed complex during early spring – a critical period where water table is elevated from snowmelt – using extended hydrological records from MEF (e.g., snow and frost depths, water table elevations, streamflow). The importance of thermal and wintertime processes within the Energy Exascale Earth System Model (E3SM) was determined through sensitivity analyses. Results from data analyses show that (i) streamflow has decreased over decades likely due to increased evapotranspiration rates (despite no detectable trends in precipitation), (ii) spring streamflow generation is controlled by a sequence of fill-and-spill pathways from the snowpack to the stream through peatland water stores, and (iii) frost depth is a key explanatory variable for the timing and magnitude of streamflow. These results suggest that frost plays an important role in connecting surface water storage in peatlands to stream outlets, and that the hydrological connectivity across the peatland watershed complex can mediate the sensitivity of hydrological responses to climate variations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Ship-based Observations and Climate Model Simulations of Cloud Phase over the Southern Ocean

The Southern Ocean (SO) clouds exert a significant influence on the Earth’s radiation budget. Here we analyzed ship-based remote sensing observations of SO clouds over a five-month long DOE ARM Measurements of Aerosols, Radiation and Clouds over the Southern Oceans (MARCUS) field campaign to better understand cloud phase variability. We developed a method to classify eight categories of hydrometeors (ice, liquid, mixed phase, rain, drizzle, snow, aerosols, and clear sky) based on measurements of lidar, radar and radiosondes. Cloud thermodynamic phases (liquid, ice and mixed phase) at coarser scales were further derived to compare with the DOE Energy Exascale Earth System Model version 1 (E3SMv1) simulation and the Earth Model Column Collaboratory (EMC 2 ) instrument simulator. For a scale-aware comparison with climate model simulations, we found that spatially averaging the raw remote sensing data (e.g., backscatter, reflectivity) results in increased cloud cover and cloud liquid, whereas, directly averaging cloud phase from higher to lower resolution maintains clear air regions and is thus recommended for in-cloud frequency comparisons. For cloud thermodynamic phases in stratiform clouds, the E3SM underestimates cloud ice and overestimates cloud liquid at temperatures between -40 and 0°C. When latitudes increase, both observations and simulations show a transition of dominant phase from liquid to ice for cloud tops as well as for the entire cloud columns, but the model underestimation of ice phase is more severe at higher latitudes. Such model bias is unlikely caused by spatial scale differences or lack of heterogeneity in cloud vertical structure in the simulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Shrub Expansion Simulations at Trail Valley Creek Tundra site using E3SM Land Model (ELM) Arctic-focused Version

The warming of the Arctic is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and densification in parts of the Arctic. Assessing the impact of these changes in vegetation composition on the Arctic’s carbon and energy budgets is important to constrain projected local and global surface-atmosphere exchanges. We conduct a sensitivity analysis of the projected surface energy fluxes, soil carbon pools, and carbon dioxide fluxes (net ecosystem exchange, gross primary production, and ecosystem respiration) between present day and 2100 to different shrub expansion rates and air temperature increases under future emission scenarios (intermediate – RCP4.5, and high – RCP8.5) using the Arctic-focused version of the Energy Exascale Earth System Model (E3SM) Land Model (ELM). We focus on Trail Valley Creek (TVC), a mineral upland tundra site located in the western Canadian Arctic, which is experiencing tall shrub densification and expansion. In this study, we run TVC under two different warming scenarios RCP4.5 and RCP 8.5 and simulate different shrubification rates projected until year 2100. In this repository, we include all the forcing, input, parameters, and output data corresponding to all the simulations performed. flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

In a computational fluid model of the atmosphere, the advective transport of trace species, or tracers, can be computationally expensive. For efficiency, models often use semi-Lagrangian advection methods. High-order interpolation semi-Lagrangian (ISL) methods, in particular, can be extremely efficient, if the problem of property preservation specific to them can be addressed. Atmosphere models often use geometrically and logically nonuniform grids for efficiency and, as a result, element-based discretizations. Such grids and discretizations make stability a particular problem for ISL methods. Generally, high-order, element-based ISL methods that use the natural polynomial interpolant associated with a nodal finite-element discretization are unstable. Here, we derive new bases having order of accuracy up to nine, with positive nodal weights, that stabilize the element-based ISL method. We use these bases to construct the linear advection operator in the property-preserving Interpolation Semi-Lagrangian Element-based Transport (Islet) method. Then we discuss key software implementation details. Finally, we show performance results for the Energy Exascale Earth System Model's atmosphere dynamical core, comparing the original and new transport methods. These simulations used up to 27,600 Graphical Processing Units (GPU) on the Oak Ridge Leadership Computing Facility's Summit supercomputer.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗

NASA NTRS2023

Changing Characteristics of Tropical Extreme Precipitation-Cloud Regimes in Warmer Climates

In this paper, we have investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). Results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by a) control experiment with present-day sea surface temperature (SST) and CO 2 concentration, b) P4K, same as in a) but with uniform increase by 4K in SST globally, and c) same as in a), but with imposed SST and CO 2 concentration from outputs of coupled E3SM forced by 4xCO 2 concentration. We find that as the surface warms under P4K and 4xCO 2 , both convective and stratiform rain increase. Importantly, there is increasing fractional contribution of stratiform rain as a function of precipitation intensity, with the most extreme but rare events occurring preferentially over land more than ocean, more so under 4xCO 2 than P4K. Extreme precipitation is facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation-cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in vertical profiles of clouds, condensation heating and vertical motions indicate increasing precipitation-cloud-circulation organization from control, P4K to 4xCO 2 . Results suggest that large-scale ocean warming, i.e., P4K, is the primary cause contributing to organization structure resembling the well-known Mesoscale Convective System (MCS) for increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating, and dynamically consistent anomalous SST further amplifies the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increase in surface moisture (temperature) under P4K and 4xCO 2 is the key driver leading to enhanced convective instability over tropical ocean (land). However, the fast and large increase in land surface temperature and lack of available local moisture result in strong reduction in land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, more so under 4xCO 2 , are likely due to the delayed onset of deep convection, i.e., the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when eventually occur, due to the release of large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

climate-scale extreme tropical precipitation↗

DOE OSTI2023

Changing Characteristics of Tropical Extreme Precipitation–Cloud Regimes in Warmer Climates

In this study, we investigated the changing characteristics of climatic scale (monthly) tropical extreme precipitation in warming climates using the Energy Exascale Earth System Model (E3SM). The results are from Atmospheric Model Intercomparison Project (AMIP)-type simulations driven by (a) a control experiment with the present-day sea surface temperature (SST) and CO 2 concentration, (b) P4K, the same as in (a) but with a uniform increase of 4K in the SST globally, and (c) the same as in (a), but with an imposed SST and CO 2 concentration from the outputs of the coupled E3SM forced by a 4xCO 2 concentration. We found that as the surface warmed under P4K and 4xCO 2 , both convective and stratiform rain increased. Importantly, there was an increasing fractional contribution of stratiform rain as a function of the precipitation intensity, with the most extreme but rare events occurring preferentially over land more than the ocean, and more so under 4xCO 2 than P4K. Extreme precipitation was facilitated by increased precipitation efficiency, reflecting accelerated rates of recycling of precipitation cloud water (both liquid and ice phases) in regions with colder anvil cloud tops. Changes in the vertical profiles of clouds, condensation heating, and vertical motions indicate increasing precipitation–cloud–circulation organization from the control and P4K to 4xCO 2 . The results suggest that large-scale ocean warming, that is, P4K, was the primary cause contributing to an organization structure resembling the well-known mesoscale convective system (MCS), with increased extreme precipitation on shorter (hourly to daily) time scales. Additional 4xCO 2 atmospheric radiative heating and dynamically consistent anomalous SST further amplified the MCS organization under P4K. Analyses of the surface moist static energy distribution show that increases in the surface moisture (temperature) under P4K and 4xCO 2 was the key driver leading to enhanced convective instability over tropical ocean (land). However, a fast and large increase in the land surface temperature and lack of available local moisture resulted in a strong reduction in the land surface relative humidity, reflecting severe drying and enhanced convective inhibition (CIN). It is argued that very extreme and rare “record-breaking” precipitation events found over land under P4K, and more so under 4xCO 2 , are likely due to the delayed onset of deep convection, that is, the longer the suppression of deep convection by CIN, the more severe the extreme precipitation when it eventually occurs, due to the release of a large amount of stored surplus convective available potential energy in the lower troposphere during prolonged CIN.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Quantifying the long-term changes of terrestrial water storage and their driving factors

Global warming is expected to cause changes in terrestrial water storage (TWS) across the land surface, with widespread impacts on ecosystems and society. Although extensive research has been performed to analyze TWS changes and possible drivers during the post-2000 period, longer-term evolution of TWS and associated environmental forcings remain relatively unexplored. In this study, we evaluated the performance of the Energy Exascale Earth System model (E3SM) land model ELM version 1 (ELM v1) in simulating global TWS, and used factorial simulations of ELMv1 to quantify global TWS changes and their drivers during 1948–2012. We found that ELM’s agreed best with existing satellites and reconstruction datasets in temperate regions unaffected by irrigation. Biome- and climate zone-averaged TWS mainly increased at rates between 0 and 10 mm/year over 1948–2012, but the second half of that period saw smaller positive trends than the first half or even negative trends. Climate change explained >80 % of the TWS trends across most biomes and climate zones, followed by land use and land cover change. The physiological and phenological effects of CO2 primarily induced noticeable TWS trends in the more humid biomes and climate zones across different latitudes. In contrast, nitrogen deposition and aerosol deposition generally had smaller and negative impacts across the biomes and climate regions. Among the meteorological drivers analyzed, the long-term average imbalance between precipitation (P), evapotranspiration (E), and runoff (Q) contributed >50 % of the TWS trends in most biomes and climate zones, with nonlinearity being induced by spatially heterogenous changes in E/P and Q/P ratios. The accumulated detrended anomalies in P, E, and Q also often contributed substantially, while the trends difference between P, E, and Q contributed little. Together, these findings unveiled an intensification of the global TWS and its diverse patterns of climate change and different non-withdrawal human-induced alterations, contributing to a more comprehensive understanding and projection of the global water cycle.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Storylines for the 1997 New Year’s Flood: The role of watershed antecedent conditions and future warming in shaping discharge in the Truckee River watershed

The 1997 New Year’s flood was among the most devastating floods in the Truckee River watershed located in western Nevada. This event resulted from complex interactions of flood drivers, such as extreme precipitation, wet antecedent watershed conditions, warm temperatures and rapid snowmelt. We leveraged simulated forcings from the regionally refined mesh capabilities of the Energy Exascale Earth System Model (RRM-E3SM) and a process-based hydrological model to recreate the 1997 New Year’s flood for the Truckee River watershed across four climate warming levels ranging from the current temperatures to + 4° C. For each scenario, we conducted ensemble simulations with the same forcing but with 100 different seasonal watershed antecedent conditions, which were randomly sampled from long-term hydrological simulations. The results show that the 1997 New Year’s flood can be reproduced or exceeded consistently only when the antecedent watershed conditions are wet, specifically when streamflows are above the 75th percentile of the climatological value. There is negligible change in ensemble mean peakflows for Truckee River near Reno; however, there are increases of 18% and 14% under the warming levels of + 3° C and + 4° C, respectively. The increases in peakflows under future climate warming are attributed to wetter antecedent watershed conditions and enhanced snowmelt. Furthermore, the largest increases in peakflows occur at small, high-elevation headwater basins along the Sierra Nevada crest. This study highlights that changes in extreme flood events will result from the complex interplay of multiple flood drivers. It also demonstrates the potential of storyline approaches to analyze future realizations of these extreme events under different climate scenarios.

Climate change impact study↗

DOE OSTI2021

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

The Oceanic Barrier Layer in the Eastern Indian Ocean as a Predictor for Rainfall Over Indonesia and Australia

Barrier layers in the tropics trap heat in a shallow and stable near-surface layer and limit entrainment of cooler water from below. Both processes act to increase sea surface temperature and enhance atmospheric convection. The high resolution fully coupled pre-industrial Energy Exascale Earth System Model version 0 (E3SMv0) is used to investigate the relationship between barrier layers in the eastern Indian Ocean during the wet season with local atmospheric convection and remote rainfall. A partial least squares regression reveals a significant relationship between Australasian rainfall and the barrier layer thickness (BLT) west of Sumatra, occurring one month earlier. Furthermore, the largest positive regression coefficients are over northern Australia. The region west of Sumatra is strategically located where the East-Asian monsoon moisture flows toward northern Australia. Thickening of the west Sumatra BLT intensifies evaporation and local convection and amplifies the moisture transported to Australia acting to increase the terrestrial rainfall.

58 GEOSCIENCES↗

DOE OSTI2021

Attribution of Snowpack Errors to Simulated Temperature and Precipitation in E3SMv1 Over the Contiguous United States

Abstract Snow water equivalent (SWE), temperature, and precipitation biases and trends are evaluated in the atmosphere‐land simulations of the Energy Exascale Earth System Model version 1 (E3SMv1) in comparison to the Community Earth System Model version 2 (CESM2) and two other models using the ground measurement‐based University of Arizona (UA) snow product. SWE, temperature, and precipitation biases are highest in magnitude in the Western contiguous United States (CONUS). SWE errors are attributed to temperature and precipitation through multiple linear regressions of normalized errors, the coefficients of which represent the sensitivities to temperature and precipitation errors. SWE errors are more sensitive to temperature errors throughout the CONUS. Model SWE and temperature trends are generally opposite from UA product trends in the Western CONUS. SWE trend errors are also attributed to temperature and precipitation trend errors using multiple linear regressions of normalized trend errors. SWE trend errors are more sensitive to those of precipitation at higher elevations (>1,500 m) in the Western CONUS in these simulations. Thus, the sensitivity to temperature and precipitation differ for SWE errors and its trend errors. Furthermore, the SWE trend errors are more sensitive to temperature and precipitation in the atmosphere‐ocean coupled simulations in which the atmosphere‐land is coupled to an active ocean model. These results suggest that both errors in simulated temperature and precipitation contribute to SWE errors.

Meteorology & Atmospheric Sciences↗

DOE OSTI2022

Cloud Phase Simulation at High Latitudes in EAMv2: Evaluation using CALIPSO Observations and Comparison with EAMv1

This study performs a comprehensive evaluation of the simulated cloud phase in the U.S. Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) atmosphere model version 2 (EAMv2) and version 1 (EAMv1). Enabled by the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) simulator, EAMv2 and EAMv1 predicted cloud phase is compared against the GCM-Oriented CALIPSO Cloud Product (CALIPSO-GOCCP) at high latitudes where mixed-phase clouds are prevalent. Our results indicate that the underestimation of cloud ice in simulated high-latitude mixed-phase clouds in EAMv1 has been significantly reduced in EAMv2. The increased ice clouds in the Arctic mainly result from the modification on the WBF (Wegner-Bergeron-Findeisen) process in EAMv2. The impact of the modified WBF process is moderately compensated by the low limit of cloud droplet number concentration (CDNC) in cloud microphysics and the new dCAPE_ULL trigger used in deep convection in EAMv2. Moreover, it is found that the new trigger largely contributes to the better cloud phase simulation over the Norwegian Sea and Barents Sea in the Arctic and the Southern Ocean where large errors are found in EAMv1. However, errors in simulated cloud phase in EAMv1, such as the overestimation of supercooled liquid clouds near the surface in both hemispheres and the underestimation of ice clouds over Antarctica, persist in EAMv2. In conclusion, this study highlights the impact of deep convection parameterizations, which has received little attention, on high-latitude mixed-phase clouds, and the importance of continuous improvement of cloud microphysics in climate models for accurately representing mixed-phase clouds.

54 ENVIRONMENTAL SCIENCES↗