1. 引言
湿地，指天然或人工、长久或暂时之沼泽地、泥炭地或水域地带，带有或静止或流动或为淡水、半咸水或咸水的水体，包括低潮时水深不超过6 m的水域。湿地生态系统是地球表面三大生态系统之一，位于陆地与水体之间，兼有水、陆两者的生态功能。湿地不仅为众多野生动植物提供栖息地，而且有巨大的环境效益，在抵御洪水、调节径流、蓄洪防旱、调节气候、控制土壤侵蚀、促淤造陆、降解环境污染等方面起着极其重要的作用 [1]。湿地被人们称为“地球之肾”、“生命的摇篮”、“文明的发源地”，受到全世界范围的广泛关注。近年来，随着人们对湿地认识的不断加深，对湿地的开发与保护活动不断开展。
城市下垫面是最具变化的下垫面。人工形成的建筑物、树林、草地、水体、沥青马路等各种类型下垫面的空间形态构成了具有强烈非均匀性的城市特有的格局，同时各种不同下垫面的热容量、含水率、反射率、粗糙度等物理特性不同，对能量平衡、水分交换和局地环流等大气过程影响也不相同，从而形成了特殊的城市大气环境 [2] [3]。研究表明，在自然植被改变为城市地表的区域，城区湿地覆盖面积减小引起地表湿度减小 [4]；建筑物的增加导致地表粗糙度增大，反照率加大 [5]。上述这些变化都直接影响地表能量平衡，感热通量加热大气的比例增加导致近地面气温升高，城市热岛效应加强 [6]。与此同时，气流的拖曳能力被加强，从而导致城区内风速的减小并影响地表的湍流动能，最终改变大气污染的输送扩散条件 [7]。但是，如果在城市区域增加水体、绿地等自然地表覆盖面积，则有利于城市的增湿减温 [8]，促使局地流场发生变化，从而改善局地微气象条件和大气物理环境 [9]。
大气数值模式，是通过数值方法对大气动力及热力过程进行求解，从而研究大气状况的重要工具，被广泛应用于天气预报、区域气候、大气化学、城市规划等领域。Chen等 [10] 详细阐述了近年来WRF模式在城市气象模拟中的应用及可行性。佟华等 [11] 使用城市边界层数值模式研究楔形绿地规划对北京城市气候的影响。冯娴慧等 [12] 利用GRAPES模式对广州城市区域的绿地布局进行规划评估，认为有利于改善局地气候的绿地布局方式是保留并建设位于城市上风方向的大规模绿地，而非进行点小面广的星星点点式布局。朱浩楠等 [13] 利用WRF模式对重庆梁平区待建湿地公园进行规划评估，认为将湿地公园分为两部分且东北部面积较大时，能产生新的风道，对当地的局地气候舒适性有较好改善。
三垟湿地有温州市“绿肾”之称，位于瓯江入海口南岸的温州生态园内，地处温州市区东南角，紧邻大罗山，360度向城市开放，是温州大都市的重要核心区域。随着温州城市化进程的加速，三垟湿地作为城市快速增长的前沿地带，正面临着城市建设开发占用自然湿地的问题。本文利用WRF模式研究三垟湿地对温州城市大气环境的影响，以科学评价三垟湿地对改善温州城市生态环境的重要性，这对政府决策和湿地保护均具有重要的理论和现实意义。
2. 研究资料和方法
2.1. 研究区域概况
温州市位于浙江省东南部，东濒东海，介于北纬27˚03'~28˚36'，东经119˚37'~121˚18'，陆地总面积11,784 km2。市区呈东西向走向，下辖鹿城、龙湾、瓯海、洞头4区，三垟湿地位于瓯海区。温州市为亚热带海洋性季风气候区，冬夏季风交替显著，温度适中，冬无严寒，夏无酷暑，四季分明，雨量充沛，年平均气温17.3℃~19.4℃，年降水量在1113~2494 mm之间，年日照数在1442~2264 h之间。正因为独特的气候优势，温州获封“中国气候宜居城市”荣誉称号。
2.2. 资料选取
本文用于检验模拟方案的观测资料主要由4个围绕三垟湿地的自动气象站提供，站号和站址分别为58659 (28.03˚N, 120.67˚E)、K3088 (28˚N, 120.67˚E)、K3166 (28.03˚N, 120.58˚E)、K3911 (28.02˚N, 120.6˚E)。
2.3. 模式和试验设计
WRF (Weather Research and Forecasting)模式是由美国国家大气研究中心(NCAR)、美国国家环境预报中心(NCEP)和美国国家海洋和大气管理局(NOAA)联合开发的新一代中尺度天气预报模式和同化系统。具有全球各地设置使用、实时输入资料、采用先进的物理过程参数化方案、可在不同的计算机平台上运行等诸多特点。本文采用WRF-V3.2版本，使用三重嵌套网格，从外到内模式网格分辨率分别为9 km、3 km、1 km。网格采用兰伯特投影坐标系，最外层模拟中心为28˚N，120.5˚E，参考纬度为30˚N和60˚N，参考经度为120˚E。图1所示为模拟区域和嵌套网格设置。NCEP的全球1˚ × 1˚再分析资料(FNL)经处理后作为WRF模式的初始场和边界条件，微物理参数化方案采用WSM3方案，长波辐射参数化采用RRTM方案、短波辐射参数化采用Dudhia方案，陆面过程采用Noah-UCM方案，近地层采用Monin-Obukhov方案，行星边界层方案采用YSU方案，城市下垫面采用UCM方案。所有嵌套网格都未采用积云参数化方案。
Figure 1. Grid settings of the simulation area
图1. 模拟区域网格设置
模式使用NCEP提供的全球30 s地形数据及MODIS下垫面分类资料。模拟时间段为2017年9月25日20:00~26日20:00 (北京时，下同)共计24 h，每隔6 h更新一次边界条件，每隔1 h输出一次结果。25日20:00~23:00作为模式启动时间而舍去。模拟试验设计成两组，一组是不改变MODIS下垫面土地利用资料，这个模拟方案称为“控制试验”；另一组是改变MODIS下垫面土地利用资料，把三垟湿地区域的土地利用资料全部修改成城市用地，即不保护三垟湿地，让三垟湿地全部被城市开发占用，从而达到抹掉三垟湿地的效果，这个模拟方案称为“敏感试验”。图2所示为控制试验和敏感试验最内层网格的土地利用分布。
(a)
(b)
Figure 2. Land use distribution of control experiment (a) and sensitive experiment (b), (the red color represents urban area and the green color represents wetland)
图2. 控制试验(a)和敏感试验(b)的土地利用(红色代表城市，绿色代表湿地)
2.4. 研究方法
敏感试验中，修改下垫面土地利用资料的方法是利用NCL编程语言，对WRF输入文件wrfinput中的LU_INDEX进行修改，将三垟湿地位置的土地利用类型由湿地修改成城市用地，而其他位置的资料不发生改变。控制试验是真实土地利用情况的数值模拟，敏感试验是没有三垟湿地的数值模拟，将控制试验模拟结果减去敏感试验模拟结果则得到三垟湿地的存在对周围气象环境的影响情况，结合分析近地面气温、湿度和风场变化特征，来说明三垟湿地对局地气象环境的影响效应。
3. 模式结果分析
图3是4个自动气象站在模拟期间近地面气温模拟值和实测值的比较。由图可见，虽然模拟值系统性都要比观测值偏低，这与模式在边界层的物理方案误差有关，但总体模拟效果尚可，WRF模式的模拟结果与观测值在变化趋势上基本一致，相关系数分别达到0.67、0.75、0.63和0.56，说明模式较好地模拟出该地的气象环境，利用模式输出的高分辨率资料用于分析该地的局地气象环境变化情况是可行且可信的。
(a) 
(b) 
(c) 
(d)
Figure 3. Comparisons of observed and simulated values of surface temperature during 00:00~20:00 26 September 2017 at automatic weather stations of (a) 58659, (b) K3088, (c) K3166 and (d) K3911 (Unit: ˚C)
图3. 2017年9月26日00:00~20:00自动气象站58659站(a)、K3088站(b)、K3166站(c)和K3911站(d)近地面气温模拟值与观测值比较(单位：℃)
图4所示为白天(9:00~17:00)与夜间(0:00~8:00)控制试验2 m温度相对敏感试验的差异。可明显看出，无论是白天还是夜间，三垟湿地对湿地附近区域呈现明显的降温作用；白天降温面积增大，对七都岛区域也有2℃左右的降温作用，降温的幅度也更大，最大降温幅度达到了7.1℃；夜间降温作用减弱，降温面积缩小，降温的幅度也减小，最大降温幅度为5.2℃。三垟湿地造成降温的原因主要是由于湿地的比热容大，在白天能吸收更多热量。湿地在吸收太阳能后获得热量，通过水面蒸发、水面有效辐射和水面与大气的对流热交换等失去热量而导致降温。白天水面蒸发、水面与大气的对流热交换作用明显，降温作用变大；夜间水面蒸发、水面与大气的对流热交换作用减弱，降温作用变小。
(a) 
(b)
Figure 4. Surface temperature differences in the daytime (a) and nighttime (b) between cases of control simulation and sensitive simulation (Unit: ˚C)
图4. 控制试验与敏感试验的白天(a)和夜间(b)地面气温差值(单位：℃)
图5所示为白天与夜间控制试验2 m相对湿度相对于敏感试验的差异。由图可见，三垟湿度对研究区域的增湿效果不明显，最大增湿幅度仅为2.5%，几乎无增湿作用。
(a) 
(b)
Figure 5. Surface relative humidity differences in the daytime (a) and nighttime (b) between cases of control simulation and sensitive simulation (Unit: %)
图5. 控制试验与敏感试验的白天(a)和夜间(b)地面相对湿度差值(单位：%)
图6为白天与夜间控制试验10 m平均风速相对敏感试验的差异，棕色线条代表平均风流场。从图中可知，三垟湿地对湿地附近区域的风环境影响较大，正因三垟湿地的存在，湿地附近区域的风速明显增大，白天和夜间风速增加幅度分别达到了1.6 m/s和1.9 m/s；白天风速增加区域变大，分别向湿地的东北方向和西南方向扩大，夜间风速增加区域则缩小至湿地附近区域；无论是白天还是夜间，三垟湿地均增大了东南风风道，对研究区域的通风环境有较好改善效果。由于风速大小的改变主要是受水陆粗糙度差异的动力作用和水陆热力差异的热力作用的影响。因三垟湿地的存在，使得下垫面粗糙度减小，气流的拖曳能力被减弱，从而导致研究区域内风速的增加。
(a) 
(b)
Figure 6. Surface wind differences in the daytime (a) and nighttime (b) between cases of control simulation and sensitive simulation (Unit: m/s)
图6. 控制试验与敏感试验的白天(a)和夜间(b)地面风场差值(单位：m/s)
由上述分析可知，三垟湿地有明显的降温效果，无论是白天还是夜间，对城市热岛效应都有明显的减缓作用。三垟湿度几乎无增湿作用。10 m风场的分析表明，三垟湿地有利于东南风风道的加强，从而改善通风环境，通风环境改善之后将改善局地空气扩散能力。
4. 结论和讨论
本文利用WRF模式下垫面敏感性试验对温州三垟湿地对局地气象环境的影响进行分析研究，主要得到以下结论：
1) 湿地对局地温度影响较大，能较为明显地对研究区域起到降温作用，最大降温幅度达到了7.1℃。
2) 湿地有利于东南风风道的加强，最大风速增加幅度达到了1.9 m/s，增强了研究区域的通风能力，并进而引起研究区域内的温度调节。
3) 三垟湿地起到了有效改善城市气候、减弱城市热岛效应和保护大气环境的作用。
本文所得到的结论仅是基于一个个例的模拟结果，尚需进一步的观测和模拟研究加以验证。使用的MODIS下垫面资料与实际地表类型和植被覆盖还有一些差异，并且WRF模式也不完全适用于局地高分辨率数值模拟。这一切都需要在今后的研究工作中解决。
基金项目
温州市气象局项目(2018ZD01)资助。