Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vegetation model”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Time Series Vegetation Aerodynamic Roughness Fields Estimated from MODIS Observations

Most land surface models used today require estimates of aerodynamic roughness length in order to characterize momentum transfer between the surface and atmosphere. The most common method of prescribing roughness is through the use of empirical look-up tables based solely on land cover class. Theoretical approaches that employ satellite-based estimates of canopy density present an attractive alternative to current look-up table approaches based on vegetation cover type that do not account for within-class variability and are oftentimes simplistic with respect to temporal variability. The current research applies Raupach s formulation of momentum aerodynamic roughness to MODIS data on a regional scale in order to estimate seasonally variable roughness and zero-plane displacement height fields using bulk land cover parameters estimated by [Jasinski, M.F., Borak, J., Crago, R., 2005. Bulk surface momentum parameters for satellite-derived vegetation fields. Agric. For. Meteorol. 133, 55-68]. Results indicate promising advances over look-up approaches with respect to characterization of vegetation roughness variability in land surface and atmospheric circulation models.

Borak, Jordan S.↗
NASA NTRS2015

Advances in Remote Sensing for Vegetation Dynamics and Agricultural Management

Spaceborne remote sensing has led to great advances in the global monitoring of vegetation. For example, the NASA Global Inventory Modeling and Mapping Studies (GIMMS) group has developed widely used datasets from the Advanced Very High Resolution Radiometer (AVHRR) sensors as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) map imagery and normalized difference vegetation index datasets. These data are valuable for analyzing vegetation trends and variability at the regional and global levels. Numerous studies have investigated such trends and variability for both natural vegetation (e.g., re-greening of the Sahel, shifts in the Eurasian boreal forest, Amazonian drought sensitivity) and crops (e.g., impacts of extremes on agricultural production). Here, a critical overview is presented on recent developments and opportunities in the use of remote sensing for monitoring vegetation and crop dynamics.

variability↗
NASA NTRS2023

Maldives Climate II: Evaluating the Potential Impacts of Sea Level Rise on Human Development and Coastal Infrastructure

The Republic of the Maldives is a low-lying island nation in the Indian Ocean which has experienced rapid urbanization, landcover changes, and sea level rise over recent years. The growth of tourism, coastal erosion, and urbanization have all driven land reclamation efforts across many islands. As in-situ landcover change monitoring has proven difficult across the vast archipelago, the NASA DEVELOP team collaborated with the Maldives Ministry of Environment, Climate Change, and Technology; USAID; and the U.S. Department of State, to utilize Earth observations to predict sea level rise impacts on coastal infrastructure. The team used a supervised classification algorithm within Google Earth Engine to create land use maps and time series analyses of nine islands and atolls using imagery from Landsat 7 Enhanced Thematic Mapper Plus, Landsat 8 Operational Land Imager, Sentinel-2 Multispectral Instrument, and PlanetScope, covering a combined period of 2000 through 2023. Additionally, the team projected coastal inundation with a modified deterministic Bathtub model utilizing elevation data from CoastalDEM and 2050-2100 Shared Socioeconomic Pathway scenarios identified in the NASA Sea Level Rise Projection Tool. The team found that islands undergoing urban growth experienced a 23% decrease in vegetation between 2014 and 2022. Furthermore, the model predicted that 51–57% of the study area’s-built environment has a chance of inundation by 2100 under the low and high sea level scenarios. These analyses demonstrate how remote sensing can be used to both track land use changes over time and project how coastlines will be affected by sea level rise.

remote sensing↗
NASA NTRS1983

Monitoring tropical vegetation succession with LANDSAT data

The shadowing problem, which is endemic to the use of LANDSAT in tropical areas, and the ability to model changes over space and through time are problems to be addressed when monitoring tropical vegetation succession. Application of a trend surface analysis model to major land cover classes in a mountainous region of the Phillipines shows that the spatial modeling of radiance values can provide a useful approach to tropical rain forest succession monitoring. Results indicate shadowing effects may be due primarily to local variations in the spectral responses. These variations can be compensated for through the decomposition of the spatial variation in both elevation and MSS data. Using the model to estimate both elevation and spectral terrain surface as a posteriori inputs in the classification process leads to improved classification accuracy for vegetation of cover of this type. Spatial patterns depicted by the MSS data reflect the measurement of responses to spatial processes acting at several scales.

Robinson, V. B.↗
NASA NTRS2022

Interactive Biogenic Emissions and Drought Stress Effects on Atmospheric Composition in NASA GISS ModelE

Drought is a hydroclimatic extreme that causes perturbations to the terrestrial biosphere and acts as a stressor on vegetation, affecting emissions patterns. During severe drought, isoprene emissions are reduced. In this paper, we focus on capturing this reduction signal by implementing a new percentile isoprene drought stress (y d ) algorithm in NASA GISS ModelE based on the MEGAN3 (Model of Emissions of Gases and Aerosols from Nature Version 3) approach as a function of a photosynthetic parameter (V c,max ) and water stress (β). Four global transient simulations from 2003–2013 are used to demonstrate the effect without y d (Default_ModelE) and with online y d (DroughtStress_ModelE). DroughtStress_ModelE is evaluated against the observed isoprene measurements at the Missouri Ozarks AmeriFlux (MOFLUX) site during the 2012 severe drought where improvements in the correlation coefficient indicate it is a suitable drought stress parameterization to capture the reduction signal during severe drought. The application of y d globally leads to a decadal average reduction of ∼ 2.7 %, which is equivalent to ∼14.6 Tg yr -1 of isoprene. The changes have larger impacts in regions such as the southeastern US. DroughtStress_ModelE is validated using the satellite ΩHCHO column from the Ozone Monitoring Instrument (OMI) and surface O 3 observations across regions of the US to examine the effect of drought on atmospheric composition. It was found that the inclusion of isoprene drought stress reduced the overestimation of ΩHCHO in Default_ModelE during the 2007 and 2011 southeastern US droughts and led to improvements in simulated O 3 during drought periods. We conclude that isoprene drought stress should be tuned on a model-by-model basis because the variables used in the parameterization responses are relative to the land surface model hydrology scheme (LSM) and the effects of y d application could be larger than seen here due to ModelE not having large biases of isoprene during severe drought.

drought stress effects↗
NASA NTRS1985

1985 International Geoscience and Remote Sensing Symposium (IGARSS '85), University of Massachusetts, Amherst, October 7-9, 1985, Digest. Volumes 1 & 2

The present conference on remote sensing instrumentation considers topics in water resources research, planetary remote sensing and mathematical geophysics, sea ice behavior in view of the Cold Regions Research and Engineering Laboratories pond measurements, Shuttle Imaging Radar-B (SIR-B) system performance and calibration results, geological applications of remote sensing, and microwave scattering from vegetation. Further attention is given to atmospheric probing with lasers, image processing for remote sensing, forest inventory and condition assessment, atmospheric remote sensing, SIR-B geological results, the reflectance, emission, and scattering characteristics of vegetation, sea ice sensing, multipolarization SAR results, the theoretical modeling of surfaces and volumes, crop condition assessment and productivity estimates, and SIR-B results for vegetation cover. Also considered are the microwave remote sensing of soil moisture, advanced sensors, advanced information extraction, wind and wave remote observations, results from the ESA Remote Sensing satellite, large area crop and land cover statistics, SIR-B results for directional ocean wave spectra, seasonal snow cover, ice sheets, and lake ice, fundamental research in terrain remote sensing, image classification methods, vegetation stress detection, remote sensing of coastal processes, SAR systems and calibration, electromagnetic geophysical methods and signal processing, image processing techniques, and SAR internal wave measurements.

Carver, K. R.↗
NASA NTRS2014

Interannual Variations and Trends in Global Land Surface Phenology Derived from Enhanced Vegetation Index During 1982-2010

Land swiace phenology is widely retrieved from satellite observations at regional and global scales, and its long-term record has been demonstmted to be a valuable tool for reconstructing past climate variations, monitoring the dynamics of terrestrial ecosystems in response to climate impacts, and predicting biological responses to future climate scenarios. This srudy detected global land surface phenology from the advanced very high resolution radiometer (AVHRR) and the Moderate Resolution Imaging Spectroradiometer (MODIS) data from 1982 to 2010. Based on daily enhanced vegetation index at a spatial resolution of 0.05 degrees, we simulated the seasonal vegetative trajectory for each individual pixel using piecewise logistic models, which was then used to detect the onset of greenness increase (OGI) and the length of vegetation growing season (GSL). Further, both overall interannual variations and pixel-based trends were examIned across Koeppen's climate regions for the periods of 1982-1999 and 2000-2010, respectively. The results show that OGI and OSL varied considerably during 1982-2010 across the globe. Generally, the interarmual variation could be more than a month in precipitation-controlled tropical and dry climates while it was mainly less than 15 days in temperature-controlled temperate, cold, and polar climates. OGI, overall, shifted early, and GSL was prolonged from 1982 to 2010 in most climate regions in North America and Asia while the consistently significant trends only occurred in cold climate and polar climate in North America. The overall trends in Europe were generally insignificant. Over South America, late OGI was consistent (particularly from 1982 to 1999) while either positive or negative OSL trends in a climate region were mostly reversed between the periods of 1982-1999 and 2000-2010. In the Northern Hemisphere of Africa, OGI trends were mostly insignificant, but prolonged GSL was evident over individual climate regions during the last 3 decades. OGI mainly showed late trends in the Southern Hemisphere of Africa while GSL was reversed from reduced GSL trends (1982-1999) to prolonged trends (2000-2010). In Australia, GSL exhibited considerable interannual variation, but the consistent trend lacked presence in most regions. Finally, the proportion of pixels with significant trends was less than I% in most of climate regions although it could be as large as 10%.

remote sensing↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
NASA NTRS2001

Mapping Arid Vegetation Species Distributions in the White Mountains, Eastern California, Using AVIRIS, Topography, and Geology

Our challenge is to model plant species distributions in complex montane environments using disparate sources of data, including topography, geology, and hyperspectral data. From an ecologist's point of view, species distributions are determined by local environment and disturbance history, while spectral data are 'ancillary.' However, a remote sensor's perspective says that spectral data provide picture of what vegetation is there, topographic and geologic data are ancillary. In order to bridge the gap, all available data should be used to get the best possible prediction of species distributions using complex multivariate techniques implemented on a GIS. Vegetation reflects local climatic and nutrient conditions, both of which can be modeled, allowing predictive mapping of vegetation distributions. Geologic substrate strongly affects chemical, thermal, and physical properties of soils, while climatic conditions are determined by local topography. As elevation increases, precipitation increases and temperature decreases. Aspect, slope, and surrounding topography determine potential insolation, so that south-facing slopes are warmer and north-facing slopes cooler at a given elevation. Topographic position (ridge, slope, canyon, or meadow) and slope angle affect sediment accumulation and soil depth. These factors combine as complex environmental gradients, and underlie many features of plant distributions. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data, digital elevation models, digitized geologic maps, and 378 ground control points were used to predictively map species distributions in the central and southern White Mountains, along the western boundary of the Basin and Range province. Minimum Noise Fraction (MNF) bands were calculated from the visible and near-infrared AVIRIS bands, and combined with digitized geologic maps and topographic variables using Canonical Correspondence Analysis (CCA). CCA allows for modeling species 'envelopes' in multidimensional environmental space, which can then be projected across entire landscapes.

VandeVen, C.↗
DOE OSTI2022

Vegetation type is an important predictor of the arctic summer land surface energy budget

Despite the importance of high-latitude surface energy budgets (SEBs) for land-climate interactions in the rapidly changing Arctic, uncertainties in their prediction persist. Here, we harmonize SEB observations across a network of vegetated and glaciated sites at circumpolar scale (1994–2021). Our variance-partitioning analysis identifies vegetation type as an important predictor for SEB-components during Arctic summer (June-August), compared to other SEB-drivers including climate, latitude and permafrost characteristics. Differences among vegetation types can be of similar magnitude as between vegetation and glacier surfaces and are especially high for summer sensible and latent heat fluxes. The timing of SEB-flux summer-regimes (when daily mean values exceed 0 Wm –2 ) relative to snow-free and -onset dates varies substantially depending on vegetation type, implying vegetation controls on snow-cover and SEB-flux seasonality. Our results indicate complex shifts in surface energy fluxes with land-cover transitions and a lengthening summer season, and highlight the potential for improving future Earth system models via a refined representation of Arctic vegetation types.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

Carbon Emissions from Deforestation in the Brazilian Amazon Region

A simulation model based on satellite observations of monthly vegetation greenness from the Moderate Resolution Imaging Spectroradiometer (MODIS) was used to estimate monthly carbon fluxes in terrestrial ecosystems of Brazilian Amazon and Cerrado regions over the period 2000-2002. The NASA-CASA (Carnegie Ames Stanford Approach) model estimates of annual forest production were used for the first time as the basis to generate a prediction for the standing pool of carbon in above-ground biomass (AGB; gC/sq m) for forested areas of the Brazilian Amazon region. Plot-level measurements of the residence time of carbon in wood in Amazon forest from Malhi et al. (2006) were interpolated by inverse distance weighting algorithms and used with CASA to generate a new regional map of AGB. Data from the Brazilian PRODES (Estimativa do Desflorestamento da Amazonia) project were used to map deforested areas. Results show that net primary production (NPP) sinks for carbon varied between 4.25 Pg C/yr (1 Pg=10(exp 15)g) and 4.34 Pg C for the region and were highest across the eastern and northern Amazon areas, whereas deforestation sources of CO2 flux from decomposition of residual woody debris were higher and less seasonal in the central Amazon than in the eastern and southern areas. Increased woody debris from past deforestation events was predicted to alter the net ecosystem carbon balance of the Amazon region to generate annual CO2 source fluxes at least two times higher than previously predicted by CASA modeling studies. Variations in climate, land cover, and forest burning were predicted to release carbon at rates of 0.5 to 1 Pg C/yr from the Brazilian Amazon. When direct deforestation emissions of CO2 from forest burning of between 0.2 and 0.6 Pg C/yr in the Legal Amazon are overlooked in regional budgets, the year-to-year variations in this net biome flux may appear to be large, whereas our model results implies net biome fluxes had actually been relatively consistent from year to year during the period 2000-2002. This is the first study to use MODIS data to model all carbon pools (wood, leaf, root) dynamically in simulations of Amazon forest deforestation from clearing and burning of all kinds.

Potter, C.↗
DOE OSTI2024

Improving North American Wildfire Prediction by Integrating a Machine-Learning Fire Model in a Land Surface Model

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simulated wildfire burned area over the CONUS during 2001-2020

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM). A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

Liu, Ye↗
NASA NTRS1981

Effect of sun and sensor geometry, canopy structure and density, and atmospheric condition on the spectral response of vegetation, with particular emphasis on across-track pointing

A computer modeling and simulation study carried out to assess the effects of various sun and sensor geometries and atmospheric conditions on the directional reflected radiance of several vegetated targets is described. Spectral responses at two wavelengths, 0.68 micron and 0.80 micron, are simulated at nine sensor zenith angles, five sensor azimuths, and nine solar zenith angles for six vegetation canopies under three atmspheric conditions. The six canopies comprise two different geometries of grass canopies at low, medium, and high leaf density. The results suggest that off-nadir viewing effects are more pronounced in the red than in the IR. However, the use of such transformations as the normalized difference index is found to reduce much of the variability seen in the bands. The magnitude of off-nadir viewing effects is found to be a function of canopy geometry.

Schnetzler, C. C.↗
NASA NTRS2011

Effective Tree Scattering and Opacity at L-Band

This paper investigates vegetation effects at L-band by using a first-order radiative transfer (RT) model and truck-based microwave measurements over natural conifer stands to assess the applicability of the tau-omega) model over trees. The tau-omega model is a zero-order RT solution that accounts for vegetation effects with effective vegetation parameters (vegetation opacity and single-scattering albedo), which represent the canopy as a whole. This approach inherently ignores multiple-scattering effects and, therefore, has a limited validity depending on the level of scattering within the canopy. The fact that the scattering from large forest components such as branches and trunks is significant at L-band requires that zero-order vegetation parameters be evaluated (compared) along with their theoretical definitions to provide a better understanding of these parameters in the retrieval algorithms as applied to trees. This paper compares the effective vegetation opacities, computed from multi-angular pine tree brightness temperature data, against the results of two independent approaches that provide theoretical and measured optical depths. These two techniques are based on forward scattering theory and radar corner reflector measurements, respectively. The results indicate that the effective vegetation opacity values are smaller than but of similar magnitude to both radar and theoretical estimates. The effective opacity of the zero-order model is thus set equal to the theoretical opacity and an explicit expression for the effective albedo is then obtained from the zero- and first- order RT model comparison. The resultant albedo is found to have a similar magnitude as the effective albedo value obtained from brightness temperature measurements. However, it is less than half of that estimated using the theoretical calculations (0.5 - 0.6 for tree canopies at L-band). This lower observed albedo balances the scattering darkening effect of the large theoretical albedo with a first-order multiple-scattering contribution. The retrieved effective albedo is different from theoretical definitions and not the albedo of single forest elements anymore, but it becomes a global parameter, which depends on all the processes taking place within the canopy, including multiple-scattering.

Kurum, Mehmet↗
DOE OSTI2025

Artificial neural networks estimate evapotranspiration for Miscanthus × giganteus as effectively as empirical model but with fewer inputs

Estimating actual evapotranspiration (ET) is particularly crucial for addressing how vegetation affects the water balance of ecosystems. ET estimation can be complex with empirical models due to their many parameters and reliance on aridity. In contrast, artificial neural networks (ANNs) could potentially estimate ET with fewer and more common meteorological parameters. In this study, we trained two ANNs, one using a feed-forward approach (FFN) and the other a nonlinear auto-regressive network (NARX), to predict ET and compared them to the commonly used empirical model Granger and Gray (GG). We trained our models on a nine-year eddy covariance (EC) dataset for Miscanthu s × giganteus ( M . × giganteus ) from Illinois (UIEF), then tested them using out-of-sample data from both UIEF and a different location in Iowa (SABR) to compare the accuracy of FFN, NARX, and GG models in estimating daily ET. A combination of air temperature (T a ) and solar radiation (R s ) was chosen as inputs due to the highest R 2 for FFN (R 2 = 0.79, 0.81, and 0.79 for training, testing, and validation, respectively) and only T a for NARX (R 2 = 0.70 for out-of-sample validation). The predictive power of the FFN model was superior to the NARX and GG models at the UIEF site (R 2 = 0.84, 0.70, and 0.83 for out-of-sample validation, respectively). Our analysis showed that ANN approaches are as accurate as empirical approaches for estimating ET but use fewer inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Median Isn’t the Message: soil nutrient hot spots have a disproportionate influence on biogeochemical structure across years, seasons, and depths

Abstract Soil nutrient distribution is heterogeneous in space and time, potentially altering nutrient acquisition by trees and microorganisms. Ecologists have distinguished “hot spots” (HSs) as areas with enhanced and sustained rates of nutrient fluxes relative to the surrounding soil matrix. We evaluated the spatial and temporal patterns in nutrient flux HSs in two mixed-conifer forest soils by repeatedly sampling the soil solution at the same spatial locations (horizontally and vertically) over multiple seasons and years using ion exchange resins incubated in situ. The climate of these forests is Mediterranean, with intense fall rains occurring following summers with little precipitation, and highly variable winter snowfall. Hot spots formed most often for NO 3 − and Na + . Although nutrient HSs often occurred in the same spatial location multiple times, HSs persisted more often for PO 4 3− NH 4 + , and NO 3 − , and were more transient for Ca 2+ , Mg 2+ , and Na + . Sampling year (annual precipitation ranged from 558 to 1223 mm) impacted the occurrence of HSs for most nutrients, but season was only significant for PO 4 3− , NH 4 + , NO 3 − , and Na + , with HSs forming more often after fall rains than after spring snowmelt. The frequency of HSs significantly decreased with soil depth for all nutrients, forming most commonly immediately below the surficial organic horizon. Although HSs accounted for less than 17% of the sampling volume, they were responsible for 56–88% of PO 4 3− , NH 4 + , and NO 3 − resin fluxes. Our results suggest that macronutrient HSs have a disproportional contribution to soil biogeochemical structure, with implications for vegetation nutrient acquisition strategies and biogeochemical models. Graphical abstract

54 ENVIRONMENTAL SCIENCES↗