# Advances in altimetric snow depth estimates using bi-frequency SARAL and CryoSat-2 Ka-Ku measurements

Florent Garnier, Sara Fleury, Gilles Garric, Jérôme Bouffard, Michel

Tsamados, Antoine Laforge, Marion Bocquet, Renée Mie Fredensborg Hansen,
Frédérique Remy

<figure><img src="image_1.png" /></figure>

## To cite this version:

Florent Garnier, Sara Fleury, Gilles Garric, Jérôme Bouffard, Michel Tsamados, et al.. Advances in altimetric snow depth estimates using bi-frequency SARAL and CryoSat-2 Ka-Ku measurements. The Cryosphere, 2021, 15, pp.5483-5512. 10.5194/tc-15-5483-2021. insu-03671315

HAL Id: insu-03671315

https://insu.hal.science/insu-03671315

Submitted on 18 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

# Advances in altimetric snow depth estimates using bi-frequency SARAL and CryoSat-2 Ka–Ku measurements

Florent Garnier<fnref n="1" />, Sara Fleury<fnref n="1" />, Gilles Garric<fnref n="2" />, Jérôme Bouffard<fnref n="3" />, Michel Tsamados<fnref n="4" />, Antoine Laforge<fnref n="5" />,
Marion Bocquet<fnref n="1" />, Renée Mie Fredensborg Hansen<fnref n="3" />, and Frédérique Remy<fnref n="1" />

**Correspondence:** Florent Garnier (florent.garnier@legosobs-mip.fr)

Received: 3 March 2021 – Discussion started: 13 April 2021

Revised: 24 September 2021 – Accepted: 25 October 2021 – Published: 10 December 2021

**Abstract.** Although snow depth on sea ice is a key parameter for sea ice thickness (SIT) retrieval, there currently does not exist reliable estimations. In the Arctic, nearly all SIT products use a snow depth climatology (the modified Warren-99 climatology, W99m) constructed from in situ data obtained prior to the first significant impacts of climate change. In the Antarctic, the lack of information on snow depth remains a major obstacle in the development of reliable SIT products.

In this study, we present the latest version of the altimetric snow depth (ASD) product computed over both hemispheres from the difference of the radar penetration into the snow pack between the Ka-band frequency SARAL/Altika and the Ku-band frequency CryoSat-2. The ASD solution is compared against a wide range of snow depth products including model data (Pan-Arctic Ice-Ocean Modelling and Assimilation System (PIOMAS) or its equivalent in the Antarctic the Global Ice-Ocean Modeling and Assimilation System (GIOMAS), the MERCATOR model, and NASA’s Eulerian Snow On Sea Ice Model (NESOSIM, only in the Arctic)), the Advanced Microwave Scanning Radiometer-2 (AMSR2) passive radiometer data, and the Dual-altimeter Snow Thickness (DuST) Ka–Ku product (only in the Arctic). The ASD product is further validated in the Arctic against the ice mass balance (IMB) buoys, the CryoSat Validation Experiment (CryoVEx) and Operation Ice Bridge’s (OIB) airborne measurements. These comparisons demonstrate that ASD is a rel-

evant snow depth solution, with spatiotemporal patterns consistent with those of the alternative Ka–Ku DuST product but with a mean bias of about 6.5 cm. We also demonstrate that ASD is consistent with the validation data: comparisons with OIB’s airborne snow radar in the Arctic during the period of 2014–2018 show a correlation of 0.66 and a RMSE of about 6 cm. Furthermore, a first-guess monthly climatology has been constructed in the Arctic from the ASD product, which shows a good agreement with OIB during 2009–2012. This climatology is shown to provide a better solution than the W99m climatology when compared with validation data. Finally, we have characterised the SIT uncertainty due to the snow depth from an ensemble of SIT solutions computed for the Arctic by using the different snow depth products previously used in the comparison with the ASD product. During the period of 2013–2019, we found a spatially averaged SIT mean standard deviation of 20 cm. Deviations between SIT estimations due to snow depths can reach up to 77 cm. Using the ASD data instead of W99m to estimate SIT over this time period leads to a reduction in the average SIT of about 30 cm.

<footnote n="1">¹Laboratoire d'Etudes en Géophysique et Océanographie Spatiales (LEGOS), CNRS/UMR5566,
Université Paul Sabbatier, 31400 Toulouse, France</footnote>

<footnote n="2">²Mercator Ocean, 31520 Ramonville Saint Agne, France</footnote>

<footnote n="3">³Earth Observation Directorate, ESA (European Space Agency), Via Galileo Galilei, 2-00044 Frascati, Italy</footnote>

<footnote n="4">⁴Centre for Polar Observation and Modelling, Department of Earth Sciences,
University College London, London, WC1E 6BT, UK</footnote>

<footnote n="5">⁵Serco c/o ESA, Earth Observation Directorate, Via Galileo Galilei, 2-00044 Frascati, Italy</footnote>
# 1 Introduction

Since the launch of CryoSat-2 (CS-2) in 2010 (Wingham et al., 2006; Parrinello et al., 2018), sea ice thickness (SIT) observations are routinely derived from altimetric measurements. The principle is to measure the fraction of the sea ice above the sea level, called the sea ice freeboard, from differences between the heights in leads (cracks in the ice referring to the local sea level) and the heights of the ice floes (Laxon et al., 2003). By integrating such sea ice freeboard estimations in the hydrostatic equilibrium equation, several SIT products have been computed (e.g. Laxon et al., 2013; Kwok and Cunningham, 2015; Guerreiro et al., 2017; Paul et al., 2018; Landy et al., 2019; Laforge et al., 2020).

Among the parameters involved in the SIT calculation, snow depth (sd) over sea ice is one of the most significant contributors adding to the overall SIT uncertainty (e.g. Giles et al., 2007; Zygmuntowska et al., 2014; Guerreiro et al., 2016). For example, it is necessary to account for the snow loading (Laxon et al., 2013) and for the decrease in altimetric radar speed as it penetrates into the snow pack (Kwok and Cunningham, 2015; Mallett et al., 2020). In fact, abilities in snow properties affect radar signals and then the snow depth measurements. More generally, snow cover has strong impacts on the sea ice (e.g. Massom et al., 2001; Powell et al., 2005; Bin et al., 2008; Sturm and Massom, 2009; Ricker et al., 2014) that affects the entire climate system (e.g. Ingram et al., 1989; Ledley, 1991; Eicken et al., 1995; Sin Garayer et al., 2006). Because of its high albedo and a low thermal conductivity, the snow regulates the transfer of solar heat energy penetration across the ice-ocean interface (e.g. Grenfell and Maykut, 1977; Sturm et al., 1997). It acts as an insulator, slowing down sea ice melt in summer and slowing down sea ice growth in winter (e.g. Perovich et al., 2003; Sturm and Massom, 2016). Such processes of sea ice formation and melting govern sea ice physical and chemical properties that impact the biological processes in sea ice (Van Leeuwe et al., 2018). The vertical distribution of light under sea ice that controls biological processes and biogeochemical fluxes is also strongly linked with the snow depth (e.g. Perovich, 2007; Arndt and Nicolaus, 2014; Arndt et al., 2017). In addition snow accumulated over sea ice from precipitation represents a tank of freshwater likely to be carried into the ocean. Recent increases in seasonal ice have promoted the amount of snow water discharge in the ocean, which impacts the freshwater budget (Andersen et al., 2019; Overland et al., 2019). Snow cover also modifies surface roughness that impacts the air-ice drag coefficient and transfer coefficients of latent and sensible heat fluxes (Andreas et al., 2005).

Despite evidence of its major importance, snow cover over sea ice is still not sufficiently understood. Until now, nearly all SIT products in the Arctic are computed using the snow depth of the modified Warren-99 climatology (W99m; Warren et al., 1999). W99m is mainly constructed using snow

depth observations obtained from the Soviet North Pole drifting stations based in the central Arctic, collected during the last century (1950–1980s). Considering the ongoing rapid modifications of the Arctic due to climate change, several studies have shown these data to be outdated, even with the modification of snow depth reduction (by 50%) applied over the seasonal ice zone (Kurtz and Farrell, 2011; Kern et al., 2015; Kwok et al., 2011). The climatology of Forsström et al. (2011) provides snow depths in areas outside the central Arctic where the W99m climatology is not working properly. More recently, Shalina and Sandven (2018) have created an improved snow depth climatology taking into account snow depth observations obtained during campaigns over seasonal ice. However, it is still based on data collected mainly in the 1960s, 1970s and 1980s. Apart from the data used to construct these climatologies, several campaigns have provided snow depth measurements in the Arctic. Between 1997 and 1998, the SHEBA (Surface HEat Budget of the Arctic Ocean) project has highlighted the complex temporal and spatial snow depth variations (Perovich et al., 1999; Sturm et al., 2002). The ice mass balance (IMB) buoys, originally deployed by the US Army CRREL-Dartmouth Mass Balance Buoy Program during the SHEBA project have measured snow depths since the 2000s. It allows for monitoring changes in the sea ice volume in key areas of the Arctic (Richter-Menge et al., 2006). Since 2003, the CryoSat Validation Experiment (CryoVEx) campaigns (e.g. Haas et al., 2006; Helm et al., 2006; Skourup et al., 2013; Hvidegaard et al., 2006) have provided data with the main goal of investigating radar penetrations into ice and snow cover. Their measurements include bi-frequency altimetry snow depth estimations as of 2017 (see Sect. 3.3). In 2003, the AMSR-E03 campaign (Sturm et al., 2006) was carried out in the Beaufort Sea region to validate the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E, from 2002 to 2011; Kelly, 2009) passive microwave snow depth data. The same year, an equivalent mission (ARISE) took place in East Antarctica (Massom et al., 2006). Since 2009, and until their final campaign in 2020, the Operation Ice Bridge (OIB; see Sect. 3.3) campaigns have provided airborne snow depth measurements each year between March and April in the Arctic (Kurtz et al., 2013; Koenig et al., 2010) using a frequency modulated continuous wave (FMCW) snow radar (Kurtz and Farrell, 2011). Note also that Webster et al. (2014) assess spring snow depth distribution on the Arctic sea ice from 2009–2013 airborne OIB radar observations. In parallel, in situ snow thickness data were collected in 2009 at the Danish GreenArc 2009 ice camp in the north of Greenland (Farrell et al., 2011). Between February and March 2015, snow depth was measured in the Nansen Basin during the N-ICE expedition (Granskog et al., 2016), demonstrating among others the impact of heavy precipitation on thin ice growth (Merkouriadi et al., 2017). From September 2019 to October 2020, the international MOSAIC expedition monitored the central Arctic (Shupe et al., 2020) provided
ing, among others, data on snow precipitation, snow water equivalent and Magnaprobe snow depth measurements (e.g. Munoz-Martin et al., 2020). Hence, several campaigns have been conducted during the last two decades in the Arctic, all of which have provided crucial information and proved invaluable as reference data.

Compared to the Arctic, snow depth measurements, and more generally sea ice in situ data, are very limited in the Antarctic. Actually, several CryoVEx and OIB missions have been conducted in the Antarctic, but the data are still only available in raw L1b level and are difficult to interpret. To the best of our knowledge, only the meereis data portal (https://seaiceportal.de/en/, last access: March 2021) provides snow buoy data of snow accumulation (Grosfeld et al., 2016) (for both hemispheres), but the comparison with snow depth is not direct. It is, for instance, limited by the flooding of ice floes due to the heavy snow loading occurring in the Antarctic. Currently, the main expectation providing snow depth data in the Southern Hemisphere is the ASPeCt programme (Worby et al., 2008a), established in 1996 to model the role of the Antarctic sea ice in the coupled atmosphere-ice-ocean system. Recent ASPeCt standardised data (Kern, 2020) cover the period 2002–2019 with an amount of information that is too sparse for a reliable assessment. Future validations in the Antarctic will also benefit from the recent AWI IceBird campaigns carrying a snow radar since 2019 (https://www.awi.de/en/science/climate-sciences/sea-ice-physics/projects/ice-bird.html, last access: March 2021).

At a basin scale, several sea ice models produce snow accumulation estimations from atmospheric reanalyses. These are converted into snow depth over sea ice by considering sea ice drifting and thermodynamic transfers of the coupled sea ice and snow system (Blanchard-Wrigglesworth et al., 2015; Lecomte et al., 2011). Likely due to biases in precipitation inputs and omission of snow processes, these global sea ice models hardly reproduce consistent on-ice snow depths (e.g. Holland et al., 1993; Serreze et al., 2000; Déry and Tremblay, 2004; Leonard and Maksym, 2011; Blazey et al., 2013). Notz (2012) states that an improved representation of the snow on sea ice in models is the most urgent task, and recently, Kaminski et al. (2018) argued that assimilating snow depth products could considerably increase sea ice model performances.

Therefore, a solution allowing us to provide large-scale information on the spatial and temporal variation in snow depth and to improve the performance of seasonal sea ice forecasts is to use satellite data. For now, only a few datasets exist. The first snow depth estimates come from passive radiometric data. Among the various sensors that have been deployed (e.g. the SMMR, Chang et al., 1987, or the SSM/I, Grody, 1991), the most relevant estimations of snow depth over sea ice have been calculated from spectral vertical gradients of brightness temperatures (Comiso et al., 2003; Worby et al., 2008b) of AMSR-E and its successor AMSR2 (Lee

et al., 2015). In spite of the relevant long-term availability of data, the commonly used snow depth retrieval algorithm is not adapted for multi-year ice (MYI) zones (Brucker and Markus, 2013). Consequently, AMSR's official data in the Arctic are only calculated in first year ice (FYI) regions. As far as we know, the AMSR2 product version of the University of Bremen (AMSR2B; Rostosky et al., 2018) and very recently the SPICES project snow depth retrieval method (Mäkynen et al., 2020) are the only attempts to compute snow depths in the Arctic MYI zones. The dataset developed from the work of Rostosky et al. (2018) is available but has the inconvenience of being re-calibrated from OIB data, which limits their availability to the months of March and April. With the same disadvantage of depending on OIB, Maaß et al. (2013) and Zhou et al. (2018) have proposed algorithms to retrieve snow depth from SMOS (Soil Moisture Ocean Satellite), with mitigated results. Recently, Braakmann-Folgmann and Donlon (2019) proposed to compute snow depth by combining AMSR2 and SMOS data with a neural network approach.

Studies have also estimated snow depth on sea ice in the Antarctic (e.g. Kacimi and Kwok, 2020; Maksym and Markus, 2008). AMSR-E and AMSR2 also provide snow depths in the Antarctic. Compare to the Arctic it covers the whole domain since all sea ice is considered as seasonal ice. The recent work of Kacimi and Kwok (2020) estimated the Antarctic snow depth from the difference between the total freeboard obtained from ICESat-2 lidar measurements and sea ice freeboards of CS-2. They show a very low variability in CS-2 freeboards compared to ICESat-2 and that ICESat-2 freeboards explained >90% of the variance in snow depth. They also highlight that the validation sea ice parameters in the Antarctic remains a challenge since no seasonally or regionally diverse datasets from field records can be used to assess the large-scale satellite retrievals. Based on this, they urge that more sustained and extensive field measurements must be conducted. Note that such lidar minus Ku-band snow depths can also be computed in the Arctic. Finally, we note that a recent study of Fons et al. (2021), an update from the former approach experimented by Fons and Kurtz (2019), directly retrieves snow freeboard from CS-2 using a two-layer scattering model and shows promise in estimations of SIT in the Antarctic.

An alternative to passive microwave radiometry is to derive snow depth estimations from bi-frequency altimetric measurements. Recently, Armitage and Ridout (2015) have demonstrated this possibility by considering the difference of penetration between the Ku-band frequency CS-2 (13.5 GHz), assuming it is reflected near the snow-ice interface, and the Ka-band frequency SARAL/Altika (35.7 GHz), assuming it is reflected near the top of the snow pack, i.e the air-snow interface. Thereafter, a preliminary altimetric snow depth (ASD) version covering the 2013–2016 winter period has been developed at the LEGOS during the ESA CryoSeaNICE project (Guerreiro et al., 2016). Meanwhile,
Lawrence et al. (2018) also developed a bi-frequency altimeter snow depth product (DuST) available for the Arctic. Similar to AMSR2B, the main difference between ASD and DuST is that DuST relies on a re-calibration of the Ka-Ku freeboards using OIB data, in this case to account for the difference of operating mode between CS-2 and SARAL (synthetic aperture radar (SAR) versus low-resolution mode (LRM) modes). These two products are the only existing publicly available snow depth products from altimetry.

A recent study of Zhou et al. (2021) presented an inter-comparison of available snow depth products from re-analyses, passive radiometry and altimetry (DuST). Similar, this paper reviews the state of the art by comparing current main snow depth estimations. Yet, the main objective is to present and assess the upgraded version of the ASD product (see Sect. 2), covering the 2013–2019 period in both hemispheres. The article fits in with the upcoming HPCM CRISTAL mission (Kern et al., 2020) and aims to demonstrate the potential of such snow depth data to further specific studies on, for instance, improved sea ice volume estimations, freshwater budgets, snow properties or data assimilation. Except from an analysis of the impact of the snow depth uncertainty on SIT retrieval, it does not explicitly address these questions.

The paper is organised as follows:

- First, we detail the methodology to process the ASD product and present all the datasets used in this study.

- Section 4 compares ASD with the other main existing snow depth satellite and model data in both hemispheres.

- The datasets are then assessed against OIB, CryoVEx and IMB validation data in the Arctic in Sect. 5.

- To circumvent the temporal limitation induced by SARAL, we propose in Sect. 6 a preliminary snow depth climatology based on ASD.

- Section 7 aims to quantify the SIT level of uncertainty due to the snow depth from an ensemble of SIT estimations calculated from the satellite and model snow depth datasets presented in the previous sections.

- We finally discuss and conclude, emphasising current needs for snow depth data in sea ice studies.

## 2 Data processing of ASD

The Altimetric Snow Depth (ASD) product presented here is an upgraded version of the data presented in Guerreiro et al. (2016). It has been developed at the LEGOS laboratory in the framework of the ESA CSAO+ (http://cryosat.mssl.ucl.ac.uk/csao/index.html, last access: March 2021) and Polar+ Snow on Ice projects. It has been freely available

from the AVISO+/ODATIS national data centres since mid-2020 (http://ctoh.legosobs.mip.fr/data/sea-ice-products, last access: March 2021). The ASD product includes data from March 2013 to October 2019 provided on a monthly basis for the 6 winter months (from November to April in the Arctic and from May to October in the Antarctic). It is projected onto a 500 × 500 EASE2 grid with a 12.5 km pixel size resolution. Snow depth calculation is based on the difference of penetration between Altika, the Ka-band range altimeter of SARAL (which is assumed to be reflected at the air-snow interface), and SIRAL, the Ku-band range altimeter of CS-2 (which is assumed to be reflected at the snow-ice interface). The main assumption is that the difference between these two surface elevations is only due to the penetration of the Ku radar in the snow pack and that the Ku radar penetrates fully to the snow-ice interface. Note that the validity of this hypothesis is strongly linked with the space and time variability in snow properties such as snow density, grain sizes, the liquid water content or the surface roughness. We use the CS-2 Ku-band waveforms in pseudo LRM (PLRM) extracted from geophysical ocean product (GOP) (https://earth.esa.int/documents/10174/125272/CryoSat-Baseline-C-Ocean-Product-Handbook, last access: March 2021) generated by the ESA CryoSat Ocean Processor (Bouffard et al., 2018a). The Ka-band SARAL/Altika waveforms are extracted from the Centre national d'études spatiales (CNES) SGDR-T (Sensor Geophysical Data Record) official products. Note that SARAL's orbit is limited to 81.5° in latitude, which reduces the coverage of snow depth data in the Arctic. The use of a degraded version of the CS-2 SAR waveforms (the PLRM mode) allows us to have a footprint in accordance with the LRM mode of SARAL/Altika waveforms, which avoids having different impact of the surface roughness. Since the latest Baseline-C PLRM GOP was only available from 2017 at the time we computed the ASD data, we have used the Baseline-B for the period 2013–2016. It does not impact ASD since we use only the L1b product levels which have identical waveforms on both baselines. However, the Baseline-B does not include the SARIN (interferometric synthetic aperture radar) data. The next version of the ASD product will be produced with only the Baseline-C PLRM GOP to include all SARIN mode zones.

The first step to compute ASD data consists in extracting heights, H, from both Ku-band and Ka-band satellite waveforms following Eq. (1):

$$H = \text{alt} - (\text{range} - (\text{tropo}_{\text{dry}} + \text{tropo}_{\text{wet}} + \text{iono} + \text{MSS} + \text{tide}_{\text{ocean}} + \text{tide}_{\text{bar}})), \quad (1)$$

where alt is the satellite altitude, range is the altimetric range retracked from the waveforms, tropo<sub>dry</sub> is the dry tropospheric correction, tropo<sub>wet</sub> is the wet tropospheric correction, iono is the ionospheric correction, MSS is the DTU15 mean sea surface correction (Andersen and Knudsen, 2015),
tide_ocean is the oceanic tide corrections, and tide\_bar is the barometer tide correction.

The surface classification between heights corresponding to ocean surfaces (leads) and heights of sea ice surfaces (ice floes) is performed using the waveform pulse-peakiness (PP) criteria calculated from the waveforms (Eq. 2):

$$PP = \frac{\max(\text{WF})}{\sum_{i=1}^{N_{\text{WF}}} \text{WF}_i}, \qquad (2)$$

where WF is the discrete waveform echoes, and N\_WF is the number of along-track measurements. Following Guerreiro et al. (2017), surfaces with a PP larger than 0.3 are considered as leads, and surfaces with a PP smaller than 0.1 are considered as sea ice floes. Observations with a PP between 0.1 and 0.3 are considered as mixed echoes and are discarded. This criteria is the same for SARAL and CS-2.

Altimetric ranges (range parameter in Eq. 1) of leads and floes are thereafter calculated using the TFMRA (Helm et al., 2014) with a 50% threshold. A 25 km median filter is then applied to the 20 Hz along-track surface height of sea ice floes and sea ice leads. For each observation, the sea level anomalies (SLAs) under floes are calculated as the median value of all leads in a 25 km radius around the considered along-track point. Freeboard heights are computed by fb = H\_floe - H\_lead, where H\_floe is the heights of the ice floes and H\_leads the heights of the leads. Finally, a 25 km radius median smoothing is applied to the retrieved freeboards. The reason is that we assume that the ASD product should not be able to provide relevant information at smaller scales. However, additional analyses and comparisons with validation data would be necessary to properly characterise this point. Note that the ability to consistently observe small scales would certainly be significantly improved in the future dual-frequency snow depth products from the CRISTAL mission. In the first ASD data version (Guerreiro et al., 2016), snow depths were only calculated at CS-2 and SARAL crossing track points. For this version, snow depths are calculated from the difference between monthly EASE2 gridded freeboards as sdr = fb\_Ka - fb\_Ku, where sdr is a radar snow depth. Furthermore, to compute snow depth (sd), we need to take into account the decrease in the Ku radar echo velocity as it penetrates into the snow pack (Ulaby et al., 1986; Kwok and Cunningham, 2015):

$$sd = sdr \times (1 + 0.51 \times \rho_s)^{(-1.5)}, \qquad (3)$$

where ρ_s is the snow density, set to 300 kg/m³. Note this formulation agrees with the conventional interpretation to correct for the slower wave propagation speed through snow, as recently shown in Mallett et al. (2020).

To estimate snow depth uncertainties, we assume that errors are unbiased, uncorrelated and follow the Gaussian propagation law. Snow depth uncertainty is then given by

$$\delta_{sd} = \sqrt{(\delta_{sd_r} \times C)^2 + (sd_r \times B \times \delta_{\rho_s})^2}, \qquad (4)$$

with

$$C = (1 + 0.51 \times \rho_s)^{-1.5}, \\ B = -1.5 \times 0.51 \times (1 + 0.51 \times \rho_s)^{-2.5}, \\ \delta_{sd_r} = \sqrt{\delta_{fb\_Ka}^2 + \delta_{fb\_Kur}^2} \quad \text{and} \quad \delta_{\rho_s} = 3.2.$$

The freeboard uncertainties δ\_fb = √(δ\_H\_leads² + δ\_H\_floes²) are calculated along-track from the statistical dispersion of heights in a 12.5 km radius. Since sea ice topography can significantly vary within this area, we assume that the statistical dispersion of floes is the same as that of the leads.

$$\delta_{H\_floes}^2 = \frac{\sigma_{H\_leads}}{N_{floes}} \quad \text{and} \quad \delta_{H\_leads}^2 = \frac{\sigma_{H\_leads}}{N_{leads}} \qquad (5)$$

Using this methodology, Fig. 1 presents an example of ASD snow depth estimations and its associated uncertainties in the two hemispheres. In the Antarctic, thicker snow is located in the Weddell Sea and, to a lesser extent, in the coastal zones of the Bellingshausen and Amundsen seas. It is also relevant to identify the very low snow values associated with the Ross Ice Shelf. In the Arctic, the snow distribution follows well the dynamics of sea ice, the most characteristic of which is the export of MYI in the Beaufort Gyre. This figure also shows that the ASD data are very different from the W99m climatology, in which the W99m climatology tends to exhibit thicker snow layers over sea ice. The ASD mean level of uncertainties (of about 4 cm on average with standard deviation of 1 cm) are smaller than the deviations between these two products (on average about 6 cm with standard deviation of 7 cm).

## 3 External datasets

This section presents the various snow depth datasets used in this study. For both hemispheres, the time period of model and satellite products ranges from 2014–2019 for the Arctic and 2013–2019 for the Antarctic (with limitations for some data products for both hemispheres), as explained in Fig. 2. Table A4 specifies the time period of the different satellite and model data.

### 3.1 Satellite data

#### 3.1.1 DuST

The DuST (Lawrence et al., 2018) data are provided in the Arctic on a 1.5° longitude × 0.5° latitude grid by the CPOM/UCL (http://www.cpom.ucl.ac.uk/DuST/, last access: March 2021). They are available for the 6 winter month (November to April) until 2018. Similar to ASD data, DuST relies on the difference of penetration between the Ka-band SARAL/Altika and the Ku-band CS-2 radar altimeters. However, unlike the ASD product (which uses the CS-2 PLRM
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Maps of the altimetric snow depth (ASD) annual mean snow depth (a) and its uncertainties (d) in the Antarctic and the Arctic (b and e) for 2015. The third column presents the annual mean snow depth map of the modified Warren-99 climatology (c) and its difference with the ASD product (ASD-W99m) in the Arctic (f) in 2015.</figcaption></figure>

data; see Sect. 2), the DuST product uses SAR mode freeboard estimates of CS-2. The CS-2 and SARAL (LRM data) data are then calibrated against data from National Aeronautics and Space Administration's (NASA's) OIB airborne measurements to align the freeboard observations from SARAL/Altika with the snow-air interface and the freeboard observations from CS-2 with the snow-ice interface. Although this method provides a generic approach adaptable to different footprints (e.g. ICESat-1 and 2), it has the disadvantage of being dependent on the OIB data which are mainly located near the Canadian Archipelago and the Beaufort Sea and are only obtained during March and April.

### 3.1.2 AMSR2

The Advanced Microwave Scanning Radiometer-2 (AMSR2) is a passive radiometer on board the JAXA/GCOM-W (Japan Aerospace Exploration Agency's Global Change Observation Mission - Water) satellite launched in 2012. Snow depth on sea ice retrieval is based on the gradient ratio of vertically polarised brightness temperature at 19 and 37 GHz (Markus and Cavalieri, 1998), and two snow depth products are available. The first product (AMSR2-NSIDC), described in Meier et al. (2018), is accessible on the NSIDC (https://nsidc.org/data/AU_S112/versions/1, last access: March 2021). It provides daily L3 snow depth data in

12.5 km × 12.5 km stereopolar grids constructed as 5d running averages. Data are only calculated over FYI (with MYI concentration of less than 20 %). This is due to the fact that MYI has a spectral signature similar to snow cover over FYI. Because of this limitation, we mainly focus on the Bremen AMSR2B product (Rostosky et al., 2018) developed at the Institute of Environmental Physics of the University of Bremen. This product is available on a daily basis from November 2012 to April 2018 on a polar stereographic grid with a 25 km × 25 km resolution. Unlike in AMSR2-NSIDC, March and April snow depth data are also calculated over MYI for AMSR2B. Other months are only available over FYI. Note that the algorithm used to calculate snow depth over MYI is also calibrated with observations from OIB.

## 3.2 Model data

### 3.2.1 PIOMAS

The pan-Arctic Ice-Ocean Modelling and Assimilation System (PIOMAS; Global Ice-Ocean Modeling and Assimilation System, GIOMAS, in the Antarctic) is a pan-Arctic coupled ocean and sea ice model developed for climate applications (Zhang and Rothrock, 2003). Snow depth data are provided on a daily basis in a 360 × 120 generalised curvilinear coordinate system from the PIOMAS model version 2.1 (http://psc.apl.uw.edu/research/projects/
arctic-sea-ice-volume-anomaly/data/model_grid, last access: March 2021). The atmospheric surface forcing fields are issued from NCEP-NCAR reanalysis. Ocean and Sea Ice Satellite Application Facility (OSI SAF) sea ice concentration from EUMETSAT is assimilated in near real time. A complete description can be found in Schweiger et al. (2011).

### 3.2.2 MERCATOR

In the framework of the Copernicus Marine Environment Monitoring Service (CMEMS), Mercator Ocean implements a real-time global analysis and forecasting system at 1/12° resolution. This system is based on the NEMO ocean platform (Madec et al., 2015) and is driven by atmospheric conditions from the ECMWF Integrated Forecast System (IFS) analysis and forecasting system. The sea ice model component is the Louvain-la-Neuve sea ice model (LIM-2) (Fichefet and Maqueda, 1997; Vancoppenolle et al., 2012). The unique sea ice quantity assimilated within the system is the near-real-time sea ice concentration from OSI SAF provided by CMEMS. The accumulation of the snow depth onto sea ice originates from ECMWF snowfall forcing fields. A complete description can be found in Lellouche et al. (2018). Snow depth data used in this study have been provided on a monthly basis on the ORCA tripolar grid. This LIM-2 model configuration will be referred to as the MERCATOR model.

### 3.2.3 NESOSIM

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a 3D snow budget model configured to produce snow depth and density over sea ice in the Arctic Ocean. Data used in this paper are snow depth monthly mean maps provided on a 100 km × 100 km stereographic polar grid issued from the NESOSIM 1.0 configuration (Petty et al., 2018). It is freely available at https://earth.gsfc.nasa.gov/cryo/data/nasa-eulerian-snow-sea-ice-model-nesosim (last access: March 2021). Snowfall data forcing fields are the ECMWF ERA-Interim global reanalysis and a median of three reanalyses: the 55-year reanalysis of the JRA, NASA's MERRA and ASR, version 1.

## 3.3 Validation data

The validation data presented in this section are only compared with ASD in the Arctic. As mentioned in the introduction CryoVEx and OIB data in the Antarctic are still only L1b. IMB buoys are only in the Arctic. The geographical locations of the validation data are shown in Sect. 4.3.

### 3.3.1 CryoSat Validation EXperiment (CryoVEx)

These airborne data are a joint effort of the DTU National Space Institute and ESA (in cooperation with AWI). The main goal of the project is to quantify and vali-

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

**Figure 2.** Times series of annual mean snow depth of the different products in both hemispheres. Annual means are calculated as the average of the spatial means of all monthly grid maps. Only the 6 winter months are considered in the calculation (November–April in the Arctic and May–October in the Antarctic). Note that the AMSR2 data are not represented in the Arctic since the AMSR2-NSIDC product is only available over FYI, and the AMSR2B product is only available with its full spatial coverage (FYI and MYI) in March and April.

date CS-2's ability to measure SIT. Since 2003, the CryoVEx missions have provided both laser (ALS, airborne laser scanner) and Ku-band radar (ASIRAS, ESA airborne Ku-band interferometric radar) altimetric measurements. In addition, in 2017, the data include for the first time along-track total freeboards from the Ka-band KAREN altimeter (Haas et al., 2017), allowing us to compute Ka-Ku snow depth estimations. These data have been devel-
oped by AWI in the context of the ESA CryoSeaNICE project. For further information, please see the technical report at https://earth.esa.int/eogateway/documents/20142/1526226/CryoVEx2017-final-report.pdf (last access: March 2021).

### 3.3.2 Operation Ice Bridge (OIB)

OIB is one of the largest airborne missions in polar regions aiming to determine sea ice properties. Among others, it carries a snow radar measuring the snow-air and snow-ice interfaces of the signal scattered by the area illuminated beneath the aircraft. All information concerning the different campaigns and instruments can be found in the literature (e.g. Newman et al., 2014; Kurtz et al., 2013; Armitage and Ridout, 2015). OIB snow depth data presented in this paper are the NSIDC OIB Quicklook version (Kurtz et al., 2012; King et al., 2015) validated from in situ data of the BROMEX campaign (Nghiem et al., 2013). These data are available at https://daacdata.apps.nsidc.org/pub/DATASETS/ICEBRIDGE/Evaluation_Products/ (last access: March 2021). For comparison with the ASD dataset, we use all spring campaigns during the period of 2014–2018. For comparison with the ASD climatology, we use all OIB spring campaigns during 2009–2012.

### 3.3.3 Ice mass balance (IMB) buoys

IMB data (Perovich et al., 2021) are issued from the US Army CRREL-Dartmouth Mass Balance Buoy Program. They provide point measurements of the ice mass balance characterised by a combination of snow and ice conditions (Richter-Menge et al., 2006; Perovich and Richter-Menge, 2015). Snow depth data are collected at 4 h intervals from acoustic sounders. In this study we analyse the data of three winters (2013–2014, 2014–2015 and 2015–2016). Note that we do not use the data of 2017 and 2018 because they are still provisional and subject to revision.

# 4 Comparison between snow depth data

## 4.1 Methodology

Model and satellite snow depth estimations are projected onto a 12.5 km pixel size EASE2 grid (similar to ASD) using a linear 2D multivariate interpolation. For each product we compute monthly mean maps for the 6 winter months (November–April in the Arctic and May–October in the Antarctic) for both hemispheres from March 2013 to October 2019. These maps are provided for 2015 in the Supplement (Figs. S1 to S13). Note that to improve the consistency of the analyses we always use the larger common time period for all the compared products. Thus, the covered time periods may vary depending on the availability of the snow depth

products. Three statistical diagnoses are used to characterise the differences between snow depth datasets:

1. The climatic annual mean is simply the average of all monthly snow depth maps from all years.

2. The mean annual variability (MAV) is the average, for the years y, of the annual snow depth standard deviation (SDy(sd)) maps computed from the six monthly mean maps. The MAV is calculated at each grid point following Eq. (6):

$$MAV = \frac{1}{N_y} \sum_y (SD_y(sd)) \quad \text{with} \\ SD_y(sd) = \sqrt{\frac{1}{N_m} \sum_m (sd_{y, m} - \overline{sd_y})^2}, \quad (6)$$

where sd_{y,m} is the monthly snow depth map of the month m (m = 11, 12, 01, 02, 03, 04 in the Arctic and m = 05, 06, 07, 08, 09, 10 in the Antarctic) for the year y, and $\overline{sd_y}$ is the snow depth annual mean map of the year y. N_m and N_y are the number of months (six in our case) and the number of year in the considered time period, respectively.

3. The mean interannual variability (MIV), which is the average of the interannual standard deviation SD_m(sd) maps for the 6 winter months m over the 2013–2019 time period. The MIV is calculated at each grid point following Eq. (7):

$$MIV = \frac{1}{N_m} \sum_y (SD_m(sd)) \quad \text{with} \\ SD_m(sd) = \sqrt{\frac{1}{N_y} \sum_y (sd_{y, m} - \overline{sd_m})^2}, \quad (7)$$

where $\overline{sd_m}$ is average snow depth map of the month m.

For now, the lack of available relevant snow depth data in the Antarctic limits these comparisons to the Northern Hemisphere. Comparisons with all in situ and validation data are performed with a comparable methodology: snow depth model and satellite gridded maps are projected along the aircraft trajectories using a bilinear interpolation. Airborne and in situ are generally daily data, and the comparisons are performed with the mean maps of the month to which that day belongs. In order to achieve similar spatial scales in the comparisons, we have applied a 25 km window rolling mean to smooth the external data. Due to the non-Gaussianities within 25 km sections, we have noticed that the smoothing has a mean effect which slightly tends to reduce the mean value of data. It is important to consider that this allows us to
correct potential bias that would be induced by a difference in resolution. Note that the space and time consistency in resolution between all model and satellite data is also ensured by projection onto similar grids. Since snow depth products and validation data are available over different time periods (see table A4), the comparisons are also performed for different time periods. Our approach is to take, for each comparison, the time period common to all the compared data. The aim is to provide reliable statistics. For each comparison, a label specifies the considered time period.

## 4.2 Comparison between satellite and model estimations

### 4.2.1 Results in the Arctic

Here, we compare the ASD, DuST and AMSR2B snow depth data in the Arctic by considering the months of March and April over the 2014–2018 period (Fig. 3). These climatic “seasonal” means and interannual variability presented in Fig. 3 are the mean and the standard deviation of all monthly maps of March and April over this time period. Statistical diagnoses are summarised in Table A1.

On average, the ASD data are about 3 cm lower than the AMSR2B data and 6 cm lower than the DuST data. ASD data present the stronger snow depth gradients (associated with a higher spatial standard deviation) with clear distinctions between low snow depths over thin sea ice regions (e.g. Laptev Sea region) and higher snow depths over regions of thicker sea ice (e.g. Queen Elizabeth Islands region). Except in the MYI zones, AMSR2B data are quite smooth and characterised by a weak variability compared to that of the two altimetric products. In MYI zones (towards Queen Elizabeth islands) AMSR2B has the highest snow depths. We assume that it might be related to the snow depth retrieval algorithm, inadequate over thick ice. Indeed, the statistical results which are presented in the Appendix (Table A1) show that the two Ka-Ku products (ASD and DuST) are in good agreement in terms of spatial variability and the magnitude of their annual cycle considering the entire set of data over the 6 winter months. Note that DuST presents the highest climatic “seasonal” mean solution, even higher than W99m (see also Fig. 2). In addition, DuST and ASD spatial distributions are comparable, showcasing that the difference between the two Ka-Ku altimetric products is likely a bias resulting from the re-calibration of the DuST data with OIB (this point will also be discussed in the next sections). This underlines the consistency of this methodology used to compute bi-frequency snow depth estimations in the absence of external data (ASD) since DuST calibration is limited by the spatial and temporal availability of OIB.

Figure 4 compares the ASD snow depth estimations with the outputs of three sea ice models in the Arctic (NESOSIM, MERCATOR and PIOMAS). The statistical diagnoses presented in Sect. 4.1 are used. They are summarised in Table A1.

Amongst the three models, the MERCATOR model provides lower snow depth estimations throughout with a climatic mean of 7.4 cm and lower mean interannual and annual variability of about 3 cm. On the contrary, climatic means of PIOMAS and NESOSIM are similar both in terms of spatial distribution and pan-Arctic mean (with values of respectively 15.5 and 16.2 cm). Considering annual and interannual variability, PIOMAS presents high values that are more evenly spread. In contrast, NESOSIM is characterised by lower values everywhere except in the Greenland and the Kara seas.

In spite of some comparable large-scale spatial patterns, these results highlight the large deviations between the ASD data and the model solutions. One striking feature is the transport of MYI along the Canadian coast by the Beaufort Gyre (see also Fig. 1), which is nearly non-existent in model solutions. The patterns of deep snow over MYI are well represented in the models but do not extend along the Canadian coastline towards the Beaufort Sea, as is the case for ASD. Overall, the models tend to overestimate snow depths nearly everywhere compared to the ASD data. The ASD data climatic mean is lower (12.3 cm), and zones of thin snow layers (e.g. around 120° E) are less pronounced in models (except by MERCATOR, where snow depth is low everywhere). Models also present higher variability. However, in all datasets the maximum variability occurs between 0 and 90° E, which is relevant because of the proximity to the open ocean and the existence of strong meteorological events (e.g. Semenov et al., 2019; Dong et al., 2019). Although investigating what causes these discrepancies is beyond the scope of this article, the strong sensitivity of the models to the reanalysis snowfall forcing data (e.g. Boisvert et al., 2018; Petty et al., 2018) could play a predominant role here.

### 4.2.2 Results in the Antarctic

ASD is currently the only publicly available altimetric snow depth product in the Antarctic. One major advantage compared to the Arctic is that the SARAL orbit does not affect the data coverage since all the Antarctic sea ice is below 81.5° S. Similarly to what has been done in the previous section for the Arctic, Fig. 5 compares ASD with the AMSR2-NSIDC data and the GIOMAS and the MERCATOR sea ice model. The statistics are summarised in Table A2.

In spite of a comparable climatic mean and some coherent spatial distributions such as the transport of snow in the Weddell Gyre, ASD and AMSR2-NSIDC data show significant discrepancies. Unlike what we have observed in the Arctic, ASD data are relatively smooth with a weak spatial variability. In comparison, AMSR2-NSIDC data are characterised by strong contrasts between large snow depth patterns in the Weddell Sea and to a lesser extent in the western part (from 150° E to 90° W) and areas with nearly no snow (eastern part, between 0 and 150° E). As in the Arctic, AMSR2-NSIDC annual and interannual variability is lower and much more localised than in the ASD data. ASD and AMSR2-NSIDC re-
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Comparisons of satellite snow depth estimations in the Arctic. (a, b, c) Bi-monthly March and April “climatological” snow depth mean maps calculated over the period 2014–2018 for the ASD data (a, d), the DuST data (b, e) and the AMSR2B data (c, f). (d, e, f) The standard deviation between all the data used to calculate the mean snow depth maps (equivalent to the mean interannual variability computed with March and April over the same 2014–2018 time period).</figcaption></figure>

produce well the thin snow pattern of the Ross Ice Shelf. One other relevant difference between these two datasets is the systematic decrease in snow depth in October in the AMSR2-NSIDC data (see Table A2).

The presence of thicker snow in East Antarctica for ASD is consistent with Worby et al. (2008b) in that comparisons with ARISE in situ data have shown radiometric measurement snow depth underestimations. Wet-snow conditions due to flooding might be responsible for radiometric brightness temperature contrasts. The low variabilities in the Weddell Sea for ASD were not expected since winter snow properties are extremely variable (Massom et al., 1997). In addition the large bias with AMSR2-NSIDC raises questions. Since the Ku band is supposed to better penetrate in cold and dry snow (Willatt et al., 2011), one hypothesis is that the presence of saline and warm moist basal snow layers, even in the absence of flooding (Massom et al., 2001; Perovich and Richter-Menge, 1994), lead to ASD underestimations. For AMSR2-NSIDC, the snow depth retrieval algorithm is probably not well adapted for rougher snow that can be compared to MYI in the Arctic.

With a climatic spatial mean of about 34.6 cm, GIOMAS simulates everywhere very high snow depths compared to the other solutions. Considering what is provided in the literature (e.g. Worby et al., 2011), these values clearly seem to

be overestimated. With the exception of the Weddell Sea, the MERCATOR model and ASD present comparable snow patterns with substantial coastal values not seen in the AMSR2-NSIDC data. However, the MERCATOR model displays thicker snow everywhere but in the Weddell Sea, while the presence of deep snow is a well observed feature, which is supported by ASD (e.g. Massom et al., 1997; Eicken et al., 1994). Because of highly dynamic weather conditions, with persistent strong winds in the Antarctic, snow thickness distribution is not directly related to snowfalls (Massom et al., 2001). Conversion from precipitation to snow depth is then very different to the Arctic, and Antarctic snow pack is not a uniform slab. These features should partly explain model difficulties and differences between the two hemispheres.

## 4.3 Comparison with in situ and airborne data

### 4.3.1 OIB

OIB data are compared to all the datasets presented in the previous sections. We only consider OIB observations below 81.5° N from the 2014–2018 campaigns. Figure 6 shows an along-track example from the OIB flight on 29 March 2015, and Fig. 7 presents scatter plots comparing the various snow depth products with the OIB data over the entire 2014–2018
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Comparisons between the ASD data and model snow depth products (PIOMAS, NESOSIM, MERCATOR) in the Arctic. The first column presents the “climatic mean” snow depth maps, the second column the “annual variability” and the third column the “interannual variability” as presented in Sect. 4.1. These diagnoses are calculated from March 2013 to April 2019 for ASD (first row), PIOMAS (second row) and MERCATOR (fourth row) and from March 2013 to April 2017 for NESOSIM (third row).</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 5.** Comparisons between ASD data and other snow depth estimations in the Antarctic (AMSR2-NSIDC, MERCATOR, GIOMAS). The first column presents the “climatic mean” snow depth maps, the second column the “annual variability” and the third column the “interannual variability” as presented in Sect. 4.1. These diagnoses are calculated from May 2013 to October 2019 for the ASD data (first row), the AMSR2-NSIDC (second row) and the MERCATOR model (third row) and from May 2013 to October 2016 for the GIOMAS model data (fourth row).</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Along-track comparison between the snow depth products and the OIB snow radar data acquired on 29 March 2015 in the Arctic. The red envelope refers to the ASD uncertainties. The geographical map specifies the geographical location of this particular 1 d OIB track.</figcaption></figure>

time period (including all spring OIB campaigns within this period). Statistical results are summarised in Table 1.

One striking feature is the good consistency between ASD and OIB in terms of magnitude and spatial variability. The OIB data are almost within the ASD envelope of uncertainties (in shaded red) at all times. These consistencies between OIB and ASD data are also demonstrated in several other OIB tracks provided in the Supplement (Figs. S14 to S17). The DuST Ka–Ku data present an along-track variability similar to ASD but tend to overestimate snow depth compared to OIB. This reinforces the hypothesis that the main difference between ASD and DuST is a bias. The AMSR2B data also reproduce the large-scale snow depth variability quite well but overestimate at a level comparable with DuST when compared to OIB. The DuST and AMSR2B overestimations are quite consistent with Kwok et al. (2017), who found that this OIB product tends to have thinner snow than the NOAA Wavelet Airborne Snow Radar dataset, which has been used for re-calibrations. In fact, it is important to be aware that various OIB datasets have been produced, leading to a variety of snow estimations that can reach 7 cm on FYI and 12 cm on MYI (Kwok et al., 2017). The variabilities in snow properties from OIB daily basis data to the monthly means of the other datasets certainly explain some of the deviations.

Models hardly represent the along-track variability and do not reproduce small-scale patterns. A resolution that is too coarse could partly explain this feature. As expected, the MERCATOR model snow thicknesses are far below the others, whereas the PIOMAS and NESOSIM solutions provide comparable results (see also Table 1). As expected, the W99m climatology strongly overestimates snow depths and is clearly a less optimal solution in this case.

### 4.3.2 IMB

The comparison between snow depth estimations and the IMB data for the winters of 2013–2014, 2014–2015 and 2015–2016 is presented in Fig. 8. Statistical results are summarised in Table 2. Except for AMSR2B, which is only available for the months of March and April, all comparisons are performed using the exact same spatial coverage considering only data below 81.5° N. Note that only very limited IMB data are available south of 81.5° N for the winter of 2015–2016.

Figure 8 and Table 2 show how the satellite-derived mean values are closer to the in situ observations than the model. However, the seasonal changes captured in the IMBs are not reproduced well by the satellite products or the models. This could reflect sampling differences as the IMB observations consist of hourly point data that are not representative of the wider averages of the gridded satellite or model products. DuST and ASD along-track variability is once again very similar, with DuST being overall higher, highlighting the very likely bias between these two products.

### 4.3.3 CryoVEx

Figure 9 presents the snow depth products projected along the CryoVEx flight track measured on 31 March 2017. This figure gathers the CryoVEx snow depth estimations calculated from (a) the difference between the Ku-band ASIRAS radar freeboard and the ALS laser (ALS/ASR) total freeboard and from (b) the Ku-band ASIRAS radar freeboard and the Karen (Ka/ASR) total freeboard. Note that because of SARAL's spatial coverage limitation (latitude < 81.5°) and the absence of the SARIN mode in the CS-2 ESA Baseline-B GOP (used to construct the ASD data; see Sect. 2), only the track of the 31 March 2017 CryoVEx survey of Baffin Bay can be compared with the ASD data.
<figure><img src="image_12.png" /><figcaption>**Figure 7.** Scatter plots between the snow depth products and OIB snow radar data considering the 2014 to 2018 missions. The red line is the linear fitting line. For NESOSIM, 2018 is not included.</figcaption></figure>

<table><caption>**Table 1.** Statistics of the comparison between OIB snow radar data and the various snow depth products considering the OIB spring campaigns of 2014–2018.</caption><thead><tr><th></th><th>OIB</th><th>ASD</th><th>DuST</th><th>PIOMAS</th><th>MERCATOR</th><th>NESOSIM</th><th>W99m</th><th>AMSR2B</th></tr></thead><tbody><tr><td>Mean (cm)</td><td>18.6</td><td>18.5</td><td>24.8</td><td>22.9</td><td>11.5</td><td>21.1</td><td>33.9</td><td>22.6</td></tr><tr><td>SD (cm)</td><td>7.5</td><td>6.9</td><td>6.5</td><td>10.4</td><td>4.4</td><td>8.5</td><td>1.5</td><td>5.8</td></tr><tr><td>Mean bias (cm)</td><td>—</td><td>≈0.0</td><td>-6.2</td><td>-4.3</td><td>7.1</td><td>-2.6</td><td>-15.03</td><td>-4.0</td></tr><tr><td>RMSE (cm)</td><td>—</td><td>5.9</td><td>8.7</td><td>10.2</td><td>9.8</td><td>8.4</td><td>17.0</td><td>7.1</td></tr><tr><td>R</td><td>—</td><td>0.66</td><td>0.62</td><td>0.50</td><td>0.45</td><td>0.51</td><td>0.16</td><td>0.63</td></tr></tbody></table>

Ka/ASR estimation exhibits very thin snow thickness (<10 cm). Although it might be expected in this area (Landy et al., 2017), this solution still contains unrealistic negative values due to the fact that ASIRAS and KAREN freeboards are nearly equivalent over FYI (without negative values). The main difference between snow depths calculated from the laser and the Ka-band-radar snow depth estimations seems to be a bias. Except around 71.3° N, where a specific event may have occurred, the ASD snow product tracks the magnitude of the CryoVEx airborne data, with ALS/ASR and Ka/ASR acting respectively as an upper and lower bound and being within ASD's range of uncertainty (in red shading). However, the along-track variability is not well reproduced by any dataset. The MERCATOR and NESOSIM models also show

similar magnitudes as the CryoVEx snow depth estimations (ALS/ASR and Ka/ASR), while PIOMAS and the W99m climatology clearly overestimate. AMSR2B and DuST overestimate snow depth, while the DuST along-track variability is comparable to ASD.

## 5 Towards an ASD snow depth climatology

In spite of the consistency of the ASD data highlighted in the previous sections from comparisons with airborne and in situ data, the limited temporal coverage (only available post-2013) remains an important limitation to the use of these data. Considering the low interannual snow precipita-
<figure><img src="image_13.png" /><figcaption>**Figure 8.** Along-track comparison between the snow depth products and the IMB in situ data. (a) IMB buoys compared with satellite data, both passive and active, and (b) IMB buoys compared with model data. The right column shows geographical maps from data for all three periods (2013–2014, 2014–2015, 2015–2016). The maps present the geographical location of the three IMB buoys.</figcaption></figure>

<table><caption>**Table 2.** Statistics of the comparison with IMB data and the various snow depth products considering data acquired during 2013–2016.</caption><thead><tr><th></th><th>IMB</th><th>ASD</th><th>DuST</th><th>PIOMAS</th><th>MERCATOR</th><th>NSIM</th><th>W99m</th><th>AMSR2B</th></tr></thead><tbody><tr><td>Mean (cm)</td><td>21.6</td><td>19.14</td><td>25.51</td><td>19.37</td><td>7.61</td><td>18.32</td><td>28.03</td><td>24.88</td></tr><tr><td>SD (cm)</td><td>2.2</td><td>3.17</td><td>2.90</td><td>4.68</td><td>2.17</td><td>3.68</td><td>5.16</td><td>1.05</td></tr><tr><td>RMSE (cm)</td><td>—</td><td>4.73</td><td>5.50</td><td>5.11</td><td>14.32</td><td>4.96</td><td>8.2</td><td>3.71</td></tr></tbody></table>
<figure><img src="image_14.png" /><figcaption>**Figure 9.** Along-track comparison, in the Arctic, between the snow depth products (models and satellites) and the CryoVEx snow depth data calculated from the Ka-band KAREN radar (blue) and the ALS laser (dark blue) with the Ku-band ASIRAS radar during the 31 March 2017 campaign. The map specifies the geographical location of the CryoVEx track.</figcaption></figure>

tion variability compared to that of the spatial variability (e.g. Figs. 2 and 4), one solution would be to develop an ASD climatology. For that purpose, we have constructed a preliminary altimetric snow depth climatology in the Arctic by averaging all the ASD snow depth maps of each month during the 2013–2019 period (designated as ASD-clim). An equivalent climatology could also be constructed for the Antarctic, but the lack of validation prevents its validity to be demonstrated at present. To demonstrate the relevance of this climatology in the Arctic for the years prior to 2013, the ASD-clim data are projected on all the tracks of the four OIB missions occurring between 2009 and 2012 and are presented in Fig. 10.

Figure 10 shows that the ASD climatology would be a better solution than W99m, at least for the CS-2 time period. This is highlighted with the lower bias (–0.1 m compared with –0.14 m), the higher correlation coefficients (0.71 compared with 0.06) and the lower RMSE (0.06 m compared with 0.17 m). Furthermore, even as a climatology, we obtained a correlation between the ASD-clim and OIB of the same order of magnitude as obtained by direct correlation between ASD and OIB measurements (Fig. 7). However, as future work, a more refined computation should be investigated. In particular, an important point would be to provide a relevant snow product representing the Envisat era. Since the W99m can be considered as relatively consistent until the year 2000, a mixed ASD and W99m product (ASD-W99m), including, for instance, a time dependency, could be considered. Also, it could be relevant to take into account the sea ice type interannual variability as is done for W99m.

Another critical point of the ASD data (for instance for sea ice volume estimations) is the reduced spatial coverage due to the SARAL orbit. While during the time period of Envisat, this does not become an issue since its orbit is equivalent to that of SARAL; the absence of data north of 81.5°N is a

major limitation to SIT estimations. A combination of ASD with W99m and/or model data could be done to extrapolate the data. Potentially, one could also investigate a combination of ASD with W99m and/or models to extrapolate the data. In the Antarctic, we do not encounter the same issue since the orbit does not reduce the spatial coverage, and FYI is generally considered to be the only sea ice type present in the Southern Ocean. Then, such a climatology would be a solution to compute SIT estimations over the Envisat and CS-2 time period provided that the temporal variability in precipitation is investigated. Note that experimental SIT estimations (Hendricks et al., 2018) have already been done in the Antarctic in the framework of the Sea Ice Climate Change Initiative (SI-CCI). These data used a snow climatology elaborated from AMSR-E and AMSR2 data in a nearly equivalent manner to what has been presented here.

## 6 Impacts of snow depth on SIT estimation

Since there is still no validated freeboard product in the Antarctic, this analysis has only been performed in the Arctic. From the CS-2 radar freeboard product presented in Guerreiro et al. (2017), extended to 2019, we have computed SIT estimates using the various snow depth datasets presented in this study. They are calculated from the CS-2 L1b ESA Baseline-C SAR waveforms (Bouffard et al., 2018b) with the methodology presented in Sect. 2. Freeboard heights are converted into SIT assuming hydrostatic equilibrium between snow-covered sea ice and the ocean (Eq. 8).

$$SIT = \frac{\rho_{wfbice} + \rho_s sd}{\rho_w - \rho_i} \qquad (8)$$
<figure><img src="image_15.png" /><figcaption>**Figure 10.** Scatterplots comparing the ASD and the W99m climatologies with OIB by considering all the tracks of the 2009, 2010, 2011 and 2012 OIB missions.</figcaption></figure>

As in Ricker et al. (2014), we assume the seawater density ($\rho_w$) is set to 1024 kg/m$^3$. Consistent with the approach of Laxon et al. (2013), sea ice density ($\rho_i$) is set to 882 kg/m$^3$ for MYI and 917 kg/m$^3$ for FYI (Alexandrov et al., 2010). Furthermore, the snow density ($\rho_s$) is set to 300 kg/m$^3$. A distinction between FYI and MYI follows the EUMETSAT OSI SAF sea ice type classification. The sea ice freeboard fb$_{ice}$ is calculated from the radar freeboards, taking into account the reduction in radar velocity within the snow layer, which depends on snow depth data following Eq. (3). SIT monthly maps are calculated for each snow depth product presented in previous sections. These monthly maps are provided for the winter of 2015 in the Supplement (Figs. S18 to S24). Figure 11 presents the annual mean time series of the spatial mean of all these SIT solutions (except for AMSR2B which is only calculated in March and April).

Except for the MERCATOR model (which has shown to be overall significantly low), using the ASD data provide the lowest pan-Arctic mean SIT estimations. The highest SIT estimations are obtained with the DuST product (and W99m), highlighting the impact of the bias previously identified between these two solutions (ASD and DuST). Considering the various snow products together (except for the MERCATOR model), the spatial mean SIT annual mean values are distributed between 1.2 and 1.6 m. This leads to a pan-Arctic SIT ensemble mean of 1.4 m (mean of the dotted blue line in Fig. 11) associated with a mean maximum deviation of 29 cm (mean of differences between the maximum and minimum dotted blue lines). This mean maximum deviation, which is approximately the mean difference between ASD and DuST, represents about 20 % of the mean SIT. Considering the good consistency between ASD and the validation snow depth data (demonstrated in Sect. 4.3), actual SIT products computed using W99m should tend to significantly overestimate sea ice volume in the Arctic. For instance, taking a deviation of

<figure><img src="image_16.png" /><figcaption>**Figure 11.** Arctic annual mean time series (2014–2019) of the averaged SIT solutions computed using various snow depth products (models and active/passive satellite-based). Spatial coverage is defined by the smallest sea ice extent. The dotted blue lines define the annual mean, minimum and maximum values between all solutions (except for MERCATOR). The AMSR2 data are not represented since they are not available over the period.</figcaption></figure>

about 0.3 m (≈ the mean deviation between ASD and W99m in 2015) and a mean Arctic sea ice extent of 13 × 10⁶ km² (approximately the mean sea ice extent of 2015), the resulting SIT overestimation would lead to a sea ice volume excess of about 4 × 10³ km³, which represents about 3.5 × 10¹⁵ L of freshwater. Such an uncertainty is not negligible compared to the 8000 ± 2000 km³ change in the western Arctic freshwater content between 1995–1996 and 2009–2010 relayed in Giles et al. (2012). To give an illustrating example, this uncertainty of 3.5 × 10¹⁵ L of freshwater would be roughly
equivalent to 6 months of the Amazon water discharge (considering a debit of about $2 \times 10^8$ L/s). Deviations of only a few centimetres of snow depth ($\approx 6$ cm of mean deviation between ASD and W99m snow depths in 2015; see Fig. 1) have strong impacts on the freshwater budget.

Although this provides a relevant overview, such large-scale pan-Arctic means do not show the spatial and temporal variability in the various SIT solutions. For instance, while DuST and W99m provide comparable SIT pan-Arctic mean estimations, the spatial and temporal information included in these two datasets are very different (e.g. Figs. 2 and 3). To better characterise the impacts of snow depth on SIT, we propose to consider the deviations between these SIT estimates for each month of each year and at each grid point. Such an ensemble of SIT solutions will allow us to characterise the level of uncertainty due to discrepancies between the different snow depths. For the month m of year y, the standard deviation ($SD_{y, m}(SIT)$) and the maximum deviation ($maxdev_{y, m}(SIT)$) maps are calculated with Eqs. (9) and (10). $SIT_{y, m, p}$ is the SIT solution provided by the snow product p, and $SIT_{y, m}$ is the mean map between all these SIT estimations. $N_p$ is the number of snow depth products, and $max_p(SIT_{y, m, p})$ and $min_p(SIT_{y, m, p})$ are respectively the maximum and minimum monthly maps between the SIT solution of the snow depth products.

$$SD_{y, m}(SIT) = \sqrt{\frac{1}{N_p} \sum_{p} (SIT_{y, m, p} - \overline{SIT_{y, m}})^2} \quad (9)$$

$$maxdev_{y, m}(SIT) = \max_{p}(SIT_{y, m, p}) - \min_{p}(SIT_{y, m, p}) \quad (10)$$

For each winter month, interannual mean, minimum and maximum maps of these two statistics (standard deviation $SD_{y, m}(SIT)$ and the maximum deviation $maxdev_{y, m}(SIT)$) can then be computed from Eqs. (11) to (16), where $N_y$ is the number of years.

$$\overline{SD_m} = \frac{1}{N_y} \sum_{y} (SD_{y, m}(SIT)) \quad (11)$$

$$\overline{maxdev_m} = \frac{1}{N_y} \sum_{y} (maxdev_{y, m}(SIT)) \quad (12)$$

$$\max(SD_m) = \max_{y}(SD_{y, m}(SIT)) \quad (13)$$

$$\max(maxdev_m) = \max_{y}(maxdev_{y, m}(SIT)) \quad (14)$$

$$\min(SD_m) = \min_{y}(SD_{y, m}(SIT)) \quad (15)$$

$$\min(maxdev_m) = \min_{y}(maxdev_{y, m}(SIT)) \quad (16)$$

These statistics quantify the SIT uncertainty, caused only for snow depth, assuming that deviations between snow depth solutions refer to the snow depth level of uncertainty. We define two cases: (1) “obs snow products”, which considers ASD, DuST, W99m and AMSR2B, and (2) “all snow

products”, which includes in addition the PIOMAS, MERCATOR and NESOSIM model solutions. Note that ASD-clim is not used in calculations in order to avoid the redundancy with ASD data. Considering “obs snow products”, Fig. 12 presents maps of the metrics presented in Eqs. (11) to (16) for April. In addition, the monthly minimum, mean and maximum maps over the period 2014–2019 are presented for April in the Supplement (Fig. S25). All the results are summarised in Table A3 which presents, for each winter month, the spatial mean value of these maps. The spatial mean of the monthly mean SIT maps ($\overline{SIT}_{y, m, p}$) are also included.

The highest impacts on SIT are generally located near coastal zones and marginal ice zones. In these areas, maximum standard deviation can reach up to 50 cm and maximum deviation to 80 cm when considering “obs snow products”. Due to larger differences between model and observations, in particular MERCATOR, these values are higher with “all snow products” (see Table A3). Except for these specific locations that mostly account for maximum values, the standard deviation maps are rather uniform. This suggests that spatial mean values (presented in Table A3) are a relevant indicator.

The large deviations between minimum and maximum (standard deviation and/or maximum deviation) maps indicate significant variations in the SIT level uncertainty due to snow depth variability from one year to another. On the contrary, spatial mean standard deviations do not vary significantly in between months. Considering “obs snow products” and all winter months, we obtained a SIT mean standard deviation of $\approx 20$ cm (14%) $\pm 11$ cm (31–9 cm interval; see Table A3) and mean maximum deviation of $\approx 49$ cm (35%) $\pm 25$ cm. These values increase to $\approx 27$ cm (20%) $\pm 9$ cm for the standard deviation and $\approx 76$ cm (58%) $\pm 25$ cm for the maximum deviation when we take into account “all snow products”. Therefore, compared to the deviation (roughly 30 cm) that can be extracted from the difference between ASD (red line) and the SIT max line (dotted blue) in Fig. 11, we obtain higher maximum deviations using this methodology. It suggests that considering large-scale global means could tend to underestimate impacts of snow depth on SIT and therefore the SIT uncertainty. Note that in addition to a more refined methodology, taking into account AMSR2B, which can reach high values, might have a small impact here as well.

Besides providing a relevant quantification of the impacts of snow depth on SIT, such a “multi-observation” approach also provides a space and time variability in uncertainties. It also avoids the necessity of relying on the strong, and usual, assumption of Gaussian propagation of de-correlated errors. These uncertainty estimates are very important for data assimilation, in which the characterisation of observation error is of main importance to constrain models (e.g. Stewart et al., 2008; Bunzel et al., 2016; Janjić et al., 2018). Note that there are a growing number of research programmes and international initiatives which aim at coordinating efforts
<figure><img src="image_17.png" /><figcaption>**Figure 12.** Interannual minimum (min($\overline{SD_m}$), first column), mean ($\overline{SD_m}$, second column) and maximum (max($\overline{SD_m}$), third column) maps of the standard deviation (first row) and the maximum deviation (second row) for the month of April (m = 04). The snow products used to compute the maps are ASD, DuST, AMSR2B and W99m. These maps correspond to the case “obs snow products”.</figcaption></figure>

for better characterising uncertainties and provide traceable quality indicators to sea ice earth observation (EO) products (e.g. World Meteorological Organisation (WMO) workshop at AWI, European Union/European Commission (EU/EC) Polar Science Week).

## 7 Discussions and conclusions

In this study, we have presented a Ka–Ku snow depth product computed, in both hemispheres, from the data obtained by altimeters on board the satellites SARAL (Ka-band) and CS-2 (Ku-band). This product has been compared over the period 2013–2019 with the other Ka–Ku snow depth product (DuST, only in the Arctic), AMSR2 passive radiometer snow depth estimations, W99m (only in the Arctic) and snow depth solutions of the models PIOMAS (GIOMAS in the Antarctic), MERCATOR and NESOSIM (only in the Arctic). Thereafter, all these products have been evaluated, in the Arctic, against OIB airborne data, IMB snow buoys and CryoVEx airborne campaign data, including Ka–Ku, as well as ALS and Ku, snow depths. In the Antarctic, the lack of validation data remains a major obstacle to the assessment of sea ice products. It is currently a primary focus of the ESA CSAO+ project.

In the Arctic, we observed a good agreement between the two Ka–Ku products (ASD and DuST) in terms of spatial dis-

tributions and annual and interannual variability, but DuST is always high biased by about 6.5 ± 0.5 cm (see Table A1). Although this is likely due to the re-calibration of the DuST product with OIB data mainly measuring thick snow, this bias probably also indicates variabilities in Ku-band penetration related to snow properties. Only the surface roughness is considered in the ASD calculation, but the impact of other snow properties should be analysed and considered.

The AMSR2B (and also AMSR2-NSIDC over FYI) data also tend to overestimate snow depths compared to ASD in the Arctic. This overestimation is lower than DuST, on average of about 3 cm, and does not seem to only be a bias since AMSR2B spatial patterns present noticeable differences. We reiterate that comparisons with AMSR2B have only been conducted in March and April because the data are not available over MYI for the other months. This absence of variability except in the MYI zone where AMSR2B has been re-calibrated is a surprising feature that needs to be considered when using these data.

Comparisons with models in the Arctic have highlighted that the MERCATOR model is always biased low when compared to the other products and simulates limited snow depth variability. This is mainly due to the use of the ERA-Interim forcing fields that underestimate precipitation. The impacts of other parameters such as the albedo or the amount of precipitation falling into the ocean should be investigated. In
contrast, NESOSIM and PIOMAS spatial patterns are con-
sistent with ASD and have a bias on average of a few cen-
timetres high. We found comparable representations in PI-
OMAS and NESOSIM, although NESOSIM presents higher
snow depths over MYI and higher variability. In spite of con-
sistencies, these two models hardly represent the MYI sea ice
drift in the Beaufort Gyre that can be observed in satellite ob-
servations. Note that we have not considered in this study the
recent SnowModel-LG model that should consistently cap-
ture the Arctic snow depth over sea ice spatial variability
(Stroeve et al., 2020; Liston et al., 2020). We refer to Zhou
et al. (2021) for a comparison between the SnowModel-LG
and various snow depth observation products.

In evaluation and for comparison of the snow depth prod-
ucts, we have compared OIB, CryoVEx and in situ IMB
snow depths over the Arctic. Comparisons with OIB have
shown satisfying consistencies with the ASD estimations.
We found a RMSE of about 6 cm with a correlation coeffi-
cient of 0.66. By comparison, the DuST product has a RMSE
of 8.7 cm and AMSR2B of 7 cm. Note that this higher RMSE
for DuST, in spite of a re-calibration, probably comes from
the OIB version used for re-calibrating. Indeed, DuST is re-
calibrated against the OIB NOAA Wavelet Airborne Snow
Radar dataset which is expected to present higher snow depth
values. In addition, further analysis could be done in order to
link the space and time variability in these biases with snow
properties The various snow depth products are less consis-
tent when comparing with IMB, although ASD mean val-
ues present a consistent order of magnitude. Due to the data
coverage, validation with CryoVEx has been limited to one
single track (where both Ku-band radar (ASIRAS), Ka-band
radar (KAREN) and lidar scanner (ALS) were available). We
observed that KAREN and ASIRAS snow depth estimations
are low, including a non-negligible amount of negative val-
ues. ASD data always remain in between the laser-Ku (ALS-
ASR) and Ka-Ku (Ka/ASR) CryoVEx snow depth estima-
tions.

Since the results of such comparisons with validation data
can vary from one methodology to another (e.g. grid sizes,
smoothing kernel, dataset versions), they do not aim to as-
sess the best snow depth product. However, they demonstrate
how ASD provides a relevant snow depth solution, in good
agreement with several validation data, and that ASD allows
for the characterisation of the deviations between the differ-
ent snow depth products. A more refined comparison with
the CryoVEx data (including various tracks) is mandatory to
understand the relative impacts of roughness and penetration
and their link with snow properties. For instance, Willatt et al.
(2011) show that the Ku-band dominant scattering surface
can significantly vary with snow temperature, with a reduced
penetration when temperatures increase. This feature could
enhance an underestimation of ASD data with the warming
of temperatures due to climate change. It also points out the
difficulty to retrieve snow depth from altimetry beyond the

6 months of winter as variations in snow properties are most
important in the summer period.

In the Antarctic, for satellite-based observations, spatial
mean snow depth values of AMSR2-NSIDC and ASD are
quite close. However, we found that AMSR2-NSIDC esti-
mations are very low in zones of thin snow, while it tends to
strongly overestimate snow depths in areas of thick ice. Such
patterns of very thin snow have not been observed in the Arc-
tic. Also patterns of thick snow are very likely due to the re-
trieval algorithm that is not adapted for thick ice. Regarding
models, over the Antarctic, the MERCATOR model simu-
lates mean values in better agreement with ASD. However,
thicker snow patterns and the clockwise drifts in the Weddell
Sea are not captured by the model. The GIOMAS model sim-
ulates very high snow depths nearly everywhere, likely lead-
ing to unrealistic snow depth representations. Sea ice mod-
els are generally very dependent on atmospheric forcings,
and differences in snow accumulation in models can lead to
very different snow depths. Model biases are largely due to
the lack of consistent observations, which has, so far, limited
the tuning of snow parameters such as albedo or snow diffu-
sion coefficients (Chevallier et al., 2017; Uotila et al., 2019).
These comparisons highlight snow properties in the Antarc-
tic that are very different than in the Arctic. Such opposite
model behaviours in the two hemispheres argue that a single
model tuning of parameters may not be adapted for a global
configuration. In the context of data assimilation, the positive
impacts of the integration of SIT (Fiedler et al., 2021; Mas-
sonnet et al., 2015) could also be enhanced using ASD either
by improving the SIT observations or by enabling the simul-
taneous use of freeboard and snow depth data to constrain the
models.

An important limitation of the ASD snow depth data is its
temporal coverage, imposed by SARAL which has only been
available since 2013. To circumvent this limitation, we have
shown that a simple monthly climatology, constructed as the
average of ASD snow depth maps over the 2013–2019 pe-
riod, could be a relevant option. More specifically, we have
demonstrated that such a climatology could provide a more
reliable solution than the W99m climatology, which strongly
overestimates snow depths and is useless outside the central
Arctic due to extrapolation. Indeed, comparison between the
ASD climatology and the 2009–2012 OIB missions shows
good consistency (comparable results of similar magnitude
as when comparing the actual ASD product with OIB cam-
paigns from 2014–2018), while W99m is on average biased
about 14 cm high with an associated RMSE of 17 cm (when
comparing with OIB). Such a climatology should be refined
in order to provide relevant snow depth representations for
at least the Envisat and CS-2 time period. For instance, the
possibility of constructing a climatology using both ASD,
AMSR and W99m information should be investigated.

A better snow depth representation significantly improves
SIT estimates. In this context, this study has investigated the
impacts of different snow depth products on SIT. Meanwhile,
we have also presented a methodology to characterise the SIT
level of uncertainty due to the snow depth deviations between
products. The approach constructed an ensemble of SIT solu-
tions using the different snow depth products over the 2013–
2019 period, thus providing an estimation of the SIT level
of uncertainty due to snow depth discrepancies and its space
and time variability. Note that equivalent methodologies are
frequently used in multi-model approaches. Taking only into
account the satellite products AMSR2B, ASD, DuST and the
W99m climatology over the 2013–2019 period, we found a
spatially averaged mean standard deviation of 20 cm (14 %
of the pan-Arctic mean SIT) and mean maximum deviation
of 49 cm (35 % of the pan-Arctic mean SIT) between the
different SIT estimations. Deviations between SIT estima-
tions reached up to 77 cm (55 % of pan-Arctic mean SIT) in
marginal and coastal zones.

We have highlighted in this study that the DuST and ASD
products, albeit computed from data acquired by the same
sensors at the same time, present some significant differences
due to the methodology applied. The main difference is that
for ASD we compare LRM data (PLRM for CS-2), while
DuST compares LRM data with SAR data. The advantage
of comparing LRM data together is that they have compara-
ble footprint sizes, by which the difference reduces the ef-
fects of surface roughness to highlight the effects of volume
backscatter and penetration. DuST, on the other hand, relies
on the calibration of the measurements on OIB without tak-
ing into account effects of the surface roughness. In fact, dif-
ferences in surface roughness and the effect it has on radar
penetration and retrieval of surface elevation is currently one
of the other main contributors to the SIT uncertainty budget
(e.g. Hendricks et al., 2010; Landy et al., 2020), which still
needs to be further investigated. We have only considered
the roughness as mainly depending on the size of the foot-
print, but the backscattering coefficient variability is not as
trivial and depends on several snow physical properties such
as density, grain size, salinity and/or moisture (e.g. Lund-
berg et al., 2006; Adodo et al., 2018; Nandan et al., 2020).
That is, specular surfaces have stronger nadir scattering of
the radar pulse, affecting the shape of the radar echo and,
in turn, the retrieval of the surface elevation. Empirical re-
trackers like TFMRA have been shown to be sensitive over
rougher surfaces (e.g. Ricker et al., 2014) In particular, com-
pared to physical retracking they have been shown to overes-
timate sea ice freeboards over rougher ice while also underes-
timating sea ice freeboard over smooth, thin ice (e.g. Laforge
et al., 2020; Landy et al., 2020). In turn, when assuming
that snow depth can be retrieved from the difference between
the Ku-band and Ka-band freeboards, the effect that surface
roughness, and particularly the Ku-band volume retrodiffu-
sion/penetration, has on the retrieval of the surface elevation
warrants further investigation. An analysis of these depen-
dencies is out the scope of this study, but it is necessary to
further use the ASD data to improve our understanding of
snow properties.

The use of two satellites, operating in different modes,
with different orbits, also leads to important uncertainties
and differences that cannot be neglected. The high-priority
candidate mission, CRISTAL, which shall require a dual-
frequency Ka–Ku SAR altimeter, could address some of
these issues. The fact of having coincident and collocated
SAR measurements obtained from the same platform, with
similar ground footprints for both frequencies, will enable
direct comparisons and a better understanding of the rela-
tive impacts of the surface roughness, radar penetration and
volume backscattering. In this context, ASD demonstrates
the possibilities of dual-frequency snow depth estimations,
which could be further improved using physical retracking
(as opposed to empirical re-trackers which have their own
limitations, e.g. Ricker et al., 2014). In addition, CRISTAL
will also mitigate the issue of the spatial coverage limitation
of 81.5° N imposed by SARAL in the Arctic. It is also impor-
tant to note that the CRISTAL mission is the only altimetric
mission planned to succeed CS-2 and avoid a gap of observa-
tions. Note that this Copernicus mission is not yet officially
selected. If it is not confirmed, sea ice thickness measure-
ment by radar altimetry above 81.5° N will end with the CS-2
mission. Prior to potential CRISTAL data, comparisons with
the very recent ICESat-2 and CS-2 snow depth estimations
(Kwok et al., 2020) computed from the difference between
laser and Ku-band measurements certainly provide impor-
tant first results. Moreover, the recent change in CS-2 orbit
to align with ICESat-2 (CRYO2ICE), providing collocated
measurements with ≈ 3 h temporal delay, also allows for the
investigation of Ku-band penetration and the discrepancies
that different radar footprints can introduce.
Appendix A: Appendices

<table><caption>**Table A1.** Statistical comparison of the different snow depth estimations in the Arctic over the 6 winters of the 2013–2019 period. AMSR2’s data are not indicated since only the months of March and April were available.</caption><thead><tr><th>Arctic</th><th>ASD</th><th>MERCATOR</th><th>PIOMAS</th><th>NESOSIM</th><th>W99m</th><th>DuST</th></tr></thead><tbody><tr><td rowspan="6">Interannual mean (cm)</td><td>11</td><td>10.3</td><td>2.5</td><td>7.8</td><td>9.7</td><td>14.9</td><td>17.1</td></tr><tr><td>12</td><td>10.6</td><td>4</td><td>10</td><td>11.9</td><td>15</td><td>17.5</td></tr><tr><td>01</td><td>11.4</td><td>6.5</td><td>14</td><td>14.8</td><td>17.1</td><td>18.1</td></tr><tr><td>02</td><td>12.4</td><td>8.8</td><td>17.5</td><td>18</td><td>18.8</td><td>19.1</td></tr><tr><td>03</td><td>13.8</td><td>10.8</td><td>20.9</td><td>20.7</td><td>20.4</td><td>20.3</td></tr><tr><td>04</td><td>15.3</td><td>11.6</td><td>23</td><td>22.38</td><td>19</td><td>20.9</td></tr><tr><td>Climatic mean</td><td>12.3</td><td>7.4</td><td>15.5</td><td>16.2</td><td>17.5</td><td>18.8</td></tr><tr><td rowspan="6">Interannual variability (cm)</td><td>11</td><td>4.8</td><td>1.4</td><td>3</td><td>2.4</td><td>2</td><td>4.3</td></tr><tr><td>12</td><td>4.6</td><td>2</td><td>3.5</td><td>3</td><td>1.8</td><td>4</td></tr><tr><td>01</td><td>4.5</td><td>2.6</td><td>4</td><td>3.1</td><td>2.3</td><td>4.3</td></tr><tr><td>02</td><td>4.4</td><td>3</td><td>4.4</td><td>3.5</td><td>2.6</td><td>4.3</td></tr><tr><td>03</td><td>4.7</td><td>3.9</td><td>5.5</td><td>4.3</td><td>2.6</td><td>4.5</td></tr><tr><td>04</td><td>5.2</td><td>4.3</td><td>5.8</td><td>4.7</td><td>3</td><td>4.9</td></tr><tr><td>Mean ≡ MIV</td><td>4.7</td><td>2.9</td><td>4.4</td><td>3.5</td><td>2.4</td><td>4.4</td></tr><tr><td rowspan="6">Annual variability (cm)</td><td>2013–2014</td><td>4.2</td><td>3.3</td><td>5.5</td><td>5</td><td>4.3</td><td>4.3</td></tr><tr><td>2014–2015</td><td>4</td><td>3.1</td><td>5.8</td><td>5</td><td>4.1</td><td>3.9</td></tr><tr><td>2015–2016</td><td>4.1</td><td>2.7</td><td>4.7</td><td>4.1</td><td>3.8</td><td>4.2</td></tr><tr><td>2016–2017</td><td>4.2</td><td>3.9</td><td>6.1</td><td>4.9</td><td>4.4</td><td>4.2</td></tr><tr><td>2017–2018</td><td>4.2</td><td>4</td><td>6.9</td><td>x</td><td>3.6</td><td>3.8</td></tr><tr><td>2018–2019</td><td>4.2</td><td>3.9</td><td>6.5</td><td>–</td><td>3.2</td><td>–</td></tr><tr><td>Mean ≡ MAV</td><td>4.18</td><td>3.5</td><td>5.9</td><td>4.8</td><td>3.9</td><td>4.1</td></tr></tbody></table>

<table><caption>**Table A2.** Statistical comparison of the different snow depth estimations in the Antarctic over the seven winters of the 2013–2019 period. MIV signifies mean interannual variability, and MAV signifies mean annual variability.</caption><thead><tr><th>Antarctic</th><th>ASD</th><th>MERCATOR</th><th>GIOMAS</th><th>AMSR2</th></tr></thead><tbody><tr><td rowspan="6">Interannual mean (cm)</td><td>05</td><td>14.9</td><td>10.2</td><td>26.1</td><td>17.3</td></tr><tr><td>06</td><td>15.2</td><td>12.8</td><td>29</td><td>18</td></tr><tr><td>07</td><td>16</td><td>16.9</td><td>33.5</td><td>18.7</td></tr><tr><td>08</td><td>16.7</td><td>21.3</td><td>37.3</td><td>19.2</td></tr><tr><td>09</td><td>17.8</td><td>25.2</td><td>40.5</td><td>17.4</td></tr><tr><td>10</td><td>19.7</td><td>26.5</td><td>41.2</td><td>15.2</td></tr><tr><td>Climatic mean</td><td>16.7</td><td>18.8</td><td>34.6</td><td>17.6</td></tr><tr><td rowspan="6">Interannual variability (cm)</td><td>05</td><td>6</td><td>3.5</td><td>8.4</td><td>5.6</td></tr><tr><td>06</td><td>5.9</td><td>4.5</td><td>7.8</td><td>5.5</td></tr><tr><td>07</td><td>5.9</td><td>5.2</td><td>7.1</td><td>5.5</td></tr><tr><td>08</td><td>5.6</td><td>5.6</td><td>6.6</td><td>5.3</td></tr><tr><td>09</td><td>6.2</td><td>6.1</td><td>6.5</td><td>4.9</td></tr><tr><td>10</td><td>7.4</td><td>6.7</td><td>7</td><td>5.3</td></tr><tr><td>Mean ≡ MIV</td><td>6.2</td><td>5.3</td><td>7.2</td><td>5.3</td></tr><tr><td rowspan="6">Annual variability (cm)</td><td>2013</td><td>5.5</td><td>6.7</td><td>7.5</td><td>4.1</td></tr><tr><td>2014</td><td>5.2</td><td>6.6</td><td>6.6</td><td>4.3</td></tr><tr><td>2015</td><td>5.1</td><td>6.8</td><td>6.4</td><td>3.7</td></tr><tr><td>2016</td><td>5.5</td><td>5.4</td><td>9.2</td><td>4.2</td></tr><tr><td>2017</td><td>5.6</td><td>6.4</td><td>–</td><td>3.7</td></tr><tr><td>2018</td><td>5.5</td><td>6.7</td><td>–</td><td>3.8</td></tr><tr><td>2019</td><td>5.4</td><td>6.0</td><td>–</td><td>4.1</td></tr><tr><td>Mean ≡ MAV</td><td>5.4</td><td>6.4</td><td>7.4</td><td>4</td></tr></tbody></table>
<table><caption>**Table A3.** Statistical impact of the various snow depths on SIT in the Arctic over the 2013–2019 time period. The upper part considers all the snow depth products while the upper part consider only the products based on observations (ASD, DuST, AMSR2B and W99m). Note that AMSR2B SIT estimations on MYI are only taken into account for the months of March and April.</caption><thead><tr><th></th><th></th><th>m11</th><th>m12</th><th>m01</th><th>m02</th><th>m03</th><th>m04</th><th>Mean</th></tr></thead><tbody><tr><td rowspan="3">All snow products</td><td>Mean</td><td>0.94</td><td>1.08</td><td>1.25</td><td>1.42</td><td>1.59</td><td>1.65</td><td>1.32</td></tr><tr><td rowspan="3">SD</td><td>mean</td><td>0.27</td><td>0.26</td><td>0.25</td><td>0.26</td><td>0.27</td><td>0.28</td><td>0.27</td></tr><tr><td>min</td><td>0.20</td><td>0.19</td><td>0.18</td><td>0.18</td><td>0.18</td><td>0.18</td><td>0.18</td></tr><tr><td>max</td><td>0.35</td><td>0.34</td><td>0.34</td><td>0.34</td><td>0.37</td><td>0.41</td><td>0.36</td></tr><tr><td rowspan="3">Max dev</td><td>mean</td><td>0.74</td><td>0.74</td><td>0.71</td><td>0.72</td><td>0.78</td><td>0.84</td><td>0.76</td></tr><tr><td>min</td><td>0.55</td><td>0.54</td><td>0.50</td><td>0.49</td><td>0.51</td><td>0.53</td><td>0.52</td></tr><tr><td>max</td><td>0.98</td><td>0.97</td><td>0.95</td><td>0.97</td><td>1.09</td><td>1.23</td><td>1.03</td></tr><tr><td rowspan="3">Obs snow products</td><td>Mean</td><td>1.09</td><td>1.21</td><td>1.34</td><td>1.48</td><td>1.62</td><td>1.65</td><td>1.40</td></tr><tr><td rowspan="3">SD</td><td>mean</td><td>0.21</td><td>0.21</td><td>0.21</td><td>0.20</td><td>0.19</td><td>0.19</td><td>0.20</td></tr><tr><td>min</td><td>0.10</td><td>0.11</td><td>0.10</td><td>0.10</td><td>0.09</td><td>0.08</td><td>0.09</td></tr><tr><td>max</td><td>0.32</td><td>0.32</td><td>0.31</td><td>0.30</td><td>0.31</td><td>0.30</td><td>0.31</td></tr><tr><td rowspan="3">Max dev</td><td>mean</td><td>0.51</td><td>0.52</td><td>0.50</td><td>0.48</td><td>0.47</td><td>0.47</td><td>0.49</td></tr><tr><td>min</td><td>0.24</td><td>0.24</td><td>0.24</td><td>0.23</td><td>0.22</td><td>0.20</td><td>0.23</td></tr><tr><td>max</td><td>0.78</td><td>0.80</td><td>0.78</td><td>0.75</td><td>0.76</td><td>0.77</td><td>0.77</td></tr></tbody></table>

<table><caption>**Table A4.** Time period of snow depth products used in this article. Note that the current data availability may have changed since this article was written.</caption><thead><tr><th>Long name</th><th>Acronym</th><th>Reference</th><th>Data coverage</th></tr></thead><tbody><tr><td>Altimetric snow depth</td><td>ASD</td><td>Present paper</td><td>March 2013–April 2019</td></tr><tr><td>Dual-altimeter Snow Thickness</td><td>DuST</td><td>Lawrence et al. (2018)</td><td>March 2013–April 2018.</td></tr><tr><td>AMSR2 data of the University of Bremen</td><td>AMSR2B</td><td>Rostosky et al. (2018)</td><td>March 2013–April 2018</td></tr><tr><td>AMSR2 data from the National Snow and Ice Data Center</td><td>AMSR2-NSIDC</td><td>Meier et al. (2018)</td><td>March 2013–present</td></tr><tr><td>The Pan-Arctic (Global) Ice-Ocean Modelling and Assimilation System (v2.1)</td><td>P(G)IOMAS</td><td>Schweiger et al. (2011)</td><td>March 2013–April 2019</td></tr><tr><td>LIM-2 model in the PSY4V3 MERCATOR configuration</td><td>MERCATOR</td><td>Lellouche et al. (2018)</td><td>2013–2019</td></tr><tr><td>The NASA Eulerian Snow On Sea Ice Model (v1.0)</td><td>NESOSIM</td><td>Petty et al. (2018)</td><td>2013–2017</td></tr></tbody></table>
Appendix B: List of acronyms

<table><thead><tr><td>Acronym</td><td>Long name</td></tr></thead><tbody><tr><td>ALS</td><td>Airborne laser scanner</td></tr><tr><td>ASD</td><td>Altimetric snow depth</td></tr><tr><td>AMSR</td><td>Advanced Microwave Scanning Radiometer</td></tr><tr><td>ARISE</td><td>Antarctic Remote Ice Sensing Experiment</td></tr><tr><td>ASIRAS</td><td>ESA Airborne Ku-band Interferometric Radar</td></tr><tr><td>ASPeCt</td><td>Antarctic Sea Ice Process and Climate</td></tr><tr><td>ASR</td><td>Arctic System Reanalysis</td></tr><tr><td>AVISO+</td><td>Archiving, Validation and Interpretation of Satellite Oceanographic data</td></tr><tr><td>AWI</td><td>Alfred Wegener Institute</td></tr><tr><td>BROMEX</td><td>Bromine, Ozone and Mercury Experiment</td></tr><tr><td>CMEMS</td><td>Copernicus Marine Environment Monitoring Service</td></tr><tr><td>CNES</td><td>Centre national d'études spatiales</td></tr><tr><td>CPOM</td><td>Centre for Polar Observation and Modelling</td></tr><tr><td>CRISTAL</td><td>Copernicus Polar Ice and Snow Topography Altimeter</td></tr><tr><td>CRREL</td><td>Cold Regions Research and Engineering Laboratory</td></tr><tr><td>CryoSeaNICE</td><td>CryoSat SciEnce-oriented data ANalysis over Sea-ICE areas</td></tr><tr><td>CryoVEx</td><td>CRYOsat Validation Experiment</td></tr><tr><td>CS-2</td><td>CryoSat-2</td></tr><tr><td>CSAO+</td><td>CryoSat+ Antarctic Ocean</td></tr><tr><td>DTU</td><td>Technical University of Denmark</td></tr><tr><td>DuST</td><td>Dual-altimeter Snow Thickness</td></tr><tr><td>EASE</td><td>Equal-Area Scalable Earth</td></tr><tr><td>ECMWF</td><td>European Centre for Medium-Range Weather Forecasts</td></tr><tr><td>ERA-Interim</td><td>Re-Analysis Interim</td></tr><tr><td>ESA</td><td>European Space Agency</td></tr><tr><td>EUMETSAT</td><td>European Organisation for the Exploitation of Meteorological Satellites</td></tr><tr><td>FMCW</td><td>Frequency modulated continuous wave</td></tr><tr><td>FYI</td><td>First year ice</td></tr><tr><td>GCOM-W</td><td>Global Change Observation Mission – Water</td></tr><tr><td>GOP</td><td>Geophysical ocean product</td></tr><tr><td>HPCM</td><td>High Priority Copernicus Candidate Mission</td></tr><tr><td>ICESat-2</td><td>Ice, Cloud, and Land Elevation Satellite-2</td></tr><tr><td>JRA</td><td>Japanese Meteorological Agency</td></tr><tr><td>JAXA</td><td>Japan Aerospace Exploration Agency</td></tr><tr><td>LEGOS</td><td>Laboratoire d'Etudes en Géophysique et Océanographie Spatiales</td></tr><tr><td>LIM</td><td>Louvain-la-Neuve Sea Ice Model</td></tr><tr><td>LRM</td><td>Low-resolution mode</td></tr><tr><td>MERRA</td><td>Modern-Era Retrospective analysis for Research and Applications</td></tr><tr><td>MYI</td><td>Multi-year ice</td></tr><tr><td>N-ICE</td><td>Nansen Basin during the Norwegian young sea ice</td></tr><tr><td>NASA</td><td>National Aeronautics and Space Administration</td></tr><tr><td>MOSAIC</td><td>Multidisciplinary drifting Observatory for the Study of Arctic Climate</td></tr><tr><td>NCAR</td><td>National Center for Atmospheric Research</td></tr><tr><td>NCEP</td><td>National Center for Environmental Prediction</td></tr><tr><td>NEMO</td><td>Nucleus for European Modelling of the Ocean</td></tr><tr><td>NSIDC</td><td>National Snow and Ice Data Center</td></tr><tr><td>SMMR</td><td>Scanning Multichannel Microwave Radiometer</td></tr><tr><td>SMOS</td><td>Soil Moisture Ocean Satellite</td></tr><tr><td>SLA</td><td>Sea level anomaly</td></tr></tbody></table>
<table><tr><td>SSM/I</td><td>Special Sensor Microwave/Imager</td></tr><tr><td>SPICES</td><td>Space-borne Observation for Detecting and Forecasting Sea Ice Cover Extremes project</td></tr><tr><td>IFS</td><td>Integrated Forecast System</td></tr><tr><td>IMB</td><td>Ice mass balance</td></tr><tr><td>ODATIS</td><td>Online Data Extraction Service</td></tr><tr><td>OIB</td><td>Operation Ice Bridge</td></tr><tr><td>OSISAF</td><td>Ocean and Sea Ice Satellite Application Facility</td></tr><tr><td>PLRM</td><td>Pseudo low-resolution mode</td></tr><tr><td>PP</td><td>Pulse peakiness</td></tr><tr><td>SAR</td><td>Synthetic aperture radar</td></tr><tr><td>sd</td><td>Snow depth</td></tr><tr><td>SGDR</td><td>Sensor Geophysical Data Record</td></tr><tr><td>SIT</td><td>Sea ice thickness</td></tr><tr><td>SHEBA</td><td>Surface HEat Budget of the Arctic Ocean</td></tr><tr><td>TFMRA</td><td>Threshold first maximum re-tracker algorithm</td></tr><tr><td>UCL</td><td>University College of London</td></tr><tr><td>W99m</td><td>Warren-99 modified climatology</td></tr></table>

*Data availability.* The ASD data used in this study are freely available at http://ctoh.legos.obs-mip.fr/data/sea-ice-products (CTOH, 2021) In the case of unavailability please contact the authors Florent Garnier (FG) or Sara Fleury (SF) to access the data.

*Supplement.* The supplement related to this article is available online at: https://doi.org/10.5194/tc-15-5483-2021-supplement.

*Author contributions.* FG and SF have conceived the study and wrote the paper. All named authors have participated in the present article within the named projects and have brought contributions to the elaboration of its final version.

*Competing interests.* The authors declare that they have no conflict of interest.

*Disclaimer.* Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Acknowledgements.* This research is supported by the ESA POLAR+ Snow on Sea Ice, the CryoSat+ Antarctic Ocean projects (CSAO+) and the CryoSat SciEnce-oriented data ANalysis over Sea-ICE areas (CryoSeaNICE). It has benefited from the support of the CNES TOSCA CASSIS project and the ESA living planet fellowship. The work has also been supported by the Programme Nationale de Télédétection Spatiale (PNTS, http://programmes.insu.cnrs.fr/pnts/, last access: March 2021), grant no. PNTS-2020-11. We also thank the Center for Topographic studies of the Oceans and Hydrosphere (CTOH) at LEGOS.

*Financial support.* This article has benefited from the financial support of the CNES TOSCA CASSIS, the ESA Living Planet Fellowship, the ESA POLAR+ Snow on Sea Ice project and the Programme Nationale de Télédétection Spatiale (PNTS, http://programmes.insu.cnrs.fr/pnts/), grant nPNTS-2020-11.

*Review statement.* This paper was edited by Stephen Howell and reviewed by three anonymous referees.

## References

Adodo, F. I., Remy, F., and Picard, G.: Seasonal variations of the backscattering coefficient measured by radar altimeters over the Antarctic Ice Sheet, The Cryosphere, 12, 1767–1778, https://doi.org/10.5194/tc-12-1767-2018, 2018.

Alexandrov, V., Sandven, S., Wahlin, J., and Johannessen, O. M.: The relation between sea ice thickness and freeboard in the Arctic, The Cryosphere, 4, 373–380, https://doi.org/10.5194/tc-4-373-2010, 2010.

Andersen, O. B. and Knudsen, P.: The DTU15 Mean Sea Surface and Mean Dynamic Topography, in: Oral Presentation in the 2015 OSTST Meeting, Reston, USA, 2015.

Andersen, O. B., Nilsen, K., Sørensen, L. S., Skourup, H., Andersen, N. H., Nagler, T., Wuite, J., Kouraev, A., Zakharova, E., and Fernandez, D.: Arctic freshwater fluxes from earth observation data, in: Fiducial Reference Measurements for Altimetry, 97–103, Springer, Cham, 2019.

Andreas, E. L., Jordan, R. E., and Makshtas, A. P.: Parameterizing turbulent exchange over sea ice: The Ice Station Weddell results, Bound.-Lay. Meteorol., 114, 439–460, 2005.

Armitage, T. W. K. and Ridout, A. L.: Arctic sea ice freeboard from AltiKa and comparison with CryoSat2 and Operation IceBridge, Geophys. Res. Lett., 42, 6724–6731, https://doi.org/10.1002/2015GL064823, 2015.

Arndt, S. and Nicolaus, M.: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice, The Cryosphere, 8, 2219–2233, https://doi.org/10.5194/tc-8-2219-2014, 2014.
Arndt, S., Meiners, K. M., Ricker, R., Krumpen, T., Katlein, C., and Nicolaus, M.: Influence of snow depth and surface flooding on light transmission through A ntarctic pack ice, J. Geophys. Res.-Oceans, 122, 2108–2119, 2017.

Bin, C., Vihma, T., Zhanhai, Z., Zhijun, L., and Huiding, W.: Snow and sea ice thermodynamics in the Arctic: Model validation and sensitivity study against SHEBA data, Advances in Polar Science, 19, 108–122, 2008.

Blanchard-Wrigglesworth, E., Farrell, S., Newman, T., and Bitz, C.: Snow cover on Arctic sea ice in observations and an Earth System Model, Geophys. Res. Lett., 42, 10–342, 2015.

Blazey, B. A., Holland, M. M., and Hunke, E. C.: Arctic Ocean sea ice snow depth evaluation and bias sensitivity in CCSM, The Cryosphere, 7, 1887–1900, https://doi.org/10.5194/tc-7-1887-2013, 2013.

Boisvert, L. N., Webster, M. A., Petty, A. A., Markus, T., Bromwich, D. H., and Cullather, R. I.: Intercomparison of Precipitation Estimates over the Arctic Ocean and Its Peripheral Seas from Reanalyses, J. Climate, 31, 8441–8462, https://doi.org/10.1175/JCLI-D-18-0125.1, 2018.

Bouffard, J., Naeije, M., Banks, C. J., Calafat, F. M., Cipollini, P., Snaith, H. M., Webb, E., Hall, A., Mannan, R., Féménias, P., and Parrinello, T.: CryoSat ocean product quality status and future evolution, Adv. Space Res., 62, 1549–1563, https://doi.org/10.1016/j.asr.2017.11.043, 2018a.

Bouffard, J., Webb, E., Scagliola, M., Garcia-Mondéjar, A., Baker, S., Brockley, D., Gaudelli, J., Muir, A., Hall, A., Mannan, R., Roca, M., Fornari, M., Féménias, P., and Parrinello, T.: CryoSat instrument performance and ice product quality status, Adv. Space Res., 62, 1526–1548, https://doi.org/10.1016/j.asr.2017.11.024, 2018b.

Braakmann-Folgmann, A. and Donlon, C.: Estimating snow depth on Arctic sea ice using satellite microwave radiometry and a neural network, The Cryosphere, 13, 2421–2438, https://doi.org/10.5194/tc-13-2421-2019, 2019.

Brucker, L. and Markus, T.: Arctic-scale assessment of satellite passive microwave-derived snow depth on sea ice using Operation IceBridge airborne data, J. Geophys. Res.-Oceans, 118, 2892–2905, 2013.

Bunzel, F., Notz, D., Baehr, J., Müller, W. A., and Fröhlich, K.: Seasonal climate forecasts significantly affected by observational uncertainty of Arctic sea ice concentration, Geophys. Res. Lett., 43, 852–859, https://doi.org/10.1002/2015GL066928, 2016.

Chang, A., Foster, J., and Hall, D. K.: Nimbus-7 SMMR derived global snow cover parameters, Ann. Glaciol., 9, 39–44, 1987.

Chevallier, M., Smith, G. C., Dupont, F., Lemieux, J.-F., Forget, G., Fujii, Y., Hernandez, F., Msadek, R., Peterson, K. A., Storto, A., Toyoda, T., Valdivieso, M., Vernieres, G., Zuo, H., Balmaseda, M., Chang, Y.-S., Ferry, N., Garric, G., Haines, K., Keeley, S., Kovach, R. M., Kuragano, T., Masina, S., Tang, Y., Tsujino H., and Wang, X.: Intercomparison of the Arctic sea ice cover in global ocean–sea ice reanalyses from the ORA-IP project, Clim. Dynam., 49, 1107–1136, 2017.

Comiso, J. C., Cavalieri, D. J., and Markus, T.: Sea ice concentration, ice temperature, and snow depth using AMSR-E data, IEEE T. Geosci. Remote Sens., 41, 243–252, 2003.

CTOH: ASD data, available at: http://ctoh.legos.obs-mip.fr/data/sea-ice-products, last access: November 2021.

Déry, S. J. and Tremblay, L.: Modeling the effects of wind redistribution on the snow mass budget of polar sea ice, J. Phys. Oceanogr., 34, 258–271, 2004.

Dong, C., Gao, X., Zhang, Y., Yang, J., Zhang, H., and Chao, Y.: Multiple-scale variations of sea ice and ocean circulation in the Bering Sea using remote sensing observations and numerical modeling, Remote Sens., 11, 1484, 2019.

Eicken, H., Lange, M. A., and Wadhams, P.: Characteristics and distribution patterns of snow and meteoric ice in the Weddell Sea and their contribution to the mass balance of sea ice, Ann. Geophys., 12, 80–93, https://doi.org/10.1007/s00585-994-0080-x, 1994.

Eicken, H., Fischer, H., and Lemke, P.: Effects of the snow cover on Antarctic sea ice and potential modulation of its response to climate change, Ann. Glaciol., 21, 369–376, 1995.

Farrell, S. L., Kurtz, N., Connor, L. N., Elder, B. C., Leuschen, C., Markus, T., McAdoo, D. C., Panzer, B., Richter-Menge, J., and Sonntag, J. G.: A first assessment of IceBridge snow and ice thickness data over Arctic sea ice, IEEE T. Geosci. Remote Sens., 50, 2098–2111, 2011.

Fichefet, T. and Maqueda, M. M.: Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics, J. Geophys. Res.-Oceans, 102, 12609–12646, 1997.

Fiedler, E. K., Martin, M., Blockley, E., Mignac, D., Fournier, N., Ridout, A., Shepherd, A., and Tilling, R.: Assimilation of sea ice thickness derived from CryoSat-2 along-track freeboard measurements into the Met Office’s Forecast Ocean Assimilation Model (FOAM), The Cryosphere Discuss. [preprint], https://doi.org/10.5194/tc-2021-127, in review, 2021.

Fons, S., Kurtz, N. T., Bagnardi, M., Petty, A. A., and Tilling, R.: Assessing CryoSat-2 Antarctic snow freeboard retrievals using data from ICESat-2, Earth and Space Science Open Archive, 23, https://doi.org/10.1002/essoar.10506473.1, 2021.

Fons, S. W. and Kurtz, N. T.: Retrieval of snow freeboard of Antarctic sea ice using waveform fitting of CryoSat-2 returns, The Cryosphere, 13, 861–878, https://doi.org/10.5194/tc-13-861-2019, 2019.

Forsström, S., Gerland, S., and Pedersen, C. A.: Thickness and density of snow-covered sea ice and hydrostatic equilibrium assumption from in situ measurements in Fram Strait, the Barents Sea and the Svalbard coast, Ann. Glaciol., 52, 261–270, 2011.

Giles, K., Laxon, S., Wingham, D., Wallis, D., Krabill, W., Leuschen, C., McAdoo, D., Manizade, S., and Raney, R.: Combined airborne laser and radar altimeter measurements over the Fram Strait in May 2002, Remote Sens. Environ., 111, 182–194, 2007.

Giles, K. A., Laxon, S. W., Ridout, A. L., Wingham, D. J., and Bacon, S.: Western Arctic Ocean freshwater storage increased by wind-driven spin-up of the Beaufort Gyre, Nat. Geosci., 5, 194–197, 2012.

Granskog, M. A., Assmy, P., Gerland, S., Sreen, G., Steen, H., and Smedsrud, L. H.: Arctic research on thin ice: Consequences of Arctic sea ice loss, Eos Trans. AGU, 97, 22–26, 2016.

Grenfell, T. C. and Maykut, G. A.: The optical properties of ice and snow in the Arctic Basin, J. Glaciol., 18, 445–463, 1977.

Grody, N. C.: Classification of snow cover and precipitation using the Special Sensor Microwave Imager, J. Geophys. Res.-Atmos., 96, 7423–7435, 1991.
Grosfeld, K., Treffeisen, R., Asseng, J., Bartsch, A., Bräuer, B., Fritzsch, B., Gerdes, R., Hendricks, S., Hiller, W., Heygster, G., Krumpen, T., Lemke, P., Melsheimer, C., Nicolaus, M., Ricker, R., and Weigelt, M.: Online sea-ice knowledge and data platform, available at: https://www.meereisportal.de/ (last access: March 2021), Polarforschung, 85, 143–155, 2016.

Guerreiro, K., Fleury, S., Zakharova, E., Rémy, F., and Kouraev, A.: Potential for estimation of snow depth on Arctic sea ice from CryoSat-2 and SARAL/AltiKa missions, Remote Sens. Environ., 186, 339–349, 2016.

Guerreiro, K., Fleury, S., Zakharova, E., Kouraev, A., Rémy, F., and Maisongrande, P.: Comparison of CryoSat-2 and ENVISAT radar freeboard over Arctic sea ice: toward an improved Envisat freeboard retrieval, The Cryosphere, 11, 2059–2073, https://doi.org/10.5194/tc-11-2059-2017, 2017.

Haas, C., Haapala, J., Hanson, S., Rabenstein, L., Rinne, E., and Wilkinson, J.: CryoVEx 2006: field report, ESA/ESTEC contract 18677/04/NL/GS, CCN 2, Bremerhaven, Alfred Wegener Institute for Polar and Marine Research, 2006.

Haas, C., Beckers, J., King, J., Silis, A., Stroeve, J., Wilkinson, J., Notenboom, B., Schweiger, A., and Hendricks, S.: Ice and snow thickness variability and change in the high Arctic Ocean observed by in situ measurements, Geophys. Res. Lett., 44, 10–462, 2017.

Helm, V., Hendricks, S., Göbell, S., Rack, W., Haas, C., Nixdorf, U., and Boebel, T.: CryoVex 2004 and 2005 (BoB) data acquisition and final report, Alfred Wegener Institute, Bremerhaven, Germany, 2006.

Helm, V., Humbert, A., and Miller, H.: Elevation and elevation change of Greenland and Antarctica derived from CryoSat-2, The Cryosphere, 8, 1539–1559, https://doi.org/10.5194/tc-8-1539-2014, 2014.

Hendricks, S., Stenseng, L., Helm, V., and Haas, C.: Effects of surface roughness on sea ice freeboard retrieval with an Airborne Ku-Band SAR radar altimeter, in: 2010 IEEE International Geoscience and Remote Sensing Symposium, 3126–3129, IEEE, Honolulu, 2010.

Hendricks, S., Paul, S., and Rinne, E.: Southern hemisphere sea ice thickness from the CryoSat-2 satellite on a monthly grid (L3C), v2.0e thickness and volume, Centre for Environmental Data Analysis [data set], https://doi.org/10.5285/48fc3d1e8ada405c8486ada522dae9e8, 2018.

Holland, D. M., Mysak, L. A., Manak, D. K., and Oberhuber, J. M.: Sensitivity study of a dynamic thermodynamic sea ice model, J. Geophys. Res.-Oceans, 98, 2561–2586, 1993.

Hvidegaard, S. M., Forsberg, R., and Skourup, H.: Sea ice thickness estimates from airborne laser scanning, Sea Ice Thickness: Past, Present and Future, 193–206, 2006.

Ingram, W., Wilson, C., and Mitchell, J.: Modeling climate change: An assessment of sea ice and surface albedo feedbacks, J. Geophys. Res.-Atmos., 94, 8609–8622, 1989.

Janjić, T., Bormann, N., Bocquet, M., Carton, J. A., Cohn, S. E., Dance, S. L., Losa, S. N., Nichols, N. K., Potthast, R., Waller, J. A., and Weston, P.: On the representation error in data assimilation, Q. J. Roy. Meteor. Soc., 144, 1257–1278, https://doi.org/10.1002/qj.3130, 2018.

Kacimi, S. and Kwok, R.: The Antarctic sea ice cover from ICESat-2 and CryoSat-2: freeboard, snow depth, and ice thickness, The

Cryosphere, 14, 4453–4474, https://doi.org/10.5194/tc-14-4453-2020, 2020.

Kaminski, T., Kauker, F., Toudal Pedersen, L., Voßbeck, M., Haak, H., Niederdrenk, L., Hendricks, S., Ricker, R., Karcher, M., Eicken, H., and Gråbak, O.: Arctic Mission Benefit Analysis: impact of sea ice thickness, freeboard, and snow depth products on sea ice forecast performance, The Cryosphere, 12, 2569–2594, https://doi.org/10.5194/tc-12-2569-2018, 2018.

Kelly, R.: The AMSR-E snow depth algorithm: Description and initial results, Journal of the Remote Sensing Society of Japan, 29, 307–317, 2009.

Kern, M., Cullen, R., Berruti, B., Bouffard, J., Casal, T., Drinkwater, M. R., Gabriele, A., Lecuyot, A., Ludwig, M., Midthassel, R., Navas Traver, I., Parrinello, T., Ressler, G., Andersson, E., Martin-Puig, C., Andersen, O., Bartsch, A., Farrell, S., Fleury, S., Gascoin, S., Guillot, A., Humbert, A., Rinne, E., Shepherd, A., van den Broeke, M. R., and Yackel, J.: The Copernicus Polar Ice and Snow Topography Altimeter (CRISTAL) high-priority candidate mission, The Cryosphere, 14, 2235–2251, https://doi.org/10.5194/tc-14-2235-2020, 2020.

Kern, S.: ESA-CCI_Phase2_Standardized_Manual_Visual_Ship-Based-SealceObservations_v02, https://doi.org/10.26050/WDCC/ESACCIPSMVSBSIOV2, WDC climate, DKRZ German climate computing center [data set], 2020.

Kern, S., Khvorostovsky, K., Skourup, H., Rinne, E., Parsakhoo, Z. S., Djepa, V., Wadhams, P., and Sandven, S.: The impact of snow depth, snow density and ice density on sea ice thickness retrieval from satellite radar altimetry: results from the ESA-CCI Sea Ice ECV Project Round Robin Exercise, The Cryosphere, 9, 37–52, https://doi.org/10.5194/tc-9-37-2015, 2015.

King, J., Howell, S., Derksen, C., Rutter, N., Toose, P., Beckers, J. F., Haas, C., Kurtz, N., and Richter-Menge, J.: Evaluation of Operation IceBridge quick-look snow depth estimates on sea ice, Geophys. Res. Lett., 42, 9302–9310, https://doi.org/10.1002/2015GL066389, 2015.

Koenig, L., Martin, S., Studinger, M., and Sonntag, J.: Polar Airborne Observations Fill Gap in Satellite Data, Eos, Transactions American Geophysical Union, 91, 333–334, https://doi.org/10.1029/2010EO380002, 2010.

Kurtz, N., Studinger, M., Harbeck, J., Onana, V., and Farrell, S.: IceBridge sea ice freeboard, snow depth, and thickness, Digital media, NASA Distributed Active Archive Center at the National Snow and Ice Data Center, Boulder, Colorado, USA, available at: http://nsidc.org/data/idcsi2 (last access: March 2021), 2012.

Kurtz, N. T., Farrell, S. L., Studinger, M., Galin, N., Harbeck, J. P., Lindsay, R., Onana, V. D., Panzer, B., and Sonntag, J. G.: Sea ice thickness, freeboard, and snow depth products from Operation IceBridge airborne data, The Cryosphere, 7, 1035–1056, https://doi.org/10.5194/tc-7-1035-2013, 2013.

Kurtz, N. T. and Farrell, S. L.: Large-scale surveys of snow depth on Arctic sea ice from Operation IceBridge, Geophys. Res. Lett., 38, https://doi.org/10.1029/2011GL049216, 2011.

Kwok, R. and Cunningham, G.: Variability of Arctic sea ice thickness and volume from CryoSat-2, Phil. Trans. R. Soc. A, 373, 20140157, https://doi.org/10.1098/rsta.2014.0157, 2015.

Kwok, R., Panzer, B., Leuschen, C., Pang, S., Markus, T., Holt, B., and Gogineni, S.: Airborne surveys of snow
depth over Arctic sea ice, J. Geophys. Res.-Oceans, 116, https://doi.org/10.1029/2011JC007371, 2011.

Kwok, R., Kurtz, N. T., Brucker, L., Ivanoff, A., Newman, T., Farrell, S. L., King, J., Howell, S., Webster, M. A., Paden, J., Leuschen, C., MacGregor, J. A., Richter-Menge, J., Harbeck, J., and Tschudi, M.: Intercomparison of snow depth retrievals over Arctic sea ice from radar data acquired by Operation IceBridge, The Cryosphere, 11, 2571–2593, https://doi.org/10.5194/tc-11-2571-2017, 2017.

Kwok, R., Kacimi, S., Webster, M., Kurtz, N., and Petty, A.: Arctic Snow Depth and Sea Ice Thickness From ICESat-2 and CryoSat-2 Freeboards: A First Examination, J. Geophys. Res.-Oceans, 125, e2019JC016008, 2020.

Laforge, A., Fleury, S., Dinardo, S., Garnier, F., Remy, F., Benveniste, J., Bouffard, J., and Verley, J.: Toward improved sea ice freeboard observation with SAR altimetry using the physical tracker SAMOSA+, Adv. Space Res., 68, 2020.

Landy, J. C., Ehn, J. K., Babb, D. G., Thériault, N., and Barber, D. G.: Sea ice thickness in the Eastern Canadian Arctic: Hudson Bay Complex & Baffin Bay, Remote Sens. Environ., 200, 281–294, 2017.

Landy, J. C., Tsamados, M., and Scharien, R. K.: A facet-based numerical model for simulating SAR altimeter echoes from heterogeneous sea ice surfaces, IEEE T. Geosci. Remote Sens., 57, 4164–4180, 2019.

Landy, J. C., Petty, A. A., Tsamados, M., and Stroeve, J. C.: Sea ice roughness overlooked as a key source of uncertainty in CryoSat-2 ice freeboard retrievals, J. Geophys. Res.-Oceans, 125, e2019JC015820, 2020.

Lawrence, I. R., Tsamados, M. C., Stroeve, J. C., Armitage, T. W. K., and Ridout, A. L.: Estimating snow depth over Arctic sea ice from calibrated dual-frequency radar freeboards, The Cryosphere, 12, 3551–3564, https://doi.org/10.5194/tc-12-3551-2018, 2018.

Laxon, S., Peacock, N., and Smith, D.: High interannual variability of sea ice thickness in the Arctic region, Nature, 425, 947–950, 2003.

Laxon, S. W., Giles, K. A., Ridout, A. L., Wingham, D. J., Willatt, R., Cullen, R., Kwok, R., Schweiger, A., Zhang, J., Haas, C., Hendricks, S., Krishfield, R., Kurtz, N., Farrell, S., and Davidson, M.: CryoSat-2 estimates of Arctic sea ice thickness and volume, Geophys. Res. Lett., 40, 732–737, 2013.

Lecomte, O., Fichefet, T., Vancoppenolle, M., and Nicolaus, M.: A new snow thermodynamic scheme for large-scale sea-ice models, Ann. Glaciol., 52, 337–346, 2011.

Ledley, T. S.: Snow on sea ice: Competing effects in shaping climate, J. Geophys. Res.-Atmos., 96, 17195–17208, 1991.

Lee, Y.-K., Kongoli, C., and Key, J.: An in-depth evaluation of heritage algorithms for snow cover and snow depth using AMSRE and AMSR2 measurements, J. Atmos. Ocean. Technol., 32, 2319–2336, 2015.

Lellouche, J.-M., Greiner, E., Le Galloudec, O., Garric, G., Regnier, C., Drevillon, M., Benkiran, M., Testut, C.-E., Bourdalle-Badie, R., Gasparin, F., Hernandez, O., Levier, B., Drillet, Y., Remy, E., and Le Traon, P.-Y.: Recent updates to the Copernicus Marine Service global ocean monitoring and forecasting real-time 1/12° high-resolution system, Ocean Sci., 14, 1093–1126, https://doi.org/10.5194/os-14-1093-2018, 2018.

Leonard, K. C. and Maksym, T.: The importance of wind-blown snow redistribution to snow accumulation on Bellingshausen Sea ice, Ann. Glaciol., 52, 271–278, 2011.

Liston, G. E., Itkin, P., Stroeve, J., Tschudi, M., Stewart, J. S., Pedersen, S. H., Reinking, A. K., and Elder, K.: A Lagrangian Snow-Evolution System for Sea-Ice Applications (SnowModel-LG): Part I – Model Description, J. Geophys. Res.-Oceans, 125, e2019JC015913, https://doi.org/10.1029/2019JC015913, 2020.

Lundberg, A., Richardson-Näslund, C., and Andersson, C.: Snow density variations: consequences for ground-penetrating radar, Hydrol. Process., 20, 1483–1495, https://doi.org/10.1002/hyp.5944, 2006.

Maaß, N., Kaleschke, L., Tian-Kunze, X., and Drusch, M.: Snow thickness retrieval over thick Arctic sea ice using SMOS satellite data, The Cryosphere, 7, 1971–1989, https://doi.org/10.5194/tc-7-1971-2013, 2013.

Madec, G., Gurvan, M., Bourdallé-Badie, R., Bouttier, P.-A., Bricaud, C., Bruciaferri, D., Calvert, D.,nut, J., Océan, M., Clementi, E., Coward, A., Delrosso, D., Ethé, C., Flavoni, S., Graham, T., Harle, J., Iovino, D., Lea, D., Lévy, C., Lovato, T., Martin, N., Masson, S., Mocavero, S., Paul, J., Rousset, C., Storkey, D., Storto, A., and Vancoppenolle, M.: NEMO ocean engine, Notes du Pôle de modélisation de l'Institut Pierre-Simon Laplace (IPSL): 27, ISSN 1288-1619, https://doi.org/10.5281/zenodo.3248739, 2015.

Maksym, T. and Markus, T.: Antarctic sea ice thickness and snow-to-ice conversion from atmospheric reanalysis and passive microwave snow depth, J. Geophys. Res.-Oceans, 113, https://doi.org/10.1029/2006JC004085, 2008.

Mäkynen, M., Haapala, J., Aulicino, G., Balan-Sarojini, B., Balmaseda, M., Gegiuc, A., Girard-Ardhuin, F., Hendricks, S., Heygster, G., Istomina, L., Kaleschke, L., Karvonen, J., Krumpen, T., Lensu, M., Mayer, M., Parmiggiani, F., Ricker, R., Rinne, E., Schmitt, A., Similä, M., Tietsche, S., Tonboe, R., Wadhams, P., Winstrup, M., and Zuo, H.: Satellite Observations for Detecting and Forecasting Sea-Ice Conditions: A Summary of Advances Made in the SPICES Project by the EU's Horizon 2020 Programme, Remote Sens., 12, 1214, https://doi.org/10.3390/rs12071214, 2020.

Mallett, R. D. C., Lawrence, I. R., Stroeve, J. C., Landy, J. C., and Tsamados, M.: Brief communication: Conventional assumptions involving the speed of radar waves in snow introduce systematic underestimates to sea ice thickness and seasonal growth rate estimates, The Cryosphere, 14, 251–260, https://doi.org/10.5194/tc-14-251-2020, 2020.

Markus, T. and Cavalieri, D. J.: Snow depth distribution over sea ice in the Southern Ocean from satellite passive microwave data, Antarctic sea ice: physical processes, interactions and variability, 19–39, https://doi.org/10.1029/AR074p0019, 1998.

Massom, R. A., Drinkwater, M. R., and Haas, C.: Winter snow cover on sea ice in the Weddell Sea, J. Geophys. Res.-Oceans, 102, 1101–1117, https://doi.org/10.1029/96JC02992, 1997.

Massom, R. A., Eicken, H., Hass, C., Jeffries, M. O., Drinkwater, M. R., Sturm, M., Worby, A. P., Wu, X., Lytle, V. I., Ushio, S., Morris, K., Reid, P. A., Warren, S. G., and Allison, I.: Snow on Antarctic sea ice, Rev. Geophys., 39, 413–445, 2001.

Massom, R. A., Worby, A., Lytle, V., Markus, T., Allison, I., Scambos, T., Enomoto, H., Tamura, T., Tateyama, K., Haran, T., Comiso, J. C., Pfaffling, A., Muto, A., Kanagaratnam, P., Giles,
B., Young, N., Hyland, G., and Key, E.: ARISE (Antarctic Remote Ice Sensing Experiment) in the East 2003: Validation of satellite-derived sea-ice data products, Ann. Glaciol., 44, 288–296, 2006.

Massonnet, F., Fichefet, T., and Goosse, H.: Prospects for improved seasonal Arctic sea ice predictions from multivariate data assimilation, Ocean Modell., 88, 16–25, 2015.

Meier, W., Markus, T., and Comiso, J.: AMSR-E/AMSR2 unified L3 Daily 12.5 km Brightness Temperatures, Sea Ice Concentration, Motion and Snow Depth Polar Grids, Version 1., NASA National Snow and Ice Data Center Distributed archive Center [data set], https://doi.org/10.5067/RA1MIJOYPK3P, 2018.

Merkouriadi, I., Cheng, B., Graham, R. M., Rösel, A., and Granskog, M. A.: Critical Role of Snow on Sea Ice Growth in the Atlantic Sector of the Arctic Ocean, Geophys. Res. Lett., 44, 10,479–10,485, https://doi.org/10.1002/2017GL075494, 2017.

Munoz-Martin, J. F., Perez, A., Camps, A., Ribó, S., Cardellach, E., Stroeve, J., Nandan, V., Itkin, P., Tonboe, R., Hendricks, S., Huntemann, M., Sreen, G., and Pastena, M.: Snow and Ice Thickness Retrievals Using GNSS-R: Preliminary Results of the MOSAiC Experiment, Remote Sens., 12, 4038, 2020.

Nandan, V., Scharien, R. K., Geldsetzer, T., Kwok, R., Yackel, J. J., Mahmud, M. S., Rösel, A., Tonboe, R., Granskog, M., Willatt, R., Stroeve, J., Nomura, D., and Frey, M.: Snow Property Controls on Modeled Ku-Band Altimeter Estimates of First-Year Sea Ice Thickness: Case Studies From the Canadian and Norwegian Arctic, IEEE J. Sel. Top. Appl., 13, 1082–1096, https://doi.org/10.1109/JSTARS.2020.2966432, 2020.

Newman, T., Farrell, S. L., Richter-Menge, J., Connor, L. N., Kurtz, N. T., Elder, B. C., and McAdoo, D.: Assessment of radar-derived snow depth over Arctic sea ice, J. Geophys. Res.-Oceans, 119, 8578–8602, 2014.

Nghiem, S. V., Clemente-Colón, P., Douglas, T., Moore, C., Obrist, D., Perovich, D. K., Pratt, K. A., Rigor, I. G., Simpson, W., Shepson, P. B., Steffen, A., and Woods, J.: Studying bromine, ozone, and mercury chemistry in the Arctic, Eos, Transactions American Geophysical Union, 94, 289–291, 2013.

Notz, D.: Challenges in simulating sea ice in Earth System Models, WIRES Clim. Change, 3, 509–526, 2012.

Overland, J., Dunlea, E., Box, J. E., Corell, R., Forsius, M., Kattsov, V., Olsen, M. S., Pawlak, J., Reiersen, L.-O., and Wang, M.: The urgency of Arctic change, Polar Sci., 21, 6–13, 2019.

Parrinello, T., Shepherd, A., Bouffard, J., Badessi, S., Casal, T., Davidson, M., Fornari, M., Maestroni, E., and Scagliola, M.: CryoSat: ESA's ice mission – Eight years in space, Adv. Space Res., 62, 1178–1190, https://doi.org/10.1016/j.asr.2018.04.014, 2018.

Paul, S., Hendricks, S., Ricker, R., Kern, S., and Rinne, E.: Empirical parametrization of Envisat freeboard retrieval of Arctic and Antarctic sea ice based on CryoSat-2: progress in the ESA Climate Change Initiative, The Cryosphere, 12, 2437–2460, https://doi.org/10.5194/tc-12-2437-2018, 2018.

Perovich, D., Richter-Menge, J., and Polashenski, C.: Observing and understanding climate change: Monitoring the mass balance, motion, and thickness of Arctic sea ice, available at: http://imb-crrel-dartmouth.org (last access: March 2021), 2021.

Perovich, D. K.: Light reflection and transmission by a temperate snow cover, J. Glaciol., 53, 201–210, 2007.

Perovich, D.-K. and Richter-Menge: Regional variability in sea ice melt in a changing Arctic, Mathematical and physical and engineering sciences, 373, 20140165, https://doi.org/10.1098/rsta.2014.0165, 2015.

Perovich, D. K. and Richter-Menge, J. A.: Surface characteristics of lead ice, J. Geophys. Res.-Oceans, 99, 16341–16350, 1994.

Perovich, D. K., Andreas, E. L., Curry, J. A., Eiken, H., Fairall, C. W., Grenfell, T. C., Guest, P., Intrieri, J., Kadko, D., Lindsay, R. W., McPhee, M. G., Morison, J., Moritz, R. E., Paulson, C. A., Pegau, W. S., Persson, P., Pinkel, R., Richter-Menge, J. A., Stanton, T., Stern, H., Sturm, M., Tucker III, W., and Uttal, T.: Year on ice gives climate insights, Eos, Transactions American Geophysical Union, 80, 481–486, https://doi.org/10.1029/EO080i041p00481-01, 1999.

Perovich, D. K., Grenfell, T. C., Richter-Menge, J. A., Light, B., Tucker III, W. B., and Eicken, H.: Thin and thinner: Sea ice mass balance measurements during SHEBA, J. Geophys. Res.-Oceans, 108, https://doi.org/10.1029/2001JC001079, 2003.

Petty, A. A., Webster, M., Boisvert, L., and Markus, T.: The NASA Eulerian Snow on Sea Ice Model (NESOSIM) v1.0: initial model development and analysis, Geosci. Model Dev., 11, 4577–4602, https://doi.org/10.5194/gmd-11-4577-2018, 2018.

Powell, D. C., Markus, T., and Stössel, A.: Effects of snow depth forcing on Southern Ocean sea ice simulations, J. Geophys. Res.-Oceans, 110, https://doi.org/10.1029/2003JC002212, 2005.

Richter-Menge, J. A., Perovich, D. K., Elder, B. C., Claffey, K., Rigor, I., and Ortmeyer, M.: Ice mass-balance buoys: a tool for measuring and attributing changes in the thickness of the Arctic sea-ice cover, Ann. Glaciol., 44, 205–210, 2006.

Ricker, R., Hendricks, S., Helm, V., Skourup, H., and Davidson, M.: Sensitivity of CryoSat-2 Arctic sea-ice freeboard and thickness on radar-waveform interpretation, The Cryosphere, 8, 1607–1622, https://doi.org/10.5194/tc-8-1607-2014, 2014.

Rostosky, P., Sreen, G., Farrell, S. L., Frost, T., Heygster, G., and Melsheimer, C.: Snow Depth Retrieval on Arctic Sea Ice From Passive Microwave Radiometers – Improvements and Extensions to Multiyear Ice Using Lower Frequencies, J. Geophys. Res.-Oceans, 123, 7120–7138, 2018.

Schweiger, A., Lindsay, R., Zhang, J., Steele, M., Stern, H., and Kwok, R.: Uncertainty in modeled Arctic sea ice volume, J. Geophys. Res.-Oceans, 116, https://doi.org/10.1029/2011JC007084, 2011.

Semenov, A., Zhang, X., Rinke, A., Dorn, W., and Dethloff, K.: Arctic intense summer storms and their impacts on sea ice – A regional climate modeling study, Atmosphere, 10, 218, https://doi.org/10.3390/atmos10040218, 2019.

Serreze, M., Walsh, J., Chapin, F. S., Osterkamp, T., Dyurgerov, M., Romanovsky, V., Oechel, W., Morison, J., Zhang, T., and Barry, R.: Observational evidence of recent change in the northern high-latitude environment, Climatic Change, 46, 159–207, 2000.

Shalina, E. V. and Sandven, S.: Snow depth on Arctic sea ice from historical in situ data, The Cryosphere, 12, 1867–1886, https://doi.org/10.5194/tc-12-1867-2018, 2018.

Shupe, M., Rex, M., Dethloff, K., Damm, E., Fong, A., Gradinger, R., Heuze, C., Loose, B., Makarov, A., Maslowski, W., Nicolaus, M., Perovich, D., Rabe, B., Rinke, A., Sokolov, V., and Sommerfeld, A: The MOSAiC Expedition: A Year Drifting with the Arctic Sea Ice, Arctic report card, NFS Public access repository, https://doi.org/10.25923/9g3v-xh92, 2020.
Singarayer, J. S., Bamber, J. L., and Valdes, P. J.: Twenty-first-century climate impacts from a declining Arctic sea ice cover, J. Climate, 19, 1109–1125, 2006.

Skourup, H., Hvidegaard, S. M., Forsberg, R., Einarsson, I., Olesen, A. V., Sornsen, L. S., Stenseng, L., Hendricks, S., Helm, V., and Davidson, M.: CryoVEX 2011-12 Airborne Campaigns for CryoSat Validation, 20 Years of Progress in Radar Altimetry, edited by: Ouwehand, L., 710, 98, 24–29 September 2012, Venice, Italy, ISBN 978-92-9221-274-2, 2013.

Stewart, L. M., Dance, S. L., and Nichols, N. K.: Correlated observation errors in data assimilation, Int. J. Numer. Meth. Fl., 56, 1521–1527, https://doi.org/10.1002/fld.1636, 2008.

Stroeve, J., Liston, G. E., Buzzard, S., Zhou, L., Mallett, R., Barrett, A., Tschudi, M., Tsamados, M., Itkin, P., and Stewart, J. S.: A Lagrangian Snow Evolution System for Sea Ice Applications (SnowModel-LG): Part II – Analyses, J. Geophys. Res.-Oceans, 125, e2019JC015900, https://doi.org/10.1029/2019JC015900, 2020.

Sturm, M. and Massom, R.: Snow in the sea ice system: Friend or foe?, 65–109, Wiley publisher, in: Sea ice, https://doi.org/10.1002/9781118778371.ch3, 2016.

Sturm, M. and Massom, R. A.: Snow and sea ice, Sea ice, 2, 153–204, 2009.

Sturm, M., Holmgren, J., König, M., and Morris, K.: The thermal conductivity of seasonal snow, J. Glaciol., 43, 26–41, 1997.

Sturm, M., Holmgren, J., and Perovich, D. K.: Winter snow cover on the sea ice of the Arctic Ocean at the Surface Heat Budget of the Arctic Ocean (SHEBA): Temporal evolution and spatial variability, J. Geophys. Res.-Oceans, 107, SHE 23–1–SHE 23–17, https://doi.org/10.1029/2000JC000400, 2002.

Sturm, M., Maslanik, J. A., Perovich, D., Stroeve, J. C., Richter-Menge, J., Markus, T., Holmgren, J., Heinrichs, J. F., and Tape, K.: Snow depth and ice thickness measurements from the Beaufort and Chukchi Seas collected during the AMSR-Ice03 campaign, IEEE T. Geosci. Remote, 44, 3009–3020, 2006.

Ulaby, F., Moore, R. K., and Fung, A. K.: Microwave remote sensing: Active and passive. Volume 3 – From theory to applications, Addison-Wesley, Theory to Applications, 997 pp., 1986.

Uotila, P., Goosse, H., Haines, K., Chevallier, M., Barthélemy, A., Bricaud, C., Carton, J., Fučkar, N., Garric, G., Iovino, D., Kauker, F., Korhonen, M., Lien, V. S., Marnela, M., Massonnet, F., Mignac, D., Peterson, K. A., Sadikni, R., Shi, L., Tietsche, S., Toyoda, T., Xie, J., and Zhang, Z.: An assessment of ten ocean reanalyses in the polar regions, Clim. Dynam., 52, 1613–1650, 2019.

Van Leeuwe, M. A., Tedesco, L., Arrigo, K. R., Assmy, P., Campbell, K., Meiners, K. M., Rintala, J.-M., Selz, V., Thomas, D. N., and Stefels, J.: Microalgal community structure and primary production in Arctic and Antarctic sea ice: A synthesis, Elementa: Science of the Anthropocene, 6, https://doi.org/10.1525/elementa.267, 2018.

Vancoppenolle, M., Bouillon, S., Fichefet, T., Goosse, H., Lecomte, O., Morales Maqueda, M., and Madec, G.: The Louvain-la-Neuve sea ice model, Notes du pole de modélisation, Institut Pierre-Simon Laplace (IPSL), Paris, France, 2012.

Warren, S. G., Rigor, I. G., Untersteiner, N., Radionov, V. F., Bryazgin, N. N., Aleksandrov, Y. I., and Colony, R.: Snow Depth on Arctic Sea Ice, J. Climate, 12, 1814–1829, https://doi.org/10.1175/1520-0442(1999)012<1814:SDOASI>2.0.CO;2, 1999.

Webster, M. A., Rigor, I. G., Nghiem, S. V., Kurtz, N. T., Farrell, S. L., Perovich, D. K., and Sturm, M.: Interdecadal changes in snow depth on Arctic sea ice, J. Geophys. Res.-Oceans, 119, 5395–5406, https://doi.org/10.1002/2014JC009985, 2014.

Willatt, R., Laxon, S., Giles, K., Cullen, R., Haas, C., and Helm, V.: Ku-band radar penetration into snow cover on Arctic sea ice using airborne data, Ann. Glaciol., 52, 197–205, 2011.

Wingham, D., Francis, C., Baker, S., Bouzinac, C., Brockley, D., Cullen, R., de Chateau-Thierry, P., Laxon, S., Mallow, U., Mavrocordatos, C., Phalippou, L., Ratier, G., Rey, L., Rostand, F., Viau, P., and Wallis, D. W.: CryoSat: A mission to determine the fluctuations in Earth’s land and marine ice fields, Adv. Space Res., 37, 841–871, 2006.

Worby, A. P., Geiger, C. A., Paget, M. J., Van Woert, M. L., Ackley, S. F., and DeLiberty, T. L.: Thickness distribution of Antarctic sea ice, J. Geophys. Res.-Oceans, 113, 2008a.

Worby, A. P., Markus, T., Steer, A. D., Lytle, V. I., and Massom, R. A.: Evaluation of AMSR-E snow depth product over East Antarctic sea ice using in situ measurements and aerial photography, J. Geophys. Res.-Oceans, 113, https://doi.org/10.1029/2007JC004181, 2008b.

Worby, A. P., Steer, A., Lieser, J. L., Heil, P., Yi, D., Markus, T., Allison, I., Massom, R. A., Galin, N., and Zwally, J.: Regional-scale sea-ice and snow thickness distributions from in situ and satellite measurements over East Antarctica during SIPEX 2007, Deep-Sea. Res. Pt. II, 58, 1125–1136, 2011.

Zhang, J. and Rothrock, D. A.: Modeling Global Sea Ice with a Thickness and Enthalpy Distribution Model in Generalized Curvilinear Coordinates, Mon. Weather Rev., 131, 845–861, https://doi.org/10.1175/1520-0493(2003)131<0845:MGSIWA>2.0.CO;2, 2003.

Zhou, L., Xu, S., Liu, J., and Wang, B.: On the retrieval of sea ice thickness and snow depth using concurrent laser altimetry and L-band remote sensing data, The Cryosphere, 12, 993–1012, https://doi.org/10.5194/tc-12-993-2018, 2018.

Zhou, L., Stroeve, J., Xu, S., Petty, A., Tilling, R., Winstrup, M., Rostosky, P., Lawrence, I. R., Liston, G. E., Ridout, A., Tsamados, M., and Nandan, V.: Inter-comparison of snow depth over Arctic sea ice from reanalysis reconstructions and satellite retrieval, The Cryosphere, 15, 345–367, https://doi.org/10.5194/tc-15-345-2021, 2021.

Zygmuntowska, M., Rampal, P., Ivanova, N., and Smedsrud, L. H.: Uncertainties in Arctic sea ice thickness and volume: new estimates and implications for trends, The Cryosphere, 8, 705–720, https://doi.org/10.5194/tc-8-705-2014, 2014.