Votre recherche

Droughts are increasingly recognized as a significant global challenge, with severe impacts observed in Canada's Prairie provinces. While less frequent in Eastern Canada, prolonged precipitation deficits, particularly during summer, can lead to severe drought conditions. This study investigates the causes and consequences of droughts in New Brunswick (NB) by employing two drought indices: the Palmer Drought Severity Index (PDSI) and Standardized Evapotranspiration Deficit Index (SEDI)– at ten weather stations across NB from 1971 to 2020. Additionally, the Canadian Gridded Temperature and Precipitation Anomalies (CANGRD) dataset (1979–2014) was utilized to examine spatial and temporal drought variability and its alignment with station-based observations. Statistical analyses, including the Mann–Kendall test and Sen's slope estimator, were applied to assess trends in drought indices on annual and seasonal timescales using both station and gridded data. The results identified the most drought-vulnerable regions in NB and revealed significant spatial and temporal variability in drought severity over the 1971–2020 period. Trend analyses further highlighted the intensification of extreme drought events during specific years. Coastal areas in southern NB were found to be particularly susceptible to severe drought conditions compared to inland regions, consistent with observed declines in both the frequency of rainy days and daily precipitation amounts in these areas. These findings underscore the need for targeted drought mitigation strategies particularly in NB’s coastal zones, to address the region’s increasing vulnerability to extreme drought events.

Abstract The present study analyses the impacts of past and future climate change on extreme weather events for southern parts of Canada from 1981 to 2100. A set of precipitation and temperature‐based indices were computed using the downscaled Coupled Model Intercomparison Project Phase 5 (CMIP5) multi‐model ensemble projections at 8 km resolution over the 21st Century for two representative concentration pathway (RCP) scenarios: RCP4.5 and RCP8.5. The results show that this region is expected to experience stronger warming and a higher increase in precipitation extremes in future. Generally, projected changes in minimum temperature will be greater than changes in maximum temperature, as shown by respective indices. A decrease in frost days and an increase in warm nights will be expected. By 2100 there will be no cool nights and cool days. Daily minimum and maximum temperatures will increase by 12 and 7°C, respectively, under the RCP8.5 scenario, when compared with the reference period 1981–2000. The highest warming in minimum temperature and decrease in cool nights and days will occur in Ontario and Quebec provinces close to the Great Lakes and Hudson Bay. The highest warming in maximum temperature will occur in the southern parts of Alberta and Saskatchewan. Annual total precipitation is expected to increase by about 16% and the occurrence of heavy precipitation events by five days. The highest increase in annual total precipitation will occur in the northern parts of Ontario and Quebec and in western British Columbia.

ABSTRACTTrends in indices based on daily temperature and precipitation are examined for two periods: 1948–2016 for all stations in Canada and 1900–2016 for stations in the south of Canada. These in...

Quantile estimates are generally interpreted in association with the return period concept in practical engineering. To do so with the peaks‐over‐threshold (POT) approach, combined Poisson‐generalized Pareto distributions (referred to as PD‐GPD model) must be considered. In this article, we evaluate the incorporation of non‐stationarity in the generalized Pareto distribution (GPD) and the Poisson distribution (PD) using, respectively, the smoothing‐based B‐spline functions and the logarithmic link function. Two models are proposed, a stationary PD combined to a non‐stationary GPD (referred to as PD0‐GPD1) and a combined non‐stationary PD and GPD (referred to as PD1‐GPD1). The teleconnections between hydro‐climatological variables and a number of large‐scale climate patterns allow using these climate indices as covariates in the development of non‐stationary extreme value models. The case study is made with daily precipitation amount time series from southeastern Canada and two climatic covariates, the Arctic Oscillation (AO) and the Pacific North American (PNA) indices. A comparison of PD0‐GPD1 and PD1‐GPD1 models showed that the incorporation of non‐stationarity in both POT models instead of solely in the GPD has an effect on the estimated quantiles. The use of the B‐spline function as link function between the GPD parameters and the considered climatic covariates provided flexible non‐stationary PD‐GPD models. Indeed, linear and nonlinear conditional quantiles are observed at various stations in the case study, opening an interesting perspective for further research on the physical mechanism behind these simple and complex interactions.

ABSTRACT Two composite sedimentary sequences sampled in the ice‐proximal (12CS) and ice‐distal (02CS) areas of Coronation Fjord (Baffin Island, Nunavut, Canada) were investigated in order to reconstruct the effect of climate variability on 600 years of changes in sediment transfer from the eastern Penny Ice Cap (PIC). Detrital proxies, and physical and sedimentological analyses revealed that glacial meltwater discharges led to frequent rapidly deposited layers (RDLs) in ice‐proximal settings. RDLs in ice‐distal settings involved the sudden release of a large quantity of sediment‐laden water during floods probably originating from adjacent fjords with large sandur deltas. Laminated sediments with ice‐rafted debris throughout the Little Ice Age interval in the ice‐proximal environment suggest that colder conditions promoted glacier growth, leading to successive episodes of turbid hyperpycnal meltwater plumes and iceberg calving in Coronation Fjord. Since 1850 ce , the accelerated Coronation retreat in response to modern warming has led to increased sedimentation rates, abrupt mineralogical and grain size proxy variations and more frequent RDLs. Similar trends between the detrital proxies of the ice‐proximal core and Atlantic Multidecadal Oscillation record and Arctic surface air temperature suggest high connectivity between atmospheric and sea surface temperature variations and PIC dynamics over the last 600 years.

A four-dimensional variational data assimilation (4DVAR) scheme has recently been implemented in the medium-range weather forecast system of the Meteorological Service of Canada (MSC). The new scheme is now composed of several additional and improved features as compared with the three-dimensional variational data assimilation (3DVAR): the first guess at the appropriate time from the full-resolution model trajectory is used to calculate the misfit to the observations; the tangent linear of the forecast model and its adjoint are employed to propagate the analysis increment and the gradient of the cost function over the 6-h assimilation window; a comprehensive set of simplified physical parameterizations is used during the final minimization process; and the number of frequently reported data, in particular satellite data, has substantially increased. The impact of these 4DVAR components on the forecast skill is reported in this article. This is achieved by comparing data assimilation configurations that range in complexity from the former 3DVAR with the implemented 4DVAR over a 1-month period. It is shown that the implementation of the tangent-linear model and its adjoint as well as the increased number of observations are the two features of the new 4DVAR that contribute the most to the forecast improvement. All the other components provide marginal though positive impact. 4DVAR does not improve the medium-range forecast of tropical storms in general and tends to amplify the existing, too early extratropical transition often observed in the MSC global forecast system with 3DVAR. It is shown that this recurrent problem is, however, more sensitive to the forecast model than the data assimilation scheme employed in this system. Finally, the impact of using a shorter cutoff time for the reception of observations, as the one used in the operational context for the 0000 and 1200 UTC forecasts, is more detrimental with 4DVAR. This result indicates that 4DVAR is more sensitive to observations at the end of the assimilation window than 3DVAR.

The Canadian Sea Ice and Snow Evolution (CanSISE) Network is a climate research network focused on developing and applying state-of-the-art observational data to advance dynamical prediction, projections, and understanding of seasonal snow cover and sea ice in Canada and the circumpolar Arctic. This study presents an assessment from the CanSISE Network of the ability of the second-generation Canadian Earth System Model (CanESM2) and the Canadian Seasonal to Interannual Prediction System (CanSIPS) to simulate and predict snow and sea ice from seasonal to multi-decadal timescales, with a focus on the Canadian sector. To account for observational uncertainty, model structural uncertainty, and internal climate variability, the analysis uses multi-source observations, multiple Earth system models (ESMs) in Phase5 of the Coupled Model Intercomparison Project (CMIP5), and large initial-condition ensembles of CanESM2 and other models. It is found that the ability of the CanESM2 simulation to capture snow-related climate parameters, such as cold-region surface temperature and precipitation, lies within the range of currently available international models. Accounting for the considerable disagreement among satellite-era observational datasets on the distribution of snow water equivalent, CanESM2 has too much springtime snow mass over Canada, reflecting a broader northern hemispheric positive bias. Biases in seasonal snow cover extent are generally less pronounced. CanESM2 also exhibits retreat of springtime snow generally greater than observational estimates, after accounting for observational uncertainty and internal variability. Sea ice is biased low in the Canadian Arctic, which makes it difficult to assess the realism of long-term sea ice trends there. The strengths and weaknesses of the modelling system need to be understood as a practical tradeoff: the Canadian models are relatively inexpensive computationally because of their moderate resolution, thus enabling their use in operational seasonal prediction and for generating large ensembles of multidecadal simulations. Improvements in climate-prediction systems like CanSIPS rely not just on simulation quality but also on using novel observational constraints and the ready transfer of research to an operational setting. Improvements in seasonal forecasting practice arising from recent research include accurate initialization of snow and frozen soil, accounting for observational uncertainty in forecast verification, and sea ice thickness initialization using statistical predictors available in real time.

Abstract. Glacier mass balance models are needed at sites with scarce long-term observations to reconstruct past glacier mass balance and assess its sensitivity to future climate change. In this study, North American Regional Reanalysis (NARR) data were used to force a physically based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess its sensitivity to climate change. A 2-year record (2014–2016) from an on-glacier automatic weather station (AWS) and historical precipitation records from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the NARR to the station sites. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multitemporal digital elevation model dataset. The model showed a good performance against recent (2012–2016) direct glaciological mass balance observations as well as with cumulative geodetic mass balance estimates. The simulated mass balance was not very sensitive to the NARR spatial interpolation method, as long as station data were used for bias correction. The simulated mass balance was however sensitive to the biases in NARR precipitation and air temperature, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these two parameters with ancillary data. The glacier-wide simulated energy balance regime showed a large contribution (57 %) of turbulent (sensible and latent) heat fluxes to melting in summer, higher than typical mid-latitude glaciers in continental climates, which reflects the local humid “icefield weather” of the Columbia Icefield. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 and 7 ∘C and precipitation between −20 % and +20 %, which comprise the spread of ensemble Representative Concentration Pathway (RCP) climate scenarios for the mid (2041–2070) and late (2071–2100) 21st century. The climate sensitivity experiments showed that future changes in precipitation would have a small impact on glacier mass balance, while the temperature sensitivity increases with warming, from −0.65 to −0.93 m w.e. a−1 ∘C−1. The mass balance response to warming was driven by a positive albedo feedback (44 %), followed by direct atmospheric warming impacts (24 %), a positive air humidity feedback (22 %) and a positive precipitation phase feedback (10 %). Our study underlines the key role of albedo and air humidity in modulating the response of winter-accumulation type mountain glaciers and upland icefield-outlet glacier settings to climate.

This study evaluates predictive uncertainties in the snow hydrology of the Fraser River Basin(FRB) of British Columbia(BC), Canada, using the Variable Infiltration Capacity(VIC) model forced with several high-resolution gridded climate datasets. These datasets include the Canadian Precipitation Analysis and the thin-plate smoothing splines(ANUSPLIN), North American Regional Reanalysis(NARR), University of Washington(UW) and Pacific Climate Impacts Consortium(PCIC) gridded products. Uncertainties are evaluated at different stages of the VIC implementation, starting with the driving datasets, optimization of model parameters, and model calibration during cool and warm phases of the Pacific Decadal Oscillation(PDO). The inter-comparison of the forcing datasets (precipitation and air temperature) and their VIC simulations (snow water equivalent – SWE – and runoff) reveals widespread differences over the FRB, especially in mountainous regions. The ANUSPLIN precipitation shows a considerable dry bias in the Rocky Mountains, whereas the NARR winter air temperature is 2°C warmer than the other datasets over most of the FRB. In the VIC simulations, the elevation-dependent changes in the maximum SWE(maxSWE) are more prominent at higher elevations of the Rocky Mountains, where the PCIC-VIC simulation accumulates too much SWE and ANUSPLIN-VIC yields an underestimation. Additionally, at each elevation range, the day of maxSWE varies from 10to 20days between the VIC simulations. The snow melting season begins early in the NARR-VIC simulation, whereas the PCIC-VIC simulation delays the melting, indicating seasonal uncertainty in SWE simulations. When compared with the observed runoff for the Fraser River main stem at Hope, BC, the ANUSPLIN-VIC simulation shows considerable underestimation of runoff throughout the water year owing to reduced precipitation in the ANUSPLIN forcing dataset. The NARR-VIC simulation yields more winter and spring runoff and earlier decline of flows in summer due to a nearly 15-day earlier onset of the FRB springtime snowmelt. Analysis of the parametric uncertainty in the VIC calibration process shows that the choice of the initial parameter range plays a crucial role in defining the model hydrological response for the FRB. Furthermore, the VIC calibration process is biased toward cool and warm phases of the PDO and the choice of proper calibration and validation time periods is important for the experimental setup. Overall the VIC hydrological response is prominently influenced by the uncertainties involved in the forcing datasets rather than those in its parameter optimization and experimental setups.

ABSTRACTStatistical relationships between weather conditions and the release of snow avalanches in the low-elevation coastal valleys of the northern Gaspe Peninsula are still poorly validated. As s...

The Appalachian Mountains of Eastern Canada are prone to several mass-wasting processes related to the geology and the nearby presence of large water bodies that influence the climate. Superimposed on this rugged terrain is the impacts of ongoing climate change, which may increase the magnitude, frequency, and duration of an array of hillslope phenomena. In this regard, the quantification of sediment fluxes at various spatiotemporal scales is prerequisite to reducing the exposure of infrastructure and communities, as well as to better understanding the mountain landscape evolution. Here, we report the quantitative modeling of sediment fluxes of several hillslope processes, mainly based on radiocarbon dating, which in turn improves understanding of how sediment has been eroded and transported through these mountain catchments since deglaciation. The results show a variable pattern of paraglacial effects at local and regional scales, highlighting the importance of ecological and hydroclimatic conditions in controlling the duration of glacially conditioned sedimentary stock exhaustion, and therefore the delay of paraglacial responses by geomorphic land systems. Current active scree slopes under the cold-temperate climate are characterized by sedimentation rates slightly lower than those calculated for the periglacial period following deglaciation, and even the sporadic remobilization of the primary stock by alluvial fan dynamics appears to be significant, testifying to a duration of paraglacial processes of more than 10,000 years.

On 15 March 2005, the Meteorological Service of Canada (MSC) proceeded to the implementation of a four-dimensional variational data assimilation (4DVAR) system, which led to significant improvements in the quality of global forecasts. This paper describes the different elements of MSC’s 4DVAR assimilation system, discusses some issues encountered during the development, and reports on the overall results from the 4DVAR implementation tests. The 4DVAR system adopted an incremental approach with two outer iterations. The simplified model used in the minimization has a horizontal resolution of 170 km and its simplified physics includes vertical diffusion, surface drag, orographic blocking, stratiform condensation, and convection. One important element of the design is its modularity, which has permitted continued progress on the three-dimensional variational data assimilation (3DVAR) component (e.g., addition of new observation types) and the model (e.g., computational and numerical changes). This paper discusses some numerical problems that occur in the vicinity of the Poles where the semi-Lagrangian scheme becomes unstable when there is a simultaneous occurrence of converging meridians and strong wind gradients. These could be removed by filtering the winds in the zonal direction before they are used to estimate the upstream position in the semi-Lagrangian scheme. The results show improvements in all aspects of the forecasts over all regions. The impact is particularly significant in the Southern Hemisphere where 4DVAR is able to extract more information from satellite data. In the Northern Hemisphere, 4DVAR accepts more asynoptic data, in particular coming from profilers and aircrafts. The impact noted is also positive and the short-term forecasts are particularly improved over the west coast of North America. Finally, the dynamical consistency of the 4DVAR global analyses leads to a significant impact on regional forecasts. Experimentation has shown that regional forecasts initiated directly from a 4DVAR global analysis are improved with respect to the regional forecasts resulting from the regional 3DVAR analysis.

Abstract River ice breakup has extensive implications on cold-region hydrological, ecological and river morphological systems. However, spatial and temporal breakup patterns under the changing climate are not well explored on large scale. This study discusses the spatial-temporal variations of breakup timing over terrestrial ecozones and five selected river basins of Canada based on long-term (1950–2016) data record. The link between the discovered patterns and climatic drivers (including air temperature, snowfall and rainfall), as well as elevation and anthropogenic activities are analyzed. An overall earlier breakup trend is observed across Canada and the spring air temperature is found to be the main driver behind it. However, the most pronounced warming trends across Canada is observed in winter. Spring warming trend is not as strong as winter warming and even becomes weak as period changes from 1950–2016 to 1970–2016, resulting in more stations showing later and significant later breakup during 1970–2016. Breakup pattern also displays evident spatial differences. Significant earlier breakup trends are mainly seen in western Canada (e.g. the Nelson River basin) and Arctic where spring warming trends are evident. Later and mixed breakup trends are generally identified in regions with weak warming or even cooling trends, such as Atlantic Canada and the St. Lawrence River basin. Spring snowfall generally delays breakup. Spring rainfall usually advances breakup dates while winter-rainfall can also delay breakup through refreezing. The increased snowfall in the north and increased rainfall in the south may be the reason why breakup timing is more sensitive to climatic warming in lower latitude regions than in higher latitude regions. Additionally, breakup timing in main streams and large rivers appears to be less sensitive to the warming trend than the headwaters and small tributaries. Elevation and flow regulation are also found to be contributing factors to the changes in breakup timing.

Abstract The database of the Quebec Ministry of Transport allowed us to analyze the occurrence of ice-block falls and snow avalanches for the past decades along national road 132. The results show that ice structure collapse may be categorized into three distinct phases by using daily temperatures (minimum, maximum, and average) and the cumulative degree day (temperatures above 0°C) since the March 1 st , corresponding to the beginning of the ice wall melting period: 1) a short and intense period of ice-block falls from the mid-April to the beginning of May; 2) a period of constant activity, mainly during the two first weeks of May; and 3) isolated residual activity, with a low frequency of ice-block falls until the month of June. The snow avalanche days were mainly characterized by significant snowfalls or rain-on-snow events with temperature>0°C. The multi-hazard probability was then evaluated based on the timing and relative frequency of ice-block fall and the modeling of sufficient snowpack for avalanching. This simple method to assess the synergistic effect of hillslope processes allows a better understanding of the spring avalanche regime related to the collapse of ice structures. These findings are expected to assist in the management of natural hazards and to improve our knowledge of spatiotemporal dynamics of mass-wasting events on highways.