A method and system for quantitative simulation and evaluation of regional vegetation carbon storage

Abstract

The present invention discloses a method and system for quantitative simulation and assessment of regional vegetation carbon storage. First, the growth of plant potential biomass is calculated based on the leaf area index of the plant. Then, the actual biomass of the plant is calculated according to the growth of plant potential biomass and the growth stress conditions during the growth of the plant. The carbon density of the specified vegetation type is calculated based on the actual biomass of the specified plant type using the biomass conversion factor method. Finally, the area of ​​different vegetation types is obtained and combined with the carbon density of different vegetation types to obtain the regional vegetation carbon storage. The method and system for quantitative simulation and assessment of regional vegetation carbon storage provided by the embodiment of the present invention have a strong physical mechanism and can realize rapid simulation and assessment of surface vegetation carbon storage on a daily scale, in a large area, and for a long time.

Description

Technical Field

[0001] The present invention relates to the technical field of terrestrial carbon cycle, and in particular to a method and system for quantitative simulation and evaluation of regional vegetation carbon reserves. Background Art

[0002] The carbon storage of surface vegetation has a profound impact on the carbon cycle of terrestrial ecosystems, and thus plays an important role in mitigating global climate change. Traditional regional carbon storage assessment methods mainly include sample plot inventory method, micrometeorological method, biomass inventory method, stock method, etc. These methods often require a lot of manpower, material and financial resources, and it is difficult to achieve long-term, daily-scale simulation assessment of regional vegetation carbon storage; existing model assessment methods mainly include carbon balance models, vegetation-climate relationship models, etc. These methods generally lack clear physical mechanisms and cannot accurately and quickly assess the carbon storage of surface vegetation. Summary of the invention

[0003] Therefore, the embodiments of the present invention provide a method and system for quantitative simulation and evaluation of regional vegetation carbon storage, which can solve the problems existing in the prior art.

[0004] In order to achieve the above object, the present invention provides the following technical solutions:

[0005] In a first aspect, an embodiment of the present invention provides a method for quantitative simulation and assessment of regional vegetation carbon storage, comprising:

[0006] Calculate the growth of plant potential biomass based on the plant leaf area index;

[0007] Calculate the actual biomass of plants according to the growth of potential biomass of plants and the growth stress conditions during plant growth;

[0008] Based on the actual biomass of the specified plant type, the carbon density of the specified vegetation type is calculated using the biomass conversion factor method;

[0009] The area of ​​different vegetation types was obtained and combined with the carbon density of different vegetation types to obtain the regional vegetation carbon stock.

[0010] Further, the leaf area index of the plant is related to the ratio of the plant's accumulated heat units to the plant's potential heat units.

[0011] Furthermore, the leaf area index of the plant is calculated by the following formula:

[0012]

[0013] LAI n =f(pro n )

[0014] PHU represents the total cumulative heat unit required for plant maturity; n It indicates the ratio of plant potential heat units on the nth day, which refers to the ratio of the sum of daily cumulative heat units within n days of plant growth to the total cumulative heat units required for plant maturity; LAI n represents the leaf area index of the plant on day n, and f(·) represents the correlation function between the leaf area index of the plant and the ratio of the plant's potential heat unit.

[0015] Furthermore, the growth of the plant potential biomass is calculated by the following formula:

[0016]

[0017] Among them, P bio represents the total potential biomass growth of plants during the calculation period; N represents the total number of simulation days; △P bio，i represents the potential biomass growth of the plant on day i, calculated by the following formula:

[0018] ΔP bio =R·H e

[0019] Among them, △P bio It indicates the potential biomass growth of the plant on that day; R indicates the light energy utilization coefficient of the plant, which is related to the carbon dioxide concentration in the air and is calculated based on the carbon dioxide concentration in the air; H e It represents the daily effective photosynthesis solar radiation, which is calculated by the following formula:

[0020] H e =0.5·H d ·(1-exp·(-k l ·LAI+0.05))

[0021] Among them, H d Represents the total daily solar radiation; K l is the extinction coefficient, which indicates the ability of plants to absorb light; LAI indicates the leaf area index of plants on simulated days.

[0022] Furthermore, the actual growth of the plant is calculated by the following formula:

[0023]

[0024] Among them, P bio,act Represents the actual biomass of plants; P bio,ori represents the initial biomass of the plant; r cfrepresents the plant growth stress coefficient; wst represents the water stress factor, which is calculated by the actual plant evapotranspiration and the plant potential evapotranspiration; tst represents the temperature stress factor, which is a function of the daily average temperature and the optimal growth temperature of the plant.

[0025] Furthermore, the carbon density of vegetation types was calculated by the following formula:

[0026] C den,i =a·W i

[0027] Among them, C den,i represents the vegetation carbon density in the i-th land use type; W i represents the actual biomass of the i-th land use type; a represents the conversion factor.

[0028] Furthermore, the carbon storage of vegetation in the region is calculated by the following formula:

[0029]

[0030] Among them, C stor Represents the regional vegetation carbon storage; A i represents the area of ​​the i-th vegetation type; N represents the number of vegetation types.

[0031] In a second aspect, an embodiment of the present invention provides a regional vegetation carbon storage quantitative simulation and evaluation system, comprising:

[0032] A plant potential biomass growth calculation module is used to calculate the plant potential biomass growth based on the plant's leaf area index;

[0033] The plant actual biomass calculation module is used to calculate the plant actual biomass according to the growth of the plant potential biomass and the growth stress conditions during the plant growth process;

[0034] The carbon density calculation module of vegetation type is used to calculate the carbon density of the specified vegetation type based on the actual biomass of the specified plant type and using the biomass conversion factor method;

[0035] The regional vegetation carbon storage calculation module is used to obtain the area of ​​different vegetation types and combine the carbon density of different vegetation types to obtain the regional vegetation carbon storage.

[0036] In a third aspect, an embodiment of the present invention provides a computer device, comprising: at least one processor, and a memory communicatively connected to the at least one processor, wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor executes the method for quantitative simulation and assessment of regional vegetation carbon stocks as described in any embodiment of the first aspect of the present invention.

[0037] In a fourth aspect, an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method for quantitative simulation and assessment of regional vegetation carbon stocks as described in any embodiment of the first aspect of the embodiment of the present invention.

[0038] The technical solution of the present invention has the following advantages:

[0039] The present invention provides a method and system for quantitative simulation and assessment of regional vegetation carbon storage. First, the growth of plant potential biomass is calculated based on the leaf area index of the plant. Then, the actual biomass of the plant is calculated based on the growth of plant potential biomass and the growth stress conditions during plant growth. The carbon density of the specified vegetation type is calculated using the biomass conversion factor method based on the actual biomass of the specified plant type. Finally, the area of ​​different vegetation types is obtained and combined with the carbon density of different vegetation types to obtain the regional vegetation carbon storage. The method and system for quantitative simulation and assessment of regional vegetation carbon storage provided by the embodiment of the present invention have a strong physical mechanism and can realize rapid simulation and assessment of surface vegetation carbon storage on a daily scale, in a large area, and for a long time. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] In order to more clearly illustrate the specific implementation methods of the present invention or the technical solutions in the prior art, the drawings required for use in the specific implementation methods or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some implementation methods of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0041] Figure 1 A flowchart of a specific example of a method for quantitative simulation and assessment of regional vegetation carbon storage provided in an embodiment of the present invention;

[0042] Figure 2 FIG. 1 is a module composition diagram of an example of a regional vegetation carbon storage quantitative simulation and evaluation system provided in an embodiment of the present invention;

[0043] Figure 3 A composition diagram of a specific example of a computer device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0044] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0045] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0046] Example

[0047] The embodiment of the present invention provides a method for quantitative simulation and assessment of regional vegetation carbon storage, such as Figure 1 As shown, the method comprises the following steps:

[0048] Step S1: Calculate the growth of plant potential biomass based on the leaf area index of the plant.

[0049] Specifically, leaf area index refers to the ratio of the total area of ​​plant leaves per unit land area to the land area. It is an important indicator reflecting the growth status of plant populations and an important parameter for calculating plant biomass growth.

[0050] Temperature is one of the important factors controlling plant growth. Each plant has its minimum, maximum and optimal growth temperature. When the average temperature is lower than the base temperature for crop growth, the crop will not grow. Only the part of the daily average temperature that exceeds the minimum temperature (base temperature) for plant growth contributes to plant growth. If the temperature exceeds the maximum, the plant will stop growing. The leaf area index of a plant is closely related to the heat units accumulated by the plant.

[0051] When the daily average temperature is higher than the minimum temperature for plant growth, the heat units accumulated by the plant per day are calculated as follows:

[0052] HU i =T av,i -T base

[0053] Among them, HU i represents the cumulative heat unit of the i-th day; T av,i represents the daily average temperature of the i-th day; T base Indicates the minimum temperature for plant growth (basal temperature).

[0054] The leaf area index of a plant is related to the plant's potential heat unit ratio:

[0055]

[0056] LAI n =f(pro n )

[0057] PHU represents the total cumulative heat unit required for plant maturity; n It indicates the ratio of plant potential heat units on the nth day, which refers to the ratio of the sum of daily cumulative heat units within n days of plant growth to the total cumulative heat units required for plant maturity; LAI nrepresents the leaf area index of the plant on day n, and f(·) represents the correlation function between the leaf area index of the plant and the ratio of the plant's potential heat unit.

[0058] It should be noted that the leaf area index of a plant can be calculated using other embodiments, for example, by measuring radiation transmittance or by remote sensing estimation. The embodiment of the present invention uses the correlation between the leaf area index of a plant and the ratio of the accumulated heat units of the plant and the potential heat units of the plant as an example, but is not limited to this.

[0059] Step S2: Calculate the actual biomass of the plant according to the growth of the plant's potential biomass and the growth stress conditions during the plant's growth process.

[0060] In the embodiment of the present invention, the growth amount of plant potential biomass is calculated by the following formula:

[0061]

[0062] Among them, P bio represents the total potential biomass growth of plants during the calculation period; N represents the total number of simulation days; △P bio，i represents the potential biomass growth of the plant on day i, calculated by the following formula:

[0063] ΔP bio =R·H e

[0064] Among them, △P bio It indicates the potential biomass growth of the plant on that day; R indicates the light energy utilization coefficient of the plant, which is related to the carbon dioxide concentration in the air and is calculated based on the carbon dioxide concentration in the air; H e It represents the daily effective photosynthesis solar radiation, which is calculated by the following formula:

[0065] H e =0.5·H d ·(1-exp·(-k l ·LAI+0.05))

[0066] Among them, H d Represents the total daily solar radiation; K l It is the extinction coefficient, which indicates the light absorption capacity of plants, and its value is generally around 0.65; LAI indicates the leaf area index of plants on simulated days.

[0067] In the actual production process, plants are affected by insufficient water and extreme temperatures, which ultimately leads to the difference between actual growth and potential biomass growth. The actual growth of plants in the embodiment of the present invention is calculated using the following formula:

[0068]

[0069] Among them, P bio,act Represents the actual biomass of plants; P bio,ori represents the initial biomass of the plant; r cf represents the plant growth stress coefficient; wst represents the water stress factor, which is calculated by the actual plant evapotranspiration and the plant potential evapotranspiration; tst represents the temperature stress factor, which is a function of the daily average temperature and the optimal growth temperature of the plant.

[0070] It should be noted that the plant growth stress coefficient in the embodiment of the present invention is calculated by the water stress factor and temperature stress factor that have a greater impact on plant growth. However, in practice, plants are also affected by growth conditions such as light and oxygen content. Therefore, according to the actual plant type, the plant growth stress coefficient can also be calculated based on more growth condition stress factors, and no specific restrictions are made here.

[0071] Step S3: Based on the actual biomass of the specified plant type, the carbon density of the specified vegetation type is calculated using the biomass conversion factor method.

[0072] The carbon density of the vegetation type in the embodiment of the present invention is calculated by the following formula:

[0073] C den,i =a·W i

[0074] Among them, C den,i represents the vegetation carbon density in the i-th land use type; W i represents the actual biomass of the i-th land use type; a represents the conversion factor, which can be determined by direct measurement in the laboratory or literature research.

[0075] Step S4: Obtain the area of ​​different vegetation types and combine the carbon density of different vegetation types to obtain the regional vegetation carbon storage.

[0076] In practical applications, the area of ​​different vegetation types can be obtained by combining remote sensing and field research, and then the regional vegetation carbon storage can be obtained based on the carbon density of different vegetation types. The regional vegetation carbon storage is calculated by the following formula:

[0077]

[0078] Where: C stor Represents the regional vegetation carbon storage; A i represents the area of ​​the i-th vegetation type; N represents the number of vegetation types.

[0079] The quantitative simulation and assessment method for regional vegetation carbon storage provided by the embodiment of the present invention first calculates the potential biomass growth of plants based on the plant leaf area index, further calculates the actual biomass growth of plants by considering the stress of conditions such as water and temperature, and then calculates the vegetation carbon density using the biomass conversion factor method. Finally, based on a method combining remote sensing and field research, the area of ​​different vegetation types is obtained, and then the regional vegetation carbon storage is obtained based on the carbon density of different vegetation types. This method has a strong physical mechanism and can realize rapid simulation and assessment of surface vegetation carbon storage on a daily scale, in a large area, and for a long time.

[0080] Example 2

[0081] The embodiment of the present invention provides a regional vegetation carbon storage quantitative simulation evaluation system, such as Figure 2 As shown, including:

[0082] The plant potential biomass growth calculation module 1 is used to calculate the plant potential biomass growth based on the plant leaf area index; this module executes the method described in step S1 of embodiment 1, which will not be repeated here.

[0083] The plant actual biomass calculation module 2 is used to calculate the plant actual biomass according to the growth of the plant potential biomass and the growth stress conditions during the plant growth process; the module executes the method described in step S2 of embodiment 1, which will not be repeated here.

[0084] The vegetation type carbon density calculation module 3 is used to calculate the carbon density of the specified vegetation type based on the actual biomass of the specified plant type using the biomass conversion factor method. This module executes the method described in step S3 of embodiment 1, which will not be repeated here.

[0085] The regional vegetation carbon storage calculation module 4 is used to obtain the area of ​​different vegetation types and the carbon density of different vegetation types to obtain the regional vegetation carbon storage. This module executes the method described in step S4 of embodiment 1, which will not be described in detail here.

[0086] Example 3

[0087] An embodiment of the present invention provides a computer device, such as Figure 3As shown, it includes: at least one processor 401, such as a CPU (Central Processing Unit), at least one communication interface 403, a memory 404, and at least one communication bus 402. Among them, the communication bus 402 is used to realize the connection and communication between these components. Among them, the communication interface 403 may include a display screen (Display), a keyboard (Keyboard), and the optional communication interface 403 may also include a standard wired interface and a wireless interface. The memory 404 may be a high-speed RAM memory (Ramdom Access Memory, volatile random access memory) or a non-volatile memory (non-volatile memory), such as at least one disk memory. The memory 404 may also be at least one storage device located away from the aforementioned processor 401. Among them, the processor 401 can execute the quantitative simulation and assessment method of regional vegetation carbon storage in Example 1. A set of program codes are stored in the memory 404, and the processor 401 calls the program code stored in the memory 404 to execute the quantitative simulation and assessment method of regional vegetation carbon storage in Example 1.

[0088] The communication bus 402 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus. The communication bus 402 may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 3 Although only one line is used in the diagram, it does not mean that there is only one bus or only one type of bus.

[0089] Among them, the memory 404 may include a volatile memory (English: volatile memory), such as a random access memory (English: random-access memory, abbreviated: RAM); the memory may also include a non-volatile memory (English: non-volatile memory), such as a flash memory (English: flash memory), a hard disk drive (English: hard disk drive, abbreviated: HDD) or a solid-state drive (English: solid-state drive, abbreviated: SSD); the memory 404 may also include a combination of the above types of memory.

[0090] The processor 401 may be a central processing unit (CPU), a network processor (NP), or a combination of a CPU and a NP.

[0091] The processor 401 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The PLD may be a complex programmable logic device (CPLD), a field-programmable gate array (FPGA), a generic array logic (GAL), or any combination thereof.

[0092] Optionally, the memory 404 is also used to store program instructions. The processor 401 can call the program instructions to implement the quantitative simulation and assessment method for regional vegetation carbon storage in Example 1 of the present application.

[0093] The embodiment of the present invention further provides a computer-readable storage medium, on which computer-executable instructions are stored, and the computer-executable instructions can execute the quantitative simulation assessment method of regional vegetation carbon storage in embodiment 1. The storage medium can be a magnetic disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk (HDD) or a solid-state drive (SSD), etc.; the storage medium can also include a combination of the above-mentioned types of memory.

[0094] Obviously, the above embodiments are merely examples for the purpose of clear explanation, and are not intended to limit the implementation methods. For those skilled in the art, other different forms of changes or modifications can be made based on the above description. It is not necessary and impossible to list all the implementation methods here. The obvious changes or modifications derived from these are still within the scope of protection of the invention.

Claims

1. A quantitative simulation and assessment method for regional vegetation carbon storage, characterized in that: include: The growth amount of the plant potential biomass was calculated based on the leaf area index of the plant, and the growth amount of the plant potential biomass was calculated by the following formula: Among them, P bio represents the total potential biomass growth of plants during the calculation period; N represents the total number of simulation days; △P bio，i represents the potential biomass growth of the plant on day i, calculated by the following formula: ΔP bio =R·H e Among them, △P bio It indicates the potential biomass growth of the plant on that day; R indicates the light energy utilization coefficient of the plant, which is related to the carbon dioxide concentration in the air and is calculated based on the carbon dioxide concentration in the air; H e It represents the daily effective photosynthesis solar radiation, which is calculated by the following formula: H e =0.5·H d ·(1-exp·(-k l ·LAI+0.05)) Among them, H d Represents the total daily solar radiation; K l is the extinction coefficient, which indicates the light absorption capacity of plants; LAI indicates the leaf area index of the plant on the simulated day; Calculate the actual biomass of plants according to the growth of potential biomass of plants and the growth stress conditions during plant growth; Based on the actual biomass of the specified plant type, the carbon density of the specified vegetation type is calculated using the biomass conversion factor method; The area of ​​different vegetation types was obtained and combined with the carbon density of different vegetation types to obtain the regional vegetation carbon stock.

2. The method for quantitative simulation and assessment of regional vegetation carbon storage according to claim 1, characterized in that: The leaf area index of a plant is related to the ratio of the plant's accumulated heat units to the plant's potential heat units.

3. The method for quantitative simulation and assessment of regional vegetation carbon storage according to claim 2, characterized in that: The leaf area index of the plant is calculated by the following formula: LAI n =f(pro n , PHU represents the total cumulative heat unit required for plant maturity; n It indicates the ratio of plant potential heat units on the nth day, which refers to the ratio of the sum of daily cumulative heat units within n days of plant growth to the total cumulative heat units required for plant maturity; LAI n represents the leaf area index of the plant on day n, and f(·) represents the correlation function between the leaf area index of the plant and the ratio of the plant's potential heat unit.

4. The method for quantitative simulation and assessment of regional vegetation carbon storage according to claim 1, characterized in that: The actual growth of the plant is calculated by the following formula: Among them, P bio,act Represents the actual biomass of plants; P bio,ori represents the initial biomass of the plant; r cf represents the plant growth stress coefficient; wst represents the water stress factor, which is calculated by the actual plant evapotranspiration and the plant potential evapotranspiration; tst represents the temperature stress factor, which is a function of the daily average temperature and the optimal growth temperature of the plant.

5. The method for quantitative simulation and assessment of regional vegetation carbon storage according to claim 1, characterized in that: The carbon density of vegetation types is calculated by the following formula: C den,i =a·W i Among them, C den,i represents the vegetation carbon density in the i-th land use type; W i represents the actual biomass of the i-th land use type; a represents the conversion factor.

6. The method for quantitative simulation and assessment of regional vegetation carbon storage according to claim 5, characterized in that: The carbon storage of vegetation in the region is calculated using the following formula: Among them, C stor Represents the regional vegetation carbon storage; A i represents the area of ​​the i-th vegetation type; N represents the number of vegetation types.

7. A regional vegetation carbon storage quantitative simulation and assessment system, characterized in that: include: The plant potential biomass growth calculation module is used to calculate the plant potential biomass growth based on the plant leaf area index. The plant potential biomass growth is calculated by the following formula: Among them, P bio represents the total potential biomass growth of plants during the calculation period; N represents the total number of simulation days; △P bio，i represents the potential biomass growth of the plant on day i, calculated by the following formula: ΔP bio =R·H e Among them, △P bio It indicates the potential biomass growth of the plant on that day; R indicates the light energy utilization coefficient of the plant, which is related to the carbon dioxide concentration in the air and is calculated based on the carbon dioxide concentration in the air; H e It represents the daily effective photosynthesis solar radiation, which is calculated by the following formula: H e =0.5·H d ·(1-exp·(-k l ·LAI+0.05)) Among them, H d Represents the total daily solar radiation; K l is the extinction coefficient, which indicates the light absorption capacity of plants; LAI indicates the leaf area index of the plant on the simulated day; The plant actual biomass calculation module is used to calculate the plant actual biomass according to the growth of the plant potential biomass and the growth stress conditions during the plant growth process; The carbon density calculation module of vegetation type is used to calculate the carbon density of the specified vegetation type based on the actual biomass of the specified plant type and using the biomass conversion factor method; The regional vegetation carbon storage calculation module is used to obtain the area of ​​different vegetation types and combine the carbon density of different vegetation types to obtain the regional vegetation carbon storage.

8. A computer device, characterized in that: include: At least one processor, and a memory communicatively connected to the at least one processor, wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor executes the method for quantitative simulation and assessment of regional vegetation carbon storage as described in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable the computer to execute the method for quantitative simulation and assessment of regional vegetation carbon storage as described in any one of claims 1-6.

Images