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

Coupling Remote Sensing With a Process Model for the Simulation of Rangeland Carbon Dynamics

Rangelands provide significant environmental benefits through many ecosystem services, which may include soil organic carbon (SOC) sequestration. However, quantifying SOC stocks and monitoring carbon (C) fluxes in rangelands are challenging due to the considerable spatial and temporal variability tied to rangeland C dynamics as well as limited data availability. We developed the Rangeland Carbon Tracking and Management (RCTM) system to track long-term changes in SOC and ecosystem C fluxes by leveraging remote sensing inputs and environmental variable data sets with algorithms representing terrestrial C-cycle processes. Bayesian calibration was conducted using quality-controlled C flux data sets obtained from 61 Ameriflux and NEON flux tower sites from Western and Midwestern US rangelands to parameterize the model according to dominant vegetation classes (perennial and/or annual grass, grass-shrub mixture, and grass-tree mixture). The resulting RCTM system produced higher model accuracy for estimating annual cumulative gross primary productivity (GPP) (R 2 > 0.6, RMSE <390 g C m -2 ) relative to net ecosystem exchange of CO 2 (NEE) (R 2 > 0.4, RMSE <180 g C m -2 ). Model performance in estimating rangeland C fluxes varied by season and vegetation type. The RCTM captured the spatial variability of SOC stocks with R 2 = 0.6 when validated against SOC measurements across 13 NEON sites. Model simulations indicated slightly enhanced SOC stocks for the flux tower sites during the past decade, which is mainly driven by an increase in precipitation. Future efforts to refine the RCTM system will benefit from long-term network-based monitoring of vegetation biomass, C fluxes, and SOC stocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Hyperspectral Reflectance for Measuring Canopy‐Level Nutrients and Photosynthesis in a Salt Marsh

Abstract Salt marsh ecosystems are underrepresented in process‐based models due to their unique location across the terrestrial–aquatic interface. Particularly, the role of leaf nutrients on canopy photosynthesis (F A ) remains unclear, despite their relevance for regulating vegetation growth. We combined multiyear information of canopy‐level nutrients and eddy covariance measurements with canopy surface hyperspectral remote sensing (CSHRS) to quantify the spatial and temporal variability of F A in a temperate salt marsh. We found that F A showed a positive relationship with canopy‐level N at the ecosystem scale and for areas dominated by Spartina cynosuroides , but not for areas dominated by short S. alterniflora . F A showed a positive relationship with canopy‐level P, K, and Na, but a negative relationship with Fe, for areas associated with S. cynosuroides , S. alterniflora , and at the ecosystem scale. We used partial least squares regression (PLSR) with CSHRS and found statistically significant data–model agreements to predict canopy‐level nutrients and F A . The red‐edge electromagnetic region and ∼770 nm showed the highest contribution of variance in PLSR models for canopy‐level nutrients and F A , but we propose that underlying sediment biogeochemistry can complicate interpretation of reflectance measurements. Our findings highlight the relevance of spatial variability in salt marshes vegetation and the promising application of CSHRS for linking information of canopy‐level nutrients with F A . We call for further development of canopy surface hyperspectral methods and analyses across salt marshes to improve our understanding of how these ecosystems will respond to global environmental change.

Environmental Sciences & Ecology↗
DOE OSTI2021

Spatial–Temporal Patterns of Historical, Near-Term, and Projected Drought in the Conterminous United States

Major droughts in the United States have heavily impacted the hydrologic system, negatively effecting energy and food production. Improved understanding of historical drought is critical for accurate forecasts. Data from global climate models (GCMs), commonly used to assess drought, cannot effectively evaluate local patterns because of their low spatial scale. This research leverages downscaled (~4 km grid spacing) temperature and precipitation estimates from nine GCMs’ data under the business-as-usual scenario (Representative Concentration Pathway 8.5) to examine drought patterns. Drought severity is estimated using the Palmer Drought Severity Index (PDSI) with the Thornthwaite evapotranspiration method. The specific objectives were (1) To reproduce historical (1966–2005) drought and calculate near-term to future (2011–2050) drought patterns over the conterminous USA. (2) To uncover the local variability of spatial drought patterns in California between 2012 and 2018 using a network-based approach. Our estimates of land proportions affected by drought agree with the known historical drought events of the mid-1960s, late 1970s to early 1980s, early 2000s, and between 2012 and 2015. Network analysis showed heterogeneity in spatial drought patterns in California, indicating local variability of drought occurrence. The high spatial scale at which the analysis was performed allowed us to uncover significant local differences in drought patterns. This is critical for highlighting possible weak systems that could inform adaptation strategies such as in the energy and agricultural sectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear-and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute time-scale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

irradiance↗
DOE OSTI2021

Simulation of PV Variability as a Function of PV Generation and Plant Size: Preprint

The deployment of photovoltaic (PV) systems continues to show significant expansion; however, this growth has brought added attention to issues around the variability of the solar resource. Both spatial and temporal variability exist. Temporal scales can range from the sub-second to multiyear, whereas spatial scales can range from a few meters to tens of kilometers. There are multiple methods described in the literature to quantify PV variability at various spatial and temporal scales. This study focuses on short-term temporal variability and uses similar approaches with the addition of PV plant size a parameter to quantify variability. The method employed here incorporates the normalization of clear- and cloudy-sky conditions and PV plant size to quantify nominal variability metrics. The distribution and fluctuations of these metrics provide relevant information that is useful for system operations. The National Solar Radiation Database (NSRDB) is used to simulate PV variability as a function of PV generation and plant size. Hypothetical but realistic system information at 33 locations is used to model PV generation by feeding NSRDB solar irradiance data to the National Renewable Energy Laboratory’s System Advisor Model (SAM). Over the selected region, it is found that the aggregated ramp rates for the 1-minute data are associated with standard deviations ranging from 0.002–0.055 on a daily basis; however, hourly intervals induce higher aggregated ramp rates than the other timescales. Even though minute-to-minute variations are significant for the 1-minute timescale, the standard deviation aggregated into a daily metric is smaller because of the cancellation of values.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

An Automated Method for Mapping Giant Kelp Canopy Dynamics from UAV

Satellite and aerial imagery have been used extensively for mapping the abundance and distribution of giant kelp ( Macrocystis pyrifera ) in southern California. There is now great potential for using unoccupied aerial vehicles (UAVs) to map kelp canopy at very high resolutions. However, tides and currents have been shown to affect the amount of floating kelp canopy on the water surface, and the impacts of these processes on remotely sensed kelp estimates in this region have not been fully quantified. UAVs were used to map fine-scale changes in canopy area due to tidal height and current speed at kelp forests off the coast of Palos Verdes, CA and Santa Barbara, CA. An automated method for detecting kelp canopy was developed that was 67% accurate using red-green-blue (RGB) UAV imagery and 93% accurate using multispectral UAV imagery across a range of weather, ocean, and illumination conditions. Increases in tidal height of 1 m reduced the amount of floating kelp canopy by 15% in Santa Barbara and by over 30% in Palos Verdes. The effect of current speed on visible kelp canopy was inconclusive, but there was a trend towards lower canopy area with increased current speed. Therefore, while tidal height and current speed can introduce significant variability to estimates of kelp abundance, the magnitude of this variability is site specific. Still, UAVs are a valuable tool for mapping of kelp canopy and can provide greater spatial resolution and temporal coverage than is possible from many satellite sensors. This data can provide insight into the patterns and drivers of high frequency fluctuations in kelp abundance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Field Campaign Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) atmospheric observatory. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120-km 2 domain in northern and central Oklahoma. The field campaign targeted summertime convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the ARM SGP site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds that penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short period of 11 days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 that resulted in variability of 10/7Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height, and depth information. We also ran Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is ongoing. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data. Lastly, we note that BIRDIE coincided with two other deep convection field campaigns in the area: The National Aeronautics and Space Administration (NASA)’s Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) experiment (https://dcotss.org/) centered in Salina, Kansas, and the Targeted Observation by Radars and (uncrewed aerial systems) UAS of Supercells (TORUS) experiment (https://catalog.eol.ucar.edu/torus_2022) across the U.S. Central Plains. Additionally, a Lawrence Livermore National Laboratory ( LLNL) team at DCOTSS attempted to measure 10Be/7Be ratios on aerosols collected in situ in the lower and mid-stratosphere over Kansas using high-altitude balloons. Precipitation samples were also collected on the ground in Salina, Kansas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Final Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the DOE Southern Great Plains (SGP) Atmospheric Radiation Measurement (ARM) site. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120 km 2 domain in northern and central Oklahoma. The field campaign targeted summer-time convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the SGP ARM site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds which penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short time period of eleven days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 which resulted in variability of 10/7 Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height and depth information. We also ran Hysplit back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is still on-going. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data.

58 GEOSCIENCES↗
DOE OSTI2023

Machine Learning of Key Variables Impacting Extreme Precipitation in Various Regions of the Contiguous United States

Abstract Amplification in extreme precipitation intensity and frequency can cause severe flooding and impose significant social and economic consequences. Variations in extreme precipitation intensity, frequencies, and return periods can be attributed to many physical variables across spatial and temporal scales. Here we employ ensemble machine learning (ML) methods, namely random forest (RF), eXtreme Gradient Boosting (XGB), and artificial neural networks (ANN), to explore key contributing variables to monthly extreme precipitation intensity and frequency in six regions over the United States. We further establish emulators for return periods. Results show that the ML models for intensity perform better in regions with obvious seasonality (i.e., Northern Great Plains, Southern Great Plains, and West Coast) than the other three regions (Northeast, Southwest, and Rocky Mountains), while for frequency the models perform well for most regions. The Shapley additive explanation is used to help explain the relationships between extreme precipitation characteristics and identify top variables for RF and XGB. We find that latent heat flux, relative humidity, soil moisture, and large‐scale subsidence are key common variables across the regions for both monthly intensity and frequency, and their compound effects are non‐negligible. The developed ML models capture the probability and return period of extreme precipitation well for all regions and may be used for decision making (e.g., infrastructure planning and design).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High resolution variability in wet deposition in the southeastern United States

Rainwater chemistry is determined by atmospheric pollutants and particles which vary spatially and temporally. Industrial and agricultural activities and meteorological events (e.g. sea breezes, severe weather, blowing dust) alter atmospheric particle and trace gas compositions. These gases and particles are scavenged by cloud and rain droplets that drive wet deposition. During an Intensive Operation Period (IOP) from April to October 2021, rainwater was collected at higher frequency intervals, usually daily, after precipitation events at three locations on the Savannah River Site (SRS). The farthest locations were separated by approximately 20 km. The mean concentration (μeq/L) of seven ions followed the Cl⁻ > SO 4 2− > Na⁺ > NO 3 ⁻ > K⁺ > Mg 2+ > Ca 2+ downward trend. Ion concentrations were compared to background ion concentrations from the National Atmospheric Deposition Program (NADP). The high frequency monthly averaged SRS data compared well with the monthly averaged NADP background but demonstrated extensive variability. In some months in 2021, the high frequency data compared better to the NADP site near the coast while in other months inland sites compared better. Strong spatial variability for ion concentrations was observed across SRS which was attributed to localized impacts in rainfall spatial variability. High frequency measurements allowed for comparison to regional weather patterns indicating influences from the Atlantic Ocean, Gulf of Mexico, and cities. This can account for spatial variability in the wet deposition flux. Sea breezes, Saharan dust, and anthropogenic sources were shown to impact wet deposition flux variability. Higher frequency precipitation chemistry sampling at numerous locations better captures ion concentration variability and improves measurement representativeness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

The biogeography of soil and airborne fungi in the Southwestern USA in relation to climate and vegetation

To assess how fungal dispersal might respond to climate change, we examined how climate and geography influence the regional distribution of fungi in soil and air. Specifically, we hypothesized that neighboring fungal communities should be more similar than distant communities (i.e. spatially autocorrelated) and that fungal dispersal should be more limited in soil than in air. We collected soil and air samples from 60 sites across five states in the Southwestern USA. Then, we sequenced the ITS2 region to identify fungal taxa in each sample. Next, we used distance-based redundancy analysis to partition variation in fungal community composition between climate variables and spatial structure. Fungi were indeed spatially autocorrelated. Moreover, precipitation, maximum vapor pressure deficit, and soil moisture were significantly related to fungal community composition in soils. In comparison, only precipitation was significantly related to community composition in the air. After accounting for climate, the strength of spatial autocorrelation did not differ significantly in soilborne versus airborne fungi. Dispersal limitation was evident in soilborne fungi at short distances (<100 km) and was not observed at any distance in airborne fungi. Altogether, climate may influence which fungal taxa are present in soil and air, and fungi could feasibly wind disperse over regional scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Data from: Understanding the biogeochemical and spatial drivers of methane and carbon dioxide fluxes in a large temperate reservoir

This dataset contains spatially resolved measurements of CO₂ and CH₄ fluxes and associated environmental variables collected across 200 sites in Douglas Reservoir (Tennessee, USA) between July 29-August 2, 2024. Measurements include diffusive fluxes of CO₂ and CH₄, CH₄ ebullition, and biogeochemical and spatial variables such as dissolved oxygen, temperature, conductivity, chlorophyll-a, pH, water depth, and distance from the dam. Sampling was conducted using a spatially balanced design to capture longitudinal and depth-related gradients throughout the reservoir. The dataset is structured to support analyses of spatial variability, flux pathway comparisons, and modeling approaches (e.g., spatial statistics and machine learning) aimed at understanding controls on reservoir CO₂ and CH₄ fluxes and improving upscaling to whole-reservoir and regional estimates.

Neeper, Jamie [ORNL] (ORCID:0009000842342101)↗
DOE OSTI2023

Robust global detection of forced changes in mean and extreme precipitation despite observational disagreement on the magnitude of change

Detection and attribution (D&A) of forced precipitation change are challenging due to internal variability, limited spatial, and temporal coverage of observational records and model uncertainty. These factors result in a low signal-to-noise ratio of potential regional and even global trends. Here, we use a statistical method – ridge regression – to create physically interpretable fingerprints for the detection of forced changes in mean and extreme precipitation with a high signal-to-noise ratio. The fingerprints are constructed using Coupled Model Intercomparison Project phase 6 (CMIP6) multi-model output masked to match coverage of three gridded precipitation observational datasets – GHCNDEX, HadEX3, and GPCC – and are then applied to these observational datasets to assess the degree of forced change detectable in the real-world climate in the period 1951–2020. We show that the signature of forced change is detected in all three observational datasets for global metrics of mean and extreme precipitation. Forced changes are still detectable from changes in the spatial patterns of precipitation even if the global mean trend is removed from the data. This shows the detection of forced change in mean and extreme precipitation beyond a global mean trend is robust and increases confidence in the detection method's power as well as in climate models' ability to capture the relevant processes that contribute to large-scale patterns of change. We also find, however, that detectability depends on the observational dataset used. Not only coverage differences but also observational uncertainty contribute to dataset disagreement, exemplified by the times of emergence of forced change from internal variability ranging from 1998 to 2004 among datasets. Furthermore, different choices for the period over which the forced trend is computed result in different levels of agreement between observations and model projections. These sensitivities may explain apparent contradictions in recent studies on whether models under- or overestimate the observed forced increase in mean and extreme precipitation. Lastly, the detection fingerprints are found to rely primarily on the signal in the extratropical Northern Hemisphere, which is at least partly due to observational coverage but potentially also due to the presence of a more robust signal in the Northern Hemisphere in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

MSW Variability Mapping and Conversion to Biofuel

MSW (Municipal Solid Waste) is a form of biomass which consists of categorized components of waste/trash. The general categories are paper, yard trash, construction & debris, appliances, tires, glass, metals, aluminum & steel cans, plastics, organics, inorganics, and HHW (Household Hazardous Waste). This project focuses on the factors within a region or population that contribute to variability in the composition of MSW and in turn MSW’s convertibility to biofuel. A list of contributors was determined (Social Vulnerability Index, Access to Public Transportation, Racial Distribution, GDP, Personal Income) and then JMP was used to perform a Multivariate analysis to determine correlations and a Partial Least-Squares regression to determine Variable Importance Plots for each MSW category. In addition to data analysis, the convertibility of MSW to biofuel was studied via microwave pyrolysis system in order to separate and characterize the various gaseous and bio-oil products.

09 BIOMASS FUELS↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers, part II: Regional simulations of building response

The existing observational database of the regional-scale distribution of strong ground motions and measured building response for major earthquakes continues to be quite sparse. As a result, details of the regional variability and spatial distribution of ground motions, and the corresponding distribution of risk to buildings and other infrastructure, are not comprehensively understood. Utilizing high-performance computing platforms, emerging high-resolution, physics-based ground motion simulations can now resolve frequencies of engineering interest and provide detailed synthetic ground motions at high spatial density. This provides an opportunity for new insight into the distribution of infrastructure seismic demands and risk. In the work presented herein, the EQSIM fault-to-structure computational framework described in a companion paper, McCallen et al., is employed to investigate the regional-scale response of buildings to large earthquakes. A representative M = 7.0 strike-slip event is used to explore the distribution and amplitude of building demand, and comparisons are made between building response computed with fault-to-structure simulations and building response computed with existing measured near-fault earthquake records. New information on the distribution and variability of building response from high-performance parallel simulations is described and analyzed, and favorable first comparisons between building response predicted with both fault-to-structure simulations and real ground motions records are presented.

58 GEOSCIENCES↗
DOE OSTI2024

Printed Potentiometric Ammonium Sensors for Agriculture Applications

Ammonium (NH 4 + ) concentration is critical to both nutrient availability and nitrogen (N) loss in soil ecosystems but can be highly variable across spatial and temporal scales. For this reason, effectively informing agricultural practices such as fertilizer management and understanding of mechanisms of soil N loss require sensor technologies to monitor ammonium concentrations in real time. Our work investigates the performance of fully printed ammonium ion-selective sensors used in diverse soil environments. Ammonium sensors consisting of a printed ammonium ion-selective electrode and a printed Ag/AgCl reference were fabricated and characterized in aqueous solutions and three different soil types (sand, peat, and clay) under the range of ion concentrations likely to be present in soil (0.01–100 mM). The response of ammonium sensors was further evaluated under variable gravimetric moisture content in the soil to reflect their reliability under field conditions. Ammonium sensors demonstrated a sensitivity of 53.6 ± 5.1 mV/decade when tested in aqueous solution, and a sensitivity of 55.7 ± 11 mV/dec, 57.5 ± 4.1 mV/dec, and 43.7 ± 4 mV/dec was measured in sand, clay, and peat soils, respectively.

60 APPLIED LIFE SCIENCES↗