Data Release: FTLOADDS (combined SWIFT2D surface-water model and SEAWAT groundwater model) simulator used to assess proposed sea-level rise response and water-resource management plans for the hydrologic system of the South Florida Peninsula for the Biscayne and Southern Everglades Coastal Transport (BISECT) Model http://dx.doi.org/10.5066/P9MDUQPK Publication: Swain, E.D., Lohmann, M.A., and Goodwin, C. R., 2019, The Hydrologic System of the South Florida Peninsula: Development and Application of the Biscayne and Southern Everglades Coastal Transport (BISECT) Model : U.S. Geological Survey Scientific Investigations Report 2017-5045 (http://dx.doi.org/10.3133/sir20195045) ------------------------------------------------------------------------------ 11/20/2017 Model archive created ------------------------------------------------------------------------------ BACKGROUND INFORMATION AND PREVIOUS PUBLICATIONS DOCUMENTING SOME OF THE MODEL DEVELOPMENT AND CALIBRATION USED FOR THESE SCENARIOS: This application utilizes the numerical application FTLOADDS, with some advancements, originally described and developed in Langevin, C.D., Swain, E.D., and Wolfert, M.A. 2005, Simulation of integrated surface-water/ ground-water flow and salinity for a coastal wetland and adjacent estuary: Journal of Hydrology 314, 212-234. The numerical model domain BISECT is created by combining two models and their domains, TIME and Biscayne, which were developed in: Wang, J, D., Swain, E.D., Wolfert, M. A., Langevin, C.D., James, D. E., and Telis, P. A., 2007 - Applications of Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) to Simulate Flow, Salinity, and Surface-Water Stage in the Southern Everglades, Florida: U.S. Geological Survey Scientific Investigations Report 2007-5010. Lohmann, Melinda A., Swain, Eric D., Wang, John D., and Dixon, Joann, 2012. Evaluation of Effects of Changes in Canal Management and Precipitation Patterns on Salinity in Biscayne Bay, Florida, Using an Integrated Surface-Water/Groundwater Model: U.S. Geological Survey Scientific Investigations Report 2012-5099. ------------------------------------------------------------------------------ Several modifications and improvements to the model representation and the FTLOADDS code are developed for this report, including: 1. The extension of the western model boundary farther offshore, 2. The vertical aquifer discretization was changed from 10 layers in TIME and 20 layers in Biscayne to 15 layers in BISECT, 3. Spatially variable albedo was included in the computation of the solar radiation term. Both the TIME and Biscayne models underwent limited calibration in their original developments. After combining into BISECT, further calibration was undertaken to adjust coastal streamflow through modifications of frictional resistance and land elevations. Groundwater evapotranspiration was reduced from initially assumed values to account for coastal vegetation which reduces its transpiration rate by additional stomatal resistance. The underlying directories contain all of the input and output files for the FTLOADDS simulations describe in the report. 22 simulations were completed during this study using BISECT. Descriptions of the data in each subdirectory are given to facilitate understanding of this model archive. Files descriptions are provided for select files to provide additional information that may be of use for understanding this model archive. 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 complete model report (https://doi.org/10.3133/SIR20195045) to understand the purpose, construction, and limitations of this model. System requirements: -------------------- The models contained in this data release were run using the run.bat batch files in the appropriate subdirectories of model\ which access bin\FTLOADDS.exe in this data release. The executable was compiled to be used on a personal computer with various Windows operating systems. For this study, simulations were run on a personal computer with 16MB Random Access Memory (RAM) and 64-bit Windows 7 OS (Service Pack 1). Reconstructing the data release from the online data release: ------------------------------------------------------------- This data release is available from: http://dx.doi.org/10.5066/P9MDUQPK The models, along with post-processing tools, 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 directory level. For example, the files in the "georef" subdirectory are zipped into a zip file named "georef.zip". For the “output” subdirectory there are zip files for each model run. For example, for the calibration run there is the “Output.Calibration.zip” which contains the individual output files. All zip files should be unzipped into a directory with the same name as the zip file name without the .zip extension. These unzipped directories are all at the same top level as the "georef" subdirectory created from "georef.zip" except for output files, which should be placed under the folder called "output" with subdirectory folders for each run, such as a folder for "Calibration". Running the model(s): --------------------- In each of the subdirectories of model\ the model is executed with the batch file run.bat, which runs the executable bin\FTLOADDS.exe with the input file ftlrunfiles.ftl. Double click on the run.bat file to initiate the run. The model\ExternalFiles subdirectory contains all the shared input files for all of the model runs. The calibrated model and the other subdirectories under model\ contain only the input files that change for the particular model run. Model output are found in the corresponding subdirectories of \output. After the model run is completed the additional model code provided in the ancillary\ directory, GWFlow.exe, should be run in order to display the flux data from the timeBB.ufl file in ASCII format. Details on GWFlow.exe output are described below. For all these simulations, the model normally terminates at 3288 days but displays an end-of-file error message due to an undefined stop time. For scenarios including CERP management, the model terminates at 1827 days. The highest-level directory structure of the original data release is: ancillary\ bin\ georef\ model\ output\ source\ Files: readme.txt: This file documents the structure of the model archive, defines the model input and output files, and provides instructions on how to run the simulations. modelgeoref.txt - ASCII file with the four corners of the model domain in decimal degrees.. Ancillary\ The directory Ancillary\ contains the raw input data collected for the input model files and three post processing codes to allow for easier viewing of the model output. InputData The input for the BISECT model was a combination of data collected for two prior models (Biscayne Bay model and TIME model) that were combined to create the new model. There are two folders within this directory that contain this data \Biscayne and \Time. Additionally, the reports for both of these models are included in pdf form. PostProcessingCodes FlowVectors – Extracts groundwater fluxes from timeBB.ufl and outputs timeseries and temporal average grids for recharge, storage, flux to constant head, aquifer flow in all three directions, well flux, drain flux, ET flux, and head-dependent boundary flux. Should be executed in the associated output\ subdirectory for the model run. The executable for this code is GWFlow.exe Newformat – Extracts surface-water output from FTLBIN1.DAT file. Outputs grids to stage.dat, depth.dat, uvel.dat, vvel.dat, Salinity.dat, Temperature.dat, and files names after individual field site locations and data type. Should be accompanied by the date.csv file and executed in the associated output\ subdirectory for the model run. Salextractor – Extracts groundwater salinity from MT3D001.UCN file and outputs grids to file salinity.dat. Should be executed in the associated output\ subdirectory for the model run. bin\ The directory bin\ has the compiled 64-bit executable FTLOADDS model software for PC: The executable for FTLOADDS Version 2 is FTLOADDS.exe. This is accessed from the run.bat batch file in each of the model\ subdirectories. Georef\ The directory Georef\ contains a shapefile defining the active and inactive areas of the model documented in the report. model\ The zipped directories model\ contains 23 subdirectories with shared input and input for the different sensitivity-run and scenario models: Calibration – Input and run files for calibrated model ExternalFiles- Input files used for all model simulations Scenario Runs: Scenario.BarrLay5 - Input and run files for simulation with aquifer barrier to one-half depth. Scenario.BarrLay10 - Input and run files for simulation with aquifer barrier to full depth. Scenario.CERP - Input and run files for simulation with Comprehensive Everglades Restoration Plan (CERP) inflows Scenario.CEP60SLB10 - - Input and run files for simulation with sea-level rise of 60 cm, aquifer barrier to full depth, and CERP inflows. Scenario.SLR30 - Input and run files for simulation with sea-level rise of 30 cm. Scenario.SLR60 - Input and run files for simulation with sea-level rise of 60 cm. Sensitivity Runs: Sens.ConstAlbedo - Input and run files for simulation with constant albedo of 16.9 percent over model domain. Sens.HighCanalCond - Input and run files for simulation with canal conductance 10 times higher. Sens.HighHorizCond - Input and run files for simulation with aquifer hydraulic conductivity doubled. Sens.HighManning - Input and run files for simulation with surface-water frictional resistance 10 percent higher. Sens.HighTidalAmp - Input and run files for simulation with the amplitudes of the tidal frequencies raised 10 percent. Sens.HighTranLayTh - Input and run files for simulation with the transition-layer thickness 0.05 meters. Sens.HighVertCond - Input and run files for simulation with vertical hydraulic conductivity doubled. Sens.HighWind - Input and run files for simulation with the wind frictional coefficient doubled. Sens.LowCanalCond - Input and run files for simulation with canal conductance lowered to 1/10 the values. Sens.LowHorizCond - Input and run files for simulation with aquifer hydraulic conductivity reduced to one half. Sens.LowManning - Input and run files for simulation with surface-water frictional resistance 10 percent lower. Sens.LowTidalAmp - Input and run files for simulation with the amplitudes of the tidal frequencies lowered 10 percent. Sens.LowTranLayTh - Input and run files for simulation with the transition-layer thickness 0.0001 meters. Sens.LowVertCond - Input and run files for simulation with vertical hydraulic conductivity reduced to one half. Sens.LowWind - Input and run files for simulation with the wind frictional coefficient halved. The directory Output\ needs to be created to contain 22 subdirectories with output from the different sensitivity-run and scenario models using the 22 zipped files (each model simulation is named based on the run name, with one folder for each simulation, below lists the correct naming convention for this folders.) Output.Calibration – Model output for calibrated simulation. Scenario Runs: Output.Scenario.BarrLay5 - Model output for simulation with aquifer barrier to full depth. Output.Scenario.BarrLay10 - Model output for simulation with aquifer barrier to half depth. Output.Scenario.CERP - Model output for simulation with CERP inflows. Output.Scenario.CP60SLB10 - Model output for simulation with sea-level rise of 60 cm, aquifer barrier to full depth, and CERP inflows. Output.Scenario.SLR30 - Model output for simulation with sea-level rise of 30 cm. Output.Scenario.SLR60 - Model output for simulation with sea-level rise of 60 cm. Sensitivity Runs: Output.Sens.ConstAlbedo - Model output for simulation with constant albedo of 16.9 percent over model domain. Output.Sens.HighCanalCond - Model output for simulation with canal conductance 10 times higher. Output.Sens.HighHorizCond - Model output for simulation with aquifer hydraulic conductivity doubled. Output.Sens.HighManning - Model output for simulation with surface-water frictional resistance 10 percent higher. Output.Sens.HighTidalAmp - Model output for simulation with the amplitudes of the tidal frequencies raised 10 percent. Output.Sens.HighTranLayTh - Model output for simulation with the transition-layer thickness 0.05 meters. Output.Sens.HighVertCond - Model output for simulation with vertical hydraulic conductivity doubled. Output.Sens.HighWind - Model output for simulation with the wind frictional coefficient doubled. Output.Sens.LowCanalCond - Model output for simulation with canal conductance lowered to 1/10 the value. Output.Sens.LowHorizCond - Model output for simulation with aquifer hydraulic conductivity halved. Output.Sens.LowManning - Model output for simulation with surface-water frictional resistance 10 percent lower. Output.Sens.LowTidalAmp - Model output for simulation with the amplitudes of the tidal frequencies lowered 10 percent. Output.Sens.LowTranLayTh - Model output for simulation with the transition-layer thickness 0.0001 meters. Output.Sens.LowVertCond - Model output for simulation with vertical hydraulic conductivity halved. Output.Sens.LowWind - Model output for simulation with the wind frictional coefficient halved. The model outputs the file timeBB.ufl. This file is close to 84 GB and is restrictive to downloads. Therefore, this file has been processed out using the GWFlow.exe code and then removed from the output folders. Fifteen files that contain the data from this file are created: GWconstH.dat - Constant head cell fluxes GWET.dat - Groundwater evapotranspiration flux GWFFF.dat - Aquifer flux front cell face GWFRF.dat - Aquifer flux right cell face GWrch.dat - Groundwater recharge flux GWRivers.dat - River-aquifer flux GWStor.dat - Groundwater storage flux aveconstH.dat - Simulation averaged constant head cell fluxes aveDrains.dat - Simulation averaged river-aquifer flux aveET.dat - Simulation averaged groundwater evapotranspiration flux aveFFF.dat - Simulation averaged aquifer flux front cell face aveFRF.dat - Simulation averaged aquifer flux right cell face aveGWrch.dat - Simulation averaged groundwater recharge flux aveGWStor.dat - Simulation averaged groundwater storage flux source\ The directory source\ contains the model code subdirectories for FTLOADDS version 2.0 FTLOADDS\ Code and executable for FTLOADDS with subdirectories swift\ containing the SWIFT2D code, Updates\ containing code updates for version 2.0, and Release\ where the compiled object files and executable are written. -------------------------------------------------------------