Climate and climate change in western canada as simulated by the Canadian regional climate model

Abstract

A þrst climate simulation performed with the novel Canadian Regional Climate Model (CRCM) is presented. The CRCM is based on fully elastic non‐hydrostatic þeld equations, which are solved with an efþcient semi‐implicit semi‐Lagrangian (SISL) marching algorithm, and on the parametrization package of subgrid‐scale physical effects of the second‐generation Canadian Global Climate Model (GCMII). Two 5‐year integrations of the CRCM nested with GCMII simulated data as lateral boundary conditions are made for conditions corresponding to current and doubled CO₂ scenarios. For these simulations the CRCM used a grid size of 45 km on a polar‐stereographic projection, 20 scaled‐height levels and a time step of 15 min; the nesting GCMII has a spectral truncation of T32, 10 hybrid‐pressure levels and a time step of 20 min. These simulations serve to document: (1) the suitability of the SISL numerical scheme for regional climate modelling, (2) the use of GCMII physics at much higher resolution than in the nesting model, (3) the ability of the CRCM to add realistic regional‐scale climate information to global model simulations, and (4) the climate of the CRCM compared to that of GCMII under two greenhouse gases (GHG) scenarios.     

Application of the Canadian regional climate model to the Laurentian great lakes region: Implementation of a lake model

Abstract

This study reports on the implementation of an interactive mixed‐layer/thermodynamic‐ice lake model coupled with the Canadian Regional Climate Model (CRCM). For this application the CRCM, which uses a grid mesh of 45 km on a polar stereographic projection, 10 vertical levels, and a timestep of 15 min, is nested with the second generation Canadian General Circulation Model (GCM) simulated output. A numerical simulation of the climate of eastern North America, including the Laurentian Great Lakes, is then performed in order to evaluate the coupled model. The lakes are represented by a “mixed layer” model to simulate the evolution of the surface water temperature, and a thermodynamic ice model to simulate evolution of the ice cover. The mixed‐layer depth is allowed to vary spatially. Lake‐ice leads are parametrized as a function of ice thickness based on observations. Results from a 5‐year integration show that the coupled CRCM/lake model is capable of simulating the seasonal evolution of surface temperature and ice cover in the Great Lakes. When compared with lake climatology, the simulated mean surface water temperature agrees within 0.12°C on average. The seasonal evolution of the lake‐ice cover is realistic but the model tends to underestimate the monthly mean ice concentration on average. The simulated winter lake‐induced precipitation is also shown, and snow accumulation patterns on downwind shores of the lakes are found to be realistic when compared with observations.     

Modulation of soil moisture–precipitation interactions over France by large scale circulation

How soil moisture affects precipitation is an important question—with far reaching consequences, from weather prediction to centennial climate change—, albeit a poorly understood one. In this paper, an analysis of soil moisture–precipitation interactions over France based on observations is presented. A first objective of this paper is to investigate how large scale circulation modulates soil moisture–precipitation interactions, thanks to a weather regime approach. A second objective is to study the influence of soil moisture not only on precipitation but also on the difference between precipitation and evapotranspiration. Indeed, to have a total positive soil moisture–precipitation feedback, the potential decrease in precipitation associated with drier soils should be larger than the decrease in evapotranspiration that drier soils may also cause. A potential limited impact of soil moisture on precipitation is found for some weather regimes, but its sign depends on large scale circulation. Indeed, antecedent dry soil conditions tend to lead to smaller precipitation for the negative phase of the North Atlantic Oscillation (NAO) regime but to larger precipitation for the Atlantic Low regime. This differential response of precipitation to soil moisture anomalies depending on large scale circulation is traced back to different responses of atmospheric stability. For all circulation regimes, dry soils tend to increase the lifted condensation level, which is unfavorable to precipitation. But for the negative phase of the NAO, low soil moisture tends to lead to an increase of atmospheric stability while it tends to lead to a decrease of stability for Atlantic Low. Even if the impact of soil moisture anomalies varies depending on large scale circulation (it is larger for Atlantic low and the positive phase of the NAO), dry soils always lead to a decrease in evapotranspiration. As the absolute effect of antecedent soil moisture on evapotranspiration is always much larger than its effects on precipitation, for all circulation regimes dry soil anomalies subsequently lead to positive precipitation minus evapotranspiration anomalies i.e. the total soil moisture feedback is found to be negative. This negative feedback is stronger for the Atlantic Low and the positive phase of the NAO regimes.     

Nonparametric analysis of high wind speed data

In this paper, nonparametric curve estimation methods are applied to analyze time series of wind speeds, focusing on the extreme events exceeding a chosen threshold. Classical parametric statistical approaches in this context consist in fitting a generalized Pareto distribution (GPD) to the tail of the empirical cumulative distribution, using maximum likelihood or the method of the moments to estimate the parameters of this distribution. Additionally, confidence intervals are usually computed to assess the uncertainty of the estimates. Nonparametric methods to estimate directly some quantities of interest, such as the probability of exceedance, the quantiles or return levels, or the return periods, are proposed. Moreover, bootstrap techniques are used to develop pointwise and simultaneous confidence intervals for these functions. The proposed models are applied to wind speed data in the Gulf Coast of US, comparing the results with those using the GPD approach, by means of a split-sample test. Results show that nonparametric methods are competitive with respect to the standard GPD approximations. The study is completed generating synthetic data sets and comparing the behavior of the parametric and the nonparametric estimates in this framework.     

An approximate energy cycle for inter-member variability in ensemble simulations of a regional climate model

The presence of internal variability (IV) in ensembles of nested regional climate model (RCM) simulations is now widely acknowledged in the community working on dynamical downscaling. IV is defined as the inter-member spread between members in an ensemble of simulations performed by a given RCM driven by identical lateral boundary conditions (LBC), where different members are being initialised at different times. The physical mechanisms responsible for the time variations and structure of such IV have only recently begun to receive attention. Recent studies have shown empirical evidence of a close parallel between the energy conversions associated with the time fluctuations of IV in ensemble simulations of RCM and the energy conversions taking place in weather systems. Inspired by the classical work on global energetics of weather systems, we sought a formulation of an energy cycle for IV that would be applicable for limited-area domain. We develop here a novel formalism based on local energetics that can be applied to further our understanding IV. Prognostic equations for ensemble-mean kinetic energy and available enthalpy are decomposed into contributions due to ensemble-mean variables (EM) and those due to deviations from the ensemble mean (IV). Together these equations constitute an energy cycle for IV in ensemble simulations of RCM. Although the energy cycle for IV was developed in a context entirely different from that of energetics of weather systems, the exchange terms between the various reservoirs have a rather similar mathematical form, which facilitates some interpretations of their physical meaning.     

Wake Response to an Ocean-Feedback Mechanism: Madeira Island Case Study

We focus on an island wake episode that occurred in the Madeira Archipelago region of the north-east Atlantic at $32.5^{\circ }\mathrm{N}, 17^{\circ }\mathrm{W}$ . The Weather Research and Forecasting numerical model was used in a (one-way) downscaling mode, considering initial and boundary conditions from the European Centre for Medium-range Weather Forecasts system. The current literature emphasizes adiabatic effects on the dynamical aspects of atmospheric wakes. Changes in mountain height and consequently its relation to the atmospheric inversion layer should explain the shift in wake regimes, from a ‘strong-wake’ to ‘weak-wake’ scenario. Nevertheless, changes in sea-surface temperature variability in the lee of an island can induce similar regime shifts because of exposure to stronger solar radiation. Increase in evaporation contributes to the enhancement of convection and thus to the uplift of the stratified atmospheric layer above the critical height, with subsequent internal gravity wave activity.     

Evolution of the Parisian urban climate under a global changing climate

The evolution of the Parisian urban climate under a changing climate is analyzed from long-term offline numerical integrations including a specific urban parameterization. This system is forced by meteorological conditions based on present-climate reanalyses (1970–2007), and climate projections (2071–2099) provided by global climate model simulations following two emission scenarios (A1B and A2). This study aims at quantifying the impact of climate change on air temperature within the city and in the surroundings. A systematic increase of 2-meter air temperature is found. In average according to the two scenarios, it reaches +?2.0/2.4°C in winter and +?3.5/5.0°C in summer for the minimum and maximum daily temperatures, respectively. During summer, the warming trend is more pronounced in the surrounding countryside than in Paris and suburbs due to the soil dryness. As a result, a substantial decrease of the strong urban heat islands is noted at nighttime, and numerous events with negative urban heat islands appear at daytime. Finally, a 30% decrease of the heating degree days is quantified in winter between present and future climates. Inversely, the summertime cooling degree days significantly increase in future climate whereas they are negligible in present climate. However, in terms of accumulated degree days, the increase of the demand in cooling remains smaller than the decrease of the demand in heating.     

Impact of climate change on the hydrogeology of two basins in northern France

This study presents an analysis of climate-change impacts on the water resources of two basins located in northern France, by integrating four sources of uncertainty: climate modelling, hydrological modelling, downscaling methods, and emission scenarios. The analysis focused on the evolution of the water budget, the river discharges and piezometric heads. Seven hydrological models were used, from lumped rainfall-discharge to distributed hydrogeological models, and led to quite different estimates of the water-balance components. One of the hydrological models, CLSM, was found to be unable to simulate the increased water stress and was, thus, considered as an outlier even though it gave fair results for the present day compared to observations. Although there were large differences in the results between the models, there was a marked tendency towards a decrease of the water resource in the rivers and aquifers (on average in 2050 about ?14 % and ?2.5 m, respectively), associated with global warming and a reduction in annual precipitation (on average in 2050 +2.1 K and ?3 %, respectively). The uncertainty associated to climate models was shown to clearly dominate, while the three others were about the same order of magnitude and 3–4 times lower. In terms of impact, the results found in this work are rather different from those obtained in a previous study, even though two of the hydrological models and one of the climate models were used in both studies. This emphasizes the need for a survey of the climatic-change impact on the water resource.     

The Indo-Australian monsoon and its relationship to ENSO and IOD in reanalysis data and the CMIP3/CMIP5 simulations

A large spread exists in both Indian and Australian average monsoon rainfall and in their interannual variations diagnosed from various observational and reanalysis products. While the multi model mean monsoon rainfall from 59 models taking part in the Coupled Model Intercomparison Project (CMIP3 and CMIP5) fall within the observational uncertainty, considerable model spread exists. Rainfall seasonality is consistent across observations and reanalyses, but most CMIP models produce either a too peaked or a too flat seasonal cycle, with CMIP5 models generally performing better than CMIP3. Considering all North-Australia rainfall, most models reproduce the observed Australian monsoon-El Niño Southern Oscillation (ENSO) teleconnection, with the strength of the relationship dependent on the strength of the simulated ENSO. However, over the Maritime Continent, the simulated monsoon-ENSO connection is generally weaker than observed, depending on the ability of each model to realistically reproduce the ENSO signature in the Warm Pool region. A large part of this bias comes from the contribution of Papua, where moisture convergence seems to be particularly affected by this SST bias. The Indian summer monsoon-ENSO relationship is affected by overly persistent ENSO events in many CMIP models. Despite significant wind anomalies in the Indian Ocean related to Indian Ocean Dipole (IOD) events, the monsoon-IOD relationship remains relatively weak both in the observations and in the CMIP models. Based on model fidelity in reproducing realistic monsoon characteristics and ENSO teleconnections, we objectively select 12 “best” models to analyze projections in the rcp8.5 scenario. Eleven of these models are from the CMIP5 ensemble. In India and Australia, most of these models produce 5–20 % more monsoon rainfall over the second half of the twentieth century than during the late nineteenth century. By contrast, there is no clear model consensus over the Maritime Continent.