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 325 records · Page 18

Do Solar Coronal Holes Affect the Properties of Solar Energetic Particle Events?

The intensities and timescales of gradual solar energetic particle (SEP) events at 1 AU may depend not only on the characteristics of shocks driven by coronal mass ejections (CMEs), but also on large-scale coronal and interplanetary structures. It has long been suspected that the presence of coronal holes (CHs) near the CMEs or near the 1-AU magnetic footpoints may be an important factor in SEP events. We used a group of 41 E (is) approx. 20 MeV SEP events with origins near the solar central meridian to search for such effects. First we investigated whether the presence of a CH directly between the sources of the CME and of the magnetic connection at 1 AU is an important factor. Then we searched for variations of the SEP events among different solar wind (SW) stream types: slow, fast, and transient. Finally, we considered the separations between CME sources and CH footpoint connections from 1 AU determined from four-day forecast maps based on Mount Wilson Observatory and the National Solar Observatory synoptic magnetic-field maps and the Wang-Sheeley-Arge model of SW propagation. The observed in-situ magnetic-field polarities and SW speeds at SEP event onsets tested the forecast accuracies employed to select the best SEP/CH connection events for that analysis. Within our limited sample and the three analytical treatments, we found no statistical evidence for an effect of CHs on SEP event peak intensities, onset times, or rise times. The only exception is a possible enhancement of SEP peak intensities in magnetic clouds.

Coronal mass ejections - low coronal signatures↗

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗

Mark Twain National Forest Ecological Forecasting - Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest

The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.

Kaitlyn Bretz↗

Mark Twain National Forest Ecological Forecasting: Utilizing NASA Earth Observations to Classify Ground Cover Types in the Mark Twain National Forest

The Mark Twain National Forest (MTNF) encompasses 1.5 million acres of public land in the Ozarks region of southeastern Missouri. The industrial boom between the 1880s and 1920s had devastating effects on the shortleaf pine (Pinus echinata),Missouri’s only native pine species. The combination of fire suppression on this fire-dependent species and timber harvest of mature pine stands inhibited the development of pine seedlings and promoted the establishment of hardwood stands. Partners at the US Forest Service’s MTNF are currently involved in restoration efforts in two ranger districts, which includes removing invasive eastern red cedar (Juniperus virginiana) and prescribed burning. To expand spatial coverage for the MTNF beyond in situ observation sites, the NASA DEVELOP team analyzed land cover change from 1986 through 2019 and forecasted changes based on a ‘business-as-usual' scenario out to 2040. The team incorporated remotely sensed data from Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) into the random trees supervised classification tool in ArcGIS Pro. This tool spectrally separated pixels into five distinct land cover classes and produced classifications for 1986 and 2019, with kappa statistics of 0.87 and 0.81, respectively. Overall, there was a net decrease in conifer and meadow land cover between 1986 and 2019 along with a net increase in water, developed, and deciduous land cover. The team used TerrSet’s Land Change Modeler to forecast land cover through 2040. Results showed an increase in coniferous land cover and a decrease in deciduous cover, indicating a high probability that current restoration efforts will produce the intended effect.

DEVELOP Tech Paper↗

Evaluation of the Relative Contribution of Observing Systems in Reanalyses: Aircraft Temperature Bias and Analysis Innovations

Reanalyses have become important sources of data in weather and climate research. While observations are the most crucial component of the systems, few research projects consider carefully the multitudes of assimilated observations and their impact on the results. This is partly due to the diversity of observations and their individual complexity, but also due to the unfriendly nature of the data formats. Here, we discuss the NASA Modern-Era Retrospective analysis for Research and Applications (MERRA) and a companion dataset, the Gridded Innovations and Observations (GIO). GIO is simply a post-processing of the assimilated observations and their innovations (forecast error and analysis error) to a common spatio-temporal grid, following that of the MERRA analysis fields. This data includes in situ, retrieved and radiance observations that are assimilated and used in the reanalysis. While all these disparate observations and statistics are in a uniform easily accessible format, there are some limitations. Similar observations are binned to the grid, so that multiple observations are combined in the gridding process. The data is then implicitly thinned. Some details in the meta data may also be lost (e.g. aircraft or station ID). Nonetheless, the gridded observations should provide easy access to all the observations input to the reanalysis. To provide an example of the GIO data, a case study evaluating observing systems over the United States and statistics is presented, and demonstrates the evaluation of the observations and the data assimilation. The GIO data is used to collocate 200mb Radiosonde and Aircraft temperature measurements from 1979-2009. A known warm bias of the aircraft measurements is apparent compared to the radiosonde data. However, when larger quantities of aircraft data are available, they dominate the analysis and the radiosonde data become biased against the forecast. When AMSU radiances become available the radiosonde and aircraft analysis and forecast error take on an annual cycle. While this supports results of previous work that recommend bias corrections for the aircraft measurements, the interactions with AMSU radiances will also require further investigation. This also provides an example for reanalysis users in examining the available observations and their impact on the analysis. GIO data is presently available alongside the MERRA reanalysis.

Bosilovich, Michael G.↗

Optimizing the use of remotely sensed atmospheric data through four-dimensional data assimilation

Results with the GLA interactive forecast-retrieval-assimilation system, which directly uses satellite radiance data in the assimilation cycle, are presented. Forecasts made from initial conditions derived with this system are shown to be more skillful than those whose initial conditions came from a comparable assimilation using statistically based temperature soundings, or ones using no sounding data at all. Examples are also given of additional parameters which are produced through the physical inversion of the satellite radiance data during the interactive assimilation cycle. These are potentially useful for analysis purposes.

Pfaendtner, J.↗

Application of an Ensemble Smoother to Precipitation Assimilation

Assimilation of precipitation in a global modeling system poses a special challenge in that the observation operators for precipitation processes are highly nonlinear. In the variational approach, substantial development work and model simplifications are required to include precipitation-related physical processes in the tangent linear model and its adjoint. An ensemble based data assimilation algorithm "Maximum Likelihood Ensemble Smoother (MLES)" has been developed to explore the ensemble representation of the precipitation observation operator with nonlinear convection and large-scale moist physics. An ensemble assimilation system based on the NASA GEOS-5 GCM has been constructed to assimilate satellite precipitation data within the MLES framework. The configuration of the smoother takes the time dimension into account for the relationship between state variables and observable rainfall. The full nonlinear forward model ensembles are used to represent components involving the observation operator and its transpose. Several assimilation experiments using satellite precipitation observations have been carried out to investigate the effectiveness of the ensemble representation of the nonlinear observation operator and the data impact of assimilating rain retrievals from the TMI and SSM/I sensors. Preliminary results show that this ensemble assimilation approach is capable of extracting information from nonlinear observations to improve the analysis and forecast if ensemble size is adequate, and a suitable localization scheme is applied. In addition to a dynamically consistent precipitation analysis, the assimilation system produces a statistical estimate of the analysis uncertainty.

Zhang, Sara↗

Dynamics and Predictability of Hurricane Humberto (2007) Revealed from Ensemble Analysis and Forecasting

This study uses short-range ensemble forecasts initialized with an Ensemble-Kalman filter to study the dynamics and predictability of Hurricane Humberto, which made landfall along the Texas coast in 2007. Statistical correlation is used to determine why some ensemble members strengthen the incipient low into a hurricane and others do not. It is found that deep moisture and high convective available potential energy (CAPE) are two of the most important factors for the genesis of Humberto. Variations in CAPE result in as much difference (ensemble spread) in the final hurricane intensity as do variations in deep moisture. CAPE differences here are related to the interaction between the cyclone and a nearby front, which tends to stabilize the lower troposphere in the vicinity of the circulation center. This subsequently weakens convection and slows genesis. Eventually the wind-induced surface heat exchange mechanism and differences in landfall time result in even larger ensemble spread. 1

Sippel, Jason A.↗

Suboptimal schemes for atmospheric data assimilation based on the Kalman filter

This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

Todling, Ricardo↗

Intraseasonal oscillation of convective activity in the tropical Southern Hemisphere - May 1984-April 1986

This paper studies the intraseasonal oscillation of convection and related variables in the Southern Hemisphere tropics utilizing European Centre for Medium-range Weather Forecasts analyses of 1984 to 1986. One aspect of the research was that the original, unfiltered, time series of various variables was first analyzed to determine if a statistically significant signal occurred on the intraseasonal time scale. The analysis showed that two variables, outgoing longwave radiation and velocity potential at 200 hPa, provided the best evidence of an intraseasonal oscillation. The oscillation propagated eastward, and the convective activity for both years was more intense over the Indian Ocean-Indonesia-western Pacific region than elsewhere.

Vincent, Dayton G.↗

Perspectives on Forced Predictability

This talk addresses the relationships between weather and El Nino/Southern Oscillation (ENSO) and various other low frequency atmospheric variations such as the Pacific/North American pattern. We focus, in particular, on the predictability of extreme wintertime precipitation events over the continental United States. We first examine fifty years of daily precipitation observations and quantify the leading patterns of winter weather variability. These weather systems, familiar to operational weather forecasters, account for many of the major flooding events on the west coast as well as the major winter snowstorms along the east coast. We next examine the statistics of these storm systems with a particular focus on the occurrence of extreme events. The leading weather systems show varying degrees of linkages to low frequency atmospheric variability. We show, for example, that on seasonal time scales the probability of occurrence of certain types of west coast storms is strongly enhanced by El Nino, while it is reduced during winters with an anomalous trough to the west of the dateline. ENSO also has a strong impact on storms that develop over the Gulf of Mexico and affect much of the eastern United States, with enhanced storminess during El Nino, and reduced storminess during La Nina conditions. We conclude with an example of how well current climate models are able to reproduce the basic linkages between weather and low frequency variability.

Schubert, S.↗

Developing Empirical Lightning Cessation Forecast Guidance for the Cape Canaveral Air Force Station and Kennedy Space Center

This research addresses the 45th Weather Squadron's (45WS) need for improved guidance regarding lightning cessation at Cape Canaveral Air Force Station and Kennedy Space Center (KSC). KSC's Lightning Detection and Ranging (LDAR) network was the primary observational tool to investigate both cloud-to-ground and intracloud lightning. Five statistical and empirical schemes were created from LDAR, sounding, and radar parameters derived from 116 storms. Four of the five schemes were unsuitable for operational use since lightning advisories would be canceled prematurely, leading to safety risks to personnel. These include a correlation and regression tree analysis, three variants of multiple linear regression, event time trending, and the time delay between the greatest height of the maximum dBZ value to the last flash. These schemes failed to adequately forecast the maximum interval, the greatest time between any two flashes in the storm. The majority of storms had a maximum interval less than 10 min, which biased the schemes toward small values. Success was achieved with the percentile method (PM) by separating the maximum interval into percentiles for the 100 dependent storms.

LDAR (LIGHTNING DETECTION AND RANGING)↗

The Influence of Synoptic Scales on Low-frequency Modes of Variability

Model analysis and predictablity studies using an equivalent barotropic empirical orthogonal function (EOF) model described in Schubert (1984b) are dicussed. An extensive set of 10-day forecasts was produced for the 12 component version of the model. Comparisons with observed data show apparently useful forecasts for individual winters. However, when averaged over 10 winters, the improvement over persistence was marginal and only one of the EOF's (resembling North Atlantic blocking) showed a statistically significant improvement. These results suggest the need for a more realistic model which includes shorter time scale baroclinic processes and spatial modes capable of efficiently representing the characteristic structures associated with baroclinic instability.

Schubert, S. D.↗

A Synopsis of AIRS Global-Mean Clear-Sky Radiance Trends From 2003 to 2020

Atmospheric Infrared Sounder (AIRS) aboard the National Aeronautics and Space Administration (NASA) Aqua satellite has been operating since September 2002. Its information content, superb instrument performance, and dense sampling pattern make the AIRS radiances an invaluable data set for climate studies. The trends of global-mean, nadir-view, clear-sky AIRS radiances from 2003 to 2020 are studied here, together with the counterparts of synthetic radiances based on two reanalyzes, European Centre for Medium-Range Weather Forecasts Reanalysis V5 (ECMWF ERA5) and NASA Goddard Earth Observing System V5.4.1 (GEOS-5.4.1; a reanalysis product without assimilation of hyperspectral radiances such as AIRS). The AIRS observation shows statistically significant negative trends in most of its CO2 channels, positive but non-significant trends in the channels over the window regions, and statistically significant positive trends in some of its H2O channels. The best agreements between observed and simulated radiance trends are seen over the CO2 tropospheric channels, while the observed and simulated trends over the CO2 stratospheric channels are opposite. ERA5 results largely agree with the AIRS observation over the H2O channels. The comparison in the H2O channels helps reveal a data continuity issue in the GEOS-5.4.1. Contributions from individual variables to the radiance trends are also assessed by performing separate simulations. This study provides the first synopsis of the global-mean trend of AIRS radiances over all its thermal-IR channels.

CO2↗

Long-Range Solar Activity Predictions: A Reprieve from Cycle #24's Activity

We discuss the field of long-range solar activity predictions and provide an outlook into future solar activity. Orbital predictions for satellites in Low Earth Orbit (LEO) depend strongly on exospheric densities. Solar activity forecasting is important in this regard, as the solar ultra-violet (UV) and extreme ultraviolet (EUV) radiations inflate the upper atmospheric layers of the Earth, forming the exosphere in which satellites orbit. Rather than concentrate on statistical, or numerical methods, we utilize a class of techniques (precursor methods) which is founded in physical theory. The geomagnetic precursor method was originally developed by the Russian geophysicist, Ohl, using geomagnetic observations to predict future solar activity. It was later extended to solar observations, and placed within the context of physical theory, namely the workings of the Sun s Babcock dynamo. We later expanded the prediction methods with a SOlar Dynamo Amplitude (SODA) index. The SODA index is a measure of the buried solar magnetic flux, using toroidal and poloidal field components. It allows one to predict future solar activity during any phase of the solar cycle, whereas previously, one was restricted to making predictions only at solar minimum. We are encouraged that solar cycle #23's behavior fell closely along our predicted curve, peaking near 192, comparable to the Schatten, Myers and Sofia (1996) forecast of 182+/-30. Cycle #23 extends from 1996 through approximately 2006 or 2007, with cycle #24 starting thereafter. We discuss the current forecast of solar cycle #24, (2006-2016), with a predicted smoothed F10.7 radio flux of 142+/-28 (1-sigma errors). This, we believe, represents a reprieve, in terms of reduced fuel costs, etc., for new satellites to be launched or old satellites (requiring reboosting) which have been placed in LEO. By monitoring the Sun s most deeply rooted magnetic fields; long-range solar activity can be predicted. Although a degree of uncertainty in the long-range predictions remains, requiring future monitoring, we do not expect the next cycle's + 2-sigma value will rise significantly above solar cycle #23's activity level.

Richon, K.↗

On the dynamic estimation of relative weights for observation and forecast in numerical weather prediction

The problem of merging direct and remotely sensed (indirect) data with forecast data to get an estimate of the present state of the atmosphere for the purpose of numerical weather prediction is examined. To carry out this merging optimally, it is necessary to provide an estimate of the relative weights to be given to the observations and forecast. It is possible to do this dynamically from the information to be merged, if the correlation structure of the errors from the various sources is sufficiently different. Some new statistical approaches to doing this are described, and conditions quantified in which such estimates are likely to be good.

Wahba, Grace↗

Global-scale, intraseasonal fluctuations of diabatic forcing of the atmosphere

Fields of diabatic heating rate estimates (H) for 5-day periods were calculated from the European Center for Medium Weather Forecasts (ECMWF) analyses since 1985 as the residual of the dry thermodynamic equation. Included in these fields are the horizontal and vertical divergences of heat for both mean and eddy statistics. Previous work dealt with 4-day periods, however, with the emphasis in the Global Precipitation Climatology Project (GPCP) for 5-day period totals the change was made to accommodate the GPCP product. H has long been associated with cold tropical cloud-top temperatures as measured by polar orbiting outgoing longwave radiation (OLR) sensors. Correlations between H and OLR fields on three time scales indicate a moderate amount of agreement. For periods less than 90 days, significant negative correlations are found between H and OLR for (1) tropical and NH midlatitude oceanic areas, and (2) for zonal and hemispheric mean values. Positive correlations are seen in NH mean and continental areas of N. Africa, N. America, N. Asia and Antarctica. These latter results reflect seasonal heating and cooling. Comparisons have been made between H as H as calculated from the ECMWF analyses and output from the CCM1 T42 simulations. The CCM1 tends to have a more cellular structure with more heating (precipitation) over land versus that observed (ECMWF) over oceans.

Christy, John R.↗

Data-driven Expectations for Electromagnetic Counterpart Searches Based on LIGO/Virgo Public Alerts

Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic after glow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo’s third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics ofO3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers

Polina Petrov↗