Remotely sensed data, with high spatial and temporal resolutions, can hardly be provided by only one sensor due to the tradeoff in sensor designs that balance spatial resolutions and temporal coverage. However, they are urgently needed for improving the ability of monitoring rapid landscape changes at fine scales (e.g., 30 m). One approach to acquire them is by fusing observations from sensors with different characteristics (e.g., Enhanced Thematic Mapper Plus (ETM+) and Moderate Resolution Imaging Spectroradiometer (MODIS)). The existing data fusion algorithms, such as the Spatial and Temporal Data Fusion Model (STDFM), have achieved some significant progress in this field. This paper puts forward an Enhanced Spatial and Temporal Data Fusion Model (ESTDFM) based on the STDFM algorithm, by introducing a patch-based ISODATA classification method, the sliding window technology, and the temporal-weight concept. Time-series ETM+ and MODIS surface reflectance are used as test data for comparing the two algorithms. Results show that the prediction ability of the ESTDFM algorithm has been significantly improved, and is even more satisfactory in the near-infrared band (the contrasting average absolute difference [AAD]: 0.0167 vs. 0.0265). The enhanced algorithm will support subsequent research on monitoring land surface dynamic changes at finer scales.

Land cover mapping in mountainous areas is a notoriously challenging task due to the rugged terrain and high spatial heterogeneity of land surfaces as well as the frequent cloud contamination of satellite imagery. Taking Southwestern China (a typical mountainous region) as an example, this paper established a new HC-MMK approach (Hierarchical Classification based on Multi-source and Multi-temporal data and geo-Knowledge), which was especially designed for land cover mapping in mountainous areas. This approach was taken in order to generate a 30 m-resolution land cover product in Southwestern China in 2010 (hereinafter referred to as CLC-SW2010). The multi-temporal native HJ (HuanJing, small satellite constellation for disaster and environmental monitoring) CCD (Charge-Coupled Device) images, Landsat TM (Thematic Mapper) images and topographical data (including elevation, aspect, slope, etc.) were taken as the main input data sources. Hierarchical classification tree construction and a five-step knowledge-based interactive quality control were the major components of this proposed approach. The CLC-SW2010 product contained six primary categories and 38 secondary categories, which covered about 2.33 million km² (accounting for about a quarter of the land area of China). The accuracies of primary and secondary categories for CLC-SW2010 reached 95.09% and 87.14%, respectively, which were assessed independently by a third-party group. This product has so far been used to estimate the terrestrial carbon stocks and assess the quality of the ecological environments. The proposed HC-MMK approach could be used not only in mountainous areas, but also for plains, hills and other regions. Meanwhile, this study could also be used as a reference for other land cover mapping projects over large areas or even the entire globe.

The Three-River Source Region (TRSR) is an important area for the ecological security of China. Vegetation growth has been affected by the climate change, topography, and human activities in this area. However, few studies have focused on analyzing time series tendencies of vegetation change in various terrain conditions. To address this issue in the TRSR, this study explored vegetation stability, tendency, and sustainability with multiple methods (e.g., coefficient of variation, Theil-Sen median trend analysis, Mann-Kendall test, and Hurst index) based on the 2000–2017 Global LAnd Surface Satellite Leaf Area Index (GLASS LAI) product. The differentiation patterns of LAI variations and multiyear mean LAI value under different topographic factors were also investigated in combination with digital elevation model (DEM). The results showed that (1) the mean LAI value in the study area increased, with a linear tendency of 0.013·10 a⁻¹; (2) LAI values decreased from southeast to northwest in terms of spatial distribution and the CV indicated LAI variations were relatively stable; (3) the trend analysis revealed that the improved area of LAI accounted for 62.72% which was larger than the degraded area (37.28%), and hurst index revealed a weak anti-sustaining effect of the current tendencies; and (4) the increasing trend was found in multiyear mean LAI value as relief amplitude and slope increased, while LAI stability improved with increasing slope. They exhibited a clear regular pattern. Moreover, significant improvement in LAI generally occurred in low-altitude and flat areas. Finally, the overall improvement and sustainability of LAI improved when moving from sunny aspects to shady aspects, but the LAI stability decreased. Note that vegetation degradation was observed in some high slope areas and was further aggravated. This study is beneficial for revealing the spatial and temporal changes of LAI and their changing rules as a function of different topographic factors in the TRSR. Meanwhile, the results of this study provide theoretical support for sustainable development of this area.

Due to the spatial heterogeneity of land surfaces, downscaling is an important issue in the development of carbon cycle models when evaluating the role of ecosystems in the global carbon cycle. In this study, a downscaling algorithm was developed to model gross primary productivity (GPP) at 500 m in a time series over rugged terrain, which considered the effects of spatial heterogeneity on carbon flux simulations. This work was carried out for a mountainous area with an altitude ranging from 2606 to 4744 m over the Gongga Mountain (Sichuan Province, China). In addition, the Moderate Resolution Imaging Spectroradiometer (MODIS) GPP product at 1 km served as the primary dataset for the downscaling algorithm, and the 500 m MODIS GPP product was used as the reference dataset to evaluate the downscaled GPP results. Moreover, in order to illustrate the advantages and benefits of the proposed downscaling method, the downscaled results in this work, along with ordinary kriging downscaled results, spline downscaled results and inverse distance weighted (IDW) downscaled results, were compared to the MODIS GPP at 500 m. The results showed that (1) the GPP difference between the 500 m MODIS GPP and the proposed downscaled GPP results was primarily in the range of [−1, 1], showing that both vegetation heterogeneity factors (i.e., LAI) and topographic factors (i.e., altitude, slope and aspect) were useful for GPP downscaling; (2) the proposed downscaled results (R² = 0.89, RMSE = 1.03) had a stronger consistency with the 500 m MODIS GPP than those of the ordinary kriging downscaled results (R² = 0.43, RMSE = 1.36), the spline downscaled results (R² = 0.40, RMSE = 1.50) and the IDW downscaled results (R² = 0.42, RMSE = 1.10) for all Julian days; and (3) the inconsistency between MODIS GPP at 500 m and 1 km increased with the increase in altitude and slope. The proposed downscaling algorithm could provide a reference when considering the effects of spatial heterogeneity on carbon flux simulations and retrieving other fine resolution ecological-physiology parameters (e.g., net primary productivity and evaporation) over topographically complex terrains.

Land surface is normally considered as a mixture of soil and vegetation. Many applications, such as drought monitoring and crop-yield estimation, benefit from accurate retrieval of both soil and vegetation temperatures through satellite observation. A preliminary study has been conducted in this study on the estimation of land surface soil and vegetation component temperature using the geostationary satellite data acquired by Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard the Meteosat Second Generation (MSG) and TERRA-MODIS data. A synergetic algorithm is proposed to derive soil and vegetation temperatures by using the temporal and spatial information in SEVIRI and MODIS products. The approach is applied to both simulation data and satellite data. For simulation data, the component temperatures are well estimated with root mean squared error (RMSE) close to 0 K. For satellite data application, reasonable spatial distributions of the soil and vegetation temperatures are derived for eight cloud-free days in the Iberian Peninsula from June to August 2009. An evaluation is performed for the estimated vegetation temperature against the near surface air temperature. The correlation analysis between two datasets is found that the R-squareds are from 0.074 to 0.423 and RMSEs are within 4 K. Considering the impact of fraction of vegetation cover (FVC) on the validation, the pixels with FVC less than 30% are excluded in the total data comparison, and an obvious improvement is achieved with R-squared from 0.231 to 0.417 and RMSE from 2.9 K to 2.58 K. The validation indicates that the proposed algorithm is able to provide reasonable estimations of soil and vegetation temperatures. It is a potential way to map soil and vegetation temperature for large areas.

Increasing attention is being paid on leaf area index (LAI) retrieval in mountainous areas. Mountainous areas present extreme topographic variability, and are characterized by more spatial heterogeneity and inaccessibility compared with flat terrain. It is difficult to collect representative ground-truth measurements, and the validation of LAI in mountainous areas is still problematic. A cost-constrained sampling strategy (CSS) in support of LAI validation was presented in this study. To account for the influence of rugged terrain on implementation cost, a cost-objective function was incorporated to traditional conditioned Latin hypercube (CLH) sampling strategy. A case study in Hailuogou, Sichuan province, China was used to assess the efficiency of CSS. Normalized difference vegetation index (NDVI), land cover type, and slope were selected as auxiliary variables to present the variability of LAI in the study area. Results show that CSS can satisfactorily capture the variability across the site extent, while minimizing field efforts. One appealing feature of CSS is that the compromise between representativeness and implementation cost can be regulated according to actual surface heterogeneity and budget constraints, and this makes CSS flexible. Although the proposed method was only validated for the auxiliary variables rather than the LAI measurements, it serves as a starting point for establishing the locations of field plots and facilitates the preparation of field campaigns in mountainous areas.