SEARCH · Engineering Papers
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.
We demonstrate current capabilities of the NASA finite-volume General Circulation Model an high-resolution global weather prediction, and discuss its development path in the foreseeable future. This model can be regarded as a prototype of a future NASA Earth modeling system intended to unify development activities cutting across various disciplines within the NASA Earth Science Enterprise.
Permafrost currently stores more than a fourth of global soil carbon. A warming climate makes this carbon increasingly vulnerable to decomposition and release into the atmosphere in the form of greenhouse gases. The resulting climate feedback can be estimated using land surface models, but the high complexity and computational cost of these models make it challenging to use them for estimating uncertainty, exploring novel scenarios, and coupling with other models. We have added a representation of permafrost to the simple, open-source global carbon–climate model Hector, calibrated to be consistent with both historical data and 21st century Earth system model projections of permafrost thaw. We include permafrost as a separate land carbon pool that becomes available for decomposition into both methane (CH4) and carbon dioxide (CO2) once thawed; the thaw rate is controlled by region-specific air temperature increases from a preindustrial baseline. We found that by 2100 thawed permafrost carbon emissions increased Hector’s atmospheric CO2 concentration by 5 %–7 % and the atmospheric CH4 concentration by 7 %–12 %, depending on the future scenario, resulting in 0.2–0.25 ∘C of additional warming over the 21st century. The fraction of thawed permafrost carbon available for decomposition was the most significant parameter controlling the end-of-century temperature change in the model, explaining around 70 % of the temperature variance, and was distantly followed by the initial stock of permafrost carbon, which contributed to about 10 % of the temperature variance. The addition of permafrost in Hector provides a basis for the exploration of a suite of science questions, as Hector can be cheaply run over a wide range of parameter values to explore uncertainty and can be easily coupled with integrated assessment and other human system models to explore the economic consequences of warming from this feedback.
Historical emission estimates for anthropogenic aerosol and precursor compounds are key data needed for Earth system models, climate models, and atmospheric chemistry and transport models; both for general analysis and assessment and also for model validation through comparisons with observations. Current global emission data sets have a number of shortcomings, including timeliness and transparency. Satellite and other earth-system data are increasingly available in near real-time, but global emission estimates lag by 5-10 years. The CEDS project will construct a data-driven, open source framework to produce annually updated emission estimates. The basic methodologies to be used for this system have been used for SO2 (Smith et al. 2011, Klimont, Smith and Cofala 2013), and are designed to complement existing inventory efforts. The goal of this system is to consistently extend current emission estimates both forward in time to recent years and also back over the entire industrial era. The project will produce improved datasets for global and (potentially) regional model, allow analysis of trends across time, countries, and sectors of emissions and emission factors, and facilitate improved scientific analysis in general. Consistent estimation of uncertainty will be an integral part of this system. This effort will facilitate community evaluation of emissions and further emission-related research more generally.
Historical emission estimates for anthropogenic aerosol and precursor compounds are key data needed for Earth system models, climate models, and atmospheric chemistry and transport models; both for general analysis and assessment and also for model validation through comparisons with observations. Current global emission data sets have a number of shortcomings, including timeliness and transparency. Satellite and other earth-system data are increasingly available in near real-time, but global emission estimates lag by 5-10 years. The CEDS project will construct a data-driven, open source framework to produce annually updated emission estimates. The basic methodologies to be used for this system have been used for SO2 (Smith et al. 2011, Klimont, Smith and Cofala 2013), and are designed to complement existing inventory efforts. The goal of this system is to consistently extend current emission estimates both forward in time to recent years and also back over the entire industrial era. The project will produce improved datasets for global and (potentially) regional model, allow analysis of trends across time, countries, and sectors of emissions and emission factors, and facilitate improved scientific analysis in general. Consistent estimation of uncertainty will be an integral part of this system. This effort will facilitate community evaluation of emissions and further emission-related research more generally.
Abstract We here present, document and validate a new atmospheric component of the Community Earth System Model (CESM1): the Specified Chemistry Whole Atmosphere Community Climate Model (SC‐WACCM). As the name implies, SC‐WACCM is a middle atmosphere‐resolving model with prescribed, rather than interactive chemistry. Ozone concentrations are specified throughout the atmosphere, using zonal and monthly mean climatologies computed by a companion integration with the Whole Atmosphere Community Climate Model (WACCM). Above 65 km, in addition, the climatological chemical and shortwave heating, nitrogen oxide, atomic and molecular oxygen and carbon dioxide are also prescribed from the companion WACCM integration. We carefully compare the climatology and the climate variability of preindustrial integrations of SC‐WACCM and WACCM each coupled with active land, ocean and sea‐ice components. We note some differences in upper stratospheric and lower mesospheric temperature, just below the 65 km transition level, due to the diurnal ozone cycle that is not captured when monthly mean ozone is used. Nonetheless, we find that the climatology and variability of the stratosphere, the troposphere and surface climate are nearly identical in SC‐WACCM and WACCM. Notably, the frequency and amplitude of Northern Hemisphere stratospheric sudden warmings in the two integrations are not significantly different. Also, we compare WACCM and SC‐WACCM to CCSM4, the “low‐top” version of CESM1, and we find very significant differences in the stratospheric climatology and variability. The removal of the chemistry reduces the computational cost of SC‐WACCM to approximately one half of WACCM: in fact, SC‐WACCM is only 2.5 times more expensive than CCSM4 at the same horizontal resolution. This considerable reduction in computational cost makes the new SC‐WACCM component of CESM1 ideally suited for studies of stratosphere‐troposphere dynamical coupling and, more generally, the role of the stratosphere in the climate system.
In this paper, the recently developed variable-resolution option within the Community Earth System Model (VR-CESM) is assessed for long-term regional climate modeling of California at 0.25° (~ 28 km) and 0.125° (~ 14 km) horizontal resolutions. The mean climatology of near-surface temperature and precipitation is analyzed and contrasted with reanalysis, gridded observational data sets, and a traditional regional climate model (RCM)—the Weather Research and Forecasting (WRF) model. Statistical metrics for model evaluation and tests for differential significance have been extensively applied. VR-CESM tended to produce a warmer summer (by about 1–3°C) and overestimated overall winter precipitation (about 25%–35%) compared to reference data sets when sea surface temperatures were prescribed. Increasing resolution from 0.25° to 0.125° did not produce a statistically significant improvement in the model results. By comparison, the analogous WRF climatology (constrained laterally and at the sea surface by ERA-Interim reanalysis) was ~1–3°C colder than the reference data sets, underestimated precipitation by ~20%–30% at 27 km resolution, and overestimated precipitation by ~ 65–85% at 9 km. Overall, VR-CESM produced comparable statistical biases to WRF in key climatological quantities. Moreover, this assessment highlights the value of variable-resolution global climate models (VRGCMs) in capturing fine-scale atmospheric processes, projecting future regional climate, and addressing the computational expense of uniform-resolution global climate models.
This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.
Explicit representation of soil microbial processes and interactions with biotic and abiotic processes in Earth System Models (ESMs) remains limited, despite their importance in biogeochemical cycles. To address this gap, which hinders prediction of global biogeochemial cycling and responses to atmospheric conditions, we integrated a microbe- and mineral-surface-explicit model, the Reaction-network-based model of soil organic matter and Microbes (ReSOM), into the Energy Exascale ESM (E3SM) land model (ELM). Here, we describe ELM-ReSOM and show a case study at a conifer forest in California. ELM-ReSOM accurately simulated surface CO 2 fluxes and SOM stocks, demonstrating improved representations of microbial and mineral interactions compared to the default ELM. We examined ELM-ReSOM sensitivity to microbial traits, enzyme properties, and organo-mineral interactions. Microbial traits such as the maximum mortality rate, transporter-density scaling factor, and maximum monomer assimilation rate were strong controllers of heterotrophic respiration, while these microbial traits and enzyme-related properties collectively influenced SOM stocks. Mineral surfaces primarily affected SOM stocks by adsorbing enzymes, thereby limiting depolymerization. Synergies among processes led to stronger impacts of parameters when evaluated together versus separately (i.e., most parameters had greater indirect than direct effects). For example, due to interactions of microbial necromass with mineral surface adsorption, the indirect effect of the maximum microbial mortality rate was 33% larger than its direct effect on SOM stock. Thus, microbial and enzyme dynamics and their interactions with mineral surfaces play critical roles in SOM cycling. Tackling the challenges of microbe-explicit models will advance understanding and modeling of SOM dynamics.
The transport of water, heat, and momentum from the surface to the atmosphere is dependent in part on the 10 characteristics of the land surface. Together with the model physics, parameterization schemes, and parameters employed, land datasets determine the spatial variability in land surface states (i.e., soil moisture and temperature) and fluxes. Despite the importance of these datasets, they are often chosen out of convenience or regional limitations without due assessment of their impacts on model results. Irrigation is an anthropogenic form of land heterogeneity that has been shown to alter the land surface energy balance, ambient weather, and local circulations. As such, irrigation schemes are becoming more 15 prevalent in weather and climate models with rapid developments in dataset availability and parameterization scheme complexity. Thus, to address pragmatic issues related to modeling irrigation, this study uses a high-resolution, regional coupled modeling system to investigate the impacts of irrigation dataset selection on land-atmosphere (L-A) coupling using a case study from the Great Plains Irrigation Experiment (GRAINEX) field campaign. The simulations are assessed in the context of irrigated versus non-irrigated regions, subregions across the irrigation gradient, and sub-grid scale process 20 representation in coarser scale models. The results show that L-A coupling is sensitive to the choice of irrigation dataset and resolution and that the irrigation impact on surface fluxes and near surface meteorology can be dominant, conditioned on the details of the irrigation map (i.e., boundaries, heterogeneity, etc), or minimal. A consistent finding across several analyses was that even a low percentage of irrigation fraction (i.e., 4-16%) can have significant local and downstream atmospheric impacts (e.g., lower PBL height), suggesting that representation of boundaries and heterogeneous areas within irrigated 25 regions is particularly important for the modeling of irrigation impacts on the atmosphere in this model. When viewing the simulations presented here as a proxy for ‘ideal’ tiling in a Earth System Model scale gridbox, the results show that some ‘tiles’ will reach critical nonlinear moisture and planetary boundary layer (PBL) thresholds that could be important for clouds and convection, implying that heterogeneity resulting from irrigation should be taken into consideration in new sub-grid land-atmosphere exchange parameterizations.
In this study, accurately estimating methane (CH 4 ) flux is critically important for investigating and predicting the biogeochemistry-climate feedback. Better simulating CH 4 flux requires explicit representations of microbial processes on CH 4 dynamics because all processes for CH 4 production and consumption are actually carried out by microbes. A microbial functional group based module was developed and tested against an incubation experiment. The module considers four key mechanisms for CH 4 production and consumption: methanogenesis from acetate or single-carbon compounds and CH 4 oxidation using molecular oxygen or other inorganic electron acceptors. These four processes were carried out by four microbial functional groups: acetoclastic methanogens, hydrogenotrophic methanogens, aerobic methanotrophs, and anaerobic methanotrophs. This module was then linked with the decomposition subroutine of the Community Land Model , and was further used to simulate dynamics of carbon dioxide (CO 2 ) and CH 4 concentrations from an incubation experiment with permafrost soils. The results show that the model could capture the dynamics of CO 2 and CH 4 concentrations in microcosms with top soils, mineral layer soils and permafrost soils under natural and saturated moisture conditions and a temperature gradient of -2°C, 3°C, and 5°C. Sensitivity analysis confirmed the importance of acetic acid's direct contribution as substrate and indirect effects through pH feedback on CO 2 and CH 4 production and consumption. This study suggests that representing the microbial mechanisms is critical for modeling CH 4 production and consumption; it is urgent to incorporate microbial mechanisms into Earth system models for better predicting the behavior of the climate system.
In this two–part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100 km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part 1, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode—with prescribed sea surface temperatures (SSTs) and sea–ice distribution—is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. Here, the model's Cess sensitivity (response in the top–of–atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part 2, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.
Here, in Part 2 of this two–part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part 1. Part 2 provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.
In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.
The ARM Climate Research Facility, a DOE scientific user facility, provides the climate research community with strategically located in situ and remote sensing observatories designed to improve the understanding and representation, in climate and earth system models, of clouds and aerosols as well as their interactions and coupling with the Earth’s surface.
The study of climate change and its impacts depends on generating projections of future temperature and other climate variables. For detailed studies, these projections usually require some combination of numerical simulation and observations, given that simulations of even the current climate do not perfectly reproduce local conditions. We present a methodology for generating future climate projections that takes advantage of the emergence of climate model ensembles, whose large amounts of data allow for detailed modeling of the probability distribution of temperature or other climate variables. The procedure gives us estimated changes in model distributions that are then applied to observations to yield projections that preserve the spatiotemporal dependence in the observations. We use quantile regression to estimate a discrete set of quantiles of daily temperature as a function of seasonality and long-term change, with smooth spline functions of season, long-term trends, and their interactions used as basis functions for the quantile regression. A particular innovation is that more extreme quantiles are modeled as exceedances above less extreme quantiles in a nested fashion, so that the complexity of the model for exceedances decreases the further out into the tail of the distribution one goes. We apply this method to two large ensembles of model runs using the same forcing scenario, both based on versions of the Community Earth System Model (CESM), run at different resolutions. The approach generates observation-based future simulations with no processing or modeling of the observed climate needed other than a simple linear rescaling. The resulting quantile maps illuminate substantial differences between the climate model ensembles, including differences in warming in the Pacific Northwest that are particularly large in the lower quantiles during winter. We show how the availability of two ensembles allows the efficacy of the method to be tested with a “perfect model” approach, in which we estimate transformations using the lower-resolution ensemble and then apply the estimated transformations to single runs from the high-resolution ensemble. Finally, we describe and implement a simple method for adjusting a transformation estimated from a large ensemble of one climate model using only a single run of a second, but hopefully more realistic, climate model.
Abstract. Human land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25∘×0.25∘) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.
There is a growing need for high-resolution land surface parameters as land surface models are being applied at increasingly higher spatial resolution offline as well as in regional and global models. The default land surface parameters for the most recent version of the Community Land Model (i.e. CLM 4.0) are at 0.5° or coarser resolutions, released with the Community Earth System Model (CESM). Plant Functional Types (PFTs), vegetation properties such as Leaf Area Index (LAI), Stem Area Index (SAI), and non-vegetated land covers were developed using remotely sensed datasets retrieved in late 1990’s and the beginning of this century. In this study, we developed new land surface parameters for CLM 4.0, specifically PFTs, LAI, SAI and non-vegetated land cover composition, at 0.05° resolution globally based on the most recent MODIS land cover and improved MODIS LAI products. Compared to the current CLM 4.0 parameters, the new parameters produced a decreased coverage by bare soil and trees, but an increased coverage by shrub, grass, and cropland. The new parameters result in a decrease in global seasonal LAI, with the biggest decrease in boreal forests; however, the new parameters also show a large increase in LAI in tropical forest. Differences between the new and the current parameters are mainly caused by changes in the sources of remotely sensed data and the representation of land cover in the source data. Advantages and disadvantages of each dataset were discussed in order to provide guidance on the use of the data. The new high-resolution land surface parameters have been used in a coupled land-atmosphere model (WRF-CLM) applied to the western U.S. to demonstrate their use in high-resolution modeling. A remapping method from the latitude/longitude grid of the CLM data to the WRF grids with map projection was also demonstrated. Future work will include global offline CLM simulations to examine the impacts of source data resolution and subsequent land parameter changes on simulated land surface processes.
The Climate Model Output Rewriter (CMOR) software was first developed by LLNL’s PCMDI program in early 2000s and was formally released with v1.0 (July 2006), v2.0 (January 2011), and v3.1(June 2016). CMOR is used to produce Climate and Forecast Convention (http://cfconventions.org/) CF-compliant netCDF files, in the standard format required to satisfy the World Climate Research Program (WCRP) Coupled Model Intercomparison Project (CMIP). The software has been used across multiple phases of the Earth System Modeling (ESM) project CMIP (CMIP3, CMIP5, CMIP6, and planned use in CMIP7) along with numerous parallel projects focused on preparation observations for use in model evaluation (obs4MIPs) and forcing datasets (input4MIPs) to guide ESM simulations to meet strict experimental protocols. More information can be obtained from the CMOR website and code repositories: https://cmor.llnl.gov/; https://github.com/pcmdi/cmor; https://github.com/PCMDI/cmor3_documentation The ESM variable definitions used as input for CMOR can also be viewed in code repositories: https://github.com/PCMDI/cmip3-cmor-tables/; https://github.com/PCMDI/cmip5-cmor-tables/; https://github.com/PCMDI/cmip6-cmor-tables/