Dynamic Model for Stormwater Treatment Areas
W. Walker & R. Kadlec for U.S. Dept. of the Interior Model Version: 4/12/2002
Enter Input File Name --------> test_data.xls Help
Click Here to View Errors
Retrieve Inflow Data
      Select Output Sheet :
Select Case:
Retrieve Case
Edit Input Values
Run Model
Save Case
Save Output
Delete Case
Sensitivity Analysis
Go to Sheet
Selected Case: E25-NE75 25% Emerg - > 75% NonEmerg
Current Series: STA-2 Input Dates: 01/01/65 thru 12/31/95
Current Case: E25-NE75 Output Dates: 01/01/65 thru 12/31/95
Description: 25% Emerg - > 75% NonEmerg
Click here to view DMSTA web site (documentation etc.) press Ctrl-m to return to menu
DMSTA User Manual:
Input File Structure:
Input data for a given treatment area are contained in a separate excel workbook (e.g., "sta2_data.xls").
The file contains daily time series ("InputSeries" sheet) and up to 50 sets of input parameters defining alternative designs.
Each set of input parameters is stored in a separate worksheet.  The sheet name is used to identify the design case.
The case worksheet also contains a summary of results from previous model runs.
The design case name corresponds to the entry in cell D4 on the DMSTA parameters sheet.
The sheet formats must be identical to those in the sample input file (sta2_data.xls) and the DMSTA parameters sheet.
Input files can be created using the procedures outlined below.
They can also be created by editing the sample input workbook (sta2_data.xls) & saving the file with a different name.
When a new input file is retrieved from the DMSTA menu:
The inflow time series is copied to the DMSTA "input series" sheet.
An index of design cases contained in the input file is presented on the DMSTA menu.
A summary of results for all cases in the input file is copied to the DMSTA cases sheet.
Sample Input Files Provided with Software:
File Description
sta2_data.xls Sample Design Cases for STA-2
sta2_sensit.xls Sensitivity Analyses for STA-2
To retrieve an existing design case from an existing input file:
Enter input file name in cell C4  (e.g."sta2_data.xls"); this file must be stored in same directory as DMSTA.XLS
Click "Retrieve Inflow Series" to retrieve inflow time series data & a  list of design cases from the selected file.
Select a design case from the menu
Click "Retrieve Case"
To add a new design case to an existing file:
Load an existing design case (instructions listed above)
Click "Edit Input Parameters" to reflect new design case
Enter new case name in cell D4 of parameters sheet
Click "Run Model"
Click "Save Case" on Menu Sheet
To create a new case & input file:
Enter file name to store input times series & design cases in Cell C4 of "Menu" sheet (replace "sta2.data.xls")
Enter or paste input time series into "InputSeries" sheet of DMSTA.XLS
Edit input parameter values on "Parameters" sheet
Specify up to 6 treatment cells in series and/or parallel.  The cells must be in downstream order.
Specify vegetation type for each cell in row 19; refer to parameter sets in "Calibrations" sheet or use selection tool
Specify a vegetation type of "none" to enter input parameters directly in rows 42-50
Enter design case label in Cell D4
Click "Run Model"
Return to Menu sheet
Click "Save Case".  This saves the input time series & design case in the specified file.
You can retrieve the input case data later by clicking "Retrieve Case"
Click "Save Outputs" to save copy of input & output tables in a separate file.
Note:
Do not delete rows or columns from any worksheet in DMSTA.XLS
To remove cell entries, select range & use "Edit Clearcontents" (from excel menu) or the "Del" key
DMSTA is a modeling tool with a constrained range of applicability. It has been developed and calibrated to information specific to South Florida. It is intended for use in evaluating Everglades Protection Projects.  It is intended for use by individuals with experience in hydrologic & water quality modeling.  It should not be exercised in any situation without careful examination of all features, assumptions and calibrations, as they relate to a given application and to the supporting research upon which the calibrations are based.  The user assumes all risks associated with using the model for engineering design or other purposes.
Click here to view additional documentation at the DMSTA Web Site
Notes on April 17, 2002 Release:        
Added sensitivity analysis & uncertainty analysis routines.  Click "Sensitivity Analysis" button on menu.  
Preliminary estimates of model & parameter standard errors are built into the calibrated parameter sets.  
P cycling parameter "Ks" replaced by "K" = net settling rate at steady state in K/C* model;  
This parameter change is transparent to user when the default calibration sets are used.  
Added parameter set  "SAV_C4" :  Caibrated to ENRP Cell 4, 1998-1999  (IGNORES DATA FROM OTHER SAV PLATFORMS !)
At this point, the "NEWS" calibration is recommended for submergents because it performs better over a range of platforms.
Attempts to reconcile differences between SAV/NEWS  & SAV_C4 calibration are underway.    
Input files from older versions should work, but each case will have to update the "Case" sheet.  
6/8/2002
Run Model
DMSTA Input Values
  Menu
 
            User Input Cells Are Red                  
                                         
Input Variable Units Value Case Description: Filename: test_data.xls   Description                    
Design Case Name  - E25-NE75 25% Emerg - > 75% NonEmerg       for labeling output & reference to input files                
Starting Date for Simulation  - 01/00/00             0 = use first date in input series                  
Ending Date for Simulation  - 01/00/00             0 = use last date in input series                  
Starting Date for Output  - 01/00/00             0 = use first date in simulation                  
Steps Per Day  - 4   Output Variable   Units Value   integration steps per day (decrease to speed up simulation, increase if mass balance errors are encountered)    
Number of Iterations  - 2   Water Balance Error   % 0.0%   number of passes through the dataset to flush out initial values; 2 normally sufficient to simulate system after growin period;
Output Averaging Interval days 7   Mass Balance Error   % 0.0%   averaging interval for time series output (e.g., 7 = weekly, 30 = monthly, 365 = yearly)        
Reservoir H2O Residence Time days 0   Flow-Wtd Conc - With Bypass ppb 27.8   rule for regulating outflow from reservoir ( 0 = residence time not constrained )          
Max Inflow / Mean Inflow  - 0   Flow-Wtd Conc - Without Bypass ppb 27.8   ratio of maximum to mean daily inflow to treatment area; operating objective for reservoir ( 0 = no peak flow control)  
Max Reservoir Storage hm3 0   Geometric Mean Conc   ppb 15.4   maximum reservoir volume (0 = no constraint)                
Reservoir P Decay Rate 1/yr/ppb 0   95th Percentile Conc   ppb 42.4   second order p removal rate in storage reservoir, nominal value ~ .15 1/yr/ppb; not tested on south florida systems  
Rainfall P Conc ppb 10   Freq Cell Outflow > 10 ppb % 68%   rainfall p concentration, wet deposition only (10 ppb assumed in calibrations)          
Atmospheric P Load (Dry) mg/m2-yr 20   Bypass Load   % 0.0%   dry deposition rate (20 mg/mg-yr assumed in calibrations)              
Cell Number -->   1 2 3 4 5 6                          
Cell Label
- EM75 NE75                                  
Vegetation Type -------> EMERG NEWS           assigns pre-calibrated p cycling parameters based upon vegetation type; move cursor to cell in range D:19-I19 & select parameter set from list boxl; select"none" to enter parameters directly in rows 42-51
Inflow Fraction - 1 0           fraction of basin flow entering cell ( 0 = receives inflow from upstream cell only)          
Downstream Cell Number  - 2 0           cell number receiving discharge from this cell ( 0 = outflow from treatment area); cells must be in downstream order  
Surface Area km2 6.500 19.500           effective treatment area (wet at normal operating level)              
Mean Width of Flow Path km 4.82 4.82           width = surface area / flow path length                
Number of Tanks in Series  - 3 6           number of cstr's in series  (can include a fraction)              
Outflow Control Depth cm 60 60           no outflow below this water level (or constant depth)  Zc            
Outflow Coefficient - Exponent  - 3.5 3.5           for computing outflow from depth (b) q / w = a  Z ^b  for Z >= Zc  typicallly ~  3.5      
Outflow Coefficient - Intercept  - 0.8 0.8           flow/width at water depth of 1 m (a),  = 0  for constant depth    typically  0.4 - 1.2      
Bypass Depth cm 140 140           depth at which bypass begins ( 0 = no limit )                
Maximum Inflow hm3/day 0 0           inflow capacity (triggers bypass)  ( 0 = no limit )              
Maximum Outflow hm3/day 0 0           outflow capacity (0 = no limit)                    
Inflow Seepage Rate (cm/d) / cm 0 0           centimeters per day per centimeter of head, reflects transmissivity of soils          
Inflow Seepage Control Elev cm 0 0           drives inflow seepage rate, depth relative to mean ground surface elev          
Inflow Seepage Conc ppb 20 20           p conc of inflow seepage                  
Outflow Seepage Rate (cm/d) / cm 0.008 0.008           centimeters per day per centimeter of head, reflects transmissivity of soils          
Outflow Seepage Control Elev cm 0 0           drives outflow seepage rate, relative to mean ground surface elev, can be < 0          
Max Outflow Seepage Conc ppb 20 20           maximum p conc of recycled & discharged seepage; seepage conc constrained to <= water column conc    
Seepage Recycle Fraction  - 0.5 0.5           fraction of outflow seepage recycled to inflow of cell              
Seepage Discharge Fraction  - 0 0           fraction of outflow seepage discharged to outlet of cell              
Initial Water Column Conc ppb 30 30           initial water-column p concentration                
Initial P Storage Per Unit Area mg/m2 500 500           initial biomass p storage                  
Initial Water Column Depth cm 50 50           initial water column depth relative to mean ground elevation              
C0 =  WC Conc at  0 g/m2 P Storage ppb 4 12           water col umn concentration @ storage = 0 mg/m2 assigned automatically based upon selected vegetation type  
C1 = WC Conc at 1 g/m2 P storage ppb 22 22           water col umn concentration @ storage =1000 mg/m2 " "            
K = Net Settling Rate at Steady State m/yr 16 129           net settling rate in steady state in K/C* model   " "            
Zx = Depth Scale Factor cm 60 60           water depth at maximum areal uptake    " "            
C0 - Periphyton ppb 0 4           C0 for periphyton       " "            
C1 - Periphyton ppb 0 22           C1  for periphyton       " "            
K  -  Periphyton 1/yr 0.00 23.80           K for periphyton       " "            
Zx - Periphyton cm 0 0           Zx for periphyton       " "            
Sm = Transition Storage Midpoint mg/m2 0 400           midpoint of transition from macrophyte to periphyton " "            
Sb = Transition Storage Bandwidth mg/m2 0 80           bandwith of transition from macrophyte to periphyton " "            
                                         
Output Variables Units 1 2 3 4 5 6 Overall Overall = combined outflow from all cells, excluding bypass            
Execution Time seconds/yr 1.58 4.39         4.39 cumulative run time in seconds per simulated year              
Run Date  - 04/12/02 04/12/02         04/12/02 model run date                    
Starting Date for Simulation - 01/01/65 01/01/65         01/01/65 start simulation                    
Starting Date for Output  - 01/01/65 01/01/65         01/01/65 after startup period                    
Ending Date  - 12/31/95 12/31/95         12/31/95 last date                      
Output Duration days 11322 11322         11322 number of simulated days  (excluding startup period)              
Cell Label   EM75 NE75         Total Outflow                        
Downstream Cell Label   NE75 Outflow          -                        
Surface Area km2 6.500 19.5         26.0                        
Mean Water Load cm/d 12.1 3.9         3.0 mean water load to treatment cell (flow / area)                
Max Water Load cm/d 122.6 41.7         30.6 maximum water load to treatment cell (flow/area)              
Inflow Volume hm3/yr 288.2 281.3         288.2 external inflow volume to treatment cell (excludes bypass)              
Inflow Load kg/yr 28833.0 21407.0         28832.1 external inflow load to treatment cell (excludes bypass)              
Inflow Conc ppb 100.0 76.1         100.0 external inflow concentration                  
Treated Outflow Volume hm3/yr 281.3 261.9         261.9 cell outflow + discharged seepage (excludes bypass)              
Treated Outflow Load kg/yr 21407.9 7269.1         7269.1 cell outflow + discharged seepage (excludes bypass)              
Treated FWM Outflow Conc ppb 76.1 27.8         27.8 cell outflow + discharged seepage (excludes bypass)              
Total Outflow Volume hm3/yr 281.3 261.9         261.9 cell outflow + discharged seepage + bypass                
Total Outflow Load kg/yr 21407.9 7270.1         7270.1 cell outflow + discharged seepage + bypass                
Total FWM Outflow Conc ppb 76.1 27.8         27.8 treatment area outflow + bypass                
Bypass Volume hm3/yr 0.00 0.01         0.01                        
Bypass Load kg/yr 0.00 0.97         0.97                        
Bypass Conc ppb 0.0 74.9         74.9                        
Bypass Load % 0% 0%         0% percent of basin load bypassed around treatment area              
Surface Outflow Load Reduc % 25.8% 66.0%         74.8% between outflow from reservoir and outflow from treatment area; ignores loads to groundwater      
Outflow Geometric Mean - Daily ppb 68.7 16.3         16.3 geometric mean of daily outflow concs on days with flow > 0            
Outflow Geo Mean  - Composites ppb 69.2 15.4         15.4 geometric mean of n-day flow-weighted-means, where n = averaging period for results        
Frequency Outflow Conc > 10 ppb % 100% #DIV/0!         68% percent of composite intervals with flow-weighted-mean outflow conc > 10 ppb, excluded bypass       
95th Percentile Outflow Conc ppb 72.3 #NUM!         42.4 95th percentile of  composite outflow concentrations              
Mean Biomass P Storage mg/m2 3733 648         1419                        
Max Reservoir Storage hm3 0.00 0.00         0.00 maximum volume of storage reservoir                
Reservoir Load Reduction % #N/A #N/A         #N/A phosphorus load reduction within reservoir                
Mean Depth cm 62.2 58.3         59.3                        
Minimum Depth cm 15.4 1.0         4.6                        
Maximum Depth cm 123.7 120.8         121.5                        
Frequency Depth < 5 cm % 0.0% 0.0%         0.0%                        
Flow/Width m2/day 160 149         151.5                        
Max 1 Yr Flow-Wtd Conc ppb #N/A #N/A         42.8 rolling average conc, computed for combined outflows only                  
Max 5 Yr Flow-Wtd Conc ppb #N/A #N/A         31.7 ""                      
Max 1 Yr Geometric Mn Conc ppb #N/A #N/A         26.3 ""                      
Max 5 Yr Geometric Mn Conc ppb #N/A #N/A         17.2 ""                      
Depth Percentile vs. Calib Set % 55% 49%         49% approximate rank relative to calibration dataset, overall = last cell            
Q/W Percentile vs. Calib Set % 46% 69%         69% approximate rank relative to calibration dataset, overall = last cell            
Cin Percentile vs. Calib Set % 82% 83%         82% approximate rank relative to calibration dataset, overall = last cell            
Cout Percentile vs. Calib Set % 72% #N/A         #N/A approximate rank relative to calibration dataset, overall = last cell            
Depth Range Flag  - OK OK         0.0 indicates depth out of calibration range                
Q/W Range Flag  - OK OK         0.0 indicates flow/width out of calibration range                
Cin Range Flag  - OK OK         0.0 indicates inflow conc out of calibration range                
Cout Range Flag  - OK           0.0 indicates outflow conc out of calibration range                
Water Balance Error % 0.00% 0.00%         0.00% increase number of steps per day if this is large (e.g.,  > 1 %)            
Mass Balance Error % -0.01% 0.08%         0.05% increase number of steps per day if this is large  (e.g.,  >1%)