Quantifying the uncertainty in passive microwave snow water equivalent observations

Passive microwave sensors (PM) onboard satellites have the capability to provide global snow observations which are not affected by cloudiness and night condition (except when precipitating events are occurring). Furthermore, they provide information on snow mass, i.e., snow water equivalent (SWE), which is critically important for hydrological modeling and water resource management. However, the errors associated with the passive microwave measurements of SWE are well known but have not been adequately quantified thus far. Understanding these errors is important for correct interpretation of remotely sensed SWE and successful assimilation of such observations into numerical models.This study uses a novel approach to quantify these errors by taking into account various factors that impact passive microwave responses from snow in various climatic/geographic regions. Among these factors are vegetation cover (particularly forest cover), snow morphology (crystal size), and errors related to brightness temperature calibration. A time-evolving retrieval algorithm that considers the evolution of snow crystals is formulated. An error model is developed based on the standard error estimation theory. This new algorithm and error estimation method is applied to the passive microwave data from Special Sensor Microwave/Imager (SSM/I) during the 1990-1991 snow season to produce annotated error maps for North America. The algorithm has been validated for seven snow seasons (from 1988 to 1995) in taiga, tundra, alpine, prairie, and maritime regions of Canada using in situ SWE data from the Meteorological Service of Canada (MSC) and satellite passive microwave observations. An ongoing study is applying this methodology to passive microwave measurements from Scanning Multichannel Microwave Radiometer (SMMR); future study will further refine and extend the analysis globally, and produce an improved SWE dataset of more than 25 years in length by combining SSMR and SSM/I measurements.     

Snowpack and runoff generation using AMSR-E passive microwave observations in the Upper Helmand Watershed, Afghanistan

Passive microwave estimates of snow water equivalent (SWE) were examined to determine their usefulness for evaluating water resources in the remote Upper Helmand Watershed, central Afghanistan. SWE estimates from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and the Special Sensor Microwave/Imager (SSM/I) passive microwave data were analyzed for six winter seasons, 2004-2009. A second, independent estimate of SWE was calculated for these same time periods using a hydrologic model of the watershed with a temperature index snow model driven using the Tropical Rainfall Measuring Mission (TRMM) gridded estimates of precipitation. The results demonstrate that passive microwave SWE values from SSM/I and AMSR-E are comparable. The AMSR-E sensor had improved performance in the early winter and late spring, which suggests that AMSR-E is better at detecting shallow snowpacks than SSM/I. The timing and magnitude of SWE values from the snow model and the passive microwave observations were sometimes similar with a correlation of 0.53 and accuracy between 55 and 62%. However, the modeled SWE was much lower than the AMSR-E SWE during two winter seasons in which TRMM data estimated lower than normal precipitation. Modeled runoff and reservoir storage predictions improved significantly when peak AMSR-E SWE values were used to update the snow model state during these periods. Rapid decreases in passive microwave SWE during precipitation events were also well aligned with flood flows that increased base flows by 170 and 940%. This finding supports previous northern latitude studies which indicate that the passive microwave signal's lack of scattering can be used to detect snow melt. The current study's extension to rain on snow events suggests an opportunity for added value for flood forecasting.     

Correcting for the influence of frozen lakes in satellite microwave radiometer observations through application of a microwave emission model

The spatial resolution of passive microwave observations from space is of the order of tens of kilometers with currently available instruments, such as the Special Sensor Microwave/Imager (SSM/I) and Advanced Microwave Scanning Radiometer (AMSR-E). The large field of view of these instruments dictates that the observed brightness temperature can originate from heterogeneous land cover, with different vegetation and surface properties.In this study, we assess the influence of freshwater lakes on the observed brightness temperature of AMSR-E in winter conditions. The study focuses on the geographic region of Finland, where lakes account for 10% of the total terrestrial area. We present a method to mitigate for the influence of lakes through forward modeling of snow covered lakes, as a part of a microwave emission simulation scheme of space-borne observations. We apply a forward model to predict brightness temperatures of snow covered sceneries over several winter seasons, using available data on snow cover, vegetation and lake ice cover to set the forward model input parameters. Comparison of model estimates with space-borne observations shows that the modeling accuracy improves in the majority of examined cases when lakes are accounted for, with respect to the case where lakes are not included in the simulation. Moreover, we present a method for applying the correction to the retrieval of Snow Water Equivalent (SWE) in lake-rich areas, using a numerical inversion method of the forward model. In a comparison to available independent validation data on SWE, also the retrieval accuracy is seen to improve when applying the influence of snow covered lakes in the emission model.     

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.     

Identification of snow cover regimes through spatial and temporal clustering of satellite microwave brightness temperatures

Areas of similar ecology are often delineated based on homogenous topography, temperature, and land cover. Once delineated, these zones become the basis for a wide variety of scientific research and management activities. For instance, in Canada, ecozones are commonly utilized ecological management units delineated using geographic, topographic, and climatic information aided by spring and summer vegetation conditions. Snow cover has an influence on local and regional hydrological conditions and climate, as well as on animal habitats. As such, we posit that inclusion of winter conditions, incorporating spatial- and temporal-variation in snow cover is an additional element for consideration when delineating areas with homogenous conditions. In our analysis we use satellite passive microwave brightness temperatures from 19 years of Special Sensor Microwave/Imager (SSM/I) measurements to produce a daily time-series on snow cover, and demonstrate how these data can be used to delineate areas of similar winter conditions. We use splines and curve fitting to generalize the dense time-series (of over 6900 days) to a set of metrics, and select three for use in cluster-based generalization of snow cover regimes: annual maximum difference between 37 and 19 GHz SSM/I measurements (with differences in magnitudes indicative of snow accumulation), variation of 37-19 GHz brightness temperatures (indicative of snow cover variability), and variation in the rate of brightness temperature change during the snow melt season (indicative of seasonal change). Our results indicate that these metrics produce spatial units that are unique, and not captured by conventional ecological management units, while also producing spatial units that cohere to those generated from summer conditions. Spatial units that are found to have spatial cohesion between summer and winter data sources are located in regions where the amount of snow tends to be low, and snow cover variability minimal. We propose that snow cover regimes may be used to augment typical vegetation-based ecological zonations or to provide insights on hydrology and animal habitat conditions. Inclusion of winter conditions is especially important when areal delineations are used to monitor impacts of climate change, and as a baseline for monitoring changes in snow cover amount, extent, and/or distribution.     

Estimation of snow water equivalent over first-year sea ice using AMSR-E and surface observations

A SWE retrieval algorithm developed in-situ using passive microwave surface based radiometer data is applied to the Advanced Microwave Scanning Radiometer for Earth Observation System (AMSR-E). Snow water equivalent is predicted from two pixels located in Canadian Arctic Shelf Exchange Study (CASES) overwintering study area in Franklin Bay, N.W.T., Canada. Results show that the satellite SWE predictions are statistically valid with measured in-situ snow thickness data in both smooth and rough ice environments where predicted values range from 15 to 25 mm. Stronger correlation between measured and predicted data is found over smooth ice with R² value of 0.75 and 0.73 for both pixels respectively. Furthermore, a qualitative study of sea ice roughness using both passive and active microwave satellite data shows that the two pixels are rougher than the surrounding areas, but the SWE predictions do not seem to be affected significantly.     

Factors affecting remotely sensed snow water equivalent uncertainty

State-of-the-art passive microwave remote sensing-based snow water equivalent (SWE) algorithms correct for factors believed to most significantly affect retrieved SWE bias and uncertainty. For example, a recently developed semi-empirical SWE retrieval algorithm accounts for systematic and random error caused by forest cover and snow morphology (crystal size — a function of location and time of year). However, we have found that climate and land surface complexities lead to significant systematic and random error uncertainties in remotely sensed SWE retrievals that are not included in current SWE estimation algorithms. Joint analysis of independent meteorological records, ground SWE measurements, remotely sensed SWE estimates, and land surface characteristics have provided a unique look at the error structure of these recently developed satellite SWE products. We considered satellite-derived SWE errors associated with the snow pack mass itself, the distance to significant open water bodies, liquid water in the snow pack and/or morphology change due to melt and refreeze, forest cover, snow class, and topographic factors such as large scale root mean square roughness and dominant aspect. Analysis of the nine-year Scanning Multichannel Microwave Radiometer (SMMR) SWE data set was undertaken for Canada where many in-situ measurements are available. It was found that for SMMR pixels with 5 or more ground stations available, the remote sensing product was generally unbiased with a seasonal maximum 20 mm average root mean square error for SWE values less than 100 mm. For snow packs above 100 mm, the SWE estimate bias was linearly related to the snow pack mass and the root mean square error increased to around 150 mm. Both the distance to open water and average monthly mean air temperature were found to significantly influence the retrieved SWE product uncertainty. Apart from maritime snow class, which had the greatest snow class affect on root mean square error and bias, all other factors showed little relation to observed uncertainties. Eliminating the drop-in-the-bucket averaged gridded remote sensing SWE data within 200 km of open water bodies, for monthly mean temperatures greater than − 2 °C, and for snow packs greater than 100 mm, has resulted in a remotely sensed SWE product that is useful for practical applications.     

Seasonal snow extent and snow mass in South America using SMMR and SSM/I passive microwave data (1979-2006)

Seasonal snow cover in South America was examined in this study using passive microwave satellite data from the Scanning Multichannel Microwave Radiometer (SMMR) on board the Nimbus-7 satellite and the Special Sensor Microwave Imagers (SSM/I) on board Defense Meteorological Satellite Program (DMSP) satellites. For the period from 1979-2006, both snow cover extent and snow water equivalent (snow mass) were investigated during the coldest months (May-September), primarily in the Patagonia area of Argentina and in the Andes of Chile, Argentina and Bolivia, where most of the seasonal snow is found. Since winter temperatures in this region are often above freezing, the coldest winter month was found to be the month having the most extensive snow cover and usually the month having the deepest snow cover as well. Sharp year-to-year differences were recorded using the passive microwave observations. The average snow cover extent for July, the month with the greatest average extent during the 28-year period of record, is 321,674 km². In July of 1984, the average monthly snow cover extent was 701,250 km² — the most extensive coverage observed between 1979 and 2006. However, in July of 1989, snow cover extent was only 120,000 km². The 28-year period of record shows a sinusoidal like pattern for both snow cover and snow mass, though neither trend is significant at the 95% level.     

Passive and active airborne microwave remote sensing of snow cover

The presence of snow cover affects the regional energy and water balance, thus having a significant impact on the global climate system. Temporal knowledge of the onset of snow melt and snow water equivalent (SWE) values are important variables in the prediction of flooding, as well as water resource applications such as reservoir management and agricultural activities. Microwave remote sensing techniques have been effective for monitoring snow pack parameters (snow extent, depth, water equivalent, wet/dry state). Coincident ground data, airborne polarimetric C-band (5.3 GHz) Synthetic Aperture Radar (SAR) and passive microwave radiometer data (19, 37 and 85 GHz) were collected on four dates (1 December 1997, 6 March 1998, 12 March 1998 and 9 March 1999) over two flight lines in Eastern Ontario, Canada. The multitemporal, multi-sensor data were analysed for changes in SAR polarimetric signatures and microwave brightness temperatures as a function of changing snow pack parameters. Results indicate that certain parameters such as linear polarizations and pedestal height are sensitive to changes in snow pack parameters, and respond differently to various snow conditions. SWE values derived from the passive microwave brightness temperatures compare well with ground measurements, with the exception of low snow volume and in the presence of significant ice layers.     

Intercomparison of microwave remote-sensing soil moisture data sets based on distributed eco-hydrological model simulation and in situ measurements over the North China Plain

Intercomparisons of microwave-based soil moisture products from active ASCAT (Advanced Scatterometer) and passive AMSR-E (Advanced Microwave Scanning Radiometer for the Earth Observing System) is conducted based on surface soil moisture (SSM) simulations from the eco-hydrological model, Vegetation Interface Processes (VIP), after it is carefully validated with in situ measurements over the North China Plain. Correlations with VIP SSM simulation are generally satisfactory with average values of 0.71 for ASCAT and 0.47 for AMSR-E during 2007–2009. ASCAT and AMSR-E present unbiased errors of 0.044 and 0.053 m³ m^?3 on average, with respect to model simulation. The empirical orthogonal functions (EOF) analysis results illustrate that AMSR-E provides more consistent SSM spatial structure with VIP than ASCAT; while ASCAT is more capable of capturing SSM temporal dynamics. This is supported by the facts that ASCAT has more consistent expansion coefficients corresponding to primary EOF mode with VIP (R = 0.825, p < 0.1). However, comparison based on SSM anomaly demonstrates that AMSR-E and ASCAT have similar skill in capturing SSM short-term variability. Temporal analysis of SSM anomaly time series shows that AMSR-E provides best performance in autumn, while ASCAT provides lower anomaly bias during highly-vegetated summer with vegetation optical depth of 0.61. Moreover, ASCAT retrieval accuracy is less influenced by vegetation cover, as it is in relatively better agreement with VIP simulation in forest than in other land-use types and exhibits smaller interannual fluctuation than AMSR-E. Identification of the error characteristics of these two microwave soil moisture data sets will be helpful for correctly interpreting the data products and also facilitate optimal specification of the error matrix in data assimilation at a regional scale.     

The impact of simulated soil temperatures on the estimation of snow depth over Siberia from SSM/I compared to a multi-model climatology

The snowpack is a key variable of the hydrological cycle. In recent years, numerous studies have demonstrated the importance of long-term monitoring of the Siberian snowpack on large spatial scales owing to evidence of increased river discharge, changes in snow fall amount and alterations with respect to the timing of ablation. This can currently only be accomplished using remote sensing methods. The main objective of this study is to take advantage of a new land surface forcing and simulation database in order to both improve and evaluate the snow depths retrieved using a dynamic snow depth retrieval algorithm. The dynamic algorithm attempts to account for the spatial and temporal internal properties of the snow cover. The passive microwave radiances used to derive snow depth were measured by the Special Sensor Microwave/ Imager (SSM/I) data between July 1987 and July 1995.The evaluation of remotely sensed algorithms is especially difficult over regions such as Siberia which are characterized by relatively sparse surface measurement networks. In addition, existing gridded climatological snow depth databases do not necessarily correspond to the same time period as the available satellite data. In order to evaluate the retrieval algorithm over Siberia for a recent multi-year period at a relatively large spatial scale, a land surface scheme reanalysis product from the Global Soil Wetness Project-Phase 2 (GSWP-2) is used in the current study. First, the high quality GSWP-2 input forcing data were used to drive a land surface scheme (LSS) in order to derive a climatological near-surface soil temperature. Four different snow depth retrieval methods are compared, two of which use the new soil temperature climatology as input. Second, a GSWP-2 snow water equivalent (SWE) climatology is computed from 12 state-of-the-art LSS over the same time period covered by the SSM/I data. This climatology was compared to the corresponding fields from the retrievals. This study reaffirmed the results of recent studies which showed that the inclusion of ancillary data into a satellite data-based snow retrieval algorithm, such as soil temperatures, can significantly improve the results. The current study also goes a step further and reveals the importance of including the monthly soil temperature variation into the retrieval, which improves results in terms of the spatial distribution of the snowpack. Finally, it is shown that further improved predictions of SWE are obtained when spatial and temporal variations in snow density are accounted for.     

Fraction images derived from NOAA AVHRR data for studying the deforestation in the Brazilian Amazon

Considerable variation in the performance of passive microwave global rainfall algorithms, both spatially and temporally,was revealed by the first WetNet Precipitation Intercomparison Project, PIP-1, with no one algorithm achieving the best results, in all locations, and all the time. In this paper a Composite Algorithm Procedure is described for the Special Sensor Microwave/Imager (SSM/I) algorithms submitted to PIP-1, that attempts through combining the best algorithm results from different regions of the globe to achieve better overall global rainfall estimates than are possible from any individual algorithm alone. The Composite Algorithm Procedure (CAP) involves the segmentation of the globe into homogeneous regions, the production of validation statistics for the various algorithm results in the different regions, and the identification of combinations of algorithms which perform best globally. The segmentations were based on aspects of the spatial and temporal variability of rainfall, or the microwave properties of the surfaces of the Earth. Initial results for the Composite Algorithm Procedure are presented for a sample month (October 1987): these confirm that improved global rainfall products can be produced in this way. Code detailing a selected Composite Algorithm based on the segmentation method of the microwave properties of the Earth has been supplied to the WetNet Support Group at the Marshall Space Flight Center, Huntsville, Alabama, for experimental, regular production of global rainfall data sets on a near real-time basis.     

Comparison of snow water equivalent retrieved from SSM/I passive microwave data using artificial neural network, projection pursuit and nonlinear regressions

The snow water equivalent (SWE) for the Red River basin of North Dakota and Minnesota was retrieved from data acquired by passive microwave SSM/I (Special Sensor Microwave Imager) sensors mounted on the US Defense Meteorological Satellite Program (DMSP) satellites, physiographic and atmospheric data by an artificial neural network called Modified Counter Propagation Network (MCPN), a Projection Pursuit Regression (PPR) and a nonlinear regression. The airborne gamma-ray measurements of SWE for 1989 and 1997 were used as observed SWE, and SSM/I data of 19 and 37 GHz frequencies, in both horizontal and vertical polarization, were used for the calibration (1989 data from DMSP-F8) and validation (1997 data from DMSP-F10 and F13 of both ascending and descending overpass times were combined) of the models. The SSM/I data were screened for the presence of wet snow, large water bodies like lakes and rivers, and depth-hoar. The MCPN model produced encouraging results in both calibration and validation stages (R² was about 0.9 for both calibration (C) and validation (V)), better than PPR (R² was 0.86 for C and 0.62 for V), which in turn was better than the multivariate nonlinear regression at the calibration stage (R² was 0.78 for C and 0.71 for V). MCPN is probably better than the linear and nonlinear regression counterparts because of its parallel computing structure resulted from neurons interconnected by a parallel network and its ability to learn and generalize information from complex relationships such as the SWE-SSM/I or other relationships encountered in geosciences.     

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

Simultaneous MSMR and SSM/I observations and analysis of sea‐ice characteristics over the Antarctic region

Polar sea ice has been monitored quasi‐continuously over the last 30 years using passive microwave radiometers onboard three satellites in polar orbit, namely Nimbus‐5, Nimbus‐7 and Defense Meteorological Satellite Program (DMSP) series. A short overlap between Scanning Multichannel Microwave Radiometer (SMMR) on Nimbus‐7 and Special Sensor Microwave Imager (SSM/I) onboard DMSP allowed inter‐calibration of the two sensors leading to a consistent series of long‐term sea‐ice measurements since 1978. With the launch of Multifrequency Scanning Microwave Radiometer (MSMR) onboard OCEANSAT‐1 in the polar sun‐synchronous orbit during 1999, India developed the capability to monitor the polar sea ice on a regular basis. The concurrent availability of SSM/I and MSMR over the last few years presents a valuable opportunity to attempt an inter‐comparison of MSMR with SSM/I measurements and derived sea‐ice parameters.

In this paper, we present an indirect validation of the brightness temperatures (T _b) observed by MSMR with near‐simultaneous measurements from SSM/I over the Antarctic and Southern Polar Ocean regions. Simultaneous MSMR and SSM/I data from two contrasting seasons—summer and winter—for the 1999–2000 period have been used. Analysis includes a comparison of T _b scatterograms to achieve confidence in the quantitative use of the T _b data to derive various geophysical parameters, e.g. sea‐ice concentration and extent. Additionally, the T _b images produced by the two sensors are compared to establish the capability of MSMR in reliable two‐dimensional portrayal of all the sea and continental ice features over the Antarctic Region. Based on a regression analysis between MSMR observed T _b in different frequency channels and polarizations, and SSM/I‐derived sea‐ice concentration (SIC) values, we have developed algorithms to estimate SIC over the Southern Polar Ocean from MSMR data. The MSMR algorithms allow estimation of SIC with better than 10% rms error. MSMR SIC images faithfully capture the observed distribution of sea ice in all the sectors of the Southern Ocean both during summer and winter periods. Using the MSMR‐derived SIC, we have also derived monthly sea‐ice extent (SIE) estimates for a period extending for about 20 months from the beginning of the launch of MSMR. These estimates show excellent agreement with values derived from SSM/I. These analyses bring out the very high level of compatibility in the measurements and derived sea‐ice parameters produced by the two sensors.     

A prototype MODIS-SSM/I snow-mapping algorithm

Data in the wavelength range 0.545-1.652 w m from the Moderate Resolution Imaging Spectroradiometer (MODIS), launched aboard the Earth Observing System (EOS) Terra in December 1999, are used to map daily global snow cover at 500 m resolution. However, during darkness, or when the satellite's view of the surface is obscured by cloud, snow cover cannot be mapped using MODIS data. We show that during these conditions, it is possible to supplement the MODIS product by mapping the snow cover using passive microwave data from the Special Sensor Microwave Imager (SSM/I), albeit with much poorer resolution. For a 7-day time period in March 1999, a prototype MODIS snow-cover product was compared with a prototype MODIS-SSM/I product for the same area in the mid-western USA. The combined MODIS-SSM/I product mapped 9% more snow cover than the MODIS-only product.     

Surface melting observations in Antarctica by microwave radiometers: Correcting 26-year time series from changes in acquisition hours

Surface melting duration and extent of the Antarctic coasts and ice-shelves is a climatic indicator related to the summer temperature and radiative budget. Surface melting is easily detectable by remote sensing using passive microwave observations. The preliminary goal of this study is to extend to 26 years an existing data set of surface melting [Torinesi, O., Fily, M., Genthon, C. (2003), Interannual variability and trend of the Antarctic summer melting period from 20 years of spaceborne microwave data, J. Climate, 16(7), pp. 1047-1060] by including the most recent years of observation. These data come from 4 microwave sensors (the Scanning Multichannel Microwave Radiometer (SMMR) and three Special Sensor Microwave Imager (SSM/I)) observing the surface at different hours of the day. Since surface melting varies throughout the day as the air temperature or the radiation, the interannual melting extent and duration time series are biased by sensor changes. Using all the sensors simultaneously available since 2002, we were able to model the diurnal variations of melting and use this hourly model to correct the long-term time series. This results in an unbiased 26-year long time series better suited for climate analysis. The cooling trend found by Torinesi et al. using uncorrected time series for the 1980-1999 period is confirmed but the decreasing rate is weaker after correction. Furthermore, extending the series up to summer 2004-2005 reveals recent changes: the last 2 summers have been particularly warmer over all the East Antarctica compared to the 10 previous years, thus ending the cold period of the 1990s. The trend over 1980-2005 is no longer toward cooling but complex climatic variations appear. The time series are available at http://www.lgge.obs.ujf-grenoble.fr/~picard/melting/.     

Monitoring sandstorms and desertification in northern China using SSM/I data and Getis statistics

A massive sandstorm enveloped most of northern China during Spring 2002. Monitoring the evolution of sandstorm and desertification has become one of the most serious problems for China's environment. Since 1989, one of the most advanced and operational passive microwave sensors is the Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave Imager (SSM/I) operated at seven channels (19, 37, 85?GHz with vertical and horizontal polarization and 22?GHz with vertical polarization only). In this paper, the sandstorm and desertification indexes, SDI and DI, are derived from the radiative transfer equation, and are employed with multi-channel measurements of the DMSP SSM/I for monitoring the sandstorm and desertification in northern China. Some SSM/I data in 1997 and 2001 are employed. The algorithm of the Getis statistics is developed to categorize the spatial correlation and its evolution during these days. It is demonstrated that the SSM/I indexes, SDI and DI, and its Getis statistics are well applicable for monitoring the sandstorm and desertification.     

Estimation of snow depth in the UK using the HUT snow emission model

Snow cover characteristics have significant effects on upwelling naturally emitted microwave radiation through processes of forward scattering. This study simulates numerically the electromagnetic responses from snow in the UK using the radiative transfer‐based semiempirical model developed at the Helsinki University of Technology (HUT), which takes into account the influence of soil surface, forest canopy and atmosphere on space‐borne observed brightness temperature by using empirical and semiempirical formulas. A sensitivity analysis of the HUT model was conducted to determine the most sensitive parameter affecting upwelling radiation from snow in the UK. The model‐based results were compared with observed Special Sensor Microwave Imager (SSM/I) brightness temperatures to better understand the SSM/I response to snow. The available ensemble of data required for input to the HUT model comprise surface physical temperature, ground level pressure and water vapour content, forest stem volume and land cover water fraction. Based on the sensitivity analyses, numerical parameters representing physical snow pack quantities (e.g. snow grain size, snow moisture and snow depth (SD)) were varied and the method of root mean square error (RMSE) minimization was used to invert the SD. The HUT model was applied to different days in 3 months (23–31 January, 1–5 and 26–27 February and 1–7 March 1995) of records of daily SD and SSM/I observations. The results show that the HUT model both underestimates and overestimates SD prediction. For the month of January 1995, the HUT model underestimated SD with a bias of ?0.59 cm, whereas for February and March 1995 the HUT model overestimated the SD with a bias of 1.89 cm and 1.64 cm, respectively. This study demonstrates that microwave remote sensing of snow can be used successfully in the UK, where most research on snow cover is conducted by using a visible and infrared radiometer. It is also evident from this work that application of algorithms to snow pack monitoring needs local calibration for effective and reasonable results.