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 289 records · Page 16

Effect of Cross-Correlation on Geomagnetic Forecast Accuracies

Surface geomagnetic observation can determine up to degree L = 14 time-varying spherical harmonic coefficients of the poloidal magnetic field. Assimilation of these coefficients to numerical dynamo simulation could help us understand better the dynamical processes in the Earth's outer core, and to provide more accurate forecast of geomagnetic secular variations (SV). In our previous assimilation studies, only the poloidal magnetic field in the core is corrected by the observations in the analysis. Unobservable core state variables (the toroidal magnetic field and the core velocity field) are corrected via the dynamical equations of the geodynamo. Our assimilation experiments show that the assimilated core state converges near the CMB, implying that the dynamo state is strongly constrained by surface geomagnetic observations, and is pulled closer to the truth by the data. We are now carrying out an ensemble of assimilation runs with 1000 years of geomagnetic and archeo/paleo magnetic record. In these runs the cross correlation between the toroidal and the poloidal magnetic fields is incorporated into the analysis. This correlation is derived from the physical boundary conditions of the toroidal field at the core-mantle boundary (CMB). The assimilation results are then compared with those of the ensemble runs without the cross-correlation, aiming at understanding two fundamental issues: the effect of the crosscorrelation on (1) the convergence of the core state, and (2) the SV prediction accuracies. The constrained dynamo solutions will provide valuable insights on interpreting the observed SV, e.g. the near-equator magnetic flux patches, the core-mantle interactions, and possibly other geodynamic observables.

Kuang, Weijia↗

Application of Aviation Turbulence Information to Air-Traffic Management (ATM): Chapter - 24

Unexpected turbulence especially in the upper troposphere and lower stratosphere where cabin crews and passengers in cruising aircraft are likely to unbuckle causes in-flight injuries, structural damage, and flight delay. Therefore, turbulence information can be used to improve safety while pursuing efficiency in Air-Traffic Management (ATM). In this chapter, simple modeling of aircraft trajectories combined with wind and turbulence predictions can suggest the optimal solution of flight plans that minimizes both total flight time (e.g., fuel consumption) and potential encounters of turbulence from departure to arrival airports. Also, probabilistic ensemble turbulence forecasts are applied to suggest an optimal strategic and tactical ATM route planning in a given weather and turbulence condition in the United States which are evaluated against in situ Eddy Dissipation Rate observations from commercial aircraft. Finally, variations of long-haul trans-Oceanic flight routes and their turbulence potentials are investigated using a global reanalysis data to understand how the upper-level large-scale flow patterns can affect the long-term ATM planning through the changes of winds and turbulence conditions.

Wind Optimal↗

Impact of GPM-retrieved surface meteorology condition on simulations of two winter storms during ICE-POP 2018

Led by the KMA as a component of the WMO's World Weather Research Program (WWRP) Research and Development and Forecast Demonstration Projects (RDP/FDP). Taken place during the Winter Olympics (February-March) of 2018. Goals of ICE-POP: To improve understanding on severe weathers (snowfalls, visibility, rapid wind changes and gusts,) over complex terrain; To improve the predictability of nowcasting and very-short range forecasting with a few kilometer horizontal resolution - Development of NWP-based nowcasting, multi-scale data assimilation and time-lagged ensemble for VSRF, and radar reflectivity and visibility data assimilation; To improve verification for high resolution model considering complex terrain.

Xuanli Li↗

Sub-Seasonal Forecasting of the Stratospheric Wave Events, Sudden Stratospheric Warmings, and Their Influence on the Troposphere

Stratospheric wave events and major sudden stratospheric warming (SSW) events are well captured in high-resolution global forecasts out to 10 days. Tropospheric influences of SSW include statistically significant shifts in the storm tracks and associated surface temperature and precipitation pattern changes. Since these events can in turn influence the troposphere on time scales of 30-60 days the ability to predict these events and the subsequent long-term response on time scales beyond 10 days is of interest. Here we examine the prediction of stratospheric wave events and their evolution using the NASA GMAO (Global Modeling and Assimilation Office) Sub-Seasonal to Seasonal (S2S) system. This is a recently released subseasonal to seasonal forecast system, GEOS-S2S version 2.1. Compared to GMAO's previous system, the new version runs at higher atmospheric resolution (approximately 1/2 degree globally), contains a substantially improved model of the cryosphere, includes additional interactive aerosol model components, and the ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter. Results are based on a comprehensive series of hindcasts starting from the year 2000. They show that while the S2S system is not as accurate at 10 days in forecasting major SSW events as the NASA GMAO Forward Processing system, it can usefully predict stratospheric anomalies out to 20 days and the subsequent stratospheric/tropospheric evolution beyond 30 days.

Coy, L.↗

Comparison of the Performance of Noise Metrics as Predictions of the Annoyance of Stage 2 and Stage 3 Aircraft Overflights

Thirty audiometrically screened test participants judged the relative annoyance of two comparison (variable level) and thirty-four standard (fixed level) signals in an adaptive paired comparison psychoacoustic study. The signal ensemble included both FAR Part 36 Stage 2 and 3 aircraft overflights, as well as synthesized aircraft noise signatures and other non-aircraft signals. All test signals were presented for judgment as heard indoors, in the presence of continuous background noise, under free-field listening conditions in an anechoic chamber. Analyses of the performance of 30 noise metrics as predictors of these annoyance judgments confirmed that the more complex metrics were generally more accurate and precise predictors than the simpler methods. EPNL was somewhat less accurate and precise as a predictor of the annoyance judgments than a duration-adjusted variant of Zwicker's Loudness Level.

Pearsons, Karl S.↗

Progression-regression effects in tracking repeated patterns

Subjects used a position control system to perform compensatory tracking of a repeated input pattern. The input pattern was 20 seconds in duration and was either an arctangent function or the sum of two sine waves. Tracking error decreased with practice and increased with the addition of a concurrent memory task. The shape of the ensemble-average tracking error resembled the shape of the input velocity signal throughout these changes in performance. Regression analyses were used to parameterize these effects and compare these results with the predictions of several conceptualizations of perceptual-motor learning.

Jagacinski, R. J.↗

On the use of the noncentral chi-square density function for the distribution of helicopter spectral estimates

A probability density function for the variability of ensemble averaged spectral estimates from helicopter acoustic signals in Gaussian background noise was evaluated. Numerical methods for calculating the density function and for determining confidence limits were explored. Density functions were predicted for both synthesized and experimental data and compared with observed spectral estimate variability.

Garber, Donald P.↗

Disturbance Dynamics in Transitional and Turbulent Flows

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. In laminar flows, such disturbances are linear superpositions of modes governed by the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport equations for representative mean quantities. The link between a description based on a deterministic evolution equation and a probability based mean transport equation is established. Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, a probability distribution must be defined even in the laminar case. Using this probability distribution, it is shown that the exponential decay of the linear modes in the laminar regime can be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and the dissipation rate. The evolution of these mean disturbance quantities is then described by transport equations similar to those for the corresponding turbulent decaying flow. Second, homogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), is studied. The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined.

Grosch, Chester E.↗

Impacts of Marine Surface Pressure Observations From A Spaceborne Differential Absorption Radar Investigated With an Observing System Simulation Experiment

A new instrument has been proposed for measuring surface air pressure over the marine surface with a combined active/passive scanning multi-channel differential absorption radar (DAR) to provide an estimate of the total atmospheric column oxygen content. A demonstrator instrument, the Microwave Barometric Radar and Sounder (MBARS), has been funded by the National Aeronautics and Space Administration (NASA) for airborne test missions. Here, a proof-of-concept study to evaluate the potential impact of spaceborne surface pressure data on numerical weather prediction is performed using the Goddard Modeling and Assimilation Office global observing system simulation experiment (OSSE) framework. This OSSE framework employs the Goddard Earth Observing System model and the hybrid 4D ensemble variational Gridpoint Statistical Interpolation data assimilation system. Multiple flight and scanning configurations of potential spaceborne orbits are examined. Swath width and observation spacing for the surface pressure data are varied to explore a range of sampling strategies. For wider swaths, the addition of surface pressures reduces the root mean square surface pressure analysis error by as much as 20% over some ocean regions. The forecast sensitivity observation impact tool estimates impacts on the Pacific Ocean basin boundary layer 24-hour forecast temperatures for spaceborne surface pressures on par with rawinsondes and aircraft, and greater impacts than the current network of ships and buoys. The largest forecast impacts are found in the southern hemisphere extratropics.

N.C. Privé↗

Using Results from the PlioMIP Ensemble to Investigate the Greenland Ice Sheet During the Mid-Pliocene Warm Period

During an interval of the Late Pliocene, referred to here as the mid-Pliocene Warm Period (mPWP; 3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the end of this century, and atmospheric carbon dioxide concentrations were higher than preindustrial levels. Sea level was also higher than today, implying a significant reduction in the extent of the ice sheets. Thus, the mPWP provides a natural laboratory in which to investigate the long-term response of the Earth's ice sheets and sea level in a warmer-than-present-day world. At present, our understanding of the Greenland ice sheet during the mPWP is generally based upon predictions using single climate and ice sheet models. Therefore, it is essential that the model dependency of these results is assessed. The Pliocene Model Intercomparison Project (PlioMIP) has brought together nine international modeling groups to simulate the warm climate of the Pliocene. Here we use the climatological fields derived from the results of the 15 PlioMIP climate models to force an offline ice sheet model. We show that mPWP ice sheet reconstructions are highly dependent upon the forcing climatology used, with Greenland reconstructions ranging from an ice-free state to a near modern ice sheet. An analysis of the surface albedo variability between the climate models over Greenland offers insights into the drivers of inter-model differences. As we demonstrate that the climate model dependency of our results is high, we highlight the necessity of data-based constraints of ice extent in developing our understanding of the mPWP Greenland ice sheet.

sea level↗

Comparing Physics Scheme Performance for a Lake Effect Snowfall Event in Northern Lower Michigan

High resolution forecast models, such as those used to predict severe convective storms, can also be applied to predictions of lake effect snowfall. A high resolution WRF model forecast model is provided to support operations at NWS WFO Gaylord, Michigan, using a 12 ]km and 4 ]km nested configuration. This is comparable to the simulations performed by other NWS WFOs adjacent to the Great Lakes, including offices in the NWS Eastern Region who participate in regional ensemble efforts. Ensemble efforts require diversity in initial conditions and physics configurations to emulate the plausible range of events in order to ascertain the likelihood of different forecast scenarios. In addition to providing probabilistic guidance, individual members can be evaluated to determine whether they appear to be biased in some way, or to better understand how certain physics configurations may impact the resulting forecast. On January 20 ]21, 2011, a lake effect snow event occurred in Northern Lower Michigan, with cooperative observing and CoCoRaHS stations reporting new snow accumulations between 2 and 8 inches and liquid equivalents of 0.1 ]0.25 h. The event of January 21, 2011 was particularly well observed, with numerous surface reports available. It was also well represented by the WRF configuration operated at NWS Gaylord. Given that the default configuration produced a reasonable prediction, it is used here to evaluate the impacts of other physics configurations on the resulting prediction of the primary lake effect band and resulting QPF. Emphasis here is on differences in planetary boundary layer and cloud microphysics parameterizations, given their likely role in determining the evolution of shallow convection and precipitation processes. Results from an ensemble of seven microphysics schemes and three planetary boundary layer schemes are presented to demonstrate variability in forecast evolution, with results used in an attempt to improve the forecasts in the 2011 ]2012 lake effect season.

Molthan, Andrew↗

Simultaneous Estimation of Model State Variables and Observation and Forecast Biases Using a Two-Stage Hybrid Kalman Filter

In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Data Assimilation↗

Radiance and Jacobian Intercomparison of Radiative Transfer Models Applied to HIRS and AMSU Channels

The goals of this study are the evaluation of current fast radiative transfer models (RTMs) and line-by-line (LBL) models. The intercomparison focuses on the modeling of 11 representative sounding channels routinely used at numerical weather prediction centers: 7 HIRS (High-resolution Infrared Sounder) and 4 AMSU (Advanced Microwave Sounding Unit) channels. Interest in this topic was evidenced by the participation of 24 scientists from 16 institutions. An ensemble of 42 diverse atmospheres was used and results compiled for 19 infrared models and 10 microwave models, including several LBL RTMs. For the first time, not only radiances, but also Jacobians (of temperature, water vapor and ozone) were compared to various LBL models for many channels. In the infrared, LBL models typically agree to within 0.05-0.15 K (standard deviation) in terms of top-of-the-atmosphere brightness temperature (BT). Individual differences up to 0.5 K still exist, systematic in some channels, and linked to the type of atmosphere in others. The best fast models emulate LBL BTs to within 0.25 K, but no model achieves this desirable level of success for all channels. The ozone modeling is particularly challenging, In the microwave, fast models generally do quite well against the LBL model to which they were tuned. However significant differences were noted among LBL models, Extending the intercomparison to the Jacobians proved very useful in detecting subtle and more obvious modeling errors. In addition, total and single gas optical depths were calculated, which provided additional insight on the nature of differences. Recommendations for future intercomparisons are suggested.

Garand, L.↗

Radiance and Jacobian Intercomparison of Radiative Transfer Models Applied to HIRS and AMSU Channels

The goals of this study are the evaluation of current fast radiative transfer models (RTMs) and line-by-line (LBL) models. The intercomparison focuses on the modeling of 11 representative sounding channels routinely used at numerical weather prediction centers: seven HIRS (High-resolution Infrared Sounder) and four AMSU (Advanced Microwave Sounding Unit) channels. Interest in this topic was evidenced by the participation of 24 scientists from 16 institutions. An ensemble of 42 diverse atmospheres was used and results compiled for 19 infrared models and 10 microwave models, including several LBL RTMs. For the first time, not only radiances, but also Jacobians (of temperature, water vapor, and ozone) were compared to various LBL models for many channels. In the infrared, LBL models typically agree to within 0.05-0.15 K (standard deviation) in terms of top-of-the-atmosphere brightness temperature (BT). Individual differences up to 0.5 K still exist, systematic in some channels, and linked to the type of atmosphere in others. The best fast models emulate LBL BTs to within 0.25 K, but no model achieves this desirable level of success for all channels. The ozone modeling is particularly challenging. In the microwave, fast models generally do quite well against the LBL model to which they were tuned. However significant differences were noted among LBL models. Extending the intercomparison to the Jacobians proved very useful in detecting subtle and more obvious modeling errors. In addition, total and single gas optical depths were calculated, which provided additional insight on the nature of differences. Recommendations for future intercomparisons are suggested.

Garand, L.↗

GEOS S2S-2_1: GMAO's New High Resolution Seasonal Prediction System

A new version of the modeling and analysis system used to produce sub-seasonal to seasonal forecasts has just been released by the NASA Goddard Global Modeling and Assimilation Office. The new version runs at higher atmospheric resolution (approximately 12 degree globally), contains a substantially improved model description of the cryosphere, and includes additional interactive earth system model components (aerosol model). In addition, the Ocean data assimilation system has been replaced with a Local Ensemble Transform Kalman Filter. Here will describe the new system, along with the plans for the future (GEOS S2S-3_0) which will include a higher resolution ocean model and more interactive earth system model components (interactive vegetation, biomass burning from fires). We will also present results from a free-running coupled simulation with the new system and results from a series of retrospective seasonal forecasts. Results from retrospective forecasts show significant improvements in surface temperatures over much of the northern hemisphere and a much improved prediction of sea ice extent in both hemispheres. The precipitation forecast skill is comparable to previous S2S systems, and the only trade off is an increased double ITCZ, which is expected as we go to higher atmospheric resolution.

Molod, Andrea↗

Scaling Observation Error for Optimal Assimilation of CCI SST Data into a Regional HYCOM EnOI System

South Africa currently possesses no operational ocean forecasting system for the purpose of predicting ocean state variables including temperature,salinity and velocity. Substantial initial efforts towards this goal have been made and resulted in a system using a regional Hybrid Coordinate Ocean Model (HYCOM) along with the Ensemble Optimal Interpolation (EnOI)assimilation scheme. Assimilating only sea surface temperature (SST) observations from the Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) product into the system resulted in a degraded forecast. Aiming to address this, Climate Change Initiative (CCI) SSTs are assimilated into the system in an effort to improve the forecast skill. Observation errors in the assimilated product are used in the EnOI to determine whether more confidence should be placed in the model or observations in producing the analysis, but overconfidence in observations can shock the model and result in failure. To tweak the impact of the assimilation, a scaling factor is applied in the assimilation code. A scaling factor of 25 was found to produce a favourable result with lowest mean root mean square error (RMSE;1.098°C) between the model and observations over time. Postulating the error to be overconfident, a floor value is introduced in order to set a minimum value for the observation error thereby reducing confidence in the observations. These experiments fared less favourably with a floor value of 0.5 and a scaling factor of 15 producing the best mean RMSE (1.118°C).

Luyt, Hermann↗

Physical Models of Layered Polar Firn Brightness Temperatures from 0.5 to 2 GHz

We investigate physical effects influencing 0.5-2 GHz brightness temperatures of layered polar firn to support the Ultra Wide Band Software Defined Radiometer (UWBRAD) experiment to be conducted in Greenland and in Antarctica. We find that because ice particle grain sizes are very small compared to the 0.5-2 GHz wavelengths, volume scattering effects are small. Variations in firn density over cm- to m-length scales, however, cause significant effects. Both incoherent and coherent models are used to examine these effects. Incoherent models include a 'cloud model' that neglects any reflections internal to the ice sheet, and the DMRT-ML and MEMLS radiative transfer codes that are publicly available. The coherent model is based on the layered medium implementation of the fluctuation dissipation theorem for thermal microwave radiation from a medium having a nonuniform temperature. Density profiles are modeled using a stochastic approach, and model predictions are averaged over a large number of realizations to take into account an averaging over the radiometer footprint. Density profiles are described by combining a smooth average density profile with a spatially correlated random process to model density fluctuations. It is shown that coherent model results after ensemble averaging depend on the correlation lengths of the vertical density fluctuations. If the correlation length is moderate or long compared with the wavelength (approximately 0.6x longer or greater for Gaussian correlation function without regard for layer thinning due to compaction), coherent and incoherent model results are similar (within approximately 1 K). However, when the correlation length is short compared to the wavelength, coherent model results are significantly different from the incoherent model by several tens of kelvins. For a 10-cm correlation length, the differences are significant between 0.5 and 1.1 GHz, and less for 1.1-2 GHz. Model results are shown to be able to match the v-pol SMOS data closely and predict the h-pol data for small observation angles.

layered polar firn emission↗

On the application of a hairpin vortex model of wall turbulence to trailing edge noise prediction

The goal is to develop a technique via a hairpin vortex model of the turbulent boundary layer, which would lead to the estimation of the aerodynamic input for use in trailing edge noise prediction theories. The work described represents an initial step in reaching this goal. The hairpin vortex is considered as the underlying structure of the wall turbulence and the turbulent boundary layer is viewed as an ensemble of typical hairpin vortices of different sizes. A synthesis technique is examined which links the mean flow and various turbulence quantities via these typical vortices. The distribution of turbulence quantities among vortices of different scales follows directly from the probability distribution needed to give the measured mean flow vorticity. The main features of individual representative hairpin vortices are discussed in detail and a preliminary assessment of the synthesis approach is made.

Liu, N. S.↗