Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical forecasting”

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 379 records · Page 21

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO Seasonal Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite sea surface salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (North American Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF (Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive satellite) V4.1 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by satellite at approximately 1-centimeter depth) and in situ measurements (typically measured by Argo floats at 3-meters depth). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from GMAO S2S spring reanalyses that assimilate/withhold along-track (L2) SSS. In particular, we contrast forecasts of the big 2015 El Nino, the 2017 La Nina and the 2018 weak El Nino. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrade the mixed layer depth leading to more accurate coupled air/sea interaction. From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of gridded satellite sea surface salinity observations on dynamical ENSO forecasts for the big 2015 El Nino.

Hackert, E.↗

The WRF-Solar Ensemble Prediction System to Provide Solar Irradiance Probabilistic Forecasts

In this study, we introduce the recently developed WRF-solar ensemble prediction system and a calibration method. The performances of forecast models are evaluated using the National Solar Radiation Database observational analysis for day-ahead solar irradiance predictions. The results demonstrate that the ensemble forecast improves the quality of the forecasts by considering the uncertainty of each ensemble member. The analog ensemble calibration contributed to the reduction of positive bias and an overall improvement in the probabilistic attributes, such as reliability and statistical consistency.

14 SOLAR ENERGY↗

A comparison of limited-area energetic processes between observations and primitive equation model predictions

Energetic analyses of the NMC initial conditions and NMC six-layer primitive equation operational prediction model 12-hr forecast for a developing cyclone are presented. Consideration is given to the total kinetic energy, the energetics of the divergent and nondivergent flows and the baroclinic (vertical shear flow) and barotropic (vertical mean flow) components of the kinetic energy. It is found that the model initial conditions lose 10-15% of the kinetic energy at various levels compared to a limited-area multivariate statistical analysis of the observational data, leading to a decrease in the horizontal kinetic energy flux, a misrepresentation of the synoptic scale wave system in the 12-hr forecast. Similar results are obtained for the nondivergent flow, while the divergent flow energetics are not reproduced accurately by the model. The horizontal flux terms of the vertical mean and vertical shear energetics are also not found to be reproduced in the upper levels, although horizontal flux contributions to the baroclinic component are improved at middle and lower levels. Finally, vertical shear kinetic energy generation is found to be well represented in the model prediction, however kinetic energy conversion between vertical shear and mean flow is not reproduced in the lower layer.

Alpert, J. C.↗

Evaluation of the Regional Atmospheric Modeling System in the Eastern Range Dispersion Assessment System

The Applied Meteorology Unit is conducting an evaluation of the Regional Atmospheric Modeling System (RAMS) contained within the Eastern Range Dispersion Assessment System (ERDAS). ERDAS provides emergency response guidance for operations at the Cape Canaveral Air Force Station and the Kennedy Space Center in the event of an accidental hazardous material release or aborted vehicle launch. The prognostic data from RAMS is available to ERDAS for display and is used to initialize the 45th Range Safety (45 SW/SE) dispersion model. Thus, the accuracy of the 45 SW/SE dispersion model is dependent upon the accuracy of RAMS forecasts. The RAMS evaluation task consists of an objective and subjective component for the Florida warm and cool seasons of 1999-2000. The objective evaluation includes gridded and point error statistics at surface and upper-level observational sites, a comparison of the model errors to a coarser grid configuration of RAMS, and a benchmark of RAMS against the widely accepted Eta model. The warm-season subjective evaluation involves a verification of the onset and movement of the Florida east coast sea breeze and RAMS forecast precipitation. This interim report provides a summary of the RAMS objective and subjective evaluation for the 1999 Florida warm season only.

Case, Jonathan↗

Spatial and Temporal scales of time-averaged 700 MB height anomalies

The monthly and seasonal forecasting technique is based to a large extent on the extrapolation of trends in the positions of the centers of time averaged geopotential height anomalies. The complete forecasted height pattern is subsequently drawn around the forecasted anomaly centers. The efficacy of this technique was tested and time series of observed monthly mean and 5 day mean 700 mb geopotential heights were examined. Autocorrelation statistics are generated to document the tendency for persistence of anomalies. These statistics are compared to a red noise hypothesis to check for evidence of possible preferred time scales of persistence. Space-time spectral analyses at middle latitudes are checked for evidence of periodicities which could be associated with predictable month-to-month trends. A local measure of the average spatial scale of anomalies is devised for guidance in the completion of the anomaly pattern around the forecasted centers.

Gutzler, D.↗

Assessing the Impact of Pre-gpm Microwave Precipitation Observations in the Goddard WRF Ensemble Data Assimilation System

The forthcoming Global Precipitation Measurement (GPM) Mission will provide next generation precipitation observations from a constellation of satellites. Since precipitation by nature has large variability and low predictability at cloud-resolving scales, the impact of precipitation data on the skills of mesoscale numerical weather prediction (NWP) is largely affected by the characterization of background and observation errors and the representation of nonlinear cloud/precipitation physics in an NWP data assimilation system. We present a data impact study on the assimilation of precipitation-affected microwave (MW) radiances from a pre-GPM satellite constellation using the Goddard WRF Ensemble Data Assimilation System (Goddard WRF-EDAS). A series of assimilation experiments are carried out in a Weather Research Forecast (WRF) model domain of 9 km resolution in western Europe. Sensitivities to observation error specifications, background error covariance estimated from ensemble forecasts with different ensemble sizes, and MW channel selections are examined through single-observation assimilation experiments. An empirical bias correction for precipitation-affected MW radiances is developed based on the statistics of radiance innovations in rainy areas. The data impact is assessed by full data assimilation cycling experiments for a storm event that occurred in France in September 2010. Results show that the assimilation of MW precipitation observations from a satellite constellation mimicking GPM has a positive impact on the accumulated rain forecasts verified with surface radar rain estimates. The case-study on a convective storm also reveals that the accuracy of ensemble-based background error covariance is limited by sampling errors and model errors such as precipitation displacement and unresolved convective scale instability.

microwave observations↗

The MSFC Solar Activity Future Estimation (MSAFE) Model

The MSAFE model provides forecasts for the solar indices SSN, F10.7, and Ap. These solar indices are used as inputs to many space environment models used in orbital spacecraft operations and space mission analysis. Forecasts from the MSAFE model are provided on the MSFC Natural Environments Branch's solar webpage and are updated as new monthly observations come available. The MSAFE prediction routine employs a statistical technique that calculates deviations of past solar cycles from the mean cycle and performs a regression analysis to predict the deviation from the mean cycle of the solar index at the next future time interval. The prediction algorithm is applied recursively to produce monthly smoothed solar index values for the remaining of the cycle. The forecasts are initiated for a given cycle after about 8 to 12 months of observations are collected. A forecast made at the beginning of cycle 24 using the MSAFE program captured the cycle fairly well with some difficulty in discerning the double peak that occurred at solar cycle maximum.

Suggs, Ronnie J.↗

Applied Meteorology Unit (AMU) Quarterly Report Third Quarter FY · 13

The AMU team worked on seven tasks for their customers: (1) Ms. Crawford completed the objective lightning forecast tool for east -central Florida airports and delivered the tool and the final report to the customers. (2) Ms. Shafer continued work for Vandenberg Air Force Base on an automated tool to relate pressure gradients to peak winds. (3) Dr. Huddleston updated and delivered the tool that shows statistics on the timing of the first lightning strike of the day in the Kennedy Space Center (KSC)/Cape Canaveral Air Force Station (CCAFS) area. (4) Dr. Bauman continued work on a severe weather forecast tool focused on the Eastern Range (ER). (5) Ms. Crawford acquired the software and radar data needed to create a dual-Doppler analysis over the east-central Florida and KSC/CCAFS areas. (6) Mr. Decker continued developing a wind pairs database for the Launch Services Program to use when evaluating upper-level winds for launch vehicles. (7) Dr. Watson continued work to assimilate observational data into the high-resolution model configurations she created for Wallops Flight Facility and the ER.

Bauman, William↗

Improving the Timeliness of Winter Wheat Production Forecast in the United States of America, Ukraine and China Using MODIS Data and NCAR Growing Degree Day Information

Wheat is themost important cereal crop traded on international markets and winter wheat constitutes approximately 80% of global wheat production. Thus, accurate and timely production forecasts are critical for making informed agricultural policies and investments, as well as increasing market efficiency and stability. Becker-Reshef et al. (2010) developed an empirical generalized model for forecasting winter wheat production. Their approach combined BRDF-corrected daily surface reflectance from Moderate resolution Imaging Spectroradiometer (MODIS) Climate Modeling Grid (CMG) with detailed official crop statistics and crop typemasks. It is based on the relationship between the Normalized Difference Vegetation Index (NDVI) at the peak of the growing season, percent wheat within the CMG pixel (area within the CMG pixel occupied by wheat), and the final yields. This method predicts the yield approximately one month to six weeks prior to harvest. In this study, we include Growing Degree Day (GDD) information extracted from NCEP/NCAR reanalysis data in order to improve the winter wheat production forecast by increasing the timeliness of the forecasts while conserving the accuracy of the original model. We apply this modified model to three major wheat-producing countries: the Unites States (US), Ukraine and China from 2001 to 2012. We show that a reliable forecast can be made between one month to a month and a half prior to the peak NDVI (meaning two months to two and a half months prior to harvest), while conserving an accuracy of 10% in the production forecast.

Production↗

Research in satellite-aided crop forecasting

Evaluations of remote sensing procedures developed specifically to estimate non-U.S. spring small grains area show accuracies of less than 10 percent relative difference to reference statistics for North Dakota in 1978 and good comparison with 9000 square miles of observations over four states and Saskatchewan, Canada during the years 1976-79. Processing a 5 x 6-nautical-mile sample site requires a few minutes manual time and a few minutes central processing unit time on an AS-3000 computer. Evaluations of summer crop, corn, and soybeans area estimates show unbiased summer crops estimates in the U.S. central corn belt but significant bias in one of two years for area estimates of corn and soybeans. Based on results to date, a highly automated corn/sorghum/soybean area estimation procedure should be achieved that is applicable to Argentina.

Erickson, J. D.↗

Distortion Representation of Forecast Errors for Model Skill Assessment and Objective Analysis

We completed the formulation of the smoothness penalty functional this past quarter. We used a simplified procedure for estimating the statistics of the FCA solution spectral coefficients from the results of the unconstrained, low-truncation FCA (stopping criterion) solutions. During the current reporting period we have completed the calculation of GEOS-2 model-equivalent brightness temperatures for the 6.7 micron and 11 micron window channels used in the GOES imagery for all 10 cases from August 1999. These were simulated using the AER-developed Optimal Spectral Sampling (OSS) model.

Hoffman, Ross N.↗

Time-continuous assimilation of remote-sounding data and its effect on weather forecasting

Methods were developed for the time-continuous assimilation of satellite-sounding temperature data: direct insertion method, asynoptic successive correction method, and local linear regression method. These methods were applied to DST-6 data from the operational and experimental temperature sounders aboard the NOAA 4 and Nimbus 6 satellites. Attention is given to a comparison of these methods and their effect on the accuracy of the initial states obtained, as well as the resulting forecasts generated from these initial states. The results suggest that (1) satellite-derived temperature data can have a modest though statistically significant positive impact on numerical weather prediction in the 2-3 day range; (2) the impact is highly sensitive to the quantity of data available; and (3) the method of satellite-data assimilation can substantially influence the magnitude of the impact obtained for the same data.

Ghil, M.↗

The Use of a Satellite Climatological Data Set to Infer Large Scale Three Dimensional Flow Characteristics

Ever since the first satellite image loops from the 6.3 micron water vapor channel on the METEOSAT-1 in 1978, there have been numerous efforts (many to a great degree of success) to relate the water vapor radiance patterns to familiar atmospheric dynamic quantities. The realization of these efforts is becoming evident with the merging of satellite derived winds into predictive models (Velden et al., 1997; Swadley and Goerss, 1989). Another parameter that has been quantified from satellite water vapor channel measurements is upper tropospheric relative humidity (UTH) (e.g., Soden and Bretherton, 1996; Schmetz and Turpeinen, 1988). These humidity measurements, in turn, can be used to quantify upper tropospheric water vapor and its transport to more accurately diagnose climate changes (Lerner et al., 1998; Schmetz et al. 1995a) and quantify radiative processes in the upper troposphere. Also apparent in water vapor imagery animations are regions of subsiding and ascending air flow. Indeed, a component of the translated motions we observe are due to vertical velocities. The few attempts at exploiting this information have been met with a fair degree of success. Picon and Desbois (1990) statistically related Meteosat monthly mean water vapor radiances to six standard pressure levels of the European Centre for Medium Range Weather Forecast (ECMWF) model vertical velocities and found correlation coefficients of about 0.50 or less. This paper presents some preliminary results of viewing climatological satellite water vapor data in a different fashion. Specifically, we attempt to infer the three dimensional flow characteristics of the mid- to upper troposphere as portrayed by GOES VAS during the warm ENSO event (1987) and a subsequent cold period in 1998.

Lerner, Jeffrey A.↗

Sensitivity of the Assimilated Ozone in the UTLS to Model and Data Selection Changes

This presentation will discuss the sensitivity of assimilated ozone fields in the upper troposphere and lower stratosphere (UTLS) to a number of factors, focusing mainly on aspects of data selection and the prediction model. This is important, because assimilation represents an attempt to construct our best estimates of the true ozone field; however, inaccuracies in the UTLS ozone distribution translate into an uncertainty in factors such as the calculated radiative forcing of climate or the inferred stratosphere-troposphere exchange (STE) of ozone. The 3D ozone data assimilation system, from NASA's Global Modeling and Assimilation Office (GMAO), combines observations of total ozone column and stratospheric profiles with predictions from an off-line, parameterized chemistry and transport model (pCTM) to produce six-hourly, global analyses. The first experiments discussed assimilate ozone retrievals from the Earth-Probe Total Ozone Mapping Spectrometer (EPTOMS) and stratospheric profiles from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The SBUV/2 ozone data have a coarse vertical resolution, with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the UTLS region are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Sensitivity studies have been conducted to examine the responses to TOMS and SBUV/2 data selection, modifications of the forecast and observation error covariance models, and the model formulation (turning off chemistry or using different wind analyses in the pCTM). The second set of experiments includes an additional data type: ozone retrieved from infrared limb-emission by MIPAS on Envisat. These data offer not only improved vertical resolution in the stratosphere, but also give measurements in the polar night. Comparisons of the assimilated ozone fields from both sets of experiments with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes. It is shown that many of the changes have a significant impact on the UTLS ozone estimates. Implications for interpretation of STE and radiative forcing of climate are discussed.

Pawson, Steven↗

Introduction to Data Assimilation

Atmospheric data assimilation is a class of techniques used for producing descriptions of fields of air temperature, pressure, humidity, wind, etc. on a spatial grid or in terms of a finite functional representation. These are then used to initialize numerical weather forecasts or to analyze the atmosphere for other purposes. The techniques combine past, present, and even future observations in an approximate statistically optimal way. Various types of statistical or physically-based models and their corresponding adjoints are employed to relate diverse fields in both time and space and to relate what is observed to what is being analyzed. Computationally, the problem is very demanding and onstraining on the techniques that can be employed on a routine basis.

Atmospheric data assimilation↗

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗

The Atmospheric Infrared Sounder (AIRS) on the NASA Aqua Spacecraft: A General Remote Sensing Tool for Understanding Atmospheric Structure, Dynamics and Composition

The Atmospheric Infrared Sounder (AIRS) on the EOS Aqua Spacecraft was launched on May 4, 2002. Early in the mission, the AIRS instrument demonstrated its value to the weather forecasting community with better than 6 hours of improvement on the 5 day forecast. Now with over eight years of consistent and stable data from AIRS, scientists are able to examine processes governing weather and climate and look at seasonal and interannual trends from the AIRSdata with high statistical confidence. Naturally, long-term climate trends require a longer data set, but indications are that the Aqua spacecraft and the AIRS instrument should last beyond 2018. This paper briefly describes the AIRS data products and presents some of the most significant findings involving the use of AIRS data in the areas of weather forecast improvement, climate processes and model validation, cloud and polar processes, and atmospheric composition (chemistry and dust).

Temperature↗

A Global Climatology of Extratropical Transition Part II: Statistical Performance of the Cyclone Phase Space

This study analyzes the differences between an objective, automated identification of tropical cyclones (TCs) that undergo extratropical transition (ET), and the designation of ET determined subjectively by human forecasters in best-track data in all basins globally. The objective identification of ET is based on the cyclone phase space (CPS), calculated from the Japanese 55-year Reanalysis (JRA-55) or the ECMWF Interim Reanalysis (ERA-Interim). The resulting classification into “ET storms” and “non-ET storms” underlies the global climatology of ET presented in Part I of this study. Here, the authors investigate how well the CPS classifications agree with those in the best-track records calculated from JRA-55 or from ERA-Interim data. According to F1 scores and Matthews correlation coefficients (MCCs), the classification of ET storms in the CPS agrees best with the best-track classification in the Western North Pacific (MCC > 0.7) and the North Atlantic (MCC > 0.5). In other basins, the correlation between the CPS classification and the best-track classification is only slightly higher than that of a random classification. The JRA-55 classification achieves higher performance scores than does the ERA-Interim classification, and the differences are statistically significant in all basins. The lower performance of ERA-Interim is mainly due to a higher false alarm rate, particularly in the Eastern North Pacific. Overall, the results show that while the CPS-based classifications are good enough to be useful for many purposes, there is almost certainly room for improvement – in the representation of the storms in reanalyses, in our objective metrics of ET, and in our scientific understanding of the ET process.

Melanie Bieli↗