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

    北京城市與西北遠(yuǎn)郊地表臭氧濃度梯度移動監(jiān)測研究

    2019-11-07 06:43:40張紅星韓立建任玉芬姚余輝孫旭王效科周偉奇鄭華
    生態(tài)學(xué)報 2019年18期
    關(guān)鍵詞:臭氧濃度西北臭氧

    張紅星,韓立建,任玉芬,姚余輝,孫旭,王效科,周偉奇,鄭華

    中國科學(xué)院生態(tài)環(huán)境研究中心,城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室,北京城市生態(tài)系統(tǒng)研究站, 北京 100085

    臭氧(O3)是空氣的重要微量成分,約90%的臭氧存在于平流層,吸收紫外線,保護(hù)地球生物與環(huán)境,對人類是有利的;約10%的臭氧存在于對流層,具有消菌殺毒的功效,但臭氧濃度持續(xù)增加可能會構(gòu)成大氣污染,對生態(tài)系統(tǒng)和生物健康造成不利影響。20世紀(jì)70年代以來,對流層中臭氧的背景濃度以每年0.5%—2%的速率增長[1]。以北京為例,2004年以來北京市臭氧濃度整體呈現(xiàn)上升趨勢[2-3],城市與區(qū)域臭氧濃度升高以及空氣氧化性增強(qiáng)引發(fā)的環(huán)境問題日益引起關(guān)注[2- 11]。

    從空間分布看,點(diǎn)位監(jiān)測表明北京城區(qū)地表臭氧濃度相對較低,周邊區(qū)縣相對較高[12]。然而飛機(jī)監(jiān)測發(fā)現(xiàn)北京及周邊區(qū)域1.5km以下受近地面人為活動影響較大,在水平方向上,0—2km范圍內(nèi)城區(qū)四環(huán)內(nèi)臭氧濃度高[13]。北京城市與區(qū)域空間臭氧格局需要進(jìn)一步監(jiān)測研究。

    相比空中臭氧,地表臭氧濃度對人與各類生態(tài)系統(tǒng)影響更直接,但是人們尚不清楚在城市與區(qū)域空間上臭氧的濃度分布規(guī)律。此前關(guān)于地表臭氧的研究多集中于點(diǎn)位觀測或聯(lián)網(wǎng)觀測[2,4,10,12,14-15]。然而,因?yàn)辄c(diǎn)位監(jiān)測代表的空間范圍不確定,臭氧在地表的空間格局仍然不清晰。盡管有研究發(fā)現(xiàn)北京西部和北部臭氧平均濃度高,但是這些研究也是基于點(diǎn)位的監(jiān)測數(shù)據(jù),目前缺少臭氧在這些區(qū)域空間上連續(xù)變化規(guī)律的報道。此外有研究表明,植物排放的BVOC在臭氧生成過程中發(fā)揮作用[16],可以使白天的臭氧濃度升高[17],然而植被與臭氧的關(guān)系的觀測研究較少。

    本研究采用北京城市生態(tài)系統(tǒng)研究站移動監(jiān)測的數(shù)據(jù)資料,選擇夏季典型臭氧污染天氣,沿道路多點(diǎn)位水平方向移動監(jiān)測北京城市到西北郊區(qū)、遠(yuǎn)郊臭氧濃度變化特征,用遙感技術(shù)獲取道路沿線植被歸一化指數(shù)NDVI,研究二者之間的關(guān)系,結(jié)合天氣形勢和氣團(tuán)軌跡分析,揭示北京城市與遠(yuǎn)郊地表臭氧遞變規(guī)律及植被機(jī)制。

    1 材料與方法

    1.1 監(jiān)測儀器

    采用臭氧濃度分析儀(49i,美國熱電公司)測定臭氧濃度, 檢測方法為紫外光度法,根據(jù)紫外光經(jīng)過樣品后被吸收的程度計(jì)算出臭氧的濃度,零點(diǎn)漂移<2.14μg/m3/24h,跨度漂移1個月小于1%,最低可檢測限μg/m3,最大量程4285.71μg/m3。每次監(jiān)測前都用Thermo 49i-PS校正49i,監(jiān)測頻率為60s。

    1.2 監(jiān)測方法

    把臭氧分析儀49i搭載在移動監(jiān)測車上,用氣泵通過采樣桿從車頂部室外抽入空氣,從車底部排出空氣,臭氧分析儀49i通過其內(nèi)置氣泵從采樣桿分流口采集氣體樣品,分析完畢通過廢氣通道從車底部排出,采樣管從車頂伸出車外40cm,離地面約3.3m。工控機(jī)實(shí)時記錄監(jiān)測數(shù)據(jù),每60s返回1次監(jiān)測數(shù)值;用開發(fā)的移動監(jiān)測車供電系統(tǒng)給儀器供電;用基于GPS的車輛定位跟蹤系統(tǒng)實(shí)時記錄車輛的經(jīng)度、緯度、速度。經(jīng)緯度的精度約10 m。

    選擇夏季典型無風(fēng)晴天,沿固定線路從北京城內(nèi)到西北遠(yuǎn)郊,行進(jìn)的同時監(jiān)測。從林業(yè)大學(xué)北路出發(fā),經(jīng)由城市支路,三環(huán)路、然后通過京藏高速(G6)出城,往北京西北遠(yuǎn)郊行駛,經(jīng)由昌赤城路到延慶大莊科鄉(xiāng)(圖1)然后返程,經(jīng)由昌赤路、京藏高速、林萃路等,回到出發(fā)地點(diǎn)林業(yè)大學(xué)北路。移動測定的時間是2012年5月10日14:16到19:21;7月24日13:25到18:30;8月24日12:08到18:05。這3天是典型的臭氧污染天氣。表1中詳細(xì)介紹了移動監(jiān)測車具體運(yùn)行時間、道路類型、名稱及代表區(qū)域。為了便于控制時間變化對臭氧濃度變化的影響,把移動監(jiān)測過程依據(jù)監(jiān)測車運(yùn)行的地段和相應(yīng)時間分段(表1)。

    在監(jiān)測車移動監(jiān)測臭氧濃度的同時,對照點(diǎn)北京城市生態(tài)系統(tǒng)研究站的教學(xué)植物園監(jiān)測平臺(116°25′39.36″,39°52′22.48″)同步監(jiān)測臭氧濃度和氣象要素,測定頻率60s。每次監(jiān)測前,都將移動監(jiān)測車的臭氧分析儀、GPS、對照點(diǎn)臭氧分析儀時間同步。

    表1 監(jiān)測線路及測定時間

    圖1 生態(tài)監(jiān)測車運(yùn)行路線圖Fig.1 Monitoring route of mobile vehicle

    1.3 數(shù)據(jù)處理

    用Origin Pro8繪制觀測期間的溫度、濕度、風(fēng)矢量變化圖;用SPSS13.0軟件,運(yùn)行ANOVA程序,分析不同路段間臭氧濃度的差異顯著性,運(yùn)行Independent samples T-test分析對照點(diǎn)和同時間范圍內(nèi)移動監(jiān)測數(shù)據(jù)之間差異顯著性;用ArcGIS10.0軟件把監(jiān)測車軌跡和臭氧濃度相對高低標(biāo)示在地圖上;用HYSPLIT模型做氣團(tuán)質(zhì)點(diǎn)后向軌跡分析,以沙塘溝(116°10′11.35″,40°24′12.78″)、大莊科(116°14′59.23″,40°23′3.23″)、十三陵(116°12′26.82″,40°14′0.10″)為氣團(tuán)后向軌跡的終點(diǎn),氣團(tuán)運(yùn)行時間為6 h,設(shè)置氣團(tuán)高度為700.00m;用GrADS 2.1分析典型臭氧污染情況下的天氣形勢,氣象數(shù)據(jù)選取NOAA數(shù)據(jù)庫2012年5月10日、7月24日、8月24日北京時間下午16時的數(shù)據(jù)。

    利用移動監(jiān)測獲取的點(diǎn)位經(jīng)緯度數(shù)據(jù),在ArcGIS支持下生成50m緩沖區(qū),基于LandSAT7衛(wèi)星遙感數(shù)據(jù),計(jì)算緩沖區(qū)內(nèi)NDVI的平均值,7月份衛(wèi)星數(shù)據(jù)不可用,鑒于植被在生長盛期變化不明顯,故用8月份的數(shù)據(jù)替代進(jìn)行NDVI計(jì)算。

    2 結(jié)果與分析

    2.1 監(jiān)測期間的環(huán)境因子變化

    2.1.1溫度日變化

    本研究中,對照點(diǎn)北京教學(xué)植物園的最高溫度都出現(xiàn)在14:00到15:00之間,且都大于30℃,在監(jiān)測車運(yùn)行時段的平均溫度都大于28.5℃。2012年5月10日12:00—20:00,對照點(diǎn)的平均氣溫是29.05℃,最高氣溫30.35℃,出現(xiàn)在14:50;7月24日的平均氣溫是31.09℃,最高氣溫32.47℃,出現(xiàn)在14:20;8月24日平均氣溫是28.52℃,最高氣溫30.66℃,出現(xiàn)在14:40(圖2)。

    圖2 監(jiān)測期間環(huán)境因子變化Fig.2 The environmental factors in monitoring days20120510:2012年5月10日,May 10, 2012;20120724:2012年7月24日,July 24, 2012;20120824:2012年8月24日,August 24, 2012

    2.1.2相對濕度日變化

    3次臭氧典型污染天氣監(jiān)測期間的平均濕度都低于50%,最低濕度都出現(xiàn)在14:00—14:30之間,和最高氣溫的出現(xiàn)時間段基本重合。2012年5月10日12:00—20:00對照點(diǎn)的平均濕度是36.42%,最低是32.11%,出現(xiàn)在14:00;7月24日的平均濕度是49.91%,最低是43.22%,出現(xiàn)在14:10;8月24日的平均濕度是33.32%,最低是27.08%,出現(xiàn)在14:20。

    2.1.3風(fēng)矢量日變化

    5月10日的平均風(fēng)速是0.52m/s,7月24日是0.52m/s,8月24日是0.69m/s。這3次典型臭氧污染天氣時,晚上幾乎是靜風(fēng)狀態(tài),10:00后整體風(fēng)向?yàn)槟巷L(fēng)。12:00后平均風(fēng)速都小于1.2 m/s(圖2)。2012年5月10日12:00—20:00對照點(diǎn)風(fēng)向?yàn)槲髂巷L(fēng),平均風(fēng)速為0.97m/s,最大風(fēng)速1.72 m/s,出現(xiàn)在15:10,最小風(fēng)速0.42 m/s,出現(xiàn)在16:50;7月24日12:00—20:00的風(fēng)向?yàn)槲髂虾湍巷L(fēng),平均風(fēng)速為0.93 m/s,最大1.70 m/s,出現(xiàn)在14:50,最小0.39 m/s,出現(xiàn)在12:20;8月24日風(fēng)向?yàn)槟虾臀髂?12:00—20:00的平均風(fēng)速為1.19 m/s,最大為2.05 m/s,出現(xiàn)在16:30,最小為0.38 m/s,出現(xiàn)在20:00。

    2.1.4光強(qiáng)日變化

    圖3 典型臭氧污染天氣臭氧濃度日變化過程 Fig.3 The diurnal change of ozone concentrations in typical pollution days

    5月10日、7月24日、8月24日都是晴天,光強(qiáng)呈現(xiàn)典型的單峰日變化特征,地表接受的光合有效輻射分別是30.93mol/m2.d、31.95mol/m2.d、40.35mol/m2.d。

    2.1.5臭氧濃度日變化

    臭氧濃度的日變化特征整體表現(xiàn)為“S”型,從零時起濃度保持低水平或者平緩下降,甚至在夜間可以回到0.00μg/m3。但是,從 7:00開始臭氧濃度開始快速上升。在5月10日,于16:00達(dá)到一天中的最高峰280.18μg/m3,然后開始下降;7月14日,13:00到達(dá)了較高水平,高峰從13:00持續(xù)到17:00,約200.00—210.00μg/m3;8月24日,12:00達(dá)到較高水平,高峰從12:00持續(xù)到17:00然后開始下降,約190.00—205.00μg/m3。高濃度的臭氧甚至可以持續(xù)到晚上21:00以后。

    2.2 臭氧濃度和氣象因子的關(guān)系

    本研究中,相關(guān)分析(Pearson Correlation)表明,臭氧濃度在典型污染天氣和空氣溫度成顯著正相關(guān),和空氣濕度成顯著負(fù)相關(guān),和光強(qiáng)成顯著正相關(guān)。

    表2 臭氧濃度和氣象因子的關(guān)系

    ** 極顯著(P<0.0001)

    2.3 典型污染期間地表臭氧濃度城市郊區(qū)梯度特征

    整體平均而言,北京西北山區(qū)植被覆蓋良好區(qū)域的臭氧濃度高于北京城區(qū)和北京西北部山前的平原區(qū)域。在選擇的固定監(jiān)測線路上,從十三陵往西北山區(qū)方向,沿昌赤路,到延慶大莊科鄉(xiāng),山區(qū)道路上監(jiān)測到的臭氧濃度明顯高于在京藏高速和北京城區(qū)三環(huán)主路的監(jiān)測值(圖4)。在3次監(jiān)測過程中有一個共同現(xiàn)象,十三陵是地表臭氧濃度的分界點(diǎn)。從十三陵開始,往植被茂密的山區(qū)方向臭氧濃度陡然升高,黃色變成紅色(圖4),并保持在高水平。返程時往城區(qū)方向,車輛駛出十三陵區(qū)域后臭氧濃度顯著降低(圖4),紅色變成黃色。西北山區(qū)(Ⅳ、Ⅴ)臭氧濃度約是高速公路上(Ⅲ、Ⅵ)的2.28到3.33倍,是平原城市區(qū)域平均水平約2.00倍。

    圖4 北京城市和郊區(qū)地表臭氧濃度連續(xù)變化特征Fig.4 Surface ozone concentrations gradient between urban and rural region of Beijing

    2.4 不同時間段,城市和郊區(qū)不同地點(diǎn)臭氧濃度的差別

    在三次移動監(jiān)測中(表3,圖5),西北山區(qū)相鄰約兩個小時內(nèi)監(jiān)測到的兩組數(shù)值(Ⅳ、Ⅴ)差別不顯著(P<0.001),平均為254.68μg/m3,但顯著高于北京城區(qū)監(jiān)測到的臭氧濃度127.08μg/m3(P<0.001)。西北山區(qū)移動監(jiān)測到的平均值是城區(qū)平均值的2.00倍,最大值是城區(qū)監(jiān)測到的最小值的11.14倍。西北山區(qū)的臭氧濃度平均是城區(qū)支路(Ⅰ)的1.41—1.76倍,是城區(qū)交通主干道(Ⅱ)的2.09倍,是高速路(Ⅲ)的2.29—3.32倍。

    在移動監(jiān)測中,除了西北山區(qū)(Ⅳ、Ⅴ)的臭氧濃度外,對照點(diǎn)在Ⅰ、Ⅱ、Ⅲ、Ⅵ、Ⅶ等各地段測得的臭氧濃度均值都顯著高于(P<0.001)各時段在不同路段監(jiān)測的臭氧濃度(表3)。西北山區(qū)臭氧和對照點(diǎn)濃度差別最小,是城區(qū)對照點(diǎn)的1.11倍。考慮到對照點(diǎn),臭氧濃度從高到低的排序是:西北山區(qū)(Ⅳ、Ⅴ)>對照點(diǎn)(教學(xué)植物園)>城市支路(Ⅰ、Ⅶ)>城市快速路和高速路(Ⅱ、Ⅲ、Ⅵ)。盡管如此,值得注意的是監(jiān)測過程中發(fā)現(xiàn),監(jiān)測車靠近道旁樹時監(jiān)測到臭氧高值,沒有道旁樹的路段臭氧明顯降低。此外發(fā)現(xiàn),城市支路上臭氧濃度也有高于對照點(diǎn)的情況,對照點(diǎn)濃度也有高于山區(qū)臭氧的情況(圖5)。

    同一類型監(jiān)測地段,臭氧濃度隨時間變化。城市支路(Ⅰ)12:08—13:55時間段監(jiān)測到的臭氧濃度平均為181.06μg/m3,顯著高于(Ⅶ)17:47—19:21時間段監(jiān)測到的臭氧濃度144.60μg/m3。高速公路上(Ⅲ)14:48—16:36時間段監(jiān)測到的臭氧濃度111.13μg/m3,顯著高于(Ⅵ)17:20—19:05時間段監(jiān)測到的值。這種現(xiàn)象符合臭氧從日間到黑夜?jié)舛戎饾u降低的日變化規(guī)律。

    相近時間段內(nèi),城市快速路和高速路臭氧濃度差別不顯著。13:14—16:08時間段在城市快速路三環(huán)上(Ⅱ)測得的城市平均臭氧濃度121.99μg/m3,和14:48—16:36時間段在京藏高速公路出京方向(Ⅲ)測得的平均臭氧濃度111.13μg/m3差異不顯著(P<0.001)。

    表3 北京城市和郊區(qū)不同地段監(jiān)測到的地表臭氧濃度/(μg/m3)

    圖5 移動監(jiān)測和對照點(diǎn)臭氧濃度值Fig.5 The ozone concentrations monitored using mobile vehicle contrast to control pointM-O3:移動監(jiān)測車測得臭氧濃度;C-O3:對照點(diǎn)北京教學(xué)植物園同時間測得臭氧濃度

    2.5 臭氧濃度和植被的關(guān)系

    研究發(fā)現(xiàn),臭氧濃度和監(jiān)測點(diǎn)周圍50米范圍內(nèi)的NDVI呈顯著正相關(guān)關(guān)系,隨著NDVI的增大,監(jiān)測到的臭氧濃度呈現(xiàn)出Logistic增長(圖6)。當(dāng)監(jiān)測車途經(jīng)區(qū)域的植被NDVI指數(shù)較大,植物茂密時,監(jiān)測到高濃度臭氧概率較大。對照點(diǎn)教學(xué)植物園的臭氧濃度和NDVI的關(guān)系同樣能夠通過Logistic方程得到較好擬合(圖6)。

    圖6 監(jiān)測到的臭氧濃度和采樣點(diǎn)周圍植被歸一化指數(shù)(NDVI)的關(guān)系Fig.6 The relationship between ozone concentrations and NDVI around monitoring sites

    2.6 城鄉(xiāng)臭氧濃度空間變異及對照點(diǎn)時間變異

    從北京城區(qū)到西北遠(yuǎn)郊,山區(qū)地表臭氧濃度的空間變異系數(shù)最小,高速公路、城市快速路、支路空間變異較大(表4)。不同時段對應(yīng)路段監(jiān)測的臭氧濃度的空間變異系數(shù)都明顯大于同時間段對照點(diǎn)的時間變異系數(shù)(表4)。在3次監(jiān)測中,城市與區(qū)域監(jiān)測的所有臭氧濃度的空間變異系數(shù)平均為0.41,監(jiān)測期間對照點(diǎn)臭氧濃度時間變異系數(shù)為0.06,說明監(jiān)測到的臭氧濃度地點(diǎn)間的差異主要是由于空間不同造成,時間變化不是引起地點(diǎn)間差異的主要原因。

    表4 移動監(jiān)測臭氧濃度空間變異系數(shù)和同時間范圍對照點(diǎn)臭氧濃度時間變異系數(shù)

    Table 4 The coefficients of spatial variation of ozone monitored using mobile vehicle and the coefficients of variation in control point at the same time

    區(qū)域Region地段Sites時間Time5-107-248-24平均MeanCK5-107-248-24平均Mean 城區(qū) Urban regionⅠ0.320.290.120.24 0.030.060.020.04城區(qū)Urban regionⅡ0.330.210.230.26 0.040.050.010.03高速路Free wayⅢ0.390.280.270.31 0.020.040.010.02西北山區(qū)North-west mountain regionⅣ0.140.120.070.11 0.030.090.020.05西北山區(qū) North-west mountain regionⅤ0.100.110.110.11 0.030.050.030.04高速路 Free wayⅥ0.400.400.430.41 0.020.050.020.03城區(qū)Urban regionⅦ0.340.260.370.32 0.010.040.020.02全部0.440.370.430.410.060.080.040.06

    西北山區(qū)的空間變異最小,約為0.11;城市支路的空間變異系數(shù)為0.24—0.32;高速路為0.31—0.41;城市快速路為0.26。

    2.7 典型污染期間郊區(qū)臭氧后向軌跡分析及天氣形勢

    在5月10日,7月24日,8月24日的3次典型污染天氣時,從16時開始,700m海拔高度上6 h氣團(tuán)后向軌跡結(jié)果表明,西北部山區(qū)的氣團(tuán)來至北京城區(qū)方向(圖7)。5月10日和7月24日的氣團(tuán)主要來至北京城區(qū),8月24日的氣團(tuán)來至城區(qū)西部,沿太行山東麓向北遷移。

    在2012年5月10典型臭氧污染天氣發(fā)生時,天氣形勢為大范圍均壓場(圖7),北京大部分位于均壓場內(nèi);7月24日和8月24日典型臭氧污染發(fā)生時,北京城市與區(qū)域位于高壓后部,低壓前部,由弱高壓控制(圖7)。

    圖7 典型臭氧污染天氣時的氣團(tuán)后向軌跡和天氣形勢圖Fig.7 The synoptic types and backward trajectories in typical ozone pollution weather圖中每條軌跡上的點(diǎn)是后向3 h氣團(tuán)質(zhì)點(diǎn)所在位置

    3 討論

    3.1 臭氧污染的氣象原因

    溫度、濕度、光照共同影響臭氧的生成和轉(zhuǎn)化。我們研究發(fā)現(xiàn),臭氧濃度和溫度正相關(guān),相關(guān)系數(shù)0.75。嚴(yán)茹莎等[10]研究發(fā)現(xiàn),臭氧濃度和溫度顯著正相關(guān),相關(guān)系數(shù)達(dá)到0.74。研究結(jié)論基本一致。同時,我們發(fā)現(xiàn)臭氧濃度和空氣濕度顯著負(fù)相關(guān),相關(guān)系數(shù)0.80,這和程念亮等[7]的發(fā)現(xiàn)相似。本研究中臭氧濃度和光強(qiáng)呈現(xiàn)顯著正相關(guān),相關(guān)系數(shù)為0.40。

    盡管有研究認(rèn)為臭氧濃度和風(fēng)速正相關(guān)[7],但更多證據(jù)支持適宜的風(fēng)速(≤3.00m/s)是造成近地層臭氧濃度高的必要條件。在風(fēng)速較低的情況下,大氣層結(jié)構(gòu)穩(wěn)定,上下混合均勻,臭氧濃度垂直廓線較平穩(wěn)[18]。本研究中,3次典型污染天氣當(dāng)天12時到20時,地表監(jiān)測到的都為平均風(fēng)速低于1.20m/s的南風(fēng)(圖2)。5月10日的氣流方向?yàn)闁|南向,7月24日的為南向氣流,8月24日的為西南向氣流(圖7)。低于1.2m/s的弱南風(fēng)是三次典型臭氧污染天氣的共同氣象條件。

    本研究中,5月10日天氣形勢為均壓場控制,7月24日和8月24日為高壓后部,低壓前部。研究認(rèn)為,北京處于低壓前部控制時,邊界層穩(wěn)定[10]本地排放的污染物不容易擴(kuò)散,高溫低濕的氣象條件有利于本地臭氧的生成[7,9]。程念亮等[7]研究發(fā)現(xiàn)臭氧超標(biāo)日地面氣壓形勢場高壓類、 低壓類、 均壓類等各占16%,36%,48%,以均壓、低壓前部、弱高壓為主。本研究選擇的天氣型和前人研究比較一致,都是典型污染天氣形勢。不同的是,本研究發(fā)現(xiàn)在均壓場控制情況下,北京城市與區(qū)域的整體臭氧污染水平相對更高(圖5、圖7)。原因可能是在均壓控制下,大氣層結(jié)更穩(wěn)定,風(fēng)速較低,不利于各種污染物的擴(kuò)散,空氣濕度小,溫度高,光化學(xué)反應(yīng)劇烈,造成城市與區(qū)域整體臭氧濃度偏高。

    3.2 北京城市與區(qū)域臭氧空間分布特征

    研究表明位于北京周邊平谷、密云、昌平、延慶、房山南部、通州東南部周邊區(qū)縣的臭氧濃度高于城區(qū)的臭氧濃度,北京周邊對照點(diǎn)臭氧濃度高于郊區(qū),郊區(qū)高于城區(qū),城區(qū)高于交通污染監(jiān)測點(diǎn)[12]。本研究發(fā)現(xiàn),在西北遠(yuǎn)郊山區(qū)森林區(qū)域的臭氧濃度從十三陵開始臭氧濃度陡然升高,從高到低的順序依次是西北山區(qū)(Ⅳ、Ⅴ)>對照點(diǎn)教學(xué)植物園>城市支路(Ⅰ、Ⅶ)>城市快速路和高速路(Ⅱ、Ⅲ、Ⅵ)。從平均水平而言,山區(qū)大于城區(qū),城區(qū)大于交通沿線,這和前人研究比較一致。不同的是,我們發(fā)現(xiàn)城市與區(qū)域臭氧空間變化比較復(fù)雜,在不同地點(diǎn)間存在較大變異。在城市支路上臭氧濃度也有高于對照點(diǎn)的情況,對照點(diǎn)濃度也有高于山區(qū)臭氧的情況(圖5)。

    以前的研究發(fā)現(xiàn)北京西部和北部山區(qū)的臭氧平均濃度高于平原地區(qū);在平原地區(qū),公園中的臭氧濃度高于道旁綠化帶[15],十三陵的定陵處臭氧濃度高于北京城區(qū)的平均臭氧濃度,是北京城區(qū)的1.01到1.56倍[2]。但點(diǎn)位監(jiān)測的方法不能發(fā)現(xiàn)大范圍內(nèi)臭氧空間變化的特征。我們用移動監(jiān)測的辦法發(fā)現(xiàn)西北山區(qū)的平均值是城區(qū)平均值的2.00倍,是城區(qū)對照點(diǎn)北京教學(xué)植物園的1.10倍,山區(qū)最大值是城區(qū)監(jiān)測到的最小值的11.14倍。此外,發(fā)現(xiàn)進(jìn)入西北植被茂密區(qū)后臭氧濃度陡然升高到平原區(qū)域約2.00倍并在山區(qū)保持在高位的現(xiàn)象,這是前人用傳統(tǒng)固定點(diǎn)位監(jiān)測無法發(fā)現(xiàn)的科學(xué)現(xiàn)象,其機(jī)制值得深入研究。

    3.3 遠(yuǎn)程輸送對臭氧污染的影響

    很多科學(xué)家都發(fā)現(xiàn)了遠(yuǎn)程輸送對臭氧污染的影響。蘇福慶等[19]將影響北京外來污染物輸入的風(fēng)帶分為3 種類型,即東風(fēng)帶、西南風(fēng)帶及東南風(fēng)帶輸送通道。認(rèn)為3 個輸送風(fēng)帶的輸入, 山前輸送匯系統(tǒng)對大范圍污染物的匯聚,是造成北京外來污染物輸入的主要原因。唐貴謙等認(rèn)為復(fù)合氣流易將河北、天津乃至山東、山西等地的污染物傳輸?shù)奖本?加之西部和北部山體的阻擋極易造成污染物的積累[14]。王占山等觀測到了一次西南風(fēng)作用下的明顯的臭氧輸送過程,榆垡、豐臺花園、奧體中心、懷柔監(jiān)測站臭氧峰值出現(xiàn)的時間從南到北滯后[12]。徐敬等[11]研究發(fā)現(xiàn)北京地區(qū)受西南氣流影響時,會導(dǎo)致下風(fēng)向的北部山區(qū)出現(xiàn)極高的臭氧濃度。安俊琳等[9]研究表明北京盛行東南、 偏南和偏西氣流時,容易造成高臭氧濃度。本研究中,5月10日典型臭氧污染天氣時,發(fā)生了均壓場控制下的東南向弱氣流,氣團(tuán)從東南向西北輸送(圖7);7月24日和8月24日天氣形勢是高壓后部,低壓前部控制下的南向弱氣流。輸送通道和前人的研究基本相一致,即發(fā)生典型臭氧污染時,均伴隨著南向弱氣流,氣流沿西山前從南向北輸送(圖7),但沒有發(fā)現(xiàn)東風(fēng)帶輸入氣流。然而,即便存在西南風(fēng)和東南風(fēng)輸送,本研究中后向軌跡分析表明,監(jiān)測到的高濃度臭氧可能不是當(dāng)天遠(yuǎn)程輸送而來。按照后向運(yùn)動軌跡,16:00時到19:00時西部山區(qū)的氣團(tuán)質(zhì)點(diǎn)來源于城區(qū)上空10:00時到13:00時(圖7)。以城區(qū)臭氧相對較高的對照點(diǎn)為例,該時間段內(nèi)城市臭氧濃度還沒有達(dá)到一天中最高的程度(圖3),且低于西北山區(qū)監(jiān)測到的平均值。即便城市中臭氧濃度達(dá)到了一天中的最高值,其仍然比郊區(qū)監(jiān)測到的平均值低(圖5)。所以,本研究中監(jiān)測到的高濃度臭氧可能不是從城區(qū)輸送來的。

    但是,北京對西北遠(yuǎn)郊森林區(qū)域確實(shí)有物質(zhì)輸送(圖7)。城市地區(qū)的氣團(tuán)在向下風(fēng)向傳輸?shù)倪^程中經(jīng)過“老化”可導(dǎo)致下風(fēng)向的臭氧濃度水平高于市區(qū)[20]。城市道路和生產(chǎn)生活排放的臭氧前體物NOx和VOCs遠(yuǎn)程運(yùn)動到北京北部山區(qū),在溫度較高,濕度較小的氣象條件下,可能因?yàn)樯种参锱欧诺腣OCs的參與[21],光化學(xué)反應(yīng)產(chǎn)生臭氧并大量富集在該區(qū)域。此外,由于遠(yuǎn)郊相對缺乏NO污染物,臭氧消耗較少,加上地形的影響可能存在累及效應(yīng)。

    3.4 高濃度臭氧和植被的可能關(guān)系

    通過移動中多點(diǎn)連續(xù)觀測發(fā)現(xiàn)北京西北山區(qū)地表臭氧濃度平均值高且空間變異小,而城區(qū)臭氧濃度平均值相比西北山區(qū)小且空間變異比較大。盡管西北山區(qū)平均臭氧水平高于北京城區(qū),但也觀測到了城區(qū)教學(xué)植物園臭氧濃度高于西北山區(qū)的情況,特別是在城市支路上距離道旁樹比較近的地方也觀測到了臭氧濃度高值(圖5)。前人的觀測多基于單位時間平均值,這種統(tǒng)計(jì)方法的好處是便于宏觀展現(xiàn)格局和過程,弊端是對發(fā)生在植物空氣界面附近的生理生化以及光化學(xué)過程監(jiān)測不到。比如,盡管從平均情況看道旁綠化帶的臭氧濃度低于公園[15],但是,在距離道旁綠化帶近的地方,確實(shí)觀測到了多組高于公園的值(圖5)。

    本研究中,除了西北森林區(qū)域外,在對照點(diǎn)教學(xué)植物園和城市道路支線等植物茂密區(qū)域也監(jiān)測到了臭氧高值,在一定時間范圍內(nèi),臭氧濃度和植被NDVI指數(shù)呈現(xiàn)正相關(guān),臭氧濃度隨NDVI指數(shù)Logistic增長。這一結(jié)果證明了在植被覆蓋良好的區(qū)域臭氧濃度高[12]。植物排放的BVOCs在大氣光化學(xué)反應(yīng)中發(fā)揮重要作用可能是植被良好地區(qū)臭氧濃度高的原因[22]。Trainer 等[21]認(rèn)為大氣污染 80%來自植物釋放的碳?xì)浠衔?。事?shí)上,在美國 亞特蘭大地區(qū),政府用10年的時間控制了人為污染物排放后,大氣環(huán)境質(zhì)量并沒有明顯改善,原因是大量松樹(主要釋放單萜烯)和橡樹(主要釋放異戊二烯)是BVOCs的主要來源,在合適的光照和濕度下,它們與人為排放的 NOx 共同作用產(chǎn)生了臭氧,從而大大增加城市發(fā)生光化學(xué)煙霧,最終導(dǎo)致控制臭氧水平失敗[22]。本研究中,西北山區(qū)植被覆蓋率達(dá)到95%以上,遠(yuǎn)郊森林區(qū)域有473種植物,他們放出大量的BVOCs[23];對照點(diǎn)北京教學(xué)植物園,景觀以植被為主,距離南部二環(huán)主路約200 m,在南風(fēng)的作用下(圖2),道路上機(jī)動車排放的NOx和VOCs進(jìn)入教學(xué)植物園。NO2光解導(dǎo)致 O3的生成,植物排放的BVOCs氧化生成活性自由基,尤其是HO2·、RO2·等,HO2·、RO2·引起NO向NO2轉(zhuǎn)化,一方面提供了生成O3的NO2源,更重要的是破壞NO2-NO-O3的光解循環(huán),使臭氧累積[24],這可能是對照點(diǎn)臭氧濃度高的原因之一。在北京城市道路的支路上有長勢良好的樹木白蠟和國槐,可能參與了臭氧生成的光化學(xué)反應(yīng)過程,造成植物附近臭氧濃度高。植物源VOCs、化石源VOCs、NOx等在植物參與下,在植物空氣界面進(jìn)行的生理學(xué)和光化學(xué)反應(yīng)和臭氧生成過程的關(guān)系,需要加強(qiáng)模擬和觀測研究。

    3.5 臭氧濃度空間變異的可能原因

    本研究中地表臭氧在城鄉(xiāng)地表的平均變異系數(shù)約為0.41,在山區(qū)平均為0.11,在高速路上為0.36,在城市支路上平均為0.27。道路交通狀況和道旁植被多寡可能會影響臭氧生成和消耗過程而影響空間變異。車輛排放出NO和VOCs,道路等級越高,車流量越大,還原性氣體越多,越容易造成臭氧消耗,城市中臭氧晚上甚至可以消耗完(圖3);另一方面,道路等級越低,監(jiān)測車離行道樹越近,緩沖范圍內(nèi)植被NDVI越高,越容易監(jiān)測到高濃度臭氧??赡苁且?yàn)橹参锱欧诺漠愇於┗蜉葡╊怴OCs在光作用下參與了臭氧生成,促使更多的NO2釋放出帶電氧原子,帶電氧原子和氧分子生成臭氧[24]。

    地表景觀均勻度和過程差異可能會影響臭氧的空間格局。山區(qū)植被茂密,地表特征相對均勻,車輛稀少,臭氧的生成和消耗空間差別不大。但是,高速路、城市道路經(jīng)過的區(qū)域,車流變化大,建筑高度和分布各異,擴(kuò)散條件差別大,道旁植被格局及距離車道遠(yuǎn)近不同,這些都可能對臭氧的生成和遷移轉(zhuǎn)化產(chǎn)生影響。城市不同景觀和干擾對臭氧的積累耗散,遷移轉(zhuǎn)化需要更多研究加以證明。

    4 結(jié)論

    (1)臭氧濃度和監(jiān)測路線周邊植被NDVI指數(shù)正相關(guān),隨其Logistic增長。地表臭氧濃度在山區(qū)植被茂密地區(qū)陡然升高到平原城市區(qū)域的約2.00倍且在山區(qū)保持此高水平,界限明顯。但在城區(qū)道旁樹或公園綠地也監(jiān)測到了高濃度臭氧,其比遠(yuǎn)郊植被茂密地區(qū)的臭氧略低??臻g上,北京臭氧濃度從高到低的順序是西北山區(qū)(Ⅳ、Ⅴ)>對照點(diǎn)教學(xué)植物園(城區(qū))>城市支路(Ⅰ、Ⅶ)>城市快速路和高速路(Ⅱ、Ⅲ、Ⅵ)。

    (2)山區(qū)地表臭氧濃度的空間變異系數(shù)最小,為0.11;高速公路最大,為0.31—0.41;、城市快速路約為0.26,城市支路空間變異較大0.24—0.32。

    (3)夏季大范圍均壓場、高壓后部、低壓前部等穩(wěn)定天氣型易造成及南向弱氣流,能夠把空氣污染物輸送到西北遠(yuǎn)郊,但本研究中監(jiān)測到的西北山區(qū)高濃度臭氧并非當(dāng)天從城區(qū)輸送而來。城市與區(qū)域產(chǎn)生的臭氧前體物輸送到西北山區(qū)與植被共同作用可能是北京西北部山區(qū)夏季典型天氣臭氧濃度高的原因之一。

    猜你喜歡
    臭氧濃度西北臭氧
    基于梯度提升算法的近地面臭氧濃度估算比較
    不同品種小麥灌漿期旗葉光合特性及光合基因表達(dá)對臭氧濃度升高的響應(yīng)
    文印室內(nèi)臭氧散發(fā)實(shí)測分析及模擬
    一座西北小城
    西北不惑
    黔西北行吟
    散文詩(2017年15期)2018-01-19 03:07:58
    看不見的污染源——臭氧
    沈陽市臭氧濃度時空分布與變化特征分析
    利用臭氧水防治韭菜遲眼蕈蚊
    中國蔬菜(2016年8期)2017-01-15 14:23:43
    臭氧分子如是說
    午夜精品一区二区三区免费看| 国产三级在线视频| 国国产精品蜜臀av免费| 99精品久久久久人妻精品| 99热这里只有精品一区| 91在线观看av| 免费看a级黄色片| 我的女老师完整版在线观看| 给我免费播放毛片高清在线观看| 校园春色视频在线观看| 天天一区二区日本电影三级| 自拍偷自拍亚洲精品老妇| 亚洲国产欧美人成| 国内精品久久久久久久电影| netflix在线观看网站| 亚洲中文字幕一区二区三区有码在线看| 在线观看一区二区三区| 精品免费久久久久久久清纯| 国产精品伦人一区二区| 久久精品国产亚洲av涩爱 | 免费看光身美女| 日韩精品中文字幕看吧| 九九在线视频观看精品| 午夜福利在线观看免费完整高清在 | 最近视频中文字幕2019在线8| 色播亚洲综合网| 一级黄片播放器| 日本黄色片子视频| 99久国产av精品| 亚洲av不卡在线观看| 国内久久婷婷六月综合欲色啪| 天堂网av新在线| 91狼人影院| av福利片在线观看| 22中文网久久字幕| 国产伦精品一区二区三区视频9| 精品午夜福利在线看| 嫁个100分男人电影在线观看| 日本三级黄在线观看| 成人国产综合亚洲| 国产伦精品一区二区三区视频9| 少妇猛男粗大的猛烈进出视频 | 三级毛片av免费| 99精品在免费线老司机午夜| 嫩草影视91久久| 男女边吃奶边做爰视频| 欧美色欧美亚洲另类二区| 久久久久久久久中文| 久久久久久九九精品二区国产| 亚洲熟妇中文字幕五十中出| 亚洲人与动物交配视频| 99热这里只有是精品50| 性欧美人与动物交配| 日韩欧美国产一区二区入口| 欧美日韩黄片免| 日日撸夜夜添| 亚洲欧美日韩无卡精品| 99久久精品热视频| 日韩欧美一区二区三区在线观看| 国产免费av片在线观看野外av| 欧美一区二区国产精品久久精品| 桃红色精品国产亚洲av| 在线免费观看不下载黄p国产 | 不卡一级毛片| 国产精品无大码| 国产精品日韩av在线免费观看| 成人性生交大片免费视频hd| 亚洲国产日韩欧美精品在线观看| 国产男靠女视频免费网站| 最近视频中文字幕2019在线8| 日本黄色片子视频| 波野结衣二区三区在线| 久久精品国产亚洲av香蕉五月| 中亚洲国语对白在线视频| 欧美国产日韩亚洲一区| 国产一级毛片七仙女欲春2| 色综合亚洲欧美另类图片| 国模一区二区三区四区视频| 精品国内亚洲2022精品成人| 亚洲真实伦在线观看| 最近视频中文字幕2019在线8| 麻豆av噜噜一区二区三区| 久久久久免费精品人妻一区二区| 一a级毛片在线观看| а√天堂www在线а√下载| 91麻豆av在线| 日本成人三级电影网站| 久久久久久国产a免费观看| 亚洲人成伊人成综合网2020| 欧美日本视频| 亚洲成人中文字幕在线播放| 日韩av在线大香蕉| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区高清视频在线| 嫩草影院精品99| 国产亚洲欧美98| 国产探花极品一区二区| 精品国产三级普通话版| 1024手机看黄色片| 亚洲无线观看免费| avwww免费| 午夜福利在线观看吧| 尤物成人国产欧美一区二区三区| 一区二区三区激情视频| 精品久久久久久久久av| 老司机午夜福利在线观看视频| 国产毛片a区久久久久| 久久这里只有精品中国| 国产淫片久久久久久久久| 真人做人爱边吃奶动态| 成人三级黄色视频| 日韩强制内射视频| 精品欧美国产一区二区三| 在线看三级毛片| 国产精品永久免费网站| 国产精品人妻久久久久久| 在线a可以看的网站| av在线观看视频网站免费| 黄色女人牲交| 偷拍熟女少妇极品色| 精品一区二区三区视频在线| 国产精品永久免费网站| 免费看光身美女| 波多野结衣高清无吗| 久久中文看片网| 69av精品久久久久久| 美女 人体艺术 gogo| 精品一区二区三区av网在线观看| 特级一级黄色大片| 三级国产精品欧美在线观看| 3wmmmm亚洲av在线观看| 狂野欧美激情性xxxx在线观看| 成年免费大片在线观看| 小说图片视频综合网站| 亚洲精品色激情综合| 别揉我奶头~嗯~啊~动态视频| 天堂av国产一区二区熟女人妻| 亚洲精品一区av在线观看| 日日干狠狠操夜夜爽| 九九爱精品视频在线观看| 精品不卡国产一区二区三区| 亚洲内射少妇av| aaaaa片日本免费| 久久久久国产精品人妻aⅴ院| 国产免费av片在线观看野外av| 麻豆国产av国片精品| 在线播放无遮挡| avwww免费| 免费观看的影片在线观看| av专区在线播放| 村上凉子中文字幕在线| 99久久精品国产国产毛片| 国产伦精品一区二区三区四那| 国产三级在线视频| 亚洲最大成人手机在线| 国产探花极品一区二区| 国产亚洲91精品色在线| 又爽又黄无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 亚洲18禁久久av| 非洲黑人性xxxx精品又粗又长| 黄色一级大片看看| 国产91精品成人一区二区三区| 日韩 亚洲 欧美在线| 欧美人与善性xxx| 婷婷六月久久综合丁香| xxxwww97欧美| 日本黄色视频三级网站网址| 最近中文字幕高清免费大全6 | 2021天堂中文幕一二区在线观| 欧美在线一区亚洲| 少妇的逼水好多| 精品久久久久久久久亚洲 | 午夜福利在线在线| 99riav亚洲国产免费| 啦啦啦韩国在线观看视频| 无遮挡黄片免费观看| 久久天躁狠狠躁夜夜2o2o| 日日夜夜操网爽| 亚洲美女黄片视频| 久久久色成人| 黄色配什么色好看| 丰满乱子伦码专区| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区 | 国产成人福利小说| 男人狂女人下面高潮的视频| 亚洲一区二区三区色噜噜| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女| 一个人免费在线观看电影| 中文字幕人妻熟人妻熟丝袜美| 成人毛片a级毛片在线播放| 最近中文字幕高清免费大全6 | 热99在线观看视频| 久久久国产成人免费| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 麻豆久久精品国产亚洲av| 偷拍熟女少妇极品色| 观看免费一级毛片| 色噜噜av男人的天堂激情| 亚洲欧美日韩高清专用| 欧美激情国产日韩精品一区| 国产成年人精品一区二区| 色视频www国产| 日韩中文字幕欧美一区二区| 亚洲av熟女| 免费在线观看日本一区| 国产伦精品一区二区三区视频9| 亚洲成人久久爱视频| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 老司机福利观看| 麻豆久久精品国产亚洲av| 日本爱情动作片www.在线观看 | 我的老师免费观看完整版| 天美传媒精品一区二区| 小蜜桃在线观看免费完整版高清| 美女高潮喷水抽搐中文字幕| 精品99又大又爽又粗少妇毛片 | 亚洲av电影不卡..在线观看| 22中文网久久字幕| 给我免费播放毛片高清在线观看| 国产精品1区2区在线观看.| 欧美性猛交黑人性爽| 亚洲美女搞黄在线观看 | 噜噜噜噜噜久久久久久91| 精品人妻一区二区三区麻豆 | 丰满的人妻完整版| 精品日产1卡2卡| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 色5月婷婷丁香| 91麻豆精品激情在线观看国产| 夜夜夜夜夜久久久久| 永久网站在线| 日韩高清综合在线| 白带黄色成豆腐渣| 国产精品女同一区二区软件 | 亚洲自偷自拍三级| 成人特级av手机在线观看| 久久久午夜欧美精品| 男女做爰动态图高潮gif福利片| 中文字幕熟女人妻在线| 无人区码免费观看不卡| 色综合婷婷激情| av中文乱码字幕在线| 国产精品精品国产色婷婷| 成年免费大片在线观看| 国产毛片a区久久久久| 国产免费av片在线观看野外av| 精品久久久久久久末码| 亚洲精品乱码久久久v下载方式| 国产精品1区2区在线观看.| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 日韩一区二区视频免费看| 国产精品福利在线免费观看| 69人妻影院| 毛片女人毛片| 亚洲三级黄色毛片| 22中文网久久字幕| 最近视频中文字幕2019在线8| 性欧美人与动物交配| 国产精品国产高清国产av| 日韩欧美一区二区三区在线观看| 在线观看午夜福利视频| 日韩欧美在线二视频| 免费无遮挡裸体视频| 精品一区二区免费观看| 美女cb高潮喷水在线观看| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 美女xxoo啪啪120秒动态图| xxxwww97欧美| 欧美日韩国产亚洲二区| 日韩精品有码人妻一区| 日韩欧美免费精品| 女的被弄到高潮叫床怎么办 | www日本黄色视频网| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 老熟妇仑乱视频hdxx| 天天躁日日操中文字幕| 熟妇人妻久久中文字幕3abv| 听说在线观看完整版免费高清| 国内精品久久久久精免费| 国产伦精品一区二区三区视频9| 又爽又黄无遮挡网站| 又粗又爽又猛毛片免费看| 俺也久久电影网| 免费看光身美女| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 波多野结衣高清作品| 色尼玛亚洲综合影院| 精品人妻视频免费看| 热99re8久久精品国产| 午夜福利在线观看免费完整高清在 | 国产真实伦视频高清在线观看 | 热99在线观看视频| 亚洲人成网站在线播| 亚洲美女黄片视频| 麻豆av噜噜一区二区三区| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 夜夜夜夜夜久久久久| 欧美激情久久久久久爽电影| 欧美zozozo另类| 国产午夜精品论理片| 成年女人永久免费观看视频| 国产成人a区在线观看| 色哟哟哟哟哟哟| 97超视频在线观看视频| 免费av不卡在线播放| 一区福利在线观看| 一本一本综合久久| 免费看光身美女| 日韩中字成人| 欧美在线一区亚洲| 女人被狂操c到高潮| 国产黄片美女视频| 国产男人的电影天堂91| 桃色一区二区三区在线观看| 欧美日韩瑟瑟在线播放| 88av欧美| 极品教师在线视频| 欧美bdsm另类| 亚洲av电影不卡..在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品国产成人久久av| 国产精品一区二区免费欧美| 欧美日韩精品成人综合77777| 乱系列少妇在线播放| 国产毛片a区久久久久| 国产三级中文精品| 欧洲精品卡2卡3卡4卡5卡区| 91狼人影院| 99热网站在线观看| 如何舔出高潮| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 国产亚洲精品综合一区在线观看| 日韩欧美一区二区三区在线观看| 国产精品嫩草影院av在线观看 | 欧美最新免费一区二区三区| 有码 亚洲区| 欧美另类亚洲清纯唯美| 国产精品一区二区性色av| 午夜亚洲福利在线播放| 麻豆一二三区av精品| 精品久久国产蜜桃| 色综合婷婷激情| a级毛片a级免费在线| 免费人成在线观看视频色| 亚洲av美国av| 欧美一级a爱片免费观看看| 美女高潮喷水抽搐中文字幕| 成人一区二区视频在线观看| 精品久久国产蜜桃| 免费av观看视频| 欧美一区二区国产精品久久精品| 在线播放无遮挡| 嫩草影院精品99| 在线观看舔阴道视频| 99九九线精品视频在线观看视频| 国产伦在线观看视频一区| 免费看美女性在线毛片视频| 欧美激情久久久久久爽电影| 精品久久久久久久末码| av.在线天堂| 国产精华一区二区三区| 1000部很黄的大片| 99热6这里只有精品| 22中文网久久字幕| 麻豆精品久久久久久蜜桃| 2021天堂中文幕一二区在线观| 日本黄色片子视频| 国产不卡一卡二| 亚洲男人的天堂狠狠| avwww免费| 别揉我奶头 嗯啊视频| 国产蜜桃级精品一区二区三区| 在线播放国产精品三级| 97超视频在线观看视频| 日日夜夜操网爽| 琪琪午夜伦伦电影理论片6080| 久久久久久久久中文| 九九热线精品视视频播放| 赤兔流量卡办理| 我要看日韩黄色一级片| 欧美精品啪啪一区二区三区| 一区二区三区免费毛片| 亚洲人成网站在线播放欧美日韩| 1024手机看黄色片| 亚洲专区中文字幕在线| 女人被狂操c到高潮| 日韩中字成人| 亚洲天堂国产精品一区在线| 国产av不卡久久| 成人欧美大片| 乱人视频在线观看| 亚洲国产精品合色在线| 性色avwww在线观看| 国产极品精品免费视频能看的| 亚洲四区av| 欧美一区二区精品小视频在线| 简卡轻食公司| 国产午夜精品久久久久久一区二区三区 | www.www免费av| 在线观看免费视频日本深夜| 春色校园在线视频观看| 成人av一区二区三区在线看| 热99re8久久精品国产| 国产三级在线视频| 国产精品久久电影中文字幕| 我的女老师完整版在线观看| 日韩一区二区视频免费看| 嫩草影院精品99| 不卡视频在线观看欧美| 亚洲精华国产精华精| 在线观看免费视频日本深夜| 成人性生交大片免费视频hd| 国产午夜精品久久久久久一区二区三区 | 人妻少妇偷人精品九色| 日韩精品青青久久久久久| 免费在线观看成人毛片| 国产一区二区三区av在线 | 别揉我奶头~嗯~啊~动态视频| 男人舔奶头视频| 欧美中文日本在线观看视频| 日本 av在线| 精品人妻一区二区三区麻豆 | 搡老妇女老女人老熟妇| 波野结衣二区三区在线| 欧美色视频一区免费| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av在线| 蜜桃亚洲精品一区二区三区| 日日撸夜夜添| 欧美潮喷喷水| 偷拍熟女少妇极品色| 免费无遮挡裸体视频| 午夜福利视频1000在线观看| 97碰自拍视频| 精品午夜福利视频在线观看一区| 精品乱码久久久久久99久播| 日韩欧美在线乱码| 久久6这里有精品| 亚洲av二区三区四区| 高清毛片免费观看视频网站| 欧美一区二区精品小视频在线| 色视频www国产| 日韩精品有码人妻一区| 免费一级毛片在线播放高清视频| 综合色av麻豆| 天天一区二区日本电影三级| 亚洲欧美清纯卡通| 波多野结衣高清无吗| 亚洲欧美精品综合久久99| 香蕉av资源在线| 日韩一区二区视频免费看| 亚洲av一区综合| 国语自产精品视频在线第100页| 亚洲精华国产精华精| 黄色日韩在线| 欧美成人a在线观看| 十八禁国产超污无遮挡网站| 白带黄色成豆腐渣| 日本免费a在线| 欧美日本视频| 成年女人永久免费观看视频| 国产精品一区二区免费欧美| 色综合亚洲欧美另类图片| 国产v大片淫在线免费观看| 亚洲精品在线观看二区| 久9热在线精品视频| av在线亚洲专区| 国产一区二区激情短视频| 国产在线精品亚洲第一网站| 丰满的人妻完整版| 99热6这里只有精品| 国产综合懂色| 少妇丰满av| 麻豆av噜噜一区二区三区| 日韩欧美国产一区二区入口| 在线免费十八禁| 一区二区三区激情视频| 精品久久久久久久末码| 在线观看66精品国产| 国产精品久久久久久av不卡| 精品福利观看| 嫩草影院入口| 黄色日韩在线| 一区二区三区免费毛片| 成年版毛片免费区| 狠狠狠狠99中文字幕| 国产免费男女视频| 国产精品三级大全| 看十八女毛片水多多多| 欧美日韩乱码在线| 国产精品人妻久久久久久| 精品久久久久久久久av| 啦啦啦观看免费观看视频高清| 97超视频在线观看视频| 又爽又黄a免费视频| www.色视频.com| 内地一区二区视频在线| 国产伦精品一区二区三区视频9| 国产v大片淫在线免费观看| 午夜福利18| 极品教师在线视频| 国产精品免费一区二区三区在线| 亚洲中文日韩欧美视频| 久久午夜福利片| 男人舔奶头视频| 亚洲不卡免费看| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区三区四区免费观看 | 国产精品精品国产色婷婷| 国产高清不卡午夜福利| x7x7x7水蜜桃| 九九久久精品国产亚洲av麻豆| 免费一级毛片在线播放高清视频| 日韩av在线大香蕉| 日本黄大片高清| 久久久色成人| 日韩一本色道免费dvd| 在线看三级毛片| 22中文网久久字幕| 18禁黄网站禁片午夜丰满| 亚洲第一区二区三区不卡| 精品一区二区三区视频在线| 三级国产精品欧美在线观看| 男人的好看免费观看在线视频| 在线观看66精品国产| 亚洲黑人精品在线| 国产白丝娇喘喷水9色精品| 97热精品久久久久久| 日韩中文字幕欧美一区二区| 午夜精品在线福利| 成人毛片a级毛片在线播放| 51国产日韩欧美| 自拍偷自拍亚洲精品老妇| 99久国产av精品| 成人国产麻豆网| 国产av一区在线观看免费| 88av欧美| 亚洲av中文字字幕乱码综合| 偷拍熟女少妇极品色| 成人国产综合亚洲| 最近最新免费中文字幕在线| 日韩中文字幕欧美一区二区| 最近在线观看免费完整版| 又粗又爽又猛毛片免费看| 国产亚洲欧美98| 欧美bdsm另类| 国产在视频线在精品| 亚洲熟妇中文字幕五十中出| 午夜免费男女啪啪视频观看 | 欧美xxxx黑人xx丫x性爽| 神马国产精品三级电影在线观看| 国产av一区在线观看免费| 国产精品自产拍在线观看55亚洲| 中文亚洲av片在线观看爽| 久久精品国产99精品国产亚洲性色| 亚洲av美国av| 国内精品一区二区在线观看| 在线免费十八禁| 国产伦精品一区二区三区视频9| 国产亚洲精品综合一区在线观看| 桃红色精品国产亚洲av| 日本色播在线视频| 国产精品野战在线观看| 欧美成人一区二区免费高清观看| 亚洲美女黄片视频| 精品国内亚洲2022精品成人| 亚洲五月天丁香| av天堂在线播放| 国产成人a区在线观看| 久久久久久久精品吃奶| 男人舔奶头视频| 男女那种视频在线观看| 日日干狠狠操夜夜爽| 久久久久久九九精品二区国产| 久久精品国产鲁丝片午夜精品 | 成年版毛片免费区| 国内精品久久久久久久电影| 美女高潮喷水抽搐中文字幕| 少妇高潮的动态图| 中文字幕精品亚洲无线码一区| 在线免费观看的www视频| 看免费成人av毛片| 在线免费观看不下载黄p国产 | 91狼人影院| 久久久国产成人免费| 久久99热这里只有精品18| 亚洲无线观看免费| 麻豆成人午夜福利视频| 日韩欧美一区二区三区在线观看| 欧美人与善性xxx| 麻豆国产av国片精品| www.色视频.com| 日本熟妇午夜| 能在线免费观看的黄片| 国产91精品成人一区二区三区| 亚洲最大成人中文| 非洲黑人性xxxx精品又粗又长| 亚洲自偷自拍三级| 国产黄片美女视频| xxxwww97欧美| 精品久久久久久久久av|