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

    東亞邊界層臭氧時空分布的數(shù)值模擬研究

    2014-05-13 03:03:42楊文夷陳煥盛王自發(fā)李健軍中國科學院大氣物理研究所大氣邊界層物理和大氣化學國家重點實驗室北京00029中國科學院大學北京00049中國環(huán)境監(jiān)測總站北京0002
    中國環(huán)境科學 2014年7期
    關(guān)鍵詞:東亞地區(qū)邊界層東亞

    楊文夷,李 杰,陳煥盛,王自發(fā),胡 波,宋 濤,李健軍(.中國科學院大氣物理研究所,大氣邊界層物理和大氣化學國家重點實驗室,北京 00029;2.中國科學院大學,北京 00049;.中國環(huán)境監(jiān)測總站,北京0002)

    東亞邊界層臭氧時空分布的數(shù)值模擬研究

    楊文夷1,2,李 杰1*,陳煥盛1,王自發(fā)1,胡 波1,宋 濤1,李健軍3(1.中國科學院大氣物理研究所,大氣邊界層物理和大氣化學國家重點實驗室,北京 100029;2.中國科學院大學,北京 100049;3.中國環(huán)境監(jiān)測總站,北京100012)

    利用嵌套網(wǎng)格空氣質(zhì)量預報模式系統(tǒng)(NAQPMS)對2010年東亞地區(qū)邊界層臭氧(O3)的時空分布進行了數(shù)值模擬,并評估了東亞邊界層光化學反應的活性.結(jié)果表明,NAQPMS模式與觀測結(jié)果較為一致,站點觀測與模擬的日均值(月均值)相關(guān)系數(shù)達到 0.56~0.91,模式能合理再現(xiàn)東亞地區(qū)地面O3的時空分布特征.東亞地區(qū)冬季邊界層O3低值區(qū)出現(xiàn)在中國東部;春季O3濃度增加,西北太平洋沿岸地區(qū)O3濃度達60μL/m3左右;夏季東亞中緯度35°N附近大陸地區(qū)O3由于強烈的光化學反應呈現(xiàn)出一濃度高值帶,濃度達60μL/m3以上;秋季東亞大部分地區(qū)O3濃度維持在40~45μL/m3左右.夏季中國京津冀和長江三角洲部分地區(qū)光化學凈生成率已超過30×10-9/d.

    嵌套網(wǎng)格空氣質(zhì)量模式;臭氧;數(shù)值模擬

    隨著城鎮(zhèn)化經(jīng)濟的快速發(fā)展,我國大氣污染呈現(xiàn)出復合型和區(qū)域型的新特點,以灰霾和光化學煙霧為代表的二次污染交替、耦合發(fā)生[1-2],對我國人體健康和農(nóng)作物生產(chǎn)帶來了潛在威脅[3-4].大氣氧化能力則是大氣復合污染的核心因素,其強弱不僅決定大氣一次污染物(如二氧化硫、氮氧化物和一氧化碳等)的清除速率,同時明顯影響大氣細顆粒物(PM2.5)主要組分(如硫酸鹽、硝酸鹽和二次有機顆粒物)的質(zhì)量濃度[5].

    對流層臭氧(O3)是重要的大氣氧化劑,同時也是OH自由基的主要前體物[6],了解其在我國的時空分布特征是準確評估我國大氣氧化能力的前提之一,數(shù)值模式是研究該問題的重要途徑之一.針對中國特別是京津冀、長江三角洲和珠江三角洲臭氧的時空分布及其來源已做了大量數(shù)值模擬研究[7-9],然而這些研究多集中在城市夏季臭氧的形成機理以及單一區(qū)域的個例模擬,對我國邊界層臭氧長期時空分布,特別是我國灰霾頻發(fā)的秋冬兩季的模擬還遠遠不足,使得對我國灰霾期間大氣二次顆粒物生消過程的認識嚴重受限[10].

    缺乏臭氧長期數(shù)值模擬也是我國臭氧生態(tài)環(huán)境效應評估的限制因素之一.研究表明,O3對植物的影響主要由于長期累積效應引起,高濃度的O3暴露會影響農(nóng)作物的生長發(fā)育和產(chǎn)量[11].

    本文利用中國科學院大氣物理研究所自主研發(fā)的嵌套網(wǎng)格空氣質(zhì)量預報模式系統(tǒng)(NAQPMS)對東亞地區(qū)2010年O3濃度進行了連續(xù)一年模擬,通過與觀測數(shù)據(jù)進行對比驗證了模式對 O3長期變化特征的模擬能力,并探討了東亞地區(qū)邊界層O3及其光化學反應速率的時空分布特征.

    1 模式設(shè)置與觀測數(shù)據(jù)

    1.1 模式介紹

    NAQPMS是由中國科學院大氣物理研究所自主研發(fā)的多物種、多尺度空氣質(zhì)量模擬系統(tǒng)[12],已被用于區(qū)域和城市尺度的沙塵輸送、酸雨、臭氧等污染物的形成機制和跨界輸送等科學研究[13];同時在業(yè)務管理方面也得到了廣泛應用,已在全國多個省市(北京、上海、廣東等)投入業(yè)務化運行.

    NAQPMS模式由氣象模塊和嵌套化學輸送模塊構(gòu)成,氣象模塊采用由美國環(huán)境預測中心(NCEP)和美國國家大氣研究中心(NCAR)等機構(gòu)研制的中尺度氣象模式 WRFv3.0;嵌套化學輸送模塊由中國科學院大氣物理研究所自主開發(fā),考慮了污染物的平流、擴散和干、濕沉降等過程,氣相化學機制采用CBM-Z機制.CBM-Z是基于CBM-IV發(fā)展的按結(jié)構(gòu)分類的物種更多、反應更完善、適用性更廣的一個歸納化學機理,反應物種增加至71個,化學反應增加至176個,模式系統(tǒng)的詳細介紹參見文獻[14].

    1.2 模式設(shè)置

    模擬區(qū)域包括了整個東亞地區(qū)(圖1),水平網(wǎng)格分辨率為 80km,網(wǎng)格數(shù)為 97(東西)×77(南北).垂直方向模式采用地形追隨坐標系,從地面至20km高空不均勻地分為20層,其中最下面7層位于距地面1~2km范圍內(nèi).

    模擬時段為2010年1月1日~2010年12月31日,初始條件與邊界條件來源于全球大氣化學輸送模式MOZARTv2.4的模擬結(jié)果.為減少因初始條件擾動對模擬結(jié)果的影響,NAQPMS的啟動時間設(shè)定為2009年12月1日.

    人為排放源取自清華大學提供的基準年為2009年的排放清單,其分辨率為 0.25°×0.25°,排放物種包括二氧化硫(SO2)、氮氧化物(NOχ)、揮發(fā)性有機物(VOCs)、黑碳(BC)、有機物(OM)和直接排放的顆粒物等[15].生物質(zhì)燃燒取自于曹國良等[16]的研究,并利用MODIS衛(wèi)星的火點資料,將污染物排放量細化為 0.25°×0.25°的空間分辨率,并將年均值細化為日均值.此外,碳氫化合物的自然源排放取自全球散發(fā)物質(zhì)名錄行動(GEIA)[17].圖1中給出了研究區(qū)域內(nèi)NOχ的排放速率.

    圖1 NAQPMS模擬區(qū)域、觀測站點分布以及NOχ的排放速率[μg/(m2·s)]Fig.1 Distribution of stations in the study and modeling domain of NAQPMS and the NOχemission rate [μg/(m2·s)]

    1.3 觀測數(shù)據(jù)

    為了檢驗模擬結(jié)果的合理性,將模式輸出的O3濃度與站點觀測結(jié)果進行對比.觀測站點的地理位置和其他相關(guān)信息如圖1和表1所示.觀測資料主要來源于:(1)中國科學院大氣物理研究所建立的中國科學院京津冀大氣環(huán)境監(jiān)測網(wǎng)[18];(2)中國環(huán)境監(jiān)測總站建立的大氣顆粒物全國觀測網(wǎng)絡(luò);(3)東亞酸雨觀測網(wǎng)(EANET)(http://www.eanet. cc/product.html);(4)其他數(shù)據(jù)來源,包括相關(guān)參考文獻.需要指出的是,個別站點的觀測年份并非2010年,以下針對這些站點的比較建立在月均值的基礎(chǔ)上,對這些站點數(shù)據(jù)的驗證用于評估模式對于季節(jié)變化規(guī)律的模擬能力.

    表1 觀測站點信息Table 1 The information of stations

    2 模擬結(jié)果分析

    2.1 模式結(jié)果驗證

    圖2展示了2010年東亞地區(qū)觀測和模擬得到的地面O3年均濃度.模擬值選取與觀測臺站同一海拔高度的結(jié)果.如圖2所示,模擬結(jié)果與觀測較為一致,基本再現(xiàn)了 O3的空間分布特征,如中國華北和珠江三角洲地區(qū)的低值區(qū)(20~30μL/m3)和日本北海道的高值區(qū)(30~40μL/m3).定量上講,模擬和觀測結(jié)果的相關(guān)系數(shù)達到 0.94,所有臺站的模擬值均在觀測值2倍范圍內(nèi).

    圖2 2010年東亞地區(qū)地面O3年均濃度水平分布以及觀測與模擬年均濃度散點圖Fig.2 The shaded figure shows the distribution of annual mixing ratios of surface O3and the scatter plot shows the relationship between observed and simulated annual O3

    圖3 觀測站點地面O3月均濃度時間序列Fig.3 Monthly mean O3mixing ratios with standard deviation of observation and model

    圖3展示了所有站點月均O3濃度及其標準偏差的觀測和模擬結(jié)果.如圖3所示,模式能較好地模擬不同地區(qū)O3濃度的季節(jié)變化特征.如我國京津冀地區(qū)的地面臭氧呈單峰型分布特征,其最大值出現(xiàn)在 6月,最小值在冬季;華南地區(qū)(如廣州)的最大值出現(xiàn)在秋季(10~11月),東北地區(qū)(龍鳳山)則出現(xiàn)在春季.與上述地區(qū)不同,日本和西太平洋地區(qū)的臭氧呈雙峰型分布特征,即春季最高、夏季最低,秋季出現(xiàn)一個次峰.模式均很好再現(xiàn)這些特征.需要指出的是,模式雖然合理再現(xiàn)了承德臭氧的季節(jié)變化特征,但是模擬濃度整體偏高,這可能與模式低估承德城區(qū)NOχ對O3的滴定效應有關(guān).

    表2給出了各觀測站點的相關(guān)統(tǒng)計參數(shù),如表 2所示,大部分站點觀測與模擬的相關(guān)系數(shù)達到 0.56~0.91,所有站點均通過了 1%的顯著性檢驗,證明模擬值與觀測值之間的相關(guān)系數(shù)是顯著的.承德站模擬結(jié)果普遍高于觀測值,其NMB為57%,其他站點的NMB在-9%~33%之間.

    圖4展示了不同季節(jié)O3日均(或月均)濃度的觀測與模擬對比散點圖,不同季節(jié)O3的模擬值大部分落在觀測值的 2倍范圍內(nèi),日均值線性回歸得到的直線斜率在 0.9左右,而月均值線性回歸斜率在0.6~0.9范圍內(nèi).不同季節(jié)觀測與模擬濃度的相關(guān)系數(shù)在0.65~0.85之間,均通過了顯著性水平為1%的t檢驗.

    圖4 不同季節(jié)O3日均濃度(a)、月均濃度(b)觀測與模擬散點圖Fig.4 Scatter plot showing the relationship between observed and simulated daily averaged O3(a) and monthly average O3(b) in different seasons

    表2 O3濃度日均值模擬結(jié)果統(tǒng)計參數(shù)Table 2 Statistical parameters of model performance with daily concentration of O3

    總體而言,NAQPMS能較好地模擬東亞地區(qū)地面O3濃度,并能基本再現(xiàn)東亞不同地區(qū)O3濃度的季節(jié)變化特征等.

    2.2 東亞地區(qū)邊界層O3時空分布特征

    圖5展示了2010年東亞地區(qū)邊界層O3濃度及輸送通量的季節(jié)分布.如圖5所示,東亞邊界層O3呈現(xiàn)明顯的季節(jié)變化特征.在冬季,東亞大部分地區(qū) O3濃度維持在 20~50μL/m3,低值區(qū)出現(xiàn)在中國東部(30~40μL/m3),山東、河南和河北等地區(qū)的濃度甚至低至 30μL/m3以下.西太平洋、蒙古和東西伯利亞等地區(qū)O3濃度較中國東部略高,達到40~50μL/m3.這種分布特征與東亞的季風氣候和污染源分布有關(guān).冬季,整個東亞地區(qū)受到蒙古冷高壓控制,30°N以北地區(qū)盛行西北季風,以南盛行東北季風[23].這些氣團將高緯地區(qū)背景臭氧輸送至整個東亞,使得該地區(qū)臭氧濃度較為均一.然而,中國東部特別是華北地區(qū)NOχ排放量較大(圖 1),在冬季太陽輻射較低的情況下,過量的NOχ與O3發(fā)生反應大量消耗O3,造成該地區(qū)O3濃度較低.受到 O3的分布特征和氣象條件的影響,O3高水平輸送通道位于西伯利亞-蒙古-東北-韓國-日本一線,其輸送通量達到 2×10-5mol/ (m2·s).

    與冬季相比,春季東亞北部地區(qū)(25°N以北)的O3濃度大幅度增加,在西伯利亞和日本、韓國等地區(qū)甚至達到 60μL/m3以上,形成東亞的兩個高值中心.東亞南部(臺灣以南)地區(qū)的O3濃度則降低5~10μL/m3,維持在40μL/m3以下.春季東亞北部地區(qū)太陽輻射增強,中國華北、華東,韓國以及日本地區(qū)排放的高強度O3前體物在邊界層內(nèi)的光化學反應明顯增強,生成大量 O3,成為東亞地區(qū)O3高值區(qū)之一;同時控制東亞大陸的蒙古高壓進一步減弱,亞洲大陸污染物向西北太平洋的O3輸送通道稍向北偏移,移至35~50°N一帶,受上游地區(qū)生成的高體積分數(shù)O3影響,黃海、東海以及日本海海域的O3濃度維持在較高值.中國南部受到來自印度洋和中國南海的海洋性氣團影響,低濃度的O3由海洋輸送至中國南部,造成東亞南部地區(qū)O3濃度明顯降低.

    夏季東亞大陸35°N中緯度地區(qū)附近形成一臭氧高值帶,中國京津冀地區(qū)、山東、山西以及河南省的邊界層O3濃度可達60μL/m3左右.太平洋及沿岸地區(qū)O3濃度大幅降低,洋面上濃度維持在 20μL/m3以下,中國臺灣和日本部分地區(qū)濃度降至 30μL/m3左右.夏季北半球太陽輻射進一步加強,東亞地區(qū)充足的O3前體物在合適的條件下發(fā)生光化學反應生成大量O3,污染物局地累積造成35°N大陸地區(qū)形成O3高值區(qū).中國南部和日本等地在東亞夏季風影響下盛行偏東或偏南氣流,來自海洋的清潔氣團使得這些地區(qū)邊界層O3濃度較低,維持在30μL/m3以下.

    圖5 2010年東亞邊界層O3濃度(μL/m3)及水平輸送通量[1×10-5mol/(m2·s)]的季節(jié)分布Fig.5 Seasonal mean ozone mixing ratios (μL/m3) and ozone horizontal fluxes [1×10-5mol/(m2·s)] in the boundary layer in 2010

    秋季東亞大部分地區(qū)邊界層O3濃度維持在40~45μL/m3左右.此時光化學反應強度隨著太陽輻射的減弱而降低,夏季35°N大陸地區(qū)O3高值區(qū)濃度大幅度下降,中國華北等地區(qū)邊界層O3濃度在秋季降低至35μL/m3左右,成為低值區(qū).由于東亞季風由夏季風向冬季風轉(zhuǎn)變,中國南部地區(qū)、日本以東的西北太平洋地區(qū)受到大陸污染氣團的影響,O3濃度回升至40μL/m3左右,并在整個秋冬季都維持在較高濃度.

    2.3 臭氧生成速率的時空分布特征

    圖6展示了2010年東亞地區(qū)邊界層O3光化學凈生成率(ONP,單位時間O3凈生成量)的季節(jié)分布.如圖6所示,ONP呈現(xiàn)明顯的季節(jié)變化特征.冬季,太陽輻射較弱,整個東亞地區(qū)ONP維持在較低水平,基本在 10×10-9/d以下.春季北半球太陽輻射增強,東亞大陸地區(qū)人為污染源排放較強的地方ONP增加5×10-9/d左右,中國35°N以南地區(qū)維持在 10×10-9/d以上.夏季光化學反應強烈,東亞大陸大部分地區(qū)ONP較春季有明顯增加,中國華北、華東等地基本維持在 20×10-9/d以上,中國京津冀和長三角部分地區(qū)ONP甚至超過30×10-9/d.秋季整個東亞地區(qū) ONP的分布與春季較為接近.總體而言,ONP與排放源分布(圖1)有較好的一致性,說明污染源排放是影響ONP的主要因子之一.

    圖6 2010年東亞邊界層O3光化學凈生成率(10-9/d)的季節(jié)分布Fig.6 Seasonal mean net ozone photochemical production (10-9/d) in the boundary layer in 2010

    隨著污染物排放強度的演變,東亞地區(qū) ONP強度也呈現(xiàn)明顯的變化趨勢.中國及韓國部分地區(qū)夏季ONP較20世紀90年代增長(5~10)×10-9/d以上

    [22,24],中國東部地區(qū)增長幅度最大,河北北部、河南南部以及江蘇、安徽、浙江三省交界處,增長幅度達到(15~20)×10-9/d(約為100%).與之相反的是,夏季日本東京地區(qū) ONP呈現(xiàn)出下降趨勢,降低至15×10-9/d以下,這與污染物的演變趨勢一致[25-26].

    大量研究表明,歐洲特別是其中部地區(qū)目前嚴重的臭氧污染,臭氧超標天數(shù)達到10%~20%[27].基于觀測數(shù)據(jù),Connor等[28]計算2000年7~8月歐洲中部地區(qū)白天邊界層 ONP為(13~45.5)×10-9/d.與之相比,我國東部大部分地區(qū)夏季全天邊界層ONP歐洲中部地區(qū)相當(10×10-9/d以上),個別地區(qū)甚至超過歐洲中部(超過40×10-9/d).

    3 結(jié)論

    3.1 利用NAQPMS對東亞地區(qū)2010年O3濃度進行了模擬,為檢驗模擬質(zhì)量,將模擬結(jié)果與觀測數(shù)據(jù)進行對比,27個觀測站點觀測與模擬的O3濃度日均(月均)值均有較好相關(guān)性.絕大部分站點的相關(guān)系數(shù)在0.6~0.9之間,所有站點均通過了顯著性水平為1%的t檢驗,表明模式可以較好地再現(xiàn)東亞地區(qū)地面O3的濃度分布和季節(jié)變化特征.

    3.2 東亞地區(qū)不同季節(jié)邊界層 O3濃度呈現(xiàn)出各自的空間分布特征,主要影響因素有對流層光化學反應生成,東亞季風環(huán)流和平流層向?qū)α鲗拥妮斔偷?冬季邊界層O3低值區(qū)位于中國東部地區(qū),山東、河南和河北等地區(qū)的濃度低至30μL/m3;春季東亞25°N以北地區(qū)O3濃度大幅度增加,西伯利亞和日本、韓國等地O3濃度高達60μL/m3左右;夏季東亞中緯度 35°N附近大陸地區(qū)邊界層O3濃度由于強烈的光化學反應呈現(xiàn)出一臭氧高值帶;秋季中國南部地區(qū)、日本以東的西北太平洋地區(qū)由夏季風向冬季風轉(zhuǎn)變,在大陸污染氣團的影響下邊界層O3濃度回升至40μL/m3左右.

    3.3 東亞地區(qū)人為活動強烈的地區(qū) ONP較大,因此污染源排放是影響ONP的主要因子.冬季整個東亞地區(qū) ONP維持在 10×10-9/d以下.春季ONP有所增加,中國 35°N以南地區(qū)維持在10× 10-9/d以上.夏季光化學反應強烈,人為污染物排放較大地區(qū)ONP大幅度增加,中國京津冀和長三角部分地區(qū)ONP甚至超過30×10-9/d.秋季ONP的分布與春季較為接近.

    [1] 張遠航,邵可聲,唐孝炎,等.中國城市光化學煙霧污染研究 [J].北京大學學報(自然科學版), 1998,34(2/3):392-400.

    [2] Wang X, Zhang Y, Hu Y, et al. Process analysis and sensitivity study of regional ozone formation over the Pearl River Delta, China, during the PRIDE-PRD2004 campaign using the Community Multiscale Air Quality modeling system [J]. Atmospheric Chemistry and Physics, 2010,10:4423-4437.

    [3] Brauer M, Brook J R. Ozone personal exposures and health effects for selected groups residing in the Fraser valley [J]. Atmospheric Environment, 1997,31(14):2113-2121.

    [4] 劉宏舉,鄭有飛,吳榮軍,等.地表臭氧濃度增加對南京地區(qū)冬小麥生長和產(chǎn)量的影響 [J]. 中國農(nóng)業(yè)氣象, 2009,30(2):195-200.

    [5] 程艷麗,王雪松,劉兆榮,等.大氣氧化性定量表征方法的建立及其應用 [J]. 中國科學B輯, 2008,38(10):938-946.

    [6] 陸克定,張遠航.HOx自由基的實地測量及其化學機制解析 [J].化學機制, 2010,22(2/3):500-514.

    [7] 馬志強,王躍思,張小玲,等.北京城區(qū)與下游地區(qū)臭氧對比研究[J]. 環(huán)境科學, 2011,32(4):924-929.

    [8] 耿福海,毛曉琴,鐵學熙,等.2006-2008年上海地區(qū)臭氧污染特征與評價指標研究 [J]. 熱帶氣象學報, 2010,26(5):584-590.

    [9] 鄧 濤,吳 兌,鄧雪嬌,等.珠江三角洲一次典型復合型污染過程的模擬研究 [J]. 中國環(huán)境科學, 2012,32(2):193-199.

    [10] 賀 泓,王新明,王躍思,等.大氣灰霾追因與控制 [J]. 中國科學院院刊, 2013,28(3):344-352.

    [11] Calatayud A, Iglesias D J, Talon M, et al. Effects of 2-month ozone exposure in spinach leaves on photosynthesis, antioxidant systems and lipid peroxidation [J]. Plant Physiology and Biochemistry, 2003,41:839-845.

    [12] 王自發(fā),謝付瑩,王喜全,等.嵌套網(wǎng)格空氣質(zhì)量預報模式系統(tǒng)的發(fā)展與應用 [J]. 大氣科學2006,30(5):778-790.

    [13] Wang Zifa, Li Jie, Wang Xiquan, et al. Modeling of regional high ozone episode observed at two mountain sites (Mt. Tai and Huang) in East China [J]. Journal of Atmospheric Chemistry, 2006,55(3): 253-272.

    [14] Li J, Wang Z F, Zhuang G, et al. Mixing of Asian mineral dust with anthropogenic pollutants over East Asia: a model cast study of a super-duststorm in March 2010 [J]. Atmospheric Chemistry and Physics, 2012,12:7591-7607.

    [15] Wang S X, Xing J, Chatani S, et al. Verification of anthropogenic emissions of China by satellite and ground observations [J]. Atmospheric Environment, 2011,45:6347-6358.

    [16] 曹國良,張小曳,王 丹,等.中國大陸生物質(zhì)燃燒排放的污染物清單 [J]. 中國環(huán)境科學, 2005,25(4):389-393.

    [17] Guenther A, Hewitt N, Erickson D, et al. A global model ofnatural volatile organic compound emissions [J]. Journal of Geophysics Research, 1995,100(D5):8873-8892.

    [18] 辛金元,王躍思,唐貴謙,等.2008年奧運期間北京及周邊地區(qū)污染物消減變化 [J]. 科學通報, 2010,55(15):1510-1519.

    [19] 徐曉斌,丁國安,李興生,等.龍鳳山大氣近地層O3濃度變化及與其它因素的關(guān)系 [J]. 氣象學報, 1998,(5):49-50,52-61.

    [20] Ge B Z, Xu X B, Lin W L, et al. Impact of the regional transport of urban Beijing pollutants on downwind areas in summer: ozone production efficiency analysis [J]. Tellus. B, 2012,64:17348.

    [21] Wang T, Wong C H, Cheung T F, et al. Relationships of trace gases and aerosols and the emission characteristics at Lin'an, a rural site in eastern China, during spring 2001 [J]. Journal of Geophysical Research-Atmospheres, 2004,109(D19), doi: 10.1029/2003JD004119.

    [22] Li Jie, Wang Zifa, Akimoto H, et al. Modeling study of ozone seasonal cycle in lower troposphere over east Aisa [J]. Journal of Geophysical Research, 2007,112(D22s25), doi: 10.1029/ 2006JD008209.

    [23] 朱 彬.東亞太平洋地區(qū)近地面臭氧的季節(jié)和年際變化特征及其與東亞季風的關(guān)系 [J]. 大氣科學學報, 2012,35(5):513-523.

    [24] Mauzerall D L, Narita D, Akimoto H, et al. Seasonal characteristics of tropospheric ozone production and mixing ratios over East Asia: A global three-dimensional chemical transport model analysis [J]. Journal of Geophysical Research, 2000, 105(D14):17895-17910.

    [25] Ohara T, Akimoto H, Kurokawa J, et al. An Asian emission inventory of anthropogenic emission sources for the period 1980-2020 [J]. Atmospheric Chemistry and Physics, 2007,7: 4419-4444.

    [26] Kurokawa J, Ohara T, Morikawa T, et al. Emissions of air pollutants and greenhouse gases over Asian regions during 2000-2008: Regional Emission inventory in Asia (REAS) version 2 [J]. Atmospheric Chemistry and Physics Discussions, 2013,13: 10049-10123.

    [27] Varotsos K V, Tombrou M, Giannakopoulos C. Statistical estimations of the number of future ozone exceedances due to climate change in Europe [J]. Journal of Geophysical Research: Atmospheres, 2013,118:6080-6099.

    [28] O’Connor F M, Law K S, Pyle J A, et al. Tropospheric ozone budget: regional and global calculations [J]. Atmospheric Chemistry and Physics Discussions, 2004,4:991-1036.

    致謝:感謝清華大學王書肖研究員提供2009年中國大陸人為排放源和日本全球變動研究所Yugo Kanaya博士提供日本Fukue臭氧的觀測結(jié)果.

    Modeling analysis of boundary layer ozone distributions over East Asia.

    YANG Wen-yi1,2, LI Jie1*, CHEN

    Huan-sheng1, WANG Zi-fa1, HU Bo1, SONG Tao1, LI Jian-jun3(1.The state Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China;2.University of Chinese Academy of Sciences, Beijing 100049, China;3.China National Environmental Monitoring Centre, Beijing 100012, China). China Environmental Science, 2014,34(7):1633~1641

    The nested air quality prediction model system (NAQPMS) was applied to investigate the temporal and spatial variations of boundary layer ozone over China in 2010, and the photochemical activity in the boundary layer was assessed. The model performed well with the observed daily (monthly) mixing ratio of ozone, and the correlation coefficients ranging from 0.56 to 0.91. Boundary O3distribution over East Asia showed a significant seasonal variability with the meteorological and photochemical condition. In winter, the seasonal averaged O3concentrations over East Asia were among 20~50μL/m3with low O3concentrations cover the East China. The major outflow transport pathway for O3to the western Pacific in spring was at 25~40°N, where the ozone mixing ratios in the lower troposphere might reach up to 60μL/m3. In summer, high ozone (about 60μL/m3) appeared in the East China around 35°N, where pollutant emissions were strong and ozone net production was high. In autumn, the distribution of O3was similar to that in spring, but O3concentrations were around 40~45μL/m3lower than those in spring. The net photochemical production of ozone in the North China and the Yangtze River Delta has exceeded 30×10-9/d in summertime.

    nested air quality prediction model system;O3;numerical simulation

    X515

    A

    1000-6923(2014)07-1633-09

    楊文夷(1989-),女,江蘇蘇州人,中國科學院大氣物理研究所博士研究生,主要從事大氣物理與大氣環(huán)境方向研究.發(fā)表論文4篇.

    《中國環(huán)境科學》2011年度引證指標

    《中國環(huán)境科學》編輯部

    2013-09-15

    中國科學院戰(zhàn)略性先導項目(XDB05030101,XDA05100501);環(huán)保行業(yè)公益性科研專項(201009002);國家自然基金(40805051)

    * 責任作者, 副研究員, lijie8074@mail.iap.ac.cn

    根據(jù)《2012年版中國科技期刊引證報告(核心版)》,《中國環(huán)境科學》2011年度引證指標繼續(xù)位居環(huán)境科學技術(shù)、安全科學技術(shù)類科技期刊前列,核心影響因子1.523,學科排名第1,綜合評價總分79.2,學科排名第2;在被統(tǒng)計的1998種核心期刊中影響因子列第18位,綜合評價總分列第52位.《中國科技期刊引證報告》每年由中國科學技術(shù)信息研究所編制,統(tǒng)計結(jié)果被科技管理部門和學術(shù)界廣泛采用.

    猜你喜歡
    東亞地區(qū)邊界層東亞
    “東亞漢詩史(多卷本)”簡介
    東疆學刊(2021年4期)2021-02-12 01:50:18
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    我校東亞研究院一行應邀訪問韓國東亞大學
    如何看待我國的高投資率
    會“珰”洞鑒以昭然——古代東亞地區(qū)冠飾“珰”之探究
    未定的東亞多邊體制與能動的中國東亞政策
    教學與研究(2016年1期)2016-02-28 01:59:27
    鼎盛與危機:明清東亞宗藩體系嬗變
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    非特征邊界的MHD方程的邊界層
    淺議東亞地區(qū)人民幣國際化的優(yōu)勢
    国产在线视频一区二区| 老司机福利观看| 青春草视频在线免费观看| www日本在线高清视频| 美女主播在线视频| av网站免费在线观看视频| 岛国毛片在线播放| 老司机午夜福利在线观看视频 | 狂野欧美激情性xxxx| 日韩三级视频一区二区三区| 老司机影院毛片| 久久精品国产亚洲av高清一级| 男女边摸边吃奶| 香蕉国产在线看| 中文精品一卡2卡3卡4更新| 国产野战对白在线观看| 两人在一起打扑克的视频| 久久人人97超碰香蕉20202| 在线av久久热| 精品国产一区二区三区久久久樱花| 超碰97精品在线观看| 久久精品久久久久久噜噜老黄| 搡老岳熟女国产| 午夜91福利影院| 欧美日韩视频精品一区| 亚洲精品自拍成人| 男女免费视频国产| 亚洲精品久久午夜乱码| 一级片免费观看大全| 亚洲精品成人av观看孕妇| 在线天堂中文资源库| 亚洲国产日韩一区二区| 99久久综合免费| 久久精品国产亚洲av高清一级| 国产视频一区二区在线看| 久久久久国产一级毛片高清牌| 欧美黄色片欧美黄色片| 亚洲专区字幕在线| 一进一出抽搐动态| 日本精品一区二区三区蜜桃| 免费日韩欧美在线观看| 亚洲国产av新网站| 国产片内射在线| 国产黄频视频在线观看| 老司机深夜福利视频在线观看 | 美女高潮喷水抽搐中文字幕| 国产成人a∨麻豆精品| 香蕉丝袜av| 亚洲成人免费电影在线观看| 一级毛片电影观看| 人成视频在线观看免费观看| 欧美+亚洲+日韩+国产| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品日韩在线中文字幕| 欧美 日韩 精品 国产| 日韩有码中文字幕| 亚洲精品成人av观看孕妇| 婷婷丁香在线五月| 精品国产一区二区三区久久久樱花| 在线亚洲精品国产二区图片欧美| 国产97色在线日韩免费| 国产主播在线观看一区二区| 亚洲欧美日韩高清在线视频 | 免费观看av网站的网址| 男女床上黄色一级片免费看| 久久国产精品影院| 国产精品一区二区在线不卡| 亚洲第一av免费看| 国产精品 国内视频| 波多野结衣av一区二区av| 伦理电影免费视频| 久久人妻熟女aⅴ| 午夜两性在线视频| 成年美女黄网站色视频大全免费| 午夜精品国产一区二区电影| 超色免费av| 色老头精品视频在线观看| 99久久国产精品久久久| 无遮挡黄片免费观看| 黑人欧美特级aaaaaa片| 又黄又粗又硬又大视频| 两人在一起打扑克的视频| 中文字幕人妻丝袜一区二区| 欧美97在线视频| tube8黄色片| 无遮挡黄片免费观看| 国产真人三级小视频在线观看| 国产欧美日韩一区二区三区在线| 成人av一区二区三区在线看 | 国产成人精品久久二区二区91| av线在线观看网站| 国产成人欧美| 国产精品一区二区精品视频观看| 欧美亚洲日本最大视频资源| 国产精品免费视频内射| 精品高清国产在线一区| 精品卡一卡二卡四卡免费| 搡老熟女国产l中国老女人| 正在播放国产对白刺激| 国产一卡二卡三卡精品| 亚洲第一青青草原| 国产欧美日韩一区二区三 | 三上悠亚av全集在线观看| 国产精品香港三级国产av潘金莲| 国产高清视频在线播放一区 | 母亲3免费完整高清在线观看| 亚洲欧美色中文字幕在线| 亚洲欧美色中文字幕在线| 另类亚洲欧美激情| a级片在线免费高清观看视频| 精品视频人人做人人爽| 99国产精品免费福利视频| 97在线人人人人妻| 亚洲欧美色中文字幕在线| 1024香蕉在线观看| 99国产精品一区二区蜜桃av | 亚洲综合色网址| 久久久久精品国产欧美久久久 | 高清欧美精品videossex| 咕卡用的链子| 男人爽女人下面视频在线观看| 在线观看免费高清a一片| 咕卡用的链子| 国产精品免费视频内射| 午夜视频精品福利| 亚洲国产毛片av蜜桃av| 两个人免费观看高清视频| 日韩电影二区| 考比视频在线观看| 欧美午夜高清在线| 黄色视频在线播放观看不卡| 欧美久久黑人一区二区| 91大片在线观看| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| svipshipincom国产片| 制服人妻中文乱码| 一级片'在线观看视频| 久久久久久亚洲精品国产蜜桃av| 麻豆国产av国片精品| 国产精品久久久久久人妻精品电影 | 婷婷色av中文字幕| 三上悠亚av全集在线观看| 亚洲欧美色中文字幕在线| 嫩草影视91久久| 午夜免费成人在线视频| 91字幕亚洲| av线在线观看网站| 久久精品熟女亚洲av麻豆精品| 精品亚洲成国产av| 老熟妇乱子伦视频在线观看 | 18禁国产床啪视频网站| 十八禁网站免费在线| 在线天堂中文资源库| 男女免费视频国产| 少妇的丰满在线观看| 多毛熟女@视频| 亚洲第一欧美日韩一区二区三区 | 18在线观看网站| 高清黄色对白视频在线免费看| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久人妻精品电影 | 国产在线观看jvid| 十八禁网站网址无遮挡| 人人妻,人人澡人人爽秒播| 一边摸一边做爽爽视频免费| 日本黄色日本黄色录像| 91成人精品电影| 精品国产一区二区三区四区第35| 午夜福利一区二区在线看| 日本91视频免费播放| 中文字幕人妻丝袜制服| 最近最新中文字幕大全免费视频| 成年人黄色毛片网站| 国产成人免费观看mmmm| 国产男人的电影天堂91| 欧美另类亚洲清纯唯美| 黄色a级毛片大全视频| 国产精品1区2区在线观看. | 欧美+亚洲+日韩+国产| 欧美+亚洲+日韩+国产| 久久精品国产a三级三级三级| 色婷婷av一区二区三区视频| 久久ye,这里只有精品| 老司机深夜福利视频在线观看 | 精品久久久久久久毛片微露脸 | 亚洲视频免费观看视频| 国产精品国产av在线观看| 国产亚洲欧美在线一区二区| 欧美激情高清一区二区三区| 天天添夜夜摸| 五月开心婷婷网| 婷婷丁香在线五月| 国产成人免费观看mmmm| 国产亚洲午夜精品一区二区久久| 午夜日韩欧美国产| 999精品在线视频| 99精品欧美一区二区三区四区| a级毛片在线看网站| 一级片'在线观看视频| 无限看片的www在线观看| 成人国产一区最新在线观看| 99国产精品免费福利视频| 国产精品偷伦视频观看了| 国产片内射在线| av免费在线观看网站| 美女中出高潮动态图| 天堂中文最新版在线下载| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕av电影在线播放| 老司机在亚洲福利影院| 女人被躁到高潮嗷嗷叫费观| av片东京热男人的天堂| 精品欧美一区二区三区在线| 亚洲成人手机| 日韩中文字幕欧美一区二区| 国产精品久久久人人做人人爽| 亚洲一码二码三码区别大吗| 久久九九热精品免费| 中文精品一卡2卡3卡4更新| 男男h啪啪无遮挡| 99国产精品一区二区三区| 好男人电影高清在线观看| 亚洲成av片中文字幕在线观看| 免费高清在线观看视频在线观看| 久久综合国产亚洲精品| 狂野欧美激情性bbbbbb| 天天躁狠狠躁夜夜躁狠狠躁| 国产有黄有色有爽视频| 亚洲国产av新网站| 国产免费一区二区三区四区乱码| 欧美日韩亚洲高清精品| 欧美日韩中文字幕国产精品一区二区三区 | 色综合欧美亚洲国产小说| a级毛片黄视频| 少妇精品久久久久久久| 欧美日韩成人在线一区二区| 不卡av一区二区三区| 十八禁网站免费在线| 欧美一级毛片孕妇| 成人av一区二区三区在线看 | 亚洲精品在线美女| 天天操日日干夜夜撸| 日本av手机在线免费观看| 伊人久久大香线蕉亚洲五| 久久天躁狠狠躁夜夜2o2o| 王馨瑶露胸无遮挡在线观看| 国产在线视频一区二区| 一级毛片电影观看| 国产xxxxx性猛交| 国产一区二区三区av在线| 老司机亚洲免费影院| 国产在线视频一区二区| 91大片在线观看| 精品亚洲成a人片在线观看| 久久性视频一级片| 久久久久国产精品人妻一区二区| 在线看a的网站| 国产成人av激情在线播放| 丁香六月天网| 午夜老司机福利片| 亚洲专区国产一区二区| 成年人免费黄色播放视频| 一本—道久久a久久精品蜜桃钙片| 建设人人有责人人尽责人人享有的| 国产精品av久久久久免费| 国产成人av激情在线播放| 青春草视频在线免费观看| 亚洲熟女精品中文字幕| 亚洲中文字幕日韩| av网站免费在线观看视频| 亚洲国产欧美网| 亚洲第一青青草原| 丰满人妻熟妇乱又伦精品不卡| 99久久人妻综合| 国产精品成人在线| 99热全是精品| 久久免费观看电影| 我要看黄色一级片免费的| 亚洲第一青青草原| 伦理电影免费视频| 国产主播在线观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 9191精品国产免费久久| 大片电影免费在线观看免费| 午夜91福利影院| 国产亚洲av高清不卡| h视频一区二区三区| av网站免费在线观看视频| svipshipincom国产片| 久久久精品区二区三区| 性高湖久久久久久久久免费观看| 99热网站在线观看| 欧美激情久久久久久爽电影 | 在线观看免费高清a一片| 欧美+亚洲+日韩+国产| 天天添夜夜摸| 色综合欧美亚洲国产小说| av不卡在线播放| 亚洲男人天堂网一区| 国产精品免费大片| 国产在线观看jvid| 亚洲成国产人片在线观看| 午夜两性在线视频| 最新的欧美精品一区二区| 超碰成人久久| 国产伦人伦偷精品视频| 一本久久精品| 国产免费现黄频在线看| 国产av国产精品国产| 精品国产一区二区三区四区第35| 日韩中文字幕视频在线看片| 黑人巨大精品欧美一区二区蜜桃| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩福利视频一区二区| 激情视频va一区二区三区| 69av精品久久久久久 | 在线 av 中文字幕| 久久久精品免费免费高清| 亚洲色图 男人天堂 中文字幕| 欧美另类亚洲清纯唯美| 中文字幕人妻丝袜制服| 精品欧美一区二区三区在线| 国产av国产精品国产| 亚洲欧洲日产国产| 9191精品国产免费久久| 黑人欧美特级aaaaaa片| 亚洲欧美激情在线| 岛国毛片在线播放| 日韩有码中文字幕| 一本综合久久免费| 国产欧美日韩一区二区三 | 大码成人一级视频| 亚洲欧美日韩另类电影网站| 女性生殖器流出的白浆| 婷婷丁香在线五月| 久久人人爽人人片av| 19禁男女啪啪无遮挡网站| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 男女下面插进去视频免费观看| 国产一区二区 视频在线| 亚洲精品国产区一区二| 亚洲av日韩精品久久久久久密| 精品久久蜜臀av无| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区三 | 悠悠久久av| 久久精品成人免费网站| 亚洲专区国产一区二区| av超薄肉色丝袜交足视频| av电影中文网址| 国产又色又爽无遮挡免| 中国国产av一级| 自线自在国产av| 精品熟女少妇八av免费久了| 国产精品久久久久久精品古装| 亚洲精华国产精华精| 咕卡用的链子| 黄频高清免费视频| 99国产精品一区二区三区| 日韩欧美免费精品| 精品卡一卡二卡四卡免费| videosex国产| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频| 欧美另类亚洲清纯唯美| 欧美黑人精品巨大| 国产片内射在线| 国产精品香港三级国产av潘金莲| √禁漫天堂资源中文www| 亚洲国产毛片av蜜桃av| 777米奇影视久久| tube8黄色片| 18禁裸乳无遮挡动漫免费视频| 伊人久久大香线蕉亚洲五| 操出白浆在线播放| 中文字幕制服av| 嫁个100分男人电影在线观看| av电影中文网址| 香蕉丝袜av| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 色视频在线一区二区三区| 日韩欧美一区二区三区在线观看 | 国内毛片毛片毛片毛片毛片| 亚洲专区字幕在线| 精品免费久久久久久久清纯 | 99热国产这里只有精品6| 曰老女人黄片| 人妻 亚洲 视频| bbb黄色大片| 黄色毛片三级朝国网站| 美女中出高潮动态图| 国产精品 国内视频| 91麻豆精品激情在线观看国产 | 亚洲一区二区三区欧美精品| 99久久综合免费| 精品一区二区三卡| 99香蕉大伊视频| 亚洲精品av麻豆狂野| 久久久久久免费高清国产稀缺| 亚洲精品中文字幕一二三四区 | 久久99热这里只频精品6学生| 国产日韩欧美视频二区| 亚洲av美国av| 三级毛片av免费| 免费一级毛片在线播放高清视频 | 亚洲va日本ⅴa欧美va伊人久久 | 性少妇av在线| 男女无遮挡免费网站观看| 18禁黄网站禁片午夜丰满| 欧美亚洲日本最大视频资源| 母亲3免费完整高清在线观看| 久久国产亚洲av麻豆专区| 午夜福利乱码中文字幕| 国产极品粉嫩免费观看在线| 啦啦啦视频在线资源免费观看| 午夜免费观看性视频| 久久久国产成人免费| 一级a爱视频在线免费观看| 国产亚洲欧美在线一区二区| 日本猛色少妇xxxxx猛交久久| 水蜜桃什么品种好| 亚洲精品美女久久久久99蜜臀| 操美女的视频在线观看| 一二三四在线观看免费中文在| 大陆偷拍与自拍| 久久这里只有精品19| 不卡一级毛片| 男女免费视频国产| 9色porny在线观看| av天堂在线播放| 亚洲人成电影免费在线| 五月天丁香电影| 免费高清在线观看日韩| svipshipincom国产片| 久久国产精品男人的天堂亚洲| 少妇精品久久久久久久| 一边摸一边抽搐一进一出视频| 极品少妇高潮喷水抽搐| 麻豆av在线久日| 在线十欧美十亚洲十日本专区| 国产欧美日韩一区二区三 | 99九九在线精品视频| 99re6热这里在线精品视频| 大码成人一级视频| 91字幕亚洲| 亚洲av国产av综合av卡| 在线观看舔阴道视频| 国产精品二区激情视频| 午夜激情久久久久久久| 亚洲成人免费av在线播放| 欧美日韩视频精品一区| 一二三四在线观看免费中文在| 国产伦人伦偷精品视频| 一级毛片精品| 国产精品二区激情视频| √禁漫天堂资源中文www| 女人被躁到高潮嗷嗷叫费观| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 1024香蕉在线观看| 最近最新免费中文字幕在线| 欧美黄色淫秽网站| 亚洲成人免费电影在线观看| 岛国毛片在线播放| 狂野欧美激情性xxxx| 国产精品免费视频内射| 俄罗斯特黄特色一大片| 欧美久久黑人一区二区| 免费少妇av软件| 热99re8久久精品国产| 在线观看www视频免费| 国产一区二区三区综合在线观看| 高清av免费在线| 十八禁高潮呻吟视频| 中国美女看黄片| 青青草视频在线视频观看| www.自偷自拍.com| 国产黄色免费在线视频| 成人国产av品久久久| 久久女婷五月综合色啪小说| 69av精品久久久久久 | 久久久精品94久久精品| 999久久久国产精品视频| 精品免费久久久久久久清纯 | 好男人电影高清在线观看| 满18在线观看网站| 搡老乐熟女国产| 波多野结衣av一区二区av| 黄片大片在线免费观看| 不卡av一区二区三区| videosex国产| 91老司机精品| 久久中文看片网| 色综合欧美亚洲国产小说| 少妇粗大呻吟视频| 亚洲国产欧美一区二区综合| 超碰成人久久| 叶爱在线成人免费视频播放| 色播在线永久视频| 精品国产一区二区三区久久久樱花| 日日摸夜夜添夜夜添小说| 中文字幕人妻丝袜一区二区| 男女高潮啪啪啪动态图| 啦啦啦免费观看视频1| 日韩大片免费观看网站| 自线自在国产av| 亚洲国产欧美一区二区综合| 美国免费a级毛片| www.av在线官网国产| 黄片播放在线免费| 亚洲国产毛片av蜜桃av| tocl精华| 国产无遮挡羞羞视频在线观看| 啪啪无遮挡十八禁网站| 午夜福利影视在线免费观看| 国产一区二区三区在线臀色熟女 | 免费久久久久久久精品成人欧美视频| 老司机在亚洲福利影院| 少妇精品久久久久久久| 天堂8中文在线网| avwww免费| av福利片在线| 我要看黄色一级片免费的| 国内毛片毛片毛片毛片毛片| 真人做人爱边吃奶动态| 在线观看免费视频网站a站| 亚洲成国产人片在线观看| 人人妻,人人澡人人爽秒播| 欧美成狂野欧美在线观看| 亚洲五月色婷婷综合| 中亚洲国语对白在线视频| 午夜激情av网站| 亚洲久久久国产精品| 高潮久久久久久久久久久不卡| 精品国产乱子伦一区二区三区 | cao死你这个sao货| 亚洲欧美一区二区三区久久| 视频区欧美日本亚洲| 热re99久久国产66热| 国产亚洲欧美在线一区二区| 一级毛片女人18水好多| 精品人妻熟女毛片av久久网站| 亚洲第一欧美日韩一区二区三区 | 他把我摸到了高潮在线观看 | 国产欧美日韩综合在线一区二区| 天堂俺去俺来也www色官网| 日韩大码丰满熟妇| a级毛片黄视频| 69av精品久久久久久 | 99国产精品99久久久久| 黄色视频不卡| 亚洲精品av麻豆狂野| 亚洲精品粉嫩美女一区| 亚洲中文日韩欧美视频| 91精品三级在线观看| 黄色a级毛片大全视频| 国产在线视频一区二区| 国产精品一区二区在线不卡| 日韩欧美一区二区三区在线观看 | 国产一级毛片在线| 国产欧美日韩一区二区精品| 12—13女人毛片做爰片一| 国产成人精品无人区| 母亲3免费完整高清在线观看| 欧美精品啪啪一区二区三区 | 在线观看免费日韩欧美大片| 一级黄色大片毛片| 一区二区三区精品91| 激情视频va一区二区三区| 两性夫妻黄色片| 亚洲精品国产av蜜桃| 亚洲人成电影免费在线| 欧美日韩亚洲高清精品| 免费在线观看黄色视频的| 免费在线观看影片大全网站| 99精品久久久久人妻精品| 欧美日韩黄片免| 别揉我奶头~嗯~啊~动态视频 | 国产一区二区三区av在线| 男人舔女人的私密视频| 日本vs欧美在线观看视频| 在线观看一区二区三区激情| 大码成人一级视频| 久久天堂一区二区三区四区| 欧美久久黑人一区二区| www.999成人在线观看| 一级毛片精品| 成人国产一区最新在线观看| 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| 亚洲性夜色夜夜综合| 国产成+人综合+亚洲专区| 天堂中文最新版在线下载| 女性被躁到高潮视频| 欧美激情高清一区二区三区| 欧美 亚洲 国产 日韩一| 欧美xxⅹ黑人| 超碰成人久久| 精品福利永久在线观看| 亚洲情色 制服丝袜| 亚洲久久久国产精品| 丰满迷人的少妇在线观看| 啦啦啦 在线观看视频| 亚洲avbb在线观看| 在线观看人妻少妇| 最新的欧美精品一区二区| 色综合欧美亚洲国产小说| 他把我摸到了高潮在线观看 | 正在播放国产对白刺激| a级毛片黄视频| 18在线观看网站|