An Introduction to the E3SM Special Collection: Goals, Science Drivers, Development, and Analysis
Not Available
54 ENVIRONMENTAL SCIENCES↗
SEARCH · Engineering Papers
Results for “Earth Systems”
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.
Numerically Relevant Timescales in the MG2 Microphysics Model
Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.
58 GEOSCIENCES↗
Improving Time Step Convergence in an Atmosphere Model With Simplified Physics: The Impacts of Closure Assumption and Process Coupling
Convergence testing is a common practice in the development of dynamical cores of atmospheric models but is not as often exercised for the parameterization of sub-grid physics. An earlier study revealed that the stratiform cloud parameterizations in several predecessors of the Energy Exascale Earth System Model (E3SM) showed strong time-step sensitivity and slower-than-expected convergence when the model's time step was systematically refined. In this work, a simplified atmosphere model is configured that consists of the spectral-element dynamical core of the E3SM atmosphere model coupled with a large-scale condensation parameterization based on commonly used assumptions. This simplified model also resembles E3SM and its predecessors in the numerical implementation of process coupling and shows poor time-step convergence in short ensemble tests. We present a formal error analysis to reveal the expected time-step convergence rate and the conditions for obtaining such convergence. Numerical experiments are conducted to investigate the root causes of convergence problems. We show that revisions in the process coupling and closure assumption help to improve convergence in short simulations using the simplified model; the same revisions applied to a full atmosphere model lead to significant changes in the simulated long-term climate. This work demonstrates that causes of convergence issues in atmospheric simulations can be understood by combining analyses from physical and mathematical perspectives. Addressing convergence issues can help to obtain a discrete model that is more consistent with the intended representation of the physical phenomena.
54 ENVIRONMENTAL SCIENCES↗
Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM
21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.
Zhang, Meng↗
Different Spatiotemporal Patterns in Global Human Population and Built‐Up Land
Abstract Population concentration and built‐up land expansion are two prominent features of contemporary urbanization. Existing literature on the population aspect of urbanization has mostly focused on national and regional aggregates, and literature on the land development aspect has often relied on spatial case studies of individual cities or their meta‐analyses. Using newly‐available data, here we conduct the first global‐coverage, spatial analysis of the relationship between (changes in) population and built‐up land at multiple spatial scales, and compare to existing common beliefs about urbanization based on individual city studies. We find that population and built‐up land show distinctly different spatial and temporal patterns (with a global correlation coefficient around 0.6). Contrary to common impressions, our results show that during recent decades, developed and developing regions across the world experienced comparable amounts of built‐up land expansion. While meta‐analyses have reported that built‐up land in urban areas expands globally on average twice as fast as population grows, our results show the global change rates of built‐up land and population are similar. Also, most global population, including what national statistics agencies call urban population, reside in areas with low land development levels (which are frequently less than 5% built up). These changes in perspective suggest that urbanization's potential large‐scale impacts may need to be re‐evaluated, and lead to best‐practice recommendations for urbanization modeling and analysis. Especially, the common practice in large‐scale earth system modeling of assuming demographically‐defined urban population resides in areas with medium to high built‐up land development levels should change.
Gao, Jing↗
InterCarb: A Community Effort to Improve Interlaboratory Standardization of the Carbonate Clumped Isotope Thermometer Using Carbonate Standards
Abstract Increased use and improved methodology of carbonate clumped isotope thermometry has greatly enhanced our ability to interrogate a suite of Earth‐system processes. However, interlaboratory discrepancies in quantifying carbonate clumped isotope (Δ 47 ) measurements persist, and their specific sources remain unclear. To address interlaboratory differences, we first provide consensus values from the clumped isotope community for four carbonate standards relative to heated and equilibrated gases with 1,819 individual analyses from 10 laboratories. Then we analyzed the four carbonate standards along with three additional standards, spanning a broad range of δ 47 and Δ 47 values, for a total of 5,329 analyses on 25 individual mass spectrometers from 22 different laboratories. Treating three of the materials as known standards and the other four as unknowns, we find that the use of carbonate reference materials is a robust method for standardization that yields interlaboratory discrepancies entirely consistent with intralaboratory analytical uncertainties. Carbonate reference materials, along with measurement and data processing practices described herein, provide the carbonate clumped isotope community with a robust approach to achieve interlaboratory agreement as we continue to use and improve this powerful geochemical tool. We propose that carbonate clumped isotope data normalized to the carbonate reference materials described in this publication should be reported as Δ 47 (I‐CDES) values for Intercarb‐Carbon Dioxide Equilibrium Scale.
58 GEOSCIENCES↗
Where Are White Roofs More Effective in Cooling the Surface?
Solar reflective roofs (also called white or cool roofs) are a common heat mitigation strategy. By reflecting more sunlight than regular roofs, cool roofs experience lower surface temperatures and produce less heating of the surrounding air. The effectiveness of cool roofs, quantified by the surface temperature difference between cool roofs and regular roofs ??????, is known to vary spatially. A common perception is that ?????? is controlled by solar radiation that reaches the roof surface. Here we use an Earth System Model and a surface energy balance model to show that the spatial variability of ??????, when normalized by the albedo difference between cool roofs and regular roofs (????), is in fact more controlled by an energy distribution factor that encodes the efficiencies of surface energy balance components in dissipating heat. Our results suggest that painting the roof white is more effective in reducing the roof surface temperature when the roof has less water holding capacity and smaller thermal admittance and is located in places with more solar radiation, less precipitation, and lower wind speed.
Wang, Linying↗
The Shifting Scales of Western U.S. Landfalling Atmospheric Rivers Under Climate Change
Atmospheric rivers (ARs) can be a boon and bane to water resource managers as they have the ability to replenish water reserves, but they can also generate million-to-billion-dollar flood damages. To investigate how anthropogenic climate change may influence AR characteristics in the coastal western United States by end century, we employ a suite of novel tools such as variable resolution in the Community Earth System Model (VR-CESM), the TempestExtremes AR detection algorithm, and the Ralph, Rutz, et al. (2019, https://doi.org/10.1175/BAMS-D-18-0023.1) AR category scale.We show that end-century ARs primarily shift from being “mostly or primarily beneficial” to “mostly or primarily hazardous” with a concomitant sharpening and intensification of winter season precipitation totals. Changes in precipitation totals are due to a significant increase in AR (+260%) rather than non-AR (+7%) precipitation, largely through increases in the most intense category of AR events and a decrease in the interval between landfalling ARs.
54 ENVIRONMENTAL SCIENCES↗
The Role of the Mean State on MJO Simulation in CESM2 Ensemble Simulation
Here we examine the role of the mean state in the propagation of the Madden-Julian oscillation (MJO) over the Maritime Continent (MC). We use an ensemble of simulations made with a single model—the Community Earth System Model version 2—to assess the effect of the mean state that is unaffected by that of model components such as parameterization schemes. Results show that the background meridional moisture gradient is much steeper over the MC region in the periods with a stronger MJO propagation. Column water vapor budget of the MJO strongly suggests that the simulated mean state affects MJO via its impacts on moisture dynamics—a greater advection of mean moisture by MJO wind in the MC region is responsible for the anomalous MJO activity.
58 GEOSCIENCES↗
Seasonal Dependent Impact of Ice Cloud Longwave Scattering on the Polar Climate
Most climate models neglect cloud longwave (LW) scattering because scattering is considered negligible compared to strong LW absorption by clouds and greenhouse gases. While this rationale is valid for simulating extrapolar regions, it is questionable for the polar regions, where the atmosphere is dry and hence has weak absorption, and ice clouds that have strong scattering capability frequently occur. Using the slab-ocean Community Earth System Model, we show that ice cloud LW scattering can warm winter surface air temperature by 0.8–1.8 K in the Arctic and 1.3–1.9 K in the Antarctic, while this warming becomes much weaker in polar summer. Such scattering effect cannot be correctly assessed when sea surface temperature and sea ice are prescribed as this effect is manifested through a surface-atmosphere coupling. Furthermore, cloud LW scattering is a necessity for the correct simulation of polar climate and surface radiation budget, especially in the winter.
54 ENVIRONMENTAL SCIENCES↗
Parameterizing Subgrid Variations of Land Surface Heat Fluxes to the Atmosphere Improves Boreal Summer Land Precipitation Simulation With the NCAR CESM1.2
Subgrid horizontal variations of land surface heat fluxes to the atmosphere resulting from subgrid land cover heterogeneity are important in land-atmosphere interaction in global climate models (GCMs). To incorporate it, a parameterization using stochastic sampling based on truncated normal distributions diagnosed from the land model and internal ensemble mean of multiple calls to the planetary boundary layer and deep convection schemes is developed. After its implementation in the National Center for Atmospheric Research (NCAR) Community Earth System Model version 1.2 (CESM1.2), large changes of simulated land precipitation are found to occur in regions where large subgrid variations of surface heat fluxes exist. The simulated precipitation is improved in boreal summer, especially over eastern China and the coastal areas of the Bay of Bengal. The improved precipitation is mainly a result of improved large-scale moisture convergence and advection by altered vertical diffusion and convection.
58 GEOSCIENCES↗
The Climate Response to Emissions Reductions Due to COVID-19: Initial Results From CovidMIP
Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.
54 ENVIRONMENTAL SCIENCES↗
A Lagrangian perspective on tropical anvil cloud lifecycle in present and future climate
The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.
climate, weather, mesocale cellular convection↗
Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5
Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.
58 GEOSCIENCES↗
Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales
Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.
54 ENVIRONMENTAL SCIENCES↗
Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model
Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.
54 ENVIRONMENTAL SCIENCES↗
Cloud Process Coupling and Time Integration in the E3SM Atmosphere Model
Abstract In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 s. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step and, in particular, on the frequency with which cloud‐related processes are coupled with each other and with the model dynamics.
54 ENVIRONMENTAL SCIENCES↗
Integrating Arctic Plant Functional Types in a Land Surface Model Using Above‐ and Belowground Field Observations
Abstract Accurate simulations of high‐latitude ecosystems are critical for confident Earth system model (ESM) projections of carbon cycle feedbacks to global climate change. Land surface model components of ESMs, including the E3SM Land Model (ELM), simulate vegetation growth and ecosystem responses to changing climate and atmospheric CO 2 concentrations by grouping heterogeneous vegetation into like sets of plant functional types (PFTs). Many such models represent high‐latitude vegetation using only two PFTs (shrub and grass), thereby missing the diversity of vegetation growth forms and functional traits in the Arctic. Here, we use field observations of biomass and leaf traits across a gradient of plant communities on the Seward Peninsula in northwest Alaska to replace the original ELM configuration for the first time with nine Arctic‐specific PFTs. The newly developed PFTs include: (1) nonvascular mosses and lichens, (2) deciduous and evergreen shrubs of various height classes, including an alder PFT, (3) graminoids, and (4) forbs. Improvements relative to the original model configuration included greater belowground biomass allocation, persistent fine roots and rhizomes of nonwoody plants, and better representation of variability in total plant biomass across sites with varying plant communities and depth to bedrock. Simulations through 2100 using the RCP8.5 climate scenario and constant PFT fractional areas showed alder‐dominated plant communities gaining more biomass and lichen‐dominated communities gaining less biomass compared to default PFTs. Our results highlight how representing the diversity of arctic vegetation and confronting models with measurements from varied plant communities improves the representation of arctic vegetation in terrestrial ecosystem models.
54 ENVIRONMENTAL SCIENCES↗