Assessment of Annual Physico-Chemical Variability via High-Temporal Resolution Monitoring in an Antarctic Shallow Coastal Site (Terra Nova Bay, Ross Sea)

: The Southern Ocean is an important atmospheric carbon sink, and potential changes in the carbon ﬂux in this region will affect the ocean as a whole. Thus, to monitor the variability of its physico-chemical parameters is becoming a priority. This study provides the ﬁrst high-resolution all-year-round record of observed and computed physico-chemical data from a shallow coastal site in Terra Nova Bay (Ross Sea). From November 2018 to November 2019, an underwater observatory deployed at a 25 m depth under an ice pack recorded pressure ( p ), temperature ( t ), electrical conductivity ( C ), dissolved oxygen ( DO ), pH in total scale (pH T ), and illuminance ( E v ). Practical salinity ( S P ), density ( ρ ), tidal constituents, carbonate system parameters (total alkalinity (TA), carbon dioxide partial pressure ( p CO 2 ), calcite, and aragonite ( Ω Ca , Ω Ar )), together with sea ice concentration ( SIC ) and chlorophyll-a ( Chl -a), were derived from measured and satellite data. t , DO , and pH T displayed the lowest values between July and November (–1.95 ◦ C, 6.61 mL L − 1 , 7.97) whereas the highest in January (+1.08 ◦ C, 10.61 mL L − 1 , 8.35). S P had the lowest values (33.72 PSU) in February and the highest (34.87 PSU) in September. E v peaked in March (201 lux), with the highest values (>50 lux) in correspondence to the lowest values of SIC and a delayed trend, between December and March, with respect to Chl -a values (0.2–1.1 mg m − 3 ). Ω Ca and Ω Ar showed their highest average monthly values ( ± s.d.) in January ( Ω Ca : 3.41 ± 0.27; Ω Ar : 2.14 ± 0.17), when DO had maximum values. The lowest Ω occurred in September ( Ω Ca : 2.11 ± 0.02; Ω Ar : 1.32 ± 0.02), at the end of phytoplankton activity. No undersaturation for both calcite and aragonite was recorded during the study period. This study highlights that biological activities and physico-chemical variables of the investigated shallow coastal site are coupled and, in many cases, inﬂuence each other.


Introduction
Antarctica is a key component of Earth's climate system, interacting with the rest of the planet through shared ocean, atmosphere, and ecological systems. Its global footprint is valuable, considering the ocean surface occupied, as well as for its extensive permafrost and glacier areas [1,2]. Despite the continuous increase in atmospheric carbon dioxide (CO 2 ) concentration, Antarctica and the Southern Ocean have experienced less uniform temperature changes and a general slowdown in the rate of warming across the past 30-50 years [3,4]. This might be due to a heat redistribution within the atmospherecryosphere system [5], with the reduction in atmospheric heating almost equating in terms of energy to the contemporaneous increases in ice melting. Over the next century, the effects of climate change on the continent and in the ocean are expected to exacerbate [6]. In the marine realm, ecosystems are expected to be affected by altered environmental factors, which are not operating in isolation but rather synergistically, either additively or antagonistically. All these factors have potentially large impacts on coastal environments and comprise increased temperature and altered sea ice coverage, iceberg scour in benthic habitats, ocean acidification, salinity/freshening, and low oxygen levels [7][8][9][10]. Considering the global importance of the Southern Ocean as an atmospheric carbon sink [1], any changes in its ecosystem characteristics will determine the future states not only of the Southern Ocean and its life but also of the oceans as a whole [11]. Thus, there is an urgent need to implement data acquisition in the marine realm, both in time and space, for understanding and monitoring the marine systems [12].
Current rapid technological advances in developing modern sensors and probes for autonomous physico-chemical monitoring have enabled marine scientists to improve the exploration of natural environmental variations in deep and coastal ocean systems [13]. For example, monitoring the annual variability of ocean physico-chemical parameters is especially important since marine systems have defined seasonal cycles that change over the year [14]. However, due to logistic constraints of the extremely adverse marine environment, the continuity of data acquisition and production of big data are still major issues. This is especially relevant for remote Antarctic coastal sites, where conditions are unpredictable, owing to seasonal or year-round ice cover, iceberg scouring, and extreme storms.
The Ross Sea is one of the most biologically productive regions of the Southern Ocean [15], which determines global biogeochemical cycles and the sequestration of anthropogenic CO 2 into the deep ocean [16,17]. In addition, the formation of the Antarctic bottom water (AABW), which originates from the dense shelf waters (DSW) that escaped from the continental shelf and mixed with ambient warm water, influences the climate on regional and global scales [18,19]. The variability of surface carbonate system properties, primarily controlled by biological activities during summer months, can be also influenced by the Southern Annular Mode and extreme El Niño oscillations [20]. These factors have been responsible for changes in the salinity, density, and thickness of the AABW due to increased sea ice formation on the continental shelf triggered by anomalous wind forces [20]. Given the complexity and importance of this highly productive system, the implementation of data through the deployment of in situ monitoring stations to record fast occurring variations under climate change pressures is becoming a priority.
Between November 2018 and November 2019, an underwater observatory was established at 25 m of depth under the ice pack in Terra Nova Bay (TNB, western Ross Sea) with the aim to collect in situ 12-month data of physico-chemical parameters. This project was realized during XXXIV and XXXV Italian Expeditions to Antarctica and had the overall aim of directly correlating the physico-chemical characteristics of the seawater with skeletal physiological features of marine calcifying ecosystems to validate the organism's role as a proxy of environmental conditions. The functions of these ecosystems, especially in polar areas, are extremely relevant even though still poorly known. Calcifying filter feeders are major players as carbon sinks [21,22] by eating phytoplankton and calcifying skeletal structures (i.e., growth) [23][24][25]. This high-temporal-resolution monitoring allowed physico-chemical processes to be described that affected biological communities at a coastal site, including these ecosystems. In this paper, we present and discuss the first high-resolution all-year-round record of observed and derived physical-chemical data set resulting from the Ross Sea. In detail, the goals of the investigation were: 1.
to describe the variability of physico-chemical parameters, both measured (pressure (p), temperature (t), electrical conductivity (C), dissolved oxygen (DO), total hydrogen ion scale concentration (pH T ), and illuminance (E v )) and derived (practical salinity (S P ), density (ρ), and carbonate parameters), as a result of an in situ high-temporalresolution monitoring;

2.
to interpret the measured physical-chemical parameters as a function of sea ice concentration (SIC) and chlorophyll-a (Chl-a) mass concentration obtained from satellite data; 3.
to estimate the sensitivity of the carbonate system (total alkalinity (TA), partial pressure of carbon dioxide in seawater (pCO 2 ), and calcite and aragonite saturation states (Ω Ca , Ω Ar )), derived by direct measurements of physical and biogeochemical drivers at the coastal site.

Study Area: Location and Characteristics
The study was carried out in a coastal site of TNB, which is located on the western margin of the Ross Sea ( Figure 1 [26] in the south. The bay comprises a tortuous continental shelf with numerous banks and trenches. The mean depth of the shelf is about 450 m, with greatest depths close to the coast and up to 1000 m depths in the adjacent basin. The orography of the region is responsible for determining its climate. Local katabatic winds often keep the area free of ice, causing temperature drops (surface seawater temperature range: −1.9, +1 • C) [27] and a delay in summer seawater stratification, and are also responsible for the formation and preservation of polynyas persisting during wintertime [28]. TNB is one of the source areas for DSW, the precursor for the AABW, which is produced on the Antarctic continental shelf, where winter cooling and salinification by sea ice formation produce shelf waters of sufficient density to sink to the deep ocean [20]. The strongest heat loss and salt injection occur in coastal polynyas, where persistent wind-driven export of sea ice allows for continuous sea ice formation [29]. Sea ice formation in the TNB polynya represents the local inputs contributing to the observed high salinity (absolute salinity values greater than 34.9 g kg −1 ) of the DSW on the western Ross Sea continental shelf [19].
The seasonal sea ice coverage in TNB strongly impacts the cycling of chemical species and biological processes by restricting the CO 2 exchange between the air and the sea during winter and thus limits the equilibration of the Antarctic surface water with the overlying atmosphere. As a result of the heterotrophic activities, this surface water becomes supersaturated with CO 2 and undersaturated in terms of DO [30]. The carbonate system in the Ross Sea and the variability associated with different water masses were investigated in a summer campaign in 2008 [31], showing that the highest pH T values in TNB (about 8.03 at 200 m of depth) are measured where waters are oversaturated (more than 1) in terms of both calcite and aragonite [32,33].
TNB is entirely covered by sea ice for at least 9 months of the year [26]. When the ice melts, phytoplankton blooms occur, which lead to big changes in the carbonate systems, also affecting local productivity of benthic ecosystems. The seasonality of phytoplankton production is generally characterized by a main peak occurring from the second half of December through to the beginning of January and eventually by a second late-summer bloom in February [34]. Phytoplankton bloom normally occurs when maximum irradiance is reached and there is enough energy to cause the ice to melt as well as to warm up and stabilize the upper layer of seawater. Irradiance values range from a maximum of 36 MJ m −2 day −l (equivalent to 70 E m −2 day −l of PAR) near the solstice of December to a minimum below 5 MJ m −2 day −l in the second half of February, in clear-sky conditions [34]. Under such conditions, thermal stratification occurs and nutrient depletion is very pronounced, with nutrients sometimes reaching concentration levels that might be limiting for phytoplankton growth [35]. Moreover, physical factors, such as the strong winds, play an important role in developing a phytoplankton system that, for TNB, shows a patchy distribution, causing a shift from an autotrophic community in areas close to the glaciers toward a heterotrophic community in the shallower layers of free ice areas [36].

The Underwater Observatory
The underwater observatory ( For the environmental monitoring, three types of probes were used: 1-one pH meter (SeaFET V2, SEA-BIRD, Bellevue, WA, USA) and a solid-state (ion-sensitive field-effect transistor (ISFET)) probe for long-lasting time series of pH T measurements; 2-one conductivity-temperature-depth probe (CTD, SBE37-SMP-ODO, SEA-BIRD, Bellevue, WA, USA) for direct measurement of p, t, C (from which indirect measures of S P and ρ were derived), and DO; 3-six light and temperature sensors (HOBO UA-002-64, Pendant, Bourne, MA, USA), essentially dedicated to Ev measurement.
The pH meter and CTD probes were mounted directly on the metal frame and then hydraulically connected by a pumped duct in order to simultaneously analyze the same portion of water; light and temperature loggers were mounted on the top of the cages containing the biological samples, not analyzed in this study ( Figure 2). Prior to deployment, a preliminary testing activity in the laboratory (operational and mechanical checks) was performed on all probes and sensors. General metrological features of deployed probes are shown in Table 1, whereas values of accuracy and stability are reported as declared on the company datasheet. CTD and pH meters were both produced in 2018. Thus, for year-round acquisition, their values of accuracy and stability were assumed to fall within the validity of their calibration documents. The extreme environmental conditions at the study site led us to assess the stability of these instruments by using ad hoc methods: when possible, their stability was verified to the best of the available experimental capabilities (i.e., by applying pre-and post-deployment metrological checks performed in the laboratory; Figures S1-S3, Supplementary Materials); otherwise the calibration documents were used. It should be noted that the information provided by these documents is not in accordance with the current rules and standards governing the expression of measurement uncertainties [37] and its reliability, as a combined standard uncertainty associated with the quantity of interest, is partially neglected because of the lack of the associated coverage factor to the accuracy or error. Given that the accuracy of SeaFET has to be intended as a standard uncertainty (company's private communication to G. Raiteri and referring to previous in-field metrological characterization of a comparable CTD probe [38], the combined standard uncertainties u c of the main quantities of interest are reported in Table 1 (values should be considered as conservative, including the time stability contribution calculated at the end of a 12-month deployment). Before being deployed, all probes were synchronized to UTC time and programmed to acquire one measure per hour (at the beginning of each hour). Simultaneous data acquisition started on 25 November 2018 (04:00 UTC) and ended 361.7 days later, on 21 November 2019 (21:00 UTC), allowing a continuous hourly time series of 8682 recordings to be logged for each parameter measured by CTD and SeaFET probes (light-intensity measurements lasted about 4 days less, due to the early recovery of the cages housing the HOBO sensors). All the acquired data satisfied quality control by the instruments employed. For pH T measurements, 2 weeks prior to the deployment in the laboratory, the SeaFET measuring cell was filled with seawater samples collected from the field and filtered for preliminary checks to be carried out. Once the probes were deployed, data acquisition started after the recommended conditioning period (24 h) needed by the instrument for recording optimal data (see [39] and Supplementary Materials for further details on probe preparation).
Once recovered, all probes were inspected to check that no hardware anomalies had occurred during data acquisition. Then, data were downloaded to check the presence of spikes, anomalous noise, or trends before being analyzed. For pH data, the pH T signal measured by one of the two different SeaFET reference electrodes (the internal (INT) one) showed very low and spiky values, so failing quality control analyses. This was probably due to the potassium chloride electrolyte escaping from the INT reference, causing a failure in the pH T internal output (company personal comm. to GR). Thus, data from the SeaFET_INT sensor are not reported in this paper. The SeaFET external (EXT) sensor, however, was completely unaffected by the INT sensor behavior and acquired pH T values free of anomalies. To provide the most accurate and stable pH T data [39], the EXT sensor data were post-processed by applying a salinity-temperature correction based on CTD data through the SBE-UCI software (v. 2.0.2) CTD data were processed by using company software (SBE-Seaterm V2, v.2-7-0-108). Relationships and charts were realized by using Ocean Data View [40]. In particular, for pressure measures, fast Fourier transform (FFT) analysis was applied to calculate principal tidal constituents (principal solar semidiurnal (S 2 ), principal lunar semidiurnal (M 2 ), larger lunar elliptic semidiurnal (N 2 ), principal lunisolar diurnal (K 1 ), principal solar diurnal (P 1 ), principal lunar diurnal (O 1 ), larger lunar elliptic (Q 1 )) and their relative periods.
TA was estimated by applying the formula from [41] based on measured sea surface temperature and practical salinity. pCO 2 , Ω Ca , and Ω Ar were calculated by TA and pH T using CO2Calc software (v. 4.0.9, U.S. Geological Survey [42]; see Supplementary Materials for methodologies). Because of the lack of proper analytical instruments at the Italian base and the very short time spent at MZS (2 months per year) during the whole study period (12 months), spectrophotometric measurements for the pH and Winkler for dissolved oxygen were not used to validate data from the sensors of pH and DO.

Satellite Data
Data of sea ice concentration (SIC) were downloaded from the Copernicus Marine Environment Monitoring Service data base (CMEMS, https://marine.copernicus.eu) to complement in situ physico-chemical parameters (Figure 3). SIC is reported as a fraction from 0 to 1, where 0 corresponds to the absence of ice and 1 is the maximum ice coverage for the region. The CMEMS product used for the analyses had the following features: product level: L4, i.e., a product for which a temporal averaging method or an interpolation procedure is applied to fill out missing data values. Temporal averaging is performed on a monthly basis. Satellite chlorophyll-a data were also obtained from the CMEMS database to complement observed parameters (DO and E V ). Data are expressed as mass concentration of Chl-a in seawater (mg m −3 ) (Figure 3). Analogously, the CMEMS product used for the analyses had the following features: 1.
type of data: daily observations of mass concentration of Chl-a in seawater (mg m −3 ); daily mean reprocessed from 1997; 3. spatial resolution of data: 4 km × 4 km (surface only); 4. product level: L3, i.e., data are the daily composite products as obtained by merging all the ocean satellite passages.

Physical and Biological Environment
Monthly and annual values for measured (pH T , p, t, S P , ρ, DO, E V ) variables for the period November 2018-November 2019 are displayed in Table 2. Table 2. Mean annual and monthly values, with their associated standard deviation (s.d.) (i.e., the measure of dispersion of the values over time), of total scale pH (pH T ), temperature (t), practical salinity (S P ), density (ρ), dissolved oxygen (DO), illuminance (E V ), total alkalinity (TA), pCO 2 , Ω Ca , and Ω Ar at the study site between November 2018 and November 2019. Values indicated by * and ** were calculated using CO2Calc from temperature and salinity measures [42] K1 and K2 dissociation constants estimated by [43], total boron from [44], and KHSO4 dissociation constant by [45]. Regarding t (Figure 4a), values showed an annual variation of 3.03 • C, with a mean of −1.63 ± 0.56 • C. Warmer temperatures occurred in January, with variations of −0.13 ± 0.37 • C and values above 0 • C (up to 1.08 • C) during part of the day when the maximum solar irradiance was reached, between 8 a.m. and 2 p.m. On the other hand, the coldest months were July and August, with a very stable temperature of −1.93 • C (minimum equal to −1.95 • C).  Regarding relationships among t, S P , and density (expressed as sigma-t density σ t , i.e., σ t − 1000 kg m −3 , data showed that, in the study site, seawater density was mainly driven by practical salinity rather than temperature differences (for t = −1.   (Table 2).
Finally, no significative relationship was found between pH T and tidal height for the study period.

Calcification Environment
To explore how the calcification environment varied over time in the shallow coastal site of TNB and to evaluate the influence of biological activities, carbon parameters (TA, pCO 2 , Ω Ca , Ω Ar ) were calculated for the period of investigation. In addition, to validate the data, our estimations from the study site were plotted with data from nearby coastal locations in McMurdo Sound (Jetty and Cape Evans) collected in 2011-2012 and 2012-2013, respectively. All values are reported in Table 2 and plotted in Figure 9. The trend of pCO 2 was comparable to that of TA (Figure 9a), with minimum peaks recorded in February 2019 (pCO 2 : 240 ± 21 µatm; TA: 2305 ± 5 µmol kg SW −1 ), but was in contrast to the measured pH T trend (max pH T : 8.24 ± 0.03; Figure 4a). Monthly mean calculated pCO 2 values ranged between 240 and 435 µatm. The highest values (greater than 400 µatm) were recorded during 42% of the study period (July-November 2019), whereas values in the range of 300-400 µatm and 200-300 µatm were recorded during 33.3% and the remaining during 25.1% of the study period (January-March 2019).
Saturation states of both calcite and aragonite showed their highest values in January 2019 (Ω Ca : 3.41 ± 0.27; Ω Ar : 2.14 ± 0.17), in correspondence to the maximum value of DO (8.85 ± 0.549 mL L −1 ), and the minimum values in September 2019 (Ω Ca : 2.11 ± 0.02; Ω Ar : 1.32 ± 0.02), at the end of phytoplankton activity (Figure 9b). Both Ω Ca and Ω Ar showed reasonable annual trends, with a significative decline from summer to autumn/winter. Experimental results are also qualitatively comparable, with measurements performed at McMurdo, especially with those during 2011-2012. No undersaturation estimates were recorded for the 12-month period for both calcite and aragonite saturation states.

Discussion
Mooring sites are well spread in Antarctica, recording physico-chemical seawater data in deep sites as well as describing Antarctic current dynamics [47][48][49] of sites, including Ross Sea and Terra Nova Bay [50][51][52]. With the exception of a few sites, such as Arthur Harbor (Palmer Station Antarctic Peninsula) [53], King George Island [54], and McMurdo Sound (western Ross Sea) [14,46,55], to date, there are not many underwater annual incontinuum measurements in coastal shallow sites, representing very unique and dynamic environments.
Prohibitive and unpredictable environmental conditions put instruments at risk, but such observatories are extremely important because of the need for data acquisition in polar areas due to the rapidity of the occurring changes [1,2]. The data presented here are from the underwater observatory (Figure 2), anchored on the sea bottom at 25 m of depth under the ice pack, which allowed in situ seawater parameters to be recorded over 12 months (Figure 4, Table 2). Despite the limitation of this study to a single station, the first successful all-year-round underwater monitoring allowed high-temporal-resolution data to be produced describing physico-chemical and biological processes occurring in a shallow coastal site at TNB. The area has been studied for the past 30 years by Italian Antarctic researchers [26,56], but in situ 12-month physico-chemical variability has never been previously monitored.
Our findings reveal that, between November 2018 and November 2019, the trends of physico-chemical parameters reflected the seasonality of the environment. Major changes occurred between December and March, probably driven by two main processes: ice melting and biological activities.
The temperature remained stable at −1.9 • C for 7 months, between April and November, and displayed a variability of 3.03 • C across a short time frame between December and March. Temperature represents one of the major drivers of biological processes by inducing the activation of physiological activities in most of the local ecosystems [57,58]. Furthermore, the temperature increase, determined by the sunray angles, represented the trigger for sea ice melting, even though it did not represent the major driver for density, confirming an important role of the polynyas and dense shelf waters (DSW) originating in coastal areas and detectable on a very local scale [20].
Density was mainly driven by practical salinity (i.e., lowest salinity = lowest seawater density), reaching the minimum value of 33.72 PSU in February 2019, in correspondence to minimum SIC (sea ice area fraction equal to 0) and low temperatures (−1.35 • C). Interestingly, the most saline waters (S P greater than 34.8 PSU) occurred between September and November 2019, when the sea ice coverage was the highest (SIC greater than 0.8).
Given the position of the study site (the observatory is located at 25 m of depth), this saline water can be attributed to brine rejection from sea ice formation in the TNB polynya and also to remote inputs, as salt content might be advected toward TNB by coastal currents (DSW) [20]. These increases in salinity and density, associated with the sea ice coverage, also drive convection and mixing and influence sea ice motion.
In fact, the sea ice motion was clearly visible between December 2018 and February 2019, with SIC values oscillating from 0 to a maximum peak of 0.4 and back to 0 again, as a result of a combination of stresses from the atmosphere and the ocean [59,60]. This motion was determined by forces applied both on the ice surface by katabatic wind stress as well as at the base by ocean currents, including tides [61]. Gradients in tidal currents, causing periodic divergence and convergence of the ice pack [62], were responsible in determining the biggest SIC range (0-0.4) during summer (January 2019). Noticeably, when local tidal effects on sea ice are coupled with the seasonal cycle (autumn-winter net ice production and spring-summer net ice loss) and with the larger-scale advection of water masses (determined by the DSW), the turnover between sea ice production and loss has a substantial effect on the large-scale ocean and sea ice state [63]. Tide-induced signals provide insight into the processes by which the oceans can affect ice sheet mass balance and dynamics [64]; thus high-temporal-resolution data from underwater observatories represent an important tool for improving understanding of the role of tides in ice sheet dynamics and for predicting the future ice sheet mass budget and sea-level rise.
The seasonal sea ice coverage is accompanied by the rise or decline in biological activities [6,31] which, also for the site, were intensified at ice melting [65]. In December, the sea ice melting exerted a possible influence on controlling the carbonate system chemistry of the site both directly through dilution processes and indirectly by favoring the development of phytoplankton blooms, as observed in McMurdo Sound [46]. Differently from other studies considering mixed layers [66], at the shallow coastal site in TNB, the increase in pH T recorded at the end of the year (the maximum peak during this period with a value of 8.35) was accompaigned by maximum values for both estimated saturation states (Ω Ca = 4.42; Ω Ar = 2.78). While pCO 2 displayed the minimum value, with 176 µatm, and a dilution of TA was recorded. In accordance with processes observed in McMurdo Sound (western Ross Sea), the variation in the thickness and attenuation of sea ice, starting in December, might have doubled the effect on the measured pH T due to photosynthesis [14,52]. This would explain the highest pH T value recorded at the beginning of the season. In the meantime, the DO peak (oversaturation: 133%) occurring in December was due to the first phytoplankton bloom (i.e., sympagic algae) of the season, with a Chl-a mass concentration of 0.3 mg m −3 and the start of its photosynthetic activity. In fact, sympagic algae are released in small quantities from sea ice before its breakup and, in a short time frame, as a strong bloom during and soon after the melting of sea ice [67][68][69]. The photosynthetic performance of phytoplankton reveals a general adaptation to low-light regimes, characterized by high values of photosynthetic efficiency and also by low-temperature adaptation of the growth performance [70].
The greatest change in the carbonate system occurred in January, confirming the effect of melting freshwater (SIC range: 0-0.1), which also affected TA, resulting in a minimum of 2292 µmol kg SW −1 [66,71]. pH T showed a second positive peak with 8.33, and lowestimated pCO 2 concentrations of approx. 200 µatm were observed. Probably due to the highest quality of light available at 25 m (approx. 50 lux), because of the free ice coverage (SIC = 0), Chl-a showed the maximum mass concentration (approx. 1 mg m −3 ). Thus, the maximum photosynthetic activities would be expected, but, interestingly, the DO trend declined (mean saturation: 78%) as an indication of a post-bloom effect. Both derived saturation states also declined from December to February, stabilizing at minimum constant values (greater than 1 mL L −1 ) in July 2019. During summer months, the variability of the surface carbonate system properties at this shallow site was primarily controlled by biological activities [20]. These include not only phytoplankton blooms, changes in species composition depending on sea ice influence, current, and water mixing [68], but also physiological processes of bacteria [72] and, on a macro scale, of benthic ecosystems (i.e., respiration, photosynthesis, calcification). Benthic ecosystems of the site were primarily characterized by calcifying filter feeders [65,73], sessile or vagile, such as bryozoans, corals, sponges, molluscs, and echinoderms. After the ice melt, an abundant quantity of food becomes available in a few hours to these consumers, which adopt a more specialized diet and increased feeding rate [74,75]. This was observed both at the level of the entire food web and in species analyzed from the same site during the two-year investigation [63]. Thus, considering the abundance of califiers at the site, DO, pCO 2 , and saturation state reduction can be due to intensification of their physiological activities, such as growth (i.e., calcification causing carbonate consumption) other than reproductive activities during summer months [58,76].
Under the rapidly progressing climate change, the complexity of processes and their variability occurring at local scales in a very short time frame call for more data collection. Even the local character of our study, TNB, is where the DSW, the precursor for the AABW, orginate and bring information on water masses from local coastal systems [20]. The increase in anthropogenic CO 2 storage in the AABW recorded in the past 40 years was accompained by a long-term change in total carbon concentration due to anthropogenic CO 2 uptake of the different formation regions [33]. All these changes were modulated by significant interannual to multi-annual variability associated with variations in physicochemical properties (i.e., temperature, practical salinity, DO, TA, pCO 2 ) orginating at local other than regional scales involving the ecosystems. Thus, we hope in the increase of in situ observatories, representing powerful tools for recording and monitoring rapidly occurring changes in shallow and deep polar environments [2,77] as well as their cascading effects on polar biodiversity [9,58,78].