城市化对河南省极端降水空间分布的影响

　　摘要：针对城市化对极端降水空间分布的影响问题，利用城市冠层模型结合 WRF (weather research and forecasting)中尺度天气预报模式对河南省 2021 年“7.20”特大暴雨进行数值模拟，并结合国家级地面观测站实测降水数据对该模型的模拟精度进行验证。结果表明，考虑城市冠层影响的耦合模型能更好地模拟出区域极端降水的强度和落区，与实测降水的空间分布更为接近，耦合模拟得到的平均降水量比未考虑城市冠层影响的模拟结果高 12.1 mm;人工耗水改变了区域的水热耦合平衡，促进了城市区域对流性降水的形成。

　　关键词：城市化;极端降水;城市冠层模型;WRF 模式;绝对湿度;河南省

　　刘家宏; 骆卓然; 张永祥; 周晋军; 邵薇薇 水资源保护2021-12-03

　　极端降水是诸多极端天气气候事件中，对人们生活影响最为直接和显著的灾害之一[1]。近年来由于受到全球气候变化和人类活动的影响，全球降水分布不均性明显增加，区域性极端降水呈现显著增多的态势[2-3]。城市热岛现象、下垫面条件改变等因素都会对水汽蒸发、空气对流产生明显影响，从而影响到降水强度、频率和落区[4-5]。城市化尤其使对流层气象活动增加，引起局部区域降水量增多，城市暴雨内涝问题愈发明显，不仅破坏了生态环境，更是对人类的生命、财产等构成了严重威胁[6-10]。国内外专家对城市地区降水进行了大量的模拟研究，但早期研究主要集中于城乡降水观测资料的对比分析。Chen 等[11]基于 CMORPH 数据评估了珠江三角洲地区城市化对降水的影响，发现与周围农村地区相比，城区短历时暴雨频率显著增加，发生下午时段的极端降水明显高于农村地区;Song 等[12]利用 1950—2012 年的降水观测数据，对北京地区降水的时间趋势和空间分布进行分析，发现降水量显著减少，由于城市化的影响，城区的降水量略大于郊区，没有发现降水在空间上的显著变化。现阶段，数值模拟已被广泛应用于城市化降水研究中[13-14]。Thielen 等[15]采用 γ 尺度数值模型分析了城市化对对流性降水影响的程度，发现城市表面潜热通量和粗糙度的增加会引起城市降水量的增加;侯爱中等[16]采用单层城市冠层模型与 WRF (weather research and forecasting)天气预报模式耦合来模拟北京市极端降水，发现城区的降水有增多的趋势。与有限的降水观测数据集相比，数值模拟可以提供具有更丰富空间分布信息，更有助于揭示降水的空间分布，并且从气象学入手，基于城市冠层模型的中尺度气象模拟，对城市的布局、形态进行量化研究，目前还是比较新的研究思路[17]。本研究基于城市冠层模型耦合 WRF 中尺度天气预报模式对河南省 2021 年“7.20”特大暴雨进行模拟，并对此次暴雨的特征、空间分布和成因进行系统分析，结合国家级地面观测站数据对 WRF 模式的模拟性能进行验证，通过对比有无城市冠层模型时降水空间分布的差异，探讨城市化对极端降水的影响，以期能够加深对城市化与极端降水内在关系的理解，在一定程度上提高应对极端降水事件的能力，并为保障国家水安全提供参考依据。

　　1 研究区概况与数据资料

　　1.1 研究区概况

　　河南省地处北纬 31°23'~36°22'，东经 110°21'~116°39'之间，全省总面积 16.7 万 km2，2020 年常住人口为 9937 万，是全国重要的综合交通枢纽，如图 1 所示，河南省地势西高东低，北、西、南三面由太行山、伏牛山、桐柏山、大别山沿省界呈半环形分布，气候具有四季分明、雨热同期、复杂多样和气象灾害频繁的特点[18]。2021年7月17—22日，河南省遭遇极端强降雨事件，强降雨中心主要分布在河南省西北部和中部地区，郑州、焦作、新乡、洛阳、许昌、平顶山等多个城市出现特大暴雨，造成重大人员伤亡和财产损失。郑州市是本次特大暴雨过程中降水量最大、受灾最严重的地区，7月20日16—17时的1h内，郑州市降水量更是高达201.9 mm，创历史之最。

　　1.2 数据资料

　　模式运行所需的地形数据采用MODIS_30s 下垫面数据 (https://www2.mmm.ucar.edu/wrf/users/do wnload/get_sources_wps_geog.html)，气象驱动数据采用国家环境预报中心 (National Centers for Environmental Prediction) 提供的再分析资料 (https://rda.ucar.edu/datasets/ds083.3/)，选择了其中时间分辨率为 6h，空间分辨率为 0.25°×0.25°的再分析数据。相比 2.5°×2.5°和 1°×1°的再分析数据，该数据集精度更高，模拟结果更为可靠。为了验证模型模拟的准确性，采用中国气象数据网的中国地面气象站逐小时观测资料对模拟结果进行验证 (http://data.cma.cn/site/index.html)。将位于地级市、市辖区、县级市(城镇化率较高的地区)的气象站点作为城市站点，位于县、自治县、乡镇(城镇化率较低的地区)的气象站点作为郊区站点。位于河南省内的国家级地面观测站共有 116 个，其中城区站共 35 个，郊区站共 81 个。

　　2 模型设置

　　WRF 模式是目前比较先进与灵活的中尺度天气预报模式，具有易携带高效且可并行运算的特性。应用场景包括：理想化模拟、参数化研究、数据同化研究、预报研究、区域气候研究以及耦合模式应用等领域。WRF 模式主要由前处理、主程序、后处理 3 部分构成，模式包含了最新发展的一些物理过程参数化方案，如微物理、积云对流、辐射、陆面、以及边界层等过程的参数化方案。

　　为了获取高分辨率的计算结果，模式采用三层嵌套方案来提升区域的分辨率[19]，如图 2 所示， D01、D02、D03 各层空间分辨率分别为 27 km、9 km、3 km，对应的网格数分别为 160×160、184×184、 229×220，根据河南省的经纬度范围确定模拟区域的中心点坐标是为 33º53′N, 113º.30′E，垂直分层为 34 层，顶层气压设置为 50 hPa，考虑城市冠层模型的稳定性及与 WRF 模式的耦合情况，本研究选用单层城市冠层模型[20]。为了确保模式运行的稳定，积分时间步长设置为 135 s (最外层区域空间分辨率的 5 倍)。由于 WRF 模式中不同的参数化方案组合在不同地区的适用性存在较大差异，因此，首先在模拟前对应用最广的几种参数化方案进行测试，对模型物理参数化方案进行区域适用性评估，并且根据前人的研究对河南省的物理过程参数化方案进行设置和调整[18]，选取其中最优的参数化方案组合，微物理过程采用 Thompson 方案，长波辐射采用 RRTM 方案，短波辐射采用 Dudhia 方案，陆面过程采用 Noah 方案，行星边界层方案采用 YSU 方案，城市冠层采用 UCM 方案，模式第一层积云过程采用 KF 方案，，由于模式第二层和第三层的空间分辨率均小于 10 km，因此未选用积云过程参数化方案[21-22]。设置的具体参数来源于 Wang 等[23]以及全球城市与建筑物特征数据集，其他参数主要采用 UCM 的默认设置。模型共设置了 2 组模拟对照试验：不考虑城市冠层模型 (记为 NON) 和考虑城市冠层模型 (记为 UCM)，城市冠层模型对城市的建筑物高度、密度、几何形状、建筑材料、街道宽度和走向、绿化面积等因素做了细致的刻画[24]。

　　3 结果与分析

　　3.1 观测结果

　　本研究选取的模拟时段为 2021 年 7 月 18—20 日，其中 18 日为预热期，19—20 日为用于分析的时段。为了进一步分析观测降水量的空间分布情况，采用地统计学中常用的克里金插值法来对 116 个站点进行插值计算，克里金插值算法是利用数据点间的空间相关性，可以自动识别采样点的空间分布，消除采样点分布不均匀带来的误差，可信度较高，是地统计学常用的空间插值方法。插值得到的降水空间分布如图 3 所示，可以看出郑州市为模拟时段降水量最高的地区，两天的累积降水量超过 500 mm，东部和南部降水量较少，累积降水量小于 50 mm，西部地区累积降水量小于 100mm。所有观测站点平均降水量为 152.1 mm。对照图 1 城区站点在河南省的分布，发现城区站点多分布在郑州市周围，可以看出城区的降水量明显高于郊区。本次河南省极端降水的原因，主要是受 1000km 外的台风“烟花”和副热带高压气流的影响，大量的水汽通过偏东风源源不断从海上输送到陆地，再加上河南省太行山区、伏牛山区特殊地形对偏东气流起到抬升辐合效应，从而导致了如此强劲的降水。

　　3.2 模拟结果

　　由于模拟是对 2021 年 7 月 19—20 日的 48h 累积降水量的展示，本身具有时间属性，而且本文对时间维的考虑主要是对比城市化前后 (有无城市冠层) 对极端降水的影响，因此在这里没有详细分析逐小时的降水序列特征。引起城市降水量变化除有城市化因素外，还有地形和区域气候的变化因素[25-27]。因此，本研究模型模拟均采用相同的地形数据和气象资料来剥离其他的影响因素，进而分析城市化对极端降水的影响。由图 4 可以看出，两个方案均能够较好地模拟出河南省降水强度和落区，模型能够用于河南省极端降水的模拟。由于河南省为典型的农业大省，两个方案模拟的主要差别分布在郑州市附近。由于考虑了城市冠层物理过程，UCM 方案更好地模拟出了强降水中心，范围相比 NON 方案更大一些，而且在东南部的低值区域，UCM 方案的模拟精度也会更高一些。UCM 方案模拟的平均降水量为 144 mm，NON 方案模拟的平均降水量 131.9 mm，高出了 12.1 mm，城市化对降水起到了增强作用，其原因可能是以下两个方面：①城市阻碍效应，城市冠层对空气具有摩擦和阻碍作用，其粗糙度高于郊区。当其与地表产生的热量和水汽量相遇时，容易形成上升的热气流;此外，建筑物的阻挡容易降低降水云团的移动速度，增加降水系统在城市上空的停留时间，导致城市降水量的增加[28]。②城市热岛效应，城市拥有自然侧植被的蒸散发量以及城市人工用水产生的蒸发量，为降水对流系统提供了充足的水汽，再加上城市热岛效应使得城市产生了高温区，促进了陆气的水汽交换，对于降水的形成起到促进作用。

　　图 5 为模拟时段 UCM 方案和 NON 方案距离地面 2m 处平均比湿和温度的差值。由图 5 可见， UCM 方案模拟的比湿大约比 NON 方案高 1g/kg，说明城市化增加了河南省近地面的空气湿度，值得注意的是，增加最多的是郑州市附近(城市下垫面)。郑州市人口密集，人工耗水随城市化进程大大增加，促进了陆-气的水汽交换，为降水对流系统提供了水汽来源。此外，城市化降低了城市植被覆盖度和自然蒸散发量的同时，空气含水率更高，说明人工耗水增加的蒸发量高于由于城市化导致的自然蒸散的减少量。UCM 方案模拟的温度比 NON 方案模拟的温度高约 1.5℃，表明城市冠层物理过程增加了城市地表的温度，水汽垂直输送的速率加快，导致对流性降水增加。气温的升高也将增强大气的保水能力，从而增加暴雨的频率。这些影响相互关联，共同影响城市地区的降水特征，突出了城市降水机制的复杂性和区域差异。总的来说，城市化使河南省城区的空气湿度和气温均升高，城市化效应显著。

　　进一步选取了位于河南省城市地区共 10 个国家级地面观测站 1960—1969 年(城市化缓慢时期) 和 2010—2019 年 (城市化加快时期) 水汽压和温度的逐月气象观测资料，并结合城市人工耗水过程，探讨城市化对河南省城市地区水汽量和温度的变化规律的影响，进而探讨城市化对降水的影响[29]。需要说明的是，位于河南省的气象站共有 15 个，其中有 10 个是城区站，5 个是郊区站，由于本研究主要是针对城市化对极端降水空间分布的影响，因此主要对位于城区的 10 个站点进行分析，站点编号为 53898、53986、57051、57083、57089、57091、57178、57290、57297、58005。绝对湿度是单位体积空气中水蒸汽的质量，反映了空气中的实际水汽含量，既包括由于天然降水所引起的植被蒸散发，也包括建筑物内人工用水及道路洒水等人类活动引起的水分蒸发[30]。绝对湿度的计算公式为 ρ vap = mvap V = pvap RvapT (1) 式中：