11. Advanced Packages with MVR

The Advanced MODFLOW 6 Packages example problem is designed to demonstrate the combined use of the UZF, SFR and LAK Packages in MODFLOW 6. It originally appeared in the “Getting Started” pdf document in the MODFLOW 6-docs.git repository, and is included here for completeness of the examples.

Flows exchanged between the advanced packages (i.e., UZF, SFR, and LAK) and the WEL package are made with the MVR Package. The problem was adapted from a previously developed test problem documented in (Prudic et al., 2004) and (Richard G. Niswonger & Prudic, 2005) and represents a developed basinfill aquifer in the northern Great Basin. Streamflow enters the valley from three perennial streams. Because the valley is much lower in altitude than the surrounding uplands, recharge in the valley floor is relatively low as compared to recharge from streams and lakes that receive flow from precipitation falling in the mountains. In the lowland areas, native plants use the ground water. Two lakes were added to the original test problem to illustrate the capabilities of the Lake Package (Figure 11.1).

../_images/ex-gwf-sfr-p01b.png

Figure 11.1 Plan view of lakes and stream reaches in the MODFLOW 6 model that demonstrates usage of the MVR package with the advanced packages (i.e., UZF, SFR, LAK) and includes a linkage with the WEL package.

The topography, and the locations of streams and lakes within the basin are shown in Figure 11.1. Note that unlike the streamflow routing package for MODFLOW-2005, the SFR Package for MODFLOW 6 does not support designation of stream segments. Stream networks now consist solely of reaches that are numbered consecutively from 1 to the total number of reaches, and all stream property information is input on a reach basis. The simulation also includes and general head boundary conditions (GHB Package) and 10 agricultural wells. Aquifer properties are provided in Table 11.1.

Table 11.1 Model parameters for example ex-gwf-sfr-p01b.

Parameter

Value

Number of periods

24

Number of layers

2

Number of rows

15

Number of columns

10

Column width (\(ft\))

5000.0

Row width (\(ft\))

5000.0

Starting head (\(ft\))

varies

Hydraulic conductivity near the stream (\(ft/s\))

0.002

Hydraulic conductivity in the basin (\(ft/s\))

varies

Lakebed leakance (\(1/s\))

2e-9

Specific storage (\(1/s\))

0.1e-5

Specific yield near the stream (unitless)

0.2

Specific yield in the basin (unitless)

0.1

Evapotranspiration rate (\(ft/s\))

9.5e-8

Evapotranspiration extinction depth (\(ft\))

15.0

11.1. Example Results

Unlike the original UZF1 Package (R. G. Niswonger et al., 2006) for MODFLOW-2005, the UZF Package for MODFLOW 6 can simulate unsaturated flow separately for each layer. Thus, rather than simulating homogenous unsaturated flow between the water table and land surface, the unsaturated zone in cells in different MODFLOW layers can be specified with different unsaturated hydraulic properties to represent vertical heterogeneity. Results from the UZF Package are shown below for the cell at row 5, column 2, layers 1 and 2 (the infiltration is shown for the cell in layer 1 while the recharge occurs in layer 2. The moisture content is for the cell located in layer 1). The unsaturated-zone package is ideal for simulating the delay between when water infiltrates at land-surface and recharges the saturated zone.

../_images/ex-gwf-sfr-p01b-uz.png

Figure 11.2 Simulated infiltration and recharge to and from the unsaturated-zone (left y-axis), respectively, for the cells located at row 5 column 2, including the moisture content of the cell in layer 1 (right y-axis).

Additionally, after calculating any rejected infiltration based on the user-specified infiltration rate, the unsaturated-zone package will partition the net infiltration into evapotranspiration (ET), changes in soil moisture storage, and recharge as mentioned above. For this problem, the MVR package is used to transfer pumped groundwater (WEL package) to land surface for simulating irrigation with groundwater. Figure 11.3 shows the monthly totals of pumped water used for irrigation for each month of the simulation. Figure 11.3 also shows the total runoff, primarily groundwater discharge to land surface (spring flow) that is transferred to the stream network (or lakes) using the MVR package.

../_images/ex-gwf-sfr-p01b-mvr.png

Figure 11.3 Monthly totals of pumped water from all 10 simulated wells transferred to UZF cells via the MVR package. Total monthly runoff amounts transferred to the surface water network are also shown.

11.2. References Cited

Niswonger, Richard G., & Prudic, D. E. (2005). Documentation of the Streamflow-Routing (SFR2) package to include unsaturated flow beneath streams—a modification to SFR1. Retrieved from https://pubs.er.usgs.gov/publication/tm6A13

Niswonger, R. G., Prudic, D. E., & Regan, R. S. (2006). Documentation of the Unsaturated-Zone Flow (UZF1) package for modeling unsaturated flow between the land surface and the water table with MODFLOW-2005. Retrieved from https://pubs.usgs.gov/tm/2006/tm6a19/

Prudic, D. E., Konikow, L. F., & Banta, E. R. (2004). A new streamflow-routing (SFR1) package to simulate stream-aquifer interaction with MODFLOW-2000. Retrieved from https://pubs.er.usgs.gov/publication/ofr20041042

11.3. Jupyter Notebook

The Jupyter notebook used to create the MODFLOW 6 input files for this example and post-process the results is: