An analysis of spatio-temporal characteristics and influencing factors of surface evapotranspiration in the Yinchuan Plain based on MOD16 data
-
摘要:
地表蒸散发是陆地水文循环的重要组成部分，分析蒸散量时空变化特征是深入了解干旱区水文过程的基础。由于银川平原缺乏区域尺度实际蒸散量的长期观测，很难得到长时间序列蒸散量的时空变化特征。基于MOD16A3地表蒸散量数据及研究区内气象站点实测数据，采用Theil Sen Median趋势度分析、MK突变检验及CA-Markov模型等方法，从时间与空间的角度分析2004—2019年银川平原地表蒸散量的变化特征及影响因素，预测2024年地表蒸散量的发展趋势。研究结果表明：2004—2019年银川平原蒸散量年际波动总体是增加趋势，MK突变检验结果显示2010年是蒸散量时序数据的突变点；银川平原实际蒸散量与潜在蒸散量空间分布格局、变化趋势均存在明显的差异性，蒸散量在近16年呈增加趋势，潜在蒸散量呈减少趋势，符合干旱区蒸散发互补相关理论。采用CA-Markov模型对2024年银川平原地表蒸散量未来发展趋势进行预测，模拟结果显示在未来5年银川平原蒸散量仍呈增加趋势；蒸散量的时空变化受气候与人类活动的共同影响，蒸散量与气温、降水、日照时数呈正相关，与相对湿度呈负相关，土地利用结构影响年蒸散量的空间格局，呈现出水田>旱田>林地>草地>荒漠的规律。
-
关键词:
- MOD16 /
- 蒸散发 /
- Theil-Sen Median趋势度 /
- CA-Markov模型 /
- 银川平原
Abstract:Surface evapotranspiration is an important part of terrestrial hydrological cycle, and the analysis of spatio-temporal variation characteristics of evapotranspiration is the basis of in-depth understanding of hydrological processes in arid areas. Due to the lack of long-term observations of actual regional evapotranspiration in the Yinchuan Plain, it is difficult to obtain the spatio-temporal variation characteristics of long-term series evapotranspiration. Based on the MOD16A3 surface evapotranspiration data and the measured data of meteorological stations in the study area, this paper analyzes the variation characteristics and influencing factors of surface evapotranspiration in the Yinchuan Plain from the perspective of time and space by using the methods of Theil Sen median trend analysis, MK mutation test and CA-Markov model, and predicts the development trend of surface evapotranspiration in the next five years. The results show that from 2004 to 2019, the interannual variation trend of evapotranspiration in the Yinchuan Plain is increasing, and the MK mutation test results show that the year of 2010 is the mutation point of evapotranspiration time series data. There are obvious differences in spatial distribution pattern and change trend between the actual evapotranspiration (ET) and potential evapotranspiration (PET) in the Yinchuan Plain. Actual evapotranspiration shows an increasing trend in recent 16 years, and potential evapotranspiration shows a decreasing trend, which conforms to the theory of complementary evapotranspiration in arid areas. The future development trend of surface evapotranspiration in the Yinchuan Plain in 2024 is predicted with the CA-Markov model. The simulation results show that the evapotranspiration will still increase in the next five years. The spatio-temporal variation of evapotranspiration is affected by climate and human activities. Evapotranspiration is positively correlated with temperature, precipitation and sunshine hours, and is negatively correlated with relative humidity. Land use structure affects annual evapotranspiration. The results also show the rule: Evapotranspiration of paddy field > dry land > woodland > grassland > desert.
-
Key words:
- MOD16 /
- evapotranspiration /
- Theil Sen median trend /
- CA-Markov model /
- Yinchuan Plain
-
蒸散发是水文循环过程中极其重要的环节，同时也是自然界水资源管理中最基本、最不确定的要素。蒸散发将地表水分与地表能量平衡联系起来[1]，由同时发生的两个过程控制，即土壤蒸发和叶面蒸腾。潜在蒸散发是指一定气象条件下，充分供水地表（植被、裸土或开阔水域）的蒸散发，是表征区域干旱程度的重要指标之一。
由于地表的高度异质性，大多数传统方法通常只提供点尺度的估算，不能有效进行区域尺度的计算。相较传统的监测方法，遥感技术可以针对地表的异质性和传热过程的动态性进行区域尺度的蒸散量监测[2]，被认为是量化蒸散发时空变化与空间尺度扩展的重要工具。目前采用遥感技术进行区域尺度蒸散量反演的方法主要有经验统计公式、能量余项模型、特征空间法等[3]。基于区域尺度蒸散量计算方法，众多学者及团队研发了蒸散量遥感数据产品，如NASA基于陆面模型提供的全球陆面通量数据GLDAS[4]；英国布里斯托大学Miralles等[5]基于Priestley-Taylor方程发布的GLEAM产品；美国地质调查局（USGS）基于简单化的陆面能量平衡模型（SSEB），以Penman-Monteith公式为理论支撑从而获取实际蒸散量的SSEBop/ET产品[6]；美国航空航天局（NASA）基于MODIS卫星遥感数据，考虑多种阻力算法制作的全球陆地蒸散量产品MOD16等[7]。MOD16产品因具有高时空分辨率、获取方式简单、对蒸散量遥感监测有较好的反演精度等特点而广受国内外学者认可，Sebinasi Dzikiti等[8]以南非干旱区为研究区，比较了MOD16和PT-JPL两种蒸散量估算模型的精度，结果表明与PT-JPL模型相比，MOD16具有准确分析干旱环境下蒸散量的潜力；Hassan等[9]采用MOD16数据与SEBS模型反演的埃及月蒸散量进行比较，两者存在高度相关性；何慧娟等[10]探讨了陕西省地表蒸散量的时空变化特征，证实MOD16数据具有较高的精度，满足陕西省蒸散发研究要求。
由于降水量低、年平均气温高以及延长的干旱期等不利条件，银川平原地区长期存在水资源空间分布不均和对传统灌溉技术过于依赖等问题，因此需要对区域水资源进行精确地估算及管理[11]。以往蒸散发研究多针对场地尺度，然而受黄河水系分布、降水等因素影响，银川平原地表蒸散量空间差异性极大，场地尺度研究难以反映区域尺度上的时空变化规律。
本文以MOD16蒸散量数据、地面气象实测数据为基础，对2004—2019年银川平原地表蒸散量的时空变化特征进行分析、模拟，并预测了2024年地表蒸散量的发展趋势。研究结果可为掌握银川平原水循环规律、合理开发利用水资源以及保护生态环境提供科学依据。
1. 研究区和数据来源
1.1 研究区概况
银川平原位于宁夏回族自治区中部黄河两岸，东依鄂尔多斯高原，西至贺兰山，南接黄土高原，北止石嘴山，南北长度约为165 km，东西长10～50 km（图1）。银川平原属于温带干旱区，年均日照时数约3000 h，年平均降水量仅约200 mm，是我国典型的引黄灌溉区，黄河过境水量年均达300×108 m3。与此同时，银川平原年蒸发量高达1600 mm，强烈的蒸发浓缩作用加剧了当地盐渍化的发展，成为制约农业发展的一个重要因素[12]。
1.2 数据源及数据预处理
本文主要采用MOD16数据进行银川平原蒸散量时空特征分析，并结合气象实测数据、Landsat数据探讨影响因素。
MOD16数据是NASA基于Penman-Monteith模型[7]，综合考虑植被蒸腾、土壤表面蒸发等条件提供的全球蒸散量数据集。本文选用2004—2019年MOD16A3数据，空间分辨率500 m，时间分辨率为年（
https://search.earthdata.nasa.gov ）。利用MRT工具对该数据集进行投影转换、波段提取等预处理，并根据NASA提供的MOD16A3产品说明对水体、建筑用地等无植被覆盖区进行无值化处理，将其转换成WGS-84基准面、UTM投影的蒸散量、潜在蒸散量初始数据。气象数据源自中国气象数据共享网（
https://data.cma.cn/ ），筛选银川平原内9个站点数据，主要包括气温、降水、日照时数、相对湿度及蒸发量气象数据。统计分析2004—2019年实测蒸发皿数据并处理成年值数据，用于MOD16蒸散量在银川平原的适用性验证。Landsat数据源自地理空间数据云网站
http://www.gscloud.cn/ ，选取2004年、2014年7月份遥感影像，根据《土地利用现状分类》国家标准并结合银川平原实际土地利用情况，采用支持向量机及面向对象的分类方法，将银川平原土地划分为水田、旱田、林地、草地、水体、建设用地及荒漠化土地7类，并通过实测调查取点验证，分类精度满足研究需求。2. 研究方法
2.1 Theil-Sen Median趋势度
Theil-Sen Median趋势度方法主要用于蒸散量时序数据在像元尺度上的趋势分析[13-14]。该方法与常用的一元线性回归方法不同，其采用时序数据点的中值进行趋势分析，减小了异常值造成的误差。非参数Sen斜率度结果可用来评估研究时段内该时序数据单位时间的变化：
β=Median(Sj−Sij−i),∀j>i (1) 式中：Si、Sj——研究时段中第i、j年对应的数据。
趋势度结果β通过正负反映该时序数据增加、下降的变化趋势。
2.2 Mann-Kendall检验
Mann-Kendall检验是一种以非参数检验为基础，进行时间序列数据趋势分析的检验方法。世界气象组织（WMO）推荐使用这种方法来评估水文气象时间序列是否存在单调趋势[15]。该方法将非参数检验与非线性趋势检验相结合，使结果合理化，主要用于蒸散量时序数据趋势变化显著性分析及突变检验[16]。
Mann-Kandall方法通过检验统计量Z值评估蒸散量时序数据显著趋势的存在：
Z={S−1√Var(S),S>00,S=0S+1√Var(S),S<0 (2) S=∑n−1i=1∑nj=i+1sgn(xj−xi) (3) sgn(xj−xi)={1,xj>xi0,xj=xi−1,xj<xi (4) Var(S)=n(n−1)(2n+5)18 (5) 式中：S——Kendall检验统计量；
xi、xj——第i、j年对应的蒸散量；
n——时序数据的长度，本文为16。
在特定α显著性水平下，采用双边趋势检验，根据正态分布表中
Z1−2/α 的值判断蒸散量趋势变化是否显著。若|Z|>|Z1−2/α| ，该时序数据趋势变化通过了α显著性水平。本文选取0.01，0.05，0.10的α水平进行显著性检验，Z1−2/α 的值分别为2.576，1.960，1.645，根据3个断点将蒸散量变化趋势分为极显著、显著、弱显著以及不显著4种状态[17]。Mann-Kandall可通过对蒸散量时序数据构造一组秩序列，用于总体趋势分析中的突变点检验。秩序列计算公式为：
sk=∑ki=1ri (6) ri={1,xi>xj0,xi⩽xj(j=1,2,⋯,i) (7) 此时设定蒸散量时序数据是独立、互不干扰的，UF检验统计量为：
UFk=sk−E(sk)√Var(sk)(k=1,2,⋯,n) (8) E(sk)=n(n+1)4 (9) Var(sk)=n(n−1)(2n+5)72 (10) 通过计算，得到蒸散量UFi曲线。通过UFi曲线与α显著性水平下对应的正态分布表临界值
Uα 对比：若|UFi|>Uα ，意味着蒸散量时序数据存在显著的趋势变化。将逆序的蒸散量时序数据X\left({x}_{n},\cdots , {x}_{2},{x}_{1}\right) 重复以上计算步骤，并将计算结果乘以−1，得到统计量曲线UBi。通过两条统计量曲线在α置信区间内的交点个数及位置，可以判断该时序数据突变发生的次数及相应的时间。2.3 CA-Markov模型
CA-Markov模型可用于蒸散量时序数据的模拟预测。元胞自动机（CA）设定元胞在时间、空间均处于独立、离散的状态，是一种能够模拟复杂系统空间演化过程的模型。该模型基本思想为：复杂系统内的任一元胞互不干扰，执行各自对应的局部转换规则。CA模型计算公式为：
St+1=f(St,N) (11) 式中：St、St+1——复杂系统内t、t+1时间的元胞状态；
N——设定的元胞邻域范围；
f——从t到t+1时刻元胞执行的转换规则。
CA-Markov模型在元胞自动机对空间格局精确模拟的基础上，结合Markov模型对时序数据预测的显著优势，实现了模拟结果在时间、空间上的合理化[18]。本文根据2004年、2009年、2014年、2019年4期年蒸散量数据统计特征进行分级，采用Markov模型计算相应时段的转移概率矩阵，以5 a为反演时段，对2024年银川平原的蒸散量进行模拟预测。
3. 结果与分析
3.1 MOD16数据适用性分析
在采用MOD16数据对银川平原蒸散量时空特征分析之前，需对数据适用性进行检验。本文采用实测蒸发量与潜在蒸散量的相关分析进行适用性检验。由于MOD16数据是年尺度，因此，将气象站点蒸发皿日值数据处理成年值，解决了气象数据中异常值及不同口径大小对结果的影响，实现两种数据在时间尺度上的统一。提取该站点所在像元点的潜在蒸散量，在点尺度上进行各个气象站点蒸发量实测值与潜在蒸散量的相关分析。结果（图2）显示气象站实测蒸发量与潜在蒸散量在0.01水平上显著相关，相关系数为0.88，表明MOD16数据可用于银川平原地区蒸散量的相关研究。
3.2 蒸散量动态特征
2004—2019年银川平原蒸散量年际变化见图3。蒸散量整体上波动较大，增加趋势明显。年际蒸散量为270.3～396.1 mm，多年平均蒸散量为330.7 mm。研究时段内蒸散量距平相对变化率波动强烈，为−18.3%～19.8%。2004—2009年蒸散量低于研究时段内平均水平，2010—2019年（2013年除外）蒸散量均高于多年平均水平。其中2005年蒸散量最小，为270.3 mm，2016年蒸散量最大，为396.1 mm，差值达到了125.8 mm。
图 2. MOD16潜在蒸散量与蒸发皿实测年蒸发量拟合关系Figure 2. Correlation between potential evapotranspiration of MOD16 and measured annual evapotranspiration of pan
图 3. 2004—2019年银川平原蒸散量年际变化Figure 3. Interannual variation of evapotranspiration in the Yinchuan Plain from 2004 to 2019在分析时序数据时，除分析蒸散量总体变化趋势外，还需要判断并检验该时序数据在研究时段内突变发生的时间以及频率。在Matlab中采用Mann-Kendall法对2004—2019年银川平原年蒸散量序列的突变情况进行检验（图4）。结果表明：（1）近16 年来，年平均蒸散量整体呈现波动增加的趋势；（2）2010—2019年UF曲线超过α=0.05显著性水平临界值（1.96），说明该时段年蒸散量增加趋势显著；（3）UF与UB曲线在置信区间（α=0.05）内有1个交点，位于2010年，说明银川平原年平均蒸散量在2010年左右发生了明显突变，突变后年蒸散量增加了45.2 mm。突变前2004—2009年时间段年蒸散量均低于多年平均水平，2010年突变后高于多年平均水平。
3.3 蒸散量空间变化特征
分析银川平原多年平均蒸散量、潜在蒸散量基本空间特征（图5）可知，2004—2019年银川平原多年平均蒸散量、潜在蒸散量空间分布具有明显差异。多年平均蒸散量变化范围为175.7～577.8 mm，整体呈现出南高北低、中间高两侧低的空间格局。从各市行政单元看，吴忠市蒸散量最高，银川市次之，石嘴山市蒸散量最低。多年平均潜在蒸散量为1429.3～1929.1 mm，整体呈现南低北高、中间低两侧高与蒸散量分布特征相反的差异特性。Bouchet等[19]认为受区域能量的影响，实际蒸散量减少，剩余能量将用于增加潜在蒸散量。这一区域蒸散互补相关理论可以解释蒸散量、潜在蒸散量空间分布的差异。多年平均潜在蒸散量以银川平原西侧荒漠带及石嘴山市最高，低值位于吴忠市、银川市内。蒸散干旱指数（EDI）[20]，即整数1与蒸散量/潜在蒸散量之差，可初步描述银川平原的干旱程度，该指数理论上量值范围为0～1，EDI值越大，说明该地区地表水分缺失越严重。图5（c）显示，银川平原地区整体干旱程度严重，EDI波动范围在0.66～0.90，石嘴山市最为干旱，银川市次之。吴忠市内黄河支系发达，整体地表水分较为充足，受水分胁迫程度较轻。
图 5. 银川平原2004—2019年年平均蒸散量、潜在蒸散量空间分布Figure 5. Spatial distribution of annual average evapotranspiration and potential evapotranspiration in the Yinchuan Plain from 2004 to 2019从银川平原2004—2019年蒸散量和潜在蒸散量趋势显著性结果分析可知（图6），银川平原蒸散量整体呈增加趋势，银川、吴忠、青铜峡城市周边区域蒸散量呈下降趋势。潜在蒸散量变化趋势与蒸散量趋势相反，石嘴山市零星区域潜在蒸散量呈上升趋势。整体上近16年来银川平原潜在蒸散量下降的同时，蒸散量在相对增加，符合区域蒸散互补相关理论，并从蒸散干旱指数方向分析，该地区干旱情况较之前有明显好转。
图 6. 2004—2019年银川平原蒸散量、潜在蒸散量变化趋势空间分布Figure 6. Spatial distribution of evapotranspiration and potential evapotranspiration in the Yinchuan Plain from 2004 to 2019银川平原绝大部分区域蒸散量呈增加趋势，占比达92.96%，极显著增加的范围最大，达75.49%；其次分别是不显著增加、显著增加及弱显著增加（表1）。蒸散量呈减少趋势的区域中以不显著减少为主，面积占比为4.41%；与蒸散量变化趋势相反，银川平原99.87%的区域潜在蒸散量呈下降趋势，极显著下降变化的面积占比最大，达53.68%，潜在蒸散量呈上升趋势的区域面积占比仅0.13%。
表 1. 银川平原蒸散量、潜在蒸散量变化趋势及显著性统计Table 1. Change trend and significance of evapotranspiration and potential evapotranspiration in the Yinchuan Plain变化趋势 显著性面积占比/% 总计/% 不显著 弱显著 显著 极显著 蒸散量 增加 7.69 2.69 7.09 75.49 92.96 减少 4.41 0.81 1.09 0.73 7.04 潜在蒸散量 增加 0.13 0 0 0 0.13 减少 3.98 6.1 36.12 53.68 99.87 | Show Table
DownLoad:
CSV
3.4 基于CA-Markov模型的蒸散量预测
基于CA-Markov模型，以2004年、2009年、2014年、2019年4期年蒸散量为基础数据，按照银川平原实际蒸散量统计特征，以200，300，400，500，600 mm为断点分为6个等级。通过Markov模型计算2004—2009年、2009—2014年、2014—2019年转移概率矩阵，选取合适的滤波器，进行预测次数的迭代，完成银川平原2009年、2014年、2019年蒸散量的预测。对模型模拟结果与实际蒸散量数据进行Kappa系数精度验证，精度检验结果显示2009年、2014年、2019年总体模拟精度分别为0.84，0.78，0.82，模拟精度较好，CA-Markov模型预测结果可信。在此基础上，以2019年蒸散量为初始数据，对银川平原2024年蒸散量空间分布进行模拟。
从2024年银川平原蒸散量模拟结果看（图7），低值主要分布在银川市、石嘴山市；与2019年相比，蒸散量整体上由低等级向高等级转移，其中银川平原北部以200～400 mm低蒸散量等级向400～500 mm较高蒸散量等级转移为主，南部主要由较高蒸散量等级向500～600 mm高蒸散量等级转移。
图 7. 2024年银川平原蒸散量 CA-Markov预测图Figure 7. CA-Markov prediction of evapotranspiration in the Yinchuan Plain in 2024银川平原蒸散量等级面积占比统计结果见表2，100～300 mm蒸散量较低值区面积逐渐减少，300～400 mm蒸散量中等区面积基本持平，400～700 mm 蒸散量高值区面积持续增加。2024年CA-Markov模型模拟结果表明蒸散量整体呈增加的趋势。
表 2. 银川平原蒸散量等级面积占比Table 2. Proportion of the evapotranspiration grade area in the Yinchuan Plain蒸散量等级/
mm面积占比/% 2004年
实际值2009年
实际值2014年
实际值2019年
实际值2024年
预测值100～200 27.83 19.79 7.05 0.09 0 200～300 33.56 32.21 31.33 29.45 19.03 300～400 28.36 33.49 32.22 31.36 31.65 400～500 9.94 14.07 23.19 28.25 32.75 500～600 0.31 0.42 6.12 10.02 14.69 600～700 0 0 0 0.80 1.78 | Show TableDownLoad:
CSV
3.5 蒸散量时空变化影响因素分析
蒸散量的时空变化受自然因素与人类活动的共同影响，气候变化影响区域水热传输过程，土地利用格局影响地表水分的供给能力。
2004—2019年年均蒸散量与各气象因子年均值进行相关分析（表3），在时间尺度上探究银川平原蒸散量与各气象因子的关系。研究结果表明：蒸散量与气温、降水、日照时数呈正相关关系，相关系数分别为0.910，0.905，0.768；与相对湿度呈负相关关系，相关系数为−0.746，相关系数均通过了P=0.01的显著性检验。相关系数表明银川平原蒸散量对气温、降水具有更高的敏感性。银川平原日照充足，年均气温呈现逐渐增加趋势，这为蒸散发过程提供了良好的热力环境，而相对湿度增加限制了该地区的蒸散需求，由于蒸散量对气温具有高敏感性，银川平原蒸散量呈增加的趋势。银川平原年降水量较少，分布不均，除南部及黄河沿岸地区外，降水是蒸散发过程主要的水分来源，受降水控制，蒸散量呈现与降水分布一致的南高北低的空间格局。
表 3. 银川平原蒸散量与气象因子相关关系Table 3. Correlation between evapotranspiration and meteorological factors in the Yinchuan Plain气象因子 气温 降水 日照时数 相对湿度 相关系数 0.910 0.905 0.768 −0.746 | Show TableDownLoad:
CSV
为进一步掌握土地利用结构对蒸散量空间分布格局的影响程度，利用2004年、2014年银川平原土地利用遥感解译图分析不同土地利用类型的蒸散量特征。
从2004年、2014年银川平原土地利用类型空间分布图（图8）来看，研究区土地类型以水田、旱田为主，呈现出南部水田、北部旱田的空间格局，这与多年平均蒸散量南高北低的空间分布特征相同。耕地的蒸散量是银川平原区域蒸散量的主要组成部分[21]，其中水田受地表水分胁迫程度较小，蒸散量高于旱田。荒漠化土地主要为裸土地、沙地等，地表裸露，土壤水分含量较少，蒸散量相对较低。草地多分布在荒漠化土地与其他土地类型过渡地带，降水是该类型地表水分的主要来源，水分缺乏，蒸散作用受限制。采用空间分析法统计各土地类型的平均年蒸散量（图9）。由于MOD16数据在水体及城镇用地的蒸散量无数据，为避免统计误差，在对各土地类型进行统计时不考虑这两类。两期土地利用类型的年蒸散量均呈现出水田>旱田>林地>草地>荒漠的规律，且2014年各土地类型的蒸散量均高于2004年，这一现象与银川平原近16年蒸散量整体变化趋势相同。
图 8. 2004年、2014年银川平原土地利用类型Figure 8. Land use types of the Yinchuan Plain in 2004 and 20144. 结论
（1）银川平原2004—2019年蒸散量年际变化较大，整体呈增加趋势，年均蒸散量为270.3～396.1 mm，多年平均蒸散量为330.7 mm，距平相对变化率波动较大，其中2010年是银川平原蒸散量年际变化的突变年。
（2）2004—2019年银川平原蒸散量、潜在蒸散量的空间格局、变化趋势存在明显的空间差异。多年平均蒸散量南高北低、东高西低，潜在蒸散量分布特征与之相反；从变化趋势看，蒸散量近16年呈增加趋势，潜在蒸散量呈减少趋势。CA-Markov模型预测结果显示未来5年银川平原蒸散量仍呈增加趋势。
（3）银川平原蒸散量与各气象因子的相关分析结果表明：蒸散量与气温、降水、日照时数呈正相关，相关系数分别为0.910，0.905，0.768，与相对湿度呈负相关，相关系数为−0.746。其中，蒸散量对气温、降水具有更高的敏感性。
（4）2014年各土地利用类型的平均蒸散量明显高于2004年，这一结果与蒸散量总体趋势一致。两期土地利用类型的年蒸散量均呈现出水田>旱田>林地>草地>荒漠的规律。
MOD16数据能够对银川平原多年蒸散量进行监测，但由于大尺度陆面蒸散发过程复杂，仍可能存在偏差。为进一步深入研究陆面蒸散发过程及定量分析各影响因素对蒸散量的影响，在该地区建立蒸散量实测系统，提高区域蒸散量遥感监测水平十分必要。
-
表 1 银川平原蒸散量、潜在蒸散量变化趋势及显著性统计
Table 1. Change trend and significance of evapotranspiration and potential evapotranspiration in the Yinchuan Plain
变化趋势 显著性面积占比/% 总计/% 不显著 弱显著 显著 极显著 蒸散量 增加 7.69 2.69 7.09 75.49 92.96 减少 4.41 0.81 1.09 0.73 7.04 潜在蒸散量 增加 0.13 0 0 0 0.13 减少 3.98 6.1 36.12 53.68 99.87
下载: 导出CSV
表 2 银川平原蒸散量等级面积占比
Table 2. Proportion of the evapotranspiration grade area in the Yinchuan Plain
蒸散量等级/
mm面积占比/% 2004年
实际值2009年
实际值2014年
实际值2019年
实际值2024年
预测值100～200 27.83 19.79 7.05 0.09 0 200～300 33.56 32.21 31.33 29.45 19.03 300～400 28.36 33.49 32.22 31.36 31.65 400～500 9.94 14.07 23.19 28.25 32.75 500～600 0.31 0.42 6.12 10.02 14.69 600～700 0 0 0 0.80 1.78
下载: 导出CSV
表 3 银川平原蒸散量与气象因子相关关系
Table 3. Correlation between evapotranspiration and meteorological factors in the Yinchuan Plain
气象因子 气温 降水 日照时数 相对湿度 相关系数 0.910 0.905 0.768 −0.746
下载: 导出CSV
-
[1] TIAN F, QIU G Y, YANG Y H, et al. Estimation of evapotranspiration and its partition based on an extended three-temperature model and MODIS products[J]. Journal of Hydrology,2013,498:210 − 220. doi: 10.1016/j.jhydrol.2013.06.038
[2] CHEN Y, XIA J Z, LIANG S L, et al. Comparison of satellite-based evapotranspiration models over terrestrial ecosystems in China[J]. Remote Sensing of Environment,2014,140(1):279 − 293.
[3] 王军, 李和平, 鹿海员. 基于遥感技术的区域蒸散发计算方法综述[J]. 节水灌溉,2016(8):195 − 199. [WANG Jun, LI Heping, LU Haiyuan. Review of regional evapotranspiration estimation models basing on the remote sensing[J]. Water Saving Irrigation,2016(8):195 − 199. (in Chinese with English abstract)
[4] RODELL M, HOUSER P R, JAMBOR U, et al. The global land data assimilation system[J]. Bulletin of the American Meteorological Society,2004,85(3):381 − 394. doi: 10.1175/BAMS-85-3-381
[5] MIRALLES D, HOLMES T, DE JEU R, et al.Global land-surface evaporation estimated from satellite-based observations[M].Geneva: World Meteorological Orgnization, 2008.
[6] SENAY G B, BOHMS S, SINGH R K, et al. Operational evapotranspiration mapping using remote sensing and weather datasets: a new parameterization for the SSEB approach[J]. JAWRA Journal of the American Water Resources Association,2013,49(3):577 − 591. doi: 10.1111/jawr.12057
[7] MU Q Z, ZHAO M S, RUNNING S W. Improvements to a MODIS global terrestrial evapotranspiration algorithm[J]. Remote Sensing of Environment,2011,115(8):1781 − 1800. doi: 10.1016/j.rse.2011.02.019
[8] DZIKITI S, JOVANOVIC N Z, BUGAN R D, et al. Comparison of two remote sensing models for estimating evapotranspiration: algorithm evaluation and application in seasonally arid ecosystems in South Africa[J]. Journal of Arid Land,2019,11(4):495 − 512. doi: 10.1007/s40333-019-0098-2
[9] HASSAN A, ISMAIL S S, ELMOUSTAFA A, et al. Evaluating evaporation rate from high Aswan Dam Reservoir using RS and GIS techniques[J]. The Egyptian Journal of Remote Sensing and Space Science,2018,21(3):285 − 293. doi: 10.1016/j.ejrs.2017.10.001
[10] 何慧娟, 卓静, 董金芳, 等. 基于MOD16监测陕西省地表蒸散变化[J]. 干旱区地理,2015,38(5):960 − 967. [HE Huijuan, ZHUO Jing, DONG Jinfang, et al. Surveying variations of evapotranspiration in Shaanxi Province Using MOD16 products[J]. Arid Land Geography,2015,38(5):960 − 967. (in Chinese with English abstract)
[11] WOOD E F, SU H B, MCCABE M, et al. Estimating evaporation from satellite remote sensing[C]//2003 IEEE International Geoscience and Remote Sensing Symposium. Proceedings (IEEE Cat. No.03CH37477). Toulouse, France. IEEE, 2003: 1163-1165.
[12] 孙玉芳. 基于遥感监测指数模型的银川平原土壤盐渍化动态研究[J]. 地下水,2019,41(5):80 − 82. [SUN Yufang. Study on soil salinization dynamics in Yinchuan plain based on remote sensing monitoring index model[J]. Ground Water,2019,41(5):80 − 82. (in Chinese with English abstract)
[13] 袁丽华, 蒋卫国, 申文明, 等. 2000—2010年黄河流域植被覆盖的时空变化[J]. 生态学报,2013,33(24):7798 − 7806. [YUAN Lihua, JIANG Weiguo, SHEN Wenming, et al. The spatio-temporal variations of vegetation cover in the Yellow River Basin from 2000 to 2010[J]. Acta Ecologica Sinica,2013,33(24):7798 − 7806. (in Chinese with English abstract)
[14] THEIL H. A rank-invariant method of linear and polynomial regression analysis[C]//Advanced Studies in Theoretical and Applied Econometrics. Dordrecht: Springer Netherlands, 1992: 345-381.
[15] TIAN Q, WANG Q, ZHAN C, et al. Analysis of climate change in the coastal zone of Eastern China, against the background of global climate change over the last fifty years: case study of Shandong peninsula, China[J]. International Journal of Geosciences,2012,3(2):379 − 390. doi: 10.4236/ijg.2012.32042
[16] 张明明. 2000—2015年中国干旱半干旱区蒸散发时空变化及其影响因素分析[D]. 西安: 长安大学, 2019.
ZHANG Mingming. Analysis of the temporal and spatial variation of evapotranspiration and its driving factors in arid and semi-arid region of China from 2000 to 2015[D]. Xi’an: Chang’an University, 2019. (in Chinese with English abstract)
[17] 温媛媛, 赵军, 王炎强, 等. 基于MOD16的山西省地表蒸散发时空变化特征分析[J]. 地理科学进展,2020,39(2):255 − 264. [WEN Yuanyuan, ZHAO Jun, WANG Yanqiang, et al. Spatiotemporal variation characteristics of surface evapotranspiration in Shanxi Province based on MOD16[J]. Progress in Geography,2020,39(2):255 − 264. (in Chinese with English abstract) doi: 10.18306/dlkxjz.2020.02.007
[18] 王丽霞, 张珈玮, 张双成, 等. 基于CA-Markov模型的陕西省植被覆盖模拟及预测[J]. 安徽农业科学,2020,48(4):53 − 56. [WANG Lixia, ZHANG Jiawei, ZHANG Shuangcheng, et al. Simulation and prediction of vegetation coverage in Shaanxi Province based on CA-Markov model[J]. Journal of Anhui Agricultural Sciences,2020,48(4):53 − 56. (in Chinese with English abstract)
[19] BOUCHET R J. Évapotranspiration potentielle et évaporation sous abri[C]//Biometeorology. Amsterdam: Elsevier, 1962: 540-545.
[20] 邵小路, 姚凤梅, 张佳华, 等. 基于蒸散干旱指数的华北地区干旱研究[J]. 气象,2013,39(9):1154 − 1162. [SHAO Xiaolu, YAO Fengmei, ZHANG Jiahua, et al. Analysis of drought in North China based on evapotranspiration drought index[J]. Meteorological Monthly,2013,39(9):1154 − 1162. (in Chinese with English abstract)
[21] 薛阳, 金晓媚, 朱晓倩. 宁夏沿黄经济区蒸散量变化特征及水均衡方法验证[J]. 水文地质工程地质,2017,44(3):27 − 32. [Xue Yang, Jin Xiaomei, Zhu Xiaoqian. Variation of evapotranspiration of Ningxia Yellow River economic zone and the validation using water budget method[J]. Hydrogeology & Engineering Geology,2017,44(3):27 − 32. (in Chinese with English abstract)
-
图(9)
表(3)
计量
- 文章访问数: 1874
- PDF下载数: 77
- 施引文献: 0