Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatiotemporal modeling”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Constraining Black Carbon Aging in Global Models to Reflect Timescales for Internal Mixing

The radiative effects of black carbon depend critically on its atmospheric lifetime, which is controlled by the rate at which freshly emitted combustion particles become internally mixed with other aerosol components. Global aerosol models strive to represent this process, but the timescale for aerosol mixing is not easily constrained using observations. In this study, we apply a timescale parameterization derived from particle‐resolved simulations to quantify, in a global aerosol model, the timescale for internal mixing. We show that, while highly variable, the average timescale for internal mixing is approximately 3 hr, which is much shorter than the 24‐hr aging timescale traditionally applied in bulk aerosol models. We then use the mixing timescale to constrain the aging criterion in the Modal Aerosol Module. Our analysis reveals that, to best reflect timescales for internal mixing, modal models should assume that particles transition from the hydrophobic (fresh) to the hydrophilic (aged) class once they accumulate a coating thickness equal to four monolayers of sulfuric acid, as opposed to the model's current aging criterion of eight monolayers. We show that, in remote regions like the Arctic and Antarctic, predictions of black carbon loading and its seasonal variation are particularly sensitive to the model representation of aging. By constraining aging in global models to reflect mixing timescales simulated by the particle‐resolved model, we eliminate one of the free parameters governing black carbon's long‐range transport and spatiotemporal distribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗
NASA NTRS2002

Tracking Avian Reservoirs of Arboviruses using Remote Sensing and Radiotelemetry

Encephalitis is caused by a virus that is transmitted by mosquitoes between mammalian hosts. The virus is closely related to the West Nile virus (WNV), which started in New York in 1999, and has since spread to 25 states. Like encephalitis, WNV is vectored by mosquitoes, and the primary hosts are birds; humans are accidental, or'dead-end' hosts. Very little is understood about the behavior of these bird populations, and how they intersect - both in time and in space - with mosquito populations. Exploring these relationships is the first step in developing models for encephalitis and WNV transmission risk. This project combines remotely sensed data with radiotelemetry to create a spatiotemporal map of encephalitis viral activity in bird and mosquito populations in the Sacramento Valley of California. Specifically, remote sensing (RS) and geographic information system (GIS) technologies were used to characterize habitats utilized by both avian viral reservoirs and the mosquitoes that vector encephalitis. Radiotelemetry and serosurveys (blood) were then used to spatially and temporally track the patterns of infection. The project uses Landsat ETM+ multitemporal satellite data to characterize habitats utilized by both birds and the mosquito vectors. Mist nets were used to sample members of individual flocks of blackbirds and cowbirds over a period of several months; these birds were then bled to assess their viral status, banded, and fitted with transmitters. Radiotelemetry was used to spatially and temporally track the distribution of banded birds and their associated flocks. The movement of these indicator flocks were compared with the location of remotely sensed (adult and larval) mosquito habitats to determine the intersection of bird's and vectors; this is key in understanding where and when transmission occurs from bird to bird, as well as from bird to mammal, via mosquito. The relationships found during the project are being used to generate a model of encephalitis transmission risk in California.

Beck, L.↗
NASA NTRS2018

Climatespark: an In-Memory Distributed Computing Framework for Big Climate Data Analytics

The unprecedented growth of climate data creates new opportunities for climate studies, and yet big climate data pose a grand challenge to climatologists to efficiently manage and analyze big data. The complexity of climate data content and analytical algorithms increases the difficulty of implementing algorithms on high performance computing systems. This paper proposes an in-memory, distributed computing framework, ClimateSpark, to facilitate complex big data analytics and time-consuming computational tasks. Chunking data structure improves parallel I/O efficiency, while a spatiotemporal index is built for the chunks to avoid unnecessary data reading and preprocessing. An integrated, multi-dimensional, array-based data model (ClimateRDD) and ETL operations are developed to address big climate data variety by integrating the processing components of the climate data lifecycle. ClimateSpark utilizes Spark SQL and Apache Zeppelin to develop a web portal to facilitate the interaction among climatologists, climate data, analytic operations and computing resources (e.g., using SQL query and Scala/Python notebook). Experimental results show that ClimateSpark conducts different spatiotemporal data queries/analytics with high efficiency and data locality. ClimateSpark is easily adaptable to other big multiple- dimensional, array-based datasets in various geoscience domains.

Hu, Fei↗
DOE OSTI2026

Integrating Maximum Entropy Production Theory and Machine Learning to Improve Global Evapotranspiration Modeling

Accurate estimation of terrestrial evapotranspiration (ET) is vital for understanding global water and energy cycles. However, current global ET estimations are not well constrained. This study introduces an integrated framework combining the Maximum Entropy Production (MEP) theory with Random Forest (RF) model to improve global ET estimation. Specifically, in contrast to direct ET estimation by the RF model, the integrated framework (MEP‐RF) trains to predict error of MEP‐simulated ET. MEP‐RF outperforms RF in spatiotemporal extrapolation. Attribution analysis with in situ observations reveals that the inputs of MEP are the most critical variables for the ET process, including net radiation, vegetated area, soil moisture, and surface temperature. We further drive MEP‐RF with global reanalysis and satellite data sets of these four inputs, yielding a global mean terrestrial ET of 548 mm/year, with 77% attributed to transpiration. The global ET increased at a rate of 0.85 mm/year per year during 2003–2021, primarily due to vegetation greening rather than rising temperature, while decreasing soil moisture led to decreasing regional ET. The integrated framework provides a novel approach for the estimation of global ET without the need for hard‐to‐obtain and thus uncertain inputs, such as wind speed, surface roughness, aerodynamic and canopy stomatal resistance. Therefore, MEP‐RF offers an independent method on existing global ET products. It represents a promising physically based approach that can be incorporated into Earth System Models to enhance water and energy cycle simulations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2019

Evaluation Analysis of NASA SMAP L3 and L4 and SPoRT-LIS Soil Moisture Data in the United States

Soil moisture has a critical role in the development, frequency and persistence of climatic and hydrologic extremes such as drought, heat wave and flooding events. In situ soil moisture data are uneven and sparse in time and space. This highlights the need to utilize other soil moisture sources to fill this spatiotemporal gap. The goal of this study is to validate one satellite-based and two model-based soil moisture datasets with in situ data across the United States. Soil moisture information from the Soil Moisture Active Passive (SMAP) enhanced level 3 (L3) (SMAP L3) and modeled level 4 (L4) (SMAP L4) data at 9-km resolution and Short-term Prediction Research and Transition-Land Information System (SPoRT-LIS) modeled at 3-km resolution were selected for evaluation. SPoRT-LIS is a near real-time, high resolution operational land analysis data. Ground-based data were obtained from the North American Soil Moisture Database (NASMD) for 362 stations. Seven statistical indicators including anomaly, Spearman, and Pearson correlation coefficients, the systematic error (Bias), root mean square error (RMSE), unbiased root mean square error (ubRMSE), and normalized standard deviation (SDV) were used to evaluate the satellite- and model-derived soil moisture data. In addition, the triple collocation (TC) error model was used to measure the error among SMAP L4, SPoRT-LIS and ground-based data. This study assesses which satellite or modeled dataset is most appropriate for specific times and locations to use as a surrogate for in situ observations. Temporal and spatial analysis demonstrated that, overall, SMAP L4 performed better than SMAP L3 and SPoRT-LIS. Strong agreement was observed between SMAP L4 and in situ observations (ρ = 0.53, Bias = −0.006) in all seasons and most regions with various land covers, especially in winter and in the central regions of the United States. For croplands, SMAP L4 presented the best agreement with in situ data, analyzing all period (ρ = 0.60) and non-winter period (ρ = 0.61) separately.

Tavakol, Ameneh↗
NASA NTRSDate not supplied

Accelerating TEMPO Air Quality Science Through STAQS

Soon after the launch of Tropospheric Emissions: Monitoring of Pollution (TEMPO) mission, NASA is supporting the Synergistic TEMPO Air Quality Science field study (STAQS) in summer 2023. This study’s main objective is to accelerate science with geostationary air quality observations from TEMPO to better understand its use in air pollution research and applications. Science objectives include, but are not limited to, TEMPO L2 product evaluation, the interpretation of the spatiotemporal evolution of TEMPO data during air quality events, assessment of anthropogenic emissions, contributions to chemical transport modeling, and environmental justice applications. These objectives can be accomplished through the integration of satellite data with systematically repeated high-resolution aircraft- and ground-based measurements in multiple urban environments. The primary urban cities considered include Los Angeles, New York City, and Chicago in June-August 2023. Measurements consist of airborne remote sensing observations of air quality constituents of nitrogen dioxide (NO2), formaldehyde (HCHO), ozone, and aerosols using the GeoCape Airborne Simulator and High-Spectral Resolution Lidar 2/Differential Absorption Lidar (HSRL2/DIAL) on the NASA JSC G-V aircraft, greenhouse gas observations from Airborne Visible InfraRed Imaging Spectrometer - Next Generation (AVIRIS-NG) and High Altitude Lidar Observatory (HALO) on the NASA LaRC G-III aircraft, and ground-based remote sensing and in situ observations of ozone, NO2, and HCHO from the Tropospheric Ozone Lidar Network (TOLNet), Pandora spectrometers, and ground-based monitoring networks. Key partnership studies providing in situ airborne observations include the NOAA Atmospheric Emissions and Reactions Observed from Megacities to Marine Areas (AEROMMA) field study on the NASA DC-8 in the same primary target areas as STAQS and Greater New York Oxidant, Tropospheric Halogens, and Aerosol Measurements and Modeling (GOTHAMM) field study on the NSF C-130 near New York City. This presentation will include a current status update of the STAQS mission and an overview of its measurement strategies and science objectives with the goal of promoting continued discussions for building and strengthening collaborations prior to the mission.

Laura Judd↗
DOE OSTI2024

Exploring impacts of electricity tariff on charging infrastructure planning: An activity-based approach

In the past decade, electric vehicles (EVs) have gained popularity for their efficiency and environmental benefits. Advances in battery technology and charging equipment have yielded long-range EVs and fast-charging. However, many major cities lack adequate charging infrastructure for daily EV use. This study addresses this gap by integrating activity-based modeling, charging behavior simulation, and charging infrastructure optimization. The research utilizes the POLARIS agent-based transportation model to accurately capture user activities, trip patterns, and traffic flows. Additionally, the study investigates the impact of fixed and spatiotemporal electricity rate distributions on optimal charging infrastructure deployment. The framework is applied to the Chicago regional area network and analyzed under various EV ownership scenarios. Further, the results reveal significant impacts of the charging pricing strategy on user decision-making and charging demand distribution. There is also a need for consistent pricing policies in charging infrastructure planning and operational phases to avoid drops in service quality.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Electrical Resistivity Changes During Heating Experiments Unravel Heterogeneous Thermal‐Hydrological‐Mechanical Processes in Salt Formations

Abstract Rock salt is considered a suitable medium for the permanent disposal of heat‐generating radioactive waste due to its isolation properties. However, excavation damage and heating induce complex and heterogeneous thermal‐hydrological‐mechanical (THM) processes across different zones. Quantifying this heterogeneity is crucial for accurate long‐term performance assessment models, but traditional methods lack the necessary resolution. This study employs 4D electrical resistivity tomography (ERT) monitoring during controlled heating experiments in a salt formation to unravel the spatiotemporal dynamics of THM processes. Advanced time‐lapse inversion and clustering analysis quantify subsurface properties and map the heterogeneity of THM dynamics. The ERT results can estimate subsurface properties and delineate the damaged and intact zones, enabling appropriate parameterization and representation of processes for long‐term modeling. This approach may be used in further improving the predictive models and ensuring the safe long‐term disposal of radioactive waste in rock salt.

58 GEOSCIENCES↗
DOE OSTI2026

AGFormer: Adaptive Spatiotemporal graph informed transformer for multi-reservoir inflow forecasting

Accurate reservoir inflow forecasting is crucial for effective water resource management, yet most machine learning models focus on single-reservoir prediction and overlook spatial dependencies among hydrologically connected reservoirs. Here, we propose AGFormer (Adaptive Graph-Informed Transformer), an end-to-end framework that integrates adaptive graph learning with temporal sequence modeling for multi-reservoir inflow forecasting. A shared encoder and graph attention mechanism generate reservoir-specific embeddings, which are then processed by the Transformer-based encoder–decoder for multi-step inflow forecasting. We also introduce a pretraining paradigm to learn robust temporal embeddings from misaligned historical records. Evaluated on 30 reservoirs in the Upper Colorado River Basin, AGFormer achieves superior seven-day-ahead forecasts, with NSE > 0.75 for 20 reservoirs—outperforming Encoder–Decoder LSTM, GCN+LSTM, and Transformer baselines. Adaptive graph learning captures dynamic inter-reservoir dependencies, and feature attribution aligns with snowmelt-driven hydrology. Incorporating forecasted meteorological inputs further enhances accuracy, demonstrating AGFormer’s potential to support reservoir management under dynamic hydrological conditions.

Adaptive graph learning↗
NASA NTRS2014

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data: Part III. Using Combined PCA to Compare Spatiotemporal Variability of MODIS, MISR and OMI Aerosol Optical Depth

Satellite measurements of global aerosol properties are very useful in constraining aerosol parameterization in climate models. The reliability of different data sets in representing global and regional aerosol variability becomes an essential question. In this study, we present the results of a comparison using combined principal component analysis (CPCA), applied to monthly mean, mapped (Level 3) aerosol optical depth (AOD) product from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Ozone Monitoring Instrument (OMI). This technique effectively finds the common space-time variability in the multiple data sets by decomposing the combined AOD field. The results suggest that all of the sensors capture the globally important aerosol regimes, including dust, biomass burning, pollution, and mixed aerosol types. Nonetheless, differences are also noted. Specifically, compared with MISR and OMI, MODIS variability is significantly higher over South America, India, and the Sahel. MODIS deep blue AOD has a lower seasonal variability in North Africa, accompanied by a decreasing trend that is not found in either MISR or OMI AOD data. The narrow swath of MISR results in an underestimation of dust variability over the Taklamakan Desert. The MISR AOD data also exhibit overall lower variability in South America and the Sahel. OMI does not capture the Russian wild fire in 2010 nor the phase shift in biomass burning over East South America compared to Central South America, likely due to cloud contamination and the OMI row anomaly. OMI also indicates a much stronger (boreal) winter peak in South Africa compared with MODIS and MISR.

satellite observation↗
DOE OSTI2024

Deep reinforcement learning for optimal control of induction welding process

Optimizing induction welding (IW) process parameters for the application of joining thermoplastic composites is challenging as it requires achieving complex spatiotemporal thermal characteristics along the weld-line to obtain desired weld quality. We formulate an optimal control problem which captures these requirements and seeks to optimize the IW coil speed using a fast-acting dynamic IW process model. We develop a novel Deep Reinforcement Learning (DRL) framework to solve this computationally challenging control problem and demonstrate via simulation study that the learned DRL feedback control policy results in better spatiotemporal thermal characteristics as compared to the current state-of-the-art.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Carbon Monitoring System Flux Estimation and Attribution: Impact of ACOS-GOSAT X(CO2) Sampling on the Inference of Terrestrial Biospheric Sources and Sinks

Using an Observing System Simulation Experiment (OSSE), we investigate the impact of JAXA Greenhouse gases Observing SATellite 'IBUKI' (GOSAT) sampling on the estimation of terrestrial biospheric flux with the NASA Carbon Monitoring System Flux (CMS-Flux) estimation and attribution strategy. The simulated observations in the OSSE use the actual column carbon dioxide (X(CO2)) b2.9 retrieval sensitivity and quality control for the year 2010 processed through the Atmospheric CO2 Observations from Space algorithm. CMS-Flux is a variational inversion system that uses the GEOS-Chem forward and adjoint model forced by a suite of observationally constrained fluxes from ocean, land and anthropogenic models. We investigate the impact of GOSAT sampling on flux estimation in two aspects: 1) random error uncertainty reduction and 2) the global and regional bias in posterior flux resulted from the spatiotemporally biased GOSAT sampling. Based on Monte Carlo calculations, we find that global average flux uncertainty reduction ranges from 25% in September to 60% in July. When aggregated to the 11 land regions designated by the phase 3 of the Atmospheric Tracer Transport Model Intercomparison Project, the annual mean uncertainty reduction ranges from 10% over North American boreal to 38% over South American temperate, which is driven by observational coverage and the magnitude of prior flux uncertainty. The uncertainty reduction over the South American tropical region is 30%, even with sparse observation coverage. We show that this reduction results from the large prior flux uncertainty and the impact of non-local observations. Given the assumed prior error statistics, the degree of freedom for signal is approx.1132 for 1-yr of the 74 055 GOSAT X(CO2) observations, which indicates that GOSAT provides approx.1132 independent pieces of information about surface fluxes. We quantify the impact of GOSAT's spatiotemporally sampling on the posterior flux, and find that a 0.7 gigatons of carbon bias in the global annual posterior flux resulted from the seasonally and diurnally biased sampling when using a diagonal prior flux error covariance.

biased sampling↗
DOE OSTI2025

Using feature importance as an exploratory data analysis tool on Earth system models

Abstract. Machine learning (ML) models are commonly used to generate predictions, but these models can also support the discovery of new science. Generating accurate predictions necessitates that a model captures the structure of the underlying data. If the structure is properly extracted, ML could be a useful exploratory and evidential tool. In this paper, we present a case study that demonstrates the use of ML for exploratory data analysis (EDA) in the climate space. We apply the ML explainability method of spatiotemporal zeroed feature importance (stZFI) to understand how climate-variable associations evolve over space and time. Our analyses focus on data from ensembles of Earth system models (ESMs) which provide data on different climate states and conditions. We elect to work with ESM ensembles since they allow us to compare feature importance across alternative scenarios not available with observed data. The ensembles also account for natural variability so that we can distinguish between signal and noise due to natural climate variability when computing feature importance. The use of perturbed initial condition ensembles introduces variability mimicking the natural variability in the atmosphere; thus the signals emerging using feature importance (FI) can be evaluated against the natural variability in the climate system. For our analyses, we consider the 1991 volcanic eruption of Mount Pinatubo, which was a large stratospheric aerosol injection. We explore the climate pathway associated with the eruption from aerosols to radiation to temperature at both the near-surface and stratospheric levels. In addition to applying the method to data generated from two different ESMs, we apply stZFI to reanalysis data to compare the associations identified by stZFI. We show how stZFI tracks the importance of aerosol optical depth over time on forecasting temperatures. This case study illustrates usefulness of an ML tool (stZFI) for EDA on a well-studied climate exemplar.

Ries, Daniel (ORCID:0000000250294647)↗
DOE OSTI2025

A Multivariate Space‐Time Dynamic Model for Characterizing the Atmospheric Impacts Following the Mt. Pinatubo Eruption

The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.

Dynamic Linear Model↗
DOE OSTI2024

Heterogeneity in Permeability and Particulate Organic Carbon Content Controls the Redox Condition of Riverbed Sediments at Different Timescales

Abstract The hydrological and biogeochemical properties of the hyporheic zone in stream and riverine ecosystems have been extensively studied over the past two decades. Although it is widely acknowledged that sediment heterogeneity can influence biogeochemical reactions, little effort has been made to understand the role of heterogeneity on the spatiotemporal variability of riverbed redox conditions under changing flow dynamics at different timescales. Here we integrate a mechanistic model and field data to demonstrate that heterogeneity in permeability plays a vital role in modulating sediment redox conditions at both seasonal (annual) and event (daily‐to‐weekly) timescales, whereas heterogeneity in particulate organic carbon (POC) content only has a comparable influence on redox conditions at the seasonal timescale. These findings underscore the importance of accurately characterizing sediment heterogeneity, in terms of permeability and POC content, in quantifying biogeochemical dynamics in the riverbed and hyporheic zones of riverine ecosystems.

Geology↗
DOE OSTI2024

UFNet: Joint U-Net and Fully Connected Neural Network to Bias Correct Precipitation Predictions from

Paper information. Shuang Yu, Indrasis Chakraborty, Gemma J. Anderson, Donald D. Lucas, Yannic Lops, and Daniel Galea. UFNet: Joint U-Net and fully connected neural network to bias correct precipitation predictions from climate models. Artificial Intelligence for the Earth Systems, 2024. Overview. This work develops the UFNet methodology to correct E3SM historical precipitation projection bias. The UFNet deep learning framework consists of a two-part architecture: a U-Net convolutional network to capture the spatiotemporal distribution of precipitation and a fully connected network to capture the distribution of higher-order statistics. The joint network, termed UFNet, can simultaneously improve the spatial structure of the modeled precipitation and capture the distribution of extreme precipitation values. Below we provide guidance for applying UFNet to correct the Energy Exascale Earth System Model (E3SM; Golaz et al. 2019) daily precipitation projection over the contiguous United States (CONUS). Getting started 1. Obtain the historical climate simulation and observation data. The E3SM historical simulation data are available through https://aims2.llnl.gov/search/cmip6/. The CPC unified gauge-based analysis of daily precipitation can be found through https://psl.noaa.gov/data/gridded/data.cpc.globalprecip.html. The ECMWF atmospheric reanalysis of the 20th century (ERA-20C) data are available through https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era-20c. The spatial resolution of E3SM and observed datasets are both regridded to a common 1° resolution grid using conservative interpolation. The regridded E3SM, CPC and ERA-20C with 1° resolution can be found throught ./data/. 2. Train the fully connected network (DNN) Python train_dnn.py 3. Train the UFNet Python train_ufnet.py 4. Evaluation and compared with the baseline Python evaluation.py

Lucas, DonaldD↗
NASA NTRSDate not supplied

Characterizing Spatiotemporal Uncertainty in Interpolated Meteorological Data

Interpolated meteorological data invariably contain errors. These errors have structure in time and space, particularly autocorrelation, which can cause the effects of errors to compound when model outputs are aggregated temporally or spatially. One way to account for this uncertainty is with a probabilistic model from which samples can be drawn that are coherent with respect to underlying spatial and temporal covariance structure. This work describes a probabilistic method for spatial interpolation of point-wise meteorological time series. Observational data from weather stations are generally sparse in space and dense in time (but sometimes missing). The method works by projecting time series onto orthogonal basis vectors and spatially interpolating each resulting component independently. Under suitable assumptions, and data transformations to better satisfy those assumptions, Gaussian process regression provides a complete description of the joint predictive distribution over a Gaussian random field. Spatiotemporally coherent realizations are generated as the sum of conditional (spatial) simulations of each orthogonal (temporal) component. Data-derived and generic orthogonal bases are considered. In addition to spatial interpolation, imputation of missing observational data is examined. The method is applied using near-surface air temperature over the Western United States and validated by comparing theoretical versus actual coverage of predictive distributions and analyzing the degree to which spatial and temporal covariance structure is reproduced. Computational considerations, relating to conditional simulation of random fields, are also addressed.

Conor T Doherty↗