第33卷,第lO期20l3年10月
光谱学与光谱分析
SpectroscopyandSpectralAnalysis
Voi.33,No.10,pp2771—2776
October,2013
土壤湿润条件下基于光谱对称度的盐渍土盐分含量预测
刘
娅1’2,潘贤章¨,王昌昆1’2,李燕丽1’2,石荣杰1’2,周
睿3,解宪丽1
I.中国科学院土壤环境与污染修复重点实验室(南京土壤研究所),江苏南京2100082.中国科学院大学,北京100049
3.中国科学院南京分院,江苏南京210008
摘要近年来光谱技术以其经济、高效的优势在土壤盐渍化监测研究中得到重视,但是由于土壤水分对反射光谱影响很大,土壤湿润条件下监测精度难以满足农业生产需求。通过对盐土土柱室内模拟蒸发过程中的反射光谱和水分、盐分变化的连续监测,利用多元逐步回归方法,建立了1370~1
610
nlTl光谱对称度
与土壤表层含盐量、含水量之间的线性关系模型,r为0.863;用该模型反演表层土壤含盐量,实测值与预测值之间线性关系的r为0.656(n=54),RMSE为2.059g・kg~。利用光谱对称度可以实现土壤湿润条件下土壤盐分含量预测。
关键词土壤盐渍化;土壤湿润条件;光谱对称度;土壤含盐量;预测中圈分类号:S127。TP79
文献标识码:A
DOI:10.3964/j.issn.1000-0593(2013)10—2771—06
降,1400和1
900
Iun水分子的振动特征变宽且更不对称[1]。
引言
土壤盐渍化是最常见的土地退化问题之一[1],不仅影响
还有研究显示,土壤饱和浸提液电导率与800~810啪处吸
收峰深度显著正相关,而与宽度和不对称度呈显著负相关[7]。Farifteh对不同含盐量的土壤光谱特征进行了研究,认为盐分含量的增加会导致大于1
300am(1
农业生产,对区域生态环境的稳定性也构成一定威胁。盐渍400和1
900
土中含有大量可溶性盐,对作物生长有较强的抑制作用,使
得该类土地生产力很低[z]。由于极易受气候、地质、地形、海岸等影响,中国盐渍土具有分布广泛、类型复杂、季节性明显等特点,因而,给盐渍化监测带来了一定困扰[3]。传统
nm左右)波段反射率变化[8]。刘焕军也指出,盐碱土在400
~2500nnl范围内有五个明显的吸收谷(500,670,1418,
1
915和2
210
nlTl附近),盐渍土由于富含盐分,在上述五个
nlTl去
波段的吸收特征均显著强于其他土壤[9J。Weng对黄河三角洲的盐渍土光谱特征进行了研究,并用2052和2演,取得了较好结果[103。
上述研究表明了采用光谱技术进行土壤盐渍化研究的可行性,但是已有研究大都以风干土为研究对象,尚未考虑到
203
的盐渍土监测、制图主要依靠野外采样或者野外快速测定,
这就使精细的盐渍化分区制图浪费大量人力、物力、财力和
包络线后反射率构建的盐分特征指数对土壤含盐量进行了反
时间[4’5],且易受到采样时问和采样范围的限制,难以实现
盐渍化程度的实时监测、区域制图的快速更新;而一些新型测定方法,比如磁感大地电导率仪的应用还不够广泛[6]。因此,目前的土壤盐渍化信息获取技术满足不了日益提高的需求,需要一种高效、精准的方法。
光谱技术以其高效、经济的优势在土壤盐渍化监测研究中得到重视,国内外学者已对盐渍化土壤的光谱特征进行了大量研究,Dehann等研究表明,高度盐渍化的土壤在680,
l
水分对盐分及光谱的影响。而水分又是土壤的重要组成部
分,亦对土壤的光谱反射率有很大影响。尤其对盐渍化土壤而言,水是土壤中盐分的溶剂和载体,土壤含水量与土壤含
盐量有着密切关系,因而水盐变化对表层光谱反射率必然有
强烈的影响,结合土壤水分研究盐渍土光谱变化特征意义重大。Wang等突破了以往只针对风干土做研究的局限,通过人为添加水分,制作不同含水量条件下各种盐渍化程度的土壤样品,并对其光谱特征进行了研究[11‘,但是没有考虑到自
180和l
780
nIn相对于非盐渍化土壤有更清晰的吸收特
200
征,且随盐渍化程度的增加,2nnl处一OH的吸收特征
变弱,800nrn处增高,在800~1
300啪光谱曲线的斜率下
收稿日期:2012—12—31,修订日期:2013—03—24
基金项目:国家自然科学基金项目(41071140),中国科学院战略先导性项目(XDA0505050902)资助作者简介:刘娅,1986年生,中国科学院南京土壤研究所博士研究生
*通讯联系人
e-mail:pan=@issa&aG
cl-i
e-mail:liuya@issas.ac.cn
2772
光谱学与光谱分析
第33卷
然条件下水盐的相互作用、相互影响,忽略了整个动态的过程。因此,本工作选择在蒸发过程中,对不同水分条件下盐渍土的光谱变化特征及水盐变化进行研究,以探索盐渍土光谱反射率随土壤水盐变化的规律,为更准确地获取盐渍土信息,以及提高区域盐渍化制图精度提供依据,也为今后土壤盐渍化信息的原位、快速获取提供依据。
直方向夹角为45。。为了消除阴影影响,每个样品采集4个方向的光谱,按90。旋转样品,每个方向采集5条光谱曲线,共20条曲线,以20条光谱曲线算术平均值代表该样品的光谱。每次测量前先做暗电流校正,然后用25射标准参考板进行校正[1
2|。
cm×25
cm的漫反
土壤含盐量数据由中国科学院南京土壤研究所研制的小型钠电极和小型氯电极测定,电极经校正后,参照供试土壤离子组成配制一系列盐分浓度的标准溶液,绘制出盐分传感器的读数值与盐分浓度之问的关系曲线,再由每次测量的土壤含水量推算出水土比,根据水土比计算出土壤含盐量,单位为g・kg~。
土壤水分数据由英国Delta-T公司生产的SM300水分速测仪测得,该仪器测量精度可达±2.5%。根据各个模拟土柱的土壤容重,将体积含水量转化为质量含水量,单位为
g・g一。
1实验部分
l-1样品与方法
实验以不同初始含盐量的土壤为研究对象,定期测定蒸发过程中表层土壤水分及盐分数据,同时采集表层土壤反射光谱,研究从湿润到蒸干过程中,三者之间的关系。供试土壤采自江苏省盐城市滨海地区,采样区内土壤为滨海盐土,滨海盐土是盐土的一个亚类,由沿海地区的盐渍淤泥发育而成,该盐土含盐量较高,尤其是Cl和Na+含量较多,pH比较高。对研究区内23个土壤样品分析发现,Cl占阴离子总量55.35%,Na+占阳离子总量53.41%,可见主要阴阳离子含量符合滨海盐土的特征(表1)。
Table1
Descriptivestatisticsofeightionscontentinthesoilsamples
1.3光谱特征吸收参数构建
将平滑后的光谱曲线进行去包络线处理[13【。包络线去除后,光谱反射率归一化到o~1之间,那些峰值点上的反射率值为1,而非峰值点上的反射率则小于1,因此,特征吸收在去包络线之后的曲线上可以明显地看出来,有利于分析光谱的吸收特征。故光谱特征吸收参数均针对去包络线后的光谱曲线而言。
光谱特征吸收参数(图1)一般为吸收波长的位置(P)、深度(H)、宽度(W)、面积(A)和对称度(S)[14|。用于构建土壤含盐量反演模型的光谱特征吸收参数为对称度,计算公式为
S一粤
Arealeft
(1)
式中,S为对称度,Area蛐。为吸收峰右边面积,Arealeft为吸收峰左边面积。
左右吸收峰面积采用积分法计算而得,首先找到吸收波长位置(P),按照逐波长积分的方法累加求得,由Visual
实验所用基础土壤为已基本脱盐土壤,风干磨细过2n1Ⅱ1筛后,装入深25cm直径10cm底部有孔的PVC土柱中,在距表层5cm处埋设盐分传感器和水分传感器,表面用直尺刮平。共制作模拟土柱5个,分别放人盛有蒸馏水,0.3%,0.6%,0.9%以及1.2%NaCl溶液的浅水容器中,从底部充分润湿,模拟5个盐分梯度的湿润土壤。然后开始定时进行光谱、水分和盐分的测量。蒸发过程在实验室条件下进行,实验室内保持恒温恒湿。整个实验过程共采集光谱数据、水盐数据65组,因有些数据为某个土柱表面出现裂隙后所获,光谱数据有些异常,未将其列入分析范围,实际分析所用数据为54组。1.2数据获取
光谱数据采集采用美国ASD
FieldSpec3
Fig.1
P
Basic编程实现。
Thecharacteristicparametersofspectroscopyabsorptioncharacteristics
Hi-Res地物光
由于350~380与2400~2500[trnt噪音较大,不能反映土壤真实的光谱信息,本研究去除上述波段,以380~2Din为光谱研究范围。1.4模型建立与检验
随着水分的蒸发,土体中的盐分必然被运移到表层,并
400
谱仪在暗室中进行,该光谱仪光谱范围为350~2500和350
~l050nin采样间隔为1.4an2,1000~2500nm采样间隔
为2nm,重采样间隔为1nm。光谱测量采用8。视场角探头,探头到土样表面的距离为40cm,光源为两盏50w的卤素灯,其到样品表面中心的距离为40cm,光源照射方向与垂
第10期
光谱学与光谱分析
程度升高而增加,而2210和2光谱在这波段的吸收深度变浅[1
243
2773
在表层积聚,水盐含量的一系列变化势必引起表层反射光谱的变化,且表层反射光谱的变化与水盐变化过程密切相关。由于表层土壤含盐量变化幅度较大,且分布不符合正态分
nIn则反之,这是由于土
壤盐渍化使土壤粘土矿物中的结晶程度减弱,从而导致土壤
0】。
布,所以,对表层含盐量进行了对数转换。以代表光谱反射率变化的光谱对称度为因变量,表层土壤含盐量和表层土壤含水量为自变量,采用多元逐步回归(stepwisemultiple
ar
2.1.2蒸发过程中盐渍土光谱变化特征
图3以其中一个土柱为例来说明蒸发过程中土壤光谱反
line—
射率与表层含水量的关系。该土柱观测过程中含水量范围为
4.32%~32.68%。由图3可知,不同含水量的土壤光谱曲线之间形态相似,各个波段光谱反射率均随含水量的增加而降低,含水量越高,光谱曲线1
414,1
regression,SMLR)方法建立三者之间的定量关系。然后
用所建立的回归模型,结合水分含量,根据光谱对称度反推土壤表层含盐量,并对反演结果进行验证。决定系数(eoeffi—
cientof
943和2
210
m附近水
determination,R2)用于判定模型拟和优度,均方根
meansquare
分的吸收谷越明显。将该土壤5个盐分梯度的观测数据按表层含水量从高到低逐条排序,可以发现,随着含水量的减少,各样品去包络线后的光谱曲线在1370~1
610
误差(root度。
error,RMSE)用于评价预测模型的精
nin表现
出极大的差异性,且该波段范围的光谱曲线因土壤含盐量的
2结果与讨论
2.1盐渍化土壤光谱特征2.1.1盐渍化土壤光谱吸收特征
按照中国盐渍土的等级划分标准[3],以在非盐渍化、轻度盐渍化、中度盐渍化、强度盐渍化和盐土五个级别内的5个风干土壤的光谱曲线为例进行说明(图2)。由5个等级土
差异也呈现出不同的形状(图4所示)。当含水量在30%以上
时,该波段范围的光谱曲线平滑,且吸收峰没有表现出明显
的不对称性;当含水量降低到10%~30%时,该处吸收深度变浅,不对称性有所表现;当含水量降低到3%~10%时,吸收深度继续变浅,且该波段范围开始呈现比较尖锐的形状,吸收峰的不对称性表现明显;含水量继续降低,达到3%以下时,吸收深度继续变浅,该处的不对称性更加明显,且因含盐量的不同导致吸收峰的对称度不同,含盐量越高不对称性越强烈。因此,选用1370~1
610
壤的原始光谱曲线[图2(a)]可知,盐渍化等级越高的土壤,其光谱反射率越低,这是由于采样区内土壤为氯化物型滨海
盐土,富含NaCl和MgClz,而氯化物极易吸收空气中的水分,因此,盐渍化程度越高的土壤其光谱反射率越低。
0.5
nm波段范围的对称度
(简称S・Ⅲ~・s・o)来预测土壤含盐量。
0.40O.350.30
o.4
褂0.3呕0.2
0.1
斛0.25亲020
1喇
“015
0.10O.05
O
0
750
l250
l750
2250
750
l250
l750
Wavelength/nm
Wavelength/nm
№3
疼一
2250
Reflectancespectraof
a
simulatedsoil
1.00
columnwithdifferentsoilmoistures
0.95
2.2多元逐步回归模型
当土壤由湿润状态蒸发至风干状态,这个过程中发生变化的不仅是土壤水分,土壤光谱和含盐量也发生着显著的变化,也就是说,光谱的变化不仅由含水量的变化导致,与含盐量的变化也有密切关系,因此,利用多元逐步回归模型建立S。啪~㈨o、含盐量和含水量三者之间的定量关系。结果显示,含盐量和含水量这两个自变量均能进入模型,并能很好地解释S。。,。~㈨。的变化。由表2可知,该模型的各系数均达到显著(n=54),多元逐步回归模型为
S1
370~1
0.90
0.85
0.80
750
1250
l750
2250
WavelengflYnm
吩2
Reflectancespectra(a)andcontinuumremoved
spectra(b)offivesalinitygrades
610一1.574+0.139109SSC十而u.ull(2)
n
n1々
dyjU
本采样区土壤,主要吸收波段的位置[图2(b)]为:489,
668,902,1143,1414,1809,l943,2210和2243
式中,S・s,o~・s,o为光谱对称度值,SSC为表层含盐量,SMC为表层含水量。
图4中浓度1~浓度10为10f盐分梯度,且盐化程度
nrn。其
中1414和1
943
nIn附近的吸收谷,其吸收深度均随盐渍化
2774
光谱学与光谱分析
第33卷
1.OO0.95
1.00
O95
O.90085
0.90
瓣
尝0.80嶙0.75
0.700.65O.60
1350
1400
1450
1500
l550
l600
1650
惭
箧0.85
O.80
O75
1600
1350l40014501500l550
Wavelength/rim
Wavelength/nm
Wavelength/nm
Fig.4
Wavelength/nm
lull
Continuumremovedspectrainthebandof1、,370~l610ofsampleswithdifferentsoilmoistureandsaltcontent
依次增大,百分数为表层土壤含水量范围。
该模型决定系数R2为0.873,其统计学意义表明,所建模型可以很好地表示自变量和因变量之间的关系,而其物理学意义则表示土壤表面光谱的变化与表层土壤含盐量及表层土壤含水量具有定量关系。
Thbk2
利用上述模型的反推模型,将表层土壤含盐量作为因变量,以光谱对称度和表层土壤含水量为自变量进行表层土壤含盐量的预测。结果表明,实测含盐量与反演含盐量具有较好的线性关系,决定系数R2为0.656,RMSE为2.059kg,结果如图5所示。
g・
模型
器警t
Coefficients
同时也采用Weng和Mashimbyrl0H5]提出的特征光谱指
数,对数据进行分析,反演结果均差于本工作基于光谱对称度方法所得反演结果,可能是由于研究区不同,土壤盐分离
si昏
intheSM匝Rmodel
子组成和母质的差异造成了这些光谱特征指数在本区反演效果较差。
3讨论
4
近年来,国内外学者用光谱方法成功进行了土壤有机质、颗粒组成、氧化铁含量、土壤含水量、pH、土壤含盐量、土壤矿物组成等土壤基本理化性质的预测。而在诸多土壤属性中,土壤水分是影响土壤光谱反射率最重要的因素之一,土壤水分增加会导致光谱反射率迅速降低,这一特性可能会掩盖其他土壤属性的光谱特征,因此,土壤水分对各土壤属性光谱预测均有重要影响。对盐渍土盐分光谱特征分析而言,由于水分的变化和盐分的运移有着密切的关系,很难找到不受土壤水分干扰的盐分敏感波段。既然水分的干扰不可避免,且土壤含水量相对土壤含盐量而言较容易获取,将水
,_
O8642O0
2
4
6
8
10
12
14
含盐鞋实测值/g・kg’1
Fig.5
Scatterdiagramofmeasuredsoilsalt
content
分信息引入反演模型应该是一种可以提高反演精度行之有效
的方法。
研究指出,盐渍化程度越高土壤光谱反射率越低[图2(a)],这种现象是对无盐结壳出现的风干研磨土样而言,而
andpredictedsoilsaltcontent
第10期光谱学与光谱分析
2775
有学者研究认为盐渍化程度越高土壤光谱反射率越高[11’“],其研究对象一般为有盐结壳出现的土壤或野外原状土壤,与本研究对象有所差别。图2(b)中盐渍土去包络线后的特殊光谱吸收波段,如489,668,902,1
2
143,1414,1809,1943,
湿润条件下土壤盐渍化程度的预测方法,得出了以下结论:
(1)利用多元逐步回归模型建立光谱对称度S。。,¨s。。、
表层土壤含盐量和土壤含水量三者之间的定量关系,模型决定系数为0.863,结果表明,土壤表面光谱的变化与表层土壤含盐量及含水量具有定量关系。
(2)以回归模型的反推模型预测土壤表层含盐量,取得了较好的预测结果,R2为0.656,RMSE为2.059
g・kg~,
210和2
343
nm,与已有研究结论类似[1…。其中某些波段
有出入,是因供试土壤的盐分类型、母质或所含矿物类型不同导致。
说明该模型具有较高的预测精度,从而改善了土壤湿润条件
4结论
以盐土土柱室内模拟蒸发实验为依据,以室内光谱数据及所测水盐数据为研究对象,重点研究了土壤湿润条件下反射光谱的变化特征及其与水盐变化的定量关系,探索了土壤
下光谱预测模型的预测精度。
由于本工作主要以江苏滨海地区土壤为基础土壤进行的模拟实验,该模型可能仅适用于该区域或类似区域的盐渍土含盐量预测,其他区域能否运用需要进一步的实验验证。
[1]Dehaan[2]Abrol
zation
R
L,TaylorGRRemoteSensingofEnvironment,2002,80(3):406.
IP,Yadav
JSP,MassoudFI.Salt-AffectedSoilsandTheirManagement:FAO,SoilsBulletin,39.FoodandAgricultureOrgani—
oftheUnitedNations,Rome,1988.
in
[3]WANGZun-qin,ZHUShou-quan,YURen-pei(T遵亲,祝寿泉,俞仁培).SalineSoils
京:科学出版社),1993.
[4]Milton
E
China(中国盐渍土).Beijing:SciencePress(北
J.InternationalJournalofRemoteSensing,1987,8(12):1807.
R,Dematte
[5]NanniMJA
M
soilScience
SocietyofAmericaJournal,2006,70(2):393.
Engineering
[6]YANGJin-song,YAORong-jiang(杨劲松,姚荣江).Hydrogeology8L[7]Shrestha
D
Geology(水文地质工程地质),2007,34(5):67.
P,MargateDE,MeerF,eta1.InternationalJournalofAppliedEarthObservationandGeoinformation,2005,7(2):85.
[8]FariftehJ,MeerF,MeijdeM,eta1.Geodei-ma,2008,145(3—4):196.
[9]LIUHuan-jun,ZHANGBai,WANGZong-ming,et
毫米波学报),2008,27(2):138.
[103[113
WengYongling,GongWangQ,Li
Peng,Zhu
al(刘焕军,张柏,王宗明,等).JournalofInfraredandMillimeterWaves(红外与
Zhiliang.Pedosphere,2010,20(3):378.
PH,ChenXGeoderma,2012,170:103.
Chang-kun,etal(刘娅,潘贤章,王昌昆,等).ActaPedologi∞Sinica(土壤学报),2012,49(4):824.
[12]LIU
Ya,PANR
Xiarrzhang,WANG
[13]Clark[14]Green
N,RoushTL
GeophysRes,1984,89(B7):6329.
ofthe
AA,Craig
MD.Proceedings
AirborneImagingSpectrometerDataAnalysisWorkshop.NASA-JPLPublication,1985,85
(41):111.
[15]Mashimbye
Z
E,ChoMA,NellJP,eta1.Pedosphere,2012,22(5):640.
D
[16]EverittJH,Escobar
E,GerbermannAH,eta1.PhotogrammetryEngineeringandRemoteSensing,1988,54(9):1283.
2776
光谱学与光谱分析
第33卷
PredictingSoilSalinityBased
on
SpectralSymmetryunder
Wet
Soil
Condition
LIUYal~,PAN
Xian-zhang¨,WANG
Chang-kunl~,LIYan-lil~,SHIRong-jiel一,ZHOU
Rui3,XIEXian-lil
1.KeyLaboratoryofSoilEnvironmentandPollutionRemediation,Instituteof
SoilScience,ChineseAcademyofSciences,
Nanjing3.Nanjing
Abstract
210008,China
2.UniversityofChineseAcademyofSciences,Beijing
100049,China
210008,China
BranchoftheChineseAcademyofSciences,Nanjing
Therehasbeen
a
growinginterestinusingspectralreflectance
a
asa
rapidandinexpensivetoolforsoilsalinitymonito—
on
ringinrecentyears.However,sincesoilmoistureoftenexerts
cy
tremendousinfluencesoilreflectance,themonitoring
accura—
undervariousmoistureconditionscannotfullysatisfytherequirementsofagriculturalpractice.Inthepresentpaper,alinear
nlTl
modelwasbuilttOrelatethespectralsymmetryinthebandof1370~1610salinesoilbased
on
withthesalt
content
content
andmoisturecontentofthe
regularlymeasureddataofreflectance.soilmoistureandsaltofthesurfaceof5soilcolumnsduring
r
thesimulatedevaporationprocessinlaboratory.The
mean
resultsshowedthatthemodelwasgoodwithgreaterthan0.8.Byinver—
singthemodel。soilsaltcontentthenwaspredictedaftermoisturecontentwasdetermined.Theresultsshowedthatthepredic—tionaccuracywasacceptablewith
a
root
square
error(RMSE)of2.059
g・k-1
and
anr
of0.656.Theresultsdemon—
stratedthefeasibilityofusingspectralsymmetrytopredictsoilsaltcontentundervariousmoistureconditions.
Keywotds
Soil
salinity;Wetcondition;Spectralsymmetry;Soilsaltcontent;Prediction
*Corresponding
author
(ReceivedDec.31,2012;acceptedMar.24,2013)
(上接2765页)重要时间
论文截稿日期:2014年4月30日第三轮会议通知:2014年10月上旬产品展示
第二轮会议通知:2014年8月上旬会议召开期:2014年10月31—11月4日
sinospectroscopy.org.en.
请您经常浏览光谱网主页,了解会议筹备情况和会议具体安排。网址:http://w、^M
会议热忱邀请国内外仪器厂商前来展示红外、近红外、拉曼、荧光更各类光谱仪以及质谱、核磁等其他相关仪器设备。我
们将在本次会议的网站和会议现场提供展出场所,希望各厂商充分利用这次机会展示自己的最新产品。联系方式
大家若有对召开本次会议的建议和想法,可与组委会联系,若您决定参加会议,请填写回执寄回(联系方式见回执表格),或以电子邮件传到组委会电子邮箱jlyao@suda.edthen,组委会将适时开通网上参会注册系统,请注意组委会的通知。
联系人:姚建林:0512—65880359,13912612398
徐敏敏:0512-65882890,13814802585
jlyao@.suda.edu.cn
x呦@suda.edu.en
中国光学学会中国化学会苏州大学
苏州市化学化工学会
2013年8月
第33卷,第lO期20l3年10月
光谱学与光谱分析
SpectroscopyandSpectralAnalysis
Voi.33,No.10,pp2771—2776
October,2013
土壤湿润条件下基于光谱对称度的盐渍土盐分含量预测
刘
娅1’2,潘贤章¨,王昌昆1’2,李燕丽1’2,石荣杰1’2,周
睿3,解宪丽1
I.中国科学院土壤环境与污染修复重点实验室(南京土壤研究所),江苏南京2100082.中国科学院大学,北京100049
3.中国科学院南京分院,江苏南京210008
摘要近年来光谱技术以其经济、高效的优势在土壤盐渍化监测研究中得到重视,但是由于土壤水分对反射光谱影响很大,土壤湿润条件下监测精度难以满足农业生产需求。通过对盐土土柱室内模拟蒸发过程中的反射光谱和水分、盐分变化的连续监测,利用多元逐步回归方法,建立了1370~1
610
nlTl光谱对称度
与土壤表层含盐量、含水量之间的线性关系模型,r为0.863;用该模型反演表层土壤含盐量,实测值与预测值之间线性关系的r为0.656(n=54),RMSE为2.059g・kg~。利用光谱对称度可以实现土壤湿润条件下土壤盐分含量预测。
关键词土壤盐渍化;土壤湿润条件;光谱对称度;土壤含盐量;预测中圈分类号:S127。TP79
文献标识码:A
DOI:10.3964/j.issn.1000-0593(2013)10—2771—06
降,1400和1
900
Iun水分子的振动特征变宽且更不对称[1]。
引言
土壤盐渍化是最常见的土地退化问题之一[1],不仅影响
还有研究显示,土壤饱和浸提液电导率与800~810啪处吸
收峰深度显著正相关,而与宽度和不对称度呈显著负相关[7]。Farifteh对不同含盐量的土壤光谱特征进行了研究,认为盐分含量的增加会导致大于1
300am(1
农业生产,对区域生态环境的稳定性也构成一定威胁。盐渍400和1
900
土中含有大量可溶性盐,对作物生长有较强的抑制作用,使
得该类土地生产力很低[z]。由于极易受气候、地质、地形、海岸等影响,中国盐渍土具有分布广泛、类型复杂、季节性明显等特点,因而,给盐渍化监测带来了一定困扰[3]。传统
nm左右)波段反射率变化[8]。刘焕军也指出,盐碱土在400
~2500nnl范围内有五个明显的吸收谷(500,670,1418,
1
915和2
210
nlTl附近),盐渍土由于富含盐分,在上述五个
nlTl去
波段的吸收特征均显著强于其他土壤[9J。Weng对黄河三角洲的盐渍土光谱特征进行了研究,并用2052和2演,取得了较好结果[103。
上述研究表明了采用光谱技术进行土壤盐渍化研究的可行性,但是已有研究大都以风干土为研究对象,尚未考虑到
203
的盐渍土监测、制图主要依靠野外采样或者野外快速测定,
这就使精细的盐渍化分区制图浪费大量人力、物力、财力和
包络线后反射率构建的盐分特征指数对土壤含盐量进行了反
时间[4’5],且易受到采样时问和采样范围的限制,难以实现
盐渍化程度的实时监测、区域制图的快速更新;而一些新型测定方法,比如磁感大地电导率仪的应用还不够广泛[6]。因此,目前的土壤盐渍化信息获取技术满足不了日益提高的需求,需要一种高效、精准的方法。
光谱技术以其高效、经济的优势在土壤盐渍化监测研究中得到重视,国内外学者已对盐渍化土壤的光谱特征进行了大量研究,Dehann等研究表明,高度盐渍化的土壤在680,
l
水分对盐分及光谱的影响。而水分又是土壤的重要组成部
分,亦对土壤的光谱反射率有很大影响。尤其对盐渍化土壤而言,水是土壤中盐分的溶剂和载体,土壤含水量与土壤含
盐量有着密切关系,因而水盐变化对表层光谱反射率必然有
强烈的影响,结合土壤水分研究盐渍土光谱变化特征意义重大。Wang等突破了以往只针对风干土做研究的局限,通过人为添加水分,制作不同含水量条件下各种盐渍化程度的土壤样品,并对其光谱特征进行了研究[11‘,但是没有考虑到自
180和l
780
nIn相对于非盐渍化土壤有更清晰的吸收特
200
征,且随盐渍化程度的增加,2nnl处一OH的吸收特征
变弱,800nrn处增高,在800~1
300啪光谱曲线的斜率下
收稿日期:2012—12—31,修订日期:2013—03—24
基金项目:国家自然科学基金项目(41071140),中国科学院战略先导性项目(XDA0505050902)资助作者简介:刘娅,1986年生,中国科学院南京土壤研究所博士研究生
*通讯联系人
e-mail:pan=@issa&aG
cl-i
e-mail:liuya@issas.ac.cn
2772
光谱学与光谱分析
第33卷
然条件下水盐的相互作用、相互影响,忽略了整个动态的过程。因此,本工作选择在蒸发过程中,对不同水分条件下盐渍土的光谱变化特征及水盐变化进行研究,以探索盐渍土光谱反射率随土壤水盐变化的规律,为更准确地获取盐渍土信息,以及提高区域盐渍化制图精度提供依据,也为今后土壤盐渍化信息的原位、快速获取提供依据。
直方向夹角为45。。为了消除阴影影响,每个样品采集4个方向的光谱,按90。旋转样品,每个方向采集5条光谱曲线,共20条曲线,以20条光谱曲线算术平均值代表该样品的光谱。每次测量前先做暗电流校正,然后用25射标准参考板进行校正[1
2|。
cm×25
cm的漫反
土壤含盐量数据由中国科学院南京土壤研究所研制的小型钠电极和小型氯电极测定,电极经校正后,参照供试土壤离子组成配制一系列盐分浓度的标准溶液,绘制出盐分传感器的读数值与盐分浓度之问的关系曲线,再由每次测量的土壤含水量推算出水土比,根据水土比计算出土壤含盐量,单位为g・kg~。
土壤水分数据由英国Delta-T公司生产的SM300水分速测仪测得,该仪器测量精度可达±2.5%。根据各个模拟土柱的土壤容重,将体积含水量转化为质量含水量,单位为
g・g一。
1实验部分
l-1样品与方法
实验以不同初始含盐量的土壤为研究对象,定期测定蒸发过程中表层土壤水分及盐分数据,同时采集表层土壤反射光谱,研究从湿润到蒸干过程中,三者之间的关系。供试土壤采自江苏省盐城市滨海地区,采样区内土壤为滨海盐土,滨海盐土是盐土的一个亚类,由沿海地区的盐渍淤泥发育而成,该盐土含盐量较高,尤其是Cl和Na+含量较多,pH比较高。对研究区内23个土壤样品分析发现,Cl占阴离子总量55.35%,Na+占阳离子总量53.41%,可见主要阴阳离子含量符合滨海盐土的特征(表1)。
Table1
Descriptivestatisticsofeightionscontentinthesoilsamples
1.3光谱特征吸收参数构建
将平滑后的光谱曲线进行去包络线处理[13【。包络线去除后,光谱反射率归一化到o~1之间,那些峰值点上的反射率值为1,而非峰值点上的反射率则小于1,因此,特征吸收在去包络线之后的曲线上可以明显地看出来,有利于分析光谱的吸收特征。故光谱特征吸收参数均针对去包络线后的光谱曲线而言。
光谱特征吸收参数(图1)一般为吸收波长的位置(P)、深度(H)、宽度(W)、面积(A)和对称度(S)[14|。用于构建土壤含盐量反演模型的光谱特征吸收参数为对称度,计算公式为
S一粤
Arealeft
(1)
式中,S为对称度,Area蛐。为吸收峰右边面积,Arealeft为吸收峰左边面积。
左右吸收峰面积采用积分法计算而得,首先找到吸收波长位置(P),按照逐波长积分的方法累加求得,由Visual
实验所用基础土壤为已基本脱盐土壤,风干磨细过2n1Ⅱ1筛后,装入深25cm直径10cm底部有孔的PVC土柱中,在距表层5cm处埋设盐分传感器和水分传感器,表面用直尺刮平。共制作模拟土柱5个,分别放人盛有蒸馏水,0.3%,0.6%,0.9%以及1.2%NaCl溶液的浅水容器中,从底部充分润湿,模拟5个盐分梯度的湿润土壤。然后开始定时进行光谱、水分和盐分的测量。蒸发过程在实验室条件下进行,实验室内保持恒温恒湿。整个实验过程共采集光谱数据、水盐数据65组,因有些数据为某个土柱表面出现裂隙后所获,光谱数据有些异常,未将其列入分析范围,实际分析所用数据为54组。1.2数据获取
光谱数据采集采用美国ASD
FieldSpec3
Fig.1
P
Basic编程实现。
Thecharacteristicparametersofspectroscopyabsorptioncharacteristics
Hi-Res地物光
由于350~380与2400~2500[trnt噪音较大,不能反映土壤真实的光谱信息,本研究去除上述波段,以380~2Din为光谱研究范围。1.4模型建立与检验
随着水分的蒸发,土体中的盐分必然被运移到表层,并
400
谱仪在暗室中进行,该光谱仪光谱范围为350~2500和350
~l050nin采样间隔为1.4an2,1000~2500nm采样间隔
为2nm,重采样间隔为1nm。光谱测量采用8。视场角探头,探头到土样表面的距离为40cm,光源为两盏50w的卤素灯,其到样品表面中心的距离为40cm,光源照射方向与垂
第10期
光谱学与光谱分析
程度升高而增加,而2210和2光谱在这波段的吸收深度变浅[1
243
2773
在表层积聚,水盐含量的一系列变化势必引起表层反射光谱的变化,且表层反射光谱的变化与水盐变化过程密切相关。由于表层土壤含盐量变化幅度较大,且分布不符合正态分
nIn则反之,这是由于土
壤盐渍化使土壤粘土矿物中的结晶程度减弱,从而导致土壤
0】。
布,所以,对表层含盐量进行了对数转换。以代表光谱反射率变化的光谱对称度为因变量,表层土壤含盐量和表层土壤含水量为自变量,采用多元逐步回归(stepwisemultiple
ar
2.1.2蒸发过程中盐渍土光谱变化特征
图3以其中一个土柱为例来说明蒸发过程中土壤光谱反
line—
射率与表层含水量的关系。该土柱观测过程中含水量范围为
4.32%~32.68%。由图3可知,不同含水量的土壤光谱曲线之间形态相似,各个波段光谱反射率均随含水量的增加而降低,含水量越高,光谱曲线1
414,1
regression,SMLR)方法建立三者之间的定量关系。然后
用所建立的回归模型,结合水分含量,根据光谱对称度反推土壤表层含盐量,并对反演结果进行验证。决定系数(eoeffi—
cientof
943和2
210
m附近水
determination,R2)用于判定模型拟和优度,均方根
meansquare
分的吸收谷越明显。将该土壤5个盐分梯度的观测数据按表层含水量从高到低逐条排序,可以发现,随着含水量的减少,各样品去包络线后的光谱曲线在1370~1
610
误差(root度。
error,RMSE)用于评价预测模型的精
nin表现
出极大的差异性,且该波段范围的光谱曲线因土壤含盐量的
2结果与讨论
2.1盐渍化土壤光谱特征2.1.1盐渍化土壤光谱吸收特征
按照中国盐渍土的等级划分标准[3],以在非盐渍化、轻度盐渍化、中度盐渍化、强度盐渍化和盐土五个级别内的5个风干土壤的光谱曲线为例进行说明(图2)。由5个等级土
差异也呈现出不同的形状(图4所示)。当含水量在30%以上
时,该波段范围的光谱曲线平滑,且吸收峰没有表现出明显
的不对称性;当含水量降低到10%~30%时,该处吸收深度变浅,不对称性有所表现;当含水量降低到3%~10%时,吸收深度继续变浅,且该波段范围开始呈现比较尖锐的形状,吸收峰的不对称性表现明显;含水量继续降低,达到3%以下时,吸收深度继续变浅,该处的不对称性更加明显,且因含盐量的不同导致吸收峰的对称度不同,含盐量越高不对称性越强烈。因此,选用1370~1
610
壤的原始光谱曲线[图2(a)]可知,盐渍化等级越高的土壤,其光谱反射率越低,这是由于采样区内土壤为氯化物型滨海
盐土,富含NaCl和MgClz,而氯化物极易吸收空气中的水分,因此,盐渍化程度越高的土壤其光谱反射率越低。
0.5
nm波段范围的对称度
(简称S・Ⅲ~・s・o)来预测土壤含盐量。
0.40O.350.30
o.4
褂0.3呕0.2
0.1
斛0.25亲020
1喇
“015
0.10O.05
O
0
750
l250
l750
2250
750
l250
l750
Wavelength/nm
Wavelength/nm
№3
疼一
2250
Reflectancespectraof
a
simulatedsoil
1.00
columnwithdifferentsoilmoistures
0.95
2.2多元逐步回归模型
当土壤由湿润状态蒸发至风干状态,这个过程中发生变化的不仅是土壤水分,土壤光谱和含盐量也发生着显著的变化,也就是说,光谱的变化不仅由含水量的变化导致,与含盐量的变化也有密切关系,因此,利用多元逐步回归模型建立S。啪~㈨o、含盐量和含水量三者之间的定量关系。结果显示,含盐量和含水量这两个自变量均能进入模型,并能很好地解释S。。,。~㈨。的变化。由表2可知,该模型的各系数均达到显著(n=54),多元逐步回归模型为
S1
370~1
0.90
0.85
0.80
750
1250
l750
2250
WavelengflYnm
吩2
Reflectancespectra(a)andcontinuumremoved
spectra(b)offivesalinitygrades
610一1.574+0.139109SSC十而u.ull(2)
n
n1々
dyjU
本采样区土壤,主要吸收波段的位置[图2(b)]为:489,
668,902,1143,1414,1809,l943,2210和2243
式中,S・s,o~・s,o为光谱对称度值,SSC为表层含盐量,SMC为表层含水量。
图4中浓度1~浓度10为10f盐分梯度,且盐化程度
nrn。其
中1414和1
943
nIn附近的吸收谷,其吸收深度均随盐渍化
2774
光谱学与光谱分析
第33卷
1.OO0.95
1.00
O95
O.90085
0.90
瓣
尝0.80嶙0.75
0.700.65O.60
1350
1400
1450
1500
l550
l600
1650
惭
箧0.85
O.80
O75
1600
1350l40014501500l550
Wavelength/rim
Wavelength/nm
Wavelength/nm
Fig.4
Wavelength/nm
lull
Continuumremovedspectrainthebandof1、,370~l610ofsampleswithdifferentsoilmoistureandsaltcontent
依次增大,百分数为表层土壤含水量范围。
该模型决定系数R2为0.873,其统计学意义表明,所建模型可以很好地表示自变量和因变量之间的关系,而其物理学意义则表示土壤表面光谱的变化与表层土壤含盐量及表层土壤含水量具有定量关系。
Thbk2
利用上述模型的反推模型,将表层土壤含盐量作为因变量,以光谱对称度和表层土壤含水量为自变量进行表层土壤含盐量的预测。结果表明,实测含盐量与反演含盐量具有较好的线性关系,决定系数R2为0.656,RMSE为2.059kg,结果如图5所示。
g・
模型
器警t
Coefficients
同时也采用Weng和Mashimbyrl0H5]提出的特征光谱指
数,对数据进行分析,反演结果均差于本工作基于光谱对称度方法所得反演结果,可能是由于研究区不同,土壤盐分离
si昏
intheSM匝Rmodel
子组成和母质的差异造成了这些光谱特征指数在本区反演效果较差。
3讨论
4
近年来,国内外学者用光谱方法成功进行了土壤有机质、颗粒组成、氧化铁含量、土壤含水量、pH、土壤含盐量、土壤矿物组成等土壤基本理化性质的预测。而在诸多土壤属性中,土壤水分是影响土壤光谱反射率最重要的因素之一,土壤水分增加会导致光谱反射率迅速降低,这一特性可能会掩盖其他土壤属性的光谱特征,因此,土壤水分对各土壤属性光谱预测均有重要影响。对盐渍土盐分光谱特征分析而言,由于水分的变化和盐分的运移有着密切的关系,很难找到不受土壤水分干扰的盐分敏感波段。既然水分的干扰不可避免,且土壤含水量相对土壤含盐量而言较容易获取,将水
,_
O8642O0
2
4
6
8
10
12
14
含盐鞋实测值/g・kg’1
Fig.5
Scatterdiagramofmeasuredsoilsalt
content
分信息引入反演模型应该是一种可以提高反演精度行之有效
的方法。
研究指出,盐渍化程度越高土壤光谱反射率越低[图2(a)],这种现象是对无盐结壳出现的风干研磨土样而言,而
andpredictedsoilsaltcontent
第10期光谱学与光谱分析
2775
有学者研究认为盐渍化程度越高土壤光谱反射率越高[11’“],其研究对象一般为有盐结壳出现的土壤或野外原状土壤,与本研究对象有所差别。图2(b)中盐渍土去包络线后的特殊光谱吸收波段,如489,668,902,1
2
143,1414,1809,1943,
湿润条件下土壤盐渍化程度的预测方法,得出了以下结论:
(1)利用多元逐步回归模型建立光谱对称度S。。,¨s。。、
表层土壤含盐量和土壤含水量三者之间的定量关系,模型决定系数为0.863,结果表明,土壤表面光谱的变化与表层土壤含盐量及含水量具有定量关系。
(2)以回归模型的反推模型预测土壤表层含盐量,取得了较好的预测结果,R2为0.656,RMSE为2.059
g・kg~,
210和2
343
nm,与已有研究结论类似[1…。其中某些波段
有出入,是因供试土壤的盐分类型、母质或所含矿物类型不同导致。
说明该模型具有较高的预测精度,从而改善了土壤湿润条件
4结论
以盐土土柱室内模拟蒸发实验为依据,以室内光谱数据及所测水盐数据为研究对象,重点研究了土壤湿润条件下反射光谱的变化特征及其与水盐变化的定量关系,探索了土壤
下光谱预测模型的预测精度。
由于本工作主要以江苏滨海地区土壤为基础土壤进行的模拟实验,该模型可能仅适用于该区域或类似区域的盐渍土含盐量预测,其他区域能否运用需要进一步的实验验证。
[1]Dehaan[2]Abrol
zation
R
L,TaylorGRRemoteSensingofEnvironment,2002,80(3):406.
IP,Yadav
JSP,MassoudFI.Salt-AffectedSoilsandTheirManagement:FAO,SoilsBulletin,39.FoodandAgricultureOrgani—
oftheUnitedNations,Rome,1988.
in
[3]WANGZun-qin,ZHUShou-quan,YURen-pei(T遵亲,祝寿泉,俞仁培).SalineSoils
京:科学出版社),1993.
[4]Milton
E
China(中国盐渍土).Beijing:SciencePress(北
J.InternationalJournalofRemoteSensing,1987,8(12):1807.
R,Dematte
[5]NanniMJA
M
soilScience
SocietyofAmericaJournal,2006,70(2):393.
Engineering
[6]YANGJin-song,YAORong-jiang(杨劲松,姚荣江).Hydrogeology8L[7]Shrestha
D
Geology(水文地质工程地质),2007,34(5):67.
P,MargateDE,MeerF,eta1.InternationalJournalofAppliedEarthObservationandGeoinformation,2005,7(2):85.
[8]FariftehJ,MeerF,MeijdeM,eta1.Geodei-ma,2008,145(3—4):196.
[9]LIUHuan-jun,ZHANGBai,WANGZong-ming,et
毫米波学报),2008,27(2):138.
[103[113
WengYongling,GongWangQ,Li
Peng,Zhu
al(刘焕军,张柏,王宗明,等).JournalofInfraredandMillimeterWaves(红外与
Zhiliang.Pedosphere,2010,20(3):378.
PH,ChenXGeoderma,2012,170:103.
Chang-kun,etal(刘娅,潘贤章,王昌昆,等).ActaPedologi∞Sinica(土壤学报),2012,49(4):824.
[12]LIU
Ya,PANR
Xiarrzhang,WANG
[13]Clark[14]Green
N,RoushTL
GeophysRes,1984,89(B7):6329.
ofthe
AA,Craig
MD.Proceedings
AirborneImagingSpectrometerDataAnalysisWorkshop.NASA-JPLPublication,1985,85
(41):111.
[15]Mashimbye
Z
E,ChoMA,NellJP,eta1.Pedosphere,2012,22(5):640.
D
[16]EverittJH,Escobar
E,GerbermannAH,eta1.PhotogrammetryEngineeringandRemoteSensing,1988,54(9):1283.
2776
光谱学与光谱分析
第33卷
PredictingSoilSalinityBased
on
SpectralSymmetryunder
Wet
Soil
Condition
LIUYal~,PAN
Xian-zhang¨,WANG
Chang-kunl~,LIYan-lil~,SHIRong-jiel一,ZHOU
Rui3,XIEXian-lil
1.KeyLaboratoryofSoilEnvironmentandPollutionRemediation,Instituteof
SoilScience,ChineseAcademyofSciences,
Nanjing3.Nanjing
Abstract
210008,China
2.UniversityofChineseAcademyofSciences,Beijing
100049,China
210008,China
BranchoftheChineseAcademyofSciences,Nanjing
Therehasbeen
a
growinginterestinusingspectralreflectance
a
asa
rapidandinexpensivetoolforsoilsalinitymonito—
on
ringinrecentyears.However,sincesoilmoistureoftenexerts
cy
tremendousinfluencesoilreflectance,themonitoring
accura—
undervariousmoistureconditionscannotfullysatisfytherequirementsofagriculturalpractice.Inthepresentpaper,alinear
nlTl
modelwasbuilttOrelatethespectralsymmetryinthebandof1370~1610salinesoilbased
on
withthesalt
content
content
andmoisturecontentofthe
regularlymeasureddataofreflectance.soilmoistureandsaltofthesurfaceof5soilcolumnsduring
r
thesimulatedevaporationprocessinlaboratory.The
mean
resultsshowedthatthemodelwasgoodwithgreaterthan0.8.Byinver—
singthemodel。soilsaltcontentthenwaspredictedaftermoisturecontentwasdetermined.Theresultsshowedthatthepredic—tionaccuracywasacceptablewith
a
root
square
error(RMSE)of2.059
g・k-1
and
anr
of0.656.Theresultsdemon—
stratedthefeasibilityofusingspectralsymmetrytopredictsoilsaltcontentundervariousmoistureconditions.
Keywotds
Soil
salinity;Wetcondition;Spectralsymmetry;Soilsaltcontent;Prediction
*Corresponding
author
(ReceivedDec.31,2012;acceptedMar.24,2013)
(上接2765页)重要时间
论文截稿日期:2014年4月30日第三轮会议通知:2014年10月上旬产品展示
第二轮会议通知:2014年8月上旬会议召开期:2014年10月31—11月4日
sinospectroscopy.org.en.
请您经常浏览光谱网主页,了解会议筹备情况和会议具体安排。网址:http://w、^M
会议热忱邀请国内外仪器厂商前来展示红外、近红外、拉曼、荧光更各类光谱仪以及质谱、核磁等其他相关仪器设备。我
们将在本次会议的网站和会议现场提供展出场所,希望各厂商充分利用这次机会展示自己的最新产品。联系方式
大家若有对召开本次会议的建议和想法,可与组委会联系,若您决定参加会议,请填写回执寄回(联系方式见回执表格),或以电子邮件传到组委会电子邮箱jlyao@suda.edthen,组委会将适时开通网上参会注册系统,请注意组委会的通知。
联系人:姚建林:0512—65880359,13912612398
徐敏敏:0512-65882890,13814802585
jlyao@.suda.edu.cn
x呦@suda.edu.en
中国光学学会中国化学会苏州大学
苏州市化学化工学会
2013年8月