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

    中國東部氣溶膠在天氣尺度上的輻射強迫和對地面氣溫的影響

    2015-12-05 07:49:03廖禮漏嗣佳符瑜常文淵廖宏
    大氣科學(xué) 2015年1期
    關(guān)鍵詞:效應(yīng)區(qū)域

    廖禮 漏嗣佳 符瑜 常文淵 廖宏

    ?

    中國東部氣溶膠在天氣尺度上的輻射強迫和對地面氣溫的影響

    廖禮1, 2漏嗣佳1, 2符瑜3常文淵1廖宏1

    1中國科學(xué)院大氣物理研究所大氣邊界層物理和大氣化學(xué)國家重點實驗室(LAPC),北京100029;2中國科學(xué)院大學(xué),北京100049;3中國科學(xué)院大氣物理研究所氣候變化研究中心(CCRC),北京100029

    本文應(yīng)用WRF-Chem(Weather Research and Forecasting—Chemistry)模式研究中國東部地區(qū)氣溶膠及其部分組分(硫酸鹽、硝酸鹽和黑碳氣溶膠)在天氣尺度下的輻射強迫和對地面氣溫的影響。5個無明顯降水時間段(2006年8月23~25日、2008年11月10~12日、2008年12月16~18日、2009年1月15~17日和2009年4月27~29日)的模擬顯示,氣溶膠濃度呈現(xiàn)顯著的白天低,夜間高的日變化特征,且北方區(qū)域(29.8°~42.6°N,110.2°~120.3°E)平均PM2.5近地面濃度(40~80 μg m?3)高于南方區(qū)域(22.3°~29.9°N,109.7°~120.2°E,30~47 μg m?3)。氣溶膠對地面2 m溫度(地面氣溫)有明顯的降溫效果,在早上08:00(北京時,下同)和下午17:00左右最為顯著,最高可降低約0.2~1 K,同時氣溶膠的參與改善了模式對地面氣溫的模擬。本文還通過對2006年8月23~25日一次個例的模擬,定量分析了氣溶膠及其部分組分(硫酸鹽、硝酸鹽和黑碳氣溶膠)的總天氣效應(yīng)(直接效應(yīng)+間接效應(yīng))、直接效應(yīng)和間接效應(yīng)分別對到達地面的短波輻射和地面氣溫的影響。北方區(qū)域平均氣溶膠直接效應(yīng)所造成的短波輻射強迫要高于南方區(qū)域,分別為-11.3 W m?2和-5.8 W m?2,導(dǎo)致地面氣溫分別降低了0.074 K和0.039 K。南方區(qū)域平均氣溶膠間接效應(yīng)所產(chǎn)的短波輻射強迫高于北方區(qū)域,分別為-14.4 W m?2和-12.4 W m?2,引起的地面氣溫的改變分別為-0.094 K和-0.035 K。對于氣溶膠組分,硫酸鹽氣溶膠的直接效應(yīng)和間接效應(yīng)的作用相當,其總效應(yīng)在北方和南方區(qū)域平均短波輻射強迫分別為-7.0 W m?2和-10.5 W m?2,對地面氣溫的影響為-0.062 K和-0.074 K,而硝酸鹽氣溶膠的作用略小。黑碳氣溶膠使得北方和南方區(qū)域平均到達地表的太陽短波輻射分別減少了6.5 W m?2和5.8 W m?2,而地表氣溫則分別增加了0.053 K和0.017 K,相比于間接效應(yīng),黑碳氣溶膠的直接效應(yīng)的影響更加顯著。

    中國東部 氣溶膠 輻射強迫 地面氣溫

    1 引言

    氣溶膠通過其直接和間接效應(yīng)影響地氣系統(tǒng)的能量平衡:氣溶膠吸收和散射太陽輻射,從而減少到達地面的太陽輻射,引起地面氣溫的降低,稱為氣溶膠的直接效應(yīng);氣溶膠還可以形成云凝結(jié)核,從而改變云的輻射特性和生命周期,稱為氣溶膠的間接效應(yīng)。IPCC(2007)指出,氣溶膠直接和間接輻射強迫總作用是減少到達近地面的太陽輻射,形成冷卻效應(yīng)。

    目前已有很多應(yīng)用數(shù)值模擬研究氣溶膠通過輻射強迫對氣候的影響。Takemura et al.(2005)應(yīng)用全球三維氣溶膠傳輸輻射模式SPRINTARS(Spectral Radiation-Transport Model for Aerosol Species)模擬了1850和2000年氣溶膠輻射強迫,指出在對流層頂人為氣溶膠全球平均直接和間接輻射強迫分別為-0.1 W m?2和-0.9 W m?2,氣溶膠緩解了約40%由于人為活動產(chǎn)生的溫室氣體所引起的地面空氣溫度的升高。Zhang et al.(2010a)應(yīng)用WRF-Chem(Weather Research and Forecasting— Chemistry)模式模擬2001年美國大陸地區(qū)氣溶膠的氣候效應(yīng),并定量分析了氣溶膠分別能夠減少1月和7月的太陽輻射多達9%和16%,進而使得1月和7月地面2 m溫度降低了多達0.16 K和0.37 K。不同種類的氣溶膠由于其散射和吸收特性的不同而引起不同的輻射效應(yīng),對太陽輻射和地面氣溫的影響也存在差異。硫酸鹽和硝酸鹽氣溶膠屬于散射型氣溶膠,能夠吸收和散射太陽輻射,從而產(chǎn)生降溫效應(yīng)。Giorgi et al.(2002)應(yīng)用RegCM(Regional Climate Model)模式模擬中國四川盆地硫酸鹽氣溶膠表明,1993~1997年間人為硫酸鹽氣溶膠的輻射強迫作用可使地面空氣溫度降低0.1~0.7 K。Li et al.(2009)利用RegCM和TACM(Tropospheric Atmosphere Chemistry Model)相耦合的模式模擬發(fā)現(xiàn),2003年1月和7月中國地區(qū)人為硝酸鹽氣溶膠在對流層頂?shù)钠骄g接輻射強迫分別為-1.63 W m?2和-2.65 W m?2,尤其在一些區(qū)域甚至達到-10 W m?2,對1月和7月地面空氣溫度分別降低了0.13 K和0.09 K。與硫酸鹽和硝酸鹽氣溶膠(散射性氣溶膠)不同,懸浮于大氣中的黑碳氣溶膠(吸收性氣溶膠)能夠吸收太陽輻射從而有很強的輻射加熱效應(yīng)(Jacobson,2001)。Wang(2004)應(yīng)用CCM3(Community Climate Model version 3)模式研究發(fā)現(xiàn)黑碳氣溶膠對太陽輻射的強烈吸收作用,導(dǎo)致黑碳氣溶膠在大氣頂?shù)?0年全球平均直接輻射強迫為0.16~0.80 W m?2,改變地面空氣溫度 0.09±0.07 K。由以上這些研究可以看出,目前的大多數(shù)相關(guān)研究主要集中在氣候研究方面,進行定量分析研究氣溶膠的輻射效應(yīng)及其對溫度的影響。然而在天氣尺度上,氣溶膠的影響和反饋作用也是十分重要的。

    目前氣溶膠對太陽輻射和地面氣溫在天氣尺度上影響的研究相對較少。一些研究中將氣溶膠考慮到天氣預(yù)報模式中評估氣溶膠對天氣預(yù)報的影響。Pérez et al.(2006)在DREAM(Dust Regional Atmospheric Modeling System)模式中加入沙塵的影響,使模式對大氣溫度和平均海表氣壓的預(yù)測有很大的提高。Rodwell and Jung(2008)在歐洲中心中尺度天氣預(yù)報模式ECMWF(European Centre for Medium-range Weather Forecasts)中加入氣溶膠的模擬,結(jié)果顯示加入氣溶膠的模擬能夠改善模式對氣溫等氣象參數(shù)的預(yù)報。此外,Zhang(2008)應(yīng)用WRF-Chem模式研究2000年8月德克薩斯州地區(qū)氣溶膠反饋作用,5天模擬結(jié)果顯示,氣溶膠的存在導(dǎo)致近地面氣溫降低最高達0.18 K。其后續(xù)的研究進一步表明大氣中的PM2.5(大氣中直徑小于或等于2.5 μm的顆粒物)的存在,可減少到達地面的太陽輻射11.2~14.4 W m?2(4.1%~5.6%),從而使近地面氣溫降低0.06~0.14 K(0.2%~0.4%)(Zhang et al., 2010b)。Wang et al.(2009)應(yīng)用SBDART(Santa Barbara Discrete Ordinate Radiative Transfere)模式研究北京地區(qū)在不同天氣狀況下的氣溶膠輻射強迫,即干凈天(2004年10月26日)、灰霾天(2004年10月29日)、霧天(2004年11月30日)和沙塵天(2006年04月17日),近地面PM10濃度分別為:38.1 μg m?3、320.8 μg m?3、229.2 μg m?3和197.7μg m?3,對應(yīng)產(chǎn)生的氣溶膠輻射強迫分別為:-30 W m?2、-76 W m?2、-72 W m?2和-136 W m?2,這也表明了氣溶膠濃度越高,產(chǎn)生的輻射強迫作用更大。從全球來看,中國地區(qū)是氣溶膠濃度高值區(qū),因此所產(chǎn)生的輻射強迫以及對地面氣溫的影響要高于美國地區(qū)。Zhang et al.(2009)利用WRF-Chem模式模擬2008年1月26~30日的中國南方暴雪過程中黑碳氣溶膠的影響,結(jié)果表明黑碳氣溶膠的沉降降低雪表面反照率,使研究區(qū)域平均地面氣溫最高增加0.033 K。然而以往的氣溶膠在天氣尺度上的輻射強迫和對地面氣溫的影響的研究多集中在對單個時間段個例的研究,對中國東部地區(qū)氣溶膠及其組分的總效應(yīng)、直接效應(yīng)和間接效應(yīng)所產(chǎn)生的輻射強迫及其對溫度的影響缺少全面系統(tǒng)的定量分析。

    本文應(yīng)用WRF-Chem模式研究中國東部地區(qū)氣溶膠及其部分組分(硫酸鹽、硝酸鹽和黑碳氣溶膠)在天氣尺度下的輻射強迫和對地面氣溫的影響。通過選取了多個時間段的個例(2006年08月23~25日、2008年11月10~12日、2008年12月16~18日、2009年01月15~17日和2009年04月27~29日)進行研究,從統(tǒng)計的角度進行定量分析,并具體分析其中的一個時間段的個例(2006年8月23~25日)中氣溶膠及其部分組分(硫酸鹽、硝酸鹽和黑碳氣溶膠)的總效應(yīng)、直接效應(yīng)和間接效應(yīng)的輻射強迫和對地面氣溫的影響。本文第二部分介紹模式的設(shè)置及使用數(shù)據(jù)的來源;第三部分簡要介紹了模式試驗的設(shè)計;模式驗證及模擬結(jié)果分析分別在本文的第四和第五部分給出;第六部分為結(jié)論。

    2 模式設(shè)置和數(shù)據(jù)

    WRF(Weather Research and Forecasting model)模式是NCAR(National Center for Atmospheric Research)及NCEP(National Centers for Environmental Prediction)等研究開發(fā)的中尺度數(shù)值天氣預(yù)報模式及同化系統(tǒng),加入化學(xué)模塊后,對于大氣化學(xué)方面的研究該模式又有了很大的發(fā)展,并使其能夠更好對于氣象—大氣氣溶膠—輻射傳輸—云之間的反饋作用進行模擬(Grell et al., 2005; Fast et al., 2006; Gustafson et al., 2007; Chapman et al., 2009)。

    本文采用WRF-Chem 3.3版本,該版本發(fā)布于2011年4月。模式所選的主要的物理參數(shù)化方案包括以下:改進的Lin et al.(1983)微物理方案, RRTM(The Rapid Radiative Transfer Model)長波輻射方案和Goddard短波輻射方案(Mlawer et al., 1997),YSU(Yonsei University)邊界層方案(Hong et al., 2006),NOAH(National Center for Environmental Prediction,Oregon State University,Air Force,and Hydrologic Research Lab’s)陸面模塊(Chen and Dudhia, 2001; Ek et al., 2003),Grell-Devenyi積云參數(shù)化方案(Grell and Dévényi, 2002)。所選的主要的化學(xué)參數(shù)化方案包括以下:光解率采用Fast-J方案(Wild et al., 2000),氣相化學(xué)反應(yīng)機制采用CBM-Z(Carbon-Bond Mechanism version Z)方案(Zaveri and Peters,1999),氣溶膠模塊采用MOSIC(Model for Simulating Aerosol Interactions and Chemistry)方案(Zaveri et al., 2008),MOSIC粒徑分為8個檔(0.039~0.078、0.078~0.156、0.156~0.3125、0.3125~0.625、0.625~1.25、1.25~2.5、2.5~5.0、5.0~10 μm),液相化學(xué)反應(yīng)方案采用Fahey and Pandis(2001)。

    本文中WRF-Chem模式模擬區(qū)域覆蓋中國大部分區(qū)域,蒙古、日本、朝鮮半島以及亞洲中部和俄羅斯的部分地區(qū),中心位于北京附近,采用Lambert投影,模式水平分辨率為60 km×60 km。從10 hPa高度開始模式垂直方向分為30層。初始氣象場和邊界條件都取自于NCEP-FNL(National Centers for Environmental Prediction–Final Analysis)全球分析資料(http://www.cgrer.uiowa.edu/EMISSION_ DATA_new/index_16.html[2013-10-05]),此資料時次間隔為6小時,水平分辨率為1°×1°,使用這些資料可以使得WRF-Chem模式模擬出更加合理和準確的氣象條件。化學(xué)模塊中沒有自帶的排放 源,只有一些背景值,因此在使用模式時要添加模擬需要的排放清單。本文使用David Streets 2006年的亞洲排放源清單INTEX-B(Intercontinental Chemical Transport Experiment-Phase B,http://mic. greenresource. cn/intex-b2006 [2013-10-08])。該源的分辨率為0.5°×0.5°,包含的物種有黑碳(BC)、一氧化碳(CO)、氮氧化物(NO)、有機碳(OC)、二氧化硫(SO2)、揮發(fā)性有機物(VOC)。其中,VOC包括烷烴,烯烴、芳香烴等30個物種。氨氣(NH3)排放選用David Streets 2000年亞洲排放清單,分辨率為0.5°×0.5°。VOC的地表植被排放根據(jù)模擬的氣象條件,使用Gunther方案在線計算獲得(Guenther et al., 1993;Guenther et al., 1994),沒有加入海鹽和沙塵排放。表1給出了中國東部和中國區(qū)域排放總量,由表1可見50%以上排放集中在了中國東部地區(qū)。氣溶膠和氣體邊界條件來自于MOZART-4/GEOS-5(Model for Ozone and Related chemical Tracers-4/Goddard Earth Observing System-5)。

    表1 中國東部和中國區(qū)域排放源年排放總量

    地面氣溫觀測數(shù)據(jù)來自于中國氣象局(http:// www.cma.gov.cn/ [2013-10-10]),該數(shù)據(jù)基于中國的751個觀測站點,為了方便與模擬值對比分析,將觀測數(shù)據(jù)插值成與模式相同的網(wǎng)格,Xu et al. (2009)對插值方法進行了說明和驗證。氣溶膠光學(xué)厚度(Aerosol Optical Depth,AOD)和單次散射反照率(Single Scattering Albedo,SSA)觀測數(shù)據(jù)都來自于AERONET(Aerosol Robotic Network,http://aeronet.gsfc.nasa.gov/new_web/index.html[2013- 10-12])。為了更好的定量描述氣溶膠對天氣影響在不同區(qū)域的特征,本文劃分兩個重點區(qū)域:中國北方區(qū)域(29.8°~42.6°N,110.2°~120.3°E)和中國南方區(qū)域(22.3°~29.9°N,109.7°~120.2°E)進行分析研究。本文中日期和時間均為北京時。

    3 試驗設(shè)計

    表2給出了兩個試驗的試驗設(shè)計。我們選用模式平衡一天之后的結(jié)果來進行分析。試驗1中,選取不同季節(jié)且無明顯降水的5個時段進行模擬分析。每個時間段模擬均包括一個控制實驗(CTRL)和一個無氣溶膠敏感性試驗(NOAER),CTRL減NOAER即可得到氣溶膠對輻射和地面氣溫的影響。通過CTRL和NOAER試驗?zāi)M的地面氣溫分別與觀測值對比,分析氣溶膠對模式模擬地面氣溫的作用。WRF-Chem模式中氣溶膠通過直接和間接效應(yīng)影響天氣。直接效應(yīng)是將氣溶膠光學(xué)厚度計入大氣輻射傳輸模塊,影響對太陽輻射的收支計算(Fast et al., 2006);間接效應(yīng)是氣溶膠通過激活(活化)過程,形成云凝結(jié)核(Chapman et al., 2009)。因此,設(shè)所研究氣溶膠AOD為0,即關(guān)閉了其直接效應(yīng),關(guān)閉所研究氣溶膠的激活(活化)過程,即關(guān)閉了其間接效應(yīng)。試驗2中,分別模擬分析2006年8月23~25日期間氣溶膠及其部分組分(硫酸鹽、硝酸鹽和黑碳氣溶膠)的直接效應(yīng),間接效應(yīng)和總效應(yīng)對輻射和地面氣溫的影響,其中均分別包括一個控制試驗(CTRL)和兩個敏感性試驗(NODIR和NODIRIND)。CTRL減NODIRIND即為所研究氣溶膠的總效應(yīng),CTRL減NODIR即為直接效應(yīng),NODIR減NODIRIND即為間接效應(yīng)。

    表2 試驗設(shè)計

    4 模擬性能分析

    此部分選區(qū)多個時段中2006年8月23~25日時段,將模式CTRL模擬試驗的氣象參數(shù)(地面2 m溫度)、氣溶膠濃度和氣溶膠光學(xué)特性與觀測值進行了對比分析,從而對模式模擬性能進行了初步驗證。由于本文中觀測和模擬均采用的是地面2 m溫度,因此一致簡稱為地面氣溫。

    4.1 地面氣溫驗證

    圖1顯示了模擬和觀測的平均地面氣溫,模式能較好的再現(xiàn)地面氣溫分布。氣溫高值區(qū)主要分布在中國中東部,南方和四川盆地地區(qū),這些區(qū)域氣溫均達25°C以上。低值區(qū)主要分布在中國東北(10~20°C)以及高原地區(qū)(0~15°C)。模式在黑龍江和華南地區(qū)模擬氣溫偏低(最大差5°C)。分析模擬驅(qū)動氣象場NCEP數(shù)據(jù)顯示,這種模擬氣溫差異與驅(qū)動場誤差有關(guān)。

    圖1 2006年8月23~25日(a)模擬和(b)觀測的平均地面氣溫(單位:°C)

    圖2 2006年8月23~25日模式模擬的(a)PM2.5、(b)SO42-(硫酸鹽)、(c)NO3(硝酸鹽)和(d)BC(黑碳)平均近地面質(zhì)量濃度(單位:μg m?3)

    4.2 氣溶膠濃度驗證

    圖2顯示了模擬的細粒子氣溶膠PM2.5、硫酸鹽(SO42-),硝酸鹽(NO3-)和黑碳氣溶膠的地表濃度分布(北方和南方區(qū)域見圖2a紅框)??傮w上,氣溶膠濃度高值位于中國東部,顯示出華北平原,長江三角洲和四川盆地幾個高值中心。這與當?shù)亟?jīng)濟發(fā)展快所導(dǎo)致的人為氣溶膠及其前體物的排放總量高有關(guān)。華北地區(qū)PM2.5濃度最高值達70 μg m?3,在川渝、湖南及長江中下游及地區(qū)次之,濃度約為30~50 μg m?3。硫酸鹽氣溶膠的濃度高值區(qū)位于華北和四川東部地區(qū),約10~14 μg m?3,中國中東部其余地區(qū)約為4~10 μg m?3。硝酸鹽氣溶膠濃度高值位于河北南部、河南和湖北北部,達20~28 μg m?3,在湖南、山東和長江下游次之約8~20 μg m?3。黑碳氣溶膠濃度最低,主要位于河北南部和河南,約8~14 μg m?3,四川地區(qū)次之,約4~6 μg m?3。

    表3為模式模擬的PM2.5及其組分:硫酸鹽、硝酸鹽、黑碳、銨鹽(NH4+)和有機碳氣溶膠和臭氧(O3),二氧化硫和二氧化氮(NO2)氣體近地面濃度與文獻搜集所獲得的四個站點的觀測值進行對比。四個站點分別為:榆垡(39°51'N,116°30'E)、興?。?0°24'N,117°30'E)、西安(34°15'N, 108°56'E)、上海(31°13'N,121°28'E)。由表3可以看出,模式模擬的PM2.5及其組分和氣體濃度與觀測值較為吻合,但在西安站點模式模擬的PM2.5和硫酸鹽氣溶膠濃度略低于觀測值,分別低估 了40 μg m?3(33%)和18 μg m?3(56%),這是由于模式模擬中對西安城市地區(qū)排放量的低估所致。

    表3 2006年8月23~25日模擬和觀測近地面PM2.5及其組分濃度( SO42-,NO3-,BC,NH4+和OC,單位:μg m?3)和氣體濃度(O3,SO2和NO2,單位:10?9)

    注:榆垡、興隆、西安和上海的觀測數(shù)據(jù)分別來自于文獻:Matsui et al.(2009),Wu et al.(2011),Zhang et al.(2011)和Hou et al.(2011)。

    4.3 光學(xué)特性驗證

    氣溶膠光學(xué)特性的驗證中,模式模擬的AOD和SSA與四個站點的觀測值進行對比,四個站點分別為:榆垡(39°18'N,116°11'E)、太湖(31°25'N,120°12'E)、蘭州(35°56'N,104°08'E)和香港(22°18'N,114°10'E)。模式模擬的AOD與觀 測值的對比如圖3所示,從總體上,模式模擬的AOD在趨勢和量級上與觀測值有較好的吻合,但模擬值均存在不同程度的低估,尤其在蘭州和 香港這兩個城市站點,分別低估了0.15和0.4左右,這可能由于模式模擬的分辨率較高,模式 低估了這兩個城市站點的氣溶膠濃度所造成 的。表4給出了模式模擬的SSA與觀測值的對比,由表4可以看出,模式模擬的SSA略低于觀測值(最大在榆垡站低估了0.21),但總體上吻合較好。

    圖3 2006年8月23~25日模擬(紅線)和觀測(黑點)的氣溶膠光學(xué)厚度(AOD,波長:500 nm)

    表4 2006年8月23~25日模擬和觀測的單次散射反照率(SSA,波長:500 nm)

    5 模擬結(jié)果

    5.1 氣溶膠對模式模擬地面氣溫的影響

    試驗1中五個不同時間段模式模擬的區(qū)域平均PM2.5近地面質(zhì)量濃度和對地面氣溫的影響如圖4所示,在此我們選用北方區(qū)域和南方區(qū)域進行對照分析。由圖4中可見,PM2.5近地面濃度在南北區(qū)域均有明顯的日變化特征。氣溶膠濃度在夜間由于大氣對流運動減弱逐漸累積,早上06:00~08:00達到最高,隨著日出后太陽輻射逐漸增強,空氣對流活動加強,有利于氣溶膠擴散,氣溶膠濃度逐漸降低,下午17:00左右達到最低。五個時段中,北方區(qū)域的氣溶膠濃度明顯高于南方區(qū)域,PM2.5近地面濃度的高值在北方區(qū)域約為40~80 μg m?3,在南方區(qū)域約為30~47 μg m?3,最高值(80 μg m?3)出現(xiàn)在2009年冬季的模擬中(圖4d),這是由于冬季大氣對流運動較弱造成的。

    圖4 模式模擬的北方(實線)和南方(虛線)區(qū)域平均PM2.5近地面質(zhì)量濃度(右坐標軸,黑線,單位:μg m?3)和對地面氣溫的影響(左坐標軸,紅線,單位:K)

    在五個時段的模擬中,氣溶膠均對地面氣溫有明顯的降溫效果(個別時間的輕微增溫是由黑碳氣溶膠的加熱作用所致),并且存在明顯的日變化。由于夜間沒有太陽輻射,氣溶膠對夜間氣溫影響較小,其影響主要集中在白天。氣溶膠對地面氣溫影響的第一個高值出現(xiàn)在早上08:00左右,雖然此時太陽輻射并未達到最高,但由于對應(yīng)氣溶膠濃度日變化最高值,因此降溫較為明顯。此后,隨著氣溶膠濃度的降低,對地面氣溫的影響逐漸減小,直至正午12:00左右,隨著太陽輻射達到最高,氣溶膠降溫作用再次增強,并在16:00~18:00左右達到另一個峰值。五個時段的模擬結(jié)果中,氣溶膠最高可使地面氣溫降低約0.2~1 K,降溫效果最強為2009年春季的模擬,氣溶膠使北方區(qū)域平均降溫最高將近1 K(圖4e),雖然此時段的氣溶膠近地面濃度不是最高,但是可能由于此時段中間接效應(yīng)起到了很強的作用,導(dǎo)致降溫劇烈。本文模擬得到的氣溶膠的冷卻作用比Zhang(2008)應(yīng)用WRF-Chem模式了德克薩斯州地區(qū)2000年8月的5天中,氣溶膠使近地面氣溫降低作用(最高0.18 K)更大,這可能的原因在于中國地區(qū)氣溶膠濃度遠高于美國地區(qū)(PM2.5質(zhì)量濃度僅為20 μg m?3),因此其冷卻效果更加顯著。

    表5統(tǒng)計了北方和南方區(qū)域模式模擬的有氣溶膠(CTRL)和無氣溶膠(NOAER)的地面氣溫與觀測值的平均絕對誤差和均方根誤差(標準誤差)。由表5可見,在有氣溶膠的地面氣溫的平均絕對誤差和均方根誤差均小于無氣溶膠,因此在有氣溶膠的情況下模擬的地面氣溫更接近于實際觀測,氣溶膠的參與改善了模式對地面氣溫的模擬。

    表5 北方(南方)區(qū)域模式模擬的有氣溶膠(CTRL)和無氣溶膠(NOAER)的地面氣溫與觀測值的平均絕對誤差(單位:°C)和均方根誤差(單位:°C)

    注:括號為南方區(qū)域。平均絕對誤差:,均方根誤差:,為網(wǎng)格數(shù)量,為網(wǎng)格,f為網(wǎng)格的模擬值,r為網(wǎng)格的觀測值(Talor,2001)。

    5.2 氣溶膠及其部分組分的輻射強迫和對地面氣溫的影響

    依據(jù)試驗2的模擬結(jié)果,此部分定量分析2006年8月23~25日期間氣溶膠及其部分組分(SO42-、NO3-和BC)的總效應(yīng)、直接效應(yīng)和間接效應(yīng)分別對到達地面短波輻射和地面氣溫的影響。

    圖5給出了北方和南方區(qū)域平均PM2.5及其部分組分(SO42-、NO3-和BC)近地面質(zhì)量濃度日變化特征。整體上看,北方區(qū)域PM2.5,SO42-,NO3-和BC的近地面質(zhì)量濃度均大于南方區(qū)域。但北方區(qū)域和南方區(qū)域PM2.5濃度均呈現(xiàn)出明顯的日變化特征,在早上07時左右達到最高值,北方區(qū)域和南方區(qū)域PM2.5濃度分別約為51 μg m?3和25 μg m?3;在傍晚17時左右濃度達到最低,分別約為23 μg m?3和16 μg m?3。硫酸鹽氣溶膠主要通過SO2氧化生成,但其濃度沒有明顯的日變化特征,北方區(qū)域和南方區(qū)域濃度最高值均分別約為8 μg m?3和5 μg m?3。而硝酸鹽氣溶膠濃度具有明顯的夜間高,白天低的趨勢,北方區(qū)域和南方區(qū)域硝酸鹽氣溶膠濃度均在早上07時左右達到最高值,分別約為16 μg m?3和8 μg m?3,此后呈不斷下降趨勢,在傍晚18時左右達到最低值,分別約為5 μg m?3和4 μg m?3。本文中模擬得到的硝酸鹽氣溶膠濃度日變化特征與Morin et al.(2011)的研究結(jié)果相一致,HNO3氣體與硝酸鹽氣溶膠間的氣粒轉(zhuǎn)化對溫度及其敏感,當溫度較低時主要以硝酸鹽存在,而溫度較高時主要以氣態(tài)HNO3形式存在,因此硝酸氣溶膠濃度呈現(xiàn)出白天低夜晚高的日變化特征。黑碳氣溶膠的濃度亦呈現(xiàn)出在早上7時出現(xiàn)峰值隨后逐漸減小至下午16~18時達到低值的日變化特征,這一日變化特征在北方地區(qū)尤為顯著。這是由于夜間對流活動較弱,氣溶膠易于累積,在早間07時左右黑碳氣溶膠濃度達到最高值,北方和南方區(qū)域黑碳氣溶膠濃度分別約為:6 μg m?3和2 μg m?3,之后隨著太陽輻射的增加,對流運動的加強,濃度逐漸減小,在下午達到最低值,北方和南方區(qū)域黑碳氣溶膠濃度分別約為:2 μg m?3和1 μg m?3。本文模擬黑碳氣溶膠濃度的日變化特點與李楊等(2005)和蔡子穎等(2011)結(jié)果較為一致。

    圖6給出了氣溶膠及其部分組分(SO42-、NO3-和BC)的地面短波輻射強迫以及對地面氣溫的影響的區(qū)域平均值??傮w上,氣溶膠及硫酸鹽和硝酸鹽氣溶膠能夠吸收和散射太陽輻射,并能夠形成云凝結(jié)核,增加云的光學(xué)厚度和云層反射率,因此氣溶膠的直接與間接效應(yīng)均降低了到達地面的短波輻射,導(dǎo)致地面氣溫降低。從區(qū)域上看,北方區(qū)域平均氣溶膠直接效應(yīng)所造成的負的短波輻射強迫作用(-11.3 W m?2)要高于南方區(qū)域(-5.8 W m?2)(圖6a、b),引起北方和南方地區(qū)地面氣溫分別降低了0.074 K和0.039 K(圖6c、d),這主要是由于北方區(qū)域氣溶膠濃度高于南方區(qū)域所致(圖5)。但是由于南方區(qū)域相對濕度高(80%~90%,圖略),能夠促進氣溶膠的激活,從而有利于云凝結(jié)核的形成和增長(Jia and Guo, 2012)。南方區(qū)域氣溶膠間接效應(yīng)所產(chǎn)的的短波輻射強迫高于北方區(qū)域,分別為-14.4 W m?2和-12.4 W m?2,引起的地面氣溫的改變分別為-0.094 K和-0.035 K。

    圖5 2006年8月23~25日模式模擬的(a)北方和(b)南方區(qū)域平均PM2.5(左坐標軸)及其部分組分(SO42-、NO3-和BC,右坐標軸)近地面質(zhì)量濃度(單位:μg m?3)的日變化

    圖6 2006年8月23~25日北方(左列)和南方(右列)區(qū)域平均氣溶膠(AER)及其部分組分(SO42-、NO3-和BC)(a、b)地面短波輻射強迫以及對(c、d)地面氣溫的影響

    對于硫酸鹽氣溶膠而言,其在南方區(qū)域的直接效應(yīng)和間接效應(yīng)的輻射強迫和對地面氣溫的影響略高于北方區(qū)域。這主要的原因在于南方區(qū)域的具有更強的太陽輻射和濕度。硫酸鹽氣溶膠在北方和南方區(qū)域的短波輻射強迫分別為-7.0 W m?2和 -10.5 W m?2(直接和間接效應(yīng)之和),對地面氣溫的影響分別為-0.062 K和-0.074 K。其中硫酸鹽氣溶膠的直接效應(yīng)和間接效應(yīng)的作用相當。相比于Giorgi et al.(2002)應(yīng)用RegCM模式模擬1993~1997年間中國四川盆地人為硫酸鹽氣溶膠使地面氣溫降低0.1~0.7 K,本文模擬數(shù)值略小,可能是由于本文計算的輻射強迫和地面氣溫變化為區(qū)域平均值,包括了部分氣溶膠低值區(qū)域,而四川盆地是硫酸鹽氣溶膠高值區(qū)(圖2b),因而其對地表氣溫的影響更強。

    硝酸鹽氣溶膠的短波輻射強迫以及對地面氣溫的影響略小于硫酸鹽氣溶膠,北方和南方區(qū)域平均地面短波輻射分別減少4.1 W m?2和6.9 W m?2,因而使得地面氣溫亦分別減少0.021 K和0.039 K,略低于Li et al.(2009)利用模式研究得到的2003年7月中國地區(qū)人為硝酸鹽氣溶膠對地面氣溫的影響(-0.09 K)。從區(qū)域來看,北方區(qū)域的硝酸鹽氣溶膠的直接效應(yīng)的作用小于間接效應(yīng),換言之,即硝酸鹽氣溶膠引起的北方區(qū)域地面氣溫的降低主要是由硝酸鹽氣溶膠的間接效應(yīng)所致。而南方區(qū)域則是相反,即由硝酸鹽的直接效應(yīng)起主要作用。這可能是由于雖然南方區(qū)域硝酸鹽氣溶膠濃度小 于北方區(qū)域,但是南方區(qū)域相對濕度大(AOD 高),且有更強的太陽輻射,因此其直接效應(yīng)更加顯著。北方和南方區(qū)域均有一定的云量(圖略),雖然南方區(qū)域相對濕度很高,有利于硝酸鹽氣溶膠的激活形成云凝結(jié)核,但是其云量的增加(0.004)沒有北方區(qū)域(0.03)多,其間接效應(yīng)相對較弱。

    與硫酸鹽和硝酸鹽氣溶膠(散射型氣溶膠)不同,黑碳氣溶膠作為重要的吸收性氣溶膠,能夠強烈地吸收太陽輻射,從而加熱大氣,因此其直接效應(yīng)雖然同樣減少了到達地表的短波輻射,但是地面氣溫卻增加了;黑碳氣溶膠間接輻射強迫很小并且有正有負,具有很強的不確定性,這主要是由于云的反饋過程比較復(fù)雜造成的(Ackerman et al., 2000; Kristjánsson, 2002; Highwood and Kinnersley, 2006; Bond et al., 2013)。由圖6可見,黑碳氣溶膠使得北方和南方區(qū)域平均到達地表的太陽短波輻射分別減少了6.5 W m?2和5.8 W m?2,而地表氣溫則分別增加了0.053 K和0.017 K。雖然在南方區(qū)域黑碳氣溶膠濃度低,但黑碳氣溶膠的消光系數(shù)(12.5 m2g?1,波長:550 nm)要強于硫酸鹽氣溶膠(4.18 m2g?1,波長:550 nm),并且其吸收性很強(Liao et al., 2004),所以黑碳氣溶膠在南方的輻射強迫的要高于硫酸鹽氣溶膠。

    此模擬的地面氣溫變化與Wang(2004)應(yīng)用CCM3模式研究發(fā)現(xiàn)黑碳氣溶膠20年全球平均直接效應(yīng)增加地面氣溫0.09±0.07 K較為相近。

    圖7進一步展示了兩個研究區(qū)域平均氣溶膠及其部分組分(SO42-、NO3-和BC)總效應(yīng)、直接效應(yīng)和間接效應(yīng)的地面短波輻射強迫及其對地面氣溫的影響的日變化結(jié)果。整體上看,氣溶膠對到達地面短波輻射的影響與對地面氣溫的影響有很好的對應(yīng)關(guān)系。

    由圖7a、d、g、j可見,氣溶膠總效應(yīng)導(dǎo)致的地面短波輻射強迫分別在上午08~10時左右和下午16~18時左右呈現(xiàn)出最大值,約為-43~-39 W m?2。氣溶膠的直接和間接效應(yīng)均對上午08~10時左右輻射強迫峰值有貢獻。上午08~10時左右氣溶膠濃度高(圖5),并且此時太陽天頂角大,因此造成了很強的直接輻射強迫。如Liao and Seinfeld(1998)指出,散射性氣溶膠(硫酸鹽)輻射強迫隨太陽天頂角增加而增加。由于此時南方區(qū)域的云量的增加達到了峰值(圖未給),所以此時的間接輻射強迫也非常顯著。下午16~18時左右雖然太陽天頂角大,但是氣溶膠濃度低,因此直接效應(yīng)較弱,此時輻射強迫峰值主要受間接效應(yīng)影響,在下午16~18時左右北方區(qū)域云量的增加達到峰值,南方區(qū)域(峰值相對較小)云量在此階段也在不斷上升。

    受到輻射強迫高值影響,地面氣溫在這2個時刻下降幅度最大,約為-0.27~-0.21 K。其中,硫酸鹽氣溶膠總效應(yīng)對短波輻射的影響在中午12時和下午16時左右分別達到最高約為-28(南方區(qū)域)~-25 W m?2(北方區(qū)域),對地面氣溫的影響最大為在-0.13 K(南方區(qū)域)??傮w而言,硝酸鹽氣溶膠對短波輻射和氣溫的影響比硫酸鹽氣溶膠小,這是由硝酸鹽氣溶膠濃度白天低夜間高的日變化特征所決定的。硝酸鹽氣溶膠對短波輻射和地面氣溫的影響在南方區(qū)域更加顯著,最大分別為-33 W m?2和-0.09 K(14:00左右)。受太陽輻射強度的影響,黑碳氣溶膠對短波輻射的影響最大值出現(xiàn)在南方區(qū)域,約為-38 W m?2,高于北方區(qū)域的-17 W m?2,其升溫作用明顯,北方和南方區(qū)域均在中午14時左右達到最高,分別為0.21 K和0.13 K。

    北方區(qū)域氣溶膠直接效應(yīng)導(dǎo)致的地面短波輻射強迫在上午09時左右達到最高值:-32 W m?2,從而使地面氣溫降低0.2 K,而隨著黑碳氣溶膠強烈的加熱作用,15時前后氣溫上升,最高升溫0.08 K。南方區(qū)域氣溶膠直接效應(yīng)導(dǎo)致的輻射強迫和降溫最大為-26 W m?2和-0.08 K,同樣受黑碳氣溶膠增溫的影響,午后出現(xiàn)增溫(最高0.05 K)。硫酸鹽和硝酸鹽氣溶膠的直接效應(yīng)對輻射和氣溫影響的日變化比較類似,均在中午時段顯現(xiàn)高值。其中,南方區(qū)域硝酸鹽的直接效應(yīng)的輻射強迫最為明顯,在14時達到最高值-45 W m?2,相對應(yīng)的使氣溫降低0.12 K。黑碳氣溶膠的直接效應(yīng)對短波輻射的影響十分顯著,在北方和南方區(qū)域短波輻射強迫最高分別約為-18 W m?2和-27 W m?2,地面氣溫在14時左右升高最為明顯,分別為0.15 K和0.04 K(圖7b、e、h、k)。

    圖7 2006年8月23~25日模式模擬的北方和南方區(qū)域平均氣溶膠(AER)及其部分組分(SO42-、NO3-和BC)總效應(yīng)、直接效應(yīng)和間接效應(yīng)的地面短波輻射強迫(單位:W m?2)及其對地面氣溫的影響的日變化(單位:K)

    北方區(qū)域氣溶膠間接效應(yīng)導(dǎo)致的地面短波輻射強迫在午后達到最高值-33 W m?2,其貢獻主要來自于硫酸鹽和硝酸鹽氣溶膠,與之相對應(yīng),使地面氣溫降低0.19 K。而在南方區(qū)域,氣溶膠間接效應(yīng)對短波輻射和地面氣溫的影響均表現(xiàn)為雙峰型,高值分別在上午10~11時左右和下午16~17時左右,約為:-33 W m?2和-0.23 K。硫酸鹽和硝酸鹽氣溶膠的間接效應(yīng)對輻射和地表氣溫的影響并沒有很明顯的日變化特征,但總體上是負強迫和降溫的效果。硫酸鹽氣溶膠間接效應(yīng)在北方區(qū)域的短波輻射強迫約-15~3 W m?2,對地面氣溫的影響為-0.07~0.02 K;在南方區(qū)域的短波輻射強迫約-13~1 W/m2,對地面氣溫的影響為-0.08~0 K。硝酸鹽氣溶膠間接效應(yīng)的短波輻射強迫和對地表氣溫的影響在北方區(qū)域更加顯著,最高分別為-16 W m?2和-0.06 K。相比于直接效應(yīng)的影響,黑碳氣溶膠間接效應(yīng)對地面短波輻射和氣溫的影響要小很多,并且其間接效應(yīng)具有很強的不確定性,有待進一步研究(圖7c、f、i、l)。

    5 結(jié)論

    本文應(yīng)用WRF-Chem模式研究中國東部地區(qū)氣溶膠及其部分組分(SO42-、NO3-和BC)在天氣尺度下的輻射強迫和對地面氣溫的影響,主要結(jié)論如下:

    模式能夠較好地模擬中國東部地區(qū)的地面氣溫,氣溶膠及其部分組分(SO42-、NO3-和BC)的分布。通過不同季節(jié)且無明顯降水的5個時間段的模擬研究氣溶膠對地面氣溫的影響,結(jié)果顯示,北方區(qū)域平均PM2.5近地面濃度(40~80 μg m?3)高于南方(30~47 μg m?3),并且均體現(xiàn)出白天低夜間高的日變化特征。氣溶膠對地面氣溫有明顯的降溫效果,在早上08時和下午17時左右最為顯著,最高可降低約0.2~1 K。模擬的地表氣溫與觀測數(shù)據(jù)的對比顯示,有氣溶膠的情況下模擬的地面氣溫更接近于實際觀測。

    本文定量分析了2006年8月23~25日期間氣溶膠及其部分組分(SO42-、NO3-和BC)的總效應(yīng),直接效應(yīng)和間接效應(yīng)分別對到達地面短波輻射和地面氣溫的影響。北方區(qū)域平均氣溶膠直接效應(yīng)所造成的負的短波輻射強迫作用(-11.3 W m?2)要高于南方區(qū)域(-5.8 W m?2),引起北方和南方地區(qū)地面氣溫分別降低了0.074 K和0.039 K,其中最高值分別約-32 W m?2和0.2 K(上午09時)。南方區(qū)域平均氣溶膠間接效應(yīng)所產(chǎn)的的短波輻射強迫高于北方區(qū)域,分別為-14.4 W m?2和-12.4 W m?2,引起的地面氣溫的改變分別為-0.094 K和-0.035 K,其中最高值分別約-33 W m?2和-0.23 K。

    對氣溶膠組分(SO42-、NO3-和BC)的研究中顯示,硫酸鹽和硝酸鹽氣溶膠做為散射型氣溶膠能夠散射和反射太陽輻射,并且形成云凝結(jié)核,增加云量,減少到達地面的太陽輻射,使地面氣溫降低。硫酸鹽氣溶膠的直接效應(yīng)和間接效應(yīng)的作用相當,其總效應(yīng)在北方和南方區(qū)域平均的短波輻射強迫分別為-7.0 W m?2和-10.5 W m?2,對地面氣溫的影響分別為-0.062 K和-0.074 K。硝酸鹽氣溶膠的短波輻射強迫以及對地面氣溫的影響略小于硫酸鹽氣溶膠,北方和南方區(qū)域平均的地面短波輻射分別減少4.1 W m?2和6.9 W m?2,因而使得地面氣溫亦分別減少0.021 K和0.039 K。硫酸鹽和硝酸鹽氣溶膠的直接效應(yīng)對輻射和氣溫影響的日變化比較類似,均在中午時段顯現(xiàn)高值,但它們的間接效應(yīng)并沒有很明顯的日變化特征,但總體上也是負強迫和降溫的效果。其中,南方區(qū)域硝酸鹽的直接效應(yīng)的輻射強迫最為明顯,在14時達到最高值-45 W m?2,相對應(yīng)的使氣溫降低0.12 K。黑碳氣溶膠使得北方和南方區(qū)域平均到達地表的太陽短波輻射分別減少了6.5 W m?2和5.8 W m?2,而地表氣溫則分別增加了0.053 K和0.017 K。黑碳氣溶膠能夠吸收太陽輻射,并放出紅外輻射,從而加熱大氣,相比于間接效應(yīng),黑碳氣溶膠的直接效應(yīng)對短波輻射和氣溫的影響更加顯著,在北方和南方區(qū)域短波輻射強迫最高分別約為-18 W m?2和-27 W m?2,增溫作用在14時左右最為明顯,分別為0.15 K和0.04 K。

    (References:)

    Ackerman A S, Toon O B, Stevens D E, et al. 2000. Reduction of tropical cloudiness by soot [J]. Science, 288 (5468): 1042?1047.

    Bond T, Doherty S, Fahey D, et al. 2013. Bounding the role of black carbon in the climate system: A scientific assessment [J]. J. Geophys. Res., 118 (11): 5380?5552.

    蔡子穎, 韓素芹, 黃鶴, 等. 2011. 天津夏季黑碳氣溶膠及其吸收特性的觀測研究[J]. 中國環(huán)境科學(xué), 31 (5): 719?723. Cai Ziying, Han Suqin, Huang He, et al. 2011. Observational study on black carbon aerosols and their absorption properties in summer in Tianjin [J]. China Environmental Science (in Chinese), 31 (5): 719?723.

    Chapman E G, Gustafson W I, Easter R C, et al. 2009. Coupling aerosol-cloud-radiative processes in the WRF-Chem model: Investigating the radiative impact of elevated point sources [J]. Atmos. Chem. Phys., 9 (3): 945?964.

    Chen F, Dudhia J. 2001. Coupling an advanced land surface-hydrology model with the Penn State-NCAR MM5 modeling system. Part I: Model implementation and sensitivity [J]. Mon. Wea. Rev., 129 (4): 569?585.

    Ek M B, Mitchell K E, Lin Y, et al. 2003. Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model [J]. J. Geophys. Res., 108 (D22), doi:10.1029/2002JD003296.

    Fahey K M, Pandis S N. 2001. Optimizing model performance: Variable size resolution in cloud chemistry modeling [J]. Atmos. Environ., 35 (26): 4471?4478.

    Fast J D, Gustafson W I Jr, Easter R C, et al. 2006. Evolution of ozone, particulates, and aerosol direct radiative forcing in the vicinity of Houston using a fully coupled meteorology-chemistry-aerosol model [J]. J. Geophys. Res., 111 (D21), doi:10.1029/2005JD006721.

    Giorgi F, Bi X Q, Qian Y. 2002. Direct radiative forcing and regional climatic effects of anthropogenic aerosols over East Asia: A regional coupled climate-chemistry/aerosol model study [J]. J. Geophys. Res., 107 (D20), doi:10.1029/2001JD001066.

    Grell G A, Dévényi D. 2002. A generalized approach to parameterizing convection combining ensemble and data assimilation techniques [J]. Geophys. Res. Lett., 29 (14): 38-1–38-4.

    Grell G A, Peckham S E, Schmitz R, et al. 2005. Fully coupled “online” chemistry within the WRF model [J]. Atmos. Environ., 39 (37): 6957?6975.

    Guenther A, Zimmerman P, Wildermuth M. 1994. Natural volatile organic compound emission rate estimates for US woodland landscapes [J]. Atmos. Environ., 28 (6): 1197–1210.

    Guenther A B, Zimmerman P R, Harley P C, et al. 1993. Isoprene and monoterpene emission rate variability-model evaluations and sensitivity analyses [J]. J. Geophys. Res.-Atmospheres, 98 (D7): 12609–12617.

    Gustafson W I, Chapman E G, Ghan S J, et al. 2007. Impact on modeled cloud characteristics due to simplified treatment of uniform cloud condensation nuclei during NEAQS 2004 [J]. Geophys. Res. Lett., 34 (19), doi:10.1029/2007GL030021.

    Highwood E J, Kinnersley R P. 2006. When smoke gets in our eyes: The multiple impacts of atmospheric black carbon on climate, air quality and health [J]. Environ. Inter., 32 (4): 560?566.

    Hong S Y, Noh Y, Dudhia J. 2006. A new vertical diffusion package with an explicit treatment of entrainment processes [J]. Mon. Wea. Rev., 134 (9): 2318?2341.

    Hou B, Zhuang G S, Zhang R, et al. 2011. The implication of carbonaceous aerosol to the formation of haze: Revealed from the characteristics and sources of OC/EC over a mega-city in China [J]. J. Hazard. Mater., 190 (1?3): 529?536.

    IPCC. 2007. The Physical Science Basis of Climate Change: Changes in Atmospheric Constituents and in Radiative Forcing [M]. New York: Cambridge University Press.

    Jacobson M Z. 2001. Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols [J]. Nature, 409 (6821): 695?697.

    Jia X C, Guo X L. 2012. Impacts of secondary aerosols on a persistent fog event in northern China [J]. Atmos. Oceanic Sci. Lett., 5 (5): 401?407.

    Kristjánsson J E. 2002. Studies of the aerosol indirect effect from sulfate and black carbon aerosols [J]. J. Geophys. Res., 107 (D15), doi:10.1029/ 2001JD000887.

    Li S, Wang T J, Zhuang B L, et al. 2009. Indirect radiative forcing and climatic effect of the anthropogenic nitrate aerosol on regional climate of China [J]. Adv. Atmos. Sci., 26 (3): 543?552.

    李楊, 曹軍驥, 張小曳, 等. 2005. 2003 年秋季西安大氣中黑碳氣溶膠的演化特征及其來源解析[J]. 氣候與環(huán)境研究, 10 (2): 229?237. Li Yang, Cao Junji, Zhang Xiaoye, et al. 2005. The variability and source apportionment of black carbon aerosol in Xi’an atmosphere during the autumn of 2003 [J]. Climatic and Environmental Research (in Chinese), 10(2): 229?237.

    Liao H, Seinfeld J H. 1998. Effect of clouds on direct aerosol radiative forcing of climate [J]. J. Geophys. Res., 103 (D4): 3781?3788.

    Liao H, Seinfeld J H, Adams P J, et al. 2004. Global radiative forcing of coupled tropospheric ozone and aerosols in a unified general circulation model [J]. J. Geophys. Res., 109 (D24), doi:10.1029/2003JD004456.

    Lin Y L, Farley R D, Orville H D. 1983. Bulk parameterization of the snow field in a cloud model [J]. J. Climate Appl. Meteor., 22 (6): 1065?1092.

    Matsui H, Koike M, Kondo Y, et al. 2009. Spatial and temporal variations of aerosols around Beijing in summer 2006: Model evaluation and source apportionment [J]. J. Geophys. Res., 114, doi:10.1029/2008JD010906.

    Mlawer E J, Taubman S J, Brown P D, et al. 1997. Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-model for the longwave [J]. J. Geophys. Res., 102 (D14): 16663?16682.

    Morin S, Sander R, Savarino J. 2011. Simulation of the diurnal variations of the oxygen isotope anomaly (△17O) of reactive atmospheric species [J]. Atmos. Chem. Phys., 11 (8): 3653?3671.

    Pérez C, Nickovic S, Pejanovic G, et al. 2006. Interactive dust-radiation modeling: A step to improve weather forecasts [J]. J. Geophys. Res., 111 (D16), doi:10.1029/2005JD006717.

    Rodwell M J, Jung T. 2008. Understanding the local and global impacts of model physics changes: An aerosol example [J]. Quart. J. Roy. Meteor. Soc., 134 (635): 1479?1497.

    Takemura T, Nozawa T, Emori S, et al. 2005. Simulation of climate response to aerosol direct and indirect effects with aerosol transport- radiation model [J]. J. Geophys. Res., 110 (D2), doi:10.1029/ 2004JD005029.

    Taylor K E. 2001. Summarizing multiple aspects of model performance in a single diagram. [J]. J. Geophys. Res., 106 (D7): 7183?7192.

    Wang C E. 2004. A modeling study on the climate impacts of black carbon aerosols [J]. J. Geophys. Res., 109 (D3), doi:10.1029/2003JD004084.

    Wang Y, Che H Z, Ma J Z, et al. 2009. Aerosol radiative forcing under clear, hazy, foggy, and dusty weather conditions over Beijing, China [J]. Geophys. Res. Lett., 36, doi:10.1029/2009GL037181.

    Wild O, Zhu X, Prather M J. 2000. Fast-j: Accurate simulation of in- and below-cloud photolysis in tropospheric chemical models [J]. J. Atmos. Chem., 37 (3): 245?282.

    Wu Q Z, Wang Z F, Gbaguidi A, et al. 2011. A numerical study of contributions to air pollution in Beijing during CAREBeijing-2006 [J]. Atmos. Chem. Phys., 11 (12): 5997?6011.

    Xu Y, Gao X J, Shen Y, et al. 2009. A daily temperature dataset over China and its application in validating a RCM simulation [J]. Adv. Atmos. Sci., 26 (4): 763?772.

    Zaveri R A, Peters L K. 1999. A new lumped structure photochemical mechanism for large-scale applications [J]. J. Geophys. Res., 104 (D23): 30387?30415.

    Zaveri R A, Easter R C, Fast J D, et al. 2008. Model for simulating aerosol interactions and chemistry (MOSAIC) [J]. J. Geophys. Res., 113 (D13), doi:10.1029/2007JD008782.

    Zhang T, Cao J J, Tie X X, et al. 2011. Water-soluble ions in atmospheric aerosols measured in Xi'an, China: Seasonal variations and sources [J]. Atmos. Res., 102 (1?2): 110?119.

    Zhang Y. 2008. Online-coupled meteorology and chemistry models: History, current status, and outlook [J]. Atmos. Chem. Phys., 8 (11): 2895?2932.

    Zhang Y, Liao H, Zhu K F, et al. 2009. Role of black carbon-induced changes in snow albedo in predictions of temperature and precipitation during a snowstorm [J]. Atmos. Oceanic Sci. Lett., 2 (4): 230?236.

    Zhang Y, Wen X Y, Jang C J. 2010a. Simulating chemistry-aerosol-cloud- radiation-climate feedbacks over the continental U. S. using the online-coupled Weather Research Forecasting Model with chemistry (WRF/Chem) [J]. Atmos. Environ., 44 (29): 3568?3582.

    Zhang Y, Pan Y, Wang K, et al. 2010b. WRF/Chem-MADRID: Incorporation of an aerosol module into WRF/Chem and its initial application to the TexAQS2000 episode [J]. J. Geophys. Res., 115 (D18), doi:10.1029/2009JD013443.

    廖禮,漏嗣佳,符瑜,等. 2015. 中國東部氣溶膠在天氣尺度上的輻射強迫和對地面氣溫的影響[J]. 大氣科學(xué), 39 (1): 68?82, doi: 10.3878/j.issn. 1006-9895.1402.13302. Liao Li, Lou Sijia, Fu Yu, et al. 2015. Radiative forcing of aerosols and its impact on surface air temperature on the synoptic scale in eastern China[J]. Chinese Journal of Atmospheric Sciences (in Chinese), 39 (1): 68?82.

    Radiative Forcing of Aerosols and Its Impact on Surface Air Temperature on Synoptic Scale in Eastern China

    LIAO Li1, 2, LOU Sijia1, 2, FU Yu3, CHANG Wenyuan1, and LIAO Hong1

    1(),,100029;2,100049;3(),,100029

    Using the WRF-Chem (Weather Research and Forecasting—Chemistry) model, we investigate aerosol radiative forcing and its impacts on surface air temperature in Eastern China on the synoptic scale. Simulations of five selected periods without precipitation (August 23–25, 2006; November 10–12, 2008; December 16–18, 2008; January 15–17, 2009; and April 27–29, 2009) indicate that aerosol concentrations are generally lower in the daytime and higher at night. Simulated surface–layer concentrations of PM2.5averaged over Northern China (29.8°–42.6°N, 110.2°–120.3°E,) are 40–80 μg m?3, higher than the simulated values of 30–47 μg m?3 in Southern China (22.3°–29.9°N,109.7°–120.2°E). Aerosols exerted a large cooling effect on surface air temperature (SAT, 2 m above ground level), and the effect is the strongest at 08:00 and 17:00 with decreases in temperature by 0.2–1 K. We also find that including aerosols in simulations improves the simulation of SAT. Based on the simulation of a weather event during August 23 to 25 in 2006, we quantitatively analyze the total (direct + indirect), direct, and indirect radiative effects of aerosols on SAT. Shortwave radiative forcing at the surface induced by direct effects of aerosols is stronger in Northern China than that in Southern China, with average forcings of -11.3 W m?2and -5.8 W m?2, respectively, and corresponding coolings of SAT by 0.074 K and 0.039 K. However, shortwave radiative forcing induced by indirect effects of aerosols is weaker in Northern China than that in Southern China, with average forcings of -12.4 W/m2and -14.4 W m?2, respectively, and corresponding coolings of SAT by 0.035 K and 0.094 K. The direct and indirect effects of sulfate are of similar magnitude. The total effect of sulfate on shortwave radiative forcing in Northern and Southern China is -7. 0 W m?2and -10.5 W m?2, respectively, causing a cooling of SAT by 0.062 K and 0.074 K, while the effect of nitrate is relative weak. Due to black carbon aerosols, the solar shortwave radiation reaching the surface is decreased by 6.5 W m?2and 5.8 W m?2and the SAT increases by 0.053 K and 0.017 K averaged over Northern and Southern China, respectively. The direct effect of black carbon on shortwave radiation is much larger than the indirect effect.

    Eastern China, Aerosols, Radiative forcing, Surface air temperature

    1006–9895(2015)01–0068–15

    P422

    A

    10.3878/j.issn.1006-9895.1402.13302

    2013–11–06;網(wǎng)絡(luò)預(yù)出版日期2014–02–17

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)項目2014CB441200

    廖禮,男,1984年出生,博士研究生,從事氣溶膠數(shù)值模擬研究。E-mail: liaoli@mail.iap.ac.cn

    廖宏,E-mail: hongliao@mail.iap.ac.cn

    猜你喜歡
    效應(yīng)區(qū)域
    鈾對大型溞的急性毒性效應(yīng)
    永久基本農(nóng)田集中區(qū)域“禁廢”
    分割區(qū)域
    懶馬效應(yīng)
    場景效應(yīng)
    應(yīng)變效應(yīng)及其應(yīng)用
    關(guān)于四色猜想
    分區(qū)域
    基于嚴重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    偶像效應(yīng)
    插逼视频在线观看| 91国产中文字幕| 97在线人人人人妻| 亚洲欧美一区二区三区国产| av天堂久久9| av免费在线看不卡| 两个人的视频大全免费| 欧美精品一区二区大全| 精品熟女少妇av免费看| 国产熟女欧美一区二区| 大香蕉久久网| 日韩熟女老妇一区二区性免费视频| 80岁老熟妇乱子伦牲交| 午夜激情福利司机影院| 狂野欧美激情性bbbbbb| 免费日韩欧美在线观看| 国产精品熟女久久久久浪| 老司机影院成人| 国产精品熟女久久久久浪| 中文乱码字字幕精品一区二区三区| 少妇猛男粗大的猛烈进出视频| 中文字幕av电影在线播放| 黄色毛片三级朝国网站| 国产免费视频播放在线视频| 一区二区三区乱码不卡18| 大香蕉97超碰在线| 日本欧美国产在线视频| 人妻少妇偷人精品九色| 日韩一区二区三区影片| 18禁动态无遮挡网站| 91精品国产九色| 精品亚洲乱码少妇综合久久| 亚洲精品乱码久久久v下载方式| 乱码一卡2卡4卡精品| 久久久国产欧美日韩av| 精品一区二区三卡| 国产一区有黄有色的免费视频| 成年女人在线观看亚洲视频| 91精品一卡2卡3卡4卡| 2021少妇久久久久久久久久久| 多毛熟女@视频| 久久久精品94久久精品| 日日爽夜夜爽网站| 国产在视频线精品| 国产成人av激情在线播放 | 午夜精品国产一区二区电影| 高清不卡的av网站| 69精品国产乱码久久久| 黄色配什么色好看| videos熟女内射| 狂野欧美白嫩少妇大欣赏| 久久国产亚洲av麻豆专区| 麻豆精品久久久久久蜜桃| 丝袜脚勾引网站| 午夜福利视频在线观看免费| 美女xxoo啪啪120秒动态图| 亚洲精品色激情综合| 在线观看一区二区三区激情| 好男人视频免费观看在线| 你懂的网址亚洲精品在线观看| 日本欧美视频一区| 国产乱来视频区| 日日摸夜夜添夜夜爱| 大码成人一级视频| 一本久久精品| 亚洲欧洲日产国产| 日日撸夜夜添| 我的女老师完整版在线观看| 人人妻人人澡人人看| 九九久久精品国产亚洲av麻豆| 亚洲精品亚洲一区二区| 国产成人a∨麻豆精品| 18禁在线播放成人免费| 国产av国产精品国产| 黄色怎么调成土黄色| 最新中文字幕久久久久| 午夜福利在线观看免费完整高清在| 亚洲欧美精品自产自拍| 欧美日韩一区二区视频在线观看视频在线| 国产免费一级a男人的天堂| 少妇的逼水好多| 国产欧美日韩一区二区三区在线 | 国内精品宾馆在线| 桃花免费在线播放| 国产成人精品在线电影| 亚洲av成人精品一区久久| 日韩电影二区| av视频免费观看在线观看| 赤兔流量卡办理| 自线自在国产av| 最近2019中文字幕mv第一页| 国产乱来视频区| 国产精品麻豆人妻色哟哟久久| 亚洲欧美中文字幕日韩二区| 欧美少妇被猛烈插入视频| 成人影院久久| 国产av国产精品国产| 国产免费一区二区三区四区乱码| 丝袜美足系列| 汤姆久久久久久久影院中文字幕| 韩国高清视频一区二区三区| 亚洲欧美日韩卡通动漫| 视频中文字幕在线观看| 日韩电影二区| 伦理电影免费视频| 亚洲av成人精品一二三区| 中文欧美无线码| 久久人妻熟女aⅴ| 最近中文字幕2019免费版| 久久鲁丝午夜福利片| 老司机亚洲免费影院| 一区二区三区四区激情视频| 高清不卡的av网站| 五月玫瑰六月丁香| 日韩人妻高清精品专区| 成人亚洲欧美一区二区av| 97超碰精品成人国产| 日韩精品免费视频一区二区三区 | 色婷婷av一区二区三区视频| 国产精品偷伦视频观看了| 久久热精品热| av.在线天堂| 免费av不卡在线播放| 人妻一区二区av| 国内精品宾馆在线| 校园人妻丝袜中文字幕| 欧美性感艳星| 日日摸夜夜添夜夜爱| 制服丝袜香蕉在线| 男女无遮挡免费网站观看| 久久久国产欧美日韩av| 国产精品国产av在线观看| 日韩,欧美,国产一区二区三区| 91精品三级在线观看| 午夜影院在线不卡| 特大巨黑吊av在线直播| 亚洲欧美一区二区三区国产| 蜜桃国产av成人99| 久久精品国产鲁丝片午夜精品| 欧美丝袜亚洲另类| 久久午夜福利片| 亚洲精品久久久久久婷婷小说| 国产男女超爽视频在线观看| 欧美日韩一区二区视频在线观看视频在线| 精品熟女少妇av免费看| 国产精品成人在线| 五月伊人婷婷丁香| 22中文网久久字幕| 视频在线观看一区二区三区| 美女国产视频在线观看| 天堂8中文在线网| 卡戴珊不雅视频在线播放| 精品少妇久久久久久888优播| 高清欧美精品videossex| 国产一级毛片在线| 日本欧美视频一区| 交换朋友夫妻互换小说| 亚洲国产精品999| av线在线观看网站| 最近2019中文字幕mv第一页| 丝袜在线中文字幕| 老女人水多毛片| 极品少妇高潮喷水抽搐| 在线观看国产h片| 大片电影免费在线观看免费| 在线观看三级黄色| 高清av免费在线| a级毛片免费高清观看在线播放| 亚洲精品国产色婷婷电影| 亚洲av免费高清在线观看| 91精品国产国语对白视频| 久久人人爽av亚洲精品天堂| 亚洲精品av麻豆狂野| av在线播放精品| 欧美少妇被猛烈插入视频| 精品少妇内射三级| 91久久精品国产一区二区三区| 精品久久蜜臀av无| 亚洲国产欧美在线一区| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久国产蜜桃| 日韩大片免费观看网站| 国产国拍精品亚洲av在线观看| 18禁在线播放成人免费| 免费观看的影片在线观看| 欧美精品亚洲一区二区| 啦啦啦中文免费视频观看日本| av又黄又爽大尺度在线免费看| 国产欧美另类精品又又久久亚洲欧美| xxxhd国产人妻xxx| 久久久久网色| 国产精品久久久久久精品古装| 蜜桃国产av成人99| 97在线人人人人妻| 国产一区二区三区综合在线观看 | 日日啪夜夜爽| 啦啦啦啦在线视频资源| 热99国产精品久久久久久7| 午夜老司机福利剧场| 亚洲欧美中文字幕日韩二区| 日本欧美视频一区| 国产白丝娇喘喷水9色精品| 一区二区三区精品91| 久久久精品94久久精品| 激情五月婷婷亚洲| 精品一区二区免费观看| 国产欧美另类精品又又久久亚洲欧美| 午夜91福利影院| 曰老女人黄片| 日韩免费高清中文字幕av| 99热6这里只有精品| 免费黄网站久久成人精品| 亚洲性久久影院| 国产免费一级a男人的天堂| 国产在线免费精品| 日韩精品免费视频一区二区三区 | 91精品伊人久久大香线蕉| 大香蕉97超碰在线| 免费观看性生交大片5| 国产精品一区二区在线不卡| 日本av手机在线免费观看| 午夜91福利影院| 97精品久久久久久久久久精品| 久热久热在线精品观看| 亚洲欧美色中文字幕在线| 午夜视频国产福利| 欧美日本中文国产一区发布| 国产精品一区二区在线观看99| 亚洲精品日韩在线中文字幕| 亚州av有码| 成人18禁高潮啪啪吃奶动态图 | 婷婷色综合www| 欧美日韩av久久| 国产无遮挡羞羞视频在线观看| 18禁观看日本| 新久久久久国产一级毛片| 两个人免费观看高清视频| 欧美成人午夜免费资源| 国产乱人偷精品视频| 天天影视国产精品| 免费黄网站久久成人精品| 91精品国产九色| 曰老女人黄片| 亚洲欧美日韩另类电影网站| 国产爽快片一区二区三区| 蜜桃国产av成人99| 新久久久久国产一级毛片| 国产成人精品福利久久| 国产精品偷伦视频观看了| 免费av中文字幕在线| 国产精品女同一区二区软件| 黄色怎么调成土黄色| 好男人视频免费观看在线| 亚洲精品中文字幕在线视频| a级片在线免费高清观看视频| 你懂的网址亚洲精品在线观看| 午夜av观看不卡| 少妇人妻久久综合中文| 亚洲av福利一区| 久久久久久久精品精品| 国产精品99久久99久久久不卡 | 亚洲精品一区蜜桃| 久久久久久久亚洲中文字幕| 国产不卡av网站在线观看| 在线播放无遮挡| 国产成人精品久久久久久| 久久久久视频综合| 亚洲精品av麻豆狂野| 亚洲国产成人一精品久久久| kizo精华| 在线天堂最新版资源| av免费在线看不卡| 99久久中文字幕三级久久日本| 综合色丁香网| 久久99热6这里只有精品| 夜夜爽夜夜爽视频| 国产有黄有色有爽视频| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 亚洲国产欧美日韩在线播放| 国产毛片在线视频| 亚洲欧美一区二区三区黑人 | 成人18禁高潮啪啪吃奶动态图 | 亚洲三级黄色毛片| 天天躁夜夜躁狠狠久久av| 九色成人免费人妻av| 少妇的逼好多水| 草草在线视频免费看| 欧美性感艳星| 老司机影院成人| 十八禁高潮呻吟视频| 不卡视频在线观看欧美| 22中文网久久字幕| 亚洲av.av天堂| 日日爽夜夜爽网站| 久久精品久久久久久久性| 在线观看免费视频网站a站| 国产成人精品久久久久久| 777米奇影视久久| 一二三四中文在线观看免费高清| 亚洲国产色片| 亚洲精品乱码久久久v下载方式| 人妻 亚洲 视频| 免费高清在线观看日韩| 超色免费av| 777米奇影视久久| 热re99久久国产66热| 天堂俺去俺来也www色官网| 成人综合一区亚洲| 国产成人freesex在线| av网站免费在线观看视频| 日韩制服骚丝袜av| 亚洲情色 制服丝袜| 国产精品久久久久久精品电影小说| 亚洲不卡免费看| 中国美白少妇内射xxxbb| 成人毛片60女人毛片免费| 七月丁香在线播放| 综合色丁香网| 久久精品久久久久久久性| 亚洲性久久影院| 最新的欧美精品一区二区| 国产精品一区二区在线不卡| 亚洲精品中文字幕在线视频| 中文字幕久久专区| 国产一区二区三区综合在线观看 | 十八禁网站网址无遮挡| 婷婷色综合大香蕉| 91精品国产国语对白视频| 亚洲国产精品国产精品| 九九在线视频观看精品| 国产一区二区在线观看日韩| 人人妻人人添人人爽欧美一区卜| 国产精品国产三级国产专区5o| 岛国毛片在线播放| 精品久久蜜臀av无| 在线免费观看不下载黄p国产| 精品人妻熟女毛片av久久网站| 国产欧美日韩综合在线一区二区| 欧美日韩视频高清一区二区三区二| 人妻少妇偷人精品九色| 在线亚洲精品国产二区图片欧美 | 看非洲黑人一级黄片| 在线观看免费高清a一片| 18禁动态无遮挡网站| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 我的女老师完整版在线观看| 黑人欧美特级aaaaaa片| 简卡轻食公司| 色5月婷婷丁香| 久久精品久久久久久噜噜老黄| 国产黄片视频在线免费观看| 亚洲精品成人av观看孕妇| 成人午夜精彩视频在线观看| 欧美亚洲 丝袜 人妻 在线| av在线播放精品| 国产深夜福利视频在线观看| 中文字幕人妻丝袜制服| 亚洲国产av新网站| 成人影院久久| 女性被躁到高潮视频| 免费少妇av软件| 超色免费av| 99久久综合免费| 日韩电影二区| 久久免费观看电影| 全区人妻精品视频| 女人精品久久久久毛片| 亚洲高清免费不卡视频| 久久久久久久国产电影| 女性生殖器流出的白浆| 亚洲国产av影院在线观看| 午夜老司机福利剧场| 亚洲三级黄色毛片| 大话2 男鬼变身卡| av一本久久久久| 国产日韩欧美视频二区| 欧美老熟妇乱子伦牲交| av免费观看日本| 亚洲精品国产av成人精品| 日本欧美视频一区| 久久国产精品大桥未久av| 日本vs欧美在线观看视频| 能在线免费看毛片的网站| 新久久久久国产一级毛片| 精品国产露脸久久av麻豆| 久久精品国产自在天天线| 日本爱情动作片www.在线观看| 91精品国产九色| 少妇高潮的动态图| av福利片在线| 国产成人a∨麻豆精品| 久久久精品免费免费高清| 人妻夜夜爽99麻豆av| 免费久久久久久久精品成人欧美视频 | 一级毛片aaaaaa免费看小| 永久网站在线| 午夜免费鲁丝| 一级爰片在线观看| 汤姆久久久久久久影院中文字幕| 一区二区三区乱码不卡18| 国产又色又爽无遮挡免| 汤姆久久久久久久影院中文字幕| 男女高潮啪啪啪动态图| 亚洲四区av| 另类亚洲欧美激情| 久久精品国产自在天天线| 国产精品熟女久久久久浪| 在线精品无人区一区二区三| 97超视频在线观看视频| 美女视频免费永久观看网站| 成人亚洲欧美一区二区av| 国产成人91sexporn| 最近手机中文字幕大全| 赤兔流量卡办理| 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 亚洲经典国产精华液单| 国产片内射在线| 精品亚洲乱码少妇综合久久| 婷婷色麻豆天堂久久| 久久 成人 亚洲| 亚洲国产毛片av蜜桃av| 亚洲国产日韩一区二区| 亚洲伊人久久精品综合| 99久国产av精品国产电影| 欧美一级a爱片免费观看看| 午夜av观看不卡| 大香蕉97超碰在线| 国产在线视频一区二区| 久久久国产一区二区| 国产爽快片一区二区三区| 欧美日韩综合久久久久久| 久久久国产一区二区| 天堂俺去俺来也www色官网| 亚洲欧美一区二区三区黑人 | 欧美 日韩 精品 国产| 制服丝袜香蕉在线| 久久97久久精品| 一区二区日韩欧美中文字幕 | 插逼视频在线观看| 国产成人精品婷婷| 秋霞在线观看毛片| 久久鲁丝午夜福利片| 亚洲国产欧美在线一区| 韩国高清视频一区二区三区| 熟女电影av网| 亚洲av电影在线观看一区二区三区| 乱人伦中国视频| 成人毛片60女人毛片免费| 成人漫画全彩无遮挡| 中文字幕免费在线视频6| 三级国产精品欧美在线观看| 99九九线精品视频在线观看视频| 久久久久久久国产电影| 欧美三级亚洲精品| 欧美变态另类bdsm刘玥| 国产在线免费精品| 十分钟在线观看高清视频www| 国产精品99久久99久久久不卡 | 99re6热这里在线精品视频| 国产成人精品婷婷| 蜜桃久久精品国产亚洲av| 午夜福利在线观看免费完整高清在| av在线播放精品| 国产又色又爽无遮挡免| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 精品熟女少妇av免费看| 美女内射精品一级片tv| 欧美激情 高清一区二区三区| 国产亚洲午夜精品一区二区久久| 国产深夜福利视频在线观看| 婷婷色综合大香蕉| 日韩一区二区三区影片| 两个人的视频大全免费| 一级毛片我不卡| 精品人妻熟女毛片av久久网站| 大片电影免费在线观看免费| 国产欧美日韩一区二区三区在线 | 久久青草综合色| 日韩视频在线欧美| 人人妻人人澡人人看| 另类精品久久| 国产精品国产三级专区第一集| 久久久久精品性色| 99精国产麻豆久久婷婷| 国产成人a∨麻豆精品| 91精品三级在线观看| a级毛片免费高清观看在线播放| 国产成人精品一,二区| 成人手机av| 曰老女人黄片| 精品国产一区二区三区久久久樱花| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 免费播放大片免费观看视频在线观看| 亚洲欧美一区二区三区国产| 天堂中文最新版在线下载| 久久韩国三级中文字幕| 99国产精品免费福利视频| 欧美性感艳星| 日韩欧美一区视频在线观看| 色婷婷av一区二区三区视频| 精品99又大又爽又粗少妇毛片| 一区在线观看完整版| 亚洲国产色片| 你懂的网址亚洲精品在线观看| 亚洲欧美日韩卡通动漫| 国产av一区二区精品久久| xxx大片免费视频| 国产色爽女视频免费观看| 国产成人91sexporn| 久久久久久人妻| 日本色播在线视频| 国产乱来视频区| 91午夜精品亚洲一区二区三区| 国产亚洲精品久久久com| 在线天堂最新版资源| 99久久综合免费| www.av在线官网国产| 欧美三级亚洲精品| 伦理电影免费视频| 亚洲精品久久成人aⅴ小说 | 国产片内射在线| 成年女人在线观看亚洲视频| 一本色道久久久久久精品综合| 18禁观看日本| 男人操女人黄网站| 亚洲精品456在线播放app| 黄片无遮挡物在线观看| 日韩亚洲欧美综合| 99国产综合亚洲精品| 久久99热这里只频精品6学生| 久久久欧美国产精品| 国产片内射在线| 免费看不卡的av| 欧美精品国产亚洲| 国产亚洲精品第一综合不卡 | 免费高清在线观看视频在线观看| 丰满少妇做爰视频| 久久精品人人爽人人爽视色| 免费观看性生交大片5| 欧美激情 高清一区二区三区| 欧美另类一区| 国产av一区二区精品久久| 啦啦啦在线观看免费高清www| 国产深夜福利视频在线观看| 26uuu在线亚洲综合色| 久久久久久人妻| 视频区图区小说| 欧美成人午夜免费资源| 最新中文字幕久久久久| 18禁动态无遮挡网站| 又大又黄又爽视频免费| 如何舔出高潮| av福利片在线| 91精品国产九色| 亚洲天堂av无毛| 晚上一个人看的免费电影| 国产熟女午夜一区二区三区 | 日产精品乱码卡一卡2卡三| 51国产日韩欧美| 伊人亚洲综合成人网| 五月天丁香电影| av网站免费在线观看视频| 在线 av 中文字幕| 熟女电影av网| 老司机亚洲免费影院| 亚洲av欧美aⅴ国产| 国产精品国产三级国产专区5o| 亚洲精品,欧美精品| 人妻系列 视频| 中文欧美无线码| 丝袜在线中文字幕| av免费观看日本| 纯流量卡能插随身wifi吗| 成人国产麻豆网| 久久精品国产a三级三级三级| 国产成人免费观看mmmm| 国产高清三级在线| 我的女老师完整版在线观看| 乱人伦中国视频| 色吧在线观看| 久久久久久久亚洲中文字幕| 一本—道久久a久久精品蜜桃钙片| 国产色婷婷99| 色哟哟·www| 国产不卡av网站在线观看| 黑人欧美特级aaaaaa片| 国产一区二区三区综合在线观看 | 国精品久久久久久国模美| 女性生殖器流出的白浆| 一区二区三区免费毛片| 免费观看av网站的网址| 女性被躁到高潮视频| 欧美精品人与动牲交sv欧美| 久久 成人 亚洲| 黑丝袜美女国产一区| 国语对白做爰xxxⅹ性视频网站| 夫妻性生交免费视频一级片| 久久午夜福利片| 欧美xxⅹ黑人| 午夜激情av网站| 男女啪啪激烈高潮av片| 秋霞在线观看毛片| 91aial.com中文字幕在线观看| 999精品在线视频| 国产极品粉嫩免费观看在线 | 国产伦理片在线播放av一区| 精品少妇久久久久久888优播| 色94色欧美一区二区| 中文字幕精品免费在线观看视频 | 满18在线观看网站|