利用龙泉驿区气象局1980~2019年近40年逐日降水观测资料与2009~2019年19个自动站的逐时降水观测资料,采用多项式插值和反距离权重插值等方法,分析了龙泉驿区暴雨的时空变化特征,得到以下结果:1) 近40年来暴雨量、暴雨日数、暴雨强度存在明显的年际变化,其旬、候变化呈现单峰型。2) 近40年来夏季暴雨量、暴雨日数在多年来暴雨中占比最高,暴雨主要集中在7、8月份。3) 在空间上,暴雨量、暴雨日数整体呈现南北高中西低,暴雨强度整体呈现东高西低。 Based on the daily precipitation observation data of Longquanyi Meteorological Bureau from 1980 to 2019 and hourly precipitation observation data of 19 automatic stations from 2009 to 2019, the methods of polynomial interpolation and inverse distance weight interpolation are adopted, and the following results were obtained: 1) For 40 years, the precipitation of rainstorm, the days of rainstorm and the intensity of rainstorm have obvious inter-annual variation, and its ten-day and pentad variation presents a single peak pattern. 2) For 40 years, the precipitation of rainstorm and the days of rainstorm in summer accounted for the highest proportion of the rainstorm in many years, and the rainstorm mainly concentrated in July and August. 3) From the spatial point of view, the precipitation and days of rainstorm were higher in the north and south, lower in the middle and west, while the intensity of rainstorm was higher in the east and lower in the west.
利用龙泉驿区气象局1980~2019年近40年逐日降水观测资料与2009~2019年19个自动站的逐时降水观测资料,采用多项式插值和反距离权重插值等方法,分析了龙泉驿区暴雨的时空变化特征,得到以下结果:1) 近40年来暴雨量、暴雨日数、暴雨强度存在明显的年际变化,其旬、候变化呈现单峰型。2) 近40年来夏季暴雨量、暴雨日数在多年来暴雨中占比最高,暴雨主要集中在7、8月份。3) 在空间上,暴雨量、暴雨日数整体呈现南北高中西低,暴雨强度整体呈现东高西低。
龙泉驿区,暴雨,变化特征
Kang Xu, Wenshu Mao, Haibo Zhao
School of Atmospheric Sciences, Chengdu University of Information Technology, Chengdu Sichuan
Received: Jun. 8th, 2021; accepted: Jul. 14th, 2021; published: Jul. 21st, 2021
Based on the daily precipitation observation data of Longquanyi Meteorological Bureau from 1980 to 2019 and hourly precipitation observation data of 19 automatic stations from 2009 to 2019, the methods of polynomial interpolation and inverse distance weight interpolation are adopted, and the following results were obtained: 1) For 40 years, the precipitation of rainstorm, the days of rainstorm and the intensity of rainstorm have obvious inter-annual variation, and its ten-day and pentad variation presents a single peak pattern. 2) For 40 years, the precipitation of rainstorm and the days of rainstorm in summer accounted for the highest proportion of the rainstorm in many years, and the rainstorm mainly concentrated in July and August. 3) From the spatial point of view, the precipitation and days of rainstorm were higher in the north and south, lower in the middle and west, while the intensity of rainstorm was higher in the east and lower in the west.
Keywords:Longquanyi District, Rainstorm, Variation Characteristics
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
根据联合国政府间的气候变化委员会在第五次科学评估的发现,在以往的几十年中,极端天气如强降水、热浪等呈增多趋势,未来极端天气出现的次数将会更多 [
在四川省,暴雨洪涝是发生频率最高、造成危害最重的自然灾害之一,其特点有范围广、强度大等,影响四川暴雨的系统有青藏高原低压、西南涡、西太平洋副热带高压、低空急流等 [
龙泉驿区位于四川省成都市平原东部,由图1可知:龙泉驿处于北纬30˚27'52~30˚43'23与东经104˚08'19~104˚27'09之间,境内东西长29.8 km、南北宽28.75 km,最高海拔为1050 m左右,最低海拔450 m左右。其总面积达556.98 km2,境内龙泉山位于东南部,呈东北–西南走向,平坝、丘陵、低山等地貌在龙泉都有表现,各地貌面积占比平坝最大、低山与丘陵次之。
本文选用1980~2018年龙泉驿区气象局的逐日降水量观测资料,以及如图1所示的黄土、柏合、西河、万兴、大面、茶店、洪安、洛带、桃花故里、七里、红花、照壁、天公山、梨园村、同安、十陵、经开区、大佛村、成龙一号共19个自动站的2009~2019年逐时降水量观测资料,得到暴雨量、暴雨日数、暴雨强度的时间演变和空间演变的数据。
图1. 龙泉驿区地形和自动站分布
通过沈剑华 [
为了研究上述函数的性质,我们希望做出一个能反映原函数特性,且便于计算的简单函数 P ( x ) 来近似代替 f ( x ) ,对于函数 P ( x ) 要满足 P ( x i ) = y i ( i = 0 , 1 , 2 , ⋯ , n ) 的条件,称这类问题为插值问题, P ( x ) 为插值函数。设函数 y = f ( x ) 在区间 [ a , b ] 上有定义,并且它在 n + 1 个不同点 a < = x 0 < x 1 < x 2 < ⋯ < x n < = b 上的函数值 y 0 , y 1 , ⋯ , y n 都已知,若存在一个次数不超过 n 次的多项式(其中 a i 为实数):
P n ( x ) = a 0 + a 1 x + a 2 x 2 + ⋯ + a n x n (1)
满足条件:
P n ( x i ) = y i , ( i = 0 , 1 , 2 , ⋯ , n ) (2)
称 P n ( x ) 为函数 f ( x ) 的 n 次代数插值多项式。本文通过1980~2018年近40年的龙泉驿区气象局逐日降水量观测资料,得到暴雨的时间演变数据,利用绘图工具绘制暴雨量、暴雨日数、暴雨强度的时间演变图形,并用多项式插值来叠加曲线,分析暴雨随时间变化的特征,最后得到龙泉驿区暴雨的时间演变概况。
通过已知样品点采集的数据来推算未知点的数据称为空间插值,这种方法常用于气象领域,反距离权重法作为一种常用于气象要素的空间插值方法,其原理简单、使用方便,在气象等领域被广泛使用。反距离权重插值其插值结果好坏与搜索半径以及样品点数量等因素有关系,搜索半径越大差值结果越光滑,效果越好。若样品点数量少,则反距离权重插值结果随机性大,会出现较大差异 [
Z o = ( ∑ i = 1 n z i d i r ) / ( ∑ i = 1 n 1 d i r ) (3)
式子中 z o 为 O 点估算值, z i 控制点 i 的值, d i 为控制点 i 与点 O 的距离, n 为估算时所用控制点数量, r 为指定的幂次数 [
本文通过2009~2019年近10年的龙泉驿区内19个自动站点逐时降水量观测资料,得到暴雨的空间演变数据,利用GIS等软件的反距离权重插值,绘制出暴雨量、暴雨日数、暴雨强度的空间分布图,通过分析不同时间下暴雨的空间变化特征,最后得到龙泉驿区暴雨的空间演变概况。
根据图2(a)~(c)龙泉驿区春、夏、秋各季暴雨量的时间演变图来看:春季暴雨量最高值为1997年的91.0 mm,最低值为1983年的50.2 mm,平均春季暴雨量为5.2 mm。39年间只有1983、1986、1997年这3年存在春季暴雨,因此春季暴雨的变化趋势参考性不大;夏季暴雨量最高值为1984年的364.0 mm,最低值为1986年的51.8 mm,平均夏季暴雨量为113.4 mm。夏季暴雨量在1980~1989年左右呈上升趋势,在1990~2002年左右呈下降趋势,2003~2018年夏季暴雨量波动较平稳;秋季暴雨量最高值为2001年的116.6 mm,最低值为2012年的50.9 mm,平均秋季暴雨量为24.8 mm。39年间秋季暴雨量呈先增再减趋势;冬季无暴雨,故不讨论。从图2(d)龙泉驿区年暴雨量的时间演变曲线来看:龙泉驿区暴雨量最高值为1984年的364.0 mm,最低值为1993年的59.3 mm,平均暴雨量为143.4 mm。差值曲线上,暴雨量在1985~1990年左右呈上升趋势在1991~2002年左右呈下降趋势,2005~2015年暴雨量升降趋势较为平稳。从图2(e)~(f)暴雨量的旬、候变化来看:8月中旬的暴雨量最高,8月第4候的暴雨量最高,暴雨量的旬、候变化都呈现单峰型,且只在4~10月有暴雨量其他月份无暴雨量。
综合来看,暴雨量存在明显的年际变化,且暴雨量的旬、候变化都呈现单峰型,夏季暴雨量约占总暴雨量的79%,而秋季暴雨约占17%,春季暴雨约占4%,冬季无暴雨出现。近40年间1~12月中7、8月份暴雨量占比最高。
图2. 龙泉驿区暴雨量的时间演变。(a)~(c):春、夏、秋季暴雨量(单位:mm),(d):年暴雨量(单位:mm),(e):旬暴雨量(单位:mm),(f):候暴雨量(单位:mm)
从图3(a)~(c)龙泉驿区春、夏、秋各季暴雨日数的时间演变曲线来看:春季暴雨日数最高、最低值都为1 d,平均春季暴雨日数为0.1 d。由于春季暴雨日数较少,不讨论其变化趋势;夏季暴雨日数最高值为1984年的5 d,最低值为1 d,平均夏季暴雨日数为1.7 d。在1980~1989年左右夏季暴雨日数呈上升趋势,1990~2003年左右呈下降趋势,2004年~2018年暴雨日数波动不明显;秋季暴雨日数最高值为2 d,最低值为1 d,平均秋季暴雨日数为0.4 d。插值曲线显示,39年来秋季暴雨日数变化波动不大。冬季无暴雨,故不讨论。从图3(d)龙泉驿区年暴雨日数的时间演变曲线来看:暴雨日数最高值为1984、2018年的5 d,最低值为1 d,平均暴雨日数为2.2 d。在1985~1989年左右暴雨日数呈上升趋势,1990~2000年左右呈下降趋势,2001~2015年左右暴雨日数波动平稳。从图3(e)、图3(f)暴雨日数的旬、候变化来看:8月中旬的暴雨日数最高,8月第4候的暴雨日数最高,暴雨日数的旬、候变化都呈现单峰型,且只在4~10月有暴雨日数其他月份无暴雨日数。
综合来看,暴雨日数存在明显的年际变化,且暴雨日数的旬、候变化都呈现单峰型,夏季暴雨日数约占总暴雨降水日数的77%,秋季暴雨日数约占18%,春季暴雨日数约占5%,而冬季无暴雨出现。近40年间1~12月份中只有4~10月份存在暴雨日数,其中7、8月份暴雨日数占比最高。
图3. 龙泉驿区暴雨日数的时间演变。(a)~(c):春、夏、秋季暴雨日数(单位:d),(d):年暴雨日数(单位:d),(e):旬暴雨日数(单位:d),(f):候暴雨日数(单位:d)
从图4(a)~(c)龙泉驿区季暴雨强度的时间演变曲线来看:春季暴雨强度最高值为1997年的91.0 mm/d,最低值为1983年的50.2 mm/d,平均春季暴雨强度为68.2 mm/d。由于39年间春季暴雨数量少,因此不分析其降水强度的升降趋势;夏季暴雨强度最高值为1987年的91.4 mm/d,最低值为1986年的51.8 mm/d,平均夏季暴雨强度为65.2 mm/d。差值曲线显示,在1983、2002年左右夏季暴雨强度存在波动谷值,1991、2011年左右夏季强度存在波动峰值;秋季暴雨强度最高值为1995年的96.9 mm/d,最低值为2015年的50.0 mm/d,平均秋季暴雨降水强度65.6 mm/d。差值曲线显示,39年来秋季暴雨强度呈多次波动状态。冬季无暴雨,故不讨论。从图4(d)龙泉驿区年暴雨强度的时间演变曲线来看:近40年来,龙泉驿区年暴雨强度最高值为1987年的91.4 mm/d,最低值为1994年的53.8 mm/d,平均暴雨强度为57.8 mm/d。差值曲线显示,1980~1994年左右暴雨强度呈下降趋势,1985~1993年左右暴雨强度呈上升趋势,1994~2018年左右暴雨强度主要呈下降趋势。从图4(e)、图4(f)暴雨强度的旬、候变化来看:4月下旬暴雨降水强度最高,5月下旬暴雨降水强度最低,6月第2候的暴雨降水强度最高,5月第5候的暴雨降水强度最低。
综合来看,暴雨强度存在明显的年际变化,且暴雨强度的旬、候变化都呈现单峰型,龙泉驿区暴雨主要集中在7、8月份,近40年间平均暴雨强度为66.3 mm/d。春、夏、秋季平均暴雨强度都超过了60.0 mm/d。
图4. 龙泉驿区暴雨强度的时间演变。(a)~(c):春、夏、秋季暴雨强度(单位:mm/d),(d):年暴雨强度(单位:mm/d),(e):旬暴雨强度(单位:mm/d),(f):候暴雨强度(单位:mm/d)
从图5(a)~(c)龙泉驿区春、夏、秋季总暴雨量的空间分布来看:春季暴雨量呈现三点分布,这三点为西北部的西河和龙泉山北部的梨园村南部的红花,这三处附近暴雨量在49.61~74.39 mm之间,其余地区降水量基本在0.00~12.40 mm间,无暴雨;夏季暴雨量呈现南北高中西低,在龙泉北部黄土、洛带和龙泉山南部的红花附近暴雨量最高达到1393.38 mm,南北大部分地区暴雨量在830.57 mm以上,龙泉西部十陵暴雨量最低,在267.77~455.37 mm间,龙泉中西部暴雨量基本在830.58 mm以下;秋季暴雨量呈现东北高西南低,龙泉北部洪安、洛带附近暴雨量最高达到302.07 mm,龙泉东北半区暴雨量基本在151.09 mm以上,西南部最低暴雨量在0.10~50.43 mm间,龙泉西南半区暴雨量基本在151.10 mm以下。冬季无暴雨,故不讨论。从图5(d)龙泉驿区总暴雨量的空间分布来看:龙泉驿区总暴雨量呈现南北高中西低,近10年间龙泉驿区北部黄土、洛带、洪安和龙泉山南部红花的暴雨量最高,在1371.31~1563.94 mm之间,中西部地区暴雨量偏低,西部十陵暴雨两最低达到408.10 mm,南北大部分地区暴雨量在986.02 mm以上,中西大部分地区暴雨量在986.03 mm以下。
综合来看,暴雨量在龙泉驿区的分布特点为北部、东南部多,西部和中部少。在北部黄土、洪安、洛带以及东南部的七里、红花、茶店附近暴雨量较高,西部十陵中部成龙一号附近暴雨量较低。
图5. 龙泉驿区暴雨量的空间分布。(a)~(c):春、夏、秋季总暴雨量(单位:mm),(d):总暴雨量(单位:mm)
从图6(a)~(c)龙泉驿区春、夏、秋季总暴雨日数的空间分布来看:春季暴雨日数主要分布在龙泉西北部的西河以及龙泉山北部的梨园村和南部的七里、红花,暴雨日数最高达到1 d,其余地区基本上无暴雨存在;夏季暴雨日数呈现南北高中西低,在龙泉北部的黄土和龙泉山南部的红花附近暴雨日数最高达到22 d,南北大部分地区暴雨日数在13 d以上,中西大部分地区暴雨日数在14 d以下,龙泉西部的十陵暴雨日数最低达到4 d;秋季暴雨日数呈现东北高西南低,在龙泉北部的洛带附近暴雨日数最高达到5 d,东北半区的大部分地区暴雨日数在3 d以上,西南半区暴雨日数基本在3 d以下,西南部的柏合等地无暴暴雨;冬季无暴雨,故不讨论。从图6(d)龙泉驿区总暴雨日数的空间分布来看:龙泉驿区暴雨日数整体呈现南北高中西低,暴雨日数最高的地区为龙泉北部的黄土、洛带附近,其暴雨日数在23~25 d之间,暴雨日数最低的地区为龙泉西部的十陵附近,在6~9 d之间,南北大部分地区暴雨日数在16 d以上,中心部地区暴雨日数基本在17 d以下。
综合来看,暴雨日数在龙泉驿区的分布特点为北部、东南部多,西部和中部少。在北部的黄土、洛带附近以及龙泉山南部的红花、七里、茶店附近较高,在西部的十陵和中部成龙一号附近较低。
图6. 龙泉驿区暴雨日数的空间分布。(a)~(c):春、夏、秋季总暴雨日数(单位:d),(d):总暴雨日数(单位:d)
从图7(a)~(c)龙泉驿区春、夏、秋季暴雨强度的空间分布来看:春季暴雨强度在西北部的西河以及龙泉山北部的梨园村、南部的红花附近最高,在62.00~74.39 mm/d之间,其余地区基本无暴雨;夏季暴雨强度整体呈现东高西低,龙泉山南部的桃花故里附近暴雨降水强度最高,在71.18~73.45 mm/d之间,东部大部分地区暴雨强度在66.62 mm/d以上,西北部西河的暴雨强度最低,在59.80~62.08 mm/d之间,西部大部分地区暴雨强度在66.63 mm/d以下;秋季暴雨强度整体呈现东北高西南低,在东北大部分地区暴雨强度在63.97~76.75 mm/d之间,西南部柏合附近降水强度最低在0.05~12.84 mm/d之间,其余地区暴雨强度基本在51.19~63.96 mm/d间。冬季无暴雨,故不讨论。从图7(d)龙泉驿区暴雨强度的空间分布来看:龙泉驿区暴雨强度整体呈现东高西低,在北部万兴和龙泉山南部大佛村、桃花故里、茶店附近暴雨强度最高,在68.08~69.68 mm/d之间,龙泉山一带暴雨强度在66.48 mm/d以上,龙泉西北西河附近暴雨强度最低,在60.05~61.66 mm/d之间,西部大部分地区暴雨强度在64.87 mm/d以下。
综合来看,5~9月暴雨强度在龙泉山北部的梨园村和南部的茶店、桃花故里附近最高。4月和10月暴雨强度在北部的西河、黄土、万兴附近最高,其余地区基本无暴雨。
图7. 龙泉驿区暴雨强度的空间分布。(a)~(c):春、夏、秋季暴雨强度(单位:mm/d),(d):暴雨强度(单位:mm/d)
综合上述对龙泉驿区暴雨的时间和空间演变的分析,得到以下结论:
1) 近40年间,暴雨量、日数、强度存在明显的年际变化,其旬、候变化表现为单峰型,在四季中,夏季暴雨量、日数最高,春秋两季较低。近40年间的各月中,4~10月存在暴雨,且暴雨主要集中在7、8月份。
2) 近10年间龙泉驿区暴雨量、日数呈现南北高中西低,强度呈现东高西低。龙泉驿区北部和龙泉山南部暴雨量、日数较高。北部和龙泉山南部附近暴雨强度较高,西北部的暴雨强度较低。
本文主要从暴雨量、日数、强度三个方面来分析龙泉驿区的暴雨时空变化特征,探究了龙泉驿区暴雨的时空变化规律。在暴雨的形成、发展等方面还未进行探讨,对于这些方面还应采用更加多样的分析方法进行详细研究,这是我下一步的研究目标。
许 康,毛文书,赵海波. 成都龙泉驿区暴雨变化特征Characteristics of Rainstorm Variation in Longquanyi District of Chengdu[J]. 自然科学, 2021, 09(04): 443-453. https://doi.org/10.12677/OJNS.2021.94050