Decadal acidification in Atlantic and Mediterranean water masses exchanging at the Strait of Gibraltar

Seawater pH is undergoing a decreasing trend due to the absorption of atmospheric CO2, a phenomenon known as ocean acidification (OA). Biogeochemical processes occurring naturally in the ocean also change pH and hence, for an accurate assessment of OA, the contribution of the natural component to the total pH variation must be quantified. In this work, we used 11 years (2005–2015) of biogeochemical measurements collected at the Strait of Gibraltar to estimate decadal trends of pH in two major Mediterranean water masses, the Western Mediterranean Deep Water (WMDW) and the Levantine Intermediate Water (LIW) and assess the magnitude of natural and anthropogenic components on the total pH change. The assessment was also performed in the North Atlantic Central Water (NACW) feeding the Mediterranean Sea. Our analysis revealed a significant human impact on all water masses in terms of accumulation of anthropogenic CO2. However, the decadal pH decline found in the WMDW and the NACW was markedly affected by natural processes, which accounted for by nearly 60% and 40% of the total pH decrease, respectively. The LIW did not exhibit a significant pH temporal trend although data indicated natural and anthropogenic perturbations on its biogeochemical signatures.

The atmospheric concentration of carbon dioxide (CO 2 ) has markedly increased from the beginning of the industrial era due to the release of carbon from combustion of fossil fuels, deforestation and large changes in land-use activities. Anthropogenic emissions in particular, became the dominant source of CO 2 to the atmosphere around the mid XX century and continue rising at present 1 . The global ocean plays a relevant role mitigating this rise through the absorption of around a quarter of the emissions 1 . Recently, the oceanic sink for anthropogenic CO 2 (C ANT ) over the period 1994-2007 has been estimated to represent 31% of the global emissions 2 . Hence, the oceanic withdrawal of atmospheric CO 2 alleviates the greenhouse effect, but there is also a downside, as the CO 2 absorbed by the ocean changes the chemistry of seawater by lowering its pH and the carbonate ion (CO 2 3− ) levels. The chemical reactions derived from the uptake of CO 2 are generally termed ocean acidification 3 and threaten the overall structure of marine ecosystems at a planetary scale 4 .
Ocean acidification trends have been already documented through multi-decadal measurements in open ocean marine observatories 5 . It has been also suggested that certain basins, such as marginal seas, will be more impacted by the phenomenon than others 6 . This is the case of the Mediterranean Sea (MedSea), whose high sensitivity to acidification has been attributed to particular biogeochemical features and water circulation patterns 7,8 . The high alkalinity of Mediterranean waters and an active overturning circulation result in both an increased absorption of atmospheric CO 2 and an intensified carbon transport from the surface to the ocean interior 6,9,10 . Moreover, C ANT is imported continuously from the North Atlantic through the Strait of Gibraltar (SoG, Fig. 1A) 10,11 , leading to its accumulation in the basin at a long term.
The MedSea is indeed experiencing a decrease in pH although different declining rates have been reported 9,10,[12][13][14] . Some works have provided a total pH decline between −0.055 and −0.156 pH units 9,12 whereas other studies reported values ranging from −0.005 to −0.06 pH 10,13 . These discrepancies are based on the simulation of the quantity of C ANT that is taken up by Mediterranean water masses. This amount cannot be measured directly because it cannot be distinguished from the much larger natural carbon background and it has to be 1 calculated through indirect approaches. The estimation of C ANT using methods that combine observable physical and biogeochemical variables or models result then in different acidification rates.
However, seawater pH also varies "naturally" as a consequence of physical and biological processes, such as ocean circulation mechanisms and/or changes in ecosystem metabolism [15][16][17] . In fact, the pH decline due to biology in Atlantic water masses is of the same order of magnitude than the pH decrease caused by C ANT absorption 15 . Similarly, in the Northeast Pacific coastal ocean, long-term ocean acidification is modulated by the pH variability driven by ecosystem metabolism 17 . Hence, the accurate assessment of the actual extent of ocean acidification requires the identification of the fraction of the total pH change that is purely due to natural processes and discriminate it from the anthropogenic fraction. This analysis is particularly relevant in the MedSea where biogeochemistry is closely linked to water circulation patterns 18,19 .
In the current study we examined the relative contribution of both C ANT uptake and biology on temporal pH trends in two major Mediterranean water masses, the Western Mediterranean Deep Water (WMDW) and the Levantine Intermediate Water (LIW). They are originated in distant regions, such as the Gulf of Lions located in the Western Mediterranean and the Levantine basin in the Eastern Mediterranean, respectively. Nevertheless, both water masses can be easily tracked in the SoG by their characteristic thermohaline signatures 13,20 .
More specifically, ocean circulation in the MedSea starts with the entry of surface water from the Atlantic Ocean through the SoG, which flows towards the east while gaining salinity and density. In the Levantine basin the density has reached such an extent that the water sinks to form the LIW, which flows back into the opposite direction in a few 100 m depth. In the Gulf of Lions, high salinity waters cool down regularly in winter reaching enhanced densities and forming the WMDW, with the highly saline LIW acting as a preconditioner for deep water formation 21,22 . Both water masses leave jointly the MedSea as part of a dense (salinity >38) Mediterranean Outflow Water (MOW) that occupies the bottom layer in the SoG (Fig. 1C). The MOW is compensated by the surface eastward warmer (and less salty) Atlantic Inflow (AI, Fig. 1C) that is mainly constituted by the North Atlantic Central Water (NACW). The MOW and AI are clearly distinguishable within the channel by their different physical and biogeochemical properties, such as potential temperature, salinity, alkalinity, dissolved oxygen, inorganic nutrients and trace gases, as shown by measurements at the Gibraltar Fixed Time series (GIFT, Fig. 1B) 11,23,24 .
Here, we present the first (to our knowledge) measurement-based temporal evolution of the carbon system parameters in water masses exchanging through the SoG using observations collected between the years 2005-2015, with special emphasis in decadal trends of pH. The anthropogenic (pH nobio ) and natural (pH bio ) components of the decadal pH change were subsequently obtained in order to identify the relative contribution of both drivers on the pH variability observed in each water mass.
Data acquisition procedure was identical for all campaigns. A temperature and salinity profile was obtained with a Seabird 911Plus CTD probe. Seawater was subsequently collected for biogeochemical analysis using Niskin bottles immersed in an oceanographic rosette platform at variable depths (from 5 to 8 levels) depending on the instant position of the interface between the AI and MOW that was identified by the CTD profiles ). pH T25 data were obtained by the spectrophotometric method with m-cresol purple as indicator 25 . Samples were taken directly from the oceanographic bottles in 10 cm path-length optical glass cells and measurements were carried out with a Shimadzu UV-2401PC spectrophotometer containing a 25 °C-thermostated cells holder 11 . Precision and accuracy of the pH measurements were ± 0.0052 and ± 0.0055 respectively, which were determined from measurements of certified reference material (CRM batches#94, #97 and #136 provided by Prof. Andrew Dickson, Scripps Institution of Oceanography, La Jolla, CA, USA).
For further analysis of the influence of the age of water masses on the pH evolution, the theoretical seawater pH values that would result exclusively from the uptake of atmospheric CO 2 were also calculated with the CO2SYS program 26 using temperature, pressure, A T , pH T25 and inorganic nutrients as inputs parameters. Atmospheric CO 2 levels were taken from the Izaña Observatory (Meteorological State Agency of Spain, AEMET, https://www.esrl.noaa.gov/) Samples for A T analysis were collected in 500-ml borosilicate bottles, and poisoned with 100 μl of HgCl 2 -saturated aqueous solution and stored until measurement in the laboratory. A T was measured by potential titration 27 with a Titroprocessor (model Metrohm 794). Precision and accuracy of A T measurements were ±0.9 and ±5.4 μmol kg −1 respectively, which were determined from measurements of 54 CRMs of the above mentioned batches.
Concentration of dissolved inorganic carbon c(DIC) and partial pressure of dissolved CO 2 (pCO 2 ) were calculated using CO 2 SYS with the dissociation constants for carbon 28,29 and sulphate 30 .
DO concentration [c(DO)] was determined through automated potentiometric modification of the original Winkler method using the Titroprocessor. Upon collection, flasks were sealed, stored in darkness and measured within 24 h. The error of measurements was ±2 μmol kg −1 (n = 115). Apparent Oxygen Utilization (AOU) values were calculated with the solubility equation 31 .
Water samples (5 mL, two replicates) for inorganic nutrients determination were taken, filtered immediately (Whatman GF/F, 0.7 μm) and stored frozen 32 for later analyses in the shore-based laboratory. Nutrients concentrations were measured with a continuous flow auto-analyzer (SkalarSan ++215) using standard colorimetric techniques 33 . Analytical precisions were always better than ±3%.
Anthropogenic carbon estimation. C ANT concentration was estimated using the back-calculation technique 34 by applying the following equation: where the stoichiometric coefficients R C (−ΔO 2 /ΔC) = 1.45 and R N (−ΔO 2 /ΔN)=10.6 were used 35 . c(AOU)/ R C corresponds to the DIC increase due to organic matter oxidation, 1/2(ΔA T + c(AOU)/R N ) accounts for the DIC change due to CaCO 3 dissolution-precipitation and ΔA T = A T − A T° is the total alkalinity variation since the water mass was formed. Preformed alkalinity (A T°) and the disequilibrium term that stands for the air-sea CO 2 difference (ΔC dis ) were calculated for each water sample from the mixing proportion of the different water masses. This was carried out by an extended optimum multiparameter analysis (eOMP) 36,37 . The A T° type values for the NACW were calculated using the approach proposed by 38 while those proposed by 39 and 40 were used in the case of the MOW, by considering existing CFC and CO 2 exchange data from the Mediterranean basin, following a previous study 37 . C°T 278 represents the c(DIC) in equilibrium with the preindustrial atmospheric CO 2 molar fraction of 278 ppm and was calculated with the proper dissociation constants [28][29][30] . (ΔC dis ) for both NACW and MOW were considered to be −12 ± 5 μmol kg −1 and 0 ± 5 μmol kg −1 , respectively 11,41 . Uncertainty associated to C ANT calculation was estimated through an error propagation analysis, resulting in ±5.0 μmol kg −1 . An additional OMP analysis was performed to resolve the water masses structure within the MOW and obtain biogeochemical parameters in the LIW and the WMDW. Potential temperature and salinity chosen as end members for water masses characterization were 13.22, 12.8 °C and 38.56, 38.45 for the LIW and the WMDW respectively 13 .
where N j is the concentration of N in sample j, and x ij is the fraction of the water mass i in the sample j.
Here, x ij of each water mass present in the SoG was obtained by the OMP analysis. Standard deviation (SD) of archetypal concentrations was calculated through the formula: Anthropogenic and natural components of pH change. The contribution of anthropogenic and natural drivers affecting total pH variation in each water mass was estimated by calculating the effect of non-biological processes on pH T25 changes (pH nobio ), extracting the biological contribution as follows: is the stoichiometric ratio 35 and ΔpH was given an annual average constant value of 0.002 5 , assuming that surface seawater is in equilibrium with the global mean rate of atmospheric CO 2 increase.

Determination of temporal trends and Statistics. Considering the annual archetypal concentrations
of each parameter during the 11 years period, annual variation trends were determined by ordinary linear regression and standard errors and 95% confidence intervals were calculated for the slopes of the regressions.

Results
Biogeochemistry in waters of the SoG: decadal spatial distribution. Water circulation in the SoG is illustrated in Fig. 2(A-C) where the decadal averaged proportions of the water masses exchanging through the channel have been plotted. The eastward Atlantic inflow (represented here by the fraction of the NACW) invariable occupied the upper layer whereas Mediterranean waters (LIW and WMDW) flew towards the west in depth ( Fig. 2A-C). However, the depth and thickness of each water mass varied along the Strait, which is due to the influence of physical mechanisms that act at different temporal scales, such as tidal currents, winds, atmospheric pressure variations and circulation processes forced by bathymetry that ultimately determine the vertical position of water masses 20 . Thus, the NACW moves upwards toward the easternmost part of the channel, shifting from 200 m at station G1 to around 70 m at station G3 ( Fig. 2A). Worth mentioning is the remarkable variability in water types proportion between the Espartel sill (358 m depth) and the Camarinal sill (285 m depth) (ES and CS in Fig. 1B), which is the consequence of the enhanced mixing and turbulence related to the internal hydraulic jump formed during most tidal cycles downstream of Camarinal [43][44][45] . Moreover, the Atlantic water accelerates in the CS (station G2) and entrains the Mediterranean layer 43 , which results in appreciable changes in the water masses proportion around this area that can be effectively resolved by the OMP analysis ( Fig. 2A-C). Even though horizontal variations in the vertical position of the LIW and the WMDW can be also observed along the channel in response to hydrodynamic features, the former was present at intermediate depths and the latter occupied the bottom layer, which is particularly clear at station G3 (Fig. 2B,C).
Water masses allocation in the SoG clearly conditions the spatial distribution of the carbon system parameters within the channel, as shown by the vertical gradients of the decadal averaged concentrations (Fig. 2D-I) The NACW exhibited the lowest averaged c(DIC) and a marked vertical variability (Fig. 2D), with a mean value of 2127 ± 33 μmol kg −1 . In contrast, the LIW and the WMDW were characterized by higher and similar mean decadal DIC levels, equivalent to 2322.6 ± 12.7 μmol kg −1 and 2311.7 ± 18.8 μmol kg −1 respectively (Fig. 2D), which match previous discrete measurements reported in different regions of the Mediterranean basin 7,9 . A T also increased gradually with depth ( Fig. 2E), with lower decadal averaged concentrations being found in the upper NACW (2310 ± 22 μmol kg −1 ) in relation to those in the LIW and WMDW (2570 ± 14 and 2567 ± 21 μmol kg −1 , respectively). The similar pattern followed by A T and c(DIC) with salinity (Figs 1C and 2D,E) is in agreement with the conservative behaviour of these parameters described in the past 11,40 .
In contrast, pH T25 decreased within the water column (Fig. 2G), and decadal averaged pH T25 values were 7.950 ± 0.043, 7.877 ± 0.015 and 7.891 ± 0.011 in the NACW, LIW and WMDW respectively (Fig. 2G). The mean values obtained in the Mediterranean water masses coincide with earlier observations in the SoG 13 and Western Mediterranean 7 whereas the higher pH variability observed in the NACW has been associated to the influence of biology into the photic zone and changes in evaporation and mixing 46 .
As expected according to the pH T25 distribution, pCO 2 levels increased in the LIW (Fig. 2H) with respect to the rest of water masses, with decadal averaged pCO 2 values being 430 ± 17, 407 ± 13 and 376 ± 30 µatm in the LIW, WMDW and NACW respectively. During the monitoring period, surface Atlantic waters were hence undersaturated with respect to atmospheric CO 2 levels, which gradually rose from 379.8 µatm in 2005 to 400.8 µatm in 2015.
The relatively low pH T25 values characterizing Mediterranean waters are taken as an indication of the active organic matter remineralization occurring in the basin 47 . Our AOU data suggest indeed a remarkable oxidation of organic compounds in Mediterranean water masses, as decadal averaged AOU levels were 65 ± 14 μmol kg −1 in the LIW and 66 ± 7 μmol kg −1 in the WMDW (Fig. 2F). In the case of the NACW, a decadal AOU concentration of 12 ± 22 μmol kg −1 was obtained.
C ANT distribution in the SoG (Fig. 2I) was also influenced by the pattern of water exchange and circulation. The highest C ANT concentrations were attained in the NACW, particularly in the first 100 m depth at the westernmost side of the channel where mixing with the MOW is insignificant (Fig. 2A). Decadal averaged C ANT concentrations were 75 ± 16, 71 ± 11 and 52 ± 8 µmol kg −1 in the NACW, WMDW and LIW respectively.
It is worth noting the appreciable changes in the vertical profiles of the carbon system parameters westwards of station G2 (Fig. 2D-I), which can be attributed to variation in the water masses proportion brought about by the energetic hydrodynamic feature caused the tidal cycle in the CS area ( Fig. 2A-C).
Archetypal concentrations trends. The temporal evolution of archetypal concentrations of the carbon system parameters in the NACW revealed a clear decadal trend (Fig. 3). In particular, pH T25 in this water mass showed a marked decrease over the decade 2005-2015 (Fig. 3A), at an annual rate of −0.0036 ± 0.0005 a −1 ( Table 2). Moreover, AOU and DIC concentrations in Atlantic waters increased annually during the period of study (Fig. 3B,C), at rates of 1.1 ± 0.4 and 3.5 ± 0.6 µmol kg −1 a −1 , respectively (Table 2). Similarly, C ANT concentration in the NACW also rose gradually with time, at an annual rate of 1.5 ± 0.6 µmol kg −1 .
Archetypal concentrations of pH T25 , AOU and DIC in the LIW were quite stable during the monitoring period and a clear decadal tendency was not observed in any case (Fig. 3). In fact, c(AOU) and c(DIC) presented the highest annual archetypal levels in relation to those found in the rest of water masses but with negligible annual variation rates (Table 2). On the other hand, a slight temporal increase in the levels of C ANT could be detected in the LIW, at a rate of 0.3 ± 0.7 µmol kg −1 a −1 .
Annual archetypal concentrations of pCO 2 were in opposition to the archetypal pH T25 levels observed in each water mass (not shown). The upper NACW that was characterized by the lowest decadal average pCO 2 level (Fig. 2H) presented the highest pCO 2 rise during the monitoring period, at a rate of 3.2 ± 0.9 µatm a −1 ( Table 2). In the WMDW, pCO 2 increased at a rate of 1.2 ± 0.7 µatm a −1 whereas the temporal evolution of pCO 2 in the LIW was very smooth ( Table 2).
The highest and lowest contribution of the anthropogenic component to the total pH decrease (∆pH ) nobio was obtained in the NACW and the LIW, respectively ( Table 2). In terms of percentage, the ∆pH nobio in the NACW had a weight of approximately 60% of the total annual decreasing pH rate and around 44% in the case of WMDW. Because of the negligible ΔpH T25 measured in the LIW, pH nobio in this water mass did not provide any significant information.

Discussion
Data collected in the SoG during the 2005-2015 decade confirm previous findings indicating that water exchange in the region results in a C ANT import from the Atlantic towards the Mediterranean basin 10,11 as the NACW was always C ANT enriched with respect to the MOW. Our periodic observations are also in agreement with the finding that the North Atlantic is the largest ocean sink for anthropogenic carbon 2,48 . The absorption of C ANT by the NACW contributed to a noticeable decadal reduction of its pH, with the contribution of the anthropogenic component representing 60% of the total pH decline.
The rise in C ANT content over time in the NACW spotted at the SoG occurred at a rate of 1.5 µmol kg −1 a −1 ( Table 2), which matches quite well the decadal rate estimated previously in this water mass 49 . This increment in the quantity of anthropogenic carbon was accompanied by parallel increases in the concentration of DIC and AOU, also suggesting an enhanced organic matter oxidation during the 2005-2015 decade. Before crossing the Strait, the NACW transits through the Gulf of Cádiz (GC in Fig. 1A), a basin characterized by a high primary production 50 . Therefore, regional inputs of organic matter are likely to occur, which could well lead to an intensification in the ecosystem respiration and inorganic carbon accumulation with time. In addition, a considerable discharge of inorganic carbon from the Guadalquivir river estuary to the continental shelf of the GC has been measured 51,52 , which could be accounted as an extra DIC source for the NACW. These external inputs would contribute to the natural component affecting the decadal pH change in this water mass, which represented 40% of the total pH reduction. It is also worth noting that pCO 2 in the NACW rose at a higher rate (3.2 ± 0.9) than the atmospheric CO 2 annual trend (approximately 2 ppm) 53 . This divergence can be related to the air-sea CO 2 disequilibrium driven by surface heat fluxes and the balance between the biological carbon uptake and CO 2 outgassing 46 that becomes larger with time 54 . Taking into account the atmospheric pCO 2 levels during the monitoring period and an average 90% equilibrium with the anthropogenic fraction 54 , the theoretical pH values in the NACW that would result exclusively from atmospheric CO 2 uptake were calculated and compared with the estimated pH nobio values. Comparison performed by a Mean Squared Error (MSE) revealed a difference of 0.0001 between both estimates and a mismatch between dissolved and atmospheric pCO 2 concentrations equivalent to 8 years. This finding provides more evidence for local ecosystem drivers of pH fluctuations, as described elsewhere, particularly in coastal regions 17,53 .
Multi-decadal declines of pH in surface Atlantic waters have been observed in the open ocean 5 yet few studies have quantified the relative importance of biological processes to modulate the pH change over time 15 . The acidification rate calculated here in the NACW from 2005 to 2015 (0.0036 pH a −1 Table 2) is of the same order of magnitude than others obtained in distant Atlantic basins although higher 5 (e.g. 0.0026 and 0.0025 pH a −1 in the Irminger Sea and CARIACO basin, respectively). Interestingly, the rate of change of the anthropogenic component (pH nobio , 0.0021 pH a −1 , Table 2) is quite similar to the total rates attained in the Atlantic time series where the rates of dissolved pCO 2 rise resembled the current atmospheric CO 2 increase 5,53 . The influence of regional www.nature.com/scientificreports www.nature.com/scientificreports/ biogeochemical processes on the carbon content of the NACW possibly amplified the temporal pH change and accounted for by the increased acidification trend found with our sustained observations, although this hypothesis remains to be tested.
On the other hand, the decadal rate of pH decrease in the WMDW is comparable to that reported for the Mediterranean finger print in the North Atlantic 55,56 and lies between the limits established in Mediterranean bottom waters 10 . It can be argued that the pH tendency observed here is not fully supported by the regression statistics. However, it has been recognized that the statistical significance of seawater carbonate chemistry trends in time-series can be weak due to several reasons, such as the irregularity of sampling, non-uniform time intervals between cruises and seasonality in the measurements 5,53 . A further analysis conducted to remove seasonality and identify potential bias of sampling in our databse revealed no seasonal effect on the estimated trends (not shown).
Previous works have already concluded that the MedSea is experiencing an acidification process 9,10,12-14,57 although the pattern of pH change is not uniform at a basin scale 9,10 , which can be attributable to both the west to east decreasing productivity gradient 58 and water circulation mechanisms 9,12 .
Each winter, deep waters are formed in both the Eastern and Western Mediterranean. The recently formed waters mix with older and resident water masses, which change their chemical and physical properties; hence, it is not simple to anticipate the effect of water formation events on pH changes. In the particular case of the WMDW, several mechanisms impacting pH at distinct temporal scales must be taken into account, which ultimately complicates the identification of a significant temporal tendency. The WMDW is formed regularly by two processes, deep convection and dense shelf water cascading. The former induces direct downward transport of newly formed organic matter from the surface layer, which is enriched in labile and easily oxidizable material whose decomposition leads to significant CO 2 release 59,60 . Hence, the WMDW receives a significant amount of pCO 2 from organic matter degradation that is transported down during the formation event. Furthermore, the active ecosystem metabolism occurring in the Alboran sub-basin (AS in Fig. 1) where this water mass resides before leaving the basin 61 is likely to increase its pCO 2 bulk. In the Alboran Sea, a high primary production in the upper layer 62 is coupled to a significant export to deep of organic material 58 and a strong remineralization in the bottom layer 59,63 , which is indeed occupied by the WMDW 61 (Fig. 2C). Therefore, the WMDW must gain a significant amount of pCO 2 from natural processes before crossing the SoG. Both mechanisms, decomposition of labile organic matter at formation site and the cascading effect during residence time in the Alboran sub-basin would explain the strong biological component (around 60%) of the total pH decline found. The high AOU values measured in the WMDW here and concordant with previous measurements 18 support this assumption.
Nevertheless, the process of the deep water mass formation in the Gulf of Lions is subject to a marked interannual variability, and thus, pH temporal evolution in the WMDW may not be regular, as evidenced in our time series. Our database identifies a clear pH decline during 2013 (Fig. 3A), which is also evident in the LIW. This particular year marks the occurrence of the Western Mediterranean Transition (WMT) at the SoG 64 . This phenomenon started with a strong WMDW production in the Gulf of Lions during the extreme cold winters of 2005 and 2006 that originated a remarkable deep convection. This deep water enriched in labile organic carbon 18 arrived to the Alboran Sea and remained in the bottom layer. A subsequent cycle of deep convection during the harsh winters of 2012 and 2013 resulted in newly formed WMDW that was warmer, saltier and denser than the old one 64 and caused its uplifting and propagation towards the Strait upon arrival in Alboran. The lower pH values observed in the WMDW in 2013 can be interpreted as the footprint of the older (and more acidified) fraction of this water mass after being almost entirely replaced by the more recent fraction. Moreover, it was inferred that the younger WMDW spilled over the SoG in 2015 64 , which would explain the rise in pH observed in this water mass that year (Fig. 3A), corresponding to the signature of a less acidified fraction. The decrease observed in both C ANT and c(AOU) in 2015 (Fig. 3B,D) is consistent with the presence of a water mass that was rapidly formed and sunk and remained in the AS a lower period of time. This interannual variability in the WMDW time series possibly blurs the identification of statistically significant long term tendencies. Changes in the proportion of fractions of WMDW carrying different levels of both natural and anthropogenic CO 2 would also account for the discrepancy found between the acidification rate provided here and the one we observed previously in the SoG using three years (from 2012 to 2015) of continuous pH measurements within the MOW 13 . In fact, in the former study considering a shorter period of time, we reported a higher annual pH decline and drastic changes in the proportion of the WMDW during 3 years, which can be related to the influence of the WMT and the subsequent cycle of deep formation events. A thorough investigation of the hydrographic profiles available in our database would certainly contribute to unravel the observed pH variability in the WMDW and strengthen our conclusions but it is not the scope of the present study.
The WMT caused a near-complete renewal of the WMDW 64 and possibly also some mixing with the LIW, thereby lowering its pH. Blurring of the LIW signal within the MOW has been identified in the SoG and diagnosed through salinity changes 20 , thus mixing between both water masses in the area is a common phenomenon.
Interestingly, during the monitoring period, the LIW exhibited similar pH T25 values than those in the WMDW although the temporal evolution of its pH was smoother, as observed in our earlier study 13 . The LIW originates in a distant sub-basin and crosses the entire Mediterranean basin before arriving at the SoG. Therefore, during transit, changes in pH, AOU and pCO 2 are expected to proceed due to physical and biological processes. The LIW route from formation site is complicated as it follows different pathways 65,66 : first it moves westward from the Levantine basin forming two veins, one that flows to the Adriatic Sea and a second one that crosses the Strait of Sicily and enters in the Tyrrhenian. This vein eventually emerges from the Sardinian Channel, reaches the Ligurian Sea, the Gulf of Lions and finally flows towards the Algerian basin 65 . During this journey, changes in biogeochemical parameters, mainly dissolved organic carbon (DOC) and AOU, have been indeed found in response to the different trophic status of the sub-basins crossed 18 . In particular, the LIW passes by areas of strong organic matter remineralization that cause minimum contents of DOC and maximum levels of AOU 18 . The c(AOU) measured in the LIW during the 2005-2010 decade (within the range 60-70 µmol kg −1 ) are possibly the final www.nature.com/scientificreports www.nature.com/scientificreports/ result of all the biological changes experienced during transit, which modified original pH as well. Furthermore, as the LIW remains at intermediate depths in the SoG (Fig. 2B), Algerian and Alboran sub-basins 61 the additional contribution of cascading processes to its pCO 2 bulk cannot be totally neglected. Therefore, re-ventilation mechanisms suffered by the LIW, vertical transport of carbon by diffusion and mixing with the WMDW must alter its original biogeochemical characteristics and consequently lower its pH.
In fact, when the ages of the WMDW and the LIW are estimated, similar values are obtained and equivalent to 32 ± 8 and 34 ± 8 years, respectively (MSE < 0.0001). Alternatively, age calculation based on available CFC data 39 yields estimates of 18 years for the WMDW and 21 years for the LIW. Although there is no consensus about the exact ages of Mediterranean water masses and discrepant values have been provided 67,68 , it is clear that our calculated ages for the LIW are lower than that expected for a water mass that has crossed the entire basin. Actually, the age of the LIW in the Western basin has been proposed to be in the range of 82-120 years 67 .
These findings and the fact that appreciable C ANT levels were measured in the LIW support the occurrence of re-ventilaton and/or mixing with younger water masses. If we consider that the LIW currently spotted at the SoG corresponds to a water mass that sank originally in the Eastern basin and lost contact with the atmosphere roughly one century ago, it is plausible to assume that its C ANT content would be negligible, as it was exposed to a very low atmospheric CO 2 concentration upon formation (resembling the CO 2 pre-industrial level). During transit, its position within the water column prevented any new contact with the atmosphere, thus hampering the absorption of the rising atmospheric CO 2 . Therefore, the LIW must gain anthropogenic carbon by other processes, and the most likely candidates are re-ventilation and mixing with more recently formed water masses, such as the Cretan Intermediate Water (CIW) in the Ionian Sea 18 and the WMDW in the SoG 20 . This would alter the original biogeochemical characteristics of this water mass and in turn, estimations of its age and anthropogenic signals.
As the drivers of anthropogenic forcing and natural variability are not entirely constrained and their signals cannot be discerned at a statistically significant level, identification of acidification trend in the LIW would require a longer time period of periodic observations, as it has been recently shown in a number of time series 53 .

Conclusions
Our study evidences the presence of an anthropogenic impact on two major Mediterranean water masses, as appreciable C ANT levels were measured in the WMWD and the LIW using 11 years of sustained observations at the Strait of Gibraltar. However, a clear acidification trend could only be identified in the former, with the contribution of the anthropogenic driver to the total pH change representing 40% of the decadal decline. As the LIW is an old water mass with no contact with the current atmosphere, its C ANT content may be attributable to further contamination during transit across the basin from formation site. Longer monitoring period than that encompassed by our study and observations in the Eastern Mediterranean are required to detect long term pH variations in the LIW and constrains the relative importance of the natural and anthropogenic signals on its pH evolution.
It is plausible that the acidification phenomenon in the MedSea could be exacerbated in the future, as our assessment also shows that the incoming Atlantic water suffered a marked acidification trend over the decade. In this case, reduction in the pH of the NACW was dominated by the anthropogenic component, which accounted for by 60% of the total decrease. Regional biogeochemical processes might amplify natural pH variability in this water mass, resulting in higher acidification rates in relation to others observed in surface Atlantic waters.
Regardless of the relative contribution of the drivers affecting pH decline in Atlantic and Mediterranean waters, our findings are in support of the current concerns regarding the fate of marine ecosystems due to the increase in acidity of the seawater. Particularly, iconic habitats of cold water corals of the North Atlantic 69 and Mediterranean Sea 70 are expected to be highly impacted by the phenomenon, which is now corroborated by our assessment.
Hence, even though further work should be directed at improving upon our estimates by increasing the length of sustained observations, our decadal records still provide a starting point for future calculations of ocean acidification trends in Mediterranean waters and reveal the importance of long time series in this marginal sea.

Data availability
The dataset generated during and/or analysed during the current study are available from the corresponding author upon reasonable request.