图3.6显示了大沽河流域地表蒸散量的逐年变化过程和年际波动情况。年均蒸散量超出多年ET 平均值的年份出现在2005 年、2007年、2008年、2010年、2012年和2013年,其中2008年和2013年尤为突出,分别超出多年平均值50.29mm 和22.34mm,其他年份均低于多年平均值。总体上来看,大沽河流域地表蒸散量年内分布呈现先增大后减小的单峰型分布趋势,蒸散量月均变化曲线见图3.7。图3.7大沽河流域2005—2014年不同月份地表蒸散量分布......
2023-09-17
(1)验证与分析。为了验证SEBAL模型反演大沽河流域日蒸散量结果的可信度,将研究区内各气象台站测得的蒸发量进行相应计算得到实际的蒸散量,以其作为真值,将反演结果与蒸散量实测值进行对比,对比结果见图3.1。根据比较分析,总体来说,经过SEBAL 模型估算的日蒸散量均值与实测的日蒸散量均值结果较为相近,平均绝对误差(MAE)为0.733mm,均方根误差(RMSE)为1.0078mm。图3.2为蒸散量估算值与实测值之间的散点图,拟合优度R2为0.8988。
图3.1 遥感估算日蒸散量与气象台实测值对比
图3.2 ET 估算值与实测值的散点图
通过上述对比结果,可以验证模型反演的正确性,因此可将SEBAL模型应用于大沽河流域的蒸散量研究中。利用遥感方法估算出来的日蒸散量和实测数据之间难免存在一定的误差,3月和6月误差相对稍大。造成误差的原因可能是模型在用瞬时蒸散量推算日蒸散量的过程中,采用的是蒸散比率不变法,此方法的基本假设是能量通量的分配在24h内服从相同的比例,而实际上能量通量的分配并不为一个固定的数值。如果有外界热量平流输入的影响,如风速和云的变化,则会破坏这种稳定的状态,使能量构成比例发生变化。此外,虽然选取的MODIS数据产品都经过严格筛选,数据质量较好,但少量地区仍然会受到云覆盖等因素的影响,导致反演蒸散量时会有较小比例的无值区,给反演结果带来一定误差。
(2)日蒸散量空间分布。选取大沽河流域2013年数据质量较好的12幅影像数据进行蒸散量的反演,当日蒸散量(mm)分布情况见图3.3。
从蒸散量的时间分布上来看,大沽河流域在冬季蒸散量最少,12月、1月、2月日蒸散均值分别在1.2mm、0.95mm、0.97mm。春季蒸散量次之,日蒸散量均在3mm 左右。夏秋两季蒸散量最多,日均蒸散量均超过3mm,在8月达到峰值,日蒸散量达到了6.02mm。
从蒸散量的空间分布上来看,流域主要的植被覆盖类型为耕地,面积约为4761km2,占整个流域面积的77.5%;其次为城镇及建筑用地,占整个流域面积的13.1%;林地占流域面积的6.8%,水体占流域面积的2.1%,未利用土地占流域面积的0.5%。在枯水期,由于地表土壤水分含量较低,流域总体蒸散量较小,流域东部及西南部部分地区蒸散量较小。在丰水期,由于降水较多和植被覆盖状况良好,流域整体蒸散量较高,且流域北部以林地为主,蒸散量高于其他地区。在全年中,胶州市所处地段日蒸散量相对较高,这主要是由于胶州市以耕地为主,植被覆盖率占总面积的71%,且境内包含大沽河、胶莱河、洋河等骨干河流。河流周围冲积平原和涝洼地区分布着大量湿生植被,南部丘陵区则分布着梯田等旱生植被。植被和水体所占较大比例决定了胶州地区具有相对较高的蒸散量。在1年内的各个月中,各土地利用类型的平均蒸散量之间的大小关系不完全是固定不变的,对于同一个月的不同土地利用类型,其平均蒸散量也不相同。以2013年7月24日的反演结果为例,研究大沽河流域不同土地利用类型日蒸散量的大小关系。水体的平均蒸散量最大,为4.71mm,林地蒸散量仅次于水体,为4.61mm,主要分布在流域北部的林地区。城镇及建筑用地蒸散量最小,为4.15mm。该结果符合地物蒸散的一般规律。研究区各种土地覆盖类型的日蒸散量的统计特征见表3.1。
图3.3(一) 2013年大沽河流域日蒸散量遥感估算图
图3.3(二) 2013年大沽河流域日蒸散量遥感估算图
图3.3(三) 2013年大沽河流域日蒸散量遥感估算图
表3.1 大沽河流域各类土地利用类型日蒸散量遥感估算结果 单位:mm
(3)年内蒸散量及变化规律。图3.4为大沽河流域2013年内逐月的蒸散量变化曲线图。1月、2月气温较低,蒸发能力较弱,且土壤含水量低,耕地多处于裸土状态,因而流域蒸散量为一年当中的最低值,月均蒸散量分别为20.01mm、22.83mm;4月、5月为多种农作物生长期,土壤供水充足,蒸发旺盛,地表蒸散量明显增加,月均蒸散量分别为69.40mm、70.01mm;蒸发量最大值出现在8月,夏季降雨较多,热量条件充足,蒸发强烈,月均值达到了141.62mm。从10月以后,随着农作物的收获、植被的枯萎凋落等原因,流域蒸散量逐步降低,到12月蒸散量降低至30.34mm。
图3.4 2013年大沽河流域逐月遥感估算蒸散量曲线图
以3—5月的蒸发量作为春季蒸散量,6—8月的蒸散量作为夏季蒸散量,9—11月的蒸散量作为秋季蒸散量,12月至次年2月的蒸散量作为冬季蒸散量,则大沽河流域各季的蒸散量的大小关系柱状图见图3.5。从图3.5中可以看出大沽河流域蒸散量年内变化基本符合如下规律:夏季大于秋季大于春季大于冬季。蒸发量最大值出现在夏季,总蒸散量达到308.85mm;其次为春季和秋季,各季度蒸散量分别为197.98mm、241.64mm;蒸发量最小值出现在冬季,总蒸散量仅为72.88mm。这是由于冬季太阳辐射时间最短,夏季太阳辐射时间最长,且夏季植被生长旺盛,气温较高,因此日蒸散量在冬季较低,在夏季较高。该结果较好地体现出了流域的自然气候条件特征和地表蒸散的一般规律。
图3.5 2013年大沽河流域各季度遥感估算蒸散量柱状图
有关大沽河流域水文要素监测体系建设与实践的文章
图3.6显示了大沽河流域地表蒸散量的逐年变化过程和年际波动情况。年均蒸散量超出多年ET 平均值的年份出现在2005 年、2007年、2008年、2010年、2012年和2013年,其中2008年和2013年尤为突出,分别超出多年平均值50.29mm 和22.34mm,其他年份均低于多年平均值。总体上来看,大沽河流域地表蒸散量年内分布呈现先增大后减小的单峰型分布趋势,蒸散量月均变化曲线见图3.7。图3.7大沽河流域2005—2014年不同月份地表蒸散量分布......
2023-09-17
黄同水库的总库容6027 万m3, 兴利库容2450 万m3,主要为农田灌溉服务。南胶莱河在北王珠镇刘家花园处流入胶州市,经店口乡沽河汇入大沽河,干流全长30km,流域面积1500km2。主要支流有胶河、墨水河及清水河等。大沽河流域水资源丰富,多年平均河川径流量为6.311 亿m3,现已建成中型水库8座,小型水库90座,塘坝拦河闸1223座,地表总拦蓄能力达4.143亿m3。表2.1大沽河流域水库一览表......
2023-09-17
选取土地利用/土地覆盖变化范围较大的4 个年份的数据做监测分析。表7.10大沽河流域2001年与2004年变化监测结果统计表7.11大沽河流域2004年与2007年变化监测结果统计表7.12大沽河流域2007年与2012年变化监测结果统计根据上述统计数据可以看出,在2001—2012年间青岛大沽河流域的7种土地利用类型的面积变化情况。......
2023-09-17
大沽河流域蒸散量估算选用了2013年数据质量较好的12 幅MODIS 的Landsat-8遥感数据,解译获取土地利用信息,采用SEBAL模型进行蒸散量估算,得到大沽河流域蒸散量情况;将遥感监测估算的大沽河流域蒸散发量与研究区内各气象台站实测蒸散量进行对比,验证SEBAL 模型反演大沽河流域日蒸散量结果的可信度。......
2023-09-17
传统的蒸散发估算方法局限于局部尺度,而较大空间尺度内陆面特征和水热传输的非均匀性导致传统的估算方法难以获取区域尺度的蒸散发。遥感以少量的地面数据估算大范围区域蒸散的优势被认为是可以经济、有效地提供区域蒸散发消耗量的唯一方法。该公式被证明在致密冠层的蒸散发估算上有良好的效果。除Penman类模型外,R.J.Bouchet于1963年提出了陆面实际蒸散与可能蒸散之间的互补关系原理,为估算实际蒸散发开辟了新的途径。......
2023-09-17
使用Matlab编程,将2个时相8d合成的地表温度数据合成为1个时相16d合成的地表温度数据。利用上一步骤中提取的数据,在Excel中对NDVI和最大和最小陆地表面温度进行线性拟合,获得干边和湿边方程的系数a1、b1、a2和b2。使用Matlab编程,根据式计算图像上每个像元的温度植被干旱指数值,获取流域TVDI的分布图,根据TVDI等级划分形成大沽河流域土壤湿度分布图。图5.4显示了部分Matlab程序文件列表。......
2023-09-17
在Matlab中编写程序,提取大沽河流域相同NDVI 下的不同像元对应的最大陆地表面温度和最小陆地表面温度,每期数据中NDVI 均匀划分了256个等级,获得研究区每16d的Ts-NDVI 特征空间,共23期结果。图5.52013年每16d的Ts-NDVI 空间散点图(部分)干湿边方程确定。图5.7(一)大沽河流域2013年每16d的土壤湿度等级分布图图5.7(二)大沽河流域2013年每16d的土壤湿度等级分布图从图5.7中可以看出,大沽河流域土壤湿度时间、空间分布均具有不均匀性。......
2023-09-17
MODIS每一个仪器的设计寿命为5年,将计划发射4颗卫星。由此估计,利用MODIS仪器至少将获得5年36个光谱波段的地球综合信息。NASA 提供MODIS全球数据产品,共有44种标准产品,具有不同的时 间和空间分辨率,均由DAAC 存储和发布。该产品中的地表温度是通过建立31通道、32通道亮温线性组合的劈窗算法计算获取的,其中通道亮温值是根据辐射度与0.1K 步长亮温的查找表来确定,在计算地表温度过程中需要的发射率是根据MODIS土地覆盖产品确定的。......
2023-09-17
相关推荐