Geology topics

Lixia Liao

Publications and source records attributed to Lixia Liao.

5 recordsLinked to original sources

USGS2021

Machine learning predictions of mean ages of shallow well samples in the Great Lakes Basin, USA

The travel time or “age” of groundwater affects catchment responses to hydrologic changes , geochemical reactions, and time lags between management actions and responses at down-gradient streams and wells. Use of atmospheric tracers has facilitated the characterization of groundwater ages, but most wells lack such measurements. This study applied machine learning to predict ages in wells across a large region around the Great Lakes Basin using well, chemistry, and landscape characteristics. For a dataset of age tracers in 961 samples, the travel time from the land surface to the sample location was estimated for each sample using parametric functions. The mean travel times were then modeled using a gradient boosting machine (GBM) algorithm with cross validation tuning of model metaparameters. The GBM approach was able to closely match estimated ages for the training data (RMSE = 0.26 natural-log scale years) and provided a reasonable match to testing data (RMSE = 0.84). Of the variables tested, well characteristics (e.g. depth), land use, hydrologic indicators (e.g. topographic wetness index), and water chemistry (e.g. nitrate, fluoride, and pH), substantially affected the predictions of age. GBM prediction was applied to 14,335 groundwater samples with median sample depth of 5.4 m, indicating for the Great Lakes Basin a broad distribution of ages among wells with a median of 32.9 years. Lag times of decades are likely for these wells to respond to changing solute fluxes near land surface. While depth variables most strongly affected predicted mean ages, chemical constituents exhibited smooth trends with age, consistent with prevailing conceptual models of evolving sources and geochemistry flowpaths. The results provide proof of concept for use of readily available variables of well, landscape, and chemical characteristics to improve groundwater age estimates across large regions.

Great Lakes basin↗

USGS2018

Regional variability of nitrate fluxes in the unsaturated zone and groundwater, Wisconsin, USA

Process-based modeling of regional NO3− fluxes to groundwater is critical for understanding and managing water quality, but the complexity of NO3− reactive transport processes make implementation a challenge. This study introduces a regional vertical flux method (VFM) for efficient estimation of reactive transport of NO3− in the vadose zone and groundwater. The regional VFM was applied to 443 well samples in central-eastern Wisconsin. Chemical measurements included O2, NO3−, N2 from denitrification, and atmospheric tracers of groundwater age including carbon-14, chlorofluorocarbons, tritium, and tritiogenic helium. VFM results were consistent with observed chemistry, and calibrated parameters were in-line with estimates from previous studies. Results indicated that (1) unsaturated zone travel times were a substantial portion of the transit time to wells and streams (2) since 1945 fractions of applied N leached to groundwater have increased for manure-N, possibly due to increased injection of liquid manure, and decreased for fertilizer-N, and (3) under current practices and conditions, approximately 60% of the shallow aquifer will eventually be affected by downward migration of NO3−, with denitrification protecting the remaining 40%. Recharge variability strongly affected the unsaturated zone lag times and the eventual depth of the NO3− front. Principal components regression demonstrated that VFM parameters and predictions were significantly correlated with hydrogeochemical landscape features. The diverse and sometimes conflicting aspects of N management (e.g. limiting N volatilization versus limiting N losses to groundwater) warrant continued development of large-scale holistic strategies to manage water quality and quantity.

Wisconsin↗

USGS2018

Metamodeling and mapping of nitrate flux in the unsaturated zone and groundwater, Wisconsin, USA

Nitrate contamination of groundwater in agricultural areas poses a major challenge to the sustainability of water resources. Aquifer vulnerability models are useful tools that can help resource managers identify areas of concern, but quantifying nitrogen (N) inputs in such models is challenging, especially at large spatial scales. We sought to improve regional nitrate (NO 3 − ) input functions by characterizing unsaturated zone NO 3 − transport to groundwater through use of surrogate, machine-learning metamodels of a process-based N flux model. The metamodels used boosted regression trees (BRTs) to relate mappable landscape variables to parameters and outputs of a previous “vertical flux method” (VFM) applied at sampled wells in the Fox, Wolf, and Peshtigo (FWP) river basins in northeastern Wisconsin. In this context, the metamodels upscaled the VFM results throughout the region, and the VFM parameters and outputs are the metamodel response variables. The study area encompassed the domain of a detailed numerical model that provided additional predictor variables, including groundwater recharge, to the metamodels. We used a statistical learning framework to test a range of model complexities to identify suitable hyperparameters of the six BRT metamodels corresponding to each response variable of interest: NO 3 − source concentration factor (which determines the local NO 3 − input concentration); unsaturated zone travel time; NO 3 − concentration at the water table in 1980, 2000, and 2020 (three separate metamodels); and NO 3 − “extinction depth”, the eventual steady state depth of the NO 3 − front. The final metamodels were trained to 129 wells within the active numerical flow model area, and considered 58 mappable predictor variables compiled in a geographic information system (GIS). These metamodels had training and cross-validation testing R 2 values of 0.52 – 0.86 and 0.22 – 0.38, respectively, and predictions were compiled as maps of the above response variables. Testing performance was reasonable, considering that we limited the metamodel predictor variables to mappable factors as opposed to using all available VFM input variables. Relationships between metamodel predictor variables and mapped outputs were generally consistent with expectations, e.g. with greater source concentrations and NO 3 − at the groundwater table in areas of intensive crop use and well drained soils. Shorter unsaturated zone travel times in poorly drained areas likely indicated preferential flow through clay soils, and a tendency for fine grained deposits to collocate with areas of shallower water table. Numerical estimates of groundwater recharge were important in the metamodels and may have been a proxy for N input and redox conditions in the northern FWP, which had shallow predicted NO 3 − extinction depth. The metamodel results provide proof-of-concept for regional characterization of unsaturated zone NO 3 − transport processes in a statistical framework based on readily mappable GIS input variables.

Wisconsin↗

USGS2014

Decadal surface water quality trends under variable climate, land use, and hydrogeochemical setting in Iowa, USA

Understanding how nitrogen fluxes respond to changes in agriculture and climate is important for improving water quality. In the midwestern United States, expansion of corn cropping for ethanol production led to increasing N application rates in the 2000s during a period of extreme variability of annual precipitation. To examine the effects of these changes, surface water quality was analyzed in 10 major Iowa Rivers. Several decades of concentration and flow data were analyzed with a statistical method that provides internally consistent estimates of the concentration history and reveals flow-normalized trends that are independent of year-to-year streamflow variations. Flow-normalized concentrations of nitrate+nitrite-N decreased from 2000 to 2012 in all basins. To evaluate effects of annual discharge and N loading on these trends, multiple conceptual models were developed and calibrated to flow-weighted annual concentrations. The recent declining concentration trends can be attributed to both very high and very low discharge in the 2000s and to the long (e.g., 8 year) subsurface residence times in some basins. Dilution of N and depletion of stored N occurs in years with high discharge. Reduced N transport and increased N storage occurs in low-discharge years. Central Iowa basins showed the greatest reduction in flow-normalized concentrations, likely because of smaller storage volumes and shorter residence times. Effects of land-use changes on the water quality of major Iowa Rivers may not be noticeable for years or decades in peripheral basins of Iowa, and may be obscured in the central basins where extreme flows strongly affect annual concentration trends.

Iowa↗

USGS2012

Factors controlling nitrate fluxes in groundwater in agricultural areas

The impact of agricultural chemicals on groundwater quality depends on the interactions of biogeochemical and hydrologic factors. To identify key processes affecting distribution of agricultural nitrate in groundwater, a parsimonious transport model was applied at 14 sites across the U.S. Simulated vertical profiles of NO 3 - , N 2 from denitrification, O 2 , Cl - , and environmental tracers of groundwater age were matched to observations by adjusting the parameters for recharge rate, unsaturated zone travel time, fractions of N and Cl - inputs leached to groundwater, O 2 reduction rate, O 2 threshold for denitrification, and denitrification rate. Model results revealed important interactions among biogeochemical and physical factors. Chloride fluxes decreased between the land surface and water table possibly because of Cl - exports in harvested crops (averaging 22% of land-surface Cl - inputs). Modeled zero-order rates of O 2 reduction and denitrification were correlated. Denitrification rates at depth commonly exceeded overlying O 2 reduction rates, likely because shallow geologic sources of reactive electron donors had been depleted. Projections indicated continued downward migration of NO 3 - fronts at sites with denitrification rates <0.25 mg-N L -1 yr -1 . The steady state depth of NO 3 - depended to a similar degree on application rate, leaching fraction, recharge, and NO 3 - and O 2 reaction rates. Steady state total mass in each aquifer depended primarily on the N application rate. In addition to managing application rates at land surface, efficient water use may reduce the depth and mass of N in groundwater because lower recharge was associated with lower N fraction leached. Management actions to reduce N leaching could be targeted over aquifers with high-recharge and low-denitrification rates.

Water Resources Research↗