A copula-based multivariate analysis of Canadian RCM projected changes to flood characteristics for northeastern Canada

In the present work, climate change impacts on three spring (March–June) flood characteristics, i.e. peak, volume and duration, for 21 northeast Canadian basins are evaluated, based on Canadian regional climate model (CRCM) simulations. Conventional univariate frequency analysis for each flood characteristic and copula based bivariate frequency analysis for mutually correlated pairs of flood characteristics (i.e. peak–volume, peak–duration and volume–duration) are carried out. While univariate analysis is focused on return levels of selected return periods (5-, 20- and 50-year), the bivariate analysis is focused on the joint occurrence probabilities P1 and P2 of the three pairs of flood characteristics, where P1 is the probability of any one characteristic in a pair exceeding its threshold and P2 is the probability of both characteristics in a pair exceeding their respective thresholds at the same time. The performance of CRCM is assessed by comparing ERA40 (the European Centre for Medium-Range Weather Forecasts 40-year reanalysis) driven CRCM simulated flood statistics and univariate and bivariate frequency analysis results for the current 1970–1999 period with those observed at selected 16 gauging stations for the same time period. The Generalized Extreme Value distribution is selected as the marginal distribution for flood characteristics and the Clayton copula for developing bivariate distribution functions. The CRCM performs well in simulating mean, standard deviation, and 5-, 20- and 50-year return levels of flood characteristics. The joint occurrence probabilities are also simulated well by the CRCM. A five-member ensemble of the CRCM simulated streamflow for the current (1970–1999) and future (2041–2070) periods, driven by five different members of a Canadian Global Climate Model ensemble, are used in the assessment of projected changes, where future simulations correspond to A2 scenario. The results of projected changes, in general, indicate increases in the marginal values, i.e. return levels of flood characteristics, and the joint occurrence probabilities P1 and P2. It is found that the future marginal values of flood characteristics and P1 and P2 values corresponding to longer return periods will be affected more by anthropogenic climate change than those corresponding to shorter return periods but the former ones are subjected to higher uncertainties.