CAEP Study Site: Cannonsville Reservoir Watershed

Transcription

1 Tammo Steenhuis, Zach Easton, Eric White, C. Shoemaker Joshua Woodbury and many others of the Soil and Water Laboratory, Farming Community and the Watershed Agriculture Program in Walton This project is supported by National Research Initiative Competitive Grant from the USDA Cooperative State Research, Education, and Extension Service

2 CAEP Study Site: Cannonsville Reservoir Watershed

3 Best Management The Catskills (or the New York City Source Practices Watersheds) is were one of the done few watersheds right! where the (dissolved) phosphorus concentration is significantly decreasing after BMP installation Why is that?

4 This is why! In the temperate northeast hydrology is, to a great extent, driven by topology Near stream areas, and areas with shallow, soils over a slowly permeable glacial till or bedrock produce the majority of runoff

5 Biological and Environmental Engineering Runoff location Variable Source Areas concept included in BMP installation from early on Soil & Water Research Group

6 BMP s BMP effectiveness is determined by the three L s Location Location Location

7 The Science of Reducing P Keep the manure out of the runoff water Requires knowledge of saturated areas Reduce P loading Reducing P in feed import To optimize P reductions, models are needed simulating hydrology of the region and are in agreement with farmer s perceptions

8 The Case of SWAT

9 Biological and Environmental Engineering Soil Water Assessment Tool (SWAT) Basin scale hydrological and water quality.6 ha continental Runs on easily available data Process based chemistry vs export coefficient type models Does not require much calibration for hydrology More for chemistry and sediment Soil & Water Research Group

10 Soils Landuse SWAT defines HRUs as the coincidence of soil type and landuse Hydrological/chemical properties are defined at the HRU HRUs So runoff/p loss is the same here (lowland pasture) As here (upland pasture) We know this is not the case

11 hydrology model approach For each hydrologic response unit rainfall ET If rainfall > ET, soil wets up Q P e P 2 e S

12 Biological and Environmental Engineering USDA-NRCS Curve Number Model Runoff =P e2 /(Pe+S) That is a big problem S=254/CN-254 Tables link CN to land use and soil infiltration capacity Infiltration capacity in excess of rainfall intensity! Soil & Water Research Group

13 hydrology model approach Total runoff and the saturated protion of the watershed is a function of amount of precipitation rainfall ET If rainfall > ET, soil wets up If rainfall > soil storage, excess overland flow Q P e P 2 e S

14 Why not use a distributed storage directly in the water balance? If done wisely, it should give same overall S and runoff distribution

15 Available Storage Defining the distribution of available storage or water deficit The difference between saturation and field capacity is the availability storage or water deficit σ S( 1/(1 As ) 1) It is a function of the fractional area, A s and the total storage, S in the watershed Area (fraction)

16 Rainfall Effective Total

17 SWAT-VSA or SWAT-WB Thus by assigning field capacities according to FldCap Saturon ` S( 1/(1 As ) 1) and not according to soil type, runoff can be determined as saturation excess. Expectsame hydrograph at outlet

18 terrain analysis Distribute runoff areas in the landscape DEM Upstrean area, A Slope β soil depth, D soil hydraulic conductivity K The location of the saturated areas is mainly determined by the topography ln tan A K s D Soil topographic wetness index

19 STI SSURGO Landuse SWAT-VSA defines HRUs as the coincidence of soil topographic index (and soil) and landuse Weighted average of soil properties nested with in an area weighted index class HRUs So runoff/p loss is now not the same here (lowland pasture) As here (upland pasture)

20 Soil Topographic Index 1=f(S) =f(s) =f(s) =f(s) Soil Topographic Index 33.1 High : Low : =f(s) =f(s) =f(s) =f(s) =f(s) f(s) Lyon et al. 24. Hydrol. Proc. 18(15): Schneiderman et al. 26. Hydrol. Proc. (in press)

21 Application of two SWATs to the Cannonsville watershed

22 Cannonsville Watershed - Modeling Options Modeling Option SWAT VSA ICRK: crack flow Turned off Turned off IRTE: water routing method Variable Variable IDEG: channel degradation code Turned off Turned on IWQ: in-stream water quality Turned off Turned off

23 Output at watershed outlet Flow (cms) Flow Measured SWAT VSA Jan-94 Mar-94 May-94 Jul-94 Sep-94 Nov-94 Jan-95 Mar-95 May-95 Jul-95 Sep-95 Nov-95 Jan-96 Mar-96 May-96 Jul-96 Sep-96 Nov-96 Jan-97 Mar-97 May-97 Jul-97 Sep-97 Nov-97 Jan-98 Mar-98 May-98 Jul-98 Sep-98 Nov-98 Jan-99 Mar-99 May-99 Jul-99 Sep-99 Nov-99 TDP (kg) TDP Sediment (tons) Sediment PP (kg) PP

24 Cannonsville Watershed - Modeling Options Modeling Option SWAT VSA ICRK: crack flow Turned off Turned off IRTE: water routing method Variable Variable IDEG: channel degradation code Turned off Turned on IWQ: in-stream water quality Turned off Turned off

25 Output at Subbasin Flow (cms) Flow SWAT VSA Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma Ma TDP (kg) TDP Sediment (tons) Sediment PP (kg) PP

26 Cannonsville Watershed - Modeling Options Modeling Option SWAT VSA ICRK: crack flow Turned off Turned off IRTE: water routing method Variable Variable IDEG: channel degradation code Turned off Turned on IWQ: in-stream water quality Turned off Turned off

27 Output at Subbasin Flow (cms) Flow S W TDP (kg) TDP Sediment (tons) Sediment PP (kg) PP

28 Average Values - Values at subbasin 15 Flow PP TDP Sediment SWAT VSA Values at subbasin 4 Flow PP TDP Sediment SWAT VSA

29 Application to Townbrook.

30 Streamflow (cm) SWAT Streamflow (cm) SWAT a r 2 =.76 E =.82 Measured SWAT-VSA Measured Oct-98 Mar-99 Aug-99 Jan- Jun- Nov- Apr-1 Sep-1 15 Feb-2 Jul-2 Dec-2 May-3 Oct-3 Mar-4 Aug-4 14 b Measured r 2 =.74 E =.83 SWAT Measured 15 2 Oct-98 Mar-99 Aug-99 Jan- Jun- Nov- Apr-1 Sep-1 Feb-2 Jul-2 Dec-2 May-3 Oct-3 Mar-4 Aug-4

31 Runoff SWAT-VSA swat_vsa_hrus SWAT_VSA SURQ Runoff [mm] (mm) swat_hrus SWAT_Standard SURQ Runoff [mm] (mm) SWAT-Standard

32 Soil Water SWAT-VSA SWAT_VSA SW_END Soil Water [mm] (mm) SWAT-Standard SWAT_Standard Soil SW_END Water [mm] (mm)

33 Cumulative Runoff (mm) SWAT-VSA Runoff Index1 Index9 Index8 Index7 Index6 Index5 Index4 Index3 Index2 Index1 Pasture 5 Jan-96 Jan-97 Jan-98 Jan-99 Jan- Jan-1 Jan-2 Jan-3 Jan-4 Jan-5

34 Cumulative Runoff (mm) SWAT-Standard Runoff 25 2 Pastures Jan-96 Jan-97 Jan-98 Jan-99 Jan- Jan-1 Jan-2 Jan-3 Jan-4 Jan-5

35 Phosphorus SWAT-VSA swat_vsa_hrus SWAT_VSA ORGP Dissolved KG/HA P (kg ha -1 ) SWAT-Standard swat_hrus SWAT_Standard ORGP Dissolved KG/HA P (kg ha -1 )

36 Dissolved P (kg) Dissolved P (kg) SWAT SWAT 8, 6, a Measured SWAT-VSA r 2 =.76 E =.71 3 Measured 3 4, 2, Base Calendar Year Cum dis nut[p] : ld_townbrook_measured Cum dis nut[p] : ld_townbrook_swat_standard Cum dis nut[p] : ld_townbrook_swat_vsa 5 25 Oct-98 Apr-99 Oct-99 Apr- Oct- Apr-1 Oct-13Apr-2 Oct-2 Apr-3 Oct-3 Apr-4 b Measured 2 SWAT r 2 =.68 E =.47 Measured 3 Oct-98 Apr-99 Oct-99 Apr- Oct- Apr-1 Oct-1 Apr-2 Oct-2 Apr-3 Oct-3 Apr-4

37 Conclusions Both SWAT models give similar watershed outlet predictions For effective BMP implementation, hotspots should be identified correctly Validating distributed output is difficult