The spatio-temporal pattern of surface ecological status affects urban thermal environment distribution significantly. Urban thermal pattern, however, is a complicated physical phenomenon involving a series of terrestrial environmental parameters. Thus, it is insufficient to employ only one ecological parameter for depicting the variation of land surface temperature (LST). This paper begins with the analysis of four ecological parameters' influence on LST using regression analysis, based on 24 Landsat images which cover Xi’an of China from 1992 to 2014. These four parameters include greenness degree (i.e. the soil adjusted vegetation index, SAVI), soil moisture degree (i.e. the normalized difference moisture index, NDMI), dryness degree (i.e. the normalized difference soil index, NDSI) and resident aggregation degree (i.e. the normalized difference build-up index, NDBI). Besides, contribution intensity index was introduced to investigate the contribution effect of four ecological parameters on LST, and a new ecological index, integrated ecological index (IEI), was founded using the principal component analysis technique to integratedly represent its spatial and mathematical correlations with LST. Results indicate that four ecological parameters all possessed pronounced performance in impacting LST pattern in all dates: SAVI and NDMI were found to be correlated negatively with LST, whereas NDBI and NDSI correlated positively with LST. Additionally, SAVI had a profound impact on LST distribution compared with the other three parameters, and there was the biggest heating contribution in the lowest SAVI category. Further finding suggests that IEI as a new ecological index can be used to integratedly estimate the spatio-temporal change of LST, manifesting a negative correlation with LST. Our study thinks that the comprehensive characterization of surface ecological status is conducive to benefit us to better understand the spatio-temporal mechanism of thermal environment and ecosystem and to help urban decision-makers to execute effective conservation policies for the ecosystem.

Spatiotemporal fusion methods provide an effective way to generate both high temporal and high spatial resolution data for monitoring dynamic changes of land surface. But existing fusion methods face two main challenges of monitoring the abrupt change events and accurately preserving the spatial details of objects. The Flexible Spatiotemporal DAta Fusion method (FSDAF) was proposed, which can monitor the abrupt change events, but its predicted images lacked intra-class variability and spatial details. To overcome the above limitations, this study proposed a comprehensive and automated fusion method, the Enhanced FSDAF (EFSDAF) method and tested it for Landsat–MODIS image fusion. Compared with FSDAF, the EFSDAF has the following strengths: (1) it considers the mixed pixels phenomenon of a Landsat image, and the predicted images by EFSDAF have more intra-class variability and spatial details; (2) it adjusts the differences between Landsat images and MODIS images; and (3) it improves the fusion accuracy in the abrupt change area by introducing a new residual index (RI). Vegetation phenology and flood events were selected to evaluate the performance of EFSDAF. Its performance was compared with the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM), the Spatial and Temporal Reflectance Unmixing Model (STRUM), and FSDAF. Results show that EFSDAF can monitor the changes of vegetation (gradual change) and flood (abrupt change), and the fusion images by EFSDAF are the best from both visual and quantitative evaluations. More importantly, EFSDAF can accurately generate the spatial details of the object and has strong robustness. Due to the above advantages of EFSDAF, it has great potential to monitor long-term dynamic changes of land surface.

Land surface temperature (LST) is a crucial biophysical parameter related closely to the land–atmosphere interface. Satellite thermal infrared measurement provides an effective method to derive LST on regional and global scales, but it is very hard to acquire simultaneously high spatiotemporal resolution LST due to its limitation in the sensor design. Recently, many LST downscaling and spatiotemporal image fusion methods have been widely proposed to solve this problem. However, most methods ignored the spatial heterogeneity of LST distribution, and there are inconsistent image textures and LST values over heterogeneous regions. Thus, this study aims to propose one framework to derive high spatiotemporal resolution LSTs in heterogeneous areas by considering the optimal selection of LST predictors, the downscaling of MODIS LST, and the spatiotemporal fusion of Landsat 8 LST. A total of eight periods of MODIS and Landsat 8 data were used to predict the 100-m resolution LST at prediction time t_p in Zhangye and Beijing of China. Further, the predicted LST at t_p was quantitatively contrasted with the LSTs predicted by the regression-then-fusion strategy, STARFM-based fusion, and random forest-based regression, and was validated with the actual Landsat 8 LST product at t_p. Results indicated that the proposed framework performed better in characterizing LST texture than the referenced three methods, and the root mean square error (RMSE) varied from 0.85 K to 2.29 K, and relative RMSE varied from 0.18 K to 0.69 K, where the correlation coefficients were all greater than 0.84. Furthermore, the distribution error analysis indicated the proposed new framework generated the most area proportion at 0~1 K in some heterogeneous regions, especially in artificial impermeable surfaces and bare lands. This means that this framework can provide a set of LST dataset with reasonable accuracy and a high spatiotemporal resolution over heterogeneous areas.

The Hindu Kush Himalayan (HKH) region is one of the most ecologically vulnerable regions in the world. Several studies have been conducted on the dynamic changes of grassland in the HKH region, but few have considered grassland net ecosystem productivity (NEP). In this study, we quantitatively analyzed the temporal and spatial changes of NEP magnitude and the influence of climate factors on the HKH region from 2001 to 2018. The NEP magnitude was obtained by calculating the difference between the net primary production (NPP) estimated by the Carnegie–Ames Stanford Approach (CASA) model and the heterotrophic respiration (Rh) estimated by the geostatistical model. The results showed that the grassland ecosystem in the HKH region exhibited weak net carbon uptake with NEP values of 42.03 gC∙m⁻²∙yr⁻¹, and the total net carbon sequestration was 0.077 Pg C. The distribution of NEP gradually increased from west to east, and in the Qinghai–Tibet Plateau, it gradually increased from northwest to southeast. The grassland carbon sources and sinks differed at different altitudes. The grassland was a carbon sink at 3000–5000 m, while grasslands below 3000 m and above 5000 m were carbon sources. Grassland NEP exhibited the strongest correlation with precipitation, and it had a lagging effect on precipitation. The correlation between NEP and the precipitation of the previous year was stronger than that of the current year. NEP was negatively correlated with temperature but not with solar radiation. The study of the temporal and spatial dynamics of NEP in the HKH region can provide a theoretical basis to help herders balance grazing and forage.

The soil dielectric constant model is essential for retrieving soil properties based on microwave remote sensing. However, the existing saline soil dielectric constant models perform poorly in simulating the dielectric constant of soil with high water content and salinity. In this study, the Wang Yueru (WYR) saline soil dielectric constant model, which was demonstrated to perform well in describing the effect of salinity and moisture on the dielectric constant, was validated based on experimental measurements of soil samples under different water content and salinity degrees. Furthermore, we adjusted the model form, refitted the empirical coefficient in the model, and finally acquired a two-stage model for simulating the soil dielectric constant. The enhanced model was validated under different soil moisture and salinity ranges using experimental measurements of soil samples. Compared to the original model, the proposed model exhibits a larger improvement in simulating the soil dielectric constant, and the RMSE of the simulated results dramatically decreased from 7.3 to 1.6, especially for soil with high salinity and water content. On this basis, a model suitable for L-band microwave was established. This model is of great significance for studying soil dielectric characteristics and retrieving soil parameters based on L-band data. Furthermore, this model can be used to retrieve soil salinity and water content using microwave remote sensing under a broadened application situation, such as in saline-alkali soils, wetlands, and salt marshes.

Soil moisture (SM), as a crucial input variable of land surface processes, plays a pivotal role in the global hydrological cycle. The aim of this paper is to examine the spatiotemporal variability in SM in the Heihe River Basin using all-weather land surface temperature (LST) and reanalysis land surface data. Initially, we downscaled and generated daily 1 km all-weather SM data (2020) for the Heihe River Basin. Subsequently, we investigated the spatial and temporal patterns of SM using geostatistical and time stability methods. The driving forces of the monthly SM were studied using the optimal parameter-based geographical detector (OPGD) model. The results indicate that the monthly mean values of the downscaled SM data range from 0.115 to 0.146, with a consistently lower SM content and suitable temporal stability throughout the year. Geostatistical analysis revealed that months with a higher SM level exhibit larger random errors and higher variability. Driving analysis based on the factor detector demonstrated that in months with a lower SM level, the q values of each driving factor are relatively small, and the primary driving factors are land cover and elevation. Conversely, in months with a higher SM level, the q values for each driving factor are larger, and the primary driving factors are the normalized difference vegetation index and LST. Furthermore, interaction detector analysis suggested that the spatiotemporal variation in SM is not influenced by a single driving factor but is the result of the interaction among multiple driving factors, with most interactions enhancing the combined effect of two factors.