Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC) are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The global net uptake of CO2 by the ocean (SOCEAN, called the ocean sink) is estimated with global ocean biogeochemistry models and observation-based fCO2 products (fCO2 is the fugacity of CO2). The global net uptake of CO2 by the land (SLAND, called the land sink) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The sum of all sources and sinks results in the carbon budget imbalance (BIM), a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2023, EFOS increased by 1.3 % relative to 2022, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (10.3 ± 0.5 GtC yr−1 when the cement carbonation sink is not included), and ELUC was 1.0 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1 ± 0.9 GtC yr−1 (40.6 ± 3.2 GtCO2 yr−1). Also, for 2023, GATM was 5.9 ± 0.2 GtC yr−1 (2.79 ± 0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 2.3 ± 1.0 GtC yr−1, with a near-zero BIM (−0.02 GtC yr−1). The global atmospheric CO2 concentration averaged over 2023 reached 419.31 ± 0.1 ppm. Preliminary data for 2024 suggest an increase in EFOS relative to 2023 of +0.8 % (−0.2 % to 1.7 %) globally and an atmospheric CO2 concentration increase by 2.87 ppm, reaching 422.45 ppm, 52 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2023, with a near-zero overall budget imbalance, although discrepancies of up to around 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the mean ocean sink. This living-data update documents changes in methods and datasets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2024 (Friedlingstein et al., 2024).

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.

Abstract The ocean will ultimately store most of the CO 2 emitted to the atmosphere by human activities. Despite its importance, estimates of the 2000−2022 trend in the ocean CO 2 sink differ by a factor of two between observation-based products and process-based models. Here we address this discrepancy using a hybrid approach that preserves the consistency of known processes but constrains the outcome using observations. We show that the hybrid approach reproduces the stagnation of the ocean CO 2 sink in the 1990s and its reinvigoration in the 2000s suggested by observation-based products and matches their amplitude. It suggests that process-based models underestimate the amplitude of the decadal variability in the ocean CO 2 sink, but that observation-based products on average overestimate the decadal trend in the 2010s. The hybrid approach constrains the 2000−2022 trend in the ocean CO 2 sink to 0.42 ± 0.06 Pg C yr −1 decade −1 , and by inference the total land CO 2 sink to 0.28 ± 0.13 Pg C yr −1 decade −1 .

The Southern Ocean is a major sink of atmospheric CO 2, but the nature and magnitude of its variability remains uncertain and debated. Estimates based on observations suggest substantial variability that is not reproduced by process-based ocean models, with increasingly divergent estimates over the past decade. We examine potential constraints on the nature and magnitude of climate-driven variability of the Southern Ocean CO 2 sink from observation-based air-sea O 2 fluxes. On interannual time scales, the variability in the air-sea fluxes of CO 2 and O 2 estimated from observations is consistent across the two species and positively correlated with the variability simulated by ocean models. Our analysis suggests that variations in ocean ventilation related to the Southern Annular Mode are responsible for this interannual variability. On decadal time scales, the existence of significant variability in the air-sea CO 2 flux estimated from observations also tends to be supported by observation-based estimates of O 2 flux variability. However, the large decadal variability in air-sea CO 2 flux is absent from ocean models. Our analysis suggests that issues in representing the balance between the thermal and non-thermal components of the CO 2 sink and/or insufficient variability in mode water formation might contribute to the lack of decadal variability in the current generation of ocean models. This article is part of a discussion meeting issue 'Heat and carbon uptake in the Southern Ocean: the state of the art and future priorities'.

The Southern Ocean plays a critical role in the uptake, transport, and storage of carbon by the global oceans. It is the ocean's largest sink of CO 2 , yet it is also among the regions with the lowest storage of anthropogenic carbon. This behavior results from a unique combination of high winds driving the upwelling of deep waters and the subduction and northward transport of surface carbon. Here we isolate the direct effect of increasing anthropogenic CO 2 in the atmosphere from the indirect effect of climate variability and climate change on the reorganization of carbon in the Southern Ocean interior using a combination of modeling and observations. We show that the effect of climate variability and climate change on the storage of carbon in the Southern Ocean is nearly as large as the effect of anthropogenic CO 2 during the period 1998–2018 compared with the climatology around the year 1995. We identify a distinct climate fingerprint in dissolved inorganic carbon (DIC), with elevated DIC concentration in the ocean at 300–600 m that reinforces the anthropogenic CO 2 signal, and reduced DIC concentration in the ocean around 2,000 m that offsets the anthropogenic CO 2 signal. The fingerprint is strongest at lower latitudes (30°–55°S). This fingerprint could serve to monitor the highly uncertain evolution of carbon within this critical ocean basin, and better identify its drivers.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land-use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based fCO2 products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. Additional lines of evidence on land and ocean sinks are provided by atmospheric inversions, atmospheric oxygen measurements, and Earth system models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and incomplete understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2022, EFOS increased by 0.9 % relative to 2021, with fossil emissions at 9.9±0.5 Gt C yr−1 (10.2±0.5 Gt C yr−1 when the cement carbonation sink is not included), and ELUC was 1.2±0.7 Gt C yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 11.1±0.8 Gt C yr−1 (40.7±3.2 Gt CO2 yr−1). Also, for 2022, GATM was 4.6±0.2 Gt C yr−1 (2.18±0.1 ppm yr−1; ppm denotes parts per million), SOCEAN was 2.8±0.4 Gt C yr−1, and SLAND was 3.8±0.8 Gt C yr−1, with a BIM of −0.1 Gt C yr−1 (i.e. total estimated sources marginally too low or sinks marginally too high). The global atmospheric CO2 concentration averaged over 2022 reached 417.1±0.1 ppm. Preliminary data for 2023 suggest an increase in EFOS relative to 2022 of +1.1 % (0.0 % to 2.1 %) globally and atmospheric CO2 concentration reaching 419.3 ppm, 51 % above the pre-industrial level (around 278 ppm in 1750). Overall, the mean of and trend in the components of the global carbon budget are consistently estimated over the period 1959–2022, with a near-zero overall budget imbalance, although discrepancies of up to around 1 Gt C yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows the following: (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living-data update documents changes in methods and data sets applied to this most recent global carbon budget as well as evolving community understanding of the global carbon cycle. The data presented in this work are available at https://doi.org/10.18160/GCP-2023 (Friedlingstein et al., 2023).

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1 (40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with a BIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

The current average state of Red Sea phytoplankton phenology needs to be resolved in order to study future variations that could be induced by climate change. Moreover, a regionalization of the Red Sea could help to identify areas of interest and guide in situ sampling strategies. Here, a clustering method used 21 years of satellite surface chlorophyll‐a concentration observations to characterize similar regions of the Red Sea. Four relevant phytoplankton spatiotemporal patterns (i.e., bio‐regions) were found and linked to biophysical interactions occurring in their respective areas. Two of them, located in the northern part the Red Sea, were characterized by a distinct winter‐time phytoplankton bloom induced by mixing events or associated with a convergence zone. The other two, located in the southern regions, were characterized by phytoplankton blooms in summer and winter which might be under the influence of water advected into the Red Sea from the Gulf of Aden in response to the seasonal monsoon winds. Some observed inter‐annual variabilities in these bio‐regions suggested that physical mechanisms could be highly variable in response to variations in air‐sea heat fluxes and ENSO phases in the northern and southern half of the Red Sea, respectively. This study reveals the importance of sustaining in situ measurements in the Red Sea to build a full understanding about the physical processes that contribute to phytoplankton production in this basin.

The oceanographic and ecological characteristics of phytoplankton dynamics in the Ligurian Sea had been described in the 1960s from transects across the area and related to some physical–biogeochemical processes. This chapter provides an overview of the micro‐, nano‐ and picophytoplankton species. Recurrent estimations of annual primary production cycles in the euphotic layer of the central Ligurian Sea were obtained at the DYFAMED sampling station. Associated with general annual cycle of phytoplankton dynamics, small spatiotemporal processes occur and may play a lead role at local scales. The biogeochemistry of the Ligurian Sea shares important characteristics with other temperate regions of the ocean, such as the North Atlantic. Modeling studies of this complex system has paved the way for the understanding of the coupling between the physics of the ocean and the response of the biogeochemistry in environments that are prone to strong seasonal variations of the surface mixed layer and to strong mesoscale dynamics.

The growth of phytoplankton at high latitudes was generally thought to begin in open waters of the marginal ice zone once the highly reflective sea ice retreats in spring, solar elevation increases, and surface waters become stratified by the addition of sea-ice melt water. In fact, virtually all recent large-scale estimates of primary production in the Arctic Ocean (AO) assume that phytoplankton production in the water column under sea ice is negligible. However, over the past two decades, an emerging literature showing significant under-ice phytoplankton production on a pan-Arctic scale has challenged our paradigms of Arctic phytoplankton ecology and phenology. This evidence, which builds on previous, but scarce reports, requires the Arctic scientific community to change its perception of traditional AO phenology and urgently revise it. In particular, it is essential to better comprehend, on small and large scales, the changing and variable icescapes, the under-ice light field and biogeochemical cycles during the transition from sea-ice covered to ice-free Arctic waters. Here, we provide a baseline of our current knowledge of under-ice blooms (UIBs), by defining their ecology and their environmental setting, but also their regional peculiarities (in terms of occurrence, magnitude, and assemblages), which is shaped by a complex AO. To this end, a multidisciplinary approach, i.e., combining expeditions and modern autonomous technologies, satellite, and modeling analyses, has been used to provide an overview of this pan-Arctic phenological feature, which will become increasingly important in future marine Arctic biogeochemical cycles.

Climate model projections suggest a substantial decrease of sea ice export into the outflow areas of the Arctic Ocean over the 21st century. Fram Strait, located in the Greenland Sea sector, is the principal gateway for ice export from the Arctic Ocean. The consequences of lower sea ice flux through Fram Strait on ocean dynamics and primary production in the Greenland Sea remain unknown. By using the most recent 16 years (2003-2018) of satellite imagery available and hydrographic in situ observations, the role of exported Arctic sea ice on water column stratification and phytoplankton production in the Greenland Sea is evaluated. Years with high Arctic sea ice flux through Fram Strait resulted in high sea ice concentration in the Greenland Sea, stronger water column stratification, and an earlier spring phytoplankton bloom associated with high primary production levels. Similarly, years with low Fram Strait ice flux were associated with a weak water column stratification and a delayed phytoplankton spring bloom. This work emphasizes that sea ice and phytoplankton production in subarctic "outflow seas" can be strongly influenced by changes occurring in the Arctic Ocean.

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.

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.

In seasonal ice zones (SIZs), such as the one of the Greenland Sea, the sea ice growth in winter and subsequent melting in summer influence the phytoplankton activity. However, studies assessing phytoplankton activities over complete annual cycles and at a fine temporal resolution are lacking in this environment. Biogeochemical‐Argo floats, which are able to sample under the ice, were used to collect physical and biogeochemical data along vertical profiles and at 5‐day resolution during two complete annual cycles in the Greenland Sea SIZ. Three phytoplankton activity phases were distinct within an annual cycle: one under ice, a second at the ice edge, and a third one around an open‐water subsurface chlorophyll maximum. As expected, the light and nitrate availabilities controlled the phytoplankton activity and the establishment of these phases. On average, most of the annual net community production occurred equally under ice and at the ice edge. The open‐water subsurface chlorophyll maximum phase contribution, on the other hand, was much smaller. Phytoplankton biomass accumulation and production thus occur over a longer period than might be assumed if under ice blooms were neglected. This also means that satellite‐based estimates of phytoplankton biomass and production in this SIZ are likely underestimated. Simulations with the Arctic‐based physical‐biologically coupled SINMOD model suggest that most of the annual net community production in this SIZ results from local processes rather than due to advection of nitrate from the East Greenland and Jan Mayen Currents.

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.

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.

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.

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.

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.

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.

Le phytoplancton est un élément primordial dans les réseaux trophiques marins et il est un acteur principal dans les cycles biogéochimiques de la planète. Cependant, des incertitudes subsistent autour des facteurs environnementaux influençant sa saisonnalité ainsi que sa capacité à se développer. L’objectif majeur de cette thèse est d’étudier la réponse du phytoplancton à la variabilité interannuelle des facteurs environnementaux en Mer Méditerranée. Plus précisément, il s’agit de déterminer l’influence de ces derniers sur la saisonnalité du phytoplancton.Dans un premier temps, la variabilité interannuelle des cycles annuels de biomasses phytoplanctoniques observables en Méditerranée a été analysée. Certaines régions, tel que les zones de formation d’eau dense, présentent une variabilité interannuelle importante. L’une des régions les plus variables est la zone de formation d’eau dense en Méditerranée Nord-Occidentale. Une approche multi-outils basée sur des observations a été mise en place pour l’étude des variations spatiale et temporelle de la saisonnalité du phytoplancton dans cette région. Le rôle crucial du mélange vertical et de la disponibilité en lumière sur la saisonnalité du phytoplancton a été évalué. Il est démontré qu’une couche de mélange profonde pendant l’hiver augmente l’intensité du bloom phytoplanctonique printanier, due à une présence plus importante dans la communauté phytoplanctonique de micro-phytoplancton. En conséquence, le taux de production primaire printanier augmente. Enfin, ces modifications de la communauté phytoplanctonique et de la production provoquent une augmentation du stock de carbone organique produit au printemps.

Planktonic organisms play crucial roles in oceanic food webs and global biogeochemical cycles1, 2. Most of our knowledge about the ecological impact of large zooplankton stems from research on abundant and robust crustaceans, and in particular copepods3, 4. A number of the other organisms that comprise planktonic communities are fragile, and therefore hard to sample and quantify, meaning that their abundances and effects on oceanic ecosystems are poorly understood. Here, using data from a worldwide in situ imaging survey of plankton larger than 600 μm, we show that a substantial part of the biomass of this size fraction consists of giant protists belonging to the Rhizaria, a super-group of mostly fragile unicellular marine organisms that includes the taxa Phaeodaria and Radiolaria (for example, orders Collodaria and Acantharia). Globally, we estimate that rhizarians in the top 200 m of world oceans represent a standing stock of 0.089 Pg carbon, equivalent to 5.2% of the total oceanic biota carbon reservoir5. In the vast oligotrophic intertropical open oceans, rhizarian biomass is estimated to be equivalent to that of all other mesozooplankton (plankton in the size range 0.2–20 mm). The photosymbiotic association of many rhizarians with microalgae may be an important factor in explaining their distribution. The previously overlooked importance of these giant protists across the widest ecosystem on the planet6 changes our understanding of marine planktonic ecosystems.