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

    半干旱地區(qū)黑碳氣溶膠和含碳氣體特征及來源

    2012-12-26 02:22:28曹賢潔黃建平張北斗蘭州大學大氣科學學院半干旱氣候變化教育部重點實驗室甘肅蘭州730000中國人民解放軍92543部隊河南濟源454650
    中國環(huán)境科學 2012年8期
    關鍵詞:觀測站氣溶膠氣流

    陳 霖 ,張 鐳 *,張 磊 ,2,曹賢潔 ,黃建平 ,張 武 ,張北斗 (.蘭州大學大氣科學學院,半干旱氣候變化教育部重點實驗室,甘肅 蘭州 730000;2.中國人民解放軍92543部隊,河南 濟源 454650)

    半干旱地區(qū)黑碳氣溶膠和含碳氣體特征及來源

    陳 霖1,張 鐳1*,張 磊1,2,曹賢潔1,黃建平1,張 武1,張北斗1(1.蘭州大學大氣科學學院,半干旱氣候變化教育部重點實驗室,甘肅 蘭州 730000;2.中國人民解放軍92543部隊,河南 濟源 454650)

    利用2009年9月~2010年6月蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)多角度吸收分光光度計(MAAP-5012)觀測數(shù)據(jù)、CO和CO2氣體成分混合比數(shù)據(jù),分析了西北半干旱地區(qū)黑碳氣溶膠和含碳氣體特征、影響源地,以及影響黑碳濃度的排放物類型.結果表明BC、CO、CO2平均濃度分別為1.75μg/m3、601.71×10-9、387.78×10-6.利用后向軌跡模式將從觀測站西部和東部輸送過來的氣流區(qū)分開,氣流從東部來時,BC、CO、CO2濃度分別為1.38μg/m3、462.79×10-9、383.03×10-6;氣流從西部來時,BC、CO、CO2濃度分別為2.2μg/m3、768.38×10-9、393.47×10-6.對500m、1500m、3000m高度氣流來向的發(fā)源地進行聚類分析,發(fā)現(xiàn)3個高度上氣流從中東、中亞及歐洲區(qū)域傳輸過來時,BC、 CO、CO2濃度較高,△B C/△C O、△ CO/△ CO2值較大,說明燃料燃燒效率較低;氣流從我國華北華中地區(qū)傳輸過來時,BC、CO、CO2濃度較低,△B C/△ CO、△ CO/△ CO2值較小,表明燃料燃燒效率較高.

    半干旱地區(qū);黑碳氣溶膠;后向軌跡;聚類分析

    黑碳氣溶膠作為吸收太陽輻射的主要氣溶膠成分,近年來引起國內(nèi)外很多學者關注[1-6].黑碳主要是由于物質(zhì)的不完全燃燒產(chǎn)生的,如工業(yè)生產(chǎn)、汽車尾氣、居民燃火、生物質(zhì)燃燒等.黑碳氣溶膠的吸收性能占大氣氣溶膠總吸收性能的90%以上[7],所以一般也用黑碳氣溶膠的吸收系數(shù)作為大氣氣溶膠的吸收系數(shù).黑碳氣溶膠能吸收太陽短波輻射,對大氣及大氣層頂造成正輻射強迫, 同時由于其吸收和散射作用,使得到達地面的太陽短波輻射減少,對地表造成負的輻射強迫作用[8].大氣中黑碳氣溶膠可以通過云-氣溶膠相互作用,對云層的厚度、云層高度、云量和云頂反照率造成影響,從而造成半直接輻射強迫效應和間接輻射強迫效應[9].近年研究發(fā)現(xiàn),在南極[10]、北極[11]地區(qū)也觀測到了黑碳氣溶膠,且濃度較大.

    大氣中的黑碳氣溶膠主要來源于火山噴發(fā)和化石燃料燃燒.IPCC第 4次評估報告[12]指出,當前全球平均每年排放黑碳約8.0Tg/a,其中化石燃料燃燒排放為 4.6Tg/a,生物質(zhì)燃燒排放為3.3Tg/a.我國大城市的日常生活中的化石燃料消耗量較大,城市工業(yè)、汽車尾氣和居民活動產(chǎn)生的黑碳較多,鄉(xiāng)村燃燒秸稈較多[13-14],也對大氣中的黑碳濃度有較大的影響.

    目前對人為排放氣溶膠的研究多集中于OC(有機碳)和 BC(黑碳)之間的關系及其季節(jié)變化[15-16]、BC的混合態(tài)[17]及硫酸鹽氣溶膠和有機氣溶膠的形成[18]方面,而對 BC與燃燒示蹤物(CO、CO2等)的研究較少,BC與CO、CO2的相互關系可以反映物質(zhì)的燃燒效率及 BC在大氣中得傳輸效率,可較好地表征不同區(qū)域的燃燒物燃燒效率及黑碳的傳輸效率.

    我國半干旱地區(qū)由于沙塵天氣及重工業(yè)污染的雙重影響,使得黑碳氣溶膠研究變得相對復雜.本課題組前期工作分析了蘭州遠郊區(qū)2007年1月1日至2009年8月28日黑碳氣溶膠濃度特征,黑碳濃度日變化及季節(jié)變化特征,并分析了沙塵過程對黑碳濃度的影響,但沒有分析周邊區(qū)域?qū)τ^測站黑碳濃度的影響,且沒有結合其他物質(zhì)進行分析.本文在此基礎上,利用蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)多角度吸收分光光度計(MAAP-5012)觀測數(shù)據(jù)進行半干旱地區(qū)黑碳氣溶膠變化特征的研究,同時結合CO和CO2數(shù)據(jù),通過后向軌跡分析影響觀測站的氣流的起源地、起源高度和燃燒物質(zhì)情況,并將觀測站西部和東部輸送來的氣流區(qū)分開,研究觀測站兩個主導風向?qū)谔細馊苣z及含碳氣體的影響.

    1 觀測儀器與資料

    1.1 觀測儀器

    黑碳氣溶膠的觀測使用Thermo Scientific公司生產(chǎn)的多角度吸收分光光度計(Model-5012).該儀器利用光束通過氣溶膠樣品后的衰減量與濾膜上黑碳氣溶膠質(zhì)量濃度之間的線性關系,在測量光透過濾膜后的衰減量后,通過線性關系反算出氣溶膠中黑碳的濃度.測量儀器腔體內(nèi)增加了光散射的測量,可用于校正顆粒物累積在濾帶上造成的多次散射的影響.MAAP-5012型黑碳儀使用濾膜玻璃纖維收集大氣中氣溶膠,其觀測波長為670nm.

    采用 ML9830型一氧化碳分析儀觀測 CO,儀器由微處理器控制,采用非色散紅外(NDIR)氣體過濾相關光度測定法,可較為精確、穩(wěn)定地測定大氣中CO的含量.用ML9820型二氧化碳分析儀進行CO2觀測,ML9820型二氧化碳分析儀由計算機控制,采用非分散紅外相關(GFC)光學測量技術,使CO和H2O的干擾降到最小,可較為精確、穩(wěn)定地測量CO2濃度.

    1.2 資料處理與質(zhì)量控制

    2009年9月~2010年6月在蘭州大學半干旱氣候與環(huán)境觀測站連續(xù)觀測了BC濃度和CO、CO2氣體混合比,BC濃度數(shù)據(jù)時間步長為1min,CO、CO2數(shù)據(jù)時間步長為 5min,實驗室每天記錄天氣情況、云量及周邊局地人為活動的變化.

    黑碳數(shù)據(jù)使用范圍檢查和時變檢查[19]的方法進行質(zhì)量控制.

    MAAP-5012黑碳儀通過測量濾膜上的黑碳對光的衰減量來反算黑碳的質(zhì)量濃度,其反演公式為:

    式中:MBC為黑碳質(zhì)量;ω0為單次散射反照率;X為透過率;A為濾膜上黑碳斑點面積;σBC為黑碳的質(zhì)量吸收系數(shù),m2/g;VOL為氣流流量;BC為儀器輸出的黑碳質(zhì)量濃度.式中X為觀測值,A、VOL值為固定值,因此訂正時需要考慮σBC的取值,儀器設定值為 6.6m2/g,但不同地區(qū) σBC的值不同,數(shù)據(jù)處理時參考歐洲西南部與非洲撒哈拉沙漠臨近地區(qū)用 MAAP-5012黑碳儀得出黑碳質(zhì)量吸收系數(shù)(約 10.3m2/g[21-23]),這一區(qū)域氣溶膠的吸收特性既受北部歐洲地區(qū)工業(yè)排放的影響,又有南部撒哈拉沙漠的沙塵影響,與我國西北半干旱地區(qū)環(huán)境較為相似.本文中σBC取10.3 m2/g.訂正時用下式:

    CBCC為訂正后的黑碳濃度.

    CO、CO2觀測儀器使用零/跨檢查、多點較準和定期對儀器做檢查與維護進行質(zhì)量控制,同時使用數(shù)據(jù)時通過校準數(shù)據(jù)所得到的多點校準方程進行訂正[24].

    通過BC與CO、CO2數(shù)據(jù)質(zhì)量控制和訂正后,BC與CO、CO2數(shù)據(jù)同時完好的共有253d,占總天數(shù)303d的83.5%,若剔除降雨和沙塵暴等特殊天氣,總天數(shù)為280d.

    2 結果分析

    2.1 兩種黑碳儀MAAP-5012與AE-31觀測結果比較分析

    黑碳儀MAAP-5012和AE-31都可用于測量大氣中的黑碳濃度,但MAAP在觀測時同時接收濾膜上氣溶膠的散射信號,并將散射信號剔除.利用蘭州大學半干旱氣候與環(huán)境觀測站2009年9月~2010年 3月 MAAP-5012(670nm)和AE-31(660nm波段)的同步觀測數(shù)據(jù)進行對比分析.所選波段差異很小,觀測波長造成的差異可以忽略不計.

    由圖1可見,MAAP-5012與AE-31測定的BC濃度5min數(shù)據(jù)的均值相關系數(shù)為0.93,日變化趨勢相關度略差;而日均值相關系數(shù)為0.98,說明長時間平均結果基本一致.

    2.2 黑碳、一氧化碳及二氧化碳變化特征

    觀測站地面風速在秋、冬季節(jié)較小,春、夏季較大.BC和 CO變化也基本為秋冬季較大,春夏季相對較小,這與冬季采暖燃煤量較大有關;而BC與 CO2變化聯(lián)系較弱,CO2冬季較大,與冬季燃煤有關,但日變化特征不明顯,可能與天氣過程有關(圖2).

    圖1 MAAP-5012與AE-31黑碳濃度相關性Fig.1 Correlation of BC between MAAP-5012 and AE-31

    圖2 地面風速風向,BC和CO及CO2日均值變化Fig.2 Daily mean of wind speed, wind direction, BC-CO and BC-CO2 concentration

    由圖3可見,BC平均濃度在一周中星期一最小,為1.53μg/m3;星期六最大,為1.77μg/m3.BC濃度星期一~星期四較星期五~星期日小.可能與周末居民活動黑碳排放增多有關,觀測站地處蘭州市遠郊區(qū),星期一~星期四附近經(jīng)過車輛相對較少,周末出行車輛增多,使得附近黑碳排放相對星期一~星期四增多.

    圖3 黑碳濃度周變化Fig.3 Weekly variation of BC Conentralion

    圖4 BC、CO、CO2濃度及地面風速風向日變化基本情況Fig.4 Diurnal variation of BC,CO,CO2 and surface wind speed, wind direction

    由圖4可見,觀測站BC日變化呈現(xiàn)雙峰結構,10:00和 20:00為峰值,對應濃度分別為2.18,1.88μg/m3,這與上下班高峰期及早晨逆溫層發(fā)展和地面風速較低有關;CO和CO2日變化呈現(xiàn)單峰結構,13:00為峰值,最大值分別為790.72× 10-9和 391.86×10-6,可能與附近村莊人為活動有關;地面風速在08:00達最小值,為2.02m/s.

    2.3 氣流后向軌跡分析

    考慮到觀測站近地面主導風向為東南方向,且觀測站以東至東南沿海地區(qū)和觀測站以西至西亞地區(qū)無論地形、植被還是工業(yè)程度都有較大差異,所以將氣團按發(fā)源地經(jīng)度分為東、西兩部分,利用 HYSPLIT-4模式[25]研究近地面氣流來向不同對BC、CO及CO2的影響.一般研究經(jīng)過城市的氣流造成的影響時,選取50~300m較為合適,由于SACOL站西部蘭州市屬河谷地形城市,需選取較低高度進行研究,所以通過后向軌跡模式給出50m高度處氣流96h后向軌跡,每d 1條軌跡,每條后向軌跡起始時間取北京時間(與儀器觀測時間同步)00:00.2009年9月~2010年6月共303條后向軌跡,其中來自觀測站東部的氣流占比略多,為 51%,西部傳輸過來的氣流占 49%,圖略.

    由表 1可見,近地面氣流從西部傳輸過來時,BC、CO、CO2濃度分別為 2.20μg/m3、768.38×10-9、393.47×10-6,比氣流從東部來時高,這與觀測站西部區(qū)域排放黑碳濃度較高有關.由圖5可見,氣流來自不同方向時,濃度日變化有較大差異:氣流從西部傳輸過來時,BC濃度峰值和次峰值出現(xiàn)在10:00和20:00,濃度分別為2.75,2.38μg/m3,這是由于氣流從西部傳輸過來時,經(jīng)過蘭州市市區(qū)近地層,城市上下班高峰期對 BC濃度峰值出現(xiàn)時間有一定影響,而氣流從東部傳輸過來時,BC濃度峰值和次峰值出現(xiàn)在 09:00和 21:00,濃度分別為 1.71, 1.55μg/m3;氣流從西部傳輸過來時,CO白天維持較高水平,峰值出現(xiàn)在 14:00,濃度為 955× 10-9,這與城市白天工業(yè)污染較重有關,而從東部傳輸過來時 CO白天平均濃度較低,但變化范圍較大,峰值出現(xiàn)在 13:00,濃度為 732×10-9;氣流從西部、東部傳輸過來時 CO2峰值都出現(xiàn)在 13:00,濃度分別為 402×10-6和 383×10-6,濃度值差別較大.

    表1 不同氣流來向BC、CO、CO2濃度Table 1 BC, CO and CO2 concentration of different airflow direction

    圖5 不同氣流來向BC、CO、CO2濃度日變化Fig.5 Diurnal variation of BC,CO,CO2 conentralion of different airflow direction

    考慮到SACOL周圍地區(qū)離地500m高度層氣溶膠混合程度較好,1500m高度一般為混合層頂,3000m高度適合分析高空氣團大尺度運行特征,所以選取500m、1500m、3000m做后向軌跡的模擬.每d 1條后向軌跡,每條后向軌跡起始時間取北京時(與儀器觀測時間同步)00:00.分析中使用馬氏距離[26]的 K-均值聚類方法[27],將后向軌跡路徑按起點分為6個類別,并對每個類別的后向軌跡做“convhull”[28](即按軌跡的邊界做出包裹線)處理.由于單一的后向軌跡有不確定性且直觀性差[26],本文作出不同高度后向軌跡圖,相對單一高度層后向軌跡圖更為直觀.圖6為各個高度上每個類別后向軌跡平面圖,3個高度上后向軌跡起源基本一致,可大致劃分為中東、中亞、歐洲、俄羅斯西部和北部、中國華北華中等區(qū)域.

    在500m、1500m高度層上,從中亞地區(qū)、西西伯利亞平原和我國華北華中地區(qū)傳輸過來的氣流所占比例較大;而3000m高度層上,主要是從觀測站西部傳輸過來的氣流,東部傳輸過來的氣流所占比例很少,這主要是由于高空受西風氣流主導所致.由表2可見,BC、CO、CO2平均濃度的變化一致,BC高值對應的CO、CO2也較高.在500m、1500m高度層上,氣流從歐洲東部、歐洲中部、中東及中亞地區(qū)傳輸過來時,觀測的 BC及CO濃度明顯較高,高于2.0μg/m3,而從我國華北華中地區(qū)傳輸過來時,BC和CO濃度較低;在3000m高度層上,氣流從歐洲地區(qū)傳輸過來時,觀測的BC及CO濃度明顯較高,而氣流從西西伯利亞平原傳輸過來時,BC和CO濃度較低.

    △CO/△CO2值(即CO-CO2離散點對比的斜率值)可反映排放源的燃燒效率,值越小說明燃燒越充分,一般化石燃料燃燒(如石油)△CO/△CO2值較小,而生物質(zhì)燃燒△CO/△CO2值較大[29],通過與排放清單的對比,也可以作為確定區(qū)域主要燃燒物質(zhì)的指標[30],△BC/△CO值可以很好地反映黑碳排放后傳輸?shù)竭吔鐚拥男室约皬倪吔鐚觽鬏數(shù)阶杂纱髿獾男蔥31].

    從3個高度層上看,氣流從歐洲、中東及中亞地區(qū)傳輸過來時,BC濃度較高,同時△BC/△CO、△CO/△CO2值明顯較高,說明這些地區(qū)物質(zhì)燃燒效率較低,物質(zhì)燃燒不完全,導致黑碳排放較多,可能由于當?shù)厣镔|(zhì)燃燒所占比例較多,長距離傳輸混合作用也有一定影響;而從西西伯利亞平原、我國華北華中地區(qū)傳輸過來時,BC濃度較低,△BC/△CO、△CO/△CO2值較低,說明這些區(qū)域物質(zhì)燃燒效率較高,尤其我國華北華中地區(qū),雖然能源消耗較大,但從結果來看產(chǎn)生的黑碳較少.

    3 結論

    3.1 BC、CO和CO2平均濃度分別為1.75μg/m3、601.71×10-9、387.78×10-6.BC濃度日變化呈現(xiàn)雙峰結構,10:00和 20:00為最大值,峰值分別為2.18μg/m3和1.88μg/m3,CO和CO2濃度日變化呈現(xiàn)單峰結構,13:00為最大值,分別為 790.72×10-9和 391.86×10-6.BC平均濃度在一周中星期一達到最小值,為1.53μg/m3;星期六較大,為1.77μg/m3,這與周末居民活動所致黑碳排放增多有關.

    圖6 500m、1500m、3000m高度后向軌跡聚類分析Fig.6 Clustering analysis of backward trajectory in 500m, 1500m and 3000m

    表2 各氣流束比例和對應BC、CO、CO2濃度和 △B C/△ CO、 △B C/△ CO2、△ CO/△ CO2比值Table 2 The percentage of airflow and concentration of BC, CO, CO2 and amount of △B C/△ C O、 △ BC/△ CO2、△C O/△C O2

    3.2 氣流從觀測站東部傳輸過來時,BC、CO、CO2濃度分別為 1.38μg/m3、462.79×10-9、383.03×10-6,BC 濃度峰值和次峰值出現(xiàn)在09:00和21:00,CO、CO2濃度峰值出現(xiàn)在13:00;氣流從西部傳輸過來時,BC、CO、CO2濃度分別為 2.20μg/m3、768.38×10-9、393.47×10-6,濃度較高, BC濃度峰值和次峰值出現(xiàn)在10:00和20:00,CO濃度白天維持較高水平,峰值出現(xiàn)在14:00,CO2峰值出現(xiàn)在13:00.

    3.3 500m、1500m、3000m高度上氣流從中東、歐洲及中亞地區(qū)上空傳輸過來時,BC濃度較大,且△BC/△CO、△CO/△CO2值較高,說明這些地區(qū)燃燒物燃燒效率較低;而氣流從我國華北華中地區(qū)傳輸過來時 BC濃度較低,且△BC/△CO、△CO/△CO2值較小,說明東部地區(qū)燃燒物燃燒效率較高.

    [1] 蔡子穎,韓素芹,黃 鶴,等.天津夏季黑碳氣溶膠及其吸收特性的觀測研究 [J]. 中國環(huán)境科學, 2011,31(5):719-723.

    [2] Huang J P, Fu Q, Zhang W, et al. Dust and black carbon in seasonal snow across Northern China [J]. Bull. Amer. Meteor. Soc., 2011,92:175-181.

    [3] Judith C C, John G W, Douglas H L, et al. PM2.5 source profiles for black and organic carbon emission inventories [J]. Atmos. Environ., 2011,45:5407-5414.

    [4] Gilardoni S, Vignati E, Wilson J. Using measurements for evaluation of black carbon modeling [J]. Atmos. Chem. Phys., 2011,11:439-455.

    [5] Jones G S, Christidis N, Stott P A. Detecting the influence of fossil fuel and bio-fuel black carbon aerosols on near surface temperature changes [J]. Atmos. Chem. Phys., 2011,11:799-816.

    [6] 張 磊,張 鐳,張丁玲,等.蘭州遠郊區(qū)黑碳氣溶膠濃度特征[J]. 中國環(huán)境科學, 2011,31(8):1248-1255.

    [7] Hansen A D A, Kapustin V N, Kopeikin V M, et a1. Optical absorption by aerosol black carbon and dust in desert region of central Asia [J]. Atmos. Environ., 1993,27(16):2527-2531.

    [8] 張 華,馬井會,鄭有飛.黑碳氣溶膠輻射強迫全球分布的模擬研究 [J]. 大氣科學, 2008,32(5):1147-1158.

    [9] Koch D, Del G A D. Black carbon semi-direct effects on cloud cover: review and synthesis [J]. Atmos. Chem. Phys., 2010,10: 7685-7696.

    [10] Graf H F, Shirsat S V, Oppenheimer C, et al. Continental scale Antarctic deposition of sulphur and black carbon from anthropogenic and volcanic sources [J]. Atmos. Chem. Phys., 2010,10:2457-2465.

    [11] Huang L, Gong S L, Sharma S, et al. A trajectory analysis of atmospheric transport of black carbon aerosols to Canadian high Arctic in winter and spring (1990–2005) [J]. Atmos. Chem. Phys., 2010,10:5065-5073.

    [12] Christensen J H, Hewitson B, Busuioc A, et a1. IPCC climate change 2007: The physical scientific basis [M]. London: Cambridge University Press, 2007,131-217.

    [13] 曹國良,鄭方成,王亞強.中國大陸生物質(zhì)燃燒排放的TSP, PM10, PM2.5清單 [J]. 過程工程學報, 2004.4:700-704.

    [14] Streets D G, Bond T C, Carmichael G R, et al. An inventory of gaseous and primary aerosol emissions in Asia in the year 2000 [J]. J. Geophys. Res., 2003,108:8809.

    [15] Rehman I H, Ahmed T, Praveen P S, et al. Black carbon emissions from biomass and fossil fuels in rural India [J]. Atmos. Chem. Phys., 2011,11:7289-7299.

    [16] Verma R L, Sahu L K, Kondo Y, et al. Temporal variations of black carbon in Guangzhou, China, in summer 2006 [J]. Atmos. Chem. Phys., 2010,10:6471-6485.

    [17] McMeeking G R, Morgan W T, Flynn M, et al. Black carbon aerosol mixing state, organic aerosols and aerosol optical properties over the United Kingdom [J]. Atmos. Chem. Phys., 2011,11:9037-9052.

    [18] Xiao R, Takegawa N, Kondo Y, et al. Formation of submicron sulfate and organic aerosols in the outflow from the urban region of the Pear River Delta in China [J]. Atmos. Environ., 2009, 43:3754-3763.

    [19] 熊安元.北歐氣象觀測資料的質(zhì)量控制 [J]. 氣象科技, 2003, 31(5):314-320.

    [20] 楊 溯,張 武,史晉森,等.半干旱地區(qū)黑碳氣溶膠特征初步分析 [J]. 氣候與環(huán)境研究, 2010,15(6):756-764.

    [21] Fernandez C R, Rodr?guez S, de la Rosa J, et al. Ultrafine particle formation in the inland sea breeze airflow in Southwest Europe [J]. Atmos. Chem. Phys., 2010, 10: 9615-9630.

    [22] Reche C, Querol X, Alastuey A, et al. New considerations for PM, black carbon and particle number concentration for air quality monitoring across different European cities [J]. Atmos. Chem. Phys., 2011,11:6207-6227.

    [23] Pandolfi M, Cusack M, Alastuey A, et al. Variability of aerosol optical properties in the Western Mediterranean Basin [J]. Atmos. Chem. Phys., 2011,11:8189-8203.

    [24] 林偉立,徐曉斌,于大江,等.龍鳳山區(qū)域大氣本底臺站反應性氣體觀測質(zhì)量控制 [J]. 氣象, 2009,35(11):93-100.

    [25] Draxler R R. Description of the HYSPLIT-4 Modeling System [M]. NOAA Technical Memorandum, 1997,24:1-28.

    [26] Roy D M, Delphine J R, Desire L M. The Mahalanobis distance [J]. Chemometrics and Intelligent Laboratory Systems, 2000,50:1-18.

    [27] Anderberg M R. Cluster Analysis for Applications [M]. New York: Academic Press, 1973:105-108.

    [28] Matlab 7 Mathematics [M]. MathWorks, 2010:296-302.

    [29] Yokelson R J, Bertschi I T, Christian T J, et al. Trace gas measurements in nascent, aged, and cloud processed smoke from African savanna fires by airborne Fourier transform infrared spectroscopy (AFTIR) [J]. J. Geophys. Res., 2003,108:8478.

    [30] Takegawa N, Kondo Y, Koike M, et al. Removal of NOxand NOyin Asian outflow plumes: Aircraft measurements over the western Pacific in January 2002 [J]. J. Geophys. Res., 2004, 109 (D23S04):1-16.

    [31] Park S S, Bae M S, Schauer J J, et al. Evaluation of the TMO and TOT methods for OC and EC measurement and their characteristics in PM2.5 at an urban site of Korea during ACE-Asia [J]. Atmos. Environ., 2005,39:5101-5112.

    Characteristics of black carbon aerosol and carbonaceous gases and their emission sources in semi-arid region.

    CHEN Lin1, ZHANG Lei1*, ZHANG Lei1,2, CAO Xian-jie1, HUANG Jian-ping1, ZHANG Wu1, ZHANG Bei-dou1(1.Key Laboratory for Semi-Arid Climate Change, Ministry of Education, College of Atmospheric Science, Lanzhou University, Lanzhou 730000, China;2.Unit of 92543 of People’s Liberation Army, Jiyuan 454650, China). China Environmental Science, 2012,32(8):1345~1352

    Adopting observational data of multi angle absorption photometer (MAAP-5012) and the gas composition mixing ratio data of CO and CO2of Semi-Arid Climate Observatory and Laboratory (SACOL) from September 2009 to June 2010, the characteristics and effects of source of black carbon aerosol and carbon gas in semi-arid areas of Northwest China, and the types of emissions which have impact on the concentration of black carbon were analyzed. The result indicated that average concentrations of BC, CO and CO2were 1.75μg/m3, 601.71×10-9and 387.78×10-6, respectively. Hybrid Single Particle Lagrangian Integrated model (HYSPLIT-4) was used to distinguish the flow areas from eastern and western of the observation station. When the airstreams came from the east, the concentrations of BC, CO, CO2were 1.38μg/m3, 462.79×10-9and 383.03×10-6, respectively. When the airstreams came from the west, the concentrations of BC, CO, CO2are 2.2μg/m3, 768.38×10-9and 393.47×10-6, respectively. Taking the clustering analysis on the origin places of the airstreams in the height of 500m, 1500m and 3000m, get the results that when the airstreams of different height flow respective from the Middle East, Central Asia, and Europe, the concentrations of BC, CO, CO2were higher and the values of △B C/△C O, △ CO/△C O2were larger, indicating the efficiency of fuel combustion was lower; when the airstreams transmits from Central and Northern China, the concentrations of BC, CO, CO2were lower, and the values of △B C/△ C O,△ CO/△ CO2were smaller, indicating the efficiency of fuel combustion is higher.

    semi-arid areas;black carbon aerosol;backward trajectory;clustering analysis

    2011-11-21

    國家自然科學基金資助項目(41075104),國家“973”項目(2010CB428604)

    * 責任作者, 教授, zhanglei@lzu.edu.cn

    X131.1, P421

    A

    1000-6923(2012)08-1345-08

    致謝:本文使用了蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)的觀測資料及NOAA提供的HYSPLIT-4模式,謹致謝忱.

    陳 霖(1988-),男,安徽六安人,蘭州大學大氣科學學院碩士研究生,主要從事大氣物理與大氣環(huán)境研究.

    猜你喜歡
    觀測站氣溶膠氣流
    GPS導航對抗數(shù)據(jù)質(zhì)量特征實例分析
    四川省甘孜州:航拍四川稻城高海拔宇宙線觀測站
    科學導報(2022年39期)2022-07-04 23:15:46
    氣流的威力
    科學大眾(2021年21期)2022-01-18 05:53:46
    氣溶膠傳播之謎
    氣溶膠中210Po測定的不確定度評定
    去中心化時差頻差直接定位方法
    航空學報(2017年5期)2017-11-20 03:44:46
    四川盆地秋季氣溶膠與云的相關分析
    固體運載火箭變軌發(fā)動機噴管氣流分離研究
    飛片下的空氣形成的“超強高速氣流刀”
    基于停留時間分布的氣流床氣化爐通用網(wǎng)絡模型
    丝袜喷水一区| 晚上一个人看的免费电影| 大又大粗又爽又黄少妇毛片口| 五月开心婷婷网| 亚洲av国产av综合av卡| 亚洲精品乱久久久久久| 成人无遮挡网站| 18禁动态无遮挡网站| 国产又色又爽无遮挡免| 午夜福利网站1000一区二区三区| 各种免费的搞黄视频| 女人十人毛片免费观看3o分钟| 国产老妇伦熟女老妇高清| 乱码一卡2卡4卡精品| 免费大片黄手机在线观看| 亚洲精品,欧美精品| 国产一区二区亚洲精品在线观看| 五月天丁香电影| 国产色爽女视频免费观看| 网址你懂的国产日韩在线| 日韩欧美精品v在线| 国产视频首页在线观看| 久久精品夜色国产| 久久久精品免费免费高清| 两个人的视频大全免费| 五月玫瑰六月丁香| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 成年版毛片免费区| 亚洲欧美成人综合另类久久久| 日本爱情动作片www.在线观看| 久久99热6这里只有精品| 国产淫语在线视频| 亚洲精品aⅴ在线观看| 亚洲av免费在线观看| 欧美日韩国产mv在线观看视频 | 内地一区二区视频在线| 伦理电影大哥的女人| 国产 一区精品| 男女国产视频网站| 99热6这里只有精品| 国产在线男女| 国产爱豆传媒在线观看| 亚洲国产欧美人成| 亚洲人与动物交配视频| a级一级毛片免费在线观看| 欧美精品一区二区大全| 国产一区二区亚洲精品在线观看| 最近的中文字幕免费完整| 欧美潮喷喷水| 22中文网久久字幕| 亚洲精品国产成人久久av| 乱系列少妇在线播放| 国产探花在线观看一区二区| 交换朋友夫妻互换小说| 大香蕉久久网| 日本wwww免费看| 国产精品久久久久久av不卡| 亚洲图色成人| 国产av码专区亚洲av| 高清日韩中文字幕在线| 亚洲四区av| 国产精品麻豆人妻色哟哟久久| 欧美日韩综合久久久久久| 亚洲欧美日韩无卡精品| 日本色播在线视频| 欧美激情国产日韩精品一区| 国产精品一区二区三区四区免费观看| 黄色日韩在线| 国产精品一区二区三区四区免费观看| 肉色欧美久久久久久久蜜桃 | 国产成人a∨麻豆精品| 国产成人a∨麻豆精品| 精品久久久精品久久久| 一本一本综合久久| 麻豆乱淫一区二区| 亚洲不卡免费看| 免费人成在线观看视频色| 十八禁网站网址无遮挡 | 久久精品夜色国产| 高清在线视频一区二区三区| 天天一区二区日本电影三级| 成人毛片a级毛片在线播放| 午夜激情久久久久久久| 国产淫片久久久久久久久| 欧美极品一区二区三区四区| 高清欧美精品videossex| 亚洲av成人精品一二三区| 日本色播在线视频| 在线观看国产h片| 男男h啪啪无遮挡| 久久午夜福利片| 在线 av 中文字幕| 国产精品一区www在线观看| 亚洲成人一二三区av| 一级黄片播放器| 直男gayav资源| 国产日韩欧美在线精品| 欧美高清成人免费视频www| 十八禁网站网址无遮挡 | 亚洲最大成人手机在线| 国产乱来视频区| 国产 一区 欧美 日韩| 夜夜爽夜夜爽视频| 免费看光身美女| 永久网站在线| 欧美少妇被猛烈插入视频| 青春草国产在线视频| 一个人看视频在线观看www免费| 26uuu在线亚洲综合色| 毛片一级片免费看久久久久| 在线观看一区二区三区激情| 亚洲精品视频女| 青春草亚洲视频在线观看| 亚洲成人一二三区av| 少妇高潮的动态图| 一本色道久久久久久精品综合| 男人和女人高潮做爰伦理| 亚洲av在线观看美女高潮| 最近最新中文字幕大全电影3| 女人十人毛片免费观看3o分钟| 亚洲va在线va天堂va国产| 亚洲成人一二三区av| 男人添女人高潮全过程视频| 国产高清不卡午夜福利| 自拍偷自拍亚洲精品老妇| 欧美日韩精品成人综合77777| 99久久精品热视频| 视频区图区小说| 亚洲精品国产色婷婷电影| 日本一二三区视频观看| 国产精品国产三级国产专区5o| 免费黄网站久久成人精品| 男女那种视频在线观看| 岛国毛片在线播放| 两个人的视频大全免费| 免费观看a级毛片全部| 3wmmmm亚洲av在线观看| 丰满乱子伦码专区| 国产在视频线精品| 内地一区二区视频在线| 纵有疾风起免费观看全集完整版| 国产亚洲最大av| 一级av片app| 看十八女毛片水多多多| 97超视频在线观看视频| 亚洲精品成人av观看孕妇| av在线亚洲专区| 搡女人真爽免费视频火全软件| 婷婷色综合www| 久久这里有精品视频免费| 一区二区三区精品91| 女人十人毛片免费观看3o分钟| 涩涩av久久男人的天堂| 久久精品人妻少妇| 一个人观看的视频www高清免费观看| 国产日韩欧美亚洲二区| 男女边摸边吃奶| 久久99热这里只频精品6学生| 别揉我奶头 嗯啊视频| 国产精品一区二区性色av| 听说在线观看完整版免费高清| 狂野欧美白嫩少妇大欣赏| 久久国内精品自在自线图片| 久久久久久久久大av| av又黄又爽大尺度在线免费看| 午夜免费鲁丝| 一本一本综合久久| 高清毛片免费看| 亚洲图色成人| 麻豆久久精品国产亚洲av| 国产精品偷伦视频观看了| 五月玫瑰六月丁香| 亚洲国产精品成人综合色| 国产人妻一区二区三区在| 啦啦啦啦在线视频资源| 丝袜喷水一区| 夜夜看夜夜爽夜夜摸| 亚洲精品成人久久久久久| 国产精品精品国产色婷婷| 嫩草影院精品99| 日日撸夜夜添| 啦啦啦中文免费视频观看日本| 欧美日韩国产mv在线观看视频 | 一本久久精品| 蜜桃亚洲精品一区二区三区| 日日啪夜夜爽| 欧美+日韩+精品| 亚洲欧美精品自产自拍| av网站免费在线观看视频| 日韩大片免费观看网站| 亚洲aⅴ乱码一区二区在线播放| 特大巨黑吊av在线直播| 精品国产一区二区三区久久久樱花 | 亚洲,欧美,日韩| 七月丁香在线播放| 成人鲁丝片一二三区免费| 一级毛片黄色毛片免费观看视频| 在线观看美女被高潮喷水网站| 高清在线视频一区二区三区| av播播在线观看一区| 最近手机中文字幕大全| 欧美激情国产日韩精品一区| 大码成人一级视频| 如何舔出高潮| 中文天堂在线官网| 美女内射精品一级片tv| 欧美老熟妇乱子伦牲交| 一本色道久久久久久精品综合| 成年免费大片在线观看| 国产精品一区二区性色av| 嫩草影院新地址| 午夜精品一区二区三区免费看| 日本色播在线视频| 亚洲欧美一区二区三区国产| 国产熟女欧美一区二区| 亚洲精品久久午夜乱码| 成人一区二区视频在线观看| 少妇被粗大猛烈的视频| 99热全是精品| 国内精品美女久久久久久| 日韩成人av中文字幕在线观看| 真实男女啪啪啪动态图| av国产免费在线观看| 在线免费十八禁| 九草在线视频观看| 在线观看国产h片| 老女人水多毛片| 日本与韩国留学比较| 一级毛片我不卡| 中文资源天堂在线| 少妇的逼好多水| 欧美成人午夜免费资源| 日本与韩国留学比较| 2021少妇久久久久久久久久久| 久久精品国产亚洲av涩爱| 国产成人精品婷婷| 日日啪夜夜爽| 3wmmmm亚洲av在线观看| 国产乱来视频区| 精品熟女少妇av免费看| 亚洲精品影视一区二区三区av| 国产午夜精品一二区理论片| 国产欧美另类精品又又久久亚洲欧美| 中文字幕亚洲精品专区| 国产精品偷伦视频观看了| 国产亚洲精品久久久com| 国产欧美亚洲国产| 一级二级三级毛片免费看| 精品久久久噜噜| 婷婷色综合大香蕉| 一级黄片播放器| 免费观看av网站的网址| 色播亚洲综合网| 免费看不卡的av| 亚洲一级一片aⅴ在线观看| 国产一级毛片在线| 久久精品人妻少妇| 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 精品一区在线观看国产| 亚洲在线观看片| 久久精品久久久久久噜噜老黄| 亚洲精品成人av观看孕妇| 肉色欧美久久久久久久蜜桃 | 日本-黄色视频高清免费观看| 嫩草影院新地址| 国产在线男女| 亚洲内射少妇av| 伊人久久国产一区二区| 久久久久精品久久久久真实原创| 韩国高清视频一区二区三区| 日韩av在线免费看完整版不卡| 春色校园在线视频观看| 国产男女超爽视频在线观看| av国产精品久久久久影院| 人人妻人人看人人澡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 伦精品一区二区三区| 夫妻午夜视频| 日韩视频在线欧美| 久久久久久伊人网av| 欧美日韩亚洲高清精品| 国产一区二区三区av在线| 黄色视频在线播放观看不卡| 欧美xxⅹ黑人| 婷婷色综合www| 国产白丝娇喘喷水9色精品| 亚洲av一区综合| 欧美变态另类bdsm刘玥| 色哟哟·www| 久久久a久久爽久久v久久| 日韩中字成人| 男人狂女人下面高潮的视频| 国产精品成人在线| 99热6这里只有精品| 人妻夜夜爽99麻豆av| 亚洲婷婷狠狠爱综合网| 伦理电影大哥的女人| 久久综合国产亚洲精品| 亚洲国产精品999| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 成人国产麻豆网| 国产精品久久久久久精品古装| 久久精品人妻少妇| 国产黄色视频一区二区在线观看| 1000部很黄的大片| 欧美成人一区二区免费高清观看| 日韩一区二区视频免费看| 岛国毛片在线播放| 91精品一卡2卡3卡4卡| 亚洲自拍偷在线| 嫩草影院新地址| 尤物成人国产欧美一区二区三区| 亚洲av日韩在线播放| 欧美zozozo另类| 国产日韩欧美在线精品| 久久精品人妻少妇| 啦啦啦在线观看免费高清www| 精品人妻熟女av久视频| 国产色婷婷99| 你懂的网址亚洲精品在线观看| 亚洲内射少妇av| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频 | 国产成人a∨麻豆精品| 一级毛片aaaaaa免费看小| 亚洲av中文av极速乱| 久久韩国三级中文字幕| 日本一本二区三区精品| av播播在线观看一区| 午夜老司机福利剧场| 91精品一卡2卡3卡4卡| 亚洲综合色惰| 国产视频内射| 99视频精品全部免费 在线| 欧美性猛交╳xxx乱大交人| 晚上一个人看的免费电影| 一个人看视频在线观看www免费| 亚洲,欧美,日韩| 最近中文字幕2019免费版| 日本爱情动作片www.在线观看| 一级av片app| 大陆偷拍与自拍| 免费观看性生交大片5| 精品久久久久久久人妻蜜臀av| 国产精品av视频在线免费观看| 韩国av在线不卡| 亚洲色图av天堂| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 亚洲自拍偷在线| 亚洲成人久久爱视频| 大香蕉久久网| 国国产精品蜜臀av免费| 成人鲁丝片一二三区免费| 女人久久www免费人成看片| 777米奇影视久久| 欧美3d第一页| 久久久久久久久久人人人人人人| 国产一区亚洲一区在线观看| 国产69精品久久久久777片| 成人免费观看视频高清| 国产精品国产三级国产av玫瑰| 热99国产精品久久久久久7| 人妻一区二区av| 99热网站在线观看| 亚洲欧美中文字幕日韩二区| 特级一级黄色大片| 免费av毛片视频| 国产精品久久久久久精品电影小说 | 精品久久久精品久久久| 天堂中文最新版在线下载 | 国产精品.久久久| 亚洲,一卡二卡三卡| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 亚洲电影在线观看av| 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 亚洲欧美中文字幕日韩二区| 日日啪夜夜撸| 免费看日本二区| 日本色播在线视频| 别揉我奶头 嗯啊视频| 日韩制服骚丝袜av| 亚洲av成人精品一区久久| 大码成人一级视频| a级一级毛片免费在线观看| 男女国产视频网站| 欧美激情国产日韩精品一区| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美中文字幕日韩二区| 午夜福利在线在线| 香蕉精品网在线| 美女主播在线视频| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 在线观看一区二区三区激情| 夫妻午夜视频| 亚洲精品一二三| av在线天堂中文字幕| 内射极品少妇av片p| 久久6这里有精品| 伦理电影大哥的女人| 一区二区三区免费毛片| 亚洲欧美一区二区三区黑人 | 91精品伊人久久大香线蕉| 免费av不卡在线播放| 亚洲va在线va天堂va国产| 欧美三级亚洲精品| 最近2019中文字幕mv第一页| 亚洲欧美日韩无卡精品| 直男gayav资源| 18禁裸乳无遮挡动漫免费视频 | 久久人人爽av亚洲精品天堂 | 国产一区二区在线观看日韩| 高清在线视频一区二区三区| 国产男人的电影天堂91| av专区在线播放| av一本久久久久| 国产伦在线观看视频一区| 久久99热这里只频精品6学生| 亚洲色图av天堂| 插阴视频在线观看视频| 久久热精品热| 干丝袜人妻中文字幕| 久久精品国产自在天天线| 日日啪夜夜爽| 有码 亚洲区| 欧美日本视频| 欧美激情在线99| 色网站视频免费| 好男人在线观看高清免费视频| 中文字幕av成人在线电影| 日本黄色片子视频| 九九久久精品国产亚洲av麻豆| 国产精品麻豆人妻色哟哟久久| 日本猛色少妇xxxxx猛交久久| 日韩在线高清观看一区二区三区| 精品少妇久久久久久888优播| 成年版毛片免费区| 在现免费观看毛片| 蜜臀久久99精品久久宅男| 18禁动态无遮挡网站| 亚洲av.av天堂| 精品一区二区三卡| 欧美一区二区亚洲| 夫妻性生交免费视频一级片| av女优亚洲男人天堂| 免费看a级黄色片| 国产91av在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 国产成人午夜福利电影在线观看| 欧美激情在线99| 高清午夜精品一区二区三区| 水蜜桃什么品种好| 欧美日韩在线观看h| 亚洲人成网站在线播| 日韩大片免费观看网站| 国产又色又爽无遮挡免| 少妇猛男粗大的猛烈进出视频 | 成年免费大片在线观看| 99久久精品热视频| 秋霞伦理黄片| 91在线精品国自产拍蜜月| 三级男女做爰猛烈吃奶摸视频| 亚洲无线观看免费| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 亚洲欧美日韩另类电影网站 | 18禁裸乳无遮挡动漫免费视频 | 在线观看免费高清a一片| 欧美三级亚洲精品| 天堂俺去俺来也www色官网| 亚洲欧洲国产日韩| 日韩视频在线欧美| 午夜免费鲁丝| 欧美日韩一区二区视频在线观看视频在线 | 久久精品久久精品一区二区三区| 全区人妻精品视频| 99热全是精品| 99久久精品热视频| 亚洲丝袜综合中文字幕| 久久精品久久精品一区二区三区| 欧美三级亚洲精品| 成年免费大片在线观看| 美女国产视频在线观看| 小蜜桃在线观看免费完整版高清| 亚洲国产精品成人综合色| 永久网站在线| 干丝袜人妻中文字幕| 亚洲av男天堂| 国产精品国产av在线观看| 久久6这里有精品| 寂寞人妻少妇视频99o| 久久99热6这里只有精品| 九九久久精品国产亚洲av麻豆| 国产淫语在线视频| 亚洲色图综合在线观看| 三级国产精品欧美在线观看| 日本黄色片子视频| 最近最新中文字幕免费大全7| 精品久久久久久久人妻蜜臀av| 久久6这里有精品| 大话2 男鬼变身卡| 91精品国产九色| 成人国产av品久久久| 免费不卡的大黄色大毛片视频在线观看| 国产伦精品一区二区三区四那| 久久久久久久精品精品| 国产成人精品久久久久久| 欧美3d第一页| 日韩三级伦理在线观看| 特级一级黄色大片| 国产精品熟女久久久久浪| 亚洲成人精品中文字幕电影| 老师上课跳d突然被开到最大视频| 91久久精品国产一区二区三区| 看十八女毛片水多多多| 日韩 亚洲 欧美在线| 美女主播在线视频| 久久精品熟女亚洲av麻豆精品| 国产综合精华液| 免费观看在线日韩| 亚洲精品中文字幕在线视频 | 国内精品美女久久久久久| 女人被狂操c到高潮| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品一二三| 午夜老司机福利剧场| 汤姆久久久久久久影院中文字幕| 久久精品国产a三级三级三级| 夫妻午夜视频| 人人妻人人澡人人爽人人夜夜| 精品人妻熟女av久视频| 国产精品av视频在线免费观看| 成人午夜精彩视频在线观看| 69av精品久久久久久| 国产成人精品久久久久久| 国产精品伦人一区二区| 精品午夜福利在线看| av专区在线播放| 舔av片在线| 看十八女毛片水多多多| 午夜激情福利司机影院| 久久精品久久久久久噜噜老黄| 国产乱人偷精品视频| 国产成年人精品一区二区| 国产伦在线观看视频一区| 国产亚洲精品久久久com| 亚洲熟女精品中文字幕| 国产色婷婷99| 卡戴珊不雅视频在线播放| 久久久亚洲精品成人影院| 丝瓜视频免费看黄片| 国产视频首页在线观看| 人体艺术视频欧美日本| 美女高潮的动态| 国产一区亚洲一区在线观看| 久久鲁丝午夜福利片| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 国产精品人妻久久久久久| 久久久久九九精品影院| 成人黄色视频免费在线看| 日韩视频在线欧美| 亚洲国产精品成人综合色| 亚洲精品视频女| 国产亚洲av片在线观看秒播厂| 青青草视频在线视频观看| 大话2 男鬼变身卡| 亚洲人成网站在线播| 国产视频首页在线观看| 国产精品久久久久久久电影| 嘟嘟电影网在线观看| 色哟哟·www| 在线a可以看的网站| 美女高潮的动态| 一级片'在线观看视频| 国产精品一区二区三区四区免费观看| 最近的中文字幕免费完整| av免费在线看不卡| 亚洲欧美日韩东京热| 亚洲人与动物交配视频| 色吧在线观看| 黄色配什么色好看| 久久精品夜色国产| 丰满人妻一区二区三区视频av| 天天一区二区日本电影三级| 日韩 亚洲 欧美在线| 国产精品一及| 亚洲国产日韩一区二区| 在线免费十八禁| 免费av毛片视频| 久久精品国产亚洲av涩爱| 国产成人aa在线观看| 亚洲国产精品专区欧美| 国产探花极品一区二区| 狂野欧美激情性bbbbbb| 日日撸夜夜添| 久久久久国产网址| 午夜激情福利司机影院| 国产成年人精品一区二区| 成年版毛片免费区| 建设人人有责人人尽责人人享有的 | 精品久久久久久久久亚洲| 赤兔流量卡办理| 国内精品宾馆在线| 亚洲美女视频黄频| 91狼人影院| 麻豆成人午夜福利视频|