Alain Royer

A simple retrieval method for land surface temperature and fraction of water surface determination from satellite microwave brightness temperatures in sub-arctic areas

A strong linear relationship is found between Special Sensor Microwave/Imager (SSM/I) microwave (19 and 37 GHz) surface emissivities at horizontal and vertical polarizations over snow- and ice-free land surfaces. This allows retrieving the land surface emissivity and temperature from satellite microwave brightness temperatures after atmospheric corrections. Over the Canadian sub-arctic continental area, we show that the main factor modifying the emissivity is the fraction of water surface (FWS) within a pixel. Accordingly, a map of the fraction of water surface across the Canadian landmass is derived, given a correspondence within 6% as compared to the 1 km 2 Canadian National Topographic Database of water-covered areas. The microwave-derived surface temperatures are compared to synchronous in situ air and ground surface temperatures and also with independent satellite IR measurements over areas without snow or ice. Root mean square differences range between 2j and 3.5j, with mean bias error of the order of 1‐3j. Better results are always obtained with the 37 GHz channel rather than with the 19 GHz channel. Over dense vegetation, the microwave-derived surface temperature is closer to the air temperature (at surface level) than to the ground temperature. The proposed simple retrieval algorithm, not sensitive to cloud cover, appears very useful for monitoring summer interannual or seasonal trends of the fraction of surface water, as well as the daily land surface temperature variation, which are very important parameters in environmental change analysis. D 2003 Elsevier Science Inc. All rights reserved.

Snow water equivalent retrieval in a Canadian boreal environment from microwave measurements using the HUT snow emission model

Snow water equivalent (SWE) is a critical parameter for climatological and hydrological studies over northern high-latitude areas. In this paper, we study the usability of the Helsinki University of Technology (HUT) snow emission model for the estimation of SWE in a Canadian boreal forest environment. The experimental data (airborne passive microwave and ground-based data) were acquired during the Boreal Ecosystem-Atmosphere Study winter field campaign held in February 1994 in Central Canada. Using the experimental dataset, surface brightness temperatures at 18 and 37 GHz (vertical polarization) were simulated with the HUT snow emission model and compared to those acquired by the airborne sensors. The results showed an important underestimation at 37 GHz (-27 K) and an overestimation at 18 GHz (10 K). In this paper, we demonstrate that the errors in the model simulations are due mainly to the extinction coefficient modeling, which is a function of snow grain size. Therefore, we propose a new semiempirical function for the extinction coefficient, based on an empirical correction to the Rayleigh scattering expression. Results presented in this paper show that the proposed function improves the HUT model accuracy to predict brightness temperature in the experimental context considered, with a mean error of /spl plusmn/5 K and /spl plusmn/9 K, respectively, at 18 and 37 GHz, and a negligible bias (less than 4 K) in both cases. These errors are comparable in magnitude to the accuracy of the radiometers used during the airborne flights. SWE was retrieved using the modified HUT snow emission model based on an iterative inversion technique. SWE was estimated with a mean error of /spl plusmn/10 mm and a negligible bias. Only a rough knowledge of mean snow grain size /spl phi/~ was required in the inversion procedure. The effects of possible errors on mean snow grain size /spl phi/~ are presented and discussed.

Analysis of Temperature Emissivity Separation (TES) algorithm applicability and sensitivity

The purpose of this paper is to assess the spectral Temperature Emissivity Separation algorithm (TES) proposed by Gillespie et al. (1998) as a simple method to retrieve surface emissivity from ground-based measurements. First, we validate different empirical relationships for the Minimum Maximum Difference module, on which the TES is based, with a large dataset (about 500 surfaces from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) spectral library including man-made materials) for multiband data in the long wave infrared (LWIR: 7.5–14 µm), and hyperspectral data in the middle wave infrared (MWIR: 3.4–5.2 µm) and LWIR. We show the applicability of TES for hyperspectral data using a specific empirical relationship; this is confirmed by experimental measurements. For multiband data, we improve the TES for high contrast emissivity surfaces by integrating broadband 8–14 µm measurements in the iterative algorithm. We also found that metals do not confirm these empirical relationships. TES accuracy, extensively assessed by simulations, remains for multiband simulations (respectively for hyperspectral) within about 0.03 (0.02) for emissivity and within about 1.2 K (0.3 K) for temperature. However, surfaces with low maximum emissivity give higher errors. Except for these particular surfaces, the TES approach, applied on measurements from a portable multiband thermal radiometer, appears as the most efficient and accurate method for emissivity determination in the field without any a priori assumption on the surface nature.

Characterization of Atmospheric Aerosols across Canada from a Ground-based Sunphotometer Network: AEROCAN

Abstract Ground‐based sunphotometry measurements acquired under clear sky conditions can be used to investigate atmospheric aerosol optical properties. Such measurements are not only important in their own right as a technique for monitoring generic aerosol dynamics, but also represent a direct means of evaluating the contribution of aerosol induced radiative forcing in the modelling of climate change. In this paper we analyze derived aerosol optical properties using datasets from the Canadian AEROCAN (AERosol CANada) sunphotometer network. The AEROCAN network currently includes eight sunphotometers distributed across Canada at sites chosen in order to obtain a diverse sampling of continental, maritime and arctic aerosols. Some of these sites have been operational since 1993 as part of the Boreal Ecosystem‐Atmosphere Study (BOREAS). These instruments permit standard and automatic multi‐wavelength measurements of solar extinction radiance centred on the solar disk as well as sky radiance scans off the solar disk. These data yields aerosol optical depth, the Ångström exponent, aerosol particle volume size distribution, refractive index, column‐averaged single scattering albedo, and precipitable water vapour content. Spatial and temporal trends of these parameters as well as observed inter‐correlations are discussed. The results demonstrate the utility and significance of these types of measurements and illustrate the potential applications of networked sunphotometry data.

Coupling the snow thermodynamic model SNOWPACK with the microwave emission model of layered snowpacks for subarctic and arctic snow water equivalent retrievals

[1] Satellite-passive microwave remote sensing has been extensively used to estimate snow water equivalent (SWE) in northern regions. Although passive microwave sensors operate independent of solar illumination and the lower frequencies are independent of atmospheric conditions, the coarse spatial resolution introduces uncertainties to SWE retrievals due to the surface heterogeneity within individual pixels. In this article, we investigate the coupling of a thermodynamic multilayered snow model with a passive microwave emission model. Results show that the snow model itself provides poor SWE simulations when compared to field measurements from two major field campaigns. Coupling the snow and microwave emission models with successive iterations to correct the influence of snow grain size and density significantly improves SWE simulations. This method was further validated using an additional independent data set, which also showed significant improvement using the two-step iteration method compared to standalone simulations with the snow model.

Snow Microwave Emission Modeling of Ice Lenses Within a Snowpack Using the Microwave Emission Model for Layered Snowpacks

Ice lens formation, which follows rain on snow events or melt-refreeze cycles in winter and spring, is likely to become more frequent as a result of increasing mean winter temperatures at high latitudes. These ice lenses significantly affect the microwave scattering and emission properties, and hence snow brightness temperatures that are widely used to monitor snow cover properties from space. To understand and interpret the spaceborne microwave signal, the modeling of these phenomena needs improvement. This paper shows the effects and sensitivity of ice lenses on simulated brightness temperatures using the microwave emission model of layered snowpacks coupled to a soil emission model at 19 and 37 GHz in both horizontal and vertical polarizations. Results when considering pure ice lenses show an improvement of 20.5 K of the root mean square error between the simulated and measured brightness temperature (Tb) using several in situ data sets acquired during field campaigns across Canada. The modeled Tbs are found to be highly sensitive to the vertical location of ice lenses within the snowpack.

Improved Corrections of Forest Effects on Passive Microwave Satellite Remote Sensing of Snow Over Boreal and Subarctic Regions

Microwave radiometry has been extensively used in order to estimate snow water equivalent in northern regions. However, for boreal and taiga environments, the presence of forest causes important uncertainties in the estimates. Variations in snow cover and vegetation in northeastern Canada (north of the Québec province) were characterized in a transect from 50°N to 60 °N during the International Polar Year field campaign of February 2008. Forest properties show a strong latitudinal gradient in fraction and stem volume. A large database (>; 2000 points with a stem volume ranging between 0 and 700 m3 ·ha-1) showed that brightness temperatures (Tb) decrease as forest cover fraction decreases until a cover fraction of about 25% is reached. Furthermore, Tb values saturate at high stem volume, particularly at 37 GHz. We defined new relationships for the forest transmissivity as a function of stem volume and depending on the frequency/polarization. The proposed relationships give asymptotic transmissivity saturation levels of 0.51, 0.55, 0.53, and 0.53 for 19 GHz [vertical (V) polarization], 19 GHz [horizontal (H) polarization], 37 GHz (V polarization), and 37 GHz (H polarization), respectively. These relationships were used to estimate snow Tb from the Advanced Microwave Scanning Radiometer-Earth Observing System brightness temperatures at 18.7 and 36.5 GHz, and results show an estimated snow brightness temperature well correlated to the airborne snow brightness temperatures over vegetation-free areas.

Brightness Temperature Simulations of the Canadian Seasonal Snowpack Driven by Measurements of the Snow Specific Surface Area

Snow grain size is the snowpack parameter that most affects the microwave snow emission. The specific surface area (SSA) of snow is a metric that allows rapid and reproducible field measurements and that well represents the grain size. However, this metric cannot be used directly in microwave snow emission models (MSEMs). The aim of this paper is to evaluate the suitability and the adaptations required for using the SSA in two MSEMs, i.e., the Dense Media Radiative Theory-Multilayer model (DMRT-ML) and the Helsinki University of Technology model (HUT n-layer), based on in situ radiometric measurements. Measurements of the SSA, using snow reflectance in the short-wave infrared, were taken at 20 snowpits in various environments (e.g., grass, tundra, and dry fen). The results show that both models required a scaling factor for the SSA values to minimize the root-mean-square error between the measured and simulated brightness temperatures. For DMRT-ML, the need for a scaling factor is likely due to the oversimplified representation of snow as spheres of ice with a uniform radius. We hypothesize that the need for a scaling factor is related to the grain size distribution of snow and the stickiness between grains. For HUT n-layer, using the SSA underestimates the attenuation by snow, particularly for snowpacks with a significant amount of depth hoar. This paper provides a reliable description of the grain size for DMRT-ML, which is of particular interest for the assimilation of satellite passive microwave data in snow models.

Extraction of bimodal aerosol-size distribution radii from spectral and angular slope (Angstrom) coefficients

A methodology is defined whereby the effective aerosol bimodal radii of the fine (accumulation) mode and the coarse (supermicron) mode can be directly and simply extracted from slope (Angstrom) coefficients derived from ground-based solar spectral transmission data and aureole scattering data.

Evaluation of Spaceborne L-Band Radiometer Measurements for Terrestrial Freeze/Thaw Retrievals in Canada

The landscape freeze/thaw (F/T) state has an important impact on the surface energy balance, carbon fluxes, and hydrologic processes; the timing of spring melt is linked to active layer dynamics in permafrost areas. L-band (1.4 GHz) microwave emission could allow the monitoring of surface state dynamics due to its sensitivity to the pronounced permittivity difference between frozen and thawed soil. The aim of this paper is to evaluate the performance of both Aquarius and soil moisture and ocean salinity (SMOS) L-band passive microwave measurements using a polarization ratio (PR)-based algorithm for landscape F/T monitoring. Weekly L-band satellite observations are compared with a large set of reference data at 48 sites across Canada spanning three environments: 1) tundra; 2) boreal forest; and 3) prairies. The reference data include in situ measurements of soil temperature (Tsoil) and air temperature (Tair), and moderate resolution imaging spectroradiometer (MODIS) land surface temperature (LST) and snow cover area (SCA) products. Results show generally good agreement between L-band F/T detection and the surface state estimated from four reference datasets. The best apparent accuracies for all seasons are obtained using Tair as the reference. Aquarius radiometer 2 (incidence angle of 39.6°) data give the best accuracies (90.8%), while for SMOS, using the Aquarius temporal characteristics, the best results (87.8% of accuracy) are obtained at higher incidence angles (55°-60°). The F/T algorithm identifies both freeze onset and end with a delay of about 1 week in tundra and 2 weeks in forest and prairies, when compared to Tair. The analysis shows a stronger F/T signal at tundra sites due to the typically clean transitions between consistently frozen and thawed conditions (and vice versa) and the absence of surface vegetation. Results in the prairies were poorer because of the influence of vegetation growth in summer (which decreases the PR) and the high frequency of ephemeral thaw events during winter. Freeze onset and end maps created from the same algorithm applied to SMOS and Aquarius measurements characterize similar F/T patterns over Canada. This study shows the potential of using L-band spaceborne observations for F/T monitoring, but underlines some limitations due to ice crusts in the snowpack, liquid water content in snow cover during the spring freeze to thaw transition, and vegetation growth.