亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        中國(guó)近海海溫年際年代際振蕩關(guān)鍵海區(qū)分析研究

        2012-12-23 06:01:06張志華陳幸榮蔡怡
        海洋預(yù)報(bào) 2012年4期
        關(guān)鍵詞:海溫距平年際

        張志華,陳幸榮,蔡怡

        (國(guó)家海洋環(huán)境預(yù)報(bào)中心,北京 100081)

        中國(guó)近海海溫年際年代際振蕩關(guān)鍵海區(qū)分析研究

        張志華,陳幸榮,蔡怡

        (國(guó)家海洋環(huán)境預(yù)報(bào)中心,北京 100081)

        用HADLEY中心的HadISST的海溫資料、NCEP的大氣資料、國(guó)家氣象局發(fā)布的74個(gè)月平均環(huán)流指數(shù)和美國(guó)華盛頓大學(xué)的PDO指數(shù),討論研究了中國(guó)近海海溫的年際年代際振蕩。通過(guò)對(duì)我國(guó)近海海溫的EOF分析,可以發(fā)現(xiàn)我國(guó)近海30?N附近是海溫年際年代際振蕩信號(hào)最強(qiáng)的關(guān)鍵區(qū),它占了總量的58.2%,與澳大利亞?wèn)|北沿海海溫年代際振蕩有相同的周期,最明顯周期是44—45年,滯后PDO信號(hào)近15年左右。該關(guān)鍵區(qū)海溫與澳大利亞?wèn)|北沿海海溫都還存在1年、5年和15年的振蕩,它們的相關(guān)系數(shù)達(dá)到0.604,屬于同一個(gè)模態(tài)。另外,1年的海溫振蕩信號(hào)除了澳大利亞?wèn)|北沿岸以外,還沿著西太平洋海岸線分布,因此中國(guó)近海海溫年際年代際振蕩不是一種局地現(xiàn)象,而與太平洋海溫變化有關(guān)。另外,它是海氣相互作用的結(jié)果,與中高緯度的東亞大槽和低緯度太平洋印度洋的風(fēng)場(chǎng)、西太平洋副高和南海副高、大氣溫度場(chǎng)甚至南半球環(huán)流有明顯的關(guān)系。

        中國(guó)近海海溫;海溫年際年代際振蕩;關(guān)鍵海區(qū)

        1 引言

        由于觀測(cè)資料的缺乏,到目前為止關(guān)于中國(guó)近海海溫的長(zhǎng)期變化研究是比較少的,馮琳等[1]用HadISST1資料分析了1945—2006年?yáng)|中國(guó)海海表溫度(SST)的長(zhǎng)期變化趨勢(shì);張秀芝等[2]用HADLEY中心的海表溫度(SST)資料,分別討論了我國(guó)渤海、黃海、東海、南海北部和南海中部各海域近100年溫度的長(zhǎng)期變化趨勢(shì);他們討論的關(guān)注點(diǎn)主要是集中在全球變暖的背景下,近百年到半個(gè)世紀(jì)以來(lái)的中國(guó)近海海溫的長(zhǎng)期變暖趨勢(shì)。李云芳等[3]用海洋衛(wèi)星遙感資料和EOF分析方法,討論了東海海域海表溫度場(chǎng)的時(shí)空變化,但是由于該資料的水平分布范圍比較小僅僅局限在我國(guó)東海,時(shí)間序列也比較短只有最近20年的資料;白彬人等[4]用中國(guó)近海13個(gè)海洋站1959—2003年逐日日平均海表溫度,討論了中國(guó)近海水溫的短時(shí)(1d)、過(guò)程、旬、月、年際變率以及海溫的長(zhǎng)期趨勢(shì)變化,并在此基礎(chǔ)上詳細(xì)討論了影響中國(guó)近海各時(shí)間尺度水溫變化的天氣氣候因素,但是由于他們用的是單站只有13個(gè)點(diǎn)的資料,因此對(duì)中國(guó)近海的代表性比較差;張松等[5]用美國(guó)NOAA極軌衛(wèi)星中的高級(jí)甚高分辨率輻射計(jì)(AV HRR)反演的海表面溫度資料和EOF方法,分析了渤、黃、東海SST的年際變化,但是它們所用的資料年限是從1985—2007年,這對(duì)于年際變化的研究,所用的資料的年限相對(duì)有點(diǎn)短,因此到目前為止對(duì)于中國(guó)近海洋海溫的研究,人們更多的是注重長(zhǎng)期變暖的趨勢(shì)研究,而對(duì)于年際年代際振蕩的研究,所用的資料年限比較短,而且所用的資料范圍比較小或者只是站點(diǎn)資料,很難準(zhǔn)確的反映中國(guó)近海海溫的年際年代際振蕩。

        然而,隨著全球變暖的愈演愈烈,關(guān)于年際年代際振蕩的研究已越來(lái)越引起人們的關(guān)注,這是因?yàn)槿藗冊(cè)谘芯咳蜃兓瘯r(shí)發(fā)現(xiàn),在全球變暖的同時(shí),年際年代際的氣候振蕩在氣候變化中起了非常重要的作用,如ENSO就是一種熱帶太平洋的海氣相互作用年際振蕩,PDO是北太平洋的一種年代際振蕩,目前比較多研究是關(guān)于ENSO、PDO等洋盆尺度的年際和年代際振蕩[6-10],但是關(guān)于中國(guó)近海局地性的年際年代際振蕩,由于缺少長(zhǎng)時(shí)間序列的觀測(cè)資料,研究還相對(duì)比較少。如果我們對(duì)中國(guó)近海海洋的年際年代際振蕩有清楚的認(rèn)識(shí),那么對(duì)我國(guó)海溫的長(zhǎng)期預(yù)測(cè)將有極大的幫助,本文希望用HADELY資料中心的SST和EOF分析等方法,研究中國(guó)近海海溫的年際年代際振蕩。

        2 資料來(lái)源和處理方法

        我們所用的資料有來(lái)自英國(guó)氣象局HADLEY中心,時(shí)間長(zhǎng)度從1870年1月到2003年12月共134年月平均海面溫度(HadISST),網(wǎng)格分辨率為1°×1°;美國(guó)NCEP大氣500 hPa月平均高度場(chǎng)和850 hPa月平均溫度場(chǎng)和風(fēng)場(chǎng),時(shí)間從1948年1月—2003年12月,網(wǎng)格分辨率為2.5°×2.5°;另外還用到了來(lái)自國(guó)家氣象局發(fā)布的74個(gè)月平均環(huán)流指數(shù),時(shí)間長(zhǎng)度都為1948年1月—2003年12月;美國(guó)華盛頓大學(xué)的PDO指數(shù),時(shí)間為1900年1月—2003年12月。為了略去季節(jié)變化,我們對(duì)所有的原始資料進(jìn)行了氣候平均,用原始資料減去氣候平均,得到氣候平均的距平,下面只用氣候平均的距平進(jìn)行討論。

        3 結(jié)果分析

        圖1 EOF分析的第一個(gè)模態(tài)和時(shí)間系數(shù)

        為了討論中國(guó)近海海溫的年際年代際振蕩,我們?nèi)≈袊?guó)近海的范圍為100°—140°E,0°—45°N,首先用HadISST資料進(jìn)行EOF分析,得到圖1。圖1是EOF分析的第一個(gè)模態(tài)和時(shí)間系數(shù),占了總模態(tài)的58.2%,而第二個(gè)模態(tài)就只占12.1%,可見(jiàn)第一模態(tài)是非常顯著的,因此,本文只討論第一模態(tài)。

        從圖1a中可以看到我國(guó)近海長(zhǎng)江口附近是我國(guó)近海海溫年際年代際振蕩信號(hào)非常強(qiáng)的關(guān)鍵區(qū)域,這可能是因?yàn)闁|海海域位于“亞印太交匯區(qū)”,海陸氣相互作用比較復(fù)雜。該海域既受到來(lái)自北邊的東亞季風(fēng)的影響,又受到來(lái)自南邊的南亞季風(fēng)、澳洲季風(fēng)和太平洋氣流等多種因素影響,是中高緯度和低緯度相互作用的關(guān)鍵區(qū)域。從第一模態(tài)的時(shí)間序列來(lái)看,它不僅反映了近50年的全球變暖,還存在很明顯的周期變化。本文主要討論年際年代際的周期變化,我們將圖1b的時(shí)間系數(shù)進(jìn)行小波分析和功率譜分析得到圖2,從圖中可以看到該模態(tài)的時(shí)間系數(shù)存在最明顯的周期是44—45年(536月)的周期,該周期與太平洋PDO有相同的周期,但是我們分析它們的位相是不同的,在20世紀(jì)PDO在1900—1924年、1947—1976年為負(fù)位相,而1925—1946年、1977—2005年為正位相,從圖2a中可以看到第一時(shí)間序列在1915—1940年、1960—1985年為負(fù)位相,而1895—1915年、1940—1960年、1985—2012為正位相,也就是說(shuō),中國(guó)近海海溫關(guān)鍵區(qū)的年代際振蕩,比PDO滯后了近15年左右,由此我們可以推測(cè),該模態(tài)是一個(gè)獨(dú)立于PDO但是又與PDO有關(guān)的年代際振蕩,它的關(guān)鍵區(qū)域在西太平洋30°N附近的中國(guó)近海。另外,從圖1b第一分量的時(shí)間系數(shù)功率譜分析可以看到,該模態(tài)存在明顯的1年(12個(gè)月)、5年(57個(gè)月)和15年(180個(gè)月)周期。為了討論我國(guó)近海年際年代際振蕩關(guān)鍵區(qū)海溫的變化,我們將122°—123°E,29°—30?N區(qū)域的海溫距平作平均,代表我國(guó)近海年際年代際振蕩關(guān)鍵區(qū)域的海溫距平變化,用該海溫距平時(shí)間序列與整個(gè)太平洋的海溫距平作相關(guān),得到圖3,從圖中可以看到沿著海岸線,我國(guó)近海關(guān)鍵區(qū)域海溫距平與西太平洋沿岸海溫距平都有很好的相關(guān),有幾個(gè)比較大的相關(guān)區(qū)域沿著海岸線從15°S到48°—49°N分布,另外在55°N、142°—145°E處,有一個(gè)明顯的負(fù)相關(guān)區(qū)域,為討論這個(gè)分布模型與我國(guó)近海關(guān)鍵區(qū)的關(guān)系,我們把近海關(guān)鍵區(qū)定義為1區(qū),再分別選取相關(guān)系數(shù)比較大的幾個(gè)區(qū)域?yàn)?—5區(qū)見(jiàn)圖3,其中1區(qū)與2—5區(qū)的相關(guān)系數(shù)分別達(dá)到-0.28(2區(qū)),0.604(3區(qū))0.424(4區(qū)),0.521(5區(qū))。我們分別對(duì)這5個(gè)區(qū)域海溫距平的時(shí)間序列進(jìn)行功率譜分析,發(fā)現(xiàn)1區(qū)和3區(qū)有著同樣的振蕩周期,分別為1年(12個(gè)月)、5年(57個(gè)月)、15年(180個(gè)月)和44—45年(530個(gè)月)周期,而且它們之間的相關(guān)系數(shù)達(dá)到0.604,由此可以推測(cè),中國(guó)近海海溫關(guān)鍵區(qū)的年際年代際振蕩與澳大利亞?wèn)|北沿海的海溫年際年代際振蕩完全屬于同一個(gè)變化模態(tài);中國(guó)近海海溫關(guān)鍵區(qū)的年際年代際振蕩不僅僅是中國(guó)近海的一個(gè)局地現(xiàn)象,另外這5區(qū)域,還都存在明顯的12個(gè)月左右(1年)的周期,因此對(duì)于1年的周期這5個(gè)區(qū)域是由同一個(gè)模態(tài)構(gòu)成,從它們的分布來(lái)看,該變化模態(tài)應(yīng)該與西太平洋沿岸的海岸線地形有關(guān)。

        圖2 EOF分析第一時(shí)間系數(shù)的小波分析結(jié)果和功率譜分析

        圖3 中國(guó)近海海溫關(guān)鍵區(qū)(1區(qū))與太平洋海溫距平相關(guān)系數(shù)分布

        圖4 中國(guó)近海海溫關(guān)鍵區(qū)海溫距平滯后56個(gè)月500 hPa高度場(chǎng)距平相關(guān)系數(shù)分布

        圖5 中國(guó)近海海溫關(guān)鍵區(qū)海溫距平滯后36個(gè)月850 hPa氣溫場(chǎng)距平相關(guān)系數(shù)分布

        圖6 中國(guó)近海海溫關(guān)鍵區(qū)海溫距平滯后44個(gè)月850 hPa緯向風(fēng)場(chǎng)距平相關(guān)系數(shù)分布

        圖7 中國(guó)近海海溫關(guān)鍵區(qū)海溫距平滯后56個(gè)月850 hPa經(jīng)向風(fēng)場(chǎng)距平相關(guān)系數(shù)分布

        從上面的討論可以發(fā)現(xiàn)中國(guó)近海海溫關(guān)鍵區(qū)年際年代際振蕩,并不是一種局地的海溫變化現(xiàn)象,它與整個(gè)太平洋的海溫變化有聯(lián)系,特別是它的年代際振蕩應(yīng)該是一種海氣相互作用現(xiàn)象。為此,我們將1區(qū)的海溫距平與各種大氣環(huán)流指數(shù)作相關(guān),發(fā)現(xiàn)當(dāng)1區(qū)出現(xiàn)最大正的海溫距平65個(gè)月(5年左右)以后,西太平洋副高的西伸位置達(dá)到最西,相關(guān)系數(shù)達(dá)到-0.304,對(duì)應(yīng)的副高面積指數(shù)和強(qiáng)度指數(shù)也達(dá)到比較強(qiáng);另外,我們還研究了南海副高和西太平洋副高對(duì)1區(qū)海溫的影響,發(fā)現(xiàn)當(dāng)西太平洋副高強(qiáng)度最強(qiáng)70個(gè)月以后,1區(qū)的正海溫最大,相關(guān)系數(shù)相關(guān)系數(shù)達(dá)到0.35,而與南海副高的關(guān)系則是當(dāng)南海副高強(qiáng)度最強(qiáng)57個(gè)月以后,1區(qū)的正海溫距平最大,相關(guān)系數(shù)也達(dá)到0.35。由此可見(jiàn),1區(qū)的海溫變化與來(lái)自熱帶副熱帶的環(huán)流變化有很大的關(guān)系。另外,我們還將1區(qū)的海溫與500 hPa的高度距平場(chǎng)作滯后的相關(guān),發(fā)現(xiàn)500 hPa的高度距平場(chǎng)超前1區(qū)的海溫56個(gè)月相關(guān)最好得到圖4,從圖4中可以發(fā)現(xiàn)熱帶印度洋與1區(qū)有非常好的正相關(guān),相關(guān)系數(shù)超過(guò)0.5,這與前面關(guān)于副高與1區(qū)海溫變化的結(jié)論是一致的,而北半球中高緯度東亞大槽則與1區(qū)有負(fù)相關(guān),相關(guān)系數(shù)為-0.2以上,可見(jiàn)1區(qū)的海溫變化與來(lái)自中高緯度的大氣環(huán)流變化也有密切的關(guān)系。

        我們還將1區(qū)的海溫與850 hPa的氣溫距平場(chǎng)作滯后的相關(guān),發(fā)現(xiàn)850 hPa的氣溫距平場(chǎng)超前1區(qū)海溫變化36個(gè)月相關(guān)最好得到圖5。從圖5中可以發(fā)現(xiàn)在中國(guó)近海相關(guān)系數(shù)的分布與圖3相似,可見(jiàn)該信號(hào)是由大氣向海洋傳播,在熱帶印度洋和東太平洋上850 hPa的氣溫距平場(chǎng)超前1區(qū)海溫變化36個(gè)月,相關(guān)系數(shù)都達(dá)到0.5以上;將1區(qū)的海溫與850 hPa的緯向風(fēng)距平場(chǎng)作滯后,發(fā)現(xiàn)850 hPa的緯向風(fēng)距平場(chǎng)超前1區(qū)海溫變化44個(gè)月相關(guān)最好得到圖6。從圖6中可以發(fā)現(xiàn)中國(guó)近海海溫關(guān)鍵區(qū),不但與熱帶東太平洋的緯向風(fēng)有很好的正相關(guān),還與熱帶印度洋的緯向風(fēng)有很好的負(fù)相關(guān),它們的相關(guān)系數(shù)都超0.5,它與東亞大陸中高緯度的緯向風(fēng)也有很好的負(fù)相關(guān)。將1區(qū)的海溫與850 hPa的經(jīng)向風(fēng)距平場(chǎng)作滯后的相關(guān),發(fā)現(xiàn)850 hPa的經(jīng)向風(fēng)距平場(chǎng)超前1區(qū)海溫變化56個(gè)月相關(guān)最好,得到圖7。從圖7可以發(fā)現(xiàn)最大的負(fù)相關(guān)在南半球中緯度澳大利亞附近,相關(guān)系數(shù)超過(guò)-0.5。目前雖然我們還不能確定大氣環(huán)流與中國(guó)近海海溫年際年代際振蕩關(guān)鍵區(qū)之間存在怎樣的聯(lián)系機(jī)制,但是如此高的相關(guān)系數(shù),可以讓我們確定中國(guó)近海海溫關(guān)鍵區(qū)年際年代際振蕩信號(hào),不但與整個(gè)太平洋海溫有關(guān),而且與中高緯度的風(fēng)場(chǎng)、低緯度的副高和風(fēng)場(chǎng)、大氣溫度場(chǎng)甚至南半球環(huán)流有明顯的關(guān)系,這種年際年代際振蕩是海氣相互作用的結(jié)果。

        4 結(jié)論

        通過(guò)對(duì)我國(guó)近海海溫的EOF分析,可以發(fā)現(xiàn)我國(guó)近海30?N附近是海溫年際年代際振蕩信號(hào)最強(qiáng)的關(guān)鍵區(qū),占了總量的58.2%,與澳大利亞?wèn)|北沿海海溫年代際振蕩有相同的周期,最明顯周期是44—45年,滯后PDO信號(hào)近15年左右,該關(guān)鍵區(qū)海溫與澳大利亞?wèn)|北沿海海溫都還存在1年、5年和15年的振蕩,它們的相關(guān)系數(shù)達(dá)到0.604,屬于同一個(gè)模態(tài);另外,1年的海溫振蕩信號(hào)除了澳大利亞?wèn)|北沿岸以外,還沿著西太平洋海岸線分布,因此中國(guó)近海海溫年際年代際振蕩不是一種局地現(xiàn)象,而與太平洋海溫變化有關(guān)。另外,它是海氣相互作用的結(jié)果,與中高緯度的東亞大槽和低緯度太平洋印度洋的風(fēng)場(chǎng)、西太平洋副高和南海副高、大氣溫度場(chǎng)甚至南半球環(huán)流有明顯的關(guān)系。但是關(guān)于中高緯度和低緯度的風(fēng)場(chǎng)、低緯度的副高、大氣溫度場(chǎng)甚至南半球環(huán)流對(duì)我國(guó)近海海溫關(guān)鍵區(qū)年際年代際振蕩的影響機(jī)制還需要作進(jìn)一步的研究。

        [1]馮琳,林霄沛.1945—2006年?yáng)|中國(guó)海海表溫度的長(zhǎng)期變化趨勢(shì)[J].中國(guó)海洋大學(xué)學(xué)報(bào),2009,39(1):13-18.

        [2]張秀芝,裘越芳,吳迅英.近百年中國(guó)近海海溫變化[J].氣候與環(huán)境研究,2005,4:799-807.

        [3]李云芳,郭佩芳,李廣云.東海海域海表溫度場(chǎng)的EOF分析[J].海洋湖沼通報(bào),2007,3:1-4.

        [4]白彬人,宋家喜.中國(guó)近海沿岸海溫多時(shí)間尺度變率及影響其變化的天氣氣候因素[J].海洋預(yù)報(bào),2005,22(4):78-88.

        [5]張松,于非,刁新源,等.渤、黃、東海海表面溫度年際變化特征分析[J].海洋科學(xué),2009,33(8):76-81.

        [6]鐘姍姍,何金海,劉宣飛.太平洋次表層海溫年代際變率及其突變特征[J].南京氣象學(xué)院學(xué)報(bào),2002,25(5):595-602.

        [7]陳幸榮,王彰貴,巢紀(jì)平,等.CCSM對(duì)太平洋年代際振蕩的敏感性試驗(yàn)[J].海洋學(xué)報(bào),2011,33(5):23-31.

        [8]陳幸榮,王彰貴,巢紀(jì)平,等.PDO的三維空間結(jié)構(gòu)和實(shí)際演變特征[J].海洋學(xué)報(bào),2011,33(6):19-28.

        [9]蔡怡,李海,張人禾.熱帶印度洋海溫異常與ENSO關(guān)系的進(jìn)一步研究[J].氣象學(xué)報(bào),2008,66(1):120-124.

        [10]蔡怡,張建華,于衛(wèi)東.赤道太平洋-印度洋海洋上層海溫分析[J].熱帶海洋學(xué)報(bào),2005,24(4):60-66。

        A study of SST decadal oscillation in China Sea

        ZHANG Zhi-hua,CHEN Xing-rong,CAI Yi
        (National Marine Environmental Forecasting Center,Beijing 100081 China)

        Based on the HadISST data from HADLEY Center,geopotential height at 500 hPa and air temperature and wind at 850hPa from NCEP,74-month averaged circulation index from Chinese Meteorology Administration and PDO index from Washington University,the SST decadal oscillation in China Sea is analyzed by EOF method.The results show that the key area is close to 30°N in the East of the China Sea,where showing strongest SST decadal oscillation and.It takes account about 58.2%of total amount.The period of SST decadal oscillation in the key area of the China Sea is same as that in the northeast coast of the Australia Ocean.The period of the strongest signal SST decadal oscillation is 44—45 years,showing a lag of 15 years with PDO.The SST oscillation in the key area of the China Sea and in the northeast coast of the Australia Ocean share the same periods with 1 year,5 years and 15 years and the correlation coefficient is 0.604,suggesting that they are in the same mode. There is also 1 year period in the western Pacific coast.The results indicate that the SST decadal oscillation in the key area of the China Sea is not a local event,but connects to the sea temperature in the Pacific Ocean.It is the results of the air—sea interaction,being related with the wind from mid-high and low latitude,the subtropical high in the west Pacific Ocean and the South China Sea,the air temperature and the circulation in the southern hemisphere.

        SST;decadal oscillation;key sea area

        P732

        A

        1003-0239(2012)04-0001-06

        2012-03-01

        國(guó)家自然科學(xué)基金(40976015);國(guó)家“973”項(xiàng)目(2010CB950302)

        張志華(1976-),男,助理研究員,主要從事工作:ENSO動(dòng)力學(xué)和海洋數(shù)值模擬研究。E-mail∶zhangzhihua@nmefc.gov.cn

        猜你喜歡
        海溫距平年際
        颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對(duì)比分析
        基于深度學(xué)習(xí)的海溫觀測(cè)數(shù)據(jù)質(zhì)量控制應(yīng)用研究
        基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
        北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
        南印度洋偶極型海溫與中國(guó)西南地區(qū)初秋降水的關(guān)系
        近40年阿里地區(qū)云量和氣溫的年際變化
        西藏科技(2018年9期)2018-10-17 05:51:30
        甘肅省降水和冰雹天氣氣候分析
        亞洲夏季風(fēng)的年際和年代際變化及其未來(lái)預(yù)測(cè)
        與北大西洋接壤的北極海冰和年際氣候變化
        基于距平百分率的萬(wàn)源市氣象干旱灰色預(yù)測(cè)
        亚洲欧美日韩国产精品一区| 国产超碰女人任你爽| 天堂草原电视剧在线观看图片高清| 免费的一级毛片| 亚洲蜜桃视频在线观看| 中文字幕午夜精品久久久| 久久www免费人成精品| 大地资源网更新免费播放视频| 亚洲美女国产精品久久久久久久久| 日韩精品中文字幕第二页 | 亚洲精品乱码久久久久久中文字幕| 国产手机在线αⅴ片无码观看| 国产精品久久无码免费看 | 国产三级av大全在线爽| 日韩中文字幕乱码在线| 久久红精品一区二区三区| 少妇性饥渴无码a区免费| 91免费永久国产在线观看| 久久精品这里就是精品| 亚洲av日韩一区二区| 久久久av精品波多野结衣| 亚洲AV无码精品色欲av| 求网址日韩精品在线你懂的| 日日碰日日摸日日澡视频播放| 美女又色又爽视频免费| 久久国产影视免费精品| 久久久精品少妇—二区| 中国美女a级毛片| 国产顶级熟妇高潮xxxxx| 亚洲精品有码在线观看| 国产毛片精品一区二区色| 欧美性白人极品1819hd| 国产成人无码区免费内射一片色欲| 一区二区三区内射视频在线观看 | 免费av网站大全亚洲一区| av人摸人人人澡人人超碰妓女| 激情久久无码天堂| 国产在线一区二区三区不卡| 国产不卡视频一区二区三区| 日本高清www午色夜高清视频| 日韩在线精品在线观看|