Air-Sea CO2 Exchange in the Strait of Gibraltar

The seasonal and spatial variability of the CO₂ system and air-sea fluxes were studied in surface waters of the Strait of Gibraltar between February 2019 and March 2021. High-resolution data was collected by a surface ocean observation platform aboard a volunteer observing ship. The CO₂ system was strongly influenced by temperature and salinity fluctuations forced by the seasonal and spatial variability in the depth of the Atlantic–Mediterranean Interface layer and by the tidal and wind-induced upwelling. The changes in seawater CO₂ fugacity (fCO_2,sw) and fluxes were mainly driven by temperature despite the significant influence of non-thermal processes in the southernmost part. The thermal to non-thermal effect ratio (T/B) reached maximum values in the northern section (>1.8) and minimum values in the southern section (<1.30). The fCO_2,sw increased with temperature by 9.02 ± 1.99 μatm °C^–1 (r² = 0.86 and ρ = 0.93) and 4.51 ± 1.66 μatm °C^–1 (r² = 0.48 and ρ = 0.69) in the northern and southern sections, respectively. The annual cycle of total inorganic carbon normalized to a constant salinity of 36.7 (NC_T) was assessed. Net community production processes described 93.5–95.6% of the total NC_T change, while air-sea exchange and horizontal and vertical advection accounted for <4.6%. The fCO_2,sw in the Strait of Gibraltar since 1999 has been fitted to an equation with an interannual trend of 2.35 ± 0.06 μatm year^–1 and a standard error of estimate of ±12.8 μatm. The seasonality of the air-sea CO₂ fluxes reported the behavior as a strong CO₂ sink during the cold months and as a weak CO₂ source during the warm months. Both the northern and the southern sections acted as a net CO₂ sink of −0.82 and −1.01 mol C m^–2 year^–1, respectively. The calculated average CO₂ flux for the entire area was −7.12 Gg CO₂ year^–1 (−1.94 Gg C year^–1).

Introduction

Coastal and transition to open-ocean areas are biogeochemically active regions which play a key role in the global carbon cycle by assembling multiple diverse ecosystems and linking the terrestrial, oceanic, and atmospheric carbon reservoirs. These regions are characterized by primary production, carbon fixation, and remineralization ratios significantly higher than in the open oceans due to higher supplies of organic and inorganic carbon (Walsh, 1991; Mackenzie et al., 1998; Borges et al., 2005; Muller-Karger et al., 2005) and show a high spatio-temporal variability in terms of air-sea CO₂ exchange, often acting as a CO₂ sink and source during cold and warm months, respectively (e.g., Frankignoulle and Borges, 2001; Cai et al., 2006; Shadwick et al., 2010, 2011). However, the air-sea CO₂ fluxes have been difficult to quantify in these areas due to the lack of underway field data and the subsequent limited number of local and regional studies. Hence, continuous high-frequency monitoring is required to assess the ocean CO₂ uptake with a higher certainty. It is important to remark and understand the role of regional areas in the climate change. In fact, the recent IPCC report 2021 (IPCC, 2021) together with the 2007 report (IPCC, 2007), highlighted the need of study of marginal seas, continental shelves, and semi-enclosed seas due to the variability on the biogeochemical cycles as well as the human pressure.

In this context, the Strait of Gibraltar (Figure 1) is a key region in physical and biogeochemical terms influenced by several different processes and therefore characterized by its high variability. The physical processes that govern the hydrodynamics in the Strait of Gibraltar (e.g., Lacombe and Richez, 1982; Gascard and Richez, 1985; Armi and Farmer, 1988; Bryden and Kinder, 1991) are involved in the variation of biogeochemical properties and have been analyzed in terms of physical–biological coupling in previous studies (Echevarría et al., 2002; Macias et al., 2006; Macías et al., 2007; Ramírez-Romero et al., 2014). However, the knowledge about the surface CO₂ system and air-sea fluxes is scarce due to insufficient data.

FIGURE 1

Figure 1. (A) Bathymetric map of the Strait of Gibraltar and location of the topographic features of interest with west-to-east distribution: Majuan Bank (MB), Espartel Sill (ES), Camarinal Sill (CS), Tarifa Narrow (TN), and Gibraltar Central (GC). (B) CanOA-VOS northern (blue) and southern (red) tracks and location of the stations of interest from S1 to S5.

Previous studies have explored the distribution of the CO₂ system parameters in the water column (Dafner et al., 2001; Santana-Casiano et al., 2002), reported an increase in total inorganic carbon from surface to the bottom of 2068–2150 and 2307–2325 μmol kg^–1 in surface Atlantic and Mediterranean waters, respectively. The annual cycle and spatial distribution of the surface CO₂ system variables was analyzed by de la Paz et al. (2009), which accounted a seasonal change in CO₂ fugacity of ∼60 μatm mainly controlled by temperature fluctuations. The Atlantic–Mediterranean interaction was also assessed in terms of natural and anthropogenic carbon exchange (Huertas et al., 2009) and decadal acidification (Flecha et al., 2019). The air-sea CO₂ exchange evaluations have determined that the Strait of Gibraltar behaves as a net atmospheric CO₂ sink on an annual scale (Santana-Casiano et al., 2002; de la Paz et al., 2009, 2011). Nevertheless, these studies are based on temporarily limited observations at certain locations along the channel and new strategies for continuous monitoring of oceanic CO₂ are required to develop high spatio-temporal resolution studies of the CO₂ system in highly variables areas such as the Strait of Gibraltar.

The aim of this study was to analyze the spatio-temporal variability of the surface CO₂ system and air-sea exchange along the Strait of Gibraltar. The processes that control the distribution of the CO₂ system parameters and fluxes were attended. The latitudinal variability related to the influence of different biogeochemical processes near both the Iberian and African coasts, the interaction between the Northeast Atlantic and the Mediterranean Sea, the effect of west-to-east surface circulation and other physical processes acting on the Strait were analyzed through the northernmost and southernmost part of the channel. These data will improve our knowledge about the role of coastal regions, marginal and semi-enclosed seas and coastal to open-ocean transitional areas in terms of air-sea CO₂ exchange.

Study Area

The Strait of Gibraltar (Figure 1) is a narrow (∼20 km) and shallow (∼600 m) channel with a west-to-east orientation that connects the Mediterranean Sea with the Northeast Atlantic. The minimum width at the Tarifa Narrow (TN, ∼14 km) and minimum depth at the Camarinal Sill (CS, ∼290 m) and at the Espartel Sill (ES, ∼360 m) play an important role in the exchange of water through the channel (e.g., Sánchez-Garrido et al., 2008, 2011; García-Lafuente et al., 2011; Sammartino et al., 2015).

The circulation pattern in the Strait has been described as a two-layer system: a surface Atlantic water inflow toward the east and a deep Mediterranean water outflow toward the west (e.g., Lacombe and Richez, 1982; Gascard and Richez, 1985; Armi and Farmer, 1988; Bryden and Kinder, 1991). The fresher and nutrient-depleted Atlantic water is separated from the saltier and nutrient-rich Mediterranean water by the Atlantic–Mediterranean Interface (AMI) layer (Lacombe and Richez, 1982; Bray et al., 1995). The AMI is deeper, thicker and colder in the westernmost part of the channel and slopes up toward the northeast from ∼200 m depth in ES to ∼75 m depth eastward of Gibraltar Central (GC) (Bray et al., 1995; Huertas et al., 2009). The shallow position of the AMI promotes the deep-water upwelling in the eastern and north-eastern region of the Strait (Minas et al., 1991; Echevarría et al., 2002; Gómez-Jakobsen et al., 2019).

The main physical components of the circulation and flow through the channel were detailed by Candela et al. (1990): a barotropic tidal flow, a barotropic sub-inertial component driven by atmospheric pressure fluctuations in the Mediterranean Sea, a long-term baroclinic component term driven by differences in density between Mediterranean Sea and Northeast Atlantic and short-period currents associated with large-amplitude internal waves. The formation of internal waves mainly occurs in the shallower area around CS and is induced by the interaction of the tidal flow with the change in bathymetry and vertical stratification through the channel (e.g., Armi and Farmer, 1988; La Violette and Arnone, 1988; Richez, 1994; Bruno et al., 2002; Alonso del Rosario et al., 2003; Vázquez-Escobar et al., 2008). Meteorological sub-inertial forcing also influences the release of these internal waves (Candela et al., 1989; Pistek and La Violette, 1999). They are propagated toward the Mediterranean on a diurnal time scale controlled by the tides and at highly variable speeds (Watson and Robinson, 1990; Richez, 1994; Sánchez-Garrido et al., 2008), while being trapped on the lee side of CS (Bruno et al., 2002). The internal waves are an important contributor to the mixing between the Mediterranean and Atlantic layers by intermittently injecting deep and nutrient-rich water into the upper layers (Wesson and Gregg, 1994; Macias et al., 2006).

The unique hydrology of the channel influenced by its circulation pattern and the generation of physical features, as well as the local climatology, are involved in the variability of the physical and biogeochemical properties in the Strait of Gibraltar. The tidal-induced upwelling phenomenon is combined with wind-induced upwelling events along the African coast during easterly winds (Stanichny et al., 2005) and both represent a source of high spatio-temporal variability of the surface CO₂.

Materials and Methods

The CO₂ system and air-sea exchange was assessed in the Strait of Gibraltar based on high spatio-temporal frequency data obtained from February 2019 to March 2021 through continuous autonomous monitoring carried out by a surface ocean observation platform (SOOP) installed on a volunteer observing ship (VOS), SOOP CanOA-VOS. The SOOP CanOA-VOS was the container ship RENATE P (IMO: 9144718) operated by Nisa Marítima, which usual route runs between the Canary Islands and Barcelona, through the eastern Canary archipelago waters and northward along the northwest African and east-southeast coast of the Iberian Peninsula. Data collection and maintenance of autonomous measurement systems took place biweekly (time required to perform a round trip) at the port of Las Palmas de Gran Canaria (28.1319°N, 15.4185°W). In total, 52 routes crossed the Strait of Gibraltar: 34 upward routes through the southernmost part (southern section) and 18 downward routes through the northernmost part (northern section). The SOOP CanOA-VOS line will be part of the Spanish contribution to the Integrated Carbon Observation System (ICOS) international program in 2021. Thus, the measurement equipment and data collection process verify their quality requirements and methodological recommendations to produces comparable and high-quality dataset.

Underway Measurements and Data Collection

Autonomous monitoring was carried out by an automated underway seawater and low atmospheric CO₂ molar fraction (xCO₂) measurement system installed inside the engine room of the RENATE P container ship (Curbelo-Hernández et al., 2021). The xCO₂ measurement system, developed by Craig Neill and commercialized by General Oceanics^TM, combines an air and seawater equilibrator with a non-dispersive infrared analyzer for gas detection (Pierrot et al., 2009).

The system was checked automatically on departure and arrival at each port, and periodically every 3 h during the ship’s journey. Four standard gasses (in the order of 0, 250, 400, and 540 ppm, with a ±0.02 ppm accuracy) provided by the National Ocean and Atmospheric Administration (NOAA) and traceable to the World Meteorological Organization (WMO) were used for system check and xCO₂ measurement corrections. Periodic calibrations were automatically performed every 12 h using standard gasses 1 and 4 to adjust the, respectively, zero and span of the infrared analyzer.

The system measures xCO₂ (ppm) in seawater with a frequency of 2–3 min and normally collect 130–150 data in the Strait of Gibraltar. The system also performs three xCO₂ measurements in low atmosphere after each calibration. In both cases, a non-dispersive infrared analyzer built by LICOR^® (initially the 6262 model and after October 2019, a 7000 model) were used for the detection of xCO₂.

The sea surface temperature (SST, in °C) was measured using an SBE38 thermometer with an estimated error of 0.01°C placed at the main seawater intake. In addition, due to the high sensitivity of xCO₂ to temperature changes, an SBE45 thermosalinograph and a Hart Scientific HT1523 Handheld Thermometer, with an estimated error of ±0.01°C each, were used to monitor the temperature near the xCO₂ system and inside the equilibrator, respectively.

Sea surface salinity (SSS) was measured with a manually calibrated SBE45 thermosalinograph, with an estimated error of ±0.005. Atmospheric pressure was measured by the transducer on the deck (these pressure records being different in the order of millibars with the pressure near the xCO₂ system, inside the ship).

Due to technical problems with the measurement equipment aboard the SOOP CanOA-VOS, no measurements were obtained during March 2019 (problems with water intake), during the second half of August 2019 and March–April 2020 (problems with LICOR flow that could not be resolved due to COVID-19 limitations) and on a limited number of return journeys from the port of Barcelona.

Determination of Variables

This study followed the data collection methodology, quality control, and calculation procedures as published in the updated version of the DOE method manual for ocean CO₂ analysis (Dickson et al., 2007). The fugacity of CO₂ was calculated in seawater (fCO_2,sw) and in the lower atmosphere (fCO_2,atm) from measured and corrected xCO₂ values following Pierrot et al. (2009).

The thermal and non-thermal contribution on fCO_2,sw and variations was studied using the equations presented by Takahashi et al. (2002), with the temperature effects on fCO_2,sw for isochemical seawater of 0.0423°C^–1 determined experimentally by Takahashi et al. (1993). The non-thermal effect (fCO_{2,non–thermal}) was obtained from the seasonal amplitude of the normalized fCO_2,sw to the average temperature (Eqs 1, 2). The thermal effect (fCO_2,thermal) was calculated by applying the observed temperature effect to the average value of fCO_2,sw (Eq. 3) and determining its seasonal amplitude (Eq. 4). The relative importance of thermal and non-thermal effects was expressed by the T/B ratio (ΔfCO_2,thermal/ΔfCO_{2,non–thermal}), with values greater than 1 indicating that the temperature effect dominated over non-physics effects.

$f{\text{CO}}_2\left({\text{T}}_{\text{mean} }\right)={\left(f{\text{CO}}_2\right)}_{\text{obs} }\cdot \exp [0.0423\left({\text{T}}_{\text{mean} }-{\text{T}}_{\text{obs} }\right)](1)$

${{\left(\mathrm{\Delta }f{\text{CO}}_2\right)}_{\text{bio} }=\text{[}f{\text{CO}}_2\left({\text{T}}_{\text{mean} }\right)]}_{\text{max} }-\text{[}f{\text{CO}}_2\left({\text{T}}_{\text{mean} }\right)]{}_{\text{min} }(2)$

$f{\text{CO}}_2\left({\text{T}}_{\text{obs} }\right)={\left(f{\text{CO}}_2\right)}_{\text{mean} }\cdot \exp [0.0423\left({\text{T}}_{\text{obs} }-{\text{T}}_{\text{mean} }\right)](3)$

${{\left(\mathrm{\Delta }f{\text{CO}}_2\right)}_{\text{temp} }=\text{[}f{\text{CO}}_2\left({\text{T}}_{\text{obs} }\right)]}_{\text{max} }-\text{[}f{\text{CO}}_2\left({\text{T}}_{\text{obs} }\right)]{}_{\text{min} }(4)$

The CO₂ fluxes (FCO₂) were determined using Eq. 5 (Broecker and Peng, 1982) with a conversion factor of 0.24 mmol m^–2 d^–1. The solubility (S, mol L^–1 atm^–1) and the difference between seawater and low atmosphere fCO₂ (ΔfCO₂ = fCO_2,sw − fCO_2,atm) were considered. Negative fluxes indicate that the ocean acts as an atmospheric CO₂ sink, while the positive ones indicate that it behaves as a source.

${\text{FCO}}_2=0.24\cdot \text{S}\cdot k\cdot \mathrm{\Delta }f{\text{CO}}_2(5)$

The updated Wanninkhof (2014) parameterization was used in this study due to it considers the most recently advances in the quantification of the input parameters and improves the wind speed products, with k being the gas transfer rate expressed in Eq. 6:

$k=0.251\cdot {\text{w}}^2\cdot {\left(\frac{\text{Sc}}{660}\right)}^{-0.5}(6)$

where w is the wind speed (m s^–1) and Sc is Schmidt number (cinematic viscosity of seawater, divided by the gas diffusion coefficient). Daily averages of satellite ocean surface wind speed derived from the Advanced Scatterometer (ASCAT) and obtained from the Satellite Research and Exploitation Center (CERSAT) at IFREMER (Plouzane, France) were considered together with daily measured wind speed averages provided by the Agencia Estatal de Meteorología (AEMET) from the meteorological station located at Tarifa (Cádiz) and corrected to 10 m height (Allen et al., 1998). These wind speed data were interpolated based on the latitude, longitude, and time of the underway xCO₂ measurements. The average ocean surface wind speed for the entire study period calculated from satellite and measured data was 7.77 ± 2.66 and 8.18 ± 2.61 m s^–1, respectively. Due to these wind speed values were consistent, the FCO₂ was calculated in this study by considering the wind speed data measured at the Tarifa meteorological station as was done previously by de la Paz et al. (2011) to assess the air-sea CO₂ exchange in the Strait of Gibraltar.

Autonomous monitoring was combined with 25 discrete surface seawater samples taken manually from the intake seawater line for total alkalinity (A_T, μmol kg^–1) and total dissolved inorganic carbon (C_T, μmol kg^–1) in February 2020 and March 2021 with in situ measurements of SST and SSS. Samples were taken with a temporal frequency of 1–2 h in borosilicate glass bottles, overfilled, preserved with 100 μl of saturated HgCl₂, kept in darkness and analyzed just after arriving at port, in a period less than 2 weeks. The A_T was determined by potentiometrically titration in a VINDTA 3C system following Mintrop et al. (2000) and using CRMs (provided by A. Dickson at Scripps Institution of Oceanography) to test the performance of the titration system and correct the determined values. The CRMs titration gave values with and accuracy of ±1.5 μmol kg^–1.

The determined A_T values were consistent with those estimated according to the global relationship of A_T with SSS and SST in global surface waters described by Lee et al. (2006) for the North Atlantic (presented an average difference of 3.08 ± 6.80 μmol kg^–1). Due to the similarity of both values, A_T was calculated for the time, longitude, and latitude of the CanOA-VOS data using the A_T-SSS relationship obtained from the discrete samples (Eq. 7; r² = 0.87 and Pearson correlation coefficient given hereinafter as ρ is equal to 0.93) (temperature was not found to improve the fitting of the experimental data). It was assumed that the change in A_T with SSS was constant throughout the year in tropical and subtropical latitudes (Lee et al., 2006).

$A_T=2501-241.9\cdot \left(SSS-35\right)-112.2\cdot {\left(SSS-35\right)}^2(7)$

The C_T and pH in total scale (pH_T) were computed with the Excel program CO₂sys, using the carbonic acid dissociation constants of Lueker et al. (2000), the HSO₄ dissociation constant of Dickson (1990) and the value of [B]_I determined by Lee et al. (2010). The longitudinal change in the seasonality of the variables measured in situ (SST and SSS), determined (fCO_2,sw and FCO₂) and computed (C_T and pH_T) was studied in the entire region along the northern and southern sections. The C_T was normalized (NC_T) to a constant SSS of 36.4 (the annual average SSS) to remove the effect of evaporation and precipitation. Similarly, the pH_T values were normalized (pH_T,18) to a constant SST of 18.0°C (annual average SST) to remove the temperature dependence on the variation of pH_T (Santana-Casiano et al., 2001; González-Dávila et al., 2003).

The temporal variation of C_T was studied using the equation provided by Sarmiento and Gruber (2006) (Eq. 8). This equation describes the change in C_T over time produced by air-sea CO₂ exchange (EX), horizontal advection (TRSP), vertical mixing (MIX), and net community production (NCP).

$\frac{{\text{dNC}}_{\text{T} }}{\text{dt}}={\frac{{\text{dNC}}_{\text{T} }}{\text{dt}}|}_{\text{EX} }+{\frac{{\text{dNC}}_{\text{T} }}{\text{dt}}|}_{\text{TRSP} }+{\frac{{\text{dNC}}_{\text{T} }}{\text{dt}}|}_{\text{MIX} }+{\frac{{\text{dNC}}_{\text{T} }}{\text{dt}}|}_{\text{NCP} }(8)$

The increase and depletion of NC_T in surface waters was studied during a whole year (February–January) using monthly average NC_T values between February 2019 and February 2021 and considering the variation of the mixed layer depth (MLD). The temporal change in NC_T was calculated as the difference between two consecutive months and was expressed in mmol m^–3, considering the seawater density. Negative values indicate that NC_T decreases in surface waters. The relative contribution of CO₂ exchange in the variation of NC_T (EX term) is given by the relationship between FCO₂ and MLD. The monthly MLD values were obtained with a spatial resolution of 0.028° × 0.028° from the operational Iberian–Biscay–Irish (IBI) Ocean Analysis and Forecasting system based on a (eddy-resolving) NEMO model application and are available at E.U. Copernicus Marine Service Information.¹ The annual MLD cycle was calculated using the harmonic fitting equation and parameters shown in Supplementary Table 1.

The relative contribution of horizontal transport (TRSP term) was calculated by considering a surface transport from the Northeast Atlantic to the Mediterranean Sea. The spatial change in NC_T between westernmost and easternmost waters [ΔNC_T = NC_{T, [=est)} − NC_{T, (East)}] was considered. Data available at the Surface Ocean CO₂ Atlas (SOCAT v2020²) referenced to 2019 have been used for the calculation of NC_T at a reference station to the west of the Strait of Gibraltar (6.75 ± 0.25°W) in a similar latitudinal range. The spatial change in salinity was removed by calculating the NC_T at all the selected locations in the Strait of Gibraltar at the constant salinity of the reference station. The horizontal transport contribution to NC_T fluctuations was obtained from the product of ΔNC_T and the change in salinity with time at the locations of interest in the Strait of Gibraltar. Positive values indicate that horizontal transport occurred in a west-to-east direction, while negative values indicated that the transport direction was reversed.

The relative contribution of vertical mixing (MIX term) was estimated following Bégovic and Copin-Montégut (2002) by considering the sum of changes in vertical diffusion across the bottom of the mixed layer and vertical entrainment into the surface layer (Eq. 9). This term was determined using the monthly values of total dissolved inorganic carbon in the MLD obtained from IBI MFC high-resolution biogeochemical forecasts and generated from the PISCES biogeochemical model and the NEMO ocean circulation model. Additional monthly values of temperature and salinity obtained from the operational IBI Ocean Analysis and Forecasting system based on a NEMO model were used to estimate the seawater density. These forecast data have a spatial resolution of 0.028° × 0.028° and are available in E.U. Copernicus Marine Service Information.¹

${\mathrm{△}N{C}_T|}_{MIX}=\frac{1}{MLD}{K}_Z\frac{\partial N{C}_T}{\partial z}+\frac{1}{h}\mathrm{\theta }\left(\frac{\partial MLD}{\partial t}\right)(N{C}_{T,f}-N{C}_T)(9)$

The difference between NC_T concentration below (NC_T,f) and at the base (NC_T) of the mixed layer was considered to calculate the term of the vertical entrainment. The function θ is equal to 0 when the mixed layer decreases (∂⁡MLD/∂⁡t) ≤ 0 due to the fact that only the deepening of the mixed layer induces mixing with the underlying waters (Peng et al., 1987; Fasham et al., 1990; Gruber et al., 1998). The vertical gradient of NC_T at the base of the mixed layer (∂⁡NC_T/∂⁡z) was considered for the calculation of vertical diffusion. The vertical diffusion coefficient (K_z, m^–2 d^–1) at the base of the mixed layer was determined according to Denman and Gargett (1983) using Eq. 10:

$K_z=\frac{0.25\mathrm{\epsilon }\mathrm{\rho }}{g\left(\frac{\partial \mathrm{\rho }}{\partial z}\right)}(10)$

where ρ is the seawater density at the base of the mixed layer, g is the acceleration of gravity and ε is the rate of turbulent energy dissipation. In this study we considered the vertical density gradient (∂⁡ρ/∂⁡z) in the pycnocline and a representative value of 2.0 × 10^–8 m² s^–3 for ε (Gruber et al., 1998).

Lastly, the change in NC_T due to net community production processes (NCP term), considering total changes in photosynthesis/respiration processes, was obtained directly from Eq. 8. Positive contributions of NCP indicated that the remineralization processes increase the NC_T in surface waters while negative contributions shows that biological production decrease it surface concentration.

Computational Methods

The raw output data was initially filtered removing data affected by the automatic sampler such as samples measured at low water rates (<2.5 L min^–1). The measured xCO₂ values were corrected following Pierrot et al. (2009) by linearly interpolating the certified standard values (section 3.2) with the xCO₂ measurement times. The three fCO_2,atm values obtained after each calibration were averaged and harmonically interpolated with the times of each fCO_2,sw. The ΔfCO₂ was computed from the difference between the fCO_2,sw and the fCO_2,atm values and was used in the FCO₂ calculations (section 2.2). Additional daily average satellite data for Chlorophyll a obtained from the Operational Mercator Ocean biogeochemical global ocean analysis and forecast system (United States Copernicus Marine Service Information) were used to improve the understanding of the biological influence on the CO₂ system variability. Chlorophyll a data was interpolated with latitude, longitude, and time of the CanOA-VOS line data.

The average physical and biogeochemical variables (y) for each selected station were fitted to a harmonic Eq. 11 as a function of time (x) in order to study the seasonal variability and the thermal and non-thermal effect on fCO_2,sw (Takahashi et al., 2002; Lüger et al., 2004). The a–e coefficients are shown in Supplementary Table 2.

$\mathrm{y}=\mathrm{a}+\mathrm{b}\text{⋅}\cos \left(2\mathrm{\pi }\text{x}\right)+\mathrm{c}\text{⋅sin}\left(2\mathrm{\pi }\text{x}\right)+\mathrm{d}\text{⋅cos}\left(4\mathrm{\pi }\text{x}\right)+\mathrm{e}\text{⋅sin}\left(4\mathrm{\pi }\text{x}\right)(11)$

The results were compared to surface data collected by other VOS and research vessels in the oceanographic environment of the Strait of Gibraltar and are available at the SOCAT v2020 database. The SOCAT data between 1999 and 2019 were used to study the monthly frequency and the annual cycle of fCO_2,sw in 2019 by considering an interannual rate of increase of 1.8 μatm per elapsed year (Bates et al., 2014).

Results

The present study provides a high-resolution evaluation of the spatio-temporal variability of surface CO₂ in the Strait of Gibraltar. The annual and seasonal averages of the carbon variables in both the northern and the southern routes are presented in Table 1 with maximum values of SST, fCO_2,sw, and FCO₂ and minimum values of C_T and pH_T during summer. The greatest seasonal change occurs from winter to summer: the SST, fCO_2,sw, and FCO₂ increased in the northern section by 4.79 ± 1.97°C, 40.75 ± 28.83 μatm, and 5.52 ± 4.19 mmol m^–2 d^–1, respectively, leading to a pH_T decrease of 0.039 ± 0.028 U. The seasonal change along the southern section was slightly lower on average except for FCO₂, C_T, and NC_T (5.67 ± 5.84 mmol m^–2 d^–1, 24.72 ± 19.1 μmol m^–3, and 23.20 ± 18.15 μmol m^–3, respectively). The seasonal change in pH along the southern section was 0.026 ± 0.014 U. In contrast, spring and autumn averages were similar for all the considered variables in both sections.

TABLE 1

Table 1. Annual and seasonal average of SST, SSS, and CO₂ system variables and fluxes with their respective standard deviations along the northern and southern sections of the Strait of Gibraltar.

Spatial Distribution of Sea Surface Temperature, Sea Surface Salinity, and CO₂ System Parameters and Fluxes

The longitudinal west-to-east distribution of the variables of interest was studied (Supplementary Figures 1, 2) and showed differences, especially during summer. The longitudinal variation was analyzed separately in the northern and southern sections: the SST, SSS, and CO₂ system variables and fluxes were annually and seasonally averaged every 0.05° longitude (Figure 2) by considering the following periods: January–March (winter), April–June (spring), July–September (summer), and October–December (autumn). The average values show that SST drives the fCO_2,sw and FCO₂ variations. The greatest latitudinal differences were found toward the east of GC. The average SST, fCO_2,sw, FCO₂, and pH_T values follow an opposite west-to-east trend between the northern and southern sections. The seasonality of SST and fCO_2,sw was maximum (5.29 ± 2.97°C and 45.51 ± 28.11 μatm, respectively) in the westernmost part (6.15–5.95°W) and minimum in the eastern area closest to the Mediterranean Sea (5.30–5.10°W) in both the northern (4.47 ± 2.68°C and 19.11 ± 24.86 μatm, respectively) and the southern (3.39 ± 2.11°C and 23.81 ± 14.96 μatm, respectively) sections. The seasonal change in pH_T (northern and southern sections, respectively) was 0.053 ± 0.028 and 0.033 ± 0.024 U in the westernmost part and 0.019 ± 0.025 and 0.025 ± 0.015 U in the easternmost part.

FIGURE 2

Figure 2. Annual and seasonal average of (A) SST, (B) SSS, (C) fCO_2,sw, (D) C_T, (E) FCO₂, (F) NC_T, (G) pH_T, and (H) pH_T,18 every 0.05° of longitude along the (1) northern and (2) southern sections in the Strait of Gibraltar.

An annual average west-to-east decrease through the northern and southern routes was obtained for SST [1.70°C (r² = 0.64 and ρ = −0.80) and 1.48°C (r² = 0.83 and ρ = −0.91) per degree of longitude, respectively] and fCO_2,sw [14.21 μatm (r² = 0.65 and ρ = −0.81) and 2.71 μatm (r² = 0.22 and ρ = −0.47) per degree of longitude, respectively] (Figures 2A,E). The seasonality of pH_T increased toward the east in both sections (Figure 2G) and correlated with the decrease in fCO_2,sw (−0.001 U μatm^–1, r² > 0.94, and ρ < −0.97). These longitudinal changes were maximum during summer in the northern section: SST and fCO_2,sw decreased by 2.95°C (r² = 0.66 and ρ = −0.81) and 47.47 μatm (r² = 0.87 and ρ = −0.93) per degree of longitude, respectively, while the pH_T increased by 0.047 U (r² = 0.86 and ρ = 0.93) per degree of longitude.

The annual average FCO₂ (Figure 2E) shows that the CO₂ sink increased toward the east in both the northern (−1.89 mmol m^–2 d^–1 per degree of longitude, r² = 0.68, and ρ = 0.82) and the southern (−1.27 mmol m^–2 d^–1 per degree of longitude, r² = 0.57, and ρ = 0.75) sections. The west-to-east increasing sink gradient was most intense during spring [−2.23 (r² = 0.34 and ρ = 0.58) and −1.38 (r² = 0.64 and ρ = 0.80) mmol m^–2 d^–1 per degree of longitude, respectively] and summer [−6.71 (r² = 0.89 and ρ = 0.94) and −2.69 (r² = 0.60 and ρ = 0.77) mmol m^–2 d^–1 per degree of longitude, respectively]. During summer, both the northern and southern sections behaved as a CO₂ source (3.11 ± 1.76 and 1.24 ± 3.84 mmol m^–2 d^–1, respectively) in the westernmost part (6.15–5.5°W) and as a weak CO₂ sink (−0.87 ± 2.62 and −0.57 ± 2.48 mmol m^–2 d^–1) in the easternmost part (5.5–5.1°W). Maximum outgassing was found in both the northern and southern sections during summer around ES (4.49 ± 2.45 and 2.28 ± 3.91 mmol m^–2 d^–1, respectively) and TN (2.99 ± 0.74 and 2.20 ± 2.59 mmol m^–2 d^–1, respectively). In contrast, maximum ingassing was obtained along the northern section between ES and CS during autumn (−6.01 ± 0.05 mmol m^–2 d^–1) and toward the east of TN along the southern section during winter (−6.57 ± 4.11 mmol m^–2 d^–1).

The longitudinal change in SSS (Figure 2B) was significant in this transition area between Northeast Atlantic and Mediterranean Sea waters and influenced the spatial distribution of C_T. Big differences in SSS and C_T values were observed between the northern and southern sections. In the northern section, the annual average of SSS slightly decreased from the most Atlantic part to GC (0.12 U per longitude, r² = 0.66, and ρ = −0.81) and increased rapidly toward the east of GC (1.40 U per longitude, r² = 0.95, and ρ = 0.97). These longitudinal changes were also observed in the southern section, but both the decrease (0.08 U per degree of longitude, r² = 0.67, and ρ = −0.82) and the increase (0.50 U per degree of longitude, r² = 0.86, and ρ = 0.93) were lower.

The longitudinal distribution of C_T (Figure 2D) was correlated with that of SSS in the northern section (r² = 0.91 and ρ = 0.95), while this control was not observed in the southern section (r² = 0.13 and ρ = 0.36). The average annual NC_T (Figure 2F) increased similarly toward the east in both the northern (11.42 μmol kg^–1 per degree of longitude, r² = 0.64, and ρ = 0.80) and the southern sections (12.64 μmol kg^–1 per degree of longitude, r² = 0.80, and ρ = 0.89). The average NC_T values were maximum during the cold months (January–June) and minimum during the warm months (July–December) in the entire region, except in the northern section during autumn to the east of CS and during summer around TN. The decrease in fCO_2,sw between CS and GC in summer and the effects of remineralization and vertical mixing processes accounted for these changes. The absolute minimum values of NC_T in summer along each route should be highlighted: the minimum value eastward of GC in the northern section (2090.60 ± 4.99 μmol kg^–1) was followed by a reduction of the oceanic CO₂ sink and high values of fCO_2,sw, while that of the southern route in TN (2090.03 ± 2.95 μmol kg^–1) coincided with maximum values of SST, fCO_2,sw, and FCO₂ and may be related to an increase in biological production.

Seasonality of the CO₂ System and Fluxes

The seasonality of SST, SSS, and CO₂ system variables and fluxes were studied on both the northern and the southern routes (Supplementary Figures 3, 4) by averaging the values obtained on each trip and fitting them to a harmonic Eq. 11. The average values of fCO_2,sw and pH ranged between ∼360 and 430 μatm and between ∼8.03 and 8.09, respectively, in phase with a seasonal temperature change of ∼7.5°C (between ∼15.5 and 23.0°C). Several differences were observed between the northern and the southern sections: the fCO_2,sw and pH trends (Supplementary Figure 3) show a lower seasonality along the southern route despite the fact that the seasonal change in SST is maximum in this section. The fCO_2,sw values were compared with surface fCO_2,sw data available at the SOCAT v2020 database for the Strait of Gibraltar between 1999 and 2019 (Supplementary Figure 5). The CanOA-VOS lines provided a total of 6440 measurements during the 2 years of study (2019–2021), almost the same number of observations (6063 data) as in the SOCAT database for a period of 20 years (1999–2019). The CanOA-VOS line collected data at times of the year when there was no previous data available (February, May–July, and December). The annual cycle of fCO_2,sw provided from both the CanOA-VOS and the SOCAT databases is described by Eqs 12, 13 (r² = 0.67 and r² = 0.64, respectively) in which SST fluctuations were considered (SSS was not found to improve the fitting). Their standard errors of estimate was ±10.2 and ±13.8 μatm, respectively.

$fC{O}_{2,sw}^{CanOA-VOS}=1118.00+0.16{}^{*}\cos \left(0.53{}^{*}Month\right)-15.14{}^{*}\sin \left(0.53{}^{*}Month\right)+2.10{}^{*}SS{T}^2-78.81{}^{*}SST(12)$

$fC{O}_{2,sw}^{SOCAT}=760.60+2.29{}^{*}\cos \left(0.55{}^{*}Month\right)-12.74{}^{*}\sin \left(0.55{}^{*}Month\right)+1.18{}^{*}SS{T}^2-42.54{}^{*}SST(13)$

The seasonal change of fCO_2,sw with SST obtained through the northern and southern sections of the CanOA-VOS line is given by Eqs 14, 15 (r² = 0.78 and r² = 0.66, respectively) with a standard error of estimate of ±10.3 and ±8.9 μatm, respectively.

$fC{O}_{2,sw}^{CanOA-VOS-North}=914.40+4.38{}^{*}\cos \left(0.54{}^{*}Month\right)-11.36{}^{*}\sin \left(0.54{}^{*}Month\right)+1.73{}^{*}SS{T}^2-60.88{}^{*}SST(14)$

$fC{O}_{2,sw}^{CanOA-VOS-South}=1118.00-2.00{}^{*}\cos \left(0.53{}^{*}Month\right)-14.84{}^{*}\sin \left(0.53{}^{*}Month\right)+2.08{}^{*}SS{T}^2-78.82{}^{*}SST(15)$

When the full data set (SOCAT and CanOA) was considered (more than 12,000 data), Eq. 16 (r² = 0.62), which present a standard error of estimate of ±12.8 μatm, provides an estimation of the interannual trend (2.33 ± 0.06 μatm year^–1), the seasonal variability and the influence of SST as the variable controlling the physical characteristics of the Strait of Gibraltar.

$fC{O}_{2,sw}=779.6+2.35{}^{*}(year-1999)+3.80{}^{*}\cos \left(0.56{}^{*}Month\right)-11.70{}^{*}\text{sin}(0.56{}^{*}Month)+1.34{}^{*}SS{T}^2-48.89{}^{*}SST(16)$

The FCO₂ seasonal trend (Supplementary Figure 3) ranged between −6 and 3 mmol m^–2 d^–1 in the northern section and between −8 and 3 mmol m^–2 d^–1 in the southern section, with a strong behavior as a sink during the cold months and as a weak source during the warm months. The maximum outgassing (2.99 ± 2.70 and 3.60 ± 1.76 mmol m^–2 d^–1 in the northern and the southern section, respectively) occurred in August 2020 and coincided with high values of SST and fCO_2,sw, while the maximum ingassing (−11.48 ± 1.36 mmol m^–2 d^–1) was observed along the southern section in February 2021 and was related with wind speed values higher than 11 m s^–1. The change in FCO₂ was linked to the change in ΔfCO₂, which is controlled by temperature fluctuations: the average values of FCO₂ increased by 1.22 mmol m^–2 d^–1 °C^–1 (r² = 0.71 and ρ = 0.84) in the northern section and 1.03 mmol m^–2 d^–1 °C^–1 (r² = 0.50 and ρ = 0.71) in the southern section. A smoothing spline with a parameter of 0.98 was applied to the daily distribution of SSS (r² > 0.27), SST (r² > 0.82), fCO_2,atm (r² > 0.94), fCO_2,sw (r² > 0.56), and wind speed (r² > 0.49). Daily values of FCO₂ were calculated from these computed variables and fitted to the observed FCO₂ values to improve the evaluation of its temporal distribution. These computed FCO₂ values differ in average by −0.25 mmol m^–2 d^–1 in the northern section and by 0.15 mmol m^–2 d^–1 in the southern section compared to those calculated from the observed data.

The west-to-east differences were assessed by selecting five stations equidistantly spaced every 0.2° of longitude along the channel (6.0–5.2°W) (Figure 3 and Supplementary Figures 6, 7). The westernmost station at 6.0 ± 0.05°W (S1) is located near ES and was used to study the seasonality in the most Atlantic part over the main gateway of the Mediterranean outflow. The station at 5.8 ± 0.05°W (S2) was selected near the shallower area of CS. The stations at 5.6 ± 0.05°W (S3) and 5.4 ± 0.05°W (S4) were in the longitudinal range around TN and GC, respectively, over the deepest area of the channel. Lastly, the easternmost station at 5.2 ± 0.05°W (S5) was used to study the most Mediterranean part. To consider the observed latitudinal variation, the northern and southern parts of each station were represented separately at all stations except S1 due to upward and backward routes passing through the same latitudinal interval in the westernmost part. A lower temporal resolution was obtained along the northern part of the channel and at S5 (Supplementary Figure 6) due to technical problems with the measurement equipment and subsequent scarcity of data during the first year of observation. Hence, the harmonic fitting is useful to understand the temporal change during periods of lack of data, but the seasonal trend shows a higher uncertainty in these positions.

FIGURE 3

Figure 3. Seasonal variability of average SST, fCO_2,sw, fCO_2,thermal, fCO_{2,non–thermal}, FCO₂, pH_T, and pH_T,18 at the northern (blue) and southern (red) part of the station S1–S4 in the Strait of Gibraltar.

The latitudinal differences are greater in the most Mediterranean part of the channel. The seasonality of the SST (Figure 3 and Supplementary Figure 6) decreased toward the east through the northern route (in the range between 16–23°C at S1 and 16–20°C at S4) and remained constant with longitude through the southern route (16–23°C). The fCO_2,sw values (Figure 3 and Supplementary Figure 6) were maximum at S1 throughout the year (ranging between 370 μatm during the cold months and 428 μatm during the warm months), and its seasonality remained constant along the southern section between 370 and 415 μatm in phase with temperature and influenced by non-thermal processes. The west-to-east pattern in fCO_2,sw seasonality along the northern section was difficult to characterize due to the lack of values at certain times of the year. However, a slight increase in the seasonal signal was observed from the Atlantic part at S1 (∼57 μatm, between 371 and 428 μatm) to TN at S3 (61 μatm, between 357 and 418 μatm) and a drastic decrease toward the east in GC at S4 (34 μatm, between 362 and 396 μatm).

The pH_T values (Figure 3 and Supplementary Figure 6) were correlated with fCO_2,sw (r² > 0.98 and ρ = −0.99) throughout the region and decreased by 0.001 U μatm^–1 at all stations. The maximum seasonality was found at S1 (8.028–8.080) and in the northern sections of S3 (8.037–8.095). The greatest differences at S2 were observed during the cold months between January and May and could be related with mixing, biological and other non-thermal processes. In GC at S4, the pH_T showed minimum pH_T seasonality ranging from 8.056 to 8.089 in the northern section and from 8.039 to 8.081 in the southern section.

The FCO₂ values (Figure 3 and Supplementary Figure 6) show a slight maximum outgassing and minimum ingassing in the northern part of each station. The outgassing weakened and the ingassing reinforced toward the east, reaching negative FCO₂ values along both the northern and the southern section during the entire first year of monitoring and in the northern part of S4 throughout the study period. A strong CO₂ sink was obtained at all stations during the cold periods, while the system behaved as a weak source during the warm months.

The influence of Mediterranean Sea waters on the distribution of SSS was only observed to the east of GC at S5 through the northern section (Supplementary Figure 7), where the values were maximum throughout the year (∼36.4–36.7). The west-to-east surface circulation from the Northeast Atlantic is reflected in the SSS between ES (S1) and GC (S4): the observed values and seasonality were minimum in the westernmost part at S1 (∼36.3–36.5) and slightly increased through the southern section between S2 and S4 (∼36.3–36.6).

The NC_T (Supplementary Figure 7) decreases between February and September and increases between October and January. The seasonality of C_T and NC_T in ES (S1) (∼30 μmol kg^–1, between 2090 and 2120 μmol kg^–1) increased through the southern section (∼40 μmol kg^–1, between 2090 and 2130 μmol kg^–1) and decreased through the northern section (∼15 μmol kg^–1, between 2100 and 2115 μmol kg^–1) in CS (S2) and TN (S3). In GC at S4, the seasonality decreased in the southern section (∼30 μmol kg^–1, between 2093 and 2123 μmol kg^–1) and increased in the northern section (∼22 μmol kg^–1, between 2098 and 2120 μmol kg^–1). Considering the temporal distribution of Chlorophyll a (Supplementary Figure 7), the seasonal change was reduced and the latitudinal differences increased toward the interior of the channel.

Discussion

High-Resolution Study of the Spatio-Temporal Variability of the CO₂ System

The current study provides a high spatio-temporal resolution evaluation of the surface CO₂ distribution based on high-frequency underway measurements obtained during 2 entire years of observation (February 2019–March 2021). More than 6000 data obtained during those 2 years is the same amount of data available in the SOCAT database for the last 20 years.

According with the high-resolution measurements carried out in this investigation, the seasonal changes in the CO₂ system and fluxes in the Strait of Gibraltar were controlled by the annual cycle of the Atlantic inflow and the position of the AMI. These processes has been shown to control other variables in the region (e.g., Bryden et al., 1994; Bray et al., 1995; Gómez et al., 2000; Skliris and Beckers, 2009). The shallower depth of the AMI favors the deep, cold and CO₂-rich Mediterranean water injection into the upper layers (de la Paz et al., 2009). The AMI elevates toward the northeast throughout the year and reaches minimum depth during late winter (February–March), coinciding with an observed decrease in SST (Figure 2A and Supplementary Figure 1) and increase in SSS and Chlorophyll a (Figure 2B and Supplementary Figure 2). The position of the AMI is also influenced by the tidal amplitude variation and internal wave generation (Gascard and Richez, 1985; Gómez et al., 2001), which favor intermittent vertical mixing events and increase the surface CO₂ during high tides (Wesson and Gregg, 1994; Macias et al., 2006) and especially during the most energetic spring tides (García-Lafuente et al., 2000).

The spatial variability of the CO₂ system during summer is also influenced by mixing processes driven by local climatology factors, as the wind-induced upwelling in the Gulf of Cádiz, the Strait of Gibraltar and the Alboran Sea at local (e.g., Peliz et al., 2009; Gómez-Jakobsen et al., 2019; Bolado-Penagos et al., 2020) and regional scale (e.g., Richez and Kergomard, 1990; Folkard et al., 1997; Stanichny et al., 2005). The signal of the wind-induced upwelling during the warm months was observed in minimum values of SST around CS and between TN and GC (Figure 2). The lowest fCO_2,sw values obtained during summer along the southern section in the westernmost part are explained by the enhancement of the biological uptake in this area due to wind-induced upwelling events at this time of the year (Stanichny et al., 2005). The opposite occurs to the east of TN, where the wind-induced upwelling weakened and higher fCO_2,sw were encountered along the southern route during summer (Figure 2C). The influence of the physical variability in the seasonality of fCO_2,sw in the area were described in term of a second degree function of SST in Eqs 12–16. The spatio-temporal variability observed in the CO₂ system and fluxes and their causes is consistent with previous studies based on datasets obtained from oceanographic surveys at certain locations and times of the year (Dafner et al., 2001; Santana-Casiano et al., 2002; de la Paz et al., 2008, 2009, 2011; Flecha et al., 2019).

The seasonality of the fCO_2,sw values obtained in this study agreed with that obtained from SOCAT data (Supplementary Figure 5) and both followed SST fluctuations. A high coincidence was found in the annual cycle, with an average difference in fCO_2,sw values between both databases of 6.94 ± 16.93 μatm. The fitting of the CanOA-VOS data and the standard error of estimate in the northern and southern section (Eqs 14, 15) and in the entire region (Eq. 12) were slightly enhanced with respect to the SOCAT data (Eq. 13) due to the higher temporal resolution in the entire annual cycle. The combination of the CanOA-VOS and SOCAT data provided an Eq. 16 that can be used to estimate fCO_2,sw in the surface waters of the Strait of Gibraltar. The fCO_2,sw is increasing in the Strait of Gibraltar at a rate of 2.35 ± 0.06 μatm year^–1, which is close to the Northeast Atlantic ESTOC site (1.92 ± 0.24 μatm year^–1) and similar to the Irminger Sea site (2.37 ± 0.49 μatm year^–1), where also vertical mixing processes contributed to the increased observed trend (Bates et al., 2014).

Thermal and Non-thermal Control Over the CO₂ System Variability

The coupling between thermal and non-thermal processes drives the seasonal variation of the CO₂ system throughout the region. These non-thermal processes are mainly biological production (including the organic matter production and remineralization) and physical processes that are strongly coupled affecting the fCO₂ distribution with both positive and negative effects (Takahashi et al., 1993). The T/B ratios (Table 2) show that fCO_2,sw seasonality in the Strait of Gibraltar was mainly controlled by thermal processes in the entire region, coinciding with previous studies in the Strait of Gibraltar (de la Paz et al., 2009) and along the northwest African coast (Curbelo-Hernández et al., 2021). However, the seasonality of fCO_2,thermal and fCO_{2,non–thermal} (Figure 3) and the T/B ratios obtained in the northern and the southern sections of each station show a high spatial variability in thermal and non-thermal coupling. The relevance of non-thermal processes in the temporal variation of fCO_2,sw (Figure 3) was greater along the southern section due to greater influence and intensity of upwelling throughout the African coast (Stanichny et al., 2005). The non-thermal effect was maximum and close to canceling the thermal effect to the south of CS, where the tidal dynamics is enhanced and vertical mixing processes are intense (Echevarría et al., 2002). The non-thermal effect decreased from S2 to S5 through the southern section and increased from S3 to S5 through the northern section, which is related to the high-intensity upwelling in the Northwest Alboran Sea and suppression of vertical mixing to the south (Minas et al., 1991; Echevarría et al., 2002; Gómez-Jakobsen et al., 2019).

TABLE 2

Table 2. Ratios T/B and linear relationship of fCO_2,sw-SST and NC_T-fCO_{2,non–thermal} at the northern and southern part of the selected stations along the Strait of Gibraltar.

The influence of thermal processes on the temporal variation of fCO_2,sw was also analyzed through the fCO_2,sw-SST relationship (Table 2). The increase in fCO_2,sw with temperature through the northern route (9.02 ± 1.99 μatm °C^–1, r² = 0.86, and ρ = 0.93) doubled that obtained through the southern route (4.51 ± 1.66 μatm °C^–1, r² = 0.48, and ρ = 0.69). Thus, the control of thermal processes was even higher in the northern section, while in the southern section the relevance of non-thermodynamic processes increased due to the enhancement of vertical mixing and biological activity (Santana-Casiano et al., 2002). The fCO_2,sw-SST relationship is consistent with the T/B ratios obtained at each station (Table 2). The highest correlation and change of fCO_2,sw with SST was obtained in the westernmost part at S1 and in the northern sections of S2 and S3 (>8.5 μatm °C^–1, r² > 0.67, and ρ > 0.82). In the southern sections of CS and TN, the change and correlation of fCO_2,sw-SST was found to be minimal (r² < 0.28 and ρ < 0.53).

The influence of non-thermal processes on the NC_T change was studied through the relationship NC_T-fCO_{2,non–thermal} (Table 2). The NC_T increased with respect to fCO_{2,non–thermal} by 0.50 ± 0.03 μmol kg^–1 μatm^–1 (r² > 0.92 and ρ > 0.96) throughout the region, but decreased with respect to fCO_2,thermal by 0.22 ± 0.06 (r² = 0.82 and ρ = −0.91) and 0.34 ± 0.06 (r² = 0.83 and ρ = −0.91) μmol kg^–1 μatm^–1 through the northern and southern routes, respectively. The NC_T showed a direct relationship and correlated with fCO_{2,non–thermal} (between 0.46 and 0.52 μmol kg^–1 μatm^–1; r² > 0.85 and ρ > 0.92) at all stations. Hence, the NC_T variation was mainly driven by non-thermal processes (Santana-Casiano et al., 2002; de la Paz et al., 2008, 2009).

Seasonal Variability of NC_T and Relative Influence of Air-Sea Exchange, Biological Activity, Horizontal Transport, and Vertical Mixing

The annual cycle of C_T was studied at each of the selected stations by considering its temporal change between consecutive months. The C_T was normalized to a constant salinity of 36.4 to remove the effect of precipitation and evaporation and study the other non-conservative processes involved in the surface variation of C_T (i.e., primary production, oxidation of organic matter, precipitation/dissolution of CaCO₃, air-sea CO₂ fluxes, and physical processes) (e.g., Chen and Pytkowicz, 1979; Takahashi et al., 1993; Wanninkhof and Feely, 1998; Lee et al., 2000).

The results allow to establish the relative contribution of biological activity, air-sea exchange, horizontal transport, and vertical mixing (Eq. 8) on the increase and depletion of NC_T in surface waters (Figure 4 and Table 3). The greatest seasonal changes of NC_T were found at S1 and throughout the southern section, especially at S2 and S3 (Table 3). A high spatial variability was found in the seasonality of NC_T. The depletion of NC_T in the southwestern part occurred from February–March to September, while it was difficult to establish a seasonal pattern through the northern section and in the easternmost part due to the high variability of the surface inflow in the tidal time-scale (de la Paz et al., 2009).

FIGURE 4

Figure 4. Monthly change in NC_T (ΔNC_T/dt; black) and relative contribution of net community production (ΔNC_T/dt| _NCP; green), air-sea exchange (ΔNC_T/dt| _EX; blue), horizontal transport (ΔNC_T/dt| _TRSP; purple), and vertical mixing (ΔNC_T/dt| _MIX; orange) in an entire annual cycle. The temporal variation of NC_T was evaluated in S1 (A) and in the northernmost and southernmost part of S2 (B), S3 (C), S4 (D), and S5 (E). The net contribution of each term in the increase and depletion of NC_T is shown in Table 4.

TABLE 3

Table 3. Seasonal change of NC_T (ΔNC_T/dt) and net contribution of the net community production (ΔNC_T/dt| _NCP; green), air-sea exchange (ΔNC_T/dt| _EX; blue), horizontal transport (ΔNC_T/dt| _TRSP; purple), and vertical mixing (ΔNC_T/dt| _MIX; orange) in the increase and depletion of NC_T at stations S1–S5.

The seasonal change of NC_T was controlled by NCP processes at most of the stations except in the northern part at S2 (0.4–10.6%; Figure 4B), where the air-sea exchange and vertical mixing dominated (62.7–68.4 and 29.7–36.7%, respectively) and horizontal transport influenced (8.0–15.6%). It is related to the generation of internal bores at CS (Boyce, 1975; Armi and Farmer, 1988; Bruno et al., 2002) and the minimum intensity of the wind-induced upwelling toward the northwest of the channel (Peliz et al., 2009). The biological effect was especially significant along the African coastal upwelling (Stanichny et al., 2005) in the southern part of S2 (81.8–83.9%) and S3 (94.2–98.4%). The upwelling effect in the Northwest Alboran Sea (Minas et al., 1991; Echevarría et al., 2002; Gómez-Jakobsen et al., 2019) was detected at the easternmost stations (Figures 4D,E). The change in NC_T explained by NCP exceeded the total change of NC_T by 24.4–30.0 and 21.9–69.3% in the northern part at S4 and S5, respectively. The biological effect in these locations was counteracted by the contribution of air-sea exchange (12.3–33.6% in the northern part at S4 and S5 and 69.04–69.2% in the southern part at S5) and the influence of vertical mixing and horizontal transport (<12 and <10%, respectively). The contribution of NCP is strongly linked to the vertical mixing processes in this region, which inject deeper and remineralized CO₂-rich water in surface layers. Positive contributions of vertical mixing indicated that the NC_T increase in the surface layer due to the deep-water injection, while negative values indicate that the concentration of NC_T on the surface is higher than in depths layers due to a minimal biological production on the surface. In some periods, the amount of injected CO₂ to the surface is not counterbalanced by primary production and positive values contributed to the observed monthly variability. Accordingly, the seasonal and spatial variation of the AMI and the deep water injection processes in the surface layers drove by the wind (wind-induced upwelling) through the southernmost part and by tides (tidal-induced upwelling) through the northernmost part could account for the observed positive and negative contributions of both NCP and vertical mixing. The contribution of vertical mixing was especially significant in certain locations due to the advection of C_T in upper layers toward the east of CS during tidal-induced upwelling events (de la Paz et al., 2008). The diffusion processes of C_T through the pycnocline increased its surface concentration by an average of 3.55 ± 2.69 mmol m^–3 during the entire annual cycle throughout the region. The vertical entrainment increased the surface C_T by an average of 0.58 ± 0.52 mmol m^–3 between August and January and was canceled between February and July due to the elevation of the mixed layer and subsequent increase in stratification. About 5.7–5.8% of the seasonal change in NC_T at S1 (Figure 4C) was explained by vertical mixing processes. The vertical mixing inversely affected the change of NC_T at S3 and S5. During the periods of increase in NC_T, between 0.02 and 9.10% of the NC_T was removed from the surface by vertical mixing processes, while during periods of depletion, between 2.90 and 8.20% of the surface NC_T was injected from deep waters. The contribution of vertical transport was maximum in the northern part of S2 (∼29.7–36.7%) and at S4 (∼10.3–11.9%) (Figures 4B,D) due to the effect of internal waves in CS (Bruno et al., 2002; Macias et al., 2006) and the elevation of the pycnocline in the most Mediterranean part (Echevarría et al., 2002). The influence of the seasonality of the wind-induced upwelling on the change of NC_T by vertical mixing was detected in S2, where a difference of 1.34 mmol m^–3 between increase and depletion periods was found.

The horizontal transport in the west-to-east direction decreased the surface NC_T (0.8–15.6%) in the entire annual cycle due to the lower A_T and C_T content of the Atlantic inflow compared to the Mediterranean outflow (de la Paz et al., 2008). An increase in NC_T by horizontal transport was only observed at S1 (0.4%) and in the northern part of S4 (2.6%) and S5 (9.6%).

The average seasonal change of NC_T for all selected locations across the Strait of Gibraltar was 29.0 ± 9.96 mmol m^–3. Biological processes controlled the seasonal change of surface NC_T (93.5–95.6%), increasing its concentration by 27.72 mmol m^–3 in periods of higher remineralization influence and decreasing it by 27.11 mmol m^–3 when production was the controlling process. The contribution of air-sea CO₂ exchange (0.9–3.3%) was higher during the increasing periods due to the strong behavior as a sink of the channel (0.94 and 0.27 mmol m^–3 during the increase and depletion periods, respectively). Vertical mixing increased the NC_T by 1.32 mmol m^–3 and decreased it by 1.02 mmol m^–3, with a total contribution of 4.1%. Finally, the effect of horizontal transport was minimal (<2.8%) and decreased the NC_T by 1.08 and 0.52 mmol m^–3 in both total increase and depletion periods, respectively.

CO₂ Fluxes in the Strait of Gibraltar

The temporal distribution of FCO₂ calculated from observed and computed data allowed the characterization of the air-sea CO₂ exchange in the Strait of Gibraltar throughout the study period. The monthly average FCO₂ values (Supplementary Figure 8) proved the strong behavior of the region as a CO₂ sink from November to July and as a weak source between August and October. This seasonal pattern is consistent with previous studies in the Strait of Gibraltar (de la Paz et al., 2009, 2011) and along the coastal to open-ocean transitional area of the Northeast Atlantic (Curbelo-Hernández et al., 2021). In the northern section, the CO₂ sink behavior remained throughout the first year of observation (February 2019–February 2020) even in the warm months, while a late-summer CO₂ source behavior was found in the second year of observation (March 2020–March 2021). In contrast, the southern section behaved as a late-summer CO₂ source in both entire years considered.

The FCO₂ values show a direct relationship with ΔfCO₂ in both the northern and the southern sections [0.11 (r² = 0.87 and ρ = 0.93) and 0.14 (r² = 0.72 and ρ = 0.84) mmol m^–2 d^–1 μatm^–1, respectively]. Despite the important role of the wind speed and direction controlling the seasonal and spatial distribution and intensity of upwelling events across coastal areas in the Strait of Gibraltar (Stanichny et al., 2005), its short-term variability is not significant on an annual time scale in this region (de la Paz et al., 2011) and specially along the northern section. The difference between the averages FCO₂ calculated with satellite and measured wind speed data was 0.14 mol C m^–2 year^–1. The changes in wind speed influence the intensity of CO₂ exchange between low atmosphere and surface ocean and could modified the physical properties of the upper layers, but do not determine the formation of CO₂ source/sink. Northeasterly winds bring on the deep-water injection in upper layers and subsequent cooling of the surface seawater in the southwest area, favoring the solubility of atmospheric CO₂. Accordingly, a more intense annual CO₂ sink was found in this area compared to the northwest area (Figure 5). The assessment of the spatio-temporal changes of FCO₂ in this study show that the formation of the seasonal CO₂ sink and source is controlled by the ΔfCO₂ variations influenced by changes in the physical properties of the surface ocean, which were driven by the local climatology and the unique hydrology of the channel (Supplementary Figure 8). Therefore, Eq. 16 which provides fCO_2,sw data for the area using SST as unknown variable (available from satellite data) together with satellite winds, can be used to estimate long term CO₂ fluxes in the Strait of Gibraltar.

FIGURE 5

Figure 5. Annual average of FCO₂ and ΔfCO₂ every 0.05° longitude calculated along the northern and southern sections and in the entire region.

The annual negative values of FCO₂ obtained along both the northern and the southern section show that the Strait of Gibraltar behaved as a net CO₂ sink (Figure 5). The annual average FCO₂, ΔfCO₂, and wind speed were calculated at all the selected stations, along the northern and southern sections and in the entire region (Table 4). The average ingassing was found to be maximum in the eastern area (Figure 5). The total annual FCO₂ in most of the Northeast Atlantic and Mediterranean Sea parts of the channel was −1.02 and −1.30 mol C m^–2 year^–1, respectively, which is related with lower SST values observed throughout the year and the enhancement of the biological uptake toward the northeastern part of the channel (Echevarría et al., 2002) that increased the ΔfCO₂.

TABLE 4

Table 4. Annual average FCO₂, ΔfCO₂, and wind speed obtained during the first and second year of observation (February 2019–February 2020 and March 2020–March 2021, respectively) and throughout the study period at stations S1–S5, along the northern and southern sections and in the entire region.

The net CO₂ sink obtained in this study between 2019 and 2021 for the entire region was −0.98 mol C m^–2 year^–1, which falls in the range of the first annual FCO₂ provided by Santana-Casiano et al. (2002) (−2.5 mol C m^–2 year^–1) and the annual FCO₂ reported in most recent studies of de la Paz et al. (2009) (−0.28 mol C m^–2 year^–1) and de la Paz et al. (2011) (−0.02 ± 0.13 mol C m^–2 year^–1). The annual FCO₂ was calculated in these studies from limited datasets and using the parameterization of Wanninkhof (1992) (W92). When applying this parameterization to the data collected, an annual value of FCO₂ of −1.21 mol C m^–2 year^–1 was obtained. This average ingassing is more than two times higher than the values reported with W92 for the Gulf of Cadiz (−0.4 mol C m^–2 year^–1; Borges et al., 2006) and for the Northwest African Atlantic coast (−0.59 mol C m^–2 year^–1, using W92; Curbelo-Hernández et al., 2021). However, it falls within the annual FCO₂ estimated by Chen et al. (2013) for the Northeast Atlantic (−1.04 mol C m^–2 year^–1) and for the Iberian upwelling (−1.33 mol C m^–2 year^–1). The net CO₂ sink reported in this study is consistent with that obtained from SOCAT data referenced to 2019 (−0.97 mol C m^–2 year^–1; using W92).

The net FCO₂ for the entire area (1650 km²) was −7.12 Gg CO₂ year^–1 (1 Gg = 10⁹ g) (−1.94 Gg C year^–1). This result is in agreement with previous studies of the behavior of the coastal regions and continental shelves as a net CO₂ sink (e.g., Borges et al., 2006, 2005; Chen and Borges, 2009; Laruelle et al., 2010, 2013, 2018; Cai et al., 2006; Cai, 2011; Chen et al., 2013) and can be integrated into the estimated net CO₂ sink for the European coastal regions (−68.1 × 10³ Gg C year^–1; Borges et al., 2006) and global continental shelves (−0.4 × 10⁶ Gg C year^–1; Chen et al., 2013).

Conclusion

The spatio-temporal variability of the CO₂ system and fluxes in the Strait of Gibraltar was assessed with a high-resolution dataset in this study. The results obtained between February 2019 and March 2021 provided the seasonal and spatial variability of the surface CO₂. Although seasonal temperature fluctuation was the main driver of the observed changes in the CO₂ system in this region, the temporal and spatial distribution of the variables was strongly influenced by physical and biological factors oscillations due to the seasonal change in the depth of the AMI layer, the wind-induced upwelling during northeasterly winds and the tidal-induced pulsating upwelling throughout the year and especially during spring tides. The elevation of the AMI layer toward the northeast and the wind-induced upwelling on the African coast toward the west of TN favors the vertical mixing and enhances the physical–biological coupling in the northeast and southwest part of the channel. In these areas of lower SST and higher SSS, the change in fCO_2,sw with temperature was lower due to the increase in the importance of the non-thermal effects. The annual cycle of surface fCO_2,sw was described in this study (Eqs 14, 15) by considering the seasonality of SST and its spatial changes controlled by the local climatology and the unique hydrology of the channel. Equation 16 allows to compute past, present, and future surface fCO_2,sw values using satellite SST values. The fCO_2,sw was found to increase with temperature by 9.02 ± 1.99 μatm °C (r² = 0.86 and ρ = 0.93) and 4.51 ± 1.66 μatm °C (r² = 0.48 and ρ = 0.69) in the northern and southern sections, respectively.

The seasonality of NC_T was mainly controlled by NCP processes (93.5–95.6%), while the contribution of air-sea exchange (0.9–3.3%), horizontal transport (1.8–3.7%), and vertical mixing (3.5–4.6%) was minimal. The seasonality of NC_T was higher in the region most influenced by wind-induced upwelling to the southwest of the channel (42.65 and 45.44 mmol m^–3 calculated in the southern part of S2 and S3, respectively). In the areas of maximum upwelling, the contribution of NCP compared to that of air-sea exchange was found to be maximum (>81.0 and <13.1%, respectively).

In terms of air-sea CO₂ exchange, the Strait of Gibraltar behaved as a strong CO₂ sink during the cold months (between −7 and −2 mmol m^–2 year^–1) and as a weak source during the warm months (<3 mmol m^–2 year^–1). Considering the high variability in the CO₂ sink and integrating the entire area (6.1–5.1°W, 1650 km²) an average CO₂ flux of −7.12 Gg CO₂ year^–1 (−1.94 Gg C year^–1) was determined. The net CO₂ sink was higher toward the east of TN (−1.30 mol C m^–2 year^–1), where the depth of the AMI layer is minimum and both solubility and biological uptake were enhanced by the injection of deep, cold, and nutrient-rich water throughout the year. The net CO₂ sink weakened toward the west of TN (−1.02 mol C m^–2 year^–1), where higher annual SST values were obtained due to vertical mixing processes are limited to periods of wind-induced upwelling.

The present study shows that VOS lines are a powerful tool to study the carbon cycle in global and regional surface ocean. The continuous monitoring of the surface CO₂ by VOS lines is especially important in poorly known and highly variable areas such as the Strait of Gibraltar, where the observational data from oceanographic surveys are limited and continuous monitoring through other observational strategies is difficult to perform due to maritime traffic and socio-political factors. We provide in this study clear evidence of the complex processes related to the unique hydrology and local climatology of the transition region of the Strait of Gibraltar involved in the surface distribution of the CO₂ system variables. The air-sea CO₂ exchange and the surface inorganic carbon concentration, as well as the processes that control them, were assessed and quantified with an improved degree of certainty. These findings enhance the understanding of the surface CO₂ system and air-sea exchange in marginal and semi-enclosed seas and can be used to further compute how much of this atmospheric CO₂ entering in the surface ocean is finally transported into the ocean interior.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.socat.info/.

Author Contributions

MG-D, JMS-C, and AGG installed and maintained the equipment in the VOS line. DC-H performed the data treatments and MATLAB^® routines. All authors made significant contributions toward the writing of the manuscript.

Funding

This study was supported by the Canary Islands Government and the Loro Parque Foundation through the CanBIO project, CanOA subproject (2019–2022), and the CARBOCAN agreement (Consejería de Transición Ecológica, Lucha contra el Cambio Climático y Planificación Territorial, and Gobierno de Canarias).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We would like to thank the RENATE P ship owner, the NISA-Marítima company and the captains and crew members for the support during this collaboration. Special thanks to the technician Adrian Castro-Alamo for biweekly equipment maintenance and discrete sampling of total alkalinity aboard the ship. We thank the editor and reviewers for their constructive comments and suggestions, which helped to improve the quality of the manuscript. The VOS line will be part of the Spanish contribution to the Integrated Carbon Observation System (ICOS), European Research Infrastructure starting in 2021. We also thank the Surface Ocean CO₂ Atlas (SOCAT), an international effort, endorsed by the International Ocean Carbon Coordination Project (IOCCP), the Surface Ocean Lower Atmosphere Study (SOLAS), and the Integrated Marine Biosphere Research (IMBeR) program, to deliver a uniformly quality-controlled surface ocean CO₂ database. The many researchers and funding agencies responsible for the collection of data and quality control are thanked for their contributions to SOCAT.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2021.745304/full#supplementary-material

Footnotes

1. ^ https://marine.copernicus.eu/
2. ^ https://www.socat.info/

References

Allen, R. G., Pereira, L. S., Raes, D., and Smith, M. (1998). Crop Evapotranspiration -Guidelines for Computing Crop Water Requirements. FAO Irrigation and Drainage Paper 56. Rome: Food and Agriculture Organization.

Google Scholar

Alonso del Rosario, J. J., Mejías, M. B., and Vázquez-Escobar, A. (2003). The influence of tidal hydrodynamic conditions on the generation of lee waves at the main sill of the Strait of Gibraltar. Deep Sea Res. 1 Oceanogr. Res. Pap. 50, 1005–1021. doi: 10.1016/S0967-0637(03)00097-9

CrossRef Full Text | Google Scholar

Armi, L., and Farmer, D. M. (1988). The flow of Mediterranean water through the Strait of Gibraltar. the flow of Atlantic water through the Strait of Gibraltar. Prog. Oceanogr. 21, 1–105.

Google Scholar

Bates, N. R., Astor, Y. M., Church, M. J., Currie, K., Dore, J. E., González-Dávila, M., et al. (2014). A time-series view of changing surface ocean chemistry due to ocean uptake of anthropogenic CO2 and ocean acidification. Oceanography 27, 126–141. doi: 10.5670/oceanog.2014.16

CrossRef Full Text | Google Scholar

Bégovic, M., and Copin-Montégut, C. (2002). Processes controlling annual variations in the partial pressure of CO2 in surface waters of the central northwestern Mediterranean Sea (Dyfamed site). Deep Sea Res. 2 Top. Stud. Oceanogr. 49, 2031–2047. doi: 10.1016/S0967-0645(02)00026-7

CrossRef Full Text | Google Scholar

Bolado-Penagos, M., González, C. J., Chioua, J., Sala, I., Gomiz-Pascual, J. J., Vázquez-Escobar, A., et al. (2020). Submesoscale processes in the coastal margins of the Strait of Gibraltar. The Trafalgar–Alboran connection. Prog. Oceanogr. 181:102219. doi: 10.1016/j.pocean.2019.102219

CrossRef Full Text | Google Scholar

Borges, A. V., Delille, B., and Frankignoulle, M. (2005). Budgeting sinks and sources of CO2 in the coastal ocean: diversity of ecosystem counts. Geophys. Res. Lett. 32, 1–4. doi: 10.1029/2005GL023053

CrossRef Full Text | Google Scholar

Borges, A. V., Schiettecatte, L. S., Abril, G., Delille, B., and Gazeau, F. (2006). Carbon dioxide in European coastal waters. Estuar. Coast. Shelf Sci. 70, 375–387. doi: 10.1016/j.ecss.2006.05.046

CrossRef Full Text | Google Scholar

Boyce, F. M. (1975). Internal waves in the Straits of Gibraltar. Deep. Res. Oceanogr. Abstr. 22, 597–610. doi: 10.1016/0011-7471(75)90047-9

CrossRef Full Text | Google Scholar

Bray, N. A., Ochoa, J., and Kinder, T. H. (1995). The role of the interface in exchange through the Strait of Gibraltar. J. Geophys. Res. 100, 10755–10776. doi: 10.1029/95JC00381

CrossRef Full Text | Google Scholar

Broecker, W. S., and Peng, T. H. (1982). Tracers in the Sea. Palisades, NY: Eldigio Press, 1–690.

Google Scholar

Bruno, M., Alonso, J. J., Cózar, A., Vidal, J., Ruiz-Cañavate, A., Echevarría, F., et al. (2002). The boiling-water phenomena at Camarinal Sill, the Strait of Gibraltar. Deep. Sea Res. 2 Top. Stud. Oceanogr. 49, 4097–4113. doi: 10.1016/S0967-0645(02)00144-3

CrossRef Full Text | Google Scholar

Bryden, H. L., Candela, J., and Kinder, T. H. (1994). Exchange through the Strait of Gibraltar. Prog. Oceanogr. 33, 201–248. doi: 10.1016/0079-6611(94)90028-0

CrossRef Full Text | Google Scholar

Bryden, H. L., and Kinder, T. H. (1991). Steady two-layer exchange through the Strait of Gibraltar. Deep Sea Res. A Oceanogr. Res. Pap. 38, S445–S463. doi: 10.1016/s0198-0149(12)80020-3

CrossRef Full Text | Google Scholar

Cai, W. J. (2011). Estuarine and coastal ocean carbon paradox: CO₂ sinks or sites of terrestrial carbon incineration? Ann. Rev. Mar. Sci. 3, 123–145. doi: 10.1146/annurev-marine-120709-142723

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, W. J., Dai, M., and Wang, Y. (2006). Air-sea exchange of carbon dioxide in ocean margins: a province-based synthesis. Geophys. Res. Lett. 33:L12603. doi: 10.1029/2006GL026219

CrossRef Full Text | Google Scholar

Candela, J., Winant, C., and Ruiz, A. (1990). Tides in the Strait of Gibraltar. J. Geophys. Res. 95, 7313–7335. doi: 10.1029/JC095iC05p07313

CrossRef Full Text | Google Scholar

Candela, J., Winant, C. D., and Bryden, H. L. (1989). Meteorologically forced subinertial flows through the Strait of Gibraltar. J. Geophys. Res. Ocean. 94, 12667–12679. doi: 10.1029/JC094iC09p12667

CrossRef Full Text | Google Scholar

Chen, C. T. A., and Borges, A. V. (2009). 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. 2: Top. Stud. Oceanogr. 56, 578–590. doi: 10.1016/j.dsr2.2009.01.001

CrossRef Full Text | Google Scholar

Chen, C. T. A., Huang, T. H., Chen, Y. C., Bai, Y., He, X., and Kang, Y. (2013). Air-sea exchanges of coin the world’s coastal seas. Biogeosciences 10, 6509–6544. doi: 10.5194/bg-10-6509-2013

CrossRef Full Text | Google Scholar

Chen, C. T. A., and Pytkowicz, R. M. (1979). On the total CO2-titration alkalinity-Oxygen system in the Pacific Ocean [5]. Nature 281, 362–365. doi: 10.1038/281362a0

CrossRef Full Text | Google Scholar

Curbelo-Hernández, D., González-Dávila, M., González, A. G., González-Santana, D., and Santana-Casiano, J. M. (2021). CO2 fluxes in the Northeast Atlantic Ocean based on measurements from a surface ocean observation platform. Sci. Total Environ. 775:145804. doi: 10.1016/j.scitotenv.2021.145804

PubMed Abstract | CrossRef Full Text | Google Scholar

Dafner, E., González-Dávila, M., Santana-Casiano, J. M., and Sempéré, R. (2001). Total organic and inorganic carbon exchange through the Strait of Gibraltar in September 1997. Deep Sea Res. 1 Oceanogr. Res. Pap. 48, 1217–1235. doi: 10.1016/S0967-0637(00)00064-9

CrossRef Full Text | Google Scholar

de la Paz, M., Debelius, B., Macías, D., Vázquez-Escobar, A., Gómez-Parra, A., and Forja, J. M. (2008). Tidal-induced inorganic carbon dynamics in the Strait of Gibraltar. Cont. Shelf Res. 28, 1827–1837. doi: 10.1016/j.csr.2008.04.012

CrossRef Full Text | Google Scholar

de la Paz, M., Gómez-Parra, A., and Forja, J. (2009). Seasonal variability of surface fCO2 in the Strait of Gibraltar. Aquat. Sci. 71, 55–64. doi: 10.1007/s00027-008-8060-y

CrossRef Full Text | Google Scholar

de la Paz, M., Huertas, E. M., Padín, X. A., Gónzalez-Dávila, M., Santana-Casiano, J. M., Forja, J. M., et al. (2011). Reconstruction of the seasonal cycle of air-sea CO2 fluxes in the Strait of Gibraltar. Mar. Chem. 126, 155–162. doi: 10.1016/j.marchem.2011.05.004

CrossRef Full Text | Google Scholar

Denman, K. L., and Gargett, A. E. (1983). Time and space scales of vertical mixing and advection of phytoplankton in the upper ocean. Limnol. Oceanogr. 28, 801–815.

Google Scholar

Dickson, A. G. (1990). Standard potential of the reaction: AgCl(s) + 1 2H2(g) = Ag(s) + HCl(aq), 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. doi: 10.1016/0021-9614(90)90074-Z

CrossRef Full Text | Google Scholar

Dickson, A. G., Sabine, C. L., and Chistian, J. R. (2007). Guide to best practices for ocean CO2 measurements. PICES Special Publ. 3:191.

Google Scholar

Echevarría, F., García-Lafuente, J., Bruno, M., Gorsky, G., Goutx, M., González, N., et al. (2002). Physical-biological coupling in the Strait of Gibraltar. Deep Sea Res. 2 Top. Stud. Oceanogr. 49, 4115–4130. doi: 10.1016/S0967-0645(02)00145-5

CrossRef Full Text | Google Scholar

Fasham, M. J. R., Ducklow, H. W., and McKelvie, S. M. (1990). A nitrogen-based model of plankton dynamics in the oceanic mixed layer. J. Mar. Res. 48, 591–639. doi: 10.1357/002224090784984678

CrossRef Full Text | Google Scholar

Flecha, S., Pérez, F. F., Murata, A., Makaoui, A., and Huertas, I. E. (2019). Decadal acidification in Atlantic and Mediterranean water masses exchanging at the Strait of Gibraltar. Sci. Rep. 9:15533. doi: 10.1038/s41598-019-52084-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Folkard, A. M., Davies, P. A., Fiúza, A. F. G., and Ambar, I. (1997). Remotely sensed sea surface thermal patterns in the Gulf of Cadiz and the Strait of Gibraltar: variability, correlations and relationships with the surface wind field. J. Geophys. Res. Ocean. 102, 5669–5683. doi: 10.1029/96JC02505

CrossRef Full Text | Google Scholar

Frankignoulle, M., and Borges, A. V. (2001). European continental shelf as a significant sink for atmospheric carbon dioxide. Glob. Biogeochem. Cycles 15, 569–576. doi: 10.1029/2000GB001307

CrossRef Full Text | Google Scholar

García-Lafuente, J., Sánchez-Román, A., Naranjo, C., and Sánchez-Garrido, J. C. (2011). The very first transformation of the Mediterranean outflow in the Strait of Gibraltar. J. Geophys. Res. Ocean. 116:7010. doi: 10.1029/2011JC006967

CrossRef Full Text | Google Scholar

García-Lafuente, J., Vargas, J. M., Plaza, F., Sarhan, T., Candela, J., and Bascheck, B. (2000). Tide at the eastern section of the Strait of Gibraltar. J. Geophys. Res. Ocean. 105, 14197–14213. doi: 10.1029/2000jc900007

CrossRef Full Text | Google Scholar

Gascard, J. C., and Richez, C. (1985). Water masses and circulation in the Western Alboran Sea and in the Strait of Gibraltar. Prog. Oceanogr. 15, 157–216. doi: 10.1016/0079-6611(85)90031-X

CrossRef Full Text | Google Scholar

Gómez, F., Echevarría, F., García, C. M., Prieto, L., Ruiz, J., Reul, A., et al. (2000). Microplankton distribution in the Strait of Gibraltar: coupling between organisms and hydrodynamic structures. J. Plankton Res. 22, 603–617. doi: 10.1093/plankt/22.4.603

CrossRef Full Text | Google Scholar

Gómez, F., Gorsky, G., Striby, L., Vargas, J. M., Gonzalez, N., Picheral, M., et al. (2001). Small-scale temporal variations in biogeochemical features in the Strait of Gibraltar, Mediterranean side–the role of NACW and the interface oscillation. J. Mar. Syst. 30, 207–220. doi: 10.1016/S0924-7963(01)00059-8

CrossRef Full Text | Google Scholar

Gómez-Jakobsen, F. J., Mercado, J. M., Cortés, D., Yebra, L., and Salles, S. (2019). A first description of the summer upwelling off the Bay of Algeciras and its role in the northwestern Alboran Sea. Estuar. Coast. Shelf Sci. 225:106230. doi: 10.1016/j.ecss.2019.05.012

CrossRef Full Text | Google Scholar

González-Dávila, M., Santana-Casiano, J. M., Rueda, M. J., Llinás, O., and González-Dávila, E. F. (2003). Seasonal and interannual variability of sea-surface carbon dioxide species at the European Station for Time Series in the Ocean at the Canary Islands (ESTOC) between 1996 and 2000. Glob. Biogeochem. Cycles 17:1076. doi: 10.1029/2002gb001993

CrossRef Full Text | Google Scholar

Gruber, N., Keeling, C. D., and Stocker, T. F. (1998). Carbon-13 constraints on the seasonal inorganic carbon budget at the BATS site in the northwestern Sargasso Sea. Deep Sea Res. 1 Oceanogr. Res. Pap. 45, 673–717. doi: 10.1016/S0967-0637(97)00098-8

CrossRef Full Text | Google Scholar

Huertas, I. E., Ríos, A. F., García-Lafuente, J., Makaoui, A., Rodríguez-Gálvez, S., Sánchez-Román, A., et al. (2009). Anthropogenic and natural CO2 exchange through the Strait of Gibraltar. Biogeosciences 6, 647–662. doi: 10.5194/bg-6-647-2009

CrossRef Full Text | Google Scholar

IPCC (2007). Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, eds S. Solomon, D. Qin, M. Manning, Z. Chen, M. Marquis, K. B. Averyt, et al. Cambridge: Cambridge University Press.

Google Scholar

IPCC (2021). Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, eds V. Masson-Delmotte, P. Zhai, A. Pirani, S. L. Connors, C. Péan, S. Berger, et al. Cambridge: Cambridge University Press.

Google Scholar

La Violette, P. E., and Arnone, R. A. (1988). A tide-generated internal waveform in the western approaches to the Strait of Gibraltar. J. Geophys. Res. Ocean. 93, 15653–15667. doi: 10.1029/JC093iC12p15653

CrossRef Full Text | Google Scholar

Lacombe, H., and Richez, C. (1982). The regime of the Strait of Gibraltar. Oceanogr. Ser. 34, 13–73. doi: 10.1016/S0422-9894(08)71237-6

CrossRef Full Text | Google Scholar

Laruelle, G. G., Cai, W. J., Hu, X., Gruber, N., Mackenzie, F. T., and Regnier, P. (2018). Continental shelves as a variable but increasing global sink for atmospheric carbon dioxide. Nat. Commun. 9:454. doi: 10.1038/s41467-017-02738-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Laruelle, G. G., Dürr, H. H., Lauerwald, R., Hartmann, J., Slomp, C. P., Goossens, N., et al. (2013). Global multi-scale segmentation of continental and coastal waters from the watersheds to the continental margins. Hydrol. Earth Syst. Sci. 17, 2029–2051. doi: 10.5194/hess-17-2029-2013

CrossRef Full Text | Google Scholar

Laruelle, G. G., Dürr, H. H., Slomp, C. P., and Borges, A. V. (2010). Evaluation of sinks and sources of CO2 in the global coastal ocean using a spatially-explicit typology of estuaries and continental shelves. Geophys. Res. Lett. 37:L15607. doi: 10.1029/2010GL043691

CrossRef Full Text | Google Scholar

Lee, K., Kim, T. W., Byrne, R. H., Millero, F. J., Feely, R. A., and Liu, Y. M. (2010). The universal ratio of boron to chlorinity for the North Pacific and North Atlantic oceans. Geochim. Cosmochim. 74, 1801–1811. doi: 10.1016/j.gca.2009.12.027

CrossRef Full Text | Google Scholar

Lee, K., Tong, L. T., Millero, F. J., Sabine, C. L., Dickson, A. G., Goyet, C., et al. (2006). Global relationships of total alkalinity with salinity and temperature in surface waters of the world’s oceans. Geophys. Res. Lett. 33, 1–5. doi: 10.1029/2006GL027207

CrossRef Full Text | Google Scholar

Lee, K., Wanninkhof, R., Feely, R. A., Millero, F. J., and Peng, T. H. (2000). Global relationships of total inorganic carbon with temperature and nitrate in surface seawater. Glob. Biogeochem. Cycles 14, 979–994. doi: 10.1029/1998GB001087

CrossRef Full Text | Google Scholar

Lueker, T. J., Dickson, A. G., and Keeling, C. D. (2000). 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. doi: 10.1016/S0304-4203(00)00022-0

CrossRef Full Text | Google Scholar

Lüger, H., Wallace, D. W. R., Körtzinger, A., and Nojiri, Y. (2004). The pCO2 variability in the midlatitude North Atlantic Ocean during a full annual cycle. Glob. Biogeochem. Cycles 18:GB3023. doi: 10.1029/2003GB002200

CrossRef Full Text | Google Scholar

Macias, D., García, C. M., Echevarría, F., Vázquez-Escobar, A., and Mejías, M. B. (2006). Tidal induced variability of mixing processes on Camarinal Sill (Strait of Gibraltar): a pulsating event. J. Mar. Syst. 60, 177–192. doi: 10.1016/j.jmarsys.2005.12.003

CrossRef Full Text | Google Scholar

Macías, D., Martin, A. P., García-Lafuente, J., García, C. M., Yool, A., Bruno, M., et al. (2007). Analysis of mixing and biogeochemical effects induced by tides on the Atlantic-Mediterranean flow in the Strait of Gibraltar through a physical-biological coupled model. Prog. Oceanogr. 74, 252–272. doi: 10.1016/j.pocean.2007.04.006

CrossRef Full Text | Google Scholar

Mackenzie, F. T., Lerman, A., and Ver, L. M. B. (1998). Role of the continental margin in the global carbon balance during the past three centuries. Geology 26, 423–426. doi: 10.1130/0091-76131998026<0423:ROTCMI<2.3.CO;2

CrossRef Full Text | Google Scholar

Minas, H. J., Coste, B., Le Corre, P., Minas, M., and Raimbault, P. (1991). Biological and geochemical signatures associated with the water circulation through the Strait of Gibraltar and in the western Alboran Sea. J. Geophys. Res. 96, 8755–8771. doi: 10.1029/91JC00360

CrossRef Full Text | Google Scholar

Mintrop, L., Pérez, F. F., González-Dávila, M., Santana-Casiano, J. M., and Körtzinger, A. (2000). Alkalinity determination by potentiometry: intercalibration using three different methods. Ciencias Mar. 26, 23–37. doi: 10.7773/cm.v26i1.573

CrossRef Full Text | Google Scholar

Muller-Karger, F. E., Varela, R., Thunell, R., Luerssen, R., Hu, C., and Walsh, J. J. (2005). The importance of continental margins in the global carbon cycle. Geophys. Res. Lett. 32, 1–4. doi: 10.1029/2004GL021346

CrossRef Full Text | Google Scholar

Peliz, A., Teles-Machado, A., Marchesiello, P., Dubert, J., and García-Lafuente, J. (2009). Filament generation off the Strait of Gibraltar in response to Gap winds. Dyn. Atmos. Ocean. 46, 36–45. doi: 10.1016/j.dynatmoce.2008.08.002

CrossRef Full Text | Google Scholar

Peng, T. H., Takahashi, T., Broecker, W. S., and Olafsson, J. (1987). Seasonal variability of carbon dioxide, nutrients and oxygen in the northern North Atlantic surface water: observations and a model. Tellus B Chem. Phys. Meteorol. 39, 439–458. doi: 10.3402/tellusb.v39i5.15361

CrossRef Full Text | Google Scholar

Pierrot, D., Neill, C., Sullivan, K., Castle, R., Wanninkhof, R., Lüger, H., et al. (2009). Recommendations for autonomous underway pCO2 measuring systems and data-reduction routines. Deep Sea Res. 2 Top. Stud. Oceanogr. 56, 512–522. doi: 10.1016/j.dsr2.2008.12.005

CrossRef Full Text | Google Scholar

Pistek, P., and La Violette, P. E. (1999). Observations of the suppression of tide-generated nonlinear internal wave packets in the Strait of Gibraltar. J. Mar. Syst. 20, 113–128. doi: 10.1016/S0924-7963(98)00073-6

CrossRef Full Text | Google Scholar

Ramírez-Romero, E., Macías, D., García, C. M., and Bruno, M. (2014). Biogeochemical patterns in the Atlantic Inflow through the Strait of Gibraltar. Deep Sea Res. 1 Oceanogr. Res. Pap. 85, 88–100. doi: 10.1016/j.dsr.2013.12.004

CrossRef Full Text | Google Scholar

Richez, C. (1994). Airborne synthetic aperture radar tracking of internal waves in the Strait of Gibraltar. Prog. Oceanogr. 33, 93–159. doi: 10.1016/0079-6611(94)90023-X

CrossRef Full Text | Google Scholar

Richez, C., and Kergomard, C. (1990). “Characteristic features occurring in the Strait of Gibraltar as seen through remote sensing data,” in The Physical Oceanography of Sea Straits, ed. L. J. Pratt (Dordrecht: Springer), 441–455. doi: 10.1007/978-94-009-0677-8_21

CrossRef Full Text | Google Scholar

Sammartino, S., García-Lafuente, J., Naranjo, C., Sánchez-Garrido, J. C., Sánchez-Leal, R., and Sánchez-Román, A. (2015). Ten years of marine current measurements in Espartel Sill, Strait of Gibraltar. J. Geophys. Res. Ocean. 120, 6309–6328. doi: 10.1002/2014JC010674

CrossRef Full Text | Google Scholar

Sánchez-Garrido, J. C., García-Lafuente, J., Criado-Aldeanueva, F., Baquerizo, A., and Sannino, G. (2008). Time-spatial variability observed in velocity of propagation of the internal bore in the Strait of Gibraltar. J. Geophys. Res. Ocean. 113, 1–6. doi: 10.1029/2007JC004624

CrossRef Full Text | Google Scholar

Sánchez-Garrido, J. C., Sannino, G., Liberti, L., García-Lafuente, J., and Pratt, L. (2011). Numerical modeling of three-dimensional stratified tidal flow over Camarinal Sill, Strait of Gibraltar. J. Geophys. Res. Ocean. 116:12026. doi: 10.1029/2011JC007093

CrossRef Full Text | Google Scholar

Santana-Casiano, J. M., Gonzalez-Davila, M., and Laglera, L. M. (2002). The carbon dioxide system in the Strait of Gibraltar. Deep. Res. 2 Top. Stud. Oceanogr. 49, 4145–4161. doi: 10.1016/S0967-0645(02)00147-9

CrossRef Full Text | Google Scholar

Santana-Casiano, J. M., González-Dávila, M., Laglera-Baquer, L. M., and Rodriguez-Somoza, M. J. (2001). Carbon dioxide system in the Canary region during October 1995. Sci. Mar. 65, 41–49. doi: 10.3989/scimar.2001.65s141

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarmiento, J. L., and Gruber, N. (2006). Ocean Biogeochemical Dynamics. Geological Magazine. Princeton: Princeton University Press, 144:1034. doi: 10.1017/S0016756807003755

CrossRef Full Text | Google Scholar

Shadwick, E. H., Thomas, H., Azetsu-Scott, K., Greenan, B. J. W., Head, E., and Horne, E. (2011). Seasonal variability of dissolved inorganic carbon and surface water pCO2 in the Scotian Shelf region of the Northwestern Atlantic. Mar. Chem. 124, 23–37. doi: 10.1016/j.marchem.2010.11.004

CrossRef Full Text | Google Scholar

Shadwick, E. H., Thomas, H., Comeau, A., Craig, S. E., Hunt, C. W., and Salisbury, J. E. (2010). Air-sea CO2 fluxes on the Scotian Shelf: seasonal to multi-annual variability. Biogeosciences 7, 3851–3867. doi: 10.5194/bg-7-3851-2010

CrossRef Full Text | Google Scholar

Skliris, N., and Beckers, J. (2009). Modelling the Gibraltar Strait/Western Alboran Sea ecohydrodynamics. Dtsch. Hydrographische Z. 59, 489–508. doi: 10.1007/s10236-009-0185-6

CrossRef Full Text | Google Scholar

Stanichny, S., Tigny, V., Stanichnaya, R., and Djenidi, S. (2005). Wind driven upwelling along the African coast of the Strait of Gibraltar. Geophys. Res. Lett. 32, 1–4. doi: 10.1029/2004GL021760

CrossRef Full Text | Google Scholar

Takahashi, T., Sutherland, S. C., Sweeney, C., Poisson, A., Metzl, N., Tilbrook, B., et al. (2002). Global air-sea flux of CO₂ based on surface ocean pCO₂, and seasonal biological and temperature effects. Deep. Res. Part II 49, 1601– 1622.

Google Scholar

Takahashi, T., Olafsson, J., Goddard, J. G., Chipman, D. W., and Sutherland, S. C. (1993). Seasonal variation of CO2 and nutrients in the high-latitude surface oceans: a comparative study. Glob. Biogeochem. Cycles 7, 843–878. doi: 10.1029/93GB02263

CrossRef Full Text | Google Scholar

Vázquez-Escobar, A., Bruno, M., Izquierdo, A., Macías, D., and Ruiz-Cañavate, A. (2008). Meteorologically forced subinertial flows and internal wave generation at the main sill of the Strait of Gibraltar. Deep Sea Res. 2 Oceanogr. Res. Pap. 55, 1277–1283. doi: 10.1016/j.dsr.2008.05.008

CrossRef Full Text | Google Scholar

Walsh, J. J. (1991). Importance of continental margins in the marine biogeochemical cycling of carbon and nitrogen. Nature 350, 53–55. doi: 10.1038/350053a0

CrossRef Full Text | Google Scholar

Wanninkhof, R. (1992). Relationship between wind speed and gas exchange over the ocean. J. Geophys. Res. 97, 7373–7382. doi: 10.1029/92JC00188

CrossRef Full Text | Google Scholar

Wanninkhof, R. (2014). Relationship between wind speed and gas exchange over the ocean revisited. Limnol. Oceanogr. Methods 12, 351–362. doi: 10.4319/lom.2014.12.351

CrossRef Full Text | Google Scholar

Wanninkhof, R., and Feely, R. A. (1998). fCO2 dynamics in the Atlantic, South Pacific and South Indian oceans. Mar. Chem. 60, 15–31. doi: 10.1016/S0304-4203(97)00088-1

CrossRef Full Text | Google Scholar

Watson, G., and Robinson, S. (1990). A study of internal waves propagation in the Strait of Gibraltar using shore-based marine radar images. J. Phys. Oceanogr. 20, 374–395. doi: 10.1175/1520-0485(1990)020<0374:ASOIWP>2.0.CO;2

CrossRef Full Text | Google Scholar

Wesson, J. C., and Gregg, M. C. (1994). Mixing at camarinal sill in the Strait of Gibraltar. J. Geophys. Res. 99, 9847–9878. doi: 10.1029/94JC00256

CrossRef Full Text | Google Scholar