MODFLOW-NWT, a Newton formulation for MODFLOW-2005

Abstract

For more information on the USGS-the Federal source for science about the Earth, its natural and living resources, natural hazards, and the environment, visit http://www.usgs.gov or call 1-888-ASK-USGS.For an overview of USGS information products, including maps, imagery, and publications, visit http://www.usgs.gov/pubprodTo order this and other USGS information products, visit http://store.usgs.gov Any use of trade, product, or firm names is for descriptive purposes only and does not imply endorsement by the U.S. Government.Although this report is in the public domain, permission must be secured from the individual copyright owners to reproduce any copyrighted materials contained within this report. Altitude, as used in this report, refers to distance above the vertical datum.*Transmissivity: The standard unit for transmissivity is cubic foot per day per square foot times foot of aquifer thickness [(ft 3 /d)/ft 2 ]ft. In this report, the mathematically reduced form, foot squared per day (ft 2 /d), is used for convenience. The UPW Package treats nonlinearities of cell drying and rewetting by use of a continuous function of groundwater head, rather than the discrete approach of drying and rewetting that is used by the BCF, LPF, and HUF Packages. This further enables application of the Newton formulation for unconfined groundwater-flow problems because conductance derivatives required by the Newton method are smooth over the full range of head for a model cell. TablesThe NWT linearization approach generates an asymmetric matrix, which is different from the standard MODFLOW formulation that generates a symmetric matrix. Because all linear solvers presently available for use with MODFLOW-2005 solve only symmetric matrices, MODFLOW-NWT includes two previously developed asymmetric matrix-solver options. The matrix-solver options include a generalized-minimum-residual (GMRES) Solver and an Orthomin / stabilized conjugate-gradient (CGSTAB) Solver. The GMRES Solver is documented in a previously published report, such that only a brief description and input instructions are provided in this report. However, the CGSTAB Solver (called χMD) is documented in this report.Flow-property input for the UPW Package is designed based on the LPF Package and material-property input is identical to that for the LPF Package except that the rewetting and vertical-conductance correction options of the LPF Package are not available with the UPW Package. Input files constructed for the LPF Package can be used with slight modification as input for the UPW Package. This report presents the theory and methods used by MODFLOW-NWT, including the UPW Package. Additionally, this report provides comparisons of the new methodology to analytical solutions of groundwater flow and to standard MODFLOW-2005 results by use of an unconfined aquifer MODFLOW example problem. The standard MODFLOW-2005 simulation uses the LPF Package with the wet/dry option active. A new example problem also is presented to demonstrate MODFLOW-NWT's ability to provide a sol...

What controls the partitioning between baseflow and mountain block recharge in the Qinghai‐Tibet Plateau?

“…The groundwater model was created using MODFLOW-NWT [Niswonger et al, 2011] in steady state mode, which assumes that the groundwater flow system is at long-term equilibrium conditions. The 90 m Shuttle Radar Topography Mission digital elevation model data for the HRB was used to define the land surface topography and hydrostratigraphy.…”

Modeling Groundwater Flow in Heterogeneous Porous Media with YAGMod

“…They give rise to numerical difficulties in the presence of dry cells or elements, e.g., under the influence of an extraction source term. When a cell becomes dry, i.e., its calculated water level falls below the bottom of the cell, two main problems arise [3][4][5][6]. First, the dry cell cannot receive external water, if it is declared as "inactive"; neither can it contain any extraction source term, unless it is rewetted, i.e., the water level rises above a prescribed threshold.…”

“…The innovative idea of Keating and Zyvoloski [4] is a weak scaling for vertical connectivity, from partially-saturated to dry conditions. On the other hand, Niswonger et al [6] use a quadratic approximation of the function that relates horizontal conductance to hydraulic head, over small intervals close to the fully-dry and fully-saturated limits.…”

GSFLOW-GRASS v1.0.0: GIS-enabled hydrologic modeling of coupled groundwater–surface-water systems

“…Other modifications are noted in the GSFLOW manual (see sections on "Changes to PRMS" and "Changes to MODFLOW-2005"). An additional feature starting in version 1.2.0 that is not described in the original manual is the inclusion of MODFLOW-NWT (Niswonger et al, 2011), a more numerically robust update to MODFLOW-2005 for groundwater flow.…”