Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble simulation”

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 199 records · Page 11

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

Soil Moisture Active Passive (SMAP) mission brightness temperature (T(b) ) observations are assimilated into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update simulations of surface and root-zone soil moisture. Different time series components of the T(b) observations are assimilated including anomalies, inter-annual variations, and high frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the uncertainties of modeled and observed T(b) is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. In a benchmark experiment, T(b) observations are assimilated using a spatially constant 4 Kelvin (K) observation uncertainty, as in the operational SMAP Level-4 algorithm. All DA experiments exhibit notable skill improvements in most regions. Improvements are largest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for surface soil moisture simulations (0.06). No notable difference in skill - neither improvement nor deterioration - is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. This may be explained by the presence of large observation operator errors, which are analytically shown to have the potential to render post-update uncertainty insensitive to inaccuracies in estimates of the Kalman gain. These results have important implications for the design of soil moisture DA systems, in particular for parameterizing model and observation uncertainties.

Hydrology↗

Impacts of Stratospheric Ozone Recovery on Southern Ocean Temperature and Heat Budget

The impacts of stratospheric ozone recovery on Southern Ocean surface and interior temperature, heat content, heat uptake, and heat transport are investigated by contrasting two ensemble chemistry-climate model simulations in 2005–2099: one with fixed ozone depleting substances (ODSs) and another with decreasing ODSs. In our simulations ozone recovery significantly affects Southern Ocean temperature, with large latitudinal and vertical variations. Ozone recovery causes a dipole change of the full-depth ocean heat content (OHC) with an increase south of 60°S and a decrease between 45°S and 60°S. Integrated over latitudes south of 40°S, OHC decreases in response to ozone recovery. This ocean heat loss is shown to be driven by weakened poleward ocean heat transport (OHT) across 40°S, which is partly canceled by enhanced heat uptake. The weakening of poleward OHT into the Southern Ocean is caused by the ozone-induced equatorward shift of the meridional overturning circulation.

stratospheric ozone↗

Evidence for Limited Indirect Aerosol Forcing in Stratocumulus

Increases in cloud cover and condensed water contribute more than half of the indirect aerosol effect in an ensemble of general circulation model (GCM) simulations estimating the global radiative forcing of anthropogenic aerosols. We use detailed simulations of marine stratocumulus clouds and airborne observations of ship tracks to show that increases in cloud cover and condensed water in reality are far less than represented by the GCM ensemble. Our results offer an explanation for recent simplified inverse climate calculations indicating that indirect aerosol effects are greatly exaggerated in GCMs.

Ackerman, Andrew S.↗

Tropical impacts of SST forcing: A case study for 1987 versus 1988

The response of the NASA/Goddard Institute for Space Studies general circulation model (GCM) to large tropical sea surface temperature (SST) anomalies is investigated by evaluating model simulations of the particularly contrasting summer monsoon seasons 1987 and 1988. These years are representative of the warm and cold phases, respectively, of a recent El Nino-Southern Oscillation (ENSO) event. An ensemble averaging the results of three simulations was considered for each season, using monthly mean observed SST anomalies for June-August 1987 and 1988 as lower boundary forcing. Consistent with the European Center for Medium Weather Forecast (ECMWF)-analyzed winds, the simulations based on 1988 as compared to 1987 SST exhibit stronger upper-tropospheric irrotational circulation between the monsoon regions and the Southern Hemispheric sub-tropical anticyclones, a stronger Pacific Walker cell and a weaker subtropical westerly jet over the South Pacific. In the same vein, the modeled precipitation, indicating a more northerly position of the Pacific Inter-Tropical Convergence Zone (ITCZ) in 1988 compared with 1987, is supported by satellite observations of outgoing longwave radiation and highly reflective clouds.

Druyan, Leonard M.↗

Additional Contributions to the Development of the New Snow-Physics Scheme for SSiB

The Simplified Simple Biosphere Model (SSiB) had a well-documented problem with snowmelt timing and infiltration. A new snow-physics scheme was developed for use in SSIB. In this, the snow layer is separated from the soil, with its own energy budget and temperature. Solar energy reaching the top of the snowpack is divided into three parts: one, reflected by the snow; two, absorbed by the snow; and three, transmitted to the ground following a simple extinction relation. Heat is exchanged between the ground and snow by conduction and by radiation through an arbitrary air-gap between them. In the GSWP exercise using the GEWEX ISLSCP Initiative I forcing data (hereafter "offline"), it was found that the new snow scheme ameliorated a significant fraction of snowmelt time-delay as compared to observations from satellite. It also produced warmer ground temperatures under the snowpack, which allowed realistic meltwater infiltration, resulting in better simulated spring soil moisture recharge and peak runoff amount as compared to observations. An ensemble of six June-July-August (JJA) simulations for 1987 and 1988 were performed with the NASA Goddard GEOS II GCM coupled with the new snow-physics SSIB using new initial soil moisture (ISM) from the offline simulations. The GCM produced more realistic precipitation in northern regions that had large snowmelt and wetter ISM in response to better snow-physics, as compared to simulations with ISM without the new snow scheme. The new SSiB-GCM also increased the interannual precipitation signal in the Indian monsoon region, resulting from changes in ISM in the Himalayas and central Asia.

Mocko, David M.↗

Uncertainty in Twenty-First-Century CMIP5 Sea Level Projections

The representative concentration pathway (RCP) simulations included in phase 5 of the Coupled Model Intercomparison Project (CMIP5) quantify the response of the climate system to different natural and anthropogenic forcing scenarios. These simulations differ because of 1) forcing, 2) the representation of the climate system in atmosphere-ocean general circulation models (AOGCMs), and 3) the presence of unforced (internal) variability. Global and local sea level rise projections derived from these simulations, and the emergence of distinct responses to the four RCPs depend on the relative magnitude of these sources of uncertainty at different lead times. Here, the uncertainty in CMIP5 projections of sea level is partitioned at global and local scales, using a 164-member ensemble of twenty-first-century simulations. Local projections at New York City (NYSL) are highlighted. The partition between model uncertainty, scenario uncertainty, and internal variability in global mean sea level (GMSL) is qualitatively consistent with that of surface air temperature, with model uncertainty dominant for most of the twenty-first century. Locally, model uncertainty is dominant through 2100, with maxima in the North Atlantic and the Arctic Ocean. The model spread is driven largely by 4 of the 16 AOGCMs in the ensemble; these models exhibit outlying behavior in all RCPs and in both GMSL and NYSL. The magnitude of internal variability varies widely by location and across models, leading to differences of several decades in the local emergence of RCPs. The AOGCM spread, and its sensitivity to model exclusion and/or weighting, has important implications for sea level assessments, especially if a local risk management approach is utilized.

Variability↗

Causes of Long-Term Drought in the United States Great Plains

The United States Great Plains (USGP) experienced a number of multi-year droughts during the last century, most notably the droughts of the 1930s and 1950s. This study examines the causes of such droughts using ensembles of long term (1930-1999) simulations carried out with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed sea surface temperatures (SSTs). The results show that the model produces long-term (multi-year) variations in the USGP precipitation that are similar to those observed. A correlative analysis suggests that the ensemble mean low frequency (time scales longer than about 6 years) rainfall variations in the USGP are linked to a pan-Pacific pattern of SST variability that is the leading empirical orthogonal function (EOF) in the low frequency SST data. The link between the SST and the Great Plains precipitation is confirmed in idealized AGCM simulations, in which the model is forced by the 2 polarities of the pan-Pacific SST pattern. The idealized simulations further show that it is primarily the tropical part of the SST anomalies that influence the USGP. As such, the USGP tend to have above normal precipitation when the tropical Pacific SSTs are above normal, while there is a tendency for drought when the tropical SSTs are cold. The upper tropospheric response to the pan-Pacific SST EOF shows a global-scale pattern with a strong wave response in the Pacific and a substantial zonally-symmetric component in which USGP pluvial (drought) conditions are associated with reduced (enhanced) heights throughout the extra-tropics. The potential predictability of rainfall in the USGP associated with SSTs is rather modest, with on average about 1/3 of the total low frequency rainfall variance forced by SST anomalies. Further idealized experiments with climatological SST, suggest that the remaining low frequency variance in the USGP precipitation is the result of interactions with soil moisture. In particular, simulations with soil moisture feedback show a six-fold increase in the variance in annual USGP precipitation compared with simulations in which the soil feedback is excluded. In addition to increasing variance, the interactions with the soil introduce year-to-year memory in the hydrological cycle that is consistent with a red noise process, in which the low frequencies in the deep soil are the result of integrating a net forcing (precipitation-evaporation-runoff) that is white noise on interannual time scales. As such, the role of low frequency SST variability is to introduce a bias to the net forcing on the soil moisture that drives the random process preferentially to either wet or dry conditions.

Schubert, Siegfried D.↗

Ensemble averaged dynamic modeling

The possibility of using the information from simultaneous equivalent Large Eddy Simulations (LAS) for improving the subgrid scale modeling is investigated. An ensemble average dynamic model is proposed as an alternative to the usual spatial average versions. It is shown to be suitable independently of the existence of any homogeneity directions, and its formulation is thus universal. The ensemble average dynamic model is shown to give very encouraging results for as few as 16 simultaneous LES's.

Carati, D.↗

Precipitation Processes Developed During ARM (1997), TOGA COARE (1992) GATE (1974), SCSMEX (1998), and KWAJEX (1999): Consistent 3D, Semi-3D and 3D Cloud Resolving Model Simulations

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D) have been used to study the response of clouds to large-scale forcing. IN these 3D simulators, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical clouds systems with large horizontal domains at the National Center of Atmospheric Research (NCAR) and at NASA Goddard Space Center. At Goddard, a 3D cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE, GATE, SCSMEX, ARM, and KWAJEX using a 512 by 512 km domain (with 2-km resolution). The result indicate that surface precipitation and latent heating profiles are very similar between the 2D and 3D GCE model simulation. The major objective of this paper are: (1) to assess the performance of the super-parametrization technique, (2) calculate and examine the surface energy (especially radiation) and water budget, and (3) identify the differences and similarities in the organization and entrainment rates of convection between simulated 2D and 3D cloud systems.

Tao, W.-K.↗

Precipitation Processes Developed During ARM (1997), TOGA COARE (1992), GATE (1974), SCSMEX (1998), and KWAJEX (1999): Consistent 2D, Semi-3D and 3D Cloud Resolving Model Simulations

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research (NACAR) and at NASA Goddard Space Flight Center . At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE, SCSMEX and KWAJEX using 512 by 512 km domain (with 2 km resolution). The results indicate that surface precipitation and latent heating profiles are very similar between the 2D and 3D GCE model simulations. The reason for the strong similarity between the 2D and 3D CRM simulations is that the same observed large-scale advective tendencies of potential temperature, water vapor mixing ratio, and horizontal momentum were used as the main focusing in both the 2D and 3D models. Interestingly, the 2D and 3D versions of the CRM used at CSU showed significant differences in the rainfall and cloud statistics for three ARM cases. The major objectives of this paper are: (1) to assess the performance of the super-parameterization technique, (2) calculate and examine the surface energy (especially radiation) and water budgets, and (3) identify the differences and similarities in the organization and entrainment rates of convection between simulated 2D and 3D cloud systems.

Tao, W-K.↗

The proton cyclotron instability and the anisotropy/B inverse correlation

Spacecraft observations in the strongly compressed subsolar magnetosheath show an inverse correlation between the proton temperature anistrophy (T(sub perpendicular to p)/T(sub parallel to p) greater than 1 where perpendicular to and parallel to denote directions perpendicular and parallel to the background magnetic field) and the parallel proton beta (beta(sub parallel to p). This manuscript uses one-dimensional hybrid simulations of the proton cyclotron anisotropy instability in homogeneous electron-proton plasmas to study this correlation which may represent a limited closure relation for fluid theories of anisotropic space plasmas. The emphasis is on driven simulation which increase the temperature anisotropy by periodically reducing the magnetic-field-aligned velocities of the protons. The late-time states from ensembles of both initial value and driven simulations yield very similar expressions for the proton anisotropy/beta(sub parallel to P) inverse correlation, and provide a basis for explaining differences between sheath observations from different spacecraft. The driven simulations also yield expressions for the maximum instability growth rate and the fluctuating field energy as functions of beta (sub parallel to p) and a parameter characterizing the anisotropy driver.

Gary, S. Peter↗

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗

Simulation of pitch angle diffusion of charged particles in a disordered magnetic field

Results are reported for computer simulation experiments in which a statistical ensemble of random magnetic field realizations is generated, orbits of charged particles in the random fields are followed, and a pitch-angle diffusion coefficient is derived from the temporal evolution of the orbits. Diffusion coefficients predicted by three nonlinear theories are compared with the derived coefficients for the standard quasilinear theory of velocity diffusion, and the goals of future simulations are outlined.

Kaiser, T. B.↗

Application of the CloudSat and NEXRAD Radars Toward Improvements in High Resolution Operational Forecasts

As computational power increases, operational forecast models are performing simulations with higher spatial resolution allowing for the transition from sub-grid scale cloud parameterizations to an explicit forecast of cloud characteristics and precipitation through the use of single- or multi-moment bulk water microphysics schemes. investments in space-borne and terrestrial remote sensing have developed the NASA CloudSat Cloud Profiling Radar and the NOAA National Weather Service NEXRAD system, each providing observations related to the bulk properties of clouds and precipitation through measurements of reflectivity. CloudSat and NEXRAD system radars observed light to moderate snowfall in association with a cold-season, midlatitude cyclone traversing the Central United States in February 2007. These systems are responsible for widespread cloud cover and various types of precipitation, are of economic consequence, and pose a challenge to operational forecasters. This event is simulated with the Weather Research and Forecast (WRF) Model, utilizing the NASA Goddard Cumulus Ensemble microphysics scheme. Comparisons are made between WRF-simulated and observed reflectivity available from the CloudSat and NEXRAD systems. The application of CloudSat reflectivity is made possible through the QuickBeam radiative transfer model, with cautious application applied in light of single scattering characteristics and spherical target assumptions. Significant differences are noted within modeled and observed cloud profiles, based upon simulated reflectivity, and modifications to the single-moment scheme are tested through a supplemental WRF forecast that incorporates a temperature dependent snow crystal size distribution.

Molthan, A. L.↗

Risk analysis simulation of rover operations for Mars surface exploration

Risk management advocates have long sought to directly influence the early stages of the systems engineering process through a more effective role in system design trade studies. The principal obstacle to this has been the lack of credible ways to represent and quantify mission risk—that is, a combination of the probability of mission success (“system safety”) and science value—for the project manager and the rest of the design team. If it were possible to quantify mission risk, then the effects of proposed mission and system design changes could be calculated, and along with life-cycle costs, could be used to explore the design space more extensively and select better designs. JPL has been working to build the capability to quantify the probability of mission success using a federation of diverse simulations and models, each of which contributes some vital piece of the puzzle. The initial institutional focus has been on Mars surface operations. This ensemble computing framework enables the diverse models and simulations to work together seamlessly. Recent work at JPL has demonstrated the capability to exercise this ensemble from end-to-end using an Oracle-based database to automatically move results from one model/simulation to the next stage in the analysis.

Shishko, Robert↗

The Impact of Model Configuration and Large-Scale, Upper-Level Forcing on CRM- Simulated Convective Systems

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D, see a brief review by Tao 2003). Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research, at NOAA GFDL, at the U. K. Met. Office, at Colorado State University and at NASA Goddard Space Flight Center (Tao 2003). At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (September 1-7, 1974), SCSMEX (June 2-11, 1998), ARM (June 26-30, 1997) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 by 512 km domain (with 2-km resolution). The results indicate that surface precipitation and latent heating profiles are similar between the 2D and 3D GCE model simulations. However, there are difference in radiation, surface fluxes and precipitation characteristics. The 2D GCE model was used to perform a long-term integration on ARM/GCSS case 4 (22 days at the ARM Southern Great Plains site in March 2000). Preliminary results showed a large temperature bias in the upper troposphere that had not been seen in previous tropical cases. The major objectives of this paper are: (1) to determine the sensitivities to model configuration (i.e., 2D in west-east, south-north or 3D), (2) to identify the differences and similarities in the organization and entrainment rates of convection between 2D- and 3D-simulated ARM cloud systems, and (3) assess the impact of upper tropospheric forcing on tropical and ARM case 4 cases.

Tao, W.-K.↗

The Impact of Model Configuration and Large-Scale, Upper-Level Forcing on CRM-Simulated Convective Systems

Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D, see a brief review by Tao 2003). Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research, at NOAA GFDL, at the U. K. Met. Office, at Colorado State University and at NASA Goddard Space Flight Center (Tao 2003). At Goddard, a 3D Goddard Cumulus Ensemble (GCE) model was used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (September 1-7, 1974), SCSMEX (June 2-11, 1998), ARM (June 26-30, 1997) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 km domain (with 2-kilometer resolution). The results indicate that surface precipitation and latent heating profiles are similar between the 2D and 3D GCE model simulations. However, there are difference in radiation, surface fluxes and precipitation characteristics. The 2D GCE model was used to perform a long-term integration on ARM/GCSS case 4 (22 days at the ARM southern Great Plains site in March 2000). Preliminary results showed a large temperature bias in the upper troposphere that had not been seen in previous tropical cases. The major objectives of this paper are: (1) to determine the sensitivities to model configuration (ie., 2D in west-east, south-north or 3D), (2) to identify the differences and similarities in the organization and entrainment rates of convection between 2D- and 3D-simulated ARM cloud systems, and (3) assess the impact of upper tropospheric forcing on tropical and ARM case 4 cases.

Tao, W.-K.↗