Increased frequency of extreme La Niña events under greenhouse warming

The El Niño/Southern Oscillation is Earth’s most prominent source of interannual climate variability, alternating irregularly between El Niño and La Niña, and resulting in global disruption of weather patterns, ecosystems, fisheries and agriculture1, 2, 3, 4, 5. The 1998–1999 extreme La Niña event that followed the 1997–1998 extreme El Niño event6 switched extreme El Niño-induced severe droughts to devastating floods in western Pacific countries, and vice versa in the southwestern United States4, 7. During extreme La Niña events, cold sea surface conditions develop in the central Pacific8, 9, creating an enhanced temperature gradient from the Maritime continent to the central Pacific. Recent studies have revealed robust changes in El Niño characteristics in response to simulated future greenhouse warming10, 11, 12, but how La Niña will change remains unclear. Here we present climate modelling evidence, from simulations conducted for the Coupled Model Intercomparison Project phase 5 (ref. 13), for a near doubling in the frequency of future extreme La Niña events, from one in every 23 years to one in every 13 years. This occurs because projected faster mean warming of the Maritime continent than the central Pacific, enhanced upper ocean vertical temperature gradients, and increased frequency of extreme El Niño events are conducive to development of the extreme La Niña events. Approximately 75% of the increase occurs in years following extreme El Niño events, thus projecting more frequent swings between opposite extremes from one year to the next.

Response of El Niño sea surface temperature variability to greenhouse warming

The destructive environmental and socio-economic impacts of the El Niño/Southern Oscillation1, 2 (ENSO) demand an improved understanding of how ENSO will change under future greenhouse warming. Robust projected changes in certain aspects of ENSO have been recently established3, 4, 5. However, there is as yet no consensus on the change in the magnitude of the associated sea surface temperature (SST) variability6, 7, 8, commonly used to represent ENSO amplitude1, 6, despite its strong effects on marine ecosystems and rainfall worldwide1, 2, 3, 4, 9. Here we show that the response of ENSO SST amplitude is time-varying, with an increasing trend in ENSO amplitude before 2040, followed by a decreasing trend thereafter. We attribute the previous lack of consensus to an expectation that the trend in ENSO amplitude over the entire twenty-first century is unidirectional, and to unrealistic model dynamics of tropical Pacific SST variability. We examine these complex processes across 22 models in the Coupled Model Intercomparison Project phase 5 (CMIP5) database10, forced under historical and greenhouse warming conditions. The nine most realistic models identified show a strong consensus on the time-varying response and reveal that the non-unidirectional behaviour is linked to a longitudinal difference in the surface warming rate across the Indo-Pacific basin. Our results carry important implications for climate projections and climate adaptation pathways.

Increasing frequency of extreme El Niño events due to greenhouse warming

El Niño events are a prominent feature of climate variability with global climatic impacts. The 1997/98 episode, often referred to as ‘the climate event of the twentieth century’1, 2, and the 1982/83 extreme El Niño3, featured a pronounced eastward extension of the west Pacific warm pool and development of atmospheric convection, and hence a huge rainfall increase, in the usually cold and dry equatorial eastern Pacific. Such a massive reorganization of atmospheric convection, which we define as an extreme El Niño, severely disrupted global weather patterns, affecting ecosystems4, 5, agriculture6, tropical cyclones, drought, bushfires, floods and other extreme weather events worldwide3, 7, 8, 9. Potential future changes in such extreme El Niño occurrences could have profound socio-economic consequences. Here we present climate modelling evidence for a doubling in the occurrences in the future in response to greenhouse warming. We estimate the change by aggregating results from climate models in the Coupled Model Intercomparison Project phases 3 (CMIP3; ref. 10) and 5 (CMIP5; ref. 11) multi-model databases, and a perturbed physics ensemble12. The increased frequency arises from a projected surface warming over the eastern equatorial Pacific that occurs faster than in the surrounding ocean waters13, 14, facilitating more occurrences of atmospheric convection in the eastern equatorial region.

Urban warming in villages

Long term meteorological records (> 100 years) from stations associated with villages are generally classified as rural and assumed to have no urban influence. Using networks installed in two European villages, the local and microclimatic variations around two of these rural-village sites are examined. An annual average temperature difference ($\Delta{T}$) of 0.6 and 0.4 K was observed between the built-up village area and the current meteorological station in Geisenheim (Germany) and Haparanda (Sweden), respectively. Considerably larger values were recorded for the minimum temperatures and during summer. The spatial variations in temperature within the villages are of the same order as recorded over the past 100+ years in these villages (0.06 to 0.17 K/10 years). This suggests that the potential biases in the long records of rural-villages also warrant careful consideration like those of the more commonly studied large urban areas effects.

Examining the predictability of the Stratospheric Sudden Warming of January 2013 using multiple NWP systems

The first multi-model study to estimate the predictability of a boreal Sudden Stratospheric Warming (SSW) is performed using five NWP systems. During the 2012-2013 boreal winter, anomalous upward propagating planetary wave activity was observed towards the end of December, which followed by a rapid deceleration of the westerly circulation around 2 January 2013, and on 7 January 2013 the zonal mean zonal wind at 60°N and 10 hPa reversed to easterly. This stratospheric dynamical activity was followed by an equatorward shift of the tropospheric jet stream and by a high pressure anomaly over the North Atlantic, which resulted in severe cold conditions in the UK and Northern Europe. In most of the five models, the SSW event was predicted 10 days in advance. However, only some ensemble members in most of the models predicted weakening of westerly wind when the models were initialized 15 days in advance of the SSW. Further dynamical analysis of the SSW shows that this event was characterized by the anomalous planetary wave-1 amplification followed by the anomalous wave-2 amplification in the stratosphere, which resulted in a split vortex occurring between 6 January 2013 and 8 January 2013. The models have some success in reproducing wave-1 activity when initialized 15 days in advance, they but generally failed to produce the wave-2 activity during the final days of the event. Detailed analysis shows that models have reasonably good skill in forecasting tropospheric blocking features that stimulate wave-2 amplification in the troposphere, but they have limited skill in reproducing wave-2 amplification in the stratosphere.

New use of global warming potentials to compare cumulative and short-lived climate pollutants

Parties to the United Nations Framework Convention on Climate Change (UNFCCC) have requested guidance on common greenhouse gas metrics in accounting for Nationally determined contributions (NDCs) to emission reductions1. Metric choice can affect the relative emphasis placed on reductions of ‘cumulative climate pollutants’ such as carbon dioxide versus ‘short-lived climate pollutants’ (SLCPs), including methane and black carbon2, 3, 4, 5, 6. Here we show that the widely used 100-year global warming potential (GWP100) effectively measures the relative impact of both cumulative pollutants and SLCPs on realized warming 20–40 years after the time of emission. If the overall goal of climate policy is to limit peak warming, GWP100 therefore overstates the importance of current SLCP emissions unless stringent and immediate reductions of all climate pollutants result in temperatures nearing their peak soon after mid-century7, 8, 9, 10, which may be necessary to limit warming to “well below 2 °C” (ref. 1). The GWP100 can be used to approximately equate a one-off pulse emission of a cumulative pollutant and an indefinitely sustained change in the rate of emission of an SLCP11, 12, 13. The climate implications of traditional CO2-equivalent targets are ambiguous unless contributions from cumulative pollutants and SLCPs are specified separately.

Understanding the rapid summer warming and changes in temperature extremes since the mid-1990s over Western Europe

Analysis of observations indicates that there was a rapid increase in summer (June-August, JJA) mean surface air temperature (SAT) since the mid-1990s over Western Europe. Accompanying this rapid warming are significant increases in summer mean daily maximum temperature, daily minimum temperature, annual hottest day temperature and warmest night temperature, and an increase in frequency of summer days and tropical nights, while the change in the diurnal temperature range (DTR) is small. This study focuses on understanding causes of the rapid summer warming and associated temperature extreme changes. A set of experiments using the atmospheric component of the state-of-the-art HadGEM3 global climate model have been carried out to quantify relative roles of changes in sea surface temperature (SST)/sea ice extent (SIE), anthropogenic greenhouse gases (GHGs), and anthropogenic aerosols (AAer). Results indicate that the model forced by changes in all forcings reproduces many of the observed changes since the mid-1990s over Western Europe. Changes in SST/SIE explain 62.2% ± 13.0% of the area averaged seasonal mean warming signal over Western Europe, with the remaining 37.8% ± 13.6% of the warming explained by the direct impact of changes in GHGs and AAer. Results further indicate that the direct impact of the reduction of AAer precursor emissions over Europe, mainly through aerosol-radiation interaction with additional contributions from aerosol-cloud interaction and coupled atmosphere-land surface feedbacks, is a key factor for increases in annual hottest day temperature and in frequency of summer days. It explains 45.5% ± 17.6% and 40.9% ± 18.4% of area averaged signals for these temperature extremes. The direct impact of the reduction of AAer precursor emissions over Europe acts to increase DTR locally, but the change in DTR is countered by the direct impact of GHGs forcing. In the next few decades, greenhouse gas concentrations will continue to rise and AAer precursor emissions over Europe and North America will continue to decline. Our results suggest that the changes in summer seasonal mean SAT and temperature extremes over Western Europe since the mid-1990s are most likely to be sustained or amplified in the near term, unless other factors intervene.

Elevation-dependent warming in mountain regions of the world

There is growing evidence that the rate of warming is amplified with elevation, such that high-mountain environments experience more rapid changes in temperature than environments at lower elevations. Elevation-dependent warming (EDW) can accelerate the rate of change in mountain ecosystems, cryospheric systems, hydrological regimes and biodiversity. Here we review important mechanisms that contribute towards EDW: snow albedo and surface-based feedbacks; water vapour changes and latent heat release; surface water vapour and radiative flux changes; surface heat loss and temperature change; and aerosols. All lead to enhanced warming with elevation (or at a critical elevation), and it is believed that combinations of these mechanisms may account for contrasting regional patterns of EDW. We discuss future needs to increase knowledge of mountain temperature trends and their controlling mechanisms through improved observations, satellite-based remote sensing and model simulations.

Making sense of the early-2000s warming slowdown

It has been claimed that the early-2000s global warming slowdown or hiatus, characterized by a reduced rate of global surface warming, has been overstated, lacks sound scientific basis, or is unsupported by observations. The evidence presented here contradicts these claims.

Can climate trends be calculated from reanalysis data?

Several global quantities are computed from the ERA40 reanalysis for the period 1958-2001 and explored for trends. These are discussed in the context of changes to the global observing system. Temperature, integrated water vapor (IWV), and kinetic energy are considered. The ERA40 global mean temperature in the lower troposphere has a trend of +0.11 K per decade over the period of 1979-2001, which is slightly higher than the MSU measurements, but within the estimated error limit. For the period 1958 2001 the warming trend is 0.14 K per decade but this is likely to be an artifact of changes in the observing system. When this is corrected for, the warming trend is reduced to 0.10 K per decade. The global trend in IWV for the period 1979-2001 is +0.36 mm per decade. This is about twice as high as the trend determined from the Clausius-Clapeyron relation assuming conservation of relative humidity. It is also larger than results from free climate model integrations driven by the same observed sea surface temperature as used in ERA40. It is suggested that the large trend in IWV does not represent a genuine climate trend but an artifact caused by changes in the global observing system such as the use of SSM/I and more satellite soundings in later years. Recent results are in good agreement with GPS measurements. The IWV trend for the period 1958-2001 is still higher but reduced to +0.16 mm per decade when corrected for changes in the observing systems. Total kinetic energy shows an increasing global trend. Results from data assimilation experiments strongly suggest that this trend is also incorrect and mainly caused by the huge changes in the global observing system in 1979. When this is corrected for, no significant change in global kinetic energy from 1958 onward can be found.

Storm tracks and climate change

Extratropical and tropical transient storm tracks are investigated from the perspective of feature tracking in the ECHAM5 coupled climate model for the current and a future climate scenario. The atmosphere-only part of the model, forced by observed boundary conditions, produces results that agree well with analyses from the 40-yr ECMWF Re-Analysis (ERA-40), including the distribution of storms as a function of maximum intensity. This provides the authors with confidence in the use of the model for the climate change experiments. The statistical distribution of storm intensities is virtually preserved under climate change using the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenarios (SRES) A1B scenario until the end of this century. There are no indications in this study of more intense storms in the future climate, either in the Tropics or extratropics, but rather a minor reduction in the number of weaker storms. However, significant changes occur on a regional basis in the location and intensity of storm tracks. There is a clear poleward shift in the Southern Hemisphere with consequences of reduced precipitation for several areas, including southern Australia. Changes in the Northern Hemisphere are less distinct, but there are also indications of a poleward shift, a weakening of the Mediterranean storm track, and a strengthening of the storm track north of the British Isles. The tropical storm tracks undergo considerable changes including a weakening in the Atlantic sector and a strengthening and equatorward shift in the eastern Pacific. It is suggested that some of the changes, in particular the tropical ones, are due to an SST warming maximum in the eastern Pacific. The shift in the extratropical storm tracks is shown to be associated with changes in the zonal SST gradient in particular for the Southern Hemisphere.

A new feedback on climate change from the hydrological cycle

An intensification of the hydrological cycle is a likely consequence of global warming. But changes in the hydrological cycle could affect sea-surface temperature by modifying diffusive ocean heat transports. We investigate this mechanism by studying a coupled general circulation model sensitivity experiment in which the hydrological cycle is artificially amplified. We find that the amplified hydrological cycle depresses sea-surface temperature by enhancing ocean heat uptake in low latitudes. We estimate that a 10% increase in the hydrological cycle will contribute a basin-scale sea-surface temperature decrease of around 0.1°C away from high latitudes, with larger decreases locally. We conclude that an intensified hydrological cycle is likely to contribute a weak negative feedback to anthropogenic climate change.

Enhancement of ENSO variability by a weakened Atlantic thermohaline circulation in a coupled GCM

A coupled ocean–atmosphere general circulation model is used to investigate the modulation of El Niño–Southern Oscillation (ENSO) variability due to a weakened Atlantic thermohaline circulation (THC). The THC weakening is induced by freshwater perturbations in the North Atlantic, and leads to a well-known sea surface temperature dipole and a southward shift of the intertropical convergence zone (ITCZ) in the tropical Atlantic. Through atmospheric teleconnections and local coupled air–sea feedbacks, a meridionally asymmetric mean state change is generated in the eastern equatorial Pacific, corresponding to a weakened annual cycle, and westerly anomalies develop over the central Pacific. The westerly anomalies are associated with anomalous warming of SST, causing an eastward extension of the west Pacific warm pool particularly in August–February, and enhanced precipitation. These and other changes in the mean state lead in turn to an eastward shift of the zonal wind anomalies associated with El Niño events, and a significant increase in ENSO variability. In response to a 1-Sv (1 Sv ≡ 106 m3 s−1) freshwater input in the North Atlantic, the THC slows down rapidly and it weakens by 86% over years 50–100. The Niño-3 index standard deviation increases by 36% during the first 100-yr simulation relative to the control simulation. Further analysis indicates that the weakened THC not only leads to a stronger ENSO variability, but also leads to a stronger asymmetry between El Niño and La Niña events. This study suggests a role for an atmospheric bridge that rapidly conveys the influence of the Atlantic Ocean to the tropical Pacific and indicates that fluctuations of the THC can mediate not only mean climate globally but also modulate interannual variability. The results may contribute to understanding both the multidecadal variability of ENSO activity during the twentieth century and longer time-scale variability of ENSO, as suggested by some paleoclimate records.

Monitoring present day changes in water vapour and the radiative energy balance using satellite data, reanalyses and models

A combination of satellite data, reanalysis products and climate models are combined to monitor changes in water vapour, clear-sky radiative cooling of the atmosphere and precipitation over the period 1979-2006. Climate models are able to simulate observed increases in column integrated water vapour (CWV) with surface temperature (Ts) over the ocean. Changes in the observing system lead to spurious variability in water vapour and clear-sky longwave radiation in reanalysis products. Nevertheless all products considered exhibit a robust increase in clear-sky longwave radiative cooling from the atmosphere to the surface; clear-sky longwave radiative cooling of the atmosphere is found to increase with Ts at the rate of ~4 Wm-2 K-1 over tropical ocean regions of mean descending vertical motion. Precipitation (P) is tightly coupled to atmospheric radiative cooling rates and this implies an increase in P with warming at a slower rate than the observed increases in CWV. Since convective precipitation depends on moisture convergence, the above implies enhanced precipitation over convective regions and reduced precipitation over convectively suppressed regimes. To quantify this response, observed and simulated changes in precipitation rate are analysed separately over regions of mean ascending and descending vertical motion over the tropics. The observed response is found to be substantially larger than the model simulations and climate change projections. It is currently not clear whether this is due to deficiencies in model parametrizations or errors in satellite retrievals.

Historical reconstruction of the Atlantic Meridional Overturning Circulation from the ECMWF operational ocean reanalysis

A reconstruction of the Atlantic Meridional Overturning Circulation (MOC) for the period 1959–2006 has been derived from the ECMWF operational ocean reanalysis. The reconstruction shows a wide range of time-variability, including a downward trend. At 26N, both the MOC intensity and changes in its vertical structure are in good agreement with previous estimates based on trans-Atlantic surveys. At 50N, the MOC and strength of the subpolar gyre are correlated at interannual time scales, but show opposite secular trends. Heat transport variability is highly correlated with the MOC but shows a smaller trend due to the warming of the upper ocean, which partially compensates for the weakening of the circulation. Results from sensitivity experiments show that although the time-varying upper boundary forcing provides useful MOC information, the sequential assimilation of ocean data further improves the MOC estimation by increasing both the mean and the time variability.