Search

For many applications, raw satellite observations need to be converted to high-level products of various essential environmental variables. While numerous products are available at kilometer spatial resolutions, there are few global products at high spatial resolutions (10–30 m), which are also referred to fine or medium resolutions in the literature. To facilitate the development of more high spatial resolution products, this paper systematically reviews the state-of-the-art progress on inversion algorithms and publicly available regional and global products. We begin with an inventory of available high-resolution satellite data, and then present different algorithms for determining cloud masks, estimating aerosol optical depth, and performing atmospheric correction and topographic correction for land surface reflectance retrieval. The majority of this paper reviews the inversion algorithms and existing regional to global products of 18 variables in four major categories: 1) Land surface radiation, including broadband albedo, land surface temperature, and all-wave net radiation; 2) Terrestrial ecosystem variables, including leaf area index, fraction of absorbed photosynthetically active radiation, fractional vegetation cover, fractional forest cover, tree height, forest above-ground biomass gross primary production, net primary production, and agricultural crop yield; 3) Water cycle and cryosphere, including soil moisture, evapotranspiration, and snow cover; and 4) Land surface types, such as global land cover, impervious surface, inland water, crop type, and fire. Since the existing products over large regions are usually spatially discontinuous due to cloud contamination, different data fusion and data assimilation algorithms and some products for producing spatially seamless and temporally continuous products are presented. In the end, we discuss a variety of challenges in generating global high spatial resolution satellite products.

China has made remarkable progress in controlling particulate matter, while O3 pollution over China has become increasingly severe in recent years according to ground observations. Continuous monitoring of dynamic changes in O3 concentrations on regional and national scales can provide valuable insights for pollution control policies. Therefore, an improved similarity distance-based space-time random forest (SDSTRF) model was developed to estimate the near-surface O3 concentration using the surface measurements, satellite O3 precursors, meteorological variables, and other auxiliary information. The O3 concentration data over China were generated based on the developed model with a spatial resolution of 10 km and a temporal resolution of 1 day from 2016 to 2022. The validation results against the ground measurements indicate that the developed SDSTRF model effectively captures O3 variations, achieving a coefficient of determination of 0.83 and a root mean square error of 20.37 μg/m3. The spatiotemporal variations of O3 concentrations were investigated using the generated dataset. A significant increasing trend of 1.243 μg/m3/yr in O3 concentrations was observed in eastern China during the COVID-19 pandemic, which was attributed to changes in NOx concentrations. In this study, the possible reasons for the increase in O3 concentrations are also discussed. Overall, the improved SDSTRF model and the comprehensive analysis of the spatiotemporal variations of near-surface O3 will significantly contribute to achieving clean air in China.

Land surface all-wave net radiation (Rn) is crucial in determining Earth's climate by contributing to the surface radiation budget. This study evaluated seven satellite and three reanalysis long-term land surface Rn products under different spatial scales, spatial and temporal variations, and different conditions. The results showed that during 2000-2018, Global Land Surface Satellite Product (GLASS)-Moderate Resolution Imaging Spectroradiometer (MODIS) performed the best (RMSE=25.54 Wm-2, bias=−1.26 Wm-2), followed by ERA5 (the fifth-generation of European Centre for Medium-Range Weather Forecast Reanalysis) (RMSE=32.17 Wm-2, bias=−4.88 Wm-2) and GLASS-AVHRR (Advanced Very-High-Resolution Radiometer) (RMSE=33.10 Wm-2, bias=4.03 Wm-2). During 1983-2018, GLASS-AVHRR and ERA5 ranked top and performed similarly, with RMSE values of 31.70 and 33.08 Wm-2 and biases of −4.56 and 3.48 Wm-2, respectively. The averaged multi-annual mean Rn over the global land surface of satellite products was higher than that of reanalysis products by about 10∼30 Wm-2. These products differed remarkably in long-term trends variations, particularly pre-2000, but no significant trends were observed. Discrepancies were more frequent in satellite data, while reanalysis products showed smoother variations. Large discrepancies were found in regions with high latitudes, reflectance, and elevation which could be attributed to input radiative components, meteorological variables (e.g., cloud properties, aerosol optical thickness), and applicability of the algorithms used. While further research is needed for detailed insights.

Surface incident shortwave radiation (Rs) can promote the circulation of substance and energy, and the accuracy of its estimation is of great significance for climate studies. The Rs can be acquired from satellite retrievals, reanalysis predictions and general circulation model (GCM) simulations. Although Rs estimates have been evaluated and compared in previous studies, most of them focus on evaluating the Rs estimates over specific regions using ground measurements from limited stations. Therefore, it is essential to comprehensively validate Rs estimates from multiple data sources. In this study, ground measurements of 690 stations from BSRN, GEBA, CMA, GC-NET and buoys were employed to validate the Rs estimates from seven representative products (GLASS, GEWEX-SRB, CERES-EBAF, ERA5, MERRA2, CFSR and CMIP6). The validation results indicated that the selected products overestimated Rs globally, with biases ranged from 0.48 to 21.27 W/m2. The satellite retrievals showed relatively better accuracy among seven datasets compared to ground measurements at the selected stations. Moreover, the selected seven products were all in poor accuracy at high-latitude regions with RMSEs greater than 50 W/m2. The long-term variation trends were also analyzed in this study.

As a major agricultural hazard, drought frequently occurs due to a reduction in precipitation resulting in a continuously propagating soil moisture (SM) deficit. Assessment of the high spatial-resolution SM-derived drought index is crucial for monitoring agricultural drought. In this study, we generated a downscaled random forest SM dataset (RF-SM) and calculated the soil water deficit index (RF-SM-SWDI) at 30 m for agricultural drought monitoring. The results showed that the RF-SM dataset exhibited better consistency with in situ SM observations in the detection of extremes than did the SM products, including SMAP, SMOS, NCA-LDAS, and ESA CCI, for different land cover types in the U.S. and yielded a satisfactory performance, with the lowest root mean square error (RMSE, below 0.055 m3/m3) and the highest coefficient of determination (R2, above 0.8) for most observation networks, based on the number of sites. A vegetation health index (VHI), derived from a Landsat 8 optical remote sensing dataset, was also generated for comparison. The results illustrated that the RF-SM-SWDI and VHI exhibited high correlations (R ≥ 0.5) at approximately 70% of the stations. Furthermore, we mapped spatiotemporal drought monitoring indices in California. The RF-SM-SWDI provided drought conditions with more detailed spatial information than did the short-term drought blend (STDB) released by the U.S. Drought Monitor, which demonstrated the expected response of seasonal drought trends, while differences from the VHI were observed mainly in forest areas. Therefore, downscaled SM and SWDI, with a spatial resolution of 30 m, are promising for monitoring agricultural field drought within different contexts, and additional reliable factors could be incorporated to better guide agricultural management practices.

The Priestley–Taylor model of the Jet Propulsion Laboratory (PT-JPL) evapotranspiration (ET) model is relatively simple and has been widely used based on meteorological and satellite data. However, soil moisture (SM) constraints include a vapor pressure deficit (VPD) that causes large uncertainty. In this study, we proposed a PT-SinRH model by introducing a sine function of air relative humidity (RH) to replace RHVPD to characterize SM constraints, which can improve the accuracy of ET estimations. The PT-SinRH model is validated by eddy covariance (EC) data from 2000–2020. These data were collected by AmeriFlux at 28 sites on the conterminous United States (CONUS), and the land cover types of the sites vary from croplands to wetlands, grasslands, shrub lands and forests. The validation results from daily scale-based on-site and satellite data inputs showed that the PT-SinRH model estimates fit the observations with a coefficient of determination (R2) of 0.55, root-mean-square error (RMSE) of 17.5 W/m2, bias of −1.2 W/m2 and Kling–Gupta efficiency (KGE) of 0.70. Additionally, the PT-SinRH model based on reanalysis and satellite data inputs has an R2 of 0.49, an RMSE of 20.3 W/m2, a bias of −8.6 W/m2 and a KGE of 0.55. The PT-SinRH model showed better accuracy when using the site-measured meteorological data than when using reanalysis meteorological data as inputs. Additionally, compared with the PT-JPL model, the results demonstrate that our approach, i.e., PT-SinRH, improved ET estimates, increasing the R2 and KGE by 0.02 and decreasing the RMSE by about 0.6 W/m2. This simple but accurate method permits us to investigate the decadal variation in regional ET over the land.

Involving the effect of atmospheric CO2 fertilization is effective for improving the accuracy of estimating gross primary production (GPP) using light use efficiency (LUE) models. However, the widely used LUE model, the remote sensing-driven Carnegie–Ames–Stanford Approach (CASA) model, scarcely considers the effects of atmospheric CO2 fertilization, which causes GPP estimation uncertainties. Therefore, this study proposed an improved method for estimating GPP by integrating the atmospheric CO2 concentration into the CASA model and generated a long time series GPP dataset with high precision for the Qinghai–Tibet Plateau. The CASA model was improved by considering the impact of atmospheric CO2 on vegetation productivity and discerning variations in CO2 gradients within the canopy and leaves. A 500 m monthly GPP dataset for the Qinghai–Tibet Plateau from 2003 to 2020 was generated. The results showed that the improved GPP estimation model achieved better performances on estimating GPP (R2 = 0.68, RMSE = 406 g C/m2/year) than the original model (R2 = 0.67, RMSE = 499.32 g C/m2/year) and MODIS GPP products (R2 = 0.49, RMSE = 522.56 g C/m2/year). The GPP on the Qinghai–Tibet Plateau increased significantly with the increase in atmospheric CO2 concentration and the gradual accumulation of dry matter. The improved method can also be used for other regions and the generated GPP dataset is valuable for further understanding the ecosystem carbon cycles on the Qinghai–Tibet Plateau.

The detection and monitoring of burned areas is crucial for vegetation recovery, loss assessment, and anomaly analysis. Although vegetation indices (VIs) have been widely used, accurate vegetation detection is challenging due to potential confusion in the spectra of different types of land cover and the interference of shadow effects caused by terrain. In this work, a novel Vegetation Anomaly Spectral Texture Index (VASTI) is proposed, which leverages the merits of both spectral and spatial texture features to identify abnormal pixels for extracting burned vegetation areas. The performance of the VASTI and its components, the Global Environmental Monitoring Index (GEMI), the Enhanced Vegetation Index (EVI), and the texture feature Autocorrelation (AC) were assessed based on a global dataset previously established, which contains 1774 pairs of samples from 10 different sites. The results illustrated that, compared with the GEMI and EVI, the VASTI improved the user’s accuracy (UA), producer’s accuracy (PA), and kappa coefficient across the ten study areas by approximately 5% to 10%. Compared to AC, the VASTI improved the accuracy of abnormal vegetation detection by 13% to 25%. The improvements were mainly caused by the fact that the incorporation of texture features can reduce spectral confusion between pixels. The innovation of the VASTI is that it considers the relationship between anomalous pixels and surrounding pixels by explicitly integrating spatial texture features with traditional spectral features.

The detection of Martian polar ice cap change patterns is important for understanding their effects on driving Mars’s global water cycle and for regulating atmospheric circulation. However, current Martian ice cap identification using optical remote sensing data mainly relies on visual interpretation, which makes it difficult to quickly extract ice caps from multiple images and analyze their fine-scale spatiotemporal variation characteristics. Therefore, this study proposes an automatic Martian polar ice cap extraction algorithm for remote sensing data and analyzes the dynamic change characteristics of the Martian North Pole ice cap using time-series data. First, the automatic Martian ice cap segmentation algorithm was developed based on the ice cap features of high reflectance in the blue band and low saturation in the RGB band. Second, the Martian North Pole ice cap was extracted for the three Martian years MY25, 26, and 28 using Mars Orbiter Camera (MOC) Mars Daily Global Maps (MDGMs) data, which had better spatiotemporal continuity to analyze its variation characteristics. Lastly, the spatiotemporal variation characteristics of the ice cap and the driving factors of ice cap ablation were explored for the three aforementioned Martian years. The results indicated that the proposed automatic ice cap extraction algorithm had good performance, and the classification accuracy exceeded 93%. The ice cap ablation boundary retreat rates and spatiotemporal distributions were similar for the three years, with approximately 105 km2 of ice cap ablation for every one degree of areocentric longitude of the Sun (Ls). The main driving factor of ice cap ablation was solar radiation, which was mainly related to Ls. In addition, elevation had a different effect on ice cap ablation at different Ls in the same latitude area near the ablation boundary of the ice cap.

The all-wave net radiation (Rn) at the land surface represents surface radiation budget and plays an important role in the Earth's energy and water cycles. Many studies have been conducted to estimate from satellite top-of-atmosphere (TOA) data using various methods, particularly the application of machine learning (ML) and deep learning (DL). However, few studies have been conducted to provide a comprehensive evaluation about various ML and DL methods in retrieving. Based on extensive in situ measurements distributed at mid-low latitudes, the corresponding Moderate Resolution Imaging Spectroradiometer (MODIS) TOA observations, and the daily from the fifth generation of European Centre for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) used as a priori knowledge, this study assessed nine models for daily estimation, including six classic ML methods (random forest -RF, adaptive boosting - Adaboost, extreme gradient boosting -XGBoost, multilayer perceptron -MLP, radial basis function neural network -RBF, and support vector machine -SVM) and three DL methods (multilayer perceptron neural network with stacked autoencoders -SAE, deep belief network -DBN and residual neural network -ResNet). The validation results showed that the three DL methods were generally better than the six ML methods except XGBoost, although they all performed poorly in certain conditions such as winter days, rugged terrain, and high elevation. ResNet had the most robust performance across different land cover types, elevations, seasons, and latitude zones, but it has disadvantages in practice because of its highly configurable implementation environment and low computational efficiency. The estimated daily values from all nine models were more accurate than the corresponding Global LAnd Surface Satellite (GLASS) product.

It is a difficult undertaking to reliably estimate global terrestrial evapotranspiration (ET) using the Visible Infrared Imaging Radiometer Suite (VIIRS) at high spatial and temporal scales. We employ deep neural networks (DNN) to enhance the estimation of terrestrial ET on a global scale using satellite data. We accomplish this by merging five algorithms that are process-based and that make use of VIIRS data. These include the Shuttleworth–Wallace dual-source ET method (SW), the Priestley–Taylor-based ET algorithm (PT-JPL), the MOD16 ET product algorithm (MOD16), the modified satellite-based Priestley–Taylor ET algorithm (MS-PT), and the simple hybrid ET algorithm (SIM). We used 278 eddy covariance (EC) tower sites from 2012 to 2022 to validate the DNN approach, comparing it to Bayesian model averaging (BMA), gradient boosting regression tree (GBRT) and random forest (RF). The validation results demonstrate that the DNN significantly improves the accuracy of daily ET estimates when compared to three other merging methods, resulting in the highest average determination coefficients (R2, 0.71), RMSE (21.9 W/m2) and Kling–Gupta efficiency (KGE, 0.83). Utilizing the DNN, we generated a VIIRS ET product with a 500 m spatial resolution for the years 2012–2020. The DNN method serves as a foundational approach in the development of a sustained and comprehensive global terrestrial ET dataset. The basis for characterizing and analyzing global hydrological dynamics and carbon cycling is provided by this dataset.

With the rapid development of economic situation, global warming and population growing, drought happened in many areas. Based on linear dual source remote sensing model we establish the modified hybrid linear dual-source’s remote sensing evaporation model to simulate the Sichuan Province actual evaporation. Then we calculate the potential evaporation by using the Hargreaves formula. In addition, we analyze evaporation and drought trend of Sichuan province. This study showed that: (1) the actual variation trend of evaporation in Sichuan Province is decreasing from 2001 to 2015. The actual evaporation in the east is higher than that in the west regions, and in the south is higher than that in the north regions. (2) The value of potential evaporation fluctuates is not obviously in every season. According to the spatial distribution, the potential evaporation distribution is decreasing from south to north. (3) The evaporation of drought index shows an upward trend, but the fluctuation phenomenon of Sichuan Province in the past seasons is not obvious. This climate research provide the basis for the prediction and guidance of agricultural activities in the region.

Measurement(s) The all-wave net radiation Technology Type(s) fusion method

Fractional vegetation cover (FVC) is a critical land surface parameter, and several large-scale FVC products have been generated based on remote sensing data. Among these existing products, the global land surface satellite (GLASS) FVC product, derived from moderate resolution imaging spectroradiometer (MODIS) 500 m reflectance data (MOD09A1), has achieved complete spatial-temporal continuity and satisfying accuracy. To further improve the spatial resolution of GLASS FVC product, this study developed a novel FVC estimation algorithm for MODIS 250 m reflectance data based on a recurrent neural network with the long short-term memory unit (RNN-LSTM). The RNN-LSTM was established using sequence training samples derived from the MODIS 250 m reflectance and GLASS FVC products, which were conducted over three vegetation types in mid-West China. Additionally, two machine learning methods, including the back propagation neural network (BPNN) and multivariate adaptive regression splines (MARS), were used to compare with the proposed method. The evaluation results showed that RNN-LSTM derived FVC had reliable spatial-temporal continuity and good consistency with the GLASS FVC product. Furthermore, the smooth temporal profiles of the RNN-LSTM FVC estimation indicated that the proposed method was capable of capturing the temporal characteristics of vegetation growth and reducing the uncertainties from the atmosphere and radiation. Finally, an independent validation case in the Heihe area indicated that the RNN-LSTM algorithm achieved the best accuracy (R2 = 0.8081, rmse = 0.0951) compared with the BPNN (R2 = 0.7320, rmse = 0.1127) and MARS (R2 = 0.7361, rmse = 0.1117). This study provides a new approach by showing the potential of the RNN-LSTM method for land surface parameter estimation and related research.

A fractional vegetation cover (FVC) estimation method incorporating a vegetation growth model and a radiative transfer model was previously developed, which was suitable for FVC estimation in homogeneous areas because the finer-resolution pixels corresponding to one coarse-resolution FVC pixel were all assumed to have the same vegetation growth model. However, this assumption does not hold over heterogeneous areas, meaning that the method cannot be applied to large regions. Therefore, this study proposes a finer spatial resolution FVC estimation method applicable to heterogeneous areas using Landsat 8 Operational Land Imager reflectance data and Global LAnd Surface Satellite (GLASS) FVC product. The FVC product was first decomposed according to the normalized difference vegetation index from the Landsat 8 OLI data. Then, independent dynamic vegetation models were built for each finer-resolution pixel. Finally, the dynamic vegetation model and a radiative transfer model were combined to estimate FVC at the Landsat 8 scale. Validation results indicated that the proposed method (R2 = 0.7757, RMSE = 0.0881) performed better than either the previous method (R2 = 0.7038, RMSE = 0.1125) or a commonly used method involving look-up table inversions of the PROSAIL model (R2 = 0.7457, RMSE = 0.1249).

Spatiotemporal fusion has provided a feasible way to generate fractional vegetation cover (FVC) data with high spatial and temporal resolution. However, when the currently available spatiotemporal fusion methods are applied over agricultural regions, they usually underestimate high FVC values at the peak vegetation growth stage with medium FVC values as base data. This mainly results from inconsistencies in the temporal variations between fine- and coarse-resolution data if substantial temporal changes occur in vegetation. Therefore, a Spatial and Temporal Fusion method combining with Vegetation Growth Models (STF-VGM) was proposed to address this problem in this study, which incorporates vegetation growth models into the fusion process. Unlike other spatiotemporal fusion methods that mainly rely on changes in coarse-resolution data for prediction, STF-VGM fully utilizes available coarse- and fine-resolution time series data, including uncontaminated information in cloud/cloud shadow contaminated images. By establishing vegetation growth models with time series data, a conversion relationship between coarse- and fine-resolution FVC that changes along with the nonlinear vegetation change process can be extracted. STF-VGM makes prediction based on this variable relationship. A typical agricultural region located in the North China Plain was selected as the study area. The validation results indicated that the prediction accuracy for high FVC values was significantly improved using STF-VGM compared to the commonly used Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM) and Flexible Spatiotemporal DAta Fusion (FSDAF) methods (STF-VGM: coefficient of determination (R2) = 0.9491, root mean square error (RMSE) = 0.0650, average difference (AD) = -0.0092; ESTARFM: R2 = 0.9341, RMSE = 0.1127, AD = -0.0631; FSDAF: R2 = 0.9224, RMSE = 0.1110, AD = -0.0599). The satisfactory performance of STF-VGM was also achieved in predicting FVC values at other vegetation growth stages (early growth stage: R2 = 0.8277, RMSE = 0.0440, AD = 0.0027; rapid growth stage: R2 = 0.9183, RMSE = 0.0844, AD = 0.0500). In addition, STF-VGM also has the potential to improve the spatiotemporal fusion accuracy of other vegetation parameters and vegetation indices, which will be evaluated in the future.

Advanced remote sensing technology has provided spatially distributed variables for estimating land–ocean heat fluxes, allowing for practical applications in drought monitoring, water resources management, and climate assessment. This Special Issue includes several research studies using state-of-the-art algorithms for estimating downward longwave radiation, surface net radiation, latent heat flux, columnar atmospheric water vapor, fractional vegetation cover, and grassland aboveground biomass. This Special Issue intends to help scientists involved in global change research and practices better comprehend the strengths and disadvantages of the application of remote sensing for monitoring surface energy, water, and carbon budgets. The studies published in this Special Issue can be applied by natural resource management communities to enhance the characterization and assessment of land–ocean biophysical variables, as well as for more accurately partitioning heat flux into soil and vegetation based on the existing and forthcoming remote sensing data.

Validation of remotely sensed evapotranspiration (RS_ET) products is important because their accuracy is critical for various scientific applications. In this study, an integrated validation framework was proposed for evaluating RS_ET products with coarse spatial resolution extending from homogenous to heterogeneous land surfaces. This framework was applied at the pixel and river basin scales, using direct and indirect validation methods with multisource validation datasets, which solved the spatial mismatch between ground measurements and remotely sensed products. The accuracy, rationality of spatiotemporal variations, and error sources of RS_ET products and uncertainties during the validation process were the focuses in the framework. The application of this framework is exemplified by validating five widely used RS_ET products (i.e., GLEAM, DTD, MOD16, ETMonitor, and GLASS) in the Heihe River Basin from 2012 to 2016. Combined with the results from direct (as the priority method) and indirect validation (as the auxiliary method), DTD showed the highest accuracy (1-MAPE) in the vegetation growing season (75%), followed by ETMonitor (71%), GLASS (68%), GLEAM (54%), and MOD16 (44%). Each product reasonably reflected the spatiotemporal variations in the validation dataset. ETMonitor exhibited the highest consistency with the ground truth ET at the basin scale (ETMap) (R = 0.69), followed by GLASS (0.65), DTD (0.63), MOD16 (0.62), and GLEAM (0.57). Error sources of these RS_ET products were mainly due to the limitations of the algorithms and the coarse spatial resolution of the input data, while the uncertainties in the validation process amounted to 15–28%. This work is proposed to effectively validate and improve the RS_ET products over heterogeneous land surfaces.

Abstract Surface downward longwave radiation (SDLR) is a major component of the energy budget. Although studies have reported the spatiotemporal variations of SDLR in China, the spatiotemporal coverage of the situ measurements used is always limited. In this study, the gradient boosting regression tree (GBRT) was developed to reconstruct SDLR based on air temperature (Ta), relative humidity (RH), and downward shortwave radiation (DSR). Ground measurements collected at the Baseline Surface Radiation Network (BSRN) and the Arid and Semi‐arid Region Collaborative Observation Project (ASRCOP) were used to build and evaluate the GBRT model. The evaluation results showed that the daily SDLR estimates are correlated well with the SDLR in situ, with an overall root mean square errors (RMSE) of 16.5 Wm−2 and a correlation coefficient (R) value of 0.91 for the validation data set. Comparison with existing SDLR products showed that accuracy and trends of the SDLR estimates based on the GBRT method are reasonable. To obtain long‐term SDLR data for spatiotemporal analysis over China, densely distributed reconstructed DSR and ground measured Ta and RH collected at 756 Chinese Meteorological Administration (CMA) stations were used as input to estimate the SDLR based on the GBRT method over China during 1958–2015. The long‐term estimated SDLRs at the selected 563 stations showed that SDLR increased at an average rate of 1.3 Wm−2 per decade over China from 1958 to 2015. The trend of SDLR is positively correlated with the trend in Ta and water vapor pressure, whereas negatively correlated with the trend in DSR.

Abstract Soil moisture is a key factor affecting crop growth and crop yield. Information on soil moisture is critical for crop growth monitoring and yield estimation. This study evaluated the OPtical TRApezoid Model (OPTRAM) on soil moisture estimates in the Sanjiang Plain, China, using Moderate Resolution Imaging Spectroradiometer (MODIS) data and in situ soil moisture data acquired from May to September in 2016 and 2017. This model was empirically calculated by parameterizing the relationship between the normalized difference vegetation index (NDVI) and shortwave infrared transformed reflectance (STR). In addition, the precipitation data were collected and used for verification of the results. According to the scatterplots of STR and NDVI, OPTRAM can estimate soil moisture from July to September, while it's inapplicable in May and June. Evaluation results indicated that OPTRAM‐based soil moisture estimates provide overall RMSE from 0.05 to 0.13 cm3/cm3, bias from −0.11 to 0.06 cm3/cm3, and R2 from 0.10 to 0.50, respectively, for all investigated sites. The performance of OPTRAM compared with the precipitation data showed good agreement. In conclusion, OPTRAM could be used to estimate soil moisture with reasonable accuracy in most areas from July to September in Sanjiang Plain.

Accurate estimation of satellite-derived ocean latent heat flux (LHF) at high spatial resolution remains a major challenge. Here, we estimate monthly ocean LHF at 4 km spatial resolution over 5 years using bulk algorithm COARE 3.0, driven by satellite data and meteorological variables from reanalysis. We validated the estimated ocean LHF by multiyear observations and by comparison with seven ocean LHF products. Validation results from monthly observations at 96 widely distributed buoy sites from three buoy site arrays (TAO, PIRATA, and RAMA) indicated a bias of less than 7 W/m2 with R2 of more than 0.80 (p

Metropolis’ contribution (anthropogenic aerosols) to solar dimming and brightening remains a hot topic of special concern over the past several decades. However, urbanization effects on surface incident solar radiation ( R _s ) have not been comprehensively investigated. In this study, the urbanization effects on solar dimming and brightening were addressed using the densely distributed reconstructed R _s data at 375 stations and 92 urban–rural station pairs over the time period of 1960–2019 in China. The results indicate that the impacts of urbanization on the monthly mean R _s is 0.86 ± 7.99 W m ^−2 during the study period, while the impact is 0.90 ± 8.30 W m ^−2 and 0.82 ± 8.26 W m ^−2 for the solar dimming (1960–1992) and brightening (1992–2019) phase, respectively. The urbanization effects on the trend of R _s is −0.39 and 0.16 W m ^−2 per decade during dimming and brightening, respectively. It also found that urbanization effects on R _s trend differs strikingly in magnitudes for specific regions in China. Generally, urbanization speeds up China’s dimming in the dimming phase and slows down China’s brightening in the brightening phase.

The surface all-wave net radiation (Rn) plays an important role in the energy and water cycles, and most studies of Rn estimations have been conducted using satellite data. As one of the most commonly used satellite data sets, Moderate Resolution Imaging Spectroradiometer (MODIS) data have not been widely used for radiation calculations at mid-low latitudes because of its very low revisit frequency. To improve the daily Rn estimation at mid-low latitudes with MODIS data, four models, including three models built with random forest (RF) and different temporal expansion models and one model built with the look-up-table (LUT) method, are used based on comprehensive in situ radiation measurements collected from 340 globally distributed sites, MODIS top-of-atmosphere (TOA) data, and the fifth generation of European Centre for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) data from 2000 to 2017. After validation against the in situ measurements, it was found that the RF model based on the constraint of the daily Rn from ERA5 (an RF-based model with ERA5) performed the best among the four proposed models, with an overall validated root-mean-square error (RMSE) of 21.83 Wm−2, R2 of 0.89, and a bias of 0.2 Wm−2. It also had the best accuracy compared to four existing products (Global LAnd Surface Satellite Data (GLASS), Clouds and the Earth’s Radiant Energy System Edition 4A (CERES4A), ERA5, and FLUXCOM_RS) across various land cover types and different elevation zones. Further analyses illustrated the effectiveness of the model by introducing the daily Rn from ERA5 into a “black box” RF-based model for Rn estimation at the daily scale, which is used as a physical constraint when the available satellite observations are too limited to provide sufficient information (i.e., when the overpass time is less than twice per day) or the sky is overcast. Overall, the newly-proposed RF-based model with ERA5 in this study shows satisfactory performance and has strong potential to be used for long-term accurate daily Rn global mapping at finer spatial resolutions (e.g., 1 km) at mid-low latitudes.

Surface upward longwave radiation (SULR) is an indicator of thermal conditions over the Earth’s surface. In this study, we validated the simulated SULR from 51 Coupled Model Intercomparison Project (CMIP6) general circulation models (GCMs) through a comparison with ground measurements and satellite-retrieved SULR from the Clouds and the Earth’s Radiant Energy System, Energy Balanced and Filled (CERES EBAF). Moreover, we improved the SULR estimations by a fusion of multiple CMIP6 GCMs using multimodel ensemble (MME) methods. Large variations were found in the monthly mean SULR among the 51 CMIP6 GCMs; the bias and root mean squared error (RMSE) of the individual CMIP6 GCMs at 133 sites ranged from −3 to 24 W m−2 and 22 to 38 W m−2, respectively, which were higher than those found between the CERES EBAF and GCMs. The CMIP6 GCMs did not improve the overestimation of SULR compared to the CMIP5 GCMs. The Bayesian model averaging (BMA) method showed better performance in simulating SULR than the individual GCMs and simple model averaging (SMA) method, with a bias of 0 W m−2 and an RMSE of 19.29 W m−2 for the 133 sites. In terms of the global annual mean SULR, our best estimation for the CMIP6 GCMs using the BMA method was 392 W m−2 during 2000–2014. We found that the SULR varied between 386 and 393 W m−2 from 1850 to 2014, exhibiting an increasing tendency of 0.2 W m−2 per decade (p < 0.05).

Ocean surface net radiation (Rn) is significant in research on the Earth’s heat balance systems, air–sea interactions, and other applications. However, there have been few studies on Rn until now. Based on radiative and meteorological measurements collected from 66 globally distributed moored buoys, it was found that Rn was dominated by downward shortwave radiation (Rg↓) when the length ratio of daytime (LRD) was greater than 0.4 but dominated by downward longwave radiation (Rl↓) for the other cases (LRD ≤ 0.4). Therefore, an empirical scheme that includes two conditional models named Case 1 (LRD > 0.4) utilizing Rg↓ as a major input and Case 2 (LRD ≤ 0.4) utilizing Rl↓ as a major input for Rn estimation was successfully developed. After validation against in situ Rn, the performance of the empirical scheme was satisfactory with an overall R2 value of 0.972, an RMSE of 9.768 Wm−2, and a bias of −0.092 Wm−2. Specifically, the accuracies of the two conditional models were also very good, with RMSEs of 9.805 and 2.824 Wm−2 and biases of −0.095 and 0.346 Wm−2 for the Case 1 and Case 2 models, respectively. However, due to the limited number of available samples, the performances of these new models were poor in coastal and high-latitude areas, and the models did not work when the LRD was too small (i.e., LRD < 0.3). Overall, the newly developed empirical scheme for Rn estimation has strong potential to be widely used in practical use because of its simple format and high accuracy.

Coarse spatial resolution sensors play a major role in capturing temporal variation, as satellite images that capture fine spatial scales have a relatively long revisit cycle. The trade-off between the revisit cycle and spatial resolution hinders the access of terrestrial latent heat flux (LE) data with both fine spatial and temporal resolution. In this paper, we firstly investigated the capability of an Extremely Randomized Trees Fusion Model (ERTFM) to reconstruct high spatiotemporal resolution reflectance data from a fusion of the Chinese GaoFen-1 (GF-1) and the Moderate Resolution Imaging Spectroradiometer (MODIS) products. Then, based on the merged reflectance data, we used a Modified-Satellite Priestley–Taylor (MS–PT) algorithm to generate LE products at high spatial and temporal resolutions. Our results illustrated that the ERTFM-based reflectance estimates showed close similarity with observed GF-1 images and the predicted NDVI agreed well with observed NDVI at two corresponding dates (r = 0.76 and 0.86, respectively). In comparison with other four fusion methods, including the widely used spatial and temporal adaptive reflectance fusion model (STARFM) and the enhanced STARFM, ERTFM had the best performance in terms of predicting reflectance (SSIM = 0.91; r = 0.77). Further analysis revealed that LE estimates using ERTFM-based data presented more detailed spatiotemporal characteristics and provided close agreement with site-level LE observations, with an R2 of 0.81 and an RMSE of 19.18 W/m2. Our findings suggest that the ERTFM can be used to improve LE estimation with high frequency and high spatial resolution, meaning that it has great potential to support agricultural monitoring and irrigation management.

The long-term estimation of grassland aboveground biomass (AGB) is important for grassland resource management in the Three-River Headwaters Region (TRHR) of China. Due to the lack of reliable grassland AGB datasets since the 1980s, the long-term spatiotemporal variation in grassland AGB in the TRHR remains unclear. In this study, we estimated AGB in the grassland of 209,897 km2 using advanced very high resolution radiometer (AVHRR), MODerate-resolution Imaging Spectroradiometer (MODIS), meteorological, ancillary data during 1982–2018, and 75 AGB ground observations in the growth period of 2009 in the TRHR. To enhance the spatial representativeness of ground observations, we firstly upscaled the grassland AGB using a gradient boosting regression tree (GBRT) model from ground observations to a 1 km spatial resolution via MODIS normalized difference vegetation index (NDVI), meteorological and ancillary data, and the model produced validation results with a coefficient of determination (R2) equal to 0.76, a relative mean square error (RMSE) equal to 88.8 g C m−2, and a bias equal to −1.6 g C m−2 between the ground-observed and MODIS-derived upscaled AGB. Then, we upscaled grassland AGB using the same model from a 1 km to 5 km spatial resolution via AVHRR NDVI and the same data as previously mentioned with the validation accuracy (R2 = 0.74, RMSE = 57.8 g C m−2, and bias = −0.1 g C m−2) between the MODIS-derived reference and AVHRR-derived upscaled AGB. The annual trend of grassland AGB in the TRHR increased by 0.37 g C m−2 (p < 0.05) on average per year during 1982–2018, which was mainly caused by vegetation greening and increased precipitation. This study provided reliable long-term (1982–2018) grassland AGB datasets to monitor the spatiotemporal variation in grassland AGB in the TRHR.

Fire is one of the most widespread and destructive disasters, which causes property losses, casualties, and disruption of the balance of ecosystems. Therefore, it is highly necessary for firefighting to study the variations in fire and its climatic attributions. This study analyzed the characteristics of fire-burned area (BA) and its response to climatic factors in seven subregions of China from 2001 to 2018 using satellite remote sensing BA products. The results show that the BA in China and most of its subregions shows a decreasing trend. In general, it is negatively correlated with precipitation and positively correlated with air temperature and wind speed based on the regression and correlation analyses. Based on Pearson correlation and random forest methods, it is also found that the temperature is commonly an important factor contributing to BA in China, except for R2 (Inner Mongolia region), where wind speed is more important, and R5 (South China), where precipitation is more important, which coexists at annual and seasonal scales. Besides temperature, precipitation in spring and summer is the main driving factor, such as in R1 (Northeast China), R5, R6 (Northwest China) and R7 (Qinghai–Tibet Plateau) in spring and R4 (Central China), R5 and R7 in summer; and wind speed in autumn and winter is the main driving factor, such as in R2 and R4 in autumn and R2, R3, R5, R6 and R7 in winter. Finally, the distributions of BA with respect to each climatic factor were also analyzed to quantify the range of climatic factors with maximum BA occurrence.

An accurate and spatially continuous estimation of terrestrial latent heat flux (LE) is fundamental and crucial for the rational utilization of water resources in the Haihe River Basin (HRB). However, the sparsity of flux observation sites hinders the accurate characterization of spatiotemporal LE patterns over the HRB. In this study, we estimated the daily LE across the HRB using the gradient boosting regression tree (GBRT) from global land surface satellite NDVI data, reanalysis data and eddy covariance data. Compared with the random forests (RF) and extra tree regressor (ETR) methods, the GBRT obtains the best results, with R2 = 0.86 and root mean square error (RMSE = 18.1 W/m2. Then, we applied the GBRT algorithm to map the average annual terrestrial LE of the HRB from 2016 to 2018 with a spatial resolution of 0.05°. When compared with the Global Land Surface Satellite (GLASS) and Moderate Resolution Imaging Spectroradiometer (MODIS) LE products, the difference between the terrestrial LE estimated by the GBRT algorithm and the GLASS and MODIS products was less than 20 W/m2 in most areas; thus, the GBRT algorithm was reliable and reasonable for estimating the long-term LE estimation over the HRB.

Extreme hydrologic events are getting more frequent under a changing climate, and a reliable hydrological modeling framework is important to understand their mechanism. However, existing hydrological modeling frameworks are mostly constrained to a relatively coarse resolution, unrealistic input information, and insufficient evaluations, especially for the large domain, and they are, therefore, unable to address and reconstruct many of the water-related issues (e.g., flooding and drought). In this study, a 0.0625-degree (~6 km) resolution variable infiltration capacity (VIC) model developed for China from 1970 to 2016 was extensively evaluated against remote sensing and ground-based observations. A unique feature in this modeling framework is the incorporation of new remotely sensed vegetation and soil parameter dataset. To our knowledge, this constitutes the first application of VIC with such a long-term and fine resolution over a large domain, and more importantly, with a holistic system-evaluation leveraging the best available earth data. The evaluations using in-situ observations of streamflow, evapotranspiration (ET), and soil moisture (SM) indicate a great improvement. The simulations are also consistent with satellite remote sensing products of ET and SM, because the mean differences between the VIC ET and the remote sensing ET range from −2 to 2 mm/day, and the differences for SM of the top thin layer range from −2 to 3 mm. Therefore, this continental-scale hydrological modeling framework is reliable and accurate, which can be used for various applications including extreme hydrological event detections.

An operational and accurate model for estimating global or regional terrestrial latent heat of evapotranspiration (ET) across different land-cover types from satellite data is crucial. Here, a simplified Priestley–Taylor (SPT) model was developed without surface net radiation (Rn) by combining incident shortwave radiation (Rs), satellite vegetation index, and air relative humidity (RH). Ground-measured ET for 2000–2009 collected by 100 global FLUXNET eddy covariance (EC) sites was used to calibrate and evaluate the SPT model. A series of cross-validations demonstrated the reasonable performance of the SPT model to estimate seasonal and spatial ET variability. The coefficients of determination (R2) of the estimated versus observed daily (monthly) ET ranged from 0.42 (0.58) (p < 0.01) at shrubland (SHR) flux sites to 0.81 (0.86) (p < 0.01) at evergreen broadleaf forest (EBF) flux sites. The SPT model was applied to estimate agricultural ET at high spatial resolution (16 m) from Chinese Gaofen (GF)-1 data and monitor long-term (1982–2018) ET variations in the Three-River Headwaters Region (TRHR) of mainland China using the Global LAnd-Surface Satellite (GLASS) normalized difference vegetation index (NDVI) product. The proposed SPT model without Rn provides an alternative model for estimating regional terrestrial ET across different land-cover types.

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 (R2) ranging from 0.6 to 0.8 and root mean squared error (RMSE) of approximately 20 W/m2. 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 R2 ranging from 0.5 to 0.7 and RMSE ranging from 20 to 35 W/m2 for MOD16, RRS and UMD. Sensitivity analysis showed that radiation and vegetation terms were the dominating variables affecting LE Fluxes in global semiarid ecosystem.

The downward longwave radiation (Ld, 4–100 μm) is a major component of research for the surface radiation energy budget and balance. In this study, we applied five machine learning methods, namely artificial neural network (ANN), support vector regression (SVR), gradient boosting regression tree (GBRT), random forest (RF), and multivariate adaptive regression spline (MARS), to estimate Ld using ground measurements collected from 27 Baseline Surface Radiation Network (BSRN) stations. Ld measurements in situ were used to validate the accuracy of Ld estimation models on daily and monthly time scales. A comparison of the results demonstrated that the estimates on the basis of the GBRT method had the highest accuracy, with an overall root-mean-square error (RMSE) of 17.50 W m−2 and an R value of 0.96 for the test dataset on a daily time scale. These values were 11.19 W m−2 and 0.98, respectively, on a monthly time scale. The effects of land cover and elevation were further studied to comprehensively evaluate the performance of each machine learning method. All machine learning methods achieved better results over the grass land cover type but relatively worse results over the tundra. GBRT, RF, and MARS methods were found to show good performance at both the high- and low-altitude sites.

Downward surface solar radiation (Rs) plays a dominant role in determining the climate and environment on the Earth. However, the densely distributed ground observations of Rs are usually insufficient to meet the increasing demand of the climate diagnosis and analysis well, so it is essential to build a long-term accurate Rs dataset. The extremely randomized trees (ERT) algorithm was used to generate Rs using routine meteorological observations (2000–2015) from the Climate Data Center of the Chinese Meteorological Administration (CDC/CMA). The estimated Rs values were validated against ground measurements at the national scale with an overall correlation coefficient value of 0.97, a mean bias of 0.04 Wm−2, a root-mean-square-error value of 23.12 Wm−2, and a mean relative error of 9.81%. It indicates that the estimated Rs from the ERT-based model is reasonably accurate. Moreover, the ERT-based model was used to generate a new daily Rs dataset at 756 CDC/CMA stations from 1958 to 2015. The long-term variation trends of Rs at 454 stations covering 46 consecutive years (1970–2015) were also analyzed. The Rs in China showed a significant decline trend (−1.1 Wm−2 per decade) during 1970–2015. A decreasing trend (−2.8 Wm−2 per decade) in Rs during 1970–1992 was observed, followed by a recovery trend (0.23 Wm−2 per decade) during 1992–2015. The recovery trends at individual stations were found at 233 out of 454 stations during 1970–2015, which were mainly located in southern and northern China. The new Rs dataset would substantially provide basic data for the related studies in agriculture, ecology, and meteorology.

Surface all-wave net radiation (Rn) is a crucial variable driving many terrestrial latent heat (LE) models that estimate global LE. However, the differences between different Rn products and their impact on global LE estimates still remain unclear. In this study, we evaluated two Rn products, Global LAnd Surface Satellite (GLASS) beta version Rn and Modern-Era Retrospective Analysis for Research and Applications-version 2 (MERRA-2) Rn, from 2007–2017 using ground-measured data from 240 globally distributed in-situ radiation measurements provided by FLUXNET projects. The GLASS Rn product had higher accuracy (R2 increased by 0.04–0.26, and RMSE decreased by 2–13.3 W/m2) than the MERRA-2 Rn product for all land cover types on a daily scale, and the two Rn products differed greatly in spatial distribution and variations. We then determined the resulting discrepancies in simulated annual global LE using a simple averaging model by merging five diagnostic LE models: RS-PM model, SW model, PT-JPL model, MS-PT model, and SIM model. The validation results showed that the estimated LE from the GLASS Rn had higher accuracy (R2 increased by 0.04–0.14, and RMSE decreased by 3–8.4 W/m2) than that from the MERRA-2 Rn for different land cover types at daily scale. Importantly, the mean annual global terrestrial LE from GLASS Rn was 2.1% lower than that from the MERRA-2 Rn. Our study showed that large differences in satellite and reanalysis Rn products could lead to substantial uncertainties in estimating global terrestrial LE.

Ocean latent heat flux (LHF) is an essential variable for air–sea interactions, which establishes the link between energy balance, water and carbon cycle. The low-latitude ocean is the main heat source of the global ocean and has a great influence on global climate change and energy transmission. Thus, an accuracy estimation of high-resolution ocean LHF over low-latitude area is vital to the understanding of energy and water cycle, and it remains a challenge. To reduce the uncertainties of individual LHF products over low-latitude areas, four machine learning (ML) methods (Artificial Neutral Network (ANN), Random forest (RF), Bayesian Ridge regression and Random Sample Consensus (RANSAC) regression) were applied to estimate low-latitude monthly ocean LHF by using two satellite products (JOFURO-3 and GSSTF-3) and two reanalysis products (MERRA-2 and ERA-I). We validated the estimated ocean LHF using 115 widely distributed buoy sites from three buoy site arrays (TAO, PIRATA and RAMA). The validation results demonstrate that the performance of LHF estimations derived from the ML methods (including ANN, RF, BR and RANSAC) were significantly better than individual LHF products, indicated by R2 increasing by 3.7–46.4%. Among them, the LHF estimation using the ANN method increased the R2 of the four-individual ocean LHF products (ranging from 0.56 to 0.79) to 0.88 and decreased the RMSE (ranging from 19.1 to 37.5) to 11 W m−2. Compared to three other ML methods (RF, BR and RANSAC), ANN method exhibited the best performance according to the validation results. The results of relative uncertainty analysis using the triangle cornered hat (TCH) method show that the ensemble LHF product using ML methods has lower relative uncertainty than individual LHF product in most area. The ANN was employed to implement the mapping of annual average ocean LHF over low-latitude at a spatial resolution of 0.25° during 2003–2007. The ocean LHF fusion products estimated from ANN methods were 10–30 W m−2 lower than those of the four original ocean products (MERRA-2, JOFURO-3, ERA-I and GSSTF-3) and were more similar to observations.

Leaf area index (LAI) is an essential vegetation parameter that represents the light energy utilization and vegetation canopy structure. As the only in-operation hyperspectral satellite launched by China, GF-5 is potentially useful for accurate LAI estimation. However, there is no research focus on evaluating GF-5 data for LAI estimation. Hyperspectral remote sensing data contains abundant information about the reflective characteristics of vegetation canopies, but these abound data also easily result in a dimensionality curse. Therefore, feature selection (FS) is necessary to reduce data redundancy to achieve more reliable estimations. Currently, machine learning (ML) algorithms have been widely used for FS. Moreover, the same ML algorithm is usually conducted for both FS and regression in LAI estimation. However, no evidence suggests that this is the optimal solution. Therefore, this study focuses on evaluating the capacity of GF-5 spectral reflectance for estimating LAI and the performances of different combination of FS and ML algorithms. Firstly, the PROSAIL model, which coupled leaf optical properties model PROSPECT and the scattering by arbitrarily inclined leaves (SAIL) model, was used to generate simulated GF-5 reflectance data under different vegetation and soil conditions, and then three FS methods, including random forest (RF), K-means clustering (K-means) and mean impact value (MIV), and three ML algorithms, including random forest regression (RFR), back propagation neural network (BPNN) and K-nearest neighbor (KNN) were used to develop nine LAI estimation models. The FS process was conducted twice using different strategies: Firstly, three FS methods were conducted to search the lowest dimension number, which maintained the estimation accuracy of all bands. Then, the sequential backward selection (SBS) method was used to eliminate the bands having minimal impact on LAI estimation accuracy. Finally, three best estimation models were selected and evaluated using reference LAI. The results showed that although the RF_RFR model (RF used for feature selection and RFR used for regression) achieved reliable LAI estimates (coefficient of determination (R2) = 0.828, root mean square error (RMSE) = 0.839), the poor performance (R2 = 0.763, RMSE = 0.987) of the MIV_BPNN model (MIV used for feature selection and BPNN used for regression) suggested using feature selection and regression conducted by the same ML algorithm could not always ensure an optimal estimation. Moreover, RF selection preserved the most informative bands for LAI estimation so that each ML regression method could achieve satisfactory estimation results. Finally, the results indicated that the RF_KNN model (RF used as feature selection and KNN used for regression) with seven GF-5 spectral band reflectance achieved the better estimation results than others when validated by simulated data (R2 = 0.834, RMSE = 0.824) and actual reference LAI (R2 = 0.659, RMSE = 0.697).

The methods for accurately fusing medium- and high-spatial-resolution satellite reflectance are vital for monitoring vegetation biomass, agricultural irrigation, ecological processes and climate change. However, the currently existing fusion methods cannot accurately capture the temporal variation in reflectance for heterogeneous landscapes. In this study, we proposed a new method, the spatial and temporal reflectance fusion method based on the unmixing theory and a fuzzy C-clustering model (FCMSTRFM), to generate Landsat-like time-series surface reflectance. Unlike other data fusion models, the FCMSTRFM improved the similarity of pixels grouped together by combining land cover maps and time-series data cluster algorithms to define endmembers. The proposed method was tested over a 2000 km2 study area in Heilongjiang Provence, China, in 2017 and 2018 using ten images. The results show that the accuracy of the FCMSTRFM is better than that of the popular enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) (correlation coefficient (R): 0.8413 vs. 0.7589; root mean square error (RMSE): 0.0267 vs. 0.0401) and the spatial-temporal data fusion approach (STDFA) (R: 0.8413 vs. 0.7666; RMSE: 0.0267 vs. 0.0307). Importantly, the FCMSTRFM was able to maintain the details of temporal variations in complicated landscapes. The proposed method provides an alternative method to monitor the dynamics of land surface variables over complicated heterogeneous regions.

Reliable estimates of terrestrial latent heat flux (LE) at high spatial and temporal resolutions are of vital importance for energy balance and water resource management. However, currently available LE products derived from satellite data generally have high revisit frequency or fine spatial resolution. In this study, we explored the feasibility of the high spatiotemporal resolution LE fusion framework to take advantage of the Moderate Resolution Imaging Spectroradiometer (MODIS) and Chinese GaoFen-1 Wide Field View (GF-1 WFV) data. In particular, three-fold fusion schemes based on Enhanced Spatial and Temporal Adaptive Reflectance Fusion Model (ESTARFM) were employed, including fusion of surface reflectance (Scheme 1), vegetation indices (Scheme 2) and high order LE products (Scheme 3). Our results showed that the fusion of vegetation indices and further computing LE (Scheme 2) achieved better accuracy and captured more detailed information of terrestrial LE, where the determination coefficient (R2) varies from 0.86 to 0.98, the root-mean-square error (RMSE) ranges from 1.25 to 9.77 W/m2 and the relative RSME (rRMSE) varies from 2% to 23%. The time series of merged LE in 2017 using the optimal Scheme 2 also showed a relatively good agreement with eddy covariance (EC) measurements and MODIS LE products. The fusion approach provides spatiotemporal continuous LE estimates and also reduces the uncertainties in LE estimation, with an increment in R2 by 0.06 and a decrease in RMSE by 23.4% on average. The proposed high spatiotemporal resolution LE estimation framework using multi-source data showed great promise in monitoring LE variation at field scale, and may have value in planning irrigation schemes and providing water management decisions over agroecosystems.

An accurate estimation of spatially and temporally continuous latent heat flux (LE) is essential in the assessment of surface water and energy balance. Various satellite-derived LE products have been generated to enhance the simulation of terrestrial LE, yet each individual LE product shows large discrepancies and uncertainties. Our study used Extremely Randomized Trees (ETR) to fuse five satellite-derived terrestrial LE products to reduce uncertainties from the individual products and improve terrestrial LE estimations over Europe. The validation results demonstrated that the estimation using the ETR fusion method increased the R2 of five individual LE products (ranging from 0.53 to 0.61) to 0.97 and decreased the RMSE (ranging from 26.37 to 33.17 W/m2) to 5.85 W/m2. Compared with three other machine learning fusion models, Gradient Boosting Regression Tree (GBRT), Random Forest (RF), and Gaussian Process Regression (GPR), ETR exhibited the best performance in terms of both training and validation accuracy. We also applied the ETR fusion method to implement the mapping of average annual terrestrial LE over Europe at a resolution of 0.05 ◦ in the period from 2002 to 2005. When compared with global LE products such as the Global Land Surface Satellite (GLASS) and the Moderate Resolution Imaging Spectroradiometer (MODIS), the fusion LE using ETR exhibited a relatively small gap, which confirmed that it is reasonable and reliable for the estimation of the terrestrial LE over Europe.

Downward shortwave radiation (RS) drives many processes related to atmosphere−surface interactions and has great influence on the earth’s climate system. However, ground-measured RS is still insufficient to represent the land surface, so it is still critical to generate high accuracy and spatially continuous RS data. This study tries to apply the random forest (RF) method to estimate the RS from the Himawari-8 Advanced Himawari Imager (AHI) data from February to May 2016 with a two-km spatial resolution and a one-day temporal resolution. The ground-measured RS at 86 stations of the Climate Data Center of the Chinese Meteorological Administration (CDC/CMA) are collected to evaluate the estimated RS data from the RF method. The evaluation results indicate that the RF method is capable of estimating the RS well at both the daily and monthly time scales. For the daily time scale, the evaluation results based on validation data show an overall R value of 0.92, a root mean square error (RMSE) value of 35.38 (18.40%) Wm−2, and a mean bias error (MBE) value of 0.01 (0.01%) Wm−2. For the estimated monthly RS, the overall R was 0.99, the RMSE was 7.74 (4.09%) Wm−2, and the MBE was 0.03 (0.02%) Wm−2 at the selected stations. The comparison between the estimated RS data over China and the Clouds and Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) RS dataset was also conducted in this study. The comparison results indicate that the RS estimates from the RF method have comparable accuracy with the CERES-EBAF RS data over China but provide higher spatial and temporal resolution.

Water use efficiency (WUE) measures the tradeoff between carbon uptake and water consumption in terrestrial ecosystems. It remains unclear how the responses of WUE to drought vary with drought severity. We assessed the spatio-temporal variations of ecosystem WUE and its responses to drought for terrestrial ecosystems in Southwest China over the period 2000−2017. The annual WUE values varied with vegetation type in the region: Forests (3.25 gC kg−1H2O) > shrublands (2.00 gC kg−1H2O) > croplands (1.76 gC kg−1H2O) > grasslands (1.04 gC kg−1H2O). During the period 2000−2017, frequent droughts occurred in Southwest China, and overall, drought had an enhancement effect on WUE. However, the effects of drought on WUE varied with vegetation type and drought severity. Croplands were the most sensitive to drought, and slight water deficiency led to the decline of cropland WUE. Over grasslands, mild drought increased its WUE while moderate and severe drought reduced its WUE. For forests and shrublands, mild and moderate drought increased their WUE, and only severe drought reduce their WUE, indicating that these ecosystems had stronger resistance to drought. Assessing the patterns and trends of ecosystem WUE and its responses to drought are essential for understanding plant water use strategy and informing ecosystem water management.

The temporal and spatial characteristics of meteorological drought have been investigated to provide a framework of methodologies for the analysis of drought in the Beijing-Tianjin-Hebei metropolitan area (BTHMA) in China. Using the Reconnaissance Drought Index (RDI) as an indicator of drought severity, the characteristics of droughts have been examined. The Beijing-Tianjin-Hebei metropolitan area was divided into 253 grid-cells of 27 × 27km and monthly precipitation data for the period of 1960–2010 from 33 meteorological stations were used for global interpolation of precipitation using spatial co-ordinate data. Drought severity was assessed from the estimated gridded RDI values at multiple time scales. Firstly, the temporal and spatial characteristics of droughts were analyzed, and then drought severity-areal extent-frequency (SAF) annual curves were developed. The analysis indicated that the frequency of moderate and severe droughts was about 9.10% in the BTHMA. Using the SAF curves, the return period of selected severe drought events was assessed. The identification of the temporal and spatial characteristics of droughts in the BTHMA will be useful for the development of a drought preparedness plan in the region.

Estimating cropland latent heat flux (LE) from continental to global scales is vital to modeling crop production and managing water resources. Over the past several decades, numerous LE models were developed, such as 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) and the modified satellite-based Priestley-Taylor LE algorithm (MS-PT). However, these LE models have not been directly compared over the global cropland ecosystem using various algorithms. In this study, we evaluated the performances of these four LE models using 34 eddy covariance (EC) sites. The results showed that mean annual LE for cropland varied from 33.49 to 58.97 W/m2 among the four models. The interannual LE slightly increased during 1982-2009 across the global cropland ecosystem. All models had acceptable performances with the coefficient of determination (R2) ranging from 0.4 to 0.7 and a root mean squared error (RMSE) of approximately 35 W/m2. MS-PT had good overall performance across the cropland ecosystem with the highest R2, lowest RMSE and a relatively low bias. The reduced performances of MOD16 and RRS, with R2 ranging from 0.4 to 0.6 and RMSEs from 30 to 39 W/m2, might be attributed to empirical parameters in the structure algorithms and calibrated coefficients.

Temporal-related features are important for improving land cover classification accuracy using remote sensing data. This study investigated the efficacy of phenological features extracted from time series MODIS Normalized Difference Vegetation Index (NDVI) data in improving the land cover classification accuracy of Landsat data. The MODIS NDVI data were first fused with Landsat data via the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) algorithm to obtain NDVI data at the Landsat spatial resolution. Next, phenological features, including the beginning and ending dates of the growing season, the length of the growing season, seasonal amplitude, and the maximum fitted NDVI value, were extracted from the fused time series NDVI data using the TIMESAT tool. The extracted data were integrated with the spectral data of the Landsat data to improve classification accuracy using a maximum likelihood classifier (MLC) and support vector machine (SVM) classifier. The results indicated that phenological features had a statistically significant effect on improving the land cover classification accuracy of single Landsat data (an approximately 3% increase in overall classification accuracy), especially for vegetation type discrimination. However, the phenological features did not improve on statistical measures including the maximum, the minimum, the mean, and the standard deviation values of the time series NDVI dataset, especially for human-managed vegetation types. Regarding different classifiers, SVM could achieve better classification accuracy than the traditional MLC classifier, but the improvement in accuracy obtained using advanced classifiers was inferior to that achieved by involving the temporally derived features for land cover classification.

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.

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/m2 to 87.6 W/m2, and with the square of correlation coefficient (R2) 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 R2 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.

As an important indicator to characterize the surface vegetation, fractional vegetation cover (FVC) with high spatio-temporal resolution is essential for earth surface process simulation. However, due to technical limitations and the influence of weather, it is difficult to generate temporally continuous FVC with high spatio-temporal resolution based on a single remote-sensing data source. Therefore, the objective of this study is to explore the feasibility of generating high spatio-temporal resolution FVC based on the fusion of GaoFen-1 Wide Field View (GF-1 WFV) data and Moderate-resolution Imaging Spectroradiometer (MODIS) data. Two fusion strategies were employed to identify a suitable fusion method: (i) fusing reflectance data from GF-1 WFV and MODIS firstly and then estimating FVC from the reflectance fusion result (strategy FC, Fusion_then_FVC). (ii) fusing the FVC estimated from GF-1 WFV and MODIS reflectance data directly (strategy CF, FVC_then_Fusion). The FVC generated using strategies FC and CF were evaluated based on FVC estimated from the real GF-1 WFV data and the field survey FVC, respectively. The results indicated that strategy CF achieved higher accuracies with less computational cost than those of strategy FC both in the comparisons with FVC estimated from the real GF-1 WFV (CF:R2 = 0.9580, RMSE = 0.0576; FC: R2 = 0.9345, RMSE = 0.0719) and the field survey FVC data (CF: R2 = 0.8138, RMSE = 0.0985; FC: R2 = 0.7173, RMSE = 0.1214). Strategy CF preserved spatial details more accurately than strategy FC and had a lower probability of generating abnormal values. It could be concluded that fusing GF-1 WFV and MODIS data for generating high spatio-temporal resolution FVC with good quality was feasible, and strategy CF was more suitable for generating FVC given its advantages in estimation accuracy and computational efficiency.

Fractional vegetation cover (FVC) is an important parameter for many environmental and ecological models. Large-scale and long-term FVC products are critical for various applications. Currently, several global-scale FVC products have been generated with remote sensing data, such as VGT bioGEOphysical product Version 2 (GEOV2), PROBA-V bioGEOphysical product Version 3 (GEOV3) and Global LAnd Surface Satellite (GLASS) FVC products. However, studies comparing and validating these global-scale FVC products are rare. Therefore, in this study, the performances of three global-scale time series FVC products, including the GEOV2, GEOV3, and GLASS FVC products, are investigated to assess their spatial and temporal consistencies. Furthermore, reference FVC data generated from high-spatial-resolution data are used to directly evaluate the accuracy of these FVC products. The results show that these three FVC products achieve general agreements in terms of spatiotemporal consistencies over most regions. In addition, the GLASS and GEOV2 FVC products have reliable spatial and temporal completeness, whereas the GEOV3 FVC product contains much missing data over high-latitude regions, especially during wintertime. Furthermore, the GEOV3 FVC product presents higher FVC values than GEOV2 and GLASS FVC products over the equator. The main differences between the GEOV2 and GLASS FVC products occur over deciduous forests, for which the GLASS product presents slightly higher FVC values than the GEOV2 product during wintertime. Finally, temporal profiles of the GEOV2 and GLASS FVC products show better consistency than the GEOV3 FVC product, and the GLASS FVC product presents more reliable accuracy (R2 = 0.7878, RMSE = 0.1212) compared with the GEOV2 (R2 = 0.5798, RMSE = 0.1921) and GEOV3 (R2 = 0.7744, RMSE = 0.2224) FVC products over these reference FVC data.

Many existing satellite evapotranspiration (ET)-based drought indices have characterized regional drought condition successfully, but the relatively short time span of ET products limits their use in long-term climatological drought assessment. In this study, we assess Evaporative Drought Index (EDI) as a drought monitoring indicator over Northeast China through a retrospective comparison with drought-related indicators. After verifying its utility for detecting documented regional drought events and impacts of drought on crop production, we apply it to improve our understanding of the variation in dryness over Northeast China from 1982 to 2015. Our results illustrate that EDI is generally effective for characterizing terrestrial moisture condition and its standardized formula, namely, Standardized Evaporative Drought Index (sEDI) corresponds well with historical drought events and inter-annual grain crop yields over Northeast China. Although the calculation of sEDI does not directly incorporate precipitation and soil moisture, statistical analyses indicate sEDI can detect drought in accordance with the Standardized Precipitation Index (SPI) and Palmer Drought Severity Index (PDSI), with the highest correlations found in the west part of Northeast China (R < −0.7). Further analysis illustrates sEDI is more related to commonly-used drought metrics over areas with short canopy vegetation (R < −0.5) than woodland (R < −0.2), which suggests precipitation may not be a good representative of drought condition over areas with deep-rooted vegetation. Then, we find 56.5% of Northeast China shows an upward dry trend from 1982 to 2015, which mainly concentrates in the west part of the study area. Conversely, 14.4% of Northeast China shows a significant wetted trend and most of them locate at cropland areas, due to the improved water management. This study suggests that EDI is a feasible method to monitor spatially distributed drought condition and can provide unique drought information not reflected by rainfall deficits, which also can be used to evaluate traditional precipitation-based indicators.