Ross Sea Dissolved Organic Matter Optical Properties During an Austral Summer: Biophysical Influences

The Ross Sea, one of the most productive regions in the Southern Ocean, plays a significant role in deep water formation and carbon cycling. Dissolved organic carbon (DOC) concentrations and chromophoric dissolved organic matter (CDOM) absorption and fluorescence (FDOM) properties were studied in conjunction with biophysical properties during austral summer. Elevated values of both DOC (mean 47.82 ± 5.70 μM) and CDOM (absorption coefficient at 325 nm, acdom325: mean 0.31 ± 0.18 m–1) observed in the upper shelf waters in the southwest (SW), north of the Ross Ice Shelf (RIS), the northwest and along a transect inward of the shelf break, suggested in situ production and accumulation linked to the productive spring/summer season. However, regional differences were observed in CDOM with acdom325 higher (0.63 ± 0.19 m–1) and its spectral slope S275–295 lower (24.06 ± 2.93 μm–1) in the SW compared to other regions (0.25 ± 0.08 m–1 and 28.92 ± 2.67 μm–1, respectively). Similarly, the specific UV absorption coefficient or SUVA254 determined at 254 nm was greater (1.85 ± 0.55 m2 mg–1 C) compared to other regions (1.07 ± 0.24 m2 mg–1 C), indicating CDOM of greater molecular weight and aromaticity in the SW. Phytoplankton absorption spectra indicated the shallow mixed layer of SW Ross Sea to be dominated by diatoms (e.g., Fragilariopsis spp.), a preferential food source for grazers such as the Antarctic krill, which in large numbers have been shown to enhance CDOM absorption, a likely source in the SW. Excitation-emission matrix (EEM) fluorescence combined with parallel factor analysis (PARAFAC) retrieved one protein-like and two humic-like FDOM fractions commonly observed in the global ocean. In contrast to acdom325 which was uncorrelated to DOC, we observed weak but significant positive correlations between the humic-like FDOM with salinity and DOC, high value of the biological index parameter BIX and an instance of increasing FDOM with depth at a location with sinking organic matter, suggesting autochthonous production of FDOM. The absorption budget showed a relatively higher contribution by CDOM (70.7 ± 18.3%) compared to phytoplankton (22.5 ± 15.2%) absorption coefficients at 443 nm with implications to ocean color remote sensing. This first study of DOM optical properties provides additional insights on carbon cycling in the Ross Sea.


INTRODUCTION
The Ross Sea, a highly productive region of the Southern Ocean, accounts for 25-30% of the annual Southern Ocean primary production and ∼25% of bottom waters formed globally (Orsi et al., 2002;Arrigo et al., 2008;Smith et al., 2014), thus playing an important role in the marine carbon cycle (DiTullio and Arrigo et al., 1999Arrigo et al., , 2008. Located north of the largest ice shelf in the world, the Ross Ice Shelf (RIS), is the largest polynya (a region of open water surrounded by sea ice), the Ross Sea polynya, associated with the formation of the high salinity and cold shelf water or the dense shelf water (DSW) owing to extensive brine release in the polynya during winter (Jacobs et al., 1970;Jacobs and Giulivi, 1998). The DSW sinks and eventually flows off the shelf contributing to the Antarctic Bottom Water (AABW), deep ocean circulation and carbon cycling (Budillon et al., 2003;Orsi and Wiederwohl, 2009;Bercovici et al., 2017). During austral spring and summer as the Ross Sea polynya expands due to increasing solar insolation and warming, extensive algal blooms occur that are dominated by the haptophyte Phaeocystis antarctica followed by increased phytoplankton growth of diatoms, mainly the Fragilariopsis species (Asper and Smith, 1999;Arrigo et al., 1999;Smith et al., 2014). Many insights into these bloom dynamics, which are often spatially and temporally distinct, have been obtained from numerous field (DiTullio and Smith et al., 1996;DiTullio et al., 2000;Sweeney et al., 2000) and satellite ocean color remote sensing studies (Arrigo et al., 1998a(Arrigo et al., , 1999Arrigo and Van Dijken, 2003;Park et al., 2019). These studies reported on two regions within the Ross Sea that exhibit physical and biogeochemical differences: (i) the weakly stratified and deeply mixed central Ross Sea polynya surface waters generally dominated by Phaeocystis antarctica, and (ii) the southwest (SW) marginal ice zone with intense surface stratification and shallow mixed layer depths dominated by diatoms. These recurrent and spatially variable distributions of phytoplankton biomass and taxa have important implications to food-web dynamics (Ainley et al., 2006;Davis et al., 2017), and the fate of dissolved (DOM) and particulate organic matter in the Ross Sea (Asper and Smith, 1999;Carlson et al., 2000;DiTullio et al., 2000;Carlson, 2002;Smith et al., 2011;Smith and Jones, 2015).
Phytoplankton production contributes to significant accumulation of newly produced dissolved (DOC) and particulate organic carbon (POC) in the surface ocean and its subsequent export to the deep ocean Carlson, 2002). DOC derived from primary production may originate from several biological processes including direct phytoplankton exudation, grazing interactions and viral lysis, while DOM removal processes include microbial mineralization and photodegradation (Carlson, 2002;Jiao et al., 2010;Ruiz-Halpern et al., 2011;Carlson and Hansell, 2015). In the Ross Sea, DOC and POC concentrations in late winter have been observed to be at background levels of ∼42 and 3 µM, respectively . An increase by as much as 30 and 107 µM of DOC and POC in excess of that observed in late winter has been linked to the high Ross Sea spring/summer production. However, the fraction of carbon fixed as DOC was found to be qualitatively more labile and enriched with nitrogen compared to the deeper refractory DOM (Carlson et al., 1998. In addition to the bulk analysis of DOC and dissolved organic nitrogen, other DOM constituents such as dissolved combined neutral sugars have been found to play an important role during the summer phytoplankton bloom in the Ross Sea . Although studies have hypothesized on the role of plankton community structure on DOC production and accumulation, the absence of large DOC accumulation within the Ross Sea has been attributed to low bacterial activity . More recently, Bercovici et al. (2017) showed that DOC produced in the Ross Sea enriches the dense shelf water (DSW) by ∼7 µM and contributes to the export of ∼4 Tg C yr −1 of DOC off the Ross Sea and into the abyssal Southern Ocean.
Chromophoric dissolved organic matter (CDOM), an optically active constituent of the DOM pool with absorption and fluorescent properties, could provide additional insights into DOM composition, source and dynamics in the Ross Sea. CDOM, a fraction of the DOM pool that absorbs UV and visible light, influences light penetration and primary productivity in aquatic ecosystems, plays a key role in photochemically induced transformations in surface waters, and acts as a 'sunscreen' protecting organisms from UV damage (Mopper and Kieber, 2002;Coble, 2007). Further, its absorption in the visible spectrum, especially in the blue band, tends to interfere with satellite ocean color estimates of chlorophyll a (D'Sa and Miller, 2003;D'Sa et al., 2006;Liu et al., 2021). CDOM is produced by a variety of autotrophic and heterotrophic organisms, influenced by physical processes such as circulation, upwelling and mixing and removed by photochemical degradation and microbial consumption (Blough and Del Vecchio, 2002;Nelson et al., 2004;Steinberg et al., 2004;Ortega-Retuerta et al., 2009;Romera-Castillo et al., 2010). CDOM absorption properties such as absorption coefficients at a specific wavelength and the slopes across specific regions of the optical absorption spectrum have been widely used to infer DOM sources and composition (Blough and Del Vecchio, 2002;D'Sa, 2008;Helms et al., 2008;D'Sa and DiMarco, 2009;Stedmon and Nelson, 2015;Hansen et al., 2016). CDOM abundance represented by absorption coefficient at 325 nm (a cdom 325) has been used to identify CDOM sources , as a tracer of biochemical processes in the global ocean (Nelson et al., 2010;Catala et al., 2015) and to study its distribution and reactivity in the Southern Ocean (Ortega-Retuerta et al., 2010a;D'Sa and Kim, 2017). Absorption spectral indices such as the spectral slope S over narrow wavelength intervals (e.g., 275-295 nm; 350-400 nm) or their ratio (S R ) provide information on CDOM photo-oxidative state, molecular size distribution, and microbial activity (Helms et al., 2008(Helms et al., , 2013. The specific UV absorbance (SUVA) or the DOC normalized absorbance at 254 nm (SUVA 254 ) is a useful parameter for estimating the dissolved aromatic content in aquatic systems (Weishaar et al., 2003). Although CDOM optical properties have been widely used to characterize DOM in aquatic systems, only a few studies have been reported for the Southern Ocean, including regions in the Antarctic Peninsula and Australasian sector (Clementson et al., 2001;Ortega-Retuerta et al., 2009, 2010aDel Castillo and Miller, 2010) with limited measurements in the Ross Sea (Kieber et al., 2009;D'Sa and Kim, 2017). Experimental studies in the Southern Ocean waters around the Antarctic Peninsula have documented the photoreactive nature (photobleaching and photohumification) of CDOM (Ortega-Retuerta et al., 2010a), and contribution by bacterioplankton and Antarctic krill to CDOM and DOC (Ortega-Retuerta et al., 2009;Ruiz-Halpern et al., 2011), suggesting the important role of these processes in carbon cycling. Insights on the role of seasonal sea ice dynamics and related growth of different biological communities in regulating CDOM and particulate matter or CDM in waters around the Antarctic Peninsula have also been obtained at larger spatiotemporal scales using satellite ocean color (Ortega-Retuerta et al., 2010b). However, only limited information on the absorption characteristics of CDOM, phytoplankton and non-algal particles and their relative contributions to the total light absorption budget that determines the optical variability of oceanic waters and thus ocean color have been reported for the Southern Ocean (Arrigo et al., 1998b;Reynolds et al., 2001;Ortega-Retuerta et al., 2010b;Hirawake et al., 2011;Lee et al., 2011;D'Sa and Kim, 2017). Knowledge of absorption properties of CDOM and other absorbing constituents could enhance understanding of DOM composition and improve ocean color algorithms for the Ross Sea.
Fluorescent DOM (FDOM), a fraction of CDOM that fluoresces, has also provided insights into DOM composition (Coble, 1996(Coble, , 2007Stedmon and Nelson, 2015). Excitation-emission matrix spectroscopy (EEMs), wherein threedimensional fluorescence measurements are performed across a range of excitation and emission wavelengths, are characterized by fluorescence in the UVA and visible spectral regions that have been associated with protein-like and humic-like fluorescent material (Coble, 2007;Stedmon and Nelson, 2015). Visible fluorescence is ubiquitous in all aquatic environments and spectral peaks observed in the EEMs spectra have been commonly referred to as A, C, and M (Murphy et al., 2018). Fluorescence observed with narrower peaks in the UVA spectral region has been associated with protein-like fluorescent material with characteristics similar to tryptophan-or T-and tyrosine-like amino acids likely derived from marine planktonic organisms or bacterial activity (Yamashita and Tanoue, 2003;Coble, 2007;Romera-Castillo et al., 2011). FDOM has also been investigated through molecular characterization of size fractions that suggest that DOM fluorescence characteristics may be controlled by molecular assemblies with similar optical properties and distributed along the molecular weight continuum (Romera- Castillo et al., 2014). The application of a multivariate tool such as parallel factor analysis (PARAFAC) to decompose a large EEM dataset into least-squares sum of several mathematically independent components allows for the efficient resolution and identification of different classes of fluorophores in an environment (Stedmon et al., 2003;Stedmon and Bro, 2008). Other fluorescent DOM compositional indicators obtained from the EEM spectra such as biological index (BIX) can also provide insights into the nature of the DOM pool (Huguet et al., 2009;Hansen et al., 2016). Although DOM fluorescence studies have been conducted in various regions of the Southern Ocean (Wedborg et al., 2007;Yamashita et al., 2007;D'Sa and Kim, 2017;Chen et al., 2019;Jeon et al., 2021), no such study has been reported for the Ross Sea. Thus, DOM fluorescence properties could further aid the characterization of the DOM pool in the Ross Sea.
In this study, we examine water column DOC concentrations, CDOM absorption (abundance, spectral slopes and ratios, UV specific absorption coefficient) and fluorescence (EEMs, PARAFAC components, BIX) properties in conjunction with biophysical (e.g., salinity, temperature, and chlorophyll) properties obtained during a field campaign in the Ross Sea conducted onboard the Korean ice breaker and research vessel Araon in the austral summer of 2014-2015. Spectral absorption properties of phytoplankton were also used to examine phytoplankton characteristics and to assess CDOM contribution to the total absorption budget, an important factor that influences ocean color algorithms. In addition, satellitederived sea ice and chlorophyll estimates were used to assess synoptic large scale features to aid the interpretation of DOM optical data in the Ross Sea.

Study Area, Water Masses and Sampling
The Ross Sea lying to the north of the RIS in the Antarctic continental shelf has an average depth of ∼500 m with an irregular bottom topography characterized by several banks and troughs and a shelf break running NW-SE that corresponds to the 700 m isobath ( Figure 1A). The hydrographic structure of the Ross Sea has been described in great detail in earlier studies (Budillon et al., 2003;Orsi and Wiederwohl, 2009). Briefly, meridional overturning brings the Circumpolar Deep Water (CDW: a relatively warm, salty and oxygen-poor layer within the Antarctic Circumpolar Current) poleward between the much colder Antarctic Surface Water (AASW) above and the AABW below (Orsi and Wiederwohl, 2009). The intrusion of CDW onto the Ross continental shelf (a source of heat and nutrients) and its interactions with the shelf waters results in the modified CDW (MCDW) which is characterized by a subsurface maximum temperature and minimum dissolved oxygen (Jacobs et al., 1985;Budillon et al., 2003). Sea-ice formation in the shelf transforms upper waters (either the near-freezing AASW or the shoaling MCDW) into cold, saltier and denser shelf water (DSW) which fills the bottom layer of the Ross Sea including under the RIS (Orsi and Wiederwohl, 2009).
Stations were occupied onboard the ice breaker RV Araon in the Ross Sea ( Figure 1A) from January 11 to February 04, 2015 (austral summer). Transect stations (sts) were sampled east of Victoria Land from the northwest (NW: sts 9, 10, and 11) to southwest (SW: sts 12, 15, and 25), north of the Ross Ice Shelf (hereafter denoted as off-RIS: sts 28 and 29) and a transect T1 (sts 30-38) inshore of the shelf break. Two additional stations (st 63 in the NW and st 71 north of the Ross Island and RIS) were also sampled toward the end of the cruise. Physical data plotted from the CTD casts obtained during the approximately month-long field campaign show the presence of the main water Potential temperature-salinity diagram with isopycnals (σ θ ) obtained from CTD profiles at the sampling stations. Water masses CDW, MCDW, AASW and DSW as defined by Jacobs et al. (1985) and Orsi and Wiederwohl (2009). masses (CDW, MCDW, DSW and AASW) at our sampling stations ( Figure 1B).

Hydrography Data
Conductivity-temperature-depth (CTD) measurements in the water column were recorded using a CTD-Rosette system (Sea-Bird SBE-911plus). An oxygen sensor (SBE 43) and a chlorophyll fluorescence sensor (WET Labs ECO-AFL/FL) were connected to the CTD system for dissolved oxygen and chlorophyll fluorescence measurements. The WET labs fluorescence sensor was factory calibrated for chlorophyll a. All data were averaged into 1-m depth intervals to reduce noise. For measurements of DOC concentrations and optical measurements (CDOM absorption and fluorescence; particulate absorption) samples were collected at multiple depths in the water column (usually at surface, 10/20, 50, 100, and 150 m; and few stations up to 200 m) using CTD-mounted Niskin water samplers. All samples were processed in the laboratory within ∼4 months of sample collection. However, due to the long duration of the field campaign, samples for CDOM absorption and fluorescence were processed in the laboratory between 3.0 to 3.6 months after sampling. While this strategy ensured any bio-alteration of the optical signals to be consistent across the samples, it could likely cause preferential removal of the bio-labile components.

Dissolved Organic Carbon Measurements
Samples for DOC were filtered through pre-rinsed Whatman GF/F filters and stored in acid cleaned, pre-combusted amber bottles with Teflon-lined caps. Laboratory measurements of DOC were made on a Shimadzu TOC 5000A (with ASI-5000A autosampler) using a high temperature combustion method to convert carbon compounds to carbon dioxide, including using the Consensus Reference Material (CRM) for QA/QC (Benner and Strom, 1993;Hansell and Carlson, 2001).

Chromophoric Dissolved Organic Matter Absorption Measurements and Analysis
For CDOM spectral absorption measurements, water samples were filtered through 0.2 µm pre-rinsed Nuclepore filters on the same day and stored in acid cleaned, pre-combusted amber bottles with Teflon-lined caps at 4 • C in the dark and processed for CDOM spectral absorption using a 10 cm quartz cuvette on a dual beam Perkin-Elmer Lambda 850 spectrophotometer equipped with a 150 mm integrating sphere. We note that while effects of storage times (∼21 days) on CDOM absorbance have been reported in wastewater effluents (Sgroi et al., 2020), a longer hold time of 6 months of filtered seawater kept at 4 • C indicated no change in absorption (Mannino et al., 2019). Details on the measurements and calculations of the absorption coefficient (a cdom , m −1 ) and spectral slopes S 275−295 and S 350−400 (µm −1 ) are given in D' Sa et al. (2014) and Joshi et al. (2017). Absorption coefficients at various wavelengths (e.g., 254, 300, 325, 355, and 443 nm) were obtained and used for comparisons to other study regions. Spectral slope ratio (S R ) was calculated as S 275−295 divided by S 350−400 (Helms et al., 2008). SUVA 254 (m 2 mg −1 C) was calculated by dividing the absorbance at 254 nm in inverse meters (m −1 ) by the DOC concentrations.

Phytoplankton and Non-algal Matter Absorption Measurements
For determination of light absorption by phytoplankton and nonalgal matter, generally 500 ml volume of water were collected from the Niskin bottles and filtered onto 25 mm Whatman GF/F filters with nominal pore size of 0.7 µm. The filters were immediately stored in a freezer (−80 • C) during the cruise and later transferred to liquid nitrogen before laboratory spectroscopic measurements. Particulate (total -a p and nonalgal -a nap ) absorption were measured inside the integrating sphere on a Perkin-Elmer Lambda 850 spectrophotometer (details in Naik and D'Sa, 2012;Roesler et al., 2018;Liu et al., 2021). Phytoplankton spectral absorption coefficient a phy were then determined as the difference between a p and a nap . Chlorophyll specific phytoplankton absorption at 443 nm, a * phy 443 (m 2 (mg chl) −1 ) was obtained by dividing phytoplankton absorption at 443 nm, a phy 443 by chlorophyll.

Excitation-Emission Matrix Spectroscopy Fluorescence Measurements and PARAFAC Analysis
Filtered seawater samples used for CDOM absorption measurements were also used to record EEMs using a FluoroMax-4 (Jobin Yvon Horiba) fluorometer. EEMs were recorded on the fluorometer by scanning the emission spectra from 290-550 nm at 5 nm intervals over excitation wavelengths between 250-450 nm at 5 nm increments. The EEMs spectra were obtained after correction of the fluorescence spectra for instrument bias, and the water Raman normalization of the fluorescence intensity (Singh et al., 2010;D'Sa et al., 2014). Due to low absorbance of the samples (0.004 to 0.015 cm −1 at 254 nm), there was no requirement for inner filter correction of the fluorescence data. The final EEMs fluorescence values are reported in equivalent Raman units (RU). The resulting EEM fluorescence observations were evaluated by PARAFAC analysis using the DOM-Fluor toolbox (Stedmon and Bro, 2008) with the model constrained by non-negativity, and run with three to seven components. Model validation was carried out using split-half analysis, residual analysis and random starts (Stedmon et al., 2003;Stedmon and Bro, 2008;D'Sa et al., 2016). A few samples were not used (e.g., Supplementary Figure 2) due to the high leverage (sample deviation from average distribution; Stedmon and Bro, 2008) and large residual EEMs (difference between sample and model derived EEMs) that improved the results of split-half analysis in the final model construction. Residual EEMs also indicated the presence of small signals in the UV corresponding to the tyrosine fluorescence but were not resolved by PARAFAC analysis due to the relatively strong tryptophan-like fluorescent signature; we note that the effects of sample hold time on the tyrosine-like FDOM cannot be ruled out due to its very labile nature (Chen and Jaffé, 2014). A fluorescence index BIX (index of recent autochthonous contribution) was calculated at 310 nm excitation, as the ratio of emission intensities at 380 nm divided by 430 nm (Huguet et al., 2009).

Satellite Data
Satellite sea ice concentration maps of the Ross Sea for the months of December 2014 and January 2015 were generated from data downloaded from the NOAA National Snow and Ice Data Center Climate Data Record of Passive Microwave Sea Ice Concentration, Version 4 1 . The monthly average of chlorophyll imagery was downloaded from NASA 2 . MODIS-Aqua level-2 products for the Ross Sea region (e.g., chlorophyll a concentration estimated by the standard OCx empirical algorithm and the phytoplankton absorption coefficient a phy 443 from the NASA GIOP semi-analytical algorithm) were downloaded for the month of January 2015 from the NASA website for comparisons to field data.

Ross Sea Polynya and Surface Chlorophyll Conditions in Summer
The Ross Sea polynya which is maintained during winter by strong katabatic winds blowing offshore off the RIS, starts expanding during spring and into summer with increasing solar insolation and warming (Park et al., 2018). For the summer of 2014-2015, we observe the polynya extent for December and its northwards expansion in January (Figures 2A,B), along with a marginal ice zone to the west and east of the polynya. A large phytoplankton bloom was observed in December north of the RIS and the SW, and a smaller bloom in the Terra Nova Bay polynya located off Victoria Land ( Figure 2C). By January, surface chlorophyll decreased considerably, but remained elevated in the SW and TNB, with lower concentrations in the central and northern part of the shelf ( Figure 2D).

Water Column Characteristics
Water column properties (temperature, salinity and chlorophyll fluorescence) and mean values in the upper water column (∼150 m) at discrete sampling depths are shown for the different regions of the Ross Sea (Figure 3 and Table 1). Transect T1 (Figures 1A, 3A,B) reveals typical water masses and its characteristics just inshore of the shelf break. The AASW with its low salinity is observed within the upper 100 m. The transect plot includes st 61 ( Figure 1A; no water sampling) at the northwestern end of the transect located just inshore of the shelf break, to demonstrate the subsurface intrusion of the warm, high salinity and low oxygen (Supplementary Figure 1A) CDW waters. However, along the main T1 transect, inflow of the MCDW (warm and low oxygen waters) into the shelf is observed at several locations generally below the low salinity AASW and above the cold and high salinity DSW, with DSW dominating over the troughs and the MCDW over the banks. Greater mixing of the MCDW with surface waters and increased depth of the surface mixed layer eastward appear to contribute to elevated subsurface phytoplankton biomass (Figures 3A-C).
Vertical profiles of the same properties at stations in the SW (sts 12, 15, and 25), off-RIS (sts 28 and 29) and station 71 reveal differences in biophysical properties (Figures 3D,E). A shallow mixed-layer depth is observed at the SW stations with strong thermocline/halocline within the upper 30 m, with subthermocline temperatures (< −1.85 • C) and salinity (> 34.65) at > 50 m associated with the cold, high salinity shelf waters (DSW). Shallow subsurface peaks in chlorophyll of 11.6, 3.8, and 8.0 mg m −3 were recorded at stations 12, 15, and 25, respectively. However, at the nearshore station north of the Ross Island in the SW (st 71, Figure 1A) water column temperature, salinity, chlorophyll and dissolved oxygen differed from stations in the SW and off-RIS (Figures 3D,E; Supplementary Figure 1B). At this station, a broad chlorophyll peak of ∼6 mg m −3 present at ∼30 m depth decreased to levels of ∼1.8 mg m −3 at 100 m depth, and thereafter remained elevated albeit decreasing with depth. A similar elevated chlorophyll profile with depth was also observed in the NW nearshore station 63 (Figure 3E). At the two RIS stations (sts 28 and 29) with deeper mixed layers, surface temperatures were relatively warmer (0.65 and −0.35 • C) and waters fresher (34.27 and 34.18). Chlorophyll increased from a low of 0.76 and 4.3 mg m −3 at the surface, to peak values of 15.6 mg m −3 at a depth of 68 m for station 28 and a broad subsurface bloom of ∼16 mg m −3 extending down to a depth of ∼120 m at station 29 ( Figure 3E).

Dissolved Organic Carbon Concentrations and Chromophoric Dissolved Organic Matter Absorption Properties
Distribution of DOC in the upper water column (depth range: 0 to ∼150 m) at each of the sampling stations corresponding to the three regions of the shelf are shown in Figure 4A. DOC concentrations over the whole shelf (mean: 47.82 ± 5.70 µM) ranged between 38.75-67.58 µM (Figure 4A), with means of 47.69 ± 7.99, 48.73 ±.0.5.71, and 49.46 ± 8.04 µM, respectively, in the NW + T1, SW and off-RIS ( Figure 1A and Table 1). Corresponding upper water column distribution of CDOM absorption coefficient a cdom 325 in the shelf ranged from 0.11 to 0.81 m −1 (Figure 4B) with a mean of 0.31 ± 0.18 m −1 , while spectral slope S 275−295 ranged from 21.20 to 33.72 µm −1 (Figure 4C), with a mean of 28.12 ± 3.24 µm −1 , respectively. The highest a cdom 325 and lowest S 275−295 was in the SW region (Figures 4B,C), with mean values of 0.63 ± 0.19 m −1 and 24.06 ± 2.93 µm −1 , respectively ( Table 1). With the exception of the high values in the SW, a cdom 325 for the other regions (i.e., NW + T1, off-RIS, sts 71 + 63) ranged between 0.11 -0.47 m −1 , with a mean of 0.25 ± 0.08 m −1 . Mean CDOM absorption coefficients at 300, 355, and 443 nm are also shown (   for all the sampling locations in the shelf show a mean value of 1.22 ± 0.44 m 2 g −1 C, with highest value in the SW region (1.85 ± 0.55 m 2 g −1 C). Example a cdom spectra at station 33 (NW + T1) and stations 15 and 25 (SW) are shown in Figure 5.
Relationships between a cdom 325 and salinity, chlorophyll and DOC were also examined. a cdom 325 showed a slight positive trend with salinity, but was uncorrelated to chlorophyll or to DOC (Figures 6A-C). The relationship between a cdom 325 and spectral slope S 275−295 showed a general decreasing trend with increasing a cdom 325 ( Figure 6D).

Fluorescence Dissolved Organic Matter Properties
An example of EEMs spectra obtained at three depths (surface, 100 m and 200 m) at station 71 is shown in Figures 7A-C.
In the visible, the humic-like peaks were similar to the A and C peaks (Coble, 1996) at excitation/emission (ex/em) of ∼ < 260/480 nm and 370/480 nm, while a blue-shifted peak was observed at ∼305/400 nm corresponding to the M-peak. In the UVA, peaks were observed at ∼275/335 nm that have been associated with tryptophan (or T-peak; Coble, 1996); although the T-peak was the main signal observed in most samples (Figure 7), comparatively weaker emission signals were also observed at shorter wavelengths that suggest the presence of other fluorophores such as the tyrosine protein-like fluorophore (e.g., Figure 7C). However, since tyrosine and tryptophan do not typically occur as pure dissolved amino acids in the environment (Yamashita and Tanoue, 2003), the broad proteinlike fluorescence with peak ex/em similar to T-peak can be considered as likely derived from a mixture of dissolved amino acids and other organic materials with similar fluorescence characteristics (Murphy et al., 2018). The PARAFAC analysis of DOM excitation-emission scans identified three major fluorescent fractions (C1, C2 and C3) and are shown along with their excitation and emission loadings (Figure 8). C1 with ex/em of 275/335 nm is associated with tryptophan-like fluorophore (Yamashita and Tanoue, 2003;Romera-Castillo et al., 2011). The C2 component with ex/em peaks at < 260(370)/480 nm corresponds to previously noted A/C humic-like material generally ubiquitous in the marine environment (Murphy et al., 2018). The C3 component with ex/em peak at 305/400 nm is similar to the "M" or marinehumic like material of biological and/or microbial origin (Coble, 1996;Stedmon et al., 2003;Romera-Castillo et al., 2011). The humic-like C2 was on average greater than the C3 component, while the protein-like C3 component showed greatest variability ( Table 2). The DOC normalized C2 and C3 components showed weak but significant positive correlations to salinity (Figure 9A), while the sum of the two humic-like components (C2 + C3) also showed weak but significant correlation to DOC concentrations ( Figure 9B). The average biological index BIX was 1.0 ± 0.1 for the whole study region with little variability between the different sampled regions ( Table 2).

Phytoplankton Absorption and Absorption Budget
The spectral characteristics of phytoplankton absorption coefficient a phy show differences in both magnitude and spectral  shape in the Ross Sea (Figure 10). In the UV region, a peak at 320 nm attributed to photoprotective pigments was observed for most surface samples with the largest signal observed for the waters in the SW (Figure 10A). These peaks were mostly low or absent in subsurface samples, but remained elevated at stations 28 and 29 off the RIS (Figure 10B). In the visible (400-700 nm), absorption spectra exhibited a clear absorption peak at 443 nm associated with chlorophyll a pigment. Absorption spectra of surface samples ( Figure 10A) were highest for the SW (sts 12, 15, 25 and 71) and lowest along the transect T1, and intermediate for the NW (sts 9, 10, and 11) and off the RIS (sts 28 and 29). Absorption spectra in the upper water column (4 depth intervals, generally up to 100 m or greater at a few at stations) show clear spectral differences (Figure 10B) that can be attributed to differences in phytoplankton taxa. Stations in the SW exhibit spectra with a single peak at 443 nm that show similarity to the diatom Fragilariopsis cylindrus (Arrigo et al., 2010), while at the two off-RIS stations (sts 28 and 29), a strong secondary absorption peak between 450-470 nm is observed for all the sampled depths that show similarity to the haptophyte Phaeocystis antarctica (Arrigo et al., 2010). At station 71, a chlorophyll spectral signature is clearly observed even at deeper depths (100 and 200 m). It is important to note, this species identification is based on visual comparisons of phytoplankton absorption spectra and is not confirmed by any explicit biological characterization.
Chromophoric dissolved organic matter absorption along with particulate matter absorption which includes phytoplankton a phy and non-algal particles a nap are examined at 443 nm as it influences the light field and ocean color. a nap 443 showed good correlation to a phy 443 (r 2 = 0.72) (Figure 11A), while a cdom 443 was uncorrelated to both a phy 443 and a nap 443. We examined the relative contribution of a phy 443, a nap 443 and a cdom 443 to total non-water absorption at 443 nm (a tnw 443 = a phy 443 + a nap 443 + a cdom 443) on a ternary plot ( Figure 11B). The ternary plot shows the relatively strong contribution by CDOM in the visible with a cdom 443 contributing between 40-90% (mean 70.7 ± 18.3%), a phy 443 about 22.5 ± 15.2% and a nap 443 only 6.1 ± 15.2% to the absorption budget; a caveat in this analysis is the unaccounted absorbance for particles between 0.2 µm and 0.7 µm due to filters used for CDOM and particulates, respectively.
We also examined the specific phytoplankton absorption spectra which are important in bio-optical models (Bricaud et al., 1998;Morel and Maritorena, 2001). A plot between chlorophyll and specific phytoplankton absorption at 443 nm (a * phy 443) for surface waters (0 to ∼20 m) shows decreasing a * phy 443 with increasing chlorophyll (Figure 11C), with SW stations showing greater values of a * phy 443 for similar range in chlorophyll concentrations. For comparison, we show the Bricaud et al. (1995) relationship which was derived from an extensive data set collected from six oceanic regions, but not including high latitude regions. Comparisons between in situ and average MODIS-derived chlorophyll using OCx algorithm and a phy 443 using the GIOP algorithm are shown for January 2015 (Figure 11D).

The Ross Sea Biophysical Properties
The Ross Sea polynya which increases with solar insolation during spring-summer, expanded rapidly northward from Dec into Jan 2014-15 austral summer (Figures 2A,B), with accompanying changes in surface algal biomass (Figures 2C,D). Such evolution of the Ross Sea polynya, which is also influenced by the variability in wind intensity and climate feedbacks (Park et al., 2018), plays an important role in the control of phytoplankton blooms and the biogeochemical cycling in the Ross Sea (Arrigo et al., 1998a;Dunbar et al., 1998;  Relationships between (A) salinity and DOC normalized C2 (C2*) and C3 (C3*), (B) DOC and C2 + C3 fluorescent components for samples obtained in the upper water column (e.g., Figure 4). Asper and Smith, 1999;DiTullio et al., 2000). Our study area which encompassed most of the significant water masses in the Ross Sea (Figures 1A,B) included: (i) the northwest and a transect in central Ross Sea inshore of the shelf break (NW + T1), (ii) the southwest (SW), and (iii) off the RIS in the Ross Sea polynya (off-RIS). These three areas have been previously characterized into distinct biogeochemical regimes based on large gradients in biomass and phytoplankton community composition and differences in water column structure and stability, namely (i) the low-productivity NW and central shelf break region with diatom-or mixed diatom and Phaeocystisdominated region, (ii) the diatom-dominated ice edge bloom region in the SW, and (iii) the Phaeocystis-dominated region north of the RIS (Goffart et al., 2000;Sweeney et al., 2000). Along the northwest and inshore of the inner shelf transect (NW + T1), we observe the subsurface presence of the MCDW (warmer, saltier, and lower dissolved oxygen) at various locations linked to CDW incursions and mixing, and the inflow of nutrients and heat into the shelf (Jacobs et al., 1985). Average chlorophyll concentration in the NW + T1 section was the lowest compared to the SW and off-RIS (Figure 3 and Table 1). In the SW region with its shallow mixed layer, the phytoplankton absorption spectra (Figure 10) show similarities to Fragilariopsis cylindrus (Arrigo et al., 2010) and is consistent with past studies that indicate this region to be dominated by diatoms, especially the Fragilariopsis spp. Goffart et al., 2000).
At station 71 north of the Ross Island, elevated chlorophyll at subsurface and deeper depths (Figure 3) suggest sinking of the phytoplankton biomass, with phytoplankton absorption spectra showing similarities to diatom Fragilariopsis cylindrus (Figure 10; Arrigo et al., 2010). Rapid organic matter flux at settling rates of 60 to > 400 m d −1 have been reported in the SW region north of the Ross Island (Dunbar et al., 1998;Asper and Smith, 2003). Sustained elevated levels of phytoplankton biomass in the SW region especially north of Ross Island observed in satellite imagery during Dec and Jan (Figures 2C,D) could be a factor in the export of organic material to depths. North of the RIS, we observe deeper mixed layer depths with high subsurface chlorophyll concentrations (Figures 3D,E). Phytoplankton absorption spectra at these stations (Figure 10) show similarity to spectra of P. Antarctica with two peaks in the blue wavelengths (Arrigo et al., 2010), including a marked shoulder at ∼475 nm usually associated with alloxanthin and 19 -hexanoyloxyfucoxanthin pigments present in haptophyte phaeocystis; this is consistent with the dominance of P. antarctica generally reported in the region north of the RIS in the central Ross Sea polynya (Arrigo et al., 1999Goffart et al., 2000;Sweeney et al., 2000).

Absorption Properties and Implications for Ocean Color
The spectral absorption properties of CDOM a cdom λ, phytoplankton a phy λ and non-algal particles a nap λ are important parameters to understand the behavior of underwater light field and in bio-optical models. The surface distributions of a cdom , a phy and a nap at 443 nm revealed distinct trends with higher absorption in the SW region of the Ross Sea ( Figure 11A). Particulate matter absorption in the Ross Sea was dominated by phytoplankton, with a phy 443 and a nap 443 well correlated. a cdom 443, which was the dominant contributor to the absorption budget (Table 2 and Figure 11B), was decoupled from a phy 443, revealing that CDOM varied independently of phytoplankton in spite of its dominant role as a source of excess CDOM in the Ross Sea; this knowledge can be useful in improving satellite ocean color algorithm performance in the Ross Sea (e.g., Figure 11D).
The specific phytoplankton absorption a * phy λ showed variability in the blue (443 nm) with a * phy 443 decreasing from 0.120 to 0.005 m 2 (mg chl a) −1 with increasing chlorophyll (0.05 to 12 mg m −3 ). However, distinct trends in the SW suggest differences in pigment packaging and/or change in pigment composition (Bricaud et al., 1995) that is consistent with the phytoplankton taxa reported in past work Arrigo et al., 1999Arrigo et al., , 2000Smith et al., 2014), and as observed in the phytoplankton absorption spectra (Figure 10); further, higher a * phy 443 also suggest better light adaptation by diatoms in the SW. The a * phy 443 versus chlorophyll relationship shows similarity to other high latitude regions (e.g., Naik et al., 2013), but is lower than the fit reported by Bricaud et al. (1995), which was obtained for middle and lower latitude waters ( Figure 11C). These differences in the a * phy λ in the blue and other wavelengths including the green band has been shown to lead to underestimations in satellite derived chlorophyll (Arrigo et al., 1998b;Dierssen and Smith, 2000), as also demonstrated in our study ( Figure 11D). We also observe that the satellite imagery misses the large sub-surface blooms observed north of the RIS in the Ross Sea polynya, thus likely underestimating satellite-derived productivity in the Ross Sea. This study demonstrates the need to better characterize linkages between bio-optical and biophysical systems in this highly productive region of the Southern Ocean.

Dissolved Organic Carbon and Dissolved Organic Matter Optical Properties and Sources
In spite of the various water masses, algal blooms and dominance of different phytoplankton species within different regions of the Ross Sea, DOC concentrations (mean: 47.82 ± 5.70 µM) during late summer varied over a relatively narrow range (38.75-67.58 µM; Figure 4A) consistent with previous studies (Carlson et al., 1998 and was ∼6.02 µM greater than the mean concentrations (41.8 ± 0.80 µM) obtained during late winter/early spring when waters are well mixed over the entire water column of the Ross Sea . This excess summer DOC accumulation is slightly lower than the ∼7 µM enriched DOC in the DSW relative to the incoming MCDW reported by Bercovici et al. (2017). The relatively small magnitude of DOC production and accumulation in the Ross Sea was found to be nitrogen-rich compared to background DOM and referred to as "semi-labile" DOM , with overall mean DOC concentrations at similar levels in the different regions of the Ross Sea ( Table 1).
As in the case of DOC, the distribution of the three FDOM components showed no discernable regional trends in the Ross Sea (Figure 9 and Table 2). The three identified PARAFAC components (Figure 8 and Table 2) include C1 protein-like (predominantly Coble, 1996 peak T), C2 humic-like (peaks A and C) and C3 humic-like (peak M); these components are similar to other widely derived FDOM fractions in the global ocean (Yamashita and Tanoue, 2003;Coble, 2007;Jorgensen et al., 2011;D'Sa et al., 2014;Murphy et al., 2018) including the Southern Ocean (Wedborg et al., 2007;D'Sa and Kim, 2017;Chen et al., 2019;Jeon et al., 2021). However, while the three retrieved PARAFAC components in the Ross Sea were similar to that reported in the Amundsen Sea, Antarctica (Jeon et al., 2021), these FDOM components differed from another study in the same region by the absence of the tyrosine-like component linked to high abundance of low molecular weight DOC (Chen et al., 2019). An approach that combines both FDOM and DOM molecular size fractions (e.g., Chen et al., 2019) may help elucidate similarities or differences in the DOM pool within the same or different regions of the Southern Ocean.
Fluorescence dissolved organic matter production of Tand M-like fractions have been shown to be exuded by marine phytoplankton (Romera-Castillo et al., 2010), while the formation of more complex molecules associated with peaks A and C could be mediated by the bacterial degradation of phytoplankton-derived organic matter (Kinsey et al., 2018;Goto et al., 2020), suggesting that in remote locations such as the Ross Sea, autochthonous FDOM production associated with the spring/summer phytoplankton blooms would be expected. As an example, Figure 7 demonstrates increasing levels of both protein-and humic-like FDOM with increasing depths (surface to 200 m) at a station with sinking algal biomass (Figures 3, 10) undergoing remineralization (as indicated by the relatively low water column dissolved oxygen (Supplementary Figure 1B). Similar increases of both protein-and humic-like FDOM were observed at the subsurface P. antarctica bloom stations off-RIS ( Supplementary Figure 2A) that show similarity to EEMs spectra described in Goto et al. (2020). The relatively high mean BIX value (> 1.0), the biological index parameter for which values > 0.8 suggest DOM with an obvious component of freshly produced organic matter of biological origin (Huguet et al., 2009), also support FDOM production in the Ross Sea. The generally increasing trend of DOC normalized C2 and C3 component intensities ( Figure 9A) with increasing salinity suggest both effects of FDOM photooxidation and/or production in the upper water column of the Ross Sea. The weak, but significant correlation between DOC concentrations and the C2 + C3 humic-like components (Figure 9B) also indicate a close linkage between DOM and FDOM in the Ross Sea. The C1 protein-like component was however, more variable ( Table 2) and varied over a greater range extending from below the detection limit to a high of 0.7 RU; such variability could be due to both production and consumption of colored materials by phytoplankton and bacteria (Romera-Castillo et al., 2011).
In contrast, CDOM showed significant differences between SW (a cdom 325: 0.63 ± 0.19 m −1 and spectral slope S 275−295 : 24.06 ± 2.93 µm −1 ) throughout the upper water column (Figures 4B,C) compared to the rest of the study area (a cdom 325: 0.25 ± 0.08 m −1 ; slope S 275−295 : 28.55 ± 1.90 µm −1 ). This compares to an increase in S 275−295 and a cdom 325 of 0.142 ± 0.01 m −1 reported before a Pheocystis antarctica bloom in Nov 2005 within the Ross Sea polynya (Kieber et al., 2009), suggesting an approximate doubling of CDOM accumulation (∼0.11 m −1 ) over the summer (we note that measurements were obtained 10-years apart on different instruments). These a cdom 325 values are higher than those reported for the AABW (∼0.19 ± 0.01 m −1 ; Catala et al., 2015) or the surface waters of the New Zealand sector of the Southern Ocean (∼0.10 ± 0.06 m −1 ; D'Sa and Kim, 2017). A decreasing trend of S 275−295 with increasing a cdom 325 ( Figure 6D) suggest greater photobleaching in the stratified waters of the NW and the central Ross Sea polynya. CDOM values however, were on average highest in the SW with lowest S 275−295 , showing generally similar properties within the upper 100 m including the shallow surface mixed layer (Figures 4, 6) suggesting greater molecular weight and aromaticity of DOM in the SW; SUVA 254 , a useful proxy for DOM aromatic content (Weishaar et al., 2003) also indicated DOM in the SW to be more aromatic suggesting an additional source of CDOM in the SW. The slope ratio S R which has been linked to shifts in DOM molecular weight and photobleaching (higher slope ratios reflective of greater amounts of low molecular weight compounds; Helms et al., 2008), were within the range of 1.60-3.17, with mean values of 2.31 ± 0.37 for the NW + T1, 2.48 ± 0.33 off-RIS, and 2.62 ± 0.58 for station 71, but lowest value of 2.23 ± 0.33 for the SW ( Table 1), suggesting relatively higher molecular weight DOM in the SW compared to other regions of the Ross Sea. Similar values have been observed in the subsurface subarctic waters and the deep ocean (Yamashita et al., 2013;Catala et al., 2015). Interestingly, an increase in the mean CDOM optical properties (a cdom 325, S 275−295 and S R : 0.21 ± 0.01 m −1 , 19.55 ± 1.0 µm −1 and 1.8 ± 0.1; Catala et al., 2015) in the CDW (the source waters for the Ross Sea) to (0.31 ± 0.18 m −1 , 28.12 ± 3.24 µm −1 , and 2.33 ± 0.38; Table 1) summer Ross Sea values strongly suggest DOM/CDOM enrichment in the Ross Sea with lower molecular weight fraction associated with the spring/summer phytoplankton bloom events. Further, lower values of CDOM optical properties in the AABW (linked to the outflowing DSW) which are observed to be similar to the CDW (Catala et al., 2015), suggest the eventual microbial degradation (Yamashita et al., 2013) of the Ross Sea DOM/CDOM. A similar process of microbial remineralization of the enriched Ross Sea DOM that is exported off the shelf has been proposed (Bercovici et al., 2017).

Chromophoric Dissolved Organic Matter Sources in the Southwest Ross Sea
While the Ross Sea DOC and CDOM enrichment (with reference to background values in late winter/early spring) can be attributed to the spring/summer phytoplankton blooms, the relatively high CDOM absorption and low spectral slope values throughout the upper water column in the SW region (Figures 4B,C and Table 1) appear to be linked to the dominance of diatoms in the highly stratified waters of the SW. We examine two factors that could have contributed to these distinct CDOM optical properties: (i) Phytoplankton absorption spectra of surface and some subsurface samples showed elevated values in the UV (300-400 nm) with peaks at ∼320 nm (Figure 10) mainly in the diatom-and haptophyte-dominated waters of the SW and off-RIS; these peaks have been attributed to mycosporine like amino acids (MAAs), a group of UV absorbing compounds that act as a sunscreen to reduce UV induced damage (Riegger and Robinson, 1997;Moisan and Mitchell, 2001). The amplitude of these absorption peaks was greater in the shallow mixed layer surface waters of the SW (Figure 10A) where, microphytoplankton such as diatoms have been shown to biosynthesize MAAs more efficiently in comparison to smaller nano-and picophytoplankton (Ha et al., 2018). Although relatively high levels of MAAs have been shown to influence CDOM absorption, especially in the UV region (D'Sa, 2008), elevated CDOM absorption in the surface as well as subsurface waters of the SW (Figure 4B) where phytoplankton absorption in the UV was low or negligible ( Figure 10B), suggest MAAs were an unlikely contributor to the CDOM anomaly in the SW.
(ii) Antarctic krill (Euphausia superba) and ice krill (E. crystallorophias) are the two dominant species found in the Ross Sea (Ainley et al., 2006;Smith et al., 2007) with Antarctic krill dominant in the northern and northwestern areas of the Ross Sea and ice krill predominant in the southwestern areas (Davis et al., 2017). Our sampling region in the SW which is close to Terra Nova Bay (e.g., station 15, Figure 1) has been shown to have high krill population (Sala et al., 2002), with krill aggregations observed to be distributed in discrete layers in the upper water column (Kang et al., 2020). Other dominant grazers such as copepods and Limacina helecina that perform diel migration in the upper 200 m have also been identified in the diatom dominated western Ross Sea (Goffart et al., 2000). Antarctic krill has been confirmed as a direct source of CDOM in the peninsular region of the Southern Ocean (Ortega-Retuerta et al., 2009), while CDOM production by zooplankton along with active DOC transport through zooplankton vertical migration has also been reported (Steinberg et al., 2002. In addition, we observe similarities of the SW CDOM absorption spectra (Figure 5, station 15) to that produced by krill (Ortega-Retuerta et al., 2009) and zooplankton . These combined factors strongly support grazer contribution to the CDOM anomaly observed in the SW.
The anomaly in the SW CDOM optical properties (low S 275−295 and high SUVA 254 ) also suggest grazer generated DOM of greater aromaticity and molecular weight (Weishaar et al., 2003;Helms et al., 2008) that demonstrate the important role that grazers such as Antarctic krill could play in the biogeochemical cycling in some regions of the Southern Ocean due to their extraordinary high biomass (Alcaraz et al., 2014;Cavan et al., 2019). Owing to potential for long-term changes in krill stock within the Southern Ocean (Atkinson et al., 2004), future field campaigns should integrate studies of krill and DOC/CDOM measurements in this region.

CONCLUSION
The Ross Sea, one of the most remote and productive shelf regions in the Southern Ocean, undergoes seasonal spring/summer phytoplankton blooms that play an important role in the marine carbon cycle. In this study, we examined upper water column DOC concentrations, CDOM absorption and fluorescent properties in conjunction with biophysical properties that indicated DOC and CDOM accumulation and FDOM production during austral summer in the Ross Sea. Consistent with past studies that identified major biogeochemical regimes in the Ross Sea, we observed differences in phytoplankton absorption properties (spectral absorption and specific absorption coefficients) that were characteristic of diatom-and phaeocystis-dominated regions in the SW and off the RIS. These properties, including the dominant contribution by CDOM to the total absorption, have implications for ocean color and demonstrate the need to better characterize bio-optical properties in the Ross Sea.
We observed no discernible regional trends in DOC concentrations in the Ross Sea, with the DOC accumulation during summer associated with the spring/summer phytoplankton blooms consistent with previously reported values. Similarly, FDOM fractions showed no discernible regional trends, with the DOC normalized humic-like C2 and C3 components and sum of humic-like FDOM components showing weak, but significant positive correlation to salinity and DOC, respectively. Further, high BIX value and an instance of increasing FDOM with depth at a location with sinking organic matter strongly suggest FDOM production in the Ross Sea.
In contrast, we observed regional differences in the CDOM absorption properties in the Ross Sea. High a cdom 325, low S 275−295 and high SUVA 254 in the SW suggest CDOM of greater aromaticity and molecular weight that appear to be linked to biogeneration by krill, the dominant planktonic grazer in the region. Overall, mean CDOM optical properties (a cdom 325, S 275−295 and S R ) in the Ross Sea revealed DOM of lower molecular weight during summer in comparison to the inflowing CDW source waters and the AABW linked to the outflowing DSW. This study enhanced our understanding of DOC, CDOM and FDOM characteristics and dynamics in the Ross Sea; it also provided insights on the potentially important contribution of Antarctic krill to carbon cycling that should be considered in future studies due to their extraordinary high biomass in the Southern Ocean.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
EJD and H-CK designed the study. EJD, H-CK, S-YH, and IJ contributed to this work through sample and data collection or processing of the field data. EJD analyzed the data and drafted the manuscript with contribution of all authors. All authors reviewed and approved the final version of the manuscript.