Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
SPI-Based Analyses of Drought Changes over the Past 60 Years in China’s Major Crop-Growing Areas
Next Article in Special Issue
Soil Moisture from Fusion of Scatterometer and SAR: Closing the Scale Gap with Temporal Filtering
Previous Article in Journal
The Impacts of the Ionospheric Observable and Mathematical Model on the Global Ionosphere Model
Previous Article in Special Issue
Ku-, X- and C-Band Microwave Backscatter Indices from Saline Snow Covers on Arctic First-Year Sea Ice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Retrieval of Effective Correlation Length and Snow Water Equivalent from Radar and Passive Microwave Measurements

1
Finnish Meteorological Institute, Erik Palménin aukio 1, FI-00560 Helsinki, Finland
2
Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, No.9 Dengzhuang South Road, Haidian District, Beijing 100094, China
3
Environment and Climate Change Canada, 4905 Dufferin Street, Toronto, ON M3H 5T4, Canada
4
ENVEO IT GmbH, Fürstenweg 176, A-6020 Innsbruck, Austria
5
Institute of Atmospheric and Cryospheric Sciences, University of Innsbruck, A-6020 Innsbruck, Austria
6
Institute of Applied Physics “Nello Carrara”, Via Madonna del Piano, 10-50019 Sesto Fiorentino (FI), Italy
7
WSL Institute for Snow and Avalanche Research SLF, Flüelastrasse 11, CH-7260 Davos Dorf, Switzerland
8
GAMMA Remote Sensing Research and Consulting AG, Worbstr. 225, CH-3073 Gümligen, Switzerland
*
Author to whom correspondence should be addressed.
Remote Sens. 2018, 10(2), 170; https://doi.org/10.3390/rs10020170
Submission received: 18 January 2018 / Accepted: 21 January 2018 / Published: 25 January 2018

Abstract

:
Current methods for retrieving SWE (snow water equivalent) from space rely on passive microwave sensors. Observations are limited by poor spatial resolution, ambiguities related to separation of snow microstructural properties from the total snow mass, and signal saturation when snow is deep (~>80 cm). The use of SAR (Synthetic Aperture Radar) at suitable frequencies has been suggested as a potential observation method to overcome the coarse resolution of passive microwave sensors. Nevertheless, suitable sensors operating from space are, up to now, unavailable. Active microwave retrievals suffer, however, from the same difficulties as the passive case in separating impacts of scattering efficiency from those of snow mass. In this study, we explore the potential of applying active (radar) and passive (radiometer) microwave observations in tandem, by using a dataset of co-incident tower-based active and passive microwave observations and detailed in situ data from a test site in Northern Finland. The dataset spans four winter seasons with daily coverage. In order to quantify the temporal variability of snow microstructure, we derive an effective correlation length for the snowpack (treated as a single layer), which matches the simulated microwave response of a semi-empirical radiative transfer model to observations. This effective parameter is derived from radiometer and radar observations at different frequencies and frequency combinations (10.2, 13.3 and 16.7 GHz for radar; 10.65, 18.7 and 37 GHz for radiometer). Under dry snow conditions, correlations are found between the effective correlation length retrieved from active and passive measurements. Consequently, the derived effective correlation length from passive microwave observations is applied to parameterize the retrieval of SWE using radar, improving retrieval skill compared to a case with no prior knowledge of snow-scattering efficiency. The same concept can be applied to future radar satellite mission concepts focused on retrieving SWE, exploiting existing methods for retrieval of snow microstructural parameters, as employed within the ESA (European Space Agency) GlobSnow SWE product. Using radar alone, a seasonally optimized value of effective correlation length to parameterize retrievals of SWE was sufficient to provide an accuracy of <25 mm (unbiased) Root-Mean Square Error using certain frequency combinations. A temporally dynamic value, derived from e.g., physical snow models, is necessary to further improve retrieval skill, in particular for snow regimes with larger temporal variability in snow microstructure and a more pronounced layered structure.

Graphical Abstract

1. Introduction

The mass of seasonal snow cover, or snow water equivalent (SWE) remains difficult to estimate on a global scale. Observational needs in terms of spatial resolution and product accuracy cannot be met by present satellite, in situ, or model-based data products at the global or regional-scale [1,2]. Global scale EO (Earth-Observation)-based products [3,4] rely on passive microwave sensors, while watershed-scale SWE has been successfully tracked with airborne LiDAR by relating observed snow height to the snow free DEM (Digital Elevation Model), and inferring SWE from the observations by modeling snow density [5]. The cost of timely airborne LiDAR surveys, however, is prohibitive for large-scale applications beyond individual watersheds or regions, so continuing efforts are made to apply Earth Observing satellite sensors for this purpose. Applying passive microwave observations from space for snow cover detection is appealing due to the availability of a long time series of daily observations with near global coverage, extending back almost 40 years. However, estimation of SWE has proved challenging despite several decades of efforts in developing retrieval approaches [6,7]. The main challenges hampering retrieval accuracy are related to the separation of the effect of increasing snow mass from other varying microstructural properties of the snowpack (density, layering, snow structural properties), and mitigating mixed pixel effects in the coarse scale passive microwave observations over heterogeneous landscapes. Existing active microwave sensors are unable to estimate SWE at the global scale and within user requirements because of the lack of current sensors at frequencies higher than X-band. In order to overcome these limits, a dual-band (X- and Ku) SAR mission called CoReH2O (Cold Regions Hydrology High-resolution Observatory, [8]), was proposed as a candidate for the 7th Earth Explorer mission of the European Agency (ESA), with the objective to provide SWE products at a spatial resolution of 200 m, exceeding that of current passive microwave methods. However, following phase-A CoReH2O was not selected for further development.
A priori characterization of snow structural parameters determining the scattering efficiency of microwaves in snow is of primary importance in the accuracy of SWE retrieval algorithms based on radiometer measurements, as this knowledge is required to resolve the total snow mass from observed signal changes [9,10,11]. A key parameter defining the scattering of microwaves has conventionally been the snow grain size, an estimate of the average size of snow grains in the snowpack [12]. The snow grain size has been used to empirically define the rate of microwave extinction in snow [13], which in turn has been applied in a radiative transfer model simulating emission from snow covered ground [14]. An effective grain size can be determined directly from passive microwave measurements, using widely available measurements from weather stations to fix the snow depth for the grain size inversion [15]; an effective grain size determined in this fashion can be related to grain sizes observed in the field [16]. Grain size estimates for application in SWE retrievals can also be obtained by means of applying a model to estimate snow grain metamorphism during the snow season [17].
The snow structural parameter itself remains difficult to measure, and the conversion from the 3D structure to effective model grain size is not unique due to the complex nature of snow grain metamorphism [18,19]. Theoretically-based emission and backscattering models based on the Dense Medium Radiative Transfer (DMRT) theory have assumed snow as a collection of spherical particles, introducing a stickiness parameter to emulate the sintering and clustering of snow grains [20,21,22]. However, it remains difficult to assign properties of snow observed in nature directly to these formulations of DMRT. While measurements of snow specific surface area (SSA) can be applied to estimate the snow optical grain diameter, empirical scaling is required to translate this value to one explaining the observed microwave response [23]. Recent efforts have focused on simulating snow as a bicontinuous medium, simulating the resulting active and passive microwave response with some success [24]. Arguably, statistical parameters such as autocorrelation length of the snow structure in different axial directions are able to describe the snow microstructure with higher fidelity than the conventional measure of grain size [25].
An effective grain size can be used to approximate the scattering behavior of the snowpack using a forward model based on that parameter [15,16]. The term ‘effective’ refers to the fact that the grain size may compensate for simplifications in the model setup (namely the aggregation of a multiple-layered snowpack to a single layer, a practical implementation for operational SWE retrieval schemes, including the one envisaged for CoReH2O), as well as other deficiencies in model input data or model physics. In this study, we present the retrieval of an effective correlation length, which similarly describes the radiative transfer properties of snow by a single parameter. The effective correlation length is retrieved from active and passive microwave observations of naturally accumulated snow over four winter seasons at a test site in northern Finland. We examine the interchangeability of the retrieved correlation length (derived independently from active and passive measurements at different frequencies) for the purpose of initializing the retrieval of SWE from radar observations. Specifically, the study has the following objectives;
  • Retrieve an effective snow correlation length by matching emission and backscattering model predictions to the radiometer and radar measurements, respectively. Examine the seasonal behavior in the observed changes and relate these to physical properties of the snowpack
  • Examine the interchangeability of the active and passive microwave effective correlation length, and determine the sensitivity of these estimates to observation frequency and polarization.
  • Apply an effective correlation length derived from one sensor type (passive) to initialize the retrieval of SWE using the other (active). Compare the impact on SWE retrievals of applying temporally dynamic effective correlation length versus a seasonally constant value, optimized separately for each winter season.
The third objective is meant as a first feasibility demonstration of synergistic active/passive retrieval methods for potential implementation in proposed satellite mission concepts exploiting co-incident radar and radiometer measurements [26].
The study makes use of data acquired from the NoSREx (Nordic Snow Radar Experiment, [27]). The experiment was initiated in November 2009 in support of CoReH2O geophysical algorithm development and included observations of a boreal forest snowpack using both active and passive microwave instruments over four winter seasons.
Section 2 of the study describes the forward model and the methods applied for retrieval of correlation length, and SWE. Section 3 presents the NoSREx datasets. Section 4 presents the main results of the study. These are discussed and the conclusions of the study are given in Section 5.

2. Forward Model and Retrieval Method

2.1. MEMLS3&a Model

The Microwave Emission Model for Layered Snowpacks (MEMLS, [28,29]) is a semi-empirical model utilizing snow correlation length (i.e., the autocorrelation of snow structural variations in a spatial dimension) to define snow microstructure. This can be related to snow-scattering parameters better than the ambiguous snow grain size [25], and is also a parameter that can be objectively estimated from field and laboratory observations. A recent modification to MEMLS expanded the model to also simulate microwave backscattering from snow (MEMLS3&a, [30]). MEMLS3&a allows, by means of model inversion, determination of an effective correlation length describing the average scattering properties of a snowpack using a unified physical approach for both active and passive microwave observations.
The original MEMLS model considers snow as a vertically stacked system of homogeneous snow layers, with each layer characterized by temperature, thickness, liquid water content, salinity, density and correlation length. Radiative transfer in individual layers is calculated using a two-flux approach (fluxes propagating in forward and backward directions), while scattering and absorption coefficients are functions of a six-flux model (including fluxes in perpendicular directions). The absorption coefficient is obtained from density, frequency, temperature, moisture and salinity; the scattering coefficient is determined based on frequency, density, and correlation length. Several empirical formulations have been derived to relate correlation length with the six-flux scattering coefficient [28]; MEMLS also includes an option to use the Improved Born Approximation (IBA) to estimate correlation length, expanding the model for coarse-grained snow [29]. For multiple layers, a sandwich model is applied to determine the effects of internal scattering and reflections at layer interfaces. In the MEMLS3&a model, a prediction of backscattering from snow is calculated based on the specular and diffuse components of snow reflectivity. The specular component is determined from a mean-square slope of slight surface undulations, assuming a Gaussian distribution. The diffuse component is determined from the diffuse component of the total snow reflectivity, which in turn is related to the total emissivity of the snowpack [30].

2.2. Retrieval of Effective Correlation Length

MEMLS3&a was applied in a one-layer configuration to retrieve effective exponential correlation length pex,eff (i.e., a value of correlation length matching model predictions to observations). The approach is analogous to the method introduced by Pulliainen [15] for large-scale passive microwave data using conventional snow grain size. The cost function for iterative inversion can be formulated as
C F ( p ex , eff ) = k = 1 P [ Φ k s ( p ex , eff , S D , x 1 , , x n ) ϕ k s ] 2 σ k s
where SD is the snow depth, pex,eff the effective (exponential) correlation length, P the number of observation channels k, Φ k s the forward model estimate for sensor types (active or passive), xn the model parameters excluding SD and pex,eff, ϕ k s the observed microwave response and σ k s the variance of the combined error of the sensor and the forward model. The same a priori parameters x were applied in both active and passive forward model simulations. The Improved Born Approximation (IBA) was applied to calculate the snow-scattering coefficient [29]. Ground reflectivity (s0v, s0h) at different frequencies corresponds to bare soil with a surface roughness (hrms) of 1 cm and a ground permittivity (εgnd) of 4, calculated using an empirical model [31]. The specular part of reflectivity in MEMLS was considered to be 0.9 across all frequencies. Downwelling sky brightness temperature was estimated using a 55% fractile atmosphere transmissivity model [32]. Snow density (ρsnow) was assigned a constant value of 200 kg m−3, which closely corresponds to the approximate bulk density over the four seasons at the test site [27], and which is also very close to the typical taiga snow density [33]. SD was taken from measured in situ snow depth, and was in baseline retrievals the only temporally variable ancillary input. All other parameters were kept as constants, while pex,eff was the sole free parameter in the iterative optimization (minimization) of F. It should be noted that this approach was chosen as the baseline, although measured values for factors such as air temperature (Tair), ground temperature (Tgnd), ground permittivity and snow density would have been available from the NoSREx experiment. However, our purpose was to emulate potential satellite scale retrievals, where such values are likely not widely available. The retrieved effective correlation length thus accounts also for potential errors arising from the use of these constant ancillary parameters.
The retrievals of effective correlation length using active and passive microwave observations are henceforth labeled pactiveex,eff and ppassiveex,eff, respectively. The ancillary parameters x, as well as other setting of MEMLS3&a, are summarized in Table 1.
To analyze the effect of uncertainties of the approach, retrievals of pactiveex,eff and ppassiveex,eff were repeated using measured values for available ancillary parameters for Tair, Tgnd, εgnd and ρsnow. Tair, Tgnd and εgnd were taken from automated sensors. The ground reflectivity was recalculated based on the measured εgnd, considering the measured value valid across all frequencies, which is a reasonable assumption for frozen ground [34]. The measured density was taken from weekly manual observations of bulk snow density (measurements were made twice a week for the first season of NoSREx). To reduce random variability due to variable location of the manual measurements, and to obtain an estimate for each day, a third order fit was performed on the density data, using the fitted values in retrievals.

2.3. Retrieval of SWE

SWE was determined iteratively from active and passive microwave measurements by inverting the MEMLS3&a model. The cost function used in iterative minimization is analogous to Equation (1), so that
C F ( SWE ) = k = 1 P [ Φ k s ( p ex , eff , SWE , x 1 , , x n ) ϕ k s ] 2 σ k s
where SWE is the sole free parameter and pex,eff is given a priori. The SWE retrieval tests in this study were made using radar backscattering observations; four configurations (labeled Configurations 1 to 4) were used to initialize pex,eff in the retrieval:
  • Configuration 1: An overall average of optimized daily values <pactiveex,eff> was calculated from all retrievals of pactiveex,eff under dry snow conditions for all four seasons. Averages were calculated separately for each channel and combination. These average values of optimizations were used to initialize the respective retrievals of SWE at time t (pex,eff (t) = <pactiveex,eff>).
  • Configuration 2: As in Configuration 1, but the average of optimized daily values <pactiveex,eff> was calculated and applied in SWE retrieval individually for each of the four winter seasons, thus applying seasonal optimization to the retrieval.
  • Configuration 3: For each radar retrieval of SWE at time t, the effective correlation length was acquired from the temporally closest passive microwave retrieval (pex,eff (t) = ppassiveex,eff (t)). As a default, ppassiveex,eff was obtained from 18.7–37 GHz, V-pol, radiometer retrievals.
  • Configuration 4: As in Configuration 3, but the value of ppassiveex,eff was scaled so that pex,eff (t) = β ppassiveex,eff (t). A constant scaling value β was applied in SWE retrieval across all seasons.
The first of the above configurations represents a baseline retrieval where pex,eff is optimized for the site in question based on a climatological average, but without any seasonal sensitivity to potential changes. The second configuration emulates a retrieval where pactiveex,eff has been seasonally optimized, but otherwise kept temporally constant. In practice, this could be achieved e.g., by means of applying e.g., a physical snow model. The contrast with Configuration 1 thus demonstrates the effect of using a seasonally sensitive microstructural indicator, compared to a static climatological average. The third and fourth configurations emulate a scheme where the active microwave retrieval is initialized using effective correlation length derived from passive microwave observations, with the last configuration applying an additional scaling factor, to account for the mean difference in the retrieved pactiveex,eff and ppassiveex,eff. All other input variables of the MEMLS3&a model were kept constant in the retrieval, following Table 1, as in the retrieval of pex,eff.
As with retrievals of pex,eff (see Section 2.2), SWE retrievals were repeated using measured ancillary data for Tair, Tgnd, εgnd and ρsnow. Furthermore, comparative retrievals using passive microwave observations were made adapting Configurations 1 and 2 for SodRad radiometer observations, as the seasonal and overall average of ppassiveex,eff to initialize the retrievals.

3. The NoSREx Campaign

3.1. Microwave Observations

The objective of the NoSREx campaign [27] was to provide a continuous time series of active and passive microwave observations of snow cover in a representative location of boreal forest, covering the complete snow accumulation and ablation period. The campaign lasted four full snow seasons, providing near-continuous observations of snow cover microwave signatures. The campaign hosted two tower-based microwave instruments; SnowScat, a stepped-frequency, fully polarized radar operating on frequencies from X- to Ku band, and the SodRad radiometer, a microwave dual-polarization radiometer system operating from X- to W bands. SnowScat was installed on a tower at the height of 9.6 m overlooking a large forest clearing. SodRad was installed on an adjacent platform at the height of 4 m, overlooking the same general test area. However, the footprints of the instruments were not entirely co-incident (see e.g., [27,30], for a detailed description).
The standard measurement of SnowScat during the NoSREx experiment took place every three to four hours (the measurement sequence was changed to four hours after the first season to allow for scanning of a secondary observation area); co- and cross-polarized backscattering was measured at four incidence angles (30°, 40°, 50°, and 60°). Each elevation scan in the main observation area contained 17 azimuth looks at 6° intervals. Observations at 50° angle of incidence were chosen for this study. The data for each incidence angle were averaged over all measured azimuth directions in order to reduce effects of random speckle and to increase the number of independent looks. The stepped-frequency scan was integrated over three bands of 2 GHz, with center frequencies 10.2, 13.3 and 16.7 GHz. Only the co-polarized vertical measurements (VV-polarization) were used in this study. The inherent reason was that the MEMLS3&a model essentially calculates cross-polarized backscattering from co-polarization using a fixed ratio; thus applying cross-polarized observations was not seen to provide additional benefit.
Brightness temperatures from SodRad at 10.65, 18.7, 21, 37, 90 and 150 GHz (H and V pol) were available at incidence angles ranging from 40 to 60°; however, only the vertically polarized 10.65, 18.7 and 37 GHz channels at 50° were used for this study, as initial retrieval tests with H polarization resulted in larger temporal variability and decreased correlation with pactiveex,eff and so were omitted from further study. SodRad measurements took place in-between SnowScat scans to avoid RFI (radio frequency interference) from the radar signal.

3.2. In Situ Data

The microwave observations were supported by regular manual snow profile measurements as well as extensive array of automated measurements on snow, ground and meteorological parameters [27]. SWE measured from snow profile measurements served as the main reference for SWE retrievals. For the third and fourth season, collected data included vertical profiles of snow specific surface area (SSA), using an infrared laser and integrating sphere [35]. The commercial IceCube instrument was used; the SSA measurements are described by Leppänen et al. [36]. In addition, the SnowMicroPen (SMP; [37]) was used on several occasions—these measurements allowed the quantitative extraction of snow correlation length from the measured vertical profile of penetration resistance [38]. Field data included also conventional measurements of snow grain size [39]; however, these data were not applied here.

3.3. Campaign Summary

Figure 1 depicts the data measured during the four seasons of NoSREx. The data are presented as channel differences to reveal trends associated with snow cover changes. For SnowScat, the channel differences of VV-pol backscattering for 13.3–10.2 GHz, 16.7–10.2 GHz, and 16.7–13.3 GHz are depicted (the data are shown in linear units). For SodRad brightness temperatures, channel differences of 10.65–18.7 GHz, 10.65–37 GHz, and 18.7–37 GHz are shown for the vertical polarization (V-pol). The main reason for applying a channel difference in passive microwave retrievals is to reduce sensitivity to variations in physical temperature affecting microwave emission; in the case of microwave backscattering, the potential benefits of examining channel differences in backscattering are less obvious, but are presented here to determine any potential improvement in retrieval results compared to individual channels.
Also depicted in Figure 1 are in situ measured air temperature (Tair), as well as the 2 cm depth ground temperature (Tgnd) and ground permittivity (εgnd) measured at the NoSREx test site. Manual snow profile measurements were conducted weekly (twice a week during the first season of NoSREx) at the test site. Bulk snow density (ρsnow), measured using a manual snow scale, is shown in the lowest panel, together with the measured snow depth (SD) from an automated sensor. A third order fit to measured ρsnow is shown, which can be used to reduce the effect of random errors in the manual measurements, which in the case of density arises from local-scale variability in density conditions (each density sample was taken from a different location due to the destructive nature of the measurement). Bulk SWE values calculated from the density measurements were used as reference to retrieved SWE in this study.
Clear differences in backscattering signatures between the four NoSREx seasons are apparent; the second season (Figure 1b) shows a large dynamic range in the channel differences involving the 16.7 GHz channel, while the third season (Figure 1c) shows comparatively less temporal dynamics; the first and fourth seasons exhibit a dynamic range between these two extremes. The channel difference of the two lower frequencies (13.3–10.2 GHz), on the other hand, shows stable behavior over all four seasons with very low dynamics relative to changing SWE. It should be noted that for the third season, radar measurements were available only after January 23. The microwave emission shows similarly large dynamics for the second season applying the frequency combinations involving the 37 GHz channel, while the overall response during the third season in particular is muted. The channel difference of the two lower frequencies (10.65 and 18.7 GHz) shows notably low dynamic behavior during the observed winter periods.
It is difficult to directly associate observed increases in microwave backscattering or decreases in microwave emission to snowfall events; rather, significant increases in backscattering at Ku-bands (shown as an increase of the difference of the 16.7 and 10.2 GHz channels, as well as the difference between 16.7 and 13.3 GHz channels) occur often in periods with no significant precipitation, indicating that such increases are due to other changes in snow or soil backscattering characteristics (e.g., in March 2011, Figure 1b). Similarly, significant precipitation events, such as those occurring in February 2013, induce little or no increases in backscattering or dampening of microwave emission, indicating low sensitivity of the microwave response to newly fallen snow. The new snow may in fact work as a dampening layer, reducing backscattering from lower layers with coarser snow microstructure (see e.g., [40]). This may also partly account for the early season decrease in the backscattering signal.

4. Results

4.1. Retrieved Active and Passive Correlation Length

The SnowScat and SodRad measurements were applied in Equation (1) to retrieve an effective correlation length pex,eff for all observation cases in dry snow conditions. These were determined from air temperatures, discarding all cases where air temperature was above −1 °C to allow for margin for unclear cases close to melting point. The measured in situ snow depth was used to initialize the retrieval; as described in Section 2.2, other forward model parameters were kept constant following Table 1. Individual channels, and combinations of two channels were used in the retrieval, mimicking possible satellite sensor configurations. In the case of SnowScat, pex,eff was retrieved using 10.2, 13.3 and 16.7 GHz observations as individual channels, using the combinations 10.2 & 13.3, 10.2 & 16.7, and 13.3 & 16.7 both independently and as a channel difference in the cost function (total of nine different retrieval combinations). VV polarized backscattering observations for each frequency channel or channel combination were used in all cases. Similarly, in the case of SodRad V-pol radiometer measurements, both individual channels, combinations of channels and channel differences were applied (total nine combinations).
Figure 2 depicts the time series of retrieved pactiveex,eff and ppassiveex,eff during the four seasons of NoSREx. Retrievals of pactiveex,eff using the three combinations of channel differences between the 10.2, 13.3 and 16.7 GHz are shown, as well as retrieved ppassiveex,eff using the channel difference combinations between the 10.65, 18.7 and 37 GHz channels. The three time series of pactiveex,eff correspond largely to one another in terms of temporal behavior, in particular in the beginning and end of the snow season, but exhibit clear differences in magnitude; also the channel difference of 13.3–10.2 GHz produces a less dynamic time series of pactiveex,eff . With the exception of 2012–2013, as well as February 2012, the channel difference of 16.7–13.3 GHz gives the largest value of pactiveex,eff, while the channel difference of 13.3–10.2 GHz results in the smallest values. Interestingly, for the last season these relations are reversed. The passive-microwave derived ppassiveex,eff follows some of the temporal behavior of the three pactiveex,eff time series, with notable exceptions in the beginnings of the first and second seasons. This may be related to early-season snow melt-refreeze and consequent metamorphism, which caused early-season scattering behavior that was not captured by the MEMLS3&a model when driven directly using in situ snow measurements [27]. It should be noted that radar observations were missing from the beginning of the third snow season of NoSREx. For the last season, both pactiveex,eff and ppassiveex,eff follow a similar temporal trend; the overall spread in the magnitude of retrievals is also the smallest for this season.
Figure 2 also shows in situ measured values of pex using the IceCube instrument and the SMP. SMP data were available from only a few days throughout the season, while IceCube measurements were conducted weekly during the third and fourth seasons. The measurements of IceCube SSA, made at approximately 5 cm vertical intervals, were averaged throughout the entire snowpack. Depth-weighted averaging was applied, estimating the layer depth represented by the sample from the recorded snow heights. The average values were converted to optical equivalent grain size D0, so that [25]
D 0 = 6 ρ i · SSA
The following empirical linear relations between D0 and pex have been derived [29]
p ex = 0.4 · D 0   for   dendritic   grains
p ex = 0.3 · D 0   for   non dendritic   grains
Similar scaling values, ranging from 0.25 to 0.4, have been presented elsewhere in literature [41]. Here, the scaling value of 0.4 for dendritic grains (Equation (4)) was applied in all cases. The resulting comparison of pex measured using IceCube for the third seasons shows relatively good agreement with ppassiveex,eff (significant R2 values for the channel differences depicted in Figure 2c ranged from 0.48 to 0.61), while pactiveex,eff generally overestimates measured values and showed a lower coefficient of determination (significant R2 values ranged from 0.37 to 0.46). Nevertheless, both ppassiveex,eff and pactiveex,eff as well as in situ measurements show decreasing values of pex in mid-February; this can be explained by snow accumulation in that period (Figure 1c), which decreased the measured average pex. For the fourth season, both pactiveex,eff and ppassiveex,eff overestimate measured values until mid-March 2013 by 20–30% using the channel difference depicted in Figure 2d. However, again the general trend of measured pex is closely replicated by the retrievals, with a generally decreasing trend in pex, pactiveex,eff and ppassiveex,eff until February 2013. Significant R2 values against in situ pex for the fourth season ranged from 0.30 to 0.61 and 0.27 to 0.34, for pactiveex,eff and ppassiveex,eff respectively.
Measurements of pex using IceCube generally agree well with SMP measurements of pex for the third and fourth snow seasons, although the IceCube values show more variability, as depicted by error bars in Figure 2. The results provide some confirmation of the physical representativeness of the retrieved pex, eff; it should, however, be noted the method for active/passive SWE retrieval presented in this study (Section 2.3) does not rely on matching of pex, eff with in situ measurements.
One objective of this study was to examine the interchangeability of the retrieved correlation length derived independently from active and passive observations. Consequently, the relationships between daily averages of pactiveex,eff to ppassiveex,eff were analyzed by calculating basic statistics between parameters retrieved using different channel configurations. Table 2 summarizes the coefficient of determination of pactiveex,eff relative to ppassiveex,eff. Coefficients of determination which are not statistically significant (p > 0.05) are indicated with asterisks. Results are shown as a matrix of the different active and passive channels and channel combinations. When calculating the daily average, obvious erroneous retrievals were ruled out from the analysis (pex,eff < 0.05 mm or pex,eff > 0.5 mm, which correspond to typical minimum and maximum values measured for natural snow [37]). Passive microwave V-pol and active microwave VV pol combinations are shown; retrievals using in situ measured values for model ancillary parameters are shown in parenthesis. The retrievals were also performed using passive microwave H-pol observations as well as a combination of H- and V-pol observations, with largely similar results, the main difference being that H-pol retrievals of ppassiveex,eff exhibited increased temporal variability compared to V-pol retrievals, arising directly from the larger variability of the observed brightness temperature at H-pol. This resulted in most cases in decreased correlation and increased bias and uRMSE against pactiveex,eff. Therefore, only the results for V-pol are presented and discussed further.
Overall, the results show that retrievals derived from single channels show low correlation between frequencies, with the exception of the highest frequencies of 16.7 GHz for active and 37 GHz for passive. Moderate association is apparent also between retrievals at 13.3 (active) and 37 (passive) GHz. Increasingly better results in terms of the coefficient of determination are obtained using multiple channels, and in particular channel differences. The highest overall R2 for most active microwave channels were obtained against the channel differences of 18.7–37 GHz; this combination also typically presented the smallest bias or uRMSE when compared to other passive microwave retrievals. Analysis for individual seasons (not shown) provided higher coefficients of determination between pactiveex,eff and ppassiveex,eff in particular for the third and fourth season of NoSREx, with R2 values exceeding 0.8. The use of measured in situ data in place of constants increased correlations moderately for most frequencies and combinations (values given in parenthesis in Table 2); a notable exception are correlations of pactiveex,eff with pactiveex,eff derived from the 10.6–18.7 GHz channel difference, where almost all correlations were deteriorated by use of in situ data.
Figure 3 shows the retrieved ppassiveex,eff against pactiveex,eff separately for the four seasons of NoSREx, using the channel differences of 16.7–10.2 GHz (active) and 18.7–37 GHz (passive). This comparison gave the best overall coefficient of determination over the four seasons (0.79). However, a consistent underestimation of 0.03–0.06 mm for ppassiveex,eff relative to pactiveex,eff was apparent. Against retrievals with 18.7–37 GHz, pactiveex,eff was on average 15% larger. This value was applied in SWE retrieval Configuration 4 (see Section 2.3) to scale between ppassiveex,eff and pactiveex,eff.
It is notable in Figure 3 that retrievals from both sensor types similarly resolve the differing magnitudes of pex,eff between the seasons; the second season gave on average the largest values, while the third season indicated the lowest values of pex,eff. This confirms that snow microstructure strongly depends on meteorologically driven seasonal variability in snowpack temporal evolution, so a temporal (or at least seasonal) dynamic retrieval of correlation length is desirable. Incidentally, these results concur with the inter-seasonal relations of the average size of snow grains derived from conventional grain size observations (Table 4 in [27]). The second season, with the largest average pex,eff, also exhibits the largest scatter and lowest correlation coefficient between active and passive retrievals.
The average of retrieved pex,eff for the four NoSREx seasons are summarized in Table 3 and Table 4 for active and passive retrievals, respectively, using the applied frequencies and combinations. For most cases, the second season exhibits the largest, and the third season the smallest pex,eff for both active and passive cases; notable exceptions are retrievals involving the lower X-band frequencies (10.2 and 10.65 GHz). The values represent retrievals with constant ancillary model parameters; the use of measured in situ data increased the value of pactiveex,eff on average by 15%, 7% and 5% for 10.2, 13.3 and 16.7 GHz, respectively. For ppassiveex,eff, the induced differences were on average 21%, 10% and −3.0% for 10.65, 18.7 and 37 GHz, respectively. For channel combinations, the effect was typically smaller (2 to 8% increase for pactiveex,eff, and a change of −3 to 16% for ppassiveex,eff). In most cases, the use of ancillary data increased the absolute difference between pactiveex,eff and ppassiveex,eff; the average difference between all retrievals increased from 22 to 40%. The difference between pactiveex,eff and ppassiveex,eff retrieved with the 18.7–37 GHz channel difference increased from 15 to 30%, which was used to scale ppassiveex,eff in Configuration 4, when performing comparative SWE retrievals with measured ancillary in situ information.

4.2. SWE Retrieval Using Radar Observations

SWE was retrieved from SnowScat observations by applying the four retrieval configurations described in Section 2. Retrievals for each season were compared to manual in situ measurements of bulk SWE made weekly at the test site. The retrievals were again constrained to the dry snow season. As described in Section 2, in the first and second configurations, pex,eff in Equation (2) was set to be the overall average, or the seasonal average, of retrieved <pactiveex,eff>, respectively. In the third configuration, pex,eff was taken from the temporally closest retrieval of ppassiveex,eff, thus testing the suitability of passive- microwave derived pex,eff in informing the SWE retrieval from radar observations. The fourth configuration applied this same method, but using also an additional linear scaling factor (20%, as described in Section 4.1), so that pex,eff = β*ppassiveex,eff. A constant value of β = 1.2 was applied over all seasons, based on the average difference between ppassiveex,eff and pactiveex,eff.
Similarly to the retrieval of pex,eff, SWE retrievals were performed for several channels and channel combinations. In Configurations 1 and 2, <pactiveex,eff> was always chosen from the corresponding retrieval at the same channel or combination of radar channels, while in Configurations 3 and 4, ppassiveex,eff was always obtained from the retrieval using the 18.7–37 GHz, V-pol channel difference (retrievals from H-pol were not applied). Retrievals using this passive microwave channel combination showed the best overall agreement with pactiveex,eff retrievals (Table 2).
Figure 4 shows the time series of retrieved SWE for the four seasons of the experiment using retrieval Configurations 1 and 3 (left panels) and 2 and 4 (right panels), using the 13.3 and 16.7 GHz channel combinations. Red lines represent retrievals with Configurations 1 and 2, applying either the overall, or the seasonal average of <pactiveex,eff>; green lines represent retrievals using ppassiveex,eff (Configurations 3 and 4). In comparison to in situ data, retrievals using Configurations 1 and 3 show large biases for all seasons, in particular for the third season (Figure 4e). Applying the seasonally optimized < pactiveex,eff> (Configuration 2, right panels) produces SWE estimates which closely match the in situ measured values, with the most notable discrepancies apparent for the fourth season. Even the weak signal dynamics observed during the third season (see Figure 1c) is sufficient to result in an increase of retrieved SWE which closely matches that of in situ observations (Figure 4f). Similarly, without scaling, the temporally dynamic ppassiveex,eff produces estimates which overestimate in situ measured SWE for all seasons (green lines in left panels. Applying the constant scaling factor β corrects for the overestimation, but actually results in underestimation of SWE in particular for the first and second seasons. This can be explained by the scaling factor β being generally applicable to scale ppassiveex,eff for all frequency combinations of pactiveex,eff; a more optimized fit could be achieved by defining β individually for each channel or combination of channels, as well as for each season.
A comparison of SWE retrievals against in situ measured SWE using the three individual SnowScat frequency bands of 10.2, 13.3 and 16.7 GHz is shown in Figure 5. The data represent all four seasons of NoSREx; the red markers represent SWE retrievals with using the seasonally optimized <pactiveex,eff> (Configuration 2), while green markers denote retrievals with a temporally dynamic but scaled ppassiveex,eff (Configuration 4). Error metrics in terms of the coefficient of determination, bias and uRMS errors are also displayed. Retrievals using 10.2 GHz (Figure 5a) show little or no correlation with measured in situ SWE, denoted by a low coefficient of determination for both configurations. With the higher frequency bands, correlation as well as bias and uRMSE are improved for both configurations. Overall, Configuration 2, using a seasonally optimized, but otherwise temporally static pex,eff = <pactiveex,eff> yields better error metrics compared to Configuration 4. The good performance of Configuration 2 can be expected due to the seasonal optimization of pex,eff; this is reflected in particular by the low retrieval bias. However, the correlation and uRMSE values also indicate good agreement, suggesting that such a static value, if optimized for a given season, is sufficient to account for most of the seasonal dynamics in snow microstructure at the test site.
Figure 6 shows a comparison of retrieved vs. measured SWE retrievals using retrieval Configurations 2 and 4 for different channel combinations, again for all seasons of NoSREx. The combination of the two lower bands (Figure 6a) yields the largest discrepancies against in situ measured SWE; the use of the 10.2 GHz channel even slightly degrades the error metrics when compared to using 13.3 GHz alone (Figure 5b). Of the combinations, the use of 10.2 and 16.7 GHz (Figure 6b) gives for both configurations the best overall error metrics. However, the results against using the 16.7 GHz channel alone are practically the same, thus any benefit of applying the 10.2 GHz channel cannot be shown.
Error metrics for SWE retrieval tests using Configurations 1 and 2 are summarized in Table 5. Using a single value for pex,eff (Configuration 1) results in bias errors for several frequency combinations, particularly notable for the third season. R2 values against the whole dataset of four seasons were below 0.2 for any frequency combination. For individual seasons using the 16.7 GHz channel, however, reasonable correlations were found (namely for the second and third seasons). Using Configuration 2, obtained bias errors are on average below 10 mm for all frequencies and combinations, demonstrating the sensitivity of the retrieval to pex,eff and the necessity for seasonal optimization. All frequency combinations show both reduced uRMSE and an increased R2 value when using Configuration 2, compared to Configuration 1. Furthermore, all retrievals with Configuration 2 using the 16.7 GHz channel give an R2 value exceeding 0.7 on average for the whole dataset, with values exceeding 0.9 for the second season. The fourth season provides the worst results in terms of R2. It can be also observed that, in this case, the addition of a secondary channel to 16.7 provides an improvement in R2 and uRMSE.
Table 6 summarizes similarly the retrievals using Configurations 3 and 4, i.e., by applying pex,eff from passive microwave observations in the retrieval. Relatively large bias errors are obtained when using ppassiveex,eff directly without scaling (Configuration 3); however, R2 values in the order of 0.6 are still obtained for several frequency combinations on average, with R2 values exceeding 0.8 for several frequency combinations for the second season, and the 16.7–13.3 GHz combination for the fourth season. This provides an indication that the passive microwave derived pex,eff is able to improve radar retrievals, at least when compared to applying a single value for the whole dataset (i.e., Configuration 1, Table 5). Applying a scaling factor to ppassiveex,eff (Configuration 4) further improves the error metrics against in situ SWE, in particular in terms of reduced bias and uRMSE. R2 values for some frequency combinations also show improvement (e.g., combinations of 10.2 & 16.7 GHz and 13.3 & 16.7 GHz), while e.g., retrievals using 16.7 GHz and the channel difference of 16.7–10.2 GHz shows slight deterioration.
Retrievals are generally overestimated during the early winter, in particular for the first and fourth seasons (see Figure 4). This is a direct result of increased backscattering observed in early winter, possibly associated with formation of basal crusts due to snow melt-refreeze events (Figure 1; see also [27,42]). Applying the channel difference mitigates this feature for some seasons and channel combinations; as an example, Figure 7 demonstrates the retrieval for the fourth season of NoSREx using the channel difference of 16.7–13.3 GHz. The overestimation in early winter retrievals for Configuration 2 (red line) is reduced compared to applying the channels individually in the cost function (compare to Figure 4h). While helpful in this case, applying the channel difference in place of individual channels in the retrieval was not found to universally improve retrieval error statistics (see Table 5).
SWE retrievals were repeated using available measured in situ data (Tair, Tgnd, εgnd, and ρsnow). Values of pex,eff applied in the various SWE retrieval configurations were also derived from retrievals involving measured in situ data. For Configurations 1, 2 and 4 retrievals involving the 16.7 GHz channel (including channel combinations) yielded on average a 27 to 75% reduction in bias errors, a reduction in uRMSE from 5 to 14%, and a 7 to 15% improvement in the coefficient of determination. For Configuration 3, however, the use of in situ data resulted in an increase of retrieval bias and uRMSE by 33 and 16%, respectively, while the coefficient of determination was marginally reduced, compared to retrievals with static ancillary parameters.
Retrievals of SWE using the SodRad radiometer observations generally showed less skill compared to radar retrievals. Figure 8 shows retrieved against measured SWE from the 10.65–36.5 GHz channel difference, which yielded the best overall error metrics. Results from Configurations 1 and 2, adapted for radiometer observations, are shown. Use of the seasonal average of ppassiveex (Configuration 2) can be seen to reduce both the bias and uRMSE, and increase the coefficient of determination compared to Configuration 1, where the overall average of ppassiveex was used. However, the error metrics are inferior compared to any radar channel combination involving the 16.7 GHz channel and using the seasonal average of pactiveex (Configuration 2 in Table 5), with the exception of the larger bias error of the 13.3–16.7 GHz channel difference retrievals.

5. Discussion

Results in Section 4.1 demonstrated that at certain frequencies and combinations of frequencies, active and passive -derived pex,eff were shown to be temporally correlated (Table 2). This can be expected, as similar physics drive the backscattering and emission signal behavior, linked to varying scattering and absorption properties induced by precipitation events, and snowpack metamorphism over time. On average, obtained correlations were higher between higher frequencies which are generally more sensitive to the volume properties of snow cover, while lower X-band active and passive frequencies were found to be mostly uncorrelated, affirming that the underlying soil and not dry snow-scattering properties govern the signal at these frequencies and snow depths. However, even at higher frequencies a seasonally consistent bias between active and passive -derived pex,eff was found (Figure 3), indicating deficiencies in either the model physics or the applied ancillary data. Possible reasons for the bias include the omission of snow layering in the applied one-layer model configuration which affects in particular the highest 37 GHz passive microwave channel simulations, as well as deficiencies in the applied ground reflectivity model [31], which is not verified for active microwave observations beyond this paper and tests conducted by Proksch et al. [30].
Considering practical satellite scale retrievals, obtaining seasonal characteristics for snow microstructure and thus scattering properties could be realized by assigning e.g., seasonally or regionally dynamic values for correlation length, derived from snow climatology or physical snow models. However, as demonstrated in this study, passive microwave observations may also be of use in providing a priori information on snow microstructure for radar retrievals. This implies a method analogous to the one introduced in [15], characterizing snow-scattering efficiency over regions where in situ data is available and extending this over larger areas by means of spatio-temporal interpolation. In such a scheme, an effective correlation length derived at large scales with daily coverage from passive microwave observations could initialize active microwave retrievals of SWE from backscattering observations at suitable frequencies. This naturally implies assuming that differences in the scale and geometry of observations do not result in detrimental differences affecting the interchangeability of ppassiveex,eff and pactiveex,eff.; this matter was not accounted for in this study and should be investigated before the method proposed here could be applied at the satellite scale. Forest cover in particular poses a challenge in relating coarse scale passive microwave observations to high resolution SAR. However, new methods for mitigating forest canopy effects in snow parameter retrieval from both radar and passive microwave measurements have been developed using recent experimental data [43].
Taiga snow typically has a poorly defined internal layer structure compared to other classes of seasonal snow: snow grains near the base are larger than snow grains near the surface (with depth hoar by the end of the season) but the internal boundaries between layers are weak. This means a single snow layer was sufficient for the MEMLS3&a simulations in this study, in particular for late winter. This may not be the case for strongly layered snowpacks that are found in open tundra and prairie environments where wind processes strongly influence snow distribution. Tundra snow is typically composed of fine-grained wind slab layers overlying depth hoar composed of large faceted grains. Under these conditions, it is likely that two layer snow simulations would be necessary [44]. However, also in the case of NoSREx data, early season radar observations, and thus retrievals of both pactiveex,eff and SWE, were assumed to be influenced by melt-refreeze crusts in the base of the snowpack; causing the high backscattering values seen in the early season (Figure 1). Accounting for these would require use of a multi-layer model to properly represent the snow-scattering properties.
Reasonable estimates of SWE (Section 4.2) were achieved by inverting active microwave backscattering observations, provided that the mean level of snowpack scattering efficiency for each season was known (i.e., retrieval Configuration 2, see Table 5). This is partly incidental, as increases in snow mass and evolution of snow microstructure both occur over the season, while the direct impact of newly precipitated snow on the microwave signature was in many cases very small (see Figure 1); the new snow may also dampen the backscattering signature from lower snow layers thus compensating for any potential increase in total volume scattering from the new snow layer [40]. Snow regimes with more short-term temporal changes in microstructure may necessitate temporally dynamic estimates of the microstructure to achieve the same retrieval skill. In an operational context, assigning seasonally constant estimates of scattering efficiency is unrealistic if short latency is required (a seasonal average can be assigned only after the snow season is over). Tests with Configuration 2 nevertheless demonstrate that temporal dynamics may not in all cases be critical for SWE retrieval performance, if the average magnitude of scattering efficiency is correctly estimated. Retrievals of SWE using passive microwave observations produced generally inferior results (Figure 8), indicating that radar retrievals may have the potential to improve the absolute accuracy of current SWE products from space.
The use of temporally variant, measured ancillary information for snow and soil properties was shown to generally improve the coefficient of determination between correlation length derived active and passive measurements (Table 2), as well as error metrics of SWE retrievals, compared to retrievals using constant values. The exception was retrieval Configuration 3, where ppassiveex was used without scaling to initialize SWE retrievals from radar. However, the average improvement in e.g., the uRMSE were less than 14% in any configuration when applying the 16.7 GHz channel. The results show that while the use of accurate ancillary information is preferable, retrieval skill may still be acceptable when using constant values derived from e.g., climatology, when measured or modeled data is unavailable.
Similar retrieval results for SWE were obtained by using an absorption-loss-based methodology using the same NoSREx dataset [45]; that study also showed that a priori definition of snow-scattering characteristics (in the form of single scattering albedo) was essential to achieve reasonable retrievals. Nevertheless, further improvement may be achieved by applying more sophisticated retrieval schemes, including balancing or constraining the parameterization of the retrieval (following e.g., [45]), as well as by estimating the microstructural evolution of snow based on a physical model. Direct coupling with physical snow models has been previously demonstrated using the same dataset [46], including retrieval of SWE using a coupled model. A demonstration of different microstructural parameterizations and their use in emission models using a suite of different emission and physical models is given in [47]. Similar schemes have been demonstrated previously for SWE estimation [48] and for forward model simulations [41,49] using other datasets.
Information on SWE can be coupled with land surface models by assimilation of snow mass estimates [50] or by direct radiance assimilation [51,52,53]. With respect to direct assimilation of radiances, this study indicates that the MEMLS3&a model may be applied to simulate both active and passive microwave data streams, provided a suitable proxy of snow microstructure (i.e., correlation length) is produced by the land surface model (e.g., [54]).

6. Conclusions

This study investigated the seasonal behavior of snow bulk scattering properties by means of retrieving a proxy variable, the effective snow correlation length, from active and passive microwave observations. Tower-based measurements from the NoSREx campaign, representing four years of differing snow conditions in a taiga environment, were used. A numerical inversion of the MEMLS3&a backscattering and emission model was used to retrieve an effective correlation length matching model estimates to observations. Snow depth was used as the only temporally variable ancillary input, mimicking a similar scheme applied for satellite passive microwave SWE retrievals, in which conventional snow grain size is used as a proxy [4]. The effective snow-scattering efficiency derived in this fashion accounts also for deficiencies in e.g., ancillary data supplied to the forward model used in the inversion, as well as physics of the forward model itself. If these deficiencies are severe, the retrieved effective snow parameters (either effective grain size or correlation length) may lose any physical relation to actual snow microstructural conditions. Nevertheless, for the last two seasons when weekly observations of SSA profiles in snow were available, the retrieved values showed similar temporal trends with MEMLS3&a retrievals of pex,eff, (Figure 2c,d). A bias was apparent for the fourth NoSREx season (Figure 2d), but these discrepancies can be due to uncertainties in the SSA measurement, the empirical relation applied to convert SSA to pex, and the use of a one-layer model for the snowpack to retrieve pex,eff, which complicates the relation to physical snow microstructure in the presence of a vertically inhomogeneous snowpack. Similar results for the second season of NoSREx were obtained using the Helsinki University of Technology (HUT) snow emission model [14] and conventional snow grain size measured in situ (see Figure 9 in [16]). This indicates that effective snow microstructural parameters may under certain conditions bear relation to measured physical snow characteristics.
This study also demonstrated SWE retrievals from radar observations, parameterizing the retrieval by a temporally static effective value for snow-scattering efficiency (Configurations 1 and 2), and by a temporally dynamic estimate, derived from passive microwave observations (Configurations 3 and 4). Applying seasonal scaling to the average scattering efficiency (Configuration 2), retrievals yielded on average highly correlated estimates against in situ observations, when higher Ku band radar (16.7 GHz) observations were applied either in a single frequency retrieval (R2 values between 0.59 and 0.91 for different seasons), or in combination with lower frequencies (R2 values between 0.48 and 0.92). It is notable that high correlations were achieved despite the value of snow correlation length being static for each season. However, when attempting to apply a single value for the whole dataset (Configuration 1), error metrics for the whole dataset deteriorated with the coefficient of determination falling below 0.2 regardless of the frequencies applied in the retrieval, highlighting the need for at least seasonal characterization of the scattering efficiency. For individual seasons, R2 values of up to 0.92 were still found, but retrieval skill was hampered by large biases (Table 5). This is consistent with different bulk and microstructural properties of the snowpack influencing the radar backscatter response each season. Applying an effective correlation length derived from passive microwave measurements to initialize the active microwave retrievals (Configurations 3 and 4) was shown to provide improvement compared to using a single constant value (i.e., Configuration 1). This indicates that passive microwave observations may, in some cases, serve to initialize retrievals using radar backscattering, in particular at Ku band frequencies.
In summary, the presented results suggest that synergistic active-passive microwave observations, such as those provided by a new radar mission concept [26], may enable the retrieval of terrestrial SWE via the method of assigning an effective parameter for snow-scattering efficiency based on passive microwave retrievals. Such retrieval products are available already from current passive microwave satellite estimates [4]. Further work includes demonstrating the synergistic retrieval method beyond plot scale observations. Being sensitive to snow volumetric properties, synergistic passive microwave and Ku-band radar measurements also support land surface modeling within operational forecast systems through radiance-based assimilation of backscatter.

Acknowledgments

The staff at FMI-ARC are acknowledged for collection of in situ data and operation of microwave instruments. Staff from WSL Institute for Snow and Avalanche Research SLF participated in field measurements and performed analysis of SMP measurements. The work was supported by the European Space Agency (ESA ESTEC contracts 22671/09/NL/JA/ef and 22829/09/NL/JC). We thank the three anonymous reviewers for providing constructive comments which helped to improve the manuscript.

Author Contributions

Juha Lemmetyinen, Helmut Rott, Andreas Wiesmann, Martin Schneebeli, Leena Leppänen, Anna Kontu and Jouni Pulliainen designed and performed the experiments; Chris Derksen, Giovanni Macelloni and Josh King contributed to data analysis; Juha Lemmetyinen wrote the paper with contributions from all authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mudryk, L.R.; Derksen, C.; Kushner, P.J.; Brown, R. Characterization of Northern Hemisphere snow water equivalent datasets, 1981–2010. J. Clim. 2015, 28, 8037–8051. [Google Scholar] [CrossRef]
  2. Larue, F.; Royer, A.; De Sève, D.; Langlois, A.; Roy, A.; Brucker, L. Validation analysis of the GlobSnow-2 database over an eco-climatic latitudinal gradient in Eastern Canada. Remote Sens. Environ. 2016, 194, 264–277. [Google Scholar] [CrossRef]
  3. Kelly, R.E.J. The AMSR-E Snow depth algorithm: Description and initial results. J. Remote Sens. Soc. Jpn. 2009, 29, 307–317. [Google Scholar]
  4. Takala, M.; Luojus, K.; Pulliainen, J.; Derksen, C.; Lemmetyinen, J.; Kärnä, J.-P.; Koskinen, J.; Bojkov, B. Estimating northern hemisphere snow water equivalent for climate research through assimilation of space-borne radiometer data and ground-based measurements. Remote Sens. Environ. 2011, 115, 3517–3529. [Google Scholar] [CrossRef]
  5. Painter, T.H.; Berisford, D.F.; Boardman, J.W.; Bormann, K.J.; Deems, J.S.; Gehrke, F.; Hedrick, A.; Joyce, M.; Laidlaw, R.; Marks, D.; et al. The Airborne Snow Observatory: Fusion of scanning LiDAR, imaging spectrometer, and physically-based modeling for mapping snow water equivalent and snow albedo. Remote Sens. Environ. 2016, 184, 139–152. [Google Scholar] [CrossRef]
  6. Foster, J.; Sun, C.; Walker, J.; Kelly, R.; Chang, A.; Dong, J.; Powell, H. Quantifying uncertainty in passive microwave snow water equivalent observations. Remote Sens. Environ. 2005, 94, 187–203. [Google Scholar] [CrossRef]
  7. Clifford, D. Global estimates of snow water equivalent from passive microwave instruments: History, challenges and future developments. Int. J. Remote Sens. 2010, 31, 3707–3726. [Google Scholar] [CrossRef]
  8. Rott, H.; Yueh, S.H.; Cline, D.W.; Duguay, C.; Essery, R.; Haas, C.; Hélière, F.; Kern, M.; Macelloni, G.; Malnes, E.; et al. Cold regions hydrology high-resolution observatory for snow and cold land processes. Proc. IEEE 2010, 98, 752–765. [Google Scholar] [CrossRef] [Green Version]
  9. Dai, L.Y.; Che, T.; Wang, J.; Zhang, P. Snow depth and snow water equivalent estimation from AMSR-E data based on a priori snow characteristics in Xinjiang, China. Remote Sens. Environ. 2012, 127, 14–29. [Google Scholar] [CrossRef]
  10. Davenport, I.J.; Sandells, M.J.; Gurney, R.J. The effects of variation in snow properties on passive microwave snow mass estimation. Remote Sens. Environ. 2012, 118, 168–175. [Google Scholar] [CrossRef]
  11. Durand, M.; Liu, D. The need for prior information in characterizing snow water equivalent from microwave brightness temperatures. Remote Sens. Environ. 2012, 126, 248–257. [Google Scholar] [CrossRef]
  12. Fierz, C.; Armstrong, R.L.; Durand, Y.; Etchevers, P.; Greene, E.; McClung, D.M.; Nishimura, K.; Satyawali, P.K.; Sokratov, S.A. The International Classification for Seasonal Snow on the Ground; IACS Contribution no. 1; UNESCOIHP, I HP-VII Technical Documents in Hydrology 83; UNESCO: Paris, France, 2009. [Google Scholar]
  13. Hallikainen, M.T.; Ulaby, F.; Deventer, T. Extinction behavior of dry snow in the 18- to 90-GHz range. IEEE Trans. Geosci. Remote Sens. 1987, GE-25, 737–745. [Google Scholar] [CrossRef]
  14. Pulliainen, J.T.; Grandell, J.; Hallikainen, M.T. HUT snow emission model and its applicability to snow water equivalent retrieval. IEEE Trans. Geosci. Remote Sens. 1999, 37, 1378–1390. [Google Scholar] [CrossRef]
  15. Pulliainen, J. Mapping of snow water equivalent and snow depth in boreal and sub-arctic zones by assimilating space-borne microwave radiometer data and ground-based observations. Remote Sens. Environ. 2006, 101, 257–269. [Google Scholar] [CrossRef]
  16. Lemmetyinen, J.; Derksen, C.; Toose, P.; Proksch, M.; Pulliainen, J.; Kontu, A.; Rautiainen, K.; Seppänen, J.; Hallikainen, M. Simulating seasonally and spatially varying snow cover brightness temperature using HUT snow emission model and retrieval of a microwave effective grain size. Remote Sens. Environ. 2015, 156, 71–95. [Google Scholar] [CrossRef]
  17. Kelly, R.; Chang, A.; Tsang, L.; Foster, J. A prototype AMSR-E global snow area and snow depth algorithm. IEEE Trans. Geosci. Remote Sens. 2003, 41, 230–242. [Google Scholar] [CrossRef]
  18. Pinzer, B.; Schneebeli, M. Snow metamorphism under alternating temperature gradients: Morphology and recrystallization in surface snow. Geophys. Res. Lett. 2009, 36, L23503. [Google Scholar] [CrossRef]
  19. Pinzer, B.; Schneebeli, M.; Kaempfer, T. Vapor flux and recrystallization during dry snow metamorphism under a steady temperature gradient as observed by time-lapse micro-tomography. Cryosphere 2012, 6, 1141–1155. [Google Scholar] [CrossRef] [Green Version]
  20. Tsang, L.; Chen, C.-T.; Chang, A.T.C.; Guo, J.; Ding, K.-H. Dense media radiative transfer theory based on quasi-crystalline approximation with applications to microwave remote sensing of snow. Radio Sci. 2000, 35, 731–749. [Google Scholar] [CrossRef]
  21. Tsang, L.; Pan, J.; Liand, D.; Li, Z.; Cline, D.W.; Tan, Y. Modeling active microwave remote sensing of snow using dense medium radiative transfer (DMRT) theory with multiple-scattering effects. IEEE Trans. Geosci. Remote Sens. 2007, 45, 990–1004. [Google Scholar] [CrossRef]
  22. Picard, G.; Brucker, L.; Roy, A.; Dupont, F.; Fily, M.; Royer, A. Simulation of the microwave emission of multi-layered snowpacks using the dense media radiative transfer theory: The DMRT-ML model. Geosci. Model Devel. 2013, 6, 1061–1078. [Google Scholar] [CrossRef]
  23. Löwe, H.; Picard, G. Microwave scattering coefficient of snow in MEMLS and DMRT-ML revisited: The relevance of sticky hard spheres and tomography-based estimates of stickiness. Cryosphere 2015, 9, 2101–2117. [Google Scholar] [CrossRef] [Green Version]
  24. Tan, S.; Chang, W.; Tsang, L.; Lemmetyinen, J.; Proksch, M. Modeling both active and passive microwave remote sensing of snow using dense media radiative transfer (DMRT) theory with multiple scattering and backscattering enhancement. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4418–4430. [Google Scholar] [CrossRef]
  25. Mätzler, C. Relation between grain-size and correlation length of snow. J. Glaciol. 2002, 162, 461–466. [Google Scholar] [CrossRef]
  26. Derksen, C.; Lemmetyinen, J.; King, J.; Garnaud, C.; Belair, S.; Lapointe, M.; Crevier, Y.; Girard, R.; Burbidge, G.; Marquez-Martinez, J.; et al. A new dual-frequency Ku-band radar mission concept for cryosphere applications. Proc. EUSAR 2018. submitted. [Google Scholar]
  27. Lemmetyinen, J.; Kontu, A.; Pulliainen, J.; Vehviläinen, J.; Rautiainen, K.; Wiesmann, A.; Mätzler, C.; Werner, C.; Rott, H.; Nagler, T.; et al. Nordic snow radar experiment. Geosci. Instrum. Methods Data Syst. 2016, 5, 403–415. [Google Scholar] [CrossRef]
  28. Wiesmann, A.; Mätzler, C. Microwave emission model of layered snowpacks. Remote Sens. Environ. 1999, 70, 307–316. [Google Scholar] [CrossRef]
  29. Mätzler, C.; Wiesmann, A. Extension of the microwave emission model of layered snowpacks to coarse grained snow. Remote Sens. Environ. 1999, 70, 318–326. [Google Scholar] [CrossRef]
  30. Proksch, M.; Mätzler, C.; Wiesmann, A.; Lemmetyinen, J.; Schwank, M.; Löwe, H.; Schneebeli, M. MEMLS3&a: Microwave emission model of layered snowpacks adapted to include backscattering. Geosci. Model Dev. 2015, 8, 2611–2626. [Google Scholar] [CrossRef]
  31. Wegmüller, U.; Mätzler, C. Rough bare soil reflectivity model. IEEE Trans. Geosci. Remote Sens. 1999, 37, 1391–1395. [Google Scholar] [CrossRef]
  32. Pulliainen, J.; Kärnä, J.-P.; Hallikainen, M.T. Development of geophysical retrieval algorithms for the MIMR. IEEE Trans. Geosci. Remote Sens. 1993, 31, 268–277. [Google Scholar] [CrossRef]
  33. Sturm, M.; Taras, B.; Liston, G.E.; Derksen, C.; Jonas, T.; Lea, J. Estimating snow water equivalent using snow depth data and climate classes. J. Hydrometeorol. 2010, 11, 1380–1394. [Google Scholar] [CrossRef]
  34. Hallikainen, M.T.; Ulaby, F.T.; Dobson, M.C.; El-Rayes, M.A.; Wu, L.-K. Microwave dielectric behavior of wet Soil—Part I: Empirical models and experimental observations. IEEE Trans. Geosci. Remote Sens. 1985, 23, 25–33. [Google Scholar] [CrossRef]
  35. Gallet, J.-C.; Domine, F.; Zender, C.; Picard, G. Measurement of the specific surface area of snow using infrared reflectance in an integrating sphere at 1310 and 1550 nm. Cryosphere 2009, 3, 167–182. [Google Scholar] [CrossRef] [Green Version]
  36. Leppänen, L.; Kontu, A.; Vehviläinen, J.; Lemmetyinen, J.; Pulliainen, J. Comparison of traditional and optical grain size field measurements with SNOWPACK simulations in a taiga environment. J. Glaciol. 2015, 61, 151–162. [Google Scholar] [CrossRef]
  37. Schneebeli, M.; Pielmeier, C.; Johnson, J.B. Measuring snow microstructure and hardness using a high resolution penetrometer. Cold Reg. Sci. Technol. 1999, 30, 101–114. [Google Scholar] [CrossRef]
  38. Proksch, M.; Löwe, H.; Schneebeli, M. Density, specific surface area, and correlation length of snow measured by high-resolution penetrometry. J. Geophys. Res. Earth Surf. 2015, 120, 346–362. [Google Scholar] [CrossRef]
  39. Leppänen, L.; Kontu, A.; Hannula, H.-R.; Sjöblom, H.; Pulliainen, J. Sodankylä manual snow survey program. Geosci. Inst. Methods Data Syst. Discuss. 2016, 5, 163–179. [Google Scholar] [CrossRef]
  40. Paloscia, S.; Pettinato, S.; Santi, E.; Valt, M. COSMO-SkyMed image investigation of snow features in alpine environment. Sensors 2017, 17, 84. [Google Scholar] [CrossRef] [PubMed]
  41. Brucker, L.; Royer, A.; Picard, G.; Langlois, A.; Fily, M. Hourly simulations of the microwave brightness temperature of seasonal snow in Quebec, Canada, using a coupled snow evolution-emission model. Remote Sens. Environ. 2011, 115, 1966–1977. [Google Scholar] [CrossRef]
  42. Lin, C.; Rommen, B.; Floury, N.; Schûttemeyer, D.; Davidson, M.W.; Kern, M.; Kontu, A.; Lemmetyinen, J.; Pulliainen, J.; Wiesmann, A.; et al. Active Microwave Scattering Signature of Snowpack-Continuous Multiyear SnowScat Observation Experiments. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 3849–3869. [Google Scholar] [CrossRef]
  43. Cohen, J.; Lemmetyinen, J.; Pulliainen, J.; Heinilä, K.; Montomoli, F.; Seppänen, J.; Hallikainen, M.T. The effect of boreal forest canopy in satellite snow mapping—A multisensor analysis. IEEE Trans. Geosci. Remote Sens. 2015, 52, 3275–3288. [Google Scholar] [CrossRef]
  44. King, J.; Derksen, C.; Toose, P.; Langlois, A.; Larsen, C.; Lemmetyinen, J.; Marsh, J.; Montpetit, B.; Roy, A.; Rutter, N.; et al. The influence of snow microstructure on dual-frequency radar measurements in a tundra environment. Remote Sens. Environ. 2018. submitted. [Google Scholar]
  45. Cui, Y.; Xiong, C.; Lemmetyinen, J.; Shi, J.; Jiang, L.; Peng, B.; Li, H.; Zhao, T.; Ji, D.; Hu, T. Estimating Snow Water Equivalent with backscattering at X and Ku band based on absorption loss. Remote Sens. 2016, 8, 505. [Google Scholar] [CrossRef]
  46. Kontu, A.; Lemmetyinen, J.; Vehviläinen, J.; Leppänen, L.; Pulliainen, J. Coupling SNOWPACK-modeled grain size parameters with the HUT snow emission model. Remote Sens. Environ. 2017, 194, 33–47. [Google Scholar] [CrossRef]
  47. Sandells, M.; Essery, R.; Rutter, N.; Wake, L.; Leppänen, L.; Lemmetyinen, J. Microstructure representation of snow in coupled snowpack and microwave emission models. Cryosphere 2017, 11, 229–246. [Google Scholar] [CrossRef]
  48. Langlois, A.; Royer, A.; Derksen, C.; Montpetit, B.; Dupont, F.; Goïta, K. Coupling the snow thermodynamic model SNOWPACK with the microwave emission model of layered snowpacks for subarctic and arctic snow water equivalent retrievals. Water Resour. Res. 2012, 48, 1–14. [Google Scholar] [CrossRef]
  49. Picard, G.; Brucker, L.; Fily, M.; Gallée, H.; Krinner, G. Modelling time series of microwave brightness temperature in Antarctica. J. Glaciol. 2009, 55, 537–551. [Google Scholar] [CrossRef]
  50. De Lannoy, G.J.M.; Reichle, R.H.; Houser, P.R.; Arsenault, K.R.; Verhoest, N.E.C.; Pauwels, V.N.R. Satellite-scale snow water equivalent assimilation into a high-resolution land surface model. J. Hydrometeorol. 2010, 11, 352–369. [Google Scholar] [CrossRef]
  51. Durand, M.; Kim, E.J.; Marguilis, S.A. Radiance assimilation shows promise for snowpack characterization. Geophys. Res. Lett. 2009, 36, L02503. [Google Scholar] [CrossRef]
  52. Kwon, Y.; Yang, Z.-L.; Hoar, T.J.; Toure, A.M. Improving the radiance assimilation performance in estimating snow water storage across snow and land-cover types in North America. J. Hydrometeorol. 2017, 18, 651–668. [Google Scholar] [CrossRef]
  53. Andreadis, K.M.; Lettenmaier, D.P. Implications of representing snowpack stratigraphy for the assimilation of passive microwave satellite observations. J. Hydrometeorol. 2012, 13, 1493–1506. [Google Scholar] [CrossRef]
  54. Roy, A.; Royer, A.; Montpetit, B.; Bartlett, P.A.; Langlois, A. Snow specific surface area simulation using the one-layer snow model in the Canadian LAnd Surface Scheme (CLASS). Cryosphere 2013, 7, 961–975. [Google Scholar] [CrossRef]
Figure 1. Data from four NoSREx campaigns in 2009–2010 (a); 2010–2011 (b); 2011–2012 (c); 2012–2013 (d). Panels from top to bottom: SnowScat channel differences for backscattering measured between 10.2, 13.3 and 16.7 GHz, VV-pol, 50° incidence angle; SodRad brightness temperature channel differences between 10.65, 18.7 and 37 GHz channels, V-pol, 50° incidence angle; air temperature (Tair), ground temperature (Tgnd) and ground permittivity (εgnd) at 2 cm depth; snow depth and snow density with 3rd order fit to manual density measurements shown adapted from [27].
Figure 1. Data from four NoSREx campaigns in 2009–2010 (a); 2010–2011 (b); 2011–2012 (c); 2012–2013 (d). Panels from top to bottom: SnowScat channel differences for backscattering measured between 10.2, 13.3 and 16.7 GHz, VV-pol, 50° incidence angle; SodRad brightness temperature channel differences between 10.65, 18.7 and 37 GHz channels, V-pol, 50° incidence angle; air temperature (Tair), ground temperature (Tgnd) and ground permittivity (εgnd) at 2 cm depth; snow depth and snow density with 3rd order fit to manual density measurements shown adapted from [27].
Remotesensing 10 00170 g001aRemotesensing 10 00170 g001b
Figure 2. Time series of effective correlation length pex,eff retrieved for NoSREx campaigns in 2009–2010 (a); 2010–2011 (b); 2011–2012 (c); 2012–2013 (d). Retrieved active microwave pactiveex,eff for channel differences 13.3–10.2, 16.7–10.2, and 16.7–13.3 GHz (VV-polarization). Retrieved passive microwave ppassiveex,eff using the channel differences of 10.65–18.7 GHz, 10.65–37 GHz, and 18.7–37 GHz (V-polarization). In situ measured pex, estimated from either SMP (black) or IceCube SSA measurements (blue), averaged over whole snowpack (see [36]). Error bars reflect standard deviation of all measured samples.
Figure 2. Time series of effective correlation length pex,eff retrieved for NoSREx campaigns in 2009–2010 (a); 2010–2011 (b); 2011–2012 (c); 2012–2013 (d). Retrieved active microwave pactiveex,eff for channel differences 13.3–10.2, 16.7–10.2, and 16.7–13.3 GHz (VV-polarization). Retrieved passive microwave ppassiveex,eff using the channel differences of 10.65–18.7 GHz, 10.65–37 GHz, and 18.7–37 GHz (V-polarization). In situ measured pex, estimated from either SMP (black) or IceCube SSA measurements (blue), averaged over whole snowpack (see [36]). Error bars reflect standard deviation of all measured samples.
Remotesensing 10 00170 g002aRemotesensing 10 00170 g002b
Figure 3. Daily ppassiveex,eff retrieved using channel difference of 18.7–37 GHz, V-pol, compared to pactiveex,eff retrieved using channel difference of 16.7–10.2 GHz, VV pol. Constant values (Table 1) used for ancillary data. Color codes represent the four different seasons of NoSREx (1 to 4).
Figure 3. Daily ppassiveex,eff retrieved using channel difference of 18.7–37 GHz, V-pol, compared to pactiveex,eff retrieved using channel difference of 16.7–10.2 GHz, VV pol. Constant values (Table 1) used for ancillary data. Color codes represent the four different seasons of NoSREx (1 to 4).
Remotesensing 10 00170 g003
Figure 4. Retrieved SWE from 13.3 and 16.7 GHz, VV pol channels of SnowScat for NoSREx campaigns in 2009–2010 (a,b); 2010–2011 (c,d); 2011–2012 (e,f); 2012–2013 (g,h). Left panels: retrievals with Configurations 1 and 3; right panels: retrievals with Configurations 2 and 4.
Figure 4. Retrieved SWE from 13.3 and 16.7 GHz, VV pol channels of SnowScat for NoSREx campaigns in 2009–2010 (a,b); 2010–2011 (c,d); 2011–2012 (e,f); 2012–2013 (g,h). Left panels: retrievals with Configurations 1 and 3; right panels: retrievals with Configurations 2 and 4.
Remotesensing 10 00170 g004aRemotesensing 10 00170 g004b
Figure 5. Comparison of retrieved SWE against manual measurements for all four NoSREx seasons. SnowScat channels 10.2 (a) 13.3 GHz (b) and 16.7 GHz (c). Red: retrieval using Configuration 2 with pex,eff = <pactiveex,eff> (seasonal optimization). Green: retrieval using Configuration 4 with pex,eff = β*ppassiveex,eff. The diagonal lines represent a 1:1 fit +/− 30 mm SWE; an RMSE of 30 mm was defined as a threshold for retrieval accuracy.
Figure 5. Comparison of retrieved SWE against manual measurements for all four NoSREx seasons. SnowScat channels 10.2 (a) 13.3 GHz (b) and 16.7 GHz (c). Red: retrieval using Configuration 2 with pex,eff = <pactiveex,eff> (seasonal optimization). Green: retrieval using Configuration 4 with pex,eff = β*ppassiveex,eff. The diagonal lines represent a 1:1 fit +/− 30 mm SWE; an RMSE of 30 mm was defined as a threshold for retrieval accuracy.
Remotesensing 10 00170 g005
Figure 6. Same as Figure 5 but for channel combinations of 10.2 and 13.3 GHz (a); 10.2 and 16.7 GHz (b); 13.3 and 16.7 GHz (c).
Figure 6. Same as Figure 5 but for channel combinations of 10.2 and 13.3 GHz (a); 10.2 and 16.7 GHz (b); 13.3 and 16.7 GHz (c).
Remotesensing 10 00170 g006
Figure 7. SWE retrieved with channel difference of 16.7–13.3 GHz for NoSREx season 4.
Figure 7. SWE retrieved with channel difference of 16.7–13.3 GHz for NoSREx season 4.
Remotesensing 10 00170 g007
Figure 8. SWE retrieved with channel difference of 10.65–37 GHz using Configurations 1 and 2 adapted for SodRad passive microwave observations from all four NoSREx seasons.
Figure 8. SWE retrieved with channel difference of 10.65–37 GHz using Configurations 1 and 2 adapted for SodRad passive microwave observations from all four NoSREx seasons.
Remotesensing 10 00170 g008
Table 1. Ancillary input parameters and MEMLS3&a model settings.
Table 1. Ancillary input parameters and MEMLS3&a model settings.
ParameterDescriptionValue
FCentre frequencyactive: 10.2, 13.3, 16.7 GHz
passive: 18.7, 37 GHz
θkIncidence angle50°
s0vV-pol reflectivity of snow-ground interface0.03(@37 GHz)…0.06(@10.2 GHz) using εgnd = 4 and hrms = 1 cm [31]
s0hH-pol reflectivity of snow-ground interface0.05(@37 GHz)…0.07(@10.2 GHz) using εgnd = 4 and hrms = 1 cm [31]
ss0v, ss0hSpecular part of reflectivity0.9
qFraction of cross-polarized scattering(cross polarization not used)
TskyDownwelling sky brightness temperature5 (@10.2 GHz) to 35 K (@37 GHz)
[32]
TgndGround temperature−5 °C
TsnowSnow temperature−5 °C
ρsnowSnow density200 kg m−3
Scattering modelIBA
Number of layers1
Table 2. Coefficient of determination between daily averages of pactiveex,eff and ppassiveex,eff retrieved using different channels and channel combinations. All four seasons of retrievals from NoSREx summarized. Statistically insignificant R2 values indicated with asterisks (*). Results from retrievals with measured ancillary data given in parenthesis.
Table 2. Coefficient of determination between daily averages of pactiveex,eff and ppassiveex,eff retrieved using different channels and channel combinations. All four seasons of retrievals from NoSREx summarized. Statistically insignificant R2 values indicated with asterisks (*). Results from retrievals with measured ancillary data given in parenthesis.
R210.2 GHz13.3 GHz16.7 GHz10.2 & 13.3 GHz10.2 & 16.7 GHz13.3 & 16.7 GHz13.3–10.2 GHz16.7–10.2 GHz16.7–13.3 GHz
10.65 GHz0.07 (0.11)0.13 (0.07)0.00 (0.20)0.12 (0.08)0.00 (0.21)0.00 (0.20)0.15 (0.01)0.01 (0.15)0.03 * (0.17)
18.7 GHz0.00 (0.08)0.10 (0.39)0.04 (0.02)0.07 (0.35)0.04 (0.03)0.02 (0.06) *0.27 (0.53)0.04 (0.01)0.10 (0.01)
37 GHz0.00 (0.09)0.40 (0.47)0.74 (0.79)0.33 (0.41)0.73 (0.79)0.75 (0.80)0.59 (0.65)0.74 (0.79)0.59 (0.63)
10.65 & 18.7 GHz0.01 (0.09)0.13 * (0.41)0.03 (0.03)0.10 (0.37)0.03 * (0.03)0.01 * (0.06) *0.29 (0.53)0.04 (0.02)0.10 (0.01)
10.65 & 37 GHz0.00 (0.04)0.40 (0.47)0.73 (0.79)0.33 (0.41)0.73 (0.78)0.75 (0.80)0.59 (0.65)0.74 (0.79)0.58 (0.63)
18.7 & 37 GHz0.00 (0.04)0.42 (0.49)0.72 (0.77)0.34 (0.42)0.72 (0.77)0.74 (0.79)0.62 (0.68)0.73 (0.78)0.56 (0.61)
10.65–18.7 GHz0.03 (0.05)0.27 (0.10)0.19 (0.07)0.19 (0.05)0.19 (0.06)0.23 (0.08)0.62 (0.40)0.23 (0.09)0.09 (0.02)
10.65–37 GHz0.01 (0.01)0.34 (0.38)0.74 (0.79)0.27 (0.32)0.74 (0.78)0.75 (0.78)0.59 (0.60)0.78 (0.82)0.62 (0.68)
18.7–37 GHz0.00 (0.02)0.29 (0.34)0.78 (0.82)0.24 (0.29)0.78 (0.82)0.77 (0.80)0.44 (0.47)0.79 (0.84)0.69 (0.76)
Table 3. The average of pactiveex,eff (in mm) across seasons for different channels and channel combinations.
Table 3. The average of pactiveex,eff (in mm) across seasons for different channels and channel combinations.
Season10.2 GHz13.3 GHz 16.7 GHz 10.2 & 13.3 GHz10.2 & 16.7 GHz 13.3 & 16.7 GHz 13.3–10.2 GHz16.7–10.2 GHz16.7–13.3 GHz
NoSREx 10.280.270.290.270.290.290.260.290.31
NoSREx 20.210.260.310.250.310.300.280.330.36
NoSREx 30.230.210.210.210.210.210.190.210.22
NoSREx 40.230.250.240.250.240.250.270.250.23
Table 4. The average of ppassiveex,eff (in mm) across seasons for different channels and channel combinations. Retrievals at 10.65 and 18.7 GHz or with the channel difference of these frequencies for NoSREx 3 did not produce realistic values of ppassiveex,eff and were discarded.
Table 4. The average of ppassiveex,eff (in mm) across seasons for different channels and channel combinations. Retrievals at 10.65 and 18.7 GHz or with the channel difference of these frequencies for NoSREx 3 did not produce realistic values of ppassiveex,eff and were discarded.
Season10.65 GHz18.7 GHz37 GHz10.65 & 18.7 GHz10.65 & 37 GHz18.7 & 37 GHz10.65–18.7 GHz10.65–37 GHz18.7–37 GHz
NoSREx 10.250.160.230.170.230.230.180.240.25
NoSREx 20.310.180.240.180.240.240.260.260.26
NoSREx 3--0.17-0.170.160.110.170.18
NoSREx 40.340.220.220.220.220.220.270.230.22
Table 5. Bias, uRMSE and R2 of retrieved SWE against in situ observations using retrieval Configurations 1 and 2. Error metrics shown for NoSREx individual seasons as well as for the whole dataset.
Table 5. Bias, uRMSE and R2 of retrieved SWE against in situ observations using retrieval Configurations 1 and 2. Error metrics shown for NoSREx individual seasons as well as for the whole dataset.
Configuration 1NoSREx 1NoSREx 2NoSREx 3NoSREx 4All
FrequencyBias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2
10.2 GHz68630.01−22370.13−26420.40−12490.0218680.05
13.3 GHz40350.4414140.89−72320.1914480.0316480.18
16.7 GHz53320.7078450.91−80190.86−31320.5927640.19
10.2 & 13.3 GHz43370.389170.86−67330.2410480.0315480.16
10.2 & 16.7 GHz51310.6971380.91−78190.86−31320.5725610.19
13.3 & 16.7 GHz47290.6760290.92−78200.84−24340.4822550.20
13.3–10.2 GHz31290.6141150.92−92290.0036510.0221510.16
16.7–10.2 GHz56400.7472440.70−87160.87−34290.7024670.19
16.7–13.3 GHz53420.3693640.71−85130.81−59270.7418800.02
Configuration 2NoSREx 1NoSREx 2NoSREx 3NoSREx 4All
FrequencyBias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2
10.2 GHz−6510.016390.13−12440.400500.02−2480.09
13.3 GHz2340.443170.88−11370.190470.030350.45
16.7 GHz3250.707130.91−3110.86−2290.592230.77
10.2 & 13.3 GHz0360.383180.86−11380.240470.030360.42
10.2 & 16.7 GHz3260.686120.91−4110.85−2300.572230.76
13.3 & 16.7 GHz2260.675110.92−6140.84−2330.481240.74
13.3–10.2 GHz9290.614120.92−11330.001480.024320.54
16.7–10.2 GHz8250.7410160.911130.87−2250.706220.78
16.7–13.3 GHz16380.7317280.8713320.802310.7013340.75
Table 6. Same as Table 5 but for Configurations 3 and 4.
Table 6. Same as Table 5 but for Configurations 3 and 4.
Configuration 3NoSREx 1NoSREx 2NoSREx 3NoSREx 4All
FrequencyBias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2
10.2 GHz41430.29−34400.07111380.3414460.1221590.29
13.3 GHz41260.716240.6561490.2486570.1543460.49
16.7 GHz81410.5759280.81108210.6159350.6273380.69
10.2 & 13.3 GHz41270.670250.6170480.2775540.1640460.48
10.2 & 16.7 GHz79400.5668380.84109200.5857350.6174390.68
13.3 & 16.7 GHz70360.5858350.81104190.5665390.5269370.70
13.3–10.2 GHz41290.7836330.7532560.26112650.0551510.44
16.7–10.2 GHz82430.3491450.82110300.5667360.6385430.61
16.7–13.3 GHz113520.5387360.67100420.0135330.8784550.43
Configuration 4NoSREx 1NoSREx 2NoSREx 3NoSREx 4All
FrequencyBias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2Bias (mm)uRMSE (mm)R2
10.2 GHz−12380.31−51380.0829430.38−32420.12−23450.28
13.3 GHz−14250.73−25250.621570.4312450.16−10370.44
16.7 GHz23320.7732340.8823270.20−8280.6219340.67
10.2 & 13.3 GHz−14270.69−28260.588610.455440.16−11380.41
10.2 & 16.7 GHz21310.7726280.8824280.13−9280.6116320.68
13.3 & 16.7 GHz13260.7716220.8718320.02−4310.5211280.71
13.3–10.2 GHz−15210.78−7190.76−32430.2546550.15−2420.41
16.7–10.2 GHz19270.6447410.8416190.60−1260.6721350.57
16.7–13.3 GHz41380.4759380.8955390.73−25160.8831470.45

Share and Cite

MDPI and ACS Style

Lemmetyinen, J.; Derksen, C.; Rott, H.; Macelloni, G.; King, J.; Schneebeli, M.; Wiesmann, A.; Leppänen, L.; Kontu, A.; Pulliainen, J. Retrieval of Effective Correlation Length and Snow Water Equivalent from Radar and Passive Microwave Measurements. Remote Sens. 2018, 10, 170. https://doi.org/10.3390/rs10020170

AMA Style

Lemmetyinen J, Derksen C, Rott H, Macelloni G, King J, Schneebeli M, Wiesmann A, Leppänen L, Kontu A, Pulliainen J. Retrieval of Effective Correlation Length and Snow Water Equivalent from Radar and Passive Microwave Measurements. Remote Sensing. 2018; 10(2):170. https://doi.org/10.3390/rs10020170

Chicago/Turabian Style

Lemmetyinen, Juha, Chris Derksen, Helmut Rott, Giovanni Macelloni, Josh King, Martin Schneebeli, Andreas Wiesmann, Leena Leppänen, Anna Kontu, and Jouni Pulliainen. 2018. "Retrieval of Effective Correlation Length and Snow Water Equivalent from Radar and Passive Microwave Measurements" Remote Sensing 10, no. 2: 170. https://doi.org/10.3390/rs10020170

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop