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

    全新世氣候?qū)η瓿叨忍柣顒拥目赡茼憫?yīng)

    2014-03-02 05:25:05尹志強韓延本王博
    關(guān)鍵詞:代用太陽活動小波

    尹志強 韓延本 王博

    (中國科學(xué)院國家天文臺,北京 100012)

    全新世氣候?qū)η瓿叨忍柣顒拥目赡茼憫?yīng)

    尹志強 韓延本 王博

    (中國科學(xué)院國家天文臺,北京 100012)

    利用熟知的山洞石筍資料及重建的太陽黑子數(shù)獲取的時間序列,討論了長時間尺度上氣候?qū)μ柣顒幼兓目赡茼憫?yīng)。小波分析顯示氣候變化與太陽活動變化在約1000及約2000年周期上具有一致性,這可能反映了全新世時期太陽活動的調(diào)制作用。另外,小波功率譜顯示在公元前2000年附近缺失千年周期,這種突變及氣候條件的長期變化可能給許多地區(qū)帶來社會性影響。

    古氣候,太陽黑子數(shù),全新世

    1 引言

    眾所周知,氣候系統(tǒng)受許多因素的影響,例如地球軌道變化、太陽活動變化、大型熱帶火山噴發(fā)、地面覆蓋差異及溫室氣體等。其中,太陽活動的變化是關(guān)鍵因素之一。目前有許多有關(guān)太陽活動變化與各種氣候記錄相關(guān)性的研究成果,但涉及千年時間尺度的太陽活動驅(qū)動氣候變化的機制與該驅(qū)動下氣候系統(tǒng)響應(yīng)之間關(guān)系的研究相對較少[1-4]。

    由于缺少更長時間的直接觀測或儀器記錄數(shù)據(jù),我們只能依靠代用指標(biāo)來研究全新世時期的古氣候。例如,樹木年輪、珊瑚、鐘乳石及冰芯等代用指標(biāo)通常能很好地反映太陽活動、氣溫、濕度及季風(fēng)強度等[5-7]。也有許多學(xué)者曾嘗試從陸地及海洋代用指標(biāo)中獲取十年至百萬年時間尺度的氣候變化,研究太陽活動對氣候變化的驅(qū)動[3,7]。

    對于千年周期尺度的氣候變化,格陵蘭冰芯及北大西洋的巖芯很好地記錄了全新世期間氣候變化[8-9]。許多學(xué)者也在古氣候資料,例如在Big Lake沉積物記錄[10]、加拿大[11]、北美花粉記錄[12]中找到了這些變化周期。Clemens等[13]認(rèn)為千年尺度的氣候變化是由太陽活動世紀(jì)周期外差得來的。Bütikofer等[14]研究表明,全新世氣候記錄的典型特征是寬波段的準(zhǔn)周期變化,而不具有某一顯著的譜峰。

    本文對收集到的一些數(shù)據(jù)做了小波分析,以進一步研究全新世千年尺度的氣候變化,并討論存在的周期振蕩及氣候?qū)μ柣顒涌赡茼憫?yīng)的時變特征。

    2 數(shù)據(jù)及分析方法

    2.1 重建的太陽黑子數(shù)

    我們可以使用宇宙成因核素(如14C,10Be等)含量變化的高頻分量作為太陽活動變化的代用指標(biāo)。宇宙成因核素是由宇宙線輻射與大氣分子相互作用形成的,其含量受太陽活動強度變化的調(diào)制。Solanki等[15]根據(jù)中緯度地區(qū)樹木年輪10年采樣間隔的高精度14C資料分析獲取的δ14C記錄,給出了覆蓋全新世的太陽黑字?jǐn)?shù)(SSN)的重建結(jié)果。該數(shù)據(jù)在公元1600年之后的結(jié)果與實際觀測有很好的一致性。雖然10年采樣間

    隔的數(shù)據(jù)不能用來分析太陽活動的高頻波動,但足以分析百年及更長的周期性振蕩。

    2.2δ18O記錄

    從冰芯、有孔蟲殼或石筍中獲取的穩(wěn)定的氧同位素比值(18O/16O)記錄被廣泛用于古氣候的研究。千年尺度變化的分析對數(shù)據(jù)的要求更為嚴(yán)格,我們需要更長時間跨度及更高分辨率的信息。根據(jù)Fleitmann等[16]的分析報告,選取了兩個高分辨率及較長時間跨度的低緯度δ18O記錄。其中,董哥(Dongge)洞(25°17′N,108°5′E)的數(shù)據(jù)被認(rèn)為記錄了過去9千年亞洲季風(fēng)的信息,數(shù)據(jù)采樣間隔從1年至10年不等,平均采樣間隔為4年。該記錄與大氣14C記錄的十年至百年尺度的互相關(guān)分析顯示,該頻率的一些季風(fēng)變化被認(rèn)為是由太陽活動的變化引起的[1,17]。為了便于進一步分析,我們利用分段3次Hermite方法對該數(shù)據(jù)做了內(nèi)插處理,得到5年間隔的數(shù)據(jù)序列。另一個從阿曼北部與南部及也門的4個山洞的全新世石筍中獲取的高分辨率δ18O數(shù)據(jù),其提供了沿12°—23°N方向的沉積作用的詳細波動信息[16,18-19]。我們同樣利用分段3次Hermite方法以10年分辨率對該數(shù)據(jù)做了內(nèi)插處理。下文中這兩個代用指標(biāo)分別簡稱為Dongge及Qunf。圖1繪出了所有數(shù)據(jù)的曲線,以展示強或弱的突變事件。

    2.3 小波分析

    為了估計一些數(shù)據(jù)頻率的非平穩(wěn)性,Morlet提出了小波變換。小波是一類用于對一給定的函數(shù)在空間及時間域上做局部化分析的函數(shù)。Morlet小波將三角函數(shù)與高斯函數(shù)相結(jié)合,通過時頻域上的小波濾波,將時間序列的信號與噪聲分離開。這樣可以有效消除噪聲的影響[20-22]。我們在研究中利用小波變換的方法對上述代用指標(biāo)的時變信息在時間及頻率域上分別進行了估計。

    圖1 SSN(紅)、Dongge(藍)及Qunf(黑)時間序列的變化曲線

    3 結(jié)果及討論

    圖2給出了SSN,Dongge及Qunf數(shù)據(jù)序列的小波變換結(jié)果。由于本文主要對百年及千年尺度長期波動的變化特征做研究,故把小波分析的周期范圍選為64~4096年。小波圖(圖2a,2c,2e)中,背景置為白色,由藍色至紅色表示能量由低到高的變化。利用時間序列的協(xié)方差對小波功率做了歸一化處理,以給出功率相對白噪聲的大小。圖2b,2d,2f分別給出了3個序列的全局小波譜。

    圖2a顯示全新世期間SSN時間序列具有明顯的長周期波動(包括約220,517,1000及2000年等)。雙世紀(jì)周期的變化在190~220年之間,千年周期長度的變化在900~1100年之間,2000年周期長度的變化幅度為2040~2400年。圖2b的全局小波變換顯示所有的周期都在置信水平以上。圖2c—d顯示Dongge時間數(shù)據(jù)具有類似的特征。Dongge序列的主要周期約為210,567,1000,2320年。Qunf序列的分析顯示其幾個顯著的周期約為:1000,1820及2400年(如圖2e—f)。幾乎所有的δ18O的主要變化周期與SSN的變化周期接近。

    全新世的十年與百年尺度上的10Be及14C中的約90年的Gleisberg周期及約200年的Suess周期以及其他太陽活動的記錄,同樣很好地在環(huán)境代用指標(biāo)記錄中展現(xiàn)出來[7,23-28]。在百年及千年尺度上的小波分析顯示SSN的約200年周期幅度較大的時段主要為公元前約8000—6000年間、公元前約4000—2700年間、公元前約700—50年間,以及公元1500年前后的幾個時段內(nèi)。

    Dongge記錄同樣顯示出約200年的周期(圖2c)。雖然該記錄的譜分析顯示出顯著的208年及86年的太陽周期,且認(rèn)為其受太陽驅(qū)動活動的影響[17],但局部小波功率譜分析顯示該周期在Dongge與SSN數(shù)據(jù)中時間不一致(圖2a與2c)。圖2e顯示了Qunf序列的時變特征,我們發(fā)現(xiàn)在置信水平之上并不存在約200年周期(圖2f)。

    上述的分析表明,約1000年及2000年的周期波動是SSN,Dongge及Qunf序列的更長期的變化特征,其都在置信水平之上。與其他波動相比,這些振蕩似乎更加穩(wěn)定。尤其是3個序列的約1000年周期幾乎在整個時間跨度上一直存在。另外,Dongge的周期長度在公元前600年附近較短,在公元前2000年之后稍微變長。對Qunf記錄,在公元前2000年之前,該周期長度更為穩(wěn)定,隨后稍微變短。這3個時間序列的局部小波功率譜顯示,在公元前2000年附近,約1000年周期有明顯的下降(圖2)。Debret等[29]給出了14C含量[10]

    及10Be[30]的類似分析結(jié)果,主要的約1000年周期同樣在5000-0 B.P.呈現(xiàn)衰減。

    為了展示這些序列的千年尺度波動,我們對每個時間序列做了歸一化尺度平均小波功率譜,圖3繪出了相應(yīng)的結(jié)果(最大值為1,最小值為0)。千年尺度變化的小波尺度約為900~1100年。

    由圖3可以看到,Dongge及Qunf在千年尺度上有明顯的相關(guān)性。公元前7000—1000年間,Dongge序列的波動幅度與Qunf非常一致。公元前8500—7000年間,SSN的千年準(zhǔn)周期幅度相對較大,隨后變小。雖然Dongge及Qunf的小波譜的平均能量似乎與SSN的相關(guān)性較差,它們在公元前2000年附近的最小值可能與上述的SSN千年尺度缺失相關(guān)。

    連續(xù)小波變換結(jié)果表明,在過去的11000年期間,連續(xù)的千年尺度振蕩周期在1000~2000年左右。它們與在北大西洋發(fā)現(xiàn)的全新世約1500年周期不同[9-10]。另外,除Dongge與Qunf記錄的夏季日照的趨勢性變化外,兩數(shù)據(jù)記錄的許多突變事件也與14C序列相關(guān)[1,16,31],我們的分析認(rèn)為顯著的約1000年波動以及由這兩個記錄獲取的全新世約2000年波動可能來自于太陽活動的驅(qū)動。Debret等[29]同樣給出了北大西洋記錄的類似結(jié)果。

    圖2 SSN、Dongge及Qunf數(shù)據(jù)序列的小波變換。(a)、(c)、(e)分別為SSN、Dongge、Qunf時間序列的局部小波功率譜;黑色虛線為Morlet小波變換端部效應(yīng)的邊界,該邊界以下的小波譜的可靠性較差;黑色實線為95%置信水平的置信線。(b)、(d)、(f)分別為SSN、Dongge、Qunf時間序列的小波全局功率譜,虛線為以白噪聲為背景的95%置信水平的置信線

    圖3 SSN、Dongge及Qunf時間序列的歸一化尺度平均小波功率譜

    4 結(jié)論

    通過再分析的中低緯度山洞石筍的δ18O數(shù)據(jù)序列,小波分析顯示全新世氣候的千年尺度的變化至少存在約一千年和約兩千年的兩個周期波動。我們的結(jié)果表明這兩個周期可能由太陽活動驅(qū)動。雖然譜分析發(fā)現(xiàn)δ18O記錄呈現(xiàn)出顯著的太陽活動百年峰值[17],然而小波分析顯示它們在同一時間內(nèi)與太陽活動周期并沒有明顯的相關(guān)性。從這一結(jié)果來看,全新世期間在千年尺度上,反應(yīng)亞洲季風(fēng)的δ18O代用指標(biāo),ITCZ或印度夏季季風(fēng)的平均緯度位置是由太陽活動驅(qū)動的。與百年尺度的變化相比,千年尺度可能是氣候系統(tǒng)變

    化的更基本的周期,這種波動很少受如復(fù)雜地形,地面覆蓋,火山驅(qū)動及溫室氣體等其他因素的影響[4]??紤]到在中緯度[32]及高緯度[29]各種氣候代用指標(biāo)發(fā)現(xiàn)的千年尺度太陽驅(qū)動[33],我們認(rèn)為這種影響可能不是一種局部特征,而是廣泛地存在于整個氣候系統(tǒng)。

    致謝:感謝Solanki教授提供重建的太陽黑子數(shù)序列,同時感謝NOAA古氣候項目提供替代數(shù)據(jù)。小波變換程序由C. Torrence及G. Comp開發(fā),可通過如下鏈接獲?。篽ttp://paos.colorado.edu/research/wavelets/。

    [1]Wang Y, Cheng H, Edwards R L, et al. The Holocene Asian monsoon: links to solar changes and North Atlantic climate. Science, 2005, 308(5723): 854-857.

    [2] Rind D. The Sun’s role in climate variations. Science, 2002, 296(5568): 673-667.

    [3]Bard E, Frank M. Climate change and solar variability: what’s new under the Sun. Earth and Planetary Science Letters, 2006, 248(1-2): 1-14.

    [4]Wanner H, Beer J, Bütikofer J, et al. Mid- to Late Holocene climate change: an overview. Quat Sci Rev, 2008, 27(19-20): 1791-1828.

    [5]Beer J, Mende W, Stellmacher R. The role of the sun in climate forcing. Quat Sci Rev, 2000, 19(1-5): 403-415.

    [6]Mayewski P A, Rohlingb E E, Stager J C, et al. Holocene climate variability. Quaternary Research, 2004, 62(3): 243-255.

    [7]Versteegh G J M. Solar forcing of climate 2: evidence from the past. Space Science Reviews, 2005, 120(3-4): 243-286.

    [8] O’Brien S R, Mayewski P A, Meeker L D, et al. Complexity of Holocene climate as reconstructed from a Greenland ice core. Science, 1995, 270(5244): 1962-1964.

    [9]Bond G, Showers W, Cheseby M, et al. A Pervasive millennialscale cycle in North Atlantic Holocene and glacial climates. Science, 1997, 278(5341): 1257-1266.

    [10]Bond G, Kromer B, Beer J, et al. Persistent solar in fl uence on North Atlantic climate during the Holocene. Science, 2001, 294(5549): 2130-2136.

    [11]Cumming B F, Laird K R, Bennett J R, et al. Persistent millennialscale shifts in moisture regimes in western Canada during the past six millennia. Proceedings of the National Academy of Sciences, 2002, 99(25): 16117-16121.

    [12]Viau A E, Gajewski K, Fines P, et al. Widespread evidence of 1500 yr climate variability in North America during the past 14000 yr. Geology, 2002, 30: 455-458.

    [13]Clemens S C. Millennial-band climate spectrum resolved and linked to centennial-scale solar cycles. Quat Sci Rev, 2005, 24(5-6): 521-531.

    [14] Bütikofer J. [Online] Available:http://www.giub.unibe.ch/klimet/ docs/diplom_jbuetikofer.pdf (2014/1/1).

    [15]Solanki S K, Usoskin I G, Kromer B, et al. Unusual activity of the Sun during recent decades compared to the previous 11,000 years. Nature, 2004, 431: 1084-1087.

    [16]Fleitmann D, Burns S. J, Mangini A, et al. Holocene ITCZ and Indian monsoon dynamics recorded in stalagmites from Oman and Yemen (Socotra). Quat Sci Rev, 2007, 26(1-2): 170-188.

    [17]Dykoski C A, Edwards R L, Cheng H, et al. A high-resolution, absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave, China. Earth and Planetary Science Letters, 2005, 233(1-2): 71-86.

    [18]Neff U, Burns S J, Mangini A, et al. Strong coherence between solar variability and the monsoon in Oman between 9 and 6 kyr ago. Nature, 2001, 411: 290-293.

    [19]Fleitmann D, Burns S J, Mudelsee M, et al. Holocene Forcing of the Indian Monsoon Recorded in a Stalagmite from Southern Oman. Science, 2003, 300(5626): 1737-1739.

    [20]Kumar P, Foufoula-Georgiou E. Wavelet analysis for geophysical applications. Rev Geophys, 1997, 35(4): 385-412.

    [21]Torrence C, Compo G. A Practical guide to wavelet analysis. Bull Amer Meteor Soc, 1998, 79(1): 61-78.

    [22]Yin Z Q, Ma L H, Han Y B, et al. Long-term variations of solar activity. Chin Sci Bull, 2007, 52(20): 2737-2741.

    [23]Mergentaler J. On the long-period changes of solar activity. Acta Astronomica, 1984, 34(2): 263-267.

    [24]Ogurtsov M G, Nagovitsyn Y A, Kocharov G E, et al. Long-period cycles of the Sun’s activity recorded in direct solar data and proxies. Sol Phys, 2002, 211(1-2), 371-394.

    [25]Ogurtsov M G, Jungner H, Kocharov G E, et al. On the link between northern fennoscandian climate and length of the quasieleven-year cycle in galactic cosmic-ray flux. Sol Phys, 2003, 218(1-2): 345-357.

    [26]Raspopov O M, DergachevV A, Kolstr?m T. Periodicity of climate conditions and solar variability derived from dendrochronological and other palaeoclimatic data in high latitudes. Palaeogeogr Palaeoclimatol Palaeoecol, 2004, 209(1-4): 127-139.

    [27]Velasco V M, Mendoza B. Assessing the relationship between solar activity and some large scale climatic phenomena. Advances in Space Research, 2008, 42(5): 866-878.

    [28]Ma L H, Vaquero J M. Is the Suess cycle present in historical naked-eye observations of sunspots? New Astronomy, 2009, 14(3): 307-310.

    [29]Debret M, Bout-Roumazeilles V, Grousset F, et al. The origin of the 1500-year climate cycles in Holocene North-Atlantic records. Clim Past, 2007, 3(4): 569-575.

    [30]Vonmoos M, Beer J, Muscheler R J. Large variations in Holocene solar activity: constraints from 10Be in the Greenland Ice Core Project ice core. Journal of Geophysical Research: Space Physics, 2006, 111(A10): A10105.

    [31]Wang Y, Cheng H, Edwards R L, et al. Millennial- and orbitalscale changes in the East Asian monsoon over the past 224,000 years. Nature, 2008, 451: 1090-1093.

    [32]Nederbragt A J, Thurow J. Geographic coherence of millennialscale climate cycles during the Holocene. Palaeogeography Palaeoclimatology Palaeoecology, 2005, 221(3-4): 313-324.

    [33]Usoskin I G, Solanki S K, Kovaltsov G A. Grand minima and maxima of solar activity: new observational constraints. A&A, 2007, 471(1): 301-309.

    Possible Response of Climate to Solar Activity on Millennial Scale during the Holocene

    Yin Zhiqiang, Han Yanben, Wang bo
    (National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012)

    Based on well-known series from land-based caves and reconstructed sunspot numbers, the nature of climatic response to solar variability is discussed over a long-time scale. Wavelet analysis reveals coherence between the climatic and solar oscillations with periods of ~1000- and ~2000-yr that may ref l ect a modulation of solar activity during the Holocene. Moreover, the wavelet power spectra present that ~1000-yr cycle of these series is absent around ~2000 BC and this abrupt change and prolonged changes in climatic conditions might have brought about major societal challenges to many regions.

    paleoclimate, Sun sunspots, Holocene

    10.3969/j.issn.2095-1973.2014.04.002

    2014年1月1日;

    2014年3月14日

    尹志強(1979—),E-mail: yinzhq@bao.ac.cn

    資助信息:國家重大科學(xué)研究計劃項目(2012CB957801)

    猜你喜歡
    代用太陽活動小波
    構(gòu)造Daubechies小波的一些注記
    第24太陽活動周中國地區(qū)電離層閃爍統(tǒng)計特性研究
    第23和24太陽活動周高緯地磁感應(yīng)電流分布特性
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    平糖仙草代用茶 降糖降脂降血壓
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    點燃式發(fā)動機氣體代用燃料的數(shù)值模擬分析
    河南科技(2015年2期)2015-02-27 14:20:28
    淺析壓力容器制造過程中的材料代用
    基于FPGA小波變換核的設(shè)計
    電測與儀表(2014年8期)2014-04-04 09:19:38
    榜單
    男插女下体视频免费在线播放| 美女 人体艺术 gogo| 欧美在线一区亚洲| 成人av在线播放网站| 久久午夜亚洲精品久久| 精品一区二区三区四区五区乱码| 波多野结衣巨乳人妻| 国产午夜福利久久久久久| 久久这里只有精品中国| www.www免费av| 亚洲成人久久爱视频| 精品无人区乱码1区二区| 久久国产精品人妻蜜桃| 亚洲精品在线美女| 亚洲国产欧洲综合997久久,| 十八禁人妻一区二区| 国产精品久久电影中文字幕| 精品久久久久久久久久免费视频| 久久久久久免费高清国产稀缺| 男人舔女人下体高潮全视频| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 日日干狠狠操夜夜爽| 成年免费大片在线观看| 欧美极品一区二区三区四区| 色综合婷婷激情| 国产高清视频在线播放一区| 中文字幕人成人乱码亚洲影| 亚洲欧美一区二区三区黑人| 成人亚洲精品av一区二区| 欧美日韩国产亚洲二区| 日韩免费av在线播放| 久久精品影院6| 99国产精品一区二区蜜桃av| 狂野欧美激情性xxxx| 人人妻人人看人人澡| 一进一出抽搐gif免费好疼| 亚洲18禁久久av| 亚洲精华国产精华精| 免费看日本二区| 女人爽到高潮嗷嗷叫在线视频| 波多野结衣高清无吗| 亚洲精品粉嫩美女一区| 精品一区二区三区四区五区乱码| 一区二区三区国产精品乱码| 日本五十路高清| 90打野战视频偷拍视频| 欧美黄色淫秽网站| 婷婷亚洲欧美| 亚洲自拍偷在线| 久久中文字幕人妻熟女| 精品欧美一区二区三区在线| 午夜激情av网站| 亚洲精品国产一区二区精华液| 久久午夜综合久久蜜桃| 全区人妻精品视频| 天堂动漫精品| 99久久国产精品久久久| 一本综合久久免费| 一a级毛片在线观看| 国产成人精品久久二区二区91| 国产精品亚洲av一区麻豆| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| cao死你这个sao货| 欧美大码av| 一边摸一边做爽爽视频免费| 日韩三级视频一区二区三区| 亚洲欧美激情综合另类| 婷婷亚洲欧美| 国产成年人精品一区二区| 桃红色精品国产亚洲av| 久久 成人 亚洲| 日本免费一区二区三区高清不卡| 久久人妻福利社区极品人妻图片| 亚洲av电影不卡..在线观看| 在线观看免费视频日本深夜| 变态另类丝袜制服| 国产区一区二久久| 亚洲色图av天堂| 久久亚洲精品不卡| 成年人黄色毛片网站| 夜夜躁狠狠躁天天躁| 身体一侧抽搐| 亚洲精品粉嫩美女一区| 欧美色欧美亚洲另类二区| 在线观看舔阴道视频| 国产男靠女视频免费网站| 舔av片在线| 最近视频中文字幕2019在线8| 1024手机看黄色片| 亚洲成av人片在线播放无| 久久久久亚洲av毛片大全| 久久久精品欧美日韩精品| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲美女久久久| 日韩欧美精品v在线| 搡老妇女老女人老熟妇| 最新美女视频免费是黄的| 国产一级毛片七仙女欲春2| 国产欧美日韩精品亚洲av| 一级作爱视频免费观看| 特级一级黄色大片| 精品国产美女av久久久久小说| 俺也久久电影网| 精品一区二区三区av网在线观看| 婷婷亚洲欧美| 好男人电影高清在线观看| 少妇粗大呻吟视频| 精品电影一区二区在线| 国内揄拍国产精品人妻在线| 亚洲精品美女久久久久99蜜臀| 午夜亚洲福利在线播放| 一级作爱视频免费观看| 国产精品永久免费网站| 日韩国内少妇激情av| 老熟妇仑乱视频hdxx| 亚洲片人在线观看| 午夜福利免费观看在线| 1024视频免费在线观看| 国产单亲对白刺激| 欧美最黄视频在线播放免费| 欧美日韩精品网址| 欧美又色又爽又黄视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成网站在线播放欧美日韩| 欧美日韩中文字幕国产精品一区二区三区| 精品日产1卡2卡| 在线观看66精品国产| 中文字幕人妻丝袜一区二区| 可以免费在线观看a视频的电影网站| 欧美日韩瑟瑟在线播放| 久久精品国产亚洲av香蕉五月| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 亚洲国产看品久久| av福利片在线| 亚洲欧美日韩东京热| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 欧美乱码精品一区二区三区| 久久精品国产清高在天天线| 久久久久性生活片| 在线a可以看的网站| 哪里可以看免费的av片| 国产97色在线日韩免费| 一级毛片精品| 又黄又爽又免费观看的视频| 亚洲欧美精品综合久久99| 中文字幕高清在线视频| 波多野结衣高清无吗| 老汉色av国产亚洲站长工具| 天天添夜夜摸| 制服诱惑二区| 国产高清videossex| 国产精品亚洲美女久久久| 久久国产乱子伦精品免费另类| 欧美一区二区国产精品久久精品 | 国产成人av激情在线播放| 黄色成人免费大全| 久久国产乱子伦精品免费另类| 国产蜜桃级精品一区二区三区| 久久亚洲真实| 亚洲av成人精品一区久久| 999久久久精品免费观看国产| 国产av一区二区精品久久| 免费观看人在逋| 一级片免费观看大全| 国内精品久久久久久久电影| 伦理电影免费视频| 亚洲欧美日韩东京热| 村上凉子中文字幕在线| 国产精品自产拍在线观看55亚洲| av福利片在线| 99精品在免费线老司机午夜| 久久精品91蜜桃| 精品国产美女av久久久久小说| 一级片免费观看大全| 蜜桃久久精品国产亚洲av| 看黄色毛片网站| 99久久久亚洲精品蜜臀av| 天堂动漫精品| 久99久视频精品免费| 夜夜躁狠狠躁天天躁| 久久欧美精品欧美久久欧美| 日韩高清综合在线| 日韩国内少妇激情av| 99热只有精品国产| 国产视频内射| 欧美三级亚洲精品| 一夜夜www| 精品第一国产精品| 五月伊人婷婷丁香| 久久九九热精品免费| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看| 亚洲人与动物交配视频| 色老头精品视频在线观看| 欧美黑人欧美精品刺激| www国产在线视频色| 久热爱精品视频在线9| 国产真人三级小视频在线观看| 亚洲一码二码三码区别大吗| 亚洲熟妇熟女久久| 欧美日韩亚洲国产一区二区在线观看| 不卡一级毛片| a级毛片在线看网站| 久久 成人 亚洲| 在线观看一区二区三区| e午夜精品久久久久久久| 亚洲色图av天堂| 国产一区二区激情短视频| 午夜福利成人在线免费观看| 亚洲av成人不卡在线观看播放网| 桃红色精品国产亚洲av| 亚洲色图av天堂| 国产三级黄色录像| 午夜久久久久精精品| 波多野结衣高清无吗| 欧美成人一区二区免费高清观看 | 精品少妇一区二区三区视频日本电影| 成在线人永久免费视频| 99精品在免费线老司机午夜| 久久中文看片网| av片东京热男人的天堂| 老汉色∧v一级毛片| 51午夜福利影视在线观看| 黄色a级毛片大全视频| 精品久久久久久,| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 国产一区二区在线av高清观看| 色综合欧美亚洲国产小说| 一进一出抽搐gif免费好疼| 一个人观看的视频www高清免费观看 | 国产精品综合久久久久久久免费| 亚洲国产欧美人成| 蜜桃久久精品国产亚洲av| 国产三级中文精品| 香蕉国产在线看| 亚洲一区中文字幕在线| 国产精品精品国产色婷婷| 久久国产精品人妻蜜桃| 男女做爰动态图高潮gif福利片| 午夜老司机福利片| 欧美日韩黄片免| 最近最新中文字幕大全免费视频| 国产精品野战在线观看| 亚洲,欧美精品.| 亚洲欧美精品综合久久99| 日韩欧美 国产精品| 亚洲国产欧美一区二区综合| 色哟哟哟哟哟哟| 亚洲aⅴ乱码一区二区在线播放 | 午夜免费激情av| 欧美绝顶高潮抽搐喷水| 天天添夜夜摸| 99国产精品99久久久久| www.999成人在线观看| 亚洲男人的天堂狠狠| 成人国产一区最新在线观看| 国产av又大| 激情在线观看视频在线高清| 亚洲九九香蕉| 夜夜躁狠狠躁天天躁| 欧美在线黄色| 国产亚洲欧美98| 精品一区二区三区视频在线观看免费| 1024视频免费在线观看| 老司机在亚洲福利影院| 国产成人精品久久二区二区91| 日韩免费av在线播放| 亚洲精品国产精品久久久不卡| 99热这里只有精品一区 | 日韩精品免费视频一区二区三区| 久久久久久亚洲精品国产蜜桃av| 久久精品人妻少妇| 精品欧美国产一区二区三| 成年版毛片免费区| 亚洲一码二码三码区别大吗| 免费在线观看视频国产中文字幕亚洲| 亚洲狠狠婷婷综合久久图片| 日日摸夜夜添夜夜添小说| 日韩 欧美 亚洲 中文字幕| 热99re8久久精品国产| 午夜两性在线视频| 18美女黄网站色大片免费观看| 一本综合久久免费| 国产精品精品国产色婷婷| 啦啦啦韩国在线观看视频| 香蕉av资源在线| 国产单亲对白刺激| 亚洲美女黄片视频| 两人在一起打扑克的视频| 人妻夜夜爽99麻豆av| av超薄肉色丝袜交足视频| 亚洲一区二区三区不卡视频| 国产精品精品国产色婷婷| 国产精品电影一区二区三区| 久久久水蜜桃国产精品网| 成人18禁在线播放| 国产免费av片在线观看野外av| 亚洲免费av在线视频| 男人舔女人的私密视频| 亚洲成人久久爱视频| 国产午夜福利久久久久久| 午夜福利在线观看吧| 欧美绝顶高潮抽搐喷水| 制服丝袜大香蕉在线| 亚洲国产欧美网| 久久香蕉国产精品| 黄色毛片三级朝国网站| 老汉色av国产亚洲站长工具| 国产精品永久免费网站| 少妇被粗大的猛进出69影院| 欧美国产日韩亚洲一区| 97人妻精品一区二区三区麻豆| 99久久国产精品久久久| 嫩草影视91久久| 精品欧美国产一区二区三| 欧美成人午夜精品| 俺也久久电影网| 黄色成人免费大全| 久久精品91无色码中文字幕| 国产精品99久久99久久久不卡| 久久精品国产亚洲av高清一级| 亚洲国产日韩欧美精品在线观看 | 黄片小视频在线播放| 美女扒开内裤让男人捅视频| 午夜福利欧美成人| 校园春色视频在线观看| 日本 av在线| 老司机福利观看| 亚洲成人精品中文字幕电影| 久久草成人影院| 精品人妻1区二区| 国产成人一区二区三区免费视频网站| 波多野结衣高清无吗| 亚洲成人久久爱视频| 少妇粗大呻吟视频| 一进一出抽搐动态| 久久精品人妻少妇| 老司机午夜十八禁免费视频| 女人高潮潮喷娇喘18禁视频| 欧美色欧美亚洲另类二区| 国产三级中文精品| xxxwww97欧美| 亚洲国产精品999在线| av片东京热男人的天堂| 日韩欧美精品v在线| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 一本久久中文字幕| 丰满人妻熟妇乱又伦精品不卡| 一区二区三区高清视频在线| 舔av片在线| 欧美成人免费av一区二区三区| 全区人妻精品视频| 香蕉丝袜av| 脱女人内裤的视频| 久久久国产欧美日韩av| 欧美日韩一级在线毛片| 久热爱精品视频在线9| 亚洲九九香蕉| 99热这里只有是精品50| 亚洲美女黄片视频| 免费人成视频x8x8入口观看| 黄色a级毛片大全视频| 中文字幕人成人乱码亚洲影| 女警被强在线播放| 亚洲午夜精品一区,二区,三区| 90打野战视频偷拍视频| 欧美最黄视频在线播放免费| 男人舔女人下体高潮全视频| 亚洲免费av在线视频| 精品一区二区三区视频在线观看免费| 18禁美女被吸乳视频| 国产激情欧美一区二区| 18禁黄网站禁片免费观看直播| tocl精华| 亚洲午夜精品一区,二区,三区| 久久精品综合一区二区三区| 久久国产精品人妻蜜桃| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 99久久久亚洲精品蜜臀av| 免费搜索国产男女视频| 欧美日韩国产亚洲二区| 国产黄a三级三级三级人| 精品久久久久久久久久久久久| 在线视频色国产色| 免费在线观看成人毛片| 日韩欧美 国产精品| 午夜精品在线福利| 国产精品99久久99久久久不卡| 在线观看美女被高潮喷水网站 | 亚洲专区字幕在线| 老汉色av国产亚洲站长工具| 久久伊人香网站| 啦啦啦免费观看视频1| 日本一区二区免费在线视频| 亚洲成人免费电影在线观看| 午夜福利在线观看吧| 欧美日本视频| 亚洲国产看品久久| 成人av在线播放网站| 日本熟妇午夜| 久久欧美精品欧美久久欧美| 欧美乱妇无乱码| 黄色成人免费大全| 母亲3免费完整高清在线观看| 一本综合久久免费| 国产一区二区三区在线臀色熟女| 淫妇啪啪啪对白视频| 久久婷婷人人爽人人干人人爱| 黄色女人牲交| 免费在线观看黄色视频的| 国语自产精品视频在线第100页| 色综合站精品国产| 国产成人av激情在线播放| 亚洲片人在线观看| 欧美人与性动交α欧美精品济南到| 国产伦一二天堂av在线观看| 亚洲成人国产一区在线观看| 性欧美人与动物交配| 久久久久亚洲av毛片大全| 国产熟女午夜一区二区三区| 全区人妻精品视频| 欧美黑人巨大hd| 少妇熟女aⅴ在线视频| 最近视频中文字幕2019在线8| 国产99久久九九免费精品| 日韩 欧美 亚洲 中文字幕| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 中文字幕人妻丝袜一区二区| 麻豆av在线久日| 国产激情欧美一区二区| 久久中文看片网| 亚洲av电影不卡..在线观看| 久久亚洲精品不卡| 男女那种视频在线观看| 欧美性猛交黑人性爽| 亚洲一区二区三区不卡视频| 中亚洲国语对白在线视频| 亚洲精品在线美女| 精品一区二区三区四区五区乱码| 国产精品久久久av美女十八| 国产人伦9x9x在线观看| 亚洲精品美女久久av网站| 亚洲午夜理论影院| 夜夜爽天天搞| 日本黄大片高清| x7x7x7水蜜桃| 日本 av在线| 国产精品一区二区免费欧美| 亚洲无线在线观看| 国内精品一区二区在线观看| 十八禁人妻一区二区| 亚洲成人久久性| 亚洲成人久久爱视频| 久久国产乱子伦精品免费另类| 欧洲精品卡2卡3卡4卡5卡区| 哪里可以看免费的av片| 熟妇人妻久久中文字幕3abv| 18禁裸乳无遮挡免费网站照片| 悠悠久久av| 国产区一区二久久| 色综合亚洲欧美另类图片| 久久人妻av系列| videosex国产| 亚洲一区中文字幕在线| 一本一本综合久久| 久久人人精品亚洲av| 国产成人啪精品午夜网站| 日本精品一区二区三区蜜桃| 色综合婷婷激情| 99国产精品一区二区蜜桃av| 亚洲精品中文字幕在线视频| 女人被狂操c到高潮| 久久99热这里只有精品18| 亚洲欧美日韩高清在线视频| 少妇裸体淫交视频免费看高清 | 19禁男女啪啪无遮挡网站| 在线观看舔阴道视频| 一级片免费观看大全| 黄色 视频免费看| 免费在线观看亚洲国产| 黄色女人牲交| 99国产精品99久久久久| 免费人成视频x8x8入口观看| 国产高清视频在线播放一区| 婷婷丁香在线五月| 欧美黄色片欧美黄色片| 一级毛片高清免费大全| 欧美在线一区亚洲| 少妇粗大呻吟视频| 老司机福利观看| 国产精品一区二区三区四区久久| 韩国av一区二区三区四区| 色播亚洲综合网| 久久婷婷成人综合色麻豆| 免费在线观看成人毛片| 国产精品久久久人人做人人爽| 丁香六月欧美| 国产精品久久久久久人妻精品电影| 精品国产亚洲在线| 欧美丝袜亚洲另类 | 欧美日韩亚洲综合一区二区三区_| 欧美精品亚洲一区二区| 欧美日韩一级在线毛片| avwww免费| 成人欧美大片| 人成视频在线观看免费观看| 久久久久久人人人人人| 熟女少妇亚洲综合色aaa.| 老汉色av国产亚洲站长工具| 成人高潮视频无遮挡免费网站| 午夜精品在线福利| 亚洲国产精品sss在线观看| 久久国产精品人妻蜜桃| 我要搜黄色片| 夜夜看夜夜爽夜夜摸| 美女 人体艺术 gogo| 国产精品av视频在线免费观看| 黄色毛片三级朝国网站| 欧美黑人巨大hd| 精品一区二区三区视频在线观看免费| 欧美中文日本在线观看视频| 久久亚洲精品不卡| 麻豆一二三区av精品| 啪啪无遮挡十八禁网站| 伊人久久大香线蕉亚洲五| 亚洲avbb在线观看| 国产精品九九99| 日韩精品中文字幕看吧| 操出白浆在线播放| 婷婷六月久久综合丁香| 亚洲国产欧美网| 一边摸一边抽搐一进一小说| 人人妻人人澡欧美一区二区| 黑人欧美特级aaaaaa片| 香蕉丝袜av| 999久久久精品免费观看国产| 国产单亲对白刺激| 免费高清视频大片| 久久久久国产精品人妻aⅴ院| 后天国语完整版免费观看| 婷婷精品国产亚洲av在线| 欧美中文综合在线视频| 90打野战视频偷拍视频| 欧美成人午夜精品| 性欧美人与动物交配| 久久久久久久久久黄片| 在线十欧美十亚洲十日本专区| 久9热在线精品视频| av在线播放免费不卡| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 亚洲成人久久性| 欧美成人免费av一区二区三区| 亚洲av成人av| 亚洲色图av天堂| 精品无人区乱码1区二区| 香蕉av资源在线| 好看av亚洲va欧美ⅴa在| 久久九九热精品免费| 啦啦啦免费观看视频1| 免费在线观看日本一区| 亚洲国产精品成人综合色| 在线a可以看的网站| 亚洲成a人片在线一区二区| 99精品欧美一区二区三区四区| 真人做人爱边吃奶动态| 久久天躁狠狠躁夜夜2o2o| 久久久久亚洲av毛片大全| √禁漫天堂资源中文www| 午夜福利免费观看在线| 夜夜看夜夜爽夜夜摸| 中出人妻视频一区二区| 变态另类丝袜制服| 老司机在亚洲福利影院| 日本撒尿小便嘘嘘汇集6| 亚洲成人国产一区在线观看| 欧美 亚洲 国产 日韩一| 在线观看日韩欧美| 国产又色又爽无遮挡免费看| 成人一区二区视频在线观看| 亚洲精华国产精华精| 国产三级在线视频| 中文字幕最新亚洲高清| 久久人妻福利社区极品人妻图片| 精品国产超薄肉色丝袜足j| 亚洲免费av在线视频| www.999成人在线观看| 久久久久国产精品人妻aⅴ院| 美女 人体艺术 gogo| 1024香蕉在线观看| 亚洲欧美日韩东京热| 午夜激情福利司机影院| 欧美黑人欧美精品刺激| 黑人操中国人逼视频| 欧美zozozo另类| 亚洲中文字幕日韩| 18禁黄网站禁片午夜丰满| 男女午夜视频在线观看| 99久久综合精品五月天人人| 黄色女人牲交| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 国产精品99久久99久久久不卡| 久久九九热精品免费| 2021天堂中文幕一二区在线观|