Estimation of Inputs to Florida Bay

Inputs from Florida Mainland - Flows

July 1998

Project Index

Table of Contents

Maps:

Measured Flows

Direct measurements made by the U.S. Geological Survey provide one resource for estimating Florida Bay inflows.   Station locations are shown on the attached map.  Periods of record are listed in a subsequent section.  Given the diffuse nature of the inflows, it not likely that these sites reflect all freshwater inflows to Florida Bay from the mainland.   Measured flows available for only a portion of the period of interest.  For these reasons, the modeling approaches described below are as the primary basis for estimating inflows from the mainland.   Comparisons of observed and predicted inflows are presented in a subsequent section.

Modeled Flows

The South Florida Water Management Model  predicts water balances for South Florida on a uniform grid (2 x 2 mile cells).  The model produces output at a daily time step for Calendar Years 1965-1995.  Output from the 95Base run of the Everglades Restudy most closely approximates current watershed and water-management conditions.   The grid size and topographic data used in the model may not be not sufficient to characterize the spatial distribution of Bay inflows.  The basic water-budget calculations performed by the model, however, should be reasonably accurate.  The model provides useful estimates of spatially-distributed rainfall and evapotranspiration rates over most of ENP for the 1965-1995 period.  Linkage of Bay inflow estimates to the SFWMM would be desirable for evaluating the effects of future changes in water management.

Modeling efforts are underway at WES and USGS (Swain, 1998) to developed a detailed hydrologic model that will serve as a linkage between the SFWMM and Florida Bay water quality models.  Given that these results are not yet available and given the immediate need for estimates of Bay inflows that can be used in preliminary calibrations of the Bay water quality models, a simple water-balance approach is used below to estimate Bay inflows from various ENP watersheds.  The resulting estimates can be refined or replaced with output from the WES/USGS watershed models, once this output is available. 

One advantage of the water budget approach is that it can be applied to produce inflow estimates for recent years (1996-1997), which are important for calibrating and testing Bay water quality models.  Applications of the SFWMM are currently restricted to the 1965-1995 period and are not likely to be extended soon. 

Water Balance Equations

The following water balance model is used below  to generate outflow time series for each ENP basin:

dV/dt     =     Qin   +    Rainfall   -   Evapotranspiration  -  Qout

where,

V = volume stored in marsh (kac-ft)
Qin = basin inflow (kac-ft/month)
Qout = basin outflows (kac-ft/month)  
t = time (months)

An empirical function of the following form is used to represent the basin stage/discharge relationship:

Qout  / A    =     a Z b

Z   =  V / A

where,

a,b = empirical coefficients
A = marsh area (1000 acres)
Z = mean water depth (feet)

Refinements to the above model would consider variations in marsh surface area as a function of elevation. The equations can be combined and integrated numerically to produce an outflow time series, given inflow, rainfall, and evapotranspiration time series, the basin area, and estimates of the empirical coefficients a and b.  Kadlec and Knight (1996) show that discharge is expected to vary with the cubic power of mean depth in a wetland with constant surface area.   This provides an initial estimate of b = 3 for a wetland with constant area.   Given that the wetland area (portion inundated) is not constant, however, calibrated values of b may differ from 3.   The coefficients a and b are calibrated to simulate mean monthly water depths at a representative marsh monitoringe site in each basin.    It should be noted that the coefficients a and b influence only the timing of the predicted outflows.  The total outflow volume is constrained by the overall water blance (inflows + rainfall - evapotranspiration).

Application to SFWMM Output

The feasibility of predicting basin outflows using this approach is demonstrated by applying the model to the monthly water budget data produced by the South Florida Management Model for the 824-mi2 ENP Basin.  Approximate ENP basin boundaries defined in SFWMM are shown below:

map_sfwm.gif (22782 bytes)

Structure inflows (blue squares) include the S12's, S333, and S174 (combined discharge through S332 & S175).  Monthly total inflows and total outflows for the 1965-1995 simulation period are plotted in the following figure:

enpb_q.gif (15873 bytes)

Total inflows primarily reflect structure inflows.  Total outflows primarily reflect overland flows.  Both terms are influenced to a minor extent by groundwater flows and levee seepage.  A strong correlation between inflows and outflows is evident.  Differences between inflows and outflows primarily reflect rainfall, evapotranspiration, and change in storage.

Average water depths have been calculated from the change-in-storage term provided in SFWMM output; depths assume a constant basin area of 824 mi2 (206 model grid cells).  Coefficients of the depth/discharge relationship (a,b) have been calibrated to match the simulated mean depth time series:

enpb_z.gif (14439 bytes)

The following figure compares "observed" and predicted monthly total outflows from the SFWMM ENP basin:

enpb_qo.gif (13836 bytes)

The "observed" values are those generated by the SFWMM from water-balance calculations on the 2 x 2 mile  grid.   The "predicted" flows are derived from the calibrated basin-wide water balance model driven by inflow, rainfall, and evapotranspiration time series.   These estimates cannot be translated directly to Bay inflows, since a portion of the outflows from this basin are towards the west and the basin does not include all of the ENP watershed.   A more detailed cell-by-cell polling of SFWMM output would be required to extract Bay inflows.  Results demonstrate that the basin-wide water-balance model produces outflow estimates which are similar to those produced by the SFWMM.

The long-term record represented by the SFWMM provides important perspectives on temporal variations in flow driven by climatologic factors.  The following figure shows 12-month rolling totals of Rainfall, ET, Inflows, and Outflows for the SFWMM ENP Basin over the 1965-1995 period, expressed as unit area flows.  The apparent correlation between rainfall and ET reflects the dependence of ET on soil moisture levels as water levels drop.  Over this period, annual outflow ranged from 10 to 120 inches/yr.

enpb_12.gif (14794 bytes)

The following figure plots net runoff from ENP watersheds (Rainfall -ET).  Over the 1965-1995 period, net runoff averaged 7.7 inches/yr and ranged from -7.8 to +36 inches/yr.  Compilation of measured rainfall and ET values will be necessary in order to extend the period of record beyond 1995.   Consideration of spatial variations in these values will be required to develop more accurate estimates of net runoff from each watershed.

enpb_ret.gif (11136 bytes)

These results demonstrate that inflows to Florida Bay are strongly dependent on water deliveries to ENP from the WCA's.    Although local contributions from the ENP watershed (Rainfall - ET) are negative in ~50% of the 12-month periods (range, -7. to +36 inches/yr), basin total outflows are positive in every 12-month period (range 10 to 120 inches/yr). 

Applications to ENP Watersheds

Using the above methodology, water-balance calculations have been performed for the ENP watersheds defined in the following figure:

map_wat.gif (22782 bytes)

Basin boundaries have been estimated from GIS layers representing land contours and hydrography.  Boundaries are approximate at best.  Monitoring station categories shown on the map include:

  • Coastal - Coastal/Mangrove fringe stations with water quality measurements
  • Marsh - ENP marsh monitoring points with stage and water quality measurements
  • Inflow - Points of Discharge into ENP (structures, pump stations, or slough monitoring stations)
  • Outflow - USGS flow and water quality monitoring stations at basin outflow points to Florida Bay or Southwest Coast
  • Structure - other regional structures with monitored flows and water quality

Click here to view the watershed map overlayed on landsat image.

Click here to view the watershed map and Florida Bay water quality model domain.

Development of a monthly rainfall time series (October 1985 - December 1997) for each watershed based upon data from the ENP hydrologic monitoring network is described in a separate section of this report.  To provide a slightly longer period for testing the water budget models, the period of record has been extended back to October 1983 using average rainfall values derived from SFWMM input data files for grid cells in each basin.   Click here to see SFWMM grid overlayed on ENP watershed map.

For October 1985-December 1995, evapotranspiration (ET) rates have been extracted from SFWMM output (95Base run) and averaged over grid cells in each basin.   Average monthly ET rates computed from 1986-1995 output have been used for 1996-1997.  Considering the importance of 1996-1997 results to the water quality   modeling effort, refinements to the ET estimates for this period might be derived from  recent measurements by the USGS, once results become available.

The SFWMM grid does not cover the entire area of interest.    Rainfall and ET values for excluded areas have been estimated from neighboring grid cells.

The equations do not directly account for inter-basin flows or levee seepage which is known to occur along the eastern boundary of the Park (into L31N, L31W, and C111 canals,   VanLent et al, 1993; ENP, 1994;  Walker, 1997).  Although S332 and S175 are generally considered inflow points to Taylor Slough, substantial portions of these inflows seep back into the L31N and/or C111 canals because of strong west-to-east hydraulic gradients and permeable soils.   Appreciable portions of flows passing through S18C, S332, and S175 consist of recycled seepage from ENP.

As indicated on the watershed map, northern portions of Taylor Slough and the Rocky Glades are excluded from the water budget calculations.  These areas are assumed to drain primarily to the East into L31N and 31W canals (VanLent et al, 1993; ENP, 1994).   The net inflows into southern Taylor Slough from northern Taylor Slough, the Rocky Glades, and discharges from S332 are assumed to be reflected in the measured flows at the Taylor Slough Bridge.   Jones & Boyer (1997) used the Taylor Slough bridge flows exclusively to represent inflows to southern Taylor Slough.  Including the S175 flows seems more consistent with the topography and provides a slightly better fit of the P37 stage hydrograph (see below), although this assumption needs further examination.  Results of more detailed hydrologic modeling may provide better resolution of the flows in this area.

The southern portion of the Rocky Glades basin does not have inflows from water management structures or direct outflows to Florida Bay.  Land slope indicates that this basin discharges primarily to Shark Slough.  Accordingly, this basin is considered part of Shark Slough for the purposes of generating Bay inflow estimates.   

Seepage losses from Northeast Shark Slough into L31N are not reflected in current calculations; these were reduced after 1991 with operation of structure G211, which raised canal stage in the northern portion of L31N by approximately 1 foot.  

The estimates do not account for potential interactions among the basins.  Under certain conditions, for example, Shark Slough has overflowed into Taylor Slough through the Rocky Glades area  (Van Lent et al., 1993).  These interactions may be important impacts on the spatial distributions of outflows to Florida Bay and may therefore need to be considered in Phase II of this project.

The "Barnes" basin drains into Barnes Sound which is part of the water quality model domain, but not part of Florida Bay.  Outflows from S197 are assumed to occur without marsh interaction.  This represents the only direct structure outflow into the model domain.  Discharges through this structure are infrequent.

Measurements of inflow through S343A and S343B (Broad River Basin) are available starting in January 1995.  Prior to this date, inflows have been estimated using a correlation with Shark Slough inflows (S12ABCD+S333).

The following table gives monitoring stations, preliminary calibration results and links to monthly time series plots for each basin:

Area

Model Coefficients

Basin

mi2

a

b

Inflows

Marsh Stations / Comments

SFWMM ENP

824.0

0.50

2.0

S12s+S333+S174

Depths Predicted by SFWMM

West

101.2

0.12

2.0

None Use Shark Coef.

Broad

232.7

0.12

2.0

S12A + S343AB

P34

Shark

462.5

0.13

2.0

S12BCD+S333

P35+P33

Rocky Glades

99.9

-

-

None

Assumed Part of Shark Slough

Taylor

300.6

0.19

2.5

Taylor  S Bridge + S175

P37

Cape Sable

206.5

0.12

2.0

None (Use Taylor Coef)

Use Shark Coef.

Coastal

70.4

0.47

2.0

S18C - S197

EP

Barnes

23.2

0.47

2.0

None Use Coastal Coef

Total

1497.0

Summary of Predicted Flows

The following table summarizes average water-balance terms for each basin, expressed as total volumes:

Average Annual Flows (kac-ft/yr)

Inflows

Rainfall  

      ET 

Rain-ET

Outflows

SFWMM ENP Basin, Calendar Years 1986-1995

SFWMM ENP (a)

1515.8

2310.5

2015.8

294.7

1779.9

SFWMM ENP (b)

1124.4

2310.5

2015.8

294.7

1589.2

ENP Watersheds, Calendar Years 1986-1997

West

0.0

280.6

275.1

5.5

43.5

Broad

206.2

649.1

596.1

53.0

262.3

Shark & Rocky Glades

731.7

1576.4

1421.7

154.8

843.7

Taylor

99.3

833.1

764.5

68.6

200.7

Cape Sable

0.0

540.8

561.2

-20.4

66.8

Coastal

143.8

183.7

203.1

-19.4

140.3

Barnes

0.0

60.3

68.4

-8.2

8.4

Total

1181.2

4124.0

3890.1

233.9

1565.7

Structure S197 (c)

39.9

39.9

Total

1220.9

4124.0

3890.1

233.9

1605.5

a - Total Inflows & Outflows (including levee seepage, overland inflows, groundwater inflows, etc.)
b - Structure Inflows, Overland + Groundwater Outflows
c - Direct Inflow to Barnes Sound from Structure S197

The following table summarizes average water-balance terms for each basin, expressed as unit area flows:

Unit Area Flows (inches/yr)

Inflows

Rainfal l

    ET  

   Rain-ET

Outflows

SFWMM   ENP Basin, Calendar Years 1986-1995

SFWMM ENP (a)

34.49

52.57

45.87

6.70

40.50

SFWMM ENP (b)

25.59

52.57

45.87

6.70

36.16

ENP Watersheds, Calendary Years 1986-1997

West

0.00

51.98

50.96

1.02

8.05

Broad

16.62

52.30

48.03

4.27

21.13

Shark & Rocky Glades

24.39

52.56

47.40

5.16

28.13

Taylor

6.19

51.96

47.68

4.28

12.52

Cape Sable

0.00

49.10

50.96

-1.86

6.07

Coastal

38.31

48.93

54.10

-5.17

37.39

Barnes

0.00

48.77

55.37

-6.60

6.76

Total

14.79

51.65

48.72

2.93

19.61

a - Total Inflows & Outflows (including levee seepage, overland inflows, groundwater inflows, etc.)
b - Structure Inflows, Overland + Groundwater Outflows

Rainfall and evapotranspiration rates dominate the water balances of each basin. Because rainfall and evapotranspiration rates are similar, basin outflows are similar to structure inflows.   Net runoff (Rainfall - ET) from areas within the park averaged ~3 inches/year during this period. 

Uncertainty in the computed outflows is probably large, especially on short time steps.  If the measured Inflow, Rainfall, and Evapotranspiration (ET) terms for the total basin are each assumed to have a relative standard error of 0.05 (standard error / mean), the computed relative standard error of the outflow term would be 0.18 (calculated by adding variances of individual terms).   Uncertainty is amplified because the outflow computation is dependent upon the difference of two relatively large terms (rainfall - ET). 

Spatial Distribution of Bay Inflows

Comparisons with flows measured by the USGS may provide a basis for refining the above flow estimates and for spatially allocating the outflows from each basin to inflow points represented in the Florida Bay Water Quality Model.   The following inflow points are currently included in the model, as shown on the attached map:

  • Broad River
  • Harney River
  • Shark River
  • Rodgers River
  • Taylor River
  • Mud+Trout Creeks

Addition of new inflow points may be needed to reflect the entire inflow distribution.  Allocation of the predicted basin outflows to specific model input puts is left to the water quality modeling team.

Additional Tasks

The following additional tasks may provide improved estimates of Bay inflows from the Florida Mainland:

  • Development of refined estimates of ET rates for 1996-1997.

  • Refinements to the watershed delineations, as suggested by reviewers or additional hydrologic information.

  • Use of more detailed hydrologic models in place of simple water budget computations (Swain, 1998)

  • Further investigation of methods for representing the effects of seepage losses to the L-31N and C-111 canals along the eastern boundary of the Park

  • Further comparisons and reconciliations of predicted outflows with values measured by the USGS.

Provide Comments

Table of Contents

Project Index - Estimation of Inputs to Florida Bay

 

Broad River Basin

broad_z.gif (9558 bytes)

broad_q.gif (11377 bytes)

Shark River Slough Basin

shark_z.gif (10186 bytes)

shark_q.gif (11456 bytes)

Taylor Slough Basin

taylor_z.gif (11567 bytes)

taylor_q.gif (10529 bytes)

Coastal Basin

coast_z.gif (12248 bytes)

coast_q.gif (14319 bytes)

Cape Sable / Whitewater Bay

sable_z.gif (9805 bytes)

sable_q.gif (10187 bytes)

Hydraulic coefficients (a,b) for the Cape Sable/ Whitewater Bay region are assumed to equal those calibrated for the Broad River Basin.

Project Index

http://www.wwwalker.net/flabay/enpflows.htm      Updated:  03/30/02