Automatic classification-based generation of thermal infrared land surface emissivity maps using AATSR data over Europe

This is the author’s version of a work that was accepted for publication in Remote Sensing of Environment. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Remote Sensing of Environment, 124, 321-333.DOI :10.1016/j.rse.2012.05.024.The remote sensing measurement of land surface temperature from satellites provides a monitoring of this magnitude on a continuous and regular basis, which is a critical factor in many research fields such as weather forecasting, detection of forest fires or climate change studies, for instance. The main problem of measuring temperature from space is the need to correct for the effects of the atmosphere and the surface emissivity. In this work an automatic procedure based on the Vegetation Cover Method, combined with the GLOBCOVER land surface type classification, is proposed. The algorithm combines this land cover classification with remote sensing information on the vegetation cover fraction to obtain land surface emissivity maps for AATSR split-window bands. The emissivity estimates have been compared with ground measurements in two validation cases in the area of rice fields of Valencia, Spain, and they have also been compared to the classification-based emissivity product provided by MODIS (MOD11_L2). The results show that the error in emissivity of the proposed methodology is of the order of ±0.01 for most of the land surface classes considered, which will contribute to improve the operational land surface temperature measurements provided by the AATSR instrument. © 2012 Elsevier Inc. All rights reserved.This work was funded by the Generalitat Valenciana (project PRO-METEO/2009/086, and contract of Eduardo Caselles) and the Spanish Ministerio de Ciencia e Innovacion (projects CGL2007-64666/CLI, CGL2010-17577/CLI and CGL2007-29819-E, co-financed with FEDER funds). AATSR data were provided by European Space Agency (ESA) under Cat-1 project 3466. We also thank ESA and the ESA GLOBCOVER Project, led by MEDIAS-France, for the GLOBCOVER classification data. The comments and suggestions of three anonymous reviewers that improved the paper are also acknowledged.Caselles, E.; Valor, E.; Abad Cerdá, FJ.; Caselles, V. (2012). Automatic classification-based generation of thermal infrared land surface emissivity maps using AATSR data over Europe. Remote Sensing of Environment. 124:321-333. https://doi.org/10.1016/j.rse.2012.05.024S32133312

Harmonization of remote sensing land surface products : correction of clear-sky bias and characterization of directional effects

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Deteção Remota), Universidade de Lisboa, Faculdade de Ciências, 2018Land surface temperature (LST) is the mean radiative skin temperature of an area of land resulting from the mean energy balance at the surface. LST is an important climatological variable and a diagnostic parameter of land surface conditions, since it is the primary variable determining the upward thermal radiation and one of the main controllers of sensible and latent heat fluxes between the surface and the atmosphere. The reliable and long-term estimation of LST is therefore highly relevant for a wide range of applications, including, amongst others: (i) land surface model validation and monitoring; (ii) data assimilation; (iii) hydrological applications; and (iv) climate monitoring. Remote sensing constitutes the most effective method to observe LST over large areas and on a regular basis. Satellite LST products generally rely on measurements in the thermal infrared (IR) atmospheric window, i.e., within the 8-13 micrometer range. Beside the relatively weak atmospheric attenuation under clear sky conditions, this band includes the peak of the Earth’s spectral radiance, considering surface temperature of the order of 300K (leading to maximum emission at approximately 9.6 micrometer, according to Wien’s Displacement Law). The estimation of LST from remote sensing instruments operating in the IR is being routinely performed for nearly 3 decades. Nevertheless, there is still a long list of open issues, some of them to be addressed in this PhD thesis. First, the viewing position of the different remote sensing platforms may lead to variability of the retrieved surface temperature that depends on the surface heterogeneity of the pixel – dominant land cover, orography. This effect introduces significant discrepancies among LST estimations from different sensors, overlapping in space and time, that are not related to uncertainties in the methodologies or input data used. Furthermore, these directional effects deviate LST products from an ideally defined LST, which should correspond to the ensemble directional radiometric temperature of all surface elements within the FOV. In this thesis, a geometric model is presented that allows the upscaling of in situ measurements to the any viewing configuration. This model allowed generating a synthetic database of directional LST that was used consistently to evaluate different parametric models of directional LST. Ultimately, a methodology is proposed that allows the operational use of such parametric models to correct angular effects on the retrieved LST. Second, the use of infrared data limits the retrieval of LST to clear sky conditions, since clouds “close” the atmospheric window. This effect introduces a clear-sky bias in IR LST datasets that is difficult to quantify since it varies in space and time. In addition, the cloud clearing requirement severely limits the space-time sampling of IR measurements. Passive microwave (MW) measurements are much less affected by clouds than IR observations. LST estimates can in principle be derived from MW measurements, regardless of the cloud conditions. However, retrieving LST from MW and matching those estimations with IR-derived values is challenging and there have been only a few attempts so far. In this thesis, a methodology is presented to retrieve LST from passive MW observations. The MW LST dataset is examined comprehensively against in situ measurements and multiple IR LST products. Finally, the MW LST data is used to assess the spatial-temporal patterns of the clear-sky bias at global scale.Fundação para a Ciência e a Tecnologia, SFRH/BD/9646

A Comparison of National Water Model Retrospective Analysis Snow Outputs at Snow Telemetry Sites Across the Western United States

This study compares the US National Water Model (NWM) reanalysis snow outputs to observed snow water equivalent (SWE) and snow-covered area fraction (SCAF) at snow telemetry (SNOTEL) sites across the Western United States. SWE was obtained from SNOTEL sites, while SCAF was obtained from moderate resolution imaging spectroradiometer (MODIS) observations at a nominal 500 m grid scale. Retrospective NWM results were at a 1000 m grid scale. We compared results for SNOTEL sites to gridded NWM and MODIS outputs for the grid cells encompassing each SNOTEL site. Differences between modelled and observed SWE were attributed to both model errors, as well as errors in inputs, notably precipitation and temperature. The NWM generally under-predicted SWE, partly due to precipitation input differences. There was also a slight general bias for model input temperature to be cooler than observed, counter to the direction expected to lead to under-modelling of SWE. There was also under-modelling of SWE for a subset of sites where precipitation inputs were good. Furthermore, the NWM generally tends to melt snow early. There was considerable variability between modelled and observed SCAF as well as the binary comparison of snow cover presence that hampered useful interpretation of SCAF comparisons. This is in part due to the shortcomings associated with both model SCAF parameterization and MODIS observations, particularly in vegetated regions. However, when SCAF was aggregated across all sites and years, modelled SCAF tended to be more than observed using MODIS. These differences are regional with generally better SWE and SCAF results in the Central Basin and Range and differences tending to become larger the further away regions are from this region. These findings identify areas where predictions from the NWM involving snow may be better or worse, and suggest opportunities for research directed towards model improvements

Generating global products of LAI and FPAR from SNPP-VIIRS data: theoretical background and implementation

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation have been successfully generated from the Moderate Resolution Imaging Spectroradiometer (MODIS) data since early 2000. As the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard, the Suomi National Polar-orbiting Partnership (SNPP) has inherited the scientific role of MODIS, and the development of a continuous, consistent, and well-characterized VIIRS LAI/FPAR data set is critical to continue the MODIS time series. In this paper, we build the radiative transfer-based VIIRS-specific lookup tables by achieving minimal difference with the MODIS data set and maximal spatial coverage of retrievals from the main algorithm. The theory of spectral invariants provides the configurable physical parameters, i.e., single scattering albedos (SSAs) that are optimized for VIIRS-specific characteristics. The effort finds a set of smaller red-band SSA and larger near-infraredband SSA for VIIRS compared with the MODIS heritage. The VIIRS LAI/FPAR is evaluated through comparisons with one year of MODIS product in terms of both spatial and temporal patterns. Further validation efforts are still necessary to ensure the product quality. Current results, however, imbue confidence in the VIIRS data set and suggest that the efforts described here meet the goal of achieving the operationally consistent multisensor LAI/FPAR data sets. Moreover, the strategies of parametric adjustment and LAI/FPAR evaluation applied to SNPP-VIIRS can also be employed to the subsequent Joint Polar Satellite System VIIRS or other instruments.Accepted manuscrip

Remote sensing phenology at European northern latitudes - From ground spectral towers to satellites

Plant phenology exerts major influences on carbon, water, and energy exchanges between atmosphere and ecosystems, provides feedbacks to climate, and affects ecosystem functioning and services. Great efforts have been spent in studying plant phenology over the past decades, but there are still large uncertainties and disputations in phenology estimation, trends, and its climate sensitivities. This thesis aims to reduce these uncertainties through analyzing ground spectral sampling, developing methods for in situ light sensor calibration, and exploring a new spectral index for reliable retrieval of remote sensing phenology and climate sensitivity estimation at European northern latitudes. The ground spectral towers use light sensors of either nadir or off-nadir viewing to measure reflected radiation, yet how plants in the sensor view contribute differently to the measured signals, and necessary in situ calibrations are often overlooked, leading to great uncertainties in ground spectral sampling of vegetation. It was found that the ground sampling points in the sensor view follow a Cauchy distribution, which is further modulated by the sensor directional response function. We proposed in situ light sensor calibration methods and showed that the user in situ calibration is more reliable than manufacturer’s lab calibration when our proposed calibration procedures are followed. By taking the full advantages of more reliable and standardized reflectance, we proposed a plant phenology vegetation index (PPI), which is derived from a radiative transfer equation and uses red and near infrared reflectance. PPI shows good linearity with canopy green leaf area index, and is correlated with gross primary productivity, better than other vegetation indices in our test. With suppressed snow influences, PPI shows great potentials for retrieving phenology over coniferous-dominated boreal forests. PPI was used to retrieve plant phenology from MODIS nadir BRDF-adjusted reflectance at European northern latitudes for the period 2000-2014. We estimated the trend of start of growing season (SOS), end of growing season (EOS), length of growing season (LOS), and the PPI integral for the time span, and found significant changes in most part of the region, with an average rate of -0.39 days·year-1 in SOS, 0.48 days·year-1 in EOS, 0.87 days·year-1 in LOS, and 0.79%·year-1 in the PPI integral over the past 15 years. We found that the plant phenology was significantly affected by climate in most part of the region, with an average sensitivity to temperature: SOS at -3.43 days·°C-1, EOS at 1.27 days·°C-1, LOS at 3.16 days·°C-1, and PPI integral at 2.29 %·°C-1, and to precipitation: SOS at 0.28 days∙cm-1, EOS at 0.05 days∙cm-1, LOS at 0.04 days∙cm-1, and PPI integral at -0.07%∙cm-1. These phenology variations were significantly related to decadal variations of atmospheric circulations, including the North Atlantic Oscillation and the Arctic Oscillation. The methods developed in this thesis can help to improve the reliability of long-term field spectral measurements and to reduce uncertainties in remote sensing phenology retrieval and climate sensitivity estimation

Downscaling Aerosol Optical Thickness from Satellite Observations: Physics and Machine Learning Approaches

In recent years, the satellite observation of aerosol properties has been greatly improved. As a result, the derivation of Aerosol Optical Thickness (AOT), one of the most popular atmospheric parameters used in air pollution monitoring, over ocean and continents from satellite observations shows comparable quality to ground-based measurements. Satellite AOT products is often applied for monitoring at global scale because of its coarse spatial resolution. However, monitoring at local scale such as over cities requires more detailed AOT information. The increase spatial resolution to suitable level has potential for applications of air pollution monitoring at global-to-local scale, detecting emission sources, deciding pollution management strategies, localizing aerosol estimation, etc. In this thesis, we investigated, proposed, implemented and validated algorithms to derive AOT maps with spatial resolution increased up to 1×1 km2 from MODerate resolution Imaging Spectrometer (MODIS) observations provided by National Aeronautics and Space Administration (NASA), while MODIS standard aerosol products provide maps at 10×10 km2 of spatial resolution. The solutions are considered on two perspectives: dynamical downscaling by improving the algorithm for remote sensing of tropospheric aerosol from MODIS and statistical downscaling using Support Vector Regression

Estimating snow cover from publicly available images

In this paper we study the problem of estimating snow cover in mountainous regions, that is, the spatial extent of the earth surface covered by snow. We argue that publicly available visual content, in the form of user generated photographs and image feeds from outdoor webcams, can both be leveraged as additional measurement sources, complementing existing ground, satellite and airborne sensor data. To this end, we describe two content acquisition and processing pipelines that are tailored to such sources, addressing the specific challenges posed by each of them, e.g., identifying the mountain peaks, filtering out images taken in bad weather conditions, handling varying illumination conditions. The final outcome is summarized in a snow cover index, which indicates for a specific mountain and day of the year, the fraction of visible area covered by snow, possibly at different elevations. We created a manually labelled dataset to assess the accuracy of the image snow covered area estimation, achieving 90.0% precision at 91.1% recall. In addition, we show that seasonal trends related to air temperature are captured by the snow cover index.Comment: submitted to IEEE Transactions on Multimedi

Estimating Regional Snow Line Elevation Using Public Webcam Images

Snow cover is of high relevance for the Earth’s climate system, and its variability plays a key role in alpine hydrology, ecology, and socioeconomic systems. Measurements obtained by optical satellite remote sensing are an essential source for quantifying snow cover variability from a local to global scale. However, the temporal resolution of such measurements is often affected by persistent cloud coverage, limiting the application of high resolution snow cover mapping. In this study, we derive the regional snow line elevation in an alpine catchment area using public webcams. We compare our results to the snow line information derived from the Moderate-Resolution Imaging Spectroradiometer (MODIS) and Sentinel-2 snow cover products and find our results to be in good agreement therewith. Between October 2017 and the end of June 2018, snow lines derived from webcams lie on average 55.8 m below and 33.7 m above MODIS snow lines using a normalized-difference snow index (NDSI) of 0.4 and 0.1, respectively, and are on average 53.1 m below snow lines derived from Sentinel-2. We further analyze the superior temporal resolution of webcam-based snow cover information and demonstrate its effectiveness in filling temporal gaps in satellite-based measurements caused by cloud cover. Our findings show the ability of webcam-based snow line elevation retrieval to complement and improve satellite-based measurements

Cloud property datasets retrieved from AVHRR, MODIS, AATSR and MERIS in the framework of the Cloud_cci project

New cloud property datasets based on measurements from the passive imaging satellite sensors AVHRR, MODIS, ATSR2, AATSR and MERIS are presented. Two retrieval systems were developed that include components for cloud detection and cloud typing followed by cloud property retrievals based on the optimal estimation (OE) technique. The OE-based retrievals are applied to simultaneously retrieve cloud-top pressure, cloud particle effective radius and cloud optical thickness using measurements at visible, near-infrared and thermal infrared wavelengths, which ensures spectral consistency. The retrieved cloud properties are further processed to derive cloud-top height, cloud-top temperature, cloud liquid water path, cloud ice water path and spectral cloud albedo. The Cloud_cci products are pixel-based retrievals, daily composites of those on a global equal-angle latitude–longitude grid, and monthly cloud properties such as averages, standard deviations and histograms, also on a global grid. All products include rigorous propagation of the retrieval and sampling uncertainties. Grouping the orbital properties of the sensor families, six datasets have been defined, which are named AVHRR-AM, AVHRR-PM, MODIS-Terra, MODIS-Aqua, ATSR2-AATSR and MERIS+AATSR, each comprising a specific subset of all available sensors. The individual characteristics of the datasets are presented together with a summary of the retrieval systems and measurement records on which the dataset generation were based. Example validation results are given, based on comparisons to well- established reference observations, which demonstrate the good quality of the data. In particular the ensured spectral consistency and the rigorous uncertainty propagation through all processing levels can be considered as new features of the Cloud_cci datasets compared to existing datasets. In addition, the consistency among the individual datasets allows for a potential combination of them as well as facilitates studies on the impact of temporal sampling and spatial resolution on cloud climatologies