Validation of Collection 6 MODIS land surface temperature product using in situ measurements

Land surface temperature (LST) is an important physical quantity at the land-atmosphere interface. Since 2016 the Collection 6 (C6) MODIS LST product is publicly available, which includes three refinements over bare soil surfaces compared to the Collection 5 (C5) MODIS LST product. To encourage the use of the C6 MODIS LST product in a wide range of applications, it is necessary to evaluate the accuracy of the C6 MODIS LST product. In this study, we validated the C6 MODIS LST product using temperature-based method over various land cover types, including grasslands, croplands, cropland/natural vegetation mosaic, open shrublands, woody savannas, and barren/sparsely vegetated. In situ measurements were collected from various sites under different atmospheric and surface conditions, including seven SURFRAD sites (BND, TBL, DRA, FPK, GCM, PSU, and SXF) in the United States, three KIT sites (EVO, KAL, and GBB) in Portugal and Namibia, and three HiWATER sites (GBZ, HZZ, and HMZ) in China. The spatial representativeness of the in situ measurements at each site was separately evaluated during daytime and nighttime using all available clear-sky ASTER LST products at 90m spatial resolution. Only six sites during daytime are selected as sufficiently homogeneous sites despite the usually high spatial thermal heterogeneity, whereas during nighttime most sites can be considered to be thermally homogeneous and have similar LST and air temperature. The C6 MODIS LST product was validated using in situ measurements from the selected homogeneous sites during daytime and nighttime: except for the GBB site, large RMSE values (> 2 K) were obtained during daytime. However, if only satellite LST with a high spatial thermal homogeneity on the MODIS pixel scale are used for LST validation, the best daytime accuracy (RMSE<1.3 K) for the C6 MODIS LST product is achieved over the BND and DRA sites. Except for the DRA site, the RMSE values during nighttime are<2 K at the selected homogeneous sites. Furthermore, the accuracy of the C6 MODIS LST product was compared with that of the C5 MODIS LST product during nighttime at the selected homogeneous sites. Except for the GBB site, there are only small differences (< 0.4 K) between the RMSE values for the C5 and C6 MODIS LST products

Land Surface Temperature Product Validation Best Practice Protocol Version 1.0 - October, 2017

The Global Climate Observing System (GCOS) has specified the need to systematically generate andvalidate Land Surface Temperature (LST) products. This document provides recommendations on goodpractices for the validation of LST products. Internationally accepted definitions of LST, emissivity andassociated quantities are provided to ensure the compatibility across products and reference data sets. Asurvey of current validation capabilities indicates that progress is being made in terms of up-scaling and insitu measurement methods, but there is insufficient standardization with respect to performing andreporting statistically robust comparisons.Four LST validation approaches are identified: (1) Ground-based validation, which involvescomparisons with LST obtained from ground-based radiance measurements; (2) Scene-based intercomparisonof current satellite LST products with a heritage LST products; (3) Radiance-based validation,which is based on radiative transfer calculations for known atmospheric profiles and land surface emissivity;(4) Time series comparisons, which are particularly useful for detecting problems that can occur during aninstrument's life, e.g. calibration drift or unrealistic outliers due to undetected clouds. Finally, the need foran open access facility for performing LST product validation as well as accessing reference LST datasets isidentified

Detection of vegetation drying signals using diurnal variation of land surface temperature: Application to the 2018 East Asia heatwave

Satellite-based vegetation monitoring provides important insights regarding spatiotemporal variations in vegetation growth from a regional to continental scale. Most current vegetation monitoring methodologies rely on spectral vegetation indices (VIs) observed by polar-orbiting satellites, which provide one or a few observations per day. This study proposes a new methodology based on diurnal changes in land surface temperatures (LSTs) using Japan's geostationary satellite, Himawari-8/Advanced Himawari Imager (AHI). AHI thermal infrared observation provides LSTs at 10-min frequencies and ∼ 2 km spatial resolution. The DTC parameters that summarize the diurnal cycle waveform were obtained by fitting a diurnal temperature cycle (DTC) model to the time-series LST information for each day. To clarify the applicability of DTC parameters in detecting vegetation drying under humid climates, DTC parameters from in situ LSTs observed at vegetation sites, as well as those from Himawari-8 LSTs, were evaluated for East Asia. Utilizing the record-breaking heat wave that occurred in East Asia in 2018 as a case study, the anomalies of DTC parameters from the Himawari-8 LSTs were compared with the drying signals indicated by VIs, latent heat fluxes (LE), and surface soil moisture (SM). The results of site-based and satellite-based analyses revealed that DTR (diurnal temperature range) correlates with the evaporative fraction (EF) and SM, whereas Tmax (daily maximum LST) correlates with LE and VIs. Regarding other temperature-related parameters, T0 (LST around sunrise), Ta (temperature rise during daytime), and δT (temperature fall during nighttime) are unstable in quantification by DTC model. Moreover, time-related parameters, such as tm (time reaching Tmax), are more sensitive to topographic slope and geometric conditions than surface thermal properties at humid sites in East Asia, although they correlate with EF and SM at a semi-arid site in Australia. Additionally, the spatial distribution of the DTR anomaly during the 2018 heat wave corresponds with the drying signals indicated as negative SM anomalies. Regions with large positive anomalies in Tmax and DTR correspond to area with visible damage to vegetation, as indicated by negative VI anomalies. Hence, combined Tmax and DTR potentially detects vegetation drying indetectable by VIs, thereby providing earlier and more detailed vegetation monitoring in both humid and semi-arid climates

Empirical approach to satellite snow detection

Lumipeitteellä on huomattava vaikutus säähän, ilmastoon, luontoon ja yhteiskuntaan. Pelkästään sääasemilla tehtävät lumihavainnot (lumen syvyys ja maanpinnan laatu) eivät anna kattavaa kuvaa lumen peittävyydestä tai muista lumipeitteen ominaisuuksista. Sääasemien tuottamia havaintoja voidaan täydentää satelliiteista tehtävillä havainnoilla. Geostationaariset sääsatelliitit tuottavat havaintoja tihein välein, mutta havaintoresoluutio on heikko monilla alueilla, joilla esiintyy kausittaista lunta. Polaariradoilla sääsatelliittien havaintoresoluutio on napa-alueiden läheisyydessä huomattavasti parempi, mutta silloinkaan satelliitit eivät tuota jatkuvaa havaintopeittoa. Tiheimmän havaintoresoluution tuottavat sääsatelliittiradiometrit, jotka toimivat optisilla aallonpituuksilla (näkyvä valo ja infrapuna). Lumipeitteen kaukokartoitusta satelliiteista vaikeuttavat lumipeitteen oman vaihtelun lisäksi pinnan ominaisuuksien vaihtelu (kasvillisuus, vesistöt, topografia) ja valaistusolojen vaihtelu. Epävarma ja osittain puutteellinen tieto pinnan ja kasvipeitteen ominaisuuksista vaikeuttaa luotettavan automaattisen analyyttisen lumentunnistusmenetelmän kehittämistä ja siksi empiirinen lähestymistapa saattaa olla toimivin vaihtoehto automaattista lumentunnistusmenetelmää kehitettäessä. Tässä työssä esitellään kaksi EUMETSATin osittain rahoittamassa H SAFissa kehitettyä lumituotetta ja niissä käytetyt empiiristä lähestymistapaa soveltaen kehitetyt algoritmit. Geostationaarinen MSG/SEVIRI H31 lumituote on saatavilla vuodesta 2008 alkaen ja polaarituote Metop/AVHRR H32 vuodesta 2015 alkaen. Lisäksi esitellään pintahavaintoihin perustuvat validointitulokset, jotka osoittavat tuotteiden saavuttavan määritellyt tavoitteet.Snow cover plays a significant role in the weather and climate system, ecosystems and many human activities, such as traffic. Weather station snow observations (snow depth and state of the ground) do not provide highresolution continental or global snow coverage data. The satellite observations complement in situ observations from weather stations. Geostationary weather satellites provide observations at high temporal resolution, but the spatial resolution is low, especially in polar regions. Polarorbiting weather satellites provide better spatial resolution in polar regions with limited temporal resolution. The best detection resolution is provided by optical and infra-red radiometers onboard weather satellites. Snow cover in itself is highly variable. Also, the variability of the surface properties (such as vegetation, water bodies, topography) and changing light conditions make satellite snow detection challenging. Much of this variability is in subpixel scales, and this uncertainty creates additional challenges for the development of snow detection methods. Thus, an empirical approach may be the most practical option when developing algorithms for automatic snow detection. In this work, which is a part of the EUMETSAT-funded H SAF project, two new empirically developed snow extent products for the EUMETSAT weather satellites are presented. The geostationary MSG/SEVIRI H32 snow product has been in operational production since 2008. The polar product Metop/AVHRR H32 is available since 2015. In addition, validation results based on weather station snow observations between 2015 and 2019 are presented. The results show that both products achieve the requirements set by the H SAF

Remote Sensing Monitoring of Land Surface Temperature (LST)

This book is a collection of recent developments, methodologies, calibration and validation techniques, and applications of thermal remote sensing data and derived products from UAV-based, aerial, and satellite remote sensing. A set of 15 papers written by a total of 70 authors was selected for this book. The published papers cover a wide range of topics, which can be classified in five groups: algorithms, calibration and validation techniques, improvements in long-term consistency in satellite LST, downscaling of LST, and LST applications and land surface emissivity research

Land surface temperature and evapotranspiration estimation in the Amazon evergreen forests using remote sensing data

Amazonian tropical forests play a significant role in global water, carbon and energy cycles. Considering the relevance of this biome and the climate change projections which predict a hotter and drier climate for the region, the monitoring of the vegetation status of these forests becomes of significant importance. In this context, vegetation temperature and evapotranspiration (ET) can be considered as key variables. Vegetation temperature is directly linked with plant physiology. In addition, some studies have shown the existing relationship between this variable and the CO2 absorption capacity and biomass loss of these forests. Evapotranspiration resulting from the combined processes of transpiration and evaporation links the terrestrial water, carbon and surface energy exchanges of these forests. How this variable will response to the changing climate is critical to understand the stability of these forests. Satellite remote sensing is presented as a feasible means in order to provide accurate spatially-distributed estimates of these variables. Nevertheless, the use of satellite passive imagery for analysing this region still has some limitations being of special importance the proper cloud masking of the satellite data which becomes a difficult task due to the continuous cloud cover of the region. Under the light of the aforementioned issues, the present doctoral thesis is aimed at estimating the land surface temperature and evapotranspiration of the Amazonian tropical forests using remote sensing data. In addition, as cloud screening of satellite imagery is a critical step in the processing chain of the previous magnitudes and becomes of special importance for the study region this topic has also been included in this thesis. We have mainly focused on the use of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) which is amongst major tools for studying this region. Regarding the cloud detection topic, the potential of supervised learning algorithms for cloud masking was studied in order to overcome the cloud contamination issue of the current satellite products. Models considered were: Gaussian Naïve Bayes (GNB), Linear Discriminant Analysis (LDA), Quadratic Discriminant Analysis (QDA), Random Forests (RF), Support Vector Machine (SVM) and Multilayer Perceptron (MLP). These algorithms are able to provide a continuous measure of cloud masking uncertainty (i.e. a probability estimate of each pixel belonging to clear and cloudy class) and therefore can be used under the light of a probabilistic approach. Reference dataset (a priori knowledge) requirement was satisfied by considering the collocation of Cloud Profiling Radar (CPR) and Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) observations with MODIS sensor. Model performance was tested using three independent datasets: 1) collocated CPR/CALIOP and MODIS data, 2) MODIS manually classified images and 3) in-situ ground data. For the case of satellite image and in-situ testing, results were additionally compared to current operative MYD35 (version 6.1) and Multi-Angle Implementation of the Atmospheric Correction (MAIAC) cloud masking algorithms. These results showed that machine learning algorithms were able to improve MODIS operative cloud masking performance over the region. MYD35 and MAIAC tended to underestimate and overestimate the cloud cover, respectively. Amongst the models considered, LDA stood out as the best candidate because of its maximum accuracy (difference in Kappa coefficient of 0.293/0.155 (MYD35 /MAIAC respectively)) and minimum computational associated. Regarding the estimation of land surface temperature (LST), the aim of this study was to generate specific LST products for the Amazonian tropical forests. This goal was accomplished by using a tuned split-window (SW) equation. Validation of the LST products was obtained by direct comparison between LST estimates as derived from the algorithms and two types of different LST observations: in-situ LST (T-based validation) and LST derived from the R-based method. In addition, LST algorithms were validated using independent simulated data. In-situ LST was retrieved from two infrared radiometers (SI-100 and IR-120) and a CNR4 net radiometer, situated at Tambopata test site (12.832 S, 62.282 W) in the Peruvian Amazon. Apart from this, current satellite LST products were also validated and compared to the tuned split-window. Although we have mainly focus on MODIS LST products which derive from three different LST algorithms: split-window, day and night (DN) and Temperature Emissivity Separation (TES), we have also considered the inclusion of the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor. In addition, a first assessment of the Sea and Land Surface Temperature Radiometer (SLSTR) is presented. Validation was performed separately for daytime and nighttime conditions. For MODIS sensor, current LST products showed Root Mean Square Errors (RMSE) in LST estimations between 2 K and 3K for daytime and 1 K and 2 K for nighttime. In the best case (with a restrictive cloud screening) RMSE errors decrease to values below 2K and around 1 K, respectively. The proposed LST showed RMSE values of 1K to approximately 2 K and 0.7-1.5 K (below 1.5 K and below 1 K in the best case) for daytime and nighttime conditions, thus improving current LST MODIS products. This is also in agreement with the R-based validation results, which show a RMSE reduction of 0.7 K to 1.7 K in comparison to MODIS LST products. For the case of VIIRS sensor daytime conditions, VIIRS-TES algorithm provides the best performance with a difference of 0.2 K to around 0.3 K in RMSE regarding the split window algorithm (in the best case it reduces to 0.2 K). All VIIRS LST products considered have RMSE values between 2 K and 3K. At nighttime, however VIIRS-TES is not able to outperform the SW algorithm. A difference of 0.7 K to 0.8 K in RMSE is obtained. Contrary to MODIS and the SW LST products, VIIRS-TES tends to overestimate in-situ LST values. Regarding SLSTR sensor, the L2 product provides a better agreement with in-situ observations than the proposed algorithm (daytime difference in RMSE around 0.6 K and up 0.07 K at nighttime). In the estimation of the ET, we focused on the evaluation of four commonly used remote-sensing based ET models. These were: i) Priestley-Taylor Jet Propulsion Laboratory (PT-JPL), ii) Penman-Monteith MODIS operative parametrization (PM-Mu), iii) Surface Energy Balance System (SEBS), and iv) Satellite Application Facility on Land Surface Analysis (LSASAF). These models were forced using remote-sensing data from MODIS and two ancillary meteorological data sources: i) in-situ data extracted from Large-Scale Biosphere-Atmosphere Experiment (LBA) stations (scenario I), and ii) three reanalysis datasets (scenario II), including Modern-Era Retrospective analysis for Research and Application (MERRA-2), European Centre for Medium-range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERA-Interim), and Global Land Assimilation System (GLDAS-2.1). Performance of algorithms under the two scenarios was validated using in-situ eddy-covariance measurements. For scenario I, PT-JPL provided the best agreement with in-situ ET observations (RMSE = 0.55 mm/day, R = 0.88). Neglecting water canopy evaporation resulted in an underestimation of ET measurements for LSASAF. SEBS performance was similar to that of PT-JPL, nevertheless SEBS estimates were limited by the continuous cloud cover of the region. A physically-based ET gap-filling method was used in order to alleviate this issue. PM-Mu also with a similar performance to PT-JPL tended to overestimate in-situ ET observations. For scenario II, quality assessment of reanalysis input data demonstrated that MERRA-2, ERA-Interim and GLDAS-2.1 contain biases that impact model performance. In particular, biases in radiation inputs were found the main responsible of the observed biases in ET estimates. For the region, MERRA-2 tends to overestimate daily net radiation and incoming solar radiation. ERA-Interim tends to underestimate both variables, and GLDAS-2.1 tends to overestimate daily radiation while underestimating incoming solar radiation. Discrepancies amongst these inputs resulted in large absolute deviations in spatial patterns (deviations greater than 500 mm/year) and temporal patterns

CIRA annual report FY 2016/2017

Reporting period April 1, 2016-March 31, 2017

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