• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    蘇州—無(wú)錫—常州城市帶熱島效應(yīng)個(gè)例研究

    2014-09-09 07:30:11康漢青朱彬朱彤高晉徽孫佳麗蘇繼峰
    大氣科學(xué)學(xué)報(bào) 2014年4期

    康漢青,朱彬,朱彤,高晉徽,孫佳麗,蘇繼峰

    (1.南京信息工程大學(xué) 氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,江蘇 南京 210044;2.南京信息工程大學(xué) 中國(guó)氣象局氣溶膠與云降水重點(diǎn)開(kāi)放實(shí)驗(yàn)室,江蘇 南京 210044;3.科羅拉多州立大學(xué),美國(guó) 科羅拉多州 80523-1375;4.NOAA/NESDIS/STAR/JCSDA,美國(guó) 馬里蘭州 20746;5.江蘇省氣候中心,江蘇 南京 210009;6.中國(guó)人民解放軍94857部隊(duì)61分隊(duì),安徽 蕪湖 241000)

    蘇州—無(wú)錫—常州城市帶熱島效應(yīng)個(gè)例研究

    康漢青1,2,朱彬1,2,朱彤3,4,高晉徽1,2,孫佳麗5,蘇繼峰6

    (1.南京信息工程大學(xué) 氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心,江蘇 南京 210044;2.南京信息工程大學(xué) 中國(guó)氣象局氣溶膠與云降水重點(diǎn)開(kāi)放實(shí)驗(yàn)室,江蘇 南京 210044;3.科羅拉多州立大學(xué),美國(guó) 科羅拉多州 80523-1375;4.NOAA/NESDIS/STAR/JCSDA,美國(guó) 馬里蘭州 20746;5.江蘇省氣候中心,江蘇 南京 210009;6.中國(guó)人民解放軍94857部隊(duì)61分隊(duì),安徽 蕪湖 241000)

    應(yīng)用WRF(weather research and forecasting)及其耦合的多層城市冠層模式BEP(building energy parameterization),對(duì)2013年8月13日長(zhǎng)江三角洲地區(qū)一次高溫天氣過(guò)程進(jìn)行了模擬。此次過(guò)程盛行東南風(fēng),風(fēng)向與蘇州—無(wú)錫—常州城市帶走向一致。模擬結(jié)果表明:蘇州—無(wú)錫—常州城市帶熱島效應(yīng)明顯,熱島強(qiáng)度向下游城市逐漸增加;在東南風(fēng)作用下,三座城市的熱島連成一片,形成了一個(gè)更強(qiáng)大的熱島環(huán)流。夜晚,邊界層逐漸趨于穩(wěn)定,熱島環(huán)流減弱,有利于熱島溫度向下游地區(qū)輸送。熱島效應(yīng)導(dǎo)致城市邊界層高度明顯上升。白天太湖產(chǎn)生強(qiáng)盛的湖風(fēng)對(duì)其周邊城市影響顯著,來(lái)自太湖的冷氣團(tuán)導(dǎo)致無(wú)錫和常州邊界層內(nèi)熱島強(qiáng)度明顯下降,抑制城市熱島向上發(fā)展,削弱了無(wú)錫與常州兩城市熱島間的聯(lián)系。白天太湖使得無(wú)錫和常州邊界層高度明顯下降。

    城市熱島;數(shù)值模擬;太湖;大氣邊界層

    0 引言

    20世紀(jì)以來(lái),隨著經(jīng)濟(jì)發(fā)展和城市人口迅速上升,城市化進(jìn)入了高速發(fā)展階段。城市化過(guò)程不僅能通過(guò)建筑物和公路等設(shè)施改變地表的幾何結(jié)構(gòu),而且能通過(guò)建筑材料的不同熱力特性以及人類活動(dòng)改變地表的熱力結(jié)構(gòu)。這種改變會(huì)使地表的反照率減小,熱容量增大,吸收存儲(chǔ)更多的熱量,進(jìn)而導(dǎo)致城市地區(qū)增溫顯著,氣溫明顯高于周圍郊區(qū),這種現(xiàn)象就是城市熱島效應(yīng)(Oke,1982;Allwine et al.,2002;Rotach et al.,2005)。

    雖然城市熱島對(duì)全球氣溫變化的影響并不明顯(Houghton et al.,2001),但它會(huì)影響用于氣候變化評(píng)估的局地氣溫(Peterson,2003;Van Weverberg et al.,2008)。在夏季,城市熱島引起的升溫會(huì)影響居民的生活環(huán)境和生活質(zhì)量。有些城市熱島效應(yīng)是正面的,如延長(zhǎng)植物生長(zhǎng)季節(jié)等,但是大多數(shù)城市熱島效應(yīng)是負(fù)面的,包括加劇能源消耗、增加污染物和溫室氣體的排放、危害人類健康、影響舒適度、造成水質(zhì)惡化等。

    數(shù)值模式由于能夠提供城市熱島的時(shí)空分布特征,一直是研究城市熱島最重要的工具之一。Masson(2000)、Kusaka et al.(2001)及Kusaka and Kimura(2004)提出并建立了單層城市冠層模式。單層城市冠層模式考慮了城市的幾何特征,考慮了建筑物對(duì)長(zhǎng)短波輻射收支的影響,能夠較好地反映出城市熱力學(xué)和動(dòng)力學(xué)過(guò)程(Rozoff et al.,2003;Chin et al.,2005;Niyogi et al.,2006;壽亦萱和張大林,2012)。多層城市冠層模式由Martilli et al.(2002)和Dupont et al.(2004)提出,相比于單層城市冠層模式,多層城市冠層模式將城市冠層劃分為若干層,分層次計(jì)算墻壁、道路、屋頂?shù)哪芰渴罩?它考慮了更多的因素,對(duì)城市特征的描述更加準(zhǔn)確。

    長(zhǎng)江三角洲地區(qū)是我國(guó)經(jīng)濟(jì)發(fā)展最迅速的地區(qū)之一,同時(shí)也是城市化進(jìn)程最快的地區(qū)之一。城市化進(jìn)程導(dǎo)致該地區(qū)城市大氣溫度呈逐年上升的趨勢(shì),1999—2005年,城市平均氣溫每年上升0.047 ℃,最高氣溫每年上升0.162 ℃(Du et al.,2007;Zhang et al.,2010)。雖然針對(duì)長(zhǎng)江三角地區(qū)城市熱島效應(yīng)的研究有很多(朱家其等,2006;邱新法等,2008;談建國(guó)等,2008;江志紅和葉麗梅,2010),但是多以單個(gè)城市為主要研究對(duì)象,而城市熱島之間以及城市熱島與海陸風(fēng)、湖陸風(fēng)之間相互影響的研究相對(duì)較少。李維亮等(2003)利用MM5V2研究了長(zhǎng)江三角洲城市熱島與太湖對(duì)局地環(huán)流的影響以及海風(fēng)與湖風(fēng)的相互作用。然而該研究沒(méi)有采用城市冠層,模式對(duì)熱島強(qiáng)度和熱島環(huán)流的表征可能存在一定誤差,且該研究討論的是1995年個(gè)例,經(jīng)過(guò)十幾年的發(fā)展,長(zhǎng)江三角洲城市熱島面積和熱島強(qiáng)度均發(fā)生了顯著變化。目前僅N. Zhang et al.(2011)討論過(guò)蘇州的城市化進(jìn)程對(duì)太湖湖陸風(fēng)的作用。蘇錫常地區(qū)是長(zhǎng)江三角洲內(nèi)經(jīng)濟(jì)發(fā)展迅速、城市化進(jìn)程較快、城市集中且鄰近太湖的區(qū)域,其城市之間熱島效應(yīng)的相互影響以及太湖對(duì)城市熱島環(huán)流的作用值得研究。本文利用WRF3.3.1耦合一個(gè)多層城市冠層模式BEP(building energy parameterization),對(duì)2013年8月13日長(zhǎng)江三角洲地區(qū)一次持續(xù)高溫天氣過(guò)程進(jìn)行模擬,以討論蘇錫常的城市熱島效應(yīng)及其受太湖的影響。

    1 數(shù)值模式和實(shí)驗(yàn)方案設(shè)計(jì)

    本研究采用的數(shù)值模式為中尺度天氣預(yù)報(bào)模式WRF3.3.1,并耦合一個(gè)多層城市冠層模式BEP。模擬區(qū)域如圖1a所示,區(qū)域中心緯度為32.5°N、經(jīng)度為118.0°E。模式采用雙向反饋四重嵌套,水平網(wǎng)格數(shù)分別為180×150、180×150、180×180、240×270,水平網(wǎng)格距分別為13.5、4.5、1.5和0.5 km。模式最外層覆蓋了大半個(gè)中國(guó)以及部分日本島和太平洋地區(qū),最內(nèi)層涵蓋了蘇州、無(wú)錫、常州和太湖地區(qū)。模式垂直方向分為30個(gè)σ層,其中約20層位于2 km高度以下。模式頂氣壓為50 hPa。模式第一層積分步長(zhǎng)為27 s,最內(nèi)層為1 s。模擬起始時(shí)間為2013年8月12日00時(shí)(北京時(shí)間,下同),結(jié)束時(shí)間為2013年8月14日00時(shí)。

    圖1 WRF模擬區(qū)域 a.水平分辨率分別為13.5、4.5、1.5和0.5 km的四重嵌套;b.模擬區(qū)域最內(nèi)層下墊面類型(圖1b中AB線為圖5和圖8垂直剖面所在位置)Fig.1 The WRF model domain a.the four nested meshes with horizontal resolution of 13.5,4.5,1.5 and 0.5 km;b.the land-use categories over the innermost domain(The line AB denotes the location of vertical cross-section in Fig.5 and Fig.8)

    需要注意的是,近十幾年來(lái)長(zhǎng)江三角洲地區(qū)經(jīng)濟(jì)發(fā)展迅速,土地覆蓋變化顯著,尤以城市擴(kuò)展為主。而WRF默認(rèn)使用的USGS數(shù)據(jù)為20世紀(jì)90年代的產(chǎn)品,已經(jīng)不能準(zhǔn)確地表示當(dāng)前的土地類型,尤其是城市區(qū)域。WRF提供的2005年30 s水平分辨率MODIS土地覆蓋產(chǎn)品能基本代表當(dāng)前城市化現(xiàn)狀,因此采用2005 MODIS土地覆蓋資料作為本研究下墊面類型的輸入。

    無(wú)云或少云天氣,風(fēng)速較小,有利于城市熱島的發(fā)生和發(fā)展。雖然在各個(gè)季節(jié)城市熱島效應(yīng)均會(huì)發(fā)生,但是只有夏季的極端高溫天氣過(guò)程會(huì)產(chǎn)生嚴(yán)重危害,且夏季盛行東南季風(fēng),風(fēng)向與蘇錫常城市帶走向一致,有利于討論城市間熱島高溫的輸送。因此,綜合考慮上述因素,本文選取2013年8月13日一次高溫過(guò)程進(jìn)行模擬試驗(yàn)。由MICAPS天氣圖(圖略)可以發(fā)現(xiàn),長(zhǎng)江三角洲地區(qū)位于入發(fā)副熱帶高壓后部,在高壓系統(tǒng)影響下,蘇錫常地區(qū)天氣晴朗,盛行東南風(fēng),有利于城市熱島的形成及其沿蘇州—無(wú)錫—常州城市帶的輸送。為了評(píng)估城市下墊面和太湖對(duì)城市熱島效應(yīng)的影響,本研究的數(shù)值試驗(yàn)設(shè)計(jì)如下:

    1)參照試驗(yàn)(CTL),采用2005 MODIS土地覆蓋資料(圖1b);

    2)無(wú)城市敏感性試驗(yàn)(NUB),將蘇州、無(wú)錫、常州的城市下墊面替換為周圍覆蓋面積最廣的農(nóng)田,其他保持不變;

    3)無(wú)太湖敏感性試驗(yàn)(NTL),將太湖替換為農(nóng)田,其他保持不變。

    2 模擬結(jié)果評(píng)估

    通過(guò)比較模擬與觀測(cè)結(jié)果來(lái)驗(yàn)證耦合的WRF-Noah-BEP模式對(duì)城市氣溫和風(fēng)場(chǎng)的模擬效果。圖2為2013年8月13日蘇州、無(wú)錫、常州三個(gè)城市氣溫、風(fēng)速和風(fēng)向模擬值與觀測(cè)值的比較。模擬的氣溫日變化與實(shí)況較一致,但是夜間模擬結(jié)果偏低,可能與模式?jīng)]有考慮人為熱釋放有關(guān)。研究表明,人為熱可導(dǎo)致近地面氣溫上升1~3 ℃(Ichinose et al.,1999;Ohashi et al.,2007;Ryu et al.,2013)。雖然白天人為熱釋放量大于夜晚,但人為熱對(duì)夜晚的增溫效應(yīng)顯著高于白天(Fan and Sailor,2005;Chen et al.,2009)。模式對(duì)早晚風(fēng)速的模擬值與觀測(cè)值較一致,中午模擬值略微偏高。模擬的風(fēng)向與觀測(cè)值吻合得很好,整個(gè)試驗(yàn)期間盛行東南風(fēng),無(wú)錫和常州中午受太湖湖風(fēng)的影響,風(fēng)向向南風(fēng)轉(zhuǎn)變。

    圖2 8月13日蘇州(SZ;a,b,c)、無(wú)錫(WX;d,e,f)和常州(CZ;g,h,i)近地面氣溫(a,d,g;單位:℃)、風(fēng)速(b,e,h;單位:m·s-1)和風(fēng)向(c,f,i;單位:(°))模擬結(jié)果(實(shí)線)與觀測(cè)值(+)的對(duì)比Fig.2 Comparison of simulated(solid lines) near surface (a,d,g)temperature(units:℃),(b,e,h)wind speed(units:m·s-1) and (c,f,i)wind direction(units:(°)) with the observed results(+) at (a,b,c)Suzhou(SZ),(d,e,f)Wuxi(WX) and (g,h,i)Changzhou(CZ) stations on 13 Auguest

    為了定量描述模型結(jié)果的準(zhǔn)確性,本研究利用均方根誤差(RMSE,ERMS)、平均標(biāo)準(zhǔn)化偏差(MNB,BMN)和平均標(biāo)準(zhǔn)誤差(MNGE,EMNG)等統(tǒng)計(jì)方法對(duì)模擬結(jié)果做進(jìn)一步的量化驗(yàn)證。計(jì)算方法如下:

    (1)

    (2)

    (3)

    其中:pi表示模式預(yù)測(cè)值;oi表示觀測(cè)值;N為樣本數(shù)。

    氣溫、風(fēng)速和風(fēng)向的ERMS分別為1.8 ℃、1.1 m·s-1和17°。三者的BMN分別為-3%、14%和5%,滿足美國(guó)國(guó)家環(huán)保局(EPA)提出的±15%的標(biāo)準(zhǔn)。Tesche et al.(2001)提出了一系列評(píng)價(jià)模式對(duì)氣象場(chǎng)模擬效果的指標(biāo):氣溫的EMNG≤2.0 ℃、風(fēng)速的EMNG≤2.0 m·s-1、風(fēng)向的EMNG≤30°,如果滿足該條件則說(shuō)明模式對(duì)氣象場(chǎng)的模擬較為精確。本次試驗(yàn)氣溫、風(fēng)速和風(fēng)向的EMNG分別為1.5 ℃、0.9 m·s-1和15°??傮w來(lái)說(shuō),模式對(duì)于氣象場(chǎng)的模擬效果較好,能夠真實(shí)反映實(shí)況氣象場(chǎng)的時(shí)空分布。

    3 結(jié)果與討論

    當(dāng)下墊面土地利用類型發(fā)生變化時(shí),其對(duì)短波、長(zhǎng)波輻射的吸收和反射作用都會(huì)隨之改變,從而引起地表感熱通量和潛熱通量的變化,進(jìn)而影響大氣邊界層結(jié)構(gòu)。本文通過(guò)控制試驗(yàn)與敏感性試驗(yàn)之間的對(duì)比,分析城市化水平和太湖的湖風(fēng)效應(yīng)對(duì)城市熱島的影響。為了研究城市化和湖陸風(fēng)對(duì)城市邊界層結(jié)構(gòu)的影響以及不同城市熱島之間的相互作用,本文沿圖1b中直線AB做垂直剖面,剖面貫穿常州(CZ)、無(wú)錫(WX)和蘇州(SZ)三座城市。

    3.1 城市化對(duì)城市熱島效應(yīng)的影響

    圖4 CTL試驗(yàn)與NUB試驗(yàn)白天(a)和夜晚(b)近地面溫度(陰影:單位:℃)和風(fēng)場(chǎng)(箭矢;單位:m·s-1)的差值Fig.4 Differences in averaged near surface temperature(shadings;units:℃) and wind field(arrows;units:m·s-1) between CTL and NUB experiments during (a)daytime and (b)nighttime

    城市由于其熱力特性和幾何特性,能吸收和儲(chǔ)存更多的熱量,加上城市人類活動(dòng)頻繁,釋放出大量的人為熱,造成了城市升溫明顯高于郊區(qū)。圖3給出了CTL試驗(yàn)8月13日白天平均(12:00—17:00)的近地面溫度和風(fēng)場(chǎng)。市區(qū)平均氣溫比郊區(qū)高約2 ℃,區(qū)域盛行東南風(fēng)。風(fēng)場(chǎng)遭遇城市粗糙下墊面阻擋,風(fēng)速下降明顯,且在城市周圍形成繞流。白天,太湖與周圍陸地之間巨大的溫差導(dǎo)致湖風(fēng)效應(yīng)明顯,一定程度上改變了無(wú)錫和常州的風(fēng)場(chǎng)結(jié)構(gòu),進(jìn)而對(duì)這兩座城市的熱島強(qiáng)度和熱島環(huán)流產(chǎn)生影響。圖4a、4b分別為白天和夜晚(18:00—23:00)CTL試驗(yàn)與NUB試驗(yàn)近地面溫度和風(fēng)場(chǎng)的差值,它反映了由城市化引起的溫度和風(fēng)場(chǎng)變化,即城市熱島強(qiáng)度和熱島環(huán)流??梢?jiàn),常州熱島強(qiáng)度最強(qiáng),達(dá)3 ℃,無(wú)錫次之,蘇州最弱,不足2 ℃。受城市熱島的影響,三座城市下風(fēng)向均出現(xiàn)不同程度增溫,增溫幅度與上游城市的熱島強(qiáng)度有關(guān),城市熱島強(qiáng)度越強(qiáng),下游增溫越顯著。D. L. Zhang et al.(2009,2011)提出,上游地區(qū)城市化過(guò)程會(huì)使下游城市的熱島強(qiáng)度加劇。由圖4a可見(jiàn),白天CTL-NUB的風(fēng)場(chǎng)差值呈現(xiàn)向市區(qū)輻合趨勢(shì),與城市熱島環(huán)流特征一致。市區(qū)的風(fēng)場(chǎng)差值呈現(xiàn)為西北風(fēng),說(shuō)明剔除城市后盛行的東南風(fēng)得到了顯著加強(qiáng)。風(fēng)場(chǎng)的輻合以及粗糙的城市地表會(huì)嚴(yán)重削弱城市地區(qū)的風(fēng)速,進(jìn)而抑制城市熱島的高溫向下游地區(qū)的輸送。夜間,由于大氣層結(jié)趨于穩(wěn)定,熱島環(huán)流減弱,風(fēng)場(chǎng)向市區(qū)的輻合被削弱。

    圖3 CTL試驗(yàn)白天平均地表氣溫(陰影:單位:℃)和10 m處風(fēng)場(chǎng)(箭矢;單位:m·s-1)Fig.3 Daytime averaged near surface temperature(shadings;units:℃) and 10 m wind fields(arrows;units:m·s-1) from CTL experiment

    圖5 AB剖面上CTL試驗(yàn)與NUB試驗(yàn)白天(a)和夜晚(b)平均氣溫(彩色陰影;單位:℃)及風(fēng)場(chǎng)(箭矢;單位:m·s-1)的差值(CZ、WX和SZ分別代表常州、無(wú)錫和蘇州)Fig.5 Vertical cross-sections of differences in averaged temperature(color shadings;units:℃) and in-plane flow vectors(arrows;units:m·s-1) between CTL and NUB experiments alone line AB during (a)daytime and (b)nighttime(CZ,WX and SZ represent Changzhou,Wuxi and Suzhou,respectively)

    為了更直觀地反映出城市熱島的垂直結(jié)構(gòu)以及城市之間熱島效應(yīng)的相互影響,本文沿蘇州—無(wú)錫—常州這一城市帶(直線AB)做垂直剖面。圖5a、5b分別為白天和夜晚剖面AB上CTL試驗(yàn)與NUB試驗(yàn)氣溫以及風(fēng)場(chǎng)的差值,它反映出城市熱島強(qiáng)度和熱島環(huán)流的垂直結(jié)構(gòu)??梢?jiàn),白天常州熱島強(qiáng)度最強(qiáng),無(wú)錫次之,蘇州最弱,但是蘇州城市熱島的垂直發(fā)展高度明顯高于無(wú)錫和常州。白天上游城市的熱島效應(yīng)對(duì)下游城市近地面溫度影響較小,但是在混合層中上層,上游城市的熱島對(duì)下游城市影響顯著。白天的風(fēng)場(chǎng)差值表明,在東南風(fēng)作用下,蘇州、無(wú)錫和常州的城市熱島連成一片,構(gòu)成了一個(gè)更強(qiáng)大的熱島環(huán)流。在近地層,風(fēng)場(chǎng)的強(qiáng)度受兩個(gè)因素的影響:城市下墊面粗糙度和城市熱島環(huán)流特征。一方面,城市地區(qū)建筑物密布,粗糙度明顯高于郊區(qū),風(fēng)場(chǎng)受城市地形的阻擋,風(fēng)速下降明顯;另一方面,城市熱島使得風(fēng)場(chǎng)在近地層向市區(qū)輻合,從而抑制熱島高溫向城市下游的輸送。受這兩種因素的共同影響,近地層的東南風(fēng)被顯著削弱。而在高空1.5 km上下,風(fēng)場(chǎng)強(qiáng)度受熱島環(huán)流的影響較大,熱島環(huán)流在高空表現(xiàn)為向城市周圍輻散,從而加劇了高空的東南風(fēng)。夜晚近地面城市熱島強(qiáng)度及其對(duì)城市下游的溫度平流作用均明顯強(qiáng)于白天,但夜晚大氣層結(jié)趨于穩(wěn)定,城市熱島平均高度不足200 m。夜間風(fēng)場(chǎng)差值表明,城市下墊面使近地面風(fēng)速顯著下降,熱島環(huán)流特征不明顯。白天,蘇州和常州城市熱島的上方均出現(xiàn)一個(gè)明顯的高空冷層,即城市化過(guò)程在導(dǎo)致近地面氣溫上升的同時(shí)還引起了高空氣溫的下降。劉壽東等(2014)在對(duì)南京市區(qū)和郊區(qū)同步進(jìn)行的邊界層探測(cè)中也發(fā)現(xiàn)了這種高空冷層現(xiàn)象。夜晚這種高空冷層特征不明顯。

    圖6為蘇州、無(wú)錫和常州區(qū)域平均溫度垂直廓線差值(CTL-NUB)及平均邊界層高度的日變化,反映了這三座城市熱島強(qiáng)度的演變過(guò)程。熱島強(qiáng)度在早晨最弱,午后和傍晚較強(qiáng),熱島的垂直發(fā)展高度在午后可達(dá)1 km以上。在無(wú)錫和常州的傍晚均出現(xiàn)了熱島強(qiáng)度和垂直發(fā)展高度突然增加的現(xiàn)象,這顯然不是受當(dāng)?shù)責(zé)釐u強(qiáng)度日變化的影響,而是與上游地區(qū)城市熱島的溫度輸送有關(guān)。通過(guò)比較CTL試驗(yàn)與NUB試驗(yàn)平均城市邊界層高度可以發(fā)現(xiàn):城市熱島效應(yīng)會(huì)引起城市邊界層高度的顯著增加,城市熱島造成蘇州、無(wú)錫和常州的平均邊界層高度分別上升115、155和168 m,上升幅度與熱島強(qiáng)度有關(guān),熱島越強(qiáng)則邊界層高度增加越明顯。

    圖6 蘇州(a)、無(wú)錫(b)和常州(c)區(qū)域平均CTL試驗(yàn)與NUB試驗(yàn)邊界層高度(實(shí)線:CTL,虛線:NUB;單位:km)及其溫度垂直廓線差值(陰影;單位:℃)的時(shí)間變化Fig.6 Variations of area-averaged PBL height(solid line:CTL,dashed line:NUB;units:km) from CTL and NUB experiments,and profiles of area-averaged temperature difference(shadings;units:℃) between CTL and NUB experiments over (a)Suzhou,(b)Wuxi and (c)Changzhou

    3.2 太湖對(duì)其周圍城市熱島的影響

    圖7 白天CTL試驗(yàn)與NTL試驗(yàn)近地面溫度(陰影;單位:℃)和風(fēng)場(chǎng)(箭矢;單位:m·s-1)的差值Fig.7 Differences in averaged near surface temperature(shadings;units:℃) and wind field(arrows;units:m·s-1) between CTL and NTL experiments during daytime

    由于水體比陸地具有更大的比熱容,升溫較慢,在水體與陸地之間形成較大溫差,進(jìn)而引起局地氣流的變化。蘇州、無(wú)錫和常州均位于太湖沿岸,本研究通過(guò)敏感性試驗(yàn)討論太湖的湖風(fēng)效應(yīng)對(duì)這三座城市熱島強(qiáng)度和熱島環(huán)流結(jié)構(gòu)的影響。圖7為CTL試驗(yàn)與NTL試驗(yàn)白天近地面氣溫以及風(fēng)場(chǎng)的差值,它反映了太湖對(duì)局地溫度和風(fēng)場(chǎng)結(jié)構(gòu)的影響。白天太湖表面氣溫比周圍陸地低5~8 ℃,巨大的溫差會(huì)引起強(qiáng)烈的湖風(fēng)效應(yīng),風(fēng)速差值(即湖風(fēng)強(qiáng)度)可達(dá)5 m·s-1。在背景風(fēng)場(chǎng)和湖風(fēng)的共同作用下,太湖對(duì)其周邊和下風(fēng)向地區(qū)的氣溫產(chǎn)生顯著影響,導(dǎo)致下游地區(qū)氣溫下降約2 ℃,其影響范圍可達(dá)40 km左右。蘇州位于太湖上風(fēng)向,背景風(fēng)與湖風(fēng)相互抵消,使得太湖對(duì)蘇州影響很小。而無(wú)錫和常州均位于下風(fēng)向,太湖使得這兩座城市白天平均近地面氣溫分別下降了0.49 ℃和0.52 ℃。

    圖8a為AB剖面上NTL試驗(yàn)與NUB試驗(yàn)白天平均溫度和風(fēng)場(chǎng)的差值。可見(jiàn),將太湖替換為農(nóng)田后,無(wú)錫和常州白天城市熱島強(qiáng)度明顯增加。太湖消失使得無(wú)錫與常州的城市熱島連成一片;而太湖存在時(shí),其白天強(qiáng)盛的湖風(fēng)會(huì)削弱無(wú)錫與常州之間近地面城市熱島的聯(lián)系(圖5a)。圖8b同圖8a,但為CTL試驗(yàn)與NTL試驗(yàn)的差值,反映了太湖對(duì)蘇州—無(wú)錫—常州城市帶上氣溫和風(fēng)場(chǎng)垂直結(jié)構(gòu)的影響。白天,湖風(fēng)帶來(lái)的冷氣團(tuán)會(huì)導(dǎo)致無(wú)錫和常州近地面溫度低于NTL試驗(yàn)。兩座城市上空的風(fēng)場(chǎng)差值以下沉運(yùn)動(dòng)為主,這是因?yàn)槔錃鈭F(tuán)進(jìn)入城市上空產(chǎn)生下沉,削弱了城市熱島效應(yīng)所引起的上升運(yùn)動(dòng)。總而言之,太湖的存在削弱了附近城市的熱島強(qiáng)度以及城市之間熱島效應(yīng)的相互作用,抑制了城市熱島的垂直發(fā)展。

    由無(wú)錫與常州區(qū)域平均溫度廓線差值日變化(圖9a、9b)可見(jiàn),太湖引起這兩座城市白天邊界層內(nèi)熱島強(qiáng)度明顯降低,最大下降幅度接近1 ℃。該過(guò)程主要發(fā)生在中午前后,因?yàn)榇藭r(shí)太湖湖風(fēng)最為強(qiáng)盛。由圖9b可見(jiàn),太湖導(dǎo)致傍晚常州邊界層內(nèi)出現(xiàn)第二次熱島強(qiáng)度下降過(guò)程。前述提及,傍晚該溫度變化特征主要受上游地區(qū)城市熱島溫度平流的影響。太湖削弱了無(wú)錫的城市熱島強(qiáng)度,且湖風(fēng)的風(fēng)向與主導(dǎo)風(fēng)向不一致,抑制了無(wú)錫對(duì)常州的溫度平流,從而導(dǎo)致傍晚時(shí)分CTL試驗(yàn)的熱島溫度平流強(qiáng)度弱于NTL試驗(yàn)。由無(wú)錫和常州CTL試驗(yàn)與NTL試驗(yàn)邊界層高度日變化(圖9a、9b)可以發(fā)現(xiàn),太湖消失會(huì)使無(wú)錫和常州白天邊界層高度顯著上升,無(wú)錫的邊界層高度最大增加約450 m,常州最大增加約670 m。NTL試驗(yàn)引起白天無(wú)錫和常州城市邊界層高度的增加甚至超過(guò)了城市熱島效應(yīng)所引起的邊界層高度增加。

    圖8 AB剖面上NTL試驗(yàn)與NUB試驗(yàn)(a)以及CTL試驗(yàn)與NTL試驗(yàn)(b)白天平均氣溫(陰影;單位:℃)和風(fēng)場(chǎng)(箭矢;單位:m·s-1)的差值(CZ、WX和SZ分別代表常州、無(wú)錫和蘇州)Fig.8 Vertical cross-sections of differences in averaged temperature(color shadings;units:℃) and in-plane flow vectors(arrows;units:m·s-1) (a)between NTL and NUB experiments,and (b)between CTL and NTL experiments alone line AB during daytime(CZ,WX and SZ represent Changzhou,Wuxi and Suzhou,respectively)

    圖9 無(wú)錫(a)和常州(b)區(qū)域平均CTL試驗(yàn)與NTL試驗(yàn)邊界層高度(實(shí)線:CTL,虛線:NTL;單位:km)及溫度垂直廓線差值(陰影;單位:℃)的時(shí)間變化Fig.9 Variations of area-averaged PBL height(solid line:CTL,dashed line:NTL;units:km) from CTL and NTL experiments,and profiles of area-averaged temperature difference(shadings;units:℃) between CTL and NTL experiments over (a)Wuxi and (b)Changzhou

    4 結(jié)論

    本文采用中尺度氣象模式WRF耦合一個(gè)多層城市冠層模式BEP,針對(duì)2013年8月13日一次高溫過(guò)程進(jìn)行模擬;設(shè)計(jì)3套試驗(yàn)方案(CTL、NUB、NTL),討論城市化過(guò)程及太湖的湖風(fēng)效應(yīng)對(duì)蘇州—無(wú)錫—常州這一城市帶熱島效應(yīng)的影響,得到如下結(jié)論:

    蘇州—無(wú)錫—常州城市帶熱島效應(yīng)明顯,部分地區(qū)熱島強(qiáng)度達(dá)3 ℃。在東南風(fēng)作用下,城市熱島強(qiáng)度向下游逐漸增加,蘇州熱島強(qiáng)度最弱,常州最強(qiáng)。城市熱島對(duì)其下風(fēng)向地區(qū)有顯著增溫效應(yīng)。通過(guò)比較AB剖面上CTL試驗(yàn)與NUB試驗(yàn)氣溫以及風(fēng)場(chǎng)的差值發(fā)現(xiàn),白天,在東南風(fēng)影響下,這三座城市的熱島效應(yīng)連成一片,形成了一個(gè)更強(qiáng)大的熱島環(huán)流;夜晚,邊界層逐漸趨于穩(wěn)定,熱島環(huán)流減弱,有利于熱島高溫向下游地區(qū)的輸送。受上游城市的影響,無(wú)錫和常州傍晚的城市熱島強(qiáng)度和垂直發(fā)展高度都有顯著增加。城市化使得蘇州、無(wú)錫和常州的平均邊界層高度分別升高了115、155和168 m。

    白天太湖表面氣溫比周圍陸地低5~8 ℃,形成強(qiáng)盛的湖風(fēng),進(jìn)而影響周邊城市白天的熱島效應(yīng)。太湖導(dǎo)致無(wú)錫和常州白天近地面平均氣溫分別下降0.49和0.52 ℃,削弱了其白天城市熱島強(qiáng)度以及這兩座城市之間城市熱島的聯(lián)系,且抑制了城市熱島的發(fā)展高度。太湖使無(wú)錫白天邊界層高度最多下降約450 m,常州下降約670 m。

    江志紅,葉麗梅.2010.近十年南京城市熱島演變的遙感研究[J].南京信息工程大學(xué)學(xué)報(bào):自然科學(xué)版,2(2):148-154.

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

    劉壽東,姜潤(rùn),王成剛,等.2014.南京夏季城市熱島時(shí)空分布特征的觀測(cè)分析[J].大氣科學(xué)學(xué)報(bào),37(1):19-27.

    邱新法,顧麗華,曾燕,等.2008.南京城市熱島效應(yīng)研究[J].氣候與環(huán)境研究,13(6):807-814.

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

    談建國(guó),鄭有飛,彭麗,等.2008.城市熱島對(duì)上海夏季高溫?zé)崂说挠绊慬J].高原氣象,27(增刊):144-149.

    朱家其,湯緒,江灝.2006.上海市城區(qū)氣溫變化及城市熱島[J].高原氣象,25(6):1154-1160.

    Allwine K J,Shinn J H,Streit G E,et al.2002.Overview of URBAN 2000:A multiscale field study of dispersion through an urban environment[J].Bull Amer Meteor Soc,83:521-536.

    Chen Y,Jiang W M,Zhang N,et al.2009.Numerical simulation of the anthropogenic heat effect on urban boundary layer structure[J].Theor Appl Climatol,97:123-134.

    Chin H N S,Leach M J,Sugiyama G A.2005.Evaluation of an urban canopy parameterization in a mesoscale model using VTMX and URBAN 2000 Data[J].Mon Wea Rev,133(7):2043-2068.

    Du Y,Xie Z Q,Zeng Y,et al.2007.Impact of urban expansion on regional temperature change in the Yangtze River Delta[J].J Geogr Sci,17(4):387-398.doi:10.1007/sll442-007-0387-0.

    Dupont S,Otte T L,Ching J K S.2004.Simulation of meteorological fields within and above urban and rural canopies with a mesoscale model[J].Bound-Layer Meteor,113:111-158.

    Fan H,Sailor D J.2005.Modeling the impacts of anthropogenic heating on the urban climate of Philadelphia:A comparison of implementations in two PBL schemes[J].Atmos Environ,39:73-84.

    Houghton J T,Ding Y,Griggs D J,et al.2001.Climate change 2001:The scientific basis[M].Cambridge:Cambridge University Press.

    Ichinose T,Shimodozono K,Hanaki K.1999.Impact of anthropogenic heat on urban climate in Tokyo[J].Atmos Environ,33:3897-3909.

    Kusaka H,Kimura F.2004.Thermal effects of urban canyon structure on the nocturnal heat island:Numerical experiment using a mesoscale model coupled with an urban canopy model[J].J Appl Meteor,43:1899-1910.

    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.

    Martilli A,Clappier A,Rotach M W.2002.An urban surface exchange parameterization for mesoscale models[J].Bound-Layer Meteor,104:261-304.

    Masson V.2000.A physically-based scheme for the urban energy budget in atmospheric models[J].Bound-Layer Meteor,94:357-397.

    Niyogi D,Holt T,Zhong S.2006.Urban and land surface effects on the 30 July 2003 MCS event observed in the Southern Great Plains[J].J Geophys Res,111,D19107.doi:10.1029/2005JD006746.

    Ohashi Y,Genchi Y,Kondo H,et al.2007.Influence of air-conditioning waste heat on air temperature in Tokyo during summer:Numerical experiments using an urban canopy model coupled with a building energy model [J].J Appl Meteor Climatol,46:66-81.

    Oke T A.1982.The energetic basis of the urban heat island[J].Quart J Roy Meteor Soc,108:1-24.

    Peterson T C.2003.Assessment of urban versus rural in situ surface temperatures in the contiguous United States:No difference found[J].J Climate,16(18):2941-2959.

    Rotach M W,Vogt R,Bernhofer C,et al.2005.BUBBLE——An urban boundary layer meteorology project[J].Theor Appl Climatol,81:231-261.

    Rozoff C,Cotton W R,Adegoke J O.2003.Simulation of St Louis,Missouri,land use impacts on thunderstorms[J].J Appl Meteor,42:716-738.

    Ryu Y H,Baik J J,Lee S H.2013.Effects of anthropogenic heat on ozone air quality in a megacity[J].Atmos Environ,80:20-30.

    Tesche T W,Mcnally D E,Emery E T.2001.Evaluation of the MM5 model over the midwestern U.S. for three 8-hour oxidant episodes[R]//Prepared for the Kansas City ozone technical workgroup,by Alpine Geophysics,LLC,Ft. Wright,KY,and ENVIRON International Corp.,Novato,CA,USA.

    Van Weverberg K,De Ridder K,Van Rompaey A.2008.Modeling the contribution of the Brussels heat island to a long temperature time series[J].J Climate Appl Meteor,47:976-990.

    Zhang D L,Shou Y X,Dickerson R R.2009.Upstream urbanization exacerbates urban heat island effects[J].Geophys Res Lett,36,L24401.doi:10.1029/2009GL041082.

    Zhang D L,Shou Y X,Dickerson R R,et al. 2011.Impact of upstream urbanization on the urban heat island effects along the Washington-Baltimore corridor[J].J Appl Meteor Climatol,50:2012-2029.

    Zhang K X,Wang R,Shen C C,et al.2010.Temporal and spatial characteristics of the urban heat island during rapid urbanization in Shanghai,China[J].Environ Monit Assess,169:101-112.

    Zhang N,Zhu L F,Zhu Y.2011.Urban heat island and boundary layer structures under hot weather synoptic conditions:A case study of Suzhou City,China[J].Adv Atmos Sci,28(4):855-865.

    (責(zé)任編輯:張福穎)

    InvestigationofanurbanheatislandepisodealongSuzhou-Wuxi-Changzhouurbancluster

    KANG Han-qing1,2,ZHU Bin1,2,ZHU Tong3,4,GAO Jin-hui1,2,SUN Jia-li5,SU Ji-feng6

    (1.Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,NUIST,Nanjing 210044,China;2.Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration,NUIST,Nanjing 210044,China;3.CIRA/Colorado State University,Fort Collins,CO,USA 80523-1375;4.NOAA/NESDIS/STAR/JCSDA,College Park,MD,USA 20746;5.Jiangsu Climate Center,Nanjing 210009,China;6.The 61 Squad of the 94857 Unit of People’s Liberation Army,Wuhu 241000,China)

    The weather research and forecasting(WRF) model,coupled with a multi-layer urban canopy model BEP(building energy parameterization),is used to investigate a high temperature episode occurred in the Yangtze River Delta region on 13 August 2013.A southeasterly breeze,consistent with the direction of Suzhou-Wuxi-Changzhou urban cluster,is prevailed over the Yangtze River Delta region during the episode.Simulating results show that significant urban heat island(UHI) effects occur in the urban cluster,and the UHI intensity increases from Suzhou to Changzhou.Under the impact of the southeasterly breeze,the UHI effects of the three individual cities merge together,forming a stronger UHI circulation.In the night,the UHI circulation is weakened by the stable planetary boundary layer,which is favorable for the downstream transport of the UHI effect.The urban planetary boundary layer height is significantly increased by the UHI effect.During the daytime,a strong lake breeze between Lake Taihu and land has great impacts on the UHI of the nearby cities.Air temperature over Wuxi and Changzhou is significantly decreased by the cold air from Lake Taihu,and the upward extension of the UHIs over the two cities are suppressed.Connections of the UHIs between Wuxi and Changzhou are weakened by Lake Taihu.The planetary boundary layer height over Wuxi and Changzhou is more significantly decreased by the lake breeze during daytime.

    urban heat island;numerical simulation;Lake Taihu;atmospheric boundary layer

    2014-03-20;改回日期2014-05-28

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41275143);公益性行業(yè)(氣象)科研專項(xiàng)(201206011);江蘇省高校自然科學(xué)研究重大基礎(chǔ)研究項(xiàng)目(12KJA170003);江蘇省“333”高層次人才培養(yǎng)工程項(xiàng)目;江蘇省“六大人才高峰”計(jì)劃項(xiàng)目

    朱彬,博士,教授,博士生導(dǎo)師,研究方向?yàn)榇髿馕锢韺W(xué)與大氣環(huán)境,binzhu@nuist.edu.cn.

    10.13878/j.cnki.dqkxxb.20140320001.

    1674-7097(2014)04-0432-09

    P461.2

    A

    10.13878/j.cnki.dqkxxb.20140320001

    康漢青,朱彬,朱彤,等.2014.蘇州—無(wú)錫—常州城市帶熱島效應(yīng)個(gè)例研究[J].大氣科學(xué)學(xué)報(bào),37(4):432-440.

    Kang Han-qing,Zhu Bin,Zhu Tong,et al.2014.Investigation of an urban heat island episode along Suzhou-Wuxi-Changzhou urban cluster[J].Trans Atmos Sci,37(4):432-440.(in Chinese)

    人人妻人人澡人人爽人人夜夜| 精品久久国产蜜桃| 内射极品少妇av片p| 十八禁高潮呻吟视频 | 久久婷婷青草| 简卡轻食公司| 高清不卡的av网站| 亚洲,一卡二卡三卡| 亚洲色图综合在线观看| 国产在视频线精品| 男男h啪啪无遮挡| 国产成人91sexporn| 日本午夜av视频| 中文欧美无线码| 国产av一区二区精品久久| 免费在线观看成人毛片| 狂野欧美激情性bbbbbb| 精品人妻熟女av久视频| 一级二级三级毛片免费看| 国产亚洲91精品色在线| 日韩制服骚丝袜av| 日日撸夜夜添| 69精品国产乱码久久久| 一级毛片黄色毛片免费观看视频| 色网站视频免费| 99视频精品全部免费 在线| 成人亚洲欧美一区二区av| 国产精品99久久久久久久久| 欧美老熟妇乱子伦牲交| 91精品一卡2卡3卡4卡| 精品视频人人做人人爽| 久久亚洲国产成人精品v| 久久久欧美国产精品| 伦理电影大哥的女人| 亚洲在久久综合| 99久久精品国产国产毛片| 各种免费的搞黄视频| xxx大片免费视频| 菩萨蛮人人尽说江南好唐韦庄| 赤兔流量卡办理| 一级毛片 在线播放| 69精品国产乱码久久久| 少妇人妻久久综合中文| xxx大片免费视频| 欧美精品人与动牲交sv欧美| 熟女人妻精品中文字幕| 最新中文字幕久久久久| 国产色婷婷99| 狂野欧美激情性bbbbbb| 免费观看性生交大片5| 亚洲精品久久午夜乱码| 日韩成人av中文字幕在线观看| 亚洲色图综合在线观看| 久久久久人妻精品一区果冻| 一区二区三区四区激情视频| 国产亚洲精品久久久com| 国产亚洲一区二区精品| 自线自在国产av| 99久久精品一区二区三区| av天堂中文字幕网| 国产免费福利视频在线观看| 2022亚洲国产成人精品| 亚洲精品久久久久久婷婷小说| 日本vs欧美在线观看视频 | 日韩 亚洲 欧美在线| 热re99久久国产66热| 久久久久久久精品精品| 亚洲综合精品二区| 国产视频内射| 性色av一级| 另类精品久久| 乱码一卡2卡4卡精品| 亚洲国产精品999| 中文字幕人妻丝袜制服| 老司机影院毛片| 久久久国产精品麻豆| 天堂中文最新版在线下载| 午夜av观看不卡| 欧美精品一区二区免费开放| 国产一区二区三区av在线| av免费在线看不卡| 午夜免费鲁丝| 大香蕉97超碰在线| 久久99热6这里只有精品| 亚洲人成网站在线播| 久久婷婷青草| 岛国毛片在线播放| 人妻夜夜爽99麻豆av| 久久国产乱子免费精品| 国产精品秋霞免费鲁丝片| 久久国产乱子免费精品| 亚洲精品日韩av片在线观看| 国内揄拍国产精品人妻在线| 只有这里有精品99| 日本午夜av视频| 伦理电影免费视频| 天堂8中文在线网| 在线观看国产h片| 中国国产av一级| 大片免费播放器 马上看| 日韩中文字幕视频在线看片| 这个男人来自地球电影免费观看 | 狠狠精品人妻久久久久久综合| 亚洲精品国产av蜜桃| 91精品一卡2卡3卡4卡| h日本视频在线播放| 两个人的视频大全免费| 一区在线观看完整版| 久久国产乱子免费精品| 视频区图区小说| 亚洲欧洲日产国产| 欧美性感艳星| 国产成人aa在线观看| 中文在线观看免费www的网站| 精品亚洲乱码少妇综合久久| 久久毛片免费看一区二区三区| 国产精品.久久久| av在线观看视频网站免费| 精品午夜福利在线看| 高清av免费在线| 成人漫画全彩无遮挡| 最近手机中文字幕大全| 麻豆乱淫一区二区| tube8黄色片| 高清视频免费观看一区二区| 人体艺术视频欧美日本| 一级,二级,三级黄色视频| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| 大片电影免费在线观看免费| 亚洲精品久久午夜乱码| 精品一区二区三区视频在线| 国产亚洲一区二区精品| 国产色爽女视频免费观看| 久久精品国产亚洲av涩爱| 在线看a的网站| 国产熟女午夜一区二区三区 | 欧美日韩精品成人综合77777| 日韩中文字幕视频在线看片| 男女国产视频网站| 在线观看免费高清a一片| 热99国产精品久久久久久7| 国产精品久久久久久久久免| av在线app专区| 在线观看免费视频网站a站| 一个人看视频在线观看www免费| 99视频精品全部免费 在线| 久久99精品国语久久久| 亚洲欧洲国产日韩| 精品久久国产蜜桃| 国产av国产精品国产| 男人和女人高潮做爰伦理| 亚洲国产成人一精品久久久| 亚洲av成人精品一区久久| 不卡视频在线观看欧美| 老司机影院毛片| 22中文网久久字幕| 视频区图区小说| 日本av手机在线免费观看| 久久久久国产精品人妻一区二区| 久久久久久久久久久丰满| 夜夜爽夜夜爽视频| 亚洲无线观看免费| 亚洲第一区二区三区不卡| 日日啪夜夜撸| 国产伦精品一区二区三区视频9| 一级毛片电影观看| 少妇猛男粗大的猛烈进出视频| 黄色视频在线播放观看不卡| 亚洲精品乱码久久久v下载方式| 欧美人与善性xxx| 建设人人有责人人尽责人人享有的| 欧美xxxx性猛交bbbb| 99久国产av精品国产电影| 如何舔出高潮| 日韩在线高清观看一区二区三区| 在线精品无人区一区二区三| 久久久久久久久久久久大奶| 亚州av有码| 成人18禁高潮啪啪吃奶动态图 | 永久免费av网站大全| 国产一区二区在线观看日韩| av天堂久久9| 久久人人爽人人片av| av黄色大香蕉| 亚洲av二区三区四区| 欧美日韩视频精品一区| 亚洲自偷自拍三级| 久久97久久精品| 亚洲国产精品成人久久小说| 国产精品久久久久久久久免| 在线观看免费视频网站a站| 人人妻人人澡人人爽人人夜夜| 亚洲真实伦在线观看| 国产成人精品婷婷| 卡戴珊不雅视频在线播放| 一本久久精品| a级毛色黄片| av在线老鸭窝| 99精国产麻豆久久婷婷| 亚洲美女视频黄频| 我要看日韩黄色一级片| 免费观看的影片在线观看| 亚洲精品成人av观看孕妇| 一级av片app| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 一级毛片 在线播放| 欧美日韩国产mv在线观看视频| 国产精品国产三级国产专区5o| 大又大粗又爽又黄少妇毛片口| 久久久久国产网址| 哪个播放器可以免费观看大片| 久久人人爽av亚洲精品天堂| 亚洲av免费高清在线观看| 新久久久久国产一级毛片| 精品少妇黑人巨大在线播放| 男女免费视频国产| 菩萨蛮人人尽说江南好唐韦庄| 久久午夜综合久久蜜桃| 国产黄色视频一区二区在线观看| 精品久久久噜噜| 久久久国产精品麻豆| 亚洲av男天堂| 春色校园在线视频观看| 91成人精品电影| 黄色日韩在线| 日韩中字成人| 校园人妻丝袜中文字幕| 亚洲国产欧美日韩在线播放 | 欧美激情极品国产一区二区三区 | 久久国产精品大桥未久av | 一级毛片久久久久久久久女| 少妇裸体淫交视频免费看高清| 一二三四中文在线观看免费高清| 精品国产国语对白av| 韩国高清视频一区二区三区| 又黄又爽又刺激的免费视频.| 久久久久久久久久人人人人人人| 免费播放大片免费观看视频在线观看| 嫩草影院入口| 丝瓜视频免费看黄片| 国产伦精品一区二区三区四那| 中文字幕人妻丝袜制服| 欧美日韩视频精品一区| 精品午夜福利在线看| 中文字幕av电影在线播放| 国产精品一区www在线观看| 亚洲精品乱码久久久v下载方式| 国产午夜精品一二区理论片| 精品熟女少妇av免费看| 国产 一区精品| av天堂久久9| 亚洲欧美中文字幕日韩二区| 亚洲第一av免费看| 少妇高潮的动态图| 久久久午夜欧美精品| 日本猛色少妇xxxxx猛交久久| 精品久久久精品久久久| 亚洲国产毛片av蜜桃av| 日韩熟女老妇一区二区性免费视频| 久久久午夜欧美精品| 夫妻午夜视频| 久久精品国产亚洲网站| 高清毛片免费看| 丝袜在线中文字幕| 99视频精品全部免费 在线| 黑人猛操日本美女一级片| 欧美丝袜亚洲另类| 99精国产麻豆久久婷婷| 午夜久久久在线观看| 秋霞在线观看毛片| 亚洲性久久影院| av福利片在线| 精品久久久精品久久久| 国产精品国产av在线观看| 久久99精品国语久久久| 18禁在线无遮挡免费观看视频| 成年av动漫网址| 大话2 男鬼变身卡| 亚洲精品久久午夜乱码| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 中文字幕av电影在线播放| 精品一区二区三卡| 免费观看a级毛片全部| 2021少妇久久久久久久久久久| 妹子高潮喷水视频| 国产伦精品一区二区三区四那| 自拍偷自拍亚洲精品老妇| 大片免费播放器 马上看| h视频一区二区三区| 国产精品成人在线| 亚洲成人手机| 亚洲图色成人| 久久久久久久久久久久大奶| 久久人人爽人人爽人人片va| 国产视频内射| av福利片在线观看| 男女国产视频网站| 亚洲图色成人| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 亚洲精品456在线播放app| 99视频精品全部免费 在线| 欧美国产精品一级二级三级 | 日本黄色日本黄色录像| 日本wwww免费看| 亚洲va在线va天堂va国产| av视频免费观看在线观看| 午夜av观看不卡| 日韩免费高清中文字幕av| 欧美精品高潮呻吟av久久| 成人毛片a级毛片在线播放| 国产一区二区在线观看av| 精品少妇久久久久久888优播| 精品99又大又爽又粗少妇毛片| 国产伦在线观看视频一区| 一区在线观看完整版| 一区二区三区乱码不卡18| 国产一区有黄有色的免费视频| 亚洲激情五月婷婷啪啪| 午夜免费鲁丝| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| 日产精品乱码卡一卡2卡三| 久久人人爽av亚洲精品天堂| 中文字幕制服av| 欧美+日韩+精品| 欧美亚洲 丝袜 人妻 在线| 国产综合精华液| 久热久热在线精品观看| av网站免费在线观看视频| 成人影院久久| 乱码一卡2卡4卡精品| 亚洲av国产av综合av卡| 偷拍熟女少妇极品色| 亚洲av不卡在线观看| 看十八女毛片水多多多| 国产真实伦视频高清在线观看| 22中文网久久字幕| 日日摸夜夜添夜夜爱| 免费黄频网站在线观看国产| 免费av不卡在线播放| 亚洲av福利一区| 99热6这里只有精品| 国产在线免费精品| 久久久久久伊人网av| 久久精品夜色国产| 久久久久国产网址| 亚洲国产欧美在线一区| 99热这里只有精品一区| 永久免费av网站大全| 黄色毛片三级朝国网站 | 亚洲精品aⅴ在线观看| 国产综合精华液| 观看免费一级毛片| 亚洲精品乱码久久久久久按摩| 人人澡人人妻人| 中文字幕免费在线视频6| 99视频精品全部免费 在线| 久久久久人妻精品一区果冻| 特大巨黑吊av在线直播| 精品国产一区二区三区久久久樱花| 又粗又硬又长又爽又黄的视频| 另类精品久久| 熟女电影av网| 午夜久久久在线观看| 亚洲av综合色区一区| 欧美最新免费一区二区三区| 又粗又硬又长又爽又黄的视频| 亚洲精品第二区| 岛国毛片在线播放| 亚州av有码| 超碰97精品在线观看| 夫妻性生交免费视频一级片| 国产亚洲91精品色在线| 性色av一级| .国产精品久久| 精品久久久久久久久av| 一级av片app| 成人毛片a级毛片在线播放| 美女xxoo啪啪120秒动态图| 丝瓜视频免费看黄片| 中文天堂在线官网| 久久国产亚洲av麻豆专区| av专区在线播放| 成人黄色视频免费在线看| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看| 久久99蜜桃精品久久| 一本久久精品| 啦啦啦啦在线视频资源| 日韩一区二区三区影片| 高清视频免费观看一区二区| 多毛熟女@视频| av线在线观看网站| 亚洲精品亚洲一区二区| 亚洲国产精品一区三区| av在线app专区| 亚洲精品久久久久久婷婷小说| 日本爱情动作片www.在线观看| 久久久久久久亚洲中文字幕| 黄片无遮挡物在线观看| 午夜91福利影院| 精品视频人人做人人爽| 亚洲欧美成人综合另类久久久| 亚洲精品第二区| 精品国产露脸久久av麻豆| 久久久久久久精品精品| 精品久久久噜噜| 日本vs欧美在线观看视频 | 中文字幕久久专区| 日韩欧美一区视频在线观看 | 欧美日韩国产mv在线观看视频| 秋霞伦理黄片| 日韩伦理黄色片| av卡一久久| 中文字幕久久专区| 国产成人精品无人区| 成人18禁高潮啪啪吃奶动态图 | 国产在视频线精品| a级片在线免费高清观看视频| 日本wwww免费看| 曰老女人黄片| 高清视频免费观看一区二区| 久久久国产一区二区| 97在线视频观看| 国产精品久久久久久精品古装| 人体艺术视频欧美日本| 欧美另类一区| 亚洲国产成人一精品久久久| 午夜福利在线观看免费完整高清在| 成人18禁高潮啪啪吃奶动态图 | 国产成人免费无遮挡视频| videos熟女内射| 亚洲精品视频女| 性色avwww在线观看| 人人妻人人看人人澡| 97超视频在线观看视频| 中文字幕制服av| 精品人妻偷拍中文字幕| 国精品久久久久久国模美| 亚洲第一av免费看| 少妇丰满av| 国产老妇伦熟女老妇高清| 国产一区二区三区综合在线观看 | 亚洲精品乱码久久久v下载方式| 亚洲av中文av极速乱| 一本久久精品| 国产毛片在线视频| 国产av一区二区精品久久| 午夜视频国产福利| 在线观看三级黄色| 亚洲av在线观看美女高潮| 91aial.com中文字幕在线观看| 色婷婷久久久亚洲欧美| 欧美最新免费一区二区三区| 久久精品久久精品一区二区三区| 国产成人午夜福利电影在线观看| 777米奇影视久久| 黄色欧美视频在线观看| 一级二级三级毛片免费看| 亚洲国产日韩一区二区| 午夜福利影视在线免费观看| 一级,二级,三级黄色视频| 美女视频免费永久观看网站| 色婷婷av一区二区三区视频| 午夜免费观看性视频| 三级国产精品欧美在线观看| 亚洲国产欧美在线一区| 乱人伦中国视频| 美女主播在线视频| 超碰97精品在线观看| 这个男人来自地球电影免费观看 | 桃花免费在线播放| 伊人久久国产一区二区| 最新中文字幕久久久久| 极品教师在线视频| 日本av免费视频播放| 久久精品久久久久久久性| 亚洲精品久久久久久婷婷小说| 精品人妻熟女av久视频| 人妻夜夜爽99麻豆av| 夜夜骑夜夜射夜夜干| 成年美女黄网站色视频大全免费 | 蜜桃在线观看..| 一二三四中文在线观看免费高清| 国产精品一区二区性色av| 国产亚洲一区二区精品| 99热网站在线观看| 成人午夜精彩视频在线观看| 啦啦啦视频在线资源免费观看| 肉色欧美久久久久久久蜜桃| 午夜激情福利司机影院| 草草在线视频免费看| 日本欧美视频一区| 人人妻人人添人人爽欧美一区卜| 精品国产露脸久久av麻豆| 高清欧美精品videossex| 全区人妻精品视频| 人妻人人澡人人爽人人| 黄色视频在线播放观看不卡| 99国产精品免费福利视频| 成人漫画全彩无遮挡| 中国美白少妇内射xxxbb| 久久久久人妻精品一区果冻| 国产一区有黄有色的免费视频| 国产一区二区在线观看日韩| 韩国高清视频一区二区三区| 九九在线视频观看精品| 精品卡一卡二卡四卡免费| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看| 中国美白少妇内射xxxbb| 中文精品一卡2卡3卡4更新| 欧美成人午夜免费资源| 亚洲欧洲国产日韩| 丰满饥渴人妻一区二区三| 午夜福利在线观看免费完整高清在| 亚洲高清免费不卡视频| 免费看不卡的av| 一区二区av电影网| 亚州av有码| 日韩强制内射视频| 午夜免费鲁丝| 久久久久久久亚洲中文字幕| 精品少妇黑人巨大在线播放| 国产成人精品婷婷| 老女人水多毛片| 亚洲精品日韩av片在线观看| 国产成人午夜福利电影在线观看| av在线播放精品| 亚洲av在线观看美女高潮| 久久久亚洲精品成人影院| 午夜福利影视在线免费观看| 久久av网站| 精品久久国产蜜桃| 精品亚洲乱码少妇综合久久| 久久99精品国语久久久| 亚洲四区av| 国产精品不卡视频一区二区| 18禁裸乳无遮挡动漫免费视频| 亚洲美女搞黄在线观看| 2022亚洲国产成人精品| 国产老妇伦熟女老妇高清| 国产精品不卡视频一区二区| 日韩av免费高清视频| 欧美成人午夜免费资源| 国产一区二区三区av在线| 国产淫片久久久久久久久| 老司机影院毛片| 久久99热这里只频精品6学生| 免费人成在线观看视频色| 国产成人免费观看mmmm| 久久久久久久大尺度免费视频| 一级,二级,三级黄色视频| 午夜老司机福利剧场| 伦精品一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 麻豆成人午夜福利视频| 国产有黄有色有爽视频| 国产成人a∨麻豆精品| 午夜福利,免费看| av一本久久久久| 亚洲天堂av无毛| 韩国高清视频一区二区三区| 边亲边吃奶的免费视频| 国产 精品1| 亚洲精品日韩在线中文字幕| 男人和女人高潮做爰伦理| 免费av不卡在线播放| 91在线精品国自产拍蜜月| 久久亚洲国产成人精品v| 亚洲欧洲日产国产| 男人添女人高潮全过程视频| 国产精品国产三级国产av玫瑰| 亚洲不卡免费看| 亚洲欧美成人综合另类久久久| 丝袜在线中文字幕| 如日韩欧美国产精品一区二区三区 | 亚洲成人一二三区av| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影小说| 国产av精品麻豆| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看| 久久 成人 亚洲| 亚洲中文av在线| 91在线精品国自产拍蜜月| 欧美激情极品国产一区二区三区 | 80岁老熟妇乱子伦牲交| 国产精品免费大片| 精品视频人人做人人爽| 夫妻午夜视频| 国产伦精品一区二区三区四那| 国产精品免费大片| 精品视频人人做人人爽| 一级a做视频免费观看| 男人添女人高潮全过程视频| 国产亚洲91精品色在线| 国产淫语在线视频| 免费av中文字幕在线| 国产在线男女| 国产精品欧美亚洲77777| 蜜桃在线观看..| 日韩强制内射视频| 伦精品一区二区三区| 亚洲色图综合在线观看| 日韩熟女老妇一区二区性免费视频| 久久 成人 亚洲| av天堂中文字幕网| 一区二区三区免费毛片| 国产av码专区亚洲av| 春色校园在线视频观看|