A novel method to improve temperature simulations of general circulation models based on ensemble empirical mode decomposition and its application to multi-model ensembles

A novel method based on the ensemble empirical mode decomposition (EEMD) method was developed to improve model performance. This method was evaluated by applying it to global surface air temperatures, which were simulated by eight general circulation models from the Coupled Model Intercomparison Project Phase 5 (CMIP5). The temperature simulations of the eight models were separated into their different components by EEMD. The model's performance improved after the first high-frequency component was removed from the original simulations by EEMD for each model, on both the global and continental scale. Moreover, EEMD was more effective in improving the model's performance compared to the wavelet transform method. The multi-model ensembles (MMEs) were calculated based on the EEMD-improved model simulations using the Average Ensemble Mean, Multiple Linear Regression, Singular Value Decomposition and Bayesian Model Averaging methods. The results showed that the MME forecasts performed better when the calculations were based on the EEMD-improved simulations as opposed to the original simulations on both the global and continental scale. Therefore, the results of the MME were further improved by using the EEMD-improved model simulations. This new method provides a simple way to improve model performance and can be easily applied to further improve MME simulations

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)

Process-conditioned bias correction for seasonal forecasting: a case-study with ENSO in Peru

This work assesses the suitability of a first simple attempt for process-conditioned bias correction in the context of seasonal forecasting. To do this, we focus on the northwestern part of Peru and bias correct 1- and 4-month lead seasonal predictions of boreal winter (DJF) precipitation from the ECMWF System4 forecasting system for the period 1981–2010. In order to include information about the underlying large-scale circulation which may help to discriminate between precipitation affected by different processes, we introduce here an empirical quantile–quantile mapping method which runs conditioned on the state of the Southern Oscillation Index (SOI), which is accurately predicted by System4 and is known to affect the local climate. Beyond the reduction of model biases, our results show that the SOI-conditioned method yields better ROC skill scores and reliability than the raw model output over the entire region of study, whereas the standard unconditioned implementation provides no added value for any of these metrics. This suggests that conditioning the bias correction on simple but well-simulated large-scale processes relevant to the local climate may be a suitable approach for seasonal forecasting. Yet, further research on the suitability of the application of similar approaches to the one considered here for other regions, seasons and/or variables is needed.This work has received funding from the MULTI-SDM project (MINECO/FEDER, CGL2015-66583-R). The authors are grateful to SENAMHI for the observational data, which are publicly available from http://www.senamhi.gob.pe/?p=data-historica, and to the European Center for Medium-Range Weather Forecast (ECMWF), for the access to the System4 seasonal forecasting hindcast

Identifying causes of Western Pacific ITCZ drift in ECMWF System 4 hindcasts

The development of systematic biases in climate models used in operational seasonal forecasting adversely affects the quality of forecasts they produce. In this study, we examine the initial evolution of systematic biases in the ECMWF System 4 forecast model, and isolate aspects of the model simulations that lead to the development of these biases. We focus on the tendency of the simulated intertropical convergence zone in the western equatorial Pacific to drift northwards by between 0.5° and 3° of latitude depending on season. Comparing observations with both fully coupled atmosphere–ocean hindcasts and atmosphere-only hindcasts (driven by observed sea-surface temperatures), we show that the northward drift is caused by a cooling of the sea-surface temperature on the Equator. The cooling is associated with anomalous easterly wind stress and excessive evaporation during the first twenty days of hindcast, both of which occur whether air-sea interactions are permitted or not. The easterly wind bias develops immediately after initialisation throughout the lower troposphere; a westerly bias develops in the upper troposphere after about ten days of hindcast. At this point, the baroclinic structure of the wind bias suggests coupling with errors in convective heating, although the initial wind bias is barotropic in structure and appears to have an alternative origin

Predicting the seasonal evolution of southern African summer precipitation in the DePreSys3 prediction system

We assess the ability of the DePreSys3 prediction system to predict austral summer precipitation (DJF) over southern Africa, defined as the African continent south of 15°S. DePresys3 is a high resolution prediction system (at a horizontal resolution of ~ 60 km in the atmosphere in mid-latitudes and of the quarter degree in the Ocean) and spans the long period 1959–2016. We find skill in predicting interannual precipitation variability, relative to a long-term trend; the anomaly correlation skill score over southern Africa is greater than 0.45 for the first summer (i.e. lead month 2–4), and 0.37 over Mozambique, Zimbabwe and Zambia for the second summer (i.e. lead month 14–16). The skill is related to the successful prediction of the El-Nino Southern Oscillation (ENSO), and the successful simulation of ENSO teleconnections to southern Africa. However, overall skill is sensitive to the inclusion of strong La-Nina events and also appears to change with forecast epoch. For example, the skill in predicting precipitation over Mozambique is significantly larger for the first summer in the 1990–2016 period, compared to the 1959–1985 period. The difference in skill in predicting interannual precipitation variability over southern Africa in different epochs is consistent with a change in the strength of the observed teleconnections of ENSO. After 1990, and consistent with the increased skill, the observed impact of ENSO appears to strengthen over west Mozambique, in association with changes in ENSO related atmospheric convergence anomalies. However, these apparent changes in teleconnections are not captured by the ensemble-mean predictions using DePreSys3. The changes in the ENSO teleconnection are consistent with a warming over the Indian Ocean and modulation of ENSO properties between the different epochs, but may also be associated with unpredictable atmospheric variability

Mediterranean winter rainfall in phase with African monsoons during the past 1.36 million years

Mediterranean climates are characterized by strong seasonal contrasts between dry summers and wet winters. Changes in winter rainfall are critical for regional socioeconomic development, but are difficult to simulate accurately1 and reconstruct on Quaternary timescales. This is partly because regional hydroclimate records that cover multiple glacial–interglacial cycles2,3 with different orbital geometries, global ice volume and atmospheric greenhouse gas concentrations are scarce. Moreover, the underlying mechanisms of change and their persistence remain unexplored. Here we show that, over the past 1.36 million years, wet winters in the northcentral Mediterranean tend to occur with high contrasts in local, seasonal insolation and a vigorous African summer monsoon. Our proxy time series from Lake Ohrid on the Balkan Peninsula, together with a 784,000-year transient climate model hindcast, suggest that increased sea surface temperatures amplify local cyclone development and refuel North Atlantic low-pressure systems that enter the Mediterranean during phases of low continental ice volume and high concentrations of atmospheric greenhouse gases. A comparison with modern reanalysis data shows that current drivers of the amount of rainfall in the Mediterranean share some similarities to those that drive the reconstructed increases in precipitation. Our data cover multiple insolation maxima and are therefore an important benchmark for testing climate model performance

Can bias correction and statistical downscaling methods improve the skill of seasonal precipitation forecasts?

Statistical downscaling methods are popular post-processing tools which are widely used in many sectors to adapt the coarse-resolution biased outputs from global climate simulations to the regional-to-local scale typically required by users. They range from simple and pragmatic Bias Correction (BC) methods, which directly adjust the model outputs of interest (e.g. precipitation) according to the available local observations, to more complex Perfect Prognosis (PP) ones, which indirectly derive local predictions (e.g. precipitation) from appropriate upper-air large-scale model variables (predictors). Statistical downscaling methods have been extensively used and critically assessed in climate change applications; however, their advantages and limitations in seasonal forecasting are not well understood yet. In particular, a key problem in this context is whether they serve to improve the forecast quality/skill of raw model outputs beyond the adjustment of their systematic biases. In this paper we analyze this issue by applying two state-of-the-art BC and two PP methods to downscale precipitation from a multimodel seasonal hindcast in a challenging tropical region, the Philippines. To properly assess the potential added value beyond the reduction of model biases, we consider two validation scores which are not sensitive to changes in the mean (correlation and reliability categories). Our results show that, whereas BC methods maintain or worsen the skill of the raw model forecasts, PP methods can yield significant skill improvement (worsening) in cases for which the large-scale predictor variables considered are better (worse) predicted by the model than precipitation. For instance, PP methods are found to increase (decrease) model reliability in nearly 40% of the stations considered in boreal summer (autumn). Therefore, the choice of a convenient downscaling approach (either BC or PP) depends on the region and the season.This study was partially supported by the SPECS and EUPORIAS projects, funded by the European Commission through the Seventh Framework Programme for Research under grant agreements 308378 and 308291, respectively. JMG acknowledges partial support from the project MULTI-SDM (CGL2015-66583-R, MINECO/FEDER)

Current and emerging developments in subseasonal to decadal prediction

Weather and climate variations of subseasonal to decadal timescales can have enormous social, economic and environmental impacts, making skillful predictions on these timescales a valuable tool for decision makers. As such, there is a growing interest in the scientific, operational and applications communities in developing forecasts to improve our foreknowledge of extreme events. On subseasonal to seasonal (S2S) timescales, these include high-impact meteorological events such as tropical cyclones, extratropical storms, floods, droughts, and heat and cold waves. On seasonal to decadal (S2D) timescales, while the focus remains broadly similar (e.g., on precipitation, surface and upper ocean temperatures and their effects on the probabilities of high-impact meteorological events), understanding the roles of internal and externally-forced variability such as anthropogenic warming in forecasts also becomes important. The S2S and S2D communities share common scientific and technical challenges. These include forecast initialization and ensemble generation; initialization shock and drift; understanding the onset of model systematic errors; bias correct, calibration and forecast quality assessment; model resolution; atmosphere-ocean coupling; sources and expectations for predictability; and linking research, operational forecasting, and end user needs. In September 2018 a coordinated pair of international conferences, framed by the above challenges, was organized jointly by the World Climate Research Programme (WCRP) and the World Weather Research Prograame (WWRP). These conferences surveyed the state of S2S and S2D prediction, ongoing research, and future needs, providing an ideal basis for synthesizing current and emerging developments in these areas that promise to enhance future operational services. This article provides such a synthesis