Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

The LANDSAT-1 multispectral scanner as a tool in the classification of inland lakes

Relationships between LANDSAT-1 multispectral scanner (MSS) data and the trophic status of a group of lakes in the north-northeastern part of the United States were studied by predicting the magnitudes of two trophic state indicators, estimating lake position on a multivariate trophic scale, and automatically classifying lakes according to their trophic state. Initially, the principal component ordination was employed with 100 lakes. MSS data for some 20 lakes was then extracted from computer-compatible tapes (CCT) using a binary marking technique. The output was in the form of descriptive statistics and photographic concatenations. Color ratios were incorporated into regression models for the prediction of Secchi disc transparency, chlorophyll a, and lake position on the tropic scale. Results indicate that the LANDSAT-1 system, although handicapped by low spectral and spatial resolutions as well as excessive cloud cover, can be used as a supplemental data source in lake survey programs.

Boland, D. H. P.↗
NASA NTRS1994

Statistical study of ULF wave occurrence in the dayside magnetosphere

Ultralow-frequency (ULF) waves are observed almost everywhere in the dayside magnetosphere. The mechanism by which these waves are generated and transformed in the dayside magnetosphere is still not understood. Here we report a statistical study of these waves based on magnetic field data from the International Sun-Earth Explorer 1 (ISEE 1) spacecraft. Data from the first traversal of the spacecraft through the entire dayside magnetosphere have been examined to determine the spatial distribution of wave occurrence. Successive 20-min segments of data were transformed to a field-aligned coordinate system. The parallel component was detrended and all three components of the field spectrally analyzed. Wave occurrence was defined by the presence of significant peaks in the power spectra. Wave events were categorized by three wave frequency bands: Pc 3 with T approximately 10-45 s; Pc 4 with T approximately 45-150 s; the short-period part of the Pc 5 wave band with T approximately 150-324 s. Properties of the spectral peaks were then entered into a data base. The data base was next sorted to determine the spatial occurrence pattern for the waves. Our results show that Pc 3 waves most frequently occur just outside synchronous orbit and are approximately centered on local noon. Pc 4 waves have a similar distribution with its peak further out. Pc 5 waves have high occurrence rate at the two flanks of the magnetosphere. Peaks in spectra obtained near the magnetopause are less clearly defined than those deeper in the magnetosphere.

Cao, M.↗
NASA NTRS2004

The Ultra-Luminous X-ray Source Population from the Chandra Archive of Galaxies

One hundred fifty-four discrete non-nuclear Ultra-Luminous X-ray (ULX) sources, with spectroscopically-determined intrinsic X-ray luminosities greater than 1 e39 ergs/s, are identified in 82 galaxies observed with Chandra's Advanced CCD Imaging Spectrometer. Source positions, X-ray luminosities, and spectral and timing characteristics are tabulated. Statistical comparisons between these X-ray properties and those of the weaker discrete sources in the same fields (mainly neutron star and stellar-mass black hole binaries) are made. Sources above approximately le38 ergs per second display similar spatial, spectral, color, and variability distributions. In particular, there is no compelling evidence in the sample for a new and distinct class of X-ray object such as the intermediate-mass black holes. 83% of ULX candidates have spectra that can be described as absorbed power laws with index = 1.74 and column density = 2.24e21 l per square centimeter, or approximately 5 times the average Galactic column. About 20% of the ULX's have much steeper indices indicative of a soft, and likely thermal, spectrum. The locations of ULXs in their host galaxies are strongly peaked towards their galaxy centers. The deprojected radial distribution of the ULX candidates is somewhat steeper than an exponential disk, indistinguishable from that of the weaker sources. About 5--15% of ULX candidates are variable during the Chandra observations (which average 39.5 ks). Comparison of the cumulative X-ray luminosity functions of the ULXs to Chandra Deep Field results suggests approximately 25% of the sources may be background objects including 14% of the ULX candidates in the sample of spiral galaxies and 44% of those in elliptical galaxies implying the elliptical galaxy ULX population is severely compromised by background active galactic nuclei. Correlations with host galaxy properties confirm the number and total X-ray luminosity of the ULXs are associated with recent star formation and with galaxy merging and interactions. The preponderance of ULXs in star-forming galaxies as well as their similarities to less-luminous sources suggest they originate in a young but short-lived population such as the high-mass X-ray binaries with a smaller contribution (based on spectral slope) from recent supernovae. The number of ULXs in elliptical galaxies scales with host galaxy mass and can be explained most simply as the high-luminosity end of the low-mass X-ray binary population.

Swartz, Douglas A.↗
NASA NTRS2014

Prediction of Acoustic Environments from Horizontal Rocket Firings

In recent years, advances in research and engineering have led to more powerful launch vehicles which can reach areas of space not yet explored. These more powerful vehicles yield acoustic environments potentially destructive to the vehicle or surrounding structures. Therefore, it has become increasingly important to be able to predict the acoustic environments created by these vehicles in order to avoid structural and/or competent failure. The current industry standard technique for predicting launch-induced acoustic environments was developed by Eldred in the early 1970's and is published in NASA SP-80721. Recent work2 has shown Eldred's technique to be inaccurate for current state-of-the-art launch vehicles. Due to the high cost of full-scale and even sub-scale rocket experiments, very little rocket noise data is available. Furthermore, much of the work thought to be applicable to rocket noise has been done with heated jets. Tam3,4 has done an extensive amount of research on jets of different nozzle exit shape, diameter, velocity, and temperature. Though the values of these parameters, especially exit velocity and temperature, are often very low compared to these values in rockets, a lot can be learned about rocket noise from jet noise literature. The turbulent nature of jet and rocket exhausts is quite similar. Both exhausts contain turbulent structures of varying scale-termed the fine and large scale turbulence by Tam. The finescale turbulence is due to small eddies from the jet plume interacting with the ambient atmosphere. According to Tam et al., the noise radiated by this envelope of small-scale turbulence is statistically isotropic. Hence, one would expect the noise from the small scale turbulence of the jet to be nearly omni-directional. The coherent nature of the large-scale turbulence results in interference of the noise radiated from different spatial locations within the jet. This interference-whether it is constructive or destructive-results in highly directional noise radiation. Tam3 has proposed a model to predict the acoustic environment due to jets and while it works extremely well for jets, it was found to be inappropriate for rockets8. A model to predict the acoustic environment due to a launch vehicle in the far-field which incorporates concepts from both Eldred and Tam was created. This was done using five sets of horizontally fired rocket data, obtained between 2008 and 2012. Three of these rockets use solid propellant and two use liquid propellant. Through scaling analysis, it is shown that liquid and solid rocket motors exhibit similar spectra at similar amplitudes. This model is accurate for these five data sets within 5 dB of the measured data for receiver angles of 30deg to 160deg (with respect to the downstream exhaust centerline). The model uses the following vehicle parameters: nozzle exit diameter and velocity, radial distance from source to receiver, receiver angle, mass flow rate, and acoustic efficiency.

Giacomoni, Clothilde↗
NASA NTRS2015

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗
DOE OSTI2021

Application of electron backscatter diffraction techniques to quantify effects of aging on sub-grain and spatial heterogeneity in NMC cathodes

Identification and evaluation of structural heterogeneity in spent cathode materials is crucial to developing appropriate remediation strategies for novel recycling processes. Native heterogeneities may be exacerbated during the cell's operational lifetime, as sub-particle-scale variations induce anisotropic expansion and contraction upon cycling. Structural transformations resulting from repeated cycling and calendar aging predominantly occur at the secondary particle surface and at the grain boundaries (GBs) between primary particles. However, the diffusion and stress build up around and across GBs are poorly understood. In this study, electron backscatter diffraction (EBSD) is employed to track sub-grain lattice structure across a statistically relevant number of Li(Ni0.33Mn0.33Co0.33)O2 (NMC-111) particles. Specifically, differences in lattice misorientation – measured as the deviation from the grain's average orientation – are tracked as a function of position within the electrode (near current collector, middle, and near separator) and as a function of electrochemical cycling. Further, a novel method of structural analysis is developed, offering insight into sub-grain diffusion behavior by comparing lattice misorientation near the grain boundary versus in the grain bulk. The present results suggest that electrode-scale spatial heterogeneity in lattice structure is induced by initial manufacturing conditions, and that radial gradients in lattice misorientation evolve at the primary particle scale with repeated electrochemical cycling.

25 ENERGY STORAGE↗
NASA NTRS1986

A statistical study of plasma sheet dynamics using ISEE 1 and 2 energetic particle flux data

Plasma sheet dynamics during substorms are studied by analyzing 461 cases of transient dropout events of the 1.5 and 6-keV particle fluxes detected by ISEE 1 and 2 satellites. The instruments for detecting low- and high-energy particles are described. The spatial distribution of flux dropout events, and the events' relationship to magnetospheric activity level are examined. Substorm events without observed flux dropout events are investigated. The data reveal that the flux dropout distribution is isotropic, between 12-23 earth radii, and is present in the entire nightside plasma sheet; and the substorms without flux dropout are more frequent near earth and magnetospheric flanks. It is observed that tailward of 12 earth radii the flux dropout events and substorms without flux dropout are similar. The Chao et al. (1977) MHD rarefaction wave propagation model and the Hones (1973, 1980) near-tail, X-type magnetic neutral line formation model are discussed and compared to the experimental data. It is noted that neither model explains the plasma sheet dynamics observed.

Dandouras, J.↗
DOE OSTI2024

Imaging and structure analysis of ferroelectric domains, domain walls, and vortices by scanning electron diffraction

Direct electron detectors in scanning transmission electron microscopy give unprecedented possibilities for structure analysis at the nanoscale. In electronic and quantum materials, this new capability gives access to, for example, emergent chiral structures and symmetry-breaking distortions that underpin functional properties. Quantifying nanoscale structural features with statistical significance, however, is complicated by the subtleties of dynamic diffraction and coexisting contrast mechanisms, which often results in a low signal-to-noise ratio and the superposition of multiple signals that are challenging to deconvolute. Here we apply scanning electron diffraction to explore local polar distortions in the uniaxial ferroelectric Er(Mn,Ti)O 3 . Using a custom-designed convolutional autoencoder with bespoke regularization, we demonstrate that subtle variations in the scattering signatures of ferroelectric domains, domain walls, and vortex textures can readily be disentangled with statistical significance and separated from extrinsic contributions due to, e.g., variations in specimen thickness or bending. The work demonstrates a pathway to quantitatively measure symmetry-breaking distortions across large areas, mapping structural changes at interfaces and topological structures with nanoscale spatial resolution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Operando characterization and regulation of metal dissolution and redeposition dynamics near battery electrode surface

Mn dissolution has been a longstanding, ubiquitous issue that negatively impacts the performance of Mn-based battery materials. Mn dissolution involves complex chemical and structural transformations at the electrode-electrolyte interface. The continuously evolving electrodeelectrolyte interface has posed great challenges for characterizing the dynamic interfacial process and quantitatively establishing the correlation with battery performance. In this study, we visualize and quantify the temporally and spatially resolved Mn dissolution/redeposition (D/R) dynamics of electrochemically operating Mn-containing cathodes. The particle-level and electrode-level analyses reveal that the D/R dynamics is associated with distinct interfacial degradation mechanisms at different states of charge. Our results statistically differentiate the contributions of surface reconstruction and Jahn-Teller distortion to the Mn dissolution at different operating voltages. Introducing sulfonated polymers (Nafion) into composite electrodes can modulate the D/R dynamics through trapping the dissolved Mn species and rapidly establishing the local Mn D/R equilibrium. This work represents an inaugural effort to pinpoint the chemical and structural transformations responsible for Mn dissolution via an operando synchrotron study and develops an effective method to regulate Mn interfacial dynamics for improving battery performance.

25 ENERGY STORAGE↗
DOE OSTI2024

Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach

Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Data and scripts associated with “Moisture content modulates DOM thermodynamic regulation of oxygen consumption in drying streambed sediments”

This data package is associated with the publication “Moisture content modulates DOM thermodynamic regulation of oxygen consumption in drying streambed sediments” published in Scientific Reports (Garayburu-Caruso et al., 2026). The package contains processed data products and scripts used to quantify how drying and re-inundation of riverbed sediments influence dissolved organic matter (DOM) thermodynamic properties and their relationship with sediment oxygen (O₂) consumption across 33 stream sites in the contiguous United States. The data package contains DOM thermodynamic metrics (e.g., Gibbs free energy of carbon oxidation and thermodynamic efficiency), and O₂ consumption along with watershed-scale climate and land-cover metrics used as explanatory variables in the analyses. Underlying unprocessed and processed ultrahigh-resolution mass spectrometry data, oxygen consumption rates from laboratory moisture-manipulation experiments, within-sample environmental properties, sediment moisture content and contextual field measurements are archived separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2428003 (Laan et al., 2024) and https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1923689 (Forbes et al.,2023). A preliminary version of this data package was published in February 2026 at the time of manuscript submission. It was updated in June 2026, at the time of manuscript acceptance, to include the finalized data and additional metadata (readme, data dictionary, and file level metadata). For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. At the top level, the data package is organized into five main folders: (1) Data, (2)Figures, (3) Map, (4) GAM_Reulsts, and (5) src. The Data folder contains analysis-ready tabular files with oxygen consumption rates, DOM thermodynamic properties by site and treatment, site-level environmental variables, watershed-scale metrics, and other derived variables referenced in the manuscript. The Figures folder contains static image files associated with the main text and supplemental figures, while the Map folder includes spatial data and map-layer files used to create the sampling-location map. The GAM results folder contains the results for each of the general additive model (GAM).The src folder contains R scripts used to perform data processing, statistical analyses (including clustering, generalized additive models, and threshold analysis), and figure generation. This data package is associated with a GitHub repository found at https://github.com/WHONDRS-Hub/ECA_DOM_Thermodynamics.

Dissolved organic matter↗
NASA NTRS1974

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. A diffraction pattern analysis of MSS images led to the development of spatial signatures for farm land, urban areas and mountains. Four spatial features are employed to describe the spatial characteristics of image cells in the digital data. Three spectral features are combined with the spatial features to form a seven dimensional vector describing each cell. Then, the classification of the feature vectors is accomplished by using the maximum likelihood criterion. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. It was also determined that for a given geographic area the statistics of the classes remain invariable for a period of a month, but vary substantially between seasons. Three ERTS-1 images from the Phoenix, Arizona area were processed, and recognition rates between 85% and 100% were obtained for the terrain classes of desert, farms, mountains, and urban areas. To eliminate the need for training data, a new clustering algorithm has been developed. Seven ERTS-1 images from four test sites have been processed through the clustering algorithm, and high recognition rates have been achieved for all terrain classes.

Gramenopoulos, N.↗
NASA NTRS1976

High-beta turbulence in two-dimensional magnetohydrodynamics

Equations of ideal magnetohydrodynamics are used to study incompressible turbulent flows in a specified geometry where all the field quantities vary with only two spatial dimensions. The procedures adopted are basically those of Kraichnan (1967), in which classical equilibrium ensembles are built around constants of the motion identified from the Fourier-transformed equations of motion. Once the constants of the motion are identified, the statistical formulation of the problem is presented in a phase space whose coordinates are the real and imaginary parts of the Fourier coefficients. Canonical ensembles are constructed in this phase space by classical arguments. The theory developed permits isolation of some qualitatively new gross physical effects which have so far not been calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long-wavelength, self-generated, magnetic fields for a wide range of initial parameters.

Fyfe, D.↗
NASA NTRS1982

The possible modifications of the Hisse model for pure LANDSAT agricultural data

An idea, due to A. Feiveson, is presented for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field). Theoretical arguments are given which show that any significant refinement of the model beyond Feiveson's proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Guassian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗
NASA NTRS1981

The possible modifications of the HISSE model for pure LANDSAT agricultural data

A method for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field) is discussed. Theoretical arguments are given which show that any significant refinement of the model beyond this proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Gaussian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗
NASA NTRS2011

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.↗
NASA NTRS2011

Recent Advances on INSAR Temporal Decorrelation: Theory and Observations Using UAVSAR

We review our recent advances in understanding the role of temporal decorrelation in SAR interferometry and polarimetric SAR interferometry. We developed a physical model of temporal decorrelation based on Gaussian-statistic motion that varies along the vertical direction in forest canopies. Temporal decorrelation depends on structural parameters such as forest height, is sensitive to polarization and affects coherence amplitude and phase. A model of temporal-volume decorrelation valid for arbitrary spatial baseline is discussed. We tested the inversion of this model to estimate forest height from model simulations supported by JPL/UAVSAR data and lidar LVIS data. We found a general good agreement between forest height estimated from radar data and forest height estimated from lidar data.

polarimetric SAR interferometry↗
DOE OSTI2020

Evaluating proxy influence in assimilated paleoclimate reconstructions – Testing the exchangeability of two ensembles of spatial processes

Climate field reconstructions (CFR) attempt to estimate spatiotemporal fields of climate variables in the past using climate proxies such as tree rings, ice cores, and corals. Data Assimilation (DA) methods are a recent and promising new means of deriving CFRs that optimally fuse climate proxies with climate model output. Despite the growing application of DA-based CFRs, little is understood about how much the assimilated proxies change the statistical properties of the climate model data. To address this question, we propose a robust and computationally efficient method, based on functional data depth, to evaluate differences in the distributions of two spatiotemporal processes. In this study, we apply our test to study global and regional proxy influence in DA-based CFRs by comparing the background and analysis states, which are treated as two samples of spatiotemporal fields. We find that the analysis states are significantly altered from the climate-model-based background states due to the assimilation of proxies. Moreover, the difference between the analysis and background states increases with the number of proxies, even in regions far beyond proxy collection sites. Our approach allows us to characterize the added value of proxies, indicating where and when the analysis states are distinct from the background states.

54 ENVIRONMENTAL SCIENCES↗