Application of SWAT to Quantify Internal Renewable Water Resources in Iran

Transcription

1 SWAT 2007,UNESCO IHE, the Netherlands, 5-Jul-07 Application of SWAT to Quantify Internal Renewable Water Resources in Iran Monireh Faramarzi Eawag: Swiss Federal Institute of Aquatic Science and Technology Picture: Doosti dam, central Khorasan province

2 Content General background Objectives Methodology model area model input model setup Results

3 Background Objectives Methodology Results and outlook General overview of the Project The proposed research project aims to assess the feasibility of applying the virtual water strategy to alleviate water stress in Iran. Regional crop structure adjustment Inter-provincial food trade Socio-economic factors improving regional/provincial and national water productivity and Water use efficiency

4 Background Objectives Methodology Results and outlook Start with SWAT: Quantification of regional water resources Quantification of provincial water resources

5 Turkey Azerbaijan Caspian Sea Area:,648,000 km 2 (65 million hectare) Altitudes: -80 m to 5670 m Average annual precipitation: 252 mm Precipitation range: mm Temperature: -44 to 56 degree C. Iraq Precipitation Afghanistan Pakistan Persian Gulf

6 Aras Dam Karaj Dam Model area and main Zayandeh Rud Dam Max. Volume: 960 MCM Max volume: MCMHydrologic regions Max volume: MCM Operation: 97 Operation: 96 Operation: 970 Karun Dam Max. volume: MCM Caspian Sea Operation: 977 Karkheh Dam Max. Volume: MCM Operation: Rajaee Dam Max volume: 96.0 MCMO peration: Karun Dam Operation: Latyan Dam 962 Persian Gulf Max volume: 99 MCM Sefid Rud Dam, Max volume: 840 MCM 967 Operation:

7 Model input (Global data) Landuse (Extracted from global USGS landuse/land cover) DEM (Extracted from global USGS DEM map) High: 5670 m Low: -80 m : URMD 2: CRDY 3: CRIR 5: CRGR 6: CRWO 7: GRAS 8: SHRB 9: MIGS 0: SAVA : FODB 2: FODN 4: FOEN 5: FOMI 6: WATB 7: WEHB 8: WEWO 9: BSVG 2: TUWO 22: TUMI Soil (Extracted from global FAO soil map, 995)

8 Synoptic stations (providing daily precipitation, max. and min. temperature) Model input (Local data) Hydrometric stations (providing river daily discharge data) Lake-inland water & River map Large reservoirs map River.shp Inland_water.shp River Lake_inland water N Meters

9 Model setup. Arc-GIS (Olivera et al., 2006) interface was used to parameterize whole the area 2. Based on DEM and stream network whole the area was divided into 506 sub-basins (threshold drainage area was set to 600 km 2 ) 3. Dominant soil and land use was selected to provide soil and land use data in the model large reservoirs/dams were contributed to the model 5. Considered simulation period for calibration was 6 years from considering 3 years warm up period

10 Background Objectives Methodology Results and outlook Definition of three approaches used in calibration processes using SUFI-2:. Global approach Global parameterization 2. Scaling approach Global parameterization 3. Regional approach Global parameterization Change of each parameter was applied globally Change of each parameter was applied based on landuse and soil WR_ WR_2 WR_3 WR_4 WR_5 WR_6 WR_7 WR_ global g opt = φi ; 8 i= i= 6 g = φ i 6 i = 0 g 2 = φ i 0.. i = 8 g 8 = φ i 8 i= parameters scaled g opt = φ i 8 ; parameters (Krause et al., 2005): 2 = br b R φ 2 for for < 0 b > b

11 Background Objectives Methodology Results Average weighted coefficient of determination (g) Global approach Regional approach Scaling approach HR_ HR_2 HR_3 HR_6 HR_5 HR_4 Hydrologic region HR_7 HR_8 Country Average weighted coefficient of determination for the eight hydrologic regions and the whole country

12 Background Objectives Methodology Results φ Weighted coefficient of determination (φ) at 8 stations across the Country Paye-Pol watershed area: km PPU 600 Best_Simulation Observed φ:0.65 P-factore: 0.99 R-factor: Date Pole-Kaleh watershed area: 5677 km 2 Pole-Choom watershed area: 80 km PPU Best_Sim Observed φ:0.32 P-factore: 0.33 R-factor: PPU Best_Sim Observed φ:0.60 P-factore: 0.70 R-factor: Date Stream flow (m 3 s - ) Stream flow (m 3 s - ) Stream flow (m 3 s - )

13 Background Objectives Methodology Results 95 PPU of regional renewable water resources (First round of calibration) Jamab data 95 PPU 95 PPU Jamab estimation Renewable water resources IRWR (km (km 3 /year) 3 /year) 0 HR_8 WR_8 HR_7 WR_7 HR_6 WR_6 WR_5 HR_5 HR_4 WR_4 HR_3 WR_3 WR_2 HR_2 WR_ HR_ Hydrologic regions

14 Background Objectives Methodology Results Final results from the final parameter set (Regional approach):. Internal renewable blue water resources and actual evapotranspiration at provinces IRWR, ET (mm year - ) Blue w ater (IRWR) Evapotranspiration (ET) 0 Ardebil Azarbayjan_Gharbi Azarbayjan_Sharghi Boshehr Charmahal Esfahan Province Fars Ghazvin Gilan Golestan

15 Background Objectives Methodology Results IRWR, ET (km 3 month - ) Soil water (km 3 ) Jan Province area: km 2 Mar May Jul Sep Nov Jan Month Province area: km 2 Mar 95 PPU - Blue water 95 PPU - Evapotransiration May Jul Sep Nov 95 PPU - Soil water 3. Average ( ) monthly internal renewable blue water (IRWR), Actual evapotranspiration (ET) and Soil water for Khozestan province 3 0 Jan Mar Mai Jul Sep Nov Jan Mar Mai Jul Sep Nov Month

16 Background Objectives Methodology Results 2. Blue water resources at sub-province level (was constructed based on final parameter set resulted from Regional approach ) Thanks for your attention