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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Cosmology with persistent homology: a Fisher forecast

Abstract Persistent homology naturally addresses the multi-scale topological characteristics of the large-scale structure as a distribution of clusters, loops, and voids. We apply this tool to the dark matter halo catalogs from theQuijotesimulations, and build a summary statistic for comparison with the joint power spectrum and bispectrum statistic regarding their information content on cosmological parameters and primordial non-Gaussianity. Through a Fisher analysis, we find that constraints from persistent homology are tighter for 8 out of the 10 parameters by margins of 13–50%. The complementarity of the two statistics breaks parameter degeneracies, allowing for a further gain in constraining power when combined. We run a series of consistency checks to consolidate our results, and conclude that our findings motivate incorporating persistent homology into inference pipelines for cosmological survey data.

Astronomy & Astrophysics↗

NASA NTRS2003

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

NASA NTRS1988

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↗

NASA NTRS1990

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

NASA NTRS2022

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↗

NASA NTRS2015

High-Resolution Mesoscale Model Setup for the Eastern Range and Wallops Flight Facility

Mesoscale weather conditions can have an adverse effect on space launch, landing, ground processing, and weather advisories, watches, and warnings at the Eastern Range (ER) in Florida and Wallops Flight Facility (WFF) in Virginia. During summer, land-sea interactions across Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) lead to sea breeze front formation, which can spawn deep convection that can hinder operations and endanger personnel and resources. Many other weak locally-driven low-level boundaries and their interactions with the sea breeze front and each other can also initiate deep convection in the KSC/CCAFS area. These convective processes often last 60 minutes or less and pose a significant challenge to the local forecasters. Surface winds during the transition seasons (spring and fall) pose the most difficulties for the forecasters at WFF. They also encounter problems forecasting convective activity and temperature during those seasons. Therefore, accurate mesoscale model forecasts are needed to better forecast a variety of unique weather phenomena. Global and national scale models cannot properly resolve important local-scale weather features at each location due to their horizontal resolutions being much too coarse. Therefore, a properly tuned local data assimilation (DA) and forecast model at a high resolution is needed to provide improved capability. To accomplish this, a number of sensitivity tests were performed using the Weather Research and Forecasting (WRF) model in order to determine the best DA/model configuration for operational use at each of the space launch ranges to best predict winds, precipitation, and temperature. A set of Perl scripts to run the Gridpoint Statistical Interpolation (GSI)/WRF in real-time were provided by NASA's Short-term Prediction Research and Transition Center (SPoRT). The GSI can analyze many types of observational data including satellite, radar, and conventional data. The GSI/WRF scripts use a cycled GSI system similar to the operational North American Mesoscale (NAM) model. The scripts run a 12-hour pre-cycle in which data are assimilated from 12 hours prior up to the model initialization time. A number of different model configurations were tested for both the ER and WFF by varying the horizontal resolution on which the data assimilation was done. Three different grid configurations were run for the ER and two configurations were run for WFF for archive cases from 27 Aug 2013 through 10 Nov 2013. To quantify model performance, standard model output will be compared to the Meteorological Assimilation Data Ingest System (MADIS) data. The MADIS observation data will be compared to the WRF forecasts using the Model Evaluation Tools (MET) verification package. In addition, the National Centers for Environmental Prediction's Stage IV precipitation data will be used to validate the WRF precipitation forecasts. The author will summarize the relative skill of the various WRF configurations and how each configuration behaves relative to the others, as well as determine the best model configuration for each space launch range.

Watson, Leela R.↗

DOE OSTI2021

A Statistical Evaluation of WRF-LES Trace Gas Dispersion Using Project Prairie Grass Measurements

In recent years, new measurement systems have been deployed to monitor and quantify methane emissions from the natural gas sector. Large-eddy simulation (LES) has complemented measurement campaigns by serving as a controlled environment in which to study plume dynamics and sampling strategies. However, with few comparisons with controlled-release experiments, the accuracy of LES for modeling natural gas emissions is poorly characterized. In this paper, we evaluate LES from the Weather Research and Forecasting (WRF) Model against Project Prairie Grass campaign measurements and surface layer similarity theory. Using WRF-LES, we simulate continuous emissions from 30 near-surface trace gas sources in two stability regimes: strong convection and weak convection. We examine the impact of grid resolutions ranging from 6.25 to 52 m in the horizontal dimension on model results. We evaluate performance in a statistical framework, calculating fractional bias and conducting Welch’s t tests. WRF-LES accurately simulates observed surface concentrations at 100 m and beyond under strong convection; simulated concentrations pass t tests in this region irrespective of grid resolution. However, in weakly convective conditions with strong winds, WRF-LES substantially overpredicts concentrations—the magnitude of fractional bias often exceeds 30%, and all but one t test fails. The good performance of WRF-LES under strong convection correlates with agreement with local free convection theory and a minimal amount of parameterized turbulent kinetic energy. The poor performance under weak convection corresponds to misalignment with Monin–Obukhov similarity theory and a significant amount of parameterized turbulent kinetic energy.

17 WIND ENERGY↗

DOE OSTI2025

Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by f NL ⁠, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure f NL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter haloes. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of f NL ⁠. In LSS surveys, the effect of imaging systematics on two-point statistics is often degenerate with the PNG signal. Our proposed model employs three-point statistics primarily to break this degeneracy. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our method. Our forecast indicates the ability to constrain f NL to a precision of σf NL ≈ 22 with one year of DESI survey data, as well as the ability to constrain the imaging systematic weights in situ.

early Universe↗

DOE OSTI2026

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗

DOE OSTI2024

Anticipating Technical Expertise and Capability Evolution in Research Communities Using Dynamic Graph Transformers

The ability to anticipate global technical expertise and capability evolution trends is essential for national and global security, especially in safety-critical domains such as nuclear nonproliferation (NN) and rapidly emerging fields like artificial intelligence (AI). Here, in this work, we extend traditional statistical relational learning approaches (e.g., link prediction in collaboration networks) and formulate a problem of anticipating technical expertise and capability evolution using dynamic heterogeneous graph representations. We develop novel capabilities to forecast collaboration patterns, authorship behavior, and technical capability evolution at different granularities (e.g., scientist and institution levels) in two distinct research fields. We implement a dynamic graph transformer (DGT) neural architecture, which pushes the state-of-the-art graph neural network models by: 1) forecasting heterogeneous (rather than homogeneous) nodes and edges; and 2) relying on both discrete- and continuous-time inputs. We demonstrate that our DGT models predict collaboration, partnership, and expertise patterns with 0.26, 0.73, and 0.53 mean reciprocal rank values for AI and 0.48, 0.93, and 0.22 for NN domains. DGT model performance exceeds the best-performing static graph baseline models by 30%–80% across AI and NN domains. Our findings demonstrate that DGT models boost inductive task performance when previously unseen nodes appear in the test data for the domains with emerging collaboration patterns (e.g., AI). Specifically, models accurately predict which established scientists will collaborate with early career scientists and vice versa in the AI domain.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Machine learning of factors for improving oyster hatchery production

Oyster aquaculture and restoration in the Chesapeake Bay are vital, yet hatcheries frequently struggle with inconsistent larval growth and sudden mass mortality events. Unpredictable disruptions in larval production cause large economic losses, represent a perceived risk to growers, and impede industry expansion. To better understand associations between production yield and its potential predictors, we applied machine learning (random forest, and neural network) and statistical (generalized additive model) models to a comprehensive dataset of environmental, water quality, and operational parameters from a Maryland oyster hatchery, aiming to identify key yield predictors and develop a robust forecasting tool. We used recursive Boruta algorithm for variable selection, pinpointing critical predictors, and employed cross-validation to fine-tune model settings. Shapley value analysis offered crucial insights into model interpretations, highlighting week number, Normalized Difference Vegetation Index, salinity, turbidity, and fecundity as primary drivers of yield variability. For low-yield cases, salinity-related variables were particularly important. Our findings provide an early warning system for potential production downturns, empowering hatchery operators to make data-driven decisions for optimizing water conditions, feeding schedules, and broodstock management. By boosting predictability and efficiency, this research directly supports economic stability of the oyster industry and ecological health of the Chesapeake Bay.

Vishwakarma, Srishti [Oak Ridge National Laborator↗

NASA NTRS1985

Synoptic scale forecast skill and systematic errors in the MASS 2.0 model

The synoptic scale performance characteristics of MASS 2.0 are determined by comparing filtered 12-24 hr model forecasts to same-case forecasts made by the National Meteorological Center's synoptic-scale Limited-area Fine Mesh model. Characteristics of the two systems are contrasted, and the analysis methodology used to determine statistical skill scores and systematic errors is described. The overall relative performance of the two models in the sample is documented, and important systematic errors uncovered are presented.

Koch, S. E.↗

NASA NTRS2004

Specification of ISS Plasma Environment Variability

Quantifying spacecraft charging risks and associated hazards for the International Space Station (ISS) requires a plasma environment specification for the natural variability of ionospheric temperature (Te) and density (Ne). Empirical ionospheric specification and forecast models such as the International Reference Ionosphere (IRI) model typically only provide long term (seasonal) mean Te and Ne values for the low Earth orbit environment. This paper describes a statistical analysis of historical ionospheric low Earth orbit plasma measurements from the AE-C, AE-D, and DE-2 satellites used to derive a model of deviations of observed data values from IRI-2001 estimates of Ne, Te parameters for each data point to provide a statistical basis for modeling the deviations of the plasma environment from the IRI model output. Application of the deviation model with the IRI-2001 output yields a method for estimating extreme environments for the ISS spacecraft charging analysis.

Minow, Joseph I.↗

NASA NTRS2009

Decadal Variation of the Number of El Nino Onsets and El Nino-Related Months and Estimating the Likelihood of El Nino Onset in a Warming World

Examination of the decadal variation of the number of El Nino onsets and El Nino-related months for the interval 1950-2008 clearly shows that the variation is better explained as one expressing normal fluctuation and not one related to global warming. Comparison of the recurrence periods for El Nino onsets against event durations for moderate/strong El Nino events results in a statistically important relationship that allows for the possible prediction of the onset for the next anticipated El Nino event. Because the last known El Nino was a moderate event of short duration (6 months), having onset in August 2006, unless it is a statistical outlier, one expects the next onset of El Nino probably in the latter half of 2009, with peak following in November 2009-January 2010. If true, then initial early extended forecasts of frequencies of tropical cyclones for the 2009 North Atlantic basin hurricane season probably should be revised slightly downward from near average-to-above average numbers to near average-to-below average numbers of tropical cyclones in 2009, especially as compared to averages since 1995, the beginning of the current high-activity interval for tropical cyclone activity.

Wilson, Robert M.↗

NASA NTRS2012

P161 Improved Impact of Atmospheric Infrared Sounder (AIRS) Radiance Assimilation in Numerical Weather Prediction

For over 6 years, AIRS radiances have been assimilated operationally into National (e.g. Environmental Modeling Center (EMC)) and International (e.g. European Centre for Medium-Range Weather Forecasts (ECMWF)), operational centers; assimilated in the North American Mesoscale (NAM) since 2008. Due partly to data latency and operational constraints, hyperspectral radiance assimilation has had less impact on the Gridpoint Statistical Interpolation (GSI) system used in the NAM and GFS. Objective of this project is to use AIRS retrieved profiles as a proxy for the AIRS radiances in situations where AIRS radiances are unable to be assimilated in the current operational system by evaluating location and magnitude of analysis increments.

Zavodsky, Bradley T.↗

DOE OSTI2026

Simulation budgeting for hybrid effective field theories

In this work, we forecast the number of, and requirements on, N-body simulations needed to train hybrid effective field theory (HEFT) emulators for a range of use cases, using a hybrid of HMcode and perturbation theory as a surrogate model. Our accuracy goals, determined with careful consideration of statistical and systematic uncertainties, are 1% accurate in the high-likelihood range of cosmological parameters, and 2% accurate over a broader parameter space volume for k < 1 h Mpc -1 and z < 3. Focusing in part on the 8-parameter w 0 w a CDM+ m ν cosmological model, we find that < 225 simulations are required to meet our error goals over our wide parameter space, including models with rapidly evolving dark energy, given our simulation and emulator recommendations. For a more restricted parameter space volume, as few as 80 simulations are sufficient. We additionally present simulation forecasts for example use cases, and make the code used in our analyses publicly available. These results offer practical guidance for efficient emulator design and simulation budgeting in future cosmological analyses.

cosmological parameters from LSS↗

NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗

NASA NTRS2013

Analysis of Runway Incursion Data

A statistical analysis of runway incursion (RI) events was conducted to ascertain relevance to the top ten challenges of the National Aeronautics and Space Administration Aviation Safety Program (AvSP). The information contained in the RI database was found to contain data that may be relevant to several of the AvSP top ten challenges. When combined with other data from the FAA documenting air traffic volume from calendar year 2000 through 2011, the structure of a predictive model emerges that can be used to forecast the frequency of RI events at various airports for various classes of aircraft and under various environmental conditions.

Runway incursions↗