999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

晝夜增溫速率的不對稱性對新疆地區植被動態的影響

2017-06-28 12:36:36劉雪佳杜自強武志濤徐小明山西大學黃土高原研究所山西太原030006山西大學環境與資源學院山西太原030006
中國環境科學 2017年6期
關鍵詞:生長影響

趙 杰,劉雪佳,杜自強*,武志濤,徐小明(.山西大學黃土高原研究所,山西 太原 030006;.山西大學環境與資源學院,山西 太原 030006)

晝夜增溫速率的不對稱性對新疆地區植被動態的影響

趙 杰1,劉雪佳2,杜自強1*,武志濤1,徐小明1(1.山西大學黃土高原研究所,山西 太原 030006;2.山西大學環境與資源學院,山西 太原 030006)

結合1982~2013年歸一化植被指數(NDVI)、植被類型和氣象數據,分別利用一元線性回歸和二階偏相關分析方法,分析了新疆地區生長季白天最高氣溫(Tmax)和生長季夜間最低氣溫(Tmin)的變化趨勢以及晝夜增溫的不對稱性對各類型植被活動的影響.結果表明:近 32年來新疆地區Tmax及Tmin都存在極為顯著的上升趨勢,生長季Tmin的增長速率約為Tmax增長速率的1.5倍,晝夜增溫速率呈現出不對稱性;晝夜增溫對植被活動的影響呈現出明顯的分異特征:相比于白天增溫,夜間溫度的變化對新疆植被的影響區域更為廣泛,而且對植被的影響程度更大;各植被類型 NDVI對晝夜增溫速率的不對稱性也產生了不同的響應:白天增溫對針葉林產生顯著的積極影響,而夜間增溫對針葉林、農業植被、草原及草甸植被表現出顯著的積極影響.

不對稱性增溫;植被動態;白天增溫;夜間增溫;新疆地區

植被作為陸地生態系統的重要組成部分,是連接大氣、土壤、水體的自然紐帶,在全球變化與陸地生態系統的研究中充當著“指示器”的作用[1-2].氣候變化是植被活動的重要影響因素,溫度、降水等作為植物生長發育的必要環境因子,對植被的生長和物候等具有重要影響[3].目前,國內外在植被動態變化及其與氣候因子的相互關系方面進行了大量的研究,得到不同的研究區域和不同的時空尺度下植被生長對氣溫和降水的時空響應特征[3-9].

過去50年的氣溫數據表明,全球地表晚間升溫增速約為白天的 1.4倍[10].由于大部分地表植物的光合作用在白天進行,然而植物呼吸作用貫穿全天.因此,這種不均衡的變暖速率必然會影響植被的碳吸收和碳消耗[11].目前,國內外在植被動態變化及其與氣候因子相互關系方面的研究大多集中于植被對平均氣溫變化的響應;考慮晝夜增溫速率不一致對植被生長變化影響的研究相對較少.Wan等[12]通過控制實驗研究了內蒙古草原日最高氣溫、夜間最低氣溫與歸一化植被指數(NDVI)的關系,發現日最高溫的增加抑制了該地區草地樣地NDVI的升高,而夜間最低溫的升高則促進了草地樣地NDVI的增加.Peng等[10]利用遙感數據、大氣CO2濃度觀測數據以及氣象數據,系統地研究了白天和晚上溫度上升對北半球植被生產力和生態系統碳源匯功能影響,發現白天溫度的升高有利于大部分寒帶和溫帶濕潤地區植被生長及其生態系統碳匯功能,但并不利于溫帶干旱和半干旱地區植被生長.而在晚上,溫度上升對植被生長的影響相反.梁春玲等[13]利用氣象數據和遙感數據計算了近40年來晝夜增溫對南四湖濕地植被 NDVI 的影響差異,發現NDVI與白天氣溫呈現極顯著的偏相關關系,而夜間氣溫與 NDVI 的偏相關關系在任何時期內均不顯著.可見,不同地區或者不同類型植被對晝夜溫度變化的響應具有明顯的差異.

新疆(E73°40′~96°33',N34°25′~49°10′)地處中國中溫帶極端干旱的荒漠地帶,自然環境復雜、氣候類型各異、水熱分布懸殊,這為從局地尺度驗證更大尺度的植被與氣候變化相互關系提供了良好的實驗場所[14].

本文利用近32年新一代NDVI數據集以及中國植被分類數據、白天和夜間極端氣溫數據,分析新疆地區晝夜氣溫、植被NDVI生長季變化及對晝夜不對稱增溫響應的差異.以期豐富全球氣候變化背景下,晝夜溫度上升對中國內陸干旱、半干旱地帶生態系統植被活動影響的認識,為不同植被類型對氣候變化的響應特征研究提供借鑒.

1 數據和方法

1.1 數據

1.1.1 遙感數據 1982~2013年 NDVI數據是由美國國家航空航天局的全球觀測模擬與制圖研究組(Global Inventory Modeling and Mapping Studies, GIMMS)提供的第3代NOAA/AVHRR遙 感 數 據(NDVI3g)(http://www.landcover.org/ data/gimms),其空間分辨率為8km,時間分辨率為15d.該數據集消除了太陽高度角、傳感器靈敏度隨時間變化等影響,并結合交叉輻射定標的方法,增強了數據的精度,已被廣泛應用于大區域與全球尺度植被覆蓋研究中[15-17].本文采用最大值合成法重建生長季各月最大化NDVI數據集[18],并利用《中國植被編碼1:400萬植被圖集》,排除無植被地區.

1.1.2 氣象數據 氣象數據來源于中國氣象科學數據共享服務網(http://cdc.cma.gov.cn)的中國地面氣候資料月值數據集,包括極端最高氣溫、極端最低氣溫、降水量等.采用反距離權重法對各個氣象因子進行空間插值[19],得到逐像元的氣象數據.將極端最高溫度作為Tmax,將極端最低氣溫作為Tmin來處理[10].

1.1.3 植被類型數據 植被類型數據來自于寒區旱區科學數據中心(http://westdc.westgis.ac.cn)《中國植被編碼1:400萬植被圖集》.該數據共有4個大類,13個綱組,50個群系綱,110個群系組.本文將新疆地區植被類型劃分為7種類型:針葉林、闊葉林、灌從(含萌生矮林)、荒漠、草原(含稀疏灌木草原)、草甸(含草木沼澤)和農業植被.

1.2 研究方法

1.2.1 氣候要素趨勢分析 采用最小二乘法擬合氣候要素(Tmax、Tmin)隨時間的變化速率,即氣候要素y隨時間t變化的線性回歸系數,可采用一元一次方程表示:

式中:t為時間;常量a和系數b由最小二乘法求得;b為回歸常數項;a為回歸系數,a值的符號及數值大小反映 y上升或下降的趨勢.線性回歸法可以清楚地顯示所給出的序列是否存在某種明顯的變化.

1.2.2 偏相關分析法 采用二階偏相關分析法,消除其他變量的干擾來研究晝夜增溫的不對稱性對植被NDVI的影響.例如,通過限制Tmin和降水量計算生長季NDVI與Tmax的偏相關系數.同樣,限制Tmax和降水量計算生長季NDVI和Tmin的偏相關系數.二階偏相關系數通過一階偏相關系數計算得到,而計算一階偏相關系數需要首先計算相關系數.相關系數的計算公式為[20]:

一階偏相關系數的計算公式為:

二階偏相關系數的計算公式為:

式中:x,y為需要進行偏相關系數計算的要素;1,2為控制變量.

偏相關系數的顯著性檢驗,一般采用 t 檢驗法.t檢驗法的統計量計算公式為[21]:

式中:r為偏相關系數,n為樣本數,q為自由度個數.

2 結果與分析

2.1 新疆地區晝夜增溫趨勢分析

1982 年以來,新疆地區Tmin和Tmax都有著不同程度的升高趨勢(圖 1).從擬合方程來看,近 32年來新疆地區Tmin和Tmax平均每10年分別增加0.642℃(P<0.01)和 0.458℃(P<0.01),晝夜溫差每10年減小0.18℃(P<0.1).Tmax和Tmin二者之間存在極為顯著的正相關性(P<0.01).這說明,新疆地區生長季 Tmin的增加趨勢遠快于 Tmax的增加趨勢,晝夜增溫速率呈現不對稱性.晝夜溫差(TmaxTmin)呈現較為顯著的減小態勢(P<0.1).

圖1 1982~2013年生長季新疆地區Tmax和Tmin變化Fig.1 Variations of Tmaxand Tminin growing seasons in Xinjiang during 1982~2013

2.2 新疆地區植被NDVI與晝夜增溫的偏相關分析

由于新疆屬于干旱、半干旱地區,降水為植被生長的主要限制因子之一.因此,在進行偏相關分析時,將降水作為控制變量處理.即通過控制降水和Tmin計算NDVI與Tmax的偏相關系數;通過控制降水和Tmax計算NDVI與Tmin的偏相關系數.通過逐像元計算NDVI與氣溫變量的偏相關系數后,得到偏相關系數的空間格局(圖 2).從像元數量上來看,去除 Tmin、降水的影響后,約有55.14%的地區NDVI同Tmax之間呈現正相關關系,其中7.39%的地區表現為顯著正相關(P<0.05).同樣,約有44.86%的地區呈現負相關,其中3.99%的地區表現為顯著負相關(P<0.05).然而,當去除Tmax、降水量的影響后,約有 58.78%的地區植被生長季NDVI和Tmin呈現負相關,其中顯著負相關的地區達到44.00%(P<0.05).同樣,約有41.22%的地區,生長季 NDVI與 Tmin呈現正相關,其中,約有 25.03%的地區表現為顯著正相關(P<0.05).從通過顯著性檢驗的區域占新疆總面積的比例來看,生長季NDVI與Tmin的偏相關性通過顯著性檢驗的地區占新疆總面積的36.18%,而生長季NDVI與Tmax的偏相關性通過顯著性檢驗的比例僅為 5.87%.可見,夜間溫度的變化對植被影響的區域更為廣泛.從受其影響地域的空間分布來看,除新疆西北部的塔城地區、克拉瑪依市、博爾塔拉蒙古自治州和哈密地區及新疆南部的巴音郭楞蒙古自治州中部與和田地區西南部外,大部分區域植被 NDVI與 Tmax呈現正相關性;生長季NDVI與Tmin偏相關性在新疆地區西北部至西南部大多呈現為正相關,東北部至東南部大多呈現為負相關.

圖2 1982~2013年生長季新疆地區NDVI與Tmax和Tmin的偏相關系數Fig.2 Partial correlation coefficients between growing-season NDVI and Tmax/ Tminin Xingjiang during 1982~2013

2.3 新疆地區各植被類型NDVI與晝夜增溫的偏相關分析

表1 不同植被類型NDVI與Tmax及Tmin的偏相關系數Table 1 Partial correlation coefficients between NDVI and Tmax/ Tminfor different vegetation types

生長季不同植被類型NDVI與Tmax和Tmin的偏相關系數(表 1)顯示,除農業植被外,其他各類植被NDVI與Tmax整體上皆表現為偏正相關性,但僅有針葉林呈現為顯著相關(P<0.05).而各類型植被對 Tmin的偏相關性差異較為明顯,除荒漠外,其它各類植被NDVI與Tmin皆呈現為偏正相關性.針葉林、農業植被、草原及草甸與 Tmin的偏相關系數最高,其中前三類植被的NDVI與Tmin表現為極顯著相關(P<0.01),后者亦表現為顯著相關(P<0.05).與 Tmax偏負相關的農業植被和與 Tmin偏負相關的荒漠植被類型,其相關關系皆未通過顯著性檢驗(P>0.05).總體上,晝夜氣溫的升高對新疆地區各類型植被產生了積極作用.

為研究不同類型的植被 NDVI與晝夜氣溫偏相關性的空間分布差異,得到偏相關性通過統計顯著性檢驗(P<0.05)的像元占該植被類型總像元的比例(表 2).可以看出,各類植被 NDVI與Tmax的偏相關性通過顯著性檢驗的比例普遍較低,其中針葉林的比例最高,但僅為 13.75%.相比之下,NDVI與 Tmin的偏相關性通過顯著性檢驗的比例普遍相對較高.其中,農業植被(38.26%)和針葉林(33.75%)NDVI與Tmin顯著正相關的比例最高;荒漠和闊葉林NDVI與Tmin顯著負相關的比例最高,分別為34.29%和30.49%.

綜合來看(表 1,表 2),除灌叢和荒漠植被外,相比于白天增溫,夜間增溫對針葉林、闊葉林、草原、草甸和農業植被的影響程度更大,影響區域更廣.

表2 通過統計學顯著性檢驗(P<0.05)的像元占該植被類型總像元的比例(%)Table 2 Percentages of statistically significant (P<0.05) pixels in different vegetation types (%)

3 討論

本文顯示近32年間新疆地區生長季夜間氣溫升溫速率約為白天氣溫升溫速率的 1.5倍.Solomon等[22]在研究過去50年全球地表溫度的晝夜增溫時發現,夜間最低氣溫的升溫速率約為日最高氣溫升溫速率的 1.4倍.這說明新疆地區與全球具有一致的晝夜變暖趨勢,晝夜增溫存在不對稱性.樸世龍等[10]在研究北半球植被對晝夜增溫速率不對稱性的響應時發現,在北方干旱、半干旱地區生長季植被NDVI與Tmin表現為明顯的正相關關系.本文同樣認為位于干旱、半干旱地區的新疆植被 NDVI對夜間氣溫變化的響應更為顯著.梁春玲等[13]在研究晝夜增溫對南四湖濕地植被 NDVI的影響差異時發現,NDVI與 Tmax呈現極顯著的偏相關關系,而本文發現Tmax僅對針葉林的影響顯著,對其他類型植被的影響并不顯著;梁春玲等[13]認為Tmin與NDVI的偏相關關系并不顯著,而本文發現 Tmin對多數植被類型的影響顯著.發生這種差異的原因可能與溫度對水份條件充足的濕地植被和干旱、半干旱的荒漠地帶植被的影響機理不同有關.

根據植物學理論,植物的光合作用通常在白天發生,在夜間一般不能進行光合作用.因此,晝夜溫度對植被的作用機理不同.夜間氣溫的升高引起植物呼吸作用的增強和體內養分代謝速率的增加,造成植物體內有機質損耗增加[13],降低灌漿期持續時間,減小胚乳細胞體積,從而會降低植被生產力[23-24].這可能會使某些植被 NDVI與Tmin產生負的偏相關關系,比如荒漠植被.另一方面,夜間氣溫的升高也可能會通過降低霜凍發生的頻率[25]以及夜間自養呼吸加劇而產生的補償作用[12]來提高植被生產能力.這可能使得某些植被NDVI與Tmin之間呈現出正的偏相關關系,比如針葉林、草原、草甸等大多數類型植被.相對于 Tmin,Tmax的上升對新疆地區植被的影響作用并不十分顯著;研究區地處干旱、半干旱地區,植被生長主要受制于降水是其可能的原因之一.

植被 NDVI的時空演變是自然因素和人為因素綜合作用的結果,本文只分析了晝夜溫度、降水量與植被變化的相關性,而未能將太陽輻射量等其他的氣候因子和人為因素考慮在內.此外,NDVI在稀疏植被區的低敏感性可能會對結果產生一定的影響.所以,需要在以后的研究中綜合各種因素,采用控制試驗或模型模擬研究等多種方法相互補充,以進一步闡釋晝夜增溫對不同植被的影響機理.

4 結論

4.1 新疆地區1982~2013年間晝夜增溫趨勢顯著;晝夜增溫速率具有不對稱性;晝夜溫差呈現減小的趨勢.

4.2 新疆地區晝夜增溫對植被的影響呈現出明顯的空間差異:總體上,相比于白天增溫,夜間增溫對植被活動的影響程度更大,影響區域更為廣泛.

4.3 新疆地區晝夜增溫對不同類型植被會產生不同的影響.從影響的程度來看,白天增溫對針葉林的影響顯著;夜間增溫對針葉林、農業植被、草原及草甸的影響顯著.從影響的空間分布來看,夜間溫度對不同類型植被有顯著影響的面積比例均高于白天溫度.

[1]袁麗華,蔣衛國,申文明,等. 2000~2010黃河流域植被覆蓋的時空變化 [J]. 生態學報, 2013,33(24):7798-7806.

[2]李 霞,李曉兵,王 宏,等.氣候變化對中國北方溫帶草原植被的影響 [J]. 北京師范大學學報(自然科學版), 2006,42(6):618-623.

[3]Zhu Z, Piao S, Myneni R B, et al. Greening of the Earth and its drivers[J]. Nature climate change, 2016,6(8):791-795.

[4]毛德華,王宗明,宋開山,等.東北多年凍土區植被 NDVI變化及其對氣候變化和土地覆被變化的響應 [J]. 中國環境科學, 2011,31(2):283-292.

[5]Fensholt R, Proud S R. Evaluation of earth observation based global long term vegetation trends-Comparing GIMMS and MODIS global NDVI time series [J]. Remote Sensing of Environment, 2012,119:131-147.

[6]崔林麗,史 軍,楊引明,等.中國東部植被 NDVI對氣溫和降水的旬響應特征 [J]. 地理學報, 2009,64(7):850-860.

[7]Jiapaer G., Liang S, Yi Q, et al. Vegetation dynamics and responses to recent climate change in Xinjiang using leaf area index as an indicator [J]. Ecological Indicators, 2015,58:64-76.

[8]王 強,張 勃,戴聲佩,等.三北防護林工程區植被覆蓋變化與影響因子分析 [J]. 中國環境科學, 2012,32(7):1302-1308.

[9]趙安周,劉憲鋒,朱秀芳,等. 2000~2014年黃土高原植被覆蓋時空變化特征及其歸因 [J]. 中國環境科學, 2016,36(5):1568-1578.

[10]Peng S S, Piao S L,Ciais P,et al. Asymmetric effects of daytime and nighttime warming on Northern Hemisphere vegetation [J]. Nature, 2013,501:88-94.

[11]Atkin O K, Turnbull M H, Zaragoza-Castells J, et al. Light inhibition of leaf respiration as soil fertility declines along a post-glacial chronosequence in New Zealand: an analysis using the Kok method [J]. Plant and Soil, 2013,367(1):163-182.

[12]Wan S Q, Xia J, Liu W, et al. Photosynthetic overcompensation under nocturnal warming enhances grassland carbon sequestration [J]. Ecology, 2009,90:2700-2710.

[13]梁春玲,于泉洲,劉煜杰,等.晝夜增溫對南四湖濕地植被 NDVI的影響差異 [J]. 應用生態學報, 2010,12(3):315-321.

[14]龐 靜,杜自強,張霄羽.新疆地區植被對水熱條件的時滯響應[J]. 中國農業資源與區劃, 2015,36(7):82-88.

[15]張景華,封志明,姜魯光,等.瀾滄江流域植被 NDVI與氣候因子的相關性分析 [J]. 自然資源學報, 2015,30(9):1425-1435.

[16]Kobayashi H, Dye D G. Atmospheric conditions for monitoring the long-term vegetation dynamics in the Amazon using normalized difference vegetation index [J]. Remote Sensing of Environment, 2005,97(4):519-525.

[17]張生軍,王 濤,王天明,等.新疆不同植被 NDVI的變化及其與氣候因子的關系 [J]. 草業科學, 2009,26(5):26-31.

[18]郭 鈮,朱燕君,王介民,等.近22年來西北不同類型植被NDVI變化與氣候因子的關系 [J]. 植物生態學報, 2008,32(2):319-327.

[19]何 月,樊高峰,張小偉,等.浙江省植被 NDVI動態及其對氣候的響應 [J]. 生態學報, 2012,32(14):4352-4362.

[20]張勇榮,周忠發,馬士彬,等.基于 NDVI的喀斯特地區植被對氣候變化的響應研究—以貴州省六盤水市為例 [J]. 2014,34(4):114-117.

[21]侯光磊,張洪巖,郭 聃,等.長白山植被生長季 NDVI時空變化及其對氣候因子敏感性 [J]. 地理科學進展, 2012,31(3):285-292.

[22]Solomon S, Qin D, Manning M, et al. Climate Change 2007: The physical science basis. Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change (Cambridge University Press, 2007).

[23]Tan J G, Piao S L, Chen A P, et al. Seasonally different response of photosynthetic activity to daytime and night-time warming in the Northern Hemisphere [J]. Global Change Biology, 2015,21(1):377-87.

[24]Welch J R, Vincent J R, Auffhammer M, et al. Rice yields in tropical/subtropical Asia exhibit large but opposing sensitivities to minimum and maximum temperatures [J]. Proceedings of the National Academy of Sciences, 2010,107(33):14562-14567.

[25]Lobell D B. Changes in diurnal temperature range and national cereal yields [J]. Agricultural and forest meteorology, 2007,145(3):229-238.

Effects of the asymmetric diurnal-warming on vegetation dynamics in Xinjiang.

ZHAO Jie1, LIU Xue-jia2, DUZi-qiang1*, WU Zhi-tao1, XU Xiao-ming1(1.Institute of Loess Plateau, Shanxi University, Taiyuan 030006, China;2.College of Environmental & Resource Science, Shanxi University, Taiyuan 030006, China). China Environment Science, 2017,37(6):2316~2321

Based on normalized difference vegetation index (NDVI) data, vegetation type data and meteorological data, this study revealed the warming trends of the day- and nighttime in growing seasons using the unitary linear regression analysis, and analyzed different effects of the asymmetry diurnal warming on vegetation activities using the two order partial correlation analysis method. Results indicated that: there was an extremely significant trend of diurnal warming from 1982 to 2013 in Xinjiang areas. The nighttime warming was approximately 1.5times that of daytime during growing seasons. The warming rate showed an asymmetry between day- and nighttime; the effects of day- and nighttime warming on vegetation activities showed obvious differentiation. Nighttime warming had significantly stronger and more extensive effects on vegetation than daytime warming; and different vegetation types responded differently to the asymmetry warming. Daytime warming had significantly positive effects on the coniferous forest, and nighttime warming had significant positive effects on the coniferous forest, agricultural vegetation, grassland, and meadow.

asymmetric warming;vegetation dynamics;daytime warming;nighttime warming;Xinjiang areas

X173

A

1000-6923(2017)06-2316-06

趙 杰(1993-),男,山東臨沂人,碩士研究生,主要從事全球變化生態學、植被與生態遙感方面的研究工作.發表論文5篇.

2016-11-17

國家自然科學基金資助項目(41161066,41401643, 41401053)

* 責任作者, 副教授, duzq@sxu.edu.cn

猜你喜歡
生長影響
是什么影響了滑動摩擦力的大小
哪些顧慮影響擔當?
當代陜西(2021年2期)2021-03-29 07:41:24
碗蓮生長記
小讀者(2021年2期)2021-03-29 05:03:48
共享出行不再“野蠻生長”
生長在哪里的啟示
華人時刊(2019年13期)2019-11-17 14:59:54
野蠻生長
NBA特刊(2018年21期)2018-11-24 02:48:04
生長
文苑(2018年22期)2018-11-19 02:54:14
沒錯,痛經有時也會影響懷孕
媽媽寶寶(2017年3期)2017-02-21 01:22:28
《生長在春天》
擴鏈劑聯用對PETG擴鏈反應與流變性能的影響
中國塑料(2016年3期)2016-06-15 20:30:00
主站蜘蛛池模板: 中日韩一区二区三区中文免费视频 | 一级毛片在线播放| 久久国产免费观看| 婷婷色一区二区三区| 久草美女视频| 无码专区国产精品第一页| 午夜欧美在线| 亚洲国产欧美国产综合久久 | 成人自拍视频在线观看| 久久福利网| 日本草草视频在线观看| 美女高潮全身流白浆福利区| 色九九视频| 91福利国产成人精品导航| 国产在线自在拍91精品黑人| 亚洲欧美日韩中文字幕一区二区三区| 青青草原国产精品啪啪视频| 亚洲成A人V欧美综合天堂| 日本亚洲欧美在线| 不卡网亚洲无码| 国产精品视频a| 国产一区亚洲一区| 成人精品午夜福利在线播放| 日本在线欧美在线| 欧美一级一级做性视频| 亚洲视频a| 中日韩欧亚无码视频| 亚洲国产精品日韩av专区| 欧美一区福利| 福利小视频在线播放| 九九久久精品免费观看| 波多野结衣无码视频在线观看| 亚洲天堂精品在线观看| 久久综合亚洲色一区二区三区| 国内精品视频| 日本一区二区三区精品AⅤ| 国产成人综合日韩精品无码首页| 欧美亚洲日韩不卡在线在线观看| 在线色国产| 久久鸭综合久久国产| 18黑白丝水手服自慰喷水网站| 国产亚洲欧美日韩在线一区二区三区| 国产无码精品在线播放| 亚洲精品视频免费观看| 成人在线视频一区| 色婷婷在线影院| 国内嫩模私拍精品视频| 伊人大杳蕉中文无码| 一本大道视频精品人妻 | 免费大黄网站在线观看| 日日拍夜夜操| 人妻无码一区二区视频| 啦啦啦网站在线观看a毛片| 欧美一级大片在线观看| 国产91久久久久久| 伊人福利视频| 伊人五月丁香综合AⅤ| 亚洲av无码久久无遮挡| 国产日韩AV高潮在线| 狼友av永久网站免费观看| 国产成人无码AV在线播放动漫| 永久免费精品视频| 国产成人亚洲综合a∨婷婷| 18禁黄无遮挡免费动漫网站| 97在线观看视频免费| 91亚洲免费| 亚洲无线一二三四区男男| 亚洲成人福利网站| 国产一级毛片在线| 国产白浆视频| 欧美69视频在线| 国产精品任我爽爆在线播放6080| 欧美一级在线播放| 四虎国产永久在线观看| 91视频区| 精品色综合| 亚洲av无码牛牛影视在线二区| 久久黄色毛片| 国产精品刺激对白在线| 欧美精品另类| 久草视频中文| 国产欧美另类|