Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble prediction”

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 181 records · Page 10

Operational Dust Prediction

Over the last few years, numerical prediction of dust aerosol concentration has become prominent at several research and operational weather centres due to growing interest from diverse stakeholders, such as solar energy plant managers, health professionals, aviation and military authorities and policymakers. Dust prediction in numerical weather prediction-type models faces a number of challenges owing to the complexity of the system. At the centre of the problem is the vast range of scales required to fully account for all of the physical processes related to dust. Another limiting factor is the paucity of suitable dust observations available for model, evaluation and assimilation. This chapter discusses in detail numerical prediction of dust with examples from systems that are currently providing dust forecasts in near real-time or are part of international efforts to establish daily provision of dust forecasts based on multi-model ensembles. The various models are introduced and described along with an overview on the importance of dust prediction activities and a historical perspective. Assimilation and evaluation aspects in dust prediction are also discussed.

assimilation↗

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.↗

Mass Conservation and Positivity Preservation with Ensemble-type Kalman Filter Algorithms

Maintaining conservative physical laws numerically has long been recognized as being important in the development of numerical weather prediction (NWP) models. In the broader context of data assimilation, concerted efforts to maintain conservation laws numerically and to understand the significance of doing so have begun only recently. In order to enforce physically based conservation laws of total mass and positivity in the ensemble Kalman filter, we incorporate constraints to ensure that the filter ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. We show that the analysis steps of ensemble transform Kalman filter (ETKF) algorithm and ensemble Kalman filter algorithm (EnKF) can conserve the mass integral, but do not preserve positivity. Further, if localization is applied or if negative values are simply set to zero, then the total mass is not conserved either. In order to ensure mass conservation, a projection matrix that corrects for localization effects is constructed. In order to maintain both mass conservation and positivity preservation through the analysis step, we construct a data assimilation algorithms based on quadratic programming and ensemble Kalman filtering. Mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate constraints. Some simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. The results show clear improvements in both analyses and forecasts, particularly in the presence of localized features. Behavior of the algorithm is also tested in presence of model error.

EnKR↗

Numerical and experimental investigation of nonswirling and swirling confined jets

An investigation of the influence of large scale structures on the flow development for coaxial jets with sudden expansion (with and without swirl) is presented. Both an experimental study and numerical predictions were performed for a configuration corresponding to that considered by Johnson and Bennett and Roback and Johnson. The effects of large scale structures on the swirling and nonswirling flows are documented, in particular their influence on turbulence modeling and the numerical simulation. The ensemble-averaged, time-dependent Navier-Stokes equations are solved by an LBI procedure to predict the turbulent flow field. Effects of artificial dissipation and placement of the upstream boundary in the numerical computation are also discussed.

Brondum, D. C.↗

Two dimensional turbulence in inviscid fluids or guiding center plasmas

Analytic theory for two-dimensional turbulent equilibria for the inviscid Navier-Stokes equations is examined mathematically. Application of the technique to electrostatic guiding center plasma is discussed. A good fit is demonstrated for the approach to a predicted energy per Fourier mode obtained from a two-temperature canonical ensemble. Negative as well as positive temperature regimes are explored. Fluctuations about the mean energy per mode also compare well with theory. In the regime of alpha less than zero, beta greater than zero, with the minimum value of alpha plus beta times k squared near zero, contour plots of the stream function reveal macroscopic vortex structures similar to those seen previously in discrete vortex simulations. Eulerian direct interaction equations, which can be used to follow the approach to inviscid equilibrium, are derived.

Seyler, C. E., Jr.↗

MHD solution of interplanetary disturbances generated by simulated velocity perturbations

An MHD time-dependent numerical simulation, restricted to the solar equatorial plane, is used to demonstrate the interplanetary disturbances caused by several simplified coronal holes. Each 'hole' is assumed to have a configuration such that the higher solar wind velocity produced within their longitudinal extent is Gaussian over a 7-day period at the inner boundary (0.3 AU) of the numerical simulation. A second, twin coronal hole is assumed to rotate on the solar disk behind its predecessor. It is shown that the first coronal hole-produced interplanetary shock ensemble is overtaken by the second ensemble because of the higher velocity, lower density environment into which the latter propagates. A number of features predicted by MHD similarity theory are confirmed by the numerical simulation. These features include (1) strong azimuthal magnetic and plasma density compression, accompanied with average temperature depression, at the contact surface between forward and reverse shock ensembles, and (2) increasing spatial separation distance between forward and reverse shocks.

Dryer, M.↗

Stationarity of magnetohydrodynamic fluctuations in the solar wind

Solar wind research and studies of charged particle propagation often assume that the interplanetary magnetic field represents a stationary random process. The extent to which ensemble averages of the solar wind magnetic fields follow the asymptotic behavior predicted by the ergodic theorem was investigated. Several time periods, including a span of nearly two years, are analyzed. Data intervals which span many solar rotations satisfy the conditions of weak stationarity if the effects of solar rotation are included in the asymptotic analysis. Shorter intervals which include a small integral number of interplanetary sectors also satisfy weak stationarity. The results are illustrated using magnetometer data from the ISEE-3, Voyager and IMP spacecraft.

Matthaeus, W. H.↗

On the Causes of the 1930s Dust Bowl

During the 1930s the United States experienced one of the most devastating droughts of the last century. The drought affected almost 2/3 of the country and parts of Mexico and Canada and was infamous for the numerous dust storms that lead to the characterization of much of the Great Plains as the "Dust Bowl". Results from an ensemble of 100-year simulations with the NASA Seasonal- to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed SSTs show that the model reproduces the basic features of the 1930s drought. In t h s study we exploit this realism in the model simulation to examine in more detail the role of the SST and soil moisture in the development and maintenance the 1930s drought. Results will be presented from simulations in which the SST anomalies are confined to the separate ocean basins, as well as from runs in which soil moisture feedback is turn off.

Schubert, Siegfried↗

On the Causes of the 1930s Dust Bowl

During the 1930s the United States experienced one of the most devastating droughts of the last century. The drought affected almost 2/3 of the country and parts of Mexico and Canada and was infamous for the numerous dust storms that lead to the characterization of much of the Great Plains as the 'Dust Bowl'. Results from an ensemble of 100-year simulations with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed SSTs show that the model reproduces the basic features of the 1930s drought. In this study we exploit this realism in the model simulation to examine in more detail the role of the SST and soil moisture in the development and maintenance the 1930s drought. Results will be presented from simulations in which the SST anomalies are confined to the separate ocean basins, as well as from runs in which soil moisture feedback is turn off.

Schubert, Siegfried↗

DYAMOND: The DYnamics of the Atmospheric General Circulation Modeled on Non-Hydrostatic Domains

A review of the experimental protocol and motivation for DYAMOND, the first intercomparison project of global storm-resolving models, is presented. Nine models submitted simulation output for a 40-day (1 August–10 September 2016) intercomparison period. Eight of these employed a tiling of the sphere that was uniformly less than 5 km. By resolving the transient dynamics of convective storms in the tropics, global storm-resolving models remove the need to parameterize tropical deep convection, providing a fundamentally more sound representation of the climate system and a more natural link to commensurately high-resolution data from satellite-borne sensors. The models and some basic characteristics of their output are described in more detail, as is the availability and planned use of this output for future scientific study. Tropically and zonally averaged energy budgets, precipitable water distributions, and precipitation from the model ensemble are evaluated, as is their representation of tropical cyclones and the predictability of column water vapor, the latter being important for tropical weather.

Bjorn Stevens↗

Upper Limits of Predictability in Long-Range Climate/Hydrologic Forecasts

The accurate forecasting of el nino or la nina conditions in the tropical Pacific can potentially lead to valuable predictions of hydrological anomalies over land at seasonal to interannual timescales. Even with highly accurate earth system models, though, our ability to generate these continental forecasts will always be limited by the chaotic nature of the atmospheric circulation. The nature of this fundamental limitation is explored through the use of 16-member ensembles of multi-decade GCM simulations. In each simulation of the first ensemble, sea surface temperatures (SSTs) are given the same realistic interannual variations over a 45-year period, and land surface state is allowed to evolve with that of the atmosphere. Analysis of the results shows that the SSTs control the temporal organization of continental precipitation anomalies to a significant extent in the tropics and to a much smaller extent in midlatitudes. In each simulation of the second ensemble, we prescribe SSTs as before, but we also prescribe interannual variations in the low frequency component of evaporation efficiency over land. Thus, in the second ensemble, we effectively make the extreme assumption that surface boundary conditions across the globe are perfectly predictable, and we quantify the consistency with which the atmosphere (particularly precipitation) responds to these boundary conditions. The resulting "absolute upper limit" on the predictability of precipitation is found to be quite high in the tropics yet only moderate in many midlatitude regions.

Koster, R. D.↗

The Experimental Regional Ensemble Forecast System (ExREF): Its Use in NWS Forecast Operations and Preliminary Verification

The Experimental Regional Ensemble Forecast (ExREF) system is a tool for the development and testing of new Numerical Weather Prediction (NWP) methodologies. ExREF is run in near‐realtime by the Global Systems Division (GSD) of the NOAA Earth System Research Laboratory (ESRL) and its products are made available through a website, an ftp site, and via the Unidata Local Data Manager (LDM). The ExREF domain covers most of North America and has 9‐km horizontal grid spacing. The ensemble has eight members, all employing WRF‐ARW. The ensemble uses a variety of initial conditions from LAPS and the Global Forecasting System (GFS) and multiple boundary conditions from the GFS ensemble. Additionally, a diversity of physical parameterizations is used to increase ensemble spread and to account for the uncertainty in forecasting extreme precipitation events. ExREF has been a component of the Hydrometeorology Testbed (HMT) NWP suite in the 2012‐2013 and 2013‐2014 winters. A smaller domain covering just the West Coast was created to minimize band‐width consumption for the NWS. This smaller domain has and is being distributed to the National Weather Service (NWS) Weather Forecast Office and California Nevada River Forecast Center in Sacramento, California, where it is ingested into the Advanced Weather Interactive Processing System (AWIPS I and II) to provide guidance on the forecasting of extreme precipitation events. This paper will review the cooperative effort employed by NOAA ESRL, NASA SPoRT (Short‐term Prediction Research and Transition Center), and the NWS to facilitate the ingest and display of ExREF data utilizing the AWIPS I and II D2D and GFE (Graphical Software Editor) software. Within GFE is a very useful verification software package called BoiVer that allows the NWS to utilize the River Forecast Center's 4 km gridded QPE to compare with all operational NWP models 6‐hr QPF along with the ExREF mean 6‐hr QPF so the forecasters can build confidence in the use of the ExREF in preparing their rainfall forecasts. Preliminary results will be presented.

Reynolds, David↗

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗

An Agile Beam Transmit Array Using Coupled Oscillator Phase Control

A few years ago York and colleagues suggested that injection locking of voltage controlled oscillators could be used to implement beam steering in a phased array [I]. The scheme makes use of the fact that when an oscillator is injection locked to an external signal, the phase difference between the output of the oscillator and the injection signal is governed by the difference between the injection frequency and the free running frequency of the oscillator (the frequency to which the oscillator is tuned). Thus, if voltage controlled oscillators (VCOs) are used, this phase difference is controlled by an applied voltage. Now, if a set of such oscillators are coupled to nearest neighbors, they can be made to mutually injection lock and oscillate as an ensemble. If they are all tuned to the same frequency, they will all oscillate in phase. Thus, if the outputs are connected to radiating elements forming a linear array, the antenna will radiate normal to the line of elements. Scanning is accomplished by antisymmetrically detuning the end oscillators in the array by application of a pair of appropriate voltages to their tuning ports. This results in a linear phase progression across the array which is just the phasing required to scan the beam. The scan angle is determined by the degree of detuning. We have constructed a seven element one dimensional agile beam array at S-band based on the above principle. Although, a few such arrays have been built in the past, this array possesses two unique features. First, the VCO MMICs have buffer amplifiers which isolate the output from the tuning circuit, and second, the oscillators are weakly coupled to each other at their resonant circuits rather than their outputs. This results in a convenient isolation between the oscillator array design and the radiating aperture design. An important parameter in the design is the so called coupling phase which determines the phase shift of the signals passing from one oscillator to its neighbors. Using this array, we have been able to verify the theoretical predictions concerning the effect of this phase on both the locking range and ensemble frequency of the array. However, the scan range achieved fell somewhat short of the theoretical value because of the amplitude variation of the oscillator outputs with tuning.

Pogorzelski, Ronald S.↗

The Subseasonal Experiment (SubX): A Multi-Model Subseasonal Prediction Experiment

SubX is a multi-model subseasonal prediction experiment designed around operational requirements with the goal of improving subseasonal forecasts. Seven global models have produced seventeen years of retrospective (re-) forecasts and more than a year of weekly real-time forecasts. The re-forecasts and forecasts are archived at the Data Library of the International Research Institute for Climate and Society, Columbia University, providing a comprehensive database for research on subseasonal to seasonal predictability and predictions. The SubX models show skill for temperature and precipitation three weeks ahead of time in specific regions. The SubX multi-model ensemble mean is more skillful than any individual model overall. Skill in simulating the Madden-Julian Oscillation (MJO) and the North Atlantic Oscillation (NAO), two sources of subseasonal predictability, is also evaluated with skillful predictions of the MJO four weeks in advance and of the NAO 2 weeks in advance. SubX is also able to make useful contributions to operational forecast guidance at the Climate Prediction Center. Additionally, SubX provides information on the potential for extreme precipitation associated with tropical cyclones which can help emergency management and aid organizations to plan for disasters. (Capsule Summary) A research to operations project in service of developing better operational subseasonal forecasts.

Precipitation↗

A comparison of CO(J = 1 yields 0) and CO(J = 2 yields 1) emission in the Milky Way molecular ring

We have carried out a CO (J = 2 right arrow 1) survey of the Scutum Arm region of the Milky Way molecular ring. Our goals are to compare CO (J = 2 right arrow 1) maps of individual Galactic clouds with the large-scale CO (J = 2 right arrow 1) emission from the Galactic plane, and to predict the CO (J = 2 right arrow 1) appearance of a Galactic cloud ensemble in an external galaxy. The angular resolution and spatial coverage of our survey are compatible with the existing CO (J = 1 right arrow 0) survey of this region by Sanders et al., which we use for comparison. We identify 34 molecular clouds in our map region; their relationships between size and line width and between virial mass and luminosity are consistent with the relationships seen in CO (J = 1 right arrow 0) emission. However, we note that previous studies have shown considerable variation in these relationships; we attribute much of this variation to differences in their cloud definition algorithm and statistical method. We find that the median ratio of integrated emission in the two lowest CO transitions for the clouds in our sample is I(2 right arrow 1)/ I(1 right arrow 0) = 0.69, implying that the typical emitting region in the line of sight contains cold gas that is not actively star-forming. Our conclusion that the molecular ring emission is not dominated by star-forming regions is consistent with other large-scale studies of the Milky Way. Our distributions of cloud size and temperature also imply that relatively massive molecular clouds that lac k star formation, such as Maddalena's cloud, are not rare in the inner Galaxy.

Chiar, Jean E.↗

Meridional Propagation of the MJO/ISO and Asian Monsoon Variability

In this study we examine the links between tropical heating, the Madden Julian Oscillation (MJO)/Intraseasonal Oscillation (ISO), and the Asian monsoon. We are particularly interested in isolating the nature of the poleward propagation of the ISO/MJO in the monsoon region. We examine both observations and idealized "MJO heating" experiments employing the NASA Seasonal-Interannual Prediction Project (NSIPP) atmospheric general circulation model (AGCM). In the idealized 10-member ensemble simulations, the model is forced by climatological SST and an idealized eastward propagating heating profile that is meant to mimic the canonical heating associated with the MJO in the Indian Ocean and western Pacific. In order to understand the impact of SST on the off equatorial convection (or Rossby-wave response), a second set of 10-member ensemble simulations is carried out with the climatological SSTs shifted in time by 6-months. The observational analysis highlights the strong link between the Indian summer monsoon and the tropical ISO/MJO activity and heating. This includes the well-known meridional propagation that affects the summer monsoons of both hemispheres. The AGCM experiments with the idealized eastward propagating MJO-like heating reproduce the observed meridional propagation including the observed seasonal differences. The impact of the SSTs are to enhance the magnitude of the propagation into the summer hemispheres. The results suggest that the winter/summer differences associated with the MJO/ISO are auxiliary features that depend on the MJO's environment (basic state and boundary conditions) and are not the result of fundamental differences in the MJO itself.

Wu, Man Li C.↗

A Single-column Model Ensemble Approach Applied to the TWP-ICE Experiment

Single-column models (SCM) are useful test beds for investigating the parameterization schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best estimate large-scale observations prescribed. Errors estimating the observations will result in uncertainty in modeled simulations. One method to address the modeled uncertainty is to simulate an ensemble where the ensemble members span observational uncertainty. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best estimate product. These data are then used to carry out simulations with 11 SCM and two cloud-resolving models (CRM). Best estimate simulations are also performed. All models show that moisture-related variables are close to observations and there are limited differences between the best estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the surface evaporation term of the moisture budget between the SCM and CRM. Differences are also apparent between the models in the ensemble mean vertical structure of cloud variables, while for each model, cloud properties are relatively insensitive to forcing. The ensemble is further used to investigate cloud variables and precipitation and identifies differences between CRM and SCM particularly for relationships involving ice. This study highlights the additional analysis that can be performed using ensemble simulations and hence enables a more complete model investigation compared to using the more traditional single best estimate simulation only.

numerical weather forecasting↗