• <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
    精品国产一区二区久久| 亚洲国产精品成人久久小说| 国产69精品久久久久777片| 国产精品麻豆人妻色哟哟久久| av视频免费观看在线观看| 欧美日本中文国产一区发布| 大香蕉久久成人网| 亚洲综合色网址| 久久 成人 亚洲| 少妇 在线观看| 老女人水多毛片| 人人妻人人爽人人添夜夜欢视频| 久久99一区二区三区| 视频区图区小说| a级毛片免费高清观看在线播放| 毛片一级片免费看久久久久| 超碰97精品在线观看| 日韩电影二区| 久久99一区二区三区| 亚洲av成人精品一二三区| 成人影院久久| 亚洲精品第二区| av在线老鸭窝| 欧美日韩一区二区视频在线观看视频在线| av福利片在线| 成人影院久久| 亚洲精品av麻豆狂野| 色视频在线一区二区三区| 久久久久人妻精品一区果冻| 亚洲色图综合在线观看| 丰满少妇做爰视频| 国产片特级美女逼逼视频| 一本大道久久a久久精品| videossex国产| 亚洲精品av麻豆狂野| 免费大片18禁| 国产精品女同一区二区软件| 天堂8中文在线网| 成人漫画全彩无遮挡| 亚洲欧洲精品一区二区精品久久久 | 欧美性感艳星| 午夜免费鲁丝| videossex国产| 精品卡一卡二卡四卡免费| 91久久精品国产一区二区成人| 国产成人av激情在线播放 | 国产精品99久久99久久久不卡 | 国产av精品麻豆| 亚洲精品日韩在线中文字幕| 午夜福利在线观看免费完整高清在| www.av在线官网国产| 一级毛片aaaaaa免费看小| 国产乱人偷精品视频| 在线观看一区二区三区激情| tube8黄色片| 2021少妇久久久久久久久久久| 久久久久久人妻| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 91精品国产国语对白视频| 丝袜在线中文字幕| 亚洲成色77777| 一区在线观看完整版| 天美传媒精品一区二区| 香蕉精品网在线| 十八禁网站网址无遮挡| 少妇人妻 视频| 99国产精品免费福利视频| 国产av精品麻豆| 99热网站在线观看| 一区二区三区乱码不卡18| 热99国产精品久久久久久7| 91久久精品国产一区二区三区| 久久久久人妻精品一区果冻| 国产精品一区二区在线不卡| 午夜福利影视在线免费观看| 少妇熟女欧美另类| 精品亚洲成国产av| 伊人久久精品亚洲午夜| 伊人亚洲综合成人网| av福利片在线| 亚洲精品日韩av片在线观看| 欧美亚洲日本最大视频资源| 国产成人91sexporn| 永久免费av网站大全| 国产成人精品在线电影| 热99久久久久精品小说推荐| 欧美精品高潮呻吟av久久| 亚洲精品456在线播放app| 成年人午夜在线观看视频| 一本久久精品| 精品少妇内射三级| 日本爱情动作片www.在线观看| 欧美精品国产亚洲| 99热网站在线观看| 日日爽夜夜爽网站| 日本黄色片子视频| 一区二区日韩欧美中文字幕 | 男女国产视频网站| xxxhd国产人妻xxx| 热re99久久国产66热| 美女国产高潮福利片在线看| 国产精品熟女久久久久浪| 一区二区三区免费毛片| 国产免费一区二区三区四区乱码| 99久久综合免费| 亚洲情色 制服丝袜| 有码 亚洲区| 亚洲精品乱久久久久久| 18禁裸乳无遮挡动漫免费视频| 国产精品 国内视频| 午夜福利在线观看免费完整高清在| 插逼视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲精品自拍成人| freevideosex欧美| 又粗又硬又长又爽又黄的视频| 日韩中文字幕视频在线看片| 国产精品一区二区在线不卡| 久久精品国产亚洲网站| 久久婷婷青草| 日本-黄色视频高清免费观看| 十分钟在线观看高清视频www| 黄色毛片三级朝国网站| 成人二区视频| 99re6热这里在线精品视频| 日本黄大片高清| 国国产精品蜜臀av免费| 国产精品久久久久成人av| 亚洲精品国产av蜜桃| 97超碰精品成人国产| 国产精品无大码| √禁漫天堂资源中文www| 久久国内精品自在自线图片| 人人妻人人爽人人添夜夜欢视频| 成人手机av| 精品久久国产蜜桃| 人妻人人澡人人爽人人| 国产午夜精品一二区理论片| 边亲边吃奶的免费视频| 国产精品三级大全| 日韩 亚洲 欧美在线| tube8黄色片| 男女边摸边吃奶| 高清av免费在线| 人妻人人澡人人爽人人| 欧美变态另类bdsm刘玥| 边亲边吃奶的免费视频| 少妇人妻久久综合中文| 极品人妻少妇av视频| 这个男人来自地球电影免费观看 | 久久久精品区二区三区| 日韩免费高清中文字幕av| 久久 成人 亚洲| 午夜91福利影院| 日日摸夜夜添夜夜爱| 人人澡人人妻人| 五月玫瑰六月丁香| 一区二区av电影网| 日本vs欧美在线观看视频| 亚洲av福利一区| 日韩不卡一区二区三区视频在线| 欧美日韩一区二区视频在线观看视频在线| 国产白丝娇喘喷水9色精品| 九草在线视频观看| 欧美激情极品国产一区二区三区 | 亚洲av不卡在线观看| 丝袜在线中文字幕| 亚洲av不卡在线观看| 蜜桃久久精品国产亚洲av| 丁香六月天网| 亚洲第一av免费看| 国产乱人偷精品视频| 国产一区有黄有色的免费视频| 成年av动漫网址| 国产日韩欧美视频二区| 久久热精品热| 精品国产一区二区三区久久久樱花| 97超碰精品成人国产| 精品酒店卫生间| 夜夜看夜夜爽夜夜摸| freevideosex欧美| 精品国产一区二区久久| 又粗又硬又长又爽又黄的视频| √禁漫天堂资源中文www| 中文字幕最新亚洲高清| 一二三四中文在线观看免费高清| 丰满少妇做爰视频| a级毛片黄视频| 亚洲四区av| 色哟哟·www| 中文字幕人妻丝袜制服| 人妻一区二区av| 久久精品熟女亚洲av麻豆精品| 人妻 亚洲 视频| 免费看不卡的av| 老司机影院毛片| 狠狠婷婷综合久久久久久88av| 国产免费一级a男人的天堂| 在线看a的网站| 22中文网久久字幕| 亚洲性久久影院| 2021少妇久久久久久久久久久| 曰老女人黄片| 性色av一级| 国产黄色视频一区二区在线观看| 五月天丁香电影| 欧美日韩国产mv在线观看视频| 夜夜骑夜夜射夜夜干| 亚洲精品第二区| 国产在线一区二区三区精| 亚洲精品中文字幕在线视频| 日韩亚洲欧美综合| 久久久久网色| 99热这里只有是精品在线观看| 99久久人妻综合| 在线观看人妻少妇| 欧美变态另类bdsm刘玥| 国产精品一区二区在线观看99| 美女大奶头黄色视频| 国产精品偷伦视频观看了| 全区人妻精品视频| 亚洲av欧美aⅴ国产| 久久久久久久久久久免费av| 国产一区二区三区综合在线观看 | 亚洲一区二区三区欧美精品| 欧美激情国产日韩精品一区| 欧美激情 高清一区二区三区| av在线app专区| 赤兔流量卡办理| 久久这里有精品视频免费| 免费人成在线观看视频色| 国产在视频线精品| 免费日韩欧美在线观看| 亚洲国产精品一区三区| 国产男女内射视频| 不卡视频在线观看欧美| 国产在线视频一区二区| 夜夜看夜夜爽夜夜摸| 久久久精品免费免费高清| 国产精品麻豆人妻色哟哟久久| 久久狼人影院| 国产不卡av网站在线观看| 十八禁网站网址无遮挡| 欧美激情国产日韩精品一区| 一个人看视频在线观看www免费| 乱人伦中国视频| 97精品久久久久久久久久精品| 王馨瑶露胸无遮挡在线观看| 日韩欧美精品免费久久| 自线自在国产av| 午夜激情av网站| 亚洲人与动物交配视频| 午夜福利视频精品| 视频区图区小说| 美女主播在线视频| 免费人妻精品一区二区三区视频| 精品国产露脸久久av麻豆| 人妻一区二区av| 欧美激情国产日韩精品一区| 在线精品无人区一区二区三| 黄色配什么色好看| 久久精品国产亚洲av涩爱| 欧美激情 高清一区二区三区| 男人添女人高潮全过程视频| 成年女人在线观看亚洲视频| 久久精品熟女亚洲av麻豆精品| 91成人精品电影| 精品亚洲成a人片在线观看| 亚洲丝袜综合中文字幕| 18禁观看日本| 精品一区二区三区视频在线| av在线播放精品| 美女xxoo啪啪120秒动态图| 2018国产大陆天天弄谢| 成人国产麻豆网| 久久久精品94久久精品| 菩萨蛮人人尽说江南好唐韦庄| 尾随美女入室| .国产精品久久| 两个人的视频大全免费| 五月伊人婷婷丁香| 国产精品三级大全| 国产综合精华液| 久久午夜福利片| 国产乱来视频区| 久久久国产欧美日韩av| 国产精品不卡视频一区二区| 人妻一区二区av| 成人免费观看视频高清| 十八禁高潮呻吟视频| 女的被弄到高潮叫床怎么办| 亚洲成人一二三区av| 在线亚洲精品国产二区图片欧美 | 在线观看人妻少妇| 欧美激情国产日韩精品一区| 国产精品.久久久| 亚洲av中文av极速乱| 啦啦啦在线观看免费高清www| 嫩草影院入口| h视频一区二区三区| 最新中文字幕久久久久| 老司机亚洲免费影院| 精品人妻熟女毛片av久久网站| 国产精品秋霞免费鲁丝片| 国产精品99久久99久久久不卡 | 一本一本综合久久| 晚上一个人看的免费电影| 午夜影院在线不卡| 日韩精品免费视频一区二区三区 | 久久99热这里只频精品6学生| 日本av手机在线免费观看| 精品久久国产蜜桃| 成人漫画全彩无遮挡| 观看美女的网站| 一本久久精品| 91精品伊人久久大香线蕉| 一本一本综合久久| 少妇 在线观看| 如何舔出高潮| 欧美97在线视频| 久久精品国产鲁丝片午夜精品| 久久99热6这里只有精品| 国产精品嫩草影院av在线观看| 王馨瑶露胸无遮挡在线观看| 国产黄片视频在线免费观看| 日韩亚洲欧美综合| 国产乱来视频区| 午夜激情福利司机影院| 亚洲四区av| 国产精品久久久久久久久免| 国产一区有黄有色的免费视频| 国产探花极品一区二区| 最近中文字幕高清免费大全6| 国产精品免费大片| 亚洲,一卡二卡三卡| 欧美 亚洲 国产 日韩一| 九色亚洲精品在线播放| 又黄又爽又刺激的免费视频.| 亚洲经典国产精华液单| 亚洲av中文av极速乱| 亚洲av.av天堂| 国产精品三级大全| a级毛片黄视频| 亚洲国产精品成人久久小说| 国产老妇伦熟女老妇高清| 纯流量卡能插随身wifi吗| 又黄又爽又刺激的免费视频.| 在线观看免费视频网站a站| 亚洲经典国产精华液单| 91久久精品电影网| 亚洲国产欧美在线一区| 久久人人爽人人爽人人片va| 丰满乱子伦码专区| 黄色毛片三级朝国网站| 我要看黄色一级片免费的| 黄色视频在线播放观看不卡| 久久精品国产亚洲av天美| 亚洲av免费高清在线观看| 特大巨黑吊av在线直播| 日韩一区二区视频免费看| 亚洲丝袜综合中文字幕| 国产在线一区二区三区精| 亚洲精华国产精华液的使用体验| 九九在线视频观看精品| 国产探花极品一区二区| 国产精品国产av在线观看| 你懂的网址亚洲精品在线观看| 91午夜精品亚洲一区二区三区| 亚洲怡红院男人天堂| 亚洲国产欧美在线一区| 中文字幕最新亚洲高清| 熟女人妻精品中文字幕| 少妇 在线观看| 国产淫语在线视频| 国产精品国产av在线观看| 永久免费av网站大全| 777米奇影视久久| 极品人妻少妇av视频| 国产伦理片在线播放av一区| 欧美精品人与动牲交sv欧美| 99热国产这里只有精品6| 欧美bdsm另类| 国产精品人妻久久久影院| 午夜av观看不卡| 97超视频在线观看视频| 有码 亚洲区| 在现免费观看毛片| 亚洲精品456在线播放app| 涩涩av久久男人的天堂| 欧美老熟妇乱子伦牲交| 欧美日韩一区二区视频在线观看视频在线| 又大又黄又爽视频免费| 99九九线精品视频在线观看视频| 国产老妇伦熟女老妇高清| 天堂俺去俺来也www色官网| 精品久久久噜噜| 寂寞人妻少妇视频99o| 制服诱惑二区| 日本猛色少妇xxxxx猛交久久| 亚洲精品久久午夜乱码| 18禁在线无遮挡免费观看视频| 一级二级三级毛片免费看| 亚洲av不卡在线观看| 激情五月婷婷亚洲| 3wmmmm亚洲av在线观看| 亚洲伊人久久精品综合| 全区人妻精品视频| 亚洲av男天堂| av天堂久久9| 色视频在线一区二区三区| www.色视频.com| 亚洲精品日韩av片在线观看| 99精国产麻豆久久婷婷| 一个人免费看片子| 午夜福利在线观看免费完整高清在| 九草在线视频观看| 最黄视频免费看| 国产精品99久久久久久久久| 免费看av在线观看网站| 自线自在国产av| 亚洲在久久综合| 国产一区二区三区av在线| 久久ye,这里只有精品| 久久国内精品自在自线图片| 亚洲国产毛片av蜜桃av| a级片在线免费高清观看视频| 国产av码专区亚洲av| 搡老乐熟女国产| 久久99热6这里只有精品| 亚洲成色77777| 国产欧美另类精品又又久久亚洲欧美| av不卡在线播放| 午夜福利,免费看| √禁漫天堂资源中文www| 久久国产精品男人的天堂亚洲 | 丝瓜视频免费看黄片| 欧美亚洲 丝袜 人妻 在线| 五月天丁香电影| 色5月婷婷丁香| 精品久久国产蜜桃| 自线自在国产av| 18在线观看网站| 麻豆成人av视频| 777米奇影视久久| 久久久久久久久久久丰满| 大香蕉久久成人网| 久久久久精品性色| 国产精品一区二区在线观看99| 欧美三级亚洲精品| 超色免费av| 天天操日日干夜夜撸| 国产亚洲午夜精品一区二区久久| 国产黄色视频一区二区在线观看| 色5月婷婷丁香| 男女啪啪激烈高潮av片| 色94色欧美一区二区| 99热国产这里只有精品6| 大片免费播放器 马上看| 亚洲精品成人av观看孕妇| 少妇的逼水好多| 91国产中文字幕| av在线播放精品| 亚洲国产精品国产精品| 国产熟女欧美一区二区| 草草在线视频免费看| 久久久国产欧美日韩av| 制服丝袜香蕉在线| 亚洲av.av天堂| 秋霞在线观看毛片| 亚洲成色77777| 18在线观看网站| 91在线精品国自产拍蜜月| 久久国内精品自在自线图片| 熟妇人妻不卡中文字幕| 国产精品久久久久久久电影| 九色成人免费人妻av| 欧美 亚洲 国产 日韩一| 涩涩av久久男人的天堂| 精品人妻熟女av久视频| 亚洲欧美一区二区三区国产| 黄色欧美视频在线观看| 日韩熟女老妇一区二区性免费视频| 成人手机av| 七月丁香在线播放| 亚洲综合精品二区| 老女人水多毛片| 免费高清在线观看日韩| 国产成人91sexporn| 插阴视频在线观看视频| 久久精品国产a三级三级三级| 精品人妻熟女毛片av久久网站| 只有这里有精品99| 国产伦精品一区二区三区视频9| 男男h啪啪无遮挡| 国产欧美日韩一区二区三区在线 | 最后的刺客免费高清国语| av线在线观看网站| 狠狠婷婷综合久久久久久88av| 日韩大片免费观看网站| 国产精品久久久久久av不卡| 亚洲国产成人一精品久久久| 久久久久久久久久人人人人人人| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| 国产伦理片在线播放av一区| 亚洲美女黄色视频免费看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av成人精品| 日本欧美视频一区| 日韩av在线免费看完整版不卡| 哪个播放器可以免费观看大片| 一区二区三区四区激情视频| 国产精品久久久久久精品古装| 国产色婷婷99| 边亲边吃奶的免费视频| 桃花免费在线播放| 下体分泌物呈黄色| 成人午夜精彩视频在线观看| 久久久亚洲精品成人影院| 成人午夜精彩视频在线观看| 我的老师免费观看完整版| 国产男女内射视频| 2021少妇久久久久久久久久久| 秋霞伦理黄片| 亚洲av二区三区四区| 啦啦啦在线观看免费高清www| 97精品久久久久久久久久精品| 少妇被粗大的猛进出69影院 | 999精品在线视频| 久久久久久久大尺度免费视频| 午夜福利,免费看| 久久久亚洲精品成人影院| 久久影院123| 天天影视国产精品| 桃花免费在线播放| av在线app专区| 日韩大片免费观看网站| 丁香六月天网| 日本av免费视频播放| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 亚洲人与动物交配视频| 熟女电影av网| 午夜福利影视在线免费观看| 欧美日韩一区二区视频在线观看视频在线| videossex国产| 全区人妻精品视频| 国产免费一区二区三区四区乱码| 久久久久久久久久久免费av| 久久午夜综合久久蜜桃| 成人午夜精彩视频在线观看| 国产成人精品婷婷| 久久久a久久爽久久v久久| 啦啦啦在线观看免费高清www| 久热这里只有精品99| av天堂久久9| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| a级毛片免费高清观看在线播放| 黄色欧美视频在线观看| 又粗又硬又长又爽又黄的视频| 国产日韩欧美在线精品| 亚洲av福利一区| 亚洲第一av免费看| 午夜精品国产一区二区电影| 视频区图区小说| 丁香六月天网| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 国产成人精品婷婷| 久久婷婷青草| 日本猛色少妇xxxxx猛交久久| 午夜影院在线不卡| 亚洲精品色激情综合| 国产免费又黄又爽又色| 亚洲精品色激情综合| 日本欧美视频一区| 久久久久视频综合| 另类精品久久| 波野结衣二区三区在线| 免费av中文字幕在线| 丝瓜视频免费看黄片| 国产精品国产三级专区第一集| 嫩草影院入口| 亚洲性久久影院| 久久久久久久久大av| 在线观看www视频免费| 老司机影院毛片| 校园人妻丝袜中文字幕| 99久久精品一区二区三区| 制服丝袜香蕉在线| 久久精品久久精品一区二区三区| 日韩精品免费视频一区二区三区 | 国精品久久久久久国模美| av视频免费观看在线观看| 黑人高潮一二区| 精品99又大又爽又粗少妇毛片| 视频在线观看一区二区三区| 日日爽夜夜爽网站| 国产不卡av网站在线观看| 亚洲精品乱久久久久久| 国产精品蜜桃在线观看| 2022亚洲国产成人精品| 老司机影院成人| 精品人妻熟女毛片av久久网站| 色吧在线观看| 国产精品国产三级国产专区5o| 永久网站在线| 卡戴珊不雅视频在线播放| kizo精华| 精品久久国产蜜桃| 97在线人人人人妻| 国产黄频视频在线观看|