European land CO2 sink influenced by NAO and East-Atlantic Pattern coupling

Abstract

Large-scale climate patterns control variability in the global carbon sink. In Europe, the North-Atlantic Oscillation (NAO) influences vegetation activity, however the East-Atlantic (EA) pattern is known to modulate NAO strength and location. Using observation-driven and modelled data sets, we show that multi-annual variability patterns of European Net Biome Productivity (NBP) are linked to anomalies in heat and water transport controlled by the NAO-EA interplay. Enhanced NBP occurs when NAO and EA are both in negative phase, associated with cool summers with wet soils which enhance photosynthesis. During anti-phase periods, NBP is reduced through distinct impacts of climate anomalies in photosynthesis and respiration. The predominance of anti-phase years in the early 2000s may explain the European-wide reduction of carbon uptake during this period, reported in previous studies. Results show that improving the capability of simulating atmospheric circulation patterns may better constrain regional carbon sink variability in coupled carbon-climate models.

ARTICLE

Received 4 Mar 2015 | Accepted 27 Nov 2015 | Published 18 Jan 2016

Ana Bastos1,2, Ivan A. Janssens3, Clia M. Gouveia2, Ricardo M. Trigo2, Philippe Ciais1, Frdric Chevallier1, Josep Penuelas4,5, Christian Rdenbeck6, Shilong Piao7, Pierre Friedlingstein8 & Steven W. Running9

Large-scale climate patterns control variability in the global carbon sink. In Europe, the North-Atlantic Oscillation (NAO) inuences vegetation activity, however the East-Atlantic (EA) pattern is known to modulate NAO strength and location. Using observation-driven and modelled data sets, we show that multi-annual variability patterns of European Net Biome Productivity (NBP) are linked to anomalies in heat and water transport controlled by the NAOEA interplay. Enhanced NBP occurs when NAO and EA are both in negative phase, associated with cool summers with wet soils which enhance photosynthesis. During anti-phase periods, NBP is reduced through distinct impacts of climate anomalies in photosynthesis and respiration. The predominance of anti-phase years in the early 2000s may explain the European-wide reduction of carbon uptake during this period, reported in previous studies. Results show that improving the capability of simulating atmospheric circulation patterns may better constrain regional carbon sink variability in coupled carbon-climate models.

1 Laboratoire des Sciences du Climat et de lEnvironnement, LSCE/IPSL, CEA-CNRS-UVSQ, Universit Paris-Saclay, F-91191 Gif-sur-Yvette, France. 2 Instituto Dom Luiz, IDL, Faculdade de Cincias, Universidade de Lisboa, Lisboa 1749-016, Portugal. 3 Department of Biology, University of Antwerp, Universiteitsplein 1, 2610 Wilrijk, Belgium. 4 CREAF, Cerdanyola del Valls, Catalonia, 08193 Barcelona, Spain. 5 CSIC, Global Ecology Unit CREAF-CSIC-UAB, Cerdanyola del Valls, Catalonia, 08193 Barcelona, Spain. 6 Max Planck Institute for Biogeochemistry, Jena 07701, Germany. 7 Department of Ecology, College of Urban and Environmental Sciences, Peking University 5 Yiheyuan Road, Haidian District, Beijing 100871, China. 8 College of Engineering, Mathematics and Physical Sciences, University of Exeter, Exeter EX4 4QF, UK. 9 Numerical Terradynamic Simulation Group, University of Montana, Missoula, Montana 59812, USA. Correspondence and requests for materials should be addressed to A.B. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms10315 OPEN

European land CO2 sink inuenced by NAO and East-Atlantic Pattern coupling

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315

Fundamental patterns of atmosphere-ocean variability are explained by teleconnections at the global1 and regional24 scales. At the global scale, the El-Nino/Southern-Oscillation

has been shown to inuence inter-annual variability in the global land carbon sink5,6 with an evident ngerprint in atmospheric CO2 growth rate uctuations7, highlighting their impact on ecosystem functioning. Hallet et al.8 proposed that teleconnections explain ecological processes even better than single climate variables, because they inuence simultaneously the range of weather variables that elicit interacting, and sometimes opposing, responses by ecosystems912.

European and North-American climate is particularly inuenced by the North-Atlantic Oscillation (NAO), which corresponds to a meridional dipole of sea-level pressure (SLP) variability between Iceland and the Azores that controls the location and strength of storm tracks in the North-Atlantic3 (Fig. 1a). NAO drives patterns and extremes in temperature, precipitation, snow cover and wind in Europe, especially during winter2,4, although its effects may propagate through the following seasons13. NAO is usually associated with warmer and wetter (colder and drier) winters in northern (southern) Europe during positive phase and approximately the reverse pattern in negative phase2. Some studies have analysed impacts of NAO on European ecosystems14,15. In the Iberian Peninsula, Carnicer et al.16 found an increase in tree defoliation in response to a period of drought (19901995) that coincided with the persistence of positive NAO. However, other teleconnections are known to also inuence the European climate, particularly the East-Atlantic (EA) pattern that has a similar conguration as NAO, albeit displaced southwards (Fig. 1a). Although EA and NAO indices are computed as independent variables (Fig. 1b), the EA pattern has been proposed to modulate the multi-decadal variability of the location and strength of the NAO dipole17. Comas-Bru and McDermott4 showed that the combined analysis of NAO and EA more appropriately describes winter climate variability in Europe.

Northern Hemisphere ecosystems are important carbon sinks18,19 that partially offset anthropogenic CO2 emissions20.

The inter-annual variability in the terrestrial carbon sink is poorly understood21, hindering better future land-sink strength projections7,22 and making regional trends in sinks difcult to detect and hence controversial23,24. Various estimates of European Net Biome Productivity (NBP) have shown that European ecosystems have been a CO2 sink during the past decades (19802005)25,26, although signs of saturation in temperate forest carbon uptake in the early 2000s have been reported27. Moreover, Piao et al.28 have found a weakening relationship between temperature and growing season vegetation greenness in Northern ecosystems, especially in Europe, during 19822011, attributable to an increase in drought conditions. Their results highlight the importance of evaluating inter-annual variability of NBP within a framework that takes into account the co-variations of the climate variables that drive ecosystem dynamics.

Here we analyse the combined impact of NAO and EA on NBP over the European area extending to the Ural Mountains (Supplementary Fig. 1), using three state-of-the-art atmospheric CO2 inversions: two versions of the Monitoring Atmospheric

Composition and ClimateInterim Implementation (MACCII) inversion29,30 and the Jena s81 v3.6 inversion31 (Supplementary Table 1). Atmospheric inversions are known to provide consistent estimates of inter-annual variability of the land sink on large aggregated regions19. The uncertainty of the annual anomalies is expected to be lower than that of the long-term mean, since these errors on annual uxes are positively correlated from 1 year to the next, and Europe is one of the regions with most robust inter-annual variability estimates32. As inversions present increasing uncertainty for smaller scale regions33 we use a complementary bottom-up estimation of NBP from eleven Dynamic Global Vegetation Models (DGVMs)20 (Supplementary Table 2). The DGVMs were forced with climate observations and include land-use changes (LUCs), but differ on a number of characteristics, which makes them independent estimations of NBP. These data sets now report more than 30 years of CO2 ux variability, a timely opportunity to revisit the relationship between carbon uxes and climate signals. Our results show that the combined NAOEA variability explains variations in the European CO2 sink from 1982 to 2012, through their control on heat and water-vapour transport towards the continent and consequent climate anomalies. We nd that the European sink is enhanced only when the two modes are in negative phase due to enhanced photosynthesis, while for both anti-phase combinations NBP is kept below average through different responses of photosynthesis and respiration.

ResultsNet biome production variability. European NBP from inversions conrms that European ecosystems were on average a net CO2 sink between 1982 and 2012 (Supplementary Fig. 2), in line with previous estimates19,25, with strong inter-annual and multi-year variability (Fig. 2a). The late 1980s and early 1990s were characterized by generally positive NBP anomalies, that is, CO2 uptake above average. A strong sink was observed during the years 19961997, followed by a decade of mostly below-average NBP, consistent with the satellite observations of Normalized Difference Vegetation Index (NDVI) data, a proxy for Gross Primary Productivity (GPP)27. From 2009 onwards, results from inversions suggest that European ecosystems became a stronger C sink again. The set of DGVMs suggests a weaker and less-variable European sink than the inversions (Supplementary Fig. 2), although the average NBP anomalies from these models

a

NAO

30.2%

EA

80 N

60 N

40 N

20 N

80 N

60 N

40 N

20 N

18.6%

70 W

30 E

70 W

30 E

50 W

10 E

50 W

10 E

30 W

10 W

30 W

10 W

0.04 0.04

0.02 0.02

0

b

1

0

1

2

NAOw EAw

1985 1990 1995 2000 2005 2010

Figure 1 | Main large-scale atmospheric circulation patterns in Europe. (a) Spatial patterns of the two rst components of 500 mb geopotential height (from NCEP/DOEII Reanalysis) variability over the North-Atlantic: NAO and EA teleconnections. (b) Time series of the winter composites of NAO (blue line) and EA (red line) from NOAA; circle markers indicate the positive and negative phases of each mode, dened by the upper and lower tercile thresholds (Supplementary Table 3).

2 NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315 ARTICLE

a

0.8

0.6

0.4

0.4

0.6

0.8

PgC per year

0.2

0.2

0

0.2

Inversions TRENDY Weaker sink

1985 1990 1995 2000 2005 2010

b

Stronger sink

Stronger sink

MACCII v11.2 MACCII v13.1 Jena s81 v3.6

CLM4.5 ISAM JSBACH JULES3.2 LPJ LPJ-GUESS LPX-Bern OCN ORCHIDEE VEGAS VISIT

0.15

0.15

0.1

0.1

PgC per year

0.05

0.05

0

Weaker sink

0.2

NAO+EA+ NAOEA NAO+EA NAOEA+

Figure 2 | Impact of NAOEA variability on NBP. (a) Annual NBP anomalies integrated over Europe (bold lines) and corresponding spread, that is, the minimum and maximum NBP estimated by each set (shaded area), for the atmospheric transport model inversions (black) and the DGVMs from TRENDY project (green). (b) Average European NBP anomalies for the four NAOEA composites (Supplementary Table 3), assessed by inversions (greyscale lines) and DGVMs (colour lines).

remain in most years within the uncertainty range of the inversions and also present no signicant trend (Fig. 2a).

Despite the lack of long-term NBP trend at continental scale, Supplementary Fig. 3 shows distinct temporal patterns in annual NBP anomalies in the four large regions encompassing most of continental Europe. The continental integral is mainly inuenced by the dynamics in central Europe (33% and 34% of variance explained on average for inversions and DGVMs, respectively) and western Russia (29%, inversions). A principal component analysis (PCA) applied on NBP elds from inversions and DGVMs conrms the importance of central Europe and western Russia dynamics, by identifying a predominant variability pattern (explaining about 30% of spatio-temporal variance) with two centres located in these regions (Supplementary Fig. 4). Particularly, anomalies in central Europe dominate the multi-year pattern observed at the continental scale, while in western Russia NBP presents an increasing trend in the rst decade followed by a peak in 19961997 and a stalling onwards (Supplementary Fig. 3). However, it must be kept in mind that atmospheric inversions do not perform as well at the regional scale as at the continental scale due to the limited, albeit still greater than in most regions, number of observation sites in Europe and to limitations of transport models.

Impact of NAO and EA on NBP. A large spread in NBP response to climate is expected from the natural variability of the climate conditions (for example, NAO years differ in NAO strength and therefore in its impacts), as well as from carryover effects from previous years climate conditions. For example, warm and moist conditions in 1 year promote organic matter decomposition, thereby increasing nutrient availability with impact on plant growth during subsequent years12. Because of

this large spread of the NBP response to current years climate, we here analyse the multi-year average NBP values within the four NAOEA phase composites (NAO EA , NAO EA ,

NAO EA and NAO EA ), rather than individual years (Fig. 2b; Supplementary Table 3).

Our analysis reveals that the sole use of NAO does not sufce to explain NBP variability, because its impacts are strongly modulated by the sign of the EA. Inversions and DGVMs agree on signicant differences in NBP composite anomalies during NAO years (Supplementary Tables 4 and 5), depending on the state of the EA. NBP is enhanced when EA is also negative, and is weaker than average when EA is positive (Fig. 2a; compare NAO EA with NAO EA ). In fact, both sets of NBP estimates identify NAO EA as the single combination with consistent estimates of NBP enhancement (except two DGVMs), even for an extended record from 1950 to 2012 (Supplementary Fig. 5). During the other phase combinations, inversions and DGVMs report NBP anomalies of similar magnitude for the two anti-phase states (NAO EA and NAO EA ), which are signicantly lower than for NAO EA . For NAO EA ,

MACCII inversions disagree with DGVMs and Jena inversion on the sign of the anomalies. However, the limited number of years (only four) in this composite highlights the need for caution in the interpretation of NBP anomalies for this combination. It is worth pointing that inversions present different sensitivities to NAOEA variability: while the sign of NBP anomaly from the two MACCII versions tends to depend on the phase/anti-phase dichotomy, NBP anomalies from the Jena inversion appear to respond more strongly to EA phases.

Mechanism. Having shown the need to consider both NAO and EA when assessing European CO2 uptake variability, we now propose the main components of the physical mechanisms associated, that is, linking variations in atmospheric circulation and climate variables to the anomalies observed in NBP for each NAOEA phase combination.

Variability in the land sink is linked to that in the teleconnections through the control of the latter on one or more of the meteorological variables inuencing NBP. The four NAOEA phases are characterized by distinct winter SLP anomaly patterns (DecemberApril (DecApr); Fig. 3a, black contours) and 500 mb geopotential height anomalies (Fig. 3b, black contours) during late winter/early spring. Although NAO clearly controls the sign of the North/South SLP anomaly dipole, with NAO (NAO ) imposing negative (positive)

anomalies in the North and positive (negative) in the south, the EA modulates their strength and conguration, consistent with (ref. 4). During in-phase combinations, a sharp North/South pressure gradient is observed at the surface and in altitude, while during anti-phase combinations anomalies tend to be less intense and present a meandering conguration.

By inuencing the strength and shape of the North-Atlantic jet, the interplay between NAO and EA controls the atmospheric transport of heat and humidity over the region (Fig. 3a,b; colours). The stronger pressure gradients during in-phase combinations create a channel of enhanced (weaker) eastward transport of heat and water-vapour spanning from the UK to eastern Europe and including Scandinavia, and reduced (higher) heat and moisture advection over the Mediterranean during NAO (NAO ). In contrast, during anti-phase periods, the weaker pressure gradients at the surface and in altitude leads to a northward displacement and attenuation of eastward heat ux, and water-vapour transport is close to average values in most of the continental Europe. The impact of these uxes on cloud cover was evaluated as a means to estimate the combined variability of

NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315

a

NAO+EA+ NAOEA NAO+EA NAOEA+

20

20

10

10

0

VT (g.km1 s1 )

HT (MW km1 )

52.5 0 2.5 5

52.5 0 2.5 5

Cloud cover (%)

b

c

Figure 3 | Impact of NAOEA variability on atmospheric circulation. Anomaly elds during extended winter (DecApr) for each of the four NAOEA phases of (a) sea-level pressure (black contours, in mb) and eastward water-vapour transport (colours); (b) 500-mb geopotential height (black contours, in g.p.m.) and heat transport (colours); (c) Cloud cover fraction.

precipitation and sunlight, both of which are relevant (with opposing effects) for vegetation growth during these months9. The patterns of winter cloud cover anomalies roughly match the ones of heat and vapour transport, with increased cloudiness in the regions affected by stronger heat and water ux anomalies. It must be noted that although precipitation anomalies are expected to be related to cloud cover, the dependency is not direct, since precipitation also depends on vertical instability of the atmosphere.

The differences in regional advection of heat and humidity transport associated with a given phase impose distinct and spatially heterogeneous climate conditions during late winter/early spring, which may propagate until summer (MaySeptember) due to the memory effects of snow cover13 and soil moisture34 (Fig. 4). During NAO EA , the inhibition of heat and water-vapour transport in winter establishes very cold conditions over most of Europe, which lead to more precipitation falling as snow rather than rain, explaining the negative soil moisture anomalies and the generally higher-than-average snow depth registered in most areas. In the Mediterranean region, soil moisture is enhanced during DecemberApril due to the stronger advection of water vapour and prevalence of cloudiness (Fig. 3). In contrast to the lower soil moisture during winter/spring, soil water was enhanced during spring/summer in the NAO EA

composite in most of central Europe (Fig. 4; Supplementary Table 6). This shift in water availability may be related to later snow melt and lower evaporative losses during the colder winter/spring periods. Higher soil moisture levels support stronger latent heat exchange and thus evaporative cooling in most regions, except the Iberian Peninsula, Scandinavia and the southern section of western Russia, in which dry conditions promote higher summer temperatures. For NAO EA these patterns are approximately the opposite.

The meandering SLP pattern and the corresponding energy transport during anti-phase periods impose a north-east/southwest gradient (rather than north/south as during in-phase combinations) in temperature anomalies with warmer-than-average (colder) winter/spring conditions over Scandinavia and western Russia and cooler (warmer) temperatures in the western-Mediterranean sector during NAO EA (NAO EA ).

The patterns of soil moisture for NAO EA (NAO EA )

reect the northward displacement and attenuation of water-vapour transport anomalies: most regions register close to average winter soil moisture, with drier (wetter) winter conditions being registered only in western Russia and Iberian Peninsula. Most of central Europe registers lower (higher) levels of snow cover, followed by dry (wet) soil conditions in summer. As for NAO EA , the snow-depth winter patterns generally match the ones observed in summer soil moisture and temperature anomalies, although the moisturetemperature coupling is most evident for NAO EA . The Iberian peninsula is also affected by very low winter cloud cover during NAO

EA , which promotes higher evaporation rates, explaining the strong water decits in winter.

The anomalies in NBP are related to two largely offsetting uxes, photosynthesis and respiration, which respond differently to climate anomalies12. It is worth comparing the continental NBP response to NAOEA phases with observations of vegetation greenness anomalies assessed by annually integrated NDVI elds (Fig. 5a) from the Global Inventory Modeling and Mapping Studies (GIMMS)35, as well as the corresponding seasonal evolution (Fig. 5b, black lines). Since DGVMs allow partitioning the NBP anomalies into GPP and respiration responses, the seasonal dynamics of these two components of NBP is also evaluated (Fig. 5b, coloured lines). It must be noted that since the climate anomalies are spatially heterogeneous, and Europe is composed of diverse biomes, these responses must also be assessed at the regional scale (Supplementary Figs 6 and 7, Supplementary Table 4). Regional NBP anomalies from inversions and DGVMs are compared, although inversions present higher uncertainty for smaller regions and in most cases do not present signicant results. At the regional scale, DGVMs are expected to better capture ecological variability patterns in response to climate anomalies (see Methods).

Inversions and most DGVMs consistently estimate enhanced NBP during NAO EA years in all regions, except in the

Iberian Peninsula, where inversions and DGVMs diverge. The clear positive NBP anomaly during NAO EA cannot be explained solely by positive NDVI anomalies, since large areas in Europe present NDVI anomalies of small magnitude, and in many regions even slight browning. On the seasonal scale, the cold and dry winter and spring conditions affecting most of

4 NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315 ARTICLE

a

NAO+EA+ NAOEA NAO+EA NAOEA+

SW anom (%) T anom (C) T anom (C)

DecApr

MaySep

2 1 0 1 2

b

1

0

0.5

0.5

1

c

DecApr MaySep

d

2 1 0 1 2

e

2 1 0 1 2

DecApr

SD anom (cm)

Figure 4 | Impact of NAOEA variability on climate variables. Average anomalies in late winter/early spring (DecApr) and summer (MaySep) of temperature (T, (a,b)), soil water content (SW, (c,d)) and snow depth (SD, only DecApr, (e)) for the four NAOEA composites.

Europe impose negative NDVI anomalies and lower-than-average GPP (especially in western Russia; Fig. 5b; Supplementary Fig. 6). High water availability due to retarded snow melt in spring and reduced evaporative soil-water loss during the colder spring sustains increased GPP from late spring to early autumn, while respiration is kept close to average values due to the mild summer temperatures, especially in parts of central Europe and western Russia. Indeed, most DGVMs estimate a strong positive dependence of summer GPP in soil water in all regions except Scandinavia (Supplementary Table 7; Supplementary Fig. 8), while only four are able to capture the known dependence of GPP in winter water availability in the Iberian Peninsula36. The only two DGVMs that estimate negative NBP anomalies at continental scale for NAO EA (ISAM and VEGAS) present very weak control of water availability on GPP. The enhanced CO2 uptake during the NAO EA phase combination is thus mostly due to enhanced GPP during the growing season, not to reduced respiration.

The highest positive integrated NDVI anomalies are observed for NAO EA in most of Europe (except Iberia). However, the higher annual NDVI does not coincide with increased

GPP, mainly because positive NDVI anomalies occur when radiation is sub-optimal (that is, late winter and early spring) and hence photosynthesis is weakly stimulated. DVGMs report higher GPP during spring, which offsets the increase in respiration during these warmer conditions. In summer, NDVI is about average, however, DGVMs estimate strong negative GPP and respiration anomalies, the former being stronger than the latter, and leading the negative NBP response (Fig. 2b). This is mainly due to the response observed in central Europe, where very dry

conditions lead to very low GPP, despite the low temperatures inducing lower-than-average respiration rates (Supplementary Figs 6 and 7).

Despite both anti-phase composites resulting in negative NBP anomalies, they diverge in the corresponding climate patterns (Fig. 4). For NAO EA , DGVMs (inversions are consistent, but not signicant) report a sink reduction in all regions except western Russia, which matches the marked vegetation browning and decrease in GPP during the year (Fig. 5) imposed by the cold and dry winter/spring in central Europe and the Mediterranean, especially in the Iberian Peninsula, and the drier summers in central Europe. This decrease in photosynthesis is accompanied by a smaller reduction in respiration, explaining the reduced NBP during this phase combination (Fig. 2; Supplementary Fig. 7). On the contrary, NAO EA is characterized by overall high annual

NDVI elds, associated with increased photosynthesis during spring (Fig. 5), promoted by warm and wet conditions over most of Europe, except in western Russia, where very cold late winter and spring inhibits vegetation activity (Supplementary Fig. 7; Supplementary Table 6). However, the GPP enhancement during spring is offset by greater increase in respiration during summer due to widespread positive temperature anomalies that reinforce, and are amplied by, soil dryness. This is especially evident in western Russia, where below average summer GPP is accompanied by enhanced respiration.

DiscussionChanges in NBP are the result of variation in land cover, in land management and/or in ecosystem activity. Changes in land cover

NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315

a b

0.1

0.1

NAO+ EA+

NAO EA

NAO+ EA

NAO EA+

g m2 per month

g m2 per month

g m2 per month

g m2 per month

0.05

0.05

0

0.01 0 0.01

NDVI

J J

M M S N

J J

0.05

0.05

0

0.1

0.1

0.01 0 0.01

NDVI

M M S N

J J

M M S N

J J

M M S N

0.05

0.05

0

0.1

0.1

0.01 0 0.01

NDVI

0.05

0.05

0

0.1

0.1

0.01 0 0.01

NDVI

0.2 0.2

0.1 0.1

0

NDVIint

GPP Ra + Rh NDVI

Figure 5 | Ecological response to NAOEA phases. (a) Anomaly elds of annually integrated NDVI from GIMMS for each of the four NAOEA composites; (b) seasonal anomalies of NDVI (black line, right yy-axis), GPP (green, left yy-axis) and aggregated autotrophic and heterotrophic respiration (red, left yy-axis). For GPP and respiration, bold line corresponds to the average of the eleven DGVMs and the shaded area indicates the model spread.

and land management are slow and not very intense in Europe, being more likely to induce long-term trends rather than short-term variation in NBP. On the inter-annual scale, NBP responds to a number of physical drivers which present co-variability patterns that are ultimately controlled by large-scale atmospheric circulation. Here, a mechanism driving part of the inter-annual and multi-year variability in the European sink is proposed, based on the interplay between the two main patterns of atmospheric circulation in the North-Atlantic sector (NAO and EA) and their control on heat and water-vapour transport across Europe.

The NBP enhancement during NAO EA conditions suggests that the highest NBP (largest uptake) should occur during periods with predominantly negative phases of both teleconnections. In fact, the largest NBP peak observed in 19961997 in the inversions corresponds to two consecutive years of NAO EA . The strong increase in NBP during 1985 also occurred during this combination of NAO and EA. In contrast, the years with lowest NBP are associated with years of opposing phases of NAO and EA: the rst 3 years of the period (19821984) and 19981999, when a strong decrease in NBP is observed. Anti-phase combinations dominated the late 1990s and the 2000s, with some exceptions, with a period of 6 consecutive years of anti-phase lasting from 1998 to 2003, coinciding with a weaker sink over the European continent11 (Fig. 2).

When computing future trends in NBP, the presence of multi-annual variability and trends in the most prominent modes of atmospheric circulation should be taken into account. In this context, anthropogenic forcing has been shown to lead to more frequent positive phases of NAO37, associated with lower NBP,

especially when EA is negative. However CMIP3 and CMIP5 models predict opposite trends in the NAO sign over the 21st century under similar scenarios38. Shepherd39 pointed out that atmospheric dynamics still constitutes an important source of uncertainty in the state-of-the-art CMIP5 models, and showed that the climate-change response of the wintertime North-Atlantic jet largely differs among models, both in magnitude and sign. Other works have shown the inuence of other teleconnections such as the Scandinavian Pattern17 or the Pacic Decadal Oscillation40 on North-Atlantic sea-level pressure, adding further complexity to the study of long-term NBP variability.

This study shows that the non-stationary relationship between teleconnections may strongly affect the relationships found between climate and ecosystem activity. This highlights the importance of correctly reproducing the mid-latitude atmospheric circulation patterns in Earth System Models to account for inter-annual and decadal variability in the land carbon sink. More long-term high-quality observation-driven data sets are needed to further improve our understanding of the multi-year variability of the land sink and its relation to climate variability.

Methods

Inversions. Atmospheric CO2 inversions estimate spatial distribution of surface CO2 uxes using a top-down approach, in which surface uxes are estimated from atmospheric CO2 concentration measurements across the globe. An atmospheric transport model that represents the global atmospheric circulation is required, as well as prior-information about the uxes. Flux estimates over land include the carbon exchange in ecosystems, LUC related uxes and forest re CO2 emissions as

6 NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315 ARTICLE

well as fossil-fuel emissions. To estimate NBP, the latter are subtracted from the ux estimates.

We use monthly CO2 uxes of two versions of the atmospheric inversion of the MACCII from the LSCE29,30. Version 11.2 covers the period from 1979 until 2011 and is used here from January 1982 onwards. It has a 2.5 3.75

(latitude, longitude) spatial resolution. Version 13.1 has a slightly higher latitudinal resolution (1.9) and covers year 2012 as well (Supplementary Table 1). Apart from the archive length, the main difference between v11.2 and v13.1 is the number of vertical layers in the underlying atmospheric transport model, the LMD General Circulation Model (LMDZ)29, that was doubled in the latter (39 versus 19), thus allowing a better representation of the variability of CO2 in the planetary boundary layer, where the assimilated measurements are located.

Since MACC inversions are based on a large but variable number of stations during the study period (134 in total), with fewer observations in the beginning of the time series, part of the variability observed on the estimated uxes may be affected by the changes in the number of sites41. Thus, we use the Jena s81 (v3.6) inversion that is based on fewer stations and is provided at 4 5 spatial resolution

(Supplementary Table 1), but uses a consistent set of sites for the 19822012 period. Jena inversion computes NBP on each grid cell using the atmospheric transport model TM3, as described in ref. 31.

To subtract fossil fuel emissions, MACCII inversions use the EDGAR3.2 FastTrack 2000 emission database42 with scaling factor to account for trends: MACC uses the annual global totals of the Carbon Dioxide Information Analysis Center and prescribes an increase of 1.4% per year after 2001. Jena inversion uses EDGAR 4.2, with FastTrack 2010 for 2009 and 2010, extrapolation based on British Petrol global totals for 2011 and 2012, and 2% year increase afterwards.

Uncertainties of annual NBP at the continental scale are given, for each model, in Supplementary Table 1. Uncertainties of the annual uxes account for errors in the prior uxes, network conguration (number of sites, Supplementary Fig. 1) and in the transport model. The uncertainty of the annual anomalies is expected to be lower than that of the absolute uxes, since these errors are positively correlated from 1 year to the next32. We use two inversions computed with the same model but variable number of sites, and a second with a xed number of observations,and use the spread between inversions as an indicator of the uncertainty in the anomalies.

In order to test the available signals from the smaller number of sites in Europe used by the Jena s81 inversion, we performed three synthetic-data runs (1982 2010) by transporting net ecosystem exchange (NEE) elds as simulated by the Biome-BGC DGVM in the atmosphere using the same transport model as in the Jena inversion (v3.7), though on coarser 10 8 spatial resolution. To create pseudo

data, modelled atmospheric CO2 is sampled at the same locations and times as the real data for three different station sets: as in Jena s81 (14 sites) and s90 (25 sites) as well as mostly the same sites used for MACCII inversion (Supplementary Table 2). The CO2 concentrations sampled in this way were inverted back using the Jena scheme (setting ocean and fossil fuel priors to zero). As the same transport model is used for synthetic data and inversion, results are not contaminated by model errors. By comparing with the original uxes (Supplementary Fig. 9), the sensitivity of the inversion skill to the site network can then be evaluated. Despite the presence of discrepancies, the run with fewer stations (Jena s81) is able to capture the variability patterns of the CO2 uxes and differences in the magnitude of the uxes are on average 0.040.05 PgC per year for each of the three runs.

Dynamic global vegetation models. DGVMs simulate the main processes of vegetation dynamics and decomposition associated with energy, water and carbon balances at the ecosystem scale and provide a bottom-up approach to evaluate NBP variability, as well as the corresponding GPP and respiration components.

The TRENDY project43 compiles outputs from a group of DGVMs to evaluate trends and drivers in landatmosphere carbon exchange20. We use 11 DGVMs (Supplementary Table 2) from simulation S3, in which all models are forced by the same values of changing atmospheric CO2 concentrations from ice core data and observations, historical climate observations from the CRU-NCEP v4 data set and LUC from the History Database of the Global Environment (HYDE) database of human-induced land-cover changes44.

Supplementary Table 2 provides a summary of the characteristics of the 11 DGVMs used in simulation S3, which is performed over the period 18602012. All models include deforestation, afforestation and regrowth, but differ in the way they represent disturbances (re), nutrient limitation and do not realistically simulate land-management practices or crop seasonality, which strongly interplays with climate in determining NBP of European ecosystems. Nevertheless, DGVMs provide an independent data set to evaluate the variability of the land carbon sink. NBP from DGVMs was selected for two periods: the common period with the inversions and satellite data (19822012) for the main analysis, and the common period with the NAO and EA indices (19502012), except for JSBACH whose record ends in 2005. To partition NBP, continental and regional GPP and respiration (computed as the sum of autotrophic and heterotrophic respiration components) anomalies from the DGVMs were also evaluated between 1982 and 2012.

NDVI. Biweekly NDVI elds over 19822012 from the GIMMS NDVI data set, provided at 8 km spatial resolution35. NDVI values were integrated over each year on a pixel basis from the biweekly elds. NDVI anomalies were then computed as

the departure from the average annual integral values, and used as a proxy for GPP anomalies28.

Monthly NBP values from inversions and DGVMs were rst deseasonalized (mean seasonal-cycle removed) on a pixel basis to compute monthly anomaly and integrated over the year for annual NBP anomaly elds, for the European continent and for the regional boxes (Supplementary Fig. 3). A positive (negative) sign of NBP anomalies correspond to either an enhanced (weaker) sink or a smaller (increased) source. To evaluate the main pattern of NBP variability, a PCA was performed on NBP elds for inversions and DGVMs. Since results from PCA depend on the resolution of the data set, all elds were resampled to a common 1 spatial resolution. The rst principal component of NBP variability (PC1), the corresponding empirical orthogonal function (EOF1) and explained variance were selected for each inversion and model. EOF1 presented in all cases dipolar pattern. The consistency of these patterns was compared by locating the corresponding centres of action for each data set (Supplementary Fig. 4).

All climate variables used for this analysis were extracted from ERA-Interim reanalysis45: vertical integral of eastward transport water vapour (VT, ing km 1 s 1) and heat (HT, in MW km 1), fraction of cloud cover (%), SLP(in mb), 500-mb geopotential height (z500, in geopotential metres, g.p.m.), average air temperature at 2 m (T, in C), volumetric soil-water content (SW, in % of volume) in the top layer46 and average snow depth (SD, in cm). The data comprise monthly average elds from 1982 to 2012, at 0.75 spatial resolution.

As complementary information about soil-water conditions, monthly values of the self-calibrated Palmer Drought Severity Index (PDSI)47, provided by the NOAA/OAR/ESRL PSD, Boulder, CO, USA at 2.5 spatial resolution was used. Since PDSI presented regional dependence (negative for Iberian Peninsula and central Europe and positive for western Russia and Scandinavia), we estimated drought conditions as the PDSI departure from the regional average (Supplementary Table 6).

Atmospheric circulation. We focus our analysis on the impacts of winter state Dec-Feb (DJF) of two teleconnections, using the monthly indices from NOAAs Climate Prediction Centre, which are calculated by rotated PCA of monthly means of 500-mb geopotential height anomalies from NCEP/DOE II (ref. 48) as described in (ref. 49). This procedure allows the calculation of orthogonal (that is, non-correlated) indices for each month. The rst mode corresponds to the NAO pattern, the leading mode of SLP variability in the North Atlantic and the main climate teleconnection affecting European weather. The second mode corresponds to EA circulation pattern. These modes affect more signicantly European weather during winter2,50. For visualization of the spatial patterns associated with the two teleconnections used in this work, we performed a PCA on standardized monthly means of 500-mb geopotential height from NCEP/DOEII Reanalysis over the North-Atlantic region between 20 to 80 N and 90 W to 50 E for the winter months (DJF) as in ref. 49.

The phases of an index are dened as those exceeding the lower (negative) or upper (positive) terciles, and identied in superscript (for example, NAO for a positive phase of winter NAO). To evaluate the combined impacts of both teleconnections, we identify four subsets of years corresponding to the four possible NAOEA phase combinations, together with a neutral composite (Supplementary Table 3).

Regional uxes. Since data sets are provided at very different spatial resolutions (Supplementary Tables 1 and 2), and because inversions perform better on larger scales, we consider regionally integrated NBP from inversions and DGVMs for the European continent dened in the TransCom inter-comparison project51 (Supplementary Fig. 1). We additionally dene four large regions encompassing most of Europe: Iberian Peninsula extending from 11 to 3.5 E and 34 to 44

N; continental central Europe encompassing the region between 5 to 25 E and

44 to 53 N, but excluding Great Britain; western Russia and eastern Europe extending from 29 to 60 E and 46 to 62 N; Scandinavia covering the region between 4.5 to 29 E and 56 to 71 N (Supplementary Fig. 1).

To evaluate the skill of inversions in capturing regional uxes, we used the same synthetic runs from Jena inversion to evaluate: (i) the ability of the inversion to reproduce uxes in large regions over Europe; (ii) the inuence of the observational network. Despite the very coarse resolution, results (Supplementary Fig. 10) indicate that inversions do have moderate skill in distinguishing sub-regions within Europe. As expected from the better observational constrain, the skill increases with increasing number of stations. Nevertheless, even the s81 run (with few sites) is able to capture part of the regional differences, especially in eastern, central and western Europe. In south-western Europe (corresponding to the Iberian Peninsula), the inversion is not able to reproduce the uxes, as expected from the very low observational constraint (only two sites in the complete MACCII set, located in the Pyrenees).

NBP response to NAOEA. All data elds were selected for the study region (Supplementary Fig. 1). Continental and regional NBP anomalies were calculated for each inversion and DGVM and evaluated for each of the four NAOEA combinations (Fig. 2). In the case of MACC v11.2 inversion and JSBACH, the data sets do not encompass the whole period, therefore, the composites were calculated

NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315

with fewer years. To partition annual NBP anomalies, continental and regional GPP and total respiration (autotrophic plus heterotrophic respiration components) were calculated at monthly scale, and compared with the corresponding NDVI anomalies (Supplementary Fig. 7).

For each composite, a one-sided analysis of variance (ANOVA) was performed to test the signicance of the average anomaly value on the continental and regional scales (Supplementary Table 4), separately for inversions and DGVMs. Absolute NBP anomalies for each composite are easier to associate with NBP variations from year to year. However, it is worth assessing whether the inversions (or DGVMs) agree on the relative response, that is, the difference in NBP anomalies between two NAOEA combinations. Therefore, a two-sided ANOVA was carried out to test the difference between NBP anomaly estimates for pairs of the four combinations (Supplementary Table 5), on the continental scale. Supplementary Tables 4 and 5 also include results of the ANOVA analyses for the neutral composite. To assess whether relationships hold for longer periods, the same analysis was extended for the period 19502012 using DGVMs (Supplementary Fig. 5).

Since the synthetic runs of Jena inversion present some dependence of absolute uxes in the observational network, it is worth assessing how the number of sites may inuence the results of NBP anomalies in response to NAOEA variations. As the observational network has increased with time, Jena inversion is provided for different periods using an increasing number (always constant for each data set) of sites. Here we compare the other two inversions from Jena v3.6 that still encompass at least 20 years and use the same model but more sites: s85 (19852012, 19 sites), s90 (19902012, 25 sites). We also performed a run with the newer version of Jena inversion (v3.7) using the MACCII sites, at the same resolution as Jena s81 v3.6. The corresponding anomalies for each NAOEA combination are presented in Supplementary Fig. 11 and show that, despite results depending on the size of the observational network, the response of continental NBP to the phases of NAO and EA for the data sets with more observations is consistent with the anomalies estimated by the two MACCII inversions and Jena s81 v3.6. In most cases, the inversion run using the MACCII sites is closer to the anomalies estimated by MACCII, highlighting the dependence of the anomalies on the observational constrains, but also the importance of other sources of uncertainty, for example, differences in the transport models.

References

1. McPhaden, M. J., Zebiak, S. E. & Glantz, M. H. ENSO as an integrating concept in earth science. Science 314, 17401745 (2006).

2. Trigo, R. M., Osborn, T. J. & Corte-Real, J. M. The North Atlantic Oscillation inuence on Europe: climate impacts and associated physical mechanisms. Clim. Res 20, 917 (2002).

3. Scaife, A. a., Folland, C. K., Alexander, L. V., Moberg, A. & Knight, J. R. European climate extremes and the North Atlantic Oscillation. J. Clim. 21, 7283 (2008).

4. Comas-Bru, L. & McDermott, F. Impacts of the EA and SCA patterns on the European twentieth century NAO-winter climate relationship. Quart. J. R. Meteorol. Soc. 140, 354363 (2013).

5. Behrenfeld, M. J. et al. Biospheric primary production during an ENSO transition. Science 291, 25942597 (2001).

6. Bastos, A., Running, S. W., Gouveia, C. & Trigo, R. M. The global NPP dependence on ENSO: La Nina and the extraordinary year of 2011. J. Geophys. Res. Biogeosci. 118, 12471255 (2013).

7. Friedlingstein, P. & Prentice, I. Carbon-climate feedbacks: a review of model and observation based estimates. Curr. Opin. Environ. Sustainability 2, 251257 (2010).

8. Hallett, T. B. et al. Why large-scale climate indices seem to predict ecological processes better than local weather. Nature 430, 7175 (2004).

9. Nemani, R. R. et al. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science 300, 15601563 (2003).

10. Angert, A. et al. Drier summers cancel out the CO2 uptake enhancement induced by warmer springs. Proc. Natl Acad. Sci. USA 102, 1082310827 (2005).

11. Piao, S. et al. Net carbon dioxide losses of northern ecosystems in response to autumn warming. Nature 451, 4952 (2008).

12. Xia, J. et al. Terrestrial carbon cycle affected by non-uniform climate warming. Nat. Geosci. 7, 173180 (2014).

13. Barriopedro, D., Garca-Herrera, R. & Hernndez, E. The role of snow cover in the Northern Hemisphere winter to summer transition. Geophys. Res. Lett. 33, L14708 (2006).

14. Gouveia, C., Trigo, R. M., DaCamara, C. C., Libonati, R. & Pereira, J. M. C. The North Atlantic Oscillation and European vegetation dynamics. Int. J. Clim. 28, 18351847 (2008).

15. Peters, W. et al. Seven years of recent European net terrestrial carbon dioxide exchange constrained by atmospheric observations. Global Change Biol. 16, 13171337 (2010).

16. Carnicer, J. et al. Widespread crown condition decline, food web disruption, and amplied tree mortality with increased climate change-type drought. Proc. Natl Acad. Sci. USA 108, 14741478 (2011).

17. Moore, G. W. K., Renfrew, I. A. & Pickart, R. S. Multidecadal mobility of the North Atlantic Oscillation. J. Clim. 26, 24532466 (2013).

18. Pan, Y. et al. A large and persistent carbon sink in the worlds forests. Science 333, 988993 (2011).

19. Peylin, P. et al. Global atmospheric carbon budget: results from an ensemble of atmospheric CO2 inversions. Biogeosciences 10, 66996720 (2013).

20. Le Qur, C. et al. The global carbon budget 1959-2011. Earth Syst. Sci. Data 5, 165185 (2013).

21. Heimann, M. & Reichstein, M. Terrestrial ecosystem carbon dynamics and climate feedbacks. Nature 451, 289292 (2008).

22. Cox, P. M. et al. Sensitivity of tropical carbon to climate change constrained by carbon dioxide variability. Nature 494, 341344 (2013).

23. Zhao, M. & Running, S. W. Drought-induced reduction in global terrestrial net primary production from 2000 through 2009. Science 329, 940943 (2010).

24. Ahlstrm, A., Miller, P. A. & Smith, B. Too early to infer a global NPP decline since 2000. Geophys. Res. Lett. 39, L15403 (2012).

25. Janssens, I. A. et al. Europes terrestrial biosphere absorbs 7 to 12% of European anthropogenic CO2 emissions. Science 300, 15381542 (2003).

26. Schulze, E. D. et al. Importance of methane and nitrous oxide for Europes terrestrial greenhouse-gas balance. Nat. Geosci. 2, 842850 (2009).

27. Piao, S. et al. Changes in satellite-derived vegetation growth trend in temperate and boreal Eurasia from 1982 to 2006. Global Change Biol. 17, 32283239 (2011).

28. Piao, J. et al. Evidence for a weakening relationship between interannual temperature variability and northern vegetation activity. Nat. Commun. 5, 5018 (2014).

29. Chevallier, F. et al. CO2 surface uxes at grid point scale estimated from a global 21 year reanalysis of atmospheric measurements. J. Geophys. Res. 115, D21307 (2010).

30. Chevallier, F. et al. Toward robust and consistent regional CO2 ux estimates from in situ and spaceborne measurements of atmospheric CO2. Geophys. Res. Lett. 41, L058772 (2014).

31. Rdenbeck, C. Estimating CO2 sources and sinks from atmospheric mixing ratio measurements using a global inversion of atmospheric transport. Technical Report No. 6 (Max Planck Institute for Biogeochemistry, Jena, 2005) Available at http://www.bgc-jena.mpg.de/mpg/websiteBiogeochemie/Publikationen/Technical_Reports/tech_report6.pdf

Web End =http://www.bgc-jena.mpg.de/mpg/websiteBiogeochemie/Publikationen/ http://www.bgc-jena.mpg.de/mpg/websiteBiogeochemie/Publikationen/Technical_Reports/tech_report6.pdf

Web End =Technical_Reports/tech_report6.pdf .

32. Baker, D. F. et al. TransCom 3 inversion intercomparison: impact of transport model errors on the interannual variability of regional CO2 uxes, 1988-2003.

Global Biogeochem. Cycles 20, GB1002 (2006).33. Kaminski, T. & Heimann, M. Inverse modeling of atmospheric carbon dioxide uxes. Science 294, 259259 (2001).

34. Quesada, B., Vautard, R., Yiou, P., Hirschi, M. & Seneviratne, S. I. Asymmetric European summer heat predictability from wet and dry southern winters and springs. Nat. Clim. Change 2, 736741 (2012).

35. Tucker, C. et al. An extended AVHRR 8km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int. J. Remote Sens. 26, 44854498 (2005).

36. Penelas, J. et al. Complex spatiotemporal phenological shifts as a response to rainfall changes. New Phytol. 161, 837846 (2004).

37. Gillett, N. P., Zwiers, F. W., Weaver, A. J. & Stott, P. A. Detection of human inuence on sea-level pressure. Nature 422, 292294 (2003).

38. Cattiaux, J., Douville, H. & Peings, Y. European temperatures in CMIP5: origins of present-day biases and future uncertainties. Clim. Dyn. 41, 28892907 (2013).

39. Shepherd, T. G. Atmospheric circulation as a source of uncertainty in climate change projections. Nat. Geosci 7, 703708 (2014).

40. Trenberth, K. E., Fasullo, J. T., Branstator, G. & Phillips, A. S. Seasonal aspects of the recent pause in surface warming. Nat. Clim. Change 4, 911916 (2014).

41. Rdenbeck, C., Houweling, S., Gloor, M. & Heimann, M. CO2 ux history 1982-2001 inferred from atmospheric data using a global inversion of atmospheric transport. Atmos. Chem. Phys. 3, 19191964 (2003).

42. Olivier, J. G. J. & Berdowski, J. J. M. The Climate System 3378 (Swets & Zeitlinger Publishers, 2001).

43. Sitch, S. et al. Trends and drivers of regional sources and sinks of carbon dioxide over the past two decades. Biogeosci. Discuss. 10, 2011320177 (2013).

44. Klein Goldewijk, K., Beusen, A., van Drecht, G. & de Vos, M. The HYDE 3.1 spatially explicit database of human-induced global land-use change over the past 12.000 years. Global Ecol. Biogeogr. 20, 7386 (2011).

45. Dee, D. P. et al. The ERA-Interim reanalysis: conguration and performance of the data assimilation system. Quart. J. R. Meteorol. Soc. 137, 553597 (2009).

46. Balsamo, G. et al. A revised hydrology for the ECMWF model: verication from eld site to terrestrial water storage and impact in the integrated forecast system. J. Hydrometeorol. 10, 623643 (2009).

47. Dai, A. Characteristics and trends in various forms of the Palmer Drought Severity Index during 1900-2008. J. Geophys. Res. Atmospheres 116, D12115 (2011).

48. Kanamitsu, M. et al. NCEP/DOE AMIP-II Reanalysis (R-2). Bull. Am. Meteorol. Soc. 83, 16311643 (2002).

49. Barnston, A. G. & Livezey, R. E. Classication, seasonality and persistence of low-frequency atmospheric circulation patterns. Mon. Weather Rev. 115, 10831126 (1987).

8 NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10315 ARTICLE

50. Hurrell, J. W. Decadal Trends in the North Atlantic Oscillation: regional temperatures and precipitation. Science 269, 676679 (1995).

51. Gurney, K. R. et al. Towards robust regional estimates of CO2 sources and sinks using atmospheric transport models. Nature 415, 626630 (2002).

Acknowledgements

Ana Bastos was partially funded by the Portuguese Foundation for Science and Tech

nology (SFRH/BD/78068/2011). Clia M. Gouveia and Ricardo M. Trigo were supported

by QSECA (PTDC/AAGGLO/4155/2012). This research was also supported by the

European Research Council Synergy grant ERC-2013-SyG 610028 IMBALANCE-P.

Author contributions

A.B. conceived the study, conducted the analysis and wrote the paper; I.A.J., R.M.T. and

C.M.G. conceived the study and supervised the project; S.P., F.C. and C.R. developed

the data sets and provided technical and conceptual advice; C.R. also performed the

sensitivity tests for Jena inversion. P.C, J.P., P.F. and S.W.R. provided expert advice on

the different topics. All authors contributed to the revision of the paper.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/

http://www.nature.com/naturecommunications

Web End =naturecommunications Competing nancial interests: The authors declare no competing nancial interests. Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/

http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Bastos, A. et al. European land CO2 sink inuenced by NAO

and East-Atlantic Pattern coupling. Nat. Commun. 7:10315 doi: 10.1038/ncomms10315

(2016).

This work is licensed under a Creative Commons Attribution 4.0

International License. The images or other third party material in this

article are included in the articles Creative Commons license, unless indicated otherwise

in the credit line; if the material is not included under the Creative Commons license,

users will need to obtain permission from the license holder to reproduce the material.

To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 7:10315 | DOI: 10.1038/ncomms10315 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9