Comparison of permafrost mean annual ground temperature derived from two different satellite-based schemes: land surface temperature based (ESA CCI+ Permafrost) versus surface status (Metop ASCAT)

Author(s):  
Dan Jakober ◽  
Helena Bergstedt ◽  
Christine Kroisleitner ◽  
Annett Bartsch

<p>Different approaches exist for a satellite-based estimation of mean annual ground temperature (MAGT). Landsurface temperature can be ingested by transient models. Surface status information (frozen/unfrozen days) has been shown to be applicable for the estimation of ground temperature as well. Such approaches are based on an empirically defined relationship. Both approaches have been evaluated with in situ bore hole measurements, but not yet compared with each other.</p><p>A comparison between yearly arctic mean temperatures, derived from the advanced scatterometer (ASCAT) and data from ESA’s CCI+ Permafrost project was carried out. The used ASCAT record is available from 2008 (first full year) onwards while the latest CCI+ Permafrost data is available from 1997 to 2018. The ASCAT data was recorded by satellites whose measurements are only intermittently available as one flyover over the whole arctic north of 60°N takes two days on average. To fill in the missing values exponentially weighted moving averages (EWMA) were used. From the number of frozen days an expected average temperature was derived based on Kroisleitner et al. (2018).</p><p>The CCI+ Permafrost data incorporates modelled MAGT for depths between the surface down to a depth of 10 meters. These data points were extracted from the raster files (~1km resolution) and averaged over polygons representing an approximation of the ASCAT grid (footprint approximation). Single polygon areas range from 150-160 km². Only footprints for which data is available in both records (and thus permafrost presence) have been eventually compared.</p><p>The CCI+ Permafrost data shows an average surface temperature of -1.42 °C for the areas analyzed between 2008 and 2018 while the statistically padded ASCAT data suggests a mean temperature of -1.18 °C over the same time period. The ASCAT retrieval corresponds to a general MAGT whereas CCI+ Permafrost values are available for certain depths. Water fraction within ASCAT footprint also affect the quality of the derivation of frozen days. New calibration considering certain depths and water fraction is suggested.</p>

2021 ◽  
Author(s):  
Alejandro Corbea-Pérez ◽  
Gonçalo Vieira ◽  
Carmen Recondo ◽  
Joana Baptista ◽  
Javier F.Calleja ◽  
...  

<p>Land surface temperature is an important factor for permafrost modelling as well as for understanding the dynamics of Antarctic terrestrial ecosystems (Bockheim et al. 2008). In the South Shetland Islands the distribution of permafrost is complex (Vieira et al. 2010) and the use of remote sensing data is essential since the installation and maintenance of an extensive network of ground-based stations are impossible. Therefore, it is important to evaluate the applicability of satellites and sensors by comparing data with in-situ observations. In this work, we present the results from the analysis of land surface temperatures from Barton Peninsula, an ice-free area in King George Island (South Shetlands). We have studied the period from March 1, 2019 to January 31, 2020 using data from the Moderate Resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature (LST) and in-situ data from 6 ground temperature loggers. MOD11A1 and MYD11A1 products, from TERRA and AQUA satellites, respectively, were used, following the application of MODIS quality filters. Given the scarce number of high-quality data as defined by MODIS, all average LST with error ≤ 2K were included. Dates with surface temperature below -20ºC, which are rare in the study area, and dates when the difference between MODIS and in-situ data exceeded 10ºC were also examined. In both cases, those days on which MOD09GA/MYD09GA products showed cloud cover were eliminated. Eight in-situ ground temperature measurements per day were available, from which the one nearest to the time of satellite overpass was selected for comparison with MODIS-LST. The results obtained show a better correlation with daytime data than with nighttime data. Specifically, the best results are obtained with daytime data from AQUA (R<sup>2</sup> between 0.55 and 0.81). With daytime data, correlation between MODIS-LST and in-situ data was verified with relative humidity (RH) values provided by King Sejong weather station, located in the study area. When RH is lower, the correlation between LST and in-situ data improves: we obtained correlation coefficients between 0.6 - 0.7 for TERRA data and 0.8 - 0.9 for AQUA data with RH values lower than 80%. The results suggest that MODIS can be used for temperature estimation in the ice-free areas of the Maritime Antarctic.</p><p>References:</p><p>Bockheim, J. G., Campbell, I. B., Guglielmin, M., and López- Martınez, J.: Distribution of permafrost types and buried ice in ice free areas of Antarctica, in: 9th International Conference on Permafrost, 28 June–3 July 2008, Proceedings, University of Alaska Press, Fairbanks, USA, 2008, 125–130.</p><p>Vieira, G.; Bockheim, J.; Guglielmin, M.; Balks, M.; Abramov, A. A.; Boelhouwers, J.; Cannone, N.; Ganzert, L.; Gilichinsky, D. A.; Goryachkin, S.; López-Martínez, J.; Meiklejohn, I.; Raffi, R.; Ramos, M.; Schaefer, C.; Serrano, E.; Simas, F.; Sletten, R.; Wagner, D. Thermal State of Permafrost and Active-layer Monitoring in the Antarctic: Advances During the International Polar Year 2007-2009. Permafr. Periglac. Process. 2010, 21, 182–197.</p><p> </p><p>Acknowledgements</p><p>This work was made possible by an internship at the IGOT, University of Lisbon, Portugal, funded by the Principality of Asturias (code EB20-16).</p><p> </p>


2021 ◽  
Author(s):  
Gitanjali Thakur ◽  
Stan Schymanski ◽  
Kaniska Mallick ◽  
Ivonne Trebs

<p>The surface energy balance (SEB) is defined as the balance between incoming energy from the sun and outgoing energy from the Earth’s surface. All components of the SEB depend on land surface temperature (LST). Therefore, LST is an important state variable that controls the energy and water exchange between the Earth’s surface and the atmosphere. LST can be estimated radiometrically, based on the infrared radiance emanating from the surface. At the landscape scale, LST is derived from thermal radiation measured using  satellites.  At the plot scale, eddy covariance flux towers commonly record downwelling and upwelling longwave radiation, which can be inverted to retrieve LST  using the grey body equation :<br>             R<sub>lup</sub> = εσ T<sub>s</sub><sup>4</sup> + (1 − ε) R<sub> ldw         </sub>(1)<br>where R<sub>lup</sub> is the upwelling longwave radiation, R<sub>ldw</sub> is the downwelling longwave radiation, ε is the surface emissivity, <em>T<sub>s</sub>  </em>is the surface temperature and σ  is the Stefan-Boltzmann constant. The first term is the temperature-dependent part, while the second represents reflected longwave radiation. Since in the past downwelling longwave radiation was not measured routinely using flux towers, it is an established practice to only use upwelling longwave radiation for the retrieval of plot-scale LST, essentially neglecting the reflected part and shortening Eq. 1 to:<br>               R<sub>lup</sub> = εσ T<sub>s</sub><sup>4 </sup>                       (2)<br>Despite  widespread availability of downwelling longwave radiation measurements, it is still common to use the short equation (Eq. 2) for in-situ LST retrieval. This prompts the question if ignoring the downwelling longwave radiation introduces a bias in LST estimations from tower measurements. Another associated question is how to obtain the correct ε needed for in-situ LST retrievals using tower-based measurements.<br>The current work addresses these two important science questions using observed fluxes at eddy covariance towers for different land cover types. Additionally, uncertainty in retrieved LST and emissivity due to uncertainty in input fluxes was quantified using SOBOL-based uncertainty analysis (SALib). Using landscape-scale emissivity obtained from satellite data (MODIS), we found that the LST  obtained using the complete equation (Eq. 1) is 0.5 to 1.5 K lower than the short equation (Eq. 2). Also, plot-scale emissivity was estimated using observed sensible heat flux and surface-air temperature differences. Plot-scale emissivity obtained using the complete equation was generally between 0.8 to 0.98 while the short equation gave values between 0.9 to 0.98, for all land cover types. Despite additional input data for the complete equation, the uncertainty in plot-scale LST was not greater than if the short equation was used. Landscape-scale daytime LST obtained from satellite data (MODIS TERRA) were strongly correlated with our plot-scale estimates, but on average higher by 0.5 to 9 K, regardless of the equation used. However, for most sites, the correspondence between MODIS TERRA LST and retrieved plot-scale LST estimates increased significantly if plot-scale emissivity was used instead of the landscape-scale emissivity obtained from satellite data.</p>


2013 ◽  
Vol 785-786 ◽  
pp. 1333-1336
Author(s):  
Xiao Feng Yang ◽  
Xing Ping Wen

Land surface temperature (LST) is important factor in global climate change studies, radiation budgets estimating, city heat and others. In this paper, land surface temperature of Guangzhou metropolis was retrieved from two MODIS imageries obtained at night and during the day respectively. Firstly, pixel values were calibrated to spectral radiances according to parameters from header files. Then, the brightness temperature was calculated using Planck function. Finally, The brightness temperature retrieval maps were projected and output. Comparing two brightness temperature retrieval maps, it is concluded that the brightness temperature retrieval are more accurate at night than during the day. Comparing the profile line of brightness temperature from north to south, the brightness temperature increases from north to south. Temperature different from north to south is larger at night than during the day. The average temperature nears 18°C at night and the average temperature nears 26°C during the day, which is consistent with the surface temperature observed by automatic weather stations.


2016 ◽  
Author(s):  
H. S. Benavides Pinjosovsky ◽  
S. Thiria ◽  
C. Ottlé ◽  
J. Brajard ◽  
F. Badran ◽  
...  

Abstract. The SECHIBA module of the ORCHIDEE land surface model describes the exchanges of water and energy between the surface and the atmosphere. In the present paper, the adjoint semi-generator software denoted YAO was used as a framework to implement a 4D-VAR assimilation method. The objective was to deliver the adjoint model of SECHIBA (SECHIBA-YAO) obtained with YAO to provide an opportunity for scientists and end users to perform their own assimilation. SECHIBA-YAO allows the control of the eleven most influent internal parameters of SECHIBA or of the initial conditions of the soil water content by observing the land surface temperature measured in situ or as it could be observed by remote sensing as brightness temperature. The paper presents the fundamental principles of the 4D-Var assimilation, the semi-generator software YAO and some experiments showing the accuracy of the adjoint code distributed. In addition, a distributed version is available when only the land surface temperature is observed.


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.


2019 ◽  
Vol 225 ◽  
pp. 16-29 ◽  
Author(s):  
Si-Bo Duan ◽  
Zhao-Liang Li ◽  
Hua Li ◽  
Frank-M. Göttsche ◽  
Hua Wu ◽  
...  

Sensors ◽  
2020 ◽  
Vol 20 (6) ◽  
pp. 1778 ◽  
Author(s):  
Md Qutub Uddin Sajib ◽  
Tao Wang

The presence of two thermal bands in Landsat 8 brings the opportunity to use either one or both of these bands to retrieve Land Surface Temperature (LST). In order to compare the performances of existing algorithms, we used four methods to retrieve LST from Landsat 8 and made an intercomparison among them. Apart from the direct use of the Radiative Transfer Equation (RTE), Single-Channel Algorithm and two Split-Window Algorithms were used taking an agricultural region in Bangladesh as the study area. The LSTs retrieved in the four methods were validated in two ways: first, an indirect validation against reference LST, which was obtained in the Atmospheric and Topographic CORection (ATCOR) software module; second, cross-validation with Terra MODerate Resolution Imaging Spectroradiometer (MODIS) daily LSTs that were obtained from the Application for Extracting and Exploring Analysis Ready Samples (A ρ ρ EEARS) online tool. Due to the absence of LST-monitoring radiosounding instruments surrounding the study area, in situ LSTs were not available; hence, validation of satellite retrieved LSTs against in situ LSTs was not performed. The atmospheric parameters necessary for the RTE-based method, as well as for other methods, were calculated from the National Centers for Environmental Prediction (NCEP) database using an online atmospheric correction calculator with MODerate resolution atmospheric TRANsmission (MODTRAN) codes. Root-mean-squared-error (RMSE) against reference LST, as well as mean bias error against both reference and MODIS daily LSTs, was used to interpret the relative accuracy of LST results. All four methods were found to result in acceptable LST products, leaving atmospheric water vapor content (w) as the important determinant for the precision result. Considering a set of several Landsat 8 images of different dates, Jiménez-Muñoz et al.’s (2014) Split-Window algorithm was found to result in the lowest mean RMSE of 1.19 ° C . Du et al.’s (2015) Split-Window algorithm resulted in mean RMSE of 1.50 ° C . The RTE-based direct method and the Single-Channel algorithm provided the mean RMSE of 2.47 ° C and 4.11 ° C , respectively. For Du et al.’s algorithm, the w range of 0.0 to 6.3 g cm−2 was considered, whereas for the other three methods, w values as retrieved from the NCEP database were considered for corresponding images. Land surface emissivity was retrieved through the Normalized Difference Vegetation Index (NDVI)-threshold method. This intercomparison study provides an LST retrieval methodology for Landsat 8 that involves four algorithms. It proves that (i) better LST results can be obtained using both thermal bands of Landsat 8; (ii) the NCEP database can be used to determine atmospheric parameters using the online calculator; (iii) MODIS daily LSTs from A ρ ρ EEARS can be used efficiently in cross-validation and intercomparison of Landsat 8 LST algorithms; and (iv) when in situ LST data are not available, the ATCOR-derived LSTs can be used for indirect verification and intercomparison of Landsat 8 LST algorithms.


2018 ◽  
Vol 10 (11) ◽  
pp. 1852 ◽  
Author(s):  
Lei Lu ◽  
Tingjun Zhang ◽  
Tiejun Wang ◽  
Xiaoming Zhou

Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) products are widely used in ecology, hydrology, vegetation monitoring, and global circulation models. Compared to the collection-5 (C5) LST products, the newly released collection-6 (C6) LST products have been refined over bare soil pixels. This study aims to evaluate the C6 MODIS 1-km LST product using multi-year in situ data covering barren surfaces. Evaluation using all in situ data shows that the MODIS C6 LSTs are underestimated with a root-mean-square error (RMSE) of 2.59 K for the site in the Gobi area, 3.05 K for the site in the sand desert area, and 2.86 K for the site in the desert steppe area at daytime. For nighttime LSTs, the RMSEs are 2.01 K, 2.88 K, and 1.80 K for the three sites, respectively. Both biases and RMSEs also show strong seasonal signals. Compared to the error of C5 1-km LSTs, the RMSE of C6 1-km LST product is smaller, especially for daytime LSTs, with a value of 2.24 K compared to 3.51 K. The large errors in the sand desert region are presumably due to the lack of global representativeness of the magnitude of emissivity adjustment and misclassification for the barren surface causing error in emissivities. It indicates that the accuracy of the MODIS C6 LST product might be further improved through emissivity adjustment with globally representative magnitude and accurate land cover classification. From this study, the MODIS C6 1-km LST product is recommended for applications.


Sign in / Sign up

Export Citation Format

Share Document