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

Feasibility of Using the Two-Source Energy Balance Model (TSEB) with Sentinel-2 and Sentinel-3 Images to Analyze the Spatio-Temporal Variability of Vine Water Status in a Vineyard

In viticulture, detailed spatial information about actual evapotranspiration (ETa) and vine water status within a vineyard may be of particular utility when applying site-specific, precision irrigation management. Over recent decades, extensive research has been carried out in the use of remote sensing energy balance models to estimate and monitor ETa at the field level. However, one of the major limitations remains the coarse spatial resolution in the thermal infrared (TIR) domain. In this context, the recent advent of the Sentinel missions of the European Space Agency (ESA) has greatly improved the possibility of monitoring crop parameters and estimating ETa at higher temporal and spatial resolutions. In order to bridge the gap between the coarse-resolution Sentinel-3 thermal and the fine-resolution Sentinel-2 shortwave data, sharpening techniques have been used to downscale the Sentinel-3 land surface temperature (LST) from 1 km to 20 m. However, the accurate estimates of high-resolution LST through sharpening techniques are still unclear, particularly when intended to be used for detecting crop water stress. The goal of this study was to assess the feasibility of the two-source energy balance model (TSEB) using sharpened LST images from Sentinel-2 and Sentinel-3 (TSEB-PTS2+3) to estimate the spatio-temporal variability of actual transpiration (T) and water stress in a vineyard. T and crop water stress index (CWSI) estimates were evaluated against a vine water consumption model and regressed with in situ stem water potential (Ψstem). Two different TSEB approaches, using very high-resolution airborne thermal imagery, were also included in the analysis as benchmarks for TSEB-PTS2+3. One of them uses aggregated TIR data at the vine+inter-row level (TSEB-PTairb), while the other is based on a contextual method that directly, although separately, retrieves soil and canopy temperatures (TSEB-2T). The results obtained demonstrated that when comparing airborne Trad and sharpened S2+3 LST, the latter tend to be underestimated. This complicates the use of TSEB-PTS2+3 to detect crop water stress. TSEB-2T appeared to outperform all the other methods. This was shown by a higher R2 and slightly lower RMSD when compared with modelled T. In addition, regressions between T and CWSI-2T with Ψstem also produced the highest R2.info:eu-repo/semantics/publishedVersio

Mapping debris-covered glaciers in the Cordillera Blanca, Peru : an object-based image analysis approach.

Accurate remote-sensing based inventories of glacial ice are often hindered by the presence of supraglacial debris cover. Attempts at automated mapping of debris-covered glacier areas from remotely-sensed multispectral data have met with limited success due to the spectral similarity of supraglacial debris to nearby bedrock, moraines, and fluvial deposition features. Data-fusion approaches leveraging terrain and/or thermal data with multispectral data have yielded improved results in certain geographic regions, but remain unproven in others. This research builds on the data-fusion approaches from the literature and explores the efficacy of object-based image analysis (OBIA) and tree-based machine learning classifiers using Landsat OLI imagery and SRTM elevation data, in effort to map debris-covered glaciers in the Cordillera Blanca range of Peru. Results suggest that the OBIA and machine learning methods render advantages over traditional methods given the unique morphological settings associated with debris-covered glaciers. Accurate inventories of glacial mass and debris-covered glaciers in the Cordillera Blanca are important for understanding the unique water resource, natural hazards, and climate change implications associated with these tropical mountain glaciers

Multispectral pansharpening with radiative transfer-based detail-injection modeling for preserving changes in vegetation cover

Whenever vegetated areas are monitored over time, phenological changes in land cover should be decoupled from changes in acquisition conditions, like atmospheric components, Sun and satellite heights and imaging instrument. This especially holds when the multispectral (MS) bands are sharpened for spatial resolution enhancement by means of a panchromatic (Pan) image of higher resolution, a process referred to as pansharpening. In this paper, we provide evidence that pansharpening of visible/near-infrared (VNIR) bands takes advantage of a correction of the path radiance term introduced by the atmosphere, during the fusion process. This holds whenever the fusion mechanism emulates the radiative transfer model ruling the acquisition of the Earth's surface from space, that is for methods exploiting a multiplicative, or contrast-based, injection model of spatial details extracted from the panchromatic (Pan) image into the interpolated multispectral (MS) bands. The path radiance should be estimated and subtracted from each band before the product by Pan is accomplished. Both empirical and model-based estimation techniques of MS path radiances are compared within the framework of optimized algorithms. Simulations carried out on two GeoEye-1 observations of the same agricultural landscape on different dates highlight that the de-hazing of MS before fusion is beneficial to an accurate detection of seasonal changes in the scene, as measured by the normalized differential vegetation index (NDVI)

Downscaling Coarse Resolution Satellite Passive Microwave SWE Estimates

The spatio-temporal heterogeneity of seasonal snow and its impact on socio-economic and environmental functionality make accurate, real-time estimates of snow water equivalent (SWE) important for hydrological and climatological predictions. Passive microwave remote sensing offers a cost effective, temporally and spatially consistent approach to SWE monitoring at the global to regional scale. However, local scale estimates are subject to large errors given the coarse spatial resolution of passive microwave observations (25 x 25 km). Regression downscaling techniques can be implemented to increase the spatial resolution of gridded datasets with the use of related auxiliary datasets at a finer spatial resolution. These techniques have been successfully implemented to remote sensing datasets such as soil moisture estimates, however, limited work has applied such techniques to snow-related datasets. This thesis focuses on assessing the feasibility of using regression downscaling to increase the spatial resolution of the European Space Agency’s (ESA) Globsnow SWE product in the Red River basin, an agriculturally important region of the northern United States that is widely recognized as a suitable location for passive microwave remote sensing research. Multiple Linear (MLR), Random Forest (RFR) and Geographically Weighted (GWR) regression downscaling techniques were assessed in a closed loop experiment using Snow Data Assimilation System (SNODAS) SWE estimates at a 1 x 1 km spatial resolution. SNODAS SWE data for a 5-year period between 2013-2018 was aggregated to a 25 x 25 km spatial resolution to match Globsnow. The three regression techniques were applied using correlative datasets to downscale the aggregated SNODAS data back to the original 1 x 1 km spatial resolution. By comparing the downscaled SNODAS estimates to the original SNODAS data, it was found that RFR downscaling produced much less variation in downscaled results, and lower RMSE values throughout the study period when compared to MLR and GWR downscaling techniques, indicating it was the optimal downscaling method. RFR downscaling was then implemented on daily Globsnow SWE estimates for the same time period. The downscaled SWE results were evaluated using SNODAS SWE as well as in situ derived SWE estimates from weather stations within the study region. Spatial and temporal errors were assessed using both the SNODAS and in situ reference datasets and overall RMSEs of 21 mm and 37 mm were found, respectively. It was observed that the southern regions of the basin and seasons with higher downscaled SWE estimates were associated with higher errors with overestimation being the most common bias throughout the region. A major contribution of this study is the illustration that RFR downscaling of Globsnow SWE estimates is a feasible approach to understanding the seasonal dynamics of SWE in the Red River basin. This is extremely beneficial for local communities within the basin for flood management and mitigation and water resource management

Relevance of UAV and sentinel-2 data fusion for estimating topsoil organic carbon after forest fire

[EN] The evaluation at detailed spatial scale of soil status after severe fires may provide useful information on the recovery of burned forest ecosystems. Here, we aimed to assess the potential of combining multispectral imagery at different spectral and spatial resolutions to estimate soil indicators of burn severity. The study was conducted in a burned area located at the northwest of the Iberian Peninsula (Spain). One month after fire, we measured soil burn severity in the field using an adapted protocol of the Composite Burn Index (CBI). Then, we performed soil sampling to analyze three soil properties potentially indicatives of fire-induced changes: mean weight diameter (MWD), soil moisture content (SMC) and soil organic carbon (SOC). Additionally, we collected post-fire imagery from the Sentinel-2A MSI satellite sensor (10–20 m of spatial resolution), as well as from a Parrot Sequoia camera on board an unmanned aerial vehicle (UAV; 0.50 m of spatial resolution). A Gram-Schmidt (GS) image sharpening technique was used to increase the spatial resolution of Sentinel-2 bands and to fuse these data with UAV information. The performance of soil parameters as indicators of soil burn severity was determined trough a machine learning decision tree, and the relationship between soil indicators and reflectance values (UAV, Sentinel-2 and fused UAV-Sentinel-2 images) was analyzed by means of support vector machine (SVM) regression models. All the considered soil parameters decreased their value with burn severity, but soil moisture content, and, to a lesser extent, soil organic carbon discriminated at best among soil burn severity classes (accuracy = 91.18 %; Kappa = 0.82). The performance of reflectance values derived from the fused UAV-Sentinel-2 image to monitor the effects of wildfire on soil characteristics was outstanding, particularly for the case of soil organic carbon content (R2 = 0.52; RPD = 1.47). This study highlights the advantages of combining satellite and UAV images to produce spatially and spectrally enhanced images, which may be relevant for estimating main impacts on soil properties in burned forest areas where emergency actions need to be applied.S

Subpixel temperature estimation from single-band thermal infrared imagery

Target temperature estimation from thermal infrared (TIR) imagery is a complex task that becomes increasingly more difficult as the target size approaches the size of a projected pixel. At that point the assumption of pixel homogeneity is invalid as the radiance value recorded at the sensor is the result of energy contributions from the target material and any other background material that falls within a pixel boundary. More often than not, thermal infrared pixels are heterogeneous and therefore subpixel temperature extraction becomes an important capability. Typical subpixel estimation approaches make use of multispectral or hyperspectral sensors. These technologies are expensive and multispectral or hyperspectral thermal imagery might not be readily available for a target of interest. A methodology was developed to retrieve the temperature of an object that is smaller than a projected pixel of a single-band TIR image using physics-based modeling. Physics-based refers to the utilization of the Multi-Service Electro-optic Signature (MuSES) heat transfer model, the MODerate spectral resolution atmospheric TRANsmission (MODTRAN) atmospheric propagation algorithm, and the Digital Imaging and Remote Sensing Image Generation (DIRSIG) synthetic image generation model to reproduce a collected thermal image under a number of user-supplied conditions. A target space is created and searched to determine the temperature of the subpixel target of interest from a collected TIR image. The methodology was tested by applying it to single-band thermal imagery collected during an airborne campaign. The emissivity of the targets of interest ranged from 0.02 to 0.91 and the temperature extraction error for the high emissivity targets were similar to the temperature extraction errors found in published papers that employed multi-band techniques

Evaluation of two applications of spectral mixing models to image fusion

None provided

A PCA-OLS Model for Assessing the Impact of Surface Biophysical Parameters on Land Surface Temperature Variations

Analysis of land surface temperature (LST) spatiotemporal variations and characterization of the factors affecting these variations are of great importance in various environmental studies and applications. The aim of this study is to propose an integrated model for characterizing LST spatiotemporal variations and for assessing the impact of surface biophysical parameters on the LST variations. For this purpose, a case study was conducted in Babol City, Iran, during the period of 1985 to 2018. We used 122 images of Landsat 5, 7, and 8, and products of water vapor (MOD07) and daily LST (MOD11A1) from the MODIS sensor of the Terra satellite, as well as soil and air temperature and relative humidity data measured at the local meteorological station over 112 dates for the study. First, a single-channel algorithm was applied to estimate LST, while various spectral indices were computed to represent surface biophysical parameters, which included the normalized difference vegetation index (NDVI), soil-adjusted vegetation index (SAVI), normalized difference water index (NDWI), normalized difference built-up index (NDBI), albedo, brightness, greenness, and wetness from tasseled cap transformation. Next, a principal component analysis (PCA) was conducted to determine the degree of LST variation and the surface biophysical parameters in the temporal dimension at the pixel scale based on Landsat imagery. Finally, the relationship between the first component of the PCA of LST and each surface biophysical parameter was investigated by using the ordinary least squares (OLS) regression with both regional and local optimizations. The results indicated that among the surface biophysical parameters, variations of NDBI, wetness, and greenness had the highest impact on the LST variations with a correlation coefficient of 0.75, −0.70, and −0.44, and RMSE of 0.71, 1.03, and 1.06, respectively. The impact of NDBI, wetness, and greenness varied geographically, but their variations accounted for 43%, 38%, and 19% of the LST variation, respectively. Furthermore, the correlation coefficient and RMSE between the observed LST variation and modeled LST variation, based on the most influential biophysical factors (NDBI, wetness, and greenness) yielded 0.85 and 1.06 for the regional approach and 0.93 and 0.26 for the local approach, respectively. The results of this study indicated the use of an integrated PCA–OLS model was effective for modeling of various environmental parameters and their relationship with LST. In addition, the PCA–OLS with the local optimization was found to be more efficient than the one with the regional optimization

Enhancing spatial resolution of remotely sensed data for mapping freshwater environments

Freshwater environments are important for ecosystem services and biodiversity. These environments are subject to many natural and anthropogenic changes, which influence their quality; therefore, regular monitoring is required for their effective management. High biotic heterogeneity, elongated land/water interaction zones, and logistic difficulties with access make field based monitoring on a large scale expensive, inconsistent and often impractical. Remote sensing (RS) is an established mapping tool that overcomes these barriers. However, complex and heterogeneous vegetation and spectral variability due to water make freshwater environments challenging to map using remote sensing technology. Satellite images available for New Zealand were reviewed, in terms of cost, and spectral and spatial resolution. Particularly promising image data sets for freshwater mapping include the QuickBird and SPOT-5. However, for mapping freshwater environments a combination of images is required to obtain high spatial, spectral, radiometric, and temporal resolution. Data fusion (DF) is a framework of data processing tools and algorithms that combines images to improve spectral and spatial qualities. A range of DF techniques were reviewed and tested for performance using panchromatic and multispectral QB images of a semi-aquatic environment, on the southern shores of Lake Taupo, New Zealand. In order to discuss the mechanics of different DF techniques a classification consisting of three groups was used - (i) spatially-centric (ii) spectrally-centric and (iii) hybrid. Subtract resolution merge (SRM) is a hybrid technique and this research demonstrated that for a semi aquatic QuickBird image it out performed Brovey transformation (BT), principal component substitution (PCS), local mean and variance matching (LMVM), and optimised high pass filter addition (OHPFA). However some limitations were identified with SRM, which included the requirement for predetermined band weights, and the over-representation of the spatial edges in the NIR bands due to their high spectral variance. This research developed three modifications to the SRM technique that addressed these limitations. These were tested on QuickBird (QB), SPOT-5, and Vexcel aerial digital images, as well as a scanned coloured aerial photograph. A visual qualitative assessment and a range of spectral and spatial quantitative metrics were used to evaluate these modifications. These included spectral correlation and root mean squared error (RMSE), Sobel filter based spatial edges RMSE, and unsupervised classification. The first modification addressed the issue of predetermined spectral weights and explored two alternative regression methods (Least Absolute Deviation, and Ordinary Least Squares) to derive image-specific band weights for use in SRM. Both methods were found equally effective; however, OLS was preferred as it was more efficient in processing band weights compared to LAD. The second modification used a pixel block averaging function on high resolution panchromatic images to derive spatial edges for data fusion. This eliminated the need for spectral band weights, minimised spectral infidelity, and enabled the fusion of multi-platform data. The third modification addressed the issue of over-represented spatial edges by introducing a sophisticated contrast and luminance index to develop a new normalising function. This improved the spatial representation of the NIR band, which is particularly important for mapping vegetation. A combination of the second and third modification of SRM was effective in simultaneously minimising the overall spectral infidelity and undesired spatial errors for the NIR band of the fused image. This new method has been labelled Contrast and Luminance Normalised (CLN) data fusion, and has been demonstrated to make a significant contribution in fusing multi-platform, multi-sensor, multi-resolution, and multi-temporal data. This contributes to improvements in the classification and monitoring of fresh water environments using remote sensing