the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-avatar-thumb150.png)
Global sensitivity analysis of simulated remote sensing polarimetric observations over snow
Matteo Ottaviani
Gabriel Harris Myers
Nan Chen
This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors.
The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum.
The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance.
The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration.
The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.
- Article
(7631 KB) - Full-text XML
- BibTeX
- EndNote
The quantification of the surface energy balance of snow-covered regions is of extreme importance for Earth-system model simulations aimed at global climate studies (Hansen and Nazarenko, 2004; Fettweis et al., 2008; van den Broeke et al., 2011; Rae et al., 2012; van Angelen et al., 2012; Tedesco et al., 2013; Alexander et al., 2014; Colgan et al., 2014). Since snow albedo fundamentally depends on the optical and microphysical properties of ice crystals (Wiscombe and Warren, 1980; Aoki et al., 2000; Flanner and Zender, 2006; Bougamont et al., 2007; Dang et al., 2016; He et al., 2018) and light-absorbing impurities potentially present in the snowpack (Warren and Wiscombe, 1980; Hansen and Nazarenko, 2004; Dumont et al., 2014), better knowledge of the properties of such components and their evolution is a high-priority objective for the modeling of the cryosphere (Tedesco et al., 2013; Dumont et al., 2014).
One fundamental source of uncertainty in the remote sensing of these properties is the treatment of snow as a collection of spherical grains (Tanikawa et al., 2020). Although useful in some contexts, such approximation can underestimate the albedo by a few percent (Xie et al., 2006; Tedesco and Kokhanovsky, 2007; Libois et al., 2013; Tedesco et al., 2013; Dumont et al., 2014; Dang et al., 2016; Räisänen et al., 2017), a discrepancy that can be exaggerated by snow albedo feedback processes (Thackeray et al., 2018; Colman, 2013; Hansen and Nazarenko, 2004). To avoid this assumption we employ hexagonal prisms, which have been demonstrated (van Diedenhoven et al., 2012) to serve well as radiative proxies for more complex shapes, while having the advantage of being characterized only by their aspect ratio (AR; with AR > 1 for columns and AR < 1 for plates) and the microscale roughness (D) of the crystal facets. The implementation of this methodology in advanced radiative transfer (RT) models has produced successful retrievals for parameters that describe the crystals forming ice clouds (van Diedenhoven et al., 2014b) and the reflectance properties of snow-covered surfaces (Ottaviani et al., 2012, 2015) from data collected with the NASA Goddard Institute for Space Studies (GISS) airborne Research Scanning Polarimeter (RSP; Cairns et al., 1999).
Another major challenge is the determination of light-absorbing particulate (LAP) content and its partitioning between LAPs deposited in snow versus those suspended above in the form of atmospheric aerosols (Warren, 2013). Because the polarization state of light is also sensitive to this partitioning (Ottaviani, 2022), remote sensors like the RSP; the Second-generation Global Imager (SGLI; Tanaka et al., 2018); the polarimeters launched aboard the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission (Werdell et al., 2019; Hasekamp et al., 2018); and the upcoming Multi-viewing, Multi-channel, Multi-polarization imaging mission (3MI) (Biron et al., 2013; Marbach et al., 2013) offer augmented retrieval capabilities.
Zhang et al. (2023) have recently evaluated the performance of a snow kernel introduced into an inverse algorithm in retrieving the microphysics of aerosols above snow, based on observations of the Polarization and Directionality of the Earth's Reflectances (POLDER) spaceborne sensor, which flew from 2004 to 2013. However, their study does not address the retrieval of the microphysical properties of the snowpack and the distribution of LAPs between the snow and the atmosphere or information content aspects. This paper extends the studies presented in Ottaviani (2022), examining these details via a global sensitivity analysis (GSA) of simulated top-of-the-atmosphere (TOA) polarimetric observations.
Section 2 explains the setup of the RT calculations needed to generate the look-up table (LUT) of the Stokes vectors at the TOA produced for random combinations of the input parameters and outlines the adaptation of the GSA formalism to the scopes of this analysis. The results of the GSA are presented in Sect. 3 along with inverse retrievals of spaceborne observations simulated for different observational configurations. We first consider pure-snow scenes to highlight the sensitivity to the ice crystal properties. We then address more realistic remote sensing scenarios where the atmosphere is allowed to contain a layer of absorbing aerosols, and the snowpack contains impurities. The paper concludes with some recommendations for operational retrievals.
2.1 Radiative transfer simulations
The plane-parallel RT code employed to generate the LUT is based on the general doubling–adding formalism described by De Haan et al. (1987). It features a consistent treatment of the radiative effects deriving from atmospheric molecular scattering, aerosols and clouds, and any surface whose reflectance is known in its analytical form or in terms of its bidirectional reflectance distribution function (BRDF) properties and the polarization counterpart (BPDF). The code has been used for decades to model measurements from the RSP over a variety of Earth scenes, including those containing ice crystals in clouds (van Diedenhoven et al., 2013) and ground snow (Ottaviani et al., 2012, 2015).
Recognizing the similarities to the polarimetric signatures of ice crystals in cirrus clouds (Ottaviani et al., 2012, 2015), the snowpack is modeled as an optically semi-infinite collection of non-spherical ice crystals at the bottom of the atmosphere so that the reflectance of the actual underlying surface is irrelevant. As with other media, the photon penetration depth in snow depends on wavelength (Kokhanovsky, 2022; Libois et al., 2013), so different instrumental channels effectively probe different depths, and this fact must be taken into account during multi-spectral retrievals (Li et al., 2001). In contrast to retrievals of grain size for mono-layer snowpacks (Nolin and Dozier, 1993; Painter et al., 2003), such an approach has been exploited to retrieve grain size in both a thin surface layer and a thick layer below using measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS; Aoki et al., 2007; Painter et al., 2009). For this reason, the snowpack in our model is vertically resolved in a thick bottom layer capped by a thin top layer (see Fig. 1).
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f01](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f01-thumb.png)
Figure 1Model for the snow–atmosphere system. The snowpack is described as a thin layer sitting on an optically semi-infinite layer, for a total thickness of 1 m. The parameters in each layer are varied independently. The atmosphere is allowed to contain light-absorbing aerosols () in the lowest layer (first 2 km above the surface). See Table 1 for the complete list of model parameters and their ranges of variability.
To span a wide spectral range as required by remote sensing applications, we consider channels at 411, 469, 555, 670, 864, 1589, and 2266 nm, which in the case of the RSP sensor are all equipped with polarization capabilities. Except for minimal differences in the precise centerband values, the channels in this set are also available from heritage instruments, such as the Moderate Resolution Imaging Spectroradiometer (MODIS), to favor atmospheric correction both over land and over the ocean. However, in many cases these sensors provide only total reflectance and at a single view per pixel. We include the MODIS band at 2112 nm in the list of channels despite it being very close to 2266 nm because these wavelengths lie on the shoulder of a major absorption band (Chen et al., 2017) where radiative differences can arise very quickly. Moreover, this channel will be available from 3MI with polarization capabilities.
Absorbing aerosols are climatologically relevant because deposition events can cause large variations in albedo (Dumont et al., 2014; Hansen and Nazarenko, 2004; Khan et al., 2023; Warren and Wiscombe, 1980). To examine the capability of different observational configurations (see Sect. 3) to distinguish them from impurities in the snow, both the lowest atmospheric layer (located within the first 2 km above the snowpack) and the snow are allowed to contain variable amounts of the same spherical LAP with properties typical of soot (n = 1.80–0.6i, reff = 0.11 µm, veff = 0.38) (Dubovik et al., 2002).
Finally, the presence of exclusively absorbing gases (H2O, O3, etc.) in the background atmosphere is neglected because it does not affect the conclusions drawn from the sensitivity study presented below. The complete list of the descriptive parameters and their bounds in the LUT can be found in Table 1.
Table 1Model parameters and their lower and upper bounds. The parameters describing the top and bottom snow layer are denoted with superscripts T and B and are varied independently. Soot-like LAPs are characterized by n = 1.80–0.6i, reff = 0.11 µm, and veff = 0.38 (Dubovik et al., 2002; Torres et al., 2017).
![](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-t01-thumb.png)
The optical properties of the hexagonal prisms are produced via a geometric optics (GO) code (Macke et al., 1996; van Diedenhoven et al., 2012) as a function of the aspect ratio (AR), microscale roughness (D), and effective radius (reff) and are integrated over a power-law size distribution (Geogdzhayev and van Diedenhoven, 2016). The microscale roughness represents the standard deviation of the distribution of angles used to randomly perturb the orientation of the ice crystal facet encountered by the incident beam in the GO calculations (van Diedenhoven et al., 2014a). Previous attempts to fit the surface contribution to the signal measured by airborne polarimeters have shown that D ≳ 0.25 (Ottaviani, 2012, 2015). Such roughness is sufficient to extinguish the halo peaks characteristic of more pristine crystals (van Diedenhoven et al., 2012).
The asymmetry parameters of columns and plates of reciprocal ARs are very similar (Ottaviani et al., 2015; van Diedenhoven et al., 2014a). To enable retrievals of representative crystal shapes and test the sensitivity to the mixing proportion, we therefore assume that the population of grains is composed of a fraction f of plates with aspect ratio ARP and a fraction (1−f) of columns with aspect ratio ARC=1 ARP. The fraction is allowed to vary independently in each layer. For ice crystals of a given AR, the extinction and scattering efficiencies are
where σext and σsca are the extinction and scattering coefficient, and A is the projected area of the hexagonal prism (van Diedenhoven et al., 2014a). Using again superscripts “P” and “C” for “plates” and “columns”, the extinction and scattering efficiencies of the mixture are
and the corresponding phase function and asymmetry parameter are
The projected area of the mixture is assumed to be that of the column crystals so that the extinction cross section, the scattering cross section, and the single-scattering albedo of the mixture are
The impurities in the snowpack are externally mixed with the snow grains (Tanikawa et al., 2020). Their optical properties are calculated by Mie calculations inherent to the RT code, as is done for the aerosols, and are assumed to follow lognormal aerosol size distributions (Hansen and Travis, 1974). Using the layer-resolved inherent optical properties above, together with the optical depths, the RT code simulates the TOA reflectances (RI, RQ, RU) corresponding to the first three parameters (I, Q, U) of the Stokes vector, describing the linear state of light polarization. Circular polarization is represented by the fourth element (V), which generally has negligible relevance to remote sensing applications (Kawata, 1978), so it is omitted from the analysis. The modeled reflectances can be output for any viewing geometry. We choose observations along the principal plane for which the collection of scattering angles is maximized.
The GSA considers the total reflectance (RI); the polarized reflectance (); and also the degree of linear polarization (DoLP = ), which in RSP-like instruments is measured at a much higher accuracy than Rp (Knobelspiesse et al., 2012; Cairns et al., 1999).
2.2 Global sensitivity analysis formalism
To illustrate the pitfalls of conventional sensitivity studies applied to hyperdimensional state spaces, Fig. 2 shows the sensitivity of RI, Rp, and DoLP to ART in the SWIR (see also Ottaviani, 2022), for two different values of . For 1280 µm (solid lines) the signals are essentially unaffected by variations in ART; conversely, detectable differences arise for 56 µm (dashed lines). Similarly, the DoLP grows with , but the curves for the two ARs are distinguishable only for small radii. Because the sensitivity to ART depends on and since the sensitivity of any parameter can similarly depend on the set of all other values kept fixed, it is necessary to use a method like the GSA to properly quantify the information content of the model.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f02](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f02-thumb.png)
Figure 2Sensitivity of RI, Rp, and DoLP (columns) to ART, for 1280 µm (solid lines) and 56 µm (dashed lines). The remaining parameters are DT = 0.35, ρT = 0.1 g cm−3, and ZT = 3 cm (top layer), as well as = 320 µm, ARB = 1.0, DB = 0.35, and ρB = 0.3 g cm−3 (bottom layer). The three rows are for the three SWIR wavelengths (1589, 2112, and 2266 nm).
The GSA framework relies on the computation of the so-called Sobol indices (Sobol, 1990). For a model function g(X) of the n state variables (in our case, the parameters in Table 1), which is square-integrable over a parameter space Kn, there exists a functional decomposition in terms of a Haar wavelet basis given by
The Haar wavelets form a basis for the space of all square-integrable functions (L2). Because the simulated reflectances are bounded and smooth over the sample space and the sample space is compact, the model functions for our simulations are square-integrable (i.e., the decomposition in Eq. 10 exists). Squaring both sides of the equation, integrating over the whole parameter space, and using the orthogonality properties of the basis, one obtains
If each variable Xi is uniformly distributed over the parameter space Kn, the left-hand side exactly defines the variance of the model function g(X):
whereas the integrals
are the variances of the functions . Combining the decomposition in Eq. (11) with Eqs. (12) and (13) gives the decomposition of the total variance:
Each Vi term in the first sum corresponds to the “main-effect” contribution of the variable Xi to the model output. The Vi,j terms quantify the pairwise interactions between Xi and Xj, the triplet-wise interactions among Xi, Xj, and Xk, and so on (Saltelli et al., 2008). Dividing both sides of Eq. (14) by the total variance, one obtains
where the “total-effect” Sobol index for the parameter Xj,
quantifies the complete contribution of Xj to the total variance over the entire parameter space, both directly and through interactions among parameters. The Sobol indices are calculated at each angle with the Python software package SALib (Herman and Usher, 2017; Iwanaga et al., 2022), which uses the quasi-Monte Carlo estimators presented in Saltelli et al. (2010). To compare the indices at different wavelengths, they are converted into absolute quantities by multiplying by the total variance at that angle:
A visual interpretation of the Sobol indices is given in Fig. 3, where the gray curves show the total reflectance at 555 nm output by the radiative transfer model for 20 random combinations of input parameters, computed for a solar zenith angle (SZA) of 65°. The region shaded in blue represents the total variance of the model curves and the area in red the absolute total-effect Sobol index relative to . A comparison between and reveals that a significant portion of the total variance of the model is due to variations in .
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f03](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f03-thumb.png)
Figure 3Total reflectance RI (gray lines) at 555 nm simulated for 20 random combinations of the state parameters, as a function of the viewing zenith angle along the principal plane (SZA = 65°). The total standard deviation for this ensemble of curves is shown by the blue-shaded region, and the red-shaded region corresponds to the absolute total-effect Sobol index associated with variations in . The fact that covers a large portion of indicates that a large amount of the variation in RI is due to .
The sensitivities are evaluated against the 1σ uncertainty corresponding to the square root of the diagonal elements of the measurement error covariance matrix (Knobelspiesse et al., 2012). Full measurement covariance matrices, which include off-diagonal elements expressing cross-correlation effects, can be hard to assess (Gao et al., 2023). In line with other studies that consider the uncertainty in multi-angle polarimeters (Hasekamp, 2010; Lebsock et al., 2007; Stamnes et al., 2018; Ottaviani et al., 2012), we also neglect such elements that in any case do not affect our main conclusions.
The light gray areas in Fig. 4 (see also discussion in the next section) correspond to a 3 % (0.5 %) radiometric (polarimetric) accuracy, nominally achievable by modern spaceborne sensors like the polarimeters launched on board PACE. The dark gray areas correspond instead to the higher accuracy of RSP-like sensors (1.5 % radiometric and 0.2 % polarimetric). Note that the uncertainty in Rp is rather different in the two cases, since the error model used for the areas in light gray includes a term proportional to . This term is a result of a filter-wheel-type design, where the measurements needed to compile the Stokes vector of any given scene are not acquired simultaneously (Dubovik et al., 2019; Knobelspiesse et al., 2012), and is large over bright surfaces, as is the case for snow in the VIS–NIR.
If is less than the threshold at all angles, the associated parameter is ruled out as a meaningful contributor to model variance and excluded from the plot. The SALib package also gives relative confidence intervals for each Sobol index, which are converted to absolute confidence intervals around . The chosen sample size is m = 216 (specified as a power of 2 as required by the software package), for a total of (n+2)m = 1 048 576 runs, where n=14 is the number of parameters (Saltelli et al., 2008). This sample size ensures that for all σT curves which lie above the uncertainty threshold, the confidence interval also lies above the uncertainty threshold.
As discussed below, values of above the uncertainty threshold do not necessarily guarantee retrievability, which is impacted by model uncertainty and other sources of error unaccounted for in the covariance matrices used here. In this respect, the retrievals can be considered to be a best-case scenario (Rodgers, 2000; Knobelspiesse et al., 2012).
3.1 Pure snow
As a first example of the application of the GSA, in Fig. 4 we show the angular variations in the absolute total Sobol indices computed for a pure-snow scene. Although such an idealized case is not commonly encountered in real-world scenarios, this exercise is useful for isolating and understanding the sensitivity to the parameters describing the ice crystals; the effects of LAPs are then discussed in the next section. In both cases, we consider optically semi-infinite snowpacks since the focus of this paper is on the remote sensing of snow; heterogeneous pixels constitute an added layer of complexity and will be the subject of future studies.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f04](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f04-thumb.png)
Figure 4Absolute total Sobol index σT for RI, Rp, and DoLP (different columns) as a function of the viewing zenith angle along the principal plane and for a SZA of 65°, computed for a pure snowpack under a clear atmosphere ( 0 ppmw, 0). The indices with their 90 % confidence intervals are only given for parameters with variance above the instrument uncertainty thresholds.
The results are for principal-plane observations and SZA = 65°. Note that the top x axis in the first row reports the conversion to scattering angles, which is useful for connecting scanning viewing zenith angles (VZAs) to forward-scattering and backscattering directions in the discussions that follow. The different wavelengths are in different rows, and only state parameters exhibiting sensitivity over the 1σ instrument uncertainty threshold are included. The GSA correctly captures the known dependence of the reflectance on the top-layer grain size in the SWIR (increasing with viewing zenith angle in the forward-scattering half-plane; see also Fig. 2) due to the large absorption by ice at these wavelengths (Wiscombe and Warren, 1980; Dang et al., 2016), which is exploited in retrieval schemes (Stamnes et al., 2007; Painter et al., 2009, 2003; Nolin and Dozier, 1993). The lack of sensitivity of RI to grain size in the VIS is explained by the fact that pure snow is highly reflective (i.e. non-absorbing) at visible wavelengths, which therefore are normally not used for size retrievals.
The exact penetration depth of each wavelength in the snowpack depends on grain size, snow density, and impurity content. For pure snow of typical densities and grain sizes, the VIS wavelengths penetrate deeper (10–25 cm) than NIR (3–20 cm) and SWIR (0–2 cm) wavelengths (Kokhanovsky, 2022; Libois et al., 2013; Li et al., 2001). Multiple scattering quickly randomizes polarization as the incident light penetrates into a dense medium (van Diedenhoven et al., 2013), explaining why no detectable sensitivity to any of the parameters is found for Rp and DoLP in the VIS within the considered uncertainties. The situation is different in the SWIR, again because of the strong ice absorption. The very limited penetration depth allows the polarization signatures determined by the single-scattering properties of the ice crystals at the very top of the snowpack to be preserved, especially for observations of the DoLP, which are typically achieved with higher measurement accuracy. Besides the evident sensitivity to and DT in this wavelength regime, an interesting result concerns the 2266 nm RSP channel, which seems to access detectable sensitivity to ART not present for the nearby MODIS channel at 2112 nm. Furthermore, the reduced uncertainty in RSP-like sensors reveals sensitivity of the polarized reflectance measured at 864 nm to DT and ART (Ottaviani et al., 2015).
All simulated measurements are insensitive to the mixing proportion (f) of columns and plates in both layers, which can be explained by the similarity of the scattering properties for crystals with reciprocal ARs (see Fig. A3). Because of its large physical thickness, the bottom layer is optically semi-infinite regardless of ρB. Directional changes in the light scattered downward in response to variations in D and AR do not prevent its fast extinction, and the fraction of upwelling photons supplied by the bottom layer stays pretty constant so that RI is insensitive to ρB.
Finally, we note that the LUT includes random selections for the thickness of the top layer and for the snow density in both layers. Other than the minimal information contained in the DoLP at 2266 nm (blue and cyan curves, nearly overlapped), the Sobol indices reveal that these parameters cannot be independently retrieved. The optical depth of the top layer largely determines the observed signal but is proportional to the product of ρT and ZT, which is invariant when one parameter is divided by the same factor used to multiply the other. The same argument is even more valid for the semi-infinite bottom layer, which prohibits passive optical measurements from accessing information on its thickness or density.
The results of the GSA are particularly useful for informing the choice of parameters to be included in the state vector of inverse retrievals. As an example, we generated synthetic TOA observations with the RT code, including random noise added according to the specifications of different sensors. As explained in Sect. 2, the snowpack consists of a mixture of crystals ( 0.5). Fresher snow (smaller grains) is simulated in the top layer ( 150 µm, ZT= 3 cm, ρT= 0.2 g cm−3, and ART= 0.05 for plates; corresponding 20 for columns; and DT= 0.3 as found by Ottaviani et al., 2015). More compact, larger, and rounder grains are located in the bottom layer ( 250 µm, ρB= 0.30 g cm−3, ARB= 0.15 for plates and 6.67 for columns, DB= 0.40), which is optically semi-infinite (τ ≈ 2000). For reasons explained in the Appendix, the physical thicknesses and densities of both layers were not included in the set of retrievables and, together with the other parameters excluded from the state vector, are constrained to the values used to generate the synthetic observations.
Using the LMFit Python library (Newville et al., 2014), a Levenberg–Marquardt nonlinear least-squares optimal estimation scheme (Levenberg, 1944; Marquardt, 1960) was then implemented to retrieve the input parameters. We first consider the configuration of an RSP-like instrument. In satellite imagery, every pixel is characterized by its own set of viewing zenith and azimuth angles. The RSP is instead a scanner, and we chose the principal plane as a scanning direction because it guarantees that the viewing geometries span the largest range of scattering angles.
To highlight the merits of polarization, we compare retrievals that consider the measurement vectors to be the simulated (i) total reflectance, (ii) total reflectance and polarized reflectance, and (iii) total reflectance and DoLP. These retrievals are repeated considering VIS, VIS–NIR, or VIS–NIR–SWIR wavelengths. The SWIR combination consists of 1589 and 2266 nm. All available viewing angles (150 measurements for RSP, roughly between ±60°) are used for the total reflectance. Figure 4 shows that the DoLP and Rp are largely unaffected by variations in any of the parameters (except ) for angles in the backscattering half-plane. If Rp or DoLP for these angles is included in the retrieval, the minimization algorithm attempts to fit very small changes in Rp or the DoLP, which cannot be distinguished from the noise, and the retrieval quality for all parameters (other than ) degrades. Consequently, we subsampled the measurements of Rp and DoLP to positive viewing zenith angles, which in our reference system correspond roughly to the forward-scattering half-plane.
Figure 5 summarizes the values of the state parameters and their uncertainty obtained from the inversion. The solid lines represent the “true” values used in the forward simulations. In contrast, the dashed lines represent the initial guess for each parameter, randomly sampled within the bounds listed in Table 1. The retrievals were repeated a few times to test the stability of the results against the different initial guesses.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f05](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f05-thumb.png)
Figure 5State parameters (different panels) retrieved from the inversion of RSP-like TOA observations, generated for a SZA of 65° along the principal plane. The scene consists of a pure snowpack and a clean atmosphere (see text). The inversions are repeated with and without the inclusion of polarization and at different wavelength combinations. The solid and dashed gray lines are the true values and the initial guesses for each parameter.
In this simplistic scenario of a pure snowpack and a clear atmosphere, Fig. 4 indicates that RI, Rp, and DoLP in the VIS are insensitive to all parameters (except for RI to at 670 nm), and retrievals attempted with these data alone are unsuccessful. The addition of NIR measurements of RI and Rp (if the accuracy of the latter matches RSP levels) gives access to information on grain shape and microscale roughness, as confirmed with real data (Ottaviani et al., 2012, 2015).
The retrieval becomes optimal when polarization capabilities in the SWIR are also included. Using RI+Rp or RI+ DoLP leads to similar performances, with the inversions converging to the true values within the error bars for all parameters except for , which is in any case retrieved moderately well.
We next turn our attention to retrievals simulated for instruments with spectral and angular configurations different from those of the RSP. A MODIS-like sensor is mimicked by considering mono-angle measurements of total reflectance in the VIS–NIR–SWIR. The addition of 16-angle polarimetric measurements in the VIS–NIR mimics POLDER. The viewing geometries for the simulations are taken from actual data collected over Greenland for a pixel near Summit Station (72° N, 39° W) on 8 April 2007 and for SZA ≈ 65°, which corresponds to the SZA used to generate Fig. 4. All available angular measurements are considered for the total reflectance in the VIS–NIR, subsampled to the same range as in the RSP-like case.
The results obtained using these two different instrument configurations are compared to those of the RSP in Fig. 6. Retrievals using MODIS-like data can recover the top-layer grain size (Stamnes et al., 2007; Hori et al., 2007; Painter et al., 2009, 2003; Nolin and Dozier, 1993) but fail to resolve all other parameters in a satisfactory manner. Information on grain shape is still accessible to the MODIS NIR channel but for the viewing zenith angle away from nadir, which exemplifies the utility of multi-angular measurements. Observations of Rp contain information on DT in the NIR, although the larger uncertainty assigned to the simulated spaceborne measurements limits the retrieval quality compared to the RSP-like case. Finally, we note that all these retrievals are robust against different choices of the initial guess for each parameter.
3.2 Scenes containing snow impurities and aerosols
In this section, we expand the analysis to include more realistic scenes characterized by the presence of LAPs both in the snowpack and as atmospheric aerosols. We also consider different SZAs.
Light-absorbing impurities in the snow are often found in very significant amounts in North America, China, and the Arctic (Warren, 2019). In Greenland, especially on the plateau, the concentrations are typically much smaller and therefore difficult to detect via remote sensing (Warren, 2013). Since they anyway have a significant impact on snow visible albedo (Warren and Wiscombe, 1980; Dang et al., 2016), their accurate determination is especially important for climate modeling (Antwerpen et al., 2022; Wang et al., 2020; Alexander et al., 2014; Ryan et al., 2019). To target these challenging retrievals, the results of the GSA are computed for maximum LAP loads of 0.4 for and 1 ppmw for and (Fig. 7), where the subscript “BC” indicates the specific type of black-carbon LAP considered in this paper, with fixed microphysical and optical properties. More sporadic events like thick burning plumes or exceptionally dirty snow are addressed in Appendix A, where the same calculations are repeated with extended ranges of (up to 1.2) and and (up to 10 ppmw). At these higher LAP amounts, the sensitivity of polarimetric measurements in the SWIR to aerosol optical depth is even more pronounced (see Fig. A1 and related discussion).
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f07](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f07-thumb.png)
Figure 7Same as Fig. 4 but for variable amounts of LAPs within the snowpack and in the atmosphere above.
The sensitivity to top- and bottom-layer parameters can be interpreted in terms of the penetration depths discussed in Sect. 3.1. The top-layer thickness ranges between 1 and 5 cm in the LUT, explaining why the DoLP is (weakly) sensitive to in the SWIR at 2266 nm but not at 1589 and 2112 nm, for which the ice absorption is slightly larger.
The total reflectance in the VIS–NIR shows the expected sensitivity to LAPs. The added sensitivity to reflects the additional information provided by polarization. The fact that the polarized reflectance and the DoLP measured in the SWIR are insensitive to the concentration of LAPs in the snowpack implies that polarimetry can also be exploited to (i) extend heritage aerosol retrievals performed over other land surfaces to snow surfaces and (ii) inform the vertical partitioning of LAPs between the atmosphere and the surface (Ottaviani, 2022).
To examine the dependence on solar illumination, Fig. 8 was produced for a SZA of 45°. At this smaller SZA, the angles of maximum shift toward larger viewing angles. Additionally, in the VIS–NIR increases for all parameters. These differences are only minor; the lack of appreciable changes with SZA, at least in this range typical of high latitudes, is therefore attractive in the context of remote sensing applications.
In the VIS–NIR, exhibits an essentially flat behavior well above the detection thresholds at all viewing zenith angles for many of the parameters, with shallow maxima at around nadir except for ART and . In the SWIR, and for and peak at the largest viewing zenith angles. The DoLP now includes sensitivity to DT, still occurring in the forward-scattering half-plane but with peaks at smaller angles. Multi-angle polarization measurements can therefore greatly supplement those of total reflectance, especially when retrieving parameters that express marked angular differences in σT.
Figure 9 provides an alternative display of the information contained in Figs. 7 and 8, for a SZA of 65° (top panel) and SZA of 45° (bottom panel). These heatmaps can aid in the choice of appropriate viewing geometries and channel combinations when designing retrieval algorithms and observational strategies. The intensity of each cell's color is proportional to the maximum value of (left columns), (middle columns), and (right columns) across all VZAs, and the number reports the angular location of these maxima. Numbers close to 0 represent nadir-looking directions, and large positive angles correspond to the forward-scattering directions (see top x axis in Figs. 7 and 8). It is evident that measurements in the forward-scattering half-plane are sensitive to the properties of aerosols and the top snow layer, while nadir-looking geometries favor the determination of parameters deeper in the snowpack. It is also clear that the addition of accurate polarimetric measurements in the VIS–SWIR benefits the retrieval of aerosol and surface properties, especially if they are present at multiple angles.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f09](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f09-thumb.png)
Figure 9Heatmaps of the maximum value of , , and across all viewing angles, for each parameter (row) and wavelength (column) combination. The top panel is for a SZA of 65° and the bottom panel for a SZA of 45°. The number in each box is the viewing zenith angle at which the maximum occurs. Only parameters with sensitivity over the uncertainty thresholds are included.
The simulated retrievals in Sect. 3.1 were duplicated with the inclusion of LAPs. The impurity concentration in the snowpack was set to ppmw, typical of the Greenland plateau (Warren, 2019). The atmosphere contains aerosols with 0.10. Figure 10 summarizes the state parameters and their uncertainty for each type of retrieval, at different measurement vector and wavelength combinations.
Retrievals which use RI + DoLP from only VIS and NIR channels are in general not successful and fail to determine the vertical distribution of impurities because they are challenged by the simultaneous sensitivity to multiple parameters. Improvements are observed when measurements in the SWIR are included due to selective sensitivity to and .
Unsurprisingly, the best performance is achieved by including both total reflectance and DoLP in the VIS–NIR–SWIR. The uncertainties decrease by an order of magnitude when compared to measurements of total reflectance only, confirming that polarimetric measurements in the SWIR are valuable for determining the vertical partitioning of LAPs. Initializing the inversion with first guesses close to the true values helps to decrease the uncertainty in and retrieved from the VIS–NIR–SWIR combination of RI and Rp, but the best performance is still obtained using RI + DoLP.
We next repeat the retrievals using MODIS and MODIS + POLDER-like measurements. For the latter, we use the DoLP in the VIS–NIR in place of Rp because the DoLP manifests detectable sensitivity to LAPs. Figure 11 shows that in the MODIS-like case the inversion struggles to retrieve all parameters except , and the considerable uncertainties can severely impact the accuracy of derived albedo estimates. The addition of multi-angle polarimetric data enables a better determination of all parameters, even at lower angular resolution than the RSP's (Wu et al., 2015), since the angular radiative behavior of the different system components is rather smooth. One exception concerns the vertical profile of impurities in the snow: only when assuming a uniform concentration in the two layers is the retrieval successful (not shown).
The information content of polarimetric simulations over snow scenes was evaluated using a global sensitivity analysis (GSA) method, which accounts for the correlated sensitivity to model parameters across the entire parameter space. A comprehensive look-up table (LUT) was created with an advanced vector radiative transfer model, spanning wavelengths from the VIS to the SWIR. The snow–atmosphere system is vertically resolved and accounts for the presence of light-absorbing particulates (LAPs) both embedded in the snow and aloft in the atmosphere above the snowpack. The Sobol indices computed from the LUT are the primary metrics for the GSA and show the expected sensitivity of total reflectance in the VIS–NIR to LAPs and in the SWIR to snow grain size. In contrast to measurements of total reflectance only, polarimetric measurements inform the vertical distribution of LAIs in the system thanks to differential sensitivity present especially in the SWIR. Retrievals of grain shape from polarimetric measurements in the NIR can be improved by the addition of SWIR channels, leading to better estimates of the asymmetry parameter and, in turn, of the albedo in climate models. The angular dependence of the sensitivity (especially of the DoLP) emphasizes the advantages of exploiting sensors with multi-angular capabilities. The findings are largely independent of the solar zenith angle at least at high latitudes, an additional advantage for remote sensing applications.
The information content analysis was used to inform the choice of state parameters to be retrieved in sample Levenberg–Marquardt inversions, which were tested on synthetically generated polar scenes for different instrument configurations. The retrievals indeed confirm that mono-angle measurements of total reflectance in the VIS–SWIR (i.e., MODIS) can adequately resolve the grain size in the top layer, while access to more complex descriptors for the snow grains (in our case the aspect ratio and microscale roughness of hexagonal prisms) is achieved by the addition of multi-angle, polarimetric measurements in the NIR–SWIR. Such observations also make it possible to differentiate LAPs in the snow from absorbing aerosol layers, a task that can improve the characterization of processes like aerosol deposition in climate models and, again, albedo simulations.
The findings generally promote the use of the DoLP over the polarized reflectance and indicate that observations from the VIS all the way to the SWIR minimize the uncertainties when attempting to distinguish impurities in snow from absorbing aerosols. The GSA can be extended to LUTs that consider a whole suite of aerosol optical properties, region-specific impurity amounts, and more elaborate mixing schemes (Tanikawa et al., 2020) or optically thin snowpacks with different underlying land cover types. However, the methods and results outlined in this paper provide cryospheric scientists with guidelines for selecting appropriate viewing geometries during data collection and for developing advanced retrieval algorithms applied to airborne and spaceborne data over snow. This perspective is particularly exciting considering the higher accuracies enabled by recent technological progress, like those of the polarimeters on PACE and 3MI.
This Appendix discusses a few aspects that could not be included directly in the main text of the paper without unnecessarily interrupting the flow.
Figure A1 is the same as Fig. 7 but considers a larger range of aerosol optical depth (up to 1.2) and impurity density (up to 10 ppmw), which can occur as a result of burning biomass plumes or extremely “dirty” snow. With higher concentrations of LAPs, the absolute Sobol indices for the visible and NIR wavelengths increase (see Eq. 17) because BC absorption causes large variations in the total reflectance and, to a more limited extent, in the polarized reflectance, especially at larger viewing zenith angles (see also Fig. 3 in Ottaviani, 2022). Larger DoLP signals also lead to a minor increase in the uncertainty threshold, but the list of parameters with identified sensitivity in Fig. 7 remains the same.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f12](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f12-thumb.png)
Figure A1Same as Fig. 7 but for a larger range of aerosol optical depth (up to 1.2) and impurity density (up to 10 ppmw).
A curious aspect concerns the growing sensitivity to reff, as the concentration of snow impurities increases, as shown with the sensitivity studies performed over specific “slices” of the LUT in Fig. A2. Size measurements are achieved by exploiting absorption. Pure snow is highly reflective in the VIS–NIR, and these wavelengths cannot be exploited to retrieve the effective radius of the snow grains, as shown by the solid lines. However, for high concentrations of impurities in the snowpack, the total reflectance and DoLP in the VIS–NIR do show dependence on (dashed lines). The reason for this behavior is that the impurities in the model are externally mixed with snow (Tanikawa et al., 2020) and occupy the empty spaces between grains; the absorption occurring in these negative spaces is therefore partly correlated with the dimension of the crystals.
The differences in the Sobol indices compared to Fig. 7 are only minor. Polarimetric measurements in the SWIR show selective sensitivity to , to DT, and (even more prominently) to . The general conclusions drawn in the main text therefore remain unaffected if a larger range of LAPs is considered. As reported in Sect. 3, the GSA finds no sensitivity to the column-to-plate fraction for mixtures of columns and plates with an aspect ratio (AR) of AR. This fact is explained by the very similar asymmetry parameters of grains with reciprocal aspect ratios (van Diedenhoven et al., 2014a), as shown in Fig. A3 for columns with ART= 19.553 (black) and plates with AR (red).
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f13](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f13-thumb.png)
Figure A2Sensitivity of RI, Rp, and DoLP (columns) to in the VIS–NIR (rows) for pure snow (solid lines) and for snow containing impurities in the top layer with 5 ppmw (dashed lines). The remaining parameters are DT= 0.35, ρT= 0.1 g cm−3, ZT= 3 cm, 320 µm, ARB= 1.0, DB= 0.35, ρB= 0.3 g cm−3, and 0 ppmw. Calculations are for the principal plane and a SZA of 65°.
![https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f14](https://amt.copernicus.org/articles/17/4737/2024/amt-17-4737-2024-f14-thumb.png)
Figure A3Comparison of RI, Rp, and DoLP (columns) for a pure snowpack consisting of column crystals with ART= 19.553 (red) and plate crystals with AR 0.051 (blue) at different wavelengths (rows). The remaining parameters are fixed at 1280 µm, 0.35, ρT= 0.26 g cm−3, ZT= 3 cm, 320 µm, ARB= 1.0, and ρB= 0.4 g cm−3. Again, calculations are for the principal plane and a SZA of 65°.
The radiative transfer code and the code used to perform the GSA, based on the SALib Python package (https://doi.org/10.5281/zenodo.160164, Usher et al., 2016), can be obtained from the corresponding author upon reasonable request.
The complete look-up table of the forward simulations is stored on a dedicated NASA server at GISS in the form of a Python xarray object and can be obtained upon request.
MO led the study. GHM conceived the adaptation of the GSA formalism to the snow dataset and carried out all calculations, assisted by NC. MO wrote the main draft of the paper; all authors contributed to the editing.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Support from the NASA Office of STEM Engagement at the Goddard Institute for Space Studies and the Minority University Research and Education Project (MUREP) to Gabriel Harris Myers is gratefully acknowledged, with additional funding provided by the NASA Remote Sensing Theory program (Lucia Tsaoussi, program manager) which was indispensable to finalize the study. All authors wish to thank Shenglong Wang and Valerio Luccio at the New York University for their patience and guidance in configuring the supercomputing architecture, Igor Geogdzhayev for the management of the GO database of the ice crystal optical properties, and Aaditya Rangan for helpful discussions. MO and NC would like to thank their PhD advisor, Knut Stamnes, for his always precious mentorship.
This research has been funded through the NASA ROSES Remote Sensing Theory program (grant no. 80NSSC21K0569). Gabriel Harris Myers was supported by the NASA internship program and additional funding provided by the Remote Sensing Theory for Earth Science program.
This paper was edited by Alexander Kokhanovsky and reviewed by two anonymous referees.
Alexander, P. M., Tedesco, M., Fettweis, X., van de Wal, R. S. W., Smeets, C. J. P. P., and van den Broeke, M. R.: Assessing spatio-temporal variability and trends in modelled and measured Greenland Ice Sheet albedo (2000–2013), The Cryosphere, 8, 2293–2312, https://doi.org/10.5194/tc-8-2293-2014, 2014.
Antwerpen, R. M., Tedesco, M., Fettweis, X., Alexander, P., and van de Berg, W. J.: Assessing bare-ice albedo simulated by MAR over the Greenland ice sheet (2000–2021) and implications for meltwater production estimates, The Cryosphere, 16, 4185–4199, https://doi.org/10.5194/tc-16-4185-2022, 2022.
Aoki, T., Fukabori, M., Hachikubo, A., Tachibana, Y., and Nishio, F.: Effects of Snow Physical Parameters on Spectral Albedo and Bidirectional Reflectance of Snow Surface, J. Geophys. Res., 105, 10219–10236, 2000.
Aoki, T., Hori, M., Motoyoshi, H., Tanikawa, T., Hachikubo, A., Sugiura, K., Yasunari, T., Storvold, R., Eide, H. A., Stamnes, K., Li, W., Nieke, J., Nakajima, Y., and Takahashi, F.: ADEOS-II/GLI snow/ice products – Part II: Validation results using GLI and MODIS data, Remote Sens. Environ., 111, 274–290, https://doi.org/10.1016/j.rse.2007.02.035, 2007.
Biron, D., Lupi, G., Montini, G., Labate, D., Bruno, U., Melfi, D., Sist, M., Zauli, F., de Leonibus, L.: METOP-SG 3MI (Multi-viewing Multi-channel Multi-polarization Imaging), a powerful observing mission for future operational applications, in: Proc. EUMETSAT Meteorological Satellite Conference, Vienna, Austria, 16–20 September 2013.
Bougamont, M., Bamber, J. L., Ridley, J. K., Gladstone, R. M., Greuell, W., Hanna, E., Payne, A. J., and Rutt, I.: Impact of Model Physics on Estimating the Surface Mass Balance of the Greenland Ice Sheet, Geophys. Res. Lett., 34, L17501, https://doi.org/10.1029/2007GL030700, 2007.
Cairns, B., Russell, E. E., and Travis, L. D.: The Research Scanning Polarimeter: Calibration and ground-based measurements, in: Polarization: Measurement, Analysis, and Remote Sensing II, 18 July 1999, Denver, Col., Proc. SPIE, vol. 3754, 186–196, https://doi.org/10.1117/12.366329, 1999.
Cairns, B., Waquet, F., Knobelspiesse, K., Chowdhary, J., and Deuzé, J.-L.: Polarimetric remote sensing of aerosols over land surfaces, Satellite Aerosol Remote Sensing over Land, edited by: Kokhanovsky, A. A. and De Leeuw, G., Springer-Praxis Books in Environmental Sciences, Springer, 295–325, https://doi.org/10.1007/978-3-540-69397-0_10, 2009.
Chen, N., Li, W., Tanikawa, T., Hori, M., Shimada, R., Aoki, T., and Stamnes, K.: Fast yet accurate computation of radiances in shortwave infrared satellite remote sensing channels, Opt. Express, 25, A649–A664, 2017.
Colgan, W., Box, J. E., Fausto, R. S., van As, D., Barletta, V. R., and Forsberg, R.: Surface Albedo as a Proxy for the Mass Balance of Greenland's Terrestrial Ice, Geol. Surv. Den. Greenl., 31, 93–96, 2014.
Colman, R. A.: Surface albedo feedbacks from climate variability and change, J. Geophys. Res.-Atmos., 118, 2827–2834, 2013.
Dang, C., Fu, Q., and Warren, S. G.: Effect of Snow Grain Shape on Snow Albedo, J. Atmos. Sci., 73, 3573–3583, 2016.
De Haan, J., Bosma, P., and Hovenier, J.: The Adding Method for Multiple Scattering Calculations of Polarized Light, Astron. Astrophys., 183, 371–391, 1987.
Dubovik, O., Holben, B., Eck, T. F., Smirnov, A., Kaufman, Y. J., King, M. D., Tanré, D., and Slutsker, I.: Variability of Absorption and Optical Properties of Key Aerosol Types Observed in Worldwide Locations, J. Atmos. Sci., 59, 590–608, 2002.
Dubovik, O., Li, Z., Mishchenko, M. I., Tanré, D., Karol, Y., Bojkov, B., Cairns, B., Diner, D. J., Espinosa, W. R., Goloub, P., Gu, X., Hasekamp, O., Hong, J., Hou, W., Knobelspiesse, K. D., Landgraf, J., Li, L., Litvinov, P., Liu, Y., Lopatin, A., Marbach, T., Maring, H., Martins, V., Meijer, Y., Milinevsky, G., Mukai, S., Parol, F., Qiao, Y., Remer, L., Rietjens, J., Sano, I., Stammes, P., Stamnes, S., Sun, X., Tabary, P., Travis, L. D., Waquet, F., Xu, F., Yan, C., and Yin, D.: Polarimetric remote sensing of atmospheric aerosols: Instruments, methodologies, results, and perspectives, J. Quant. Spectrosc. Ra., 224, 474–511, https://doi.org/10.1016/j.jqsrt.2018.11.024, 2019.
Dumont, M., Brun, E., Picard, G., Michou, M., Libois, Q., Petit, J.-R., Geyer, M., Morin, S., and Josse, B.: Contribution of Light-Absorbing Impurities in Snow to Greenland's Darkening Since 2009, Nat. Geosci., 7, 509–512, 2014.
Fettweis, X., Hanna, E., Gallée, H., Huybrechts, P., and Erpicum, M.: Estimation of the Greenland ice sheet surface mass balance for the 20th and 21st centuries, The Cryosphere, 2, 117–129, https://doi.org/10.5194/tc-2-117-2008, 2008.
Flanner, M. G. and Zender, C. S.: Linking Snowpack Microphysics and Albedo Evolution, J. Geophys. Res.-Atmos., 111, D12208, https://doi.org/10.1029/2005JD006834, 2006.
Gao, M., Franz, B. A., Knobelspiesse, K., Zhai, P.-W., Martins, V., Burton, S., Cairns, B., Ferrare, R., Gales, J., Hasekamp, O., Hu, Y., Ibrahim, A., McBride, B., Puthukkudy, A., Werdell, P. J., and Xu, X.: Efficient multi-angle polarimetric inversion of aerosols and ocean color powered by a deep neural network forward model, Atmos. Meas. Tech., 14, 4083–4110, https://doi.org/10.5194/amt-14-4083-2021, 2021.
Gao, M., Knobelspiesse, K., Franz, B. A., Zhai, P.-W., Cairns, B., Xu, X., and Martins, J. V.: The impact and estimation of uncertainty correlation for multi-angle polarimetric remote sensing of aerosols and ocean color, Atmos. Meas. Tech., 16, 2067–2087, https://doi.org/10.5194/amt-16-2067-2023, 2023.
Geogdzhayev, I. V. and van Diedenhoven, B.: The effect of roughness model on scattering properties of ice crystals, J. Quant. Spectrosc. Ra., 178, 134–141, 2016.
Hansen, J. and Nazarenko, L.: Soot Climate Forcing via Snow and Ice Albedos, P. Natl. Acad. Sci. USA, 101, 423–428, 2004.
Hansen, J. and Travis, L.: Light Scattering in Planetary Atmospheres, Space Sci. Rev., 16, 527–610, https://doi.org/10.1007/BF00168069, 1974.
Hasekamp, O. P.: Capability of multi-viewing-angle photo-polarimetric measurements for the simultaneous retrieval of aerosol and cloud properties, Atmos. Meas. Tech., 3, 839–851, https://doi.org/10.5194/amt-3-839-2010, 2010.
Hasekamp, O. P., Fu, G., Rusli, S. P., Wu, L., Di Noia, A., aan de Brugh, J., Landgraf, J., Smit, J. M., Rietjens, J., and van Amerongen, A.: Aerosol measurements by SPEXone on the NASA PACE mission: expected retrieval capabilities, J. Quant. Spectrosc. Ra., 227, 170–184, 2018.
He, C., Liou, K.-N., Takano, Y., Yang, P., Qi, L., and Chen, F.: Impact of grain shape and multiple black carbon internal mixing on snow albedo: Parameterization and radiative effect analysis, J. Geophys. Res.-Atmos., 123, 1253–1268, 2018.
Herman, J. and Usher, W.: SALib: An open-source Python library for sensitivity analysis, Journal of Open Source Software, 2, 97, https://doi.org/10.21105/joss.00097, 2017.
Hori, M., Aoki, T., Stamnes, K., and Li, W.: ADEOS-II/GLI snow/ice products – Part III: Retrieved results, Remote Sens. Environ., 111, 291–336, https://doi.org/10.1016/j.rse.2007.01.025, 2007.
Iwanaga, T., Usher, W., and Herman, J.: Toward SALib 2.0: Advancing the accessibility and interpretability of global sensitivity analyses, Socio-Environmental Systems Modelling, 4, 18155, https://doi.org/10.18174/sesmo.18155, 2022.
Jin, Z., Charlock, T., Yang, P., Xie, Y., and Miller, W.: Snow Optical Properties for Different Particle Shapes with Application to Snow Grain Size Retrieval and MODIS/CERES Radiance Comparison Over Antarctica, Remote Sens. Environ., 112, 3563–3581, 2008.
Kawata, Y.: Circular polarization of sunlight reflected by planetary atmospheres, Icarus, 33, 217–232, 1978.
Khan, A. L., Xian, P., and Schwarz, J. P.: Black carbon concentrations and modeled smoke deposition fluxes to the bare-ice dark zone of the Greenland Ice Sheet, The Cryosphere, 17, 2909–2918, https://doi.org/10.5194/tc-17-2909-2023, 2023.
Knobelspiesse, K., Cairns, B., Mishchenko, M., Chowdhary, J., Tsigaridis, K., van Diedenhoven, B., Martin, W., Ottaviani, M., and Alexandrov, M.: Analysis of fine-mode aerosol retrieval capabilities by different passive remote sensing instrument designs, Opt. Express, 20, 21457–21484, 2012.
Kokhanovsky, A.: Light penetration in snow layers, J. Quant. Spectrosc. Ra., 278, 108040, https://doi.org/10.1016/j.jqsrt.2021.108040, 2022.
Lebsock, M. D., L'Ecuyer, T. S., and Stephens, G. L.: Information content of near-infrared spaceborne multiangular polarization measurements for aerosol retrievals, J. Geophys. Res., 112, D14206, https://doi.org/10.1029/2007JD008535, 2007.
Levenberg, K.: A Method for the Solution of Certain Non-Linear Problems in Least Squares, Q. Appl. Math., 2, 164–168, 1944.
Li, W., Stamnes, K., Chen, B., and Xiong, X.: Retrieval of the depth dependence of snow grain size from near-infrared radiances at multiple wavelengths, Geophys. Res. Lett., 28, 1699–1702, 2001.
Libois, Q., Picard, G., France, J. L., Arnaud, L., Dumont, M., Carmagnola, C. M., and King, M. D.: Influence of grain shape on light penetration in snow, The Cryosphere, 7, 1803–1818, https://doi.org/10.5194/tc-7-1803-2013, 2013.
Macke, A., Mueller, J., and Raschke, E.: Single Scattering Properties of Atmospheric Ice Crystals, J. Atmos. Sci., 53, 2813–2825, 1996.
Marbach, T., Phillips, P., Lacan, A., and Schlüssel, P.: The Multi-Viewing, -Channel, -Polarisation Imager (3MI) of the EUMETSAT Polar System - Second Generation (EPS-SG) dedicated to aerosol characterization, in: Proc. SPIE 8889, Sensors, Systems, and Next-Generation Satellites XVII, Dresden, Germany, 23–26 September 2013, 88890I, https://doi.org/10.1117/12.2028221, 2013.
Marquardt, D.: An Algorithm for Least-Squares Estimation of Nonlinear Parameters, SIAM J. Appl. Math., 11, 431–441, 1963.
Newville, M., Stensitzki, T., Allen, D. B., and Ingargiola, A.: LMFIT: Non-Linear Least-Square Minimization and Curve-Fitting for Python, Version 0.8.0, Zenodo [code], https://doi.org/10.5281/zenodo.11813, 2014.
Nolin, A. W. and Dozier, J.: Estimating snow grain size using AVIRIS data, Remote Sens. Environ., 44, 231–238, https://doi.org/10.1016/0034-4257(93)90018-S, 1993.
Ottaviani, M.: Polarization as a Discriminator of Light-Absorbing Impurities in or Above Snow, Front. Remote Sens., 3, 863239, https://doi.org/10.3389/frsen.2022.863239, 2022.
Ottaviani, M., Cairns, B., Rogers, R. R., and Ferrare, R.: Iterative atmospheric correction scheme and the polarization color of alpine snow, J. Quant. Spectrosc. Ra., 113, 789–804, 2012.
Ottaviani, M., van Diedenhoven, B., and Cairns, B.: Photopolarimetric retrievals of snow properties, The Cryosphere, 9, 1933–1942, https://doi.org/10.5194/tc-9-1933-2015, 2015.
Painter, T. H., Dozier, J., Roberts, D. A., Davis, R. E., and Green, R. O.: Retrieval of subpixel snow-covered area and grain size from imaging spectrometer data, Remote Sens. Environ., 85, 64–77, https://doi.org/10.1016/S0034-4257(02)00187-6, 2003.
Painter, T. H., Rittger, K., McKenzie, C., Slaughter, P., Davis, R. E., and Dozier, J.: Retrieval of subpixel snow covered area, grain size, and albedo from MODIS, Remote Sens. Environ., 113, 868–879, https://doi.org/10.1016/j.rse.2009.01.001, 2009.
Rae, J. G. L., Aðalgeirsdóttir, G., Edwards, T. L., Fettweis, X., Gregory, J. M., Hewitt, H. T., Lowe, J. A., Lucas-Picher, P., Mottram, R. H., Payne, A. J., Ridley, J. K., Shannon, S. R., van de Berg, W. J., van de Wal, R. S. W., and van den Broeke, M. R.: Greenland ice sheet surface mass balance: evaluating simulations and making projections with regional climate models, The Cryosphere, 6, 1275–1294, https://doi.org/10.5194/tc-6-1275-2012, 2012.
Räisänen, P., Makkonen, R., Kirkevåg, A., and Debernard, J. B.: Effects of snow grain shape on climate simulations: sensitivity tests with the Norwegian Earth System Model, The Cryosphere, 11, 2919–2942, https://doi.org/10.5194/tc-11-2919-2017, 2017.
Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice. World Scientific, 240 pp., https://doi.org/10.1142/3171, 2000.
Ryan, J. C., Smith, L. C., van As, D., Cooley, S. W., Cooper, M. G., Pitcher, L. H., and Hubbard, A.: Greenland ice sheet surface melt amplified by snowline migration and bare ice exposure, Science Advances, 5, eaav3738, https://doi.org/10.1126/sciadv.aav3738, 2019.
Saltelli, A., Ratto, M., Andres, T., Campolongo, F., Cariboni, J., Gatelli, D., and Tarantola, S.: Global Sensitivity Analysis: The Primer, JohnWiley, Chichester, UK, https://doi.org/10.1002/9780470725184, 2008.
Saltelli, A., Annoni, P., Azzini, I., Campolongo, F., Ratto, M., and Tarantola, S.: Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index, Comput. Phys. Commun., 181, 259–270, https://doi.org/10.1016/j.cpc.2009.09.018, 2010.
Sobol, I. M.: Sensitivity estimates for nonlinear mathematical models, Matematicheskoe Modelirovanie, 2, 112–118, 1990.
Stamnes, K., Li, W., Eide, H., Aoki, T., Hori, M., and Storvold, R.: ADEOS-II/GLI snow/ice products – Part I: Scientific basis, Remote Sens. Environ., 111, 258–273, https://doi.org/10.1016/j.rse.2007.03.023, 2007.
Stamnes, S., Hostetler, C., Ferrare, R., Burton, S., Liu, X., Hair, J., Hu, Y., Wasilewski, A., Martin, W., van Diedenhoven, B., Chowdhary, J., Cetinić, I., Berg, L. K., Stamnes, K., and Cairns, B.: Simultaneous polarimeter retrievals of microphysical aerosol and ocean color parameters from the “MAPP” algorithm with comparison to high-spectral-resolution lidar aerosol and ocean products, Appl. Optics, 57, 2394–2413, 2018.
Tanaka, K., Okamura, Y., Mokuno, M., Amano, T., Yoshida, J.: First year on-orbit calibration activities of SGLI on GCOM-C satellite, in: Proceedings of SPIE AsiaPacific Remote Sensing, Honolulu, Hawaii, USA, 24–26 September 2018, SPIE, 107810Q, https://doi.org/10.1117/12.2324703, 2018.
Tanikawa, T., Kuchiki, K., Aoki, T., Ishimoto, H., Hachikubo, A., Niwano, M., Hosaka, M., Matoba, S., Kodama, Y., Iwata, Y., and Stamnes, K.: Effects of snow grain shape and mixing state of snow impurity on retrieval of snow physical parameters from ground-based optical instrument, J. Geophys. Res.-Atmos., 125, e2019JD031858, https://doi.org/10.1029/2019JD031858, 2020.
Tedesco, M. and Kokhanovsky, A.: The Semi-Analytical Snow Retrieval Algorithm and its Application to Modis Data, Remote Sens. Environ., 111, 228–241, https://doi.org/10.1016/j.rse.2007.02.036, 2007.
Tedesco, M., Fettweis, X., Mote, T., Wahr, J., Alexander, P., Box, J. E., and Wouters, B.: Evidence and analysis of 2012 Greenland records from spaceborne observations, a regional climate model and reanalysis data, The Cryosphere, 7, 615–630, https://doi.org/10.5194/tc-7-615-2013, 2013.
Thackeray, C. W., Qu, X., and Hall, A.: Why do models produce spread in snow albedo feedback?, Geophys. Res. Lett., 45, 6223–6231, https://doi.org/10.1029/2018GL078493, 2018.
Torres, B., Dubovik, O., Fuertes, D., Schuster, G., Cachorro, V. E., Lapyonok, T., Goloub, P., Blarel, L., Barreto, A., Mallet, M., Toledano, C., and Tanré, D.: Advanced characterisation of aerosol size properties from measurements of spectral optical depth using the GRASP algorithm, Atmos. Meas. Tech., 10, 3743–3781, https://doi.org/10.5194/amt-10-3743-2017, 2017.
Usher, W., Herman, J., Whealton, C., Hadka, D., xantares, Rios, F., bernardoct, Mutel, C., and van Engelen, J.: SALib/SALib: Launch! (v1.0.0), Zenodo [code], https://doi.org/10.5281/zenodo.160164, 2016.
van Angelen, J. H., Lenaerts, J. T. M., Lhermitte, S., Fettweis, X., Kuipers Munneke, P., van den Broeke, M. R., van Meijgaard, E., and Smeets, C. J. P. P.: Sensitivity of Greenland Ice Sheet surface mass balance to surface albedo parameterization: a study with a regional climate model, The Cryosphere, 6, 1175–1186, https://doi.org/10.5194/tc-6-1175-2012, 2012.
van den Broeke, M. R., Smeets, C. J. P. P., and van de Wal, R. S. W.: The seasonal cycle and interannual variability of surface energy balance and melt in the ablation zone of the west Greenland ice sheet, The Cryosphere, 5, 377–390, https://doi.org/10.5194/tc-5-377-2011, 2011.
van Diedenhoven, B., Cairns, B., Geogdzhayev, I. V., Fridlind, A. M., Ackerman, A. S., Yang, P., and Baum, B. A.: Remote sensing of ice crystal asymmetry parameter using multi-directional polarization measurements – Part 1: Methodology and evaluation with simulated measurements, Atmos. Meas. Tech., 5, 2361–2374, https://doi.org/10.5194/amt-5-2361-2012, 2012.
van Diedenhoven, B., Cairns, B., Fridlind, A. M., Ackerman, A. S., and Garrett, T. J.: Remote sensing of ice crystal asymmetry parameter using multi-directional polarization measurements – Part 2: Application to the Research Scanning Polarimeter, Atmos. Chem. Phys., 13, 3185–3203, https://doi.org/10.5194/acp-13-3185-2013, 2013.
van Diedenhoven, B., Ackerman, A. S., Cairns, B., and Fridlind, A. M.: A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals, J. Atmos. Sci., 71, 1763–1782, https://doi.org/10.1175/jas-d-13-0205.1, 2014a.
van Diedenhoven, B., Fridlind, A. M., Cairns, B., and Ackerman, A. S.: Variation of Ice Crystal Size, Shape, and Asymmetry Parameter in Tops of Tropical Deep Convective Clouds, J. Geophys. Res.-Atmos., 119, 11809–11825, https://doi.org/10.1002/2014JD022385, 2014b.
Wang, S., Tedesco, M., Alexander, P., Xu, M., and Fettweis, X.: Quantifying spatiotemporal variability of glacier algal blooms and the impact on surface albedo in southwestern Greenland, The Cryosphere, 14, 2687–2713, https://doi.org/10.5194/tc-14-2687-2020, 2020.
Warren, S. G.: Can Black Carbon in Snow Be Detected by Remote Sensing?, J. Geophys. Res.-Atmos., 118, D018476, https://doi.org/10.1029/2012jd018476, 2013.
Warren S. G.: Optical properties of ice and snow, Philos. T. Roy. Soc. A, 377, 20180161, https://doi.org/10.1098/rsta.2018.0161, 2019.
Warren, S. G. and Wiscombe, W.: A Model for the Spectral Albedo of Snow. II: Snow Containing Atmospheric Aerosols, J. Atmos. Sci., 37, 2734–2745, https://doi.org/10.1175/1520-0469(1980)037<2734:AMFTSA>2.0.CO;2, 1980.
Werdell, P. J., Behrenfeld, M. J., Bontempi, P. S., Boss, E., Cairns, B., Davis, G. T., Franz, B. A., Gliese, U. B., Gorman, E. T., Hasekamp, O., Knobelspiesse, K. D., Mannino, A., Martins, J. V., McClain, C. R., Meister, G., and Remer, L. A.: The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission: Status, science, advances, B. Am. Meteorol. Soc., 100, 1775–179, https://doi.org/10.1175/BAMS-D-18-0056.1, 2019.
Wiscombe, W. J. and Warren, S. G.: A Model for the Spectral Albedo of Snow. I: Pure Snow, J. Atmos. Sci., 37, 2712–2733, https://doi.org/10.1175/1520-0469(1980)037<2712:AMFTSA>2.0.CO;2, 1980.
Wu, L., Hasekamp, O., van Diedenhoven, B., and Cairns, B.: Aerosol retrieval from multiangle, multispectral photopolarimetric measurements: importance of spectral range and angular resolution, Atmos. Meas. Tech., 8, 2625–2638, https://doi.org/10.5194/amt-8-2625-2015, 2015.
Xie, Y., Yang, P., Gao, B.-C., Kattawar, G. W., and Mishchenko, M. I.: Effect of Ice Crystal Shape and Effective Size on Snow Bidirectional Reflectance, J. Quant. Spectrosc. Ra., 100, 457–469, https://doi.org/10.1016/j.jqsrt.2005.11.056, 2006.
Zhang, Z., Fu, G., and Hasekamp, O.: Aerosol retrieval over snow using the RemoTAP algorithm, Atmos. Meas. Tech., 16, 6051–6063, https://doi.org/10.5194/amt-16-6051-2023, 2023.