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

    基于地基雷達探究和田市一次沙塵污染過程

    2024-01-06 01:33:50陰璐璐李京龍付光祥吳彩云肖鴻丹范旭雨晨
    中國環(huán)境科學 2023年12期
    關(guān)鍵詞:散射系數(shù)沙塵沙塵暴

    陰璐璐,何 清*,李京龍,孟 露,付光祥,吳彩云,肖鴻丹,范旭雨晨

    基于地基雷達探究和田市一次沙塵污染過程

    陰璐璐1,2,何 清1,2*,李京龍2,3,孟 露2,3,付光祥1,2,吳彩云2,3,肖鴻丹2,范旭雨晨2

    (1.新疆師范大學地理科學與旅游學院,新疆 烏魯木齊 830054;2.中國氣象局烏魯木齊沙漠氣象研究所/新疆塔克拉瑪干沙漠氣象國家野外科學觀測研究站/中國氣象局塔克拉瑪干沙漠氣象野外科學試驗基地/新疆維吾爾自治區(qū)沙漠氣象與沙塵暴重點實驗室,新疆 烏魯木齊 830002;3.新疆大學地理與遙感科學學院,新疆 烏魯木齊 830017)

    2022年5月新疆和田市遭遇了近5a來最強的一次沙塵暴天氣侵襲,沙塵污染過程持續(xù)了一周.利用大氣氣溶膠激光雷達、環(huán)境監(jiān)測站污染物濃度和近地面氣象觀測資料,分析了此次沙塵污染過程中顆粒物濃度變化、氣溶膠光學信息垂直分布特征,并使用全球再分析資料得出本次污染天氣的氣象匹配形態(tài),最后結(jié)合HYSPLIT后向軌跡模型解析過境期間氣團的輸送及潛在來源.結(jié)果表明:沙塵污染期間,PM2.5和PM10變化趨勢較為一致且在垂直高度上存在非均一性,PM2.5/PM10均值為0.25,表明主要以粗顆粒物為主,風速與溫度的增加、濕度的降低有利于促成污染天氣的發(fā)生;沙塵污染期間出現(xiàn)了明顯的氣溶膠層,在垂直方向上有顯著的波動性變化,沙塵發(fā)生當日消光系數(shù)、后向散射系數(shù)和退偏振比達到了本次過程最大值,分別為3.5km-1、0.07km-1·sr-1和0.17;和田市在春季處于熱低壓場中,在西風背景環(huán)境與塔里木盆地周邊地形阻擋下污染物難以稀釋與擴散,沙塵氣團長時間在此盤踞,是本次污染持續(xù)時間較長的原因;本次沙塵天氣過程污染物主要來自于新疆北部、新疆南部的塔克拉瑪干沙漠以及周邊的斯坦國家且污染源主要以內(nèi)源性沙塵氣溶膠的排放貢獻為主.

    和田市;沙塵氣溶膠;垂直分布;光學信息

    大氣氣溶膠影響著大氣輻射傳輸過程、云的形成、人體健康和出行安全[1].沙塵作為氣溶膠的重要來源之一,在大氣氣溶膠中有著舉足輕重的地位,它作為大氣氣溶膠的主要組成部分[2-4],不但可以影響天氣和氣候[5-7],而且對海洋生態(tài)及全球物質(zhì)循環(huán)也具有重要的作用[8-11],是當前地球科學及環(huán)境變化領(lǐng)域研究的熱點問題[12-14].我國沙塵天氣高發(fā)區(qū)主要位于民豐與和田地區(qū)的南疆盆地和民勤到吉蘭泰為中心的河西地帶,其中沙塵暴與浮塵天氣的頻發(fā)區(qū)位于我國的南疆地區(qū).目前,沙塵暴發(fā)生的頻率與發(fā)展的趨勢沒有一個全球模式,在多數(shù)情況下,自然過程決定了任何一年中它的發(fā)生頻率,并且人類活動也決定著沙塵暴頻率的變異性[15].在全球干旱半干旱地區(qū)的相關(guān)研究領(lǐng)域中,沙塵暴作為一種特殊且影響顯著的災害性天氣,至今仍然是國內(nèi)外學者所探究的重點[16-18].同時,氣象背景條件是影響沙塵天氣發(fā)展變化的關(guān)鍵因素,有學者研究發(fā)現(xiàn)氣象條件直接作用大氣顆粒物的積累、傳輸和擴散,其中風作為起沙最重要的氣象要素,當其超過10m/s時就會產(chǎn)生傳輸沉淀物運動,使大氣中顆粒物濃度迅速上升[19],而污染天氣的出現(xiàn)也往往與弱氣壓場有關(guān)[20-21].塔克拉瑪干沙漠作為我國及東亞地區(qū)重要的沙塵源區(qū),對區(qū)域甚至全球的能量收支、生態(tài)系統(tǒng)及水循環(huán)過程都有非常重要的影響.和田市位于塔克拉瑪干沙漠南緣,其靠近沙漠獨特的地理位置與大氣環(huán)流特征決定了南疆的和田市成為北方沙塵天氣高發(fā)區(qū)之一[22-25].

    大氣氣溶膠垂直結(jié)構(gòu)的觀測手段一般包括了現(xiàn)場直接觀測和基于遙感技術(shù)的觀測.現(xiàn)場直接觀測是指利用無人機、飛機、塔臺以及探空氣球等設(shè)備主動地獲取實驗數(shù)據(jù),基于遙感技術(shù)的觀測主要包括了地基激光雷達、微波輻射計以及星載激光雷達等.地基雷達是當今科研界用來觀測大氣氣溶膠的有效手段之一,它具有分辨率高、體積小、可以連續(xù)且長時間觀測的優(yōu)勢,被廣泛地用于探測AOD(氣溶膠光學厚度)、氣溶膠的垂直結(jié)構(gòu)以及大氣邊界層高度等領(lǐng)域中.目前,國內(nèi)外相關(guān)學者圍繞塔克拉瑪干沙漠沙塵氣溶膠的垂直分布[26-27]、時空分布特征[28]以及質(zhì)量濃度[29]的研究較多且成果豐碩,但是目前以沙漠為背景且基于地基雷達對沙塵污染過程中顆粒物垂直變化”點”的個例分析研究還相對較少.

    根據(jù)地面氣象觀測數(shù)據(jù),顯示2022年5月中旬,新疆維吾爾自治區(qū)南疆盆地與東疆的11個縣市遭遇了近5a以來最強沙塵暴天氣,據(jù)中國氣象局相關(guān)數(shù)據(jù)統(tǒng)計,本次沙塵天氣過程主要集中在喀什地區(qū)、和田地區(qū)以及哈密市,過程的綜合強度可達6.76,其中和田市能見度明顯低于500m,空氣污染指數(shù)嚴重超標,給當?shù)厝嗣竦纳a(chǎn)與生活造成了嚴重的影響.和田市空氣質(zhì)量報告顯示2022年5月12~18日和田市均處于污染狀態(tài),這是近些年以來少見的一次典型沙塵污染天氣過程.鑒于此,本研究基于環(huán)境監(jiān)測站污染物濃度數(shù)據(jù),分析了和田市PM10、PM2.5濃度與近地面氣象要素之間的關(guān)系,并利用大氣氣溶膠激光雷達數(shù)據(jù)探究沙塵天氣過程沙塵氣溶膠垂直變化以及光學信息特征,最后結(jié)合全球再分析資料(ERA5)和HYSPLIT(混合單粒子拉格朗日綜合軌跡模型)后向軌跡模型追蹤氣候變化背景下沙塵可能傳輸?shù)穆窂?并對沙塵源地進行判斷,以期為大氣污染防治提供有效的參考.

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

    1.1 產(chǎn)品簡介

    表1 產(chǎn)品技術(shù)參數(shù)指標

    本文所使用的設(shè)備為光博量子有限公司生產(chǎn)的大氣氣溶膠激光雷達(GBQL-01型),它是由激光器、控制單元板、光學單元、高速信號采集卡、工業(yè)主板和通信模塊等部分構(gòu)成.通過探測激光和大氣中各種分子與氣溶膠相互作用的后向散射信號來反演大氣特性.GBQL-01型氣溶膠激光雷達采用了高能量、高重頻的532nm波長的線偏振激光,克服了傳統(tǒng)米散射激光雷達白天信噪比低的問題,做到了晝夜高信噪比全天時工作.另外,它依賴其較高的脈沖頻率能實現(xiàn)較高的時間分辨率,并通過搭配3D掃描云臺實現(xiàn)大范圍空間掃描探測,具體參數(shù)如表1所示.

    1.2 實驗地點與實驗設(shè)計

    和田市(圖1)位于塔克拉瑪干沙漠南緣,南部毗鄰青藏高原,是南疆地區(qū)政治、經(jīng)濟、文化重要的交流陣地.據(jù)統(tǒng)計,在極端氣候條件下其年平均降塵量約為1400t/km2,年沙塵暴天數(shù)為32.9d,揚沙天數(shù)為53.ld,浮塵天數(shù)甚至可達202.4d,和田年內(nèi)沙塵天氣主要集中在春季和夏季,占全年沙塵天氣總數(shù)的70.23%,其中5月份浮塵和揚沙天氣發(fā)生頻次最高,是我國北方地區(qū)遭受沙塵天氣危害最嚴重的城市之一[30],作為沙塵暴天氣的頻發(fā)地區(qū),探究該地沙塵氣溶膠的垂直分布與光學信息垂直結(jié)構(gòu)特征有著重要意義.雷達觀測場位于新疆和田市氣象局辦公樓樓頂(37.12°N,79.92°E),地面平整、視野開闊,周圍無任何高大的障礙物,滿足實驗觀測的條件.

    圖1 研究區(qū)概況

    1.3 地基雷達的反演與原理

    本研究利用了地基雷達信號反演的要素包括了消光系數(shù)、退偏振比和后向散射系數(shù).在進行反演前需要確定回波信號具有較高的信噪比且信號可靠,隨后對數(shù)據(jù)進行噪聲訂正、重疊因子訂正等預處理,基于米散射-偏振技術(shù),激光雷達信號的反演采用了當今主流的Fernald算法,激光雷達反演方程式如下所示:

    式中:()為雷達接受到的大氣后向散射回波功率;為雷達系統(tǒng)常數(shù);為幾何因子;2為探測距離;m為大氣分子的后向散射系數(shù);m為大氣分子的消光系數(shù);a為氣溶膠消光系數(shù).

    氣溶膠消光系數(shù)反映了不同的氣溶膠粒子在大氣中散射和吸收作用,是對某一波段太陽輻射衰減的綜合描述,主要貢獻來自各種形式的降水、顆粒物和氣態(tài)污染物,而顆粒物的消光系數(shù)可以占到總消光系數(shù)的90%,通常比值越大說明其能見度較低,污染也就越嚴重;后向散射系數(shù)是用來表征氣溶膠粒子后向散射強度的物理量,其數(shù)值越大則代表粒子的散射能力強;而退偏振比值則反映的是大氣氣溶膠和云粒子的非球形特征,波動范圍為0~1,即比值越接近0則代表其粒子形狀越規(guī)則;越接近1表示非球形程度很強,而沙塵多數(shù)為不規(guī)則粒子,粗顆粒物占比較高,因此退偏振比值也較高.相關(guān)研究認為沙塵氣溶膠退偏振比一般大于0.1,最大則可以達到0.4[31-33].因此,一般結(jié)合氣溶膠激光雷達的消光系數(shù)、退偏振比和后向散射系數(shù)等光學參數(shù)來判斷分析沙塵.

    1.4 數(shù)據(jù)與模式

    其他數(shù)據(jù)包括大氣顆粒物數(shù)據(jù),來自中國環(huán)境監(jiān)測總站的全國城市空氣質(zhì)量實時發(fā)布平臺;近地面氣象資料來自于和田市氣象局常規(guī)氣象觀測,包括風速、溫度、濕度、氣壓、能見度、極大風速、最大風速數(shù)據(jù);地基雷達數(shù)據(jù)來自于GBQL-01型大氣氣溶膠激光雷達的實時觀測;不同高度溫度、氣壓、風場與位勢高度等數(shù)據(jù)來自于ERA5,即ECMWF(歐洲中期天氣預報中心)全球氣候的第五代大氣再分析數(shù)據(jù)集;混合單粒子拉格朗日綜合軌跡模式(HYSPLIT)所用氣象資料來自于美國氣象環(huán)境預報中心(NCEP)提供的全球資料同化系統(tǒng)(GDAS)氣象再分析資料.

    HYSPILT模型是由澳大利亞氣象局(BOM)和美國國家海洋和大氣管理局(NOAA)所共同研發(fā),用以計算和分析大氣污染物傳輸與擴散軌跡[34-35].由于該模型能夠較好地模擬和預測出沙塵或其他類型污染物的輸送、擴散以及沉降過程,對多種的氣象要素輸入場、物理過程的污染物排放源進行很好地處理,現(xiàn)如今已被廣泛地應(yīng)用于大氣污染物的傳輸過程分析中[36-37].本文所使用的是HYSPLIT-4后向軌跡模式,采用了NCEP提供的全球資料同化系統(tǒng)數(shù)據(jù)(GDAS,1°×1°)提供的三維氣象要素場作為模式的輸入場進行同步模擬.本文以和田市為擬受點,時間為2022年5月12日22:00和2月18日22:00進行氣團的72h后向軌跡模擬,時間分辨率為6h,進而判斷本次沙塵污染天氣中氣溶膠的來源與傳輸過程.

    依據(jù)《地面氣象觀測規(guī)范》[38]中的《沙塵暴天氣等級》,結(jié)合沙塵天氣地面能見度水平來劃分天氣狀況,分別為浮塵、揚沙、沙塵暴、強沙塵暴和特強沙塵暴.當處于無風或平均風速小于3m/s的天氣條件下,能見度低于10km被定義為浮塵;本地或近地面沙塵被風吹起使得空氣混濁,能見度在1~10km被定義為揚沙;大風將地面沙塵吹起,能見度低于1km定義為沙塵暴;最大風速大于或等于20m/s,能見度低于500m被定義為強沙塵暴;當最大風速大于或等于25m/s且能見度小于50m時則被定義為特強沙塵暴.

    2 結(jié)果與討論

    2.1 污染過程顆粒物濃度變化及地面氣象要素分析

    圖2為本次觀測期間顆粒物濃度與近地面氣象要素的時間變化.和田市在本次污染期間PM10與PM2.5變化趨勢較為一致且在垂直高度上存在很強的非均一性.污染期間顆粒物質(zhì)量濃度均高于GB 3095-2012《環(huán)境空氣質(zhì)量標準》[39]日平均二級標準限值(PM10和PM2.5日平均濃度閾值分別為150和75μg/m3).5月12日2:00~6:00,風速從2.4m/s迅速增加到10.5m/s,顆粒物濃度也呈現(xiàn)出強烈的上升趨勢,PM2.5從73μg/m3增至118μg/m3,PM10從205μg/m3增至661μg/m3,PM2.5/PM10的比值由0.3降至0.18,能見度明顯下降,從10039m降至1520m,開始由晴天變?yōu)檫B續(xù)的揚沙天氣;16:00至次日14:00,沙塵暴天氣發(fā)生,期間風速與濕度迅速增加,分別為15.6m/s和54%,而溫度持續(xù)降低,PM2.5與PM10陡升并出現(xiàn)了此次污染過程的最大值,分別為1174和4708μg/m3,且PM2.5/PM10的比值為0.24,呈現(xiàn)波動遞減趨勢,能見度達到此次沙塵污染天氣的最低值,僅為373m;伴隨著風速的降低,沙塵暴天氣結(jié)束,又開始了長達4d的揚沙天氣,能見度持續(xù)在10km內(nèi),PM10和PM2.5濃度較之前有所回落,但空氣污染程度依舊嚴重超標;5月17日7:00至次日12:00,風速明顯降低,最低達0.9m/s,夜間溫度較低但此時伴有降水天氣濕度較高,最高值為63%,天氣開始由揚沙轉(zhuǎn)變?yōu)楦m;隨后天氣變化較為劇烈,揚沙與浮塵天氣相間出現(xiàn),直至5月18日1:00,能見度高達12687m,天氣狀況變好,一次持續(xù)的沙塵污染天氣結(jié)束.污染期間PM2.5/PM10的比值整體處于0~0.49之間,均值為0.25,表明粗顆粒物在本次污染期間占據(jù)主導地位.由此看來,PM10和PM2.5濃度的變化與沙塵天氣過程發(fā)展具有較好的時間一致性,風速與濕度的增大、溫度的降低有利于促成沙塵污染天氣的發(fā)生,這與先前所研究的結(jié)果一致[40-41].而在之后的天氣發(fā)展過程中,雖然有風速高值出現(xiàn),但顆粒物的污染濃度卻呈現(xiàn)下降趨勢,這主要是因為風速增加會促進污染物的消散[42].

    2.2 氣溶膠光學信息垂直分布特征

    圖3 沙塵過境期間大氣氣溶膠消光系數(shù)、退偏振比和后向散射系數(shù)變化

    地基激光雷達對于地面以上以及云層以下氣溶膠垂直探測具有明顯的優(yōu)勢,能夠很好地分析沙塵與天氣污染過程[43-44].圖3為地基激光雷達反演得到的2022年5月12~18日沙塵過境前后氣溶膠消光系數(shù)、退偏振比和后向散射系數(shù)的垂直分布變化情況.由圖3可知,沙塵污染期間出現(xiàn)了明顯的氣溶膠層,并且在垂直方向上有顯著的波動性變化.5月12日為沙塵事件始發(fā)當日,局部高空和近地面PM10與PM2.5質(zhì)量的濃度均有所上升,使得消光系數(shù)、退偏振比與后向散射系數(shù)的變化都非常明顯,在強大上升氣流與風的剪切作用下卷起地面沙塵從而引起了氣溶膠的觀測值發(fā)生變化.自5月12日起,近地面消光系數(shù)明顯升高并持續(xù)至5月16日,此期間消光系數(shù)越大反映出沙塵氣溶膠含量越高,污染開始持續(xù)增強,這與圖2中顆粒物濃度在5月12~16日呈現(xiàn)高值所得出的結(jié)果一致.在5:00之前大氣邊界層比較稀薄,但在2.5~6km的高空出現(xiàn)了明顯的氣溶膠層,消光系數(shù)在0.5~0.9km-1之間,且退偏振比大于0.1,判斷該時間段的高空為形狀不規(guī)則的沙塵氣溶膠因風速迅速增大被吹揚至高空,5:00開始近地面消光系數(shù)高達1km-1左右,后向散射系數(shù)和退偏振比同時也呈現(xiàn)增大趨勢.在3km左右的高空至地面,由于氣流的下沉效應(yīng)并伴隨著風速的降低,邊界層上大量的高空沙塵開始逐漸向下混合并延展至地表沉降.結(jié)合圖2可知此時大氣中粗顆粒物所占比例明顯大于細顆粒物,和田市在此期間近地層污染十分嚴重.整個沙塵影響期間,除個別區(qū)域高度為監(jiān)測盲區(qū)導致數(shù)據(jù)缺測之外,退偏振比的整體在0.13左右,這表明在6km高度以下沙粒子分布比較均勻.在5月16日10:00前后,消光系數(shù)、退偏振比和后向散射系數(shù)隨著高度增加開始降低,沙塵逐漸消散,濃度開始有所減弱,直至18日1:00沙塵污染過程結(jié)束,消光系數(shù)、后向散射系數(shù)和退偏振比也呈現(xiàn)出明顯降低的趨勢.

    圖4為2022年5月12日02:00、5月12日15:00、5月15日02:00和5月18日03:00四個典型時刻的消光系數(shù)、后向散射系數(shù)和退偏振比的廓線圖,分別代表著本次沙塵污染過程前、污染中、污染積累和污染消散時段內(nèi)氣溶膠光學信息垂直分布情況.從圖上可以直觀地看出在5月12日02:00,沙塵氣溶膠主要分布在5km以下的高空范圍內(nèi),此時為沙塵污染的前期,污染程度較低,消光系數(shù)在1km-1左右,后向散射系數(shù)最高到0.02km-1·sr-1,退偏振比大于0.1;12日15:00為沙塵天氣污染嚴重時刻,沙塵氣溶膠主要集中分布在1.3km高空以下,消光系數(shù)、后向散射系數(shù)與退偏振比均呈現(xiàn)出大幅上漲趨勢并達到本次沙塵污染天氣過程的最大值,此時分別為3.5km-1、0.07km-1·sr-1和0.17,5月15日02:00至18日03:00,沙塵天氣污染過程有所緩和,消光系數(shù)與后向散射系數(shù)也明顯降低,沙塵氣溶膠所處的高度也在逐漸下降.

    圖4 激光雷達反演得到污染期間內(nèi)四個典型時刻消光系數(shù)、后向散射系數(shù)和退偏振比廓線圖

    2.3 沙塵污染天氣氣象背景分析

    強大的風力、充足的沙源以及不穩(wěn)定的大氣層結(jié)是導致沙塵天氣發(fā)生的原因[45-47],而蒙古氣旋與冷鋒是造成本次和田市沙塵天氣發(fā)生的主要影響系統(tǒng).此外,天氣形勢發(fā)展決定著沙塵顆粒物濃度的變化,這對深入剖析沙塵污染天氣過程起著重要的指示作用.為研究本次沙塵過境期間影響沙塵起沙的天氣狀況,分別對和田市2022年5月12日08:00(污染前)、12日20:00 (污染中)、15日(污染累積)以及18日(污染消散)的天氣形勢做出了分析,如圖5~7所示.

    結(jié)合500hPa位勢高度和溫度場圖(圖5)與500hPa風場圖(圖6)分析可知,新疆位于北半球,因此受西風環(huán)流影響較大[48],并且新疆北部地區(qū)夏季也主要受到蒙古西伯利亞高壓的控制.沙塵天氣發(fā)生前(圖5a),在里海以北、貝加爾湖以西存在豎槽,高層為偏西風且和田地區(qū)風力較為微弱,使得沙塵發(fā)生前日浮塵天氣持續(xù)時間較長,和田市屬于槽前,溫度槽落后于高度槽并配合有冷平流,這在未來有利于高度槽發(fā)展;在沙塵發(fā)生當日圖5(b),可見低壓槽位于蒙古國地區(qū)及中國東北地區(qū)附近并不斷地向東移動,而盆地東側(cè)出現(xiàn)了明顯的冷平流,加之高空西風急流加強且風力強勁,盆地內(nèi)上升氣流明顯,這使得和田在5月12日開始發(fā)生了強沙塵天氣;隨后到沙塵天氣發(fā)展的后期圖5(c)和圖5(d),塔里木盆地內(nèi)氣流較之前而言有所穩(wěn)定,風力明顯下降.

    圖5 研究區(qū)沙塵過境前后500hPa位勢高度(單位:gpm;等值線)和溫度(單位:℃;顏色)的空間分布

    圖6 研究區(qū)沙塵過境前后500hPa風場圖

    圖中箭頭表示水平風場,顏色表示垂直風場,其中紅色代表下沉氣流且藍色表示上升氣流

    為了更加直觀了解到沙塵污染過程和田市低空天氣形勢狀況,結(jié)合圖7中850hPa風場示意圖可知,5月12日塔里木盆地南端上升氣流明顯且在當日達到了最強,強烈的西風受到了天山山脈阻擋與青藏高原北坡的大地形的影響,使其轉(zhuǎn)變?yōu)槠珫|風并通過了天山東部缺口,伴隨東灌大風氣流襲來,冷氣團中又包裹了大量的本地沙塵在偏北風與偏西風的作用下于羅布泊附近匯合并向塔里木盆地西南側(cè)輸送,使得局地溫度有所降低,盆地東西側(cè)溫差加大,而盆地內(nèi)部地形閉塞,阻擋環(huán)流,在上升氣流的作用下促進起沙,當天發(fā)生沙塵暴天氣,能見度迅速降低至373m,5月15日以后風力有所下降,天氣形勢逐漸穩(wěn)定,污染濃度水平也逐漸有所回落,并于5月18日結(jié)束.這次沙塵污染天氣過程中,和田市在春季處于熱低壓場中[49],并在西風背景環(huán)境與塔里木盆地周邊地形阻擋的影響下,使得污染物難以稀釋與擴散,沙塵氣團長時間在此地盤踞,繼而沙塵污染天氣過程持續(xù)時間較長.

    圖7 研究區(qū)沙塵過境前后850hPa風場圖

    圖中箭頭表示水平風場,顏色表示垂直風場,其中紅色代表下沉氣流且藍色表示上升氣流

    2.4 沙塵污染傳輸路徑分析

    為了解本次沙塵污染過程中顆粒物的來源以及傳輸路徑,本文采用HYSPLIT后向軌跡模型進行72h后向軌跡氣流的模擬,以和田市(37.12°N, 79.92°E)為擬受點,分別計算2022年5月12日(污染期間)與5月18日(污染結(jié)束)22:00 100m、1000m、2000m不同高度氣團中72h的后向軌跡,如圖8所示.從圖8a中的模擬結(jié)果可以看出,沙塵污染天氣發(fā)生前天氣處于靜穩(wěn)狀態(tài),從5月12日開始沙塵開始從地面迅速被揚至2000m的高空,表示沙塵天氣開始.由圖8a可見,100m、1000m以及2000m高度的氣團傳輸軌跡基本一致,本次天氣過境期間源自北疆的氣團隨著西風急流與偏東風在天山缺口處相遇并一同經(jīng)過吐哈盆地,近地面的沙塵被吹散至高空,在此期間攜帶了大量的污染物并伴隨著氣流的抬升運動被輸送至塔里木盆地的西南地區(qū),使得和田地區(qū)于5月12日開始發(fā)生沙塵天氣.根據(jù)圖8b可知, 1000m高度的氣團來源于阿富汗東北部,途徑了塔吉克斯坦;100m高度的氣團來源于準噶爾盆地,移動速度較慢,隨偏東風一同從天山缺口處涌入塔里木盆地,攜帶了大量來自于古爾班通古特沙漠與塔克拉瑪干沙漠的沙粒在和田發(fā)生沉降;而2000m高度的氣團來源于土庫曼斯坦的東部,移動速度較快,途徑烏茲別克斯坦、塔吉克斯坦與吉爾吉斯斯坦,最終從新疆的西部經(jīng)過喀什地區(qū)抵達盆地內(nèi)部,三種不同高度的氣團在5月18日呈現(xiàn)出了明顯下降的趨勢.HYSPLIT后向軌跡中不同高度氣團的傳輸路徑也與天氣形勢所得結(jié)果基本吻合.此次天氣污染過程中,污染物主要來自于新疆北部、新疆南部的塔克拉瑪干沙漠以及周邊的斯坦國家,這些地區(qū)工農(nóng)業(yè)活動較多,所排放出的污染物與沙塵氣溶膠是這次天氣污染過程的主要貢獻源[50],這也說明本次天氣污染過程不僅是污染物的遠距離輸送,也包含了大量的內(nèi)源性沙塵氣溶膠的排放貢獻.還可以發(fā)現(xiàn),模擬后的后向軌跡圖中不同高度氣團的傳輸路徑也與天氣形勢所得結(jié)果基本吻合.

    圖8 2022年5月12日22:00與5月18日22:00和田72h后向軌跡圖

    3 結(jié)論

    3.1 和田市在本次沙塵過境天氣期間PM10和PM2.5的變化趨勢較為一致且在垂直高度上存在非均一性,PM2.5/PM10均值為0.25,表明主要以粗顆粒物為主.污染物濃度的變化與沙塵天氣發(fā)展過程發(fā)展具有較好的時間一致性,風速與濕度的增加、溫度的降低有利于促成污染天氣的發(fā)生.

    3.2 沙塵污染期間出現(xiàn)了明顯的氣溶膠層,并且在垂直方向上有顯著的波動性變化,5月12日消光系數(shù)、后向散射系數(shù)和退偏振比達到了此次沙塵污染天氣中的最大值,分別為3.5km-1、0.07km-1·sr-1和0.17,在5月18日濃度開始降低,沙塵氣溶膠逐漸消散,沙塵污染天氣結(jié)束.

    3.3 西風受到天山山脈阻擋與青藏高原大地形的影響轉(zhuǎn)變?yōu)槠珫|風,在盆地缺口處與東灌氣流匯合,向西輸送過程中攜帶了大量的本地沙塵匯聚到了盆地西南部,和田市在春季處于熱低壓場中,并在西風背景環(huán)境與盆地周邊地形阻擋的影響下致使污染物難以稀釋與擴散,沙塵氣團長時間在此地盤踞,是本次沙塵污染持續(xù)較長的原因.

    3.4 后向軌跡模擬結(jié)果表明,本次沙塵污染天氣過程污染物主要來自于新疆北部、新疆南部的塔克拉瑪干沙漠以及周邊的斯坦國家,這些地區(qū)工農(nóng)業(yè)活動較多,所排放出的污染物與沙塵氣溶膠是這次天氣污染過程的主要貢獻源,表明本次天氣污染過程不僅是污染物的遠距離輸送,也包含了大量的內(nèi)源性沙塵氣溶膠的排放貢獻.

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

    Mao J T, Qian J H, Wang M H. Summary comment on research of atmospheric aerosl in CHINA [J]. Journal of Meteorology, 2002,5:625-634.

    [2] Miller R L,Tegen I, Perlwitz J. Surface radiative forcing by soil dust aerosols and the hydrologic cycle [J]. Journal of Geophysical Research, 2004,109(109):361-375.

    [3] Han Y, Dai X, Fang X, et al. Dust aerosols: a possible accelerant for an increasingly arid climate in North China [J]. Journal of arid environments, 2008,72(8):1476-1489.

    [4] Liu X, Huneeus N, Schulz M, et al. Global dust model intercomparison in AeroCom phase1 [J]. Atmospheric Chemistry and Physics, 2011,11(15):77-81.

    [5] Jemmett-Smith B C, Marsham J H, Knippertz P, et al. Quantifying global dust devil occurrence from meteorological analyses [J]. Geophysical Research Letters, 2015,42(4):1275-1282.

    [6] 成天濤,沈志寶.中國西北大氣沙塵光學特性的數(shù)值試驗 [J]. 高原氣象, 2001,20(3):291-297.

    Cheng T T, Shen Z B. A numerical simulation of optical characteristics for atmospheric dust aerosols in northwest China [J]. Journal of Meteorology, 2001,20(3):291-297.

    [7] Creamean J M, Suski K J, Rosenfeld D, et al. Dust and biological aerosols from the Sahara and Asia influence precipitation in the Western US [J]. Science, 2013,339(6127):1572-1578.

    [8] Carslaw K, Boucher O, Spracklen D, et al. A review of natural aerosol interactions and feedbacks within the earch system [J]. Atmospheric Chemistry and Physics, 2010,10(4):1-37.

    [9] Liu M, Westphal D L, Wang S G, et al. A high-resolution numerical study of the Asian dust storms of April 2001 [J]. Journal of Geophysical Research Atmospheres, 2003,34(6):618-622.

    [10] 韓永翔,奚曉霞,方小敏,等.亞洲大陸沙塵過程與北太平洋地區(qū)生物環(huán)境效應(yīng):以2001年4月中旬中亞特大沙塵暴為例 [J]. 科學通報, 2005,50(23):2649-2655.

    Han Y X, Xi X X, Fang X M, et al. Dust and sand processes in the Asian continent and bioenvironmental effects in the North Pacific. An example of the Central Asian Dust Storm of mid-April 2001 [J]. Science Bulletin, 2005,50(23):2649-2655.

    [11] 韓永翔,宋連春,趙天良,等.北太平洋地區(qū)沙塵沉降與海洋生物興衰的關(guān)系 [J]. 中國環(huán)境科學, 2006,26(2):157-160.

    Han Y X, Song L C, Zhao T L, et al. The relationship between continental dust and marine phytoplankton in the North Pacific [J]. China Environmental Science, 2006,26(2):157-160.

    [12] Ridgwell A J. Implications of the glacial CO2“iron hypothesis” for Quaternary climate change [J]. Geochemistry Geophysics Geosystems, 2013,4(9):130-145.

    [13] 欒兆鵬,趙天良,韓永翔,等.干旱半干旱地區(qū)塵卷風研究進展 [J].沙漠與綠洲氣象, 2016,56(2):5-12.

    Luan Z P, Zhao T L, Han Y X, et al. Advances in study of dust devils over arid and semi-arid regions [J]. Desert and Oasis Meteorology, 2016,56(2):5-12.

    [14] Goudie A S, Middleton N J. The changing frequency of dust storms through time [J]. Climate Change, 1992,20:197-225.

    [15] Hankin E H. On dust raising winds and descending currents [J]. India Met Memoirs, 1921,22(1):80-84.

    [16] Li H L, He Q, Liu X C. Identification of long-range transport pathways and potential source regions of PM2.5and PM10at Akedala Station, Central Asia [J]. Atmosphere, 2020,11:1183-1200.

    [17] Saunois M, Stavert A R, Poulter B, et al. The global methane budget 2000~2017 [J]. Earth System Science Data, 2020,12(3):1561-1623.

    [18] Raipal D K, Deka S N. ANDHI, the convective dust storm of northwest India [J]. Mausam, 1980,31:431-442.

    [19] Huang J, Minnis P, Chen B, et al. Long-range transport and vertical structure of Asian dust from CALIPSO and surface measurements during PACDEX [J]. Journal of Geophysical Research: Atmospheres, 2008,113:D23212.

    [20] 毛卓成,許建明,楊丹丹,等.上海地區(qū)PM2.5-O3復合污染特征及氣象成因分析 [J]. 中國環(huán)境科學, 2019,39(7):2730-2738.

    Mao Z C, Xu J M, Yang D D, et al. Analysis of characteristics and meteorological causes of PM2.5-O3compound pollution in Shanghai [J]. China Environmental Science, 2019,39(7):2730-2738.

    [21] 王躍思,姚 利,王莉莉,等.2013年元月我國中東部地區(qū)強霾污染成因分析 [J]. 中國科學:地球科學, 44(1):15-26.

    Wang Y S, Yao L, Wang L L, et al. 2014. Mechanism for the formation of the January 2013 heavy haze pollution episode over central and eastern China. Science China: Earth Sciences, 57:14–25.

    [22] 周自江,王錫穩(wěn),牛若蕓.近47年中國沙塵暴氣候特征研究 [J]. 應(yīng)用氣象學報, 2002,13(2):193-200.

    Zhou Z J, Wang X W, Niu R Y. Climate characteristics of sandstorm in CHINA in recent 47years [J]. Journal of Applied Meteorology, 2002,13(2):193-200.

    [23] 錢正安,蔡 英,劉景濤,等.中國北方沙塵暴研究的若干進展和問題 [C]//中國氣象學會.新世紀氣象科技創(chuàng)新與大氣科學發(fā)展——中國氣象學會2003年年會“氣候系統(tǒng)與氣候變化”分會論文集.氣象出版社, 2003:265-269.

    Qian Z A, Cai Y, Liu J T, et al. Some progress and problems in the study of sand and dust storms in northern China [C]//Chinese Meteorological Society. Meteorological Science and Technology Innovation and Atmospheric Science Development in the New Century - Proceedings of the 2003Annual Meeting of the Chinese Meteorological Society on "Climate System and Climate Change". Meteorological Publishing House, 2003:265-269.

    [24] 劉曉東,田 良,張小曳.塔克拉瑪干沙塵活動對下游大氣PM10濃度的影響 [J]. 中國環(huán)境科學, 2004,(5):17-21.

    Liu X D, Tian L, Zhang X Y. Influence of spring dust activities over the Taklimakan Desert area on concentrations of atmospheric PM10in east of Northwest China [J]. China Environmental Science, 2004,(5):17-21.

    [25] 艾沙江·艾力,徐至遠,徐海量,等.基于軌跡模型的和田綠洲沙塵移動軌跡及氣象特征分析[J]. 環(huán)境工程技術(shù)學報, 2022,12(4):1007-1014.

    Aishajiang A L, Xu Z Y, Xu H L, et al. Analysis of the origin and meteorological characteristics of dust-storm in Hetian Oasis based on a trajectory model approach [J]. Journal of Environmental Engineering Technology, 2022,12(4):1007-1014.

    [26] Li J L, He Q, Ge X. Spatiotemporal distribution of aerosols over the Tibet Plateau and Tarim Basin (1980~2020) [J]. Journal of Cleaner Production, 2022,374:133958.

    [27] 畢道金,何 清,李京龍,等.基于無人機垂直觀測塔克拉瑪干沙漠一次沙塵污染過程研究 [J]. 環(huán)境科學學報, 2022,42(12):298-310.

    Bi D J, He Q, Li J L, et al. Study on a dust pollution process in Taklamakan Desert based on vertical observation by UAV [J]. Acta Scientiae Circumstantiae, 2022,42(12):298-310.

    [28] Li J L, He Q, Jin L, et al. Three-dimensional distribution of dust aerosols over the Tarim Basin and the Tibet Plateau during 2007~2021 derived from CALIPSO lidar observations [J]. Journal of Cleaner Production, 2023,400:136746.

    [29] 劉新春,鐘玉婷,何 清,等.塔克拉瑪干沙漠腹地及周邊地區(qū)PM10時空變化特征及影響因素分析 [J]. 中國沙漠, 2011,31(2):323-330.

    Liu X C, Zhong Y T, He Q, et al. Spatio-temporal pattern of PM10concentration and impact factors in the hinterland and surrounding area of Taklimakan desert [J]. China Desert, 2011,31(2):323-330.

    [30] 張小嘯,陳 曦,王自發(fā),等.新疆和田綠洲大氣降塵和PM10濃度變化特征分析 [J]. 干旱區(qū)地理, 2015,38(3):454-462.

    Zhang X X, Chen x, Wang Z F, et al. Variation characteristics of atmospheric dustfall and PM10concentration in Hotan oasis, Xinjiang [J]. ARID LAND GEOGRAPHY, 2015,38(3):454-462.

    [31] 劉文彬,劉 濤,黃祖照,等.利用偏振-米散射激光雷達研究廣州一次浮塵天氣過程 [J]. 中國環(huán)境科學, 2013,33(10):1751-1757.

    Liu W B, Liu T, Huang Z Z, Study of a dust case in Guangzhou using polarization mie-scattering lidar. [J]. China Environmental Science, 2013,33(10):1751-1757.

    [32] 劉文彬,黃祖照,陳彥寧,等.廣州市春季一次沙塵天氣過程綜合觀測 [J]. 中國環(huán)境監(jiān)測, 2017,33(5):42-48.

    Liu W B, Huang Z X, Chen Y N, Observation Studies on a Dust Case in the Spring of Guangzho [J]. China Environmental Monitoring, 2017,33(5):42-48.

    [33] 黃忠偉.氣溶膠物理光學特性的激光雷達遙感研究[D]. 蘭州:蘭州大學, 2012:19-70.

    Huang Z W. Lidar remote sensing study of aerosol physico-optical properties [D]. Lanzhou: Lanzhou University, 2012:19-70.

    [34] Liao T T, Wang S, Ai J, et al. Heavy pollution episodes, transport pathways and potential sources of PM2.5, during the winter of 2013 in Chengdu (China) [J]. Science of the Total Environment, 2017,584-585:1056-1065.

    [35] 王芳龍,李忠勤,尤曉妮,等.2015~2017年天水市大氣污染物變化特征及來源分析 [J]. 環(huán)境科學學報, 2018,38(12):4592-4604.

    Wang F L, Li Z Q, You X N, et al. Variation characteristics and source analysis of atmospheric pollutants in Tianshui from 2015 to 2017 [J]. Acta Scientiae Circumstantiae, 2018,38(12):4592-4604.

    [36] Alam K, Qureshi S, Blaschke T. Monitoring spatio-temporal aerosol patterns over Pakistan based on MODIS, TOMS and MISR satellite data and a HYSPLIT model [J]. Atmospheric environment, 2011,45 (27):4641-4651.

    [37] Rolph G, Stein A, Stunder B. Real-time environmental applications and display system:READY [J]. Environmental Modelling & Software, 2017,95:210-228.

    [38] 中央氣象局.地面氣象觀測規(guī)范 [M]. 北京:氣象出版社, 1979:21-27.

    Central Meteorological Administration. Specification for surface meteorological observations [M]. Beijing: Meteorological Press, 1979,21-27.

    [39] GB 3095-2012 環(huán)境空氣質(zhì)量標準 [S]. GB 3095-2012 Ambient air quality standards [S].

    [40] Zhou C, Liu Y, He Q, et al. Dust Characteristics Observed by Unmanned Aerial Vehicle over the Taklimakan Desert. Remote Sensing. 2022:14(4):990.

    [41] 曹云擎,王體健,高麗波,等.基于無人機垂直觀測的南京PM2.5污染個例研究 [J]. 氣候與環(huán)境研究, 2020,25(3):292-304.

    Cao Y Q, Wang T J, Gao L B, et al. A Case Study of PM2.5Pollution in Nanjing Based on Unmanned Aerial Vehicle Vertical Observations [J]. Climatic and Environmental Research, 2020,25(3):292?304.

    [42] 張 喆,丁建麗,王瑾杰.中亞沙塵氣溶膠時空分布特征及潛在擴散特性分析[J]. 地理學報, 2017,72(3):507-520.

    Zhang Z, Ding J L, Wang J J. Spatio-temporal variations and potential diffusion characteristics of dust aerosol originating from Central Asia [J]. Journal of Geography, 2017,72(3):507-520.

    [43] 楊富燕,張 寧,朱蓮芳,等.基于激光雷達和微波輻射計觀測確定混合層高度方法的比較 [J]. 高原氣象, 2016,35(4):1102-1111.

    Yang F Y, Zhang N, Zhu L F, et al. Comparison of the mixing layer height determination methods using Lidar and Microwave Radiometer [J]. Plateau Meteorology, 2016,35(4):1102-1111.

    [44] 郝巨飛,張功文,王曉娟,等.一次環(huán)境大氣重污染過程的監(jiān)測分析 [J]. 高原氣象, 2017,36(5):1404-1411.

    Hao J F, Zhang G W, Wang X J, et al. Analysis of meteorology detection data during a heavy pollution event [J]. Plateau Meteorology, 2017,36(5):1404-1411.

    [45] 高慶先,李令軍,張運剛,等.我國春季沙塵暴研究 [J]. 中國環(huán)境科學, 2000,20(6):495-500.

    Gao Q X, Li L J, Zhang Y G, Studies on the springtime dust storm of China [J]. 中國環(huán)境科學, 2000,20(6):495-500.

    [46] 云靜波,姜學恭,孟雪峰,等.冷鋒型和蒙古氣旋型沙塵暴過程若干統(tǒng)計特征的對比分析[J]. 高原氣象, 2013,32(2):423-434.

    Yun J B, Jiang X G, Meng X F, et al. Comparative Analyses on Some Statistic Characteristics between Cold Front and Mongolia Cyclone Duststorm Processes, [J]. Plateau Meteorology, 2013,32(2):423-434.

    [47] 宗志平,張恒德,馬 杰.2009年4月下旬蒙古氣旋型大范圍沙塵暴天氣過程的診斷分析[J]. 沙漠與綠洲氣象, 2012,6(1):1-9.

    Zong Z P, Zhang Z D, Ma J. Diagnosis and Analysis of a Large-scale Sand-dust Storm Process Resulted from Mongolia Cyclone during Late April 2009 [J]. Desert and Oasis Meteorology, 2012,6(1):1-9.

    [48] 康 林,季明霞,黃建平.歐亞大氣環(huán)流對中國北方夏季沙塵天氣的影響 [J]. 中國沙漠, 2014,28(5):968-975.

    Kang L, Ji M X, Huang J P. lmpact of eurasian atmospheric circulation on aust events in Summer over Northern China , journa of Desert Research, 2014,34(1):197-205.

    [49] 韓永翔,方小敏,宋連春,等.塔里木盆地中的大氣環(huán)流及沙塵暴成因探討——根據(jù)沙漠風積地貌和氣象觀測重建的風場 [J]. 大氣科學, 2005,(4):627-635.

    Han Y X, Fang X M, Song L C, et al. A study of At mospheric circulation and dust storm causes of for mation in the Tarim Basin -the restructured wind field by shapes of dune and observed prevailing wind [J]. Atmospheric Sciences, 2005,(4):627-635.

    [50] 薛福民,李 娟,黃 侃,等.塔克拉瑪干沙漠黑碳氣溶膠的特性及來源[J]. 中國科學:化學, 2010,40(5):556-566.

    Xue F M, Li J, Huang K, et al. Characteristics and source of black carbon over Taklimakan Desert [J]. Chinese Science: Chemistry, 2010, 40(5):556-566.

    Ground-based ridar-based investigation of a sand and dust pollution process in Hotan City.

    YIN Lu-lu1,2, HE Qing1,2*, LI Jing-long2,3, MENG Lu2,3, FU Guang-xiang1,2, WU Cai-yun2,3, XIAO Hong-dan2, FAN Xu-yu-chen2

    (1.College of Geographic Science and Tourism, Xinjiang Normal University, Urumqi 830054,China;2.Institute of Desert Meteorology, China Meteorological Administration/National Observation and Research Station of Desert Meteorology, Taklimakan Desert of Xinjiang/Taklimakan Desert Meteorology Field Experiment Station of China Meteorological Administration/Xinjiang Key Laboratory of Desert Meteorology and Sandstorm, Urumqi 830002, China;3.College of Geography and Remote Sensing Sciences, Xinjiang University, Urumqi 800017, China)., 2023,43(12):6290~6300

    In May 2022, Hotan City,Xinjiang suffered the strongest sandstorm in nearly five years, and the sand and dust pollution process lasted for a week. In this paper, the changes of particulate matter concentration and the vertical distribution characteristics of aerosol optical information in the process of sand and dust pollution were analyzed by using atmospheric aerosol lidar, pollutant concentration of environmental monitoring station and near-ground meteorological observation data, and the meteorological matching pattern of this polluted weather was obtained by using global reanalysis data, and finally the transport and potential sources of air mass during transit were analyzed by combining HYSPLIT backward trajectory model. The results showed that during the period of sand and dust pollution, the trend of PM2.5and PM10was consistent and there was heterogeneity in vertical height, and the average PM2.5/PM10was 0.25, indicating that coarse particulate matter was the main thing, and the increase of wind speed, temperature and humidity were conducive to the occurrence of polluted weather. During the dust pollution, there was an obvious aerosol layer, and there were significant fluctuations in the vertical direction, and the extinction coefficient, backscatter coefficient and depolarization ratio reached the maximum values of the process, which were 3.5km-1, 0.07km-1·sr-1and 0.17, respectively. In the spring, Hotan City is in a hot and low-pressure field, and it is difficult to dilute and diffuse pollutants under the background environment of the westerly wind and the surrounding terrain of the Tarim Basin, and the sand and dust air masses have been entrenched here for a long time, which is the reason for the long duration of this pollution. The pollutants in this sand and dust weather process mainly came from the Taklamakan Desert in northern Xinjiang, southern Xinjiang and the surrounding countries of stan, and the pollution sources were mainly the emission contribution of endogenous sand and dust aerosols.

    Hotan City;dust aerosols;vertical distribution;optical information

    X513

    A

    1000-6923(2023)12-6290-11

    陰璐璐,何 清,李京龍,等.基于地基雷達探究和田市一次沙塵污染過程 [J]. 中國環(huán)境科學, 2023,43(12):6290-6300.

    Yin L L, He Q, Li J L, et al. Ground-based ridar-based investigation of a sand and dust pollution process in Hotan City [J]. China Environmental Science, 2023,43(12):6290-6300.

    2023-04-14

    第三次新疆綜合科學考察項目(2021xjkk030501)

    * 責任作者, 研究員, qinghe@idm.cn

    陰璐璐(1996-),女,新疆吐魯番人,新疆師范大學碩士研究生,主要從事干旱區(qū)環(huán)境演變與災害防治研究.1498240373@qq.com.

    感謝中國氣象局烏魯木齊沙漠氣象研究所提供的雷達實測數(shù)據(jù),以及氣象局相關(guān)工作人員對于儀器的悉心維護.感謝歐洲中期天氣預報中心(ECMWF)提供的全球氣候的第五代大氣再分析數(shù)據(jù)集.同時也感謝美國國家環(huán)境預報中心(NCEP)提供的同期全球資料同化系統(tǒng)(GDAS)數(shù)據(jù).

    猜你喜歡
    散射系數(shù)沙塵沙塵暴
    什么是沙塵天氣?
    國土綠化(2024年3期)2024-04-17 01:02:22
    等離子體層嘶聲波對輻射帶電子投擲角散射系數(shù)的多維建模*
    物理學報(2022年22期)2022-12-05 11:16:04
    可怕的沙塵暴
    北部灣后向散射系數(shù)的時空分布與變化分析
    可怕的沙塵天氣
    學生天地(2017年12期)2017-05-17 05:51:25
    大話西游之沙塵暴
    天外來客:火星沙塵暴
    一維帶限Weierstrass分形粗糙面電磁散射的微擾法研究
    基于PM譜的二維各向異性海面電磁散射的微擾法研究
    如何看待沙塵暴的利與弊
    河南科技(2014年1期)2014-02-27 14:04:33
    欧美中文日本在线观看视频| 人人妻人人澡欧美一区二区| 久久中文看片网| 亚洲熟妇熟女久久| www.999成人在线观看| 国产免费av片在线观看野外av| 一进一出好大好爽视频| 日韩精品免费视频一区二区三区| 国产aⅴ精品一区二区三区波| 午夜福利成人在线免费观看| 国产片内射在线| 久久久久久久久免费视频了| 亚洲国产日韩欧美精品在线观看 | 美女午夜性视频免费| 国产精品精品国产色婷婷| 日韩有码中文字幕| 亚洲av片天天在线观看| 女人高潮潮喷娇喘18禁视频| 性色av乱码一区二区三区2| 亚洲成人久久爱视频| 淫妇啪啪啪对白视频| 变态另类成人亚洲欧美熟女| 一区二区日韩欧美中文字幕| 高清毛片免费观看视频网站| 国产成人av激情在线播放| 亚洲人成网站在线播放欧美日韩| 深夜精品福利| 亚洲成人久久爱视频| 99精品在免费线老司机午夜| 国产高清视频在线播放一区| 成人亚洲精品一区在线观看| ponron亚洲| 男女做爰动态图高潮gif福利片| 男女之事视频高清在线观看| 成人三级做爰电影| 久久久久免费精品人妻一区二区 | 久久久久久国产a免费观看| 中亚洲国语对白在线视频| 禁无遮挡网站| 国产高清videossex| 婷婷亚洲欧美| 久久欧美精品欧美久久欧美| 精品少妇一区二区三区视频日本电影| 久久九九热精品免费| 国产成人影院久久av| 麻豆成人午夜福利视频| 国产极品粉嫩免费观看在线| 狂野欧美激情性xxxx| 老司机福利观看| 91老司机精品| 国产熟女xx| 久热爱精品视频在线9| 欧美亚洲日本最大视频资源| 91成人精品电影| 欧美激情 高清一区二区三区| 黄色视频不卡| 国产av不卡久久| 久久久久久国产a免费观看| 波多野结衣高清作品| 欧美性长视频在线观看| 色老头精品视频在线观看| 999久久久精品免费观看国产| 18禁观看日本| 久热这里只有精品99| 美女高潮到喷水免费观看| 最近最新中文字幕大全免费视频| 久久性视频一级片| 性欧美人与动物交配| 日韩欧美国产在线观看| 波多野结衣巨乳人妻| 国产黄片美女视频| 国内少妇人妻偷人精品xxx网站 | 久久精品成人免费网站| 男女之事视频高清在线观看| 999久久久精品免费观看国产| 亚洲av片天天在线观看| 欧美一级a爱片免费观看看 | 亚洲欧洲精品一区二区精品久久久| 看片在线看免费视频| 欧美国产精品va在线观看不卡| 99国产极品粉嫩在线观看| or卡值多少钱| 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 亚洲,欧美精品.| 欧美成狂野欧美在线观看| xxx96com| 黄色毛片三级朝国网站| 欧美乱色亚洲激情| 中文亚洲av片在线观看爽| 神马国产精品三级电影在线观看 | 长腿黑丝高跟| 中亚洲国语对白在线视频| 午夜久久久久精精品| 成人特级黄色片久久久久久久| 88av欧美| 免费在线观看成人毛片| 国产熟女xx| 亚洲真实伦在线观看| 日本五十路高清| 久久午夜综合久久蜜桃| 人人澡人人妻人| 两性午夜刺激爽爽歪歪视频在线观看 | 麻豆成人午夜福利视频| 国产精品久久电影中文字幕| 9191精品国产免费久久| 黄色成人免费大全| 老汉色∧v一级毛片| 一级a爱片免费观看的视频| 国产欧美日韩精品亚洲av| 两个人视频免费观看高清| 国产精品美女特级片免费视频播放器 | 黄网站色视频无遮挡免费观看| 午夜影院日韩av| 国产一区二区激情短视频| 禁无遮挡网站| 美女扒开内裤让男人捅视频| 免费一级毛片在线播放高清视频| 变态另类成人亚洲欧美熟女| 日韩成人在线观看一区二区三区| 久久精品影院6| 嫩草影院精品99| 国产成人影院久久av| 免费看日本二区| 真人做人爱边吃奶动态| 日韩中文字幕欧美一区二区| 亚洲国产看品久久| 91成年电影在线观看| 国产色视频综合| 国产v大片淫在线免费观看| 国产精品电影一区二区三区| 麻豆成人av在线观看| 97人妻精品一区二区三区麻豆 | 51午夜福利影视在线观看| 久久久久国内视频| 亚洲精品在线观看二区| 亚洲国产精品999在线| 国产av一区二区精品久久| 一本久久中文字幕| 人人妻人人看人人澡| 亚洲美女黄片视频| 国产亚洲av嫩草精品影院| 亚洲av中文字字幕乱码综合 | 色尼玛亚洲综合影院| 亚洲一区二区三区色噜噜| 欧美性猛交黑人性爽| 国产v大片淫在线免费观看| 变态另类成人亚洲欧美熟女| 国产一区二区三区在线臀色熟女| 午夜久久久久精精品| 国产精品综合久久久久久久免费| 18禁美女被吸乳视频| 色尼玛亚洲综合影院| 99久久久亚洲精品蜜臀av| 黄色视频,在线免费观看| 午夜福利成人在线免费观看| 在线av久久热| www国产在线视频色| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看 | 免费高清在线观看日韩| 国产av一区在线观看免费| 丝袜在线中文字幕| 中文字幕另类日韩欧美亚洲嫩草| 中文在线观看免费www的网站 | 亚洲av日韩精品久久久久久密| 性欧美人与动物交配| 国产亚洲av高清不卡| 琪琪午夜伦伦电影理论片6080| 国产单亲对白刺激| 亚洲成人久久性| 欧美激情高清一区二区三区| 美女扒开内裤让男人捅视频| 男女那种视频在线观看| 欧美国产精品va在线观看不卡| 午夜两性在线视频| 国产在线观看jvid| 欧美性猛交黑人性爽| 国产熟女午夜一区二区三区| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 熟女电影av网| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 欧美人与性动交α欧美精品济南到| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av| 一a级毛片在线观看| 国产亚洲精品久久久久5区| 亚洲人成伊人成综合网2020| 国产1区2区3区精品| 人妻丰满熟妇av一区二区三区| 级片在线观看| 在线视频色国产色| 国产欧美日韩一区二区精品| 曰老女人黄片| 亚洲男人天堂网一区| 免费在线观看成人毛片| 久久久久久免费高清国产稀缺| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看 | 国产一区二区三区在线臀色熟女| 亚洲人成网站高清观看| 亚洲欧美日韩无卡精品| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色成人免费大全| 免费高清在线观看日韩| 精品国产乱子伦一区二区三区| 美女大奶头视频| a在线观看视频网站| netflix在线观看网站| 久久99热这里只有精品18| 国产国语露脸激情在线看| 成年版毛片免费区| 日韩视频一区二区在线观看| 淫妇啪啪啪对白视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久香蕉激情| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文av在线| 午夜久久久在线观看| 国内少妇人妻偷人精品xxx网站 | 精品国产乱码久久久久久男人| 观看免费一级毛片| 一区福利在线观看| 久久久国产成人精品二区| 久久精品91蜜桃| 少妇裸体淫交视频免费看高清 | 精品卡一卡二卡四卡免费| 少妇被粗大的猛进出69影院| 精品午夜福利视频在线观看一区| 成人国产一区最新在线观看| 亚洲色图 男人天堂 中文字幕| 欧美 亚洲 国产 日韩一| 欧美黑人欧美精品刺激| 在线免费观看不下载黄p国产| 黄色欧美视频在线观看| 小说图片视频综合网站| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 亚洲美女搞黄在线观看 | 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 男人的好看免费观看在线视频| 亚洲,欧美,日韩| 精品午夜福利在线看| 特大巨黑吊av在线直播| 国产av在哪里看| 国产亚洲av嫩草精品影院| 亚洲最大成人中文| 黄色视频,在线免费观看| 乱码一卡2卡4卡精品| 精品人妻视频免费看| 亚洲欧美日韩高清专用| 国产高清不卡午夜福利| 日韩av在线大香蕉| 欧美激情国产日韩精品一区| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品一区二区| 亚洲精品色激情综合| 午夜精品国产一区二区电影 | 国产一区二区在线观看日韩| 亚洲一区高清亚洲精品| 亚洲成人久久性| 久久精品久久久久久噜噜老黄 | 亚洲av熟女| 天堂av国产一区二区熟女人妻| 日本 av在线| 97超级碰碰碰精品色视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国内精品一区二区在线观看| eeuss影院久久| 嫩草影院入口| 国产高清三级在线| 99久久久亚洲精品蜜臀av| a级毛片免费高清观看在线播放| 在线看三级毛片| 日本熟妇午夜| 久久精品国产亚洲av香蕉五月| 久久九九热精品免费| 午夜免费男女啪啪视频观看 | 免费看日本二区| 婷婷精品国产亚洲av| 欧美精品国产亚洲| 日韩 亚洲 欧美在线| 不卡视频在线观看欧美| 久久这里只有精品中国| 国产精品嫩草影院av在线观看| 亚洲人成网站在线观看播放| 久久人人精品亚洲av| 亚洲七黄色美女视频| 尤物成人国产欧美一区二区三区| 人妻夜夜爽99麻豆av| 日韩强制内射视频| 精品午夜福利视频在线观看一区| 国产高清不卡午夜福利| 成年女人永久免费观看视频| 国产亚洲欧美98| 精品福利观看| 亚洲av免费在线观看| 最近的中文字幕免费完整| 天天一区二区日本电影三级| 99久久精品热视频| 最近在线观看免费完整版| 国产亚洲精品综合一区在线观看| 成人三级黄色视频| 日韩高清综合在线| 国产精品野战在线观看| 色在线成人网| 最后的刺客免费高清国语| 欧美潮喷喷水| 大型黄色视频在线免费观看| 我要搜黄色片| 晚上一个人看的免费电影| 婷婷精品国产亚洲av| 亚洲自偷自拍三级| 日本撒尿小便嘘嘘汇集6| 欧美激情久久久久久爽电影| 在线免费观看的www视频| 午夜激情福利司机影院| 又黄又爽又免费观看的视频| 国产私拍福利视频在线观看| 久久精品夜色国产| 大香蕉久久网| 国产男人的电影天堂91| 国产单亲对白刺激| 日本黄色视频三级网站网址| 乱人视频在线观看| 欧美国产日韩亚洲一区| 亚洲成人久久性| 少妇裸体淫交视频免费看高清| 尤物成人国产欧美一区二区三区| 国内少妇人妻偷人精品xxx网站| 狠狠狠狠99中文字幕| 又爽又黄a免费视频| 人妻夜夜爽99麻豆av| 一夜夜www| 国产色爽女视频免费观看| 男人舔奶头视频| 欧美成人a在线观看| 大又大粗又爽又黄少妇毛片口| 欧美日韩在线观看h| 一级毛片我不卡| 日韩在线高清观看一区二区三区| 精品国产三级普通话版| 99在线视频只有这里精品首页| 特级一级黄色大片| av免费在线看不卡| 国产黄色视频一区二区在线观看 | 亚洲欧美中文字幕日韩二区| 日韩精品中文字幕看吧| 黑人高潮一二区| 色视频www国产| 日本黄色视频三级网站网址| 男女边吃奶边做爰视频| 欧美极品一区二区三区四区| 国产精品一区二区性色av| 日韩av在线大香蕉| 99久久九九国产精品国产免费| av在线播放精品| 国产精品乱码一区二三区的特点| 亚洲av熟女| 国内精品一区二区在线观看| av卡一久久| 大香蕉久久网| 97热精品久久久久久| av福利片在线观看| 日本一二三区视频观看| 国产免费一级a男人的天堂| 老司机影院成人| av在线蜜桃| 欧美色视频一区免费| 国产极品精品免费视频能看的| 中国美女看黄片| 尾随美女入室| 男人舔奶头视频| 国产高清视频在线观看网站| 干丝袜人妻中文字幕| 国产精品久久久久久亚洲av鲁大| 桃色一区二区三区在线观看| 亚洲精品日韩av片在线观看| 亚洲美女搞黄在线观看 | 精品久久久久久久久久免费视频| 日本精品一区二区三区蜜桃| 成人欧美大片| 噜噜噜噜噜久久久久久91| 亚洲国产日韩欧美精品在线观看| 国产黄色小视频在线观看| 日韩国内少妇激情av| 一夜夜www| 国产毛片a区久久久久| 久久韩国三级中文字幕| 欧美区成人在线视频| 最后的刺客免费高清国语| 在线天堂最新版资源| 亚洲欧美清纯卡通| 12—13女人毛片做爰片一| 午夜爱爱视频在线播放| 人人妻,人人澡人人爽秒播| 欧美色欧美亚洲另类二区| av中文乱码字幕在线| 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 久久草成人影院| www日本黄色视频网| 99久国产av精品| 久久热精品热| 国产亚洲精品久久久com| 日本-黄色视频高清免费观看| 91久久精品国产一区二区成人| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 国产精品不卡视频一区二区| 69人妻影院| 婷婷亚洲欧美| 男人的好看免费观看在线视频| 三级国产精品欧美在线观看| 秋霞在线观看毛片| 亚洲成人中文字幕在线播放| 国产乱人偷精品视频| 免费电影在线观看免费观看| 精品国内亚洲2022精品成人| 成人漫画全彩无遮挡| 一个人看视频在线观看www免费| 欧美zozozo另类| 在线a可以看的网站| 一本久久中文字幕| 精品人妻偷拍中文字幕| 夜夜夜夜夜久久久久| 3wmmmm亚洲av在线观看| 老师上课跳d突然被开到最大视频| 国产高清视频在线观看网站| 国产精品人妻久久久久久| 午夜亚洲福利在线播放| 一级黄色大片毛片| 麻豆成人午夜福利视频| 亚洲乱码一区二区免费版| 欧美成人免费av一区二区三区| 久久精品久久久久久噜噜老黄 | 在线免费观看不下载黄p国产| 最好的美女福利视频网| 性插视频无遮挡在线免费观看| 99热这里只有精品一区| 国产一区二区三区av在线 | 成年女人毛片免费观看观看9| 22中文网久久字幕| 免费观看的影片在线观看| 少妇高潮的动态图| 欧美日韩乱码在线| 赤兔流量卡办理| 亚洲av成人av| 精品久久久久久久久久免费视频| 久久精品国产亚洲av香蕉五月| 亚洲精品456在线播放app| ponron亚洲| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| 成人午夜高清在线视频| 欧美日韩国产亚洲二区| 综合色av麻豆| 国产激情偷乱视频一区二区| 99热只有精品国产| 亚洲图色成人| 国产不卡一卡二| 成人鲁丝片一二三区免费| 亚洲av中文av极速乱| 国产人妻一区二区三区在| 精品一区二区三区人妻视频| 一边摸一边抽搐一进一小说| 99久久九九国产精品国产免费| 午夜亚洲福利在线播放| 国产精品不卡视频一区二区| 成年女人永久免费观看视频| 国内精品久久久久精免费| av免费在线看不卡| 日本三级黄在线观看| 深夜a级毛片| 最新在线观看一区二区三区| 久久精品人妻少妇| 国产真实乱freesex| 亚洲中文字幕一区二区三区有码在线看| 久久天躁狠狠躁夜夜2o2o| 国产精品一区二区三区四区免费观看 | 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 免费观看的影片在线观看| 国产视频一区二区在线看| 亚洲五月天丁香| 波野结衣二区三区在线| videossex国产| av福利片在线观看| 高清毛片免费观看视频网站| 尾随美女入室| 中国国产av一级| 成人精品一区二区免费| 国产一区二区在线av高清观看| 美女被艹到高潮喷水动态| 亚洲熟妇中文字幕五十中出| 最近最新中文字幕大全电影3| 国产蜜桃级精品一区二区三区| 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 精品午夜福利视频在线观看一区| 亚洲成人中文字幕在线播放| 欧美成人精品欧美一级黄| 搡老岳熟女国产| 免费看日本二区| 美女大奶头视频| 老司机午夜福利在线观看视频| 久久国产乱子免费精品| 老司机午夜福利在线观看视频| 亚洲国产精品sss在线观看| 欧美色欧美亚洲另类二区| 亚洲乱码一区二区免费版| 麻豆av噜噜一区二区三区| 亚洲美女黄片视频| 成人鲁丝片一二三区免费| 色在线成人网| 性欧美人与动物交配| 美女大奶头视频| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩精品成人综合77777| 99久国产av精品| 欧美bdsm另类| 亚洲国产精品久久男人天堂| 菩萨蛮人人尽说江南好唐韦庄 | 国产午夜精品久久久久久一区二区三区 | 精品午夜福利视频在线观看一区| www日本黄色视频网| 成年av动漫网址| 99热网站在线观看| 丝袜喷水一区| 变态另类丝袜制服| 真人做人爱边吃奶动态| 中国国产av一级| 久久精品国产亚洲av涩爱 | 毛片女人毛片| 国产高清视频在线播放一区| 51国产日韩欧美| 蜜臀久久99精品久久宅男| a级毛色黄片| 国产一区二区三区av在线 | 尤物成人国产欧美一区二区三区| 成人精品一区二区免费| 久久中文看片网| 成人美女网站在线观看视频| АⅤ资源中文在线天堂| 免费一级毛片在线播放高清视频| 成人特级av手机在线观看| 97超碰精品成人国产| 免费观看人在逋| 国产精品一区二区三区四区免费观看 | 欧美另类亚洲清纯唯美| 久久精品国产亚洲av香蕉五月| 级片在线观看| 亚洲熟妇中文字幕五十中出| 日韩一区二区视频免费看| 中文字幕熟女人妻在线| 内射极品少妇av片p| 日韩av在线大香蕉| 小说图片视频综合网站| 久久九九热精品免费| 午夜福利高清视频| 又爽又黄无遮挡网站| 搡女人真爽免费视频火全软件 | 人妻夜夜爽99麻豆av| 三级毛片av免费| 国产精品无大码| 成年av动漫网址| av免费在线看不卡| 啦啦啦韩国在线观看视频| 久久久久久大精品| 久久久久久久久久成人| 91在线精品国自产拍蜜月| 少妇熟女aⅴ在线视频| 69av精品久久久久久| 欧美日韩综合久久久久久| 亚洲四区av| 两个人的视频大全免费| 99热全是精品| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕| 亚洲一级一片aⅴ在线观看| 大又大粗又爽又黄少妇毛片口| 丰满乱子伦码专区| 国产探花极品一区二区| 国产精品久久久久久av不卡| 成年女人看的毛片在线观看| 欧美xxxx性猛交bbbb| 搡老岳熟女国产| 乱码一卡2卡4卡精品| 成人特级av手机在线观看| 亚洲,欧美,日韩| 天天躁日日操中文字幕| 夜夜看夜夜爽夜夜摸| 好男人在线观看高清免费视频| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久久久免| eeuss影院久久| 老女人水多毛片| 麻豆国产97在线/欧美| 精品少妇黑人巨大在线播放 | 国产视频一区二区在线看| 国产人妻一区二区三区在| 国产真实伦视频高清在线观看| 国产精品久久电影中文字幕| 美女黄网站色视频| 亚洲三级黄色毛片| 亚洲精品国产成人久久av| 伊人久久精品亚洲午夜| 一级毛片我不卡| 大型黄色视频在线免费观看| 中出人妻视频一区二区| 99热网站在线观看| 真人做人爱边吃奶动态|