Publication: Knight, J.E., 2020, Simulation of Groundwater-level Changes from Projected Groundwater Withdrawals in the Truxton Basin, Northwestern Arizona: U.S. Geological Survey Scientific Investigations Report 2020–5017E, 52 p., https://doi.org/10.3133/sir20205017E. Data release: Knight, J.E., 2020, MODFLOW-NWT groundwater model used for simulating potential future pumping scenarios and forecasting associated groundwater- level changes in the Truxton aquifer on the Hualapai Reservation and adjacent areas, Mohave County, Arizona: U.S. Geological Survey data release https://doi.org/10.5066/P9O2WGLS -------------------------------------------------------------------------------- SIR2020-5017E/ Description: ----------- The underlying directories contain all of the input and output files for the simulations described in the report, MODFLOW-NWT (v 1.1.4) source code, and PEST++ (v 4.2.5) source code. Descriptions of the data in each subdirectory are given to facilitate understanding of this data release. Files descriptions are provided for select files to provide additional information that may be of use for understanding this data release. See MODFLOW-NWT or PEST++ documentation in the appropriate sub-folder of the source directory for information on file types not explicitly described herein. 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/sir20205017E) to understand the purpose, construction, and limitations of this model. Reconstructing the data release from the online data release: ------------------------------------------------------------- This data release is available from: https://doi.org/10.5066/P9O2WGLS The models will run successfully only if the original directory structure is correctly restored. The data release is broken into several pieces to reduce the likelihood of download timeouts. Small files (readme.txt and modelgeoref.txt) are available as uncompressed files. All other files are zipped at the subdirectory level. For example, the files in the "georef" subdirectory are zipped into a zip file named "georef.zip". All zip files should be unzipped into a directory with the same name as the zip file name without the .zip extension. The highest-level directory structure of the original data release is: SIR2020-5017E/ ancillary/ bin/ georef/ model/ output/ source/ The full directory structure of the data release and the zip files associated with each subdirectory are listed below. System requirements: -------------------- The models contained in the model subdirectory of this data release were run using the MODFLOW-NWT_64.exe executable in the bin directory of this data release. These models have been run successfully on computers running the following Windows operating systems (OS): o 64-bit Windows 7 OS o 64-bit Windows 10 OS For Monte Carlo forecasts, the model input files of the transient scenario models were adapted to run in parallel on the linux-based USGS "Yeti" high performance computer. Each ensemble of 10,000 models took approximately 24 hours to complete using 240 parallel processing nodes. External file pointers were modified to include a leading "./" and the model files were run by executing the MODFLOW-NWT executable compiled for linux located in the subdirectory ancillary\tbhm_mc_runs\mc_bin Optional system requirements: ----------------------------- Python 3 and the following packages need to be installed in order to use the run-model and post-processing scripts included in the tbhm_mc_runs\ and tbhm_ss_calibration\ subdirectories in the ancillary directory: Numpy (v 1.13.1) Flopy (v 3.2.11) Pyemu (v 0.6) Previous versions of these codes may possibly work, later versions will probably work. Use of these scripts are only necessary to reproduce the steady-state calibration results or Monte Carlo forecast results. The calibrated steady-state model and 9 transient scenario models located in the model\ directory can be run and post-processed using standard MODFLOW or user- preferred post-processing utilities. Running the model(s): --------------------- Each simulation has a name, listed below, which has an associated name file used to run the corresponding simulation. When MODFLOW-NWT is executed, a single prompt is issued for the name file that specifies the input files to be used in the simulation, as well as, the output files that are created. Responding to this prompt with the name file listed below will initiate the corresponding simulation. The simulations that are documented can be run by opening a command line/terminal in the model subdirectory and typing: ..\..\bin\MODFLOW-NWT_64.exe .nam For example, the calibration/verification simulation can be run by opening a terminal in the model/tbhm_ss/ subdirectory and typing: ..\..\bin\MODFLOW-NWT_64.exe tbhm_ss.nam Summary of MODFLOW-NWT simulations for the area-wide model ---------------------------------------------------------------------- MODEL SIMULATION SUBDIRECTORY DESCRIPTION OF AND NAME FILE NAME SIMULATION ----------------------------- --------------------------------- tbhs_ss Steady-State, primary bedrock and no pumping tbhm_tr_Xhigh_Qhigh Transient, high bedrock and high pumping tbhm_tr_Xhigh_Qlow Transient, high bedrock and low pumping tbhm_tr_Xhigh_Qmedium Transient, high bedrock and medium pumping tbhm_tr_Xlow_Qhigh Transient, low bedrock and high pumping tbhm_tr_Xlow_Qlow Transient, low bedrock and low pumping tbhm_tr_Xlow_Qmedium Transient, low bedrock and medium pumping tbhm_tr_Xprimary_Qhigh Transient, primary bedrock and high pumping tbhm_tr_Xprimary_Qlow Transient, primary bedrock and low pumping tbhm_tr_Xprimary_Qmedium Transient, primary bedrock and medium pumping NOTE1: The outputs from each are created in the appropriate output directory. For example, tbhm_ss model results are saved to the output/output.tbhm_ss directory. For the model to run successfully as written, you must either download and extract the output zip file, or create an empty output folder directory of the same name and structure. Note that running the model will write output files to the associated output folder and will overwrite any existing output files of the same name. NOTE2: The models must be run using the MODFLOW-NWT_64.exe executable in the bin directory created using the bin.zip file in this data release. Files: ----- readme.txt : This file documents the structure, directories and files, and instructions on how to run the model for this data release. modelgeoref.txt : ASCII file with the four corners of the model domain in decimal degrees. ancillary/ Description: ----------- This directory contains ancillary data including a directory with python scripts for post-processing model input and output along with resulting post-processed results. tbhm_ss_calibration Description: ----------- This directory contains the model input files, template and instruction files, and pest control file for PEST++ GLM calibration of the steady-state model documented in the report. tbhm_mc_runs/ Description: ----------- This directory contains the model input files, template and instruction files, and pest control file for running Monte Carlo forecasts on each of the transient scenarios using PEST++ SWP as documented in the report. Output files from Monte Carlo runs are also included. NOTE: Monte Carlo forecasts were run on a linux-based high performance computer. MODFLOW input files in the tbhm_mc_runs/ directory are written to work with linux machines, and may not work on other operating systems. mc_bin/ Description: ----------- This directory contains versions of MODFLOW-NWT and PEST++ SWP compiled for linux operating systems and used to produce the Monte Carlo forecast results. mc_sweep_out Description: ----------- This directory contains Monte Carlo forecast output files. One file in .csv format for each of the 9 transient scenarios. See PEST++ documentation in directory source/pestpp/documentation for more information. Description: ----------- The rest of the subdirectories under tbhm_mc_runs contain the model input files, template and instruction files, and pest control file for running Monte Carlo forecasts for each of the 9 scenario. The corresponds to the name in the table above. Files: ------ .bas BAS package input file .dis DIS package input file .drn DRN package input file .hyd HYDMOD package input file .nam NAM file .nwt NWT package input file .oc OC package input file .rch RCH package input file .upw UPW package input file .wel WEL package input file botm_layer_.ref Bottom elevation of delc.ref Cell lengths in column direction delr.ref Cell lengths in row direction model_top.ref Top elevation of model layer 1 usgs.model.reference Spatial and temporal setting info WEL_.dat WEL package data by stress period DRN_.dat Drain package data by stress period hani_.ref Horizontal anisotropy by layer hk_.ref Horizontal hydraulic conductivity by layer ibound_.ref Active/Inactive flag and aquifer property zone by layer rech_.ref Recharge flux by stress period ss_.ref Specific storage by layer strt_.ref Starting heads by layer sy_.ref Specific yield by layer vani.ref Vertical anisotropy by layer factors.dat interpolation factors used to interpolate estimated properties at pilot points to rest of zone multipliers.dat multiplier parameters selected by pestpp for a particular model run points_.dat XY location and parameter value at pilot point locations ZB_array.dat Zone budget array for processing simulated model fluxes zb_out.dat output file of model fluxes summed by observation zone for.ddns model output of simulated drawdowns at forecast locations hydout.heads formatted output from .hyd.bin pot.heads simulated heads tracked for potential data-worth analysis *.ins pestpp instruction files mc_sweep_in.csv Monte Carlo ensemble parameter values tbhm_mc_sweep.pst pestpp control file for running Monte Carlo with pestpp-swp run_model_tbhm_mc_sweep.py Forward-run python script tbhm_mc.rec pestpp rec file trux_zone.ref Array used for locating potential head observation sites within Truxton Basin ZonePropArray_.ref Lookup array created in forward-run script to store property values in zonal arrays, which are subsequently pulled into model layers based on zone arrays of each layer Zones_.ref Property zone coverages by model layer, for assigning interpolated zone properties tbhm_mc_sweep.rns pestpp rns file tbhm_mc_sweep.rst pestpp rst file *.tpl pestpp template files tbhm_mc_sweep.pfm pestpp pfm file bin/ Description: ----------- This directory contains the executable files used to run the calibrated steady-state model and all transient scenarios documented in the report. These executables were compiled for the 64-bit Windows Operating Systems. Files: ------ MODFLOW-NWT_64.exe: MODFLOW-NWT executable (v 1.1.4) used to run steady-state and transient models pestpp-glm.exe: PEST++ GLM executable (v 4.2.5) used to calibrate steady-state model pestpp-swp.exe: PEST++ SWP executable (v 4.2.5) used to run Monte Carlo ensembles of transient models georef/ Description: ----------- This directory contains a shapefile defining the active and inactive areas of the area-wide model documented in the report. Files: ----- SIR2020_5017E.dbf: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.prj: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.shp: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.shx: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.sbx: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.cpg: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.sbn: part of the SIR2020_5017E.shp shapefile SIR2020_5017E.shp.xml: part of the SIR2020_5017E.shp shapefile model/ Description: ----------- This directory contains all of the input files for the calibration/verification run and the 9 scenario runs documented in the report. doi: https://doi.org/10.3133/sir20205017E externalfiles/ Description: ----------- This directory contains common or shared files used by the 9 transient models documented in the report. Steady-state model does not use these files. Files: ------ DRN_.dat Drain package data by stress period hani_.ref Horizontal anisotropy by layer hk_.ref Horizontal hydraulic conductivity by layer ibound_.ref Active/Inactive flag and aquifer property zone by layer rech_.ref Recharge flux by stress period ss_.ref Specific storage by layer strt_.ref Starting heads by layer sy_.ref Specific yield by layer vani.ref Vertical anisotropy by layer / Description: ----------- The rest of the subdirectories under model directory contain the MODFLOW-NWT model input files for running the calibrated steady- state model and each of the 9 scenario. The corresponds to the name in the table above. Files: ------ .bas BAS package input file .dis DIS package input file .drn DRN package input file .hyd HYDMOD package input file .nam NAM file .nwt NWT package input file .oc OC package input file .rch RCH package input file .upw UPW package input file .wel WEL package input file botm_layer_.ref Bottom elevation of delc.ref Cell lengths in column direction delr.ref Cell lengths in row direction model_top.ref Top elevation of model layer 1 usgs.model.reference Spatial and temporal setting info WEL_.dat WEL package data by stress period output/ Description: ----------- This directory contains all of the output files for the calibrated steady-state model run and the 9 transient scenario model runs documented in the report. output./ Description: ----------- The rest of the subdirectories under output directory contain the MODFLOW-NWT model output files for running the calibrated steady-state model and each of the 9 transient scenarios. The corresponds to the name in the table above. Files: ------ .cbc: Cell-by-cell flow binary output file .hds: Heads binary output file .hyd.bin: HYDMOD binary output file .list: Readable model run log & output file source/ Description: ----------- This directory contains the source files for code versions used to run model simulations included in the model data release. MODFLOW-NWT_1.1.4: Source code of MODFLOW-NWT (v 1.1.4), the software used to run steady-state and transient models pestpp: Source code of the PEST++ software suite, of which two programs are used by this model: PEST++ GLM executable is used to calibrate the steady-state model, and PEST++ SWP is used to run Monte Carlo ensembles of transient models