艾比湖流域盐渍土含水量光谱特征分析与建模

张海威,张飞,3†,李哲,井云清

(1.新疆大学资源与环境科学学院,830046,乌鲁木齐;2.新疆大学绿洲生态教育部重点实验室,830046,乌鲁木齐; 3.新疆智慧城市与环境建模普通高校重点实验室,830046,乌鲁木齐)

艾比湖流域盐渍土含水量光谱特征分析与建模

张海威1,2,张飞1,2,3†,李哲1,2,井云清1,2

(1.新疆大学资源与环境科学学院,830046,乌鲁木齐;2.新疆大学绿洲生态教育部重点实验室,830046,乌鲁木齐; 3.新疆智慧城市与环境建模普通高校重点实验室,830046,乌鲁木齐)

盐渍土表层盐分累积与土壤含水量有着直接的关系。为了建立干旱区盐渍土含水量高光谱遥感监测模型,本文以艾比湖流域不同含水量的盐渍土为研究对象,采用光谱反射率变换和多元统计分析(MSAM)方法,对土壤含水量的光谱特征进行分析和建模。结果表明:随着土壤含水量的增加,土壤反射率呈下降趋势;在一定范围内,波长越长,土壤光谱反射率与含水量的相关性越高,其中1 937 nm处的土壤光谱反射率与含水量具有最高的相关性(r=-0.636)。对土壤光谱反射率进行8种光谱反射率变换后,在此基础上利用多元统计方法分析盐渍土的不同含水量与光谱之间的相关性,筛选敏感波段,建立关系模型。得出对数一阶微分(Logarithm First Order Differential)在波长2 024和2 357 nm建立的模型以及均方根一阶微分(Root Mean Square First Order Differential)在波长1 972和2 357 nm建立的模型最优,相关系数r分别为0.894和0.865。基于上述模型作者构建了一种耦合模型,其相关性r=0.926比对数一阶微分模型提高了0.032,比均方根一阶模型提高了0.061;因此,所构建的盐渍土含水量估算模型是可行的,可以为遥感反演提供理论参考,对高光谱遥感反演具有一定意义。

光谱特征;含水量;模型;艾比湖流域

干旱区盐渍土研究一直是土壤学科的热门问题之一。土壤含水量不仅是微观气象学中一个重要的能量平衡参数,而且还是旱情监测的重要参数之一[1],因此土壤含水量对土壤盐渍化的研究具有特殊意义。土壤含水量不仅与次生盐渍化有关,还与盐渍土的盐分运移有着密切的关系,是影响绿洲生态环境稳定性的重要因素之一[2]。国内外一些学者建立光谱与土壤含水量模型,通常利用微波遥感对地表层湿度直接进行数据收集,或利用热量距平植被指数、作物缺水指数等方法,获取地表能量与作物生长状况建立模型[35]。当前分析土壤理化性质与光谱反射率之间关系,主要使用主成分分析、逐步多元线性回归分析、偏最小二乘法回归分析及支持向量机和人工神经网络方法[67]。Lu Ning等[8]采用PHI高光谱图像数据,利用多元回归分析与最小二乘法等方法对土壤盐渍化信息进行了可视化。孙小艳等[9]使用手持便携式光谱仪对小麦的光谱反射率进行研究。夏学齐等[10]对不同理化性质的土壤进行了光谱分析与建模。吴亚坤等[11]通过变换光谱反射率与土壤EC之间构建回归模型,分析变化光谱反射率与EC关系。以上研究均是在野外进行的实验研究,另外也有很多学者开展了室内光谱研究,如刘焕军等[12]对黑土土壤水分进行了光谱响应特征分析与建模,何挺等[13]对土壤水分光谱特征进行对比分析,翁永玲等[14]基于PLSR方法对青海茶卡-共和盆地的土壤盐分进行高光谱遥感反演。他们都是在室内进行土壤含水量研究,均取得一定的成果。目前,利用实测光谱与土壤理化性质建立估算模型,已成为盐渍土研究的主流方法。

笔者以室内实测光谱数据和盐渍土含水量为数据源,提出基于不同模型下的耦合建模思想,弥补单一性建模方法的不足。由于针对新疆典型生态脆弱区土壤含水量的光谱特性及建模研究还比较缺乏;因此,选择新疆艾比湖流域为示范区,以盐渍土为研究对象,选用多元统计分析等方法,建立盐渍土不同含水量估算模型,为治理与监测盐渍土提供技术支持。

1 研究区概况

艾比湖流域地处欧亚大陆腹地,位于E 79°53＇～85°02＇,N 43°38＇～45°52＇,三面环山,远离海洋,位于托里县和乌苏南部。研究区属于典型的中温带干旱大陆气候,以干旱少雨、气温变化剧烈为特点。年平均气温6.6～7.8℃,年降水量为116.0～169.2 mm。近40年来,由于入湖水量的减少,艾比湖湖面常年萎缩,湖面曾缩至不足500 km2,水位下降近23 m,干涸湖地也早已沦为盐漠,成为浮尘天气发源地。植被也依次形成旱生、超旱生、沙生、盐生、水生等多种类型,主要植物类型有梭梭(Haloxylon ammodendron)、胡杨(Populus euphratica)、柽柳(Tamar-rix chinensis)和芦苇(Phragmites australis)等[1516]。

2 材料与方法

2.1 数据的采集与处理

2015年10月,对艾比湖流域进行综合调查,整个流域共布设有53个样点,用GPS对其进行坐标定位同时采集距地面0～10 cm的土壤样品,采集土壤时用铝盒装土样,这样可以防止水分的散失。对采集回来的土样,立即使用便携式地物光谱(Field Spec Hi res,ASD,美国)的人工光源采集每个土样光谱曲线,该仪器光谱范围为350～2 500 nm,光谱分辨率3 nm。容器半径为2.5 cm,垂直探头,探头距离为5 cm。测量步骤为:首先开机预热2～3 min、去除暗电流、利用白板定标获取绝对反射率、容器土样表面刮平、收集光谱曲线5条;每采集5个样点光谱要利用白板定标;后利用光谱仪的后处理软件ASD View Spec Pro对采集的光谱曲线分析,对每个点号的5条光谱数据在Process下的Statistics内的Mean命令进行平均值计算,最后光谱数据的归一化;最后在Origin 8.0软件中进行微分变换与多元统计分析。

将土样放入烘箱,在105℃下烘干12 h,放于干燥器中冷却,然后进行干土样称量,根据土壤含水量测定方法[17]计算出土壤含水量

式中:msc为土壤含水量,g/kg,m为铝盒质量,g;m1为土壤烘干前质量,g;m2为土壤烘干后质量,g。

2.2 光谱数据分析

为更好地分析与展示光谱数据与土壤含水量之间的相关性,对土壤光谱反射率分别进行倒数、对数和均方根等数学变换,然后再进行一阶、二阶微分变换消除干扰。一阶、二阶微分方程[18]分别为:

式中:λi+1、λi-1和λi为波长;R＇λi为波长λi的一阶微分光谱;R″λi为波长λi的二阶微分光谱。将变换后的反射率一阶微分、反射率二阶微分、倒数一阶微分、对数一阶微分、均方根一阶微分与原始反射率作为盐渍土的高光谱指数。

2.3 相关性分析

通过每个波段计算盐渍土的高光谱指数与土壤含水量的相关系数,筛选出盐渍土含水量的敏感波段[19]即

式中:rj为高光谱指数与土壤含水量的相关系数;j为波段;Rij为第i个样品第j波段的反射率值;个土样在第j波_段反射率的平均值;Wi为第i个土样的水分含量,g;为土壤样本含水量的平均值,g;n为盐渍土样品个数。

3 结果与分析

3.1 盐渍土含水量光谱特征分析

通过室内分析发现艾比湖流域盐渍土含水量的范围为2.6%～25.6%。把含水量的范围分为5个层次:0～5%、5～10%、10～15%、15～20%、20～25%,分别从各个层次随机选取一条来分析不同含水量的光谱特征。不同含水量下盐渍土光谱特征曲线趋势和态势很相似,但反射率高低不同。如图1所示,总体来看,光谱曲线从400到1350 nm呈平缓上升趋势,但不同含水量倾斜度不同。其中含水量为10.2%的光谱曲线在380～600 nm范围内,土壤光谱反射率最高。含水量为9.7%的光谱曲线在380～1 350 nm范围内,不符合随着含水量的增加而反射率降低的规律,这可能是由于土壤中的含盐量过高所导致的此现象,在600～1450 nm附近反射率曲线斜率变缓,在1 470和1 950 nm左右各有一个水分反射谷,其中在1 950 nm附近的反射谷比较明显。不同土壤含水量反射谷程度不同,反射率呈先升高后降低的趋势。22.1%、10.2%和9.7%3条不同含水量的光谱曲线,在1 950 nm附近出现很明显的反射谷,且这3条曲线反射率差异很小,而4.6%含水量的曲线差异较高。从整体看来随着含水量的增加反射率降低,不同含水量增减速率不同,这与顾燕等[20]和刘曾媛等[21]的研究一致。

图1 不同含水量盐渍土光谱反射率曲线Fig.1 Spectral reflectance curves of salinized soil with different levels of water content

3.2 敏感波段的选取与分析

在Origin 8.0软件中进行8种数学变换:一阶微分、二阶微分、倒数、倒数一阶微分、对数、对数一阶微分、均方根、均方根一阶微分,之后分别与土壤含水量建立相关关系,可以看出在不同波段上变化差异比较大(图2)。

图2 土壤含水量和光谱指数的相关系数(R为反射率)Fig.2 Correlation coefficients between soil water content and spectral indices(R represents reflectance)

由图2可得,对光谱反射率进行微分变换后,其与盐渍土含水量的相关性要明显高于原始反射率,既有明显正相关,又有负相关,且均通过显著性检验。这由于经过数学变换后,放大某些原始光谱数据中细小的波段信息,所以微分变换后的波段与盐渍土含水量相关性比较高。据此可以筛选出盐渍土含水量的敏感波段,确定相关系数最高值对应的波段(表1)。可见,出各种数学变换筛选出的敏感波段波长均＞1 000 nm;一阶微分、倒数一阶、对数一阶和均方根一阶都有相同的敏感波段2 357 nm,且此敏感波段的相关系数较高。

3.3 建模与分析

根据相关分析结果得到盐渍土含水量的敏感波段,建立相应的土壤含水量的估算模型。由于实验误差或数据异常,53个样本中只有44个样点进行了建模与验证,然后从44个样本中随机选取34个样品构建土壤含水量的估算模型,以敏感波段作为自变量,实测的土壤含水量作为因变量,并选择相关系数r较大的敏感波段,进行多元线性回归分析,建立盐渍土含水量估算模型(表2)。

表1 光谱反射率及其变换形式与土壤含水量敏感波段Tab.1 Spectral reflectance and its transformations,and sensitive bands of soil water content

表2 不同自变量回归模型Tab.2 Prediction models based on different characteristic bands

可见,土壤含水量与所选敏感波段呈良好的线性关系,原始反射率和经各数学变换后构建的模型均通过了0.05显著性检验。其中原始反射率的回归模型回归系数最低(r=0.618)。一阶微分回归模型回归系数最高(r=0.858)。依据回归系数和线性回归方程显著性检验统计量F值最大、均方根误差最小的原则,本文选取一阶微分,对数一阶,均方根一阶进行构建模型。

3.3.1 模型验证

模型的验证主要考虑估算模型的稳定性和效果。利用剩下的10个真实值作为验证数据,对上述所选模型对其进行验证,对反演模型的估算效果进行检验其中一阶微分模型精度r=0.554,对数一阶微分模型精度r=0.894,均方根一阶微分精度r= 0.865。根据以上结果剔除一阶微分模型,选取对数一阶微分和均方根一阶微分模型进行盐渍土含水量的反演。从模型的估算值和真实值的线性拟合(图3)可以看出,对数一阶模型样点基本上聚集在对角线附近,样本相关系数r=0.894。相比之下均方根一阶样点没有对数一阶样点那样紧凑,样点相关系数r=0.865。除个别异常点,估算值和实测值的线性拟合程度均较好。本文建立的盐渍土含水量遥感监测模型估算能力通过了检验,在一定程度上可以用来估算盐渍土含水量。

3.3.2 多元线性回归与构建的耦合估算模型

与原始反射率相比,不同变换后的光谱波段与含水量相关性有着不同程度的提高。经过以往研究模型的验证可知:构建的模型比较单一性,都是一种数学变换下的模型构建,有一定的不足;因此本文基于以上模型选取不同变换下的敏感波段构建模型。为了保证模型数据的一致性,模型选取的样本和上述模型一致。构建模型为y=-6.147x1+53.535x2-181.344x3+0.044,式中:y为盐渍土含水量,x1为倒数一阶变换下2 032 nm处的光谱反射率,x2为对数一阶变换2 357 nm处的光谱反射率,x3为均方根一阶变换2 357 nm处的光谱反射率。通过10个样点对模型进行精度验证,可以得到,模型的精度有所提高,相关系数r=0.926,比对数一阶微分模型相关性r提高了0.032,比均方根一阶微分模型相关系数如提高了0.061(表3)。

图3 模型估算值与实测值线性拟合图(N为样本数,下同)Fig.3 Linearfit of measured values and predicted values with the model(N:number of samples,the same below)

表3 不同回归模型对比Tab.3 Comparison of different regression models

相比之前的模型,此模型基于对数一阶微分模型与比均方根一阶微分模型的相关性最优波段,且经过验证发现耦合模型对盐渍土含水量的估算是可行的(图4)。

图4 模型估算值与实测值线性拟合图Fig.4 Linearfit of measured values and predicted values with the model

4 结论与讨论

1)盐渍土不同土含水量光谱曲线的形态与趋势基本一致;在可见光波段范围内,不同含水量的盐渍土反射率呈规律性变化,即随着土壤含水量的增加,土壤反射率呈下降走势。

2)对盐渍土含水量与光谱建模与验证得:经建模与验证后得出对数一阶的模型最好(r=0.804)其次是均方根一阶(r=0.798)。对数一阶微分模型的估算值与实测值的相关系数最高(r=0.894),其次是均方根一阶微分模型,相关系数为r=0.865,即盐渍土含水量模型估算良好。

3)基于均方根一阶微分模型与对数一阶微分模型建立耦合模型,其方程:y=-6.147x1+ 53.535x2-181.344x3+0.044,其模型效果良好(r= 0.926)。证明了选取不同变换后的敏感波段进行建模,对盐渍土含水量估算是可行的。

笔者提出的基于均方根一阶微分模型与对数一阶微分模型下的耦合模型,对艾比湖流域盐渍土的治理与监测有着重要意义和作用;但本研究只是在室内的一种尝试,并没有在室外进行验证,因此耦合模型有待于进一步检验与改进。

[1] 李美婷,武红旗,蒋平安,等.利用土壤的近红外光谱特征测定土壤含水量[J].光谱学与光谱分析,2012, 32(8):2117. LI Meiting,WU Hongqi,JIANG Pingan,et al.Measuring soil water content by using near infrared spectral characteristics of soil[J].Spectroscopy and Spectral Analysis, 2012,32(8):2117.

[2] 中国土壤学会.中国土壤学在前进[M].北京:中国农业科技出版社,1995:9. Soil Science Society of China.Progress of soil science in China[M].Beijing:China's Agricultural Science and Technology Press,1995:9.

[3] ROBINOVE C J,CHAVEZ P S,GEHRING D,et al. Arid land monitoring using Landsat albedo difference images[J].Remote Sensing of Environment,1981,11:133.

[4] HENRICKSEN B L.Reflections on drought:Ethiopia 1983-1984[J].International Journal of Remote Sensing,1986,7(11):1447.

[5] 刘志明,张柏,晏明,等.土壤水分与干旱遥感研究的进展与趋势[J].地球科学进展,2003,18(4):576. LIU Zhiming,ZHANG Bai,YAN Ming,et al.Some research advances and trends on soil moisture and drought monitoring by remote sensing[J].Advance in Earth Sciences,2003,18(4):576.

[6] LIU Huanjun,ZHANG Yuanzhi,ZHANG Bai.Novel hyper spectral reflectance models for estimating black-soil organic matter in Northeast China[J].Environmental Monitoring and Assessment,2009,154(1/2/3/4): 147.

[7] LI Y,DEMETRIADES-SHAH T H,KANEMASU E T, et al.Use of second derivatives of canopy reflectance for monitoring prairie vegetation over different soil backgrounds[J].Remote Sensing of Environment,1993,44 (1):81.

[8] LU Ning,ZHANG Zhi,GAO Yang.Recognition and mapping of soil salinization in arid environment with hyper spectral data[C].Proceedings.2005 IEEE International Geo science and Remote Sensing Symposium,2005. IGArSS'05.IEEE,2005,6:4520.

[9] 孙小艳,常学礼,张宁,等.不同取样单元对干旱区绿洲小麦地上生物量光谱估算模型的影响[J].中国沙漠,2012,32(2):568. SUN Xiaoyan,CHANG Xueli,ZHANG Ning,et al.Impact of different sampling units on ground spectral modelsfor estimating wheat aboveground biomass in oases of arid zone[J].Journal of Desert Research,2012,32(2):568.

[10] 夏学齐,季峻峰,陈骏,等.土壤理化参数的反射光谱分析[J].地学前缘,2009,16(4):354. XIA Xueqi,JI Junfeng,CHEN Jun,et al.Analysis of soil physical and chemical properties by reflectance spectroscopy[J].Earth Science Frontiers,2009,16 (4):354.

[11] 吴亚坤,杨劲松,李晓明.基于光谱指数与EM38的土壤盐分空间变异性研究[J].光谱学与光谱分析, 2009(4):1023. WU Yakun,YANG Jinsong,LI Xiaoming.Study on spatial variability of soil salinity based on spectral indices and EM38 readings[J].Spectroscopy and Spectral Analysis,2009,29(4):1023.

[12] 刘焕军,张柏,宋开山,等.黑土土壤水分光谱响应特征与模型[J].中国科学院研究生院学报,2008,25 (4):503. LIU Huanjun,ZHANG Bai,SONG Kaishan,et al.Soil moisture spectral analysis and its spectral model[J]. Journal of the Graduate School of the Chinese Academy of Sciences,2008,25(4):503.

[13] 何挺,王静,程烨,等.土壤水分光谱特征研究[J].土壤学报,2006,43(6):1027. HE Ting,WANG Jing,CHENG Ye,et al.Soil water spectrum characteristic research[J].Acta Pedologica Sinica,2006,43(6):1027.

[14] 翁永玲,戚浩平,方洪宾,等.基于PLSR方法的青海茶卡-共和盆地土壤盐分高光谱遥感反演[J].土壤学报,2010,47(6):1255. WENG Yongling,QIE Haoping,FANG Hongbin,et al. PLSR based hyper spectral remote sensing retrieval of soil salinity of Chaka·Gonghe basin in Qinghai province[J].Acta Pedologica Sinica,2010,47(6):1255.

[15] 姜红涛,特依拜,克里木,等.艾比湖流域NDWI变化及其与降水、温度的关系[J].中国沙漠,2014,34 (6):1678. JIANG Hongtao,TIYIP,KELIMU,et al.Responses of NDVI to the variation of precipitation and temperature the Ebinur Lake Basin[J].Journal of Desert Research, 2014,34(6):1678.

[16] 张小龙.艾比湖流域气候变化及其径流响应[J].盐湖研究,2011,19(2):11. ZHANG Xiaolong Climatic variation and runoff response in Ebinur Lake Catchment[J].Journal of Salt Lake ResearcH,2011,19(2):11.

[17] 王昌佐,王纪华,王锦地,等.裸土表层含水量高光谱遥感的最佳波段选择[J].遥感信息,2003(4):33. WANG Changzuo,WANG Jihua,WANG Jindi,et al. The choice of best detecting band for hyper spectral remote sensing on surface water content of bare soil[J]. Remote Sensing Information,2003.(4):33.

[18] 关红,贾科利,张至椅,等.盐渍化土壤光谱特征分析与建模[J].国土资源遥感,2015,27(2)100. GUAN Hong,JIA Keli,ZHANG Zhinan.Research on remote sensing monitoring model of soil salinization based on spectrum characteristic analysis[J].Remote Sensing for Land and resources,2015,27(2):100.

[19] 赵振亮,特依拜,丁建丽,等.新疆典型绿洲土壤电导率和PH值的光谱响应特征[J].中国沙漠,2013,33 (5):1413. ZHAO Zhenliang,TIYIP,DING Jianli,et al.Characteristics of spectral responding to soil electrical conductivity PH in the typical Oasis of Xinjiang[J].Journal of Desert Research,2013,33(5):1413.

[20] 顾燕,张鹰,李欢.基于实测光谱的潮滩土壤含水量遥感反演模型研究[J].湿地科学,2013,11(2): 167. GU Yan,ZHANG Ying,LI Huan,et al.Remote sensing retrieval model on soil moisture content of tidal flat based on measured spectra[J].Wetland Science, 2012,11(2):168.

[21] 刘曾媛,贡璐,张雪妮,等.克里雅河流域于田绿洲土壤酶活性c理化因子相关性分析[J].中国土壤与肥料,2014,4(7):35. LIU Zeng yuan,GONG Lu,ZHANG Xueni,et al.Soil enzyme activities and their relationship with physicochemical factors in Yutian oasis of the Keriya river watershed[J].Soil and Fertilizer Sciences in China,2014, 4(7):35.

Modeling and analysis of spectral characteristic of soil water content in the salinized soil of Ebinur Lake Watershed

ZHANG Haiwei1,2,ZHANG Fei1,2,3,LI Zhe1,2,JING Yunqing1,2
(1.College of Resources&Environmental Science,Xinjiang University,830046,Urumqi,China; 2.Key Laboratory of Oasis Ecology,830046,Urumqi,China;3.General Institutes of Higher Learning Key Laboratory of Smart City and Environmental Modeling,Xinjiang University,830046,Urumqi,China)

[Background]With the development of spatial information science,the hyperspectral remote sensing becomes more and more important in nowadays.Studying on spectral characteristic of soil water content is an important work,for it is the base of monitoring remote sensing.This study aims to investigate the spectral characteristic of soil water content and the relationship between the hyperspectral data and the soil water content,to explore a rapid and accurate method for estimating soil water contents, and to establish the hyperspectral remote sensing monitoring model for saline soil water in the arid area of the Ebinur Lake Watershed.[Methods]This study took saline soil with different water contents in Ebinur Lake Watershed as the research object,used the spectral reflectance transformation and multivariate statistical analysis methods(MSAM)to analyze the spectral characteristic of soil water content,and built models.[Results]The result showed that with the increase of soil water content,thereflectance of soil declined.In a certain range,the longer the wavelength,the higher the correlation between the soil spectral reflectance and soil water content;the soil spectral reflectance at the wavelength of 1 937 nm(r=-0.636)had the highest correlation with water content.After 8 transformation of soil spectral reflectance,the correlation of logarithmic first order differential sensitive band of 2 357 nm was the best(r=-0.808 6).Subsequently,MSAM were applied to analyze the correlation between spectrum and salinized soil with different water contents,then the spectral sensitive band of the soil was screened,and the correlation model was established.The result indicated that the model set up by Logarithm First Order Differential at the wavelength of 2 024 nm and 2 357 nm and the Root Mean Square First Order Differential at the wavelength of 1 972 nm and 2 357 nm was the best,the correlation coefficientrwas 0.894 and 0.865 respectively.Based on the above established models,the authors constructed a new model,and the correlationrof which was 0.926,increased 0.032 against the Logarithm First Order Differential model,and increased 0.061 compared to the Root Mean Square First Order Differential model.[Conclusions]Therefore,the estimation model of soil water content is feasible.In addition,this study provides a new model for the indoor hyperspectral estimation of soil water content in Ebinur Lake Watershed,and it also provides a theoretical and technical reference for the hyperspectral quantitative estimation of soil water content,which has certain guiding significance for hyperspectral remote sensing.

spectrum characteristic;water content;modeling;Ebinur Lake Watershed

O433.4

:A

:2096-2673(2017)01-0008-07

10.16843/j.sswc.2017.01.002

2016- 07- 03

2016- 11- 27

项目名称:国家自然科学基金“新疆联合基金本地优秀青年人才培养专项”(U1503302)

张海威(1990—),男,硕士。主要研究方向:干旱区遥感应用。E-mail:1522634426@qq.com

†通信作者简介:张飞(1980—),男,博士,副教授。主要研究方向:干旱区资源环境遥感。E-mail:zhangfei3s@163.com