Studies on the Relationships Between Land Surface Temperature and Environmental Factors in an Inland River Catchment Based on Geographically Weighted Regression and MODIS Data

Abstract

Despite growing concerns in Land Surface Temperature (LST) and its related environmental factors (geographical, climate, and atmospheric conditions), little attention was about the spatial variation that consider above conditions together. Our purpose is to analyze and quantify LST and related environmental factors, using Geographically Weighted Regression (GWR), and Moderate Resolution Imaging Spectroradiometer (MODIS) data in a typical inland river catchment, named Heihe River catchment, China. Considering thirteen environmental factors (altitude, latitude, Topographic Wetness Index, Cos(aspect), temperature, precipitation, humidity, wind speed, radiation, albedo, the normalized difference vegetation index (NDVI), water vapor, COT), 18 GWR models were set up. Results showed that yearly averaged LST changed from 264 K to 309 K, with the highest value recorded in the downstream desert region. LST has the same variable trend and seasonality with NDVI, precipitable water vapor, and cloud optical thickness (COT), but has an inverse relationship with albedo. All GWR models indicated better simulation with smaller Akaike Information Criterion (AICc), and higher coefficient of determination (R <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> ), compared with Ordinary Least Squares method (OLS). Furthermore, performance of multi-factor analysis was better than single-factor analysis, with model 18 showing the best performance achieving higher R <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> (0.94) and lower AICc (7760). For all GWR model, 86.4% of R <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> was higher than 0.60, most values distributed in the range of 0.80-0.99, and 86.59% of residual values were within the range of ±2 K. Different parameters resulted in different slope distribution, which indicated that altitude is the major driving factor, followed by NDVI, and albedo.