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

        ?

        土壤濕度對(duì)中國(guó)南部熱帶森林土壤甲烷吸收的影響

        2022-08-03 01:54:54李君怡趙俊福
        生態(tài)學(xué)報(bào) 2022年12期
        關(guān)鍵詞:森林土壤山腳土壤濕度

        李君怡,席 毅,趙俊福

        1 北京大學(xué)城市與環(huán)境學(xué)院, 北京 100871 2 海南省生態(tài)環(huán)境監(jiān)測(cè)中心, 海口 571126

        甲烷(CH4)是僅次于二氧化碳(CO2)的第二大溫室氣體,盡管它在大氣中的濃度只有CO2的0.45%,其增溫潛勢(shì)卻是CO2的28—34倍[1]。在過去的幾十年,大氣CH4的濃度不斷上升,到2020年,已經(jīng)達(dá)到1873 mm3/m3[2]。大氣CH4濃度的增長(zhǎng)是CH4排放與吸收平衡之后的結(jié)果。根據(jù)最新的全球甲烷排放收支(2008—2017)[3],全球CH4年排放量約為576(550—594)Tg/a,主要來源于水稻田、濕地、海洋以及化石燃料的開采與燃燒。全球CH4年吸收量約為556(501—574)Tg/a,主要通過在對(duì)流層中與羥基自由基發(fā)生氧化還原反應(yīng)、向平流層傳輸、在通氣條件下被土壤中的甲烷氧化菌氧化消除[4]。土壤產(chǎn)甲烷過程是缺氧環(huán)境中有機(jī)質(zhì)分解的核心環(huán)節(jié),而甲烷氧化是缺氧-有氧界面的重要微生物過程[5]。既往研究顯示,土壤甲烷通量存在較大的時(shí)空差異[6],氣候[7]、植被類型[8—9]、底質(zhì)可利用性[10]等影響CH4產(chǎn)生、氧化、傳輸過程的因子都將引起CH4通量的變化。氣體擴(kuò)散率隨土壤充氣孔隙的增加而增加,影響土壤的通氣狀況且控制土壤與大氣間的氣體交換。在溫帶森林[11—15]和溫帶草地[16—18]的實(shí)地研究中,觀測(cè)到CH4通量與土壤溫度、土壤濕度的季節(jié)動(dòng)態(tài)呈現(xiàn)一致性。然而,對(duì)于熱帶地區(qū)的土壤CH4通量與環(huán)境因子之間的關(guān)系卻鮮有研究。

        熱帶森林土壤是大氣CH4一個(gè)重要的匯[19],在全球甲烷收支中發(fā)揮著重要的作用。一項(xiàng)基于318個(gè)站點(diǎn)(包含62個(gè)熱帶森林站點(diǎn))的全球土壤甲烷吸收的薈萃分析表明,熱帶森林土壤的CH4吸收量約為4.65 Tg CH4-C/a,約占全球土壤CH4吸收量的28%[20]。另有研究表明,熱帶森林土壤CH4的吸收量在近幾十年正在逐步增加[21]。然而,由于氣候類型的多樣性及森林結(jié)構(gòu)的復(fù)雜性,有關(guān)熱帶森林土壤CH4吸收對(duì)全球變化的響應(yīng)及機(jī)制的認(rèn)識(shí)仍十分有限。相較于寒帶森林[22—24]和溫帶森林[25—26],目前熱帶尤其是東南亞的CH4通量觀測(cè)較少。東南亞的熱帶森林土壤是甲烷的源還是匯尚不清楚,其土壤CH4通量的量級(jí)大小、季節(jié)動(dòng)態(tài)、驅(qū)動(dòng)因子亟需進(jìn)一步的研究。

        吊羅山國(guó)家森林公園位于我國(guó)海南省東南部,是我國(guó)重要的熱帶森林保護(hù)區(qū)之一。本研究通過對(duì)吊羅山熱帶森林土壤原位甲烷通量及相關(guān)環(huán)境因子連續(xù)兩年(2016年9月至2018年9月)的觀測(cè),旨在分析熱帶森林土壤甲烷通量的季節(jié)變化規(guī)律及其與環(huán)境因子的關(guān)系,研究結(jié)果將為東南亞地區(qū)長(zhǎng)時(shí)間甲烷通量觀測(cè)及熱帶森林土壤甲烷匯的估算提供數(shù)據(jù)支持。

        1 材料與方法

        1.1 研究區(qū)概況

        本研究選取海南省吊羅山國(guó)家森林公園(18°43′—18°58′N,109°43′—110°03′E)為研究區(qū)域(圖1的紅色五角星處)。該公園位于海南省陵水縣、保亭縣和瓊中縣交界處,總面積3.8萬hm2,是我國(guó)重要的熱帶森林保護(hù)區(qū)之一。氣候?qū)贌釒ШQ笮詺夂?全年暖熱,降水豐沛。年降水量1870—2760 mm,年均氣溫24.4 ℃,最冷月平均氣溫15.4 ℃,月均相對(duì)空氣濕度80%—85%[27]。

        圖1 研究區(qū)地理位置Fig.1 The geographical location of the study area

        由于吊羅山公園地形起伏較大,山頂和山腳分別是典型的山地森林和低地森林,因此本研究分別在山頂(海拔為970 m)和山腳(海拔為260 m)布設(shè)兩個(gè)樣地。山腳的母巖為黑云母花崗巖,土壤類型以山地赤紅壤為主,主要組成樹種有青梅(Vaticamangachapoi)、鈍葉新木姜子(Neolitseaobtusifolia)、光葉巴豆(Crotonlaevigatus)、子凌蒲桃(Syzygiumtephrodes)、瓊欖(Gonocaryumlobbianum)和九節(jié)(Psychotriarubra)等。山頂?shù)哪笌r為花崗巖,土壤類型主要為山地黃壤,主要組成樹種有陸均松(Dacrydiumpierrei)、叢花山礬(Symplocospoilanei)、楔葉柃(Euryacuneate)、五列木(Pentaphylaxeuryoides)、紅鱗蒲桃(Syzygiumhancei)和脈葉虎皮楠(Daphniphyllumpaxianum)等[27]。位于山頂?shù)臒釒降厣趾臀挥谏侥_的熱帶低地森林均為20世紀(jì)70—80年代擇伐后天然更新的次生林,林下基本無草本層,兩類森林的土壤基本理化性質(zhì)如表1所示。地表常年有凋落物覆蓋,山頂樣地的凋落物多于山腳樣地。凋落物的主要組成是葉片,樹枝及果實(shí)較少。

        表1 熱帶山地森林和熱帶低地森林的土壤基本理化性質(zhì)Table 1 Soil physical and chemical properties in tropical lowland and montane forest

        1.2 樣品采集與處理

        在研究區(qū)山頂和山腳進(jìn)行CH4通量、CO2通量及主要環(huán)境因子包括土壤溫度和土壤濕度的測(cè)量,于2016年9月至2018年9月每月測(cè)量1次。受地形因素的影響,設(shè)置山頂樣地大小為50 m × 50 m,山腳樣地的大小為30 m × 30 m。樣地的地理位置如圖1所示。在每個(gè)樣地中,分別沿對(duì)角線選取三個(gè)樣方(10 m × 10 m),布設(shè)土壤呼吸環(huán)以測(cè)量土壤的甲烷通量和二氧化碳通量。在每一個(gè)樣方中,隨機(jī)布設(shè)五個(gè)土壤呼吸環(huán)(土壤呼吸環(huán)插入土壤中的深度為10 cm)。土壤呼吸環(huán)的外直徑20 cm,長(zhǎng)10 cm。五個(gè)土壤呼吸環(huán)測(cè)得通量的平均值代表該樣方的數(shù)值。

        甲烷通量和二氧化碳通量均采用便攜式溫室氣體分析儀(LGR Inc., San Jose, USA)測(cè)量,其測(cè)量原理為離軸積分輸出腔光譜技術(shù)。本研究采取原位測(cè)量,測(cè)量前不移除凋落物層。距土壤表面約5 cm處的土壤溫度和土壤濕度(全土體體積含水量)由Decagon 5 TM傳感器(Decagon Devices Inc., USA)測(cè)量。為了盡可能減小人為測(cè)量造成的數(shù)據(jù)偏差,每次溫室氣體通量和環(huán)境因子的測(cè)量工作均在上午進(jìn)行,此時(shí)的土壤甲烷通量接近全天的平均值[28]。每個(gè)呼吸環(huán)測(cè)量1次的時(shí)間是170 s,一共測(cè)量三次,取其平均值作為當(dāng)次測(cè)量的結(jié)果。在本研究中,從大氣到土壤的通量記為負(fù)值,反之記為正值。

        在升尺度估算全球熱帶森林土壤甲烷通量時(shí),使用的土壤水?dāng)?shù)據(jù)為SMOS衛(wèi)星遙感數(shù)據(jù),來源于SMOS數(shù)據(jù)處理中心(CATDS,https://www.catds.fr/),其有效載荷為使用孔徑合成(MIRAS)儀器的微波成像輻射計(jì),是全球第一個(gè)采用L波段的遙感器。由于L波段能夠穿透大氣和植被,因此有效地排除了其對(duì)于土壤水分反演的干擾[29]。

        1.3 數(shù)據(jù)處理

        采用MATLAB R2018b和RStudio軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。采用Pearson相關(guān)分析對(duì)CH4通量、CO2通量和土壤溫度、土壤濕度進(jìn)行相關(guān)分析。

        相關(guān)系數(shù)的計(jì)算公式如下:

        (1)

        由于土壤中CH4吸收隨著土壤濕度的增加并非呈現(xiàn)出單調(diào)遞減的關(guān)系,而是當(dāng)土壤濕度非常大甚至土壤趨于飽和時(shí),土壤由CH4的匯變成CH4的源。因此,本研究使用分段線性回歸[30]擬合土壤CH4吸收和土壤濕度的關(guān)系:

        (2)

        式中,y是甲烷通量,x是土壤濕度,α是CH4通量變化趨勢(shì)開始改變的斷點(diǎn),β0,β1,β2是回歸系數(shù),ε是殘差。利用RStudio(http://rstudio.com/)中的segmented包進(jìn)行分段回歸,并進(jìn)行斷點(diǎn)估計(jì)。

        運(yùn)用MATLAB R2018b軟件進(jìn)行數(shù)據(jù)制圖。

        2 結(jié)果與分析

        2.1 土壤溫度和土壤濕度的季節(jié)變化

        圖2展示了2016年9月至2018年9月研究區(qū)表層土壤(5 cm)的溫度和濕度的季節(jié)變化。與氣溫的季節(jié)變化趨勢(shì)一致,土壤溫度的最大值出現(xiàn)在夏季(山頂樣地為24.9 ℃,山腳樣地為27.5 ℃),最小值出現(xiàn)在冬季(山頂樣地為11.5 ℃,山腳樣地為15.4 ℃)。山頂樣地的年平均土壤溫度比山腳樣地低3.3 ℃。對(duì)于土壤濕度,山頂樣地年平均土壤濕度為19.2%,年內(nèi)的波動(dòng)較小(2.8%)。而山腳樣地的年平均濕度相對(duì)較低,為12.7%,且年內(nèi)波動(dòng)大(5.4%)。

        圖2 2016—2018年土壤溫度和土壤濕度的季節(jié)動(dòng)態(tài)Fig.2 Seasonal variations of soil temperature and soil moisture from 2016 to 2018誤差線表示標(biāo)準(zhǔn)誤

        2.2 二氧化碳通量和甲烷通量的季節(jié)變化及相關(guān)性

        圖3展示了研究時(shí)段研究區(qū)CO2通量和CH4通量的季節(jié)變化。山頂樣地和山腳樣地均為CO2的源(即土壤呼吸過程),山頂樣地年平均CO2排放通量約為1.72 μmol m-2s-1,山腳樣地約為2.36 μmol m-2s-1,且在4—10月CO2的排放量更多。對(duì)于CH4通量,山頂樣地和山腳樣地在研究時(shí)段內(nèi)基本都為CH4的匯(山頂樣地年平均CH4吸收通量約為0.25 nmol m-2s-1,山腳樣地約為0.51 nmol m-2s-1)。相比于CO2,山頂和山腳CH4吸收的季節(jié)波動(dòng)表現(xiàn)出更明顯的分歧,山頂CH4吸收的季節(jié)波動(dòng)比山腳低0.15 nmol m-2s-1。

        圖3 2016—2018年二氧化碳通量和甲烷通量的季節(jié)動(dòng)態(tài)Fig.3 Seasonal variations of CO2 flux and CH4 flux from 2016 to 2018誤差線表示標(biāo)準(zhǔn)誤

        圖4展示了研究時(shí)段CH4吸收通量與CO2排放通量的關(guān)系,CO2排放通量?jī)H能解釋CH4吸收通量月際變化的0.2%,二者之間的相關(guān)性并不顯著(P=0.76)。

        圖4 2016—2018年甲烷通量對(duì)二氧化碳通量的響應(yīng)Fig.4 Response of CH4 flux to CO2 flux from 2016 to 2018

        2.3 二氧化碳通量和甲烷通量對(duì)土壤溫度和土壤濕度的響應(yīng)

        圖5展示了研究時(shí)段CO2排放和CH4吸收與土壤溫度和濕度的關(guān)系。對(duì)于CO2的排放,研究時(shí)段CO2的排放與土壤溫度呈現(xiàn)出顯著的正相關(guān)(相關(guān)系數(shù)R=0.80,P<0.01),土壤溫度可以解釋研究區(qū)CO2排放月際變化的64%,而土壤濕度和CO2的排放之間沒有顯著的相關(guān)關(guān)系。相比之下,CH4的吸收與環(huán)境因子的關(guān)系與之相反。研究時(shí)段CH4的吸收和土壤濕度之間呈現(xiàn)出顯著的負(fù)相關(guān)關(guān)系(相關(guān)系數(shù)R=0.88,P<0.01),但當(dāng)土壤濕度超過20%時(shí),土壤CH4的吸收基本接近于0,并存在轉(zhuǎn)化為CH4的源的現(xiàn)象。因此,對(duì)于土壤濕度和CH4吸收的關(guān)系,我們采用分段函數(shù)進(jìn)行擬合:

        圖5 2016—2018年二氧化碳通量和甲烷通量對(duì)土壤溫度和土壤濕度的響應(yīng)Fig.5 Responses of CO2 flux and CH4 flux to soil temperature and soil moisture from 2016 to 2018

        (3)

        2.4 熱帶森林土壤甲烷吸收量升尺度估算

        由于本研究中發(fā)現(xiàn)熱帶森林土壤濕度對(duì)于甲烷吸收有明顯的主導(dǎo)作用,使用SMOS衛(wèi)星觀測(cè)的土壤濕度數(shù)據(jù)對(duì)熱帶森林土壤甲烷通量進(jìn)行了升尺度估算。圖6展示了2017年熱帶地區(qū)土壤濕度和甲烷吸收通量的空間格局。綜合本研究及表2中研究的結(jié)果,得到熱帶森林土壤濕度和CH4吸收的關(guān)系:

        圖6 2017年熱帶地區(qū)土壤濕度以及甲烷吸收通量的空間格局圖Fig.6 Spatial pattern of soil moisture and CH4 uptake flux in tropics in 2017

        (4)

        另外,由于研究站點(diǎn)測(cè)量的土壤濕度受到儀器和土壤質(zhì)地的影響,測(cè)量值比衛(wèi)星觀測(cè)偏低。因此在升尺度估算前,對(duì)測(cè)量數(shù)據(jù)進(jìn)行了偏差校正。經(jīng)過偏差校正之后升尺度估算到的2017年熱帶森林土壤(1.35×109hm2,來源于MODIS MCD12Q1土地覆被數(shù)據(jù))的甲烷吸收通量約為 1.72 Tg CH4-C/a。

        3 討論

        3.1 熱帶森林土壤作為甲烷的匯

        土壤甲烷通量是甲烷產(chǎn)生過程(主要在厭氧條件下發(fā)生)和甲烷氧化過程(主要在有氧條件下發(fā)生)平衡之后的結(jié)果[31]。由于熱帶森林土壤一般為不飽和的土壤,因此更有利于發(fā)生甲烷的氧化作用。本研究中,無論是山頂樣地還是山腳樣地,均表現(xiàn)為甲烷的匯。熱帶森林土壤表現(xiàn)為甲烷的匯這一結(jié)論已經(jīng)被之前的大量研究證實(shí)(表2)。王全成等在中國(guó)福建省(27°03′N)的研究發(fā)現(xiàn)杉木人工林土壤全年均表現(xiàn)為甲烷匯(甲烷年通量為1.86 kg CH4-C hm-2a-1)[32]。Gong等在中國(guó)廣東省(21°27′N)持續(xù)一年的通量觀測(cè)表明,熱帶森林土壤為甲烷匯(甲烷年通量為0.33 kg CH4-C hm-2a-1)[33]。Davidson等在巴西(02°54′S)的研究也發(fā)現(xiàn),熱帶森林土壤表現(xiàn)為甲烷的吸收,年甲烷吸收通量為0.83 kg CH4-C hm-2a-1[34]。Zhao等搜集了全球54項(xiàng)熱帶森林土壤研究中的甲烷通量數(shù)據(jù),其中94%的甲烷通量為負(fù)值,也證實(shí)了熱帶森林土壤是較強(qiáng)的甲烷吸收匯[35]。

        表2 熱帶不同區(qū)域森林土壤甲烷年通量Table 2 Annual CH4 flux of forest soils in different regions of the tropics

        季節(jié)上,本研究發(fā)現(xiàn)甲烷吸收通量在干季(11月—次年4月)明顯高于濕季(5—10月)(約0.28 nmol m-2s-1)。Keller和Reiners的研究也發(fā)現(xiàn),哥斯達(dá)黎加(10°26′N)的老齡林中,甲烷吸收通量在干季更高(約0.61 nmol m-2s-1)[36]。賈高輝等在海南島尖峰嶺(18°23′N)熱帶山地雨林的研究也表明,土壤甲烷通量整體呈現(xiàn)出旱季大于雨季(約0.26 nmol m-2s-1)[28]。莫江明等在鼎湖山(23°10′N)的研究也發(fā)現(xiàn),土壤甲烷的吸收速率在研究期間出現(xiàn)夏季最低(約0.09 nmol m-2s-1),冬季最高(約0.93 nmol m-2s-1)[37]。

        3.2 土壤溫度和濕度對(duì)熱帶森林土壤甲烷吸收的影響

        土壤濕度決定了土壤的通氣情況,是影響甲烷吸收和排放的重要環(huán)境因子[40]。本研究中,無論是在熱帶山地森林還是熱帶低地森林,甲烷吸收和土壤濕度均呈現(xiàn)出顯著的負(fù)相關(guān),且土壤濕度可以解釋甲烷吸收變化的至少76%(P<0.01),顯示出熱帶森林土壤濕度對(duì)甲烷吸收的主導(dǎo)作用。一致的結(jié)論在熱帶其他地區(qū)的森林土壤也有證實(shí)(表2)。Fang等在中國(guó)云南省西雙版納(21°56′N)的研究發(fā)現(xiàn),土壤濕度可以解釋甲烷吸收變化的37%(P<0.01)[38]。Jones等在秘魯(13°11′S)的研究同樣發(fā)現(xiàn),甲烷吸收通量與土壤濕度呈顯著負(fù)相關(guān)(R=0.79,P<0.05)[39]。Rowlings等在澳大利亞(26°45′S)的熱帶森林的實(shí)驗(yàn)表明,土壤濕度控制甲烷吸收(R=0.52,P<0.01)[41]。這些實(shí)驗(yàn)都表明,低的土壤濕度決定的有氧環(huán)境是甲烷氧化過程強(qiáng)于甲烷吸收過程的重要前提條件。

        土壤溫度的變化可以影響甲烷產(chǎn)生和氧化反應(yīng)中酶的活性,進(jìn)而影響甲烷的排放和吸收[19]。然而,本研究中,無論是山頂樣地還是山腳樣地,土壤溫度和甲烷吸收通量之間均沒有顯著的相關(guān)關(guān)系,但是土壤溫度明顯影響了土壤呼吸的過程。這意味著,在研究區(qū),土壤濕度對(duì)于甲烷通量的變化是占主導(dǎo)作用的,土壤溫度的影響則非常輕微。土壤溫度對(duì)于熱帶森林土壤甲烷通量的弱的影響在其他研究中也有體現(xiàn)。例如,Werner等在中國(guó)云南省西雙版納(21°57′N)原始林的研究發(fā)現(xiàn),土壤溫度僅能解釋甲烷吸收變化的5%(P<0.05)[42]。Werner等在肯尼亞西部(00°16′N)熱帶森林實(shí)驗(yàn)中也發(fā)現(xiàn),土壤溫度與甲烷吸收通量無顯著相關(guān)性(R2<0.05,P<0.01)[43]。

        3.3 熱帶森林土壤甲烷吸收量的估算

        本研究經(jīng)過偏差校正之后升尺度估算到的2017年熱帶(30°N—30°S)森林土壤(1.35×109hm2)的甲烷吸收通量約為 1.72 Tg CH4-C/a,而Dutaur和Verchot運(yùn)用薈萃分析搜集了熱帶森林土壤甲烷吸收的58個(gè)測(cè)量值,估算出熱帶森林(1.93×109hm2)土壤的甲烷吸收通量為4.65 Tg CH4-C/a[20]。由于Dutaur和Verchot的研究使用了比本研究中更大的森林面積,且不同站點(diǎn)直接平均受到站點(diǎn)不具備代表性的影響(其熱帶站點(diǎn)大多集中在南美洲北部,缺乏非洲、亞洲的實(shí)測(cè)數(shù)據(jù)),因此其估算的結(jié)果可能存在較大的不確定性。除了基于站點(diǎn)測(cè)量的估算結(jié)果, Curry使用CLASS(Canadian Land Surface Scheme)模型模擬的熱帶森林(1.55×109hm2)土壤的甲烷吸收通量為3.11 Tg CH4-C/a[21],由于該模型使用了更大的熱帶森林面積(2×108hm2),且模型的關(guān)鍵參數(shù)β和k0的驗(yàn)證使用的觀測(cè)站點(diǎn)較少,因此模型模擬結(jié)果也存在較大的不確定性??傮w上,大尺度的熱帶森林土壤甲烷吸收的估算仍然需要更多實(shí)地觀測(cè)的結(jié)果。

        4 結(jié)論

        本研究通過測(cè)量2016—2018年吊羅山熱帶森林土壤甲烷通量和相關(guān)的環(huán)境因子,發(fā)現(xiàn)我國(guó)海南省吊羅山森林公園的熱帶森林土壤表現(xiàn)為甲烷的匯,且干季的甲烷吸收通量明顯高于濕季。土壤濕度主導(dǎo)了研究區(qū)甲烷吸收的月際變化。由于熱帶地區(qū)的甲烷地面觀測(cè)站點(diǎn)分布稀疏,且大部分站點(diǎn)觀測(cè)時(shí)長(zhǎng)有限,更多長(zhǎng)期的實(shí)地觀測(cè)研究可以幫助了解熱帶森林土壤甲烷通量的季節(jié)變化特征,提升對(duì)其在全球甲烷收支中作用的認(rèn)識(shí)。

        猜你喜歡
        森林土壤山腳土壤濕度
        天元智能干選機(jī)在山腳樹礦平田矸石山的應(yīng)用
        桂源鋪
        土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
        基于51單片機(jī)控制花盆土壤濕度
        電子制作(2019年15期)2019-08-27 01:12:12
        濃霧
        抵達(dá)山腳
        四川盆地土壤濕度時(shí)空分布及影響因子分析
        中國(guó)不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
        兩種亞熱帶森林土壤甲烷氧化活性的垂直分布特征
        龍灣自然保護(hù)區(qū)森林土壤動(dòng)物群落多樣性及功能類群
        一区二区三区蜜桃在线视频| 亚洲精品无码久久久久牙蜜区| 婷婷丁香五月中文字幕| 狠狠丁香激情久久综合| 亚洲黑寡妇黄色一级片| 亚洲国产高清精品在线| 人妻少妇精品视频无码专区| 国产精品视频久久久久| 人妻少妇无乱码中文字幕| 久久精品av在线观看| gv天堂gv无码男同在线观看| 亚洲va在线va天堂va手机| 日韩精品人妻少妇一区二区| 在线观看亚洲av每日更新影片| 极品嫩模高潮叫床| 亚洲日本天堂| 黄片午夜免费观看视频国产 | 亚洲AV色欲色欲WWW| 亚洲禁区一区二区三区天美| 综合色就爱涩涩涩综合婷婷| 国产乱沈阳女人高潮乱叫老| 国产免费午夜福利蜜芽无码| 国产黄久色一区2区三区| 成人无码网www在线观看| 欧美日韩成人在线| 综合激情中文字幕一区二区| 日本二一三区免费在线| 国产麻豆md传媒视频| 色噜噜狠狠色综合中文字幕| 亚洲女同性恋第二区av| 在线观看的网站| 欧美日本日韩aⅴ在线视频| 蜜桃网站在线免费观看视频| 国产精品专区第一页天堂2019| 日韩一卡2卡3卡4卡新区亚洲| 国产精品无码久久久一区蜜臀 | 高清中文字幕一区二区三区| 岛国av无码免费无禁网站| 在线观看免费午夜大片| av在线网站手机播放| 麻豆资源在线观看视频|