Search USGS⌕ Search
Home/Archive/USGS

USGS · 70026804

Three-dimensional local grid refinement for block-centered finite-difference groundwater models using iteratively coupled shared nodes: A new method of interpolation and analysis of errors

S. Mehl·M. C. Hill
View original source ↗Export metadata ↓

Abstract

This paper describes work that extends to three dimensions the two-dimensional local-grid refinement method for block-centered finite-difference groundwater models of Mehl and Hill [Development and evaluation of a local grid refinement method for block-centered finite-difference groundwater models using shared nodes. Adv Water Resour 2002;25(5):497-511]. In this approach, the (parent) finite-difference grid is discretized more finely within a (child) sub-region. The grid refinement method sequentially solves each grid and uses specified flux (parent) and specified head (child) boundary conditions to couple the grids. Iteration achieves convergence between heads and fluxes of both grids. Of most concern is how to interpolate heads onto the boundary of the child grid such that the physics of the parent-grid flow is retained in three dimensions. We develop a new two-step, "cage-shell" interpolation method based on the solution of the flow equation on the boundary of the child between nodes shared with the parent grid. Error analysis using a test case indicates that the shared-node local grid refinement method with cage-shell boundary head interpolation is accurate and robust, and the resulting code is used to investigate three-dimensional local grid refinement of stream-aquifer interactions. Results reveal that (1) the parent and child grids interact to shift the true head and flux solution to a different solution where the heads and fluxes of both grids are in equilibrium, (2) the locally refined model provided a solution for both heads and fluxes in the region of the refinement that was more accurate than a model without refinement only if iterations are performed so that both heads and fluxes are in equilibrium, and (3) the accuracy of the coupling is limited by the parent-grid size - A coarse parent grid limits correct representation of the hydraulics in the feedback from the child grid.

Explore related subjects

Advances in Water ResourcesArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

S. Mehl, M. C. Hill. 2004. Three-dimensional local grid refinement for block-centered finite-difference groundwater models using iteratively coupled shared nodes: A new method of interpolation and analysis of errors. https://doi.org/10.1016/j.advwatres.2004.06.004

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

KEEP EXPLORING

Related USGS reports

USGS2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen–Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen–Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Advances in Water Resources↗
USGS2025

The acoustic-Doppler current profiler (ADCP): A comprehensive tool for river-reach hydromorphodynamics

This paper introduces the use of acoustic Doppler current profiler (ADCP) measurements as input for the Acoustic Mapping Velocimetry (AMV) method, a technique for characterizing the dynamics of riverine bedforms. The performance of this new approach, ADCP-AMV, is compared with input from a multibeam echosounder through a field study conducted on the Mississippi River (USA). A virtual ADCP tool has been created to support the ADCP-AMV measurements with optimal data density predictions. To the authors’ knowledge, this is the first time ADCP measurements have been used in conjunction with the AMV dune-tracking method. Subsequently, the paper discusses the coupling of ADCP-AMV measurements with ancillary data extracted from the ADCP. These ancillary data are processed using previously developed protocols to characterize hydrodynamics and the suspended sediment distribution in the water column. This paper emphasizes the capability of ADCPs to characterize open-channel river hydromorphodynamic parameters with high spatiotemporal resolution. Recommendations to accurately and efficiently acquire these multi-variable measurements and derived datasets are discussed.

Tennessee↗
USGS2024

Precipitation uncertainty estimation and rainfall-runoff model calibration using iterative ensemble smoothers

The introduction of iterative ensemble smoothers (IES) for parameter calibration opens avenues for expanding parameter space in surface water hydrologic modeling. Here, we have introduced independent parameters into a model calibration experiment to estimate errors in rainfall forcing data. This approach has the potential to estimate rainfall errors using other hydrological observations and to improve model calibration. Using high-resolution rain gauge data, we estimated “real” rainfall errors across the Turkey River watershed at storm and daily scales. Tests on synthetic and real-world scenarios successfully estimated errors correlated with observed values – even at daily scales. However, a bias remained from model parameter compensation, and identifying errors was challenging for low precipitation and snowfall. Despite synthetic results showing good error correlation, the biases in parameter identification masked potential improvements in hydrological calibration. This study highlights the potential of IES to provide additional information on rainfall errors, even only using streamflow observations.

Advances in Water Resources↗