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

Ocean primary production is a key process that regulates marine ecosystems and the global climate, but its estimation is still affected by multiple uncertainties. Typically, the chlorophyll-a concentration (CHL) is used to characterise this process, as it is considered as a proxy of phytoplankton biomass. To date, the most common observing systems for studying CHL are ocean colour satellites and Biogeochemical-Argo (BGC-Argo) floats. These are complementary systems: satellite observations provide global coverage but are limited to the ocean surface, while BGC-Argo floats provide punctual observations along the whole water column. Quantitative matching of these two observing systems has been obtained only at regional or single-float scales, while at a global scale the relatively low and irregular BGC-Argo coverage results in large uncertainties. Here, we propose a different method, by comparing satellite and BGC-Argo climatological annual time series within seven different bioregions, each characterised by a homogeneous phytoplankton phenology, allowing us to smooth the uncertainties. By comparing the mean values, amplitudes, and shapes of the two time series, we identify regions (a) where they agree (58%-61% of the ocean surface area); (b) regions undersampled by the BGC-Argo float network (particularly in the Arabian Sea and near the Amazon delta); (c) where the discrepancy may stem from satellite or (d) BGC-Argo performance (mainly found at subtropical and high latitudes, respectively). Caution is required when using BGC-Argo and satellite data in regions b-d, and, for each region, we provide suggestions on which system could be affected by the largest uncertainties.

<div><p>Inversion models, in the context of oceanography, relate the observed ocean color to the concentrations of the different biogeochemical components present in the water of the ocean. However, building accurate inversion models can be quite complex due to the many factors that can influence the observed ocean color, such as variations in the composition or the optical properties of biogeochemical products. Here we assess the feasibility of the inversion approach, by implementing the three-stream light inversion model in a one-dimensional water column configuration, represented at the BOUSSOLE site in the northwestern Mediterranean Sea. Moreover, we provide a comprehensive sensitivity analysis of the model's skill by perturbing the parameters of the bio-optical properties and phytoplankton physiology. Analysis of the inversion indicates that the model is able to reconstruct the variability of the optical constituents. Results indicate that chlorophyll-a and coloured dissolved organic matter play a major role in light modulation. The sensitivity analysis shows that the parameterization of the ratio of chlorophyll-a to carbon is important for the performance of the inversion model. A coherent inversion model, as presented, can be used as an observational operator to assimilate remote sensing reflectance.</p><p>Satellites sensors provide useful data at different temporal and spatial scales to reconstruct the variability of marine ecosystems. Ocean color sensors in Earth orbit allow the derivation of spectral water leaving radiance and inference of physical and biogeochemical properties of the water masses 1 . In simple terms, the color of the sea is related to dissolved and particulate matter present in the water that for most of the ocean are of biological origin. In recent years, suitable algorithms have been developed to estimate the biogeochemical state of the oceans from measured radiance 2 . In particular, the derivation of water biogeochemical properties can be formulated in terms of the inversion of a forward problem. Mathematically, the so-called forward description resolves light propagation according to the properties of the medium in which the light propagates; the corresponding mathematical framework is well established 3 . In parallel, the inversion algorithms use the available information about the light field to retrieve the Inherent Optical Properties (IOPs) of the medium in which the light propagates. The major optical constituents of seawater in open oceans are phytoplankton, chromophoric dissolved organic matter (CDOM), and non algal particles (NAP) such as organic and mineral particles, bacteria, viruses and air bubbles. These biogeochemical components are important indicators of the ecosystem trophic regime and carbon pool formation 4 and are among the most common products derived from inversion algorithms 2 . In addition, these indicators are extremely useful to validate surface dynamics of coupled hydrodynamic biogeochemical ocean models, by comparing the distribution of satellite-derived biogeochemical products with corresponding distributions derived from model simulations 5 .</p><p>Traditionally, biogeochemical models used a simple Beer-Lambert formulation to describe photosynthetically available radiation (PAR) propagation along the water column, i.e. the downwelling irradiance integrated over visible wavelengths (from 400 to 700 nm) that decays exponentially due to attenuation. In recent years, a</p></div>

Assessment of water quality indicators, such as Suspended Particulate Matter (SPM) concentration and water turbidity (TUR) is essential for marine ecosystems health evaluations. This work focuses on the estimation of these two variables from water-leaving reflectance measurements, either derived from in-situ observations or satellite data, for the western Black Sea basin. The regional characteristics of the water optically active constituents can have important impact on the quality of such estimations if not properly accounted for. We test several existing inversion algorithms and quantify their accuracy. New, regionally adapted methods are then proposed, together with a new formulation for a multi-conditional switching mechanism. Calibration of these models is performed based on SPM and TUR in-situ measurements. It is shown that the improvements achieved through regional adaptation can be significant. Application of these algorithms to Sentinel-3 OLCI satellite data reveals the consistency of the proposed methodology. Also, it shows the degree of uncertainty if improper formulations are used, with major impact on the absolute values of retrieved SPM or TUR.

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.

Chromophoric dissolved organic matter (CDOM) significantly contributes to the non-water absorption budget in the Mediterranean Sea. The absorption coefficient of CDOM, aCDOM(λ), is measurable in situ and can be retrieved remotely, although ocean-colour algorithms do not distinguish it from the absorption of detritus. These observations can be used as indicators for the concentration of other relevant biogeochemical variables in the ocean, e.g. dissolved organic carbon. However, our ability to model the biogeochemical processes that determine CDOM concentrations is still limited. Here we propose a novel parameterization of the CDOM cycle that accounts for the interplay between the light- and nutrient-dependent dynamics of local CDOM production and degradation, as well as its vertical transport. The parameterization is included in a one-dimensional (1D) configuration of the Biogeochemical Flux Model (BFM), which is here coupled to the General Ocean Turbulence Model (GOTM) through the Framework for Aquatic Biogeochemical Models (FABM). Here the BFM is augmented with a bio-optical component that resolves spectrally the underwater light transmission. We run this new GOTM-(FABM)-BFM configuration to simulate the seasonal aCDOM(λ) cycle at the deep-water site of the Bouée pour l'acquisition de Séries Optiques à Long Terme (BOUSSOLE) project in the northwestern Mediterranean Sea. Our results show that accounting for both nutrient and light dependence of CDOM production improves the simulation of the seasonal and vertical dynamics of aCDOM(λ), including a subsurface maximum that forms in spring and progressively intensifies in summer. Furthermore, the model consistently reproduces the higher-than-average concentrations of CDOM per unit chlorophyll concentration observed at BOUSSOLE. The configuration, outputs, and sensitivity analyses from this 1D model application will be instrumental for future applications of BFM to the entire Mediterranean Sea in a three-dimensional configuration.

Synthèse des ateliers du colloque de synthèse de prospective du domaine Océan-Atmosphère (OA) des 10-13 janvier 2023 à Autrans et des groupes de travail préparatoire. Sous la coordination de Jean-Francois Doussin – Directeur adjoint scientifique Océan-Atmosphère de l’INSU et Thibault de Garidel-Thoron – Président de la CSOA.

The Société Civile Mercator Océan International (MOi) carries out a large range of activities developed over its 25-years history, increasing its ocean expertise, and enlarging its offer in terms of services for multiple users and stakeholders of the society. A consultation of the members of the MOi management board based on seven fields of activities has been carried out to collect shareholders’ requirements for the potential framework, objectives, and goals of the future Intergovernmental Organisation (IGO) MERCATOR to be developed from the transition of the current organization. The present document aims to provide the Mercator Ocean shareholders with a comprehensive view of their respective expectations and with a common background to interact with their national delegations preparing the IGO. A comprehensive report of the collected requirements is proposed for these seven areas, as well as annexes with the detailed wording provided by each shareholder. A summary short is given below, split in the 4 categories of activities (operations / R&D / international / capacity development) adopted by the IGO Board of Delegates to describe the organisation’s tasks

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.

Using Argo profiling floats, cruises and mooring data, we reconstructed the dissolved oxygen (O2) dynamics in the Gulf of Lion and the Ligurian Sea, with a focus on the intermediate waters. By applying the CANYON-MED neural network-based method on the large network of O2-equipped Argo floats we derived nutrients and carbonate system variables in the Gulf of Lion and the Ligurian Sea at different depths in the water column and derived trends over the 2012-2020 period. In these waters, the O2 minimum is strongly affected by the intermittent convection process, and the two areas show dissimilar responses to the mixing events. In the absence of deep convection events, the O2-depleted layer tends to spread vertically and intensify even more so in the Ligurian than in the Gulf of Lion. In both areas, over the 2012-2020 period, nutrients increase overall in deep layers, with a concomitant impact on nutrient molar ratios tending toward an increase in P-limitation. Acidification estimates derived in different layers of the water column show an overall increase in dissolved inorganic carbon and a concurrent pH decrease. These trends were strongly affected by convection events slowing down the overall acidification trend.

We simulate and analyze the effects of a high CO 2 emission scenario on the Mediterranean Sea biogeochemical state at the end of the XXI century, with a focus on carbon cycling, budgets and fluxes, within and between the Mediterranean subbasins, and on ocean acidification. As a result of the overall warming of surface water and exchanges at the boundaries, the model results project an increment in both the plankton primary production and the system total respiration. However, productivity increases less than respiration, so these changes yield to a decreament in the concentrations of total living carbon, chlorophyll, particulate organic carbon and oxygen in the epipelagic layer, and to an increment in the DIC pool all over the basin. In terms of mass budgets, the large increment in the dissolution of atmospheric CO 2 results in an increment of most carbon fluxes, including the horizontal exchanges between eastern and western sub-basins, in a reduction of the organic carbon component, and in an increament of the inorganic one. The eastern sub-basin accumulates more than 85% of the absorbed atmospheric CO 2. A clear ocean acidification signal is observed all over the basin, quantitatively similar to those projected in most oceans, and well detectable also down to the mesopelagic and bathypelagic layers.

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.

The Mediterranean Sea has been identified as a hotspot for climate change. Furthermore, its very diverse trophic regimes, in such a little area, make it an extremely interesting region from a biogeochemical perspective. Numerous studies aim at better understanding and representing the Mediterrenean dynamics and biogeochemistry through modeling. This is a crucial step in order to predict the future anthropogenic impacts on the Mediterranean Sea and their possible effects on its biogeochemistry,  and all what depends on it. The number of models that simulate the Mediterranean biogeochemistry, and the data available to compare with are now sufficient to draw an overall picture of the Mediterranean Sea biogeochemical models state of the art.</p><p>In this study, we gathered 10 biogeochemical simulations of the Mediterranean Sea, including 8 regional and 2 high-resolution global configurations. The simulations are compared with surface chlorophyll estimates derived from satellite observations; chlorophyll, nitrate, oxygen, and particulate organic carbon concentrations derived from BGC-Argo floats, and  phytoplankton group-specific  primary production estimated from ocean color satellite observations. </p><p>Our first aim is to describe and compare all known Mediterranean biogeochemical models, and to highlight their specificity. This should give an insight into the current achievements, and expose what biogeochemical model products are hence available for further ecological analysis. </p><p>Furthermore, a specific attention is given to how well each model performs in selected regions of the Mediterranean Sea, in order to understand which specific process is needed to adequately represent the different trophic regimes of the Mediterranean Sea.</p><p> </p>

The Mediterranean Monitoring and Forecasting Center (Med-MFC) is part of the Copernicus Marine Environment Monitoring Service (CMEMS) and operationally produces analysis, forecast and reanalysis products for the Mediterranean Sea hydrodynamics, waves and biogeochemistry. The modelling systems are based on state-of-the-art community models, assimilate observational in situ and satellite observations and are forced by high resolution atmospheric fields. Improvements and functioning of the Med-MFC systems are based on the full consistency among the three components which are jointly upgraded and include a continuous amelioration of the accuracy of the products. The focus of this work is to present the Med-MFC modelling systems and the available products, their skill assessment, main recent achievements and future upgrades.

MOOSE is a multidisciplinary integrated Ocean observing system part of the French national Research Infrastructure for coastal ocean and seashore observations (ILICO-RI). It was established in 2010 to monitor the Northwestern Mediterranean Sea in the context of rapid climate change and its impacts on marine ecosystems.

This study investigates the spatial and temporal variability of chromophoric-dissolved organic matter (CDOM) in the Mediterranean Sea. The analysis is carried out using a state-of-the-art 3D biogeochemical model. The model describes the plankton dynamics, the cycles of the most important limiting nutrients, and the particulate and dissolved pools of carbon. The source of CDOM is directly correlated to the dynamics of dissolved organic carbon (DOC) by a fixed production quota. Then CDOM degrades by photobleaching and remineralization. The main innovation of the system is the inclusion of a bio-optical radiative transfer model that computes surface upwelling irradiance, and therefore simulates remotely sensed reflectance (Rrs). Simulation results of three model configurations are evaluated using satellite Rrs, particularly at 412 nm, 443 nm, and 490 nm. All simulations show a winter minimum in Rrs for the considered bands. However, different parameterizations of DOC-release induce a different accumulation of CDOM, especially in the eastern Mediterranean, and a different Rrs signature: a more active microbial loop during summer implies a decrease of Rrs at 412 nm. We demonstrate how the usage of a bio-optical model allows us to corroborate hypotheses on CDOM-cycling based on blue–violet Rrs data, supporting the importance of this complementary data stream with respect to satellite-derived chlorophyll.

A radiative transfer model was parameterized and validated using Biogeochemical Argo float data acquired between 2012 and 2017 across the Mediterranean Sea. Fluorescence-derived chlorophyll urn:x-wiley:21699275:media:jgrc24745:jgrc24745-math-0001 concentration, particulate optical backscattering at 700 nm, and fluorescence of chromophoric dissolved organic matter (CDOM) were used to parametrize the light absorption and scattering coefficients of the optically significant water constituents (such as pure water, non-algal particles, CDOM, and phytoplankton). The model was validated with in situ downwelling irradiance profiles and apparent optical properties derived both from irradiance profiles and satellite data, such as the diffuse attenuation coefficients and remote sensing reflectance. Results showed that by using regional parameterizations that are not only related to chlorophyll concentration and vertical distribution, the model was able to capture a more accurate spectral response in the examined wavelength range compared to chlorophyll-related (or Case 1) optical models. When using alternative models that incorporated also measurements of CDOM fluorescence or particulate optical backscattering, the model skill increased at all examined wavelengths. Finally, using a multi-spectral optical configuration also enabled the estimation of the relative contribution of separate water constituents in the examined spectral range. Simulations including non-algal particles and CDOM performed up to 61% and 79% better than when considering the optical properties of pure seawater alone. Moreover, a simulation including phytoplankton light absorption resulted in an error reduction of up to 42%, especially at 412 nm and with a more uniform response at the wavelengths considered.

Stratified oceanic systems are characterized by the presence of a so-called Deep Chlorophyll a Maximum (DCM) not detectable by ocean color satellites. A DCM can either be a phytoplankton (carbon) biomass maximum (Deep Biomass Maximum, DBM), or the consequence of photoacclimation processes (Deep photoAcclimation Maximum, DAM) resulting in the increase of chlorophyll a per phytoplankton carbon. Even though these DCM (further qualified as either DBMs or DAMs) have long been studied, no global-scale assessment has yet been undertaken and large knowledge gaps still remain in relation to the environmental drivers responsible for their formation and maintenance. In order to investigate their spatial and temporal variability in the open ocean, we use a global data set acquired by more than 500 Biogeochemical-Argo floats given that DCMs can be detected from the comparative vertical distribution of chlorophyll a concentrations and particulate backscattering coefficients. Our findings show that the seasonal dynamics of the DCMs are clearly region-dependent. High-latitude environments are characterized by a low occurrence of intense DBMs, restricted to summer. Meanwhile, oligotrophic regions host permanent DAMs, occasionally replaced by DBMs in summer, while subequatorial waters are characterized by permanent DBMs benefiting from favorable conditions in terms of both light and nutrients. Overall, the appearance and depth of DCMs are primarily driven by light attenuation in the upper layer. Our present assessment of DCM occurrence and of environmental conditions prevailing in their development lay the basis for a better understanding and quantification of their role in carbon budgets (primary production and 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.

Measuring the underwater light field is a key mission of the international Biogeochemical-Argo program. Since 2012, 0–250 dbar profiles of downwelling irradiance at 380, 412 and 490 nm besides photosynthetically available radiation (PAR) have been acquired across the globe every 1 to 10 days. The resulting unprecedented amount of radiometric data has been previously quality-controlled for real-time distribution and ocean optics applications, yet some issues affecting the accuracy of measurements at depth have been identified such as changes in sensor dark responsiveness to ambient temperature, with time and according to the material used to build the instrument components. Here, we propose a quality-control procedure to solve these sensor issues to make Argo radiometry data available for delayed-mode distribution, with associated error estimation. The presented protocol requires the acquisition of ancillary radiometric measurements at the 1000 dbar parking depth and night-time profiles. A test on >10,000 profiles from across the world revealed a quality-control success rate >90% for each band. The procedure shows similar performance in re-qualifying low radiometry values across diverse oceanic regions. We finally recommend, for future deployments, acquiring daily 1000 dbar measurements and one night profile per year, preferably during moonless nights and when the temperature range between the surface and 1000 dbar is the largest.

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.

A novel pan-European marine model ensemble was established, covering nearly all seas under the regulation of the Marine Strategy Framework Directive (MSFD), with the aim of providing a consistent assessment of the potential impacts of riverine nutrient reduction scenarios on marine eutrophication indicators. For each sea region, up to five coupled biogeochemical models from institutes all over Europe were brought together for the first time. All model systems followed a harmonised scenario approach and ran two simulations, which varied only in the riverine nutrient inputs. The load reductions were evaluated with the catchment model GREEN and represented the impacts due to improved management of agriculture and wastewater treatment in all European river systems. The model ensemble, comprising 15 members, was used to assess changes to the core eutrophication indicators as defined within MSFD Descriptor 5. In nearly all marine regions, riverine load reductions led to reduced nutrient concentrations in the marine environment. However, regionally the nutrient input reductions led to an increase in the non-limiting nutrient in the water, especially in the case of phosphate concentrations in the Black Sea. Further core eutrophication indicators, such as chlorophyll-a, bottom oxygen and the Trophic Index TRIX, improved nearly everywhere, but the changes were less pronounced than for the inorganic nutrients. The model ensemble displayed strong consistency and robustness, as most if not all models indicated improvements in the same areas. There were substantial differences between the individual seas in the speed of response to the reduced nutrient loads. In the North Sea ensemble, a stable plateau was reached after only three years, while the simulation period of eight years was too short to obtain steady model results in the Baltic Sea. The ensemble exercise confirmed the importance of improved management of agriculture and wastewater treatments in the river catchments to reduce marine eutrophication. Several shortcomings were identified, the outcome of different approaches to compute the mean change was estimated and potential improvements are discussed to enhance policy support. Applying a model ensemble enabled us to obtain highly robust and consistent model results, substantially decreasing uncertainties in the scenario outcome.

Biogeochemical Argo floats (BGC-Argo) provide an unprecedented availability of high-resolution biogeochemical and optical vertical profiles at near real time. The integration of biogeochemical and optical observations with marine ecosystem models allows to improve the model capability to describe marine ecosystem dynamics at different spatial and temporal scales, also resulting in an increase of the model skill. Recent advancements and future upgrades of the Copernicus Marine Service (CMS) Mediterranean Sea biogeochemical modelling system include the use of BGC-Argo data for validation and assimilation of both biogeochemical and optical components. The focus of this work is to present the upgrade of the BGC-Argo data stream quality check for the CMS Mediterranean Sea biogeochemical modelling system workflow and to discuss a novel skill assessment framework oriented to evaluate key biogeochemical processes and ecosystem dynamics (e.g. deep chlorophyll maximum depth, nitracline depth, minimum oxygen depth) and optical characteristics (bbp700 converted data) that benefit from the particularly rich and high quality level of the BGC-Argo network in this semi-enclosed sea.

A multiplatform assessment of the Ocean–Atmosphere Spectral Irradiance Model (OASIM) radiative model focussed on the Mediterranean Sea for the period 2004–2017 is presented. The BOUée pour l'acquiSition d'une Série Optique à Long termE (BOUSSOLE) mooring and biogeochemical Argo (BGC-Argo) float optical sensor observations are combined with model outputs to analyse the spatial and temporal variabilities in the downward planar irradiance at the ocean–atmosphere interface. The correlations between the data and model are always higher than 0.6. With the exception of downward photosynthetic active radiation and the 670 nm channel, correlation values are always higher than 0.8 and, when removing the inter-daily variability, they are higher than 0.9. At the scale of the BOUSSOLE sampling (15 min temporal resolution), the root mean square difference oscillates at approximately 30 %–40 % of the averaged model output and is reduced to approximately 10 % when the variability between days is filtered out. Both BOUSSOLE and BGC-Argo indicate that bias is up to 20 % for the irradiance at 380 and 412 nm and for wavelengths above 670 nm, whereas it decreases to less than 5 % at the other wavelengths. Analysis of atmospheric input data indicates that the model skill is strongly affected by cloud dynamics. High skills are observed during summer when the cloud cover is low.

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.

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.

A regional neural network-based method, "CANYON-MED" is developed to estimate nutrients and carbonate system variables specifically in the Mediterranean Sea over the water column from pressure, temperature, salinity, and oxygen together with geolocation and date of sampling. Six neural network ensembles were developed, one for each variable (i.e., three macronutrients: nitrates (NO − 3), phosphates (PO 3− 4) and silicates (SiOH 4), and three carbonate system variables: pH on the total scale (pH T), total alkalinity (A T), and dissolved inorganic carbon or total carbon (C T), trained using a specific quality-controlled dataset of reference "bottle" data in the Mediterranean Sea. This dataset is representative of the peculiar conditions of this semi-enclosed sea, as opposed to the global ocean. For each variable, the neural networks were trained on 80% of the data chosen randomly and validated using the remaining 20%. CANYON-MED retrieved the variables with good accuracies (Root Mean Squared Error): 0.73 µmol.kg −1 for NO − 3 , 0.045 µmol.kg −1 for PO 3− 4 and 0.70 µmol.kg −1 for Si(OH) 4 , 0.016 units for pH T , 11 µmol.kg −1 for A T and 10 µmol.kg −1 for C T. A second validation on the ANTARES independent time series confirmed the method's applicability in the Mediterranean Sea. After comparison to other existing methods to estimate nutrients and carbonate system variables, CANYON-MED stood out as the most robust, using the aforementioned inputs. The application of CANYON-MED on the Mediterranean Sea data from autonomous observing systems (integrated network of Biogeochemical-Argo floats, Eulerian moorings and ocean gliders measuring hydrological properties together with oxygen concentration) could have a wide range of applications. These include data quality control or filling gaps in time series, as well as biogeochemical data assimilation and/or the initialization and validation of regional biogeochemical models still lacking crucial reference data. Matlab and R code are available at https:// github.com/MarineFou/CANYON-MED/.

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.

Linear regression is widely used in applied sciences and, in particular, in satellite optical oceanography, to relate dependent to independent variables. It is often adopted to establish empirical algorithms based on a finite set of measurements, which are later applied to observations on a larger scale from platforms such as autonomous profiling floats equipped with optical instruments (e.g., Biogeochemical Argo floats; BGC-Argo floats) and satellite ocean colour sensors (e.g., SeaWiFS, VIIRS, OLCI). However, different methods can be applied to a given pair of variables to determine the coefficients of the linear equation fitting the data, which are therefore not unique. In this work, we quantify the impact of the choice of "regression method" (i.e., either type-I or type-II) to derive bio-optical relationships, both from theoretical perspectives and by using specific examples. We have applied usual regression methods to an in situ data set of particulate organic carbon (POC), total chlorophyll-a (TChla), optical particulate backscattering coefficient (b bp), and 19 years of monthly TChla and b bp ocean colour data. Results of the regression analysis have been used to calculate phytoplankton carbon biomass (C phyto) and POC from: i) BGC-Argo float observations; ii) oceanographic cruises, and iii) satellite data. These applications enable highlighting the differences in C phyto and POC estimates relative to the choice of the method. An analysis of the statistical properties of the dataset and a detailed description of the hypothesis of the work drive the selection of the linear regression method.

The emerging availability of Biogeochemical-Argo (BGC-Argo) float data creates new opportunities to combine models and observations in investigations of the interior structures and dynamics of marine ecosystems. An existing variational data assimilation scheme (3DVarBio) has been upgraded and coupled with the Copernicus Marine Environment Monitoring Service biogeochemical model of the Mediterranean Sea to assimilate BGC-Argo chlorophyll profile observations. Our results show that the assimilation of BGC-Argo float data is feasible. Moreover, the proposed data assimilation framework provides significant corrections to the chlorophyll concentrations and is able to consistently readjust the shapes of chlorophyll profiles during surface blooms occurring in winter vertically mixed conditions, and in the case of the summer deep chlorophyll maxima. Sensitivity analysis and diagnostic metrics have been applied to evaluate the impact of assimilation and the relevance of different factors of the 3DVarBio method. A key factor is the observation error that has been tuned on a monthly basis to incorporate the representation error. Different frequencies of the assimilation cycle have been tested: daily or 3-day assimilation fosters the highest skill performances despite the reduced impacts and the increase of computational burden. Considering the present size of the BGC-Argo Mediterranean network (about 15 floats) and the estimated non-homogeneous correlation radius length scale (15.2 km on average), the chlorophyll assimilation can constrain the phytoplankton dynamics along the whole water column over an area up to 10% of the Mediterranean Sea.

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.

The Mediterranean community represented in this paper is the result of more than 30 years of EU and nationally funded coordination, which has led to key contributions in science concepts and operational initiatives. Together with the establishment of operational services, the community has coordinated with universities, research centers, research infrastructures and private companies to implement advanced multi-platform and integrated observing and forecasting systems that facilitate the advancement of operational services, scientific achievements and mission-oriented innovation. Thus, the community can respond to societal challenges and stakeholders needs, developing a variety of fit-for-purpose services such as the Copernicus Marine Service. The combination of state-of-the-art observations and forecasting provides new opportunities for downstream services in response to the needs of the heavily populated Mediterranean coastal areas and to climate change. The challenge over the next decade is to sustain ocean observations within the research community, to monitor the variability at small scales, e.g., the mesoscale/submesoscale, to resolve the sub-basin/seasonal and inter-annual variability in the circulation, and thus establish the decadal variability, understand and correct the model-associated biases and to enhance model-data integration and ensemble forecasting for uncertainty estimation. Better knowledge and understanding of the level of Mediterranean variability will enable a subsequent evaluation of the impacts and mitigation of the effect of human activities and climate change on the biodiversity and the ecosystem, which will support environmental assessments and decisions. Further challenges include extending the science-based added-value products into societal relevant downstream services and engaging with communities to build initiatives that will contribute to the 2030 Agenda and more specifically to SDG14 and the UN's Decade of Ocean Science for sustainable development, by this contributing to bridge the science-policy gap. The Mediterranean observing and forecasting capacity was built on the basis of community best practices in monitoring and modeling, and can serve as a basis for the development of an integrated global ocean observing system.

New autonomous robotic platforms for observing the ocean, i.e. Biogeochemical-Argo (BGC-Argo) floats, have drastically increased the number of vertical profiles of irradiance, photosynthetically available radiation (PAR), and algal chlorophyll concentrations around the globe independent of the season. Such data may therefore be a fruitful resource to improve performances of numerical models for marine biogeochemistry. Here we present a work that integrates 1314 vertical profiles of PAR acquired by 31 BGC-Argo floats operated in the Mediterranean Sea between 2012 and 2016 into a one-dimensional model to simulate the vertical and temporal variability of algal chlorophyll concentrations. The model was initially forced with PAR measurements to assess its skill when using quality-controlled light profiles, and subsequently with a number of alternative bio-optical models to analyse the model capability when light observations are not available. Model outputs were evaluated against co-located chlorophyll profiles measured by BGC-Argo floats. Results highlight that the data-driven model is able to reproduce the spatial and temporal variability of deep chlorophyll maxima depth observed at a number of Mediter-ranean sites well. Further, we illustrate the key role of PAR and vertical mixing in shaping the vertical dynamics of primary producers in the Mediterranean Sea. The comparison of alternative bio-optical models identifies the best simple one to be used, and suggests that model simulations benefit from considering the diel cycle.

The Argo Program has been implemented and sustained for almost two decades, as a global array of about 4000 profiling floats. Argo provides continuous observations of ocean temperature and salinity versus pressure, from the sea surface to 2000 dbar. The successful installation of the Argo array and its innovative data management system arose opportunistically from the combination of great scientific need and technological innovation. Through the data system, Argo provides fundamental physical observations with broad societally-valuable applications, built on the cost-efficient and robust technologies of autonomous profiling floats. Following recent advances in platform and sensor technologies, even greater opportunity exists now than 20 years ago to (i) improve Argo’s global coverage and value beyond the original design, (ii) extend Argo to span the full ocean depth, (iii) add biogeochemical sensors for improved understanding of oceanic cycles of carbon, nutrients, and ecosystems, and (iv) consider experimental sensors that might be included in the future, for example to document the spatial and temporal patterns of ocean mixing. For Core Argo and each of these enhancements, the past, present, and future progression along a path from experimental deployments to regional pilot arrays to global implementation is described. The objective is to create a fully global, top-to-bottom, dynamically complete, and multidisciplinary Argo Program that will integrate seamlessly with satellite and with other in situ elements of the Global Ocean Observing System (Legler et al., 2015). The integrated system will deliver operational reanalysis and forecasting capability, and assessment of the state and variability of the climate system with respect to physical, biogeochemical, and ecosystems parameters. It will enable basic research of unprecedented breadth and magnitude, and a wealth of ocean-education and outreach opportunities.

Understanding spatial and temporal dynamics of non-algal particles (NAP) in open ocean is of the utmost importance to improve estimations of carbon export and sequestration. These particles covary with phytoplankton abundance but also accumulate independently of algal dynamics. The latter likely represents an important fraction of organic carbon but it is largely overlooked. A possible way to study these particles is via their optical backscattering properties (b bp) and relationship with chlorophyll-a (Chl). To this aim, we estimate the fraction of b bp associated with the NAP portion () that does not covary with Chl by using a global Biogeochemical-Argo dataset. We quantify the spatial, temporal and vertical variability of. In the northern productive areas, is a small fraction of b bp and shows a clear seasonal cycle. In the Southern Ocean, b k bp is a major fraction of total b bp. In oligotrophic areas, has a smooth annual cycle.

Biogeochemical (BGC)-Argo floats observations are becoming a major data source for assimilation into and constraining of ocean biogeochemical models. An important prerequisite for a successful synthesis between models and observations is the characterization of the observational errors in BGC-Argo float data. The root-mean-square error and multiplicative and additive biases in qualitycontrolled data sets of oxygen, nitrate, and chlorophyll a concentrations collected with 17 BGC-Argo floats in the Mediterranean Sea between 2013 and 2017 are assessed using the triple collocation analysis. The analysis suggests that BGC-Argo float oxygen, nitrate and chlorophyll a data suffer from an additive bias of 2.9 ± 5.5 μmol/kg, 0.46 ± 0.07 μmol/kg, and −0.06 ± 0.02 mg/m 3 , respectively. The root-mean-square error is evaluated at 5.1 ± 0.8 μmol/kg, 0.25 ± 0.07 μmol/kg, and 0.03 ± 0.01 mg/m 3. Additional studies should determine whether these values are applicable to the global ocean. Plain Language Summary The Biogeochemical-Argo program is a network of ocean robots whose sensors monitor oxygen, nitrate, and chlorophyll a concentration information that is needed to detect decadal changes in biological carbon production, ocean acidification, ocean carbon uptake, and hypoxia in the world ocean. One of the goals of the Biogeochemical-Argo program is to incorporate these observations into ocean models to understand and forecast the changing state of the carbon cycle. The successful integration of the float data into numerical models, however, requires the specification of the observational errors. This study provides, for the first time, the biases and errors of the three cores variables of the Biogeochemical-Argo floats network: oxygen, nitrate, and chlorophyll a concentrations.

A wavelet analysis has been applied, for the first time, to 3-year high-frequency field observations of bio-optical properties (i.e. chlorophyll-fluorescence, beam attenuation and backscattering coefficients) in the northwestern Mediterranean Sea (BOUSSOLE site), in order to identify their dominant temporal patterns and evolution. A cross-wavelet and coherence analysis has also been applied to paired bio-optical coefficients time-series at the BOUSSOLE site, which allows identifying the temporal relationship between the cycles of the bio-optical properties. Annual, six-and four-month, intra-seasonal (i.e., mid-and short-terms) cycles are identified from the time-series analysis. The periodicities of chlorophyll-fluorescence, beam attenuation and particulate back-scattering coefficients correlate well at different temporal scales and specific seasons. At annual, six-and four-month scales, different bio-optical properties follow rather similar patterns, likely driven by physical forcing. Intra-seasonal variability consists in both mid-and short-term variations. The former dominates during the winter and are related to episodic bloom events, while the latter variations (i.e., diel) prevail during summer, in a stratified water column.

Characterizing phytoplankton distribution and dynamics in the world's open oceans requires in situ observations over a broad range of space and time scales. In addition to temperature/salinity measurements , Biogeochemical-Argo (BGC-Argo) profiling floats are capable of autonomously observing at high-frequency bio-optical properties such as the chlorophyll fluorescence, a proxy of the chlorophyll a concentration (Chla), the particulate backscattering coefficient (b bp), a proxy of the stock of particulate organic carbon , and the light available for photosynthesis. We analyzed an unprecedented BGC-Argo database of more than 8,500 multivariable profiles collected in various oceanic conditions, from subpolar waters to subtropical gyres. Our objective is to refine previously established Chla versus b bp relationships and gain insights into the sources of vertical, seasonal, and regional variability in this relationship. Despite some regional, seasonal and vertical variations, a general covariation occurs at a global scale. We distinguish two main contrasted situations: (1) concomitant changes in Chla and b bp that correspond to actual variations in phytoplankton biomass, e.g., in subpolar regimes; (2) a decoupling between the two variables attributed to photoacclima-tion or changes in the relative abundance of nonalgal particles, e.g., in subtropical regimes. The variability in the b bp :Chla ratio in the surface layer appears to be essentially influenced by the type of particles and by photoacclimation processes. The large BGC-Argo database helps identifying the spatial and temporal scales at which this ratio is predominantly driven by one or the other of these two factors.

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.

The northwestern Mediterranean Sea is one of the few sites of open-sea deep convection and dense water formation. The area is characterized by intense air-sea exchanges favored by the succession of strong northerly and northwesterly winds during autumn and winter that eventually induce deep convection episodes and the formation of the Western Mediterranean Deep Water. This region exhibits a significant spring phytoplankton bloom, which appears to be largely influenced in intensity and diversity by the winter mixing. To understand and resolve the interplay between the atmosphere and ocean, and the impact of ocean circulation and mixing on biogeochemistry, we carried out an unprecedented observational ffort during a large experiment between July 2012 and July 2013. A multiplatform approach—combining aircraft, balloons, ships, moorings,floats, and gliders - aimed both at characterizing the dynamics of the atmosphere and at quantifying the physical, biogeochemical, and biological properties of the water masses.Beyond a better understanding of the wind dynamics, the interannual variability of the deep convection, and the seasonal variability of the nutrient distribution and plankton structure in the pelagic ecosystem, the experiment provided a reference data set that was used as a benchmark for advancing the modeling of the surface fluxes, convective processes, dense water formation rates, and physical-biogeochemical coupling processes. It also represented an opportunity for complementary investigations such as evaluating model parameterizations or studying the role of submesoscale eddies in dense water spreading and biogeochemistry.

A new special issue of JGR: Oceans and JGR: Atmospheres presents new insights into the dynamics of dense water formation in the western Mediterranean Sea and its biogeochemical consequences.

Currently, observations from low-Earth orbit (LEO) ocean color sensors represent one of the most used tools to study surface optical and biogeochemical properties of the ocean. LEO observations are available at daily temporal resolution, and are often combined into weekly, monthly, seasonal, and annual averages in order to obtain sufficient spatial coverage. Indeed, daily satellite maps of the main oceanic variables (e.g., surface phytoplankton chlorophyll-a) generally have many data gaps, mainly due to clouds, which can be filled using either Optimal Interpolation or the Empirical Orthogonal Functions approach. Such interpolations, however, may introduce large uncertainties in the final product. Here, our goal is to quantify the potential benefits of having high-temporal resolution observations from a geostationary (GEO) ocean color sensor to reduce interpolation errors in the reconstructed hourly and daily chlorophyll-a products. To this aim, we used modeled chlorophyll-a fields from the Copernicus Marine Environment Monitoring Service's (CMEMS) Baltic Monitoring and Forecasting Centre (BAL MFC) and satellite cloud observations from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor (on board the geostationary satellite METEOSAT). The sampling of a GEO was thus simulated by combining the hourly chlorophyll fields and clouds masks, then hourly and daily chlorophyll-a products were generated after interpolation from neighboring valid data using the Multi-Channel Singular Spectral Analysis (M-SSA). Two cases are discussed: (i) A reconstruction based on the typical sampling of a LEO and, (ii) a simulation of a GEO sampling with hourly observations. The results show that the root mean square and interpolation bias errors are significantly reduced using hourly observations.

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.

An efficient system to produce in situ high quality radiometric measurements is compulsory to rigorously perform the vicarious calibration of satellite sensors dedicated to Ocean Color Radiometry (OCR) and to validate their derived products. This requirement is especially needed during the early stages of an OCR satellite activity or for remote areas poorly covered by oceanographic cruises with possible bio-optical anomalies. Taking advantage of Argo's profiling float technology, we present a new autonomous profiling float dedicated to in situ radiometric measurements. The float is based on the Provor CTS5 (manufacturer NKE) with an added novel two protruding arm design allowing for sensor redundancies, shading mitigation and near-surface data. Equipped with two identical radiometers on each arm that measure downward irradiance and upwelling radiance at seven wavelengths, the ProVal float generates both redundant radiometric profiles as well as an estimate of Remote Sensing Reflectance. Results from 449 profiles obtained in the NW Mediterranean Sea and in the Indian sector of the Southern Ocean are presented to illustrate the ProVal float technical maturity. Analysis of the behavior of the profiling float, including tilting and ascent speeds is presented. The vertical stability of the ProVal exhibits 85% of surface data of the Mediterranean Sea with a tilt smaller than 10 degrees. This percentage is 40% in the Southern Ocean due to rougher seas. Redundant sensors provide a characterization of the relative drift between sensors over the deployment which is found to be <0.15% per month over a year. Post-cruise calibration of a recovered float revealed no significant drift. As an example of the utility of ProVal floats, a match-up of Remote Sensing Reflectance measured with the European Space Agency Ocean and Land Color Imager (OLCI onboard Sentinel-3A) is shown. It follows that profiling floats, such as ProVal, could provide a significant contribution to an upcoming global System Vicarious Calibration of space-based radiometers.

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.

Regionalisation aims at delimiting provinces within which physical conditions, chemical properties, and biological communities are reasonably homogeneous. This article proposes a synthesis of the many recent regionalisations of the open-sea regions of the Mediterranean Sea. The nine studies considered here defined regions based on different, and sometimes complementary, criteria: dynamics of surface chlorophyll concentration, ocean currents, three-dimensional hydrological and biogeochemical properties, or the distribution of organisms. Although they identified different numbers and patterns of homogeneous regions, their compilation in the epipelagic zone identifies nine consensus frontiers, eleven consensus regions with relatively homogeneous conditions, and four heterogeneous regions with highly dynamical conditions. The consensus frontiers and regions are in agreement with well-known hydrodynamical features of the Mediterranean Sea, which constrain the distribution of hydrological and ecological variables. The heterogeneous regions are rather defined by intense mesoscale activity. The synthesis proposed here could constitute a reference step for management actions and spatial planning, such as the application of the European Marine Strategy Framework Directive, and for future biogeochemical and ecological studies in the Mediterranean Sea.

Chlorophyll fluorometers provide the largest in situ global data set for estimating phytoplankton biomass because of their ease of use, size, power consumption, and relatively low price. While in situ chlorophyll a (Chl) fluorescence is proxy for Chl a concentration, and hence phytoplankton biomass, there exist large natural variations in the relationship between in situ fluorescence and extracted Chl a concentration. Despite this large natural variability, we present here a global validation data set for the WET Labs Environmental Characterization Optics (ECO) series chlorophyll fluorometers that suggests a factor of 2 overestimation in the factory calibrated Chl a estimates for this specific manufacturer and series of sensors. We base these results on paired High Pressure Liquid Chromatography (HPLC) and in situ fluorescence match ups for which non-photochemically quenched fluorescence observations were removed. Additionally, we examined match-ups between the factory-calibrated in situ fluorescence and estimates of chlorophyll concentration determined from in situ radiometry, absorption line height, NASA's standard ocean color algorithm as well as laboratory calibrations with phytoplankton monocultures spanning diverse species that support the factor of 2 bias. We therefore recommend the factor of 2 global bias correction be applied for the WET Labs ECO sensors , at the user level, to improve the global accuracy of chlorophyll concentration estimates and products derived from them. We recommend that other fluorometer makes and models should likewise undergo global analyses to identify potential bias in factory calibration.

Satellite ocean color observations revealed that unusually deep convection events in 2005, 2006, 2010, and 2013 led to an increased phytoplankton biomass during the spring bloom over a large area of the northwestern Mediterranean Sea (NWM). Here we investigate the effects of these events on the seasonal phytoplankton community structure, we quantify their influence on primary production, and we discuss the potential biogeochemical impact. For this purpose, we compiled in situ phytoplankton pigment data from five ship surveys performed in the NWM and from monthly cruises at a fixed station in the Ligurian Sea. We derived primary production rates from a light photosynthesis model applied to these in situ data. Our results confirm that the maximum phytoplankton biomass during the spring bloom is larger in years associated with intense deep convection events (151%). During these enhanced spring blooms, the contribution of diatoms to total phytoplankton biomass increased (133%), as well as the primary production rate (1115%). The occurrence of a highly productive bloom is also related to an increase in the phytoplankton bloom area (1155%) and in the relative contribution of diatoms to primary production (163%). Therefore, assuming that deep convection in the NWM could be significantly weakened by future climate changes, substantial decreases in the spring production of organic carbon and of its export to deep waters can be expected.

In June 2013, a glider equipped with oxygen and fluorescence sensors has been used to extensively sample an anticyclonic Submesoscale Coherent Vortex (SCV) in the Ligurian Sea (NW Mediterranean Sea). Those measurements are complemented by full‐depth CTD casts (T, S, and oxygen) and water samples documenting nutrients and phytoplankton pigments within the SCV and outside. The SCV has a very homogeneous core of oxygenated waters between 300 and 1200 m formed 4.5 months earlier during the winter deep convection event. It has a strong dynamical signature with peak velocities at 700 m depth of 13.9 cm s⁻¹ in cyclogeostrophic balance. The eddy has a small radius of 6.2 km corresponding to high Rossby number of −0.45. The vorticity at the eddy center reaches 0.8f. Cross‐stream isopycnic diffusion of tracers between the eddy core and the surroundings is found to be very limited due to dynamical barriers set by the SCV associated with a diffusivity coefficient of about 0.2 m² s⁻¹. The deep core is nutrients‐depleted with concentrations of nitrate, phosphate, and silicate, 13–18% lower than the rich surrounding waters. However, the nutriclines are shifted of about 20–50 m toward the surface thus increasing the nutrients availability for phytoplankton. Chlorophyll‐a concentrations at the deep chlorophyll maximum are subsequently about twice bigger as compared to outside. Pigments further reveal the predominance of nanophytoplankton inside the eddy and an enhancement of the primary productivity. This study demonstrates the important impact of postconvective SCVs on nutrients distribution and phytoplankton community, as well as on the subsequent primary production and carbon sequestration.

The Southern Ocean (SO), an area highly sensitive to climate change, is currently experiencing rapid warming and freshening. Such drastic physical changes might significantly alter the SO's biological pump. For more accurate predictions of the possible evolution of this pump, a better understanding of the environmental factors controlling SO phytoplankton dynamics is needed. Here we present a satellite-based study deciphering the complex environmental control of phytoplankton biomass (PB) and phenology (PH; timing and magnitude of phytoplankton blooms) in the SO. We reveal that PH and PB are mostly organized in the SO at two scales: a large latitudinal scale and a regional scale. Latitudinally, a clear gradient in the timing of bloom occurrence appears tightly linked to the seasonal cycle in irradiance, with some exceptions in specific light-limited regimes (i.e., well-mixed areas). Superimposed on this latitudinal scale, zonal asymmetries, up to 3 orders of magnitude, in regional-scale PB are mainly driven by local advective and iron supply processes. These findings provide a global understanding of PB and PH in the SO, which is of fundamental interest for identifying and explaining ongoing changes as well as predicting future changes in the SO biological pump.

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.

Since 2012, an array of 105 Biogeochemical-Argo (BGC-Argo) floats has been deployed across the world's oceans to assist in filling observational gaps that are required for characterizing open-ocean environments. Profiles of biogeochemical (chlorophyll and dissolved organic matter) and optical (single-wavelength particulate optical backscattering, downward irradiance at three wavelengths, and photosynthetically available radiation) variables are collected in the upper 1000m every 1 to 10 days. The database of 9837 vertical profiles collected up to January 2016 is presented and its spatial and temporal coverage is discussed. Each variable is quality controlled with specifically developed procedures and its time series is quality-assessed to identify issues related to biofouling and/or instrument drift. A second database of 5748 profile-derived products within the first optical depth (i.e., the layer of interest for satellite remote sensing) is also presented and its spatiotemporal distribution discussed. This database, devoted to field and remote ocean color applications, includes diffuse attenuation coefficients for downward irradiance at three narrow wavebands and one broad waveband (photosynthetically available radiation), calibrated chlorophyll and fluorescent dissolved organic matter concentrations, and single-wavelength particulate optical backscattering. To demonstrate the applicability of these databases, data within the first optical depth are compared with previously established bio-optical models and used to validate remotely derived bio-optical products. The quality-controlled databases are publicly available from the SEANOE (SEA scieNtific Open data Edition) publisher at https://doi.org/10.17882/49388 and https://doi.org/10.17882/47142 for vertical profiles and products within the first optical depth, respectively.

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.

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.

In the 80s, the POEM (Physical Oceanography of the Eastern Mediterranean) cruises in the Levantine Basin first revealed the presence of a very pronounced dynamical structure off Cyprus: The Cyprus Warm Core Eddy. Since then, a large amount of data have been collected thanks to the use of autonomous oceanic gliders (+8000 profiles since 2009). Part of those profiles were carried out in the upper layers down to 200 m, and we take benefit of a novel approach named ITCOMP SOM that uses a statistical approach to extend them down to 1000 m (see [1] for more details). This dataset have a particularly good spatio-temporal coverage in 2009 for about a month, thanks to simultaneous deployments of several gliders (up to 6). In this study, we present a set of 3D reconstruction of the dynamical and hydrographical characteristics of the Warm Core Cyprus Eddy between 2009 and 2015. Moreover, chlorophyll-a fluorescence data measured by the gliders give evidence to strong vertical velocities at the edge of the eddy. We discuss possible mechanisms (frontogenesis, symmetric instability) that could generate such signals and provide an assessment of the role of this peculiar circulation feature on the circulation and biogeochemistry of the Levantine basin. Reference: [1] Charantonis, A., P. Testor, L. Mortier, F. D'Ortenzio, S. Thiria (2015): Completion of a sparse GLIDER database using multi-iterative Self-Organizing Maps (ITCOMP SOM), Procedia Computer Science, 51(1):2198-2206. DOI: 10.1016/j.procs.2015.05.496

D'Ortenzio and Ribera d'Alcalà (2009, DR09 hereafter) divided the Mediterranean Sea into " bioregions " based on the climatological seasonality (phenology) of phy-toplankton. Here we investigate the interannual variability of this bioregionalization. Using 16 years of available ocean color observations (i.e., SeaWiFS and MODIS), we analyzed the spatial distribution of the DR09 trophic regimes on an annual basis. Additionally, we identified new trophic regimes, exhibiting seasonal cycles of phytoplankton biomass different from the DR09 climatological description and named " Anomalous ". Overall, the classification of the Mediter-ranean phytoplankton phenology proposed by DR09 (i.e., " No Bloom " , " Intermittently " , " Bloom " and " Coastal "), is confirmed to be representative of most of the Mediterranean phytoplankton phenologies. The mean spatial distribution of these trophic regimes (i.e., bioregions) over the 16 years studied is also similar to the one proposed by DR09, although some annual variations were observed at regional scale. Discrepancies with the DR09 study were related to interannual variability in the sub-basin forcing: winter deep convection events, frontal instabilities, inflow of Atlantic or Black Sea Waters and river runoff. The large assortment of phytoplank-ton phenologies identified in the Mediterranean Sea is thus verified at the interannual scale, further supporting the " sen-tinel " role of this basin for detecting the impact of climate changes on the pelagic environment.

The evolution of the stratification of the north-western Mediterranean between summer 2012 and the end of winter 2013 was simulated and compared with different sets of observations. A summer cruise and profiler observations were used to improve the initial conditions of the simulation. This improvement was crucial to simulate winter convection. Variations of some parameters involved in air - sea exchanges (wind, coefficient of transfer used in the latent heat flux formulation, and constant additive heat flux) showed that the characteristics of water masses and the volume of dense water formed during convection cannot be simply related to the time-integrated buoyancy budget over the autumn - winter period. The volume of dense water formed in winter was estimated to be about 50,000 km³ with a density anomaly larger than 29.113 kg m^-3. The effect of advection and air/sea fluxes on the heat and salt budget of the convection zone was quantified during the preconditioning phase and the mixing period. Destratification of the surface layer in autumn occurs through an interaction of surface and Ekman buoyancy fluxes associated with displacements of the North Balearic front bounding the convection zone to the south. During winter convection, advection stratifies the convection zone: from December to March, the absolute value of advection represents 58 % of the effect of surface buoyancy fluxes.

We present here a unique oceanographic and meteorological data set focus on the deep convection processes. Our results are essentially based on in situ data (mooring, research vessel, glider, and profiling float) collected from a multiplatform and integrated monitoring system (MOOSE: Mediterranean Ocean Observing System on Environment), which monitored continuously the northwestern Mediterranean Sea since 2007, and in particular high-frequency potential temperature, salinity, and current measurements from the mooring LION located within the convection region. From 2009 to 2013, the mixed layer depth reaches the seabed, at a depth of 2330m, in February. Then, the violent vertical mixing of the whole water column lasts between 9 and 12 days setting up the characteristics of the newly formed deep water. Each deep convection winter formed a new warmer and saltier “vintage” of deep water. These sudden inputs of salt and heat in the deep ocean are responsible for trends in salinity (3.3 ± 0.2 × 10−3/yr) and potential temperature (3.2 ± 0.5 × 10−3 C/yr) observed from 2009 to 2013 for the 600–2300 m layer. For the first time, the overlapping of the three “phases” of deep convection can be observed, with secondary vertical mixing events (2–4 days) after the beginning of the restratification phase, and the restratification/spreading phase still active at the beginning of the following deep convection event.

An array of Bio-Argo floats equipped with radiometric sensors has been recently deployed in various open ocean areas representative of the diversity of trophic and bio-optical conditions prevailing in the so-called case 1 waters. Around solar noon and almost every day, each float acquires 0-250-m vertical profiles of photosynthetically available radiation and downward irradiance at three wavelengths (380, 412, and 490 nm). Up until now, more than 6500 profiles for each radiometric channel have been acquired. As these radiometric data are collected out of an operator's control and regardless of meteorological conditions, specific and automatic data processing protocols have to be developed. This paper presents a data quality-control procedure aimed at verifying profile shapes and providing near-real-time data distribution. This procedure is specifically developed to 1) identify main issues of measurements (i.e., dark signal, atmospheric clouds, spikes, and wavefocusing occurrences) and 2) validate the final data with a hierarchy of tests to ensure a scientific utilization. The procedure, adapted to each of the four radiometric channels, is designed to flag each profile in a way compliant with the data management procedure used by the Argo program. Main perturbations in the light field are identified by the new protocols with good performances over the whole dataset. This highlights its potential applicability at the global scale. Finally, the comparison with modeled surface irradiances allows for assessing the accuracy of quality-controlled measured irradiance values and identifying any possible evolution over the float lifetime due to biofouling and instrumental drift.

The HYdrological cycle in the Mediterranean Experiment (HyMeX) Special Observing Period 2 (SOP2, January 27–March 15, 2013) was dedicated to the study of dense water formation in the Gulf of Lion in the northwestern Mediterranean. This paper outlines the deep convection of winter 2012–2013 and the meteorological conditions that produced it. Alternating phases of mixing and restratification are related to periods of high and low heat losses, respectively. High-resolution, realistic, three-dimensional models are essential for assessing the intricacy of buoyancy fluxes, horizontal advection, and convective processes. At the submesoscale, vertical velocities resulting from symmetric instabilities of the density front bounding the convection zone are crucial for ventilating the deep ocean. Finally, concomitant atmospheric and oceanic data extracted from the comprehensive SOP2 data set highlight the rapid, coupled evolution of oceanic and atmospheric boundary layer characteristics during a strong wind event.

The present study proposes a novel method that merges satellite ocean color bio-optical products with Argo temperature-salinity profiles to infer the vertical distribution of the particulate backscattering coefficient (bbp). This neural network-based method (SOCA-BBP for Satellite Ocean-Color merged with Argo data to infer the vertical distribution of the Particulate Backscattering coefficient) uses three main input components: (1) satellite-based surface estimates of bbp and chlorophyll a concentration matched up in space and time with (2) depth-resolved physical properties derived from temperature-salinity profiles measured by Argo profiling floats and (3) the day of the year of the considered satellite-Argo matchup. The neural network is trained and validated using a database including 4725 simultaneous profiles of temperature-salinity and bio-optical properties collected by Bio-Argo floats, with concomitant satellite-derived products. The Bio-Argo profiles are representative of the global open-ocean in terms of oceanographic conditions, making the proposed method applicable to most open-ocean environments. SOCA-BBP is validated using 20% of the entire database (global error of 21%). We present additional validation results based on two other independent data sets acquired (1) by four Bio-Argo floats deployed in major oceanic basins, not represented in the database used to train the method; and (2) during an AMT (Atlantic Meridional Transect) field cruise in 2009. These validation tests based on two fully independent data sets indicate the robustness of the predicted vertical distribution of bbp. To illustrate the potential of the method, we merged monthly climatological Argo profiles with ocean color products to produce a depth-resolved climatology of bbp for the global ocean.

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.

Agulhas rings provide the principal route for ocean waters to circulate from the Indo-Pacific to the Atlantic basin. Their influence on global ocean circulation is well known, but their role in plankton transport is largely unexplored. We show that, although the coarse taxonomic structure of plankton communities is continuous across the Agulhas choke point, South Atlantic plankton diversity is altered compared with Indian Ocean source populations. Modeling and in situ sampling of a young Agulhas ring indicate that strong vertical mixing drives complex nitrogen cycling, shaping community metabolism and biogeochemical signatures as the ring and associated plankton transit westward. The peculiar local environment inside Agulhas rings may provide a selective mechanism contributing to the limited dispersal of Indian Ocean plankton populations into the Atlantic.

The distribution of the chlorophyll a concentration ([Chl a]) in the Mediterranean Sea, mainly obtained from satellite surface observations or from scattered in situ experiments, is updated by analyzing a database of fluorescence profiles converted into [Chl a]. The database, which includes 6790 fluorescence profiles from various origins, was processed with a specific quality control procedure. To ensure homogeneity between the different data sources, 65 % of fluorescence profiles have been intercalibrated on the basis of their concomitant satellite [Chl a] estimation. The climatological pattern of [Chl a] vertical profiles in four key sites of the Mediterranean Sea has been analyzed. Climatological results confirm previous findings over the range of existing [Chl a] values and throughout the principal Mediterranean trophic regimes. They also provide new insights into the seasonal variability in the shape of the vertical [Chl a] profile, inaccessible through remote-sensing observations. An analysis based on the recognition of the general shape of the fluorescence profile was also performed. Although the shape of [Chl a] vertical distribution characterized by a deep chlorophyll maximum (DCM) is ubiquitous during summer, different forms are observed during winter, thus suggesting that factors affecting the vertical distribution of the biomass are complex and highly variable. The [Chl a] spatial distribution in the Mediterranean Sea mimics, on smaller scales, what is observed in the global ocean. As already evidenced by analyzing satellite surface observations, midlatitude- and subtropical-like phytoplankton dynamics coexist in the Mediterranean Sea. Moreover, the Mediterranean DCM variability appears to be characterized by patterns already observed on the global scale.

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.

Atmospheric fluxes of bio-available inorganic nitrogen (DIN, i.e. nitrate+ammonium) and phosphorus (DIP, i.e. phosphate) were measured in 2010, 2011 and 2013 at the sampling station of Cap Ferrat (Ligurian Sea). Wet and dry fluxes of DIN, averaged over three years, were 35 and 19 mmol m−2 yr−1, respectively. Most of the nitrate was deposited under dry form, whilst ammonia was twice more found in wet deposition. Wet and dry fluxes of DIP, averaged over three years, were 0.11 and 0.64 mmol m−2 yr-1, respectively. Atmospheric fluxes of DIN and DIP were compared with other photic zone nutrient input sources, physical and biological, i.e. winter convection, N2 fixation, and upward diffusion. Even if convection is by far the most important nutrient input for surface waters, atmospheric sources may be the second one, supplying more nutrients than diazotrophy and diffusion, particularly in conditions of water column stratification.

Abstract. In vivo chlorophyll a fluorescence is a proxy of chlorophyll a concentration, and is one of the most frequently measured biogeochemical properties in the ocean. Thousands of profiles are available from historical databases and the integration of fluorescence sensors to autonomous platforms has led to a significant increase of chlorophyll fluorescence profile acquisition. To our knowledge, this important source of environmental data has not yet been included in global analyses. A total of 268 127 chlorophyll fluorescence profiles from several databases as well as published and unpublished individual sources were compiled. Following a robust quality control procedure detailed in the present paper, about 49 000 chlorophyll fluorescence profiles were converted into phytoplankton biomass (i.e., chlorophyll a concentration) and size-based community composition (i.e., microphytoplankton, nanophytoplankton and picophytoplankton), using a method specifically developed to harmonize fluorescence profiles from diverse sources. The data span over 5 decades from 1958 to 2015, including observations from all major oceanic basins and all seasons, and depths ranging from the surface to a median maximum sampling depth of around 700 m. Global maps of chlorophyll a concentration and phytoplankton community composition are presented here for the first time. Monthly climatologies were computed for three of Longhurst's ecological provinces in order to exemplify the potential use of the data product. Original data sets (raw fluorescence profiles) as well as calibrated profiles of phytoplankton biomass and community composition are available on open access at PANGAEA, Data Publisher for Earth and Environmental Science.

An analysis of seasonal variations in climatological surface chlorophyll points to distinct biogeographical zones in the North Atlantic subpolar gyre. In particular, the Labrador Sea appears well delineated into two regions on either side of the 60°N parallel, with very different climatological phytoplankton biomass cycles. Indeed, north of 60°N, an early and short spring bloom occurs in late April, while south of 60°N, the bloom gradually develops 1 month later and significant biomass persists all summer long. Nevertheless, at climatological scale, the first-order mechanism that controls the bloom is identical for both bioregions. The light-mixing regime can explain the bloom onset in both bioregions. In the Labrador Sea, the blooms seem to rely on a mean community compensation irradiance threshold value of 2.5 mol photon m−2 d−1 over the mixed layer.

Thanks to the new genera.on of profiling floats and in par.cular to iridium telemetry, the acquisi.on frequency can be drama.cally increased. We inves.gate here the possibility to extract informa.on related to sea-state from the analysis of high-resolu.on measurements of the pressure data. We par.cularly focus on the study of the speed anomaly as compared to a nominal speed expected for a calm sea-state. By comparison between speed anomaly of a float in the Med Sea and concurrent seastate measurements by a weather buoy in the same area, we suggest that float behaviour can be an indicator of sea-state. In the context of remOcean and NOAS projects, we set up a high frequency mode (every 2 s) for the sub-surface layer and for more than fourty floats deployed in various open ocean areas, we present a preliminary analysis of the speed anomaly. observa.on of float behavior : Characteris.cs of the NKE CTS4 float in the upper layer:-when the float passes 10 dbars, it does not ac.vate its pump for a dura.on of 600 sec-ader 600 sec, the float starts to pump 360 sec to emerge-under standard (calm) weather condi.ons, it takes 90-100 sec for the float to rise the surface from 10m.-data acquisi.on from 1m to 0 m: mean value-data acquisi.on from 10 to 1m : raw data @ 2 sec (0.5Hz)-data acquisi.on from 350 to 10m : raw data @ 10 sec (0.1Hz) CTS4, NKE 600 s 4.5 m 80 s 0 m Standard behavior sta.s.c & climagraphs Non Standard behavior MOORED AUTOMATIC WEATHER STATIONS valida.on with weather mooring Time-serie of behavior indicators of two floats compared to meteorological buoy data in the NorthWestern Mediterranean Sea. Even if the distance of the floats to the mooring are greater than 100 km, the comparison of both time series confirms that float behavior indicators can be used with a certain confidence to track sea-state. These indicators can thus be extracted from all the CTS4 floats we manage to date. We now envisage to extract similar float behavior indicators from standard Argo profiles and in particular the information related to big event.. big events > 10m re-dive 0.4% what next Characteris.cs of the NKE CTS4 float in Indicators of float behavior in the upper layer In this poster we explore the idea that the float behavior is related to the sea state. Time spend between 10m and 1m (s) proxy of speed anomaly Maximum re−dive between 10m and 1m [m],

We present a Mediterranean climatology (1° × 1° × 12 months) of the mixed layer and of the seasonal thermocline, based on a comprehensive collection of temperature profiles spanning 44 years (1969–2012). The database includes more than 190,000 profiles, merging CTD, MBT/XBT, profiling floats, and gliders observations. This data set is first used to describe the seasonal cycle of the mixed layer depth and temperature, together with the seasonal thermocline depth and averaged temperature, on the whole Mediterranean on a monthly climatological basis. Our analysis discriminates several regions with coherent behaviors, in particular the deep water formation sites, characterized by significant differences in the winter mixing intensity. Heat Storage Rate (HSR) is calculated as the time rate of change of the heat content due to variations in the temperature integrated from the surface down to the base of the seasonal thermocline. For the first time the quantification of heat storage rate in the upper-ocean, based only on in situ oceanographic data, is made for the whole Mediterranean. The spatial and temporal variability of the HSR in the Mediterranean Sea and its link with dynamic structures like oceanic gyres are also discussed.

A neural network-based method is developed to assess the vertical distribution of (1) chlorophyll a concentration ([Chl]) and (2) phytoplankton community size indices (i.e., microphytoplankton, nanophytoplankton, and picophytoplankton) from in situ vertical profiles of chlorophyll fluorescence. This method (FLAVOR for Fluorescence to Algal communities Vertical distribution in the Oceanic Realm) uses as input only the shape of the fluorescence profile associated with its acquisition date and geo-location. The neural network is trained and validated using a large database including 896 concomitant in situ vertical profiles of High-Performance Liquid Chromatography (HPLC) pigments and fluorescence. These profiles were collected during 22 oceanographic cruises representative of the global ocean in terms of trophic and oceanographic conditions, making our method applicable to most oceanic waters. FLAVOR is validated with respect to the retrieval of both [Chl] and phytoplankton size indices using an independent in situ data set and appears to be relatively robust spatially and temporally. To illustrate the potential of the method, we applied it to in situ measurements of the BATS (Bermuda Atlantic Time Series Study) site and produce monthly climatologies of [Chl] and associated phytoplankton size indices. The resulting climatologies appear very promising compared to climatologies based on available in situ HPLC data. With the increasing availability of spatially and temporally well-resolved data sets of chlorophyll fluorescence, one possible global-scale application of FLAVOR could be to develop 3-D and even 4-D climatologies of [Chl] and associated composition of phytoplankton communities. The Matlab and R codes of the proposed algorithm are provided as supporting information.

Nitrate, phosphate, and silicate concentration profiles were measured at monthly frequency at the DYFAMED time-series station (central Ligurian Sea) between 1991 and 2011. The resulting data set, which constitutes the longest open-ocean time-series in the Mediterranean Sea, underwent quality control. A reproducible climatological pattern was observed with an unprecedented resolution, confirming the typical seasonal cycle of mid-latitudes. In summer and autumn, when the water mass is well stratified, i.e. the mixed layer depth (MLD) is shallow, nutrient concentrations in surface are very low or under the detection limit. In winter, as a result of the MLD extent, nutrients are supplied to the surface layer. Then, nutrient concentrations progressively decrease during spring. MLD appears to play a key role in controlling nutrient availability in the surface layer, but a direct, quantitative relationship between MLD and nutrient concentrations is difficult to establish due to undersampling. Regarding nutrient molar ratios (N:P, Si:N, and Si:P), results show anomalous values compared to those of other oceanic regions, presumably due to strong influence of external sources. As a consequence, nutrient molar ratios exhibit a seasonal pattern, with, in particular, an increase of the N:P ratio in condition of stratification. Over the period 1991–2011, the DYFAMED data set reveals decadal trends in nitrate and phosphate concentrations in deep waters (+0.23% and –0.62%, respectively) resulting in increasing N:P and Si:P ratios (+1.14% and +0.85% per year, respectively). Such a long-term variability is presumably related to changes in water mass and/or changes in external sources, even if it is difficult to assess due to not enough concomitant data from atmospheric and riverine inputs.

Two profiling floats, equipped with nitrate concentration sensors were deployed in the northwestern Mediterranean from summer 2012 to summer 2013. Satellite ocean color data were extracted to evaluate surface chlorophyll concentration at float locations. Time series of mixed layer depths and nitrate and chlorophyll concentrations were analyzed to characterize the interplay between the physical-chemical and biological dynamics in the area. Deep convection (mixed layer depth > 1000 m) was observed in January-February, although high-nitrate surface concentrations could be already observed in December. Chlorophyll increase is observed since December, although high values were observed only in March. The early nitrate availability in subsurface layers, which is likely due to the permanent cyclonic circulation of the area, appears to drive the bloom onset. The additional nitrate supply associated to the deep convection events, although strengthening the overall nitrate uptake, seems decoupled of the December increase of chlorophyll.

We deployed four Bio-Argo profiling floats in various oligotrophic locations of the Pacific subtropical gyres and Mediterranean Sea to address the seasonal phytoplankton dynamics in the euphotic layer and explore its dependence on light regime dynamics. Results show that there is a similar phytoplankton biomass seasonal pattern in the four observed oceanic regions. In the lower part of the euphotic layer, the seasonal displacement of the deep chlorophyll maximum (DCM) is light driven. During winter, the chlorophyll a concentration ([Chl a]) always increases in the upper euphotic mixed layer. This increase always results from a photoacclimation to the reduced irradiance. Depending on the location, however, the concentration can also be associated with an actual increase in biomass. The winter increase in [Chl a] results in an increase in irradiance attenuation that impacts the position of the isolume (level where the daily integrated photon flux is constant) and DCM, which becomes shallower. In summer when the [Chl a] in the upper layer decreases along with light attenuation, the DCM deepens and becomes closer to (and sometimes reaches) the nitracline, which enhances the phytoplankton biomass at the DCM. The bio-optical mechanisms and their relationship to light regimes that are revealed by the time series appear to be generic and potentially characteristic of all of the areas where a DCM forms, which is 50% of the open ocean.

nazionale per le nuove tecnologie, l'energia e lo sviluppo economico sostenibile,

The temporal evolution of the vertical export flux at the DYFAMED time-series station (Ligurian Sea) over the last 20 years reveals a strong interannual variability. Winter convection allows particulate (and dissolved) matter to be vertically exported (“flush-down” effect). The efficiency of this process determines also the concentration of nutrients brought to surface waters and, therefore, the intensity of the subsequent phytoplankton bloom. The sequence “convection-bloom” is the main driving force of vertical export flux in this region. The present work attempts to better identify the parameters that control vertical export flux dynamics by observing a 20 year time-series in relation with the temporal variability of mixed layer depth and surface primary production. The consequences of a more stratified water column in the future on biological productivity and vertical export flux are pointed out.

The winter of 2012 experienced peculiar atmospheric conditions that triggered a massive formation of dense water on the continental shelf and in the deep basin of the Gulf of Lions. Multiplatforms observations enabled a synoptic view of dense water formation and spreading at basin scale. Five months after its formation, the dense water of coastal origin created a distinct bottom layer up to a few hundreds of meters thick over the central part of the NW Mediterranean basin, which was overlaid by a layer of newly formed deep water produced by open-sea convection. These new observations highlight the role of intense episodes of both dense shelf water cascading and open-sea convection to the progressive modification of the NW Mediterranean deep waters.

Phytoplankton phenology is primarily affected by physical forcing. However, its quantification is far from being completely understood. Among the physical forcing factors, the mixed layer depth (MLD) is considered to have the strongest impact on phytoplankton dynamics, and consequently, on their phenology. The role of MLD variations in shaping the phytoplankton phenology was explored in the Mediterranean Sea, a basin displaying contrasting phenological regimes. A database of MLD estimations was merged with ocean color chlorophyll concentrations ([Chl]SAT) to generate concomitant annual MLD and [Chl]SAT cycles. Several indices were calculated to quantitatively analyze these cycles. The relevance of indices summarizing the temporal difference between main characteristics of MLD and [Chl]SAT cycles was emphasized. As previously observed, two dominant phenological regimes coexist in the Mediterranean Sea. The first is marked by a typical spring bloom, as in temperate regions. The second displays a low seasonality and an absence of an intense [Chl]SAT peak as in subtropical areas. The MLD is shown to play a key role in determining the dominant phenological regime in a given area. Results also show that regions having low seasonality display concomitant MLD and [Chl]SAT maxima, whereas [Chl]SAT peaks are generally observed 30 days after MLD peaks in regions with strongest seasonality. Over the whole basin, [Chl]SAT increase starts 1 month after the initiation of MLD deepening. Finally, after examining the impact of MLD on light and nutrient availability for phytoplankton, mechanisms were proposed to explain the time lags between MLD and [Chl]SAT increase and MLD and [Chl]SAT maxima.

The winter of 2012 experienced peculiar atmospheric conditions that triggered a massive formation of dense water on the continental shelf and in the deep basin of the Gulf of Lions. Multi-platforms observations enabled, with an unprecedented resolution, a synoptic view of dense water formation and spreading at basin scale. Five months after its formation, the dense water of coastal origin created a distinct bottom layer up to few hundreds of meters thick over the central part of the NW Mediterranean basin, which was overlaid by a layer of newly formed deep water produced by open-sea convection. These observations highlight the role of intense episodes of both dense shelf water cascading and open-sea convection to the alteration of the characteristics of the NW Mediterranean deep waters.

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.

Terrestrial inputs (natural and anthropogenic) from rivers, the atmosphere and physical processes strongly impact the functioning of coastal pelagic ecosystems. The objective of this study was to develop a tool for the examination of these impacts on the Marseille coastal area, which experiences inputs from the Rhone River and high rates of atmospheric deposition. Therefore, a new 3D coupled physical/biogeochemical model was developed. Two versions of the biogeochemical model were tested, one model considering only the carbon (C) and nitrogen (N) cycles and a second model that also considers the phosphorus (P) cycle. Realistic simulations were performed for a period of 5 years (2007– 2011). The model accuracy assessment showed that both versions of the model were able of capturing the seasonal changes and spatial characteristics of the ecosystem. The model also reproduced upwelling events and the intrusion of Rhone River water into the Bay of Marseille well. Those processes appeared to greatly impact this coastal oligotrophic area because they induced strong increases in chlorophyll-a concentrations in the surface layer. The model with the C, N and P cycles better reproduced the chlorophyll-a concentrations at the surface than did the model without the P cycle, especially for the Rhone River water. Nevertheless, the chlorophyll-a concentrations at depth were better represented by the model without the P cycle. Therefore, the complexity of the biogeochemical model introduced errors into the model results, but it also improved model results during specific events. Finally, this study suggested that in coastal oligotrophic areas, improvements in the description and quantification of the hydrodynamics and the terrestrial inputs should be preferred over increasing the complexity of the biogeochemical model. Citation: Fraysse M, Pinazo C, Faure VM, Fuchs R, Lazzari P, et al. (2013) Development of a 3D Coupled Physical-Biogeochemical Model for the Marseille Coastal Area (NW Mediterranean Sea): What Complexity Is Required in the Coastal Zone? PLoS ONE 8(12): e80012.

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.

In situ observation of the marine environment has traditionally relied on ship-based platforms. The obvious consequence is that physical and biogeochemical properties have been dramatically undersampled, especially in the remote Southern Ocean (SO). The difficulty in obtaining in situ data represents the major limitations to our understanding, and interpretation of the coupling between physical forcing and the biogeochemical response. Southern elephant seals (Mirounga leonina) equipped with a new generation of oceanographic sensors can measure ocean structure in regions and seasons rarely observed with traditional oceanographic platforms. Over the last few years, seals have allowed for a considerable increase in temperature and salinity profiles from the SO, but we were still lacking information on the spatiotemporal variation of phytoplankton concentration. This information is critical to assess how the biological productivity of the SO, with direct consequences on the amount of CO2 "fixed" by the biological pump, will respond to global warming. In this research programme, we use an innovative sampling fluorescence approach to quantify phytoplankton concentration at sea. For the first time, a low energy consumption fluorometer was added to Argos CTD-SRDL tags, and these novel instruments were deployed on 27 southern elephant seals between 25 December 2007 and the 4 February 2011. As many as 3388 fluorescence profiles associated with temperature and salinity measurements were thereby collected from a vast sector of the Southern Indian Ocean. This paper addresses the calibration issue of the fluorometer before being deployed on elephant seals and presents the first results obtained for the Indian sector of the Southern Ocean. This in situ system is implemented in synergy with satellite ocean colour radiometry. Satellite-derived data is limited to the surface layer and is restricted over the SO by extensive cloud cover. However, with the addition of these new tags, we are able to assess the 3-dimension distribution of phytoplankton concentration by foraging southern elephant seals. This approach reveals that for the Indian sector of the SO, the surface chlorophyll a (chl a) concentrations provided by MODIS were underestimated by a factor 2 compared to chl a concentrations estimated from HPLC corrected in situ fluorescence measurements. The scientific outcomes of this programme include an improved understanding of both the present state and variability in ocean biology, and the accompanying biogeochemistry, as well as the delivery of real-time and open-access data to scientists (doi:10.7491/MEMO.1x)

Heat fluxes across the ocean-atmosphere interface play a crucial role in the upper turbulent mixing. The depth reached by this turbulent mixing is indicated by an homogenization of seawater properties in the surface layer, and is defined as the Mixed Layer Depth (MLD). The thickness of the mixed layer determines also the heat content of the layer that directly interacts with the atmosphere. The seasonal variability of these air-sea fluxes is crucial in the calculation of heat budget. An improvement in the estimate of these fluxes is needed for a better understanding of the Mediterranean ocean circulation and climate, in particular in Regional Climate Models. There are few estimations of surface heat fluxes based on oceanic observations in the Mediterranean, and none of them are based on mixed layer observations. So, we proposed here new estimations of these upper-ocean heat fluxes based on mixed layer. We present high resolution Mediterranean climatology (0.5°) of the mean MLD based on a comprehensive collection of temperature profiles of last 43 years (1969-2012). The database includes more than 150,000 profiles, merging CTD, XBT, ARGO Profiling floats, and gliders observations. This dataset is first used to describe the seasonal cycle of the mixed layer depth on the whole Mediterranean on a monthly climatological basis. Our analysis discriminates several regions with coherent behaviors, in particular the deep water formation sites, characterized by significant differences in the winter mixing intensity. Heat storage rates (HSR) were calculated as the time rate of change of the heat content integrated from the surface down to a specific depth that is defined as the MLD plus an integration constant. Monthly climatology of net heat flux (NHF) from ERA-Interim reanalysis was balanced by the 1°x1° resolution heat storage rate climatology. Local heat budget balance and seasonal variability in the horizontal heat flux are then discussed by taking into account uncertainties, due to errors in monthly value estimation and to intra-annual and inter-annual variability.

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.

The deep ocean is the largest and least known ecosystem on Earth. It hosts numerous pelagic organisms, most of which are able to emit light. Here we present a unique data set consisting of a 2.5-year long record of light emission by deep-sea pelagic organisms, measured from December 2007 to June 2010 at the ANTARES underwater neutrino telescope in the deep NW Mediterranean Sea, jointly with synchronous hydrological records. This is the longest continuous time-series of deep-sea bioluminescence ever recorded. Our record reveals several weeks long, seasonal bioluminescence blooms with light intensity up to two orders of magnitude higher than background values, which correlate to changes in the properties of deep waters. Such changes are triggered by the winter cooling and evaporation experienced by the upper ocean layer in the Gulf of Lion that leads to the formation and subsequent sinking of dense water through a process known as "open-sea convection". It episodically renews the deep water of the study area and conveys fresh organic matter that fuels the deep ecosystems. Luminous bacteria most likely are the main contributors to the observed deep-sea bioluminescence blooms. Our observations demonstrate a consistent and rapid connection between deep open-sea convection and bathypelagic biological activity, as expressed by bioluminescence. In a setting where dense water formation events are likely to decline under global warming scenarios enhancing ocean stratification, in situ observatories become essential as environmental sentinels for the monitoring and understanding of deep-sea ecosystem shifts.

As the proxy for Chlorophyll a (Chl a) concentration, thousands of fluorescence profiles were measured by instrumented elephant seals in the Kerguelen region (Southern Ocean). For accurate retrieval of Chl a concentrations acquired by in vivo fluorometer, a two-step procedure is applied: 1) A predeployment intercalibration with accurate determination by high performance liquid chromatography (HPLC) analysis, which not only calibrates fluorescence in appropriate Chl a concentration units, but also strongly reduces variability between fluorometers, and 2) a profile-by-profile quenching correction analysis, which effectively eliminates the fluorescence quenching issue at surface around noon, and results in consistent profiles between day and night. The quenching correction is conducted through an extrapolation of the deep fluorescence value toward surface. As proved by a validation procedure in the Western Mediterranean Sea, the correction method is practical and relatively reliable when there is no credible reference, especially for deep mixed waters, as in the Southern Ocean. Even in the shallow mixed waters, the method is also effective in reducing the influence of quenching.

Eight autonomous profiling "Bio-Argo" floats were deployed offshore during about 2 years (2008-2010) in Pacific, Atlantic, and Mediterranean zones. They were equipped with miniaturized bio-optical sensors, namely a radiometer measuring within the upper layer the downward irradiance at 412, 490, and 555 nm, and two fluorometers for detection of chlorophyll-a (Chla) and colored dissolved organic matter (CDOM; profiles from 400 m to surface). A first study dealt with the interpretation of the Chla fluorescence signal in terms of concentration, using for this purpose the diffuse attenuation coefficient for irradiance at 490 nm, K d (490), taken as a proxy for the Chla absorption. The present study examines the possibility of similarly using the K d (412) values combined with retrieved Chla profiles to convert the CDOM fluorometric qualitative information into a CDOM absorption coefficient (a y). The rationale is to take advantage of the fact that K d is more sensitive to CDOM presence at 412 nm than at 490 nm. A validation of this method is tested through its application to field data, collected from a ship over a wide range of trophic conditions (Biogeochemistry and Optics South Pacific Experiment (BIOSOPE) cruise); these data include both in situ fluorescence profiles and CDOM absorption as measured on discrete samples. In addition, near-surface a y values retrieved from the floats agree with those derivable from ocean color imagery (Moderate Resolution Imaging Spectroradiometer (MODIS-A)). The low sensitivity of commercially available CDOM fluorometers presently raises difficulties when applying this technique to open ocean waters. It was nevertheless possible to derive from the floats records meaningful time series of CDOM vertical distribution.

Understanding the ocean carbon cycle requires a precise assessment of phytoplankton biomass in the oceans. In terms of numbers of observations, satellite data represent the largest available data set. However, as they are limited to surface waters, they have to be merged with in situ observations. Amongst the in situ data, fluorescence profiles constitute the greatest data set available, because fluorometers have operated routinely on oceanographic cruises since the 1970s. Nevertheless, fluorescence is only a proxy of the total chlorophyll a concentration and a data calibration is required. Calibration issues are, however, sources of uncertainty, and they have prevented a systematic and wide range exploitation of the fluorescence data set. In particular, very few attempts to standardize the fluorescence databases have been made. Consequently , merged estimations with other data sources (e.g. satellite) are lacking. We propose a merging method to fill this gap. It consists firstly in adjusting the fluorescence profile to impose a zero chlorophyll a concentration at depth. Secondly, each point of the fluorescence profile is then multiplied by a correction coefficient, which forces the chlorophyll a integrated content measured on the fluorescence profile to be consistent with the concomitant ocean colour observation. The method is close to the approach proposed by Boss et al. (2008) to correct fluorescence data of a profiling float, although important differences do exist. To develop and test our approach , in situ data from three open ocean stations (BATS, HOT and DYFAMED) were used. Comparison of the so-called "satellite-corrected" fluorescence profiles with con-comitant bottle-derived estimations of chlorophyll a concentration was performed to evaluate the final error (estimated at 31 %). Comparison with the Boss et al. (2008) method, using a subset of the DYFAMED data set, demonstrated that the methods have similar accuracy. The method was applied to two different data sets to demonstrate its utility. Using fluo-rescence profiles at BATS, we show that the integration of "satellite-corrected" fluorescence profiles in chlorophyll a climatologies could improve both the statistical relevance of chlorophyll a averages and the vertical structure of the chlorophyll a field. We also show that our method could be efficiently used to process, within near-real time, profiles obtained by a fluorometer deployed on autonomous platforms, in our case a bio-optical profiling float. The application of the proposed method should provide a first step towards the generation of a merged satellite/fluorescence chlorophyll a product, as the "satellite-corrected" profiles should then be consistent with satellite observations. Improved climatolo-gies with more consistent satellite and in situ data are likely to enhance the performance of present biogeochemical models .

An approach that combines a recently developed procedure for improved estimation of surface chlorophyll a concentration (Chl surf) from ocean color and a phytoplankton class-specific bio-optical model was used to examine primary production in the Mediterranean Sea. Specifically, this approach was applied to the 10 year time series of satellite Chl surf data from the Sea-viewing Wide Field-of-view Sensor. We estimated the primary production associated with three major phytoplankton classes (micro, nano, and picophytoplankton), which also yielded new estimates of the total primary production (P tot). These estimates of P tot (e.g., 68 g C m À2 yr À1 for the entire Mediterranean basin) are lower by a factor of $2 and show a different seasonal cycle when compared with results from conventional approaches based on standard ocean color chlorophyll algorithm and a non-class-specific primary production model. Nanophytoplankton are found to be dominant contributors to P tot (43-50%) throughout the year and entire basin. Micro and picophytoplankton exhibit variable contributions to P tot depending on the season and ecological regime. In the most oligotrophic regime, these contributions are relatively stable all year long with picophytoplankton ($32%) playing a larger role than microphytoplankton ($22%). In the blooming regime, picophytoplankton dominate over microphytoplankton most of the year, except during the spring bloom when microphytoplankton (27-38%) are considerably more important than picophytoplankton (20-27%).

We investigated the phenology of oceanic phytoplankton at large scales over two 5-year time periods: 1979-1983 and 1998-2002. Two ocean-color satellite data archives (Coastal Zone Color Scanner (CZCS) and Sea-viewing Wide Field-of-view Sensor (SeaWiFS)) were used to investigate changes in seasonal patterns of concentration- normalized chlorophyll. The geographic coverage was constrained by the CZCS data distribution. It was best for the Northern Hemisphere and also encompassed large areas of the Indian, South Pacific, and Equatorial Atlantic regions. For each 2° pixel, monthly climatologies were developed for satellite-derived chlorophyll, and the resulting seasonal cycles were statistically grouped using cluster analysis. Five distinct groups of mean seasonal cycles were identified for each half-decade period. Four types were common to both time periods and correspond to previously identified phytoplankton regimes: Bloom, Tropical, Subtropical North, and Subtropical South. Two other mean seasonal cycles, one in each of the two compared 5-year periods, were related to transitional or intermediate states (Transitional Tropical and Transitional Bloom). Five mean seasonal cycles (Bloom, Tropical, Subtropical North, and Subtropical South, Transitional Bloom) were further confirmed when the whole SeaWiFS data set (1998-2010) was analyzed. For ~35% of the pixels analyzed, characteristic seasonal cycles of the 1979-1983 years differed little from those of the 1998-2002 period. For ~65% of the pixels, however, phytoplankton seasonality patterns changed markedly, especially in the Northern Hemisphere. Subtropical regions of the North Pacific and Atlantic experienced a widespread expansion of the Transitional Bloom regime, which appeared further enhanced in the climatology based on the full SeaWiFS record (1998-2010), and, as showed by a more detailed analysis, is associated to La Niña years. This spatial pattern of Transitional Bloom regime reflects a general smoothing of seasonality at macroscale, coming into an apparent greater temporal synchrony of the Northern Hemisphere. The Transitional Bloom regime is also the result of a higher variability, both in space and time. The observed change in phytoplankton dynamics may be related not only to biological interactions but also to large-scale changes in the coupled atmosphere-ocean system. Some connections are indeed found with climate indices. Changes were observed among years belonging to opposite phases of ENSO, though discernible from the change among the two periods and within the SeaWiFS era (1998-2010). These linkages are considered preliminary at present and are worthy of further investigation.

Eight autonomous profiling floats equipped with miniaturized radiometers and fluorimeters have collected data in Pacific, Atlantic, and Mediterranean offshore zones. They measured in particular 0-400 m vertical profiles of the downward irradiance at three wavelengths (412, 490, and 555 nm) and of the chlorophyll a fluorescence. Such autonomous sensors collect radiometric data regardless of sky conditions and collect essentially uncalibrated fluorescence data. Usual processing and calibration techniques are no longer usable in such remote conditions and have to be adapted. The proposition here is an interwoven processing by which missing parts of irradiance profiles (due to intermittent cloud occurrence) are interpolated by accounting for possible changes in optical properties (detected by the fluorescence signal) and by which the attenuation coefficient for downward irradiance, used as proxy for [Chl a] (the chlorophyll a concentration), allows the fluorescence signal to be calibrated in absolute units (mg m −3). This method is successfully applied to about 600 irradiance and fluorescence profiles. Validation of the results in terms of [Chl a] is made by matchup with satellite (MODIS-A) chlorophyll (24.3% RMSE, N = 358). Validation of the method is obtained by applying it on similar field data acquired from ships, which, in addition to irradiance and fluorescence profiles, include the [Chl a] HPLC determination, used for final verification.

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.

Phytoplankton chlorophyll-a (Chl) seasonal cycles of the North Atlantic are described using satellite ocean color observations covering the 1980s and the 2000s. The study region is where warmer SST and higher Chl in the 2000s as compared to the 1980s have been reported. It covers latitudes from 30°N to 50°N and longitudes from 60°W to 0°W, where two phytoplankton blooms take place: a spring bloom that follows stratification of upper layers, and a fall bloom due to nutrient entrainment through deepening of the mixed layer. In the 1980s, spring and fall blooms were of similar amplitude over the entire study region. In the 2000s, the fall bloom was weaker in the eastern Atlantic (east of 40°W), because of a delayed deepening of the mixed layer at the end of summer (mixed-layer depths - MLD - determined from in situ data). Conversely, the spring bloom of the eastern Atlantic was stronger in the 2000s than it was in the 1980s, because of a deeper MLD and stronger winds in winter. In the North Western Atlantic (northwest of 38°N-40°W), little differences are observed for spring and fall blooms, and for the wintertime MLD. Our results show that the links between upper-layer stratification, SST changes, and biological responses are more complex than the simple paradigm that sequentially relates higher stratification with warmer SST and an enhanced growth of the phytoplankton population.

In vivo fluorescence of Chlorophyll-a (Chl-a) is a potentially useful property to study the vertical distribution of phytoplankton biomass. However the technique is presently not fully exploited as it should be, essentially because of the difficulties in converting the fluorescence signal into an accurate Chl-a concentration. These difficulties arise noticeably from natural variations in the Chl-a fluores-cence relationship, which is under the control of community composition as well as of their nutrient and light status. As a consequence, although vertical profiles of fluores-cence are likely the most recorded biological property in the open ocean, the corresponding large databases are underex-ploited. Here with the aim to convert a fluorescence profile into a Chl-a concentration profile, we test the hypothesis that the Chl-a concentration can be gathered from the sole knowledge of the shape of the fluorescence profile. We analyze a large dataset from 18 oceanographic cruises conducted in case-1 waters from the highly stratified hyperoligotrophic waters (surface Chl-a = 0.02 mg m −3) of the South Pacific Gyre to the eutrophic waters of the Benguela upwelling (sur-face Chl-a = 32 mg m −3) and including the very deep mixed waters in the North Atlantic (Mixed Layer Depth = 690 m). This dataset encompasses more than 700 vertical profiles of Chl-a fluorescence as well as accurate estimations of Chl-a by High Performance Liquid Chromatography (HPLC). Two typical fluorescence profiles are identified, the uniform profile , characterized by a homogeneous layer roughly corresponding to the mixed layer, and the non-uniform profile, characterized by the presence of a Deep Chlorophyll Maximum. Using appropriate mathematical parameterizations, a fluorescence profile is subsequently represented by 3 or 5 shape parameters for uniform or non-uniform profiles, respectively. For both situations, an empirical model is de-Correspondence to: A. Mignot () veloped to predict the "true" Chl-a concentration from these shape parameters. This model is then used to calibrate a flu-orescence profile in Chl-a units. The validation of the approach provides satisfactory results with a median absolute percent deviation of 33 % when comparing the HPLC Chl-a profiles to the Chl-a-calibrated fluorescence. The proposed approach thus opens the possibility to produce Chl-a clima-tologies from uncalibrated fluorescence profile databases that have been acquired in the past and to which numerous new profiles will be added, thanks to the recent availability of autonomous platforms (profiling floats, gliders and animals) in-strumented with miniature fluorometers.

The subtropical gyres represent 40% of the total ocean. They have been considered, since a long time, as steady state biological desert due to low nutrients and phytoplankton biomass throughout the year. Primary production (carbon uptake by phytoplankton) is low too in these regions, but as these structures are huge, it makes the total contribution to ocean primary production important. The South Pacific gyre is considered to be the place where the phytoplankton biomass is the lowest in the global ocean. This area remains poorly sampled and studied, that's why the seasonal variability of phytoplankton, primary production and physical forcing remains unclear. By using in-situ bio-optical profiling float capable to measure physical parameters (temperature and salinity) and proxy of the phytoplankton biomass in the water column, we acquired one year of data near the eastern islands (in the middle of the South Pacific gyre). We found a seasonal variability of the phytoplankton biomass and the carbon flux. This variability is driven by the vertical structure density of the ocean, from very low phytoplankton biomass and carbon uptake in winter (due to the mixing of the water column) to low phytoplankton biomass and carbon uptake in summer (associated with strong stratification). Studying the interaction between the physical forcing and phytoplankton biomass in the subtropical gyres is essential. Any increase in sea water temperature will change the stratification of the water column and thus affect the annual carbon fluxes budget.

The ten years of the SeaWiFS satellite surface chlorophyll concentration observations, presently available, were used to characterize the biogeography of the Mediterranean Sea and the seasonal cycle of the surface biomass in different areas of the basin. The K-means cluster analysis was applied on the satellite time-series of chlorophyll concentration. The resulting coherent patterns were then explained on the basis of the present knowledge of the basin's functioning. Winter biomass enhancements were shown to occur in most of the basin and last for 2-3 months depending on the region. Classical spring bloom regimes were also observed, regularly in the North Western Mediterranean, and intermittently in four other specific areas. The geographical correspondence between specific clusters and regions showing high values of mean chlorophyll concentration indicates that, at least in the Mediterranean Sea, accumulations of phytoplankton are observed only where specific temporal trends are present.

The match-up of satellite-derived reflectances with in situ observations is crucial to evaluate their quality and temporal stability. To contribute to this effort, a project has been set up to collect a data set of in situ radiometric and bio-optical quantities, in support to satellite ocean color calibration and validation. The project has been named ''BOUSSOLE'', and one of its key elements is a deep-sea optics mooring collecting data on a near-continuous basis since September 2003. This buoy is deployed in the deep clear waters of the northwestern Mediterranean Sea, and is visited on a monthly basis for servicing and acquisition of complementary data. The characteristics of the work area establish the site as a satisfactory location for validating satellite ocean color observations. A description of the data processing protocols is provided, followed by an analysis of the uncertainty of the buoy radiometry measurements. The results of a match-up analysis of the marine reflectances, diffuse attenuation coefficients, and chlorophyll concentrations for three major missions, i.e., MERIS, SeaWiFS, and MODIS-A, are then analyzed. They show poor performances for the bluest band (412 nm) of the three sensors, and performances within requirements at 443 and 490 nm for SeaWiFS and MODIS-A. These results suggest that a vicarious calibration should be introduced for the MERIS sensor. This analysis also demonstrates that a major effort is still required to improve atmospheric correction procedures whatever the mission.

A new 0.5°resolution Mediterranean climatology of the mixed layer depth based on individual profiles of temperature and salinity has been constructed. The criterion selected is a threshold value of temperature from a near-surface value at 10 m depth, mainly derived by a method applied on the global (de Boyer Montégut et al., 2004 dBM04). With respect to dBM04, the main differences reside in the absence of spatial interpolation of the final fields and in the improved spatial resolution. These changes to the method are necessary to reproduce the Mediterranean mixed layer's behavior. In the derived climatological maps, the most relevant features of the basin surface circulation are reproduced, as well as the areas prone of the deep water formation are clearly identified. Finally, the role of density in the definition of the mixed layer's differing behaviors between the oriental and the occidental regions of the basin is presented.

The year-to-year variability in the timing, duration, and spatial extent of the surface phytoplankton bloom over the winter-spring period is examined in the southern Adriatic Sea using Sea-viewing Wide Field-of-view Sensor (SeaWiFS)-derived chlorophyll images for three years (1998, 1999, and 2000). Each year's image time series shows that blooms were intermittent and differed in onset, duration, and intensity with relatively low values observed in 2000. The relation between atmospheric forcing and interannual variability of the bloom timing and intensity is investigated using a coupled physical-biological model. The simulations focus on the effect of cumulative buoyancy loss on convective depths and its implications on surface nutrient availability, chlorophyll concentrations, and other ecosystem components during the study period. We test the hypothesis that the south Adriatic bloom is essentially controlled by the local winter climatic conditions (i.e., maximum convective depth), as suggested by recent findings, rather than the available nutrient pool at the intermediate depths (200-800 m), which also varies from year to year. For all three years the simulations produced convective depths that were in good agreement with in situ observations. However, the fluctuations in SeaWiFS phytoplankton biomass could be reproduced only if the particular year's nutrient pool was also taken into account. Thus the most probable explanation for the low SeaWiFS phytoplankton biomass observed in 2000 is the reduced nutrient pool because of the return from the transient phase to the pretransient regime of the Mediterranean Sea. Our results indicate that the south Adriatic bloom is a complex phenomenon and cannot simply be explained by interannual changes in convective depth.

Two years and six months of night-time Advanced Very High Resolution Radiometer (AVHRR) sea surface temperature (SST) and daytime Sea viewing Wide Field of view Sensor (SeaWiFS) data collected during the MFSPP have been used to examine spatial and temporal variability of SST and chlorophyll (Chl) in the Adriatic Sea. Flows along the Albanian and the Italian coasts can be distinguished year-round in the monthly averaged Chl but only in the colder months in the monthly averaged SST's. The Chl monthly-averaged fields supply less information on circulation features away from coastal boundaries and where conditions are generally oligotrophic, except for the early spring bloom in the Southern Adriatic Gyre. To better characterise the year-to-year and seasonal variability, exploratory data analysis techniques, particularly the plotting of multiple Chl-SST histograms, are employed to make joint quantitative use of monthly-averaged fields. Modal water mass (MW), corresponding to the Chl-SST pairs in the neighbourhood of the maximum of each monthly histogram, are chosen to represent the temporal and spatial evolution of the prevalent processes and their variability in the Adriatic Sea. Over an annual cycle, the MW followed a triangular path with the most pronounced seasonal and interannual variations in both Chl-SST properties and spatial distributions of the MW in the colder part of the year. The winter of 1999 is the colder (by at least 0.5°C) and most eutrophic (by 0.2 mg/m³). The fall of the year 2000 is characterised by the lack of cooling in the month of November that was observed in the previous year. In addition to characterising the MW, the two-dimensional histogram technique allows a distinction to be made between different months in terms of the spread of SST values at a given Chl concentration. During spring and summer, the spread is minimal indicating surface homothermal conditions. In fall and winter, on the other hand, a spread of points suggesting a linear negative correlation between SST and Chl is found. This behaviour is related to the high nutrient content of cooler water associated with upwelling or the Po River fresh water outflow. <b>Key words. </b>Oceanography: general (diurnal, seasonal and annual cycles; marginal and semi-enclosed seas; water masses)