The spatially heterogeneous nature and geographical scale of surface urban heat island (SUHI) driving mechanisms remain largely unknown, as most previous studies have focused solely on their global performance and impact strength. This paper analyzes diurnal and nocturnal SUHIs in China based on the multiscale geographically weighted regression (MGWR) model for 2005, 2010, 2015, and 2018. Compared to results obtained using the ordinary least square (OLS) model, the MGWR model has a lower corrected Akaike information criterion value and significantly improves the model’s coefficient of determination (OLS: 0.087–0.666, MGWR: 0.616–0.894). The normalized difference vegetation index (NDVI) and nighttime light (NTL) are the most critical drivers of daytime and nighttime SUHIs, respectively. In terms of model bandwidth, population and Δfine particulate matter are typically global variables, while ΔNDVI, intercept (i.e., spatial context), and NTL are local variables. The nighttime coefficient of ΔNDVI is significantly negative in the more economically developed southern coastal region, while it is significantly positive in northwestern China. Our study not only improves the understanding of the complex drivers of SUHIs from a multiscale perspective but also provides a basis for urban heat island mitigation by more precisely identifying the heterogeneity of drivers.

Satellite-based vegetation indices (VIs) and Apparent Thermal Inertia (ATI) derived from temperature change provide valuable information for estimating evapotranspiration (LE) and detecting the onset and severity of drought. The modified satellite-based Priestley-Taylor (MS-PT) algorithm that we developed earlier, coupling both VI and ATI, is validated based on observed data from 40 flux towers distributed across the world on all continents. The validation results illustrate that the daily LE can be estimated with the Root Mean Square Error (RMSE) varying from 10.7 W/m² to 87.6 W/m², and with the square of correlation coefficient (R²) from 0.41 to 0.89 (p < 0.01). Compared with the Priestley-Taylor-based LE (PT-JPL) algorithm, the MS-PT algorithm improves the LE estimates at most flux tower sites. Importantly, the MS-PT algorithm is also satisfactory in reproducing the inter-annual variability at flux tower sites with at least five years of data. The R² between measured and predicted annual LE anomalies is 0.42 (p = 0.02). The MS-PT algorithm is then applied to detect the variations of long-term terrestrial LE over Three-North Shelter Forest Region of China and to monitor global land surface drought. The MS-PT algorithm described here demonstrates the ability to map regional terrestrial LE and identify global soil moisture stress, without requiring precipitation information.

Potential evapotranspiration (PET) is an important indicator of atmospheric evaporation demand and has been widely used to characterize hydrological change. However, sparse observations of pan evaporation (EP) prohibit the accurate characterization of the spatial and temporal patterns of PET over large spatial scales. In this study, we have estimated PET of China using the Penman-Monteith (PM) method driven by gridded reanalysis datasets to analyze the spatial and decadal variations of PET in China during 1982–2010. The results show that the estimated PET has decreased on average by 3.3 mm per year (p < 0.05) over China during 1982–1993, while PET began to increase since 1994 by 3.4 mm per year (p < 0.05). The spatial pattern of the linear trend in PET of China illustrates that a widely significant increasing trend in PET appears during 1982–2010 in Northwest China, Central China, Northeast China and South China while there are no obvious variations of PET in other regions. Our findings illustrate that incident solar radiation (Rs) is the largest contributor to the variation of PET in China, followed by vapor pressure deficit (VPD), air temperature (Tair) and wind speed (WS). However, WS is the primary factor controlling inter-annual variation of PET over Northwest China.

Accurate estimation of latent heat flux (LE) is critical in characterizing semiarid ecosystems. Many LE algorithms have been developed during the past few decades. However, the algorithms have not been directly compared, particularly over global semiarid ecosystems. In this paper, we evaluated the performance of five LE models over semiarid ecosystems such as grassland, shrub, and savanna using the Fluxnet dataset of 68 eddy covariance (EC) sites during the period 2000–2009. We also used a modern-era retrospective analysis for research and applications (MERRA) dataset, the Normalized Difference Vegetation Index (NDVI) and Fractional Photosynthetically Active Radiation (FPAR) from the moderate resolution imaging spectroradiometer (MODIS) products; the leaf area index (LAI) from the global land surface satellite (GLASS) products; and the digital elevation model (DEM) from shuttle radar topography mission (SRTM30) dataset to generate LE at region scale during the period 2003–2006. The models were the moderate resolution imaging spectroradiometer LE (MOD16) algorithm, revised remote sensing based Penman–Monteith LE algorithm (RRS), the Priestley–Taylor LE algorithm of the Jet Propulsion Laboratory (PT-JPL), the modified satellite-based Priestley–Taylor LE algorithm (MS-PT), and the semi-empirical Penman LE algorithm (UMD). Direct comparison with ground measured LE showed the PT-JPL and MS-PT algorithms had relative high performance over semiarid ecosystems with the coefficient of determination (R²) ranging from 0.6 to 0.8 and root mean squared error (RMSE) of approximately 20 W/m². Empirical parameters in the structure algorithms of MOD16 and RRS, and calibrated coefficients of the UMD algorithm may be the cause of the reduced performance of these LE algorithms with R² ranging from 0.5 to 0.7 and RMSE ranging from 20 to 35 W/m² for MOD16, RRS and UMD. Sensitivity analysis showed that radiation and vegetation terms were the dominating variables affecting LE Fluxes in global semiarid ecosystem.

Continuous high spatio-temporal resolution monitoring of evapotranspiration (ET) is critical for water resource management and the quantification of irrigation water efficiency at both global and local scales. However, available remote sensing satellites cannot generally provide ET data at both high spatial and temporal resolutions. Data fusion methods have been widely applied to estimate ET at a high spatio-temporal resolution. Nevertheless, most fusion methods applied to ET are initially used to integrate land surface reflectance, the spectral index and land surface temperature, and few studies completely consider the influencing factor of ET. To overcome this limitation, this paper presents an improved ET fusion method, namely, the spatio-temporal adaptive data fusion algorithm for evapotranspiration mapping (SADFAET), by introducing critical surface temperature (the corresponding temperature to decide soil moisture), importing the weights of surface ET-indicative similarity (the influencing factor of ET, which is estimated from remote sensing data) and modifying the spectral similarity (the differences in spectral characteristics of different spatial resolution images) for the enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM). We fused daily Moderate Resolution Imaging Spectroradiometer (MODIS) and periodic Landsat 8 ET data in the SADFAET for the experimental area downstream of the Heihe River basin from April to October 2015. The validation results, based on ground-based ET measurements, indicated that the SADFAET could successfully fuse MODIS and Landsat 8 ET data (mean percent error: −5%), with a root mean square error of 45.7 W/m², whereas the ESTARFM performed slightly worse, with a root mean square error of 50.6 W/m². The more physically explainable SADFAET could be a better alternative to the ESTARFM for producing ET at a high spatio-temporal resolution.

Evapotranspiration (ET) is the combination process of the surface evaporation and plant transpiration, which occur simultaneously, and it links the terrestrial water cycles, carbon cycles, and energy exchange. In this study, based on the observations from 242 global FLUXnet sites, with daily mean temperature, relative humidity, net radiation, wind speed, incoming shortwave radiation, maximum temperature, minimum temperature, normalized difference vegetation index, altitude, difference in temperature, and observed ET as input data, we used a support vector machine and a semiempirical algorithm to estimate the land surface daily ET at nine different vegetation-type sites. Subsequently, based on the meteorological reanalysis data combined with remote sensing data, we estimated regional land surface ET of China during 1982-2010. The results showed that, for all vegetation-type sites, when the predicted ET was validated with the eddy covariance measurements, the support vector machine algorithm undervalued ET while the semiempirical algorithm overvalued ET. When five indicators and the second classification method were selected, the semiempirical algorithm probably could explain 56%-76% of the land surface ET change, whereas the support vector machine algorithm probably could explain 71%-85%. The regional values of annual daily average ET varied from 5.8 to 110.5 W/m ² , and the land surface ET overall trend decreased from the southeast to the northwest in China.