Method and system for detecting self-adaptive recovery of ecosystem carbon sequestration capacity interference

Abstract

The application provides an ecosystem carbon fixation capacity interference recovery adaptive detection method and system, and relates to the field of ecological environment data processing.The method comprises the following steps: obtaining type-related information and carbon fixation capacity interference recovery detection data of a plurality of sample ecosystems based on a plurality of ecosystem type factors, interference identification factors, carbon fixation capacity indexes and interference recovery indexes, determining a plurality of ecosystem types and key interference identification factors, key carbon fixation capacity indexes and key interference recovery indexes of each ecosystem type; determining the ecosystem type to which the to-be-detected ecosystem belongs; obtaining interference identification data and carbon fixation capacity data based on the key interference identification factors and the key carbon fixation capacity indexes, determining an interference time period, obtaining interference recovery data in combination with the key interference recovery indexes, and performing carbon fixation capacity interference recovery progress evaluation, which has the advantages of improving the efficiency and accuracy of the carbon fixation capacity interference recovery detection of the ecosystem.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of ecological environment data processing, and in particular to an ecological system carbon fixation capacity interference recovery adaptive detection method and system. BACKGROUND

[0002] The ecological system itself is in dynamic balance and has certain elasticity and recovery force. Within the range of elasticity, interference fluctuation, i.e. moderate interference, can bring certain external stimulation and vitality to the ecological system, promoting species diversity and healthy development of the ecological system. However, when the interference becomes frequent or high-intensity, the recovery force of the ecological system may be challenged, leading to the ecological system failing to recover to the original balance state in time. Frequent or high-intensity ecological system interference may lead to the ecological system failing to recover its carbon fixation capacity in time, and long-term interference may make the ecological system change from a carbon sink to a carbon source. Ecological system interference has multi-scale characteristics in space and time, involving changes from micro to macro and from short-term to long-term. The multi-scale nature increases the difficulty of integrating different data and results.

[0003] The existing ecological system carbon fixation capacity interference recovery detection technology uses unified indicators and weights, which cannot accurately reflect the differences in carbon fixation mechanisms of different ecological systems. For example, forest ecosystems mainly rely on aboveground biomass carbon storage, while wetlands mainly rely on soil organic carbon accumulation. Unified indicators may mask the key carbon fixation processes of each system, leading to evaluation results deviating from reality. Secondly, data acquisition and integration is another challenge. The existing ecological system carbon fixation capacity interference recovery detection technology requires high-quality biological, climatic, hydrological and soil data, and the workload of invalid data is large, the cost is high and the detection efficiency is low.

[0004] Therefore, it is necessary to provide an ecological system carbon fixation capacity interference recovery adaptive detection method and system to improve the efficiency and accuracy of ecological system carbon fixation capacity interference recovery detection. SUMMARY

[0005] The present application provides an ecological system carbon sequestration capacity interference recovery adaptive detection method, comprising: determining a plurality of ecological system type factors, interference identification factors, carbon sequestration capacity indicators and interference recovery indicators; based on the plurality of ecological system type factors, interference identification factors, carbon sequestration capacity indicators and interference recovery indicators, obtaining type-related information and carbon sequestration capacity interference recovery detection data of a plurality of sample ecosystems; determining a plurality of ecological system types according to the type-related information of the plurality of sample ecosystems; determining key interference identification factors, key carbon sequestration capacity indicators and key interference recovery indicators of each ecological system type according to the carbon sequestration capacity interference recovery detection data of the plurality of sample ecosystems; determining the ecological system type to which the to-be-detected ecosystem belongs according to the type-related information of the to-be-detected ecosystem; obtaining interference identification data and carbon sequestration capacity data of the to-be-detected ecosystem based on the key interference identification factors and key carbon sequestration capacity indicators of the ecological system type to which the to-be-detected ecosystem belongs; determining the interference time period of the to-be-detected ecosystem according to the interference identification data and carbon sequestration capacity data of the to-be-detected ecosystem; obtaining interference recovery data of the to-be-detected ecosystem based on the interference time period of the to-be-detected ecosystem and the key interference recovery indicators of the ecological system type to which the to-be-detected ecosystem belongs, and performing carbon sequestration capacity interference recovery progress evaluation.

[0006] Further, the type-related information of the sample ecosystem at least includes the factor value of each ecological system type factor corresponding to the sample ecosystem; determining a plurality of ecological system types according to the type-related information of the plurality of sample ecosystems comprises: calculating the ecological system type similarity of any two sample ecosystems according to the factor value of each ecological system type factor corresponding to the two sample ecosystems; and dividing the plurality of sample ecosystems into a plurality of sample ecosystem groups according to the ecological system type similarity of any two sample ecosystems, wherein one sample ecosystem group corresponds to one ecological system type.

[0007] Further, the interference recovery detection data of the carbon fixation capacity of the sample ecosystem includes interference identification data, carbon fixation capacity data, and interference recovery data, wherein the interference identification data includes factor values of each interference identification factor at multiple sample time points, the carbon fixation capacity data includes index values of each carbon fixation capacity index at multiple sample time points, and the interference recovery data includes index values of each interference recovery index at multiple sample interference recovery time points; and the key interference identification factors of each ecosystem type are determined according to the interference recovery detection data of the carbon fixation capacity of the multiple sample ecosystems, including: for each ecosystem type, based on the index values of each carbon fixation capacity index corresponding to each sample ecosystem included in the ecosystem type at multiple sample time points, calculating a dynamic difference coefficient of each carbon fixation capacity index corresponding to the ecosystem type, and based on the dynamic difference coefficient of each carbon fixation capacity index corresponding to the ecosystem type, determining the key carbon fixation capacity index of the ecosystem type and the weight of each key carbon fixation capacity index.

[0008] Further, the key interference identification factors and the key interference recovery indexes of each ecosystem type are determined according to the interference recovery detection data of the carbon fixation capacity of the multiple sample ecosystems, including: for each ecosystem type, based on the key carbon fixation capacity index of the ecosystem type, determining candidate interference recovery indexes of the ecosystem type from the multiple interference recovery indexes, based on the index values of each candidate interference recovery index corresponding to each sample ecosystem included in the ecosystem type at multiple sample interference recovery time points, calculating a dynamic difference coefficient of each candidate interference recovery index corresponding to the ecosystem type, and according to the dynamic difference coefficient of each candidate interference recovery index corresponding to the ecosystem type, determining the key interference recovery index of the ecosystem type and the weight of each key interference recovery index; and for each ecosystem type, based on the index values of each key carbon fixation capacity index corresponding to each sample ecosystem included in the ecosystem type at multiple sample time points and the factor values of each interference identification factor at multiple sample time points, calculating an influence coefficient of each interference identification factor on each key carbon fixation capacity index, and based on the influence coefficient of each interference identification factor on each key carbon fixation capacity index, determining the key interference identification factor of the ecosystem type.

[0009] Further, based on the dynamic difference coefficient of each carbon fixation capacity index corresponding to the ecosystem type, the weight of each key carbon fixation capacity index of the ecosystem type is determined, including: calculating the first weight of each key carbon fixation capacity index of the ecosystem type by the analytic hierarchy process; determining the second weight of each key carbon fixation capacity index of the ecosystem type based on the dynamic difference coefficient of each carbon fixation capacity index corresponding to the ecosystem type; calculating the comprehensive weight of each key carbon fixation capacity index of the ecosystem type based on the first weight and the second weight of each key carbon fixation capacity index of the ecosystem type; determining the weight of each key carbon fixation capacity index of the ecosystem type based on the comprehensive weight of each key carbon fixation capacity index of the ecosystem type.

[0010] Further, based on the influence coefficient of each disturbance identification factor on each key carbon fixation capacity index, the key disturbance identification factor of the ecosystem type is determined, including: calculating the comprehensive key coefficient of the disturbance identification factor based on the influence coefficient of each disturbance identification factor on each key carbon fixation capacity index and the dynamic difference coefficient of each key carbon fixation capacity index corresponding to the ecosystem type; determining the key disturbance identification factor of the ecosystem type based on the comprehensive key coefficient of each disturbance identification factor.

[0011] Further, according to the type related information of the to-be-detected ecosystem, the ecosystem type to which the to-be-detected ecosystem belongs is determined, including: for each ecosystem type, determining the type related information corresponding to the ecosystem type based on the factor value of each sample ecosystem included in the ecosystem type corresponding to each ecosystem type factor; calculating the belonging coefficient of the to-be-detected ecosystem to each ecosystem type based on the type related information of the to-be-detected ecosystem and the type related information corresponding to each ecosystem type; determining the ecosystem type to which the to-be-detected ecosystem belongs based on the belonging coefficient of the to-be-detected ecosystem to each ecosystem type.

[0012] Further, the disturbance identification data of the to-be-detected ecosystem includes the index value of each key disturbance identification factor of the ecosystem type to which the to-be-detected ecosystem belongs at a plurality of historical time points, and the carbon fixation capacity data of the to-be-detected ecosystem includes the index value of each key disturbance identification factor of the ecosystem type to which the to-be-detected ecosystem belongs at a plurality of historical time points; according to the disturbance identification data and the carbon fixation capacity data of the to-be-detected ecosystem, the disturbance time period of the to-be-detected ecosystem is determined, including: determining the carbon fixation capacity baseline of each ecosystem type; determining the disturbance time period of the to-be-detected ecosystem based on the disturbance identification data, the carbon fixation capacity data of the to-be-detected ecosystem, and the carbon fixation capacity baseline of the ecosystem type to which the to-be-detected ecosystem belongs.

[0013] Further, the carbon fixation capacity index at least includes a water-vegetation coupling index, and the calculation formula is:

[0014]

[0015]

[0016]

[0017] wherein, is a water-vegetation coupling index, is an enhanced vegetation index, is a comprehensive reflectance ratio, is a near-infrared reflectance of the ecosystem, is a red light band reflectance of the ecosystem, is a blue light band reflectance of the ecosystem, is a short-wave infrared reflectance of the ecosystem, and is a coefficient, and is greater than 0, is a correction parameter.

[0018] The present application provides an ecosystem carbon sequestration capacity interference recovery adaptive detection system, which applies the above-mentioned ecosystem carbon sequestration capacity interference recovery adaptive detection method, comprising: an index determination module for determining a plurality of ecosystem type factors, interference identification factors, carbon sequestration capacity indexes and interference recovery indexes; a sample analysis module for obtaining type-related information and carbon sequestration capacity interference recovery detection data of a plurality of sample ecosystems based on the plurality of ecosystem type factors, interference identification factors, carbon sequestration capacity indexes and interference recovery indexes; determining a plurality of ecosystem types according to the type-related information of the plurality of sample ecosystems; determining key interference identification factors, key carbon sequestration capacity indexes and key interference recovery indexes of each ecosystem type according to the carbon sequestration capacity interference recovery detection data of the plurality of sample ecosystems; a recovery detection module for determining the ecosystem type to which the to-be-detected ecosystem belongs according to the type-related information of the to-be-detected ecosystem; obtaining interference identification data and carbon sequestration capacity data of the to-be-detected ecosystem based on the key interference identification factors and key carbon sequestration capacity indexes of the ecosystem type to which the to-be-detected ecosystem belongs; determining the interference time period of the to-be-detected ecosystem according to the interference identification data and carbon sequestration capacity data of the to-be-detected ecosystem; obtaining interference recovery data of the to-be-detected ecosystem based on the interference time period of the to-be-detected ecosystem and the key interference recovery indexes of the ecosystem type to which the to-be-detected ecosystem belongs, and performing carbon sequestration capacity interference recovery evaluation.

[0019] Compared with the prior art, the present application provides an ecosystem carbon sequestration capacity interference recovery adaptive detection method and system, which at least has the following beneficial effects:

[0020] 1. Existing technologies use a uniform indicator system, ignoring the differences in carbon sequestration mechanisms among different ecosystems (e.g., different carbon sequestration pathways in forests, wetlands, and farmland). This invention categorizes ecosystems by type factors and selects key disturbance identification factors and key carbon sequestration capacity indicators for each type, achieving precise detection with a "one type, one solution" approach. This reduces redundant data interference and improves the accuracy of disturbance identification and recovery assessment. For example, in arid forests, using "soil moisture" rather than "total precipitation" as the key factor is more consistent with the actual changes in carbon sequestration capacity.

[0021] 2. By calculating the dynamic difference coefficients of carbon sequestration capacity indicators and disturbance recovery indicators, key indicators sensitive to ecosystem changes are screened out, redundant data is eliminated, and the detection indicators are ensured to be strongly correlated with carbon sequestration capacity recovery. Key indicators are screened separately for different ecosystem types (e.g., forests, wetlands, farmland) to avoid a "one-size-fits-all" approach and improve the specificity of the detection results. Weights are assigned based on dynamic difference coefficients or influence coefficients to avoid subjective experience bias, allowing high-weight indicators (e.g., biomass, soil carbon density) to play a greater role in the comprehensive assessment and improve the accuracy of the assessment. By calculating the influence coefficients of disturbance identification factors on key carbon sequestration capacity indicators, the direct correlation between disturbance factors and changes in carbon sequestration capacity is clarified, avoiding misjudgments (e.g., misjudging natural fluctuations as disturbances). Key disturbance factors are determined based on the magnitude of the influence coefficients to improve the accuracy of disturbance identification.

[0022] 3. The carbon sequestration capacity baseline represents the carbon sequestration capacity level of an ecosystem under undisturbed or stable conditions. By comparing the historical carbon sequestration capacity data of the ecosystem under test with the baseline, the degree and duration of deviation in carbon sequestration capacity caused by disturbance can be quantified. Existing technologies may misjudge the disturbance period due to manually set thresholds that are too high or too low. In contrast, the determination of the disturbance period in this invention is based entirely on the comparison of historical data with the baseline, rather than human experience or threshold settings, thus avoiding human bias and improving the accuracy of assessing the recovery progress of carbon sequestration capacity disturbance. Attached Figure Description

[0023] This specification will be further described by way of exemplary embodiments, which will be described in detail with reference to the accompanying drawings. These embodiments are not limiting; in these embodiments, the same reference numerals denote the same structures, wherein:

[0024] Figure 1 This is a flowchart illustrating an adaptive detection method for the restoration of ecosystem carbon sequestration capacity under disturbance, as shown in some embodiments of this specification.

[0025] Figure 2 This is a schematic diagram of a module of an adaptive detection system for the disturbance recovery of ecosystem carbon sequestration capacity, as shown in some embodiments of this specification. Detailed Implementation

[0026] In order to more clearly illustrate the technical solutions of the embodiments of the present specification, the drawings needed to be used in the embodiment description will be briefly introduced as follows. Obviously, the drawings in the following description are only some examples or embodiments of the present specification, and for those skilled in the art, the present specification can also be applied to other similar scenarios without creative labor on the basis of these drawings. Unless it is obvious from the language environment or otherwise stated, the same reference numbers in the drawings represent the same structure or operation.

[0027] Figure 1 is a flowchart of an ecosystem carbon sequestration capacity interference recovery adaptive detection method according to some embodiments of the present specification, as shown in Figure 1 An ecosystem carbon sequestration capacity interference recovery adaptive detection method can include the following steps.

[0028] Step 110, determine a plurality of ecosystem type factors, interference identification factors, carbon sequestration capacity indicators and interference recovery indicators.

[0029] The ecosystem type factor reflects the characteristic parameters of the basic type, structure and function of the ecosystem, and is used to distinguish the core attributes of different ecosystems. For example, the plurality of ecosystem type factors can at least include vegetation type, biological community structure, climate zone, topographic feature, etc.

[0030] The factor value of the ecosystem in the plurality of ecosystem type factors can be determined in an encoded manner.

[0031] For example, the factor value of the ecosystem in the plurality of ecosystem type factors is determined in a three-level encoding manner of "main class-subclass-factor value".

[0032] For example, the main class of the vegetation type is valued at 2, and the factor value of the vegetation type in different ecosystems can be as shown in Table 1.

[0033] Table 1

[0034] Subclass code Main class Factor value 210 Single-layer structure 211 - Herb-dominated (coverage ≥ 80%); 212 - Shrub-dominated (coverage ≥ 60%) 220 Double-layer structure 221 - Shrub + Herb; 222 - Tree + Herb 230 Multi-layer structure 231 - Tree + Shrub + Herb (tropical rainforest); 232 - Tree + Herb (temperate forest)

[0035] The interference identification factor is used to identify and quantify the external pressure or internal change of the ecosystem. For example, the plurality of interference identification factors can at least include precipitation, temperature, etc.

[0036] The carbon sequestration capacity indicator is the core parameter for evaluating the carbon sink function of the ecosystem. For example, the plurality of carbon sequestration capacity indicators can be as shown in Table 2.

[0037] Table 2

[0038] Index category Specific index Vegetation carbon fixation Aboveground biomass carbon storage, net primary productivity Soil carbon fixation Soil organic carbon content, carbon turnover time Total carbon fixation Carbon storage per unit area (tons / hectare) Carbon fixation efficiency Carbon fixation per unit of photosynthetically active radiation

[0039] In some embodiments, the carbon sequestration capacity indicator comprises at least a water-vegetation coupling indicator, and the formula is:

[0040]

[0041]

[0042]

[0043] wherein, is the water-vegetation coupling indicator, is the enhanced vegetation index, is the comprehensive reflectance ratio, is the near-infrared reflectance of the ecosystem, is the red light band reflectance of the ecosystem, is the blue light band reflectance of the ecosystem, is the short-wave infrared reflectance of the ecosystem, and is a coefficient for correcting the interference of aerosol scattering existing in the atmosphere, and is greater than 0, for example, = 6 and = 7.5, is a correction parameter for correcting the interference of soil and tree crown background, for example, L is 1.

[0044] Specifically, the enhanced vegetation index (EVI) is a vegetation spectral index proposed on the basis of improving NDVI. EVI captures the photosynthetic activity of an area, i.e., its “greenness”, and can be used to estimate the ecosystem function of an area through its relationship with aboveground biomass production. In order to reduce the interference of water, clouds, heavy aerosols, etc., EVI uses the near-infrared and red light bands, and introduces the blue light band, and corrects the aerosol scattering existing in the atmosphere and the soil and tree crown background through the coefficient and the correction parameter.

[0045] Shortwave infrared / near infrared (SWIR / NIR) is sensitive to soil moisture, vegetation water status, and fire disturbance, while the near infrared band mainly reflects photosynthesis and chlorophyll content of vegetation. The ratio of SWIR to NIR can be introduced into the analysis to better understand the dynamic changes of the ecosystem under different wet and dry conditions. It can complement the shortcomings of EVI and more comprehensively capture changes in water, carbon cycling, and other factors in the ecosystem. Time-series remote sensing image data of Landsat series (Landsat 5, 7, 8), Sentinel-2, and MODIS can be obtained through Google Earth Engine (GEE), among which Landsat and Sentinel-2 provide high-resolution SWIR and NIR band data, and MODIS provides long-time series of large-scale EVI data. The cloud and shadow are removed using the cloud mask function (such as Fmask) built-in GEE, the image is standardized using the automatic atmospheric correction algorithm, and finally the Savitzky-Golay filter is used to smooth the time series to ensure data quality and reduce noise, to determine the near-infrared reflectance, red band reflectance, and blue band reflectance of the ecosystem.

[0046] CVI reflects the photosynthetic capacity (productivity) and water status (drought / wetness) of vegetation by integrating the information of vegetation greenness (EVI) and water-sensitive bands (SWIR / NIR). This coupling relationship is crucial for assessing the carbon sequestration capacity of ecosystems, as water conditions directly affect the photosynthetic efficiency of vegetation and the stability of soil carbon pools. CVI can be used to monitor the functional changes of ecosystems under different environmental conditions. For example, under drought stress, the EVI of vegetation may decrease due to water limitation, while the SWIR / NIR ratio may increase due to water deficit or soil drying, leading to a decrease in CVI value. This change can sensitively reflect the weakening of the carbon sequestration capacity of the ecosystem. CVI has the ability to respond to disturbances such as fire, drought, and land degradation. For example, after a fire, vegetation cover decreases, EVI decreases, and the SWIR / NIR ratio may significantly increase due to soil exposure or vegetation death, leading to a sharp decrease in CVI value. This change can be used to quickly identify disturbance events and their impact on carbon sequestration capacity.

[0047] The disturbance recovery index is used to assess the ability of the ecosystem to recover its carbon sequestration function after being disturbed.

[0048] Each carbon sequestration capacity index can correspond to at least one disturbance recovery index, and the disturbance recovery index corresponding to the carbon sequestration capacity index is used to reflect the changes of the carbon sequestration capacity index in the process of recovering the carbon sequestration function after disturbance.

[0049] For example, for the water-vegetation coupling index, the corresponding disturbance recovery index can include a water-vegetation coupling index growth rate, a water-vegetation coupling index recovery index, and the like.

[0050] For example, the water-vegetation coupling index growth rate and the water-vegetation coupling index recovery index can be calculated according to the following formula:

[0051]

[0052]

[0053] wherein, is the water-vegetation coupling index growth rate, is the index value of the current water-vegetation coupling index, is the index value of the initial water-vegetation coupling index after disturbance, is the recovery time, is the water-vegetation coupling index recovery index.

[0054] Step 120, based on the multiple ecosystem type factors, the disturbance identification factors, the carbon sequestration capacity index and the disturbance recovery index, obtaining the type related information of the multiple sample ecosystems and the carbon sequestration capacity disturbance recovery detection data.

[0055] Specifically, the type related information of the sample ecosystem at least includes the factor value of each ecosystem type factor corresponding to the sample ecosystem.

[0056] After determining the multiple ecosystem type factors, the disturbance identification factors, the carbon sequestration capacity index and the disturbance recovery index, different disturbances can be applied to the sample ecosystem, and the carbon sequestration capacity disturbance recovery detection data corresponding to the time period before and after the disturbance can be obtained.

[0057] The obtained type related information of the multiple sample ecosystems and the carbon sequestration capacity disturbance recovery detection data are subjected to standard processing, including data accuracy, format, metadata, projection, framing, integrity, spatial position processing and multi-scale conversion, to realize effective compilation of multi-source heterogeneous data, and the data is subjected to quality inspection to ensure the consistency and availability of the same data. The compiled data is subjected to database processing, and through the establishment of a database, users can conveniently query and search in multiple ways, and can conveniently browse the database data.

[0058] The carbon sequestration capacity disturbance recovery detection data of the sample ecosystem includes disturbance identification data, carbon sequestration capacity data and disturbance recovery data, wherein the disturbance identification data includes the factor value of each disturbance identification factor at multiple sample time points, the carbon sequestration capacity data includes the index value of each carbon sequestration capacity index at multiple sample time points, and the disturbance recovery data includes the index value of each disturbance recovery index at multiple sample disturbance recovery time points.

[0059] Step 130, according to the type of multiple sample ecosystem related information, determine the type of multiple ecosystem.

[0060] Specifically includes:

[0061] According to the factor value of each ecosystem type factor corresponding to any two sample ecosystems, the ecosystem type similarity of any two sample ecosystems is calculated, specifically, for each ecosystem type factor, if the factor values of the two sample ecosystems corresponding to the ecosystem type factor are consistent, the similarity of the two ecosystems corresponding to the ecosystem type factor is 1; if different, the similarity of the two ecosystems corresponding to the ecosystem type factor is 0, and the similarity of the two ecosystems corresponding to each ecosystem type factor is summed up to obtain the ecosystem type similarity of the two sample ecosystems;

[0062] According to the ecosystem type similarity of any two sample ecosystems, the multiple sample ecosystems are divided into multiple sample ecosystem groups, wherein a sample ecosystem group corresponds to an ecosystem type, that is, the sample ecosystems included in the sample ecosystem group belong to the same ecosystem type, for example, through clustering algorithm (such as K-Means clustering algorithm, hierarchical clustering algorithm, etc.), according to the ecosystem type similarity of any two sample ecosystems, the multiple sample ecosystems are divided into multiple sample ecosystem groups.

[0063] Step 140, according to the carbon sequestration capacity disturbance recovery detection data of multiple sample ecosystems, determine the key disturbance identification factor, key carbon sequestration capacity index and key disturbance recovery index of each ecosystem type.

[0064] Specifically includes:

[0065] For each ecosystem type, based on the index value of each carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type at multiple sample time points, the dynamic difference coefficient of each carbon sequestration capacity index corresponding to the ecosystem type is calculated, and based on the dynamic difference coefficient of each carbon sequestration capacity index corresponding to the ecosystem type, the key carbon sequestration capacity index of the ecosystem type and the weight of each key carbon sequestration capacity index are determined, specifically, for each sample ecosystem included in the ecosystem type, the standard deviation of the index value of a certain carbon sequestration capacity index corresponding to the sample ecosystem at multiple sample time points can be calculated, and the standard deviation of the index value of the carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type at multiple sample time points is averaged as the dynamic difference coefficient of the carbon sequestration capacity index corresponding to the ecosystem type, and the key carbon sequestration capacity index of the ecosystem type with the dynamic difference coefficient greater than the dynamic difference coefficient threshold;

[0066] For each ecosystem type, based on the key carbon sequestration capacity index of the ecosystem type, determine the candidate disturbance recovery index of the ecosystem type from the plurality of disturbance recovery indexes, calculate the dynamic difference coefficient of each candidate disturbance recovery index corresponding to the ecosystem type based on the index value of each candidate disturbance recovery index corresponding to each sample ecosystem included in the ecosystem type at the plurality of sample disturbance recovery time points, determine the key disturbance recovery index of the ecosystem type and the weight of each key disturbance recovery index according to the dynamic difference coefficient of each candidate disturbance recovery index corresponding to the ecosystem type. Specifically, the disturbance recovery index corresponding to the key carbon sequestration capacity index of the ecosystem type can be taken as the candidate disturbance recovery index, the dynamic difference coefficient of the candidate disturbance recovery index corresponding to the ecosystem type is calculated in a similar manner to the dynamic difference coefficient of the carbon sequestration capacity index corresponding to the ecosystem type, and the weight of the key disturbance recovery index is determined in a similar manner to the weight of the key carbon sequestration capacity index. Here, no further description is given.

[0067] For each ecosystem type, based on the index value of each key carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type at the plurality of sample time points and the factor value of each disturbance identification factor at the plurality of sample time points, calculate the influence coefficient of each disturbance identification factor on each key carbon sequestration capacity index, determine the key disturbance identification factor of the ecosystem type based on the influence coefficient of each disturbance identification factor on each key carbon sequestration capacity index. Specifically, for each sample ecosystem included in the ecosystem type, the correlation coefficient between the disturbance identification factor and the key carbon sequestration capacity index corresponding to the sample ecosystem can be calculated according to the correlation coefficient calculation algorithm (e.g., Pearson correlation coefficient, Spearman rank correlation coefficient, Kendall rank correlation coefficient, etc.) according to the index value of the key carbon sequestration capacity index corresponding to the sample ecosystem at the plurality of sample time points and the factor value of the disturbance identification factor at the plurality of sample time points. The correlation coefficient between the disturbance identification factor and the key carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type is averaged to obtain the influence coefficient of the disturbance identification factor on the key carbon sequestration capacity index.

[0068] In some embodiments, determining the weight of each key carbon sequestration capacity index of the ecosystem type based on the dynamic difference coefficient of the carbon sequestration capacity index corresponding to the ecosystem type comprises:

[0069] calculating the first weight of each key carbon sequestration capacity index of the ecosystem type by the analytic hierarchy process;

[0070] determining the second weight of each key carbon sequestration capacity index of the ecosystem type based on the dynamic difference coefficient of the carbon sequestration capacity index corresponding to the ecosystem type;

[0071] based on the first weight and the second weight of each key carbon fixation capacity indicator of the ecosystem type, calculating a comprehensive weight of each key disturbance identification factor of the ecosystem type, for example, performing a weighted sum of the first weight and the second weight of the key carbon fixation capacity indicator as the comprehensive weight of the key disturbance identification factor;

[0072] based on the comprehensive weight of each key carbon fixation capacity indicator of the ecosystem type, determining the weight of each key carbon fixation capacity indicator of the ecosystem type.

[0073] The weight of the key carbon fixation capacity indicator can be calculated according to the following formula:

[0074]

[0075] wherein, wi is the weight of the i th key carbon fixation capacity indicator of the ecosystem type, wi is the comprehensive weight of the i th key carbon fixation capacity indicator of the ecosystem type, wn is the comprehensive weight of the n th key carbon fixation capacity indicator of the ecosystem type, N is the total number of the key carbon fixation capacity indicators of the ecosystem type.

[0076] In some embodiments, based on the influence coefficient of each disturbance identification factor on each key carbon fixation capacity indicator, the key disturbance identification factor of the ecosystem type is determined, including:

[0077] based on the influence coefficient of each disturbance identification factor on each key carbon fixation capacity indicator and the dynamic difference coefficient of the corresponding ecosystem type of each key carbon fixation capacity indicator, calculating a comprehensive key coefficient of the disturbance identification factor, specifically, the dynamic difference coefficient of the corresponding ecosystem type of the key carbon fixation capacity indicator can be taken as a weight, and the influence coefficient of a certain disturbance identification factor on each key carbon fixation capacity indicator is weighted and summed as the comprehensive key coefficient of the disturbance identification factor;

[0078] based on the comprehensive key coefficient of each disturbance identification factor, the key disturbance identification factor of the ecosystem type is determined, for example, the disturbance identification factor with a comprehensive key coefficient greater than a comprehensive key coefficient threshold can be taken as the key disturbance identification factor of the ecosystem type.

[0079] Step 150, determining the ecosystem type to which the to-be-detected ecosystem belongs according to the type-related information of the to-be-detected ecosystem.

[0080] Specifically, it includes:

[0081] For each ecosystem type, based on the factor values ​​of each ecosystem type factor corresponding to each sample ecosystem included in the ecosystem type, the type-related information corresponding to the ecosystem type is determined. Specifically, the mean of the factor values ​​of each ecosystem type factor corresponding to each sample ecosystem included in the ecosystem type can be calculated to obtain the mean of the factor values ​​of each ecosystem type factor, which is used as the type-related information corresponding to the ecosystem type.

[0082] Based on the type-related information of the ecosystem to be detected and the type-related information corresponding to each ecosystem type, the attribution coefficient of the ecosystem to be detected to each ecosystem type is calculated.

[0083] Based on the association coefficient of the ecosystem to be tested with each ecosystem type, the ecosystem type to which the ecosystem to be tested belongs can be determined. For example, the ecosystem type with the largest association coefficient can be taken as the ecosystem type to which the ecosystem to be tested belongs.

[0084] For example, the classification coefficient of the ecosystem to be tested for each ecosystem type can be calculated using the following formula:

[0085]

[0086] in, Let be the coefficient of the ecosystem to be tested for the j-th ecosystem type. Let j be the factor value of the ecosystem type factor corresponding to the ecosystem to be detected. Let be the mean factor value of the factor corresponding to the j-th ecosystem type. The total number of ecosystem type factors. It is a positive integer.

[0087] Step 160: Based on the key disturbance identification factors and key carbon sequestration capacity indicators of the ecosystem type to which the ecosystem to be tested belongs, obtain the disturbance identification data and carbon sequestration capacity data of the ecosystem to be tested.

[0088] Specifically, the disturbance identification data of the ecosystem to be tested includes the index values ​​of each key disturbance identification factor of the ecosystem type to which the ecosystem to be tested belongs at multiple historical time points, and the carbon sequestration capacity data of the ecosystem to be tested includes the index values ​​of each key disturbance identification factor of the ecosystem type to which the ecosystem to be tested belongs at multiple historical time points.

[0089] Step 170: Determine the disturbance time period of the ecosystem to be tested based on the disturbance identification data and carbon sequestration capacity data of the ecosystem to be tested.

[0090] Specifically, it includes:

[0091] determining a baseline of carbon sequestration capacity for each type of ecosystem;

[0092] Based on the interference identification data of the to-be-detected ecosystem, the carbon sequestration capacity data, and the baseline of carbon sequestration capacity of the type of ecosystem to which the to-be-detected ecosystem belongs, the interference time period of the to-be-detected ecosystem is determined.

[0093] Specifically, for each type of ecosystem, using linear regression or moving average method, the trend of each key carbon sequestration capacity index over time is analyzed according to the index values of each key carbon sequestration capacity index of each sample ecosystem included in the type of ecosystem at the sample time points before the interference occurs, and the baseline of carbon sequestration capacity of the type of ecosystem is determined based on the trend, wherein the baseline of carbon sequestration capacity can include a baseline value corresponding to each key carbon sequestration capacity index.

[0094] The target time period can be determined based on the interference identification data of the to-be-detected ecosystem through an interference prediction model, wherein the interference prediction model can be a long short-term memory network model, and the target time period is a time period in which interference is likely to occur. The standard deviation of the index value of each key carbon sequestration capacity index is calculated according to the index values of each key carbon sequestration capacity index of the type of ecosystem to which the to-be-detected ecosystem belongs at multiple historical time points in the target time period, and the index difference between the mean value of the index value of each key carbon sequestration capacity index and the baseline value of the key carbon sequestration capacity index of the type of ecosystem to which the to-be-detected ecosystem belongs is calculated. Whether the target time period is an interference time period is determined according to the standard deviation of the index value of each key carbon sequestration capacity index and the index difference of each key carbon sequestration capacity index. For example, if the predicted value of a certain key carbon sequestration capacity index deviates from the historical mean value by more than 1.5 times the standard deviation, and / or if the difference between the predicted value of a certain key carbon sequestration capacity index and the baseline value exceeds 10%, it is determined that the target time period is an interference time period.

[0095] In step 180, interference recovery data of the to-be-detected ecosystem is obtained based on the interference time period of the to-be-detected ecosystem and the key interference recovery indicators of the type of ecosystem to which the to-be-detected ecosystem belongs, and carbon sequestration capacity interference recovery progress evaluation is performed.

[0096] Specifically, the interference recovery data of the to-be-detected ecosystem can include the index values of each key interference recovery indicator of the type of ecosystem to which the to-be-detected ecosystem belongs. The comprehensive evaluation model corresponding to the type of ecosystem to which the to-be-detected ecosystem belongs can be used to obtain the comprehensive recovery index of the to-be-detected ecosystem based on the comprehensive evaluation model, wherein the comprehensive evaluation model can be a multivariate linear regression model.

[0097] The evaluation result of the carbon sequestration capacity interference recovery progress evaluation can be:

[0098] Recovery success: the comprehensive recovery index reaches or exceeds the target value (such as 0.8), and the key interference recovery indicators are stable.

[0099] Recovery lag: the comprehensive recovery index is below the target value, and the key interference recovery indicators have not improved significantly.

[0100] Recovery failure: the comprehensive recovery index continues to decline or shows no significant improvement, and there may be secondary interference or ineffective management measures.

[0101] An ecosystem carbon sequestration capacity interference recovery adaptive detection method can be applied to forest ecosystem management, grassland ecosystem protection, wetland ecosystem restoration, and urban ecosystem optimization scenarios.

[0102] Taking forest ecosystem management as an example, in tropical rainforest or temperate forest areas, due to illegal logging, fires or pest damage, forest area reduction or vegetation degradation, real-time monitoring of carbon sequestration capacity changes and recovery progress evaluation are needed. Through the interference identification factor, the interference event (such as the time of logging) is quickly located. Based on the key carbon sequestration capacity indicators, the carbon sequestration capacity loss is quantified to guide the forest restoration priority. The interference recovery indicators are used to evaluate the effect of ecological restoration projects (such as afforestation).

[0103] Figure 2 According to some embodiments of the present specification, a module schematic diagram of an ecosystem carbon sequestration capacity interference recovery adaptive detection system is shown as Figure 2 The ecosystem carbon sequestration capacity interference recovery adaptive detection system can include an indicator determination module, a sample analysis module, and a recovery detection module.

[0104] The indicator determination module is used to determine multiple ecosystem type factors, interference identification factors, carbon sequestration capacity indicators, and interference recovery indicators.

[0105] The sample analysis module is used to obtain type-related information and carbon sequestration capacity interference recovery detection data of multiple sample ecosystems based on multiple ecosystem type factors, interference identification factors, carbon sequestration capacity indicators, and interference recovery indicators; determine multiple ecosystem types according to the type-related information of the multiple sample ecosystems; determine the key interference identification factor, the key carbon sequestration capacity indicator, and the key interference recovery indicator of each ecosystem type according to the carbon sequestration capacity interference recovery detection data of the multiple sample ecosystems.

[0106] The recovery detection module is configured to determine an ecosystem type to which the to-be-detected ecosystem belongs according to type-related information of the to-be-detected ecosystem; obtain interference identification data and carbon fixation capacity data of the to-be-detected ecosystem based on key interference identification factors and key carbon fixation capacity indexes of the ecosystem type to which the to-be-detected ecosystem belongs; determine an interference time period of the to-be-detected ecosystem according to the interference identification data and the carbon fixation capacity data of the to-be-detected ecosystem; and obtain interference recovery data of the to-be-detected ecosystem based on the interference time period of the to-be-detected ecosystem and key interference recovery indexes of the ecosystem type to which the to-be-detected ecosystem belongs, and perform carbon fixation capacity interference recovery evaluation.

[0107] An ecosystem carbon fixation capacity interference recovery adaptive detection system can be used to perform an ecosystem carbon fixation capacity interference recovery adaptive detection method, which will not be described herein again.

[0108] Finally, it should be understood that the embodiments described herein are merely for the purpose of illustrating the principles of the embodiments described herein. Other variations can also be within the scope of the embodiments described herein. Therefore, alternative configurations of the embodiments described herein can be considered as consistent with the teachings of the embodiments described herein, as examples rather than limitations. Accordingly, the embodiments described herein are not limited to the embodiments explicitly introduced and described herein.

Claims

1. An ecosystem carbon sequestration capacity interference recovery adaptive detection method, characterized in that, The method comprises the following steps: determining a plurality of ecosystem type factors, disturbance identification factors, carbon sequestration capacity indicators and disturbance recovery indicators, wherein the plurality of ecosystem type factors at least include vegetation type, biological community structure, climate zone and topographic feature, each carbon sequestration capacity indicator corresponds to at least one disturbance recovery indicator, and the disturbance recovery indicator corresponding to the carbon sequestration capacity indicator is used to reflect the change of the carbon sequestration capacity indicator in the process of recovering the carbon sequestration function after disturbance; based on the plurality of ecosystem type factors, disturbance identification factors, carbon sequestration capacity indicators and disturbance recovery indicators, obtaining type-related information and carbon sequestration capacity disturbance recovery detection data of a plurality of sample ecosystems, wherein the type-related information of the sample ecosystem at least includes the factor value of each ecosystem type factor corresponding to the sample ecosystem, different disturbances are applied to the sample ecosystem, and carbon sequestration capacity disturbance recovery detection data corresponding to the time period before and after the disturbance is obtained, the carbon sequestration capacity disturbance recovery detection data of the sample ecosystem includes disturbance identification data, carbon sequestration capacity data and disturbance recovery data, the disturbance identification data includes the factor value of each disturbance identification factor at a plurality of sample time points, the carbon sequestration capacity data includes the indicator value of each carbon sequestration capacity indicator at a plurality of sample time points, and the disturbance recovery data includes the indicator value of each disturbance recovery indicator at a plurality of sample disturbance recovery time points; determining a plurality of ecosystem types according to the type-related information of the plurality of sample ecosystems, specifically comprising: calculating the ecosystem type similarity of any two sample ecosystems according to the factor value of each ecosystem type factor corresponding to the two sample ecosystems, and dividing the plurality of sample ecosystems into a plurality of sample ecosystem groups according to the ecosystem type similarity of any two sample ecosystems, wherein one sample ecosystem group corresponds to one ecosystem type; determining the key disturbance identification factor, the key carbon sequestration capacity indicator and the key disturbance recovery indicator of each ecosystem type according to the carbon sequestration capacity disturbance recovery detection data of the plurality of sample ecosystems; determining the ecosystem type to which the to-be-detected ecosystem belongs according to the type-related information of the to-be-detected ecosystem; obtaining the disturbance identification data and the carbon sequestration capacity data of the to-be-detected ecosystem based on the key disturbance identification factor and the key carbon sequestration capacity indicator of the ecosystem type to which the to-be-detected ecosystem belongs; determining the disturbance time period of the to-be-detected ecosystem according to the disturbance identification data and the carbon sequestration capacity data of the to-be-detected ecosystem; obtaining the disturbance recovery data of the to-be-detected ecosystem based on the disturbance time period of the to-be-detected ecosystem and the key disturbance recovery indicator of the ecosystem type to which the to-be-detected ecosystem belongs, and evaluating the progress of carbon sequestration capacity disturbance recovery; wherein, determining the key disturbance identification factor, the key carbon sequestration capacity indicator and the key disturbance recovery indicator of each ecosystem type according to the carbon sequestration capacity disturbance recovery detection data of the plurality of sample ecosystems comprises: For each sample ecosystem included in the ecosystem type, a standard deviation of the index values of the carbon sequestration capacity index corresponding to the sample ecosystem at multiple sample time points is calculated, and a mean of the standard deviations of the index values of the carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type at multiple sample time points is calculated as a dynamic difference coefficient of the carbon sequestration capacity index corresponding to the ecosystem type, and a key carbon sequestration capacity index of the ecosystem type whose dynamic difference coefficient is greater than a dynamic difference coefficient threshold value is determined; The disturbance recovery index corresponding to the key carbon sequestration capacity index of the ecosystem type is taken as a candidate disturbance recovery index, and based on the index values of each candidate disturbance recovery index corresponding to each sample ecosystem included in the ecosystem type at multiple sample disturbance recovery time points, a dynamic difference coefficient of each candidate disturbance recovery index corresponding to the ecosystem type is calculated, and according to the dynamic difference coefficient of each candidate disturbance recovery index corresponding to the ecosystem type, a key disturbance recovery index of the ecosystem type is determined; For each sample ecosystem included in the ecosystem type, based on a correlation coefficient calculation algorithm, according to the index values of the key carbon sequestration capacity index corresponding to the sample ecosystem at multiple sample time points and the factor values of the disturbance identification factor at multiple sample time points, a correlation coefficient of the disturbance identification factor and the key carbon sequestration capacity index corresponding to the sample ecosystem is calculated, and a mean of the correlation coefficients of the disturbance identification factor and the key carbon sequestration capacity index corresponding to each sample ecosystem included in the ecosystem type is calculated as an influence coefficient of the disturbance identification factor on the key carbon sequestration capacity index; The dynamic difference coefficient of the ecosystem type corresponding to the key carbon sequestration capacity index is taken as a weight, and the influence coefficients of the disturbance identification factor on each key carbon sequestration capacity index are weighted and summed as a comprehensive key coefficient of the disturbance identification factor; The disturbance identification factor whose comprehensive key coefficient is greater than a comprehensive key coefficient threshold value is taken as a key disturbance identification factor of the ecosystem type; The carbon sequestration capacity index at least includes a water-vegetation coupling index, and the calculation formula is: wherein is a water-vegetation coupling index, is an enhanced vegetation index, is a composite reflectance ratio, is a near-infrared reflectance of the ecosystem, is a red light band reflectance of the ecosystem, is a blue light band reflectance of the ecosystem, is a shortwave infrared reflectance of the ecosystem, and is a coefficient, and is greater than 0, is a correction parameter.

2. The method according to claim 1, wherein, According to the type-related information of the to-be-detected ecosystem, the ecosystem type to which the to-be-detected ecosystem belongs is determined, including: For each ecosystem type, based on the factor values of each ecosystem type factor corresponding to each sample ecosystem included in the ecosystem type, type-related information corresponding to the ecosystem type is determined; Based on the type-related information of the to-be-detected ecosystem and the type-related information corresponding to each ecosystem type, an attribution coefficient of the to-be-detected ecosystem to each ecosystem type is calculated; Based on the attribution coefficient of the to-be-detected ecosystem to each ecosystem type, the ecosystem type to which the to-be-detected ecosystem belongs is determined.

3. The method of claim 1, wherein the method is characterized by, The disturbance identification data of the to-be-detected ecosystem includes the index values of each key disturbance identification factor of the ecosystem type to which the to-be-detected ecosystem belongs at multiple historical time points, and the carbon sequestration capacity data of the to-be-detected ecosystem includes the index values of each key disturbance identification factor of the ecosystem type to which the to-be-detected ecosystem belongs at multiple historical time points; According to the interference identification data and the carbon fixation capacity data of the to-be-detected ecosystem, a disturbance time period of the to-be-detected ecosystem is determined, including: determining a carbon fixation capacity baseline of each ecosystem type; based on the interference identification data, the carbon fixation capacity data of the to-be-detected ecosystem, and the carbon fixation capacity baseline of the ecosystem type to which the to-be-detected ecosystem belongs, determining the disturbance time period of the to-be-detected ecosystem.

4. An ecosystem carbon sequestration capacity interference recovery adaptive detection system, characterized in that, The application of the adaptive detection method for disturbance recovery of carbon fixation capacity of an ecosystem according to any one of claims 1-3, including: an index determination module for determining a plurality of ecosystem type factors, interference identification factors, carbon fixation capacity indexes, and disturbance recovery indexes; a sample analysis module for obtaining type-related information and carbon fixation capacity disturbance recovery detection data of a plurality of sample ecosystems based on the plurality of ecosystem type factors, interference identification factors, carbon fixation capacity indexes, and disturbance recovery indexes; determining a plurality of ecosystem types according to the type-related information of the plurality of sample ecosystems; determining key interference identification factors, key carbon fixation capacity indexes, and key disturbance recovery indexes of each ecosystem type according to the carbon fixation capacity disturbance recovery detection data of the plurality of sample ecosystems; a recovery detection module for determining the ecosystem type to which the to-be-detected ecosystem belongs according to the type-related information of the to-be-detected ecosystem; obtaining interference identification data and carbon fixation capacity data of the to-be-detected ecosystem based on the key interference identification factors and the key carbon fixation capacity indexes of the ecosystem type to which the to-be-detected ecosystem belongs; determining a disturbance time period of the to-be-detected ecosystem according to the interference identification data and the carbon fixation capacity data of the to-be-detected ecosystem; obtaining disturbance recovery data of the to-be-detected ecosystem based on the disturbance time period of the to-be-detected ecosystem and the key disturbance recovery indexes of the ecosystem type to which the to-be-detected ecosystem belongs, and performing carbon fixation capacity disturbance recovery evaluation.

Images