Quantifying climate change impacts on streamflow ...

5 downloads 0 Views 1MB Size Report
Grand river at Brantford: -PBIAS = 2% (very good). -NSE = 0.77 (good). -R2 = 0.88 (very good). Themes river at. Thamesville: -PBIAS = 1% (very good).
Quantifying climate change impacts on streamflow dynamics of Grand and Thames river basins, Canada Binbin Zhang, Narayan Kumar Shrestha, Prasad Daggupati, Ramesh Rudra, Rituraj Shukla, Baljeet Kaur, Jun Hou School of Engineering, University of Guelph, ON, Canada

The Great Lakes and the Grand and Thames river basins



The Northern Lake Erie (NLE) drains 21750 km2



The Grand river drains 6764 km2



The Thames drains 5883 km2



Agriculture is the dominant land use type (44% and 69%) respectively for the Grand and Thames river basin

Climate Change in Cold Climate Regions and in Great Lakes     

The Great Lakes provides fresh water resources for consumption, transportation, tourism, power generation, among others The Lake Erie is home to approximately 14 million In the last few decades, serious blue-green algae blooms has been observed in the Lake Erie Growing evidences that climate change is a real In cold climate region: o o o o

o o o

Increases in winter temperature decreases in snowfall proportion on total precipitation changes in lake level changes in streamflow input to the lake changes in freeze-thaw frequencies increases in drought and flood frequency changes in ice cover, among others

The Grand and Thames river basins  



Soil and Water Assessment Tool (SWAT) Precipitation and Temperature data: Gridded Climate Data for Canada (GCDC) Base period: 1980-1993 (1983-1988: Calibration, 1989-1993: Validation)



Streamflow at four stations used  Bias correction tool: Climate Model data for Hydrologic modeling (CMhyd)  Environment Canada (EC) station: Glen Allan is used for assessing performance of different bias correction methods

Different Bias Correction Methods and their Evaluation 



Methods o

Linear Scaling of Precipitation and Temperature

o

Local Intensity Scaling (LOCI) of Precipitation

o

Power Transformation (PT) of Precipitation

o

Variance Scaling (VS) of Temperature

o

Distribution Mapping (DM) of Precipitation and Temperature

Evaluation o

Frequency based statistics

o

Time Series based statistics

o

Coefficient of Determination (R2)

o

Nash-Sutcliffe Efficiency (NSE)

o

Percentage of Bias (PBIAS)

The Climate Change Data 

Canadian regional climate model (CanRCM4)



Resolution: 0.22°, equivalent to 25 km



The projections from the CanRCM4 are tailored made for applications in Canada



CanRCM4 is one of RCMs used in the Coordinated Regional Climate Downscaling Experiment (CORDEX) project implemented by the World Climate Research Programme (WCRP)



Two emission scenarios: the Representative Concentration Pathway RCP 4.5 and RCP 8.5



Two future periods: a mid-century (2043-2053) and an end-century (2083-2093)

SWAT Model Results for Streamflow during Base Period Simulated (Calibration)

0 100

300 200 200

100 0 500 400

Streamflow (m3/s)

Simulated (Validation)

300

-PBIAS = 2% (very good) -NSE = 0.77 (good) -R2 = 0.88 (very good)

0 100

200 200

0

Grand river at Brantford:

400

300

100

Precipitation (mm)

Streamflow (m3/s)

400

Observed Flow

300

Precipitation (mm)

Precipitation

500

Themes river at Thamesville: -PBIAS = 1% (very good)

-NSE = 0.88 ( very good) -R2 = 0.91 (very good)

400

Bias Correction Results for Precipitation 

Distribution Mapping (DM) was found to be the best performing method

Bias Correction Results for Temperature 

Distribution Mapping (DM) was again found to be the best performing method

Different Combination of Bias Correction Methods for Streamflow Simulation at four Stations Simulations/ Combinations

Grand river near Marsville PBIAS (%) NSE R2

Grand river at Brantford PBIAS (%) NSE R2

Thames river at Ingersoll PBIAS (%) NSE R2

0.79 6.42 0.81 0.86 -4.25 0.94 0.91 -10.9 0.93 Baseline 0.59 -2.01 0.69 0.64 4.77 0.68 0.22 11.49 0.61 Raw 0.3 -36.69 0.62 0.52 -36.66 0.77 0.46 -44 0.68 P(LS)+T(LS) 0.25 -35.99 0.62 0.51 -36.6 0.76 0.4 -43.36 0.67 P(LS)+T(VS) 0.28 -36.72 0.6 0.52 -36.78 0.77 0.43 -43.24 0.67 P(LS)+T(DM) 0.48 -19.95 0.67 0.73 -18.86 0.83 0.6 -26.7 0.73 P(LOCI)+T(LS) 0.43 -19.03 0.67 0.72 -18.62 0.81 0.54 -25.91 0.71 P(LOCI)+T(VS) 0.46 -19.97 0.65 0.73 -18.93 0.83 0.58 -25.95 0.73 P(LOCI)+T(DM) 0.43 -41.55 0.72 0.49 -43.77 0.77 0.44 -52.83 0.7 P(PT)+T(LS) 0.41 -41.11 0.74 0.48 -43.92 0.77 0.39 -52.39 0.69 P(PT)+T(VS) 0.42 -41.63 0.72 0.49 -43.97 0.77 0.42 -52.23 0.7 P(PT)+T(DM) 0.57 -7.04 0.63 0.83 -6.2 0.84 0.72 -14 0.78 P(DM)+T(LS) 0.56 -6.13 0.64 0.83 -6.19 0.85 0.69 -13.45 0.78 P(DM)+T(VS) 0.58 -7.07 0.63 0.84 -6.48 0.85 0.72 -13.44 0.79 P(DM)+T(DM) P: Precipitation; T: Temperature LS: Linear Scaling; VS: Variance Scaling; DM: Distribution Map; LOCI: Local Intensity Scaling; PT: Power Transformation



P(DM)+T(DM) was the best and used for further processing

Thames river at Thamesville PBIAS (%) NSE R2 0.98 0.83 0.78 0.83 0.76 0.86 0.86 0.85 0.78 0.77 0.76 0.92 0.92 0.91

1.63 -4.01 -29.81 -22.07 -29.57 -15.35 -14.68 -15.16 -33.81 -33.37 -33.75 -4.95 -4.34 -4.9

0.97 0.81 0.52 0.56 0.5 0.75 0.73 0.73 0.56 0.54 0.54 0.9 0.9 0.89

Future Changes in Precipitation

-40%

+55%

Future Changes in Mean Temperature

up to +12.5oC

Future Changes in Green Water Flow (Evapotranspiration: ET)

Future Changes in Green Water Storage (Soil Water: SW)

water stress conditions to growing crops

Future Changes in Streamflow Dynamics Historical

RCP4.5

RCP4.5

100

0

-50

120

-100

80

-150 -200

40

-250

0

Flow Discharge (m3/s)

50

160

-300

RCP4.5

100 50 0

-50

120

-100

80

-150 -200

40

-250

-300

RCP4.5

RCP8.5

200

200

160

100

0 120 -100

80 -200 40

-300

0

-400 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

RCP8.5

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Flow Discharge (m3/s)

Historical

RCP4.5

0

Percentage Changes (%)

Flow Discharge (m3/s)

RCP8.5

Historical

160

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

RCP4.5

RCP8.5

200

Percentage Changes (%)

Flow Discharge (m3/s)

200

RCP8.5

Percentage Changes (%)

RCP8.5

RCP8.5

Historical

RCP4.5

RCP8.5

200

200

160

100

0 120 -100

80 -200 40

-300

0

-400

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

Percentage Changes (%)

RCP4.5

Uncertainty Future Changes in Streamflow Dynamics Grand river near Marsville

Thames river at Ingersoll

Grand river at Brantford

Thames river at Thamesville

Key Findings 

We found the combination of distribution mapping (DM) methods for both precipitation and maximum temperature to be the best method



Results showed marked temporal and spatial variability on future changes in precipitation and temperature



Also reflected in evapotranspiration and soil water storage changes, leading to heterogeneity in streamflow changes at four considered stations



On average, increases in winter (+11%) and decreases in spring (-33%), summer (-23%) and autumn (-15%) streamflow were expected in future



Region-specific temporal variability in water resources availability, as expected in future in these river basins, pose considerable challenges to water planners and managers and to integrated water resources management of these two river basins