三阈值法干旱识别及特征值提取见图1,步骤如下。
步骤1 逐时段干旱状态初步识别
的时段为干旱,识别出a、b、c、d共4次干旱。
步骤2 干旱合并与剔除
对历时1个时段的干旱(如a、d),若其(如a),则最终被识别为1次干旱,否则,剔除该时段(如d)。对于间隔短于L个时段的相邻干旱,且间隔期
,合并为1次干旱(如L=3,b、c可合并为一次干旱b+c);间隔期为R>R1的时段,则为独立的干旱(如c、d)。
步骤3 干旱特征值提取
记干旱烈度为干旱起始至结束过程中R≤R0的累积水分亏缺量,峰值强度I为过程中水分亏缺最大值。、
分别取《气象干旱等级》中轻旱所对应的降水距平百分率的下、上限值,即R0=-40%、R2=-60%;R1为干旱解除阈值,据降水量超过蒸发量和土壤耕作层有效含水率的时段内雨养农业不干旱,取对应R1的月降水量为月水面蒸发量与耕作层有效含水率W之和,即:
W=(θf-θw)h
式中,θf、θw分别为田间持水量、凋萎含水率(体积含水率),%;h为耕作层厚度,mm。
3 合肥市干旱识别及特征值重现期分析
3.1 干旱识别结果
采用1953~2019年合肥站逐月降水和水面蒸发资料识别干旱。其中1953~2019年合肥站逐月降水取自中国气象数据网(http://data.cma.cn),1963~2019年月水面蒸发量取自合肥水文水资源局设在董铺的蒸发量观测资料,设站前1953~1962年的水面蒸发量以1963~2019年多年平均的月水面蒸发量代替。根据安徽省土壤含水率田间持水量数据和太湖地区典型水稻土土壤水分特征,取合肥市土壤耕作层有效含水率为22mm(太湖地区水稻土耕作层土壤有效体积含水率为15.7%,厚度为14cm[8])。合肥市1953~2019年各月平均蒸发量及R0、R2所对应的降水量见表1。
对于S<10mm干旱,在其过程中遭遇1场中雨量级的降水,所识别的干旱即不存在,故剔除S<10mm的干旱。干旱识别结果受小干旱合并的间隔参数L影响,L=1、3时识别的干旱分别为83、80场。比较识别的干旱和实际旱情,发现L=3识别的干旱更符合实情,80场干旱的S、I的分布函数的拟合见图2。