USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 577 records · Page 32Linked to original sources
USGS2012

Monitoring groundwater-surface water interaction using time-series and time-frequency analysis of transient three-dimensional electrical resistivity changes

Time-lapse resistivity imaging is increasingly used to monitor hydrologic processes. Compared to conventional hydrologic measurements, surface time-lapse resistivity provides superior spatial coverage in two or three dimensions, potentially high-resolution information in time, and information in the absence of wells. However, interpretation of time-lapse electrical tomograms is complicated by the ever-increasing size and complexity of long-term, three-dimensional (3-D) time series conductivity data sets. Here we use 3-D surface time-lapse electrical imaging to monitor subsurface electrical conductivity variations associated with stage-driven groundwater-surface water interactions along a stretch of the Columbia River adjacent to the Hanford 300 near Richland, Washington, USA. We reduce the resulting 3-D conductivity time series using both time-series and time-frequency analyses to isolate a paleochannel causing enhanced groundwater-surface water interactions. Correlation analysis on the time-lapse imaging results concisely represents enhanced groundwater-surface water interactions within the paleochannel, and provides information concerning groundwater flow velocities. Time-frequency analysis using the Stockwell (S) transform provides additional information by identifying the stage periodicities driving groundwater-surface water interactions due to upstream dam operations, and identifying segments in time-frequency space when these interactions are most active. These results provide new insight into the distribution and timing of river water intrusion into the Hanford 300 Area, which has a governing influence on the behavior of a uranium plume left over from historical nuclear fuel processing operations.

Washington↗
USGS2012

Practical estimates of field-saturated hydraulic conductivity of bedrock outcrops using a modified bottomless bucket method

The bottomless bucket (BB) approach (Nimmo et al., 2009a) is a cost-effective method for rapidly characterizing field-saturated hydraulic conductivity Kfs of soils and alluvial deposits. This practical approach is of particular value for quantifying infiltration rates in remote areas with limited accessibility. A similar approach for bedrock outcrops is also of great value for improving quantitative understanding of infiltration and recharge in rugged terrain. We develop a simple modification to the BB method for application to bedrock outcrops, which uses a non-toxic, quick-drying silicone gel to seal the BB to the bedrock. These modifications to the field method require only minor changes to the analytical solution for calculating Kfs on soils. We investigate the reproducibility of the method with laboratory experiments on a previously studied calcarenite rock and conduct a sensitivity analysis to quantify uncertainty in our predictions. We apply the BB method on both bedrock and soil for sites on Pahute Mesa, which is located in a remote area of the Nevada National Security Site. The bedrock BB tests may require monitoring over several hours to days, depending on infiltration rates, which necessitates a cover to prevent evaporative losses. Our field and laboratory results compare well to Kfs values inferred from independent reports, which suggests the modified BB method can provide useful estimates and facilitate simple hypothesis testing. The ease with which the bedrock BB method can be deployed should facilitate more rapid in-situ data collection than is possible with alternative methods for quantitative characterization of infiltration into bedrock.

Water Resources Research↗
USGS2012

The influence of reservoirs, climate, land use and hydrologic conditions on loads and chemical quality of dissolved organic carbon in the Colorado River

Longitudinal patterns in dissolved organic carbon (DOC) loads and chemical quality were identified in the Colorado River from the headwaters in the Rocky Mountains to the United States-Mexico border from 1994 to 2011. Watershed- and reach-scale climate, land use, river discharge and hydrologic modification conditions that contribute to patterns in DOC were also identified. Principal components analysis (PCA) identified site-specific precipitation and reach-scale discharge as being correlated with sites in the upper basin, where there were increases in DOC load from the upstream to downstream direction. In the lower basin, where DOC load decreased from upstream to downstream, sites were correlated with site-specific temperature and reach-scale population, urban land use and hydrologic modification. In the reaches containing Lakes Powell and Mead, the two largest reservoirs in the United States, DOC quantity decreased, terrestrially derived aromatic DOC was degraded and/or autochthonous less aromatic DOC was produced. Taken together, these results suggest that longitudinal patterns in the relatively unregulated upper basin are influenced by watershed inputs of water and DOC, whereas DOC patterns in the lower basin are reflective of a balance between watershed contribution of water and DOC to the river and loss of water and DOC due to hydrologic modification and/or biogeochemical processes. These findings suggest that alteration of constituent fluxes in rivers that are highly regulated may overshadow watershed processes that would control fluxes in comparable unregulated rivers. Further, these results provide a foundation for detailed assessments of factors controlling the transport and chemical quality of DOC in the Colorado River.

Arizona, California, Colorado, Nevada, New Mexico,↗
USGS2012

Role of surface-water and groundwater interactions on projected summertime streamflow in snow dominated regions : An integrated modeling approach

Previous studies indicate predominantly increasing trends in precipitation across the Western United States, while at the same time, historical streamflow records indicate decreasing summertime streamflow and 25th percentile annual flows. These opposing trends could be viewed as paradoxical, given that several studies suggest that increased annual precipitation will equate to increased annual groundwater recharge, and therefore increased summertime flow. To gain insight on mechanisms behind these potential changes, we rely on a calibrated, integrated surface and groundwater model to simulate climate impacts on surface water/groundwater interactions using 12 general circulation model projections of temperature and precipitation from 2010 to 2100, and evaluate the interplay between snowmelt timing and other hydrologic variables, including streamflow, groundwater recharge, storage, groundwater discharge, and evapotranspiration. Hydrologic simulations show that the timing of peak groundwater discharge to the stream is inversely correlated to snowmelt runoff and groundwater recharge due to the bank storage effect and reversal of hydraulic gradients between the stream and underlying groundwater. That is, groundwater flow to streams peaks following the decrease in stream depth caused by snowmelt recession, and the shift in snowmelt causes a corresponding shift in groundwater discharge to streams. Our results show that groundwater discharge to streams is depleted during the summer due to earlier drainage of shallow aquifers adjacent to streams even if projected annual precipitation and groundwater recharge increases. These projected changes in surface water/groundwater interactions result in more than a 30% decrease in the projected ensemble summertime streamflow. Our findings clarify causality of observed decreasing summertime flow, highlight important aspects of potential climate change impacts on groundwater resources, and underscore the need for integrated hydrologic models in climate change studies.

California;Nevada↗
USGS2012

Wildfire impacts on soil-water retention in the Colorado Front Range, United States

This work examined the plot-scale differences in soil-water retention caused by wildfire in the area of the 2010 Fourmile Canyon Fire in the Colorado Front Range, United States. We measured soil-water retention curves on intact cores and repacked samples, soil particle-size distributions, and organic matter content. Estimates were also made of plant-available water based on the soil-water retention curves. Parameters for use in soil-hydraulic property models were estimated; these parameters can be used in unsaturated flow modeling for comparing burned and unburned watersheds. The primary driver for measured differences in soil-water retention in burned and unburned soils was organic matter content and not soil-particle size distribution. The tendency for unburned south-facing soils to have greater organic matter content than unburned north-facing soils in this field area may explain why unburned south-facing soils had greater soil-water retention than unburned north-facing soils. Our results suggest that high-severity wildfire can “homogenize” soil-water retention across the landscape by erasing soil-water retention differences resulting from organic matter content, which for this site may be affected by slope aspect. This homogenization could have important implications for ecohydrology and plant succession/recovery in burned areas, which could be a factor in dictating the window of vulnerability of the landscape to flash floods and erosion that are a common consequence of wildfire.

Colorado↗
USGS2012

Direct geoelectrical evidence of mass transfer at the laboratory scale

Previous field-scale experimental data and numerical modeling suggest that the dual-domain mass transfer (DDMT) of electrolytic tracers has an observable geoelectrical signature. Here we present controlled laboratory experiments confirming the electrical signature of DDMT and demonstrate the use of time-lapse electrical measurements in conjunction with concentration measurements to estimate the parameters controlling DDMT, i.e., the mobile and immobile porosity and rate at which solute exchanges between mobile and immobile domains. We conducted column tracer tests on unconsolidated quartz sand and a material with a high secondary porosity: the zeolite clinoptilolite. During NaCl tracer tests we collected nearly colocated bulk direct-current electrical conductivity (σ b ) and fluid conductivity (σ f ) measurements. Our results for the zeolite show (1) extensive tailing and (2) a hysteretic relation between σ f and σ b , thus providing evidence of mass transfer not observed within the quartz sand. To identify best-fit parameters and evaluate parameter sensitivity, we performed over 2700 simulations of σ f , varying the immobile and mobile domain and mass transfer rate. We emphasized the fit to late-time tailing by minimizing the Box-Cox power transformed root-mean square error between the observed and simulated σ f . Low-field proton nuclear magnetic resonance (NMR) measurements provide an independent quantification of the volumes of the mobile and immobile domains. The best-fit parameters based on σ f match the NMR measurements of the immobile and mobile domain porosities and provide the first direct electrical evidence for DDMT. Our results underscore the potential of using electrical measurements for DDMT parameter inference.

Water Resources Research↗
USGS2018

Regionalization of groundwater residence time using metamodeling

Groundwater residence-time distributions (RTDs) are critical for assessing susceptibility of water resources to degradation. A novel combination of numerical modeling and statistical methods allows estimation of regional RTDs with unprecedented speed. In this method, particle RTDs are generated in 30 type locales in the northeastern glaciated U.S using automated generalized finite-difference groundwater flow and advective transport models. Targets for statistical learning were created from particle RTDs by fitting Weibull, gamma, and inverse Gaussian distributions. Whole-basin flux-weighted RTDs were well fit by one-component Weibull distributions. Flux-weighted RTDs at stressed receptors such as wells often produced more complicated RTDs that required a two-component mixture to fit. A Multitask Lasso regression was trained on the parametric RTDs using hydrogeographic features of the modeled areas as explanatory features. In this way, RTDs are regionalized using mappable physical features such as recharge and aquifer volume. The shape, location, and scale parameters of the parametric RTDs are strongly related to the mean exponential age. The shape parameter of the distribution, which controls deviation from exponential, is additionally a function of aquifer heterogeneity and hydrologic features. Regionalized RTDs provide useful metrics with respect to groundwater lag times and solute loading to streams. The lag time between input and output contained in the RTD is critical to understanding the relation between the land surface and human and ecological receptors.

Water Resources Research↗
USGS2018

Direct channel precipitation and storm type influence short-term fallout radionuclide assessment of sediment source

Fallout radionuclides (FRNs) and their ratios, such as Beryllium‐7 ( 7 Be) and excess Lead‐210 ( 210 Pb xs ), have been used to determine suspended sediment source and age in catchments. These models are based on numerous assumptions, for example, that channel deposition of FRNs from precipitation is negligible in comparison to their delivery to the channel from land surface erosion during individual storm events. We test this assumption using a mass balance approach during eight storms from summer 2011 to fall 2012 in a mid‐Atlantic United States piedmont region watershed with mixed land use. Event peak discharge and storm type corresponded to the importance of direct channel FRN deposition from precipitation. During relatively low discharge summer thunderstorms, with minimal overland flow, less than 1% of 7 Be and 210 Pb xs flux deposited on the watershed exits the watershed associated with suspended sediment. The majority but not all deposited on the stream channel exits the watershed associated with suspended sediment (60% of 7 Be and 80% of 210 Pb xs ). Here precipitation and throughfall onto the wetted channel area can be responsible for any FRN newly associated with suspended sediment, as opposed to landscape surface erosion. Furthermore, FRNs can be stored with sediments in the channel between events. Events with higher discharges, including hurricanes, show the opposite pattern—FRN flux associated with suspended sediment exported from the reach is greater than channel FRN wet deposition, suggesting net erosion from the watershed landscape and/or stored material during these types of storms.

Water Resources Research↗
USGS2018

What makes a first‐magnitude spring?: Global sensitivity analysis of a speleogenesis model to gain insight into karst network and spring genesis

Often, karstic conduit network geometry is unknown. This lack of knowledge represents a significant limitation when modeling flow and solute transport in karst systems. In this study, we apply Morris Method Global Sensitivity Analysis to a speleogenesis model to identify model input parameters, and combinations thereof, that most significantly influence evolution of karst conduit networks, development of first‐magnitude springs, and resulting flow and solute transport pulse responses. Based on an idealized model of the Silver Springshed in Central Florida USA, results showed that porous matrix hydraulic conductivity and parameters that govern connectivity of vertical and horizontal preferential flow paths (proto‐conduits) are the most influential parameters. In particular, a lower porous matrix conductivity is more likely to produce a first‐order magnitude spring. For the boundary conditions assumed in this application, conduits tend to develop in low topographic regions that drained nearby high regions. Morris ensemble realizations that generated first‐magnitude springs exhibit similar flow and solute transport pulse responses at the spring vent, despite differences in network configuration. However, distributed head fields are highly spatially variable, implying substantial spatial variability among solute flow paths and travel times from the land surface to the spring across realizations.

Water Resources Research↗
USGS2018

Mechanisms of earthquake‐induced chemical and fluid transport to carbonate groundwater springs after earthquakes

Mechanisms by which hydrochemical changes occur after earthquakes are not well documented. We use the 2016–2017 central Italy seismic sequence, which caused notable hydrochemical transient variations in groundwater springs to address this topic, with special reference to effects on fractured carbonate aquifers. Hydrochemistry measured before and after the earthquakes at four springs at varying distances from the epicenters all showed immediate postmainshock peaks in trace element concentrations but little change in major elements. Most parameters returned to preearthquake values before the last events of the seismic sequence. The source of solutes, particularly trace elements, is longer residence time pore water stored in slow‐moving fractures or abandoned karstic flow paths. These fluids were expelled into the main flow paths after an increase in pore pressure, hydraulic conductivity, and shaking from coseismic aquifer stress. The weak response to the later earthquakes is explained by progressive depletion of high solute fluids as earlier shocks flushed out the stored fluids in the fractures. Spring δ 13 C DIC values closest to a deep magma source to the west became enriched relative to preearthquake values following the 24 August event. This enrichment indicates input from deeply sourced dissolved CO 2 gas after dilation of specific fault conduits. Differences in carbon isotopic responses between springs are attributed to proximity to the deep gaseous CO 2 source. Most of the transient chemical changes seen in the three fractured carbonate aquifers are attributed to local shaking and emptying of isolated pores and fractures and are not from rapid upward movement of deep fluids.

Water Resources Research↗
USGS2018

Rethinking the use of seabed sediment temperature profiles to trace submarine groundwater flow

Submarine groundwater fluxes across the seafloor facilitate important hydrological and biogeochemical exchanges between oceans and seabed sediment, yet few studies have investigated spatially distributed groundwater fluxes in deep‐ocean environments such as continental slopes. Heat has been previously applied as a submarine groundwater tracer using an analytical solution to a heat flow equation assuming steady state conditions and homogeneous thermal conductivity. These assumptions are often violated in shallow seabeds due to ocean bottom temperature changes or sediment property variations. Here heat tracing analysis techniques recently developed for terrestrial settings are applied in concert to examine the influences of groundwater flow, ocean temperature changes, and seabed thermal conductivity variations on deep‐ocean sediment temperature profiles. Temperature observations from the sediment and bottom ocean water on the Scotian Slope off eastern Canada are used to demonstrate how simple thermal methods for tracing groundwater can be employed if more comprehensive techniques indicate that the simplifying assumptions are valid. The spatial distribution of the inferred groundwater fluxes on the slope suggests a downward groundwater flow system with recharge occurring over the upper‐middle slope and discharge on the lower slope. We speculate that the downward groundwater flow inferred on the Scotian Slope is due to density‐driven processes arising from underlying salt domes, in contrast with upward slope systems driven by geothermal convection. Improvements in the design of future submarine hydrogeological studies are proposed for thermal data collection and groundwater flow analysis, including new equations that quantify the minimum detectable flux magnitude for a given sensor accuracy and profile length.

Water Resources Research↗
USGS2018

Metamodeling for groundwater age forecasting in the Lake Michigan Basin

Groundwater age is an important indicator of groundwater susceptibility to anthropogenic contamination and a key input to statistical models for forecasting water quality. Numerical models can provide estimates of groundwater age, enabling interpretation of measured age tracers. However, to extend to national‐scale groundwater systems where numerical models are not routinely available, a more efficient metamodeling approach can provide a less precise but widely applicable estimate of groundwater age, trained to make forecasts based on predictor variables that can be measured independent of numerical models. We trained gradient‐boosted regression tree statistical metamodels to MODFLOW/MODPATH‐derived groundwater age estimates in five inset models in the Lake Michigan Basin, USA. Using high‐throughput computing, we explored an exhaustive range of tuning parameters and tested metamodels through cross validation, a 20% holdout, and a round robin approach among the five inset models withholding each inset model from training and testing on the held‐out inset model. Forecast skill—measured by Nash Sutcliffe efficiency—was high for age‐related responses in the 20% hold‐out case (ranging from 0.73 to 0.84). The round robin analysis provided the opportunity to explore extending to unmodeled areas and a greater range of skill indicated the need to evaluate when it is appropriate to apply a metamodel from one region to another. We further explored the ramifications of metamodel simplification achieved through removing predictor variables based on their estimated importance. We found that similar metamodel performance was achievable with a fraction of the candidate set of predictor variables with well construction variables being most important.

Water Resources Research↗
USGS2018

Exploring drivers of regional water-quality change using differential spatially referenced regression – A pilot study in the Chesapeake Bay watershed

An understanding of riverine water-quality dynamics in regional mixed-land use watersheds is the foundation for advances in landscape biogeochemistry and informed land management. A differential implementation of the statistical/process-based model SPAtially Referenced Regressions on Watershed attributes (SPARROW; Smith et al., https://doi.org/10.1029/97wr02171 ) is proposed to empirically relate a regional pattern of changes in flow-normalized constituent flux, over a multiyear period, to contemporaneous changes in spatially referenced explanatory variables. In a pilot application, the differential model, called Spatiotemporal Watershed Accumulation of Net effects (SWAN), is used to explore factors influencing changes in flow-normalized flux of total nitrogen over the period 1990–2010 at 43 sites in the nontidal Chesapeake Bay watershed. A seven-parameter model explains 80% of the transformed variability in independently estimated flux changes, indicating that storage effects having characteristic time scales greater than 20 years had a small influence, relative to changes in inputs, on regional water-quality response. Results suggest that 1990–2010 changes in total-nitrogen flux are largely the outcome of increased nonpoint-source pollution associated with urban and suburban development, modulated to the point of negation by terrestrial losses stemming from widespread increases in air temperature and precipitation. The loss mechanism is qualitatively consistent with denitrification; however, increases in aboveground biomass, agricultural nitrogen exports, or hydrologic flushing are also plausible contributors. Although qualified by a small sample size and constraints on explanatory data availability, the pilot suggests that SWAN is a promising approach for broadening scientific understanding of factors driving regional water-quality change and for supporting evidence-based land-management decisions.

Chesapeake Bay watershed↗
USGS2018

A physically based daily simulation of the glacier-dominated hydrology of the Copper River basin, Alaska

The large, highly glacierized Copper River basin is an important water resource for the south-central region of Alaska. Thus, information is needed on the reaction of its hydrologic timing and streamflow volumes to historical changes in climate, in order to assess the possible impact of future changes. However, the basin is remote, and therefore, it has proved difficult to collect field data in a frequent temporal and spatial manner. An extension of the distributed-parameter, physical-process code Precipitation Runoff Modeling System, PRMSglacier, has been specifically developed to simulate daily hydrology without requiring extensive input data. In this study, PRMSglacier was used to characterize the hydrology of the Copper River basin from 1959 to 2015. The basin was split into subbasins for specific regional climatic calibrations and finer resolution characterization. The model was calibrated and performed well against data of glacier mass balance, glacier area change, snow cover, gaged streamflow, evapotranspiration, and solar radiation. Ice melt contributed 26% of the total basin streamflow, with differences temporally from climate oscillations. Furthermore, differences were seen geographically in subbasins depending on the state of the glaciers in each subbasin. Decreasing trends in ice melt volume were mostly seen on smaller steeper glaciers responding to a critical level of glacier recession, while increasing trends in ice melt volume were mostly seen on larger valley glaciers responding to increasing temperature. The areas with substantially decreasing ice melt had decreasing streamflow, possibly indicating health concerns for the ecosystems therein.

Alaska↗
USGS2018

Stream‐centric methods for determining groundwater contributions in karst mountain watersheds

Climate change influences on mountain hydrology are uncertain, but likely to be mediated through changes in subsurface hydrologic residence times and flowpaths. The heterogeneity of karst aquifers add complexity in assessing the resiliency of these water sources to perturbation, suggesting a clear need to quantify contributions from and losses to these aquifers. Here we develop a stream centric method that combines mass and flow balances to quantify net and gross gains and losses at different spatial scales. We then extend these methods to differentiate between karst conduit and matrix contributions from the aquifer. In the Logan River watershed in Northern Utah we found significant amounts of the river water repeatedly gained and then lost through a 35 km study reach. Further, the direction and amount of water exchanged varied over space, time, and discharge. Streamflow was dominated by discharge of karst conduit groundwater after runoff with increasing, yet still small, fractions of matrix water later in the summer. These findings were combined with geologic information, prior subsurface dye tracing, and chemical sampling to provide additional lines of evidence that repeated groundwater exchanges are likely occurring and river flows are highly dependent on karst aquifer recharge and discharge. Given the large population dependent on karst aquifers throughout the world, there is a continued need to develop simple methods, like those presented here, for determining the resiliency of karst groundwater resources.

Water Resources Research↗
USGS2018

A comparison of methods for streamflow uncertainty estimation

Streamflow time series are commonly derived from stage‐discharge rating curves, but the uncertainty of the rating curve and resulting streamflow series are poorly understood. While different methods to quantify uncertainty in the stage‐discharge relationship exist, there is limited understanding of how uncertainty estimates differ between methods due to different assumptions and methodological choices. We compared uncertainty estimates and stage‐discharge rating curves from seven methods at three river locations of varying hydraulic complexity. Comparison of the estimated uncertainties revealed a wide range of estimates, particularly for high and low flows. At the simplest site on the Isère River (France), full width 95% uncertainties for the different methods ranged from 3 to 17% for median flows. In contrast, uncertainties were much higher and ranged from 41 to 200% for high flows in an extrapolated section of the rating curve at the Mahurangi River (New Zealand) and 28 to 101% for low flows at the Taf River (United Kingdom), where the hydraulic control is unstable at low flows. Differences between methods result from differences in the sources of uncertainty considered, differences in the handling of the time‐varying nature of rating curves, differences in the extent of hydraulic knowledge assumed, and differences in assumptions when extrapolating rating curves above or below the observed gaugings. Ultimately, the selection of an uncertainty method requires a match between user requirements and the assumptions made by the uncertainty method. Given the significant differences in uncertainty estimates between methods, we suggest that a clear statement of uncertainty assumptions be presented alongside streamflow uncertainty estimates.

Water Resources Research↗
USGS2018

Tidal response of groundwater in a leaky aquifer—Application to Oklahoma

Quantitative interpretation of the tidal response of water levels measured in wells has long been made either with a model for perfectly confined aquifers or with a model for purely unconfined aquifers. However, many aquifers may be neither totally confined nor purely unconfined at the frequencies of tidal loading but behave somewhere between the two end‐members. Here we present a more general model for the tidal response of groundwater in aquifers with both horizontal flow and vertical leakage. The model has three independent parameters: the transmissivity ( T ) and storativity ( S ) of the aquifer and the specific leakage ( K′/b′ ) of the leaking aquitard, where K′ and b′ are the hydraulic conductivity and the thickness of the aquitard, respectively. If T and S are known independently, this model may be used to estimate aquitard leakage from the phase shift and amplitude ratio of water level in wells obtained from tidal analysis. We apply the model to interpret the tidal response of water level in a US Geological Survey (USGS) deep monitoring well installed in the Arbuckle aquifer in Oklahoma, into which massive amount of wastewater coproduced from hydrocarbon exploration has been injected. The analysis shows that the Arbuckle aquifer is leaking significantly at this site. We suggest that the present method may be effective and economical for monitoring leakage in groundwater systems, which bears on the safety of water resources, the security of underground waste repositories, and the outflow of wastewater during deep injection and hydrocarbon extraction.

Oklahoma↗
USGS2019

POLARIS properties: 30-meter probabilistic maps of soil properties over the contiguous United States

Soils play a critical role in the cycling of water, energy, and carbon in the Earth system. Until recently, due primarily to a lack of soil property maps of a sufficiently high‐quality and spatial detail, a minor emphasis has been placed on providing high‐resolution measured soil parameter estimates for land surface models and hydrologic models. This study introduces Probabilistic Remapping of SSURGO (POLARIS) soil properties—a database of 30‐m probabilistic soil property maps over the contiguous United States (CONUS). The mapped variables over CONUS include soil texture, organic matter, pH, saturated hydraulic conductivity, Brooks‐Corey and Van Genuchten water retention curve parameters, bulk density, and saturated water content. POLARIS soil properties was assembled by (1) depth harmonizing and aggregating the pedons in the National Cooperative Soil Survey Soil Characterization Database and the components in Soil Survey Geographic Database into a database of 21,481 different soil series, each soil series having its own vertical profiles of different soil properties, (2) pruning the original POLARIS soil series maps using conventional soil maps to improve soil series prediction accuracy, and (3) merging the assembled soil series databases with the pruned POLARIS soil series maps to construct the soil property maps over CONUS. POLARIS soil properties includes 100‐bin histograms for each layer and variable per grid cell and a series of summary statistics at 30‐, 300‐, and 3,000‐m spatial resolution. Evaluation of POLARIS soil properties using in situ measurements shows an average R2 of 0.41, normalized root‐mean‐square error of 12%, and a normalized mean absolute error of 8.8%.

Water Resources Research↗