Nolwenn Lemaitre
email: nolwenn.lemaitre@erdw.ethz.ch

Hélène Planquette

Frédéric Planchon

Géraldine Sarthou

Stéphanie Jacquet

Maribel I García-Ibáñez

Arthur Gourain

Marie Cheize

Laurence Monin

Luc André

Priya Laha

Herman Terryn

Frank Dehairs

Particulate barium tracing of significant mesopelagic carbon remineralisation in the North Atlantic

Abstract. The remineralisation of sinking particles by prokaryotic heterotrophic activity is important for controlling oceanic carbon sequestration. Here, we report mesopelagic particulate organic carbon (POC) remineralisation fluxes in the North Atlantic along the GEOTRACES-GA01 section (GEOVIDE cruise; May-June 2014) using the particulate biogenic barium (excess barium; Ba xs ) proxy. Important mesopelagic (100-1000 m) Ba xs differences were observed along the transect depending on the intensity of past blooms, the phytoplankton community structure, and the physical forcing, including downwelling. The subpolar province was characterized by the highest mesopelagic Ba xs content (up to 727 pmol L -1 ), which was attributed to an intense bloom averaging 6 mg chl a m -3 between January and June 2014 and by an intense 1500 m deep convection in the central Labrador Sea during the winter preceding the sampling. This downwelling could have promoted a deepening of the prokaryotic heterotrophic activity, increasing the Ba xs content. In comparison, the temperate province, characterized by the lowest Ba xs content (391 pmol L -1 ), was sampled during the bloom period and phytoplankton appear to be dominated by small and calcifying species, such as coccolithophorids. The Ba xs content, related to oxygen consumption, was converted into a remineralisation flux using an updated relationship, proposed for the first time in the North Atlantic. The estimated fluxes were of the same order of magnitude as other fluxes obtained using independent methods (moored sediment traps, incubations) in the North Atlantic. Interestingly, in the subpolar and subtropical provinces, mesopelagic POC remineralisation fluxes (up to 13 and 4.6 mmol C m -2 d -1 , respectively) were equalling and occasionally even exceeding upper-ocean POC export fluxes, deduced using the 234 Th method. These results highlight the important impact of the mesopelagic remineralisation on the biological carbon pump of the studied area with a near-zero, deep (> 1000 m) carbon sequestration efficiency in spring 2014.

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

Introduction

The ocean represents the largest active CO 2 sink [START_REF] Sabine | The Oceanic Sink for Anthropogenic CO 2[END_REF] partly materialised by the oceanic biological carbon pump (BCP), which controls the export of carbon and nutrients to the deep ocean through the production of biogenic sinking particles [START_REF] Boyd | Understanding the export of biogenic particles in oceanic waters: Is there consensus?[END_REF][START_REF] Sigman | Glacial/interglacial variations in atmospheric carbon dioxide[END_REF][START_REF] Volk | Ocean carbon pumps: analysis of relative strengths and efficiencies in ocean-driven atmospheric CO 2 changes, in: The carbon cycle and atmospheric CO 2 : Natural variations Archean to Present[END_REF]. The North Atlantic sustains one of the most productive spring phytoplankton blooms of the world ocean [START_REF] Esaias | Monthly satellite-derived phytoplankton pigment distribution for the North Atlantic basin[END_REF][START_REF] Henson | Decadal variability in North Atlantic phytoplankton blooms[END_REF][START_REF] Longhurst | Ecological geography of the sea[END_REF]Pommier et al., 2009). The high primary productivity in combination with the water mass formation there as part of the thermohaline circulation [START_REF] Seager | Is the Gulf Stream responsible for Europe's mild winters?[END_REF] results in a particularly efficient BCP in the North Atlantic [START_REF] Buesseler | Carbon and nitrogen export during the JGOFS North Atlantic Bloom experiment estimated from 234 Th : 238 U disequilibria[END_REF][START_REF] Buesseler | Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean[END_REF]Herndl and Reinthaler, 2013;[START_REF] Honjo | Annual biogenic particle fluxes to the interior of the North Atlantic Ocean; studied at 34 • N 21 • W and 48 • N 21[END_REF]Le Moigne et al., 2013b), estimated to contribute up to 18 % of the global oceanic BCP [START_REF] Sanders | The Biological Carbon Pump in the North Atlantic[END_REF]. However, the magnitude of the carbon transfer to the deep ocean depends on many factors including the efficiency of bacterial remineralisation within the mesopelagic layer (100-1000 m depth layer). In this layer, most of the particulate organic carbon (POC) exported from the upper mixed layer is respired or released to the dissolved phase as dissolved organic carbon (DOC; [START_REF] Buesseler | Revisiting Carbon Flux Through the Ocean's Twilight Zone[END_REF][START_REF] Buesseler | Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean[END_REF][START_REF] Burd | Towards a transformative understanding of the ocean' s biological pump: Priorities for future research[END_REF]Herndl and Reinthaler, 2013;[START_REF] Lampitt | Particle flux in deep seas: Regional characteristics and temporal variability[END_REF][START_REF] Martin | VERTEX: carbon cycling in the northeast Pacific[END_REF]. Mesopelagic remineralisation has often been reported to balance or even exceed the carbon supply from the surface (i.e. POC and DOC; [START_REF] Aristegui | Microbial oceanography of the dark ocean' s pelagic realm[END_REF][START_REF] Baltar | Evidence of prokaryotic metabolism on suspended particulate organic matter in the dark waters of the subtropical North Atlantic[END_REF][START_REF] Burd | Assessing the apparent imbalance between geochemical and biochemical indicators of meso-and bathypelagic biological activity: What the @$#! is wrong with present calculations of carbon budgets?[END_REF][START_REF] Collins | The multiple fates of sinking particles in the North Atlantic Ocean[END_REF][START_REF] Fernández-Castro | Mesopelagic respiration near the ESTOC (European Station for Time-Series in the Ocean, 15.5 • W, 29.1 • N) site inferred from a tracer conservation model[END_REF]Giering et al., 2014;[START_REF] Reinthaler | Prokaryotic respiration and production in the meso-and bathypelagic realm of the eastern and western North Atlantic basin[END_REF], highlighting the impact of mesopelagic processes on bathypelagic carbon sequestration. Unfortunately, studies focusing on the mesopelagic layer are scarce, and the remineralisation process in this part of the water column remains poorly constrained. A variety of methods have been used to assess deep remineralisation. The attenuation of the particulate organic matter concentration with depth can be deduced from POC fluxes recorded by bottom-tethered or free-floating neutrally buoyant sediment traps (e.g. [START_REF] Buesseler | Revisiting Carbon Flux Through the Ocean's Twilight Zone[END_REF][START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF][START_REF] Martin | VERTEX: carbon cycling in the northeast Pacific[END_REF] deployed at different depths. Bacterial respiration can be assessed by measuring the rate of dissolved oxygen consumption, but this approach is usually limited to the upper 200 m of depth because of sensitivity issues [START_REF] Arístegui | Respiration in the mesopelagic and bathypelagic zones of the oceans[END_REF]Christaki et al., 2014;[START_REF] Lefèvre | The metabolic balance at two contrasting sites in the Southern Ocean: The ironfertilized Kerguelen area and HNLC waters[END_REF]. However, sediment traps and direct respiration measurements are insufficiently reliable for depths exceeding 200 to 500 m (i.e. the lower mesopelagic area). Earlier work has revealed that the accumulation of particulate biogenic barium (excess barium; Ba xs ) in the mesopelagic water column (100-1000 m) is related with organic carbon remineralisation. This biogenic Ba is essentially carried by barite (BaSO 4 ) microcrystals, which form inside oversaturated micro-environments, mostly aggregates of organic material where prokaryotic activity is intense [START_REF] Bishop | The barite-opal-organic carbon association in oceanic particulate matter[END_REF][START_REF] Collier | The trace element geochemistry of marine biogenic particulate matter[END_REF][START_REF] Dehairs | Discrete suspended particles of barite and the barium cycle in the open ocean[END_REF][START_REF] Ganeshram | An experimental investigation of barite formation in seawater[END_REF][START_REF] Gonzalez-Munoz | Precipitation of barite by Myxococcus xanthus: possible implications for the biogeochemical cycle of barium[END_REF]. Bacterial activity will result in the disruption of these aggregates, thereby releasing the barite crystals in the ambient water. As a result, the concentration of Ba xs relates with oxygen consumption rate [START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF][START_REF] Shopova | A simple model of biogeochemical element distribution in the oceanic water column[END_REF] and can be converted into a remineralisation rate of POC in the mesopelagic layer [START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF]. Ba xs has been successfully used as a proxy of POC remineralisation flux in the Southern Ocean [START_REF] Cardinal | Variations of carbon remineralisation in the Southern Ocean illustrated by the Baxs proxy[END_REF]Jacquet et al., 2008aJacquet et al., , b, 2011aJacquet et al., , 2015;;Planchon et al., 2013) and Pacific Ocean [START_REF] Dehairs | Barium in twilight zone suspended matter as a potential proxy for particulate organic carbon remineralization: Results for the North Pacific[END_REF].

We examined mesopelagic POC remineralisation along the GEOTRACES-GA01 section during the GEOVIDE cruise (15 May-30 June 2014; R/V Pourquoi Pas?) by assessing particulate biogenic barium (excess barium; Ba xs ) contents. This study is the first one to report the use of the Ba xs proxy in the North Atlantic. Regional variations in the Ba xs distributions along the crossed biogeochemical provinces are discussed regarding the stage and intensity of the bloom, the phytoplankton community structure, and the physical forcing. We reassessed the algorithm between Ba xs content and oxygen consumption developed for the Southern Ocean, adapting it for the North Atlantic. We compared the remineralisation fluxes resulting from this new North Atlantic-specific algorithm with those obtained using other methods in the same area. This comparison, in combination with surface primary production (PP) and POC export estimates [START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue), allowed us to evaluate the fate of POC to the deep ocean and to constrain the BCP in the North Atlantic.

Methods

Study area

The GEOVIDE section (15 May-30 June 2014; R/V Pourquoi pas?) crossed different biogeochemical provinces in the North Atlantic including the North Atlantic subtropical gyre (NAST; Stations 1 and 13), the North Atlantic drift (NADR) covering the West European (Stations 21 and 26) and Icelandic (Stations 32 and 38) basins, and the Atlantic Arctic (ARCT) divided between the Irminger (Stations 44 and 51) and Labrador (Stations 64, 69, and 77) seas [START_REF] Longhurst | Seasonal cycles of pelagic production and consumption[END_REF]Figs. 1, 2).

The evolution of chlorophyll a (chl a) concentrations from satellite imagery (Fig. 1) revealed the decline of the bloom in the NAST and the Labrador Sea and the bloom period within the NADR province and the Irminger Sea. Indeed, the highest daily PP rates were measured in the NADR and in the Irminger Sea (> 150 mmol C m -2 d -1 ; [START_REF] Fonseca-Batista | Significant N 2 fixation in productive waters of the temperate Northeast Atlantic[END_REF]this issue;[START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue). The phytoplankton community structure also varied regionally, with diatoms dominating the ARCT province and the West Euro- pean basin of the NADR, coccolithophorids dominating the Icelandic basin of the NADR, and cyanobacteria dominating the NAST province [START_REF] Tonnard | Phytoplankton assemblages along the GEOVIDE section (GEOTRACES section GA01) using CHEMTAX[END_REF]. Finally, as described elsewhere [START_REF] Daniault | The northern North Atlantic Ocean mean circulation in the early 21st century[END_REF][START_REF] García-Ibáñez | The water mass structure and transports in the Atlantic Subpolar Gyre[END_REF][START_REF] Kieke | Studies of Labrador Sea Water formation and variability in the subpolar North Atlantic in the light of international partnership and collaboration[END_REF]Zunino et al., 2017; this issue), these provinces also differ in terms of their hydrographic features. The NADR province is crossed by the subarctic front (SAF), which was located near Station 26 during GEOVIDE (Fig. 2). Strong currents were observed near the Greenland margin (probably influencing Stations 51 and 64), and an intense 1500 m deep convection happened during the winter preceding GEOVIDE in the central Labrador Sea (Station 69) due to the formation of the Labrador Sea Water (LSW) in winter (Fig. 2). These features influenced the magnitude of the carbon export fluxes, as well as the export and transfer efficiencies along the transect [START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue). The highest POC export fluxes from the upper-ocean (calculated at the depth "z" ranging from 40 to 130 m at Stations 44 and 32, respectively) were observed in the NADR province and in the Labrador Sea and reached up to 10 mmol C m -2 d -1 at Station 69 [START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue). Export efficiency (i.e. the ratio of the POC export to the PP) was generally low (≈ 10 %), except at Sta-tions 1 and 69 where it reached 35 %. The transfer efficiency (defined as the ratio of the POC export at z + 100 m to the POC export at z) was more variable, ranging from 30 % at Station 69 to 78 % at Station 32 [START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue).

Sampling and analyses

For different objectives, during GEOVIDE, suspended particles were collected using different sampling techniques. The main goal of the Niskin sampling was to derive Ba xs concentrations and, thus, carbon remineralisation fluxes in the mesopelagic zone (high resolution in the 100-1000 m layer) at stations where PP data and carbon export fluxes were also determined. The goal of the Go-Flo sampling was, at first, dedicated to the determination of all dissolved and particulate trace elements and their isotopes. Since particulate Ba and Al were determined on samples collected by both sampling techniques, we took the opportunity to compare both datasets in order to assess the quality of our data. [START_REF] Schlitzer | Ocean Data View[END_REF].

pled between the surface and 1500 m in order to cover a high vertical resolution in the mesopelagic layer (Table S1 in the Supplement). An amount of 4 to 8 L of seawater was filtered on acid-cleaned polycarbonate membranes of 0.4 µm porosity (Nuclepore ® , 47 or 90 mm diameter). Filter membranes were rinsed with Milli-Qgrade water (18.2 M cm; ≤ 5 mL) to remove sea salt, dried at ambient temperature under a laminar flow hood, and finally stored in clean petri slides until analysis in the home-based laboratory.

Filters were totally digested overnight with a concentrated tri-acid mixture (1.5 mL HCl / 1 mL HNO 3 / 0.5 mL HF; all Merck Suprapur grades) using clean Teflon vials (Savillex ® ) on a hot plate at 90 • C. The acid solution was then evaporated at 110 • C until near dryness and the residue dissolved in 13 mL 0.32 M HNO 3 (Merck; distilled Normapur). The solutions were transferred to polypropylene tubes (VWR) and analysed for barium (Ba), aluminium (Al), and other major and minor elements using an inductively coupled plasma quadrupole mass spectrometer (ICP-QMS; X Series 2 Thermo Fisher) equipped with a collision cell technology (CCT). We used a concentric quartz nebulizer (1 mL min -1 ) and nickel sample and skimmer cones. During the analyses, internal standards (Ru, In, Re, and Bi) were added to samples in order to monitor and correct the instrumental drift and matrix-dependent sensitivity variations.

Two multi-element artificial standard solutions were prepared for external calibration. The first contained major elements (Na, Mg, Al, Ca, and Ti) and the second was prepared with minor elements (Sr, Ba, REEs, Th, and U). Standards were prepared with dilution of the multi-element mixed standard stock solutions to span the expected range of sample concentrations, with concentrations in the standard curve spaced to cover potential sample variations.

The accuracy and precision of our analyses were assessed using the following Certified Reference Materials (CRM): BHVO-1, JB-3, JGb-1, and SLRS-5 (Table 1).

The detailed procedure for sample preparation and analysis is given in [START_REF] Cardinal | Geochemistry of suspended particles in the Subantarctic and Polar Frontal zones south of Australia: Constraints on export and advection processes[END_REF].

2. Ba xs concentrations measured in samples collected using the trace metal clean rosette, equipped with 22 12 L Go-Flo bottles at higher spatial resolution (31 stations) but with lower vertical resolution in the mesopelagic layer. Details about filtration, sample processing, and analyses can be found in [START_REF] Gourain | Particulate trace metals along the GEOVIDE section[END_REF]this issue). Briefly, at each depth, two size fractions were investigated: 0.45-5 µm using polysulfone filters (Supor ® ) and > 5 µm using mixed ester cellulose filters (Millipore ® ). Between 2 and 5 L of seawater was filtered for the upper water column (surface to 100 m) and 10 L for depths exceeding 100 m. Excess seawater from the filters was drawn off and then the filters were frozen in acid-cleaned petri dishes until home analysis. In the laboratory, filters were digested with a solution of 8 M HNO 3 (Ultrapur grade, Merck) and 2.3 M HF (Suprapur grade, Merck). Vials were then refluxed at 130 • C on a hotplate for 4 h. After gentle evaporation, the residue was redissolved with approximately 2 mL of 0.32 M HNO 3 spiked with 1 µg L -1 of indium. Solutions were analysed using a sector field inductively coupled plasma mass spectrometer (Element 2, Thermo) following the method of [START_REF] Planquette | Sampling for particulate trace element determination using water sampling bottles: methodology and comparison to in situ pumps[END_REF]. Total Ba and Al concentrations were calculated by summing the two size fractions. The accuracy Table 1. Particulate barium (Ba) and aluminium (Al) concentrations and resulting recoveries of the certified reference materials SLRS-5 (river water), BHVO-1 (basalt powder), JB-3 (basalt powder), and JGb-1 (gabbro powder).

Ba Al

SLRS-5 (µg kg -1 ) 13 ± 1 47 ± 2 n = 4 95 % 95 % BHVO-1 (µg g -1 ) 129 ± 1 70118 ± 984 n = 4 93 % 96 % JB-3 (µg g -1 ) 229 ± 13 92144 ± 1620 n = 4 94 % 101 % JGb-1 (µg g -1 ) 68 ± 15 91491 ± 732 n = 4 106 % 99 %
and precision of these analyses were assessed using the BCR-414 CRM (see [START_REF] Gourain | Particulate trace metals along the GEOVIDE section[END_REF]this issue).

For both Niskin and Go-Flo samples, the Ba xs concentrations were calculated by subtracting the particulate lithogenic barium (pBa-litho) from the total particulate barium (pBa). The pBa-litho was determined by multiplying the particulate aluminium (pAl) concentration by the upper continental crust (UCC) Ba : Al molar ratio (0.00135 mol mol -1 ; [START_REF] Taylor | The continental crust: its composition and evolution[END_REF].

Along the GEOVIDE section, the pBa-litho fraction represented less than 7 % of total barium, except at Stations 1 and 53, close to the Iberian and Greenland margin where pBa-litho accounted for 28 and 44 % of total Ba, respectively. Because of the rather large uncertainty associated with the UCC Ba : Al ratio and because of the high lithogenic particle loads at Stations 1 and 53, those stations were not considered further in this study. Uncertainties on Ba xs concentrations were estimated using error propagation and ranged between 6 and 25 %.

For stations where total pBa and pAl concentrations were available at similar depths, the regression of Ba xs concentrations (100-1000 m layer) from the Go-Flo samples vs. those of the Niskin samples was significant (regression slope: 0.87; R 2 : 0.61; p<0.01; n = 66; Fig. S1 in the Supplement) despite some discrepancies, especially in the higher concentration domain. Such discrepancies could have resulted from differences in the chemical protocols and most likely the filters used. The Niskin samples collected on 0.4 µm polycarbonate filters were digested using a tri-acid mix (50 % HCl / 33 % HNO 3 / 17 % HF), while the Go-Flo samples collected on paired 0.45 µm polyestersulfone / 5 µm mixed ester cellulose filters were digested using a 50 % HNO 3 / 10 % HF acid mix. The use of different filter types has been shown to lead to different concentrations, depending on the element of consideration, despite using the same digestion technique [START_REF] Planquette | Sampling for particulate trace element determination using water sampling bottles: methodology and comparison to in situ pumps[END_REF]. The addition of HCl has been shown to not improve elemental recoveries of marine particles [START_REF] Ohnemus | Cycling of Lithogenic Marine Particulates in the US GEOTRACES North Atlantic Transect[END_REF] but the larger HF con-centration of the tri-acid mix used for digesting the Niskin samples, likely, dissolved more of the refractory particles, explaining the slightly higher concentrations obtained for these samples.

In addition, filtered suspended matter was also analysed using a field emission scanning electron microscope (FE-SEM; JEOL JSM-7100F) to detect the presence of barite particles (Fig. 3). Because of time-consuming analyses, seven filters from the different basins were scanned: six samples with high mesopelagic Ba For each sample, a filter surface of 0.5 cm 2 was analysed.

To verify the relationship between Ba xs and barite particles (see Sect. 4.1), we evaluated the contribution of the barite particles to Ba xs concentrations for the sample collected at 600 m at Station 69 and which has a high mesopelagic Ba xs content (see Sect. 3). Using the FE-SEM, 0.003 % of the total filter surface was scanned and size and volume of all detected barite particles present in this surface area were assessed. To this aim, each barite particle was pictured using a magnification setting between 12 000 and 15 000×. Images were then analysed with the software ImageJ and, for each barite particle, the longest and shortest axes were measured and pixels were converted to nanometres. Barite particles were then assimilated to ellipses to calculate their volume. Finally, the concentration of pBa of each barite particle was calculated using Eq. (1):

pBa in barite = V × µ BaSO 4 × M Ba /M BaSO 4 /V SW , ( 1
)
where V is the volume of the barite particle (between 0.01 and 3.96 µm 3 ), µ BaSO 4 is the density of barite (4.45 g cm -3 ), M Ba /M BaSO 4 is the molar proportion of barium in BaSO 4 (0.59), and V SW is the volume of seawater filtered (equivalent to 0.2 mL for the portion of filter surface analysed).

3 Results

Particulate biogenic Ba xs distribution

Section overview

The Ba xs longitudinal section of concentrations (Fig. 4) shows elevated concentrations between 100 and 1000 m, in the mesopelagic layer (333

±

Individual profiles

In this section, we only present Ba xs concentrations obtained from Niskin bottles. This is because (i) both data sets converge (regression slope: 0.87; R 2 : 0.61; p<0.01; n = 66);

(ii) Niskin casts had a better sampling resolution in the 100-1000 m layer; and (iii) Niskin casts were also used for the determination of POC export using the 234-Th deficit method (see [START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF]this issue).

All the vertical Ba xs profiles (Fig. 5) show increased concentrations between 100 and 1000 m, followed by lower concentrations deeper that tend to return to a background level of 180 ± 54 pmol L -1 (n = 10) as the average along the GEO-VIDE transect. This background value is quite characteristic for the deep ocean (> 1000 m) and is considered to represent the residual Ba xs left over after partial dissolution and sedimentation of Ba xs produced during previous phytoplankton growth events [START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF].

In the NAST province (Station 13), the Ba xs concentrations steadily increased from the surface to 400 m, reaching 961 pmol L -1 , then decreased with depth, reaching the background level of 180 pmol L -1 at 1500 m.

In the West European basin of the NADR province, vertical profiles of Ba xs were similar, yet concentrations in the mesopelagic layer were smaller at Station 21 with Ba xs peaking only at 524 pmol L -1 . Ba xs concentrations returned to the background value at 1200 m. Ba xs concentration in surface waters of Station 26 were the highest of the entire section reaching 1888 pmol L -1 at 50 m (note that the value for the Go-Flo sample at 35 m reaches 24 643 pmol L -1 ; Sect. 3.1.1). Below this depth, Ba xs concentrations decreased back to the background level at 100 m, then increased again, with a second peak of 451 pmol L -1 at 200 m. In the Icelandic basin of the NADR province, Ba xs concentrations were relatively high, reaching 646 and 711 pmol L -1 at 200 and 300 m at stations 32 and 38, respectively. Station 38 was characterized by a double Ba xs peak at 300 and 700 m. Below this second maximum, Ba xs concentrations decreased to the background level at 1000 m for both stations.

In the ARCT province, a similar double peak profile was observed at Station 44, in the Irminger Sea, with Ba

Mesopelagic Ba xs

The Ba xs concentrations were integrated (trapezoidal integration) over two depth intervals of the mesopelagic layer (100-500 m and 100-1000 m; Table 2) to obtain depth-weighted average (DWA) Ba xs values.

The DWA Ba xs values between 100 and 500 m ranged from 399 to 672 pmol L -1 and from 315 to 727 pmol L -1 between 100 and 1000 m (Stations 51 and 69, respectively). The DWA Ba xs values varied by less than a factor of 1.4 between both modes of integration. Only for the Labrador Sea (Stations 64, 69, and 77) were the DWA Ba xs values for the 100-1000 m larger than for the 100-500 m interval. For the latter stations, the Ba xs inventories for the interval between 100 m and the depths where concentrations decreased to background level (1300, 1700, and et al., 1997). The approximate depth range of the major water masses is also indicated in blue shading. Similarly, low median DWA Ba xs contents were observed within the NADR province (403 ± 34 pmol L -1 , n = 4), with the lowest DWA Ba xs observed at Station 26 (391 ± 58 pmol L -1 ).

The highest median DWA Ba xs value was observed in the ARCT province (566 ± 155 pmol L -1 , n = 5). There, the DWA Ba xs contents were more variable between stations, ranging from 315 pmol L -1 at Station 51 to 727 pmol L -1 at Station 69, with a high DWA Ba xs also observed at Station 44 (633 pmol L -1 ).

Discussion

Barite is the main carrier of Ba xs

Several barite particles were observed associated with or in close proximity to biogenic fragments such as coccoliths (Fig. 3), suggesting they were originally formed inside biogenic micro-environments as proposed by others [START_REF] Bishop | The barite-opal-organic carbon association in oceanic particulate matter[END_REF][START_REF] Dehairs | Discrete suspended particles of barite and the barium cycle in the open ocean[END_REF][START_REF] Stroobants | Barite formation in the Southern Ocean water[END_REF]. However, no barite crystals were observed in the surface waters at Station 26 where very high Ba xs concentrations were recorded (up to 1888 pmol L -1 ), most likely the result of Ba uptake and/or adsorption by biota, as reported by [START_REF] Sternberg | Barium uptake and adsorption in diatoms[END_REF] for culture experiments. This result was expected as it fits in the concept of barite formation proposed by [START_REF] Stroobants | Barite formation in the Southern Ocean water[END_REF], showing that the barium sulfate in biogenic aggregates of surface waters is not crystallised, whereas below this surface layer, when organic matter degradation occurs, barite is present as discrete micron-sized particles.

Regarding the contribution of the barite particles to the Ba xs concentration for the sample collected at 600 m of Station 69, we assumed that the small filter portion (only 0.003 % of the total filter surface was analysed by FE-SEM) is representative of the whole filter. The Ba concentration deduced from the FE-SEM particle sizing analysis and using Eq. ( 1) is 1260 pmol L -1 . This is of the same order of magnitude, although 1.5 times larger, as the concentration of total Ba xs obtained by ICP-MS (831 pmol L -1 ) after whole filter digestion. The similarity between both values is remarkable considering the limitations of the FE-SEM procedure (the very small fraction of filter analysed). This also confirms that Ba xs in the mesopelagic layer is carried mostly by barite particles, as observed earlier [START_REF] Dehairs | Discrete suspended particles of barite and the barium cycle in the open ocean[END_REF].

Factors influencing the DWA Ba xs in the North Atlantic

Influence of the intensity and stage of the bloom

We compared our Ba xs inventories with the average biomass development from January to June 2014 (Fig. 6), which covers the entire productive period in the North Atlantic, starting from winter till the date of our sampling.

Along the GEOVIDE transect, the most productive area was clearly the Labrador Sea of the ARCT province, where chl a concentrations averaged 6 mg m -3 (Fig. 6). This basin was sampled during the decline of the bloom (Fig. 1; chl a concentration was > 3 mg m -3 1 month before the sampling, and PP and nutrient concentrations were low during sampling). The high DWA Ba xs observed in this area (Table 2) likely results from the large biological activity during the period preceding sampling. The West European basin of the NADR province, and in particular the area around Station 21, was also characterized by a high phytoplankton biomass between January and June (Fig. 6), though lower than in the Labrador Sea. Here, the bloom started in May (Fig. 1; chl a concentration ≈ 3 mg m -3 , 1 month before the sampling) and was still in progress during the sampling, as indicated by the high PP (135 mmol C m -2 d -1 ). These features can explain the lower DWA Ba xs observed at Station 21 (Table 2) compared to the Labrador Sea. The other stations of the NADR (Stations 26, 32, and 38) were sampled during the bloom development (Fig. 1) with high PP reaching 174 mmol m -2 d -1 at Station 26. The latter stations were characterized by lower DWA Ba xs values compared to other stations, pointing to a time lag between the phytoplankton bloom and the build-up of the Ba xs signal. However, this was not the case for Station 44, in the Irminger Sea of the ARCT province, which was sampled during the bloom (high PP, high chl a, and high nutrient concentrations during the sampling period) and characterized by one of the highest DWA Ba xs values, possibly reflecting an important past bloom.

As also indicated by others, the mesopelagic Ba xs signal builds up during the growth season and therefore integrates effects of past surface production [START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF][START_REF] Cardinal | Geochemistry of suspended particles in the Subantarctic and Polar Frontal zones south of Australia: Constraints on export and advection processes[END_REF][START_REF] Cardinal | Variations of carbon remineralisation in the Southern Ocean illustrated by the Baxs proxy[END_REF]. The large regional and temporal variability in the bloom development stage thus results in a large variability in the mesopelagic Ba xs signal in the North Atlantic.

Influence of water masses and physical forcing

The Labrador Sea (Stations 64, 69, and 77) had the largest Ba xs inventory coinciding with the presence of the Labrador Sea Water (LSW; potential temperature between 2.7 and 3.8 • C and salinity below 34.9; [START_REF] Harvey | 0-S relationships and water masses in the eastern North Atlantic[END_REF][START_REF] Yashayaev | Hydrographic changes in the Labrador Sea, 1960-2005[END_REF] in the upper 1500 m. The LSW formation takes place in the central Labrador Sea, where convection reached ∼ 1700 m during the winter preceding GEOVIDE (Fig. 2; [START_REF] Kieke | Studies of Labrador Sea Water formation and variability in the subpolar North Atlantic in the light of international partnership and collaboration[END_REF]. The deepening of the mixed layer has been recently shown to represent a major mechanism to convey organic carbon to the mesopelagic zone (from 23 to > 100 % in high-latitude regions; Dall 'Olmo et al., 2016), supporting the carbon demand of the mesopelagic food web [START_REF] Burd | Assessing the apparent imbalance between geochemical and biochemical indicators of meso-and bathypelagic biological activity: What the @$#! is wrong with present calculations of carbon budgets?[END_REF][START_REF] Aristegui | Microbial oceanography of the dark ocean' s pelagic realm[END_REF]. Moreover, the highest mesopelagic prokaryotic heterotrophic abundance during GEOVIDE was observed in the central Labrador Sea (Station 69), reaching 896 cells µL -1 at 500 m, while the median values at the other stations for which bacterial cell numbers were available for the mesopelagic zone (Stations 13,21,26,32,and 38), reached only 258 ± 60 cells µL -1 at similar depth (Julie Laroche, Jennie-Marie Ratten and Ryan Barkhouse, personal communication, 2017). Therefore, the LSW subduction area appears to reinforce the microbial loop by increasing the layer in which the bacteria can thrive feeding on increased food supplies. This condition appears to increase the Ba xs inventory.

The LSW was also present in the Irminger Sea between 500 and 1000 m at Station 44 (Fig. 2). In the temperaturesalinity plot, the high Ba xs concentrations of the second peak (823 pmol L -1 at 700 m; Fig. 5) are clearly associated with the presence of LSW (Fig. 7a), suggesting that this second deeper Ba xs maximum represents an advected signal. We cal-culated the DWA Ba xs without taking into account the second peak (100-600 m depth interval) and subtracted it from the total DWA Ba xs (100-1000 m depth interval) to estimate the advected signal. At Station 44, the advected Ba xs signal represents some 89 pmol L -1 , or 14 %, of the total signal, which is within the uncertainty of the remineralisation flux calculation (see later below). Similarly, at Station 32 the temperature-salinity plot (Fig. 7b) points out that the second Ba xs peak (at 450 m; Fig. 5) was related to the presence of the Subarctic Intermediate Water (SAIW; temperature of 5.6 ± 0.1 • C and salinity of 34.70 ± 0.02; [START_REF] Alvarez | Physical and biogeochemical transports structure in the North Atlantic subpolar gyre[END_REF], which contributes 14 pmol L -1 (3 %) of the total signal. Association of Ba xs maxima with water masses is not always clear, as is evident from the case of Station 38 where the second Ba xs maximum at 700 m (Fig. 5) does not coincide with a specific water mass (Fig. 7). In this case, the deep Ba xs maximum may possibly result from remineralisation generated by larger or heavier organic aggregates reaching greater depths. At the remaining stations there was no evidence of water mass influence. Overall, lateral transport influencing the local Ba xs distributions was observed at two stations during GEOVIDE but did not significantly modify the magnitude of the local mesopelagic Ba xs inventory. However, the subduction occurring in the Labrador Sea resulted in larger mesopelagic DWA Ba xs , probably due to high organic export and associated prokaryotic heterotrophic activity in these areas.

Influence of the phytoplankton community structure

The different Ba xs inventories may also be influenced by the differences in phytoplankton community compositions.

The ARCT province was dominated by diatoms (median value: 63 ± 19 % of the total phytoplankton community taxa; [START_REF] Tonnard | Phytoplankton assemblages along the GEOVIDE section (GEOTRACES section GA01) using CHEMTAX[END_REF] and was characterized by the highest DWA Ba xs values while the NAST and NADR provinces were characterized by higher abundance of haptophytes (median value: 43 ± 16 % of the total phytoplankton community taxa; [START_REF] Tonnard | Phytoplankton assemblages along the GEOVIDE section (GEOTRACES section GA01) using CHEMTAX[END_REF] and by lower Ba xs inventories. Coccolithophorids are part of the haptophyte family and their dominance was confirmed by visual observations on filters (surface down to 400 m) using a FE-SEM. Calcifiers, such as coccolithophorids, have been shown to be more efficient in transferring carbon to the deep ocean compared to diatoms [START_REF] Francois | Factors controlling the flux of organic carbon to the bathypelagic zone of the ocean[END_REF][START_REF] Klaas | Association of sinking organic matter with various types of mineral ballast in the deep sea: Implications for the rain ratio[END_REF][START_REF] Lam | The dynamic ocean biological pump: Insights from a global compilation of particulate organic carbon, CaCO 3 , and opal concentration profiles from the mesopelagic[END_REF]. This difference could result from the low compaction or the high fluffiness of diatom aggregates, the high degree of degradability of organic compounds within diatom aggregates, the greater density of calcite, the resistance of calcite to grazing, and the more refractory nature of the exported organic matter associated with calcite [START_REF] Bach | Influence of plankton community structure on the sinking velocity of marine aggregates[END_REF][START_REF] Francois | Factors controlling the flux of organic carbon to the bathypelagic zone of the ocean[END_REF][START_REF] Klaas | Association of sinking organic matter with various types of mineral ballast in the deep sea: Implications for the rain ratio[END_REF][START_REF] Lam | The dynamic ocean biological pump: Insights from a global compilation of particulate organic carbon, CaCO 3 , and opal concentration profiles from the mesopelagic[END_REF]Le Moigne et al., 2013a;[START_REF] Ragueneau | Si and C interactions in the world ocean: Importance of ecological processes and implications for the role of diatoms in the biological pump[END_REF]. Therefore, enhanced particle degradation when diatoms are predominant seems to increase the mesopelagic DWA Ba xs .

Relationship between Ba xs and carbon remineralisation in the North Atlantic

In previous studies focusing on the Southern Ocean, Ba xsbased mesopelagic carbon remineralisation fluxes were estimated using Eq. ( 2), which relates the accumulated mesopelagic Ba xs inventory to the rate of oxygen consumption [START_REF] Shopova | A simple model of biogeochemical element distribution in the oceanic water column[END_REF][START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF]:

mesopelagic Ba xs = 17 200 × J O 2 + Ba residual , (2)
where mesopelagic Ba xs is the DWA in the mesopelagic layer (pmol L -1 ), J O 2 is the rate of oxygen consumption (µmol L -1 d -1 ), and Ba residual is the deep-ocean Ba xs value observed at zero oxygen consumption (or Ba xs background signal), which was determined to reach 180 pmol L -1 (Dehairs et al., 1997).

The oxygen consumption JO 2 can be converted into a C remineralisation flux through Eq. ( 3):

POC mesopelagic remineralisation = Z × J O 2 × (C : O 2 ) Redfield Ratio , (3)
where the POC mesopelagic remineralisation is in mmol C m -2 d -1 , Z is the thickness of the layer in which the mesopelagic Ba xs is calculated, J O 2 is the rate of oxygen consumption given by Eq. (2), and (C : O 2 ) Redfield Ratio is the stoichiometric molar ratio of carbon to dioxygen (127/175; [START_REF] Broecker | Source and flow patterns of deep ocean[END_REF]. However, it is of interest to investigate if this relationship can be applied in the North Atlantic. Therefore, we determined the oxygen utilisation rate (OUR; µmol kg -1 yr -1 ), which is obtained by dividing the apparent oxygen utilisation (AOU, µmol kg -1 ) by the water mass age (Table S1). From the Iberian coast to Greenland, the age calculation was based on the CFC-12 distribution (when available, otherwise CFC-11) determined in 2012 (OVIDE CARINA cruise; de la [START_REF] De La Paz | Ventilation versus biology: What is the controlling mechanism of nitrous oxide distribution in the North Atlantic?[END_REF]. For the Labrador Sea, the mean age of LSW has been estimated by Rhein et al. (2015) based on a 25-year record of CFC contents. The OUR was then integrated over the 100-1000 m layer. The resulting regression between DWA Ba xs and OUR is as follows (see Fig. 8): mesopelagic Ba xs = 23391(±6368)

× J O 2 + 247(±61), (4)
where mesopelagic Ba xs and J O 2 are defined in Eq. ( 2).

Here, Ba residual is 247 pmol L -1 . This regression is significant (R 2 = 0.63; p value = 0.006) when Station 44 is excluded. This latter station was located in the Irminger Gyre (Zunino et al., 2017; this issue) and it is possible that the gyre system induced an accumulation and retention of mesopelagic Ba xs , which then no longer reflects remineralisation associated with the present growth season. Figure 8 also shows the oxygen consumption related to the GEOVIDE Ba xs values using the Southern Ocean regression (Eq. 2). It appears that for a given mesopelagic Ba xs inventory the oxygen consumption is smaller when using the Southern Ocean regression. However, both regressions are not significantly different when taking into account the errors associated with the slope and intercept of the regression in Eq. ( 4). The Southern Ocean regression appears to represent a lower limit that seems to overestimate the remineralisation fluxes. Furthermore, the relationship here deduced for the North Atlantic is sensitive to potential errors. Indeed, calculation of OUR has been shown to underestimate the ocean respiration because of the non-proportional diffusive mixing of AOU and water mass age resulting in an excess loss of AOU versus age [START_REF] Koeve | Oxygen utilization rate (OUR) underestimates ocean respiration: a model study[END_REF]. This would decrease the mismatch between the Southern Ocean and North Atlantic regressions. Errors can also be directly associated with the CFC-based age values of the water masses, which would appear especially critical for LSW. Indeed, the severe winter preceding the cruise (2013/2014) appeared to have strongly ventilated LSW with a mixed-layer depth exceeding 1700 m [START_REF] Kieke | Studies of Labrador Sea Water formation and variability in the subpolar North Atlantic in the light of international partnership and collaboration[END_REF], indicating that the mean age (4 years) estimated by Rhein et al. (2015) may have overestimated the real LSW age (Pascale Lherminier, personal communication, 2017). Moreover, in the Labrador Sea, the residence time of LSW strongly varies between the central Labrador Sea (4-5 years) and the boundary currents off the Greenland and Newfoundland coasts (a few months; [START_REF] Deshayes | Formation and export of deep water in the Labrador and Irminger Seas in a GCM[END_REF][START_REF] Straneo | Spreading of Labrador sea water: an advective-diffusive study based on Lagrangian data[END_REF]. An overestimation of these ages leads to underestimating OUR, resulting in reducing the apparent discrepancy between the both North Atlantic and Southern Ocean regressions.

In the following discussion, carbon remineralisation fluxes are estimated for the North Atlantic (GEOVIDE and GEOSECS cruises) using Eqs. ( 4) and (3). The GEOVIDE remineralisation fluxes are compared with values reported for the World Ocean and also based on Ba xs inventories (Table 3; Fig. 9). In the North Atlantic, the fluxes obtained during the GEOVIDE and GEOSECS (symbolised by stars in Fig. 9) cruises are of the same order of magnitude, highlighting a relatively constant remineralisation over the last 44 years. The remineralisation fluxes reported for the Southern and Pacific oceans are similar to those in the NAST and NADR provinces of the North Atlantic. However, the fluxes in the ARCT province are clearly higher, highlighting an important remineralisation in the northern part of the North Atlantic compared to other oceans.

Remineralisation from direct measurements

In the North Atlantic, carbon respiration rates were also deduced from surface drifting sediment traps and associated shipboard incubations (Fig. 9). [START_REF] Collins | The multiple fates of sinking particles in the North Atlantic Ocean[END_REF] determined very high respiration rates reaching 39 and 72 mmol C m -2 d -1 at sites located in the NADR and in the ARCT provinces, respectively. Nevertheless, these high fluxes were deduced in the upper mesopelagic layer (50-150 m) where respiration is larger compared to the lower mesopelagic layer (100-1000 m). This difference in depth interval could thus explain the lower remineralisation rates in our study. Also using surface drifting sediment traps and associated-shipboard incubations but supplemented by measurements of zooplankton respiration, Giering et al. (2014) determined respiration rates in the NADR province (PAP site) reaching 7.1 mmol C m -2 d -1 during summer. This flux, determined over the 50-1000 m depth interval, is of the same order of magnitude as our estimates for the NADR province.

Remineralisation from the deep sediment traps

The remineralisation flux in the mesopelagic layer can also be derived from the difference between a deep POC export flux and a surface POC export flux. [START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF] compiled deep POC fluxes from bottom-tethered sediment traps and calculated the corresponding export production (upperocean POC export flux) using an ecosystem model [START_REF] Laws | Temperature effects on export production in the open ocean[END_REF] for most world ocean provinces. Then, by difference, the authors estimated an annual average of carbon remineralisation fluxes in the mesopelagic layer, which were converted into daily average fluxes. Remineralisation fluxes reached values of 34 mmol C m -2 d -1 in the ARCT province, 9 mmol C m -2 d -1 in the NADR province, and 4 mmol C m -2 d -1 in the NAST province (Fig. 9). Note that the flux in the ARCT province was one of the highest mesopelagic remineralisation fluxes estimated worldwide, confirming the occurrence of important remineralisation in the northern part of the North Atlantic as compared to other [START_REF] Dehairs | Barium in twilight zone suspended matter as a potential proxy for particulate organic carbon remineralization: Results for the North Pacific[END_REF] and the Southern Ocean [START_REF] Cardinal | Variations of carbon remineralisation in the Southern Ocean illustrated by the Baxs proxy[END_REF]Jacquet et al., 2008aJacquet et al., , b, 2011bJacquet et al., , 2015;;Planchon et al., 2013) were calculated based on the Ba xs inventories. The insert shows data for the North Atlantic: sites indicated by circles lined in black are from the present study; at sites labelled with # symbols remineralisation was deduced from POC fluxes recorded by moored sediment traps [START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF]; at sites labelled by • remineralisation was obtained from on-board incubations [START_REF] Collins | The multiple fates of sinking particles in the North Atlantic Ocean[END_REF]Giering et al., 2014); sites labelled with * are GEOSECS sites for which we calculated remineralisation from existing Ba xs profiles (Brewer et al., unpublished results). Data were plotted using the ODV software [START_REF] Schlitzer | Ocean Data View[END_REF].

oceans. The values published by [START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF] for the North Atlantic are quite similar to our median values obtained during GEOVIDE. Indeed, mesopelagic remineralisation fluxes based on the Ba xs proxy were similar to the value reported by [START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF] for the NAST province, while they were respectively 2-and 4-fold lower in the NADR and in the ARCT provinces.

Overall, the remineralisation fluxes deduced from the Ba xs proxy are in concordance with those obtained using the other methods, confirming the order of magnitude of the mesopelagic remineralisation fluxes determined in this study of the North Atlantic (Fig. 9).

The biological carbon pump in the North Atlantic

In order to investigate the efficiency of the biological carbon pump in the North Atlantic, we examined the daily PP [START_REF] Fonseca-Batista | Significant N 2 fixation in productive waters of the temperate Northeast Atlantic[END_REF]this issue;[START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF] this issue), the upper-ocean POC export (Lemaitre et al., 2008;this issue), and the POC remineralisation in the mesopelagic layer (Table 4; Fig. 10).

During GEOVIDE, low (≤ 14 %) export efficiencies (i.e. the ratio between PP and POC export) were observed at most stations, indicating an accumulation of biomass in surface waters or a strong turnover of the exported organic matter due to important remineralisation occurring in the upper water column (< 100 m). Furthermore, mesopelagic POC remineralisation fluxes were relatively high, equalling or exceed-ing the POC export fluxes at some stations. This highlights a strong mesopelagic remineralisation with little or no material left for export to the deep ocean, but above all, it involves an imbalance between carbon supplies and mesopelagic remineralisation.

This imbalance can result from the differences between the time windows over which the PP, POC export, and POC remineralisation fluxes are integrated. Indeed, the measurements of PP represent a snapshot (24 h incubations) while measurements of export ( 234 Th) integrate several weeks [START_REF] Benitez-Nelson | A time-series study of particulate matter export in the North Pacific Subtropical Gyre based on 234 Th : 238 U disequilibrium[END_REF][START_REF] Buesseler | Carbon and nitrogen export during the JGOFS North Atlantic Bloom experiment estimated from 234 Th : 238 U disequilibria[END_REF] and remineralisation (from the Ba xs proxy) probably integrates much longer timescales. Moreover, previous studies in the Southern Ocean showed that mesopelagic processing of exported organic carbon, as reflected by Ba xs , has a phase lag relative to the upper-ocean processes [START_REF] Dehairs | Particulate barium stocks and oxygen consumption in the Southern Ocean mesopelagic water column during spring and early summer: Relationship with export production[END_REF][START_REF] Cardinal | Variations of carbon remineralisation in the Southern Ocean illustrated by the Baxs proxy[END_REF]. Thus, we do not expect mesopelagic Ba xs to be in phase with coinciding amplitude of PP and subsequent export. Because of the observed high remineralisation fluxes relative to the export fluxes, particularly in the ARCT province, it is likely that particulate organic matter sank out of the surface waters and became subject to mineralisation in the mesopelagic layer during the period preceding the specific time windows for POC export and PP. Such discrepancies between fluxes can be amplified by the spatial and temporal variability in the phytoplankton blooms in this province, generating sudden high export events and associated remineralisation. In contrast to the above, at Station 32 in the NADR province, a large fraction of exported POC (50 %; Table 4) appears to escape remineralisation. The more efficient POC transfer through the mesopelagic layer of this province may be explained by the fact that sampling took place in an early stage of the bloom and/or by the presence of calcified phytoplankton species, ballasting aggregates thereby increasing their settling velocity (see Sect. 4.2.3).

Overall, the remineralisation in the mesopelagic layer is an important process that needs to be taken into account as our results point to the poor capacity of specific areas within the North Atlantic to sequester carbon at depths below 1000 m in spring 2014.

Conclusions

We investigated mesopelagic carbon remineralisation fluxes in the North Atlantic during spring 2014 (GEOVIDE section) using, for the first time, the particulate biogenic barium inventories measured for this area. The excess biogenic barium (Ba xs ) inventories in the mesopelagic layer varied between the different provinces of the North Atlantic. The largest Ba xs inventory was observed in the ARCT province, where high carbon production rates were also observed earlier in the season. The regional variations in the Ba xs inventory may also result from differences with phytoplankton community composition encountered along this transatlantic section. Lower mesopelagic Ba xs contents occurred where smaller calcified phytoplankton species dominated, such as in the NADR province. Finally, the ARCT province was also characterized by important water mass subduction, generating a large transport of organic matter to the deep ocean, which might have resulted into an important Ba xs accumulation in the mesopelagic layer.

Using the OUR method, we confirmed that the mesopelagic Ba xs inventory can be related to the oxygen utilisation rate, but the relationship between both parameters is slightly different compared to the relationship proposed elsewhere for the Southern Ocean. A new relationship is thus proposed for the North Atlantic. This proxy approach provided estimations of mesopelagic remineralisation fluxes of similar magnitude to those obtained by others using independent methods (free-floating and moored sediment traps, incubations) in the North Atlantic.

Overall, in spring 2014, the mesopelagic remineralisation balanced or exceeded POC export in the subtropical and subpolar provinces of the North Atlantic, highlighting the important impact of the mesopelagic remineralisation on the biological carbon pump and indicating that little to no organic matter was transferred below 1000 m in this region.

Figure 1 .

 1 Figure 1. Satellite-derived chlorophyll a concentrations (MODIS Aqua from http://giovanni.sci.gsfc.nasa.gov/, last access: September 2017), in milligrams per cubic metre, during the GEOVIDE cruise (May and June 2014). The GEOVIDE transect (grey line) and the main crossed provinces are indicated. NAST: North Atlantic subtropical gyre; NADR: North Atlantic drift; ARCT: Atlantic Arctic. Coloured circles indicate stations sampled at the corresponding month.

Figure 2 .

 2 Figure 2. (a) Schematic of the circulation features, adapted from García-Ibáñez et al. (2015). Bathymetry is plotted in colour with interval boundaries at 100 m, at 1000 m, and every 1000 m below 1000 m. The red and green arrows represent the main surface currents, the pink and orange arrows represent currents at intermediate depths, and the blue and purple arrows represent the deep currents. Diamonds indicate station positions. The approximate locations of the subarctic front (SAF; black bar crossing Station 26) and the formation site of the Labrador Sea Water (LSW form.) are indicated. (b) Salinity along the GEOVIDE section and associated water masses are as follows: LSW: Labrador Sea Water; ISOW: Iceland-Scotland Overflow Water; IcSPMW: Iceland Subpolar Mode Water; SAIW: Subarctic Intermediate Water; NACW: North Atlantic Central Water; MW: Mediterranean water; DSOW: Denmark Strait Overflow Water; NEADW: North East Atlantic Deep Water. Station labels in red indicate sites where Niskin casts were sampled. NAST: North Atlantic subtropical gyre; NADR: North Atlantic drift; ARCT: Atlantic Arctic. Data were plotted using ODV software[START_REF] Schlitzer | Ocean Data View[END_REF].

  xs concentrations (Station 13 at 400 m; Station 38 at 300 m; Station 44 at 300 and 700 m; Station 69 at 600 m, and Station 77 at 300 m) and one sample with high surface Ba xs concentrations (Station 26 at 50 m).

Figure 3 .

 3 Figure 3. Barite particles observed with the FE-SEM at (a) Station 38 (300 m), (b, c) Station 44 (700 m), and (d) Station 69 (600 m). Panel (c) shows the backscattered electron image of the aggregate in (b) highlighting the shape of the partly hidden barite crystal. White arrows indicate the position of barite crystals.

Figure 4 .

 4 Figure 4. Section of the particulate biogenic barium (Ba xs ) (pmol L -1) determined in samples collected with the Go-Flo bottles. Stations labelled in red are those where profiles were also obtained from Niskin casts. Data were plotted using the ODV software[START_REF] Schlitzer | Ocean Data View[END_REF].

Figure 5 .

 5 Figure 5. Vertical profiles of Ba xs concentrations (pmol L -1 ) determined from Niskin casts during the GEOVIDE (squares) and GEOSECS (circles) cruises. The vertical black dashed line (at 180 pmol L -1 ) represents the deep-ocean Ba xs value (or Ba xs background signal; Dehairset al., 1997). The approximate depth range of the major water masses is also indicated in blue shading.

Figure 6 .Figure 7 .

 67 Figure 6. Map of time-averaged chlorophyll a concentrations (mg m -3 ) for the period from January to June 2014 (monthly 4 km MODIS Aqua model; http://giovanni.sci.gsfc.nasa.gov/, last access: September 2017).

Figure 8 .

 8 Figure 8. Regression of DWA mesopelagic Ba xs (pmol L -1 ) versus O 2 consumption rate (µmol L -1 d -1 ) using the Southern Ocean transfer function from Dehairs et al. (1997; red circles) and the transfer function obtained here for the North Atlantic (black circles). Station 44 (triangle) was excluded from the regression. If station 44 is included, R 2 = 0.33 and the p value = 0.07.

4. 4

 4 Comparison of remineralisation fluxes 4.4.1 Remineralisation from the Ba xs proxy

Figure 9 .

 9 Figure 9. Summary of published POC remineralisation fluxes (mmol C m -2 d -1 ) in the world ocean. The remineralisation fluxes for the Pacific Ocean[START_REF] Dehairs | Barium in twilight zone suspended matter as a potential proxy for particulate organic carbon remineralization: Results for the North Pacific[END_REF] and the Southern Ocean[START_REF] Cardinal | Variations of carbon remineralisation in the Southern Ocean illustrated by the Baxs proxy[END_REF] Jacquet et al., 2008a Jacquet et al., , b, 2011b Jacquet et al., , 2015;; Planchon et al., 2013) were calculated based on the Ba xs inventories. The insert shows data for the North Atlantic: sites indicated by circles lined in black are from the present study; at sites labelled with # symbols remineralisation was deduced from POC fluxes recorded by moored sediment traps[START_REF] Honjo | Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983[END_REF]; at sites labelled by • remineralisation was obtained from on-board incubations[START_REF] Collins | The multiple fates of sinking particles in the North Atlantic Ocean[END_REF] Giering et al., 2014); sites labelled with * are GEOSECS sites for which we calculated remineralisation from existing Ba xs profiles (Brewer et al., unpublished results). Data were plotted using the ODV software[START_REF] Schlitzer | Ocean Data View[END_REF].

Figure 10 .

 10 Figure 10. General schematic of the biological carbon pump in the NAST, NADR, and ARCT provinces during GEOVIDE. Primary production (PP) data from Fonseca-Batista et al. (2018; this issue) and Lemaitre et al. (2018; this issue); particulate organic carbon (POC) export fluxes from Lemaitre et al. (2018; this issue) and POC remineralisation fluxes from this study. The dominating phytoplankton communities and the stage of the bloom are also indicated.

  [START_REF] Gourain | Particulate trace metals along the GEOVIDE section[END_REF] this issue). Concentrations ranged from 4 (Station 11, 55 m) to 24 643 (Station 26, 35 m) pmol L -1 in surface waters and from 7 (Station 71, 350 m) to 1388 (Station 15, 300 m) pmol L -1 in the mesopelagic layer (100-1000 m). For the mesopelagic layer, where the maximum Ba xs concentrations were generally observed, the highest Ba xs concentration was observed in the NAST province, reaching 1388 pmol L -1 at 300 m of Station 15. These maxima occurred between 200 and 600 m but were spread over larger depth intervals in the ARCT province, where high Ba xs values occurred until 1200 m depth at Station 69.

	224 pmol L -1 ; median ± 1 SD;
	n = 209). In comparison, the surface ocean (depths <100 m)
	and the deep ocean (depths > 1000 m) are characterized by

lower median values (94 and 114 pmol L -1 , n = 113 and 199, respectively). Exceptions can be observed for the upper waters at Stations 25 and 26 and bottom waters at

 32, 36, 38, and 71

, where high Ba xs concentrations may be attributed to Ba assimilation by phytoplank-ton and sediment resuspension, respectively (

Table 2 .

 2 Depth-weighted average (DWA) values of mesopelagic Ba

	Province Station Latitude Longitude DWA Ba xs 100-500 m DWA Ba xs 100-1000 m
			( • N)	( • E)	(pmol L -1 )	(pmol L -1 )
	NAST	13	41.4	-13.9	578 ± 89	419 ± 71
	NADR	21	46.5	-19.7	428 ± 69	394 ± 64
		26	50.3	-22.6	405 ± 59	391 ± 58
		32	55.5	-26.7	522 ± 81	413 ± 66
		38	58.8	-31.3	572 ± 86	465 ± 78
	ARCT	44	59.6	-38.9	678 ± 104	633 ± 98
		51	59.8	-42	399 ± 72	315 ± 58
		64	59.1	-46.1	464 ± 95	566 ± 99
		69	55.8	-48.1	672 ± 111	727 ± 118
		77	53	-51.1	472 ± 80	505 ± 83

xs (pmol L -1 ) for the 100-500 m and 100-1000 m depth intervals. The biogeochemical provinces defined by

[START_REF] Longhurst | Seasonal cycles of pelagic production and consumption[END_REF]

are also indicated: NAST: North Atlantic subtropical gyre; NADR: North Atlantic drift; ARCT: Atlantic Arctic. casts at Stations 64, 69, and 77, respectively) were somewhat smaller than for the inventories between 100 and 1000 m (up to 1.5 times in the case of Station 77). To facilitate intercomparison among stations, we consistently considered Ba xs inventories over the 100-1000 m depth interval in the following discussion.

Within the NAST province, Station 13 was characterized by a relatively low DWA Ba xs value of 419 pmol L -1 .

Table 3 .

 3 Comparison of the Ba xs inventory (pmol L -1 ) and related carbon mesopelagic remineralisation fluxes (mmol C m -2 d -1 ) obtained in the world ocean. Fluxes are calculated with the new North Atlantic regression (Eq. 4) for the GEOVIDE and GEOSECS cruises and with the Southern Ocean regression (Eq. 2) for the other studies. HNLC: high nutrient-low chlorophyll; art. Fe fertilised: artificially Fe fertilised; nat. Fe fertilised: naturally Fe fertilised; PF: Polar Front; NAST: North Atlantic subtropical gyre; NADR: North Atlantic drift; ARCT: Atlantic Arctic.

	Cruise (season)
	Location
	Features
	Depth DWA Ba xs

Table 4 .

 4 Comparison of the mesopelagic POC remineralisation fluxes (Remineralisation) with primary production (PP) and POC export fluxes in the upper water column (Export). All fluxes are expressed in mmol C m -2 d -1 .

		ARCT -		ARCT -		NADR		NAST
		Labrador Sea	Irminger Sea
	Station	77 69 64	51	44 38	32	26	21	13
	PP 1	80 27 54 166	137 68 105 174 135	79
	Export 2	6 10	8	3	1	5	8	7	5	2
	Remineralisation	7 13	9	2	10	6	4	4	4	5

Biogeosciences, 15, 2289-2307, 2018 www.biogeosciences.net/15/2289/2018/

www.biogeosciences.net/15/2289/2018/ Biogeosciences, 15, 2289-2307, 2018

PP data from Fonseca-Batista et al. (2018; this issue) and[START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF] this issue).

Export data from[START_REF] Lemaitre | High variability of export fluxes along the North Atlantic GEOTRACES section GA01 -Part I: Particulate organic carbon export deduced from the 234 Th method[END_REF] this issue).

Acknowledgements. We would like to thank the captain and the crew of the R/V Pourquoi Pas?, as well as Fabien Perault and Emmanuel De Saint Léger from the CNRS DT-INSU for their help during the CTD deployments. Pierre Branellec, Floriane Desprez de Gésincourt, Michel Hamon, Catherine Kermabon, Philippe Le Bot, Stéphane Leizour, and Olivier Ménage are also acknowledged for their technical help during the cruise. We acknowledge Lorna Foliot, Raphaëlle Sauzède, Joséphine Ras, Hervé Claustre, and Céline Dimier for the sampling and analysis of pigments. We would like to thank the co-chief scientists Pascale Lherminier and Géraldine Sarthou. We also express our thanks to Pascale Lherminier, Herlé Mercier, Monika Rhein, Julie Deshayes, and Claude Talandier for providing useful help on the characteristics of the Labrador Sea Water. Satellite chlorophyll a data and visualisations used in this study were produced with the Giovanni and the Ocean Color (Ocean Biology Processing Group; OBPG) online data system, developed and maintained by NASA.

This work was funded by the Flanders Research Foundation (project G071512N), the Vrije Universiteit Brussel (strategy research program: project SRP-2), the French National Research Agency (ANR-13-BS06-0014 and ANR-12-PDOC-0025-01), the French National Center for Scientific Research (CNRS-LEFE-CYBER), IFREMER, and the "Laboratoire d'Excellence" Labex-Mer (ANR-10-LABX-19). For this work, Maribel I. Garcia-Ibáñez was supported by the Centre for Climate Dynamics at the Bjerknes Centre and by the Spanish Ministry of Economy and Competitiveness through the BOCATS (CTM2013-41048-P) project co-funded by the Fondo Europeo de Desarrollo Regional 2014-2020 (FEDER).

Edited by: Gilles Reverdin Reviewed by: two anonymous referees

Data availability. Underlying data are reported in the Supplement.

The Supplement related to this article is available online at https://doi.org/10.5194/bg-15-2289-2018-supplement.

Competing interests. The authors declare that they have no conflict of interest. Special issue statement. This article is part of the special issue "GEOVIDE, an international GEOTRACES study along the OVIDE section in the North Atlantic and in the Labrador Sea (GA01)". It is not associated with a conference.