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

    浙江省氣溶膠垂直分布及區(qū)域傳輸分析

    2021-05-28 16:08:46張亦舒李冬會
    中國環(huán)境科學 2021年5期
    關鍵詞:氣溶膠激光雷達高空

    張亦舒,徐 達,何 秦,李冬會,秦 凱*

    浙江省氣溶膠垂直分布及區(qū)域傳輸分析

    張亦舒1,徐 達2,何 秦1,李冬會1,秦 凱1*

    (1.中國礦業(yè)大學環(huán)境與測繪學院,江蘇 徐州 221116;2.浙江省生態(tài)環(huán)境監(jiān)測中心,浙江 杭州 310012)

    利用6個激光雷達站點數(shù)據(jù)從年際和季節(jié)變化尺度對浙江省2017~2018年氣溶膠垂直分布特征進行研究,并結合污染物濃度數(shù)據(jù)分析了浙江省主要顆粒物類型的季節(jié)差異. 結果表明,大部分站點2018年氣溶膠消光系數(shù)整體小于2017年; 邊界層內消光系數(shù)空間上呈中部高東西低的分布形態(tài); 消光系數(shù)季節(jié)變化表現(xiàn)為冬季最高,秋季其次,而春夏較低;氣溶膠在春、夏季集中在地面上方1~2km以內,冬季主要分布在1km以下,秋季受高空氣溶膠影響,6個站點呈現(xiàn)出不同的垂直分布模式. 進一步, 選取發(fā)生在2018年4月和7月兩次氣溶膠跨區(qū)域傳輸過程,綜合CALIPSO和MODIS衛(wèi)星數(shù)據(jù)、以及后向軌跡?PSCF與CWT模式對污染物來源與傳輸機制進行探索, 發(fā)現(xiàn)來自AOD高值區(qū)的氣團顯著加重本地污染.

    微脈沖激光雷達;氣溶膠;消光系數(shù);垂直分布;浙江

    大氣氣溶膠可直接或間接影響地球氣候系統(tǒng),一方面通過吸收和散射太陽輻射直接影響地氣系統(tǒng)的輻射收支平衡;還可作為云凝結核和冰核引起間接輻射強迫,繼而影響氣候變化.對流層氣溶膠垂直分布和光學特性對于研究氣溶膠對氣候的影響非常重要[1],其垂直分布特征更是評估氣溶膠輻射效應的關鍵因素之一.長江三角洲地區(qū)作為中國重要的經(jīng)濟文化中心,正面臨著嚴重的空氣污染問題[2-3],而浙江作為長江三角洲的一部分,研究其氣溶膠垂直結構對于探究長江三角洲的污染源和傳輸機制具有重要價值.

    微脈沖激光雷達(MPL)作為一種典型的主動遙感手段,可連續(xù)獲取氣溶膠在垂直方向上的信息[4].近年來,已有眾多學者利用微脈沖激光雷達在氣溶膠垂直分布方面開展了科學研究.丁輝[5]利用2008年5月14日~12月28日微脈沖激光數(shù)據(jù)研究了安徽壽縣地區(qū)氣溶膠消光系數(shù)垂直廓線的分布特征;Fan等[6]利用3a的激光雷達探測數(shù)據(jù)分析了長江三角洲地區(qū)氣溶膠的垂直分布特征,并結合潛在源貢獻函數(shù)和濃度加權軌跡模型探討了氣溶膠的來源;Sun等[7]利用微脈沖激光雷達和CE-318太陽光度計對2013~2015年長江三角洲上空氣溶膠直接輻射強迫(ADRF)與氣溶膠垂直結構的月變化及相互作用進行了研究,并結合中分辨率成像光譜儀(MODIS)與云氣溶膠激光雷達紅外探索觀測系統(tǒng)(CALIPSO)衛(wèi)星產(chǎn)品討論了2015年8月的一起生物質燃燒事件.然而,以往研究多局限于單個地面激光雷達站點,代表區(qū)域有限,難以獲得大范圍氣溶膠的垂直分布特征.

    氣溶膠來源與形成機制錯綜復雜,不同地區(qū)的氣溶膠受地理位置、氣象條件和人類活動等因素影響,其水平與垂直分布呈顯著差異[8-10].因此,本文選取位于浙江省的6個激光雷達站點數(shù)據(jù),首先對2017~2018年垂直方向上的氣溶膠消光系數(shù)進行多站點、長時間序列的連續(xù)觀測與數(shù)據(jù)反演,旨在探究不同站點的空間差異性,為大尺度范圍的氣溶膠垂直分布研究提供重要依據(jù).其次,為進一步了解長江三角洲地區(qū)的污染傳輸機制,選取2018年4月24~28日和7月25~29日浙江省兩次典型高空氣溶膠傳輸過程,結合地基激光雷達數(shù)據(jù),MODIS氣溶膠光學厚度(AOD)產(chǎn)品和火點資料以及HYSPLIT4后向軌跡模型對兩例事件中氣溶膠空間分布情況與污染來源進行綜合分析.最后,利用CALIPSO L2產(chǎn)品對垂直方向上的氣溶膠進行識別與分類,并對污染輸送路線進行探討,有助于加深了解浙江省顆粒物來源,以期為長江三角洲跨區(qū)域污染物的聯(lián)合防控提供一定借鑒與參考.

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)域

    圖1 浙江省激光雷達站點的地理位置

    本文使用的激光雷達分別位于浙江省的嘉興市南湖區(qū)善西站點(30.82°N,120.87°E);金華市塘雅中學站點(29.16°N,119.80°E);杭州市富陽區(qū)漁山鄉(xiāng)(30.08°N,120.09°E);杭州市淳安縣千島湖站點(29.61°N,119.03 °E);寧波市奉化區(qū)滕頭村站點(29.704°N,121.38°E);諸暨市陳蔡水庫站點(29.58°N, 120.39°E).

    浙江省位于我國的東南沿海?長江三角洲南翼,氣候較為溫和濕潤.地形自西南向東北方向呈現(xiàn)階梯狀變化,西南多為山地,中部以丘陵為主,東北部是沖積平原.圖1給出了激光雷達站點的具體空間位置分布.其中千島湖與陳蔡水庫站點眥鄰水庫,生態(tài)環(huán)境良好,車輛較少且樹木繁多,無明顯排放源;而塘雅站點位于塘雅中學附近,人口密集,人類活動較為頻繁.

    1.2 地基激光雷達數(shù)據(jù)

    本文所使用的氣溶膠垂直結構觀測儀器為微脈沖激光雷達(Micro Pulse LiDAR, MPL),其激光發(fā)射器的工作波長為532nm,空間垂直分辨率為30m,時間分辨率為1min,最大探測高度可達45km,探測盲區(qū)約為270m.選取數(shù)據(jù)時段為2017~2018年,對原始測量數(shù)據(jù)取小時平均值以提高信噪比,根據(jù)消光系數(shù)大于2km-1和退偏比的大小刪除云?雨時段[11-12].

    由于激光雷達接收的信號是能量值,通常通過求解雷達方程獲得光學特性值. 在對激光雷達信號處理前應先對相應的雷達因子進行修正,修正后的雷達方程如下:

    式中:()是激光雷達接收到高度(km)處的大氣后向散射回波功率, W;是激光雷達發(fā)射的能量, μJ;代表激光雷達系統(tǒng)常數(shù), W·km3·sr;是高度變量,范圍從0到r;()和()分別是在高度(km)處的大氣總后向散射系數(shù)(SR/km)和消光系數(shù)(km?1).

    由于散射系數(shù)和消光系數(shù)都為未知數(shù),因此需要選用一定的反演方法.本文采用 Fernald 算法[13]反演激光雷達信號.Fernald方法考慮了氣溶膠和空氣分子,激光雷達方程因此可以表示為:

    式中:下標a和m分別代表氣溶膠和大氣分子.定義大氣氣溶膠的消光后散射比即激光雷達比a()=a()/a(),空氣分子的消光后散射比m()=m()/m().基于Fernald算法,假設m(km)高度不存在氣溶膠,得到處的氣溶膠消光系數(shù)為:

    激光雷達比a假設為不隨高度變化的常數(shù)[14],變化范圍通常在20~70之間,這里取為50;空氣分子的消光后散射比m根據(jù)瑞利散射理論得出為(8π/3)sr.

    氣溶膠光學厚度(AOD)可由消光系數(shù)積分得到,將激光達的探測盲區(qū)(270m)等同于地表,計算公式如下:

    1.3 星載激光雷達數(shù)據(jù)

    2006年發(fā)射的云氣溶膠激光雷達紅外探索觀測系統(tǒng)(CALIPSO)是A-Train衛(wèi)星編隊的成員之一.其搭載的主要探測設備CALIOP是一臺偏振雙波長米散射激光雷達,可輸出波長為532和1064nm的脈沖.CALIOP通過接收大氣中不同層的后向散射,能夠較為準確地測量氣溶膠和云的光學性質及其形態(tài)的垂直分布廓線.本文使用的數(shù)據(jù)為 CALIPSO L2產(chǎn)品中4.20版本的垂直結構特性(VFM)和氣溶膠分類(AS)圖像,這些產(chǎn)品可從美國蘭利研究中心(NASA Langley Research Center https://www-CALIPSO.larc.nasa. gov/)獲得.

    1.4 MODIS衛(wèi)星數(shù)據(jù)

    中分辨率成像光譜儀(MODIS)是搭載在Terra和Aqua衛(wèi)星上的重要傳感器之一.MODIS衛(wèi)星每隔1~2d實現(xiàn)對整個地球表面的重復觀測,可收集范圍從0.4~14.4μm之間共36個波段的光譜輻照度數(shù)據(jù),以分辨率高?覆蓋范圍廣等優(yōu)點被廣泛運用于氣溶膠監(jiān)測.為了獲取氣溶膠在中國東部的空間分布,選擇MODIS L2中Terra 和Aqua C6.1版本中融合暗像元(dark target) 和深藍算法(deep blue)的10km×10km分辨率的550nm波段AOD數(shù)據(jù)[15-16],時間范圍分別是2018年4月25~ 29日和7月24~28日.與此同時,選用MODIS C6的火點產(chǎn)品(https://firms.modaps.eosdis.nasa. gov/)了解2018年4月25~29日和7月24~28日中國東部地區(qū)生物質燃燒的分布情況.

    1.5 污染物濃度數(shù)據(jù)

    本文使用的PM2.5、PM10以及空氣質量指數(shù)(AQI)來自中國國家環(huán)境監(jiān)測中心(http://106.37.208. 233:20035/).選取與各激光雷達站點距離最近的環(huán)境監(jiān)測站點數(shù)據(jù),以此代表激光雷達站點的污染物濃度水平.

    1.6 后向軌跡與潛在源分析

    為了明確浙江省污染氣團的來源與傳輸機制,通過基于美國國家海洋和大氣管理局(NOAA) HYSPLIT4模型的TrajStat軟件(http://www.arl.noaa. gov/ready/hysplit4.html),以6個激光雷達站點的位置作為起點,起始高度選擇兩次污染事件中高空氣溶膠層中心所在高度(1400和1500m),以NECP-FNL (國家環(huán)境預測中心,Final Analysis)數(shù)據(jù)作為氣象數(shù)據(jù)輸入,每隔6h (00:00,06:00,12:00,18:00UTC)對浙江省2次外來傳輸事件的高空氣團來源進行72h后向軌跡分析.基于后向軌跡結果[17],結合PM2.5小時濃度數(shù)據(jù),進一步利用潛在源貢獻函數(shù)(PSCF)和濃度加權軌跡(CWT)模型確定污染物的來源與位置,將PM2.5濃度等于35μg/m3作為污染閾值,具體算法與模型參數(shù)介紹可參考文獻[6].

    2 結果與分析

    2.1 年均氣溶膠垂直分布特征

    氣溶膠消光系數(shù)表示氣溶膠粒子通過散射和吸收作用,對某一波段太陽輻射單位面積?單位長度的削減程度.消光系數(shù)的大小可間接反映大氣污染程度:通常,排除云等因素影響,氣溶膠消光系數(shù)越小,表明污染越輕,反之,污染重.區(qū)別于先前研究只采用單個激光雷達代表整個研究區(qū)域氣溶膠垂直分布的局限,采用同型號多站點激光雷達細化浙江省不同區(qū)域的氣溶膠垂直分布的空間差異性,如圖2、圖3所示.

    圖2 由微脈沖激光雷達(MPL)反演得到的年均氣溶膠消光系數(shù)垂直廓線

    水平誤差棒對應每個站點的標準差

    圖3 由微脈沖激光雷達(MPL)反演得到的2017~2018年年均氣溶膠消光系數(shù)垂直廓線

    圖2表示浙江省激光雷達站點2017, 2018年年均氣溶膠消光系數(shù)的垂直分布特征.不難看出,大部分站點消光系數(shù)在2018年整體小于2017年且呈相似的垂直廓線形態(tài).而在漁山站,0.5km以下2018年消光系數(shù)整體較高;0.5~3km,2017年消光系數(shù)整體較高;以1.3km為界,陳蔡水庫站點呈現(xiàn)出與之相反的變化趨勢.此外,塘雅和漁山站消光系數(shù)年均值最高可達0.58km-1和0.61km-1,明顯大于其它站點,忽略相對濕度對顆粒物消光的影響,可間接反映出塘雅和漁山站點近地面污染情況較為嚴重.圖3展現(xiàn)的是各激光雷達站點2017~2018氣溶膠消光系數(shù)年均值的垂直分布情況.邊界層內(大概在1km以下)氣溶膠消光系數(shù)年均值由高到低依次為:塘雅>漁山>滕頭村>善西>陳蔡水庫>千島湖,空間分布整體為中部高,東西低,地形地貌和經(jīng)濟活動強度是近地面氣溶膠消光系數(shù)存在顯著空間差異的主要原因:塘雅站點位于塘雅中學附近,人口密集,人為氣溶膠的大量排放造成近地面消光系數(shù)較高;漁山站位于杭州市,為浙江省經(jīng)濟、金融中心,人類活動頻繁,污染排放較為嚴重;陳蔡水庫和千島湖站點生態(tài)環(huán)境良好并臨近水域,溫和的氣候條件更有利于去除大氣顆粒物,且無明顯排放源,因此底層消光較低.

    2.2 季均氣溶膠垂直分布特征

    為進一步研究浙江省氣溶膠的季節(jié)垂直分布特征,以月均值作為時間變化規(guī)律基礎,按照3~5月,6~8月,9~11月,12月~次年2月劃分為春季、夏季、秋季和冬季.

    激光雷達站點氣溶膠消光系數(shù)垂直分布和大氣邊界層高度的季節(jié)特征如圖4所示,其中大氣邊界層高度采用梯度法[18]反演得到.夏季邊界層高度達到峰值,可超過1.2km,春秋次之,冬季邊界層高度僅維持在0.5km左右,太陽輻射減少加之靜穩(wěn)天氣居多抑制了大氣邊界層發(fā)展.對比各站點的四季消光系數(shù)垂直廓線,發(fā)現(xiàn)消光系數(shù)大多集中在1km以下,與之前的研究結論一致[7,19].這可能是由于長江三角洲地區(qū)工業(yè)發(fā)達,加之強烈的人為活動導致氣溶膠主要積聚于對流層低空大氣中[19].氣溶膠消光系數(shù)季節(jié)變化范圍存在明顯差異,冬季近地面的消光系數(shù)顯著高于其他季節(jié),消光系數(shù)整體表現(xiàn)為冬季最高,秋季其次,而春夏較低.秋冬季節(jié)由于氣溫較低,相對穩(wěn)定的天氣形勢導致低層的氣溶膠層消光更強,氣溶膠更易集中在近地面[20].在1~3km高度上,氣溶膠消光系數(shù)在冬季快速衰減并趨近于0;夏季,觀察到各站點消光廓線在約1.5km高度出現(xiàn)明顯凸起,春秋季峰值稍弱,其中以滕頭村和陳蔡水庫站尤為明顯.沈吉等[19]研究發(fā)現(xiàn)長江三角洲地區(qū)6~9月大氣垂直分布不均,大氣垂直運動劇烈使夏秋季對流層高空出現(xiàn)多處較大峰值.綜上所述,出現(xiàn)高空氣溶膠峰值可能是本地排放或其他污染地區(qū)輸送所致,具體原因值得進一步探究.

    總體而言,氣溶膠集中分布在2km以下(圖5),最低占78.79%,最高達94.80%.具體而言,在春?夏兩季,除善西站點,氣溶膠多集中在地面上方1~2km處,占比超40%,最高至50.25%,這與夏季對流旺盛,強烈的垂直混合使更多的氣溶膠從邊界層被帶到自由對流層有關[21-22];其次分布在1km以下,占比達42.17%.秋季,6個站點呈現(xiàn)出2種不同的空間分布模式:對于善西?塘雅和漁山站點,氣溶膠在1km以下比率明顯超過1~2km,可達60.82%,氣溶膠在近地面出現(xiàn)聚集情況;相反,對于千島湖,滕頭村和陳蔡水庫站點,氣溶膠在1~2km的占比大于1km以下.在冬季, 1km以下的氣溶膠比例遠超1~2和2~3km,這可能是冬季多靜穩(wěn)天氣,較低的大氣邊界層和太陽輻射加之大氣湍流活動減弱,導致污染物被困在邊界層內難以擴散和清除

    圖4 由微脈沖激光雷達(MPL)反演得到的2017~2018年季均氣溶膠消光系數(shù)垂直廓線及邊界層高度

    水平誤差棒對應每個站點的標準差

    圖5 各季節(jié)激光雷達站點不同高度內AOD與總AOD的比值

    從圖6可見,各站點顆粒物濃度具有如下特征:首先,PM2.5和PM10質量濃度呈現(xiàn)一致的季節(jié)變化規(guī)律,從高到低依次為冬季,春季,秋季,夏季.冬季各站點顆粒物濃度水平顯著高于其他季節(jié),這主要歸因于供暖期燃煤等能源消耗的增加使人為顆粒物排放顯著上升;春季受風沙等污染源因素影響使顆粒物尤其是PM10濃度普遍較高;夏季由于大氣邊界層較高且溫度升高,垂直方向的空氣對流較強,為顆粒在垂直方向上擴散和稀釋提供有利條件,PM濃度水平較低.

    PM2.5/PM10比值較PM10和PM2.5更能準確地表示二次污染程度,可作為判斷大氣污染程度和污染源的依據(jù):PM2.5在PM10中所占比重越高,污染越嚴重[23].許多研究將PM2.5/PM10比值0.5作為區(qū)分自然污染和人為污染的分界線[24-25].圖6各站點PM2.5/ PM10始終保持在51%以上,最高可達76.50%,反映出浙江省局地細顆粒污染嚴重.冬季,PM2.5/PM10在65.71%~76.50%間變化,表明污染氣團以細顆粒為主導,且最高值出現(xiàn)在塘雅站點,該地人口密集,人為污染嚴重,故而細顆粒物富集.此外,漁山、滕頭村、塘雅、善西站點PM2.5/PM10在春夏秋變化較為平穩(wěn)且以夏季為谷值;千島湖、陳蔡水庫的變化趨勢與之不同,PM2.5/PM10在夏季轉為峰值,達63%,細顆粒物占比僅次于冬季,推測產(chǎn)生這種空間變化的原因為:一方面,氣象因素通過一系列相互矛盾的過程影響PM濃度,最終取決于當?shù)貧夂?地形和PM成分的綜合作用結果[26-27].考慮到浙江省各地區(qū)經(jīng)濟發(fā)展程度不均,污染排放類型也存在相應差異.另一方面,有研究表明PM10與氣象因子的相關系數(shù)高于PM2.5與氣象因子的相關系數(shù),即氣象條件對PM10的影響更顯著:當氣象條件轉好時,有助于顆粒物沉降或擴散且粗粒子更易被清除[28].結合圖6發(fā)現(xiàn)PM10在夏季的下降幅度明顯大于PM2.5,這也可能是千島湖?陳蔡水庫在夏季細顆粒物占比上升的成因之一.

    圖6 2017~2018年各站點PM2.5和PM10的質量濃度及其比值的季均變化

    2.3 氣溶膠跨區(qū)域傳輸事件分析

    在區(qū)域污染源排放相對穩(wěn)定的情況下,污染物濃度的高低主要取決于大氣擴散能力,區(qū)域氣象條件的變化會直接影響污染物累積以及重污染天氣的形成、發(fā)展和消散[29].大量研究表明,長江三角洲區(qū)域性污染時常發(fā)生,除了受到區(qū)域內部污染的交互影響外,京津冀魯豫等長江三角洲以北城市群的跨區(qū)域污染輸送也不容忽視[30-31].

    2018年4月24~28日以及7月25~29日在浙江省上空發(fā)生了兩起氣溶膠跨區(qū)域傳輸事件,多個激光雷達站點同時觀測到高空氣溶膠層存在,如圖7所示.滕頭村站在4月25日早晨06:00左右最先觀測到高空氣溶膠;除漁山站由于噪聲影響導致數(shù)據(jù)缺失外,其余站點于25日08:00前后出現(xiàn)高空氣溶膠層.對比各站點時間-高度的消光系數(shù)變化可知,該氣溶膠層位于1000~1800m高度且明顯與地面氣溶膠層相脫離,可知為高空氣溶膠層;善西和滕頭村對應高度上退偏比大于0.12,說明高空顆粒物不規(guī)則程度較高,可能是混入部分沙塵所致;其余站點退偏比在0.06~0.1之間,不規(guī)則程度相對較低.隨時間推移,1000m以上氣溶膠消光系數(shù)顯著增大,超過2km-1.近地面層消光系數(shù)在25日白天始終高于0.5km-1,鑒于此時高空氣溶膠還未發(fā)生沉降,因此推測為局地排放引發(fā)污染所致;25日20:00過后,各站點近地面消光系數(shù)逐漸減小,污染程度有所緩解;27日高空氣溶膠層高度不斷下降,19:00開始降至1000m以下并逐步與大氣邊界層內氣溶膠相混合; 28日00:00前后,隨著空中氣溶膠層與近地面氣溶膠層混合均勻,底層消光系數(shù)急劇上升至1km-1以上,與之對應的高空消光系數(shù)明顯減小,近地面退偏比受沉降顆粒物影響略有上升,污染傳輸過程至此趨于結束.

    圖8描述了2018年4月24~28日期間PM2.5, PM10質量濃度和空氣質量指數(shù)(AQI)的變化情況.依據(jù)主要污染來源不同可劃分為兩階段:24~25日為第一階段,以本地污染排放為主引發(fā)重度污染;26日為短暫間歇期;27~28日為第二階段,高空氣溶膠沉降導致污染程度再次加重并達到嚴重污染.4月24~25日上午,浙江省各站點PM質量濃度持續(xù)攀升并于25日陸續(xù)達到峰值,其中善西站點的PM2.5和PM10濃度高達232和297μg/m3;各站點AQI也不斷上升并超過200,最高達282,造成重度污染.25日12:00后PM2.5和PM10濃度開始快速回落,26日PM2.5和PM10濃度在50~100μg/m3之間震蕩,表明污染進入間歇期,27日凌晨前后,各站點PM濃度開始大幅增加并在27日夜間至28日上午陸續(xù)達到峰值,這與氣溶膠沉降時間吻合.值得注意的是,在善西、滕頭村和陳蔡水庫站點,由跨區(qū)域傳輸引發(fā)的第二次污染強度明顯超過第一階段因區(qū)域內部污染累積而導致的污染,PM2.5和PM10濃度甚至高達271和322μg/m3,說明高空氣溶膠沉降引起污染物再次累積,是導致污染程度更為惡化的主要成因.

    激光雷達在2018年7月25~29日觀測到浙江上空有一段連續(xù)的污染物輸送過程(漁山站對應觀測時段數(shù)據(jù)缺失),高空氣溶膠所在高度為1000~ 2000m,退偏比在0.04~0.08之間,顆粒物較為規(guī)則. 27~29日高空氣溶膠層消光系數(shù)明顯上升,污染輸送強度呈增加趨勢,其中千島湖站點的傳輸強度明顯弱于其余站點.得益于夏季大氣邊界層較高,良好的氣象條件利于污染物擴散,盡管高空氣溶膠發(fā)生了沉降,但環(huán)保監(jiān)測站點數(shù)據(jù)顯示在此期間PM2.5和PM10濃度始終保持在75μg/m3以下,高空粒子的輸送并未明顯引發(fā)當?shù)匚廴?

    2018年4月24~28日及7月25~29日MODIS氣溶膠光學厚度以及火點位置如圖9a和圖9e所示.4月24~28日中國東部地區(qū)多為AOD高值區(qū)且與火點分布位置大致吻合,其中江蘇,安徽和浙江北部地區(qū)火點分布較為密集(圖9a),說明這些區(qū)域存在大量的生物質燃燒氣溶膠.7月25~29日浙江,安徽東部,山東和河南東北部表現(xiàn)為AOD高值區(qū);火點集中分布在江蘇,浙江和安徽北部地區(qū)的同時,還廣泛分布在浙江以南地區(qū).

    圖8 4月24~28日PM2.5、PM10質量濃度和AQI的逐小時變化

    圖9 2018年4月24~28日和7月25~29日AOD和火點位置、72h后向軌跡聚類分析、WPSCF以及WCWT的空間分布

    為進一步明確兩次氣溶膠跨區(qū)域傳輸事件的污染來源,利用HYSPLIT4模式分別計算2018年4月24~28日以及7月25~29日的72h后向軌跡,時間間隔為6h,模擬起始高度選取為海拔高度1400和1500m,分別對應兩次傳輸事件中高空氣溶膠層中心位置所在高度.依據(jù)聚類分析結果,4月的污染輸送事件主要被分為4軌跡,如圖9b所示.從水平輸送方向來看,來自山西、河北方向的軌跡1占總軌跡比重最大,為44.17%,來自湖北、江西方向的軌跡3次之,占23.33%,說明浙江省春季主要受大陸氣團控制;從垂直方向看,軌跡1和4從2800m以上的高度不斷沉降至1400m,為下沉氣團,軌跡2和3的氣團表現(xiàn)出先上升后沉降的過程.為降低計算誤差,將權重函數(shù)與潛在源貢獻(PSCF)模型和濃度加權軌跡模型(WCWT)相結合,得出加權后的WPSCF和WCWT結果[32-35].如圖9c~d所示,潛在源貢獻區(qū)域與PM2.5濃度加權軌跡貢獻呈現(xiàn)出較好的一致性,潛在源區(qū)域集中在江浙滬以及湖北和江西的交界地帶, WCWT值均高于75μg/m3,表明相應地區(qū)對浙江省此次污染輸送PM2.5貢獻最大.7月傳輸事件的72h后向軌跡聚類成4類,如圖9f所示.水平輸送方向看,軌跡大部分來自浙江西南方向:江西安徽方向的軌跡2占據(jù)比重最高,達40%,其次是來自福建江西方向的軌跡1;從垂直方向看,3類軌跡均呈先升后降態(tài)勢.由圖9g~h可見,潛在源貢獻與PM2.5濃度加權軌跡貢獻區(qū)域存在一定差異,對污染物輸送貢獻最高的是江西北部以及安徽和江蘇南部地區(qū),PM2.5貢獻超40μg/m3.

    圖10 氣溶膠類型隨高度的分布情況

    上圖代表中國范圍內氣溶膠類型隨高度的分布情況,虛線框為浙江省范圍,下圖為浙江省局部放大圖;1.潔凈海洋型;2.沙塵;3.污染大陸型/煙;4.潔凈大陸型;5.污染沙塵型;6.煙塵型;7.污染海洋型

    為研究兩次跨區(qū)域傳輸氣溶膠在垂直方向上的分布特征,選取中國范圍內CALIPSO衛(wèi)星穿過浙江上空時的L2產(chǎn)品,對應北京時間分別為2018年4月28日13:07~13:21,7月26日02:02~02:16和7月28日01:50~02:04.CALIPSO L2產(chǎn)品通過不同高度的后向散射數(shù)據(jù)推導出532nm后向散射系數(shù)、體積退偏比、區(qū)域類型以及抬升現(xiàn)象等特征雷達信號,對大氣層內不同高度的氣溶膠類型進行區(qū)分[36],利用CALIPSO的AS產(chǎn)品可清晰地揭示氣溶膠在不同垂直高度上的類別.浙江省AS分布如圖10虛線框及下圖所示,結果顯示在2018年4月28日(圖10a),浙江以北地區(qū)諸如內蒙、遼寧地區(qū)出現(xiàn)大范圍沙塵,高度從1km蔓延至6km以上;浙江區(qū)域緯度內,在1.4km以下表現(xiàn)為潔凈海洋型與污染海洋型混雜,1.4~2km以污染大陸型/煙為主,其次為污染海洋型和沙漠沙塵型;浙江以南地區(qū)在1km以下多為潔凈海洋型而高空夾帶部分煙塵型氣溶膠.結合HYSPLIT模型的72h后向軌跡及AOD和火點空間分布(圖9),可看出2018年4月24~28日以西北方向為主的氣流(軌跡1)經(jīng)過河北、山東和浙江等AOD高值區(qū),易將這些地區(qū)的大氣污染物諸如生物質燃燒氣溶膠輸送至浙江上空,加重本地污染.春季,東亞沙塵事件頻發(fā),中國西北地區(qū)由于戈壁沙漠和黃土分布廣泛,是亞洲揚塵的主要來源地區(qū)之一,起源于塔克拉瑪干和戈壁沙漠的沙塵氣溶膠可長距離運輸?shù)街袊鴸|部,韓國和日本[37-39].因此,外來沙塵極可能在高空西北氣流的輸送下到達浙江上空,沉降后對浙江PM10濃度增加具有重要貢獻.7月26~28日,蘇魯?shù)貛С霈F(xiàn)大量煙塵型氣溶膠(圖10b),與AOD高值區(qū)和火點密集區(qū)分布相吻合,隨時間推移同時向南北擴張至浙江、吉林和遼寧.26日,浙江省內1~2km的氣溶膠為污染大陸型/煙,28日轉變?yōu)闊焿m型,其次為潔凈海洋型和污染大陸型/煙(圖10c).值得注意的是,7月25~29日氣流均來自浙江的西南方向,盡管這些地區(qū)屬于AOD低值區(qū)但存在大量火點,生物質燃燒氣溶膠隨西南氣流不斷輸送至研究區(qū)域,導致其上空出現(xiàn)大規(guī)模煙塵型氣溶膠.生物質燃燒作為大氣污染的主要來源之一,不僅能夠惡化區(qū)域空氣質量、降低大氣能見度,還可顯著影響生態(tài)系統(tǒng)的循環(huán),對人類健康和環(huán)境產(chǎn)生不利效應.

    3 結論

    3.1 大部分站點2018年氣溶膠消光系數(shù)整體小于2017年.邊界層內消光系數(shù)年均值表現(xiàn)為塘雅>漁山>滕頭村>善西>陳蔡水庫>千島湖,在空間上呈現(xiàn)中部高而東西低的分布形態(tài),這主要歸因于地形地貌以及經(jīng)濟活動強度差異.

    3.2 氣溶膠消光系數(shù)存在明顯季節(jié)變化,整體表現(xiàn)為冬季最高,秋季其次,而春夏較低.從垂直分布上看,春季和夏季大部分站點的氣溶膠集中在地面上方1~2km處,冬季則分布在1km以下.秋季,千島湖,滕頭村和陳蔡水庫站點氣溶膠在1~2km占比更多,高空氣溶膠的存在使垂直分布形態(tài)發(fā)生改變.各站點PM2.5和PM10質量濃度的季節(jié)變化規(guī)律一致表現(xiàn)為冬季>春季>秋季>夏季,浙江省細粒子污染嚴重且受人為排放影響較高.

    3.3 2018年4月24~28日,高空氣溶膠沉降促使局地排放引發(fā)污染加重,人為污染產(chǎn)生的細粒子占主要地位,西北方向氣流經(jīng)AOD高值區(qū)時易將人為排放氣溶膠攜帶至浙江本地,污染強度明顯增加.7月25~29日,西南方向氣流途經(jīng)江西、福建等廣泛存在火點的地區(qū),導致浙江上空出現(xiàn)大規(guī)模煙塵型氣溶膠.

    [1] Kaufman Y J, Tanre D, Gordon H R, et al. Passive remote sensing of tropospheric aerosol and atmospheric correction for the aerosol effect [J]. Journal of Geophysical Research, 1997,102(14):16815-16830.

    [2] Cai Q L, Dai X R, Li J R, et al. The characteristics and mixing states of PM2.5during a winter dust storm in Ningbo of the Yangtze River Delta, China [J]. Science of the Total Environment, 2019,709:136- 146.

    [3] Du W, Zhang Y, Chen Y. Chemical characterization and source apportionment of PM2.5during spring and winter in the Yangtze River Delta, China [J]. Aerosol and Air Quality Research, 2017,17(9):2165- 2180.

    [4] 賀千山,毛節(jié)泰.微脈沖激光雷達及其應用研究進展 [J]. 氣象科技, 2004,32(4):219-224. He Q S, Mao J T. Micro-pulse lidar and its applications [J]. Meteorological Science and Technology, 2004,32(4):219-224.

    [5] 丁 輝.利用微脈沖激光雷達(MPL)探測氣溶膠消光系數(shù)廓線和大氣混合層高度的初步研究 [D]. 南京:南京信息工程大學, 2012. Ding H. Measurements of aerosol vertical profiles and the mixed layer height using a micro pulse lidar [D]. Nanjing: Nanjing University of Information Science & Technology, 2012.

    [6] Fan W, Qin K, Xu J, et al. Aerosol vertical distribution and sources estimation at a site of the Yangtze River Delta region of China [J]. Atmospheric Research, 2019,217:128-136.

    [7] Sun T, Che H, Qi B, et al. Characterization of vertical distribution and radiative forcing of ambient aerosol over the Yangtze River Delta during 2013~2015 [J]. Science of the Total Environment, 2019,650(2): 1846-1857.

    [8] 李正強,謝一凇,張 瑩,等.大氣氣溶膠成分遙感研究進展 [J]. 遙感學報, 2019,23(3):359-373. Li Z Q, Xie Y S, Zhang Y, et al. Advance in the remote sensing of atmospheric aerosol composition [J]. Journal of Remote Sensing, 2019,23(3):359-373.

    [9] 邱 昀,安欣欣,劉保獻,等.北京市氣溶膠消光系數(shù)垂直特征及影響因子探討[J]. 環(huán)境科學研究, 2020,33(3):519-525. Qiu Y, An X X, Liu B X, et al. Vertical distribution of aerosol extinction coefficient and its influencing factor in Beijing [J]. Research of Environmental Sciences, 2020,33(3):519-525.

    [10] 吳 儀,鄧孺孺,秦 雁,等.秋冬季節(jié)珠三角典型區(qū)域氣溶膠濃度垂直分布的激光雷達反演 [J]. 遙感學報, 2020,24(3):302-318. Wu Y, Deng R R, Qin Y, et al. Vertical distribution of aerosol mass concentration over Pearl River Delta observed by Lidar during autumn and winter [J]. Journal of Remote Sensing, 2020,24(3):302-318.

    [11] Zhao C, Wang Y, Wang Q, et al. A new cloud and aerosol layer detection method based on micropulse lidar measurements [J]. Journal of Geophysical Research: Atmospheres, 2014,119(11):6788-6802.

    [12] Zhang J L, Liu P, Zhang F, et al. CloudNet: Ground‐based cloud classification with deep convolutional neural network [J]. Geophysical Research Letters, 2018,45(16):8665-8672.

    [13] Fernald F G. Analysis of atmospheric lidar observations: some comments [J]. Applied Optics, 1984,23(5):652-652.

    [14] Bucci S, Cristofanelli P, Decesari S, et al. Vertical distribution of aerosol optical properties in the Po Valley during the 2012 summer campaigns [J]. Atmospheric Chemistry & Physics, 2018,18(8):5371- 5389.

    [15] Levy R C, Remer L A, Mattoo S, et al. Second-generation operational algorithm: Retrieval of aerosol properties over land from inversion of Moderate Resolution Imaging Spectroradiometer spectral reflectance [J]. Journal of Geophysical Research-Atmospheres, 2007,112(D13): 319-321.

    [16] Levy R C, Remer L A, Dubovik O. Gobal aerosol optical models and lookup tables for the new MODIS aerosol retrieval over land [Z]. Journal of Geophysical Research-Atmospheres, 2006.

    [17] Wang Y Q, Zhang X Y, Draxler R R. TrajStat: GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data [J]. Environmental Modelling & Software, 2009,24(8):938-939.

    [18] 張艷品,陳 靜,鈐偉妙,等.石家莊冬季典型污染過程氣溶膠激光雷達觀測 [J]. 中國環(huán)境科學, 2020,40(10):4205-4215. Zhang Y P, Chen J, Qian W M, et al. Aerosol lidar observation of typical pollution process in Shijiazhuang in Winter [J]. China Environmental Science, 2020,40(10):4205-4215.

    [19] 沈 吉,曹念文.2018年中國長江三角洲地區(qū)氣溶膠的垂直分布特征 [J]. 環(huán)境科學, 2019,40(11):4743-4754. Shen J, Cao N W. Characteristics of aerosol vertical distribution over the Yangtze River Delta region of China in 2018 [J]. Environmental Science, 2019,40(11):4743-4754.

    [20] Leng C, Duan J, Xu C, et al. Insights into a historic severe haze weather in Shanghai: synoptic situation, boundary layer and pollutants [J]. Atmospheric Chemistry and Physics Discussions, 2015,15(22): 32561-32605.

    [21] 高星星,陳 艷,張 鐳,等.華北地區(qū)氣溶膠的季節(jié)垂直分布特征及其光學特性 [J]. 蘭州大學學報(自然科學版), 2018,54(3):395-403. Gao X X, Chen Y, Zhang L, et al. Vertical distribution of seasonal aerosols and their optical properties over Northern China [J]. Journal of Lanzhou University: Natural Sciences, 2018,54(3):395-403.

    [22] 唐志偉,許瀟鋒,楊曉玥,等.基于衛(wèi)星資料的華東地區(qū)氣溶膠三維分布特征研究[J]. 中國環(huán)境科學, 2019,39(9):3624-3634. Tang Z W, Xu X F, Yang X Y, et al. Characteristics of three- dimensional aerosol distribution in Eastern China based on satellite data [J]. China Environmental Science, 2019,39(9):3624-3634.

    [23] 楊文濤,姚詩琪,鄧 敏,等.北京市PM2.5時空分布特征及其與PM10關系的時空變異特征[J]. 環(huán)境科學, 2018,39(2):684-690. Yang W T, Yao S Q, Deng M, et al. Spatio-temporal distribution characteristics of PM2.5and spatio-temporal variation characteristics of the relationship between PM2.5and PM10in Beijing [J]. Environmental Science, 2018,39(2):684-690.

    [24] 鄧發(fā)榮,康 娜,Kanike R K,等.長江三角洲地區(qū)大氣污染過程分析 [J]. 中國環(huán)境科學, 2018,38(2):401-411. Deng F R, Kang N, Kanike R K, et al. Analysis of air pollution episodes over different cites in the Yangtze River Delta [J]. China Environmental Science, 2018,38(2):401-411.

    [25] 李 敏,何紅弟,郝楊楊.上海市大氣環(huán)境中PM2.5/PM10時空分布特征[J]. 云南大學學報(自然科學版), 2019,41(2):323-332. Gehrig R, Buchmann B. Characterising seasonal variations and spatial distribution of ambient PM10and PM2.5concentrations based on long-term Swiss monitoring data [J]. Atmospheric Environment, 2003, 37(19):2571-2580.

    [26] Yang Q Q, Yuan Q Q, Li T W, et al. The relationships between PM2.5and meteorological factors in China: Seasonal and regional variations [J]. arXiv e-prints, 1708.06072.

    [27] 李 瑞,李 清,徐 健,等.秋冬季區(qū)域性大氣污染過程對長三角北部典型城市的影響 [J]. 環(huán)境科學, 2020,41(4):1520-1534. Li R, Li Q, Xu J, et al. Regional air pollution process in winter over the Yangtze River Delta and its influence on typical northern cities [J]. Environmental Science, 2020,41(4):1520-1534.

    [28] Hua Y, Cheng Z, Wang S, et al. Characteristics and source apportionment of PM2.5during a fall heavy haze episode in the Yangtze River Delta of China [J]. Atmospheric Environment, 2015,123:380- 391.

    [29] Tang L, Yu H, Ding A, et al. Regional contribution to PM1pollution during winter haze in Yangtze River Delta, China [J]. Science of the Total Environment, 2016,541:161-166.

    [30] Wang M, Cao C, Li G, et al. Analysis of a severe prolonged regional haze episode in the Yangtze River Delta, China [J]. Atmospheric Environment, 2015,102:112-121.

    [31] Seibert P, Kromp-Kolb H, Baltensperger U, et al. Trajectory analysis of high-alpine air pollution data [J]. Air Pollution Modeling and Its Application X, 1994,18:595-596.

    [32] Polissar A. The aerosol at Barrow, Alaska: long-term trends and source locations [J]. Atmospheric Environment, 1999,33(16):2441- 2458.

    [33] Hsu Y K, Holsen T M, Hopke P K. Comparison of hybrid receptor models to locate PCB sources in Chicago [J]. Atmospheric Environment, 2003,37(4):545-562.

    [34] 孫天澤.杭州地區(qū)氣溶膠光學—輻射特性垂直觀測研究 [D]. 北京:中國氣象科學研究院, 2018. Sun T Z. Vertical observation and analysis of aerosol optical properties and direct radiative forcing over Hangzhou [D]. Beijing: Chinese Academy of Meteorological Sciences, 2018.

    [35] 于彩霞,楊元建,鄧學良,等.基于CALIOP探測的合肥氣溶膠垂直分布特征[J]. 中國環(huán)境科學, 2017,37(5):1677-1683. Yu C X, Yang Y J, Deng X L, et al. Vertical distribution characteristics of aerosol optical properties on haze and clear day in Hefei based on CALIOP satellite measuring [J]. China Environmental Science, 2017, 37(5):1677-1683.

    [36] Uematsu M, Yoshikawa A, Muraki H, et al. Transport of mineral and anthropogenic aerosols during a Kosa event over East Asia [J]. Journal of Geophysical Research-Atmospheres, 2002,107(D7):AAC3-1- AAC3-7.

    [37] Sheng Z, Ramesh S. Aerosol and meteorological parameters associated with the intense dust event of 15 April 2015 over Beijing, China [J]. Remote Sensing, 2018,10(6):957.

    [38] Zhang Z, Huang J, Chen B, et al. Three-year continuous observation of pure and polluted dust aerosols over northwest China using the ground-based lidar and sun photometer data [J]. Journal of Geophysical Research Atmospheres, 2019,124(2):1118-1131.

    Vertical distribution and regional transport of aerosols in Zhejiang Province.

    ZHANG Yi-shu1, XU Da2, HE Qin1, LI Dong-hui1, QIN Kai1*

    (1.School of Environment and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China;2.Zhejiang Ecological and Environmental Monitoring Center, Hangzhou 310012, China)., 2021,41(5):2097~2107

    The interannual and seasonal changes of aerosol vertical distributions in Zhejiang Province from 2017 to 2018 were studied using Micro Pulse Lidar (MPL) measurements in six sites. Air quality data were combined to analyse the dominant particulate matter types in terms of seasonal differences. Aerosol extinction coefficients of most sites in 2018 were lower than that in 2017. Spatially, the extinction coefficient within the boundary layer was high in central Zhejiang and low in the eastern and western Zhejiang. Seasonally, the extinction coefficient was the highest in winter, second in autumn, and lower in spring and summer; aerosols were concentrated within 1~2km above the ground in spring and summer, and mainly distributed below 1km in winter; in autumn, affected by high-altitude aerosols, the six stations showed different vertical distribution patterns. Furthermore, the aerosol sources and transmission mechanisms of two trans-regional transport episodes in April and July 2018 were revealed, by using the CALIPSO and MODIS satellite data, as well as the backward trajectory, PSCF and CWT models. It was found thatair masses from high AOD areas can significantly increase local pollution in Zhejiang.

    micro pulse LiDAR;aerosol;extinction coefficient;vertical distribution;Zhejiang

    X513

    A

    1000-6923(2021)05-2097-11

    張亦舒(1995-),女,河南焦作人,中國礦業(yè)大學環(huán)境與測繪學院碩士研究生,主要從事大氣環(huán)境遙感研究.

    2020-10-09

    徐州市重點研發(fā)計劃(社會發(fā)展)項目(KC18225);國家自然科學基金資助項目(41975041)

    * 責任作者, 教授, qinkai@cumt.edu.cn

    猜你喜歡
    氣溶膠激光雷達高空
    手持激光雷達應用解決方案
    北京測繪(2022年5期)2022-11-22 06:57:43
    高空走繩
    法雷奧第二代SCALA?激光雷達
    汽車觀察(2021年8期)2021-09-01 10:12:41
    氣溶膠傳播之謎
    基于激光雷達通信的地面特征識別技術
    氣溶膠中210Po測定的不確定度評定
    高空纜車
    不要高空拋物!
    高空莫拋物
    基于激光雷達的多旋翼無人機室內定位與避障研究
    電子制作(2018年16期)2018-09-26 03:27:00
    最近最新免费中文字幕在线| 久久久精品欧美日韩精品| 男女之事视频高清在线观看| 日韩有码中文字幕| 黄片播放在线免费| 亚洲中文字幕一区二区三区有码在线看 | 又黄又粗又硬又大视频| 亚洲欧美日韩无卡精品| 亚洲无线在线观看| 日本vs欧美在线观看视频| 黄片播放在线免费| 成人av一区二区三区在线看| 成人三级做爰电影| 亚洲色图综合在线观看| 动漫黄色视频在线观看| 亚洲,欧美精品.| 国产精品 国内视频| 桃色一区二区三区在线观看| 丁香六月欧美| 男人舔女人的私密视频| 一个人免费在线观看的高清视频| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 亚洲中文字幕日韩| www.999成人在线观看| 亚洲精品久久国产高清桃花| 久久人妻av系列| 99国产精品一区二区三区| 乱人伦中国视频| 视频区欧美日本亚洲| 中文字幕人妻熟女乱码| 亚洲午夜理论影院| 91成人精品电影| 露出奶头的视频| bbb黄色大片| 国产99久久九九免费精品| 亚洲精品在线美女| www.999成人在线观看| 91国产中文字幕| av天堂久久9| 麻豆成人av在线观看| avwww免费| 亚洲精品中文字幕一二三四区| 99久久精品国产亚洲精品| 看片在线看免费视频| 999精品在线视频| 人人妻人人爽人人添夜夜欢视频| 精品熟女少妇八av免费久了| 好看av亚洲va欧美ⅴa在| 99国产极品粉嫩在线观看| 亚洲无线在线观看| 女人被狂操c到高潮| 亚洲一区二区三区不卡视频| 久久狼人影院| 99riav亚洲国产免费| 国产精品一区二区在线不卡| 老司机靠b影院| 免费高清视频大片| 精品国产一区二区久久| 搡老熟女国产l中国老女人| 亚洲伊人色综图| 国产一区二区在线av高清观看| 国产欧美日韩一区二区精品| 久久精品国产综合久久久| 亚洲av成人不卡在线观看播放网| 露出奶头的视频| 很黄的视频免费| 男男h啪啪无遮挡| 精品一区二区三区视频在线观看免费| 亚洲自拍偷在线| 免费在线观看完整版高清| 亚洲av第一区精品v没综合| 日韩有码中文字幕| 午夜日韩欧美国产| 一a级毛片在线观看| 69av精品久久久久久| www.精华液| 免费无遮挡裸体视频| 国产主播在线观看一区二区| 这个男人来自地球电影免费观看| 啦啦啦韩国在线观看视频| 欧美精品啪啪一区二区三区| 亚洲av成人一区二区三| 精品国产国语对白av| 不卡一级毛片| 国产高清激情床上av| or卡值多少钱| 一边摸一边抽搐一进一小说| 精品一区二区三区四区五区乱码| 午夜久久久在线观看| 怎么达到女性高潮| 久热爱精品视频在线9| 亚洲精品在线观看二区| 1024香蕉在线观看| 满18在线观看网站| 97超级碰碰碰精品色视频在线观看| 中文字幕久久专区| 亚洲欧美日韩另类电影网站| 好男人电影高清在线观看| 宅男免费午夜| 欧美激情高清一区二区三区| 人人妻人人澡人人看| 久久热在线av| 99国产精品一区二区三区| 999久久久国产精品视频| 日韩成人在线观看一区二区三区| 午夜福利,免费看| 亚洲精品中文字幕一二三四区| 好男人在线观看高清免费视频 | 狂野欧美激情性xxxx| 免费观看精品视频网站| 777久久人妻少妇嫩草av网站| 久久人人97超碰香蕉20202| 一进一出好大好爽视频| 国产成人精品久久二区二区91| 亚洲欧美激情综合另类| 91大片在线观看| 国产在线精品亚洲第一网站| 真人做人爱边吃奶动态| 成人特级黄色片久久久久久久| 黄色女人牲交| 69av精品久久久久久| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 国产一卡二卡三卡精品| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 中国美女看黄片| 欧美另类亚洲清纯唯美| 色尼玛亚洲综合影院| 亚洲av美国av| 免费在线观看影片大全网站| 视频在线观看一区二区三区| 国产精品一区二区在线不卡| 国产成人av教育| 亚洲国产毛片av蜜桃av| 韩国精品一区二区三区| 日韩成人在线观看一区二区三区| 亚洲在线自拍视频| 9色porny在线观看| 亚洲色图综合在线观看| 人妻丰满熟妇av一区二区三区| 91老司机精品| 欧美乱码精品一区二区三区| 丝袜美足系列| 欧美国产精品va在线观看不卡| 成人欧美大片| 曰老女人黄片| 精品少妇一区二区三区视频日本电影| 亚洲色图av天堂| 身体一侧抽搐| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品sss在线观看| 国产精品久久久av美女十八| 黄片小视频在线播放| 亚洲国产精品999在线| 少妇被粗大的猛进出69影院| 一进一出好大好爽视频| 亚洲九九香蕉| 一级片免费观看大全| 一区二区三区精品91| 国产亚洲欧美98| 国产亚洲欧美精品永久| 制服人妻中文乱码| 国产99久久九九免费精品| 丝袜美足系列| 欧美日韩瑟瑟在线播放| 亚洲国产精品sss在线观看| 亚洲欧美激情在线| 精品人妻1区二区| 男女做爰动态图高潮gif福利片 | 女生性感内裤真人,穿戴方法视频| 国内毛片毛片毛片毛片毛片| 嫁个100分男人电影在线观看| 后天国语完整版免费观看| 欧美黄色淫秽网站| 91大片在线观看| 亚洲最大成人中文| 一二三四在线观看免费中文在| 国产精品综合久久久久久久免费 | 18禁国产床啪视频网站| 欧美老熟妇乱子伦牲交| 亚洲中文av在线| 亚洲精品久久国产高清桃花| 美女 人体艺术 gogo| www.精华液| 久久 成人 亚洲| 天天一区二区日本电影三级 | 日韩高清综合在线| 国产熟女午夜一区二区三区| 国产区一区二久久| 午夜久久久久精精品| 91精品国产国语对白视频| 色婷婷久久久亚洲欧美| 亚洲三区欧美一区| 国产av在哪里看| 1024香蕉在线观看| 国内毛片毛片毛片毛片毛片| 两性夫妻黄色片| av片东京热男人的天堂| 中文字幕最新亚洲高清| 欧美一级a爱片免费观看看 | 天堂√8在线中文| 大陆偷拍与自拍| 又大又爽又粗| 美女 人体艺术 gogo| 亚洲欧洲精品一区二区精品久久久| 精品高清国产在线一区| 亚洲午夜精品一区,二区,三区| 欧美一区二区精品小视频在线| 亚洲成人精品中文字幕电影| 欧美激情 高清一区二区三区| 免费女性裸体啪啪无遮挡网站| 亚洲五月婷婷丁香| 午夜老司机福利片| 级片在线观看| 无限看片的www在线观看| 妹子高潮喷水视频| 一a级毛片在线观看| 99国产精品免费福利视频| 99国产综合亚洲精品| 日韩大码丰满熟妇| 狂野欧美激情性xxxx| 日韩大尺度精品在线看网址 | 淫秽高清视频在线观看| 一边摸一边抽搐一进一小说| 可以在线观看的亚洲视频| 亚洲欧美精品综合久久99| 窝窝影院91人妻| 国产精品免费视频内射| 亚洲成人精品中文字幕电影| 欧美日本中文国产一区发布| 亚洲精华国产精华精| 精品欧美一区二区三区在线| 日日摸夜夜添夜夜添小说| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲精品第一综合不卡| 国产一区二区激情短视频| 午夜福利视频1000在线观看 | 亚洲美女黄片视频| 色综合站精品国产| 十八禁人妻一区二区| 中出人妻视频一区二区| 精品久久久精品久久久| 日本免费一区二区三区高清不卡 | 91麻豆av在线| 搡老熟女国产l中国老女人| 亚洲全国av大片| 人人妻人人爽人人添夜夜欢视频| 久久久久久久精品吃奶| 中文字幕色久视频| 国产精品电影一区二区三区| 夜夜看夜夜爽夜夜摸| 久久精品国产亚洲av香蕉五月| 香蕉丝袜av| 99久久精品国产亚洲精品| 日日摸夜夜添夜夜添小说| 变态另类成人亚洲欧美熟女 | tocl精华| 国产精品 欧美亚洲| 国产成+人综合+亚洲专区| 久久欧美精品欧美久久欧美| 国产黄a三级三级三级人| 亚洲一卡2卡3卡4卡5卡精品中文| 在线播放国产精品三级| 欧美性长视频在线观看| 大型黄色视频在线免费观看| 久久香蕉国产精品| 黄片大片在线免费观看| 国产高清videossex| 精品久久久精品久久久| 18禁观看日本| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 十分钟在线观看高清视频www| 大陆偷拍与自拍| avwww免费| 日韩欧美免费精品| 18禁黄网站禁片午夜丰满| 夜夜躁狠狠躁天天躁| 欧美精品亚洲一区二区| 免费少妇av软件| 天天添夜夜摸| 久久这里只有精品19| 亚洲一区二区三区不卡视频| 亚洲精品在线观看二区| 俄罗斯特黄特色一大片| bbb黄色大片| 99国产精品99久久久久| 妹子高潮喷水视频| 亚洲熟妇中文字幕五十中出| 这个男人来自地球电影免费观看| 欧洲精品卡2卡3卡4卡5卡区| 天天一区二区日本电影三级 | 午夜福利视频1000在线观看 | 精品国产一区二区三区四区第35| 日本a在线网址| 国产精品秋霞免费鲁丝片| 丰满人妻熟妇乱又伦精品不卡| 少妇熟女aⅴ在线视频| 久久国产乱子伦精品免费另类| 久久人妻福利社区极品人妻图片| 成年女人毛片免费观看观看9| 欧美乱码精品一区二区三区| 妹子高潮喷水视频| 亚洲精品美女久久久久99蜜臀| 欧美成人一区二区免费高清观看 | 国产精品永久免费网站| 手机成人av网站| 欧美黑人精品巨大| 亚洲五月天丁香| 午夜久久久久精精品| 亚洲国产欧美日韩在线播放| 十八禁网站免费在线| 久久这里只有精品19| 精品久久久精品久久久| 亚洲精品一区av在线观看| 欧美激情高清一区二区三区| e午夜精品久久久久久久| 成人精品一区二区免费| 亚洲自拍偷在线| 免费高清在线观看日韩| 国产极品粉嫩免费观看在线| 精品国内亚洲2022精品成人| 亚洲 国产 在线| videosex国产| 国产区一区二久久| 好看av亚洲va欧美ⅴa在| 人人妻人人澡欧美一区二区 | 夜夜躁狠狠躁天天躁| 国产欧美日韩一区二区三区在线| 十八禁网站免费在线| 国产色视频综合| 国产精品 欧美亚洲| 成年版毛片免费区| 午夜福利一区二区在线看| 一边摸一边抽搐一进一小说| svipshipincom国产片| 变态另类成人亚洲欧美熟女 | 老司机福利观看| 亚洲色图av天堂| 少妇被粗大的猛进出69影院| 欧美在线一区亚洲| 丰满人妻熟妇乱又伦精品不卡| 精品少妇一区二区三区视频日本电影| 色综合亚洲欧美另类图片| 亚洲精品久久国产高清桃花| 淫妇啪啪啪对白视频| 欧美日韩黄片免| 99国产综合亚洲精品| 日韩有码中文字幕| 亚洲国产欧美网| 亚洲成人精品中文字幕电影| 电影成人av| 50天的宝宝边吃奶边哭怎么回事| 999精品在线视频| 免费高清在线观看日韩| 精品午夜福利视频在线观看一区| 国产精品1区2区在线观看.| 欧美午夜高清在线| 好男人在线观看高清免费视频 | 国产激情久久老熟女| 亚洲精品在线美女| 亚洲熟妇熟女久久| 亚洲成av人片免费观看| 精品久久久久久久人妻蜜臀av | 19禁男女啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 成人永久免费在线观看视频| 国产免费av片在线观看野外av| 一区在线观看完整版| 日韩大尺度精品在线看网址 | 欧美乱妇无乱码| 免费一级毛片在线播放高清视频 | 嫁个100分男人电影在线观看| 怎么达到女性高潮| 午夜视频精品福利| 丝袜美腿诱惑在线| 国产激情久久老熟女| 老司机靠b影院| 精品国产乱码久久久久久男人| 国内久久婷婷六月综合欲色啪| 久久香蕉精品热| 日韩视频一区二区在线观看| 亚洲av电影不卡..在线观看| 亚洲欧美激情在线| 国产精品九九99| 日本vs欧美在线观看视频| 欧美性长视频在线观看| 国产精品久久久人人做人人爽| 好男人在线观看高清免费视频 | 50天的宝宝边吃奶边哭怎么回事| 此物有八面人人有两片| 国产高清有码在线观看视频 | 激情视频va一区二区三区| 久久精品人人爽人人爽视色| 久久久久久久午夜电影| 中文字幕人妻丝袜一区二区| 这个男人来自地球电影免费观看| 精品久久蜜臀av无| 一边摸一边抽搐一进一出视频| 欧美亚洲日本最大视频资源| 亚洲狠狠婷婷综合久久图片| 成年女人毛片免费观看观看9| 在线观看免费视频日本深夜| 欧美黄色片欧美黄色片| 可以在线观看的亚洲视频| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| 国产欧美日韩综合在线一区二区| 欧美 亚洲 国产 日韩一| 久久精品国产清高在天天线| 国产av又大| 亚洲国产欧美网| 一区二区日韩欧美中文字幕| 18禁黄网站禁片午夜丰满| 亚洲av五月六月丁香网| 亚洲精品美女久久久久99蜜臀| 成人国语在线视频| 亚洲 欧美 日韩 在线 免费| 国产精品精品国产色婷婷| 色综合欧美亚洲国产小说| 九色亚洲精品在线播放| 亚洲精品美女久久av网站| 校园春色视频在线观看| 国产一级毛片七仙女欲春2 | 国产一区二区三区综合在线观看| 日韩一卡2卡3卡4卡2021年| 精品不卡国产一区二区三区| 精品久久久久久久毛片微露脸| 成人18禁在线播放| 乱人伦中国视频| 嫁个100分男人电影在线观看| 又黄又粗又硬又大视频| 黄片小视频在线播放| 在线av久久热| 亚洲电影在线观看av| 亚洲中文字幕一区二区三区有码在线看 | 亚洲久久久国产精品| 老司机福利观看| 午夜精品国产一区二区电影| 啪啪无遮挡十八禁网站| 亚洲中文日韩欧美视频| 日韩欧美一区二区三区在线观看| 精品国产一区二区久久| 又黄又爽又免费观看的视频| 19禁男女啪啪无遮挡网站| 精品日产1卡2卡| 韩国精品一区二区三区| 大型av网站在线播放| 国产男靠女视频免费网站| 久久人妻熟女aⅴ| 丰满人妻熟妇乱又伦精品不卡| 欧美一区二区精品小视频在线| 日韩av在线大香蕉| 精品国产乱子伦一区二区三区| av欧美777| 国产亚洲av嫩草精品影院| 一级毛片女人18水好多| 男人舔女人下体高潮全视频| 无人区码免费观看不卡| 精品少妇一区二区三区视频日本电影| 午夜福利免费观看在线| 亚洲欧洲精品一区二区精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 一个人观看的视频www高清免费观看 | 色综合站精品国产| 一级毛片女人18水好多| 一区福利在线观看| 伦理电影免费视频| 黄色 视频免费看| 色精品久久人妻99蜜桃| 欧美日本视频| 国产真人三级小视频在线观看| aaaaa片日本免费| 亚洲欧美精品综合一区二区三区| 99国产综合亚洲精品| 亚洲欧美日韩高清在线视频| 免费av毛片视频| 在线av久久热| 老熟妇仑乱视频hdxx| 禁无遮挡网站| 曰老女人黄片| 嫩草影院精品99| 日日爽夜夜爽网站| 国产精品免费一区二区三区在线| 成年女人毛片免费观看观看9| 欧美+亚洲+日韩+国产| 最近最新免费中文字幕在线| 欧美日韩瑟瑟在线播放| 露出奶头的视频| 亚洲色图综合在线观看| 国产91精品成人一区二区三区| 久久这里只有精品19| 国产精品综合久久久久久久免费 | 午夜成年电影在线免费观看| 18禁观看日本| 成人手机av| 亚洲欧美激情在线| 亚洲国产精品成人综合色| 99国产精品一区二区三区| 丝袜美腿诱惑在线| 黄色毛片三级朝国网站| 欧美激情 高清一区二区三区| 亚洲男人的天堂狠狠| 成人三级黄色视频| 国产成人一区二区三区免费视频网站| 麻豆成人av在线观看| 免费看美女性在线毛片视频| 波多野结衣巨乳人妻| 亚洲熟妇熟女久久| 黑人欧美特级aaaaaa片| 午夜福利一区二区在线看| 亚洲国产中文字幕在线视频| 美女午夜性视频免费| 精品久久蜜臀av无| 午夜影院日韩av| 99精品久久久久人妻精品| 久久久久久免费高清国产稀缺| 一个人观看的视频www高清免费观看 | 亚洲美女黄片视频| 好男人在线观看高清免费视频 | 精品少妇一区二区三区视频日本电影| 免费人成视频x8x8入口观看| www日本在线高清视频| 日本a在线网址| 又黄又粗又硬又大视频| 成在线人永久免费视频| 日日摸夜夜添夜夜添小说| 岛国视频午夜一区免费看| 老司机午夜福利在线观看视频| 国产亚洲精品第一综合不卡| 国产亚洲欧美98| 久久久久久人人人人人| 精品熟女少妇八av免费久了| 日本免费一区二区三区高清不卡 | 亚洲精品国产区一区二| 国产精华一区二区三区| 欧美精品亚洲一区二区| 69av精品久久久久久| 欧美性长视频在线观看| 亚洲九九香蕉| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久午夜电影| 久久 成人 亚洲| 亚洲男人天堂网一区| 午夜久久久久精精品| 国产免费av片在线观看野外av| 人妻久久中文字幕网| 亚洲精品久久成人aⅴ小说| 欧美成人免费av一区二区三区| 97人妻精品一区二区三区麻豆 | 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 亚洲九九香蕉| 在线观看舔阴道视频| 亚洲性夜色夜夜综合| 9191精品国产免费久久| 国产精品九九99| 国产成人免费无遮挡视频| 国产精品九九99| 精品国内亚洲2022精品成人| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩高清在线视频| 深夜精品福利| 人人澡人人妻人| 一二三四社区在线视频社区8| 亚洲 欧美一区二区三区| 一二三四社区在线视频社区8| 亚洲在线自拍视频| 日韩欧美一区视频在线观看| 亚洲一区高清亚洲精品| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 国产精品久久电影中文字幕| 午夜两性在线视频| 国产一区二区三区综合在线观看| 国产成人一区二区三区免费视频网站| 超碰成人久久| 国产精品二区激情视频| 精品无人区乱码1区二区| 亚洲成人免费电影在线观看| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| 一夜夜www| 一二三四在线观看免费中文在| 国产精品 国内视频| 免费无遮挡裸体视频| 视频在线观看一区二区三区| 久久香蕉精品热| 午夜免费激情av| 国产精品综合久久久久久久免费 | 中文字幕av电影在线播放| 90打野战视频偷拍视频| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 久热爱精品视频在线9| 亚洲成av人片免费观看| 欧美一区二区精品小视频在线| aaaaa片日本免费| 欧美在线黄色| 高清在线国产一区| 欧美黑人欧美精品刺激| 国产精品电影一区二区三区| 国产成人精品在线电影| 欧美绝顶高潮抽搐喷水| 啪啪无遮挡十八禁网站| 亚洲电影在线观看av| 男人操女人黄网站| 黑人巨大精品欧美一区二区蜜桃| 欧美中文日本在线观看视频| 97人妻精品一区二区三区麻豆 | 久久婷婷人人爽人人干人人爱 | 亚洲精品在线美女|