苗峻峰

(南京信息工程大學(xué) 大氣科學(xué)學(xué)院, 江蘇 南京 210044)

城市熱島和海風環(huán)流相互作用的數(shù)值模擬研究進展

苗峻峰

(南京信息工程大學(xué) 大氣科學(xué)學(xué)院, 江蘇 南京 210044)

在沿海城市地區(qū),城市熱島環(huán)流和海風環(huán)流往往同時存在;它們在空氣污染物的傳輸和擴散中均起至關(guān)重要的作用,對強對流天氣亦有重要的觸發(fā)作用。隨著城市化的發(fā)展,空氣污染等環(huán)境問題變得越來越嚴峻,突然性強對流天氣所造成的災(zāi)害也越來越嚴重。因此,沿海城市地區(qū)城市熱島和海風環(huán)流相互作用的研究受到日益廣泛的關(guān)注。本文回顧了近30 a來國內(nèi)外關(guān)于城市熱島和海風環(huán)流相互作用數(shù)值模擬的研究歷史,分析了目前的研究現(xiàn)狀及存在的一些問題,概述了城市化、城市熱島對海風環(huán)流的影響,海風對城市熱島、城市熱島環(huán)流的影響,以及城市熱島環(huán)流和海風環(huán)流的耦合效應(yīng)。最后,本文提出了一些有待于研究或需深入研究的問題;這些問題的研究將有助于進一步揭示沿海城市地區(qū)空氣污染動力學(xué)機制、強對流天氣觸發(fā)機制。

熱島環(huán)流;海陸風;湖陸風;局地環(huán)流;局地對流;城市效應(yīng);空氣污染動力學(xué)

0 引言

21世紀以來,隨著城市化的發(fā)展,城市氣象學(xué)(urban meteorology)和城市氣候?qū)W(urban climatology)的研究已越來越受到國內(nèi)外廣大氣象學(xué)家的重視和政府部門的強有力支持(Oke,2006;Grimmond et al.,2010;壽亦萱和張大林,2012)。國際城市氣候協(xié)會(International Association for Urban Climate,IAUC)的創(chuàng)立(2000年)及每3 a一次的國際城市氣候會議(International Conference for Urban Climate,ICUC)就是一個重要標志。

城市熱島(urban heat island,UHI)是城市氣侯最顯著的特征之一(Kanda,2007;黃利萍等,2012),其形成原因主要是城市下墊面性質(zhì)的改變與城市人為熱源的排放。大量的研究表明:城市化、城市熱島對空氣污染(徐祥德等,2006)、大氣邊界層(陳燕和蔣維楣,2007)、(強)對流天氣和降水(Shepherd,2005;Guo et al.,2006;蒙偉光等,2007;Chen et al.,2007;Zhang et al.,2009;Lin et al.,2011;Miao et al.,2011;徐蓉等,2013)均有重要的影響;城市熱島強度與城市特征(城市規(guī)模、建筑物平均高度及密度、地表物理屬性、街區(qū)走向及天空可視因子、人口密度、人為熱源的排放量等)和氣象條件密切相關(guān)(Atkinson,2003;蔣維楣等,2009,2010)。

海風環(huán)流(sea breeze circulation,SBC)是沿海地區(qū)重要的中尺度氣象現(xiàn)象之一(張立鳳等,1999;Lin et al.,2001;Miao et al.,2003,2009;Miller et al.,2003;Crosman and Horel,2010;許啟慧等,2013a),主要是由于海陸之間的熱力差異所造成的。海風環(huán)流在空氣污染物的傳輸和擴散中起重要的作用(陳訓(xùn)來等,2008),海風鋒能誘發(fā)(強)對流天氣的發(fā)生(王樹芬,1990;Fovell,2005;尹東屏等,2010;盧煥珍等,2012;易笑園等,2012;汪雅等,2013)。大量的研究表明:海風非常敏感于實際地理環(huán)境的復(fù)雜性、地面溫度和大尺度背景風(Crosman and Horel,2010)。因此,地面的不均勻性對海風環(huán)流的發(fā)生和發(fā)展以及環(huán)流強度有至關(guān)重要的影響。

眾所周知,許多大城市位于沿海地區(qū),因而城市冠層(urban canopy)在局地氣象學(xué)(local meteorology)中通常起關(guān)鍵的作用,尤其是當城市熱島效應(yīng)產(chǎn)生城市熱島環(huán)流(urban heat island circulation,UHIC)(Vukovich et al.,1979;周明煜等,1980;Savij?rvi,1985;Eliasson and Holmer,1990;Han and Baik,2008)時更是如此。城市熱島是城市熱島環(huán)流的源,城市熱島環(huán)流是大氣對與感熱通量梯度有關(guān)的溫度水平變化的中尺度響應(yīng)(Hidalgo et al.,2010)。城市熱島環(huán)流的形成原因類似于海風環(huán)流,后者主要是由于海陸之間的熱力差異所造成的,而前者主要是由于城鄉(xiāng)之間的熱力差異(和城市大小)所造成的。城市熱島環(huán)流的范圍可以是城市大小的2～3倍(Hidalgo et al.,2008),屬城市尺度(urban scale)或中小尺度環(huán)流,而海風環(huán)流是局地或區(qū)域中尺度環(huán)流。

在沿海城市地區(qū),尤其是大城市地區(qū),城市熱島環(huán)流和海風環(huán)流往往同時存在(黃利萍等,2013;許啟慧等,2013b;Meir et al,2013),因而它們是相互作用的(Cenedese and Monti,2003;Childs and Raman,2005)。也就是說,城區(qū)的存在影響海風環(huán)流的結(jié)構(gòu)和發(fā)展;反過來,海風環(huán)流控制城區(qū)及周圍地區(qū)溫度、濕度和空氣污染物濃度的時空變化。隨著城市化的發(fā)展,空氣污染等環(huán)境問題變得越來越嚴峻,城市熱島和(或)海風鋒誘發(fā)的強對流天氣所造成的災(zāi)害也越來越嚴重(東高紅等,2013)。因此,城市熱島和海風環(huán)流的相互作用受到科學(xué)家們越來越多的關(guān)注。本文在分析了20世紀80年代以來國內(nèi)外相關(guān)研究成果的基礎(chǔ)上,概述了城市熱島和海風環(huán)流相互作用的數(shù)值模擬研究進展,并提出了一些有待于深入研究的問題。

1 歷史回顧與現(xiàn)狀

早在20世紀80年代,科學(xué)家們就開始關(guān)注城市熱島效應(yīng)與海風環(huán)流之間的相互作用問題。30多年來,取得了不少研究成果。研究方法以數(shù)值試驗(numerical experiment)和數(shù)值模擬(numerical simulation、numerical modelling)為主,研究內(nèi)容涉及城市化、城市熱島對海風環(huán)流的影響(Schultz and Warner,1982;Savij?rvi,1985;Yoshikado,1992,1994;張雷鳴等,1994;Sarkar et al.,1998;苗曼倩和唐有華,1998;Kusaka et al.,2000;Ohashi and Kida,2002a,2002b,2002c;李維亮等,2003;Ferretti et al.,2003;Martilli,2003;Lemonsu et al.,2006;Freitas et al.,2007;Lo et al.,2007;Lin et al.,2008;Simpson et al.,2008;陸希和壽紹文,2009;文偉俊等,2009;Dandou et al.,2009;Lu et al.,2010;Chen et al.,2011a;梁釗明等,2013;張亦洲等,2013;Ryu and Baik,2013),海風對城市熱島、城市熱島環(huán)流的影響(Schultz and Warner,1982;Yoshikado,1992;Ohashi and Kida,2002a,2002b;Chemel and Sokhi,2012),以及城市熱島環(huán)流—海風環(huán)流的耦合效應(yīng)(Childs and Raman,2005;Freitas et al.,2007;劉樹華等,2009)。數(shù)值模擬研究區(qū)域涉及日本東京、大阪—京都,意大利羅馬—奧斯蒂亞,美國紐約、休斯頓,法國馬賽,巴西圣保羅,印度欽奈,希臘雅典,英國倫敦,韓國首爾,中國京津冀、長江三角洲、珠江三角洲、廣西北部灣、臺灣地區(qū)等。

20世紀80年代,Schultz and Warner(1982)使用二維原始方程模式模擬了美國洛杉磯盆地夏天的局地環(huán)流,通過數(shù)值試驗探討了海風、城市熱島和山谷風對“氣流型”(airflow pattern)的貢獻。Savij?rvi(1985)使用二維干靜力模式研究了海陸風、城市熱島環(huán)流及其它們的相互作用,發(fā)現(xiàn)熱的城市大大加強了海風,并把海風“鎖”在海岸線附近。

20世紀90年代,Yoshikado(1992)使用二維靜力邊界層模式(boundary-layer model),針對日本東京的情況,模擬了城市效應(yīng)與海風的相互作用,發(fā)現(xiàn)沿海區(qū)域的城市熱島在海風發(fā)展階段起加強海風速度的作用;海風鋒傳播的停滯是由海風環(huán)流和熱島環(huán)流的相互作用引起的,因此在海風鋒前部產(chǎn)生了一個弱風區(qū)。Yoshikado(1994)使用同樣的二維邊界層模式,研究了海風和城市熱島相互作用的程度對城區(qū)大小(寬度)、城區(qū)離開海岸線的距離以及熱島強度的依賴性,發(fā)現(xiàn)當城區(qū)寬度超過10 km時,海風和城市熱島的相互作用變得更為重要;對于較小的城市區(qū)域而言,熱島發(fā)展得不夠強,相互作用不明顯;當城區(qū)位于離海岸線的一定距離內(nèi)時,海風和熱島的相互作用變得比較強,這是由于在海風到達城區(qū)之前,熱島在很大程度上已經(jīng)發(fā)展起來了。Sarkar et al.(1998)用二維數(shù)值模式研究了沿海城市熱島對海風環(huán)流特征的影響,討論了城區(qū)寬度和人為熱通量的效應(yīng),結(jié)果表明:在下午早些時候,城市熱島的存在加強了海風鋒輻合的強度,減少了海風鋒向內(nèi)陸傳播的距離;較大的城市寬度能導(dǎo)致較大的垂直速度,較大的熱島強度會導(dǎo)致較強的輻合。苗曼倩和唐有華(1998)應(yīng)用三維中尺度模式研究了長江三角洲夏季海陸風與熱島環(huán)流的相互作用,發(fā)現(xiàn)白天由于東海(East China Sea)海風和太湖湖風環(huán)流與上海市城市熱島相互增強,最大垂直速度可達6.2 cm/s;上海市因海風的調(diào)節(jié)作用,致使城市熱島環(huán)流強度比內(nèi)陸城市弱1～2 ℃。

進入21世紀以來,隨著城市化的發(fā)展,城市熱島和海風環(huán)流相互作用的研究得到更廣泛的關(guān)注。Kusaka et al.(2000)使用三維局地環(huán)流模式對日本東京都市區(qū)(metropolitan area)進行的土地利用(land use)變化敏感性試驗研究發(fā)現(xiàn)城市發(fā)展和下墊面改變使得海風鋒抵達內(nèi)陸所需的時間增加。Ohashi and Kida(2002a)使用三維中尺度大氣模式研究了兩個理想相鄰的城區(qū)(一個在沿海、一個在內(nèi)陸)對局地環(huán)流的影響,發(fā)現(xiàn)當海風鋒向內(nèi)陸移動時,會形成一個“鏈流”(chain flow),它從沿海城區(qū)的高層向內(nèi)陸城區(qū)的低層流動。如果內(nèi)陸城區(qū)存在的話,所引起的海風—熱島環(huán)流系統(tǒng)比較強,維持時間也比較長;當這兩個城區(qū)之間的距離加寬時,這樣的特征會更明顯。它是由沿海城區(qū)海風—熱島環(huán)流與內(nèi)陸城區(qū)熱島環(huán)流的合并而產(chǎn)生的;Ohashi and Kida(2002b)使用Ohashi and Kida(2002a)中所用的模式的二維版本研究了沿海城區(qū)海風鋒前部形成的“弱風區(qū)”特征及其形成機制,結(jié)果表明:“弱風區(qū)”是由城區(qū)熱島環(huán)流的發(fā)展所引起的,它的空間尺度極大地依賴于城市的大小(寬度),即依賴于城區(qū)發(fā)展起來的熱島環(huán)流的空間尺度。Ohashi and Kida(2002c)使用Ohashi and Kida(2002a)中所用的模式再次研究了日本大阪—京都平原地區(qū)城區(qū)對白天局地環(huán)流的影響,發(fā)現(xiàn)海風環(huán)流前部的上升運動和海風環(huán)流的垂直尺度大于沒有城區(qū)時的情況。Martilli(2003)利用一個包含詳細城市參數(shù)化的二維中尺度模式研究了城市對海風的影響,發(fā)現(xiàn)城市的存在加速了海風的形成,但減慢了海風鋒的傳播。

Freitas et al.(2007)利用區(qū)域大氣模式RAMS(regional atmospheric modelling system)模擬了巴西圣保羅都市區(qū)超大城市環(huán)境中城市熱島和海風環(huán)流的相互作用,結(jié)果表明:1)城市熱島使得在城市中心形成了強輻合區(qū),因此加速了海風鋒向城市中心的傳播;2)與沒有城市的情況相比,城區(qū)的存在使海風鋒傳播的平均速度增加了約0.32 m/s;3)在海風鋒到達城區(qū)之后,海風鋒在城市中心滯留約2 h,而熱島消失后,海風鋒穿過城市繼續(xù)前進。Dandou et al.(2009)利用中尺度模式MM5(fifth-generation Penn State/NCAR mesoscale model)模擬了希臘大型城市雅典對海風鋒演變的影響,發(fā)現(xiàn)在白天海風鋒的移動被推遲,城市熱島加強了海風,而海風使城市熱島向內(nèi)陸產(chǎn)生了位移。

在我國,李維亮等(2003)利用中尺度模式MM5研究了長江三角洲地區(qū)城市熱島對局地環(huán)流的影響,指出在有城市熱島的情況下,城市熱島環(huán)流與海陸風之間存在相互作用,明顯地反映在局地環(huán)流的改變上。Lo et al.(2007)利用MM5研究了珠江三角洲城市化及與城市化相關(guān)的城市熱島效應(yīng)對局地、區(qū)域大氣環(huán)流的影響,發(fā)現(xiàn)在下午城市化加強了海風環(huán)流,在夜間和清晨,城市熱島環(huán)流效應(yīng)弱化了陸風;珠江三角洲地區(qū)較強的城市熱島增加了城區(qū)和附近海面之間的溫度梯度差異,進而增強了海風環(huán)流;城市化加強了白天的海風環(huán)流,也增加了大氣層最下部2 km內(nèi)的氣溫。陸希和壽紹文(2009)利用MM5對珠江三角洲進行的城市化敏感性試驗發(fā)現(xiàn)城市化對該地區(qū)海陸風的結(jié)構(gòu)和強度有很大的影響,中午增強的海風使得珠江三角洲東岸的低層輻合加強。Lu et al.(2010)使用MM5模擬了冬季珠江三角洲地區(qū)城市化對海陸風環(huán)流的影響,也發(fā)現(xiàn)珠江三角洲地區(qū)城市化對海陸風的分布和環(huán)流特征有顯著的影響。劉樹華等(2009)使用MM5模擬了京津冀地區(qū)不同季節(jié)溫度場和風場等邊界層特征量及其變化特征,結(jié)果表明:在弱天氣系統(tǒng)控制下,京津冀地區(qū)大氣邊界層中可同時存在海陸風、山谷風和城市熱島環(huán)流,同時三者之間還存在明顯的耦合效應(yīng);海陸風環(huán)流極盛時可深入陸地200 km左右,山谷風環(huán)流的最大影響可覆蓋北京區(qū)域內(nèi)的平原地區(qū),而城市熱島環(huán)流則發(fā)生在城市中心幾十公里范圍內(nèi),并對海陸風環(huán)流和山谷風環(huán)流起明顯的削弱作用。

近些年來,中尺度數(shù)值模式WRF(weather research and forecasting model)的不斷發(fā)展和完善,以及城市冠層參數(shù)化(urban canopy parameterization,UCP)的深入研究極大推動了城市熱島和海風環(huán)流相互作用的研究。Lin et al.(2008)利用WRF模式模擬了臺灣北部城市熱島效應(yīng)對邊界層的發(fā)展以及海陸風環(huán)流的影響,探討了人為熱源和土地利用的影響,結(jié)果表明:人為熱源對臺北地區(qū)邊界層的發(fā)展和城市熱島強度起重要的作用,城市熱島效應(yīng)對海陸風環(huán)流有顯著的影響,主要表現(xiàn)在白天加強海風,夜間削弱陸風。文偉俊等(2009)利用WRF模式模擬了廣西北部灣地區(qū)海陸風的結(jié)構(gòu),發(fā)現(xiàn)使用“修正后的城市下墊面試驗”模擬的海陸風強度和分布更接近實測,城市化加強了海風的發(fā)展,但削弱了陸風。Chen et al.(2011a)利用耦合了單層(single-layer)城市冠層模式(urban canopy model,UCM)的WRF模式研究了美國休斯頓地區(qū)局地陸面和城市強迫對海陸風環(huán)流演變的影響,發(fā)現(xiàn)城市的存在有利于地面風的停滯。Chemel and Sokhi(2012)利用耦合了城市冠層模式(UCM)、“建筑物效應(yīng)參數(shù)化”(building effect parameterization,BEP)、“建筑物能量模型”(building energy model,BEM)的WRF模式系統(tǒng)研究了英國北海(North Sea)地區(qū)海風的侵入對倫敦城市熱島的影響,發(fā)現(xiàn)海風減少了城市熱島的強度,并改變了城市熱島的結(jié)構(gòu),使城市熱島中心發(fā)生了位移。Ryu and Baik(2013)利用耦合了韓國首爾國立大學(xué)城市冠層模式的中尺度模式WRF研究了首爾都市區(qū)局地環(huán)流的相互作用,發(fā)現(xiàn)城市風(urban breeze)扮演著抑制海風向內(nèi)陸縱深發(fā)展的角色,使海風停滯在首爾。張亦洲等(2013)利用耦合了Noah陸面過程模式(land surface model,LSM)和單層城市冠層模式的中尺度模式WRF模擬系統(tǒng)(WRF-Noah-UCM)模擬了京津城市下墊面對海風的影響,發(fā)現(xiàn)京津城市對海風的發(fā)展和推進過程有明顯影響,能夠阻礙海風的推進,北京城市下墊面還能在海風到達前增加其強度和向內(nèi)陸傳播速度,并在海風經(jīng)過后延緩其消亡、增加其向內(nèi)陸傳播距離。梁釗明等(2013)利用耦合了城市冠層參數(shù)化(UCP)和建筑物能量模型(BEM)的中尺度模式WRF模擬系統(tǒng)研究了城市(天津)下墊面對渤海灣海風鋒特征的影響,發(fā)現(xiàn)城市下墊面高粗糙度對低層海風風速的明顯削弱造成海風鋒往內(nèi)陸推進距離稍減,高粗糙度的城市下墊面對海風環(huán)流的摩擦力效應(yīng)使得海風得到抬升。

從上述的研究進展來看,城市熱島和海風環(huán)流相互作用的數(shù)值模擬研究已經(jīng)在不少地區(qū)開展,并得到科學(xué)家們較為廣泛的重視。美國科學(xué)家Schultz and Warner(1982)和芬蘭科學(xué)家Savij?rvi(1985)是這個研究領(lǐng)域的開拓者之一。日本科學(xué)家Yoshikado(1992,1994)、Kusaka et al.(2000)、Ohashi and Kida(2002a,2002b,2002c)、印度科學(xué)家Sarkar et al.(1998)、瑞士科學(xué)家Martilli(2003)、巴西科學(xué)家Freitas et al.(2007)以及希臘科學(xué)家Dandou et al.(2009)的工作比較有代表性。從研究內(nèi)容上來看,近30 a來城市化、城市熱島效應(yīng)對海風環(huán)流的影響研究較多,而海風對城市熱島、城市熱島環(huán)流的影響研究則較少,城市熱島和海風環(huán)流相互作用機理的研究也不多見。

2 評述與展望

早期的數(shù)值研究以二維(2D)數(shù)值試驗(Schultz and Warner,1982;Savij?rvi,1985;Yoshikado,1992,1994;Ohashi and Kida,2002b;Martilli,2003)和三維(3D)理想數(shù)值模擬(張雷鳴等,1994;苗曼倩和唐有華,1998;Kusaka et al.,2000;Ohashi and Kida,2002a,2002c)為主,實際個例的三維數(shù)值模擬研究則是進入21世紀才開展起來的。早期研究所使用的數(shù)值模式基本上是理論分析模式,物理過程過于簡單,也不完整,缺乏對實際大氣的模擬能力,所做的模擬是“虛擬的”(理想狀態(tài)下的模擬)。

近一、二十年以來,中尺度數(shù)值模式的發(fā)展和應(yīng)用對城市熱島和邊界層的研究起到了極大的推動作用(Miao et al.,2008;蔣維楣等,2010;壽亦萱和張大林,2012)。21世紀以來,包含全物理過程的三維中尺度數(shù)值模式RAMS、MM5、WRF在城市熱島和海風環(huán)流相互作用的研究中得到了較為廣泛的應(yīng)用。比如,李維亮等(2003)采用中尺度模式MM5,Lo et al.(2007)、劉樹華等(2009)以及Lu et al.(2010)采用耦合了陸面過程模式Noah的中尺度模式MM5模擬系統(tǒng)(MM5-Noah)對實例進行了三維數(shù)值模擬研究,但MM5-Noah耦合中尺度模式只是簡單地在陸面模式Noah中修改了植被類型“城市”的特征參數(shù)值(粗糙度長度、地面反照率、熱容量、土壤熱傳導(dǎo)率、綠色植被覆蓋度和土壤濕度)(Liu et al.,2006;Miao et al.,2007),并不包括真正意義上的城市冠層參數(shù)化(UCP)方案。Freitas et al.(2007)采用區(qū)域大氣模擬系統(tǒng)(RAMS)對實例進行了三維數(shù)值模擬研究,并在大氣模式中耦合了“城市地表能量平衡”(town energy balance,TEB)模式,這使得他們的研究工作較他們之前的同類研究在方法學(xué)上有了很大的改進。然而,在Freitas et al.(2007)的研究中,為了識別城市熱島環(huán)流和海風的相互作用,在數(shù)值試驗中沒有考慮“地形”(topography),這顯然是不真實的。Dandou et al.(2009)利用中尺度模式MM5的“城市修改版”通過對實際個例的三維數(shù)值模擬研究了大型城市雅典對海風鋒演變的影響,他們的研究較以前使用MM5所做的同類研究在方法學(xué)上占有一定優(yōu)勢。

近年來,隨著城市冠層參數(shù)化(UCP)、建筑物效應(yīng)參數(shù)化(BEP)、建筑物能量模型(BEM)研究的不斷深入(Holt and Pullen,2007;Chen et al.2011b;Lee et al.2011;Salamanca et al.,2011;王詠薇等,2013;張艷霞等,2013),以及中尺度數(shù)值模式的不斷發(fā)展和完善,WRF-Noah-UCM、WRF-BEP、WRF-UCP-BEM、WRF-BEP-BEM等模擬系統(tǒng)在城市熱島與海風環(huán)流相互作用“實際個例”三維數(shù)值模擬研究中得到了初步應(yīng)用(Lin et al.,2008;Chen et al.,2011a;Chemel and Sokhi,2012;梁釗明等,2013;張亦洲等,2013;Ryu and Baik,2013)。從方法論的角度來看,使用耦合UCM、BEP、UCP、BEM的中尺度模式模擬系統(tǒng)研究“真實大氣”中城市熱島和海風環(huán)流的相互作用是未來發(fā)展的方向。

值得注意的是:盡管不少研究已表明大尺度背景風不僅對海風環(huán)流有顯著的影響(Crosman and Horel,2010),而且對城市熱島也有顯著的影響(霍飛和陳海山,2011;Wu et al.,2011),但以前的研究很少關(guān)注大尺度背景風對城市熱島和海風環(huán)流相互作用或耦合效應(yīng)的影響,絕大部分研究中的大尺度背景風是靜風或弱的定常風。此外,以前絕大多數(shù)城市熱島和海風環(huán)流相互作用的數(shù)值模擬研究并不考慮云過程,基本上是干模式結(jié)構(gòu),但多云天氣條件下城市熱島和海風環(huán)流相互作用的研究會更有意義。另外,隨著大渦模擬(large eddy simulation,LES)技術(shù)的不斷發(fā)展和完善,科學(xué)家們開始關(guān)注城市熱島、城市熱島環(huán)流的大渦模擬(Zhang et al.,2014)和海風的大渦模擬(Antonelli and Rotunno,2007)研究,但目前城市熱島和海風環(huán)流相互作用的大渦模擬研究還極為少見。

綜上所述,城市熱島和海風環(huán)流相互作用的研究在近30 a來已有了長足的進展,但研究的深度和廣度還遠遠不夠,尚有許多問題值得進一步探討,譬如:

1)中尺度數(shù)值模式中城市冠層參數(shù)化的適應(yīng)性與評估;

2)沿海城市地區(qū)非均勻復(fù)雜下墊面特征(地形、土地利用、人為熱源等)對局地環(huán)流的影響;

3)海風影響下,城市熱島、城市熱島環(huán)流的演變;

4)大尺度背景風對城市熱島和海風環(huán)流相互作用的影響;

5)大氣穩(wěn)定度或天氣條件(晴天和多云)對城市熱島和海風環(huán)流相互作用的影響;

6)城市熱島和海風環(huán)流耦合效應(yīng)對沿海城市地區(qū)局地強對流天氣的觸發(fā)作用;

7)城市熱島環(huán)流和海風環(huán)流在局地對流(云降水)中的具體作用;

8)城市熱島和海風環(huán)流相互作用的大渦模擬;

9)城市熱島和海風環(huán)流相互作用的物理機制、概念模型。

此外,湖陸風是一種形成機制類似于海陸風的中小尺度環(huán)流(Laird et al.,2001;Crosman and Horel,2010),對空氣質(zhì)量(air quality)、天氣和氣候均有重要的影響(Harris and Kotamarthi,2005;楊薇等,2014)。因此,沿湖城市熱島與湖陸風相互作用的研究也已引起科學(xué)家們的關(guān)注(李維亮等,2003;Crosman and Horel,2010;Keeler and Kristovich,2012;楊健博等,2013),其相互作用機理本質(zhì)上與沿海城市熱島—海陸風相互作用機理一樣。然而,相對于城市熱島、城市化對海陸風的影響研究而言,城市熱島效應(yīng)對湖陸風的影響研究還不夠廣泛和深入,需要進一步加強。

另外,值得一提的是:盡管數(shù)值試驗是研究城市熱島—海風環(huán)流(湖風環(huán)流)相互作用最常使用的方法,但觀測試驗也是一種重要的研究手段,具有不可替代性。迄今為止,相對于數(shù)值試驗研究而言,城市熱島—海風環(huán)流(湖風環(huán)流)相互作用的觀測試驗研究(Barbato,1978;Bornstein and Thompson,1981;于恩洪等,1987;Yoshikado and Kondo,1989;Yoshikado,1990;Ohashi and Kida,2001;Gedzelman et al.,2003;Childs and Raman,2005;Keeler and Kristovich,2012)還不夠廣泛,這在很大程度上制約于研究成本與邊界層觀測手段。隨著科研經(jīng)費投入的不斷增長、邊界層氣象觀測手段的日益豐富,城市熱島—海風環(huán)流(湖風環(huán)流)相互作用的觀測試驗研究也應(yīng)是未來發(fā)展的重要方向之一。

致謝:兩位匿名審稿專家、邱新法教授對本文提出了建設(shè)性的修改意見,2013級氣象學(xué)碩士研究生韓芙蓉協(xié)助整理了參考文獻,在此表示衷心的感謝!

陳訓(xùn)來,馮業(yè)榮,范紹佳,等.2008.離岸型背景風和海陸風對珠江三角洲地區(qū)灰霾天氣的影響[J].大氣科學(xué),32(3):530-542.

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

東高紅,劉一瑋,孫蜜娜,等.2013.城市熱島與海風鋒疊加作用對一次局地強降水的影響[J].氣象,39(11):1422-1430.

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

黃利萍,苗峻峰,劉月琨,等.2013.天津地區(qū)夏季海陸風對城市熱島日變化特征影響的觀測分析[J].大氣科學(xué)學(xué)報,36(4):417-425.

霍飛,陳海山.2011.大尺度環(huán)境風場對城市熱島效應(yīng)影響的數(shù)值模擬試驗[J].氣候與環(huán)境研究,16(6):679-689.

蔣維楣,王詠薇,張寧.2009.城市陸面過程與邊界層結(jié)構(gòu)研究[J].地球科學(xué)進展,24(4):411-419.

蔣維楣,苗世光,張寧,等.2010.城市氣象與邊界層數(shù)值模擬研究[J].地球科學(xué)進展,25(5):463-473.

李維亮,劉洪利,周秀驥,等.2003.長江三角洲城市熱島與太湖對局地環(huán)流影響的分析研究[J].中國科學(xué)D輯,33(2):97-104.

梁釗明,高守亭,王東海,等.2013.城市下墊面對渤海灣海風鋒特征影響的一次數(shù)值試驗[J].大氣科學(xué),37(5):1013-1024.

劉樹華,劉振鑫,李炬,等.2009.京津冀地區(qū)大氣局地環(huán)流耦合效應(yīng)的數(shù)值模擬[J].中國科學(xué)D輯,39(1):88-98.

盧煥珍,劉一瑋,劉愛霞,等.2012.海風鋒導(dǎo)致雷暴生成和加強規(guī)律研究[J].氣象,38(9):1078-1086.

陸希,壽紹文.2009.珠三角地區(qū)海陸風的數(shù)值模擬及城市化對其影響研究[J].安徽農(nóng)業(yè)科學(xué),37(1):312-315.

蒙偉光,閆敬華,扈海波.2007.城市化對珠江三角洲強雷暴天氣的可能影響[J].大氣科學(xué),31(2):364-376.

苗曼倩,唐有華.1998.長江三角洲夏季海陸風與熱島環(huán)流的相互作用及城市化的影響[J].高原氣象,17(3):280-289.

壽亦萱,張大林.2012.城市熱島效應(yīng)的研究進展與展望[J].氣象學(xué)報,70(3):338-353.

王樹芬.1990.一次由海風鋒觸發(fā)的強對流天氣分析[J].大氣科學(xué),14(4):504-507.

汪雅,苗峻峰,談?wù)苊?2013.寧波地區(qū)海陸下墊面差異對雷暴過程影響的數(shù)值模擬[J].氣象學(xué)報,71(6):1146-1159.

王詠薇,伍見軍,杜欽,等.2013.不同城市冠層參數(shù)化方案對重慶髙密度建筑物環(huán)境的數(shù)值模擬研究[J].氣象學(xué)報,71(6):1130-1145.

文偉俊,沈桐立,丁治英,等.2009.城市化對廣西夏季海陸風影響的數(shù)值試驗[J].熱帶氣象學(xué)報,25(3):350-356.

許啟慧,苗峻峰,劉月琨,等.2013a.渤海灣西岸海風時空演變特征觀測分析[J].海洋預(yù)報,30(1):9-19.

許啟慧,苗峻峰,劉月琨,等.2013b.渤海灣西岸海陸風特征對城市熱島響應(yīng)的觀測分析[J].氣象科學(xué),33(4):408-417.

徐蓉,苗峻峰,談?wù)苊?2013.南京地區(qū)城市下墊面特征對雷暴過程影響的數(shù)值模擬[J].大氣科學(xué),37(6):1235-1246.

徐祥德,丁國安,卞林根.2006.北京城市大氣環(huán)境污染機理與調(diào)控原理[J].應(yīng)用氣象學(xué)報,17(6):815-828.

楊健博,劉紅年,費松,等.2013.太湖湖陸風背景下的蘇州城市化對城市熱島特征的影響[J].氣象科學(xué),33(5):473-484.

楊薇,苗峻峰,談?wù)苊?2014.太湖地區(qū)湖陸風對雷暴過程影響的數(shù)值模擬[J].應(yīng)用氣象學(xué)報,25(1):59-70.

易笑園,張義軍,沈永海,等.2012.一次海風鋒觸發(fā)的多單體雹暴及合并過程的觀測分析[J].氣象學(xué)報,70(5):974-985.

尹東屏,吳海英,張備,等.2010.一次海風鋒觸發(fā)的強對流天氣分析[J].高原氣象,29(5):1261-1269.

于恩洪,陳彬,白玉榮.1987.渤海灣西部海陸風的空間結(jié)構(gòu)[J].氣象學(xué)報,45(3):379-381.

張雷鳴,苗曼倩,洪鐘祥,等.1994.城市發(fā)展對夜間海陸風環(huán)流影響的預(yù)測模擬[J].大氣科學(xué),18(3):366-372.

張立鳳,張銘,林宏源.1999.珠江口地區(qū)海陸風系的研究[J].大氣科學(xué),23(5):581-589.

張艷霞,蒙偉光,戴光豐,等.2013.WRF耦合城市冠層模式對珠三角城市群天氣模擬影響的評估[J].熱帶氣象學(xué)報,29(6):935-946.

張亦洲,苗世光,戴永久,等.2013.北京夏季晴天邊界層特征及城市下墊面對海風影響的數(shù)值模擬[J].地球物理學(xué)報,56(8):2558-2573.

周明煜,曲紹厚,李玉英,等.1980.北京地區(qū)熱島和熱島環(huán)流特征[J].環(huán)境科學(xué),1(5):12-18.

Antonelli M,Rotunno R.2007.Large-eddy simulation of the onset of the sea breeze[J].J Atmos Sci,64:4445-4457.

Atkinson B W.2003.Numerical modelling of urban heat-island intensity[J].Bound-Layer Meteor,109:285-310.

Barbato J P.1978.Areal parameters of the sea breeze and its vertical structure in the Boston basin[J].Bull Amer Meteor Soc,59:1420-1431.

Bornstein R D,Thompson W T.1981.Effects of frictionally retarded sea breeze and synoptic frontal passages on sulfur dioxide concentrations in New York City[J].J Appl Meteor,20:843-858.

Cenedese A,Monti P.2003.Interaction between an inland urban heat island and a sea-breeze flow:A laboratory study[J].J Appl Meteor,42:1569-1583.

Chemel C,Sokhi R S.2012.Response of London’s urban heat island to a marine air intrusion in an easterly wind regime[J].Bound-Layer Meteor,144:65-81.

Chen F,Miao S,Tewari M,et al.2011a.A numerical study of interactions between surface forcing and sea breeze circulations and their effects on stagnation in the greater Houston area[J].J Geophys Res,116,D12105.doi:10.1029/2010JD015533.

Chen F,Kusaka H,Bornstein R,et al.2011b.The integrated WRF/urban modelling system:Development,evaluation,and applications to urban environmental problems[J].Int J Climatol,31:273-288.

Chen T C,Wang S Y,Yen M C.2007.Enhancement of afternoon thunderstorm activity by urbanization in a valley:Taipei[J].J Appl Meteor Climatol,46:1324-1340.

Childs P P,Raman S.2005.Observations and numerical simulations of urban heat island and sea breeze circulations over New York City[J].Pure Appl Geophy,162:1955-1980.

Crosman E T,Horel J D.2010.Sea and lake breezes:A review of numerical studies[J].Bound-Layer Meteor,137:1-29.

Dandou A,Tombrou M,Soulakellis N.2009.The influence of the City of Athens on the evolution of the sea-breeze front[J].Bound-Layer Meteor,131:35-51.

Eliasson I,Holmer B.1990.Urban heat island circulation in G?teborg,Sweden[J].Theor Appl Climatol,42:187-196.

Ferretti R,Mastrantonio G,Argentini S,et al.2003.A model-aided investigation of winter thermally driven circulation on the Italian Tyrrhenian coast:A case study[J].J Geophys Res,108(D24),4777.doi:10.1029/2003JD003424.

Fovell R G.2005.Convective initiation ahead of the sea-breeze front[J].Mon Wea Rev,133:264-278.

Freitas E D,Rozoff C M,Cotton W R,et al.2007.Interactions of an urban heat island and sea-breeze circulations during winter over the metropolitan area of S?o Paulo,Brazil[J].Bound-Layer Meteor,122:43-65.

Gedzelman S D,Austin S,Cermak R,et al.2003.Mesoscale aspects of the urban heat island around New York City[J].Theor Appl Climatol,75:29-42.

Grimmond C S B,Blackett M,Best M J,et al.2010.The international urban energy balance models comparison project:First results from phase 1[J].J Appl Meteor Climatol,49:1268-1292.

Guo X,Fu D,Wang J.2006.Mesoscale convective precipitation system modified by urbanization in Beijing City[J].Atmos Res,82:112-126.

Han J Y,Baik J J.2008.A theoretical and numerical study of urban heat island-induced circulation and convection[J].J Atmos Sci,65:1859-1877.

Harris L,Kotamarthi V R.2005.The characteristics of the Chicago lake breeze and its effects on trace particle transport:Results from an episodic event simulation[J].J Appl Meteor,44:1637-1654.

Hidalgo J,Pigeon G,Masson V.2008.Urban-breeze circulation during the CAPITOUL experiment:Observational data analysis approach[J].Meteor Atmos Phys,102:223-241.

Hidalgo J,Masson V,Gimeno L.2010.Scaling the daytime urban heat island and urban-breeze circulation[J].J Appl Meteor Climatol,49:889-901.

Holt T,Pullen J.2007.Urban canopy modeling of the New York City metropolitan area:A comparison and validation of single-and multilayer parameterizations[J].Mon Wea Rev,135:1906-1930.

Kanda M.2007.Progress in urban meteorology:A review[J].J Meteor Soc Japan,85B:363-383.

Keeler J M,Kristovich D A R.2012.Observations of urban heat island influence on lake-breeze frontal movement[J].J Appl Meteor Climatol,51:702-710.

Kusaka H,Kimura F,Hirakuchi H,et al.2000.The effects of land-use alteration on the sea breeze and daytime heat island in the Tokyo metropolitan area[J].J Meteor Soc Japan,78:405-420.

Laird N F,Kristovich D A R,Liang X Z,et al.2001.Lake Michigan lake breezes:Climatology,local forcing,and synoptic environment[J].J Appl Meteor,40:409-424.

Lee S H,Kim S W,Angevine W M,et al.2011.Evaluation of urban surface parameterizations in the WRF model using measurements during the Texas Air Quality Study 2006 field campaign[J].Atmos Chem Phys,11:2127-2143.

Lemonsu A,Pigeon G,Masson V,et al.2006.Sea-town interactions over Marseille:3D urban boundary layer and thermodynamic fields near the surface[J].Theor Appl Climatol,84:171-178.

Lin C Y,Chen F,Huang J C,et al.2008.Urban heat island effect and its impact on boundary layer development and land-sea circulation over northern Taiwan[J].Atmos Environ,42:5635-5649.

Lin C Y,Chen W C,Chang P L,et al.2011.Impact of the urban heat island effect on precipitation over a complex geographic environment in Northern Taiwan[J].J Appl Meteor Climatol,50:339-353.

Lin W S,Wang A Y,Wu C S,et al.2001.A case modeling of sea-land breeze in Macao and its neighborhood[J].Adv Atmos Sci,18:1231-1240.

Liu Y B,Chen F,Warner T,et al.2006.Verification of a mesoscale data-assimilation and forecasting system for the Oklahoma City area during the Joint Urban 2003 field project[J].J Appl Meteor Climatol,45:912-929.

Lo J C F,Lau A K H,Chen F,et al.2007.Urban modification in a mesoscale model and the effects on the local circulation in the Pearl River Delta Region[J].J Appl Meteor Climatol,46:457-476.

Lu X,Chow K C,Yao T,et al.2010.Effects of urbanization on the land sea breeze circulation over the Pearl River Delta region in winter[J].Int J Climatol,30:1089-1104.

Martilli A.2003.A two-dimensional numerical study of the impact of a city on atmospheric circulation and pollutant dispersion in a coastal environment[J].Bound-Layer Meteor,108:91-119.

Meir T,Orton P M,Pullen J,et al.2013.Forecasting the New York City urban heat island and sea breeze during extreme heat events[J].Wea Forecasting,28:1460-1477.

Miao J F,Kroon L J M,Vil-Guerau de Arellano J,et al.2003.Impacts of topography and land degradation on the sea breeze over eastern Spain[J].Meteor Atmos Phys,84:157-170.

Miao J F,Chen D,Borne K.2007.Evaluation and comparison of Noah and Pleim-Xiu land surface models in MM5 using G?TE2001 data:Spatial and temporal variations in near-surface air temperature[J].J Appl Meteor Climatol,46:1587-1605.

Miao J F,Chen D,Wyser K,et al.2008.Evaluation of MM5 mesoscale model at local scale for air quality applications over the Swedish west coast:Influence of PBL and LSM parameterizations[J].Meteor Atmos Phys,99:77-103.

Miao J F,Wyser K,Chen D,et al.2009.Impacts of boundary layer turbulence and land surface process parameterizations on simulated sea breeze characteristics[J].Ann Geophys,27:2303-2320.

Miao S G,Chen F,Li Q C,et al.2011.Impacts of urban processes and urbanization on summer precipitation:A case study of heavy rainfall in Beijing on 1 August 2006[J].J Appl Meteor Climatol,50:806-825.

Miller S T K,Keim B D,Talbot R W,et al.2003.Sea breeze:Structure,forecasting,and impacts[J].Rev Geophys,41:1011.doi:10.1029/2003RG000124.

Ohashi Y,Kida H.2001.Observational results of the sea breeze with a weak wind region over the northern Osaka urban area[J].J Meteor Soc Japan,79:949-955.

Ohashi Y,Kida H.2002a.Local circulations developed in the vicinity of both coastal and inland urban areas:A numerical study with a mesoscale atmospheric model[J].J Appl Meteor,41:30-45.

Ohashi Y,Kida H.2002b.Numerical experiments on the weak-wind region formed ahead of the sea-breeze front[J].J Meteor Soc Japan,80:519-527.

Ohashi Y,Kida H.2002c.Effects of mountain and urban areas on daytime local-circulations in the Osaka and Kyoto regions[J].J Meteor Soc Japan,80:539-560.

Oke T R.2006.Towards better scientific communication in urban climate[J].Theor Appl Climatol,84:179-190.

Ryu Y H,Baik J J.2013.Daytime Local Circulations and Their Interactions in the Seoul Metropolitan Area[J].J Appl Meteor Climatol,52:784-801.

Salamanca F,Martilli A,Tewari M,et al.2011.A study of the urban boundary layer using different urban parameterizations and high-resolution urban canopy parameters with WRF[J].J Appl Meteor Climatol,50:1107-1128.

Sarkar A,Saraswat R,Chandrasekar A.1998.Numerical study of the effects of urban heat island on the characteristic features of the sea breeze circulation[J].Proc Indian Acad Sci-Earth Planet Sci,107:127-137.

Savij?rvi H.1985.The sea breeze and urban heat island circulation in a numerical model[J].Geophysica,21:115-126.

Schultz P,Warner T T.1982.Characteristics of summertime circulations and pollutant ventilation in the Los Angeles basin[J].J Appl Meteor,21:672-682.

Shepherd J M.2005.A review of current investigations of urban-induced rainfall and recommendations for the future[J].Earth Interact,9(12):1-27.

Simpson M,Raman S,Suresh R,et al.2008.Urban effects of Chennai on sea breeze induced convection and precipitation[J].J Earth Syst Sci,117:897-909.

Vukovich F M,King W J,Dunn J W,et al.1979.Observations and simulations of the diurnal variation of the urban heat island circulation and associated variations of the ozone distribution:a case study[J].J Appl Meteor,18:836-854.

Wu J B,Chow K C,Fung J C H,et al.2011.Urban heat island effects of the Pearl River Delta city clusters-their interactions and seasonal variation[J].Theor Appl Climatol,103:489-499.

Yoshikado H.1990.Vertical structure of the sea breeze penetrating through a large urban complex[J].J Appl Meteor,29:878-891.

Yoshikado H.1992.Numerical study of the daytime urban effect and its interaction with the sea breeze[J].J Appl Meteor,31:1146-1164.

Yoshikado H.1994.Interaction of the sea breeze with urban heat islands of different sizes and locations[J].J Meteor Soc Japan,72:139-143.

Yoshikado H,Kondo H.1989.Inland penetration of the sea breeze over the suburban area of Tokyo[J].Bound-Layer Meteor,48:389-407.

Zhang C L,Chen F,Miao S G,et al.2009.Impacts of urban expansion and future green planting on summer precipitation in the Beijing metropolitan area[J].J Geophys Res,114,D02116.doi:10.1029/2008JD010328.

Zhang N,Wang X Y,Peng Z.2014.Large-eddy simulation of mesoscale circulations forced by inhomogeneous urban heat island[J].Bound-Layer Meteor,151:179-194.

(責任編輯：張福穎)

Anoverviewofnumericalstudiesofinteractionofurbanheatislandandseabreezecirculations

MIAO Jun-feng

(School of Atmospheric Sciences,NUIST,Nanjing 210044,China)

Urban heat island circulation(UHIC) and sea breeze circulation(SBC) usually coexist in coastal urban areas;they play a crucial role in transport and dispersion of air pollutants,and have a significant effect on convection initiation.With the development of urbanization,environmental problems such as air pollution become more and more serious,and disaster induced by severe convective weather becomes more and more serious too.Therefore,the study of UHIC and SBC interaction receives increasing and extensive attention.In this paper,the history of numerical studies of UHIC and SBC interaction in last 30 years is overviewed;The current research status and some problems are analyzed;Some important results are outlined;Some issues and topics are addressed for further study.

heat island circulation;sea and land breezes;lake and land breezes;local circulation;local convection;urban effect;air pollution dynamics

2014-05-03;改回日期2014-06-02

公益性行業(yè)(氣象)科研專項經(jīng)費項目(GYHY201306009);中荷科學(xué)合作“中國交流項目”(530-4CDP09);江蘇省高?！扒嗨{工程”科技創(chuàng)新團隊項目(東亞季風與區(qū)域氣候變化)

苗峻峰,博士,教授,博士生導(dǎo)師,研究方向為中尺度數(shù)值模擬,miaoj@nuist.edu.cn.

10.13878/j.cnki.dqkxxb.20140503001.

1674-7097(2014)04-0521-08

P435.2

A

10.13878/j.cnki.dqkxxb.20140503001

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

Miao Jun-feng.2014.An overview of numerical studies of interaction of urban heat island and sea breeze circulations[J].Trans Atmos Sci,37(4):521-528.(in Chinese)