USGS Science⌕ Search

USGS · 70217587

Machine learning predictions of pH in the Glacial Aquifer System, Northern USA

Abstract

A boosted regression tree model was developed to predict pH conditions in three dimensions throughout the glacial aquifer system of the contiguous United States using pH measurements in samples from 18,386 wells and predictor variables that represent aspects of the hydrogeologic setting. Model results indicate that the carbonate content of soils and aquifer materials strongly controls pH and, when coupled with long flowpaths, results in the most alkaline conditions. Conversely, in areas where glacial sediments are thin and carbonate‐poor, pH conditions remain acidic. At depths typical of drinking‐water supplies, predicted pH >7.5—which is associated with arsenic mobilization—occurs more frequently than predicted pH <6—which is associated with water corrosivity and the mobilization of other trace elements. A novel aspect of this model was the inclusion of numerically based estimates of groundwater flow characteristics (age and flowpath length) as predictor variables. The sensitivity of pH predictions to these variables was consistent with hydrologic understanding of groundwater flow systems and the geochemical evolution of groundwater quality. The model was not developed to provide precise estimates of pH at any given location. Rather, it can be used to more generally identify areas where contaminants may be mobilized into groundwater and where corrosivity issues may be of concern to prioritize areas for future groundwater monitoring.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Paul E. Stackelberg, Kenneth Belitz, Craig J. Brown, Melinda L. Erickson, Sarah M. Elliott, Leon J. Kauffman, Katherine Marie Ransom, James E. Reddy. 2020-12-31. Machine learning predictions of pH in the Glacial Aquifer System, Northern USA. https://doi.org/10.1111/gwat.13063

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Computing flow-field distortion coefficients from well-construction and formation properties

Direct measurements of groundwater velocity made with borehole flowmeters in screened wells must be compensated for the effects of flow-field distortion (also known as borehole acceleration). A theoretical equation developed by Drost et al. (1968) and simple inputs describing hydraulic properties of well construction and geologic formation were programmed into an Excel workbook to facilitate computation by groundwater-flowmeter users. Tables describing the physical and hydraulic properties for well constructions and gravel pack media are provided with an example to facilitate use of the workbook. Groundwater flowlines converge or diverge as they pass from a geologic formation, through a gravel pack and well screen. The extent of flowline convergence or divergence and the value of the flow-field distortion coefficient is related to the relative changes in hydraulic conductivity of the well screen, gravel pack, and geologic formation. Convergence or divergence is accompanied by acceleration or deceleration of groundwater. Direct measurements of groundwater velocity at the center of the monitoring well can be adjusted to provide a more accurate estimate of velocity in the formation by applying a correction for flow-field distortion. Variables required to compute the flow-field distortion coefficient include the hydraulic conductivity of the gravel pack, well screen, and the geologic formation surrounding the well screen; the borehole radius, and the inside radius and outside radius of the well screen.

Groundwater↗

Predicting groundwater hydrochemical facies in three dimensions with random forest classification, USA

A random forest classification (RFC) model was developed to predict hydrochemical facies (HCFs) of groundwater in three dimensions across the conterminous United States (CONUS). Major-ion data from 152,673 sites were used to categorize groundwater into one of six HCFs (CaMg-HCO 3 , NaK-HCO 3 , CaMg-SO 4 , NaK-SO 4 , Cl, or Mixed). These six HCFs were used as targets for RFC modeling. Model features that represent relevant geochemical processes and/or physical conditions were derived from previously published data. Additional model features were specifically engineered to support this analysis: elevation of the bottom of a well relative to the base of drinking water (ERDW) and six flags that relate geologic units to HCFs. The most important model feature was ERDW. The model was used to map HCFs at a 1-km 2 resolution across CONUS and to depths of 400 m below the base of drinking water (which varies from 22 m to 2 km). Model predictions are consistent with expectations. CaMg-HCO 3 is predicted to occur near the water table in more humid settings, and areas underlain by carbonate or crystalline rocks. At depths below the base of drinking-water supplies, the model predicts a rapid transition from HCO 3 HCFs to Cl. Model predictions are accurate based on point data, and data averaged across hydrogeologic regions and with depth. Model predictions of HCFs could be used for multiple purposes, including the mapping of salinity and other groundwater characteristics.

Groundwater↗