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.

Using high-resolution satellite data to perform routine (i.e., daily to weekly) monitoring of surface evapotranspiration, evapotranspiration (ET) (or LE, i.e., latent heat flux) has not been feasible because of the low frequency of satellite coverage over regions of interest (i.e., approximately every two weeks). Cloud cover further reduces the number of useable observations, and the utility of these data for routine ET or LE monitoring is limited. Moderate-resolution satellite imagery is available multiple times per day; however, the spatial resolution of these data is too coarse to enable the estimation of ET from individual agricultural fields or variations in ET or LE. The objective of this study is to combine high-resolution satellite data collected in the visible and near-infrared (VNIR) bands with data from the MODIS thermal-infrared (TIR) bands to estimate subpixel surface LE. Two temperature-sharpening methods, the disaggregation procedure for radiometric surface temperature (DisTrad) and the geographically-weighted regression (GWR)-based downscaling algorithm, were used to obtain accurate subpixel land surface temperature (LST) within the Zhangye oasis in China, where the surface is heterogeneous. The downscaled LSTs were validated using observations collected during the HiWATER-MUSOEXE (Multi-Scale Observation Experiment on Evapotranspiration) project. In addition, a remote sensing-based energy balance model was used to compare subpixel MODIS LST-based turbulent heat fluxes estimates with those obtained using the two LST downscaling approaches. The footprint validation results showed that the direct use of the MODIS LST approach does not consider LST heterogeneity at all, leading to significant errors (i.e., the root mean square error is 73.15 W·m⁻²) in LE, whereas the errors in the LE estimates obtained using DisTrad and GWR were 45.84 W·m⁻² and 47.38 W·m⁻², respectively. Furthermore, additional analysis showed that the ability of DisTrad and GWR to capture subpixel LST variations depends on the value of Shannon’s diversity index (SHDI) and the surface type within the flux contribution source area.

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.

Surface downward longwave radiation (DLR) is a crucial component in Earth’s surface energy balance. Yu et al. (2013) developed a parameterization for retrieving clear-sky DLR at high spatial resolution by combined use of satellite thermal infrared (TIR) data and column integrated water vapor (IWV). We extended the Yu2013 parameterization to Moderate Resolution Imaging Spectroradiometer (MODIS) data based on atmospheric radiative simulation, and we modified the parameterization to decrease the systematic negative biases at large IWVs. The new parameterization improved DLR accuracy by 1.9 to 3.1 W/m² for IWV ≥3 cm compared to the Yu2013 algorithm. We also compared the new parameterization with four algorithms, including two based on Top-of-Atmosphere (TOA) radiance and two using near-surface meteorological parameters and water vapor. The algorithms were first evaluated using simulated data and then applied to MODIS data and validated using surface measurements at 14 stations around the globe. The results suggest that the new parameterization outperforms the TOA-radiance based algorithms in the regions where ground temperature is substantially different (enough that the difference between them is as large as 20 K) from skin air temperature. The parameterization also works well at high elevations where atmospheric parameter-based algorithms often have large biases. Furthermore, comparing different sources of atmospheric input data, we found that using the parameters interpolated from atmospheric reanalysis data improved the DLR estimation by 7.8 W/m² for the new parameterization and 19.1 W/m² for other algorithms at high-altitude sites, as compared to MODIS atmospheric products.