Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical Downscaling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Assessment of Data Assimilation with the Prototype High Resolution Rapid Refresh for Alaska (HRRRAK)

The Arctic Region Supercomputing Center has been running a quasi-operational prototype of a High Resolution Rapid Refresh for Alaska (HRRRAK) at 3km resolution, initialized by the 13km Rapid Refresh (RR). Although the RR assimilates a broad range of observations into its analyses, experiments with the HRRRAK suggest that there may be added value in assimilating observations into the 3km initial conditions, downscaled from the 13km RR analyses. The NASA Short-term Prediction Research and Transition (SPoRT) group has been using assimilated data from the Atmospheric Infrared Sounder (AIRS) in WRF and WRF-Var simulations since 2004 with promising results. The sounder is aboard NASA s Aqua satellite, and provides vertical profiles of temperature and humidity. The Gridpoint Statistical Interpolation (GSI) system is then used to assimilate these vertical profiles into WRF forecasts. In this work, we assess the use of AIRS data in combination with other global data assimilation products on non-assimilated HRRRAK case studies. Two separate weather events will be assessed to qualitatively and quantitatively assess the impacts of AIRS data on HRRRAK forecasts.

Harrison, Kayla↗
DOE OSTI2025

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Decline in Seasonal Snow during a Projected 20-Year Dry Spell

Snowpack loss in midlatitude mountains is ubiquitously projected by Earth system models, though the magnitudes, persistence, and time horizons of decline vary. Using daily downscaled hydroclimate and snow projections, we examine changes in snow seasonality across the U.S. Pacific Southwest region during a simulated severe 20-year dry spell in the 21st century (2051–2070) developed as part of the 4th California Climate Change Assessment to provide a “stress test” for water resources. Across California’s mountains, substantial declines (30–100% loss) in median peak annual snow water equivalent accompany changes in snow seasonality throughout the region compared to the historic period. We find that 80% of historic seasonal snowpacks transition to ephemeral conditions. Subsetting empirical-statistical wildfire projections for California by snow seasonality transition regions indicates a two-to-four-fold increase in the area burned, consistent with recent observations of high elevation wildfires following extended drought conditions. By analyzing six of the major California snow-fed river systems, we demonstrate snowpack reductions and seasonality transitions result in concomitant declines in annual runoff (47–58% of historical values). The negative impacts to statewide water supply reliability by the projected dry spell will likely be magnified by changes in snowpack seasonality and increased wildfire activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Global Corn Heat Stress: Mean and SD of Degree Days Above 29°C based on NEX-GDDP-CMIP6 Climate Projections

Description This global dataset provides the estimated mean and standard deviation (SD) of corn heat stress (degree days above 29°C) for a set of climate models in NEX-GDDP-CMIP6 at 0.25-degree resolution. The NEX-GDDP-CMIP6 dataset is comprised of global downscaled climate scenarios derived from the General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 6 (CMIP6). The current dataset includes: Long-Term Average Degree Days Above 29°C- Historical Long-Term Average Degree Days Above 29°C- SSP245 Long-Term Standard Deviation of Degree Days Above 29°C- Historical Long-Term Standard Deviation of Degree Days Above 29°C- SSP245 The mean and SD are calculated over 1985-2014 for the historical period and over 2035-2064 for future projections. A full description of methods, including growing season, daily temperature distribution, and statistical coefficients, can be found in Haqiqi (2024). The source climate data are obtained from https://ds.nccs.nasa.gov/thredds2/catalog/catalog.html and are described in Thrasher et al (2022). The codes used to create this dataset are available at https://github.com/ihaqiqi/dd29c_nex_cmip6. Acknowledgments This work was supported by the US Department of Energy, Office of Science, Biological and Environmental Research Program, Earth and Environmental Systems Modeling, MultiSector Dynamics under Cooperative Agreement DE-SC0022141. The data processing, computation, and storage were completed on Purdue Anvil supercomputer and cyberinfrastructure supported by the National Science Foundation HDR award # 2118329: "NSF Institute for Geospatial Understanding through an Integrative Discovery Environment (I-GUIDE)". References Haqiqi. I. (2024). Trade can buffer climate-induced risks and volatilities in crop supply. Environmental Research: Food Systems. https://doi.org/10.1088/2976-601X/ad7d12 Thrasher, B., Wang, W., Michaelis, A., Melton, F., Lee, T., & Nemani, R. (2022). NASA global daily downscaled projections, CMIP6. Scientific Data, 9(1), 262. https://doi.org/10.1038/s41597-022-01393-4

Climate Change↗
NASA NTRS2014

Reply to Comment by Laprise on 'the Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: a Case Study over the Continental U.S. Using the GISS-ModelE2 and WRF Models'

In his comment, Laprise raises several points that we agree merit consideration. His primary critique is that our study [Racherla et al., 2012] tested the ability of the WRF regional climate model to reproduce historical temperature and precipitation change relative to the driving global climate model (GCM) using only a single simulation rather than an ensemble. He asserts that the observed changes are smaller than the internal variability in the climate system (i.e., not statistically significant) and that thus a single simulation should not necessarily be able to capture the observations. Laprise points out that the statistical signal is reduced for a multi-decadal trend such as the one we analyzed in comparison with mean climatology and cites two studies showing that for particular climate parameters it can take any years for a signal to be discerned over internal variability. He states that The results of theexperiment as designed were strongly influenced by the presence of internal variability and sampling errors,which masked the rather small climate changes that may have occurred as a consequence of changes inforcing during the period considered. While Laprise discusses statistics in general terms at some length, for the actual climate trends examined in our study, he offers no evidence that the forced signal was smallcompared with internal variability. The two studies he cites [de Ela et al., 2013; Maraun, 2013] do not provide convincing evidence as they concern climate variables averaged over different times and areas. One in fact examines extreme precipitation events, which by definition are rare and thus have a lower significance level. We accept the general point that it is important to consider internal variability, and as noted in our paper we agree that an ensemble of simulations is in principle an optimal, though computationally expensive, approach. While we did not present the statistical significance of the observations in our original paper, we have now evaluated those for the regional temperature trends used in our study to evaluate the added value of WRF and thus can analyze data as to the magnitude of the trends with respect to internal variability.

climate↗
NASA NTRS1993

The helical decomposition and the instability assumption

Direct numerical simulations show that the triadic transfer function T(k,p,q) peaks sharply when q (or p) is much smaller than k. The triadic transfer function T(k,p,q) gives the rate of energy input into wave number k from all interactions with modes of wave number p and q, where k, p, q form a triangle. This observation was thought to suggest that energy is cascaded downscale through non-local interactions with local transfer and that there was a strong connection between large and small scales. Both suggestions were in contradiction with the classical Kolmogorov picture of the energy cascade. The helical decomposition was found useful in distinguishing between kinematically independent interactions. That analysis has gone beyond the question of non-local interaction with local transfer. In particular, an assumption about the statistical direction of triadic energy transfer in any kinematically independent interaction was introduced (the instability assumption). That assumption is not necessary for the conclusions about non-local interactions with local transfer recalled above. In the case of turbulence under rapid rotation, the instability assumption leads to the prediction that energy is transferred in spectral space from the poles of the rotation axis toward the equator. The instability assumption is thought to be of general validity for any type of triad interactions (e.g. internal waves). The helical decomposition and the instability assumption offer detailed information about the homogeneous statistical dynamics of the Navier-Stokes equations. The objective was to explore the validity of the instability assumption and to study the contributions of the various types of helical interactions to the energy cascade and the subgrid-scale eddy-viscosity. This was done in the context of spectral closures of the Direct Interaction or Quasi-Normal type.

Waleffe, Fabian A.↗
DOE OSTI2020

A global ensemble of ocean wave climate projections from CMIP5-driven models

This dataset, produced through the Coordinated Ocean Wave Climate Project (COWCLIP) phase 2, represents the first coordinated multivariate ensemble of 21st Century global wind-wave climate projections available (henceforth COWCLIP2.0). COWCLIP2.0 comprises general and extreme statistics of significant wave height (H S ), mean wave period (T m ), and mean wave direction (θ m ) computed over time-slices 1979–2004 and 2081–2100, at different frequency resolutions (monthly, seasonally and annually). The full ensemble comprising 155 global wave climate simulations is obtained from ten CMIP5-based state-of-the-art wave climate studies and provides data derived from alternative wind-wave downscaling methods, and different climate-model forcing and future emissions scenarios. The data has been produced, and processed, under a specific framework for consistency and quality, and follows CMIP5 Data Reference Syntax, Directory structures, and Metadata requirements. Technical comparison of model skill against 26 years of global satellite measurements of significant wave height has been undertaken at global and regional scales. This new dataset provides support for future broad scale coastal hazard and vulnerability assessments and climate adaptation studies in many offshore and coastal engineering applications.

Science & Technology - Other Topics↗
DOE OSTI2023

Dynamically Downscaled Hourly Future Weather Data with 12-km Resolution Covering Most of North America

This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.

Array↗
DOE OSTI2021

Future projections in tropical cyclone activity over multiple CORDEX domains from RegCM4 CORDEX-CORE simulations

The characteristics of tropical cyclone (TC) activity over five TC basins lying within four Coordinated Regional Downscaling Experiment (CORDEX) domains are examined for present and future climate conditions using a new ensemble of simulations completed as part of the CORDEX-CORE initiative with the regional climate model RegCM4. The simulations are conducted at a 25 km horizontal grid spacing and are driven by three CMIP5 general circulation models (GCMs) under two Representative Concentration Pathways (RCP2.6 and RCP8.5). The RegCM4 captures most features of the observed TC climatology, except for the TC intensity, which is thus statistically adjusted using a bias correction procedure to account for the effect of the coarse model resolution. The RegCM4 exhibits an improved simulation of several TC statistics compared to the driving GCMs, over most basins analyzed. In future climate conditions we find significant increases in TC frequency over the North Indian Ocean, the Northwest Pacific and Eastern Pacific regions, which are consistent with an increase in mid-tropospheric relative humidity. The North Atlantic and Australasia regions show a decrease in TC frequency, mostly associated with an increase in wind shear. We also find a consistent increase in future storm rainfall rates associated with TCs and in the frequency of the most intense TCs over most domains. Finally, our study shows robust responses often, but not always, in line with previous studies, still implying the presence of significant uncertainties in the projection of TC characteristics, which need to be addressed using large ensembles of simulations with high-resolution models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Antarctic ice sheet model comparison with uncurated geological constraints shows that higher spatial resolution improves deglacial reconstructions

Accurately reconstructing past changes to the shape and volume of the Antarctic ice sheet relies on the use of physically based and thus internally consistent ice sheet modeling, benchmarked against spatially limited geologic data. The challenge in model benchmarking against geologic data is diagnosing whether model-data misfits are the result of an inadequate model, inherently noisy or biased geologic data, and/or incorrect association between modeled quantities and geologic observations. In this work we address this challenge by (i) the development and use of a new model-data evaluation framework applied to an uncurated data set of geologic constraints, and (ii) nested high-spatial-resolution modeling designed to test the hypothesis that model resolution is an important limitation in matching geologic data. While previous approaches to model benchmarking employed highly curated datasets, our approach applies an automated screening and quality control algorithm to an uncurated public dataset of geochronological observations (specifically, cosmogenic-nuclide exposure-age measurements from glacial deposits in ice-free areas). This optimizes data utilization by including more geological constraints, reduces potential interpretive bias, and allows unsupervised assimilation of new data as they are collected. We also incorporate a nested model framework in which high-resolution domains are downscaled from a continent-wide ice sheet model. We highlight the application of this framework by applying these methods to a small ensemble of deglacial ice-sheet model simulations, and demonstrate that the nested approach improves the ability of model simulations to match exposure age data collected from areas of complex topography and ice flow. We develop a range of diagnostic model-data comparison metrics to provide more insight into model performance than possible from a single-valued misfit statistic, showing that different metrics capture different aspects of ice sheet deflation.

Geosciences↗
DOE OSTI2025

Upper bounds for 21st-century surface air temperatures in the Western United States

The last decade has seen a large number of severe heatwaves that were unprecedented in the observational record, highlighting challenges associated with observationally-based statistical quantification of the likelihood and magnitude of future extreme temperatures. An alternative to such probabilistic assessments is identification of upper bounds that quantify the hottest surface air temperatures that can possibly be achieved by the end of the 21st century. Theory, simulations, and observational analyses support the existence of a finite upper bound for surface air temperature; however, estimates for future upper-bound values that are realistic and usable for planning remain unavailable. Here, we combine atmospheric theory with large ensembles of dynamically downscaled projections to estimate historical and end-of-century upper bounds for surface air temperatures. A number of physical mechanisms can influence upper bounds, and at the end of the 21st century, estimates based on mechanisms that yield more moderate upper-bounds produce values around 60∘C for much of the Western United States and in excess of 80∘C for the hottest parts of the domain. Even cooler high-altitude locations have end-of-century upper bounds over 50∘C. Although these upper-bound estimates might seem implausibly large, increases in the upper bounds over the 21st century are similar to increases in dynamically downscaled peak surface temperatures after adjusting those downscaled temperatures to eliminate the possibly biased model trends in surface specific humidity. While upper bound estimates are high relative to historical observations, they nonetheless suggest that heatwave intensity risk is bounded, with uncertainty dominated by projections of surface and upper-level humidity.

Risser, Mark D↗
DOE OSTI2024

The first ensemble of kilometer-scale simulations of a hydrological year over the third pole

An accurate understanding of the current and future water cycle over the Third Pole is of great societal importance, given the role this region plays as a water tower for densely populated areas downstream. An emerging and promising approach for skillful climate assessments over regions of complex terrain is kilometer-scale climate modeling. As a foundational step towards such simulations over the Third Pole, we present a multi-model and multi-physics ensemble of kilometer-scale regional simulations for the hydrological year of October 2019 to September 2020. The ensemble consists of 13 simulations performed by an international consortium of 10 research groups, configured with a horizontal grid spacing ranging from 2.2 to 4 km covering all of the Third Pole region. These simulations are driven by ERA5 and are part of a Coordinated Regional Climate Downscaling EXperiment Flagship Pilot Study on Convection-Permitting Third Pole. The simulations are compared against available gridded and in-situ observations and remote-sensing data, to assess the performance and spread of the model ensemble compared to the driving reanalysis during the cold and warm seasons. Although ensemble evaluation is hindered by large differences between the gridded precipitation datasets used as a reference over this region, we show that the ensemble improves on many warm-season precipitation metrics compared with ERA5, including most wet-day and hour statistics, and also adds value in the representation of wet spells in both seasons. As such, the ensemble will provide an invaluable resource for future improvements in the process understanding of the hydroclimate of this remote but important region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Examining the Impacts of Great Lakes Temperature Perturbations on Simulated Precipitation in the Northeastern United States

Most inland water bodies are not resolved by general circulation models, requiring that lake surface temperatures be estimated. Given the large spatial and temporal variability of the surface temperatures of the North American Great Lakes, such estimations can introduce errors when used as lower boundary conditions for dynamical downscaling. Lake surface temperatures (LSTs) influence moisture and heat fluxes, thus impacting precipitation within the immediate region and potentially in regions downwind of the lakes. For this study, the Advanced Research version of the Weather Research and Forecasting Model (WRF-ARW) was used to simulate precipitation over the six New England states during a 5-yr historical period. The model simulation was repeated with perturbed LSTs, ranging from 10°C below to 10°C above baseline values obtained from reanalysis data, to determine whether the inclusion of erroneous LST values has an impact on simulated precipitation and synoptic-scale features. Results show that simulated precipitation in New England is statistically correlated with LST perturbations, but this region falls on a wet–dry line of a larger bimodal distribution. Wetter conditions occur to the north and drier conditions occur to the south with increasing LSTs, particularly during the warm season. Additionally, the precipitation differences coincide with large-scale anomalous temperature, pressure, and moisture patterns. Care must therefore be taken to ensure reasonably accurate Great Lakes surface temperatures when simulating precipitation, especially in southeastern Canada, Maine, and the mid-Atlantic region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Drought Indicators Based on Model Assimilated GRACE Terrestrial Water Storage Observations

The Gravity Recovery and Climate Experiment (GRACE) twin satellites observe time variations in Earth's gravity field which yield valuable information about changes in terrestrial water storage (TWS). GRACE is characterized by low spatial (greater than 150,000 square kilometers) and temporal (greater than 10 day) resolution but has the unique ability to sense water stored at all levels (including groundwater) systematically and continuously. The GRACE Data Assimilation System (GRACE-DAS), based on the Catchment Land Surface Model (CLSM) enhances the value of the GRACE water storage data by enabling spatial and temporal downscaling and vertical decomposition into moisture 39 components (i.e. groundwater, soil moisture, snow), which individually are more useful for scientific applications. In this study, GRACE-DAS was applied to North America and GRACE-based drought indicators were developed as part of a larger effort that investigates the possibility of more comprehensive and objective identification of drought conditions by integrating spatially, temporally and vertically disaggregated GRACE data into the U.S. and North American Drought Monitors. Previously, the Drought Monitors lacked objective information on deep soil moisture and groundwater conditions, which are useful indicators of drought. Extensive datasets of groundwater storage from USGS monitoring wells and soil moisture from the Soil Climate Analysis Network (SCAN) were used to assess improvements in the hydrological modeling skill resulting from the assimilation of GRACE TWS data. The results point toward modest, but statistically significant, improvements in the hydrological modeling skill across major parts of the United States, highlighting the potential value of GRACE assimilated water storage field for improving drought detection.

Houborg, Rasmus↗
DOE OSTI2023

Climate Change Signal in Atlantic Tropical Cyclones Today and Near Future

This manuscript discusses the challenges in detecting and attributing recently observed trends in the Atlantic tropical cyclone (TC) and the epistemic uncertainty we face in assessing future risk. We use synthetic storms downscaled from five CMIP5 models by the Columbia HAZard model (CHAZ), and directly simulated storms from high-resolution climate models. We examine three aspects of recent TC activity: the upward trend and multi-decadal oscillation of the annual frequency, the increase in storm wind intensity, and the decrease in forward speed. Some data sets suggest that these trends and oscillation are forced while others suggest that they can be explained by natural variability. Projections under warming climate scenarios also show a wide range of possibilities, especially for the annual frequencies, which increase or decrease depending on the choice of moisture variable used in the CHAZ model and on the choice of climate model. The uncertainties in the annual frequency lead to epistemic uncertainties in TC risk assessment. Here, we investigate the potential for reduction of these epistemic uncertainties through a statistical practice, namely likelihood analysis. We find that historical observations are more consistent with the simulations with increasing frequency than those with decreasing frequency, but we are not able to rule out the latter. We argue that the most rational way to treat epistemic uncertainty is to consider all outcomes contained in the results. In the context of risk assessment, since the results contain possible outcomes in which TC risk is increasing, this view implies that the risk is increasing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subseasonal Clustering of Atmospheric Rivers Over the Western United States

Abstract The serial occurrence of atmospheric rivers (ARs) along the US West Coast can lead to prolonged and exacerbated hydrologic impacts, threatening flood‐control and water‐supply infrastructure due to soil saturation and diminished recovery time between storms. Here a statistical approach for quantifying subseasonal temporal clustering among extreme events is applied to a 41‐year (1979–2019) wintertime AR catalog across the western United States (US). Observed AR occurrence, compared against a randomly distributed AR timeseries with the same average event density, reveals temporal clustering at a greater‐than‐random rate across the western US with a distinct geographical pattern. Compared to the Pacific Northwest, significant AR clusters over the northern Coastal Range of California and Sierra Nevada are more frequent and occur over longer time periods. Clusters along the California Coastal Range typically persist for 2 weeks, are composed of 4–5 ARs per cluster, and account for over 85% of total AR occurrence. Across the northwest Coast‐Cascade Ranges, clusters account for ∼50% of total AR occurrence, typically last 8–10 days, and contain 3–4 individual AR events. Based on precipitation data from a high‐resolution dynamical downscaling of reanalysis, the fractions of total and extreme hourly precipitation attributable to AR clusters are largest along the northern California coast and in the Sierra Nevada. Interannual variability among clusters highlights their importance for determining whether a particular water year is anomalously wet or dry. The mechanisms behind this unusual clustering are unclear and require further research.

Meteorology & Atmospheric Sciences↗
NASA NTRS2015

Evaluation of the 7-km GEOS-5 Nature Run

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

GEOS-5↗
DOE OSTI2021

A Binomial Stochastic Framework for Efficiently Modeling Discrete Statistics of Convective Populations

Abstract Understanding the coupling between convective clouds and the general circulation, as well as addressing the gray zone problem in convective parameterization, requires insight into the genesis and maintenance of spatial patterns in cumulus cloud populations. In this study, a simple toy model for recreating populations of interacting convective objects as distributed over a two‐dimensional Eulerian grid is formulated to this purpose. Key elements at the foundation of the model include i) a fully discrete formulation for capturing discrete behavior in convective properties at small population sample sizes, ii) object age‐dependence for representing life‐cycle effects, and iii) a prognostic number budget allowing for object interactions and co‐existence of multiple species. A primary goal is to optimize the computational efficiency of this system. To this purpose the object birth rate is represented stochastically through a spatially aware Bernoulli process. The same binomial stochastic operator is applied to horizontal advection of objects, conserving discreteness in object number. The applicability to atmospheric convection as well as behavior implied by the formulation is assessed. Various simple applications of the BiOMi model (Binomial Objects on Microgrids) are explored, suggesting that important convective behavior can be captured at low computational cost. This includes i) subsampling effects and associated powerlaw scaling in the convective gray zone, ii) stochastic predator‐prey behavior, iii) the downscale turbulent energy cascade, and iv) simple forms of spatial organization and convective memory. Consequences and opportunities for convective parameterization in next‐generation weather and climate models are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗