Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AUTOCORRELATION”

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 901 records · Page 50

Dust Abrasion Damage on Martian Solar Arrays: Experimental Investigation and Opportunity Rover Performance Analysis

Here we investigate the effects of erosion and weathering that occur on III-V cover-glass interconnected cells (CICs) after exposure to Mars dust storm conditions. The durability of these materials in a Martian environment is not well characterized so we perform analogous experimentation. To replicate the dust impingement, test coupons were placed in an enclosure and sandblasted with Mars dust simulant. We show the J-V response dependency on both incident angle and exposure times. We find that the simulated Martian dust storm often results in damage to the anti-reflective coating and subsequent reduced short circuit current. Reduction in the open circuit voltage is observed, likely caused by structural damage to the crystal lattice after CIC fracture. We employ data-driven modeling to determine a performance degradation rate that is consistent with zero within uncertainty. We also quantify the soiling contribution and power degradation of the photovoltaic cells on Mars through analysis of 4.95 Martian years of report-out power conditions from the Opportunity rover. We find that atmospheric dust suspended due to a weather event does not result in instantaneous settled dust on the PV cells. We calculate via autocorrelation function that the dust settling rate is approximately 21 Sols from atmospheric dust suspension. The findings presented here deliver a realistic approximation for the insolation values and subsequent PV power expected over time on the Martian surface thus informing future dust abatement systems.

McMillon-Brown, Lyndsey↗

Dust Abrasion Damage on Martian Solar Arrays: Experimental Investigation and Opportunity Rover Performance Analysis

Here we investigate the effects of erosion and weathering that occur on III-V cover-glass interconnected cells (CICs) after exposure to Mars dust storm conditions. The durability of these materials in a Martian environment is not well characterized so we perform analogous experimentation. To replicate the dust impingement, test coupons were placed in an enclosure and sandblasted with Mars dust simulant. We show the J-V response dependency on both incident angle and exposure times. We find that the simulated Martian dust storm often results in damage to the anti-reflective coating and subsequent reduced short circuit current. Reduction in the open circuit voltage is observed, likely caused by structural damage to the crystal lattice after CIC fracture. We employ data-driven modeling to determine a performance degradation rate that is consistent with zero within uncertainty. We also quantify the soiling contribution and power degradation of the photovoltaic cells on Mars through analysis of 4.95 Martian years of report-out power conditions from the Opportunity rover. We find that atmospheric dust suspended due to a weather event does not result in instantaneous settled dust on the PV cells. We calculate via autocorrelation function that the dust settling rate is approximately 21 Sols from atmospheric dust suspension. The findings presented here deliver a realistic approximation for the insolation values and subsequent PV power expected over time on the Martian surface thus informing future dust abatement systems.

McMillon-Brown, Lyndsey↗

Evaluation and Validation of a High Spatial Resolution Satellite Soil Moisture Product over the Continental United States

The soil moisture (SM) data retrieved from the Soil Moisture Active and Passive (SMAP) satellite are available at a 9 km grid spacing since April 2015. This product can provide valuable information for research and applications in hydrology and other related fields. However, the resolution may be too coarse for applications at catchment or field scale. In this study, an established downscaling methodology, which had a major modification regarding its application on the SMAP 33 km domain, was implemented to develop a 1 km soil moisture product based on the SMAP 9 km data. The algorithm proposed here is based on the thermal inertia principle and developed by modeling the relationship between surface temperature difference and SM for different Normalized Difference Vegetation Index (NDVI) classes. The model functions were established and tuned using data from the NASA’s Land Information System (LIS) North America Land Data Assimilation System (NLDAS) and remotely sensed VISible/InfRared (VIS/IR) reflectance data from Long Term Data Record (LTDR) AVHRR (Advanced Very High Resolution Radiometer) for the growing season months of April-September 1981–2018. These were then implemented using the MODIS (Moderate Resolution Imaging Spectroradiometer) data over the Continental United States (CONUS) domain. Validation activities were carried out using in situ measurements distributed through the International Soil Moisture Network (ISMN). The validation results computed using the 1 km SM data showed that the R2, unbiased RMSE (root mean square error) and bias were improved relative to the 9 km SMAP product by 0.045, 0.018m3/m3 and 0.001m3/m3, respectively. The 1 km SM also exhibited a strong time-series autocorrelation. Further accuracy assessment analyses indicated that precipitation might contribute to the uncertainties in both the 9 km SMAP and 1 km downscaled SMAP SM products.

SMAP↗

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils from Combined Active-Passive Microwave SMAP Observations

An active-passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height (s) and horizontal correlation length (l)] is presented for bare soils. The algorithm is based on active-passive microwave covariation, including the improved Integral Equation Method (I²EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for εₛ > 10 [-]. Furthermore, the physical model basis of this approach (I²EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

correlation length↗

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils From Combined Active–Passive Microwave SMAP Observations

An active–passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height ( s ) and horizontal correlation length ( l )] is presented for bare soils. The algorithm is based on active–passive microwave covariation, including the improved Integral Equation Method (I(exp 2) EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for E(sub s) > 10 [-]. Furthermore, the physical model basis of this approach (I(exp 2) EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

Correlation length↗

Improved SIFTER v2 algorithm for long-term GOME-2A satellite retrievals of fluorescence with a correction for instrument degradation

Solar-induced fluorescence (SIF) data from satellites are increasingly used as a proxy for photosynthetic activity by vegetation, and as a constraint on gross primary production. Here we report on improvements in the algorithm to retrieve mid-morning (09:30 hrs local time) SIF estimates on the global scale from GOME-2 sensor on the Metop-A satellite (GOME-2A) for the period 2007-2019. Our new SIFTER (Sun-Induced Fluorescence of Terrestrial Ecosystems Retrieval) v2 algorithm improves over a previous version by using a narrower spectral window that avoids strong oxygen absorption and is less sensitive to water vapour absorption, by constructing stable reference spectra from a 6-year period (2007-2012) of atmospheric spectra over the Sahara, and by applying a latitude-dependent zero-level adjustment that accounts for biases in the data product. We generated stable, good-quality SIF retrievals between January 2007 and June 2013, when GOME-2A degradation in the near infrared was still limited. After the narrowing of the GOME-2A swath in July 2013, we characterized the throughput degradation of the level-1 data in order to derive reflectance corrections and apply these for the SIF retrievals between July 2013 and December 2018. SIFTER v2 data compares well with the independent NASA v2.8 data product. Especially in the evergreen tropics, SIFTER v2 no longer shows the underestimates against other satellite products that were seen in SIFTER v1. The new data product includes uncertainty estimates for individual observations, and is best used for mostly clear-sky scenes, and when spectral residuals remain below a certain spectral autocorrelation threshold. Our results support the use of SIFTER v2 data to be used as an independent constraint on photosynthetic activity on regional to global scales.

Solar-induced fluorescence (SIF)↗

Drought Relief and Reversal Over North America from 1500 to 2016

Season-to-season persistence of soil moisture drought varies across North America. Such inter-seasonal autocorrelation can have modest skill in forecasting future conditions several months in advance. Because robust instrumental observations of precipitation span less than 100 years, the temporal stability of the relationship between seasonal moisture anomalies is uncertain. The North American Seasonal Precipitation Atlas (NASPA) is a gridded network of separately reconstructed cool (December-April) and warm (May-July) season precipitation series and offers new insights on the intra-annual changes in drought for up to 2,000 years. Here, the NASPA precipitation reconstructions are rescaled to represent the long-term soil moisture balance during the cool season and three-month long atmospheric moisture during the warm season. These rescaled seasonal reconstructions are then used to quantify the frequency, magnitude, and spatial extent of cool season drought that was relieved or reversed during the following summer months. The adjusted seasonal reconstructions reproduce the general patterns of large-scale drought amelioration and termination in the instrumental record during the 20th century and are used to estimate relief and reversals for the most skillfully reconstructed past 500-years. Sub-continental to continental scale reversals of cool season drought in the following warm season have been rare, but the reconstructions display periods prior to the instrumental data of increased reversal probabilities for the Mid-Atlantic region and the US Southwest. Drought relief at the continental scale may arise in part from macroscale ocean-atmosphere processes, while the smaller scale regional reversals may reflect land-surface feedbacks and stochastic variability.

Drought↗

Spatial validation reveals poor predictive performance of large-scale ecological mapping models

Mapping aboveground forest biomass is central for assessing the global carbon balance. However, current large-scale maps show strong disparities, despite good validation statistics of their underlying models. Here, we attribute this contradiction to a flaw in the validation methods, which ignore spatial autocorrelation (SAC) in data, leading to overoptimistic assessment of model predictive power. To illustrate this issue, we reproduce the approach of large-scale mapping studies using a massive forest inventory dataset of 11.8 million trees in central Africa to train and validate a random forest model based on multispectral and environmental variables. A standard nonspatial validation method suggests that the model predicts more than half of the forest biomass variation, while spatial validation methods accounting for SAC reveal quasi-null predictive power. This study underscores how a common practice in big data mapping studies shows an apparent high predictive power, even when predictors have poor relationships with the ecological variable of interest, thus possibly leading to erroneous maps and interpretations.

biomass↗

Leveraging the ALMA Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide

We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations from the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH3OH, H2CO, and HNC (along with continuum emission from dust) probed the inner coma of C/2015 ER61, revealing asymmetric outgassing and discerning parent from daughter/distributed source species. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for each molecule and a stringent upper limit for CO. HCN is consistent with direct nucleus release in C/2015 ER61, whereas CS, H2CO, HNC, and potentially CH3OH are associated with distributed sources in the coma. Adopt- ing a Haser model, parent scale lengths determined for H2CO (Lp ∼ 2200 km) and HNC (Lp ∼ 3300 km) are consistent with previous work in comets, whereas significant extended source production (Lp ∼ 2000 km) is indicated for CS, suggesting production from an unknown parent in the coma. The con- tinuum presents a point-source distribution with a flux density implying an excessively large nucleus, inconsistent with other estimates of the nucleus size. It is best explained by the thermal emission of slowly moving outburst ejectas, with total mass 5–8 × 1010 kg. These results demonstrate the power of the ACA for revealing the abundances, spatial distributions, and locations of molecular production for volatiles in moderately bright comets such as C/2015 ER61.

Molecular spectroscopy (2095) — High resolution sp↗

Leveraging the Atacama Compact Array for Cometary Science: An Interferometric Survey of Comet C/2015 ER61 (PanSTARRS) and Evidence for a Distributed Source of Carbon Monosulfide

We report the first survey of molecular emission from cometary volatiles using standalone Atacama Compact Array (ACA) observations from the Atacama Large Millimeter/Submillimeter Array (ALMA) toward comet C/2015 ER61 (PanSTARRS) carried out on UT 2017 April 11 and 15, shortly after its April 4 outburst. These measurements of HCN, CS, CH3OH, H2CO, and HNC (along with continuum emission from dust) probed the inner coma of C/2015 ER61, revealing asymmetric outgassing and discerning parent from daughter/distributed source species. This work presents spectrally integrated flux maps, autocorrelation spectra, production rates, and parent scale lengths for each molecule and a stringent upper limit for CO. HCN is consistent with direct nucleus release in C/2015 ER61, whereas CS, H2CO, HNC, and potentially CH3OH are associated with distributed sources in the coma. Adopt- ing a Haser model, parent scale lengths determined for H2CO (Lp ∼ 2200 km) and HNC (Lp ∼ 3300 km) are consistent with previous work in comets, whereas significant extended source production (Lp ∼ 2000 km) is indicated for CS, suggesting production from an unknown parent in the coma. The continuum presents a point-source distribution with a flux density implying an excessively large nucleus, inconsistent with other estimates of the nucleus size. It is best explained by the thermal emission of slowly moving outburst ejectas, with total mass 5–8 × 1010 kg. These results demonstrate the power of the ACA for revealing the abundances, spatial distributions, and locations of molecular production for volatiles in moderately bright comets such as C/2015 ER61.

Molecular spectroscopy↗

Weekly Mapping of Sea Ice Freeboard in the Ross Sea from ICESat-2

NASA’s ICESat-2 has been providing sea ice freeboard measurements across the polar regions since October 2018. In spite of the outstanding spatial resolution and precision of ICESat-2, the spatial sparsity of the data can be a critical issue for sea ice monitoring. This study employs a geostatistical approach (i.e., ordinary kriging) to characterize the spatial autocorrelation of the ICESat-2 freeboard measurements (ATL10) to estimate weekly freeboard variations in 2019 for the entire Ross Sea area, including where ICESat-2 tracks are not directly available. Three variogram models (exponential, Gaussian, and spherical) are compared in this study. According to the cross-validation results, the kriging-estimated freeboards show correlation coefficients of 0.56–0.57, root mean square error (RMSE) of ~0.12 m, and mean absolute error (MAE) of ~0.07 m with the actual ATL10 freeboard measurements. In addition, the estimated errors of the kriging interpolation are low in autumn and high in winter to spring, and low in southern regions and high in northern regions of the Ross Sea. The effective ranges of the variograms are 5–10 km and the results from the three variogram models do not show significant differences with each other. The southwest (SW) sector of the Ross Sea shows low and consistent freeboard over the entire year because of the frequent opening of wide polynya areas generating new ice in this sector. However, the southeast (SE) sector shows large variations in freeboard, which demonstrates the advection of thick multiyear ice from the Amundsen Sea into the Ross Sea. Thus, this kriging-based interpolation of ICESat-2 freeboard can be used in the future to estimate accurate sea ice production over the Ross Sea by incorporating other remote sensing data.

Satellite altimeter↗

Markovian Statistical Model of Cloud Optical Thickness. Part I: Theory and Examples

We present a generalization of the binary-value Markovian model previously used for statistical characterization of cloud masks to a continuous-value model describing 1D fields of cloud optical thickness (COT). This model has simple functional expressions and is specified by four parameters: the cloud fraction, the autocorrelation (scale) length, and the two parameters of the normalized probability density function of (non-zero) COT values (this PDF is assumed to have gamma-distribution form). Cloud masks derived from this model by separation between the values above and below some threshold in COT appear to have the same statistical properties as in binary-value model described in our previous publications. We demonstrate the ability of our model to generate examples of various cloud-field types by using it to statistically imitate actual cloud observations made by the Research Scanning Polarimeter (RSP) during two field experiments.

binary-value Markovian model↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗

Insights Into the Hydrology of the Congo Peatlands Through Land Surface Modeling and Data Assimilation

The 16.8 million ha of peatlands in the Cuvette Centrale wetland complex in the Congo Basin is one of the largest peatland regions on Earth but still highly understudied. Understanding the hydrological functioning of these peatlands and the effects of external disturbances thereon remains a major challenge. Recent research suggested fundamental hydrological differences between the Congo peatlands and the well-studied Southeast Asian peatlands. The Congo peatlands have a doming gradient that is up to ten times smaller, and they are influenced by river hydrology to some extent. In this study, we explore the Congo peatland hydrology through land surface modeling and data assimilation. We build upon our recently developed tropical PEATCLSM module (Apers et al., 2022) that was parameterized based on data from Southeast Asian peatlands due to the lack of field data from other tropical peatland regions. In a first step, we derive Congo-specific peat hydraulic and discharge function parameters from a scalar parametrization of water level dynamics in the Congo peatlands, using observed water level data at two locations. These Congo-specific parameters differ considerably from the original literature-based parameters from Southeast Asian peatlands. In a second step, we apply our original and Congo-specific parameters in an assimilation scheme for L-band brightness temperature (Tb) data from the Soil Moisture and Ocean Salinity (SMOS) mission. The data assimilation results are used in two ways. First, the effect of these parameters on the simulated peatland hydrology and the observation-minus-forecast Tb residuals is evaluated. It is hypothesized that the new parameters reduce the previously reported modeling errors over the Congo peatlands and reduce the residuals in Tb as well. Second, we analyze the data assimilation diagnostics to learn about other model improvement possibilities. In preliminary results, we found long periods of temporally autocorrelated total water storage increments (difference of pre- and post-update) that coincided with anomalies in river stages measured upstream of the peatlands. Since PEATLCSM neglects possible river influence, this concurrence suggests that the typically used grid-based approach of land surface models should be combined with a river routing scheme over the Congo peatlands.

Sebastian Apers↗

A Spaceborne Demonstration of P-Band Signals-of-Opportunity (SoOp) Reflectometry

Land-reflected signals from a geosynchronous communication satellite broadcasting in P-band (367.5 MHz) were captured in low Earth orbit (LEO) using a simple dipole antenna. A delay-Doppler map (DDM) was generated through autocorrelation. Estimates of the specular point delay were obtained from the lag of the second peak in the DDM with a bias of 239.4 m and a standard deviation of 44 m (12 m over a frozen lake) with respect to a predicted orbit model. Relative magnitudes of the first and second DDM peaks fell within the range of values predicted using dielectric models for the frozen ground and lake. Finally, retrievals of surface reflection coefficient were generated using a range of realistic values for the transmitter link budget G/T , and these also fell within the range of possible values for the antenna gain pattern. Given the lack of calibration and the large uncertainties in the receiver orbit and attitude, this agreement is sufficient to conclude a successful demonstration of the fundamental principle of single-antenna reflectometry in P-band. P-band reflectometry may offer a new approach to remote sensing of sub-canopy and root-zone soil moisture.

Delays↗

Global Intercomparison of hyper-resolution ECOSTRESS coastal sea surface temperature measurements from the Space Station with VIIRS-N20

The ECOSTRESS multi-channel thermal radiometer on the Space Station has an unprecedented spatial resolution of 70 m and a return time of hours to 5 days. It resolves details of oceanographic features not detectable in imagery from MODIS or VIIRS, and has open-ocean coverage, unlike Landsat. We calibrated two years of ECOSTRESS sea surface temperature observations with L2 data from VIIRS-N20 (2019–2020) worldwide but especially focused on important upwelling systems currently undergoing climate change forcing. Unlike operational SST products from VIIRS-N20, the ECOSTRESS surface temperature algorithm does not use a regression approach to determine temperature, but solves a set of simultaneous equations based on first principles for both surface temperature and emissivity. We compared ECOSTRESS ocean temperatures to well-calibrated clear sky satellite measurements from VIIRS-N20. Data comparisons were constrained to those within 90 min of one another using co-located clear sky VIIRS and ECOSTRESS pixels. ECOSTRESS ocean temperatures have a consistent 1.01 ◦C negative bias relative to VIIRS-N20, although deviation in brightness temperatures within the 10.49 and 12.01 µm bands were much smaller. As an alternative, we compared the performance of NOAA, NASA, and U.S. Navy operational split-window SST regression algorithms taking into consideration the statistical limitations imposed by intrinsic SST spatial autocorrelation and applying corrections on brightness temperatures. We conclude that standard bias-correction methods using already validated and well-known algorithms can be applied to ECOSTRESS SST data, yielding highly accurate products of ultra-high spatial resolution for studies of biological and physical oceanography in a time when these are needed to properly evaluate regional and even local impacts of climate change.

Nicolas Weidberg↗

Anisotropy of Density Fluctuations in the Solar Wind at 1 au

A well-known property of solar wind plasma turbulence is the observed anisotropy of the autocorrelations, or equivalently the spectra, of velocity and magnetic field fluctuations. Here we explore the related but apparently not well-studied issue of the anisotropy of plasma density fluctuations in the energy-containing and inertial ranges of solar wind turbulence. Using 10 yr (1998–2008) of in situ data from the Advanced Composition Explorer mission, we find that for all but the fastest wind category, the density correlation scale is slightly larger in directions quasi-parallel to the large-scale mean magnetic field as compared to quasi-perpendicular directions. The correlation scale in fast wind is consistent with isotropic. The anisotropy as a function of the level of correlation is also explored. We find at small correlation levels, i.e., at energy-containing scales and larger, the density fluctuations are close to isotropy for fast wind, and slightly favor more rapid decorrelation in perpendicular directions for slow and medium winds. At relatively smaller (inertial range) scales where the correlation values are larger, the sense of anisotropy is reversed in all speed ranges, implying a more "slablike" structure, especially prominent in the fast wind samples. We contrast this finding with published results on velocity and magnetic field correlations.

Interplanetary turbulence↗

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↗