The recent roadmap IndOOS-2 has stressed the need to expand the biogeochemical-Argo observing system in the Indian Ocean. The Monaco Explorations Indian Ocean expedition offered a unique opportunity to meet this goal in the southwestern sector which was, in this regard and at that time, one of the least covered oceanic regions. We designed a deployment strategy for the biogeochemical float array grounded on past experiences, existing knowledge, and the analysis of historical datasets to cover the contrasting biophysical regimes from the Seychelles Chagos Thermocline Ridge to the subtropical gyre. Aligning with IndOOS-2 recommendations, a denser float distribution was set in the tropical band to enhance biogeochemical observations in upwelling zones. Following this strategy, a fleet of seventeen biogeochemical floats was successfully deployed during the expe- dition in October–November 2022. After two years of operations, the spatio-temporal distribution covered by the fleet confirmed that the goals of the deployment strategy have been reached, revealing seasonal modulations of the meridional trophic gradient with respect to phytoplankton biomass from tropical mesotrophy to subtropical oligotrophy

Abstract. This study examines conductivity–temperature–depth (CTD), acoustic Doppler current profiler (ADCP) and optical data collected during the PERLE-3 cruise in March 2020 between the surface and 2000 m depth over an east–west section of the Tyrrhenian Sea in the Mediterranean. The focus is on the impact of double-diffusion processes, in particular salt fingering, on the distribution and dynamics of particulate and dissolved matter. Staircases form at the interface between the warm, saline Eastern Intermediate Water (EIW) and the colder, less saline Tyrrhenian Deep Water (TDW) in the centre of the basin in a region of weak hydrodynamic activity. The results show that thermohaline staircases formed by salt fingering significantly influence particle sedimentation and biogeochemical cycling in deep ocean environments by altering vertical flux patterns. These density steps create distinct vertical layers that act as physical barriers, slowing the descent of particles and facilitating their retention and aggregation. The retention of fine particles at density gradients promotes the formation of larger aggregates, thereby affecting the particle size distribution. The staircases also affect dissolved matter by creating pronounced concentration gradients of oxygen and nutrients, which may influence microbial activity and nutrient fluxes.

During the TONGA cruise (2019), seawater samples were collected to assess the effect of volcanic eruption versus submarine hydrothermal system on the water column. For this purpose, two locations were investigated, the first one located directly under the influence of the New Late'iki island (eruption in October 2019), and the second one showing ongoing submarine hydrothermal activity. At both locations, the total strontium (TSr) and lithium (TLi) concentrations vary between 94.4 and 152.3 µmol/L and 13.2 and 203.5 µmol/L, respectively. When combined, TSr and TLi concentrations of all samples in the water column are higher than those of the oligotrophic water. Both volcanic eruption and submarine hydrothermal activity (e.g. volcanic ashes, particles, gas condensate) can deliver substantial amount of TSr and TLi to the water column. The distribution of TSr versus TLi evidences linear trends either with a negative or positive slope. The negative correlation is observed in the water column at both sites, directly under the influence of the eruption and in the vicinity of the volcano with hydrothermal activity. The positive TSr versus TLi correlation is observed at site under submarine hydrothermal influence and is in line with black smokers related hydrothermal plumes. The 87 Sr/ 86 Sr ratios vary between 0.709147 and 0.709210 and d 7 Li values vary between +10.1 and +37.6 ‰. While 92% of the measured 87 Sr/ 86 Sr ratios are in line with the mean value of oligotrophic waters, once combined with the d 7 Li values, only 20% of them remains within this field. The wide range of d 7 Li values decreases from sea-surface down to ~140 mbsl, before increasing at greater depth, while defining different linear trend according to the dissolved inorganic carbon concentrations. The variability of d 7 Li values reflect hydrothermal contribution, mineral-seawater interaction and potentially biology-environment interaction. In the particular geological setting of the study, where both hydrothermal and volcanic activities were at play, disentangling both contributions on water column implies a combined use of elemental and isotopic signatures of Sr and Li tracers.

Diazotrophs regulate marine productivity in the oligotrophic ocean by alleviating nitrogen limitation, contributing to particulate organic carbon (POC) export to the deep ocean. Yet, the characterization of particles composing the sinking POC flux has never been explored in such ecosystems. Moreover, the contribution of the direct gravitational export of diazotrophs to the overall flux is seldom assessed. Here we explore the composition of the sinking POC flux in a hot spot of N 2 fixation (the western sub-tropical South Pacific) using polyacrylamide gel-filled traps deployed at two stations (S05M and S10M) and three depths (170 m, 270 m, 1000 m) during the TONGA expedition (November-December 2019). Image analyses of particles collected in the gels was used to classify them into 5 categories (fecal aggregates, phytodetrital aggregates, mixed aggregates, cylindrical fecal pellets, and zooplankton carcasses). Fecal aggregates were the most abundant at both stations and all depths and dominated the flux (average of 56 ± 28% of the POC flux), followed by zooplankton carcasses (24 ± 19%), cylindrical fecal pellets (15 ± 14%) and mixed aggregates (5 ± 4%), whereas phytodetrital aggregates contributed less (<1%). Since N isotope budgets show that export is mainly supported by diazotrophy at these stations, these results suggest that the diazotroph-derived N has been efficiently transferred to the foodweb up to zooplankton and fecal pellets before being exported, pleading for an indirect export of diazotrophy. However, random confocal microscopy examination performed on sinking particles revealed that diazotrophs were present in several categories of exported particles, suggesting that diazotrophs are also directly exported, with a potential contribution to overall POC fluxes increasing with depth. Our results provide the first characterization of particle categories composing the sinking flux and their contribution to the overall flux in a hot spot of N 2 fixation.

The Arabian Sea is a highly productive tropical ecosystem of the Indian Ocean that supports high fluxes of particulate organic carbon to the mesopelagic zone from two distinct periods of elevated biological productivity associated with the semiannual reversals of the monsoonal wind system. There are now strong indications that the Arabian Sea's monsoonal wind patterns and hydrographic conditions are being impacted by long-term temperature increases, but the consequences of these changes on primary production and carbon export to the mesopelagic zone are unknown. This is especially true for the summer monsoon period when cloud cover obscures much of the Arabian Sea basin and therefore precludes remotely sensed ocean color measurements for estimating phytoplankton biomass and productivity. Here we overcome this limitation by using a database of bio-optical profiles from Biogeochemical Argo floats collected over the last decade to evaluate the impact of interannual temperature increases on Arabian Sea primary production and carbon export. We classify individual years of float observations based on the spatial extent of the Arabian Sea Mini Warm Pool that appears in the southeast Arabian Sea before the onset of the summer monsoon. This Mini Warm Pool, which begins to build in winter and collapses with the onset of the summer monsoon in late spring, has gained considerable interest on account of its influence on the timing of the onset of the summer monsoon. We observed a 35 percent decrease in primary production during the summer monsoon phytoplankton bloom in strong warm pool years, and a 13 percent decrease in particle stocks in the upper mesopelagic zone following the peak of the bloom. Decreases in production and export were additionally accompanied by a decrease in average particle size, indicating a shift from larger cells like diatoms that appear from fertilization of the oligotrophic waters to smaller phytoplankton size classes in response to a deepening of the thermocline and increased stratification of the water column. These results suggest changes in phytoplankton community structure and further decreases in primary production and carbon export in the Arabian Sea in response to future warming.

Abstract. Hydrolysis of dissolved organic phosphorus by marine planktonic microorganisms is a key process in the P cycle, particularly in P-depleted, oligotrophic environments. The present study assessed spatiotemporal variations in phosphomonoesterase (PME) and phosphodiesterase (PDE) activities using concentration kinetics in the eastern Mediterranean Sea in two contrasting situations: the end of winter (including a small bloom period) and autumn. The distribution and regulation of the maximum hydrolysis rate (Vm) and half-saturation constant (Km) of both ectoenzymes were assessed in relation to the vertical structure of the epipelagic layers. PME reached its maximum activities (Vm) after the addition of 1 µM MUF-P (4-methylumbelliferyl phosphate), whereas, for PDE, it was necessary to add up to 50 µM bis(4-methylumbelliferyl)phosphate (bis-MUF-P) to reach saturation state. On average, the Km of PDE was 33 ± 25 times higher than that of PME. The Vm of PME and Vm of PDE were linearly correlated. Conversely to the Km values, Vm values were on the same order of magnitude for both ectoenzymes, with their ratio (Vm PME : Vm PDE) ranging between 0.2 and 6.3. Dissolved organic phosphorus (DOP) and the phosphomonoesterase hydrolysable fraction of DOP explained most of the lack of variability in Vm PME and Vm PDE. On the contrary, Vm of both phosphohydrolase enzymes was inversely correlated to the concentration of dissolved inorganic phosphorus. The particular characteristics of concentration kinetics obtained for PDE (saturation at 50 µM, high Km, high turnover times) are discussed with respect to the possible unequal distribution of PDE and PME among the size continuum of organic material and accessibility of phosphodiesters.

The objective of the TONGA oceanographic expedition was to study the control of productivity and carbon sequestration by micronutrients of shallow hydrothermal origin in the Western Tropical South Pacific (WTSP) Ocean. The 37-day oceanographic survey took place on board the R/V L’Atalante in 2019 between Oct. 31 to Dec. 6 (Nouméa-Nouméa). Over a large area of the WTSP the team acquired numerous results on both the entire water column (up to the sediment) and the atmosphere. Specific task are represented on figure 1: (task 1) to characterize chemically and optically shallow hydrothermal fluids and to compare the source from below (shallow hydrothermal fluids) with the source from above (atmospheric deposition); (task 2) to quantify the dynamical dispersion of the fluids at small and regional scale; (task 3) to investigate the impact of the shallow hydrothermal sources on the biological activity and diversity, and the feedback to the atmosphere via the oceanic emissions of primary and secondary aerosols. (Task 4) to communicate about the campaign (see for example our Tweeter account (https://twitter.com/tongaproject) and the movie (26’) both in French (https://www.youtube.com/watch?v=e5kAd0i6Dck) and English (https://www.youtube.com/watch?v=UeABf-cVR-k). A long west to east (up to the blue waters of the gyre) transect allowed to characterize the different biogeochemical provinces crossed and a focus in the region of the Lau Basin allowed to investigate the impact of shallow hydrothermal sources. A series of short and long stations allowed to fully characterize the stocks and the fluxes in the different provinces. Short-term (up to 10 days) processes studies have been conducted (drifting moorings and minicosms experiments). Part of these results will feed into important modeling work. A fixed mooring line launched at the end of the campaign and recovered in Nov. 2020 as well as the 7 ARGO floats and 20 drifting buoys that were dropped during the campaign provide a broader temporal context of the acquisitions done during the campaign. An important focus of the campaign was the trace metal characterization of the entire water column. For this, TONGA has been labeled by the international program GEOTRACES (https://www.geotraces.org/). The impact on biological communities of fluids is supported by the international IMBER program (https://imber.info/). The TONGA project is also part of the LEFE program (funding by LEFE-CYBER and LEFE-GMMC), the ANR (Appel à projets génériques) and the Fondation A-MIDeX of the Aix-Marseille Université.

Abstract. Numerical models of ocean biogeochemistry are becoming the major tools used to detect and predict the impact of climate change on marine resources and to monitor ocean health. However, with the continuous improvement of model structure and spatial resolution, incorporation of these additional degrees of freedom into fidelity assessment has become increasingly challenging. Here, we propose a new method to provide information on the model predictive skill in a concise way. The method is based on the conjoint use of a k-means clustering technique, assessment metrics, and Biogeochemical-Argo (BGC-Argo) observations. The k-means algorithm and the assessment metrics reduce the number of model data points to be evaluated. The metrics evaluate either the model state accuracy or the skill of the model with respect to capturing emergent properties, such as the deep chlorophyll maximums and oxygen minimum zones. The use of BGC-Argo observations as the sole evaluation data set ensures the accuracy of the data, as it is a homogenous data set with strict sampling methodologies and data quality control procedures. The method is applied to the Global Ocean Biogeochemistry Analysis and Forecast system of the Copernicus Marine Service. The model performance is evaluated using the model efficiency statistical score, which compares the model–observation misfit with the variability in the observations and, thus, objectively quantifies whether the model outperforms the BGC-Argo climatology. We show that, overall, the model surpasses the BGC-Argo climatology in predicting pH, dissolved inorganic carbon, alkalinity, oxygen, nitrate, and phosphate in the mesopelagic and the mixed layers as well as silicate in the mesopelagic layer. However, there are still areas for improvement with respect to reducing the model–data misfit for certain variables such as silicate, pH, and the partial pressure of CO2 in the mixed layer as well as chlorophyll-a-related, oxygen-minimum-zone-related, and particulate-organic-carbon-related metrics. The method proposed here can also aid in refining the design of the BGC-Argo network, in particular regarding the regions in which BGC-Argo observations should be enhanced to improve the model accuracy via the assimilation of BGC-Argo data or process-oriented assessment studies. We strongly recommend increasing the number of observations in the Arctic region while maintaining the existing high-density of observations in the Southern Oceans. The model error in these regions is only slightly less than the variability observed in BGC-Argo measurements. Our study illustrates how the synergic use of modeling and BGC-Argo data can both provide information about the performance of models and improve the design of observing systems.

Iron is an essential nutrient that regulates productivity in ~30% of the ocean. Compared with deep (>2000 meter) hydrothermal activity at mid-ocean ridges that provide iron to the ocean’s interior, shallow (<500 meter) hydrothermal fluids are likely to influence the surface’s ecosystem. However, their effect is unknown. In this work, we show that fluids emitted along the Tonga volcanic arc (South Pacific) have a substantial impact on iron concentrations in the photic layer through vertical diffusion. This enrichment stimulates biological activity, resulting in an extensive patch of chlorophyll (360,000 square kilometers). Diazotroph activity is two to eight times higher and carbon export fluxes are two to three times higher in iron-enriched waters than in adjacent unfertilized waters. Such findings reveal a previously undescribed mechanism of natural iron fertilization in the ocean that fuels regional hotspot sinks for atmospheric CO 2 .

Climatic changes and interannual variability in the Mediterranean overturning circulation are crucially linked to dense water formation in the Levantine Sea, namely the Levantine Intermediate Water whose formation zone, comprising multiple and intermittent sources, extends over fluctuating pathways. To probe into the variability of this water formation and spreading, a unique dataset was collected during the winter of 2019 in the western Levantine Sea, via oceanographic cruises, profiling floats and a glider, at a spatio-temporal distribution suited to resolve mesoscale circulation features and intermittent convection events. This study highlights the competition between two source regions, the Cretan Sea and the Rhodes Cyclonic Gyre, to supply the Mediterranean overturning circulation in Levantine Intermediate Water. The Cretan source was estimated as the most abundant, supported by increasingly saltier water masses coming from the Levantine Sea under the pumping effect of a water deficit caused by strong western outflow toward the Ionian Sea.

N2 fixation rates were measured in the 0–1000 m layer at 13 stations located in the open western and central Mediterranean Sea (MS) during the PEACETIME cruise (late spring 2017). While the spatial variability of N2 fixation was not related to Fe, P nor N stocks, the surface composition of the diazotrophic community indicated a strong eastward increasing longitudinal gradient for the relative abundance of non-cyanobacterial diazotrophs (NCD) (mainly γ-Proteobacteria) and conversely eastward decreasing for UCYN-A (mainly -A1 and -A3) as did N2 fixation rates. UCYN-A4 and A3 were identified for the first time in the MS. The westernmost station influenced by Atlantic waters, and characterized by highest stocks of N and P, displayed a patchy distribution of diazotrophic activity with an exceptionally high rate in the euphotic layer of 72.1 nmol N L−1 d−1, which could support up to 19 % of primary production. At this station at 1 %PAR depth, UCYN-A4 represented up to 94 % of the diazotrophic community. These in situ observations of higher UCYN-A relative abundance in nutrient rich stations while NCD increased in the more oligotrophic stations, suggest that the nutrient conditions could determine the composition of the diazotrophic communities and in turn the N2 fixation rates. The impact of Saharan dust deposition on N2 fixation and diazotrophic communities was also investigated, under present and future projected conditions of temperature and pH during short term (3–4 days) experiments at three stations. New nutrients from simulated dust deposition triggered a significant stimulation of N2 fixation (from 41 % to 565 %). The strongest increase in N2 fixation was observed at the stations dominated by NCD and did not lead on this short time scale to change in the diazotrophic community composition. Under projected future conditions, N2 fixation was either exacerbated or unchanged, in that later case this was probably due to a too low nutrient bioavailability or an increased grazing pressure. The future warming and acidification likely benefited NCD (Pseudomonas) and UCYN-A2 while disadvantaged UCYN-A3 without knowing which effect (alone or in combination) is the driver, especially since we do not know the temperature optima of these species not yet cultivated as well as the effect of acidification.

In the Western Tropical South Pacific, a hotspot of dinitrogen-fixing organisms has been identified. The survival of these species depends on the availability of dissolved iron (DFe); however, the source of this DFe is still unclear. DFe was measured along a transect from 175°E to 166°W near 19–21°S. The distribution of DFe showed high spatial variability: low concentrations (∼0.2 nmol kg−1) in the South Pacific gyre and high concentrations (up to 50 nmol kg−1) in the west of the Tonga arc, indicating that this arc is a clear boundary between iron-poor and iron-rich waters. An optimal multiparameter analysis was used to distinguish the relative importance of physical transport relative to non-conservative processes on the observed distribution. This analysis demonstrated that the shallow hydrothermal sources present along the Tonga-Kermadec arc are responsible for the high concentrations observed in the photic layer. Nevertheless, in contrast to what has been observed for deep hydrothermal plumes, our results highlighted the rapid decrease in DFe concentrations near shallow hydrothermal sources. This is likely due to a shorter residence time of surface water masses combined with several biogeochemical processes at play (precipitation, scavenging, biological uptake, and photoreduction). This study clearly highlights the role of shallow hydrothermal sources on the DFe cycle within the Tonga-Kermadec arc where a strong link to biological activity in surface waters can be assessed, despite the small but significant fraction of DFe ultimately stabilized. It also emphasizes the need to consider the impact of these sources for a better understanding of the global iron cycle.

This study assesses marine community production based on the diel variability of bio-optical properties monitored by two BioGeoChemical-Argo (BGC-Argo) floats. Experiments were conducted in two distinct Mediterranean systems, the northwestern Ligurian Sea and the central Ionian Sea, during summer months. We derived particulate organic carbon (POC) stock and gross community production integrated within the surface, euphotic and subsurface chlorophyll maximum (SCM) layers, using an existing approach applied to diel cycle measurements of the particulate beam attenuation (c_p) and backscattering (b_bp) coefficients. The diel cycle of c_p provided a robust proxy for quantifying biological production in both systems; that of b_bp was comparatively less robust. Derived primary production estimates vary by a factor of 2 depending upon the choice of the bio-optical relationship that converts the measured optical coefficient to POC, which is thus a critical step to constrain. Our results indicate a substantial contribution to the water column production of the SCM layer (16 %-42 %), which varies largely with the considered system. In the Ligurian Sea, the SCM is a seasonal feature that behaves as a subsurface biomass maximum (SBM) with the ability to respond to episodic abiotic forcing by increasing production. In contrast, in the Ionian Sea, the SCM is permanent, primarily induced by phytoplankton photoacclimation, and contributes moderately to water column production. These results clearly demonstrate the strong potential for transmissometers deployed on BGC-Argo profiling floats to quantify non-intrusively in situ biological production of organic carbon in the water column of stratified oligotrophic systems with recurring or permanent SCMs, which are widespread features in the global ocean.

Abstract. N2 fixation rates were measured in the 0–1000 m layer at 13 stations located in the open western and central Mediterranean Sea (MS) during the PEACETIME cruise (late spring 2017). While the spatial variability in N2 fixation was not related to Fe, P nor N stocks, the surface composition of the diazotrophic community indicated a strong longitudinal gradient increasing eastward for the relative abundance of non-cyanobacterial diazotrophs (NCDs) (mainly γ-Proteobacteria) and conversely decreasing eastward for photo-heterotrophic group A (UCYN-A) (mainly UCYN-A1 and UCYN-A3), as did N2 fixation rates. UCYN-A4 and UCYN-A3 were identified for the first time in the MS. The westernmost station influenced by Atlantic waters and characterized by highest stocks of N and P displayed a patchy distribution of diazotrophic activity with an exceptionally high rate in the euphotic layer of 72.1 nmolNL-1d-1, which could support up to 19 % of primary production. At this station at 1 % PAR (photosynthetically available radiation) depth, UCYN-A4 represented up to 94 % of the diazotrophic community. These in situ observations of greater relative abundance of UCYN-A at stations with higher nutrient concentrations and dominance of NCDs at more oligotrophic stations suggest that nutrient conditions – even in the nanomolar range – may determine the composition of diazotrophic communities and in turn N2 fixation rates. The impact of Saharan dust deposition on N2 fixation and diazotrophic communities was also investigated, under present and future projected conditions of temperature and pH during short-term (3–4 d) experiments at three stations. New nutrients from simulated dust deposition triggered a significant stimulation of N2 fixation (from 41 % to 565 %). The strongest increase in N2 fixation was observed at the stations dominated by NCDs and did not lead on this short timescale to changes in the diazotrophic community composition. Under projected future conditions, N2 fixation was either increased or unchanged; in that later case this was probably due to a too-low nutrient bioavailability or an increased grazing pressure. The future warming and acidification likely benefited NCDs (Pseudomonas) and UCYN-A2, while disadvantaged UCYN-A3 without knowing which effect (alone or in combination) is the driver, especially since we do not know the temperature optima of these species not yet cultivated as well as the effect of acidification.

Abstract. Mineral dust deposition is an important supply mechanism for trace elements in the low-latitude ocean. Our understanding of the controls of such inputs has been mostly built on laboratory and surface ocean studies. The lack of direct observations and the tendency to focus on near-surface waters prevent a comprehensive evaluation of the role of dust in oceanic biogeochemical cycles. In the frame of the PEACETIME project (ProcEss studies at the Air-sEa Interface after dust deposition in the MEditerranean sea), the responses of the aluminum (Al) and iron (Fe) cycles to two dust wet deposition events over the central and western Mediterranean Sea were investigated at a timescale of hours to days using a comprehensive dataset gathering dissolved and suspended particulate concentrations, along with sinking fluxes. Dissolved Al (dAl) removal was dominant over dAl released from dust. The Fe/Al ratio of suspended and sinking particles revealed that biogenic particles, and in particular diatoms, were key in accumulating and exporting Al relative to Fe. By combining these observations with published Al/Si ratios of diatoms, we show that adsorption onto biogenic particles, rather than active uptake, represents the main sink for dAl in Mediterranean waters. In contrast, systematic dissolved Fe (dFe) accumulation occurred in subsurface waters (∼ 100–1000 m), while dFe input from dust was only transient in the surface mixed layer. The rapid transfer of dust to depth, the Fe-binding ligand pool in excess to dFe in subsurface (while nearly saturated in surface), and low scavenging rates in this particle-poor depth horizon are all important drivers of this subsurface dFe enrichment. At the annual scale, this previously overlooked mechanism may represent an additional pathway of dFe supply for the surface ocean through diapycnal diffusion and vertical mixing. However, low subsurface dFe concentrations observed at the basin scale (< 0.5 nmol kg−1) cause us to question the residence time for this dust-derived subsurface reservoir and hence its role as a supply mechanism for the surface ocean, stressing the need for further studies. Finally, these contrasting responses indicate that dAl is a poor tracer of dFe input in the Mediterranean Sea.

Abstract. The study of phosphorus cycling in phosphate-depleted oceanic regions, such as the Mediterranean Sea, has long suffered from methodological limitations, leading to a simplistic view of a homogeneous surface phosphate pool with concentrations below the detection limit of measurement above the phosphacline. During the PEACETIME (Process studies at the air-sea interface after dust deposition in the Mediterranean Sea) cruise, carried out from 10 May to 11 June 2017, we conducted co-located measurements of phosphate pools at the nanomolar level, alkaline phosphatase activities and atmospheric deposition of phosphorus, across a longitudinal gradient from the west to the central Mediterranean Sea. In the phosphate-depleted layer (PDL), between the surface and the phosphacline, nanomolar phosphate was low and showed little variability across the transect spanning from 6 ± 1 nmol L−1 in the Ionian basin to 15 ± 4 nmol L−1 in the westernmost station. The low variability in phosphate concentration contrasted with that of alkaline phosphatase activity, which varied over 1 order of magnitude across the transect. Nanomolar phosphate data revealed gradients of phosphate concentration over density inside the PDL ranging between 10.6 ± 2.2 µmol kg−1 in the westernmost station to values close to zero towards the east. Using the density gradients, we estimated diapycnal fluxes of phosphate to the PDL and compared them to atmospheric deposition, another external source of phosphate to the PDL. Phosphate supply to the PDL from dry deposition and diapycnal fluxes was comparable in the western part of the transect. This result contrasts with the longtime idea that, under stratification conditions, the upper waters of the Mediterranean Sea receive new P almost exclusively from the atmosphere. The contribution of atmospheric deposition to external P supply increased under the occurrence of rain and Saharan dust. Although this finding must be taken cautiously given the uncertainties in the estimation of diapycnal fluxes, it opens exciting questions on the biogeochemical response of the Mediterranean Sea, and more generally of marine oligotrophic regions, to expected changes in atmospheric inputs and stratification regimes. Taken together, external sources of phosphate to the PDL contributed little to total phosphate requirements which were mainly sustained by in situ hydrolysis of dissolved organic phosphorus. The results obtained in this study show a highly dynamic phosphorus pool in the upper layer of the euphotic zone, above the phosphacline, and highlight the convenience of combining highly sensitive measurements and high-resolution sampling to precisely depict the shape of phosphate profiles in the euphotic zone with still unexplored consequences on P fluxes supplying this crucial layer for biogeochemical cycles.

The surface mixed layer (ML) in the Mediterranean Sea is a well-stratified domain characterized by low macronutrients and low chlorophyll content for almost 6 months of the year. In this study we characterize the biogeochemical cycling of nitrogen (N) in the ML by analyzing simultaneous in situ measurements of atmospheric deposition, nutrients in seawater, hydrological conditions, primary production, heterotrophic prokaryotic production, N2 fixation and leucine aminopeptidase activity. Dry deposition was continuously measured across the central and western open Mediterranean Sea, and two wet deposition events were sampled, one in the Ionian Sea and one in the Algerian Basin. Along the transect, N budgets were computed to compare the sources and sinks of N in the mixed layer. In situ leucine aminopeptidase activity made up 14 % to 66 % of the heterotrophic prokaryotic N demand, and the N2 fixation rate represented 1 % to 4.5 % of the phytoplankton N demand. Dry atmospheric deposition of inorganic nitrogen, estimated from dry deposition of nitrate and ammonium in aerosols, was higher than the N2 fixation rates in the ML (on average 4.8-fold). The dry atmospheric input of inorganic N represented a highly variable proportion of biological N demand in the ML among the stations, 10 %–82 % for heterotrophic prokaryotes and 1 %–30 % for phytoplankton. As some sites were visited on several days, the evolution of biogeochemical properties in the ML and within the nutrient-depleted layers could be followed. At the Algerian Basin site, the biogeochemical consequences of a wet dust deposition event were monitored through high-frequency sampling. Notably, just after the rain, nitrate was higher in the ML than in the nutrient-depleted layer below. Estimates of nutrient transfer from the ML into the nutrient-depleted layer could explain up to a third of the nitrate loss from the ML. Phytoplankton did not benefit directly from the atmospheric inputs into the ML, probably due to high competition with heterotrophic prokaryotes, also limited by N and phosphorus (P) availability at the time of this study. Primary producers decreased their production after the rain but recovered their initial state of activity after a 2 d lag in the vicinity of the deep chlorophyll maximum layer.

During the winter from 2009 to 2013, the mixed layer reached the seafloor at about 2500min the northwestern Mediterranean Sea. Intense fronts around the deep convection area were repeatedly sampled by autonomous gliders. Subduction down to 200–300 m, sometimes deeper, below the mixed layer was regularly observed testifying of important frontal vertical movements. Potential Vorticity dynamics was diagnosed using glider observations and a high resolution realistic model at 1-km resolution. During down-front wind events in winter, remarkable layers of negative PV were observed in the upper 100m on the dense side of fronts surrounding the deep convection area and successfully reproduced by the numerical model. Under such conditions, symmetric instability can grow and overturn water along isopycnals within typically 1–5 km cross-frontal slanted cells. Two important hotpspots for the destruction of PV along the topographically-steered Northern Current undergoing frequent down-front winds have been identified in the western part of Gulf of Lion and Ligurian Sea. Fronts were there symmetrically unstable for up to 30 days per winter in the model, whereas localized instability events were found in the open sea, mostly influenced by mesoscale variability. The associated vertical circulations also had an important signature on oxygen and fluorescence, highlighting their under important role for the ventilation of intermediate layers, phytoplankton growth and carbon export.

Abstract. Numerical models of ocean biogeochemistry are becoming a major tool to detect and predict the impact of climate change on marine resources and ocean health. Classically, validation of such models relies on comparison with surface quantities from satellite (such as chlorophyll-a concentrations), climatologies, or sparse in situ data (such as cruises observations, and permanent fixed oceanic stations). However, these datasets are not fully suitable to assess how models represent many climate-relevant biogeochemical processes. These limitations now begin to be overcome with the availability of a large number of vertical profiles of light, pH, oxygen, nitrate, chlorophyll-a concentrations and particulate backscattering acquired by the Biogeochemical-Argo (BGC-Argo) floats network. Additionally, other key biogeochemical variables such as dissolved inorganic carbon and alkalinity, not measured by floats, can be predicted by machine learning-based methods applied to float oxygen concentrations. Here, we demonstrate the use of the global array of BGC-Argo floats for the validation of biogeochemical models at the global scale. We first present 18 key metrics of ocean health and biogeochemical functioning to quantify the success of BGC model simulations. These metrics are associated with the air-sea CO2 flux, the biological carbon pump, oceanic pH, oxygen levels and Oxygen Minimum Zones (OMZs). The metrics are either a depth-averaged quantity or correspond to the depth of a particular feature. We also suggest four diagnostic plots for displaying such metrics.

The seasonal variability of the carbonate system in the eastern Mediterranean Sea (EMed) was investigated based on discrete total alkalinity (A T ), total dissolved inorganic carbon (C T ), and pH measurements collected during three cruises around Crete between June 2018 and March 2019. This study presents a detailed description of this new carbonate chemistry dataset in the eastern Mediterranean Sea. We show that the North Western Levantine Basin (NWLB) is unique in terms of range of A T variation vs. C T variation in the upper water column over an annual cycle. The reasons for this singularity of the NWLB can be explained by the interplay between strong evaporation and the concomitant consumption of C T by autotrophic processes. The high range of A T variations, combined to temperature changes, has a strong impact on the variability of the seawater p CO 2 ( p CO 2 S W ). Based on Argo float data, an entire annual cycle for p CO 2 S W in the NWLB has been reconstructed in order to estimate the temporal sequence of the potential “source” and “sink” of atmospheric CO 2 . By combining this dataset with previous observations in the NWLB, this study shows a significant ocean acidification and a decrease in the oceanic surface pH T 25 of −0.0024 ± 0.0004 pH T 25 units.a –1 . The changes in the carbonate system are driven by the increase of atmospheric CO 2 but also by unexplained temporal changes in the surface A T content. If we consider that the EMed will, in the future, encounter longer, more intense and warmer summer seasons, this study proposes some perspectives on the carbonate system functioning of the “future” EMed.

The extension of Argo float coverage to the European marginal seas is one of the strategic targets of the Euro-Argo European Research Infrastructure Consortium (ERIC). Under this general framework, the Argo capability to monitor the shallow coastal shelf remains an open question. In the Euro-Argo RISE H2020 project, targeted deployments have been undertaken to investigate this potential. In this study, we present the experience and outcomes from 4 such deployments in areas of the Mediterranean Sea with intrigue coastlines and complex bathymetry (north Aegean, north Adriatic, south Palma, and Gulf of Lions). We focus on the floats’ configuration settings and the monitoring tools/software that have been utilized to follow the floats’ performance. Our results show that certain configuration parameters such as the drifting depth, and the sampling frequency, play a significant role in the floats‘ performance. Technological advances both on the floats’ characteristics and on the monitoring-controlling tools can lead to significant improvements of similar missions in the near future. The fact that all fl oats achieved successful missions and acquired data showing important hydrographic features, highlights the importance of Argo expansion in targeted shallow coastal areas where Argo can be a complementary part of an integrated oceanographic monitoring system for the Mediterranean Sea.

The Mediterranean Sea is a hotspot for climate change, and recent studies have reported its intense warming and salinification. In this study, we use an outstanding dataset relying mostly on glider endurance lines but also on other platforms to track these trends in the northwestern Mediterranean where deep convection occurs. Thanks to a high spatial coverage and a high temporal resolution over the period 2007-2017, we observed the warming (+0.06 [Formula: see text]C year[Formula: see text]) and salinification (+0.012 year[Formula: see text]) of Levantine Intermediate Water (LIW) in the Ligurian Sea. These rates are similar to those reported closer to its formation area in the Eastern Mediterranean Sea. Further downstream, in the Gulf of Lion, the intermediate heat and salt content were exported to the deep layers from 2009 to 2013 thanks to deep convection processes. In 2014, a LIW step of +0.3 [Formula: see text]C and +0.08 in salinity could be observed concomitant with a weak winter convection. Warmer and more saline LIW subsequently accumulated in the northwestern basin in the absence of intense deep convective winters until 2018. Deep stratification below the LIW thus increased, which, together with the air-sea heat fluxes intensity, constrained the depth of convection. A key prognostic indicator of the intensity of deep convective events appears to be the convection depth of the previous year.

The international array of profiling floats known as Argo is a major component of the global ocean-and climate-observing system. In 2010, the NAOS (Novel Argo Observing System) project was selected as part of France's Equipex "Investissement d'Avenir" program. The objectives of NAOS were to consolidate the French contribution to the Argo core mission (global temperature and salinity measurements down to 2,000 m) as well as to develop the future generation of French Argo profiling floats and prepare the next phase of the Argo program with an extension to the deep ocean (Deep-Argo), biogeochemistry (BGC-Argo) and polar seas. This paper summarizes the main technological advances and at-sea validations carried out as part of NAOS: development of a deep (4,000 m) float, a new BGC float for Research & Development (R&D) applications, and a BGC float for deployments in Arctic areas, assessment of a new density and Absolute Salinity optical sensor, improvement of the reliability of the standard Argo float, and upgraded satellite-transmission performance. French profiling floats developed in this way are now operational and among the most deployed worldwide, and the density sensor is the most promising of its kind for profiling floats applications.

Argo, the international array of profiling floats, is a major component of the global ocean and climate observing system. In 2010, the NAOS (Novel Argo Observing System) project was selected as part of the French "Investissements d'Avenir" Equipex program. The objectives of NAOS were to consolidate the French contribution to Argo's core mission (global temperature and salinity measurements down to 2000 m), and also to develop the future generation of French Argo profiling floats and prepare the next phase of the Argo program with an extension to the deep ocean (Deep Argo), biogeochemistry (BGC-Argo) and polar seas. This paper summarizes how NAOS has met its objectives. The project significantly boosted France's contribution to Argo's core mission by deploying more than 100 NAOS standard Argo profiling floats. In addition, NAOS deployed new-generation floats as part of three scientific experiments: biogeochemical floats in the Mediterranean Sea, biogeochemical floats in the Arctic Ocean, and deep floats with oxygen sensors in the North Atlantic. The experiment in the Mediterranean Sea, launched in 2012, implemented and maintained a network of BGC-Argo floats at basin scale for the first time. The 32 BGC-Argo floats deployed and about 4000 BGC profiles collected have vastly improved characterization of the biogeochemical and ecosystem dynamics of the Mediterranean. Meanwhile, experiments in the Arctic and in the North Atlantic, starting in 2015 and deploying 20 Arctic BGC floats and 23 deep floats, have provided unique observations on biogeochemical cycles in the Arctic and deep-water masses, as well as ocean circulation variability in the North Atlantic. NAOS has therefore paved the way to the new operational phase of the Argo program in France that includes BGC and Deep Argo extensions. The objectives and characteristics of this new phase of Argo-France are discussed in the conclusion.

The necessity of wide, global-scale observing systems for marine biogeochemistry emerged dramatically in the last decade. A global network based on Biogeochemical (BGC) Argo floats is considered to be one of the most promising approaches for reaching this goal. As a first step, pilot studies were encouraged to test the feasibility of a global BGC-Argo array, to consolidate the methods and practices under development, and to set up the array's characteristics. A pilot study in The Mediterranean Sea-deemed a suitable candidate for a test case because it combines a relatively large diversity of oceanic BGC conditions in a reduced open-ocean basin-was consequently approved as a part of the "Novel Argo ocean Observing System" (NAOS) project, a French national initiative to promote, consolidate, and develop the Argo network. We present here a first assessment of the NAOS Mediterranean array, in view of scientific choices on observing-system strategy, on implementation and statistics on network performances, and on data-quality control.

In the western Mediterranean Sea, Levantine intermediate waters (LIW), which circulate below the surface productive zone, progressively accumulate nutrients along their pathway from the Tyrrhenian Sea to the Algerian Basin. This study addresses the role played by diffusion in the nutrient enrichment of the LIW, a process particularly relevant inside step-layer structures extending down to deep waters-structures known as thermohaline staircases. Profiling float observations confirmed that staircases develop over epicentral regions confined in large-scale circulation features and maintained by saltier LIW inflows on the periphery. Thanks to a high profiling frequency over the 4-year period 2013-2017, float observations reveal the temporal continuity of the lay-ering patterns encountered during the cruise PEACETIME and document the evolution of layer properties by about + 0.06°C in temperature and +0.02 in salinity. In the Al-gerian Basin, the analysis of in situ lateral density ratios untangled double-diffusive convection as a driver of thermoha-line changes inside epicentral regions and isopycnal diffusion as a driver of heat and salt exchanges with the surrounding sources. In the Tyrrhenian Sea, the nitrate flux across ther-mohaline staircases, as opposed to the downward salt flux, contributes up to 25 % of the total nitrate pool supplied to the LIW by vertical transfer. Overall, however, the nutrient enrichment of the LIW is driven mostly by other sources, coastal or atmospheric, as well as by inputs advected from the Algerian Basin.

The societal benefits of satellite ocean colour include aiding the management of the marine ecosystem, helping understand the role of the ocean ecosystem in climate change, aquaculture, fisheries, coastal zone water quality, and the mapping and monitoring of harmful algal blooms. Ocean colour is also designated as an essential climate variable by the Global Climate Observing System (GCOS). However, in order to have confidence in earth observation data, measurements made at the surface of the Earth, with the intention of providing verification or validation of satellite mounted sensor measurements, should be trustworthy and of the same high quality as those taken with the satellite sensors themselves. In order to be trustworthy, in situ validation measurements should include an unbroken chain of SI traceable calibrations and comparisons and full uncertainty budgets for each of the in situ sensors used. This metrological traceability is beginning to be demanded by the space agencies for satellite validation measurements and, for ocean colour, should follow the guidelines and protocols of the ESA Fiducial Reference Measurements for Satellite Ocean Colour (FRM4SOC) project (www.frm4soc.org). Until now, this has not been the case for most measurements used for validation, including those taken in the Aegean and Eastern Mediterranean. Subsequently, the Hellenic Centre for Marine Research (HCMR), in cooperation with the Laboratory of Optical Metrology (LOM), has started to follow the FRM direction by ensuring that the radiometers of its optical suite underwent SI-traceable absolute radiometric calibration. This included an estimate of the radiometry calibration uncertainty budget and was performed at the marine optical laboratory of the European Commission’s Joint Research Centre prior to their deployment on the recent PERLE-2 oceanographic cruise in the Eastern Mediterranean (Feb-Mar 2019). As well as irradiance and radiance sensors, the HCMR optical suite also houses instruments for measuring inherent optical properties (IOP) of the water column. Therefore, this paper presents the in-water radiometry matchups from PERLE-2 with Sentinel-3 Ocean and Land Colour Instrument (OLCI) measurements, and investigates their validation potential. It also presents the PERLE-2 cruise profile chlorophyll and backscatter measurements that aid this effort through characterizing the light scattering and absorbing constituents that contribute to the signal detected by satellite ocean colour sensors during validation matchups.

As commonly observed in oligotrophic stratified waters, a subsurface (or deep) chlorophyll maximum (SCM) frequently characterizes the vertical distribution of phyto-plankton chlorophyll in the Mediterranean Sea. Occurring far from the surface layer "seen" by ocean colour satellites , SCMs are difficult to observe with adequate spatio-temporal resolution and their biogeochemical impact remains unknown. Biogeochemical-Argo (BGC-Argo) profiling floats represent appropriate tools for studying the dynamics of SCMs. Based on data collected from 36 BGC-Argo floats deployed in the Mediterranean Sea, our study aims to address two main questions. (1) What are the different types of SCMs in the Mediterranean Sea? (2) Which environmental factors control their occurrence and dynamics? First, we analysed the seasonal and regional variations in the chlorophyll concentration (Chl a), particulate backscattering coefficient (b bp), a proxy of the particulate organic carbon (POC) and environmental parameters (photosynthetically active radiation and nitrates) within the SCM layer over the Mediter-ranean Basin. The vertical profiles of Chl a and b bp were then statistically classified and the seasonal occurrence of each of the different types of SCMs quantified. Finally, a case study was performed on two contrasted regions and the environmental conditions at depth were further investigated to understand the main controls on the SCMs. In the eastern basin, SCMs result, at a first order, from a photoacclima-tion process. Conversely, SCMs in the western basin reflect a biomass increase at depth benefiting from both light and nitrate resources. Our results also suggest that a variety of intermediate types of SCMs are encountered between these two endmember situations.

An accurate understanding of the biogeochemistry of dissolved phosphate pool in the upper waters of P-depleted oceanic regions is constrained by the low sensitivity of routine phosphate measurements. In this study, by using the sensitive Liquid Waveguide Capillary Cell method, we report the first extensive cross-basin survey of nanomolar dissolved inorganic phosphate (DIP) and dissolved organic phosphate (DOP) concentration in P-depleted surface waters of the Mediterranean Sea during the stratification period. In the north western Mediterranean Sea (NWMS), DIP above the mixed layer depth (MLD) ranged between 4.9 and 26.5 nM. Along an E-W transect crossing Ionian and Tyrrhenian Seas (E-W transect), DIP above the MLD was lower, ranging between 0.9 and 11.4 nM. Contrarily to the traditional view of a depleted and invariant surface dissolved phosphate pool, a significant vertical variability of DIP and DOP was revealed in upper waters. A positive gradient of DIP was observed above the phosphacline, between the MLD and the deep chlorophyll maximum (DCM) depth, suggesting a potential diffusion of new phosphate to near-surface waters, even under stratified conditions. Interestingly, despite this apparent DIP availability, a significant negative gradient of DOP concentration was observed in the same layer. Finally, the positive gradient in DIP coincided with a significant increase in N:P ratio, suggesting a higher rate of increase of N than of P. The results obtained in this study indicate that acquiring nanomolar DIP data is a sine qua non condition for the comprehension and prediction of the biogeochemical functioning of P-depleted oceanic regions, such as the Mediterranean Sea.

We report on data from an oceanographic cruise, covering western, central and eastern parts of the Mediterranean Sea, on the French research vessel Tethys 2 in May 2015. This cruise was fully dedicated to the maintenance and the metrological verification of a biogeochemical observing system based on a fleet of BGC-Argo floats. During the cruise, a comprehensive data set of parameters sensed by the autonomous network was collected. The measurements include ocean currents, seawater salinity and temperature, and concentrations of inorganic nutrients, dissolved oxygen and chlorophyll pigments. The analytical protocols and data processing methods are detailed, together with a first assessment of the calibration state for all the sensors deployed during the cruise.

Dissolved oxygen (O2) is a relevant tracer to interpret variations of both water mass properties in the open ocean and biological production in the surface layer of both coastal and open waters. Deep-water formation is very active in the northwestern Mediterranean Sea, where it influences intermediate and deep waters properties, nutrients replenishment and biological production. This study analyses, for the first time, the 20-year time series of monthly O2 concentrations at the DYFAMED long-term sampling site in the Ligurian Sea. Until the winters of 2005 and 2006, a thick and strong oxygen minimum layer was present between 200 and 1300 m because dense water formation was then local, episodic and of low intensity. In 2005-2006, intense and rapid deep convection injected 24 mol O2 m-2 between 350 and 2000 m from December 2005 to March 2006. Since this event, the deep layer has been mostly ventilated during winter time by newly formed deep water spreading from the Gulf of Lion 250 km to the west and by some local deep mixing in early 2010, 2012 and 2013. In the context of climate change, it is predicted that the intensity of deep convection will become weaker in the Mediterranean, which could potentially lead to hypoxia in intermediate and deep layers with substantial impact on marine ecosystems. With the exception of winters 2005 and 2006, the O2 changes in surface waters followed a seasonal trend that reflected the balance between air-sea O2 exchanges, changes in the depth of the mixed layer and phytoplankton net photosynthesis. We used the 20-year O2 time series to estimate monthly and annual net community production. The latter was 7.1 mol C m-2 yr-1, consistent with C-14 primary production determinations and sediment-trap carbon export fluxes at DYFAMED.

During winter 2012–2013, open‐ocean deep convection which is a major driver for the thermohaline circulation and ventilation of the ocean, occurred in the Gulf of Lions (Northwestern Mediterranean Sea) and has been thoroughly documented thanks in particular to the deployment of several gliders, Argo profiling floats, several dedicated ship cruises, and a mooring array during a period of about a year. Thanks to these intense observational efforts, we show that deep convection reached the bottom in winter early in February 2013 in a area of maximum 28 ± 3 10⁹ m². We present new quantitative results with estimates of heat and salt content at the subbasin scale at different time scales (on the seasonal scale to a 10 days basis) through optimal interpolation techniques, and robust estimates of the deep water formation rate of 2.0 ± 0.2 Sv. We provide an overview of the spatiotemporal coverage that has been reached throughout the seasons this year and we highlight some results based on data analysis and numerical modeling that are presented in this special issue. They concern key circulation features for the deep convection and the subsequent bloom such as Submesoscale Coherent Vortices (SCVs), the plumes, and symmetric instability at the edge of the deep convection area.

A multi-element biogeochemical model forced by a 1 km resolution hydrodynamical model was used to gain in understanding of the biogeochemical functioning of the North-Western Mediterranean (NW Med), the only region in the whole Mediterranean Sea with a marked and recurrent spring bloom behavior related to the winter dense water formation characterizing this area. After an assessment of the simulation using satellite derived chlorophyll and Dewex project in situ nutrients observations, the nitrogen and phosphorus seasonal cycles were analyzed using model outputs on the period 2012-2013. Injections of nutrients during the wind intensification period allow the triggering of the autumn bloom. Then, convection in winter upwells large amounts of nutrients in the euphotic layer. When the conditions for phytoplankton development are gathered (reduction of vertical mixing, low grazing pressure), a bloom is triggered with a massive consumption of nutrients during more than one month resulting at the end of April in a depletion of nutrients at the surface. Nutrients consumption continues to deplete nutrients at increasing depth, increasing the nutriclines and deep chlorophyll maximum depths. That finally leads to the summer oligotrophy of the water column. Then a quantification of nitrogen and phosphorus budgets of the open-sea convection area was performed on an annual basis. The deep convection area represents a sink of nitrate and phosphate, and a source of organic nitrogen and phosphorus for the peripheric regions. Regarding the biogeochemical nitrogen cycle, the deep-nitrate based new production is responsible for 19% of the total nitrogen uptake. This new production dominates during the winter deep convection and spring bloom periods. Finally, our results suggest that the NW Med open sea convection represents a major source of nutrients for the Mediterranean surface sea.

This study was a part of the DeWEX project (Deep Water formation Experiment), designed to better understand the impact of dense water formation on the marine biogeochemical cycles. Here, nutrient and phytoplankton vertical and horizontal distributions were investigated during a deep open‐ocean convection event and during the following spring bloom in the Northwestern Mediterranean Sea (NWM). In February 2013, the deep convection event established a surface nutrient gradient from the center of the deep convection patch to the surrounding mixed and stratified areas. In the center of the convection area, a slight but significant difference of nitrate, phosphate and silicate concentrations was observed possibly due to the different volume of deep waters included in the mixing or to the sediment resuspension occurring where the mixing reached the bottom. One of this process, or a combination of both, enriched the water column in silicate and phosphate, and altered significantly the stoichiometry in the center of the deep convection area. This alteration favored the local development of microphytoplankton in spring, while nanophytoplankton dominated neighboring locations where the convection reached the deep layer but not the bottom. This study shows that the convection process influences both winter nutrients distribution and spring phytoplankton distribution and community structure. Modifications of the convection's spatial scale and intensity (i.e., convective mixing depth) are likely to have strong consequences on phytoplankton community structure and distribution in the NWM, and thus on the marine food web.

Knowledge of the relative contributions of phytoplankton size classes to zooplankton biomass is necessary to understand food-web functioning and response to climate change. During the Deep Water formation Experiment (DEWEX), conducted in the northwest Mediterranean Sea in winter (February) and spring (April) of 2013, we investigated phytoplankton-zooplankton trophic links in contrasting oligotrophic and eutrophic conditions. Size fractionated particulate matter (pico-POM, nano-POM, and micro-POM) and zooplankton (64 to >4000 lm) composition and carbon and nitrogen stable isotope ratios were measured inside and outside the nutrient-rich deep convection zone in the central Liguro-Provencal basin. In winter, phytoplankton biomass was low (0.28 mg m 23) and evenly spread among picophytoplankton, nanophyto-plankton, and microphytoplankton. Using an isotope mixing model, we estimated average contributions to zooplankton biomass by pico-POM, nano-POM, and micro-POM of 28, 59, and 15%, respectively. In spring, the nutrient poor region outside the convection zone had low phytoplankton biomass (0.58 mg m 23) and was dominated by pico/nanophytoplankton. Estimated average contributions to zooplankton biomass by pico-POM, nano-POM, and micro-POM were 64, 28 and 10%, respectively, although the model did not differentiate well between pico-POM and nano-POM in this region. In the deep convection zone, spring phyto-plankton biomass was high (1.34 mg m 23) and dominated by micro/nano phytoplankton. Estimated average contributions to zooplankton biomass by pico-POM, nano-POM, and micro-POM were 42, 42, and 20%, respectively, indicating that a large part of the microphytoplankton biomass may have remained ungrazed. Plain Language Summary The grazing of zooplankton on algal phytoplankton is a critical step in the transfer of energy through all ocean food webs. Although microscopic, phytoplankton span an enormous size range. The smallest picophytoplankton are generally thought to be too small to be directly grazed by zooplankton, resulting in less efficient energy transfer through the food web. This has implications for our future oceans where warming and lower nutrient supply are predicted to favor picophytoplankton over the larger nanosize and microsize classes. We tested the importance of phytoplankton size classes in the transfer of energy to zooplankton in the northwest Mediterranean Sea, where conditions naturally result in contrasting regions of small and large phytoplankton dominance. Contrary to expectation, biochemical tracers showed that microphytoplankton never contributed more than 20% to zooplankton biomass, even in regions where microphytoplankton were plentiful. On the other hand, picophytoplankton contributed 25-65% to zooplankton biomass. This finding indicates that there are well-established food-web pathways from picophytoplankton to zooplankton, and that these pathways play an important role even in ocean regions where microphytoplankton dominate. Accordingly, a decline in phytoplankton size classes may have a greater effect on carbon sequestration than on food-web productivity.

During the winter 2013, an intense observation and monitoring was performed in the north-western Mediterranean Sea to study deep water formation process that drives thermohaline circulation and biogeochemical processes (HYMEX SOP2 and DEWEX projects). To observe intensively and continuously the impact of deep convection on oxygen (O₂) ventilation, an observation strategy was based on the enhancement of the Argo-O₂ floats to monitor the offshore dense water formation area (DWF) in the Gulf of Lion prior to and at the end of the convective period (December 2012 to April 2013). The intense O₂ measurements performed through shipborne CTD casts and Argo-O₂ floats deployment revealed an O₂ inventory rapidly impacted by mixed layer (ML) deepening on the month scale. The open-sea convection in winter 2013 ventilated the deep waters from mid-February to the end of May 2013. The newly ventilated dense water volume, based on an Apparent Oxygen Utilization (AOU) threshold, was estimated to be about 1.5 × 10¹³ m³ during the DWF episode, increasing the deep O2₂ concentrations from 196 to 205 µmol kg⁻¹ in the north-western basin.

The North Western Mediterranean Sea exhibits recurrent and significant autumnal and spring phytoplankton blooms. The existence of these two blooms coincides with typical temperate dynamics. To determine the potential control of physical and biogeochemical factors on these phytoplankton blooms, data from a multiplatform approach (combining ships, Argo and BGC‐Argo floats, and bio‐optical gliders) were analyzed in association with satellite observations in 2012–2013. The satellite framework allowed a simultaneous analysis over the whole annual cycle of in situ observations of mixed layer depth, photosynthetical available radiation, particle backscattering, nutrients (nitrate and silicate), and chlorophyll‐a concentrations. During the year 2012–2013, satellite ocean color observations, confirmed by in situ data, have revealed the existence of two areas (or bioregions) with comparable autumnal blooms but contrasting spring blooms. In both bioregions, the ratio of the euphotic zone (defined as the isolume 0.415 mol photons m−2 d−1, Z0.415) and the MLD identified the initiation of the autumnal bloom, as well as the maximal annual increase in [Chl‐a] in spring. In fact, the autumnal phytoplankton bloom might be initiated by mixing of the summer shallowing deep chlorophyll maximum, while the spring restratification (when Z0.415/MLD ratio became >1) might induce surface phytoplankton production that largely overcomes the losses. Finally, winter deep convection events that took place in one of the bioregions induced higher net accumulation rate of phytoplankton in spring associated with a diatom‐dominated phytoplankton community principally. We suggest that very deep winter MLD lead to an increase in surface silicates availability, which favored the development of diatoms.

The Northwestern Mediterranean (NWMed) Sea includes one of the best observed ocean deep convection sites in the World. An observing system simulation experiment (OSSE) is developed to provide a methodology for estimating observing network errors. It is applied to quantify dense water volumes in the NWMed during 2012-2013 with their observation error from MOOSE network. Results from the OSSE show low spatiotemporal sampling errors, which confirms MOOSE network ability to measure dense waters. However, results are highly sensitive to instrumental stability. The dense water volume is then estimated in observations from four ship cruises between summers 2012 and 2013. A large seasonal cycle is found, maximal in spring 2013 and dominated by the area west of 6.58E. The dense water volume ( \sigma₀ 29.11 kg/m³) is stable between summer 2012 (13.3 +/- 0.6 x 10¹³ m³) and winter 2013 (13.7 +/- 1.3 x 10¹³ m³). It increases dramatically in spring 2013 (17.7 +/- 0.9 x 10¹³ m³) due to an intense convective event, and it finally decreases rapidly in summer 2013 (15.1 +/- 0.6 x 10¹³ m³) due to restratification and spreading. We estimate an open-sea dense water formation (DWF) rate of 1.4 +/- 0.3 Sv between summer 2012 and spring 2013 over the studied area, extrapolated to 2.3 +/- 0.5 Sv over the whole NWMed Sea and for the optimal timing. This is to our knowledge the highest measured DWF rate, suggesting winter 2013 was exceptionally convective. The observed restratification rate between spring and summer 2013 is -0.8 +/-0.4 Sv. This study provides robust quantifications of deep convection during an exceptional event that will allow to evaluate numerical simulations.

Since 2010, an intense effort in the collection of in situ observations has been carried out in the northwestern Mediterranean Sea thanks to gliders, profiling floats, regular cruises, and mooring lines. This integrated observing system enabled a year-to-year monitoring of the deep waters formation that occurred in the Gulf of Lions area during four consecutive winters (2010-2013). Vortical structures remnant of wintertime deep vertical mixing events were regularly sampled by the different observing platforms. These are Submesoscale Coherent Vortices (SCVs) characterized by a small radius ($5-8 km), strong depth-intensified orbital velocities ($10-20 cm s 21) with often a weak surface signature, high Rossby ($0.5) and Burger numbers O(0.5-1). Anticyclones transport convected waters resulting from intermediate ($300 m) to deep ($2000 m) vertical mixing. Cyclones are characterized by a 500-1000 m thick layer of weakly stratified deep waters (or bottom waters that cascaded from the shelf of the Gulf of Lions in 2012) extending down to the bottom of the ocean at $2500 m. The formation of cyclonic eddies seems to be favored by bottomreaching convection occurring during the study period or cascading events reaching the abyssal plain. We confirm the prominent role of anticyclonic SCVs and shed light on the important role of cyclonic SCVs in the spreading of a significant amount ($30%) of the newly formed deep waters away from the winter mixing areas. Since they can survive until the following winter, they can potentially have a great impact on the mixed layer deepening through a local preconditioning effect.

In 2013, as part of the French NAOS (Novel Argo Oceanic observing System) program, five profiling floats equipped with nitrate sensors (SUNA-V2) together with CTD and bio-optical sensors were deployed in the Mediterranean Sea. At present day, more than 500 profiles of physical and biological parameters were acquired, and significantly increased the number of available nitrate data in the Mediterra-nean Sea. Results obtained from floats confirm the general view of the basin, and the well-known west-to-east gradient of oligotrophy. At seasonal scale, the north western Mediterranean displays a clear temperate pattern sustained by both deep winter mixed layer and shallow nitracline. The other sampled areas follow a subtropical regime (nitracline depth and mixed layer depth are generally decoupled). Float data also permit to highlight the major contribution of high-frequency processes in controlling the nitrate supply during winter in the north western Mediterranean Sea and in altering the nitrate stock in subsurface in the eastern basin.

Since 2007, gliders have been regularly deployed in the northwestern Mediterranean Sea, a crucial region regarding the thermohaline circulation of the Mediterranean Sea. It revealed for the first time very warm (10.48C) and saline (10.1) submesoscale anticyclones at intermediate depth characterized by a small radius ($5 km), high Rossby ($0.3), and Burger ($0.7) numbers. They are likely order of 10 to be formed each year, have a life time order a year and certainly contribute significantly to the spreading of the Levantine Intermediate Waters (LIW) toward the whole subbasin, thus potentially impacting wintertime vertical mixing through hydrographical and dynamical preconditioning. They could be mainly formed by the combined action of turbulent mixing and flow detachment of the northward flow of LIW at the northwestern headland of Sardinia. Upwelling conditions along the western coast of Sardinia associated with a southward geostrophic flow within the upper layers seem to play a key role in their formation process.

The present data publication provides permanent links to original and updated versions of validated data files. The data files include properties of seawater, particulate matter and dissolved matter that were measured from discrete water samples collected with Niskin bottles during the 2009-2013 Tara Oceans expedition. Properties include pigment concentrations from HPLC analysis (10 depths per vertical profile, 25 pigments per depth), the carbonate system (Surface and 400m; pH (total scale), CO₂, pCO₂, <em>f</em>CO₂, HCO₃, CO₃, Total alkalinity, Total carbon, OmegaAragonite, OmegaCalcite, and dosage Flags), nutrients (10 depths per vertical profile; NO₂, PO₄, NO₂/NO₃, SI, quality Flags), DOC, CDOM, and dissolved oxygen isotopes. The Service National d'Analyse des Paramètres Océaniques du CO₂, at the Université Pierre et Marie Curie, determined CT and AT potentiometrically (Edmond 1970; DOE 1994) on samples preserved according to Dickson et al. (2007). More than 250 vertical profiles of these properties were made across the world ocean. DOC, CDOM and dissolved oxygen isotopes are available only for the Arctic Ocean and Arctic Seas (2013).

In the last 5 years, an unprecedented effort in the sampling of the Northern Current (NC) has been carried out using gliders which collected more than 50 000 profiles down to 1000m maximum along a few repeated sections perpendicular to the French coast. Based on this dataset, this study presents a very first quantitative picture of the NC on 0-1000m depth. We show its mean structure of temperature and salinity characterized by the different Water Masses of the basin (Atlantic Water, Winter Intermediate Water, Levantine Intermediate Water and Western Mediterranean Deep Water) for each season and at different location. Geostrophic currents are derived from the integration of the thermal-wind balance using the mean glider-estimate of the current during each dive as a reference. Estimates of the heat, salt, and volume transport are then computed in order to draw an heat and salt budget of the NC. The results show a strong seasonal variability due to the intense surface buoyancy loss in winter resulting in a vertical mixing offshore that makes the mixed layer depth reaching several hundreds of meters in the whole basin and in a very particular area down to the bottom of the sea-floor (deep convection area). The horizontal density gradient intensifies in winter leading to geostrophic currents that are more intense and more confined to the continental slope, and thus to the enhancement of the mesoscale activity (meandering, formation of eddies through baroclinic instability...). The mean transport estimates of the NC is found to be about 2-3Sv greater than previous spurious estimates. The heat budget of the NC also provides an estimate of the mean across shore heat/salt flux directly impacting the region in the Gulf of Lion where deep ocean convection, a key process in the thermohaline circulation of the Mediterranean Sea, can occur in Winter.

Investigating transport pathways in the ocean

Since 2008, gliders repeated transects crossing the basin of the North Western Mediterranean Sea and regularly sampled mesoscale structures with an high horizontal resolution of about 2-3 km between each profile required in that region of small internal deformation radius (<10km). By analysing more than 50 000 profiles collected by these gliders in the last 5 years, we were able to identify several types of eddies regarding the water mass composing their cores: Winter Intermediate Water (WIW), Levantine Intermediate Water (LIW), Western Mediterranean Deep Water (WMDW). Most of them are anticyclonic structures with Rossby Number greater than 0.1 and tend to be characterized by a core in the inner ocean. Some of them whose formation has been dated several months back in time can be qualified as long lived features. Of particular interest to assess the role of mesoscale eddies in the ocean circulation, a Submesoscale Coherent Vortex (SCV) composed of newly WMDW was observed nine months after its formation. We also used a 1 year run of a high resolution (1km, 40 vertical levels) numerical model of the region (SYMPHONIE) which is able to reproduce similar eddies. In this study we discuss their formation process (instability of the boundary current, or diapycnal mixing followed by geostrophic adjustment) based on comparisons between these observations and the model outputs, and try to estimate their impact on the general circulation of this basin.

From 2008 on, repeated sections crossing the Northern Current (NC) were operated by gliders as part of a global observing system (MOOSE project) of the North Western Mediterranean Sea. This work is dedicated to the analysis of the submesoscale thermohaline variability at the margin of this current observed by gliders. The mean circulation of the basin is characterized by a cyclonic gyre (whose Northern part is the so-called NC) associated with a doming of the isopycnals preconditionning the whole interior basin to great vertical mixing. The thermal and haline differences between the Atlantic Water (AW) transported by the NC and older and modified AW off the coast leads to a frontal structure. Especially in winter, when the mixed layer depth used to reach several hundreds of meters offshore, isopycnal outcropping and the role of frontal processes are enhanced leading to intense variability at scales smaller than the deformation radius. Based on diagnostics using the Potential Vorticity (PV) computed from the glider data assuming quasi-geostrophic conditions and no variation in the alongshore direction, we discuss the dynamical processes at work, with a focus on 2 typical examples: (1) the first example takes place in winter during a strong vertical mixing event. While the glider crossed the frontal region, the temperature and salinity fields exhibit vertical motions at depths about 0-400m. Frontogenesis might be at play through mesoscale strain since the glider shows an intense mesoscale activity but a weak stratification and enhanced horizontal buoyancy gradient actually make the Ertel PV reach negative values and symmetric instability is likely to be a prominent mechanism explaining the observed variability. (2) the second example takes place in spring. We identify an episode of down-front wind blowing during the glider deployment which could have extracted PV from the surface layer. However, the geostrophic turbulence is in that case likely to play a key role in the formation of the observed variability of the temperature and salinity since it is organized along slopes characterized by an aspect ratio of an order of f/N.

The Mediterranean Forecasting System (MFS) has been operational for a decade, and is continuously providing forecasts and analyses for the region. These forecasts comprise local- and basin-scale information of the environmental state of the sea and can be useful for tracking oil spills and supporting search-and-rescue missions. Data assimilation is a widely used method to improve the forecast skill of operational models and, in this study, the three-dimensional variational (OceanVar) scheme has been extended to include Argo float trajectories, with the objective of constraining and ameliorating the numerical output primarily in terms of the intermediate velocity fields at 350 m depth. When adding new datasets, it is furthermore crucial to ensure that the extended OceanVar scheme does not decrease the performance of the assimilation of other observations, e.g., sea-level anomalies, temperature, and salinity. Numerical experiments were undertaken for a 3-year period (2005-2007), and it was concluded that the Argo float trajectory assimilation improves the quality of the forecasted trajectories with similar to 15%, thus, increasing the realism of the model. Furthermore, the MFS proved to maintain the forecast quality of the sea-surface height and mass fields after the extended assimilation scheme had been introduced. A comparison between the modeled velocity fields and independent surface drifter observations suggested that assimilating trajectories at intermediate depth could yield improved forecasts of the upper ocean currents.

The semi-enclosed nature of the Mediterranean Sea, together with its smaller inertia due to the relative short residence time of its water masses, make it highly reactive to external forcings, in particular variations of water, energy and matter fluxes at the interfaces. This region, which has been identified as a “hotspot” for climate change, is therefore expected to experience environmental impacts that are considerably greater than those in many other places around the world. These natural pressures interact with the increasing demographic and economic developments occurring heterogeneously in the coastal zone, making the Mediterranean even more sensitive. This review paper aims to provide a review of the state of current functioning and responses of Mediterranean marine biogeochemical cycles and ecosystems with respect to key natural and anthropogenic drivers and to consider the ecosystems’ responses to likely changes in physical, chemical and socio-economical forcings induced by global change and by growing anthropogenic pressure at the regional scale. The current knowledge on and expected changes due to single forcing (hydrodynamics, solar radiation, temperature and acidification, chemical contaminants) and combined forcing (nutrient sources and stoichiometry, extreme events) affecting the biogeochemical fluxes and ecosystem functioning are explored. Expected changes in biodiversity resulting from the combined action of the different forcings are proposed. Finally, modeling capabilities and necessity for modeling are presented. A synthesis of our current knowledge of expected changes is proposed, highlighting relevant questions for the future of the Mediterranean ecosystems that are current research priorities for the scientific community. Finally, we discuss how these priorities can be approached by national and international multi-disciplinary research, which should be implemented on several levels, including observational studies and modeling at different temporal and spatial scales.

The impact of Argo float trajectory assimilation on the quality of ocean analyses is studied by means of an operational oceanographic model implemented in the Mediterranean Sea and a 3D-Var assimilation scheme. For the first time, both Argo trajectories and vertical profiles of temperature and salinity (TS) together with satellite altimeter data of sea level anomaly (SLA) are assimilated to produce analyses for short-term forecasts. The study period covers 3 months during winter 2005 when four Argo trajectories were present in the northwestern Mediterranean Sea. The scheme is first assessed computing the misfits between observations and model forecast and analysis. The misfit statistics appear improved for float trajectories, while they are not degraded for the other assimilated variables (TS profiles and SLA). This indicates that the trajectory integration is consistent with the other components of the assimilation system and provides new information on horizontal pressure gradients. Comparisons between analyses obtained with and without trajectory assimilation suggest that trajectory assimilation can have an impact on the description of boundary currents and their instabilities, as well as mesoscale activity at regional scales. Changes are depicted by intermediate water mass redistributions, mesoscale eddy relocations, and net transport modulations. These impacts are detailed and assessed considering historical and simultaneous in situ data sets. The results motivate the integration of Argo trajectories in the operational Mediterranean Forecasting System.

Not Available

Not Available

In this paper we present an application of a variational method for the reconstruction of the velocity field in a coastal flow in the central Adriatic Sea, using in situ data from surface drifters and outputs from the ROMS circulation model. The variational approach, previously developed and tested for mesoscale open ocean flows, has been improved and adapted to account for inhomogeneities on boundary current dynamics over complex bathymetry and coastline and for weak Lagrangian persistence in coastal flows. The velocity reconstruction is performed using nine drifter trajectories over 45 d, and a hierarchy of indirect tests is introduced to evaluate the results as the real ocean state is not known. For internal consistency and impact of the analysis, three diagnostics characterizing the particle prediction and transport, in terms of residence times in various zones and export rates from the boundary current toward the interior, show that the reconstruction is quite effective. A qualitative comparison with sea color data from the MODIS satellite images shows that the reconstruction significantly improves the description of the boundary current with respect to the ROMS model first guess, capturing its main features and its exchanges with the interior when sampled by the drifters.

La dynamique de la circulation marine en région côtière peut être appréhendée par des moyens d'observation et de modélisation. Si le domaine côtier bénéficie de réseaux de mesure adaptés en résolution et en couverture, la mise en oeuvre d'un outil de modélisation s'avère délicate, particulièrement dans la spécification de conditions initiales et de conditions aux frontières ouvertes appropriées. Cette question de l'estimation des conditions aux limites peut être considérée comme relevant d'une problématique plus générale de détermination des paramètres de fonctionnement d'un modèle à partir de l'information apportée par les observations. Ainsi posé, le problème est abordé par l'outil mathématique du contrôle optimal qui conduit à une opération d'assimilation de données par méthode d'adjoint. L'adaptation et l'application de cette formulation inverse à un problème de modélisation océanographique contrôlé par ses conditions frontières font l'objet de ce travail de thèse. Une approche 4D-variationnelle d'assimilation pour le contrôle des conditions frontières est ainsi développée et implémentée sur la base du modèle aux équations primitives OPA (version 'toit rigide') et de son adjoint existant. L'algorithme d'inversion est validé sur trois expériences d'identification: (i) les perturbations du forçage barotrope sur la frontière ouverte d'un modèle non stratifié sont détectées par l'assimilation de courants de surface pris à l'intérieur du domaine; (ii) le développement d'un courant de densité dans un canal semi infini et homogène est retrouvé à partir de profils hydrologiques distribués en aval; (iii) l'observation du " spin up " d'une circulation de plateau induit par un courant côtier permet d'identifier les conditions amonts prévues par le mécanisme de séparation d'un front de Kelvin barocline en incidence sur un talus. Cette méthode d'assimilation est alors proposée pour des applications s'inscrivant dans la problématique scientifique des échanges côte - large sur la zone du Golfe du Lion. Riche en observations et largement ouverte sur le bassin méditerranéen occidental, cette région est soumise à l'influence du courant Liguro-Provençal qui s'écoule en bord de plateau. L'algorithme d'inversion permet de contrôler cette influence régionale sur la circulation de plateau; il est donc évalué à travers des expériences réalistes menées dans le mode diagnostique. En particulier, la faculté de réseaux d'observation à détecter des épisodes d'intrusion du courant côtier sur le golfe est traitée suivant leur couverture spatiale et temporelle.