2. 中国气象局气象影视中心,北京 100081
2. Huafeng Meteorological Media Group, CMA, Beijing 100081
随着经济发展和城市化进程的加快,突发性、局地性强的异常天气事件越来越多地出现在一些大城市中,如2001年8月5日上海强降水(曹晓岗等,2011)、2004年7月10日北京局地强降水(孙继松等,2006)等。有研究表明,城市热岛效应不仅造成城市及其周边地区气温的升高,还可能通过流场的作用,对冬夏季降水过程产生影响。Daniel(2000)研究认为,城市热岛效应将造成城市下游降水增加;毕宝贵等(2004)的研究指出城市化过程造成的地表特征变化也有可能对降水过程产生较大影响;孙继松(2005)、孙继松等(2007)针对北京地区,分析了地形与城市热岛及下垫面物理属性造成的热力差异对不同天气过程的作用。Thielen等(2000)的数值试验表明,城市地表的感热通量、特别是降水发生前4 h内的感热通量变化对γ中尺度的对流降水有重要影响。Chen等(2000)通过观测研究表明, 在周边地区气温变冷的情况下,随着城市群的发展,长江三角洲地区的气温在迅速变暖;与此同时, 该地区呈现出降水量增加的趋势。苗曼倩等(1998)通过数值模拟证实, 长江三角洲地区的海(江、湖)风环流与城市热岛效应之间存在相互增强的过程。蒙伟光等(2007)则研究了城市化对珠江三角洲强雷暴天气的影响,他们的数值试验表明,模拟的与城市影响有关的低层辐合主要位于500 m以下的近地面层。
天津作为沿海大城市,随着经济发展和环渤海经济新区的开发建设,城区范围不断扩大、城市建筑群迅速增多,城市热岛也展现了扩大的趋势;同时天津地处渤海西岸,受海陆风环流的影响明显,地形地貌及下垫面不均匀性的差异使得天津降水的特征与北京等地有明显不同。已有统计表明,近十年来,天津地区出现的局地强降水次数和降水的强度明显高于北京及河北等地。针对北京地区近年来出现的局地强降水天气过程,孙继松等(2006)通过个例分析指出,城市与郊区下垫面物理属性造成的热力差异利于城市中尺度低空风场辐合线的形成、同时还可能造成边界层内中心城区风场垂直切变加强,触发暴雨的出现和持续。针对天津受海陆风影响的相关研究也很多,于恩洪等(1987)、陈彬等(1989)研究表明,渤海湾海风锋对强对流天气有触发作用;何群英等(2011)、东高红等(2011)也注意到海风锋对天津局地强降水有触发和加强作用;而刘一玮等(2010)通过个例分析指出天津地区下垫面热力不均匀会对强冰雹天气造成影响。
事实上,城市热岛、海风锋之间必然存在相互作用过程,这种相互作用可能对局地强天气的发生、发展产生影响。天津新一代天气雷达位于渤海西岸,它既能探测陆地强对流天气,同时也可获得海洋气象信息,更重要的是还可捕获海-陆交界的大气边界层特征。2010年8月16日下午天津城区出现局地强降水天气,其时间短、强度大、范围小,给我们的预报带来很大难度。本文充分利用各种观测资料、雷达观测资料与VDRAS反演资料及中尺度TJ-WRF模式输出资料,对其进行综合分析,重点分析城市热岛与海风锋二者叠加作用对局地强降水的触发机制,以期为发生在沿海城市局地强对流天气的短时临近预报提供参考依据。
1 资料由于本文重点研究的是城市热岛与海风锋叠加作用对天津城区一次局地强降水的触发机制,因此选取资料时间为2010年8月16日00—12 UTC。
所用资料:常规观测资料、NECP/NCAR再分析资料,天津逐时加密自动站观测资料,包括风向、风速、温度、相对湿度。
天津多普勒雷达低层仰角反射率资料,时间间隔为6 min;利用变分多普勒雷达分析系统(VDRAS)反演得到的对流层低层热力、动力场分析资料;考虑到北京和天津2部雷达的有效探测范围,VDRAS的运算范围设置为300 km×300 km,模式水平分辨率为3 km×3 km。由于利用VDRAS反演主要是获得低层(3 km以下)的热动力三维结构分析场,因此将模式顶高设置为5.4375 km,而且仅在3 km以下的雷达资料被同化,3 km以上作为模式顶的海绵边界层处理;考虑到时空一致性,模式的垂直分辨率为0.375 km,模式最低层为垂直分辨率的一半,即0.1875 km。模式输出资料的时间间隔为12 min(陈明轩等,2011)。
天津255 m气象塔观测资料,气象塔位于天津市城区南部,海拔高度为3.3 m,地理位置为39.06°N、117.1°E。实验站周边为典型的城市下垫面条件,气象铁塔上安装有15层风、温度与相对湿度梯度观测,涡动相关观测仪采用美国CAMPBELL公司生产的CSAT3型超声风温仪,安装在气象铁塔40 m高度的南伸臂顶端,采样频率为10 Hz。
天津中尺度TJ-WRF模式输出资料(模式简介见后文)。
2 降水实况与天气形势 2.1 降水实况特征2010年8月16日下午,天津城区附近出现局地强雷雨、冰雹天气,降水从07时(文中时间为世界时,未特别说明,下同)开始,雨区自东向西移动,到10时结束,自动站最大降水量为31.6 mm。图 1a为降水量空间分布图,图 1b为降水区自动站点降水量随时间的演变。从图上看到,此次降水的特点为:(1) 降水落区范围小、局地性强,20 mm以上的降水区域仅有几平方千米;(2) 降水时间短、强度大,单站主要降水基本集中在1~2 h内。(3) 降水首先出现在天津城区东面郊区,雨区向西移动到城区时,雨强突然加强,造成城区单站1 h 29.8 mm的短时暴雨。
从8月16日00时500 hPa图上(图略)看到,在河北经河南至湖南一线有一高空槽,天津处于槽区前部,700 hPa已处于槽后,850 hPa高度以下津京附近为一弱的反气旋性环流;对比高低层槽区位置看该槽为明显后倾槽。应用NECP/NCAR再分析资料看到,到06时天津正好处于500 hPa高空槽槽区,700 hPa为槽后西北气流、850 hPa上弱的反气旋环流明显东移,天津处于反气旋环流后部的偏南气流里(图 2)。对应地面场天津地区由之前处于鞍形场的弱气压场转为处于高压后部弱的偏南气流里。
对比分析北京探空资料看,00时925 hPa以上风随高度顺转为暖平流,925 hPa以下风随高度逆转为弱冷平流,且整层大气的湿度都不大;到06时在700~850 hPa高度有一明显垂直风切变,850 hPa以下为偏南风,700 hPa以上为一致的西北风。从两时次各物理量参数变化看到:对流有效位能由00时的164.9 J·kg-1迅速增大到990.7 J·kg-1,抬升指数由-1.39℃降为-3.32℃,对流抑制由545.2 J·kg-1减小到0.4 J·kg-1,对流温度保持为34℃。这说明降水前虽然天津处于后倾槽500 hPa槽前、低层槽后弱的反环流下,但该地区对流不稳定能量迅速增加,对流抑制迅速降低,该地如有好的触发抬升条件,只要克服较小的对流抑制,就能触发该地不稳定能量强烈释放,激发强对流的新生发展。
3 局地强降水发生前后城区下垫面特征分析 3.1 天津城区附近下垫面气象要素特征从天津逐时加密自动站观测资料看到,16日00时在天津市区东面东丽—津南区有一南北方向、西南—西北风的中尺度切变存在,市区为一个小的反气旋风场环流,对应1000~850 hPa各层天津上空均为弱的反气旋风场环流。随时间地面中尺度切变向西北方移动、范围扩大、方向逐渐转为西南—东北向,到05时中尺度切变压在宁河—东丽西北部至市区北部一线,此时市区的反气旋风场环流消失,转为一致西南风(图 3a)。随后,中尺度切变向西北方收缩、特征减弱,而市区南部逐渐出现向城区的风速辐合,并于08时前后在城区风场辐合处形成一条新的地面中尺度辐合线(图 3c)。
分析中尺度辐合线形成的原因,应和天津城区与郊区热力差异强迫有关。从上面的分析知,05时前市区为一反气旋环流风场说明存在弱的冷下沉;之后由于城市下垫面物理属性的差异,太阳辐射造成该地区温度上升速度加快,气温明显高于郊区,城区附近出现城市热岛。从天津温度距平与自动站风场分布看到,城区附近的温度正距平为最大(图 3b),而且风场在此处出现辐合,辐合线走向与最大温度梯度方向基本一致。分析天津255 m气象塔观测资料也看到,05时前31℃温度的伸展高度在70 m以下,之后其高度迅速向上伸展,到07时31℃温度的伸展高度达到近160 m(图 3d)。从对应40 m高度垂直速度分布也看到,07时之前垂直速度一直为负,而07时后垂直速度转为正值,并一直持续到降水前的08:30(图 3e)。这说明随城区附近热岛效应的出现与加强,该地相对干暖空气得到抬升,因不断有气块被迫抬升,所以导致该地出现弱上升运动,上升运动又导致该地出现辐合,从而形成地面中尺度辐合线。
3.2 渤海沿岸自动站气象要素特征从渤海沿岸自动站气象要素随时间演变看,上午沿岸地区一直维持西南风,03时后逐渐转为东南风,而且风速有所加大,并于05时在岸上离海岸线20 km处形成另一条辐合线(图 3a),辐合线前侧为西南风、后侧为东南风且温度相对较低、湿度较大,说明辐合线后侧为从东南海上吹来的相对湿冷空气。对应雷达观测,在相应位置处有一强度为15 dBz的窄带回波带随时间缓慢向西移动,所以这辐合线应是海风锋。分析沿岸站点风能塔资料也看到,在04时后沿岸(风能塔观测高度最高为100 m)风向突转,由西南风转为东南风、风速最大加大到4 m·s-1(图略),海上相对湿冷空气随东南气流不断向陆地输送堆积从而形成辐合。
3.3 局地强降水前后自动站与雷达回波特征分析海风锋形成后,形状为南北向弧状。其随时间不断向西移动,移动中其后侧东南风风速最大达到6 m·s-1,雷达观测其回波强度也略有加强。从06时开始在海风锋后侧不断有雷暴生消发展,但雷暴的生命史都很短、持续时间大约仅半个小时(图 3b)。分析原因应为海风锋本身有一辐合抬升区,在其内一直有弱的对流发展,当海风锋移动到一大气层结不稳定和水汽条件适合区域时,海风锋的辐合抬升将会得到加强,从而触发雷暴新生和发展, 这与作者之前的分析结论一致(东高红等,2011)。到07:24海风锋前端移到天津市区东面,这时开始在其后侧有雷暴新生并迅速发展,两个体扫就发展成强度为50 dBz的雷暴单体。之后海风锋加速向西移动,雷暴单体也迅速发展、范围不断扩大,并与海风锋一直保持一定距离向西移动,08:36以后回波的强度有所减弱,影响范围缩小,但当海风锋移到天津市区上空时,处于其后侧回波的强度再次迅速加强、范围扩大,造成天津城区部分地区出现29.8 mm·h-1的短时强降水。从铁塔资料也看到,降水开始前后,城区近地层出现明显的上升运动,最大垂直上升速度达到0.2 m·s-1(图 3e),而且降水开始后120 m高度以下风向由西南风转为北风、风速略增加,同时相对湿度增大,气温由34℃迅速降至30℃以下(图 3d),这从10时自动站温度分布也能明显看出。降水后市区附近出现明显负距平区、降水区温度与周围自动站温度相差最大达10℃(图略),这应是由于强降水拖曳作用造成的冷池。
将雷达1.5°仰角反射率因子图与自动站风场相叠加看到,新生雷暴单体正好位于地面辐合线与海风锋之间的区域。雷暴单体随海风锋缓慢向市区方向移动,当海风锋移到市区附近,其与位于城区的中尺度辐合线相遇,相遇处回波单体迅速加强发展,12 min时间雷暴单体回波中心强度最强就达到50 dBz(图 3c)。分析原因应是随海风锋移近,海风锋后相对湿冷的气团与位于城区的中尺度辐合线上相对暖干的气团这两种不同属性气团在天津城区附近相遇,城区附近相对暖干气团不断被抬升,致使此处辐合上升运动进一步加强,从而促发该地不稳定能量强烈释放,雷暴单体新生并迅速发展加强,而热岛效应造成城区附近热力不稳定加强,使该地辐合上升运动加强,这对新生雷暴单体的加强发展起到增强作用。
4 利用雷达反演和数值模拟分析城市热岛与海风锋叠加作用对局地强降水的触发机制 4.1 利用雷达反演资料分析局地雷暴发生地的低层热、动力特征研究证实,VDRAS能够从单(双)多普勒雷达资料反演低层大气三维风场、热力场和微物理场(许小永等,2004)。本文利用北京、天津两部多普勒雷达资料及自动站资料进行反演,得到VDRAS反演的低层热动力场资料,来分析此次发生在天津城区的局地强降水天气。
从VDRAS反演结果看到,中午前后天津城区附近出现中心值为1.2℃的扰动温度大值区,说明城区附近存在城市热岛。随时间城市热岛效应明显加强,到07:53达到最强,过市区沿39.15°N做扰动温度与流场垂直剖面看到,此时城区附近扰动温度大于1.4℃的水平范围大约为15 km、垂直伸展高度接近1.4 km,而且扰动温度大值区处对应有辐合抬升。随后城区附近扰动温度大值区范围有所减小西缩,但中心值却随时间明显增高(图 4a~4e)。配合图 4f~4h看到城区附近扰动温度大值区对应着辐合上升气流,在城区东西两侧为辐散下沉气流,城区附近的辐合上升与其两侧的辐散下沉相配合,在城区和近郊构成两个垂直环流,环流的垂直伸展高度为2.7 km、中心位置在1.5 km高度附近。另外从图上看到市区东部沿海地区在近地层为偏东气流,前沿有弱的辐合抬升,而且与扰动温度低值区相对应,对比图 3b知,近地层偏东风前沿辐合抬升区所在位置为海风锋锋区。随时间海风锋不断向市区方向推进,随其移近,城区附近扰动温度大值区范围有所缩小、但辐合抬升却有所加强(图 4)。随后海风锋前沿移到城区附近,海风锋自身的辐合抬升气流与城市热岛产生的辐合上升气流相遇,该处的辐合上升明显加强,最大上升速度达到3 m·s-1,位于2.5 km高度附近(图 4i~4j);而且强回波位置与垂直上升区相对应,大于40 dBz的强回波中心位于最大上升区中心右侧(图 4k)。这说明虽然海风锋对城市热岛效应有消弱作用,但二者相遇处的辐合却得到加强,低层辐合促使该地上升运动加强,从而触发该地不稳定能量释放导致局地强对流的新生发展,造成此次局地强降水的出现。
中尺度TJ-WRF模式是天津气象科学研究所与北京大学合作引进本地化业务运行的一套数值预报业务化系统,该系统采用ARPS-ADAS数据同化系统,背景场资料使用NECP预报模式GFS资料(分辨率为0.5°×0.5°,3 h一次)。TJ-WRF模式的积分区域中心在北京附近(40.0°N、115°E), X、Y方向格点数分别为401和369,垂直方向51层,边界层内大致有11层;水平分辨率5 km,积分步长30 s,预报时效72 h,每3 h输出一次预报结果,对于降水量资料,每1 h输出一次预报结果(WRF模式主要参数设置表略)。模式中的土地利用类型分类是依据美国地理测量土地利用系统提供的24类全球土地利用分布(戴俐卉等,2008),其是利用2008年遥测资料得到的,WRF模式中天津市区周围土地利用类型为城市建筑用地(见图 5)。
本次过程模拟起始时间为2010年8月15日12时,将TJ-WRF模式模拟的降水出现时间、降水量与落区和降水实况进行对比看到,模拟降水出现的时间和观测实况基本相同,模拟的降水量相对实况偏小而且降水落区范围偏大一些(图略);对比07时实况温度与TJ-WRF模拟的温度与风场分布看到(图 6), TJ-WRF模式模拟的温度场与观测实况大体相同, 这说明TJ-WRF模式模拟结果可以用于对此次天气过程发生发展机制的分析。
分析TJ-WRF模式模拟结果看到,16日04时开始在天津城区附近出现明显高温区,到07时城区附近温度比周围郊区高出1℃左右(图 7a),这与实况观测相一致(图 3b)。而且城区上空水汽含量较小,为一相对暖干区域,这说明天津城区附近热岛效应明显。从地面风场分布看到,城区附近为西南风,而且有弱的风速辐合。同时在天津东部有一辐合线,其前侧为暖干的西南气流、后侧为从海上来的相对湿冷的东南气流,对应前面分析知这是海风锋。过市区沿39.2°N做流场与垂直速度剖面看到,07时低层西南风的偏西分量在城区上空出现向上爬升,而且在低层形成一垂直上升区。1小时后垂直上升区范围明显扩大,最大上升速度由0.1 m·s-1增加到0.15 m·s-1,上升区的高度也从800 hPa伸展到700 hPa(图 7b)。另外海风锋上也有一垂直上升区,其中心强度略强于城区附近处的。随时间海风锋逐渐向市区方向移动,随其移近,城区附近高温区范围略有缩小,但强度没有明显变化。但从图 7a上看到,海风锋移近到城区附近时,出现明显的向后弯曲现象,其后侧的风场也出现明显绕流。
过市区沿117.3°E做流场与垂直速度剖面看到(图 7c),低层流场在城区附近出现分流绕流现象,气流被分为三支:一支气流流向北、一支气流流向南、还有一支沿垂直方向向上爬升;而在远离市区的南北两侧郊区处的剖面上则没有气流分流现象(图略),这说明城市热岛对海风锋有明显的阻挡,海风锋所带气流在经过城区附近时出现了明显的气流绕流现象。在过市区沿39.2°N的流场与垂直速度剖面图上看到(图 7d),在低层有南北两支气流在市区处相遇产生辐合,强的辐合导致强的上升运动,对比图 7b知这一强辐合上升区为海风锋自身的辐合上升区与城市热岛造成的局地辐合上升气流相遇合并加强造成的。从图 7d看到上升气流的高度接近650 hPa,明显高于之前海风锋与城市热岛相遇前各自的上升区高度;而且最大上升速度也由之前的0.25及0.15 m·s-1合并加强达到0.5 m·s-1以上。这说明由于城市热岛的存在,当海风锋移动到城区附近时,其后侧气流在城区附近出现弯曲绕流和爬升现象;同时海风锋与城市热岛自身的辐合上升区在二者相遇处得到合并加强,而加强的上升运动又为雷暴的发生发展提供了有利的动力热力条件,从而触发了该地不稳定能量的强烈释放,导致强对流天气的出现。
5 结论(1) 2010年8月16日下午出现在天津城区的局地强降水天气时间短、强度大,具有明显的γ中尺度对流降水特征。降水发生在后倾槽低层槽后弱的反环流下,降水前该地不稳定能量迅速增加,大气处于明显不稳定状态,如有触发抬升条件,只要克服较小的对流抑制,就能触发该地不稳定能量强烈释放、激发强对流的新生发展。
(2) 天津市区附近城市热岛效应明显,扰动温度大值区处的低层有辐合抬升;海风锋自身也有一辐合抬升区,海风锋由岸边向西移动过程中,当其移动到城区附近时,与城市热岛效应产生的辐合抬升气流相遇,相遇处的辐合上升运动明显加强;大于40 dBz的强回波位置正好位于强上升中心的右侧。
(3) 城市热岛对海风锋的移动有明显阻挡作用,海风锋移到热岛效应明显的城区附近,其后侧低层气流会出现明显分支绕流和气流爬升现象。两者相遇处的辐合上升运动会明显加强,从而为该地不稳定能量的释放及雷暴的发生发展提供了有利的动力热力条件。中尺度TJ-WRF模式能很好地模拟出这一现象。
致谢:非常感谢天津市气象科学研究所正研级高工解以扬的悉心指导。
毕宝贵, 李泽椿, 李晓莉, 等, 2004. 北京地区降水的特殊性及其预报方法[J]. 南京气象学院学报, 27(1): 79-89. |
曹晓岗, 王慧, 邹兰军, 等, 2011. 上海"010805"特大暴雨与"080825"大暴雨对比分析[J]. 高原气象, 30(3): 739-748. |
陈彬, 于恩洪, 1989. 渤海湾西部海陆风的天气气候特征[J]. 海洋通报, 8(1): 23-27. |
陈明轩, 王迎春, 高峰, 等, 2011. 基于雷达资料4DVar的低层热动力反演系统及其在北京奥运期间的初步应用分析[J]. 气象学报, 69(1): 64-78. DOI:10.11676/qxxb2011.006 |
戴俐卉, 洪景山, 庄秉洁, 等, 2008. WRF模式台湾地区土地利用类型之更新与个案研究[J]. 台湾大气科学, 36: 43-61. |
东高红, 何群英, 刘一玮, 等, 2011. 海风锋在渤海西岸局地暴雨过程中的作用[J]. 气象, 37(9): 1100-1107. DOI:10.7519/j.issn.1000-0526.2011.09.006 |
何群英, 解以扬, 东高红, 等, 2011. 海陆风环流在天津2009年9月26日局地暴雨过程中的作用[J]. 气象, 37(3): 291-297. DOI:10.7519/j.issn.1000-0526.2011.03.006 |
刘一玮, 寿绍文, 解以扬, 等, 2010. 热力不均匀对一次冰雹天气影响的诊断分析[J]. 高原气象, 30(1): 226-234. |
蒙伟光, 闫敬华, 扈海波, 2007. 城市化对珠江三角洲强雷暴天气的可能影响[J]. 大气科学, 31(2): 364-376. |
苗曼倩, 唐有华, 1998. 长江三角洲夏季海陆风与热岛环流的相互作用及城市化影响[J]. 高原气象, 17(3): 280-289. |
孙继松, 王华, 王令, 等, 2006. 城市边界层过程在北京2004年7月10日局地暴雨过程中的作用[J]. 大气科学, 3(2): 221-234. |
孙继松, 2005. 气流的垂直分布对地形雨落区的影响[J]. 高原气象, 24(1): 62-69. |
孙继松, 舒文军, 2007. 北京城市热岛效应对冬夏季降水的影响研究[J]. 大气科学, 31(2): 311-320. |
许小永, 郑国光, 刘黎平, 2004. 多普勒雷达资料4DVAR同化反演的模拟研究[J]. 气象学报, 62(4): 410-422. DOI:10.11676/qxxb2004.042 |
于恩洪, 陈彬, 1987. 渤海湾西部海陆风的空间结构[J]. 气象学报, 45(3): 379-381. DOI:10.11676/qxxb1987.050 |
Chen Longxun, Zhu Wenqin, Zhou Xiuji, 2000. Characteristics of environmental and climate change in Changjiang Delta and it spossible mechanism[J]. Acta Meteor Sin, 14(2): 129-140. |
Daniel Rosenfeld, 2000. Suppression of rain and snow by urban and industrial air pollution[J]. Science, 287(10): 1793-1796. |
Thielen J, Wobrock W, Gadian A, et al, 2000. The possible inf1uence of urban surfaces on rainfall development:A sensitivity study in 2D in the meso-γ-scale[J]. Atmos Res, 54(1): 15-39. DOI:10.1016/S0169-8095(00)00041-7 |