CN113591032B - Method for calculating optimal growth elevation of mangrove plant - Google Patents
Method for calculating optimal growth elevation of mangrove plant Download PDFInfo
- Publication number
- CN113591032B CN113591032B CN202110798731.5A CN202110798731A CN113591032B CN 113591032 B CN113591032 B CN 113591032B CN 202110798731 A CN202110798731 A CN 202110798731A CN 113591032 B CN113591032 B CN 113591032B
- Authority
- CN
- China
- Prior art keywords
- time
- elevation
- point
- tidal
- water level
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 240000002044 Rhizophora apiculata Species 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 20
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 65
- 230000000694 effects Effects 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 230000004083 survival effect Effects 0.000 claims abstract description 7
- 102000004190 Enzymes Human genes 0.000 claims abstract description 6
- 108090000790 Enzymes Proteins 0.000 claims abstract description 6
- 230000000243 photosynthetic effect Effects 0.000 claims abstract description 6
- 238000005259 measurement Methods 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 46
- 238000004364 calculation method Methods 0.000 claims description 34
- 239000013598 vector Substances 0.000 claims description 22
- 230000000630 rising effect Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 4
- 241000120622 Rhizophoraceae Species 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 241001635619 Avicennia marina Species 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 240000003793 Rhizophora mangle Species 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000002028 Biomass Substances 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 241001366278 Leptotes marina Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Cultivation Of Plants (AREA)
Abstract
The invention discloses a method for calculating the optimal growth elevation of mangrove plants, which comprises the following steps: receiving a user instruction, setting a flooding time gradient to perform an indoor flooding simulation test, and determining the optimal flooding time of the mangrove plant according to the survival rate, the growth state, the photosynthetic efficiency and the enzyme activity of the mangrove plant; receiving tidal level data collected by nearby tidal stations or field measurements entered by a user; after the two data are obtained, determining the flooding time under each elevation by calculating the total time when the tidal water level is higher than each elevation of the coastal zone, and comparing the flooding time with the optimal flooding time of the mangrove plant, wherein the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant. The method can be used in the mangrove ecological restoration field, has important significance for hydrologic restoration, saves time and cost compared with the method for determining the optimal growth elevation of the mangrove by the traditional field test, and has universal applicability.
Description
Technical Field
The invention relates to the field of mangrove plants, in particular to a method for calculating the optimal growth height of a mangrove plant based on the optimal flooding time and the tidal water level of the mangrove plant.
Background
Mangrove is a evergreen woody plant grown in tropical and subtropical coastal zones, and is mainly distributed in coastal beach, swamps, tidal ditches and estuary regions between 30 degrees in north and south latitude. Mangrove has important ecological system service value, such as providing habitat for coastal marine organisms, wind prevention, shore fixation, water purification, important places as marine blue carbon sink, and the like. However, the global mangrove area continues to decrease and the chinese mangrove area undergoes a rapid decrease to a slow increase. Numerous mangrove restoration practices show that hydrologic conditions (such as flooding time and the like) have important influence on the growth and propagation of mangrove plants, and too long or too short flooding time can influence the survival rate of the mangrove and influence the mangrove restoration effect. The flooding time is closely related to the elevation of the coastal zone, and mangrove repair activities can be carried out after the suitable forest areas of the mangrove are determined or the elevation is modified by determining Gao Chenglai the optimal growth of the mangrove plant.
Along with the continuous improvement of ecological environmental awareness of people and the propaganda investment of government and NGO organizations, for politics, economy, ecological protection and other reasons, mangrove forestation activities in China are numerous and mainly concentrated on the beach. The biggest problem of pond withdrawal and forest withdrawal is that the ecological environment of the culture pond is improved, particularly the hydrologic transformation, which is closely related to the elevation, so that the optimum growth elevation of mangrove plants needs to be known. At present, the optimal growth elevation of mangrove plants is mainly determined through field tests, the period is long, the cost is high, and field test results are not universal due to different tides and terrains in various places. However, the current low-cost and reproducible method for determining the optimal growth elevation of mangrove plants in China is far from mature.
Disclosure of Invention
In order to overcome the defects of long period, high cost and irreproducibility of the prior art, the invention provides a method for calculating the optimal growth elevation of mangrove plants.
The technical problems of the invention are solved by the following technical scheme:
the invention provides a method for calculating the optimal growth elevation of mangrove plants, which comprises the following steps: s1: receiving a user instruction, setting a flooding time gradient to perform an indoor flooding simulation test, and determining the optimal flooding time of the mangrove plant according to the survival rate, the growth state, the photosynthetic efficiency and the enzyme activity of the mangrove plant; s2: receiving tidal level data collected by nearby tidal stations or field measurements entered by a user; s3: after the two data are obtained, determining the flooding time under each elevation by calculating the total time when the tidal water level is higher than each elevation of the coastal zone, and comparing the flooding time with the optimal flooding time of the mangrove plant, wherein the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant.
In some embodiments, the setting the flooding time gradient in the step S1 for performing the indoor flooding simulation test includes: setting the artificial simulation wetland with the daily flooding time as a certain gradient by using a water pump and a timer, and planting corresponding mangroves.
In some embodiments, step S3 comprises the following specific steps: s31: performing function fitting on the tidal water level data for calculation, adopting cubic spline curve fitting and linear fitting according to tidal collection time and actual tidal conditions, wherein the fitted function is a piecewise function, and calculating the time point of the tidal water level at a certain elevation and the total flooding time under the elevation through the fitted function; s32: establishing a cycle taking 1cm as step length from the lowest tide level to the highest tide level, establishing a cycle body inside the cycle body, calculating the highest water level and the lowest water level of the tide piecewise function in the time interval, and determining by calculating and comparing the end point value of the function and the extreme value inside the function interval; s33: judging whether the elevation in the circulation is between the lowest water level and the highest water level of the piecewise function; if so, calculating the intersection point of the elevation and the tidal piecewise function, determining the time of the tidal water level at the elevation, solving the derivative at the intersection point, and carrying out the next tidal piecewise function cycle; if not, directly carrying out the next tide piecewise function circulation; s34: obtaining the intersection point of all elevations and the tidal function between the lowest tide level and the highest tide level and the derivative of the intersection point through the small cycle calculated by the intersection point of the tidal piecewise function and the elevations and the large cycle of the external elevations; three vectors corresponding one to one are obtained: the elevation vector, the time intersection point vector and the intersection point derivative vector are used for removing the corresponding vector element with the derivative of 0; s35: after obtaining the flooding time point of each elevation and the derivative of the flooding time point, calculating the total annual flooding time of each elevation, namely the sum of the time when the tidal water level is higher than a certain elevation; the time interval when the derivative of the intersection point is changed from positive to negative is the time when the tidal water level is higher than a certain elevation; the derivative of the first intersection point of the elevation and the tidal function has positive or negative conditions, and the total number of intersection points of the elevation and the tidal function can be divided into singular numbers and even numbers; respectively calculating according to the positive and negative of the first intersection derivative and the single and double division of the total number of the intersections; s36: comparing the calculated flooding time under all elevations with the optimal flooding time of the mangrove plant, and when the flooding time under a certain elevation is consistent with the optimal flooding time of the mangrove plant, obtaining the elevation which is the optimal growth elevation of the mangrove plant.
In some embodiments, in the step S34, the elevation vector is an elevation of the lowest level to the highest level of the tide at intervals of 1cm, the time intersection vector is a time point of the tide level at the corresponding elevation, the intersection derivative vector is a value capable of judging rise or fall of the tide level according to signs therein, and positive indicates rise of the tide level and negative indicates fall of the tide level.
In some embodiments, in the step S35, the calculating according to the positive and negative of the first intersection derivative and the single-double division of the total number of intersections includes: case 1: when the derivative of the first intersection point is smaller than 0 and the total number of intersection points is singular, the total time calculation point is the intersection point plus the initial time; case 2: when the first intersection derivative is smaller than 0 and the total number of the intersections is double, the total time calculation point is the intersection plus an initial time point and an end time point; case 3: when the first intersection derivative is greater than 0 and the total number of intersections is singular, the total time calculation point is the intersection plus the ending time point; case 4: when the first intersection derivative is greater than 0 and the total number of intersections is a double number, the total time calculation is based on the intersection calculation alone.
In some embodiments, in the case 1, when the derivative of the first intersection point is less than 0 and the total number of intersection points is singular, it is indicated that the tidal water level is reduced to a certain elevation at the initial time point, and the tidal water level is continuously reduced at a certain elevation at the end time point, and the total time calculation point is the intersection point plus the initial time point, and the total flooding time of a certain elevation is obtained according to the sum of the adjacent time points.
In some embodiments, in case 2, when the derivative of the first intersection is less than 0 and the total number of intersections is two, it is indicated that the tidal level is lowered to a certain elevation at the initial time point, the tidal level is at a certain Gao Chengshang liter at the end, and the total time calculation point is the intersection plus the initial time point and the end time point, and the total flooding time of a certain elevation is obtained according to the sum of the differences between the adjacent time points.
In some embodiments, in the case 3, when the derivative of the first intersection is greater than 0 and the total number of intersections is singular, it is indicated that the tidal level rises to a certain elevation at the initial time point, and the tidal level continues to rise at a certain elevation at the end time point, and the total time calculation point is the intersection plus the end time point, and the total flooding time of a certain elevation is obtained according to the sum of the differences between the adjacent time points.
In some embodiments, in case 4, when the first intersection derivative is greater than 0 and the total number of intersections is a double number, it is stated that the tidal level rises to a certain elevation at the initial point in time, the tidal level falls at a certain elevation at the end, the total time calculation is based on the intersection alone, and the total flooding time for a certain elevation is obtained based on the sum of the differences between adjacent points in time.
The invention also provides a terminal device of the mangrove plant optimum growth elevation, which comprises a memory, a processor and a computer program stored in the memory and capable of running on the processor, and is characterized in that the processor realizes the steps of any one of the methods when executing the computer program.
Compared with the prior art, the invention has the beneficial effects that: according to the method, a calculation method for calculating the optimal growth elevation of the mangrove plant is developed according to the optimal flooding time of the mangrove plant and tidal water level data, an indoor flooding simulation test is carried out by setting a flooding time gradient, the flooding time at each elevation is determined by calculating the total time of the tidal water level higher than each elevation of a coastal zone, and compared with the optimal flooding time of the mangrove plant, the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant; the method can be used in the mangrove ecological restoration field, has important significance in hydrologic restoration, saves time and cost compared with the method for determining the most suitable growth elevation of the mangrove by the traditional field test, and has universal applicability; the method can effectively guide the elevation hydrologic transformation when the mangrove is restored by destroying the habitat, thereby improving the mangrove restoration effect, maintaining the coastal zone ecological system and ensuring the sustainable development of the economy and society.
Drawings
FIG. 1 is a flow chart of calculation of flooding time at different elevations according to an embodiment of the present invention;
FIG. 2 is a schematic diagram I of the flooding time calculation according to an embodiment of the present invention;
FIG. 3 is a schematic diagram II of flooding time calculation according to an embodiment of the present invention;
FIG. 4 is a schematic diagram III of flooding time calculation according to an embodiment of the present invention;
fig. 5 is a schematic diagram iv of flooding time calculation according to an embodiment of the present invention.
Detailed Description
The invention will be further described with reference to the following drawings in conjunction with the preferred embodiments. It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other.
It should be noted that, in this embodiment, the terms of left, right, upper, lower, top, bottom, etc. are merely relative terms, or refer to the normal use state of the product, and should not be considered as limiting.
The invention aims to provide a low-cost, large-scale and reproducible method for determining the optimal growth elevation of mangrove plants so as to ensure the effect of mangrove restoration activities. In order to solve the technical problems, the invention adopts the following technical scheme:
the method for calculating the optimal growth elevation of the mangrove plant comprises the following steps of: s1: receiving a user instruction, setting a flooding time gradient to perform an indoor flooding simulation test, and determining the optimal flooding time of the mangrove plant according to the survival rate, the growth state, the photosynthetic efficiency and the enzyme activity of the mangrove plant; s2: receiving tidal level data collected by nearby tidal stations or field measurements entered by a user; s3: after the two data are obtained, determining the flooding time under each elevation by calculating the total time when the tidal water level is higher than each elevation of the coastal zone, and comparing the flooding time with the optimal flooding time of the mangrove plant, wherein the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant.
According to the method for calculating the optimal growth elevation of the mangrove plant, the indoor flooding simulation test is carried out by setting the flooding time gradient, for example, a water pump, a timer and the like are used for setting the flooding time of each day to be 0h, 2h, 4h, 6h, 8h, 10h, 12h, 14h, 16h, 18h, 20h, 22h and 24h, corresponding mangrove is planted in the artificial simulation wetland, and the optimal flooding time of the mangrove plant is determined according to the survival rate, growth state (plant height, base diameter, biomass and the like), photosynthetic efficiency, enzyme activity and the like of the mangrove plant, for example, when the flooding time of each day is 9-12 h, the survival rate of the plant is highest when the flooding time of each day is 6-8 h, and the growth state, photosynthetic efficiency, enzyme activity and the like of the mangrove plant are optimal. Tidal water level data is collected by nearby tide stations or field measurements.
After the two kinds of data are obtained, calculating the flooding time under each elevation of the coastal zone, determining the flooding time under each elevation by calculating the total time when the tidal water level is higher than each elevation of the coastal zone, and comparing the flooding time with the optimal flooding time of the mangrove plant, wherein the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant. The calculation flow chart of the flooding time at different elevations is shown in figure 1.
And carrying out function fitting on the tidal water level data for calculation, carrying out spline curve fitting for three times or linear fitting for an interp1 function by adopting a MATLAB built-in function spline function according to the tidal collection time and the actual tidal condition, wherein the fitted function is a piecewise function, and calculating the time point of the tidal water level at a certain elevation and the total flooding time under the elevation through the fitted function.
A cycle is established from the lowest tide level to the highest tide level in 1cm steps, then a cycle body is established in the cycle body, and the highest water level and the lowest water level of the tide segmentation function in the time interval are calculated (determined by comparing the end value of the function with the extreme value in the function interval).
Then, judging whether the elevation in the circulation is between the lowest water level and the highest water level of the piecewise function, if so, calculating the intersection point of the elevation and the tidal piecewise function, determining the time of the tidal water level at the elevation, solving the derivative at the intersection point, and carrying out the next tidal piecewise function circulation; if not, the next tidal piecewise function cycle is performed directly.
The intersection point of all elevations between the lowest tide level and the highest tide level and the tide function and the derivative of the intersection point are obtained through small cycles calculated by intersection points of the tide sectional functions and the elevations and large cycles of the external elevations, the intersection point is the time point of the tide level at a certain elevation, and the intersection point derivative can judge whether the tide level rises or falls at the time point; three vectors corresponding one by one are obtained, an elevation vector (the element is the elevation of the lowest water level to the highest water level of the tide at intervals of 1 cm), a time intersection vector (the element is the time point of the tide level at the corresponding elevation), and an intersection derivative vector (the tide level rising or falling can be judged according to the sign of the element, the element is positive and represents the tide level rising, the element is negative and represents the tide level falling), and the corresponding vector element with the derivative of 0 is removed.
After the flooding time point and the derivative of the time point of each elevation are obtained, the total annual flooding time of each elevation is calculated, namely the sum of the time when the tidal water level is higher than a certain elevation. When the derivative at the intersection point is positive, the tidal level is rising; when the derivative at the intersection is negative, it is indicated that the tidal level is falling, so the time interval in which the derivative changes from positive to negative at the intersection is the time when the tidal level is above a certain elevation. The programming is mainly divided into 4 calculation cases.
When the derivative of the first intersection point is smaller than 0 and the total number of intersection points is singular, the total time calculation point is the intersection point plus the initial time; it is explained that the tidal water level is lowered to a certain elevation at the initial time point, and the tidal water level is continuously lowered at the certain elevation at the end, so that the total time calculation point is the intersection point plus the initial time point, and then the total flooding time of the certain elevation is obtained according to the sum of the differences between the adjacent time points (even points minus odd points), and the total flooding time is shown as a shaded part of fig. 2.
When the first intersection derivative is smaller than 0 and the total number of the intersections is double, the total time calculation point is the intersection plus an initial time point and an end time point; it is explained that the tidal level drops to a certain elevation at the initial time point, the tidal level is at a certain Gao Chengshang liter at the end, the total time calculation point is the intersection point plus the initial time point and the end time point, and then the total flooding time of a certain elevation is obtained according to the sum of the differences between the adjacent time points (even points minus odd points), and the total flooding time is shown as the shaded part of fig. 3.
When the first intersection derivative is greater than 0 and the total number of intersections is singular, the total time calculation point is the intersection plus the ending time point; the tidal water level rises to a certain elevation at the initial time point, the tidal water level continues to rise at the certain elevation at the end time point, the total time calculation point is the intersection point plus the end time point, and then the total flooding time of the certain elevation is obtained according to the sum of the differences between the adjacent time points (even points minus odd points), and the total flooding time is shown in a shadow part of fig. 4.
When the first intersection derivative is greater than 0 and the total number of intersections is double, calculating the total time only according to the intersections; it is explained that at the initial time point the tidal level rises to a certain elevation and at the end the tidal level falls at a certain elevation, the total time calculation is calculated only from the intersections, and then the total flooding time for a certain elevation is obtained from the sum of the differences between adjacent time points (even points minus odd points), the total flooding time being shown in the shaded portion of fig. 5.
And finally, comparing the calculated flooding time under all elevations with the optimal flooding time of the mangrove plant, and when the flooding time under a certain elevation is consistent with the optimal flooding time of the mangrove plant, obtaining the elevation which is the optimal growth elevation of the mangrove plant.
According to the invention, the optimal flooding time of the avicennia marina is 8-12 h per day, the optimal flooding time is converted into the annual optimal flooding time, namely, the annual optimal flooding time of the avicennia marina is 2920-4280 h, then tidal water level data of a Zhanjiang harbor tide station of the peninsula, is collected, and finally, the optimal growth elevation of the avicennia marina near the Zhanjiang harbor of the peninsula, which is the Zhanzhanjiang, is 216-259 cm through MATLAB programming calculation.
According to the invention, an MATLAB program algorithm for calculating the optimal growth elevation of the mangrove plant is developed according to the optimal flooding time and tidal water level data of the mangrove plant, and the elevation hydrologic transformation during the mangrove forest restoration in the damaged habitat can be effectively guided, so that the mangrove forest restoration effect is improved, the coastal zone ecological system is maintained, and the sustainable development of the economy and society is ensured.
The foregoing is a further detailed description of the invention in connection with the preferred embodiments, and it is not intended that the invention be limited to the specific embodiments described. It will be apparent to those skilled in the art that several equivalent substitutions and obvious modifications can be made without departing from the spirit of the invention, and the same should be considered to be within the scope of the invention.
Claims (3)
1. The method for calculating the optimal growth elevation of the mangrove plant is characterized by comprising the following steps of:
s1: receiving a user instruction, setting a flooding time gradient to perform an indoor flooding simulation test, and determining the optimal flooding time of the mangrove plant according to the survival rate, the growth state, the photosynthetic efficiency and the enzyme activity of the mangrove plant;
s2: receiving tidal level data collected by nearby tidal stations or field measurements entered by a user;
s3: after the two data are obtained, determining the flooding time under each elevation by calculating the total time when the tidal water level is higher than each elevation of the coastal zone, and comparing the flooding time with the optimal flooding time of the mangrove plant, wherein the elevation meeting the optimal flooding time of the mangrove plant is taken as the optimal growth elevation of the mangrove plant;
wherein, the step S3 comprises the following specific steps:
s31: performing function fitting on the tidal water level data for calculation, adopting cubic spline curve fitting and linear fitting according to tidal collection time and actual tidal conditions, wherein the fitted function is a piecewise function, and calculating the time point of the tidal water level at a certain elevation and the total flooding time under the elevation through the fitted function;
s32: establishing a cycle taking 1cm as step length from the lowest tide level to the highest tide level, establishing a cycle body inside the cycle body, calculating the highest water level and the lowest water level of the tide piecewise function in the time interval, and determining by calculating and comparing the end point value of the function and the extreme value inside the function interval;
s33: judging whether the elevation in the circulation is between the lowest water level and the highest water level of the piecewise function;
if so, calculating the intersection point of the elevation and the tidal piecewise function, determining the time of the tidal water level at the elevation, solving the derivative at the intersection point, and carrying out the next tidal piecewise function cycle;
if not, directly carrying out the next tide piecewise function circulation;
s34: obtaining the intersection point of all elevations and the tidal function between the lowest tide level and the highest tide level and the derivative of the intersection point through the small cycle calculated by the intersection point of the tidal piecewise function and the elevations and the large cycle of the external elevations; three vectors corresponding one to one are obtained: the elevation vector, the time intersection point vector and the intersection point derivative vector are used for removing the corresponding vector element with the derivative of 0;
the elevation vector is the elevation of the lowest water level to the highest water level of the tide at intervals of 1cm, the time intersection point vector is the time point of the tide water level under the corresponding elevation, the intersection point derivative vector is the time point of the tide water level, the rising or the falling of the tide water level can be judged according to the signs in the intersection point derivative vector, the rising of the tide water level is represented positively, and the falling of the tide water level is represented negatively;
s35: after obtaining the flooding time point of each elevation and the derivative of the flooding time point, calculating the total annual flooding time of each elevation, namely the sum of the time when the tidal water level is higher than a certain elevation; the time interval when the derivative of the intersection point is changed from positive to negative is the time when the tidal water level is higher than a certain elevation;
the derivative of the first intersection point of the elevation and the tidal function has positive or negative conditions, and the total number of intersection points of the elevation and the tidal function can be divided into singular numbers and even numbers; respectively calculating according to the positive and negative of the first intersection derivative and the single and double division of the total number of the intersections;
the calculating according to the single-double division of the positive and negative of the first intersection derivative and the total number of the intersections respectively comprises the following steps:
case 1: when the derivative of the first intersection point is smaller than 0 and the total number of intersection points is singular, the total time calculation point is the intersection point plus the initial time;
when the derivative of the first intersection point is smaller than 0 and the total number of intersection points is singular, the tidal water level is indicated to drop to a certain elevation at the initial time point, the tidal water level continues to drop at a certain elevation at the end, the total time calculation point is the intersection point plus the initial time point, and the total flooding time of a certain elevation is obtained according to the sum of the differences of the adjacent time points;
case 2: when the first intersection derivative is smaller than 0 and the total number of the intersections is double, the total time calculation point is the intersection plus an initial time point and an end time point;
when the derivative of the first intersection point is smaller than 0 and the total number of the intersection points is equal to two, the tidal water level is reduced to a certain elevation at the initial time point, the tidal water level is at a certain Gao Chengshang liter when the tidal water level is ended, the total time calculation point is the intersection point plus the initial time point and the ending time point, and the total flooding time of a certain elevation is obtained according to the sum of the differences of the adjacent time points;
case 3: when the first intersection derivative is greater than 0 and the total number of intersections is singular, the total time calculation point is the intersection plus the ending time point;
when the derivative of the first intersection point is larger than 0 and the total number of intersection points is singular, the tidal water level rises to a certain elevation at an initial time point, the tidal water level continues to rise at a certain elevation at the end, the total time calculation point is the intersection point plus the end time point, and the total flooding time of a certain elevation is obtained according to the sum of the differences of adjacent time points;
case 4: when the first intersection derivative is greater than 0 and the total number of intersections is double, calculating the total time only according to the intersections;
when the derivative of the first intersection point is larger than 0 and the total number of intersection points is double, the tidal water level rises to a certain elevation at the initial time point, the tidal water level falls at a certain elevation at the end, the total time is calculated only according to the intersection points, and the total flooding time of a certain elevation is obtained according to the sum of the differences of the adjacent time points;
s36: comparing the calculated flooding time under all elevations with the optimal flooding time of the mangrove plant, and when the flooding time under a certain elevation is consistent with the optimal flooding time of the mangrove plant, obtaining the elevation which is the optimal growth elevation of the mangrove plant.
2. The method for calculating the optimal growth height of mangrove plants according to claim 1, wherein the step S1 of setting the flooding time gradient for the indoor flooding simulation test comprises: setting the artificial simulation wetland with the daily flooding time as a certain gradient by using a water pump and a timer, and planting corresponding mangroves.
3. Terminal device for the optimal growth of mangrove plants, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any one of claims 1-2 when the computer program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110798731.5A CN113591032B (en) | 2021-07-15 | 2021-07-15 | Method for calculating optimal growth elevation of mangrove plant |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110798731.5A CN113591032B (en) | 2021-07-15 | 2021-07-15 | Method for calculating optimal growth elevation of mangrove plant |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113591032A CN113591032A (en) | 2021-11-02 |
| CN113591032B true CN113591032B (en) | 2023-06-16 |
Family
ID=78247393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110798731.5A Active CN113591032B (en) | 2021-07-15 | 2021-07-15 | Method for calculating optimal growth elevation of mangrove plant |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113591032B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116168015B (en) * | 2023-03-21 | 2024-03-26 | 国家海洋信息中心 | Identification method and device for mangrove repair suitable area |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102619192A (en) * | 2012-04-06 | 2012-08-01 | 广西红树林研究中心 | Method for increasing elevation of intertidal zone to form suitable land for mangrove forest |
| CN111275254A (en) * | 2020-01-17 | 2020-06-12 | 河海大学 | Method for predicting development trend of future intertidal zone area based on intertidal zone area change model |
| CN111402169A (en) * | 2020-03-23 | 2020-07-10 | 宁波大学 | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide |
| CN112307420A (en) * | 2020-10-29 | 2021-02-02 | 大连海洋大学 | Estuary wetland vegetation habitat restoration method |
| CN113011716A (en) * | 2021-02-25 | 2021-06-22 | 长江水资源保护科学研究所 | Method and system for evaluating influence of lake water level change on functional characters of aquatic plants |
| CN113029290A (en) * | 2021-02-26 | 2021-06-25 | 澜途集思生态科技集团有限公司 | Method for measuring appropriate water level of aquatic plant |
-
2021
- 2021-07-15 CN CN202110798731.5A patent/CN113591032B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102619192A (en) * | 2012-04-06 | 2012-08-01 | 广西红树林研究中心 | Method for increasing elevation of intertidal zone to form suitable land for mangrove forest |
| CN111275254A (en) * | 2020-01-17 | 2020-06-12 | 河海大学 | Method for predicting development trend of future intertidal zone area based on intertidal zone area change model |
| CN111402169A (en) * | 2020-03-23 | 2020-07-10 | 宁波大学 | Method for repairing remote sensing vegetation index time sequence under influence of coastal tide |
| CN112307420A (en) * | 2020-10-29 | 2021-02-02 | 大连海洋大学 | Estuary wetland vegetation habitat restoration method |
| CN113011716A (en) * | 2021-02-25 | 2021-06-22 | 长江水资源保护科学研究所 | Method and system for evaluating influence of lake water level change on functional characters of aquatic plants |
| CN113029290A (en) * | 2021-02-26 | 2021-06-25 | 澜途集思生态科技集团有限公司 | Method for measuring appropriate water level of aquatic plant |
Non-Patent Citations (2)
| Title |
|---|
| 《effects of exotic and native mangrove forests plantation on soil organic catbon,nitrogen,and phosphorus contents and pools in leizhou,china》;jianxiang feng etal;《CATENA》;全文 * |
| 《mariculture pond influence on mangrove areas in south china:significantly larger nitrogen and phosphorus loadings from sediment wash-out than from tidal water exchange》;Hao Wu etal;《Aquaculture》;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113591032A (en) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113591032B (en) | Method for calculating optimal growth elevation of mangrove plant | |
| Wu et al. | Quantifying the dynamics and driving forces of the coastal wetland landscape of the Yangtze River Estuary since the 1960s | |
| Anderson et al. | Sources and fates of nitrogen in Virginia coastal bays | |
| CN206014519U (en) | Multi-functional water field planting chinampa | |
| CN204548429U (en) | Multifunctional floating island base on a kind of Green Water | |
| CN111567301B (en) | Method for planting mangrove forest in deep water bare beach | |
| CN204350790U (en) | A kind of integrated form expanded matrix stereoscopic strawberry cultivation and management devices | |
| CN204272940U (en) | A kind of staged mountain spring sturgeon culturing pond | |
| Srinivas et al. | Atmospheric forcing on the seasonal variability of sea level at Cochin, southwest coast of India | |
| CN210885524U (en) | Wind-wave-resistant ecological floating island | |
| CN210537668U (en) | Submerged plant planting module for ecological restoration | |
| CN107907627B (en) | Controlled experimental device for simulating habitat of fish on mudflat and construction method | |
| CN102919113A (en) | Float planting device and float planting field | |
| CN205922200U (en) | Aquatic plant circulating water experimental tank | |
| CN105660358A (en) | Alternate utilization method of porphyra haitanensis shell conchocelis | |
| US20080005959A1 (en) | Method for cultivating seaweed having adherence | |
| CN207046947U (en) | A kind of high-efficient purifier using emergent aquactic plant | |
| CN205556246U (en) | Water ecological remediation mobile solar power chinampa | |
| CN206042934U (en) | A kind of solid growth box of wet land protection | |
| CN112219594B (en) | Submerged plant planting device and river channel treatment method in ecological restoration system engineering | |
| CN104276662A (en) | Method for purifying eutrophic landscape water body on planted floating bed | |
| CN209594470U (en) | A kind of suspend plants the planting unit of leaf vegetables | |
| Do | Application of GIS to evaluate saltwater intrusion on agricultural production in the Mekong Delta | |
| CN215872800U (en) | Underwater planting disc | |
| CN218043016U (en) | Mangrove floating plate device of growing seedlings |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |