• <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)絡模型
    久久久久久国产a免费观看| 久久亚洲国产成人精品v| 欧美成人a在线观看| 亚洲av免费在线观看| 日韩人妻高清精品专区| 日本一二三区视频观看| 精品久久久久久久久亚洲| 1000部很黄的大片| 在线a可以看的网站| 伦理电影大哥的女人| 免费av不卡在线播放| 日本免费在线观看一区| 婷婷六月久久综合丁香| 欧美成人一区二区免费高清观看| 秋霞在线观看毛片| 国产欧美日韩精品一区二区| 天天躁夜夜躁狠狠久久av| 日日啪夜夜撸| 亚洲精品aⅴ在线观看| 欧美三级亚洲精品| 色尼玛亚洲综合影院| 久久人妻av系列| АⅤ资源中文在线天堂| 久久久午夜欧美精品| 九草在线视频观看| 精品久久久久久久久久久久久| 久久综合国产亚洲精品| 亚洲av日韩在线播放| 精品久久久噜噜| 国产麻豆成人av免费视频| 一级毛片电影观看 | 国产毛片a区久久久久| 少妇熟女aⅴ在线视频| 熟妇人妻久久中文字幕3abv| 内射极品少妇av片p| 久久精品人妻少妇| 中国国产av一级| 国产精品一区www在线观看| 国产成人a区在线观看| 波野结衣二区三区在线| 国产精品一及| 日韩亚洲欧美综合| 一级毛片aaaaaa免费看小| 色网站视频免费| 久久精品国产自在天天线| 天美传媒精品一区二区| 成年av动漫网址| 国产探花极品一区二区| 看非洲黑人一级黄片| 乱人视频在线观看| 免费观看的影片在线观看| 欧美zozozo另类| 别揉我奶头 嗯啊视频| 啦啦啦啦在线视频资源| 建设人人有责人人尽责人人享有的 | 欧美精品一区二区大全| 国产综合懂色| 99热全是精品| 久久久久久久亚洲中文字幕| 亚洲成av人片在线播放无| 久久精品91蜜桃| 最近中文字幕2019免费版| 亚洲精品aⅴ在线观看| 国产老妇女一区| 国产黄色视频一区二区在线观看 | 秋霞伦理黄片| 免费无遮挡裸体视频| 自拍偷自拍亚洲精品老妇| 日本猛色少妇xxxxx猛交久久| 国产成人a∨麻豆精品| 精品欧美国产一区二区三| 国产91av在线免费观看| 在线播放无遮挡| 波野结衣二区三区在线| 天天躁夜夜躁狠狠久久av| 国产成人精品久久久久久| 国产成人aa在线观看| 久99久视频精品免费| 色视频www国产| 欧美一级a爱片免费观看看| 国产精品99久久久久久久久| 内射极品少妇av片p| 最近中文字幕2019免费版| 22中文网久久字幕| 亚洲精品亚洲一区二区| 免费看a级黄色片| 在线免费观看的www视频| 听说在线观看完整版免费高清| 国产成人aa在线观看| 久久久久久久久中文| 日产精品乱码卡一卡2卡三| 中文字幕熟女人妻在线| 一级av片app| 看免费成人av毛片| av免费在线看不卡| 久久精品夜色国产| 男女那种视频在线观看| 日韩,欧美,国产一区二区三区 | 色5月婷婷丁香| 成人亚洲欧美一区二区av| 97人妻精品一区二区三区麻豆| 国产男人的电影天堂91| 九九热线精品视视频播放| 毛片一级片免费看久久久久| 少妇人妻精品综合一区二区| 国产精品永久免费网站| 国产精品久久电影中文字幕| 18禁在线无遮挡免费观看视频| 中文乱码字字幕精品一区二区三区 | 五月玫瑰六月丁香| 国产精品久久久久久精品电影| 舔av片在线| 国产真实乱freesex| 久久草成人影院| 国产免费又黄又爽又色| 免费无遮挡裸体视频| 免费看a级黄色片| 亚洲精品乱久久久久久| 日韩成人av中文字幕在线观看| 男人和女人高潮做爰伦理| 国产v大片淫在线免费观看| 99热精品在线国产| 免费无遮挡裸体视频| 汤姆久久久久久久影院中文字幕 | 日日啪夜夜撸| 亚洲av免费在线观看| 女人十人毛片免费观看3o分钟| 欧美xxxx性猛交bbbb| 又爽又黄a免费视频| 国产成人一区二区在线| 身体一侧抽搐| 国产精品久久久久久精品电影| 亚洲国产精品成人综合色| 久久精品夜色国产| 久久久a久久爽久久v久久| 国产成人freesex在线| 午夜福利高清视频| 草草在线视频免费看| 听说在线观看完整版免费高清| 99久国产av精品国产电影| 特大巨黑吊av在线直播| 看片在线看免费视频| 小说图片视频综合网站| 亚洲成人中文字幕在线播放| 国产一区二区在线av高清观看| 国产免费视频播放在线视频 | 午夜福利成人在线免费观看| 国产老妇伦熟女老妇高清| 午夜激情欧美在线| 国产综合懂色| 国产精品一区二区三区四区免费观看| 午夜日本视频在线| 精品午夜福利在线看| 男女视频在线观看网站免费| 建设人人有责人人尽责人人享有的 | 18禁裸乳无遮挡免费网站照片| 久久草成人影院| 免费av观看视频| 身体一侧抽搐| 亚洲av.av天堂| 女人十人毛片免费观看3o分钟| 韩国av在线不卡| 日韩视频在线欧美| 国产精品一区二区在线观看99 | 国内精品美女久久久久久| 亚洲av成人精品一二三区| 成人漫画全彩无遮挡| 国产精品美女特级片免费视频播放器| 国产高清三级在线| 国产精品.久久久| 自拍偷自拍亚洲精品老妇| 精品久久国产蜜桃| 免费观看的影片在线观看| 九九爱精品视频在线观看| 国产免费男女视频| 午夜福利高清视频| 亚洲av福利一区| 18禁在线无遮挡免费观看视频| 如何舔出高潮| 免费av观看视频| 国产乱人视频| 天堂√8在线中文| 最近手机中文字幕大全| 高清毛片免费看| 岛国毛片在线播放| 国内揄拍国产精品人妻在线| 成人漫画全彩无遮挡| 精品无人区乱码1区二区| 又黄又爽又刺激的免费视频.| 在线天堂最新版资源| 在线观看一区二区三区| 久久久久网色| 99热精品在线国产| 日韩av不卡免费在线播放| 国产成人精品久久久久久| 午夜福利在线观看免费完整高清在| 午夜精品一区二区三区免费看| 国产精品99久久久久久久久| 春色校园在线视频观看| 韩国av在线不卡| 日韩视频在线欧美| 永久网站在线| 免费不卡的大黄色大毛片视频在线观看 | 九九爱精品视频在线观看| 大香蕉久久网| 色综合亚洲欧美另类图片| 久久99热这里只频精品6学生 | 人妻制服诱惑在线中文字幕| 黄片无遮挡物在线观看| 七月丁香在线播放| 春色校园在线视频观看| 精品人妻一区二区三区麻豆| 亚洲国产高清在线一区二区三| 国产男人的电影天堂91| 大香蕉97超碰在线| 天堂网av新在线| 国产免费一级a男人的天堂| 日本wwww免费看| 亚洲一级一片aⅴ在线观看| 国产av不卡久久| 国产亚洲最大av| 中国国产av一级| 成人午夜高清在线视频| 午夜老司机福利剧场| 久久久国产成人免费| 中文字幕精品亚洲无线码一区| 国产亚洲一区二区精品| 亚洲久久久久久中文字幕| 观看免费一级毛片| 国产伦精品一区二区三区四那| 中国美白少妇内射xxxbb| 白带黄色成豆腐渣| 免费观看在线日韩| 亚洲精品自拍成人| 日本午夜av视频| 99久久中文字幕三级久久日本| 床上黄色一级片| 狂野欧美激情性xxxx在线观看| 搡老妇女老女人老熟妇| 两个人视频免费观看高清| 日韩欧美精品免费久久| 麻豆成人午夜福利视频| 亚洲一区高清亚洲精品| 亚洲人成网站在线播| 国产成人a区在线观看| 亚洲精华国产精华液的使用体验| 亚洲国产精品合色在线| 久久精品91蜜桃| 汤姆久久久久久久影院中文字幕 | 在线播放无遮挡| 国语自产精品视频在线第100页| 国内精品美女久久久久久| 久久99精品国语久久久| 亚洲国产精品专区欧美| 床上黄色一级片| 女人十人毛片免费观看3o分钟| 亚洲最大成人手机在线| 99热全是精品| 亚洲国产精品久久男人天堂| 美女被艹到高潮喷水动态| 天天一区二区日本电影三级| av.在线天堂| 亚洲乱码一区二区免费版| 精品一区二区免费观看| 一边亲一边摸免费视频| 日本av手机在线免费观看| 一本久久精品| 国产视频内射| 亚洲,欧美,日韩| 天天躁夜夜躁狠狠久久av| eeuss影院久久| 搡女人真爽免费视频火全软件| 变态另类丝袜制服| 免费观看性生交大片5| 午夜免费激情av| 国产精品一二三区在线看| 乱码一卡2卡4卡精品| 91精品一卡2卡3卡4卡| 国产老妇女一区| 久久99精品国语久久久| 国产精品国产三级国产av玫瑰| 成人漫画全彩无遮挡| videossex国产| 国产亚洲一区二区精品| 日本一本二区三区精品| 高清视频免费观看一区二区 | 国产淫语在线视频| h日本视频在线播放| 午夜福利高清视频| 3wmmmm亚洲av在线观看| 亚洲av免费在线观看| 最近中文字幕高清免费大全6| 大香蕉97超碰在线| 欧美日本视频| 国内精品美女久久久久久| 色综合站精品国产| 91精品国产九色| 色网站视频免费| 菩萨蛮人人尽说江南好唐韦庄 | 成人漫画全彩无遮挡| 色尼玛亚洲综合影院| 国产成人精品久久久久久| 免费av毛片视频| 亚洲精品国产成人久久av| 亚洲综合精品二区| 国产精品国产三级国产专区5o | 欧美bdsm另类| or卡值多少钱| 亚洲最大成人手机在线| 一卡2卡三卡四卡精品乱码亚洲| 国产淫语在线视频| 一级黄色大片毛片| 成人三级黄色视频| 国产精品日韩av在线免费观看| 精品久久久久久久末码| 欧美3d第一页| 欧美一级a爱片免费观看看| 女人久久www免费人成看片 | 老司机影院成人| 久久精品国产鲁丝片午夜精品| 精品久久久久久久人妻蜜臀av| 你懂的网址亚洲精品在线观看 | 午夜激情欧美在线| 亚洲av成人精品一区久久| 白带黄色成豆腐渣| 中文欧美无线码| 久久精品综合一区二区三区| 最后的刺客免费高清国语| 久久99热6这里只有精品| 日韩视频在线欧美| 国产老妇女一区| 男女下面进入的视频免费午夜| 乱人视频在线观看| 亚洲最大成人手机在线| 久久久久久久久久久丰满| 亚洲欧美日韩东京热| 亚洲精品乱码久久久久久按摩| 在线免费观看的www视频| 青春草亚洲视频在线观看| 亚洲最大成人手机在线| 别揉我奶头 嗯啊视频| 日本熟妇午夜| 国产极品天堂在线| АⅤ资源中文在线天堂| 国产精华一区二区三区| 亚洲不卡免费看| 搞女人的毛片| 亚洲精品日韩在线中文字幕| 国产高清视频在线观看网站| 亚洲国产精品成人综合色| 啦啦啦啦在线视频资源| 秋霞伦理黄片| 夫妻性生交免费视频一级片| av播播在线观看一区| 国国产精品蜜臀av免费| 免费观看的影片在线观看| 天堂影院成人在线观看| 白带黄色成豆腐渣| 国产不卡一卡二| 偷拍熟女少妇极品色| 人妻系列 视频| 亚洲精华国产精华液的使用体验| 久久热精品热| 国产成人精品婷婷| 99久久精品国产国产毛片| 中文在线观看免费www的网站| 国产黄片美女视频| 黄片wwwwww| 国产高清国产精品国产三级 | 久久久国产成人免费| 欧美一级a爱片免费观看看| 成年女人看的毛片在线观看| 久久精品国产自在天天线| av视频在线观看入口| 嘟嘟电影网在线观看| 搞女人的毛片| 97人妻精品一区二区三区麻豆| 校园人妻丝袜中文字幕| 国产69精品久久久久777片| 91精品一卡2卡3卡4卡| 国产精品国产三级专区第一集| 91久久精品国产一区二区三区| 男人和女人高潮做爰伦理| 最近最新中文字幕大全电影3| av天堂中文字幕网| 亚洲国产欧美人成| 精品久久久久久成人av| 中文字幕av成人在线电影| 一边亲一边摸免费视频| 在线免费十八禁| a级一级毛片免费在线观看| 91久久精品电影网| 欧美精品国产亚洲| 午夜爱爱视频在线播放| 欧美变态另类bdsm刘玥| 久久精品人妻少妇| 嫩草影院精品99| 亚洲一区高清亚洲精品| 亚洲在久久综合| 久久精品国产自在天天线| 免费搜索国产男女视频| 久久久久久久久久久丰满| 国产av码专区亚洲av| 美女黄网站色视频| 国产精品嫩草影院av在线观看| 天堂网av新在线| 午夜福利高清视频| 国产精品爽爽va在线观看网站| 青春草国产在线视频| 久久久久久久久久成人| 亚洲成人中文字幕在线播放| 久久6这里有精品| 亚洲国产精品国产精品| 色综合色国产| 欧美性感艳星| 天堂av国产一区二区熟女人妻| 亚洲av.av天堂| 女人被狂操c到高潮| 欧美不卡视频在线免费观看| 黄片wwwwww| 午夜福利在线观看吧| 欧美丝袜亚洲另类| 免费观看在线日韩| 中文精品一卡2卡3卡4更新| 久久久国产成人免费| 国产精品人妻久久久久久| 亚洲乱码一区二区免费版| 亚洲欧美日韩卡通动漫| 级片在线观看| 欧美精品一区二区大全| 99九九线精品视频在线观看视频| 女的被弄到高潮叫床怎么办| 小说图片视频综合网站| 亚洲国产欧洲综合997久久,| 非洲黑人性xxxx精品又粗又长| 两性午夜刺激爽爽歪歪视频在线观看| 少妇丰满av| 一本久久精品| 中文字幕av在线有码专区| 欧美日韩在线观看h| 亚洲婷婷狠狠爱综合网| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 嫩草影院精品99| 高清毛片免费看| 国内精品美女久久久久久| 青春草亚洲视频在线观看| 精品99又大又爽又粗少妇毛片| 日本免费在线观看一区| 色综合站精品国产| 欧美成人精品欧美一级黄| 美女脱内裤让男人舔精品视频| 男女边吃奶边做爰视频| 欧美成人精品欧美一级黄| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 少妇的逼水好多| 日韩精品青青久久久久久| 欧美激情在线99| 亚洲人成网站高清观看| 国产av在哪里看| 好男人视频免费观看在线| 亚洲三级黄色毛片| 少妇人妻一区二区三区视频| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 亚洲国产精品成人久久小说| 青春草视频在线免费观看| 国产一区二区三区av在线| 国产在视频线精品| 亚洲精华国产精华液的使用体验| 免费播放大片免费观看视频在线观看 | 又粗又硬又长又爽又黄的视频| 干丝袜人妻中文字幕| 日本黄色视频三级网站网址| 欧美性猛交╳xxx乱大交人| 国产色爽女视频免费观看| 日韩制服骚丝袜av| 精品久久久久久成人av| 欧美日韩在线观看h| 国产精华一区二区三区| 亚洲性久久影院| 亚洲av免费在线观看| 干丝袜人妻中文字幕| 国产极品精品免费视频能看的| 丰满乱子伦码专区| 国产久久久一区二区三区| 日韩av在线免费看完整版不卡| 男女视频在线观看网站免费| 国产日韩欧美在线精品| 深夜a级毛片| 大香蕉97超碰在线| 国产精品乱码一区二三区的特点| 亚洲熟妇中文字幕五十中出| 女人被狂操c到高潮| 国产成人精品久久久久久| 国产黄片视频在线免费观看| 亚洲三级黄色毛片| 26uuu在线亚洲综合色| 麻豆成人午夜福利视频| 黄色一级大片看看| 国产一区二区在线av高清观看| 成年av动漫网址| 精华霜和精华液先用哪个| 热99re8久久精品国产| 亚洲欧美日韩东京热| 亚洲最大成人av| 99久久人妻综合| 别揉我奶头 嗯啊视频| 国产不卡一卡二| 成年免费大片在线观看| 国产一级毛片七仙女欲春2| 日日摸夜夜添夜夜爱| 久久久久久久久中文| 国产成人免费观看mmmm| 女人被狂操c到高潮| 1000部很黄的大片| 亚洲婷婷狠狠爱综合网| 伦精品一区二区三区| 99在线人妻在线中文字幕| 成人二区视频| 又粗又爽又猛毛片免费看| 18禁裸乳无遮挡免费网站照片| 日韩国内少妇激情av| 精品久久久久久久久av| av黄色大香蕉| 精华霜和精华液先用哪个| 在线免费观看不下载黄p国产| av天堂中文字幕网| 校园人妻丝袜中文字幕| 小蜜桃在线观看免费完整版高清| 美女黄网站色视频| 免费黄网站久久成人精品| 日本五十路高清| 午夜爱爱视频在线播放| av福利片在线观看| 校园人妻丝袜中文字幕| 免费一级毛片在线播放高清视频| 日本wwww免费看| 久久草成人影院| 2022亚洲国产成人精品| 国产真实乱freesex| 熟女电影av网| 嫩草影院入口| 日韩,欧美,国产一区二区三区 | 在线免费十八禁| 人妻夜夜爽99麻豆av| 国产麻豆成人av免费视频| 亚洲欧美日韩卡通动漫| 亚洲成av人片在线播放无| 日本猛色少妇xxxxx猛交久久| 亚洲人成网站在线播| 日韩精品有码人妻一区| 99久久人妻综合| 高清日韩中文字幕在线| 欧美97在线视频| 丰满乱子伦码专区| 91狼人影院| 日日摸夜夜添夜夜爱| 亚洲无线观看免费| 两个人视频免费观看高清| 国产精品一区二区在线观看99 | 男女啪啪激烈高潮av片| 欧美xxxx黑人xx丫x性爽| 国产毛片a区久久久久| 一夜夜www| 欧美激情国产日韩精品一区| 久久久精品欧美日韩精品| 国产激情偷乱视频一区二区| 免费看a级黄色片| 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| 色播亚洲综合网| 亚洲无线观看免费| 听说在线观看完整版免费高清| 少妇的逼水好多| 国产av不卡久久| 啦啦啦观看免费观看视频高清| 国产欧美另类精品又又久久亚洲欧美| 精品久久久久久久久久久久久| 久久这里只有精品中国| 老司机福利观看| 亚洲成人精品中文字幕电影| 99热全是精品| 欧美又色又爽又黄视频| 三级毛片av免费| 午夜福利网站1000一区二区三区| 在现免费观看毛片| 国产一级毛片七仙女欲春2| 99久久精品国产国产毛片| 日本免费在线观看一区| 春色校园在线视频观看| 亚洲电影在线观看av| 大香蕉97超碰在线| 日本免费在线观看一区| 久久韩国三级中文字幕| 久久久久国产网址| 美女xxoo啪啪120秒动态图| 日韩在线高清观看一区二区三区| 国产av不卡久久| a级一级毛片免费在线观看| 久久久精品94久久精品| kizo精华| 偷拍熟女少妇极品色| 欧美激情国产日韩精品一区| 欧美性猛交╳xxx乱大交人| 97超碰精品成人国产| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 国产黄色小视频在线观看| 国内精品一区二区在线观看| 欧美一区二区亚洲| 国产成人freesex在线| 成人午夜高清在线视频| 久久久久久久久久久免费av|