Combining a Conceptual Hydrologic Model (SWAT) and a Hydrodynamic Model (Telemac 3D) to Simulate Reservoir Dynamics in. Eagle Creek Reservoir, Indiana

Transcription

1 Combining a Conceptual Hydrologic Model (SWAT) and a Hydrodynamic Model (Telemac 3D) to Simulate Reservoir Dynamics in Eagle Creek Reservoir, Indiana 1

2 Purpose of Modeling Effort Fate and Transport of Chemicals from Agriculture, Industry, Residences Support field work on chemical absorbtion in fish Improve protection of reservoir public water supply

3 Road Map I. Eagle Creek Watershed-Reservoir System II. Modeling Framework III. Calibration Data IV. Preliminary Results A. Calibration with Upstream USGS Station B. Calibration with Calculated Water Balance Flow V. Conclusions and Future Work 3

4 I. Eagle Creek Watershed- Reservoir System 104,000 acres Huc 10 in size Northern two-thirds mostly agricultural Southern one-third mostly urban Reservoir primary public water supply for Indianapolis, IN 4

5 II. Modeling Framework Upstream SWAT (Eagle Creek and Fishback) Side streams SWAT (School Branch) Bathymetry data Land use data Daily precipitation Soil cover data Digital Elevation Map Hydrologic Model Inflow 1 (daily discharges)... Inflow N (daily discharges) Hydrodynamic Model Stream flow velocities in 3 directions (North, East, Up) Pool elevation Outflow data (daily discharges) 5

6 SWAT Model Upstream of Reservoir

7 SWAT Model Upstream of Reservoir 7

8 2007 Cropland Data Layer 8

9 60:40 Corn:Soybean Ratio 9

10 SSURGO 10

11 Slope 11

12 Point Sources 10 Permitted Discharges in Watershed 4 Sub-basins Affected Actual Discharge Data Reported to EPA Under CWA Formatted for SWAT Conversion to Metric 12

13 SWAT Model Side Streams Entering Reservoir Reaches and drainage basins 2007 CDL SSURGO soils No USGS monitoring in these drainage basins 13

14 TELEMAC3D Numerical modeling system for free surface hydrodynamics Developed by the Laboratoire National d'hydraulique, a department of Electricité de France's Research and Development Division Solves the 3D-fluid mechanic equations(navier-stokes) and transport-diffusion equation for intrinsic values (temperature, salinity, concentration) High capacity algorithms based on finite-element method Space is discretised in the form of an unstructured grid of triangular elements Represents the physical processes observed in the real world 14

15 Telemac Reservoir Grids 7335 Triangular Grids

16 TELEMAC3D Bathymetry Data Using the depth contour lines measured by the IDNR in 1994 Rating curve Volume-Water level computed by Telemac3D depending on the bathymetry of the reservoir Voulme (m3) 3.5E E E E E E E E+00 Eagle Creek Reservoir Relationship between water level and volume Pool elevation above NGVD 1929 (m) 16

17 III. Calibration Data A.Upstream USGS Station Eagle Creek measured at Zionsville: daily discharges [USGS gage station ] B.Upstream Watershed Outlet Flow from Reservoir Water Balance Water Balance Calculation: Water intake outflow measured by Veolia Water Indianapolis: daily discharges Pool elevation in the reservoir measured: daily water level [USGS gage station ] Eagle Creek measured at Clermont (1 miles downstream the reservoir): daily discharges [USGS gage station ] 17

18 Upstream Watershed Outlet Flow from Reservoir Water Balance 1. Use Bathymetry data and Telemac 3D model design to estimate relationship between pool elevation and reservoir storage volume. 2. Use water balance calculation to estimate Total Reservoir Inflow. 3. Total Reservoir Inflow - Uncalibrated SideStreams Inflow = Upstream Watershed Outlet flow Uncalibrated side streams could introduce errors in accuracy of upstream watershed outlet flow Flow from side streams could be re-estimated by using calibrated parameter values from upstream watershed (Future Work) 18

19 WATER BALANCE CALCULATION What should be the discharge of all the inflows which reach the reservoir based on the overall watershed-reservoir outflow and the water level in the reservoir? 3.5E E+07 Eagle Creek Reservoir Relationship between water level and volume Qinflow -Qoutflow = ΔV/Δt Qinflow = Qoutflow + ΔV/Δt ΔV/Δt = [V(i+1) -V(i-1)] / [2* Δt ] V estimated via the Volume-Water level rating curve computed by TELEMAC3D using the bathymetry data Water level and Ouflow measured by the USGS and Veolia Water Indianapolis Voulme (m3) 2.5E E E E E E Pool elevation above NGVD 1929 (m) 19

20 IV. Preliminary Results A.Uncalibrated Simulation with Observed Flows from Upstream USGS Station 20

21 B. Manual Calibrated Results Altering PET Methods 21

22 C. Manual Calibrated Results Adding Value to Curve Number 22

23 Auto-Calibration a Success 23

24 Future Work Model watershed for water quality Smaller sub-basin size breaking across streams where fish are being studied Use Indiana s high resolution (1.5 m) DEMs 24

25 Team with NRCS to Incorporate Agricultural Land Uses 25

26 Conclusions TBD Field work and modeling work will be used to draw conclusions 26