Causes of the rapid warming of the North Atlantic ocean in the mid-1990s

In the mid-1990s the subpolar gyre of the North Atlantic underwent a remarkable rapid warming, with sea surface temperatures increasing by around 1C in just 2 years. This rapid warming followed a prolonged positive phase of the North Atlantic Oscillation (NAO), but also coincided with an unusually negative NAO index in the winter of 1995/96. By comparing ocean analyses and carefully designed model experiments we show that this rapid warming can be understood as a delayed response to the prolonged positive phase of the NAO, and not simply an instantaneous response to the negative NAO index of 1995/96. Furthermore, we infer that the warming was partly caused by a surge, and subsequent decline, in the Meridional Overturning Circulation and northward heat transport of the Atlantic Ocean. Our results provide persuasive evidence of significant oceanic memory on multi-annual timescales, and are therefore encouraging for the prospects of developing skillful predictions.

A tutorial on pseudospectral methods for computational optimal control

[English] This paper is a tutorial introduction to pseudospectral optimal control. With pseudospectral methods, a function is approximated as a linear combination of smooth basis functions, which are often chosen to be Legendre or Chebyshev polynomials. Collocation of the differential-algebraic equations is performed at orthogonal collocation points, which are selected to yield interpolation of high accuracy. Pseudospectral methods directly discretize the original optimal control problem to recast it into a nonlinear programming format. A numerical optimizer is then employed to find approximate local optimal solutions. The paper also briefly describes the functionality and implementation of PSOPT, an open source software package written in C++ that employs pseudospectral discretization methods to solve multi-phase optimal control problems. The software implements the Legendre and Chebyshev pseudospectral methods, and it has useful features such as automatic differentiation, sparsity detection, and automatic scaling. The use of pseudospectral methods is illustrated in two problems taken from the literature on computational optimal control. [Portuguese] Este artigo e um tutorial introdutorio sobre controle otimo pseudo-espectral. Em metodos pseudo-espectrais, uma funcao e aproximada como uma combinacao linear de funcoes de base suaves, tipicamente escolhidas como polinomios de Legendre ou Chebyshev. A colocacao de equacoes algebrico-diferenciais e realizada em pontos de colocacao ortogonal, que sao selecionados de modo a minimizar o erro de interpolacao. Metodos pseudoespectrais discretizam o problema de controle otimo original de modo a converte-lo em um problema de programa cao nao-linear. Um otimizador numerico e entao empregado para obter solucoes localmente otimas. Este artigo tambem descreve sucintamente a funcionalidade e a implementacao de um pacote computacional de codigo aberto escrito em C++ chamado PSOPT. Tal pacote emprega metodos de discretizacao pseudo-spectrais para resolver problemas de controle otimo com multiplas fase. O PSOPT permite a utilizacao de metodos de Legendre ou Chebyshev, e possui caractersticas uteis tais como diferenciacao automatica, deteccao de esparsidade e escalonamento automatico. O uso de metodos pseudo-espectrais e ilustrado em dois problemas retirados da literatura de controle otimo computacional.

European plant phenology and climate as seen in a 20-year AVHRR land-surface parameter dataset

Vegetation distribution and state have been measured since 1981 by the AVHRR (Advanced Very High Resolution Radiometer) instrument through satellite remote sensing. In this study a correction method is applied to the Pathfinder NDVI (Normalized Difference Vegetation Index) data to create a continuous European vegetation phenology dataset of a 10-day temporal and 0.1° spatial resolution; additionally, land surface parameters for use in biosphere–atmosphere modelling are derived. The analysis of time-series from this dataset reveals, for the years 1982–2001, strong seasonal and interannual variability in European land surface vegetation state. Phenological metrics indicate a late and short growing season for the years 1985–1987, in addition to early and prolonged activity in the years 1989, 1990, 1994 and 1995. These variations are in close agreement with findings from phenological measurements at the surface; spring phenology is also shown to correlate particularly well with anomalies in winter temperature and winter North Atlantic Oscillation (NAO) index. Nevertheless, phenological metrics, which display considerable regional differences, could only be determined for vegetation with a seasonal behaviour. Trends in the phenological phases reveal a general shift to earlier (−0.54 days year−1) and prolonged (0.96 days year−1) growing periods which are statistically significant, especially for central Europe.

Initialisation and predictability of the AMOC over the last 50 years in a climate model

The mechanisms involved in Atlantic meridional overturning circulation (AMOC) decadal variability and predictability over the last 50 years are analysed in the IPSL–CM5A–LR model using historical and initialised simulations. The initialisation procedure only uses nudging towards sea surface temperature anomalies with a physically based restoring coefficient. When compared to two independent AMOC reconstructions, both the historical and nudged ensemble simulations exhibit skill at reproducing AMOC variations from 1977 onwards, and in particular two maxima occurring respectively around 1978 and 1997. We argue that one source of skill is related to the large Mount Agung volcanic eruption starting in 1963, which reset an internal 20-year variability cycle in the North Atlantic in the model. This cycle involves the East Greenland Current intensity, and advection of active tracers along the subpolar gyre, which leads to an AMOC maximum around 15 years after the Mount Agung eruption. The 1997 maximum occurs approximately 20 years after the former one. The nudged simulations better reproduce this second maximum than the historical simulations. This is due to the initialisation of a cooling of the convection sites in the 1980s under the effect of a persistent North Atlantic oscillation (NAO) positive phase, a feature not captured in the historical simulations. Hence we argue that the 20-year cycle excited by the 1963 Mount Agung eruption together with the NAO forcing both contributed to the 1990s AMOC maximum. These results support the existence of a 20-year cycle in the North Atlantic in the observations. Hindcasts following the CMIP5 protocol are launched from a nudged simulation every 5 years for the 1960–2005 period. They exhibit significant correlation skill score as compared to an independent reconstruction of the AMOC from 4-year lead-time average. This encouraging result is accompanied by increased correlation skills in reproducing the observed 2-m air temperature in the bordering regions of the North Atlantic as compared to non-initialized simulations. To a lesser extent, predicted precipitation tends to correlate with the nudged simulation in the tropical Atlantic. We argue that this skill is due to the initialisation and predictability of the AMOC in the present prediction system. The mechanisms evidenced here support the idea of volcanic eruptions as a pacemaker for internal variability of the AMOC. Together with the existence of a 20-year cycle in the North Atlantic they propose a novel and complementary explanation for the AMOC variations over the last 50 years.

Arctic Oscillation and Antarctic Oscillation in internal atmospheric variability

n this study, we investigated the features of Arctic Oscillation (AO) and Antarctic Oscillation (AAO), that is, the annular modes in the extratropics, in the internal atmospheric variability attained through an ensemble of integrations by an atmospheric general circulation model (AGCM) forced with the global observed SSTs. We focused on the interannual variability of AO/AAO, which is dominated by internal atmospheric variability. In comparison with previous observed results, the AO/AAO in internal atmospheric variability bear some similar characteristics, but exhibit a much clearer spatial structure: significant correlation between the North Pacific and North Atlantic centers of action, much stronger and more significant associated precipitation anomalies, and the meridional displacement of upper-tropospheric westerly jet streams in the Northern/Southern Hemisphere. In addition, we examined the relationship between the North Atlantic Oscillation (NAO)/AO and East Asian winter monsoon (EAWM). It has been shown that in the internal atmospheric variability, the EAWM variation is significantly related to the NAO through upper-tropospheric atmospheric teleconnection patterns.

Decadal to multidecadal variability and the climate change background

Three prominent quasi-global patterns of variability and change are observed using the Met Office's sea surface temperature (SST) analysis and almost independent night marine air temperature analysis. The first is a global warming signal that is very highly correlated with global mean SST. The second is a decadal to multidecadal fluctuation with some geographical similarity to the El Niño–Southern Oscillation (ENSO). It is associated with the Pacific Decadal Oscillation (PDO), and its Pacific-wide manifestation has been termed the Interdecadal Pacific Oscillation (IPO). We present model investigations of the relationship between the IPO and ENSO. The third mode is an interhemispheric variation on multidecadal timescales which, in view of climate model experiments, is likely to be at least partly due to natural variations in the thermohaline circulation. Observed climatic impacts of this mode also appear in model simulations. Smaller-scale, regional atmospheric phenomena also affect climate on decadal to interdecadal timescales. We concentrate on one such mode, the winter North Atlantic Oscillation (NAO). This shows strong decadal to interdecadal variability and a correspondingly strong influence on surface climate variability which is largely additional to the effects of recent regional anthropogenic climate change. The winter NAO is likely influenced by both SST forcing and stratospheric variability. A full understanding of decadal changes in the NAO and European winter climate may require a detailed representation of the stratosphere that is hitherto missing in the major climate models used to study climate change.

Midlatitude Forcing Mechanisms for Glacier Mass Balance Investigated Using General Circulation Models

A process-oriented modeling approach is applied in order to simulate glacier mass balance for individual glaciers using statistically downscaled general circulation models (GCMs). Glacier-specific seasonal sensitivity characteristics based on a mass balance model of intermediate complexity are used to simulate mass balances of Nigardsbreen (Norway) and Rhonegletscher (Switzerland). Simulations using reanalyses (ECMWF) for the period 1979–93 are in good agreement with in situ mass balance measurements for Nigardsbreen. The method is applied to multicentury integrations of coupled (ECHAM4/OPYC) and mixed-layer (ECHAM4/MLO) GCMs excluding external forcing. A high correlation between decadal variations in the North Atlantic oscillation (NAO) and mass balance of the glaciers is found. The dominant factor for this relationship is the strong impact of winter precipitation associated with the NAO. A high NAO phase means enhanced (reduced) winter precipitation for Nigardsbreen (Rhonegletscher), typically leading to a higher (lower) than normal annual mass balance. This mechanism, entirely due to internal variations in the climate system, can explain observed strong positive mass balances for Nigardsbreen and other maritime Norwegian glaciers within the period 1980–95. It can also partly be responsible for recent strong negative mass balances of Alpine glaciers.

Assessing levels of uncertainty in recent temperature time series

We examine to what degree we can expect to obtain accurate temperature trends for the last two decades near the surface and in the lower troposphere. We compare temperatures obtained from surface observations and radiosondes as well as satellite-based measurements from the Microwave Soundings Units (MSU), which have been adjusted for orbital decay and non-linear instrument-body effects, and reanalyses from the European Centre for Medium-Range Weather Forecasts (ERA) and the National Centre for Environmental Prediction (NCEP). In regions with abundant conventional data coverage, where the MSU has no major influence on the reanalysis, temperature anomalies obtained from microwave sounders, radiosondes and from both reanalyses agree reasonably. Where coverage is insufficient, in particular over the tropical oceans, large differences are found between the MSU and either reanalysis. These differences apparently relate to changes in the satellite data availability and to differing satellite retrieval methodologies, to which both reanalyses are quite sensitive over the oceans. For NCEP, this results from the use of raw radiances directly incorporated into the analysis, which make the reanalysis sensitive to changes in the underlying algorithms, e.g. those introduced in August 1992. For ERA, the bias-correction of the one-dimensional variational analysis may introduce an error when the satellite relative to which the correction is calculated is biased itself or when radiances change on a time scale longer than a couple of months, e.g. due to orbit decay. ERA inhomogeneities are apparent in April 1985, October/November 1986 and April 1989. These dates can be identified with the replacements of satellites. It is possible that a negative bias in the sea surface temperatures (SSTs) used in the reanalyses may have been introduced over the period of the satellite record. This could have resulted from a decrease in the number of ship measurements, a concomitant increase in the importance of satellite-derived SSTs, and a likely cold bias in the latter. Alternately, a warm bias in SSTs could have been caused by an increase in the percentage of buoy measurements (relative to deeper ship intake measurements) in the tropical Pacific. No indications for uncorrected inhomogeneities of land surface temperatures could be found. Near-surface temperatures have biases in the boundary layer in both reanalyses, presumably due to the incorrect treatment of snow cover. The increase of near-surface compared to lower tropospheric temperatures in the last two decades may be due to a combination of several factors, including high-latitude near-surface winter warming due to an enhanced NAO and upper-tropospheric cooling due to stratospheric ozone decrease.

Mediterranean cyclones and windstorms in a changing climate

Changes in the frequency and intensity of cyclones and associated windstorms affecting the Medi-terranean region simulated under enhanced Greenhouse Gas forcing conditions are investigated. The analysis is based on 7 climate model integrations performed with two coupled global models (ECHAM5 MPIOM and INGV CMCC), comparing the end of the twentieth century and at least the first half of the twenty-first century. As one of the models has a considerably enhanced resolution of the atmosphere and the ocean, it is also investigated whether the climate change signals are influenced by the model resolution. While the higher resolved simulation is closer to reanalysis climatology, both in terms of cyclones and windstorm distributions, there is no evidence for an influence of the resolution on the sign of the climate change signal. All model simulations show a reduction in the total number of cyclones crossing the Mediterranean region under climate change conditions. Exceptions are Morocco and the Levant region, where the models predict an increase in the number of cyclones. The reduction is especially strong for intense cyclones in terms of their Laplacian of pressure. The influence of the simulated positive shift in the NAO Index on the cyclone decrease is restricted to the Western Mediterranean region, where it explains 10–50 % of the simulated trend, depending on the individual simulation. With respect to windstorms, decreases are simulated over most of the Mediterranean basin. This overall reduction is due to a decrease in the number of events associated with local cyclones, while the number of events associated with cyclones outside of the Mediterranean region slightly increases. These systems are, however, less intense in terms of their integrated severity over the Mediterranean area, as they mostly affect the fringes of the region. In spite of the general reduction in total numbers, several cyclones and windstorms of intensity unknown under current climate conditions are identified for the scenario simulations. For these events, no common trend exists in the individual simulations. Thus, they may rather be attributed to long-term (e.g. decadal) variability than to the Greenhouse Gas forcing. Nevertheless, the result indicates that high-impact weather systems will remain an important risk in the Mediterranean Basin.

Past and recent changes in the North Atlantic oscillation

The North Atlantic oscillation (NAO) is under current climate conditions the leading mode of atmospheric circulation variability over the North Atlantic region. While the pattern is present during the entire year, it is most important during winter, explaining a large part of the variability of the large-scale pressure field, being thus largely determinant for the weather conditions over the North Atlantic basin and over Western Europe. In this study, a review of recent literature on the basic understanding of the NAO, its variability on different time scales and driving physical mechanisms is presented. In particular, the observed NAO variations and long-term trends are put into a long term perspective by considering paleo-proxy evidence. A representative number of recently released NAO reconstructions are discussed. While the reconstructions agree reasonably well with observations during the instrumental overlapping period, there is a rather high uncertainty between the different reconstructions for the pre-instrumental period, which leads to partially incoherent results, that is, periods where the NAO reconstructions do not agree even in sign. Finally, we highlight the future need of a broader definition of the NAO, the assessment of the stability of the teleconnection centers over time, the analysis of the relations to other relevant variables like temperature and precipitation, as well as on the relevant processes involved

Large-scale atmospheric dynamics of the wet winter 2009–2010 and its impact on hydrology in Portugal

The anomalously wet winter of 2010 had a very important impact on the Portuguese hydrological system. Owing to the detrimental effects of reduced precipitation in Portugal on the environmental and socio-economic systems, the 2010 winter was predominantly beneficial by reversing the accumulated precipitation deficits during the previous hydrological years. The recorded anomalously high precipitation amounts have contributed to an overall increase in river runoffs and dam recharges in the 4 major river basins. In synoptic terms, the winter 2010 was characterised by an anomalously strong westerly flow component over the North Atlantic that triggered high precipitation amounts. A dynamically coherent enhancement in the frequencies of mid-latitude cyclones close to Portugal, also accompanied by significant increases in the occurrence of cyclonic, south and south-westerly circulation weather types, are noteworthy. Furthermore, the prevalence of the strong negative phase of the North Atlantic Oscillation (NAO) also emphasises the main dynamical features of the 2010 winter. A comparison of the hydrological and atmospheric conditions between the 2010 winter and the previous 2 anomalously wet winters (1996 and 2001) was also carried out to isolate not only their similarities, but also their contrasting conditions, highlighting the limitations of estimating winter precipitation amounts in Portugal using solely the NAO phase as a predictor.

Cyclones causing wind storms in the Mediterranean: characteristics, trends and links to large-scale patterns

A climatology of cyclones with a focus on their relation to wind storm tracks in the Mediterranean region (MR) is presented. Trends in the frequency of cyclones and wind storms, as well as variations associated with the North Atlantic Oscillation (NAO), the East Atlantic/West Russian (EAWR) and the Scandinavian variability pattern (SCAND) are discussed. The study is based on the ERA40 reanalysis dataset. Wind storm tracks are identified by tracking clusters of adjacent grid boxes characterised by extremely high local wind speeds. The wind track is assigned to a cyclone track independently identified with an objective scheme. Areas with high wind activity – quantified by extreme wind tracks – are typically located south of the Golf of Genoa, south of Cyprus, southeast of Sicily and west of the Iberian Peninsula. About 69% of the wind storms are caused by cyclones located in the Mediterranean region, while the remaining 31% can be attributed to North Atlantic or Northern European cyclones. The North Atlantic Oscillation, the East Atlantic/West Russian pattern and the Scandinavian pattern all influence the amount and spatial distribution of wind inducing cyclones and wind events in the MR. The strongest signals exist for the NAO and the EAWR pattern, which are both associated with an increase in the number of organised strong wind events in the eastern MR during their positive phase. On the other hand, the storm numbers decrease over the western MR for the positive phase of the NAO and over the central MR during the positive phase of the EAWR pattern. The positive phase of the Scandinavian pattern is associated with a decrease in the number of winter wind storms over most of the MR. A third of the trends in the number of wind storms and wind producing cyclones during the winter season of the ERA40 period may be attributed to the variability of the North Atlantic Oscillation.

Implications of the North Atlantic Oscillation for a UK–Norway renewable power system

UK wind-power capacity is increasing and new transmission links are proposed with Norway, where hydropower dominates the electricity mix. Weather affects both these renewable resources and the demand for electricity. The dominant large-scale pattern of Euro-Atlantic atmospheric variability is the North Atlantic Oscillation (NAO), associated with positive correlations in wind, temperature and precipitation over northern Europe. The NAO's effect on wind-power and demand in the UK and Norway is examined, focussing on March when Norwegian hydropower reserves are low and the combined power system might be most susceptible to atmospheric variations. The NCEP/NCAR meteorological reanalysis dataset (1948–2010) is used to drive simple models for demand and wind-power, and ‘demand-net-wind’ (DNW) is estimated for positive, neutral and negative NAO states. Cold, calm conditions in NAO− cause increased demand and decreased wind-power compared to other NAO states. Under a 2020 wind-power capacity scenario, the increase in DNW in NAO− relative to NAO neutral is equivalent to nearly 25% of the present-day average rate of March Norwegian hydropower usage. As the NAO varies on long timescales (months to decades), and there is potentially some skill in monthly predictions, we argue that it is important to understand its impact on European power systems.

The transient atmospheric response to midlatitude SST anomalies

To study the transient atmospheric response to midlatitude SST anomalies, a three-layer quasigeostrophic (QG) model coupled to a slab oceanic mixed layer in the North Atlantic is used. As diagnosed from a coupled run in perpetual winter conditions, the first two modes of SST variability are linked to the model North Atlantic Oscillation (NAO) and eastern Atlantic pattern (EAP), respectively, the dominant atmospheric modes in the Atlantic sector. The two SST anomaly patterns are then prescribed as fixed anomalous boundary conditions for the model atmosphere, and its transient responses are established from a large ensemble of simulations. In both cases, the tendency of the air–sea heat fluxes to damp the SST anomalies results in an anomalous diabatic heating of the atmosphere that, in turn, forces a baroclinic response, as predicted by linear theory. This initial response rapidly modifies the transient eddy activity and thus the convergence of eddy momentum and heat fluxes. The latter transforms the baroclinic response into a growing barotropic one that resembles the atmospheric mode that had created the SST anomaly in the coupled run and is thus associated with a positive feedback. The total adjustment time is as long as 3–4 months for the NAO-like response and 1–2 months for the EAP-like one. The positive feedback, in both cases, is dependent on the polarity of the SST anomaly, but is stronger in the NAO case, thereby contributing to its predominance at low frequency in the coupled system. However, the feedback is too weak to lead to an instability of the atmospheric modes and primarily results in an increase of their amplitude and persistence and a weakening of the heat flux damping of the SST anomaly.

Mechanisms for AMOC variability simulated by the NEMO model

We have investigated mechanisms for the Atlantic Meridional Overturning Circulation (AMOC) variability at 26.5° N (other than the Ekman component) that can be related to external forcings, in particular wind variability. Resolution dependence is studied using identical experiments with 1° and 1/4° NEMO model runs over 1960–2010. The analysis shows that much of the variability in the AMOC at 26° N can be related to the wind strength over the North Atlantic, through mechanisms lagged on different timescales. At ~ 1-year lag the January–June difference of mean sea level pressure between high and mid-latitudes in the North Atlantic explains 35–50% of the interannual AMOC variability (with negative correlation between wind strength and AMOC). At longer lead timescales ~ 4 years, strong (weak) winds over the northern North Atlantic (specifically linked to the NAO index) are followed by higher (lower) AMOC transport, but this mechanism only works in the 1/4° model. Analysis of the density correlations suggests an increase (decrease) in deep water formation in the North Atlantic subpolar gyre to be the cause. Therefore another 30% of the AMOC variability at 26° N can be related to density changes in the top 1000 m in the Labrador and Irminger seas occurring ~ 4 years earlier.