METHOD FOR DETERMINING DYNAMIC WETLAND BOUNDARY BASED ON HYDROLOGY, ORGANISM AND SOIL ELEMENTS

Abstract

The invention discloses a method for determining dynamic wetland boundary based on hydrology, organism and soil elements, including the following steps: step 1: extracting dynamic wetland hydrology boundary; step 2: obtaining data on wetland vegetation; step 3: obtaining data on wetland soil; step 4: simulating dynamic wetland vegetation boundary and dynamic wetland soil boundary according to the data from step 1 to step 3; step 5: determining dynamic wetland boundary; the invention can reflect the dynamic wetland boundary in multiple dimensions comprehensively and accurately.

Claims

1. A method for determining dynamic wetland boundary based on hydrology, organism and soil elements, wherein the method includes the following steps: step 1: extracting dynamic wetland hydrology boundary; step 2: obtaining data on wetland vegetation; step 3: obtaining data on wetland soil; step 4: simulating dynamic wetland vegetation boundary and dynamic wetland soil boundary according to the data from step 1 to step 3; step 5: determining dynamic wetland boundary.

2. The method for determining dynamic wetland boundary based on hydrology, organism and soil elements according to claim 1, wherein the dynamic wetland hydrology boundary in step 1 includes two boundaries: actual water surface boundary, and animal and plant remain boundary.

3. The method for determining dynamic wetland boundary based on hydrology, organism and soil elements according to claim 1, wherein the step 1 includes the following sub-steps: sub-step 1: conducting on-site investigation by high-precision handheld GPS, or using remote sensing image with a resolution of less than 2 m, so as to obtain the actual water surface boundary; searching for the animal and plant remains through on-site investigation; sub-step 2: tracking and recording the positions of animal and plant remains by high-precision handheld GPS, so as to obtain the animal and plant remain boundary, wherein the animal and plant remain boundary is the farthest hydrology boundary; overlaying the actual water surface boundary and the animal and plant remain boundary, so as to obtain the dynamic wetland hydrology boundary.

4. The method for determining dynamic wetland boundary based on hydrology, organism and soil elements according to claim 1, wherein the step 2 includes the following sub-steps: sub-step 1: laying investigation line transects in the direction perpendicular to the wetland hydrology boundary, with an interval of 100 m-500 m among the line transects, with the water surface boundary taken as the starting point of investigation line transect. sub-step 2: setting three investigation quadrats (1 m×1 m) on the line transects every 10 m, recording the vegetation coverage, the plant species and the numbers and coverage of various plants in the quadrats, and determining the species of hygrophytes and xerophytes thereof; sub-step 3: calculating the importance value of each species, wherein the importance value is the sum of relative coverage, relative density and relative frequency of each species; calculating the importance value proportion of hygrophyte.

5. The method for determining dynamic wetland boundary based on hydrology, organism and soil elements according to claim 1, wherein the step 3 includes the following sub-steps: sub-step 1: in the investigation quadrats, collecting the soil samples within the range of 0 cm-10 cm of the surface layer, weighing the soil samples, then drying them at 105° C., and finally calculating the soil moisture content by using the weight difference; sub-step 2: taking the animal and plant remain boundary, i.e., the farthest hydrology boundary, as a benchmark, and calculating the minimum perpendicular distance of quadrat away from the actual water body boundary: $D={\left[{\left(x_0-x\right)}^2+{\left(y_0-y\right)}^2\right]}^{\frac{1}{2}}{D}_m=\underset{\left(x,y\right)\in R}{\min }D$ Where, (x.sub.0, y.sub.0) and (x, y) represent the decimal latitude and longitude coordinates of pixels of each quadrat and hydrology boundary respectively, R represents the decimal latitude and longitude coordinates set of all pixels of hydrology boundary, D represents the minimum perpendicular distance of quadrat away from the actual water body boundary, D.sub.m represents the quadrat with the minimum perpendicular distance away from the actual water body boundary and with the greatest influence by water body, and the smaller the index D, the closer the hydrological connectivity between the quadrat and the wetland water body; sub-step 3: conducting regression analysis on the minimum perpendicular distance of quadrat away from the actual water body boundary, the hygrophyte proportion and the measured soil moisture content so as to obtain a regression equation, and simulating the wetland boundary by using the geographically weighted regression model in ArcGIS software so as to fit curves of the measured and estimated hygrophyte proportion or soil moisture content; sub-step 4: verifying the simulation effect of the model: when the fitted curve R.sup.2 is greater than or equal to 0.6, the simulation result of wetland boundary is reliable, so that the dynamic wetland vegetation boundary and dynamic wetland soil boundary are obtained; when the fitted curve R.sup.2 is less than 0.6, it is necessary to readjust the model or conduct field investigation again; sub-step 5: considering the wetland hydrology boundary, wetland vegetation boundary and wetland soil boundary comprehensively, and taking the actual water surface boundary as a benchmark, so as to analyze the differences of the three boundaries on the perpendicular line of each pixel of water surface boundary; taking the boundary point with the largest D among the three boundary points as the wetland boundary of this perpendicular line, so as to obtain the dynamic wetland boundary.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0030] FIG. 1 is an operation flow chart of the invention.

[0031] FIG. 2 is a hydrology boundary diagram in Embodiment 1.

[0032] FIG. 3 is an analysis chart of vegetation data dependency in Embodiment 1.

[0033] FIG. 4 is an analysis chart of soil data dependency in Embodiment 1.

[0034] FIG. 5 is a scatter diagram of vegetation boundary verification results in Embodiment 1.

[0035] FIG. 6 is a scatter diagram of soil boundary verification results in Embodiment 1.

[0036] FIG. 7 is a vegetation boundary diagram in Embodiment 1.

[0037] FIG. 8 is a soil boundary diagram in Embodiment 1.

[0038] FIG. 9 is an overlay chart of hydrology boundary, vegetation boundary and soil boundary in Embodiment 1.

[0039] FIG. 10 is a dynamic wetland boundary diagram in Embodiment 1.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0040] The detailed embodiment of the invention will be described below in order to facilitate those skilled in the technical field to understand the invention. However, it should be clear that the invention is not limited to the scope of the detailed embodiment. For those skilled in the technical field, if various changes are within the spirit and scope of the invention defined and determined by the appended claims, and these changes are obvious, then all inventions and creations using the concepts of the invention are protected.

[0041] Embodiment 1: taking West Lake of Hanshiqiao Wetland, Shunyi District, Beijing as an example to determine dynamic wetland boundary;

[0042] As shown in FIG. 1, the method for determining dynamic wetland boundary based on hydrology, organism and soil elements in the invention is carried out according to the following steps:

Step 1: Extracting Dynamic Wetland Hydrology Boundary

[0043] As shown in FIG. 2, the dynamic wetland hydrology boundary in step 1 includes two boundaries: actual water surface boundary, and animal and plant remain boundary. The actual water surface boundary is obtained through on-site investigation by high-precision handheld GPS or through remote sensing image with a resolution of less than 2 m; the animal and plant remain boundary is obtained by searching for the animal and plant remains through on-site investigation and then tracking and recording by high-precision handheld GPS, wherein the animal and plant remain boundary is the farthest hydrology boundary. The dynamic wetland hydrology boundary is obtained by overlaying the actual water surface boundary and the animal and plant remain boundary;

Step 2: Obtaining Data on Wetland Vegetation

[0044] Investigation line transects are laid in the direction perpendicular to the wetland hydrology boundary, with an interval of 100 m-500 m among the line transects, and the water surface boundary is taken as the starting point of investigation line transect; three investigation quadrats (1 m×1 m) are set on the line transects every 10 m, the vegetation coverage, the plant species and the numbers and coverage of various plants in the quadrats are recorded, and the species of hygrophytes and xerophytes thereof are determined; the importance value of each species is calculated, wherein the importance value is the sum of the relative coverage, relative density and relative frequency of each species; the importance value proportion of hygrophyte is calculated;

Step 3: Obtaining Data on Wetland Soil

[0045] In the investigation quadrats, the soil samples within the range of 0 cm-10 cm of surface layer are collected, weighed and then dried at 105° C., and finally the soil moisture content is calculated by weight difference;

Step 4: Simulating Dynamic Wetland Vegetation Boundary and Dynamic Wetland Soil Boundary

[0046] Taking the animal and plant remain boundary, i.e., the farthest hydrology boundary, as a benchmark, and calculating the minimum perpendicular distance of quadrat away from the actual water body boundary:

[00002]$D={\left[{\left(x_0-x\right)}^2+{\left(y_0-y\right)}^2\right]}^{\frac{1}{2}}{D}_m=\underset{\left(x,y\right)\in R}{\min }D$

[0047] Where, (x.sub.0, y.sub.0) and (x, y) represent the decimal latitude and longitude coordinates of pixels of each quadrat and hydrology boundary respectively, R represents the decimal latitude and longitude coordinates set of all pixels of hydrology boundary, D represents the minimum perpendicular distance of quadrat away from the actual water body boundary, D.sub.m represents the quadrat with the minimum perpendicular distance away from the actual water body boundary and with the greatest influence by water body, and the smaller the index D, the closer the hydrological connectivity between the quadrat and the wetland water body;

[0048] As shown in FIGS. 3 and 4, regression analysis is conducted on the minimum perpendicular distance of quadrat away from the actual water body boundary, the hygrophyte proportion and the measured soil moisture content, so as to obtain a regression equation;

[0049] Spatial statistics toolbox—spatial relationship modeling tool set—geographically weighted regression tool is selected in ArcGIS software, the wetland vegetation boundary or wetland soil boundary is simulated by geographically weighted regression model, curves of the measured and estimated hygrophyte proportion or soil moisture content are fitted, and the simulation effect of the model is verified. As shown in FIGS. 5 and 6, the fitted curve R.sup.2 is greater than 0.6, this shows that the simulation result of wetland boundary is reliable; as shown in FIGS. 7 and 8, when the soil moisture content is 34%-37% or the hygrophyte proportion is 45%-55%, the soil boundary or vegetation boundary is obtained respectively;

Step 5: Determining Dynamic Wetland Boundary

[0050] As shown in FIG. 9, considering the wetland hydrology boundary, wetland vegetation boundary and wetland soil boundary comprehensively, and taking the actual water surface boundary as a benchmark, so as to analyze the differences of the three boundaries on the perpendicular line of each pixel of water surface boundary; taking the boundary point with the largest D among the three boundary points as the wetland boundary of this perpendicular line, so as to obtain the dynamic wetland boundary, as shown in FIG. 10.