Comment on gmd-2023-34

Abstract

<strong class="journal-contentHeaderColor">Abstract.</strong> The Bay of Marseille (BoM), located in the north-western Mediterranean Sea, is affected by various hydrodynamic processes (e.g., Rh&ocirc;ne River intrusion and upwelling events) that result in a highly complex local carbonate system. In any complex environment, the use of models is advantageous since it allows to identify the different environmental forcings, thereby facilitating a better understanding. By combining approaches from two biogeochemical ocean models and improving the formulation of total alkalinity, we develop a more realistic representation of the carbonate system variables at high temporal resolution which enables us study air-sea CO₂ fluxes and seawater <em>p</em>CO₂ variations more reliably. We apply this new formulation to two particular scenarios, typical for the BoM: (i) summer upwelling and (ii) Rh&ocirc;ne River intrusion events. In both scenarios, our model was able to correctly reproduce the observed patterns of <em>p</em>CO₂ variability. Summer upwelling events are typically associated with <em>p</em>CO₂ decrease that mainly results from decreasing near-surface temperatures. Furthermore, Rh&ocirc;ne River intrusion events are typically associated with <em>p</em>CO₂ decrease, although in this case the <em>p</em>CO₂ decrease results from a decrease in salinity and an overall increase in total alkalinity. While our model was able to correctly represent the daily range of air-sea CO₂ fluxes, we were unable to correctly estimate the yearly total air-sea CO₂ flux. Although the model consistent with observations, predicted the BoM to be a sink of CO₂ on a yearly basis, the magnitude of this CO₂ sink was underestimated which may be an indication of the limitations inherent in dimensionless models for representing air-sea CO₂ fluxes.