SIR2018-5059 MODEL ARCHIVE Created March 22, 2018 Publication Nishikawa, Tracy, ed., 2018, Santa Barbara and Foothill groundwater basins geohydrology and optimal water resources management - developed using density dependent solute transport and optimization models: U.S. Geological Survey Scientific Investigations Report 2018-5059, 4 chap. (A-D), variously paged, https://doi.org/10.3133/sir20185059 Data Release Nishikawa, T., Paulinski S., Stanko, Z., Cromwell, G., and Boyce, S. 2018, SEAWAT model used to evaluate water management issues in the Santa Barbara and Foothill groundwater basins, California: U.S. Geological Survey data release, http://doi.org/10.5066/F74J0DF5 The highest-level directory structure of the original data release is: SIR2018-5059\ ancillary bin georef model output source SIR2018-5059\ Description: ------------ This archive contains all of model input and output files for the Santa Barbara Flow and Transport Models of the Santa Barbara and Foothill groundwater basins (SMFTM). The model simulations and optimizations were run using modified version of SEAWAT (v. 4.00.05) (Langevin and others, 2007) which is included in this model archive data release. 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/sir2018xxxx) to understand the purpose, construction, and limitations of this model. System requirements: -------------------- The models require at least 500 MB of memory and have been run successfully on 64-bit Windows 7 Enterprise OS. Python 2.7 and the Numpy package need to be installed in order to use the pre-processing and post-processing scripts included in the ancillary directory. In addition, the flopy package is necessary for the optimization scripts. The calibrated model can be post-processed using user-preferred post-processing utilities. Reconstruction of data release from online data release archive files: ---------------------------------------------------------------------- The data release archive files need to be extracted into the original directory structures for the model to run correctly. The model depends on the model, bin, and output folders all being extracted into the same top-level folder. To properly assemble the model archive from archive files: 1) Create the directory .\SIR2018-5059\ 2) Download the data release archive files (.zip files) into this directory 3) Uncompress each zip file placing the contents of the zip files into the following subdirectories: zip file name place contents in -------------------------------------------------------------- ancillary.zip .\ancillary\ bin.zip .\bin\ georef.zip .\georef\ model.zip .\model\ output.zip .\output\ source.zip .\source\ Running the models: ------------------ SBFTM consists of two models, an initial conditions model and a transient model. These models are run sequentially independent of PEST by running the RunModels.bat batch file in the .\model directory. RunModels.bat first runs the initial conditions model, then runs a post-processing program on the initial conditions model output, and lastly runs the transient model. Output for both models will be created in the .\output subdirectory. Running the optimization model is described in .\ancillary\optimization readme.txt file. Files and folders: ------------------ The files and folders contained in this archive are described below: A. .\ 1. modelgeoref.txt: File containing the four corners of the model domain. Model data files are in In Lambert Conformal Conic projection (NAD 83) UTM Zone 11N, in meters. The upper left corner of the model area corresponds to the upper left corner of the model cells at row 1 column 1. 2. readme.txt: Contains this readme file documenting the model archive. B. .\ancillary\optimization: This folder contains the instructions and files for the multi-objective optimization model. Additional details are provided in the readme.txt file within each subfolder. C. .\bin\: Contains all of the executable files used to run the model and PEST. beopest64.exe: PEST executable used for running PEST NoComment.exe: Executable used to turn a commented PEST command file into addendum PEST command file readable by PEST. par2par.exe: PEST executable for pre-processing parameters. savelast.exe: Executable used to extract data from MODFLOW/SEAWAT output files. swt_v4x64.exe: Executable to run SEAWAT. D. .\georef\: Contains a shapefile defining the active and inactive model areas. SIR2018_0000.shp: GIS file containing the active and inactive extent of the SBFTM. E. .\model\: Contains the files to run the SBFTM initial conditions and transient models. 1. .\externalfiles: Data files for the SBFTM initial conditions and transient models. _.txt Recharge data in the model area for specific month and year indicated in the file name. .txt Recharge data in the model area for specific year indicated in the file name. Average_Recharge_Zone_#.txt Recharge multipliers for four zones (# = 1-4) defined based on soil permeability and parcel land-use data. Average_Septic.txt Recharge from septic systems. BotLayer.txt Files defining the bottoms of each model layer . IBOUND.txt Files defining the IBOUND array for each model layer . LHK.txt Files defining the base horizontal hydraulic conductivity of layer in region . LVK.txt Files defining the base vertical hydraulic conductivity of layer in region . LSS.txt Files defining the base specific storage of layer . LIP.txt Immobile porosity for layer . LMP.txt Mobile porosity for layer . R_Zone_.txt Files defining the extent of the five recharge zones. StorageUnitZones.txt File defining the HK/VK zones. TopLayer01.txt File defining the top of model layer 1. V50_-_PC_and_parameter_Updates.rch.R_Mult1638 Defines locations of septic system recharge. 2. .\IC: Model files needed to run the SBFTM initial conditions model. RunModels.bat Runs the SBFTM initial conditions and transient models. SBIC.adv Defines how advective flow is simulated. SBIC.bas Defines the initial conditions and some of the boundary conditions. SBIC.btn Defines the initial conditions, some boundary conditions, and how the model is spatially and temporally broken up into model cells, time steps, and stress periods. SBIC.dis Defines how the model is spatially and temporally broken up into model cells, time steps, and stress periods. SBIC.drn Simulates a head-dependent flux boundary at specified cells. The flux from the drain drops to zero when the head in the cell falls below a threshold. SBIC.dsp Defines how dispersion and molecular diffusion is simulated. SBIC.gcg Setup and solve equations simulating a basic ground-water transport model. SBIC.ghb Head-dependent boundary condition to allow groundwater to flow into or out of the model. SBIC.hfb Simulates a groundwater barrier using a hydraulic conductance between two adjacent cells. SBIC.lpf Calculates the conductance between cell centers. SBIC.mlt Defines multiplier arrays for calculation of model properties. SBIC.nam SBFTM initial conditions name file SBIC.oc Used to output head, drawdown, and budget information for specified time periods. SBIC.pcg Solves the groundwater flow equations. SBIC.rch Simulates a specified flux distributed over the top of the model. SBIC.rct Chemical reaction package. SBIC.riv Simulates a head-dependent flux boundary at specified cells. SBIC.ssm Defines the sources of solute and the solute concentration from each source. SBIC.tob Defines the solute concentration and weight by layer(s), row, column, and time. SBIC.vdf Setup and solve the variable-density ground-water flow equation. SBIC.wel Simulates a specified flux to individual cells. SBIC.zon Defines zones to distinguish different parts of the model. slint.txt Text file with instructions for post-processing initial conditions output. 3. .\TCP: Model files needed to run the SBFTM transient model. SBModel.adv Defines how advective flow is simulated. SBModel.bas Defines the initial conditions and some ofthe boundary conditions. SBModel.btn Defines the initial conditions, some boundary conditions, and how the model is spatially and temporally broken up into model cells, time steps, and stress periods. SBModel.dis Defines how the model is spatially and temporally broken up into model cells, time steps, and stress periods. SBModel.drn Simulates a head-dependent flux boundary at specified cells. The flux from the drain drops to zero when the head in the cell falls below a threshold. SBModel.dsp Defines how dispersion and molecular diffusion are simulated. SBModel.gcg Setup and solve equations simulating a basic ground-water transport model. SBModel.ghb Head-dependent boundary condition to allow groundwater to flow into or out of the model. SBModel.hfb Simulates a groundwater barrier using a hydraulic conductance between two adjacent cells. SBModel.hyd Defines settings for hydmod output. SBModel.lpf Calculates the conductance between cell centers. SBModel.mlt Defines multiplier arrays for calculation of model properties. SBModel.nam SBFTM transient name file SBModel.ob_hob Defines the head observation and weight by layer(s), row, column, and time. SBModel.obs Defines location of observation output. SBModel.oc Used to output head, drawdown, and budget information for specified time periods. SBModel.pcg Solves the groundwater flow equations. SBModel.rch Simulates a specified flux distributed over the top of the model. SBModel.rct Chemical reaction package. SBModel.riv Simulates a head-dependent flux boundary at specified cells. SBModel.ssm Defines the sources of solute and the solute concentration from each source. SBModel.tob Defines the solute concentration and weight by layer(s), row, column, and time. SBModel.vdf Setup and solve the variable-density ground-water flow equation. SBModel.wel Simulates a specified flux to individual cells. SBModel.zon Defines zones to distinguish different parts of the model. F. .\output\: Contains all of the output for the SBFTM initial conditions and transient model runs. Due to the size of the transient cbc budget file (SBModel.cbc) it is not included here. A complete set of output files can be generated by running the model. 1. \output\output_IC: Contains the output for the SBFTM initial conditions model run ConcCl.ocn TOB package chloride concentration observation output file SBIC.bhd MODFLOW binary head output file (initial conditions) SBIC.cbc MODFLOW binary cell-by-cell output file (initial conditions) SBIC.glo SEAWAT global listing file (initial conditions) SBIC.lst MODFLOW listing file (initial conditions) SBIC.ucn Post processing file that contains the last stress period of MT3DMS unformatted concentration file (initial conditions) SBIC_Cl.mto MT3DMS observation point concentration file (initial conditions) SBICModel_Cl._mas MT3DMS mass balance file (initial conditions) SBICModel_Cl.ucn MT3DMS unformatted concentration file (initial conditions) 2. \output\output_TCP: Contains the output for the SBFTM transient model run SBModel.cbc MODFLOW binary head output file (transient) *** not included SBModel.fhd MODFLOW head output file (transient) SBModel.glo SEAWAT global listing file (transient) SBModel.hyd_out MODFLOW hydmod output file (transient) SBModel.lst MODFLOW listing file (transient) SBModel.ocn MT3DMS concentration observation output file (transient) SBModel_Cl._mas MT3DMS mass balance file (transient) SBModel_Cl.mto MT3DMS observation point concentration file (transient) SBModel_Cl.ucn MT3DMS unformatted concentration file (transient) G. .\source\: Contains the source code for USGS programs contained in this archive. 1. .\swt_v4_00_05\: Subdirectory containing all the source code for SEAWAT. Minor modifications were made to two original source code files as follows: 1. gmg1.f: To compile for 64-bit, line 399 and 400 have to be edited such that: there is a 'C' in the first column of line 399 and the 'C' is removed from the first column of line 400 2. param.inc: Line 7 was changed from: PARAMETER (MXPAR=999,MXCLST=5000,MXZON=500,MXMLT=500,MXINST=10000) to: PARAMETER (MXPAR=999,MXCLST=5000,MXZON=500, & MXMLT=5000,MXINST=50000) References: ----------- Doherty J.E., 2010, PEST, Model-independent parameter estimation: user manual, 5th edn. (and addendum to the PEST manual), 336 p. Langevin, C.D., Thorne, D.T., Jr., Dausman, A.M., Sukop, M.C., and Guo, Weixing, 2008 SEAWAT Version 4: A Computer Program for Simulation of Multi-Species Solute and Heat Transport: U.S. Geological Survey Techniques and Methods Book 6, Chapter A22, 39 p.