Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial variability”

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

COMPASS-FME Synoptic Site Characterization

This dataset contains soil biogeochemical and physicochemical characterization data for the COMPASS-FME synoptic sites.This dataset also contains data for the paper Patel et al. 2025 "Transition zones at the changing coastal terrestrial-aquatic interface", https://doi.org/10.1029/2025JG008978.Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience unique spatial and temporal variability in biogeochemical processes, driven by wetland-to-upland gradients and hydrological fluctuations. We studied drivers of coastal soil variability (a) at regional scales and (b) across transects from upland forest to wetland, in two contrasting regions — Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity-related analytes were a key driver of soil variability, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI – however, contrary to expectations, transition soils were not consistently intermediate between upland and wetland endmembers; the non-monotonic trends of carbon, phosphorus, iron along our transects suggest that these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients provide insights into which soil processes may act as precursors to ecosystem shifts. Our comprehensive soil characterization across the coastal transects provides essential data for mechanistic modeling of ecosystem dynamics.The data are provided as processed, csv files. Raw data and processing scripts can be accessed on GitHub (https://github.com/COMPASS-DOE/cmps-soil_characterization).A note on the nomenclature: the experimental design represents three points along the coastal gradient -- upland, transition, and wetland. "wetland" is referred to as "marsh" in the corresponding paper. The two terms can be used interchangeably for the sites in this study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Prestack Waveform Inversion at GC 955: Trials and sensitivity of PWI to high-resolution seismic data (Data Report)

Understanding the spatial extent and variability of hydrate deposits has important implications for potential economic production, climate change, and assessing natural hazards risks. Seismic reflection techniques are capable of determining the extent of gas hydrate deposits and can shed light on hydrate concentrations and properties. Using high resolution prestack time migrated seismic data collected during a 2013 US Geological Survey expedition and prestack waveform inversion (PWI) techniques, we produce highly resolved velocity models and compare them to co-located well logs. Coupling our PWI results with velocity-porosity relationships and nearby well control, we compare hydrate properties at site 955 in the Green Canyon and site 313 at Walker Ridge.

03 NATURAL GAS↗
DOE OSTI2023

Physics-based Induced Earthquake Forecasting: Process Understanding, and Hazards Mitigation

Disposal of saltwater co-produced with oil and gas is linked to elevated seismicity in the Central and Midwest US. There is a concern that these events may lead to widespread damage and an overall increase in seismicity. Thus an improved understanding of the spatially and temporally variable deformation and stress field associated with fluid injection operation is critically necessary for evaluating time-varying seismic hazards. Despite the improvements in seismic monitoring capacity and the resulting decrease in the magnitude detection threshold, estimates of induced earthquake probability remain elusive due to insufficient models incapable of accounting for the complex physics governing the process of induced seismicity. The proposed research effort will comprehensively analyze, integrate, and interpret geodetic, injection, and seismic data in the vicinity of the injection sites in Oklahoma to resolve the 4-dimensional distribution of pore pressure and stress in the shallow crust. This project, in particular, is focused on exploring the statistical relation between injection operation and increased earthquake hazard. The amplitude of and the extent to which pore pressure changes are determined by some factors, in particular, the hydrogeological properties of the rocks, such as diffusivity. Thus the available deformation data will be used to constrain hydrogeological properties of the medium, to accurately resolve the evolution of crustal stresses due to fluid injection. Having the time-varying models of stress changes, a statistical framework will be implemented to estimate the time-dependent probability of large earthquakes on the nearby fault systems. These data and models help to improve seismic hazard estimates and aid in constructing operational-induced earthquake forecast models. This information can also be integrated into the updated U.S. National Seismic Hazard Map, which local communities and authorities use in their earthquake risk estimates and mitigation efforts.

58 GEOSCIENCES↗
DOE OSTI2026

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Degradability of Biodegradable Soil Moisture Sensor Components and Their Effect on Maize (Zea mays L.) Growth

Inexpensive and no-maintenance biodegradable soil moisture sensors could improve existing knowledge on spatial and temporal variability of available soil water at field-scale. Such sensors can unlock the full potential of variable-rate irrigation (VRI) systems to optimize water applications in irrigated cropping systems. The objectives of this study were to assess (i) the degradation of soil moisture sensor component materials and (ii) the effects of material degradation on maize (Zea Mays L.) growth and development. This study was conducted in a greenhouse at Colorado State University, Colorado, USA, by planting maize seeds in pots filled with three growing media (field soil, silica sand, and Promix commercial potting media). The degradation rate of five candidate sensor materials (three blends of beeswax and soy wax, balsa wood, and PHBV (poly(3-hydroxybutyrate-co-3-hydroxyvalerate))) was assessed by harvesting sensor materials at four maize growth stages (30, 60, 90, and 120 days after transplanting). All materials under consideration showed stability in terms of mass and dimension except PHBV. PHBV was degraded entirely within 30 days in soil and Promix, and within 60 days in sand. Balsa wood did now show any significant reduction in mass and dimensions in all growth media. Similarly, there was no significant mass loss across wax blends (p = 0.05) at any growth stage, with a few exceptions. Among the wax blends, 3:1 (beeswax:soy wax) was the most stable blend in terms of mass and dimension with no surface cracks, making it a suitable encapsulant for soil sensor. All materials under consideration did not have any significant effect on maize growth (dry biomass, green biomass, and height) as compared to control plants. These results indicated that 3:1 beeswax:soy wax blend, PHBV, and balsa wood could be suitable candidates for various components of biodegradable soil moisture sensors.

Dahal, Subash (ORCID:0000000205489103)↗
DOE OSTI2026

Data from: Coupled machine learning-ecosystem ensemble models substantially improve predictions of nitrous oxide (N 2 O) fluxes from US croplands

Nitrous oxide (N₂O) is a potent and persistent greenhouse gas, with rising atmospheric concentrations driven in part by inefficient use of synthetic nitrogen (N) fertilizers in agriculture. Predicting soil N₂O emissions is challenging due to high spatial and temporal variability arising from complex soil biogeochemical processes. Process-based ecosystem models and standalone machine learning (ML) approaches without extensive site-specific calibration often miss high emission episodes. Here, we show how an Ensemble Modeling System (EMS) based on outputs from an ensemble of ecosystem models coupled to an ensemble of ML models can improve predictions and understanding of N2O fluxes from US cropland. Trained and validated on approximately 12,000 N2O chamber measurements at 17 U.S. Midwest sites (six crops, 35 management practices), the EMS accurately predicted daily fluxes of N2O at both training (R² = 0.84, RMSE = 16.4 g N ha⁻¹ d⁻¹) and held-out testing sites (R² = 0.84, RMSE = 6.2 g N ha⁻¹ d⁻¹). Analyses identified six dominant N₂O drivers: soil organic carbon (SOC), NH₄⁺, NO₃⁻, water-filled pore space (WFPS), soil temperature, and biomass production. Wet, warm soils produced large N₂O peaks only with sufficient SOC and mineral N; in low-SOC soils, fluxes remained low. Incorporating these drivers into process-based models might significantly improve their predictive capacity. The EMS demonstrates a strong potential to predict N₂O fluxes at unseen sites, enabling more reliable regional inventories, improved gap-filling where measurements are sparse, and enhanced understanding of mechanisms to advance targeted mitigation strategies in food, feed, and bioenergy crops.

agricultural sciences↗
DOE OSTI2020

Ship-based measurements of ice nuclei concentrations over the Arctic, Atlantic, Pacific and Southern oceans

Ambient concentrations of ice-forming particles measured during ship expeditions are collected and summarised with the aim of determining the spatial distribution and variability in ice nuclei in oceanic regions. The presented data from literature and previously unpublished data from over 23 months of ship-based measurements stretch from the Arctic to the Southern Ocean and include a circumnavigation of Antarctica. In comparison to continental observations, ship-based measurements of ambient ice nuclei show 1 to 2 orders of magnitude lower mean concentrations. To quantify the geographical variability in oceanic areas, the concentration range of potential ice nuclei in different climate zones is analysed by meridionally dividing the expedition tracks into tropical, temperate and polar climate zones. We find that concentrations of ice nuclei in these meridional zones follow temperature spectra with similar slopes but vary in absolute concentration. Typically, the frequency with which specific concentrations of ice nuclei are observed at a certain temperature follows a log-normal distribution. A consequence of the log-normal distribution is that the mean concentration is higher than the most frequently measured concentration. Finally, the potential contribution of ship exhaust to the measured ice nuclei concentration on board research vessels is analysed as function of temperature. We find a sharp onset of the influence at approximately -36 °C but none at warmer temperatures that could bias ship-based measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using automated machine learning for the upscaling of gross primary productivity

Estimating gross primary productivity (GPP) over space and time is fundamental for understanding the response of the terrestrial biosphere to climate change. Eddy covariance flux towers provide in situ estimates of GPP at the ecosystem scale, but their sparse geographical distribution limits larger-scale inference. Machine learning (ML) techniques have been used to address this problem by extrapolating local GPP measurements over space using satellite remote sensing data. However, the accuracy of the regression model can be affected by uncertainties introduced by model selection, parameterization, and choice of explanatory features, among others. Recent advances in automated ML (AutoML) provide a novel automated way to select and synthesize different ML models. In this work, we explore the potential of AutoML by training three major AutoML frameworks on eddy covariance measurements of GPP at 243 globally distributed sites. We compared their ability to predict GPP and its spatial and temporal variability based on different sets of remote sensing explanatory variables. Explanatory variables from only Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and photosynthetically active radiation explained over 70 % of the monthly variability in GPP, while satellite-derived proxies for canopy structure, photosynthetic activity, environmental stressors, and meteorological variables from reanalysis (ERA5-Land) further improved the frameworks' predictive ability. We found that the AutoML framework Auto-sklearn consistently outperformed other AutoML frameworks as well as a classical random forest regressor in predicting GPP but with small performance differences, reaching an r 2 of up to 0.75. We deployed the best-performing framework to generate global wall-to-wall maps highlighting GPP patterns in good agreement with satellite-derived reference data. This research benchmarks the application of AutoML in GPP estimation and assesses its potential and limitations in quantifying global photosynthetic activity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A restructured and updated global soil respiration database (SRDB-V5)

Field-measured soil respiration (R S , the soil-to-atmosphere CO 2 flux) observations were compiled into a global soil respiration database (SRDB) a decade ago, a resource that has been widely used by the biogeochemistry community to advance our understanding of R S dynamics. Novel carbon cycle science questions require updated and augmented global information with better interoperability among datasets. Here, we restructured and updated the global R S database to version SRDB-V5. The updated version has all previous fields revised for consistency and simplicity, and it has several new fields to include ancillary information (e.g., R S measurement time, collar insertion depth, collar area). The new SRDB-V5 includes published papers through 2017 (800 independent studies), where total observations increased from 6633 in SRDB-V4 to 10 366 in SRDB-V5. The SRDB-V5 features more R S data published in the Russian and Chinese scientific literature and has an improved global spatio-temporal coverage and improved global climate space representation. We also restructured the database so that it has stronger interoperability with other datasets related to carbon cycle science. For instance, linking SRDB-V5 with an hourly timescale global soil respiration database (HGRsD) and a community database for continuous soil respiration (COSORE) enables researchers to explore new questions. The updated SRDB-V5 aims to be a data framework for the scientific community to share seasonal to annual field R S measurements, and it provides opportunities for the biogeochemistry community to better understand the spatial and temporal variability in R S , its components, and the overall carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

CMIP5 model selection for ISMIP6 ice sheet model forcing: Greenland and Antarctica

Abstract. The ice sheet model intercomparison project for CMIP6 (ISMIP6) effort brings together the ice sheet and climate modeling communities to gain understanding of the ice sheet contribution to sea level rise. ISMIP6 conducts stand-alone ice sheet experiments that use space- and time-varying forcing derived from atmosphere–ocean coupled global climate models (AOGCMs) to reflect plausible trajectories for climate projections. The goal of this study is to recommend a subset of CMIP5 AOGCMs (three core and three targeted) to produce forcing for ISMIP6 stand-alone ice sheet simulations, based on (i) their representation of current climate near Antarctica and Greenland relative to observations and (ii) their ability to sample a diversity of projected atmosphere and ocean changes over the 21st century. The selection is performed separately for Greenland and Antarctica. Model evaluation over the historical period focuses on variables used to generate ice sheet forcing. For stage (i), we combine metrics of atmosphere and surface ocean state (annual- and seasonal-mean variables over large spatial domains) with metrics of time-mean subsurface ocean temperature biases averaged over sectors of the continental shelf. For stage (ii), we maximize the diversity of climate projections among the best-performing models. Model selection is also constrained by technical limitations, such as availability of required data from RCP2.6 and RCP8.5 projections. The selected top three CMIP5 climate models are CCSM4, MIROC-ESM-CHEM, and NorESM1-M for Antarctica and HadGEM2-ES, MIROC5, and NorESM1-M for Greenland. This model selection was designed specifically for ISMIP6 but can be adapted for other applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

BIRDIE beryllium-10 and -7 measurements from 2022 summer deep convection rain at ARM SGP

Water samples were manually collected from 11 SGP ARM sites in the early summer of 2022 to measure cosmogenic isotopes of beryllium. The water samples were brought back to the Center for Accelerator Mass Spectrometry (CAMS) at LLNL and analysed for beryllium-10 and beryllium-7 as well as for chlorine-37, chlorine-36, and tritium. A total of six storm events were measured over three weeks. At least five of those storms included overshooting cloud events. Available data include precipitation amount, beryllium-10 concentration, beryllium-7 concentration, and 10/7Be ratio. Spatial and temporal variability in the cosmogenic isotopes was relatively high but could not directly be linked to precipitation amount. Analysis is on-going to explain the measured variability and to see if it can be linked to cloud top height and overshooting depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elastic Changepoint Detection for Globally-indexed Functional Time Series Data with Climate Applications

Changepoint detection is a vital tool in the application of climate data analysis. Numerous types of climate observation data are most properly represented by functional time series, implying a need for accurate changepoint detection methods applicable to functional time series data. Such data taken at a global scale often contain both spatial heterogeneity and dependence as well as phase (time) misalignment. In this report, we present methods which can detect spatially-dependent changepoints while allowing different estimates of change time and change strength depending on location. Additionally, we provide extensions to this spatially-predicted model which controls for phase variability among observations. Our methods provide the ability to detect a single change, or control for epidemic changes (where a “return-to-normal” change is more likely to be detected than the initial change). We showcase results analyzing the June 1991 eruption of Mt. Pinatubo, where our methods demonstrate the ability to accurately detect both single and epidemic changepoints even in the presence of strong seasonal variability. We find that our spatially-predicted model improves the detection of relevant changepoints versus methods which do not take spatial information into account, and we find that controlling for phase variability helps to control the false discovery rate during the detection process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Characterizing and Mitigating Intraday Variability: Reconstructing Source Structure in Accreting Black Holes with mm-VLBI

The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87* and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse (u, v) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87*. We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87*, finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Meteoric 10Be Flux Calibration Data for the East River Watershed, Colorado, USA

This data package contains tabular and geospatial data used to quantify and model meteoric beryllium-10 fluxes in the East River watershed, Colorado, USA. The tabular component includes calibration-site data from five glacial moraine sites and includes environmental variables used to evaluate spatial controls on meteoric 10Be delivery, including elevation, mean annual precipitation (MAP), mean snow depth, and mean snow water equivalent (SWE). These site-level data were used to compare observed fluxes with environmental gradients across the watershed and to evaluate the effects of erosion correction on flux estimates. The package also includes supporting slope and curvature values used to assess topographic inputs to the erosion analysis. A second component of the data package contains updated manuscript tables and regression outputs used to summarize the relationships between meteoric 10Be flux and environmental predictors. These tables include meteoric 10Be sample information and AMS results, site-level environmental values, site-level meteoric 10Be inventory and flux values, watershed-averaged predicted fluxes, soil bulk density measurements, fine-fraction values, soil pH measurements, and regression statistics including slope, intercept, coefficient of determination, and p-value. The regression products include both standard linear regressions and regressions in which the intercept is constrained to pass through zero, and they support the analyses presented in the companion manuscript. Together, these tabular files provide the numerical basis for the manuscript tables and the regression-based interpretation of meteoric 10Be flux variability in a snow-dominated mountain watershed. The geospatial component of the package consists of GeoTIFF raster files used to generate the map products presented in Figures 2 and 6 of the companion manuscript. These rasters represent watershed-scale spatial layers for environmental variables and regression-based predictions of meteoric 10Be flux. This dataset contains comma-separated values files (.csv), Microsoft Excel files (.xlsx), GeoTIFF raster files (.tif), and upporting metadata files, including CSV data dictionaries and readme text files (.csv, .txt). The tabular files can be opened with standard spreadsheet software, and the raster files can be viewed and analyzed in GIS software such as ArcGIS Pro or QGIS. Together, these files document the numerical and spatial datasets used to calibrate and predict meteoric 10Be delivery in the East River watershed.

East River↗
DOE OSTI2021

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗
DOE OSTI2019

Arctic River Delta Morphologic Variability and Implications for Riverine Fluxes to the Coast

Arctic riverine fluxes are anticipated to increase as the Arctic warms and have a large impact on the Arctic ocean. Deltas modify the spatial and temporal distributions of riverine fluxes, but no thorough studies have been conducted to analyze Arctic delta morphologies to determine their influence on land–ocean fluxes. We performed an analysis of six high–latitude deltas (Colville, Kolyma, Lena, Mackenzie, Yenisei, and Yukon) to characterize delta morphologies and determine the influence of morphology on the distribution of fluxes to the coast. All six deltas deliver material to the coast at discrete locations across small areas despite differences in delta shoreline length. Large Arctic deltas exhibit large variability in channel width, which we hypothesize is due to a feedback with ice cover and retreat that favors the growth of large channels over geologic timescales. Spatial variability in island sizes suggests variability in channel activity, island nourishment, and susceptibility to drowning by sea level rise. Potential lake storage is highest on the Mackenzie delta, thus providing a means for reducing nutrient and sediment loading of the coastal ocean. Here, the connected lakes are also prevalent on the Colville and Yukon deltas, suggesting that these deltas can filter riverine fluxes even when the deltas are not flooded. Differences in Arctic delta morphologies can be explained by varying levels of riverine and marine influence, antecedent topography, and local channel dynamics. Ice cover also plays a large role in controlling Arctic delta morphologies and dynamics that has not been previously represented in interpretations of existing delta metrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Comparison of four learning-based methods for predicting groundwater redox status

Knowing the location where groundwater denitrification occurs, or by proxy the groundwater redox status (oxic, mixed, and anoxic), is valuable information for assessing and managing potential agricultural land-use impacts on freshwater quality. We compare the efficacy of supervised (Linear Discriminant Analysis LDA; Boosted Regression Trees, BRT; and Random Forest, RF) and unsupervised (Modified Self-Organizing Map, MSOM) learning-based methods to predict groundwater redox status in the agriculturally dominated Tasman, Waikato, and Wellington regions of New Zealand. Thresholds applied to regional groundwater-quality samples provide redox status variables and learn heuristics constrained by these variables and applied to spatial factors (climate, elevation, geologic, hydrology soils, and well depth) identify optimal sets of regional predictor variables. A split-sample approach is used to train and test the learning methods ability to predict redox status using the optimal predictor variables. Overall, the supervised methods demonstrate a prediction bias toward oxic conditions and inability to perform statistically well when using independent regional data; for example, consider kappa statistics for BRT (Tasman: 0.42, Waikato: 0.38, Wellington: 0.17), RF (Tasman: 0.42, Waikato: 0.47, Wellington: 0.17 and LDA (Tasman: 0.46, Waikato: 0.32, Wellington: 0.17). By contrast, the unsupervised method performs statistically well when predicting oxic, mixed, and anoxic conditions and corresponding depths when using independent regional data; for example, consider MSOM kappa statistics for Tasman: 0.78, Waikato: 0.80, Wellington: 0.76. The unsupervised learning method provides the added benefits of being (1) able to combine predictions into 3D regional anoxic probability plots for interpreting the spatial influence of paleosols and groundwater flowpaths on redox status, and (2) readily extended to map 3D redox status across New Zealand and other countries despite data bias and sparsity.

machine learning algorithms, Random Forest, booste↗
DOE OSTI2022

A Universal Power-law Prescription for Variability from Synthetic Images of Black Hole Accretion Flows

We present a framework for characterizing the spatiotemporal power spectrum of the variability expected from the horizon-scale emission structure around supermassive black holes, and we apply this framework to a library of general relativistic magnetohydrodynamic (GRMHD) simulations and associated general relativistic ray-traced images relevant for Event Horizon Telescope (EHT) observations of Sgr A*. We find that the variability power spectrum is generically a red-noise process in both the temporal and spatial dimensions, with the peak in power occurring on the longest timescales and largest spatial scales. When both the time-averaged source structure and the spatially integrated light-curve variability are removed, the residual power spectrum exhibits a universal broken power-law behavior. On small spatial frequencies, the residual power spectrum rises as the square of the spatial frequency and is proportional to the variance in the centroid of emission. Beyond some peak in variability power, the residual power spectrum falls as that of the time-averaged source structure, which is similar across simulations; this behavior can be naturally explained if the variability arises from a multiplicative random field that has a steeper high-frequency power-law index than that of the time-averaged source structure. We briefly explore the ability of power spectral variability studies to constrain physical parameters relevant for the GRMHD simulations, which can be scaled to provide predictions for black holes in a range of systems in the optically thin regime. We present specific expectations for the behavior of the M87* and Sgr A* accretion flows as observed by the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗