一种基于反射光谱数据的土壤有机碳密度估算方法

摘要

本发明公开了一种基于反射光谱数据的土壤有机碳密度估算方法，属于计量土壤学技术领域，包括步骤：土壤样品采集和预处理，获取土壤样品的物理化学属性，计算每个土样的单位厘米厚度土壤有机碳密度；测量预处理后的土壤样品反射率，得到反射光谱数据；反射光谱数据预处理，包括断点校正，以及对光谱曲线进行平滑处理；构建单位厘米厚度土壤有机碳密度的反射光谱估算模型；建立土壤有机碳密度与反射光谱数据的函数关系。本发明通过土壤反射光谱数据直接、快速、精确估算土壤碳密度，提高了土壤有机碳储量数据的可获取性。

说明书

技术领域

本发明涉及计量土壤学技术领域，尤其涉及一种基于反射光谱数据的土壤有机碳密度估算方法。

背景技术

土壤碳主要包括土壤有机碳和土壤无机碳，其中土壤无机碳较为稳定。土壤有机碳是动植物残体、微生物体、土壤微生物作用生成的腐殖质的合称，主要分布在1米深度的土体内。土壤有机碳对土壤的物理、化学性质有着重要影响，并直接影响着土壤质量。土壤碳库储量约是大气碳库储量与陆地生态圈碳库储量的总和，是最大的陆地生态系统碳库，也是全球第三大碳库。较小的土壤碳库变化会对大气二氧化碳产生较大影响，因此土壤碳储量数据对于减缓土壤中温室气体的排放、预测全球气候变化趋势及采取应对决策具有重要意义。

土壤碳密度是评价和衡量土壤中有机碳储量的一个极其重要的指标。土壤有机碳密度是指单位面积一定深度的土层中土壤有机碳的储量。土壤碳密度的获取以土壤有机碳含量、土壤容重、土层厚度、大于2mm的石砾体积百分比等属性数据为基础。目前对土壤有机碳的测定主要为化学方法，具有污染性、成本高、检测周期长等不足，而容重和石砾含量数据的获取较为困难(有时甚至无法测定)，这些因素共同制约了土壤有机碳密度数据的获取。因此，现有技术中急需一种无污染、低成本、快速准确的土壤有机碳密度的检测方法。

发明内容

发明目的：针对以上问题，本发明提出一种基于反射光谱数据的土壤有机碳密度估算方法，以解决目前有机碳密度和碳储量计算过程中有机碳含量测定周期长成本高、以及对容重和石砾含量数据过度依赖的问题。

技术方案：为实现本发明的目的，本发明所采用的技术方案是：一种基于反射光谱数据的土壤有机碳密度估算方法，包括如下步骤：

步骤A：土壤样品采集和预处理，获取土壤样品的物理化学属性，计算每个土样的单位厘米厚度土壤有机碳密度SOCD

步骤B：测量预处理后的土壤样品反射率，得到350-2500nm波段的反射光谱数据；

步骤C：反射光谱数据预处理，包括断点校正，以及对光谱曲线进行平滑处理；

步骤D：构建单位厘米厚度土壤有机碳密度SOCD

步骤E：基于步骤D，建立土壤有机碳密度与反射光谱数据的函数关系。

进一步地，步骤A获取的土壤样品的物理化学属性包括有机碳含量、大于2mm石砾所占的体积百分比、土壤容重、土层厚度；单位厘米厚度土壤有机碳密度是单位面积1cm厚度的土壤有机碳储量，表达式为式一：SOCD

进一步地，步骤B采集土壤光谱前首先去除暗电流，然后利用漫反射标准参照板获取相对反射率；光谱采集过程中，光源与垂直方向成一定角度，与土壤样品距离为a；探头垂直于土壤样品进行光谱采集，与土壤样品距离为b；每个土壤样品采集n条光谱曲线，算术平均后得到该土壤样品的350-2500nm波段反射光谱数据。

进一步地，根据计算得到的单位厘米厚度土壤有机碳密度SOCD

土层i的有机碳密度SOCD

设某一土体的剖面由m个土层组成，该剖面的有机碳密度SOCD表示为土壤反射光谱数据的函数关系，表达式为式四：

有益效果：与现有技术相比，本发明的技术方案具有以下有益的技术效果：

光谱分析技术具有分析速度快、成本低、无危险、无破坏、可同时估算多种成分等特点，为土壤研究提供了新的手段与方法。通过本发明提供的一种基于反射光谱数据的土壤有机碳密度估算方法，解决了土壤有机碳密度和碳储量计算过程中属性数据测定费用高周期长、以及对容重和石砾含量数据过度依赖的问题，提供了一种通过土壤反射光谱数据直接、快速、精确估算土壤碳密度的方法，提高了土壤有机碳储量数据的可获取性。本发明为了解土壤有机碳密度状况以及研究土壤碳库对全球气候变化的影响提供了快速获取基础数据的新手段，适合推广应用。

附图说明

图1是本发明整体的方法流程图；

图2是本发明实施例的单位厘米厚度土壤有机碳密度真实值与预测值散点图；

图3是本发明实施例的1m剖面土壤有机碳密度真实值与预测值散点图。

具体实施方式

为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。

在本实施例中提供了一种基于反射光谱数据的土壤有机碳密度估算方法，图1是估算方法流程图，该流程包括如下步骤：

步骤A：土壤样品采集和预处理，获取土壤样品的物理化学属性，包括有机碳含量、大于2mm石砾所占的体积百分比、土壤容重、土层厚度；计算每个土样的单位厘米厚度土壤有机碳密度SOCD

步骤B：测量预处理后的土壤样品反射率，得到350-2500nm波段的反射光谱数据。

本实施例中，土壤反射光谱数据测量在暗室内采用FieldSpec 3便携式光谱分析仪进行。采集土壤光谱前首先去除暗电流，然后利用40cm×40cm漫反射标准参照板获取相对反射率。光谱采集过程中，光源为一盏50W卤素灯，与垂直方向角度为15°，距离样品30cm；采用5度视场角探头，距土壤样品15cm处垂直进行光谱采集；每个土壤样品采集20条光谱曲线，算术平均后得到该土壤样品的350-2500nm波段反射光谱数据。

步骤C：反射光谱数据预处理；首先使用软件ViewSpecPro进行断点校正，消除仪器自身原因产生的断点，然后对光谱曲线进行Savitzky-Golay平滑处理。

步骤D：根据计算得到的单位厘米厚度土壤有机碳密度SOCD

步骤E：基于步骤D，建立土壤有机碳密度与反射光谱数据的函数关系。

土层i的有机碳密度SOCD

本实施例以江苏省东台市范公堤以东地区土壤有机碳密度反射光谱估算为例。在研究区内与海岸线垂直自东而西依次采集了11个剖面，每个剖面按照特定采样区间采集10个样品，采样深度为0-5cm、5-10cm、10-15cm、15-20cm、20-30cm、30-40cm、40-50cm、50-60cm、60-80cm和80-100cm，共采集土壤样品110个。研究区有机碳含量范围介于1.00g/kg和17.03g/kg之间，平均值为4.67g/kg。通过表达式一计算得到单位厘米厚度土壤有机碳密度SOCD

对预处理后的土壤样品采用美国ASD公司生产的FieldSpec 3便携式光谱分析仪分别测定110个土壤样品的反射率，得到350-2500nm波段反射光谱数据。需要说明的是，光谱采集过程中，盛于容器内的土样随自动旋转平台旋转的同时测反射率，每个样本得到20条光谱曲线，算术平均后得到该土样的反射光谱数据。ASD光谱仪由三个光谱仪组成，在1000nm和1800nm附近会出现由于仪器转换产生的断点，因此使用软件ViewSpecPro进行断点校正，然后对光谱曲线进行Savitzky-Golay平滑处理。

偏最小二乘回归集典型相关分析、主成分分析与多元线性回归分析于一体，可以较好地解决自变量多重共线性问题，允许在样本个数少于变量个数条件下进行回归建模。该方法广泛应用于反射光谱估算土壤属性的研究。使用Unscrambler软件，采用偏最小二乘回归分析方法构建单位厘米厚度土壤有机碳密度SOCD

图2为本发明实施例的单位厘米厚度土壤有机碳密度真实值与预测值散点图，可以发现采用偏最小二乘回归分析方法构建单位厘米厚度土壤有机碳密度SOCD

在采用偏最小二乘回归分析方法构建单位厘米厚度土壤有机碳密度SOCD

在上述实施例中虽然仅使用了偏最小二乘回归分析方法构建单位厘米厚度土壤有机碳密度SOCD

本发明方法针对土壤有机碳密度和碳储量计算过程中属性数据测定费用高周期长、以及对容重和石砾含量数据过度依赖的问题，通过单位厘米厚度土壤有机碳密度SOCD