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

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 Remote Sens. 2018, 10, 170; doi:10.3390/rs10020170 www.mdpi.com/journal/remotesensing Remote Sens. 2018, 10, 170 2 of 28 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.


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.

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].

Retrieval of Effective Correlation Length
MEMLS3&a was applied in a one-layer configuration to retrieve effective exponential correlation length p ex,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 where SD is the snow depth, p ex,eff the effective (exponential) correlation length, P the number of observation channels k, Φ s k the forward model estimate for sensor types (active or passive), x n the model parameters excluding SD and p ex,eff , φ s k the observed microwave response and σ s k 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 (h rms ) 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 p ex,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 (T air ), ground temperature (T gnd ), 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 p active ex,eff and p passive ex,eff , respectively.The ancillary parameters x, as well as other setting of MEMLS3&a, are summarized in Table 1.

Scattering model IBA Number of layers 1
To analyze the effect of uncertainties of the approach, retrievals of p active ex,eff and p passive ex,eff were repeated using measured values for available ancillary parameters for T air , T gnd , ε gnd and ρ snow .
T air , T gnd 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.

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 where SWE is the sole free parameter and p ex,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 p ex,eff in the retrieval: • Configuration 1: An overall average of optimized daily values <p active ex,eff > was calculated from all retrievals of p active ex,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 (p ex,eff (t) = <p active ex,eff >).

•
Configuration 2: As in Configuration 1, but the average of optimized daily values <p active ex,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 (p ex,eff (t) = p passive ex,eff (t)).As a default, p passive ex,eff was obtained from 18.7-37 GHz, V-pol, radiometer retrievals.

•
Configuration 4: As in Configuration 3, but the value of p passive ex,eff was scaled so that p ex,eff (t) = β p passive ex,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 p ex,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 p active ex,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 p active ex,eff and p passive ex,eff .All other input variables of the MEMLS3&a model were kept constant in the retrieval, following Table 1, as in the retrieval of p ex,eff .
As with retrievals of p ex,eff (see Section 2.2), SWE retrievals were repeated using measured ancillary data for T air , T gnd , ε 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 p passive ex,eff to initialize the retrievals.

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 p active ex,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.

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.

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.2GHz, 16.7-10.2GHz, and 16.7-13.3GHz 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 (T air ), as well as the 2 cm depth ground temperature (T gnd ) 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.2GHz), 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.

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 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 (T air ), ground temperature (T gnd ) 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].

Retrieved Active and Passive Correlation Length
The SnowScat and SodRad measurements were applied in Equation ( 1) to retrieve an effective correlation length p ex,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, p ex,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 p active ex,eff and p passive ex,eff during the four seasons of NoSREx.Retrievals of p active ex,eff using the three combinations of channel differences between the 10.2, 13.3 and 16.7 GHz are shown, as well as retrieved p passive ex,eff using the channel difference combinations between the 10.65, 18.7 and 37 GHz channels.The three time series of p active ex,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.2GHz produces a less dynamic time series of p active ex,eff .With the exception of 2012-2013, as well as February 2012, the channel difference of 16.7-13.3GHz gives the largest value of p active ex,eff , while the channel difference of 13.3-10.2GHz results in the smallest values.Interestingly, for the last season these relations are reversed.The passive-microwave derived p passive ex,eff follows some of the temporal behavior of the three p active ex,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 p active ex,eff and p passive ex,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 p ex 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 D 0 , so that [25] The following empirical linear relations between D 0 and p ex have been derived [29] p ex = 0.4•D 0 for dendritic grains (4) 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 p ex measured using IceCube for the third seasons shows relatively good agreement with p passive ex,eff (significant R 2 values for the channel differences depicted in Figure 2c ranged from 0.48 to 0.61), while p active ex,eff generally overestimates measured values and showed a lower coefficient of determination (significant R 2 values ranged from 0.37 to 0.46).Nevertheless, both p passive ex,eff and p active ex,eff as well as in situ measurements show decreasing values of p ex in mid-February; this can be explained by snow accumulation in that period (Figure 1c), which decreased the measured average p ex .For the fourth season, both p active ex,eff and p passive ex,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 p ex is closely replicated by the retrievals, with a generally decreasing trend in p ex , p active ex,eff and p passive ex,eff until February 2013.Significant R 2 values against in situ p ex for the fourth season ranged from 0.30 to 0.61 and 0.27 to 0.34, for p active ex,eff and p passive ex,eff respectively.Measurements of p ex using IceCube generally agree well with SMP measurements of p ex 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 p ex, 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 p ex, eff with in situ measurements.3 GHz (VV-polarization).Retrieved passive microwave p passive ex,eff using the channel differences of 10.65-18.7 GHz, 10.65-37 GHz, and 18.7-37 GHz (Vpolarization).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.
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 p active ex,eff to p passive ex,eff were analyzed by calculating basic statistics between parameters retrieved using different channel configurations.Table 2 summarizes the coefficient of determination of p active ex,eff relative to p passive ex,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 p passive ex,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 p active ex,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 R 2 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 p active ex,eff and p passive ex,eff in particular for the third and fourth season of NoSREx, with R 2 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 ex,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 p ex , 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.
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 p active ex,eff to p passive ex,eff were analyzed by calculating basic statistics between parameters retrieved using different channel configurations.Table 2 summarizes the coefficient of determination of p active ex,eff relative to p passive ex,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 (p ex,eff < 0.05 mm or p ex,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 p passive ex,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 p active ex,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 R 2 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 p active ex,eff and p passive ex,eff in particular for the third and fourth season of NoSREx, with R 2 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 p active ex,eff with p active ex,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 p passive ex,eff against p active ex,eff separately for the four seasons of NoSREx, using the channel differences of 16.7-10.2GHz (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.06mm for p passive ex,eff relative to p active ex,eff was apparent.Against retrievals with 18.7-37 GHz, p active ex,eff was on average 15% larger.This value was applied in SWE retrieval Configuration 4 (see Section 2.3) to scale between p passive ex,eff and p active ex,eff .It is notable in Figure 3 that retrievals from both sensor types similarly resolve the differing magnitudes of p ex,eff between the seasons; the second season gave on average the largest values, while the third season indicated the lowest values of p ex,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 p ex,eff , also exhibits the largest scatter and lowest correlation coefficient between active and passive retrievals.
The average of retrieved p ex,eff for the four NoSREx seasons are summarized in Tables 3 and 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 p ex,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 p active ex,eff on average by 15%, 7% and 5% for 10.2, 13.3 and 16.7 GHz, respectively.For p passive ex,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 p active ex,eff , and a change of −3 to 16% for p passive ex,eff ).In most cases, the use of ancillary data increased the absolute difference between p active ex,eff and p passive ex,eff ; the average difference between all retrievals increased from 22 to 40%.The difference between p active ex,eff and p passive ex,eff retrieved with the 18.7-37 GHz channel difference increased from 15 to 30%, which was used to scale p passive ex,eff in Configuration 4, when performing comparative SWE retrievals with measured ancillary in situ information.Figure 3. Daily p passive ex,eff retrieved using channel difference of 18.7-37 GHz, V-pol, compared to p active ex,eff retrieved using channel difference of 16.7-10.2GHz, VV pol.Constant values (Table 1) used for ancillary data.Color codes represent the four different seasons of NoSREx (1 to 4).

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 <p active ex,eff>, respectively.In the third configuration, pex,eff was taken from the temporally closest retrieval of p passive ex,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 = β*p passive ex,eff.A constant value of β = 1.2 was applied over all seasons, based on the average difference between p passive ex,eff and p active ex,eff.
Similarly to the retrieval of pex,eff, SWE retrievals were performed for several channels and channel combinations.In Configurations 1 and 2, <p active ex,eff> was always chosen from the corresponding retrieval at the same channel or combination of radar channels, while in Configurations 3 and 4, p passive ex,eff was always obtained from the retrieval using the 18.7-37 GHz, Vpol channel difference (retrievals from H-pol were not applied).Retrievals using this passive ex,eff retrieved using channel difference of 18.7-37 GHz, V-pol, compared to p active ex,eff retrieved using channel difference of 16.7-10.2GHz, VV pol.Constant values (Table 1) used for ancillary data.Color codes represent the four different seasons of NoSREx (1 to 4).

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, p ex,eff in Equation ( 2) was set to be the overall average, or the seasonal average, of retrieved <p active ex,eff >, respectively.In the third configuration, p ex,eff was taken from the temporally closest retrieval of p passive ex,eff , thus testing the suitability of passive-microwave derived p ex,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 p ex,eff = β*p passive ex,eff .A constant value of β = 1.2 was applied over all seasons, based on the average difference between p passive ex,eff and p active ex,eff .Similarly to the retrieval of p ex,eff , SWE retrievals were performed for several channels and channel combinations.In Configurations 1 and 2, <p active ex,eff > was always chosen from the corresponding retrieval at the same channel or combination of radar channels, while in Configurations 3 and 4, p passive ex,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 p active ex,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 <p active ex,eff >; green lines represent retrievals using p passive ex,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 < p active ex,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 p passive ex,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 p passive ex,eff for all frequency combinations of p active ex,eff ; a more optimized fit could be achieved by defining β individually for each channel or combination of channels, as well as for each season. Remote Sens. 2018, 10, x FOR PEER REVIEW 15 of 28 microwave channel combination showed the best overall agreement with p active ex,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 <p active ex,eff>; green lines represent retrievals using p passive ex,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 < p active ex,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 p passive ex,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 p passive ex,eff for all frequency combinations of p active ex,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<p active ex,eff> (Configuration 2), while green markers denote retrievals with a temporally dynamic but scaled p passive ex,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 = <p active ex,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.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 <p active ex,eff > (Configuration 2), while green markers denote retrievals with a temporally dynamic but scaled p passive ex,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 p ex,eff = <p active ex,eff > yields better error metrics compared to Configuration 4. The good performance of Configuration 2 can be expected due to the seasonal optimization of p ex,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 p ex,eff (Configuration 1) results in bias errors for several frequency combinations, particularly notable for the third season.R 2 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 p ex,eff and the necessity for seasonal optimization.All frequency combinations show both reduced uRMSE and an increased R 2 value when using Configuration 2, compared to Configuration 1. Furthermore, all retrievals with Configuration 2 using the 16.7 GHz channel give an R 2 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 R 2 .It can be also observed that, in this case, the addition of a secondary channel to 16.7 provides an improvement in R 2 and uRMSE.
Table 6 summarizes similarly the retrievals using Configurations 3 and 4, i.e., by applying p ex,eff from passive microwave observations in the retrieval.Relatively large bias errors are obtained when using p passive ex,eff directly without scaling (Configuration 3); however, R 2 values in the order of 0.6 are still obtained for several frequency combinations on average, with R 2 values exceeding 0.8 for several frequency combinations for the second season, and the 16.7-13.3GHz combination for the fourth season.This provides an indication that the passive microwave derived p ex,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 p passive ex,eff (Configuration 4) further improves the error metrics against in situ SWE, in particular in terms of reduced bias and uRMSE.R 2 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.2GHz 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.3GHz.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 (T air , T gnd , ε gnd , and ρ snow ).Values of p ex,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 p passive ex (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 p passive ex 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 p active ex (Configuration 2 in Table 5), with the exception of the larger bias error of the 13.3-16.7 GHz channel difference retrievals.
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.3GHz.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 p passive ex (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 p passive ex 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 p active ex (Configuration 2 in Table 5), with the exception of the larger bias error of the 13.3-16.7 GHz channel difference retrievals.

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 p passive ex,eff and p active ex,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

Discussion
Results in Section 4.1 demonstrated that at certain frequencies and combinations of frequencies, active and passive -derived p ex,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 p ex,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 p passive ex,eff and p active ex,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 p active ex,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 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 p passive ex 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]).

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 p ex,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 p ex , and the use of a one-layer model for the snowpack to retrieve p ex,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 (R 2 values between 0.59 and 0.91 for different seasons), or in combination with lower frequencies (R 2 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, R 2 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.

Figure 3 .
Figure 3. Daily p passiveex,eff retrieved using channel difference of 18.7-37 GHz, V-pol, compared to p active ex,eff retrieved using channel difference of 16.7-10.2GHz, VV pol.Constant values (Table1) used for ancillary data.Color codes represent the four different seasons of NoSREx (1 to 4).

Figure 5 .
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 p ex,eff = <p active ex,eff > (seasonal optimization).Green: retrieval using Configuration 4 with p ex,eff = β*p passive ex,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 8 .
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 .
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.

Table 1 .
Ancillary input parameters and MEMLS3&a model settings.

Table 2 .
Coefficient of determination between daily averages of p active ex,eff and p passive ex,eff retrieved using different channels and channel combinations.All four seasons of retrievals from NoSREx summarized.Statistically insignificant R 2 values indicated with asterisks (*).Results from retrievals with measured ancillary data given in parenthesis.

Table 3 .
The average of p active ex,eff (in mm) across seasons for different channels and channel combinations.

Table 4 .
The average of p passive ex,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 p passive ex,eff and were discarded.

Table 3 .
The average of p active ex,eff (in mm) across seasons for different channels and channel combinations.

Table 4 .
The average of p passive ex,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 p passive ex,eff and were discarded.

Table 5 .
Bias, uRMSE and R 2 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 6 .
Same as Table 5 but for Configurations 3 and 4.