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

    中國3個AERONET站點氣溶膠微物理特性分析及比較

    2014-05-09 08:22:50張志薇王宏斌田鵬飛梁捷寧王式功
    中國環(huán)境科學 2014年8期
    關(guān)鍵詞:香河氣溶膠太湖

    張志薇,王宏斌,張 鐳,胡 蝶,田鵬飛,梁捷寧,王式功

    (1.蘭州大學半干旱氣候變化教育部重點實驗室,大氣科學學院,甘肅 蘭州730000;2.江蘇省氣象服務中心,江蘇 南京210008;3.江蘇省氣象科學研究所,江蘇 南京210009;4.中國氣象局蘭州干旱氣象研究所,甘肅 蘭州730020)

    中國3個AERONET站點氣溶膠微物理特性分析及比較

    張志薇1,2,王宏斌3,張 鐳1*,胡 蝶4,田鵬飛1,梁捷寧1,王式功1

    (1.蘭州大學半干旱氣候變化教育部重點實驗室,大氣科學學院,甘肅 蘭州730000;2.江蘇省氣象服務中心,江蘇 南京210008;3.江蘇省氣象科學研究所,江蘇 南京210009;4.中國氣象局蘭州干旱氣象研究所,甘肅 蘭州730020)

    選取中國地區(qū)區(qū)域代表性較強且觀測時間序列較長的3個AERONET站點(SACOL、香河和太湖),分析了其氣溶膠微物理參數(shù)特征. 香河和太湖多年平均氣溶膠光學厚度(AOD)分別為0.67±0.66和0.72±0.44,是SACOL AOD平均值(0.38±0.27)的近2倍,且AOD變化范圍較大. SACOL春冬季AOD較大,夏秋季AOD較小;而香河和太湖夏春季較大,秋冬季較小.結(jié)合尺度分布、體積濃度等參數(shù)特征說明沙塵是SACOL春季最主要的氣溶膠類型,香河春季受沙塵的影響也較嚴重,而太湖受沙塵影響的頻率較香河要小的多;香河和太湖AOD最大值出現(xiàn)月份與細模態(tài)粒子體積濃度最大值出現(xiàn)月份一致,是由于細模態(tài)粒子的消光效率是粗模態(tài)粒子的3~4倍.細模態(tài)體積比(Vf/Vt)的年變化趨勢與?ngstr?m波長指數(shù)(α)的年變化趨勢相似,Vf/Vt和α均可以用來分析粒子尺度大小的年變化特征.但α<0.75時,3個站點Vf/Vt均小于0.4,以粗模態(tài)粒子為主; α>1.7時,Vf/Vt大于0.6,以細模態(tài)粒子為主;而0.75<α<1.7時,也是α出現(xiàn)概率最高的區(qū)間,3個站點Vf/Vt均在0.1~0.8之間變化.

    氣溶膠光學厚度;體積濃度;細模態(tài)體積比;?ngstr?m波長指數(shù)

    大氣氣溶膠作為地氣系統(tǒng)的重要組成部分,其直接散射和吸收太陽輻射影響氣候系統(tǒng)的輻射平衡(直接效應);同時作為云凝結(jié)核,可以改變云的物理和微物理特征(云反照率、云壽命、云水路徑等)進而改變云的輻射特征,影響太陽能在氣候系統(tǒng)中的分配(間接效應)[1-3];除此,吸收性氣溶膠會將其吸收的太陽輻射能轉(zhuǎn)變?yōu)闊崮苤匦箩尫?影響大氣穩(wěn)定度、地表能量收支,并可能導致云滴蒸發(fā),造成云反照率和云量的減小,進而影響氣候[4-5].除此,大氣氣溶膠對環(huán)境和人體健康有重要影響,如嚴重降低能見度[6-7]、對人體肺功能造成損傷等[8].與溫室氣體不同,氣溶膠氣候和環(huán)境效應的評估非常困難,因為氣溶膠是由不同形狀、尺度分布、化學組成和光學性質(zhì)的物質(zhì)構(gòu)成的,同時它們濃度的時空變化可達幾個數(shù)量級且其時空變化的觀測資料較少[9-12].

    政府間氣候變化專業(yè)委員會(IPCC)第四次科學評估報告[13]對大氣氣溶膠輻射效應理解和評估比之前三次報告準確.研究表明大多數(shù)氣候模式低估了人為氣溶膠導致的負強迫[14].而大氣氣溶膠光學和微物理特性的長期觀測對這些結(jié)果的得到有著至關(guān)重要的作用.長期系統(tǒng)的氣溶膠光學和微物理特性觀測資料不僅是認識區(qū)域氣溶膠特性的重要依據(jù),而且是定量研究氣溶膠環(huán)境和輻射效應的重要基礎(chǔ).目前在全球范圍內(nèi)已建立了多個氣溶膠觀測網(wǎng),如 AERONET[15]、PFR-GAW和SKYNET[16]等.

    目前對中國地區(qū)氣溶膠光學和微物理特性及其區(qū)域差異的認識還有很大的不確定性,從而限制了中國地區(qū)衛(wèi)星資料的應用和對氣溶膠環(huán)境和輻射效應的認識.氣溶膠的光學特性與其尺度分布有很強的依賴性,任何可能影響其尺度分布的因素都有可能影響其光學特性的變化,如相對濕度的變化會引起吸濕性氣溶膠單次散射反照率和不對稱因子的變化[17-18].因此對氣溶膠尺度分布、體積濃度、細模態(tài)粒子體積比和有效半徑等微物理特性的分析有助于增強對其環(huán)境和輻射效應的認識.本文選取中國地區(qū)區(qū)域代表性較強且觀測時間序列較長的3個AERONET站點(蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)、香河和太湖)分析其氣溶膠微物理參數(shù)特征及區(qū)域差異.

    1 站點與數(shù)據(jù)

    AERONET是美國國家宇航局(NASA)和法國國家科學研究中心(CNRS)共同組建的全球數(shù)百個太陽光度計組成的氣溶膠監(jiān)測網(wǎng),其觀測數(shù)據(jù)精度較高,氣溶膠光學厚度(AOD)觀測誤差為0.01~0.02.其對全球和區(qū)域氣溶膠光學特性、環(huán)境和輻射效應、衛(wèi)星遙感和數(shù)值模式產(chǎn)品的驗證和評估等方面的研究發(fā)揮了重要作用.

    AERONET反演算法[19-20]根據(jù)太陽光度計測量得到的太陽直接輻射和天空散射輻射可以提供大氣柱中的氣溶膠光學特性.利用太陽光度計測得的直接輻射數(shù)據(jù)可以計算各波段的氣溶膠光學厚度(AOD)和 ?ngstr?m 波長指數(shù)(α)[21].利用440,675,870,1020nm波長的等高度角天空掃描輻射數(shù)據(jù)結(jié)合相應波段的 AOD數(shù)據(jù),可以得到氣溶膠粒子尺度分布和復折射指數(shù).在此基礎(chǔ)上,AERONET還提供了其他一些參數(shù),如體積濃度、細粒子體積比、有效半徑、單次散射反照率、相函數(shù)等.

    選取的站點是中國地區(qū)區(qū)域代表性較強且觀測時間序列較長的SACOL、香河和太湖3個AERONET站點,站點信息見表1.SACOL位于甘肅中部的榆中縣,距蘭州市48km,地處西北黃土高原半干旱區(qū),下墊面屬典型的黃土高原地貌,源面梁峁基本為原生植被;香河站位于華北平原,距北京市70km,屬城市遠郊站點,站點周圍為農(nóng)業(yè)生產(chǎn)區(qū),地表植被覆蓋較好;太湖站位于太湖北側(cè)湖岸,屬長江三角洲地區(qū),距無錫市12km,距蘇州市40km,被上海、杭州、南京、蘇州和無錫等大城市包圍,是中國經(jīng)濟發(fā)展最快的地區(qū)之一.

    表1 AERONET站點信息及AOD數(shù)據(jù)樣本數(shù)Table1 Introduction of AERONET sites and the number of AOD data samples

    2 結(jié)果分析

    2.1 氣溶膠光學厚度AOD和?ngstr?m波長指數(shù)α

    圖1為SACOL、香河和太湖3個AERONET站點 AOD(550nm)、?ngstr?m 波長指數(shù)α440-870和大氣柱可降水量的年變化.AERONET公布的數(shù)據(jù)沒有550nm的AOD產(chǎn)品,利用其觀測的其他波段 AOD數(shù)據(jù),采用二次多項式擬合算法[22]求得550nm的 AOD.α440-870是利用波長在440到870nm范圍內(nèi)的可用的AOD數(shù)據(jù)與對應波長對?ngstr?m波長指數(shù)公式進行線性擬合得到的.

    由圖1a可見,SACOL站AOD550具有明顯的年變化,年均值為0.38±0.27,4月 AOD550達到最大值,為0.52±0.35,9月達到最小值,為0.25±0.12;從季節(jié)變化上看,春季最大(0.46),冬季次之(0.45),秋季最小(0.27).α440-870也具有明顯的年變化,年平均值為0.88±0.34,4月達到最小值,為0.54±0.33,7月達到最大值為1.14±0.29;從季節(jié)變化上看,春季(0.60)明顯小于其他季節(jié),冬季較大(0.87),夏秋季分別為1.06和1.04.SACOL春季受沙塵影響非常嚴重,導致 AOD在春季達到最大,沙塵主要為大粒子氣溶膠,使得α在春季達到最小值.SACOL冬季AOD較大是由中國北方冬季燃煤取暖生成的煤煙型氣溶膠導致的.大氣柱可降水量同樣具有顯著的年變化,6~9月平均值超過1.2cm,其他月份小于1.0cm.

    香河站 AOD550的年平均值為0.67±0.66,最大值出現(xiàn)在7月,為1.13±0.94,最小值出現(xiàn)在12月,為0.41±0.40.從季節(jié)變化上看,夏季最大(1.03),春季次之(0.68),冬季最小(0.47).α440-870年平均值為1.09±0.32,只有4月和5月略小于1,其他月份均大于1.可見香河站春季,雖然受到來自中國西部和蒙古戈壁沙漠沙塵的影響,但仍以工業(yè)燃燒和城市污染產(chǎn)生的細粒子為主.香河 AOD夏季最大,是由于夏季風的作用阻礙了北方沙塵對該地區(qū)的影響,導致人為氣溶膠的排放成為其夏季AOD值過高的主要原因.由于華北沿海地區(qū)夏季溫濕度較高,有利于氣粒轉(zhuǎn)化,同時利于增大吸濕性氣溶膠的平均粒徑,從而增加粒子消光效率,最終導致AOD值增大.大氣柱可降水量6~9月平均值超過1.8cm,其他月份小于1.4cm.

    太湖站 AOD550年變化不明顯,平均值為0.72±0.44,最大值出現(xiàn)在6月,為1.12±0.64,最小值出現(xiàn)在12月,為0.56±0.29.α440-870在4月達到最小值,為0.99,9月份最大,為1.34.由此可見,太湖站在春季也會受到沙塵的影響,但沙塵并不是其主要的氣溶膠來源.夏季出現(xiàn)的主要是由于氣粒轉(zhuǎn)化過程生成的細模態(tài)氣溶膠粒子.太湖站,大氣柱可降水量也具有明顯的年變化,6~9月平均值超過3cm,其他月份小于2.0cm,高溫高濕時氣溶膠的吸濕性增長導致的消光效率的增加也可能是AOD增大的原因之一.

    圖2為3個AERONET站點AOD、?ngstr?m波長指數(shù)α和柱可降水量的頻率分布,豎實線表示平均值,豎虛線表示中值.位于中國東部的香河和太湖2個站點,AOD大值出現(xiàn)的頻率明顯高于位于中國西北的SACOL站,而SACOL站α小值出現(xiàn)的頻率明顯高于其他兩個站點,香河站α小值出現(xiàn)的頻率高于太湖站.這是因為中國東部經(jīng)濟發(fā)達,工業(yè)、交通排放的人為氣溶膠較多,從而導致了較大的AOD和較大的α值(細粒子).而位于中國西北黃土高原半干旱區(qū)的 SACOL站,受人類活動的影響相對較小,春季主要受沙塵氣溶膠的影響,冬季由于中國北方燃煤取暖,受到煤煙型氣溶膠的影響,從而導致 SACOL站春冬季AOD出現(xiàn)大值的頻率高于夏秋季,但不同的是春季,SACOL站α出現(xiàn)小值的頻率明顯增大,而冬季α值較大,這一特征在季節(jié)頻率圖上更明顯(未顯示).如SACOL春、夏、秋、冬α小于0.75的頻率分別為70%、20%、14%、35%.

    圖1 AERONET站點觀測的AOD(550nm)、波長指數(shù)α440-870和柱可降水量的年變化Fig.1 Monthly means of AOD (550nm), Angstrom exponent (α440-870), and precipitable water computed from multiyear monitoring at AERONET sites located in SACOL, Xianghe and Taihu

    2.2 氣溶膠粒子尺度譜分布

    AERONET反演算法可以得到0.05~15.0μm范圍內(nèi)22個對數(shù)等間距的半徑處的粒子體積濃度,從而得到氣溶膠粒子體積尺度譜分布(dV/dlnR).

    圖3為3個AERONET站點多年平均的氣溶膠粒子體積尺度譜分布.3個站點粒子體積尺度譜分布均為雙峰型,SACOL站細模態(tài)的峰值較小,以粗模態(tài)粒子為主;香河站粗模態(tài)峰值最大,且粗模態(tài)峰值大于其細模態(tài)峰值;太湖站細模態(tài)峰值最大,且其細模態(tài)的峰值略大于粗模態(tài)的峰值.

    圖4a是SACOL站各月多年平均的氣溶膠粒子尺度分布.SACOL站粗模態(tài)粒子的濃度年變化遠大于細模態(tài)粒子,如粗模態(tài)的峰值從0.028μm3/ μm2(9月)變化到0.211μm3/μm2(4月),變化幅度為0.183μm3/μm2;細模態(tài)的峰值從0.020μm3/μm2(9月)變化到0.029μm3/μm2(2月),變化幅度為0.009μm3/ μm2.春季(3、4和5月)粗模態(tài)粒子濃度顯著大于細模態(tài)粒子濃度,冬季(12、1和2月)粗模態(tài)的粒子濃度小于春季粗模態(tài)的粒子濃度,但明顯大于其對應月份的細模態(tài)粒子濃度.粗模態(tài)和細模態(tài)的粒子體積尺度分布最小值均出現(xiàn)在9月,這與SACOL站AOD的最小值出現(xiàn)月份一致.粗模態(tài)的最大值出現(xiàn)在4月為0.211μm3/μm2,而4月細模態(tài)的峰值僅為0.027μm3/μm2,此時粗模態(tài)對AOD的貢獻占主導地位,而4月的平均α最小,為0.54±0.33,也說明了這一點.

    香河站(圖4b),粗模態(tài)粒子濃度峰值在春季均較大,最大值出現(xiàn)在4月,為0.169μm3/μm2,而7月的粗模態(tài)濃度峰值最小,為0.039μm3/μm2;細模態(tài)粒子濃度峰值在夏季均較大,最大值出現(xiàn)在7月,為0.099μm3/μm2,而1月的細模態(tài)濃度峰值最小,為0.040μm3/μm2.值得注意的是,4月份的粗模態(tài)粒子濃度大于7月份的細模態(tài)粒子濃度,但AOD的最大值出現(xiàn)在7月份.4月份高的粗模態(tài)粒子濃度與α的最小值出現(xiàn)在4月份對應.同時還注意到,各月份粗模態(tài)粒子濃度峰值對應的粒子半徑(模態(tài)半徑)基本不變;而細模態(tài)粒子濃度峰值對應的半徑則從0.15μm(1月)變化到0.33μm(7月),細模態(tài)粒子濃度峰值隨著細模態(tài)半徑的增加而增加.

    圖2 三個AERONET站點AOD、?ngstr?m波長指數(shù)和柱可降水量的頻率分布Fig.2 Frequency of AOD at550nm wavelength, ?ngstr?m exponent, and precipitable water

    圖3 三個AERONET站點多年平均的氣溶膠粒子體積尺度譜分布Fig.3 Multi-annual average of volume particle size distribution at SACOL, Xianghe, and Tiahu sites

    太湖站(圖4c),粗細模態(tài)粒子濃度的峰值相當,粗模態(tài)粒子濃度峰值同樣在春季較大,最大值出現(xiàn)在4月,為0.121μm3/μm2,最小值出現(xiàn)在7月,為0.026μm3/μm2;細模態(tài)粒子濃度峰值在6月最大,為0.126μm3/μm2,在7月最小,為0.057μm3/μm2.太湖站AOD最大值出現(xiàn)月份(6月)與細模態(tài)粒子濃度峰值的最大值出現(xiàn)月份一致.

    2.3 氣溶膠粒子體積濃度

    氣溶膠粒子體積濃度根據(jù)(1)式計算.表2中列出了AOD550、α440-870、總體積濃度VolCon-T、細模態(tài)(粒子半徑r為0.05~0.6μm)體積濃度VolCon-F、粗模態(tài)(粒子半徑 r為0.6~15μm)體積濃度VolCon-C和細模態(tài)體積比(Vf/Vt)等參數(shù)統(tǒng)計值.其中細模態(tài)體積比定義為細模態(tài)體積濃度與總體積濃度之比,即Vf/Vt=VolCon-F/VolCon-T.

    SACOL站,總、細模態(tài)和粗模態(tài)粒子體積濃度分別為0.18、0.041和0.13 μm3/μm2,且其標準差均達到自身的量級,說明數(shù)據(jù)變化范圍較大.總粒子體積濃度的平均值和中值分別為0.18和0.13μm3/μm2,變化達到0.05μm3/μm2(28%).中值小于平均值,說明總粒子體積濃度小于平均值的情況出現(xiàn)的次數(shù)比較多.總粒子體積濃度大于平均值的情況(粗模態(tài)粒子)出現(xiàn)次數(shù)雖然較少,但其對總體積濃度的貢獻達到72%,粗模態(tài)粒子(沙塵)對SACOL站的影響非常重要,甚至可以決定該站氣溶膠總體特征.

    圖4 三個AERONET站點各月多年平均的氣溶膠粒子尺度分布Fig.4 Multi-annual monthly average of volume particle size distribution at SACOL, Xianghe, and Tiahu sites

    表2 三個AERONET站點AOD550、α440-870、體積濃度和細模態(tài)體積比等參數(shù)統(tǒng)計值Table2 Statistical parameters of AOD at550nm, α440-870, volume particle concentration (VolCon) of total, fine and coarse mode, and fine mode volume fraction (Vf/Vt) at SACOL, Xianghe and Taihu sites

    圖5a是SACOL站氣溶膠體積濃度和細模態(tài)粒子體積比的年變化.SACOL總體積濃度和粗模態(tài)體積濃度具有明顯的年變化,且變化趨勢一致,均是9月最小(0.09和0.06μm3/μm2),自9月起增加到4月到達最大(分別為0.37和0.33μm3/ μm2),之后不斷減小.而細模態(tài)體積濃度年變化非常小,最小值出現(xiàn)在9月(0.030μm3/μm2),最大值出現(xiàn)在1月(0.052μm3/μm2).

    香河站,總粒子、細模態(tài)和粗模態(tài)粒子體積濃度分別為0.22、0.082和0.14μm3/μm2.圖5b顯示,總粒子、粗模態(tài)和細模態(tài)粒子體積濃度均具有明顯的年變化.總粒子和粗模態(tài)粒子體積濃度最大值出現(xiàn)在4月(0.33和0.25μm3/μm2),總粒子體積濃度最小值出現(xiàn)在1月(0.16μm3/μm2),粗模態(tài)粒子體積濃度最小值出現(xiàn)在7月(0.066μm3/ μm2).細模態(tài)粒子體積濃度最大值出現(xiàn)在7月(0.16μm3/μm2),最小值出現(xiàn)在12月(0.064μm3/ μm2).

    圖5 三個AERONET站點氣溶膠體積濃度和細模態(tài)體積比的年變化Fig.5 Multi-annual monthly average of the volume particle concentration total, coarse , fine mode and fine mode volume fraction

    太湖站,總粒子、細模態(tài)和粗模態(tài)粒子體積濃度分別為0.23、0.11和0.12μm3/μm2,粗、細模態(tài)的體積濃度相當.從圖5c可以看出,總粒子、粗模態(tài)和細模態(tài)粒子體積濃度也均具有較明顯的年變化.總體積濃度在4月和6月分別出現(xiàn)極大值 (0.30μm3/μm2),最 小 值 出 現(xiàn) 在 7 月(0.16μm3/μm2).粗模態(tài)體積濃度最大值出現(xiàn)在4月(0.20μm3/μm2),最小值出現(xiàn)在8月(0.052μm3/ μm2).細模態(tài)體積濃度最大值出現(xiàn)在6月(0.19μm3/μm2),最小值出現(xiàn)在3月(0.090μm3/ μm2).

    圖6 三個AERONET站點不同細模態(tài)體積比的AOD550與總體積濃度的關(guān)系Fig.6 Scatter plot between AOD550and total volume particle concentration as a function fine mode volume fraction

    圖6是3個AERONET站點在不同細模態(tài)粒子體積比的 AOD550與總體積濃度的散點圖.根據(jù)細模態(tài)體積比的分檔,從圖6可以區(qū)分出右下部(深色)代表的沙塵為主粗粒子氣溶膠和左上部(淺色)代表的細粒子氣溶膠,如大陸型、污染型和生物質(zhì)燃燒型氣溶膠等.每一檔 Vf/Vt對應的AOD550與總體積濃度的擬合直線的斜率即為這一檔對應氣溶膠的消光效率,即單位體積濃度氣溶膠對應的消光(AOD).由此得到SACOL、香河和太湖站的粗粒子平均消光效率分別為:1.44、1.43和1.57μm2/μm3;細粒子消光效率為分別為:4.64、4.97和4.29μm2/μm3.可見細粒子的消光效率均高于粗粒子的消光效率,即同樣體積濃度的粗細粒子,細粒子氣溶膠的AOD要大3~4倍.2.4 細模態(tài)體積比

    細模態(tài)體積比定義為細模態(tài)體積濃度與總體積濃度之比,即 Vf/Vt=VolCon-F/VolCon-T;與?ngstr?m波長指數(shù)α相似,Vf/Vt也可以表征氣溶膠以粗、細何種模態(tài)為主導.α和Vf/Vt的小值代表以粗模態(tài)氣溶膠為主,大值代表以細模態(tài)氣溶膠為主.

    由圖5可見,3個站點Vf/Vt在春季均出現(xiàn)最小值,這主要是由沙塵的影響造成的.同時可以看出Vf/Vt的年變化趨勢與圖1中α的年變化趨勢相似,Vf/Vt和α均可以用來說明粒子尺度大小的年變化特征.圖7為細模態(tài)比的相對頻率和累積頻率的分布.但將圖7(a、b、c)分別與圖2(b、e、h)比較,顯示Vf/Vt的頻率分布特征與α的頻率分布特征并不一致.

    圖8為3個AERONET站點?ngstr?m波長指數(shù)α和細模態(tài)比的散點圖.由圖可見,當α<0.75時,3個站點Vf/Vt均小于0.4,以粗模態(tài)粒子為主,這種情況在SACOL、香河和太湖出現(xiàn)的概率分別為33%、15%和8%;此時α的大小可以較好的說明粒子的大小,如當α<0.5時,3個站點Vf/Vt均小于0.2,即粗模態(tài)粒子占了80%以上,這種情況在SACOL、香河和太湖出現(xiàn)的概率分別為16%、6%和2%.當α>1.75時,SACOL和太湖Vf/Vt大于0.6,以細模態(tài)粒子為主,但這一情況出現(xiàn)的概率小于0.5%.而當0.75<α<1.7時,也是α出現(xiàn)概率最高的區(qū)間,3個站點Vf/Vt均在0.1~0.8之間變化,此時僅利用α無法判斷何種氣溶膠粒子占主導地位.

    圖7 細模態(tài)比的相對頻率和累積頻率的分布Fig.7 Relative and cumulative frequency of fine mode volume fraction

    2.5 有效半徑

    Hansen等[23]用(2)式定義有效半徑,即以粒子的幾何截面積(πr2)為權(quán)重的平均半徑,并指出有效半徑是球形粒子的任何物理上合理的尺度分布的首要特征參數(shù).有效半徑經(jīng)常被用于涉及尺度分布的計算中,主要是因為球形粒子的散射輻射與其幾何截面積成正比,有效半徑與尺度分布的單次散射特性關(guān)系比中值半徑更密切[24-26]也就是說,雖然對數(shù)正態(tài)分布的形式并不統(tǒng)一,但只要氣溶膠尺度分布有相同的有效半徑和半徑譜寬,則具有相近的單次散射反照率.

    表3列出了SACOL、香河和太湖氣溶膠站粒子有效半徑的各統(tǒng)計值.總粒子的有效半徑的平均值在 SACOL最大(0.52μm),香河站居中(0.45μm),太湖站最小(0.34μm).粗模態(tài)粒子有效半徑平均值香河站最大,SACOL居中,太湖站最小.細模態(tài)粒子有效半徑平均值各站點基本相等.各站總粒子和粗模態(tài)粒子的有效半徑的變化幅度較大,而細模態(tài)粒子的有效半徑變化很小.總粒子、細模態(tài)和粗模態(tài)粒子有效半徑的平均值和中值均比較接近,說明氣溶膠粒子的尺度分布比較接近對數(shù)正態(tài)分布.

    圖9是總粒子、細模態(tài)和粗模態(tài)粒子有效半徑的多年月平均值的年變化.各站點總粒子有效半徑的平均值均在春冬季較大,夏秋季較小.SACOL站最大值出現(xiàn)在3、4月,為0.67μm,最小值出現(xiàn)在10月,為0.39μm;香河站最大值出現(xiàn)在4月,為0.67μm,之后減小,至8月出現(xiàn)最小值(0.32μm)后又逐漸增大;太湖站最大值出現(xiàn)在3、4月,為0.41μm,7月減小到最小值,為0.28μm,后逐漸增大.

    表3 三個AERONET氣溶膠粒子有效半徑的統(tǒng)計值Table3 Statistical values of effective radius for total, fine and coarse mode at SACOL, Xianghe, and Tiahu sites

    圖9 三個AERONET站點氣溶膠粒子有效半徑的年變化Fig.9 Monthly average of the effective radius of total, fine and coarse mode at SACOL, Xianghe, and Taihu sites

    細模態(tài)有效半徑年變化幅度較小.SACOL細模態(tài)有效半徑春季最小,為0.13μm,冬季最大,為0.17μm.香河春季最小,為0.14μm,而最大值出現(xiàn)在夏季,為0.19μm.太湖與香河類似,春季最小,為0.15μm,夏季最大,為0.19μm.可以看到3個站點細模態(tài)有效半徑均在春季最小,而3個站點春季都受到沙塵的影響,可見由于沙塵氣溶膠的出現(xiàn),會使細模態(tài)有效半徑有所減小.粗模態(tài)有效半徑年變化不明顯,但可以看到3個站點粗模態(tài)有效半徑同樣均在春季最小.可見,春季沙塵的入侵,不僅會使細模態(tài)有效半徑減小,同樣也會使粗模態(tài)的有效半徑減小.

    3 結(jié)論

    3.1 香河和太湖多年平均 AOD分別為0.67±0.66和0.72±0.44,是SACOL AOD平均值(0.38±0.27)的近2倍,且AOD變化范圍較大.SACOL 春冬季AOD較大,夏秋季AOD較小;而香河和太湖是夏春季較大,秋冬季較小.

    3.2 SACOL α明顯小于香河和太湖,且三站春季α值均小于其年平均值,特別是SACOL站要小的多,說明春季沙塵氣溶膠對三站均有影響,且在SACOL影響最為嚴重,結(jié)合尺度分布特征可以說明沙塵是SACOL春季最主要的氣溶膠類型.

    3.3 SACOL總體積濃度和粗模態(tài)體積濃度具有明顯的年變化,且趨勢一致,而細模態(tài)體積濃度全年均比較小.香河,總粒子、粗和細模態(tài)粒子體積濃度均具有較明顯的年變化,總粒子和粗模態(tài)粒子體積濃度最大值(0.33和0.25 μm3/μm2)出現(xiàn)在4月,但AOD的最大值卻出現(xiàn)在細模態(tài)體積濃度最大的7月(0.16μm3/μm2).太湖,總粒子和粗模態(tài)粒子體積濃度最大值(0.30,0.20μm3/μm2)也出現(xiàn)在4月,但AOD的最大值出現(xiàn)在細模態(tài)體積濃度最大的6月(0.19μm3/μm2).這是由于細模態(tài)粒子的消光效率是粗模態(tài)粒子的3~4倍.

    3.4 Vf/Vt的年變化趨勢與α的年變化趨勢相似,Vf/Vt和α均可以用來分析粒子尺度大小的年變化特征. 但α<0.75時,3個站點Vf/Vt均小于0.4,以粗模態(tài)粒子為主; α>1.7時,Vf/Vt大于0.6,以細模態(tài)粒子為主; 而0.75<α<1.7時,也是α出現(xiàn)概率最高的區(qū)間,3個站點Vf/Vt均在0.1~0.8之間變化,此時僅利用α無法判斷是以何種模態(tài)粒子為主.

    3.5 各站點總粒子有效半徑的平均值均在春冬季較大,夏秋季較小. 粗、細模態(tài)有效半徑的年變化不大,春季沙塵的入侵,會使粗、細模態(tài)有效半徑有所減小.

    [1] Twomey S. Atmospheric aerosols [M]. New York, USA: Elsevier Scientific Pub. Co.,1977.

    [2] Charlson R J, Lovelock J E, Andreae M O, et al. Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate [J]. Nature,1987,326(6114):655-661.

    [3] Lohmann U, Feichter J. Global indirect aerosol effects: a review [J]. Atmos. Chem. Phys.,2005,5(3):715-737.

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

    [5] Koren I, Kaufman Y J, Remer L A, et al. Measurement of the effect of Amazon smoke on inhibition of cloud formation [J]. Science,2004,303(5662):1342-1345.

    [6] B?umer D, Vogel B, Versick S, et al. Relationship of visibility, aerosol optical thickness and aerosol size distribution in an ageing air mass over South-West Germany [J]. Atmospheric Environment,2008,42(5):989-998.

    [7] Wang K, Dickinson R E, Liang S. Clear sky visibility has decreased over land globally from1973 to2007 [J]. Science,2009,323(5920):1468-1470.

    [8] Gauderman W J, Mcconnell R, Gilliland F, et al. Association between air pollution and lung function growth in southern California children [J]. American Journal of Respiratory and Critical care Medicine,2000,162(4):1383-1390.

    [9] 毛節(jié)泰,張軍華,王美華.中國大氣氣溶膠研究綜述 [J]. 氣象學報,2002,60(5):625-634.

    [10] 石廣玉,王 標,張 華,等.大氣氣溶膠的輻射與氣候效應 [J].大氣科學,2008,32(4):826-840.

    [11] 張小曳.中國大氣氣溶膠及其氣候效應的研究 [J]. 地球科學進展,2007,22(1):12-16.

    [12] 沈凡卉,王體健,莊炳亮,等.中國沙塵氣溶膠的間接輻射強迫與氣候效應 [J]. 中國環(huán)境科學,2011,31(7):1057-1063.

    [13] Climate change2007-the physical science basis: Working group I contribution to the fourth assessment report of the IPCC [M]. Cambridge University Press,2007.

    [14] Hansen J, Sato M, Kharecha P, et al. Earth’s energy imbalance and implications [J]. Atmos Chem Phys,2011,11(24):13421-13449.

    [15] Holben B N, Eck T F, Slutsker I, et al. AERONET-A federated instrument network and data archive for aerosol characterization [J]. Remote Sensing of Environment,1998,66(1):1-16.

    [16] AOKI K, Fujiyoshi Y. Sky radiometer measurements of aerosoloptical properties over Sapporo, Japan [J]. Journal of the Meteorological Society of Japan. Ser. II,2003,81(3):493-513.

    [17] Redemann J, Russell P B, Hamill P. Dependence of aerosol light absorption and single scattering albedo on ambient relative humidity for sulfate aerosols with black carbon cores [J]. J. Geophys. Res.,2001,106(D21):27485-27495.

    [18] 錢 凌,銀 燕,童堯青,等.南京北郊大氣細顆粒物的粒徑分布特征 [J]. 中國環(huán)境科學,2008,28(1):18-22.

    [19] Dubovik O, King M D. A flexible inversion algorithm for retrieval of aerosol optical properties from Sun and sky radiance measurements [J]. J. Geophys. Res.,2000,105(D16):20673-20696.

    [20] Dubovik O, Sinyuk A, Lapyonok T, et al. Application of spheroid models to account for aerosol particle nonsphericity in remote sensing of desert dust [J]. J. Geophys. Res.,2006,111(D11).

    [21] 王 靜,牛生杰,許 丹,等.南京一次典型霧霾天氣氣溶膠光學特性 [J]. 中國環(huán)境科學,2013,33(2):201-208.

    [22] Eck T F, Holben B N, Reid J S, et al. Wavelength dependence of the optical depth of biomass burning, urban, and desert dust aerosols [J]. J. Geophys. Res.,1999,104(D24):31333-31349.

    [23] Hansen J E, Travis L D. Light scattering in planetary atmospheres [J]. Space Science Reviews,1974,16(4):527-610.

    [24] Liou K N, Davies R. Radiation and Cloud Processes in the Atmosphere [J]. Physics Today,1993,46:66.

    [25] Mishchenko M I, Travis L D, Kahn R A, et al. Modeling phase functions for dustlike tropospheric aerosols using a shape mixture of randomly oriented polydisperse spheroids [J]. J. Geophys. Res.,1997,102(D14):16831-16847.

    [26] 孫賢明,哈恒旭.基于反射太陽光反演氣溶膠光學厚度和有效半徑 [J]. 物理學報,2008,57(9):5565-5570.

    致謝:本文資料由全球氣溶膠監(jiān)測網(wǎng)(AERONET)和香河站、太湖站、蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)提供,謹致謝忱!

    德國將中止頁巖氣鉆探, 因勘探技術(shù)或污染地下水

    據(jù)美國《華爾街日報》報道,德國計劃未來七年中止頁巖氣鉆探,原因是擔心頁巖氣的勘探技術(shù)可能會污染地下水.

    據(jù)報道,頁巖氣的開采需要使用水力壓裂技術(shù).這種技術(shù)是用摻入水、沙子和化學物質(zhì)的高壓混合物將巖石層壓裂,從而釋放出天然氣.德國政府計劃禁止在離地表不到3000m的鉆探作業(yè)中使用水力壓裂技術(shù),并希望在明年初制定一項法令.這一禁令將于2021年被重新評估.

    德國環(huán)境部長亨德里克斯表示,德國在可預見的未來將不會有頁巖氣壓裂.

    報道指出,在德國提出這一擬定的規(guī)定之際,德國與其主要天然氣供應國俄羅斯正陷入政治僵局中.另外,環(huán)保人士等此前也因擔心頁巖氣鉆探可能污染地下水而進行了密集的游說.

    亨德里克斯稱,保護飲用水和健康是德國最為重視的事情.

    摘自《中國環(huán)境報》

    2014-07-07

    Analysis and comparison of the aerosol microphysical properties at three AERONET sites in China.


    ZHANG Zhi-wei1,2, WANG Hong-bin3, ZHANG Lei1*, HU Die4, TIAN Peng-fei1, LIANG Jie-ning1, WANG Shi-gong1
    (1.Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, College of Atmospheric Sciences, Lanzhou University, Lanzhou730000, China;2.Jiangsu Meteorological Service Center, Nanjing210008, China;3.Jiangsu Institute of Meteorological Sciences, Nanjing210009, China;4.Lanzhou Institute of Arid Meteorology, China Meteorological Administration, Lanzhou730020, China). China Environmental Science,2014,34(8):1927~1937

    The microphysical parameter characteristics of aerosol were analyzed at the AERONET stations located in SACOL, Xianghe and Taihu. The results showed that the annual average AOD at Xianghe and Taihu stations were0.67±0.66 and0.72±0.44, respectively, which were nearly2times of the value of AOD at SACOL (0.38±0.27). The standard deviation presented a significant temporal variability. At SACOL, the higher value of AOD was found in the spring and winter, and the lower value of AOD was found in the summer and fall. However, At Xianghe and Taihu, the higher values of AOD were found in the spring and summer, and the lower values of AOD were found in the fall and winter. Combined with the size distribution, volume concentration of aerosol particles, it can be concluded that dust was the main type of aerosols at SACOL in spring. At Xianghe and Taihu stations, the types of aerosol were mainly fine mode particles in summer due to the gas-to-particle conversion process. We could observe that for the measurements with ?ngstr?m exponent (α<0.75), there was a prevalence of coarse particles, because for the lower α values, data with the fine mode volume fraction (Vf/Vt>0.4) was almost unavailable. According to analysis we could establish the limit of the fine particle predominance with α>1.7approximately. But for0.75<α<1.7, the range of Vf/Vt was predominantly between0.2and0.8, which represented a great scattering character, hence it was not possible to establish a priori predominance for the fine or coarse particles.

    t:aerosol optical depth(AOD);volume concentration;fine mode volume fraction;?ngstr?m exponent

    X513

    :A

    :1000-6923(2014)08-1927-11

    張志薇(1987-),女,廣東大埔人,博士,主要從事大氣環(huán)境和人體健康研究.

    2013-11-04

    國家重大科學研究計劃資助項目(2012CB955302);國家自然科學基金資助項目(41075104);中國氣象局氣溶膠與云降水重點開放實驗室開放課題(KDW1303);蘭州大學中央高?;究蒲袠I(yè)務費專項資金(lzujbky-2013-246)

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

    猜你喜歡
    香河氣溶膠太湖
    堅守中的回望、凝思與嬗變
    氣溶膠傳播之謎
    雁飛潮白六十年
    博覽群書(2019年3期)2019-04-16 18:34:42
    氣溶膠中210Po測定的不確定度評定
    香河旅游推介走進鳥巢
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    四川盆地秋季氣溶膠與云的相關(guān)分析
    養(yǎng)螃蟹
    金山(2016年9期)2016-10-12 14:09:07
    91久久精品国产一区二区成人| 亚洲精品国产成人久久av| 一区二区三区免费毛片| 国产精品久久久久久精品古装| 久久久精品免费免费高清| 永久网站在线| 国产成人免费观看mmmm| 午夜福利影视在线免费观看| 一级爰片在线观看| 亚洲第一区二区三区不卡| 精品久久久久久久久亚洲| 我要看日韩黄色一级片| 狂野欧美白嫩少妇大欣赏| 中文在线观看免费www的网站| 国产精品一区二区三区四区免费观看| 少妇丰满av| 激情 狠狠 欧美| 国产美女午夜福利| 成年女人在线观看亚洲视频| 国产成人精品久久久久久| 日日撸夜夜添| 大香蕉97超碰在线| av在线老鸭窝| 亚洲成人av在线免费| 黄片wwwwww| 日韩成人av中文字幕在线观看| 午夜激情福利司机影院| 另类亚洲欧美激情| 成人毛片60女人毛片免费| 黑人高潮一二区| 日韩欧美精品免费久久| 国产日韩欧美在线精品| 中文天堂在线官网| 久久久久人妻精品一区果冻| av女优亚洲男人天堂| 久久久久久久亚洲中文字幕| 成年美女黄网站色视频大全免费 | 80岁老熟妇乱子伦牲交| 自拍偷自拍亚洲精品老妇| 国产视频首页在线观看| 青春草国产在线视频| 午夜视频国产福利| 一区二区三区精品91| 国产精品一区二区在线不卡| 国产一级毛片在线| 国产白丝娇喘喷水9色精品| 精品国产乱码久久久久久小说| 免费观看性生交大片5| 热99国产精品久久久久久7| 国产精品伦人一区二区| 久久久久国产网址| 久久99热这里只有精品18| 亚洲成人手机| 国产乱人偷精品视频| 国产黄片美女视频| av又黄又爽大尺度在线免费看| 高清欧美精品videossex| 国产 一区 欧美 日韩| 美女脱内裤让男人舔精品视频| 国产成人91sexporn| 好男人视频免费观看在线| 免费黄网站久久成人精品| 国产在线免费精品| 七月丁香在线播放| 国产片特级美女逼逼视频| 国产色婷婷99| 自拍偷自拍亚洲精品老妇| 一级毛片久久久久久久久女| 蜜桃久久精品国产亚洲av| 黄色配什么色好看| 亚洲欧美成人综合另类久久久| 一二三四中文在线观看免费高清| 观看免费一级毛片| 两个人的视频大全免费| 国产高清国产精品国产三级 | 国内揄拍国产精品人妻在线| 春色校园在线视频观看| 国产在线男女| 精华霜和精华液先用哪个| 久久精品久久精品一区二区三区| 老师上课跳d突然被开到最大视频| 九九在线视频观看精品| 国产精品欧美亚洲77777| 欧美精品亚洲一区二区| a级一级毛片免费在线观看| 91午夜精品亚洲一区二区三区| 久久久精品94久久精品| 免费看日本二区| 亚洲中文av在线| 亚洲美女视频黄频| 国模一区二区三区四区视频| 高清不卡的av网站| 极品教师在线视频| 欧美日韩一区二区视频在线观看视频在线| 黑人高潮一二区| 国产人妻一区二区三区在| 日韩欧美一区视频在线观看 | 亚洲精品国产色婷婷电影| 亚洲人与动物交配视频| 成人高潮视频无遮挡免费网站| 少妇被粗大猛烈的视频| 99久久精品国产国产毛片| av黄色大香蕉| 99re6热这里在线精品视频| 两个人的视频大全免费| 国产成人精品久久久久久| 久久久亚洲精品成人影院| 亚洲熟女精品中文字幕| 国产亚洲91精品色在线| 一二三四中文在线观看免费高清| 下体分泌物呈黄色| 亚洲美女搞黄在线观看| 中国美白少妇内射xxxbb| 三级国产精品欧美在线观看| 最黄视频免费看| 国国产精品蜜臀av免费| 久久 成人 亚洲| 午夜激情福利司机影院| 国产深夜福利视频在线观看| 亚洲在久久综合| 天天躁夜夜躁狠狠久久av| 亚洲精品乱久久久久久| 91精品国产国语对白视频| 少妇熟女欧美另类| 国产成人精品福利久久| 少妇的逼好多水| 少妇 在线观看| 国产精品99久久久久久久久| 美女国产视频在线观看| 国产成人aa在线观看| 欧美日本视频| 欧美bdsm另类| 在线免费十八禁| 日日摸夜夜添夜夜添av毛片| 久久久久久伊人网av| 免费观看的影片在线观看| 国产精品麻豆人妻色哟哟久久| 在线看a的网站| 精品久久国产蜜桃| 特大巨黑吊av在线直播| 欧美三级亚洲精品| 激情 狠狠 欧美| 中文字幕免费在线视频6| 久久国产精品大桥未久av | 啦啦啦中文免费视频观看日本| 国产久久久一区二区三区| 肉色欧美久久久久久久蜜桃| 国产视频首页在线观看| 又爽又黄a免费视频| 午夜福利网站1000一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 黄色怎么调成土黄色| 国产91av在线免费观看| 纵有疾风起免费观看全集完整版| 久久综合国产亚洲精品| 精品亚洲成a人片在线观看 | 老熟女久久久| 大话2 男鬼变身卡| 最近中文字幕高清免费大全6| 直男gayav资源| 国产有黄有色有爽视频| 国产精品一区二区在线不卡| 久热久热在线精品观看| 高清午夜精品一区二区三区| 成人影院久久| 老女人水多毛片| 国产黄色免费在线视频| av国产久精品久网站免费入址| 国产精品麻豆人妻色哟哟久久| 全区人妻精品视频| 国产精品无大码| 网址你懂的国产日韩在线| 亚洲天堂av无毛| 国产成人免费无遮挡视频| 特大巨黑吊av在线直播| 岛国毛片在线播放| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 国产精品蜜桃在线观看| 女性生殖器流出的白浆| 久久久久精品久久久久真实原创| 一二三四中文在线观看免费高清| 99热6这里只有精品| 久久久a久久爽久久v久久| 黑丝袜美女国产一区| 日韩制服骚丝袜av| 美女福利国产在线 | 日韩在线高清观看一区二区三区| 男女啪啪激烈高潮av片| 又粗又硬又长又爽又黄的视频| 国产精品嫩草影院av在线观看| 免费黄频网站在线观看国产| 人人妻人人澡人人爽人人夜夜| 最新中文字幕久久久久| 国产成人免费无遮挡视频| 欧美日韩亚洲高清精品| 免费大片黄手机在线观看| 国产一区二区三区av在线| 久久人人爽av亚洲精品天堂 | 丰满乱子伦码专区| 精品亚洲成a人片在线观看 | 青春草视频在线免费观看| 亚洲天堂av无毛| 亚洲av免费高清在线观看| 欧美+日韩+精品| 国产精品国产三级国产av玫瑰| 韩国av在线不卡| 九九在线视频观看精品| 国精品久久久久久国模美| 美女中出高潮动态图| 最黄视频免费看| 丝袜脚勾引网站| 久久精品国产亚洲av涩爱| 黑人高潮一二区| 极品教师在线视频| 日韩av在线免费看完整版不卡| 男男h啪啪无遮挡| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 如何舔出高潮| 国产精品偷伦视频观看了| 高清欧美精品videossex| 欧美激情国产日韩精品一区| 久久久精品94久久精品| 欧美97在线视频| 欧美bdsm另类| 蜜桃在线观看..| 只有这里有精品99| 国产人妻一区二区三区在| 三级经典国产精品| 日韩成人av中文字幕在线观看| 亚洲第一av免费看| 久久99热这里只有精品18| 日韩一区二区三区影片| 国产高潮美女av| 天美传媒精品一区二区| 狠狠精品人妻久久久久久综合| 久久这里有精品视频免费| 热99国产精品久久久久久7| 国产精品久久久久久精品古装| 欧美高清成人免费视频www| 亚洲一区二区三区欧美精品| 亚洲精品国产色婷婷电影| 国产精品国产三级国产专区5o| 国产探花极品一区二区| 一区二区av电影网| av国产久精品久网站免费入址| 国产精品久久久久久精品古装| 91精品国产国语对白视频| 91精品国产九色| 女人久久www免费人成看片| 久久99热6这里只有精品| 国产毛片在线视频| 日韩av在线免费看完整版不卡| 激情五月婷婷亚洲| a 毛片基地| 成年女人在线观看亚洲视频| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 香蕉精品网在线| 久久久久久久久久人人人人人人| 精品午夜福利在线看| 最近的中文字幕免费完整| 99久久精品一区二区三区| av在线观看视频网站免费| 一区二区三区四区激情视频| av黄色大香蕉| 边亲边吃奶的免费视频| 成人影院久久| 特大巨黑吊av在线直播| 99视频精品全部免费 在线| 777米奇影视久久| 人妻系列 视频| 日本-黄色视频高清免费观看| 亚洲色图综合在线观看| 内射极品少妇av片p| 亚洲欧洲日产国产| 国产免费又黄又爽又色| 26uuu在线亚洲综合色| 秋霞在线观看毛片| 中文字幕av成人在线电影| 春色校园在线视频观看| 国产av码专区亚洲av| 最近中文字幕高清免费大全6| 国产午夜精品一二区理论片| 色哟哟·www| 99热这里只有精品一区| 国产精品偷伦视频观看了| 下体分泌物呈黄色| 国产高潮美女av| 插阴视频在线观看视频| 免费av不卡在线播放| 久久青草综合色| 蜜桃在线观看..| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 精品国产三级普通话版| 成人二区视频| 国产精品无大码| 亚洲精品第二区| 三级国产精品欧美在线观看| 99久国产av精品国产电影| 黄色日韩在线| 卡戴珊不雅视频在线播放| 国产精品99久久久久久久久| 天堂8中文在线网| 亚洲av中文字字幕乱码综合| 欧美国产精品一级二级三级 | 国产69精品久久久久777片| 亚洲精品aⅴ在线观看| 日韩免费高清中文字幕av| 色综合色国产| 久久6这里有精品| 久久毛片免费看一区二区三区| 啦啦啦在线观看免费高清www| 国产亚洲av片在线观看秒播厂| 99热国产这里只有精品6| 成年免费大片在线观看| 只有这里有精品99| 热re99久久精品国产66热6| 亚洲av中文av极速乱| 简卡轻食公司| 亚洲精品aⅴ在线观看| 高清欧美精品videossex| 亚洲成人av在线免费| 男男h啪啪无遮挡| 久久久久国产网址| 在线观看免费高清a一片| 亚洲欧美日韩东京热| 亚洲国产精品一区三区| 麻豆成人av视频| 国产精品爽爽va在线观看网站| 成年女人在线观看亚洲视频| 国产精品三级大全| 久久久欧美国产精品| 国产男女内射视频| av专区在线播放| 国产成人精品婷婷| 美女中出高潮动态图| av又黄又爽大尺度在线免费看| 在线观看人妻少妇| 久久久久视频综合| 国产成人精品一,二区| 18禁裸乳无遮挡动漫免费视频| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 久久精品国产a三级三级三级| 又大又黄又爽视频免费| 九九爱精品视频在线观看| 欧美日韩精品成人综合77777| 亚洲成人一二三区av| 在线观看免费日韩欧美大片 | 99久国产av精品国产电影| 夜夜看夜夜爽夜夜摸| 性色avwww在线观看| 夫妻性生交免费视频一级片| 最近手机中文字幕大全| 亚洲欧美一区二区三区黑人 | a级一级毛片免费在线观看| 国产亚洲欧美精品永久| 午夜免费观看性视频| 色哟哟·www| 日韩一区二区视频免费看| 精品久久久久久久久av| 人妻 亚洲 视频| 亚洲欧洲日产国产| 国产精品不卡视频一区二区| 久久综合国产亚洲精品| 久久久久国产网址| 久久午夜福利片| 人人妻人人澡人人爽人人夜夜| 美女内射精品一级片tv| 人体艺术视频欧美日本| av不卡在线播放| 久久人妻熟女aⅴ| 日本欧美国产在线视频| 国产成人精品久久久久久| 久久久色成人| 国产成人精品久久久久久| 亚洲精品日本国产第一区| 两个人的视频大全免费| 国产中年淑女户外野战色| 午夜福利在线在线| 日本黄色片子视频| av国产免费在线观看| 欧美 日韩 精品 国产| 黄片wwwwww| 中文字幕亚洲精品专区| 日本av免费视频播放| 99久久精品国产国产毛片| 最近中文字幕2019免费版| 亚洲综合精品二区| 交换朋友夫妻互换小说| 久久久久性生活片| 亚洲欧美精品专区久久| 中文字幕免费在线视频6| 亚洲怡红院男人天堂| videos熟女内射| 国产一区亚洲一区在线观看| 成人午夜精彩视频在线观看| 欧美一区二区亚洲| 国产成人精品福利久久| 天美传媒精品一区二区| 精品亚洲乱码少妇综合久久| 日本-黄色视频高清免费观看| 久久久久久九九精品二区国产| 免费大片18禁| 观看免费一级毛片| 日日啪夜夜撸| 午夜老司机福利剧场| 欧美+日韩+精品| 色哟哟·www| 国产 一区 欧美 日韩| 欧美日韩精品成人综合77777| 狂野欧美白嫩少妇大欣赏| 日本免费在线观看一区| 久久99热这里只频精品6学生| 亚洲av中文字字幕乱码综合| 日韩三级伦理在线观看| 亚洲av在线观看美女高潮| 中文欧美无线码| 亚洲欧美一区二区三区国产| 高清午夜精品一区二区三区| 国产免费视频播放在线视频| 国产成人一区二区在线| 国国产精品蜜臀av免费| av在线老鸭窝| 日本wwww免费看| 美女视频免费永久观看网站| 这个男人来自地球电影免费观看 | 中文字幕制服av| 日韩视频在线欧美| 亚洲精品视频女| 国产免费一区二区三区四区乱码| 久久毛片免费看一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 成年av动漫网址| 国产一区二区三区av在线| 中文乱码字字幕精品一区二区三区| 久久久久久久久久久免费av| 搡老乐熟女国产| 联通29元200g的流量卡| 成人国产av品久久久| 国产精品99久久99久久久不卡 | 国产免费视频播放在线视频| 亚洲成人av在线免费| 国产精品嫩草影院av在线观看| 亚洲精品乱码久久久久久按摩| 国产无遮挡羞羞视频在线观看| 男女下面进入的视频免费午夜| av在线app专区| 一级片'在线观看视频| 国产成人午夜福利电影在线观看| 啦啦啦中文免费视频观看日本| 成人亚洲欧美一区二区av| 老司机影院成人| 亚洲精品乱码久久久久久按摩| 国产精品伦人一区二区| 欧美日韩在线观看h| 久久久久国产精品人妻一区二区| 五月伊人婷婷丁香| 少妇精品久久久久久久| 狂野欧美激情性xxxx在线观看| 不卡视频在线观看欧美| 午夜福利影视在线免费观看| 婷婷色麻豆天堂久久| 亚洲精品中文字幕在线视频 | 久久人人爽人人爽人人片va| 国产成人一区二区在线| 丰满少妇做爰视频| 男的添女的下面高潮视频| 啦啦啦啦在线视频资源| 看非洲黑人一级黄片| 免费黄频网站在线观看国产| 久久ye,这里只有精品| 国产精品一区二区三区四区免费观看| 夜夜骑夜夜射夜夜干| 亚洲精品日韩在线中文字幕| 在线亚洲精品国产二区图片欧美 | 亚洲国产av新网站| 国产欧美日韩精品一区二区| 国产高清有码在线观看视频| 国产av码专区亚洲av| 久久精品夜色国产| 国内揄拍国产精品人妻在线| 男人添女人高潮全过程视频| 校园人妻丝袜中文字幕| 欧美三级亚洲精品| 久久久久性生活片| 久久久成人免费电影| 日韩欧美一区视频在线观看 | 亚洲人成网站在线播| 亚洲精品乱久久久久久| 中文资源天堂在线| 王馨瑶露胸无遮挡在线观看| 男的添女的下面高潮视频| 蜜桃亚洲精品一区二区三区| 各种免费的搞黄视频| 在线播放无遮挡| 大片免费播放器 马上看| av在线蜜桃| 国产精品国产三级国产专区5o| 一区二区av电影网| av视频免费观看在线观看| 啦啦啦视频在线资源免费观看| 免费观看av网站的网址| 建设人人有责人人尽责人人享有的 | 久久久欧美国产精品| 亚洲,一卡二卡三卡| 国产精品一区二区性色av| 久久99热6这里只有精品| 少妇人妻久久综合中文| 人妻 亚洲 视频| 少妇猛男粗大的猛烈进出视频| 制服丝袜香蕉在线| 欧美3d第一页| 日日撸夜夜添| 91狼人影院| 丰满少妇做爰视频| 免费人妻精品一区二区三区视频| 中文在线观看免费www的网站| 国产精品爽爽va在线观看网站| 日本黄色日本黄色录像| 国产淫片久久久久久久久| 精品人妻视频免费看| 国产高清三级在线| 视频中文字幕在线观看| 寂寞人妻少妇视频99o| 亚洲av综合色区一区| 中国三级夫妇交换| 欧美人与善性xxx| 毛片女人毛片| 久久国产乱子免费精品| 欧美97在线视频| 亚洲色图综合在线观看| 欧美精品人与动牲交sv欧美| 久久久久久人妻| 免费av中文字幕在线| 久久女婷五月综合色啪小说| 99热全是精品| 国产精品国产三级国产av玫瑰| 91狼人影院| 国产精品久久久久成人av| 亚洲精品自拍成人| 精品人妻视频免费看| 日本欧美视频一区| 亚洲色图综合在线观看| 欧美精品人与动牲交sv欧美| 国产久久久一区二区三区| 国产亚洲最大av| 人体艺术视频欧美日本| 啦啦啦啦在线视频资源| 中文精品一卡2卡3卡4更新| 久久久久国产网址| 国产黄色视频一区二区在线观看| 精品一区二区免费观看| 日韩成人伦理影院| 大香蕉97超碰在线| 国产一区二区在线观看日韩| a级一级毛片免费在线观看| 99久久综合免费| 国产精品国产三级专区第一集| 97在线人人人人妻| 亚洲国产欧美人成| 午夜免费观看性视频| 亚洲国产精品一区三区| 久久韩国三级中文字幕| 国产成人精品福利久久| 欧美xxxx黑人xx丫x性爽| 国产免费福利视频在线观看| 麻豆乱淫一区二区| 街头女战士在线观看网站| 熟女人妻精品中文字幕| 日产精品乱码卡一卡2卡三| 久久精品国产亚洲网站| 亚洲精品日韩av片在线观看| 欧美高清性xxxxhd video| 日日撸夜夜添| 精品一品国产午夜福利视频| 一级二级三级毛片免费看| 欧美激情极品国产一区二区三区 | 青春草国产在线视频| 成人黄色视频免费在线看| 久久久久久久久久久丰满| 一级片'在线观看视频| 熟女av电影| 久久 成人 亚洲| 亚洲第一区二区三区不卡| 国产中年淑女户外野战色| 五月伊人婷婷丁香| 一个人免费看片子| 国产精品爽爽va在线观看网站| 亚洲最大成人中文| 午夜福利影视在线免费观看| 日韩人妻高清精品专区| 婷婷色av中文字幕| 七月丁香在线播放| 久久毛片免费看一区二区三区| av国产精品久久久久影院| 国产av精品麻豆| 亚洲欧美一区二区三区国产| av在线app专区| 久久久久久久久大av| 黄色欧美视频在线观看| 99九九线精品视频在线观看视频| 久久久久久久久久成人| 亚洲精品久久午夜乱码| .国产精品久久| 国产成人aa在线观看| 欧美日韩精品成人综合77777| 丝袜喷水一区| 国产亚洲最大av| 岛国毛片在线播放| 国产精品久久久久久av不卡| 国产熟女欧美一区二区|