Improvement of evapotranspiration simulation study in the Hailar River basin under the influence of vegetation dynamics.

In arid and semi-arid areas with <400 mm of precipitation, evapotranspiration (ET) accounts for about 80% of precipitation and is the main water consumer in the watershed. However, vegetation greening in recent years will increase ET and exacerbate the aridity of the area by affecting soil moisture in the root system. Vegetation changes are regional and spatially heterogeneous, therefore, in order to characterize ET changes under vegetation dynamics, it is necessary to expand the spatial scale of ET simulation. However, widely used evapotranspiration simulation models, such as the Shuttleworth-Wallace model (SW model), are deficient in reflecting the direct and indirect effects of vertical (i.e., soil depths) and horizontal (i.e., vegetation dynamics) directions. Based on field sampling and constructed structural equation model (SEM), we found that vegetation dynamics affect evapotranspiration not only directly, but also indirectly by affecting soil moisture at different depths. On this basis, we defined the weighting coefficients of 0.85 and 0.15 for grassland vegetation zones, 0.3, 0.15, 0.20, 0.25, 0.10 for forest-grass interspersed zones, and 0.20, 0.55, 0.25 for forested zones, respectively, based on the SEM results. Different soil moisture weighting coefficients were defined within different vegetation type zones and the improved SW model is called S-W-α. Comparing the simulation results with the measured data, S-W-α improved the ET simulation accuracy in this region by 33.92% and the improved ET spatial trend can respond to the dynamic changes of vegetation. Replacing the ET module in the Block-wise use of TOPMODEL and Muskingum-Cunge method mode (BTOP model) with the modified S-W-α, the results show that the simulation accuracy of the improved model is increased by 25%, and the Nash is higher than 75% for both the rate period and the validation period, which realizes the extension of the model from the point scale to the basin scale. The modified model may provide technical support for simulation of evapotranspiration and management of ecosystem health in ecologically fragile areas.

Mapping Chinese annual gross primary productivity with eddy covariance measurements and machine learning.

Annual gross primary productivity (AGPP) is the basis for grain production and terrestrial carbon sequestration. Mapping regional AGPP from site measurements provides methodological support for analysing AGPP spatiotemporal variations thereby ensures regional food security and mitigates climate change. Based on 641 site-year eddy covariance measuring AGPP from China, we built an AGPP mapping scheme based on its formation and selected the optimal mapping way, which was conducted through analysing the predicting performances of divergent mapping tools, variable combinations, and mapping approaches in predicting observed AGPP variations. The reasonability of the selected optimal scheme was confirmed by assessing the consistency between its generating AGPP and previous products in spatiotemporal variations and total amount. Random forest regression tree explained 85 % of observed AGPP variations, outperforming other machine learning algorithms and classical statistical methods. Variable combinations containing climate, soil, and biological factors showed superior performance to other variable combinations. Mapping AGPP through predicting AGPP per leaf area (PAGPP) explained 86 % of AGPP variations, which was superior to other approaches. The optimal scheme was thus using a random forest regression tree, combining climate, soil, and biological variables, and predicting PAGPP. The optimal scheme generating AGPP of Chinese terrestrial ecosystems decreased from southeast to northwest, which was highly consistent with previous products. The interannual trend and interannual variation of our generating AGPP showed a decreasing trend from east to west and from southeast to northwest, respectively, which was consistent with data-oriented products. The mean total amount of generated AGPP was 7.03 ± 0.45 PgC yr-1 falling into the range of previous works. Considering the consistency between the generated AGPP and previous products, our optimal mapping way was suitable for mapping AGPP from site measurements. Our results provided a methodological support for mapping regional AGPP and other fluxes.