南京地區下墊面變化對城市熱島效應影響的數值模擬

葉麗梅,江志紅,霍飛

(1.氣象災害教育部重點實驗室(南京信息工程大學),江蘇 南京 210044;2.武漢區域氣候中心,湖北 武漢 430074;3.南京信息工程大學 大氣科學學院,江蘇 南京 210044)

南京地區下墊面變化對城市熱島效應影響的數值模擬

葉麗梅1,2,江志紅1,3,霍飛1

(1.氣象災害教育部重點實驗室(南京信息工程大學),江蘇 南京 210044;2.武漢區域氣候中心,湖北 武漢 430074;3.南京信息工程大學 大氣科學學院,江蘇 南京 210044)

利用新一代中尺度研究和預報模式(Weather Research and Forecasting Model,簡稱 WRF)分別耦合多冠層、單冠層和平板模式三種情況進行南京地區2007年8月1日的天氣過程模擬,分析不同城市冠層方案對南京氣象場的模擬效果。在此基礎上,結合模擬效果最好的城市冠層方案,研究南京城市下墊面的變化對其熱島的影響。結果表明:多冠層方案對近地面氣溫、10 m風場的模擬效果最好;城區的擴張使南京地區近地面氣溫升高,主要表現為城市區域夜間升溫顯著,并且導致熱島強度明顯增強;城市擴張后,城區白天風速大范圍地減小,同時熱島環流更加顯著,且具有明顯的城市熱島的“下游效應”。

下墊面變化;城市熱島;WRF;數值模擬

0 引言

近年來,我國東部大規模城市化帶來了土地利用/土地覆蓋的快速變化,南京作為長三角的重要城市,城市建設面貌發生了巨大的變化。2000年以前南京市區面積占南京市土地總面積不到20%,2003年后南京市區面積達到南京市土地總面積的70%之多,城市發展處于極度擴張期(呂憲軍和王梅,2006)。這種大規模下墊面類型的快速變化,必然對城市的天氣和氣候帶來重要影響,并伴隨著夏季高溫熱浪、城市積澇等災害性氣象事件頻發。深入研究快速城市化的直接氣候效應及其形成機理,可以為合理規劃城市發展規模、布局,盡可能減少城市化造成的負面影響提供科學依據。

城市熱島(Urban Heat Island,UHI)作為城市氣候效應主要特征,一直受到氣候學家的廣泛關注。以往UHI的研究主要通過觀測分析和數值模擬兩種途徑(黃利萍等,2012;苗峻峰,2014)。由于氣象站點的分布通常較為稀疏,且主要為地面觀測,觀測分析一般無法給出熱島的三維結構,無法了解具體的物理過程。目前數值試驗被廣泛地用于城市效應的研究當中,如宋靜等(2009)將WRF模式與城市冠層模式(Urban Canopy Model,UCM)耦合,結合同期外場試驗觀測結果和模式模擬結果,分析城市冠層效應對局地天氣的影響;陳燕和蔣維楣(2007)運用數值模擬手段研究城市化進程對城市邊界層結構的影響,得到城市反照率減小、植被減少、地表濕度減小,蒸發耗熱減少、感熱通量增加、城市波恩比增加等結論。雖然這些研究使用了不同時期的下墊面資料以分別反映城市化前后的陸面狀況,但是其結論一般基于某一種UCM。不同的UCM在某一地區的模擬能力是有區別的,因此有必要對多種UCM的模擬能力進行評估。近年來,國際上先后建立和發展出數十種UCM,有些已經與中尺度天氣模式進行了耦合。最近發布的中尺度WRF模式,V3中能夠耦合多冠層、單冠層和平板3種UCM,目前被廣泛地運用于城市天氣與氣候的研究中。

本文將WRF模式分別耦合多冠層、單冠層和平板3類UCM,對2007年8月1日南京地區氣象場進行模擬試驗,并比較3類UCM的模擬效果。在此基礎上,使用最優的UCM,考察不同時期土地利用變化對南京地區地面氣溫、熱島強度、熱島環流等氣象場的影響。

1 WRF模式及UCM模式簡介

WRF模式是由美國國家大氣研究中心、美國國家環境預測中心和俄克拉荷馬大學暴雨分析預報中心等多家單位聯合發展起來的新一代非靜力平衡、高分辨率、科研和業務預報統一的中尺度預報和資料同化模式。第三版較之前的版本主要增加了近期MODIS下墊面資料,同時UCM中增加了多層城市冠層方案,可以更加全面地考慮城市冠層的物理過程。下面簡單介紹本文使用的UCM。

1)平板模式(Surface-Layer Scheme),平板模式將建筑作為平坦的、具有一定厚度的小塊介質層處理,城市在模式中為裸土或平板,考慮了城市表面的物理屬性,如:熱容量、熱傳導、反射率、粗糙度等,沒有考慮城市的幾何形狀。

2)單層城市冠層模式(Kusaka et al.,2001;Kusaka and Kimura,2004;Tewari et al.,2004)(Single-Layer Urban Canopy Model),簡稱單冠層模式,單冠層模式計算了墻面、屋頂、道路產生的熱通量,在表面能量收支平衡和風切變環流中考慮城市幾何形狀的影響。

3)多層城市冠層模式(Kondo and Liu,1998;Kondo et al.,1999;Vu et al.,1999)(Multi-Layer Urban Canopy Model),簡稱多冠層模式,多冠層模式中將地面到屋頂之間分為若干層,并按這些層次分層計算屋頂、墻面、道路的能量收支。它對城市特征的描述更準確,更接近真實,甚至考慮了各個格點上建筑物的高度和分布密度的不同。

2 數值試驗方案與資料

試驗采用雙向反饋的四重嵌套網格,模擬中心點為118.192°E、30.006°N,水平格距分別為27、9、3和1 km(圖1a)。四層區域的格點數分別為66×66、147×147、180×189和111×177,最內層的區域及地形高度如圖1b所示。

初始場選用2007年7月30日20時—8月2日02時(北京時間,下同)的NCEP/NCAR 1°×1° FNL資料,6 h更新一次側邊界,Spin-up時間為4 h。500 hPa風場(圖略)分析表明,該時段南京上空盛行西風,天氣形勢穩定。模式參數化方案包含:Ferrier微物理方案,RRTM長波輻射方案,Dudhia短波輻射方案,Noah陸面過程方案,Eta MYJ (Mellor-Yamada-Janjic)TKE邊界層方案,Kain-Fritsch積云對流方案。

圖1 模擬區域及地形高度(單位:m) a.模式模擬區域;b.最內層模擬區域Fig.1 Modeling domains and terrain elevation(units:m) a.nested domains;b.inner domain

圖2 不同時期南京地區的土地利用情況(黑色矩形框表示主要城市和建筑用地;白色虛、實線分別表示計算熱島強度時選取的市區和郊區的所在位置) a.1993年;b.2001年Fig.2 Land use of Nanjing in (a)1993 and (b)2001 (black rectangles denote urban and construction land;white solid and dashed rectangles represent the location of urban and rural areas respectively in calculation of heat island intensity)

WRF模式自帶的兩種土地利用資料,分別是美國地質調查局分辨率30″的格點資料(圖2a)和Boston University/NCEP提供的MODIS 30″的格點資料(圖2b),前者是1993年土地利用狀況,可以代表大規模城市化前的下墊面狀況;后者是2001年的土地利用狀況,可代表大規模城市化后的下墊面狀況。對比該區域2007年MODIS資料的下墊面類型(圖略)與模式自帶2001年的下墊面類型,發現城市面積沒有較大變化。這說明模式自帶的2001年下墊面類型可以代表2000年后的情況。

為了研究下墊面的變化對南京UHI的影響,本文首先通過模擬結果和觀測資料的比較選擇最優城市冠層方案,然后通過耦合最優冠層方案研究城市發展對UHI所產生的影響。

城市冠層耦合方案模擬選擇試驗:使用2001年MODIS土地利用類型資料,分別耦合多冠層、單冠層、平板3種UCM,記為:1)m-multi試驗;2)m-single試驗;3)m-slab試驗。

下墊面變化影響的對比試驗:1)利用1993年USGS土地利用類型資料,耦合最優城市冠層方案進行模擬試驗,稱為USGS試驗;2)使用MODIS 2001年的土地利用類型資料,耦合最優城市冠層方案進行試驗,稱為MODIS試驗。

圖3 2007年8月1日00—24時觀測站近地面氣溫日變化的模擬值與觀測值(單位:℃) a.南京站;b.六合站;c.江浦站;d.江寧站Fig.3 Simulated and observed temperature(℃) at different stations between 00:00 BST and 24:00 BST on 1 August 2007 a.Nanjing;b.Luhe;c.Jiangpu;d.Jiangning

3 不同城市冠層方案對南京城市氣象場模擬效果分析

3.1 溫度場的模擬效果分析

根據模式自帶的2001年MODIS下墊面資料,南京站、六合站為城市站,江浦站、江寧站為非城市站,圖3給出2007年8月1日00—24時此4站觀測氣溫與模擬值的對比。可以看到城市站(圖3a、b)3種參數化方案模擬的氣溫差別較大,m-multi試驗的模擬夜間氣溫高于m-single試驗和m-slab試驗,與觀測氣溫最為接近,白天氣溫則低于m-single試驗和m-slab試驗;在非城市站(圖3c、d)3種方案模擬的氣溫比較一致,模擬氣溫都低于觀測值;但總體m-multi試驗氣溫日變化曲線與觀測值日變化曲線的形狀最為吻合。表1給出各測站模擬氣溫的均方誤差,可以看出:各測站m-multi試驗的均方誤差較小,其均方誤差的平均值也是最小的。

表1近地面氣溫模擬結果與觀測值的均方誤差

Table 1 Standard deviations of simulated and observed temperatures ℃

觀測站m-multi試驗m-single試驗m-slab試驗南京1.111.371.34六合0.940.881.03江浦2.112.262.23江寧1.341.621.56平均1.381.531.54

進一步分析南京地區地表溫度空間分布的模擬效果。圖4d是基于MODIS衛星數據劈窗算法反演得到的地表溫度(丁莉東等,2005),反映出強熱島高值區被長江割裂為南北兩部分,其中最高值出現在32°N的主城區中心。城區的地表溫度比郊區高約8 ℃。圖4a—c給出的是不同UCM模擬的地表溫度:m-multi試驗(圖4a)模擬的城郊溫差約為8 ℃,主城區中心對應最高值,熱島結構與衛星反演形態上十分接近;m-single試驗(圖4b)和m-slab試驗(圖4c)的城郊溫差約為5～8 ℃,但主城區的模擬溫度分布均一,沒有模擬出城市中心的最高值。由此可見,無論是對氣溫的日變化還是空間分布,多冠層模式的模擬效果較單冠層和平板模式好。

圖4 2007年8月1日11時m-multi試驗(a)、m-single試驗(b)、m-slab試驗(c)的地表溫度模擬值和MODIS衛星反演值(d)(單位:℃)Fig.4 Simulated surface temperature by (a)m-multi test,(b)m-single test,(c)m-slab test,and (d)MODIS surface temperature at 11:00 BST on 1 August 2007(units:℃)

3.2 風場的模擬效果分析

圖5為4個測站10 m風速在2007年8月1日00—24時日變化的模擬與觀測對比。由圖可見,城市站(圖5a、b)3種參數化方案模擬的風速差別較大,其中夜間m-multi試驗模擬的風速與觀測值接近,m-single和m-slab試驗模擬值偏高;白天3組試驗的模擬值與實際值相差都較大。對于非城市觀測站(江浦站、江寧站),3組試驗的風速模擬值與觀測值比較接近,3組試驗間的差異不明顯。

表2給出了10 m風場不同方案均方差的對比,可以看出,m-multi試驗模擬的風速值均方差最小。3組試驗風向的均方差都較大,試驗間差別不明顯。

由3組試驗在2007年8月1日12時和20時模擬的10 m風場水平分布與臺站觀測資料對比(圖6),可以發現觀測風速(圖6a、e)在32°N的城市中心表現為低值區。m-multi試驗(圖6b、f)能很好地模擬出城市中心風速的最低值,城郊交界處有著明顯的風速梯度。m-single試驗(圖6c、g)和m-slab試驗(圖6d、h)不僅沒有模擬出城市中心的低值區,而且城區風速還略大于郊區。3組試驗對于風向的模擬差異不大,與觀測值大致吻合。

表210m風場模擬結果與觀測值的均方差

Table 2 Standard deviations of simulated and observed surface wind

m-multi試驗m-single試驗m-slab試驗風速/(m·s-1)風向/(°)風速/(m·s-1)風向/(°)風速/(m·s-1)風向/(°)南京1.1826.021.3426.741.6426.22六合1.0832.201.6036.341.7534.82江浦1.0032.931.1737.091.1937.69江寧0.9858.901.1745.741.1445.44平均1.0637.511.3236.481.4336.04

圖6 2007年8月1日12時(a—d)和20時(e—h)10 m風速的觀測值與模擬值(單位:m·s-1) a,e.觀測;b,f.m-multi試驗;c,g.m-single試驗;d,h.m-slab試驗Fig.6 Observed and simulated surface wind speed at (a—d)12:00 BST and (e—h)20:00 BST on 1 August 2007 (units:m·s-1) a,e.observation;b,f.m-multi;c,g.m-single;d,h.m-slab

圖7 2007年8月1日00時—2日00時MODIS試驗與USGS試驗中近地面氣溫差值日變化(單位:℃) a.非城市站發展為城市站;b.城市站發展為非城市站Fig.7 Daily temperature difference (℃) of MODIS test and USGS test between 00:00 BST on 1 August 2007 and 00:00 BST on 2 August 2007 a.stations changing from rural to urban;b.stations changing from urban to rural

綜上比較可以看到,相比單冠層和平板模式,多冠層模式對城市氣溫的模擬效果較好,對于風向的模擬優勢不明顯,從均方差來看,多冠層模式對風速的模擬較優。

4 下墊面變化對南京UHI影響的模擬研究

前述試驗表明,WRF模式耦合多冠層模式能較好地模擬南京地區氣象場。為進一步研究城區面積擴展對南京UHI的影響,以2007年8月1日為代表,進行下墊面變化對比試驗(USGS試驗、MODIS試驗),分析不同的城市化規模對近地面氣溫、熱島強度、風場和熱島環流的影響。

4.1 對近地面氣溫的影響

根據USGS、MODIS的下墊面資料,可以得到南京站、溧水站、江浦站3個觀測站在不同時期的土地利用類型(表3),其中南京、溧水站由非城市站發展為城市站,江浦站由城市站變為非城市站。圖7給出了MODIS試驗與USGS試驗近地面氣溫差值的日變化。對于非城市站發展為城市站的觀測站(南京站、溧水站),下墊面的變化使夜間升溫明顯高于白天,這是由于城市冠層與自然下墊面不同,白天對輻射的吸收和攔截作用增大,夜間將增大的儲存熱量釋放。對于由城市站變為非城市站(江浦站)的觀測站,MODIS試驗與USGS試驗相比,溫度沒有明顯上升,在某些時刻(00—02時)MODIS試驗模擬的氣溫反而偏低,可見城市化主要表現出夜間UHI的顯著增強現象。

進一步分析溫度變化的空間分布,夜間02時(圖8a),下墊面的變化使南京城區近地面氣溫升高,城市擴張區升高2 ℃以上,而郊區沒有明顯的升溫;14時(圖8b),南京地區氣溫差異普遍不顯著,再次表明城市化主要導致夜間城區顯著升溫,而白天升溫并不明顯。

表33個觀測站在不同試驗中的土地利用類型

Table 3 Land use of 3 stations in different experiments

測站土地利用類型USGSMODIS南京灌溉農田和牧場城市和建筑用地溧水耕地城市和建筑用地江浦城市和建筑用地耕地

4.2 對熱島強度的影響

圖8 02時(a)和14時(b)MODIS試驗與USGS試驗中近地面氣溫差值的水平分布(單位:℃)Fig.8 Horizontal distribution of surface air temperature difference(℃) in MODIS test and USGS test at (a)02:00 BST and (b)14:00 BST

圖9 2007年8月1日00—24時MODIS試驗與USGS試驗中熱島強度日變化(單位:℃)Fig.9 Daily difference of Heat Island intensity(℃) in MODIS test and USGS test between 00:00 BST and 24:00 BST on 1 August 2007

選取城市地區(118.76～118.80°E,32.00～32.05°N,圖2中白色虛線框)與郊區(118.90～118.94°E,31.90～31.95°N,圖2中白色實線框),計算地面2 m高度上氣溫差值的時間序列,并將這種城鄉差異定義為熱島強度。為了保證城區與郊區的典型性,所選城區與郊區的土地利用類型在兩組試驗中是一致的。圖9給出了不同下墊面下熱島強度的日變化對比,可以看到,兩組試驗都表現出熱島強度夜間強于白天。MODIS試驗模擬的熱島強度全天大于USGS試驗。MODIS試驗的熱島強度最大值出現在02時,達3.4 ℃;USGS試驗的熱島強度最大值出現在01時,達到2.8 ℃。圖中虛線標出了兩組試驗熱島強度高于1 ℃時所對應的時間段:USGS試驗對應的是06時之前與17時之后,而MODIS試驗是08時之前與15時之后。這說明,隨著城市面積的擴展,強熱島持續時間變長了。結合上一節的研究,城市化導致南京城區夜間升溫顯著,表明城市的擴張使夜間的熱島強度顯著增強,持續時間明顯延長。

4.3 對10 m風速的影響

與自然表面不同,城市表面幾何結構復雜,粗糙度較大。圖10是MODIS試驗與USGS試驗模擬的2007年8月1日00—24時10 m風速差的日變化,可以看到,對于發展為城市的觀測站(南京、溧水),全天風速差為負值,風速差異最大值出現在12時左右;而對于由城市站變為非城市站的江浦,由于缺少城市冠層的摩擦耗散,風速是明顯增加的,差異最大值出現在15時。風速差異極值總是出現于正午至傍晚之間,說明大規模城市化對白天風速的影響最明顯。

圖11是MODIS試驗與USGS試驗中10 m水平風速差值的空間分布對比,可以清楚地看到,無論是12時還是20時,32°N所在城市中心及其周邊的區域都表現為負值區,這表明城市擴張后城區的風速會明顯減小。

圖10 2007年8月1日00—24時MODIS試驗與USGS試驗中10 m風速差值的日變化(單位:m·s-1) a.非城市站發展為城市站;b.城市站發展為非城市站Fig.10 Daily difference of surface wind speed(m·s-1) in MODIS test and USGS test between 00:00 BST and 24:00 BST on 1 August 2007 a.stations changing from rural to urban;b.stations changing from urban to rural

圖11 12時(a)和20時(b)MODIS試驗與USGS試驗中10 m水平風速差值的空間分布 (單位:m·s-1)Fig.11 Distribution of 10m surface wind speed difference in MODIS test and USGS test at (a)12:00 BST and (b)20:00 BST (units:m·s-1)

4.4 對熱島環流的影響

以往對UHI的研究,因為缺少高分辨率三維數據,很少涉及城市熱島環流特征分析。WRF模式能模擬UHI的垂直結構,可以加深了解UHI的局地環流特征。

考察城市擴張對城區垂直環流所產生的影響,考慮到模擬時間內近地面為大范圍的西風氣流,因而選取沿東西向通過城市中心的剖面分析其環流特征。圖12分別是兩組試驗8月1日13時位溫場及風場經過32.065°N沿東西向的垂直剖面。可以看到,USGS試驗模擬的熱島環流不明顯,而MODIS試驗中出現了典型的熱島環流。在MODIS試驗中,城市(紅色橫條)表面對應著高溫中心,空氣受熱上升,考慮到低層風場為西風氣流,可見上升運動最強支出現在城區的下游。

圖13給出了2007年8月1日13時MODIS試驗與USGS試驗沿32.065°N東西向的位溫差垂直剖面及MODIS試驗的垂直環流。可以看到,隨著城區下墊面擴張,大部分地區表現為氣溫上升,并且增溫的最大值出現在城市東部上空。上升運動高值區與氣溫高值區位置比較吻合(等值線為MODIS試驗的垂直速度),出現在城市的東部。結合低層西風氣流,表明城市擴張后城區下游的UHI更加明顯,并且在城市下游出現了強的熱島環流上升支。上述結果說明南京地區UHI存在明顯的“下游效應”。Zhang et al.(2009)的研究發現Baltimore上游的城市化會加劇下游的UHI效應,并且提出UHI的“下游效應”。本文的數值試驗結果表明:UHI的“下游效應”在南京也同樣存在。

圖12 2007年8月1日13時USGS試驗(a)和MODIS試驗(b)與沿32.065°N東西向的位溫(單位:℃)及風場(單位:m·s-1)垂直剖面(紅條:城市和建筑用地;藍條:長江;綠條:紫金山)Fig.12 Vertical profile of potential temperature(℃) and wind field(m·s-1) in (a)USGS test and (b)MODIS test along 32.065°N latitude at 13:00 BST on 1 August 2007 (red bar:urban and construction land;blue bar Yangtze River;green bar:Purple Mountain)

圖13 2007年8月1日13時MODIS試驗與USGS試驗沿32.065°N東西向的位溫差(單位:℃)垂直剖面及MODIS試驗的風場(單位:m·s-1)(陰影區為位溫差;等值線為MODIS試驗的垂直速度;紅條、藍條、綠條分別為MODIS下墊面下城市和建筑用地、長江、紫金山)Fig.13 Vertical profile of potential temperature difference (℃) in MODIS test and USGS test and wind field(m·s-1) in MODIS test along 32.065°N latitude at 13:00 BST on 1 August 2007 (shaded area:potential temperature difference;contour:vertical velocity in MODIS test;red bar:urban and construction land;blue bar:Yangtze River;green bar:Purple Mountain)

綜上數值試驗研究結果表明,城區的擴張使南京地區近地面氣溫升高,熱島強度明顯增強,夜間強熱島維持時間顯著延長;城市化后城區白天的風速大幅減小;熱島環流更加顯著,且具有明顯的UHI“下游效應”。

5 結論與討論

本文利用WRF模式分別耦合多冠層、單冠層和平板模式,模擬南京地區2007年8月1日的天氣過程,分析不同城市冠層方案對南京氣象場的模擬效果。在此基礎上,結合模擬效果最好的城市冠層方案,研究南京城市下墊面的變化對其熱島的影響。通過比較3種UCM的模擬效果,發現相比單冠層和平板模式,多冠層模式的氣溫日變化曲線與觀測值日變化曲線的形狀最為吻合。從均方差上看,多冠層模式對風速的模擬較優。利用多冠層模式模擬下墊面變化對UHI影響,結果表明:城區的擴張使南京地區近地面氣溫升高,且主要表現為城市區域夜間升溫明顯,熱島強度明顯增強;城區白天風速大幅減小,熱島環流更加顯著,且表現出明顯的UHI“下游效應”。

由于WRF模式的下墊面資料只有1993年的USGS資料和2001年MODIS資料,本文對于2007年各氣象要素場的模擬使用的是2001年MODIS資料。對比該區域2007年MODIS資料的下墊面類型與模式自帶2001年的下墊面類型,發現城市面積沒有較大變化。這說明模式自帶的2001年下墊面類型可以代表2000年后的情況。因為冬季熱島效應所反映的特征和夏季并不相同,所以模擬試驗只選擇了夏季某一天與實測資料對比,這造成本文的結論在一定程度上存在局限性,在以后的研究中可以通過增加對于不同年份、季節的模擬,增強結論的普遍性。

陳燕,蔣維楣.2007.南京城市化進程對大氣邊界層的影響研究[J].地球物理學報,50(1):66-73.

丁莉東,覃志豪,毛克彪.2005.基于MODIS影像數據的劈窗算法研究及其參數確定[J].遙感技術與應用,20(2):284-289.

黃利萍,苗峻峰,劉月琨.2012.天津城市熱島效應的時空變化特征.大氣科學學報,35(5):620-632.

呂憲軍,王梅.2006.行政區劃調整與城市擴張研究——以南京市為例[J].現代城市研究,6(1):67-72.

苗峻峰.2014.城市熱島和海風環流相互作用的數值模擬研究進展.大氣科學學報,37(4):521-528.

宋靜,湯劍平,孫鑒濘.2009.南京地區城市冠層效應的模擬試驗研究[J].南京大學學報:自然科學版,45(6):779-789.

Kondo H,Liu F H.1998.A study on the urban thermal environment obtained through one-dimensional urban canopy model[J].J Japan Soc Atmos Environ,33:179-192.

Kondo H,Kikegawa Y,Genchi Y,et al.1999.Heating in the urban canopy by anthropogenic energy use[C]//Proceedings of 15th International Congress of Biometeorology and International Conference on Urban Climatology(ICB-ICUC 99).Sydney,Australia.

Kusaka H,Kimura F.2004.Coupling a single-layer urban canopy model with a simple atmospheric model:Impact on urban heat island simulation for an idealized case[J].J Meteor Soc Japan,82:67-80.

Kusaka H,Kondo H,Kikegawa Y,et al.2001.A simple single-layer urban canopy model for atmospheric models:Comparison with multi-layer and slab models[J].Bound-Layer Meteor,101:329-358.

Tewari M,Chen F,Wang W,et al.2004.Implementation and verification of the unified Noah land surface model in the WRF model[C]//20th Conference on Weather Analysis and Forecasting/16th Conference on Numerical Weather Prediction.Seattle.

Vu T C,Asaeda T,Ashie Y.1999.Development of a numerical model for the evaluation of the urban thermal environment[J].J Wind Eng Ind Aerodyn,81:181-191.

Zhang D,Shou Y,Russe D.2009.Upstream urbanization exacerbates Urban Heat Island effects[J].Geophys Res Lett,36,L24401.doi:10.1029/2009GL041082.

(責任編輯：劉菲)

NumericalsimulationoftheimpactoflandcoverchangeontheurbanheatislandeffectinNanjing

YE Li-mei1,2,JIANG Zhi-hong1,3,HUO Fei1

(1.Key Laboratory of Meteorological Disaster(NUIST),Ministry of Education,Nanjing 210044,China;2.Wuhan Regional Climate Center,Wuhan 430074,China;3.School of Atmospheric Sciences,NUIST,Nanjing 210044,China)

In this paper,the Weather Research and Forecast Model(WRF) is coupled with Surface-Layer Scheme,Single-Layer Urban Canopy Model and Mingle-Layer Urban Canopy Model respectively to evaluate the simulation effect of various parameterizations on the weather conditions on 1 August 2007 in Nanjing.The best urban parameterization scheme is coupled into WRF to study the impact of land cover change on the Urban Heat Island(UHI) effect in Nanjing.Results show that the Mingle-Layer Urban Canopy Model shows the best simulation effect for surface temperature and 10m wind field.Urbanization makes surface air temperature increase over the region,especially at night and thus intensifies the UHI effect.After urbanization,the wind speed in the downtown area decreases obviously while the Urban Heat Circulation occurs more apparently.There also exists the downstream effect of UHI in Nanjing.

land cover change;Urban Heat Island;WRF;numerical simulation

2012-10-17;改回日期2013-01-13

國家重點基礎研究發展規劃項目(2010CB428505)

江志紅,教授,博士生導師,研究方向為短期氣候診斷,zhjiang@nuist.edu.cn.

10.13878/j.cnki.dqkxxb.20121017006.

1674-7097(2014)05-0642-11

P404

A

10.13878/j.cnki.dqkxxb.20121017006

葉麗梅,江志紅,霍飛.2014.南京地區下墊面變化對城市熱島效應影響的數值模擬[J].大氣科學學報,37(5):642-652.

Ye Li-mei,Jiang Zhi-hong,Huo Fei.2014.Numerical simulation of the impact of land cover change on the urban heat island effect in Nanjing[J].Trans Atmos Sci,37(5):642-652.(in Chinese)