Search USGS⌕ Search
Home/Archive/USGS

USGS · 70168426

Extending the MODPATH algorithm to rectangular unstructured grids

David W. Pollock
View original source ↗Export metadata ↓

Abstract

The recent release of MODFLOW-USG, which allows model grids to have irregular, unstructured connections, requires a modification of the particle-tracking algorithm used by MODPATH. This paper describes a modification of the semi-analytical particle-tracking algorithm used by MODPATH that allows it to be extended to rectangular-based unstructured grids by dividing grid cells with multi-cell face connections into sub-cells. The new method will be incorporated in the next version of MODPATH which is currently under development.

Explore related subjects

Ground WaterArticle

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

David W. Pollock. 2015-03-05. Extending the MODPATH algorithm to rectangular unstructured grids. https://doi.org/10.1111/gwat.12328

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

KEEP EXPLORING

Related USGS reports

USGS2018

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water↗
USGS2016

A semi-structured MODFLOW-USG model to evaluate local water sources to wells for decision support

In order to better represent the configuration of the stream network and simulate local groundwater-surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater-flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface-water depletion by wells) without recourse to labor-intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW-USG) allows application of regional models to address local problems. A “semi-structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20-layer model with uniform 5000-foot (1524-m) lateral spacing was converted to 4 layers with 500-foot (152-m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi-structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi-structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface-water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).

Illinois, Indiana, Michigan, Wisconsin↗
USGS2015

PHT3D-UZF: A reactive transport model for variably-saturated porous media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

Ground Water↗