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 451 records · Page 25

The Transition of Satellite Observations Assimilated in GEOS to JEDI

The Goddard Earth Observing System (GEOS) is used by NASA’s Global Modeling and Assimilation Office (GMAO) to produce weather, climate, and air quality forecasts and reanalysis datasets. In order to incorporate the Joint Effort for Data assimilation Integration (JEDI) system into GEOS it is necessary to provide a validated set of input observations in JEDI. The GMAO, in collaboration with the Joint Center for Satellite Data Assimilation (JCSDA), is developing the Unified Forward Operator (UFO) in JEDI and adding all the necessary features to replicate the existing capability of the Gridpoint Statistical Interpolation (GSI)–based GEOS atmospheric data assimilations system. GMAO has been adding, validating, and updating procedures to assimilate the vast array of satellite and conventional observations assimilated in GEOS, including the GEOS all-sky microwave radiance assimilation framework to assimilate those observations in the UFO. Robust tests have been conducted to ensure correct configurations of observational data bias correction, quality control, and observation error in UFO and good agreements between UFO and GSI results. Our work on satellite observations will be reported in this presentation.

Jianjun Jin↗

The Simulation and Assimilation of Doppler Wind Lidar Observations in Support of Future Instruments

With the launch of the European Space Agency's Atmospheric Dynamics Mission (ADM-Aeolus) in 2011 and the call for the 3D-Winds mission in National Research Council's decadal survey, direct spaceborne measurements of vertical wind profiles are imminent via Doppler wind lidar technology. Part of the preparedness for such missions is the development of the proper data assimilation methodology for handling such observations. Since no heritage measurements exist in space, the Joint Observing System Simulation Experiment (Joint OSSE) framework is being utilized to generate a realistic proxy dataset as a precursor to flight. These data are being used for the development of the Gridpoint Statistical Interpolation (GSI) data assimilation system utilized at a number of centers through the United States including the Global Modeling and Assimilation Office (GMAO) at NASA/Goddard Space Flight Center and at the National Centers for Environmental Prediction (NOAA/NWS/NCEP). This effort will be presented, including the methodology of proxy data generation, the handling of line-of-sight wind measurements within the GSI, and the impact on both analyses and forecasts with the addition of the new data type.

McCarty, Will↗

Forecast of an exceptionally large even-numbered solar cycle

Using the 'dynamo theory' method to predict solar activity, an accurate prediction was made for solar cycle 21 by Schatten et al. (1978). Using the same dynamo technique for solar cycle 22, a value for the smoothed sunspot number of 170 + or - 25 is obtained. This large sunspot number is expected to peak in 1990 + or - 1 year. The F(10.7) radio flux is expected to reach a smoothed value of 210 + or - 25 flux units. Since this value is larger than values obtained with prediction schemes based upon 'statistical' and 'periodicity' methods, it provides a useful test for the current methodology, based upon the strength of the sun's polar field near solar minimum. The predicted degree of solar activity is expected to enhance the density and temperature of the earth's thermosphere to values somewhat larger than those found in solar cycle 21. This will have an impact on the orbital lifetime of low altitude satellites.

Schatten, Kenneth H.↗

Comparison of the CALIPSO Level 3 Ice Cloud Product with the DARDAR and 2C-ICE products

Understanding the vertical distribution of ice clouds is crucial for climate modeling and weather forecasting. Since its launch to space in 2006, the Cloud-Aerosol Lidar with Orthogonal Polarization instrument (CALIOP) onboard the CALIPSO spacecraft has been providing unprecedented high-quality profiles of ice clouds, especially optically thin cirrus clouds, during both day and night, on a global scale. Recently the lidar science working group at the NASA Langley Research Center delivered a level 3 (L3) ice cloud product which reports monthly statistics of ice cloud extinction coefficient and ice water content on a uniform 3-dimensional spatial grid. This presentation compares the ice cloud climatology derived from this product with two other ice cloud products, namely the raDAR/liDAR(DARDAR) ice cloud product and the CloudSat-CALIPSO Ice Cloud Property Product (2C-ICE). Using processing similar to the CALIPSO L3 product, granules of DARDAR and 2C-ICE data have been processed into monthly statistics for a better comparison. Similarities and differences will be presented and discussed. The findings illustrate the strengths of the different products and help to identify which product might be most appropriate for a given research topic.

Xia Cai↗

Reconstruction of the 1997/1998 El Nino from TOPEX/POSEIDON and TOGA/TAO Data Using a Massively Parallel Pacific-Ocean Model and Ensemble Kalman Filter

Two massively parallel data assimilation systems in which the model forecast-error covariances are estimated from the distribution of an ensemble of model integrations are applied to the assimilation of 97-98 TOPEX/POSEIDON altimetry and TOGA/TAO temperature data into a Pacific basin version the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. in the first system, ensemble of model runs forced by an ensemble of atmospheric model simulations is used to calculate asymptotic error statistics. The data assimilation then occurs in the reduced phase space spanned by the corresponding leading empirical orthogonal functions. The second system is an ensemble Kalman filter in which new error statistics are computed during each assimilation cycle from the time-dependent ensemble distribution. The data assimilation experiments are conducted on NSIPP's 512-processor CRAY T3E. The two data assimilation systems are validated by withholding part of the data and quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The pros and cons of each system are discussed.

Keppenne, C. L.↗

Dynamic Hurricane Data Analysis Tool

A dynamic hurricane data analysis tool allows users of the JPL Tropical Cyclone Information System (TCIS) to analyze data over a Web medium. The TCIS software is described in the previous article, Tropical Cyclone Information System (TCIS) (NPO-45748). This tool interfaces with the TCIS database to pull in data from several different atmospheric and oceanic data sets, both observed by instruments. Users can use this information to generate histograms, maps, and profile plots for specific storms. The tool also displays statistical values for the user-selected parameter for the mean, standard deviation, median, minimum, and maximum values. There is little wait time, allowing for fast data plots over date and spatial ranges. Users may also zoom-in for a closer look at a particular spatial range. This is version 1 of the software. Researchers will use the data and tools on the TCIS to understand hurricane processes, improve hurricane forecast models and identify what types of measurements the next generation of instruments will need to collect.

Knosp, Brian W.↗

Assessing Carbon Properties in Coastal Waters with a New Observing System Testbed

Large rainfall events over land can lead to a substantial flux of carbon and nutrients to estuaries and the coastal ocean. In the mid-Atlantic on the east coast of North America (35° - 42° N), these events often happen due to tropical storm activity as well as less predictable anomalously large midlatitude storms or abrupt spring snow melt runoff. Storms can directly impact the coastal carbon cycle via export of carbon from land to sea, while also stimulating phytoplankton production due to the influx of nutrients to the coastal ocean. To assess the impact of high precipitation and river flow events on the coastal carbon cycle, we have integrated multiple observing platforms in an analytical framework to dynamically observe carbon-related properties. The new observing system testbed (NOS-T) allows for an estimate of phytoplankton and organic carbon stocks in the surface ocean, with a goal of providing near real time analytical capability. A case study of high river discharge in the mid-Atlantic in the summer and fall of 2018 and 2021 were used to examine how riverine carbon manifests along the land-estuary-ocean continuum particularly in Chesapeake Bay. Chesapeake Bay was chosen as the study site because of long term monitoring by the Chesapeake Bay Program and a robust and well-developed regional biogeochemical modeling system, the Chesapeake Bay Environmental Forecast System (CBEFS), that is publicly available. There is also new capability in Chesapeake Bay to combine hyperspectral radiometric data from the Aerosol Robotic Network site that is online in 2022 with in-water observations of optical properties. Precipitation data from the Global Precipitation Measurement IMERG data set was used to establish triggering criteria for storm carbon flux observation. USGS and in-water carbon data were used to establish statistical models to estimate the mass flux of organic carbon into Chesapeake Bay using measured river discharge. The year 2018 was examined as there was record rainfall and near record river flow in the late summer. The CBEFS modeling was compared with in-water estimates of dissolved organic carbon (DOC) to assess the model’s ability to capture storm fluxes of carbon to the ecosystem. Multiple satellite platforms were also used to assess how remote sensing using passive Earth orbiting sensors can be used to observe carbon in this complex coastal region. A set of recommendations have been established to improve sensing capability to measure aquatic carbon during storm events. A unique challenge in these dynamic inland waters is how rapidly carbon fluxes can evolve in space and time, with many sources contributing to the water leaving reflectance that satellites can observe. To improve our ability to quantify carbon stocks and fluxes in near real time, a suite of satellite sensors and high-resolution modeling capability is needed, all supplemented by in situ monitoring. Future inclusion of in-water observations that would be deployed when the system reaches triggering criteria, as well as taskable orbital instruments, will improve estimates of ocean carbon properties and provide the ability to calculate major carbon stocks and fluxes in near real time following episodic storm events.

land-ocean continuum↗

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗

Technical report series on global modeling and data assimilation. Volume 4: Documentation of the Goddard Earth Observing System (GEOS) data assimilation system, version 1

This report describes the analysis component of the Goddard Earth Observing System, Data Assimilation System, Version 1 (GEOS-1 DAS). The general features of the data assimilation system are outlined, followed by a thorough description of the statistical interpolation algorithm, including specification of error covariances and quality control of observations. We conclude with a discussion of the current status of development of the GEOS data assimilation system. The main components of GEOS-1 DAS are an atmospheric general circulation model and an Optimal Interpolation algorithm. The system is cycled using the Incremental Analysis Update (IAU) technique in which analysis increments are introduced as time independent forcing terms in a forecast model integration. The system is capable of producing dynamically balanced states without the explicit use of initialization, as well as a time-continuous representation of non- observables such as precipitation and radiational fluxes. This version of the data assimilation system was used in the five-year reanalysis project completed in April 1994 by Goddard's Data Assimilation Office (DAO) Data from this reanalysis are available from the Goddard Distributed Active Center (DAAC), which is part of NASA's Earth Observing System Data and Information System (EOSDIS). For information on how to obtain these data sets, contact the Goddard DAAC at (301) 286-3209, EMAIL daac@gsfc.nasa.gov.

Suarez, Max J.↗

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗

Grid Resolution Study over Operability Space for a Mach 1.7 Low Boom External Compression Inlet

This paper presents a statistical methodology whereby the probability limits associated with CFD grid resolution of inlet flow analysis can be determined which provide quantitative information on the distribution of that error over the specified operability range. The objectives of this investigation is to quantify the effects of both random (accuracy) and systemic (biasing) errors associated with grid resolution in the analysis of the Lockheed Martin Company (LMCO) N+2 Low Boom external compression supersonic inlet. The study covers the entire operability space as defined previously by the High Speed Civil Transport (HSCT) High Speed Research (HSR) program goals. The probability limits in terms of a 95.0% confidence interval on the analysis data were evaluated for four ARP1420 inlet metrics, namely (1) total pressure recovery (PFAIP), (2) radial hub distortion (DPH/P), (3) ) radial tip distortion (DPT/P), and (4) ) circumferential distortion (DPC/P). In general, the resulting +/-0.95 delta Y interval was unacceptably large in comparison to the stated goals of the HSCT program. Therefore, the conclusion was reached that the "standard grid" size was insufficient for this type of analysis. However, in examining the statistical data, it was determined that the CFD analysis results at the outer fringes of the operability space were the determining factor in the measure of statistical uncertainty. Adequate grids are grids that are free of biasing (systemic) errors and exhibit low random (precision) errors in comparison to their operability goals. In order to be 100% certain that the operability goals have indeed been achieved for each of the inlet metrics, the Y+/-0.95 delta Y limit must fall inside the stated operability goals. For example, if the operability goal for DPC/P circumferential distortion is ≤0.06, then the forecast Y for DPC/P plus the 95% confidence interval on DPC/P, i.e. +/-0.95 delta Y, must all be less than or equal to 0.06.

Grid↗

General-aviation's view of progress in the aviation weather system

For all its activity statistics, general-aviation is the most vulnerable to hazardous weather. Of concern to the general aviation industry are: (1) the slow pace of getting units of the Automated Weather Observation System (AWOS) to the field; (2) the efforts of the National Weather Service to withdraw from both the observation and dissemination roles of the aviation weather system; (3) the need for more observation points to improve the accuracy of terminal and area forecasts; (4) the need for improvements in all area forecasts, terminal forecasts, and winds aloft forecasts; (5) slow progress in cockpit weather displays; (6) the erosion of transcribed weather broadcasts (TWEB) and other deficiencies in weather information dissemination; (7) the need to push to make the Direct User Access Terminal (DUAT) a reality; and (7) the need to improve severe weather (thunderstorm) warning systems.

Lundgren, Douglas J.↗

Analysis/forecast experiments with a multivariate statistical analysis scheme using FGGE data

A three-dimensional, multivariate, statistical analysis method, optimal interpolation (OI) is described for modeling meteorological data from widely dispersed sites. The model was developed to analyze FGGE data at the NASA-Goddard Laboratory of Atmospherics. The model features a multivariate surface analysis over the oceans, including maintenance of the Ekman balance and a geographically dependent correlation function. Preliminary comparisons are made between the OI model and similar schemes employed at the European Center for Medium Range Weather Forecasts and the National Meteorological Center. The OI scheme is used to provide input to a GCM, and model error correlations are calculated for forecasts of 500 mb vertical water mixing ratios and the wind profiles. Comparisons are made between the predictions and measured data. The model is shown to be as accurate as a successive corrections model out to 4.5 days.

Baker, W. E.↗

Use NU-WRF and GCE Model to Simulate the Precipitation Processes During MC3E Campaign

One of major CRM approaches to studying precipitation processes is sometimes referred to as "cloud ensemble modeling". This approach allows many clouds of various sizes and stages of their lifecycles to be present at any given simulation time. Large-scale effects derived from observations are imposed into CRMs as forcing, and cyclic lateral boundaries are used. The advantage of this approach is that model results in terms of rainfall and QI and Q2 usually are in good agreement with observations. In addition, the model results provide cloud statistics that represent different types of clouds/cloud systems during their lifetime (life cycle). The large-scale forcing derived from MC3EI will be used to drive GCE model simulations. The model-simulated results will be compared with observations from MC3E. These GCE model-simulated datasets are especially valuable for LH algorithm developers. In addition, the regional scale model with very high-resolution, NASA Unified WRF is also used to real time forecast during the MC3E campaign to ensure that the precipitation and other meteorological forecasts are available to the flight planning team and to interpret the forecast results in terms of proposed flight scenarios. Post Mission simulations are conducted to examine the sensitivity of initial and lateral boundary conditions to cloud and precipitation processes and rainfall. We will compare model results in terms of precipitation and surface rainfall using GCE model and NU-WRF

Tao, Wei-Kuo↗

Forecasting Coronal Mass Ejections from Magnetograms

We report further results from our ongoing assessment of magnetogram-based measures of active-region nonpotentiality (magnetic shear and twist) and size as predictors of coronal mass ejections (CMEs). From a set of 36 vector magnetograms of bipolar active regions (Falconer, Moore, & Gary 2004, ApJ, submitted), we have found that (1) each of our 5 measures of active-region nonpotentiality has a 75-80% success rate in predicting whether an active region will produce a CME within 2 days after the magnetogram, (2) the correlation of each of these measures with CME production in this window is statistically significant (confidence level greater than 95%), (3) our measure of active-region size has a 65% success rate in predicting CMEs in this window, but (4) the correlation between active-region size and CME productivity is not statistically significant (confidence level approximately 80%). As part of the work under our pending 2003 LWS TR&T proposal, we will double our sample to approximately 70 active regions in order to demonstrate a statistically significant correlation between active-region size and CME productivity, and to determine whether this correlation is a byproduct of any Sun-produced correlation between magnetic twist and size of active regions. Since the 2002 LWS Science Workshop, we have found a measure of active-region nonpotentiality (the length of strong-gradient neutral line) that can be measured from a line-of-sight magnetogram, and we have generalized this measure and the corresponding measure from a vector magnetogram (the length of strong-shear neutral line) so that they can be applied to multiple-bipole complex active regions as well as to bipolar active regions. Preliminary results will be presented from application of these two generalized measures to our sample of bipolar active regions and to a new sample of multiple-bipole active regions.

Falconer, D. A.↗

Clear air turbulence forecasting techniques

A method to improve clear air turbulence (CAT) forecasting by more effectively using the currently operational rawinsonde (RW) system is discussed. The method is called the Diagnostic Richardson Number Tendency (DRT) technique. The technique does not attempt to use the RW as a direct detector of the turbulent motion or even of the CAT mechanism structure but rather senses the synoptic scale centers of action which provide the energy to the CAT mechanism at the mesoscale level. The DRT algorithm is deterministic rather than statistical in nature, using the hydrodynamic equations (equations of motion) relevant to the synoptic scale. However, interpretation, by necessity, is probabilistic. What is most important with respect to its operational implementation is that this method uses the same input data as currently used by the operational National Meteorological Center prognostic models.

Keller, J. L.↗

New Releases and Upgrades of Giant Planet Global Reference Atmospheric Models (GRAMs)

• GRAMs are engineering-oriented atmospheric models that estimate mean values and statistical variations of atmospheric properties for numerous planetary destinations – Outputs atmospheric density, temperature, pressure, winds, and chemical composition along a user-defined path – Provide mean values and variability for any point in an atmosphere – Includes seasonal, geographic, and altitude variations – Used by engineering community because of their ability to create realistic atmospheric dispersions; can be integrated into high fidelity flight dynamic simulations of launch, entry, descent, and landing (EDL), aerobraking, and aerocapture • GRAMs are not forecast models • GRAMs are available for: Earth, Mars, Venus, Neptune, Titan, Jupiter, and Uranus • Available through the NASA Software Catalog https://software.nasa.gov/software/MFS-33888-1

atmospheric models↗

Assessing the Usefulness of AIRS Radiance Observations in a 4D-Var Assimilation Scheme Using the Penn State/NCAR Mesoscale Model Version 5 (MM5) and a Stand Alone Radiative Transfer Algorithm (SARTA)

Short-term precipitation forecast skill remains poor despite further advances in numerical modeling and data analysis. Utilizing indirect observations, such as the radiance observations obtained from the Atmospheric Infrared Sounder (AIRS), in a 4D-Var assimilation scheme may lead to a rapid improvement in the skill of 12-h QPF. In order to investigate the usefulness of AIRS observations, the adjoint code of a radiative transfer model (SARTA) is developed and linked to MM5 4D-Var system. A statistical comparison of the AIRS radiances with those calculated using the linked MMS-SARTA model is carried out to assess errors in the radiative transfer model, MM5 and AIRS data, and to develop a quality control scheme. Outliers and any bias that may exist in radiative transfer model, MM5 and AIRS data are identified. An adjoint sensitivity study, utilizing the SARTA adjoint model, is then done to identify the relative sensitivities of each spectral channel to such parameters as atmospheric temperature, pressure, and mixing ratio. Another sensitivity study is done using the MM5 alone to determine the sensitivity of the 12-h QPF to model parameters such as atmospheric temperature, pressure, and mixing ratio under various weather regimes (i.e. summertime convection, winter snowfall). Results from these two types of adjoint sensitivity studies are used to determine which AIRS channels contribute most to 12-h QPF under differing weather regimes. This information is useful when conducting 4D-Var experiments using a subset of the 2378 AIRS spectral channels. Finally, 4D-Var experiments are done to analyze the impact of the AIRS observations on the quality of the 12-h QPF results. This project falls under the data assimilation general session topic as specified by the conference organizers.

Carrier, M.↗