Contrasting responses of mean and extreme snowfall to climate change

Snowfall is an important element of the climate system, and one that is expected to change in a warming climate. Both mean snowfall and the intensity distribution of snowfall are important, with heavy snowfall events having particularly large economic and human impacts. Simulations with climate models indicate that annual mean snowfall declines with warming in most regions but increases in regions with very low surface temperatures. The response of heavy snowfall events to a changing climate, however, is unclear. Here I show that in simulations with climate models under a scenario of high emissions of greenhouse gases, by the late twenty-first century there are smaller fractional changes in the intensities of daily snowfall extremes than in mean snowfall over many Northern Hemisphere land regions. For example, for monthly climatological temperatures just below freezing and surface elevations below 1,000 metres, the 99.99th percentile of daily snowfall decreases by 8% in the multimodel median, compared to a 65% reduction in mean snowfall. Both mean and extreme snowfall must decrease for a sufficiently large warming, but the climatological temperature above which snowfall extremes decrease with warming in the simulations is as high as −9 °C, compared to −14 °C for mean snowfall. These results are supported by a physically based theory that is consistent with the observed rain–snow transition. According to the theory, snowfall extremes occur near an optimal temperature that is insensitive to climate warming, and this results in smaller fractional changes for higher percentiles of daily snowfall. The simulated changes in snowfall that I find would influence surface snow and its hazards; these changes also suggest that it may be difficult to detect a regional climate-change signal in snowfall extremes.National Science Foundation (U.S.) (Grant AGS-1148594)United States. National Aeronautics and Space Administration (ROSES Grant 09-IDS09-0049

Very extreme seasonal precipitation in the NARCCAP ensemble: model performance and projections

Seasonal extreme daily precipitation is analyzed in the ensemble of NARCAPP regional climate models. Significant variation in these models' abilities to reproduce observed precipitation extremes over the contiguous United States is found. Model performance metrics are introduced to characterize overall biases, seasonality, spatial extent and the shape of the precipitation distribution. Comparison of the models to gridded observations that include an elevation correction is found to be better than to gridded observations without this correction. A complicated model weighting scheme based on model performance in simulating observations is found to cause significant improvements in ensemble mean skill only if some of the models are poorly performing outliers. The effect of lateral boundary conditions are explored by comparing the integrations driven by reanalysis to those driven by global climate models. Projected mid-century future changes in seasonal precipitation means and extremes are presented and discussions of the sources of uncertainty and the mechanisms causing these changes are presented. © 2012 The Author(s)

Bias adjustment and ensemble recalibration methods for seasonal forecasting: a comprehensive intercomparison using the C3S dataset

This work presents a comprehensive intercomparison of diferent alternatives for the calibration of seasonal forecasts, ranging from simple bias adjustment (BA)-e.g. quantile mapping-to more sophisticated ensemble recalibration (RC) methods- e.g. non-homogeneous Gaussian regression, which build on the temporal correspondence between the climate model and the corresponding observations to generate reliable predictions. To be as critical as possible, we validate the raw model and the calibrated forecasts in terms of a number of metrics which take into account diferent aspects of forecast quality (association, accuracy, discrimination and reliability). We focus on one-month lead forecasts of precipitation and temperature from four state-of-the-art seasonal forecasting systems, three of them included in the Copernicus Climate Change Service dataset (ECMWF-SEAS5, UK Met Ofce-GloSea5 and Météo France-System5) for boreal winter and summer over two illustrative regions with diferent skill characteristics (Europe and Southeast Asia). Our results indicate that both BA and RC methods efectively correct the large raw model biases, which is of paramount importance for users, particularly when directly using the climate model outputs to run impact models, or when computing climate indices depending on absolute values/thresholds. However, except for particular regions and/or seasons (typically with high skill), there is only marginal added value-with respect to the raw model outputs-beyond this bias removal. For those cases, RC methods can outperform BA ones, mostly due to an improvement in reliability. Finally, we also show that whereas an increase in the number of members only modestly afects the results obtained from calibration, longer hindcast periods lead to improved forecast quality, particularly for RC methods.This work has been funded by the C3S activity on Evaluation and Quality Control for seasonal forecasts. JMG was partially supported by the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER). FJDR was partially funded by the H2020 EUCP project (GA 776613)

How to create an operational multi-model of seasonal forecasts?

Seasonal forecasts of variables like near-surface temperature or precipitation are becoming increasingly important for a wide range of stakeholders. Due to the many possibilities of recalibrating, combining, and verifying ensemble forecasts, there are ambiguities of which methods are most suitable. To address this we compare approaches how to process and verify multi-model seasonal forecasts based on a scientific assessment performed within the framework of the EU Copernicus Climate Change Service (C3S) Quality Assurance for Multi-model Seasonal Forecast Products (QA4Seas) contract C3S 51 lot 3. Our results underpin the importance of processing raw ensemble forecasts differently depending on the final forecast product needed. While ensemble forecasts benefit a lot from bias correction using climate conserving recalibration, this is not the case for the intrinsically bias adjusted multi-category probability forecasts. The same applies for multi-model combination. In this paper, we apply simple, but effective, approaches for multi-model combination of both forecast formats. Further, based on existing literature we recommend to use proper scoring rules like a sample version of the continuous ranked probability score and the ranked probability score for the verification of ensemble forecasts and multi-category probability forecasts, respectively. For a detailed global visualization of calibration as well as bias and dispersion errors, using the Chi-square decomposition of rank histograms proved to be appropriate for the analysis performed within QA4Seas.The research leading to these results is part of the Copernicus Climate Change Service (C3S) (Framework Agreement number C3S_51_Lot3_BSC), a program being implemented by the European Centre for Medium-Range Weather Forecasts (ECMWF) on behalf of the European Commission. Francisco Doblas-Reyes acknowledges the support by the H2020 EUCP project (GA 776613) and the MINECO-funded CLINSA project (CGL2017-85791-R)

Decadal predictions with the HiGEM high resolution global coupled climate model: description and basic evaluation

This paper describes the development and basic evaluation of decadal predictions produced using the HiGEM coupled climate model. HiGEM is a higher resolution version of the HadGEM1 Met Office Unified Model. The horizontal resolution in HiGEM has been increased to 1.25◦ × 0.83◦ in longitude and latitude for the atmosphere, and 1/3◦ × 1/3◦ globally for the ocean. The HiGEM decadal predictions are initialised using an anomaly assimilation scheme that relaxes anomalies of ocean temperature and salinity to observed anomalies. 10 year hindcasts are produced for 10 start dates (1960, 1965,..., 2000, 2005). To determine the relative contributions to prediction skill from initial conditions and external forcing, the HiGEM decadal predictions are compared to uninitialised HiGEM transient experiments. The HiGEM decadal predictions have substantial skill for predictions of annual mean surface air temperature and 100 m upper ocean temperature. For lead times up to 10 years, anomaly correlations (ACC) over large areas of the North Atlantic Ocean, the Western Pacific Ocean and the Indian Ocean exceed values of 0.6. Initialisation of the HiGEM decadal predictions significantly increases skill over regions of the Atlantic Ocean,the Maritime Continent and regions of the subtropical North and South Pacific Ocean. In particular, HiGEM produces skillful predictions of the North Atlantic subpolar gyre for up to 4 years lead time (with ACC > 0.7), which are significantly larger than the uninitialised HiGEM transient experiments

Analysis of rainfall seasonality from observations and climate models

Two new indicators of rainfall seasonality based on information entropy, the relative entropy (RE) and the dimensionless seasonality index (DSI), together with the mean annual rainfall, are evaluated on a global scale for recently updated precipitation gridded datasets and for historical simulations from coupled atmosphere--ocean general circulation models. The RE provides a measure of the number of wet months and, for precipitation regimes featuring a distinct wet and dry season, it is directly related to the duration of the wet season. The DSI combines the rainfall intensity with its degree of seasonality and it is an indicator of the extent of the global monsoon region. We show that the RE and the DSI are fairly independent of the time resolution of the precipitation data, thereby allowing objective metrics for model intercomparison and ranking. Regions with different precipitation regimes are classified and characterized in terms of RE and DSI. Comparison of different land observational datasets reveals substantial difference in their local representation of seasonality. It is shown that two-dimensional maps of RE provide an easy way to compare rainfall seasonality from various datasets and to determine areas of interest. Models participating to the Coupled Model Intercomparison Project platform, Phase 5, consistently overestimate the RE over tropical Latin America and underestimate it in West Africa, western Mexico and East Asia. It is demonstrated that positive RE biases in a general circulation model are associated with excessively peaked monthly precipitation fractions, too large during the wet months and too small in the months preceding and following the wet season; negative biases are instead due, in most cases, to an excess of rainfall during the premonsoonal months