Early winter barium excess in the southern Indian Ocean as an annual remineralisation proxy (GEOTRACES GIPr07 cruise)

. The Southern Ocean (SO) is of global importance to the carbon cycle, and processes such as mesopelagic remineralisation that impact the efﬁciency of the biological carbon pump in this region need to be better constrained. During this study early austral winter barium excess (Ba xs ) concentrations were measured for the ﬁrst time, along 30 ◦ E in the southern Indian Ocean. Winter Ba xs concentrations of 59 to 684 pmol L −


Introduction
The Southern Ocean (SO) is a carbon sink of global significance responsible for 40 %-50 % of the global oceans' carbon uptake (Friedlingstein et al., 2019;Gregor et al., 2019;Gruber et al., 2019). Oceanic carbon uptake is regulated by various processes, including the biological carbon pump (BCP). Inorganic carbon is consumed and released by photosynthetic organisms through photosynthesis and respiration (Sarmiento and Gruber, 2006), thereby regulating the earth's carbon cycle by partially sequestering photosynthetically fixed CO 2 in the ocean interior (Honjo et al., 2014). In particular, the SO BCP is a crucial contributor to the earth's carbon cycle by exporting, from surface waters, ∼ 3 Pg C yr −1 of the ∼ 10 Pg C yr −1 of global export production (Schlitzer, 2002). The efficiency of the BCP is linked to the export and preservation of surface particulate matter and is directly linked to atmospheric CO 2 levels on glacialinterglacial timescales (Honjo et al., 2014;Sigman et al., 2010).
Published by Copernicus Publications on behalf of the European Geosciences Union. Sedimentation out of the surface layer (∼ 100 m) is defined as surface export and out of the mesopelagic zone (∼ 1000 m) as deep export (Passow and Carlow, 2012). There are large gaps in our knowledge with regard to deep carbon export, internal cycling and the seasonality of these processes (Takahashi et al., 2012). The magnitude of deep carbon export is dependent on the efficiency of mesopelagic remineralisation (Jacquet et al., 2015) which can balance or even exceed particulate organic carbon (POC) surface export, especially later in the growing season, thereby limiting deep export (Buesseler and Boyd, 2009;Jacquet et al., 2011Jacquet et al., , 2015Lemaitre et al., 2018;Planchon et al., 2013). A possible explanation for imbalances between surface export and mesopelagic processes can be the lateral advection of surface waters with lower particle export relative to the mesopelagic signal (Planchon et al., 2013). It is also possible that continued remineralisation of earlier larger export fluxes is detected in the mesopelagic signal but not in the export fluxes of in situ observations (Planchon et al., 2013). In addition to this, the efficiency of remineralisation is influenced by the size and composition of exported particles (Rosengard et al., 2015;Twining et al., 2014), as well as the pathway by which these particles are transported downwards (e.g. eddy subduction, active migration, sinking or mixing) from the surface mixed layer to the mesopelagic zone (Boyd et al., 2019;Le Moigne, 2019), creating an intricate web of processes to disentangle. Mesopelagic remineralisation has also been shown to be influenced by environmental factors, such as temperature, phytoplankton community structure and nutrient availability (Bopp et al., 2013;Buesseler and Boyd, 2009). Indeed, nutrient limitation in surface waters limits export and consequently mesopelagic remineralisation by promoting the shift to smaller phytoplankton assemblages that preferentially take up recycled nutrients in the surface mixed layer (Planchon et al., 2013). Phytoplankton community composition exerts an important control in which diatoms are more efficiently exported due to their large size and ballasting by biogenic silica compared to smaller non-diatom phytoplankton (Armstrong et al., 2009;Buesseler, 1998;Ducklow et al., 2001). Latitudinal trends in remineralisation efficiency can also be linked to temperaturedependent heterotrophs that are responsible for remineralisation (DeVries and Weber, 2017;Marsay et al., 2015). The mesopelagic layer is under-studied, especially in the high latitudes, and therefore these processes are poorly constrained despite their importance to global elemental cycles, including that of carbon (Le Moigne, 2019; Robinson et al., 2010).
Export and remineralisation tracers, such as 234 Th / 238 U and apparent oxygen utilisation (AOU), have been used to study mesopelagic POC remineralisation fluxes Planchon et al., 2013;Lemaitre et al., 2018). Surface export is set by the deficit of 234 Th activities over 238 U activities, while remineralisation processes are reflected by 234 Th / 238 U ratios larger than 1 below the surface mixed layer integrating processes over a 2-to 3-week period (Bues-seler et al., 2005;Planchon et al., 2013). AOU is the depletion of oxygen (O 2 ) in the ocean interior relative to surface saturation due to biological respiration when surface water masses are subducted. AOU is dependent on salinity and temperature and integrates remineralisation on timescales of years to decades (Ito et al., 2004). Inaccuracies have, however, been detected with AOU as a remineralisation proxy, specifically in high-latitude areas due to O 2 undersaturation as a consequence of large temperature gradients (Ito et al., 2004).
Barium excess (Ba xs ) is another proxy utilised to yield estimates of mesopelagic POC remineralisation fluxes. It is defined as the "biogenic" portion of particulate barium (pBa) as barite crystals, formed by the decay of bio-aggregates below the surface mixed layer (Bishop, 1988;Dehairs et al., 1980;Lam and Bishop, 2007;Legeleux and Reyss, 1996;van Beek et al., 2007). As these crystals are released, a Ba xs peak is formed within the mesopelagic zone which has been found to correlate to primary production (PP), O 2 consumption and POC remineralisation (Dehairs and Goeyens, 1996;Dehairs et al., 1997). Depth-integrated rates of O 2 consumption between the base of the mixed layer and 1000 m were estimated using an inverse 1-D advection-diffusionconsumption model (Shopova et al., 1995) to develop transfer functions between the Ba xs signal and the rate of surface POC export for subsequent mesopelagic remineralisation (Dehairs and Goeyens, 1996;Dehairs et al., 1997). Strong correlations have been obtained between the wellestablished export and remineralisation flux proxy 234 Th and Ba xs during studies conducted in the SO and the North Atlantic, confirming the validity of Ba xs as a remineralisation proxy Lemaitre et al., 2018;Planchon et al., 2013). Estimates of POC remineralisation fluxes, using the Ba xs proxy, are directly influenced by the background signal of Ba xs after partial dissolution and sedimentation from the previous bloom season. It can be thought of as "preformed" Ba xs , defined as the Ba residual signal at zero O 2 consumption (Jacquet et al., 2015). Because studies conducted in spring and summer suggest that the mesopelagic Ba xs signal lasts between a few days to a few weeks (Dehairs et al., 1997;Jacquet et al., 2007Jacquet et al., , 2008a, it is postulated that winter measurements should give the true SO Ba residual value (Jacquet et al., 2008b(Jacquet et al., , 2011. In this context, as part of a GEOTRACES process study (GIpr07) of a transect along 30 • E in the southern Indian Ocean (58.5 to 41.0 • S), we studied Ba xs distributions during early austral winter (July 2017) to better constrain the SO Ba residual concentrations and the timescale of this proxy. To our knowledge these are the first reported wintertime values for this proxy in the SO.
During the GEOTRACES GIpr07 cruise, which took place in early austral winter (28 June-13 July 2017) on board the R/V SA Agulhas II, seven stations were sampled along 30 • E, from 58.5 to 41.0 • S (WOCE I06S, Fig. 1a). At each station between 15 and 21 samples were collected from 25 m down to 1500 m, for shallow stations, and down to 4250 m, for deep stations, to be analysed for multiple parameters.
Positions of the fronts during the cruise were determined using the July monthly mean absolute dynamic topography data from the CLS/AVISO product (Rio et al., 2011), with boundary definitions from Swart et al. (2010). From north to south the identified fronts are the subtropical front (STF), the subantarctic front (SAF), the polar front (PF), the Southern Antarctic Circumpolar Current front (SACCf) and the southern boundary (SBdy) (Fig. 1a). The marginal ice zone, identified as the position of 30 % ice cover, was positioned at 61.7 • S, approximately 3 • (356 km) south of the southernmost station (de Jong et al., 2018). Therefore, a potential sea ice influence on our study area can be disregarded.

Temperature, salinity and dissolved O 2
Temperature ( • C), salinity and dissolved O 2 (µmol L −1 ) profiles were measured by sensors (SBE 911plus) which were calibrated by the manufacturer within a year prior to the cruise. At each cast, discrete seawater samples were collected and analysed on board for in situ calibration of sensor data for salinity (8410A Portasal salinometer, R 2 = 0.99) and dissolved O 2 concentrations (Metrohm 848 Titrino plus, R 2 = 0.83; Ehrhardt et al., 1983). Temperature and salinity measurements were used to calculate potential density (σ θ ; Gill, 1982) to characterise water masses sampled and to identify the mixed layer depth (MLD). The MLD is the depth at which there is a change of 0.03 kg m −3 in σ θ from a near-surface value at ∼ 10 m (de Boyer Montégut et al., 2004). Decreases in dissolved O 2 concentrations at intermediate depths, together with Ba xs concentrations, were used to define the mesopelagic remineralisation layer.

pBa and particulate aluminium
Profile sampling of the water column was conducted with a GEOTRACES-compliant trace-metal-clean CTD (conductivity, temperature, depth) housed on an epoxy coated aluminium frame with titanium bolts equipped with 24 × 12 L trace-metal-clean Teflon-coated GO-FLO bottles (General Oceanics). All sampling and analyses were conducted following the GEOTRACES clean sampling and analysis protocols (Cutter et al., 2017). Volumes of 2 to 7 L seawater were filtered from the GO-FLO bottles onto acid-washed polyethersulfone filters (25 mm diameter, Supor, 0.45 µm pore size) for pBa and particulate aluminium (pAl) analyses.
Filters were mounted in-line on the side spigot of each GO-FLO bottle on Swinnex filter holders. Furthermore, bottles were mixed three times before filtration, as recommended by the GEOTRACES protocols (Cutter et al., 2017), to ensure homogenous sampling. Although the large fast-sinking fraction of particles may be under-sampled by using bottles (Bishop and Edmond, 1976;Planquette and Sherrell, 2012), comparing data that were generated using the same internationally validated sampling systems and protocols (Cutter et al., 2017) as we do in this study minimises potential bias. After filtration, filters were placed in trace-metal-clean petri slides (Pall) and kept frozen at −20 • C until further processing on land. Sample processing was conducted under a class 100 HEPA-filtered laminar flow and extraction hood in a clean laboratory.
The pBa and pAl samples were processed and analysed 6 months after sample collection at LEMAR (France). Unused blank filters and filters containing the samples were acid reflux digested at 130 • C in acid-cleaned Savillex vials using a mixture of HF and HNO 3 (both Ultrapur grade, Merck) solutions (Planquette and Sherrell, 2012). Archive solutions were stored in 3 mL of 0.12 M HNO 3 (Ultrapur grade), of which 250 µL was diluted up to 2 mL for analysis by sector field inductively coupled plasma mass spectrometry (SF-ICP-MS; Element XR, Thermo Scientific). Samples were spiked with 1 µg L −1 indium as an internal standard to correct for instrument drift. The detection limits, defined as 3 times the standard deviation of the blanks (unused filter blanks), were 0.39 pmol L −1 and 0.03 nmol L −1 (n = 5) for pBa and pAl, respectively. Mean amounts (in nmol) of a given element determined in unused filter blanks were subtracted from the amounts in the sample filter and then divided by the volume filtered. Three certified reference materials (BCR 414, MESS 4, and PACS 3) were processed and analysed with the samples to assess the accuracy of the methodology. Our values were in good agreement with the certified values of the reference materials (Table 1) (Jochum et al., 2005). Percentage error of analyses was determined by the repeat analysis of random samples during each run, and the mean percentage error of sample analysis for pBa and pAl was 9.2 ± 2.5 % and 11.1 ± 4.6 % (mean ± SD, n = 6), respectively.

Ba xs as a proxy for mesopelagic POC remineralisation
The non-lithogenic fraction of pBa, Ba xs , was calculated by subtracting the lithogenic fraction of pBa from the total pBa measured using Eq. (1). The lithogenic contribution to pBa was calculated by multiplying the pAl concentration with the Ba/Al upper continental crust (UCC) ratio, 0.00135, as determined by Taylor and McLennan (1985).  Total pBa and Ba xs profiles were nearly identical with a mean percentage Ba xs to total pBa of 99 ± 1 % (mean ± SD, n = 124; Table S1 in the Supplement), indicating that pBa from lithogenic sources was negligible. This ensures the accurate estimation of Ba xs , which requires that less than 50 % of pBa should be associated with lithogenic inputs (Dymond et al., 1992). The mesopelagic POC remineralisation flux was estimated using Eq. (2) (Dehairs and Goeyens, 1996;Shopova et al., 1995): where the mesopelagic POC remineralisation flux is expressed in milligrams carbon per square metre per day (mg C m −2 d −1 ), Z is the depth range of the mesopelagic Ba xs layer (100-1000 m), C : O 2 is the stoichiometric molar ratio of carbon to O 2 consumption by remineralisation as per the Redfield ratio (127/175; Broecker et al., 1985), 12.01 is the molar mass of carbon (g mol −1 ), and J O 2 is the rate of O 2 consumption (µmol L −1 d −1 ) as estimated using Eq. (3).
Equation (3) (Dehairs and Goeyens, 1996;Shopova et al., 1995) is the linearisation of the exponential function by Dehairs et al. (1997). Mesopelagic Ba xs is the depth-weighted average Ba xs of the mesopelagic zone (pmol L −1 ), the constant value of 17 200 is the slope of the linear regression of depth-weighted average Ba xs (pmol L −1 ) versus O 2 consumption rate (µmol L −1 d −1 ), and Ba residual is the deep ocean background value of Ba xs at zero oxygen consumption. The literature value of 180 pmol L −1 was used as the Ba residual value (Jacquet et al, 2008a, b;2011Planchon et al., 2013) in our calculations. The integrated mesopelagic Ba xs stock (µmol m −2 ) over the mesopelagic layer (100-1000 m) was calculated from the depth-weighted average Ba xs in order to investigate the link between the accumulated mesopelagic signal and the corresponding integrated remotely sensed primary productivity (PP).

Integrated remotely sensed PP
The integrated remotely sensed PP (mg C m −2 d −1 ) within the surface mixed layer was calculated using the CbPM algorithm (Behrenfeld et al., 2005), which requires chlorophyll concentration (mg m −3 ), particulate backscatter (λ 443 nm, m −1 ), photosynthetically active radiation (PAR; µmol photons m −2 d −1 ) and the MLD (m). Ocean Colour Climate Change Initiative (OC-CCI) data (https: //esa-oceancolour-cci.org/, last access: 14 August 2020), which blend existing data streams into a coherent record, meeting the quality requirements for climate assessment (Sathyendranath et al., 2019), were used for chlorophyll and particulate backscatter. PAR was taken from GlobColour (http://www.globcolour.info/, last access: 14 August 2020), and the MLD was taken from the climatology of de Boyer Montegut et al. (2004). The integrated remotely sensed PP data were regridded to 0.25 • spatially, using bilinear interpolation using the Python programming package xESMF (Zhuang, 2018), and averaged monthly. The area-averaged PP was averaged over a 6 × 1 • rectangular sample area, starting 6 • upstream extending up to the sampled station longitudinally, and 1 • latitudinally centred around each sampled station (see "Discussion" for details). In order to assess the validity of the remotely sensed PP data and demonstrate no meridional bias across the SO, the percentage of valid pixels was calculated for data north (90 ± 20 %; mean ± SD, n = 370) and south (82 ± 29 %; mean ± SD, n = 488) of the PF, revealing no bias.

Integrated percentage of POC re-mineralised
The integrated re-mineralised POC (mg C m −2 ) was estimated by multiplying the POC remineralisation flux (mg C m −2 d −1 ), as estimated using the Ba xs proxy method, by the number of days over which the corresponding remotely sensed PP (mg C m −2 d −1 ) was subsampled. The percentage of POC re-mineralised was then estimated as the percentage of integrated remotely sensed PP (mg C m −2 ) remineralised, assuming that the mesopelagic Ba xs stock signal observed is due to the remineralisation of the integrated surface PP signal.

Statistical analysis
For statistical analysis, the least squares regression method was applied for the assessment of significant correlations (Barbur et al., 1994). Significant differences between regions and regressions were tested using Welch's t test with an alpha of 0.05 (95 % confidence level) (Kokoska and Zwillinger, 2000).
As can be observed on the T -S plot of stations sampled (

Dissolved O 2
The water column dissolved O 2 concentrations ranged from 159 to 333 µmol L −1 (Fig. 2). Maximum concentrations were observed in the surface mixed layer, increasing southwards along the transect, with a mean value of 287 ± 40 µmol L −1 (mean ± SD, n = 700). A decrease in concentrations below the MLD coincided with an increase in σ θ . South of the PF, the decrease in dissolved O 2 concentrations at the MLD was sharp and relatively shallow when compared to profiles NPF, which were more gradual, spanning a wider depth range. Within the mesopelagic zone concentrations decreased down to 204 ± 29 µmol L −1 (mean ± SD, n = 6373) and then remained relatively uniform below 1000 m at 192 ± 113 µmol L −1 (mean ± SD, n = 12 950). Figure 2. Ba xs (black circles) with error bars, potential density (σ θ ; red) and dissolved O 2 (blue) profiles sampled along the transect, plotted against depth, for stations TM1 to TM8 from south to north. The shaded green area is the mesopelagic zone, and the hatched area is the ocean floor.

Ba xs and estimated POC remineralisation fluxes
Along the transect, Ba xs concentrations ranged from 59 to 684 pmol L −1 . All profiles exhibited a depletion of Ba xs in the upper surface waters (59-152 pmol L −1 ) and then a rapid increase below the MLD (∼ 150 m), with concentrations ranging between 113 and 684 pmol L −1 in the mesopelagic zone (100-1000 m, Fig. 2). At the two southernmost stations (TM1 and TM2), mesopelagic Ba xs peaks spanned a narrower depth range (100-600 m) than stations further north, with concentrations reaching values of ∼ 400 pmol L −1 . Concentrations were higher in the PFZ and SAZ with a maximum of 684 pmol L −1 in the PFZ at 48 • S (TM5). The subsurface increase in Ba xs started at slightly deeper depths (150-200 m) and spanned wider depth ranges down to 1000 m at stations north of the PF. The STZ station, at 41 • S (TM8), had the lowest concentrations, only increasing up to ∼ 200 pmol L −1 . Double peaks were observed at all stations north of the PF, with a shallow and more substantial peak occurring in the upper mesopelagic zone and a second peak in the lower mesopelagic zone. Below the mesopelagic zone, Ba xs concentrations decreased down to ∼ 180 pmol L −1 and remained relatively uniform.
The mean Ba residual concentration south of the PF was 183 ± 29 pmol L −1 (mean ± SD, n = 7), whereas it was 142 ± 45 pmol L −1 (mean ± SD, n = 8) between the PF and the STF. The two regions were not, however, significantly different to each other when conducting Welch's t test (t statistic = 2.10; p value = 0.06), and when averaging all concentrations below 2000 m along the transect, the Ba residual concentration was 161 ± 43 pmol L −1 (mean ± SD, n = 15). This concentration is not statistically different from the literature value of 180 pmol L −1 (Jacquet et al, 2008a, b;2011;Planchon et al., 2013), which is widely used for estimates of POC remineralisation fluxes. For a better comparison with these previous estimates, we used 180 pmol L −1 for the Ba residual concentration in our calculations.
The estimated POC remineralisation fluxes for the study area ranged from 6 to 96 mg C m −2 d −1 (Table S3), increasing northwards from the southernmost station up to the PFZ from 32 to 92 mg C m −2 d −1 and then decreasing down to 70 mg C m −2 d −1 at the SAF. The highest flux was estimated in the SAZ, and the lowest flux was estimated in the STZ.

Early wintertime Ba xs and Ba residual concentrations
A noticeable difference between profiles sampled early in the bloom season (Dehairs et al., 1997;Jacquet et al., 2015) versus those sampled later (Cardinal et al., 2001;Planchon et al., 2013) is the contrasted Ba xs concentrations in the surface mixed layer. Dehairs et al. (1997) has shown that these concentrations of Ba xs can be as high as 9000 pmol L −1 in areas of high productivity during spring, which then become depleted to concentrations below the SO Ba residual value of ∼ 180 pmol L −1 as productivity declines and surface POC export increases (Planchon et al., 2013). These high surface concentrations are not, however, due to the same process as the one that controls the Ba xs concentrations within the mesopelagic zone (Jacquet et al., 2011). Surface water concentrations are associated with Ba adsorbed onto particles, whereas the mesopelagic Ba xs signal is due to barite crystals formed within decaying bio-aggregates Lam and Bishop, 2007;Lemaitre et al., 2018;Sternberg et al., 2005). In this study, we observed surface depletion of Ba xs at all stations, in line with the assumption that the bulk surface export from the preceding bloom had been achieved at the time of sampling and that the majority of the Ba xs had been transferred to the mesopelagic zone.
A sharp increase in σ θ observed at the MLD has previously been identified as the depth at which decaying bioaggregates are formed (Lam and Bishop, 2007). These increases coincided with an increase in Ba xs (Fig. 2), linking the subsurface Ba xs signal to decaying bio-aggregates as per previous studies Dehairs et al., 1997;Jacquet et al., 2011). Additionally, decreases observed in dis-solved O 2 profiles along the transect were also accompanied by coinciding increases in Ba xs , in line with O 2 consumption due to remineralisation within the mesopelagic zone ( Fig. 2) Jacquet et al., 2005Jacquet et al., , 2011. The observed range of mesopelagic Ba xs concentrations (113-684 pmol L −1 ) were comparable to those previously reported in SO open waters (∼ 200-1000 pmol L −1 ; Cardinal et al., 2001Jacquet et al., 2005Jacquet et al., , 2008aJacquet et al., , b, 2011Jacquet et al., , 2015Planchon et al., 2013).
Ba xs profiles exhibited similar distributions to those reported throughout bloom seasons in the SO, with distinct peaks observed within the mesopelagic zone at all stations. Earlier in the bloom season, peaks mostly occur within the upper half of the mesopelagic zone (100-500 m; Cardinal et al., 2001Jacquet et al., 2005Jacquet et al., , 2008aJacquet et al., , 2011Jacquet et al., , 2015, but as the season progresses, they deepen down towards the bottom half of the mesopelagic zone (500− > 1000 m; Jacquet et al., 2008b, Planchon et al., 2013. Deepening and widening of the remineralisation depth range can be expected as the season progresses due to continued remineralisation taking place as particles sink to the bottom of the mesopelagic zone Planchon et al., 2013). This is also what we observed during early winter at stations NPF, with a second peak in deeper waters, as observed by Jacquet et al. (2008b) during the iron (Fe) fertilisation experiment (EIFEX). The deeper peak could also be linked to relatively larger cells that sink faster as they re-mineralise, possibly resulting in a large bloom earlier in the season.
A distinct latitudinal trend in mesopelagic depth-weighted average Ba xs has generally been observed in the SO with the highest values in the PFZ, decreasing north and southwards from the PF. These latitudinal trends tend to be accompanied by a coinciding trend in in situ surface biomass measurements Dehairs et al., 1997, Jacquet et al., 2011Planchon et al., 2013). During our early winter study, we observed a similar latitudinal trend in mesopelagic Ba xs stock (µmol m −2 ), with an increase from the southernmost station up to the PF, then varying around a maximum in the SAZ, down to the lowest value in the STZ, whereas temporally integrated remotely sensed PP increased progressively northwards to a maximum in the STZ (Fig. S1 in the Supplement). Time of sampling and extended blooms, which are characteristic of the SAZ (Thomalla et al., 2011), could be contributing factors to the higher values observed in PP and mesopelagic Ba xs distributions at stations north of the PF (Fig. S1). Contrary to what was expected, the profiles observed during our early winter study still show a significant mesopelagic remineralisation signal well after the summer bloom termination, which occurred between April and May (Fig. 3), as defined by the point in time when community losses outweigh the growth rate (Thomalla et al., 2011).
In deeper waters (below 2000 m) along the transect, south of the STF, where remineralisation is minimal compared to the mesopelagic zone, our Ba xs concentrations of 161 ± 43 pmol L −1 (mean ± SD, n = 15) are not signif- icantly different from the widely used Ba residual concentration of 180 pmol L −1 , measured during early spring to late summer (e.g. Jacquet et al., 2008a, b;2011;Planchon et al., 2013). We thus did not observe a wintertime decline to an expected "true" SO background value, when PP and bacterial activity are suspected to be minimal (Jacquet et al., 2011). There are two possible explanations for this. Firstly, the decline to a winter background signal might never be achieved due to ongoing barite precipitation and remineralisation, as well as the release of labile Ba attached to phytoplankton as they decay, precipitating into barite crystals, which could possibly continue throughout winter . Secondly, the low sinking speed of suspended barite (∼ 0.3 m d −1 ; Sternberg et al., 2008), once produced in the mesopelagic layer, implies that it would take ∼ 6 years (not considering reaggregation and redissolution) to sink from 300 m (∼ peak of production) to the bottom of the mesopelagic layer (1000 m). The "true" background value may thus have to be measured at the very end of winter just before the initiation of the spring bloom. This also suggests that the Ba xs signal in the mesopelagic layer may represent remineralisation activity over more than a few days to weeks, per previous reports (e.g. Dehairs et al., 1997;Jacquet et al., 2015;Planchon et al., 2013).

Timescale of the mesopelagic Ba xs signal
The Ba xs signal that we observed in winter is in agreement with the suggestion by Dehairs et al. (1997) that there can be significant carry-over between bloom seasons. Other studies have also pointed out that the timescale of this proxy is longer than a snapshot view  and have highlighted a seasonal increase in mesopelagic Ba xs (Jacquet et al., 2011). This strongly suggests that the Ba xs signal is not directly linked to synoptic measurements of PP at the time of sampling. In order to investigate this hypothesis, for the first time, we compiled a SO mesopelagic Ba xs stock dataset with all available literature data including data from this study (Fig. 4a, Table S3). The mesopelagic Ba xs stock was integrated over the Ba xs peak depth range (as identified in each study). As can be seen on the map of the compilation dataset (Fig. 4a), these data points were collected across the three basins of the SO over ∼ 20 years. Despite this diversity in observations, a statistically significant accumulation of mesopelagic Ba xs with time is still observed SPF (Fig. 4b) and NPF (Fig. 4c). Mesopelagic Ba xs accumulates at a rate of 0.86 (±0.15) µmol m −2 d −1 SPF (R 2 = 0.43, p value 0.05, n = 43; Fig. 4b) and at 0.88 (±0.20) µmol m −2 d −1 NPF (R 2 = 0.41, p value < 0.05, n = 31; Fig. 4c), with no statistically significant difference between the two regions (Welch's t test = 0.24; p value = 0.80).
A possible link between the integrated mesopelagic Ba xs stock and the corresponding integrated remotely sensed PP was assessed for all studies conducted after September 1997, when remotely sensed PP data became available. To do so, we first estimated that sub-millimetre-sized aggregates would take ∼ 20 d to sink down to 1000 m (considered as the bottom of the mesopelagic zone in this study) using a sinking speed of 50 m d −1 that corresponds to an average literature value (50-100 m d −1 : Riebesell et al., 1991;50- (Ferrari and Nikurashin, 2010), it was estimated that these aggregates would have originated 346 km west from the station that was sampled for mesopelagic Ba xs ∼ 20 d prior. Using this distance, the dimensions of the sample area were set with the southernmost station (TM1) of this study, where degrees of longitude cover the smallest area. For the sake of consistency this sample area was applied to all sampling locations of the considered dataset. The integrated remotely sensed PP (see Sect. 2.5) was then averaged spatially, starting 6 • upstream extending up to the sampled station longitudinally, and 1 • latitudinally centred around each station, in order to capture the surface PP that is assumed to translate to the mesopelagic remineralisation and measured Ba xs stock.
The monthly averaged remotely sensed PP, at the time of sampling, was compiled for the considered dataset, and we found that the PP over the growing season (Fig. 4d, e) reaches the highest values between January and February (day 125 to 175 of the year), thereafter steadily decreasing to minimal values in July (∼ day 310 of the year, i.e. during our study). The mesopelagic Ba xs accumulation over time can not, therefore, be matched with the remotely sensed PP measured during the month of sampling. A possible relationship between mesopelagic Ba xs stock and temporally integrated remotely sensed PP was further investigated by considering longer timescales. Remotely sensed PP of the preceding bloom was temporally integrated from the preceding September prior to sampling as the start of the bloom, in general agreement with previous bloom phenology studies (Thomalla et al., 2011), and up to 1 month prior to the sampling date of the study, taking into consideration time needed for export, aggregate formation and barite crystal release through remineralisation (∼ 1 month).
Varying timescales were considered between the preceding September up to 1 month prior to sampling (Sept -T1,  Table S4), in monthly increments, that could influence the relationship between remotely sensed PP and the mesopelagic Ba xs stock (Table S4). The strongest and most significant correlation between the mesopelagic Ba xs stock and integrated remotely sensed PP, for both north and south of the PF, was obtained from the preceding September up to 1 month prior to sampling (Sept -T1, Table S4; SPF: Fig. 5a, R 2 = 0.55, p value < 0.05, n = 39; NPF: Fig. 5b, R 2 = 0.42, p value < 0.05, n = 31). When remote sensing data were limited due to cloud cover and low sunlight during winter months, specifically at the southernmost stations, all available data were used for the duration of the season. The correlation observed in the STZ is not significant at a 95 % confidence level (p value = 0.10); however, the limited number of data points (n = 6) may preclude any significance from emerging. The significant positive correlations obtained south of the STF suggest that mesopelagic Ba xs stock can be used as a remineralisation proxy on an annual timescale instead of only a few weeks. Figure 5 also reveals that for a given PP the mesopelagic Ba xs stock was 2-fold higher SPF compared to NPF (Welch's t test, t statistic = 2.24; p value < 0.05); this is further discussed below.

Environmental factors influencing mesopelagic remineralisation and carbon export efficiency
Estimated POC remineralisation fluxes along the transect (6-96 mg C m −2 d −1 ) were on the upper end of the range of fluxes reported in previous studies, with the exception of the STZ station, but within the same order of magnitude for the SO as estimated from spring to autumn (0.2-118 mg C m −2 d −1 ; Table S3; Jacquet et al., 2011Jacquet et al., , 2015Planchon et al., 2013). As the bloom season progresses, more efficient remineralisation rates have been reported in multiple studies Jacquet et al., 2011;Planchon et al., 2013). However, during late summer as the bloom declines, observations indicate an inefficient BCP due to enhanced surface nutrient recycling (Dehairs et al., 1992;Jacquet et al., 2011;Planchon et al., 2013), leading to a decrease in surface POC export (Planchon et al., 2013). Seasonal variation is reported to be more pronounced northwards within the SO with the least variation observed in the southern Antarctic circumpolar current (Dehairs et al., 1997;Planchon et al., 2013). The percentage of mesopelagic POC remineralisation as calculated from estimated POC remineralisation fluxes over integrated remotely sensed PP, for the SO compilation dataset (SPF: 19 ± 15 %, n = 39; NPF: 10 ± 10 %, n = 29; mean ± SD; t statistic = 2.75; p value < 0.05; Table S3), was ∼ 2-fold higher SPF than NPF, revealing the higher surface carbon export efficiency SPF. Our estimates of percentage of POC re-mineralised fall within the range of reported export efficiencies throughout the SO (2 %-58 %; Jacquet et al., 2011;Morris et al., 2007;Savoye et al., 2008). Our values also support the inverse relationship between export efficiency and productivity, with higher export efficiency in areas of lower production and high productivity low e-ratio (HPLE), where e-ratio refers to the ratio between export production and net primary productivity (Fan et al., 2020;Maiti et al., 2013). Estimated mesopelagic POC remineralisation has been reported to account for a significant fraction of exported carbon in the PFZ and southwards, from 31 % to 97 % from spring to summer, whereas it only accounts for ∼ 50 % in the SAZ and SAF during summer . A combination of variables can influence surface export efficiency and the magnitude of the subsequent mesopelagic remineralisation, even more so when considering longer timescales. These variables include physical dynamics and interlinked biogeochemical factors, i.e. bacterial activity, phytoplankton community structure, zooplankton grazing and nutrient availability (Bopp et al., 2013;Buesseler and Boyd, 2009;Jacquet et al., 2008b;Pyle et al., 2018). In previous studies, supply and loss via physical transport were deemed negligible relative to decay and loss via production due to minimal advection and diffusion gradients observed on the timescale of days to weeks. These processes were therefore assumed to have minimal impact on the mesopelagic signal (Dehairs et al., 1997;Figure 4. (a) Positions of Ba xs observations compiled from all known SO studies on a cylindrical equal-area projection of the SO; the three SO basin cut-offs are indicated by the dashed black lines: from left to right Pacific, Atlantic and Indian. Integrated mesopelagic Ba xs stock plotted against day of year sampled, with the 1st of September set as day 1, for all available literature data and winter data from this study. Data were split into two zones using the polar front (PF) to divide the SO (b) south of the PF (SPF) and (c) north of the PF (NPF). Monthly averaged remotely sensed PP plotted against day of year for locations and dates of the SO compilation dataset and winter data from this study: (d) SPF and (e) NPF. Open circles are data points from studies which did not use HF in the particulate sample digestion procedure; regressions did not include these data. There was, however, no significant difference when including these data points (Table S3). Planchon et al., 2013;Rutgers van der Loeff et al., 2011). It has, however, been observed that features such as mesoscale eddies can have an effect on Ba xs distribution by influencing particle patterns on a broad spatial scale, homogenising mesopelagic remineralisation signals by causing relatively flat profiles or shallower remineralisation peaks Jacquet et al., 2008b). The region of our winter study is known for being a mesoscale eddy hotspot due to the South-west Indian Ridge (Ansorge et al., 2015). In the STZ, extremely dynamic submesoscale activity due to the Agulhas return current may indeed have significantly influenced the mesopelagic signal and may help explain the absence of correlation with integrated surface PP. On the contrary, south of the STF, the significant correlations seem to indicate that physical transport variability is not the main process affecting the mesopelagic Ba xs signal and that biogeochemical factors may be dominant.
The Fe-limited SAZ  and AZ (Viljoen et al., 2018) have generally mixed and seasonally changing assemblages of pico-, nano-and microphytoplankton (Eriksen et al., 2018;Gall et al., 2001). Diatoms tend to dominate in the silicate-rich waters south of the PF (Petrou et al., 2016;Rembauville et al., 2017;Wright et al., 2010), whilst seasonally silicate-limited waters north of the PF favour smaller phytoplankton groups (Freeman et al., 2018;Nissen et al., 2018;Trull et al., 2018). HPLE regimes are indeed characteristic of large areas of the SAZ mainly due to surface POC accumulation caused by non-sinking particles, tending towards less efficient export of smaller cells (Fan et al., 2020). Even when large particles are abundant in HPLE surface layers, a complex grazing community may prevent the export of large particles (Dehairs et al., 1992;Lam and Bishop, 2007). This can explain the higher surface carbon export efficiency that we estimate SPF compared to NPF. Export efficiency has also been linked to bacterial Figure 5. Integrated mesopelagic Ba xs stock plotted against integrated remotely sensed PP from the preceding September up to 1 month prior to sampling and all available literature data and winter data from this study (a) south of the PF (SPF, black squares) and (b) north of the PF (NPF, black circles). Red open squares are data points from our winter dataset for which there were not sufficient remote sensing PP data to integrate up to 1 month prior to sampling, and available data up to 3 months prior to sampling were plotted but not included in the statistical analysis. productivity with efficient surface remineralisation limiting surface POC export when most of the water column integrated bacterial productivity is restricted to the upper mixed layer (Dehairs et al., 1992;Jacquet et al., 2011), which can be the case to varying degrees throughout the SO. In the STZ phytoplankton communities are reported to be dominated by prokaryotic picoplankton including cyanobacteria and prochlorophytes (Mendes et al., 2015). These groups utilise regenerated nutrients in the surface mixed layer tending towards diminished surface export efficiency with high concentrations of non-sinking POC (Fan et al., 2020;Planchon et al. 2013). In addition to this, the potential influence of high submesoscale activity may explain the low mesopelagic Ba xs measured at the STZ station of this study despite it being the station with the highest integrated PP (Fig. S1). Linking temporally integrated remotely sensed PP to mesopelagic Ba xs stock, coupled with the added influence of physical dynamics affecting surface export efficiencies and along longer timescales, could give better estimates of export and remineralisation signals throughout the SO on an annual and basin scale. Our estimates of percentage re-mineralised POC over remotely sensed PP may contribute to the improved modelling of the C cycle over the SO on an annual timescale.

Conclusions
Our unique early winter Ba xs data were similar in magnitude and exhibited the same relationship with σ θ and dissolved O 2 gradients as observed in summer, indicating that processes controlling this signal in summer are still driving the signal in early winter. The expected decline of the mesopelagic Ba xs signal to background values during winter was not observed in this study, supporting the hypothesis that this remineralisation proxy likely has a longer timescale than previously reported. The absolute decline might be delayed due to the cumulative behaviour of mesopelagic Ba xs , ongoing remineralisation and barite precipitation. The "true" SO background value may thus have to be measured at the very end of winter, prior to bloom initiation.
Significant positive correlations north and south of the PF between mesopelagic Ba xs stock and remotely sensed PP integrated from September to 1 month before sampling (Sept -T1, Table S4), in combination with significant Ba xs accumulation trends obtained for the SO compilation dataset, suggest an annual timescale. They may also indicate that physical processes do not dominate the mesopelagic signal on an annual scale within the SO and that biogeochemical factors are dominant. There is no significant difference in mesopelagic Ba xs and POC remineralisation north and south of the PF, but the significantly higher integrated remotely sensed PP to the north when compared to the south indicates a greater export efficiency south of the PF. This is in accordance with the phenomenon of HPLE regimes which are common throughout the SO, more so north of the PF than south of the PF (Fan et al., 2020). The longer timescale of Ba xs and the cumulative behaviour of this proxy in the mesopelagic zone make it possible to use Ba xs on an annual scale for the estimation of POC remineralisation fluxes throughout the SO and to better understand how variable environmental factors influence these processes on a basin scale. We believe that the significance of these relationships will improve as more data 3220 N. R. van Horsten et al.: Early winter barium excess in the southern Indian Ocean become available (e.g. GEOTRACES IDP2021), which will assist in better understanding and constraining the timescale of remineralisation and carbon export efficiency throughout the SO.
Data availability. Data used in this study have been published in the online open-source repository Zenodo and can be accessed at https://doi.org/10.5281/zenodo.6583338 (van Horsten et al., 2022).
Author contributions. This study was conceptualised by NRvH, HP, GS and EB. Formal analysis, investigation and validation of data was carried out by NRvH, HP, GS, EB and TJRK. NRvH and TJRK contributed towards the visualisation of the data. HP, GS, TNM, AR, NL and EB contributed towards supervision and resources. Funding was acquired by NRvH, TNM, AR and EB. All authors contributed towards writing, reviewing and editing the final manuscript.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements. This work was supported by the ISblue project, Interdisciplinary Graduate School for the Blue Planet (ANR-17-EURE-0015), and co-funded by a grant from the French government under the programme "Investissements d'Avenir". International collaboration was made possible by funding received by the French-South African National Research Foundation (NRF) Collaboration (PROTEA; FSTR180418322331), NRF funding (SNA170518231343 and UID 110715) including funding from South African Department of Science and Innovation, French Ministry of National Education, Higher Education and Research, and the French Ministry of Foreign Affairs and International Development. We would like to thank the captain and crew of the R/V SA Agulhas II for their invaluable efforts, as well as all the research participants who assisted our fieldwork. Thanks to Isabelle Ansorge, Marcel du Plessis and Esther Portela for their assistance with water mass identification and Catherine Jeandel for her invaluable expert insight. Finally, we thank the three reviewers, Frank Dehairs, Stéphanie Jacquet and James K. B. Bishop, for their insightful reviews, thereby improving the quality of our manuscript. A special thank you to Frank Dehairs for sharing the Indigo 3, EPOS 2 and ANTX/6 data.
Financial support. This research has been supported by the Council for Scientific and Industrial Research, South Africa (grant no. SNA170518231343).
Review statement. This paper was edited by Carol Robinson and reviewed by Frank Dehairs, Stéphanie Jacquet, and J. K. B. Bishop.