Abstract
Regional water conservation is the key to maintain water supply. Based on the data of 2000, 2010 and 2020, this paper takes Jingyu County in the Changbai Mountain as the study area. By land use conversion matrix, geostatistical analysis, InVEST model and scenario analysis, this paper aims to find out the variation of water conservation and the its related influencing factors in Changbai Mountain region. From 2000 to 2020, forest and cultivated land decreased, while grassland and developed land increased. The water conservation increased from 2.94 × 108 m3 to 4.83 × 108 m3, with a total increase of 64.29%. Forest has the highest average water conservation, followed by grassland, cultivated land and developed land. Climatic factors were important to 91.92% of the study area, whereas 6.22% of it was influenced by land use. When forest decrease in some areas, the main control factor of 48.71% of the area shifted into land use, that could lead to a sharp decrease in water conservation. Climate determines the overall spatial pattern of water conservation; land use determines the local characteristics by controlling runoff. The increase of water conservation is beneficial to improve the reserves and development potential of mineral water.
HIGHLIGHTS
Application of InVEST model and scenario analysis in Changbai Mountain area.
Climate change is a dominant factor impacting water conservation.
Change of water conservation is more severe in area of land use conversion.
High water conservation capacity is conducive to mineral water exploitation.
Suggestions are put forward for mineral water resource management.
Graphical Abstract
INTRODUCTION
Water conservation is an important ecological service function, playing important role in improving the climate, regulating runoff, and maintaining balance and stability of the ecosystem (Sharafatmandrad & Mashizi 2021). It refers to the amount of water resources retained after evaporation and runoff and can be preserved in the ecosystem through plants, interception of litter and soil storage in a specific spatio-temporal range (Costanza et al. 2004). The amount of water conservation can reflect the capacity of water supply, soil and water retention and runoff regulation in the region (Xu et al. 2019).
The methods of precipitation storage, comprehensive storage capacity and water budget are widely used to evaluate water conservation. However, the precipitation storage method didn't distinguish the evapotranspiration by different vegetation (Hou et al. 2018), and the comprehensive water storage capacity method required a large demand for data (Zhang et al. 2012). While the water budget method utilized precipitation, evapotranspiration and runoff to estimate the water conservation by vegetation and soil with multiple advantages such as easy data acquisition and good adaptability to a Geographic Information System (GIS), thus applicable to wide range (Sun et al. 2018; Hu et al. 2021).
The InVEST model (Integrated Valuation of Ecosystem Services and trade-offs) developed by Stanford University, The Nature Conservancy (TNC) and World Wide Fund for Nature (WWF) (https://naturalcapitalproject.stanford.edu/software/invest), aims at providing a scientific basis for decision makers to weigh the benefits and impacts of human activities by simulating changes in the amount of ecosystem goods and values under different land cover scenarios (Sharp et al. 2020). The Water Yield module of InVEST model based on water budget method is widely used currently (Donohue et al. 2012; Gao et al. 2017). Scholars at home and abroad can use it to classify water conservation levels, discuss the characteristics of its intra-annual variation (Su & Fu 2013; Kim & Jung 2020), analyse the particularity of water conservation capacity in special climate and geomorphic areas (Lang et al. 2018; Yang et al. 2020), evaluate the importance of its influencing factors (Lang et al. 2017; Dai et al. 2020; Lian et al. 2020), and predict the future development trend (Shrestha et al. 2020). In addition, the model is applied to analyse the impact of cultivated conversion and afforestation on water conservation (Kim et al. 2017; Qi et al. 2019; Daneshi et al. 2021), explore the changes in water yield of the watershed by ecological management measures (such as the developed of protected areas) (Shoyama & Yamagata 2014; Mafios & Burney 2017). In addition, its combination with other models (such as CA-Markov) and parameters correction conducted recently makes it a more comprehensive approach. (Chu et al. 2018; Yang et al. 2019).
Among the above studies, forest ecosystem has always been the major focus of scholars’ attention due to its strong water conservation capacity important role in controlling the regional water cycle (Yu et al. 2017). When the changes of forest ecosystem happen, they will have great impact on the ecological environment, especially on the drinking water source, which may further affect its water supply function (Zhang et al. 2020). In China, balancing the protection and exploitation of water resources in forest ecosystems is a critical concern for policy makers. Changbai Mountain area is an important distribution area of forest and wetland ecosystem in China with rich mineral water resources. Jingyu County in this region is one of the three major mineral spring cities in the world, owning large spring flow and excellent water quality (Yan et al. 2015, 2016). At present, studies on water conservation in Changbai Mountain are mainly static analyses of annual water conservation, investigation on changes of different periods, and analysis involving dynamic law of water conservation capacity under climate or land use change scenarios are very limited. Therefore, in order to understand and optimize the further exploitation of mineral water resources in Changbai Mountain, it is important to clarify the mechanism of water conservation changes in the region first.
In view of the current research status and gaps, as well as urgent needs for future water resources and ecological protection, this paper calculated the water conservation of the study area and analysed its spatio-temporal change rules based on InVEST model. Then, three different climate and land use scenarios were created based on the data of 2000 and 2020, so as to study the impact of climatic and land use factors on water conservation capacity. Finally, the relationship between water conservation capacity and local mineral water resources was discussed.
METHODS
Study area
Jingyu is an important source of mineral water in China. A total of 84 springs have been discovered in this county with large amount of mineral water and good water quality (Bian et al. 2019; Wang et al. 2021). Some mineral water also contains strontium and other trace elements. The southwest part of study area is the concentrated distribution area of mineral water. The number of mineral springs in the area is large, and the water quality is good. The local famous mineral springs such as Qinglong Spring and Baijing spring are exposed here.
Research methods
Calculation of water yield

Calculation of water conservation
where, WRx is the amount of water conservation of grid x; Rsx is the surface runoff of grid x; α is runoff coefficient.
Scenario analysis
Data requirements
The data needed to run InVEST model are annual precipitation grid, annual reference evapotranspiration grid, land use type grid, soil maximum root depth grid, PAWC (Plants Available Water Content) grid, watershed vector map, biophysical series table and seasonal coefficient (Z coefficient) of Jingyu from 2000 to 2020. The data sources are shown in Table 1.
Data sources
Data type . | Data sources . | Accuracy . |
---|---|---|
Annual precipitation | China Meteorological Data Service Centre (http://data.cma.cn/) | 14 stations |
Annual reference evapotranspiration | ||
Land use | GlobeLand30 (http://globeland30.org/) | 30 m |
Watershed | Geospatial Data Cloud site, Chinese Academy of Sciences. (http://www.gscloud.cn) | 30 m |
Maximum root depth | National Tibetan Plateau Data Center (http://www.tpdc.ac.cn/) (Meng & Wang 2018) | 1 km |
PAWC |
Data type . | Data sources . | Accuracy . |
---|---|---|
Annual precipitation | China Meteorological Data Service Centre (http://data.cma.cn/) | 14 stations |
Annual reference evapotranspiration | ||
Land use | GlobeLand30 (http://globeland30.org/) | 30 m |
Watershed | Geospatial Data Cloud site, Chinese Academy of Sciences. (http://www.gscloud.cn) | 30 m |
Maximum root depth | National Tibetan Plateau Data Center (http://www.tpdc.ac.cn/) (Meng & Wang 2018) | 1 km |
PAWC |
In addition, based on previous research results on water conservation capacity of forest ecosystem in Northeast China and empirical data, biophysical table and runoff coefficient are determined as shown in Table 2 (Zhang et al. 2016; Zhang 2019).
Biophysical table and runoff coefficient
Land use code . | Land use type . | LULC_vega . | Root depth (mm) . | Evapotranspiration coefficient . | Runoff coefficient . |
---|---|---|---|---|---|
10 | Cultivated | 1 | 1,500 | 0.65 | 0.6 |
20 | Forest | 1 | 3,500 | 1 | 0.2 |
30 | Grassland | 1 | 2,000 | 0.65 | 0.35 |
50 | Wetlands | 0 | 1 | 1 | 0.8 |
60 | Water | 0 | 1 | 1.1 | 0.9 |
80 | Developed | 0 | 1 | 0.3 | 0.85 |
Land use code . | Land use type . | LULC_vega . | Root depth (mm) . | Evapotranspiration coefficient . | Runoff coefficient . |
---|---|---|---|---|---|
10 | Cultivated | 1 | 1,500 | 0.65 | 0.6 |
20 | Forest | 1 | 3,500 | 1 | 0.2 |
30 | Grassland | 1 | 2,000 | 0.65 | 0.35 |
50 | Wetlands | 0 | 1 | 1 | 0.8 |
60 | Water | 0 | 1 | 1.1 | 0.9 |
80 | Developed | 0 | 1 | 0.3 | 0.85 |
aLULC_veg Values must be 1 for vegetated land use except wetlands, and 0 for all other land uses.
ArcGIS processing diagram of input data required by InVEST (a. Precipitation; b. Evapotranspiration; c. Land use; d. Watershed; e. Maximum root depth; f. PAWC).
ArcGIS processing diagram of input data required by InVEST (a. Precipitation; b. Evapotranspiration; c. Land use; d. Watershed; e. Maximum root depth; f. PAWC).
Parameter calibration and model validation
RESULTS
Land use change
The area of land use types in the study area from large to small is forest, followed by cultivated land, grassland, developed land, water and wetlands. In 2020, the forest area was 2,460.21 km2, accounting for 79.75% of the total area. The arable land area was 453.65 km2, accounting for 14.71% of the total area. From 2000 to 2020, the transition between cultivated land and forest was significant. Based on the land use conversion matrix (Tables 3 and 4), from 2000 to 2010, the forest area decreased from 2,560.58 km2 to 2,400.79 km2, accounting for 5.18%, amongst which the conversion area from forest to cultivated land and grassland was 51.36 km2 and 100.6 km2 respectively. From 2010 to 2020, the forest area has increased, with a total area of 2,400.88 km2, and the area of cultivated land and grassland converted to forest was 40.9 km2 and 19.65 km2 respectively. In addition, cultivated land area decreased by 56.26 km2, with an abandonment rate of 11.03%, Although the forest area remains smaller than in 2000, the policy of returning farmland to forest has achieved good results.
Land use conversion matrix 2000–2010 (unit: km2)
2010\2000 . | Cultivated . | Developed . | Forest . | Grassland . | Water . | Wetlands . | Total . |
---|---|---|---|---|---|---|---|
Cultivated | 458.26 | 0.68 | 2.92 | 1.42 | 0.32 | 0 | 463.61 |
Developed | 0.05 | 17.03 | 0.01 | 0 | 0 | 0 | 17.09 |
Forest | 51.36 | 1.65 | 2,397.62 | 100.6 | 7.03 | 2.32 | 2,560.58 |
Water | 0.23 | 0.04 | 0.24 | 0.05 | 42.81 | 0.01 | 43.39 |
Wetlands | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Total | 509.90 | 19.41 | 2,400.79 | 102.07 | 50.18 | 2.33 | 3,084.67 |
2010\2000 . | Cultivated . | Developed . | Forest . | Grassland . | Water . | Wetlands . | Total . |
---|---|---|---|---|---|---|---|
Cultivated | 458.26 | 0.68 | 2.92 | 1.42 | 0.32 | 0 | 463.61 |
Developed | 0.05 | 17.03 | 0.01 | 0 | 0 | 0 | 17.09 |
Forest | 51.36 | 1.65 | 2,397.62 | 100.6 | 7.03 | 2.32 | 2,560.58 |
Water | 0.23 | 0.04 | 0.24 | 0.05 | 42.81 | 0.01 | 43.39 |
Wetlands | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Total | 509.90 | 19.41 | 2,400.79 | 102.07 | 50.18 | 2.33 | 3,084.67 |
Land use conversion matrix 2010–2020 (unit: km2)
2020\2010 . | Cultivated . | Developed . | Forest . | Grassland . | Water . | Wetlands . | Total . |
---|---|---|---|---|---|---|---|
Cultivated | 422.23 | 12.97 | 63.91 | 7.64 | 3.12 | 0 | 509.88 |
Developed | 1.26 | 17.69 | 0.12 | 0.25 | 0.10 | 0 | 19.41 |
Forest | 23.01 | 1.15 | 2,349.02 | 24.38 | 3.25 | 0.07 | 2,400.88 |
Grassland | 5.79 | 0.56 | 44.03 | 51.32 | 0.35 | 0.02 | 102.08 |
Water | 1.32 | 0.25 | 2.93 | 2.00 | 43.61 | 0.02 | 50.13 |
Wetlands | 0.00 | 0 | 0.09 | 0.00 | 0.02 | 2.21 | 2.33 |
Total | 453.62 | 32.62 | 2,460.10 | 85.59 | 50.45 | 2.33 | 3,084.71 |
2020\2010 . | Cultivated . | Developed . | Forest . | Grassland . | Water . | Wetlands . | Total . |
---|---|---|---|---|---|---|---|
Cultivated | 422.23 | 12.97 | 63.91 | 7.64 | 3.12 | 0 | 509.88 |
Developed | 1.26 | 17.69 | 0.12 | 0.25 | 0.10 | 0 | 19.41 |
Forest | 23.01 | 1.15 | 2,349.02 | 24.38 | 3.25 | 0.07 | 2,400.88 |
Grassland | 5.79 | 0.56 | 44.03 | 51.32 | 0.35 | 0.02 | 102.08 |
Water | 1.32 | 0.25 | 2.93 | 2.00 | 43.61 | 0.02 | 50.13 |
Wetlands | 0.00 | 0 | 0.09 | 0.00 | 0.02 | 2.21 | 2.33 |
Total | 453.62 | 32.62 | 2,460.10 | 85.59 | 50.45 | 2.33 | 3,084.71 |
The developed land is primarily distributed in the central area of Jingyu Town. During the 20 years, the area increased from 17.09 to 32.62 km2, accounting for 90.81%, indicating that the villages and towns were developed greatly, and the development was the most rapid in the last decade. The area of cultivated and forest converted to developed land was 12.45 km2 and 2.98 km2 respectively, accounting for 80.21% and 19.22% of the increment respectively. In the study area, the land use changes indicate a decrease in forest and cultivated land and an increase in grassland and developed land.
Results of water yield and water conservation
Temporal distribution of water conservation
Distribution map of water yield and water conservation (a. Water yield; b. Water conservation; c. Water conservation in the nature reserve).
Distribution map of water yield and water conservation (a. Water yield; b. Water conservation; c. Water conservation in the nature reserve).
Data of the concentrated distribution area of mineral water were extracted from the distribution of water conservation, and the results are shown in Figure 6(c). The average water conservation in the distribution area of mineral water in the three periods was 99.79 mm, 167.91 mm and 176.37 mm, and the perennial average was 148.02 mm, which is higher than the average of the study area in the same period. Thus, this value is conducive to the supply of mineral water.
Spatial distribution of water conservation
Statistical chart of water conservation and contribution rate of each town in study area. (a. the average water conservation capacity of each township; b. Contribution rate of each township to total water conservation.).
Statistical chart of water conservation and contribution rate of each town in study area. (a. the average water conservation capacity of each township; b. Contribution rate of each township to total water conservation.).
Influencing factors of water conservation
Climatic factors
The average precipitation was 719.82 mm, 799.05 mm and 830.81 mm in 2000, 2010 and 2020 respectively, showing an increasing trend, whereas the water yield and water conservation also increased. Comparing the distribution maps of precipitation, water yield and water conservation in the three stages, the overall distribution pattern of water yield and water conservation in the whole study area is the same as that of precipitation. Meanwhile, SPSS was used to calculate the correlation amongst water yield, precipitation and PET in three stages and further analyse the influence mechanism of meteorological factors. A positive correlation was observed between water yield and precipitation, with a correlation coefficient of 0.883. By contrast, a negative correlation was found between water production and PET, with a correlation coefficient of −0.770.
Land use factors
Statistical chart of water yield and water conservation of different land use types.
Statistical chart of water yield and water conservation of different land use types.
DISCUSSION
Response of water conservation to climate and land use change
Influence mechanism of climatic and land use factors
With regard to climatic factors, compared with the precipitation map of 2000, 2010 and 2020, the distribution law of precipitation is not evident. This finding is due to the small size of the study area and the large randomness of annual precipitation distribution in different areas, which leads to the different spatial distribution of precipitation in different years and the inconsistent spatial distribution of annual water yield and water conservation in 2000, 2010 and 2020. However, the spatial distribution of water yield and water conservation is consistent with that of precipitation in the same year. Therefore, climatic factors determine the overall spatial and temporal distribution pattern of water and water conservation in the study area.
For land use factors, the water yield of cultivated land and developed land is high, but the water conservation is low. This finding is due to the low vegetation coverage and root depth of these two land use types, resulting in less precipitation intercepted and utilised by vegetation, and less transpiration loss than that of forest and grassland. Thus, most of the precipitation is converted into water yield. In addition, the infiltration condition of cultivated land and developed land is poor, and the ability of vegetation to retain water is low. Therefore, the runoff coefficient is relatively large, and the water yield is completely lost through the surface runoff out of the region, resulting in low water conservation. However, for forests, the canopy and litter are broad and thick, respectively, which result in an evident interception effect on runoff and large amount of water conserved (Julian & Gardner 2014).
For wetlands and water area, this study considers that precipitation falling into the water body will directly form runoff out of the study area, which cannot be used in determining the statistics of water conservation in the study area. Given the scattered distribution of cultivated land and developed land, local high-value or low value areas appear in the distribution map of water yield and water conservation. Therefore, land use affects water conservation by influencing runoff, which determines the local characteristics of water conservation in the study area.
Relative importance analysis
Distribution diagram of relative importance of factors and changes of water conservation.
Distribution diagram of relative importance of factors and changes of water conservation.
Statistical chart of water conservation under different factors from 2010 to 2020.
Statistical chart of water conservation under different factors from 2010 to 2020.
As shown in Figure 10, if the climatic effect is not considered, when the forest changes to cultivated, developed land and grassland, then the water conservation will decrease. The reduction is the largest when forest changes to cultivated land (172 mm), and it is the least when forest changes to grassland (49 mm). Most of the crops in northeast China, such as rice and corn, are relatively low. These crops conserve water and soil only by themselves and their shallow roots; but forest can conserve water through canopy, understory shrubs, humus and roots. Thus, the conversion of forest to cultivated leads to a sharp decline in water conservation. Although Section 3.3.2 considers that developed land has a worse water conservation capacity than cultivated land, the area of forest converted into cultivated land from 2010 to 2020 is 23.01 km2, which is larger than the area converted into construction land (1.15 km2). Therefore, the loss of water conservation caused by the conversion of forest to farmland is the largest. In addition, the Changbai Mountain is a national key construction tourism area. The newly built towns in this region are small with a large internal green area, and some areas retain the characteristics of forest or grassland. When cultivated land and developed land are converted to forest and grassland, respectively, water conservation will increase, and the increment is the largest when the cultivated land changes to forest (149 mm). When the climatic effect is considered, water conservation under different categories (Figure 10) fluctuates because of differences in precipitation and evapotranspiration in different areas, with variation ranging from −19 to 34 mm, which is significantly smaller than that under the influence of land use only −172 to 149 mm. Therefore, although climatic factors are the main factors affecting water conservation in the study area, large-scale land use change, particularly deforestation and reclamation caused by human activities, will change the amount of runoff loss. Such changes in water conservation are often more dramatic than those caused by climate. At present, the expansion of cultivated land and developed land remains the main inhibiting factor of water conservation in the area. Therefore, the policy of returning farmland to forest should be implemented for the sustainable exploitation of mineral water resources.
The effect of water conservation on mineral water resources
The Changbai Mountain is a world-famous mineral water source. It is of reference significance for other water sources in the world to clarify the influence of water conservation on mineral water resources and establish the relationship between them. Based on previous studies, mineral water in the Changbai Mountain comes from a single source and it is primarily supplied by atmospheric precipitation (Bian et al. 2019; Li et al. 2022). Precipitation initially enters into the aquifer through soil and then into the vadose zone, where it is fully exchanged with basalt to form high-quality mineral water. The mineral water in this area has strong renewability, and it is sensitive to changes in atmospheric precipitation. Water conservation in the InVEST model is directly or indirectly generated by atmospheric precipitation, and it is the residual part of water yield through runoff loss. In this paper, the mineral water concentrated distribution area is densely forested, the ecosystem is well protected, and the runoff loss is less, which is conducive to the storage of water resources. The water conservation of the mineral water concentrated distribution area is higher than the average level of the study area, which is conducive to sufficient water–rock reaction between mineral water and basalt in the basin. Consequently, the supply and storage of mineral water will increase, and the water quantity and quality of mineral water, suitable for the development and utilisation of mineral water.
Given its location in the middle temperate zone of China, the Changbai Mountain has a high altitude (average altitude of 1,640–1,826 m), and the snow and ice cover period last for 9 months every year. The local forest ecosystem can effectively hold the meltwater, allowing it to recharge mineral water slowly. The recharge time of snow meltwater for runoff in the Changbai Mountain is 1 month, which indicates also reflects the strong water conservation capacity of forest ecosystem (Feng 2021). However, in areas with low vegetation coverage, ice and snow meltwater often flows out of the ecosystem in ice flood, resulting in regional water resource loss. Therefore, the forest ecosystem can effectively ensure not only high water conservation in the Changbai Mountain area, but also mineral water supply.
Suggestions on the management of local mineral water resources
This study is can clarify the causes of gains and losses of water conservation services, which can be further extended to ecosystem protection policies in related areas. At the end of the past century, large-scale deforestation and land reclamation activities in China severely affected the ecological environment, leading to decline of river flow in many areas, reduction of forest area and weakening of ecological services. In solving the abovementioned problems, the Chinese government has vigorously implemented the ecological policy of returning farmland to forests, grasslands and lakes, which has achieved remarkable effects (Tang et al. 2020). However, with the rapid economic development and large-scale urban expansion from the 21st century, the contradiction between construction land and natural area (such as forest and grassland) has to be put on the urgent agenda. Hence, policy makers must identify the contribution of various land uses to ecological service functions and identify areas that are most in need of protection or improvement.
In this paper, the change map of water conservation can directly reflect the gains and losses of water conservation services under the influence of climate and land use, and scenario analysis provides the separate effect of two influencing factors on water conservation. Areas with reduced water conservation should be improved, that is, areas with reduced forest and grassland in areas where land use change is the main factor. Considering the actual content of the study area, most of these areas are urban and road construction land in recent years. Therefore, the relationship between urban expansion and ecological protection should be weighed.
The formulation of the ecological policy should not only consider the overall changes of water conservation in the region, but also reflect the particularity of different functional areas to make the policy more localised. The central part of the research area is the urban centre, and economic development is the main goal. Therefore, the development of construction land and cultivated land can be carried out without affecting the ecology, thereby ensuring the normal operation of the ecosystem. However, for mineral water protection areas in the southwest of the study area, farming and animal husbandry activities should be prohibited. Furthermore, mineral water sources and surrounding rivers should be protected, particularly the upper reaches of rivers, to prevent the establishment of mineral water plant and related human activities on the water conservation function caused by negative effects.
CONCLUSIONS
This paper makes a comprehensive analysis of the change mechanism of water conservation in the Changbai mountain area and evaluated the influence of water conservation on mineral water resources. The results show that from 2000 to 2020, the changes of land use in the study area decreased in forest and cultivated land and increased in grassland and developed land. Water conservation increased from 2.94 × 108m3 to 4.83 × 108m3, total conservation increased by 64.29%, and the proportion of water production into water conservation was 40%. Under natural conditions, climatic factors played a decisive role in water conservation in 91.92% of the study regions. However, when forest decreased, cultivated and construction land increase occurred in some areas, the main control factor of 48.71% of the area shifted into land use factor, and these behaviours will cause a sharp decrease in water conservation. The change of land use factors may greatly reduce water conservation, or even directly go to 0, which is more severe than the impact of climatic factors.
This study showed that climatic factors determine the overall spatial pattern of water conservation, while land use factors determine the local characteristics of water conservation by controlling runoff. In addition, forest ecosystems can effectively intercept precipitation and store snow, so that water resources can be well stored. For mineral water source, this is an important supply source, conducive to the sustainable exploitation of high-quality mineral water. Therefore, in the process of mineral water mining and economic development, the conflict between urban construction, arable land expansion, and forest protection should be weighed to avoid forest destruction. This is conducive to the steady and healthy development of the local economy dominated by mineral water. The methods and results used in this study will be helpful to better understand the change mechanism of water conservation and reveal the significance of water conservation function to mineral water resources, playing a guiding role in decision making of local policy development and ecological protection.
ACKNOWLEDGEMENTS
The study was supported by the National Key Research and Development Program of China (grant number 2019YFC0409103) and the Key Projects of Jilin Provincial Department of Science and Technology (grant number 20190303076SF). We sincerely appreciate the funding of the above projects and the work of each author.
FUNDING
The study was supported by:
- 1.
National Key Research and Development Program of China (grant number 2019YFC0409103);
- 2.
Key Projects of Jilin Provincial Department of Science and Technology (grant number 20190303076SF).
CREDIT AUTHOR STATEMENT
Sun Wenhao: Conceptualization; Formal Analysis; Methodology; Investigation; Writing – Original Draft. Bian Jianmin: Funding Acquisition; Methodology; Writing – Review & Editing. Li Yihan: Resources; Data Curation. Li Jialin: Investigation; Resources.
DECLARATION OF INTERESTS
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
AVAILABILITY OF DATA AND MATERIALS
All relevant data are available from an online repository or repositories. The URL or references of all data sources have been noted in the relevant positions in the text.
DATA AVAILABILITY STATEMENT
All relevant data are available from an online repository or repositories. China Meteorological Data Service Centre (http://data.cma.cn/) GlobeLand30 (http://globeland30.org/) Geospatial Data Cloud site, Chinese Academy of Sciences. (http://www.gscloud.cn) National Tibetan Plateau Data Center (http://www.tpdc.ac.cn/).
CONFLICT OF INTEREST
The authors declare there is no conflict.