Search USGSSearch

USGS · 70248311

Benefits and cautions in data assimilation strategies: An example of modeling groundwater recharge

Abstract

Assimilating recent observations improves model outcomes for real-time assessments of groundwater processes. This is demonstrated in estimating time-varying recharge to a shallow fractured-rock aquifer in response to precipitation. Results from estimating the time-varying water-table altitude ( h ) and recharge, and their error covariances, are compared for forecasting, filtering, and fixed-lag smoothing (FLS), which are implemented using the Kalman Filter as applied to a data-driven, mechanistic model of recharge. Forecasting uses past observations to predict future states and is the current paradigm in most groundwater modeling investigations; filtering assimilates observations up to the current time to estimate current states; and FLS estimates states following a time lag over which additional observations are collected. Results for forecasting yield a large error covariance relative to the magnitude of the expected recharge. With assimilating recent observations of h , filtering and FLS produce estimates of recharge that better represent time-varying observations of h and reduce uncertainty in comparison to forecasting. Although model outcomes from applying data assimilation through filtering or FLS reduce model uncertainty, they are not necessarily mass conservative, whereas forecasting outcomes are mass conservative. Mass conservative outcomes from forecasting are not necessarily more accurate, because process errors are inherent in any model. Improvements in estimating real-time groundwater conditions that better represent observations need to be weighed for the model application against outcomes with inherent process deficiencies. Results from data assimilation strategies discussed in this investigation are anticipated to be relevant to other groundwater processes models where system states are sensitive to system inputs.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Allen M. Shapiro, Frederick Day-Lewis. 2023-09-15. Benefits and cautions in data assimilation strategies: An example of modeling groundwater recharge. https://doi.org/10.1111/gwat.13349

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

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota