CN112782385B - Method for estimating density and reserve of soil organic carbon - Google Patents
Method for estimating density and reserve of soil organic carbon Download PDFInfo
- Publication number
- CN112782385B CN112782385B CN202011471133.9A CN202011471133A CN112782385B CN 112782385 B CN112782385 B CN 112782385B CN 202011471133 A CN202011471133 A CN 202011471133A CN 112782385 B CN112782385 B CN 112782385B
- Authority
- CN
- China
- Prior art keywords
- soil
- organic carbon
- density
- carbon content
- reserve
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention discloses a method for estimating density and reserve of soil organic carbon, belonging to the technical field of soil measurement, and comprising the following steps: setting soil sampling point positions, collecting profile samples of soil bodies, and measuring the thickness of each soil layer; testing the organic carbon content of each soil layer, and the volume percentage and volume weight index of the gravel with the thickness of more than 2 mm; calculating the organic carbon density of soil in unit centimeter depth of each soil layer sample, and establishing a linear regression model of the organic carbon density of the soil in unit centimeter depth and the organic carbon content of the soil; and respectively establishing the functional relations between the organic carbon density and the organic carbon content of the soil and between the organic carbon reserve and the organic carbon content based on the established linear regression model. According to the method, the carbon density and carbon reserve of the soil are directly, quickly and accurately estimated according to the organic carbon content, and the acquirability of the organic carbon reserve data of the soil is improved.
Description
Technical Field
The invention relates to the technical field of measurement soil science, in particular to a method for estimating organic carbon density and reserve of soil.
Background
Soil organic carbon is an important component of soil, is a source of plant nutrients, and plays an important role in forming soil aggregates, increasing soil water content, and reducing soil degradation and erosion. The soil is the third large carbon reservoir in the earth system, and the total organic carbon content of the global soil exceeds the total organic carbon content of the atmospheric carbon reservoir and the ecological carbon reservoir. The small change of the organic carbon in the soil can generate large influence on the carbon dioxide in the atmosphere, so that the climate change is influenced, and the rapid and accurate estimation of the soil organic carbon reserve is a precondition and a basis for realizing the analysis of the space-time change thereof and further evaluating the climate change.
Soil carbon density is used for evaluating and measuring the organic matter in soilAn extremely important indicator of carbon reserves. The soil organic carbon density refers to the reserve of soil organic carbon in a soil layer with a certain depth in unit area. If a section of a certain soil body consists of m layers, the organic carbon density SOCD (t/hm) of the section 2 ) The calculation formula of (c) is:
SOC in the formula i The organic carbon content (g/kg) of the soil of the i-th layer, BD i The i-th layer is soil volume weight (g/cm) 3 )，H i Thickness of soil layer of i-th layer, G i Is the volume percentage (%) of the ith layer of gravel larger than 2 mm. The calculation formula of the soil organic carbon density is not complicated, but the volume weight and the gravel content data are difficult to obtain (sometimes even cannot be measured), so that the calculation of the soil carbon storage is very inconvenient.
Disclosure of Invention
The invention aims to provide a soil organic carbon density and reserve estimation method to solve the problems that volume weight and gravel content data are difficult to obtain and over-dependent in the current soil organic carbon density and carbon reserve calculation process.
The technical scheme is as follows: in order to realize the purpose of the invention, the technical scheme adopted by the invention is that the method for estimating the density and the reserve of the soil organic carbon comprises the following steps:
step A: setting soil sampling point positions, collecting profile samples of soil bodies, and measuring the thickness of each soil layer; testing the organic carbon content of each soil layer, and the volume percentage and volume weight index of the gravel larger than 2 mm;
and B: calculating the organic carbon density of soil in unit centimeter depth of each soil layer sample, and establishing a linear regression model of the organic carbon density of the soil in unit centimeter depth and the organic carbon content of the soil;
and C: and B, respectively establishing the functional relations between the organic carbon density and the organic carbon content of the soil and between the organic carbon reserve and the organic carbon content based on the linear regression model established in the step B.
Further, the density of the organic carbon in the soil with the unit centimeter depth in the step B is the reserve of the organic carbon in the soil with the unit area of 1cm thickness, and the expression is as follows: SOCD pcm SOC × BD × H × (1-G/100)/10, where SOC is soil organic carbon content, BD is soil bulk density, H is soil thickness, and G is a volume percentage of gravel larger than 2 mm; establishing a linear regression model of the soil organic carbon density and the soil organic carbon content in unit centimeter depth, wherein the expression is as follows:
SOCD pcm and = a × SOC + b, where a and b are both calibration coefficients.
Further, the modeling precision of the linear regression model in the step B is determined by a coefficient R 2 And root mean square error RMSE, where R 2 The larger the RMSE, the smaller the RMSE, the higher the modeling accuracy.
Further, the organic carbon density SOCD of the soil layer i in the step C i Expressed as a function of the organic carbon content of the soil, the expression is: SOCD i ＝H i ×(a×SOC i + b) wherein H i The thickness of the ith layer; the section of a certain soil body consists of m soil layers, the organic carbon density SOCD of the section is expressed as a function relation of the organic carbon content of the soil, and the expression is as follows:
the soil organic carbon storage SOCS is expressed as a function of soil organic carbon content, and the expression is as follows:
wherein S is the area of the study region.
Has the advantages that: compared with the prior art, the technical scheme of the invention has the following beneficial technical effects:
the method for estimating the soil organic carbon density and the soil reserve provided by the invention avoids the problem of excessive dependence on soil volume weight and gravel content data, and shortens the detection time. The method provided by the invention has the advantages of high precision, good stability, certain universality and suitability for popularization.
Drawings
FIG. 1 is a method flow diagram of the present invention as a whole;
FIG. 2 is a plot of organic carbon density versus organic carbon content per centimeter of soil depth for example 1 of the present invention (n = 692);
fig. 3 is a scatter plot of organic carbon density versus organic carbon content per centimeter of depth of soil in example 2 of the present invention (n = 110).
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of the present invention and the above-described drawings, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
In this embodiment, a method for estimating soil organic carbon density and reserve is provided, and fig. 1 is a flow chart of a method for estimating soil organic carbon density and reserve according to an embodiment of the present invention, where the flow chart includes the following steps:
step A: setting soil sampling point positions aiming at a research area, collecting profile samples of soil bodies in the field, and measuring the thickness of each soil layer; the laboratory tests the organic carbon content of each soil layer, and the volume percentage and volume weight index of the gravel larger than 2 mm. The soil sample collection and index acquisition are implemented strictly according to the requirements of the manual of field soil description and sampling and the analysis method of the soil survey laboratory.
And B: soil organic carbon density per centimeter depth (SOCD) pcm ，t/hm 2 ) The organic carbon reserve of the soil with the unit area of 1cm thickness is expressed as a formula I: SOCD pcm (1-G/100)/10, wherein SOC is the organic carbon content (G/kg) of the soil, and BD is the volume weight (G/cm) of the soil 3 ) H is the volume percentage (%) of the gravel with the soil layer thickness of 1cm and G is more than 2 mm. Calculating the organic carbon density of soil in unit centimeter depth of each soil layer sample, and establishing a linear regression model of the organic carbon density of the soil in unit centimeter depth and the organic carbon content of the soil, wherein the expression is as follows:
SOCD pcm and = a × SOC + b, where a and b are both calibration coefficients.
Step C: and B, respectively establishing a functional relation between the Soil Organic Carbon Density (SOCD) and the organic carbon content and a functional relation between the organic carbon storage capacity (SOCS) and the organic carbon content on the basis of the linear regression model established in the step B. Organic carbon Density of soil layer i SOCD i Expressed as a function relation of soil organic carbon content SOC, the expression is three: SOCD i ＝H i ×(a×SOC i + b) wherein H i Is the thickness (cm) of the soil layer of the ith layer. The section of a certain soil body consists of m soil layers, the organic carbon density SOCD of the section is expressed as a function relation of the organic carbon content of soil, and the expression is as follows:
the soil organic carbon storage SOCS is expressed as a function relation of soil organic carbon content, and the expression is as follows:
wherein S is the area of the study region.
Example 1: (Chinese soil system classification method is adopted for dividing soil layers in soil section)
This example takes the estimation of organic carbon density and carbon reserve in soil in Jiangsu province as an example. The soil systems in the range of Jiangsu province are mainly distributed in 5 soil classes of artificial soil, eluvial soil, salinized soil, prototype soil and new soil, 173 representative soil system sections and 692 soil layer samples are obtained. The organic carbon content ranged between 0.17g/kg and 28.59g/kg with an average of 6.51g/kg. The organic carbon density (SOCD) of the soil in unit centimeter depth is obtained by calculation pcm ) Between 0.029t/hm 2 And 4.295t/hm 2 In between, the average value is 0.867t/hm 2 。
Construction of soil organic carbon Density per centimeter depth (SOCD) pcm ) And a linear regression model of the organic carbon content to obtain an expression as follows: SOCD pcm =0.1119 × SOC +0.1387. Fig. 2 is a scatter diagram of the organic carbon density and organic carbon content of soil per centimeter depth in example 1 of the present invention, and it can be seen how well the organic carbon density and organic carbon content of soil samples per centimeter depth fit. Modeling determination coefficient R between organic carbon density and organic carbon content of soil in unit centimeter depth 2 The modeling root mean square error RMSE is only 0.16t/hm as high as 0.9374 2 And the model has the capability of high-precision prediction. The result shows that the method can accurately and quickly estimate the organic carbon density of the soil in the depth of unit centimeter through the organic carbon content.
Further, the organic carbon density SOCD of the soil layer i can be obtained i The expression is: SOCD i ＝H i ×(0.1119×SOC i + 0.1387) wherein H i The thickness (cm) of the soil layer of the ith layer; the section organic carbon density SOCD expression of a certain soil body is as follows:
the expression of the organic carbon reserve of the soil is as follows:
wherein S is the area of the study region.
Example 2: (soil layer in soil profile is divided by fixed sampling depth)
This example takes the estimation of soil organic carbon density and carbon storage in east area Fan Gongdi, east tai city, jiangsu province, as an example. The research area is saline soil developed from coastal sediments for nearly millennium, 11 sections are sequentially collected from east to west in Fan Gongdi in an east area and a coastline, 10 samples are collected from each section according to a specific sampling interval, the sampling depth is 0-5cm, 5-10cm, 10-15cm, 15-20cm, 20-30cm, 30-40cm, 40-50cm, 50-60cm, 60-80cm and 80-100cm, and 110 soil samples are collected in total.
The organic carbon content in the study area ranged between 1.00g/kg and 17.03g/kg, with an average of 4.67g/kg. Calculating to obtain the Soil Organic Carbon Density (SOCD) in unit centimeter depth by an expression I pcm ) In the range of 0.13t/hm 2 And 1.86t/hm 2 In between, the average value is 0.64t/hm 2 . Soil organic carbon density per centimeter depth (SOCD) pcm ) Both are significantly correlated at the 0.01 level (two-sided), with Pearson correlation coefficients as high as 0.985. Constructing a linear regression model of the two, and obtaining an expression as follows: SOCD pcm =0.1243 × SOC +0.0605. FIG. 3 is a scatter plot of organic carbon density and organic carbon content in soil per centimeter depth in example 2 of the present invention, which shows that the fitting degree of organic carbon density and organic carbon content in soil sample per centimeter depth is very high, and that the model determination coefficient R is 2 The modeling root mean square error RMSE is only 0.08t/hm as high as 0.9699 2 . The result shows that the method can accurately and quickly estimate the organic carbon density of the soil in the depth of unit centimeter through the organic carbon content.
Further, the organic carbon density SOCD of the soil layer i can be obtained i The expression is as follows: SOCD i ＝H i ×(0.1243×SOC i + 0.0605); the section organic carbon density SOCD expression of a certain soil body is as follows:
the expression of the organic carbon reserve of the soil is as follows:
in the formula H i The thickness (cm) of the soil layer of the ith layer and S is the area of the research area.
Although in the above examples two types of soil profile division data were used for soil organic carbon density and carbon reserve estimation, the method is still applicable to other soil profile hierarchical division data, for example, the method is still applicable to soil organic carbon density and carbon reserve estimation based on the thickness of the occurrence layer.
Aiming at the problems that the volume weight and the gravel content data are difficult to obtain in the soil organic carbon density calculation, and further the soil carbon reserve calculation is extremely inconvenient, the method establishes a new method and a new model for accurately and quickly estimating the organic carbon density and the organic carbon reserve through the organic carbon content through the high fitting degree of the soil organic carbon density and the organic carbon content in the unit centimeter depth; the method solves the problems in the prior art, and can obtain the organic carbon density and carbon reserve of the soil under the condition of volume weight and gravel content data loss.
Claims (1)
1. A soil organic carbon density and reserve estimation method is characterized by comprising the following steps: the method comprises the following steps:
step A: setting soil sampling point positions, collecting profile samples of soil bodies, and measuring the thickness of each soil layer; testing the organic carbon content of each soil layer, and the volume percentage and volume weight index of the gravel with the thickness of more than 2 mm;
and B: calculating the organic carbon density of soil in unit centimeter depth of each soil layer sample, and establishing a linear regression model of the organic carbon density of the soil in unit centimeter depth and the organic carbon content of the soil;
and C: respectively establishing the functional relations between the organic carbon density and the organic carbon content of the soil and between the organic carbon reserve and the organic carbon content based on the linear regression model established in the step B;
in the step B, the density of the organic carbon in the soil with the depth of unit centimeter refers to the reserve of the organic carbon in the soil with the thickness of 1cm in unit area, and the calculation expression is as follows:
SOCD pcm ＝SOC×BD×H×(1-G/100)/10，
wherein SOC is the organic carbon content of soil, BD is the soil volume weight, H is the soil layer thickness, and G is the volume percentage of gravel larger than 2 mm; establishing a linear regression model, wherein the expression is as follows:
SOCD pcm = a × SOC + b, where a and b are both calibration coefficients;
organic carbon Density of soil layer i SOCD in step C i Expressed as a function of the organic carbon content of the soil, the expression is:
SOCD i ＝H i ×(a×SOC i + b) wherein H i The thickness of the ith layer;
the section of a certain soil body consists of m soil layers, the organic carbon density SOCD of the section is expressed as a function relation of the organic carbon content of soil, and the expression is as follows:
the soil organic carbon storage SOCS is expressed as a function of soil organic carbon content, and the expression is as follows:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011471133.9A CN112782385B (en) | 2020-12-14 | 2020-12-14 | Method for estimating density and reserve of soil organic carbon |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011471133.9A CN112782385B (en) | 2020-12-14 | 2020-12-14 | Method for estimating density and reserve of soil organic carbon |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112782385A CN112782385A (en) | 2021-05-11 |
| CN112782385B true CN112782385B (en) | 2022-12-02 |
Family
ID=75750932
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011471133.9A Active CN112782385B (en) | 2020-12-14 | 2020-12-14 | Method for estimating density and reserve of soil organic carbon |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112782385B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110174506B (en) * | 2019-06-28 | 2020-07-07 | 中国科学院地球化学研究所 | Method for estimating organic carbon in soil in karst region |
| CN114663782A (en) * | 2022-04-13 | 2022-06-24 | 中科海慧（天津）科技有限公司 | Forestry carbon source remittance method for forestry environment data |
| CN114858977A (en) * | 2022-04-29 | 2022-08-05 | 中国科学院东北地理与农业生态研究所 | Method and system for calculating field organic carbon conversion rate of straws |
| CN115480046A (en) * | 2022-10-11 | 2022-12-16 | 中国科学院地理科学与资源研究所 | Estimation method for predicting soil carbon reserves based on loess thickness |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105699624A (en) * | 2016-03-07 | 2016-06-22 | 中国科学院南京土壤研究所 | Soil organic carbon storage amount estimation method based on soil genetic horizon thickness prediction |
| CN110174506A (en) * | 2019-06-28 | 2019-08-27 | 中国科学院地球化学研究所 | A kind of Karst region soil organic matter evaluation method |
-
2020
- 2020-12-14 CN CN202011471133.9A patent/CN112782385B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105699624A (en) * | 2016-03-07 | 2016-06-22 | 中国科学院南京土壤研究所 | Soil organic carbon storage amount estimation method based on soil genetic horizon thickness prediction |
| CN110174506A (en) * | 2019-06-28 | 2019-08-27 | 中国科学院地球化学研究所 | A kind of Karst region soil organic matter evaluation method |
Non-Patent Citations (4)
| Title |
|---|
| 区域土壤有机碳密度及碳储量计算方法探讨;徐艳等;《土壤通报》;20051231;第36卷(第6期);第836-839页 * |
| 华南地区几种典型人工林土壤有机碳密度及其与土壤物理性质的关系;冯嘉仪等;《华南农业大学学报》;20180131;第39卷(第1期);第85页1.2样地设置与土壤采集分析和1.3数据处理，第86-87页2.3土壤有机碳密度与有机碳含量的关系 * |
| 我国东北地区土壤有机碳密度和储量的估算研究;孙维侠等;《土壤学报》;20040331;第41卷(第2期);第298-300页 * |
| 浅谈土壤有机碳密度及储量的估算方法;邵月红等;《土壤通报》;20061031;第37卷(第5期);第1007-1011页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112782385A (en) | 2021-05-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112782385B (en) | Method for estimating density and reserve of soil organic carbon | |
| Casson et al. | Simulation and extremal analysis of hurricane events | |
| Yakir et al. | Hydrologic response of a semi-arid watershed to spatial and temporal characteristics of convective rain cells | |
| CN109543356A (en) | Consider the ocean interior temperature-salinity structure remote sensing inversion method of Space atmosphere | |
| CN113221439B (en) | BP neural network-based drainage system real-time calibration and dynamic prediction method | |
| CN112035448A (en) | Neural network short-term rainfall forecasting method integrating water vapor and meteorological elements of foundation GNSS | |
| CN112329334B (en) | MWHTS and MWTS-II fusion inversion sea surface air pressure method based on simulated bright temperature | |
| CN115147239B (en) | Method for calculating environmental background water temperature and calculating temperature rise of newly-built coastal power plant temperature rise area | |
| CN114840616A (en) | Dynamic atmospheric natural environment modeling method based on space-time interpolation | |
| CN108663727B (en) | Method for estimating height of evaporation waveguide in world sea area by using evaporation rate | |
| CN115203625A (en) | Drought and waterlogging index data missing value interpolation method and device | |
| CN115100819A (en) | Landslide hazard early warning method and device based on big data analysis and electronic equipment | |
| CN109145258B (en) | Weibull distribution parameter confidence interval estimation method based on nonlinear fitting | |
| CN112329272A (en) | Residual life fitting prediction method integrating photovoltaic module degradation and life data | |
| CN112782096A (en) | Soil organic carbon density estimation method based on reflection spectrum data | |
| CN117172149A (en) | Evaporation waveguide prediction method based on data feature classification and neural network model | |
| CN110057990B (en) | PH correction method of multi-parameter water quality profiler | |
| CN115345075B (en) | Integrated aerosol pollution meteorological index-aerosol concentration estimation method and system | |
| Lee et al. | Joint probability distribution of significant wave height and peak wave period using gaussian copula method | |
| Connell et al. | An interpretation of radiosonde errors in the atmospheric boundary layer | |
| CN109374686A (en) | A kind of gas sensor | |
| CN111897032B (en) | Unknown precision correction method for tipping bucket type water quantity sensor based on unit time | |
| CN115236772A (en) | Data quality control system and method for drifting observation instrument | |
| CN113916181A (en) | Data processing method of surface-internal integrated deformation monitoring device | |
| WO2022217332A1 (en) | Monitoring an ambient air parameter using a trained model |
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 |