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

    臨安本底站2010~2015年P(guān)M10污染特征及影響因素分析

    2017-01-23 21:08:37李金娟馬千里
    中國(guó)環(huán)境科學(xué) 2017年8期
    關(guān)鍵詞:氣團(tuán)臨安能見(jiàn)度

    岳 毅,李金娟,馬千里

    ?

    臨安本底站2010~2015年P(guān)M10污染特征及影響因素分析

    岳 毅1,2,李金娟1*,馬千里3

    (1.貴州大學(xué)資源與環(huán)境工程學(xué)院,貴州貴陽(yáng) 550025;2.臨安市氣象局,浙江臨安 311300;3.臨安區(qū)域大氣本底站,浙江臨安 311307)

    利用長(zhǎng)江三角洲區(qū)域代表性站點(diǎn)臨安區(qū)域大氣本底站2010~2015年的PM10質(zhì)量濃度和各氣象要素資料,分析了PM10質(zhì)量濃度的變化特征及其影響因素.結(jié)果表明:與《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB3096-2012)二級(jí)標(biāo)準(zhǔn)相比,2010~2015年P(guān)M10質(zhì)量濃度達(dá)標(biāo)率為96.17%;日變化中峰值出現(xiàn)在09:00和20:00,低值出現(xiàn)在06:00和13:00,并表現(xiàn)出周末比平時(shí)高的“周末效應(yīng)”;月份和季度的變化為7月最低、1月最高,夏季最低、冬季最高,與能見(jiàn)度呈明顯負(fù)相關(guān)關(guān)系;年均值從2010年的79μg/m3到2015年的56μg/m3,總體呈下降趨勢(shì),并在全國(guó)范圍內(nèi)處于中間水平.相對(duì)濕度、日降水量和氣壓與PM10質(zhì)量濃度有更好的相關(guān)性,風(fēng)速越大對(duì)PM10質(zhì)量濃度影響越大.氣團(tuán)后向軌跡聚類分析和潛在源貢獻(xiàn)解析結(jié)果表明,臨安本底站在2010~2015年不同季節(jié)中,春季、秋季和冬季的氣團(tuán)大部分來(lái)自西北和北方,夏季更多的是來(lái)自東部海洋和沿海地帶,而不同年份里,每年的氣團(tuán)軌跡均以北方居多,其次是南方和東部海洋.總體分析得出, 2010~2015年臨安本底站的氣團(tuán)主要來(lái)自安徽、浙江、江蘇、廣東、福建和北方地區(qū),并且PM10高值主要來(lái)源于安徽、江蘇、上海和浙江等長(zhǎng)三角附近區(qū)域.

    區(qū)域本底站;PM10;氣象因素;后向軌跡;潛在源貢獻(xiàn)解析

    長(zhǎng)江三角洲是目前我國(guó)重要的大氣污染防治區(qū)之一,大氣環(huán)境污染已成為當(dāng)?shù)厝罕姺浅jP(guān)注的環(huán)境問(wèn)題[1-2].顆粒物是涉及范圍廣泛且具危害性的空氣污染物[2-4],可吸入顆粒物(PM10)可以被人體吸入危害健康[5],并對(duì)大氣能見(jiàn)度[6]和全球氣候變化[7]都有著重要影響,是影響我國(guó)環(huán)境空氣質(zhì)量的主要污染物[8].因此,臨安區(qū)域大氣本底站(以下簡(jiǎn)稱臨安本底站)作為反映長(zhǎng)三角地區(qū)大氣成分本底特征的唯一站點(diǎn),逐漸引起了相關(guān)學(xué)者們關(guān)注[9].學(xué)者們主要開展了臨安本底站氣溶膠的物理化學(xué)特征[10-11]、氣溶膠無(wú)機(jī)離子成分[12]、北方沙塵暴天氣影響下的氣溶膠演變特征[13]、細(xì)粒子的理化特性[14-16]等的研究.為了探討PM10質(zhì)量濃度長(zhǎng)時(shí)間段的變化趨勢(shì)及其影響因素,本研究通過(guò)對(duì)2010~2015年P(guān)M10質(zhì)量濃度的收集與整理,分析了長(zhǎng)三角臨安本底站PM10質(zhì)量濃度的日、周、月、季、年的變化規(guī)律,及其與大氣能見(jiàn)度等氣象要素的相關(guān)關(guān)系,并通過(guò)氣團(tuán)的后向軌跡聚類分析和PM10潛在源貢獻(xiàn)解析來(lái)研究氣團(tuán)輸送PM10的軌跡特征,可為長(zhǎng)三角地區(qū)大氣污染防治提供重要的科學(xué)依據(jù).

    1 材料與方法

    1.1 觀測(cè)站點(diǎn)介紹

    臨安本底站(119.44°E,30.18°N),海拔138. 6m,是世界氣象組織于1983年選址建成的全球大氣觀測(cè)網(wǎng)(WMO/GAW)的成員站,代表長(zhǎng)三角地區(qū)的本底狀況,2005年被選為大氣成分本底國(guó)家野外站.站址地屬于中亞熱帶季風(fēng)氣候區(qū),常年主導(dǎo)風(fēng)向?yàn)闁|北風(fēng)(如圖1,為風(fēng)向頻率和平均風(fēng)速的風(fēng)玫瑰圖,其中平均風(fēng)速擴(kuò)大了10倍),森林覆蓋良好,該站距杭州市約50km,距上海市約210km,距南京約250km處(如圖2),是長(zhǎng)三角地區(qū)的大氣本底站.

    1.2 數(shù)據(jù)采集來(lái)源和處理

    1.2.1 PM10質(zhì)量濃度 數(shù)據(jù)由美國(guó)Rupprecht & Patashick公司的TEOM1400a型大氣粒子監(jiān)測(cè)儀觀測(cè)得到,該儀器每10min采集1個(gè)數(shù)據(jù),測(cè)量范圍是0~5g/m3,分辨率為0.1μg/m3,精密度是± 1.5μg/m3(1小時(shí)濃度平均值)和±0.5μg/m3(24時(shí)濃度平均值).數(shù)據(jù)處理:先收集2010~2015年P(guān)M10的小時(shí)值,再去掉因儀器故障、進(jìn)氣管內(nèi)水滴、施工等產(chǎn)生的異常值,然后參考《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB3095-2012)[17]中PM10有效性規(guī)定,即每日至少20個(gè)小時(shí)平均濃度值或采樣時(shí)間,每年至少有分布均勻的324個(gè)日均值獲得有效數(shù)據(jù).

    1.2.2 能見(jiàn)度 數(shù)據(jù)由Vaisala公司的FD12儀對(duì)大氣光學(xué)能見(jiàn)度進(jìn)行觀測(cè), 采樣頻率為5min,探測(cè)范圍是10~50000m,在對(duì)觀測(cè)數(shù)據(jù)進(jìn)行質(zhì)量控制后平均獲得.本文所用地面氣象數(shù)據(jù)來(lái)自于臨安本底站自動(dòng)氣象觀測(cè)站(K1448號(hào))采集的逐時(shí)本站氣壓、溫度、風(fēng)向、風(fēng)速、相對(duì)濕度、日降水量等常規(guī)氣象數(shù)據(jù).

    1.2.3 后向軌跡 本文利用TrajStat軟件[18]來(lái)做氣團(tuán)的后向軌跡聚類分析,并結(jié)合PM10質(zhì)量濃度進(jìn)行潛在源貢獻(xiàn)解析(PSCF),氣象數(shù)據(jù)為GDAS數(shù)據(jù)(全球資料同化系統(tǒng)),下載網(wǎng)址為ftp://arlftp.arlhq. noaa.gov/pub/archives/gdas1/.

    2 結(jié)果與討論

    2.1 2010~2015年P(guān)M10質(zhì)量濃度變化特征

    2.1.1 PM10質(zhì)量濃度日均值 選取2010~ 2015年臨安本底站的PM10質(zhì)量濃度的小時(shí)資料(共計(jì)52515個(gè)有效小時(shí)值)進(jìn)行分析,按照《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》一共得到2192個(gè)有效日均值,其中有96.17%達(dá)到《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》的二級(jí)空氣質(zhì)量標(biāo)準(zhǔn)(24h平均濃度限值為150μg/m3).就超標(biāo)情況來(lái)看(表1),2010年超標(biāo)率最高,為8.49%,超標(biāo)倍數(shù)范圍為0.01~5.75, 2015年超標(biāo)率最低,是1.64%,超標(biāo)倍數(shù)范圍為0.01~0.24,6年日均值超標(biāo)率為3.83%.

    表1 2010~2015年臨安本底站PM10 質(zhì)量濃度日均值超標(biāo)率統(tǒng)計(jì)表

    如圖3所示,2010~2015年,臨安本底站PM10質(zhì)量濃度日均值成規(guī)律的波浪狀,其中最大日均值為1012μg/m3、473μg/m3和406μg/m3,分別出現(xiàn)在2010年3月21日、20日和22日,是由于在此期間,由蒙古氣旋強(qiáng)烈發(fā)展引發(fā)的一次強(qiáng)沙塵暴天氣,先后影響了21個(gè)省(區(qū)、市)[19];其次是在2011年4月28日至30日,在西北、華北、黃淮等地出現(xiàn)了2011年影響時(shí)間最長(zhǎng)、范圍最廣的一次沙塵天氣過(guò)程[20];另外在2013年12月1~9日,我國(guó)中東部地區(qū)集中發(fā)生了嚴(yán)重的灰霾污染過(guò)程,造成74個(gè)主要城市的重度污染天氣,其中長(zhǎng)三角地區(qū)污染最為嚴(yán)重[21].

    2.1.2 PM10質(zhì)量濃度的日變化、周變化及其與氣象要素的相關(guān)關(guān)系 大氣能見(jiàn)度與大氣顆粒物濃度密切相關(guān)[22-23],能很好地反映空氣質(zhì)量的變化趨勢(shì)[24-26].同時(shí),顆粒物濃度受天氣系統(tǒng)和氣象條件等多因素的共同影響[27],本文用2010~ 2015年臨安本底站的能見(jiàn)度(共計(jì)52077個(gè)有效小時(shí)值)、本站氣壓、空氣溫度和相對(duì)濕度的平均日變化與PM10質(zhì)量濃度日變化進(jìn)行對(duì)比分析,得到其關(guān)系,如圖4、圖5.

    研究發(fā)現(xiàn),因人類工作日和休息日中生產(chǎn)及生活方式的變化,將不同量的污染物排入大氣,所以全球很多站點(diǎn)的大氣氣溶膠質(zhì)量濃度、污染物質(zhì)量濃度、氣溫等天氣要素在周中和周末存在明顯的差異[30],這種周循環(huán)被稱為周末效應(yīng)[31].本文分析了2010~2015年P(guān)M10質(zhì)量濃度日均值的周變化及其與能見(jiàn)度的關(guān)系(如圖6),以及2010~ 2015年P(guān)M10質(zhì)量濃度、本站氣壓、空氣溫度和相對(duì)濕度的平均周變化及其關(guān)系(圖7).

    由圖6、圖7可以看出, 2010~2015年P(guān)M10質(zhì)量濃度和各氣象要素都有比較明顯的“周末效應(yīng)”.在平均周變化中:周六、周日PM10質(zhì)量濃度最高,能見(jiàn)度最差.導(dǎo)致這種現(xiàn)象的原因可能與周邊人類活動(dòng)和觀測(cè)站的位置有關(guān).該站遠(yuǎn)離工業(yè)區(qū),受工作日工業(yè)污染的影響相對(duì)較小,而現(xiàn)代人物資越來(lái)越豐富,周末活動(dòng)也越來(lái)越多,交通污染等使得能見(jiàn)度降低,空氣污染加劇.

    2.1.3 PM10質(zhì)量濃度月變化和季變化 將2010~2015年P(guān)M10質(zhì)量濃度和能見(jiàn)度的日均值按每年的月份和季度(3~5月為春季,6~8月為夏季,9~11月為秋季,12月、1月和2月為冬季)統(tǒng)計(jì),得出每年的月均值和季均值.由圖8可以看出,PM10質(zhì)量濃度和大氣能見(jiàn)度均有顯著的月變化和季變化規(guī)律:(1)PM10質(zhì)量濃度夏季最低,平均值為48μg/m3,其中每年月均值變化范圍為32~70μg/m3,平均月均值最低為7月44μg/m3;春季、秋季和冬季每年月均值變化范圍比較大,分別為52~111μg/m3、50~119μg/m3和52~ 124μg/m3,春季和秋季的平均值分別為76μg/m3和71μg/m3,冬季則最高,為81μg/m3,是夏季的1.7倍左右,其中平均月均值最高月份為1月93μg/m3.(2)大氣能見(jiàn)度剛好相反,夏季的能見(jiàn)度最好,平均值為8663m,其中全年7月能見(jiàn)度最好為11038m,其次是春季和秋季,能見(jiàn)度為7179m和6203m,而冬季最差,平均值為5267m,只有夏季能見(jiàn)度的61%,其中1、2月只有4770m和4638m.這是因?yàn)槎灸鏈貙虞^強(qiáng)較厚,空氣流動(dòng)慢,不利于污染物的擴(kuò)散,而在夏季,大氣垂直擴(kuò)散強(qiáng),溫度高,降水充沛,對(duì)大氣顆粒物有較好的清除作用,故而PM10質(zhì)量濃度最低,能見(jiàn)度也最好.

    由圖8(a)可知,大部分年份PM10質(zhì)量濃度在2月會(huì)有一個(gè)降低,3月又升高的過(guò)程,這是由于每年2月起降水量開始增加,風(fēng)速開始增大,但是由于3月左右我國(guó)開始出現(xiàn)沙塵天氣現(xiàn)象(2010~2015年首次沙塵天氣現(xiàn)象出現(xiàn)時(shí)間分別為:3月8日、3月22日、3月20日、2月24日、3月19日和2月21日),所以PM10質(zhì)量濃度又有增加.

    2.1.4 PM10質(zhì)量濃度年變化 2010~2015年,臨安本底站PM10質(zhì)量濃度年均值總體呈下降趨勢(shì)(圖9),6年的平均年均值為69μg/m3,達(dá)到環(huán)境空氣質(zhì)量二級(jí)標(biāo)準(zhǔn)(70μg/m3).PM10質(zhì)量濃度由2010年的79μg/m3到2015年的56μg/m3,降幅為28.62%.并且從2012年起, PM10質(zhì)量濃度年均值基本達(dá)標(biāo),分別為64μg/m3、71μg/m3、70μg/m3和56μg/m3, 與2002~2012年全國(guó)地級(jí)及以上城市PM10年均濃度的平均值呈下降趨勢(shì)一致[3].而2010~2015年,大氣能見(jiàn)度年均值呈波動(dòng)變化,與PM10質(zhì)量濃度的變化沒(méi)有明顯的同步性.

    根據(jù)國(guó)家環(huán)保部公布的《中國(guó)環(huán)境公報(bào)》[32]結(jié)果顯示,在全國(guó)范圍來(lái)看,如圖10所示, 長(zhǎng)三角主要城市2013~2015年P(guān)M10年均質(zhì)量濃度(103μg/m3、92μg/m3和83μg/m3)在京津冀、長(zhǎng)三角和珠三角區(qū)三大重點(diǎn)區(qū)域中處于中間水平,并且比全國(guó)74個(gè)主要城市的平均PM10質(zhì)量濃度年均值(118μg/m3、105μg/m3和87μg/m3)略低,而臨安本底站的年均值比長(zhǎng)三角區(qū)低,略高于珠三角區(qū).

    2.2 PM10質(zhì)量濃度與主要地面氣象要素的相關(guān)性分析

    在污染源一定的條件下,污染物濃度的大小主要取決于氣象條件[33-34],氣象因素是影響PM10質(zhì)量濃度的主要因素[35].本文對(duì)2010~2015年P(guān)M10質(zhì)量濃度與主要?dú)庀笠?能見(jiàn)度、相對(duì)濕度、本站氣壓、空氣溫度、日降水量等)進(jìn)行了相關(guān)性分析,如表2所示.

    由表2可知,PM10質(zhì)量濃度與能見(jiàn)度、相對(duì)濕度、本站氣壓、空氣溫度、降水量、極大風(fēng)速(小時(shí)內(nèi)3秒鐘滑動(dòng)平均風(fēng)速的最大值)和十分鐘風(fēng)速(十分鐘滑動(dòng)平均風(fēng)速)均呈顯著相關(guān).其中,相對(duì)濕度和日降水量與PM10質(zhì)量濃度呈較好的負(fù)相關(guān),相關(guān)系數(shù)分別為-0.367和-0.254,這主要是由于降水對(duì)大氣污染物的沖刷清潔作用.其次是本站氣壓與PM10質(zhì)量濃度呈較好的正相關(guān),相關(guān)系數(shù)為0.249,因?yàn)榈蛪簳r(shí),氣流上升,近地面的污染物會(huì)隨著上升氣流擴(kuò)散到高空,繼而擴(kuò)散到范圍大的其他區(qū)域,污染減輕,而高壓時(shí),氣流下沉,污染物不易擴(kuò)散.同時(shí),PM10質(zhì)量濃度與極大風(fēng)速和十分鐘風(fēng)速的相關(guān)系數(shù)分別為-0.119和-0.087,這是由于大氣湍流對(duì)空氣污染物的擴(kuò)散作用,而極大風(fēng)速的相關(guān)性比十分鐘風(fēng)速的更為顯著,說(shuō)明了大風(fēng)對(duì)PM10質(zhì)量濃度的影響更大.

    表2 2010~2015年臨安本底站PM10質(zhì)量濃度與各氣象要素的Pearson相關(guān)系數(shù)(n=2191)

    注:**表示在0.01水平上顯著相關(guān).

    2.3 氣團(tuán)后向軌跡聚類分析和潛在源貢獻(xiàn)解析

    2.3.1 氣團(tuán)后向軌跡聚類分析 氣團(tuán)后向軌跡分析可以了解氣團(tuán)輸送顆粒物的路線軌跡和輸送速度,為了探討長(zhǎng)三角臨安本底站氣團(tuán)軌跡的輸送特征,本文對(duì)長(zhǎng)三角臨安本底站2010~2015年每天當(dāng)?shù)貢r(shí)間08時(shí)的72h氣團(tuán)后向軌跡(距地500m)進(jìn)行聚類分析,根據(jù)氣團(tuán)軌跡的途經(jīng)區(qū)域、移動(dòng)速度和方向?qū)⑵浞譃?類,得到2010~2015年每個(gè)季節(jié)、每年和6年總體的后向軌跡聚類結(jié)果圖.

    如圖11所示,2010~2015年4個(gè)季節(jié)中,從輸送方向上看,春季和秋季的氣團(tuán)聚類軌跡大部分來(lái)自西北和北方向,少部分來(lái)自海洋,春季還有部分來(lái)自西南方向,夏季更多的是來(lái)自東部海洋和沿海地帶,而冬季均來(lái)自西北和北方;從出現(xiàn)概率來(lái)看,春、秋、冬季的氣團(tuán)軌跡出現(xiàn)概率最高的均來(lái)自北方向,概率分別為36.5%、43.7%和35.3%,而夏季出現(xiàn)概率最高的來(lái)自東北方向,為34.7%;其次,出現(xiàn)概率第二高的,春季是來(lái)自西南方(33.1%),夏季是來(lái)自南方(34.1%),秋季是來(lái)自東北方(26.0%),冬季是來(lái)自西北方(34.0%);從輸送速度上看, 4個(gè)季節(jié)的輸送距離和速度變化都比較大,尤其是都存在出現(xiàn)概率比較高的慢速短距離輸送.所以從4個(gè)季節(jié)氣團(tuán)軌跡聚類分析看出, 2010~2015年臨安本底站春季大氣顆粒物的外來(lái)輸送可能主要受江蘇、浙北、江西、山東、河北和內(nèi)蒙古的影響;夏季可能主要受到浙江、上海、蘇南、江西和廣東的影響;而秋季和冬季可能主要受浙北、江蘇、安徽、上海、山東、河北和內(nèi)蒙古的影響.

    由圖12可知,從每年的軌跡規(guī)律來(lái)看,總體上,2010~2015年臨安本底站每年的氣團(tuán)軌跡聚類分析結(jié)果比較相似,以北方居多,南方和東邊海洋次之.每年聚類軌跡中出現(xiàn)概率最高的(2011年55.8%、2012年53.6%、2013年51.1%、2015年42.6%)來(lái)自北方向,東北方向(2010年47.5%)和南方(2014年35.1%),且都是短距離慢速輸送氣團(tuán);而出現(xiàn)概率第二高的除了2013年(23.2%)的來(lái)自西南方,其余年度的(2010年27.4%、2011年20.1%、2012年18.5%、2014年30.3%、2015年25.1%)均來(lái)自北方或東北方,同時(shí)每年都有部分來(lái)自西北方或北方的遠(yuǎn)距離快速輸送氣團(tuán).所以,2010~2015年臨安本底站每年的氣團(tuán)輸送可能受江蘇、浙江、山東、河北和內(nèi)蒙古的影響較大,部分年份福建、江西和廣東也有一定影響.

    最后,由2010~2015年臨安本底站1961條氣團(tuán)后向軌跡線聚合成4條軌跡線,如圖13所示.總體來(lái)看,出現(xiàn)概率最高的氣團(tuán)軌跡是來(lái)自安徽的短距離慢速輸送,出現(xiàn)概率是39.0%;其次是來(lái)自東部的海洋、浙北和南部的廣東、福建和浙江西部,出現(xiàn)概率分別是30.1%和22.0%;然后是來(lái)自北部的江蘇、山東、北京、河北和內(nèi)蒙古的遠(yuǎn)距離快速輸送,出現(xiàn)概率為8.8%.

    2.3.2 潛在源貢獻(xiàn)解析 為了更好地進(jìn)行PM10的來(lái)源分析,本文采用潛在源貢獻(xiàn)的解析.

    潛在源貢獻(xiàn)解析,Potential Source Contribution Function(PSCF)是基于條件概率函數(shù)[36],通過(guò)結(jié)合氣團(tuán)軌跡和要素值來(lái)識(shí)別大氣污染物潛在來(lái)源可能方位的一種方法[37-38]. PSCF函數(shù)的定義是經(jīng)過(guò)某一區(qū)域的氣團(tuán)到達(dá)點(diǎn)所對(duì)應(yīng)某要素值(本文指PM10質(zhì)量濃度)超過(guò)設(shè)定閾值(本研究設(shè)定值為70μg/m3,總體樣本的中值)的條件概率. PSCF的算法是計(jì)算網(wǎng)格范圍內(nèi)每條軌跡的要素值,如果要素值高于閾值,則此格點(diǎn)對(duì)于受體要素高值貢獻(xiàn)的可能性較大[39],屬于這個(gè)格點(diǎn)的排放物將在氣團(tuán)的攜帶下到達(dá)最終的受體,所以通過(guò)PSCF可以給出可能造成受體地點(diǎn)某一要素高值的可能源地. 但是由于網(wǎng)格與采樣點(diǎn)的距離增加會(huì)使得PSCF的誤差增加[40],所以為了降低PSCF的不確定性本文引入Wij(權(quán)重因子)[41].將研究區(qū)域劃分為0.5°*0.5°的格點(diǎn),落在某格點(diǎn)的軌跡條數(shù)為nij,格點(diǎn)軌跡對(duì)應(yīng)要素值大于閾值的數(shù)量為mij, PSCF和Wij計(jì)算如公式(1)和(2)所示:

    (n)=(2)

    結(jié)合2010~2015年的氣團(tuán)軌跡和PM10質(zhì)量濃度, 用PSCF分析對(duì)臨安本底站PM10質(zhì)量濃度起到貢獻(xiàn)的可能來(lái)源.PSCF分析結(jié)果如圖14.

    從圖14可以看出,2010~2015年臨安本底站PM10可能主要來(lái)自于江蘇、安徽、上海和浙江,尤其是集中在經(jīng)濟(jì)發(fā)展較好、排放較多的長(zhǎng)三角地區(qū).而黃海和東海臨近海域?qū)﹂L(zhǎng)三角臨安本底站PM10質(zhì)量濃度也有一定貢獻(xiàn),一方面可能是受海洋氣溶膠的影響,另外可能是因?yàn)楹C嫔蠞穸容^大,氣溶膠吸濕增加,加劇了污染氣團(tuán)中顆粒物的二次生成.

    3 結(jié)論

    3.1 臨安本底站2010~2015年P(guān)M10質(zhì)量濃度日均值的達(dá)標(biāo)率是96.17%,從日變化上看,PM10質(zhì)量濃度小時(shí)變化呈“雙峰雙谷”分布,峰值出現(xiàn)在09:00和20:00,低值出現(xiàn)在6:00和13:00;從周變化上看,PM10質(zhì)量濃度呈現(xiàn)出平時(shí)低周末高的“周末效應(yīng)”;從月、季變化上看,PM10質(zhì)量濃度7月最低1月最高,夏季最低冬季最高,與能見(jiàn)度呈明顯負(fù)相關(guān)關(guān)系;從年變化上看,從2010~2015年, 總體呈下降趨勢(shì),降幅為28.62%,2013~2015年年均值高于珠三角地區(qū)而較其他地區(qū)低.

    3.2 臨安本底站PM10質(zhì)量濃度與相對(duì)濕度、日降水量呈現(xiàn)明顯的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)分別為-0.367和-0.254;與氣壓呈現(xiàn)明顯的正相關(guān)關(guān)系,相關(guān)系數(shù)為0.249;而與極大風(fēng)速和十分鐘風(fēng)速的相關(guān)系數(shù)分別為-0.119和-0.087,表明風(fēng)速越大對(duì)PM10質(zhì)量濃度的影響越大.

    3.3 氣團(tuán)后向軌跡聚類分析和潛在源貢獻(xiàn)解析結(jié)果表明:臨安本底站2010~2015年春季、秋季和冬季的氣團(tuán)大部分來(lái)自西北和北方向,少部分來(lái)自海洋和西南方向,夏季更多的是來(lái)自東部海洋和沿海地帶; 2010~2015年臨安本底站每年的氣團(tuán)軌跡以北方居多,其次是南方和東邊海洋,每年聚類軌跡中出現(xiàn)概率最高的(2010年47.5%、2011年55.8%、2012年53.6%、2013年51.1%、2015年42.6%)大部分都來(lái)自北方或東北方;總體來(lái)看2010~2015年,臨安本底站出現(xiàn)概率最高的氣團(tuán)軌跡是來(lái)自安徽(39.0%),其次是來(lái)自浙北、東部的海洋(30.1%),和南部的廣東、福建、浙西(22.0%),以及來(lái)自江蘇等北方向的遠(yuǎn)距離快速輸送(8.8%).2010~2015年臨安本底站可吸入顆粒物可能主要來(lái)自江蘇、安徽、上海和浙江,以及黃海和東海臨近海域.

    [1] 李 凱,劉 濤,曹廣忠.城市群空間集聚和擴(kuò)散的特征與機(jī)制—以長(zhǎng)三角城市群、武漢城市群和成渝城市群為例 [J].城市規(guī)劃, 2016,40(2):18-26.

    [2] 金 均,吳 建,蔡菊珍,等.杭州市灰霾天氣基本特征及成因分析[J]. 環(huán)境污染與防治, 2010,32(5):61-67.

    [3] 李名升,張建輝,張殷俊,等.近10年中國(guó)大氣PM10污染時(shí)空格局演變[J]. 地理學(xué)報(bào), 2013,68(11):1504-1512.

    [4] Hu M G, Jia L, Wang J F et al. Spatial and temporal characteristics of particulate matter in Beijing, China using the empirical mode decomposition method [J]. Science of the Total Environment, 2013,(458-460):70-80.

    [5] Qiu H, Yu I T, Tian L W et al. Effects of coarse particulate matter on emergency hospital admissions for respiratory diseases: A time-series analysis in Hong Kong [J]. Environmental Health Perspectives, 2012,120(4):572-576.

    [6] Noll K E, Mueller P K, Inada M. Visibility and aerosol concentration in urban air [J]. Atmospheric Environment, 1968, 2(5):465-475.

    [7] Buseck P R, Posfai M. Airborne minerals and related aerosol particles: Effects on climate and the environment [J]. Proceedings of the National Academy of Sciences of the United States of America, 1999,96(7):3372-3379.

    [8] 中華人民共和國(guó)環(huán)境保護(hù)部.中國(guó)環(huán)境質(zhì)量報(bào)告2011 [R]. 北京:中國(guó)環(huán)境科學(xué)出版社, 2012.

    [9] 丁國(guó)安,鄭向東,馬建中,等.2006.等.近30年大氣化學(xué)和大氣環(huán)境研究回顧—紀(jì)念中國(guó)氣象科學(xué)研究院成立50周年[J].應(yīng)用氣象學(xué)報(bào), 17(6):796-814.

    [10] 楊東貞,于曉嵐,李興生,等.臨安大氣污染本底站氣溶膠特征分析[J]. 大氣科學(xué), 1995,19(2):219-228.

    [11] 張養(yǎng)梅,顏 鵬,楊東貞,等.臨安大氣氣溶膠理化特性季節(jié)變化[J]. 應(yīng)用氣象學(xué)報(bào), 2007,18(5):635-644.

    [12] 楊東貞,顏 鵬,張養(yǎng)梅,等.WMO區(qū)域本底站氣溶膠特征分析[J]. 第四紀(jì)研究, 2006,26(5):733-741.

    [13] 顏 鵬,毛節(jié)泰,楊東貞,等.臨安一次沙塵暴過(guò)程影響氣溶膠物理化學(xué)特性演變的初步分析[J]. 第四紀(jì)研究, 2004,24(7): 437-446.

    [14] 孟昭陽(yáng),賈小芳,張仁建,等.2010年長(zhǎng)三角洲臨安本底站PM2.5理化特征[J].應(yīng)用氣象學(xué)報(bào), 2012,23(4):424-432.

    [15] 徐宏輝,劉 潔,王躍思,等.臨安本底站大氣氣溶膠水溶性離子濃度變化特征[J]. 環(huán)境化學(xué), 2012,31(6):796-802.

    [16] 張軼雯,張小曳,張養(yǎng)梅,等.長(zhǎng)三角臨安本底PM1中各主要化學(xué)成分質(zhì)量-粒度分布與變化[J]. 氣候變化研究進(jìn)展, 2015, 11(4):270-280.

    [17] GB3095-2012 環(huán)境空氣質(zhì)量標(biāo)準(zhǔn) [S].

    [18] Wang Y Q, Zhang X Y, Draxler 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.

    [19] 王汝佛,馮 強(qiáng),尚 可.2014.2010春季我國(guó)一次強(qiáng)沙塵暴過(guò)程分析[J]. 干旱區(qū)地理, 37(1):31-44.

    [20] 中華人民共和國(guó)環(huán)境保護(hù)部.中國(guó)環(huán)境公報(bào), 2011. http://www.zhb.gov.cn/hjzl/zghjzkgb/lnzghjzkgb/.

    [21] 劉 璇,朱 彬,侯雪偉,等.基于衛(wèi)星遙感和地面觀測(cè)資料的霾過(guò)程分析——以長(zhǎng)江三角洲地區(qū)一次嚴(yán)重霾污染過(guò)程為例[J]. 中國(guó)環(huán)境科學(xué), 2015,35(3):641-651.

    [22] 馬雁軍,左洪超,張?jiān)坪?等.遼寧中部城市群大氣能見(jiàn)度變化趨勢(shì)及影響因子分析[J]. 高原氣象, 2005,24(4):623-628.

    [23] Charlson R J. Atmospheric visibility related to aerosol mass concentration: A review [J]. Environmental Science and Technology, 1969,3(10):430- 435.

    [24] 吳 兌,畢雪巖,鄧雪嬌,等.珠江三角洲大氣灰霾導(dǎo)致能見(jiàn)度下降問(wèn)題研究[J]. 氣象學(xué)報(bào), 2006,64(4):510-517.

    [25] Noll K E, Mueller P K, Imada M. Visibility and aerosol concentration in urban air [J]. Atmospheric Environment (1967), 1968,2(5):465-475.

    [26] Charlson R J. Atmospheric visibility related to aerosol mass concentration: review [J]. Environmental Science & Technology, 1969,3(10):913-918.

    [27] 劉 寧,王雪松,胡泳濤,等.珠江三角洲秋季PM10污染模擬與形成過(guò)程分析[J]. 中國(guó)環(huán)境科學(xué), 2012,32(9):1537-1545.

    [28] 孫景群.能見(jiàn)度與相對(duì)濕度的關(guān)系 [J]. 氣象學(xué)報(bào), 1985,43(2): 230-234.

    [29] 周奕珂,朱 彬,韓志偉,等.長(zhǎng)江三角洲地區(qū)冬季能見(jiàn)度特征及影響因子分析[J]. 中國(guó)環(huán)境科學(xué), 2016,36(3):660-669.

    [30] Forster P M, Solomons. Observations of a weekend effect in diurnal temperature range [J]. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(2):11225-11230.

    [31] 龔道溢,郭 棟,羅 勇.中國(guó)夏季日降水頻次的周末效應(yīng)[J]. 氣候變化研究進(jìn)展, 2006,2(3):131-134.

    [32] 中華人民共和國(guó)環(huán)境保護(hù)部.中國(guó)環(huán)境公報(bào)2013~2015 [EB/OL] //http://www.zhb.gov.cn/hjzl/zghjzkgb/lnzghjzkgb/.

    [33] 蔣維楣,曹文俊,蔣瑞賓.空氣污染氣象學(xué)教程[M]. 北京:氣象出版社, 1993.

    [34] 吳 兌,鄧雪嬌.環(huán)境氣象學(xué)與特種氣象預(yù)報(bào)[M]. 北京:氣象出版社, 2001.

    [35] 邱玉珺,牛生杰,沈建國(guó).內(nèi)蒙古沙地和草地環(huán)境下近地面PM10分布特征及其與相關(guān)氣象要素的關(guān)系[J].氣候與環(huán)境研究, 2009,14(1):97-104.

    [36] Begum B A, Kim E, Jeong C H, et al. Evaluation of the potential source contribution function using the 2002 Quebec forest fire episode [J]. Atmospheric Environment, 2005,39(20):3719-3724.

    [37] Gao N, Cheng M D, Hopke P K. Potential source contribution function analysis and source apportionment of sulfur species measured at Rubidoux, CA during the Southern California Air Quality Study, 1987 [J]. Analytica chimica acta, 1993,277(2): 369-380.

    [38] Gao N, Hopke P, Reid N. Possible sources for some trace elements found in airborne particles and precipitation in Dorset, Ontario [J]. Journal of the air & waste management association, 1996,46(11):1035-1047.

    [39] Zeng Y, Hopke P. A study of the sources of acid precipitation in Ontario, Canada [J]. Atmospheric Environment (1967), 1989,23(7):1499-1509.

    [40] 王愛(ài)平,朱 彬,銀 燕,等.黃山頂夏季氣溶膠數(shù)濃度特征及其輸送潛在源區(qū)[J]. 中國(guó)環(huán)境科學(xué), 2014,34(4):852-861.

    [41] Polissar A V, Hopke P K, Kaufmann P P, et al. The aerosol at Barrow, Alaska: long-term trends and source location [J]. Atmospheric Environment, 1999,33:2441-2458.

    Variation characteristics of PM10and its influencing factors at Lin'an regional background station in Yangtze River Delta region during 2010~2015.

    YUE Yi1,2, LI Jin-juan1*, MA Qian-li3

    (1.Resource and Environmental Engineering College, Guizhou University, Guiyang 550025, China;2.Lin'an Meteorological Bureau, Lin'an 311300, China;3.Lin'an Regional Atmosphere Background Station, Lin'an 311307, China)., 2017,37(8):2877~2887

    PM10mass concentration variations in the Yangtze River Delta Region and the causes were investigated, using measurement of PM10mass concentration and meteorology during 2010~2015 at the regional representative Lin’an Regional Atmosphere Background Station. The results show that PM10mass concentration met the Grade 2 Standard specified in the Ambient Air Quality Standards (GB3096-2012) for 96.17% during 2010~2015. Diurnal PM10mass concentration showed the maximum value at 09:00 and 20:00, minimum at 06:00 and 13:00, and it also presented a “weekend effect” with higher value on weekend than on weekdays. Monthly variation showed the lowest PM10mass concentration in July in summer and the highest in January, which negatively correlated with visibility. Annual average of PM10mass concentration at Lin’an decreased from 79μg/m3to 56μg/m3from 2010 to 2015, which is a moderate concentration level in China. PM10mass concentration correlated well with relative humidity, daily precipitation and local air pressure. Stronger wind also had greater impact on PM10mass concentration. Back-trajectories and cluster analysis reveal that prevailing air masses were from the northwest and the north in spring, autumn, and winter, and from the sea and the coastal areas in the east in summer. Throughout the years, the dominant air masses were from the north, followed by southerly and east marine air masses. It is concluded that the majority air mass at Lin’an arrived from Anhui, Zhejiang, Jiangsu, Guangdong, Fujian and the north region. Potential source contribution function (PSCF) results suggest that PM10pollution may mainly originate from the areas around the Yangtze River Delta region, including Shanghai, Anhui, Jiangsu and Zhejiang Provinces.

    regional background station;PM10;meteorological factors;back-trajectory;potential source contribution function (PSCF)

    X513

    A

    1000-6923(2017)08-2877-11

    岳 毅(1982-),女,布依族,貴州貴陽(yáng)人,貴州大學(xué)碩士研究生,研究方向?yàn)榇髿饪茖W(xué).

    2017-01-24

    國(guó)家自然科學(xué)基金項(xiàng)目(41265008);公益性行業(yè)(氣象)科研專項(xiàng)子項(xiàng)目(GYHY201206037-06);貴州省重點(diǎn)學(xué)科建設(shè)項(xiàng)目(黔學(xué)位合字ZDXK[2016]11號(hào))

    * 責(zé)任作者, 教授, summy_lee@163.com

    猜你喜歡
    氣團(tuán)臨安能見(jiàn)度
    輸水管線啟動(dòng)填充過(guò)程含滯留氣團(tuán)瞬變流數(shù)值模擬
    基于后向軌跡模式的合肥市臭氧來(lái)源分析
    2005—2017年白云機(jī)場(chǎng)能見(jiàn)度變化特征及其與影響因子關(guān)系研究
    氣團(tuán)來(lái)源對(duì)沿海城市PM2.5中二次水溶性無(wú)機(jī)離子形成特征的影響研究
    HZ330CL橫折缺陷原因分析與改進(jìn)
    我校臨安校區(qū)簡(jiǎn)介
    低能見(jiàn)度下高速公路主動(dòng)誘導(dǎo)技術(shù)的應(yīng)用
    前向散射能見(jiàn)度儀的常見(jiàn)異?,F(xiàn)象處理及日常維護(hù)
    前向散射能見(jiàn)度儀故障實(shí)例分析
    浙江臨安郎碧村宋墓 一座平凡墓葬的啟示
    大眾考古(2015年10期)2015-06-26 08:00:04
    日本五十路高清| 久久久久精品国产欧美久久久| 女人精品久久久久毛片| 亚洲精品国产一区二区精华液| 色综合欧美亚洲国产小说| 中出人妻视频一区二区| 国产激情久久老熟女| 怎么达到女性高潮| 王馨瑶露胸无遮挡在线观看| 下体分泌物呈黄色| 老熟妇乱子伦视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 日本五十路高清| 国产高清激情床上av| 天天影视国产精品| www.自偷自拍.com| 国产精品国产高清国产av | 叶爱在线成人免费视频播放| 交换朋友夫妻互换小说| 两性夫妻黄色片| aaaaa片日本免费| 操美女的视频在线观看| 国产欧美日韩综合在线一区二区| 亚洲色图 男人天堂 中文字幕| 国产精品美女特级片免费视频播放器 | 大片电影免费在线观看免费| 久久久国产一区二区| 一级毛片女人18水好多| 欧美激情久久久久久爽电影 | 国产成人一区二区三区免费视频网站| 成年女人毛片免费观看观看9 | 午夜视频精品福利| 午夜福利一区二区在线看| av线在线观看网站| 99re在线观看精品视频| 久久久久久久午夜电影 | 久久人妻福利社区极品人妻图片| 久久99一区二区三区| 成人手机av| bbb黄色大片| 欧美黄色淫秽网站| 夫妻午夜视频| 欧美丝袜亚洲另类 | tocl精华| 无人区码免费观看不卡| 国产av又大| 久久这里只有精品19| 老司机福利观看| 夜夜躁狠狠躁天天躁| 91九色精品人成在线观看| 欧美性长视频在线观看| 中文字幕高清在线视频| 9色porny在线观看| 丰满迷人的少妇在线观看| 美女福利国产在线| 亚洲欧美激情在线| 极品少妇高潮喷水抽搐| 国产一卡二卡三卡精品| 黄频高清免费视频| 亚洲黑人精品在线| 久久香蕉激情| 最新美女视频免费是黄的| 一级毛片高清免费大全| 欧美日韩av久久| 亚洲欧美激情在线| 国产一区在线观看成人免费| 一个人免费在线观看的高清视频| 国产精品偷伦视频观看了| 久久这里只有精品19| 久久久久久久久免费视频了| 成年版毛片免费区| 国产高清视频在线播放一区| 成年动漫av网址| 午夜福利免费观看在线| 大型av网站在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲一区二区三区欧美精品| bbb黄色大片| 亚洲免费av在线视频| 啦啦啦在线免费观看视频4| 9热在线视频观看99| 天天躁夜夜躁狠狠躁躁| 国产精品国产av在线观看| 狠狠狠狠99中文字幕| 精品国产一区二区久久| 夜夜夜夜夜久久久久| 91麻豆av在线| 男男h啪啪无遮挡| 成人特级黄色片久久久久久久| 精品亚洲成国产av| 日韩欧美三级三区| 高清黄色对白视频在线免费看| 女警被强在线播放| 久久久久久久久免费视频了| 老司机深夜福利视频在线观看| 国产在线一区二区三区精| 在线天堂中文资源库| 亚洲一码二码三码区别大吗| 精品一区二区三区四区五区乱码| 制服人妻中文乱码| 美女福利国产在线| 亚洲精品国产区一区二| 国产在线观看jvid| 国产91精品成人一区二区三区| 老熟妇仑乱视频hdxx| 麻豆成人av在线观看| 国产精品99久久99久久久不卡| 99久久人妻综合| 欧美黑人精品巨大| 村上凉子中文字幕在线| 国产一卡二卡三卡精品| 欧美精品人与动牲交sv欧美| 18禁裸乳无遮挡动漫免费视频| 日日爽夜夜爽网站| 97人妻天天添夜夜摸| 亚洲中文字幕日韩| 老司机靠b影院| 1024香蕉在线观看| 99国产精品免费福利视频| 露出奶头的视频| 18在线观看网站| 亚洲一区高清亚洲精品| 啦啦啦在线免费观看视频4| 热99国产精品久久久久久7| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看 | 女警被强在线播放| 国产精品一区二区在线不卡| 成年人午夜在线观看视频| 国产高清视频在线播放一区| av在线播放免费不卡| 久久国产精品大桥未久av| 久久亚洲真实| 首页视频小说图片口味搜索| 日日摸夜夜添夜夜添小说| 一级黄色大片毛片| 国产欧美日韩一区二区三区在线| 日韩成人在线观看一区二区三区| 亚洲成人国产一区在线观看| 精品免费久久久久久久清纯 | 久久婷婷成人综合色麻豆| 搡老岳熟女国产| 免费在线观看完整版高清| 久久国产精品大桥未久av| 中国美女看黄片| 色婷婷av一区二区三区视频| 性少妇av在线| 脱女人内裤的视频| 色老头精品视频在线观看| 日本黄色视频三级网站网址 | 国产一区二区三区在线臀色熟女 | 在线看a的网站| 嫁个100分男人电影在线观看| 99久久综合精品五月天人人| 一级毛片女人18水好多| 亚洲午夜理论影院| 亚洲熟女毛片儿| 国产精品国产av在线观看| 午夜亚洲福利在线播放| 99热网站在线观看| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠躁躁| 高清视频免费观看一区二区| 亚洲精华国产精华精| 国产亚洲欧美98| 啦啦啦在线免费观看视频4| 国产亚洲一区二区精品| 欧美在线黄色| 日韩视频一区二区在线观看| 一级作爱视频免费观看| 国产精品免费大片| 99久久99久久久精品蜜桃| 美女国产高潮福利片在线看| 丝袜在线中文字幕| 亚洲欧美色中文字幕在线| 免费在线观看完整版高清| 国产精品一区二区在线观看99| 人成视频在线观看免费观看| 国产精华一区二区三区| 老司机午夜十八禁免费视频| 高清在线国产一区| 亚洲va日本ⅴa欧美va伊人久久| 丝袜在线中文字幕| 热99久久久久精品小说推荐| 中文字幕另类日韩欧美亚洲嫩草| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷丁香在线五月| 人妻丰满熟妇av一区二区三区 | 亚洲色图av天堂| 中文字幕制服av| 99国产精品免费福利视频| 亚洲五月婷婷丁香| 一级作爱视频免费观看| 黄色a级毛片大全视频| 国产亚洲av高清不卡| 精品一区二区三卡| 精品国产一区二区三区久久久樱花| 午夜福利在线观看吧| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文日韩欧美视频| 不卡一级毛片| 黄色女人牲交| 亚洲aⅴ乱码一区二区在线播放 | 国产午夜精品久久久久久| 亚洲aⅴ乱码一区二区在线播放 | 国产xxxxx性猛交| 成人亚洲精品一区在线观看| 香蕉国产在线看| 久久久久视频综合| 咕卡用的链子| 乱人伦中国视频| 国产欧美日韩精品亚洲av| 久久性视频一级片| 久久影院123| 9热在线视频观看99| 涩涩av久久男人的天堂| 欧美午夜高清在线| 桃红色精品国产亚洲av| 正在播放国产对白刺激| 777久久人妻少妇嫩草av网站| 91成人精品电影| 免费在线观看完整版高清| 久久亚洲真实| 免费av中文字幕在线| 色综合婷婷激情| 男人舔女人的私密视频| 国产一卡二卡三卡精品| 99国产综合亚洲精品| 首页视频小说图片口味搜索| 精品高清国产在线一区| 久久人人97超碰香蕉20202| 国产一区二区激情短视频| 中文字幕人妻熟女乱码| 亚洲成人手机| 久久久国产欧美日韩av| 中文字幕人妻丝袜制服| 久久精品国产99精品国产亚洲性色 | 99热网站在线观看| 夜夜躁狠狠躁天天躁| 建设人人有责人人尽责人人享有的| 一进一出好大好爽视频| 99riav亚洲国产免费| 久久青草综合色| 一a级毛片在线观看| 99久久国产精品久久久| 身体一侧抽搐| av网站免费在线观看视频| 中文字幕高清在线视频| 亚洲成人国产一区在线观看| 午夜久久久在线观看| 99国产极品粉嫩在线观看| 国产精华一区二区三区| 国产乱人伦免费视频| 国产精品一区二区在线观看99| 最新美女视频免费是黄的| 一级a爱片免费观看的视频| 国产欧美日韩综合在线一区二区| 热re99久久精品国产66热6| 在线看a的网站| 99国产精品99久久久久| 国产一区二区三区综合在线观看| av天堂久久9| 巨乳人妻的诱惑在线观看| 久99久视频精品免费| 麻豆乱淫一区二区| 狠狠狠狠99中文字幕| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 国产高清国产精品国产三级| 久久久精品免费免费高清| 亚洲人成伊人成综合网2020| 久久婷婷成人综合色麻豆| av超薄肉色丝袜交足视频| 韩国av一区二区三区四区| 91麻豆精品激情在线观看国产 | 亚洲avbb在线观看| 美女视频免费永久观看网站| 国产免费av片在线观看野外av| 丝袜人妻中文字幕| 女人被狂操c到高潮| 天天躁日日躁夜夜躁夜夜| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品一区在线观看| 两人在一起打扑克的视频| 午夜激情av网站| 男女床上黄色一级片免费看| 精品久久久久久电影网| 国产精品国产高清国产av | 女人精品久久久久毛片| 9热在线视频观看99| 国产一卡二卡三卡精品| 精品熟女少妇八av免费久了| 在线国产一区二区在线| 亚洲国产毛片av蜜桃av| a级毛片在线看网站| 久久国产精品人妻蜜桃| 老汉色∧v一级毛片| 大陆偷拍与自拍| 高清av免费在线| 国产亚洲精品久久久久久毛片 | 高潮久久久久久久久久久不卡| 亚洲欧美日韩另类电影网站| 欧美黑人欧美精品刺激| 操美女的视频在线观看| 亚洲专区国产一区二区| 亚洲,欧美精品.| 日本一区二区免费在线视频| 亚洲欧美一区二区三区久久| 久久狼人影院| 交换朋友夫妻互换小说| 成人影院久久| ponron亚洲| 亚洲欧美一区二区三区黑人| 国产精品欧美亚洲77777| 老司机靠b影院| 日韩欧美免费精品| 欧美黄色片欧美黄色片| cao死你这个sao货| 校园春色视频在线观看| 国产精品一区二区精品视频观看| 18禁国产床啪视频网站| 午夜精品久久久久久毛片777| 欧美乱码精品一区二区三区| 国产又色又爽无遮挡免费看| 香蕉丝袜av| 91麻豆av在线| 水蜜桃什么品种好| 捣出白浆h1v1| 精品无人区乱码1区二区| 热re99久久国产66热| 国产日韩欧美亚洲二区| 国产伦人伦偷精品视频| 国产成人欧美在线观看 | 久久精品人人爽人人爽视色| 欧美日韩福利视频一区二区| 精品少妇久久久久久888优播| 亚洲欧美日韩另类电影网站| 两个人免费观看高清视频| 男人的好看免费观看在线视频 | 宅男免费午夜| 久久精品亚洲av国产电影网| 人妻丰满熟妇av一区二区三区 | 午夜福利一区二区在线看| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 久久香蕉国产精品| 如日韩欧美国产精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 丝袜在线中文字幕| 欧美激情极品国产一区二区三区| 国产成人一区二区三区免费视频网站| 黑人操中国人逼视频| 乱人伦中国视频| 一夜夜www| 国产精品一区二区免费欧美| 婷婷成人精品国产| 极品教师在线免费播放| 国产精品免费大片| 老汉色∧v一级毛片| 18禁国产床啪视频网站| 日韩有码中文字幕| 国产精品乱码一区二三区的特点 | 成年人黄色毛片网站| 伦理电影免费视频| 老熟妇仑乱视频hdxx| 黄频高清免费视频| 国产精品久久电影中文字幕 | 嫩草影视91久久| 男女之事视频高清在线观看| 高清av免费在线| 女同久久另类99精品国产91| 另类亚洲欧美激情| 18在线观看网站| 中文字幕av电影在线播放| 午夜免费成人在线视频| 色精品久久人妻99蜜桃| 在线观看免费高清a一片| 啦啦啦免费观看视频1| tocl精华| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 久久久国产成人精品二区 | 欧美性长视频在线观看| 亚洲视频免费观看视频| 国产麻豆69| 色老头精品视频在线观看| 亚洲视频免费观看视频| 日本wwww免费看| 久久久精品区二区三区| 满18在线观看网站| 丁香六月欧美| 久久天堂一区二区三区四区| 一进一出抽搐动态| 免费日韩欧美在线观看| av欧美777| 午夜精品在线福利| 王馨瑶露胸无遮挡在线观看| 精品国产一区二区三区四区第35| 男女之事视频高清在线观看| xxx96com| 天天躁夜夜躁狠狠躁躁| 免费人成视频x8x8入口观看| 黄色怎么调成土黄色| 欧美av亚洲av综合av国产av| 亚洲人成电影观看| 亚洲成人手机| 国产精品香港三级国产av潘金莲| 如日韩欧美国产精品一区二区三区| 午夜日韩欧美国产| 国产主播在线观看一区二区| 久久香蕉激情| 久久精品亚洲熟妇少妇任你| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 欧美精品亚洲一区二区| 19禁男女啪啪无遮挡网站| 在线av久久热| 麻豆乱淫一区二区| 国产有黄有色有爽视频| 久久影院123| 久久久国产欧美日韩av| 午夜精品在线福利| 热99国产精品久久久久久7| 99精品欧美一区二区三区四区| 国产成人一区二区三区免费视频网站| 手机成人av网站| 黄片小视频在线播放| 国产成人av教育| 亚洲午夜精品一区,二区,三区| 深夜精品福利| 国产国语露脸激情在线看| 欧美性长视频在线观看| 天堂√8在线中文| 伦理电影免费视频| 99久久人妻综合| 成年人黄色毛片网站| 欧美激情久久久久久爽电影 | 欧美久久黑人一区二区| 亚洲国产欧美日韩在线播放| 在线观看午夜福利视频| 欧美人与性动交α欧美软件| 女人久久www免费人成看片| 制服人妻中文乱码| videos熟女内射| 久久精品成人免费网站| 香蕉久久夜色| 99re6热这里在线精品视频| 制服人妻中文乱码| 久久午夜综合久久蜜桃| 中文字幕人妻丝袜制服| 看免费av毛片| 日韩欧美一区二区三区在线观看 | 亚洲av片天天在线观看| 捣出白浆h1v1| 最新美女视频免费是黄的| 亚洲精品久久成人aⅴ小说| 母亲3免费完整高清在线观看| 欧美乱码精品一区二区三区| 后天国语完整版免费观看| 精品福利永久在线观看| 久久国产精品人妻蜜桃| 国产精品1区2区在线观看. | 天天躁日日躁夜夜躁夜夜| 看黄色毛片网站| 999精品在线视频| 国产精品.久久久| 首页视频小说图片口味搜索| 欧美成人午夜精品| 日本a在线网址| 夜夜躁狠狠躁天天躁| 亚洲精品在线观看二区| 免费看十八禁软件| av不卡在线播放| 久久精品国产亚洲av香蕉五月 | 久久人妻av系列| 女人爽到高潮嗷嗷叫在线视频| 国产乱人伦免费视频| 老司机靠b影院| 精品乱码久久久久久99久播| 99精品久久久久人妻精品| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线 | 性少妇av在线| 精品国产乱子伦一区二区三区| 好看av亚洲va欧美ⅴa在| 97人妻天天添夜夜摸| 国产深夜福利视频在线观看| 精品乱码久久久久久99久播| av网站免费在线观看视频| 99热网站在线观看| 9色porny在线观看| 曰老女人黄片| tube8黄色片| 老司机影院毛片| 超色免费av| 高清av免费在线| 一本综合久久免费| 亚洲专区中文字幕在线| 99热只有精品国产| 久久久久国内视频| 熟女少妇亚洲综合色aaa.| 黄色毛片三级朝国网站| 国产亚洲一区二区精品| 精品一区二区三卡| 叶爱在线成人免费视频播放| 女警被强在线播放| 色尼玛亚洲综合影院| 欧美精品啪啪一区二区三区| 亚洲精品在线观看二区| 国产亚洲欧美在线一区二区| 国产色视频综合| 亚洲欧洲精品一区二区精品久久久| 国产精品乱码一区二三区的特点 | 国产精品香港三级国产av潘金莲| 人人妻人人澡人人爽人人夜夜| 亚洲熟女精品中文字幕| 国产日韩一区二区三区精品不卡| 国产乱人伦免费视频| 成年动漫av网址| 俄罗斯特黄特色一大片| 久久久精品国产亚洲av高清涩受| 国产99久久九九免费精品| 动漫黄色视频在线观看| 久久精品国产清高在天天线| 在线看a的网站| 99久久综合精品五月天人人| 一本综合久久免费| 免费黄频网站在线观看国产| 欧美激情高清一区二区三区| 操美女的视频在线观看| 精品卡一卡二卡四卡免费| 18禁观看日本| av有码第一页| 人人妻人人澡人人看| 99国产精品99久久久久| 飞空精品影院首页| 日本黄色日本黄色录像| 久9热在线精品视频| 免费观看精品视频网站| 12—13女人毛片做爰片一| 成年人免费黄色播放视频| 久久香蕉激情| 真人做人爱边吃奶动态| 婷婷精品国产亚洲av在线 | 99久久综合精品五月天人人| 色老头精品视频在线观看| 久久久精品免费免费高清| 亚洲五月天丁香| av天堂久久9| 欧美日韩亚洲高清精品| 欧美激情极品国产一区二区三区| 成年人黄色毛片网站| 丁香六月欧美| 精品人妻1区二区| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性xxxx| 成年人午夜在线观看视频| 日韩欧美一区视频在线观看| 免费看十八禁软件| 精品少妇久久久久久888优播| 久久久久久人人人人人| 精品久久久久久久毛片微露脸| 十分钟在线观看高清视频www| 国产精品久久久久成人av| 国产精品香港三级国产av潘金莲| 一边摸一边抽搐一进一小说 | 天天躁夜夜躁狠狠躁躁| 精品国产一区二区三区久久久樱花| 亚洲成国产人片在线观看| 成年人午夜在线观看视频| 国产精品免费一区二区三区在线 | 免费在线观看黄色视频的| 精品国产国语对白av| 国产真人三级小视频在线观看| 丝袜人妻中文字幕| 免费在线观看日本一区| 人妻丰满熟妇av一区二区三区 | 午夜精品国产一区二区电影| 久久中文字幕一级| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟女毛片儿| 亚洲中文av在线| 欧美一级毛片孕妇| 精品少妇一区二区三区视频日本电影| 亚洲熟妇熟女久久| 欧美黄色淫秽网站| 变态另类成人亚洲欧美熟女 | 少妇裸体淫交视频免费看高清 | 国产精品二区激情视频| 亚洲人成伊人成综合网2020| www日本在线高清视频| 日韩欧美三级三区| 精品一区二区三区四区五区乱码| 欧美国产精品va在线观看不卡| 免费在线观看完整版高清| 亚洲成a人片在线一区二区| 午夜激情av网站| 成年人免费黄色播放视频| 婷婷丁香在线五月| 免费在线观看影片大全网站| 国产伦人伦偷精品视频| 最近最新中文字幕大全电影3 | 捣出白浆h1v1| 久久久久久亚洲精品国产蜜桃av| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 国产激情久久老熟女| 黑人巨大精品欧美一区二区蜜桃| cao死你这个sao货| 日日爽夜夜爽网站| 国产不卡一卡二| 精品熟女少妇八av免费久了| 在线十欧美十亚洲十日本专区| 99精品欧美一区二区三区四区|