This document describes the theoretical basis of the algorithms used to determine properties of snow and ice from the measurements of the Ocean and Land Color Instrument (OLCI) onboard Sentinel-3 satellites within the Pre-operational Sentinel-3 snow and ice products (SICE) project: http://snow.geus.dk/. The code used for the SICE retrieval and its documentation can be found at https://github.com/GEUS-SICE/pySICE. The algorithms were developed after the work from Kokhanovsky et al. (2018, 2019, 2020).

We present the update of the Snow and Ice (SICE) property retrieval algorithm proposed initially by Kokhanovsky et al. (2019). The algorithm is based on the spectral measurements of Ocean and Land Color Instrument (OLCI) onboard Sentinel-3 satellites combined with the asymptotic radiative transfer theory valid for weakly absorbing turbid media. The main improvements include the introduction of a new atmospheric correction, retrieval of snow impurity load and properties, retrievals for partially snow-covered ground and also accounting for various thresholds to be used to assess the retrieval quality. The algorithm is available as python and Fortran packages at https://github.com/GEUS-SICE/pySICE. The technique can be applied to various optical sensors (satellite and ground-based) operated in the visible and near infrared regions of electromagnetic spectra.

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing (http://step.esa.int/main/toolboxes/snap/). In this work we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400-1020nm, we derive important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on the spatial grid of 300m. The algorithm also incorporates cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long – term snow property records essential for climate monitoring and data assimilation studies - especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.