Regional climate models (RCMs) are important tools used for downscaling climate simulations from global scale models. In project CECILIA, two RCMs were used to provide climate change information for regions of Central and Eastern Europe. Models RegCM and ALADIN-Climate were employed in downscaling global simulations from ECHAM5 and ARPEGE-CLIMAT under IPCC A1B emission scenario in periods 2021-2050 and 2071-2100. Climate change signal present in these simulations is consistent with respective driving data, showing similar large-scale features: warming between 0 and 3 degrees C in the first period and 2 and 5 degrees C in the second period with the least warming in northwestern part of the domain increasing in the southeastern direction and small precipitation changes within range of +1 to -1 mm/day. Regional features are amplified by the RCMs, more so in case of the ALADIN family of models.

A multi-model ensemble of decadal prediction experiments, performed in the framework of the EU-funded COMBINE (Comprehensive Modelling of the Earth System for Better Climate Prediction and Projection) Project following the 5th Coupled Model Intercomparison Project protocol is examined. The ensemble combines a variety of dynamical models, initialization and perturbation strategies, as well as data assimilation products employed to constrain the initial state of the system. Taking advantage of the multi-model approach, several aspects of decadal climate predictions are assessed, including predictive skill, impact of the initialization strategy and the level of uncertainty characterizing the predicted fluctuations of key climate variables. The present analysis adds to the growing evidence that the current generation of climate models adequately initialized have significant skill in predicting years ahead not only the anthropogenic warming but also part of the internal variability of the climate system. An important finding is that the multi-model ensemble mean does generally outperform the individual forecasts, a well-documented result for seasonal forecasting, supporting the need to extend the multi-model framework to real-time decadal predictions in order to maximize the predictive capabilities of currently available decadal forecast systems. The multi-model perspective did also allow a more robust assessment of the impact of the initialization strategy on the quality of decadal predictions, providing hints of an improved forecast skill under full-value (with respect to anomaly) initialization in the near-term range, over the Indo-Pacific equatorial region. Finally, the consistency across the different model predictions was assessed. Specifically, different systems reveal a general agreement in predicting the near-term evolution of surface temperatures, displaying positive correlations between different decadal hindcasts over most of the global domain.

A case of daytime development of deep convection over tropical semi-arid land is used to evaluate the representation of convection in global and regional models. The case is based on observations collected during the African Monsoon Multidisciplinary Analysis (AMMA) field campaign and includes two distinct transition phases, from clear sky to shallow cumulus and from cumulus to deep convection. Different types of models, run with identical initial and boundary conditions, are intercompared: a reference large-eddy simulation (LES), single-column model (SCM) version of four different Earth system models that participated in the Coupled Model Intercomparison Project 5 exercise, the SCM version of the European Centre for Medium-range Weather Forecasts operational forecast model, the SCM version of a mesoscale model and a bulk model. Surface fluxes and radiative heating are prescribed preventing any atmosphere-surface and cloud-radiation coupling in order to simplify the analyses so that it focuses only on convective processes. New physics packages are also evaluated within this framework. As the LES correctly reproduces the observed growth of the boundary layer, the gradual development of shallow clouds, the initiation of deep convection and the development of cold pools, it provides a basis to evaluate in detail the representation of the diurnal cycle of convection by the other models and to test the hypotheses underlying convective parametrizations. Most SCMs have difficulty in representing the timing of convective initiation and rain intensity, although substantial modifications to boundary-layer and deep-convection parametrizations lead to improvements. The SCMs also fail to represent the mid-level troposphere moistening during the shallow convection phase, which we analyse further. Nevertheless, beyond differences in timing of deep convection, the SCM models reproduce the sensitivity to initial and boundary conditions simulated in the LES regarding boundary-layer characteristics, and often the timing of convection triggering.

The main purpose of this study is to investigate the influence of the Arctic Oscillation (AO), the dominant mode of natural variability over the northerly high latitudes, on the spatial (horizontal and vertical) distribution of clouds in the Arctic. To that end, we use a suite of sensors on-board NASA's A-Train satellites that provide accurate observations of the distribution of clouds along with information on atmospheric thermodynamics. Data from three independent sensors are used (AQUA-AIRS, CALIOP-CALIPSO and CPR-CloudSat) covering two time periods (winter half years, November through March, of 2002-2011 and 2006-2011, respectively) along with data from the ERA-Interim reanalysis. We show that the zonal vertical distribution of cloud fraction anomalies averaged over 67-82 degrees N to a first approximation follows a dipole structure (referred to as "Greenland cloud dipole anomaly", GCDA), such that during the positive phase of the AO, positive and negative cloud anomalies are observed eastwards and westward of Greenland respectively, while the opposite is true for the negative phase of AO. By investigating the concurrent meteorological conditions (temperature, humidity and winds), we show that differences in the meridional energy and moisture transport during the positive and negative phases of the AO and the associated thermodynamics are responsible for the conditions that are conducive for the formation of this dipole structure. All three satellite sensors broadly observe this large-scale GCDA despite differences in their sensitivities, spatio-temporal and vertical resolutions, and the available lengths of data records, indicating the robustness of the results. The present study also provides a compelling case to carry out process-based evaluation of global and regional climate models.

Multiyear climate predictions with two initialization strategies are systematically assessed in the EC-Earth V2.3 climate model. In one ensemble, an estimate of the observed climate state is used to initialize the model. The other uses estimates of observed ocean and sea ice anomalies on top of the model climatology. The ensembles show similar spatial characteristics of drift related to the biases in control simulations. As expected, the drift is less with anomaly initialization. The full field initialization overshoots to a colder state which is related to cold biases in the tropics and North Atlantic, associated with oceanic processes. Despite different amplitude of the drift, both ensembles show similar skill in multiyear global temperature predictions, but regionally differences are found. On multiyear time scales, initialization with observations enhances both deterministic and probabilistic skill scores in the North Atlantic. The probabilistic verification shows skill over the European continent.

Many theories for the Madden-Julian oscillation (MJO) focus on diabatic processes, particularly the evolution of vertical heating and moistening. Poor MJO performance in weather and climate models is often blamed on biases in these processes and their interactions with the large-scale circulation. We introduce one of the three components of a model evaluation project, which aims to connect MJO fidelity in models to their representations of several physical processes, focusing on diabatic heating and moistening. This component consists of 20day hindcasts, initialized daily during two MJO events in winter 2009-2010. The 13 models exhibit a range of skill: several have accurate forecasts to 20days lead, while others perform similarly to statistical models (8-11days). Models that maintain the observed MJO amplitude accurately predict propagation, but not vice versa. We find no link between hindcast fidelity and the precipitation-moisture relationship, in contrast to other recent studies. There is also no relationship between models' performance and the evolution of their diabatic heating profiles with rain rate. A more robust association emerges between models' fidelity and net moistening: the highest-skill models show a clear transition from low-level moistening for light rainfall to midlevel moistening at moderate rainfall and upper level moistening for heavy rainfall. The midlevel moistening, arising from both dynamics and physics, may be most important. Accurately representing many processes may be necessary but not sufficient for capturing the MJO, which suggests that models fail to predict the MJO for a broad range of reasons and limits the possibility of finding a panacea.

Decadal prediction is one focus of the upcoming 5th IPCC Assessment report. To be able to interpret the results and to further improve the decadal predictions it is important to investigate the potential predictability in the participating climate models. This study analyzes the upper limit of climate predictability on decadal time scales and its dependency on sea ice albedo parameterization by performing two perfect ensemble experiments with the global coupled climate model EC-Earth. In the first experiment, the standard albedo formulation of EC-Earth is used, in the second experiment sea ice albedo is reduced. The potential prognostic predictability is analyzed for a set of oceanic and atmospheric parameters. The decadal predictability of the atmospheric circulation is small. The highest potential predictability was found in air temperature at 2 m height over the northern North Atlantic and the southern South Atlantic. Over land, only a few areas are significantly predictable. The predictability for continental size averages of air temperature is relatively good in all northern hemisphere regions. Sea ice thickness is highly predictable along the ice edges in the North Atlantic Arctic Sector. The meridional overturning circulation is highly predictable in both experiments and governs most of the decadal climate predictability in the northern hemisphere. The experiments using reduced sea ice albedo show some important differences like a generally higher predictability of atmospheric variables in the Arctic or higher predictability of air temperature in Europe. Furthermore, decadal variations are substantially smaller in the simulations with reduced ice albedo, which can be explained by reduced sea ice thickness in these simulations.

Seasonal prediction skill of winter mid and high northern latitudes climate from sea ice variations in eight different Arctic regions is analyzed using detrended ERA-interim data and satellite sea ice data for the period 1980-2013. We find significant correlations between ice areas in both September and November and winter sea level pressure, air temperature and precipitation. The prediction skill is improved when using November sea ice conditions as predictor compared to September. This is particularly true for predicting winter NAO-like patterns and blocking situations in the Euro-Atlantic area. We find that sea ice variations in Barents Sea seem to be most important for the sign of the following winter NAO-negative after low ice-but amplitude and extension of the patterns are modulated by Greenland and Labrador Seas ice areas. November ice variability in the Greenland Sea provides the best prediction skill for central and western European temperature and ice variations in the Laptev/East Siberian Seas have the largest impact on the blocking number in the Euro-Atlantic region. Over North America, prediction skill is largest using September ice areas from the Pacific Arctic sector as predictor. Composite analyses of high and low regional autumn ice conditions reveal that the atmospheric response is not entirely linear suggesting changing predictive skill dependent on sign and amplitude of the anomaly. The results confirm the importance of realistic sea ice initial conditions for seasonal forecasts. However, correlations do seldom exceed 0.6 indicating that Arctic sea ice variations can only explain a part of winter climate variations in northern mid and high latitudes.

We investigated high-resolution simulations of regional climate models (RCMs) driven by ERA-40 reanalyses over areas of selected European countries (Austria, Czech Republic, Hungary, Slovakia and Romania) for the period 1961-1990. RCMs were run at a spatial resolution of 10 km in the framework of the CECILIA project, and their outputs were compared with the EOBS dataset of gridded observations and RCM simulations at coarser 25 km resolution from the ENSEMBLES project to identify a possible gain from the CECILIA experiments over ENSEMBLES. Cold biases of air temperature and wet biases of precipitation dominate in the CECILIA simulations. Spatial variability and distribution of the air temperature field are well captured. The precipitation field, relative to observations, often shows inadequately small spatial variability and lowered correlations but is nevertheless comparable to the ENSEMBLES model. Inter-annual variability (IAV) of air temperature is captured differently among seasons but mostly improved in CECILIA compared with ENSEMBLES. Precipitation IAV shows a similar or worse score. The detected weaknesses found within the validation of the CECILIA RCMs are attributed to the resolution dependence of the set of physical parameterizations in the models and the choice of integration domain. The gain obtained by using a high resolution over a small domain (as in CECILIA) relative to a lower resolution (25 km) over a larger domain (as in ENSEMBLES) is clear for air temperature but limited for precipitation.

We present the first climate prediction of the coming decade made with multiple models, initialized with prior observations. This prediction accrues from an international activity to exchange decadal predictions in near real-time, in order to assess differences and similarities, provide a consensus view to prevent over-confidence in forecasts from any single model, and establish current collective capability. We stress that the forecast is experimental, since the skill of the multi-model system is as yet unknown. Nevertheless, the forecast systems used here are based on models that have undergone rigorous evaluation and individually have been evaluated for forecast skill. Moreover, it is important to publish forecasts to enable open evaluation, and to provide a focus on climate change in the coming decade. Initialized forecasts of the year 2011 agree well with observations, with a pattern correlation of 0.62 compared to 0.31 for uninitialized projections. In particular, the forecast correctly predicted La Nia in the Pacific, and warm conditions in the north Atlantic and USA. A similar pattern is predicted for 2012 but with a weaker La Nia. Indices of Atlantic multi-decadal variability and Pacific decadal variability show no signal beyond climatology after 2015, while temperature in the Nio3 region is predicted to warm slightly by about 0.5 A degrees C over the coming decade. However, uncertainties are large for individual years and initialization has little impact beyond the first 4 years in most regions. Relative to uninitialized forecasts, initialized forecasts are significantly warmer in the north Atlantic sub-polar gyre and cooler in the north Pacific throughout the decade. They are also significantly cooler in the global average and over most land and ocean regions out to several years ahead. However, in the absence of volcanic eruptions, global temperature is predicted to continue to rise, with each year from 2013 onwards having a 50 % chance of exceeding the current observed record. Verification of these forecasts will provide an important opportunity to test the performance of models and our understanding and knowledge of the drivers of climate change.

An analysis of diabatic heating and moistening processes from 12 to 36h lead time forecasts from 12 Global Circulation Models are presented as part of the Vertical structure and physical processes of the Madden-Julian Oscillation (MJO) project. A lead time of 12-36h is chosen to constrain the large-scale dynamics and thermodynamics to be close to observations while avoiding being too close to the initial spin-up of the models as they adjust to being driven from the Years of Tropical Convection (YOTC) analysis. A comparison of the vertical velocity and rainfall with the observations and YOTC analysis suggests that the phases of convection associated with the MJO are constrained in most models at this lead time although the rainfall in the suppressed phase is typically overestimated. Although the large-scale dynamics is reasonably constrained, moistening and heating profiles have large intermodel spread. In particular, there are large spreads in convective heating and moistening at midlevels during the transition to active convection. Radiative heating and cloud parameters have the largest relative spread across models at upper levels during the active phase. A detailed analysis of time step behavior shows that some models show strong intermittency in rainfall and differences in the precipitation and dynamics relationship between models. The wealth of model outputs archived during this project is a very valuable resource for model developers beyond the study of the MJO. In addition, the findings of this study can inform the design of process model experiments, and inform the priorities for field experiments and future observing systems.