Observed and modelled effects of ice lens formation on passive microwave brightness temperatures over snow covered tundra

Immediately before an April 2007 snow survey and passive microwave radiometer field campaign in the Northwest Territories, Canada, a rain-on-snow event deposited a thin (~ 3 mm) continuous layer of ice on the surface of the snowpack. At eight sites the brightness temperature (T_b) of the undisturbed snow pack was measured with a multi-frequency dual polarization (6.9, 19, 37, and 89 GHz) ground based radiometer system. The ice lens was then carefully removed and the T_bs were measured again. The individual V-pol channels and the 37 V − 19 V difference were largely unaffected by the presence of the ice lens, exhibiting a systematic shift of about 3 K. In comparison, the ice lens had a considerable effect on the H-pol T_b at all frequencies, with a mean difference (ice lens present − ice lens removed) of − 9 K (± 5.3 K) at 6.9 GHz, − 40 K (± 11.3 K) at 19 GHz, − 33 K (± 7.6 K) at 37 GHz, and − 19 K (± 8.0 K) at 89 GHz. The effect of the ice lens on H-pol measurements was also observed with spaceborne data from the Advanced Microwave Scanning Radiometer (AMSR-E) satellite data.Simulations of T_b were produced for each site using a new two layer formulation of the Helsinki University of Technology (HUT) snow emission model. The ice lens was used as the top layer and the underlying snowpack considered as a homogenous second layer. The agreement between observations and simulations was variable, with agreement strongest at 19 GHz. A comparison with simulations produced using the Microwave Emission Model of Layered Snowpacks (MEMLS) suggests HUT model uncertainty is related not to the ice lens, but to difficulties in simulating emission from deep snow. Overall, the observations and simulations suggest H-pol measurements are capable of detecting new ice layers across the tundra snowpack, while V-pol measurements are more appropriate for snow water equivalent (SWE) retrievals due to their relative insensitivity to ice layers.     

AMSR-E data inversion for soil temperature estimation under snow cover

Climate warming is the focus of several studies where the soil temperature plays an essential role as a state variable for the surface energy balance of the Earth. Many methods have been developed to determine summer surface temperature, but the determination in presence of snow is an ill-conditioned problem for microwave techniques because snow changes the emissivity of the surface. This project aims to improve the estimation of soil temperature, within the top 5 cm of the ground, under the snowpack using passive microwave remote sensing. Results show the potential of the passive microwave brightness temperature inversion at 10 GHz (derived from the Advanced Microwave Scanning Radiometer—Earth Observing System, AMSR-E) for the estimation of soil temperature using a physical multilayer snow-soil model (SNTHERM) coupled with a snow emission model (HUT). The snow model is driven with meteorological measurements from ground-based stations as well as data generated from reanalysis. The proposed iterative retrieval method minimizes the difference between the simulated and measured brightness temperature using the soil temperature as a free parameter given by SNTHERM. Results are validated against ground-based measurements at several sites across Canada through several winter seasons. The overall root mean square error and bias in the retrieved soil temperature is respectively 3.29 K and 0.56 K, lower than the error derived from the snow-soil model without the use of remote sensing. The accuracy in detection of frozen/unfrozen soil under the snowpack is 78%, which is improved up to 81% if the spring melting period is not considered. This original procedure constitutes a very promising tool to characterize the soil (frozen or not) under snow cover, as well as its evolution in northern remote locations where measurements are unavailable.     

Development of a winter snow water equivalent algorithm using in situ passive microwave radiometry over snow-covered first-year sea ice

A snow water equivalent (SWE) algorithm has been developed for thin and thick snow using both in situ microwave measurements and snow thermophysical properties, collected over landfast snow covered first-year sea ice during the Canadian Arctic Shelf Exchange Study (CASES) overwintering mission from December 2003 to May 2004. Results showed that the behavior of brightness temperatures (T_bs) in thin snow covers was very different from those in a thick snowpack. Microwave SWE retrievals using the combination of T_b 19 GHz and air temperature (multiple regression) over thick snow are quite accurate, and showed very good agreement with the physical data (R² = 0.94) especially during the cooling period (i.e., from freeze up to the minimum air temperature recorded) where the snow is dry and cold. Thin snow SWE predictions also showed fairly good agreement with field data (R² = 0.70) during the cold season. The differences between retrieved and in situ SWE for both thin and thick snow cover are mainly attributable to the variations in air temperature, snow wetness and spatial heterogeneity in snow thickness.     

Hourly simulations of the microwave brightness temperature of seasonal snow in Quebec, Canada, using a coupled snow evolution-emission model

To interpret the snowpack evolution, and in particular to estimate snow water equivalent (SWE), passive microwave remote sensing has proved to be a useful tool given its sensitivity to snow properties. However, the main uncertainties using existing SWE algorithms arise from snow metamorphism which evolves during the winter season, and changes the snow emissivity. To consider the evolution in snow emissivity a coupled snow evolution-emission model can be used to simulate the brightness temperature (T_B) of the snowpack.During a dedicated campaign in the winter season, from November to April, of 2007-2008 two surface-based radiometers operating at 19 GHz and 37 GHz continuously measured the passive microwave radiation emitted through a seasonal snowpack in southern Quebec (Canada). This paper aims at modeling and interpreting this time series of T_B over the whole season, with an hourly step, using a coupled multi-layer snow evolution-emission model. The thermodynamic snow evolution model, referred as to Crocus, was driven by local meteorological measurements. Results from this model provided, in turn, the input variables to run the Microwave Emission Model of Layered Snowpacks (MEMLS) in order to predict T_B at 19 GHz and 37 GHz for both vertical (V) and horizontal (H) polarizations. The accuracy of T_B predicted by the Crocus-MEMLS coupled model was evaluated using continuous measurements from the surface-based radiometers.The weather conditions observed during the winter season were diverse, including several warm periods with melting snow and rain-on-snow events, producing very complex variations in the time series of T_B. To aid our analysis, we identified days with melting snow versus days with dry snow. The Crocus-MEMLS coupled model was able to accurately predict melt events with a success rate of 86%. The residual error was due to an overestimation of the duration of several melt events simulated by Crocus. This problem was explained by 1) a limitation of percolation, and 2) a very long-acting melt of lower layers due to geothermal flux.When the snowpack was completely dry, the global trend of T_B during the season was characterized by a decrease of T_B due to growth in the snow grain size. During most of the season, Crocus-MEMLS correctly predicted the evolution of T_B resulting from temperature gradient metamorphism; the root mean square errors ranged between 2.8 K for the 19 GHz vertical polarization (19V) and 6.9 K for the 37 GHz horizontal polarization (37H). However, during dry periods near the end of the season, the values of T_B were strongly overestimated. This overestimation was mainly due to a limitation of the growth of large snow grains in the wet snowpack simulated by Crocus. This effect was confirmed by estimating snow grain sizes from the observed T_B and the coupled model. The estimated snow grain sizes were larger and more realistic than those initially predicted by the Crocus model.     

Detection of small-scale roughness and refractive index of sea ice in passive satellite microwave remote sensing

Polar ice masses and sheets are sensitive indicators of climate change. Small-scale surface roughness significantly impacts the microwave emission of the sea ice/snow surface; however, published results of surface roughness measurements of sea ice are rare. Knowing the refractive index is important to discriminate between objects. In this study, the small-scale roughness and refractive index over sea ice are estimated with AMSR-E observations and a unique method. Consequently, the small-scale surface roughness of 0.25 cm to 0.5 cm at AMSR-E 6.9 GHz shows reasonable agreement with the results of known observations, ranging from 0.2 cm to 0.6 cm for the sea ice in the Antarctic and Arctic regions. The refractive indexes are retrieved from 1.6 to 1.8 for winter, from 1.2 to 1.4 for summer in the Arctic and the Antarctic, which are similar to those of the sea ice and results from previous studies. This research shows the physical characteristics of the sea ice edges and melting process. Accordingly, this investigation provides an effective procedure for retrieving the small-scale roughness and refractive index of sea ice and snow. Another advantage of this study is the ability to distinguish sea ice from the sea surface by their relative small-scale roughness.     

Development of a tundra-specific snow water equivalent retrieval algorithm for satellite passive microwave data

Airborne and satellite brightness temperature (T_B) measurements were combined with intensive field observations of sub-Arctic tundra snow cover to develop the framework for a new tundra-specific passive microwave snow water equivalent (SWE) retrieval algorithm. The dense snowpack and high sub-grid lake fraction across the tundra mean that conventional brightness temperature difference approaches (such as the commonly used 37 GHz-19 GHz) are not appropriate across the sub-Arctic. Airborne radiometer measurements (with footprint dimensions of approximately 70 × 120 m) acquired across sub-Arctic Canada during three field campaigns during the 2008 winter season were utilized to illustrate a slope reversal in the 37 GHz T_B versus SWE relationship. Scattering by the tundra snowpack drives a negative relationship until a threshold SWE value is reached near 130 mm at which point emission from the snowpack creates a positive but noisier relationship between 37 GHz T_B and SWE.The change from snowpack scattering to emission was also evident in the temporal evolution of 37 GHz T_B observed from satellite measurements. AMSR-E brightness temperatures (2002/03-2006/07) consistently exhibited decreases through the winter before reaching a minimum in February or March, followed by an increase for weeks or months before melt. The cumulative absolute change (Σ|Δ37V|) in vertically polarized 37 GHz T_B was computed at both monthly and pentad intervals from a January 1 start date and compared to ground measured SWE from intensive and regional snow survey campaigns, and climate station observations. A greater (lower) cumulative change in |Δ37V| was significantly related to greater (lower) ground measured SWE (r² = 0.77 with monthly averages; r² = 0.67 with pentad averages). Σ|Δ37V| was only weakly correlated with lake fraction: monthly r² values calculated for January through April 2003-2007 were largely less than 0.2. These results indicate that this is a computationally straightforward and viable algorithmic framework for producing tundra-specific SWE datasets from the complete satellite passive microwave record (1979 to present).     

Retrieval of snow grain size over Greenland from MODIS

This paper presents a new automatic algorithm to derive optical snow grain size at 1 km resolution using Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. The retrieval is conceptually based on an analytical asymptotic radiative transfer model which predicts spectral bidirectional snow reflectance as a function of the grain size and ice absorption. The snow grains are modeled as fractal rather than spherical particles in order to account for their irregular shape. The analytical form of solution leads to an explicit and fast retrieval algorithm. The time series analysis of derived grain size shows a good sensitivity to snow melting and snow precipitation events. Pre-processing is performed by a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which includes gridding MODIS data to 1 km resolution, water vapor retrieval, cloud masking and an atmospheric correction. MAIAC cloud mask is a new algorithm based on a time series of gridded MODIS measurements and an image-based rather than pixel-based processing. Extensive processing of MODIS TERRA data over Greenland shows a robust discrimination of clouds over bright snow and ice. Because in-situ grain size measurements over Greenland were not available at the time of this work, the validation was performed using data of Aoki et al. (Aoki, T., Hori, M., Motoyoshi, H., Tanikawa, T., Hachikubo, A., Sugiura, K., et al. (2007). ADEOS-II/GLI snow/ice products — Part II: Validation results using GLI and MODIS data. Remote Sensing of Environment, 111, 274-290) collected at Barrow (Alaska, USA), and Saroma, Abashiri and Nakashibetsu (Japan) in 2001-2005. The retrievals correlate well with measurements in the range of radii ~ 0.1-1 mm, although retrieved optical diameter may be about a factor of 1.5 lower than the physical measured diameter. As part of validation analysis for Greenland, the derived grain size from MODIS over selected sites in 2004 was compared to the microwave brightness temperature measurements of SSM/I radiometer which is sensitive to the amount of liquid water in the snowpack. The comparison showed a good qualitative agreement, with both datasets detecting two main periods of snowmelt. Additionally, MODIS grain size was compared with predictions of the snow model CROCUS driven by measurements of the automatic weather stations of the Greenland Climate Network. We found that the MODIS value is on average a factor of two smaller than CROCUS grain size. This result agrees with the direct validation analysis indicating that the snow reflectance model may need a “calibration” factor of ~ 1.5 for the retrieved grain size to match the physical snow grain size. Overall, the agreement between CROCUS and MODIS results was satisfactory, in particular before and during the first melting period in mid-June. Following detailed time series analysis of snow grain size for four permanent sites, the paper presents maps of this important parameter over the Greenland ice sheet for the March-September period of 2004.     

Merging complementary remote sensing datasets in the context of snow water equivalent reconstruction

Time series of snow covered area (SCA) estimates from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat Enhanced Thematic Mapper (ETM+) were merged with a spatially explicit snowmelt model to reconstruct snow water equivalent (SWE) in the Rio Grande headwaters (3419 km²). A linear optimization scheme was used to derive SCA estimates that preserve the statistical moments of the higher spatial resolution (i.e. 30 m) ETM+ data and resolve the superior temporal signal (i.e. ∼ daily) of the MODIS data. It was found that merging the two SCA products led to an 8% decrease and an 18% increase in the basinwide SWE in 2001 and 2002, respectively, compared to the SWE estimated from ETM+ only. Relative to SWE simulations using only ETM+ data, the hybrid SCA estimates reduced the mean absolute SWE error by 17 and 84% in 2001 and 2002, respectively; errors were determined using intensive snow survey data and two separate methods of scaling snow survey field measurements of SWE to the 1-km model pixel resolution. SWE bias for both years was reduced by 49% and skewness was reduced from − 0.78 to 0.49. These results indicate that the hybrid SWE was closer to being an unbiased estimate of the measured SWE and errors were distributed more normally. The accuracy of the SCA estimates is likely dependent on the vegetation fraction.     

一种基于雪粒径演化过程的积雪亮温模拟新方法

本研究采用HUT模型、DMRT模型和MEMLS模型模拟积雪雪粒子与不同波段（18.7 GHz和36.5 GHz）微波相互作用（吸收和消光）,并用于辐射传输模型。而雪粒径的获取一直是一个难点,本研究由Jordan91雪粒径演化模型演化得到雪粒径,并将其作为辐射传输模型的输入参数,基于像元内实测数据进行混合像元18.7和36.5 GHz水平极化亮温模拟。结果表明：采用HUT模型、DMRT模型和MEMLS模型的消光系数在18.7 GHz时模拟亮温的偏差分别为-3.6、-1.8和-0.7 K,在36.5 GHz时分别为4.0、10.4和14.4 K。对于18.7 GHz水平极化和36.5 GHz水平极化,基于有效雪粒径的亮温模拟与基于雪粒径演化过程的亮温模拟精度呈现出很好的线性关系。因此,基于雪粒径演化过程的方法是一种合适的获取辐射传输模型中雪粒径参数的方法。     

A New Method of Simulating Bright Temperature of Snow Cover based on Snow Grain Size Evolution Process

This research used HUT model, DMRT model and MEMLS model to simulate interactions(absorption and extinction) between snow grainsfor different wave bands (18.7 GHz and 36.5 GHz) of microwave which were used for radiative transfer model. Obtaining the snow grain size is always a difficulty. So this research used Jordan91 snow grain size evolution model to evolve snow grain size which was regarded as input parameter of radiative transfer model, and used measured data to simulate spaceborne brightness temperature for 18.7 GHz horizontal polarization and 36.5 GHz horizontal polarization in a mixed pixel. The results showed that the bias of simulation brightness temperature using extinction coefficient of HUT model, DMRT model and MEMLS model for 18.7 GHz horizontal polarization were -3.6 K、-1.8 K and -0.7 K respectively, and for 36.5 GHz horizontal polarization were 4.0 K、10.4 K and 14.4 K respectively. For 18.7 GHz horizontal polarization and 36.5 GHz horizontal polarization, the bright temperature simulation based on effective snow grain size shows a good linear relationship with the brightness temperature simulation basedon snow grain size evolution process. Therefore, the method based on the snow grain size evolution process is a suitable method for obtaining the snow grain size parameter in the radiative transfer model.     

Accuracy assessment of the MODIS 16-day albedo product for snow: comparisons with Greenland in situ measurements

The accuracy of the Moderate Resolution Imaging Spectroradiometer (MODIS) 16-day albedo product (MOD43) is assessed using ground-based albedo observations from automatic weather stations (AWS) over spatially homogeneous snow and semihomogeneous ice-covered surfaces on the Greenland ice sheet. Data from 16 AWS locations, spanning the years 2000-2003, were used for this assessment. In situ reflected shortwave data were corrected for a systematic positive spectral sensitivity bias of between 0.01 and 0.09 on a site-by-site basis using precise optical black radiometer data. Results indicate that the MOD43 albedo product retrieves snow albedo with an average root mean square error (RMSE) of ±0.07 as compared to the station measurements, which have ±0.035 RMSE uncertainty. If we eliminate all satellite retrievals that rely on the backup algorithm and consider only the highest quality results from the primary bidirectional reflectance distribution function (BRDF) algorithm, the MODIS albedo RMSE is ±0.04, slightly larger than the in situ measurement uncertainty. There is general agreement between MODIS and in situ observations for albedo <0.7, while near the upper limit, a −0.05 MODIS albedo bias is evident from the scatter of the 16-site composite.     

Semi-empirical approach for estimating broadband albedo of snow

Snow is a medium that exhibits highly anisotropic reflectance throughout the solar spectrum. The anisotropic nature of snow shows more variability in snow metamorphic processes for wavelengths beyond 1.0 μm than in the visible spectrum. This behavior poses challenges for the development of a model that can retrieve broadband albedo from reflectance measurements throughout the snow season. In this paper, a semi-empirical model is presented to estimate near infrared (0.8-2.5 μm) albedo of snow. This model estimates spectral albedo at a wavelength of 1.240 μm using only three variables: solar zenith angle, scattering angle and measured reflectance, which is used to retrieve near infrared albedo. To form a base for such a model, quantification of reflectance patterns and variability in varying snow condition, i.e. snow grain size, and sun-sensor geometry are prerequisite. In this study the DIScrete Ordinate Radiative Transfer (DISORT) model is used to simulate bi-directional reflectance. The performance of the developed model is evaluated by using DISORT simulated spectral albedo for various snow grain sizes and solar zenith angles, as well as the Moderate Resolution Imaging Spectroradiometer (MODIS) and in-situ measurements. The developed model is shown to be capable of estimating spectral albedo at 1.240 μm with acceptable accuracy. The mean error (ME), mean absolute error (MAE), and root mean squared error (RMSE) in the estimates are found to be 0.053, 0.055 and 0.064, respectively, for a wide range of sun-sensor geometries and snow grain sizes. The model shows better accuracy for spectral albedo estimates than for those computed using the Lambertian reflectance assumption for snow, reducing the error in the range and standard deviation by 75% and 65%, respectively. Applying the model to MODIS, the retrieved albedo is found to be in good quantitative agreement (r = 0.82) with in-situ measurements. These improvements in albedo estimation should allow more accurate use of remote sensing measurements in climate and hydrological models.     

Validation of the Climate-SAF surface broadband albedo product: Comparisons with in situ observations over Greenland and the ice-covered Arctic Ocean

This paper describes a validation study performed by comparing the Climate-SAF Surface Albedo Product (SAL) to ground truth observations over Greenland and the ice-covered Arctic Ocean. We compare Advanced Very High Resolution Radiometer (AVHRR)-based albedo retrievals to data from the Greenland Climate Network (GCN) weather stations and the floating ice station Tara for polar summer 2007. The AVHRR dataset consists of 2755 overpasses. The overpasses are matched to in situ observations spatially and temporally. The SAL algorithm presented here derives the surface broadband albedo from AVHRR channels 1 and 2 using an atmospheric correction, temporal sampling of an empirical Bidirectional Reflectance Distribution Function (BRDF), and a narrow-to-broadband conversion algorithm. The satellite product contains algorithms for snow, sea ice, vegetation, bare soil, and water albedo. At the Summit and DYE-2 stations on the Greenland ice sheet, instantaneous SAL RMSE is 0.073. The heterogeneous surface conditions at satellite pixel scale over the stations near the Greenland west coast increase RMSE to > 0.12. Over Tara, the instantaneous SAL RMSE is 0.069. The BRDF sampling approach reduces RMSE over the ice sheet to 0.053, and to 0.045 over Tara. Taking into account various sources of uncertainty for both satellite retrievals and in situ observations, we conclude that SAL agrees with in situ observations within their limits of accuracy and spatial representativeness.     

Artificial neural network-based techniques for the retrieval of SWE and snow depth from SSM/I data

The retrieval of snow water equivalent (SWE) and snow depth is performed by inverting Special Sensor Microwave Imager (SSM/I) brightness temperatures at 19 and 37 GHz using artificial neural network ANN-based techniques. The SSM/I used data, which consist of Pathfinder Daily EASE-Grid brightness temperatures, were supplied by the National Snow and Ice Data Centre (NSIDC). They were gathered during the period of time included between the beginning of 1996 and the end of 1999 all over Finland. A ground snow data set based on observations of the Finnish Environment Institute (SYKE) and the Finnish Meteorological Institute (FMI) was used to estimate the performances of the technique. The ANN results were confronted with those obtained using the spectral polarization difference (SPD) algorithm, the HUT model-based iterative inversion and the Chang algorithm, by comparing the RMSE, the R², and the regression coefficients. In general, it was observed that the results obtained through ANN-based technique are better than, or comparable to, those obtained through other approaches, when trained with simulated data. Performances were very good when the ANN were trained with experimental data.     

Estimating northern hemisphere snow water equivalent for climate research through assimilation of space-borne radiometer data and ground-based measurements

The key variable describing global seasonal snow cover is snow water equivalent (SWE). However, reliable information on the hemispheric scale variability of SWE is lacking because traditional methods such as interpolation of ground-based measurements and stand-alone algorithms applied to space-borne observations are highly uncertain with respect to the spatial distribution of snow mass and its evolution. In this paper, an algorithm assimilating synoptic weather station data on snow depth with satellite passive microwave radiometer data is applied to produce a 30-year-long time-series of seasonal SWE for the northern hemisphere. This data set is validated using independent SWE reference data from Russia, the former Soviet Union, Finland and Canada. The validation of SWE time-series indicates overall strong retrieval performance with root mean square errors below 40 mm for cases when SWE < 150 mm. Retrieval uncertainty increases when SWE is above this threshold. The SWE estimates are also compared with results obtained by a typical stand-alone satellite passive microwave algorithm. This comparison demonstrates the benefits of the newly developed assimilation approach. Additionally, the trends and inter-annual variability of northern hemisphere snow mass during the era of satellite passive microwave measurements are shown.     

MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland

In situ measurements of snow albedo at five stations along a north-south transect in the dry-snow facies of the interior of Greenland follow the theoretically expected dependence of snow albedo with solar zenith angle (SZA). Greenland Climate Network (GC-Net) measurements from 1997 through 2007 exhibit the trend of modest surface brightening with increasing SZA on both diurnal and seasonal timescales. SZA explains up to 50% of seasonal albedo variability. The two other environmental factors considered, temperature and cloudiness, play much less significant roles in seasonal albedo variability at the five stations studied. Compared to the 10-year record of these GC-Net measurements, the five-year record of MODIS satellite-retrieved snow albedo shows a systematic negative bias for SZA larger than about 55°. Larger bias of MODIS snow albedo exists at more northerly stations. MODIS albedos successfully capture the snow albedo dependence on SZA and have a relatively good agreement with GC-Net measurements for SZA < 55°. The discrepancy of MODIS albedo with in situ albedo and with theory is determined mainly by two related factors, SZA and retrieval quality. While the spatiotemporal structure, especially zonal features, of the MODIS-retrieved albedo may be correct for large SZA, the accuracy deteriorates for SZA > 55° and often becomes physically unrealistic for SZA > 65°. This unphysical behavior biases parameterizations of surface albedo and restricts the range of usefulness of the MODIS albedo products. Seasonal-to-interannual trends in surface brightness in Greenland, and in polar (i.e., large SZA) regions in general, and model simulations of these trends, should be evaluated in light of these limitations.     

Utilizing SAR imagery and aquatic vegetation to map fresh and brackish lakes in the Brazilian Pantanal wetland

South America's Pantanal, the world's largest tropical wetland, contains hundreds of thousands of geochemically diverse lakes, ranging from dilute to brackish to saline in composition. These lakes form the backbone of the habitat that supports the highly diverse flora and fauna of the Pantanal, yet the natural processes that create and destroy them are largely unknown. The quantities and types of lakes in the Pantanal and their spatial distribution are therefore essential, but missing information required to understand the dynamics of the Pantanal ecosystem.RADARSAT S1 and S7, and JERS-1 imagery were integrated with field measurements of water geochemistry and characteristics of emergent aquatic vegetation for fresh and brackish lakes of the Nhecolândia region of the Brazilian Pantanal. A supervised classification was used to classify forest, pasture, bare soil, and lakes. A mask is then applied to produce an image of only lakes. The radar backscattering values were found to have a strong relationship with the emergent aquatic plant assemblages of the lakes—S1 imagery was the most useful. The plant assemblages, in turn, were observed to be strongly controlled by the total dissolved solids (TDS) and pH of the lakes. The relationships between backscattering values, plant assemblages, and geochemistry were then exploited to map the type and distribution of the lakes in the study area.Threshold rules were used to perform Level 1 and Level 2 classifications of the lakes. For the Level 1 classification, the σ^o values of RADARSAT S1 effectively separated brackish (10,000 > TDS > 1000 mg/kg) from fresh water lakes (TDS < 1000 mg/kg) with a total accuracy of 91%. For the Level 2 classification, the σ^o values of RADARSAT S1 effectively separated lakes into three geochemical groups: brackish (10,000 > TDS > 1000 mg/kg), hard with only Typha (1000 > TDS > 100 mg/kg), and fresh water lakes (TDS mg/kg < 100 mg/kg) with a total accuracy of 83%. Considering that the area is very remote and the lakes are very numerous, this may be the most feasible way to map lake type in the Pantanal.     

Comparison of Envisat radar and airborne laser altimeter measurements over Arctic sea ice

Sea ice thickness is a crucial, but very undersampled cryospheric parameter of fundamental importance for climate modeling. Advances in satellite altimetry have enabled the measurement of sea ice freeboard using satellite microwave altimeters. Unfortunately, validation of these new techniques has suffered from a lack of ground truth measurements. Therefore, an airborne campaign was carried out in March 2006 using laser altimetry and photo imagery to validate sea ice elevation measurements derived from the Envisat/RA-2 microwave altimeter.We present a comparative analysis of Envisat/RA-2 sea ice elevation processing with collocated airborne measurements collected north of the Canadian Archipelago. Consistent overall relationships between block-averaged airborne laser and Envisat elevations are found, over both leads and floes, along the full 1300 km aircraft track. The fine resolution of the airborne laser altimeter data is exploited to evaluate elevation variability within the RA-2 ground footprint. Our analysis shows good agreement between RA-2 derived sea ice elevations and those measured by airborne laser altimetry, particularly over refrozen leads where the overall mean difference is about 1 cm. Notwithstanding this small 1 cm mean difference, we identify a larger elevation uncertainty (of order 10 cm) associated with the uncertain location of dominant radar targets within the particular RA-2 footprint. Sources of measurement uncertainty or ambiguity are identified, and include snow accumulation, tracking noise, and the limited coverage of airborne measurements.     

Comparison between a multi-scattering and multi-layer snow scattering model and its parameterized snow backscattering model

Snow Water Equivalent (SWE) is a crucial parameter in the study of climatology and hydrology. Active microwave remote sensing is one of the most promising techniques for estimating the distribution of SWE at high spatial resolutions in large areas. Development of reliable and accurate inversion techniques to recover SWE is one of the most important tasks in current microwave researches. However, a number of snow pack properties, including snow density, particle size, crystal shape, stratification, ground surface roughness and soil moisture, affect the microwave scattering signals and need to be properly modeled and exploited. In this paper, we developed a multi-layer, multi-scattering model for dry snow based on recent theoretical advances in snow and surface modeling. In the proposed multi-layer model, Matrix Doubling method is used to account for scattering from each snow layer; and Advanced Integral Equation Model (AIEM) is incorporated into the model to describe surface scattering. Comparisons were made between the model predictions and field observations from NASA Cold Land Processes Field Experiment (CLPX) during Third Intensive Observation Period (IOP3) and SARALPS-2007 field experiment supported by ESA. The results indicated that model predictions were in good agreement with field observations. With the confirmed confidence, the analyses on multiple scattering, scatterer shape, and snow stratification effects were further made based on the model simulations. Furthermore, a parameterized snow backscattering model with a simple form and high computational efficiency was developed using a database generated by the multiple-scattering model. For a wide range of snow and soil properties, this parameterized model agrees well with the multiple-scattering model, with the root mean square error 0.20 dB, 0.24 dB and 0.43 dB for VV, HH and VH polarizations, respectively. This simplified model can be useful for the development of SWE retrieval algorithm and for fast simulations of radar signals over snow cover in land data assimilation systems.