Biogeosciences, 11, 16831692, 2014www.biogeosciences.net/11/1683/2014/doi:10.5194/bg-11-1683-2014 Author(s) 2014. CC Attribution 3.0 License.

Biogeosciences

Open A

ccess

The seasonal cycle of 3CDIC in the North Atlantic subpolar gyre

V. Racap1, N. Metzl1, C. Pierre1, G. Reverdin1, P. D. Quay2, and S. R. Olafsdottir3

1Sorbonne Universits (UPMC, univ. Paris 6)-CNRS-IRD-MNHN, UMR7159, Laboratoire LOCEAN-IPSL,4, place Jussieu, 75005 Paris, France2School of Oceanography, University of Washington, Seattle, Washington, USA3Marine Research Institute, Skulagata 4, IS 121 Reykjavik, Iceland

Correspondence to: V. Racap (virginie.racape@locean-ipsl.upmc.fr)

Received: 12 August 2013 Published in Biogeosciences Discuss.: 29 August 2013Revised: 24 January 2014 Accepted: 11 February 2014 Published: 28 March 2014

Abstract. This study introduces for the first time the 13CDICseasonality in the North Atlantic subpolar gyre (NASPG)using 13CDIC data obtained in 20052006 and 20102012with dissolved inorganic carbon (DIC) and nutrient obser-vations. On the seasonal scale, the NASPG is characterizedby higher 13CDIC values during summer than during win-ter, with a seasonal amplitude between 0.70 0.10 (Au-gust 2010March 2011) and 0.77 0.07 (20052006).This is mainly attributed to photosynthetic activity in sum-mer and to a deep remineralization process during winterconvection, sometimes influenced by ocean dynamics andcarbonate pumps. There is also a strong and negative lin-ear relationship between 13CDIC and DIC during all sea-sons. Winter data also showed a large decrease in 13CDICassociated with an increase in DIC between 2006 and20112012, but the observed time rates (0.04 yr1and+1.7 mol kg1 yr1) are much larger than the expected an-thropogenic signal.

1 Introduction

The stable isotopic composition of dissolved inorganic car-bon (DIC) in the ocean is a useful semi-conservative tracerfor various oceanic analysis (e.g., water masses, paleo-circulation) and global carbon studies (constraint for theglobal carbon budget, interpreting anthropogenic carbon).This is reported relative to the reference V-PDB (ViennaPeeDee Belemnite) in terms of 13C in (Craig, 1957):

13C =

(

13C12C

)sample(

13C12C

)reference

1

1000. (1)The mean annual large-scale distribution of this tracer is

now relatively well documented at the sea surface and atdepth. Values range typically between 0.5 and 2.5 inthe modern surface ocean (Gruber et al., 1999). Regular ob-servations conducted since the 1970s1980s at a few timeseries stations or repetitive tracks were used to describe andinterpret its temporal variability. In subtropical oceanic re-gions, the seasonal 13C amplitude of 0.20.3 at the BATS(Bermuda Atlantic Time series) station and of 0.1 at theALOHA station (Pacific time series north of Hawaii) reflectthe local balance between airsea CO2 exchange, biologi-cal activity and vertical mixing (Gruber et al., 1998, 1999;Gruber and Keeling, 1999; Quay and Stutsman, 2003). Thisis also observed in the Southern Indian Ocean, with a sea-sonal 13C amplitude of 0.3 in the western part of thebasin (Racap et al., 2010), whereas it reaches 0.45 southof Tasmania (McNeil and Tilbrook, 2009). Studies have alsohighlighted a significant decrease in 13CDIC, attributed tothe anthropogenic CO2 uptake in the ocean. This results fromhuman-induced perturbations that reduce the 13CO2 in theatmosphere with a time rate of change of 0.024 yr1 (es-timate for the 19852008 period with the Alert station dataset; Keeling et al., 2010), as the combustion of the organiccarbon pool produces CO2 gas that is strongly depleted inthe heavy carbon isotope (13C). This signal, known as theoceanic 13C Suess effect, provides an additional constraintto estimate the anthropogenic carbon uptake in the ocean.

Published by Copernicus Publications on behalf of the European Geosciences Union.

1684 V. Racap et al.: 3CDIC seasonality in NASPG

This has been estimated as less than 0.007 yr1 in polarregions (McNeil et al., 2001; Olsen et al., 2006; Sonnerupand Quay, 2012) and as large as 0.025 yr1 in subtropi-cal regions (Bacastow et al., 1996; Gruber et al., 1999; Son-nerup and Quay, 2012). The high latitudes of the North At-lantic Ocean are considered to be one of the strongest anthro-pogenic CO2 sinks, a consequence of the large heat loss anddeep convection during winter as well as strong biologicalactivity during spring and summer. In this region, the tempo-ral variability (seasonal to decadal) of the 13CDIC still needsto be investigated. The wintersummer difference is expectedto be as large as 1 (Gruber et al., 1999), which could po-tentially mask the oceanic 13C Suess effect in surface water(McNeil et al., 2001). In this context, this study introducesfor the first time the 13CDIC seasonality in the North Atlanticsubpolar gyre using 13CDIC surface data obtained in 20052006 and 20102012 from merchant vessels Skogafoss andReykjafoss (14 cruises) along with DIC and nutrient measure-ments. Driving processes will be identified from isotopic dataand from the semi-quantitative tracer (113Cbio) of Gruberet al. (1999), which enables us to isolate the biological contri-bution to the seasonal amplitude from the physical contribu-tions (e.g., airsea gas exchange fractionation, ocean dynam-ics). The time rate of change in 13CDIC will be investigatedover the short period of 5 to 7 years and compared to otherregions.

2 Materials and methods

2.1 Data collection

Data used in this study were collected in the North Atlanticsubpolar gyre, mostly from merchant vessels Skogafoss andReykjafoss between Iceland and Newfoundland (Fig. 1) overtwo periods in 20052006 and 20102012. This is part ofsampling efforts carried from vessels of opportunity withinthe SURATLANT project initiated in 1993. The samplingperiods for each cruise and vessel are summarized in Table 1.

During all cruises, seawater was pumped from a depth of 35 m, depending on the ship. Sea surface salinity (SSS)was from discrete salinity samples. Sea surface temperature(SST) was either from the continuously recording ThermoS-alinoGraph (TSG seabird Electronics, Inc. USA, model 21)corrected for an average bias to sea surface temperature, orfrom nearly simultaneous XBT temperature profiles (at 3 m).SST and SSS accuracies were estimated at 0.2 C and 0.01,respectively. 13CDIC, DIC and nitrate concentrations (NO

3 )

were also obtained from discrete samples collected every4 h and analyzed later in land laboratories. Over the 20052006 period, DIC concentrations were estimated manometri-cally during the acid CO2 extraction procedure for 13CDICmeasurements from the helium stripping technique. This an-alytical method has been described previously by Quay etal. (1992, 2003). These measurements have an accuracy of

Fig. 1: Cruises track where 13CDIC was sampled in 2005-06 and 2010-12 (red dots). Dashed lines

indicate schematically surface main currents in the North Atlantic Ocean. The interior of the North

Atlantic SubPolar Gyre (NASPG, 53N-62N), our study area, is delimited by the black box.

Greenland Iceland

Labrador

Sea

Irminger

Sea

North Atlantic Ocean

Newfoundland

Fig. 1. Cruises track where 13CDIC was sampled in 20052006 and20102012 (red dots). Dashed lines indicate schematically surfacemain currents in the North Atlantic Ocean. The interior of the NorthAtlantic subpolar gyre (NASPG, 53 N62 N), our study area, isdelimited by the black box.

0.02 for 13CDIC based on a helium stripping tech-nique adapted from the one used by Kroopnick (1974) and 5 mol kg1 for DIC, based on a comparison to coulo-metric DIC values and to the Certified Reference Material(CRM) provided by Prof. A. Dickson (Scripps Institution ofOceanography, San Diego, USA). Over the 20102012 pe-riod during the SURATLANT cruises, 13CDIC were mea-sured via an acid CO2 extraction method in a vacuum sys-tem developed by Kroopnick (1974), whereas DIC were de-termined at the same time as total alkalinity (Alk) by a po-tentiometric titration derived from the method developed byEdmond (1970) using a closed cell and a CRM calibration.Further details on the sampling methods and analytical tech-niques are provided in Racap et al. (2010, 2013) for 13CDICand in Corbire et al. (2007) for DIC. Both parameters wereanalyzed at the LOCEAN-IPSL Laboratory (Paris, France).13CDIC values have a precision of 0.01 (Vangriesheimet al., 2009) and a reproducibility of 0.02 , whereas DICconcentrations have an accuracy of 3 mol kg1. Finally,nutrient concentrations were measured for all cruises withstandard colorometric methods at the Marine Research Insti-tute (Reykjavik, Iceland). The analytical procedure and thequality control for the nutrient analyses have been describedin detail in Olafsson et al. (2010), where the long-term ac-curacy has been estimated as 0.2 mol L1 for nitrate andsilicate, and 0.03 mol L1 for phosphate.

2.2 113Cbio and 13Cphys estimation

To determine the major process driving the seasonal 13CDICdistribution (t1 t2), we used two isotopic and semi-quantitative tracers, noted 113Cbio and 113Cphys, which

Biogeosciences, 11, 16831692, 2014 www.biogeosciences.net/11/1683/2014/

V. Racap et al.: 3CDIC seasonality in NASPG 1685

Table 1. Summary of sampling periods conducted on board the merchant vessels Skogafoss and Reykjafoss between Reykjavik (64 N,Iceland) and Newfoundland (47 N, Canada) over the 20052006 and 20102012 periods, respectively. The average values of 13CDIC, DICand salinity normalized DIC (nDIC) as well as their standard deviation (SD) were estimated for the interior of the North Atlantic subpolargyre (NASPG 53 N62 N). The number of data (n) used to calculate these mean values are indicated in the last column.

NASPG (53 N62 N)

Date Mean 13CDIC SD Mean DIC SD Mean nDIC SD n() (mol kg1) (mol kg1)

M/V Skogafoss

23 to 27 Jul 2005 1.65 0.10 2077.4 10.0 2093.6 5.3 1517 to 24 Sep 2005 1.46 0.11 2089.9 9.3 2107.9 4.3 812 to 16 Nov 2005 1.06 0.07 2120.7 10.6 2137.7 9.7 167 to 11 Feb 2006 0.87 0.05 2142.4 4.5 2144.1 7.3 2018 to 28 Apr 2006 0.90 0.09 2137.7 8.5 1225 to 27 Jun 2006 1.64 0.08 2058.7 8.4 2094.9 9.1 411 to 15 Nov 2006 0.99 0.12 2114.3 11.6 2118.6 10.3 13

M/V Reykjafoss

21 to 25 Aug 2010 1.36 0.11 2054.4 9.9 2074.0 9.2 2115 to 19 Dec 2010 0.97 0.10 2114.6 7.6 2131.2 5.3 1113 to 18 Mar 2011 0.66 0.08 2149.4 5.9 2150.8 6.3 1516 to 20 Dec 2011 0.69 0.05 2150.6 4.1 2151.0 6.9 2311 to 13 Mar 2012 0.61 0.04 2152.5 2.4 2152.1 9.5 1630 Jun to 4 Jul 2012 1.28 0.09 2089.7 9.1 2100.4 17.2 1121 to 25 Sep 2012 1.33 0.12 2095.7 13.5 2115.6 5.0 13

reveal the biological and physical signatures, respectively, asshown by Eqs. (2) and (3).

DICt2(13Ct2DIC) = DICt1(13Ct1DIC)

+

[1DIC(13C)

] t1t2bio

+

[1DIC(13C)

]t1t2phys

(2)

13Ct2DIC =

DICt1

DICt2(13C

t1DIC)

+ 113Ct1t2bio + 1

13Ct1t2phys . (3)

Described by Gruber et al. (1999), the biological pro-cesses contribution (113Cbio) results from the combinedphotosynthesis, respiration and remineralization processes(the soft-tissue pump) as well as the formationdissolutionof the calcium carbonate shells (the carbonate pump). Inthis study, the 113Cbio will be expressed according to ob-served biogeochemical conditions from one month or sea-son (t1) to another one (t2) by Eq. (4), with the stoichio-metric carbonnitrate ratio (rc:n) of 117 : 16 from Ander-son and Sarmiento (1994), the 13C value of the carbonate(13Ccarb) of 2 like Gruber et al. (1999), and the 13Cvalue of organic matter (13Corg) of 24 from Goerickeand Fry (1994). This last value depicts the average for thenorthern North Atlantic where sea surface temperature isaround 10 C in summer.

113Ct1t2bio = rc:n

13Corg

(NOt23 NO

t13

)DICt2

+ 13Ccarb

(Alkt2 Alkt1 + NOt23 NO

t13

)2DICt2

(4)

As the 13C variation associated with the carbonate pump(< 0.02 of the measurement error) does not contributemuch to the biological signature, this mechanism will be af-terwards neglected.

According to Eq. (3), the physical processes contribu-tion (airsea CO2 exchange fractionations and ocean dynam-ics; 113Cphys) corresponds to the difference between the13CDIC at t2 and the expected 13CDIC if biological activitywas the only process controlling the t1 to t2 distribution. Thisis expressed in Eq. (5).

113Ct1t2phys =

13Ct2DIC

(DICt1

DICt2(13Ct1DIC) + 1

13Ct1t2bio

). (5)

As t1 to t2 DIC change due to physical and biologicalprocesses (up to 100 mol kg1) is lower than the meanDIC concentration observed in the surface water of this re-gion (between 2050 and 2150 mol kg1), we assume thatDICt1/DICt2 equals 1.

3 Results and discussion

3.1 Surface 13CDIC distribution

During each transect, surface 13CDIC is distributed ac-cording to a northsouth gradient, with an amplitude be-tween 0.2 and 0.4 depending on the season and withhigher values close to the Newfoundland coast (up to 1.8 ;Fig. 2, bottom panel). This is opposite to the DIC distribution

www.biogeosciences.net/11/1683/2014/ Biogeosciences, 11, 16831692, 2014

1686 V. Racap et al.: 3CDIC seasonality in NASPG

Fig. 2. Surface distribution of salinity (top panel), dissolved inorganic carbon concentration (DIC, middle panel) and 13CDIC (bottompanel) between Reykjavik (Iceland, 64 N) and Newfoundland (Canada, 47 N) over the periods 20052006 (left panel) and 20102012(right panel). The interior of the North Atlantic subpolar gyre (NASPG, 53 N62 N), our study area, is delimited by the light grey box. Anysalinity measurements are realized in April 2006. August 2010 (full red square) distinguishes itself by a lower DIC concentration observedover both periods but without higher 13CDIC values.

(variation up to 100 mol kg1; Fig. 2, middle panel), whichincreases northward as observed previously by Corbire etal. (2007). These authors have also shown that the area southof 53 N was characterized by large spatial variability in SSS(Fig. 2, top panel), SST, DIC (Fig. 2, bottom panel) and nutri-ent concentrations when crossing the Labrador Current andapproaching the Newfoundland coast (Fig. 1). Influence ofcoastal waters was also identified north of 62 N near Icelandand on 13CDIC distribution (Fig. 2, bottom panel). Thus, weselected a region where both hydrological and carbon prop-erties (including 13CDIC) are much more homogeneous be-tween 53 N and 62 N, which corresponds to the interior ofthe North Atlantic subpolar gyre, denoted NASPG, west ofthe Reykjanes Ridge.

In the NASPG, the 13CDIC meridional gradient is weaker(close to 0.1 ). Highest values were observed during sum-mer (up to 1.65 0.10 in 20052006 and 1.36 0.11 in 20102012), while minima were measured during winter(down to 0.87 0.05 in 20052006 and 0.61 0.04 in

20102012, Fig. 2, bottom panel and Table 1). Comparedto 20052006, 13CDIC values collected in 20102012 arelower, suggesting a potential impact of anthropogenic carbonuptake. This would then require the correction of 13CDICobservations by this anthropogenic signal to compose theclimatological seasonal cycle. If this change in 13CDIC(decrease) represents the oceanic 13C Suess effect only, thisshould also be observed on DIC, with an expected theoreticalincrease of 57 mol kg1 in winter (Corbire et al., 2007;Metzl et al., 2010).

3.2 13CDIC versus DIC in NASPG

Figure 3a highlights a strong and negative link between13CDIC and DIC concentration in the NASPG over bothsampling periods. Lower 13CDIC in winter are associatedwith higher DIC concentrations, conversely for summer. Ac-cording to this figure, August 2010 is clearly outside the gen-eral trend observed here. This is a linear relationship fromwinter to winter, nevertheless with a difference in the values

Biogeosciences, 11, 16831692, 2014 www.biogeosciences.net/11/1683/2014/

V. Racap et al.: 3CDIC seasonality in NASPG 1687

Fig. 3: (a,b) 13CDIC versus dissolved inorganic carbon (DIC) concentration collected in the surface waters of the North Atlantic SubPolar

Gyre (NASPG, 53N-62N) in 2005-06 and 2010-12. (b) Dark grey squares represent 2005-06 data set, light grey squares symbolize the

2010-12 data set and colored symbols (like Fig. 3a legend) show mean values (and standard deviation) of 13CDIC and DIC concentration

for each sampling period, but without June 2006 because only 4 data were sampled in the south of the transect. The linear regressions

indicated with dark lines were estimated for the periods from July 2005 to June 2006 and from March 2011 to September 2012 (from the

bottom up: y -0.0108x+23.81, r=0.95 and y= -0.0108x+24.10 , r=0.94) and with dashed line was estimated for the period from August

2010 to March 2011 (y -0.007x+16.58, r=0.94).

1

3C

DIC

(

vs

V-P

DB

)

DIC concentration (mol kg -1)

a)

b)

1

3C

DIC

(

vs

V-P

DB

)

Fig. 3. (a, b) 13CDIC versus dissolved inorganic carbon (DIC) concentration collected in the surface waters of the North Atlantic subpolargyre (NASPG, 53 N62 N) in 20052006 and 20102012. (b) Dark grey squares represent the 20052006 data set, light grey squaressymbolize the 20102012 data set and colored symbols (like the Fig. 3a legend) show mean values (and standard deviation) of 13CDICand DIC concentration for each sampling period, but without June 2006 because only four data were sampled in the south of the transect.The linear regressions indicated with dark lines were estimated for the periods from July 2005 to June 2006 and from March 2011 toSeptember 2012 (from the bottom up: y 0.0108x + 23.81, r2 = 0.95 and y = 0.0108x + 24.10, r2 = 0.94), and with dashed lines wereestimated for the period from August 2010 to March 2011 (y 0.007x + 16.58, r2 = 0.94).

detected between all data collected in 20052006, except forNovember 2006 and those sampled in 20102012. If we fo-cus on the winter data set (FebruaryMarch) like Corbireet al. (2007) and Metzl et al. (2010), we notice an increasein DIC by 7 to 10 mol kg1 over the 5- to 6-year period.This is much larger than those expected from the anthro-pogenic effect ( 57 mol kg1). To remove natural vari-ation due to a potential salinity anomaly (Fig. 2, top panel),we evaluated the trend for salinity normalized DIC (nDICreported in Table 1) of 1.3 mol kg1 yr1 during this period(the 5- to 6-year period). This suggests a significant addi-tional contribution of the natural variability, thus masking theimpact of the anthropogenic carbon uptake in the NASPG re-gion. This comment is supported by the observed time ratesof change in 13CDIC (0.04 yr1) that are much largerthan those observed by Quay et al. (2007) in the NASPG(0.017 yr1 over the 19932007 period), Sonnerup andQuay (2012) in the Labrador Sea (0.014 yr1 over the19701995 period), Olsen et al. (2006) in the Nordic seas (upto 0.024 yr1 over the 19812003 period) and by Keel-ing et al. (2010) in the atmosphere (0.024 yr1 over the19852008 period). As the residence time of the surface wa-ter ( 1 yr) is shorter than the time needed for CO2 to reachisotopic equilibrium with the atmosphere ( 10 yr; Broeckerand Peng, 1974; Lynch-Stieglitz et al., 1995), this large de-crease in 13CDIC cannot be attributed to the anthropogenic

signal only. In addition, the mixed layer depth time series de-picted in Fig. 4b reveals two winter episodes, with the deep-est layers between 2005 and 2012, which should be asso-ciated with a higher DI12C input of subsurface water, thuscontributing to the 13CDIC decrease during this period (andDIC increase). As we suspect a large impact of the naturalvariability that we cannot estimate, observations are not cor-rected for an anthropogenic signal and the mean 13CDIC sea-sonal cycle will be composed for each detected period.

Study of the relationships between 13CDIC and DIC overboth sampling periods reveals, in fact, another important sig-nal to compose and interpret the 13CDIC seasonal cycle.According to Fig. 3b, observations are actually distributedwithin three periods: from July 2005 to June 2006, Au-gust 2010 to March 2011 and from March 2011 to Septem-ber 2012. These are characterized by the following linear re-gressions (Eqs. 6 to 8):

From July 2005 to June 2006,

13C20052006DIC = 0.0108 [DIC] + 24.10;

r2 = 0.94; slope error = 3.3 104. (6)

From August 2010 to March 2011,

13CAugust 2010March 2011DIC = 0.0074 [DIC] + 16.58;

r2 = 0.94; slope error = 3.1 104. (7)

www.biogeosciences.net/11/1683/2014/ Biogeosciences, 11, 16831692, 2014

1688 V. Racap et al.: 3CDIC seasonality in NASPG

Chlo

rophyll

a c

once

ntr

atio

n (

mg m

-3)

Mix

ed l

ayer

dep

th (

m)

Fig. 4: Distribution of a) Chlorophyll a concentration from MODIS-Aqua 4 km (disc.sci.gsfc.nasa.gov/giovanni, time series

generated between 44W-27W and 53N-62N) and of b) mixed layer depth (MLD, m) from Argo-float observations in

NASPG. The MLD is estimated on each profile as the depth at which temperature deviates by 0.1C from the temperature at

5m. The full line represents the median depth whereas the dashed lines correspond to the first and third quartiles

Jan.

2010

Jan.

2005

Jan.

2006

Jan.

2011

Jan.

2012

a)

b)

Fig. 4. Distribution of (a) chlorophyll a concentration from MODIS-Aqua 4 km (disc.sci.gsfc.nasa.gov/giovanni, time series generated be-tween 44 W27 W and 53 N62 N) and of (b) mixed layer depth (MLD, m) from Argo-float observations in NASPG. The MLD isestimated on each profile as the depth at which temperature deviates by 0.1 C from the temperature at 5 m. The full line represents themedian depth, whereas the dashed lines correspond to the first and third quartiles.

From March 2011 to September 2012

13C20112012DIC = 0.0108 [DIC] + 23.81;

r2 = 0.95; slope error = 2.9 104. (8)

Based on a Students t test and its associated error,there is no significant difference between the slopes of re-gression 1 (Eq. 6) and 3 (Eq. 8). These are also equiva-lent to those expected if the soft-tissue pump (photosynthe-sis/respiration/remineralization) were the only mechanismaffecting DIC and 13CDIC seasonal distributions (0.011 (mol kg1)1 with 13Corg of 24 and rc:n of 117 : 16).As a result, the slope of regression 2 (Eq. 7), estimated fromAugust 2010 to March 2011, is significantly different fromthe other two (0.007 (mol kg1)1). This is influencedby the August 2010 DIC values that are the lowest concen-trations reported in this study (Table 1, as well as for nDIC)but are not associated with the highest 13CDIC as expected.Indeed, high productivity recorded in the whole region dur-ing the summer of 2010 (Fig. 4a) would explain the lowDIC values (2054.40 9.92 mol kg1, Table 1) but, accord-ing to previous results, an enhanced photosynthetic activitywould increase 13CDIC linearly with a ratio of 0.011 (mol kg1)1 to reach a maximum value. This suggests thatadditional processes also contribute to DIC and 13CDIC sea-sonal distributions during this period, thus changing their re-lationship. The estimate of the isotopic tracers (113Cbio and

113Cphys) might be clarified to interpret this peculiar situa-tion.

3.3 13CDIC seasonal cycle and driving processes inNASPG

3.3.1 The period July 2005June 2006

Based on satellite information (the monthly SST compos-ite of MODIS Aqua 9 km, disc.sci.gsfc.nasa.gov/giovanni),the SURATLANT data set (NO3 and DIC) and Eq. (6), webuilt the climatological seasonal cycle of SST, NO3 , DICand 13CDIC for 20052006 (Fig. 5). Because DIC concen-trations have slightly changed during this period (by less thanthe climatological cycle standard deviation; Corbire et al.,2007; Ullman et al., 2009), we chose to consider the avail-able DIC (and NO3 ) data set sampled between 2001 and2006, which is much more complete for composing the sea-sonality. The observations collected on board M/V Skogafossduring this period fit well within this climatology (Fig. 5band c).

Figure 5 thus reveals the mean seasonal cycle for 13CDICin 20052006, with an amplitude of 0.77 0.07 , inphase with SST (6 C) and in opposite phase with NO3(9 mol kg1) and DIC (70 mol kg1). Minima are ob-served in FebruaryMarch, whereas maxima are estimated

Biogeosciences, 11, 16831692, 2014 www.biogeosciences.net/11/1683/2014/

disc.sci.gsfc.nasa.gov/giovannidisc.sci.gsfc.nasa.gov/giovanni

V. Racap et al.: 3CDIC seasonality in NASPG 1689

b)

a)

Fig. 5: Mean seasonal cycle of a) Sea

Surface Temperature (SST), b) nitrate

(NO3-), c) Dissolved Inorganic Carbon

(DIC), and d) 13CDIC in the surface

waters of the NASPG (53N-62N).

Light grey squares (and standard

deviation) symbolize the 2010-12 data

set whereas dark grey squares (and

standard deviation) represent the 2005-

06 data set. Black dashed curves and

their standard deviation were

established from a) MODIS-aqua data

set recorded in 2005-06 (monthly

composite with a spatial resolution of

9km, generated by NASAs Giovanni,

disc.sci.gsfc.nasa.gov/giovanni), b and

c) both NO3- and DIC measurements

collected during SURATLANT cruises

between 2002 and 2006 and c) from

2005-06 13CDIC DIC relationship

presented in Fig. 3b. Anomalies are

identified by colored symbols

following the legend of Fig. 3a.

Fig. 5. Mean seasonal cycle of (a) sea surface temperature (SST),(b) nitrate (NO3 ), (c) dissolved inorganic carbon (DIC), and (d)

13CDIC in the surface waters of the NASPG (53 N62 N). Light

grey squares (and standard deviation) symbolize the 20102012data set, whereas dark grey squares (and standard deviation) rep-resent the 20052006 data set. Black dashed curves and their stan-dard deviation were established from (a) the MODIS-aqua data setrecorded in 20052006 (monthly composite with a spatial resolu-tion of 9 km, generated by NASAs Giovanni, disc.sci.gsfc.nasa.gov/giovanni), b and c) both NO3 and DIC measurements col-lected during SURATLANT cruises between 2002 and 2006 and(c) from the 20052006 13CDICDIC relationship presented inFig. 3b. Anomalies are identified by colored symbols following thelegend of Fig. 3a.

in August during the summer bloom (Antoine et al., 2005;Corbire et al., 2007). Based on Eqs. (4) and (5), processescontrolling this seasonality could be identified. Table 2 sum-marizes the results. The standard deviation displays the

meridional gradient reported previously. Positive 113Cbiois associated with phytoplankton activity, which preferen-tially uses the light carbon isotope during the photosyn-thetic mechanism (1NOWS3 < 0 mol kg

1, nutrient con-sumption), whereas negative 113Cbio highlights the rem-ineralization process, which releases much more DI12C(1NOWS3 > 0, nutrient input). Based on 20052006 values(Fig. 5), we estimate a 113Cbio of 0.74 0.12 for thewinter to summer amplitude of 0.77 0.07 . During thisperiod, the 13CDIC seasonality observed in the NASPG sur-face waters is largely controlled by biological activity, whichincreased 13CDIC values during summer (photosynthesis)and decreased them during winter (remineralization). Thesehigh concentrations of NO3 and DIC related to remineral-ization are attributed to the winter entrainment of subsurfacewater (Takahashi et al., 1993).

According to the previous section, data from Novem-ber 2006 do not get into line with other data collectedin 20052006. Results obtained from isotopic tracers indi-cate that physical processes also contribute to controllingthe 13CDIC and DIC distributions from summer to Novem-ber 2006 (113Cphys = 0.23 0.15 ; Table 2). As noanomaly was detected in November 2006 (e.g., SST, MLD),this is difficult to assess further.

3.3.2 The period August 2010March 2011

There were only three cruises conducted between Au-gust 2010 and March 2011 (Table 1), but 13CDIC were mea-sured during winter and summer, when maximal and mini-mal values are expected. According to this comment, we can-not compose the mean 13CDIC seasonal cycle for this periodbut, we can evaluate a seasonal amplitude of 0.70 0.10 for 13CDIC and of 92.9 5.9 mol kg1 for DIC concen-tration. The biological and physical 13CDIC signature es-timated from summer 2010 to winter 2011 is much higher(0.97 0.06 and 0.27 0.13 , respectively; Table 2)than those obtained from the 20052006 data set. This is inagreement with the hypothesis of high productivity that ismasked by additional processes.

The extensive winter mixed layer observed in 2009(Fig. 4b) could have lowered 13CDIC at the sea surfaceand persisted through to August 2010. The lack of 13CDICmeasurements in the winter of 2010 precludes estimating113Cphys directly. If photosynthesis were the only pro-cess controlling the seasonal amplitude, as was observedduring the 20052006 cycle, the 13CDIC values of thewinter of 2010 would be around 0.36 (13CAugust 2010DIC +

113CWAugust 2010bio = 1.361.00). However, such low val-ues have not yet been observed in the NASPG (Quay et al.,2007; Racap et al., 2013) and are unlikely.

The NASPG also experienced a significant anomalouswarming (+2 C) during 2010. According to Zhang etal. (1995), a warming of 2 C would decrease 13CDIC by

www.biogeosciences.net/11/1683/2014/ Biogeosciences, 11, 16831692, 2014

disc.sci.gsfc.nasa.gov/giovannidisc.sci.gsfc.nasa.gov/giovanni

1690 V. Racap et al.: 3CDIC seasonality in NASPG

Table 2. Seasonal (t1 to t2) amplitude of 13CDIC (113CDIC) and their standard deviation (SD). The biological and physical contributions(113Cbio and 1

13Cphys, respectively) for the 20052006 and 20102012 periods estimated from Eqs. (4) and (5).

t1 to t2 113CDIC SD () 1

13Cbio SD () 113Cphys SD ()

20052006 period

Mean Feb to mean Aug 0.77 0.07 0.74 0.12 0.03 0.14of black dashed curves in Fig. 5Mean Aug to Nov 2006 0.65 0.14 0.42 0.10 0.23 0.15

20102012 period

Mar 2010 to Aug 2010 NA 1.00 0.10 NAAug 2010 to Dec 2010 0.44 0.09 0.61 0.11 0.18 0.11Aug 2010 to Mar 2011 0.70 0.10 0.97 0.06 0.27 0.13Mar 2012 to Jun 2012 0.68 0.14 0.76 0.11 0.08 0.12

0.2 at the isotopic equilibrium conditions, thus maskingthe biological signal on 13CDIC. In the subpolar region, thesurface 13CDIC is never in isotopic equilibrium with the at-mospheric 13CO2 (Quay et al., 2003). Nevertheless, if weassume that change in SST affects 13CDIC substantially atthe seasonal scale in this region, this is in opposite effect tothe uptake of 13C poor CO2 in surface water. Data collectedin December 2010 are relatively higher in 13CDIC and lowerin DIC than those sampled in December 2011. This suggestseither that the impact due to the change in SST in 2010 isdominated by the CO2 uptake effect, or that additional phys-ical driving processes, like horizontal transport, influencedthe 13CDIC and DIC seasonal distributions in 2010.

In the subpolar region, phytoplanktonic bloom is regularlydominated by the coccolithophore species, which are pro-moted by specific environmental conditions like high SSTand shallow MLD (Fig. 4b; Raitsos et al., 2006) observedduring this period. Coccolithophores are photosynthetic cal-cifying algal species. As we mentioned in Sect. 3.2, a small13C fractionation occurs during the formation of the calciumcarbonate coccoliths so that it is insignificant for the 13CDICvariation. Coccolithophores bloom starts generally in Juneand would induce by their soft-tissue pump a winter to sum-mer increase in 13CDIC and decrease in DIC according to aratio of 0.011 mol kg1 and by their carbonate pump,an additional decrease in DIC without change in 13CDIC(conversely from summer to winter), thus impacting the sea-sonal relationships observed between these both parametersover the period August 2010March 2011. There is no in-formation on 13CDIC in the spring of 2010, and the cloudyconditions during summer precluded observations with satel-lite imagery, so that it is difficult to confirm this hypothesis.

According to regression slope 2 (Eq. 7), the lowest DICvalues associated with low 13CDIC and the isotopic tracers(113Cbio and 113Cphys), DIC and 13CDIC seasonal distri-butions may be controlled by these combined processes.

3.3.3 The period March 2011September 2012

In 20112012, the seasonal minimum was sampled two timesbut the seasonal maximum was never reached (Table 1).As all DIC data collected in the NASPG region betweenMarch 2011 and September 2012 are presented in this study,it is impossible to complete the seasonal DIC cycle like in20052006 and to deduce that of 13CDIC (Fig. 5). Basedon 2012 values (from March 2012 to June 2012; Fig. 5),we can nevertheless observe a seasonal 13CDIC amplitudeof at least 0.68 0.14 , with a biological contribution of0.76 0.11 (Table 2). These results are comparable tothose observed in 20052006, suggesting that the soft-tissuepump was the major driving process of 13CDIC seasonalityin the NASPG region.

4 Conclusions

In this study, we present new sea surface observations of13CDIC and associated hydrological and biogeochemicalproperties (DIC, NO3 ) in the NASPG (53

N62 N) basedon 14 transects conducted in 20052006 and 20102012.During all transects, we observe a positive northsouth gra-dient coherent with model studies of Gruber et al. (1999)and Tagliabue and Bopp (2008), and opposite to DIC gra-dients (high DIC in the north, low in the south). We alsodetect a very large range of 13CDIC depending on the pe-riods. Coupling 13CDICDIC observations enables to re-veal strong parallel relationships in 20052006 and in 20112012. For both tracers and for winter (FebruaryMarch) thechange from 20052006 to 20112012 is coherent with an-thropogenic carbon uptake (increase in DIC, decrease in13CDIC), but with a time rate of change much larger than thetheoretical value due to anthropogenic effect. Deeper wintermixing recorded in the winter of 20082009 and in the win-ter of 20112012 most likely contributed to reduce 13CDICand increase DIC, thus masking the anthropogenic signal.

Biogeosciences, 11, 16831692, 2014 www.biogeosciences.net/11/1683/2014/

V. Racap et al.: 3CDIC seasonality in NASPG 1691

During these both periods, the seasonality of 13CDIC wasinvestigated taking into account the observed changes andwas totally computed for the years 20052006. Our cli-matological seasonal cycle shows a large 13CDIC varia-tion of 0.77 0.07 (larger in summer, lower in winter),which represents the largest seasonal amplitude recorded atthe sea surface in this region by Gruber et al. (1999) andbuilt by Jonkers et al. (2013). This is coherent with large DIC(70 mol kg1), nutrients and pCO2 seasonal signals previ-ously observed in the region (Corbire et al., 2007). A diag-nostic process model suggests that the biological processescontrol the 13CDIC seasonality, thus masking airsea CO2exchange and the ocean dynamics effect. The comparison be-tween the slope obtained from the 13CDICDIC relationshipand those expected if the soft-tissue pump is the only drivingprocess emphasized the importance of the biological pumpin DIC seasonality in the subpolar North Atlantic region.

The 13CDICDIC relationship has also helped to identifya third period (from August 2010 to March 2011) charac-terized by a lower regression slope, significantly differentfrom the other two. This may be due to a combination ofprocesses such as the formation of a carbonate shell asso-ciated with a coccolitophore bloom, the temperature effecton the carbon isotopic fractionation, and the ocean dynam-ics processes. Coccolithophore bloom being common in theNASPG surface water, it will be interesting to estimate therelationship between 13CDIC and DIC in such a bloom inorder to validate one of these hypotheses and to better under-stand the carbon cycle evolution in this region.

The observed 13CDIC seasonal cycle now being well es-tablished, a challenge for future analysis will be to investigatethis seasonality in ocean models. We strongly recommendmaintaining such observations in the future to better separatenatural and anthropogenic signals. The 13CDIC is also ableto estimate the oceanic 13C in disequilibrium with the atmo-sphere (13Cdes) used in the atmospheric inversion models toevaluate CO2 flux. This study has revealed a large seasonalvariability of 13CDIC that strongly influenced the 13Cdesseasonality, although we have observed a winter to summerchange in SST of 6 C. Taking into account this variabilitywould improve the CO2 flux uncertainty due to the evaluationof the 13Cdes. This highlights the need to continue measur-ing both parameters to investigate the anthropogenic carbonsignals, as well as to better constrain atmospheric inversionsand validate ocean carbon models.

Acknowledgements. This study has been supported by EU FP7project CARBOCHANGE (264879) and by CNRS/INSU in France(O.R.E. SSS and OCEANS-C13 projects). We thank the EIMSKIPCompany, the captains and crews of the M/V Skogafoss and theM/V Reykjafoss, all the very helpful embarked observers for the13C and DIC sampling, A. Naamar and the SNAPO-CO2 for theirhelp in laboratory analysis and C. Deschamps-Berger for the MLDtime series. We thank also L. Bopp for fruitful discussions and twoanonymous reviewers for their constructive comments.

Edited by: K. Suzuki

The publication of this article isfinanced by CNRS-INSU.

References

Anderson, L. A. and Sarmiento, J. L.: Redfield ratios of rem-ineralization determined by nutrient data analysis, Global Bio-geochem. Cy., 8, 6580, 1994.

Antoine, D., Morel, A., Gordon, H. R., Banzon, V. F., and Evans,R. H.: Bridging ocean color observations of the 1980s and 2000sin search of long term trends, J. Geophys. Res., 110, C06009,doi:10.1029/2004JC002620, 2005.

Bacastow, R. B., Keeling, C. D., Lueken, T. J., Wahlen, M., andMook, W. G.: The 13C Suess effect in the world surface oceansand its implications for oceanic uptake of CO2: analysis of ob-servations at Bermuda, Global Biogeochem. Cy., 10, 335346,1996.

Broecker, W. S. and Peng, T. H.: Gas exchange rates between airand sea, Tellus 26, 2135, 1974.

Corbire, A., Metzl, N., Reverdin, G., Brunet, C., and Takahashi,T.: Interannual and decadal variability of the oceanic carbon sinkin the North Atlantic SubPolar Gyre, Tellus B, 59, 168178,doi:10.1111/j.1600-0889.2006.00232.x, 2007.

Craig, H.: Isotopic standards for carbon and oxygen and correc-tion factor for mass-spectrometric analysis of carbon dioxide,Geochim. Cosmochim. Acta, 12, 133149, 1957.

Edmond, J. M.: High precision determination of titration alkalinityand total carbon dioxide content of sea water by potentiometrictitration, Deep-Sea Res., 17, 737750, 1970.

Goericke, R. and Fry, B.: Variations of marine plankton 13C withlatitude, temperature, and dissolved CO2 in the world ocean,Global Biogeochem. Cy., 8, 8590, 1994.

Gruber, N. and Keeling, C. D.: Seasonal carbon cycling in the Sar-gasso Sea near Bermuda, Scripps Inst. Oceanogr. Bull., 30, 96pp., 1999.

Gruber, N., Keeling, C. D., and Stocker, T. F.: Carbon-13 constraintson the seasonal inorganic carbon budget at the BATS site in theNorthwestern Sargasso Sea, Deep-Sea Res. Pt. I, 45, 673717,1998.

Gruber, N., Keeling, C. D., Bacastow, R. B., Guenther, P. R., Lueker,T. J., Whalen, M., Meijer, H. A. J., Mook, W. G., and Stocker,T. F.: Spatiotemporal patterns of carbon-13 in the global surface

www.biogeosciences.net/11/1683/2014/ Biogeosciences, 11, 16831692, 2014

http://dx.doi.org/10.1029/2004JC002620http://dx.doi.org/10.1111/j.1600-0889.2006.00232.x

1692 V. Racap et al.: 3CDIC seasonality in NASPG

oceans and the oceanic Suess effect, Global Biogeochem. Cy.,13, 307335, 1999.

Jonkers, L., van Heuven, S., Zahn, R., and Peeters, F. J. C.: Sea-sonal patterns of shell flux, 18O and 13C of small and large N.pachyderma(s) and G. bulloides in the subpolar North Atlantic,Paleoceanography, 28, 111, doi:10.1002/palo.20018, 2013.

Keeling, R. F., Piper, S. C., Bollenbacher, A. F., and Walker, S. J.:Monthly atmospheric 13C/12C isotopic ratios for 11 SIO stations,in: Trends: A Compendium of Data on Global Change, CarbonDioxide Information Analysis Center, Oak Ridge National Labo-ratory, US Department of Energy, Oak Ridge, Tenn., USA, 2010.

Kroopnick, P.: Correlations between 13C and 6CO2 in surface wa-ters and atmospheric CO2, Earth Planet. Sci. Lett., 22, 397403,1974.

Lynch-Stieglitz, J., Stocker, T. F., Broecker, W. S., and Fairbanks,R. G.: The influence of air-sea exchange on the isotopic com-position of oceanic carbon: Observations and modelling, GlobalBiogeochem. Cy., 9, 653665, 1995.

McNeil, B. and Tilbrook, B.: A seasonal carbon budget for the sub-antarctic Ocean, South of Australia, Mar. Chem., 115, 196210,2009.

McNeil, B., Matear, R., and Tilbrook, B.: Does carbon 13 track an-thropogenic CO2 in Southern ocean, Global Biogeochem. Cy.,15, 597613, 2001.

Metzl, N., Corbire, A., Reverdin, G., Lenton, A., Takahashi, T.,Olsen, A., Johannessen, T., Pierrot, D., Wanninkhof, R., lafs-dttir, S. R., Olafsson, J., and Ramonet, M.: Recent accelerationof the sea surface fCO2 growth rate in the North Atlantic Sub-Polar Gyre (19932008) revealed by winter observations, GlobalBiogeochem. Cy., 24, 113, doi:10.1029/2009GB003658, 2010.

Olafsson, J., Olafsdottir, S. R., Benoit-Cattin, A., and Takahashi,T.: The Irminger Sea and the Iceland Sea time series measure-ments of sea water carbon and nutrient chemistry 19832008,Earth Syst. Sci. Data, 2, 99104, doi:10.5194/essd-2-99-2010,2010.

Olsen, A., Omar, A. M., Bellerby, R. G. J., Johannessen, T., Ninne-mann, U., Brown, K. R., Olsson, K. A., Olafsson, J., Nondal, G.,Kivime, C., Kringstad, S., Neill, C., and lafsdttir, S. R.: Mag-nitude and origin of the anthropogenic CO2 increase and Suesseffect in the Nordic seas since 1981, Global Biogeochem. Cy.,20, 112, doi:10.1029/2005GB002669, 2006.

Quay, P. D. and Stutsman, J.: Surface layer carbon budget forthe subtropical N. Pacific: 13C constraints at station ALOHA,Deep-Sea Res. Pt. I, 50, 10451061, 2003.

Quay, P. D., Tilbrook, B., and Wong, C.: Oceanic uptake of fossilfuel CO2: Carbon-13 evidence, Science, 256, 7479, 1992.

Quay, P. D., Sonnerup, R., Westby, T., Stutsman, J., and McNichol,A.: Changes in the 13C/12C of dissolved inorganic carbon in theocean as a tracer of anthropogenic CO2 uptake, Global Biogeoch.Cy., 17, 1004, doi:10.1029/2001GB001817, 2003.

Quay, P. D., Sonnerup, R., Stutsman, J., Maurer, J., Krtzinger, A.,Padin, X. A., and Robinson, C.: Anthropogenic CO2 accumula-tion rates in the North Atlantic Ocean from changes in the 13C/12C of dissolved inorganic carbon, Global Biogeoch. Cy., 21, 115, doi:10.1029/2006GB002761, 2007.

Racap, V., Lo Monaco, C., Metzl, N., and Pierre, C.: Sum-mer and winter distribution of 13CDIC in surface waters ofthe South Indian Ocean (20 S60 S), Tellus B, 62, 660673,doi:10.1111/j.1600-0889.2010.00504.x, 2010.

Racap, V., Pierre, C., Metzl, N., Lo Monaco, C., Reverdin, G.,Olsen, A., Morin, P., Ros, A. F., Vzquez-Rodriguez, M., andPrez, F. F.: Anthropogenic carbon changes in the Irminger Basin(19812006): Coupling 13CDIC and DIC observations, J. Ma-rine Syst., 126, 2432, doi:10.1016/j.jmarsys.2012.12.005, 2013.

Raitsos, D. E., Lavender, S. J., Pradhan, Y., Tyrell, T., Reid, P. C.,and Edwards, M.: Coccolithophore bloom size variation in re-sponse to the regional environment of the subarctic North At-lantic, Limnol. Oceanogr., 51, 21222130, 2006.

Sonnerup, R. E. and Quay, P. D.: 13C constraints on oceancarbon cycle models, Global Biogeochem. Cy., 26, GB2014,doi:10.1029/2010GB003980, 2012.

Tagliabue, A. and Bopp, L.: Towards understanding global variabil-ity in ocean carbon-13, Global Biogeochem. Cy., 22, GB1025,doi:10.1029/2007GB003037, 2008.

Takahashi, T., Olafsson, J., Goddard, J. G., Chipman, D. W., andSutherland, S. C.: Seasonal variations of CO2 and nutrients in thehigh latitude surface oceans: a comparative study, Global Bio-geochem. Cy., 7, 843878, 1993.

Ullman, D. J., Galen, A., McKinley, G. A., Bennington, V.,and Dutkiewicz, S.: Trends in the North Atlantic carbonsink: 19922006, Global Biogeochem. Cy., 23, GB4011,doi:10.1029/2008GB003383, 2009.

Vangriesheim, A., Pierre, C., Aminot, A., Metzl, N., Baurand, F.,and Caprais, J.-C.: The influence of Congo river discharges inthe surface and deep layers of the Gulf of Guinea, Deep-Sea Res.Pt. II, 56, 21832196, doi:10.1016/j.dsr2.2009.04.002, 2009.

Zhang, J., Quay, P. D., and Wilbur, D. O.: Carbon isotope frac-tionation during gas-water exchange and dissolution of CO2,Geochim. Cosmochim. Ac., 59, 107114, 1995.

Biogeosciences, 11, 16831692, 2014 www.biogeosciences.net/11/1683/2014/

http://dx.doi.org/10.1002/palo.20018http://dx.doi.org/10.1029/2009GB003658http://dx.doi.org/10.5194/essd-2-99-2010http://dx.doi.org/10.1029/2005GB002669http://dx.doi.org/10.1029/2001GB001817http://dx.doi.org/10.1029/2006GB002761http://dx.doi.org/10.1111/j.1600-0889.2010.00504.xhttp://dx.doi.org/10.1016/j.jmarsys.2012.12.005http://dx.doi.org/10.1029/2010GB003980http://dx.doi.org/10.1029/2007GB003037http://dx.doi.org/10.1029/2008GB003383http://dx.doi.org/10.1016/j.dsr2.2009.04.002