Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Variables”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

Regression Analysis with the Directed Infusion of Data

Integrating artificial intelligence and machine learning tools into industry necessitates large-scale collaborative efforts that ensure the robust and accurate execution of downstream analytics such as time series prediction, uncertainty quantification, grid optimization, and condition monitoring. However, concerns related to data privacy pervade the nuclear industry due to the proprietary nature of its data and the possibility of data leakage. Legacy techniques such as encryption often require the explicit transmission of data to trustworthy parties, thereby inviting data leakage concerns. The ideal collaboration scenario avoids the explicit dissemination of data/code while maintaining experimental fidelity, which is currently accomplished using various techniques such as trusted execution environments, homomorphic encryption, differential privacy, and multimatrix masking. These techniques, however, often necessitate a trade-off between trust, efficiency, and utility. This article extends a previously proposed technique called the directed infusion of data (DIOD) that ensures data privacy, allows for scalable obfuscation, and combats the risk of data leakage without compromising utility. The experiments discussed in this article examine a regression-type scenario using DIOD with the goal of preserving the inferential link between two variables. Using the point-kinetics equations, regression experiments compare the performance of a model trained using the original data to that of a model trained using the obfuscated data, which produced identical results. Our claim is further strengthened by an information theoretic proof and experiment, which showed that the inferential content between variables remains the same after obfuscation, thereby avoiding the required communication of the proprietary data.

47 - OTHER INSTRUMENTATION↗

Intra-annual variations of regional total column ozone, aerosol optical depth, and water vapor from ground-based, satellite-based and model-based observations

The Microtops II ozonometer has been deployed to continuously measure the total column ozone (TCO), pre-cipitable water content (PWC) and aerosol optical thickness (AOT) at Atigre village (16.74° N, 74.37° E, 604 m above sea level, masl) located on the southeastern slope of Indian western Ghats. The Microtops II ozonometer measurements during the period from September 2017 to June 2018 were analyzed along with the retrieved products of Ozone Monitoring Instrument (OMI) and Moderate Resolution Imaging Spectroradiometer (MODIS) onboard NASA's Terra satellites and also the ERA-interim reanalysis model. We found that the TCO (or PWC) data from the OMI (or MODIS) and the ERA reanalysis model data products are in good agreement. The mea-surements showed strong seasonal variability such that TCO was lowest in the winter season but increased in the summer season, and both PWC and AOT were lowest in the post-monsoon season and increased in the summer season. In continuous ERA-interim reanalysis TCO observations, we found the three types of periodicities in the daily mean TCO as well as PWC (i) The weekly (7 days), (ii) quasi-biweekly (14 to 16 days), (iii) Madden Julian Oscillations (MJO) (30–60 days) oscillations. These oscillations depend on the periodic weather changes in the troposphere. We also observe that the average day to day variability in daily mean TCO observations in the observational period is about 1.4%, 1.1%, 2.62% in the Microtops II ozonometer, ERA-interim reanalysis, and OMI data respectively. We found the positive correlation between daily mean AOT and PWC in all seasons (except monsoon), which is maximum in the winter season, showing the hygroscopic nature of aerosols. The sources of water vapor and aerosol at our location are also studied using back-trajectory analysis.

Nade, Dada R.↗

A Review of Recent and Emerging Machine Learning Applications for Climate Variability and Weather Phenomena

Abstract Climate variability and weather phenomena can cause extremes and pose significant risk to society and ecosystems, making continued advances in our physical understanding of such events of utmost importance for regional and global security. Advances in machine learning (ML) have been leveraged for applications in climate variability and weather, empowering scientists to approach questions using big data in new ways. Growing interest across the scientific community in these areas has motivated coordination between the physical and computer science disciplines to further advance the state of the science and tackle pressing challenges. During a recently held workshop that had participants across academia, private industry, and research laboratories, it became clear that a comprehensive review of recent and emerging ML applications for climate variability and weather phenomena that can cause extremes was needed. This article aims to fulfill this need by discussing recent advances, challenges, and research priorities in the following topics: sources of predictability for modes of climate variability, feature detection, extreme weather and climate prediction and precursors, observation–model integration, downscaling, and bias correction. This article provides a review for domain scientists seeking to incorporate ML into their research. It also provides a review for those with some ML experience seeking to broaden their knowledge of ML applications for climate variability and weather.

54 ENVIRONMENTAL SCIENCES↗

Monthly Mean In Situ Surface Flux Observations Paired with Satellite-Derived and Reanalysis-Based Flux Data for the Great Lakes Region, 2001–2020

Surface radiative and turbulent heat fluxes over the Great Lakes strongly influence regional hydrological and meteorological processes, and their accurate representation is critical for numerical weather prediction and coupled atmosphere–lake modeling. However, direct flux observations are spatially sparse across the region, so gridded reanalysis and satellite-derived products are often used for climatological analyses and model evaluation despite differences in their flux representations. This dataset provides processed, quality-controlled, monthly mean surface flux observations from the Great Lakes Evaporation Network (GLEN), AmeriFlux, and the National Data Buoy Center, paired with spatiotemporally matched flux estimates from two reanalysis products, the fifth generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis dataset (ERA5) and the Modern Era Reanalysis for Research and Applications, version 2 (MERRA-2), and two satellite-derived products, the Clouds and Earth's Radiant Energy Systems Energy Balanced and Filled (CERES-EBAF) and the Cloud, Albedo and Surface Radiation dataset from AVHRR data - Edition 3 (CLARA-A3). The dataset includes sixteen observational stations with variable temporal coverage within 2001–2020. For each station, a CSV file contains monthly time series of available flux variables, including surface downwelling shortwave radiation (SW), surface downwelling longwave radiation (LW), sensible heat (SH) flux, and latent heat flux (LH), alongside matched gridded product values where available. Columns in the CSV file correspond to different variables sourced from each dataset, with column titles structured as "{dataset}_{variable}". Columns with relevant metadata are also provided in each CSV file, including station latitude and longitude, monthly timestamps, and the name of the sourced observational data. These files are structured for direct use in common analysis tools, including Microsoft Excel, Python pandas, and Python matplotlib. This dataset supports climatological analysis of the Great Lakes regional surface energy budget, evaluation of satellite-derived and reanalysis-based flux products, and development or validation of flux representations in numerical weather prediction and coupled atmosphere–lake models.

Great Lakes↗

Symmetry-mode analysis for local structure investigations using pair distribution function data

Symmetry-adapted distortion modes provide a natural way of describing distorted structures derived from higher-symmetry parent phases. Structural refinements using symmetry-mode amplitudes as fit variables have been used for at least ten years in Rietveld refinements of the average crystal structure from diffraction data; more recently, this approach has also been used for investigations of the local structure using real-space pair distribution function (PDF) data. Here, the value of performing symmetry-mode fits to PDF data is further demonstrated through the successful application of this method to two topical materials: TiSe2, where a subtle but long-range structural distortion driven by the formation of a charge-density wave is detected, and MnTe, where a large but highly localized structural distortion is characterized in terms of symmetry-lowering displacements of the Te atoms. Here, the analysis is performed using fully open-source code within the DiffPy framework via two packages developed for this work: isopydistort, which provides a scriptable interface to the ISODISTORT web application for group theoretical calculations, and isopytools, which converts the ISODISTORT output into a DiffPy-compatible format for subsequent fitting and analysis. These developments expand the potential impact of symmetry-adapted PDF analysis by enabling high-throughput analysis and removing the need for any commercial software.

36 MATERIALS SCIENCE↗

Data assimilation empowered neural network parametrizations for subgrid processes in geophysical flows

In the past couple of years, there has been a proliferation in the use of machine learning approaches to represent subgrid-scale processes in geophysical flows with an aim to improve the forecasting capability and to accelerate numerical simulations of these flows. Despite its success for different types of flow, the online deployment of a data-driven closure model can cause instabilities and biases in modeling the overall effect of subgrid-scale processes, which in turn leads to inaccurate prediction. To tackle this issue, we exploit the data assimilation technique to correct the physics-based model coupled with the neural network as a surrogate for unresolved flow dynamics in multiscale systems. In particular, we use a set of neural network architectures to learn the correlation between resolved flow variables and the parametrizations of unresolved flow dynamics and formulate a data assimilation approach to correct the hybrid model during their online deployment. We illustrate our framework in a set of applications of the multiscale Lorenz 96 system for which the parametrization model for unresolved scales is exactly known, and the two-dimensional Kraichnan turbulence system for which the parametrization model for unresolved scales is not known a priori. Our analysis, therefore, comprises a predictive dynamical core empowered by (i) a data-driven closure model for subgrid-scale processes, (ii) a data assimilation approach for forecast error correction, and (iii) both data-driven closure and data assimilation procedures. We show significant improvement in the long-term prediction of the underlying chaotic dynamics with our framework compared to using only neural network parametrizations for future prediction. Moreover, we demonstrate that these data-driven parametrization models can handle the non-Gaussian statistics of subgrid-scale processes, and effectively improve the accuracy of outer data assimilation workflow loops in a modular nonintrusive way.

42 ENGINEERING↗

Neural chaos: A spectral stochastic neural operator

Building surrogate models for operators with uncertainty quantification capabilities is essential for many engineering applications where randomness–such as variability in material properties, boundary conditions, and initial conditions–is unavoidable. Polynomial Chaos Expansion (PCE) is widely recognized as a go-to method for constructing stochastic surrogates in both intrusive and non-intrusive ways, and it has recently been used in the context of operator learning. However, its application becomes challenging for complex or high-dimensional processes, as achieving accuracy requires higher-order polynomials, which can increase computational demand and/or the risk of overfitting. Furthermore, PCE requires specialized treatments to manage random variables that are not independent, and these treatments may be problem-dependent or may fail with increasing complexity. Here, in this work, we adopt the same formalism as the spectral expansion used in PCE; however, we replace the classical polynomial basis functions with neural network (NN) basis functions to leverage their expressivity. To achieve this, we propose an algorithm that identifies NN-parameterized basis functions in a purely data-driven manner, without any prior assumptions about the joint distribution of the random variables involved, whether independent or dependent, or about their marginal distributions. The proposed algorithm identifies each NN-parameterized basis function sequentially, ensuring they are orthogonal with respect to the data distribution. The basis functions are constructed directly on the joint stochastic variables without requiring a tensor product structure or assuming independence of the random variables. This approach may offer greater flexibility for complex stochastic models, while simplifying implementation compared to the tensor product structures typically used in PCE to handle random vectors. This is particularly advantageous given the current state of open-source packages, where building and training neural networks can be done with just a few lines of code and extensive community support. We demonstrate the effectiveness of the proposed scheme through several numerical examples of varying complexity and provide comparisons with classical PCE.

Polynomial chaos expansion↗

Implementing Superresolution of Nonstationary Tides with Wavelets: An Introduction to CWT_Multi

Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.

Lobo, Matthew↗

Data-driven Detection and Early Prediction of Thermoacoustic Instability in a Multi-nozzle Combustor

Thermoacoustic instability (TAI) is a critical issue in modern lean-burn gas-turbine combustors, which is induced by a strong coupling between the resonant combustor acoustics and fluctuations in the heat release rate. This instability may lead to high-amplitude pressure waves that generate undesirable noise levels as well as fatigue stresses in mechanical structures of the combustor. The intense pressure fluctuations due to TAI may also cause large flow perturbations and possibly flow reversal that may lead to flame oscillations, flame liftoff, and even flame blow-out. Hence, there is a strong need for exercising control actions in a timely fashion to mitigate the TAI phenomena. Anomaly detection is an essential prerequisite to the design of a good controller and such a detector must be able to reliably predict a forthcoming TAI. To detect and predict the onset of a TAI from an ensemble of pressure time series, this paper investigates three data-driven methods: Fast Fourier transform (FFT), symbolic time series analysis (STSA), and hidden Markov modeling (HMM). The main focus of the paper is to make a comparative evaluation of these three anomaly detection methods for classification of the current regime of operation into stable and unstable categories as well as for real-time identification of precursors to impending instabilities with short-length time series of measured variables (e.g., pressure oscillations). Finally, the results, generated on experimental data from a multi-nozzle combustor apparatus, have been compared to evaluate the performance of FFT, STSA, and HMM methods for TAI analysis.

42 ENGINEERING↗

Data-driven wind turbine wake modeling via probabilistic machine learning

Wind farm design primarily depends on the variability of the wind turbine wake flows to the atmospheric wind conditions and the interaction between wakes. Physics-based models that capture the wake flow field with high-fidelity are computationally very expensive to perform layout optimization of wind farms, and, thus, data-driven reduced-order models can represent an efficient alternative for simulating wind farms. In this work, we use real-world light detection and ranging (LiDAR) measurements of wind-turbine wakes to construct predictive surrogate models using machine learning. Specifically, we first demonstrate the use of deep autoencoders to find a low-dimensional latent space that gives a computationally tractable approximation of the wake LiDAR measurements. Then, we learn the mapping between the parameter space and the (latent space) wake flow fields using a deep neural network. Additionally, we also demonstrate the use of a probabilistic machine learning technique, namely, Gaussian process modeling, to learn the parameter-space-latent-space mapping in addition to the epistemic and aleatoric uncertainty in the data. Finally, to cope with training large datasets, we demonstrate the use of variational Gaussian process models that provide a tractable alternative to the conventional Gaussian process models for large datasets. Furthermore, we introduce the use of active learning to adaptively build and improve a conventional Gaussian process model predictive capability. Overall, we find that our approach provides accurate approximations of the wind-turbine wake flow field that can be queried at an orders-of-magnitude cheaper cost than those generated with high-fidelity physics-based simulations.

Deep neural networks↗

Evaluation and development of flow condensation correlations using the data from low GWP refrigerants in an axial micro-fin aluminum tube

To mitigate global warming, the world is transitioning to refrigerants with low global warming potential (GWP). Supporting this shift requires a model that can accurately predict the heat transfer and pressure drop of new refrigerants, crucial for designing efficient heat exchangers. Existing models, however, are largely based on currently deployed refrigerants and primarily developed for unexpanded micro-fin tubes with spiral angles of 6° to 30°. Their applicability to new refrigerants, especially in expanded micro-fin tubes, is uncertain. This study assesses the performance of four well-known condensation models for six emerging refrigerants—R-32, R-454B, R-454C, R-455A, R-1234yf, and R-1234ze(E)—against experimental data. Initially, the Han and Lee (2005) model shows the best prediction accuracy with a mean absolute deviation (MAD) of 22.1 %. To enhance the accuracy of heat transfer models for new refrigerants and geometries with large temperature glides, two approaches are proposed. Here, the first approach applies a simple correction factor, reducing the MAD of the Cavallini et al. (2009) model from 68.2 % to 15.4 %. The second approach uses the variable metric method for minimization, fitting new constants to the data. This optimization results in the Kedzierski and Goncalves (1997) model achieving the highest accuracy, with a MAD of 13.1 %. For pressure drop models, the Cavallini et al. (1997) model is the most accurate with a MAD of 6.4 %, followed by the Haraguchi et al. (1993) model with a MAD of 9.4 %. Due to its simplicity, the Haraguchi et al. (1993) model is a practical option for predicting frictional pressure drop.

Axial micro-fin tubes↗

Analyzing Wildland Fire Smoke Emissions Data Using Compositional Data Techniques

By conservation of mass, the mass of wildland fuel that is pyrolyzed and combusted must equal the mass of smoke emissions, residual char and ash. For a given set of conditions, these amounts are fixed. This places a constraint on smoke emissions data which violates statistical assumptions for many of the methods currently used to analyze these data such as linear regression, analysis of variance, and t-tests. These data are inherently multivariate and non-negative parts of a whole. This paper introduces the field of compositional data analysis to the emissions community and provides examples of appropriate statistical treatment of emissions data. It is shown that modified combustion efficiency should not be used as a predictor variable for other smoke emissions because it is not an independent variable. An alternative method based on compositional linear trends to estimate trace gas composition using CO and CO2 is presented. The data used in this paper resulted from projects the DOD/DOE/EPA Strategic 586 Environmental Research and Development Program projects RC-1648 and 1649. The senior 587 author appreciates the guidance and R scripts provided by Prof. Girty at San Diego State 588 University to estimate linear trends by perturbation. J. P.-A. was supported by the Spanish 589 Ministry of Science, Innovation and Universities under the project CODAMET (RTI2018-590 095518-B-C21, 2019-2021). The data used in this study have been previously published and are 591 available in the original publications. DRW conceived the initial manuscript (70 percent) and 592 performed the bulk of the data analysis. JPA provided statistical guidance and compositional data 593 expertise and contributed 20 percent of the manuscript. TJJ and HJ were extensively involved in 594 the study that provided the data. TJJ provide smoke emissions expertise and HJ provided 595 combustion expertise. The authors declare that they have no conflict of interest. The use of trade 596 or firm names in this publication is for reader information and does not imply endorsement by 597 the U.S. Department of Agriculture of any product or service.

simplex, compositional data analysis, balance, log↗

On the Emergence of Human Influence on Surface Air Temperature Changes Over India

Human activities in terms of greenhouse gas (GHG) emissions and aerosols resulting from the combustion of fossil fuels have been shown to have affected the temperature of the Earth on global and continental scales. The surface air temperature (TAS) over India has also been observed to be increasing over the last 100 years. Understanding the underlying causes of regional climate change over India can help in developing appropriate mitigation and adaptation strategies. Differentiating signals of externally forced climate changes from the noise of natural internal variability generally becomes more difficult as spatial scale reduces. Therefore detecting and attributing the influence of external forcings such as GHGs and aerosols is harder at local and regional scales. In this study, we applied a detection and attribution (D&A) method to study annual and seasonal mean TAS over the Indian region. We found that the observed warming over India from 1906 to 2005 cannot be explained by natural climate variability alone. We found that the warming is largely driven by the increase in GHGs, and partially offset by regional anthropogenic emissions of aerosols. These results were confirmed for the shorter 1956–2005 period, but results were sensitive to the choice of observational data set. Furthermore, the changes cannot be explained by internal climate variability or natural external forcings alone, but are compatible with the responses to combined anthropogenic GHG and aerosol forcings.

54 ENVIRONMENTAL SCIENCES↗

Assessment of rainfall variability and future change in Brazil across multiple timescales

Abstract Rainfall variability change under global warming is a crucial issue that may have a substantial impact on society and the environment, as it can directly impact biodiversity, agriculture, and water resources. Observed precipitation trends and climate change projections over Brazil indicate that many sectors of society are potentially highly vulnerable to the impacts of climate change. The purpose of this study is to assess model projections of the change in rainfall variability at various temporal scales over sub‐regions of Brazil. For this, daily data from 30 CMIP5 models for historical (1900–2005) and future (2050–2100) experiments under a high‐emission scenario are used. We assess the change in precipitation variability, applying a band‐pass filter to isolate variability on daily, weekly, monthly, intra‐seasonal, and El Nino Southern Oscillation (ENSO) time scales. For historical climate, simulated precipitation is evaluated against observations to establish model reliability. The results show that models largely agree on increases in variability on all timescales in all sub‐regions, except on ENSO timescales where models do not agree on the sign of future change. Brazil will experience more rainfall variability in the future that is, drier or more frequent dry periods and wetter wet periods on daily, weekly, monthly, and intra‐seasonal timescales, even in sub‐regions where future changes in mean rainfall are currently uncertain. This may provide useful information for climate change adaptation across, for example, the agriculture and water resource sectors in Brazil.

Alves, Lincoln M.↗

Machine learning-enabled model-data integration for predicting subsurface water storage

Subsurface water storage (SWS) is a key variable of the climate system and a storage component for precipitation and radiation anomalies, inducing persistence in the climate system. It plays a critical role in climate-change projections and can mitigate the impacts of climate change on ecosystems. However, because of the difficult accessibility of the underground, hydrologic properties and dynamics of SWS are poorly known. Direct observations of SWS are limited, and accurate incorporation of SWS dynamics into Earth system land models remains challenging. We propose a machine learning-enabled model-data integration framework to improve the SWS prediction at local to conus scales in a changing climate by leveraging all the available observation and simulation resources, as well as to inform the model development and guide the observation collection. The accurate prediction will enable an optimal decision of water management and land use and improve the ecosystem's resilience to the climate change.

Lu, Dan↗

The MOSAiC Distributed Network: Observing the coupled Arctic system with multidisciplinary, coordinated platforms

Central Arctic properties and processes are important to the regional and global coupled climate system. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Distributed Network (DN) of autonomous ice-tethered systems aimed to bridge gaps in our understanding of temporal and spatial scales, in particular with respect to the resolution of Earth system models. By characterizing variability around local measurements made at a Central Observatory, the DN covers both the coupled system interactions involving the ocean-ice-atmosphere interfaces as well as three-dimensional processes in the ocean, sea ice, and atmosphere. The more than 200 autonomous instruments (“buoys”) were of varying complexity and set up at different sites mostly within 50 km of the Central Observatory. During an exemplary midwinter month, the DN observations captured the spatial variability of atmospheric processes on sub-monthly time scales, but less so for monthly means. They show significant variability in snow depth and ice thickness, and provide a temporally and spatially resolved characterization of ice motion and deformation, showing coherency at the DN scale but less at smaller spatial scales. Ocean data show the background gradient across the DN as well as spatially dependent time variability due to local mixed layer sub-mesoscale and mesoscale processes, influenced by a variable ice cover. The second case (May–June 2020) illustrates the utility of the DN during the absence of manually obtained data by providing continuity of physical and biological observations during this key transitional period. We show examples of synergies between the extensive MOSAiC remote sensing observations and numerical modeling, such as estimating the skill of ice drift forecasts and evaluating coupled system modeling. The MOSAiC DN has been proven to enable analysis of local to mesoscale processes in the coupled atmosphere-ice-ocean system and has the potential to improve model parameterizations of important, unresolved processes in the future.

54 ENVIRONMENTAL SCIENCES↗

Data associated with “Different methods of estimating riverbed sediment grain size diverge at the basin scale ” (v2)

This data package is associated with the publication “Different methods of estimating riverbed sediment grain size diverge at the basin scale” published in Frontiers in Earth Science (Regier et al., 2025). The distribution of sediment grain size in streams and rivers is often quantified by the median grain size (d50), a key metric for understanding and predicting hydrologic and biogeochemical function of streams and rivers. Manual methods to measure d50 are time-consuming and ignore larger grains, while model-based methods to estimate d50 often over-generalize basin characteristics, and therefore cannot accurately represent site-scale heterogeneity. Here, we apply a machine learning-enabled photogrammetry methodology (You Only Look Once, or YOLO) for estimating d50 for grains > 2 mm based on images collected from streams and rivers throughout the Yakima River Basin (YRB). To understand how such methods may help bridge the gaps in resolution and accuracy between manual and catchment characteristics model-based d50 estimates, we compared YOLO d50 values to manual and model-based estimates across the YRB. We found distinct differences among methods for d50 averages and variability, and relationships between d50 estimates and basin characteristics. Source images can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1892052. This data package was originally published in May 2023. It was updated August 2025 (v2; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. In addition to the 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. This dataset is comprised of one main data folder containing (1) file-level metadata; (2) data dictionary; (3) readme; (4) and subfolders containing data, figures, and scripts. The data folder contains datasets used for the analyses in the manuscript in image, text-delimited or geospatially-referenced formats. The figures folder contains the figures from the manuscript in different formats. The scripts folder contains all of the scripts used to complete the analyses in the manuscript. All files are .csv, .rds, .dbf, .prj, .shp, .shx, .jpg, .png, .R, .Rproj, or .pdf. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the influence of natural climate variability on in situ measurements of seasonal total and extreme daily precipitation

Abstract While various studies explore the relationship between individual sources of climate variability and extreme precipitation, there is a need for improved understanding of how these physical phenomena simultaneously influence precipitation in the observational record across the contiguous United States. In this work, we introduce a single framework for characterizing the historical signal (anthropogenic forcing) and noise (natural variability) in seasonal mean and extreme precipitation. An important aspect of our analysis is that we simultaneously isolate the individual effects of seven modes of variability while explicitly controlling for joint inter-mode relationships. Our method utilizes a spatial statistical component that uses in situ measurements to resolve relationships to their native scales; furthermore, we use a data-driven procedure to robustly determine statistical significance. In Part I of this work we focus on natural climate variability: detection is mostly limited to DJF and SON for the modes of variability considered, with the El Niño/Southern Oscillation, the Pacific–North American pattern, and the North Atlantic Oscillation exhibiting the largest influence. Across all climate indices considered, the signals are larger and can be detected more clearly for seasonal total versus extreme precipitation. We are able to detect at least some significant relationships in all seasons in spite of extremely large (> 95%) background variability in both mean and extreme precipitation. Furthermore, we specifically quantify how the spatial aspect of our analysis reduces uncertainty and increases detection of statistical significance while also discovering results that quantify the complex interconnected relationships between climate drivers and seasonal precipitation.

54 ENVIRONMENTAL SCIENCES↗