• <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级毛片在线播放| 成年女人看的毛片在线观看| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 亚洲人与动物交配视频| 成人性生交大片免费视频hd| 国产午夜精品论理片| av天堂中文字幕网| 亚洲va日本ⅴa欧美va伊人久久| 成人欧美大片| 亚洲精品国产成人久久av| 亚洲性夜色夜夜综合| 偷拍熟女少妇极品色| 亚洲五月天丁香| 91久久精品国产一区二区成人| 男女之事视频高清在线观看| 亚洲狠狠婷婷综合久久图片| 我要搜黄色片| 最后的刺客免费高清国语| 无人区码免费观看不卡| 直男gayav资源| 欧美日本视频| 小说图片视频综合网站| 亚洲精品日韩av片在线观看| 丰满人妻一区二区三区视频av| 一个人免费在线观看电影| 麻豆国产av国片精品| 免费黄网站久久成人精品| 中文字幕av在线有码专区| 国产色婷婷99| 成人高潮视频无遮挡免费网站| 日本-黄色视频高清免费观看| 久久热精品热| 一区福利在线观看| 高清日韩中文字幕在线| 亚洲av日韩精品久久久久久密| 亚洲aⅴ乱码一区二区在线播放| 毛片一级片免费看久久久久 | 天美传媒精品一区二区| 精品久久久久久久末码| 国产美女午夜福利| 亚洲人成伊人成综合网2020| 欧美色欧美亚洲另类二区| 在线看三级毛片| 99热6这里只有精品| 两人在一起打扑克的视频| 精品99又大又爽又粗少妇毛片 | 国内揄拍国产精品人妻在线| 国产精品国产三级国产av玫瑰| 国产精品一区二区三区四区久久| 国产黄a三级三级三级人| 国产精品精品国产色婷婷| 啪啪无遮挡十八禁网站| 成年人黄色毛片网站| 中文字幕人妻熟人妻熟丝袜美| av天堂中文字幕网| 国产精品,欧美在线| 美女cb高潮喷水在线观看| av在线老鸭窝| 午夜福利在线观看免费完整高清在 | 我要搜黄色片| 99国产精品一区二区蜜桃av| 男女视频在线观看网站免费| 在线观看午夜福利视频| 国产精品久久久久久av不卡| 免费观看的影片在线观看| 三级毛片av免费| 亚洲欧美精品综合久久99| 中文在线观看免费www的网站| 欧美一级a爱片免费观看看| 日韩欧美在线乱码| 人妻制服诱惑在线中文字幕| 亚洲电影在线观看av| 91av网一区二区| 偷拍熟女少妇极品色| 久久精品国产99精品国产亚洲性色| 色5月婷婷丁香| 日韩欧美精品v在线| 在线观看美女被高潮喷水网站| 国产精品久久久久久久久免| 国产精品国产高清国产av| 日本 av在线| 亚洲五月天丁香| 老女人水多毛片| av.在线天堂| 五月玫瑰六月丁香| 亚洲一区高清亚洲精品| 白带黄色成豆腐渣| 亚洲av免费在线观看| 久久九九热精品免费| 欧美最黄视频在线播放免费| 国产精品福利在线免费观看| 少妇猛男粗大的猛烈进出视频 | 欧美成人免费av一区二区三区| 国产精品乱码一区二三区的特点| 国产蜜桃级精品一区二区三区| 一区二区三区免费毛片| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 91久久精品国产一区二区成人| 99久久成人亚洲精品观看| 国产亚洲精品久久久久久毛片| 免费不卡的大黄色大毛片视频在线观看 | 久久久久九九精品影院| 一区二区三区免费毛片| 一区二区三区高清视频在线| 嫩草影视91久久| 人人妻,人人澡人人爽秒播| 午夜福利欧美成人| 精品久久国产蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 99九九线精品视频在线观看视频| 国内毛片毛片毛片毛片毛片| 亚洲aⅴ乱码一区二区在线播放| 少妇猛男粗大的猛烈进出视频 | 狠狠狠狠99中文字幕| 国产视频内射| a级一级毛片免费在线观看| 一区二区三区激情视频| av国产免费在线观看| 最新中文字幕久久久久| 欧美黑人欧美精品刺激| 最近视频中文字幕2019在线8| 少妇猛男粗大的猛烈进出视频 | 国产在线男女| 又黄又爽又免费观看的视频| av天堂中文字幕网| 亚洲精品在线观看二区| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 夜夜爽天天搞| 免费在线观看成人毛片| 亚洲第一电影网av| 欧美精品啪啪一区二区三区| 高清日韩中文字幕在线| 国产精品一区www在线观看 | 99在线视频只有这里精品首页| 欧美高清成人免费视频www| 亚洲精品成人久久久久久| 久久天躁狠狠躁夜夜2o2o| 啦啦啦观看免费观看视频高清| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月| 亚洲第一电影网av| 两个人的视频大全免费| 一区二区三区四区激情视频 | 老司机福利观看| 九九在线视频观看精品| or卡值多少钱| 国产av一区在线观看免费| 人妻夜夜爽99麻豆av| 亚洲精品粉嫩美女一区| 亚洲av电影不卡..在线观看| 一区二区三区高清视频在线| 日本a在线网址| 亚洲天堂国产精品一区在线| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久久av| 欧美中文日本在线观看视频| 欧美日韩黄片免| 我要看日韩黄色一级片| 51国产日韩欧美| 综合色av麻豆| 在线免费观看的www视频| 观看美女的网站| or卡值多少钱| 99riav亚洲国产免费| 亚洲一区高清亚洲精品| 亚洲国产日韩欧美精品在线观看| 小说图片视频综合网站| 国产真实伦视频高清在线观看 | 很黄的视频免费| 黄色欧美视频在线观看| .国产精品久久| 欧美日本亚洲视频在线播放| 久久中文看片网| 欧美黑人巨大hd| 97超级碰碰碰精品色视频在线观看| 国产综合懂色| 在线观看66精品国产| 精品国产三级普通话版| 校园人妻丝袜中文字幕| 老司机午夜福利在线观看视频| 国产高清不卡午夜福利| 噜噜噜噜噜久久久久久91| 欧美国产日韩亚洲一区| 欧美+亚洲+日韩+国产| 久久亚洲真实| 午夜免费男女啪啪视频观看 | 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 欧美日本亚洲视频在线播放| 日韩亚洲欧美综合| 免费av不卡在线播放| 国产一区二区在线av高清观看| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 毛片一级片免费看久久久久 | 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区| 大型黄色视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 国产视频内射| 久久久久久久久久成人| 亚洲精品亚洲一区二区| 97超视频在线观看视频| 一夜夜www| 两人在一起打扑克的视频| 国产激情偷乱视频一区二区| 精品一区二区三区视频在线观看免费| 日韩精品有码人妻一区| 欧美一区二区亚洲| 国产美女午夜福利| 国产在线精品亚洲第一网站| 一级黄色大片毛片| 神马国产精品三级电影在线观看| 哪里可以看免费的av片| 欧美xxxx黑人xx丫x性爽| 亚洲天堂国产精品一区在线| 97人妻精品一区二区三区麻豆| 日本欧美国产在线视频| 亚洲专区国产一区二区| 深爱激情五月婷婷| 国产精品野战在线观看| 在线免费观看的www视频| 欧美中文日本在线观看视频| 国产三级中文精品| 欧美xxxx黑人xx丫x性爽| 中文资源天堂在线| 欧美一区二区国产精品久久精品| 直男gayav资源| 亚洲精品在线观看二区| 久99久视频精品免费| 欧美又色又爽又黄视频| 国产av在哪里看| 久久人人爽人人爽人人片va| 久久久久九九精品影院| av在线老鸭窝| 亚洲中文字幕日韩| 免费av毛片视频| 啪啪无遮挡十八禁网站| 色av中文字幕| 成人国产综合亚洲| 午夜福利18| 嫩草影院精品99| 色视频www国产| 国产探花极品一区二区| 村上凉子中文字幕在线| 国内毛片毛片毛片毛片毛片| 精品午夜福利视频在线观看一区| 免费av不卡在线播放| 国产精品一及| 亚洲一区高清亚洲精品| 1000部很黄的大片| 高清毛片免费观看视频网站| 免费无遮挡裸体视频| 狂野欧美激情性xxxx在线观看| 又爽又黄a免费视频| 亚洲欧美日韩无卡精品| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 久久久久九九精品影院| 日本成人三级电影网站| 最近最新免费中文字幕在线| 少妇高潮的动态图| 日本a在线网址| 午夜视频国产福利| 女生性感内裤真人,穿戴方法视频| 又爽又黄a免费视频| 97人妻精品一区二区三区麻豆| 日韩 亚洲 欧美在线| 一级黄色大片毛片| 亚洲精华国产精华精| 欧美激情国产日韩精品一区| 深爱激情五月婷婷| 亚洲国产日韩欧美精品在线观看| 天堂√8在线中文| 内地一区二区视频在线| 少妇丰满av| 在线看三级毛片| 国产高潮美女av| 性色avwww在线观看| 国内精品一区二区在线观看| 69人妻影院| 欧美日本视频| 黄色视频,在线免费观看| 淫妇啪啪啪对白视频| 黄色女人牲交| 色尼玛亚洲综合影院| 国产 一区精品| 国内精品久久久久久久电影| 久久热精品热| 乱系列少妇在线播放| 国产精品久久久久久久电影| 久久久久九九精品影院| 久久香蕉精品热| av视频在线观看入口| 国产午夜精品论理片| 亚洲内射少妇av| 国产在线精品亚洲第一网站| 久久99热这里只有精品18| 91麻豆av在线| av国产免费在线观看| 欧美色欧美亚洲另类二区| 国产高清视频在线播放一区| 成人无遮挡网站| 欧美激情久久久久久爽电影| 69av精品久久久久久| 伦精品一区二区三区| 亚洲精品国产成人久久av| 大又大粗又爽又黄少妇毛片口| 久久中文看片网| 直男gayav资源| 欧美bdsm另类| 一进一出抽搐动态| 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av香蕉五月| 好男人在线观看高清免费视频| 国产精品不卡视频一区二区| 男人舔奶头视频| 日韩欧美免费精品| 国产精品久久电影中文字幕| 精品一区二区三区人妻视频| 亚洲av日韩精品久久久久久密| 国产单亲对白刺激| 小说图片视频综合网站| 色综合站精品国产| 亚洲男人的天堂狠狠| 免费在线观看影片大全网站| 欧美性感艳星| 日韩欧美国产在线观看| 国产三级中文精品| 久久久久久久精品吃奶| 国产精品一区二区免费欧美| 国产男靠女视频免费网站| 中国美女看黄片| 波多野结衣巨乳人妻| 草草在线视频免费看| 亚洲欧美日韩高清专用| 搡老熟女国产l中国老女人| 一级a爱片免费观看的视频| 国产精品久久视频播放| 国产 一区 欧美 日韩| 亚洲人成网站在线播放欧美日韩| 变态另类丝袜制服| 97超级碰碰碰精品色视频在线观看| 亚洲国产精品合色在线| 久久久久久伊人网av| 别揉我奶头~嗯~啊~动态视频| 天天一区二区日本电影三级| 十八禁网站免费在线| 午夜免费激情av| 99精品在免费线老司机午夜| 少妇人妻一区二区三区视频| 俺也久久电影网| 亚洲国产日韩欧美精品在线观看| 九九久久精品国产亚洲av麻豆| 别揉我奶头~嗯~啊~动态视频| 免费无遮挡裸体视频| 日韩强制内射视频| 尤物成人国产欧美一区二区三区| av天堂在线播放| 亚洲欧美激情综合另类| 亚洲欧美日韩无卡精品| 亚洲欧美精品综合久久99| av国产免费在线观看| 精品日产1卡2卡| 少妇丰满av| 蜜桃久久精品国产亚洲av| 日本 av在线| 日韩欧美在线二视频| 国产精品av视频在线免费观看| 亚洲美女搞黄在线观看 | 一个人免费在线观看电影| 内射极品少妇av片p| 观看美女的网站| 伊人久久精品亚洲午夜| 亚洲欧美日韩卡通动漫| 在线a可以看的网站| 99riav亚洲国产免费| 日本黄色视频三级网站网址| 久久午夜亚洲精品久久| 国产乱人视频| 亚洲一区高清亚洲精品| 最新中文字幕久久久久| 国产大屁股一区二区在线视频| 他把我摸到了高潮在线观看| 久久久久久久亚洲中文字幕| 国产成人福利小说| 在线观看舔阴道视频| 超碰av人人做人人爽久久| 两性午夜刺激爽爽歪歪视频在线观看| 夜夜爽天天搞| 日本-黄色视频高清免费观看| 精品欧美国产一区二区三| 国产精品综合久久久久久久免费| 久久精品国产亚洲网站| 亚洲成人中文字幕在线播放| 亚洲久久久久久中文字幕| 搡老岳熟女国产| 欧美色欧美亚洲另类二区| 欧美区成人在线视频| 欧美性感艳星| 国内精品久久久久久久电影| 免费无遮挡裸体视频| 我要看日韩黄色一级片| 欧美潮喷喷水| 丰满的人妻完整版| 男人狂女人下面高潮的视频| 身体一侧抽搐| 亚洲国产精品成人综合色| 国产激情偷乱视频一区二区| 国产精品免费一区二区三区在线| 国产精品久久久久久久久免| 91午夜精品亚洲一区二区三区 | 国产精品久久久久久久电影| 亚洲av熟女| 中文字幕高清在线视频| 国产精品国产三级国产av玫瑰| 欧美3d第一页| 久久热精品热| 国产亚洲精品久久久com| 一个人看的www免费观看视频| 最近视频中文字幕2019在线8| 美女被艹到高潮喷水动态| 丰满的人妻完整版| 亚洲美女搞黄在线观看 | 国内精品宾馆在线| 波多野结衣高清无吗| 欧美中文日本在线观看视频| 日韩欧美一区二区三区在线观看| 亚洲三级黄色毛片| 欧美一区二区亚洲| 国语自产精品视频在线第100页| 中国美女看黄片| 精品一区二区三区人妻视频| 免费在线观看成人毛片| 最新中文字幕久久久久| av在线老鸭窝| 欧美一区二区国产精品久久精品| 两性午夜刺激爽爽歪歪视频在线观看| 一个人观看的视频www高清免费观看| 国产亚洲精品久久久久久毛片| 悠悠久久av| 人妻少妇偷人精品九色| 村上凉子中文字幕在线| 在线播放国产精品三级| 亚洲熟妇熟女久久| 色5月婷婷丁香| 91久久精品国产一区二区成人| 在线天堂最新版资源| 乱系列少妇在线播放| 国产午夜精品久久久久久一区二区三区 | 成年人黄色毛片网站| 91午夜精品亚洲一区二区三区 | 久久久久久久久久黄片| 午夜日韩欧美国产| 91久久精品国产一区二区成人| 成人特级黄色片久久久久久久| 日韩 亚洲 欧美在线| 日韩欧美精品v在线| 久久久久久久久久久丰满 | 高清在线国产一区| 少妇的逼好多水| 国产伦在线观看视频一区| 欧美黑人巨大hd| 精品一区二区三区人妻视频| av女优亚洲男人天堂| 国产高潮美女av| 最近视频中文字幕2019在线8| 女人十人毛片免费观看3o分钟| 美女被艹到高潮喷水动态| 波多野结衣高清作品| 免费在线观看日本一区| 九色成人免费人妻av| 色综合婷婷激情| 观看美女的网站| 白带黄色成豆腐渣| 日韩高清综合在线| 国产精品一区二区性色av| 一进一出抽搐动态| 人人妻人人澡欧美一区二区| 桃色一区二区三区在线观看| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 日韩人妻高清精品专区| 免费观看的影片在线观看| h日本视频在线播放| 日本a在线网址| 欧美激情国产日韩精品一区| 欧美性感艳星| 久久久久免费精品人妻一区二区| 淫秽高清视频在线观看| 午夜免费男女啪啪视频观看 | 国产精品一区二区三区四区久久| 天堂网av新在线| 国产高清视频在线观看网站| 日本一本二区三区精品| 亚洲自偷自拍三级| 男女那种视频在线观看| 日韩欧美三级三区| 日本一本二区三区精品| 亚洲精品色激情综合| 国产色婷婷99| 亚洲内射少妇av| 少妇的逼水好多| 日韩欧美免费精品| 动漫黄色视频在线观看| 欧美中文日本在线观看视频| 国产成年人精品一区二区| 成人永久免费在线观看视频| 18禁裸乳无遮挡免费网站照片| 日韩欧美国产一区二区入口| 国产av一区在线观看免费| 91久久精品国产一区二区三区| 国内精品一区二区在线观看| 别揉我奶头~嗯~啊~动态视频| 日韩中字成人| 美女被艹到高潮喷水动态| 女人十人毛片免费观看3o分钟| 亚洲成人久久性| 亚洲av不卡在线观看| 五月玫瑰六月丁香| 亚洲最大成人中文| 99精品久久久久人妻精品| 国产不卡一卡二| 国产欧美日韩精品一区二区| 久久久久久大精品| 中文字幕av在线有码专区| 国产一区二区三区在线臀色熟女| 中文资源天堂在线| 制服丝袜大香蕉在线| 91狼人影院| а√天堂www在线а√下载| 日本熟妇午夜| 狠狠狠狠99中文字幕| 色综合色国产| www.www免费av| 一区二区三区激情视频| 久9热在线精品视频| 成人av一区二区三区在线看| 别揉我奶头~嗯~啊~动态视频| 国产乱人视频| 国产亚洲欧美98| 精品久久久久久久久亚洲 | 男女做爰动态图高潮gif福利片| 午夜久久久久精精品| 亚洲国产欧洲综合997久久,| 美女被艹到高潮喷水动态| 我要看日韩黄色一级片| 不卡视频在线观看欧美| 嫩草影院入口| 动漫黄色视频在线观看| 国产人妻一区二区三区在| 午夜影院日韩av| 女的被弄到高潮叫床怎么办 | 成熟少妇高潮喷水视频| 男女下面进入的视频免费午夜| 黄色丝袜av网址大全| 国产男人的电影天堂91| 99久久精品一区二区三区| 噜噜噜噜噜久久久久久91| 搞女人的毛片| 亚洲精品影视一区二区三区av| 日韩中文字幕欧美一区二区| 麻豆成人午夜福利视频| 国产男靠女视频免费网站| 精品日产1卡2卡| 国产精品一及| 国产精品野战在线观看| 日韩欧美免费精品| 亚洲av五月六月丁香网| 国产精华一区二区三区| 国产乱人视频| 久久精品国产99精品国产亚洲性色| 热99在线观看视频| 亚洲在线观看片| 91在线精品国自产拍蜜月| 大型黄色视频在线免费观看| 一本精品99久久精品77| 亚洲七黄色美女视频| 国产爱豆传媒在线观看| 精品人妻一区二区三区麻豆 | 日本一二三区视频观看| 精品人妻1区二区| 精品一区二区三区av网在线观看| 久久久久久久精品吃奶| 我的老师免费观看完整版| 午夜激情福利司机影院| 精品国产三级普通话版| 国内精品久久久久精免费| 国产一区二区亚洲精品在线观看| 搡老熟女国产l中国老女人| 一个人观看的视频www高清免费观看| 中文亚洲av片在线观看爽| 美女高潮喷水抽搐中文字幕| 欧美一级a爱片免费观看看| 噜噜噜噜噜久久久久久91| 少妇熟女aⅴ在线视频| 18禁裸乳无遮挡免费网站照片| 国内少妇人妻偷人精品xxx网站| 色视频www国产| 婷婷六月久久综合丁香| 精品一区二区三区av网在线观看| or卡值多少钱| 免费无遮挡裸体视频| 亚洲欧美日韩高清专用| 亚洲乱码一区二区免费版| 欧美高清性xxxxhd video| 99在线人妻在线中文字幕| 欧美绝顶高潮抽搐喷水|