土壤水分的遥感监测方法概述
吴黎1, 张有智1, 解文欢1, 李岩1, 宋静波2
1. 黑龙江省农业科学院遥感技术中心,哈尔滨 150086
2. 黑龙江省社会科学院应用经济所,哈尔滨 150000
张有智(1979-),男,硕士,副研究员,主要从事农业遥感研究工作。Email:zhangyouzhi68@sina.com

第一作者简介: 吴 黎 (1983-),女,硕士,助理研究员,主要从事农业遥感应用研究。Email:aromawu@163.com

摘要

回顾了目前国内外土壤水分的遥感监测方法,介绍了反射率法、植被指数法、地表温度法、温度-植被指数法、作物水分胁迫指数法、热惯量法和微波法,并对各方法的优缺点进行了详细比对; 在总结国内外土壤水分遥感监测研究方法的基础上,对目前该研究领域的重点、难点和未来的发展方向进行了评价。认为: 热惯量法和植被温度指数法是较为成熟的方法; 微波遥感因其独特的优越性,将是该领域的重点研究方向。

关键词: 遥感; 土壤水分; 地表温度; 热惯量; 微波
中图分类号:TP79 文献标志码:A 文章编号:1001-070X(2014)02-0019-08 doi: 10.6046/gtzyyg.2014.02.04
Summary of remote sensing methods for monitoring soil moisture
WU Li1, ZHANG Youzhi1, XIE Wenhuan1, LI Yan1, SONG Jingbo2
1. Remote Sensing Technique Center of Heilongjiang Academy of Agricultural Science, Harbin 150086, China
2. Institute of Applied Economic of Heilongjiang Academy of Social Science, Harbin 150000, China
Abstract

A comparison with traditional soil moisture monitoring methods shows that the remote sensing method has great superiority. This paper presents a review of the remote sensing methods currently used both in China and abroad for monitoring soil moisture, which include the reflectivity method, the vegetation index method, the surface temperature, temperature-vegetation index method, the crop water stress index method, the thermal inertia method and the microwave method, with a detailed comparative description of the advantages and disadvantages of these methods. Based on summarizing researches on remote sensing monitoring methods for soil water, this paper evaluated the focal points, difficulties and development trend of this research field. It is held that the thermal inertia method and the vegetation temperature index method are relatively mature methods for soil moisture monitoring. With the wide application of geographic information system, the microwave remote sensing will become the key research direction in this field because of its unique advantages.

Keyword: remote sensing; soil moisture; land surface temperature; thermal inertia; microwave
0 引言

土壤水分是自然界水分平衡的重要参量, 根据土壤-作物-大气连续体(soil-plant-atmosphere continuum, SPAC)理论, 水分在连续体内的运动主要由水势差决定。土壤水分与干旱存在密切关系, 而干旱是我国农业生产的最大威胁。我国农业自然灾害的近60%是干旱造成的, 每年直接经济损失达4~7亿元[1]。土壤水分的变化及区域差异不但对于区域水文, 而且通过植被、土壤等因素对气候产生很大影响。土壤水分监测在农业生产、旱情监测、农田水分管理、水资源管理、水文模拟和预报等许多方面都有重要意义。

由于影响土壤水分的因素很多, 大范围、实时准确地监测土壤水分仍是一个难题。传统监测土壤水分含量的方法有重量法、中子仪法、张力计法、时域反射法等, 这些方法只能采集样点的离散数据, 需要投入大量的人力物力, 耗费长、效率低、成本高、代表性差, 很难完成大范围的土壤水分数据的获取。

遥感监测土壤水分是通过遥感器测量得到土壤表面反射或发射的电磁能量, 然后分析、建立该信息与土壤水分之间的关系从而得到土壤水分含量。目前遥感监测土壤水分主要用到的遥感波段为可见光-近红外、热红外和微波。热红外遥感因其对土壤水分的敏感性, 在土壤水分含量监测中已日益成熟, 并得到了较广泛的应用。微波遥感因其全天时、全天候、多极化、高分辨率、穿透性及对土壤水分含量敏感性的优势, 使其成为该领域未来发展的主要方向。本文在查阅了大量相关资料的基础上, 就土壤水分的遥感监测方法进行了总结、分析、评价和展望。

1 土壤水分的遥感监测方法
1.1 反射率法

土壤在可见光-近红外波段的反射率随波长增加而增大; 对于给定类型的土壤, 同一波段的反射率随土壤水分的增加而减小, 可以建立土壤反射特性与土壤水分的关系来监测土壤水分。Bowers[2]等发现土壤水分的增加会引起土壤反射率的降低, 随后Curran[3, 4, 5, 6]等学者的研究也证实了反射率与土壤水分含量存在密切关系, 并建立了不同模型反演土壤水分。国内, 刘培君[7]等提出“ 光学植被覆盖度” 的概念, 通过分解像元从土壤中分离出植被信息, 建立了利用AVHRR可见光与近红外通道数据监测土壤水分的模型; 张仁华[8]通过地面控制试验发现, 采用不同遥感仪器研究近红外波段反射率与土壤水分之间的关系, 仍然可以得到一致的结论, 并建立了10 cm土层含水量与近红外波段反射率的线性关系模型; 詹志明[9]等提出了一个基于红光和近红外反射率光谱空间特征的土壤水分监测模型, 并进行了验证, 取得了满意的效果。

1.2 植被指数法

植被指数是反映植物生长状况的参数, 是传感器不同通道数据的线性或非线性组合。最常用的植被指数为归一化植被指数(normalized difference vegetation index, NDVI), 即

NDVI=(ρ NIR-ρ RED)/(ρ NIR+ρ RED), (1)

式中, ρ NIRρ RED分别为近红外和红光波段反射率。

在有植被覆盖地区, 土壤水分状况越好植被生长状况越好, 植被指数越高, 所以可建立土壤水分与植被指数之间的关系监测土壤水分。目前已经发展了很多用于监测土壤水分的植被指数模型, 常用的有距平植被指数(anomaly vegetation index, AVI)和条件植被指数(vegetation condition index, VCI)。

AVI定义为

AVI=NDVIi- NDVI¯, (2)

式中: NDVIi为某一特定年份某一时期(旬、月等)的NDVI值; NDVI¯为多年同一时期NDVI的平均值。AVI> 0说明该时期植被生长较一般年份好, AVI< 0则说明较一般年份差, 土壤水分亏缺, 旱情出现。

VCI定义为

VCI=(NCVIi-NDVImin)/(NDVImax-NDVImin) , (3)

式中: NDVIi为特定年份第i个时期的值; NDVImaxNDVImin分别代表所研究年内第i个时期NDVI的最大值和最小值。

国内陈维英[10]等利用NOAA卫星数据计算AVI, 对1992年中国干旱状况进行了监测; 王鹏新等[11]用VCI和AVI 对陕西关中平原的土壤水分进行了监测; 蔡斌等[12]应用 VCI结合常规资料进行综合分析对我国的干旱状况进行宏观动态监测。

植被指数法监测土壤水分方法是建立在计算植被指数的基础上进行的, 所以该方法只能应用于有植被覆盖地区; 在大区域范围应用时, 由于地表覆盖状况非常复杂, 需要先判断是不是满足该方法的应用要求, 再决定可否应用此方法进行监测。

1.3 地表温度法

地表温度(Ts)是大气-土壤-植被系统内物质(水分)和能量(热能)交换共同作用的结果。水的比热大, 接受太阳辐射后温度变化较慢, 白天下垫面温度的空间分布能间接反映土壤水分的分布, 据此建立地表温度与土壤水分的关系来监测土壤水分状况。Mcvicar T R[13]通过对地表温度的归一化处理提出了归一化温度指数(normalized difference temperature index, NDTI)的概念; Kogan[14]提出了条件温度指数(temperature condition index, TCI)的概念并进行了干旱监测。国内刘志明[15]研究发现白天地表温度与土壤水分存在较好的负相关, 而夜晚温度与土壤水分没有确定的相关关系, 并利用该方法对四川省的土壤水分进行了动态监测; 乔平林等[16]利用TM6波段反演得到的地表温度监测了甘肃石羊河流域民勤盆地区土壤水分, 结果表明在某些土壤水分范围内可以获得成功; 张树誉[17]等利用NOAA/AVHRR的热红外通道对陕西省的土壤水分进行了监测; 张仁华[8]通过对地表温度归一化处理得到了一个相对温差模型, 即

Ts=K(Tmax-Tmin)/(Tmax+Tmin), (4)

式中: Ts为地表温度; K为与天气、植被类型、季节等因素相关的系数; TmaxTmin分别为最大和最小亮度温度。研究表明, (Tmax-Tmin)可以较好地反映土壤热通量G; 而(Tmax+Tmin)则较好地表达了土壤热量收入, 并建立了该温差模型与土壤水分的线性回归方程。研究结果表明, 用该模型与土壤水分建立的回归方程比单纯的利用地表温度效果要好。

地表温度对下垫面湿度状况敏感, 反应迅速, 该方法简单易行。但影响地表温度的因素很多, 除土壤水分外, 太阳辐射、显热、潜热、风及植被等因素都对地表温度有很大影响。到目前为止, 遥感反演地表温度并没有非常完美精确的方法; 土壤水分与地表温度的关系也比较复杂, 如何找到两者合适的定量关系也是需要近一步研究。

1.4 温度-植被指数法

土壤水分、地表温度、植被指数3者之间相互影响。单一使用植被指数或者地表温度都存在某些不足, 比如作物类型、地理位置、气候类型等对NDVI有影响, 而地表温度除受土壤水分影响之外, 地表— 大气之间传导作用、混合像元分解精度、发射率等都影响地表温度状况。毛学森[18]发现, 冬小麦在受到水分胁迫时NDVI对土壤水分的反应具有一定的滞后性, 所以将地表温度和植被指数联合起来监测土壤水分效果更好。

1.4.1 温度植被干旱指数法

Goetz[19]通过分析认为NDVITs之间存在明显的负相关关系, 主要原因是植被受到水分胁迫时下垫面温度会急剧升高, 并估算了区域平均的土壤湿度条件, 且认为传感器的分辨率对两者之间的关系没有太大影响; Price等[20]分析了不同卫星传感器得到的NDVITs数据, 认为NDVI和地表辐射温度构成的散点图呈三角形(图1)。

图1 地表温度和植被指数构成的三角形空间[21]Fig.1 Triangular space of vegetation index and ground surface temperature[21]

Nemani[21]从理论上分析认为地表温度和植被指数之间应为梯形关系; Moran[22]等加入空气温度(Ta)数据, 建立植被覆盖度和植被指数之间的线性关系, 定义了植被指数-地表温度梯形图(图2)。

图2 植被覆盖和地表温度构成的梯形图[22]Fig.2 Ladder space of vegetation index and ground surface temperature

通过卫星资料得到区域植被指数和地表温度, 建立二者散点图, 确定干边、湿边和模型各个顶点坐标, 就可以得到区域土壤水分的空间分布。国内外学者对此进行了大量的研究和应用。 Sandholt[23]等在植被指数-地表温度特征空间时发现了很多条直线, 据此提出了温度植被干旱指数(temperature vegetagle dryness index, TVDI)监测地表湿度状况, 即

TVDI= Ts-Ts-minTs-max-Ts-min= Ts-(a2+b2NDVI)(a1+b2NDVI)-(a2+b2NDVI), (5)

式中: Ts为地表温度; Ts-min 为相同NDVI条件下的最小地表温度, 对应Ts-NDVI特征空间的湿边; Ts-max为相同NDVI条件下的最大地表温度, 对应Ts-NDVI特征空间的干边; a1, a2, b1, b2为回归系数。

Lambin[24]研究认为在NDVI-Ts特征空间中, 确定干、湿边后, 可以进一步确定在相同植被覆盖条件下因土壤供水能力差异而造成的温差变化范围, 并以此为条件反演得到土壤水分含量; 韩丽娟[25]等详细解释了NDVI-Ts特征空间, 并用蒸散和温度植被干旱指数解释了此特征空间的内涵; 刘良云等[26]利用两者关系对地物进行分类, 提取了植被覆盖和土壤水分的信息; 姚春生[27]利用MODIS数据得到的TVDI反演了新疆地区2个月的土壤水分; 王鹏新[11]等在NDVI-Ts特征空间和TVDI的基础上, 提出了条件植被温度指数(vegetable temperature condition index, VTCI) 模型, 监测土壤水分和干旱。

1.4.2 植被供水指数法

考虑到植被指数和地表温度与土壤水分的关系, Carlson[28]综合考虑植物受旱时在红光、近红外及热红外波段上的反应, 提出了植被供水指数(vegetation supply water index, VSWI), 即

VSWI=NDVI/Ts。 (6)

VSWI值越小, 说明植被受干旱胁迫, 土壤水分越少。该指数计算简单, 易于实现, 适用于植被覆盖区, 尤其在作物覆盖良好的情况下效果很好。刘丽等[29]应用NOAA数据建立了VSWI监测贵州干旱状况的模型, 认为VSWI方法监测干旱效果较好, 并适合常年植被覆盖较高地区; 莫伟华等[30]使用该方法根据NOAA资料对广西贵港地区进行了干旱遥感监测, 均取得了良好的效果。

1.5 作物水分胁迫指数法

Idso[31]在Monteith等人研究的基础上提出作物水分胁迫指数(crop water stress index, CWSI ), 即

CWSI= (Tc-Ta)-(Tc-Ta)LL(Tc-Ta)UL-(Tc-Ta)LL, (7)

式中: Tc为地表(灌层)温度; Ta为空气温度; (Tc-Ta)LL表示充分湿润条件下的冠气温差下限; (Tc-Ta)UL表示极端干旱条件下的冠气温差上限。

Idso认为冠气温差的上下限可以与空气饱和水气压建立线性方程来计算, 被称为CWSI的经验模式。Jackson[32]在Idso经验模式的基础上, 根据冠层热量平衡方程对冠气温差的上下限进行理论解释, 得到新的计算模型, 即

CWSI=1- ETETp= γ(1+γcγa)-γ(1+γcpγa)Δ+γ(1+γcγa)。 (8)

式中: γ 为干湿表常数; γ a为空气动力学阻力; γ cp为潜在蒸发条件下的冠层阻力, 随植被类型变化而变化; Δ 为饱和水汽压随温度变化的斜率。

Jackson的CWSI计算公式机理性较强, 物理意义明确, 适于结合冠层温度和微气象条件研究土壤水分状况。许多学者针对CWSI的计算模型和影响因素进行了深入研究, 并广泛应用于土壤水分和干旱监测。张仁华[8]通过重新定义水汽饱和活动面得到了一个新的CWSI微气象参数模型, 并对模型进行了验证和灵敏行分析; 田国良[33]等利用NOAA 资料和气象站资料, 计算了CWSI, 并与气象站观测的5~50 cm土壤水分建立关系, 对河南省进行了土壤水分估算; 申广荣[34]等在地理信息系统的支持下, 通过分析NOAA遥感数据和地面气象站资料进行了黄淮海平原土壤水分监测; 刘安麟[35]等通过能量平衡方程对潜在蒸散的计算进行简化, 得到一种计算CWSI的模型, 根据NOAA数据和气象资料计算了CWSI, 并对陕西地区春季土壤水分状况进行了监测。

作物水分胁迫指数法基于明确的物理公式推导, 意义明确, 机理性强, 精度较高, 是监测土壤水分状况的理想方法。但计算作物水分胁迫指数的方法需要大量参数, 特别是微气象参数在遥感中不易获得; 另外差值过程也会产生一定的误差。

1.6 热惯量法

土壤热惯量是度量土壤阻止其自身温度变化能力大小的物理量, 反映土壤与周围环境能量交换能力的强弱, 表征了土壤的热学特性, 即

P= λρC, (9)

式中: P为土壤热惯量; λ 为土壤热传导率; ρ 为土壤密度; C为比热容。

土壤热惯量随土壤水分增加而增大, 可以通过建立二者之间的关系来监测土壤水分。目前还不能通过遥感手段得到3个参数直接计算真实热惯量, 主要通过求解热传导方程和地表热量平衡方程来实现。Watson[36]于1973年提出了利用地表温度日较差计算热惯量的模型, 随后很多学者发展了热惯量的计算模型; Kahle[37]在热量平衡方程中加入显热和潜热通量, 计算了热惯量; Price[38]在地表能量平衡方程基础上通过引入“ 地表综合参量” 的概念, 提出了一个热惯量计算模型, 即

P=2Q(1-A)/| ω(T午后-T夜间)|-0.9B/ ω, (10)

式中: Q为到达地表的太阳总辐射; A为地表反照率; ω 为地球自转频率; T午后T夜间分别为由卫星数据反演得到的午后和夜间的地表温度; B为一地表综合参量, 需要地面实测数据计算。由于B的计算比较困难, Price[39]简化以上热惯量的计算模型, 提出表观热惯量的概念, 并认为真实热惯量在某些情况下可以用表观热惯量计算。表观热惯量(apparent thermal inertia, ATI)计算公式为

ATI=1 000π (1-A)/(T1330-T0230)。 (11)

式中, T1330T0230为由卫星数据反演得到的13:30时刻和02:30时刻的地表温度。如果只考虑研究区域的反照率和温度变化, 上式可以进一步简化为

ATI=(1-A)/T, (12)

式中△ T为地表温度日较差。

表观热惯量计算模型简单, 所需资料可以由卫星数据提供。但是根据Carlson[40]等人的研究, 当地表蒸发量很大时, 表观热惯量为无效, 在植被覆盖度较大和地表适度较大的情况下, 不能用表观热惯量代替真实热惯量。

热惯量遥感土壤水分的研究主要集中在热惯量计算模型和热惯量与土壤水分关系2个方面, 并取得了很大进步。Xue[41]和Cai[42]等对热惯量模型进行了深入研究, 提出了不同的计算模型; Sobrino[43]等通过3个时刻的遥感信息建立了一种只用卫星数据计算真实热惯量的模型。国内也有很多学者提出了不同的热惯量模式, 张仁华[44]提出了一个考虑显热和潜热通量的热惯量模型, 利用遥感影像中最干点和最湿点订正显热和潜热通量对热惯量计算的影响; 余涛[45]等通过简化地表能量平衡方程, 建立了地表综合参量与热惯量直接的关系, 发展了一种计算真实热惯量的方法; 郭茜[46]等利用表观热惯量对东北地区土壤水分进行了监测; 黄妙芬[47]等提出一个利用地表辐射温度计算表观热惯量的模型; 刘振华[48]等提出了一个利用地表最高温度结合土壤热平衡方程和显热、潜热通量的真实热惯量模型。另外很多学者研究了各种因素对热惯量影响, 如李星敏[49]等分析了植被覆盖对热惯量的影响, 陈怀亮[50]等分析了风速对热惯量的影响, 这些工作提高了热惯量的反演精度。

表观热惯量模型相对简单, 需要参数较少, 有的表观热惯量模型可以只用遥感手段计算得到。对表观热惯量的研究是目前热惯量模式反演土壤水分的重点, 但表观热惯量是表示热惯量相对大小的一个量, 只能表达理论热惯量的一部分。在裸土和低植被覆盖区热惯量模式监测土壤水分效果较好, 但不能用在浓密植被覆盖地区。

1.7 微波法

微波波段(1 mm~1 m)对云层和地表有较强的穿透力, 不依赖太阳光, 可以全天候工作, 具有光学遥感不具备的优势。土壤的介电特性与土壤水分含量有着密切的关系, 水的介电常数大约为80, 而干土仅为3, 它们之间巨大的反差使得土壤的介电常数随土壤湿度的变化而变化, 成为遥感监测土壤水分的重要切入点, 因而微波法具有很好的应用前景。微波遥感法根据工作原理差异可以分为主动微波遥感法和被动微波遥感法。

1.7.1 主动微波遥感法

主动微波遥感是通过遥感器主动发射已知的微波信号, 然后接受所发射的微波信号与地表相互作用后的回波信号, 通过分析2种信号的差异获取地表的后向散射系数(σ ), 而后向散射系数与地表介电常数(ε )直接相关。如前所述, 土壤的地表介电常数与土壤含水量密切相关, 随着含水量的增加, 介电常数几乎呈线性增加; 土壤含水量越高, 介电常数也越高, 雷达信号穿透深度便越小, 后向散射与回波强度越强, 后向则越大。根据以上原理可建立后向散射系数和土壤水分含量的关系用来反演地表土壤水分。通过实验数据的相关分析建立土壤水分与后向散射系数之间的线性回归这种方法简单实用。有的学者还从面散射理论和介电常数模型等理论上建立两者之间的定量关系。Ulady[51]等针对两者的关系进行了定量研究, 提出了一个后向散射模型, 对土壤水分进行了反演, 并估算了植被覆盖对反演精度的影响; Dobson[52]等则认为对于太干或者太湿润的土壤不能用线性模拟; Narayannan等[53]对两者进行了相关的研究, 并发展了一个非线性算法反映两者的关系。国内李杏朝[54]通过同步测量土壤水分、土壤后向散射系数和同步机载雷达图像, 建立了土壤湿度与后向散射系数的线性关系; 田国良[55]等利用合成孔径雷达对冬小麦田进行了土壤水分监测。

1.7.2 被动微波遥感法

被动微波遥感是利用微波辐射计对地表本身发射的微波辐射进行测量来监测土壤水分。在微波波段, 土壤的比辐射率从0.6(湿土)变化到0.9(干土), 通常利用微波辐射计得到土壤的亮度温度, 研究亮度温度与土壤水分的关系, 通过建立亮度温度与土壤水分的关系反演得到土壤水分。被动微波遥感用来监测地表土壤水分的历史更长, 相对来说算法更为成熟。在20世纪70年代NASA就在美国进行了航空平台的被动微波遥感试验, 1987年Schmugge[56]根据试验结果得到了亮度温度与土壤水分的线性关系。许多学者针对地表粗糙度、植被信息等对亮度温度的影响进行了研究, 提高了被动微波遥感法反演土壤水分的精度。常用的反演算法可以分为3类: 基于数理统计的经验算法、基于正向模型的算法和神经网络算法。经验算法是在统计描述和相关分析的基础上直接建立遥感参数与土壤水分之间的统计关系, 应用简单, 对于特定地区适用性较好, 但是缺乏物理基础, 普适性较差; 基于正向模型的算法是充分考虑微波由地表到达传感器的物理过程, 选择一个正向模型, 建立模型参数与土壤水分之间的关系, 其物理意义较明确, 反演精度较高, 但需要的模型参数较多, 反演过程较复杂; 神经网络方法是利用神经网络具有非常强大的模拟任意复杂的非线性关系这个优点来建立遥感数据与土壤水分的关系, 该方法不需要复杂的物理模型和算法, 但受到算法输入端参数特点的限制。

2 方法对比

1)反射率法简单易行, 操作简单快捷, 适宜地形平坦、单一土壤组成较典型的区域。但不同类型土壤发射率的差异对土壤的反射率的影响与土壤水分的作用在同一个量级甚至更大[57], 加之太阳高度、大气条件和地表状况等引起的误差, 使得利用这类方法监测土壤水分精度难以保证。

2)植被指数法需要有较长时间的NDVI积累, 对遥感数据要求较多。

3)在温度-植被指数特征空间中确定干湿边和顶点时只通过温度和NDVI散点图的拟合, 有一定的随意性, 缺乏严格的物理定义, 因此在高植被覆盖区该方法会失效。

4)作物水分胁迫指数法中概念明确、 考虑因素多 , 具有广泛的理论和实用价值 , 但在应用上受到参数的获取性困难和尺度转换问题的制约, 且在作物地表植被较少时, 可能会夸大植被的作用。

5) 微波遥感具有明确的物理基础, 具有全天候、高精度的优势; 且对云层有较强的穿透力, 在地球表面被云层覆盖时, 有其独特的优越性。但是微波遥感数据对地表粗糙度、植被覆盖、地形等参数与对土壤水分同样敏感, 如何确定各种因素对微波遥感信息的影响程度是还需要深入研究的问题。另外, 微波遥感的相对成本过高, 使其在国内的研究和应用都受到了一定的限制。

3 结论与展望

土壤水分的遥感监测方法经过几十年的发展已经取得了巨大的进步, 并且越来越显示出其优越性。然而土壤水分监测是一个非常复杂的问题, 到目前为止还没有一种方法是完美的。笔者认为今后可以从以下方面加强研究:

1)地表温度的反演。目前监测土壤水分的大多数方法中很多都依赖地表温度, 但由于地表温度影响的复杂性, 到目前为止各种遥感地表温度反演方法均存在较大的误差, 即使利用地面仪器仍无法测得精确的地表温度。如何提高地表温度反演的精度一直是遥感学科研究的难点和重点。

2)主动与被动微波遥感相结合。微波遥感具有全天时、全天候和一定的穿透性, 对土壤水分敏感。被动微波遥感方法时间分辨率较高, 但空间分辨率较低; 而主动微波方法空间分辨率较高, 对地表形态比较敏感。充分利用主动与被动微波的优点联合反演土壤水分, 可以达到提高空间分辨率, 提高土壤水分的反演精度的目的。

3)土壤表层水分与深层水分的关系研究。到目前为止, 各种研究方法只能监测到土壤表层的湿度状况, 但是对于农业生产来说, 深层土壤的水分更具有意义, 建立表层土壤水分与深层土壤水分的模型可以促进遥感监测土壤水分方法的发展。

4)模型的优化与改进。遥感信息与土壤水分的关系受到辐射传输等多种因素的影响, 两者之间的关系通常是复杂的非线性关系, 需要从物理机制等方面入手深入研究, 以建立遥感信息与土壤水分之间更加完美的科学模型。

The authors have declared that no competing interests exist.

参考文献
[1] 邓辉, 周清波. 土壤水分遥感监测方法进展[J]. 中国农业资源与区划, 2004, 25(3): 46-49.
Deng H, Zhou Q B. The progress in remote sensing method for monitoring soil moisture[J]. Journal of China Agricultural Resources and Regional Planning, 2004, 25(3): 46-49. [本文引用:1]
[2] Bowers S A, Hunks R J. Reflection of radiant energy from soils[J]. Soil Science, 1965, 100(2): 130-138. [本文引用:1] [JCR: 1.051]
[3] Curran P J. The use of polarized panchromatic and false color infrared film in the monitoring of soil surface moisture[J]. Remote Sensing of Environment, 1979, 8(3): 249-266. [本文引用:1] [JCR: 5.103]
[4] Robinove C J, Chavez P S, Gehring D, et al. Arid land monitoring using Land sat albedo difference images[J]. Remote Sensing of Environment, 1981, 11(2): 133-156. [本文引用:1] [JCR: 5.103]
[5] Henricksen B L. Reflections on drought: Ethiopia 1983—1984[J]. International Journal of Remote Sensing, 1986, 7(11): 1447-1451. [本文引用:1]
[6] Everitt J H, Escobar D E, Alaniz M A, et al. Using multispectral video imagery for detecting soil surface conditions[J]. Photogrammetric Engineering and Remote Sensing, 1989, 55(4): 467-471. [本文引用:1] [JCR: 1.802]
[7] 刘培君, 张琳. 卫星遥感估测土壤水分的一种方法[J]. 遥感学报, 1997, 1(2): 135-138, 81.
Liu P J, Zhang L. A method for monitoring soil water contents using satellite remote sensing[J]. Journal of Remote Sensing, 1997, 1(2): 135-138, 81. [本文引用:1] [CJCR: 0.992]
[8] 张仁华. 对于定量热红外遥感的一些思考[J]. 国土资源遥感, 1999, 11(1): 1-6.
Zhang R H. Some thinkings on quantitative thermal infrared remote sensing[J]. Remote Sensing for Land and Resources, 1999, 11(1): 1-6. [本文引用:3]
[9] 詹志明, 秦其明, 阿布都瓦斯提·吾拉本, 等. 基于NIR-Red光谱特征空间的土壤水分监测新方法[J]. 中国科学D辑: 地球科学, 2006, 36(11): 1020-1026.
Zhan Z M, Qin Q M, Ghulan A, et al. Based on NIR-Red spectral space based new method for soil moisture monitoring[J]. Science in China Series D: Earth Sciences, 2006, 36(11): 1020-1026. [本文引用:1]
[10] 陈维英, 肖乾广, 盛永伟. 距平植被指数在1992年特大干旱监测中的应用[J]. 环境遥感, 1994, 9(2): 106-112.
Chen W Y, Xiao Q G, Sheng Y W. Application of the anomaly vegetation index to monitoring heavy drought in 1992[J]. Remote Sensing of Environmental China, 1994, 9(2): 106-112. [本文引用:1]
[11] 王鹏新, 龚健雅, 李小文. 条件植被温度指数及其在干旱监测中的应用[J]. 武汉大学学报(信息科学版), 2001, 26(5): 412-418.
Wang P X, Gong J Y, Li X W. Vegetation-temperature condition index and its application for drought monitoring[J]. Editoral Board of Geomatics and Information Science of Wuhan University, 2001, 26(5): 412-418. [本文引用:2]
[12] 蔡斌, 陆文杰, 郑新江. 气象卫星条件植被指数监测土壤状况[J]. 国土资源遥感, 1995, 7(4): 45-50, 66.
Cai B, Lu W J, Zheng X J. Using meterorological satellite’s VIC monitor soil state[J]. Remote Sensing for Land and Resources, 1995, 7(4): 45-50, 66. [本文引用:1]
[13] Mcvicar T R, Jupp D L B, Yang X, et al. Linking regional water balance models with remote sensing[C]//Proceedings of the 13th Asian Conference on Remote Sensing. Ulaanbaatar, Mongolia: 1992. B6. 1-B6. 6. [本文引用:1]
[14] Kogan F N. Application of vegetation index and brightness temperature for drought detection[J]. Advances in Space Research, 1995, 15(11): 91-100 [本文引用:1] [JCR: 1.183]
[15] 刘志明. 利用气象卫星信息遥感土壤水分的探讨[J]. 遥感信息, 1992(1): 21-23.
Liu Z M. Discussed on using meteorological satellite remote sensing of soil moisture information[J]. Remote Sensing Information, 1992(1): 21-23. [本文引用:1] [CJCR: 0.7]
[16] 乔平林, 张继贤, 燕琴, . 利用TM6进行土壤水分的监测研究[J]. 测绘通报, 2003(7): 14-15, 18.
Qiao P L, Zhang J X, Yan Q, et al. A study of soil water contents monitoring using TM6[J]. Bulletin of Surveying and Mapping, 2003(7): 14-15, 18. [本文引用:1] [CJCR: 0.661]
[17] 张树誉, 赵杰明, 袁亚社, . NOAA/AVHRR资料在陕西省干旱动态监测中的应用[J]. 中国农业气象, 1998, 19(5): 26-29.
Zhang S Y, Zhao J M, Yuan Y S, et al. Application of NOAA/ AVHRR data to dynamic monitoring of drought in Shaanxi Province[J]. Chinese Journal of Agrow eteorology, 1998, 19(5): 26-29. [本文引用:1]
[18] 毛学森, 张永强, 沈彦俊. 水分胁迫对冬小麦植被指数NDVI影响及其动态变化特征[J]. 干旱地区农业研究, 2002, 20(1): 69-71.
Mao X S, Zhang Y Q, Shen Y J. Effects of water stress in different stages on winter wheat vegetation index and NDVI’s dynamic[J]. Agricultural Research in the Arid Areas, 2002, 20(1): 69-71. [本文引用:1] [CJCR: 0.916]
[19] Goetz S J. Multi-sensor analysis of NDVI, surface temperature and biophysical variables at a mixed grassland site[J]. International Journal of Remote Sensing, 1997, 18(1): 71-94. [本文引用:1]
[20] Price J C. Using spatial context in satellite data to infer regional scale evapotranspiration[J]. IEEE Transactions on Geoscience and Remote Sensing, 1990, 28(5): 940-948. [本文引用:1] [JCR: 3.467]
[21] Nemani R, Pierce L, Running S. Developing satellite-derived estimates of surface moisture status[J]. Journal of Applied Meteorology, 1993, 32(3): 548-557. [本文引用:1] [JCR: 1.702]
[22] Moron M S, Clarke T R, Inoue Y, et al. Estimating crop water deficit using the relation between surface air temperature and spectral vegetation index[J]. Remote Sensing of Environment, 1994, 49(3): 246-263. [本文引用:1] [JCR: 5.103]
[23] Sand holt I, Rasmussen K, Andersen J. A simple interpretation of the surface temperature-vegetation index space for assessment of surface moisture status[J]. Remote Sensing of Environment, 2002, 79(2-3): 213224. [本文引用:1] [JCR: 5.103]
[24] Lambin E F, Ehrlich D. The surface temperature-vegetation index for land cover and land -cover change analysis[J]. International Journal of Remote Sensing, 1996, 17: 463-487. [本文引用:1]
[25] 韩丽娟, 王鹏新, 王锦地, . 植被指数-地表温度构成的特征空间研究[J]. 中国科学D辑: 地球科学, 2005, 35(4): 371-377.
Han L J, Wang P X, Wang J D, et al. Surface temperature vegetation index-congstitute the characteristics of space research[J]. Science in China Series D: Earth Sciences, 2005, 35(4): 371-377. [本文引用:1]
[26] 刘良云, 张兵, 郑兰芬, . 利用温度和植被指数进行地物分类和土壤水分反演[J]. 红外与毫米波学报, 2002, 21(4): 269-273.
Liu L Y, Zhang B, Zheng L F, et al. Target classification and soil water content regression using land surface temperature and vegetation index[J]. Journal of Infrared and Millimeter Waves, 2002, 21(4): 269-273. [本文引用:1] [JCR: 0.292] [CJCR: 1.14]
[27] 姚春生, 张增祥, 汪潇. 使用温度植被干旱指数法(TVDI)反演新疆土壤湿度[J]. 遥感技术与应用, 2004, 19(6): 473-478.
Yao C S, Zhang Z X, Wang X. Evaluating soil moisture status in Xinjiang using the temperature vegetation dryness index[J]. Remote Sensing Technology and Application, 2004, 19(6): 473-478. [本文引用:1] [CJCR: 0.749]
[28] Carlson T N, Gillies R R, Perry E M. A method to make use of thermal infrared temperature and NDVI measurements to infer soil water content and fractional vegetation cover[J]. Remote Sensing Reviews, 1994, 9(1-2): 161-173. [本文引用:1]
[29] 刘丽, 周颖, 杨凤, . 用遥感植被供水指数监测贵州干旱[J]. 贵州气象, 1998, 22(6): 17-21.
Liu L, Zhou Y, Yang F, et al. With remote sensing of vegetation water supply index monitoring drought in Guizhou[J]. Journal of Guizhou Meteorology, 1998, 22(6): 17-21. [本文引用:1] [CJCR: 0.509]
[30] 莫伟华, 王阵会, 孙涵, . 基于植被供水指数的农田干旱遥感监测研究[J]. 南京气象学院学报, 2006, 29(3): 396-401.
Mo W H, Wang Z H, Sun H, et al. Remote sensing monitoring of farmland drought based on vegetation supply water index[J]. Journal of Nanjing Institute of Meteorology, 2006, 29(3): 396-401. [本文引用:1] [CJCR: 1.005]
[31] Idso S B, Jackson R D, Pinter P J, et al. Normalizing the stress degree-day parameter for environmental variability[J]. Agricultural Meteorology, 1981, 24: 45-55. [本文引用:1]
[32] Jackson R D, Jdso S B, Reginato R J, et al. Canopy temperature as a crop water stress indicator[J]. Water Resources Research, 1981, 17(4): 1133-1138 [本文引用:1] [JCR: 3.149]
[33] 田国良, 杨希华, 郑柯. 冬小麦旱情遥感监测模型研究[J]. 环境遥感, 1992, 7(2): 83-89.
Tian G L, Yang X H, Zheng K. Remote sensing model for wheat drought monitoring[J]. Remote sensing of Environmental, 1992, 7(2): 83-89. [本文引用:1]
[34] 申广荣, 田国良. 基于GIS的黄淮海平原旱灾遥感监测研究—作物缺水指数模型的实现[J]. 生态学报, 2000, 20(2): 224-228.
Shen G R, Tian G L. Remote sensing monitoring of drought in Huanghe, Huaihe and Haihe plain based on GIS the calculation of crop water stress index model[J]. Acta Ecologica Sinica, 2000, 20(2): 224-228. [本文引用:1]
[35] 刘安麟, 李星敏, 何延波, . 作物缺水指数法的简化及在干旱遥感监测中的应用[J]. 应用生态学报, 2004, 15(2): 210-214.
Liu A L, Li X M, He Y D, et al. Simplification of crop shortage water index and its application in drought remote sensing monitoring[J]. Chinese Journal of Applied Ecology, 2004, 15(2): 210-214. [本文引用:1] [CJCR: 1.742]
[36] Watson K, Rowen L C, Offield T W. Application of thermal modeling in the geologic interpretation of IR images[J]. Remote Sensing of Environment, 1971, 3: 2017-2041. [本文引用:1] [JCR: 5.103]
[37] Kahle A B. A simple thermal model of the earth’s surface for geologic mapping by remote sensing[J]. Journal of Geophysical Research, 1977, 82(11): 1673-1680. [本文引用:1]
[38] Price J C. Thermal inertia mapping: A new view of the earth[J]. Journal of Geophysical Research, 1977, 82(18): 2582-2590. [本文引用:1]
[39] Price J C. On the analysis of thermal infrared imagery: The limited utility of apparent thermal inertia[J]. Remote Sensing of Environment, 1985, 18(1): 59-73. [本文引用:1] [JCR: 5.103]
[40] Carlson T N, Boland F F. Analysis of urban-rural canopy using a surface heat flux temperature model[J]. Journal of Applied Meteorology, 1978, 17(7): 998-1013. [本文引用:1] [JCR: 1.702]
[41] Xue Y, Cracknell A P. Advanced thermal inertia modeling[J]. International Journal of Remote Sensing, 1995, 16(3): 431-446. [本文引用:1]
[42] Cai G, Xue Y, Hu Y, et al. Soil moisture retrieval from MODIS data in Northern China Plain using thermal inertia model[J]. International Journal of Remote Sensing, 2007, 28(16): 3567-3581. [本文引用:1]
[43] Sobrino J A, El Kharraz M H. Combining afternoon and morning NOAA satellites for thermal inertia estimation: 1. Algorithm and its testing with Hydrologic Atmospheric Pilot Experiment-Sahel data[J]. Journal of Geophysical Research, 1999, 104: 9445-9453. [本文引用:1]
[44] 张仁华. 改进的热惯量模式及遥感土壤水分[J]. 地理研究, 1990, 9(2): 101-112.
Zhang R H. Improved thermal inertia model and remote sensing of soil moisture[J]. Geographical Research, 1990, 9(2): 101-112. [本文引用:1]
[45] 余涛, 田国良. 热惯量法在监测土壤表层水分变化中的研究[J]. 遥感学报, 1997, 1(1): 24-31.
Yu T, Tian G L. The application of thermal inertia method the monitoring of soil moisture of north China plain based on NOAA-AVHRR data[J]. Journal of Remote Sensing, 1997, 1(1): 24-31. [本文引用:1] [CJCR: 0.992]
[46] 郭茜, 李国春. 用表观热惯量法计算土壤含水量探讨[J]. 中国农业气象, 2005, 26(4): 215-219.
Guo Q, Li G C. Monitoring of soil moisture by apparent thermal inertia method[J]. Chinese Journal of Agrometeorology, 2005, 26(4): 215-219. [本文引用:1] [CJCR: 1.474]
[47] 黄妙芬, 邢旭峰, 刘素红, . 运用热惯量估算大气下行长波辐射的遥感方法研究——以Land sat/TM为例[J]. 资源科学, 2006, 28(3): 37-43.
Huang M F, Xing X F, Liu S H, et al. Estimation of downward atmospheric long-wave radiation using thermal inertia based on remote sensing[J]. Resources Science, 2006, 28(3): 37-43. [本文引用:1] [CJCR: 2.123]
[48] 刘振华, 赵英时. 一种改进的遥感热惯量模型初探[J]. 中国科学院研究生学报, 2005, 22(3): 380-385.
Liu Z H, Zhao Y S. An improved thermal inertia model[J]. Journal of the Graduate School of the Chinese Academy of Science, 2005, 22(3): 380-385. [本文引用:1] [CJCR: 0.4174]
[49] 李星敏, 刘安麟, 张树誉, . 热惯量法在干旱遥感监测中的应用研究[J]. 干旱地区农业研究, 2005, 23(1): 54-59.
Li X M, Liu A L, Zhang S Y, et al. Use of thermal inertia approach in the monitoring of drought by remote sensing[J]. Agricultural Research in the Arid Areas, 2005, 23(1): 54-59. [本文引用:1] [CJCR: 0.916]
[50] 陈怀亮, 冯定原, 邹春辉. 麦田土壤水分NOAA/AVHRR遥感监测方法研究[J]. 遥感技术与应用, 1998, 13(4): 27-35.
Chen H L, Feng D Y, Zou C H. The study of method on monitoring soil moisture in wheat field by NOVV/AVHRR data[J]. Remote Sensing Technology and Application, 1998, 13(4): 27-35. [本文引用:1] [CJCR: 0.749]
[51] Ulaby F T, Aslam A, Dobson M C. Effects of vegetation cover on the Radar sensitivity to soil moisture[J]. IEEE Transactions on Geoscience and Remote Sensing, 1982, 20(4): 476-481. [本文引用:1] [JCR: 3.467]
[52] Ulaby F T, Batlivala P P, Dobson M C. Microwave backscatter dependence on surface roughness, soil moisture, and soil texture[J]. IEEE Transactions on Geoscience and Remote Sensing, 1978, 16(4): 286-295. [本文引用:1] [JCR: 3.467]
[53] Narayanan R M, Horner J R, St Germain K M. Simulation study of a robust algorithm for soil moisture and surface roughness estimation using L-band Radar back scatter[J]. Geocarto International, 1999, 14(1): 6-13. [本文引用:1] [JCR: 0.575]
[54] 李杏朝. 微波遥感监测土壤水分的研究初探[J]. 遥感技术与应用, 1995, 10(4): 1-8.
Li X C. A study on soil moisture monitoring by using microwave remote sensing[J]. Remote Sensing Technology and Application, 1995, 10(4): 1-8. [本文引用:1] [CJCR: 0.749]
[55] 田国良, 李长乐, 杨习荣, . 机载合成孔径雷达图像监测土壤水分的初步分析[C]//田国良. 黄河流域典型地区遥感动态研究. 北京: 科学出版社, 1990: 102-110.
Tian G L, Li C L, Yang X R, et al. Monitoring soil moisture using image of airborne synthetic aperture Radar: A case study[C]//Tian G L. Dynamic Research by Remote Sensing on Typical Region in the Yellow River Vally. Beijing: Science Press, 1990: 102-110. [本文引用:1]
[56] Schmugge T J, O’Neill P E, Wang J R. Passive microwave soil moisture research[J]. IEEE Transactions on Geoscience and Remote Sensing, 1986, 24(1): 12-20. [本文引用:1] [JCR: 3.467]
[57] 虞献平, 贺红仕. 生态与环境遥感研究[M]. 北京: 科学出版社, 1990.
Yu X P, He H S. Study on ecological and environmental and remote sensing[M]. Beijing: Science Press, 1990. [本文引用:1]