Method and system for evaluating contribution rate of driving factor to vegetation carbon sequestration change

Abstract

The invention provides a method and system for evaluating the contribution rate of a driving factor to vegetation carbon sequestration change, and relates to the field of data processing, and the method comprises the steps: calculating the vegetation carbon sequestration data of a to-be-evaluated region in a to-be-evaluated time period according to the vegetation net primary productivity data of the to-be-evaluated region in the to-be-evaluated time period; obtaining driving factor data of the to-be-evaluated area in the to-be-evaluated time period, and calculating a change trend of a plurality of driving factors according to the driving factor data of the to-be-evaluated area in the to-be-evaluated time period; according to the vegetation carbon sequestration data and the driving factor data of the to-be-evaluated region in the to-be-evaluated time period, calculating correlation between the plurality of driving factors and vegetation carbon sequestration; and based on the change trend of the plurality of driving factors and the correlation between the plurality of driving factors and the vegetation carbon sequestration, calculating the contribution rates of the plurality of driving factors to the vegetation carbon sequestration change of the to-be-evaluated area in the to-be-evaluated time period, and having the advantage of evaluating the contribution rates of different driving factors to the vegetation carbon sequestration change.

Description

Technical Field

[0001] The present invention relates to the field of data processing, and in particular to a method and system for evaluating the contribution rate of driving factors to changes in vegetation carbon fixation. Background Art

[0002] With the rapid development of social economy, global carbon dioxide emissions have shown a clear growth trend in recent years. Vegetation is the most basic component of terrestrial ecosystems. Vegetation Carbon Sequestration (VCS) has attracted widespread attention as the most direct way to reduce carbon dioxide.

[0003] At present, the research on vegetation carbon sequestration is mostly focused on the correlation analysis of temporal and spatial variation trends and influencing factors, and few scholars have proposed quantitative evaluation methods and systems for the contribution rate of influencing factors to vegetation carbon sequestration. The use of traditional correlation analysis can only explore the relationship between a single influencing factor and vegetation carbon sequestration, and the failure to exclude the influence of other meteorological factors is not conducive to extracting the impact of human activities on vegetation carbon sequestration. Some studies have used the difference method between potential vegetation net primary productivity (Net Primary Productivity, NPP) and actual vegetation net primary productivity to distinguish the contribution of human activities and meteorological factors, but this method cannot distinguish the contribution of a single influencing factor to NPP, and it is also difficult to quantitatively estimate the contribution of influencing factors.

[0004] Therefore, it is necessary to provide a method for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration, which is used to evaluate the contribution rate of different driving factors to changes in vegetation carbon sequestration. Summary of the invention

[0005] The present invention provides a method for evaluating the contribution rate of driving factors to changes in vegetation carbon fixation, comprising: obtaining vegetation net primary productivity data of an area to be evaluated during a time period to be evaluated, and calculating vegetation carbon fixation data of the area to be evaluated during the time period to be evaluated based on the vegetation net primary productivity data of the area to be evaluated during the time period to be evaluated; obtaining driving factor data of the area to be evaluated during the time period to be evaluated, and calculating change trends of multiple driving factors based on the driving factor data of the area to be evaluated during the time period to be evaluated; calculating the correlation between multiple driving factors and vegetation carbon fixation based on the vegetation carbon fixation data and driving factor data of the area to be evaluated during the time period to be evaluated; and calculating the contribution rate of multiple driving factors to the changes in vegetation carbon fixation in the area to be evaluated during the time period to be evaluated based on the change trends of the multiple driving factors and the correlations between the multiple driving factors and vegetation carbon fixation.

[0006] Furthermore, the vegetation carbon sequestration data of the area to be assessed during the period to be assessed is calculated based on the following formula according to the vegetation net primary productivity data of the area to be assessed during the period to be assessed, including: ,in, is the vegetation in the area to be evaluated per unit time t and pixel The amount of CO2 fixed at is the vegetation in the area to be evaluated per unit time t and pixel of net primary productivity of vegetation.

[0007] Furthermore, the driving factor data of the area to be evaluated during the time period to be evaluated includes meteorological data of the area to be evaluated during the time period to be evaluated; based on the driving factor data of the area to be evaluated during the time period to be evaluated, the changing trends of multiple driving factors are calculated, including: based on the meteorological data of the area to be evaluated during the time period to be evaluated, the changing trends of multiple meteorological driving factors are calculated.

[0008] Furthermore, the change trends of multiple meteorological driving factors are calculated according to the meteorological data of the area to be evaluated during the time period to be evaluated by the following formula: ,in, is the changing trend of meteorological driving factors, To take the median value, Time index k Positional variables, Time index The variable at the position, n is the total number of time index positions.

[0009] Furthermore, based on the vegetation carbon fixation data and driving factor data of the area to be assessed during the time period to be assessed, the correlation between multiple driving factors and vegetation carbon fixation is calculated, including: for each of the meteorological driving factors, based on the vegetation carbon fixation data and meteorological data of the area to be assessed during the time period to be assessed, calculating the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data; for any two of the meteorological driving factors, based on the meteorological data of the area to be assessed during the time period to be assessed, calculating the correlation coefficient between the two meteorological driving factors; for each of the meteorological driving factors, based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficient of any two meteorological driving factors, calculating the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon fixation data, wherein n is the total number of driving factors.

[0010] Further, the multiple meteorological driving factors include temperature driving factors, precipitation driving factors and radiation driving factors; for each of the meteorological driving factors, the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon fixation data is calculated according to the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficients of any two meteorological driving factors, including: for the temperature driving factor, according to the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficients of any two meteorological driving factors, the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor are calculated; for the precipitation driving factor, according to the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficients of any two meteorological driving factors, the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the radiation driving factor are calculated; for the radiation driving factor, according to the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data The correlation coefficient of any two meteorological driving factors is used to calculate the first-order partial correlation coefficient of the radiation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, and the first-order partial correlation coefficient of the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor; according to the first-order partial correlation coefficient of the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, and the first-order partial correlation coefficient of the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor, and the first-order partial correlation coefficient of the precipitation driving factor and the vegetation carbon fixation data after excluding the radiation driving factor, The first-order partial correlation coefficients between the temperature driving factor and the vegetation carbon fixation data, as well as the first-order partial correlation coefficients between the radiation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, and the first-order partial correlation coefficients between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, were calculated; the second-order partial correlation coefficients between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the radiation driving factor, the second-order partial correlation coefficients between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the radiation driving factor, and the second-order partial correlation coefficients between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the temperature driving factor.

[0011] Furthermore, based on the changing trends of the multiple driving factors and the correlations between the multiple driving factors and vegetation carbon sequestration, the contribution rates of the multiple driving factors to the vegetation carbon sequestration changes in the to-be-assessed area during the to-be-assessed time period are calculated, including: calculating the contribution rates of the multiple driving factors to the vegetation carbon sequestration changes in the to-be-assessed area during the to-be-assessed time period based on the changing trends of the temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between the temperature driving factors and vegetation carbon sequestration data after excluding the precipitation driving factors and radiation driving factors, and the second-order partial correlation coefficients between the radiation driving factors and vegetation carbon sequestration data after excluding the precipitation driving factors and temperature driving factors, The contribution rate of human activities to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed is calculated according to the contribution rate of temperature driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed; the contribution rate of human activities to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed is calculated according to the contribution rate of temperature driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factor to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed.

[0012] Furthermore, based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between temperature driving factors and vegetation carbon fixation data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between precipitation driving factors and vegetation carbon fixation data after excluding temperature driving factors and radiation driving factors, and the second-order partial correlation coefficients between radiation driving factors and vegetation carbon fixation data after excluding precipitation driving factors and temperature driving factors, the contribution rate of temperature driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed are calculated, including: calculating the contribution rate of temperature driving factors, precipitation driving factors and radiation driving factors Contribution rate sign; based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient of temperature driving factors and vegetation carbon fixation data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient of precipitation driving factors and vegetation carbon fixation data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient of radiation driving factors and vegetation carbon fixation data after excluding precipitation driving factors and temperature driving factors, and the contribution rate signs of temperature driving factors, precipitation driving factors and radiation driving factors, calculate the contribution rate of temperature driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed.

[0013] Furthermore, the contribution rate sign is calculated according to the following formula: ,in, is the symbol of contribution rate, For VCS variables, It is the time index of the variable, indicating the position of the variable in the data set, positive contribution (+, indicating an increase in VCS), and negative contribution (-, indicating a decrease in VCS).

[0014] The present invention provides a system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration, which applies the above-mentioned method for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration, including: a data acquisition module, used to acquire vegetation net primary productivity data and driving factor data of an area to be evaluated during a time period to be evaluated; a data processing module, used to calculate the vegetation carbon sequestration data of the area to be evaluated during the time period to be evaluated based on the vegetation net primary productivity data of the area to be evaluated during the time period to be evaluated, and to calculate the change trends of multiple driving factors based on the driving factor data of the area to be evaluated during the time period to be evaluated, and also used to calculate the correlation between multiple driving factors and vegetation carbon sequestration based on the vegetation carbon sequestration data and driving factor data of the area to be evaluated during the time period to be evaluated; a contribution evaluation module, used to calculate the contribution rate of multiple driving factors to the change in vegetation carbon sequestration in the area to be evaluated during the time period to be evaluated based on the change trends of the multiple driving factors and the correlations between the multiple driving factors and vegetation carbon sequestration.

[0015] Compared with the prior art, the method and system for evaluating the contribution rate of driving factors to vegetation carbon fixation changes provided by the present invention have at least the following beneficial effects: Existing technologies only focus on the spatiotemporal variation trends of VCS and the correlation of influencing factors. When studying the relationship between multiple driving factors and VCS and the relationship between multiple driving factors, traditional correlation analysis can only explore the relationship between a single driving factor and VCS, and cannot rule out the influence of other meteorological driving factors. The use of partial derivatives can effectively rule out the influence of redundant variables, thereby more accurately reflecting the relationship between a single driving factor and VCS, and is also conducive to extracting the impact of human activities on VCS. The quantitative study of the contribution of multiple factors to VCS is achieved, and the quantitative resolution is high. Combined with the complex terrain and diverse vegetation of the area to be evaluated, a localized dynamic time-effectiveness modeling method is adopted, which has a high degree of refinement and strong scalability. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] This specification will be further described in the form of exemplary embodiments, which will be described in detail by the accompanying drawings. These embodiments are not restrictive, and in these embodiments, the same number represents the same structure, wherein: Figure 1 It is a flow chart of a method for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration according to some embodiments of this specification; Figure 2 It is a schematic diagram of modules of a system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration according to some embodiments of this specification; Figure 3 is a schematic diagram of test analysis results shown in some embodiments of this specification; Figure 4ais a schematic diagram of the contribution of temperature to VCS according to some embodiments of this specification; Figure 4b is a schematic diagram of the contribution of precipitation to VCS according to some embodiments of this specification; Figure 4c is a schematic diagram of the contribution of radiation to VCS according to some embodiments of this specification; Figure 4d is a schematic diagram of the contribution of human activities to VCS according to some embodiments of the present specification. DETAILED DESCRIPTION

[0017] In order to more clearly illustrate the technical solutions of the embodiments of this specification, the following is a brief introduction to the drawings required for the description of the embodiments. Obviously, the drawings described below are only some examples or embodiments of this specification. For ordinary technicians in this field, this specification can also be applied to other similar scenarios based on these drawings without creative work. Unless it is obvious from the language environment or otherwise explained, the same reference numerals in the figures represent the same structure or operation.

[0018] Figure 1 is a flow chart of a method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to some embodiments of this specification, such as Figure 1 As shown, the method for assessing the contribution of driving factors to changes in vegetation carbon sequestration may include the following steps.

[0019] Step 110, obtaining the vegetation net primary productivity data of the area to be assessed during the time period to be assessed, and calculating the vegetation carbon fixation data of the area to be assessed during the time period to be assessed based on the vegetation net primary productivity data of the area to be assessed during the time period to be assessed.

[0020] The amount of CO2 fixed by vegetation in a certain period of time can be calculated by calculating the organic matter photosynthesized by vegetation in the study area during the corresponding period combined with CO2. There are differences in vegetation types in different regions, and the feedback mechanism of carbon fixation on the ecological environment is different due to the influence of geographical location, climatic conditions, etc.

[0021] The vegetation carbon sequestration data of the area to be assessed during the period to be assessed can be calculated based on the following formula according to the vegetation net primary productivity data of the area to be assessed during the period to be assessed, including:

[0022] in, is the vegetation in the area to be evaluated per unit time t and pixel The amount of CO2 fixed at is the vegetation in the area to be evaluated per unit time t and pixel The net primary productivity of vegetation at the site.

[0023] The total amount of CO2 fixed in the area to be evaluated is the product of the amount of CO2 fixed by the vegetation in the area to be evaluated per unit time and pixel and the area.

[0024] Step 120 , obtaining driving factor data of the area to be evaluated in the time period to be evaluated, and calculating the change trends of multiple driving factors according to the driving factor data of the area to be evaluated in the time period to be evaluated.

[0025] Specifically, multiple driving factors are divided into meteorological driving factors (precipitation driving factors, temperature driving factors, radiation driving factors) and human activity driving factors as shown in Table 1.

[0026] Table 1 Classification unit precipitation mm temperature ℃ radiation MJ Human activities / In some embodiments, the driving factor data of the area to be evaluated during the time period to be evaluated includes meteorological data of the area to be evaluated during the time period to be evaluated.

[0027] According to the driving factor data of the area to be evaluated during the period to be evaluated, the changing trends of multiple driving factors are calculated, including: According to the meteorological data of the area to be evaluated during the period to be evaluated, the changing trends of multiple meteorological driving factors are calculated.

[0028] The Theil-Sen Median method was used to analyze the temporal and spatial variation trend of VCS in the evaluated area during the evaluated period, and the Mann-Kendall method was used for significance test. The Theil-Sen Median method estimates the long-term variation trend of variables by calculating the median in a long time series. This method can avoid discrete and erroneous data and reduce the impact of errors on the overall trend and judgment.

[0029] The following formula is used to calculate the changing trends of multiple meteorological driving factors based on the meteorological data of the area to be evaluated during the period to be evaluated:

[0030] in, is the changing trend of meteorological driving factors, To take the median value, and Represent the variables at different times, and n is the total number of time index positions. >0 indicates an upward trend, whereas a downward trend indicates a downward trend. The larger the value, the more obvious the increase or decrease trend of the meteorological driving factor in the evaluated area during the evaluated period. For example, taking the temperature driving factor as an example, the sample is the temperature data of a certain area from 2001 to 2020. n That is, take 20; perform a regression calculation on the temperature for 20 years, and use the data from the last year ( X k ) and the previous year's data ( X l ) According to the formula, we can get the changing trend of the temperature driving factor in the past 20 years.

[0031] The Mann-Kendall test is a long-term nonparametric statistical method that tests samples with non-normality or certain distribution characteristics. It is relatively less affected by outliers, can have missing values, and has strong quantitative analysis capabilities. The VCS of the area to be evaluated in the time period to be evaluated is analyzed based on the pixel scale through Theil-Sen Median trend analysis and Mann-Kendall test, and the Mann-Kendall test results are reclassified.

[0032] Based on the above VCS estimation and multi-year spatiotemporal trend analysis, the coefficient of variation method was further used to evaluate the stability and sustainability of VCS spatiotemporal changes in the study area. The value of the value is used to construct a spatial distribution characteristic model of the stability of VCS in the study area based on pixels during the evaluation period:

[0033] The coefficient of variation can reveal the fluctuation and stability of the variable during the study period, and also characterize the stability of the ecosystem in the corresponding time series. The larger the value, the greater the fluctuation of the variable and the lower the stability. On the contrary, the smaller the value, the more stable the variable. , indicating that the VCS increase / decrease trend is extremely stable; , indicating that the VCS increase / decrease change trend is stable; , indicating that the VCS increase / decrease change trend is unstable; , indicating that the VCS increase / decrease trend is extremely unstable.

[0034] Based on the stability assessment results, we can determine the VCS change trend in the study area during the study period and whether the VCS maintains a stable growth or steady decrease throughout the year. In the study period, VCS maintained a very stable increasing / decreasing trend, while The area where VCS changes during the study period is sometimes increasing and sometimes decreasing, with discontinuous mutations.

[0035] The coefficient of variation was used to analyze the stability of the VCS in the test area to be evaluated during the test period. The test analysis results are as follows: Figure 3 As shown, the interannual fluctuation of VCS in the experimental area for evaluation during the experimental period of evaluation is generally small, of which about 80.58% of the areas are stable (19.21% are extremely stable), and the extremely stable areas are mostly found in some areas in the northwest of the experimental area for evaluation; while 19.42% of the areas are unstable (4.74% are extremely unstable), among which the densely populated areas in the experimental area for evaluation have a greater impact on vegetation, and the urban centers in the experimental area for evaluation are mostly unstable. In addition, some high-altitude areas and rivers in the experimental area for evaluation also show obvious instability.

[0036] Step 130, calculating the correlation between multiple driving factors and vegetation carbon sequestration based on the vegetation carbon sequestration data and driving factor data of the area to be assessed during the time period to be assessed.

[0037] Specifically, correlation analysis can explore the relationship between VCS and a single driving factor, but the relationship cannot exclude the influence of other driving factors. When studying the relationship between multiple variables, the partial correlation coefficient can effectively exclude the influence of redundant variables, thereby more accurately reflecting the relationship between two variables. The partial correlation coefficient is also widely used to explore the relationship between vegetation and influencing factors. Based on the calculation of correlation and partial correlation, the second-order partial correlation coefficient is calculated to exclude the influence of other factors and explore the relationship between a single driving factor and VCS.

[0038] In some embodiments, step 130 specifically includes: For each meteorological driving factor, the correlation coefficient between each meteorological driving factor and vegetation carbon sequestration data is calculated based on the vegetation carbon sequestration data and meteorological data of the area to be assessed during the period to be assessed; For any two meteorological driving factors, the correlation coefficient of the two meteorological driving factors is calculated based on the meteorological data of the area to be evaluated during the period to be evaluated; For each meteorological driving factor, the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon sequestration data is calculated based on the correlation coefficient between each meteorological driving factor and the vegetation carbon sequestration data and the correlation coefficient between any two meteorological driving factors, where n is the total number of driving factors.

[0039] Multiple meteorological driving factors include temperature driving factors, precipitation driving factors and radiation driving factors.

[0040] In some embodiments, for each meteorological driving factor, the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon sequestration data is calculated based on the correlation coefficient between each meteorological driving factor and the vegetation carbon sequestration data and the correlation coefficients between any two meteorological driving factors, including: For the temperature driving factor, the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor were calculated based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficient between any two meteorological driving factors. For precipitation driving factors, based on the correlation coefficient between each meteorological driving factor and vegetation carbon fixation data and the correlation coefficient between any two meteorological driving factors, the first-order partial correlation coefficient between the precipitation driving factor and vegetation carbon fixation data after excluding the temperature driving factor and the first-order partial correlation coefficient between the precipitation driving factor and vegetation carbon fixation data after excluding the radiation driving factor were calculated; For the radiation driving factors, the first-order partial correlation coefficients between the radiation driving factors and the vegetation carbon fixation data after excluding the temperature driving factors and the first-order partial correlation coefficients between the radiation driving factors and the vegetation carbon fixation data after excluding the precipitation driving factors were calculated based on the correlation coefficients between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficients between any two meteorological driving factors. According to the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor, the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, and the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, the second-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the radiation driving factor, the second-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the radiation driving factor, and the second-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the temperature driving factor are calculated.

[0041] The first-order partial correlation coefficient can be calculated according to the following formula:

[0042] in, To exclude the influence of variable k, and The partial correlation of variables, for , The correlation between two variables.

[0043] The second-order partial correlation coefficient can be calculated according to the following formula:

[0044] in, To exclude , After the two variables and The partial correlation of variables, , , The partial correlation coefficients of the two driving factors and VCS were calculated respectively. The relationship between VCS and meteorological driving factors and human activity driving factors was explored through the second-order partial correlation coefficient, and the t-test method was used to test the significance of the results.

[0045] Step 140, based on the change trends of the multiple driving factors and the correlation between the multiple driving factors and vegetation carbon sequestration, calculate the contribution of the multiple driving factors to the change in vegetation carbon sequestration in the area to be assessed during the time period to be assessed.

[0046] The quantitative formula that characterizes the contribution of multiple driving factors to the change of vegetation carbon sequestration in the area to be assessed during the period to be assessed is:

[0047] In the above formula, Indicates the change amount of VCS. , , , , , and Respectively represent the contribution of radiation, precipitation, temperature, and human activities to the change in VCS. Each component in the right formula represents the product of the second-order partial correlation coefficient between the driving factor and VCS and the change trend of the driving factor. The contribution rate of each driving factor can be obtained by the ratio of the contribution of each driving factor to the change in VCS. For example, the contribution rate of the temperature driving factor to the change in vegetation carbon fixation can be obtained by and The ratio of is obtained.

[0048] In some embodiments, step 140 specifically includes: Based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between temperature driving factors and vegetation carbon fixation data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between precipitation driving factors and vegetation carbon fixation data after excluding temperature driving factors and radiation driving factors, and the second-order partial correlation coefficients between radiation driving factors and vegetation carbon fixation data after excluding precipitation driving factors and temperature driving factors, calculate the contribution rate of temperature driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed; Based on the contribution rate of temperature driving factors to the change of vegetation carbon sequestration in the assessment area during the assessment period, the contribution rate of precipitation driving factors to the change of vegetation carbon sequestration in the assessment area during the assessment period, and the contribution rate of radiation driving factors to the change of vegetation carbon sequestration in the assessment area during the assessment period, the contribution rate of human activities to the change of vegetation carbon sequestration in the assessment area during the assessment period is calculated.

[0049] The contribution sign is calculated according to the following formula:

[0050] in, is the symbol of contribution rate, For VCS variables, It is the time index of the variable, indicating the position of the variable in the data set, positive contribution (+, indicating an increase in VCS), and negative contribution (-, indicating a decrease in VCS).

[0051] The contribution rates of multiple driving factors to the experimental area to be evaluated during the experimental period to be evaluated are shown in Table 2.

[0052] Table 2 radiation precipitation temperature Human activities The overall area to be evaluated for the test 1.07（8.98%） 2.21（18.54%） 0.002（0.02%） 8.64（72.47%） Part of the area to be evaluated for testing A 1.07（5.20%） 4.39（21.34%） 0.002（0.01%） 15.11（73.45%） Part of the area to be evaluated for testing B 0.81（15.00%） 0.24（4.44%） 0.0004（0.01%） 4.35（80.55%） Part of the area to be evaluated for testing C 1.15（30.81%） 0.23（6.16%） 0.003（0.08%） 2.35（62.95%） Depend on Figures 4a-4dIt can be seen that the impact of temperature on the VCS of the experimental area to be evaluated during the experimental evaluation period is mainly positive, with the proportion of positive contribution areas accounting for about 60.08%, and the proportion of negative contribution areas accounting for 39.92%. It is relatively concentrated in the part of area A in the experimental area to be evaluated, showing that the northern part of the part of area A in the experimental area to be evaluated has a positive contribution, and the plain has a negative contribution. The spatial distribution characteristics of precipitation contribution are more obvious. It has a positive contribution to most areas of the experimental area to be evaluated, accounting for 70.38%, and the proportion of negative contribution areas is only 28.62%. The positive contribution to the basin area is more significant, among which the positive contribution of the basin center and its southern area to VCS exceeds 10 g / m2 / a. The positive contribution of radiation to the VCS of the experimental area to be evaluated accounts for 57.12% of the area, and the negative contribution accounts for about 42.88%. The contribution of human activities is mainly positive, with positive contribution areas accounting for 71.90% and negative contribution accounting for only 28.10%. In the basin area and its surrounding mountainous areas, the positive contribution of human activities exceeds 30 g / m2 / a, but in some area C of the experimental evaluation area, the negative contribution of human activities is more obvious. Secondly, around the two large urban centers, human activities also show obvious negative contribution.

[0053] Figure 2 is a schematic diagram of a module of a system for evaluating the contribution of driving factors to changes in vegetation carbon fixation according to some embodiments of this specification, such as Figure 2 As shown, the system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration may include a data acquisition module, a data processing module and a contribution evaluation module.

[0054] The data acquisition module can be used to obtain the vegetation net primary productivity data and driving factor data of the area to be evaluated during the time period to be evaluated; The data processing module can be used to calculate the vegetation carbon sequestration data of the area to be assessed during the time period to be assessed based on the vegetation net primary productivity data of the area to be assessed during the time period to be assessed, calculate the change trend of multiple driving factors based on the driving factor data of the area to be assessed during the time period to be assessed, and calculate the correlation between multiple driving factors and vegetation carbon sequestration based on the vegetation carbon sequestration data and driving factor data of the area to be assessed during the time period to be assessed; The contribution assessment module can be used to calculate the contribution rate of multiple driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed based on the change trends of multiple driving factors and the correlation between multiple driving factors and vegetation carbon sequestration.

[0055] The system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration can be used to execute the method for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration. For more descriptions of the system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration, please refer to the relevant descriptions of the system for evaluating the contribution rate of driving factors to changes in vegetation carbon sequestration, which will not be repeated here.

[0056] Finally, it should be understood that the embodiments described in this specification are only used to illustrate the principles of the embodiments of this specification. Other variations may also fall within the scope of this specification. Therefore, as an example and not a limitation, alternative configurations of the embodiments of this specification may be considered consistent with the teachings of this specification. Accordingly, the embodiments of this specification are not limited to the embodiments explicitly introduced and described in this specification.

Claims

1. A method for assessing the contribution of driving factors to changes in vegetation carbon sequestration, characterized in that: include: Obtaining vegetation net primary productivity data of the area to be assessed during the time period to be assessed, and calculating vegetation carbon sequestration data of the area to be assessed during the time period to be assessed based on the vegetation net primary productivity data of the area to be assessed during the time period to be assessed; Acquire driving factor data of the area to be evaluated during the time period to be evaluated, and calculate the change trends of multiple driving factors according to the driving factor data of the area to be evaluated during the time period to be evaluated; Calculating the correlation between multiple driving factors and vegetation carbon sequestration based on the vegetation carbon sequestration data and driving factor data of the area to be assessed during the time period to be assessed; Based on the changing trends of the multiple driving factors and the correlation between the multiple driving factors and vegetation carbon sequestration, the contribution rates of the multiple driving factors to the changes in vegetation carbon sequestration in the area to be assessed during the time period to be assessed are calculated.

2. The method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to claim 1, characterized in that: The vegetation carbon sequestration data of the area to be assessed during the period to be assessed is calculated based on the following formula according to the vegetation net primary productivity data of the area to be assessed during the period to be assessed, including: in, is the vegetation in the area to be evaluated per unit time t and pixel The amount of CO2 fixed at is the vegetation in the area to be evaluated per unit time t and pixel The net primary productivity of vegetation at the site.

3. The method for evaluating the contribution of driving factors to vegetation carbon sequestration changes according to claim 2, characterized in that: The driving factor data of the area to be evaluated in the time period to be evaluated includes the meteorological data of the area to be evaluated in the time period to be evaluated; According to the driving factor data of the area to be evaluated in the time period to be evaluated, the changing trends of multiple driving factors are calculated, including: According to the meteorological data of the area to be evaluated during the time period to be evaluated, the changing trends of multiple meteorological driving factors are calculated.

4. The method for evaluating the contribution of driving factors to vegetation carbon sequestration changes according to claim 3, characterized in that: The change trends of multiple meteorological driving factors are calculated according to the meteorological data of the area to be evaluated during the period to be evaluated by the following formula: in, is the changing trend of meteorological driving factors, To take the median value, Time index k Positional variables, Time index The variable at the position, n is the total number of time index positions.

5. The method for evaluating the contribution of driving factors to vegetation carbon sequestration changes according to claim 3 or 4, characterized in that: According to the vegetation carbon sequestration data and driving factor data of the area to be assessed during the period to be assessed, the correlation between multiple driving factors and vegetation carbon sequestration is calculated, including: For each of the meteorological driving factors, based on the vegetation carbon sequestration data and meteorological data of the area to be assessed during the time period to be assessed, calculate the correlation coefficient between each of the meteorological driving factors and the vegetation carbon sequestration data; For any two of the meteorological driving factors, based on the meteorological data of the area to be evaluated during the time period to be evaluated, calculate the correlation coefficient of the two meteorological driving factors; For each meteorological driving factor, the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon fixation data is calculated based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficient between any two meteorological driving factors, where n is the total number of driving factors.

6. The method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to claim 5, characterized in that: Multiple meteorological drivers include temperature drivers, precipitation drivers, and radiation drivers; For each of the meteorological driving factors, the n-1 order partial correlation coefficient between the meteorological driving factor and the vegetation carbon fixation data is calculated based on the correlation coefficient between each of the meteorological driving factors and the vegetation carbon fixation data and the correlation coefficients between any two of the meteorological driving factors, including: For the temperature driving factor, based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficient between any two meteorological driving factors, calculate the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor; For the precipitation driving factor, based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficients between any two meteorological driving factors, calculate the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the radiation driving factor; For the radiation driving factor, based on the correlation coefficient between each meteorological driving factor and the vegetation carbon fixation data and the correlation coefficient between any two meteorological driving factors, calculate the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor; According to the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, the first-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the radiation driving factor, the first-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor, and the first-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor, the second-order partial correlation coefficient between the temperature driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the radiation driving factor, the second-order partial correlation coefficient between the precipitation driving factor and the vegetation carbon fixation data after excluding the temperature driving factor and the radiation driving factor, and the second-order partial correlation coefficient between the radiation driving factor and the vegetation carbon fixation data after excluding the precipitation driving factor and the temperature driving factor are calculated.

7. The method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to claim 6, characterized in that: Based on the change trends of the multiple driving factors and the correlation between the multiple driving factors and vegetation carbon sequestration, the contribution rate of the multiple driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed is calculated, including: Based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between temperature driving factors and vegetation carbon sequestration data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficients between precipitation driving factors and vegetation carbon sequestration data after excluding temperature driving factors and radiation driving factors, and the second-order partial correlation coefficients between radiation driving factors and vegetation carbon sequestration data after excluding precipitation driving factors and temperature driving factors, calculate the contribution rate of temperature driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed; Based on the contribution rate of temperature driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of human activities to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed is calculated.

8. The method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to claim 7, characterized in that: Based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient between temperature driving factors and vegetation carbon fixation data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient between precipitation driving factors and vegetation carbon fixation data after excluding temperature driving factors and radiation driving factors, and the second-order partial correlation coefficient between radiation driving factors and vegetation carbon fixation data after excluding precipitation driving factors and temperature driving factors, the contribution rate of temperature driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed are calculated, including: Calculate the contribution rate signs of temperature driving factor, precipitation driving factor and radiation driving factor; Based on the changing trends of temperature driving factors, precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient between temperature driving factors and vegetation carbon sequestration data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient between precipitation driving factors and vegetation carbon sequestration data after excluding precipitation driving factors and radiation driving factors, the second-order partial correlation coefficient between radiation driving factors and vegetation carbon sequestration data after excluding precipitation driving factors and temperature driving factors, and the contribution rate signs of temperature driving factors, precipitation driving factors and radiation driving factors, the contribution rate of temperature driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, the contribution rate of precipitation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed, and the contribution rate of radiation driving factors to the change of vegetation carbon sequestration in the area to be assessed during the time period to be assessed are calculated.

9. The method for evaluating the contribution of driving factors to changes in vegetation carbon sequestration according to claim 8, characterized in that: The contribution sign is calculated according to the following formula: in, is the symbol of contribution rate, For VCS variables, The year index of the variable, indicating the position of the variable in the data set.

10. A system for evaluating the contribution of driving factors to changes in vegetation carbon sequestration, characterized in that: The method for evaluating the contribution rate of driving factors to vegetation carbon sequestration changes according to any one of claims 1 to 9 comprises: A data acquisition module is used to obtain the net primary productivity data and driving factor data of the vegetation in the area to be evaluated during the time period to be evaluated; A data processing module, used to calculate the vegetation carbon sequestration data of the area to be assessed during the time period to be assessed based on the vegetation net primary productivity data of the area to be assessed during the time period to be assessed, calculate the change trend of multiple driving factors based on the driving factor data of the area to be assessed during the time period to be assessed, and also used to calculate the correlation between multiple driving factors and vegetation carbon sequestration based on the vegetation carbon sequestration data and driving factor data of the area to be assessed during the time period to be assessed; The contribution assessment module is used to calculate the contribution rate of multiple driving factors to the change of vegetation carbon fixation in the area to be assessed during the time period to be assessed based on the change trends of the multiple driving factors and the correlation between the multiple driving factors and vegetation carbon fixation.

Images