Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-energy emission from a magnetar giant flare in the Sculptor galaxy

Magnetars are the most highly magnetized neutron stars in the cosmos (with magnetic field 10 13 –10 15 G). Giant flares from magnetars are rare, short-duration (about 0.1 s) bursts of hard X-rays and soft γ rays. Owing to the limited sensitivity and energy coverage of previous telescopes, no magnetar giant flare has been detected at gigaelectronvolt (GeV) energies. Here, we report the discovery of GeV emission from a magnetar giant flare on 15 April 2020 (refs. 3,4 and A. J. Castro-Tirado et al., manuscript in preparation). The Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope detected GeV γ rays from 19 s until 284 s after the initial detection of a signal in the megaelectronvolt (MeV) band. Our analysis shows that these γ rays are spatially associated with the nearby (3.5 megaparsecs) Sculptor galaxy and are unlikely to originate from a cosmological γ-ray burst. Thus, we infer that the γ rays originated with the magnetar giant flare in Sculptor. Furthermore, we suggest that the GeV signal is generated by an ultra-relativistic outflow that first radiates the prompt MeV-band photons, and then deposits its energy far from the stellar magnetosphere. After a propagation delay, the outflow interacts with environmental gas and produces shock waves that accelerate electrons to very high energies; these electrons then emit GeV γ rays as optically thin synchrotron radiation. This observation implies that a relativistic outflow is associated with the magnetar giant flare, and suggests the possibility that magnetars can power some short γ-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Global monthly sectoral water use for 2010–2100 at 0.5° resolution across alternative futures

Water usage is closely linked with societal goals that are both local and global in scale, such as sustainable development and economic growth. It is therefore of value, particularly for long-term planning, to understand how future sectoral water usage could evolve on a global scale at fine resolution. Additionally, future water usage could be strongly shaped by global forces, such as socioeconomic and climate change, and the multi-sector dynamic interactions those forces create. We generate a novel global gridded monthly sectoral water withdrawal and consumption dataset at 0.5° resolution for 2010–2100 for a diverse range of 75 scenarios. The scenarios are harmonized with the five Shared Socioeconomic Pathways (SSPs) and four Representative Concentration Pathways (RCPs) scenarios to support its usage in studies evaluating the implications of uncertain human and earth system change for future global and regional dynamics. To generate the data, we couple the Global Change Analysis Model (GCAM) with a land use spatial downscaling model (Demeter), a global hydrologic framework (Xanthos), and a water withdrawal downscaling model (Tethys).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

MEP: Mobility Energy Productivity Metric

In response to the need to assess fundamental improvements in mobility as technology provides new modes, methods, and opportunities in urban areas to connect people with goods, services, and employment, a mobility-energy productivity score (MEPS) has been developed. The MEPS fundamentally measures the potential of a multi-modal system to connect a person to a variety of services and activities that define a high-quality of life. The MEPS is measured relative to the convenience and energy needed to provide connectivity. Fundamentally derived from accessibility theory, the MEPS advances practice by using readily available travel time data (either from web-based APIs or outputs from an urban transportation model) combined with established parameters that reflect the energy intensity of various travel modes, and relative frequency of activity engagement. The construction of the MEPS allows for aggregation and disaggregation to the appropriate spatial, modal, and trip purpose resolution, as analysis needs dictate. The MEPS could be used to compare alternative futures related to technology, infrastructure investment, or policy, providing a much needed tool for planners, researchers, and analysts. In other words, the metric measures how efficiently connected a place is.

Hou, Yi↗
DOE OSTI2020

A Novel and Practical Method to Quantify the Quality of Mobility: The Mobility Energy Productivity Metric: Preprint

Recent technology innovations are enabling fundamental improvements in mobility systems, including options for new travel modes, methods, and opportunities to connect people with goods, services, and employment. A desire to quantify and compare both existing and emerging transportation options motivated development of the mobility energy productivity (MEP) metric described herein. The MEP metric fundamentally measures the potential of a city's transportation system to connect a person to a variety of services and activities that define a high-quality of life, relative to the convenience, cost and energy needed to provide these connections. Fundamentally derived from accessibility theory, the MEP advances practice by using readily available travel time data (either from web-based application programming interfaces (APIs) or outputs from an urban transportation model) combined with established parameters that reflect the energy intensity and cost of various travel modes, and relative frequency of activity engagement. The construction of the MEP metric allows for aggregation and disaggregation to the appropriate spatial, modal, and trip purpose resolution, as analysis needs dictate. The MEP could be used to compare alternative futures related to technology, infrastructure investment, or policy, providing a much-needed tool for planners, researchers, and analysts.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Characterization of Nitrate, Nitrite, Ammonia, and Tritium in D0220 Cores

This report describes characterization of nitrogen species (nitrate, nitrite, ammonia) in aqueous and solid phases in cores taken in borehole D0220 (well 299-E25-245) under crib 216-A-37-1 at the Hanford Site in 2022 to evaluate (a) the types of nitrogen species currently in the vadose zone and (b) the migration of nitrogen species in the vadose zone. Ammonia and tritium from PUREX decladding condensate were sporadically discharged to the crib from 1977 to 1989. During discharge operations, the estimated travel time through the vadose zone to groundwater was 2.5 to 9 months. By 2003, a characterization borehole (C4106) showed residual pore water with elevated tritium and nitrate in the first 100 ft of the vadose zone, which is likely from the crib. However, there was significantly more nitrate present at the shallowest depth (15 to 22 ft), which may indicate a different source for nitrate, such as adsorbed ammonia slowly being oxidized or nitrogen species precipitates slowly dissolving. In this study, more extensive nitrogen species characterization was done on D0220 cores (drilled in 2022) at 40- and 262-ft depths, which included (a) aqueous nitrate, nitrite, and tritium; (b) adsorbed ammonia; (c) nitrogen in carbonates (or other minerals dissolved in acidic acid); (d) nitrogen in iron oxides (or other minerals dissolved in oxalic acid); and (e) nitrogen in hard-to-extract minerals (minerals dissolved in nitric acid). High pore water nitrate (226 to 331 mg/L) at 40-ft depth measured in D0220 (2022) compared to 60 mg/L at 40-ft depth in C4106 (2003) may indicate nitrate is migrating deeper. Tritium concentrations (pore water 132,000 to 148,000 pCi/L) measured in D0220 at 40-ft depth in 2022 were considerably higher than in C4106 at 40-ft depth (160 pCi/L). Additional nitrogen species mass was present in adsorbed and precipitated phases in D0220 cores. Low adsorbed ammonia was measured at 40.1- and 261.7-ft depths. Low concentrations of carbonate associated nitrogen and iron oxide-associated nitrogen were present at 40.1 and 40.6 ft depths. Nitrogen species in solid phase extractions indicate nitrogen precipitates or aqueous nitrate or ammonia trapped in sediment microfractures that are coated by precipitates. Overall, the nitrogen species and tritium characterization at two depths in D0220 showed that additional nitrogen species were present in the vadose zone. In addition, elevated pore water nitrate in D0220 at 40 ft depth from 2022 may indicate vertical migration compared to C4106 nitrate profile from 2003. Analysis of additional depths in D0220 and spatial variability of the nitrate plume along the length of the crib (from surface electrical resistivity) would be extremely useful for this comparison.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High Spatial Resolution Mapping of Retained Fission Gas

Fission gas isotopic analysis provides quantitative high precision determination of irradiated nuclear fuel burnup, offers diagnostic value, and informs fuel performance models. A measurement capability has been developed at Idaho National Laboratory (INL) for the release of retained fission gas using a focused laser and static noble gas mass spectrometry (MS) analysis. This high resolution (10s microns spot size) capability was demonstrated using Xe implanted metal foils.

07 - ISOTOPES AND RADIATION SOURCES↗
DOE OSTI2025

Risk Assessment in a Chemical Laboratory Following an Explosive Incident Involving a Novel Diazonium Compound: Retrospective Analysis and Lessons Learned

Diazonium compounds are synthetically useful in the production of dyes and textiles, however they are highly explosive under dry conditions. Explosion prevention becomes more difficult when new diazonium compounds are synthesized, because while some syntheses include a counterion to increase their stability, this is not always a reliable method to prevent an explosive incident. Due to the uncertainty surrounding the explosiveness of different diazonium compounds, it is important to understand how to safely clean up after an incident and how to determine when it is safe to return a laboratory to typical operational use, particularly when the incident involves a novel compound where a standard does not exist for instrument calibration. Here, an explosive event is discussed involving the synthesis of 4-bromo-benzenediazonium-2-carboxylate. Following the explosive incident and 3-step cleanup, which involved a precautionary neutralization step, samples were collected from the fume hood where the incident occurred. Because the incident involved an unstable, novel compound that is not commercially available and was deemed unsafe to resynthesize for instrument calibration, we assessed the risk of further explosion by analyzing for the stable decomposition products. Mass spectrometry analysis confirmed that the residue in the fume hood contained 5-bromosalicylic acid, a decomposition product of 4-bromo-benzenediazonium-2-carboxylate. Samples were taken from multiple points in the fume hood and analyzed to estimate the spatial distribution of the decomposition product. Based on this analysis, we inferred that the primary decomposition product was far more abundant than residual energetic, indicating the energetic had been consumed or neutralized to a trace quantity where the risk of further explosion was low. Furthermore, the steps presented here─specifically, initial neutralization and then analyzing the spatial distribution of expected decomposition products to assess risk when a novel explosive material is detonated in a confined space─were our approach to assess further risk following an explosion due to a novel diazonium compound without the need for any further handling or resynthesis of the energetic. Here, we present our approach and critically analyze these steps by discussing retrospective lessons learned and alternative analytical approaches.

Computer simulations↗
DOE OSTI2023

ENSO-Induced Teleconnection: Process-Oriented Diagnostics to Assess Rossby Wave Sources and Ambient Flow Properties in Climate Models

Climate model fidelity in representing ENSO-induced teleconnection is assessed with process-oriented diagnostics that examine a chain of processes, from equatorial Pacific precipitation to the midlatitude circulation pattern over the Pacific–North American regions. Such processes are rarely addressed during model development. Using an upper-tropospheric divergent level, local vorticity gradient of the ambient zonal flow (∂ 2 U¯/∂y 2 ) and a restoring force for Rossby waves (β * ) are estimated, the equivalent barotropic vorticity equation is solved, and an anomalous Rossby wave source (RWS') quantified. The analysis is applied to AMIP5 and AMIP6 simulations. For a realistic circulation response representation, the hypothesis that models accurately represent the strength and location of RWS', and spatial variations in β * is tested. Compared to AMIP5, in AMIP6 there are clear improvements in representing RWS' and β * . To validate the hypothesis, the analysis identifies two metrics: spatially coherent RWS' in the subtropical North Pacific, and longitudes of negative β * over the western-central North Pacific. By projecting these metrics in two and three-dimensional views, improvements or degradations in model versions are apparent. If a model’s fidelity in representing ∂ 2 U¯/∂y 2 and RWS' are compromised, then radiated Rossby waves are reflected more equatorward, resulting in zonally elongated circulation anomalies over the central North Pacific. Thus, during climate model development, applying this analysis frequently will keep a regular check on the fidelity of the modeled response to anomalous El Niño convection in conjunction with changing model ambient flow dependencies. As a result, this analysis is intended to form a process-oriented diagnostics package, a community contribution to the NOAA Model Diagnostics Task Force.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Willow Biomass Crops Are a Carbon Negative or Low-Carbon Feedstock Depending on Prior Land Use and Transportation Distances to End Users

Few life cycle assessments (LCAs) on willow biomass production have investigated the effects of key geographically specific parameters. This study uses a spatial LCA model for willow biomass production to determine spatially explicit greenhouse gas (GHG) emissions and energy return on investment (EROI), including land use conversion from pasture and cropland or grassland. There were negative GHG emissions on 92% of the land identified as suitable for willow biomass production, indicating this system’s potential for climate change mitigation. For willow planted on cropland or pasture, life cycle GHG emissions ranged from −53.2 to −176.9 kg CO2eq Mg-1. When willow was grown on grassland the projected decrease in soil organic carbon resulted in a slightly positive GHG balance. Changes in soil organic carbon (SOC) associated with land use change, transportation distance, and willow yield had the greatest impacts on GHG emissions. Results from the uncertainty analysis exhibited large variations in GHG emissions between counties arising from differences in these parameters. The average EROI across the entire region was 19.2. Willow biomass can be a carbon negative or low-carbon energy source with a high EROI in regions with similar infrastructure, transportation distances, and growing conditions such as soil characteristics, land cover types, and climate.

09 BIOMASS FUELS↗
DOE OSTI2026

Spatially resolved phase transformation mapping in 410 stainless steel during additive manufacturing visualized via real-time infrared data

Here, in this study, we demonstrate that spatially resolved cooling curves derived from real-time infrared (IR) thermography during additive manufacturing (AM) can capture spatial variations in phase transformation temperatures through cooling curve analysis (CCA). Using this approach, we show that during laser hot-wire deposition of 410 stainless steel (410SS), the martensite start temperature (M s ) evolves dynamically throughout the build. The M s temperature is spatially nonuniform, ranging from 185 °C to 348 °C, with the lowest values toward the build center and higher values toward the upper region of the deposit. In the lower portion of the build, no M s inflection is detected via CCA, consistent with transformation occurring earlier during thermal cycling followed by tempering during subsequent thermal cycles. These trends in M s are corroborated by characterizing the microstructure by electron backscatter diffraction (EBSD). Traditionally, M s is assumed to be constant, and a single interpass temperature is applied during both deposition and residual stress modeling. Our results demonstrate that IR-derived cooling curves provide a route to spatially and temporally resolved transformation temperature tracking for dynamic interpass control and improved residual-stress modeling.

Additive manufacturing↗
DOE OSTI2025

Non-stationary precipitation design standards for stormwater infrastructure modernization at USAF installations

The resilience of defense infrastructure systems to a changing climate is critical for national security. Climate induced recurrent flooding is already impacting over 20 U.S. Air Force installations, underscoring the urgency of revisiting precipitation standards and stormwater infrastructure design. Despite growing scientific knowledge and an expanding set of tools for updating outdated precipitation standards based on the assumption of climate stationarity, the adoption of climate informed analyses remain limited in practice. This study utilizes an existing framework to update Intensity (or Depth)-Duration-Frequency (DDF) curves using an ensemble of future climate projections. Change factors in precipitation estimates are derived and applied to six USAF installations across the U.S. The analysis is further extended to evaluate the implications of climate-informed DDFs on stormwater infrastructure performance and flood analysis at Tyndall AFB. Results indicate that the current design precipitation estimates are likely to become obsolete in all six USAF bases by the end of the century. The wide range of change factors across 32 GCM ensembles highlights the need to integrate uncertainty and evolving scientific data into infrastructure planning. The study also finds that the impacts of a changing climate vary spatially and temporally, emphasizing the value of localized analysis for infrastructure decision-making. The work advances ongoing DoD and societal efforts to implement adaptation strategies aimed at enhancing infrastructure resilience.

Intensity-duration-frequency curves↗
DOE OSTI2026

Scaled Implementation of Smart Charge Management for Electric Vehicles

Smart charge management (SCM) has become a critical strategy for mitigating potential grid impacts and reducing electricity costs for all customers. This study will evaluate the economics of implementing light-duty EV SCM at scale across the United States. This work will enhance distribution system analysis by estimating SCM implementation costs, exploring viable business cases, and developing a framework for national-level applications. The primary methodology involves leveraging detailed grid modeling from a specific service territory and using spatial extrapolation techniques to generalize findings to other regions. The analysis will develop key metrics to quantify the costs and benefits of SCM, including implementation costs relative to strategy and scale, the cost of distribution system upgrades with and without SCM, and the percentage of peak-load reduction. The objective is to produce a comprehensive report and a parameterized framework that enables utilities to self-assess the value of SCM in their own service territories. This will support the development of cost-effective charging strategies, accelerate the energization of new EV chargers, and facilitate the seamless integration of EVs into the nation's power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Characterization of Surface Heterogeneity-Induced Convection Using Cluster Analysis

Land-atmosphere interactions and boundary layer processes often control the formation of shallow clouds and subsequently deep convective precipitation over the Southern Great Plains. In this study, we examine the impacts of large-scale advection on the cloud populations and land-atmospheric coupling observed during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign in 2016. In this work, we performed two Large Eddy Simulations (LES) using the Weather Research and Forecasting Model for a day in which the transition from clear-sky to shallow precipitating clouds and isolated deep convection was observed. The control simulation reproduced the overall distribution of cloud populations by using realistic soil conditions with an interactive land model. In the sensitivity simulation where large-scale advection is removed, a strong relationship between the land and boundary layer is found. To study the timing, location, and intensity of convective initiation and the relationship of clouds with land surface properties, a cluster analysis of equivalent potential temperature is performed for the simulation without large-scale advection. That analysis shows that convective clouds first form over regions with higher surface sensible heat flux. Precipitation from those convective clouds likely triggers new updrafts nearby about two hours later through the lifting associated with cold pools. Furthermore, the cluster analysis also shows that in addition to the spatial pattern of soil moisture, land use and soil texture in western Oklahoma also influence the location of convective initiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Solution Irregularity Remediation for Spatial Discretization Error Estimation for S N Transport Solutions

The discrete ordinates linear Boltzmann transport equation is typically solved in its spatially discretized form, incurring spatial discretization error. Quantification of this error for purposes such as adaptive mesh refinement or error analysis requires an a posteriori estimator, which utilizes the numerical solution to the spatially discretized equation to compute an estimate. Because the quality of the numerical solution informs the error estimate, irregularities, present in the true solution for any realistic problem configuration, tend to cause the largest deviation in the error estimate vis-a-vis the true error. In this paper, an analytical partial singular characteristic tracking (pSCT) procedure for reducing the estimator’s error is implemented within our novel residual source estimator for a zeroth-order discontinuous Galerkin scheme, at the additional cost of a single inner iteration. Here, a metric-based evaluation of the pSCT scheme versus the standard residual source estimator is performed over the parameter range of a Method of Manufactured Solutions test suite. The pSCT scheme generates near-ideal accuracy in the estimate in problems where the dominant source of the estimator’s error is the solution irregularity, namely, problems where the true solution is discontinuous and problems where the true solution’s first derivative is discontinuous and the scattering ratio is low. In problems where the scattering ratio is high and the true solution is discontinuous in the first derivative, the error in the scattering source, which is not converged by the pSCT scheme, is greater than the error incurred due to the irregularity. Ultimately, a pSCT scheme is judged to be useful for error estimation in problems where the computational cost of the scheme is justified. In the presence of many irregularities, such a scheme may be intractable for general use, but in benchmarks, as an analytical tool, or in problems that have nondissipative discontinuities, the scheme may prove invaluable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Unveiling Spatial and Temporal Dynamics of Plasmon-Enhanced Localized Fields in Metallic Nanoframes through Ultrafast Electron Microscopy

Plasmonic nanomaterials, particularly noble metal nanoframes (NFs), are important for applications such as catalysis, biosensing, and energy harvesting due to their ability to enhance localized electric fields and atomic efficiency via localized surface plasmon resonance (LSPR). Yet the fundamental structure-function relationships and plasmonic dynamics of the NFS are difficult to study experimentally and thus far rely predominately on computational methodologies, limiting their utilization. This study leverages the capabilities of ultrafast electron microscopy (UEM), specifically photon-induced near-field electron microscopy (PINEM), to probe the light-matter interactions within plasmonic NF structures. Here, the effects of shape, size, and plasmonic coupling of Pt@Au core-shell NFs on spatial and temporal characteristics of plasmon-enhanced localized electric fields are explored. Importantly, time-resolved PINEM analysis reveals that the plasmonic fields around hexagonal NF prisms exhibit a spatially dependent excitation and decay rate, indicating a nuanced interplay between the spatial geometry of the NF and the temporal evolution of the localized electric field. These results and observations uncover nanophotonic energy transfer dynamics in NFs and highlight their potential for applications in biosensing and photocatalysis.

Plasmonics↗
DOE OSTI2024

xCDAT: A Python Package for Simple and Robust Analysis of Climate Data

xCDAT (Xarray Climate Data Analysis Tools) is an open-source Python package that extends Xarray (Hoyer & Hamman, 2017) for climate data analysis on structured grids. xCDAT streamlines analysis of climate data by exposing common climate analysis operations through a set of straightforward APIs. Some of xCDAT’s key features include spatial averaging, temporal averaging, and regridding. These features are inspired by the Community Data Analysis Tools (CDAT) library (Dean N. Williams et al., 2009) (D. N. Williams, 2014) (Doutriaux et al., 2019) and leverage powerful packages in the Xarray ecosystem including xESMF (Zhuang et al., 2023), xgcm (Abernathey et al., 2022), and CF xarray (Cherian et al., 2023). To ensure general compatibility across various climate models, xCDAT operates on datasets that are compliant with the Climate and Forecast (CF) metadata conventions (Hassell et al., 2017).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Benchmarking the PCMCI Causal Discovery Algorithm for Spatiotemporal Systems

Causal discovery algorithms construct hypothesized causal graphs that depict causal dependencies among variables in observational data. While powerful, the accuracy of these algorithms is highly sensitive to the underlying dynamics of the system in ways that have not been fully characterized in the literature. In this report, we benchmark the PCMCI causal discovery algorithm in its application to gridded spatiotemporal systems. Effectively computing grid-level causal graphs on large grids will enable analysis of the causal impacts of transient and mobile spatial phenomena in large systems, such as the Earth’s climate. We evaluate the performance of PCMCI with a set of structural causal models, using simulated spatial vector autoregressive processes in one- and two-dimensions. We develop computational and analytical tools for characterizing these processes and their associated causal graphs. Our findings suggest that direct application of PCMCI is not suitable for the analysis of dynamical spatiotemporal gridded systems, such as climatological data, without significant preprocessing and downscaling of the data. PCMCI requires unrealistic sample sizes to achieve acceptable performance on even modestly sized problems and suffers from a notable curse of dimensionality. This work suggests that, even under generous structural assumptions, significant additional algorithmic improvements are needed before causal discovery algorithms can be reliably applied to grid-level outputs of earth system models.

54 ENVIRONMENTAL SCIENCES↗