Recovery priority zone identification method and device based on ecological vulnerability and recovery potential

Abstract

The invention provides a restoration priority zone identification method and device based on ecological vulnerability and restoration potential, and belongs to the technical field of ecological protection. The restoration priority area identification method based on ecological vulnerability and restoration potential comprises the steps that current ecological system service data, reference ecological system service data and ecological vulnerability data of a target area are acquired, and the reference ecological system service data are ecological system service data in a reference time period; according to the ecological vulnerability data, calculating an ecological vulnerability spatial distribution map of the target area; according to the current ecosystem service data and the reference ecosystem service data, obtaining an ecological restoration potential spatial distribution diagram of the target area; and according to the ecological restoration potential spatial distribution diagram and the ecological vulnerability spatial distribution diagram, identifying a restoration priority region in the target region.

Description

Technical Field

[0001] The present disclosure relates to the technical field of ecological protection, and more specifically, to a method and device for identifying restoration priority areas based on ecological fragility and restoration potential. Background Art

[0002] With the global population growth, economic and social development and climate change, the problems of biodiversity loss, habitat degradation and ecosystem function decline are becoming increasingly serious. In order to achieve sustainable development, it is necessary not only to strengthen ecological protection, but also to curb the process of ecological degradation through ecological restoration measures. Ecological restoration improves the structure of the ecosystem through human intervention and enhances the self-sustaining capacity of the system. It has become a key direction of future global development. In order to improve the efficiency and benefits of ecological restoration, it is crucial to identify priority areas for ecological restoration.

[0003] In the process of realizing the concept of the present disclosure, it was found that there are at least the following problems in the relevant technology: in the relevant technology, the method of identifying priority areas for ecological restoration often relies on successful cases in other regions or attempts to restore the ecosystem to a certain historical state, ignoring the uniqueness of the region and the dynamic changes of the environment, and cannot effectively reflect the comprehensive requirements of current ecological functions and socio-economic needs. Summary of the invention

[0004] In view of this, the present disclosure provides a method and apparatus for identifying restoration priority areas of ecological vulnerability and restoration potential.

[0005] One aspect of the present disclosure provides a method for identifying restoration priority areas based on ecological fragility and restoration potential, comprising: obtaining current ecosystem service data, reference ecosystem service data and ecological fragility data of a target area, wherein the reference ecosystem service data is ecosystem service data within a reference time period; calculating a spatial distribution map of ecological fragility of the target area based on the ecological fragility data; obtaining a spatial distribution map of ecological restoration potential of the target area based on the current ecosystem service data and the reference ecosystem service data; identifying restoration priority areas within the target area based on the spatial distribution map of ecological restoration potential and the spatial distribution map of ecological fragility.

[0006] According to an embodiment of the present disclosure, the above-mentioned obtaining the spatial distribution map of ecological restoration potential of the above-mentioned target area based on the above-mentioned current ecosystem service data and the above-mentioned reference ecosystem service data includes: obtaining the current ecosystem service spatial distribution map of the above-mentioned target area based on the above-mentioned current ecosystem service data; obtaining the optimal recovery state spatial distribution map of the above-mentioned target area based on the above-mentioned reference ecosystem service data; comparing the above-mentioned current ecosystem service spatial distribution map with the spatial distribution map of the above-mentioned optimal recovery state to obtain the ecological restoration potential spatial distribution map of the above-mentioned target area.

[0007] According to an embodiment of the present disclosure, the above-mentioned reference time period is divided into multiple sub-reference time periods; the above-mentioned reference ecosystem service data includes average climate condition data within the reference time period, and soil and terrain data within the above-mentioned multiple sub-reference time periods.

[0008] According to an embodiment of the present disclosure, the above-mentioned reference ecosystem service data is used to obtain the optimal recovery state spatial distribution map of the above-mentioned target area, including: calculating multiple groups of ecosystem service indicators based on the above-mentioned average climate condition data and the soil and terrain data in the above-mentioned multiple sub-reference time periods; selecting the maximum value of each indicator in the above-mentioned multiple groups of ecosystem service indicators to obtain the optimal recovery state spatial distribution map of the above-mentioned target area.

[0009] According to an embodiment of the present disclosure, the above-mentioned ecosystem service indicators include carbon sequestration, habitat quality, water conservation and soil conservation.

[0010] According to an embodiment of the present disclosure, the above-mentioned carbon sequestration includes above-ground biological carbon reserves, underground biological carbon reserves, soil carbon reserves and dead organic matter.

[0011] According to an embodiment of the present disclosure, the water conservation capacity is the difference between the water production and the net surface flow.

[0012] According to an embodiment of the present disclosure, the above-mentioned soil retention amount is the difference between the potential erosion amount and the actual erosion amount of the slope.

[0013] According to an embodiment of the present disclosure, the above habitat quality is calculated using formula (2):

[0014] (2);

[0015] in, is the habitat quality of grid x in land use type j, is the habitat suitability of land use type j, is the habitat degradation index of grid x, z is the normalization parameter, k is the half-saturation coefficient, and land use type j includes at least one of the following: cultivated land, forest, grassland, shrub, water area, construction land, unused land, snow / ice.

[0016] According to an embodiment of the present disclosure, the above-mentioned ecological vulnerability data include one or more of temperature, precipitation, terrain undulation, slope, habitat degradation factor, normalized difference vegetation index, landscape fragmentation, gross domestic product of the above-mentioned target area, and human footprint.

[0017] According to an embodiment of the present disclosure, the above-mentioned calculation of the ecological fragility spatial distribution map of the above-mentioned target area based on the above-mentioned ecological fragility data includes: using the hierarchical analysis method and the principal component analysis method to calculate the weights corresponding to each indicator in the above-mentioned ecological fragility data; generating the ecological fragility spatial distribution map of the above-mentioned target area according to each indicator in the above-mentioned ecological fragility data and the weights corresponding to each of the above indicators.

[0018] According to an embodiment of the present disclosure, the above-mentioned use of the hierarchical analysis method and the principal component analysis method to calculate the weights corresponding to each indicator in the above-mentioned ecological fragility data includes: using the above-mentioned hierarchical analysis method to obtain the first weight corresponding to each indicator in the above-mentioned ecological fragility data; using the above-mentioned principal component analysis method to calculate the second weight corresponding to each indicator in the above-mentioned ecological fragility data; according to the above-mentioned first weight and the above-mentioned second weight, obtaining the weight corresponding to each indicator in the above-mentioned ecological fragility data.

[0019] According to an embodiment of the present disclosure, the above-mentioned identifying the restoration priority area within the target area according to the above-mentioned ecological restoration potential spatial distribution map and the above-mentioned ecological fragility spatial distribution map includes: identifying the first ecological restoration potential area within the target area according to the above-mentioned ecological restoration potential spatial distribution map; identifying the first ecological fragile area within the target area according to the above-mentioned ecological fragility spatial distribution map; obtaining the above-mentioned restoration priority area according to the above-mentioned first ecological restoration potential area and the above-mentioned first ecological fragile area.

[0020] Another aspect of the present disclosure provides a device for identifying restoration priority areas based on ecological fragility and restoration potential, including: an acquisition module, used to acquire ecosystem service data and ecological fragility data of a target area; a calculation module, used to calculate a spatial distribution map of ecological fragility of the target area based on the above-mentioned ecological fragility data; an acquisition module, used to obtain a spatial distribution map of ecosystem services and a spatial distribution map of an optimal restoration state of the target area based on the above-mentioned ecosystem service data; a comparison module, used to compare the spatial distribution map of ecosystem services with the spatial distribution map of the above-mentioned optimal restoration state to obtain a spatial distribution map of ecological restoration potential of the target area; and an identification module, used to identify the restoration priority areas within the target area based on the spatial distribution map of ecological restoration potential and the spatial distribution map of ecological fragility.

[0021] Another aspect of the present disclosure provides an electronic device, comprising:

[0022] one or more processors;

[0023] a memory for storing one or more programs,

[0024] When the one or more programs are executed by the one or more processors, the one or more processors implement the method as described above.

[0025] Another aspect of the present disclosure provides a computer-readable storage medium storing computer-executable instructions, which are used to implement the above method when executed.

[0026] Another aspect of the present disclosure provides a computer program product, which includes computer executable instructions, and when the instructions are executed, are used to implement the method as described above.

[0027] According to the embodiments of the present disclosure, based on the ecological vulnerability data, the ecological vulnerability spatial distribution map of the target area is calculated, and the ecological vulnerability characterizes the background situation of the ecology; based on the current ecosystem service data and the reference ecosystem service data, the ecological restoration potential spatial distribution map is obtained; and then the ecological restoration potential spatial distribution map is combined with the ecological vulnerability spatial distribution map to identify the restoration priority area in the target area. The embodiments of the present disclosure combine the ecological restoration potential with the ecological vulnerability to construct a scientific and systematic restoration priority area identification method, which can reasonably determine the priority restoration area, thereby achieving efficient restoration of the ecosystem and optimal allocation of resources on a large scale. BRIEF DESCRIPTION OF THE DRAWINGS

[0028] The above and other objects, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which:

[0029] Figure 1 An exemplary system architecture of the restoration priority area identification method and apparatus based on ecological vulnerability and restoration potential disclosed in the present invention is schematically shown;

[0030] Figure 2 A flowchart of a method for identifying priority restoration areas based on ecological vulnerability and restoration potential according to an embodiment of the present disclosure is schematically shown;

[0031] Figure 3 A technical roadmap of a restoration priority area identification method based on ecological vulnerability and restoration potential according to an embodiment of the present disclosure is schematically shown;

[0032] Figure 4 A block diagram schematically shows a device for identifying restoration priority areas based on ecological fragility and restoration potential according to an embodiment of the present disclosure; and

[0033] Figure 5 A block diagram of an electronic device suitable for implementing a robot according to an embodiment of the present disclosure is schematically shown. DETAILED DESCRIPTION

[0034] Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the present disclosure. In the following detailed description, for ease of explanation, many specific details are set forth to provide a comprehensive understanding of the embodiments of the present disclosure. However, it is apparent that one or more embodiments may also be implemented without these specific details. In addition, in the following description, descriptions of known structures and technologies are omitted to avoid unnecessary confusion of the concepts of the present disclosure.

[0035] The terms used herein are only for describing specific embodiments and are not intended to limit the present disclosure. The terms "comprise", "include", etc. used herein indicate the existence of the features, steps, operations and / or components, but do not exclude the existence or addition of one or more other features, steps, operations or components.

[0036] All terms (including technical and scientific terms) used herein have the meanings commonly understood by those skilled in the art unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification and should not be interpreted in an idealized or overly rigid manner.

[0037] When using expressions such as "at least one of A, B, and C, etc.", they should generally be interpreted according to the meaning of the expression commonly understood by those skilled in the art (for example, "a system having at least one of A, B, and C" should include but is not limited to a system having A alone, B alone, C alone, A and B, A and C, B and C, and / or A, B, C, etc.).

[0038] In the embodiments of the present disclosure, the collection, updating, analysis, processing, use, transmission, provision, disclosure, storage, etc. of the data involved (for example, including but not limited to user personal information) are in compliance with the provisions of relevant laws and regulations, are used for legitimate purposes, and do not violate public order and good morals. In particular, necessary measures are taken for user personal information to prevent illegal access to user personal information data and maintain the security of user personal information and network security.

[0039] In order to improve the efficiency and benefits of ecological restoration, it is crucial to identify priority areas for ecological restoration. Ecosystem services are indicators of the level of ecological functions. One possible way to identify them is to use ecosystem services to evaluate the effects of ecological restoration and determine priority restoration areas. In related technologies, although the identification method of priority areas for ecological restoration can objectively describe the current ecological status, it has not yet clarified the goals of ecological restoration, and lacks clear standards to judge the success of restoration effects. Therefore, it is particularly important to introduce the optimal recovery state as an evaluation method. However, traditional methods often rely on successful cases in other regions or attempt to restore the ecosystem to a certain historical state, ignoring the uniqueness of the region and the dynamic changes in the environment, and cannot effectively reflect the comprehensive requirements of current ecological functions and socio-economic needs.

[0040] In addition, relevant technologies usually fail to fully consider the potential for ecosystem recovery, making it difficult to scientifically quantify the ability of ecosystems to recover after being disturbed. This shortcoming leads to a lack of targeted resource input, limiting the sustainability of ecological restoration effects. As far as ecological vulnerability assessment is concerned, relevant technologies often focus on the vulnerability of a single ecosystem or a specific area, failing to effectively combine it with actual restoration needs. Moreover, its division method usually lists areas with high ecosystem services and vulnerability as priority restoration targets, while areas with high ecological vulnerability but low service values ​​are ranked second. This approach can easily lead to inefficient resource utilization and cannot ensure the scientific division of priority areas for ecological restoration.

[0041] In the process of realizing the concept of the present disclosure, it is found that the optimal recovery state can be used as the goal and evaluation standard of ecological restoration, and ecological vulnerability characterizes the background situation of ecology. Combining the two can better divide the ecological restoration priority areas. The construction of the optimal recovery state should comprehensively consider all factors that affect the function of the ecosystem.

[0042] In view of this, an embodiment of the present disclosure provides a method for identifying restoration priority areas based on ecological fragility and restoration potential, including: obtaining current ecosystem service data, reference ecosystem service data, and ecological fragility data of a target area, wherein the reference ecosystem service data is the ecosystem service data within a reference time period; calculating a spatial distribution map of ecological fragility of the target area based on the ecological fragility data; obtaining a spatial distribution map of ecological restoration potential of the target area based on the current ecosystem service data and the reference ecosystem service data; and identifying restoration priority areas within the target area based on the spatial distribution map of ecological restoration potential and the spatial distribution map of ecological fragility.

[0043] Figure 1 The exemplary system architecture 100 to which the restoration priority area identification method and apparatus based on ecological vulnerability and restoration potential can be applied according to an embodiment of the present disclosure is schematically shown. It should be noted that: Figure 1What is shown is merely an example of a system architecture to which the embodiments of the present disclosure can be applied, in order to help those skilled in the art understand the technical content of the present disclosure, but it does not mean that the embodiments of the present disclosure cannot be used in other devices, systems, environments or scenarios.

[0044] like Figure 1 As shown, the system architecture 100 according to this embodiment may include a first terminal device 101, a second terminal device 102, a third terminal device 103, a network 104, and a server 105. The network 104 is used to provide a medium for a communication link between the first terminal device 101, the second terminal device 102, the third terminal device 103, and the server 105. The network 104 may include various connection types, such as wired and / or wireless communication links, etc.

[0045] The user may use the first terminal device 101, the second terminal device 102, and the third terminal device 103 to interact with the server 105 through the network 104 to receive or send messages, etc. Various communication client applications may be installed on the first terminal device 101, the second terminal device 102, and the third terminal device 103, such as shopping applications, web browser applications, search applications, instant messaging tools, email clients, and / or social platform software, etc. (only for example).

[0046] The first terminal device 101, the second terminal device 102, and the third terminal device 103 may be various electronic devices having display screens and supporting web browsing, including but not limited to smart phones, tablet computers, laptop computers, desktop computers, and the like.

[0047] The server 105 may be a server that provides various services, such as a background management server (only as an example) that provides support for websites browsed by users using the first terminal device 101, the second terminal device 102, and the third terminal device 103. The background management server may analyze and process the received data such as user requests, and feed back the processing results (such as web pages, information, or data obtained or generated according to user requests) to the terminal device.

[0048] It should be noted that the restoration priority zone identification method based on ecological fragility and restoration potential provided in the embodiment of the present disclosure can generally be executed by the server 105. Accordingly, the restoration priority zone identification device based on ecological fragility and restoration potential provided in the embodiment of the present disclosure can generally be set in the server 105. The restoration priority zone identification method based on ecological fragility and restoration potential provided in the embodiment of the present disclosure can also be executed by a server or server cluster that is different from the server 105 and can communicate with the first terminal device 101, the second terminal device 102, the third terminal device 103 and / or the server 105. Accordingly, the restoration priority zone identification device based on ecological fragility and restoration potential provided in the embodiment of the present disclosure can also be set in a server or server cluster that is different from the server 105 and can communicate with the first terminal device 101, the second terminal device 102, the third terminal device 103 and / or the server 105. Alternatively, the restoration priority zone identification method based on ecological fragility and restoration potential provided in the embodiment of the present disclosure may also be performed by the first terminal device 101, the second terminal device 102, or the third terminal device 103, or may also be performed by other terminal devices different from the first terminal device 101, the second terminal device 102, or the third terminal device 103. Accordingly, the restoration priority zone identification device based on ecological fragility and restoration potential provided in the embodiment of the present disclosure may also be set in the first terminal device 101, the second terminal device 102, or the third terminal device 103, or may be set in other terminal devices different from the first terminal device 101, the second terminal device 102, or the third terminal device 103.

[0049] For example, the ecosystem service data and ecological vulnerability data of the target area may be originally stored in any one of the first terminal device 101, the second terminal device 102 or the third terminal device 103 (for example, the first terminal device 101, but not limited thereto), or stored in an external storage device and may be imported into the first terminal device 101. Then, the first terminal device 101 may locally execute the restoration priority area identification method based on ecological vulnerability and restoration potential provided in the embodiment of the present disclosure, or send the ecosystem service data and ecological vulnerability data of the target area to other terminal devices, servers, or server clusters, and the other terminal devices, servers, or server clusters that receive the ecosystem service data and ecological vulnerability data of the target area may execute the restoration priority area identification method based on ecological vulnerability and restoration potential provided in the embodiment of the present disclosure.

[0050] It should be understood that Figure 1 The number of terminal devices, networks and servers in the embodiment is only for illustration. Any number of terminal devices, networks and servers may be provided according to implementation requirements.

[0051] Figure 2 The flowchart of the method for identifying priority restoration areas based on ecological fragility and restoration potential according to an embodiment of the present disclosure is schematically shown.

[0052] like Figure 2 As shown, the method includes operations S210 to S240.

[0053] In operation S210, current ecosystem service data, reference ecosystem service data, and ecological vulnerability data of a target area are acquired, wherein the reference ecosystem service data is ecosystem service data within a reference time period.

[0054] In operation S220, a spatial distribution map of ecological vulnerability of the target area is calculated based on the ecological vulnerability data.

[0055] In operation S230, a spatial distribution map of ecological restoration potential of the target area is obtained according to the current ecosystem service data and the reference ecosystem service data.

[0056] In operation S240, restoration priority areas within the target area are identified based on the ecological restoration potential spatial distribution map and the ecological vulnerability spatial distribution map.

[0057] According to the embodiments of the present disclosure, ecosystem services refer to all benefits that humans obtain from ecosystems, including provision services (such as providing food and water), regulation services (such as controlling floods and diseases), cultural services (such as spiritual, recreational and cultural benefits), and support services (such as maintaining the nutrient cycle of the earth's living environment). Ecosystem service data may include one or more of various types of data related to ecosystem services. For example, it may include one or more of provision service data, regulation service data, cultural service data, and support service data.

[0058] According to an embodiment of the present disclosure, the current ecosystem service data may be the ecosystem service data within the current time period. The reference ecosystem service data may be the ecosystem service data within the reference time period. Exemplarily, the current ecosystem service data may be the ecosystem service data of 2020, and the reference ecosystem service data may be the ecosystem service data of 1990 to 2020.

[0059] According to the embodiments of the present disclosure, the types of ecosystem service data can be flexible and diverse. The types of ecosystem service data can be adjusted according to the restoration demands of different regions. For example, the calculation of windbreak and sand fixation service data can be increased in the northwest region, while the calculation of grain yield service data can be emphasized in the northeast region.

[0060] According to the embodiments of the present disclosure, ecological vulnerability is the sensitive response and self-recovery ability of an ecosystem to external interference at a specific time and space scale, and is an inherent attribute of the ecosystem. Ecological vulnerability data may include one or more of various types of data related to ecological vulnerability. For example, it may include one or more of climate factor data, terrain factor data, environmental factor data, and human activity data.

[0061] According to an embodiment of the present disclosure, ecological restoration potential can be the ability of an ecosystem to return to a stable state after being disturbed. A more representative ecological restoration potential is constructed based on current ecosystem service data and reference ecosystem service data. By comparing the reference state with the current state, the ecological restoration potential can be quantified. The greater the external disturbance a region is subjected to, the greater the restoration effect that can be achieved. The spatial distribution map of ecological restoration potential can characterize the distribution of the ecological restoration potential of the target area in geographic space.

[0062] According to the embodiments of the present disclosure, ecological vulnerability, as a risk characterization of ecosystem disturbance, reflects the internal and external constraints of the ecosystem, complements the ecosystem services, and can more accurately delineate the priority areas for ecological restoration. Ecological vulnerability assessment has been widely used in the vulnerability analysis of species, habitats, and socio-economics. Incorporating it into the division of priority areas helps identify the main driving factors of ecological degradation, reflects the background conditions of the ecological environment, and provides a basis for taking differentiated management measures in different regions.

[0063] According to the embodiments of the present disclosure, based on the ecological vulnerability data, the ecological vulnerability spatial distribution map of the target area is calculated, and the ecological vulnerability characterizes the background situation of the ecology; based on the current ecosystem service data and the reference ecosystem service data, the ecological restoration potential spatial distribution map is obtained; and then the ecological restoration potential spatial distribution map is combined with the ecological vulnerability spatial distribution map to identify the restoration priority area in the target area. The embodiments of the present disclosure combine the ecological restoration potential with the ecological vulnerability to construct a scientific and systematic restoration priority area identification method, which can reasonably determine the priority restoration area, thereby achieving efficient restoration of the ecosystem and optimal allocation of resources on a large scale.

[0064] Reference below Figure 3 , combined with specific embodiments Figure 2 The method shown is further explained.

[0065] Figure 3 The technical roadmap of the method for identifying priority restoration areas based on ecological fragility and restoration potential according to an embodiment of the present disclosure is schematically shown.

[0066] According to the embodiments of the present disclosure, Figure 3As shown, according to the current ecosystem service data and the reference ecosystem service data, the spatial distribution map of the ecological restoration potential of the target area is obtained, including: according to the current ecosystem service data, obtaining the current ecosystem service spatial distribution map of the target area; according to the reference ecosystem service data, obtaining the optimal recovery state spatial distribution map of the target area; comparing the current ecosystem service spatial distribution map with the spatial distribution map of the optimal recovery state to obtain the ecological restoration potential spatial distribution map of the target area.

[0067] According to an embodiment of the present disclosure, the spatial distribution map of ecosystem services in the target area may characterize the distribution of ecosystem services in the target area in the geographic space.

[0068] According to an embodiment of the present disclosure, the optimal recovery state can be a relatively healthy, stable and sustainable state determined based on ecosystem service data. The spatial distribution map of the optimal recovery state can characterize the distribution of the optimal recovery state of the target area in geographic space. According to an embodiment of the present disclosure, a more representative optimal recovery state is constructed based on reference ecosystem service data, so as to more effectively reflect ecological functions and meet current socio-economic needs. The optimal recovery state can also be used as a goal and evaluation standard for ecological restoration. By comparing the optimal recovery state with the current status, the potential for ecological restoration can be quantified. Regions with great potential indicate that the greater the external disturbance they are subjected to, the greater the restoration effect that can be achieved.

[0069] According to an embodiment of the present disclosure, the reference time period may include multiple sub-reference time periods. The reference ecosystem service data includes average climate condition data within the reference time period, and soil and terrain data within multiple sub-reference time periods. Exemplarily, the reference time period may be 2016-2020, and the sub-reference time periods may be 2016, 2017, 2018, 2019, and 2020. The reference ecosystem service data includes average climate condition data from 2016 to 2020, and soil and terrain data from 2016, 2017, 2018, 2019, and 2020. The optimal recovery state spatial distribution map of the target area can be obtained based on the ecosystem service data from 2016 to 2020. The reference time period may also be the past 30 years, and the optimal recovery state spatial distribution map of the target area can be obtained based on the ecosystem service data from the past 30 years.

[0070] According to an embodiment of the present disclosure, obtaining a spatial distribution map of the optimal recovery state of a target area based on reference ecosystem service data includes: calculating multiple groups of ecosystem service indicators based on average climate condition data and soil and terrain data in multiple sub-reference time periods; selecting the maximum value of each indicator in the multiple groups of ecosystem service indicators to obtain a spatial distribution map of the optimal recovery state of the target area.

[0071] According to an embodiment of the present disclosure, multiple groups of ecosystem service indicators are calculated based on average climate condition data and soil and terrain data in multiple sub-reference time periods, taking into account historical climate change factors and temporal dynamics, so that the established optimal recovery state can reflect the impact of climate change on the ecological environment.

[0072] According to an embodiment of the present disclosure, the reference ecosystem service data may also include future predicted ecosystem service sub-data to construct a spatial distribution map of optimal restoration states at different time scales.

[0073] According to the embodiments of the present disclosure, Figure 3 As shown in Figure 2, ecosystem service indicators can include carbon sequestration (CS), habitat quality (HQ), water conservation (WC), and soil retention (SR).

[0074] According to an embodiment of the present disclosure, a comprehensive assessment model of ecosystem services and trade-offs, such as the InVEST 3.14.0 model, can be used to calculate the spatial distribution of ecosystem services in a target area and assign the same weights for overlay analysis.

[0075] According to an embodiment of the present disclosure, carbon sequestration includes above-ground biological carbon reserves, underground biological carbon reserves, soil carbon reserves and dead organic matter.

[0076] According to the embodiments of the present disclosure, in the calculation of carbon sequestration, aboveground organic matter, underground organic matter, soil carbon content and dead organic matter can be considered as the main sources of carbon, and unstable aboveground carbon and carbon stored in wood products can be ignored. Specifically, the calculation can be performed with reference to formula (1).

[0077] (1);

[0078] in, represents carbon sequestration, represents aboveground biomass carbon storage, represents the underground biomass carbon stock, represents soil carbon storage, Represents dead organic matter.

[0079] According to the embodiments of the present disclosure, the assessment of habitat quality can comprehensively consider the sensitivity of different habitats or vegetation types to three types of threat factors, namely, cultivated land, bare land, and built-up land, and their degree of degradation. By analyzing the impact of threat factors on habitats, the sensitivity of habitats, and the distance between the two, the habitat quality of the entire watershed can be calculated.

[0080] According to an embodiment of the present disclosure, the habitat quality can be calculated using formula (2):

[0081] (2);

[0082] in, is the habitat quality of grid x in land use type j, is the habitat suitability of land use type j, is the habitat degradation index of raster x, z is the normalization parameter, k is the half-saturation coefficient, and land use type j includes at least one of the following: cultivated land, forest, grassland, shrub, water area, construction land, unused land, snow / ice.

[0083] According to an embodiment of the present disclosure, in the above formula (2), the normalization parameter z may be 2.5, and the half-saturation constant k may be 0.5.

[0084] According to the embodiments of the present disclosure, by incorporating habitat degradation indicators into the ecosystem service habitat quality calculation module, the connection between ecosystem services and vulnerability is enhanced, making the vulnerability assessment more comprehensive.

[0085] According to an embodiment of the present disclosure, the water conservation capacity may be the difference between the water production and the net surface flow.

[0086] According to the embodiments of the present disclosure, the calculation of water conservation can be divided into two steps. First, the water yield of the basin is calculated in a comprehensive assessment model of ecosystem services and trade-offs, such as the InVEST model, using the water-heat coupling balance hypothesis. Second, referring to the method in the document "Delineation of Ecological Red Lines", the water yield is subtracted from the net surface flow to obtain the water conservation of the ecosystem. Specifically, the calculation can be performed with reference to formulas (3) to (5).

[0087] (3)

[0088] (4)

[0089] (5)

[0090] in, represents the annual actual evapotranspiration of grid cell x, Represents the annual precipitation in grid cell x (mm). is the surface runoff (mm), is the surface runoff coefficient, and Y(x) represents the water yield of grid cell x.

[0091] According to an embodiment of the present disclosure, the soil retention amount may be the difference between the potential erosion amount and the actual erosion amount of the slope.

[0092] According to the embodiments of the present disclosure, the actual and potential soil erosion of the slope can be obtained based on the soil and water loss equation by calculating the influencing factors such as soil, landform, climate and vegetation, and the soil conservation amount can be obtained by subtracting the actual erosion amount from the potential erosion amount of the slope. Specifically, the calculation can be performed with reference to formulas (6) to (8).

[0093] (6)

[0094] (7)

[0095] (8)

[0096] in, is the potential soil erosion, is the actual amount of soil erosion, is the soil conservation capacity (t km -2 ·a -1 ), R is rainfall erosion factor, K is soil erodibility factor, LS is slope factor, C is vegetation management factor, and P is soil and water conservation measures factor.

[0097] According to an embodiment of the present disclosure, obtaining a spatial distribution map of the optimal recovery state of a target area based on reference ecosystem service data may include: calculating multiple groups of ecosystem service indicators based on average climate condition data and soil and terrain data in multiple sub-reference time periods, wherein each group of ecosystem service indicator data includes carbon sequestration, habitat quality, water conservation, and soil conservation; obtaining a spatial distribution map of the optimal recovery state of the target area based on the maximum carbon sequestration, the maximum habitat quality, the maximum water conservation, and the maximum soil conservation in the multiple groups of ecosystem service indicators.

[0098] According to the embodiments of the present disclosure, the optimal recovery state based on ecosystem services can be constructed through long-term series of ecosystem service data, and the ecological restoration potential can be quantified.

[0099] For example, the optimal recovery state based on ecosystem services can be constructed by combining the average climate conditions and soil and topographic data of 30 years. Specifically, the average climate conditions from 1990 to 2020, as well as the soil and topographic data of 1990, 2000, 2010 and 2020, can be combined to calculate ecosystem services, and the optimal state of the four ecosystem services can be obtained by superimposing the maximum values. The data involved and their preprocessing methods can be referred to Table 1. The optimal state of the four ecosystem services is superimposed with the same weight to obtain the optimal recovery state. By comparing the distribution of ecosystem services in 2020 with the optimal recovery state, the ecological restoration potential can be obtained.

[0100] Table 1

[0101]

[0102] According to the embodiments of the present disclosure, in terms of ecological vulnerability assessment, in order to fully describe the ecological environment background, such as Figure 3 As shown in the figure, we can select nine indicators in four categories, namely, climate factors (temperature and precipitation), terrain factors (terrain relief and slope), environmental factors (habitat degradation factors, normalized difference vegetation index (NDVI) and landscape fragmentation), and human activities (gross domestic product (GDP), human footprint), and use the analytic hierarchy process (AHP) and principal component analysis (PCA) to calculate the weights of the indicators and generate a spatial distribution map of ecological vulnerability in the target area. Combine the ecological restoration potential with the ecological vulnerability level to identify the priority areas for ecological restoration.

[0103] According to an embodiment of the present disclosure, the ecological vulnerability data includes one or more of temperature, precipitation, terrain relief, slope, habitat degradation factor, normalized difference vegetation index, landscape fragmentation, gross domestic product of the target area, and human footprint.

[0104] According to the embodiments of the present disclosure, by selecting one or more of temperature, precipitation, terrain undulation, slope, habitat degradation factor, normalized difference vegetation index, landscape fragmentation, gross domestic product of the target area, and human footprint as ecological vulnerability data, the natural environment and socio-economic factors are integrated to fully describe the ecological background conditions.

[0105] According to the embodiments of the present disclosure, due to the dimensional differences between ecosystem service data and ecological vulnerability data, the range method can be used to first standardize them before subsequent calculation and analysis.

[0106] According to an embodiment of the present disclosure, ecosystem service data includes positive indicators and negative indicators. Positive indicators can be calculated with reference to formula (9), and negative indicators can be calculated with reference to formula (10). Each ecosystem service data can be regarded as a positive indicator.

[0107] (9)

[0108] (10)

[0109] in, is the standardized value of the jth data in the i-th row, is the original value of the data, is the minimum value of the jth data in the i-th row, The maximum value of the jth data in the i-th row.

[0110] The ecological vulnerability data and its properties provided by the embodiments of the present disclosure can be referred to Table 1.

[0111] Table 2

[0112]

[0113] According to an embodiment of the present disclosure, based on ecological fragility data, calculating the ecological fragility spatial distribution map of the target area includes: using the analytic hierarchy process (AHP) and the principal component analysis (PCA) to calculate the weight corresponding to each indicator in the ecological fragility data; generating the ecological fragility spatial distribution map of the target area according to each indicator in the ecological fragility data and the weight corresponding to each indicator.

[0114] According to an embodiment of the present disclosure, the weights corresponding to each indicator in the ecological fragility data are calculated using the hierarchical analysis method and the principal component analysis method, including: using the hierarchical analysis method to obtain the first weight corresponding to each indicator in the ecological fragility data; using the principal component analysis method to calculate the second weight corresponding to each indicator in the ecological fragility data; according to the first weight and the second weight, obtaining the weight corresponding to each indicator in the ecological fragility data.

[0115] According to the embodiments of the present disclosure, the analytic hierarchy process can be used to collect expert scores, construct a judgment matrix and calculate the consistency ratio, retain the scores that meet the consistency requirements, and obtain the weights of each indicator. Secondly, relying on a geographic information system platform, such as the ArcGIS 10.8 platform, first resample the data to the same resolution, and then use the spatial analysis module to perform principal component analysis, extract the main components that explain more than 80% of the variation rate, and further calculate the precise weights. Finally, the weights obtained by the analytic hierarchy process and the principal component analysis method are averaged to obtain the final weights, which are used to calculate ecological fragility and are divided into three categories: high, medium, and low by the natural breakpoint method. Specifically, the calculation formula refers to formula (11) to formula (17).

[0116] (11)

[0117] (12)

[0118] in, is the consistency index; is the maximum eigenvalue of the judgment matrix; n is the order of the judgment matrix and the number of judgment factors. is the consistency ratio, is the stochastic index.

[0119] (13)

[0120] (14)

[0121] (15)

[0122] in, and They represent the standardized principal component loading and eigenvector of each indicator j in the input layer i respectively. is the eigenvalue of each layer, is the proportion of the feature values ​​of each extraction layer, is the weight of index j calculated by principal component analysis.

[0123] (16)

[0124] (17)

[0125] in, represents the ecological vulnerability index; represents the standardized value of the jth indicator, Its weight value, is the weight of index j calculated by analytic hierarchy process

[0126] According to the embodiments of the present disclosure, the ecological fragility of the target area can be divided into three categories, high, medium and low, or two categories, high and low, by using the natural breakpoint method, the geometric spacing method, the standard deviation classification method, etc.; the ecological restoration potential of the target area can be divided into three categories, high, medium and low, or two categories, high and low.

[0127] According to the embodiments of the present disclosure, Figure 3 As shown, the restoration priority areas in the target area can be identified based on the spatial distribution map of ecological restoration potential and the spatial distribution map of ecological fragility, including: identifying the first ecological restoration potential area in the target area based on the spatial distribution map of ecological restoration potential; identifying the first ecological fragile area in the target area based on the spatial distribution map of ecological fragility; and obtaining the restoration priority area based on the first ecological restoration potential area and the first ecological fragile area.

[0128] According to an embodiment of the present disclosure, the first ecological restoration potential area may be an area with greater ecological restoration potential. The first ecologically fragile area may be an area with higher ecological fragility.

[0129] According to the embodiments of the present disclosure, areas with high ecological restoration potential and high vulnerability can be defined as restoration priority areas (Ecological restoration priority areas, ERPAs), and restoration projects need to be carried out first in these areas. Areas with high ecological restoration potential but low vulnerability are assisted restoration areas (Assisted restoration areas, ARAs), which are the second in restoration level and require human intervention to effectively restore ecosystem functions. For areas with low ecological restoration potential and high ecological vulnerability, more attention is paid to their protection rather than restoration, and they are named natural conservation areas (Natural conservation areas, NCAs). Finally, according to the identified areas and their unique ecological characteristics, regional mergers can be carried out to further divide the target area into large restoration zones. Unlike the zoning technology in the related art, the method provided in the embodiments of the present disclosure can give priority to carrying out restoration measures in areas with high ecological restoration potential and vulnerability, followed by areas with high restoration potential and low vulnerability; for areas with high vulnerability but limited restoration potential, more emphasis is placed on protection measures to maximize the effectiveness of ecological restoration and reduce unnecessary resource input.

[0130] The restoration priority area identification method based on ecological fragility and restoration potential provided in the embodiments of the present disclosure can be comprehensively used in large-scale ecological restoration projects to identify ecological restoration priority areas.

[0131] According to the embodiments of the present disclosure, the restoration priority areas identified by the entities can propose targeted ecological restoration measures to provide scientific guidance for the overall ecological restoration of the watershed.

[0132] When implementing an ecological restoration project, restoration work can be initiated preferentially in the restoration priority areas identified in the embodiments of the present disclosure to ensure the rapid recovery of ecosystem service capabilities; secondly, work can be carried out in auxiliary restoration areas to maximize the ecological environment service capabilities. In nature reserves, emphasis is placed on maintaining their protected status during the development process. At the same time, in large-scale ecological restoration projects, the ecological restoration areas demarcated often cover multiple regions. It is recommended that multiple departments collaborate to optimize resource allocation and jointly ensure the effective advancement of ecological restoration projects. For different restoration areas, targeted restoration suggestions can be made based on the changes in ecosystem services and the different weights of vulnerability indicators within them. For example, in areas with low vegetation coverage, measures such as planting trees, optimizing tree species configuration, and thinning and tending can be taken; in areas with large terrain undulations, improvements can be made through terrain transformation, riparian restoration, and ecological engineering. These measures will effectively promote ecological restoration and improve ecosystem service functions.

[0133] Figure 4The block diagram of a device for identifying priority restoration areas based on ecological fragility and restoration potential according to an embodiment of the present disclosure is schematically shown.

[0134] like Figure 4 As shown, the restoration priority area identification device 400 based on ecological fragility and restoration potential includes an acquisition module 410 , a calculation module 420 , an acquisition module 430 , and an identification module 440 .

[0135] The acquisition module 410 is used to acquire current ecosystem service data, reference ecosystem service data and ecological vulnerability data of the target area, wherein the reference ecosystem service data is the ecosystem service data within a reference time period.

[0136] The calculation module 420 is used to calculate the ecological vulnerability spatial distribution map of the target area according to the ecological vulnerability data.

[0137] The acquisition module 430 is used to obtain a spatial distribution map of ecological restoration potential of a target area based on current ecosystem service data and reference ecosystem service data.

[0138] The identification module 440 is used to identify the restoration priority areas in the target area according to the ecological restoration potential spatial distribution map and the ecological vulnerability spatial distribution map.

[0139] According to an embodiment of the present disclosure, the acquisition module 430 includes a first acquisition submodule, a second acquisition submodule and a comparison submodule. The first acquisition submodule is used to obtain the current ecosystem service spatial distribution map of the target area based on the current ecosystem service data; the second acquisition submodule is used to obtain the optimal recovery state spatial distribution map of the target area based on the reference ecosystem service data; the comparison submodule is used to compare the current ecosystem service spatial distribution map with the spatial distribution map of the optimal recovery state to obtain the ecological restoration potential spatial distribution map of the target area.

[0140] According to an embodiment of the present disclosure, the reference time period is divided into a plurality of sub-reference time periods; the reference ecosystem service data includes average climate condition data within the reference time period, and soil and terrain data within a plurality of sub-reference time periods.

[0141] According to an embodiment of the present disclosure, the second acquisition submodule includes a first calculation unit and a selection unit. The first calculation unit is used to calculate multiple groups of ecosystem service indicators based on average climate condition data and soil and terrain data in multiple sub-reference time periods; the selection unit is used to select the maximum value of each indicator in the multiple groups of ecosystem service indicators to obtain the optimal recovery state spatial distribution map of the target area.

[0142] According to an embodiment of the present disclosure, the calculation module 420 includes a calculation submodule and a generation submodule. The calculation submodule is used to calculate the weights corresponding to each indicator in the ecological vulnerability data using the hierarchical analysis method and the principal component analysis method; the generation submodule is used to generate the ecological vulnerability spatial distribution map of the target area according to each indicator in the ecological vulnerability data and the weight corresponding to each indicator.

[0143] According to the embodiment of the present disclosure, the calculation submodule includes a first obtaining unit, a second calculating unit and a second obtaining unit. The first obtaining unit is used to obtain the first weight corresponding to each indicator in the ecological vulnerability data by using the hierarchical analysis method; the second calculating unit is used to calculate the second weight corresponding to each indicator in the ecological vulnerability data by using the principal component analysis method; the second obtaining unit is used to obtain the weight corresponding to each indicator in the ecological vulnerability data according to the first weight and the second weight.

[0144] According to an embodiment of the present disclosure, the identification module 440 includes a first identification submodule, a second identification submodule and a third acquisition submodule. The first identification submodule is used to identify the first ecological restoration potential area in the target area according to the ecological restoration potential spatial distribution map; the second identification submodule is used to identify the first ecologically fragile area in the target area according to the ecologically fragile spatial distribution map; the third acquisition submodule is used to obtain the restoration priority area according to the first ecological restoration potential area and the first ecologically fragile area.

[0145] According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits, or at least part of the functions of any one of them can be implemented in one module. According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits can be split into multiple modules for implementation. According to the embodiments of the present invention, any one or more of the modules, submodules, units, and subunits can be at least partially implemented as hardware circuits, such as field programmable gate arrays (FPGAs), programmable logic arrays (PLAs), systems on chips, systems on substrates, systems on packages, application specific integrated circuits (ASICs), or can be implemented by hardware or firmware in any other reasonable way of integrating or packaging circuits, or can be implemented in any one of the three implementation methods of software, hardware, and firmware, or in any appropriate combination of any of them. Alternatively, according to the embodiments of the present invention, one or more of the modules, submodules, units, and subunits can be at least partially implemented as computer program modules, and when the computer program modules are run, the corresponding functions can be performed.

[0146] For example, any of the acquisition module 410, the calculation module 420, the acquisition module 430, and the identification module 440 can be combined in one module / unit / sub-unit for implementation, or any of the modules / units / sub-units can be split into multiple modules / units / sub-units. Alternatively, at least part of the functions of one or more of these modules / units / sub-units can be combined with at least part of the functions of other modules / units / sub-units and implemented in one module / unit / sub-unit. According to an embodiment of the present disclosure, at least one of the acquisition module 410, the calculation module 420, the acquisition module 430, and the identification module 440 can be at least partially implemented as a hardware circuit, such as a field programmable gate array (FPGA), a programmable logic array (PLA), a system on a chip, a system on a substrate, a system on a package, an application specific integrated circuit (ASIC), or can be implemented by hardware or firmware such as any other reasonable way of integrating or packaging the circuit, or by any one of the three implementation methods of software, hardware, and firmware, or by a suitable combination of any of them. Alternatively, at least one of the acquisition module 410 , the calculation module 420 , the obtaining module 430 , and the identification module 440 may be at least partially implemented as a computer program module, and when the computer program module is executed, a corresponding function may be performed.

[0147] It should be noted that the data processing system part in the embodiments of the present disclosure corresponds to the data processing method part in the embodiments of the present disclosure. The description of the data processing system part specifically refers to the data processing method part, which will not be repeated here.

[0148] Figure 5 A block diagram of an electronic device suitable for implementing the method described above according to an embodiment of the present disclosure is schematically shown. Figure 5 The electronic device shown is only an example and should not bring any limitation to the functions and scope of use of the embodiments of the present disclosure.

[0149] like Figure 5 As shown, the electronic device 500 according to an embodiment of the present disclosure includes a processor 501, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 502 or a program loaded from a storage part 508 to a random access memory (RAM) 503. The processor 501 may include, for example, a general-purpose microprocessor (such as a CPU), an instruction set processor and / or a related chipset and / or a special-purpose microprocessor (for example, an application-specific integrated circuit (ASIC)), etc. The processor 501 may also include an onboard memory for caching purposes. The processor 501 may include a single processing unit or multiple processing units for performing different actions of the method flow according to an embodiment of the present disclosure.

[0150] In RAM 503, various programs and data required for the operation of electronic device 500 are stored. Processor 501, ROM 502 and RAM 503 are connected to each other via bus 504. Processor 501 performs various operations of the method flow according to the embodiment of the present disclosure by executing the program in ROM 502 and / or RAM 503. It should be noted that the program can also be stored in one or more memories other than ROM 502 and RAM 503. Processor 501 can also perform various operations of the method flow according to the embodiment of the present disclosure by executing the program stored in the one or more memories.

[0151] According to an embodiment of the present disclosure, the electronic device 500 may further include an input / output (I / O) interface 505, which is also connected to the bus 504. The electronic device 500 may further include one or more of the following components connected to the input / output (I / O) interface 505: an input portion 506 including a keyboard, a mouse, etc.; an output portion 507 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage portion 508 including a hard disk, etc.; and a communication portion 509 including a network interface card such as a LAN card, a modem, etc. The communication portion 509 performs communication processing via a network such as the Internet. A drive 510 is also connected to the input / output (I / O) interface 505 as needed. A removable medium 511, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 510 as needed, so that a computer program read therefrom is installed into the storage portion 508 as needed.

[0152] According to an embodiment of the present disclosure, the method flow according to an embodiment of the present disclosure can be implemented as a computer software program. For example, an embodiment of the present disclosure includes a computer program product, which includes a computer program carried on a computer-readable storage medium, and the computer program contains a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from the network through the communication part 509, and / or installed from the removable medium 511. When the computer program is executed by the processor 501, the above-mentioned functions defined in the system of the embodiment of the present disclosure are executed. According to an embodiment of the present disclosure, the system, equipment, device, module, unit, etc. described above can be implemented by a computer program module.

[0153] The present disclosure also provides a computer-readable storage medium, which may be included in the device / apparatus / system described in the above embodiments; or may exist independently without being assembled into the device / apparatus / system. The above computer-readable storage medium carries one or more programs, and when the above one or more programs are executed, the method according to the embodiment of the present disclosure is implemented.

[0154] According to an embodiment of the present disclosure, the computer-readable storage medium may be a non-volatile computer-readable storage medium. For example, it may include, but is not limited to: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present disclosure, a computer-readable storage medium may be any tangible medium containing or storing a program that may be used by or in combination with an instruction execution system, apparatus, or device.

[0155] For example, according to an embodiment of the present disclosure, the computer-readable storage medium may include the ROM 502 and / or the RAM 503 described above and / or one or more memories other than the ROM 502 and the RAM 503 .

[0156] An embodiment of the present disclosure also includes a computer program product, which includes a computer program, which contains program code for executing the method provided by the embodiment of the present disclosure. When the computer program product is run on an electronic device, the program code is used to enable the electronic device to implement the method for identifying restoration priority areas based on ecological fragility and restoration potential provided by the embodiment of the present disclosure.

[0157] When the computer program is executed by the processor 501, the above functions defined in the system / device of the embodiment of the present disclosure are executed. According to the embodiment of the present disclosure, the system, device, module, unit, etc. described above can be implemented by a computer program module.

[0158] In one embodiment, the computer program may rely on tangible storage media such as optical storage devices, magnetic storage devices, etc. In another embodiment, the computer program may also be transmitted and distributed in the form of signals on a network medium, and downloaded and installed through the communication part 509, and / or installed from the removable medium 511. The program code contained in the computer program may be transmitted using any appropriate network medium, including but not limited to: wireless, wired, etc., or any suitable combination of the above.

[0159] According to an embodiment of the present disclosure, the program code for executing the computer program provided by the embodiment of the present disclosure can be written in any combination of one or more programming languages. Specifically, these computing programs can be implemented using high-level process and / or object-oriented programming languages, and / or assembly / machine languages. Programming languages ​​include, but are not limited to, Java, C++, python, "C" language or similar programming languages. The program code can be executed entirely on the user computing device, partially on the user device, partially on the remote computing device, or entirely on the remote computing device or server. In the case of a remote computing device, the remote computing device can be connected to the user computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or can be connected to an external computing device (for example, using an Internet service provider to connect through the Internet).

[0160] The flowcharts and block diagrams in the accompanying drawings illustrate the possible implementation architecture, functions and operations of the systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each box in the flowchart or block diagram may represent a module, a program segment, or a part of a code, and the above-mentioned module, program segment, or a part of the code contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box may also occur in an order different from that marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart, can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions. It can be understood by those skilled in the art that the features recorded in the various embodiments of the present disclosure can be combined and / or combined in a variety of ways, even if such a combination or combination is not explicitly recorded in the present disclosure. In particular, without departing from the spirit and teaching of the present disclosure, the features described in the various embodiments of the present disclosure may be combined and / or combined in a variety of ways. All of these combinations and / or combinations fall within the scope of the present disclosure.

[0161] The embodiments of the present disclosure are described above. However, these embodiments are only for illustrative purposes and are not intended to limit the scope of the present disclosure. Although the embodiments are described above, this does not mean that the measures in the various embodiments cannot be used in combination to advantage. Without departing from the scope of the present disclosure, those skilled in the art may make a variety of substitutions and modifications, which should all fall within the scope of the present disclosure.

Claims

1. A method for identifying restoration priority areas based on ecological vulnerability and restoration potential, including: Obtaining current ecosystem service data, reference ecosystem service data, and ecological vulnerability data of the target area, wherein the reference ecosystem service data is the ecosystem service data within a reference time period; Calculating and obtaining a spatial distribution map of ecological vulnerability of the target area based on the ecological vulnerability data; Obtaining a spatial distribution map of ecological restoration potential of the target area based on the current ecosystem service data and the reference ecosystem service data; Based on the spatial distribution map of ecological restoration potential and the spatial distribution map of ecological vulnerability, restoration priority areas within the target area are identified.

2. The method according to claim 1, wherein: The step of obtaining the spatial distribution map of ecological restoration potential of the target area according to the current ecosystem service data and the reference ecosystem service data includes: Obtaining a current ecosystem service spatial distribution map of the target area according to the current ecosystem service data; Obtaining a spatial distribution map of the optimal restoration state of the target area based on the reference ecosystem service data; The current ecosystem service spatial distribution map and the optimal recovery state spatial distribution map are compared to obtain the ecological restoration potential spatial distribution map of the target area.

3. The method according to claim 2, wherein: The reference time period is divided into a plurality of sub-reference time periods; The reference ecosystem service data include average climate condition data within a reference time period, and soil and topography data within the plurality of sub-reference time periods; The step of obtaining the optimal restoration state spatial distribution map of the target area according to the reference ecosystem service data comprises: Calculating multiple groups of ecosystem service indicators based on the average climate condition data and the soil and terrain data in the multiple sub-reference time periods; The maximum value of each indicator in the multiple groups of ecosystem service indicators is selected to obtain the optimal recovery state spatial distribution map of the target area.

4. The method according to claim 3, wherein: The ecosystem service indicators include carbon sequestration, habitat quality, water conservation and soil conservation; The carbon sequestration includes aboveground biomass carbon storage, underground biomass carbon storage, soil carbon storage and dead organic matter; The water conservation capacity is the difference between the water production and the net surface flow; The soil conservation amount is the difference between the potential erosion amount and the actual erosion amount of the slope surface; The habitat quality is calculated using formula (2): （2）； in, is the habitat quality of grid x in land use type j, is the habitat suitability of land use type j, is the habitat degradation index of grid x, z is the normalization parameter, k is the half-saturation coefficient, and land use type j includes at least one of the following: cultivated land, forest, grassland, shrub, water area, construction land, unused land, snow / ice.

5. The method according to claim 1, wherein: The ecological vulnerability data include one or more of temperature, precipitation, terrain relief, slope, habitat degradation factor, normalized difference vegetation index, landscape fragmentation, gross domestic product of the target area, and human footprint.

6. The method according to any one of claims 1 to 5, wherein The step of calculating the ecological vulnerability spatial distribution map of the target area according to the ecological vulnerability data includes: Using the analytic hierarchy process and principal component analysis method, the weights corresponding to each indicator in the ecological vulnerability data are calculated; A spatial distribution map of ecological fragility of the target area is generated according to each indicator in the ecological fragility data and the weight corresponding to each indicator.

7. The method according to claim 6, wherein: The calculation of the weights corresponding to the various indicators in the ecological vulnerability data by using the hierarchical analysis method and the principal component analysis method includes: Using the analytic hierarchy process, obtaining the first weight corresponding to each indicator in the ecological vulnerability data; Using the principal component analysis method, calculating the second weight corresponding to each indicator in the ecological vulnerability data; The weight corresponding to each indicator in the ecological vulnerability data is obtained according to the first weight and the second weight.

8. The method according to any one of claims 1 to 5, wherein The identifying the restoration priority areas within the target area according to the ecological restoration potential spatial distribution map and the ecological vulnerability spatial distribution map comprises: identifying a first ecological restoration potential area within the target area according to the ecological restoration potential spatial distribution map; According to the ecological vulnerability spatial distribution map, identifying a first ecologically vulnerable area in the target area; The restoration priority area is obtained based on the first ecological restoration potential area and the first ecological fragile area.

9. A restoration priority area identification device based on ecological vulnerability and restoration potential, comprising: The acquisition module is used to obtain ecosystem service data and ecological vulnerability data of the target area; A calculation module, used for calculating the ecological vulnerability spatial distribution map of the target area according to the ecological vulnerability data; An acquisition module, used to obtain an ecosystem service spatial distribution map and an optimal recovery state spatial distribution map of the target area according to the ecosystem service data; A comparison module, used for comparing the ecosystem service spatial distribution map with the spatial distribution map of the optimal recovery state to obtain a spatial distribution map of ecological restoration potential of the target area; An identification module is used to identify the restoration priority areas in the target area according to the ecological restoration potential spatial distribution map and the ecological vulnerability spatial distribution map.

10. An electronic device, comprising: one or more processors; a memory for storing one or more programs, When the one or more programs are executed by the one or more processors, the one or more processors implement the method according to any one of claims 1 to 8.

Images