西北干旱区河西走廊荒漠绿洲土地覆盖类型与蒸散的关系研究——基于Landsat 8和ZY3数据融合

焦丹丹,吉喜斌,金博文,赵丽雯,张靖琳,郭 飞

1 中国科学院西北生态环境资源研究院临泽内陆河流域研究站,中国科学院生态水文与流域科学重点实验室, 兰州 730000 2 中国科学院大学, 北京 100049

蒸散是地气相互作用[1]、全球气候与水循环及农业灌溉用水[2]等领域研究的重要内容,在调节局地小气候乃至区域气候方面,蒸散可作为潜热减少净辐射向感热的转换,对改善大气环流有重要影响[3- 4]；另一方面,蒸散作为地气水热交换的主要组分,在将地表水汽带到大气中的同时消耗了大量地表热量,影响了不同下垫面地表各分量和降水的重新分配,在土壤-植被-大气系统中也发挥了重要作用[5]。Rosenberg等[6]指出降落在地面的降水有70%通过蒸散作用回到大气中,在干旱区甚至可达90%,因此提高水分利用效率对于干旱区而言至关重要,而任何提高水分利用效率的尝试都应基于对每日蒸散量的精确估计[7],Ramos等[8]曾使用遥感数据来评估灌溉用水量,依靠SEBAL模型计算了实际蒸散值并通过水分利用率指数计算了净灌溉量,发现此法是可行的；Yang等[9]发现黄河引水的减少将直接导致灌区内水量减少,为了提高水资源利用率,需要根据变化的用水策略来对不同地表类型的用水量进行分配,因此获取蒸散值是必不可少的。以上研究均涉及蒸散变化趋势和灌溉水量分配情况,表明作物用水量的计算需要蒸散量来辅助,而要达到最大水分利用效率和最高产量,灌水的时间和总量都要进行合理安排[10],由此可见蒸散的研究已经被广泛应用在农田作物生长方面[11-12],它不仅可以反映作物耗水的空间分布状况,据此来做灌溉计划,还能减少作物因水分缺失或过度灌溉而导致减产现象的发生,因此准确估算蒸散对干旱区气候动态变化及维持水资源供需平衡至关重要。

作为水文循环中最难准确估算的分量之一,自1802年Dalton提出计算蒸发的公式以来,出现了许多蒸散量的计算公式,其中传统估算法有蒸渗仪法[13]、波文比-能量平衡法[14]、涡度相关法[15]以及公式计算法[16-23],其中适用性最好的当属基于能量平衡和空气动力学的Penman法和FAO- 56法。传统估算蒸散的方法虽然能够获取相对准确的均匀下垫面的蒸散量,但大多局限于一点或同一种作物的均匀下垫面来开展,难以反映蒸散的空间异质性。因遥感技术能够提供大范围地表特征参数,它的应用使得区域尺度蒸散量的估算成为可能,是定量估算区域尺度蒸散量最有效的途径,干旱区下垫面几何结构和物理性质的水平不均匀性使得该方法在干旱区蒸散研究方面得到了广泛应用。目前应用较多的有地表能量平衡指数SEBI模型[24]、陆面能量平衡算法SEBAL模型[25]、简化的地表能量平衡指数S-SEBI模型[26]、地表能量平衡系统SEBS模型[27]和METRIC模型[28]等。以上模型均为单层模型,因其使用效果较好,近年来在国内外被广泛应用于估算不同气候条件下的地表蒸散。Elhag等[7]利用SEBS模型和MERIS、AATSR资料估算了尼罗河三角洲地区的日蒸散和蒸发比,并与均匀分布于研究区内的92个地面实测点数据进行对比验证,结果具有很强的一致性；Folhes等[29]利用METRIC模型和Landsat遥感数据估算了巴西东北部半干旱地区的蒸散量,将计算值与微气象塔实测蒸散值进行对比,研究结果为灌溉农业用水量提供了数据支撑；何延波等[30]在对SEBS模型参数进行订正的基础上,利用MODIS遥感数据对黄渤海地区地表通量进行了估算,并分别与地面观测点和LAS数据进行对比验证,结果表明精度可靠且能满足日蒸散量计算的要求。从众多研究可以看出,利用遥感计算蒸散需要反演的其中一个重要参数为地表温度,将地表温度与其他参数相结合能得到较为精确的近地表层湍流通量值和下垫面的干湿状况[31-32],对于计算区域蒸散意义重大,而温度的反演需要用到热红外波段数据,但高分辨率的热红外数据难以获取,因此尽管对区域范围内蒸散发估算的研究较多,但大部分为中低分辨率尺度获取蒸散值[33-34],为提高计算精度,有必要在该研究区开展结合高空间分辨率数据利用多尺度空间融合获取日蒸散量的研究,为小区域范围内水资源合理利用与精准分配提供借鉴。

本研究区位于西北干旱区河西走廊中段临泽绿洲北部,地表覆盖包括绿洲内部农田、绿洲-荒漠过渡带、人工-天然植被防风固沙林和荒漠。利用Landsat 8遥感影像的反演参数及气象数据估算研究区的日蒸散量,依靠临泽内陆河流域研究站的通量塔涡动相关数据和气象数据来验证遥感反演的精度,将不同季节的蒸散值与通量塔通量数据进行对比验证,分析SEBS模型在该地区的适用性,同时探讨了不同土地覆盖类型对蒸散量空间分布的影响；并利用资源三号(ZY3)影像的高分辨率植被信息来获取空间增强后的地表温度值,以提高蒸散量的估算精度,通过遥感手段来估算蒸散量为陆面过程的研究和水资源管理提供了关键参量,对于农业生态和水文学的发展意义重大[35]。

1 资料与方法

1.1 研究区概况

该荒漠绿洲位于甘肃省张掖市临泽县北部(100°06′04″E—100°09′53″E,39°19′07″N—39°24′40″N),地处河西走廊中段,属黑河流域中游,最高海拔1481 m,最低1273 m,根据中国科学院临泽内陆河流域研究站(100°07′48″E,39°20′50″N,海拔1384 m)气象站的年观测数据[36],该地区年均气温9.2℃,多年平均降水量124 mm,有77%集中在7、8月份,年均无霜期195 d,属温带干旱气候。研究区内土地覆盖类型复杂,东部和西北部多裸地,中部和南部多农田及农田外围防护林地,主要种植杨树,裸地与农田的过渡带长有梭梭、泡泡刺、沙拐枣等防风固沙灌木,在雨水较多的年份还会生长刺蓬、雾冰藜、白茎盐生草等一年生草本植物,西南有黑河流经,为典型的荒漠绿洲农田,由于农田比例较高,区内受人类活动影响较大。该研究区主要种植作物有制种玉米、大田玉米,作物生长以灌溉为主,灌区内密布各级水渠,在生长季内大部分区域通过水渠引平川水库内水进行灌溉,也有部分抽取地下水灌溉,由于该研究区地处西北干旱区,水资源短缺,蒸发旺盛,加之人类对水资源的不合理利用,导致该地区生态环境十分脆弱,沙尘暴等恶劣天气时有发生,因此合理利用水资源以提高利用率对于生态环境保护和该地区农业的发展具有重要意义。

1.2 研究方法

本文采用SEBS模型[27],该模型是一种利用遥感数据和气象数据来估算地表通量的模型,模型的计算需输入两套参数：第一套为遥感反演数据,其中包括空间分布的地表反照率α、植被归一化指数NDVI、植被覆盖度Pv、叶面积指数LAI、地表比辐射率ε、地表温度Ts、高程DEM等,第二套为研究区内卫星过境时刻的气象数据,包括气温、风速、湿度、气压、日照时数等,在输入ILWIS软件的SEBS模块时,将第一套参数中的遥感反演数据统一为30 m分辨率的TIFF格式以计算日蒸散量。它的理论基础是地表能量平衡方程：

Rn=G0+H+λE

(1)

式中,Rn为净辐射量(W/m2),G0为土壤热通量(W/m2),H为感热通量(W/m2),λE为潜热通量(W/m2),其中λ为水的蒸发潜热,取值2.49×106J/kg,E为实际蒸散量(kg m-2s-1)。

1.2.1净辐射量

(2)

1.2.2土壤热通量

G0=Rn[Γc+(1-Pv)(Γs-Γc)]

(3)

式中,Γc=0.05,为全植被覆盖条件下土壤热通量与净辐射的比率；Γs=0.315,为裸土条件下土壤热通量与净辐射的比率。Pv为植被覆盖度[39-40],由归一化植被指数NDVI计算得到：

Pv=[(NDVI-NDVIs)/(NDVIv-NDVIs)]2

(4)

式中,NDVIv、NDVIs分别为植被和裸土的NDVI值,分别取研究区内累积百分比98%和5%的值。

根据下式中叶面积指数LAI与植被覆盖度之间的关系,可以得到LAI值[41]：

LAI=[ln(1-Pv)]/(-0.5)

(5)

1.2.3感热通量

(6)

(7)

(8)

式中,u为风速(m/s),u*为摩擦风速(m/s),k为von Karman常数,取值0.4,z为参考高度(m),d0为零平面位移(m),z0m、z0h分别为动量、热量交换粗糙度(m),ψm、ψh分别为动力交换和热力交换的MOS稳定度订正函数,L为奥布霍夫长度,定义为式(8),θ0、θa分别为观测高度和参考高度的位温(K),ρ为空气密度(kg/m3),Cp为空气热容(J kg-1k-1),θv为近地表虚位温(K),g为重力加速度(m/s2)。其中,摩擦风速、奥布霍夫长度和感热通量利用MOS相似性理论通过Broyden非线性方程组进行求解,其余所需参数均可通过气象数据结合遥感影像信息获取。

1.2.4日蒸散量

由于蒸发比在一日之内可被认为是常数,利用蒸发比不变法可进行日时间尺度扩展[42],因此由卫星过境时刻的瞬时蒸散值可得到日蒸散量：

(9)

在极端干旱状态下,由于土壤水分的缺失潜热通量约为零,而此时的感热通量达最大值：Hdry=Rn-G0,式中：Hdry为干燥地表环境下的感热通量(W/m2)；在土壤供水充足的极端湿润环境下,蒸散发旺盛且达到潜在蒸发速率,此时的感热通量取最小值：Hwet=Rn-G0-λEwet,式中：Hwet为湿润地表环境下的感热通量(W/m2),λEwet为湿润地表环境下的潜热通量(W/m2)。

根据以上两种极端状态就可计算相对蒸发比Λr,在该模型中定义Λr为实际蒸发速率与潜在蒸发速率(即湿限蒸发速率)之比：

(10)

由此可知瞬时蒸散量可用下式计算：

(11)

同时可以得到蒸发比Λ,定义为实际蒸散与潜热通量与感热通量之和的比值,即可带入式(9)计算日蒸散量：

(12)

1.3 数据来源与处理

根据Landsat 8影像条带号及本研究区所在地理位置,选择133行、33列遥感数据,下载地址为https://earthexplorer.usgs.gov/,该影像全色波段空间分辨率为15 m,基本分辨率30 m,还包括了两个100 m空间分辨率的热红外波段,所下载热红外数据已经过重采样将分辨率转为30 m。本文选取了研究区2014年8景质量较好的Landsat 8遥感影像(1月13日、2月14日、3月18日、6月6日、7月24日、8月9日、9月26日、11月13日),这几日均为晴朗天气,研究区云量均在8%以下,并用ENVI 5.3软件进行辐射定标、大气校正、裁剪等处理,反演出SEBS模型需要输入的地表参数。同时在该网站下载ASTER DEM高程数据,投影为UTM/WGS84,空间分辨率为1 ″(约30 m),按照研究区边界进行拼接裁剪。在2014年6—8月份覆盖该研究区的众多高分辨率遥感影像中,资源三号卫星在2014年8月27日获取的该研究区影像晴朗无云,因此选购该景影像与2014年8月9日的Landsat 8影像进行融合,因时间距离较近,可将植被信息看作同一时相,其全色数据空间分辨率为2.1 m,多光谱数据分辨率为5.8 m。本文只选了生长季的高分辨率影像来进行植被空间信息的增强,而未选择冬季月份的,是因为该研究区冬季降水稀少,几乎无植被,蒸散量较小,即便使用高分辨率影像进行空间信息增强,变化也不甚明显。

图1 研究区土地覆盖类型分类Fig.1 Classification of land cover in study area

气象数据来自临泽内陆河流域研究站,同时收集位于该研究区内的绿洲内部农田观测场和绿洲-荒漠过渡带观测场两个不同下垫面(图1)的通量塔数据对相应位置像元内遥感反演的日蒸散值进行地面验证。

1.3.1土地覆盖类型分类

土地覆盖类型的特点影响着蒸散量的空间分布,为进一步分析不同土地覆盖类型对蒸散发的影响,运用监督分类里的支持向量机分类法将地表分为裸地、防护林、梭梭群落、白刺群落、农田、水体、其他共计7类,由于荒漠内部在几米范围内就会有多种植被覆盖,为提高分类精度,选用ZY3数据进行分类,并用混淆矩阵法检验分类精度,总体分类精度为88.96%,Kappa系数为0.83,分类精度较高,图1所示为地物类型图,以6种地表类型(在统计分析时将梭梭群落和白刺群落合并为绿洲外围防风固沙灌木丛)为掩膜图像,借助ArcGIS软件中的区域分析功能,计算各类土地覆盖在不同时间的日蒸散值并分析。

1.3.2地表温度求取

基于Kustas[43]和Anderson[44]提出的植被指数与地表温度之间的最小二乘法拟合关系,可以利用高分辨率的植被信息获取空间增强后的地表温度,该方法使用二次回归分析式而非线性公式求解回归系数是因为线性公式只能考虑到平均状态,会忽略土壤湿度对模型的影响。本研究利用空间分辨率为2.1 m的ZY3影像反演植被指数NDVI,并结合Landsat 8影像模拟高空间分辨率的地表温度信息,步骤如下：

(1)将2.1 m空间分辨率的NDVI值重采样至30 m,在不同地物类型上分别选取纯像元的NDVI30值(ZY3提供)和T30值(Landsat 8提供),计算系数a、b、c,王万同等[45]构建了一种纯像元选取需遵循的原则：区域相似度判定函数,可用来进行纯像元的选取。

T30=a+bNDVI30+cNDVI302

(13)

(14)

1.3.3蒸散量测定

蒸散量利用位于绿洲内部农田和绿洲-荒漠过渡带的开路式涡动相关系统测定的潜热通量换算得到,其中位于绿洲内部农田的开路式涡动相关系统由三维超声风速仪(HS- 50,Gill Solent Instruments,英国)和红外CO2/H2O气体分析仪(LI- 7500,LI-COR inc.,美国)组成,观测高度距地面4.5 m,风浪区范围均为灌溉玉米农田；绿洲-荒漠过渡带开路式涡动相关系统由三维超声风速仪(WindMaster Pro, Gill Solent Instruments,英国)和红外CO2/H2O气体分析仪(LI- 7500A,LI-COR inc.,美国)组成,观测高度距地面8.0 m,风浪区范围是由梭梭、沙拐枣和白刺等组成的稀疏灌木林地。两套开路式涡动相关观测系统的红外CO2/H2O气体分析仪与垂直方向夹角为10°,与三维超声风速仪的感应面相距均为20 cm。两个通量塔周围地势相对平坦。

2 结果

2.1 SEBS模型验证

绿洲农田观测场下垫面主要种植玉米,植被均一,荒漠-绿洲过渡带观测场下垫面为梭梭群落、白刺群落、裸地等,植被类型较为复杂,为保证数据的准确性,在SEBS模型计算的蒸散结果中提取两个通量塔所在经纬度的像元值与通量塔数据分别进行对比,以模型计算值为自变量,实测日蒸散值为因变量,用回归方程对二者进行拟合,以决定系数R2来表示回归方程拟合的好坏,同时计算了模型计算值与实测日蒸散值之间的均方根误差(RMSE)、平均绝对误差(MAE),二者绝对值越小,表明两种方法之间的差异越小。由图2明显可以看出：模型估算值不管是在绿洲农田还是在绿洲-荒漠过渡带都存在高估现象,且在蒸散量越小时高估现象越严重。农田区模型计算值与实测值回归方程斜率为1.07(R2=0.96,P<0.001),相比回归方程斜率为0.80(R2=0.72,P<0.001)的绿洲-荒漠过渡带来看,变化趋势较为一致且方程拟合效果较好,且农田区随着蒸散量的增加,模型值与实测值的拟合程度越好,说明用SEBS模型估算该研究区蒸散值是可行的。单独看SEBS模型计算值与绿洲农田蒸散值之间的差异(表1),RMSE、MAE分别为0.84 mm/d、0.56 mm/d,而模型计算值与绿洲-荒漠过渡带蒸散值之间的RMSE、MAE分别为1.33 mm/d、1.12 mm/d,均比与绿洲农田之间的差异要大,说明在下垫面为均一农田时,SEBS模型的估算结果更为准确。

2.2 绿洲-荒漠过渡带蒸散季节变化

SEBS模型计算的研究区日蒸散量分布如图3所示,四季可以明显看出该地区蒸散量的季节变化规律：夏季>春季>秋季>冬季。空间上整体变化特征主要受土地覆盖类型的影响,南部多农田、林地,靠近水域(河流：黑河；水库：平川水库),多植被覆盖且供水充足,蒸散量较大；北部和东部多裸地,水分不足因此蒸散值较小。

图2 实测值和地表能量平衡系统(SEBS)模型计算值之间的相关性Fig.2 Correlation between measured values and ETdaily (daily evapotranspiration) estimated by surface energy balance system (SEBS) model

Table 1 Comparison between measured values and ETdaily(daily evapotranspiration)estimated by surface energy balance system(SEBS)model

时间TimeSEBS计算值Simulated values by SEBS/(mm/d)绿洲农田值Values for oasis/(mm/d)SEBS计算值Simulated values by SEBS/(mm/d)绿洲-荒漠过渡带值Values for oasis-desert ecotone/(mm/d)1-130.440.090.960.022-140.970.101.610.093-182.010.272.250.066-62.722.741.830.777-247.166.683.622.208-95.665.382.041.299-261.760.881.580.4011-131.020.221.110.07RMSE:均方根误差Root mean squared error0.841.33MAE:平均绝对误差Mean absolute error0.561.12

由图3可知研究区的时间变化特征：1月13日、2月14日表征冬季蒸散情况,气温低且地表裸露,水库内水结冰,蒸发微弱,大部分蒸散量均在1.5 mm/d以下；比较这两日蒸散值可知,2月14日蒸散量在整个研究区内相比1月13日有小幅上升,这与温度值升高有关。

图3 SEBS模型计算日蒸散量对比Fig.3 ETdaily estimated by SEBS model for various days

3月18日代表初春,农田地区作物还未播种,但因温度升高,蒸散值也相应增大,均在2.0 mm/d左右,水库与河流解冻,蒸发加强。6月6日为春末夏初时节,农田作物已进入生长阶段,该地区为引水灌溉,且此时作物已经过第一轮浇水及施肥,生长旺盛,蒸散量明显增加。从图3可以看到北部农田蒸散值比南部高,这是由于引水灌溉并不能在一天内对研究区内所有农田进行浇水,这充分说明了水分对于蒸散的影响是极其显著的,同时水库、河流及湖面(中部水域为锁龙潭)的蒸散值大幅上升,可达7.0 mm/d。

7月24日蒸散结果代表仲夏时节,此时为该地区雨季,是荒漠地区植物的生长季,因此荒漠地区蒸散值增加明显,约为3.0 mm/d；农田作物此时进入生长关键期,需大量引水灌溉,且已经过第二轮施肥,蒸散量达全年最高,均在6.0 mm/d以上,同时蒸散差异也达最大。8月9日整个研究区蒸散值较7月24日有所回落,但仍处于较高水平。

9月26日蒸散结果对应秋季,此时农作物生长结束,进入收割季节,该地区种植的作物大部分为玉米,需在9月20日至10月1日前后进行收割[46],且此时温度也有所下降,因此蒸散量较8月9日大幅降低。11月13日为秋末,温度降低且地表又恢复裸露状态,蒸散值极小,均在1.5 mm/d以下。

2.3 土地覆盖类型对蒸散的影响

图4 Landsat 8与ZY3数据融合的蒸散结果 Fig.4 The results of evapotranspiration of Landsat 8 and ZY3 data fusion

将8月 27日ZY3数据的高分辨率影像和8月9日Landsat 8影像融合(图4)与仅使用8月9日Landsat 8影像(图3f)计算结果对比来看,在大面积农田的中部和水域中部,蒸散值有所增大,而在东部、西北部的荒漠处蒸散值整体上变小,仅在0.5 mm/d以下,在这些地区大部分为纯像元,空间增强后像元内地表类型未发生变化,因此蒸散值发生变化主要是由于植被信息的增强。南部占比例较小的居民地和道路(归在其他类)在空间增强后蒸散值变小,与周围农田、防护林的蒸散值区分度增大。由此可见,在地物类型交界的地方,如农田与荒漠交界处、居民地和道路与农田的交界处,都不再是缓慢地发生变化,区分度变得明显,变化值约在2 mm/d左右；而在面积较大的地方空间增强前后变化程度较小,约在1 mm/d。

将所选8景数据的日期以儒略日的形式表示出来,得到不同土地覆盖类型日平均蒸散值的比较(表2),可知这6种地物蒸散值的季节变化趋势一致,均在7月24日(205)取得最大值。但各土地覆盖类型的季节变化幅度不同,水体的季节变幅最大(6.20 mm/d),其次为农田(5.43 mm/d),而灌木丛的变幅最小,仅为2.66 mm/d。由表中信息还可看出,在3月18日(77)以前,各土地覆盖类型的蒸散值均在缓慢变大且差距不大,而从6月6日(157)开始蒸散量出现明显差异,在9月26日(269)之后,各土地覆盖类型的蒸散值又回到相差不大的状态。蒸散量的多少按照土地覆盖类型排序依次为：水体>农田>防护林>裸地>绿洲外围防风固沙灌木丛,其中水体蒸散值在这几日均为最大,其他类由于包含水体和农田等,蒸散值也较高,在此不列入排序。

表2 不同土地覆盖类型日平均值及变化幅度

裸地与生长在荒漠中的灌木丛群落的蒸散变化趋势极为相似,由于荒漠地区植物生长需依靠降水和地下水作为水分补给来源,而7月份为该研究区雨季,水热条件较好,植物生长迅速,所以7月24日蒸散量较大,其他时间均在2 mm/d以下,由此可见,裸地与灌木丛可代表荒漠地区蒸散量的季节变化。农田的日蒸散值除夏季的7月24日和8月9日外,其他时间均小于林地,这是由于作物在夏季供水充足,作物生长旺盛,蒸散能力极强,因此要高于林地的蒸散值。

3 讨论

从SEBS模型值和通量塔实测值所建立的回归方程来看,斜率接近于1,说明变化趋势较一致,拟合程度较好且两者之间差异不大,表明该模型适用性较强。但是该模型在估算日蒸散值时存在明显的高估现象,在缺水的荒漠地区比农田地区高估现象更严重,这是由于在干旱区植被生长过程中会遭受水分胁迫,因此会低估感热通量,高估蒸散量,Gokmen等[47]在研究中发现：SEBS模型中在计算蒸散量时并未明确考虑到土壤水分的依赖性,对此他们将利用微波遥感获取的土壤水分信息整合进SEBS模型中显热通量的计算,以获取在水分胁迫下的蒸散值,改进后的模型即成为将土壤水分考虑在计算蒸散量内的单源模型。连晋姣等[48]将利用METRIC模型计算的黑河中游生长季日蒸散值与通量数据做对比,同样也发现模型存在高估现象。Ma等[49]利用改进后的SEBS模型估算了黑河中游6—9月农田地区的日蒸散,为4.9 mm/d,与本文在7月24日利用SEBS模型估算得到的农田区日蒸散值6.0 mm/d有一定出入,但两个研究的范围有所不同,亦可能与SEBS模型存在高估现象有关。相比较异质性较强的绿洲-荒漠过渡带,下垫面为均一农田时模型的模拟效果更好,该结果与吴雪娇等[50]利用MODIS/Terra影像,基于SEBS模型在黑河中游估算日蒸散值并用不同下垫面的通量塔数据进行验证而得出的结论一致,也是在下垫面为均一条件时的估算结果更为精准。值得注意的是在将模型计算值与实测值进行对比时所选取的像元个数和位置会影响模拟的精度,Meijninger等[51]曾测试了大口径闪烁仪LAS在非均质地表的适用性并确定了用于模拟的像元数量,本文中为了精确起见,仅选择了通量塔所在经纬度的一个像元值,在实际应用中,可根据分辨率大小来确定所选像元个数。同时,由于通量塔覆盖范围与像元大小并不完全一致,因此蒸散值存在差别是确定的,这也会对验证结果的准确性产生影响,为进一步对模型计算值做精确验证,应选取与通量塔覆盖范围相当的像元大小范围,以消除卫星轨道漂移及其他原因造成的像元变形等问题。由于该研究区只有农田和荒漠两种土地覆盖类型下设有通量塔,其他几种类型的蒸散值无验证数据,仅用两组数据来验证模型的适用性显得过于单薄。

从时间变化来看,夏季热量充足,作物生长需水量大,植被覆盖度增加,且是荒漠地区的雨季,达到水热最佳耦合点,因此是蒸散量最大的季节。蒸散的多少主要是由于下垫面供水状况和植被生长差异造成,绿洲农田和防护林供水充足,植被指数较大,具有高蒸散量,荒漠则恰好相反,为低蒸散值区。为具体说明植被与蒸散之间的关系,将不同土地覆盖类型的NDVI像素点与对应的蒸散量(7月24日影像)进行对比分析,在农田取NDVI为0.74时,模型计算蒸散值为7.14 mm/d,防护林取NDVI为0.44时,蒸散量为4.91 mm/d,灌木丛取NDVI为0.22时,蒸散量为3.19 mm/d,裸地取NDVI为0.13时,蒸散量为2.98 mm/d,由此可说明随着植被覆盖度增加,蒸散量也随之增加。

利用空间增强后的地表温度值求取蒸散量相当于采用较高空间分辨率的数据来解析混合像元内部的情况,能够得到更为准确的地表蒸散值。对比分析空间增强前后的影像,受混合像元的影响,土地覆盖类型交界部分的信息会被夸大或忽略,在交界处一个像元内若存在着裸土和植被至少两种土地覆盖类型,则该部分的土壤热通量因为裸土的存在会比空间增强前增加,因此蒸散值会变小；若交界处一个像元内存在着农田、林地、道路及居民地等地表类型,面积较小的地表类型将会被忽略,而在空间增强后又被考虑进去,显热通量就会增加,蒸散值相应变小。该研究结果与辛晓洲等[52]将MODIS数据与CBERS-02提供的地表类型图相结合并用面积权重平均法得到MODIS像元的通量一文中所得的结论,以及刘雅怩等[53]用同样方法在黑河流域进行地表通量的估算所得出的结论是一致的。王万同等[45]在利用ETM+与MODIS数据对伊洛河流域地表蒸散值进行估算时也发现低分辨率影像中的混合像元容易被忽略或夸大。为提高计算的精度,所选择高分辨率植被信息的时间应尽可能地接近Landsat 8影像的获取时间,这样才可以保证地表植被信息是在相同时相下获取的,但是由于两种数据获取周期的差异,导致很难选择到时间近乎接近的两日数据。并且从整个计算过程来看,仅有植被信息部分使用了高分辨率数据,而影响蒸散的还有多种气象因素,单纯提高某个参数对计算结果的影响是极小的,以上都是空间增强法需要改进的地方。

从各土地覆盖类型的日蒸散量统计分布来看,季节变化趋势是一致的,但变化幅度不尽相同,变化幅度由大而小依次为：水体>农田>防护林>裸地>绿洲外围防风固沙灌木丛。蒸散值达最大的时间都在夏季。夏季是农作物生长的黄金时期,灌溉充足且热量条件丰富,农田和林地的蒸散值明显增加。此时为该地区雨季,裸地水热条件好且荒漠植物也进入生长季节,因此裸地和绿洲外围防风固沙灌木丛的蒸散变化特征可用来表示荒漠地区蒸散量的季节变化。其他类主要包括了居民地和道路,其具有较高蒸散值可能与实际不符,这是由于尽管使用了分辨率较高的ZY3数据来提取地物类型,但不可避免地将农田或林地划分进去,且农田和其他类的地表与近地面大气之间热力性质差异较大,地气能量交换的主要方式不同,居民地、道路等以感热交换为主,近地表形成不稳定层结,而农田、林地主要通过潜热进行交换,因此两者间容易造成平流交换,产生“冷岛效应”[54-55],使得临近居民地的林地和农田的蒸散值更高,也就解释了其他类蒸散值较高的原因。

4 结论

本文利用我国资源三号影像提供的高分辨率植被信息获取了空间增强后的地表温度值,运用SEBS模型对黑河中游绿洲-荒漠过渡带地表蒸散量进行了估算,探讨了不同土地覆盖类型对蒸散量空间分布的影响,得到的主要结论有：

(1)将高分辨率的ZY3影像与Landsat 8数据进行融合后由于植被信息的增强,在地物交界处蒸散量的区分度变得明显；通过与通量数据的观测值对比分析表明：SEBS模型对该地区日蒸散量进行估算可行,且在均一地表类型时精度更高。

(2)从季节变化来看,各土地覆盖类型的季节变化趋势一致(夏季>春季>秋季>冬季),但变幅存在差异：水体>农田>防护林>裸地>绿洲外围防风固沙灌木丛,主要与下垫面供水、热量收支和植被覆盖有关。

(3)在空间分布上,除水体外,地表植被覆盖对蒸散量影响明显,夏季绿洲内部农田和防护林蒸散量较高,最大值分别为5.95 mm/d、5.14 mm/d,是绿洲以外荒漠裸地的2—3倍。