Satellite-based leaf area index (LAI) products, such as the MODIS LAI, play an essential role in land surface and climate modeling research, from regional to global scales. However, data gaps and high-level noise can exist, thus limiting their applications to a broader scope. Our previous work has reprocessed the MODIS LAI Collection 5 (C5) product, and the reprocessed data have been widely used these years. In this study, the MODIS C6.1 LAI data were reprocessed to broaden its application as a successor. We updated the integrated two-step method that is used for MODIS C5 LAI and implemented it into the MODIS C6.1 LAI product. Comprehensive evaluations for the original and reprocessed products were conducted. The results showed that the reprocessed LAI data had better performance in validation against reference maps. In addition, the site scale time series of reprocessed data was much smoother and more consistent with adjacent values. The global scale comparison showed that, though the MODIS C6.1 LAI does have improvements in ground validation with LAI reference maps, its spatial continuity, temporal continuity, and consistency showed little improvement when compared to C5. In contrast, the reprocessed data were more spatiotemporally continuous and consistent. Based on this evaluation, some suggestions for using various MODIS LAI products were given. This study assessed the quality of these different versions of MODIS LAI products and demonstrated the improvement of the reprocessed C6.1 data, which we recommended for use as a substitute for the reprocessed C5 data in land surface and climate modeling.

Urban trees are recognized for mitigating urban thermal stress, therefore incorporating their effects is crucial for urban climate research. However, due to the limitation of remote sensing, the LAI in urban areas is generally masked (e.g., MODIS), which in turn limits its application in Urban Canopy Models (UCMs). To address this gap, we developed a high-resolution (500 m) and long-time-series (2000–2022) urban tree LAI dataset derived through the Random Forest model trained with MODIS LAI data, with the help of meteorological variables and tree height datasets. The results show that our dataset has high accuracy when validated against site reference maps, with R of 0.85 and RMSE of 1.03 m²/m². Compared to reprocessed MODIS LAI, our modeled LAI exhibits an RMSE ranging from 0.36 to 0.64 m²/m² and an R ranging from 0.89 to 0.97 globally. This dataset provides a reasonable representation of urban tree LAI in terms of magnitude and seasonal changes, thereby potentially enhancing its applications in UCMs and urban climate studies.

Transpiration represents more than 30% of the global land–atmosphere water exchange but is highly uncertain. Plant hydraulics was ignored in traditional land surface modeling, but recently plant hydraulics has been found to play an essential role in transpiration simulation. A new physical-based representation of plant hydraulic schemes (PHS) was recently developed and implemented in the Common Land Model (CoLM). However, it is unclear to what extent PHS can reduce these uncertainties. Here, we evaluated the PHS against measurements obtained at 81 FLUXNET sites. The transpiration of each site was estimated using an empirical evapotranspiration partitioning approach. The metric scores defined by the International Land Model Benchmarking Project (ILAMB) were used to evaluate the model performance and compare it with that of the CoLM default scheme (soil moisture stress (SMS)). The bias score of transpiration in PHS was higher than SMS for most sites, and more significant improvements were found in semi-arid and arid sites where transpiration was limited by soil moisture. The hydraulic redistribution in PHS optimized the soil water supply and thus improved the transpiration estimates. In humid sites, no significant improvement in seasonal or interannual variability of transpiration was simulated by PHS, which can be explained by the insensitivity of transpiration demand coupled to the photosynthesis response to precipitation. In arid and semi-arid sites, seasonal or interannual variability of transpiration was better captured by PHS than SMS, which was interpreted by the improved drought sensitivity for transpiration. Arid land is widespread and is expected to expand due to climate change, thus there is an urgent need to couple PHS in land surface models.

Land surface models (LSMs) require reliable forcing, validation, and surface attribute data as the foundation for effective model development and improvement. Eddy covariance flux tower data are widely regarded as the benchmark for LSMs. However, currently available flux tower datasets often require multiple aspects of processing to ensure data quality before application to LSMs. More importantly, these datasets frequently lack site-observed attribute data, such as fractional vegetation cover and leaf area index, which limits their utility as benchmarking data. Here, we conducted a comprehensive quality screening of the existing reprocessed flux tower dataset, including the proportion of gap-filled data, energy balance closure (EBC), and external disturbances such as irrigation and deforestation, leading to 90 high-quality sites. For these sites, we collected vegetation, soil, and topography data as well as wind speed, temperature, and humidity measurement heights from literature; regional networks; and Biological, Ancillary, Disturbance, and Metadata (BADM) files. We then compiled the final flux tower attribute dataset by filling in missing attributes with global data and classifying plant functional types (PFTs). This dataset is provided in NetCDF (Network Common Data Form) format with necessary descriptions and reference sources. Model simulations revealed substantial disparities in the output between the attribute data observed at the site and those commonly used by LSMs, underscoring the critical role of site-observed attribute data and increasing the emphasis on flux tower attribute data in the LSM community. The dataset addresses the lack of the site attribute to some extent, reduces uncertainty in LSM data source, and aids in diagnosing parameter and process deficiencies. The dataset is available at https://doi.org/10.5281/zenodo.12596218 (Shi et al., 2024).

Leaf optical properties (LOPs, i.e., leaf reflectance and transmittance), as a fundamental property of vegetation, are a key parameter in the canopy radiative transfer process. LOPs have a direct impact on the surface solar radiation partition and further affect surface flux exchanges. Recent works have provided reliable LOP data and mentioned that notable differences exist between the prescribed LOP values in current land surface models and measured LOP values, especially in the near-infrared (NIR) band. To evaluate the effects of different LOP values in land surface modeling, we ran two land surface models (the Community Land Model and the Common Land Model) with their default prescribed and measured values to examine the differences in simulated surface radiation partitions and fluxes. Our analyses show that differences in LOP values can lead to a large discrepancy in albedo, radiation partition, sensible heat flux and net radiation simulations. By using the measured LOP values, in the boreal forest zone, Southeast China, and the eastern United States, both models have a significantly increased surface albedo in the NIR band, with the difference exceeding 10% during JJA. Thus, the measured LOP values can improve the negative albedo bias in the boreal forest zone during summertime. Moreover, both models simulate less net radiation with a maximum reduction of 11 W/m2 when incorporating the measured LOP values. Therefore, the total sensible heat flux can be reduced by as much as 11 W/m2. The results of this study emphasize that different LOP values can have a considerable effect on the surface radiation budget and sensible heat flux simulations which need attention in land surface model development. However, in current offline simulations, the measured LOP values cause slight changes in land surface temperatures and gross primary productivity (GPP).

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 R² (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.