Global LAnd Surface Satellite Research Articles

Analyzing Spatiotemporal Variations and Driving Factors of Grassland in the Arid Region of Northwest China Surrounding the Tianshan Mountains

The Tianshan Mountains, the largest arid mountain range in Central Asia, feature diverse terrains and significant landscape heterogeneity. The grasslands within the Xinjiang Tianshan region are particularly sensitive to climate change and human activities. However, until recently, the patterns and mechanisms underlying grassland changes in this region have been unclear. In this study, we analyzed spatial and temporal changes in grassland fractional vegetation cover (FVC) from 2001 to 2020, analyzed spatial and temporal changes in grassland, and predicted future trends using Global Land Surface Satellite (GLASS) FVC data, trend analysis, and the Hurst index method. We also explored the driving mechanisms behind these changes through the structural equation model (SEM). The results showed that from 2001 to 2020, the grassland FVC in the Tianshan region of Xinjiang was higher in the central and western regions and lower in the northern and southern regions, showing an overall fluctuating growth trend, with a change in the growth rate of 0. 0017/a (p < 0.05), and that this change was spatially heterogeneous, with the sum of significant improvement (20.6%) and slight improvement (29.9%) being much larger than the sum of significant degradation (0.6%) and slight degradation (9.5%). However, the Hurst index (H = 0.47) suggests that this trend may not continue, and there is a risk of degradation. Our study uncovers the complex interactions between the Tianshan barrier effect and grassland ecosystems, highlighting regional differences in driving mechanisms. Although the impacts of climatic conditions in grasslands vary over time in different regions, the topography and its resulting hydrothermal conditions are still dominant, and the extent of the impact is susceptible to fluctuations of varying degrees due to extreme climatic events. Additionally, the number of livestock changes significantly affects the grasslands on the southern slopes of the Tianshan Mountains, while the effects of nighttime light are minimal. By focusing on the topographical barrier effect, this study enhances our understanding of grassland vegetation dynamics in the Tianshan Mountains of Xinjiang, contributing to improved ecosystem management strategies under climate change.

Evaluation of seven satellite-based and two reanalysis global terrestrial evapotranspiration products

Although comprehensive evaluation of different types of global terrestrial evapotranspiration (ET) products has been conducted, the satellite remote sensing techniques have prompted the development of several available global ET products, warranting a reassessment as products continue to evolve. Recently, we produced the long-term Global LAnd Surface Satellite (GLASS) ET product, but there is a lack of comparison and evaluation with other ET products and EC observations on a global scale. In this study, we evaluated the accuracy and uncertainty of seven satellite-based (GLASS-AVHRR, GLASS-MODIS, BESS, FLUXCOM, GLEAM, MOD16, and PML_V2) and two reanalysis (ERA5 and MERRA2) global terrestrial ET products at multiple scales for selecting the most suitable ET products and developing large-scale ET models. At the point scale, their accuracy was evaluated through direct comparison with in situ observations from 230 global flux towers. The results indicate that no single ET product can provide the most accurate ET estimates for all land cover types, although GLASS-MODIS [coefficient of determination (R2) of 0.51, Kling–Gupta efficiency (KGE) of 0.68] and FLUXCOM [R2 of 0.51, KGE of 0.66] outperform the seven other products [0.54 =< KGE =< 0.66, 0.35 =< R2 =< 0.50] at all sites. At the basin scale, the accuracy of ET products was assessed through 36 large river basins. The R2 values between all ET products and the water balance-derived ET (WBET) are >0.88, while the accuracies of the nine ET products differ in some sense. The three-cornered hat (TCH) method and comparison analysis are applied to assess the uncertainty of nine ET products at the pixel level. The TCH outputs reveal that GLASS-AVHRR and GLASS-MODIS are the two products with the lowest relative uncertainty, while MERRA2 has the largest relative uncertainty. Our results imply that there is no single ET product performing best in all respects. The selection of ET products for scientific research should consider their performance differences in spatial scale as well as the influence of land cover and climate conditions.

Benefit of incorporating GLASS remote sensing vegetation products in improving Noah-MP land surface temperature simulations on the Tibetan Plateau

Land Surface Temperature (LST) is important for diagnosing surface energy balance in land surface models (LSMs). However, LST simulation in current LSMs tends to show large cold biases, partially due to the reason that the model's prescribed vegetation parameters (e.g., Leaf Area Index (LAI) and Fraction of Vegetation Cover (FVC)) are misrepresented, especially in regions with complex topography and climate such as Tibetan Plateau. Recent advancements in remote sensing technologies provide a unique opportunity to improve the model's vegetation parameters at large scales. In this study, we practice two experiments to improve LST simulations in Noah-MP LSM by (1) incorporating LAI and FVC from the Global Land Surface Satellite (GLASS) remote sensing product (exp_RS); and (2) incorporating an empirical LAI and FVC parameterization scheme based on the soil temperature stress factor (exp_RL02). Results show that the effect of vegetation on simulated LST is the most significant in summer season when the model-satellite LAI and FVC differences are the largest. Compared to the default experiment that uses static LAI and FVC values from the model's look-up table (exp_CTL), the results in exp_RS and exp_RL02 show domain-wide improvement of the simulated LST. The LAI and FVC effect on LST are also well reflected in model's energy budget components (i.e., longwave emissivity, sensible and latent heat fluxes, etc). Validation of the model simulated soil temperature with in-situ observations further demonstrate the model improvements. Our study underscores the important role of vegetation in regulating surface energy transfer processes. Our study also highlights the feasibility and benefit of incorporating remote sensing data in improving land surface model simulations.

Simultaneous estimation of five temporally regular land variables at seven spatial resolutions from seven satellite data using a multi-scale and multi-depth convolutional neural network

Various satellite sensors have provided a huge amount of observations of Earth's environment at variable spatial and temporal resolutions. Many global coarse-resolution land products have been generated from single-satellite data, but global temporally regular land products at fine spatial resolutions (say 10-30 m) are scarce because of infrequent observations. An ideal inversion framework can estimate multiple global continuous land variables at different spatial resolutions by combining all sources of satellite data. This paper proposes a new framework that can estimate five land variables simultaneously from the top-of-atmosphere (TOA) reflectance acquired by seven satellite sensors based on a multi-scale and multi-depth convolutional neural network (MSDCNN). This framework enables us to estimate temporally regular land variables at as fine as 10 m spatial resolution by transforming information from satellite data at coarser spatial resolutions. The estimated land variables include Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), shortwave albedo, visible albedo, and spectral reflectance. Seven sensors that acquire data at different spatial resolutions, include Visible Infrared Imaging Radiometer Suite (VIIRS) (750 m), Moderate Resolution Imaging Spectroradiometer (MODIS) (500 m), Fengyun-3 (FY-3) Medium Resolution Spectral Imager (MERSI) (250 m), China-Brazil Earth Resources Satellite 04 (CBERS-04) Wide Field Imager (WFI) (73 m), Landsat 8 Operational Land Imager (OLI) (30 m), Gaofen-1 (GF-1) Wide-Field-of-View (16 m) and Sentinel 2 A/B Multispectral Imager (MSI) (10 m). This framework is mainly composed of four steps. First, a Shuffled Complex Evolution (SCE) optimization method is adopted to estimate these variables from VIIRS TOA reference. Second, a joint-output random forest regression (RF) method is used to link the satellite observations and the estimated values from step 1. Third, the six-type multi-resolution satellite observations are used to downscale the VIIRS TOA reflectance to six different spatial resolutions by using the MSDCNN. Finally, the downscaled six-resolution VIIRS TOA reflectance is fed into the multiple-variable RF model to estimate land variables. The framework was validated, and the results based on high-resolution reference maps from ImagineS network and time series shortwave albedo field values from Surface Radiation (SURFRAD) and Integrated Carbon Observation System (ICOS) network show that the retrieved variables had high validation accuracy, with root mean square error (RMSE) ranges of 0.361–0.489 (LAI), 0.023–0.120 (FAPAR), 0.013–0.026 (snow-free shortwave albedo), respectively. Comparison results of the retrieved multi-scale variables with the Sentinel 2 A/B (10 m), Landsat 8 (30 m), Global LAnd Surface Satellite (GLASS, 250 m, 500 m), MODIS (500 m), and VIIRS (750 m) products show that their values were close, with RMSE ranges of 0.107–0.273 (LAI), 0.015–0.027 (FAPAR), 0.003–0.007 (shortwave albedo), 0.001–0.007 (visible albedo), and 0.003–0.025 (spectral reflectance), respectively. The results of the direct validation as well as the product intercomparison show that this novel framework has the potential to be used in estimating global land variables at various spatial scales using a variety of satellite data sources.

Spatio-temporal variation in leaf area index in the Yan Mountains over the past 40 years and its relationship to hydrothermal conditions

Changes in hydrothermal conditions have significant effects on vegetation, but there is still a lack of understanding of how vegetation responds to land surface (surface temperature and soil moisture) and meteorological (temperature and precipitation) conditions in mountain regions. This study examined the trends of leaf area index (LAI) in the Yan Mountains over the last four decades using Global Land Surface Satellite (GLASS) data. The results showed a persistent increase of LAI (greening) over 20 % to 80 % of the study area in growing season, spring, summer and autumn. Anthropogenic activities caused the greening trend by crop management before 2000 and afforestation after 2000. The increasing rate of LAI varied with elevation, and the most significant increase occurred in areas between 300 and 900 m, and the lowest increase occurred in areas below 300 m. Moreover, we found that LAI was negatively correlated with land surface temperature and soil moisture, but positively correlated with precipitation and air temperature. The time-lag effect was found between hydrothermal factors and LAI in the past four decades. There was a time lag of 2–3 months between LAI changes and temperature/precipitation during the early and late stages of the growing season, and a time lag of 0–1 month during the middle stage. Specifically, there was no time lag in vegetation response to surface soil moisture, and a time lag of 2–3 months in vegetation response to land surface temperature from July to October. Our findings provide insights into how vegetation adapts to land surface and climatic hydrothermal conditions in mountain regions and can be used by governments to develop policies for ecological protection.

Improved global 250 m 8-day NDVI and EVI products from 2000–2021 using the LSTM model

Satellite vegetation index (VI) products, such as normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI), have been widely used. However, they are severely contaminated by clouds and other factors and provide false signals of the surface vegetation conditions. In this study, the new global seamless 250 m, eight-day NDVI and EVI products from 2000–2021 were developed from Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data using a long short-term memory (LSTM) neural network method. High-quality globally representative time series VI samples were constructed to train the model using a combination of the Savitzky-Golay filter (SG), Global LAnd Surface Satellite (GLASS) leaf area index (LAI) fitting and upper envelope methods. To evaluate the proposed method and the 250 m VI products, the MODIS VI product (MOD13Q1) was used for the inter-comparisons using four widely used VI reconstruction methods. Assuming that the MODIS VI data of high quality represents the true values, the root mean square error (RMSE) for NDVI and EVI generated by the LSTM model are 0.0734 and 0.0509, respectively.

Climate impact greater on vegetation NPP but human enhance benefits after the Grain for Green Program in Loess Plateau

Net primary productivity (NPP) has been the subject of numerous studies, but the qualitative and comparative study on NPP that considers human activities, climate change, and their interactions is still insufficient. The Grain for Green Program (GTGP) was established in the Loess Plateau since 1999, could be used as a reference study. The temporal and spatial patterns, change trend, and driving factors of vegetation NPP were analysed in Loess Plateau using data derived from GLASS NPP (the Global Land Surface Satellite Product), meteorology, atmospheric quality, hydrology and land use by correlation and improved residual trend analysis from 2000 to 2020. After GTGP implementation, the climate and human contributions to vegetation NPP were 62 % and 38 %. NPP increased at a rate of 6.3 g C m−2 a−1 from 2000 to 2020, which was higher in 2010–2020 (6.9 g C m−2 a−1) than it was in 2000–2010 (5.7 g C m−2 a−1). The same trend was followed by the contribution of climate change to NPP between 2010 and 2020 (70.6 %) and 2000–2010 (65 %). The GTGP core area was mostly where the vegetation NPP rising zone was located. From the original bare land in the desert area, 4442 km2 of vegetation grew and increased the NPP (68.7 g C m−2). Under the warm–humid climatic trend, the climate factors contributed to NPP more in grasslands (40 %) compared to woodlands (29 %) and precipitation (33 %) was more sensitive than temperature (12.7 %). Consequently, sediment transport, PM2.5 and dusty weather decreased significantly. In conclusion, although climate change has a greater impact on vegetation, human activities have enhanced the impact faster in the first decade, which benefit and profit from the environment finally.

Spatiotemporally consistent global dataset of the GIMMS leaf area index (GIMMS LAI4g) from 1982 to 2020

Abstract. Leaf area index (LAI) with an explicit biophysical meaning is a critical variable to characterize terrestrial ecosystems. Long-term global datasets of LAI have served as fundamental data support for monitoring vegetation dynamics and exploring its interactions with other Earth components. However, current LAI products face several limitations associated with spatiotemporal consistency. In this study, we employed the back propagation neural network (BPNN) and a data consolidation method to generate a new version of the half-month 1/12∘ Global Inventory Modeling and Mapping Studies (GIMMS) LAI product, i.e., GIMMS LAI4g, for the period 1982–2020. The significance of the GIMMS LAI4g was the use of the latest PKU GIMMS normalized difference vegetation index (NDVI) product and 3.6 million high-quality global Landsat LAI samples to remove the effects of satellite orbital drift and sensor degradation and to develop spatiotemporally consistent BPNN models. The results showed that the GIMMS LAI4g exhibited overall higher accuracy and lower underestimation than its predecessor (GIMMS LAI3g) and two mainstream LAI products (Global LAnd Surface Satellite (GLASS) LAI and Long-term Global Mapping (GLOBMAP) LAI) using field LAI measurements and Landsat LAI samples. Its validation against Landsat LAI samples revealed an R2 of 0.96, root mean square error of 0.32 m2 m−2, mean absolute error of 0.16 m2 m−2, and mean absolute percentage error of 13.6 % which meets the accuracy target proposed by the Global Climate Observation System. It outperformed other LAI products for most vegetation biomes in a majority area of the land. It efficiently eliminated the effects of satellite orbital drift and sensor degradation and presented a better temporal consistency before and after the year 2000. The consolidation with the reprocessed MODIS LAI allows the GIMMS LAI4g to extend the temporal coverage from 2015 to a recent period (2020), producing the LAI trend that maintains high consistency before and after 2000 and aligns with the reprocessed MODIS LAI trend during the MODIS era. The GIMMS LAI4g product could potentially facilitate mitigating the disagreements between studies of the long-term global vegetation changes and could also benefit the model development in earth and environmental sciences. The GIMMS LAI4g product is open access and available under Attribution 4.0 International at https://doi.org/10.5281/zenodo.7649107 (Cao et al., 2023).

Estimating daily surface downward shortwave radiation over rugged terrain without bright surface at 30 m on clear-sky days using CERES data

ABSTRACTIn this study, the authors propose a model, called the Daily Downward Shortwave Radiation Random Forest Model over Rugged Terrain (DSRMT), to accurately calculate the downward shortwave radiation over a terrain without bright surface on clear days at a daily scale (DSRdaily−rugged). It was built by using the random forest method based on the comprehensive samples from CERES4_SYN1deg_Ed4A within 17 typical mountainous regions. DSRMT could directly estimate DSRdaily-rugged from the instantaneous direct and diffuse solar radiation on a flat surface during 10:30–14:30hrs on each day by comparing with the terrain factors from a digital elevation model, broadband albedo from the Global Land Surface Satellite, and ancillary information. The in-situ validation results showed that it generally delivered superior performance in estimating DSRdaily-rugged at any time during 10:30–14:30hrs, especially at noon, yielding a validated root mean-squared error (RMSE) of 24.90–29.22 Wm−2 and mean absolute error (MAE) of 19.16–22.94 Wm−2, and the average weighted DSRdaily-rugged were usually more accurate with the RMSE and MAE of 21.63 and 17.14 Wm−2. Overall, DSRMT was found to deliver satisfactory performance because of its high accuracy, robustness, ease of implementation, and efficiency, so it has the strong potential to be widely used in practice.

Assessing the Reliability of Global Carbon Flux Dataset Compared to Existing Datasets and Their Spatiotemporal Characteristics

Land carbon fluxes play a critical role in ecosystems, and acquiring a comprehensive global database of carbon fluxes is essential for understanding the Earth’s carbon cycle. The primary methods of obtaining the spatial distribution of land carbon fluxes include utilizing machine learning models based on in situ measurements, estimating through satellite remote sensing, and simulating ecosystem models. Recently, an innovative machine learning product known as the Global Carbon Flux Dataset (GCFD) has been released. In this study, we assessed the reliability of the GCFD by comparing it with existing data products, including two machine learning products (FLUXCOM and NIES (National Institute for Environmental Studies)), two ecosystem model products (TRENDY and EC-LUE (eddy covariance–light use efficiency model)), and one remote sensing product (Global Land Surface Satellite), on both site and global scales. Our findings indicate that, in terms of average absolute difference, the spatial distribution of the GCFD is most similar to the NIES product, albeit with slightly larger discrepancies compared to the other two types of products. When using site observations as the benchmark, gross primary production (GPP), respiration of ecosystem (RECO), and net ecosystem exchange of machine learning products exhibit higher R2 (ranging from 0.57 to 0.85, 0.53–0.79, and 0.31–0.70, respectively) compared to model products and remote sensing products. Furthermore, we analyzed the spatial and temporal distribution characteristics of carbon fluxes in various regions. The results demonstrate an upward trend in both GPP and RECO over the past two decades, while NEE exhibits an opposite trend. This trend is particularly pronounced in tropical regions, where higher GPP is observed in tropical, subtropical, and oceanic climate zones. Additionally, two remote sensing variables that influence changes in carbon fluxes, i.e., fraction absorbed photosynthetically active radiation and leaf area index, exhibit relatively consistent spatial and temporal characteristics. Overall, our study can provide valuable insights into different types of carbon flux products and contribute to understanding the general features of global carbon fluxes.

Exploring the contribution of environmental factors to evapotranspiration dynamics in the Three-River-Source region, China

Evapotranspiration (ET) exerts substantial effects on terrestrial ecosystem carbon and water cycles. A deep understanding of the spatial and temporal variability in ET and its dominant controls has crucial implications for future ecological management, especially for the ecologically-fragile Three-River-Source region (TRSR) in China. The study examined the performance of three ET products including Global Land Evaporation Amsterdam Model (GLEAM), Famine Early Warning Systems Network Land Data Assimilation System (FLDAS) and Global Land Surface Satellite (GLASS) based on eddy covariance-based flux observations. On this basis, ET dynamics in the TRSR from 1982 to 2018 were analyzed at two time scales using the GLASS ET product. Spatially, the pattern in ET showed a gradual upward trend from northwest to southeast. Meanwhile, ET in the majority of regions significantly enhanced, with multiannual mean ET and summer ET (June to August) increasing by 0.609 mm/yr and 0.765 mm/yr, respectively. Especially, ET in summer dominated the interannual dynamics of ET in the TRSR. The individual contributions of five environmental factors including Leaf Area Index (LAI), solar radiation (Rs), air temperature (Ta), precipitation (PRE) and vapour pressure deficit (VPD), to the variability in ET were isolated using ridge regression. Significant differences appeared when explaining ET over different time scales. Rs (33.05 %) and LAI (25.29 %) were the dominant control variables in shaping the interannual ET, with VPD contributing the least by 7.19 %. However, the long-term dynamics in summer ET were mainly controlled by LAI (35.73 %) and PRE (31.19 %), and Ta exerted the least effect by 5.73 %. All analyses will strengthen our understanding of the relationships between climate, vegetation and ET in the TRSR area, China.

Significant discrepancies of land surface daily net radiation among ten remotely sensed and reanalysis products

ABSTRACT 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.

Dense flux observations reveal the incapability of evapotranspiration products to capture the heterogeneity of evapotranspiration

Accurate terrestrial evapotranspiration (ET) estimation over complex surfaces with spatial heterogeneity is crucial for local, regional, and global applications. Currently, a variety of global ET products with different spatiotemporal resolutions have been developed and evaluated. However, little is known about their performance in capturing the heterogeneity of ET over complex surfaces. Focusing on the Heihe River Basin (HRB), a typical arid and semi-arid region that is hydrologically vulnerable, this study compared ET from eleven global ET products (six remotely- sensed ET, two land surface model ET, one hydrological model-based ET, one reanalysis ET, and one synthesis ET) and one regional ET dataset against dense eddy covariance observations in terms of the magnitude, seasonal cycle, and spatial pattern. In general, the remotely-sensed ET and synthesis ET outperformed other categories, with the operational Simplified Surface Energy Balance (SSEBop), Penman-Monteith–Leuning (PML), Moderate Resolution Imaging Spectroradiometer (MOD16) and Global LAnd Surface Satellite (GLASS) performing relatively better (root mean square error ranging from 1.22 to 1.57 mm d-1). Across all the land cover types, ET products reproduced a relatively feasible ET over the desert steppe, meadow, and barren, but substantially underestimated ET over ecosystems with high ET values but a low land area fraction over HRB (cropland, wetland, and forest). This highlights the importance of accurately representing the heterogeneity and local climates in the complex land surfaces for ET quantification in regional water resource management.

Responses of ecosystem water use efficiency to drought in the Lancang–Mekong River Basin

The ecosystem water use efficiency (WUE), a crucial indicator of how climate change will affect terrestrial ecosystems, depicts the coupling of the carbon gain and water loss in terrestrial ecosystems. In this study, the spatiotemporal variations in the WUE and its responses to drought in the Lancang–Mekong River Basin (LMRB) from 1982 to 2018 were investigated using the gross primary productivity (GPP) and evapotranspiration (ET) data acquired from the Global Land Surface Satellite (GLASS) products. The analyses revealed that: (1) the mean yearly WUE for the LMRB was 1.63 g C kg−1 H2O, with comparatively higher values in forests and warm temperate climatic types. The interaction of temperature and leaf area index was the main factor affecting the spatial distribution of WUE. The yearly WUE for the entire region exhibited a decreasing trend with a rate of −0.0009 g C kg −1 H2O·yr−1, and the spatially significantly decreasing area accounted for 41.67% of the total area. (2) The annual WUE was positively correlated with drought in the humid regions, accounting for 66.55% of the total area, while a negative relationship mainly occurred in the high-altitude cold region. (3) The ecosystem WUE lagged behind the drought by 3 months in most regions. The lag effect was more apparent in the grassland-dominated upstream region and the cropland-dominated Mekong Delta. (4) The resilience analysis revealed that the ecosystems in forests and temperate climate types were strongly resistant to drought, while the grassland and high-altitude regions with a dry and cold climate had relatively poor resilience. The results of this study shed light on how the WUE responds to drought across diverse land use types, climate types, and elevation gradients, uncovering fresh insights into the potential mechanisms behind the impact of drought on water and carbon cycles within ecosystems.

Generation of global 1 km daily soil moisture product from 2000 to 2020 using ensemble learning

Abstract. Motivated by the lack of long-term global soil moisture products with both high spatial and temporal resolutions, a global 1 km daily spatiotemporally continuous soil moisture product (GLASS SM) was generated from 2000 to 2020 using an ensemble learning model (eXtreme Gradient Boosting – XGBoost). The model was developed by integrating multiple datasets, including albedo, land surface temperature, and leaf area index products from the Global Land Surface Satellite (GLASS) product suite, as well as the European reanalysis (ERA5-Land) soil moisture product, in situ soil moisture dataset from the International Soil Moisture Network (ISMN), and auxiliary datasets (Multi-Error-Removed Improved-Terrain (MERIT) DEM and Global gridded soil information (SoilGrids)). Given the relatively large-scale differences between point-scale in situ measurements and other datasets, the triple collocation (TC) method was adopted to select the representative soil moisture stations and their measurements for creating the training samples. To fully evaluate the model performance, three validation strategies were explored: random, site independent, and year independent. Results showed that although the XGBoost model achieved the highest accuracy on the random test samples, it was clearly a result of model overfitting. Meanwhile, training the model with representative stations selected by the TC method could considerably improve its performance for site- or year-independent test samples. The overall validation accuracy of the model trained using representative stations on the site-independent test samples, which was least likely to be overfitted, was a correlation coefficient (R) of 0.715 and root mean square error (RMSE) of 0.079 m3 m−3. Moreover, compared to the model developed without station filtering, the validation accuracies of the model trained with representative stations improved significantly for most stations, with the median R and unbiased RMSE (ubRMSE) of the model for each station increasing from 0.64 to 0.74 and decreasing from 0.055 to 0.052 m3 m−3, respectively. Further validation of the GLASS SM product across four independent soil moisture networks revealed its ability to capture the temporal dynamics of measured soil moisture (R=0.69–0.89; ubRMSE = 0.033–0.048 m3 m−3). Lastly, the intercomparison between the GLASS SM product and two global microwave soil moisture datasets – the 1 km Soil Moisture Active Passive/Sentinel-1 L2 Radiometer/Radar soil moisture product and the European Space Agency Climate Change Initiative combined soil moisture product at 0.25∘ – indicated that the derived product maintained a more complete spatial coverage and exhibited high spatiotemporal consistency with those two soil moisture products. The annual average GLASS SM dataset from 2000 to 2020 can be freely downloaded from https://doi.org/10.5281/zenodo.7172664 (Zhang et al., 2022a), and the complete product at daily scale is available at http://glass.umd.edu/soil_moisture/ (last access: 12 May 2023).

Impacts of land cover changes on biogenic emission and its contribution to ozone and secondary organic aerosol in China

Abstract. The greening impacts on China from 2000 to 2017 led to an increase in vegetated areas and thus enhanced biogenic volatile organic compound (BVOC) emissions. BVOCs are regarded as important precursors for ozone (O3) and secondary organic aerosol (SOA). As a result, accurate estimation of BVOC emissions is critical to understand their impacts on air quality. In this study, the Model of Emissions of Gases and Aerosols from Nature (MEGAN) v2.1 was used to investigate the impact of different leaf area index (LAI) and land cover (LC) datasets on BVOC emissions in China in 2016, and the effects on O3 and SOA were evaluated based on the Community Multiscale Air Quality (CMAQ) modeling system. Three LAI satellite datasets of the Global LAnd Surface Satellite (GLASS), the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD15A2H version 6 (MOD15), and the Copernicus Global Land Service (CGLS), as well as three LC satellite datasets of the MODIS MCD12Q1 LC products, the Copernicus Climate Change Service (C3S) LC products, and the CGLS LC products, were used in five parallel experiments (cases: C1–C5). Results show that changing LAI and LC datasets of the model input has an impact on BVOC estimations. BVOC emissions in China ranged from 25.42 to 37.39 Tg in 2016 and were mainly concentrated in central and southeastern China. Changing the LC inputs for the MEGAN model has a more significant difference in BVOC estimates than using different LAI datasets. The combination of C3S LC and GLASS LAI performs better in the CMAQ model, indicating that it is the better choice for BVOC estimations in China. The highest contribution of BVOCs to O3 and SOA can reach 12 ppb and 9.8 µg m−3, respectively. Changing the MEGAN inputs further impacts the concentrations of O3 and SOA, especially changing LC datasets. The relative difference between MCD12Q1 LC and C3S LC is over 52 % and 140 % in O3 and biogenic SOA (BSOA) in central and eastern China. The BSOA difference is mainly attributed to the isoprene SOA (ISOA), a major contributor to BSOA. The relative differences in ISOA between different cases are up to 160 % in eastern China. Therefore, our results suggest that the uncertainties in MEGAN inputs should be fully considered in future O3 and SOA simulations.

Shortwave radiation balance modulates potential evapotranspiration over China

ABSTRACT Shortwave radiation is an influential driver of global hydrological cycle, as its variation will alter evapotranspiration (ET). While climate change discussion extensively examined ET response to temperature, wind speed and land use/cover variations, little is known about the contribution of shortwave radiation balance on ET partly because of data availability. In this study, we applied the newly-released Global LAnd Surface Satellite (GLASS) products to detect changes in shortwave radiation and albedo for 2003–2018 over China, and quantified their contributions on potential evapotranspiration (PET) which is the upper limit of ET and generally the basis to estimate it. The results showed that downward shortwave radiation declined, together with a slight decrease in albedo, causing a decrease of net shortwave radiation. While climate change affected PET, the decreased net shortwave radiation reduced PET during the period. We identified a hotspot area in the north China plain that is one of the main agricultural production regions, as the radiation balance induced a contrasting effect on PET for the warm and the cold seasons. Our study indicates that it is important to consider shortwave radiation balance when estimating PET and ET to make reliable evaluation of climate change on hydrological cycle and ecosystem response.

Triple collocation-based merging of multi-source gridded evapotranspiration data in the Nordic Region

Accurate evapotranspiration (ET) data are required for many hydro-meteorological applications. Compared with the traditional evaluation that requires in-situ measurements, the triple collocation (TC) technique estimates geophysical product errors without the need for ground truth, which is especially suitable over large areas lacking a dense in-situ network. However, violations of the zero-error cross-correlation (ECC) assumption are found to be the dominant sources of impairing the TC robustness. This study presents the first application of a TC-based merging framework that optimally considers ECC to merge multi-source gridded ET products in the Nordic Region during 2003–2018. The ECC estimates of each ET dataset pair calculated by the quadruple collocation approach are used to select the qualified triplets from four products, including FLUXCOM, Global Land Surface Satellite (GLASS), Global Land Evaporation and Amsterdam Model (GLEAM), and Penman-Monteith-Leuning Version 2 (PML-V2). Then the ET merged datasets are generated by weighting TC-based rescaled error variances of the parent datasets through least square merging. Finally, the accuracy of both the parent and the merged datasets are assessed with the Integrated Carbon Observation System (ICOS) flux data in the Nordic Region based on multiple statistical metrics. Results demonstrate that the ECC values provide intuitive evidence for filtering unqualified TC triplets. Both the absolute and relative error variances (signal-to-noise ratio) are considered for ET dataset evaluation. Overall PML-V2 has the best performance among the evaluated four products. Two merged ET datasets with the reference climatology of FLUXCOM outperform all parent products with the lowest errors by using ICOS data as reference among all sites – indicating the feasibility of TC technique for improving ET accuracy in the Nordic Region. This study also analyses the impacts of reference climatology selection on the TC merged results.

Climate Change and Anthropogenic Activity Co-Driven Vegetation Coverage Increase in the Three-North Shelter Forest Region of China

The Three-North Shelter Forest (TNSF) program is a significant ecological safety barrier in northern China, where both climate change and anthropogenic activity contribute to the increase in vegetation coverage observed. However, comprehensive effects of these factors on vegetation have not been accurately quantified yet. This study utilized the Global Land Surface Satellite (GLASS) Advanced Very-High-Resolution Radiometer (AVHRR) Fractional Vegetation Cover (FVC) data, meteorological data, and spatial distribution of ecological engineering to analyze spatiotemporal variation of FVC and climate in the TNSF program region in China during the period 1982–2018. A partial correlation analysis and residual analysis were performed to determine the relative contribution of climate change and anthropogenic activity to the FVC and the overall effect of ecological governance. Results showed that since 1982, the average FVC in the TNSF program region was 0.201–0.253, with an average growth rate of 0.01·(10a)−1. The FVC showed a significant increase in 66.45% of the TNSF region, and will continue to increase, while only 7.02% showed a significant decrease. The coefficient of variation showed a large spatial variation, with 30.86% being in very low stability regions, mainly distributed in Inner Mongolia and the Loess Plateau. A warm and wet climate is more conducive to increasing the FVC than the warm and dry climate, and ecological engineering has the largest impact on areas with an annual accumulated precipitation greater than 300 mm. A quantitative analysis revealed that climate change and anthropogenic activity contributed to the significant increase in the FVC in 15.58% and 46.81% of the TNSF region, respectively. Therefore, ecological governance projects, such as the TNSF program, play a crucial role in enhancing the FVC in this region.

The Hi-GLASS all-wave daily net radiation product: Algorithm and product validation

The surface net radiation (Rn) represents the balance of the radiative budget on the land surface and drives many physical and biological processes. An accurate and long-term product for global daily coverage of Rn at a high spatial resolution is needed for a variety of applications at regional and local scales. This study proposes two algorithms, called the downward shortwave radiation (DSR)-based algorithm and the top-of-atmosphere (TOA)-based algorithm, to estimate Rn by using Landsat data. The DSR-based algorithm consists of three conditional models, and was developed based on the analysis of the relationship between Rn and shortwave radiation as well as ancillary information from ground measurements and various datasets. The TOA-based algorithm was developed by linking Rn to TOA observations from Landsat sensors and ancillary information. The two algorithms were developed by using the random forest method. The results of their validation against ground measurements showed that the DSR-based algorithm outperformed the TOA-based algorithm in terms of accuracy, with a determination coefficient (R2) of 0.93, root-mean-squared error (RMSE) of 17.58 Wm−2, and bias of −4.27 Wm−2. It was stable under various conditions. We then applied the DSR-based algorithm to generate a product of the global daily Rn, called the High-resolution (Hi)- Global LAnd Surface Satellite (GLASS), from 2013 to 2018 at a spatial resolution of 30 m under a clear sky based on remotely sensed products, including the DSR from GLASS, the normalized difference vegetation index (NDVI) obtained from Landsat, surface broadband albedo from Hi-GLASS, and meteorological factors based on reanalysis data from MERRA2. Following its validation using in-situ observations from 2013 to 2018, the overall accuracy of the daily Rn acquired by Hi-GLASS under clear sky was found to be satisfactory, with a value of R2 of 0.90 and an RMSE of 25.03 Wm−2. Moreover, compared with the daily Rn obtained from the GLASS product at a spatial resolution of 5 km, that obtained by Hi-GLASS can better characterize the surface by providing more details and capturing the variations in the measurements, especially large and small values. However, due to limitations of the available datasets and the algorithm, the data on Rn for most regions lacked information on cloudy skies and areas at high latitudes. This information thus cannot be provided by Hi-GLASS yet. Moreover, the influence of the topography on values of Rn has not been thoroughly considered. Nonetheless, values of Rn under clear sky obtained from Hi-GLASS offer promise for use in a wide range of areas, and efforts are underway to improve this product.