Large-Scale Variability of Physical and Biological Sea-Ice Properties in Polar Oceans

In this study, we present unique data collected with a Surface and Under-Ice Trawl (SUIT) during five campaigns between 2012 and 2017, covering the spring to summer and autumn transition in the Arctic Ocean, and the seasons of winter and summer in the Southern Ocean. The SUIT was equipped with a sensor array from which we retrieved: sea-ice thickness, the light field at the underside of sea ice, chlorophyll a concentration in the ice (in-ice chl a), and the salinity, temperature, and chl a concentration of the under-ice water. With an average trawl distance of about 2 km, and a global transect length of more than 117 km in both polar regions, the present work represents the first multi-seasonal habitat characterization based on kilometer-scale profiles. The present data highlight regional and seasonal patterns in sea-ice properties in the Polar Ocean. Light transmittance through Arctic sea ice reached almost 100% in summer, when the ice was thinner and melt ponds spread over the ice surface. However, the daily integrated amount of light under sea ice was maximum in spring. Compared to the Arctic, Antarctic sea-ice was thinner, snow depth was thicker, and sea-ice properties were more uniform between seasons. Light transmittance was low in winter with maximum transmittance of 73%. Despite thicker snow depth, the overall under-ice light was considerably higher during Antarctic summer than during Arctic summer. Spatial autocorrelation analysis shows that Arctic sea ice was characterized by larger floes compared to the Antarctic. In both Polar regions, the patch size of the transmittance followed the spatial variability of sea-ice thickness. In-ice chl a in the Arctic Ocean remained below 0.39 mg chl a m−2, whereas it exceeded 7 mg chl a m−2 during Antarctic winter, when water chl a concentrations remained below 1.5 mg chl a m−2, thus highlighting its potential as an important carbon source for overwintering organisms. The data analyzed in this study can improve large-scale physical and ecosystem models, habitat mapping studies and time series analyzed in the context of climate change effects and marine management.

In this study, we present unique data collected with a Surface and Under-Ice Trawl (SUIT) during five campaigns between 2012 and 2017, covering the spring to summer and autumn transition in the Arctic Ocean, and the seasons of winter and summer in the Southern Ocean. The SUIT was equipped with a sensor array from which we retrieved: sea-ice thickness, the light field at the underside of sea ice, chlorophyll a concentration in the ice (in-ice chl a), and the salinity, temperature, and chl a concentration of the under-ice water. With an average trawl distance of about 2 km, and a global transect length of more than 117 km in both polar regions, the present work represents the first multi-seasonal habitat characterization based on kilometer-scale profiles. The present data highlight regional and seasonal patterns in sea-ice properties in the Polar Ocean. Light transmittance through Arctic sea ice reached almost 100% in summer, when the ice was thinner and melt ponds spread over the ice surface. However, the daily integrated amount of light under sea ice was maximum in spring. Compared to the Arctic, Antarctic sea-ice was thinner, snow depth was thicker, and sea-ice properties were more uniform between seasons. Light transmittance was low in winter with maximum transmittance of 73%. Despite thicker snow depth, the overall under-ice light was considerably higher during Antarctic summer than during Arctic summer. Spatial autocorrelation analysis shows that Arctic sea ice was characterized by larger floes compared to the Antarctic. In both Polar regions, the patch size of the transmittance followed the spatial variability of sea-ice thickness. In-ice chl a in the Arctic Ocean remained below 0.39 mg chl a m −2 , whereas it exceeded 7 mg chl a m −2 during Antarctic winter, when water chl a concentrations remained below 1.5 mg chl a m −2 , thus highlighting its potential as an important carbon source for overwintering organisms. The data analyzed in this study can improve large-scale physical and ecosystem models, habitat mapping studies and time series analyzed in the context of climate change effects and marine management.

INTRODUCTION
Sea ice is one of the Earth system components most sensitive to climate change. Besides playing an essential role in global ocean circulation (e.g., Schmitz, 1995;Ferrari et al., 2014) and in regulating Earth's climate and weather (e.g., Liu, 2012;Dethloff et al., 2019), sea ice is crucial for the Arctic and Antarctic polar food webs (e.g., Eicken, 1992;McMinn et al., 2010;Meyer and Auerswald, 2014). Arctic and Antarctic pack ice serve as unique habitats for microalgae (Arrigo, 2014), which contribute to primary production and carbon transfer to higher trophic levels (Gradinger, 1999;Søreide et al., 2006;Budge et al., 2008;Fernández-Méndez et al., 2015;Wang et al., 2015Wang et al., , 2016Kohlbach et al., 2016Kohlbach et al., , 2017aKohlbach et al., ,b, 2018Schaafsma et al., 2017). Since light is the energy source of algae, the quantity and spectral composition of the light that penetrates the ice and reaches the surface ocean impacts primary productivity and biological activity at the bottom of the sea ice and below in the water column. Due to changes from complete darkness to 24 h daylight, and extreme natural variations in water temperature in high latitudes, the concentration of primary production in both the sea ice and the water column undergoes a seasonal cycle, which has a major impact on food availability and life cycle of animals living in polar regions (Swadling et al., 1997;Brierley and Thomas, 2002). Small organisms living inside or under the sea ice, and feeding on residing algae in (ice algae) or below (phytoplankton) it, transfer both ice-algae-produced and phytoplankton-produced carbon to the pelagic food web (Budge et al., 2008;Wang et al., 2015Wang et al., , 2016Kohlbach et al., 2016Kohlbach et al., , 2017aKohlbach et al., ,b, 2018. Thus, knowledge of sea ice and under-ice environmental properties as drivers of large-scale patterns in the abundance and distribution of seaice algae and phytoplankton, and, consequently, of zooplankton and nekton, is essential for understanding ecosystem functioning and predicting possible consequences of climate change for the ecosystems (Flores et al., 2011David et al., 2015David et al., , 2016Ehrlich et al., 2020).
Sampling in polar regions is limited in space and time because of the harsh weather during most of the year and because of the remote and hard to access locations. Particularly for sea-ice algae, for which chlorophyll a (in the following abbreviated as chl a) is usually used as a proxy for their biomass, estimates are often based on a small number of ice core observations, which do not fully capture the patchiness and the spatial and temporal variability in their distribution (Miller et al., 2015;Lange et al., 2016Lange et al., , 2017b. Satellite remote sensing during the past decades has vastly improved, enabling large-scale monitoring of certain sea-ice parameters, such as sea-ice concentration and thickness. However, uncertainties of global satellite sea-ice concentration and thickness retrievals on smaller scales, e.g., in the range of a few kilometers, can reach the magnitude of the measurement itself (Ivanova et al., 2014;Ricker et al., 2015). Moreover, despite new satellite advancements (i.e., ALOS-Advanced Land Observation Satellite), resolving small scale features such as ridges is still a challenge due to the relatively large footprint. Satellite snow depth retrievals are under development, but further evaluation and validation are needed (Guerreiro et al., 2016). Furthermore, satellite measurements do not provide enough information for the study of in-ice and under-ice environmental characteristics. This is particularly true for ice algae and phytoplankton in ice-covered regions that cannot be observed by satellite so that a comprehensive picture of their distribution on large scales remains difficult to obtain. Therefore, non-disruptive time series of sea-ice algae are non-existent. Only recently, technological developments allowed the estimation of under-ice properties, including light and inice chl a, on the scale of hundreds of meters in the Arctic region (Nicolaus et al., 2012;Nicolaus and Katlein, 2013;Katlein et al., 2015Katlein et al., , 2019Lange et al., 2017b) and in the Antarctic .
The effects of ice algae on the spectral distribution of light under the ice were first observed by Maykut and Grenfell (1975), followed by studies focused on investigating the wavelengths that are affected by the presence of algae (Legendre and Gosselin, 1991), and on distinguishing between effects of algae and snow on the spectral distribution of under-ice light (Perovich, 1990). Mundy et al. (2007) used, for the first time, normalized difference indices (NDI) for studying the effects of ice algae and snow on the under-ice hyperspectral measurements conducted under landfast first-year ice in the Canadian Arctic Archipelago. After comparing different methods for in-ice chl a retrieval, the NDI method was employed by Melbourne-Thomas et al. (2015) to estimate ice-algae biomass in two Antarctic regions, the Weddell Sea and off East Antarctica, and to retrieve inice chl a in the Weddell Sea in winter . Lange et al. (2016) developed and compared different algorithms for the retrieval of ice algae in the Arctic pack ice during summer 2012, and applied such algorithms to under-ice hyperspectral measurements collected with under-ice profiling platforms . For the first time, in this work these methods are applied to multi-annual data sets of under-ice hyperspectral measurements conducted over hundreds to thousands of meters, providing a unique large-scale study of the spatial variability of Arctic and Antarctic sea-ice algae biomass estimates.
In the present study, we present data collected with a Surface and Under Ice Trawl (SUIT) and used these data to characterize ice-associated environments in the Arctic Ocean and the Southern Ocean. We aim to provide a regional, and seasonal comparison of sea-ice properties on large scales in both Polar regions with a focus on: (1) under-ice water properties, (2) ice thickness, (3) light transmittance through sea ice and under-ice light, and (4) in-ice chl a estimates. Moreover, the spatial scale covered by the present data sets also allow for an (5) unprecedented investigation of the scales of variability for all the variables measured. Our account of the meter-to kilometer variability of structural and optical sea-ice properties should contribute to the validation and parameterizations of sea-ice models as well as habitat mapping in the polar oceans. and Table 1). In the Arctic Ocean, data were collected during early springtime (PS92, May-June 2015, Peeken, 2016), late springtime (PS106, June-July 2017, Macke and Flores, 2018), and during summer-autumn (PS80, August-September 2012, Boetius, 2013). In the Southern Ocean, the campaigns cover wintertime (PS81, August-October 2013, Meyer and Auerswald, 2014), and summertime (PS89, December 2014-January 2015, Boebel, 2015).
Sampling was performed with Surface and Under Ice Trawls (SUIT, van Franeker et al., 2009;Flores et al., 2012). The SUIT consists of a steel frame with a 2 × 2 m opening and two 15 m long nets: a 7 mm half-mesh commercial shrimp net, and a zooplankton net with a 0.15 mm (PS106) or 0.30 mm (PS80, PS81, PS89, PS92) mesh. Floats attached to the frame keep the net at the surface or right at the sea-ice underside. An asymmetric brindle forces the net to tow off at an angle of approximately 60 • so that it samples undisturbed sea ice away from the ship's wake. Previous studies using the SUIT nets have described in detail the ecological and biological aspects of the catch and distribution of organisms in both polar regions (Flores et al.,  , 2012David et al., 2015David et al., , 2016David et al., , 2017Schaafsma et al., 2016Schaafsma et al., , 2017. Since 2012 the SUIT is equipped with a sensors array Lange, 2017) developed in order to conduct coincident quantitative observations of the sea-ice and under-ice water environments. The sensors allow the direct observations of: (1) water inflow speed and direction, pitch and roll angles, and pressure (i.e., depth) using an Acoustic Doppler Current Profiler (ADCP; Nortek Aquadopp R Profiler) with three acoustic beams, which allows 3-dimensional measurements of current velocities, at a frequency of 2 MHZ, and a sampling interval of 1 s; (2) water temperature, water salinity (practical salinity scale PSS-78;Fofonoff, 1985), and water depth by using a Conductivity Temperature Depth (CTD) probe (Sea and Sun Technology CTD75M memory probe) with a sampling interval of 0.1 s; (3) under-ice water chl a concentrations using a fluorometer (Cyclops, Turner Designs, USA) incorporated into the CTD; (4) distance from the ice surface by using an altimeter (Tritech PA500/6-E) incorporated into the CTD; (5) under-ice  Table with listed, for each expedition, the number of profiles considered in the analysis (N hauls ), the total profile length, the mean ice concentration (Ā), mean total ice thickness (H i ) and mean snow depth (H s ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, temperature and salinity. light levels using Ramses spectral radiometers (Trios GmbH, Rastede, Germany) with a wavelength range from 350 to 920 nm and a resolution of 3.3 nm. Incident solar radiation and under-ice irradiance were measured using an irradiance sensor (RAMSES-ACC) containing a cosine receptor with a 180 • fieldof-view. Under-ice radiance measurements were acquired using a radiance sensor (RAMSES-ARC) with a 9 • field of view. As explained in detail in Lange (2017), the combined information from these different sensors allow also the retrieval of: (1) seaice draft derived by combining the CTD depth measurements with the distance from the ice (altimeter), and then corrected with pitch and roll measurements from the ADCP, as described in section 2.2; (2) ice-algal chl a derived from the Ramses spectral radiometers as described in section 2.4 (see also Table S1). Snow depth was recorded by visual observation of the ice under which the SUIT was traveling: A marked stick extending from the starboard side of the ship helps in quantifying the thickness of snow on top of the ice floes that are tilted during the ship passage. Part of these data have been used in previous works to study the relationship of environmental properties of seaice habitats with the community structure of sympagic fauna in the Arctic Ocean (PS80 and PS92, David et al., 2015David et al., , 2016Schaafsma, 2018;Flores et al., 2019;Ehrlich et al., 2020), and in the Southern Ocean (PS81, Schaafsma et al., 2016Schaafsma et al., , 2017David et al., 2017), to develop and test algorithms for estimating inice chl a (PS80, Lange et al., 2016;Lange, 2017), to retrieve Arctic primary production on large scales in the Arctic summer ( PS80 Lange, 2017;Lange et al., 2017b), and to estimate Arctic under-ice primary production based on under-ice irradiance measurements (PS92, Massicotte et al., 2019). All data collected during PS89 and PS106 are so far unpublished. Although data were collected during different years, this data set allows the investigation of differences between seasons in both the Arctic and the Antarctic, assuming that the sampling is representative of the season. The profile lengths used for calculation of mean quantities (e.g., ice and snow thickness, water temperature and salinity, chl a in water and sea ice) correspond to the distance between the start and end trawl points, excluding parts of the profiles where the sensors data were unreliable due, for example, to sensors failure or damage. Thus, the profile length does not always correspond to the trawled distance, i.e., the total distance during which the net was towed in water by the ship, used to calculate the density of animals (e.g., David et al., 2017).

2011
The sampled regions were different between expeditions: during Arctic spring (PS92 and PS106) data were collected north of Svalbard (Figure 1), in the latitudinal band between 80 • N and 85 • N. PS92 stations were located on the Yermak Plateau (lon < 10 • E, later referred to as Yermak stations), in the Sophia Basin and along the Svalbard shelf/slope (lon > 10 • E, later referred to as basin & shelf stations). Some PS106 stations were also located in the Basin & shelf area (10 • E < lon < 20 • E), the others were located between 20 • E and 30 • E. During Arctic summer (PS80), the sampling region extended toward the North Pole and to 130 • E, thus covering the Nansen Basin and the Amundsen Basin. The two Southern Ocean expeditions (Figure 2) took place in the Weddell Sea, but at different latitudinal locations. During winter (PS81), data were collected between 52 • S and 61.5 • S. During summer (PS89), the Marginal Ice Zone (MIZ) had already retreated, forming a belt relatively close to the continent, and facilitating the sampling of southern areas (between 66 • S and 69 • S). During the southward trip in December 2014, the MIZ extended to ∼66 • S, during the return trip northward in January the MIZ had retreated to ∼68 • S. We will refer to the two different MIZs as Dec-MIZ and Jan-MIZ. In the period in-between the sampling of the two MIZs, the ship traveled through the inner pack-ice area that we will refer to as the pack-ice area.

Draft Calculation
Retrieval of the sea-ice draft makes use of the combination of depth h w given by the pressure sensor of the CTD, distance from the ice h a computed from the altimeter, and pitch β and roll φ movement of the SUIT: The CTD and the altimeter are connected but mounted in different parts of the SUIT, thus they are located at different depths. h CTD is the distance between CTD and altimeter, α is the angle formed between the CTD probe, the CTD sensor, and the altimeter (see Figure 3 in Lange, 2017). In some cases, due to the failure of one instrument, the entire set of information was not available. In this case, we used the simplified version of equation (1) that does not include the pitch and roll correction: (2) FIGURE 3 | Correlation surfaces (A) of normalized difference indices for integrated chl a for PS92 and PS106. λ 1 and λ 2 are wavelengths pairs. Relationship (B) between observed normalized difference indices and integrated chl a (dots). The green line represents the fitted relationship given by the algorithm c in Table 2.
This is equivalent of assuming that the SUIT is towed perfectly parallel to the ice and, even if only an approximation, it has been proven reliable (R 2 = 0.78) in determining draft when the entire set of information was not available (Lange, 2017;Lange et al., 2017b). Handling of data gaps due to sensors malfunctioning or failure is explained in the Supplementary Material. Ice concentration for each profile was computed as the percentage of data points where ice thickness H i was larger than zero. Ice thickness was computed by using a fixed density value (ρ = 0.917 g cm −3 ) for the sea ice. By doing this, final ice thickness values include both sea ice and snow, and we will refer to it as total thickness in the following . Mean and median thicknesses were computed for H i > 0 and they include ridges. Keels of ridges were detected along each profile by using the Rayleigh criterion (Rabenstein et al., 2010;Castellani et al., 2014Castellani et al., , 2015. A detailed analysis of the distribution and properties of ridges detected along each SUIT profile will be provided in a separate study. During all expeditions, a helicopter-borne frequency-domain electromagnetic induction sounding system (EM-Bird) was used to measure the total sea-ice thickness (sea-ice thickness plus snow depth, Haas et al., 2009). The 4-m long EM-Bird was towed at a height of 10-15 m above the surface under a helicopter along mainly triangular flight tracks. With this method, the larger scale sea-ice thickness distribution can be observed also away from the ship track, which usually follows the easier seaice conditions in a region. Each triangle leg covered 15-30 nautical miles.

Under-Ice Light
Irradiance and radiance values were integrated over the Photosynthetically Active Radiation wavelength range (PAR; 400-700 nm). Since the SUIT does not always travel parallel to the ice but can encounter obstacles that make it swing, we included in the analysis only data with an inclination angle lower than 15 • . In order to reduce the effect of underice water between the actual position of the sensors mounted in the SUIT frame and the underneath side of the ice, we excluded all measurements associated with an altimeter value larger than 1.5 m following works by Katlein et al. (2015Katlein et al. ( , 2016Katlein et al. ( , 2017, Lange et al. (2016Lange et al. ( , 2017b. Transmittance was calculated as the ratio between underice and incoming radiation, the latter measured with a Ramses sensor mounted on the ship's crows nest. Due to a sensor failure during summer 2012 (PS80) we used incoming global radiation data, which are obtained with a pyanometer, also placed in the same region as the above described sensor. For comparison, a regression between the two sensors was calculated for the simultaneous measuring period and gave excellent agreement. Thus, data from this were used to calculate transmittance values in September 2012. A detailed method used to obtain radiation values to cover this period can be found in the Supplementary Material.
The quantity of light reaching the underside of the ice layer is important for sea-ice algae and phytoplankton, as well as for sympagic fauna. In order to compare the amount of light between expeditions and sampling stations independently from the time of the day when we sampled, we calculated the insolation parameter as: where S i is the insolation value (mol photons m −2 d −1 ) at position i along each profile, TR i is the transmittance value at position i, and I h is the modeled hourly incoming solar radiation (µmol photons m −2 s −1 ) above the surface (1 • x1 • spatial resolution, daily temporal resolution, interpolated hourly) based on the radiative transfer model SBDART (Ricchiazzi et al., 1998) as described in Laliberté et al. (2016). These data do not consider atmospheric parameters, such as cloudiness, which can affect the amount of radiation reaching the ship, where the incoming irradiance sensor is usually mounted. So defined, the insolation S i is the amount of light passing through the bottom of the sea ice (i.e., PAR available for ice algae), at a specific position along the profile over an entire day. It is important to notice that the absence of information on cloud cover, or any other atmospheric parameter that could dampen the amount of light reaching the ice surface, potentially lead to overestimating the actual insolation values.

Retrieval of In-Ice Chlorophyll a
Since the retrieval of extensive spatial and temporal observations of sea-ice algae is limited in large part due to destructive techniques (i.e., based on the collection of sea ice from the field and subsequent melting and filtering of the ice), many studies in the past decades focused on estimating sea-ice algae based on under-ice hyperspectral measurements (Perovich, 1996;Mundy et al., 2007;Melbourne-Thomas et al., 2015;Lange et al., 2016;Meiners et al., 2017;Wongpan et al., 2018). In the present study, we employed normalized difference indices (NDIs) and Empirical Orthogonal Functions (EOFs) of under-ice spectra to estimate the concentration of chl a present in the ice. The heterogeneity of the environment covered by the present data set required the development of a specific algorithm for each expedition. The algorithms applied to the different data sets are listed in Table 2. For Arctic summer (PS80) data we used the EOF algorithm developed by Lange et al. (2016) based, amongst others, on the PS80 data set. Lange et al. (2016) recommended the EOF algorithm, instead of NDI, because it proved to be most reliable in catching the variability of in-ice chl a in the PS80 data set. For the Southern Ocean (PS81 and PS89) we used the NDI algorithm derived by Melbourne-Thomas et al. (2015) for the Weddell Sea. This algorithm was found to provide the most robust predictions of integrated chl a in comparison to others tested by Melbourne-Thomas et al. (2015). The NDI algorithm was developed based on summer data, and it is thus applicable to Antarctic summer (PS89). Meiners et al. (2017) applied the same algorithm to under-ice hyperspectral data collected with an ROV deployed during the Antarctic winter (PS81). Results showed that the algorithm provides a good fit with ice-core data collected during ice stations, we thus applied the same algorithm to the SUIT winter data (PS81). The NDI algorithm for Arctic spring (PS92 and PS106) was developed by comparison between coincident under-ice hyperspectral and in-ice chl a measurements of in total 8 ice stations carried out during the same campaigns. The accordance in sampling time (May-June-July), sampling region (Figure 1), and similarity of sea-ice properties (see section 3.2.1) justifies merging the data from these two spring expeditions. We used a total of 19 coincident measurements of under-ice irradiance and integrated chl a retrieved by measuring autotrophic pigments with High-Performance Liquid Chromatography (HPLC, for more details, see Tran et al., 2013). The NDIs for each wavelength pair were correlated with the integrated chl a values. Values lower than 0.1 mg m −2 were excluded from the analysis. Correlation surfaces of normalized difference indices are shown in Figure 3A. We then used a linear model to explore the relationship between predicted chl a from the NDI algorithm and integrated chl a from the ice cores ( Figure 3B). The NDI algorithm developed explains 60% of the variability (R 2 = 0.60). An example of a SUIT profile including CTD data, retrieved ice draft, and retrieved in-ice chl a is shown in Figure 4. The sampling frequency of the RAMSES radiometers was lower than the one of the CTD and ADCP. In order to match ice thickness/draft with under-ice hyperspectral measurements, we computed averages of thickness and draft values falling between the beginning and the end of each light measurement.

Statistical Analysis
We applied Principal Component Analysis (PCA) to assess spatial and seasonal patterns in the variability of physical properties of the sea ice and the underlying water sampled by the SUIT separately for each hemisphere. We included in this analysis the parameters water temperature, salinity, total ice thickness from the present study, and ridge depth and ridge density (Rabenstein et al., 2010;Castellani et al., 2014. Near-normal distribution of data, as assumed by PCA, was checked by visual inspection of histograms. In order to achieve near-normal distribution of the data, total ice thickness and ridge density were square-root transformed. To test for significant differences between the frequency distributions of total ice thickness data from the EM-bird and the SUIT, we applied the Kolmogorov-Smirnov test. Before applying the test, the SUIT data were resampled to the footprint of the EM-bird (40 m) by averaging the data for each 40 m segment FIGURE 4 | Example of a SUIT profile (station PS92-47_1). Orange and red lines represent surface salinity and temperature, respectively. The gray line is the draft obtained with Equation (1). The blue line is the draft resampled at the frequency of the RAMSES radiometers sampling (section 2.4). The green circles represent in-ice chl a retrieved by applying the NDI algorithm c in Table 2. of each thickness profile. All statistical analyses were performed with the software R version 3.5.2 (R-Development-Core-Team, 2018), applying the package vegan (Oksanen et al., 2013).

Spatial Autocorrelation Analysis
Spatial autocorrelation was used to investigate the horizontal patchiness of sea-ice draft for 55 stations. Spatial autocorrelation analyses were also conducted for transmittance and in-ice chl a biomass. However, the transect nature of the surveys, the large interval spacing and small sample sizes, in comparison to draft, limited the spatial analyses of in-ice chl a to only 7 stations and of transmittance to only 10 stations. Autocorrelation was estimated using Moran's I (Moran, 1950;Legendre and Fortin, 1989;Legendre and Legendre, 1998), which was calculated for each SUIT survey at equally spaced (25 m) distance classes. Individual autocorrelation coefficients (e.g., Moran's I estimates) were plotted for each distance class as a spatial correlogram (Legendre and Fortin, 1989;Legendre and Legendre, 1998) using the R software function correlog from the pgirmess package (Giraudoux, 2018). Autocorrelation coefficients for each distance class were assigned a two-sided p-value according to Legendre and Fortin (1989) and Legendre and Legendre (1998). The presence of spatial autocorrelation (i.e., patchiness) was determined if the correlogram was considered globally significant at p < 0.05. We used the first x-intercept of globally significant correlogram lines as an indicator of the patch size for sea ice draft, Pd (Legendre and Fortin, 1989;Legendre and Legendre, 1998). This methodology is consistent with spatial autocorrelation analyses conducted on ROV gridded data from one of the same cruises (PS80, Lange et al., 2017b) and in other snow and sea ice studies (e.g., Gosselin et al., 1986;Rysgaard et al., 2001;Granskog et al., 2005;Søgaard et al., 2010).

Sea Ice and Under-Ice Water Properties
Under-ice water properties for the three Arctic expeditions are presented in detail in Tables 3-5. In spring (PS92 and PS106), under-ice water temperatures were lower and salinities were higher compared to summer (PS80). The mean under-ice water chl a was highest in spring (6.34 ± 5.94 mg m −2 in May 2015, PS92) and lowest in autumn (1.60 ± 0.64 mg m −2 in September 2012, PS80). Besides seasonal differences, the data also show a regional pattern. During spring 2015 (PS92) there was a distinct difference between Yermak stations and the stations on the Sophia Basin & shelf ( Table 3). The latter had higher under-ice water chl a (with a mean of 12.52 ± 9.14 mg m −2 compared to 1.80 ± 2.2 mg m −2 for the Yermak stations), higher mean under-ice water temperature (−1.54 ± 0.13 • C compared to −1.77 ± 0.05 • C) and lower mean salinity (33.70 ± 0.43 compared to 34.13 ± 0.23). We could not confirm if this pattern was repeated on PS106 because the Yermak Plateau was not reached. However, during this expedition higher temperatures, lower salinities and very high chl a values (up to 20 mg m −2 ) in the surface layer were associated with the shelf-slope area toward the end of the expedition, and were probably related to the beginning of ice breakup ( Table 4). During spring 2017 (PS106), time progression corresponded to a decrease in temperature and salinity and an increase in under-ice water chl a (PS106 , Table 4). During summer 2012 (PS80), the stations situated in the Nansen Basin were characterized by higher salinities (mean 31.85 ± 0.94) and lower surface chl a (mean 0.66 ± 0.28 mg m −2 ), whereas stations in the Amundsen Basin had lower mean salinities (31.02 ± 1.71) and higher mean chl a concentrations (1.46 ± 0.60 mg m −2 , Table 5, David et al., 2015).  Table for PS92 with, for each station (stn), the total profile length, the mean ice concentration (Ā) along profile, mean snow depth (H s ) and mean total ice thickness (H i ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, under-ice water temperature and salinity, ice concentration retrieved from satellite (Ā Sat ), median in-ice chl a (interquartile range), mean transmittance TR (± one standard deviation), mean Insolation S i (± one standard deviation), and the draft patch size (P d ).   (stn), the total profile length, the mean ice concentration (Ā) along profile, mean snow depth (H s ) and mean total ice thickness (H i ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, under-ice water temperature and salinity, ice concentration retrieved from satellite (Ā Sat ), and median in-ice chl a (interquartile range), mean transmittance TR (± one standard deviation), mean Insolation S i (± one standard deviation), and the draft patch size (P d ).
stn   (stn), the total profile length, the mean ice concentration (Ā) along profile, mean snow depth (H s ) and mean total ice thickness (H i ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, under-ice water temperature and salinity, ice concentration retrieved from satellite (Ā Sat ), and median in-ice chl a (interquartile range), mean transmittance TR (± one standard deviation), mean Insolation S i (± one standard deviation), and the draft patch size (P d ).
stn In gray shades are the stations located in the Nansen Basin whereas the white ones are located in the Amundsen Basin. During the springtime campaigns in 2015 (PS92) and 2017 (PS106), Arctic sea ice was covered with snow, whereas during summer 2012 (PS80) the ice surface was covered by melt ponds. Despite the difference in timing, snow depth decreased during both PS92 and PS106 and melt ponds started to form at the end of both expeditions (Tables 3, 4). Toward the end of summer 2012 (PS80), refreeze started and snow accumulated over the ice (Table 5). Modal sea-ice thickness (Figure 5) was the highest in the Arctic spring, with modes of 1.25 m and 1.65 m in 2015 (PS92) and 2017 (PS106), respectively, and lowest during the summer campaign in 2012 (PS80), with a primary mode for very thin ice (<10 cm). With an explained variability of 81.9% for the first two principal components, the outcome of the PCA (Figure 6A) confirms a high similarity of the two springtime data sets (PS92, PS106), and a distinctly different combination of physical sea-ice and under-ice properties in most stations sampled during summertime (PS80). Salinity and ridge properties were the main drivers of the variability in sea-ice and under-ice properties ( Figure 6A). The analysis of regional patterns shows that total ice thickness decreased eastwards with a decrease of modal thickness from the Yermak stations EM-bird data ( Figure 5) show a similar pattern as the SUIT data, with a decrease in mean total ice thickness between spring and summer from 1.8 and 1.44 m during spring 2015 (PS92) and 2017 (PS106), respectively, to 0.71 m in summer 2012 (PS80). The shape of the ice-thickness distributions in spring was broadly similar between EM-bird data and SUIT data. However, a statistical comparison of the two data sets indicated significant differences between the two methods in all three sampling seasons (Kolmogorov-Smirnov test, p ≪ 0.001), probably due to differences of several decimeters in the modes of the thickness distributions. Histograms of total thickness for both the SUIT and the EM-bird data are presented in the Supplementary Material.

Under-Ice Light and Light-Derived In-Ice Chlorophyll a Concentration
Values of transmittance show a mode for TR < 0.05 in both spring and summer. However, the summer mean value is larger (0.203 compared to 0.136 and 0.184 in 2015 and 2017, respectively). During summer 2012, high transmittance values were more frequent compared to spring 2015 and 2017. The similarity in the distribution of transmittance data during spring 2015 and 2017, shows that the two campaigns sampled the same variability in sea-ice and snow conditions (Figure 7). A monthly comparison shows that there was a change in mean transmittance from 0.22 ± 0.03 in May 2015 to 0.33 ± 0.3 in August 2012, then a decrease to 0.12 ± 0.19 in September 2012. The extinction coefficient k b calculated by non-linear regression from the May subset of the data showed the largest value with k b = 8.678 m −1 . In June and July the extinction coefficients were lower, with values of 2.209 and 1.679 m −1 , respectively, reaching the minimum at 1.576 m −1 in August (PS80). The extinction coefficient increased again in September to 6.552 m −1 . Insolation values in the Arctic (Figure 7)    During PS92 the stations on the Yermak Plateau had the lowest in-ice chl a values, whereas the Basin & shelf stations had the highest ( Table 3). High chl a concentrations were also found in the Basin & shelf area during PS106, but in this case the spatial variability was less pronounced ( Table 4). During summer 2012 (PS80), the Nansen Basin presented lower values of in-ice chl a compared to the Amundsen Basin. In the area of overlap between PS80 and PS106 stations, chl a values from the two campaigns were very similar, varying around 0.2 mg chl a m −2 . Distributions of in-ice chl a are similar for the two spring campaigns, whereas the distribution of values from the summer campaign shows significantly higher values ( Figure 7C). Histograms of in-ice chl a are presented in the Supplementary Material.

Sea Ice and Under-Ice Water Properties
During the Antarctic winter in 2013 (PS81), under-ice water temperatures were always close to the freezing point (Tables 1, 6) and uniform between sampling stations with a mean temperature of −1.85 ± 0.01 • C. During summertime 2014/15 (PS89 , Table 7), mean surface temperatures ranged between −1.67 ± 0.12 • C in the Dec-MIZ and −1.58 ± 0.12 • C in the Jan-MIZ. Lowest under-ice water temperatures characterized the pack-ice area sampled in December, with a mean of −1.86 ± 0.01 • C. Likewise, the mean chl a concentration in the under-ice water was lower in winter 2013 (PS81, mean 0.42 ± 0.12 mg m −2 ) than in summer 2014/15 (PS89, mean 0.88 ± 0.62 mg m −2 ), and mean salinity was higher (34.22 ± 0.12 psu) in winter 2013  Table for PS81 with, for each station (stn), the total profile length, the mean ice concentration (Ā) along profile, mean snow depth (H s ) and mean total ice thickness (H i ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, under-ice water temperature and salinity, ice concentration retrieved from satellite (Ā Sat ), and median in-ice chl a (interquartile range), mean transmittance TR (± one standard deviation), mean Insolation S i (± one standard deviation), and the draft patch size (P d ). In gray shades are the westernmost stations, in white are the easternmost stations. *For station 555_1 there is only one value for in-ice chl a and we report it here only to help data interpretation.  (stn), the total profile length, the mean ice concentration (Ā) along profile, mean snow depth (H s ) and mean total ice thickness (H i ), under-ice water chl a (chl a w ) multiplied by the vertical section of the SUIT (2 m) to obtain integrated values, under-ice water temperature and salinity, ice concentration retrieved from satellite (Ā Sat ), and median in-ice chl a (interquartile range), mean transmittance TR (± one standard deviation), mean Insolation S i (± one standard deviation), and the draft patch size (P d ).
stn In light gray shades are the stations located in the Marginal Ice Zone in December (referred to as Dec-MIZ), in white are the stations in the pack-ice area, whereas in dark gray are the stations located in the Marginal Ice Zone in January (referred to as Jan-MIZ). *Values for in-ice chl a at stations 70_2 and 71_1 were retrieved even if the altimeter values were missing (so it was not possible to filter the data) and are given to help the data interpretation in section 4.
(PS81) than in summer 2014/15 (33.63 ± 0.24 psu). Between December (PS89, Dec-MIZ) and January (PS89, Jan-MIZ), both under-ice water temperature and chl a decreased, whereas salinity did not show any pattern ( Table 7). The pack-ice area had different water properties compared to the MIZs, and it was characterized by significantly colder temperatures, more saline waters, and higher surface chl a concentrations ( Table 7). Sea ice in the Southern Ocean was covered with snow during both winter 2013 (PS81, Table 6) and summer 2014/15 (PS89, Table 7). The snow depth ranged from 5 to 60 cm in winter 2013 (PS89) and from 5 and 100 cm in summer 2014/15 (PS89). Modal ice thickness was lowest in winter 2013 (PS81, 0.45 m) and highest in summer 2014/15 (PS89, 0.75 m) in the Jan-MIZ (Figure 5). During both expeditions, the differences in total ice thickness varied from thinner ice in the northern latitudes to thicker ice in southern parts of the sampling areas. In PS89, the shape of the thickness frequency distribution in the pack-ice zone was different from the ones of the MIZs with a longer tail toward larger thickness values (data not shown). Ice thickness distribution from the EM-bird data (Figure 5) was similar in shape to those of the SUIT for both expeditions. In EM-bird data, however, modal ice thicknesses were somewhat higher compared to SUIT data (Figure 5), with a mode of 0.65 m in winter (PS81), and a mode of 0.85 m in summer (PS89). As in the Arctic Ocean, a statistical comparison of the two data sets indicated significant differences between the two methods in both sampling seasons (Kolmogorov-Smirnov test, p < 0.001). Besides the differences in modal thickness, this difference was also related to greater variations in the sampling of very thin ice between the two methods, with a larger proportion of thinner ice in the SUIT data compared to EM-bird data.
The PCA biplot of the physical sea-ice and under-ice properties during the two Southern Ocean expeditions showed a gradual separation of the two sampling seasons, indicating that most stations sampled in winter (PS81) were associated with lower under-ice water temperatures and higher salinities than most stations sampled in summer (PS89) (Figure 6B). The cumulative explained variability of the two first principal components was 95.4%. Within each expedition, the PCA ordination reflected the differences between sampling regions ( Figure 6B).  , Table 7).

Under-Ice Light and Derived Chlorophyll a Concentration
During winter 2013, in-ice chl a was significantly higher than during summer, with a mean of 7.91 ± 1.07 mg chl a m −2 in winter compared 0.88 ± 0.51 mg chl a m −2 in summer ( Figure 7C). Principal modes were at 8.5 and 1.3 mg chl a m −2 during winter and summer, respectively. The variability within expeditions was low in both seasons. During summer 2014/15 higher in-ice chl a values were measured in the Jan-MIZ, compared to the Dec-MIZ.

Spatial Autocorrelation Analysis
The summary of the spatial autocorrelation analysis and the results for the patch size in both the Arctic and the Antarctic are presented in Table 8. Spatial autocorrelation analysis resulted in significant sea ice draft patch size estimates for 13 of the 19 SUIT stations in the Southern Ocean and 42 of the 48 SUIT stations in the Arctic Ocean, with significant results for almost all cruises. Patch size, which for ice draft can be interpreted as the size of smooth sea-ice areas, varied from 50 to 750 m with a mean of 221±156 m (Table 8) for the Arctic Ocean. In the Southern Ocean, autocorrelation analysis showed a draft patch size varying from 50 to 400 m (with an average of 129±92 m, Table 8), with larger values for the summer expedition (PS89). The spatial autocorrelation analysis for transmittance (i.e., patch sizes) showed statistical significant results for SUIT in 3 of the 19 SUIT stations in Antarctica and in 7 of the 48 SUIT stations in the Arctic Ocean, with no statistically significant results during the summer Arctic cruise PS80. Our analysis showed a significant Pearson's correlation between the patch size of draft and transmittance (correlation coefficient = 0.77; N = 9; p = 0.016), indicating that sea-ice draft (i.e., total ice thickness) is the driver of variability in sea-ice transmittance. Correlograms for one SUIT haul for each expedition are shown in Figure 8 together with a comparison in patch size between the two Polar regions for both ice draft and transmittance.

Environmental Properties
This study presents the first characterization of sea ice and underice environments from multiple expeditions and seasons on the scale of kilometers in both ice-covered polar oceans. In the Arctic Ocean, ice and snow start to melt with the spring to summer progression and the increase in atmospheric and oceanic temperatures. The results showed this seasonal progression with a decrease of snow cover with subsequent formation of melt ponds, and a decrease of salinity in the under-ice water, an indicator of the presence of melt water under the ice. More light was found to penetrate the ice as time progressed, which can be attributed to a decrease in both ice thickness and snow cover (Tables 3-5). Similar patterns were seen in previous studies (e.g., Nicolaus et al., 2010). At some sampling stations, increased light penetration was accompanied by an increase of under-ice water chl a concentrations. During the two spring campaigns in 2015 and 2017 (PS92 and PS106), sampled locations were situated closely together with an overlap in the region 10 • E < lon < 20 • E, defined as the Basin & shelf region. When comparing the two expeditions, the spatial pattern appeared stronger than the inter-annual one. During the PS92 expedition, the Basin & shelf region showed a difference in water masses, characterized by more saline and warmer water, compared to the Yermak Plateau . The under-ice water chl a concentration was lower over the Yermak Plateau compared to the Basin & shelf region. In addition to these measurements, Assmy et al. (2017) found that the Basin & shelf area was characterized by large leads allowing light to penetrate through the ice which potentially explains increased chl a concentrations and a higher primary production compared to the Yermak region (Massicotte et al., 2019). The two regions also differed in the composition of under-ice fauna, and zooplankton abundances were generally very low in the under-ice water over the Yermak Plateau (Ehrlich et al., 2020). In summer 2012 (PS80), the Nansen vs. Amundsen Basin separation (Table 5) emerged clearly from the under-ice water salinity and chl a concentration, pointing to two different environmental regimes: an Atlantic influenced and nutrientrich regime in the Nansen Basin, and a nutrient-poor regime influenced from the Laptev Sea shelf in the Amundsen Basin (David et al., 2015;Flores et al., 2019). Surface temperature, however, showed a seasonal decrease toward winter. Ice thickness differences between the two basins (not shown here) remained very small pointing to similar sea-ice conditions in the two basins ( Figure 6A). However, the Amundsen Basin, which was sampled later in the season, had a principal mode at a thickness value of 0.65 m (figures not shown) and a secondary mode for very thin ice (values < 0.1), indicating that refreeze started already.
Additionally, under-ice water parameters differed between the two basins (Table 5), as well as the zooplankton community structure, which showed a dominance of copepods in the Nansen Basin, but a co-dominance of copepods and amphipods in the Amundsen Basin (David et al., 2015). The shape of the total sea-ice thickness distributions from the SUIT data agreed largely with those from the EM-measurements. In both, the Arctic Ocean and the Southern Ocean, however, median and modal total thicknesses from the EM-bird and the SUIT differed by up to several decimeters, causing significant differences between the two distributions with no systematic bias toward one instrument (Figure 5). These differences could be due to differences in the geographical area covered by the two devices, thus one of the instruments might have collected a larger amount of data points in a region with thicker/thinner ice compared to the other instrument. This was particularly the case in PS89, where the EM-measurements included areas covered by landfast ice, which were not sampled with the SUIT (data not shown). Differences between SUIT and EM-measurements were probably also a result of a balance between systematic biases of the two devices: On the one hand, SUIT-based total thickness estimates could have been negatively biased because the ship tends to travel in areas where navigation is easier (i.e., thinner ice; PS92, PS81, PS89). On the other hand, the lower spatial resolution of the EM-measurements causes an underrepresentation of deep ridges compared to the SUIT, shifting the distribution of total thicknesses toward the lower end (PS106, PS80). Moreover, the way snow is treated might introduce bias. The EM-bird directly measures snow plus ice thickness, on the other hand, for the SUIT, ice draft is converted into total ice thickness. Therefore, a change in the actual snow depth will have different impacts on SUIT and EM-bird ice thickness retrievals. When comparing with satellite data (columnĀ Sat in Tables 3-7), the SUIT can provide representative ice-concentration values. Only in the marginal ice zone, where the SUIT preferentially sampled only the ice-covered parts of the grid cell detected by satellites, greater differences occurred (David et al., 2015).
Transmittance values largely reflected sea ice and snow conditions. The greater ice thickness during spring 2017 (PS106) compared to 2015 (PS92) was associated with lower mean transmittance values (Tables 3, 4), whereas the absence of snow during summer 2012 (PS80) led to higher transmittance ( Table 5 and Figure 7). A larger spread of values during spring 2017 (PS106) compared to spring 2015 (PS92), however, points to higher variability of sea-ice and snow conditions, which was confirmed by Katlein et al. (2019), and which could be explained by the slightly different sampling period extending more toward summer during PS106. We estimated the insolation parameter (i.e., the daily integrated sunlight received) in order to obtain a representative comparison between SUIT hauls even if they were carried out at different times of the day. Mean insolation values in summer were higher than in spring, as a result of the larger transmittance values. However, the spread of values was wider in spring, which was related to incoming solar radiation that is maximum in the months of June-July compared to August-September (Arrigo et al., 2012;Arndt and Nicolaus, 2014). Extinction coefficients calculated from the present data set include the effect of both snow thickness and ice thickness, thus they represent bulk coefficients. Our summer values correspond to the ones generally known for bare ice (Grenfell and Maykut, 1977;Perovich, 1996;Light et al., 2008). The values presented here are in good agreement with those shown recently by Katlein et al. (2019), where bulk coefficients varied from ∼5 m −1 in May to ∼1 m −1 in August, and then increased again in September   Figure 8).
In the Southern Ocean, variations within winter 2013 (PS81) stations in salinity and temperature were rather small due to a quasi-homogeneous winter layer circulation in the Weddell Sea . The north-eastward sampling performed during the winter-spring transition led to a gradual decrease in sea-ice coverage, higher under-ice insolation, and an increase of under-ice water chl a by the end of the expedition, indicating that the productive season commenced . As for the Arctic expeditions, disentangling spatial and temporal patterns is challenging. Despite this difficulty, we can conclude that larger thickness values in the Jan-MIZ can be explained by the more southern locations. Total ice thickness distribution for winter 2013 (PS81) and summer 2014/15 (PS89) were similar, indicating similar sea-ice conditions. The cyclonic drift pattern of the Weddell Gyre moves sea ice from the eastern Weddell Sea/Lazarev Sea northwestward (Kwok et al., 2017) in a regime of mainly free drift. This decreases the tendency of sea ice to deform (G. Castellani, unpublished data) and it leads to more uniform sea-ice conditions in that area.
The similarity in snow cover and ice thickness between winter 2013 (PS81) and summer 2014/15 (PS89) leads to comparable transmittance values, even though the maximum transmittance values were lower in winter than in summer. Differently from the Arctic, Antarctic sea ice remained covered by relatively thick snow also during summer, explaining why transmittance values never approached 1.

In-Ice Chl a
We applied NDI and EOF algorithms to retrieve in-ice chl a from under-ice irradiance and transmittance spectra. We used one algorithm for Arctic summer , one for the Antarctic regions (Melbourne-Thomas et al., 2015), and one for the Arctic spring. In all the three cases, the retrieved values of inice chl a were in good agreement with ice-cores data collected during PS80 (Fernández-Méndez et al., 2015), and during PS92 and PS106 (Ehrlich et al., 2020, I. Peeken, unpublished data), and with derived in-ice chl a based on ROV measurements during PS81 .
In the Arctic, the progression in snowmelt is visible in the increase of in-ice chl a with seasonal progression. Spring campaigns were conducted from late May onward, when incoming light levels were already high enough to cause an algal bloom, but snow cover obviously hampered this. Ice-algae growth was limited by the snow cover and the increase of in-ice chl a in the last stations during both PS92 and PS106 reflected the better light conditions after the onset of snowmelt. In-ice chl a concentrations were higher toward the end of summer 2012 (PS80) than in the two Arctic spring studies. This indicates that the ice-algae bloom only commenced after the onset of snowmelt, as expected from modeling studies (Castellani et al., 2017). During our sampling of PS80, large parts of the ice-algae bloom had already been released to the sea floor due to basal melting , indicating that the observed in-ice chl a values had probably been considerably higher during the peak of the bloom.
In-ice chl a in the Antarctic winter was much higher than on any other expedition in the Arctic Ocean or the Southern Ocean. Our winter values agree well with values obtained from ice cores collected during the same expedition . The presence of high biomass despite very low transmittance and insolation values indicate that the algae present in winter are the result of a previous accumulation of high biomass along the entire core length, more than an active growth at the bottom (Meiners et al., 2012). Winter in-ice algal assemblages often consist mainly of diatoms (Garrison and Close, 1993;Ugalde et al., 2016). During PS81 this was confirmed by sea-ice fatty acid compositions (Kohlbach et al., 2017b), and algae found in the stomach of larval krill . Notably, in-ice chl a mean values in the Antarctic in winter were always higher than phytoplankton chl a values from the under-ice water. This was valid also for some stations during summer (i.e., stations 30_4, 37_2, 62_1, 70_1, and 71_1). Several Antarctic species are attracted to the sea ice during both summer and winter, including Antarctic krill Euphasia superba (Flores et al., 2011(Flores et al., , 2012(Flores et al., , 2014, showing the importance of sea ice-derived carbon for the Southern Ocean food web (Jia et al., 2016;Kohlbach et al., 2017bKohlbach et al., , 2018Kohlbach et al., , 2019. This basal role of sea-ice derived carbon has a major impact on the distribution, abundance, life-cycle and possibly survival of many organisms residing in ice-covered oceans.

Physical and Biological Comparison
The present data set agrees with the well-known characteristics of the two polar regions, thus proves to be representative of the environment sampled. Variability in sea-ice properties was higher in the Arctic Ocean compared to the Southern Ocean. Furthermore, our spatial autocorrelation analysis of total sea-ice thickness showed smaller patch sizes in the Southern Ocean than in the Arctic Ocean. Interpreting this as bottom roughness, we can say that the Antarctic sea ice was characterized by smaller floes, or smaller patches of level ice compared to the Arctic, pointing to the different sea-ice drift, growth, and deformation regimes (Haas, 2008;Massom and Stammerjohn, 2010). Additionally, the correlation between total thickness and transmittance patch sizes indicates that, at floe-scales of around 120 m for the Southern Ocean and around 220 m for the Arctic Ocean, total ice thickness is a good predictor for transmittance, which is in agreement with previous work from the Arctic (e.g., Katlein et al., 2015). At smaller spatial scales, typically around 10 m, Katlein et al. (2015) found that surface properties such as melt ponds (during Arctic summer) were the best predictor for transmittance, which was not possible to resolve with the present data set due to the large spacing intervals between subsequent under-ice irradiance measurements (>10 m). The snow remains thick in the Antarctic during summer, whereas during Arctic summer snow melts completely and melt ponds form on the sea ice surface (Eicken, 2003;Nicolaus et al., 2012;Webster et al., 2018). Despite higher transmittance, insolation values (i.e., the amount of daily light that an algae receives at the bottom of the sea ice) were higher in the Arctic spring than in summer. The range of transmittance values during PS80 was higher compared to values collected in the same Arctic region and in the same season, but one year earlier, with a Remotely Operated Vehicle (ROV, Nicolaus and Katlein, 2013) diving under the ice (Nicolaus et al., 2012), and compared to those collected in July 2014 with an AUV diving under an ice floe . The same holds for the transmittance values for PS81 and PS89: the range of transmittance variability measured with the SUIT was higher, and with higher mean values, than those obtained during the same expeditions with the ROV (PS81 data were published in Arndt et al., 2017, PS89 are unpublished data). Differences arise mainly from sampling methodologies. Whereas, with the ROV data are collected under a single, relatively uniform ice floe, the SUIT travels through a more heterogeneous environment, including also very thin ice or brash ice present in between ice floes. For this reason, the different methods sample differently and transmittance values can be higher when calculated from the SUIT data set compared to the ROV data set (Massicotte et al., 2019). Furthermore, a towed approach such as the SUIT is more robust compared to robotic platforms (e.g., ROVs and AUVs) and is thus able to ride directly along the underside of the sea ice at large distances with minimal concern about damage and losing the platform. However, the towed SUIT is subject to minimum speed requirements in order to maintain momentum while breaking sea ice and to representatively catch under-ice fauna, which is a main objective of SUIT deployments. Thus, ROVs and AUVs can survey the under ice environment at a more controlled and lower speed providing a greater spatial resolution, however, at the cost of spatial coverage.
As expected, insolation values were very low in winter (PS81). Despite the highest snow cover and the low transmittance values, during the Antarctic summer (PS89) the under-ice insolation was greater than in the Arctic summer (PS80). This difference reflected the different latitudinal ranges of the two sampling areas. The difference in maximum latitude (70 • S compared to 85 • N) and the difference in solar inclination lead to stronger incoming radiation during the Antarctic summer (PS89) than during the Arctic summer (PS80), leading to stronger underice insolation. Surface chl a values were higher in the Arctic (particularly in spring) than in the Antarctic. In contrast, inice chl a values were higher in Antarctic sea ice. Moreover, in the Arctic the amount of in-ice chl a was always lower than phytoplankton chl a, whereas the opposite is true in the Southern Ocean. This highlights the importance of sea icederived carbon especially for Antarctic ecosystems, both in summer and in winter.
Sea-ice represents a permanent or seasonal habitat for many species in the polar oceans. The SUIT does not only uniquely cover relatively large areas, but also enables a simultaneous sampling of environmental data and under-ice fauna (collected with the two nets attached to the SUIT frame). Collected data have, therefore, enabled the study of the relationship between the sea-ice environment and the distribution of zooplankton and nekton species at the ice-water interface at a scale of kilometers. The surface zooplankton communities of both the Arctic and Southern Ocean have been found to respond to the presence (or absence) of sea ice in all sampled seasons, and variations in under-ice zooplankton community structure have been related to sea-ice environmental properties (Flores et al., 2011(Flores et al., , 2012(Flores et al., , 2014David et al., 2015David et al., , 2016Schaafsma et al., 2017;Ehrlich et al., 2020). In the Arctic Ocean, for example, abrupt changes were observed from a dominance of ice-associated amphipods (Apherusa glacialis) at ice-covered stations to a dominance of pelagic amphipods (Themisto libellula) at nearby ice-free stations (David et al., 2015). In addition, high numbers of the copepods Calanus hyperboreus and Calanus glacialis were found in the icewater interface layer, in spite of low under-ice water chl a within the sampling area (David et al., 2015), which indicated that they were grazing on ice algae (Runge and Ingram, 1991;Kohlbach et al., 2016). Higher abundances of the polar cod (Boreogadus saida), which is regarded as a key species in the Arctic Ocean, were related to relatively low under-ice surface salinity, thick seaice and high sea-ice coverage . In the Southern Ocean, differences in sea-ice drift pathways were closely related to the condition of krill larvae and the proportion of ice algae in their diet, at the beginning of the productive season (PS81) in the Weddell Sea (Kohlbach et al., 2017b;Schaafsma et al., 2017). The community structure of under-ice fauna was associated to changes in under-ice water parameters .
The presented data set and the environmental information from this study allow for further research on the relationship between large scale distributional patterns of under-ice fauna and sea-ice structures, such as the abundance and size of pressure ridges and hummocks, which have been suggested as structures where organisms, including both algae and zooplankton, accumulate (Gradinger et al., 2010;Lange et al., 2017a,b). Studies of Gradinger (1999) and Nozais et al. (2001) analyzed the potential role of sympagic fauna in controlling algal production due to grazing, but both drew contradictory conclusions. Whereas, Nozais et al. (2001) strongly suggest that the grazing impact of sea-ice meiofauna on ice algae was negligible, Gradinger (1999) found significant positive correlations between the integrated ice-algae and sympagic fauna biomass, and potential ingestion rates. Also other studies, often using data collected by divers, found evidence that abundance and distribution of certain species were related to the structure of the sea ice (Gruzov et al., 1967;Carey, 1985;Grainger and Hsiao, 1990;Garrison, 2015;Meyer B. et al., 2017), examples being the amphipods Gammarus wilkitzkii (Beuchel and Lønne, 2002) and Apherusa glaciali (Poltermann, 1998;Gradinger et al., 2010) in the Arctic Ocean. Gaining knowledge on the relationship between structure and faunal biomass on a larger scale, potentially will allow for prediction of an oceanwide distribution based on sea-ice parameters.

CONCLUSIONS
In this study, we presented and analyzed a large data set collected during five expeditions covering different years, seasons, and regions in both the Arctic Ocean and the Southern Ocean. Data were collected with a Surface and Under-Ice Trawl (SUIT), which allows the characterization of the sympagic environment on the kilometer scale. The data proved to be representative in characterizing the physical environment, and highlighted regional and seasonal differences in both hemispheres. Analysis of spatial variability showed a significant difference in the spatial scale of the variability of sea-ice properties between the Arctic and Antarctic, and showed that the total ice thickness was a driver for variability in sea-ice transmittance. The introduction of a new parameter, the insolation (i.e., daily integrated underice PAR) showed that more light for photosynthesis is available for biological production in the Antarctic summer compared to the Arctic summer, despite the large snow cover in the Antarctic and the higher transmittance in the Arctic. For the first time, we could provide in-ice chl a estimates over large spatial scales derived with the same methodology, for both the Arctic and the Antarctic. Our results showed that ice algae biomass was generally at comparable levels in both hemispheres during summertime, it was low in Arctic spring and particularly high in the Antarctic winter. In a future study, the present data set will be used to compare results and to improve parameters of sea-ice biogeochemical models.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available at https://doi. pangaea.de/10.1594/PANGAEA.902056.

AUTHOR CONTRIBUTIONS
GC wrote the paper, processed the data to obtain sea-ice draft and in-ice chl a, and prepared most of the figures. SA processed the RAMSES data to obtain integrated irradiance and radiance values and also to extract the wavelengths. RR provided the satellite ice concentration data. PM provided modeled downwelling irradiance data. IP supervised the Arctic sea ice biology sampling and processed all HPLC samples. FS, HF, and JE contributed to develop the structure of the manuscript and FS provided information on the snow depth per each SUIT haul. TK, SH, SS, and BL collected the EM-bird data. BL processed and plotted the EM-data, provided in-ice chl a estimates for PS80, and carried out the spatial autocorrelation analysis. GC, HF, JE, BL, FS, and CD contributed to SUIT data collection. JE, FS, CD, and HF contributed to data interpretation in connection with biological data. HF supervised the construction of the SUIT and designed the sampling for all the campaigns and conducted statistical analyses. All co-authors contributed to the discussion of results and to the preparation of the manuscript.