引言

闪电是伴随雷暴天气系统的一种天气现象，研究闪电特征及其时空分布，对了解雷暴发生规律和雷电灾害的防御具有非常重要的意义^[1-2]。闪电活动除与天气系统有关外，还与地形等因素有紧密的关系，存在明显的地域性特征^[3-5]，不同地区的闪电活动受当地的气候、地形和环境条件的影响而呈现出较大差别。KILINC and BERINGER^[6]发现澳大利亚南部闪电密度和海拔成正相关，北部这一特征却不明显。BOURSCHEIDT et al.^[7]研究发现巴西南部地区地形坡度对闪电活动的影响大于海拔的影响。ZAJAC and RUTLEDGE^[8]对美国闪电密度与地形进行分析，发现地形与闪电密度和雷暴日之间存在较好的相关性。国内成鹏伟等^[9]研究发现北京市与成都市地闪密度与海拔高度呈明显的负相关。刘海兵等^[10-11]对不同尺度地闪密度与海拔进行分析发现，江西省地闪密度与海拔呈负相关，但江西省的三清山地闪密度在800 m以下随海拔的上升而减小，800 m以上随海拔上升表现出明显的增长趋势，不同地区及不同尺度地闪密度和海拔之间存在不同的相关性^[12]。目前山东地区尤其是泰山景区的地闪分布特征及与下垫面关系研究相对较少。

泰山位于中国东南部的山东省泰安市内，是山东省最高大的山脉，地势较为险峻，有山峰150多座。泰山景区地理范围为117.0°~117.2°E，36.20°~36.36°N，约300 km²，海拔高度在150~1 545 m之间。冯桂力等^[13]、刘向科等^[14]和唐巧玲等^[15]在研究山东地区闪电特征时，均发现泰山景区为山东的闪电高发区。本文利用山东省2007—2018年12 a的地闪数据与地理高程数据，分析地形、高度对地闪分布特征的影响，在此基础上，分析泰山的地闪密度、雷电流强度空间分布特征，同时对闪电的时空分布与海拔高度的关系进行分析。由此可对制定相应的雷电防御措施，减少雷击灾害造成的人员伤亡、森林火灾和财产损失提供科学指导，并为地方政府部门决策提供技术支撑。

1 资料来源

文中所用闪电资料由山东省雷电监测定位网提供，该系统由13个ADTD闪电定位仪组成(图 1)，可覆盖山东95%以上陆地区域，探测效率为90%左右，其定位精度为300 m，地闪雷电波形鉴别率达95%，测向精度理论随机误差为0.5°，时钟同步精度为0.1 μs。此系统能够提供闪电事件发生的时间、位置、强度等数据^[16-17]，根据相关专家对数据质量控制的研究成果^[18]，剔除强度值小于10 kA的正地闪。

文中所使用泰山景区地貌、海拔数据为美国国家航空航天局(National Aeronautics and Space Administration, NASA)与日本经济产业省(Ministry of Economy, Trade and Industry, METI)共同推出的最新地球电子地形数据，分辨率为30 m，资料获取地址为http://bjdl.gscloud.cn/search。

2 落雷次数与海拔高度相关性分析 2.1 年落雷次数与海拔高度关系变化特征

不同海拔高度所占面积与对应海拔内落雷次数占总闪比例分布趋势基本一致(图 2a)，即所占面积越大落雷次数相对较多。整体来看：1 000 m以下的区域占总面积的94.96%，对应的落雷次数占总闪的92.20%；同时发现海拔600 m以下区域占总面积的66.67%，对应的落雷次数占总闪的55.27%，此区域虽然占总面积较多，但落雷次数占总闪比例不高；而海拔600 m以上区域占总面积的33.33%，对应的落雷次数占总闪的44.93%，其中海拔1 300 m以上区域占总面积的0.28%，但对应的落雷次数占总闪的2.29%，这也反映了海拔较高的区域比海拔较低的区域更容易发生雷击。由年际变化趋势(图 2b)可知，2007、2010和2012年1 000 m以下落雷次数较多，其中2010年落雷次数最多。整体来看，虽然不同年份的闪电数量有所差别，但每年不同海拔高度落雷次数的分布趋势比较一致，即：100~500 m落雷次数占的比重最大；其次为500~1 000 m落雷次数较多；1 000 m以上落雷次数较少。

2.2 月落雷次数与海拔高度关系变化特征

泰山景区地闪活动的季节性分布特征明显(图 3a)，春季(3—5月)地闪活动逐渐增多，夏季(6—8月)地闪次数约占全年总地闪次数的86.86%，秋季(9—11月)地闪活动迅速减少，冬季仅2月出现过少量地闪，11月、12月和1月没有出现地闪活动。正闪比的季节变化与地闪频数的季节变化基本相反，正闪比最高出现在4月。泰山景区平均海拔在600 m左右，由图 3b可以看出落雷区域主要集中在海拔1 000 m以下，约占总闪的90%以上，其中8月海拔在100~500 m之间时落雷次数相对较多，由于夏季7、8月山东地区受西南季风和副热带高压的影响，强对流天气较多，泰山受地形热力和动力特征影响^[19-20]，海拔900 m以下落雷次数较多。

2.3 日落雷次数与海拔高度关系变化特征

由泰山景区落雷次数日变化(图 4a)可见：总地闪数与负闪数呈多峰震荡分布，在06时和14—20时出现峰值，其中19—20时最高，00—05时及08—13时为低值时段；正闪比在00—08时较高。图 4b为地闪频数随海拔高度每升高100 m的日变化分布，可以看出落雷区域主要集中在900 m以下；06时落雷高度在海拔400~800 m明显增加，14—20时在海拔200~800 m落雷次数最多。由于太阳辐射在中午达到最强，大气在午后容易形成不稳定层结，当受山体抬升作用影响时，比较容易出现强对流天气，所以该时段闪电活动最为活跃^[21-23]。

3 地形对地闪空间分布的影响

图 5a为泰山景区遥感影像与地闪密度叠加图。泰山景区地势呈现北高南低、西高东低的特征，北部为山区，西南部为低山丘陵，东南部为平原。泰山景区年均地闪密度为1.82次·km^-2·a^-1，高于周边区域年平均值，景区地闪密度高值区位于泰山主峰(A峰，1 545 m)峰顶区域，主峰西北方4 km附近C峰(964 m)峰顶的北侧和主峰西北方7 km附近D峰(886 m)峰顶区域的地闪密度也相对较高。

通过地闪密度与海拔高度等值线(图 5b)可看出，主峰(A峰)在海拔1 200 m以上区域地闪密度最高，年地闪密度在4.90次·km^-2·a^-1以上。主峰南侧和东侧3 km以上区域海拔迅速降低到500 m左右，其中主峰南10 km附近垂直落差在1 300 m以上。随着地形的平缓，地闪密度也随之减少，仅主峰东北侧山脊处地闪密度相对较高。C峰(964 m)的北侧地闪密度大于3.00次·km^-2·a^-1，D峰(886 m)峰顶地闪密度相对较高，泰山景区闪电密度高值区主要位于主峰和主峰西北方山峰的峰顶附近。根据陈善炳^[24]对1959—1996年6—8月泰山顶逐日雷暴方向统计结论可知，雷暴移动方向基本规律为自西、北半面向东、南半面移动。结合泰山景区地形特点，雷暴自西北方进入泰山景区，受地形抬升作用，首先在C峰、D峰处放电，虽然B峰(1 130 m)海拔较高，但相对C峰、D峰来说，高度差异并不明显且距离较近，所以B峰闪电密度较低，雷暴继续向东南移动，由于主峰相对陡峭且海拔较高，加剧了对流的发展，因此主峰闪电密度最高。

4 雷电流幅值随海拔高度变化规律

闪电强度是衡量闪电破坏能力的重要因素之一^[25]，通过对闪电监测资料统计分析，泰山景区地闪平均电流强度为11.91 kA，正闪的平均电流强度为45.42 kA，负闪的平均电流强度为11.01 kA。

为进一步了解地闪的电流强度分布特征，将地闪根据电流强度绝对值以5 kA为间隔分为21个等级，不同等级的地闪百分率变化曲线(图 6a)显示，电流强度在5~10 kA的地闪比例最大，约为31.28%，强度主要集中在0~20 kA之间，占总数的92.43%。正闪电流强度在15~20 kA之间的占正地闪总数的比例最大，约为20.83%，正地闪主要集中在15~40 kA之间，占正地闪总数的76.19%。电流强度在5~10 kA之间的负闪占负闪总数的比例最大，约为32.11%，负地闪主要集中在0~20 kA之间，占负闪总数的93.86%。由于正地闪占总闪的比例非常低，负地闪占绝大多数，所以总闪的电流强度分布曲线与负地闪接近。

由图 6b可以看出，在同一海拔高度处，正闪强度大于负闪强度，正闪强度均值随海拔高度升高波动较大。900 m以下变化较小，正闪强度均值为40.64 kA；在500 m时有轻微下降，正闪强度均值平均为25.98 kA；在900~1 100 m之间，随着海拔高度升高闪电强度均值明显上升，海拔高度1 100 m时为最高值89.40 kA；在1 200~1 500 m之间，随着海拔高度升高闪电强度均值波动中下降，海拔高度1 500 m时为最低值15.61 kA。负闪强度随海拔高度的变化不大，比较稳定，负闪强度均值在14.06 kA附近，最大负闪强度均值为20.65 kA，在海拔高度1 400 m处。由于正地闪强度远大于负地闪，正地闪造成的危害相对较高。

5 地闪密度与海拔高度相关性分析

将海拔高度按100 m划分为高度间隔，对泰山地区2007—2018年均地闪密度(0.01°×0.01°网格)进行分析，发现地闪密度随海拔高度上升呈现上升趋势(图 7a)，海拔高度小于800 m的区域地闪密度随海拔高度增加呈缓慢上升的趋势，平均值在1.80次·km^-2·a^-1左右，海拔高度大于1 000 m的区域地闪密度随海拔高度增加呈明显的上升趋势，利用线性回归方法拟合海拔高度(x)与地闪密度(y)的关系，得出R²=0.292 4，地闪密度与海拔高度呈现正相关，即高度越高，闪电密度越高，雷击事故发生的概率也越大。

将海拔高度和地闪密度均按0.01°×0.01°网格进行划分，将中心网格与其周围相邻网格数据按从小到大排序(共3×3=9个等级)，见图 7b，该等级可以较好地反映相对陡度与相对地闪密度之间的关系，可以看出随海拔高度的增加地闪密度和陡度均呈现增高的趋势，800 m以上尤其明显，可见泰山景区地闪密度与陡度呈现较好的正相关。

6 讨论与结论 6.1 讨论

虽然泰山景区地闪密度与海拔高度呈现较好的正相关，但紧邻主峰的B峰(1 130 m)海拔相对较高，地闪密度却低于2.5次·km^-2·a^-1。通过图 5b可以看出主峰A峰西北侧地闪密度高值区在C峰的东侧、北侧山坡和西侧D峰(886 m)附近，整体来看主峰A峰西北侧地闪密度高值区位于B、C山峰的西北方。形成这一现象的原因可能是：由于夏季泰山景区雷暴来向以西北方居多^[26]，受地形抬升作用，雷暴能量会在B、C山峰的西北方山坡处进行释放，当到达B、C山峰时能量已经释放完成；同时B峰、C峰山顶与主峰A峰距离较近且相对A峰比较平缓；另外可能是泰山主峰地形阻挡了低层的南风暖湿空气，使其沿山区西侧绕流北上到泰山北部，使得主峰西北侧低层大气增暖增湿，增强此处的位势不稳定，从而导致B、C山峰西北侧山坡地闪密度较高。

6.2 结论

1) 泰山景区海拔高度在200~500 m之间的年落雷次数呈明显的上升趋势，海拔600 m以上年闪电数量相对较少；泰山景区年均地闪密度为1.82次·km^-2·a^-1，地闪密度较高的区域集中在主峰和主峰西北方4~7 km附近的山顶区域。虽然不同年份的闪电数量有所差别，但年闪电密度分布趋势比较一致。

2) 地闪活动的季节性分布特征明显。夏季(6—8月)地闪次数约占全年总地闪次数的86.86%，海拔900 m以下落雷次数较多；落雷次数的高峰时段是14—21时，峰值时段闪电主要落区集中在200~800 m之间。

3) 正闪强度均值随海拔升高波动较大，900 m以下变化较小，在900~1 100 m之间，随着海拔升高闪电强度均值明显上升，在1 200~1 500 m之间，随着海拔升高闪电强度均值波动中下降；负闪强度随海拔的变化不大。

4) 地闪密度随海拔上升呈现上升趋势，海拔高度小于800 m的区域地闪密度随海拔增加呈缓慢上升的趋势，平均值在1.8次·km^-2·a^-1左右，海拔高度大于1 000 m的区域地闪密度随海拔增加呈明显的上升的趋势，地闪密度与海拔高度呈现正相关；同时发现随海拔高度的增加地闪密度和陡度均呈现增高的趋势，800 m以上尤其明显，可见泰山景区地闪密度与陡度呈现较好的正相关。

泰山景区在森林防火、旅游设施建设和游客安全等方面可参考以上结论，有效减少或避免因雷电导致的相关危害。