基于标准化降水指数的鲁西地区干旱时空特征分析
2020-09-21
来源:易榕旅网
山东农业科学2017,49(7):131~138 DOI:10.14083/j.issn.1001—4942.2017.07.029 Shandong Agricultural Sciences 基于标准化降水指数的鲁西地区 干旱时空特征分析 李楠。,李慧 ,李强 ,杨志勇 ,韩学蕾 (1.聊城市气象局,山东聊城252000;2.南通市气象局,江苏南通226018) 摘要:利用鲁西地区1962--2015年8个国家气象观测站逐月的降水资料,采用标准化降水指数(SPI)、主 成分分析和ArcGIS空问插值方法,分析了该地区的年度和季节干旱演变特征。结果表明:(1)鲁西地区的干 旱发生较为频繁且持续时间较长,年度和季节干旱发生频率均为30%左右,主要以轻旱和中旱为主,重旱和 特旱发生频率较低;年度和季节的干旱站次比中,春旱和冬旱呈减少趋势,夏旱和秋旱呈增加趋势。(2)主成 分分析中前两个特征向量的空间分布与干旱频率站次比的区域分布较为一致。 关键词:鲁西地区;干旱;标准化降水指数;时空特征 中图分类号:¥161.6 文献标识号:A 文章编号:1001—4942(2017)07-0131-08 Analysis of Spatial and Temporal Characteristics of Drought in Western Shandong Based on Standardized Precipitation Index Li Nan ,Li Hui ,Li Qiang ,Yang Zhiyong ,Han Xuelei (1.Liaocheng Meteorological Bureau of Shandong Province,Liaocheng 252000,China; 2.Nantong Meteorological of Jiangsu Province,Nantong 226018,China) Abstract The annual and seasonal drought evolution characteristics in Western Shandong were analyzed based on the monthly precipitation data of 8 national meteorological stations from 1962 to 2015 by the standard- ized precipitation index(SPI),principal component analysis and spatial interpolation method of ArcGIS.The results showed that in Western Shandong,the drought had rather high occurrence frequency and lasted a long time.The annual and seasonal drought frequency fluctuated around 30%.The drought degree usually was light drought and moderate drought,and the heavy drought and excessive drought had lower occurrence frequency. In the station number ratio of annual and seasonal drought,the drought in spring and winter showed a decline trend,and the drought in summer and autumn showed an increase trend.The spatial distribution of the first two feature vectors in principal component analysis was in accordance with the regional distribution of station number ratio of drought frequency. Keywords characterist cs Western Shandong;Drought;Standardized precipitation index;Temporal and spatial 人类活动对旱涝灾害的影响十分显著 ]。全 时空分布不均匀,干旱灾害频繁发生,极端事件偏 球气候变暖,大气中能量分布发生变化,导致降雨 收稿日期:2017一O1—07 多[2 3,严重威胁着人民的生命财产安全 。而城市 基金项目:山东省气象局青年基金项目(2016SDQN19);山东省气象局面上课题(2014sdqxm03) 作者简介:李楠(1985一),女,工程师,硕士,研究方向为气象信息技术。E—mail:569386124@qq.con l 32 山东农业科学 第49卷 化的加剧进一步加重了于旱灾害的发生和程度 。 假设某一时段的降水量为X,则概率函数为: f(x):—— x ~e一 0. 聊城地处中原腹地、山东省西部,地形单一,以平原 为主,故有鲁西平原之称 。该地区自然灾害时有 发生,其中旱灾危害面积最大且频繁。目前,常用 、 p F( ) ’ (1) X>0,p>0, >0 。 的=f 早评价指标有降水距平、z指数、标准化降水 指数、帕默尔干旱指数等 ,以标准化降水指数 (SPI)方法评价鲁两地区干旱情况的研究尚未见报 道。SP1是采用r函数的标准化降水累积分布来 式中:13和^y分别为尺度和形状参数,可用极大似 然估计方法求得: 一 1+ 1+4A/3 —— — 一; ,1、 z B= ; (3) 描述降水量的变化,该指数是中国气象局推荐的气 象干旱分析指标,以月为评价尺度。本文应用该指 数从气象干旱的角度分析了鲁西地区年及季节性 }二旱发牛频率和强度的时空演变规律,明确了鲁西 地区干旱发生的时空变化规律,以期为应对气候变 化和抗旱减灾提供理论依据 。 A=lg 一音。 1 g xl。 式中:xi为降水量, 为降水量气候平均值。 (4) 概率密度函数中的参数确定后,对于某一年 的降水量x (x。≠0),可求出随机变量x小于x , 事件的概率为: F(X<X ,):』0 ̄f(X)dx。 (5) (6) l研究区域概况 聊城市地处山东省西部(鲁西),北纬35。47 ~ 37。Cl3 ,东经115。16 一116。30 。共8个县(市、区)。 处于暖温带季风气候区,属半干旱大陆性气候,光热 当降水量为0时的概率为: F(X=0)=m/n 。 式中:n 是降水量为0的样本数;n是总样本数。 将(5)、(6)式求得的概率值分别代入标准化 正态分布函数,即得: 资源充足,但降水缺乏,且时空分布不均。 F(x<x0)={ e。 !dx。(7)  ̄/2,n- 2数据来源与研究方法 2.1 数据来源 利用数值积分对(7)式进行近似求解可得: 选取1962 0l5年聊城市8个气象观测站 s 嵩 蒜 。 ㈩ r—— 点(见图1)经质控后的逐月降水资料。 式中:t √ n 1,F为(5)式或(6)式求得的概率; 在公式(8)中当F>0.5时,F=I.0一F,S=1;当 F≤0.5时,S=一1;其中:c【]=2.515517,cl= 0.802853, 2=0.010328,dI=I_432788,d2= 0.189269,d =0.001308。由(8)式求得z值,即 标准化降水指数SPI。 主成分分析及ArcGIS中的反距离加权插值 分析详见参考文献[1 1]、[12]。 本文主要分析不同时间尺度——月(SPII)、季 图1 鲁西地区气象站点分布 度(SPI3)、年(SPI12)的降水量,并根据《GB/ T20481—20HD6气象干旱等级》划分干旱等 级:一1.0<SPI≤一0.5,轻旱;一1.5<SPI≤一1.0 中旱;一2.0<SPI≤一1.5;重旱;SPI≤一2.0,特旱。 2.3干旱评价指标 2.2标准化降水指数(SPI)及其计算方法 标准化降水指数(简称SPI)是基于降水量的 分布,一般不是正态分布,而是一种偏态分布 ” 。 (1)干旱频率(Pi):某站发生干旱的频繁程度。 P;=(ni/N)X 100% 。 (9) 标准化降水指数的计算步骤如下: 式中:N为某站有气象资料的总年数;n;为i气象 第7期 李楠,等:基于标准化降水指数的鲁西地区干旱时空特征分析 133 站点发生干旱的总年数;i代表不同的气象站点。 (2)干旱站次比(P;):某地区干旱发生站数 占全部站数的比列。 P =(m /M)X 100% 。 (10) 式中:m为发生干旱的站数,s 为i站发生干旱时 SPI的绝对值。当S <0.5时,干旱强度不明显; 0.5≤sjj<1.0为轻度干旱;1.0≤S <1.5时为中 度干旱;S ≥1.5时为重度干旱。 式中:M为总气象站数;mi为J年发生干旱的站 数;j代表不同的年份。利用Pi可以描述干旱发 生范围的大小:Pi≥70%为全域性干旱,50%≤ 3结果与分析 3.1 不同时间尺度的SPI值比较 利用聊城市8个气象观测站1962--2015年 P;<70%为区域性干旱,30%≤P <50%为部分 地区干旱,10%≤Pi<30%为局部地区干旱,Pi< 10%则表示无明显干旱。 逐月降水数据,计算得出54年来各站SPI1、SPI3、 SPI12值。结果(图2)表明,短时间尺度的SPI值 (SPI1)受较短时间降水影响大,没有考虑前期降 水的影响,随机性较强,时间持续很弱,旱涝转换 (3)干旱强度S 用于评价干旱的严重程度。 l 1 Is j=f ∑m_1SPIi l。 3 (11) 2 l 0 1 2 3 —… HH一-,.-4 o-1 e e e e o,l e 3 2 l 0 一O 3 昌器 2宝 3 2 曷 磊器昌 蓦罄g 苫兽罟2譬= 2 2曼2 2 0 2 2 2暑0宝2窖窖呙 0 舄 一l 9 /^\/\人/\。。入。八 ¨. 。。 八 :w V~V’ \\/V VV\\/\//V\ .。 。—3 O-I ∞∞H c。∞H ∞(o0 O,1 卜t"-卜∞卜o。卜0 00 00∞∞∞00 0 00 0①①∞00 0 ∞00 0 。)0 0 0 0 0 0 州 N 0 0 0 0 0H ∞∞∞∞0 0H ⑦ ∞0西0∞0 0 0 0 0 0 0 0 0 N N H H H H H N 往:份 图2 1962--2015年鲁西地区不同时间尺度SPI变化 l34 山东农业科学 一第49卷 很突然,SPI1值频繁地在0线上下波动,各月间 次,以轻旱和中旱为主。基于SPI12序列计算 变化比较大。¥PI3考虑了短时间内水分盈亏情 况,反映了季节干旱状况。随着时间尺度的进一 年尺度各等级干旱频率,各站点干旱总频率平均 为31%左右,轻旱的发生频率平均为10.38%,中 旱发生频率为14.56%,重旱发生频率为6.O1%。 3.2.2 干旱发生范围 基于SPI12序列计算年 尺度干旱站次比,结果(图3)表明,54年来鲁西 地区干旱站次比在0~100%之间波动变化,共24 年无明显干旱,占总年数的44.4%。 步延长,一次降水对SPI12的影响变小,只有在持 续多次降水的情况下才能使SPI12值发生波动。 3.2鲁西地区干旱发生的年度特征分析 3.2.1干旱发生频率基于l2个月时间尺度的 SPI序列表明,1962--2015年间干旱较为频繁,8 个站点平均出现17.25次干旱事件,平均每三年 ………干早强度——干早站次比 1OO 2・O 8O L 5 60 餐1.o H_ O・5 4O榴 辞 20 H- O.0 0 g 器器昌 曷 器器昌 墨墨吕 g兽蓦三 = 宝霉宝2 2 2 0 2 2 2 2 竺 年份 2 0 2 昌昌昌 图3 1962--2015年鲁西地区干旱站次比和干旱强度变化 由图3和表1可知,54年中共有12年发生 全域性干旱,主要集中在2O世纪80年代和90年 代,其中1986、1988、1995、1997、1999、2002、2006 3.2.3干旱发生强度54年来鲁西地区的干旱 强度在0~1.73之间波动,平均1.13,共有10年 干旱强度在1.0以上,其中1986年干旱强度值最 年干旱站次比达到100%;1986--1989连续四年 发生全域性或区域性干旱。20世纪60-70年代 干旱发生的区域较少,1962--1970年间未发生全 域性干旱,但有两次区域性干旱。1971--1980年 间,只有1972年发生一次全域性干旱。2000年 高,达1.73,为重度干旱;其中,鲁西地区发生的 干旱等级主要是轻旱和中旱(图3)。由表2可 知,轻旱和中旱发生次数较多,各年代均未有特旱 的发生;20世纪8O年代旱情较重,发生重旱14次, 中旱10次,轻旱11次。轻旱多集中在20世纪60 以后全域性干旱有所减少,干旱有减轻趋势。 表1 鲁西地区不同年代范围干旱 发生次数统计 年代 1962--1970 1971--1980 l981一l990 年代和80年代,中旱集中在20世纪80年代和90 表2 鲁西地区不同年代各级别干旱 发生次数统计 无明显干旱局部地区干旱部分地区干旱区域性干旱全域性干旱 l991—2000 2【x】1—2010 20l1_2015 第7期 李楠,等:基于标准化降水指数的鲁西地区干旱时空特征分析 135 年代。总体来看,20世纪8O年代和90年代是干旱 的高发期,并且干旱的量级也相对较重。 3.3 鲁西地区干旱发生的季节特征分析 基于SPI3的各气象站点干旱事件发生次数 平均为64.62次。季节干旱总频率平均为30% 左右,轻旱发生频率14.35%,中旱发生频率 8.55%,重旱发生频率4.44%,特旱发生频率 2.65%。以下从各季节干旱发生特征详细分析。 (1)春旱:鲁西地区1962--2015年共发生春 旱15次,分别出现在1962、1965--1967、1972— 1974、1976—1978、1981、1995—1996、2Oo0— 2001、2010年,呈阶段连续性分布,平均3~4 a出 现一次。其中,1981和2001年为重旱,SPI值分 另0为一1.79和一1.81,1962、1967、1973、1976、 1978、2000年为中旱,其余年份为轻旱。 图4(A)表明,春季干旱站次比在0~100% 之间变化,共21年全域未出现春旱,占总年数的 38.8%。20世纪60年代和70年代春旱较为频 繁而且发生区域较大,其中出现全域性干旱5次, 区域性干旱4次。20世纪80—90年代干旱的发 生较为缓和,只发生全域性干旱1次,区域性干旱 3次。2000年以后发生全域性干旱2次,区域性 干旱2次。春季干旱站次比的线性倾向率为 一5.67/10a,表明鲁西地区春季干旱站次比呈现 减少的趋势。 (2)夏旱:1962--2015年间夏旱亦较频繁,共 发生夏旱14次,分别出现在1965、1968、1972、 1983 1985 1989 1992 1997 1999 20011 2002 2006、2014、2015年,其中1983、2002年为重旱, 1997年为特旱,SPI值为一2.11。 图4(B)表明,夏季干旱的站次比也在0~100% 之间变化,共22年全域未出现夏旱,占总年数的 40.7%。20世纪6O年代初期到80年代初期近二十 年的时间里,夏旱发生区域较少,全域性夏旱只有1 次,部分地区夏旱1次,其余均为局部地区干旱。20 世纪80年代初到90年代末,夏旱的发生区域相对 较广,全域性干旱5次,区域性干旱2次。2000年以 后夏旱发生区域相对上个时期有所减少,其中全域 性干旱3次,部分地区干旱3次。夏旱站次比的线 性倾向率为3.77/10a,呈增加趋势。 (3)秋旱:1962--2015年问秋季干旱发生16 次,分布在1963、1966、1970、1978、1979、1981、 1986 1988 1989 1996 1998 2001 1 2002、2006 .2012、2013年。其中,重旱年份为1966、1981年, 特旱年份1988和1998年,其SPI值分别 为一2.31和一2.07。 图4(c)表明,秋季发生的干旱站次比也在 0~100%之间变化。共27年未有全域性秋旱的 发生,占总年数的50%。20世纪6O和70年代发 生全域性干旱4次,区域性干旱2次,局部地区秋 旱1次。20世纪80年代和90年代共发生6次全 域性干旱,1次部分地区干旱,4次局部地区干旱。 21世纪初的16年间,出现全域性干旱3次,区域 性干旱2次,局部地区干旱4次。秋旱站次比的 线性倾向率为1.01/10a,表明鲁西地区发生秋旱 的站次比呈现增加趋势。 (4)冬旱:1962--2015年问出现冬季干旱17 次,分别分布在1963、1966、1968、1970、1974、 1977 1983 1984 1986 1988 1991 1994 1996 1999、2006、2007、2012年。其中重旱年份有 1984、1996年;特早年份为1977年和1988年,其 SPI值分别为一2.05和一2.0O。 图4(D)表明,冬季干旱的站次比在0~ 100%之间变化。共33年全域未有明显冬旱发 生,占总年数的61.1%。全域性干旱17次,占总 年数的31.48 ̄0。20世纪60年代和7O年代全域 性干旱6次,局部干旱2次。20世纪80年代和 90年代干旱发生站次较高且量级相对较重,共出 现全域性干旱8次。21世纪以后发生全域性干 旱3次,局部干旱1次。冬旱站次比的线性倾向 率为一3.51/10a,表明聊城地区冬旱的站次比呈 现下降趋势。 综上可见,聊城地区发生的干旱以轻旱和中旱 较为频繁,且伴随有季节性连旱;春旱和夏旱发生 的频率高于秋旱及冬旱,但秋旱和冬旱的干旱程度 要重于夏旱及春旱,其中秋季干旱程度最为严重。 l36 山东农业科学 A:春季 2.0 第49卷 干旱强度——干早站次比 1.5 越骶噼 隈噼 1.0 西_【 0 _【 ∞∞一 0 H ∞∞ 嘧 H ∞ ∞_【 ∞ ∞一 0.5 00_【 0∞0H ∞∞ H O.0 0卜0'[ ∞0∞一 0卜①_【 卜∞_【 卜∞_【 0卜西一 卜 H 卜 .E 卜0一 0卜∞ ∞卜∞H 2・5 2.O 0∞ ∞①H ∞0H 0∞0H 61 ∞卜∞'E ∞卜∞_【 lOO 0∞0 B:夏季 0卜0 ∞卜∞-[ 0∞ H ∞∞一 呻∞∞H 心∞①一 ∞∞0— 0西0_【 ∞∞一 ∞ H 80 60 ∞ H 萋 1.5 1.o O.5 O.O ∞0'【 ∞∞∞ ∞∞ _[ 0 0.【 昏0— ∞ H 酱 40辟 2O 0 000 0 呻 ∞一 0∞∞H ∞∞0.【 000 ∞ _E ∞∞∞.【 ∞西∞.【 000 00 00 00 ∞∞ _【 000N 00 00 ∞00 N00N 00 000 ∞00N 0H0 _【0 【'0 80,、 1O0 C:秋季 2.5 2.O ∞00 ∞00 0_【0 '【0 一0 0—0 N 0 【_0 m ∞ ∞ ∞ ∞ O 1.5 6O丑 40 ^水v丑 橱畔 1.O O.5 20 ‘’ 0.0 0 D: 冬季 2.5 2.0 lO0 8060 ,、 1.5 1.0 0 5 40莓 20 0 H- O.O 图4 1962--2015年鲁西地区不同季节干旱站次比和干旱强度变化 第7期 李楠,等:基于标准化降水指数的鲁西地区干旱时空特征分析 l37 3.4 鲁西地区干旱发生的空间特征分析 t:SPI12值进行主成分分析,方差贡献牢 第一主成分各特征向} 的系数均为正值且卡¨ 差较小,说明各站降水埘l t 的贡献较为均衡,较 人值ff{现 临清和苇县,低值叫t心位于 谷和冠 县,表fJjJl临清、莘县的降水世对_r鲁两地区f:旱发 乍的影响程度相对较大 .第一t成分所代表的降 水 问变化特征可以视为鲁 地区干旱伞 分布 的主要形态。第二主成分各特征向量的系数值由 表3,…]=前两个特征量值的 {贡献率等于 0.866,人十0.7,所以取前两个代 性的主成分, 使川A J ( ls中的反距离加仪捕值法对前两个辛 成分特 世系数做空间分 , 【皋J 5所示, .表3 主成分分析特征向量系数及方差贡献率 【f1l尔部县 两部县市逐渐减小,最低值位于冠县 (一0.469),说明从东部县『 到西部县市干旱发 牛 增加趋势。 为了更好地r解鲁眄地 的分 情况, 利川8个气象站点54年米sPI12指数所得的各 站点卜旱频率进行反距离』J【I杖插值,得到研究区 干 频率分 图(图6) r,J‘ ,高唐、茌平、尔阿 的总r旱发生频率较高,轻 发生频率的各县分 布 总干旱发生频 的分 较为一致一 平、冠 ・ (氰站点 Ili 0盯 ・气象站』 I1igh 0 ■■l_('w:0 0 ■■l 0 l7 ;【】 l0 米 i r玳 图5聊城地区SPI12第一、二主成分各特征向量的空间分布 县、『ju谷的 }1旱发生频率均为1 8.5%。中西部县 的发,I 频率高于 部的高膺、 平、东阿 (1)鲁西地区月、季、年的SP1值在波动I}l均 现降低趋势,表明该地I F旱发展的趋势; 琏J:SPIl2的十早分析结果总体显示,鲁 地 以轻 和中旱为主,艰旱以f 发生概率相对 较低;第・. 成分各特征向城的空M分布中,冠县 币Il…谷位于低值中心,其次为筝县,说明其对降水 的贝献柑 r较小。小同级别于早发 频率的空间 分 结 也显示,冠县、刚谷、节县的中旱及重旱 发 频牢较高,与主成分分析的冠县、 谷、莘县 处于低 lIl心的结果有较好的一致惟。 年度和季节干旱较为频繁I] 持续时问较长,年度 及季 总干旱频率均在30%左 ,主要以轻旱和 【_{l 为主。1 (2)基于12个月尺度的Sf)1序列干旱rf1,发 生令域性干甲12年,主要集fl}l 20世纪80年代 币『 l90 q 代,其巾l986—1989年连续四年发牛 域性及全域性干旱、‘fi- 度I: 强度平均值为 1.73 、 (3)告季和冬季干旱站次比线性倾向率均为 负值,表现为f: 站次比下降趋势;夏季和秋季的 站次线性倾 J率为 值, t 增加趋势 、 (4)丰成分分析的特 ㈨级圳r 发生频率分 情况 问分布 各站不 较好的相父性 4讨论与结论 分析近54年来鲁两地区 度币【J季节1- 发 的频率、范围、强度干¨ 分布等特征得出以下 结沦 1 38 山东农业科学 第49卷 图6鲁西地区不同级别干旱的发生频率空间分布 f 的发q 过大气圈、水圈、岩石圈和生物圈 共同作月j『I,J结果 .本文仅从降水凶素分析f:旱 的发生利分 特征,没有考虑水文及灌溉等因素 的影响“,以后将结合其他干旱指标 进一步 深入研究份 地区:F旱发生的时空演变特征。 参 考文 献 [7 1 帏海涛,胡义趟,陈学君,等 种 L象1。 指怀的啦I{{比较 研究lJ r 地区农、 研究,2009,27(1):237—24I.247. [8I 孙钳辉.1:治亮. 梅,等.墟 协 准化降水指数的陕 靖f 高原地区1971—2O J() 1。 变化特征[J] 叶1国沙 漠,2013.33(5):1560一I567. [9] 雌学志,粟晓玲,;{』;东国.J, ’ l的陕 美r{ 地区气象f: ‘ 时 特征分析[J]下‘ 地l 农业研究,2011.29(2): 224—229. 张 雪J‘ :,孙畅.E J:l;峄,等.巾囝干旱I C湍变化时全球变暖 的I爱域II¨J越J:.I:‘ l≮铆}究.201o.27(4):592—599. 姚玉‘挺.IK rf杰.邓振铺,等.气象、衣、 指怀综述:J J [10 J 媛媛,张勃.綦于标准化降水指数的近40a陇东地 旱涝 时空特flF[J].自然资源学报.20I 2,27(12):2135—2l44. 1Ij苏时光.谱末成分分析及 多指怀评价体系中的应川 11j北京:中同农、Jk大学.2004. 卜 地I≮农、l 研究.2007,25(1):1 85一I89.21 1 112j 5盛.壤】:支持向量机的_I1f 列组台坝测模 f I)1. 林琳近30 我 卞要气象灾害影响特 分析:D]兰州: 上州人。' ’,_.21)1 3. K沙:湖南农业大学,2011 [I 3 J 绚i_】『.『,淡哲敏象,2014.40(9):1097—1 105. 大气圈碳循 0模拟研究进展[.1-.气象 ,J 法研究lI)j北京:斗,围水 怀准化降水指数的中 安翱 ,仆铂 ,李韵缝.等.近5O年华北区域 气象r旱 啊 Fff0j 】 分"il.I] 科学,2000,20(3):400—4I2. [14j蒋桂片1 驱动机制与 李橘,聊城地I《冬小麦i 旱灾害评仙【I)j.南京:南京信息 I 程火 :,2015. flJ水LU } :研究院.2013. [15 0趟晚 l匕,畅晓光,李茂松,等. (7):50—59. 张龙,越椭‘I-. I f1,等.2种怀准化 。1 指怀在H‘耵武成f:旱 监测lt 的× J比 J1_f:早 象,2013.31(2):412—418,431 力’季 ’ T。 近58a演变特 l1].农、I l:徉学报,2()lO.:6