USGS ScienceSearch

USGS · 70031127

A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids

Abstract

A flexible, robust method for linking parent (regional-scale) and child (local-scale) grids of locally refined models that use different numerical methods is developed based on a new, iterative ghost-node method. Tests are presented for two-dimensional and three-dimensional pumped systems that are homogeneous or that have simple heterogeneity. The parent and child grids are simulated using the block-centered finite-difference MODFLOW and control-volume finite-element FEHM models, respectively. The models are solved iteratively through head-dependent (child model) and specified-flow (parent model) boundary conditions. Boundary conditions for models with nonmatching grids or zones of different hydraulic conductivity are derived and tested against heads and flows from analytical or globally-refined models. Results indicate that for homogeneous two- and three-dimensional models with matched grids (integer number of child cells per parent cell), the new method is nearly as accurate as the coupling of two MODFLOW models using the shared-node method and, surprisingly, errors are slightly lower for nonmatching grids (noninteger number of child cells per parent cell). For heterogeneous three-dimensional systems, this paper compares two methods for each of the two sets of boundary conditions: external heads at head-dependent boundary conditions for the child model are calculated using bilinear interpolation or a Darcy-weighted interpolation; specified-flow boundary conditions for the parent model are calculated using model-grid or hydrogeologic-unit hydraulic conductivities. Results suggest that significantly more accurate heads and flows are produced when both Darcy-weighted interpolation and hydrogeologic-unit hydraulic conductivities are used, while the other methods produce larger errors at the boundary between the regional and local models. The tests suggest that, if posed correctly, the ghost-node method performs well. Additional testing is needed for highly heterogeneous systems. ?? 2007 Elsevier Ltd. All rights reserved.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

J.E. Dickinson, S.C. James, S. Mehl, M. C. Hill, S. A. Leake, G.A. Zyvoloski, C.C. Faunt, A.-A. Eddebbarh. 2007. A new ghost-node method for linking different models and initial investigations of heterogeneity and nonmatching grids. https://doi.org/10.1016/j.advwatres.2007.01.004

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

KEEP EXPLORING

Related USGS reports

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

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

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