CN102818885A - Method for measuring salt storage in soil - Google Patents

Method for measuring salt storage in soil Download PDF

Info

Publication number
CN102818885A
CN102818885A CN2012103390518A CN201210339051A CN102818885A CN 102818885 A CN102818885 A CN 102818885A CN 2012103390518 A CN2012103390518 A CN 2012103390518A CN 201210339051 A CN201210339051 A CN 201210339051A CN 102818885 A CN102818885 A CN 102818885A
Authority
CN
China
Prior art keywords
soil
salt
grid
reserves
total
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.)
Pending
Application number
CN2012103390518A
Other languages
Chinese (zh)
Inventor
盛建东
武红旗
张丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xinjiang Agricultural University
Original Assignee
Xinjiang Agricultural University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xinjiang Agricultural University filed Critical Xinjiang Agricultural University
Priority to CN2012103390518A priority Critical patent/CN102818885A/en
Publication of CN102818885A publication Critical patent/CN102818885A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

The invention discloses a method for measuring salt storage in soil. The method is characterized by comprising the following steps: (a) arranging sampling points in a research area; (b) performing soil sample collection at each sampling point in three layers according to the coordinates of the sampling points, and performing analytic determination of indexes, including total salt, electric conductance and pH, on the collected soil samples according to a conventional analysis method; (c) using ArcGIS to perform interpolation on the total salt of the three layers of the soil and a total soil layer through an ordinary Kriging interpolation method; and (d) converting a Kriging interpolation diagram of the salt in the soil of the three soil layers into a grip map, and then utilizing a grid calculator in the ArcGIS to calculate the salt storage in each layer of the soil, wherein the formula is as follows, in the formula, i refers to the ith grid, k refers to the total number of grids in the research area, cell refers to the area of the grids, Xi refers the salt content predicted result of the ith grid, d refers to the thickness of the soil layer, and rho i refers to the soil volume weight of the ith grid. According to the method disclosed by the invention, the result is high in precision and not affected by outliers.

Description

A kind of measuring method of soil salt reserves
Technical field
The present invention relates to soil salt computing technique field, specifically a kind of measuring method of soil salt reserves.
Background technology
The calculating of soil salt reserves is basic contents of the research of landlocked arid biogeographic zone water resource utilization and ecologic environment aspect; Be to be study area, the origin mechanism of agroecological environment of carrying out and sustainable utilization of water resource, salinization of soil, irrigated area water salt evolution process, groundwater development, economize on water and control the important foundation content of research of the aspect such as relation and water resources rational regulation and control of salt with basin or enclosed areas.This also is one of foundation about study for globe change problem simultaneously.
Owing to weather that whole world change caused day by day climate arid He Shui, soil, living resources excessively developed, reason such as utilization, cause some representative basin district deterioration of the ecological environments, problems such as soil salinization of soil, desertification of land and grassland degeneration are serious day by day.Numerous researchers have been carried out many-sided research around the ecologic environment of key area, oasis, basin.These researchs mainly concentrate on aspects such as basin yardstick water salt balance, soil utilization variation, ecological environmental protection, and technological less to the salination spatial variability rule of basin soil itself and Evolution Mechanism research.
The technology of the calculating of regional soil salt reserve still less; Methods such as evaporate to dryness mass method, conductance method and main zwitterion summation calculating are arranged traditionally; But because these methods all need great deal of experiment data; And the experimental implementation process is complicated, time-consuming, energy consumption is big and measure relatively cause such as strictness of environmental requirement, causes its data precision and uses range and also be restricted.
Seeing that all have defective in the technology of the calculating of above-mentioned existing soil salt reserves; The inventor is according to working experience for many years and the professional knowledge of enriching; Actively study and innovate measuring method, have the influence that interpolation speed is fast, computational solution precision is high, do not receive special value with the soil salt reserves of founding a kind of novelty.
Summary of the invention
In order to solve the problems referred to above that exist in the prior art, the invention provides a kind of measuring method of soil salt reserves, have the computational solution precision height, do not receive the advantages such as influence of special value.
In order to solve the problems of the technologies described above, the present invention has adopted following technical scheme:
A kind of measuring method of soil salt reserves comprises the steps:
A. arrange sampled point at study area, and make the sampled point arrangenent diagram of study area;
B. according to the coordinate of sampled point, carry out the pedotheque collection for three layers by 0-30cm, 30-60cm, 60-100cm, and according to conventional method of analysis the pedotheque of gathering is comprised that total salt, electricity are led with the index analysis of pH and measure at each sampled point;
C. with ArcGIS total salt use Ordinary Kriging Interpolation method of interpolation (Ordinary Kriging) of study area soil 0-30cm, 30-60cm, 60-100cm and 0-100cm soil layer is carried out interpolation;
D. the Krieger interpolation graphs with 0-30cm, 30-60cm and three soil layer soil salts of 60-100cm converts grid map to, utilizes the salinity reserves of every layer of soil of grid counter calculating among the ArcGIS then, and its formula is following:
V = ∑ i = 1 k cell i × X i × d × ρ i
Wherein, V is regional a certain soil layer salinity reserves, and i is an i grid, and k is total grid number of study area, and cell is the area of grid, X iThe salt content that is i grid predicts the outcome, and d is a soil thickness, ρ iBe the soil weight of i grid.
Further, wherein the equispaced of sampled point is 2km.
Further, wherein said grid size is 50m * 50m.
Further, with the salinity reserves addition of 0-30cm, 30-60cm and three soil layer soil of 60-100cm, obtain total soil salt reserves of study area.
Compared with prior art, beneficial effect of the present invention is:
The measuring method of soil salt reserves of the present invention is based on the salinity and the unit weight data of grid cell; The method that Using geographic information system spatial analysis and geo-statistic are learned is carried out the calculating of basin soil salt total reserves; As long as this method is carried out the mensuration of total salt of indoor sample and unit weight; Utilize the optimum partially Krieger interpolation method of nothing to calculate; Compare with other method and to have the influence that interpolation speed is fast, computational solution precision is high, do not receive special value, also embodied the advantage that the inventive method is easy to realize visual map simultaneously.
Description of drawings
The sampled point arrangenent diagram of Fig. 1 for generating in the specific embodiment of the inventive method;
Fig. 2, Fig. 3, Fig. 4 and Fig. 5 are respectively the total salt Ordinary Kriging Interpolation interpolation graphs of the different soil (0-30cm, 30-60cm, 60-100cm, 0-100cm) that generates in the specific embodiment of the inventive method;
Fig. 6 calculates salinity reserves density profile for the layering that generates in the specific embodiment of the inventive method;
Fig. 7 calculates salinity reserves density profile for the 0-100cm layer that generates in the specific embodiment of the inventive method.
Embodiment
Below in conjunction with accompanying drawing and specific embodiment the present invention is described in further detail, but not as to qualification of the present invention.
The sampled point arrangenent diagram of Fig. 1 for generating in the specific embodiment of the inventive method; Fig. 2, Fig. 3, Fig. 4 and Fig. 5 are respectively the total salt Ordinary Kriging Interpolation interpolation graphs of the different soil (0-30cm, 30-60cm, 60-100cm, 0-100cm) that generates in the specific embodiment of the inventive method; Fig. 6 calculates salinity reserves density profile for the layering that generates in the specific embodiment of the inventive method; Fig. 7 calculates salinity reserves density profile for the 0-100cm layer that generates in the specific embodiment of the inventive method.In conjunction with Fig. 1 to Fig. 7, a kind of measuring method of soil salt reserves comprises the steps:
A. arrange sampled point at study area, and make the sampled point arrangenent diagram of study area, wherein the equispaced of sampled point is 2km;
B. according to the coordinate of sampled point, carry out the pedotheque collection for three layers by 0-30cm, 30-60cm, 60-100cm, and according to conventional method of analysis the pedotheque of gathering is comprised that total salt, electricity are led with the index analysis of pH and measure at each sampled point;
C. with ArcGIS total salt use Ordinary Kriging Interpolation method of interpolation (Ordinary Kriging) of study area soil 0-30cm, 30-60cm, 60-100cm and 0-100cm soil layer is carried out interpolation;
D. the Krieger interpolation graphs with 0-30cm, 30-60cm and three soil layer soil salts of 60-100cm converts grid map to; Grid size is 50m * 50m; Utilize grid counter among the ArcGIS (Raster caculation) to calculate every layer of soil salt reserves then, its formula is following:
V = ∑ i = 1 k cell i × X i × d × ρ i
Wherein, V is regional a certain soil layer salinity reserves, and i is an i grid, and k is total grid number of study area, and cell is the area of grid, X iThe salt content that is i grid predicts the outcome, and d is a soil thickness, ρ iBe the soil weight of i grid.
Wherein, the interval of sampled point is more little, and the sampled point that is provided with in the same area so is just many more, and result of calculation is also accurate more.But also increase workload simultaneously, so according to the interval of the requirement of precision as a result being confirmed sampled point.General sampled point be set to about 2Km or littler, its result of calculation can satisfy the demand of all researchs.Equally, the setting of grid size also need be sought the point that is fit between computational accuracy and calculated amount, and general grid size is got about 50m * 50m and promptly met the demands.
Being the study area measuring method of using soil salt reserves of the present invention with the Yan Qi basin below calculates the salinity reserves of this study area, so that understand the present invention better.
1, field sampling and in office analysis
Utilize the preliminary interpretation atural object of remote sensing images and the soil in study area Yan Qi basin in 2000 to utilize characteristic, arrange sampled point and make the sampled point arrangenent diagram of study area that with equispaced 2km the arrangenent diagram of sampled point is as shown in Figure 1.According to the coordinate of sampled point, utilize the GPS technology to carry out the pedotheque collecting work at study area, dig section respectively at 474 sampled points, press 0-30cm, 30-60cm, 60-100cm three level sampling.The pedotheque of gathering being carried out index analysis such as total salt, electricity are led, pH according to conventional method of analysis then measures.
2, basin yardstick soil salt spatial distribution map is drawn
2.1 soil salt statistics characteristic analysis
The index of the sample that obtains according to above-mentioned steps is carried out statistical study to the study area soil salt, and The result of statistics sees the following form 1, and table 1 is the total salt statistic analysis result of basin soil over sixty years of age for here.
Table 1
From table 1, can find out, here over sixty years of age the basin soil salinity average along with the degree of depth reduces fast.Here the total salt Variation Lines number average of each layer of basin soil over sixty years of age shows strong Characteristics of spatial variability greater than 100%.
2.2 the Characteristics of spatial variability analysis of each layer of soil salt salinity
Obtain the variation function theoretical model of the total salt of basin soil over sixty years of age here and related parameter is arranged through Characteristics of spatial variability analysis, specifically see the following form 2 each layer of soil salt salinity.The variation function theoretical model of table 2 total salt of basin soil over sixty years of age and related parameter is arranged for here.
Table 2
Figure BDA00002134008500052
Semivariance calculating through the total salt of soil in the last table 2 shows that the total salt semi-variance function of Yan Qi basin each layer soil theoretical model is spherical model (Spherical model).Except the 60-100cm soil layer of Yan Qi basin; Whole each layer of study area soil salt piecemeal gold effect (Co/Sill) is all less than 25%; Show that the study area soil salt has strong spatial coherence; Its spatial variability mainly receives structural factor such as parent soil material (matrix saliferous), landform (basin), the influence of weather natural causes such as (arids), can not ignore artificial factor simultaneously.
From effective range of semivariance model, the soil salt range in Yan Qi basin has exceeded the scope of study area, shows that this regional soil salinity has strong spatial coherence in whole basin, and this maybe to belong to the basin landform relevant with this area.
2.3 the total salt Krieger interpolation analysis of each layer of soil
Use ArcGIS total salt use Ordinary Kriging Interpolation method of interpolation (Ordinary Kriging) of study area soil 0-30cm, 30-60cm, 60-100cm and 0-100cm soil layer is carried out interpolation; See Fig. 2-Fig. 5, Fig. 2, Fig. 3, Fig. 4 and Fig. 5 are respectively the total salt Ordinary Kriging Interpolation of the soil layer interpolation graphs of 0-30cm, 30-60cm, 60-100cm and the 0-100cm soil layer of generation.Soil salt divides grade standard and corresponding area sees the following form 3.Table 3 is the different salinization rank soil areas in Yan Qi basin and ratio (g/kg, km 2).
Table 3
Figure BDA00002134008500061
Can be known that by chart top layer, Yan Qi basin 0-30cm soil is maximum with solonchak area proportion, reaches 44.5%, non-salinization soil area accounts for 34.1%, and middle salinization soil and geavy saltization soil area quite and all are lower than 10%.In vertical direction; The Krieger interpolation graphs of 30-60cm and the total salt of 60-100cm soil layer is compared with study area 0-30cm soil layer; Along with the increase of soil depth, non-salinization soil area is and increases considerably, and the solonchak area then is significantly and reduces; Middle salinization soil and geavy salt soil then change not obvious, and salinity is by the trend development that heavily lightens.Simultaneously, the interpolation graphs of 0-100cm weighting salinity also shows same trend.
2.4 soil salt reserve estimate
Following table 4 for here over sixty years of age basin soil salinity calculation of reserves result.
Table 4
Figure BDA00002134008500071
Annotate: unit weight 1.45g/cm 3
The study area soil texture is main with sandy loam mainly, is generally 1.45g/cm according to the soil weight of data sandy loam 3The calculating of salinity reserves is carried out with these data in (Saxton and Rawls, 2006).Here over sixty years of age, basin soil salinity calculation of reserves result was: mean value method result of calculation is 5.10 * 10 7T (three layers of sum) and 5.23 * 10 7T (0-100cm).And utilize the result of Krieger interpolation calculation to be: 7.12x10 7T (three layers of sum), 7.65x10 7T (0-100cm).Relatively the result of calculation of two kinds of methods can know that the result of Krieger interpolation calculation is greater than mean value method.Mean value method calculates the salinity reserves with All Ranges as a computing unit, can not react the special heterogeneity of soil, especially for the actual reserves that can not reflect salinity in the strong zone of soil space variability (coefficient of variation is greater than 100%); The Ke Lige method of interpolation utilizes geo-statistic scientific principle opinion and method to carry out the extraction of soil salt spatial variability model and characteristic parameter; Obtain the salt content figure of different spatial after the interpolation; The zoning is divided into different grids; Calculate the salinity reserves of each grid respectively, add with method obtain the zone in soil nutrient reserves data, result of calculation has been reacted the actual conditions of soil salt spatial variability.
Above embodiment is merely exemplary embodiment of the present invention, is not used in restriction the present invention, and protection scope of the present invention is defined by the claims.Those skilled in the art can make various modifications or be equal to replacement the present invention in essence of the present invention and protection domain, this modification or be equal to replacement and also should be regarded as dropping in protection scope of the present invention.

Claims (4)

1. the measuring method of soil salt reserves is characterized in that, comprises the steps:
A. arrange sampled point at study area, and make the sampled point arrangenent diagram of study area;
B. according to the coordinate of sampled point, carry out the pedotheque collection for three layers by 0-30cm, 30-60cm, 60-100cm, and according to conventional method of analysis the pedotheque of gathering is comprised that total salt, electricity are led with the index analysis of pH and measure at each sampled point;
C. with ArcGIS total salt use Ordinary Kriging Interpolation method of interpolation of study area soil 0-30cm, 30-60cm, 60-100cm and 0-100cm soil layer is carried out interpolation;
D. the Krieger interpolation graphs with 0-30cm, 30-60cm and three soil layer soil salts of 60-100cm converts grid map to, utilizes the salinity reserves of every layer of soil of grid counter calculating among the ArcGIS then, and its formula is following:
V = ∑ i = 1 k cell i × X i × d × ρ i
Wherein, V is regional a certain soil layer salinity reserves, and i is an i grid, and k is total grid number of study area, and cell is the area of grid, X iThe salt content that is i grid predicts the outcome, and d is a soil thickness, ρ iBe the soil weight of i grid.
2. the measuring method of soil salt reserves according to claim 1 is characterized in that, wherein the equispaced of sampled point is 2km.
3. the measuring method of soil salt reserves according to claim 1 is characterized in that, wherein said grid size is generally 50m * 50m.
4. the measuring method of soil salt reserves according to claim 1 is characterized in that, with the salinity reserves addition of 0-30cm, 30-60cm and three soil layer soil of 60-100cm, obtains total soil salt reserves of study area.
CN2012103390518A 2012-09-13 2012-09-13 Method for measuring salt storage in soil Pending CN102818885A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2012103390518A CN102818885A (en) 2012-09-13 2012-09-13 Method for measuring salt storage in soil

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2012103390518A CN102818885A (en) 2012-09-13 2012-09-13 Method for measuring salt storage in soil

Publications (1)

Publication Number Publication Date
CN102818885A true CN102818885A (en) 2012-12-12

Family

ID=47303125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2012103390518A Pending CN102818885A (en) 2012-09-13 2012-09-13 Method for measuring salt storage in soil

Country Status (1)

Country Link
CN (1) CN102818885A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101696968A (en) * 2009-10-14 2010-04-21 四川农业大学 New method for monitoring heavy metal content in soil
CN101788553A (en) * 2010-03-02 2010-07-28 中国农业大学 Multiscale analysis method for vegetation indexes of refuse dump and soil nutrient space
KR20120068206A (en) * 2010-12-17 2012-06-27 한국지질자원연구원 Method for verifying representativeness of sample collected in contaminated soil

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101696968A (en) * 2009-10-14 2010-04-21 四川农业大学 New method for monitoring heavy metal content in soil
CN101788553A (en) * 2010-03-02 2010-07-28 中国农业大学 Multiscale analysis method for vegetation indexes of refuse dump and soil nutrient space
KR20120068206A (en) * 2010-12-17 2012-06-27 한국지질자원연구원 Method for verifying representativeness of sample collected in contaminated soil

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
谷海斌: "基于地统计学和 GIS 的灌区盐渍化土壤制图及评价", 《干旱区研究》 *
谷海斌: "灌区尺度土壤盐渍化调查与评价－以石河子灌区和玛纳斯灌区为例", 《新疆农业大学学报》 *

Similar Documents

Publication Publication Date Title
Winder et al. Temporal organization of phytoplankton communities linked to physical forcing
Hutley et al. Monsoonal influences on evapotranspiration of savanna vegetation of northern Australia
Gao et al. Dynamics of alpine grassland NPP and its response to climate change in Northern Tibet
CN106204312B (en) Rapid economic evaluation method for ion adsorption type rare earth ore
CN102508961A (en) Design method for high-resolution fully distributed hydrological model TOPX
CN107886216A (en) A kind of land carrying capacity analysis method based on Evaluation of Land Use Suitability
CN102034003A (en) Watershed hydrological model design method based on storage capacity curve and TOPMODEL
CN102915387B (en) A kind of power grid ice region distribution diagram method for drafting
CN108573302A (en) A kind of simulation of basin non-point source pollution loading and Best Management Practices optimization method
CN103413293A (en) Method for obtaining urban land use landscape pattern intrinsic observation scale
CN103081624B (en) Partitioned management method of soil nutrients
CN108446782A (en) A kind of optimization of land resource configuration method of consideration of multi-scale
CN107239615A (en) A kind of power station water purification footprint evaluates computational methods
Li et al. Study on spatio-temporal evolution of habitat quality based on land-use change in Chongming Dongtan, China
CN1651859A (en) Layered regression remote sensing sampling technique for estimating crops per unit yield and forest unit storage quantity
CN104698508A (en) Portable device and method for soil erosion field investigation
CN107705002A (en) The determination method of mining soil content of beary metal sampled point exception high level coverage
Wang et al. Remote sensing evaluation of Chinese mainland's comprehensive natural resources carrying capacity and its spatial-temporal variation characteristics
CN112347657A (en) Wind energy theoretical reserve assessment method
CN102818885A (en) Method for measuring salt storage in soil
Ren et al. Investigation into spatiotemporal characteristics of coastal winds around the Taiwan Island
Zhao et al. Evaluation of water resources carrying capacity in Shandong Province based on Fuzzy Comprehensive Evaluation
Tao et al. Drought monitoring and analysing on typical karst ecological fragile area based on GIS
CN113435631B (en) Flood forecasting method, flood forecasting system, readable storage medium and computing device
CN114897364B (en) Method for calculating liquid mineral resource quantity

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20121212