Articles | Volume 14, issue 1
Geosci. Model Dev., 14, 295–321, 2021
Geosci. Model Dev., 14, 295–321, 2021

Model description paper 19 Jan 2021

Model description paper | 19 Jan 2021

Implementation and assessment of a carbonate system model (Eco3M-CarbOx v1.1) in a highly dynamic Mediterranean coastal site (Bay of Marseille, France)

Implementation and assessment of a carbonate system model (Eco3M-CarbOx v1.1) in a highly dynamic Mediterranean coastal site (Bay of Marseille, France)
Katixa Lajaunie-Salla1, Frédéric Diaz1, Cathy Wimart-Rousseau1, Thibaut Wagener1, Dominique Lefèvre1, Christophe Yohia2, Irène Xueref-Remy3, Brian Nathan3, Alexandre Armengaud4, and Christel Pinazo1 Katixa Lajaunie-Salla et al.
  • 1Aix Marseille Univ., Université de Toulon, CNRS, IRD, MIO, UM 110, 13288, Marseille, France
  • 2Aix Marseille Univ., CNRS, IRD, OSU Institut Pythéas, 13288, Marseille, France
  • 3Aix Marseille Univ., Université d'Avignon, CNRS, IRD, IMBE, Marseille, France
  • 4AtmoSud: Observatoire de la qualité de l'air en région Sud Provence Alpes Côte d'Azur, le Noilly Paradis, 146 rue Paradis, 13294 Marseille, France

Correspondence: Katixa Lajaunie-Salla ( and Fredéric Diaz (


A carbonate chemistry balance module was implemented into a biogeochemical model of the planktonic food web. The model, named Eco3M-CarbOx, includes 22 state variables that are dispatched into 5 compartments: phytoplankton, heterotrophic bacteria, detrital particulate organic matter, labile dissolved organic, and inorganic matter. This model is applied to and evaluated in the Bay of Marseille (BoM, France), which is a coastal zone impacted by the urbanized and industrialized Aix–Marseille Metropolis, and subject to significant increases in anthropogenic emissions of CO2.

The model was evaluated over the year 2017, for which in situ data of the carbonate system are available in the study site. The biogeochemical state variables of the model only change with time, to represent the time evolution of a sea surface water cell in response to the implemented realistic forcing conditions. The model correctly simulates the value ranges and seasonal dynamics of most of the variables of the carbonate system except for the total alkalinity. Several numerical experiments were conducted to test the response of carbonate system to (i) a seawater temperature increase, (ii) wind events, (iii) Rhône River plume intrusions, and (iv) different levels of atmospheric CO2 contents. This set of numerical experiments shows that the Eco3M-CarbOx model provides expected responses in the alteration of the marine carbonate balance regarding each of the considered perturbation. When the seawater temperature changes quickly, the behavior of the BoM waters alters within a few days from a source of CO2 to the atmosphere to a sink into the ocean. Moreover, the higher the wind speed is, the higher the air–sea CO2 gas exchange fluxes are. The river intrusions with nitrate supplies lead to a decrease in the pCO2 value, favoring the conditions of a sink for atmospheric CO2 into the BoM. A scenario of high atmospheric concentrations of CO2 also favors the conditions of a sink for atmospheric CO2 into the waters of the BoM. Thus the model results suggest that external forcings have an important impact on the carbonate equilibrium in this coastal area.

1 Introduction

Current climate change mostly originates from the carbon dioxide (CO2) increase in the atmosphere at a high annual rate (+2.63 ppm from May 2018 to May 2019,, last access: September 2019). This atmospheric CO2 increase impacts the carbonate chemistry equilibrium of the oceanic water column (Allen et al., 2009; Matthews et al., 2009). Oceans are known to act as a sink for anthropogenic CO2, i.e., 30 % of emissions, which leads to a marine acidification (Gruber et al., 2019; Orr et al., 2005; Le Quéré et al., 2018). CO2 is a key molecule in the biogeochemical functioning of the marine ecosystem. Photo-autotrophic organisms, mainly phytoplankton and macro-algae, fix this gas through photosynthesis in the euphotic zone and, in turn, produce organic matter and dissolved oxygen. Heterotrophic organisms, mainly heterotrophic protists and metazoans, consume organic matter and dissolved oxygen by aerobic respiration and, in turn, produce CO2. In the ocean, the main processes regulating CO2 exchanges between the atmosphere and sea are the solubility pump and the biological pump at different timescales. Overall, the thermohaline gradients drive the solubility pump, while the metabolic processes of gross primary production and respiration set the intensity of the biological pump (Raven and Falkowski, 1999).

The coastal zones, despite their small surface area and volume compared to those of the open ocean, have a large influence upon carbon dynamics and represent 14 % to 30 % of the oceanic primary production (Gattuso et al., 1998). At the interface between open ocean and continents, these zones receive large inputs of nutrients and organic matter from rivers, groundwater discharge, and from atmospheric depositions (Cloern et al., 2014; Gattuso et al., 1998). On coasts, shorelines are subject to an increasing density of population and associated urbanization (Small and Nicholls, 2003). This rapid alteration of shorelines all over the world accelerates the emissions of greenhouses gases near the coastal ocean, and it also involves large discharges of material into the seawater by wastewater runoff and/or rivers (Cloern, 2001). These anthropogenic forcings alter the biogeochemical functioning of these zones and could lead to a growing eutrophication (Cloern, 2001). Moreover, these forcings could affect the carbonate chemistry dynamics and amplify or attenuate the acidification in coastal zones. This alteration of the marine environment may provoke further changes in the structure of the plankton community, including, in fine, consequences on the populations with high trophic levels, such as teleosts (Esbaugh et al., 2012). At the global scale, coastal zones are considered to be a significant sink for atmospheric CO2, with an estimated flux converging to 0.2 Pg C yr−1 (Roobaert et al., 2019). However, some studies highlight that the status of these areas as a net sink or source still remains uncertain due to the complexity of the interactions between biological and physical processes, and also due to the lack of in situ measurements (Borges and Abril, 2011; Chen et al., 2013; Chen and Borges, 2009). Moreover, the capacity for coastal zones to absorb atmospheric CO2 resulting from the increasing human pressure also remains poorly known. There are few works which highlight, under future atmospheric CO2 levels, whether coastal zones will become a net sink or a reduced source of CO2 (Andersson and Mackenzie, 2012; Cai, 2011).

The current increase in the atmospheric CO2 partial pressure (pCO2) is slowly shifting the marine carbonate chemistry equilibrium towards increases in the seawater pCO2 and bicarbonate ions (HCO3-) and decreases in pH and carbonate ions (CO32-) (Hoegh-Guldberg et al., 2018). These trends were already described in several coastal and open-ocean locations worldwide (Cai et al., 2011). In a coastal northwestern Mediterranean site, a 10-year time series of in situ measurements highlights a trend of pH decrease and pCO2 increase (Kapsenberg et al., 2017). Low pH values can inhibit the ability of many marine organisms to form the calcium carbonate (CaCO3) used in the making of skeletons and shells (Gattuso et al., 2015). In an extreme case, this shift may promote dissolution of CaCO3 because the water will become under-saturated with respect to CaCO3 minerals (Doney et al., 2009).

The present study is dedicated to the implementation of a carbonate system module into a preexisting biogeochemical model of the planktonic food web. This new model, named Eco3M-CarbOx (v1.1), is then evaluated in a highly dynamic coastal site, i.e., the Bay of Marseille (BoM) in the northwestern Mediterranean Sea. This evaluation is performed on the seasonal dynamics of biogeochemical and carbonate modeled variables against that of the corresponding in situ data available over the year 2017. This study is extended by a fine analysis of the variability of the marine carbonate system (stocks, fluxes) in relation to physical (e.g., wind events, river intrusions, temperature increases, changes in the atmospheric pCO2 levels) and biogeochemical processes (gross primary production (GPP) and respiration, R) in the study site. The BoM is suitable for this kind of study because this coastal area is subject to high emissions of atmospheric CO2 from the nearby urban area, and it also receives effluents from the Aix–Marseille metropolis. In addition, strong wind events (the mistral) regularly occur, which could lead to (i) strong latent heat losses at the surface (Herrmann et al., 2011) and upwelling along the coast with a common consequence of a cooling effect and (ii) Rhône River plume intrusion under specific wind conditions (Fraysse et al., 2013, 2014). In this regional context, many anthropogenic forcings can interact with the dynamics of the carbonate systems. Natural determinants of the composition of the marine planktonic community can also play a crucial role in these dynamics.

Figure 1Schematic diagram of the biogeochemical model Eco3M-CarbOx. Explicit state variables of the model are indicated in continuous-line box or circles except the implicit variable for zooplankton (dotted line box). State variables written in orange are added variables compared to the preexisting biogeochemical model of Fraysse et al. (2013). Arrows represent processes between two state variables. TA: total alkalinity; DIC: dissolved inorganic carbon; CO2: dissolved carbon dioxide; O2: dissolved oxygen; CaCO3: calcium carbonate.


2 Materials and methods

2.1 Numerical model description

The Eco3M-CarbOx biogeochemical model was developed to represent the dynamics of the seawater carbonate system and plankton food web in the BoM. The model was implemented using the Eco3M (Ecological Mechanistic and Modular Modelling) platform (Baklouti et al., 2006). The model structure used is based on an existing model of the plankton ecosystem (Fraysse et al., 2013), including a description of carbon (C), nitrogen (N), and phosphorus (P) marine biogeochemical cycles. The Eco3M-CarbOx model includes 22 prognostic state variables that are split into several compartments: phytoplankton, heterotrophic bacteria, detrital particulate organic matter, labile dissolved organic matter, nutrients (ammonia, nitrate, and phosphate), dissolved oxygen, and carbonate system variables (Fig. 1). In this study, the state variables of the Eco3M-CarbOx model only change over time (i.e., usually termed “model 0D”); they are representative of the time evolution of a sea surface water cell, but this biogeochemical model is not coupled with a hydrodynamic model.

The model presented in this study includes a set of new developments and improvements in the realism of the plankton web structure and process formulations compared to the model of Fraysse et al. (2013). In order to improve the representation of chlorophyll concentration in the Bay of Marseille the phytoplankton is divided into two groups: one with some ecological and physiological traits of the Synechococcus cyanobacteria, which is one of the major constitutive members of pico-autotrophs in the Mediterranean Sea (Mella-Flores et al., 2011), and another with traits of large diatoms, which are generally observed during spring blooms at mid-latitudes (Margalef, 1978). For both of the phytoplankton, there is a diagnostic chlorophyll a variable related to the phytoplankton C biomass, the phytoplankton N-to-C ratio, and the limiting internal ratio fQN (Faure et al., 2010; Smith and Tett, 2000; Table B2, Appendix B). The functional response of primary production was modified using another formulation of the temperature limitation function which takes into account the optimal temperature of growth for each phytoplankton group. The exudation of phytoplankton was modified taking into account the intracellular phytoplankton ratio. For the uptake of matter by bacteria and the remineralization processes the dependence on the intracellular bacteria ratio was modified. A temperature dependence of all biogeochemical processes was added to take into account the effects of rapid and strong variations of seawater temperature on plankton during episodes of upwelling, for instance, that are usually observed in the BoM. Also certain parameters in some formulations were modified owing to the alterations of some formulations (Table B4, Appendix B).

Additionally, a carbonate system module was developed and three state variables were added: dissolved inorganic carbon (DIC), total alkalinity (TA) and the calcium carbonate (CaCO3) implicitly representing calcifying organisms. The knowledge of DIC and TA allows the calculation of the pCO2 and pH (total pH scale) diagnostic variables, necessary for resolving all the equations of the carbonate system. These equations use apparent equilibrium constants, which depend on temperature, pressure, and salinity (Dickson, 1990a, b; Dickson and Riley, 1979; Lueker et al., 2000; Millero, 1995; Morris and Riley, 1966; Mucci, 1983; Riley, 1965; Riley and Tongudai, 1967; Uppström, 1974; Weiss, 1974). The details of the resolution of carbonate system module are given in Appendix A. For this module three processes were also added: the precipitation and dissolution of calcium carbonate and the gas exchange of pCO2 with the atmosphere. Based on the review of Middelburg (2019), it is considered that (i) TA decreases by 2 moles for each mole of CaCO3 precipitated, by 1 mole for each mole of ammonium nitrified, by 1 mole for each mole of ammonium assimilated by phytoplankton, and TA increases by 2 moles for each mole of CaCO3 dissolved, and by 1 mole for each mole of organic matter mineralized by bacteria in ammonium (Table B2, Appendix B); (ii) DIC is consumed during the photosynthesis and calcification processes and is produced by respiration (of phytoplankton, zooplankton, and bacteria) and the CaCO3 dissolution processes. Moreover, the dynamics of DIC are altered by the CO2 exchanges with the atmosphere (Table B2, Appendix B). The air–sea CO2 fluxes are calculated from the pCO2 gradient across the air–sea interface and the gas transfer velocity (Table B3, Appendix B) estimated from the wind speed and using the parametrization of Wanninkhof (1992).

In the Eco3M-CarbOx model, zooplankton is considered as an implicit variable. However, a closure term based on the assumption that all of the matter grazed by the zooplankton and higher trophic levels return as either organic or inorganic matter by excretion, egestion, and mortality processes is taken into account (Fraysse et al., 2013). The model considers a “non-Redfieldian” stoichiometry for phytoplankton and bacteria. All the biogeochemical model formulations, equations, and associated parameter values are detailed in Appendix B.

Figure 2Map of study area: the Region Sud (a), Aix–Marseille Metropolis (area in the red rectangle) (b), the Bay of Marseille (c). CAV = Cinq Avenues Station (urban site); OHP: Observatoire de Haute Provence station (non-urban site); Carry, Solemio, Planier: marine study sites in the Bay of Marseille.

2.2 Study area

The BoM is located in the eastern part of the Gulf of Lion, in the northwestern Mediterranean Sea (Fig. 2). The city of Marseille, located on the coast of the BoM, is the second largest city of France, with a population of ca. 1 million. The Rhône River, which flows into the Gulf of Lion, is the greatest source of freshwater and nutrients for the Mediterranean Sea, with a river mean flow of 1800 m3 s−1 (Pont et al., 2002). Several studies highlight the eastward intrusion events from the Rhône River plume in the BoM under east and southeasterly wind conditions, which favor biological productivity (Fraysse et al., 2014; Gatti et al., 2006; Para et al., 2010). The biogeochemistry of the BoM is complex and highly driven by hydrodynamics. For example, north–northwesterly winds induce upwelling events which bring cold and nutrient-rich waters upward (Fraysse et al., 2013). Moreover, the oligotrophic Northern Current occasionally intrudes into the BoM (Petrenko, 2003; Ross et al., 2016).

Despite the presence of several marine protected areas around the BoM (the Regional Park of Camargue, the Marine protected area Côte Bleue, and the National Park of Calanques), it is strongly impacted by diverse anthropogenic forcing, because industrialized and urbanized areas are located all along the coast. From the land, the BoM receives nutrients and organic matter from the urban area of the Aix–Marseille metropolis (Millet et al., 2018), the industrialized area of Fos-sur-Mer (one of the biggest oil-based industry areas in Europe), and the Berre Lagoon, which is eutrophized (Gouze et al., 2008; Fig. 2c). From the atmosphere, the BoM is subject to fine-particle deposition and greenhouse gas emissions (including CO2) from the nearby urban area, and it also receives effluents from the Aix–Marseille metropolis.

2.3 Dataset

The modeled variables of the carbonate system (DIC, TA, pH, and pCO2) and chlorophyll a are hereafter compared to observations collected at the SOLEMIO station (Figs. 2c and 3), which is a component of the French national monitoring network (Service d'Observation en Milieu Littoral – SOMLIT,, last access: January 2020). Major biogeochemical parameters have been recorded since 1994. Carbonate chemistry variables (pH, pCO2, DIC, and TA), sampled every 2 weeks, have been available since 2016.

Table 1Forcing of the different scenarios (S) simulated with the model. See Sect. 2.4 for details of the scenarios.

Download Print Version | Download XLSX

2.4 Design of numerical experiments

In the present work, the Eco3M-CarbOx model was run for the whole year of 2017. This year was chosen because in situ data of carbonate systems (DIC, TA, pH, and pCO2) are available for the whole year at the SOLEMIO station (Fig. 2c). The biogeochemical variables were initialized using in situ data from winter conditions (Table B1, Appendix B). The model was forced by time series of sea surface temperature and salinity, wind (at 10 m), light, and atmospheric CO2 concentrations. The sea temperature time series is from in situ hourly data recorded at the Planier station (Fig. 2c). For salinity, hourly in situ data from the SOLEMIO station and from the CARRY buoy were used (Fig. 2c). Wind and light hourly time series were extracted from the WRF meteorological model at the SOLEMIO station (Yohia, 2017). Finally, we used hourly atmospheric CO2 values from in situ measurements recorded at the Cinq Avenues station (CAV station; Fig. 2b) by the AtmoSud Regional Atmospheric Survey Network, France (, last access: August 2019). This simulation is the reference simulation (noted S0). As highlighted previously, Rhône River plume intrusions (due to wind-specific conditions) have an impact on the dynamics of primary production (Fraysse et al., 2014; Ross et al., 2016) and then on the seawater carbonate system. Moreover, the seawater temperature and atmospheric CO2 variations control the seawater CO2 dynamics via the solubility equilibrium and gas exchange with the atmosphere (Middelburg, 2019). In order to quantify the impact of different forcing, several simulations (hereafter noted S), which are summarized in Table 1, were conducted:

  • Impact of seawater temperature increase, S1. The forcing time series of in situ temperatures was shifted by +1.5 C (Cocco et al., 2013).

  • Impact of wind events. A first simulation S2 was run with a constant wind intensity of 7 m s−1 (2017 annual average wind speed) throughout the year and a second one (S3) with two three-day periods of strong wind speed (20 m s−1) representative of short bursts of the mistral (data not shown) starting on 15 May and 15 August, and a constant value of 7 m s−1 the rest of the year.

  • Impact of nutrient supply (nitrate) during a Rhône River plume intrusion (S4). A threshold of 37 has been chosen to identify the presence of low-salinity waters from the Rhône River plume in the forcing file of salinity. Here, concentrations in nitrate supplied by the river depend on the salinity. A relationship between these two variables was then established for the SOLEMIO point from the MARS3D-RHOMA coupled physical and biogeochemical model (Fraysse et al., 2013; Pairaud et al., 2011). This relationship has already been used successfully to reproduce realistic observed conditions in the studies of Fraysse et al. (2014) and Ross et al. (2016): NO3intrusion (mmol m−3) =-1.70×S+65.

  • Non-urban atmospheric CO2 concentrations (S5). This simulation takes into account the forcing of atmospheric CO2 values measured at the Observatoire de Haute Provence station (OHP; Fig. 2b), located outside of the Aix–Marseille metropolis, by the ICOS National Network, France (, last access: August 2019).

In this work, we calculated the daily mean values of state variables, statistical parameters, and mean fluxes of modeled processes throughout the year and over two main hydrological periods: the stratified and mixed water column periods. The stratified water column (SWC) is defined with a temperature difference between the surface and bottom of more than 0.5 C (Monterey and Levitus, 1997). For the simulated year (2017), the SWC period lasts from 10 May to 20 October. The mixed water column (MWC) period corresponds to the rest of the year.

3 Results

3.1 Model skills

Following the recommendations of Rykiel (1996), three criteria were considered to evaluate the performance of our model:

  • Does the model reproduce the timing of the observed variations of carbonate system at the seasonal timescale?

  • Does the model reproduce the observed pCO2 and pH ranges at the seasonal timescale?

  • Are the observations correctly represented according to the Willmott skill score (WSS)? This index is an objective measurement of the degree of agreement between the modeled results and the observed data. A correct representation of observations by the model is achieved when this index is higher than 0.70 (Willmott, 1982).

Figure 3Comparison of model results (red) and in situ data (blue) at the surface of the SOLEMIO station. (a) Chlorophyll a concentrations (mg m−3), (b) pCO2 (µatm), (c) pH, (d) DIC (µmol kg−1), (e) TA (µmol kg−1). The value of each state variable represents the mean around ± 5 d of the seawater sampling date. The shaded area and dotted black line delimit the SWC and MWC periods.


Table 2Statistical evaluation of observations vs. model for 2017: observed and simulated minimum and maximum values. WSS = Willmott skill score; N= number of measurements. Units of bias are those of modeled variables: chlorophyll a (Chl a, mg m−3), seawater partial pressure of CO2 (seawater pCO2, µatm), pH, dissolved inorganic carbon (DIC, µmol kg−1), and total alkalinity (TA, µmol kg−1).

* Significant value of WSS (> 0.70).

Download Print Version | Download XLSX

Over most of the studied period, the model simulates lower chlorophyll a concentrations than the in situ observations, especially during the MWC period (Fig. 3a). Two maxima of chlorophyll a concentrations are observed in situ: the first one at ca. 1.71 mg m−3 in March and the second one at ca. 0.68 mg m−3 in May. They are both linked to Rhône River plume intrusions. Several in situ maxima between 0.50 and 0.70 mg m−3 are observed between March and April (at the end of the MWC period), and they signaled the spring bloom event (Table 2 and Fig. 3a). The biogeochemical model quantitatively reproduces the spring bloom observed at the end of the MWC period (Fig. 3a) with a maximum value of ca. 0.69 mg m−3. The model does not catch the two aforementioned maxima of chlorophyll, and it contains a low WSS and a strong bias (0.37 and +0.22 mg m−3, respectively – Table 2).

On the whole, the seasonal variations of the seawater pCO2 are correctly simulated by the biogeochemical model (Fig. 3b), even if the values are rather overestimated during the MWC period. From January to February, the model reproduces the slight decrease in the observed pCO2, and from February to March it reproduces the increase in pCO2 even if the amplitude of this increase remains smaller. In mid-April, during the simulated spring bloom period, the observed drop in pCO2 and increase in pH are also spotted in the model (Fig. 3b and c). The model especially succeeds in reproducing the observed increase in relation to high temperatures during the SWC period. The reduction of the CO2 solubility due to thermal effects mostly explains the increase in pCO2 during the SWC period. The strong SD of modeled values during the SWC period can be explained by the rapid changes in temperature, probably due to upwelling usually occurring at this time of the year (Millot, 1990). The range of modeled pCO2 values (345–503 µatm) encompasses the range of observed values (358–471 µatm; Table 2). The statistical analysis provides a mean bias of +23 µatm, and a WSS of 0.69 (Table 2).

The seasonal dynamic of pH is mostly reproduced by the model, and in particular, the decrease during the SWC period (Fig. 3c). However, the modeled pH is generally underestimated throughout the year, except during the SWC period, with a mean bias of 0.015 (Table 2). The seasonal range is captured by the model with a minimum value during the SWC period (7.994 vs. 8.014 for observations; Table 2) and a maximum one during the MWC period (8.137 vs. 8.114 for observations; Table 2). The statistical analysis highlights an index of agreement between the in situ data and the model outputs higher than 0.70 (Table 2).

The seasonal variations of DIC show the highest values during the MWC period and a decrease (increase) during the beginning (the end) of the SWC period (Fig. 3d). The lowest values are observed during September. The Eco3M-CarbOx model closely matches the seasonal dynamic by reproducing the range of extreme observed values (Table 2). The mean bias is also small (8.48 µmol kg−1; Table 2). More than 70 % (0.73; Table 2) of modeled DIC concentrations are in statistical agreement with the corresponding observations.

The seasonal cycle of measured TA does not show a clear pattern (Fig. 3e). Large variations of values ranging between 2561 and 2624 µmol kg−1 (Table 2) are observed, whatever the hydrological season being considered. The biogeochemical model provides almost constant values of around 2570 µmol kg−1 all throughout the year, which is lower than in situ data. With a low WSS index of agreement and a large mean bias (Table 2), the model is not able to confidently reproduce the observed variations of TA (Fig. 3e and Table 2).

Figure 4(a–f) 2017. (g–l) Temporal focus between 1 April and 31 July 2017. In situ daily average of (a, g) temperature (C, black line) and salinity (blue line) at the SOLEMIO station (at the surface). Modeled daily average (b, h) chlorophyll a concentrations (mg m−3, black line) (c, i) DIC (µmol kg−1, black line) and pH (blue line), (d, j) seawater pCO2 (µatm, black line) and atmosphere pCO2 from OHP (µatm, red line), (e, k) air–sea CO2 fluxes (mmolm-3d-1), (f, l) gross primary production (mmolm-3d-1, green line), total respiration (mmolm-3d-1, red line) and net ecosystem production (mmolm-3d-1, black line). The shaded areas and dotted black lines delimit the SWC and MWC periods.


Table 3Mean flux values (mmolm-3d-1) and the contribution of each process to the DIC variations for the reference simulation over the year and SWC/MWC periods.

GPP: gross primary production, RA: autotrophic respiration, RH: heterotrophic respiration, NEP: net ecosystem production.

Download Print Version | Download XLSX

3.2 Carbon fluxes and budgets

For the year 2017, the values of temperature vary between 13.3 and 25.9 C (Fig. 4a). The DIC variations closely match those of temperature (correlation coefficient 0.75). For example, the spring increase in temperature leads to a decrease in DIC concentrations (Fig. 4a and c), and the minimum values are reached at the end of SWC period. Over the simulated period, the air–sea CO2 fluxes (Faera) vary between 14 and 17 mmolm-3d-1, with a weakly positive annual budget of +6 mmolm-3yr-1 (or +0.017 mmolm-3d-1, Table 3). Then, the BoM waters would act as a net source of CO2 to the atmosphere on an annual basis. However, on a seasonal basis, the BoM waters would change from a net sink during the MWC period (Faera< 0; Table 3) to a net source during the SWC one (Faera> 0; Table 3).

On an annual basis, the gross primary production (GPP) and total respiration (R) are balanced, leading to a null average net ecosystem production (NEP; NEP = GPPR) (Fig. 4f and Table 3). The intensity of autotrophic respiration (Ra) is lower than that of primary production (annual mean of 0.065 vs. 0.413 mmolm-3d-1, respectively – Table 3), while the zooplankton and bacterial respiration account for an average of 0.348 mmolm-3d-1 (Table 3). On a seasonal basis, the model highlights an ecosystem dominated by autotrophy during the MWC period (NEP > 0; Table 3) and heterotrophy during the SWC period with higher flux values (NEP < 0; Table 3). The biogeochemical fluxes show the strongest variations along the SWC period, following those of temperature (Fig. 4f). The maximum GPP occurs in April and is correlated with the maximum chlorophyll concentration. At this time, the ecosystem is autotrophic (NEP > 0; Fig. 4b and f) and is a net sink for atmospheric CO2, which explains the DIC and seawater pCO2 decreases during the bloom period (Fig. 4c–e)

When looking in detail at the 2017 temperature and salinity time series (Fig. 4a), several crucial events can be seen occurring, including freshwater intrusions (e.g., 15 March and 6 May) into the BoM and large variations of temperature in relation to upwelling events or latent heat losses due to wind bursts. The largest freshwater intrusion from the Rhône River plume occurs in mid-March, with a minimum observed salinity of ca. 32.5 at the SOLEMIO station (Fig. 4a). During this event, the seawater pCO2 decreases and pH increases concomitantly (Fig. 4c and d). Then, seawater appears to be temporarily under-saturated in CO2, and the BoM waters thus act as a sink for atmospheric CO2 at the time of intrusion (Fig. 4e).

During the SWC period, upwelling events quickly cool the surface seawater. In two days, from 25 to 27 July, the water temperature drops from 24.7 to 16.9 C (Fig. 4g). The decrease in temperature corresponds to the increase in DIC concentrations (Fig. 4i). Concomitantly, the value of seawater pCO2 decreases from 497 to 352 µatm, and pH increases from 7.99 to 8.12 (Fig. 4i and j). This event quickly changes the BoM waters from a source to a sink for atmospheric CO2 (from +17 to 14 mmolm-3d-1; Fig. 4k), and also from a net heterotrophic to a net autotrophic ecosystem (Fig. 4l).

Figure 5Modeled daily average chlorophyll a concentrations (mg m−3(a), light limitation (b), temperature limitation, and a close-up from 15 April to 5 May of temperature (c) and nutrient limitation (d) for picophytoplankton and the same set for diatoms (f–i). Modeled daily average NEP (mmolm-3d-1e) and air–sea CO2 fluxes (mmolm-3d-1j). Reference simulation (S0, blue line) and temperature-shifted simulation by 1.5 C (S2, red line). The shaded area and dotted black lines delimit the SWC and MWC periods.


3.3 Impact of external forcing on the dynamics of carbonate system

3.3.1 Temperature increase

Here we compare the reference simulation S0 with the S1 simulation (seawater temperature elevation of 1.5 C – Fig. 5). During the year, there are few changes on the carbonate system variables such as the pCO2 and pH (data not shown). The main alterations occur during the blooms of phytoplankton. The simulated bloom of phytoplankton occurs later, at beginning of May, for both diatoms and picophytoplankton, with maximum values of chlorophyll at 1.4 and 0.4 mg m−3, respectively (Fig. 5a and f).

As both the limitations due to light and nutrients remain about the same during the simulations S0 and S1, this counterintuitive occurrence of a bloom relative to changes in temperature is mainly explained by the temperature limiting function, which depends on the optimal temperature of growth (Topt). For the picophytoplankton, from January to April, the increase of 1.5 C drastically reduces the limitation by temperature (Fig. 5c), because the temperature is closer to the optimal temperature (Topt= 16 C, Table A4) during S1 than S0. In the S0 simulation, the temperature reaches Topt ca. 15 April and it induces the bloom, while at the same time in S1 the temperature moves slightly away from Topt and it does not enable the triggering of a bloom. At the time of the bloom in S1, the opposite configuration occurs. In S0, the ambient temperature is again far from Topt, explaining the absence of a bloom, while in the S1 the ambient temperature is closer to Topt, enabling the occurrence of a bloom. The picophytoplankton bloom then occurs later in the warm simulation S1 than in the reference simulation S0 (Fig. 5a). The duration and termination of a bloom is controlled by both the nutrient availability and the temperature (Fig. 5c and d). Inversely, from January to April, the diatoms' growth limitation by temperature is strengthened in the warm simulation S1 (Fig. 5h), because the resulting ambient temperature is further from the optimum temperature (Topt= 13 C, Table A4) than that in the reference simulation S0. This induces a slower growth of diatoms and a delay of the maximum concentration (Fig. 5f). Afterwards the photosynthesis is mainly limited by temperature (Fig. 5h).

The ecosystem is dominated by autotrophs at the time of blooms whatever the simulation considered (NEP > 0; Fig. 5e), and the quantity of DIC (not shown) fixed through autotrophic processes is larger than that released by heterotrophic processes. During the short period of a bloom, the seawater pCO2 decreases, leading to some negative air–sea fluxes (i.e., an oceanic sink for atmospheric CO2). In the warm simulation, the later occurrence of a bloom enables the period of the spring sink to extend by ca. 3 weeks over May relative to the reference simulation (Fig. 5j).

Figure 6Temporal evolution for May (a–c) and August (d–f) 2017 of the wind speed (m s−1, a, d); air–sea CO2 fluxes (mmolm-3d-1, b, e); seawater pCO2 (µatm, c, f). Constant wind scenario (S2, blue line) and wind event scenario (S3, red line). On (c) and (f), the dashed line represents the atmosphere pCO2 (µatm) at the CAV station.


3.3.2 Wind speed

The Bay of Marseille is periodically under the influence of strong wind events (Millot, 1990). Here we compare two simulations: one with a constant wind value (S2) and the other one with two wind events that occur in May and August (S3) (Fig. 6a and d). The result of this numerical experiment shows that the stronger the wind speed is, the higher the air–sea fluxes are, mainly owing to the increase in gas transfer velocity. Depending on the gradient of CO2 between seawater and the atmosphere, strong wind speeds will favor either the emission or uptake of CO2 (Fig. 6b and e). In May, with the air–sea CO2 flux being positive, the outgassing of CO2 to the atmosphere is enhanced, leading to a decrease in seawater pCO2 (Fig. 6c). On the contrary, in August the oceanic sink of atmospheric CO2 is amplified, which leads to an increase in the seawater pCO2 value (Fig. 6f).

Figure 7(a–f) 2017. (g–l) Temporal focus between 1 May and 1 July 2017. (a, g) In situ daily average of salinity. Modeled daily average (b, h) nitrate concentrations (mmol m−3); (c, i) chlorophyll a concentrations (mg m−3); (d, j) DIC (µmol kg−1); (e, k) seawater pCO2 (µatm); and (f, l) pH. Reference simulation (S0, blue line) and nitrate supply simulation (S4, red line). On (e) and (k), the dashed line represents the atmosphere pCO2 (µatm) at the CAV station. The shaded area and dotted black lines delimit the SWC and MWC periods.


3.3.3 Supply in nitrate by river inputs

According to the model results (Fig. 7), the occasional inputs of nitrate (S4) that are linked to Rhône River plume intrusions favor primary production, and they led to increased chlorophyll concentrations (Fig. 7b and c) five times during the SWC period. These blooms, as seen previously, lead to a decrease (increase) in the seawater pCO2 (pH) (Fig. 7e and f). It can be noted that with the strongest river supply at mid-March (Fig. 7a and b) the occurrence of the spring bloom is earlier (Fig. 7c) than that occurring in the reference simulation (S0). The time lag between river nutrient supply and bloom is due to the temperature limitation (Fig. 4c). During blooms occurring within the SWC period following intrusions, the DIC concentrations are generally lower than those of the reference simulation, as in the case of the bloom of mid-May (decrease by ca. 15 µmol kg−1; Fig. 7j), due to the autotrophic processes dominating the heterotrophic ones. In turn, the seawater pCO2 drops by ca. 30 µatm (Fig. 7k) and pH increases by ca. 0.030 (Fig. 7l). Nitrate inputs, favoring primary production, reduce the source of CO2 to the atmosphere and intensify the sink of atmospheric CO2 into the waters of BoM (Fig. 7e and k).

Figure 8(a) Temporal evolution for the year 2017 of the observed pCO2 (µatm) in the atmosphere at the CAV station, called the “urban scenario” (S0, blue line), and at the OHP station, called the “non-urban scenario” (S6, red line). (b) Temporal evolution for the year 2017 of the modeled seawater pCO2 (µatm) with forcings from the urban (S0, blue line) and non-urban (S6, red line) scenarios. The shaded area and dotted black lines delimit the SWC and MWC periods.


3.3.4 Urban air CO2 concentrations

The Aix–Marseille metropolis is strongly subject to urban emissions to the atmosphere (Xueref-Remy et al., 2018a). The seasonal variability of atmospheric CO2 concentrations at the urban site (CAV station; Fig. 2) is much higher than that observed in a non-urban area (OHP station; Fig. 2), especially during the MWC period (Fig. 8a): CO2 concentrations vary between 379 and 547 µatm at the CAV station and between 381 and 429 µatm at the OHP station. Moreover, in winter the atmospheric pCO2 is higher in the urban area than the non-urban area, whereas in summer those of both areas are quite close. These differences in the seasonal pattern and between areas are usually explained by (i) the thinner atmospheric boundary layer, (ii) the decreased fixation of CO2 by terrestrial vegetation, and (iii) the greater influence of anthropogenic activities by emissions from heating (Xueref-Remy et al., 2018b). Forcing the model by atmospheric pCO2 values from urban or non-urban sites can lead to significant differences in the values of the seawater pCO2, especially during the MWC period. The air–sea gradient of pCO2 is higher when using a forcing derived from the CO2 concentrations originating from an urban area than from a non-urban area, which strengthens the sink of atmospheric CO2 into the waters of BoM. The seawater pCO2 is then lower with non-urban area pressure (S5) than with urban area pressure (S0), because of lower CO2 solubility in the BoM (Fig. 8b).

4 Discussion

4.1 Model performance

The evaluation of model skill vs. in situ data highlights that the modeled pH, pCO2, and DIC are in acceptable agreement with observations (Fig. 3). The seasonal variations observed for the different variables are captured by the model, including for example the seasonal decrease in DIC and pH during the SWC period, in relation to the increase in pCO2, and the inverse scenario during the MWC period. The chlorophyll content variability is not well reproduced, especially during spring (Fig. 3a), even taking into account the nitrate supply from the Rhône River plume intrusion (Fig. 7c). This is due to the multiple origins of chlorophyll, organic matter, and nutrients in the BoM that are not accounted for in the Eco3M-CarbOx model: autochthonous marine production, and allochthonous origins from the Rhône and Huveaune River plumes (Fraysse et al., 2013). The observed variations and levels of TA are not correctly simulated by the model (Fig. 3f). The study of Soetaert et al. (2007) highlights that the main variations of TA in the marine coastal zones are linked to freshwater supplies and marine sediments. The present study does not take into account the inputs of TA from the Rhône River and the water–sediment interface, and it may explain why the TA variable is not correctly predicted by our model.

4.2 Contribution of physical and biogeochemical processes to the variability of carbonate system

The contribution of each biogeochemical process to the DIC variability can be assessed using the presented model: the aeration process contributes to 78 % of the DIC variations, and all biogeochemical processes together contribute to 22 % (Table 3). As mentioned by Wimart-Rousseau et al. (2020), the model suggests that the seawater pCO2 variations and associated fluxes would be mostly driven by the seawater temperature dynamics. Moreover, the seasonal variations of the air–sea CO2 flux are in agreement with some previous field studies (De Carlo et al., 2013; Wimart-Rousseau et al., 2020), which measured a weak oceanic sink for atmospheric CO2 during winter and a weak source to the atmosphere during summer.

The model results reveal that temperature would play a crucial role in controlling two counterbalanced processes: (1) the carbonate system equilibrium and (2) the phytoplankton growth. The increase in temperature during SWC leads to a higher pCO2 in seawater due to the decrease in the CO2 solubility (Middelburg, 2019) and, at the same time, the fixation of DIC by phytoplankton is favored, leading to a decrease in the pCO2 level. The imbalance between the latter two processes leads to a change in the ecosystem status (autotrophic or heterotrophic) and the corresponding behavior as a sink or source to the atmosphere. In case of a 1.5 C rise over the whole year, the temperature variation has a very small impact on the carbonate system dynamics. However, it favors the autotrophic processes and strengthens the oceanic sink of atmospheric CO2 during the bloom of phytoplankton (Fig. 5e and j).

4.3 Contribution of the external forcing to the variability of carbonate system

In line with several previous works on the northwestern Mediterranean Sea (De Carlo et al., 2013; Copin-Montégut et al., 2004; Wimart-Rousseau et al., 2020), the model also suggests that the status of the Bay of Marseille regarding sink or source for CO2 could change at high temporal frequency (i.e., hours to days). Bursts of north–northwestern winds lead to sudden and sharp decreases in seawater temperature (< 2 d; Fig. 4g) either directly by latent heat loss through evaporation at the surface (Herrmann et al., 2011) or indirectly by creating upwelling (Millot, 1990), with the consequences of a decrease in the seawater pCO2 values (Fig. 4j) and, in fine, an alteration of the CO2 air–sea fluxes. Model results suggest that the fast variations of temperature could lead to rapid changes of the sink vs. source status in this coastal zone (Fig. 4k). Moreover, Fraysse et al. (2013) highlight that upwelling in the BoM favors ephemeral blooms of phytoplankton by nutrient supplies up to the euphotic layer and would, in turn, contribute to the seawater pCO2 decrease. North and northwestern winds through latent heat losses and/or upwelling events could then enhance the sink for atmospheric CO2 due to the temperature drop and nutrients inputs. However, these results remain preliminary because in our experimental design only the cooling effect of upwelling on the carbonate balance is taken into account. But concomitantly, upwelling usually bring nutrients and DIC at the surface and these supplies could also perturb the balance of the carbonate system. A further coupling of the Eco3M-CarbOx model with a tridimensional hydrodynamic model would certainly enable the multiple effects of upwelling on the dynamics of the carbonate system in this area to be embraced and the results presented in this study to be refined.

High wind speeds (> 7 m s−1) amplified the gaseous exchange of CO2 considerably (De Carlo et al., 2013; Copin-Montégut et al., 2004; Wimart-Rousseau et al., 2020). The model highlights that a strong wind event of 3 d has a significant impact on the seawaterpCO2 values during a longer period of ca. 15 d (Fig. 6). A combination of high atmospheric pCO2 values and high wind speeds would then favor the sink for CO2 into the waters of the BoM. The aeration process also depends on the choice of the formulation of the gas transfer velocity (k600). In this study, the formulation of Wanninkhof (1992) is used and depends of the wind speed at 10 m above the water surface. However, the current velocity could favor the gas exchange, and suspended matter concentration could limit the gas exchange (Abril et al., 2009; Upstill-Goddard, 2006; Zappa et al., 2003). Due to the important heterogeneity of physical and biogeochemical forcings in coastal zones, other factors that control the air–sea gas exchange should certainly be taken into account.

The simulation with intrusions of the Rhône River plume shows that inputs of nitrate cause a drop of seawater pCO2 owing to the nutrient supply favoring the phytoplankton development (Fig. 7). In this scenario, the oceanic sink of atmospheric CO2 is enhanced. But rivers also supply TA (e.g., Gemayel et al., 2015; Schneider et al., 2007) and DIC (e.g., Sempéré et al., 2000) that shift the carbonate system equilibrium toward a pCO2 decrease and a DIC increase (Middelburg, 2019). Taking into account these further supplies may sensibly modify the modeled carbonate balance in the BoM. A next step to the present work will be to design more realistic numerical experiments to refine the results obtained in this preliminary study. The intrusions of the Rhône River plume also induce a salinity decrease in the BoM waters, which leads to a drop in the pCO2 levels in the model. This drop of pCO2 is due to the decrease in the CO2 solubility when salinity decreases (Middelburg, 2019).

In the scenario of forcing the model by using urban atmospheric pCO2 time series, the air–sea gradient increases, and then it enhances the status of the BoM as a sink for atmospheric CO2. As suggested by the in situ study of Wimart-Rousseau et al. (2020), the Eco3M-Carbox model highlights the crucial role of the coastal ocean in an urbanized area, and with an increase in atmospheric CO2, the CO2 uptake by the coastal ocean may increase. This result is in line with the studies of Andersson and Mackenzie (2004) and Cai (2011), which predict an increase in the intensity of the CO2 sink and a potential threat to coastal marine biodiversity in coastal areas owing to high atmospheric CO2 levels.

5 Conclusions

A marine carbonate chemistry module was implemented in the Eco3M-CarbOx biogeochemical model and evaluated against in situ data available in the Bay of Marseille (northwestern Mediterranean Sea) over the year 2017. The model correctly simulates the value ranges and seasonal dynamics of most of the variables of the carbonate system except for the total alkalinity. Several numerical experiments were also conducted to test the sensitivity of the carbon balance to physical processes (temperature and salinity), biogeochemical processes (GPP and respiration processes), and external forcing (wind, river intrusion, and atmospheric CO2). This set of numerical experiments shows that the Eco3M-CarbOx model provides expected responses in the alteration of the marine carbonate balance regarding each of the considered perturbation.

On the whole, the model results suggest that the carbonate system is mainly driven by the seawater temperature dynamics. At a seasonal scale, the BoM marine waters appear to be a net sink of atmospheric CO2 and a dominantly autotrophic ecosystem during the MWC period, and a net source of CO2 to the atmosphere during the SWC period, which is mainly characterized by a dominance of heterotrophic processes. However, the model results highlight that sharp seawater cooling observed within the SWC period, probably owing to upwelling events, cause the CO2 status of the BoM marine waters to change from a source to the atmosphere to a sink into the ocean within a few days. External forcing as the temperature increases leads to a delay in the bloom of phytoplankton. Strong wind events enhance the gas exchange of CO2 with the atmosphere. A Rhône River plume intrusion with input of nitrate favors pCO2 decreases, and the sink of atmospheric CO2 into the BoM waters is enhanced. The higher atmospheric pCO2 values from the urban area intensify the oceanic sink of atmospheric CO2.

The BoM biogeochemical functioning is mainly forced by wind-driven hydrodynamics (upwelling, downwelling), urban rivers, wastewater treatment plants, and atmospheric deposition (Fraysse et al., 2013). In addition, Northern Current and Rhône River plume intrusions frequently occurred (Fraysse et al., 2014; Ross et al., 2016). Moreover, the BoM harbors the second biggest metropolis of France (Marseille), which is impacted by many harbor activities. The next step of this study will be to couple the Eco3M-CarbOx biogeochemical model to a 3-D hydrodynamic model that will mirror the complexity of the BoM functioning. In this way, the contributions of hydrodynamic, atmospheric, anthropic, and biogeochemical processes to the DIC variability will be able to be determined with higher refinement and realism, and an overview of the air–sea CO2 exchange could be made at the scale of the Bay of Marseille. The main results of our study could be transposed to other coastal sites that are also impacted by urban and anthropic pressures. Moreover, in this paper we highlighted that fast and strong variations of pCO2 values occur, and thus it is essential to acquire more in situ values at high frequency (at least with an hourly resolution) to understand the rapid variations of the marine carbon system at these short spatial and temporal scales.

Appendix A: Details of resolution of carbonate system module

A1 Calculation of carbonate system constants

The concentrations in conservative elements (relative to salinity) are calculated as follows:

  • Total fluoride (TF) concentrations from Riley (1965) in mol kg−1:

  • Total sulfate (TS) concentration from Morris and Riley (1966) in mol kg−1:

  • Calcium ion concentration from Riley and Tongudai (1967) in mol kg−1:

  • Total boron (TB) concentration from Uppström (1974) in mol kg−1:

  • Ionic strength (IonS) from Millero (1982):


The constants are calculated on the total pH scale except for KS on free pH scale. If necessary, pH scale conversion factors are as follows:

  • From seawater pH scale (SWS) to total pH scale: SWStoTOT=1+TSKS1+TSKS+TFKF;

  • From free pH scale to total pH scale: FREEtoTOT=1+TSKS.

Further carbonate system constants were calculated as follows:

  • KS equilibrium constant of dissociation of HSO4- from Dickson (1990a) in mol kg−1:

  • KF equilibrium constant of dissociation of hydrogen fluoride (HF) formation from Dickson and Riley (1979) in mol kg−1:

  • KB equilibrium constant of dissociation of boric acid from Dickson (1990b) in mol kg−1:

  • K0 constant of CO2 solubility from Weiss (1974) in molkg-1atm-1:

  • Ke dissociation constant of water from Millero (1995) in (mol kg−1)2:

    Ke=exp(-13847.26T(K)+148.9802-23.6521log(T(K))+(-5.977+118.67T(K)+1.0495log(T(K)))S12-0.01615S),Ke=KeSWStoTOT, on total pH scale in molkg-1;
  • K1 and K2 from Lueker et al. (2000) in mol kg−1:

  • Kca for calcite from Mucci (1983) in (mol kg−1)2:


All the constants are corrected by the effect of hydrostatic pressure:

R=83.1451 in mlbar-1K-1mol-1;


Calculation of the fugacity factor

We suppose that the pressure is at one atmosphere or close to it (Weiss, 1974):

Patm=1.01325bar,delta=57.7-0.118T in cm3mol-1,b=-1636.75+12.0408T-0.0327957T2,+3.165280.00001T3 in cm3mol-1,FugFac=exp(b+2delta)PatmRT.

A2 Resolution of carbonate system

To resolve the carbonate system, we calculate the deltapH, which is the difference of pH between two iterations of the model. We initialize the run by imposing a pH value of 8.

if(nbiter<1)pH=8pHtol=0.001! tolerance for iterations enddeltapH=pHtol+1do while (abs(deltapH)>0.0001)H=10-pHDenom=H2+K1H+K1K2CAlk=DICK1H+2K2DenomBAlk=TBKBKB+HOH=KeHFreetoTot=1+TSKSHfree=HFreetoTotHSO4=TS1+KSHfreeHF=TF1+KFHfreeResidual=TA-CAlk-BAlk-OH+Hfree+HSO4+HFSlope=DICHK1(H2+K1K2+4HK2)Slope=SlopeDenom2+OH+H+BAlkHKB+HSlope=log10SlopedeltapH=Residual/Slope! this is Newton's methoddo while (abs(deltapH)>1)deltapH=deltapH2! to keep the jump from being too bigenddopH=pH+deltapH! Is on the same scale asK1andK2were calculated, i.e., totalpHscalepCO2=DICH2H2+K1H+K1K2106K0FugFac! in µatmCO2=DIC1061+K1H+K1K2H2HCO3=K1CO2H
Appendix B: Biogeochemical model variables and parameters

Table B1Initial conditions of the state variables of the Eco3M-CarbOx model.

* Diagnostic variables.

Download Print Version | Download XLSX

Table B2Balance equations of the Eco3M-CarbOx model.

Download Print Version | Download XLSX

Table B3Biogeochemical processes simulated by the Eco3M-CarbOx model.

Download Print Version | Download XLSX

Table B4Value of parameters.

Download XLSX

Appendix C: Short user manual

After uploading the whole archive on the Zenodo site (Lajaunie-Salla et al., 2020;, the exact version of the Eco3M-CarbOx code used in this study can be run as follows:

make !two executable files will be created: eco3M_ini.exe and eco3M.exe

Please note the following:

  • The file config.ini allows the following to be defined: the time, time step, and save time of simulation variables biogeochemical process.

  • Results files are stocked in “SORTIES” directory.

  • Boundary conditions and forcings data are stocked in “DATA” directory.

  • All subroutines of biogeochemical processes are stocked in “F_PROCESS” directory.

For further information, please contact Frédéric Diaz ( or Christel Pinazo (

Code availability

Eco3M is freely available under the CeCILL license agreement (a French equivalent to the L-GPL license; (last access: 10 February 2020). The Eco3M-CarbOx model is written in Fortran-90/95 and the plotting code is written in MATLAB®. The exact version of the model used to produce the results presented in this paper is archived on Zenodo (Lajaunie-Salla et al., 2020; A short user manual is given in Appendix C of this study.

Author contributions

KLS conceptualized this study, made code developments on the carbonate module, designed the numerical experiments, performed the final analysis to the model results (validation, statistics), developed MATLAB® software to visualize the results of model and statistics tests, wrote and edited the initial draft, revised the manuscript, and answered reviewers. FD participated in the conceptualization of the present study, supervised this research, wrote some parts of the initial draft, helped to answer reviewer comments, revised the draft after the second review, and edited the final version. CWR participated in the data acquisition of the seawater pCO2 and corresponding measurements at the laboratory and reviewed and helped with the initial draft. TW participated in the data acquisition of the seawater pCO2 and corresponding measurements at the laboratory, helped to answer reviews, and helped with the final editing. DL acquired funding, conducted the project administration, and participated in the data acquisition of the seawater pCO2. CY helped to maintain computing resources and gave expertise on the Eco3M-CarbOx code and software developments and compiling. IX-R acquired funding, conducted the project administration, and provided data on the atmospheric pCO2 at OHP. BN participated in the acquisition of pCO2 data at OHP and reviewed and edited the first submitted draft. AA provided temporal data on the atmospheric pCO2 at CAV. CP participated in the conceptualization of the present study, provided the initial Eco3M code (without the carbonate module), supervised this research, reviewed the initial draft, and helped to answer reviewer comments.

Competing interests

The authors declare that they have no conflict of interest.


We thank the National Service d'Observation en MILieu Littoral (SOMLIT) for its permission to use SOLEMIO data. We wish to thank the crew members of the RV Antedon II, operated by the DT-INSU, for making these samplings possible. We wish to acknowledge the team of the SAM platform (Service Atmosphère Mer) of the MIO institute for their help with the field work. For the collection and analyses of the seawater sample, we thank Michel Lafont and Véronique Lagadec of the PACEM (Plateforme Analytique de Chimie des Environnenments Marins) platform of the MIO institute and also the SNAPO-CO2 at LOCEAN, Paris. The SNAPO-CO2 service at LOCEAN is supported by CNRS-INSU and OSU Ecce-Terra.

We acknowledge the staff of the “Cluster de calcul intensif HPC” platform of the OSU Institut PYTHEAS (Aix–Marseille Université, INSU-CNRS) for providing the computing facilities. We gratefully acknowledge Julien Lecubin from the Service Informatique de OSU Institut PYTHEAS (SIP) for their technical assistance. Moreover, we thank Camille Mazoyer and Claire Seceh for their contribution to the Eco3M-CarbOx model development.

Financial support

This study is part of the AMC project (Aix–Marseille Carbon Pilot Study, 2016–2019), funded and performed in the framework of the Labex OT-MED (ANR-11-LABEX-0061, part of the “Investissement d'Avenir” program through the A*MIDEX project ANR-11-IDEX-0001-02), which is funded by the French National Research Agency (ANR). The project leading to this publication has received funding from the European FEDER Fund under project 1166-39417.

Review statement

This paper was edited by Sandra Arndt and reviewed by Guy Munhoven and Sebastiaan van de Velde.


Abril, G., Commarieu, M. V., Sottolichio, A., Bretel, P., and Guérin, F.: Turbidity limits gas exchange in a large macrotidal estuary, Estuar. Coast. Shelf S., 83, 342–348,, 2009. 

Allen, M. R., Frame, D. J., Huntingford, C., Jones, C., Lowe, J. A., Meinshausen, M., and Meinshausen, N.: Warming caused by cumulative carbon emissions towards the trillionth tonne, Nature, 458, 1163–1166,, 2009. 

Andersson, A. J. and Mackenzie, F. T.: Shallow-water oceans: a source or sink of atmospheric CO2?, Front. Ecol. Environ., 2, 348–353,[0348:SOASOS]2.0.CO;2, 2004. 

Andersson, A. J. and Mackenzie, F. T.: Revisiting four scientific debates in ocean acidification research, Biogeosciences, 9, 893–905,, 2012. 

Auger, P. A., Diaz, F., Ulses, C., Estournel, C., Neveux, J., Joux, F., Pujo-Pay, M., and Naudin, J. J.: Functioning of the planktonic ecosystem on the Gulf of Lions shelf (NW Mediterranean) during spring and its impact on the carbon deposition: a field data and 3-D modelling combined approach, Biogeosciences, 8, 3231–3261,, 2011. 

Baklouti, M., Faure, V., Pawlowski, L., and Sciandra, A.: Investigation and sensitivity analysis of a mechanistic phytoplankton model implemented in a new modular numerical tool (Eco3M) dedicated to biogeochemical modelling, Prog. Oceanogr., 71, 34–58,, 2006. 

Borges, A. V. and Abril, G.: Carbon dioxide and methane dynamics in estuaries, in :Treatise on Estuarine and Coastal Science, edited by: Wolanski, E. and McLusky, D., Academic Press, Waltham, 119–161, 2011. 

Cai, W.-J.: Estuarine and coastal ocean carbon paradox: CO2 sinks or sites of terrestrial carbon incineration?, Annu. Rev. Mar. Sci., 3, 123–145,, 2011. 

Cai, W.-J., Hu, X., Huang, W.-J., Murrell, M. C., Lehrter, J. C., Lohrenz, S. E., Chou, W.-C., Zhai, W., Hollibaugh, J. T., Wang, Y., Zhao, P., Guo, X., Gundersen, K., Dai, M., and Gong, G.-C.: Acidification of subsurface coastal waters enhanced by eutrophication, Nat. Geosci., 4, 766–770,, 2011. 

Campbell, R., Diaz, F., Hu, Z., Doglioli, A., Petrenko, A., and Dekeyser, I.: Nutrients and plankton spatial distributions induced by a coastal eddy in the Gulf of Lion. Insights from a numerical model, Prog. Oceanogr., 109, 47–69,, 2013. 

Chen, C.-T. A. and Borges, A. V: Reconciling opposing views on carbon cycling in the coastal ocean: Continental shelves as sinks and near-shore ecosystems as sources of atmospheric CO2, Deep-Sea Res. Pt. II, 56, 578–590,, 2009. 

Chen, C.-T. A., Huang, T.-H., Chen, Y.-C., Bai, Y., He, X., and Kang, Y.: Air–sea exchanges of CO2 in the world's coastal seas, Biogeosciences, 10, 6509–6544,, 2013. 

Cloern, J. E.: Our evolving conceptual model of the coastal eutrophication problem, Mar. Ecol. Prog. Ser., 210, 223–253,, 2001. 

Cloern, J. E., Foster, S. Q., and Kleckner, A. E.: Phytoplankton primary production in the world's estuarine-coastal ecosystems, Biogeosciences, 11, 2477–2501,, 2014. 

Cocco, V., Joos, F., Steinacher, M., Frölicher, T. L., Bopp, L., Dunne, J., Gehlen, M., Heinze, C., Orr, J., Oschlies, A., Schneider, B., Segschneider, J., and Tjiputra, J.: Oxygen and indicators of stress for marine life in multi-model global warming projections, Biogeosciences, 10, 1849–1868,, 2013. 

Copin-Montégut, C., Bégovic, M., and Merlivat, L.: Variability of the partial pressure of CO2 on diel to annual time scales in the Northwestern Mediterranean Sea, Mar. Chem., 85, 169–189,, 2004. 

De Carlo, E. H., Mousseau, L., Passafiume, O., Drupp, P. S., and Gattuso, J.-P.: Carbonate Chemistry and Air–Sea CO2 Flux in a NW Mediterranean Bay Over a Four-Year Period: 2007–2011, Aquat. Geochem., 19, 399–442,, 2013. 

Dickson, A. G.: Standard potential of the reaction: AgCl(s)+ 12H2(g)=Ag(s)+HCl(aq), and and the standard acidity constant of the ion HSO4- in synthetic sea water from 273.15 to 318.15 K, J. Chem. Thermodyn., 22, 113–127,, 1990a. 

Dickson, A. G.: Thermodynamics of the dissociation of boric acid in synthetic seawater from 273.15 to 318.15 K, Deep-Sea Res., 37, 755–766,, 1990b. 

Dickson, A. G. and Riley, J. P.: The estimation of acid dissociation constants in sea-water media from potentiometric titrations with strong base. II. The dissociation of phosphoric acid, Mar. Chem., 7, 101–109,, 1979. 

Doney, S. C., Tilbrook, B., Roy, S., Metzl, N., Le Quéré, C., Hood, M., Feely, R. A., and Bakker, D.: Surface-ocean CO2 variability and vulnerability, Deep-Sea Res. Pt. II, 56, 504–511,, 2009. 

Esbaugh, A. J., Heuer, R., and Grosell, M.: Impacts of ocean acidification on respiratory gas exchange and acid-base balance in a marine teleost, Opsanus beta, J. Comp. Physiol. B, 182, 921–934,, 2012. 

Faure, V., Pinazo, C., Torréton, J.-P., and Jacquet, S.: Modelling the spatial and temporal variability of the SW lagoon of New Caledonia I: A new biogeochemical model based on microbial loop recycling, Mar. Pollut. Bull., 61, 465–479,, 2010. 

Fraysse, M., Pinazo, C., Faure, V., Fuchs, R., Lazzari, P., Raimbault, P., and Pairaud, I.: 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, 1–18,, 2013. 

Fraysse, M., Pairaud, I., Ross, O. N., Faure, V., and Pinazo, C.: Intrusion of Rhône River diluted water into the Bay of Marseille: Generation processes and impacts on ecosystem functioning, J. Geophys. Res.-Oceans, 119, 6535–6556,, 2014. 

Fukuda, R., Ogawa, H., Nagata, T., and Koike, I.: Direct Determination of Carbon and Nitrogen Contents of Natural Bacterial Assemblages in Marine Environments, Appl. Environ. Microb., 64, 3352–3358,, 1998. 

Gatti, J., Petrenko, A., Devenon, J.-L., Leredde, Y., and Ulses, C.: The Rhône river dilution zone present in the northeastern shelf of the Gulf of Lion in December 2003, Cont. Shelf Res., 26, 1794–1805,, 2006. 

Gattuso, J.-P., Frankignoulle, M., and Wollast, R.: Carbon and carbonate metabolism in coastal aquatic ecosystems, Annu. Rev. Ecol. Syst., 29, 405–34,, 1998. 

Gattuso, J.-P., Magnan, A., Bille, R., Cheung, W. W. L., Howes, E. L., Joos, F., Allemand, D., Bopp, L., Cooley, S. R., Eakin, C. M., Hoegh-Guldberg, O., Kelly, R. P., Portner, H.-O., Rogers, A. D., Baxter, J. M., Laffoley, D., Osborn, D., Rankovic, A., Rochette, J., Sumaila, U. R., Treyer, S., and Turley, C.: Contrasting futures for ocean and society from different anthropogenic CO2 emissions scenarios, Science, 349,, 2015. 

Gehlen, M., Gangstø, R., Schneider, B., Bopp, L., Aumont, O., and Ethe, C.: The fate of pelagic CaCO3 production in a high CO2 ocean: a model study, Biogeosciences, 4, 505–519,, 2007. 

Gemayel, E., Hassoun, A. E. R., Benallal, M. A., Goyet, C., Rivaro, P., Abboud-Abi Saab, M., Krasakopoulou, E., Touratier, F., and Ziveri, P.: Climatological variations of total alkalinity and total dissolved inorganic carbon in the Mediterranean Sea surface waters, Earth Syst. Dynam., 6, 789–800,, 2015. 

Gerber, R. P. and Gerber, M. B.: Ingestion of natural particulate organic matter and subsequent assimilation, respiration and growth by tropical lagoon zooplankton, Mar. Biol., 52, 33–43,, 1979. 

Gouze, E., Raimbault, P., Garcia, N., and Picon, P.: Nutrient dynamics and primary production in the eutrophic Berre Lagoon (Mediterranean, France), Transitional Waters Bull., 2, 17–40,, 2008. 

Gruber, N., Clement, D., Carter, B. R., Feely, R. A., van Heuven, S., Hoppema, M., Ishii, M., Key, R. M., Kozyr, A., Lauvset, S. K., Lo Monaco, C., Mathis, J. T., Murata, A., Olsen, A., Perez, F. F., Sabine, C. L., Tanhua, T., and Wanninkhof, R. H.: The oceanic sink for anthropogenic CO2 from 1994 to 2007, Science, 363, 1193–1199,, 2019. 

Gutiérrez-Rodríguez, A., Latasa, M., Scharek, R., Massana, R., Vila, G., and Gasol, J. M.: Growth and grazing rate dynamics of major phytoplankton groups in an oligotrophic coastal site, Estuar. Coast. Shelf S., 95, 77–87,, 2011. 

Harrison, W. G., Harris, L. R., and Irwin, B. D.: The kinetics of nitrogen utilization in the oceanic mixed layer: Nitrate and ammonium interactions at nanomolar concentrations, Limnol. Oceanogr., 41, 16–32,, 1996. 

Herrmann, M., Somot, S., Calmanti, S., Dubois, C., and Sevault, F.: Representation of spatial and temporal variability of daily wind speed and of intense wind events over the Mediterranean Sea using dynamical downscaling: impact of the regional climate model configuration, Nat. Hazards Earth Syst. Sci., 11, 1983–2001,, 2011. 

Hoegh-Guldberg, O., Jacob, D., Taylor, M., Bindi, M., Brown, S., Camilloni, I., Diedhiou, A., Djalante, R., Ebi, K., Engelbrecht, F., Guiot, J., Hijioka, Y., Mehrotra, S., Payne, A., Seneviratne, S. I., Thomas, A., Warren, R., and Zhou, G.: Impacts of 1.5 C Global Warming on Natural and Human Systems, in Global warming of 1.5 C. An IPCC Special Report on the impacts of global warming of 1.5 C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, edited by: Masson-Delmotte, V., Zhai, P., Pörtner, H. O., Roberts, D., Skea, J., Shukla, P. R., Pirani, A., Moufouma-Okia, W., Péan, C., Pidcock, R., Connors, S., Matthews, J. B. R., Chen, Y., Zhou, X., Gomis, M. I., Lonnoy, E., Maycock, T., Tignor, M., and Waterfield, T., World Meteorological Organization Technical Document, 2018. 

Kapsenberg, L., Alliouane, S., Gazeau, F., Mousseau, L., and Gattuso, J.-P.: Coastal ocean acidification and increasing total alkalinity in the northwestern Mediterranean Sea, Ocean Sci., 13, 411–426,, 2017. 

Lacroix, G. and Grégoire, M.: Revisited ecosystem model (MODECOGeL) of the Ligurian Sea: seasonal and interannual variability due to atmospheric forcing, J. Marine Syst., 37, 229–258,, 2002. 

Lajaunie-Salla, K., Diaz, F., and Pinazo, C.: configuration files and code for Eco3M-CarbOx v1.1, Zenodo,, 2020. 

Le Borgne, R.: Zooplankton production in the eastern tropical Atlantic Ocean: Net growth efficiency and P:B in terms of carbon, nitrogen, and phosphorus, Limnol. Oceanogr., 27, 681–698,, 1982. 

Le Borgne, R. and Rodier, M.: Net zooplankton and the biological pump: a comparison between the oligotrophic and mesotrophic equatorial Pacific, Deep-Sea Res. Pt. II, 44, 2003–2023,, 1997. 

Le Quéré, C., Andrew, R. M., Friedlingstein, P., Sitch, S., Hauck, J., Pongratz, J., Pickers, P. A., Korsbakken, J. I., Peters, G. P., Canadell, J. G., Arneth, A., Arora, V. K., Barbero, L., Bastos, A., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Doney, S. C., Gkritzalis, T., Goll, D. S., Harris, I., Haverd, V., Hoffman, F. M., Hoppema, M., Houghton, R. A., Hurtt, G., Ilyina, T., Jain, A. K., Johannessen, T., Jones, C. D., Kato, E., Keeling, R. F., Goldewijk, K. K., Landschützer, P., Lefèvre, N., Lienert, S., Liu, Z., Lombardozzi, D., Metzl, N., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., Neill, C., Olsen, A., Ono, T., Patra, P., Peregon, A., Peters, W., Peylin, P., Pfeil, B., Pierrot, D., Poulter, B., Rehder, G., Resplandy, L., Robertson, E., Rocher, M., Rödenbeck, C., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Steinhoff, T., Sutton, A., Tans, P. P., Tian, H., Tilbrook, B., Tubiello, F. N., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., Wright, R., Zaehle, S., and Zheng, B.: Global Carbon Budget 2018, Earth Syst. Sci. Data, 10, 2141–2194,, 2018. 

Leblanc, K., Quéguiner, B., Diaz, F., Cornet, V., Michel-Rodriguez, M., Durrieu de Madron, X., Bowler, C., Malviya, S., Thyssen, M., Grégori, G., Rembauville, M., Grosso, O., Poulain, J., de Vargas, C., Pujo-Pay, M., and Conan, P.: Nanoplanktonic diatoms are globally overlooked but play a role in spring blooms and carbon export, Nat. Commun., 9, 953,, 2018. 

Lueker, T. J., Dickson, A. G., and Keeling, C. D.: Ocean pCO2 calculated from dissolved inorganic carbon, alkalinity, and equations for K1 and K2: Validation based on laboratory measurements of CO2 in gas and seawater at equilibrium, Mar. Chem., 70, 105–119,, 2000. 

Margalef, R.: Life-forms of phytoplankton as survival alternatives in an unstable environment, edited by: Gauthier-Villars, Oceanol. Acta, 1, 493–509, available at: (last access: December 2019), 1978. 

Marty, J.-C., Chiavérini, J., Pizay, M.-D., and Avril, B.: Seasonal and interannual dynamics of nutrients and phytoplankton pigments in the western Mediterranean Sea at the DYFAMED time-series station (1991–1999), Deep-Sea Res. Pt. II, 49, 1965–1985,, 2002. 

Matthews, H. D., Gillett, N. P., Stott, P. A., and Zickfeld, K.: The proportionality of global warming to cumulative carbon emissions, Nature, 459, 829–32,, 2009. 

Mella-Flores, D., Mazard, S., Humily, F., Partensky, F., Mahé, F., Bariat, L., Courties, C., Marie, D., Ras, J., Mauriac, R., Jeanthon, C., Mahdi Bendif, E., Ostrowski, M., Scanlan, D. J., and Garczarek, L.: Is the distribution of Prochlorococcus and Synechococcus ecotypes in the Mediterranean Sea affected by global warming?, Biogeosciences, 8, 2785–2804,, 2011. 

Middelburg, J. J.: Marine carbon biogeochemistry. A primer for earth system scientists, Springer Briefs in Earth System Sciences, Springer Nature Switzerland AG, Cham, Switzerland, 2019. 

Millero, F. J.: Thermodynamics Seawater – 1. The PVT Properties, Ocean Science and Engineering, 7, 403–460 1982. 

Millero, F. J.: Thermodynamics of the carbon dioxide system in the oceans, Geochim. Cosmochim. Ac., 59, 661–677,, 1995. 

Millet, B., Pinazo, C., Daniela, B., Remi, P., Pierre, G., and Ivane, P.: Unexpected spatial impact of treatment plant discharges induced by episodic hydrodynamic events: Modelling Lagrangian transport of fine particles by Northern Current intrusions in the bays of Marseille (France), PLoS One, 13, e0195257,, 2018. 

Millot, C.: The Gulf of Lions hydrodynamics, Cont. Shelf Res., 10, 885–894,, 1990. 

Monterey, G. and Levitus, S.: Seasonal variability of mixed layer depth for the World Ocean, NOAA Atlas NESDIS 14, U.S. Dept. of Commerce, National Oceanic and Atmospheric Administration, National Environmental Satellite, Data, and Information Service, Washington, D.C., 1997. 

Moran, M. A.: The global ocean microbiome, Sci. Am. Assoc. Adv. Sci., 350,, 2015. 

Morris, A. W. and Riley, J. P.: The bromide/chlorinity and sulphate/chlorinity ratio in sea water, Deep Sea Res. Oceanogr. Abstr., 13, 699–705,, 1966. 

Mucci, A.: The solubility of calcite and aragonite in seawater at various salinities, temperatures, and one atmosphere total pressure, Am. J. Sci., 283, 780–799,, 1983. 

Orr, J. C., Fabry, V. J., Aumont, O., Bopp, L., Doney, S. C., Feely, R. A., Gnanadesikan, A., Gruber, N., Ishida, A., Joos, F., Key, R. M., Lindsay, K., Maier-Reimer, E., Matear, R., Monfray, P., Mouchet, A., Najjar, R. G., Plattner, G.-K., Rodgers, K. B., Sabine, C. L., Sarmiento, J. L., Schlitzer, R., Slater, R. D., Totterdell, I. J., Weirig, M.-F., Yamanaka, Y., and Yool, A.: Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms, Nature, 437, 681–6,, 2005. 

Pairaud, I., Gatti, J., Bensoussan, N., Verney, R., and Garreau, P.: Hydrology and circulation in a coastal area off Marseille: Validation of a nested 3D model with observations, J. Marine Syst., 88, 20–33,, 2011. 

Para, J., Coble, P. G., Charrière, B., Tedetti, M., Fontana, C., and Sempéré, R.: Fluorescence and absorption properties of chromophoric dissolved organic matter (CDOM) in coastal surface waters of the northwestern Mediterranean Sea, influence of the Rhône River, Biogeosciences, 7, 4083–4103,, 2010. 

Petrenko, A.: Variability of circulation features in the Gulf of Lion, NW Mediterranean Sea. Importance of inertial currents, Oceanol. Acta, 26, 323–338,, 2003. 

Pont, D., Simonnet, J.-P., and Walter, A. V.: Medium-term changes in suspended sediment delivery to the Ocean: Consequences of catchment heterogeneity and river management (Rhône River, France), Estuar. Coast. Shelf S., 54, 1–18,, 2002. 

Raven, J. A. and Falkowski, P. G.: Oceanic sinks for atmospheric CO2, Plant Cell Environ., 22, 741–755,, 1999. 

Riley, J. P.: The occurrence of anomalously high fluoride concentrations in the North Atlantic, Deep Sea Research and Oceanographic Abstracts, 12, 219–220,, 1965. 

Riley, J. P. and Tongudai, M.: The major cation/chlorinity ratios in sea water, Chem. Geol., 2, 263–269,, 1967. 

Roobaert, A., Laruelle, G. G., Landschützer, P., Gruber, N., Chou, L., and Regnier, P.: The spatiotemporal dynamics of the sources and sinks of CO2 in the global coastal ocean, Global Biogeochem. Cy., 33,, 2019. 

Ross, O. N., Fraysse, M., Pinazo, C., and Pairaud, I.: Impact of an intrusion by the Northern Current on the biogeochemistry in the eastern Gulf of Lion, NW Mediterranean, Estuar. Coast. Shelf S., 170, 1–9,, 2016. 

Rykiel, E. J.: Testing ecological models: The meaning of validation, Ecol. Model., 90, 229–244,, 1996. 

Sarthou, G., Timmermans, K. R., Blain, S., and Tréguer, P.: Growth physiology and fate of diatoms in the ocean: a review, J. Sea Res., 53, 25–42,, 2005. 

Schneider, A., Wallace, D. W. R., and Körtzinger, A.: Alkalinity of the Mediterranean Sea, Geophys. Res. Lett., 34,, 2007. 

Sempéré, R., Charrìère, B., van Wambeke, F., and Cauwet, G.: Carbon inputs of the Rhône River to the Mediterranean Sea: Biogeochemical implications, Global Biogeochem. Cy., 14, 669–681,, 2000. 

Small, C. and Nicholls, R. J.: A global analysis of human settlement in coastal zones, J. Coastal Res., 19, 584–599, available at: (last access: August 2019), 2003. 

Smith, C. L. and Tett, P.: A depth-resolving numerical model of physically forced microbiology at the European shelf edge, J. Marine Syst., 26, 1–36,, 2000. 

Soetaert, K., Hofmann, A. F., Middelburg, J. J., Meysman, F. J. R., and Greenwood, J.: The effect of biogeochemical processes on pH, Mar. Chem., 105, 30–51,, 2007. 

Tett, P.: A three-layer vertical and microbiological processes model for shelf seas, Report No. 14. Proudman Oceanographic Laboratory, Birkenhead, UK, 85 pp., 1990. 

Thingstad, T. F.: Utilization of N, P, and organic C by heterotrophic bacteria. I. Outline of a chemostat theory with a consistent concept of “maintenance” metabolism, Mar. Ecol. Prog. Ser., 35, 99–109,, 1987. 

Uppström, L. R.: The boron/chlorinity ratio of deep-sea water from the Pacific Ocean, Deep Sea Res. Oceanogr. Abstr., 21, 161–162,, 1974. 

Upstill-Goddard, R. C.: Air–sea gas exchange in the coastal zone, Estuar. Coast. Shelf S., 70, 388–404,, 2006. 

Wanninkhof, R. H.: Relationship between wind speed and gas exchange, J. Geophys. Res., 97, 7373–7382, 1992. 

Weiss, R. F.: Carbon dioxide in water and seawater: the solubility of a non-ideal gas, Mar. Chem., 2, 203–215,, 1974. 

Willmott, C. J.: Some comments on the evaluation of model performance, B. Am. Meteorol. Soc., 63, 1309–1313,<1309:SCOTEO>2.0.CO;2, 1982. 

Wimart-Rousseau, C., Lajaunie-Salla, K., Marrec, P., Wagener, T., Raimbault, P., Lagadec, V., Lafont, M., Garcia, N., Diaz, F., Pinazo, C., Yohia, C., Garcia, F., Xueref-Remy, I., Blanc, P. E., Armengaud, A., and Lefèvre, D.: Temporal variability of the carbonate system and air–sea CO2 exchanges in a Mediterranean human-impacted coastal site, Estuar. Coast. Shelf S., 236,, 2020. 

Xueref-Remy, I., Milne, M., Zoghbi, N., Yohia, C., Armengaud, A., Blanc, P.-E., Delmotte, M., Piazzola, J., Nathan, B., Ramonet, M., and Lac, C.: Assessing atmospheric CO2 variability in the Aix–Marseille metropolis area (France) and its coastal Mediterranean Sea at different time scales within the AMC project, Prague, Aus, available at: (last access: February 2020), 2018a. 

Xueref-Remy, I., Dieudonné, E., Vuillemin, C., Lopez, M., Lac, C., Schmidt, M., Delmotte, M., Chevallier, F., Ravetta, F., Perrussel, O., Ciais, P., Bréon, F.-M., Broquet, G., Ramonet, M., Spain, T. G., and Ampe, C.: Diurnal, synoptic and seasonal variability of atmospheric CO2 in the Paris megacity area, Atmos. Chem. Phys., 18, 3335–3362,, 2018.  

Yohia, C.: Genèse du mistral par interaction barocline et advection du tourbillon potentiel, Climatologie, 13, 24–37,, 2017. 

Zappa, C. J., Raymond, P. A., Terray, E. A., and McGillis, W. R.: Variation in surface turbulence and the gas transfer velocity over a tidal cycle in a macro-tidal estuary, Estuaries, 26, 1401–1415, 2003. 

Short summary
A biogeochemical model of planktonic food webs including a carbonate balance module is applied in the Bay of Marseille (France) to represent the carbon marine cycle expected to change in the future owing to significant increases in anthropogenic emissions of CO2. The model correctly simulates the ranges and seasonal dynamics of most variables of the carbonate system (pH). This study shows that external physical forcings have an important impact on the carbonate equilibrium in this coastal area.