Understanding the Implications of Hydrographic Processes on the Dynamics of the Carbonate System in a Sub-Antarctic Marine-Terminating Glacier-Fjord (53°S)

The biogeochemical dynamics of fjords in the southeastern Pacific Ocean are strongly influenced by hydrological and oceanographic processes occurring at a seasonal scale. In this study, we describe the role of hydrographic forcing on the seasonal variability of the carbonate system of the Sub-Antarctic glacial fjord, Seno Ballena, in the Strait of Magellan (53°S). Biogeochemical variables were measured in 2018 during three seasonal hydrographic cruises (fall, winter and spring) and from a high-frequency pCO2-pH mooring for 10 months at 10 ± 1 m depth in the fjord. The hydrographic data showed that freshwater input from the glacier influenced the adjacent surface layer of the fjord and forced the development of undersaturated CO2 (< 400 μatm) and low aragonite saturation state (ΩAr < 1) water. During spring, the surface water had relatively low pCO2 (mean = 365, range: 167 - 471 μatm), high pH (mean = 8.1 on the total proton concentration scale, range: 8.0 - 8.3), and high ΩAr (mean = 1.6, range: 1.3 - 4.0). Concurrent measurements of phytoplankton biomass and nutrient conditions during spring indicated that the periods of lower pCO2 values corresponded to higher phytoplankton photosynthesis rates, resulting from autochthonous nutrient input and vertical mixing. In contrast, higher values of pCO2 (range: 365 – 433 μatm) and relatively lower values of pHT (range: 8.0 – 8.1) and ΩAr (range: 0.9 – 2.0) were recorded in cold surface waters during winter and fall. The naturally low freshwater carbonate ion concentrations diluted the carbonate ion concentrations in seawater and decreased the calcium carbonate saturation of the fjord. In spring, at 10 m depth, higher primary productivity caused a relative increase in ΩAr and pHT. Assuming global climate change will bring further glacier retreat and ocean acidification, this study represents important advances in our understanding of glacier meltwater processes on CO2 dynamics in glacier–fjord systems.

The biogeochemical dynamics of fjords in the southeastern Pacific Ocean are strongly influenced by hydrological and oceanographic processes occurring at a seasonal scale. In this study, we describe the role of hydrographic forcing on the seasonal variability of the carbonate system of the Sub-Antarctic glacial fjord, Seno Ballena, in the Strait of Magellan (53°S). Biogeochemical variables were measured in 2018 during three seasonal hydrographic cruises (fall, winter and spring) and from a high-frequency pCO 2 -pH mooring for 10 months at 10 ± 1 m depth in the fjord. The hydrographic data showed that freshwater input from the glacier influenced the adjacent surface layer of the fjord and forced the development of undersaturated CO 2 (< 400 matm) and low aragonite saturation state (W Ar < 1) water. During spring, the surface water had relatively low pCO 2 (mean = 365, range: 167 -471 matm), high pH (mean = 8.1 on the total proton concentration scale, range: 8.0 -8.3), and high W Ar (mean = 1.6, range: 1.3 -4.0). Concurrent measurements of phytoplankton biomass and nutrient conditions during spring indicated that the periods of lower pCO 2 values corresponded to higher phytoplankton photosynthesis rates, resulting from autochthonous nutrient input and vertical mixing. In contrast, higher values of pCO 2 (range: 365 -433 matm) and relatively lower values of pH T (range: 8.0 -8.1) and W Ar (range: 0.9 -2.0) were recorded in cold surface waters during winter and fall. The naturally low freshwater carbonate ion concentrations diluted the carbonate ion concentrations in

INTRODUCTION
As global climate change continues, there is increasing awareness of the influence of anthropogenic CO 2 on the melting of glaciers in polar and sub-polar marine ecosystems (Meredith et al., 2019). The input of glacial meltwater into fjords in these ecosystems strongly modulates their biogeochemistry with implications for sea-air CO 2 fluxes and ocean acidification Fransson et al., 2013). Seasonal ice-melting events change the water column's physical and optical properties, resulting in stratified, higher turbidity fjords. The suspended material decreases light availability for phytoplankton growth, causing a non-linear response of primary production rates (Hopwood et al., 2020) and potential feedback for the inorganic carbon cycle.
Northern Patagonian fjords (41°S) are typically CO 2 sinks during the warm and productive summer season and CO 2 sources in the cold and low productivity winter season (Torres et al., 2011b). While winter's high CO 2 levels, cold temperatures and low salinity surface waters result in strongly "corrosive conditions" for calcium carbonate (Feely et al., 2018), in the warm, less rainy and sunny season, primary productivity and high temperatures lead to higher levels of calcium carbonate saturation (W) (Alarcoń et al., 2015).
High freshwater runoff of low alkalinity in Northern Hemisphere fjord systems during spring and summer produces low aragonite saturation state, low pH and favors the flux of CO 2 to the atmosphere. In some cases, the low aragonite saturation state and low pH are not compensated by the mixing of higher alkalinity water from depth, or by the uptake of carbon by primary productivity (Chierici and Fransson, 2009). Freshwater runoff in these fjords may also modulate the availability of organic matter (Dissolved Organic Carbon, DOC), macronutrients (N, P, Si), and micronutrients (Fe) in the euphotic zone, which are all important to phytoplankton and bacterial processes (Fransson et al., 2016). Under winter conditions, the coupling between physical (ice formation and physical mixing), chemical (increase in macronutrients and pCO 2 , but decreasing availability of CO 2− 3 ), and biological processes (ratio of primary production to respiration <1) may decrease the CaCO 3 saturation state in surface waters .
In the Sub-Antarctic region (50-55°S), glacier retreat (Bown et al., 2014) along with important atmosphere-ocean interactions (Garreaud, 2018), may profoundly influence biogeochemical processes. In this southernmost region (Figure 1), seasonal changes in glacier mass may simultaneously influence biogeochemical and physical oceanographic features. First, like Northern Hemisphere fjords, stratification and mixing processes impact nutrient fluxes. Second, melting glacier ice creates low temperature and low salinity surface water, generating large vertical and horizontal gradients in density. In fjords with a very shallow sill, tidal dynamics may produce Bernoulli aspiration (Kinder and Bryden, 1990), resulting in the injection at depth of nutrients to the inner portion of the fjord (Torres et al., 2011b). Reduced winter rates of glacier melting may result in longer residence times of deep water in the inner portion of the fjord, seasonally reducing the ventilation of deep waters. All of these factors may play a role in modulating the biological carbon pump, including the chemical speciation of the carbonate system in the fjord region.
Fjord-glacier systems are singular but diverse in their characteristics and thus likely also in their dynamics. Factors such as the hydrographic characteristics of the watershed, glacier dynamics, bathymetry, wind and sun orientation, latitude, location of the ice edge (e.g., glacier-terminating fjord, glacier edge on the coast or far from the shoreline) may be critical to understand the response of fjords to climate forcing. It has been suggested that the increase in air temperature and atmospheric pCO 2 in recent years may result in increased cold freshwater input from the glaciers (Bown et al., 2014), with a large potential to perturb the carbonate system speciation (Meire et al., 2015;Hopwood et al., 2020) and thus marine biota (Orr et al., 2005;Kurihara, 2008;Grear et al., 2017;Giesecke et al., 2019). Current projections of how changes in the ocean carbonate system may influence biological processes (i.e., photosynthesis and microbial respiration) in the future still remain speculative. Hydrological and biogeochemical information for Seno Ballena, a glacier-fjord system, will be relevant for explaining biological responses to coastal acidification interacting with other climate stressors (i.e., warming and reduced salinity due to increased freshwater flow). In this study, the main objective is to examine the carbonate system's seasonal dynamics due to the influence of glacier meltwater and oceanic waters on the upper surface layer along the fjord. This study specifically presents the seasonal hydrography and time series of pCO 2 , pH T , and O 2 in the Seno Ballena glacier Fjord from austral fall through spring (March to December, 2018). Frangoṕulos et al. (2007); Figure 1). It is 18 km long from the head, where the Santa Ineś glacier is located, to the outside of the sill in the direction of Carlos III Island, Strait of Magellan ( Figure 1). The lower (southwest) area of the fjord is approximately 9 km long and less than 1 km wide. A sill reaching a depth of 2-3 m from the surface during high tide is located 7 km from the glacier (Valle-Levinson et al., 2006). The bathymetry is highly variable from the inside to the outside of the fjord. Near the glacier depths are between 20 to 50 m, followed by a deeper area of 150 m up to the sill; past the sill depths of more than 200 m are found. This system is strongly influenced by freshwater input from (1) the melting of the glacial ice located at the head of the fjord; (2) a river of glacial origin; (3) melting snow of the adjacent mountains, that forms several small tributaries along the fjord (Supplementary Figure 1); and (4) high rainfall (4500-6000 mm y -1 ; Meier et al., 2018). These hydrological processes are the primary controls of the stratification of the fjord, characterized by a thin (~10 m) surface layer of lowsalinity water (24 -28 S p ).
The fjord is within the Sub-Antarctic fjord and channel system of Patagonia (Figure 1), for which a strong influence of Modified Sub-Antarctic Waters mass (MSAAW) has been proposed (Supplementary Figure 2). This water mass occupies a significant fraction of the Strait of Magellan and Sub-Antarctic fjords, having practical salinity values (S p ) between 31 and 34 (Silva et al., 1998;Sievers et al., 2002;Valdenegro and Silva, 2003). When this water mass mixes with freshwater runoff from rivers and glaciers it forms estuarine waters, which are classified into three categories: the saline water of the estuary is more than 66% seawater (21-31 S p ); estuary-brackish water is 33-66% seawater (11-21 S p ); and freshwater in estuaries is less than 33% seawater (2-11 S p ) (Sievers and Silva, 2008). MSAAW has low dissolved silicate (DSi) compared to nitrate concentration. In contrast, given the geology of the drainage basin, the siliceous crystalline batholiths of Patagonia, the continental shelf waters of western Patagonia are characterized by high DSi, low nitrate, low alkalinity and low Ca 2+ (Hervéet al., 2007;Torres et al., 2020).

Hydrographic Measurements
Three seasonal oceanographic cruises were conducted along the fjord in March (fall), August (winter), and December (spring), 2018. Five to ten stations along the main axis of the fjord were visited during each cruise ( Figure 1). Hydrographic profiles (in situ temperature and salinity) were obtained using a CTD (SeaBird SBE model 19). The CTD salinity data in spring were not used due to CTD technical problems, and consequently salinity was determined from water samples taken at distinct depths (0, 5, 10, 25, and 50 m) and measured using a YSI-Pro30 probe at constant temperature previously calibrated with an IAPSO standard seawater (35 S p ); the nominal uncertainties on these measurements were ±1%. The water masses were identified by the mixing triangle method (Mamayev, 1975), using the information from the fall and winter campaigns to confirm the presence of MSAAW (Silva et al., 2009). Throughout practical salinity (PSS-78) is reported as S p (UNESCO, I. 1981).

Water Samples and Analysis
Samples for the determination of carbonate chemistry parameters, inorganic nutrients and autotrophic biomass were taken at 5-8 stations, at 0, 5, 10, 25, and 50 m depths. Samples were also taken from glacier ice as well as river-glacier water and small freshwater tributaries along the fjord to determine freshwater end-members of salinity and alkalinity. Water samples were collected in 10 L Go-Flo bottles and stored in the dark at low temperature (< 7°C) in 250 mL gas-tight containers for the analysis of total alkalinity (A T ) and pH T . Seawater samples for A T analysis were poisoned with 50 ml of a saturated mercuric chloride solution (HgCl 2 ; Dickson et al., 2007). pH T was analyzed soon after collection (<12 h), using purified m-cresol purple as an indicator at 25.0°C with an Ocean Optics STS-Vis (350 -800 nm) spectrophotometer (Byrne et al., 1988). When possible, pH T was also measured with an aquatrode probe (Metrohm TM) with a single point calibration (pH Tris buffer pH=8.089 at 25.0°C and salinity 35) (Dickson and Goyet, 1994). Mean differences from three replicate samples between methods were 0.01 pH units in the 28 -31 salinity range (spectrophotometer and electrode pH). Tris buffer was made for salinity 35 whereas the field salinities ranged from 27.4 to 31.2. Nevertheless, we have determined empirically (using buffers at salinity 35 and 25, and this electrode) that the uncertainties in pH due to differences in salinity between buffer and sample is lower than 0.01 pH units. The performance (i.e., Nernstian slope) of the electrode was better than 99.8% of the theoretical value.
A T analysis was performed with the automatic open-cell potentiometric titration method (Haraldsson et al., 1997). This technique allows A T to be obtained quickly and very precisely with a small sample volume (~40 ml) by adding hydrochloric acid (0.05 M HCl; Merck Titrisol ® ) in increasing volumes with a Dosimat 665. The reading was done with a combination Ross electrode (Orion 8102BN). The end-point was determined by the Gran method (Gran, 1952) according to Haraldsson et al. (1997). Each sample was analyzed twice and the average and standard deviation are reported. Seawater distributed by Dr. Andrew Dickson's laboratory was used to verify the A T estimates.
Filtered (0.7 μm glass fiber filters; Whatman GF/F) samples were frozen (-20°C) for analysis of soluble reactive phosphate (SRP), NO − 3 and dissolved silicate (DSi), using the manual techniques recommended by Parsons et al. (1984). The inorganic nutrients and A T were analyzed in the carbonate laboratory at the Centro de Investigacioń de Ecosistemas de la Patagonia (CIEP, Coyhaique). Chlorophyll-a concentration of seston larger than 0.7 mm was extracted with 90% acetone and measured using a fluorometer (Turner P700) following Parsons et al. (1984).
Estimates of in situ pH T , pCO 2 and aragonite saturation state (W Ar ) were calculated using the chemical speciation model program CO 2 SYS (Pierrot et al., 2006), using the dissociation constants (K 1 and K 2 ) estimated by Mehrbach et al. (1973), modified by Dickson and Millero (1987) for salinity = 20 -40 and temperatures 2 -35°C. Our values fit into these ranges of S P and temperatures. We used the equilibrium constants for KHSO 4 determined by Dickson (1990) and those for Total Boron determined by Uppström (1974). The dissociation constants provided by Dickson and Millero (1987) have been shown to be consistent with results of previous studies carried out in the Chilean Patagonia fjords (Torres et al., 2011b;Alarcoń et al., 2015;Vergara-Jara et al., 2019;Torres et al., 2020). Input data for CO 2 SYS used T, S P , A T , pH T @25°C and nutrients (SRP and DSi) (Kim and Lee, 2009). We measured only two parameters in discrete samples (pH T and A T ) and two in the in situ sampling (pH T and pCO 2 , see section 2.5); therefore we cannot estimate the internal consistency of the measured carbonate system parameters. Organic matter (organic acids, fulvic and humic acids) could be a significant source of non-carbonate alkalinity, produced during decomposition of organic matter in sediments carried away by inland waters (Lukawska-Matuszewska, 2016). Therefore, care was taken to interpret calculated carbonate system parameters since they may result in overestimation of DIC and pCO 2 . Organic alkalinity does not play an important role in this system, given the low average contribution to the total alkalinity of the freshwater end-members (5.8 ± 16.0 mmol kg -1 ).

Water Mass Distribution Analysis
To quantify the relative contributions of freshwater from the Santa Ines glacier and adjacent river-glacier from the oceanic-type water mass (estuarine-brackish: MSAAW), an optimum multiparameter (OMP) analysis was performed, which uses a simple model of linear mixing to calculate the fraction of the source-water types. Two conditions are taken into account to apply the method: (1) all the calculated fractions are positive and (2) the sum of all the fractions is close to 100% (conservation of mass). In our analysis, we found that only two water masses were not sufficient to reproduce the properties of the samples. Thus, we assume mixing that involves 3 water masses: (1) near-surface and (2) deep inner fjord water masses, whose property definitions are determined from the average of a set of observations made during seasonal synoptic sampling, and (3) the oceanic water mass (MSSAW) that has traditionally been identified for this study area (Sievers and Silva, 2008;Silva et al., 2009).
To characterize the inner fjord water masses, propertyproperty diagrams were made of the observations of salinity, temperature, A T and nutrients to select values that represent the source water types. The observations used for the Inner Surface Source Water (ISSW) were from the surface sampling station closest to the glacier, to represent the properties of the freshwater from the glacier-ice and river-glacier. The observations used to characterize the Inner Bottom Source Water (IBSW) were from sampling stations 2 at 50 m depth (2 km from the glacier). The oceanic water mass that enters this Sub-Antarctic zone is the Sub-Antarctic water mass (SAAW), which is modified within the Strait of Magellan by mixing with continental freshwater (runoff from glaciers and rivers), forming MSAAW and is found from the surface to 150 m depth. The properties of the MSAAW were established based on reported information (Sarmiento et al., 2004;Sievers and Silva, 2008;Silva et al., 2009;Llanillo et al., 2012;Torres et al., 2014;Forceń-Vaźquez et al., 2021). The strong salinity and alkalinity relationship reported for the Patagonian archipelago interior sea (41-56°S: Torres et al., 2011a;Torres et al., 2020) makes it possible to calculate alkalinity of the MSAAW using the equation: A Tsal (mmol kg -1 ) = 66.4 × Salinity + 1.1, (R 2 = 0.93, n=81; Figure 2; Supplementary Table 1). The regression is close to the equations reported for the Sub-Antarctic region A T (mmol kg −1 ) = 61.9 × Salinity + 151 and A T (mmol kg −1 ) = 61.8 × Salinity + 145; Torres et al. (2011b) ( Figure 2). The intercept (1.1 mmol kg −1 ) is a reasonable estimate of the freshwater end-members (glacier-ice and small tributaries along the fjord) with the value mean (A T of 5.8 ± 16.0 mmol kg −1 ; n = 7, salinity 0.0 ± 0.1 S P ).
The water masses proposed in this analysis have ranges of values, thus 5 source types of water (SWT) were used as described in Table 1 (ISSW 1 , ISSW 2, ISBW 1 , MSAAW 1, and MSAAW 2 ). The weights for each variable were chosen taking into account the extent of their conservative behavior and their spatial and temporal variability (S p = 25, Temp = 25, DSi = 5, NO 3 = 5, A T = 25). The salinity and A T data were assigned relatively high weights, because they are conservative in the temporal and spatial scales considered here. Temperature and nutrients were assigned lower weights; temperature due to its high seasonal variability caused by heating and cooling and nutrients due to the changes caused by biological processes. The linear system of mixing equations, that were solved using the classic OMP analysis (MATLAB -version 1.2.0.0; Karstensen, 2013), are the following: where T obs , S obs , A Tobs , DSi obs and NO 3obs are the values observed during the seasonal sampling and the R ∑ are their respective residuals. The values T i , S i , A Ti , DSi i and NO i (i = a, b, c, d and e) define the fixed parameter values of the five source-water types, that define the three water masses (ISSW, IBSW and MSAAW).
The five tracers plus conservation of mass leaves little redundancy for five water types (Karstensen and Tomczak, 1998), since we need to include the fifth water type (IBSW) as it is distinctly different than the other four. x i is the fraction for each data point. Equation (1f) is the mass conservation constraint.

pH-pCO 2 Mooring Sensors
High-resolution measurements of pCO 2 and pH T were recorded in situ using the autonomous submersible SAMI-CO 2 and SAMI-pH sensors (Sunburst Sensors, LLC) (DeGrandpre et al., 1995), temperature and conductivity using SBE37 sensors and dissolved oxygen (DO) using an Aanderaa 5331A sensor. SAMI-CO 2 measures the partial pressure of carbon dioxide (pCO 2 ) in water from 200 to 1000 matm (precision <1 matm, accuracy~10 μatm), while SAMI-pH measures pH T on total proton concentration scale in the marine pH range of 7 -9 for salinity ranges 25 -40 S P (Clayton and Byrne, 1993). The two sensors use a high-precision colorimetric reagent method (http://www. sunburstsensors.com/). The factory SAMI-CO 2 calibration, based on coincident measurements of pCO 2 using gas-equilibration with infrared detection, was used for this study. The SAMI-pH is not factory calibrated but was validated using a Tris buffer at 25°C (accuracy and precision~± 0.003). All of the sensors were placed at 10 ± 1 m and mounted on one single steel frame with the sensors' water intakes in the same vertical position. The sensors were deployed in the lower part of the pycnocline on the outer side of the sill ( Figure 1). Mooring measurements started early in austral fall (March, 2018) and ended in austral spring (December, 2018). All sensors recorded data every 4 hours. Discrete samples of the three seasonal oceanographic cruises (March, August, and December 2018) were used for data quality control.  , where DIC ST.10-10m is station 10 at 10 m that is near pure MSAAW water, f ISSW is fraction of ISSW at station 10 to 10 m and f MSAAW is fraction of MSAAW at station 10 to 10 m. Subscripts 1 and 2 in the water mass indicate slight variations in their physical-chemical characteristics.
The SAMI-pH and A Tsal data were used to compute the additional carbonate system parameters because this pair was shown to yield good accuracy (Cullison Gray et al., 2011). The total alkalinity (A Tsal ) values derived from salinity carry an error of 42 mmol kg -1 . This error is estimated from the values of the A Tsal time series for the days that coincide with the seasonal A T-measured samples measured in fall, winter, and spring. For this calculation, A T-measured data were collected from a sampling station near the mooring at a depth of 10 m.
All inorganic carbon parameters were calculated with CO 2 SYS software (Pierrot et al., 2006) as described above (see item 2.3). Nutrient data were not included in the computations due to the lack of continuous measurements during the study. The calculated pCO 2 changed by~1.0 matm when the highest observed levels of SRP (1.76 mM) and DSi (5.47 mM) were included in CO 2 SYS.

CO 2 Flux Determination
The air-sea CO 2 flux was determined using the diffusive boundary layer model, through the volumetric flux equation expressed in terms of CO 2 partial pressure: where F is the air-sea flux in units of moles area −1 time −1 , K is the gas transfer velocity in units of length time −1 (Wanninkhof, 1992), K O is the solubility coefficient of CO 2 (mol m −3 atm −1 ), estimated from in situ salinity and temperature according to Weiss (1974) and the pCO 2W -pCO 2a is the difference between the air and sea surface pCO 2 values in μatm. Results with positive values indicate that there is a flux from sea to atmosphere. The gas transfer velocity of CO 2 was calculated using the revised relationship recommended by Wanninkhof (2014): where <U 2 > is the mean squared wind speed (m s -1 ). The wind speed data were obtained from a meteorological station near the study area (Chile Meteorological Directorate, Alberto Hurtado School station, coordinates: -53.16694°S, -70.94528°W). The updated Wanninkhof (2014) parameterization was used in the study of the gas transfer rate estimation in Seno Ballena fjord because it considers the most recent advances in the quantification of the input parameters, improvements in the wind speed products and it provides good estimates for most insoluble gases in intermediate wind speed ranges (3 -15 m s 1 ). The mean wind speed in Seno Ballena Fjord was 4.04 m s -1 during the study period (range of 3.73 to 6.04 m s -1 ). Sc is the Schmidt number which accounts for the difference in molecular diffusivity between gases, with the leading coefficient updated for seawater (35 S P ) and freshwater (0 S P ) at temperatures ranging from -2 to 40°C by Wanninkhof (2014). Sc is calculated using: where t is temperature (°C) and A, B, C, D, and E are fitting coefficients (Wanninkhof, 2014).
The formula described by Dickson et al., (2007) was used to calculate the atmospheric pCO 2 values in humid air and the water vapor pressure values: where pCO 2 (dryair) is the air pressure at sea level taken from the Earth System Research Laboratory database (National Oceanic and Atmospheric Administration Marine Boundary Layer Reference 53.1 to 17.5°S; www.esrl.noaa.gov/gmd/ccgg/ mbl/data.php), interpolated for the year 2018 and corrected with local barometric pressure. These values assume clean marine air, and therefore a possible error is created due to terrestrial effects on pCO 2 , mainly of the export of organic and inorganic carbon from land runoff, rivers and glaciers (Lafon et al., 2014). The error is not readily quantifiable because no regional pCO 2 data are available. P w is the equilibrium water vapor for the in situ temperature (°C) and salinity (S P ) (Forstner and Gnaiger, 1983).

Temperature Effect on pCO 2
Temperature changes affect the partial pressure of seawater due to its solubility. For this reason, the values of the pCO 2 time series were adjusted to a mean temperature, using the equation of Takahashi et al. (2009) (Eqn. 6). The temperature coefficient (0.0459) for cold (-1.8 to 10°C) and less saline (30 < S P < 35) water (Ericson et al., 2018) was used, as these ranges match the observed salinity and temperature ranges of the Seno Ballena fjord.
where pCO 2obs are the values of the pCO 2 time series. The mean temperature during the study period (March to December 2018; T ave ) was 7.67°C, and T obs is the measured temperature in degrees Celsius. This equation makes it possible to determine if changes in pCO 2 in spring are related to phytoplankton bloom events rather than just to temperature changes.

Apparent Oxygen Utilization
Apparent oxygen utilization is calculated to estimate the effect of biology on oxygen concentrations, eliminating the effect of temperature and solubility. It was calculated using the following equation: This equation represents the change in oxygen since a mass of water was last in contact with the atmosphere, assuming that the oxygen concentration at the surface was at equilibrium with the atmosphere (DO sat ). DO is the measured dissolved O 2 . AOU estimates biological processes; positive values indicate aerobic remineralization processes, which consume DO, and negative values indicate photosynthetic processes, which produce DO (Pytkowicx, 1971;Ito et al., 2004;Jackson et al., 2021).

Mixing Model
We use a simplified mixing model to determine the effects of mixing versus photosynthesis/respiration in the time series.
Because we have no nutrient measurements at the mooring, we use only salinity and temperature and thus can only use three end-members. The choice of the end-members was based on a T-S diagram (Supplementary Figure 3). ISSW 1 , IBSW 2 , and MSAAW 1 , define a triangle that best contains the mooring data and were selected to estimate the fractions of each type of water at the mooring at each moment. Table 1 shows the values of end-members of temperature and salinity for each water type. Equations (8a), (8b), and (8c) were used to determine the mixing ratios f. Equations (8d) and (8e) were used to calculate the DIC and A T that would result from mixing these three water masses.
where f is the end-member mixing ratio. T°C, S p , DIC, and A T are temperature, salinity, dissolved inorganic carbonate, and total alkalinity for each end-member ( Table 1). The pCO 2 due to mixing was calculated from the following equation: where the term f (DIC mix , A Tmix , S mix , T mix ) is the calculation of pCO 2 using DIC, A T , S p and T°C to the CO 2 SYS program (Pierrot et al., 2006).

1-D Biochemical Mass Balance Model
A simple 1-D mass budget model determined the contribution of biological processes, physical mixing, air-sea CO 2 exchange, and thermodynamics to pCO 2 variability. The changes in pCO 2 due to dissolution/formation of CaCO 3 and salinity are omitted. The daily changes in the pCO 2 were estimated as previously carried out by various authors (Chierici et al., 2006;Xue et al., 2016;Fransson et al., 2017;Li et al., 2018;Gac et al., 2021). Changes in pCO 2 (DpCO 2 ) are driven by changes in temperature (DpCO 2tem ), air-sea CO 2 exchange (DpCO 2gas ), mixing (DpCO 2mix ), and biological activity (DpCO 2bio ). DpCO 2 and DDIC are calculated using: where D is the change in parameters between times (begin) t n and t n+1 (end) are the difference between the daily value at t n+1 minus the daily value a t n .

The Effect of Thermal Changes
Temperature changes affect the dynamics of the pCO 2 but not the DIC of the water. Therefore, the thermal effect on DpCO 2tem was determined according to the equation (similar to Eqn 6 above; Takahashi et al., 1993) (11) where (pCO 2 ) n is pCO 2 at time 1, DT is the temperature difference between time t n+1 and t n .

The Effect of Air-Sea Gas Exchange
Air-sea CO 2 exchange affects the dynamics of DIC and pCO 2 but not the A T . Therefore, the pCO 2 changes due to the air-sea CO 2 exchange (DpCO 2gas ) are calculated using: where F gas is the air-sea CO 2 flux (see Section "Air-sea CO 2 exchanges") t n+1 -t n is the number of days between two measurements, r is the seawater density (kg m -3 ) calculated from the TEOS-10 calculations (McDougall et al., 2011), and d is the mixed layer depth taken to be 10 m, the depth of the mooring.
DIC n+1 ð Þ gas = DIC n + DDIC gas (12b) where (DIC n+1 ) gas is DIC predicted for time t n+1 based solely on air-sea CO 2 exchange. The term f ((DIC n+ 1 ) gas , A Tn , S n T n ) is the calculation of pCO 2 using DIC, A T , S, and T to the CO 2 SYS program (Pierrot et al., 2006).

The Effect of Water Mass Changes
To determine the changes in pCO 2 due to mixing the pCO 2 was calculated from DIC mix and A Tmix derived from the fraction of each water mass using equations: where DDIC mix and DA Tmix are DIC and A T of mixing (see Section "Mixing model") t n+1 -t n is the number of days between two measurements DIC mix(n+1) and A Tmix(n+1) are predicted DIC and A T based solely on the mixing.
is the calculation of pCO 2 using DIC, A T , S and T to the CO 2 SYS program (Pierrot et al., 2006).

The Effect of Biological Processes
The biological effect on pCO 2 is estimated as the remainder in the change of DIC.
where (DIC n+1 ) bio is the predicted DIC based solely on the biological factor. f ((DIC n+1 ) bio , A Tn , S n , T n ) is the calculation of pCO 2 using DIC, A T , S, and T in the CO 2 SYS program (Pierrot et al., 2006).

Hydrography and Nutrients
The hydrography of the Seno Ballena fjord was modulated by the seasonal input of freshwater from glacier meltwater, the glacial river located at the head of the fjord, and the small tributaries along the fjord. The influence of the MSAAW was also evident given the range of salinities observed (24 -31); the fjord water is estuarine-brackish (33 -66% seawater). The contribution of freshwater from the Inner Surface Estuarine Water (ISSW) near the glacier is distributed from the head of the fjord towards the outside of the sill near the surface, with smaller fractions outside of the sill ( Figures 3A-C). In fall and spring, the deep layer contained not only MSAAW, but also another type of water (IBSW) near the glacier, that is cold and nutrient-rich ( Figures 3D, F). In fall the surface ISSW was mixed deep inside the fjord. Finally, the MSAAW dominated the deep layer of the inner fjord, which confirms that this estuarine-brackish water enters the fjord ( Figures 3G-I). The MSAAW contribution was greater in winter (from 5 -50 m depth). This mixed water mass was associated with a strong vertical density gradient (s t range from surface to depth: fall = 19.20 -23.97; winter = 20.18 -23.98; spring = 17.17 -22.82). We divided the fjord area into three zones (inner sill, outer sill, and off-shelf), based on stratification and mixing conditions. In the inner sill region, the 0 -10 m surface layer showed low salinity, between S P = 27.4 to 28.6, over a deeper layer with S P ranging from 30.2 to 31.2 ( Figures 4A-C). Salinity varied also with the season, showing high salinity in winter (range from surface to depth: S P = 28.6 -30.3) and low salinity in spring (range from surface to depth S P = 27.4 -31.2); the latter modulated by the input of glacier meltwater (Figures 4-C).
The water column had a thermal inversion and weak temperature variation in the inner sill region during fall and winter, with an increase of 0.7°C from the surface to the deep layer ( Figures 4D, E). There were high temperature values in the fall (range from surface to depth 7.2 -7.7°C) and lowtemperature values in winter (range from surface to depth 6.4 -7.3°C).
The stratification was much shallower in fall and spring and practically disappeared in winter in the outer sill section (Figure 4), coincident with a more significant influence of the oceanic water and lower influence of freshwater (S P between 28.6 and 31.3 in the surface and deep layers, respectively; Figures 3, 4). Mixed  conditions were observed in the off-shelf region, with a high S P of 30.5 throughout the water column ( Figures 4A-C). The fall and winter data showed a very stable surface layer (>10 m) in the water column of the inner section, with Brunt-Väisälä frequency squared values of 10 -30x10 -3 s -2 compared to outer sections (< 10x10 -3 s -2 ; Supplementary Figure 4). Nutrient concentrations (SRP, NO − 3 , and DSi) had a clear spatial variability from the inner to the outer sections of the fjord ( Figure 5; Table 2). Surface NO − 3 showed low to medium concentrations along the fjord, ranging from 5.7 -7.4 mM (inner sill) to 9.3 mM (off-shelf), and in the deep layer NO − 3 concentrations ranged from 6.2 -11.1 mM (inner sill) to 11.5 mM (off-shelf) ( Figure 5). Most of the time DSi concentrations were below 5 mM throughout the water column along the fjord. DSi concentrations were very low in the surface water along the fjord in fall and spring, with a range from 0.2 mM (off-shelf) to 2.7 mM (inner sill), increasing up to 4.6 mM in the deep layer in the inner sill region. During winter, DSi values increased up to 5.5 mM, especially in the middle and inner sections of the fjord ( Figure 5, Table 2). SRP concentrations were more homogeneously distributed throughout the water column with a range from 0.6 to 1.3 mM, with the highest concentrations in the inner section below 10 m during spring and fall ( Figure 5; Table 2). NO − 3 was strongly correlated with SRP in fall (R 2 = 0.9; n=28; p< 0.001) and spring (R 2 = 0.6; n=25; p< 0.001). N:P fluctuated between 7.2 in fall and 6.9 in spring, with lower ratios than the expected Redfield N:P ratio (16:1). DSi:N ratios were low in fall (0.5), winter (0.1), and spring (0.2) during the hydrographic cruises. These values are consistent with the values of 0.2 -0.5 recorded in this fjord during spring by Torres et al. (2011a).

Carbonate Parameters From Hydrography Cruises
Low pCO 2 values were observed in the upper 10 m layer during the fall and spring cruises (< 400 matm) in the inner fjord section; whereas, higher mean values were seen during winter ( Table 2). Below the 10 m surface layer, pCO 2 values were much higher close to the glacier, 500 to 589 matm (Figures 6A-C; Table 2).
Mean pH values of 8.0 and 7.9 were observed in the inner section   Table 2). The lowest pH values were observed below 10 m in the inner sections of the fjord and for most seasons were coincident with the highest pCO 2 values (Figure 6). High surface Chl-a was observed during spring in the upper layer in the whole transect, ranging from 2.3 (inner sill section) to 1.7 mg m -3 (outer sill section) (Figure 4). These high Chl-a biomass was associated with high pH values (R 2 = 0.48; n=7; p<0.001) (Figures 4,  6, respectively).
The A T was low in the surface layer (< 10 m), coincident with low-salinity waters in the inner section, with seasonal mean values of 1736 to 1993 mmol kg -1 ( Table 2), associated with low values of W Ar in the surface (R 2 = 0.42 n=50; p<0.001) (Figures 6G-I). Values of A T were over 2000 mmol kg -1 in the inner section of the fjord below the 10 m layer. The maximum W Ar and A T values were recorded in the off-shelf sections (2.1 and 2081 mmol kg -1 , respectively). These results indicate that during all seasons, freshwater inputs from the melting glacier and the glacial river were low in alkalinity (5.82 mmol kg -1 ), whereas MSAAW was high in alkalinity (2081 mmol kg -1 ).

Parameters From Mooring Sensors
The sea surface temperature values (outer sill section) also showed low seasonal variability, with amplitude not greater than 2°C. High values were recorded in fall with a mean of 8.2 ± 0.32°C (range 7.8 -8.6°C) and low values were recorded in spring with a mean of 7.2 ± 0.22°C (range 6.6 -7.7°C) (Figure 7; Supplementary Table 3). The S P values recorded at the mooring showed a range from 29.8 to 30.5 with high values occurring in fall; whereas, minimum S P occurred in spring (Figure 7). The salinity was above 28 throughout most of the studied seasons but presented a clear low amplitude seasonality (S P <1.5) (Figure 7). The seasonal variability of T and S P is related to the variable contributions of the different water masses of the fjord, with minor contribution of ISSW (0.04 -0.23) compared to IBSW (0.05-0.97) and MSAAW (0.01 -0.85) at the mooring (Figure 8). MSAAW dominated the water masses during the fall; whereas, IBSW dominated during the spring. During the spring season, we observed at the end of October, an event during which the contribution of MSAAW (0.51) and IBSW (0.45) was more similar coinciding with a noticeable change in T and S P , with respect to the values observed when a larger fraction of IBSW or MSAAW is recorded.
The carbonate system showed clear seasonality in the upper 10 m layer, with the most noticeable changes during fall and spring ( Figure 7D). During fall, pCO 2 values had a mean of 385 ± 26 matm (range 240 -431 matm). pCO 2 values decreased in spring, with a mean of 365 ± 47 matm (range 167 -471 matm) and were associated with an increase in dissolved oxygen (negative values of AOU) and a decrease in salinity ( Figure 7B). An increase of pCO 2 was observed in winter, which reached a mean of 394 ± 10 matm (range 365 -433 matm), coincident with a decrease in temperature as the system changed from warming to cooling. A negative relationship was observed between pCO 2 normalized to an average temperature (npCO 2 ) and oxygen saturation percentage (%DO) (R 2 = 0.47, p<0.001; slope -3.7 matm %DO -1 ; Figure 9A). Although the npCO 2salinity relationship had high scatter ( Figure 9B), it showed a significant positive relationship during spring (R 2 = 0.30, p<0.001). The pCO 2 -temperature relationship ( Figure 9C) also showed a significant positive relationship during winter (R 2 = 0.41, p<0.001; Supplementary Table 2).
Air-sea CO 2 fluxes varied from -3.08 mmol m -2 d -1 in spring to 0.29 mmol m -2 d -1 in winter ( Figure 10A). CO 2 uptake in fall (March-May) fluctuated from -1.62 to -0.24 mmol m -2 d -1 , while in spring (September to December) it varied from -0.57 to -3.08 mmol m -2 d -1 . Degassing occurred only during the June-July winter months, with mean monthly fluxes of 0.09 and 0.26 mmol m -2 d -1 , respectively. The data show continuous CO 2 absorption over seasons, indicating that this fjord behaves as a net sink of atmospheric CO 2 (mean = -1.01 mmol m -2 d -1 ). Water mass changes also played an important role in the change of DIC (range 41.39 -180 mmol kg -1 ) and therefore in pCO 2 (153 -454 matm). We observed large changes during the spring related to an influx of IBSW. The differences in DIC and pCO 2 at the mooring compared to a mixture of the water masses give the impact of the other factors: biology, air-sea gas exchange, and temperature ( Figures 10A-C). The difference in DIC is always negative, which suggests that the mooring water mass always has more of a production signal, compared to a remineralization signal, than the underlying water types. These water masses are mainly subsurface water masses (ISBW and MSAAW) and have larger nutrient concentrations favoring primary productivity compared to the small contribution from the surface water mass (ISSW; Figure 8). The variability of pH T was inversely correlated with pCO 2 , showing maximum pH values during spring with a mean of 8.1 (range 8.0-8.3) and low values during winter and fall with a mean of 8.0 (range 8.0-8.1) ( Figure 7E). The seasonal trend of pCO 2 and pH is followed by W Ar (Figures 7D-F).

Hydrography and Nutrients
The results of our studies reveal the hydrographic control on the seasonal variability of the carbonate system, including the variability of the saturation state of seawater relative to calcium carbonate, pCO 2, and air-sea CO 2 fluxes (Figures 7, 10). Within the marine fjord system of Patagonia, we find fjords that are: (1) simple, two-layer vertical structures, with an estuarine surface layer of low salinity and a subsurface layer of oceanic origin with higher salinity (Schneider et al., 2014;Saldıás et al., 2019); and (2) more complex three-layer structures as found in the Seno Reloncavı́Fjord . The above is modulated by freshwater contributions of glacial origin, river discharge, and rainfall, which significantly influence the salinity gradient of the surface layer in the vast coastal ocean of western Patagonia (Acha et al., 2004;Saldıás et al., 2019).
Our results indicate that the Seno Ballena Fjord has two layers; the salinity of the surface layer is modulated by the input of meltwater from the glacier, mainly through stratifying the inner sill region of the fjord. The freshwater stratification and surface water salinity varied seasonally, with a relatively deeper halocline in fall and spring due to higher freshwater input (Figures 3, 4). The conditions were generally less stratified in winter, suggesting a limited supply of freshwater and strong ocean water contribution ( Figure 3). As was observed in the water mass distribution analysis and the salinity gradients along the Seno Ballena Fjord, the influence of freshwater was higher on the inner-sill surface and reduced toward the outer-sill section. Only in spring does the freshwater plume spread past the sill (Figure 3), as has been observed in previous studies in this fjord during spring (Torres et al., 2011a). The observed seasonal variability of temperature and salinity was weaker (range Temp: 7.2 -8.3°C/range S P : 29.8 -30.2) compared to other fjord systems of northern Patagonia (Seno Reloncavı́(41°S); range Temp: 10.6 -16.9°C/range S P : 20.1 -29.8) (Vergara-Jara et al., 2019) and western Canada (Strait of Georgia (49°N); S P range: 5 -29) (Moore-Maley et al., 2018), where strong seasonality in salinity is related to larger, seasonal river freshwater inputs (Puelo River: 713 to 4000 m 3 s −1 and Fraser River: 800 to 12.000 m 3 s −1 ).
The deep layer with higher salinity was influenced by the inflow of MSAAW: S P = 31-33, which is modified upon entry into the Strait of Magellan, mixing with freshwater from rivers and melting glaciers, forming estuarine water with salinity around 28 -31 S P (Sievers and Silva, 2008). During tidal flood, the acceleration of the flow forces aspiration (Bernoulli aspiration; Valle-Levinson et al., 2007), pumping denser deep water which sinks below the less salty buoyant water in the inner section . The high salinity-deeper water that flows inward compensates for the surface layer of freshwater that flows toward the Strait of Magellan. This type of circulation also influences the temperature of the fjord, causing an inverted thermocline with low surface temperature, underlain by a deep homogeneous layer of higher temperature. The low temperature in the surface layer of the Seno Ballena Fjord was more pronounced in winter due to the formation of a thin, surface layer of ice near the glacier, ice that likely is derived from re-frozen subglacial discharges rather than sea ice.
Our results in the Seno Ballena Fjord show that the nutrient concentrations were modulated inside and outside by the sill and revealed the significant spatial variability of nutrients (N and DSi), consistent with previous reports (Frangoṕulos et al., 2007). The presence of the shallow sill dividing the fjord causes a thin plume of low salinity water to spread over the sill during ebb tide (Torres et al., 2011a), which is probably modulated by the wind. The deepwater aspiration, when the water flow accelerates during the flood tide, probably influences the nutrient concentration patterns in the inner fjord ( Figure 11; Valle-Levinson et al., 2006;Torres et al., 2011a). There is a strong influence of MSAAW, which has low DSi compared to nitrate concentrations, which explains the higher nitrate concentration in the deep water. In contrast, the continental waters of western Patagonia are characterized by high DSi and low nitrate, which explains the higher concentration of DSi in surface waters and near the glacier (Hervéet al., 2007;Torres et al., 2020). Plumes, tides and wind stress have also been found to drive circulation and renewal of the water masses of Northern Hemisphere fjords with sills, demonstrating their importance in the physical, chemical, and biological processes that occur inside and outside the fjords (Mortensen et al., 2011).
The plume of low salinity is blocked to a greater extent by the sill, particularly during low tide (Valle-Levinson et al., 2006). This blockage was observed during fall and spring. Hence, we suggest that the inner fjord residence time seems to allow the accumulation of local remineralization products such as nutrients and CO 2 in the inner fjord, allowing the IBSW to form near the glacier during fall and spring. IBSW represents the water mass observed near the glacier at subsurface levels (30 -50 m), which has different properties (Figures 3, 5, 6) than those of the ISSW and MSAAW.

Carbonate Chemistry
The seasonal variability recorded in the Seno Ballena Fjord carbonate system can be assigned to several processes, including the input of freshwater from the melting of the Santa Ineś glacier, the exchange with oceanic waters, biological production, seasonal variability in temperature, water column mixing and advection. The surface layer of the inner-sill region had undersaturated pCO 2 values with respect to the atmosphere during fall and spring, associated with the input of freshwater from the glacier (Figures 3, 4, 6). In this area, a higher fraction of ISSW is observed (range 1 to 0.3 of ISSB fraction, the surface to 10 m, respectively), which coincides with the lowest DIC mix values (range 1489 -1852 mmol kg -1 of DIC mix , surface to 10 m, respectively; Supplementary Figure 5), this indicates that the DIC input of freshwater was low and contributes to the decrease of the DIC of the mixture and directly to the pCO 2 . These observations are consistent with previous studies in fjords and channels of Patagonia in which low-salinity surface waters (e.g. salinity less than 28) are strongly undersaturated in CO 2 (Torres et al., 2011b). Subglacial environments have low pCO 2 due to the effect of temperature (cooling), photosynthesis, and dilution (Bates and Mathis, 2009;Meire et al., 2015).
Although the pCO 2 was low in surface waters in the inner section of the fjord region during fall and spring, it did not result in the increase of calcium carbonate saturation (W CaCO3 : W Ar and W cal ) whose variation seems to be dominated by dilution (i.e. the reduction of CO −2 3 or Ca 2+ ions by dilution with freshwater with very low levels of these ions). Our results were similar to conditions found in Alaska and Arctic fjords, where glacier fjords had low pCO 2 with respect to the atmosphere, and corrosive levels of aragonite (W Ar ) and calcite (W cal ) saturation near the surface (<50 m) due to the input of freshwater with low A T and Ca 2+ from glacier melting (Evans et al., 2014). The above finding suggests that the inner surface water of the Seno Ballena Fjord was under-saturated and corrosive to aragonite during the survey period, due to the intrusion of low A T and Ca 2+ (relationship W Ar -A T ; R 2 = 0.42 n=50; p< 0.001) freshwater that drains the siliceous crystalline batholiths of Patagonia (Torres et al., 2011b;Torres et al., 2020). The above demonstrates the sensitivity of this Sub-Antarctic system to global stressors such as global warming and cooling, mainly for this Sub-Antarctic area where there has been a significant retreat of the glaciers in recent years (Bown et al., 2014).
The carbonate system presented clear seasonal variability in the outer sill region in March-December period, 2018 (Figure 7). A substantial reduction in pCO 2 and an increase in pH T in spring months can be attributed to the dominance of photosynthesis (AOU: -191 mmol kg -1 ) over respiration (%DO versus npCO 2 , Figure 9; AOU Figure 7C). After the spring bloom, a signal of remineralization with higher values of pCO 2 (471 matm) related to positive values of AOU (67 mmol kg -1 ; Figure 7C) is observed, these results illustrate the key role that phytoplankton blooms and microbial respiration play in the variability of pCO 2 , and therefore the change in pH and to a lesser extent the change in W Ar , taking into account the signal that observed during the spring season that coincides in the three variables (pCO 2 , pH T, and W Ar ). In contrast, high pCO 2 values were recorded in winter, related to low pH T and W Ar , which could be attributed to subsurface water masses brought to the surface by increased water column mixing, which is more consistent between the months of July and August (IBSW and MSSAW; Figure 8) as well as the predominance of respiration over photosynthesis.
The variations recorded in the pCO 2 time-series seem to be driven by net primary productivity events due to the consistency of pCO 2 and dissolved oxygen variability (Figure 7). However, it should be borne in mind that the dynamics of coastal systems are also driven by other processes such as the mixing of different water masses; air-sea fluxes of CO 2, and variations in temperature. Therefore, in this study, these effects were also calculated (DpCO 2mix , DpCO 2gas , and DpCO 2tem , respectively). During spring, the biological processes (-4.72 to -20.40 matm day -1 ), the changes in water mass (-5.42 to -17.2 matm day -1 ), and the sea-air CO 2 flux (-5.47 to -21.81 matm day -1 ) contributed significantly to the decrease of pCO 2. The temperature effect is much smaller (-0.10 to 0.29 matm day -1 ). Similarly in fall, DpCO 2bio (22.68 to 30.13 matm day -1 ), DpCO 2mix (10.8 to 30.13 matm day -1 ), DpCO 2gas (11.25 to 30.72 matm day -1 ) dominate, which is why there are higher values of pCO 2 in fall compared to spring. The DpCO 2tem contribution is again small (Supplementary Table 4), as is observed in the linear regression of pCO 2 /T ( Figure 9). In contrast, in other studies, T is the most important factor in the variability of pCO 2 , whereas in our studies it is the factor with the least importance in modulating the daily variability of pCO 2 (Xue et al., 2016).
The sea-air CO 2 flux varies seasonally, while low temperature and low haline stratification favors the positive sea-air CO 2 flux estimated in winter, the low wind speed limits those fluxes. In contrast, stronger winds during fall and spring contribute to the greater CO 2 uptake of Patagonian fjords. The mean wind speed was 4.04 m s -1 during the study period, with a maximum monthly average of 6.04 m s -1 (November) and a minimum of 3.73 m s -1 (July; see Supplementary Figure 6). The higher CO 2 uptake in spring is also related to the decrease of pCO 2 in surface water by photosynthesis. Previous studies have also suggested that the fjords of Patagonia in warm seasons behave as a sink for CO 2 in response to high productivity (Torres et al., 2011b;Vergara-Jara et al., 2019). This behavior is also observed in fjords of the Northern Hemisphere (Meire et al., 2015).
Variation in pCO 2 is also driven by water mass changes. In fall, MASSW dominates with high values of DIC, whereas in spring, ISBW dominates with much lower DIC values. When these subsurface water masses are mixed to the surface, they bring macronutrients, products of remineralization, to the surface layer and favor primary productivity (Jones et al., 2020). The nutrient ratios (Si:N<1 and N:P<15) may indicate that the availability of silicic acid and nitrate plays a critical role in determining the extent and/or the persistence of phytoplankton blooms and the effects on the variability of the carbonate chemistry of waters in the study area, mainly in the spring season.

CONCLUSION
Frequent pCO 2 and pH T measurements allowed us to determine that the carbonate system's seasonal variability in the Seno Ballena Fjord was weaker compared to the northern Patagonia fjords. The most notable changes in pCO 2 , pH T , and W Ar were driven by spring primary productivity, likely induced by the effect of seasonal increase of solar insolation in a recently stratified layer and local nutrient fertilization processes by the subsurface water mass ISBW (through aspiration and entrainment). In our results, we also found that the changes of the water masses play an important role in the daily variability of pCO 2 throughout the time series. In fall, MSAAW dominates with greater DIC contributing to the mixture and higher values of pCO 2 are recorded, whereas in spring the IBSW dominates with lower DIC compared to MSAAW, and lower pCO 2 values are recorded. The important role of mixing of the water masses ISSW, ISBW, and MSAAW and the air-sea exchange of CO 2 in the daily variability of pCO 2 was shown, which was generally just as important as biological processes. The very low freshwater-A T end member means that meltwater significantly dilutes the concentration of A T /DIC and calcium ions and plays a dominant role in driving low W Ar values in estuarine waters ( Figure 11). This study points out the importance of subglacial discharge from Santa Ineś glacier in influencing the hydrography, nutrient, and carbonates system dynamics of the Seno Ballena Fjord, illustrating the sensitivity of this Sub-Antarctic system to global stressors such as global warming and freshening. The results highlight the important role of this high-latitude fjord as a CO 2 sink and the main processes driving this trend, including wind speed, low temperatures, and low salinity, these last two driven by surface freshwater runoff. These results suggest the importance of sustained high-resolution and long-term measurements of different physical and biogeochemical variables to understand the complex dynamics of the subantarctic fjords.

DATA AVAILABILITY STATEMENT
The data of the sensor SAMI-CO 2 , SAMI-pH, SBE37, and Aanderaa 5331A used for this study are publicly available and can be accessed at https://figshare.com/articles/dataset/ Buoy_Seno_Ballena/12888008

AUTHOR CONTRIBUTIONS
JV was responsible for data processing and analysis, designed the manuscript, and wrote the manuscript. GS and SA. Participated FIGURE 11 | Schematic of processes that occur inside and outside of the Seno Ballena glacial fjord: (1) The input of freshwater from the marine-terminating glacier and river glacier may change physical-chemical properties, such as the following: i) generate strong surface stratification in the inner fjord section, ii) bring low A T , low W AR , low nutrients (SPR and NO 3 ), and low pCO 2 , and iii) bring sediment and dissolved organic carbon (DOC). (2) Advection of subsurface nutrient-rich waters near the glacier. (3) Bernoulli aspiration caused by the fast shallow flow into the inner fjord is capable of sucking deep water coming from the Strait of Magellan up and over the sill. The density of this water mass then causes it to sink.