This study presents a model of chlorophyll and primary production in the pelagic Mediterranean Sea. A 3-D-biogeochemical model (OPATM-BFM) was adopted to explore specific system characteristics and quantify dynamics of key biogeochemical variables over a 6 yr period, from 1999 to 2004. We show that, on a basin scale, the Mediterranean Sea is characterised by a high degree of spatial and temporal variability in terms of primary production and chlorophyll concentrations. On a spatial scale, important horizontal and vertical gradients have been observed. According to the simulations over a 6 yr period, the developed model correctly simulated the climatological features of deep chlorophyll maxima and chlorophyll west-east gradients, as well as the seasonal variability in the main offshore regions that were studied. The integrated net primary production highlights north-south gradients that differ from surface net primary production gradients and illustrates the importance of resolving spatial and temporal variations to calculate basin-wide budgets and their variability. According to the model, the western Mediterranean, in particular the Alboran Sea, can be considered mesotrophic, whereas the eastern Mediterranean is oligotrophic. During summer stratified period, notable differences between surface net primary production variability and the corresponding vertically integrated production rates have been identified, suggesting that care must be taken when inferring productivity in such systems from satellite observations alone. Finally, specific simulations that were designed to explore the role of external fluxes and light penetration were performed. The subsequent results show that the effects of atmospheric and terrestrial nutrient loads on the total integrated net primary production account for less than 5 % of the its annual value, whereas an increase of 30 % in the light extinction factor impacts primary production by approximately 10 %.

We present an analysis of the spatial distributions and seasonal dynamics of phosphate and nitrate in the Mediterranean Sea. The analysis involves observations under present conditions and relies on statistical analysis of the output of a state-of-the-art corroborated numerical model. Model results are corroborated with available data from cruises and fixed monitoring stations. The results capture and quantitatively detail the main spatial patterns observed in the Mediterranean Sea, with an east–west gradient in the surface layer for phosphate and a north–south gradient for nitrate. The model results highlight that plankton growth in Mediterranean waters is usually nutrient limited and that such limitation is less frequent in the western sub-basin and is usually due to a lack of phosphate, rather than nitrate, with the exception of the areas close to the Atlantic Ocean, where nitrogen limitation is important. The model simulation illustrates the benefit of a variable stoichiometric formulation of the phytoplankton cell, showing how the ratio of the internal N:P quota adapts to different limiting conditions

Ocean and marine ecosystems provide a range of valuable services to humans, including benefits such as carbon sequestration, whose economic value are as yet poorly understood. This paper presents a novel contribution to the valuation of carbon sequestration services in marine ecosystems with an application to the Mediterranean Sea. We combine a state-of-the-art biogeochemical model with various estimates of the social cost of carbon emissions to provide a spatially explicit characterization of the current flow of values that are attributable to the various sequestration processes, including the biological component. Using conservative estimates of the social cost of carbon, we evaluate the carbon sequestration value flows over the entire basin to range between 127 and 1722 million €/year. Values per unit area range from −135 to 1000 €/km2 year, with the exclusive economic zone of some countries acting as net carbon sources. Whereas the contribution of physical processes can be either positive or negative, also depending on the properties of incoming Atlantic water, the contribution of biological processes to the marine “blue carbon” sequestration is always positive, and found to range between 100 to 1500 million €/year for the whole basin

The paper provides a basin-scale assessment of the spatiotemporal distribution of alkalinity in the Mediterranean Sea. The assessment is made by integrating the available observations into a 3-D transport–biogeochemical model. The results indicate the presence of complex spatial patterns: a marked west-to-east surface gradient of alkalinity is coupled to secondary negative gradients: (1) from marginal seas (Adriatic and Aegean Sea) to the eastern Mediterranean Sea and (2) from north to south in the western region. The west–east gradient is related to the mixing of Atlantic water entering from the Strait of Gibraltar with the high-alkaline water of the eastern sub-basins, which is correlated to the positive surface flux of evaporation minus precipitation. The north-to-south gradients are related to the terrestrial input and to the input of the Black Sea water through the Dardanelles. In the surface layers, alkalinity has a relevant seasonal cycle (up to 40 μmol kg−1) that is driven by physical processes (seasonal cycle of evaporation and vertical mixing) and, to a minor extent, by biological processes. A comparison of alkalinity vs. salinity indicates that different regions present different relationships: in regions of freshwater influence, the two quantities are negatively correlated due to riverine alkalinity input, whereas they are positively correlated in open sea areas of the Mediterranean Sea.

The projected warming, nutrient decline, changes in net primary production, deoxygenation and acidification of the global ocean will affect marine ecosystems during the 21st century. Here, the climate change-related impacts on the marine ecosystems of the Mediterranean Sea in the middle and at the end of the 21st century are assessed using high-resolution projections of the physical and biogeochemical state of the basin under Representative Concentration Pathways (RCPs) 4.5 and 8.5. In both scenarios, the analysis shows changes in the dissolved nutrient contents of the euphotic and intermediate layers of the basin, net primary production, phytoplankton respiration and carbon stock (including phytoplankton, zooplankton, bacterial biomass and particulate organic matter). The projections also show uniform surface and subsurface reductions in the oxygen concentration driven by the warming of the water column and by the increase in ecosystem respiration as well as an acidification signal in the upper water column linked to the increase in the dissolved inorganic carbon content of the water column due to CO₂ absorption from the atmosphere and the increase in respiration. The projected changes are stronger in the RCP8.5 (worst-case) scenario and, in particular, in the eastern Mediterranean due to the limited influence of the exchanges in the Strait of Gibraltar in that part of the basin. On the other hand, analysis of the projections under the RCP4.5 emission scenario shows a tendency to recover the values observed at the beginning of the 21st century for several biogeochemical variables in the second half of the period. This result supports the idea – possibly based on the existence in a system such as the Mediterranean Sea of a certain buffer capacity and renewal rate – that the implementation of policies for reducing CO₂ emission could indeed be effective and could contribute to the foundation of ocean sustainability science and policies.

The quality of the upgraded version of the Copernicus Marine Environment Monitoring Service (CMEMS) biogeochemical operational system of the Mediterranean Sea (MedBFM) is assessed in terms of consistency and forecast skill, following a mixed validation protocol that exploits different reference data from satellite, oceanographic databases, Biogeochemical Argo floats, and literature. We show that the quality of the MedBFM system has been improved in the previous 10 years. We demonstrate that a set of metrics based on the GODAE (Global Ocean Data Assimilation Experiment) paradigm can be efficiently applied to validate an operational model system for biogeochemical and ecosystem forecasts. The accuracy of the CMEMS biogeochemical products for the Mediterranean Sea can be achieved from basin-wide and seasonal scales to mesoscale and weekly scales, and its level depends on the specific variable and the availability of reference data, the latter being an important prerequisite to build robust statistics. In particular, the use of the Biogeochemical Argo floats data proved to significantly enhance the validation framework of operational biogeochemical models. New skill metrics, aimed to assess key biogeochemical processes and dynamics (e.g. deep chlorophyll maximum depth, nitracline depth), can be easily implemented to routinely monitor the quality of the products and highlight possible anomalies through the comparison of near-real-time (NRT) forecasts skill with pre-operationally defined seasonal benchmarks. Feedbacks to the observing autonomous systems in terms of quality control and deployment strategy are also discussed.