RESEARCH LETTER Seasonality of Dissolved Organic Carbon Exchange 10.1029/2020GL089601 Across the Key Points: X. A. Álvarez‐Salgado1 , J. Otero1 , S. Flecha2 , and I. E. Huertas2 • Dissolved organic carbon exchange across the Strait of Gibraltar exhibits 1CSIC Instituto de Investigacións Mariñas (IIM‐CSIC), Vigo, Spain, 2CSIC Instituto de Ciencias Marinas de Andalucía a marked seasonal variability −1 ‐ • Net DOC input of 4.2 Tg C yr (ICMAN CSIC), Puerto Real, Spain across the Strait of Gibraltar represents 53% of the external DOC input to the Abstract Exchange of mass, heat, and solutes across the Strait of Gibraltar is fundamental to understand • DOC of Atlantic origin supports 32% the circulation and biogeochemistry of the Mediterranean. Here we focus on the exchange of dissolved of the net heterotrophy of the Mediterranean Sea organic carbon (DOC) using data from 12 surveys conducted between 2008 and 2015. DOC exchange exhibits a marked bimodal distribution with minima in late June and late October and maxima in mid‐April Supporting Information: and late August. This pattern is mainly due to seasonal variation of the DOC gradient between the Atlantic • Supporting Information S1 Surface Water entering the Mediterranean and the deep opposite flow of Mediterranean Overflow Water. The gradient is controlled by the different seasonal cycles followed by the DOC in the two layers. Annual − average DOC import from the Atlantic is equivalent to 4.2 ± 1.5 Tg C yr 1, which represents 53% of the Correspondence to: X. A. Álvarez‐Salgado, external DOC inputs and contributes to support 32% of the net heterotrophy of the Mediterranean. [email protected] Plain Language Summary The Mediterranean Sea is a semienclosed basin connected with the through the Strait of Gibraltar. At this hot spot of ocean circulation, about 0.8 Sv Citation: ¼ 6 3 −1 ‐ (1 Sv 10 m s ) of dissolved organic carbon (DOC) rich Atlantic Surface Water enters the Mediterranean Álvarez Salgado, X. A., Otero, J., fl fl Flecha, S., & Huertas, I. E. (2020). Sea and the same volume of DOC poor Mediterranean Over ow Water ows oppositely to the Atlantic Seasonality of dissolved organic carbon Ocean. Both DOC concentrations and water flows are not stationary but vary seasonally. Differences in the exchange across the Strait of Gibraltar. amplitude and timing of those seasonal cycles produce a marked bimodal variation in the net DOC flux Geophysical Research Letters, 47, e2020GL089601. https://doi.org/ of Atlantic water that enters the Mediterranean Sea, with minima in late June and late October and maxima 10.1029/2020GL089601 in mid‐April and late August. This pattern has been observed for the first time and allowed us to better constrain this organic carbon flux, which represents about half of the total input of DOC in the Received 29 JUN 2020 Mediterranean Sea and supports about one third of its net organic carbon demand. Accepted 2 SEP 2020 Accepted article online 9 SEP 2020

1. Introduction The two‐layer water circulation in the Strait of Gibraltar (Figure 1) ties the overturning circulation cells of the North Atlantic and the Mediterranean Sea (Bryden & Kinder, 1991; Carracedo‐Segade et al., 2016; Reid, 1979; Volkov et al., 2019). Atlantic Surface Water (ASW) enters the Mediterranean basin across the strait and flows eastward until the Levantine basin, progressively gaining salinity due to evaporation. At the Rhodes gyre, it sinks and returns westward as Levantine Intermediate Water (LIW) in the Eastern basin and Eastern Intermediate Water (EIW) in the Western basin, constituting the shallow overturning circula- tion cell of the Mediterranean Sea (Schneider et al., 2014; Tsimplis et al., 2006). EIW mixes mostly with Western Mediterranean Deep Water (WMDW) to form the Mediterranean Overflow water (MOW) that leaves the Mediterranean Sea through the Strait of Gibraltar (Bryden & Stommel, 1982; García‐Lafuente et al., 2007; Naranjo et al., 2015). This conspicuous circulation pattern is the reason behind the extreme oli- gotrophy and net heterotrophy that dictate the carbon and nutrient biogeochemistry of the Mediterranean

(Luna et al., 2012; Reygondeau et al., 2017) and makes this basin a transient anthropogenic CO2 sink (Palmiéri et al., 2015). Given that the strait is a “hot spot” of this circulation cell, it has been repeatedly vis- ited to study its hydrography (García‐Lafuente et al., 2015; Millot, 2014; Naranjo et al., 2015) and dynamics (Bryden et al., 1994; García‐Lafuente et al., 2011; Naranjo et al., 2017; Sammartino et al., 2015). Similarly, the quantification of the exchange of ASW and MOW carrying contrasting biogeochemical properties has been ©2020. The Authors. This is an open access article under the crucial to constrain the net heterotrophy of the Mediterranean sea in terms of dissolved inorganic nutrients terms of the Creative Commons (Dafner et al., 2003; Huertas et al., 2012) and CO2 system components (Flecha et al., 2019; Huertas Attribution License, which permits use, et al., 2009). Furthermore, the MOW leaving the strait also impacts on the hydrography, circulation, and bio- distribution and reproduction in any medium, provided the original work is geochemistry of the adjacent Northeast Atlantic: It causes the massive subduction of Eastern North Atlantic properly cited. Central Water (ENACW) in the Gulf of Cadiz (GC) (Barbosa Aguiar et al., 2015; Baringer & Price, 1997;

ÁLVAREZ‐SALGADO ET AL. 1of9 Geophysical Research Letters 10.1029/2020GL089601

Figure 1. Map of the study area (a), topography and GIFT stations (b), and average salinity distribution in the water column (c). Position of ASW and MOW flows is also indicated along with the isohaline 37.5. GC, ES, and CS stand for Gulf of Cadiz, Espartel sill, and , respectively.

Rhein & Hinrichsen, 1993; Sánchez‐Leal et al., 2017) and likely affects the formation of North Atlantic Deep water (NADW) in the high latitudes of the North Atlantic, with direct implications for the Atlantic Meridional Overturning Circulation (Carracedo‐Segade et al., 2016). Particularly, the impact of the

subduction of ENACW: MOW has been quantified in terms of the downward transports of CO2 (Álvarez et al., 2005) and dissolved organic carbon (DOC) (Santana‐Falcón et al., 2017) to intermediate depths of

the Northeast Atlantic, where anthropogenic CO2 is stored and DOC is a major substrate for microbial growth. While the assessment of net inorganic carbon and nutrient fluxes between the North Atlantic and the Mediterranean Sea is based on time series of water transport records and repeated bottle measurements in ASW and MOW (Huertas et al., 2009, 2012), previous estimates of the net exchange of DOC through the Strait of Gibraltar are limited to discrete DOC measurements performed more than two decades ago (September 1997 and April 1998) and a constant value of water transport taken from the literature (Dafner, González‐Dávila, et al., 2001; Dafner, Sempéré, & Bryden, 2001). Latter, Santinelli (2015) also pro- vided an estimate based on Dafner's and her own data. These estimations have been considered as represen- tative of a stationary net import of DOC from the Atlantic toward the Mediterranean. This steady‐state assumption would not be accurate if water flows and DOC concentration profiles are not stationary throughout the year. In this regard, a marked seasonality in the thermohaline characteristics and water flow of the MOW at the Strait of Gibraltar have been reported in response to basin scale processes, such as the seasonal and interannual variability in the formation of WMDW in the Gulf of Lions (García‐ Lafuente et al., 2007; Naranjo et al., 2015, 2017; Sammartino et al., 2015; Schroeder et al., 2016). Furthermore, DOC profiles in temperate and subtropical areas are also characterized by well‐defined seaso- nal cycles (Carlson et al., 1994; Doval et al., 2016; Romera‐Castillo et al., 2013). Recently, Amaral et al. (2020) have also reported seasonal differences in the DOC of the GC. Therefore, there are solid evidences to hypothesize that the net DOC exchange fluxes along the Strait of Gibraltar are not stationary but exhibit a marked seasonal variability. To test the hypothesis we combined periodic measurements of DOC taken in the ASW and the MOW from 2008 to 2015 and water flows through the channel, to calculate net DOC

ÁLVAREZ‐SALGADO ET AL. 2of9 Geophysical Research Letters 10.1029/2020GL089601

exchange fluxes throughout the seasonal cycle. We further revisited the role of the Strait of Gibraltar in the DOC budget of the Mediterranean Sea on the light of this new estimate.

2. Materials and Methods 2.1. Sampling Program Data were acquired at the Gibraltar Fixed Time series (GIFT, Figure 1b) during 12 surveys carried out between 2008 and 2015 (Table S1 in the supporting information). Measurements were taken in the three sta- tions that form the GIFT, G1 (5°58.60′W, 35°51.68′N), G2 (5°44.75′W, 35°54.71′N) and G3 (5°22.10′W, 35°59.19′N) by following an identical procedure. In all surveys, a temperature and salinity profile was obtained with a Seabird 911Plus conductivity/temperature depth probe connected to a rosette sampler con- taining 10 L Niskin PVC bottles. Conductivity measurements were converted into practical values of the sali- nity scale with the UNESCO equation (1986). The accuracy of CTD measurements for temperature and salinity were 0.004°C and 0.005, respectively. Seawater samples for dissolved oxygen, nutrient salts, and chlorophyll measurements were subsequently collected with the rosette. Sampling occurred at variable depths (from five to eight levels) depending on the instant position of the interface between the ASW and the MOW in the water column (Figure 1c), which had been previously identified by the CTD profiles. On average, 21 (range 16 to 26) water samples were collected on each survey. Dissolved oxygen concentration was determined by automated potentiometric modification of the original Winkler method using a Titroprocessor (model Metrohm 794). Upon collection, Winkler flasks were sealed, stored in darkness, − and measured within 24 hr. The error of measurements was ±2 μmol kg 1 (n ¼ 252). Apparent oxygen uti- lization (AOU) values were calculated with the solubility equation (Benson & Krause, 1984). For nutrients (nitrate, phosphate, and silicate) determination, water samples (5 ml, two replicates) were taken from the Niskin bottles in 50 ml polyethylene bottles, filtered immediately through Whatman GF/F glass fibber filters (0.7 μm pore size), and stored frozen at −20°C for later analyses in the shore‐based laboratory. Nutrients con- centrations were measured with a continuous flow autoanalyzer using standard colorimetric techniques (Hansen & Koroleff, 1999). Analytical precisions were always better than ±3%. Chlorophyll analysis was conducted by filtering 0.5 L samples through Whatman GF/F filters, extracting in 90% acetone, and measur- ing concentration by standard fluorometric methods (Parsons et al., 1984) using a Turner Designs Model 10 fluorometer. The fluorometer was calibrated using pure chlorophyll a from the cyanobacterium Anacystis nidulans (Sigma chemical Co.) with the concentration determined spectrophotometrically.

2.2. DOC Measurements Seawater samples for the determination of DOC were collected in 0.25 L acid–cleaned glass bottles and immediately filtered through precombusted (450°C, 4 hr) Whatman GF/F filters with an acid‐cleaned all‐ glass filtration system previously rinsed with about 50 ml of the sample. Aliquots of 20 ml were collected

for DOC analysis in precombusted (450°C, 12 hr) 24 ml glass vials. After acidification with H3PO4 (85%) to pH < 2, they were sealed with Teflon‐lined caps and stored in the dark at 4°C until analyzed in the shore‐based laboratory within 1 month after collection. DOC content in samples was measured with a com- mercial Shimadzu TOC‐VCPH organic carbon analyzer working under the principle of high‐temperature catalytic oxidation according to the protocol described by Álvarez‐Salgado and Miller (1998). Potassium hydrogen phthalate (99.95–100.05%, p.a., Merck) was used to calibrate the system. The precision of the ana- − lyzer was ±0.5 μmol L 1. Accuracy was checked with consensus reference materials provided by D. A. Hansell (University of Miami, USA).

2.3. Seasonal Cycles Modeling Generalized additive models (GAMs) (Wood, 2006) were used to infer the seasonal cycles of DOC and com- panion variables in the two opposite flowing layers of the strait. The surface eastward flowing ASW was demarcated by salinity <37.5 and the bottom westward flowing MOW by salinity >37.5 (Figure 1c; Huertas et al., 2009, 2012). Data collected deeper than 500 m (n ¼ 7 cases) were excluded from the analyses. Independent GAMs were fitted to the ASW and MOW data. A GAM is a nonparametric regression technique that allows inspecting the relationship between a response variable and one (or more) explanatory variable(s) without the need to choose a particular parametric form for describing the shape of the relationship(s). In this work, we modeled the variability of DOC and other

ÁLVAREZ‐SALGADO ET AL. 3of9 Geophysical Research Letters 10.1029/2020GL089601

variables (Chl a, dissolved oxygen, and nutrients) using the day of the year (DoY), salinity (S) and potential temperature (T) as explanatory variables. DoY accounts for the seasonal cycle while salinity and temperature standardize the effect of thermohaline variability within the two layers. Salinity and temperature ranged from 37.483 to 35.876 and from 13.208 to 23.446°C in ASW, respectively. In MOW they varied from 37.504 to 38.513 and from 12.965 to 14.632°C. GAMs were formulated as follows: ÀÁÀÁ ÀÁ Y i; j ¼ α þ g DoY i; j þ h1 Si; j þ h2 Ti; j þ ϵi; j

where Y is the DOC, or any other variable of interest, measured at a day i and depth level j, α is the inter-

cept, and εi,j is the error term assumed to be normally distributed; g and hs are nonparametric smoothing functions specifying the effect of the covariates (DoY, S, and T) on the response variable Y. Smoothing functions g and h were fit by penalized cyclic cubic and cubic regression splines, respectively, and the number of knots were restricted to a maximum of six and three knots for each smoothing function type. Note that when modeling S and T we included in the model the covariate depth (in m).

3. Results and Discussion 3.1. Seasonal Variability of DOC at the Strait of Gibraltar Average DOC profiles along the Strait of Gibraltar (Figure S1a) and year round (Figure S1b) show well‐defined spatial and temporal patterns. The DOC‐rich ASW is characterized by higher concentrations at the Atlantic side of the strait (stn G1). Conversely, the DOC‐poor MOW presents lower concentration in the bottom layer of the Mediterranean side of the strait (stn G3) (Figure S1a). The annual cycle of the aver- age DOC profile (Figure S1b) shows a differential seasonal accumulation of DOC in both the ASW and MOW. Individual GAMs were tested for the three stations along the strait (Figure 1b) finding roughly equivalent patterns between stations (Figure S2 and Table S2). Therefore, GAMs were fit to data for the three stations together to obtain average seasonal distributions of DOC along the strait, which is more convenient to calculate net exchange fluxes rather than taking individual stations at both ends of the strait. These statistical models described well the seasonal cycle of DOC (Table S3) and the rest of variables in both ASW (Table S4) and MOW (Table S5) flowing layers. Regarding the DOC in ASW, although the model − explained 36.4% of the variance, the standard error of the estimate was low (about 1 μmol L 1) and the sig- nificance of DoY as explanatory variable was very high (p < 0.001) (Table S3). The time course of DOC in − ASW was characterized by a well‐defined annual cycle with minimum levels of 64 μmol L 1 by mid‐June − and maximum of 88 μmol L 1 by late August (Figure 2a). In the MOW the model produced comparable results in terms of the error of the estimate and significance of DoY and showed also a clear annual cycle although with a lower seasonal amplitude and decoupled from the ASW as the seasonal minimum of − − 49 μmol L 1 happens by early May and the maximum of 66 μmol L 1 by early August (Figure 2a). Thus, the seasonal DOC extremes occurred 20–30 days later in the ASW than in the MOW. The contrasting origin of the water masses that meet in the Strait of Gibraltar and their opposed circulation are the reasons behind the observed differences in the amplitude and timing of the seasonal cycles of DOC. ASW flowing eastward in the upper layer and MOW flowing westward in depth constitute the shallow open overturning circulation cell of the Mediterranean Sea (Schneider et al., 2014; Tsimplis et al., 2006). Thus, ASW is a well‐ventilated DOC‐rich water mass originated in the Eastern North Atlantic Ocean and trans- ported to the Strait of Gibraltar by the Portugal and Azores currents (García‐Lafuente et al., 2015; Sánchez‐Leal & Relvas, 2003; Santana‐Falcón et al., 2017). ASW also modifies further while crossing the GC (Amaral et al., 2020; Criado‐Aldeanueva et al., 2006) (Figure 1a), a basin characterized by the classical seasonal cycle of primary production of temperate regions (Navarro & Ruiz, 2006). Furthermore, a consider- able discharge of DOC from the Guadalquivir river estuary to the continental shelf of the GC has been docu- mented (Amaral et al., 2020; de la Paz et al., 2007; Flecha et al., 2015). Therefore, the ASW likely receives regional inputs of organic matter with a marked temporal pattern during transit toward the Strait. On the other hand, the MOW is largely a mixture of EIW and WMDW. EIW is an aged water mass ultimately − formed in the Levantine basin and it is characterized by its high salinity (38.84) and low DOC (51 μmol L 1) (Catalá et al., 2018; Santinelli, 2015; Santinelli et al., 2010). WMDW forms by deep convection in the Gulf of

ÁLVAREZ‐SALGADO ET AL. 4of9 Geophysical Research Letters 10.1029/2020GL089601

Lions. Several varieties of this water mass, previous and posterior to the Western Mediterranean Transient (WMT), coexist in the Western Mediterranean basin, with salinities ranging from 38.46 to 38.49 and − DOC 42–43 μmol L 1 (Catalá et al., 2018). A time series study of the seasonal variability of DOC in surface ocean waters of the temperate Eastern North Atlantic has not been reported yet. The nearest available time series observations are in the oligotrophic waters of the Bermuda Atlantic Time Series station (BATS), in the western side of the North Atlantic subtropical gyre. DOC at BATS displays a marked seasonal cycle characterized by a minimum in February‐March and a maximum in July‐August (Carlson et al., 1994). Thus, the DOC minimum overlaps with the winter convective mixing period due to DOC dilution within the winter mixed layer. The DOC maximum hap- pens during summer stratification due to “malfunctioning” of the micro- bial food web caused by nutrient limitation in warm, highly stratified waters (Thingstad et al., 1997). By contrast, the DOC minimum in ASW of the Strait of Gibraltar does not occur during winter mixing: the lowest temperature in ASW is observed in early February (Figure S3), that is, about 3 months before the DOC minimum. Conversely, the DOC maxi- mum overlaps with the highest value of temperature and lowest concen- trations of Chl a and nitrate (Figure S3), that is, the period of maximum stratification. In this sense, the seasonal cycle of primary production in the GC (Navarro & Ruiz, 2006) would be coupled to the fluctuations in DOC concentrations measured in the ASW in the same way as already described in BATS. Interestingly, a quasi‐permanent phytoplankton pro- liferation is found around the Trafalgar Cape (NW of the Strait), which has been attributed to local fertilization by the injection of nutrients from the MOW due to tidal effect, but which is particularly pronounced at the end of summer (Macías et al., 2007; Navarro & Ruiz, 2006; Sala et al., 2018). Macías et al. (2007) also showed an eastward advection of phytoplankton toward the Alboran Sea by the Atlantic jet, which could well carry the DOC formed by the local bloom occurring at the end of summer, contri- buting to the DOC maximum detected over this period of the year. In the case of the MOW, seasonal minimum levels of DOC at the Strait of Gibraltar are coincident with the DOC concentration in the EIW (Catalá et al., 2018). The remarkable amplitude of the seasonal change of DOC − in the MOW, 17 μmol L 1 (Figure 2a), suggests a vertical transference from the ASW by the strong shear and tidal induced mixing between both water masses while crossing the Strait in opposite directions, particularly at the Camarinal sill (Figures 1b and 1c) area (Dafner, González‐Dávila, et al., 2001; García‐Lafuente et al., 2007; Macías et al., 2007). Even though salinity, potential temperature, AOU, and inorganic nutrients in the MOW also experience seasonal changes (Figure S4), they are small com- pared with the temporal changes in DOC. In this regard, Santinelli et al. (2012) found a highly significant linear relationship between DOC − Figure 2. Seasonal cycles of DOC (in μmol L 1) in ASW (red) and MOW and AOU in the core of the LIW, with DOC representing 49% of the oxy- − (blue) (a), ΔDOC (in μmol L 1) (b), water transport in the MOW (in gen demand of this water mass in the Eastern Mediterranean. However, ¼ 6 3 −1 fl 9 −1 Sv 10 m s ) (c), and DOC exchange ux (in 10 g C day ) (d). Note the DOC‐AOU relationship was lost in the Western Mediterranean, as that panels (a), (b), and (d) show the mean value and the standard error, we observed in the Strait of Gibraltar. whereas panel (c) shows the daily average (orange curve) and 1σ interval (gray area), and superimposed in black a GAM fit. Vertical bars in panel (d) show the range of values provided by Dafner, González‐Dávila, 3.2. DOC Exchange Across the Strait of Gibraltar et al. (2001) for September 1997 and Dafner, Sempéré, and Bryden (2001) The decoupling between the annual cycles of DOC in the upper Atlantic for April 1998. layer and Mediterranean bottom waters in the Strait of Gibraltar results

ÁLVAREZ‐SALGADO ET AL. 5of9 Geophysical Research Letters 10.1029/2020GL089601

in a marked bimodal distribution of the annual cycle of the difference between the DOC in the ASW and the − MOW (ΔDOC; Figure 2b). Two ΔDOC maxima of 18 and 22 μmol L 1 occur by mid‐April and end of August − whereas two ΔDOC minima of 6 and 8 μmol L 1 are observed by the beginning of summer and fall at these latitudes. When ΔDOC is multiplied by the value of the Mediterranean water transport (Q) measured in the Atlantic side of the Strait (at the Espartel sill area, Figure 1b) the net exchange of DOC across the channel can be obtained (F ¼ ΔDOC × Q). Here we have used the modeled seasonal variability of the DOC gradient in the Strait (Figure 2b) and the seasonal pattern observed in the Mediterranean water transport (Figure 2c) provided by Sammartino et al. (2015). This procedure allows to calculate the seasonal cycle of the net exchange of DOC between the Atlantic and the Mediterranean (Figure 2d). As expected, the annual distribution of DOC fluxes is also bimodal, although the spring and summer max- − ima of F are of about the same magnitude (15109 g C day 1) compared with the two ΔDOC maxima, 20% higher in summer than in spring. This can be explained by the seasonal pattern of Q, which is 20% faster in spring (0.90 Sv) than in summer (0.75 Sv) (Sammartino et al., 2015) (Figure 2c). DOC exchange at the − Strait of Gibraltar varies over a wide range of values, from 5 to 15109 g C day 1 and our work suggests that the periods of minimum and maximum DOC import from the Atlantic toward the Mediterranean can be anticipated. Previous estimates of net DOC exchange through the Strait are restricted to two values obtained with dis- crete measurements collected in September 1997 (Dafner, González‐Dávila, et al., 2001) and April 1998 (Dafner, Sempéré, & Bryden, 2001). The small differences between both estimates, 0.28–0.35 and 0.30– − − 0.561012 mol C yr 1 (or 10–18 and 9–12 109 g C day 1) would not lead to infer the occurrence of a seasonal cycle of net DOC transport. However, when the estimates by Dafner and coworkers are superimposed on the seasonal cycle of the net DOC exchange provided here, they fit conveniently (Figure 2d). These authors used data taken in spring and late summer, when the fluxes between both basins are relatively high. Hence, con- sidering the large seasonal variability of DOC fluxes, any estimate based in observations gathered at a short − − time period should be expressed on daily (day 1) or monthly basis (month 1) whereas annual extrapolations should be avoided. 3.3. Contribution of Water Exchange Through the Strait of Gibraltar to the DOC Budget of the Mediterranean Sea The annual average net transport of DOC from the North Atlantic toward the Mediterranean Sea through − the Strait of Gibraltar is 4.2 ± 1.51012 gCyr 1, as calculated from the seasonal cycle of DOC exchange (Figure 2d). Note that our flux estimates cover from DoY 58 to 338, the period comprised by the surveys (Table S1). Therefore, our annual average net DOC transport does not take into account the fluxes from December to February. Extrapolation of the modeled DOC gradient (Figure 2b) to the winter months does − − not alter significantly the annual average net flux. Our estimate is half than the 7.7–9.710 12 gCyr 1 reported by Santinelli (2015) because we have used more recent and accurate water flows. In order to con- textualize the net DOC transport from the Atlantic with other DOC sources to the Mediterranean, we − obtained that the atmospheric DOC inputs represent 0.45–1.801012 gCyr1 (Copin‐Montegut, 1993; Djaoudi et al., 2018; Galletti et al., 2020; Loÿe‐Pilot et al., 1992; Violaki et al., 2018; Willey et al., 2000) and − riverine inputs account for by 1.0–1.6 1012 gCyr 1. These fluvial fluxes were calculated from the average − − − DOC concentrations in rivers Po (215 μmol C L 1) Rhone (250 μmol C L 1) and Nile (301 μmol C L 1) recently reviewed by Dai et al. (2012) and using the total freshwater flow to the Mediterranean Sea provided − by Lionello et al. (2012) (12.7–14.3; Violaki et al., 2018103 m3 s 1). Our estimates are lower than the atmo- − − spheric input of 0.4–13 1012 gCyr 1 and slightly higher than the fluvial input of 0.6–0.71012 gCyr 1 pro- vided by Santinelli (2015). The discrepancy in the range of atmospheric DOC inputs is because we have discarded outliers when extrapolating to the entire Mediterranean. In fact, Santinelli et al. (2015) state that − 1.3–4.6 g C yr 1 would be a more reasonable range. Finally, the net DOC input from the Black Sea across the − Dardanelles strait has been estimated in 1.21–1.261012 gCyr1 (Copin‐Montegut, 1993; Polat & Tugrul, 1996). Therefore, according to our assessment, the Atlantic Ocean can be considered the main source − of DOC to the Mediterranean Sea, representing 53% of the total DOC input of 7.9 ± 1.61012 gCyr 1. This percentage is very consistent with the 55% contribution of the ASW inflow to the DOC of the deep Mediterranean Sea based on δ13C‐DOC measurements by Santinelli et al. (2015) despite our disagreements in the estimates of next fluxes.

ÁLVAREZ‐SALGADO ET AL. 6of9 Geophysical Research Letters 10.1029/2020GL089601

Applying a similar two‐layer model of water mass exchange, Huertas et al. (2009) estimated the net export of

dissolved inorganic carbon (CT) from the Mediterranean Sea to the adjacent Atlantic Ocean as 12 −1 29.2 ± 0.610 gCyr , after discounting the anthropogenic CT signal. These authors also reported that 12 −1 45% of such export (13.2 ± 0.210 gCyr ) corresponded to CT produced during organic matter minerali- zation. Therefore, total DOC inputs to the Mediterranean Sea constitute 60 ± 12% of the organic carbon mineralization in this semienclosed basin, with 32 ± 12% of an Atlantic origin. The DOC transported east- ward in the ASW has an immediate impact on the heterotrophy of the adjacent Alboran Sea (Sempéré et al., 2003). The remaining 40% should come from the mineralization of organic particles suspended in flu- vial discharge, Dardanelle and Gibraltar straits flows, as well as that deposited from the atmosphere. Concerning fluvial particulate organic carbon (POC), the annual flux from river Rhone is about twice the flux of DOC (Sempéré et al., 2000). If this ratio is applicable to all Mediterranean rivers, the total fluvial flux − of POC would be 2.0–3.21012 gCyr 1, which would explain about 50% of the POC demand. The net POC − exchange across the Dardanelle strait has been estimated in 0.17–0.18 1012 gCyr1 (Polat & Tugrul, 1996). For the case of the Strait of Gibraltar, Santana‐Falcón et al. (2017) estimated a net POC − exchange of 0.6 1012 gCyr 1 from a single POC profile obtained in May 2014. Therefore, POC imported from the surface NE Atlantic Ocean and Marmara Sea would represent about 15% of the POC demand. The remaining 35% could be explained by the organic content of particles deposited from the atmosphere, likely during the frequent and massive Sahara dust events (Gallisai et al., 2016). Submarine groundwater discharge has also been suggested as a source of DOC to the Mediterranean but it would be essentially refractory (Santinelli et al., 2015). All these estimates are just based on the product of DOC or POC concentrations in marine, atmospheric, and fluvial waters to the Mediterranean multiplied by their respective water flows. Therefore, what we obtain is a net budget of inputs minus outputs, assuming steady state, and treating the Mediterranean as a “black box.” Despite all these limitations, our numbers reveal the major role played by DOC in the carbon cycle of the Mediterranean Sea and particularly by the DOC transported from the Atlantic Ocean by the ASW through the Strait of Gibraltar.

4. Conclusions Straits connecting ocean basins are “hot spots” of ocean circulation and biogeochemical cycles. Despite their importance, studies on the seasonal variability of the net exchange fluxes of essential compounds across straits are still scarce. Here we have shown that differences in the amplitude and timing of the seasonal cycles of DOC concentrations in ASW and MOW as well as water exchange fluxes across the Strait of Gibraltar produced a marked bimodal distribution of the net DOC flux to the Mediterranean basin. These fluxes are about threefold during the seasonal maxima that during the seasonal minima. Quantification of − this seasonality is crucial to better constraining this significant flux of 4.2 Tg C yr 1, which represents 53% of the external DOC input to the Mediterranean Sea and supports 32% of the net heterotrophy of this basin. Our estimate is based on discrete DOC measurements along the strait occupied on different months through several years. Using our data to validate a 3‐D hydrodynamic + biogeochemical model would provide better‐constrained results. Furthermore, application of DOC isotope, optical, or molecular proxies would allow gaining new insights on the role of the DOC transported by the Atlantic jet into the Mediterranean.

Acknowledgments We are grateful to the crews of all the Data Availability Statement research vessels involved in this study. We also thank Maria Ferrer, Antonio Data used in this study correspond to the database generated at the GIFT (Gibraltar Fixed Time Series) Moreno, and David Roque for their observatory and are available at Digital CSIC repository (https://digital.csic.es/handle/10261/205367). assistance with sample collection and analysis. Funding for this work was provided by the European Commission through the projects CARBOOCEAN References (FP6‐511176), SESAME (FP6‐036949), Álvarez, M., Pérez, F. F., Shoosmith, D. R., & Bryden, H. L. (2005). Unaccounted role of Mediterranean water in the drawdown of CARBOCHANGE (FP7‐264879), anthropogenic carbon. Journal of Geophysical Research, 110, C09S03. https://doi.org/10.1029/2004JC002633 PERSEUS (FP7‐287600), and Álvarez‐Salgado, X. A., & Miller, A. E. J. (1998). Simultaneous determination of dissolved organic carbon and total dissolved nitrogen in COMFORT (H2020‐820989). The work seawater by high temperature catalytic oxidation: Conditions for accurate shipboard measurements. Marine Chemistry, 62(3–4), reflects only the author's/authors' view; 325–333. https://doi.org/10.1016/S0304-4203(98)00037-1 the European Commission and their Amaral, V., Romera‐Castillo, C., & Forja, J. (2020). Dissolved organic matter in the Gulf of Cádiz: Distribution and drivers of chromophoric executive agency are not responsible for and fluorescent properties. Frontiers in Marine Sciences, 7, 126. https://doi.org/10.3389/fmars.2020.00126 any use that may be made of the Barbosa Aguiar, A. C., Peliz, A., Neves, F., Bashmachnikov, I., & Carton, X. (2015). Mediterranean outflow transports and entrainment information the work contains. estimates from observations and high‐resolution modelling. Progress in Oceanography, 131,33–45.

ÁLVAREZ‐SALGADO ET AL. 7of9 Geophysical Research Letters 10.1029/2020GL089601

Baringer, M. O., & Price, J. F. (1997). Mixing and spreading of the Mediterranean outflow. Journal of Physical Oceanography, 27(8), 1654–1677. https://doi.org/10.1175/1520-0485(1997)027<1654:MASOTM>2.0.CO;2 Benson, B. B., & Krause, D. (1984). The concentration and isotopic fractionation of oxygen dissolved in freshwater and seawater in equi- librium with the atmosphere. Limnology and Oceanography, 29(3), 620–632. https://doi.org/10.4319/lo.1984.29.3.0620 Bryden, H. L., Candela, J., & Kinder, T. H. (1994). Exchange through the Strait of Gibraltar. Progress in Oceanography, 33(3), 201–248. https://doi.org/10.1016/0079-6611(94)90028-0 Bryden, H. L., & Kinder, T. H. (1991). Steady two‐layer exchange through the Strait of Gibraltar. Deep‐Sea Research, 38, S445–S463. https:// doi.org/10.1016/S0198-0149(12)80020-3 Bryden, H. L., & Stommel, H. M. (1982). Origin of the Mediterranean outflow. Journal of Marine Research, 40,55–71. Carlson, C. A., Ducklow, H. W., & Michaels, A. F. (1994). Annual flux of dissolved organic carbon from the euphotic zone in the north- western Sargasso Sea. Nature, 371(6496), 405–408. https://doi.org/10.1038/371405a0 Carracedo‐Segade, L. I., Pardo, P. C., Flecha, S., & Pérez, F. F. (2016). On the Mediterranean water composition. Journal of Physical Oceanography, 46(4), 1339–1358. https://doi.org/10.1175/JPO-D-15-0095.1 Catalá, T. S., Martínez‐Pérez, A. M., Nieto‐Cid, M., Álvarez, M., Otero, J., Emelianov, M., et al. (2018). Dissolved organic matter (DOM) in the open Mediterranean Sea. I. Basin‐wide distribution and drivers of chromophoric DOM. Progress in Oceanography, 165,35–51. https://doi.org/10.1016/j.pocean.2018.05.002 Copin‐Montegut, C. (1993). Alkalinity and carbon budgets in the Mediterranean Sea. Global Biogeochemical Cycles, 7(4), 915–925. https:// doi.org/10.1029/93GB01826 Criado‐Aldeanueva, F., Garcia‐Lafuente, J., Vargas, J. M., Del Rio, J., Vazquez, A., Reul, A., & Sanchez, A. (2006). Distribution and cir- culation of water masses in the Gulf of Cadiz from in situ observations. Deep‐Sea Research II, 53(11–13), 1144–1160. https://doi.org/ 10.1016/j.dsr2.2006.04.012 Dafner, E., González‐Dávila, M., Santana‐Casiano, J. M., & Sempéré, R. (2001). Total organic and inorganic carbon exchange through the strait of Gibraltar in September 1997. Deep‐Sea Research I, 48(5), 1217–1235. https://doi.org/10.1016/S0967-0637(00)00064-9 Dafner, E. V., Boscolo, R., & Bryden, H. L. (2003). The N:Si:P molar ratio in the Strait of Gibraltar. Geophysical Research Letters, 30(10), 1506. https://doi.org/10.1029/2002GL016274 Dafner, E. V., Sempéré, R., & Bryden, H. L. (2001). Total organic carbon distribution and budget through the strait of Gibraltar in April 1998. Marine Chemistry, 73(3‐4), 233–252. https://doi.org/10.1016/S0304‐4203(00)00109‐2 Dai, M., Yin, Z., Meng, F., Liu, Q., & Cai, W.‐J. (2012). Spatial distribution of riverine DOC inputs to the ocean: An updated global synthesis. Current Opinion in Environmental Sustainability, 4(2), 170–178. https://doi.org/10.1016/j.cosust.2012.03.003 de la Paz, M., Gómez‐Parra, A., & Forja, J. (2007). Inorganic carbon dynamic and air–water CO2 exchange in the Guadalquivir Estuary (SW Iberian Peninsula). Journal of Marine Systems., 68(1–2), 265–277. https://doi.org/10.1016/j.jmarsys.2006.11.011 Djaoudi, K., Van Wambeke, F., Barani, A., Hélias‐Nunige, S., Sempéré, R., & Pulido‐Villena, E. (2018). Atmospheric fluxes of soluble organic C, N, and P to the Mediterranean Sea: Potential biogeochemical implications in the surface layer. Progress in Oceanography, 163, 59–69. https://doi.org/10.1016/j.pocean.2017.07.008 Doval, M. D., López, A., & Madriñán, M. (2016). A decade of weekly dissolved organic carbon values in the coastal upwelling of the NW Spain (Atlantic Galician rías). Marine Chemistry, 179,34–43. https://doi.org/10.1016/j.marchem.2016.01.009 Flecha, S., Huertas, I. E., Navarro, G., Morris, E. P., & Ruiz, J. (2015). Air–water CO2 fluxes in a highly heterotrophic estuary. Estuaries and Coasts, 38(6), 2295–2309. https://doi.org/10.1007/s12237-014-9923-1 Flecha, S., Pérez, F. F., Murata, A., Makaoui, A., & Huertas, I. E. (2019). Decadal acidification in Atlantic and Mediterranean water masses exchanging at the strait of Gibraltar. Scientific Reports, 9(1), 15533. https://doi.org/10.1038/s41598-019-52084-x Galletti, Y., Becagli, S., di Sarra, A., Gonnelli, M., Pulido‐Villena, E., Sferlazzo, D. M., et al. (2020). Atmospheric deposition of organic matter at a remote site in the Central Mediterranean Sea: Implications for marine ecosystem. Biogeosciences Discussions, 17(13), 3669–3684. https://doi.org/10.5194/bg-2020-14 Gallisai, R., Volpe, G., & Peters, F. (2016). Large Saharan dust storms: Implications for chlorophyll dynamics in the Mediterranean Sea. Global Biogeochemical Cycles, 30, 1725–1737. https://doi.org/10.1002/2016GB005404 García‐Lafuente, J., Naranjo, C., Sánchez‐Leal, R., Sammartino, S., Bellanco, M. J., Sánchez‐Garrido, J. C., & Soto‐Navarro, J. (2015). On the origin of the seasonal and interannual T–S variability of the inflow through the strait of Gibraltar. Deep‐Sea Research I, 101,38–53. https://doi.org/10.1016/j.dsr.2015.03.005 García‐Lafuente, J., Sánchez Román, A., Naranjo, C., & Sánchez Garrido, J. C. (2011). The very first transformation of the Mediterranean outflow in the Strait of Gibraltar. Journal of Geophysical Research, 116, C07010. https://doi.org/10.1029/2011JC006967 García‐Lafuente, J., Sánchez‐Román, A., Díaz del Río, G., Sannino, G., & Sánchez‐Garrido, J. C. (2007). Recent observations of seasonal variability of the Mediterranean outflow in the strait of Gibraltar. Journal of Geophysical Research, 112, C10005. https://doi.org/10.1029/ 2006JC003992 Hansen, H. P., & Koroleff, F. (1999). Determination of nutrients. In K. Grasshof, K. Kremling, & M. Ehrhardt (Eds.), Methods of seawater analysis (pp. 159–228). Weilheim: Wiley‐Vch. Huertas, I. E., Ríos, A. F., García‐Lafuente, J., Makaoui, A., Rodríguez‐Gálvez, S., Sánchez, A., et al. (2009). Anthropogenic and natural CO2 exchange through the strait of Gibraltar. Biogeosciences, 6(4), 647–662. https://doi.org/10.5194/bg-6-647-2009 Huertas, I. E., Ríos, A. F., García‐Lafuente, J., Navarro, G., Makaoui, A., Sánchez‐Román, A., et al. (2012). Atlantic forcing of the Mediterranean oligotrophy. Global Biochememical Cycles, 26, GB2022. https://doi.org/10.1029/2011GB004167 Lionello, P., Abrantes, F., Congedi, L., Dulac, F., Gacic, M., Gomis, D., et al. (2012). Introduction: Mediterranean climate—Background information. In P. Lionello (Ed.), The climate of the Mediterranean Region: From the past to the future (pp. 35–90). London: Elsevier. https://doi.org/10.1016/B978-0-12-416042-2.00012-4 Loÿe‐Pilot, M.‐D., Cauwet, G., Spitzy, A., & Martin, J.‐M. (1992). Preliminary results on atmospheric wet deposition of organic carbon and nitrogen in Corsica. In J.‐M. Martin, & H. Barth (Eds.), Water Pollution Research Report 28, EROS 2000, Third Workshop (pp. 519–532). Brussels: Commission of the European Communities. Luna, G. M., Bianchelli, S., Decembrini, F., De Domenico, E., Danovaro, R., & Dell'Anno, A. (2012). The dark portion of the Mediterranean Sea is a bioreactor of organic matter cycling. Global Biogeochemical Cycles, 26, GB2017. http://doi.org/10.1029/2011GB004168 Macías, D., Martin, A. P., García, C. M., Bruno, M., García‐Lafuente, J., Izquierdo, A., 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. Progress in Oceanography, 74(2–3), 252–272. https://doi.org/10.1016/j.pocean.2007.04.006 Millot, C. (2014). Heterogeneities of in‐ and out‐flows in the Mediterranean Sea. Progress in Oceanography, 120, 254–278. https://doi.org/ 10.1016/j.pocean.2013.09.007

ÁLVAREZ‐SALGADO ET AL. 8of9 Geophysical Research Letters 10.1029/2020GL089601

Naranjo, C., García‐Lafuente, J., Sammartino, S., Sánchez‐Garrido, J. C., Sánchez‐Leal, R., & Bellanco, J. M. (2017). Recent changes (2004– 2016) of temperature and salinity in the Mediterranean outflow. Geophysical Research Letters, 44, 5665–5672. https://doi.org/10.1002/ 2017GL072615 Naranjo, C., Sammartino, S., García‐Lafuente, J., Bellanco, M. J., & Taupier‐Letage, I. (2015). Mediterranean waters along and across the strait of Gibraltar, characterization and zonal modification. Deep Sea Research I, 105,41–52. https://doi.org/10.1016/j.dsr.2015.08.003 Navarro, G., & Ruiz, J. (2006). Spatial and temporal variability of phytoplankton in the Gulf of Cádiz through remote sensing images. Deep Sea Research II, 53(11–13), 1241–1260. https://doi.org/10.1016/j.dsr2.2006.04.014 Palmiéri, J., Orr, J. C., Dutay, J. C., Béranger, K., Schneider, A., Beuvier, J., & Somot, S. (2015). Simulated anthropogenic CO2 storage and acidification of the Mediterranean Sea. Biogeosciences, 12, 781–802. https://doi.org/10.5194/bg-12-781-2015 Parsons, T. R., Maita, Y., & Lalli, C. M. (1984). A manual of chemical and biological methods for seawater analysis. Oxford: Pergamon Press. Polat, C., & Tugrul, S. (1996). Chemical exchange between the Mediterranean and the Black Sea via the Turkish straits. Bulletin de l'institut Oceanographique, Monaco, 17, 167–186. Reid, J. L. (1979). On the contribution of the Mediterranean Sea outflow to the Norwegian‐Greenland Sea. Deep‐Sea Research A, 26(11), 1199–1223. https://doi.org/10.1016/0198-0149(79)90064-5 Reygondeau, G., Guieu, C., Benedetti, F., Irisson, J. O., Ayata, S. D., Gasparini, S., & Koubbi, P. (2017). Biogeochemical regions of the Mediterranean Sea: An objective multidimensional and multivariate environmental approach. Progress in Oceanography, 151, 138–148. https://doi.org/10.1016/J.POCEAN.2016.11.001 Rhein, M., & Hinrichsen, H. H. (1993). Modification of Mediterranean water in the Gulf of Cadiz, studied with hydrographic, nutrient and chlorofluoromethane data. Deep Sea Research I, 40(2), 267–291. https://doi.org/10.1016/0967-0637(93)90004-M Romera‐Castillo, C., Álvarez‐Salgado, X. A., Galí, M., Marrasé, C., & Gasol, J. M. (2013). Combined effect of light exposure and microbial activity on distinct dissolved organic matter pools. A seasonal field study in an oligotrophic coastal system (Blanes Bay, NW Mediterranean). Marine Chemistry, 148,44–51. https://doi.org/10.1016/j.marchem.2012.10.004 Sala, I., Navarro, G., Bolado‐Penagos, M., Echevarría, F., & García, C. M. (2018). High‐chlorophyll‐area assessment based on remote sen- sing observations: The case study of Cape Trafalgar. Remote Sensing, 10(2), 165. https://doi.org/10.3390/rs10020165 Sammartino, S., García‐Lafuente, J., Naranjo, C., Sánchez‐Garrido, J. C., Sánchez‐Leal, R., & Sánchez‐Román, A. (2015). Ten years of marine current measurements in Espartel Sill, Strait of Gibraltar. Journal of Geophysical Research: Oceans, 120, 6309–6328. https://doi. org/10.1002/2014JC010674 Sánchez‐Leal, R. F., Bellanco, M. J., Fernández‐Salas, L. M., García‐Lafuente, J., Gasser‐Rubinat, M., González‐Pola, C., et al. (2017). The Mediterranean overflow in the Gulf of Cadiz: A rugged journey. Science Advances, 3(11), eaao0609. https://doi.org/10.1126/sciadv. aao0609 Sánchez‐Leal, R. F., & Relvas, P. (2003). Spring–summer climatological circulation in the upper layer in the region of Cape St. Vincent, Southwest Portugal. ICES Journal of Marine Science, 60(6), 1232–1250. https://doi.org/10.1016/S1054-3139(03)00137-1 Santana‐Falcón, Y., Álvarez‐Salgado, X. A., Pérez‐Hernández, M. D., Hernández‐Guerra, A., Mason, E., & Arístegui, J. (2017). Organic carbon budget for the eastern boundary of the North Atlantic subtropical gyre: Major role of DOC in mesopelagic respiration. Scientific Reports, 7(1), 10129. https://doi.org/10.1038/s41598-017-10974-y Santinelli, C. (2015). DOC in the Mediterranean Sea. In D. A. Hansell, & C. A. Carlson (Eds.), Biogeochemistry of marine dissolved organic matter (second ed.pp. 579–608). Boston, MA: Academic Press. https://doi.org/10.1016/B978-0-12-405940-5.00013-3 Santinelli, C., Follet, C., Retelletti Brogi, S., Xu, L., & Repeta, D. (2015). Carbon isotope measurements reveal unexpected cycling of dis- solved organic matter in the deep Mediterranean Sea. Marine Chemistry, 177, 267–277. https://doi.org/10.1016/j.marchem.2015.06.018 Santinelli, C., Nannicini, L., & Seritti, A. (2010). DOC dynamics in the meso and bathypelagic layers of the Mediterranean Sea. Deep‐Sea Research II, 57(16), 1446–1459. https://doi.org/10.1016/j.dsr2.2010.02.014 Santinelli, C., Sempéré, R., Van Wambeke, F., Charriere, B., & Seritti, A. (2012). Organic carbon dynamics in the Mediterranean Sea: An integrated study. Global Biogeochemical Cycles, 26, GB4004. https://doi.org/10.1029/2011GB004151 Schneider, A., Tanhua, T., Roether, W., & Steinfeldt, R. (2014). Changes in ventilation of the Mediterranean Sea during the past 25 year. Ocean Sciences, 10(1), 1–16. https://doi.org/10.5194/os-10-1-2014 Schroeder, K., Chiggiato, J., Bryden, H. L., Borghini, M., & Ismail, S. B. (2016). Abrupt climate shift in the Western Mediterranean Sea. Scientific Reports, 6(1), 23009. https://doi.org/10.1038/srep23009 Sempéré, R., Charriere, B., Van Wambeke, F., & Cauwet, G. (2000). Carbon inputs of the Rhone River to the Mediterranean Sea: Biogeochemical implications. Global Biogeochemical Cycles, 14(2), 669–681. https://doi.org/10.1029/1999GB900069 Sempéré, R., Dafner, E., Van Wambeke, F., Lefèvre, D., Magen, C., Allègre, S., et al. (2003). Distribution and cycling of total organic carbon across the Almeria Oran Front Biogeochemical implications. Journal of Geophysical Research, 108(C11), 3361. https://doi.org/10.1029/ 2002JC001475 Thingstad, T. F., Hagström, Å., & Rassoulzadegan, F. (1997). Accumulation of degradable DOC in surface waters: Is it caused by a mal- functioning microbial loop? Limnology Oceanography, 42(2), 398–404. https://doi.org/10.4319/lo.1997.42.2.0398 Tsimplis, M. N., Zervakis, V., Josey, S. A., Peneva, E. L., Struglia, M. V., Stanev, E. V., et al. (2006). Changes in the oceanography of the Mediterranean Sea and their link to climate variability. Developments in Earth and Environmental Sciences, 4, 227–282. https://doi.org/ 10.1016/S1571-9197(06)80007-8 Violaki, K., Bourrin, F., Aubert, D., Kouvarakis, G., Delsaut, N., & Mihalopoulos, N. (2018). Organic phosphorus in atmospheric deposition over the Mediterranean Sea: An important missing piece of the phosphorus cycle. Progress in Oceanography, 163,50–58. https://doi.org/ 10.1016/j.pocean.2017.07.009 Volkov, D. L., Baringer, M., Smeed, D., Johns, W., & Landerer, F. W. (2019). Teleconnection between the Atlantic Meridional overturning circulation and sea level in the Mediterranean Sea. Journal of Climate, 32(3), 935–955. https://doi.org/10.1175/JCLI-D-18-0474.1 Willey, J. D., Kieber, R. J., Eyman, M. S., & Avery, G. B. (2000). Rainwater dissolved organic carbon: Concentrations and global flux. Global Biogeochemical Cycles, 14(1), 139–148. https://doi.org/10.1029/1999GB900036 Wood, S. (2006). Generalized additive models: An introduction. London: R. Chapman & Hall. https://doi.org/10.1201/9781420010404

ÁLVAREZ‐SALGADO ET AL. 9of9