Particulate biogenic barium tracer of mesopelagic carbon remineralization in the Mediterranean Sea (PEACETIME project)

. We report on the sub-basin variability in particulate organic carbon (POC) remineralization in the western and central Mediterranean Sea in late spring during the PEACETIME (ProcEss studies at the Air–sEa Interface after dust deposition in the MEditerranean sea) cruise. POC remineralization rates were estimated using the excess biogenic particulate barium (Ba xs ) inventories in the mesopelagic layers (100–1000 m depth) and compared with prokaryotic heterotrophic production (PHP). Ba xs -based mesopelagic remineralization rates (MRs) ranged from 25 ± 2 to 306 ± 70 mgCm − 2 d − 1 . MRs were larger in the Algero-Provençal (ALG) Basin than in the Tyrrhenian (TYR) and Ionian (ION) basins. Our Ba xs inventories and integrated PHP data also indicated that signiﬁcant mesopelagic remineralization occurred down to 1000 m depth in the ALG Basin in contrast to the ION and TYR basins, where remineralization was mainly located above 500 m depth. We propose that the higher and deeper MRs in the ALG Basin were sustained by an additional particle export event driven by deep convection. The TYR Basin (in contrast to the ALG and ION basins) presented the impact of a previous dust event, as reﬂected by our particulate Al water column concentrations. The ION and TYR basins showed small-scale heterogeneity in remineralization processes, reﬂected by our Ba xs inventories and integrated PHP data at the Tyrr long-duration station. This heterogeneity was linked to the mosaic of blooming and non-blooming patches reported in this area during the cruise. In contrast to the western Mediterranean Sea (ALG Basin), the central Mediterranean Sea (ION and TYR basins) showed lower remineralization rates restricted to the upper mesopelagic layer during the late spring PEACETIME cruise.


Introduction
In the ocean, the remineralization rate associated with sinking particles is a crucial variable for the air-sea CO 2 balance (Kwon et al., 2009). Most of the sinking particulate organic carbon (POC) conversion (i.e. remineralization) into CO 2 by heterotrophic organisms (i.e. respiration) occurs within the mesopelagic zone (100-1000 m) (Martin et al., 1987;Buesseler et al., 2007;Buesseler and Boyd, 2009). Thus, a quantitative representation of this process is crucial for future predictions of the ocean's role in the global C cycle (IPCC, 2014). Particulate biogenic barium (Ba xs ) is a geochemical tracer of POC remineralization in the mesopelagic layer. Ba xs occurs in the form of barite (BaSO 4 crystals) in the dark ocean as a by-product of prokaryotic remineralization. In a global ocean undersaturated with respect to barite (Monnin and Cividini, 2006), Ba xs precipitates inside oversaturated biogenic micro-environments during POC degradation by heterotrophic prokaryotes, through sulfate and/or barium enrichment (Dehairs et al., 1980;Stroobant et al., 1991;Bertram and Cowen, 1997;Ganeshram et al., 2003). By applying a transfer function relating Ba xs to O 2 consumption (Dehairs et al., 1997) Ba xs has been widely used since the 1990s as an estimator of mesopelagic POC remineralization rates in various sectors of the Southern Ocean, North Pacific and North Atlantic (Cardinal et al., 2001(Cardinal et al., , 2005Dehairs et al., 2008;Jacquet et al., 2008aJacquet et al., , b, 2011Jacquet et al., , 2015Planchon et al., 2013;Lemaitre et al., 2018). Jacquet et al. (2021) recently reported that such a transfer function could be applied in the Mediterranean Sea without restriction. The above-mentioned study complemented previous investigations that aimed at improving the use of Ba xs to estimate local POC remineralization processes in the Mediterranean Sea (Jacquet et al., 2016;Jullion et al., 2017). The Mediterranean Sea represents a unique case study, mainly due to unresolved issues related to the imbalance in the regional C budget, such as the coupling between surface biology and deeper remineralization, the timescales of their variability between basins, and discrepancies between mesopelagic trophic structure and respiration dynamics (Sternberg et al., 2007(Sternberg et al., , 2008Santinelli et al., 2010;Lopez-Sandoval et al., 2011;Luna et al., 2012;Tanhua et al., 2013b;Malanotte-Rissoli et al., 2014).
The present work is part of the PEACETIME (ProcEss studies at the Air-sEa Interface after dust deposition in the MEditerranean sea) project (http://peacetime-project.org/, last access: 7 April 2021), which aims to study the impact of atmospheric Saharan dust on Mediterranean biogeochemistry (Guieu et al., 2020a). Dust deposition is a major source of macro-and micronutrients and ballasted material to surface waters that likely impacts the biological carbon pump through organic matter production (i.e. primary production) and its subsequent export and remineralization in the water column (Pabortsava et al., 2017;Gazeau et al., 2021). Overall, the aims of the present contribution to the PEACETIME project were (1) to document particulate biogenic Ba xs in different ecoregions of the western and central parts of the Mediterranean Sea, as previous Ba xs data in the Mediterranean Sea are relatively scarce with limited vertical sampling resolution (Sanchez-Vidal et al., 2005) or restricted locations (Dehairs et al., 1987;Sternberg et al., 2007Sternberg et al., , 2008van Beek et al., 2009); (2) to determine the relationship between Ba xs and environmental variability, including dust deposition; (3) to estimate Ba xs -based POC remineralization rates (MRs) at mesopelagic depths using Dehairs' transfer function (Dehairs et al., 1997), which we recently validated for the Mediterranean Sea ; and (4) to assess potential differences in the remineralization length scale of POC in the various ecoregions of the Mediterranean Sea.

Study area
The PEACETIME cruise (https://doi.org/10.17600/ 17000300) was conducted during late spring, from 10 May to 11 June 2017, in the western and central Mediterranean using the French research vessel (RV) Pourquoi pas? (Fig. 1a). The Mediterranean Sea is a semi-landlocked sea with limited but crucial exchange with the Atlantic Ocean, two deep overturning cells, one shallow circulation, and a complex upper-layer circulation with several permanent and quasi-permanent eddies.
The hydrography during the PEACETIME cruise was characterized by three layers: surface, intermediate and deep waters, typical for the Mediterranean (Tamburini et al., 2013;Tanhua et al., 2013a;Hainbucher et al., 2014, Malanotte-Rizzoli et al., 2014. Briefly, the main water masses are (see potential temperature-salinity diagram in Fig. 1b) as follows: (1) from west to east, Atlantic Surface Water (SW) is gradually replaced by Ionian Surface Water (ISW) and Levantine Surface Water (LSW); (2) Winter Intermediate Water (WIW); (3) Levantine Intermediate Water (LIW), which is present at intermediate depths (from 200 to 800 m) and is characterized by a local maximum of salinity and a local minimum of dissolved oxygen concentration; and (4) Mediterranean Deep Water (MDW).
Three main ecoregions (Reygondeau et al., 2017;Ayata et al., 2018) were crossed during the cruise: the Algero-Provençal Basin (later referred to as ALG), the Tyrrhenian Basin (TYR) and the Ionian Basin (ION) (Fig. 1a). These basins displayed the typical eastward oligotrophic gradient reported in previous studies (Moutin and Raimbault, 2002;Durrieu de Madron et al., 2011;Pujo-Pay et al., 2011;Tanhua et al., 2013a;Reygondeau et al., 2017;Guieu et al., 2020a). However, this trend was not homogeneous -for instance, in the Ionian Sea (a crossroad for waters with contrasting biological histories) where a mosaic of blooming and nonblooming areas co-occurred in spring (Berline et al., 2021). A diatom-dominated deep chlorophyll maximum that coincided with a maximum in biomass and primary production (PP) was well developed and observed along the entire cruise track (Marañón et al., 2021). PP is described in detail in Van Wambeke et al. (2020). Furthermore, important dust deposition affected the TYR Basin a few days before our arrival at station nos. Tyrr and 5, whereas in the ALG Basin, dust deposition occurred few hours before our sampling at Fast station (Bressac et al., 2021). POC downward fluxes measured at 200 m depth were similar at the three long-duration stations (Fast, Tyrr and Ion).

Barium sampling and sample processing
A total of 13 stations were sampled for particulate barium from the surface to 2000 m (30 depths in total) in the ALG (station nos. 1, 2, 3, 10, Fast, 9 and 4), TYR (station nos. 5, Tyrr and 6) and ION (station nos. 8, 7 and Ion) basins ( Table 1). Three of these stations were sampled twice on different days (long-duration stations Fast, Tyrr and Ion); however, due to technical problems, no particulate barium data are available for the second visit to Ion station. A total of 3 d separate both visits to Fast station and 2 d separate both visits to Tyrr station.
For particulate barium, 4-6 L of seawater sampled using Niskin bottles was filtered onto 47 mm polycarbonate membranes (0.4 µm porosity) under slight overpressure sup-  Fig. 2. (b) Potential temperature-salinity diagram with isopynals (kg m −3 ) for PEACETIME profiles. The figure was produced using Ocean Data View (Schlitzer, 2002). plied by filtered air. Filters were rinsed with a few millilitres of Milli-Q water to remove sea salt, dried (50 • C) and stored in Petri dishes for later analysis. In the laboratory, we performed a total digestion of filters using a concentrated tri-acid (0.5 mL HF/1.5 mL HNO 3 /HCl 1 mL; all Optima grade) mixture in closed Teflon beakers overnight at 95 • C in a clean pressurized room. After evaporation close to dryness, samples were re-dissolved into 10 mL of 2 % HNO 3 . Subsequently, samples were analysed for Ba and other elements of interest (i.e. Al, Na, Sr and Ca) using high-resolution inductively coupled plasma-mass spectrometry (HR-ICP-MS; Element XR, Thermo Scientific). Based on analyses of external certified reference standards, accuracy and reproducibility were both within ±5 %. More details on sample processing and analysis are given in Cardinal et al. (2001) and Jacquet et al. (2015). The presence of sea salt was checked by analysing Na, and the sea-salt particulate Ba contribution was found to be negligible (< 0.1 % of total Ba). Particulate biogenic barium in excess (hereafter referred to as Ba xs ) was calculated as the difference between total Ba and lithogenic Ba. The lithogenic Ba concentration was determined using the Al concentration and the upper-continentalcrust (UCC) Ba : Al molar ratio (Dymond et al., 1992;Taylor and Mc Lennan, 1985). The standard uncertainty (Ellison et al., 2000) on Ba xs concentration ranges from 5.0 % to 5.5 %. The term "in excess" is used to indicate that concentrations are larger than the Ba xs background (Ba Bkg). The background (or residual value) is considered as "preformed" Ba xs at zero oxygen consumption left over after the transfer and partial dissolution of Ba xs produced during the degradation of previous particle export events. This background Ba xs value likely depends on the saturation state of the water with respect to barite (BaSO 4 , the main phase of partic-ulate biogenic barium). Saturation indexes were reported in Jacquet et al. (2016) over a high-resolution and quasi-zonal Mediterranean transect (M84/3 cruise; Tanhua et al., 2013a, b). They revealed that the water column throughout the study area is largely undersaturated, with a saturation state ranging from 0.2 to 0.6. A background Ba xs value of 130 pM was recently reported in Jacquet et al. (2021). It is close to the average Ba xs content observed at greater depth (> 1000 m) in the present study (see below).

Prokaryotic heterotrophic production
Prokaryotic heterotrophic production (PHP) estimation was measured by the L-(4,5-3 H)Leucine ( 3 H-Leu, specific activity 109 Ci mmol −1 , PerkinElmer ® ) incorporation technique (Kirchman, 1993). Details of the protocols can be found in Van Wambeke et al. (2020). Briefly, in epipelagic layers (0-200 m), 1.5 mL seawater samples were incubated at 20 nM 3 H-Leu final concentration using the microcentrifuge technique (Smith and Azam, 1992). For the mesopelagic layers, 20 mL (200-800 m depth) and 40 mL (below 800 m depth) seawater samples were incubated using 20 nM and 10 nM 3 H-Leu final concentration respectively with the filtration technique (Tamburini et al., 2002). Samples were incubated at in situ temperature. To calculate PHP, we used the empirical conversion factor of 1.5 ng C per picomole of incorporated leucine assuming that isotopic dilution was negligible under saturating concentrations of leucine as checked occasionally from concentration kinetics (Van Wambeke et al., 2020).

POC remineralization rates
We recently reported on the validity of the Dehairs's transfer function (Dehairs et al., 1997) in the Mediterranean Basin to estimate mesopelagic POC remineralization . We applied a similar approach to estimate the remineralization rate (MR): where Ba xs is the depth-weighted average Ba xs concentration (DWA; pM), i.e. the Ba xs inventory divided by the depth layer considered (Z), and RR is the Redfield C/O 2 molar ratio (127/175; Broecker et al., 1985). As reported above, a Ba Bkg concentration of 130 pM was used. MRs were then integrated over the 100-500 m (upper mesopelagic zone) and 100-1000 m (entire mesopelagic zone) depth layers.

Results
Particulate biogenic Ba xs , biogenic Ba fraction (%), and particulate Al, Sr and Ca concentrations are reported in Fig. 2 in the upper 2000 m of the water column along a zonal transect crossing the three main sub-basins. PHP rates are also reported in Fig. 2 in the upper 500 m along the same transect. Ba xs displayed a similar water column distribution to that reported in various sectors of the global Ocean: relatively low surface concentrations; a maximum in the mesopelagic layer (100-1000 m); and a decrease in concentrations back to a background level at deeper depths, usually 1000 m (Fig. 3). At station nos. 9, Tyrr, 8 and Ion, Ba xs concentrations in the upper 100 m were quite high (> 5000 pM), with values reaching up to 11 700 pM at 80 m depth at Tyrr station (Fig. 2a). Such high Ba xs concentrations in the upper layers  (Schlitzer, 2002). are quite unusual, although similar values (up to 9000 pM) were occasionally observed in earlier Southern Ocean studies (Dehairs et al., , 1997Jacquet et al., 2007bJacquet et al., , 2008a. The high Ba xs values at Tyrr and Ion stations were associated with higher Sr (up to 4267 pM at 100 m at Ion station) and Ca (> 130 nM) concentrations (Fig. 2d, e). Throughout the water column, Sr and Ca concentrations ranged from to 448 to 6938 pM and from 30 to 488 nM respectively. Sr/Ca molar ratios ranged from 7 to 45 mmol mol −1 , which is within the range of ratios reported for organic material (Martin and Knauer, 1976). The upper mesopelagic layer (100-500 m) showed the characteristic Ba excess (maximum), as illustrated in Fig. 3a. The lithogenic impact on the Ba xs signal was relatively low (< 20 %) except at station nos. 4, 5 and Tyrr where it reached up to 30 % at some depths in the water column (Figs. 2b, 3b). High Ba xs concentrations at stations in the ALG Basin and at station no. 7 in the ION Basin extended deeper than at stations in the TYR Basin (Fig. 2a). At Ion station, the Ba xs maximum coincided with higher Ca (up to 186 nM) concentrations in the upper mesopelagic layer (Fig. 2e). However, the Ba xs maximum also extended deeper. This is especially salient at stations in the ALG Basin. At the other stations, Ba xs concentrations below 500 m decreased to reach the background value of around 130 pM. Among stations sampled twice for barium during the cruise, Fast station (ALG Basin) presented similar Ba xs profiles except between 400 and 1000 m depth where lower concentrations were measured during the second visit (3 d later; Fig. 3a). The Ba xs signal was mostly biogenic and rather stable over the whole water column at this station. This was also the case at Ion station. In contrast, at Tyrr station, differences between Ba xs profiles mainly occurred in the surface layer and upper mesopelagic layer, with relatively higher Ba xs peaks during the second visit (2 d later; Fig. 3b). The biogenic Ba fraction was also more variable throughout the water column at Tyrr station.
PHP rates decreased from west to east in surface waters (Fig. 2f). At Fast station, PHP rates decreased from 49 ng C L −1 h −1 in surface waters to values between 7 and 11 ng C L −1 h −1 at 100 m depth and below 6 ng C L −1 h −1 below 200 m depth (Fig. 3d). The same trends were found at Tyrr and Ion stations with values in surface waters of around 36 and 25 ng C L −1 h −1 respectively (Fig. 3e, f). Depth-weighted average (DWA) concentrations of Ba xs are reported in Table 1 and Fig. 5 for the upper (100-500 m) and entire (100-1000 m) mesopelagic layer. As the base of the mixed layer was shallower than 100 m, this depth is taken as the upper boundary of the mesopelagic domain. DWA values ranged from 221 to 979 pM. On average, stations located in the ALG Basin presented higher DWA values than in the TYR and ION basins. DWA Ba xs values remained rather stable over the 3 d period at Fast station between 100 and 500 m depth, but decreased in deeper layers (Fig. 5). As a consequence, the DWA changed from 527 to 381 pM for the entire 100-1000 depth layer. In contrast, at Tyrr station, DWA Ba xs values for the 100-500 m and 100-1000 m depth layers increased over the 2 d period (from 284 to 542 pM and from 200 to 380 pM respectively). On average DWA Ba xs reached 577 ± 286, 378 ± 123 and 529 ± 213 pM (100-500 m), and 527 ± 288, 280 ± 82 and 358 ± 112 pM (100-1000 m) in the ALG, TYR and ION basins respectively.

Ba xs distributions across the sub-basins
The very high Ba xs concentrations reported in the surface layer at station nos. 9, Tyrr, 8 and Ion were associated with local Sr and Ca maxima, likely linked to potentially ballasted phytoplankton-derived material. Similar observations have previously been reported in the Southern Ocean, revealing that particulate Ba xs is incorporated into or adsorbed onto biogenic material in the surface water, with barite being a minor component (Dehairs et al., , 1997Jacquet et al., 2007aJacquet et al., , 2008a. In deeper layers, Ba xs presented the characteristic maximum reflecting mesopelagic remineralization processes. Mesopelagic Ba xs distributions presented here were similar to those reported in Jacquet et al. (2021) and Sternberg et al. (2008) in the northwestern Mediterranean Sea (the ANTARES and DYFAMED observatory sites respectively). The Ba xs maximum extended down to 1000 m depth in the ALG Basin, whereas it was mostly located in the upper 500 m in the TYR Basin. The lithogenic impact on the Ba xs signal was relatively very low (< 5 %), except at station nos. 4, 5 and Tyrr where it was more variable and reached up to 30 % at some depths (Fig. 2b, c). A large dust deposition event occurred over a large area including the southern Tyrrhenian Sea just before the beginning of the PEACETIME campaign. Particulate Al concentrations and the estimated lithogenic Ba fraction were sampled at these stations 5-12 d after the event and reflected the impact of this dust event at depth. These conclusions are further supported by results reported in Bressac et al. (2021), showing that Saharan dust depositions strongly impacted station nos. 4, 5, Tyrr and 6, where a significant fraction of dust particles was transferred to mesopelagic depths.

Mesopelagic Ba xs and prokaryotic heterotrophic production
Previous studies have highlighted the relationship between the mesopelagic Ba xs and the vertical distribution of prokaryotic heterotrophic production (PHP), reflecting the temporal progression of the POC remineralization processes. In mesopelagic layers, the Ba xs content is smaller when most of the PHP occurs in the upper mixed layer (indicating an efficient, close to complete remineralization within the surface), compared with situations where a significant part of the PHP is located deeper in the water column (reflecting significant deep prokaryotic activity and POC export). Figure 3 shows the PHP profiles at the Fast, Tyrr and Ion long-duration stations (see also Van Wambeke et al., 2020, for more details on PHP). Figure 4 shows the ratio of integrated surface (100 m) to integrated upper mesopelagic (500 m) PHP vs. DWA Ba xs values calculated over the 100-500 m depth interval. Results are compared to the data obtained in the Southern Ocean (Jacquet et al., 2008a(Jacquet et al., , 2015 and recently in the northeast Atlantic and northwestern Mediterranean Sea (the PAP and ANTARES observatory sites respectively) . The blue line in Fig. 4 represents the trend obtained during KEOPS2 (Kerguelen Ocean and Plateau Study 2; Jacquet et al., 2015); it does not include the encircled data points referred to as "season advancement". Results during PEACETIME followed a similar trend to that found for KEOPS2, with higher DWA Ba xs values in situations where a significant part of the column-integrated PHP is located deeper in the water column (high IntPHP100/IntPHP500 ratio, Fig. 4). Note that some data points, characterized by low DWA Ba xs values, did not follow the trend from KEOPS2 (station nos. 3, 5 and Tyrr2). During KEOPS2, the lowest DWA values were reported for stations located in a meander and reflecting different (earlier) stages of a bloom compared with the other stations (see "season advancement" in Fig. 4). Similarly, station nos. 5 and Tyrr2 reflected the temporal evolution of the establishment (or advanced stages) of mesopelagic remineralization processes in the TYR Basin compared with the other basins. Measurements carried out during the second visit to Tyrr station 2 d later corroborated this hypothesis, showing an increase in remineraliza-tion in the upper mesopelagic layer (DWA Ba xs increased from 284 to 542 pM). At the DYFAMED station, Sternberg et al. (2008) reported the seasonal evolution of Ba xs profiles on a monthly basis between February and June 2003. These authors showed the mesopelagic Ba xs build up (and barite stock increase) following the spring phytoplankton bloom development, enhanced POC fluxes and subsequent remineralization. Overall, the DWA Ba xs values reported in the present study were higher than those reported by Sternberg et al. (2008) (maximum of 463 pM; 0-600 m). The variability over the 2 d period at Tyrr station was of the same order of magnitude as the seasonal DWA Ba xs dynamics found at DYFAMED and similar to changes found over a few days to a week-long period in different sectors of the Southern Ocean (Cardinal et al., 2005;Jacquet et al., 2007aJacquet et al., , 2015. The column-integrated PHP vs. DWA Ba xs ratio at Tyrr station confirms that the second visit experienced higher remineralization rates in the upper mesopelagic layer than during the first one (Table 1).

Mesopelagic C remineralization
POC remineralization rates (MRs) estimated from DWA Ba xs values using Eq.
(1) are shown in Fig. 5 for the upper (100-500 m) and entire (100-1000 m) mesopelagic layer along with primary productivity (Van Wambecke et al., 2021). The MRs ranged from 25 ± 2 to 306 ± 70 mg C m −2 d −1 , and primary production ranged from 138 to 284 mg C m −2 d −1 . A large difference in MRs between the upper and the whole mesopelagic layers can be seen in the ALG Basin. This is more pronounced at station no. 9, which has MRs of 91 mg C m −2 d −1 in the upper (100 to 500 m depth) layer and 306 mg C m −2 d −1 in the entire mesopelagic layer (Fig. 5).
These results show that significant remineralization occurred between 500 and 1000 m in the ALG Basin in contrast to the ION and TYR basins, where remineralization mainly occurred in the mesopelagic layer between 100 and 500 m depth. A similar conclusion was reached by Jullion et al. (2017) from dissolved Ba and parametric optimum multiparameter (POMP) derived POC remineralization rates along a zonal transect between the Lebanese coast and Gibraltar (from 156 to 348 mg C m −2 d −1 ; M84/3 cruise, April 2011). Independent of any dust input considerations, Jullion et al. (2017) showed significant differences in the mesopelagic MRs between the western and eastern Mediterranean, indicating an additional organic carbon export pathway to depth. The western basin is indeed the site of deep-shelf and open-ocean convection, transferring organic matter to deeper layers (Durrieu de Madron et al., 2013;Stabholz et al., 2013). The larger MR fluxes found in the ALG Basin during PEACETIME are in line with an ecoregion with the recurrent injection of material by winter convection (the particle injection pump hypothesis; Boyd et al., 2019), sustaining higher rates of remineralization below 500 m depth. In contrast, in the TYR Basin, remineralization was mainly located in the upper mesopelagic layer. Stations in the TYR Figure 4. Ratio of surface layer integrated PHP (Int.PHPx1) to mesopelagic integrated PHP (IntPHPx2) vs. mesopelagic depth-weighted average (DWA) Ba xs (pM) during PEACETIME. The same data are reported for the KEOPS1 and KEOPS2 cruises (Southern Ocean; Jacquet et al., 2015) and at the PAP (northeastern Atlantic -"N Atlantic") and ANTARES/EMSO-LO (northwestern Mediterranean Sea -"NW Medsea") observatory sites . The blue line (R 2 = 0.88) represents the trend reported during KEOPS2 (Jacquet et al., 2005). The data points referred to as "season advancement" (encircled by the blue line) were excluded from the KEOPS2 regression analysis shown here.
Basin received dust inputs a few days before our arrival at these locations; the particulate Al concentrations and estimated lithogenic Ba fraction reflected the impact of this event ( Fig. 2; Bressac et al., 2021). At Tyrr station, the DAW Ba xs vs. column-integrated PHP increase between the two visits indicated higher MRs. MRs were mainly localized in the upper 500 m. Another atmospheric deposition event occurred on 5 June, a few hours after the first sampling at Fast station in the ALG Basin. However, Fast station does not present any evidence of an impact from particulate Al concentrations and estimated lithogenic Ba at mesopelagic depths. In contrast to conditions in the surface mixed layer, the generation of an observable signal from mesopelagic remineralization and subsequent Ba xs formation from a single dust event would require more time than the time span between atmospheric deposition and sampling at Fast station (in contrast to Tyrr station where the dust event occurred 5-12 d before sampling). In the ION Basin, where stations did not reflect the impact of any deposition event and were not subject to potential deep convection, DWA Ba xs and MR fluxes were mostly restricted to the upper mesopelagic layer. Berline et al. (2021) report small-scale heterogeneity in particle abundances at ION stations, emphasizing the spatial decoupling between particle production and particle distribution and adding complexity to estimating the time lag between the production and export of particles and, thus, C transfer at depth (Stange et al., 2017;Henson et al., 2011). Further, no significant surface production events occurred in the ION Basin However, surface particles at station no. 8 seemed related to a past production event without significant vertical export by the time the station was sampled. As reported in Van Wambeke et al. (2020), primary production fluxes were slightly higher in the ION Basin (from 158 to 208 mg C m −2 d −1 ) than in the TYR Basin (from 142 to 170 mg C m −2 d −1 ). Thus, DWA Ba xs values and MR fluxes reported in the ION Basin would generally reflect the earlier stage of export and remineralization processes. The same applies to Tyrr2 station (in contrast to Tyrr4 station) according to the DWA Ba xs vs. integrated-PHP trend.

Conclusions
The present paper expands the data coverage of the Ba xs distribution in the ALG, TYR and ION basins (western and central Mediterranean Sea) in late spring 2017. Results high-light that mesopelagic remineralization processes are mainly located in the upper 500 m horizon in the TYR and ION basins, whereas they occur in the lower mesopelagic zone (down to 1000 m) in the ALG Basin. We suggest that particle injection driven by the seasonal winter deep convection in the western basin would sustain the larger and deeper MRs that we observed in the ALG Basin. In both the TYR and ION basins, Ba xs indicated lower (intensity) and upper mesopelagic-layer-restricted remineralization processes that could be the results of a previous dust deposition event (in particular at Tyrr station) or the patchiness of time lags between production and export of particles.
Author contributions. SJ wrote the paper with contributions from all co-authors. SJ and AD managed the barium analyses, and CT, MG, SG and FVV managed the PHP analyses. MG, SG and NB performed Ba sampling during the cruise.
Competing interests. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Special issue statement. This article is part of the special issue "Atmospheric deposition in the low-nutrient-low-chlorophyll (LNLC) ocean: effects on marine life today and in the future (ACP/BG interjournal SI)". It is not associated with a conference.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.