艾比湖湿地自然保护区土壤盐分多光谱遥感反演模型
周晓红1,2,3, 张飞1,2,3,*, 张海威1,2,3, 张贤龙1,2,3, 袁婕1,2,3
1. 新疆大学资源与环境科学学院智慧城市与环境建模自治区普通高校重点实验室, 新疆 乌鲁木齐 830046
2. 新疆大学绿洲生态教育部重点实验室, 新疆 乌鲁木齐 830046
3. 中亚地理信息开发利用国家测绘地理信息局工程技术研究中心, 新疆 乌鲁木齐 830002
*通讯联系人 e-mail: zhangfei3s@163.com

作者简介: 周晓红, 女, 1994年生, 新疆大学资源与环境科学学院智慧城市与环境建模自治区普通高校重点实验室硕士研究生,e-mail: 18290622359@163.com

摘要

土壤盐分是衡量土壤质量的要素, 也是作物生长发育的基本条件。 因此, 迫切地需要一种可以快速了解土壤盐分含量(SSC)的方法。 针对艾比湖湿地自然保护区, 基于Landsat8 OLI多光谱遥感影像, 以该研究区36个土壤表层样品的盐分含量为数据源, 选择相关性较好的多光谱遥感指数分析研究区土壤盐分分布状况, 并将其分别与实测SSC构建线性、 对数、 二次函数模型, 进而优选精度最高的模型来反演该研究区SSC。 结果表明: (1)在多光谱遥感指数中, 与SSC相关性最高的是增强型植被指数(EVI), 其相关性范围为(-0.70~-0.67); 其次是传统型植被指数(TVI), 其范围为(-0.58~-0.46); 土壤盐分指数(SI)与SSC的相关性最低, 其范围为(-0.45~0.16), 其中SI3和SI4与SSC均没有相关性。 (2)将实测土壤盐分值所反演的分布图与EVI对比分析, 发现在西北、 正南方向的艾比湖湖边周围和东北方向盐池桥的SSC均较高, 其EVI的值较低, 说明通过该研究区实测土壤盐分值所反演的盐分分布图与EVI的空间分布结果较为一致, 表明EVI对该地区土壤盐分具有一定的敏感性, 能较好地反演SSC的空间分布; (3)分别将三种EVI与实测SSC建模分析比较, 发现SSC与增强型比值植被指数(ERVI)所构建的二次函数模型最好; 其验证集的决定系数( R2)为0.92, 均方根误差(RMSE)为2.48, 相对分析误差(RPD)为2.09, 模型精度较高、 稳定性较为可靠, 相比之下, 说明ERVI对该湿地自然保护区土壤盐分有更高的敏感性, 可以用来预测该区域SSC, 从而进行空间反演。 在TVI中加入Landsat8多光谱遥感影像的b6和b7波段, 得到EVI, 以此来反演SSC是可行的, 且比传统可见光和近红外波段所构建的植被指数反演效果更好。 因此该研究不仅可以为遥感反演提供理论参考, 而且对该地区SSC的定量估算和动态监测具有重要的意义, 也可作为其他区域SSC预测反演的备选方案。

关键词: 自然保护区; 增强型植被指数; 土壤盐分指数; 反演模型
中图分类号:TP79 文献标识码:A
A Study of Soil Salinity Inversion Based on Multispectral Remote Sensing Index in Ebinur Lake Wetland Nature Reserve
ZHOU Xiao-hong1,2,3, ZHANG Fei1,2,3,*, ZHANG Hai-wei1,2,3, ZHANG Xian-long1,2,3, YUAN Jie1,2,3
1. College of Resources &Environmental Science, Xinjiang University, Urumqi 830046, China;
2. Key Laboratory of Oasis Ecology, Xinjiang University, Urumqi 830046, China
3. Engineering Research Center of Central Asia Geoinformation Development and Utilization, National Administration of Surveying, Mapping and Geoinformation, Urumqi 830002, China
Abstract

Soil salinity is an important factor for measuring soil quality, and it is also a basic condition for the growth of crops. Therefore, it is urgent to find a method that can understand soil salt content quickly. This paper is based on the Landsat8 OLI multispectral remote sensing image for the Ebinur Lake Wetland Nature Reserve, and we use the salt content of 36 soil surface samples in the study area as the data source, and choose several multispectral remote sensing indices which have the superior correlation with soil salinity to analyze the soil salinity distribution in the study area. The linear, logarithmic and quadratic function models were constructed with the measured soil salinity, and optimum inversion model of soil salt content was selected. The result shows that: (1)Among these multispectral remote sensing indices, the enhanced vegetation indices show the closest correlation with soil salinity, and the correlation coefficient range is between -0.67 and -0.70. The second is the traditional vegetation indices, and the correlation coefficient range is between -0.46 and -0.58. The correlation of soil salt index is the farthest t,and its range is between 0.16 and -0.45, and there is no correlation between SI3, SI4 and soil salt content. (2)Comparing and analyzing the salt distribution map inverted by measured soil salinity values and the spatial distribution of Enhanced Vegetation indices, we found that the soil salt content around the Ebinur Lake of the northwest and south direction and tne Yan Chi Bridge in the northeast is higher, but the enhanced vegetation indices are lower. The result shows that the salt distribution map inverted by measured soil salinity values is consistent with the spatial distribution of Enhanced Vegetation indices. It indicates that the enhanced vegetation indices have a higher sensitivity to soil salinity,which can better reverse the spatial distribution of soil salinity in the study area. (3)From the comparison and analysis of those models, which builds the three enhanced vegetation indices and measured soil salt content respectively.We found that the enhanced ratio vegetation index is the best choice to construct the quadratic function model. The determination coefficient of its validation set ( R2) is 0.92, and the root mean square error (RMSE) is 2.48, and the relative analysis error (RPD) is 2.09. The data show that this model is more accurate and reliable. In summary, ERVI is more sensitive to soil salinity and predict the soil salinity content, while is more suitable for inversion of soil salinity in this study area. Therefore, the study indicates that it is feasible to invert the soil salinity by the enhanced vegetation index constructed by the b6 and b7 band of Landsat8 multipectral remote sensing imagery. And its inversion effect is better than that of traditional visible light band. Therefore, this study not only provides a theoretical reference for remote sensing inversion, but also has important implications for the quantitative estimation and dynamic monitoring of soil salinity for the study area. Otherwise, it can be used as an alternative offer for prediction of soil salt content in other regions.

Key words: Natural reserve; Enhanced vegetation index; Soil salinity index; Inversion model
引 言

土壤盐渍化是由于自然和人类不合理的灌溉造成的, 主要发生在地势低且气候干旱的地区, 会严重影响耕地质量、 致使土壤结块, 农业生产损失, 严重破坏地理生态环境, 干扰人们的正常生活[1]。 因此, 快速了解盐渍土的特征以及发展现状等对治理盐渍土尤为重要, 同时也是盐渍土被科学改良、 合理管理利用的前提条件。 实地采样是获取土壤盐分含量(soil salinity content, SSC)的传统方法, 该方法不仅实测点少、 典型性差, 且不能大范围、 实时动态地监测。 遥感技术具有空间性、 实时性、 多波段性等特点, 不仅能够全面、 快速地监测SSC, 可以为快速获取SSC提供了一种新的技术方法。

目前, 遥感技术的一个研究热点就是通过各遥感指数的构建来分析某区域盐渍土的信息, 王飞等[2]以盐渍区实测SSC与landsat7和landsat8建立各种遥感指数和原始波段来构建线性和多个非线性模型, 探讨其对新疆玛纳斯河绿洲、 阜北地区和渭-库绿洲土壤盐度的敏感性; 罗冲等[3]利用Landsat8计算植被指数(VI), 采用PLRS法对SSC与植被指数进行模型构建、 比较, 从而确定玛纳斯流域绿洲土壤盐渍化分布状况; Alexakisa等[4]针对克里特岛中南部, 希腊土壤盐渍化问题, 基于WV2和Landsat 8计算土壤盐分指数(soil salinity index, SI), 可知在该区域S5对SSC最为敏感, 在高电导率ECe区域可为监测地表表层土壤盐分提供一种新的方案。 TahaGorji[5]等针对土耳其图兹湖周边的土壤盐渍化问题, 基于Landsat8 OLI遥感影像提取出土壤盐分指数, 经过回归分析, 利用反演模型最好的SI反演出土耳其图兹湖周边土壤盐渍化的空间分布。 根据以上研究可知, 通过计算多光谱遥感影像中的b3, b4和b5波段得到的TVI, 虽然可以用来分析盐渍土性状, 但用来识别盐渍土特征的准确性还是比较低的[6]。 TVI仅仅是通过计算b3, b4和b5而得到的, 并没考虑到与土壤盐分和植被相关的短红外波段, 所以TVI并没有包括全部的土壤、 植被信息, 故造成模型精度不高, 预测能力差、 难以快速准确地监测到SSC。

本文在TVI计算中加入波长更长、 信息量更大的b6和b7波段, 来检验是否能够提高SSC建模反演的准确性, 目前类似研究还比较少。 在TVI中引入短红外波段得到EVI, 并以艾比湖湿地自然保护区为研究对象, 拟尝试用实测土壤盐分数据和OLI多光谱遥感影像, 计算各种植被指数和SI等十种多光谱遥感指数, 探讨其对SSC的敏感性, 选出最强敏感指数, 构建最优SSC反演模型, 为该地区准确获取、 监测SSC提供理论依据。

1 实验部分
1.1 研究区概况

艾比湖湿地自然保护区地处新疆维吾尔族自治区博尔塔拉蒙古自治州境内, 地理位置43° 38'— 45° 52'N, 79° 53'— 85° 02'E之间, 三面环山, 远离海洋, 位于乌苏东部和托里县, 气温变化强烈, 干旱少雨, 总面积为2 670.85 km2, 包括湖滨地带形成5 500 km2的湿地。 近年来, 由于艾比湖人口剧增、 水土的大规模开发, 致使湖面缩小、 耕地盐渍化、 草场退化严重、 土地沙漠化等问题突出[7]。 所以人们正在大力改善生态环境, 但并没有改变当地恶劣的生态环境, 研究区地下水位每年仍然在下降, 枯死的荒漠植被越来越多, 气候也越来越干燥, 致使该地区的土壤盐分也在不断地发生变化。

图1 研究区概况和采样点分布Fig.1 Location of research area and distribution of sampling points

1.2 土样采样与多光谱遥感影像

土壤样品采样时间为2016年5月中旬, 以3~10 km采样点间隔为原则并结合保护区的自然景观、 土壤质地以及道路通达等实际情况, 在各采样点30 m× 30 m样方内且具有显著特征的地表表层选取土壤, 在每个样方内按照“ 五点采样法” 采集五份土壤, 将其均匀混合在一起制成混合土样带回实验室内自然风干、 磨碎, 然后用0.5 mm孔径的筛子过滤, 按水土比5∶ 1的比例提取浸提液, 使用DDS-307A型电导率仪测定其浸出液盐度, 并转换为土壤含盐量(g· kg-1)。

选用的多光谱遥感影像是由传感器陆地成像仪(OLI)在2016年5月获取的Landsat 8影像, 获取时间与采样时间相对应。 利用Envi 5.1软件, 对该多光谱影像预处理, 包括几何校正、 辐射校正、 大气校正、 彩色合成, 最后裁剪获得研究区影像; 然后, 利用ArcGIS10.3软件, 将实测点的经纬度导入多光谱遥感影像, 得到每个实测点所对应的反射率值, 从而进行SSC的反演。

1.3 多光谱遥感指数

本研究选取VI和SI分别与SSC进行相关性分析。 TVI分别是归一化植被指数(normalized difference vegetation index, NDVI)、 比值植被指数(ratio vegetation index, RVI)、 差值植被指数(difference vegetation index, DVI); TVI就是由常用多光谱遥感波段b3, b4计算所得, 而EVI就是在b3, b4和b5基础上再加入短红外波段b6和b7计算所得, 分别是增强型归一化植被指数(enhanced normalized difference vegetation index, ENDVI)、 增强型比值植被指数(enhanced ratio vegetation index, ERVI)、 增强型差值植被指数(enhanced difference vegetation index, EDVI); 土壤盐分指数(soil salinity index, SI)是将多光谱遥感影像的b1, b2, b3, b4和b5波段相互参与计算, 重新组合成B1, B2, B3和B4波段, 以此得到四个SI。

表1 多光谱指数计算公式 Table 1 Multispectral indices formula
1.4 模型反演与检验

首先给土样随机分组, 得到建模集样本24个, 验证集样本12个。 然后利用建模集SSC和选定的多光谱遥感指数, 分别建立各种多光谱遥感指数与SSC的回归模型。 在本研究中运用SPSS 19软件完成一元线性、 二次函数、 对数函数、 多元线性回归模型, 并将其进行比较分析, 从而选出最佳模型。

表2 土壤盐分的描述性统计分析 Table 2 Descriptive statistics analysis of soil salinity

利用12个验证样本对上述构建的模型进行检验, 从两个方面评价反演模型的精度即模型的稳定性和预测能力, 稳定性用R2的大小来检验[8], 预测能力用RMSE和RPD综合评定。 R2越大, 说明模型精度越高; RMSE越小, 模型预测能力越强; 当RPD> 2时, 模型预测效果较好, 1.4< RPD< 2时, 模型预测效果一般, RPD< 1.4时, 模型无预测能力。 各检验表达式如下所示:

R2=i=1n(Xi-X̅)(Yi-Y̅)i=1n(Xi-X̅)2i=1n(Yi-Y̅)2(1)

RMSE=i=1n(Yoi-Ysi)2n(2)

RPD=stdevrmse(3)

其中, n为采样点个数; YioYis分别为土壤样本i的观测值和模拟值。

2 结果与讨论
2.1 不同区域土壤盐分含量

植被的分布主要受气候条件、 地形和土壤性状的影响, 其中SSC是影响植被分布的关键因素, 经过野外实地考察, 发现在西北方向的艾比湖湖边周围, 常年无植被生长, 可能是因为该区域盐渍土较多; 在正南方向的艾比湖湖周边, 其植被稀疏; 而在正北、 东南方向上的山地, 植被生长较为茂盛, 可能是因为在东南方向, 阿奇克苏河的支流使该区域水份充沛且山地上少有人生活, 因此植被生长较为茂盛。 根据野外实地考察和SSC值, 将该研究区实测点主要分为植被生长区、 无植被生长(稀疏)区、 裸地三种区域, 如表3

表3 土壤盐分含量分布 Table 3 The distribution of SSC

该研究区属于典型湿地自然保护区, 路广等[9]研究发现在湿地土壤盐渍化分布区, 植被覆盖变化与土壤盐分呈负相关性, 盐分较高的土壤会抑制植被的生长发育, 在土壤盐分较低的区域, 有较多的植被生长。 王飞等[10]研究发现在SSC较高区域的土壤表层在植被指数灰度图中较亮。 由表3可以看出, 在有植被生长的区域, 其SSC较低, 范围在0~3.3 g· kg-1之间, 其实测点区域的土壤表层在植被指数灰度图亮度中较暗; 在无植被生长(稀疏)区, 其SSC在3.7~13.6 g· kg-1之间, 由于该区域土壤盐分比有植被生长区高, 因此在该区域的土壤表层在植被指数灰度图中的亮度比有植被生长区亮度更亮; 在裸地区, 其SSC较高, 均大于20 g· kg-1, 该区域的土壤表层在植被指数灰度图中较暗。

2.2 多光谱指数与土壤盐分含量的相关性分析

通过Envi5.1中的Band manth计算出TVI, EVI和SI, 将其分别与SSC进行相关性分析, 从表4可以看出ERVI, ENDVI, EDVI与SSC的相关性最强, 其相关系数分别是-0.70, -0.70和-0.67, 在0.01水平上显著相关; 其次是NDVI, RVI和DVI, 其相关系数分别是-0.46, -0.46和-0.58, 在0.01水平上显著相关; SI2与SSC的相关系数为-0.45, 在0.01的水平上显著相关, SI1与SSC的相关系数为-0.40, 在0.05的水平上显著相关, 而SI3, SI4与SSC均没有相关性。

表4 多光谱指数和土壤盐分的相关性 Table 4 The correlation of Multispectral index and soil salinity
2.3 增强型植被指数特征分析

表3比较可知, ERVI, ENDVI, EDVI与SSC的相关性最强。 植物叶子内部结构会影响植被的反射光谱。 健康植物的反射光谱值在380~3 000 nm波段范围内较高, 在2 000~2 500 nm光谱范围之间, 光谱曲线存在着明显的变化, 有一个明显的峰值, 说明植被的光谱信息在2 000~2 500 nm之间具有一定的敏感性, 因此加入远红外波段的EVI, 可以用来表示植被的生长, 从而间接地了解到该地区的SSC。

基于与SSC相关性最强的EVI, 对SSC进行定量分析, 图2显示了该研究区实测点盐分反演的盐分分布图和EVI图。 从实测点盐分反演的盐分分布图目视可知, SSC较高的区域主要分布在西北、 正南方向的艾比湖湖边周围和东北方向的盐池桥, 这些区域的EVI都相对偏低。 主要是因为在艾比湖湖边周围土壤盐渍化较为严重, 因此植被生长较少, 植被覆盖率较低, 所以EVI值较低。 而在东北方向上的盐池桥, 虽然SSC较高, 但其EVI值没有很低, 可能是因为在该区域两边都生长着茂盛的植被, 该区域可能有适宜植被生长的其他条件, 所以该区域生长着少量的植被, 因此其EVI值比艾比湖周边的高; 在正北、 东南方向的山地, 其SSC较低, EVI值较高, 可能是因为在该区域SSC较少, 且少有人在此生活, 因此该区域植被生长较为茂盛, EVI值也相对较高。

图2 盐分分布图与增强型植被指数Fig.2 The distribution map of soil salinity and enhanced vegetation index

2.4 土壤盐分含量模型反演优选

在相关性分析和SSC定性分析的基础上, 将ERVI, ENDVI, EDVI与SSC分别进行回归分析, 构建简单的回归模型。 由表5可以看出, 首先分别以ERVI, ENDVI, EDVI为自变量建立多元线性回归模型, 可得到建模集R2为0.59, 验证集R2为0.75, RMSE为3.42, 其次对ENDVI, EDVI, ERVI分别建立一元线性、 对数函数、 二次函数回归模型比较分析, 可以得到以ENDVI为自变量所构建的最优模型为二次函数模型, 其建模集R2为0.51, 验证集R2和RMSE分别为0.61和3.88; 以EDVI为自变量所构建的最优模型为二次函数模型, 其建模集R2为0.53, 验证集R2和RMSE分别为0.74和3.02; 以ERVI所构建的模型, 其对数模型和二次函数模型建模集R2分别为0.60和0.62, 验证集R2均为0.92, RMSE分别为2.40和2.48, 其模型精度都相对较高, 基于R2最大, RMSE最小的原则, 表明ERVI对该研究区SSC的敏感性较强, 其与SSC所构建模型精度较高。

表5 土壤盐分含量反演模型 Table 5 Soil salt content inversion model
2.5 模型验证

用剩下的12个实测值, 对上述所选最佳模型进行验证, 其中以ERVI所建立的对数回归模型验证集R2为0.92, RMSE为2.40; 二次函数回归模型验证集R2为0.92, RMSE为2.48, 根据以上结果剔除ERVI所建立的一元线性回归模型, 然后比较对数、 二次函数回归模型建模集精度与验证集精度[见图3(a)和(b)], 可知对数函数模型的建模集R2为0.597, 验证集RPD为2.181; 二次函数模型的建模集R2为0.624, 验证集RPD为2.089, 根据R2最大, RMSE最小, RPD> 2的原则, 确定二次函数模型为SSC的最佳反演模型。

图3 模型预测值和实测值线性拟合图
(a): 对数模型; (b): 二次函数模型; (c): 预测值和实测值拟合图
Fig.3 Linear fit of measured values and predicted values with the model
(a): Logarithmic model; (b): Quadratic function model; (c): Linear fit of measured values and predicted values

从图3(b)可以看出, 二次函数模型的验证点基本上都分布在对角线附近, 该模型建模集R2为0.624, 验证集R2为0.916, RPD为2.089, 除个别异常点偏离对角线附近外, 其他验证点的线性拟合程度均较好, 因此所构建的SSC模型的预测能力通过了检验, 在一定程度上可以用来反演该地区SSC, 同时该模型为了保持建模集与验证集数据分布的大致一致性, 在建模集和验证集中都存在一个SSC相对较高的值, 使其造成建模集与验证集精度差异较大; 若将验证集的SSC最高值去掉, 见图3(c), 验证集R2将达到0.703, RPD为0.938, 由于RPD< 1.4, 表明该模型预测能力较差, 由此可以说明验证集与预测集数据分布的一致性对模型的预测是很有必要的。

3 结 论

以艾比湖湿地自然保护区为研究区, 利用野外实测地表表层SSC, 并结合从多光谱遥感影像中提取出的各种多光谱遥感指数, 选择出相关性较高的指数将其与实测点盐分值反演的盐分分布图进行定性分析, 并构建统计回归模型来预测SSC, 探究最佳多光谱遥感指数最佳模型在预测SSC的适用性, 结论如下:

(1) 从多光谱遥感影像中提取出SI, TVI, EVI, 分别将其与SSC进行相关性分析比较可知, EVI与SSC的相关性最高, 其相关系数平均比TVI高0.15左右, EVI和SSC的相关系数分别为0.70, 0.70和0.67, SI与SSC的相关性最低。

(2) 选择出与SSC相关性较好的ERVI, ENDVI和EDVI, 将其与实测点盐分值反演的盐分分布图进行定性分析, 可以发现EVI分布图与盐分分布图具有空间一致性, 表明EVI对该区域土壤盐分具有一定的敏感性。

(3) 将SSC与多光谱遥感指数进行建模, 基于R2最大, RMSE最小和RPD> 2的原则比较分析, 可知ERVI所建的二次函数模型的R2为0.92, RMSE为2.48和RPD为2.09, 该模型满足R2最大, RMSE较小和RPD> 2的原则, 表明模型精度较高, 可以用来预测SSC。

本研究选用多光谱遥感影像构建的EVI, 首次用于预测艾比湖湿地自然保护区SSC, 且该反演模型精度较高, 是因为在TVI公式中加入b6和b7波段重新组合得到EVI, 使其具有较为全面的植被、 土壤光谱信息, 从而使SSC反演模型的稳定性得到明显改进, 精度较高, 对该研究区盐渍土的治理与监测尤为重要。 前人也验证过短红外波段可用来反映植被的光谱信息, 用光谱范围(2 100~2 300 nm)之间的波段来构建植被指数, 能够反映出某区域植被的生长状况差异和覆盖度[11], 根据植被的生长状况间接地反演出该地区的土壤盐分。 而且, 许多关于土壤的研究都证明了SSC对2 100~2 400 nm范围的光谱都具有显著影响, 该光谱范围可认为是SSC的典型谱段[12, 13]。 还有研究利用2 153~2 254和1 941~2 092 nm数据来计算光谱指数、 构建土壤盐分的反演模型, 来预测某区域SSC, 取得了理想的效果[14]。 但是在该研究中出现了建模与验证精度差异较大的情况, 主要是采集的土壤盐分异常点较多造成的, 作者从三个方面来讨论异常点较多的原因: 首先, 艾比湖位于阿拉山口大风口通道, 气候干燥, 常年大风, 发生沙尘暴机率增加, 导致该研究区土壤表层盐分分布不均匀; 其次在5月中旬, 日照时间长, 水分蒸发量大, 致使土壤盐分随之向上蒸发, 聚集在地表, 导致该区域土壤表层盐分分布不均; 另外, 艾比湖土壤类型的不同, 也是造成异常点较多的一个重要因素。 砂土土壤颗粒间孔隙较大, 持水性较差, 盐分会随着水分的蒸发聚集在地表, 地表土壤盐分较多, 粘土质地粘重, 土粒之间孔隙小, 透水性差, 土壤盐分聚集在地下深层, 因此地表土壤盐分较少。 以上原因造成该模型精度较不稳定。

根据前人研究, 张添佑等[15, 16]利用多光谱遥感指数来进行干旱区SSC反演, 均取得较好的反演效果。 该研究区艾比湖湿地自然保护区属于典型干旱区, 利用多光谱遥感指数来进行SSC反演, 得到加入2100~2400nm范围的EVI对该区域SSC更为敏感, 与前人研究[13]结果相符, ERVI所构建的SSC反演模型效果最好, 因此可以用多光谱遥感指数与SSC建立的模型作为干旱区SSC预测、 反演的备选方案。

The authors have declared that no competing interests exist.

参考文献
[1] Netternicht G I, Zinck J A. Remote Sensing of Environment, 2003, 85(1): 1. [本文引用:1]
[2] WANG Fei, DING Jian-li, WEI Yang, et al(王飞, 丁建丽, 魏阳, ). Acta Ecologica Sinica(生态学报), 2017, 37(15): 5007. [本文引用:1]
[3] LUO Chong, WANG Ling, PENG Li, et al(罗冲, 王玲, 彭丽, ). Journal of Shihezi University·Natural Science(石河子大学学报·自然科学版), 2016, 34(6): 779. [本文引用:1]
[4] Alexakis D D, Daliakopoulos I N, Panagea I S, et al. Geocarto International, 2016, 33(4): 321. [本文引用:1]
[5] Gorji T, Sertel E, Tanik A. Ecological Indicators, 2017, 74: 384. [本文引用:1]
[6] Gorji T, Sertel E, Tanki A, et al. Geoderma, 2011, 162: 1. [本文引用:1]
[7] SU Xiang-ming, LIU Zhi-hui, WEI Tian-feng, et al(苏向明, 刘志辉, 魏天锋, ). Research of Soil and Water Conservation(水土保持研究), 2016, 23(3): 252. [本文引用:1]
[8] PENG Xiang, HU Dan, ZENG Wen-zhi, et al(彭翔, 胡丹, 曾文治, ). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2016, 32(11): 167. [本文引用:1]
[9] LU Guang, HAN Mei, WANG Min, et al(路广, 韩美, 王敏, ). Ecology and Environmental Sciences(生态环境学报), 2017, 26(3): 422. [本文引用:1]
[10] WANG Fei, DING Jian-li, WU Man-chun(王飞, 丁建丽, 伍漫春). Transactions of the Chinese Society of Agricultural Engineering(农业工程学报), 2010, 26(8): 168. [本文引用:1]
[11] Steig E J, Ding Q, Battisti D S, et al. Annals of Glaciology, 2012, 53(60): 19. [本文引用:1]
[12] LIU Ya, PAN Xian-zhang, WANG Chang-kun, et al(刘娅, 潘贤章, 王昌昆, ). Acta Pedologica Sinica(土壤学报), 2012, 49(4): 824. [本文引用:1]
[13] LAI Ning, LI Xin-guo, LIANG Dong(赖宁, 李新国, 梁东). Journal of Arid Land Resources and Environment(干旱区资源与环境), 2015, 29(2): 151. [本文引用:2]
[14] Weng Y L, Gong P, Zhu Z L. Pedosphere, 2010, 20(3): 378. [本文引用:1]
[15] ZHANG Tian-you, WANG Ling, WANG Hui, et al(张添佑, 王玲, 王辉, ). Acta Ecologica Sinica(生态学报), 2017, 37(9): 3009. [本文引用:1]
[16] WANG Duo-duo, JIA Wen-xiao, WANG Zhi-bao, et al(王多多, 贾文晓, 王志保, ). Journal of Agricultural Science and Technology(中国农业科技导报), 2018, 20(3): 55. [本文引用:1]