Search

Vegetation phenology and its spatiotemporal driving factors are essential to reflect global climate change, the surface carbon cycle and regional ecology, and further quantitative studies on spatiotemporal heterogeneity and its two-way driving are needed. Based on MODIS phenology, meteorology, land cover and other data from 2001 to 2019, this paper analyzes the phenology change characteristics of the Yangtze River Delta from three dimensions: time, plane space and elevation. Then, the spatiotemporal heterogeneity of phenology and its driving factors are explored with random forest and geographic detector methods. The results show that (1) the advance of start of season (SOS) is insignificant—with 0.17 days per year; the end of season (EOS) shows a significant delay—0.48 days per year. The preseason temperature has a greater contribution to SOS, while preseason precipitation is main factor in determining EOS. (2) Spatial differences of the phenological index do not strictly obey the change rules of latitude at a provincial scale. The SOS of Jiangsu and Anhui is earlier than that of Zhejiang and Shanghai, and EOS shows an obvious double-clustering phenomenon. In addition, a divergent response of EOS with elevation grades is found; the most significant changes are observed at grades below 100 m. (3) Land cover (LC) type is a major factor of the spatial heterogeneity of phenology, and its change may also be one of the insignificant factors driving the interannual change of phenology. Furthermore, nighttime land surface temperature (NLST) has a relatively larger contribution to the spatial heterogeneity in non-core urban areas, but population density (PD) contributes little. These findings could provide a new perspective on phenology and its complex interactions between natural or anthropogenic factors. [ABSTRACT FROM AUTHOR]

As one of the best indicators of the periglacial environment, ice-wedge polygons (IWPs) are important for arctic landscapes, hydrology, engineering, and ecosystems. Thus, a better understanding of the spatiotemporal dynamics and evolution of IWPs is key to evaluating the hydrothermal state and carbon budgets of the arctic permafrost environment. In this paper, the dynamics of ground surface deformation (GSD) in IWP zones (2018–2019) and their influencing factors over the last 20 years in Saskylakh, northwestern Yakutia, Russia were investigated using the Interferometric Synthetic Aperture Radar (InSAR) and Google Earth Engine (GEE). The results show an annual ground surface deformation rate (AGSDR) in Saskylakh at −49.73 to 45.97 mm/a during the period from 1 June 2018 to 3 May 2019. All the selected GSD regions indicate that the relationship between GSD and land surface temperature (LST) is positive (upheaving) for regions with larger AGSDR, and negative (subsidence) for regions with lower AGSDR. The most drastic deformation was observed at the Aeroport regions with GSDs rates of −37.06 mm/a at tower and 35.45 mm/a at runway. The GSDs are negatively correlated with the LST of most low-centered polygons (LCPs) and high-centered polygons (HCPs). Specifically, the higher the vegetation cover, the higher the LST and the thicker the active layer. An evident permafrost degradation has been observed in Saskylakh as reflected in higher ground temperatures, lusher vegetation, greater active layer thickness, and fluctuant numbers and areal extents of thermokarst lakes and ponds. [ABSTRACT FROM AUTHOR]

Incoming surface solar irradiance (SSI) is essential for calculating Earth's surface radiation budget and is a key parameter for terrestrial ecological modeling and climate change research. Remote sensing images from geostationary and polar-orbiting satellites provide an opportunity for SSI estimation through directly retrieving atmospheric and land-surface parameters. This paper presents a new scheme for estimating SSI from the visible and infrared channels of geostationary meteorological and polar-orbiting satellite data. Aerosol optical thickness and cloud microphysical parameters were retrieved from Geostationary Operational Environmental Satellite (GOES) system images by interpolating lookup tables of clear and cloudy skies, respectively. SSI was estimated using pre-calculated offline lookup tables with different atmospheric input data of clear and cloudy skies. The lookup tables were created via the comprehensive radiative transfer model, Santa Barbara Discrete Ordinate Radiative Transfer (SBDART), to balance computational efficiency and accuracy. The atmospheric attenuation effects considered in our approach were water vapor absorption and aerosol extinction for clear skies, while cloud parameters were the only atmospheric input for cloudy-sky SSI estimation. The approach was validated using one-year pyranometer measurements from seven stations in the SURFRAD (SURFace RADiation budget network). The results of the comparison for 2012 showed that the estimated SSI agreed with ground measurements with correlation coefficients of 0.94, 0.69, and 0.89 with a bias of 26.4 W/m-2, -5.9 W/m-2, and 14.9 W/m-2 for clear-sky, cloudy-sky, and all-sky conditions, respectively. The overall root mean square error (RMSE) of instantaneous SSI was 80.0 W/m-2 (16.8%), 127.6 W/m-2 (55.1%), and 99.5 W/m-2 (25.5%) for clear-sky, cloudy-sky (overcast sky and partly cloudy sky), and all-sky (clear-sky and cloudy-sky) conditions, respectively. A comparison with other state-of-the-art studies suggests that our proposed method can successfully estimate SSI with a maximum improvement of an RMSE of 24 W/m-2. The clear-sky SSI retrieval was sensitive to aerosol optical thickness, which was largely dependent on the diurnal surface reflectance accuracy. Uncertainty in the pre-defined horizontal visibility for 'clearest sky' will eventually lead to considerable SSI retrieval error. Compared to cloud effective radius, the retrieval error of cloud optical thickness was a primary factor that determined the SSI estimation accuracy for cloudy skies. Our proposed method can be used to estimate SSI for clear and one-layer cloud sky, but is not suitable for multi-layer clouds overlap conditions as a lower-level cloud cannot be detected by the optical sensor when a higher-level cloud has a higher optical thickness. [ABSTRACT FROM AUTHOR]

Satellite-based temperature measurements are an important indicator for global climate change studies over large areas. Records from Moderate Resolution Imaging Spectroradiometer (MODIS), Advanced Very High Resolution Radiometer (AVHRR) and (Advanced) Along Track Scanning Radiometer ((A)ATSR) are providing long-term time series information. Assessing the quality of remote sensing-based temperature measurements provides feedback to the climate modeling community and other users by identifying agreements and discrepancies when compared to temperature records from meteorological stations. This paper presents a comparison of state-of-the-art remote sensing-based land surface temperature data with air temperature measurements from meteorological stations on a pan-arctic scale (north of 60° latitude). Within this study, we compared land surface temperature products from (A)ATSR, MODIS and AVHRR with an in situ air temperature (Tair) database provided by the National Climate Data Center (NCDC). Despite analyzing the whole acquisition time period of each land surface temperature product, we focused on the inter-annual variability comparing land surface temperature (LST) and air temperature for the overlapping time period of the remote sensing data (2000-2005). In addition, land cover information was included in the evaluation approach by using GLC2000. MODIS has been identified as having the highest agreement in comparison to air temperature records. The time series of (A)ATSR is highly variable, whereas inconsistencies in land surface temperature data from AVHRR have been found. [ABSTRACT FROM AUTHOR]

Urban heat island (UHI), referring to higher temperatures in urban extents than its surrounding rural regions, is widely reported in terms of negative effects to both the ecological environment and human health. To propose effective mitigation measurements, spatiotemporal variations and control machines of surface UHI (SUHI) have been widely investigated, in particular based on the indicator of SUHI intensity (SUHII). However, studies on SUHI frequency (SUHIF), an important temporal indicator, are challenged by a large number of missing data in daily land surface temperature (LST). Whether there is any city with strong SUHII and low SUHIF remains unclear. Thanks to the publication of daily seamless all-weather LST, this paper is proposed to investigate spatiotemporal variations of SUHIF, to compare SUHII and SUHIF, to conduct a pattern classification, and to further explore their driving factors across 305 Chinese cities. Four main findings are summarized below: (1) SUHIF is found to be higher in the south during the day, while it is higher in the north at night. Cities within the latitude from 20° N and 40° N indicate strong intensity and high frequency at day. Climate zone-based variations of SUHII and SUHIF are different, in particular at nighttime. (2) SUHIF are observed in great diurnal and seasonal variations. Summer daytime with 3.01 K of SUHII and 80 of SUHIF, possibly coupling with heat waves, increases the risk of heat-related diseases. (3) K-means clustering is employed to conduct pattern classification of the selected cities. SUHIF is found possibly to be consistent to its SUHII in the same city, while they provide quantitative and temporal characters respectively. (4) Controls for SUHIF and SUHII are found in significant variations among temporal scales and different patterns. This paper first conducts a comparison between SUHII and SUHIF, and provides pattern classification for further research and practice on mitigation measurements. [ABSTRACT FROM AUTHOR]

Remote sensing technologies are critical for analyzing the escalating impacts of global climate change and increasing urbanization, providing vital insights into land surface temperature (LST), land use and cover (LULC) changes, and the identification of urban heat island (UHI) and surface urban heat island (SUHI) phenomena. This research focuses on the nexus between LULC alterations and variations in LST and air temperature (Tair), with a specific emphasis on the intensified SUHI effect in Kharkiv, Ukraine. Employing an integrated approach, this study analyzes time-series data from Landsat and MODIS satellites, alongside Tair climate records, utilizing machine learning techniques and linear regression analysis. Key findings indicate a statistically significant upward trend in Tair and LST during the summer months from 1984 to 2023, with a notable positive correlation between Tair and LST across both datasets. MODIS data exhibit a stronger correlation (R2 = 0.879) compared to Landsat (R2 = 0.663). The application of a supervised classification through Random Forest algorithms and vegetation indices on LULC data reveals significant alterations: a 70.3% increase in urban land and a decrement in vegetative cover comprising a 15.5% reduction in dense vegetation and a 62.9% decrease in sparse vegetation. Change detection analysis elucidates a 24.6% conversion of sparse vegetation into urban land, underscoring a pronounced trajectory towards urbanization. Temporal and seasonal LST variations across different LULC classes were analyzed using kernel density estimation (KDE) and boxplot analysis. Urban areas and sparse vegetation had the smallest average LST fluctuations, at 2.09 °C and 2.16 °C, respectively, but recorded the most extreme LST values. Water and dense vegetation classes exhibited slightly larger fluctuations of 2.30 °C and 2.24 °C, with the bare land class showing the highest fluctuation 2.46 °C, but fewer extremes. Quantitative analysis with the application of Kolmogorov-Smirnov tests across various LULC classes substantiated the normality of LST distributions p > 0.05 for both monthly and annual datasets. Conversely, the Shapiro-Wilk test validated the normal distribution hypothesis exclusively for monthly data, indicating deviations from normality in the annual data. Thresholded LST classifies urban and bare lands as the warmest classes at 39.51 °C and 38.20 °C, respectively, and classifies water at 35.96 °C, dense vegetation at 35.52 °C, and sparse vegetation 37.71 °C as the coldest, which is a trend that is consistent annually and monthly. The analysis of SUHI effects demonstrates an increasing trend in UHI intensity, with statistical trends indicating a growth in average SUHI values over time. This comprehensive study underscores the critical role of remote sensing in understanding and addressing the impacts of climate change and urbanization on local and global climates, emphasizing the need for sustainable urban planning and green infrastructure to mitigate UHI effects. [ABSTRACT FROM AUTHOR]

This is an update of Sobrino et al.'s paper, published in January 2020, which extends the calculation of the Earth's surface temperature to the period 2003–2019 and uses the new version 2019.0 for the sea surface temperature product MODIS, which is available from 15 January 2020 and replaces version 2014.0. The land surface temperature was estimated from the MCD11C1 product for the same period. The results corroborate the temperature anomalies retrieved from climate models and improve the comparison with global annual air temperatures estimated by the NOAA's National Climatic Data Center (NOAA-NCDC), with a correlation coefficient of 0.96. In addition, a trend of 0.021 ± 0.001 °C/year increase was found for the Earth's surface temperature in this 17-year period. [ABSTRACT FROM AUTHOR]

The study of the response of vegetation phenology in the Qinghai Tibet Plateau to various climatic variables is paramount to unveiling the reaction of alpine ecosystems to worldwide climate alterations. Nonetheless, the lagged and cumulative effects of various climatic variables on vegetation phenology in the Qinghai Tibet Plateau remain unclear. Therefore, based on MODIS NDVI data, we extracted vegetation phenological parameters from 2001 to 2020, including the start of the vegetation growing season (SOS) and the end of the vegetation growing season (EOS), and then analyzed the response mechanisms of vegetation phenology to pre-seasonal air temperature (T), precipitation (P), and daytime and nighttime land surface temperatures (DLST, NLST) in the Qinghai Tibet Plateau on the basis of an investigation of the lag and cumulative effects. The results showed that: (1) the multiyear mean values of the SOS mainly occurred from 120 to 160 days, accounting for 86.17% of the study area, while the multiyear mean values of the EOS were mainly concentrated between 260 and 280 days, accounting for 77.05% of the study area; (2) air temperature (T), precipitation (P), and daytime and nighttime land surface temperatures (DLST, NLST) had different degrees of lagging effects on the SOS and the EOS. Among them, the time lag effect of precipitation on vegetation phenology was more pronounced; (3) different climatic variables had distinct cumulative effects on vegetation phenology. In contrast to the insignificant cumulative effects of temperature and nighttime surface temperature on the SOS and the EOS, the cumulative effects of precipitation and daytime land surface temperature on the SOS were more pronounced than those on the EOS; (4) the SOS and air temperature, precipitation, and NLST were mainly negatively correlated, in which the proportion of the negative correlation between SOS and NLST was up to 68.80%, and SOS and DLST were mainly positively correlated with a positive correlation proportion of 73.27%, EOS and air temperature, precipitation, and NLST were positively correlated with a positive correlation proportion of EOS and precipitation of up to 71.52%, and EOS and DLST were mainly negatively correlated with a negative correlation ratio of 55.87%. [ABSTRACT FROM AUTHOR]

The mountain systems of the Himalayan regions are changing rapidly due to climatic change at a local and global scale. The Indian Western Himalaya ecosystem (between the tree line and the snow line) is an underappreciated component. Yet, knowledge of vegetation distribution, rates of change, and vegetation interactions with snow-hydroclimatic elements is lacking. The purpose of this study is to investigate the linkage between the spatiotemporal variability of vegetation (i.e., greenness and forest) and related snow-hydroclimatic parameters (i.e., snow cover, land surface temperature, Tropical Rainfall Measuring Mission (TRMM), and Evapotranspiration (ET)) in Himachal Pradesh (HP) Basins (i.e., Beas, Chandra, and Bhaga). Spatiotemporal variability in forest and grassland has been estimated from MODIS land cover product (MCD12Q1) using Google Earth Engine (GEE) for the last 19 years (2001–2019). A significant inter- and intra-annual variation in the forest, grassland, and snow-hydroclimatic factors have been observed during the data period in HP basins (i.e., Beas, Chandra, and Bhaga basin). The analysis demonstrates a significant decrease in the forest cover (214 ha/yr.) at the Beas basin; however, a significant increase in grassland cover is noted at the Beas basin (459 ha/yr.), Chandra (176.9 ha/yr.), and Bhaga basin (9.1 ha/yr.) during the data period. Spatiotemporal forest cover loss and gain in the Beas basin have been observed at ~7504 ha (6.6%) and 1819 ha (1.6%), respectively, from 2001 to 2019. However, loss and gain in grassland cover were observed in 3297 ha (2.9%) and 10,688 ha (9.4%) in the Beas basin, 1453 ha (0.59%) and 3941 ha (1.6%) in the Chandra basin, and 1185 ha (0.92%) and 773 ha (0.60%) in the Bhaga basin, respectively. Further, a strong negative correlation (r = −0.65) has been observed between forest cover and evapotranspiration (ET). However, a strong positive correlation (r = 0.99) has been recorded between grassland cover and ET as compared to other factors. The main outcome of this study in terms of spatiotemporal loss and gain in forest and grassland shows that in the Bhaga basin, very little gain and loss have been observed as compared to the Chandra and Beas basins. The present study findings may provide important aid in the protection and advancement of the knowledge gap of the natural environment and the management of water resources in the HP Basin and other high-mountain regions of the Himalayas. For the first time, this study provides a thorough examination of the spatiotemporal variability of forest and grassland and their interactions with snow-hydroclimatic factors using GEE for Western Himalaya. [ABSTRACT FROM AUTHOR]

Accurate, seamless, and long-term land surface temperature (LST) data sets are crucial for investigating climate change and agriculture production. However, factors like cloud contamination have led to invalid values in the LST product, which has restricted the application of the LST dataset. Therefore, the reconstruction of LST products is challenging, and it is attracting widespread attention. This study compared the performance of different algorithms (XGBoost, GBDT, RF, POLY, MLR) and different training sets (using only good-quality pixels or using both good-quality and other-quality pixels) in the estimation of missing pixels in the LST data, obtaining a seamless daily 1 km LST dataset of MODIS Terra-day, Aqua-day, Terra-night, and Aqua-night data for Zhejiang Province and its surrounding areas from 2000 to 2022. The results demonstrated that the performance of machine-learning models is significantly better than that of linear models, and among the five models, XGBoost performed the best, with an RMSE of less than 1 °C. The Wilcoxon test between the reconstructed LST and the true LST values revealed that including both good-quality and other-quality pixels for reconstruction resulted in a 33% increase in the number of days with non-significant differences compared with using only good-quality pixels. Moreover, the reconstructed nighttime LST has a lower RMSE compared with the reconstructed daytime LST, and the RMSE of the reconstructed LST on the Terra satellite is lower than the RMSE of the reconstructed LST on the Aqua satellite. The RMSEs for the reconstructed LSTs are 0.50 °C, 0.61 °C, 0.36 °C, and 0.39 °C, corresponding to Terra-day, Aqua-day, Terra-night, and Aqua-night for images with coverage reaching 70%, 0.65 °C, 0.83 °C, 0.49 °C, respectively, and 0.52 °C for images with coverage less than 70%. The accuracy of the reconstructed LSTs using our proposed framework outperforms the existing reconstruction methods. The 1 km daily seamless LST products can be applied in various fields, such as air temperature estimation, climate change, urban heat island, and crop temperature stress monitoring. [ABSTRACT FROM AUTHOR]

The assessment of satellite-derived land surface temperature (LST) data is essential to ensure their high quality for climate applications and research. This study intercompared seven LST products (i.e., ATSR_3, MODISA, MODIST, SLSTRA, SLSTRB, SEVIR2 and SEVIR4) of the European Space Agency's (ESA) LST Climate Change Initiative (LST_cci) project, which are retrieved for polar and geostationary orbit satellites, and three operational LST products: NASA's MODIS MOD11/MYD11 LST and ESA's AATSR LST. All data were re-gridded on to a common spatial grid of 0.05° and matched for concurrent overpasses within 5 min. The matched data were analysed over Europe and Africa for monthly and seasonally aggregated median differences and studied for their dependence on land cover class and satellite viewing geometry. For most of the data sets, the results showed an overall agreement within ±2 K for median differences and robust standard deviation (RSD). A seasonal variation of median differences between polar and geostationary orbit sensor data was observed over Europe, which showed higher differences in summer and lower in winter. Over all land cover classes, NASA's operational MODIS LST products were about 2 K colder than the LST_cci data sets. No seasonal differences were observed for the different land covers, but larger median differences between data sets were seen over bare soil land cover classes. Regarding the viewing geometry, an asymmetric increase of differences with respect to nadir view was observed for day-time data, which is mainly caused by shadow effects. For night-time data, these differences were symmetric and considerably smaller. Overall, despite the differences in the LST retrieval algorithms of the intercompared data sets, a good consistency between the LST_cci data sets was determined. [ABSTRACT FROM AUTHOR]

A good understanding of the processes of land surface temperature (LST) change is important for assessing regional climate change. In the present study, we obtained the MODIS MOD11A2 LST products over the Yangtze River Delta (YRD) from 2001 to 2020. In order to comprehensively assess the spatial and temporal variability of LST in the YRD region over the past two decades, the Theil–Sen Median trend analysis and Mann–Kendall test, BFAST01 trend decomposition, and landscape pattern analysis were used in this study. We show that the rate of linear change in LST in the YRD ranges from −0.019 °C/month to 0.046 °C/month. The BFAST01 trend decomposition identifies more details of LST change and monotonic increases, reversal increase, and interruption increase are the main warming trends. The distribution of the different trend types shows strong aggregation with high spatial heterogeneity. The LST breakpoints are mainly located in the northern and southern YRD, which frequently occurred during 2010–2013. Of the various land types, breakpoints occur most frequently in cropland and high NDVI (0.5–0.7) areas, and the intensity of most of them is within 2 °C. In addition, much stronger warming occurs in urban areas than in other land types. Our study provides a better understanding of the dynamics of LST in the YRD region over the past 20 years and highlights that breakpoints cannot be circumvented in regional temperature assessment. [ABSTRACT FROM AUTHOR]

Satellite-derived Land Surface Temperature (LST) dynamics have been increasingly used to study various geophysical processes. This review provides an extensive overview of the applications of LST in the context of global change. By filtering a selection of relevant keywords, a total of 164 articles from 14 international journals published during the last two decades were analyzed based on study location, research topic, applied sensor, spatio-temporal resolution and scale and employed analysis methods. It was revealed that China and the USA were the most studied countries and those that had the most first author affiliations. The most prominent research topic was the Surface Urban Heat Island (SUHI), while the research topics related to climate change were underrepresented. MODIS was by far the most used sensor system, followed by Landsat. A relatively small number of studies analyzed LST dynamics on a global or continental scale. The extensive use of MODIS highly determined the study periods: A majority of the studies started around the year 2000 and thus had a study period shorter than 25 years. The following suggestions were made to increase the utilization of LST time series in climate research: The prolongation of the time series by, e.g., using AVHRR LST, the better representation of LST under clouds, the comparison of LST to traditional climate change measures, such as air temperature and reanalysis variables, and the extension of the validation to heterogenous sites. [ABSTRACT FROM AUTHOR]

As an essential climate variable (ECV), land surface albedo plays an important role in the Earth surface radiation budget and regional or global climate change. The Tibetan Plateau (TP) is a sensitive environment to climate change, and understanding its albedo seasonal and inter-annual variations is thus important to help capture the climate change rules. In this paper, we analyzed the large-scale spatial patterns, temporal trends, and seasonal variability of land surface albedo overall the TP, based on the moderate resolution imaging spectroradiometer (MODIS) MCD43 albedo products from 2001 to 2019. Specifically, we assessed the correlations between the albedo anomaly and the anomalies of normalized difference vegetation index (NDVI), the fraction of snow cover (snow cover), and land surface temperature (LST). The results show that there are larger albedo variations distributed in the mountainous terrain of the TP. Approximately 10.06% of the land surface is identified to have been influenced by the significant albedo variation from the year 2001 to 2019. The yearly averaged albedo was decreased significantly at a rate of 0.0007 (Sen's slope) over the TP. Additionally, the yearly average snow cover was decreased at a rate of 0.0756. However, the yearly average NDVI and LST were increased with slopes of 0.0004 and 0.0253 over the TP, respectively. The relative radiative forcing (RRF) caused by the land cover change (LCC) is larger than that caused by gradual albedo variation in steady land cover types. Overall, the RRF due to gradual albedo variation varied from 0.0005 to 0.0170 W/m2, and the RRF due to LCC variation varied from 0.0037 to 0.0243 W/m2 during the years 2001 to 2019. The positive RRF caused by gradual albedo variation or the LCC can strengthen the warming effects in the TP. The impact of the gradual albedo variations occurring in the steady land cover types was very low between 2001 and 2019 because the time series was short, and it therefore cannot be neglected when examining radiative forcing for a long time series regarding climate change. [ABSTRACT FROM AUTHOR]

In the Tibetan Plateau (TP), the changes of lake ice phenology not only reflect regional climate change, but also impose substantial ecohydrological impacts on the local environment. Due to the limitation of ground observation, remote sensing has been used as an alternative tool to investigate recent changes of lake ice phenology. However, uncertainties exist in the remotely sensed lake ice phenology owing to both the data and methods used. In this paper, three different remotely sensed datasets are used to investigate the lake ice phenology variation in the past decade across the Tibetan Plateau, with the consideration of the underlying uncertainties. The remotely sensed data used include reflectance data, snow product, and land surface temperature (LST) data of moderate resolution imaging spectroradiometer (MODIS). The uncertainties of the three methods based on the corresponding data are assessed using the triple collocation approach. Comparatively, it is found that the method based on reflectance data outperforms the other two methods. The three methods are more consistent in determining the thawing dates rather than the freezing dates of lake ice. It is consistently shown by the three methods that the ice-covering duration in the northern part of the TP lasts longer than that in the south. Though there is no general trend of lake ice phenology across the TP for the period of 2000–2015, the warmer climate and stronger wind have led to the earlier break-up of lake ice. [ABSTRACT FROM AUTHOR]

In this work, we investigate how the seasonal hysteresis of the Surface Urban Heat Island Intensity (SUHII) differs across climates and provide a detailed typology of the daytime and nighttime SUHII hysteresis loops. Instead of the typical tropical/dry/temperate/continental grouping, we describe Earth's climate using the Köppen–Geiger system that empirically maps Earth's biome distribution into 30 climate classes. Our thesis is that aggregating multi-city data without considering the biome of each city results in temporal means that fail to reflect the actual SUHII characteristics. This is because the SUHII is a function of both urban and rural features and the phenology of the rural surroundings can differ considerably between cities, even in the same climate zone. Our investigation covers all the densely populated areas of Earth and uses 18 years (2000–2018) of land surface temperature and land cover data from the European Space Agency's Climate Change Initiative. Our findings show that, in addition to concave-up and -down shapes, the seasonal hysteresis of the SUHII also exhibits twisted, flat, and triangle-like patterns. They also suggest that, in wet climates, the daytime SUHII hysteresis is almost universally concave-up, but they paint a more complex picture for cities in dry climates. [ABSTRACT FROM AUTHOR]

Surface shortwave (SW) irradiation is the primary driving force of energy exchange in the atmosphere and land interface. The global climate is profoundly influenced by irradiation changes due to the special climatic condition in Antarctica. Remote-sensing retrieval can offer only the instantaneous values in an area, whilst daily cycle and average values are necessary for further studies and applications, including climate change, ecology, and land surface process. When considering the large values of and small diurnal changes of solar zenith angle and cloud coverage, we develop two methods for the temporal extension of remotely sensed downward SW irradiance over Antarctica. The first one is an improved sinusoidal method, and the second one is an interpolation method based on cloud fraction change. The instantaneous irradiance data and cloud products are used in both methods to extend the diurnal cycle, and obtain the daily average value. Data from South Pole and Georg von Neumayer stations are used to validate the estimated value. The coefficient of determination (R2) between the estimated daily averages and the measured values based on the first method is 0.93, and the root mean square error (RMSE) is 32.21 W/m2 (8.52%). As for the traditional sinusoidal method, the R2 and RMSE are 0.68 and 70.32 W/m2 (18.59%), respectively The R2 and RMSE of the second method are 0.96 and 25.27 W/m2 (6.98%), respectively. These values are better than those of the traditional linear interpolation (0.79 and 57.40 W/m2 (15.87%)). [ABSTRACT FROM AUTHOR]

Urban growth and its associated expansion of built-up areas are expected to continue through to the twenty second century and at a faster pace in developing countries. This has the potential to increase thermal discomfort and heat-related distress. There is thus a need to monitor growth patterns, especially in resource constrained countries such as Africa, where few studies have so far been conducted. In view of this, this study compares urban growth and temperature response patterns in Freetown and Bo town in Sierra Leone. Multispectral Landsat images obtained in 1998, 2000, 2007, and 2015 are used to quantify growth and land surface temperature responses. The contribution index (CI) is used to explain how changes per land use and land cover class (LULC) contributed to average city surface temperatures. The population size of Freetown was about eight times greater than in Bo town. Landsat data mapped urban growth patterns with a high accuracy (Overall Accuracy > 80%) for both cities. Significant changes in LULC were noted in Freetown, characterized by a 114 km2 decrease in agriculture area, 23 km2 increase in dense vegetation, and 77 km2 increase in built-up area. Between 1998 and 2015, built-up area increased by 16 km2, while dense vegetation area decreased by 14 km2 in Bo town. Average surface temperature increased from 23.7 to 25.5 °C in Freetown and from 24.9 to 28.2 °C in Bo town during the same period. Despite the larger population size and greater built-up extent, as well as expansion rate, Freetown was 2 °C cooler than Bo town in all periods. The low temperatures are attributed to proximity to sea and the very large proportion of vegetation surrounding the city. Even close to the sea and abundant vegetation, the built-up area had an elevated temperature compared to the surroundings. The findings are important for formulating heat mitigation strategies for both inland and coastal cities in developing countries. [ABSTRACT FROM AUTHOR]

Global change affects vegetation cover and processes through multiple pathways. Long time series of surface land surface properties derived from satellite remote sensing give unique abilities to observe these changes, particularly in areas with complex topography and limited research infrastructure. Here, we focus on Nepal, a biodiversity hotspot where vegetation productivity is limited by moisture availability (dominated by a summer monsoon) at lower elevations and by temperature at high elevations. We analyze the normalized difference vegetation index (NDVI) from 1981 to 2015 semimonthly, at an 8 km spatial resolution. We use a random forest (RF) of regression trees to generate a statistical model of the NDVI as a function of elevation, land use, CO2 level, temperature, and precipitation. We find that the NDVI increased over the studied period, particularly at low and middle elevations and during the fall (post-monsoon). We infer from the fitted RF model that the NDVI linear trend is primarily due to CO2 level (or another environmental parameter that is changing quasi-linearly), and not primarily due to temperature or precipitation trends. On the other hand, interannual fluctuation in the NDVI is more correlated with temperature and precipitation. The RF accurately fits the available data and shows promise for estimating trends and testing hypotheses about their causes. [ABSTRACT FROM AUTHOR]

Urban parks have been shown to form park cool islands (PCIs), which can effectively alleviate the negative influences of urban heat islands (UHI). However, few studies have examined the detailed characteristics of PCIs, the effect of urban park features on their individual temperatures, and monthly variation in PCIs. Land surface temperature (LST) retrieved from Landsat 8 TIR images between May and October were used to represent the thermal environment. Urban park characteristics were extracted from high-resolution GF-2 images. Using these datasets, the relationships between urban park characteristics and PCIs were explored in this study using Changchun, which has a snow climate, as a case study. The results showed the following: (1) the urban parks exhibited a cooling island effect, and the PCIs showed significant monthly variations with the highest intensities in the hot months; (2) the effects of composition (e.g., park size and the percentage of water area) on LSTs and PCIs showed significant monthly variability and were stronger than the configuration effects. Furthermore, an unexpected, negative correlation between PCIs and the area of park grass was also found; and (3) larger parks tended to have stronger PCI intensities and extents of influence. For parks larger than 30 ha, the cooling effects extended approximately 480 m from the park edge between June and August. For all of parks during the study duration, the rate of temperature increase was highest within 60 m from the park edge. The PCI we employ specifically in this study is characterized by LST. [ABSTRACT FROM AUTHOR]

Global climate change and increasing urbanization worldwide intensify the need for a better understanding of human heat stress dynamics in urban systems. During heat waves, which are expected to increase in number and intensity, the development of urban cool islands could be a lifesaver for many elderly and vulnerable people. The use of remote sensing data offers the unique possibility to study these dynamics with spatially distributed large datasets during all seasons of the year and including day and night-time analysis. For the city of Basel 32 high-quality Landsat 8 (L8) scenes are available since 2013, enabling comprehensive statistical analysis. Therefore, land surface temperature (LST) is calculated using L8 thermal infrared (TIR) imagery (stray light corrected) applying improved emissivity and atmospheric corrections. The data are combined with a land use/land cover (LULC) map and evaluated using administrative residential units. The observed dependence of LST on LULC is analyzed using a thermal unmixing approach based on a multiple linear regression (MLR) model, which allows for quantifying the gradual influence of different LULC types on the LST precisely. Seasonal variations due to different solar irradiance and vegetation cover indicate a higher dependence of LST on the LULC during the warmer summer months and an increasing influence of the topography and albedo during the colder seasons. Furthermore, the MLR analysis allows creating predicted LST images, which can be used to fill data gaps like in SLC-off Landsat 7 ETM+ data. [ABSTRACT FROM AUTHOR]

Air temperature (Tair) near the ground surface is a fundamental descriptor of terrestrial environment conditions and one of the most widely used climatic variables in global change studies. The main objective of this study was to explore the possibility of retrieving high-resolution Tair from the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) products, covering complex terrain in Northeast China. The All Subsets Regression (ASR) method was adopted to select the predictors and build optimal multiple linear regression models for estimating maximum (Tmax), minimum (Tmin), and mean (Tmean) air temperatures. The relative importance of predictors in these models was evaluated via the Standardized Regression Coefficients (SRCs) method. The results indicated that the optimal models could estimate the Tmax, Tmin, and Tmean with relatively high accuracies (Model Efficiency ≥ 0.90). Both LST and day length (DL) predictors were important in estimating Tmax (SRCs: daytime LST = 0.53, DL = 0.35), Tmin (SRCs: nighttime LST = 0.74, DL = 0.23), and Tmean (SRCs: nighttime LST = 0.72, DL = 0.28). Models predicting Tmin and Tmean had better performance than the one predicting Tmax. Nighttime LST was better at predicting Tmin and Tmean than daytime LST data at predicting Tmax. Land covers had noticeable influences on estimating Tair, and even seasonal vegetation greening could result in temporal variations of model performance. Air temperature could be accurately estimated using remote sensing, but the model performance was varied across different spatial and temporal scales. More predictors should be incorporated for the purpose of improving the estimation of near surface Tair from the MODIS LST production. [ABSTRACT FROM AUTHOR]

Satellite remote sensing provides a rapid and broad-scale means for monitoring vegetation phenology and its relationship with fluctuations in air temperature. Investigating the response of plant communities to climate change is needed to gain insight into the potentially detrimental effects on ecosystem processes. While many studies have used satellite-derived land surface temperature (LST) as a proxy for air temperature, few studies have attempted to create and validate models of forest understory temperature (Tust), as it is obscured from these space-borne observations. This study worked to predict instantaneous values of Tust using daily Moderate Resolution Imaging Spectroradiometer (MODIS) LST data over a 99,000 km2 study area located in the Rocky Mountains of western Alberta, Canada. Specifically, we aimed to identify the forest characteristics that improve estimates of Tust over using LST alone. Our top model predicted Tust to within a mean absolute error (MAE) of 1.4 °C with an overall model fit of R2 = 0.89 over two growing seasons. Canopy closure and the LiDAR-derived standard deviation of canopy height metric were found to significantly improve estimations of Tust over MODIS LST alone. These findings demonstrate that canopy structure and forest stand-type function to differentiate understory air temperatures from ambient canopy temperature as seen by the sensor overhead. [ABSTRACT FROM AUTHOR]

Under the circumstance of global climate change, the evolution of thermal environments has attracted more attention, for which the surface urban heat island (SUHI) is one of the major concerns. In this research, we focused on the spatiotemporal patterns for two "stove cities" in China, i.e., Wuhan and Nanchang, based on the long-term (1984–2018) and fine-scale (Landsat-like) series of satellite images. The results showed opposite spatiotemporal patterns for the two cities, even though they were both widely concerned to be the hottest cities. No matter which definition of surface urban heat island intensity (SUHII) was selected, Nanchang presented higher and more fluctuating SUHII than Wuhan, with a relatively higher land surface temperature (LST) of the urban area and lower LST of the rural area in Nanchang, especially in recent years. For the spatial pattern, the highest LST center (i.e., the SUHI) has expanded obviously for the past 35 years in Nanchang. For Wuhan, the LST in SUHI has gone through a trend of a relatively increase at first, followed by a decrease. For the temporal pattern, an increasing trend of LST could be detected in Nanchang. However, the LST in Wuhan presented a slightly decreasing trend. Moreover, the SUHII evolution in Nanchang decreased at first then increased, while Wuhan showed a slight increasing trend at first, followed by a decrease for SUHII. In addition, different SUHII definitions would not affect the spatial pattern and temporal trend of SUHI, but only controlled the exact SUHII value, especially in those years with extreme weather. [ABSTRACT FROM AUTHOR]

In the Alps, understanding how climate change is affecting evapotranspiration (ET) is relevant due to possible implications on water availability for large lowland areas of Europe. Here, changes in ET were studied based on 20 years of MODIS data. MOD16 and operational Simplified Surface Energy Balance (SSEBop) products were compared with eddy-covariance data and analyzed for trend detection. The two products showed a similar relationship with ground observations, with RMSE between 0.69 and 2 mm day−1, and a correlation coefficient between 0.6 and 0.83. A regression with the potential drivers of ET showed that, for climate variables, ground data were coherent with MOD16 at grassland sites, where r2 was 0.12 for potential ET, 0.17 for precipitation, and 0.57 for air temperature, whereas ground data agreed with SSEBop at forest sites, with an r2 of 0.46 for precipitation, no correlation with temperature, and negative correlation with potential ET. Interestingly, ground-based correlation corresponded to SSEBop for leaf area index (LAI), while it matched with MOD16 for land surface temperature (LST). Through the trend analysis, both MOD16 and SSEBop revealed positive trends in the south-west, and negative trends in the south and north-east. Moreover, in summer, positive trends prevailed at high elevations for grasslands and forests, while negative trends dominated at low elevations for croplands and grasslands. However, the Alpine area share with positive ET trends was 16.6% for MOD16 and 3.9% for SSEBop, while the share with negative trends was 1.2% for MOD16 and 15.3% for SSEBop. A regression between trends in ET and in climate variables, LST, and LAI indicated consistency, especially between ET, temperature, and LAI increase, but low correlation. Overall, the discrepancies in the trends, and the fact that none of the two products outperformed the other when compared to ground data, suggest that, in the Alps, SSEBop and MOD16 might not be accurate enough to be a robust basis to study ET changes. [ABSTRACT FROM AUTHOR]

Unsustainable development paths have reached critical levels in Europe. In recent years, in cities, urbanization has been contributing to the intensification of urban heat islands. To analyze the development of surface urban heat islands (SUHI) in Europe in the last few years, the present study combines the land surface temperature (LST) from MODIS with the urban classes of the CORINE land cover data within 617 functional urban areas (FUAs). Urban and industrial uses have significantly higher LST than green urban areas across all years (about 4 to 6 °C), as do agricultural areas within cities. Besides land cover, location also influences LST differences. While, e.g., Bolzano (Italy) shows particularly large LST differences (>6 °C) between the core and the commuting zone, this effect is hardly visible in Porto (Portugal) and Madrid (Spain) (<2.5 °C). Cities of moderate climates show increasing differences between a city and its commuting zones with rising LST (r = 0.68), i.e., less cooling effects at night. [ABSTRACT FROM AUTHOR]

The permafrost in the Qilian Mountains (QLMs), the northeastern margin of the Qinghai–Tibet Plateau, changed dramatically in the context of climate warming and increasing anthropogenic activities, which poses significant influences on the stability of the ecosystem, water resources, and greenhouse gas cycles. Yet, the characteristics of the frozen ground in the QLMs are largely unclear regarding the spatial distribution of active layer thickness (ALT), the maximum frozen soil depth (MFSD), and the temperature at the top of the permafrost or the bottom of the MFSD (TTOP). In this study, we simulated the dynamics of the ALT, TTOP, and MFSD in the QLMs in 2004–2019 in the Google Earth Engine (GEE) platform. The widely-adopted Stefan Equation and TTOP model were modified to integrate with the moderate-resolution imaging spectroradiometer (MODIS) land surface temperature (LST) in GEE. The N-factors, the ratio of near-surface air to ground surface freezing and thawing indices, were assigned to the freezing and thawing indices derived with MODIS LST in considerations of the fractional vegetation cover derived from MODIS normalized difference vegetation index (NDVI). The results showed that the GEE platform and remote sensing imagery stored in Google cloud could be quickly and effectively applied to obtain the spatial and temporal variation of permafrost distribution. The area with TTOP < 0 °C is 8.4 × 104 km2 (excluding glaciers and lakes) and accounts for 46.6% of the whole QLMs, the regional mean ALT is 2.43 ± 0.44 m, while the regional mean MFSD is 2.54 ± 0.45 m. The TTOP and ALT increase with the decrease of elevation from the sources of the sub-watersheds to middle and lower reaches. There is a strong correlation between TTOP and elevation (slope = −1.76 °C km−1, p < 0.001). During 2004–2019, the area of permafrost decreased by 20% at an average rate of 0.074 × 104 km2·yr−1. The regional mean MFSD decreased by 0.1 m at a rate of 0.63 cm·yr−1, while the regional mean ALT showed an exception of a decreasing trend from 2.61 ± 0.45 m during 2004–2005 to 2.49 ± 0.4 m during 2011–2015. Permafrost loss in the QLMs in 2004–2019 was accelerated in comparison with that in the past several decades. Compared with published permafrost maps, this study shows better calculation results of frozen ground in the QLMs. [ABSTRACT FROM AUTHOR]

Land surface temperature (LST) plays a critical role in the water cycle and energy balance at global and regional scales. Large-scale LST estimates can be obtained from satellite observations and reanalysis data. In this study, we first investigate the long-term changes of LST during 2003–2017 on a per-pixel basis using three different datasets derived from (i) the Atmospheric Infrared Sounder (AIRS) onboard Aqua satellite, (ii) the Moderate Resolution Imaging Spectroradiometer (MODIS) also aboard Aqua, and (iii) the recently released ERA5-Land reanalysis data. It was found that the spatio-temporal patterns of these data agree very well. All three products globally showed an uptrend in the annual average LST during 2003–2017 but with considerable spatial variations. The strongest increase was found over the region north of 45° N, particularly over Asian Russia, whereas a slight decrease was observed over Australia. The regression analysis indicated that precipitation (P), incoming surface solar radiation (SW↓), and incoming surface longwave radiation (LW↓) can together explain the inter-annual LST variations over most regions, except over tropical forests, where the inter-annual LST variation is low. Spatially, the LST changes during 2003–2017 over the region north of 45° N were mainly influenced by LW↓, while P and SW↓ played a more important role over other regions. A detailed look at Asian Russia and the Amazon rainforest at a monthly time scale showed that warming in Asian Russia is dominated by LST increases in February–April, which are closely related with the simultaneously increasing LW↓ and clouds. Over the southern Amazon, the most apparent LST increase is found in the dry season (August–September), primarily affected by decreasing P. In addition, increasing SW↓ associated with decreasing atmospheric aerosols was another factor found to cause LST increases. This study shows a high level of consistency among LST trends derived from satellite and reanalysis products, thus providing more robust characteristics of the spatio-temporal LST changes during 2003–2017. Furthermore, the major climatic drivers of LST changes during 2003–2017 were identified over different regions, which might help us predict the LST in response to changing climate in the future. [ABSTRACT FROM AUTHOR]

Lake phenology is essential for understanding the lake freeze-thaw cycle effects on terrestrial hydrological processes. The Qinghai-Tibetan Plateau (QTP) has the most extensive ice reserve outside of the Arctic and Antarctic poles and is a sensitive indicator of global climate changes. Qinghai Lake, the largest lake in the QTP, plays a critical role in climate change. The freeze-thaw cycles of lakes were studied using daily Moderate Resolution Imaging Spectroradiometer (MODIS) data ranging from 2000–2018 in the Google Earth Engine (GEE) platform. Surface water/ice area, coverage, critical dates, surface water, and ice cover duration were extracted. Random forest (RF) was applied with a classifier accuracy of 0.9965 and a validation accuracy of 0.8072. Compared with six common water indexes (tasseled cap wetness (TCW), normalized difference water index (NDWI), modified normalized difference water index (MNDWI), automated water extraction index (AWEI), water index 2015 (WI2015) and multiband water index (MBWI)) and ice threshold value methods, the critical freeze-up start (FUS), freeze-up end (FUE), break-up start (BUS), and break-up end (BUE) dates were extracted by RF and validated by visual interpretation. The results showed an R2 of 0.99, RMSE of 3.81 days, FUS and BUS overestimations of 2.50 days, and FUE and BUE underestimations of 0.85 days. RF performed well for lake freeze-thaw cycles. From 2000 to 2018, the FUS and FUE dates were delayed by 11.21 and 8.21 days, respectively, and the BUS and BUE dates were 8.59 and 1.26 days early, respectively. Two novel key indicators, namely date of the first negative land surface temperature (DFNLST) and date of the first positive land surface temperature (DFPLST), were proposed to comprehensively delineate lake phenology: DFNLST was approximately 37 days before FUS, and DFPLST was approximately 20 days before BUS, revealing that the first negative and first positive land surface temperatures occur increasingly earlier. [ABSTRACT FROM AUTHOR]

Forest cover change is critical in the regulation of global and regional climate change through the alteration of biophysical features across the Earth's surface. The accurate assessment of forest cover change can improve our understanding of its roles in the regulation processes of surface temperature. In spite of this, few researchers have attempted to discern the varying effects of multiple satellite-derived forest changes on local surface temperatures. In this study, we quantified the actual contributions of forest loss and gain associated with evapotranspiration (ET) and albedo to local surface temperature in Guangdong Province, China using an improved spatiotemporal change pattern analysis method, and explored the interrelationships between surface temperature and air temperature change. We specifically developed three forest change products for Guangdong, combining satellite observations from Landsat, PALSAR, and MODIS for comparison. Our results revealed that the adjusted simple change detection (SCD)-based Landsat/PALSAR forest cover data performed relatively well. We found that forest loss and gain between 2000 and 2010 had opposite effects on land surface temperature (LST), ET, and albedo. Forest gain led to a cooling of −0.12 ± 0.01 °C, while forest loss led to a warming of 0.07 ± 0.01 °C, which were opposite to the anomalous change of air temperature. A reduced warming to a considerable cooling was estimated due to the forest gain and loss across latitudes. Specifically, mid-subtropical forest gains increased LST by 0.25 ± 0.01 °C, while tropical forest loss decreased LST by −0.16 ± 0.05 °C, which can demonstrate the local differences in an overall cooling. ET induced cooling and warming effects were appropriate for most forest gain and loss. Meanwhile, the nearby temperature changes caused by no-change land cover types more or less canceled out some of the warming and cooling. Albedo exhibited negligible and complex impacts. The other two products (i.e., the GlobeLand30 and MCD12Q1) affect the magnitude of temperature response due to the discrepancies in forest definition, methodology, and data resolution. This study highlights the non-negligible contributions of high-resolution maps and a robust temperature response model in the quantification of the extent to which forest gain reverses the climate effects of forest loss under global warming. [ABSTRACT FROM AUTHOR]

The scientific community has widely reported the impacts of climate change on the Central Himalaya. To qualify and quantify these effects, long-term land surface temperature observations in both the daytime and nighttime, acquired by the Moderate Resolution Imaging Spectroradiometer from 2000 to 2017, were used in this study to investigate the spatiotemporal variations and their changing mechanism. Two periodic parameters, the mean annual surface temperature (MAST) and the annual maximum temperature (MAXT), were derived based on an annual temperature cycle model to reduce the influences from the cloud cover and were used to analyze their trend during the period. The general thermal environment represented by the average MAST indicated a significant spatial distribution pattern along with the elevation gradient. Behind the clear differences in the daytime and nighttime temperatures at different physiographical regions, the trend test conducted with the Mann-Kendall (MK) method showed that most of the areas with significant changes showed an increasing trend, and the nighttime temperatures exhibited a more significant increasing trend than the daytime temperatures, for both the MAST and MAXT, according to the changing areas. The nighttime changing areas were more widely distributed (more than 28%) than the daytime changing areas (around 10%). The average change rates of the MAST and MAXT in the daytime are 0.102 °C/yr and 0.190 °C/yr, and they are generally faster than those in the nighttime (0.048 °C/yr and 0.091 °C/yr, respectively). The driving force analysis suggested that urban expansion, shifts in the courses of lowland rivers, and the retreat of both the snow and glacier cover presented strong effects on the local thermal environment, in addition to the climatic warming effect. Moreover, the strong topographic gradient greatly influenced the change rate and evidenced a significant elevation-dependent warming effect, especially for the nighttime LST. Generally, this study suggested that the nighttime temperature was more sensitive to climate change than the daytime temperature, and this general warming trend clearly observed in the central Himalayan region could have important influences on local geophysical, hydrological, and ecological processes. [ABSTRACT FROM AUTHOR]