A mechanism of internal decadal atlantic ocean variability in a high-resolution coupled climate model

The North Atlantic Ocean subpolar gyre (NA SPG) is an important region for initialising decadal climate forecasts. Climate model simulations and palaeo climate reconstructions have indicated that this region could also exhibit large, internally generated variability on decadal timescales. Understanding these modes of variability, their consistency across models, and the conditions in which they exist, is clearly important for improving the skill of decadal predictions — particularly when these predictions are made with the same underlying climate models. Here we describe and analyse a mode of internal variability in the NA SPG in a state-of-the-art, high resolution, coupled climate model. This mode has a period of 17 years and explains 15–30% of the annual variance in related ocean indices. It arises due to the advection of heat content anomalies around the NA SPG. Anomalous circulation drives the variability in the southern half of the NA SPG, whilst mean circulation and anomalous temperatures are important in the northern half. A negative feedback between Labrador Sea temperatures/densities and those in the North Atlantic Current is identified, which allows for the phase reversal. The atmosphere is found to act as a positive feedback on to this mode via the North Atlantic Oscillation which itself exhibits a spectral peak at 17 years. Decadal ocean density changes associated with this mode are driven by variations in temperature, rather than salinity — a point which models often disagree on and which we suggest may affect the veracity of the underlying assumptions of anomaly-assimilating decadal prediction methodologies.

Serial clustering of extratropical cyclones in a multi-model ensemble of historical and future simulations

This study has investigated serial (temporal) clustering of extra-tropical cyclones simulated by 17 climate models that participated in CMIP5. Clustering was estimated by calculating the dispersion (ratio of variance to mean) of 30 December-February counts of Atlantic storm tracks passing nearby each grid point. Results from single historical simulations of 1975-2005 were compared to those from historical ERA40 reanalyses from 1958-2001 ERA40 and single future model projections of 2069-2099 under the RCP4.5 climate change scenario. Models were generally able to capture the broad features in reanalyses reported previously: underdispersion/regularity (i.e. variance less than mean) in the western core of the Atlantic storm track surrounded by overdispersion/clustering (i.e. variance greater than mean) to the north and south and over western Europe. Regression of counts onto North Atlantic Oscillation (NAO) indices revealed that much of the overdispersion in the historical reanalyses and model simulations can be accounted for by NAO variability. Future changes in dispersion were generally found to be small and not consistent across models. The overdispersion statistic, for any 30 year sample, is prone to large amounts of sampling uncertainty that obscures the climate change signal. For example, the projected increase in dispersion for storm counts near London in the CNRMCM5 model is 0.1 compared to a standard deviation of 0.25. Projected changes in the mean and variance of NAO are insufficient to create changes in overdispersion that are discernible above natural sampling variations.

Extratropical cyclones in a warmer, moister climate: a recent Atlantic analogue

Current climate model projections do not exhibit a large change in the intensity of extratropical cyclones. However, there are concerns that current models represent moist processes poorly, and this provides motivation for investigating observational evidence for how cyclones behave in warmer climates. In the North Atlantic in particular, recent decades provide a clear contrast between warm and cold climates due to Atlantic Multidecadal Variability. In this paper we investigate these periods as analogues which may provide a guide to future cyclone behavior. While temperature and moisture rise in recent warm periods as in the projections, differences in energetics and temperature gradients imply that these periods are only partial analogues. The main result from current reanalyses is that while increased cyclone-associated precipitation is seen in the recent warm periods, there is no robust evidence of an increase in cyclone intensity by other measures, such as maximum wind speed or vorticity. A set of low- and high-resolution model simulations are also studied, suggesting that changes in cyclone intensity may be different in higher-resolution reanalyses.

On the speed of the eddy-driven jet and the width of the Hadley cell in the Southern Hemisphere

A strong correlation between the speed of the eddy-driven jet and the width of the Hadley cell is found to exist in the Southern Hemisphere, both in reanalysis data and in twenty-first-century integrations from the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report multimodel archive. Analysis of the space–time spectra of eddy momentum flux reveals that variations in eddy-driven jet speed are related to changes in the mean phase speed of midlatitude eddies. An increase in eddy phase speeds induces a poleward shift of the critical latitudes and a poleward expansion of the region of subtropical wave breaking. The associated changes in eddy momentum flux convergence are balanced by anomalous meridional winds consistent with a wider Hadley cell. At the same time, faster eddies are also associated with a strengthened poleward eddy momentum flux, sustaining a stronger westerly jet in midlatitudes. The proposed mechanism is consistent with the seasonal dependence of the interannual variability of the Hadley cell width and appears to explain at least part of the projected twenty-first-century trends.

The relationship between the ITCZ and the Southern Hemispheric eddy-driven jet

We study the effect of a thermal forcing confined to the midlatitudes of one hemisphere on the eddy-driven jet in the opposite hemisphere. We demonstrate the existence of an “interhemispheric teleconnection,” whereby warming (cooling) the Northern Hemisphere causes both the intertropical convergence zone (ITCZ) and the Southern Hemispheric midlatitude jet to shift northward (southward). The interhemispheric teleconnection is effected by a change in the asymmetry of the Hadley cells: as the ITCZ shifts away from the Equator, the cross-equatorial Hadley cell intensifies, fluxing more momentum toward the subtropics and sustaining a stronger subtropical jet. Changes in subtropical jet strength, in turn, alter the propagation of extratropical waves into the tropics, affecting eddy momentum fluxes and the eddy-driven westerlies. The relevance of this mechanism is demonstrated in the context of future climate change simulations, where shifts of the ITCZ are significantly related to shifts of the Southern Hemispheric eddy-driven jet in austral winter. The possible relevance of the proposed mechanism to paleoclimates is discussed, particularly with regard to theories of ice age terminations.

Objective determination of the extratropical transition of tropical cyclones in the Northern Hemisphere

Extratropical transition (ET) has eluded objective identification since the realisation of its existence in the 1970s. Recent advances in numerical, computational models have provided data of higher resolution than previously available. In conjunction with this, an objective characterisation of the structure of a storm has now become widely accepted in the literature. Here we present a method of combining these two advances to provide an objective method for defining ET. The approach involves applying K-means clustering to isolate different life-cycle stages of cyclones and then analysing the progression through these stages. This methodology is then tested by applying it to five recent years from the European Centre of Medium-Range Weather Forecasting operational analyses. It is found that this method is able to determine the general characteristics for ET in the Northern Hemisphere. Between 2008 and 2012, 54% (±7, 32 of 59) of Northern Hemisphere tropical storms are estimated to undergo ET. There is great variability across basins and time of year. To fully capture all the instances of ET is necessary to introduce and characterise multiple pathways through transition. Only one of the three transition types needed has been previously well-studied. A brief description of the alternate types of transitions is given, along with illustrative storms, to assist with further study

Links between tropical Pacific seasonal, interannual and orbital variability during the Holocene

The El Niño/Southern Oscillation (ENSO) is the leading mode of interannual climate variability. However, it is unclear how ENSO has responded to external forcing, particularly orbitally induced changes in the amplitude of the seasonal cycle during the Holocene. Here we present a reconstruction of seasonal and interannual surface conditions in the tropical Pacific Ocean from a network of high-resolution coral and mollusc records that span discrete intervals of the Holocene. We identify several intervals of reduced variance in the 2 to 7 yr ENSO band that are not in phase with orbital changes in equatorial insolation, with a notable 64% reduction between 5,000 and 3,000 years ago. We compare the reconstructed ENSO variance and seasonal cycle with that simulated by nine climate models that include orbital forcing, and find that the models do not capture the timing or amplitude of ENSO variability, nor the mid-Holocene increase in seasonality seen in the observations; moreover, a simulated inverse relationship between the amplitude of the seasonal cycle and ENSO-related variance in sea surface temperatures is not found in our reconstructions. We conclude that the tropical Pacific climate is highly variable and subject to millennial scale quiescent periods. These periods harbour no simple link to orbital forcing, and are not adequately simulated by the current generation of models.

An assessment of Indian monsoon seasonal forecasts and mechanisms underlying monsoon interannual variability in the Met Office GloSea5-GC2 system

We assess Indian summer monsoon seasonal forecasts in GloSea5-GC2, the Met Office fully coupled subseasonal to seasonal ensemble forecasting system. Using several metrics, GloSea5-GC2 shows similar skill to other state-of-the-art forecast systems. The prediction skill of the large-scale South Asian monsoon circulation is higher than that of Indian monsoon rainfall. Using multiple linear regression analysis we evaluate relationships between Indian monsoon rainfall and five possible drivers of monsoon interannual variability. Over the time period studied (1992-2011), the El Nino-Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) are the most important of these drivers in both observations and GloSea5-GC2. Our analysis indicates that ENSO and its teleconnection with the Indian rainfall are well represented in GloSea5-GC2. However, the relationship between the IOD and Indian rainfall anomalies is too weak in GloSea5-GC2, which may be limiting the prediction skill of the local monsoon circulation and Indian rainfall. We show that this weak relationship likely results from a coupled mean state bias that limits the impact of anomalous wind forcing on SST variability, resulting in erroneous IOD SST anomalies. Known difficulties in representing convective precipitation over India may also play a role. Since Indian rainfall responds weakly to the IOD, it responds more consistently to ENSO than in observations. Our assessment identifies specific coupled biases that are likely limiting GloSea5-GC2 prediction skill, providing targets for model improvement.

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.

Steric sea level variability (1993-2010) in an ensemble of ocean reanalyses and objective analyses

Quantifying the effect of the seawater density changes on sea level variability is of crucial importance for climate change studies, as the sea level cumulative rise can be regarded as both an important climate change indicator and a possible danger for human activities in coastal areas. In this work, as part of the Ocean Reanalysis Intercomparison Project, the global and regional steric sea level changes are estimated and compared from an ensemble of 16 ocean reanalyses and 4 objective analyses. These estimates are initially compared with a satellite-derived (altimetry minus gravimetry) dataset for a short period (2003–2010). The ensemble mean exhibits a significant high correlation at both global and regional scale, and the ensemble of ocean reanalyses outperforms that of objective analyses, in particular in the Southern Ocean. The reanalysis ensemble mean thus represents a valuable tool for further analyses, although large uncertainties remain for the inter-annual trends. Within the extended intercomparison period that spans the altimetry era (1993–2010), we find that the ensemble of reanalyses and objective analyses are in good agreement, and both detect a trend of the global steric sea level of 1.0 and 1.1 ± 0.05 mm/year, respectively. However, the spread among the products of the halosteric component trend exceeds the mean trend itself, questioning the reliability of its estimate. This is related to the scarcity of salinity observations before the Argo era. Furthermore, the impact of deep ocean layers is non-negligible on the steric sea level variability (22 and 12 % for the layers below 700 and 1500 m of depth, respectively), although the small deep ocean trends are not significant with respect to the products spread.

Ocean heat content variability and change in an ensemble of ocean reanalyses

Accurate knowledge of the location and magnitude of ocean heat content (OHC) variability and change is essential for understanding the processes that govern decadal variations in surface temperature, quantifying changes in the planetary energy budget, and developing constraints on the transient climate response to external forcings. We present an overview of the temporal and spatial characteristics of OHC variability and change as represented by an ensemble of dynamical and statistical ocean reanalyses (ORAs). Spatial maps of the 0–300 m layer show large regions of the Pacific and Indian Oceans where the interannual variability of the ensemble mean exceeds ensemble spread, indicating that OHC variations are well-constrained by the available observations over the period 1993–2009. At deeper levels, the ORAs are less well-constrained by observations with the largest differences across the ensemble mostly associated with areas of high eddy kinetic energy, such as the Southern Ocean and boundary current regions. Spatial patterns of OHC change for the period 1997–2009 show good agreement in the upper 300 m and are characterized by a strong dipole pattern in the Pacific Ocean. There is less agreement in the patterns of change at deeper levels, potentially linked to differences in the representation of ocean dynamics, such as water mass formation processes. However, the Atlantic and Southern Oceans are regions in which many ORAs show widespread warming below 700 m over the period 1997–2009. Annual time series of global and hemispheric OHC change for 0–700 m show the largest spread for the data sparse Southern Hemisphere and a number of ORAs seem to be subject to large initialization ‘shock’ over the first few years. In agreement with previous studies, a number of ORAs exhibit enhanced ocean heat uptake below 300 and 700 m during the mid-1990s or early 2000s. The ORA ensemble mean (±1 standard deviation) of rolling 5-year trends in full-depth OHC shows a relatively steady heat uptake of approximately 0.9 ± 0.8 W m−2 (expressed relative to Earth’s surface area) between 1995 and 2002, which reduces to about 0.2 ± 0.6 W m−2 between 2004 and 2006, in qualitative agreement with recent analysis of Earth’s energy imbalance. There is a marked reduction in the ensemble spread of OHC trends below 300 m as the Argo profiling float observations become available in the early 2000s. In general, we suggest that ORAs should be treated with caution when employed to understand past ocean warming trends—especially when considering the deeper ocean where there is little in the way of observational constraints. The current work emphasizes the need to better observe the deep ocean, both for providing observational constraints for future ocean state estimation efforts and also to develop improved models and data assimilation methods.

Explaining patterns of avian diversity and endemicity: climate and biomes of southern Africa over the last 140,000 years

Aim Test hypotheses that present biodiversity and endemic species richness are related to climatic stability and/or biome persistence.Location Africa south of 15° S. Methods Seventy eight HadCM3 general circulation model palaeoclimate experiments spanning the last 140,000 years, plus a pre-industrial experiment,were used to calculate measures of climatic variability for 0.5° grid cells. Models were fitted relating distributions of the nine biomes of South Africa,Lesotho and Swaziland to present climate. These models were used to simulate potential past biome distribution and extent for the 78 palaeoclimate experiments, and three measures of biome persistence. Climatic response surfaces were fitted for 690 bird species regularly breeding in the region and used to simulate present species richness for cells of the 0.5° grid. Species richness was evaluated for residents, mobile species (nomadic or partially/altitudinally migrant within the region), and intra-African migrants, and also separately for endemic/near-endemic (hereafter ‘endemic’) species as a whole and those associated with each biome. Our hypotheses were tested by analysing correlations between species richness and climatic variability or biome persistence. Results The magnitude of climatic variability showed clear spatial patterns. Marked changes in biome distributions and extents were projected, although limited areas of persistence were projected for some biomes. Overall species richness was not correlated with climatic variability, although richness of mobile species showed a weak negative correlation. Endemic species richness was significantly negatively correlated with climatic variability. Strongest correlations, however, were positive correlations between biome persistence and richness of endemics associated with individual biomes. Main conclusions Low climatic variability, and especially a degree of stability enabling biome persistence, is strongly correlated with species richness of birds endemic to southern Africa. This probably principally reflects reduced extinction risk for these species where the biome to which they are adapted persisted

Latitudinal gradients in butterfly population variability are influenced by landscape heterogeneity

The variability of populations over time is positively associated with their risk of local extinction. Previous work has shown that populations at the high-latitude boundary of species’ ranges show higher inter-annual variability, consistent with increased sensitivity and exposure to adverse climatic conditions. However, patterns of population variability at both high- and low-latitude species range boundaries have not yet been concurrently examined. Here, we assess the inter-annual population variability of 28 butterfly species between 1994 and 2009 at 351 and 18 sites in the United Kingdom and Catalonia, Spain, respectively. Local population variability is examined with respect to the position of the species’ bioclimatic envelopes (i.e. whether the population falls within areas of the ‘core’ climatic suitability or is a climatically ‘marginal’ population), and in relation to local landscape heterogeneity, which may influence these range location – population dynamic relationships. We found that butterfly species consistently show latitudinal gradients in population variability, with increased variability in the more northerly UK. This pattern is even more marked for southerly distributed species with ‘marginal’ climatic suitability in the UK but ‘core’ climatic suitability in Catalonia. In addition, local landscape heterogeneity did influence these range location – population dynamic relationships. Habitat heterogeneity was associated with dampened population dynamics, especially for populations in the UK. Our results suggest that promoting habitat heterogeneity may promote the persistence of populations at high-latitude range boundaries, which may potentially aid northwards expansion under climate warming. We did not find evidence that population variability increases towards southern range boundaries. Sample sizes for this region were low, but there was tentative evidence, in line with previous ecological theory, that local landscape heterogeneity may promote persistence in these retracting low-latitude range boundary populations.

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.

Annular modes and apparent eddy feedbacks in the Southern Hemisphere

Lagged correlation analysis is often used to infer intraseasonal dynamical effects but is known to be affected by non-stationarity. We highlight a pronounced quasi-two-year peak in the anomalous zonal wind and eddy momentum flux convergence power spectra in the Southern Hemisphere, which is prima facie evidence for non-stationarity. We then investigate the consequences of this non-stationarity for the Southern Annular Mode and for eddy momentum flux convergence. We argue that positive lagged correlations previously attributed to the existence of an eddy feedback are more plausibly attributed to non-stationary interannual variability external to any potential feedback process in the mid-latitude troposphere. The findings have implications for the diagnosis of feedbacks in both models and re-analysis data as well as for understanding the mechanisms underlying variations in the zonal wind.