Dimensionality Reduction and Physics-Informed Neural Networks for Climate Land Models.
Abstract not provided.
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.
Abstract not provided.
Abstract not provided.
Explore the source record for details and available documents.
Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.
Modern scientific datasets present numerous data management and analysis challenges. State-of-the- art index and query technologies such as FastBit can significantly improve accesses to these datasets by augmenting the user data with indexes and other secondary information. However, a challenge is that the indexes assume the relational data model but the scientific data generally follows the array data model. To match the two data models, we design a generic mapping mechanism and implement an efficient input and output interface for reading and writing the data and their corresponding indexes. To take advantage of the emerging many-core architectures, we also develop a parallel strategy for indexing using threading technology. This approach complements our on-going MPI-based parallelization efforts. We demonstrate the flexibility of our software by applying it to two of the most commonly used scientific data formats, HDF5 and NetCDF. We present two case studies using data from a particle accelerator model and a global climate model. We also conducted a detailed performance study using these scientific datasets. The results show that FastQuery speeds up the query time by a factor of 2.5x to 50x, and it reduces the indexing time by a factor of 16 on 24 cores.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Despite the fact that Global Climate Model (GCM) outputs have been used to project hydrologic impacts of climate change using off-line hydrologic models for two decades, many of these efforts have been disjointed applications or at least calibrations have been focused on individual river basins and using a few of the available GCMs. This study improves upon earlier attempts by systematically projecting hydrologic impacts for the entire conterminous United States (US), using outputs from ten GCMs from the latest Coupled Model Intercomparison Project phase 5 (CMIP5) archive, with seamless hydrologic model calibration and validation techniques to produce a spatially and temporally consistent set of current hydrologic projections. The Variable Infiltration Capacity (VIC) model was forced with ten-member ensemble projections of precipitation and air temperature that were dynamically downscaled using a regional climate model (RegCM4) and bias-corrected to 1/24 (~4 km) grid resolution for the baseline (1966 2005) and future (2011 2050) periods under the Representative Concentration Pathway 8.5. Based on regional analysis, the VIC model projections indicate an increase in winter and spring total runoff due to increases in winter precipitation of up to 20% in most regions of the US. However, decreases in snow water equivalent (SWE) and snow-covered days will lead to significant decreases in summer runoff with more pronounced shifts in the time of occurrence of annual peak runoff projected over the eastern and western US. In contrast, the central US will experience year-round increases in total runoff, mostly associated with increases in both extreme high and low runoff. Furthermore, the projected hydrological changes described in this study have implications for various aspects of future water resource management, including water supply, flood and drought preparation, and reservoir operation.
An integrated chemistry-climate model (RegCM4-CHEM) simulates present-day climate, ozone and tropospheric aerosols over Egypt with a focus on Greater Cairo (GC) region. The densley populated GC region is known for its severe air quality issues driven by high levels of anthropogenic pollution in conjuction with natural sources such as dust and agricultural burning events. We find that current global emission inventories underestimate key pollutants such as nitrogen oxides and anthropogenic aerosol species. In the GC region, average-ground-based NO2 observations of 40-60 ppb are substantially higher than modeled estimates (5-10 ppb), likely due to model grid resolution, improper boundary layer representation, and poor emissions inventories. Observed ozone concentrations range from 35 ppb (winter) to 80 ppb (summer). The model reproduces the seasonal cycle fairly well, but modeled summer ozone is understimated by approximately 15 ppb and exhibits little interannual variability. For aerosols, springtime dust events dominate the seasonal aerosol cycle. The chemistry-climate model captures the springtime peak aerosol optical depth (AOD) of 0.7-1 but is slightly greater than satellite-derived AOD. Observed AOD decreases in the summer and increases again in the fall due to agricultural burning events in the Nile Delta, yet the model underestimates this fall observed AOD peak, as standard emissions inventories underestimate this burning and the resulting aerosol emissions. Our comparison of modeled gas and particulate phase atmospheric chemistry in the GC region indicates that improved emissions inventories of mobile sources and other anthropogenic activities are needed to improve air quality simulations in this region.
A primary goal of the Climate and Environmental Sciences Division (CESD) is to improve the predictive ability of regional and global climate models (GCMs). One set of tools used for evaluating and improving model predictions and projections are testbeds that combine data derived from observations and model simulations. In these testbeds, outputs of model simulations can be compared to observations of the climate system in order to identify errors and/or biases in the model simulations and determine the specific model processes that need improvements. Testbeds can also be used to provide scientific insights into dominant processes and process interactions as well as increase our understanding of the role of various physical processes involved in a particular case study or meteorological event. CESD currently funds multiple activities, which are developing testbeds to evaluate and improve atmospheric model parameterizations across a range of models, scales, and processes. For each of these testbeds, the Atmospheric Radiation Measurement (ARM) Climate Research Facility represents a primary source of observational data for evaluating the models. The Atmospheric System Research (ASR), Earth System Modeling (ESM), Regional and Global Climate Modeling (RGCM), and ARM program managers co-sponsored a workshop on atmospheric testbed activities across CESD. The workshop was held August 5-6, 2013, at the U. S. Department of Energy (DOE) Germantown facility. The principal objectives of the workshop were to survey the CESD testbed portfolios; identify the distinctive capabilities of each testbed as well as commonalities between the testbeds; identify areas for coordination and collaboration among testbeds; improve linkages between testbed activities and the ARM facility; identify gaps in testbed activities; and identify ways to improve coordination and communication among the testbeds and between the testbeds and the broader scientific community.
We use a nascent data-driven causal discovery method to find and compare causal relationships in observed data and climate model output. We consider ten different features in the Arctic climate collected from public databases on observational and Energy Exascale Earth System Model (E3SM) data. In identifying and analyzing the resulting causal networks, we make meaningful comparisons between observed and climate model interdependencies. This work demonstrates our ability to apply the PCMCI causal discovery algorithm to Arctic climate data, that there are noticeable similarities between observed and simulated Arctic climate dynamics, and that further work is needed to identify specific areas for improvement to better align models with natural observations.
Here, reproducing characteristics of observed sea ice extent remains an important climate modeling challenge. This study describes several approaches to improve how model biases in total sea ice distribution are quantified, and applies them to historically forced simulations contributed to phase 5 of the Coupled Model Intercomparison Project (CMIP5). The quantity of hemispheric total sea ice area, or some measure of its equatorward extent, is often used to evaluate model performance. A new approach is introduced that investigates additional details about the structure of model errors, with an aim to reduce the potential impact of compensating errors when gauging differences between simulated and observed sea ice. Using multiple observational datasets, several new methods are applied to evaluate the climatological spatial distribution and the annual cycle of sea ice cover in 41 CMIP5 models. It is shown that in some models, error compensation can be substantial, for example resulting from too much sea ice in one region and too little in another. Error compensation tends to be larger in models that agree more closely with the observed total sea ice area, which may result from model tuning. The results herein suggest that consideration of only the total hemispheric sea ice area or extent can be misleading when quantitatively comparing how well models agree with observations. Further work is needed to fully develop robust methods to holistically evaluate the ability of models to capture the finescale structure of sea ice characteristics; however, the “sector scale” metric used here aids in reducing the impact of compensating errors in hemispheric integrals.
We present a new version of the Community Earth System Model, version 1 (CESM1) with the Whole Atmosphere Community Climate Model (WACCM) featuring numerous improvements that are unique among earth system models. Improved horizontal resolution, dynamics, and chemistry now provide the development of an internally generated quasi-biennial oscillation, and significant improvements to temperatures and ozone throughout the stratosphere. The prognostic treatment of stratospheric sulfate aerosols is shown to represent well the evolution of stratospheric aerosol optical depth and perturbations to solar and longwave radiation following volcanic eruptions. We identify the inclusion of interactive OH chemistry as crucial to the study of aerosol formation following large inputs of SO 2 to the stratosphere. We show that depletion of OH levels within the dense SO 2 cloud in the first weeks following the June 1991 eruption of Mt. Pinatubo significantly prolonged the e-folding decay time for SO 2 oxidation to 47 days. Previous observational and model studies showing a 30-day decay time have not accounted for the large initial losses of SO 2 on ash and ice in the first 7-9 days following the eruption, and have not correctly accounted for OH depletion. The completeness of the chemistry, dynamics, and aerosol microphysics in WACCM uniquely qualify it for studies of stratospheric sulfate aerosol geoengineering.
Abstract Model evidence for the “pattern effect” assumes that global climate models (GCMs) faithfully simulate how clouds respond to varying sea surface temperature (SST) patterns and associated meteorological perturbations. We exploit time-invariant satellite-based estimates of the sensitivity of marine low clouds to meteorological perturbations to estimate how these clouds responded to time-varying SST patterns and meteorology between 1870 and 2014. GCMs and reanalyses provide estimates of the historical meteorological changes. Observations suggest that increasing estimated inversion strength (EIS) between 1980 and 2014 produced a negative low cloud feedback, opposite to the positive feedback expected from increasing CO 2 . This indicates that the processes responsible for marine cloud changes from 1980 to the near present are distinct from those associated with an increase in CO 2 . We also observationally constrain the difference between the historical near-global marine low cloud feedback, , and that arising from increasing CO 2 , . We find that this cloud feedback pattern effect depends strongly on time period and reanalysis dataset, and that varying changes in EIS and SST with warming explain much of its variability. Between 1980 and 2014, we estimate that (90% confidence) assuming meteorological changes from the Multiple Reanalysis Ensemble, implying a total pattern effect (that arising from all climate feedbacks) of 1.86 ± 0.45 W m −2 K −1 . This observational evidence corroborates previous quantitative estimates of the pattern effect, which heretofore relied largely upon GCM-based cloud changes. However, disparate historical meteorological changes across individual reanalyses contribute to considerable uncertainty in its magnitude.
This study explores the added value (AV) of a regional Earth system model (ESM) compared to an atmosphere-only regional climate model (RCM) in simulating West African monsoon (WAM) rainfall. The primary goals are to foster discussions on the suitability of coupled RCMs for WAM projections and deepen our understanding of ocean–atmosphere coupling’s influence on the WAM system. The study employs results from dynamical downscaling of the ERA-Interim reanalysis and Max Plank Institute ESM, low resolution (MPI-ESM-LR), by two RCMs, atmosphere only (REMO) and REMO coupled with Max Planck Institute Ocean Model (MPIOM) (ROM), at ~25-km horizontal resolution. Results show that in regions distant from coupling domain boundaries such as West Africa (WA), constraint conditions from ERA-Interim are more beneficial than coupling effects. REMO, reliant on oceanic sea surface temperatures (SSTs) from observations and influenced by ERA-Interim, is biased under coupling conditions, although coupling offers potential advantages in representing heat and mass fluxes. Contrastingly, as intended, coupling improves SSTs and monsoon fluxes’ relationships under ESM-forced conditions. In this latter case, the coupling features a dipole-like spatial structure of AV, improving precipitation over the Guinea Coast but degrading precipitation over half of the Sahel. Our extensive examination of physical processes and mechanisms underpinning the WAM system supports the plausibility of AV. Additionally, we found that the monsoonal dynamics over the ocean respond to convective activity, with the Sahara–Sahel surface temperature gradient serving as the maintenance mechanism. While further efforts are needed to enhance the coupled RCM, we advocate for its use in the context of WAM rainfall forecasts and projections.
Abstract Large ensembles of model simulations are frequently used to reduce the impact of internal variability when evaluating climate models and assessing climate change induced trends. However, the optimal number of ensemble members required to distinguish model biases and climate change signals from internal variability varies across models and metrics. Here we analyze the mean, variance and skewness of precipitation and sea surface temperature in the eastern equatorial Pacific region often used to describe the El Niño–Southern Oscillation (ENSO), obtained from large ensembles of Coupled model intercomparison project phase 6 climate simulations. Leveraging established statistical theory, we develop and assess equations to estimate, a priori, the ensemble size or simulation length required to limit sampling‐based uncertainties in ENSO statistics to within a desired tolerance. Our results confirm that the uncertainty of these statistics decreases with the square root of the time series length and/or ensemble size. Moreover, we demonstrate that uncertainties of these statistics are generally comparable when computed using either pre‐industrial control or historical runs. This suggests that pre‐industrial runs can sometimes be used to estimate the expected uncertainty of statistics computed from an existing historical member or ensemble, and the number of simulation years (run duration and/or ensemble size) required to adequately characterize the statistic. This advance allows us to use existing simulations (e.g., control runs that are performed during model development) to design ensembles that can sufficiently limit diagnostic uncertainties arising from simulated internal variability. These results may well be applicable to variables and regions beyond ENSO.
The land ice contribution to global mean sea level rise has not yet been predicted using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models, but primarily used previous-generation scenarios and climate models, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios using statistical emulation of the ice sheet and glacier models. In this work, we find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.