Floods in Future Climates - Process and Statistical Issues Paul H. Whitfield Meteorological Service of Canada, Environment Canada, Vancouver Department of Earth Sciences, Simon Fraser University, Burnaby
Floods in Future 4AR views What are floods? Mechanisms What causes floods? Will that change in future Measurement Issues Modelling issues Equifinality Extrapolation beyond observations Scale and resolution
4AR 2007 IPCC Increased precipitation intensity and variability are projected to increase the risks of flooding and drought in many areas. The frequency of heavy precipitation events (or proportion of total rainfall from heavy falls) will be very likely to increase over most areas during the 21st century, with consequences for the risk of rain-generated floods.
Fig 10. 9
Uncertainties arise from the incorporation of climate model results into freshwater studies for two reasons: the different spatial scales of global climate models and hydrological models biases in the long-term mean precipitation as computed by global climate models for the current climate Timing biases from GCM patterns A number of methods have been used to address the scale differences: simple interpolation of climate model results dynamic or statistical downscaling methods all such methods introduce uncertainties into the projection Biases in simulated mean precipitation are often addressed by adding modelled anomalies to the observed precipitation in order to obtain the driving dataset for hydrological models. changes in interannual or day-to-day variability of climate parameters are not taken into account in most hydrological impact studies. Hydrology has both short-term and long-term memory Over/underestimation of future floods, droughts and irrigation water requirements.
3.39
Impacts on Europe Increasing risk of winter flood in northern Europe and of flash flood in all of Europe Risk of snowmelt flood shifts from spring to winter Today s 100-year floods are projected to occur more frequently northern and north-eastern Europe (Sweden, Finland, N. Russia) Ireland, Central and E. Europe (Poland, Alpine rivers), Atlantic parts of S. Europe (Spain, Portugal); less frequently in large parts of S. Europe Globally - wet areas/periods get wetter and dry areas/periods get drier
Uncertainties of Climate Change Science Increasing GHG Concentrations Significance of Future Climate Change Magnitude/Rate of Global Temperature, Precipitation and Sea Level Change Discernible human Influence Amplified Polar/Continental warming Change in Regional Rainfall Patterns/Extreme Event Detailed Characteristics of Local Change Local Impacts High Confidence Level Low
What are floods? a rising and overflowing of a body of water especially onto normally dry land ; also : a condition of overflowing <rivers in flood. Flooding is primarily but not exclusively caused by hydro-meteorological conditions, acting either individually or in combination
Studying Future Floods Landuse Assumptions Future Climates Extreme events Pattern dependencies Uncertainty Issues Measurements Statistical Distribution return periods Equifinality Process changes Scaling effects Coupling models affected by central limit it theorem Extrapolation beyond observations Examples
10.18
Arnell et al 2004. Global Environmental Change 14:3-20
What Causes a Flood? Combination of circumstances that t results in an excess of water In a stream channel On a surface Combinations of precipitation and antecedent conditions Combinations of temperature and antecedent conditions Combination of conditions and landuse
Precipitation Unresolved in GCMs Downscaling Changes in Intensity Variety of processes Extrapolation beyond the observed
Scope and nature of precipitation p processes Process Scope Scale Orographic Windward sides on elevation features Locally persistent local Convection Differential heating of atmosphere local Km s Tropical Cyclonic Convergence in cyclones Tropical regional Extra-tropical Pineapple Express mid-latitude regional Frontal regional Monsoon Heavy precipitation alternates with hot, dry conditions on an annual basis due to seasonal reversal of winds widespread continental Unlikely that observed patterns and frequencies in the past 100 years will continue
Pineapple Express
Flood Processes Timing Storm Duration Rainfall depths Catchment Runoff Spatial state response coherence Long -rain No seasonality >1day Substantial wet due to persistent rainfall Short-rain No seasonality Hours to 1 day Moderate to substantial Flash flood Periods of <90min Small to extreme rainfall moderate Rain on snow Cold to warm transitions Moderate rainfall events Snowmelt and rainfall Slow Extent of storm Wet for large Fast Local or regional flood events Any Flashy Limited <30km Wet, snow covered Snow-melt Spring to Unimportant Minor rainfall Wet, snow summer covered Glacial outburst Coastal flood Event Event Fast to slow Medium to slow Areas of snow cover Medium Hurricanes and tropical storms can produce heavy rains, or drive ocean water onto land. Urban flood Various Rain water can not be absorbed into the ground and becomes runoff, filling parking lots, making roads into rivers, and flooding basements and businesses. Ice-jam flooding Groundwater Rise Event
Figure 7. Change in extremes by 2050, under HADCM3. This panel shows change in the magnitude of the 10-year return period maximum monthly runoff. Arnell et al 1999. Global Environmental Change 9:S31-S49
GCM Prudhomme et al 2003 J Hydrology 277:1-23
Measuring Flood Flows
Norish Creek near Dewdney 08MH058 2.5 2 >3Q >2h Stage (m) 1.5 1 0.5 0 0 20 40 60 80 100 120 140 Q(m3/s)
How are Floods Statistically Distributed? Standard Practice Log Pearson Type 3 GEV Annual series partial series Mixed populations Multiple mixtures
Alila and Mtiraoui 2002 Hydrological Processes 16:1065-1084
Alila and Mtiraoui 2002 Hydrological Processes 16:1065-1084
Equifinality Models with many parameters have equally likely solutions that fit observations equally well Within simple models parameters may lack clear physical meaning Within complex models the physical meaning is improve but not better resolved
Cameron et al. 2000 Hydrology and Earth System Sciences 4:393-405
Booji 2005 J Hydrology 303:176-198 Scaling Effects
Coupling Models Prudhomme et al. 2002; Hydrological Processes 16:1137-1150
Fig. 4. Errors of the simulated peak flows for peaks of different magnitudes during the test period for the four catchments. Both the median and the 80-per cent range of the prediction errors obtained using the 50 best parameter sets are shown. The assemblies of best parameter set were determined with regard to the model efficiency (Reff, left column), the mean of efficiency and goodness of groundwater level simulations (Reff and r2, middle), as well as the mean of efficiency and goodness of peak flow simulations (Reff and Rpeak, right). The shaded area indicates where the model was extrapolated, i.e., events that were larger than those observed during the calibration period. Seibert 2003 Nordic Hydrology 34:477-492 HBV
Seibert 2003 Nordic Hydrology 34:477-492 HBV
Streamflow Processes e m3/sec Discharg 60 Englishman 600 Capilano Snowmelt Coquitlam Fall Rains 50 Cheakamus 500 40 30 20 Elaho Lillooet Winter Rains Glacier melt Dry 10 Summer 100 400 300 200 0 1 6 11 16 21 26 31 36 41 46 51 56 61 66 71 0 After Wade et al. 2001 observations
ENSO - Streamflow Fleming et al. 2008 Statistical bootstrapping
PDO - Streamflow Fleming et al. 2008 Statistical bootstrapping
Future Changes Coquitlam River 18 Discharg ge (cms) 16 1973-1993 2013-2033 p<0.05 14 2043-2063 p<0.05 12 2073-2093 p<0.05 10 8 6 4 2 0 1 6 11 16 21 26 31 36 41 46 51 56 61 66 71 5 Day Periods Whitfield et al 2002 UBC Watershed - snowmelt
Floods and Process Change in % Coastal Rainfall Snowmelt Upper Campbell River Illecillewaet River Episodes 11.4-23.3 Flood Days 59.2-27.5 Duration 44.44 2.6 Flood Volume 93.7-38.2 Mean Flood Flow 4.1-5.2 Mean Flood Peak 14.11-7.0 70 Annual Maximum 16.1-10.7 Day of Flood Centroid -2.7-7.7 Day of Occurence of Peak -16.9-6.4 Loukas and Quick, 1999 UBC Watershed Model - snowmelt
Statistical Models Model validation- Process & scale effects Index of agreement- Decadal change Index of agreement- Mann- Whitney U Inde ex of agreeme ent (decadal ch hange) 1.00 0.90 0.80 0.70 0.60 0.50 0.40 0.30 0.20 010 0.10 0.00 0.74 0.66 b c d 0.86 0.83 0.79 0.76 0.74 0.52 0.83 0.78 0.64 0.58 0.76 0.78 0.53 0.85 081 0.81 08 0.8 420 453 1170 3080 9710 10200 11500 14600 6.22 80.5 201 230 337 1850 2430 37.3 172 181 185 355 477 Drainage Area (sq. km) 0.64 0.84 0.77 ney U) Inde ex of agreemen nt (Mann-Whit b c d 1.00 0.91 0.9 0.91 0.87 0.87 0.88 0.89 090 0.90 085084 0.81 0.85 0.86 0.84 0.83 0.81 0.83 0.86 0.86 0.83 0.78 0.78 0.80 0.75 0.70 0.60 050 0.50 0.40 0.30 0.20 0.10 0.00 0.53 420 453 1170 3080 9710 10200 11500 14600 6.22 80.5 201 230 337 1850 2430 37.3 172 181 185 355 477 Drainage Area (sq. km) Cannon and Whitfield, 2002 Neural Network Models
1.20 1.00 0.80 0.60 0.40 Statistical Models Model validation- Process & scale effects Predicted mean/obs. mean Predicted variance/obs. variance b c d 107 1.07 1.03 0.99 0.99 0.99 0.99 1 0.98 0.94 0.97 0.99 0.94 0.95 0.91 0.89 0.8 0.83 0.84 0.77 0.78 0.75 b c d 1.20 1.09 1.07 1 0.99 1 0.95 0.95 0.97 1.00 0.94 0.88 0.79 0.81 0.80 0.68 0.60 0.53 0.53 0.54 0.37 0.40 0.36 0.28 0.26 0.21 020 0.20 0.00 Predicted mean/observed mea n 420 453 1170 3080 9710 10200 11500 14600 6.22 80.5 201 230 337 1850 2430 37.3 172 181 185 355 477 420 453 1170 3080 9710 10200 11500 14600 6.22 80.5 201 230 337 1850 2430 37.3 172 181 185 355 477 Predicted variance/ /observed variance 0.20 0.00 Drainage Area (sq. km) Drainage Area (sq. km) Cannon and Whitfield, 2002 Neural Network Models
Georgia Basin Model Bias 1.2 Magnitude of 10 year Flood 0.8 0.4 Englis shman ilano Capi Coqu uitlam Chea akamus Ela aho Lillo ooet UBC Watershed - snowmelt 0 High Low Decreasing Rainfall Influence Whitfield et al, 2003 Relative to Mean Observed Magnitude of 10 year Flood Observed, Danard, NCEP, CGCM1
Flood Days Percent 16 14 12 10 8 Generating Mechanism Frequency Location Magnitude 6 4 2 0 Englishman Capilano Coquitlam Cheakamus Elaho Lillooet Modelled 2013 2043 2073 Percent of Events Greater than Q05 Whitfield et al, 2003 UBC Watershed - snowmelt
Location Flows above Q05 Modelledd d 2013 Englishman River Modelled 2013 Capilano River 2043 2043 2073 2073 1-Jan 31-Jan 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec 1-Jan 31-Jan 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec UBC Watershed - snowme elt Coquitlam River Cheakamus River Modelled Modelled 2013 2013 2043 2043 2073 2073 1Jan 1-Jan 31-Jan 1Mar 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec 1Jan 1-Jan 31-Jan 1Mar 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec Elaho River Lillooet River Modelled Modelled 2013 2013 2043 2043 2073 2073 1-Jan 31-Jan 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec 1-Jan 31-Jan 1-Mar 31-Mar 30-Apr 30-May 29-Jun 29-Jul 28-Aug 27-Sep 27-Oct 26-Nov 26-Dec Whitfield et al, 2003
120 10 Year - Englishman River 200 10 Year - Capilano River 10 Year Annual Series Gumbel; Log Pearson Type 3 80 40 160 120 80 40 0 1960 2000 2040 2080 10 Year - Coquitlam River 50 40 30 20 10 0 1960 2000 2040 2080 10 Year - Cheakamus River 160 120 80 40 UBC Watershed - snowmelt 0 1960 2000 2040 2080 0 1960 2000 2040 2080 800 10 Year - Elaho River 1000 10 Year - Lillooet River 800 600 600 400 400 200 200 Whitfield et al, 2003 0 1960 2000 2040 2080 0 1960 2000 2040 2080
10 Year - Englishman River - Partial Series 10 Year - Capilano River - Partial Series 120 160 10 Year Partial Series Gumbel; Log Pearson Type 3 80 40 120 80 40 0 1960 2000 2040 2080 10 Year - Coquitlam River - Partial Series 50 40 30 20 10 0 1960 2000 2040 2080 0 1960 2000 2040 2080 10 Year - Cheakamus River - Partial Series 200 160 120 80 40 0 1960 2000 2040 2080 UBC Watershed - snowmelt 800 10 Year - Elaho River - Partial Series 800 10 Year - Lillooet River - Partial Series 600 600 400 400 200 200 Whitfield et al, 2003 0 1960 2000 2040 2080 0 1960 2000 2040 2080
CGCM2 emissions A2 Englishman Riv er near Parksville - 08HB002 70 60 50 Mean - ENGLISH CGCM2A (1973-1993) Mean - ENGLISH CGCM2A (2013-2033) p<0.05 Mean - ENGLISH CGCM2A (2043-2063) p<0.05 Mean - ENGLISH CGCM2A (2073-2093) p<0.05 Disch harge (m3/sec) 40 30 20 10 0 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 IHACRES Rainfall model
CGCM2 emissions B2 Englishman River near Parksville - 08HB002 Disc charge (m3/sec) 60 50 40 30 20 Mean - ENGLISH CGCM2B (1973-1993) Mean - ENGLISH CGCM2B (2013-2033) p<0.05 Mean - ENGLISH CGCM2B (2043-2063) p<0.05 Mean - ENGLISH CGCM2B (2073-2093) p<0.05 10 0 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 IHACRES Rainfall model
Annual Series 10-Year Floods IHACRES Rainfall model
Partial Series 10-Year Floods IHACRES Rainfall model
10-, 100-, 200-Year Floods CGCM2B Annual series (annual maximum) IHACRES Rainfall model
10-, 100-, 200-Year Floods CGCM2B Partial series (all events above Q 05 )
Ice-induced Flood Frequency Liard River at Mouth Ice effect return period of ice-jam floods tend to be lower than for open-water events extremes can extend into flood levels not considered for open-water events rarely are ice events analyzed when determining flood risk numerous flood analyses over the last century of minimal value because analysis conducted on discharge not water levels (ignore ice)
TRENDS IN NORTHERN RIVER-ICE BREAKUP AND TIMING OF SPRING 0 C-ISOTHERMS A B 25 20 C days/49 years 15 10 5 0-5 -10 (1950-98) 98) -15-20 D -25 earlier in the west later in the east matches some W- E GCM forecasts western trends match even longer term historical breakup records Bonsal & Prowse 2003 Climatic i Change Bonsal et al., 2001
Summary Shorter duration of precipitation season coupled with overall increase in precipitation intensity Changing dominance of snow to rainfall Larger volumes of water in shorter time periods More runoff in winter Less spring runoff
Not all floods are equal Importance of precipitation process as generating g force Changes in mixture of driving forces Changes in generating processes Modelling issues Atmosphere Hydrological models Equifinality