Vegetation plays a fundamental role within terrestrial ecosystems, serving as a cornerstone of their functionality. Presently, these crucial ecosystems face a myriad of threats, including deforestation, overgrazing, wildfires, and the impact of climate change. The implementation of remote sensing for monitoring the status and dynamics of vegetation ecosystems has emerged as an indispensable tool for advancing ecological research and effective resource management. This study takes a comprehensive approach by integrating ecosystem monitoring indicators and aligning them with the objectives of SDG15. We conducted a thorough analysis by leveraging global 500 m resolution products for vegetation Leaf Area Index (LAI) and land cover classification spanning the period from 2016 to 2020. This encompassed the calculation of annual average LAI, identification of anomalies, and evaluation of change rates, thereby enabling a comprehensive assessment of the global status and transformations occurring within major vegetation ecosystems. In 2020, a discernible rise in the annual Average LAI of major vegetation ecosystems on a global scale became evident when compared to data from 2016. Notably, the ecosystems demonstrating a slight increase in area constituted the largest proportion (34.23%), while those exhibiting a significant decrease were the least prevalent (6.09%). Within various regions, such as Eastern Europe, Central Africa, and South Asia, substantial increases in both forest ecosystem area and annual Average LAI were observed. Furthermore, Eastern Europe and Central America recorded significant expansions in both grassland ecosystem area and annual average LAI. Similarly, regions experiencing notable growth in both cropland ecosystem areas and annual average LAI encompassed Southern Africa, Northern Europe, and Eastern Africa.

Cloud plays a crucial role in surface downward longwave radiation (DLR). To determine optimal algorithms for cloud‐sky DLR calculation under various climates, three types of algorithms were assessed, (i) four empirical algorithms determining cloudy DLR by simple cloud correction using the cloud fraction, (ii) three parameterized algorithms determining the cloud contribution by cloud temperature, and (iii) a semiempirical algorithm, Zhou‐Cess, parameterized by the cloud water path. A sensitivity study was conducted using a Moderate Resolution Transmittance (MODTRAN) code to first determine the sensitive cloud parameters of DLR. Then, these algorithms were validated using synthetic data simulated by MODTRAN, ground‐observed data, and satellite‐observed data. When all the input parameters were accurate, the cloud‐correction algorithms showed poor performance. Cloud‐temperature‐based algorithms showed much better results but exhibited positive systematic biases. The Zhou‐Cess algorithm performed best but could not precisely describe the effect caused by cloud variations. When these algorithms were applied to ground‐ and satellite‐based data, the accuracies of DLR calculations were affected by the uncertainty in the atmospheric and cloud parameters. The simple empirical algorithms showed the poorest results. The cloud‐temperature‐based algorithms were greatly influenced by the uncertainty in cloud base temperature and cloud fraction and showed acceptable results when cloud fractions were accurate. The Zhou‐Cess algorithm revealed the best results at most sites and was less impacted by cloud parameter uncertainties; therefore, this algorithm is suggested for cloudy‐sky DLR calculation with poor‐quality cloud parameters.

Incoming surface solar irradiance (SSI) is essential for calculating Earth’s surface radiation budget and is a key parameter for terrestrial ecological modeling and climate change research. Remote sensing images from geostationary and polar-orbiting satellites provide an opportunity for SSI estimation through directly retrieving atmospheric and land-surface parameters. This paper presents a new scheme for estimating SSI from the visible and infrared channels of geostationary meteorological and polar-orbiting satellite data. Aerosol optical thickness and cloud microphysical parameters were retrieved from Geostationary Operational Environmental Satellite (GOES) system images by interpolating lookup tables of clear and cloudy skies, respectively. SSI was estimated using pre-calculated offline lookup tables with different atmospheric input data of clear and cloudy skies. The lookup tables were created via the comprehensive radiative transfer model, Santa Barbara Discrete Ordinate Radiative Transfer (SBDART), to balance computational efficiency and accuracy. The atmospheric attenuation effects considered in our approach were water vapor absorption and aerosol extinction for clear skies, while cloud parameters were the only atmospheric input for cloudy-sky SSI estimation. The approach was validated using one-year pyranometer measurements from seven stations in the SURFRAD (SURFace RADiation budget network). The results of the comparison for 2012 showed that the estimated SSI agreed with ground measurements with correlation coefficients of 0.94, 0.69, and 0.89 with a bias of 26.4 W/m², −5.9 W/m², and 14.9 W/m² for clear-sky, cloudy-sky, and all-sky conditions, respectively. The overall root mean square error (RMSE) of instantaneous SSI was 80.0 W/m² (16.8%), 127.6 W/m² (55.1%), and 99.5 W/m² (25.5%) for clear-sky, cloudy-sky (overcast sky and partly cloudy sky), and all-sky (clear-sky and cloudy-sky) conditions, respectively. A comparison with other state-of-the-art studies suggests that our proposed method can successfully estimate SSI with a maximum improvement of an RMSE of 24 W/m². The clear-sky SSI retrieval was sensitive to aerosol optical thickness, which was largely dependent on the diurnal surface reflectance accuracy. Uncertainty in the pre-defined horizontal visibility for ‘clearest sky’ will eventually lead to considerable SSI retrieval error. Compared to cloud effective radius, the retrieval error of cloud optical thickness was a primary factor that determined the SSI estimation accuracy for cloudy skies. Our proposed method can be used to estimate SSI for clear and one-layer cloud sky, but is not suitable for multi-layer clouds overlap conditions as a lower-level cloud cannot be detected by the optical sensor when a higher-level cloud has a higher optical thickness.

The quantitative characterization of the thermal conditions in the Tibetan Plateau has long been a focal point of global research. Downward shortwave radiation, as a crucial component, plays an important role in numerous land surface processes while also serving as a significant indicator of the plateau’s thermal state. In order to gain a more comprehensive understanding of the Earth’s radiation budget in the Tibetan Plateau region, this study undertook an evaluation of six radiation products (ISCCP-FH, CERES-SYN, GLASS DSR, Himawari-8, MCD18A1, and ERA5). Two sets of ground measurements (downward shortwave radiation values from 10 CMA sites and 6 sites provided by the National Tibetan Plateau Data Center) in 2015 and 2016 were used as validation data to verify the accuracy of the remote sensing products. The results show that in the Tibetan Plateau region, CERESC products show the highest accuracy among the six data products with a bias (relative bias) of −7.57 W/m² (3.46%), RMSE (relative RMSE) of 32.77 W/m² (14.99%), and coefficient of determination of 0.80. Among all products, only the ERA5 products overestimated the value of downward shortwave radiation in the Tibetan Plateau region with a bias (relative bias) of 15.62 W/m² (7.14%). By employing a spatial resolution upscaling approach, we assessed the influence of varying spatial resolutions on the validation accuracy, with the results indicating minimal impact. Through an analysis of the impact of cloud factors and aerosol factors on the validation accuracy, it is deduced that ERA5, Himawari-8, and MCD18A1 products are significantly influenced by cloud factors, whereas the CERES-SYN product is notably affected by aerosol factors.