免费文献传递   相关文献

Applicability of established drought index in Huang-Huai-Hai region based on actual evapotranspiration.

基于实际蒸散构建的干旱指数在黄淮海地区的适用性


基于NOAH陆面模式模拟的实际蒸散产品,分析了2002—2010年黄淮海地区实际蒸散的时空分布特征.同时,结合MOD17潜在蒸散数据和MOD13 NDVI构建了2002—2010年的农业干旱指数——干旱敏感性指数(DSI),并以2002年1—12月为例,利用帕默尔干旱指数(PDSI)、冬小麦减产率以及实际旱情资料,分析了DSI在该地区干旱监测的适用性.结果表明: 黄淮海地区年均实际蒸散从西北向东南递增,最高值出现在研究区域东南部(800~900 mm),最低值出现在西北部(<300 mm);DSI、PDSI两种干旱指数的年际变化呈正相关(R2=0.61)和变化趋势的一致性,均在2002年达到最低(-0.61和-1.33),2003年最高(0.81和0.92),DSI与冬小麦减产率的相关性(R2=0.43)明显优于PDSI(R2=0.06),且表征干旱的空间分辨率较高,对区域农业干旱程度的判断和旱情的指示比较可靠.

Based on the modeled products of actual evapotranspiration with NOAH land surface model, the temporal and spatial variations of actual evapotranspiration were analyzed for the Huang-Huai-Hai region in 2002-2010. In the meantime, the agricultural drought index, namely, drought severity index (DSI) was constructed, incorporated with products of MOD17 potential evapotranspiration and MOD13 NDVI. Furthermore, the applicability of established DSI in this region in the whole year of 2002 was investigated based on the Palmer drought severity index (PDSI), the yield reduction rate of winter wheat, and drought severity data. The results showed that the annual average actual evapotranspiration within the survey region increased from the northwest to the southeast, with the maximum of 800-900 mm in the southeast and the minimum less than 300 mm in the northwest. The DSI and PDSI had positive correlation (R2=0.61) and high concordance in change trend. They all got the low point (-0.61 and -1.33) in 2002 and reached the peak (0.81 and 0.92) in 2003. The correlation between DSI and yield reduction rate of winter wheat (R2=0.43) was more significant than that between PDSI and yield reduction rate of winter wheat (R2=0.06). So, the DSI reflected a high spatial resolution of drought pattern and could reflect the region agricultural drought severity and intensity more accurately.
 


全 文 :基于实际蒸散构建的干旱指数在
黄淮海地区的适用性
王  莹1  吴荣军2,3  郭照冰1∗
( 1南京信息工程大学环境与科学学院, 南京 210044; 2南京信息工程大学应用气象学院, 南京 210044; 3江苏省农业气象重点
实验室, 南京 210044)
摘  要  基于 NOAH陆面模式模拟的实际蒸散产品,分析了 2002—2010 年黄淮海地区实际
蒸散的时空分布特征.同时,结合 MOD17 潜在蒸散数据和 MOD13 NDVI 构建了 2002—2010
年的农业干旱指数———干旱敏感性指数(DSI),并以 2002 年 1—12 月为例,利用帕默尔干旱
指数(PDSI)、冬小麦减产率以及实际旱情资料,分析了 DSI在该地区干旱监测的适用性.结果
表明: 黄淮海地区年均实际蒸散从西北向东南递增,最高值出现在研究区域东南部(800~900
mm),最低值出现在西北部(<300 mm);DSI、PDSI 两种干旱指数的年际变化呈正相关(R2 =
0.61)和变化趋势的一致性,均在 2002 年达到最低( -0.61 和-1.33),2003 年最高(0.81 和
0.92),DSI与冬小麦减产率的相关性(R2 = 0.43)明显优于 PDSI(R2 = 0.06),且表征干旱的空
间分辨率较高,对区域农业干旱程度的判断和旱情的指示比较可靠.
关键词  实际蒸散; 农业干旱指数; 黄淮海地区
Applicability of established drought index in Huang⁃Huai⁃Hai region based on actual evapo⁃
transpiration. WANG Ying1, WU Rong⁃jun2,3, GUO Zhao⁃bing1∗ ( 1 School of Environmental
Science and Engineering, Nanjing University of Information Science Technology, Nanjing 210044,
China; 2School of Applied Meteorology, Nanjing University of Information Science Technology, Nan⁃
jing 210044, China; 3Jiangsu Province Key Laboratory of Agricultural Meteorology, Nanjing
210044, China) .
Abstract: Based on the modeled products of actual evapotranspiration with NOAH land surface
model, the temporal and spatial variations of actual evapotranspiration were analyzed for the Huang⁃
Huai⁃Hai region in 2002-2010. In the meantime, the agricultural drought index, namely, drought
severity index (DSI) was constructed, incorporated with products of MOD17 potential evapotranspi⁃
ration and MOD13 NDVI. Furthermore, the applicability of established DSI in this region in the
whole year of 2002 was investigated based on the Palmer drought severity index (PDSI), the yield
reduction rate of winter wheat, and drought severity data. The results showed that the annual aver⁃
age actual evapotranspiration within the survey region increased from the northwest to the southeast,
with the maximum of 800-900 mm in the southeast and the minimum less than 300 mm in the
northwest. The DSI and PDSI had positive correlation (R2 = 0.61) and high concordance in change
trend. They all got the low point ( -0.61 and -1.33) in 2002 and reached the peak (0.81 and
0.92) in 2003. The correlation between DSI and yield reduction rate of winter wheat (R2 = 0.43)
was more significant than that between PDSI and yield reduction rate of winter wheat (R2 = 0.06).
So, the DSI reflected a high spatial resolution of drought pattern and could reflect the region agricul⁃
tural drought severity and intensity more accurately.
Key words: actual evapotranspiration; agricultural drought index; Huang⁃Huai⁃Hai region.
本文由中国气象局农业气象保障与应用技术重点实验室开放基金项目(AMF201305)、国家自然科学基金项目(41575110)和国家自然科学基
金委员会⁃广东联合基金(二期)超级计算科学应用研究专项(nsfc2015_233)资助 This work was supported by the Key Laboratory of Agrometeorolo⁃
gical Safeguard and Applied Technique of China Meteorological Administration ( AMF201305), the National Natural Science Foundation of China
(41575110) and the National Natural Science Foundation of China⁃Guangdong Province (Second) Supercomputer Science Project (nsfc2015_233).
2015⁃08⁃10 Received, 2016⁃02⁃04 Accepted.
∗通讯作者 Corresponding author. E⁃mail: guocumt@ nuist.edu.cn
应 用 生 态 学 报  2016年 5月  第 27卷  第 5期                                            http: / / www.cjae.net
Chinese Journal of Applied Ecology, May 2016, 27(5): 1603-1610                  DOI: 10.13287 / j.1001-9332.201605.012
    干旱是中国的主要自然灾害之一,其中农业干
旱是以植物生长状态以及土壤含水量为特征,生长
季节内因缺少水分供给而造成土壤缺水,使农作物
生长发育受抑,导致作物减产甚至绝收的一种农业
气象灾害.在全球气候变化的大背景下,干旱发生的
频率和强度均在迅速增长.黄淮海地区是我国重要
的粮食产地之一,属于季风气候,春季少雨、夏季降
水变率大,不同时空尺度的干旱在此区域频繁发生,
对农作物的生长和产量影响较大[1-6],因此对黄淮
海地区的干旱综合评估和监测具有重要的现实意
义[7] .为了研究干旱发生和发展的规律,国内外学者
们做了大量研究[8-10] .
干旱监测技术已有多年的发展历史,传统的干
旱监测指标包括仅依赖降水的单要素指标(如降水
距平等)以及降水和温度要素结合的干旱指标
等[11],逐渐发展到后来的多因素结合的复杂干旱指
标,其对干旱的监测精度较传统的监测指标更为全
面和精确.实际蒸散过程涉及气象、土壤、植物、地
理、水利等领域,是多学科交叉的复杂过程,在地表
能量平衡以及水分循环方面起着关键作用[12-14] .准
确模拟实际蒸散能够为干旱的监测工作提供更多有
利依据,因此基于蒸散的干旱指数具有明确的物理
意义. Palmer[15]综合考虑了实际蒸散量、潜在蒸散
量、水分需求、水分供给等要素,基于水分平衡原理
建立了一个干旱指数(Palmer drought severity index,
PDSI),该指数被广泛应用于各个地区的干旱监测、
旱情比较以及旱情的时空分布特征分析[16-19];An⁃
derson等[20-21]提出了一种基于 ET / PET的蒸散发胁
迫指数(ESI),获得了较高的干旱监测精度;邵小路
等[22]基于蒸散干旱指数(EDI)对华北地区干旱进
行研究;Mu 等[23]在 ET / PET 的基础上结合归一化
植被指数(NDVI)提出了一种干旱指数 (drought se⁃
verity index, DSI),并应用于全球尺度干旱监测,与
PDSI、植被净初级生产力(NPP)等对比得到了良好
的效果.DSI 综合考虑干旱过程中农作物的水分胁
迫信息和生长状况,在农业干旱监测方面具有较大
的应用潜力.
本研究以黄淮海地区为研究区域,采用 NOAH
陆面模式模拟 2002—2010年研究区域实际蒸散,并
利用地面实测数据对模拟结果进行验证,然后结合
潜在蒸散和归一化植被指数构建 DSI,通过与 PDSI、
冬小麦减产率以及实际旱情等资料对比分析,开展
DSI在研究区干旱监测的适用性分析.
1  研究地区与研究方法
1􀆰 1  研究区概况
黄淮海地区位于我国东部(32°—40° N,114°—
121° E),北起长城,南至桐柏山、大别山北麓,西倚
太行山和豫西伏牛山地,东濒渤海和黄海,其主体为
由黄河、淮河与海河及其支流冲积而成的黄淮海平
原,地势平坦、利于耕作.全区土地总面积 46.95×104
km2,总耕地面积占全国耕地面积的 25%,是我国耕
地面积最多的农业区,冬小麦和夏玉米轮作为本区
典型农业种植模式.黄淮海地区属于暖温带季风气
候,降水季节差异明显且变率较大,导致干旱频发.
目前干旱已经成为影响本区域农作物生长发育最主
要的农业灾害,因此选择黄淮海地区进行农业干旱
监测具有一定的典型性.
1􀆰 2  数据来源
1􀆰 2􀆰 1实际蒸散  NOAH陆面模式(NOAH land sur⁃
face model, LSM)被广泛用于陆面过程的综合模
拟[24-26] .在模式中,实际蒸散量(ET)的计算主要包
括 3大部分,即表层土壤直接蒸散(Es)、植被截留
蒸散量(Ec)和植被的蒸腾(Ev).其中,Es被认为与
潜在蒸散量间存在较好的线性关系,潜在蒸散量采
用彭曼公式计算;Ec通过计算地表覆盖率、降水量与
截留系数的乘积获得;Ev采用传统的阻抗法计算.
LSM的主要输入数据包括:气象强迫数据(GL⁃
DAS)、逐月 MODIS⁃BRDF(MCD43C3)地表反照率、
SURFRAD(U. S. Surface Radiation Budget Network)
辐射观测值、UMD⁃AVHRR 地表类型、逐周 AVHRR
植被覆盖率等,将输入数据空间分辨率统一后,调试
并运行模式得到中国区域 ET 模拟结果,然后进行
裁剪、累加、双线性插值,得到 2002—2010 年 0.1°×
0.1°月值以及年值研究区域数据集.
1􀆰 2􀆰 2实测数据   本研究采用长白山、千烟洲和禹
城站点的通量数据对 LSM模型模拟的 ET产品进行
验证.这些站点的气候、海拔和地表覆盖类型存在显
著差异.长白山站点的土地覆盖类型为温带混交林,
千烟洲站点为常绿针叶林,禹城站为农田(表 1).其
中,长白山和千烟洲站点实测数据源于文献[27];
禹城站点 2003、2004 年蒸散量数据来自源于文献
[28],2008年的日蒸散数据来自禹城站 CERN监测
项目农田生态系统农田蒸散量表(水量平衡法),剔
除异常值后与 LSM日蒸散值进行对比分析.
1􀆰 2􀆰 3潜在蒸散   本文所用潜在蒸散 ( PET)为
1 km×1 km空间分辨率的 MODIS PET月合成产品,
4061 应  用  生  态  学  报                                      27卷
表 1  观测站点信息
Table 1  Information of observed sites
站点
Site
坐标
Coordinate
土地类型
Land type
年份
Year
实测数据来源
Measured data source
长白山
CBS
42°24′00″ N,128°28′00″ E 温带混交林
Temperate mixed forest
2003、2004 [27]
千烟洲
QYZ
26°43′48″ N,115°42′00″ E 常绿针叶林
Evergreen coniferous forest
2003、2004 [27]
禹城
YC
36°49′45″ N,116°34′14″ E 农田
Cropland
2003、2004、
2008
[28]
CERN监测项目 CERN Program
CBS: Changbai Mountains; QYZ: Qianyinzhou; YC: Yucheng. 下同 The same below.
即 MOD16A2.研究区域共有 5 景数据覆盖,经过拼
接、双线性插值、裁剪、定标、提取,得到 2002—2010
年黄淮海区域 0.1°×0.1°PET 的长时间序列数据集.
1􀆰 2􀆰 4归一化植被指数  本文中的归一化植被指数
(NDVI)数据为 1 km×1 km 空间分辨率的 MODIS
NDVI 月合成产品,即 MOD13A3.统一时空分辨率
后,进行拼接、双线性插值、裁剪、定标和提取,得到
研究区 2002—2010年 0.1° ×0.1°NDVI 长时间序列
数据集.
1􀆰 2􀆰 5辅助数据   1)中国农业灾情数据.研究区
2002年 5—9 月发生持续旱灾的 12 个农业站点的
灾情数据源于中国气象共享网 ( http: / / cdc. nmic.
cn).
2)1970—2010年研究范围内 29 个县级(包括
河间市、阜城县、隆尧县、香河县、曲周县、定兴县、栾
城县、济阳县、青州市、兖州市、邹平县、陵县、茌平
县、阳谷县、滑县、濮阳县、息县、太康县、正阳县、西
平县、颍上县、涡阳县、谯城区、凤台县、兴化市、东海
县、泗洪县、高邮市)冬小麦产量数据,来自中国农
业部种植业管理司县级农作物数据库.本研究采用
滑动平均法对站点 1970—2010 年的实际产量进行
5 年滑动平均求出趋势产量,然后根据实际产量和
趋势产量求出相对气象产量.
3)帕默尔指数.帕尔默指数(PDSI)综合了水分
亏缺和持续时间因子对干旱程度的影响,且考虑了
前期天气状况,具有较好的时空可比性[29-30],在美
国的干旱事件监测、分析以及序列重建上被广泛应
用.帕默尔干旱等级划分如表 2所示.本文对 MOD17
的 0.5°×0.5°空间分辨率的 PDSI月产品进行提取与
裁剪,建立 2002—2010年 PDSI月数据集.
1􀆰 3  研究方法
基于 ET、PET和 NDVI,Mu等[23]构建了干旱敏
感度指数(drought sensitive index,DSI),用于全球尺
度的干旱监测,并对干旱等级进行划分(表 3).DSI
算法如下:
    R=ET / PET (1)
ZR =(R-RAVE) / RSTU (2)
ZNDVI = (NDVI-NDVIAVE) / NDVISTU (3)
Z=ZR+ZNDVI (4)
DSI = (Z-ZAVE) / ZSTU (5)
式中:ET是实际蒸散;PET为潜在蒸散;RAVE是 R 的
多年均值;RSTU是 R 的标准差;ZNDVI是 NDVI 标准化
的值;NDVIAVE是 NDVI的多年均值;NDVISTU是NDVI
的标准差;ZAVE是 Z 的多年均值;ZSTU是 Z 的标
准差.
2  结果与分析
2􀆰 1  实际蒸散模拟结果的验证
从表 4可以看出,模拟年 ET值与实测年 ET 值
较为接近,模拟值和实测值全年相对偏差在
-11.7%~5.7%.由于实测值是基于站点测量,与模
型模拟尺度存在差异,并且容易受到下垫面不均匀
以及空间一致性的影响,导致模拟结果与实测结果
之间存在一定差异.
表 2  帕默尔干旱等级划分标准
Table 2  Drought grade determined by Palmer index
帕默尔指数 PDSI 等级 Level     
0.99< PDSI <-0.99 正常 Normal
-1.99< PDSI <-1 轻微干旱 Incipient drought
-2.99< PDSI <-2 中等干旱 Moderate drought
-3.99< PDSI <-3 严重干旱 Severe drought
PDSI≤-4 极端干旱 Extreme drought
表 3  DSI干旱等级划分
Table 3  Drought grade determined by DSI
干旱敏感度指数 DSI 等级 Level
-0.29< DSI <0.29 无旱 Normal
-0.59< DSI <-0.29 轻旱 Incipient drought
-0.89< DSI <-0.59 中旱 Medium drought
-1.19< DSI <-0.9 中等偏重旱 Moderate drought
-1.49< DSI <-1.2 重旱 Severe drought
DSI≤-1.5 特旱 Extreme drought
50615期                        王  莹等: 基于实际蒸散构建的干旱指数在黄淮海地区的适用性         
表 4  LSM模拟 ET与观测数据的对比
Table 4   Comparison of the ET simulated by LSM with
observed ET
站点
Site

Year
实测值
Observed value
(mm·a-1)
模拟值
Simulated value
(mm·a-1)
相对偏差
RE
(%)
长白山 2003 536 534.8 -0.2
CBS 2004 523 506.8 -3.1
千烟洲 2003 708 720.0 1.7
QYZ 2004 769 679.2 -11.7
禹城 2003 541.7 572.8 5.7
YC 2004 581.9 557.8 -4.1
    由于数据获取的有限性,并考虑缺测和异常值
较少的情况,本文选择禹城观测站 2008年日蒸散数
据与 LSM模拟值进行对比(图 1),模拟值与实测值
相关系数为 0. 67,均方根误差为 0. 82.春季(3—5
月),禹城模拟蒸散值基本小于测量值,但其他月份
的蒸散值变化量与实测值基本一致,这是由于禹城
地下水位偏高,且半干旱地区在春季易受到近地层
水平平流输送的影响[31],导致实测数据偏高.此外,
模拟值与实测值的季节变化趋势一致,6 月有一个
较窄的低值区,这是因为此时冬小麦收获,玉米刚刚
播种,所以植被蒸腾较低,导致蒸散值出现低谷.说
明 LSM对蒸散的模拟效果较好.
    由图 2 可以看出,年均 ET 最高值出现在研究
区域南部(800~900 mm),最低值出现在西北(<300
mm).ET在黄淮海地区从西北到东南,阶梯状分布
较明显,呈逐渐增加的趋势.研究区域内,河北北部
以及京津部分地区的年均 ET值为 300~400 mm;河
北南部和山东北部由北向南递增,达到 400 ~ 500
mm;山东南部、江苏北部以及安徽北部的年均 ET
值从 500 mm 递增到 800 mm;江苏、安徽南部部分
地区因降水量高且植被覆盖度高,达到研究区域年
图 1  禹城观测站点 2008 年实测 ET(Ⅰ)与 LSM 模型模拟
结果(Ⅱ)对比
Fig.1   Comparison between observed ET (Ⅰ) and simulated
ET (Ⅱ) at the Yucheng Station in 2008.
图 2  2002—2010年黄淮海地区年均 ET的空间分布
Fig.2   Spatial distribution of annual mean values of ET in
Huang⁃Huai⁃Hai region in 2002-2010.
均 ET最高值,为 900 mm左右.
2􀆰 2  干旱指数的时空变化特征
2002—2010年,研究区 DSI 和 PDSI 的变化趋
势一致,两者的相关系数为 0.61,最低值均出现在
2002年,分别为- 0. 61 和- 1. 33,最高值均出现在
2003年,分别为 0.81和 0.92(图 3).
由于 DSI、PDSI均在 2002 年出现低值,说明相
较于其他年份研究区域内旱情发生概率最大. 2002
年,黄淮海地区降水较常年减少,旱灾频发,主要集
中在 1—4月和 7—12 月,为干旱偏重的一年,相较
于其他年份无旱比例也最小,为了探讨 DSI 对干旱
监测的适用性,下面以旱情显著的 2002 年为例,对
研究区域内干旱变化过程进行实例分析.
由图 4 可以看出,1—3 月,河北、山东、江苏和
安徽等地均出现大面积干旱,直到 6月,干旱空间分
布范围缩小,只有河北和山东的小范围还存在轻微
干旱.MODIS17 PDSI 产品在 1—3 月也出现了大范
围旱情,不仅河北、京津冀,以及江苏和安徽北部出
现较重干旱,河南和山西也出现大范围干旱,且山东
只显示在西部出现中等干旱,其余地区都是轻微
图 3  2002—2010年 DSI和 PDSI指数的年际变化
Fig.3  Annual variation of DSI and PDSI from 2002 to 2010.
6061 应  用  生  态  学  报                                      27卷
图 4  2002年黄淮海地区 DSI和 PDSI的空间分布
Fig.4  Spatial distribution of DSI and PDSI in Huang⁃Huai⁃Hai region in 2002.
a~ l: 1—12月 January to December.
干旱.与 DSI一样,直到 6 月降雨的增多,PDSI 才显
示旱情有所缓解,干旱面积大幅减小,干旱主要分布
在京津冀和河北部分地区以及山东西部.这是由于
2002年 1—3月,华北大部、黄淮北部等地降水量较
往年偏少一半以上,同时气温持续偏高,土壤蒸发量
较高,导致春旱发生早且发展迅速;随着 6月下旬的
大范围降雨,上述旱区的部分旱情得到缓解.从 7 月
开始由于降水逐渐减少,两种干旱指数均显示旱情
急剧发展,DSI 旱情由京津冀地区向山东地区以及
江苏和安徽北部开始持续蔓延开来,尤其 9—11 月,
旱情较重干旱面积较大.
    图 5中相对气象产量<0的站点发生减产,而相
对气象产量负值的绝对值即为该站点的减产率.可
以看出,在山东省内,邹平、青州、莱阳 3个站点减产
率明显高于其他站点,分别为 9.5%、17.8%和 5.9%.
在主要影响冬小麦生长的 1—5 月,上述 3 站点的
PDSI结果显示 3 月有轻微干旱,而其余月份无干
旱 ,与减产率不一致;DSI的结果显示,这3个站点
70615期                        王  莹等: 基于实际蒸散构建的干旱指数在黄淮海地区的适用性         
图 5  2002年县级农业站点冬小麦相对气象产量分布
Fig.5  Distribution of the relative meteorological yield of winter
wheat in 2002 at county agro⁃meteorological station in the study
area (%).
除 5月因为降雨干旱得到缓解外,其余月份均发生
严重干旱,干旱覆盖范围与站点减产率分布基本相
符.江苏省内 5 个农业监测站点(东海、泗洪、高邮、
兴化、姜堰)的减产率分别为 2. 0%、3. 6%、9. 0%、
4.7%和 2.5%,PDSI结果显示 1—5月上述 5 个站点
均未发生干旱,DSI结果显示 1—5月 5 个站点所在
区域均发生干旱,DSI 结果与减产率分布更为符合.
河北省内农业站点冬小麦产量均发生减产,两种指
数显示在相应区域同时发生大面积干旱.河南东部
站点减产率为 14.0%,DSI 产品显示 1—5月除 2 月
有干旱外,其他月份均未发生明显干旱,PDSI 产品
则与站点减产率分布更为相符.
由图 6可以看出,DSI 与减产的 19 个站的减产
率的决定系数(R2 )为 0. 43,明显优于 PDSI.由于
PDSI产品为全球产品,容易存在空间异质性;另一
方面,DSI在计算过程中考虑作物生长因素而 PDSI
主要考虑气象因素,导致 PDSI 的旱情指示与减产
率结果存在差异性.
图 6  干旱指数 DSI和 PDSI与冬小麦减产率散点图
Fig.6   Scatterplots of DSI and PDSI with the yield reduction
rate of winter wheat.
    通过灾害类型、经纬度、灾害发生时间等条件的
筛选,选出表 5 中的 11 个点作为实际灾情对比对
象.表中受害百分比代表农业站点所在区域受干旱
影响的作物面积与总的耕种区域面积的比值,受害
程度代表由干旱引起的作物损害程度.临汾站受旱
比例小且旱灾对冬小麦影响较轻,此站的 DSI 为
-0.4(轻度干旱),而 PDSI为-2.5(中等干旱),这是
由于该站 5月降水减少、温度升高,容易发生气象干
旱,但可能由于灌溉等原因该地区农作物并没有发
生重大农业干旱,PDSI更倾向于考虑气候环境因素
而并未考虑作物生长状况,所以旱情指示偏重.宿州
站受旱比例大但旱灾程度较轻,DSI 为-0.5(轻度干
旱 ) ,PDSI为-2.9(中等干旱) ,因为PDSI产品分辨
表 5  2002年黄淮海地区 11个农业气象站点的干旱灾情数据
Table 5  Drought observation data of 11 agro⁃meteorological stations in Huang⁃Huai⁃Hai region in 2002
站点
Site
经度
Longitude
(°)
纬度
Latitude
(°)
月份
Month
干旱百分比
Drought
percentage
干旱程度 
Drought 
intensity 
受灾作物
Affected
crop
DSI PDSI
黄骅 Huanghua 117.35 38.37 5 90~100 重 Severe 冬小麦 Winter wheat 重 Severe 重 Severe
临汾 Linfen 111.50 36.07 5 30~39 轻 Mild 冬小麦 Winter wheat 轻 Mild 中 Moderate
莱阳 Laiyang 120.70 36.93 6 40~49 中 Moderate 冬小麦 Winter wheat 中 Moderate 轻 Mild
南宫 Nangong 115.38 37.37 7 40~49 重 Severe 其他作物 Other crops 中 Moderate 中 Moderate
宿州 Suzhou 116.98 33.63 8 90~100 轻 Mild 大豆 Soybean 轻 Mild 中 Moderate
怀来 Huailai 115.50 40.40 8 0~9 中 Moderate 其他作物 Other crops 中 Moderate 中 Moderate
南宫 Nangong 115.38 37.37 9 60~69 重 Severe 所有作物 All crops 重 Severe 重 Severe
商丘 Shangqiu 115.67 34.45 9 30~39 中 Moderate 夏玉米 Summer corn 中 Moderate 中 Moderate
赣榆 Ganyu 119.12 34.83 9 60~69 不明 Unknown 花生 Peanut 重 Severe 重 Severe
惠民 Huimin 117.53 37.50 9 70~79 中 Moderate 所有作物 All crops 重 Severe 重 Severe
潍坊 Weifang 119.08 36.70 9 60~69 重 Severe 所有作物 All crops 重 Severe 重 Severe
8061 应  用  生  态  学  报                                      27卷
率较低,容易存在空间异质性,所以旱情指示产生偏
差.整体来讲,站点受灾程度与 DSI 和 PDSI 的干旱
等级大部分一致,个别站点 DSI 的旱情等级指示更
为准确.
3  讨    论
本文通过 LSM 模型计算了 2002—2010 年黄淮
海地区 ET的时空变化,与前人的模拟结果基本一
致[27-28,32] .其中,黄淮南部阜阳地区年蒸散量大于
900 mm,河北南部以及河南的新乡、濮阳、郑州等中
西部地区的年蒸散量少于 700 mm,与王菱等[33]以
田间试验资料为基础计算黄淮海地区在自然条件下
农田蒸散量的研究结果相符.黄淮海地区基本为冬
小麦、夏玉米一年两熟,本文结果的农田站 ET 与作
物生长进程基本一致,呈现双峰型变化,春季作物进
入生长旺季,ET 逐渐增加,5 月上旬达到高峰期,6
月小麦成熟收割,ET达到波谷,直到 8 月上旬,玉米
生长旺季 ET 达到第 2 个高峰,与邵小路等[22]对华
北地区干旱研究中栾城农田站的实际蒸散变化趋势
一致.
在气候变暖的大背景下,农业干旱所带来的一
系列影响成为焦点.DSI 作为一种新型的农业干旱
指数,重点考虑了 NDVI 与 ET / PET 两个基本要素,
即结合了农作物的生长状态和作物的水分胁迫状
况,通过两者综合反映干旱状况.由于 NDVI 在春季
和冬季植被裸露时易产生异常值[23],且存在一定的
时间滞后性,客观程度上影响了干旱监测的准确性.
而土壤湿度体现了地表降水和蒸发的综合效应,是
监测干旱的关键指标之一[34],所以在进一步研究中
将土壤湿度加入计算中将更准确可靠.本研究构建
的空间分辨率为 0.1°×0.1°的 DSI 干旱指数,进一步
提高了区域尺度监测精度,解决了黄健熙等[35]在计
算冬小麦产量与 DSI相关分析中提出的空间精度问
题,提高了区域干旱监测精度,对干旱程度的判断和
旱情的指示更为准确.计算 DSI 时,数据时间尺度的
变化对计算结果的影响有待深入研究.
4  结    论
利用陆面模式 LSM 模拟的实际蒸散数据与禹
城等站点实测值的相关性较好,说明将该数据用于
研究区域的干旱监测模拟是可靠的.研究区域内 ET
空间分布呈现东南向西北方向降低的明显阶梯趋
势,山东西北部为 300 ~ 400 mm·a-1,山东地区为
400~600 mm·a-1,长黄淮南部地区为 900 mm·a-1
左右.2002—2010年,DSI 和 PDSI 年际变化趋势基
本一致,且 DSI的农业干旱指示相较于 PDSI更为准
确,在黄淮海平原农业干旱监测应用中具有一定的
适用性.本文基于陆面模式 LSM 模拟的 ET,结合
NDVI和 PET构建了干旱指数 DSI,将分辨率提高到
0.1°×0.1°,较 MODIS17 PDSI 干旱监测产品进一步
提高了干旱区域监测精度.
参考文献
[1]  An S⁃Q (安顺清), Zhu Z⁃X (朱自玺), Wu N⁃Y (吴
乃元), et al. Studied results of CROP water stress and
drought in central region of Huang⁃Huai⁃Hai Plain.
Scientia Agricultura Sinica (中国农业科学), 1991, 24
(2): 13-18 (in Chinese)
[2]  Xu H⁃Y (徐焕颖), Jia J⁃H (贾建华), Liu L⁃Y (刘
良云), et al. Drought monitoring in Huang⁃Huai⁃Hai
Plain using the multi⁃drought indices. Remote Sensing
Technology and Application (遥感技术与应用), 2015,
30(1): 25-32 (in Chinese)
[3]  Yin H (尹  晗), Li Y⁃H (李耀辉). Summary of ad⁃
vance on drought study in southwest China. Journal of
Arid Meteorology (干旱气象), 2013, 31(1): 182-
193 (in Chinese)
[4]  Liu X⁃Y (刘晓云), Li D⁃L (李栋梁), Wang J⁃S (王
劲松). Spatiotemporal characteristics of drought over
China during 1961 - 2009. Journal of Desert Research
(中国沙漠), 2012, 32(2): 473-483 (in Chinese)
[5]  He J⁃Y (贺晋云), Zhang M⁃J (张明军), Wang P
(王   鹏), et al. Climate characters of the extreme
drought events in southwest China during recent 50
years. Acta Geographica Sinica (地理学报), 2011, 66
(9): 1179-1190 (in Chinese)
[6]  Zhang W⁃C (张万诚), Zheng J⁃M (郑建萌), Ren J⁃Z
(任菊章). Climate characteristics of extreme drought
events in Yunnan. Journal of Catastrophology (灾害
学), 2013, 28(1): 59-64 (in Chinese)
[7]  Xu J⁃W (徐建文). Potential Effects of Drought on Win⁃
ter Wheat Yield Based on DSSAT Model Over Huang⁃
Huai⁃Hai Plain. Master Thesis. Beijing: Chinese Acade⁃
my of Agricultural Sciences, 2014 (in Chinese)
[8]  Hou Y⁃Y (侯英雨), He Y⁃B (何延波), Liu Q⁃H (柳
钦火), et al. Research progress on drought index. Chi⁃
nese Journal of Ecology (生态学杂志), 2007, 26(6):
892-897 (in Chinese)
[9]  Fei Z⁃Y (费振宇), Sun H⁃W (孙宏巍), Jin J⁃L (金
菊良), et al. Temporal and spatial patterns of meteoro⁃
logical drought hazard in China for recent 50 years. Wa⁃
ter Resources and Power (水电能源科学), 2014, 32
(12): 5-10 (in Chinese)
[10]  Chen F⁃H (陈发虎), Huang W (黄   伟), Jin L⁃Y
(靳立亚), et al. Spatiotemporal precipitation variations
in the arid central Asia in the context of global warming.
Science China Earth Sciences (中国科学: 地球科学),
2011, 41(11): 1647-1657 (in Chinese)
[11]  Zhang Q (张  强), Zhang L (张  良), Cui X⁃C (崔
显成), et al. Progresses and challenges in drought as⁃
sessment and monitoring. Advances in Earth Science (地
球科学进展), 2011, 26(7): 763-778 (in Chinese)
[12]  Wu R⁃J (吴荣军), Ge Q (葛  琴), Zhan X⁃W (詹
90615期                        王  莹等: 基于实际蒸散构建的干旱指数在黄淮海地区的适用性         
习武), et al. Research of land surface evapotranspira⁃
tion based on remote sensing process coupling model
BEPS. Ecology and Environmental Sciences (生态环境
学报), 2013, 22(6): 1001-1008 (in Chinese)
[13]  Guadalupe SC, Jurgen V, Alirio A, et al. Assessment of
the Eumetsat LSA⁃SAF evapotranspiration product for
drought monitoring in Europe. International Journal of
Applied Earth Observations and Geoinformation, 2014,
30: 190-202
[14]  Jin XM, Guo RH, Xia W. Distribution of actual evapo⁃
transpiration over Qaidam basin, an arid area in China.
Remote Sensing, 2013, 5: 6976-6996
[15]  Palmer WC. Meteorological Drought. Washington DC:
US Weather Bureau, 1965
[16]  Zhou Y⁃W (周跃武), Feng J⁃Y (冯建英). The
drought index in the 20th century. Arid Meteorology (干
旱气象), 2006, 24(1): 79-89 (in Chinese)
[17]  Yang X⁃L (杨小利), Yang X⁃G (杨兴国), Ma P⁃L
(马鹏里), et al. Modification and applications of PDSI
for mid and eastern Gansu Province. Advances in Earth
Science (地球科学进展), 2005, 20(9): 1022-1028
(in Chinese)
[18]  Yan G⁃X (闫桂霞), Lu G⁃H (陆桂华), Wu Z⁃Y (吴
志勇), et al. Study on integrated meteorological drought
index based on PDSI and SPI. Water Resources and Hy⁃
dropower Engineering (水利水电技术), 2009, 40
(4): 10-13 (in Chinese)
[19]  Zhang Y (张   永), Chen F⁃H (陈发虎), Gou X⁃H
(勾晓华), et al. The temporal and spatial distribution
of seasonal dry⁃wet changes over the northwestern China
based on PDSI. Acta Geographica Sinica (地理学报),
2007, 62(11): 1142-1152 (in Chinese)
[20]  Anderson MC, Hain C, Wardlow B, et al. Evaluation of
drought indices based on thermal remote sensing of eva⁃
potranspiration over the Continental United States. Jour⁃
nal of Climate, 2011, 248: 2025-2044
[21]  Anderson MC, Hain C, Wardlow B, et al. A climatolo⁃
gical study of evapotranspiration and moisture stress
across the continental U. S. based on thermal remote
sensing: Surface moisture climatology. Journal of Geo⁃
physical Research, 2007, 112: 1029-1039
[22]  Shao X⁃L (邵小路), Yao F⁃M (姚凤梅), Zhang J⁃H
(张佳华), et al. Analysis of drought in north China
based on evapotranspiration drought index. Acta Meteoro⁃
logica Sinica (气象学报), 2013, 39(9): 1154-1162
(in Chinese)
[23]  Mu QZ, Zhao MS, Kimball JS, et al. A remotely sensed
global terrestrial drought severity index. Bulletin of the
American Meteorological Society, 2013, 213: 83-98
[24]  Chen F, Mitchell KE, Schaake J, et al. Modeling of
landsurface evaporation by four schemes and comparison
with FIFE observations. Journal of Geophysical Research,
1996, 101: 7251-7268
[25]  Chen F, Janic Z, Mitchell KE. Impact of atmospheric
surface layer parameterizations in the new land⁃surface
scheme of the NCEP mesoscale Eta model. Bound Layer
Meteorology, 1997, 85: 391-421
[26]  Chen F, Mitchell KE. Using the GEWEX / ISLSCP for⁃
cing data to simulate global soil moisture fields and
hydrological cycle for 1987 - 1988. Journal of Meteoro⁃
logical Society of Japan, 1999, 77: 167-182
[27]  Zhang F⁃M (张方敏), Ju W⁃M (居为民), Chen J⁃M
(陈镜明), et al. Study on evapotranspiration in East
Asia using the BEPS ecological model. Journal of Natu⁃
ral Resources (自然资源学报), 2010, 25(9): 1596-
1606 (in Chinese)
[28]  Liu Y, Zhou W, Ju J, et al. Changes of evapotranspira⁃
tion and water yield in China’ s terrestrial ecosystems
during the period from 2000 to 2010. Hydrology and
Earth System Sciences, 2013, 10: 5397-5456
[29]  Wei J (卫   捷), Ma Z⁃G (马柱国). Comparison of
palmer drought severity index, percentage of precipita⁃
tion anomaly and surface humid index. Acta Geographica
Sinica (地理学报), 2003, 58(suppl.1): 117-124 (in
Chinese)
[30]  Shen L⁃F (沈良芳), Gao R⁃B (高荣斌), Lyu J (吕
军), et al. Research on Palmer drought index in Jiangsu
Province. Journal of Anhui Agricultural Sciences (安徽农
业科学), 2007, 35(28): 8768-8788 (in Chinese)
[31]  Zhang Q (张  强), Zhang Z⁃X (张之贤), Wen X⁃M
(问晓梅), et al. Comparisons of observational methods
of land surface evapotranspiration and their influence
factors. Advances in Earth Science (地球科学进展),
2011, 26(5): 538-547 (in Chinese)
[32]  Zhou L (周   蕾), Wang S⁃Q (王绍强), Chen J⁃M
(陈镜明), et al. The spatial⁃temporal characteristics of
evapotranspiration of China’s terrestrial ecosystems
during 1991 - 2000. Resources Science (资源科学),
2009, 31(6): 962-972 (in Chinese)
[33]  Wang L (王  菱), Ni J⁃H (倪建华). Research of ac⁃
tual evaporanspiration in farmland: A case study in
Huang⁃Huai⁃Hai river region. Acta Meteorologica Sinica
(气象学报), 2001, 59(6): 785-793 (in Chinese)
[34]   Guo W⁃D (郭维栋), Ma Z⁃G (马柱国), Yao Y⁃H
(姚永红). Regional characteristics of soil moisture evo⁃
lution in northern China over recent 50 years. Acta Geo⁃
graphica Sinica (地理学报), 2003, 58(suppl.1): 84-
88 (in Chinese)
[35]  Huang J⁃X (黄健熙), Zhang J (张   洁), Liu J⁃M
(刘峻明), et al. Correlation analysis between drought
and winter wheat yields based on remotely sensed
drought severity index. Transactions of the Chinese Socie⁃
ty for Agricultural Machinery (农业机械学报), 2015,
46(3): 166-173 (in Chinese)
作者简介  王  莹,女,1991年生,硕士研究生. 主要从事农
业气象灾害研究. E⁃mail: 494900718@ qq.com
责任编辑  杨  弘
王莹, 吴荣军, 郭照冰. 基于实际蒸散构建的干旱指数在黄淮海地区的适用性. 应用生态学报, 2016, 27(5): 1603-1610
Wang Y, Wu R⁃J, Guo Z⁃B. Applicability of established drought index in Huang⁃Huai⁃Hai region based on actual evapotranspiration.
Chinese Journal of Applied Ecology, 2016, 27(5): 1603-1610 (in Chinese)
0161 应  用  生  态  学  报                                      27卷