Next Article in Journal
Evaluating the Applicability of Four Latest Satellite–Gauge Combined Precipitation Estimates for Extreme Precipitation and Streamflow Predictions over the Upper Yellow River Basins in China
Next Article in Special Issue
Estimating Snow Depth Using Multi-Source Data Fusion Based on the D-InSAR Method and 3DVAR Fusion Algorithm
Previous Article in Journal
Sparse Weighted Constrained Energy Minimization for Accurate Remote Sensing Image Target Detection
Previous Article in Special Issue
Mapping Radar Glacier Zones and Dry Snow Line in the Antarctic Peninsula Using Sentinel-1 Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Davos-Laret Remote Sensing Field Laboratory: 2016/2017 Winter Season L-Band Measurements Data-Processing and Analysis

1
Swiss Federal Research Institute WSL, CH-8903 Birmensdorf, Switzerland
2
Gamma Remote Sensing AG, CH-3073 Gümligen, Switzerland
*
Author to whom correspondence should be addressed.
Remote Sens. 2017, 9(11), 1185; https://doi.org/10.3390/rs9111185
Submission received: 16 October 2017 / Revised: 9 November 2017 / Accepted: 15 November 2017 / Published: 21 November 2017
(This article belongs to the Special Issue Snow Remote Sensing)

Abstract

:
The L-band radiometry data and in-situ ground and snow measurements performed during the 2016/2017 winter campaign at the Davos-Laret remote sensing field laboratory are presented and discussed. An improved version of the procedure for the computation of L-band brightness temperatures from ELBARA radiometer raw data is introduced. This procedure includes a thorough explanation of the calibration and filtering including a refined radio frequency interference (RFI) mitigation approach. This new mitigation approach not only performs better than conventional “normality” tests (kurtosis and skewness) but also allows for the quantification of measurement uncertainty introduced by non-thermal noise contributions. The brightness temperatures of natural snow covered areas and areas with a reflector beneath the snow are simulated for varying amounts of snow liquid water content distributed across the snow profile. Both measured and simulated brightness temperatures emanating from natural snow covered areas and areas with a reflector beneath the snow reveal noticeable sensitivity with respect to snow liquid water. This indicates the possibility of estimating snow liquid water using L-band radiometry. It is also shown that distinct daily increases in brightness temperatures measured over the areas with the reflector placed on the ground indicate the onset of the snow melting season, also known as “early-spring snow”.

1. Introduction

Microwave remote sensing is preeminently suitable to achieve quantitative estimates of large-scale terrestrial state parameters. This is partially due to the high transparency of the atmosphere regardless of weather conditions, particularly within the low-frequency bands of the microwave regime. Furthermore, remote microwave observations do not rely on sunlight, which is decisive to achieve continuous information, especially over the Cryosphere at high latitudes. Another key advantage of microwave remote sensing, over, for example, optical remote sensing, is the significantly larger penetration depth of microwaves in media such as snow [1], vegetation, and soils [2] with moderate amounts of liquid water. However, to retrieve information on terrestrial state-parameters from remotely-measured brightness temperatures, proper microwave emission models are needed, which relate such state parameter values to microwave brightness temperatures. Typically, retrieval approaches involve computationally-expensive iterative procedures employed to optimally fit modeled microwave signatures to the corresponding measurements. Thus, it is important to develop radiative transfer models which not only capture the observed scene’s inherent radiative processes as accurately as possible, but also keep the number of model parameters as low as possible and limited explicitly to quantities which are accessible from other sources, such as from models or auxiliary remote sensing information.
Nevertheless, the development of radiative transfer models with the aforementioned features necessitates the development and operation of field laboratories allowing for comprehensive, synchronous and collocated in-situ and near-distance remote sensing measurements. Such field laboratories, dedicated to the development of methods to optimally exploit microwave remote sensing of the Cryosphere, are still rare today. Examples of such currently operative research sites include: (i) the Pallas-Sodankylä supersite in Finland [3,4] operated by the Finnish Meteorological Institute Arctic Research Centre (FMI-ARC); (ii) the DOMEX site at Dome-C, Antarctica [5]; (iii) the Kernen Crop Research Farm (KCRF) [6] in Saskatoon operated by the University of Saskatchewan, and (iv) the Maqu regional SMST monitoring network [7]. The most notable example in the Alps is the snow laboratory built on the Weissfluhjoch (2540 m a.s.l., Switzerland) in 1936. Data from this field laboratory include daily records of snow depths and snowfall dating back more than 70 years. During the 1980s, a comprehensive long-term program of microwave-signature studies (summarized in [8]) was conducted at the Weissfluhjoch snow laboratory. This research laid the foundation for the “microwave emission model of layered snowpacks” (MEMLS) [9,10], which is one of today’s widely-used tools to simulate thermal microwave emission of snow [11,12,13,14,15,16].
Nearly three decades later, the Swiss Federal Research Institute WSL established the Davos-Laret remote sensing field laboratory in 2016 as the only currently-operating Alpine test site dedicated to the development of novel retrieval approaches for the estimation of snow properties from microwave remote sensing data. The present paper provides details on this field laboratory and presents the in-situ data collected during its first operation in Winter 2016/2017 in Section 2 and Section 3, respectively. The content is then narrowed to what is relevant for L-band radiometry in the retrieval of snow and ground parameters. In Section 4, we outline the steps taken to achieve calibrated L-band brightness temperatures from ELBARA-II [17] radiometer raw data. This section includes novel approaches to mitigate and filter non-thermal disturbances in ELBARA-II measurements, as well as a method to optimally compensate for thermal noise imposed by transmission losses. With the main purpose being its use for ground-based measurements and validation campaigns, ELBARA-II is designed to have very similar technical specifications to the on-board radiometer of the ESA’s SMOS satellite [18]. Thus, ELBARA-II (and subsequent generations) is used by many research institutes around the world (see for example [7,15,19,20,21]), indicating that refined data processing methods are beneficial to the passive L-band remote sensing scientific community.
It is worth noting that the penetration depth of L-band microwaves drops from >300 m in dry snow to only a few centimeters for wet snow with ~5% snow liquid water content [22,23]. This, consequently, results in immediate retrieval failure of snow-column properties, or even more so of the underlying ground, as shown in [15,24,25]. Accordingly, as one of the main research questions in the Davos-Laret 2016/2017 winter campaign, the impact of snow liquid water content on L-band brightness temperatures is investigated. Section 5 explains the emission model (Section 5.1) used to analyze the response of simulated brightness temperatures to snow liquid water for snowpacks over “natural ground” and over areas covered with a reflecting metal grid (“reflector areas”—Section 5.2). Afterwards, the temporal and spatial variations of the calibrated L-band brightness temperatures measured over the 2016/2017 winter campaign are presented and discussed. The specific content of Section 4 and Section 5 is, on the one hand, valuable on its own, and, on the other hand, preparative for ongoing developments in the retrieval of snow density and ground permittivity [14,15,16]. In addition, it is of value for the development of a new retrieval approach for the estimation of snow liquid water-column from passive L-band radiometry. The latter two instances are further outlined in Section 6.

2. Test-Site Description

The Davos Laret Remote Sensing Field Laboratory (48°50′53″ N, 9°52′19″ E) in Switzerland is a 50 m × 50 m area with an approximate elevation of 1450 m above sea level. The ground is mostly flat with some smooth slopes on the northwestern side of the site. The valley, including the site area, is surrounded by mountains with an average height difference of ~400 m with respect to the site. The site area is surrounded by Lake Schwarz on the northwestern side, canopy forest on the southeastern side and local buildings on the northeastern and southwestern sides. The spring and summertime vegetation cover of the site is grass.
Dielectric permittivity and temperature of the ground were measured along the two transects shown in Figure 1 using twelve SMT-100 [26] in-situ sensors. Additionally, two SMT-100 sensors were installed on a wooden stick 15- and 50-cm above the ground to measure the permittivity and temperature of the snow at the respective heights (dashed red square in Figure 1). All the sensors were driven with a DT80 data logger which recorded the measured qualities every five minutes. Section 3 presents the time series of these in-situ recordings over the 2016/2017 campaign. In the 2016/2017 Davos Laret campaign, three remote sensing instruments were used:
  • The second generation ELBARA-II of the ETH L-band radiometer (ELBARA) operating at horizontal and vertical polarization within the protected part of the L-band (1400–1427 MHz). ELBARA-II technical specifications and its performance under field conditions are outlined by Schwank et al. [17,27]. The first generation ELBARA was designed and constructed in 2001 at the Institute of Applied Physics, University of Bern [28]. A number of ELBARA-II radiometers were produced by the company Gamma Remote Sensing AG (Gümligen, Switzerland) with the main purpose of calibration and validation activities supporting ESA’s SMOS mission [18]. The ELBARA-II assembly, mounted atop an 8-meter tower, includes a tracking system to allow for automated observations of brightness temperatures T B p ( θ ) at different nadir angles θ and azimuth. For the work presented here, ELBARA-II is the key remote sensing instrument.
  • MORA is a single polarization X-band radiometer operating at 11.4 GHz [29]. It is mounted on top of the ELBARA-II antenna (as in [30]) to measure collocated and synchronous X-band brightness temperatures.
  • The SnowScat scatterometer is a fully polarimetric, coherent stepped-frequency continuous wave radar operating in the frequency range 9–18 GHz [31]. It is mounted on a rail and attached to a 10-m tower, which allows for elevation and azimuth scanning over the site. Using Synthetic Aperture Radar (SAR) techniques, 3D backscatter maps of the site are made indicating parameters such as radiation penetration depth and snow water equivalent of dry snow.

3. In-Situ Measurements

The in-situ measurements performed during the first operation in Winter 2016/2017 of the Davos-Laret field laboratory are presented here, focusing on those used to analyze the collocated and synchronous passive L-band data. Manual in-situ snow-pit characterization (height h S , mass-density ρ S , and near infrared (NIR) photography) are shown in Section 3.1; automated measurements of ground permittivities ε G and temperatures T G , as well as snow temperatures T S , and meteorological data (air temperatures T air , precipitation r) are shown in Section 3.2.

3.1. Snow Density Profiles

The first major snow fall event during the winter of 2016/2017 occurred on 3 January. During the snow covered period, regular snow-pit measurements were performed manually. Snow height ( h S ) and mass-density ( ρ S ) were measured approximately once a week with the first and last measurement on 9 January and 22 March, respectively. A snow cutter was used to measure snow density with a depth resolution of 10 cm. The green crosses in Figure 2a,b show measured h S and ρ S , where the latter represents the average density of the bottom 10 cm of the snowpack derived from the respective profile measurements. Red lines are B-splines fitted to estimate temporal variations of h S and ρ S in between the times of measurements. The reason for showing snow bottom-layer density ρ S is that the snowpack parameter is expected to be the most influential on L-band emission for conditions as explained in [13,14]. Observed h S and ρ S reveal the expected increasing trend throughout the snow season, while short-term variations are mostly explained by snow metamorphism, settlement, and melting.
With each of the quantitative snow-pit measurements, a NIR image of the snowpack profile was taken to qualitatively track the evolution of the snowpack structure. Figure 3a,b shows the NIR images taken from the newly-formed snowpack on 9 January during the “cold winter period” and the same snowpack and location on 27 February during the “early spring period”, respectively. The most distinct structural features that resulted from snowpack evolution associated with changes in, for example, density and grain size, are marked with dashed green lines. The young snowpack (Figure 3a) has a simple structure with a homogeneous layer in the lowest ~20 cm which was persistent throughout the “cold winter period”. The ~5-cm layer above is less dense, and the fine-grained crisp topmost ~5 cm represents the fresh snow fallen the day before the image was taken. Compared to the young snowpack, the mature snowpack (Figure 3b), nearly two months later, shows almost no clear layering structure. The respective increase in complexity is associated with depth hoar, rounded-snow, and faceted snow formed at the beginning of the “early spring period”. Additionally, the preferential vertical flow paths are signs of liquid water infiltration caused by, for instance, rain or wet snowfall.

3.2. Permittivity, Temperature, and Precipitation

Ground permittivities ε G and temperatures T G were measured automatically every 5 minutes using SMT-100 sensors [26] installed approximately 5 cm below the ground’s surface. Additionally, two SMT-100 sensors were installed 15 cm and 50 cm above ground to measure either air or snow temperatures T 15   cm and T 50   cm depending on snow height h S at time of measurement. The thin lines in Figure 4a,b show the time-series of ε G from the respective six sensors installed along transects one and two shown in Figure 1. The thick red line in Figure 4a,b are identical ε G resulting from averaging all twelve ground sensor readings. Likewise, Figure 4c shows T G the average of all ground sensor measurements. Figure 4d shows the temperatures T 15 cm (red), T 50 cm (green), and T air (blue) measured approximately 8 m above the ground by the PT-100 sensor mounted outside the ELBARA-II radiometer. Figure 4e shows precipitation rates r (for 10-min time windows) for the time period 28 November to 15 March, corresponding to the time period shown for the above in-situ measurements.
The first week of the in-situ time series reveals the rapid freezing of the bare ground surface. This, of course, results from the low diurnal heat input to the ground associated with air temperatures mostly below the freezing point. Figure 4c,d indicate that while during afternoons the temperatures above ground still show short-term peaks above 0 °C, ground temperatures steadily decrease until they fall below the freezing point. However, at the latest by the middle of December, ground permittivities measured by all sensors drop to 4 ε G 7 indicating that at least the top 5 cm of the ground is completely frozen. Around 26 December, a slight thawing event takes place as the result of increased time-integrated heat input and precipitation to the still bare ground. This, in turn, is reflected by T G approaching the 0 °C curtain, and finally the increased ε G recorded that indicate increased liquid water in the ground surface.
During the “snow-free period” where the ground is frozen, permittivity readings still show small diurnal variations as the result of partial surface melting caused by increased heat input during cloudless afternoons and refreezing overnight. These variations almost completely disappear by 3 January with the onset of snow due to thermal insulation of the accumulating dry snow. The recorded ground permittivities stay at low levels 4 ε G 7 between mid-December and mid-February, indicating that at least the uppermost 5 cm of the ground is frozen during this approximately two-month period.
As shown in Figure 4a,b, in-situ permittivities start increasing gradually after 31 January, indicating that the ground surface starts thawing with the beginning of the “early spring period”. By the middle of March, the ground has almost entirely thawed. Furthermore, deviations between permittivity readings of the different sensors are greater under not-deeply-frozen ground conditions, indicating that heterogeneity of ground permittivity is significantly higher under such transient ground conditions. The most marked example takes place around 21 February, when the readings of sensor d (black symbols in Figure 4b suddenly increase to ~20, while the other sensors along transect 2 react in a much less distinct manner. This distinct heterogeneity in ground permittivity is most likely caused by surface water runoff ponding locally at the position of sensor d.

4. L-Band Radiometry and Raw Data Processing

Brightness temperatures T B p ( θ ) measured at given nadir angles θ and horizontal (p = H) and vertical (p = V) polarizations are the basis for successful retrievals based on passive L-band remote sensing. However, calibrated T B p ( θ ) are not directly provided by the ELBARA-II radiometer. Instead, they are derived from the instrument’s output data using a radiometer-specific raw-data processing approach. Furthermore, radiometer measurements are sometimes affected by radiative sources other than the natural scene of interest. These perturbing radiation contributions can include man-made Radio Frequency Interferences (RFI) and thermal noise from the surrounding area such as a nearby forest or lake.
Section 4.1 outlines the processing-chain to achieve calibrated T B p ( θ ) from ELBARA-II raw-data outputs. Section 4.2 explains and demonstrates an updated RFI mitigation approach. Section 4.3 addresses the issue of computing effective transmission line losses with respect to their time-variability. Finally, the approach used to derive T B , R p ( θ ) emitted exclusively from the “reflector areas” (see Figure 1) is explained in Section 4.4.

4.1. L-Band Brightness Temperatures

ELBARA-II raw data outputs are voltage samples U i consisting of i = 1,…, nps = 2400 measurements performed at the sampling rate of 800 Hz during the integration time of 3 s [17]. An ELBARA-II measurement cycle includes the sequential acquisition of voltage samples associated with two internal calibration noise sources (the Active Cold Source (ACS) and the Resistive Source (RS)), and associated with the noise entering the horizontal (p = H) and vertical (p = V) polarization port of the radiometer assembly (RMA). Each of these voltage samples are recorded for two 11-MHz channels ch = 1, 2 within the protected part (1400–1427 MHz) of the L-band. Accordingly, an ELBARA-II measurement cycle consists of eight voltage samples U ACS , i c h , U RS ,   i c h , and U RMA ,   i p , c h ultimately used to derive calibrated T B p ( θ ) . The processing steps, illustrated in the flowchart in Figure 5 and explained in Section 4.1.1 and Section 4.1.2, exclusively use mean values of voltage samples. This is henceforth indicated by omitting the sub-index i.

4.1.1. Active Cold Source (ACS) Calibration Using Sky Measurements

The noise temperature T ACS p , c h of the ACS, implemented in the ELBARA-II radiometer, is calibrated using two noise standards. The accurately known Resistive Source (RS) noise temperature T RS = T RMA (= actual physical temperature of the interior of the radiometer assembly (RMA)) is one of the noise sources associated with corresponding raw data (voltages) U RS , i c h where ch = 1, 2 and i = 1,..., nps = 2400. The second standard noise temperatures are T RMA p = T sky p at the RMA input-ports p = H, V associated with raw data U RMA ,   i p , c h = U sky ,   i p , c h when pointing the antenna towards the sky at the nadir angle θ = 140 ° . Noise temperatures T sky p = T RMA p (including inherent noise of the transmission line (TL)) are smaller than ACS noise temperatures T ACS p , c h , which is crucial for the accurate calibration of the ACS. Considering the fact that T ACS p , c h , determined through sky measurements, slightly varies when performed at antenna-port polarization p = H, V, and also in the two frequency channels ch = 1 and 2, ACS noise temperatures T ACS H , 1 , T ACS H , 2 , T ACS V , 1 , and T ACS V , 2 are calculated separately for each of the polarization-channel pairs. This is an improvement over the corresponding original ACS calibration approach described in [17] where T ACS p , c h is assumed to be independent of p and ch. It is seen as a pragmatic solution to compensate for impacts of differences in TL losses along the H and V polarization paths (explained in Section 4.3), and differences between the RMA transfer functions of the two frequency channels (ch = 1 and 2). The rest of this section outlines the procedure used to obtain the four ACS noise temperatures T B p , c h (ch = 1, 2 and p = H, V), which are used to eventually derive calibrated T B p .
Lossy coaxial transmission lines used to connect the H and V antenna port with the corresponding RMA input ports add their own thermal noise. Consequently, reference noise temperatures T sky p = T RMA p at the RMA input p = H, V are larger than the un-polarized L-band brightness T sky 4.5   K entering the aperture of the sky-looking antenna:
T sky p = T RMA p = T sky + ( 1 t TL p ) ( T TL T sky )
Sky brightness temperature T sky entering the aperture of the sky-looking antenna is simulated using the model described in [32,33] for the following parameter values: air humidity q air = 10   g   kg 1 , rain rate r air = 0   mm   h 1 , the elevation of the Davos Laret site h site = 1450   m , the nadir angle θ = 140 ° and T air recorded by ELBARA-II at the time of the sky measurement. For T air = 273.15   K = 0   ° C , this yields T sky 4.4   K . The second summand ( 1 t TL p ) ( T TL T sky ) in Equation (1) is the noise added by the TLs (p = H, V) at common physical temperature T TL = T air . Their power transmission factor t TL p is given by the respective transmission loss L TL p [dB] as:
t TL p = 10 L TL p / 10
Finally, the four ACS noise temperatures T ACS p , c h (p = H, V and ch = 1, 2) are derived from the means U RS c h , U ACS c h , U sky p , c h of the voltage samples U RS , i c h , U ACS , i c h , U RMA , i p , c h = U sky , i p , c h and the two reference noise temperatures T RS and T sky p as:
T ACS p , c h = T RS T sky p U RS c h U sky p , c h ( U ACS c h U sky p , c h ) + T sky p

4.1.2. Antenna Brightness Temperatures

Calibrated antenna brightness temperatures T B p (p = H, V) measured over ground-areas are derived as the mean of T B p , 1 and T B p , 2 at the frequency channels ch = 1 and 2:
T B p = T B p , 1 + T B p , 2 2
The channel specific T B p , c h are expressed by means of noise temperatures T RMA p , c h associated with the mean sample voltage U RMA p , c h measured for the radiometer pointed towards ground. Again, T RMA p , c h at the RMA-input exceed brightness temperatures T B p , c h by the amount of inherent noise added by the TL. Therefore, in analogy with Equation (1), they are related via:
T RMA p , c h = T B p , c h + ( 1 t TL p ) ( T TL T B p , c h )
Solving Equation (5) for T B p , c h immediately yields the channel specific T B p , c h corrected by the TL noise as is used in Equation (4):
T B p , c h = T RMA p , c h ( 1 t TL p ) T TL t TL p
It should be noted that losses L TL p , used to represent the TLs’ inherent noise (second summands in Equations (1) and (5)), can vary in time for a number of reasons such as temperature effects, mechanical stress, and immersive moisture. This can cause both reversible changes in L TL p , as well as irreversible increases in L TL p as the result of ageing. These variations can be taken into account by computing and using the effective L TL p at any given time during the measurement campaign, the method of which is explained in Section 4.3.
Again, in close analogy to Equation (3), the noise temperatures T RMA p , c h used in Equation (5) and ultimately used in Equation (4) to derive the calibrated T B p , are computed as:
T RMA p , c h = T RS T ACS p , c h U RS c h U ACS c h S RMA p , c h ( U RMA p , c h U ACS c h ) + T ACS p , c h
The radiometer sensitivity S RMA p , c h 0.322   K · mV 1 , defined above, is almost independent of p = H, V and ch = 1, 2. It will be used as part of the novel Radio Frequency Interference (RFI) mitigation method introduced in Section 4.2 to quantify distortions Δ T B p , c h of T B p caused by non-thermal noise contributions. In addition, as part of the RFI mitigation method, the raw data sample mean U RMA p , c h used in Equation (7) is replaced with the mean U Gauss p , c h of the Gaussian-model fitted to the probability density function (PDF) of the measured voltage sample U RMA ,   i p , c h .

4.2. Radio Frequency Interference Detection and Mitigation

RFI can be defined as any electromagnetic signal, other than thermal noise emitted from within the antenna footprint, contributing to the measured antenna brightness temperature T B p . Man-made RFI sources include spurious signals from other frequency bands, spread-spectrum signals overlapping the “protected” band of operation, or out-of-band emissions not properly rejected by the pre-detection filters [34]. The statistical distribution of a measured sample U RMA , i p , c h performed on a perfectly undisturbed thermal noise source should theoretically match a Gaussian probability density function PDF Gauss , whereas its standard deviation results from the uncertainties of the individual U RMA , i p , c h determined by the time-bandwidth product (see Section 3.2 in [17]).
Using RFI corrupted T B p in remote sensing results in erroneous retrievals [35,36,37]. Thus, it is critically important to mitigate, or at least detect, RFI in passive microwave remote sensing. RFI detection methods used most often are essentially statistical “normality” tests which quantify the non-Gaussianity of measured sample distributions by means of using thresholds for associated Kurtosis and Skewness [34]. Among the more sophisticated algorithms are the “Kolmogorov-Smirnov” (KS) method [38] and the Lilliefors (L) normality test [39]. However, the former can often wrongly estimate the mean and variance of the normal distribution and the latter has limitations in that it tends to be more sensitive at the center than near the tail of the distribution.
Widely used “normality” tests (based on Kurtosis and Skewness), were initially tried on our ELBARA-II data to detect RFI. However, as also reported in [40,41,42,43], these simple approaches often failed even for physically meaningless measurements T B p , c h . Consequently, we use a more rigorous normality test which not only identifies RFI, but also quantifies the extent to which non-thermal distortion Δ T B p , c h affects T B p , c h . Quantitative estimates of distortions become valuable especially when radiometer measurements of T B p , c h are used in retrieval schemes because it allows to assign less weight to T B p , c h in accordance with its non-thermal distortion Δ T B p , c h .
The RFI mitigation approach used here consists of the following steps. Firstly, the probability density function P D F m ( U RMA , i p , c h ) associated with the measured (m) raw-data voltage sample U RMA , i p , c h is computed numerically. Secondly, the measurement based P D F m ( U RMA , i p , c h ) is approximated with a Gaussian model P D F Gauss ( U RMA , i p , c h ) including the parameters U ^ Gauss p , c h , U Gauss p , c h , σ Gauss p , c h to represent optimized values of its peak, mean, and standard deviation, respectively:
P D F Gauss ( U RMA , i p , c h ) = U ^ Gauss p , c h exp ( ( U RMA , i p , c h U Gauss p , c h ) 2 2 σ Gauss p , c h )
To achieve high sensitivity to residuals in the fit of P D F Gauss ( U RMA , i p , c h ) with respect to underlying measured P D F m ( U RMA , i p , c h ) , the following conditions are imposed on the numerical optimization:
  • The starting value of U Gauss p , c h is set to the mean U RMA p , c h of the measured sample U RMA , i p , c h .
  • The codomain of U Gauss p , c h is restricted by the lowest measurement to U Gauss p , c h min ( U RMA p , c h )
  • The starting value of σ Gauss p , c h is set to 0.02   mV corresponding to the typical standard deviation of a sample U RS , i c h measured on the resistive noise source.
  • The codomain of U ^ Gauss p , c h is restricted to U ^ Gauss p , c h 20
In the third step, the coefficient of determination R2 between P D F m ( U RMA , i p , c h ) and P D F Gauss ( U RMA , i p , c h ) is computed to quantify the Gaussianity of the statistics of a measured sample U RMA , i p , c h . Fourth, the threshold R2 < 0.95 is applied to flag the sample U RMA , i p , c h for exclusion from further analysis due to RFI. Upon successful Gaussian fitting ( R 2 0.95 ), the optimized mean value U Gauss p , c h of P D F Gauss ( U RMA , i p , c h ) is used to compute T B p , c h (using Equation (7)) instead of the mean U RMA p , c h of the measured sample U RMA , i p , c h .
The fifth step is the quantification of the distortion imposed by non-thermal RFI on the measurement of T B p , c h . The resulting measurement error Δ T B p , c h is estimated as:
Δ T B p , c h = | U Gauss p , c h U RMA p , c h | S RMA
Thereby, radiometer sensitivity S RMA p , c h 0.322   K · mV 1 (defined in Equation (7)) is used to “translate” the difference between the mean voltages U RMA p , c h and U Gauss p , c h of the measured sample U RMA ,   i p , c h and the fitted P D F Gauss to the difference Δ T B p , c h between brightness temperatures. The usefulness of the measurement uncertainty factor Δ T B p , c h lies in its application in retrievals based on L-band radiometry (such as in [14,15,16]) where T B p , c h can be weighted according to their non-thermal distortion levels. This essentially indicates that less influence on the retrieval results occurs from more disturbed measurements.
Figure 6 shows representative examples of PDFs associated with measured voltage samples U RMA , i p , c h (p = V, and ch = 1) including different levels of non-thermal RFI. Blue lines are P D F m ( U RMA , i V , 1 ) derived from measured samples U RMA , i V , 1 , red lines represent the Gaussian model P D F Gauss ( U RMA , i V , 1 ) fitted to the measurement based P D F m ( U RMA , i V , 1 ) . Figure 6a shows the case of an almost RFI-free sample, for which P D F m ( U RMA , i V , 1 ) and P D F Gauss ( U RMA , i V , 1 ) match well. Accordingly, the coefficient of determination R 2 = 0.99 is high, and the non-thermal contribution Δ T B V , 1 = 0.66   K is minor. Likewise, Kurtosis K = 2.77 and Skewness S = 0.35 indicate an almost RFI-free sample. Figure 6b shows an example of a distorted sample, which would not necessarily be recognized by simply looking at Kurtosis K = 3.06 and Skewness S = 0.46 values. However, the non-Gaussianity of P D F m ( U RMA , i V , 1 ) is obvious and becomes apparent in the reduced R 2 = 0.64 (RFI-flag is raised), and the significant RFI-distortion Δ T B V , 1 = 1.21   K . The PDFs shown in Figure 6c,d are examples of heavily RFI distorted measurements. As a consequence of the imposed conditions (indicated in the bullet points above) the Gaussian fits P D F Gauss ( U RMA , i V , 1 ) clearly fail, yielding R 2 = 0.0 in both cases. The RFI distortion shown in Figure 6c would also have been detected by means of the statistical moments K = 5.16 and S = 1.19. However, the example in Figure 6d would not necessarily be rejected as RFI-distorted based on K = 2.29 and S = 0.03 given their assumed thresholds of around 3.0 and 0.0, respectively.

4.3. Effective Transmission Losses and ACS Temperatures

As explained earlier, T B p corresponds with the noise T B , RMA p at the RMA-inputs p = H, V reduced by the transmission lines’ inherent noise (Equation (6)). Because transmission line (TL) noise depends on its physical temperature T TL and on its loss L TL p , the use of an “effective” value of L TL p is necessary to optimally minimize distortions of T B p imposed by temporal variations of T TL . The approach developed to achieve effective L TL p and optimal ACS noise-temperatures T ACS p , c h is outlined in the flowchart shown in Figure 7. This approach is based on the assumption that TL losses are optimally chosen if, for a range of TL-temperatures T TL = T air , deviations between theoretical sky-brightness and measurement-based T sky p are minimal.
Sky measurements were performed regularly throughout the campaign between 28 November and 3 March by pointing the radiometer skyward at nadir angle θ = 140 ° (= 40° relative to zenith). This resulted in N sky = 88 measurement cycles, each of which consists of 8 undisturbed voltage samples U ACS , i c h , U RS , i c h , and U RMA ,   i p , c h (i = 1,..., nps = 2400) and covers the temperature range of 20   ° C T air = T TL 14   ° C . To begin with, the first N 0 = 50 sample-means U ACS c h , U RS c h , and U RMA p , c h (28 November–11 January) with associated measurements T TL = T air were used together with initial guesses L TL p ini of transmission lines to compute the associated ACS temperatures T ACS p , c h (Section 4.1.1). From these N 0 foursomes of T ACS p , c h , the corresponding four temporal mean values 〈 T ACS p , c h 〉 were computed as illustrated in Boxs 1–3 in the flowchart (Figure 7). Next, the four mean values 〈 T ACS H , 1 〉, 〈 T ACS H , 2 〉, 〈 T ACS V , 1 〉, and 〈 T ACS V , 2 〉 were used to compute the four time-series T sky H , 1 , j , T sky H , 2 , j , T sky V , 1 , j , and T sky V , 2 , j (using Equation (6) with Equation (7)) based on the voltage sample-means recorded for the first j = 1 , , N 0 + k sky measurements. Furthermore, synchronous theoretical time series T sky j were simulated using [32] (Box 4 in the flow-chart). Naturally, the measurement-based time-series T sky p , c h , j do not perfectly agree with the simulated T sky j , and result in a Root-Mean-Square Error RMSE ( T sky p ) , defined in Box 5 of the flow chart.
The magnitude of RMSE ( T sky p ) depends on the assumption made regarding L TL p because it determines the degree to which fluctuations in T TL = T air translate into fluctuations of T sky p , c h , j . An excessively high value of L TL p causes T sky p , c h , j that are distorted by overcompensation of temperature-dependent TL-noise. An excessively low value of L TL p leads to an underestimation of this effect. Accordingly, the computation of effective TL-losses L TL p and ACS temperatures 〈 T ACS p , c h 〉 for a given time period (including j sky measurements) consists of iterating Steps 2–5 (Loop 6 in the flow chart) until RMSE ( T sky p ) becomes minimal (Box 7). As indicated above, the initial period, consisting of the first N 0 = 50 sky measurements, is used to compute L TL p and 〈 T ACS p , c h 〉 representative of this training period, and ultimately applied to derive T B p of footprint areas for this period. For later time periods, T B p are derived from effective values L TL p and 〈 T ACS p , c h 〉 estimated from the growing number j = 1 , , N 0 + k of sky measurements as indicted by the outermost Loop 8 in the flowchart shown in Figure 7.
Figure 8a shows the increasing trend of L TL p over time for both polarizations (p = H, V). This can be explained, for instance, by ageing and the accumulation of moisture in the TL-connectors as well as the λ / 4 dipoles within the Picket-horn antenna [17]. Throughout the entire campaign a total increase of about 0.08 dB and 0.21 dB is observed for L TL H and L TL V , respectively. After 23 January, the L TL H values almost saturate to 0.43 dB whilst L TL V continues increasing until the end of the campaign. The systematically larger L TL V , in comparison to L TL H , is consistent with the smaller bending radius of the V-polarization TL, compared to the H-polarization TL. This is because the λ / 4 dipole receiving vertical polarized radiance is mounted at the very back of the antenna-feed (see Figure 4 in [17]). Beyond this, the V-polarization λ / 4 dipole is more susceptible to contaminations because of its mounting position.
Figure 8b shows similar increasing trends for the four ACS noise temperatures 〈 T ACS p , c h 〉 (p = H, V; ch = 1,2) computed as the temporal mean of the first j = 1 , , N 0 + k individual ACS temperatures T ACS p , c h . This increase, however, is a direct consequence of increasing L TL p as can be seen when calculating Equations (1)–(3) for increasing L TL p . It was verified that the sample-means U ACS c h do not increase throughout the campaign. This confirms that increasing 〈 T ACS p , c h 〉 are not due to increasing noise power of the RMA internal ACS, but indeed results from ACS calibration based on sky measurements U RMA p , c h = U sky p , c h which increase with increasing TL losses.

4.4. Brightness Temperatures of “Reflector Areas”

Brightness temperatures T B , R p originating exclusively from the areas with the metal-mesh reflector beneath the snow (“reflector areas” shown in Figure 1) are useful to investigate snow emission segregated from ground emission. This is important, for instance, when investigating impacts of snow liquid water on the snowpack’s L-band emission. However, even when pointing the radiometer along the central line of the “reflector areas”, the antenna field of view still captures radiation emanating from areas surrounding the metal-mash reflector. Accordingly, it is necessary to separate T B , R p ( θ ) from T B , R , m p ( θ ) measured (m) along the azimuth direction of the reflector. The corresponding procedures, used similarly in [27,30,44], are outlined below.
Measurements (m) T B , R , m p ( θ ) taken along the azimuth direction of the “reflector area” are represented as the linear combination of radiances T B , R p ( θ ) and T B , N p ( θ ) originating from areas covered by the reflector (R) and natural (N) areas surrounding the reflector, respectively:
T B , R , m p ( θ ) = μ R p ( θ ) T B , R p ( θ ) + ( 1 μ R p ( θ ) ) T B , N p ( θ )
Here we further assume that emission from areas surrounding the reflector can be represented by measurements T B , N p ( θ ) performed on the “natural areas” shown in Figure 1. The factors μ R p ( θ ) and ( 1 μ R p ( θ ) ) in Equation (10) express, respectively, the weightings of T B , R p ( θ ) and T B , N p ( θ ) for nadir angles θ = 30°, 35°,…,65° and polarization p = H and V. Once μ R p ( θ ) is known, T B , R p ( θ ) are computed from measurements T B , R , m p ( θ ) and T B , N p ( θ ) using Equation (10) solved for T B , R p ( θ ) :
T B , R p ( θ ) = T B , R , m p ( θ ) ( 1 μ R p ( θ ) ) T B , N p ( θ ) μ R p ( θ )
The weighting factors μ R p ( θ ) in Equation (11) are derived on a daily basis using exclusively T B , R , m p ( θ ) and T B , N p ( θ ) measured quasi-simultaneously at night (00:00–07:00). It is essential to use exclusively night measurements for the computation of μ R p ( θ ) , which is representative of the corresponding day, because spatial heterogeneities across the field site are expected to be generally diminished during the night when temperature drops below freezing. This instance is in support of the assumption made that measurements T B , N p ( θ ) performed over “natural areas” (Figure 1) are also representative of areas surrounding the reflector. However, this assumption is still considered critical, especially during relatively warm and transient conditions where spatial variability of, for example, ground permittivity ε G is increased (as proven by corresponding in-situ measurements shown in Figure 4). However, throughout the main period of interest (i.e., cold winter and early spring periods), the ground is homogeneously frozen. Thus, using measurements T B , N p ( θ ) to represent emission of areas surrounding the “reflector areas” is reasonable in most cases.
Another important reason to estimate μ R p ( θ ) from night measurements of T B , R , m p ( θ ) and T B , N p ( θ ) is that a surface-thawed snowpack refreezes overnight, and thus snow emission during cold nights is almost always zero [13,22,23]. As a consequence, during the cold nights throughout our main periods of interest (cold winter and early spring), T B , R p ( θ ) correspond with downwelling sky brightness temperatures T sky ( θ ) reflected by the metal-mash reflector towards the antenna aperture. Accordingly, T B , R p ( θ ) = T sky   ( θ ) is used in Equation (10), yielding the following expression to compute μ R p ( θ ) with T sky   ( θ ) simulated by the model [32]:
μ R p ( θ ) = T B , N p ( θ ) T B , R , m p ( θ ) T B , N p ( θ ) T sky ( θ )
In summary, the weighting factors μ R p ( θ ) are computed with Equation (12) for each day using the corresponding night measurements of T B , R , m p ( θ ) ,   T B , N p ( θ ) and with T sky   ( θ ) simulated using [32]. μ R p ( θ ) is then kept constant throughout that day and used in Equation (11) to derive T B , R p ( θ ) for that day from the corresponding quasi-synchronous measurements T B , R , m p ( θ ) and T B , N p ( θ ) .

5. Analysis of Brightness Temperatures

5.1. Methodology for the Analysis of L-Band Brightness Temperatures’ Response to Snow Liquid Water

The model used to simulate L-band brightness temperatures T B p of a rough ground surface beneath a snowpack is a composite of components adopted from the L-band Microwave Emission of the Biosphere “L—MEB” model [45] and the Microwave Emission Model of Layered Snowpacks “MEMLS” [10,46] while considering simplifications applicable to L-band. The resultant L-band Specific emission model, henceforth denoted as “LS—MEMLS”, is now outlined following the model structure illustrated in Figure 9.
Block 1 of the diagram (Figure 9a) shows the “LS—MEMLS” inputs, which are the radiometer nadir angle θ, polarization p = H (for horizontal) and p = V (for vertical), and the L-band vacuum wavelength λ = 21 cm. Air temperature T air   2 m above ground, specific air humidity q air , rain rate r air , and the altitude h site of the test site above sea-level are additional inputs necessary to simulate downwelling L-band sky radiance T sky using the model described in both [32] and [33].
The “LS—MEMLS” inputs given to each of the layers j of the snowpack (S) are the layer-thickness d S , j , the physical temperature T S , j , the mass density ρ S , j , and the volumetric liquid water content W S , j . Inputs used to characterize the ground (G) beneath the snowpack are effective ground temperature T G , ground permittivity ε G , and the HQN ground roughness model [45,47] parameters h G , q G and n G p (p = H, V) used and described in [48]. Intermediate quantities computed along with the evaluation of “LS—MEMLS” are listed in Block 2 of Figure 9a, and outlined here.
As mentioned above, L-band sky radiance T sky is simulated from the “LS—MEMLS” inputs θ, T air , q air , r air , h site based on [32]. The same approach is also used in our earlier work [14] to represent T N + 1 = T sky incident on the uppermost interface j = N + 1 of the snowpack. For h site = 1450   m (specific to the Davos Laret site), T air = 273.15   K , q air = 10   g   kg 1 , r air = 0   mm   h 1 , and the considered nadir angles 30 ° θ 65 ° , sky radiance is in the range 4.2   K T sky 5.8   K , revealing that the atmosphere is largely transparent at L-band.
Complex snow permittivities ε S , j = ε S , j + i · ε S , j are computed with the Maxwell–Garnett dielectric model considering a two-phase mixture of liquid water with permittivity ε W and dry snow with permittivity ε S , dry , j :
ε S , j = ( 1 W S ) ε S , dry , j + W S ε W K 1 W S ( 1 K )   with   K = K a + K b + K c 3 and K a = ε S , dry , j ε S , dry , j + A a ( ε W ε S , dry , j ) K b = K c = ε S , dry , j ε S , dry , j + A b ( ε W ε S , dry , j )
The same approach is implemented in the latest version of “MEMLS” (pages 15–16 in the “MEMLS-III” manual [49]). Liquid snow water with volumetric content W S is taken as prolate ellipsoidal inclusions in dry snow considered as the dielectric background, and the depolarization factors are assumed as A a = 0.005 , A b = A c = 0.4975 [23,50].
Permittivity ε W = 85.82 + i · 12.64 at the frequency 1.4 GHz and at freezing temperature (273.15 K) is computed with the model [51]. Permittivity ε S , dry , j of the dry snow phase is computed from snow mass-density ρ S , j using the same model as in “MEMLS” [10,46] ( ρ S = 0.001 · ρ S , j in units of [ g   cm 3 ] ):
ε S , dry , j ( ρ S ) = { 1 + 1.5995 ρ S + 1.861 ρ S 3   for   0   g cm 3 ρ S 0 . 4   g cm 3 ( ( 1 ν ) ε h + ν ε b ) 3   for   ρ S 0 . 4   g cm 3
with υ = ρ S / 0.917 , ε h = 0.99913 and ε b = 1.4759 . For dry snow ( W S = 0.0   m 3 m 3 ) with mass density 0   kg   m 3 ρ S 400   kg   m 3 , imaginary parts 0 ε S 1.2 · 10 4 are very small and significantly smaller than real parts 1 ε S 1.759 . This explains the very large L-band penetration depth (>300 m) in dry snow [22,23], which allows a dry snow layer to be treated as transparent (transmissivity t S , j = 1 ). However, for a moist snow layer, t S , j rapidly decreases with increasing W S , j . Accordingly, for W S , j > 0.0   m 3 m 3 , transmissivity t S , j is computed from the layers nadir optical-depth d S , j · α S , j and the propagation angle θ S , j across the layer j with thickness d j and the power absorption coefficient α S , j . Following Beer’s and Snell’s law, this yields:
t S , j = exp ( d S , j α S , j / cos θ S , j )
and
θ S , j = arcsin ( sin ( θ ) / | ε S , j | )
The power absorption coefficient α S , j within the snow-layer j is computed from the respective snow permittivity ε S , j = ε S , j + i · ε S , j as:
α S , j = 4 π λ Im ( ε S , j )
In accordance with earlier work [13,14], the volume reflectivities of snow layers are assumed to be r S , j = 0 , indicating the fact that volume scattering caused by snow microstructures is negligible at L-band. Clearly, neglecting snow volume scattering is applicable only in absence of prominent macrostructures (such as ice lenses) with dimensions of the order of λ.
Specular reflectivities s S , j p (p = H, V) of the upper interface of each layer j result from the dielectric contrast between the respective layer permittivities ε S , j and ε S , j + 1 . Accordingly, s S , j p are computed with Fresnel’s equations considering ε out = ε S , j and ε in = ε S , j + 1 as the permittivities of the layers containing the propagated (out) and the incident (in) waves, respectively. Furthermore, the incident-angle θ in = θ S , j + 1 at the interface j = 1,…,N−1 is computed using Equation (16) to correspond with the propagation angle within the above layer j+1:
s S , j H = | ε in A ε out B ε in A + ε out B | s S , j V = | ε out A ε in B ε out A + ε in B | with A = cos θ in B = 1 ( 1 A 2 ) ε in / ε out ε out = ε S , j ε in = ε S , j + 1
Likewise, the reflectivities s S , N H , s S , N V of the snow-air interface (j = N) are computed via Equation (18) considering ε in = ε air = 1 , ε out = ε S , j and the radiometer nadir angle θ in = θ .
In contrast to the j = 1,…,N specular snow-layer interface reflectivities s S , j p , the effects of ground roughness are indeed considered in the computation of the reflectivities s G p = s S , 0 p of the ground-snow interface (j = 0). Thereto, first the specular (*) ground-snow interface reflectivities s G H and s G V are computed using Equation (18) with ε in = ε S , 1 , ε out = ε G and θ in = θ S , 1 . Accordingly, the dominant effect of snow refraction is incorporated into “LS—MEMELS” at this point. Next, the effect of ground surface roughness is taken into account by correcting the specular reflectivities s G H and s G V using the semi-empirical HQN roughness model [45,46,47]
s G H = exp ( h G ( cos θ S , 1 ) n G H ) ( s G H ( 1 q G ) + q G s G V ) s G V = exp ( h G ( cos θ S , 1 ) n G V ) ( s G V ( 1 q G ) + q G s G H )
Detailed information on the HQN roughness model and its typical parameter values for different types of landscapes can be found in Section 3.1.2 and Table II in [52]. However, in all of our investigations, roughness effects are assumed to be independent of the incidence angle θ S , 1 on the ground. For this reason, we set n G V = n G H = 0 . The polarization mixing parameter is assumed to be q G = 0.05 . The most influential HQN roughness model parameter is the dimensionless effective surface roughness parameter h G , which is mostly correlated with “geometric roughness”. The values ultimately selected for the four HQN model ground surface roughness parameters ( h G ,   q G ,   n G V ,   n G H ) = ( 0.1 ,   0.05 ,   0.0 ,   0.0 ) are in agreement with the values used in the current SMOS soil moisture retrieval algorithm [53].
As indicated in Block 3 in Figure 9a, final outputs T B p of “LS—MEMLS” are computed with the two-stream emission model (2S-EM) employed in “MEMLS” [10,46] using the intermediate quantities (for snow-layers j :   t S , j ,   r S , j ,   T S , j ; for the ground: s G p = s S , 0 p ; for the sky: T sky ) made available as outlined above. The 2S-EM “balances” in each layer the up- and down-welling energy streams (fluxes) propagated by the electromagnetic waves by taking into account the boundary conditions at layer interfaces and considering Kirchhoff’s law. These (incoherent) linkages between up- and down-welling layer-fluxes take into consideration multiple reflections between the layers’ interfaces. Formulating the resulting four flux-equations for each layer (up- and down-welling at the respective two layer interfaces) yields a matrix equation for the vector T B p = ( T B , 1 p , T B , 2 p , , T B , j p , , T B , N 1 p ,   T B , N p , T B p ) containing the upwelling radiances T B , j p across the j = 1,…,N snow-layers, and, as the last element N + 1, the upwelling radiance above the uppermost snow layer corresponding to the sought T B p = T B , N + 1 p above the snow covered ground. The detailed derivation of the matrix formulation of the 2S-EM can be found in Section 3 of [49].
In the alternative formulation sketched in Figure 9a, T B p is written as T B p = Σ j = 0 N + 1 a j p · T j whereas the Kirchhoff coefficients a 0 p = a G p , a j p = a S , j p (j = 1,…,N), and a j + 1 p = a sky p fulfill Σ j = 0 N + 1 a j p = 1 and weight the respective temperatures T 0 = T G , T j = T S , j (j = 1,…,N), and T j + 1 = T sky . This so-called Kirchhoff formulation is appropriate for small numbers N of snow layers. In the case of N = 1 snow layer ( s S , 1 p = s S p , t S , 1 = t S ) and assuming r S , 1 = r S = 0 , the associated 2S-EM Kirchhoff coefficients read (Appendix of [14]):
a G p = ( 1 s G p ) ( 1 s S p ) t S 1 + r S 2 s G p s S p r S ( s G p + s S p ) s G p s S p t S 2 a S p = ( 1 s G p ) ( 1 r S t S ) ( 1 r S s G p + s G p t S ) ( 1 r S s G p ) ( 1 r S s S p ) s G p s S p t S 2 a sky p = 1 a G p a S p
Brightness temperatures T B p above a homogeneous moist ( t S < 1 ) snowpack above a given ground is therefore simulated as:
T B p = a G p T G + a S p T S + a sky p T sky
This one-layer version of “LS—MEMLS” is employed in Figure 10a to simulate impacts of uniformly distributed snow liquid water on brightness temperatures. Furthermore, the simplicity of the one-layer version is the method of choice when “LS—MEMLS” is used in an iterative retrieval scheme. For more complex assumptions made, for instance, on the vertical distribution of snow liquid water (Figure 10b–d), the matrix formation of the 2S-EM is used.

5.2. Sensitivities of Brightness Temperatures to Liquid Snow Water

Following the model Equations (13) and (14), effective snow permittivities ε S = ε S + i · ε S depend on both the snow mass-density ρ S and its volumetric liquid water content W S . For W S = 0   m 3 m 3 , power absorption coefficients (Equation (17)) of dry snow ( W S = 0   m 3 m 3 ) with density ρ S 500   kg   m 3 are α S < 10 6   m 1 , implying that any seasonal snowpack can be considered transparent ( t S 1 ) as long as it is dry. However, this changes drastically when W S > 0   m 3 m 3 , mainly as a consequence of the strong increase in the imaginary part ε S of snow permittivity with increasing W S . For example, for wet snow with W S = 0.05   m 3 m 3 , or even moist snow with W S = 0.01   m 3 m 3 the respective absorption coefficients α S 7.4545   m 1 and α S 1.44034   m 1 are more than six orders of magnitude higher than for dry snow. From a practical point of view, this rapidly increasing absorption of snow with increasing W S led, in the first attempt, to the omission of L-band-based retrievals ( ρ S ,   ε G ) of snow density and ground permittivity for not entirely dry snow conditions [15]. This, however, makes it even more important to explore the sensitivities of L-band T B p ’s with respect to snow liquid water in support of (i) an in-depth understanding of the measured brightness temperatures discussed in Section 5.3; and (ii) the potential to estimate snow wetness based on L-band brightness temperatures.
Accordingly, “LS—MEMLS” is used to explore sensitivities of L-band T B , R p expected for moist snow atop a perfectly reflecting (R) ground and T B , N p expected for moist snow atop frozen natural (N) ground with respect to snow liquid water. Figure 10 shows simulated T B , R p ( W C S ) (thin lines) and T B , N p ( W C S ) (bold lines) for the nadir angles θ = 0 ° (black), θ = 30 ° (red), and θ = 60 ° (green). Snow mass-density is considered as ρ S = 300   kg   m 3 ; ground- and snow temperatures are assumed as T G = T S = 273.15   K . The parameter varied is the snow liquid water column W C S = 0 h s W S ( z ) d z , defined as the volumetric liquid water content W S , integrated over the entire snow depth h S . Figure 10a is for “uniform” W S = W C S / h S across snowpacks with 0   mm W C S 70   mm and snow-depth h S = 0.5   m . Further simulations are shown for the same range of W C S contained in snowpacks with a moist snow-layer of thickness d WS = 0.1   m either at the “top” (Figure 10b), “sandwiched” in between the dry snow (Figure 10c), or at the “bottom” of the dry snowpack (Figure 10d). These types of W S -profiles are hereafter referred to as “uniform”, “top”, “sandwiched”, and “bottom”. Corresponding simulations are also conducted for ρ S = 500   kg   m 3 yielding results fully consistent with those presented below.
For totally dry snowpacks, brightness temperatures T B , R p ( W C S = 0   mm ) necessarily agree with the downwelling sky brightness temperature assumed as T sky = 5   K and specularly reflected at the ground (black open squares in Figure 10). A further characteristic common to all the T B , R p ( W C S ) is the sharp increase with increasing W C S for small values of W C S , followed by gradually decreasing T B , R p ( W C S ) for larger W C S . The steeply increasing tail of T B , R p ( W C S ) is explained by the strong increase in snow emission with increasing snow liquid water (see Equations (13)–(17)). The gradients d T B , R p / d W C S become even larger for increasing nadir angles θ as a result of longer propagation path-lengths across the wet snow. Furthermore, the steeply increasing tail of T B , R p ( W C S ) at small W C S is largely independent of polarization because volume emission of wet snow is assumed to be isotropic. The cause of the gradually decreasing tail of T B , R p ( W C S ) seen with larger W C S , is caused by increased interface reflectivity in the upper boundary of the wet snow-layer. The latter increases with increasing W S as a result of likewise increasing permittivity (Equation (13)). Moreover, the gradually decreasing part of T B , R p ( W C S ) depends significantly on θ and polarization p, which is a consequence of the respective sensitivities of the Fresnel-like interface reflectivities (Equation (18)).
The plateau of T B , R p ( W C S ) occurs at intermediate values of W C S when increasing snow volume emission compensates for decreasing interface emission. The impact of different assumptions made on the depth-profile of snow liquid water content W S becomes apparent when comparing Figure 10a–d. For the “uniform” scenario (Figure 10a), the steeply ascending part of T B , R p ( W C S ) shows larger gradients d T B , R p / d W C S 0   K   mm 1 , whereas the respective descending part reveals smaller sensitivities d T B , R p / d W C S < 0   K   mm 1 compared to the scenarios “top” (Figure 10b), “sandwiched” (Figure 10c), and “bottom” (Figure 10d). This finding suggests that the assumption of uniformly distributed snow liquid water may lead to biased W C S -retrievals estimated from T B , R p when liquid snow water is actually, for example, mostly present below the snow surface. Regarding the relevance of each of the considered W S -scenarios, it should be mentioned that the corresponding scenario “sandwiched” is seen as the most realistic during the mid-winter period with pronounced diurnal temperature variations around 0 °C. This is because strong solar short-wave radiation penetrates the first few centimeters of dry snow and can cause partial subsurface melting—a phenomenon that can be best called the “snow greenhouse effect”.
Different snow mass-densities ρ S = 300   kg   m 3 (Figure 10) and ρ S = 500   kg   m 3 manifest primarily in moderate offsets of T B , R p ( W C S ) , while gradients d T B , R p / d W C S are not affected greatly. The predominant and very high sensitivities of T B , R p with respect to W C S imply the potential to retrieve snow liquid water column W C S from measurements T B , R p performed over the “reflector areas” with high sensitivity, especially for the lower range 0   mm W C S 20   mm .
Regarding the simulations T B , N p ( W C S ) (bold lines in Figure 10), representing the emission of snow atop a ground with ε G = 5 , the same four W S -scenarios “uniform” (Figure 10a), “top” (Figure 10b), “sandwiched” (Figure 10c), and “bottom” (Figure 10d) are considered. Polarization is indicated as p = H (solid), V (dashed); nadir angles θ = 0 ° (black), θ = 30 ° (red), and θ = 60 ° (green), all with the considered snow density ρ S = 300   kg   m 3 (Figure 10).
For the completely dry snowpack ( W C S = 0   mm ) and θ = 0 ° (black lines), the respective T B , N p ( W C S = 0   mm ) are slightly larger than the brightness temperatures above the snow-free ground (black solid circles) because of the impedance matching of the snow as explained in [13]. For θ > 0 ° and W C S = 0   mm , the T B , N p ( W C S = 0   mm ) deviate from corresponding bare-ground brightness temperatures (red circles for θ = 30 ° , green circles for θ = 60 ° and solid for p = H, open for p = V). This deviation results from not only impedance matching, but also refraction effects caused by snow interfaces [13,14]. Most of the T B , N p ( W C S ) show an increasing and a decreasing tail separated by a plateau. Again, the increasing tails of T B , N p ( W C S ) result from the increasing emission of the snow-volume with increasing snow liquid water. However, they are distinctly less pronounced than the corresponding increasing tails simulated for the T B , R p ( W C S ) .
This is the compelling consequence of the much larger emission of the frozen natural ground (with ε G = 5 ) compared to the zero-emission of the ground covered with the reflector. The decreasing tails of T B , N p ( W C S ) , caused by increased interface reflectivities, are more pronounced than the decreasing tails of T B , R p ( W C S ) . This is because snow-ground interface reflectivities are much more comparable with snow-sir interface reflectivities in the case of “natural” ( ε G = 5 ) than in the case of the ground covered with the reflector. Accordingly, T B , N p is diminished by both increased reflectivities at the snow-sir and ground-snow interfaces, while T B , R p is exclusively diminished by increased reflectivities at the snow-sir interface. Of course, snow-ground interface reflectivities are only affecting brightness temperatures as long as the snowpack is not totally opaque. Accordingly, when brightness temperatures become dominated by snow volume emission for large W C S , the T B , N p ( W C S ) and T B , R p ( W C S ) necessarily approach each other. In summary, the discussed sensitivities of T B , N p with respect to W C S seem to be sufficiently high to achieve at least reliable qualitative information on snow wetness states based on L-band brightness temperatures measured over natural frozen grounds.

5.3. Analysis of Measured Brightness Temperatures

We now present the time series of measured calibrated brightness temperatures T B , R p ( θ ) and T B , N p ( θ ) , and discuss examples of their prominent signatures in response to selected snow and ground conditions. Figure 11a shows the time series of T B , N H and T B , R H averaged over nadir angles 40 ° θ 50 ° with blue and magenta lines, respectively. Averaging over θ is performed to achieve more measured brightness temperatures after RFI filtering. The onset of snow cover on 3 January and the beginning of “early spring” on 31 January are indicated by vertical dashed lines. Thus, as labeled above Figure 11, the entire time series is broken into three periods: the “snow-free period”, the “cold winter period” with frozen ground and mostly T air < 0 °   C , and the “early spring period” with air temperatures rising above the freezing point. Figure 11b shows T B , R V ( θ ) and T B , N V ( θ ) . Figure 11c shows air temperature T air and precipitation rate, respectively.
As explained in Section 5.2 and shown with black open squares in Figure 10, T B , R p ( θ ) (p = H, V) simulated for the “reflector areas” equals downwelling sky radiance T sky = T B , R p if the reflector is snow-free and also if the reflector is covered with entirely dry snow. These model findings are corroborated by the measurements shown in Figure 11. Throughout the snowless period, T B , R H (Figure 11a) is close to 4.5 K, which is virtually equal to T sky 4.7   K for 40 ° θ 50 ° [32]. This holds true also after the onset of the snow cover throughout the “cold winter period” in January when the snowpack is dry. However, with the beginning of the “early spring period”, T B , R H occasionally rises distinctly above T sky , which is interpreted as intermittent occurrences of snow liquid water. This is, again, consistent with the corresponding response of simulated T B , R H shown in Figure 10. On the other hand, the perception of intermittent increases in measured T B , R H caused by snow liquid water is also meaningful from the environmental conditions viewpoint because time-integrated heat-input to the snowpack starts temporarily rising above the ice latent-heat during the “early spring period”. The distinct sensitivity of T B , R H with respect to snow liquid water becomes even more evident during the period 9–23 February, when T B , R H reaches its maximum every day in the afternoon and falls to T sky over night when the snowpack refreezes. Accordingly, the measured temporal characteristics of T B , R H clearly indicate the onset of “early-spring snow”.
Additionally, during the aforementioned two weeks in February, where the T B , R H shows exclusively sharp increases during the afternoon hours, the concurrent T B , N H measured over the “natural areas” show both sharply increasing (11–13 February) and decreasing (16–19 February) deflections. Although at first glance these opposite measurement responses of T B , R H and T B , N H to snow wetness may seem counter-intuitive, it is in qualitative agreement with the corresponding simulations shown in Figure 10. On the one hand, the simulated T B , R H ( W C S ) and T B , N H ( W C S ) in Figure 10, reveal sharply increasing T B , R H ( W C S ) for 0   mm W C S 20   mm explained by increased snow emission. On the other hand, T B , R H ( W C S ) is expected to decrease for W C S 20   mm , while simulated T B , N H ( W C S ) already decrease for lower W C S 5   mm . In other words, for 0   mm W C S 20   mm , the T B , R H ( W C S ) are expected to increase steadily, while T B , N H ( W C S ) increase with very little snow liquid-water ( 0   mm W C S 5   mm ) and decrease with larger values W C S 5   mm . These simulation-based findings explain why measurements T B , R H and T B , N H can indeed respond in contrast to each other with respect to low W C S values, as is expected from 9 to 23 February. The fact that “LS—MEMLS” is an incoherent emission model also implies that it is not necessary to adduct coherent effects to explain the observed oppositional responses of T B , R H and T B , N H to snow liquid water.
As mentioned earlier, Figure 11b shows similar information to Figure 11a only with “reflector areas” and “natural areas” brightness temperatures T B , R V (magenta) and T B , N V (blue) at vertical polarization. The signature in T B , R V (Figure 11b) associated with the first appearance of snow liquid water (“early-spring snow”) coincides with the corresponding response of T B , R H (Figure 11a). However, three additional distinct peaks in T B , R V are observed during measurements on 3, 13, and 18 of January. For certain, these peaks are not caused by snow liquid water because T air (Figure 11c) was well below the freezing point during the “cold winter period”. The still very shallow and homogeneous snow-depth present during this period suggests that these peaks are caused by coherent effects. This supposition is further borne out by the fact that the amplitudes of the three peaks decrease with time and, consequently, with growing snow depth (see Figure 2a). Likewise, the excessively high T B , R V measured for 9 h during 7 and 8 of February are likely the result of coherent effects. In the case of the respective measurements during the beginning of the “early spring period” the suspected coherent effect is most likely induced by coherent layers of moist snow rather than by shallow snow height h S . This is because h S > 0.5   m (see Figure 2a) during the “early spring period” is significantly larger than the maximum thickness 3 λ / 8 ε S 0.07   m of a snow layer expected to induce coherent effects (see [10,49]) at L-band ( λ 0.21   m ), while a coherent moist snow layer can develop and vanish as a result of refreezing during a day of this period.
Figure 11d shows frequent precipitation at the end of February and early March. Most of these precipitation events are rain or wet snowfall, which naturally increase snow liquid-water distinctly. Thus, the higher T B , R H and T B , R V observed over the “reflector areas” for this period are consistent with simulated responses shown in Figure 10. One should note that increased T B , R p (p = H, V) is almost exclusively a result of increased snow wetness. Thus, neither precipitation events alone nor T air 0   ° C imply increased T B , R p because dry snowfall does not cause more liquid snow water and snow melt-down does not occur for T air 0   ° C unless the snowpack has reached ice latent heat.
As outlined in Section 5.1, snow liquid-water distinctly increases snow absorption. This impact is so large that the L-band penetration depth of >300 m in dry snow drops to less than 3 m for moist snow with liquid water content of ~1% and to even less than 0.3 m for wet snow with liquid water content of ~3% [22,23]. Accordingly, a natural snowpack can eventually become largely opaque especially during late winter periods, implying that T B , R p are no longer dominated by the very low emissivity of the metal reflector. Accordingly, we define the snowpack as “opaque” if more than 63% of the total emission originates from the snowpack, whereas up to this limit the snowpack is considered “semi-transparent”. Considering that the effective snow temperature for moist snow is necessarily close to the freezing-point of water, 273.15 K, the upper limit of T B , R p indicating an “opaque” snowpack is 273.15 × e 1 172   K . The horizontal dashed lines in Figure 11a,b indicate this threshold. It can be seen that the great majority of T B , R p (p = H, V) measurements are below this threshold, and as a result snow is still classified as “semi-transparent” rather than “opaque”. This experimental result validates the theoretically predicted advantage of L-band radiometry over higher frequency bands because of moist snow semi-transparency which is an absolute precondition to estimating snowpack column properties (such as snow liquid water-column) from remote sensing data.
Finally, we briefly interpret the T B , N p (p = H, V) measurements during the “snow-free period” before 3 January. It can be seen in Figure 11a,b that T B , N p in both polarizations p = H and V demonstrate an evanescing fluctuation pattern synchronous and anti-correlated with the diurnal variations of T air . This anti-correlation implies that the T B , N p fluctuations are not a direct result of T air diurnal variations but are rather a dielectric effect. To understand these fluctuations, one should note that in the absence of a snow cover acting as a heat-insulator, daily partial thaw-up takes place at a thin surface layer of bare ground which results in slightly higher ε G and lower emissivity. This change in emissivity results in lower and higher T B , N p during days and nights, respectively. The amplitude of these fluctuations decreases with deeper freezing of soil resulting from decreasing average air temperature, shorter days, and so on. The exception to this trend takes place on 26 and 27 December when a rain precipitation event moisturizes the soil surface resulting in significantly lower T B , N p for both p = H and V, but after 27 December, the soil refreezes and T B , N p returns to the approximately same values as before the rain event.

6. Summary and Conclusions

The present work is to be seen as the scientific basis for upcoming research on L-band remote sensing applied to the retrieval of ground freeze/thaw state and snow volume properties, such as mass density and liquid water content. Accordingly, details regarding the Davos Laret field laboratory and its instrumentation during its first winter of operation (2016/2017) are here explained. The in-situ measurements of ground and snow state parameters, climatological data, as well as L-band brightness temperatures recorded during the first winter 2016/2017 campaign, were presented and discussed. These data provide insight into how ground-snow system evolution is influenced by the climatological drivers. The temporal evolution of the measured data suggests a three-phase partitioning of the entire campaign into “snow-free period”, “cold winter period”, and “early spring period”.
As part of ELBARA-II’s raw data processing to achieve calibrated brightness temperatures, a refined RFI mitigation approach is suggested based on fitting a Gaussian model to the distribution of measured raw-data voltage samples. Its better performance compared to conventional “normality” tests (kurtosis, skewness) is demonstrated, and its ability to quantify non-thermal disturbances T B p , c h in individual measurements is highlighted. It is considered prudent to employ T B p , c h in future retrieval schemes using data measured with radiometers that output raw data in time-domain, such as ELBARA-II and the SMAP radiometer [54]. A further improvement in the calibration of ELBARA-II brightness temperatures concerns the use of values L TL p of “effective” transmission losses, which optimally represent the extent of associated thermal noise. Accordingly, the introduced refinement minimizes under- and over-estimation of thermal noise caused by the lossy transmission line at different physical temperatures and ultimately imposed on the calibrated brightness temperatures.
The response of the brightness temperatures T B , R p and T B , N p , over “reflector areas” and “natural areas”, with respect to snow liquid water column W C S is analyzed and presented. These analysis results reveal the high sensitivity of T B , R p and moderate, but measureable, sensitivity of T B , N p with respect to the liquid water column W C S . These model findings were corroborated by corresponding measured time series of T B , R p and T B , N p . For example, the measured T B , R p equals the sky brightness temperature T sky 4.7   K for the “snow-free period” and the “cold winter cold” when the snowpack is dry. This is predicted by the corresponding simulations of T B , R p for dry snow. It was also shown that distinct daily increases in measured T B , R p indicate the onset of the snow melting season, otherwise called “early-spring snow”. The consistency between modeled and measured responses of T B , R p and T B , N p suggests the need to explore interfering effects of snow liquid water on estimates ( ρ S , ε G ) of snow density and ground permittivity retrieved with the approach [14], which inherently assumes snow to be dry. Even more noticeably, our analysis hints at the possibility of estimating snow liquid water column W C S using L-band radiometry.

Acknowledgments

This paper was supported in part by the Swiss National Science Foundation under Grant 200021L_156111 /1, and in part by ESA within the framework of ESA’s SMOS “Expert Support Laboratory” (ESL) for Level 2. Furthermore, we acknowledge the Helmholtz Centre Potsdam (GFZ) for providing the ELBARA-II radiometer funded via the German Helmholtz Associations TERrestrial ENvironmental Observatories (TERENO) infrastructure project, and Matthias Jaggi from the WSL Institute for Snow and Avalanche Research SLF for conducting in-situ snow characterization measurements. We thank Curtis Gautschi for editing the manuscript as well as Manfred Stähli for his helpful comments. In addition, the authors would like to hereby express their appreciation to the anonymous reviewers for their valuable comments and inputs which helped with further improvement of this work.

Author Contributions

Reza Naderpour and Mike Schwank conceived, designed, and performed the experiments in equal shares. Likewise, they contributed equally to the data analyses and the writing of the manuscript. Christian Mätzler provided valuable scientific inputs to the research presented.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mätzler, C. Cost action 712: Microwave radiometry. In Remote Sensing of Atmosphere and Ocean from Space: Models, Instruments and Techniques; Springer: Berlin, Germany, 2002; pp. 231–246. [Google Scholar]
  2. Jackson, T.J.; Schmugge, T.J. Passive microwave remote sensing system for soil moisture: Some supporting research. IEEE Trans. Geosci. Remote Sens. 1989, 27, 225–235. [Google Scholar] [CrossRef]
  3. Rautiainen, K.; Lemmetyinen, J.; Pulliainen, J.; Vehviläinen, J.; Drusch, M.; Kontu, A.; Kainulainen, J.; Seppänen, J. L-Band Radiometer Observations of Soil Processes in Boreal and Subarctic Environments. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1483–1497. [Google Scholar] [CrossRef]
  4. Rautiainen, K.; Lemmetyinen, J.; Schwank, M.; Kontu, A.; Ménard, C.B.; Mätzler, C.; Drusch, M.; Wiesmann, A.; Ikonen, J.; Pulliainen, J. Detection of soil freezing from L-band passive microwave observations. Remote Sens. Environ. 2014, 147, 206–218. [Google Scholar] [CrossRef]
  5. Macelloni, G.; Brogioni, M.; Pampaloni, P.; Cagnati, A.; Drinkwater, M.R. DOMEX 2004: An experimental campaign at Dome-C Antarctica for the calibration of spaceborne low-frequency microwave radiometers. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2642–2653. [Google Scholar] [CrossRef]
  6. Roy, A.; Toose, P.; Williamson, M.; Rowlandson, T.; Derksen, C.; Royer, A.; Berg, A.A.; Lemmetyinen, J.; Arnold, L. Response of L-Band brightness temperatures to freeze/thaw and snow dynamics in a prairie environment from ground-based radiometer measurements. Remote Sens. Environ. 2017, 191, 67–80. [Google Scholar] [CrossRef]
  7. Zheng, D.; Wang, X.; van der Velde, R.; Zeng, Y.; Wen, J.; Wang, Z.; Schwank, M.; Ferrazzoli, P.; Su, Z. L-Band Microwave Emission of Soil Freeze-Thaw Process in the Third Pole Environment. IEEE Trans. Geosci. Remote Sens. 2017, 55, 5324–5338. [Google Scholar] [CrossRef]
  8. Mätzler, C.; Hüppi, R. Review of Signature Studies for Microwave Remote Sensing of Snowpacks. Adv. Space Res. 1989, 9, 253–265. [Google Scholar] [CrossRef]
  9. Mätzler, C. Improved Born Approximation for scattering of radiation in a granular medium. J. Appl. Phys. 1998, 83, 6111–6117. [Google Scholar] [CrossRef]
  10. Wiesmann, A.; Mätzler, C. Microwave emission model of layered snowpacks. Remote Sens. Environ. 1999, 70, 307–316. [Google Scholar] [CrossRef]
  11. Liang, D.; Xu, X.; Tsang, L.; Andreadis, K.M.; Josberger, E.G. The Effects of Layers in Dry Snow on Its Passive Microwave Emissions Using Dense Media Radiative Transfer Theory Based on the Quasicrystalline Approximation (QCA/DMRT). IEEE Trans. Geosci. Remote Sens. 2008, 46, 3663–3671. [Google Scholar] [CrossRef]
  12. Harlow, R.C.; Essery, R. Tundra Snow Emissivities at MHS Frequencies: MEMLS Validation Using Airborne Microwave Data Measured During CLPX-II. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4262–4278. [Google Scholar] [CrossRef]
  13. Schwank, M.; Rautiainen, K.; Mätzler, C.; Stähli, M.; Lemmetyinen, J.; Pulliainen, J.; Vehviläinen, J.; Kontu, A.; Ikonen, J.; Ménard, C.B.; et al. Model for microwave emission of a snow-covered ground with focus on L band. Remote Sens. Environ. 2014, 154, 180–191. [Google Scholar] [CrossRef]
  14. Schwank, M.; Mätzler, C.; Wiesmann, A.; Wegmüller, U.; Pulliainen, J.; Lemmetyinen, J.; Rautiainen, K.; Derksen, C.; Toose, P.; Drusch, M. Snow Density and Ground Permittivity Retrieved from L-Band Radiometry: A Synthetic Analysis. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 3833–3845. [Google Scholar] [CrossRef]
  15. Lemmetyinen, J.; Schwank, M.; Rautiainen, K.; Kontu, A.; Parkkinen, T.; Mätzler, C.; Wiesmann, A.; Wegmüller, U.; Derksen, C.; Toose, P.; et al. Snow density and ground permittivity retrieved from L-band radiometry: Application to experimental data. Remote Sens. Environ. 2016, 180, 377–391. [Google Scholar] [CrossRef]
  16. Naderpour, R.; Schwank, M.; Mätzler, C.; Lemmetyinen, J.; Steffen, K. Snow Density and Ground Permittivity Retrieved From L-Band Radiometry: A Retrieval Sensitivity Analysis. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 3148–3161. [Google Scholar] [CrossRef]
  17. Schwank, M.; Wiesmann, A.; Werner, C.; Mätzler, C.; Weber, D.; Murk, A.; Völksch, I.; Wegmüller, U. ELBARA II, An L-Band Radiometer System for Soil Moisture Research. Sensors 2010, 10, 584–612. [Google Scholar] [CrossRef] [PubMed]
  18. Kerr, Y.H.; Waldteufel, P.; Wigneron, J.P.; Delwart, S.; Cabot, F.; Boutin, J.; Escorihuela, M.J.; Font, J.; Reul, N.; Gruhier, C.; et al. The SMOS mission: New tool for monitoring key elements of the global water cycle. Proc. IEEE 2010, 98, 666–687. [Google Scholar] [CrossRef] [Green Version]
  19. De Jeu, R.; Holmes, T.R.H.; Owe, M. Deriving land surface parameters from 3 different vegetated sites with the ELBARA 1.4 GHz passive microwave radiometer. In Proceedings of the Remote Sensing of Agriculture, Ecosystems, and Hydrology V (SPIE), Barcelona, Spain, 8–12 September 2003; Volume 5232. [Google Scholar]
  20. Saleh, K.; Wigneron, J.-P.; Waldteufel, P.; de Rosnay, P.; Schwank, M.; Calvet, J.-C.; Kerr, Y.H. Estimates of Surface Soil Moisture under Grass Covers using L-band Radiometry. Remote Sens. Environ. 2007, 109, 42–53. [Google Scholar] [CrossRef]
  21. Wigneron, J.P.; Chanzy, A.; Kerr, Y.H.; Lawrence, H.; Shi, J.; Escorihuela, M.J.; Mironov, V.; Mialon, A.; Demontoux, F.; de Rosnay, P.; et al. Evaluating an Improved Parameterization of the Soil Emission in L-MEB. IEEE Trans. Geosci. Remote Sens. 2011, 49, 1177–1189. [Google Scholar] [CrossRef]
  22. Hofer, R.; Mätzler, C. Investigations on snow parameters by radiometry in the 3- to 60-mm wavelength region. J. Geophys. Res. 1980, 85, 453. [Google Scholar] [CrossRef]
  23. Mätzler, C.; Aebischer, H.; Schanda, E. Microwave dielectric properties of surface snow. IEEE J. Ocean. Eng. 1984, 9, 366–371. [Google Scholar] [CrossRef]
  24. Pulliainen, J.; Hallikainen, M. Retrieval of Regional Snow Water Equivalent from Space-Borne Passive Microwave Observations. Remote Sens. Environ. 2001, 75, 76–85. [Google Scholar] [CrossRef]
  25. Kaleschke, L.; Tian-Kunze, X.; Maaß, N.; Mäkynen, M.; Drusch, M. Sea ice thickness retrieval from SMOS brightness temperatures during the Arctic freeze-up period. Geophys. Res. Lett. 2012, 39. [Google Scholar] [CrossRef]
  26. Bogena, H.R.; Huisman, J.A.; Schilling, B.; Weuthen, A.; Vereecken, H. Effective calibration of low-cost soil water content sensors. Sensors 2017, 17, 208. [Google Scholar] [CrossRef] [PubMed]
  27. Schwank, M.; Wigneron, J.P.; Lopez-Baeza, E.; Völksch, I.; Mätzler, C.; Kerr, Y. L-Band Radiative Properties of Vine Vegetation at the SMOS Cal/Val Site MELBEX III. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1587–1601. [Google Scholar] [CrossRef] [Green Version]
  28. Mätzler, C.; Weber, D.; Wüthrich, M.; Schneeberger, K.; Stamm, C.; Wydler, H.; Fluhler, H. ELBARA, the ETH L-band radiometer for soil-moisture research. In Proceedings of the 2003 IEEE International Geoscience and Remote Sensing Symposium (IEEE Cat. No.03CH37477), Toulouse, France, 21–25 July 2003; Volume 5, pp. 3058–3060. [Google Scholar]
  29. Institue of Applied Physics. Handbook of the MORA 11.4 GHz Radiometer; Institue of Applied Physics: Bern, Switzerland, 1991. [Google Scholar]
  30. Guglielmetti, M.; Schwank, M.; Mätzler, C.; Oberdörster, C.; Vanderborght, J.; Flühler, H. FOSMEX: Forest Soil Moisture Experiments with Microwave Radiometry. IEEE Trans. Geosci. Remote Sens. 2008, 46, 727–735. [Google Scholar] [CrossRef]
  31. Werner, C.; Wiesmann, A.; Strozzi, T.; Schneebeli, M.; Matzler, C. The snowscat ground-based polarimetric scatterometer: Calibration and initial measurements from Davos Switzerland. In Proceedings of the 2010 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 25–30 July 2010; pp. 2363–2366. [Google Scholar]
  32. Pellarin, T.; Wigneron, J.-P.; Calvet, J.-C.; Berger, M.; Douville, H.; Ferrazzoli, P.; Kerr, Y.H.; Lopez-Baeza, E.; Pulliainen, J.; Simmonds, L.P. Two-year global simulation of L-band brightness temperatures over land. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2135–2139. [Google Scholar] [CrossRef]
  33. Loth, B.; Graf, H.-F.; Oberhuber, J.M. Snow cover model for global climate simulations. J. Geophys. Res. 1993, 98, 10451–10464. [Google Scholar] [CrossRef]
  34. Tarongi, J.M.; Camps, A. Normality Analysis for RFI Detection in Microwave Radiometry. Remote Sens. 2010, 2, 191–210. [Google Scholar] [CrossRef]
  35. Le Vine, D.M. ESTAR experience with RFI at L-band and implications for future passive microwave remote sensing from space. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Toronto, ON, Canada, 24–28 June 2002; pp. 847–849. [Google Scholar]
  36. Li, L.; Njoku, E.G.; Im, E.; Chang, P.S.; Germain, K.S. A preliminary survey of radio-frequency interference over the US in Aqua AMSR-E data. IEEE Trans. Geosci. Remote Sens. 2004, 42, 380–390. [Google Scholar] [CrossRef]
  37. Njoku, E.G.; Ashcroft, P.; Chan, T.K.; Li, L. Global survey and statistics of radio-frequency interference in AMSR-E land observations. IEEE Trans. Geosci. Remote Sens. 2005, 43, 938–947. [Google Scholar] [CrossRef]
  38. Massey, F.J., Jr. The Kolmogorov-Smirnov test for goodness of fit. J. Am. Stat. Assoc. 1951, 46, 68–78. [Google Scholar] [CrossRef]
  39. Lilliefors, H.W. On the Kolmogorov-Smirnov test for normality with mean and variance unknown. J. Am. Stat. Assoc. 1967, 62, 399–402. [Google Scholar] [CrossRef]
  40. Ruf, C.S.; Gross, S.M.; Misra, S. RFI detection and mitigation for microwave radiometry with an agile digital detector. IEEE Trans. Geosci. Remote Sens. 2006, 44, 694–706. [Google Scholar] [CrossRef]
  41. Roo, R.D.D.; Misra, S.; Ruf, C.S. Sensitivity of the Kurtosis Statistic as a Detector of Pulsed Sinusoidal RFI. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1938–1946. [Google Scholar] [CrossRef]
  42. Misra, S.; Ruf, C.S. Detection of Radio-Frequency Interference for the Aquarius Radiometer. IEEE Trans. Geosci. Remote Sens. 2008, 46, 3123–3128. [Google Scholar] [CrossRef]
  43. Johnson, J.T.; Potter, L.C. Performance Study of Algorithms for Detecting Pulsed Sinusoidal Interference in Microwave Radiometry. IEEE Trans. Geosci. Remote Sens. 2009, 47, 628–636. [Google Scholar] [CrossRef]
  44. Schwank, M.; Guglielmetti, M.; Mätzler, C.; Flühler, H. Testing a New Model for the L-band Radiation of Moist Leaf Litter. IEEE Trans. Geosci. Remote Sens. 2008, 46, 1982–1994. [Google Scholar] [CrossRef]
  45. Wigneron, J.-P.; Kerr, Y.; Waldteufel, P.; Saleh, K.; Richaume, P.; Ferrazzoli, P.; Escorihuela, M.-J.; Grant, J.P.; Hornbuckle, B.; de Rosnay, P.; et al. L-band microwave emission of the biosphere (L-MEB) model: Description and calibration against experimental data sets over crop fields. Remote Sens. Environ. 2007, 107, 639–655. [Google Scholar] [CrossRef]
  46. Mätzler, C.; Wiesmann, A. Extension of the microwave emission model of layered snowpacks to coarse-grained snow. Remote Sens. Environ. 1998, 70, 317–325. [Google Scholar] [CrossRef]
  47. Wigneron, J.-P.; Laguerre, L.; Kerr, Y.H. A simple parameterization of the L-band microwave emission from rough agricultural soils. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1697–1707. [Google Scholar] [CrossRef]
  48. Flores, A.N.; Entekhabi, D.; Bras, R.L.; Ivanov, V.Y. Hillslope-Scale Controls on Remote Sensing of Soil Moisture with Microwave Radiometry. In Proceedings of the 2008 IEEE International Geoscience and Remote Sensing Symposium, Boston, MA, USA, 7–11 July 2008; pp. II-695–II-698. [Google Scholar]
  49. Mätzler, C.; Wiesmann, A. Documentation for MEMLS, Version 3 ‘Microwave Emission Model of Layered Snowpacks’. MEMLS Manual, Version: 04.05.12. 2012. Available online: http://www.iap.unibe.ch/publications/pub-detail.php?lang=de&id=236 (accessed on 17 November 2017).
  50. Hallikainen, M.; Ulaby, F.; Dobson, M.; El-Rayes, M. Dielectric measurements of soils in the 3- to 37-GHz band between −50 °C and 23 °C. In Proceedings of the IGARSS’84: Remote Sensing, from Research Towards Operational Use, Strasbourg, France, 27–30 August 1984. [Google Scholar]
  51. Liebe, H.J.; Hufford, G.A.; Manabe, T. A model for the complex permittivity of water at frequencies below 1 THz. Int. J. Infrared Millim. Waves 1991, 12, 659–675. [Google Scholar] [CrossRef]
  52. Wigneron, J.-P.; Jackson, T.; O’Neill, P.; de Lannoy, G.; de Rosnay, P.; Walker, J.; Ferrazzoli, P.; Mironov, V.; Bircher, S.; Grant, J. Modelling the passive microwave signature from land surfaces: A review of recent results and application to the L-band SMOS & SMAP soil moisture retrieval algorithms. Remote Sens. Environ. 2017, 192, 238–262. [Google Scholar]
  53. Kerr, Y.H.; Waldteufel, P.; Richaume, P.; Wigneron, J.P.; Ferrazzoli, P.; Mahmoodi, A.; Al Bitar, A.; Cabot, F.; Gruhier, C.; Juglea, S.E.; et al. The SMOS Soil Moisture Retrieval Algorithm. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1384–1403. [Google Scholar]
  54. Entekhabi, D.; Njoku, E.G.; Neill, P.E.O.; Kellogg, K.H.; Crow, W.T.; Edelstein, W.N.; Entin, J.K.; Goodman, S.D.; Jackson, T.J.; Johnson, J.; et al. The Soil Moisture Active Passive (SMAP) Mission. Proc. IEEE 2010, 98, 704–716. [Google Scholar] [CrossRef]
Figure 1. Schematics of the footprint areas and the location of the in-situ sensors. ELBARA-II was initially installed at the center (position P1) of the upper-most platform of the tower. The radiometer scaffold was moved to position P2 (bottom right corner of the tower) on 12 December 2016 for RFI improvement.
Figure 1. Schematics of the footprint areas and the location of the in-situ sensors. ELBARA-II was initially installed at the center (position P1) of the upper-most platform of the tower. The radiometer scaffold was moved to position P2 (bottom right corner of the tower) on 12 December 2016 for RFI improvement.
Remotesensing 09 01185 g001
Figure 2. (a) Measured snowpack height h S ; and (b) average bottom-layer snow density ρ S over time. The snow cover quickly melted down in the second half of March 2017 and almost disappeared within approximately the last 10 day of the measurement campaign.
Figure 2. (a) Measured snowpack height h S ; and (b) average bottom-layer snow density ρ S over time. The snow cover quickly melted down in the second half of March 2017 and almost disappeared within approximately the last 10 day of the measurement campaign.
Remotesensing 09 01185 g002
Figure 3. NIR photos of the snowpack profile taken on (a) 9 January and (b) 27 February. In addition to an increased snow height, significantly more complex layering and more variable snow grain size and types can be observed in the second profile.
Figure 3. NIR photos of the snowpack profile taken on (a) 9 January and (b) 27 February. In addition to an increased snow height, significantly more complex layering and more variable snow grain size and types can be observed in the second profile.
Remotesensing 09 01185 g003
Figure 4. Panels (a,b) show the time series of in-situ measured ε G along transects 1 and 2 (shown in Figure 1), respectively. Panel (c) shows the average ground temperature T G measured by the 12 SMT-100 sensors along transects 1 and 2. Panel (d) indicates temperatures T air ,   T 15   cm , and T 50   cm measured by ELBARA-II’s PT-100 temperature sensor and SMT-100 sensors placed 15 cm and 50 cm above ground, respectively. Panel (e) shows the recorded precipitation (both rain and snow) in units of mm/10 min. over the entire campaign.
Figure 4. Panels (a,b) show the time series of in-situ measured ε G along transects 1 and 2 (shown in Figure 1), respectively. Panel (c) shows the average ground temperature T G measured by the 12 SMT-100 sensors along transects 1 and 2. Panel (d) indicates temperatures T air ,   T 15   cm , and T 50   cm measured by ELBARA-II’s PT-100 temperature sensor and SMT-100 sensors placed 15 cm and 50 cm above ground, respectively. Panel (e) shows the recorded precipitation (both rain and snow) in units of mm/10 min. over the entire campaign.
Remotesensing 09 01185 g004
Figure 5. Flowchart of the approach used to convert ELBARA-II raw data into calibrated T B p .
Figure 5. Flowchart of the approach used to convert ELBARA-II raw data into calibrated T B p .
Remotesensing 09 01185 g005
Figure 6. Four examples of measured sample voltage distributions P D F m ( U RMA , i V , 1 ) , and the corresponding Gaussian fits P D F Gauss ( U RMA , i V , 1 ) : (a) a “healthy” measurement; and (bd) distorted measurements.
Figure 6. Four examples of measured sample voltage distributions P D F m ( U RMA , i V , 1 ) , and the corresponding Gaussian fits P D F Gauss ( U RMA , i V , 1 ) : (a) a “healthy” measurement; and (bd) distorted measurements.
Remotesensing 09 01185 g006aRemotesensing 09 01185 g006b
Figure 7. Flowchart illustrating the course of action followed to estimate effective losses of transmission line (TL) and active cold source (ACS) noise temperatures.
Figure 7. Flowchart illustrating the course of action followed to estimate effective losses of transmission line (TL) and active cold source (ACS) noise temperatures.
Remotesensing 09 01185 g007
Figure 8. (a) Effective L TL p and (b) 〈 T ACS p , c h 〉 (ch = 1, 2 and p = H, V) computed with increasing numbers of sky measurements.
Figure 8. (a) Effective L TL p and (b) 〈 T ACS p , c h 〉 (ch = 1, 2 and p = H, V) computed with increasing numbers of sky measurements.
Remotesensing 09 01185 g008
Figure 9. (a) Block diagram of the L-band Specific Microwave Emission Model of Layered Snowpacks (“LS—MEMLS”) used to simulate L-band brightness temperatures T B p over snow covered grounds; and (b) sketch of the two-stream emission model (2S-EM) employed in “LS—MEMLS”. Symbols and model components are explained in the text.
Figure 9. (a) Block diagram of the L-band Specific Microwave Emission Model of Layered Snowpacks (“LS—MEMLS”) used to simulate L-band brightness temperatures T B p over snow covered grounds; and (b) sketch of the two-stream emission model (2S-EM) employed in “LS—MEMLS”. Symbols and model components are explained in the text.
Remotesensing 09 01185 g009
Figure 10. Sensitivities of brightness temperatures with respect to snow liquid water column W C S for snow mass density ρ S = 300   kg   m 3 : (a) “uniform” snow; (b) snow on “top” of the snowpack; (c) moist snow “sandwiched” in-between dry snow; and (d) a moist “bottom” snow-layer. T B , R p (thin lines) and T B , N p (bold lines) are for snow atop the reflector (R) and atop the natural (N) frozen ground with ε G = 5 , respectively. Polarization p = H (solid), p = V (dashed), and the nadir angles θ = 0 ° (black), θ = 30 ° (red), and θ = 60 ° (green).
Figure 10. Sensitivities of brightness temperatures with respect to snow liquid water column W C S for snow mass density ρ S = 300   kg   m 3 : (a) “uniform” snow; (b) snow on “top” of the snowpack; (c) moist snow “sandwiched” in-between dry snow; and (d) a moist “bottom” snow-layer. T B , R p (thin lines) and T B , N p (bold lines) are for snow atop the reflector (R) and atop the natural (N) frozen ground with ε G = 5 , respectively. Polarization p = H (solid), p = V (dashed), and the nadir angles θ = 0 ° (black), θ = 30 ° (red), and θ = 60 ° (green).
Remotesensing 09 01185 g010
Figure 11. T B , N p (blue) and T B , R p (magenta) of brightness temperatures at: (a) horizontal (p = H); and (b) vertical (p = V) polarization measured over “natural” and “reflector” areas, respectively. Panels (c,d) show air temperature T air and precipitation rate, respectively. The vertical dashed lines delimit the “snow-free period” (before 3 January), the “cold winter period” (3–31 January), and the “early spring period” (after 31 January).
Figure 11. T B , N p (blue) and T B , R p (magenta) of brightness temperatures at: (a) horizontal (p = H); and (b) vertical (p = V) polarization measured over “natural” and “reflector” areas, respectively. Panels (c,d) show air temperature T air and precipitation rate, respectively. The vertical dashed lines delimit the “snow-free period” (before 3 January), the “cold winter period” (3–31 January), and the “early spring period” (after 31 January).
Remotesensing 09 01185 g011

Share and Cite

MDPI and ACS Style

Naderpour, R.; Schwank, M.; Mätzler, C. Davos-Laret Remote Sensing Field Laboratory: 2016/2017 Winter Season L-Band Measurements Data-Processing and Analysis. Remote Sens. 2017, 9, 1185. https://doi.org/10.3390/rs9111185

AMA Style

Naderpour R, Schwank M, Mätzler C. Davos-Laret Remote Sensing Field Laboratory: 2016/2017 Winter Season L-Band Measurements Data-Processing and Analysis. Remote Sensing. 2017; 9(11):1185. https://doi.org/10.3390/rs9111185

Chicago/Turabian Style

Naderpour, Reza, Mike Schwank, and Christian Mätzler. 2017. "Davos-Laret Remote Sensing Field Laboratory: 2016/2017 Winter Season L-Band Measurements Data-Processing and Analysis" Remote Sensing 9, no. 11: 1185. https://doi.org/10.3390/rs9111185

APA Style

Naderpour, R., Schwank, M., & Mätzler, C. (2017). Davos-Laret Remote Sensing Field Laboratory: 2016/2017 Winter Season L-Band Measurements Data-Processing and Analysis. Remote Sensing, 9(11), 1185. https://doi.org/10.3390/rs9111185

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