Downscaling sub-daily Land Surface Temperature time series for monitoring heat in urban environments

Author(s):  
Nikos Alexandris ◽  
Matteo Piccardo ◽  
Vasileios Syrris ◽  
Alessandro Cescatti ◽  
Gregory Duveiller

<p>The frequency of extreme heat related events is rising. This places the ever growing number of urban dwellers at higher risk. Quantifying these phenomena is important for the development and monitoring of climate change adaptation and mitigation policies. In this context, earth observations offer increasing opportunities to assess these phenomena with an unprecedented level of accuracy and spatial reach. Satellite thermal imaging systems acquire Land Surface Temperature (LST) which is fundamental to run models that study for example hotspots and heatwaves in urban environments.</p><p>Current instruments include TIRS on board Landsat 8 and MODIS on board of Terra satellites. These provide LST products on a monthly basis at 100m and twice per day at 1km respectively. Other sensors on board geostationary satellites, such as MSG and GOES-R, produce sub-hourly thermal images. For example the SEVIRI instrument onboard MSG, captures images every 15 minutes. However, this is done at an even coarser spatial resolution, which is 3 to 5 km in the case of SEVIRI. Nevertheless, none of the existing systems can capture LST synchronously with fine spatial resolution at a high temporal frequency, which is a prerequisite for monitoring heat stress in urban environments.</p><p>Combining LST time series of high temporal resolution (i.e. sub-daily MODIS- or SEVIRI-derived data) with products of fine spatial resolution (i.e. Landsat 8 products), and potentially other related variables (i.e. reflectance, spectral indices, land cover information, terrain parameters and local climatic variables), facilitates the downscaling of LST estimations. Nonetheless, considering the complexity of how distinct surfaces within a city heat-up differently during the course of a day, such a downscaling is meaningful for practically synchronous observations (e.g. Landsat-8 and MODIS Terra’s morning observations).</p><p>The recently launched ECOSTRESS mission provides multiple times in a day high spatial resolution thermal imagery at 70m. Albeit, recording the same locations on Earth every few days at varying times. We explore the associations between ECOSTRESS and Landsat-8 thermal data, based on the incoming radiation load and distinct surface properties characterised from other datasets. In our approach, first we upscale ECOSTRESS data to simulate Landsat-8 images at moments that coincide the acquisition times of other sensors products. In a second step, using the simulated Landsat-8 images, we downscale LST products acquired at later times, such as MODIS Aqua (ca. 13:30) or even the hourly MSG data. This composite downscaling procedure enables an enhanced LST estimation that opens the way for better diagnostics of the heat stress in urban landscapes.</p><p>In this study we discuss in detail the concepts of our approach and present preliminary results produced with the JEODPP, JRC's high throughput computing platform.</p>

2020 ◽  
Vol 12 (5) ◽  
pp. 791 ◽  
Author(s):  
Jingjing Yang ◽  
Si-Bo Duan ◽  
Xiaoyu Zhang ◽  
Penghai Wu ◽  
Cheng Huang ◽  
...  

Land surface temperature (LST) is vital for studies of hydrology, ecology, climatology, and environmental monitoring. The radiative-transfer-equation-based single-channel algorithm, in conjunction with the atmospheric profile, is regarded as the most suitable one with which to produce long-term time series LST products from Landsat thermal infrared (TIR) data. In this study, the performances of seven atmospheric profiles from different sources (the MODerate-resolution Imaging Spectroradiomete atmospheric profile product (MYD07), the Atmospheric Infrared Sounder atmospheric profile product (AIRS), the European Centre for Medium-range Weather Forecasts (ECMWF), the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA2), the National Centers for Environmental Prediction (NCEP)/Global Forecasting System (GFS), NCEP/Final Operational Global Analysis (FNL), and NCEP/Department of Energy (DOE)) were comprehensively evaluated in the single-channel algorithm for LST retrieval from Landsat 8 TIR data. Results showed that when compared with the radio sounding profile downloaded from the University of Wyoming (UWYO), the worst accuracies of atmospheric parameters were obtained for the MYD07 profile. Furthermore, the root-mean-square error (RMSE) values (approximately 0.5 K) of the retrieved LST when using the ECMWF, MERRA2, NCEP/GFS, and NCEP/FNL profiles were smaller than those but greater than 0.8 K when the MYD07, AIRS, and NCEP/DOE profiles were used. Compared with the in situ LST measurements that were collected at the Hailar, Urad Front Banner, and Wuhai sites, the RMSE values of the LST that were retrieved by using the ECMWF, MERRA2, NCEP/GFS, and NCEP/FNL profiles were approximately 1.0 K. The largest discrepancy between the retrieved and in situ LST was obtained for the NCEP/DOE profile, with an RMSE value of approximately 1.5 K. The results reveal that the ECMWF, MERRA2, NCEP/GFS, and NCEP/FNL profiles have great potential to perform accurate atmospheric correction and generate long-term time series LST products from Landsat TIR data by using a single-channel algorithm.


Author(s):  
F. Farhanj ◽  
M. Akhoondzadeh

Land surface temperature image is an important product in many lithosphere and atmosphere applications. This image is retrieved from the thermal infrared bands. These bands have lower spatial resolution than the visible and near infrared data. Therefore, the details of temperature variation can't be clearly identified in land surface temperature images. The aim of this study is to enhance spatial information in thermal infrared bands. Image fusion is one of the efficient methods that are employed to enhance spatial resolution of the thermal bands by fusing these data with high spatial resolution visible bands. Multi-resolution analysis is an effective pixel level image fusion approach. In this paper, we use contourlet, non-subsampled contourlet and sharp frequency localization contourlet transform in fusion due to their advantages, high directionality and anisotropy. The absolute average difference and RMSE values show that with small distortion in the thermal content, the spatial information of the thermal infrared and the land surface temperature images is enhanced.


2020 ◽  
Author(s):  
Maite Lezama Valdes ◽  
Marwan Katurji ◽  
Hanna Meyer

<p>Anthropogenic Climate Change is expected to take a toll on the Antarctic environment and its biodiversity, which is concentrated on the continent’s few ice-free areas, such as the McMurdo Dry Valleys (MDV). To model the current terrestrial habitat distribution and predict possible climate induced changes, high spatial and temporal resolution abiotic variables, especially land surface temperature (LST) and soil moisture are needed, but are currently unavailable.</p><p>The aim of this project is to fill this gap and create a high resolution LST dataset of the Antarctic Dry Valleys. This variable is acquired in a high temporal resolution (sub-daily) by the MODIS sensor aboard Terra and Aqua satellites. However, as LST varies greatly in space, the spatial resolution provided by this data source (1000 m) is too low to give a meaningful impression of LST and to study biodiversity patterns. Therefore, we use data from Landsat and ASTER sensors as a reference to downscale MODIS LST to a spatial resolution of 30 m. 7 year’s worth of satellite data as well as terrain-derived auxiliary variables went into the development of the model, which predicts 30 m LST for the Antarctic Dry Valleys. </p><p>To model complex relations between terrain, radiation, land cover and LST, machine learning models are used. Multiple algorithms (Random Forest, NN, SVM, Gradient Boosting) are compared to find the best approach for predicting high resolution LST based on MODIS data. Using the best performing model, a daily dataset is created that provides LST for the Antarctic Dry Valleys from 2002 on.</p>


2019 ◽  
Vol 11 (11) ◽  
pp. 1319 ◽  
Author(s):  
Paulina Bartkowiak ◽  
Mariapina Castelli ◽  
Claudia Notarnicola

In this study, we evaluated three different downscaling approaches to enhance spatial resolution of thermal imagery over Alpine vegetated areas. Due to the topographical and land-cover complexity and to the sparse distribution of meteorological stations in the region, the remotely-sensed land surface temperature (LST) at regional scale is of major area of interest for environmental applications. Even though the Moderate Resolution Imaging Spectroradiometer (MODIS) LST fills the gap regarding high temporal resolution and length of the time-series, its spatial resolution is not adequate for mountainous areas. Given this limitation, random forest algorithm for downscaling LST to 250 m spatial resolution was evaluated. This study exploits daily MODIS LST with a spatial resolution of 1 km to obtain sub-pixel information at 250 m spatial resolution. The nonlinear relationship between coarse resolution MODIS LST (CR) and fine resolution (FR) explanatory variables was performed by building three different models including: (i) all pixels (BM), (ii) only pixels with more than 90% of vegetation content (EM1) and (iii) only pixels with 75% threshold of homogeneity for vegetated land-cover classes (EM2). We considered normalized difference vegetation index (NDVI) and digital elevation model (DEM) as predictors. The performances of the thermal downscaling methods were evaluated by the Root Mean Square Error (RMSE) and the Mean Absolute Error (MAE) between the downscaled dataset and Landsat LST. Validation indicated that the error values for vegetation fraction (EM1, EM2) were smaller than for basic modelling (BM). BM model determined averaged RMSE of 2.3 K and MAE of 1.8 K. Enhanced methods (EM1 and EM2) gave slightly better results yielding 2.2 K and 1.7 K for RMSE and MAE, respectively. In contrast to the EMs, BM showed a reduction of 22% and 18% of RMSE and MAE respectively with regard to Landsat and the original MODIS LST. Despite some limitations, mainly due to cloud contamination effect and coarse resolution pixel heterogeneity, random forest downscaling exhibits a large potential for producing improved LST maps.


Sensors ◽  
2020 ◽  
Vol 20 (18) ◽  
pp. 5336
Author(s):  
Sorin Cheval ◽  
Alexandru Dumitrescu ◽  
Vlad-Alexandru Amihaesei

The Landsat 8 satellites have retrieved land surface temperature (LST) resampled at a 30-m spatial resolution since 2013, but the urban climate studies frequently use a limited number of images due to the problems related to missing data over the city of interest. This paper endorses a procedure for building a long-term gap-free LST data set in an urban area using the high-resolution Landsat 8 imagery. The study is applied on 94 images available through 2013–2018 over Bucharest (Romania). The raw images containing between 1.1% and 58.4% missing LST data were filled in using the Data INterpolating Empirical Orthogonal Functions (DINEOF) algorithm implemented in the sinkr R packages. The resulting high-spatial-resolution gap-filled land surface temperature data set was used to explore the LST climatology over Bucharest (Romania) an urban area, at a monthly, seasonal, and annual scale. The performance of the gap-filling method was checked using a cross-validation procedure, and the results pledge for the development of an LST-based urban climatology.


2017 ◽  
Vol 9 (1) ◽  
Author(s):  
Andrzej Chybicki ◽  
Zbigniew Łubniewski

AbstractSatellite imaging systems have known limitations regarding their spatial and temporal resolution. The approaches based on subpixel mapping of the Earth’s environment, which rely on combining the data retrieved from sensors of higher temporal and lower spatial resolution with the data characterized by lower temporal but higher spatial resolution, are of considerable interest. The paper presents the downscaling process of the land surface temperature (LST) derived from low resolution imagery acquired by the Advanced Very High Resolution Radiometer (AVHRR), using the inverse technique. The effective emissivity derived from another data source is used as a quantity describing thermal properties of the terrain in higher resolution, and allows the downsampling of low spatial resolution LST images. The authors propose an optimized downscaling method formulated as the inverse problem and show that the proposed approach yields better results than the use of other downsampling methods. The proposed method aims to find estimation of high spatial resolution LST data by minimizing the global error of the downscaling. In particular, for the investigated region of the Gulf of Gdansk, the RMSE between the AVHRR image downscaled by the proposed method and the Landsat 8 LST reference image was 2.255°C with correlation coefficient R equal to 0.828 and Bias = 0.557°C. For comparison, using the PBIM method, it was obtained RMSE = 2.832°C, R = 0.775 and Bias = 0.997°C for the same satellite scene. It also has been shown that the obtained results are also good in local scale and can be used for areas much smaller than the entire satellite imagery scene, depicting diverse biophysical conditions. Specifically, for the analyzed set of small sub-datasets of the whole scene, the obtained RSME between the downscaled and reference image was smaller, by approx. 0.53°C on average, in the case of applying the proposed method than in the case of using the PBIM method.


2019 ◽  
Vol 11 (21) ◽  
pp. 2588 ◽  
Author(s):  
Otgonbayar ◽  
Atzberger ◽  
Mattiuzzi ◽  
Erdenedalai

The objective of this research was to develop a robust statistical model to estimate climatologies (2002–2017) of monthly average near-surface air temperature (Ta) over Mongolia using Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) time series products and terrain parameters. Two regression models were analyzed in this study linking automatic weather station data (Ta) with Earth observation (EO) images: partial least squares (PLS) and random forest (RF). Both models were trained to predict Ta climatologies for each of the twelve months, using up to 17 variables as predictors. The models were applied to the entire land surface of Mongolia, the eighteenth largest but most sparsely populated country in the world. Twelve of the predictor variables were derived from the LST time series products of the Terra MODIS satellite. The LST MOD11A2 (collection 6) products provided thermal information at a spatial resolution of 1 km and with 8-day temporal resolution from 2002 to 2017. Three terrain variables, namely, elevation, slope, and aspect, were extracted using a Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM), and two variables describing the geographical location of weather stations were extracted from vector data. For training, a total of 8544 meteorological data points from 63 automatic weather stations were used covering the same period as MODIS LST products. The PLS regression resulted in a coefficient of determination (R2) between 0.74 and 0.87 and a root-mean-square error (RMSE) from 1.20 °C to 2.19 °C between measured and estimated monthly Ta. The non-linear RF regression yielded even more accurate results with R2 in the range from 0.82 to 0.95 and RMSE from 0.84 °C to 1.93 °C. Using RF, the two best modeled months were July and August and the two worst months were January and February. The four most predictive variables were day/nighttime LST, elevation, and latitude. Using the developed RF models, spatial maps of the monthly average Ta at a spatial resolution of 1 km were generated for Mongolia (~1566 × 106 km2). This spatial dataset might be useful for various environmental applications. The method is transparent and relatively easy to implement.


Sign in / Sign up

Export Citation Format

Share Document