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

        ?

        呼倫湖流域近55年氣候變化特征

        2016-11-14 11:50:37吳丹瑞耿欣欣
        水土保持研究 2016年6期

        吳丹瑞, 耿欣欣, 李 明

        (1.鄱陽湖濕地與流域研究教育部重點實驗室(江西師范大學), 南昌 330000;2.西北農(nóng)林科技大學 資源環(huán)境學院, 陜西 楊凌 712100; 3.江蘇南水科技有限公司, 南京 210012)

        ?

        呼倫湖流域近55年氣候變化特征

        吳丹瑞1,2, 耿欣欣3, 李 明1,2

        (1.鄱陽湖濕地與流域研究教育部重點實驗室(江西師范大學), 南昌 330000;2.西北農(nóng)林科技大學 資源環(huán)境學院, 陜西 楊凌 712100; 3.江蘇南水科技有限公司, 南京 210012)

        收集了呼倫湖流域1960—2014年的逐日氣溫、降水、蒸發(fā)和相對濕度數(shù)據(jù),提取了包括年平均、季節(jié)平均值及極端天氣等9個指標,采用Mann-Kendall趨勢檢驗、線性傾向估計法以及5 a滑動平均等方法在年和季節(jié)尺度上分析了該地區(qū)的氣候變化趨勢,并利用Pearson相關系數(shù)探究這四個氣象要素之間的相關性。結(jié)果表明:(1) 年氣溫顯著上升,降水變化不顯著,蒸發(fā)顯著減少,相對濕度顯著降低,氣候向暖干化方向發(fā)展。(2) 春夏秋季溫度、蒸發(fā)和相對濕度的變化趨勢與年相同,最大升溫速率發(fā)生在春季,每十年上升0.43℃,冬季變化均不顯著。(3) 平均降水和極端降水在冬春季節(jié)都顯著增加,然而夏季降水以最大幅度減少,并且夏季持續(xù)時間較長的極端降水減少顯著。(4) 夏季氣溫升高最顯著,相對濕度顯著降低,容易發(fā)生干旱。冬季極端降水增多則可能導致雪災加劇。

        氣候變化; 氣象要素; 呼倫湖流域; 氣溫; 干旱

        全球氣候變化帶來一系列生態(tài)環(huán)境問題,如海平面上升、干旱、極端天氣等等,近些年成為各國關注的熱點問題。政府間氣候變化專門委員會(IPCC)第五次評估報告指出[1],過去130年全球升溫0.85℃,最近110年全球平均海平面上升0.19 m,過去10 a間冰川融化的速度比20世紀90年代加快了數(shù)倍,人類活動極有可能是20世紀中期以來全球氣候變暖的主要原因。在RCP8.5情形下,預計到21世紀末全球表面氣溫變化可能超過2℃,平均降水在中緯度和副熱帶干旱區(qū)可能減少,在中緯度濕潤地區(qū)可能增加[1]。20世紀80年代以來,極端天氣事件如高溫熱浪、暴雨、干旱、冰雹等發(fā)生頻率急劇上升,各種氣象災害也隨之增加,給全球社會經(jīng)濟以及生態(tài)環(huán)境造成了嚴重的損失和破壞[2]。

        呼倫湖是我國第五大淡水湖,也是北方第一大湖,與周圍濕地被稱為呼倫貝爾“草原之腎”,對呼倫貝爾草原的經(jīng)濟發(fā)展和生態(tài)保護起著重要作用[3]。呼倫湖流域位于中高緯度受東南季風影響的邊緣區(qū)域,遠離海洋,大陸性極強,同時又是南北冷暖氣流交匯沖撞頻繁的地區(qū)[4],成為對全球氣候變化的響應最積極,最敏感的地區(qū)之一。近年來,呼倫湖流域春夏連續(xù)干旱以及冬季特大雪災發(fā)生頻率加劇。研究表明,受氣候趨向暖干化的影響,呼倫湖水源補給不足,導致水位下降,湖面萎縮,周圍蘆葦大面積減少,珍稀鳥類大量遷徙,生態(tài)環(huán)境嚴重惡化[5]。目前,針對呼倫湖流域?qū)θ驓夂蜃兓憫难芯砍蔀闊狳c,但是大多只是集中在氣候變化對呼倫湖濕地生態(tài)環(huán)境的影響[3,5-7],對呼倫湖流域的氣候變化特征分析大多以年均氣溫和降水為主要氣象要素,指標選取單一,缺乏對極端天氣的研究,也沒有分析相對濕度的變化。所以,利用連續(xù)多年觀測到的氣象要素資料提取多種指標以不同的時間尺度對呼倫湖流域的氣候變化進行分析研究,以揭示多年來當?shù)貧夂蜃兓厔?,對呼倫湖流域生態(tài)環(huán)境的恢復、重建以及可持續(xù)發(fā)展提供科學的參考依據(jù)。

        1 材料與方法

        1.1 流域概況

        呼倫湖也稱達賚湖[8],位于呼倫貝爾草原的中部,東經(jīng)117°00′10′′—117°41′40′′,北緯48°30′40′′—49°20′40′′,呈不規(guī)則斜長方形。呼倫湖的水源補給除來源于大氣降水和地下水外,主要依靠境外河流克魯倫河,以及連接其“姐妹湖”貝爾湖的烏爾遜河。呼倫湖隨著歷史的構造運動和氣候的變化,其面積也曾多次擴大或縮小。其歷史記錄最大水深8 m,蓄水量達138.5億立方米,由于氣候變的干旱等原因,湖體水位不斷下降。自2000—2009年,水位下降了約4.6 m,湖面縮小了477 km2,蓄水量比歷史最大蓄水量減少了一半。從前呼倫湖最北部與海拉爾河相通,湖水可以外流進入黑龍江,現(xiàn)在已經(jīng)斷流成為內(nèi)陸湖。

        呼倫湖流域?qū)儆谥袦貛О敫珊挡菰瓪夂颍瑲夂蚝涓稍?,冬長夏短,春季干旱多風沙,秋季霜凍早,降溫快。該地區(qū)年均氣溫1.1℃,無霜期110 d,年降雨量集中在6—9月,年均降雨量為200~300 mm左右[6],湖水封凍時間長達180 d,冰層厚可達1 m以上。

        1.2 數(shù)據(jù)獲得

        本文根據(jù)位于呼倫湖流域中部的新巴爾虎右旗的氣象觀測站1960—2014年逐日平均氣溫、最高氣溫、降水、蒸發(fā)量以及相對濕度資料,分別計算了近55 a平均氣溫、最高氣溫和相對濕度的年、季節(jié)平均值以及降水和蒸發(fā)量的年、季節(jié)的總量,另外根據(jù)降水的日數(shù)據(jù)分別提取出年、季節(jié)的極端降水各指標所需數(shù)據(jù)。季節(jié)尺度的數(shù)據(jù)采用1960年3月至2013年2月,規(guī)定春季為3—5月,夏季6—8月,秋季9—11月,冬季12月至次年2月。

        1.3 指標選取

        本文分析氣候變化趨勢的氣象要素有4種,其中氣溫有兩個指標,分別為日平均最高氣溫和日平均氣溫,降水有五個指標,分為平均降水和極端降水,極端降水又有日最大降水量(AM1R)、連續(xù)3日最大降水量(AM3R)、連續(xù)5日最大降水量(AM5R)和連續(xù)7日最大降水量(AM7R),蒸發(fā)指標為總蒸發(fā)量,相對濕度為平均相對濕度。

        1.4 分析方法

        此研究采用Mann-Kendall非參數(shù)統(tǒng)計檢驗對呼倫湖流域近55 a的氣候要素指標變化趨勢進行檢驗,同時應用線性傾向估計、5 a滑動平均分析氣候要素的變化趨勢,此外,還利用SPSS計算出Pearson相關系數(shù)來探究氣象要素間的相關關系。

        Mann-Kendall(簡稱MK)是水文序列中最受歡迎的一種趨勢檢驗方法[9-10],被世界氣象組織(WMO)廣泛推薦應用于時間序列的趨勢分析。許多研究人員利用MK趨勢檢驗方法分析降水、氣溫、泥沙和徑流等的時間序列的變化趨勢。MK檢驗中,原假設H0為各指標時間序列數(shù)據(jù)(x1,x2,…,xn)是獨立同分布的;備擇假設H1是對于所有的k,j≤n,且k≠j,xi與xj具有不同的分布,檢驗統(tǒng)計量S被定義如下:

        (1)

        (2)

        式中:n是數(shù)據(jù)記錄長度,即時間序列長度55;xi和xj是連續(xù)數(shù)據(jù)值。當n≥10時,檢驗統(tǒng)計量S近似服從正態(tài)分布,均值和方差如下:

        (3)

        (4)

        式中:m是固定組的數(shù)量;ti是第i個固定組的大小。我們通過以下計算將檢驗統(tǒng)計量S轉(zhuǎn)化為標準檢驗統(tǒng)計量Z:

        (5)

        MK方法對數(shù)據(jù)有序列獨立的要求。對時間序列進行分析時,序列自相關性的出現(xiàn)會嚴重影響趨勢檢驗的結(jié)果[11]。序列相關為正,可能高估趨勢的顯著性;相關為負,可能導致低估。

        yt=xt-φxt-1

        (6)

        式中:yt是時間序列預白化的值;xt是間隔時間為t的原始序列值;φ是滯后1個時間間隔的自相關系數(shù)。當序列相關在置信水平為5%顯著時,使用預白化處理。我們的指標變量中,大部分沒有顯著的序列自相關性。在進行MK趨勢檢驗之前,首先利用預白化去除時間序列中的自相關。

        利用已經(jīng)提取出來的各項指標序列,以時間t為自變量(年或月),各氣象要素指標為因變量,建立一元回歸方程,方程為下式,即

        y(t)=b0+b1t

        (7)

        dy(t)/dt=b1

        (8)

        式中:b0,b1通過最小二乘法獲得。將b1×10稱為變化傾向率,表示每十年的變化大小。b1為正值,表示在計算時段內(nèi)有增加(上升)趨勢,b1為負值,表示有減少(降低)趨勢,絕對值越大,變化程度越大[12]。

        5 a滑動平均我們采用前五年的平均值作為第五年的值,依次后推。

        2 結(jié)果與分析

        2.1 溫度與蒸發(fā)的年際變化

        近55 a呼倫湖流域的年平均氣溫存在年際間的微小波動,但整體上升趨勢明顯(圖1a)。其變化傾向率為0.37℃/10 a,即55 a氣溫上升了2℃。MK檢驗的標準檢驗統(tǒng)計量Z值為4.54,即年平均氣溫在0.01水平上顯著上升,與眾多前人研究結(jié)果[3,6,13]以及全球變暖的大趨勢一致。從5 a滑動平均可以看出,20世紀60年代初至80年代后期該地區(qū)年平均氣溫呈現(xiàn)微小波動上升,80年代后期至今上升速率加大,年際間波動也較之前大。年際間平均最高氣溫的變化趨勢與平均氣溫相同(圖1b),其變化傾向率為0.27℃/10 a,Z值為3.17,上升趨勢顯著。

        圖1 呼倫湖流域年際氣溫變化

        近55 a研究區(qū)年蒸發(fā)量平均值為2 113.9 mm。MK檢驗的Z值為-2.58,絕對值大于2.32,說明蒸發(fā)量在0.01水平顯著減少,與張娜等的研究結(jié)果相反[6],可能是我們的時間尺度更長的原因。該地區(qū)溫度顯著上升,飽和蒸汽壓增大,理論上蒸發(fā)應該增加,但是由圖2中的5 a滑動平均可知,年總蒸發(fā)量在波動中顯著減少??梢姡摰貐^(qū)蒸發(fā)量的減少與氣溫變化關系不大,影響蒸發(fā)量的因子除氣溫外還包括太陽輻射、風速、云量、水汽壓等[14],如沙塵、霧霾以及其他氣溶膠增加可能會對太陽短波輻射的散射增強,這種衰減效應會導致蒸發(fā)量減少。

        圖2 呼倫湖流域年際蒸發(fā)變化曲線

        2.2 降水與濕度的年際變化

        近55 a呼倫湖流域年際變化呈現(xiàn)微小減少趨勢(圖3),年均降水量243.3 mm。Z值為-0.73,沒有達到0.01的顯著性水平,說明降水量減少不顯著,與趙慧穎等的研究結(jié)果一致[3]。由5 a滑動平均可以看出,20世紀60年代初至70年代初降水量明顯減少,70年代初至90年代中期降水量緩慢波動增加,90年代中期至今降水量減少迅速。

        圖3 呼倫湖流域年際降水變化曲線

        圖4 呼倫湖流域年際濕度變化曲線

        多年來呼倫湖流域平均相對濕度為59.6%。MK檢驗的Z值為-4.19,絕對值大于2.32,所以相對濕度在0.01水平上顯著降低。近55 a相對濕度的變化傾向率為-1.19%/10 a,即55 a降低了6.5%。相對濕度的穩(wěn)定波動降低(圖4)反映了研究區(qū)近55 a來干旱正在逐漸加重,是溫度顯著上升、降水減少的結(jié)果。

        2.3 氣候要素的季節(jié)變化規(guī)律

        由表1可以看出,春夏秋季的溫度、濕度和蒸發(fā)的變化趨勢與冬季差異很大。根據(jù)Z值可知溫度的兩個指標日平均氣溫和日最高氣溫在春夏秋季都顯著上升。日平均氣溫的變化傾向率在春季最大,為0.43℃/10 a,55 a來上升了2.37℃,比年平均氣溫上升快,說明春季為溫度上升做貢獻最大。日最高氣溫的變化傾向率在夏季最大,并且Z值超過2.32,說明夏季日最高氣溫在0.01水平顯著上升,極端高溫天氣在夏季發(fā)生的強度增大,促進干旱的產(chǎn)生。冬季的溫度變化不明顯,有微小的上升趨勢,Z值沒有達到0.05的顯著水平。蒸發(fā)在春夏秋都有減少趨勢,其中夏季和秋季均達到了顯著水平,夏季蒸發(fā)減少最迅猛,變化傾向率幾乎是春季或秋季的五倍,與之相反,冬季蒸發(fā)呈不顯著的增加趨勢。相對濕度在四季都下降,除冬季外,其他季節(jié)都顯著下降,其中夏季下降最迅速,依次為秋季、春季。

        表1 不同季節(jié)溫度、濕度及蒸發(fā)的線性傾向估計及M-K分析結(jié)果

        注:“*”,“**”分別代表通過0.05,0.01水平顯著性檢驗。

        降水是呼倫湖流域的敏感氣象要素,降水減少是導致干旱的最主要原因。在年降水量變化不顯著的前提下,呼倫湖流域的冬春季降水顯著增加(表2),秋季稍微增加,夏季降水減少,雖然夏季減少趨勢沒有達到0.05的顯著性水平,但是其變化傾向率的絕對值最大,是冬春秋季降水增加的二倍還要多。極端降水在冬春季同樣顯著增加,春季除AM3R外都達到0.05的顯著性水平,冬季幾乎所有指標都達到0.01的顯著性水平,夏季持續(xù)時間較長的極端降水顯著減少,秋季沒有明顯變化。

        呼倫湖流域?qū)儆谥芯暥劝敫珊祬^(qū),季節(jié)尺度對全球氣候變化的響應是夏季溫度大幅度上升同時伴隨著降水減少,冬季增溫不明顯同時平均降水和極端降水均顯著增加。夏季越來越干旱,冬季出現(xiàn)雪災的頻率可能加劇。春季溫度升高的同時降水增多,可能對作物生長有利。秋季溫度升高、相對濕度降低的同時降水沒有變化,呈現(xiàn)出干燥的趨勢。總的來說,氣候趨向暖干化。

        2.4 各要素間的相關關系

        本文利用呼倫湖流域降水、氣溫和相對濕度的年均值以及年總蒸發(fā)量四個氣象要素的時間序列進行相關性分析。由表3可以看出,相對濕度同降水和氣溫的相關系數(shù)分別是0.509和-0.391,根據(jù)t檢驗,相對濕度同降水和氣溫的相關系數(shù)都通過了α=0.01的顯著性檢驗,說明相對濕度同降水量之間呈正相關,同氣溫呈負相關,即降水量少的年份相對濕度低,氣溫越高相對濕度越低。蒸發(fā)同降水和相對濕度的相關性系數(shù)分別為0.596,0.744,均在0.01水平顯著相關,說明蒸發(fā)量大的年份相對濕度也高,由于此地區(qū)相對濕度整體較小,因此,相對濕度不是蒸發(fā)量的限制因子,反倒隨著蒸發(fā)量的增加,相對濕度會略微增大。蒸發(fā)和氣溫的相關性系數(shù)是0.015,沒有達到顯著相關。

        表2 不同季節(jié)平均降雨和極端降雨指標的線性傾向估計及M-K分析結(jié)果

        注:“*”,“**”分別代表通過0.05,0.01水平顯著性檢驗。

        表3 各氣象要素的Pearson相關性分析結(jié)果

        注:“**”代表在0.01水平(雙側(cè))顯著性相關。

        呼倫湖流域近55 a降水、氣溫、相對濕度和蒸發(fā)的相關性與其他地區(qū)有些不同,甚至完全相反[15],特別是溫度高的年份相對濕度低,蒸發(fā)和溫度沒有顯著相關性。研究表明潛在蒸散量與平均風速、平均氣溫、日照時數(shù)呈正相關,與相對濕度呈負相關,風速的下降是蒸發(fā)量減少的主要原因[16]。而Stahill等[17]則認為是大氣中浮塵、霾等氣溶膠的增加導致了有效輻射量的減少,蒸發(fā)量減少。由此可知,此研究區(qū)溫度上升,氣溶膠粒子上升,削弱太陽輻射,減少蒸發(fā)量。研究區(qū)屬于半干旱地區(qū),水是限制因子,蒸發(fā)量隨著溫度的升高而減小,但是溫度升高使得飽和水汽壓增大,造成相對濕度下降。

        3 結(jié) 論

        (1) 近55 a(1960—2014年)呼倫湖流域平均氣溫和最高溫均在0.01水平上顯著上升,降水量變化不顯著,相對濕度顯著降低,氣候向暖干化發(fā)展。

        (2) 氣候變化在不同季節(jié)趨勢有所不同。春、夏和秋季在日均溫、日最高溫以及相對濕度這三個指標的變化趨勢相同,均表現(xiàn)為溫度顯著上升、相對濕度顯著下降,最大升溫速率發(fā)生在春季,日平均氣溫每十年上升0.43℃,其他指標夏季變化傾向率最大,這三個指標在冬季的變化趨勢都不顯著。蒸發(fā)在春夏秋季為減少趨勢,并且夏秋季達到顯著水平,冬季為不顯著的增加趨勢。無論平均降水還是極端降水冬春季節(jié)均表現(xiàn)為顯著的增多趨勢,相反,夏秋季節(jié)為不顯著的減少趨勢,其中夏季歷時較長的極端降水(AM5R和AM7R)減少顯著。

        (3) 呼倫湖流域近55 a干旱加重主要是因為氣溫升高,飽和水汽壓增大,然而全年降水卻沒有增加,蒸發(fā)不足,相對濕度顯著降低,干旱持續(xù)加重。夏季氣溫升高顯著,降水減少幅度最大,相對濕度顯著降低,相比于其他季節(jié),最易變干旱,而冬季雖然降水增多,極端降水也增多了,雪災可能也隨之增多。

        [1] Intergovernmental panel on climate change. climate change 2013: The physical science basis:Working group I contribution to the fifth assessment report of the intergovernmental panel on climate change [M]. Cambridge:Cambridge University Press, 2014.

        [2] 萬仕全.中國降水與溫度極值的時空分布規(guī)律模擬[D].蘭州:蘭州大學,2010.

        [3] 趙慧穎,李成才,趙恒和,等.呼倫湖濕地氣候變化及其對水環(huán)境的影響[J].冰川凍土,2007,29(5):795-801.

        [4] 蘭玉坤.內(nèi)蒙古地區(qū)近50年氣候變化特征研究[D].北京:中國農(nóng)業(yè)科學院,2007.

        [5] 趙慧穎,烏力吉,郝文俊.氣候變化對呼倫湖濕地及其周邊地區(qū)生態(tài)環(huán)境演變的影響[J].生態(tài)學報,2008,28(3):1064-1071.

        [6] 張娜,烏力吉,劉松濤,等.呼倫湖地區(qū)氣候變化特征及其對湖泊面積的影響[J].干旱區(qū)資源與環(huán)境,2015,29(7):192-197.[7] 白美蘭,郝潤全,沈建國.近46a氣候變化對呼倫湖區(qū)域生態(tài)環(huán)境的影響[J].中國沙漠,2008,28(1):101-107.

        [8] 王素慧,梁宏偉,楊玉生.呼倫湖濕地水環(huán)境治理對策[J].內(nèi)蒙古水利,2006,105(1):35-40.

        [9] Mann HB. Nonparametric tests against trend [J]. Econometrica:Journal of the Econometric Society, 1945,13(3):245-259.

        [10] Kendall MG. Rank [M]. Correlation Measures, London:Charles Griffin.1975.

        [11] Hu Y, Maskey S, Uhlenbrook S. Trends in temperature and rainfall extremes in the Yellow River source region, China[J]. Climatic Change, 2012,110(1/2):403-429.

        [12] 陸志華,夏自強,于嵐嵐,等.松花江流域年降水和四季降水變化特征分析[J].水文,2012,32(2):62-71.

        [13] 顧潤源,李思慧,趙慧穎,等.呼倫湖流域徑流對氣候變化的響應[J].生態(tài)學雜志,2012,31(6):1517-1524.

        [14] 鄭廣芬,陳曉光,孫銀川,等.寧夏氣溫,降水,蒸發(fā)的變化及其對氣候變暖的響應[J].氣象科學,2006,26(4):412-421.

        [15] 吉磊,何新林,劉兵,等.近50a瑪納斯河流域上游氣候變化的分析[J].長江科學院院報,2014,31(12):21-27.

        [16] 李敏敏,延軍平.“蒸發(fā)悖論”在北方農(nóng)牧交錯帶的探討[J].資源科學,2013,35(11):2298-2307.

        [17] Stanhill G, Cohen S. Global dimming:a review of the evidence for a widespread and significant reduction in global radiation with discussion of its probable causes and possible agricultural consequences[J]. Agricultural and forest meteorology, 2001,107(4):255-278.

        Characteristics of Climatic Changes in Hulun Lake Basin in Recent 55 Years

        WU Danrui1,2, GENG Xinxin3, LI Ming1,2

        (1.KeyLaboratoryofPoyangLakeWetlandandWatershedResearch,MinistryofEducation,JiangxiNormalUniversity,Nanchang330000,China; 2.CollegeofNaturalResourcesandEnvironment,NorthwestA&FUniversity,Yangling,Shaanxi712100,China; 3.JiangsuNaiwchCo.,Ltd,Nanjing210012,China)

        The changes of nine meteorological indicators on temperature, precipitation, evaporation and relative humidity from 1960 to 2014 in Hulun Lake basin were analyzed at the scales of year and season. The results show that: (1) annual temperature significantly increased but evaporation and relative humidity decreased significantly, the variation of annual precipitation was insignificant, these results indicated that the climate in this area became warmer and drier; (2) the variation of temperature, evaporation and relative humidity in spring, summer and autumn are similar to their interannual variations, the maximum rate of mean temperature rise occurred in spring with 0.43℃/decade, however, variations of the three meteorological elements were not significant in winter; (3) the mean precipitation and extreme precipitation increased significantly in both winter and spring, but in summer the precipitation decreased, especially the long-time extreme precipitation which decreased significantly; (4) summer climate in this area became much drier because of significant increase in temperature and decrease in relative humidity. The increasing extreme precipitation in winter would cause snow disaster in this area.

        climate change; meteorological element; Hulun Lake Basin; temperature; drought

        2015-11-23

        2015-12-07

        鄱陽湖濕地與流域研究教育部重點實驗室(江西師范大學)開放基金資助項目(PK2014003);國家自然科學基金(51409216)

        吳丹瑞(1992—),女,河南周口人,碩士研究生,研究方向是土壤侵蝕與水土保持。E-mail:wudanrui15@mails.ucas.ac.cn

        李明(1986—),男,甘肅慶陽人,副教授,博士,研究方向是水環(huán)境與生態(tài)水利。E-mail:lileaf@163.com

        P467

        A

        1005-3409(2016)06-0174-05

        99精品视频69v精品视频| 精品丝袜一区二区三区性色| 国产自拍在线视频观看| 后入丝袜美腿在线观看| 国产h视频在线观看| 四虎影院在线观看| 国产人禽杂交18禁网站| 99久久婷婷亚洲综合国产| 无码人妻精品一区二区三区夜夜嗨| 久久久午夜精品福利内容| 91精品国产色综合久久不卡蜜| 国产免费99久久精品| 友田真希中文字幕亚洲| 熟妇人妻av中文字幕老熟妇| 亚洲偷自拍另类图片二区| 国产av精品久久一区二区| 激情人妻另类人妻伦| 亚洲aⅴ天堂av天堂无码麻豆| 国产精品视频一区日韩丝袜| 日本精品少妇一区二区| 日本成本人片视频免费| 人妻激情偷乱一区二区三区| 狠狠色噜噜狠狠狠狠97俺也去| 亚洲丰满熟女乱一区二区三区| 国产精品爽爽久久久久久竹菊| 亚洲欧美日韩精品高清| 国产精品亚洲精品日产久久久| 国产视频一区二区在线免费观看 | 国产乱子伦一区二区三区国色天香| 性饥渴的农村熟妇| 日本一区二区三区高清千人斩| 日本人妻av在线观看| 91九色成人蝌蚪首页| 人妻无码αv中文字幕久久琪琪布| 国产日韩欧美911在线观看| 精品亚洲在线一区二区| 成人午夜福利视频后入| 亚洲精品国产成人AV| 冲田杏梨av天堂一区二区三区| 无遮挡激情视频国产在线观看| 欧美最猛黑人xxxx黑人表情|