Model documentation report: Smith, S.J., Ellis, J.H., Wagner, D.L., and Peterson, S.M., 2017, Hydrogeology and simulated groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013: U.S. Geological Survey Scientific Investigations Report 2017-5098, https://doi.org/10.3133/sir20175098 Model data release: Smith, S.J., Ellis, J.H., Wagner, D.L., and Peterson, S.M., 2017, MODFLOW-NWT model used in simulation of groundwater flow and availability in the North Fork Red River aquifer, southwest Oklahoma, 1980–2013: U.S. Geological Survey data release, https://doi.org/10.5066/F7JQ0ZXH. ------------------------------------------------------------------------------------------------------------------ GENERAL The model directory and subdirectories contain all of the input and output files for the simulations described in the report. About 33 GB of space is required to unzip the entire model archive and run all of the simulations. The steady-state model represents the average 1980-2013 conditions, and is included as stress period 1. The transient model runs 1980-2013. PYTHON All of the groundwater-flow models in this archive require Python 2.X and the numpy module to run successfully. These are available at this link in a single package: https://www.continuum.io. Once installed, python must be added to the system path variable under 'environmental variables.' To test if this is needed, open a command prompt (Winkey+R, type 'cmd') and type 'python' at the prompt. If an error results, python needs to be added to the system path variable. To run the Parameter Estimation (PEST) python script (sens_table.py), in addition to the requirements stated above, the PESTools and pandas packages are required, which requires Python 2.7. The PESTools package can be downloaded here: https://github.com/PESTools/pestools. The Pandas package can be downloaded here: http://pandas.pydata.org/getpandas.html. MODEL INFORMATION Length Unit: Meters Steady-State Stress Period length: 365 days Transient Stress Period length: 1 month (28 to 31 days) Transient Stress Period duration: January 1980-December 2013 Cell size: 270 meters uniform grid spacing Time steps per stress period: 2 Layers: 2 (Layer 1 is Quaternary-age alluvium and terrace deposits, Layer 2 is Permian-age bedrock) Output Control: Displays outputs for Time Step 2 in each Stress Period Solver: NWT System requirements The models contained in this data release were run using the executables in the bin.zip file in this data release. The simulations were run on a personal computer with 32GB Random Access Memory with 64-bit Windows 7 OS. The models will take approximately 1 GB of memory each to run. Support Support is provided for correcting errors in the data release and clarification of the modeling conducted by the U.S. Geological Survey. Users are encouraged to review the model documentation report (https://doi.org/10.3133/sir20175098) to understand the purpose, construction, and limitations of this model. ------------------------------------------------------------------------------------------------------------------ In order to run the models, all zipped folders in the model archive must be unzipped in accordance with the directory structure shown below. DIRECTORY DESCRIPTION SIR2017-5098 \ancillary \bin \georef \model \current_head \current_use_2013pump \current_use_avgpump \current_use_growthpump \current_use_nopump \drought \externalfiles \arrays \EVT \RCH \calibration_targets \heads \stagevolume \final_model \may_20_year \may_20_year_dec \may_20_year_inc \may_20_year_storage \may_40_year \may_40_year_dec \may_40_year_inc \may_50_year \may_50_year_dec \may_50_year_inc \may_head \output \output.current_head \output.current_use_2013pump \output.current_use_avgpump \output.current_use_growthpump \output.current_use_nopump \output.drought \output.final_model \output.may_20_year \output.may_20_year_dec \output.may_20_year_inc \output.may_20_year_storage \output.may_40_year \output.may_40_year_dec \output.may_40_year_inc \output.may_50_year \output.may_50_year_dec \output.may_50_year_inc \output.may_head \source \MODFLOW-NWT_1.0.8 \Zonbud.3_01 ------------------------------------------------------------------------------------------------------------------ SUBDIRECTORY AND FILE DESCRIPTIONS \ancillary.zip: This directory contains the Soil-Water-Balance (SWB) code used to compute and distribute recharge, a polygon shapefile of the MODFLOW model grid, final PEST files used to calculate the model sensitivity, a python script (Calc_storage.py, with batch files) to calculate storage, a python script (may.py) to compute equal-proportionate-share pumping rates (with initialization .ini files). Calc_storage.py: Python script used to calculate groundwater in storage Calc_storage_20EPS.bat: Batch file used to calculate groundwater storage for the 20-year equal proportionate share (EPS scenario) Calc_storage_2013pump.bat: Batch file used to calculate groundwater storage for the projected (50-year) pumping scenario with 2013 rates Calc_storage_avgpump.bat: Batch file used to calculate groundwater storage for the projected (50-year) pumping scenario with mean rates for the study period Calc_storage_growthpump.bat: Batch file used to calculate groundwater storage for the projected (50-year) pumping scenario with increased demand rates Calc_storage_nopump.bat: Batch file used to calculate groundwater storage for the projected (50-year) pumping scenario with no pumping Calc_storage_drought.bat: Batch file used to calculate groundwater storage for the drought scenario Calc_storage_final_model.bat: Batch file used to calculate groundwater storage for the final model may.py: This script calculates the number of cells at the end of each simulation to determine whether 50 percent of the cells have a saturated thickness of 1.524 meters (5 feet) or less. nfrr_bpa3.jco: Jacobian matrix file from the PEST calibration nfrr_bpa3.pst: PEST control file nfrr_bpa3.rec: PEST run record file nfrr_bpa3.rei: PEST interim observation residuals file NFRR_may_xx.ini: Contains command line information read in by the may.py python script, where 'xx' indicates the 20-, 40-, or 50-year time period. The 'inc' or 'dec' suffix indicates either recharge increased or decreased by 10 percent, respectively. pargroupsnoobs.txt: Contains a list of the parameter groups to be read in by the sens_table.py script. sens_table.py: This python script multiplies each observation weight by the derivatives of observations with respect to parameters (the JCO), strips the regularization groups, and then sums the derivatives of observations for each parameter together. It then sums the sensitivities of all the parameters in each group. sensbygroupsumout.txt: Parameter sensitivities calculated by the sens_table.py script. SIR2017_5098_ModelGrid.dbf: Attribute information for a polygon shapefile of the model grid SIR2017_5098_ModelGrid.prj: Coordinate system information for a polygon shapefile of the model grid SIR2017_5098_ModelGrid.shp: Feature geometry information for a polygon shapefile of the model grid SIR2017_5098_ModelGrid.shx: Feature geometry index information for a polygon shapefile of the model grid SIR2017_5098_NorthForkRedRiverAquiferExtent.dbf: Attribute information for a polygon shapefile of the modified North Fork Red River aquifer extent SIR2017_5098_NorthForkRedRiverAquiferExtent.prj: Coordinate system information for a polygon shapefile of the modified North Fork Red River aquifer extent SIR2017_5098_NorthForkRedRiverAquiferExtent.shp: Feature geometry information for a polygon shapefile of the modified North Fork Red River aquifer extent SIR2017_5098_NorthForkRedRiverAquiferExtent.shx: Feature geometry index information for a polygon shapefile of the modified North Fork Red River aquifer extent \bin.zip: This directory contains compiled executable codes for MODFLOW-NWT, version 1.0.8 and Zonebudget 3.01. MODFLOW-NWT_64.exe: Compiled MODFLOW-NWT executable code zonebud.exe: Compiled ZONEBUDGET executable code \georef.zip: This directory contains a polygon shapefile showing the active and inactive areas of the model domain. SIR2017_5098.dbf: Attribute information for a polygon shapefile showing the active and inactive areas of the model domain SIR2017_5098.prj: Coordinate system information for a polygon shapefile showing the active and inactive areas of the model domain SIR2017_5098.shp: Feature geometry information for a polygon shapefile showing the active and inactive areas of the model domain SIR2017_5098.shx: Feature geometry index information for a polygon shapefile showing the active and inactive areas of the model domain \model.zip: This directory contains all of the files specific to each scenario (e.g. current use, drought, final model, may) and common files used by all models (externalfiles). IMPORTANT: The calibrated model (e.g. final_model folder) needs to be run first before any of the other scenarios are run. See the 'RUNNING THE MODELS' section for more details. \current_head: Holds files to produce input (starting) heads for the projected (50-year) pumping scenarios. The input files in this folder and identical to those in the 'final_model' folder, with the exception of the 'current_head.bat', 'NFRR_current.nam', and 'NFRR_sp1.oc' files, described below. current_head.BAT: Batch file for running the calibrated model to produce input heads for the projected (50-year) pumping scenarios NFRR.lak: Lake Package input file NFRR.rch: Recharge Package input file NFRR.sfr: Streamflow-Routing Package input file NFRR.wel: Well Package input file NFRR_current.nam: Name file for running the calibrated model to produce input heads for the projected (50-year) pumping scenarios NFRR_sp1.oc: Output Control file that saves heads for only the first stress period (steady-state) usgs.model.reference: File contains data to register the model in space and time \current_use_2013pump: Holds files for the projected (50-year) pumping scenario using 2013 rates NFRR_current_50_2013pump.BAT: Batch file for running the 50-year use scenario with 2013 rates NFRR_current_50_2013pump.nam: Name file for the 50-year use scenario with 2013 rates NFRR_current_50_2013.wel: Well Package input file for 2013 rates usgs.model.reference: File contains data to register the model in space and time \current_use_avgpump: Holds files for the projected (50-year) pumping scenario using mean rates for the study period NFRR_current_50_avgpump.BAT: Batch file for running the 50-year use scenario with mean rates for the study period NFRR_current_50_avgpump.nam: Name file for the 50-year use scenario with mean rates for the study period NFRR_current_50_avg.wel: Well Package input file for mean rates usgs.model.reference: File contains data to register the model in space and time \current_use_growthpump: Holds files for the projected (50-year) pumping scenario using increased demand rates NFRR_current_50_growthpump.BAT: Batch file for running the 50-year use scenario with increased demand rates NFRR_current_50_growthpump.nam: Name file for the 50-year use scenario with increased demand rates NFRR_current_50_growth.wel: Well Package input file for increased demand rates usgs.model.reference: File contains data to register the model in space and time \current_use_nopump: Holds files for the projected (50-year) pumping scenario with no pumping NFRR_current_50_nopump.BAT: Batch file for running the 50-year use scenario with no pumping NFRR_current_50_nopump.nam: Name file for the 50-year use scenario with no pumping usgs.model.reference: File contains data to register the model in space and time \drought: Holds files for the hypothetical (10-year) drought scenario NFRR_drought.BAT: Batch file for running the hypothetical (10-year) drought scenario NFRR_drought.nam: Name file for the hypothetical (10-year) drought scenario NFRR_drought.lak: Lake Package input file NFRR_drought.rch: Recharge Package input file NFRR_drought.sfr: Streamflow-Routing Package input file NFRR_drought.wel: Well Package input file usgs.model.reference: File contains data to register the model in space and time \externalfiles: Holds common files used by all of the models Altus_stagevolume.dat: Text files containing relations among lake stage, surface area, and volume Altus_stagevolume_MAY.dat: Text files containing relations among lake stage, surface area, and volume for the EPS scenarios. Has an additional row for 0 surface area and volume to avoid the lake package failing to converge during the EPS scenarios when the lake goes dry. head1979c.asc: Starting head for the steady state stress period NFRR.ba6: Basic Package file NFRR.chd: Time-Variant Specified-Head input file NFRR.dis: Discretization file NFRR.drn: Drain Package input file NFRR.evt: Evapotranspiration Package input file NFRR.gag: Gage Package input file NFRR.ghb: General-Head Boundary Package input file NFRR.hob: Head observation package file NFRR.nwt: Newton Solver file NFRR.oc: Output Control file NFRR.upw: Upstream Weighting Package file NFRR_current_50.ba6: Basic package file NFRR_current_50.chd: Time-Variant Specified-Head input file NFRR_current_50.dis: Discretization file NFRR_current_50.drn: Drain Package input file NFRR_current_50.evt: Evapotranspiration Package input file NFRR_current_50.ghb: General-Head Boundary Package input file NFRR_current_50.lak: Lake Package input file NFRR_current_50.oc: Output Control file NFRR_current_50.rch: Recharge Package input file NFRR_current_50.sfr: Streamflow-Routing Package input file NFRR_may.ba6: Basic Package file NFRR_may_20.dis: Discretization file NFRR_may_20.oc: Output Control file NFRR_may_40.dis: Discretization file NFRR_may_40.oc: Output Control file NFRR_may_50.chd: Time-Variant Specified-Head input file NFRR_may_50.dis: Discretization file NFRR_may_50.drn: Drain Package input file NFRR_may_50.evt: Evapotranspiration Package input file NFRR_may_50.ghb: General-Head Boundary Package input file NFRR_may_50.lak: Lake Package input file NFRR_may_50.oc: Output Control file NFRR_may_50.rch: Recharge Package input file NFRR_may_50.sfr: Streamflow-Routing Package input file NFRR_may_50_dec.rch: Recharge Package input file NFRR_may_50_inc.rch: Recharge Package input file NFRR_sp1.hds: Stress period 1 output heads for the projected (50-year) pumping scenarios NFRR_sp409.hds: Stress period 409 output heads for the EPS scenarios (which calculate the maximum annual yield, or MAY) TwoArray_Kh.py: Python script that takes the horizontal hydraulic conductivity array from the results of the hydrogeologic framework (kh_1.asc), reshapes it to a single column, multiplies by a calibrated K multiplier array (kh_mult_r1_1D.asc), then reshapes the multiplied column back to an array to be read into the model(kh_1_mult.asc). \arrays: Holds common array files used by all of the models et_extdpth.asc: Evapotranspiration package extinction depth array ib_1.asc: Boundary variable array indicating active or inactive cells k1bot.asc: Bottom elevation of layer 1 k1topfil.asc: Top elevation of layer 1 k2bot.asc: Bottom elevation of layer 2 kh_1.asc: Horizontal hydraulic conductivity array from the results of the hydrogeologic framework kh_1_mult.asc: Final horizontal hydraulic conductivity array read into the model for each scenario. kh_mult_1d.asc: Calibrated horizontal hydraulic conductivity multiplier array lake.asc: Holds array for lake package data set 5 (LKARR) Sy.asc: Specific yield array for layer 1 zonebudget.zon: ZONEBUDGET Zone File input \EVT: Holds common evapotranspiration rate array files used by all of the models. This folder holds two types of arrays: averaged monthly evapotranspiration rate arrays for the projected (50-year) pumping scenario (current_use folder), and unique monthly evapotranspiration rate arrays for the calibrated model and drought scenario. xxx_EVT.asc: Averaged monthly evapotranspiration rate arrays, where 'xxx' indicates the month. pet_minus_aet_avg80_13_70rz.asc: The 1980 evapotranspiration rate array. pet_minus_aetxxx.asc: The unique monthly evapotranspiration rate arrays, where 'xxx' indicates the stress period. \RCH: Holds common recharge array files used by all of the models. This folder holds two types of arrays: averaged monthly recharge rate arrays for the projected (50-year) pumping scenario (current_use folder), and unique monthly recharge rate arrays for the calibrated model and drought scenario. xxx_RCH.asc: Averaged monthly recharge rate arrays, where 'xxx' indicates the month. rchavgss.asc: Holds the steady-state recharge rate array rchxxx.asc: Monthly recharge rate arrays, where 'xxx' indicates the stress period. \final_model: Holds files for the calibrated model NFRR.BAT: Batch file for running the calibrated model NFRR.nam: Name file for the calibrated model zonebudget.in: ZONEBUDGET input file NFRR.lak: Lake Package input file NFRR.rch: Recharge Package input file NFRR.sfr: Streamflow-Routing Package input file NFRR.wel: Well Package input file usgs.model.reference: File contains data to register the model in space and time \may_20_year: Holds files for the 20-year EPS scenario may_20.bat: Batch file for running the EPS scenario with normal (same) recharge. This batch file runs the may.py python script. NFRR_may_20.nam: Name file for the 20-year EPS scenario. NFRR_may_20.BAT: Batch file for the 20-year EPS scenario. This file is called by the may.py script. NFRR_may_20.wel: Well Package input file for normal (same) recharge EPS scenario usgs.model.reference: File contains data to register the model in space and time \may_20_year_dec: Holds files for the 20-year EPS scenario with recharge decreased 10 percent may_20_dec.bat: Batch file for running the EPS scenario with recharge decreased 10 percent. This batch file runs the may.py python script. NFRR_may_20_dec.BAT: Batch file for the 20-year EPS scenario with recharge decreased 10 percent. This file is called by the may.py script. NFRR_may_20_dec.nam: Name file for the 20-year EPS scenario with recharge decreased 10 percent. input files: The input files are structured the same as the 20-year EPS scenario with normal (same) recharge. usgs.model.reference: File contains data to register the model in space and time \may_20_year_inc: Holds files for the 20-year EPS scenario with recharge increased 10 percent may_20_inc.bat: Batch file for running the EPS scenario with recharge increased 10 percent. This batch file runs the may.py python script. NFRR_may_20_inc.BAT: Batch file for the 20-year EPS scenario with recharge increased 10 percent. This file is called by the may.py script. NFRR_may_20_inc.nam: Name file for the 20-year EPS scenario with recharge increased 10 percent. input files: The input files are structured the same as the 20-year EPS scenario with normal (same) recharge. usgs.model.reference: File contains data to register the model in space and time \may_20_year_storage: Holds files for the 20-year EPS scenario NFRR_may_20.wel: Well Package input file for normal (same) recharge EPS scenario NFRR_may_20_allSP.oc: Output Control File used for groundwater storage calculation NFRR_may_20_calcstorage.nam: Name file for the 20-year EPS scenario. NFRR_may_20_storage.BAT: Batch file for the 20-year EPS scenario. This file is called by the may.py script. usgs.model.reference: File contains data to register the model in space and time \may_40_year: Holds files for the 40-year EPS scenario with normal (same) recharge. The folder structure and input files follow the same format as the 20-year EPS scenario with normal (same) recharge \may_40_year_dec: Holds files for the 40-year EPS scenario with recharge decreased 10 percent. The folder structure and input files follow the same format as the 20-year EPS scenario with recharge decreased 10 percent \may_40_year_inc: Holds files for the 40-year EPS scenario with recharge increased 10 percent. The folder structure and input files follow the same format as the 20-year EPS scenario with recharge increased 10 percent \may_50_year: Holds files for the 50-year EPS scenario with normal (same) recharge. The folder structure and input files follow the same format as the 20-year EPS scenario with normal (same) recharge \may_50_year_dec: Holds files for the 50-year EPS scenario with recharge decreased 10 percent. The folder structure and input files follow the same format as the 20-year EPS scenario with recharge decreased 10 percent \may_50_year_inc: Holds files for the 50-year EPS scenario with recharge increased 10 percent. The folder structure and input files follow the same format as the 20-year EPS scenario with recharge increased 10 percent \may_head: Holds files to produce input (starting) heads for the EPS scenarios. The input files in this folder and identical to those in the 'final_model' folder, with the exception of the 'may_head.bat', 'NFRR_EPS.nam', and 'NFRR_sp409.oc' files, described below. may_head.BAT: Batch file for running the calibrated model to produce input heads for the projected EPS scenarios NFRR.lak: Lake Package input file NFRR_EPS.nam: Name file for running the calibrated model to produce input heads for the EPS scenarios NFRR.rch: Recharge Package input file NFRR.sfr: Streamflow-Routing Package input file NFRR.wel: Well Package input file NFRR_sp409.oc: Output Control file that saves heads for only the first stress period (steady-state) usgs.model.reference: File contains data to register the model in space and time \output.zip: This directory contains all of the output files specific to each scenario (e.g. current use, drought, final model, may) documented in the report \output.current_head: Holds files to produce input (starting) heads for the projected (50-year) pumping scenarios. The output files in this folder and identical to those in the 'final_model' folder, with the exception of the storage and ZONEBUDGET files, and the file described below. NFRR.flw: Streamflow-Routing Package flow file NFRR.hobs: Output head observations NFRR_current.out: List file of model outputs and budgets (not used) USGS07301350.sg: Flow observations for streamgage 07301350 USGS07301420.sg: Flow observations for streamgage 07301420 USGS07301481.sg: Flow observations for streamgage 07301481 USGS07301500.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow.sg: Flow observations into Lake Altus USGS07302500_lkaltus.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500.sg: Flow observations for streamgage 07303500 USGS07304500.sg: Flow observations for streamgage 07304500 USGS07305000.sg: Flow observations for streamgage 07305000 USGS07307010.sg: Flow observations for streamgage 07307010 \output.current_use_2013pump: Holds output files for the projected (50-year) pumping scenario with 2013 rates NFRR_current_50_2013pump.flw: Streamflow-Routing Package flow file NFRR_current_50_2013pump.hds: Unformatted simulated hydraulic heads NFRR_current_50_2013pump.hobs: Head observation file NFRR_current_50_2013pump.out: List file of model outputs and budgets storage_2013pump.asc: Output calculated groundwater storage for the projected (50-year) pumping scenario with 2013 rates USGS07301350_50_2013pump.sg: Flow observations for streamgage 07301350 USGS07301420_50_2013pump.sg: Flow observations for streamgage 07301420 USGS07301481_50_2013pump.sg: Flow observations for streamgage 07301481 USGS07301500_50_2013pump.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow_50_2013pump.sg: Flow observations into Lake Altus USGS07302500_lkaltus_50_2013pump.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500_50_2013pump.sg: Flow observations for streamgage 07303500 USGS07304500_50_2013pump.sg: Flow observations for streamgage 07304500 USGS07305000_50_2013pump.sg: Flow observations for streamgage 07305000 USGS07307010_50_2013pump.sg: Flow observations for streamgage 07307010 \output.current_use_avgpump: Holds output files for the projected (50-year) pumping scenario using mean rates for the study period. The file structure is identical to the 'current_use_2013pump' folder. \output.current_use_growthpump: Holds output files for the projected (50-year) pumping scenario with increased demand rates. The file structure is identical to the 'current_use_2013pump' folder. \output.current_use_nopump: Holds output files for the projected (50-year) pumping scenario with no pumping. The file structure is identical to the 'current_use_2013pump' folder. \output.drought: Holds output files for the hypothetical (10-year) drought scenario drought.asc: Output calculated groundwater storage for the hypothetical (10-year) drought scenario drought.hobs: Output head observations drought.flw: Streamflow-Routing Package flow file drought.hds: Unformatted simulated hydraulic heads drought.out: List file of model outputs and budgets USGS07301350_drought.sg: Flow observations for streamgage 07301350 USGS07301420_drought.sg: Flow observations for streamgage 07301420 USGS07301481_drought.sg: Flow observations for streamgage 07301481 USGS07301500_drought.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow_drought.sg: Flow observations into Lake Altus USGS07302500_lkaltus_drought.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500_drought.sg: Flow observations for streamgage 07303500 USGS07304500_drought.sg: Flow observations for streamgage 07304500 USGS07305000_drought.sg: Flow observations for streamgage 07305000 USGS07307010_drought.sg: Flow observations for streamgage 07307010 \output.final_model: Holds output files for the calibrated model NFRR.flw: Streamflow-Routing Package flow file NFRR.hobs: Output head observations NFRR.out: List file of model outputs and budgets NFRR_storage.asc: Output calculated groundwater storage for the calibrated model USGS07301350.sg: Flow observations for streamgage 07301350 USGS07301420.sg: Flow observations for streamgage 07301420 USGS07301481.sg: Flow observations for streamgage 07301481 USGS07301500.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow.sg: Flow observations into Lake Altus USGS07302500_lkaltus.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500.sg: Flow observations for streamgage 07303500 USGS07304500.sg: Flow observations for streamgage 07304500 USGS07305000.sg: Flow observations for streamgage 07305000 USGS07307010.sg: Flow observations for streamgage 07307010 zb.2.csv: ZONEBUDGET output file with budget terms displayed in rows zb.csv: ZONEBUDGET output file with budget terms displayed in colums zb.zblst: ZONEBUDGET output file with budgets for each zone \output.may_20_year: Holds output files for the 20-year EPS scenario storage_may_20.asc: Output calculated groundwater storage for the 20-year EPS scenario may_20.flw: Streamflow-Routing Package flow file NFRR_may_20.csv: EPS scenarios results. Each file will contain several lines of data, discussed in the 'RUNNING THE MODELS' section. NFRR_may_20.hds: Unformatted simulated hydraulic heads NFRR_may_20.out: List file of model outputs and budgets NFRR_may_20_calcstorage.hds: Unformatted simulated hydraulic heads (not used) NFRR_may_20_calcstorage.out: List file of model outputs and budgets (not used) NFRR_may_20_optWel: EPS pumping rates at 50 percent cell saturation of five feet or less (not used) storage_may_20: Output calculated groundwater storage for the 20-year EPS scenario USGS07301350_may_20.sg: Flow observations for streamgage 07301350 USGS07301420_may_20.sg: Flow observations for streamgage 07301420 USGS07301481_may_20.sg: Flow observations for streamgage 07301481 USGS07301500_may_20.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow_may_20.sg: Flow observations into Lake Altus USGS07302500_lkaltus_may_20.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500_may_20.sg: Flow observations for streamgage 07303500 USGS07304500_may_20.sg: Flow observations for streamgage 07304500 USGS07305000_may_20.sg: Flow observations for streamgage 07305000 USGS07307010_may_20.sg: Flow observations for streamgage 07307010 \output.may_20_year_dec: Holds output files for the 20-year EPS scenario with decreased recharge. The file structure is identical to the 'output.may_20_year' folder other than the removal of the 'Calc_storage' file \output.may_20_year_inc: Holds output files for the 20-year EPS scenario with increased recharge. The file structure is identical to the 'output.may_20_year' folder other than the removal of the 'Calc_storage' file \output.may_20_year_storage: Holds storage calculation for the 20-year EPS scenario may_20.flw: Streamflow-Routing Package flow file NFRR_may_20_calcstorage.cbb: Unformatted cell-by-cell flow budget file NFRR_may_20_calcstorage.ddn: Unformatted drawdown file (not used) NFRR_may_20_calcstorage.hds: Unformatted simulated hydraulic heads (not used) NFRR_may_20_calcstorage.out: List file of model outputs and budgets (not used) storage_may_20: Output calculated groundwater storage for the 20-year EPS scenario USGS07301350_may_20.sg: Flow observations for streamgage 07301350 USGS07301420_may_20.sg: Flow observations for streamgage 07301420 USGS07301481_may_20.sg: Flow observations for streamgage 07301481 USGS07301500_may_20.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow_may_20.sg: Flow observations into Lake Altus USGS07302500_lkaltus_may_20.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500_may_20.sg: Flow observations for streamgage 07303500 USGS07304500_may_20.sg: Flow observations for streamgage 07304500 USGS07305000_may_20.sg: Flow observations for streamgage 07305000 USGS07307010_may_20.sg: Flow observations for streamgage 07307010 \output.may_40_year: Holds output files for the 40-year EPS scenario. The file structure is identical to the 'output.may_20_year' folder with the only changes involving the removal of the 'Calc_storage' file. \output.may_40_year_dec: Holds output files for the 40-year EPS scenario with decreased recharge. The file structure is identical to the 'output.may_20_year_dec' \output.may_40_year_inc: Holds output files for the 40-year EPS scenario with increased recharge. The file structure is identical to the 'output.may_20_year_inc' folder \output.may_50-year: Holds output files for the 50-year EPS scenario. The file structure is identical to the 'output.may_20_year' folder with the only changes involving the removal of the 'Calc_storage' file. \output.may_50_year_dec: Holds output files for the 50-year EPS scenario with decreased recharge. The file structure is identical to the 'output.may_20_year_dec' folder. \output.may_50_year_inc: Holds output files for the 50-year EPS scenario with increased recharge. The file structure is identical to the 'output.may_20_year_inc' folder. \output.may_head: Holds files to produce input (starting) heads for the EPS scenarios. The output files in this folder are identical to those in the 'final_model' folder, with the exception of the storage and ZONEBUDGET files, and the file described below. NFRR.flw: Streamflow-Routing Package flow file NFRR.hobs: Output head observations NFRR_EPS.out: List file of model outputs and budgets USGS07301350.sg: Flow observations for streamgage 07301350 USGS07301420.sg: Flow observations for streamgage 07301420 USGS07301481.sg: Flow observations for streamgage 07301481 USGS07301500.sg: Flow observations for streamgage 07301500 USGS07302000_lkinflow.sg: Flow observations into Lake Altus USGS07302500_lkaltus.sg: File listing a variety of outputs from the Lake Package, including stage, volume, and water budget components USGS07303500.sg: Flow observations for streamgage 07303500 USGS07304500.sg: Flow observations for streamgage 07304500 USGS07305000.sg: Flow observations for streamgage 07305000 USGS07307010.sg: Flow observations for streamgage 07307010 \source.zip: Holds MODFLOW-NWT (v 1.0.8) source and compiled executable code, and Zonebudget (v 3.01) compiled executable code. \MODFLOW-NWT_1.0.8: Holds source Fortran code and documentation for MODFLOW-NWT. \Zonbud.3_01: Holds source Fortran code and documentation for Zonebudget 3.01. modelgeoref.txt: Text file of geo-reference information, including data release title and DOI, publication title and DOI, datum, latitude and longitude in decimal degrees of the corners of a rectangle outlining the model study area. readme.txt: Provides a description of the archive folder structure and of the files contained in each subfolder, instructions for running the model(s), instructions for reconstructing the model archive from the downloaded model files, and a brief description of each model run. ------------------------------------------------------------------------------------------------------------------ RUNNING THE MODELS To run the current use, drought, final model, and MAY simulations, navigate to the folders described below. IMPORTANT: The calibrated model (e.g. final_model folder) needs to be run first, followed by the batch files in the 'current_head' and 'may_head' folders before any other scenarios are run. See the 'Calibrated model' subsection below for details. Running batch files: Each simulation can be run by double clicking on the batch file (e.g.xxx.bat) in the model subdirectories. This will bring up a command prompt window which will record the simulation progress. Once the simulation has finished running, press any key to exit and close the command prompt window. Projected (50-year) pumping scenario (current_use folders): Run the four scenarios with the batch files located in the folders below. When these batch files are run, storage is also automatically calculated at the end of the scenario. An explanation of the file name is below: \current_use_2013pump NFRR_current_50_2013pump.bat: Runs the projected (50-year) pumping scenario with 2013 rates and calculates storage \current_use_avgpump NFRR_current_50_avgpump.bat: Runs the projected (50-year) pumping scenario with mean rates for the study period and calculates storage \current_use_growthpump NFRR_current_50_growthpump.bat: Runs the projected (50-year) pumping scenario with increased demand rates and calculates storage \current_use_nopump NFRR_current_50_nopump.bat: Runs the projected (50-year) pumping scenario with no pumping and calculates storage Hypothetical (10-year) drought scenario (drought folder): Run the drought scenario with the 'NFRR_drought.bat' file. When the batch file is run, storage is also automatically calculated at the end of the scenario. Calibrated model (final_model folder): To run the calibrated model and ZONEBUDGET, and calculate groundwater storage, use the 'NFRR.bat' file. Note that when the transient model finishes, ZONEBUDGET is automatically run, then the Calc_storage python script is automatically run to calculate groundwater storage. Zone 1 is the model area upgradient from Lake Altus, while Zone 2 is the model area downgradient from Lake Altus. Once complete, the batch files in the 'current_head' and 'may_head' folders need to be run. The batch files in these folders will produce the output heads for the EPS scenarios, and the output heads for the projected (50-year) pumping scenario. The output heads are copied from the respective output folder into the heads folder and replace the existing heads files. EPS scenarios (may folder): Run the EPS scenarios with the batch files located in the folders below. Note that when the 20-year EPS scenario is run, storage is also automatically calculated at the end of the scenario. The other batch files in each folder starting with 'NFRR' should not be used, as this is called by the batch files below to run the model. \may_20_year may_20.bat: Runs the 20-year EPS scenario with normal recharge \may_20_year_dec may_20_dec.bat: Runs the 20-year EPS scenario with recharge reduced 10 percent \may_20_year_inc may_20_inc.bat: Runs the 20-year EPS scenario with recharge increased 10 percent \may_40_year may_40.bat: Runs the 40-year EPS scenario with normal recharge \may_40_year_dec may_40_dec.bat: Runs the 40-year EPS scenario with recharge reduced 10 percent \may_40_year_inc may_40_inc.bat: Runs the 40-year EPS scenario with recharge increased 10 percent \may_50_year may_50.bat: Runs the 50-year EPS scenario with normal recharge \may_50_year_dec may_50_dec.bat: Runs the 50-year EPS scenario with recharge reduced 10 percent \may_50_year_inc may_50_inc.bat: Runs the 50-year EPS scenario with recharge increased 10 percent ---------------------------------------------------------------------------------------------------------------- RESULTS FOR THE EPS SCENARIOS The results for the EPS scenarios are saved in csv files located in the output locations for each scenario. Each file will contain several lines of data, with an example from the 20-year EPS shown below. there are 27,622.0 active cells pump rate (CMD) pump rate aft/ac/yr dry cell ratio dry cell count -33 0.54 0.46 12,785 -34 0.56 0.48 13,140 -35 0.58 0.49 13,485 -36 0.59 0.5 13,837 The first row states the total number of active model cells. This will be the same for each scenario. The second row lists the column headings. These are, in order of left to right, the model pumping rate in cubic meters per day, the EPS in acre-feet per acre per year (ac-ft/ac/yr), the fraction of dry cells (e.g. 0.5 is 50 percent), and the total dry cell count. The remaining rows list the outputs from the EPS scenarios. The EPS is obtained when 50 percent of cells have a saturated thickness of 1.524 meters (5 feet) or less. Thus, the row with the dry cell ratio at 0.5, or 50 percent, is the EPS for each respective scenario. The EPS value in ac-ft/ac/yr is determined as the pumping rate in CMD*((CMD->AFY)/acres in a model cell)). Thus, for this example, the EPS pump rate in ac-ft/ac/yr would be -36*(((1/0.3048)^3)*(640 acres/5280^2 feet)*365.25 days)/((270 meters/0.3048 meters/5280 feet)^2*640 acres)*-1. The values in table 15 are obtained from two minor changes to the model output csv files. First, the total number of active cells used is 27,502 versus 27,622. This accounts for a small number of cells that may be dry at the start of the simulation. Second, a linear regression was used to obtain the EPS values on table 15; thus, the final EPS values may sometimes be slightly different (up to 0.1 ac-ft/ac/yr). ------------------------------------------------------------------------------------------------------------------