一种快速识别土壤有机碳关键输出区的多参数评价方法与流程

技术特征：

1.一种快速识别土壤有机碳关键输出区的多参数评价方法，其特征在于：包括以下步骤：

S1：计算多参数评估模型所需的四个关键参数；

(a1)流域土壤有机碳输出参数(organic carbon output index,COindex)：选定试验区，利用国家公共数据平台获取覆盖试验区的土壤类型和土地利用/覆被类型空间分布矢量数据图层，结合野外实测不同土壤类型下的(0-35cm)土壤有机碳垂直迁移速率，不同土地利用/覆被类型有机碳坡面水平迁移速率，确定试验区间内、单位面积上流域内每种土壤类型、每种土地利用/覆被类型经过自然或人为干扰进而从横向和垂向两个方面综合释放有机碳的潜力值，借助ArcGIS平台将有机碳潜力值作为属性信息赋予每种土壤类型和每种土地利用/覆被类型，进而生成覆盖试验区的有机碳潜力空间分布栅格图层，该有机碳潜力空间分布栅格图层即为流域土壤有机碳输出参数；

(a2)流域土壤有机碳输出迁移路径参数(organic carbon flowpath index,CFindex)：利用覆盖试验区的数字高程模型(DEM)，借助ArcGIS平台中开放式的建模工具(Start Model Builder)获取流域土壤有机碳的汇水路径，汇水路径即流域内每个单元格在地形因子影响下，在遵循相对高程由高到低、由陆地向下游水体进行迁移的原则下，所经过的最短路径，所有最短路径的集合构成的图层即为流域土壤有机碳输出迁移路径参数；

(a3)流域土壤有机碳输出径流指数(organic carbon runoff index,CRindex)：利用覆盖试验区的土壤类型、土地利用/覆被类型、地形、土壤质地、降雨格局的空间分布图所提供的属性信息，再结合径流曲线方程(Soil Conservation Service,SCS)估算流域内每个栅格点的坡面径流值参数(CN)，即促使土壤有机碳沿陆地坡面向水体发生迁移行为的基本动力，所有坡面径流值的集合构成的图层为流域土壤有机碳输出径流指数(organic carbon runoff index,CRindex)；

(a4)流域土壤有机碳输出降水强度指数栅格图层(Precipitation index，Pindex)：利用覆盖试验区及其外围的水文气象站获取研究时段的逐日降水量数据,结合每个站点的空间地理坐标信息，并借助利用ArcGIS平台空间分析模块，将各站点逐日降水统计数据转变为的覆盖试验区的月度和年度降水矢量点图层，再借助ArcGIS平台空间分析模块中的kriging插值方法生成覆盖试验区的月度和年度降水量栅格图层，所有降水量的的集合构成的图层为流域土壤有机碳输出降水强度指数栅格图层；

S2：构建和应用多参数评估模型；

(b1)构建试验区标准化评估参数矩阵R_ij(见公式1)：

式中χ_ij为第i个景观栅格单元中的第j个评估参数，i为覆盖试验区的所有景观单位栅格点，j为(a1)-(a4)步骤所估算的四个参数；

(b2)为(a1)-(a4)步骤所估算的四个参数分别建立权重系数，采用可消除专家经验主观性的主成分分析法确定各参数所对应的权重系数(ω_j)；

(b3)构建试验区标准化后的加权决策矩阵W_ij，即为标准化后的评估参数矩阵与其权重系数的乘积(R_ij×ω_j)；

(b4)根据构建的试验区标准化后的加权决策矩阵W_ij，确定试验区土壤有机碳输出最大可能值(H_j)和最小可能值(D_j)，上述两个参数分别依据公式(2)与(3)确定；

H_j＝maxw_ij (2)

D_j＝minw_ij (3)

(b5)通过计算试验区每个景观栅格单元与有机碳输出最大可能值和最小可能值之间的距离，分别为D⁺_i和D^-_i，分别参照公式(4)和公式(5)，并据此判定每个景观栅格单元与最大输出可能值或最小输出可能值的接近程度；

式中：D⁺_i值越大，表明试验区第i个待评价景观栅格单元与最大输出可能值间的距离越远，其释放有机碳潜力越小，反之亦然，D^-_i值越大，表明试验区第i个待评价景观栅格单元与最小输出可能值间的距离越远，其释放有机碳潜力越大，反之亦然；

(b5)通过公式(6)计算试验区每个景观栅格单元与有机碳输出最大可能值和最小可能值之间的距离，评估试验区研究区每景观栅格单元中释放有机碳的潜在可能性，P_i值越大，表明第i个备选单元向附近河流水体中释放有机碳的潜力越大；

S3：有机碳风险区识别与预警；

(c1)识别有机碳释放风险区：采用滤波窗口(filter window)法识别有机碳释放风险区。选择32×32滤波窗口，滤波条件为大于试验区内P_i均值的5倍，满足滤波条件获取的区域为有机碳关键输出源区。

(c2)有机碳释放源区预警：根据试验区的气象预报数据，获取试验区未来某段时间内降水概率和土壤湿度数据，并据此获得S1中的(a3)和(a4)参数在未来某段时间内的预报值，结合S1、S2和S3中(c1)方法可预报试验区在未来某段时间内的有机碳释放源区。

2.根据权利要求1所述的快速识别土壤有机碳关键输出区的多参数评价方法，其特征在于：所述步骤S1中(a3)流域土壤有机碳输出径流指数的计算过程参照M1)和M2)：

M1)试验区土地利用数据通过对覆盖试验区的Landsat TM/ETM影像进行监督分类和结合野外调查验证工作获取，试验区内不同土壤类型及其对应的土壤质地，可通过中国科学院南京土壤所提供的中国土壤数据库和野外典型采样得到，见表1，

表1土壤渗透系数范围

M2)径流指数(CR_index)的求算可分为以下五种情形进行求算：

i.对于坡度≤5°的试验区，根据划分的土壤类别(参照表1)，结合土壤前期湿润度(AMC)查找不同土地利用对应的CN值表；其中，土壤前期湿润度可划分为以下三类：AMC I为土壤干旱，但未达到萎蔫点；AMCⅡ为土壤水分平均状况；AMCⅢ为土壤水分呈饱和状态；当试验区土壤处于AMCⅡ状态时，径流指数(CR_index)参照主要土地利用类型对应的CN值表(见表2至表4)；

表2栽培农业用地CN值表

表3其他农业用地CN值表

表4城镇用地CN值表

ii.对于坡度>5°的试验区，在土壤含水量为AMCⅡ状态下，每个栅格单位的径流CN值(CN_2s)参照公式(7)进行求算：

式中CN_2s为在土壤含水量(Ⅱ)状态θ坡度下的CN修正值，CN₃计算方法见公式(8)；在该种情形下，径流指数(CR_index)为通过公式(7)求得的CN_2s值；

iii.对于坡度≤5°，土壤含水量为(AMCⅢ)状态下的试验区，每个栅格单元的径流CN值(CN₃)参照(8)进行求算：

CN₃＝CN₂exp[0.00673(100-CN₂)] (8)

在该种情形下，径流指数(CR_index)为通过公式(8)求得的CN₃值；

iv.对于坡度≤5°的试验区，在土壤水分为(AMC I)状态下，每个栅格单元的径流CN值为0；

v.对于坡度>5°的试验区，在土壤水分为(AMC Ⅲ)状态下，每个栅格单位的径流CN值(CN_3s)参照公式(9)进行求算。

CN_3s＝CN_2sexp[0.00673(100-CN_2s)] (9)。

3.根据权利要求1所述的快速识别土壤有机碳关键输出区的多参数评价方法，其特征在于：所述步骤S1中(a4)部分中的kriging插值方法具体计算过程参照公式(10)：

式中S₀为气象站，Z为S₀处的降雨量，Z₁为距离S_φ最近的三个降水观测站的距离加权平均值；Z₂为距离S_ψ最近的周围八个点的高程均值，λ_φ和λ_ψ为S₀分别到S_φ和S_ψ的距离权重。