ENSO representation in climate models: from CMIP3 to CMIP5

We analyse the ability of CMIP3 and CMIP5 coupled ocean–atmosphere general circulation models (CGCMs) to simulate the tropical Pacific mean state and El Niño-Southern Oscillation (ENSO). The CMIP5 multi-model ensemble displays an encouraging 30 % reduction of the pervasive cold bias in the western Pacific, but no quantum leap in ENSO performance compared to CMIP3. CMIP3 and CMIP5 can thus be considered as one large ensemble (CMIP3 + CMIP5) for multi-model ENSO analysis. The too large diversity in CMIP3 ENSO amplitude is however reduced by a factor of two in CMIP5 and the ENSO life cycle (location of surface temperature anomalies, seasonal phase locking) is modestly improved. Other fundamental ENSO characteristics such as central Pacific precipitation anomalies however remain poorly represented. The sea surface temperature (SST)-latent heat flux feedback is slightly improved in the CMIP5 ensemble but the wind-SST feedback is still underestimated by 20–50 % and the shortwave-SST feedbacks remain underestimated by a factor of two. The improvement in ENSO amplitudes might therefore result from error compensations. The ability of CMIP models to simulate the SST-shortwave feedback, a major source of erroneous ENSO in CGCMs, is further detailed. In observations, this feedback is strongly nonlinear because the real atmosphere switches from subsident (positive feedback) to convective (negative feedback) regimes under the effect of seasonal and interannual variations. Only one-third of CMIP3 + CMIP5 models reproduce this regime shift, with the other models remaining locked in one of the two regimes. The modelled shortwave feedback nonlinearity increases with ENSO amplitude and the amplitude of this feedback in the spring strongly relates with the models ability to simulate ENSO phase locking. In a final stage, a subset of metrics is proposed in order to synthesize the ability of each CMIP3 and CMIP5 models to simulate ENSO main characteristics and key atmospheric feedbacks.

Mid-holocene and last glacial maximum climate simulations with the IPSL model-part I: comparing IPSLCM5A to IPSLCM4

The climates of the mid-Holocene (MH), 6,000 years ago, and of the Last Glacial Maximum (LGM), 21,000 years ago, have extensively been simulated, in particular in the framework of the Palaeoclimate Modelling Intercomparion Project. These periods are well documented by paleo-records, which can be used for evaluating model results for climates different from the present one. Here, we present new simulations of the MH and the LGM climates obtained with the IPSL_CM5A model and compare them to our previous results obtained with the IPSL_CM4 model. Compared to IPSL_CM4, IPSL_CM5A includes two new features: the interactive representation of the plant phenology and marine biogeochemistry. But one of the most important differences between these models is the latitudinal resolution and vertical domain of their atmospheric component, which have been improved in IPSL_CM5A and results in a better representation of the mid-latitude jet-streams. The Asian monsoon’s representation is also substantially improved. The global average mean annual temperature simulated for the pre-industrial (PI) period is colder in IPSL_CM5A than in IPSL_CM4 but their climate sensitivity to a CO2 doubling is similar. Here we show that these differences in the simulated PI climate have an impact on the simulated MH and LGM climatic anomalies. The larger cooling response to LGM boundary conditions in IPSL_CM5A appears to be mainly due to differences between the PMIP3 and PMIP2 boundary conditions, as shown by a short wave radiative forcing/feedback analysis based on a simplified perturbation method. It is found that the sensitivity computed from the LGM climate is lower than that computed from 2 × CO2 simulations, confirming previous studies based on different models. For the MH, the Asian monsoon, stronger in the IPSL_CM5A PI simulation, is also more sensitive to the insolation changes. The African monsoon is also further amplified in IPSL_CM5A due to the impact of the interactive phenology. Finally the changes in variability for both models and for MH and LGM are presented taking the example of the El-Niño Southern Oscillation (ENSO), which is very different in the PI simulations. ENSO variability is damped in both model versions at the MH, whereas inconsistent responses are found between the two versions for the LGM. Part 2 of this paper examines whether these differences between IPSL_CM4 and IPSL_CM5A can be distinguished when comparing those results to palaeo-climatic reconstructions and investigates new approaches for model-data comparisons made possible by the inclusion of new components in IPSL_CM5A.

The sensitivity of convective aggregation to diabatic processes in idealized radiative-convective equilibrium simulations

Idealized explicit convection simulations of the Met Office Unified Model exhibit spontaneous self-aggregation in radiative-convective equilibrium, as seen in other models in previous studies. This self-aggregation is linked to feedbacks between radiation, surface fluxes, and convection, and the organization is intimately related to the evolution of the column water vapor field. Analysis of the budget of the spatial variance of column-integrated frozen moist static energy (MSE), following Wing and Emanuel [2014], reveals that the direct radiative feedback (including significant cloud longwave effects) is dominant in both the initial development of self-aggregation and the maintenance of an aggregated state. A low-level circulation at intermediate stages of aggregation does appear to transport MSE from drier to moister regions, but this circulation is mostly balanced by other advective effects of opposite sign and is forced by horizontal anomalies of convective heating (not radiation). Sensitivity studies with either fixed prescribed radiative cooling, fixed prescribed surface fluxes, or both do not show full self-aggregation from homogeneous initial conditions, though fixed surface fluxes do not disaggregate an initialized aggregated state. A sensitivity study in which rain evaporation is turned off shows more rapid self-aggregation, while a run with this change plus fixed radiative cooling still shows strong self-aggregation, supporting a “moisture memory” effect found in Muller and Bony [2015]. Interestingly, self-aggregation occurs even in simulations with sea surface temperatures (SSTs) of 295 K and 290 K, with direct radiative feedbacks dominating the budget of MSE variance, in contrast to results in some previous studies.

Enhanced long-range forecast skill in boreal winter following stratospheric strong vortex conditions

There has been a great deal of recent interest in producing weather forecasts on the 2–6 week sub-seasonal timescale, which bridges the gap between medium-range (0–10 day) and seasonal (3–6 month) forecasts. While much of this interest is focused on the potential applications of skilful forecasts on the sub-seasonal range, understanding the potential sources of sub-seasonal forecast skill is a challenging and interesting problem, particularly because of the likely state-dependence of this skill (Hudson et al 2011). One such potential source of state-dependent skill for the Northern Hemisphere in winter is the occurrence of stratospheric sudden warming (SSW) events (Sigmond et al 2013). Here we show, by analysing a set of sub-seasonal hindcasts, that there is enhanced predictability of surface circulation not only when the stratospheric vortex is anomalously weak following SSWs but also when the vortex is extremely strong. Sub-seasonal forecasts initialized during strong vortex events are able to successfully capture the associated surface temperature and circulation anomalies. This results in an enhancement of Northern annular mode forecast skill compared to forecasts initialized during the cases when the stratospheric state is close to climatology. We demonstrate that the enhancement of skill for forecasts initialized during periods of strong vortex conditions is comparable to that achieved for forecasts initialized during weak events. This result indicates that additional confidence can be placed in sub-seasonal forecasts when the stratospheric polar vortex is significantly disturbed from its normal state.

The combined influences of westerly phase of the quasi biennial oscillation and 11-year solar maximum conditions on the Northern Hemisphere extratropical winter circulation

The combined influences of the westerly phase of the quasi-biennial oscillation (QBO-W) and solar maximum (Smax) conditions on the Northern Hemisphere extratropical winter circulation are investigated using reanalysis data and Center for Climate System Research/National Institute for Environmental Studies chemistry climate model (CCM) simulations. The composite analysis for the reanalysis data indicates strengthened polar vortex in December followed by weakened polar vortex in February–March for QBO-W during Smax (QBO-W/Smax) conditions. This relationship need not be specific to QBO-W/Smax conditions but may just require strengthened vortex in December, which is more likely under QBO-W/Smax. Both the reanalysis data and CCM simulations suggest that dynamical processes of planetary wave propagation and meridional circulation related to QBO-W around polar vortex in December are similar in character to those related to Smax; furthermore, both processes may work in concert to maintain stronger vortex during QBO-W/Smax. In the reanalysis data, the strengthened polar vortex in December is associated with the development of north–south dipole tropospheric anomaly in the Atlantic sector similar to the North Atlantic oscillation (NAO) during December–January. The structure of the north–south dipole anomaly has zonal wavenumber 1 (WN1) component, where the longitude of anomalous ridge overlaps with that of climatological ridge in the North Atlantic in January. This implies amplification of the WN1 wave and results in the enhancement of the upward WN1 propagation from troposphere into stratosphere in January, leading to the weakened polar vortex in February–March. Although WN2 waves do not play a direct role in forcing the stratospheric vortex evolution, their tropospheric response to QBO-W/Smax conditions appears to be related to the maintenance of the NAO-like anomaly in the high-latitude troposphere in January. These results may provide a possible explanation for the mechanisms underlying the seasonal evolution of wintertime polar vortex anomalies during QBO-W/Smax conditions and the role of troposphere in this evolution.

An assessment of upper ocean salinity content from the ocean reanalyses inter-comparison project (ORA-IP)

Many institutions worldwide have developed ocean reanalyses systems (ORAs) utilizing a variety of ocean models and assimilation techniques. However, the quality of salinity reanalyses arising from the various ORAs has not yet been comprehensively assessed. In this study, we assess the upper ocean salinity content (depth-averaged over 0–700 m) from 14 ORAs and 3 objective ocean analysis systems (OOAs) as part of the Ocean Reanalyses Intercomparison Project. Our results show that the best agreement between estimates of salinity from different ORAs is obtained in the tropical Pacific, likely due to relatively abundant atmospheric and oceanic observations in this region. The largest disagreement in salinity reanalyses is in the Southern Ocean along the Antarctic circumpolar current as a consequence of the sparseness of both atmospheric and oceanic observations in this region. The West Pacific warm pool is the largest region where the signal to noise ratio of reanalysed salinity anomalies is >1. Therefore, the current salinity reanalyses in the tropical Pacific Ocean may be more reliable than those in the Southern Ocean and regions along the western boundary currents. Moreover, we found that the assimilation of salinity in ocean regions with relatively strong ocean fronts is still a common problem as seen in most ORAs. The impact of the Argo data on the salinity reanalyses is visible, especially within the upper 500m, where the interannual variability is large. The increasing trend in global-averaged salinity anomalies can only be found within the top 0–300m layer, but with quite large diversity among different ORAs. Beneath the 300m depth, the global-averaged salinity anomalies from most ORAs switch their trends from a slightly growing trend before 2002 to a decreasing trend after 2002. The rapid switch in the trend is most likely an artefact of the dramatic change in the observing system due to the implementation of Argo.

Interannual-decadal variability of wintertime mixed layer depths in the North Pacific detected by an ensemble of ocean syntheses

The interannual-decadal variability of the wintertime mixed layer depths (MLDs) over the North Pacific is investigated from an empirical orthogonal function (EOF) analysis of an ensemble of global ocean reanalyses. The first leading EOF mode represents the interannual MLD anomalies centered in the eastern part of the central mode water formation region in phase opposition with those in the eastern subtropics and the central Alaskan Gyre. This first EOF mode is highly correlated with the Pacific decadal oscillation index on both the interannual and decadal time scales. The second leading EOF mode represents the MLD variability in the subtropical mode water (STMW) formation region and has a good correlation with the wintertime West Pacific (WP) index with time lag of 3 years, suggesting the importance of the oceanic dynamical response to the change in the surface wind field associated with the meridional shifts of the Aleutian Low. The above MLD variabilities are in basic agreement with previous observational and modeling findings. Moreover the reanalysis ensemble provides uncertainty estimates. The interannual MLD anomalies in the first and second EOF modes are consistently represented by the individual reanalyses and the amplitudes of the variabilities generally exceed the ensemble spread of the reanalyses. Besides, the resulting MLD variability indices, spanning the 1948–2012 period, should be helpful for characterizing the North Pacific climate variability. In particular, a 6-year oscillation including the WP teleconnection pattern in the atmosphere and the oceanic MLD variability in the STMW formation region is first detected.

Interdecadal changes on the seasonal prediction of the western North Pacific summer climate around the late 1970s and early 1990s

Identifying predictability and the corresponding sources for the western North Pacific (WNP) summer climate in the case of non-stationary teleconnections during recent decades benefits for further improvements of long-range prediction on the WNP and East Asian summers. In the past few decades, pronounced increases on the summer sea surface temperature (SST) and associated interannual variability are observed over the tropical Indian Ocean and eastern Pacific around the late 1970s and over the Maritime Continent and western–central Pacific around the early 1990s. These increases are associated with significant enhancements of the interannual variability for the lower-tropospheric wind over the WNP. In this study, we further assess interdecadal changes on the seasonal prediction of the WNP summer anomalies, using May-start retrospective forecasts from the ENSEMBLES multi-model project in the period 1960–2005. It is found that prediction of the WNP summer anomalies exhibits an interdecadal shift with higher prediction skills since the late 1970s, particularly after the early 1990s. Improvements of the prediction skills for SSTs after the late 1970s are mainly found around tropical Indian Ocean and the WNP. The better prediction of the WNP after the late 1970s may arise mainly from the improvement of the SST prediction around the tropical eastern Indian Ocean. The close teleconnections between the tropical eastern Indian Ocean and WNP summer variability work both in the model predictions and observations. After the early 1990s, on the other hand, the improvements are detected mainly around the South China Sea and Philippines for the lower-tropospheric zonal wind and precipitation anomalies, associating with a better description of the SST anomalies around the Maritime Continent. A dipole SST pattern over the Maritime Continent and the central equatorial Pacific Ocean is closely related to the WNP summer anomalies after the early 1990s. This teleconnection mode is quite predictable, which is realistically reproduced by the models, presenting more predictable signals to the WNP summer climate after the early 1990s.

Climatology and variability of trace gases in extratropical double-tropopause regions from MLS, HIRDLS, and ACE-FTS measurements

Upper tropospheric and lower stratospheric measurements from the Aura Microwave Limb Sounder (MLS), the Aura High Resolution Dynamics Limb Sounder (HIRDLS), and the Atmospheric Chemistry Experiment-Fourier transform spectrometer (ACE-FTS) are used to present the first global climatological comparison of extratropical, nonpolar trace gas distributions in double-tropopause (DT) and single-tropopause (ST) regions. Stratospheric tracers, O3, HNO3, and HCl, have lower mixing ratios ∼2–8 km above the primary (lowermost) tropopause in DT than in ST regions in all seasons, with maximum Northern Hemisphere (NH) differences near 50% in winter and 30% in summer. Southern Hemisphere winter differences are somewhat smaller, but summer differences are similar in the two hemispheres. H2O in DT regions of both hemispheres shows strong negative anomalies in November through February and positive anomalies in July through October, reflecting the strong seasonal cycle in H2O near the tropical tropopause. CO and other tropospheric tracers examined have higher DT than ST values 2–7 km above the primary tropopause, with the largest differences in winter. Large DT-ST differences extend to high NH latitudes in fall and winter, with longitudinal maxima in regions associated with enhanced wave activity and subtropical jet variations. Results for O3 and HNO3 agree closely between MLS and HIRDLS, and differences from ACE-FTS are consistent with its sparse and irregular midlatitude sampling. Consistent signatures in climatological trace gas fields provide strong evidence that transport from the tropical upper troposphere into the layer between double tropopauses is an important pathway for stratosphere-troposphere exchange.

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Variability and trends in Laptev Sea ice outflow between 1992-2011

Variability and trends in seasonal and interannual ice area export out of the Laptev Sea between 1992 and 2011 are investigated using satellite-based sea ice drift and concentration data. We found an average total winter (Octo- ber to May) ice area transport across the northern and east- ern Laptev Sea boundaries (NB and EB) of 3.48 × 10 5 km 2 . The average transport across the NB (2.87 × 10 5 km 2 ) is thereby higher than across the EB (0.61 × 10 5 km 2 ), with a less pronounced seasonal cycle. The total Laptev Sea ice area flux significantly increased over the last decades (0.85 × 10 5 km 2 decade − 1 , p> 0 . 95), dominated by increas- ing export through the EB (0.55 × 10 5 km 2 decade − 1 , p> 0 . 90), while the increase in export across the NB is smaller (0.3 × 10 5 km 2 decade − 1 ) and statistically not significant. The strong coupling between across-boundary SLP gradient and ice drift velocity indicates that monthly variations in ice area flux are primarily controlled by changes in geostrophic wind velocities, although the Laptev Sea ice circulation shows no clear relationship with large-scale atmospheric in- dices. Also there is no evidence of increasing wind velocities that could explain the overall positive trends in ice export. The increased transport rates are rather the consequence of a changing ice cover such as thinning and/or a decrease in con- centration. The use of a back-propagation method revealed that most of the ice that is incorporated into the Transpolar Drift is formed during freeze-up and originates from the cen- tral and western part of the Laptev Sea, while the exchange with the East Siberian Sea is dominated by ice coming from the central and southeastern Laptev Sea. Furthermore, our re- sults imply that years of high ice export in late winter (Febru- ary to May) have a thinning effect on the ice cover, which in turn preconditions the occurence of negative sea ice extent anomalies in summer.

Impact of the Arctic Ocean Atlantic water layer on Siberian shelf hydrography

This paper examines the role of the Arctic Ocean Atlantic water (AW) in modifying the Laptev Sea shelf bottom hydrography on the basis of historical records from 1932 to 2008, field observations carried out in April–May 2008, and 2002–2009 cross‐slope measurements. A climatology of bottom hydrography demonstrates warming that extends offshore from the 30–50 m depth contour. Bottom layer temperature‐time series constructed from historical records links the Laptev Sea outer shelf to the AW boundary current transporting warm and saline water from the North Atlantic. The AW warming of the mid‐1990s and the mid‐2000s is consistent with outer shelf bottom temperature variability. For April–May 2008 we observed on‐shelf near‐bottom warm and saline water intrusions up to the 20 m isobath. These intrusions are typically about 0.2°C warmer and 1–1.5 practical salinity units saltier than ambient water. The 2002–2009 cross‐slope observations are suggestive for the continental slope upward heat flux from the AW to the overlying low‐halocline water (LHW). The lateral on‐shelf wind‐driven transport of the LHW then results in the bottom layer thermohaline anomalies recorded over the Laptev Sea shelf. We also found that polynya‐induced vertical mixing may act as a drainage of the bottom layer, permitting a relatively small portion of the AW heat to be directly released to the atmosphere. Finally, we see no significant warming (up until now) over the Laptev Sea shelf deeper than 10–15 m in the historical record. Future climate change, however, may bring more intrusions of Atlantic‐modified waters with potentially warmer temperature onto the shelf, which could have a critical impact on the stability of offshore submarine permafrost.

Decadal predictability of regional scale wind speed and wind energy potentials over Central Europe

Decadal predictions on timescales from one year to one decade are gaining importance since this time frame falls within the planning horizon of politics, economy and society. The present study examines the decadal predictability of regional wind speed and wind energy potentials in three generations of the MiKlip (‘Mittelfristige Klimaprognosen’) decadal prediction system. The system is based on the global Max-Planck-Institute Earth System Model (MPI-ESM), and the three generations differ primarily in the ocean initialisation. Ensembles of uninitialised historical and yearly initialised hindcast experiments are used to assess the forecast skill for 10 m wind speeds and wind energy output (Eout) over Central Europe with lead times from one year to one decade. With this aim, a statistical-dynamical downscaling (SDD) approach is used for the regionalisation. Its added value is evaluated by comparison of skill scores for MPI-ESM large-scale wind speeds and SDD-simulated regional wind speeds. All three MPI-ESM ensemble generations show some forecast skill for annual mean wind speed and Eout over Central Europe on yearly and multi-yearly time scales. This forecast skill is mostly limited to the first years after initialisation. Differences between the three ensemble generations are generally small. The regionalisation preserves and sometimes increases the forecast skills of the global runs but results depend on lead time and ensemble generation. Moreover, regionalisation often improves the ensemble spread. Seasonal Eout skills are generally lower than for annual means. Skill scores are lowest during summer and persist longest in autumn. A large-scale westerly weather type with strong pressure gradients over Central Europe is identified as potential source of the skill for wind energy potentials, showing a similar forecast skill and a high correlation with Eout anomalies. These results are promising towards the establishment of a decadal prediction system for wind energy applications over Central Europe.

The teleconnection of the tropical Atlantic to Indo-Pacific sea surface temperatures on inter-annual to centennial time scales: a review of recent findings

In this paper, the teleconnections from the tropical Atlantic to the Indo-Pacific region from inter-annual to centennial time scales will be reviewed. Identified teleconnections and hypotheses on mechanisms at work are reviewed and further explored in a century-long pacemaker coupled ocean-atmosphere simulation ensemble. There is a substantial impact of the tropical Atlantic on the Pacific region at inter-annual time scales. An Atlantic Niño (Niña) event leads to rising (sinking) motion in the Atlantic region, which is compensated by sinking (rising) motion in the central-western Pacific. The sinking (rising) motion in the central-western Pacific induces easterly (westerly) surface wind anomalies just to the west, which alter the thermocline. These perturbations propagate eastward as upwelling (downwelling) Kelvin-waves, where they increase the probability for a La Niña (El Niño) event. Moreover, tropical North Atlantic sea surface temperature anomalies are also able to lead La Niña/El Niño development. At multidecadal time scales, a positive (negative) Atlantic Multidecadal Oscillation leads to a cooling (warming) of the eastern Pacific and a warming (cooling) of the western Pacific and Indian Ocean regions. The physical mechanism for this impact is similar to that at inter-annual time scales. At centennial time scales, the Atlantic warming induces a substantial reduction of the eastern Pacific warming even under CO2 increase and to a strong subsurface cooling.

Simulating decadal variability in the North Atlantic Ocean

Observations and climate models suggest significant decadal variability within the North Atlantic subpolar gyre (NA SPG), though observations are sparse and models disagree on the details of this variability. Therefore, it is important to understand 1) the mechanisms of simulated decadal variability, 2) which parts of simulated variability are more faithful representations of reality, and 3) the implications for climate predictions. Here, we investigate the decadal variability in the NA SPG in the state-of-the-art, high resolution (0.25◦ ocean resolution), climate model ‘HadGEM3’. We find a decadal mode with a period of 17 years that explains 30% of the annual variance in related indices. The mode arises due to the advection of heat content anomalies, and shows asymmetries in the timescale of phase reversal between positive and negative phases. A negative feedback from temperature-driven density anomalies in the Labrador Sea (LS) allows for the phase reversal. The North Atlantic Oscillation (NAO), which exhibits the same periodicity, amplifies the mode. The atmosphere-ocean coupling is stronger during positive rather than negative NAO states, explaining the asymmetry. Within the NA SPG, there is potential predictability arising partly from this mode for up to 5 years. There are important similarities between observed and simulated variability, such as the apparent role for the propagation of heat content anomalies. However, observations suggest interannual LS density anomalies are salinity-driven. Salinity control of density would change the temperature feedback to the south, possibly limiting real-world predictive skill in the southern NA SPG with this model. Finally, to understand the diversity of behaviours, we analyse 42 present-generation climate models. Temperature and salinity biases are found to systematically influence the driver of density variability in the LS. Resolution is a good predictor of the biases. The dependence of variability on the background state has important implications for decadal predictions.