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

    貢嘎山海螺溝冰川退縮區(qū)4種常見樹種的異速生長方程

    2019-05-14 11:22:48商宏莉孫守琴何詠梅李安迪
    西南農(nóng)業(yè)學報 2019年4期
    關鍵詞:樣木沙棘樹干

    劉 亭, 商宏莉,羅 輯,孫守琴,何詠梅,李安迪,張 軍

    (1.四川師范大學生命科學學院,四川 成都 610000; 2.中國科學院水利部成都山地災害與環(huán)境研究所山地表生過程和生態(tài)調(diào)控重點實驗室,四川 成都 610041)

    【研究意義】森林生物量是研究群落演替、森林生產(chǎn)力乃至整個生態(tài)系統(tǒng)結構、功能和元素的生物地球化學循環(huán)的重要基礎[1-2]。特別是近幾十年來,隨著氣候變暖等生態(tài)環(huán)境問題的加劇和CO2減排呼聲的高漲,對森林生物量的估算作為生態(tài)學和全球變化研究的重要內(nèi)容之一備受研究者關注[3-6]?!厩叭搜芯窟M展】傳統(tǒng)的森林生物量野外調(diào)查工作量和勞動強度大,嚴重影響了相關工作的進展。隨著數(shù)學模型研究的發(fā)展,各種生物量異速生長模型被開發(fā)并被廣泛用于估算森林生物量[7-11],例如冉飛等[7]對西藏森林的6種林木的生長方程進行了擬合;張繪芳[9]建立了新疆西伯利亞云杉的生物量方程;曾偉生[10]對3種常用方程的優(yōu)劣進行了分析。該方法主要是通過模型將樹木的總生物量及不同組分(如枝、葉、樹干、根等)生物量與胸徑和樹高等易測指標聯(lián)系起來[12],為森林生態(tài)系統(tǒng)生物量和凈生產(chǎn)力的估測提供了一種普遍且可靠的方法。該方法簡單可靠,具有一定的精度保證,同時克服了傳統(tǒng)生物量調(diào)查法費時、費力和破壞性大的缺點,因此成為森林生態(tài)系統(tǒng)生物量和生產(chǎn)力估測的一種普遍可靠的方法?!颈狙芯壳腥朦c】川西亞高山森林是我國第二大林區(qū),在長江上游生態(tài)安全屏障及生態(tài)系統(tǒng)功能維持中起著重要作用。雖然采用異速生長模型預測樹木地上部分生物量目前已有不少報道[13-19],但針對川西亞高山生態(tài)系統(tǒng)樹木生長模型的研究比較少,對地下部分生物量的異速生長模型研究也比較缺乏。例如,僅有的少量針對川西典型針葉樹種峨眉冷杉(Abiesfabri)的研究主要針對人工針葉林[20]或者自然暗針葉林中齡林和成熟林[21]進行,涉及的樹木徑級偏大,難以用于描述植被演替過程中物種的生長動態(tài);而且由于回歸方程中生物量干重數(shù)據(jù)估算的多、實測解析木數(shù)據(jù)偏少,所得的回歸方程在實際應用中還有待進一步修正[21];除此之外,針對川西植被原生演替先鋒樹種沙棘和川滇柳以及演替中期優(yōu)勢樹種冬瓜楊的生物量模型還未見報道,導致該區(qū)植被原生演替及相應植被生態(tài)功能研究進展緩慢?!緮M解決的關鍵問題】本研究以川西貢嘎山海螺溝冰川退縮區(qū)植被原生演替過程中的幾種主要樹種為研究對象,通過整株收獲法研究這些樹種地上地下生物量分配規(guī)律,并分別在整株及不同組分(干、枝、葉、根)水平上建立整株及不同器官生物量異速生長模型,旨在為川西亞高山森林生物量及碳儲量估測提供有效的研究方法,同時為異速生長模型優(yōu)化以及植被原生演替過程研究提供基礎依據(jù)。

    1 研究方法

    1.1 研究區(qū)概況

    本研究在海螺溝(29°34’21’’ N,102°59’42’’ E)冰川退縮區(qū)開展。海螺溝冰川位于四川省甘孜藏族自治州東南部,貢嘎山東坡,是青藏高原東緣的代表性冰川[22]。小冰期以來海螺溝冰川逐步退縮,形成的退縮區(qū)平均海拔約2800 m,氣候濕冷,屬于山地寒溫帶氣候類型[23-24],年平均氣溫5~8 ℃,年降水日數(shù)在260 d以上,年平均降水量約1500~1900 mm。年平均相對濕度約90 %。由于該區(qū)域水熱條件好,冰川消融速度快,近百多年來大約退縮了2 km,且中間沒有冰進過程,土壤為連續(xù)成土過程 。冰川退縮后形成了大量原生裸地,在其長達2 km的序列范圍內(nèi)形成了一個完整的從裸地到先鋒群落、再到云冷杉頂級群落的連續(xù)植被原生演替序列[24-25]。主要樹種有峨眉冷杉(A.fabri)、麥吊云杉(Piceabrachytyla)、冬瓜楊(Populuspurdomii)、糙皮樺(Betulaalbo-sinensisBurk var.septentrionalisSchneid)、沙棘(Hippophaerhamnoides)、川滇柳(Salixrehderiana)以及杜鵑(Rhododendronsimsii)、花楸(Sorbuspohuashanensis)、莢蒾(Viburnumdilatatum)、茶藨子(Ribesnigrum)等灌木。

    1.2 樣木選擇

    由于海螺溝冰川退縮區(qū)1998年以前的區(qū)域為柳樹、沙棘幼樹群落,本研究主要沿1980年到1890年近110年的演替序列進行?;谇叭藢χ脖谎萏骐A段的劃分[25-26],選取5個具有代表性演替階段的4個主要樹種(峨眉冷杉、冬瓜楊、川滇柳和沙棘)進行研究。出于生態(tài)保護的原因,本研究未能大量采集樣木進行調(diào)查;因此樣木選擇時盡量依據(jù)樣木選擇標準及不同演替階段樹木徑級變化,每個樹種基本能反應其不同生長階段的變化。每個樹種樣本的基本情況見表1。

    1.3 樣木生物量測定

    選取樣木后,首先測量其基徑、胸徑、南北向冠幅及東西向冠幅。隨后貼近地面將樹伐倒,測量樹高,然后去除枝葉,分別稱量樹干、樹根、枝和葉的重量(盡可能的收集樹木伐倒時掉落的枝葉)。將整個樹冠劃分為上、中、下3 層,用臺秤現(xiàn)場稱量各層鮮重;從每層各選取3枝不同大小的枝條,現(xiàn)場用游標卡尺及卷尺測量每個枝條基部直徑及長度,測定鮮重后帶回實驗室;同時在樹干基徑、胸徑及1/4、1/2和3/4高處分別截取圓盤,稱重后帶回實驗室。

    表1 樣木基本信息統(tǒng)計

    所有樣木地下部分生物量采取全挖法獲得。個別樹根難以全部挖出時根據(jù)根冠比來進行估算[27-28]。挖出后清除表面泥土并分為根樁、粗根(>3 cm)、中根、細根(<1 cm)4個級別,稱量鮮重后各取部分(500~1000 g)帶回實驗室,用流動水浸泡、漂洗。所有帶回樣品經(jīng)60 ℃烘干、稱重,計算出各器官含水率,根據(jù)各部分干濕比推算各器官生物量。總生物量是所有器官生物量的總和。

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

    大量文獻表明,冪函數(shù)y=axb是描述林木及其器官相對生長關系最常用的模型[29-31]。由于建立生物量模型的目的是為了以后進行更好地預估,因此在模型中一般盡可能選擇和生物量關系密切并在林木中容易獲取的測樹因子[10-12]。胸徑(D)由于是最容易準確獲取的指標,而且和生物量關系密切,因此常常被作為估算立木生物量的最重要指標在世界范圍內(nèi)得到廣泛應用[31]。目前,基于胸徑和樹高的模型主要有以下3種:①基于胸徑(D)的一元模型;②以D2H來代替D作為自變量的二元生物量模型;③以D3/H為自變量的二元生物量模型[9,32]。因此,本研究中也選用上述3個異速生長方程進行生物量模型研究,3個方程相應的表達式如下:

    W=aDb

    (1)

    W=a(D2H)b

    (2)

    W=a(D3/H)b

    (3)

    式中,W為樹木各部分(器官)生物量及總生物量,D是胸徑,H為樹高,a、b是模型參數(shù)。采用相關系數(shù)R2、估計值標準誤差SEE、總相對誤差TRE、平均系統(tǒng)誤差MSE作為模型精度的評價指標[33-34]。使用SPSS和Origin軟件進行數(shù)據(jù)分析和圖形制作。

    表2 樹木各組分生物量異速生長方程與評價指標

    續(xù)表2 Continued table 2

    物種Species組分Component變量Variable參數(shù)Parameter abR2P指標IndexSEETREMSED2H0.0240.9250.993<0.0019.300-1.6890.655D3/H0.3150.9390.978<0.00127.100-6.2592.168根RootD0.0082.4280.922<0.013.9307.1267.010D2H0.0020.9510.906<0.014.19611.8438.186D3/H0.0310.9880.935<0.014.5152.6316.073總TotalD0.0932.3460.992<0.00120.572-3.9710.732D2H0.0260.9290.996<0.0019.891-1.4040.365D3/H0.3480.9450.985<0.00131.316-6.2791.476川滇柳S. rehderiana枝BranchD0.0361.7920.983<0.0010.2871.0900.627D2H0.0160.7010.959<0.0010.3082.6681.457D3/H0.0940.7160.985<0.0010.432-0.0020.577葉FoliageD0.0141.8980.979<0.0010.216-0.4820.907D2H0.0050.7570.994<0.0010.1460.3190.245D3/H0.0410.7430.942<0.0010.307-0.6742.548樹干StemD0.1162.0690.924<0.017.1444.6024.070D2H0.0400.8220.929<0.017.5785.9283.858D3/H0.3690.8140.898<0.017.0383.9705.624地上AbovegroundD0.1652.0160.945<0.017.3973.1942.754D2H0.0590.7990.946<0.017.8684.5212.716D3/H0.5050.7950.921<0.017.4022.5184.019根RootD0.0521.7530.928<0.011.2084.4132.623D2H0.0230.6870.907<0.011.3865.9253.523D3/H0.1340.7000.928<0.011.0863.3932.631總TotalD0.2151.9730.945<0.018.4603.2752.591D2H0.0800.7810.944<0.019.1324.6312.725D3/H0.6370.7790.925<0.018.3382.5433.628沙棘H. rhamnoides枝BranchD0.0421.8010.783<0.050.4972.3882.107D2H0.0300.6450.758<0.050.5362.6132.360D3/H0.0640.8130.811<0.050.4502.1351.825葉FoliageD0.0042.6450.866<0.050.3441.1632.804D2H0.0030.9500.843<0.050.3751.5353.278D3/H0.0081.1890.891<0.050.3110.7462.284樹干StemD0.0782.3790.945<0.011.598-0.0830.811D2H0.0330.9040.999<0.011.556-0.1190.791D3/H0.1501.0510.940<0.011.750-0.0070.884地上Above groundD0.1162.3000.971<0.011.385-0.2560.386D2H0.0590.8540.978<0.011.589-0.2120.451D3/H0.2151.0210.973<0.011.311-0.2790.354根RootD0.0871.5970.820<0.050.6251.8161.290D2H0.0750.5420.692<0.050.6742.0071.496D3/H0.1270.7210.851<0.050.5641.5981.060總TotalD0.1712.1920.975<0.011.576-0.1140.294D2H0.0940.8060.965<0.011.945-0.0460.383D3/H0.3060.9750.981<0.011.251-0.1680.229

    注:a、b代表3個冪函數(shù)方程的參數(shù),R2指相關系數(shù)、SEE指估計值標準誤差、TRE指總相對誤差、MSE指平均系統(tǒng)誤差。

    圖1 樣木不同器官生物量相對分配比例Fig.1 Relative biomass allocation of the different components

    2 結果與分析

    2.1 生物量分配格局

    本研究涉及的4個樹種地上總生物量占比差別不大(圖1),均在80 %左右;從整株水平來看,峨眉冷杉、川滇柳、沙棘3個樹種不同器官生物量分配比例均表現(xiàn)為:干>根>枝>葉(P< 0.001),冬瓜楊各器官生物量分配比例表現(xiàn)為干>枝>根>葉(P< 0.01)。4個樹種的樹干生物量均占整株生物量的60 %左右;但相對而言沙棘樹干占整株生物量的比例最高,達67.9 %,枝生物量的占比最低,僅為10.9 %;而冬瓜楊枝生物量所占比例在4個樹種中最高,達16.6 %,根生物量占比最低,僅占14.7 %。4個樹種地上生物量與地下生物量的平均比值為4.72,其中冬瓜楊比值最高,為5.79,沙棘地上、地下生物量比值最低,僅3.92。

    2.2 生物量異速生長模型

    R2、SEE、TRE和MSE4項統(tǒng)計指標中,R2數(shù)值越大效果越好,其它3項指標均是數(shù)值越小(其中TRE按絕對值進行比較)效果越好[7,32]。

    分別采用(1)~(3)3個模型對峨眉冷杉、冬瓜楊、川滇柳和沙棘4個樹種不同器官生物量及總生物量進行擬合,結果如表2所示,表明3個模型總體回歸F檢驗均顯示出顯著性(P<0.05),因此均適合四種樹木各器官生物量和總生物量的估算。對于峨眉冷杉,模型(2)對樹干及地上和地下生物量的解釋量分別為97 %、94 %和94 %,對枝和葉生物量的解釋量稍低,約為76 %左右;模型(2)對峨眉冷杉各器官及總生物量的擬合效果明顯優(yōu)于模型(1)和(3)。冬瓜楊3個方程的R2值都大于0.9,樹干、地上和總生物量模型(2)擬合效果最好,枝葉和地下生物量模型(3)擬合效果最好。川滇柳3個方程的擬合效果相差不大,總的來說模型(2)對樹干和地上生物量的擬合效果最好,模型(3)對枝生物量和地下生物量擬合效果較好,模型(1)對葉生物量擬合效果較好。沙棘枝、葉、地下生物量和總生物量,模型(3)的擬合效果最優(yōu),樹干和地上生物量模型(2)擬合效果最優(yōu);但如果考慮單一方程的擬合效果,枝葉和樹根模型的解釋能力(R2)明顯低于其他器官。

    對比模型擬合效果評價指標SEE、TRE和MSE發(fā)現(xiàn),樹干和地上生物量模型在引入樹高H,并以D2H作為自變量后,SEE、TRE和MSE都有一定程度的降低,模型(2)擬合效果較好。枝、葉和根生物量各評價指標顯示,除冷杉外另外3個物種模型(3)的擬合效果最好,模型(1)次之,模型(2)最差,說明以D2H為自變量的方程不適合枝、葉和根生物量的擬合。

    3 討 論

    3.1 生物量分配格局

    生物量是植物積累能量的主要體現(xiàn),其分配方式受外界環(huán)境(光照、水分、營養(yǎng)等)以及植株年齡大小的限制[35-38]。在植被演替過程中,優(yōu)勢木的生長與森林的更新密切相關[39],針葉和闊葉受光照的影響不同,其生物量的分配也存在差異。有人得出[36]:通常地下生物量占21 %,樹干生物量占全株生物量的比例范圍在25 %~76 %,本研究得出的樹干生物量占60 %左右,但是,地下生物量占比在14 %~21 %,說明在冰川退縮區(qū)原生演替過程中,根在總生物量中的占比可能跟演替進展有關,激烈的種間和種內(nèi)斗爭,導致根系分布變化不均,這個問題有待繼續(xù)研究。

    3.2 模型評價

    生長方程的回歸分析中最重要的一步就是模型評價,經(jīng)報道過的用于模型評價的指標很多,如SEE、TRE、AEE、RE、MSE、MPE、R2、FI等,但一些指標(如FI)僅被用于特殊情況[7,40]。本研究中,構建一元和二元生物量模型時,選取最常用的R2、SEE、TRE和MSE作為模型評價的指標。前人在模型評價時使用交叉驗證[41],即將可用的樣本數(shù)據(jù)分成2組,一組用于模型的建立,一組用于模型驗證,但已有研究稱交叉驗證存在模型信息缺失的缺點[7],鑒于本次研究的樣本有限,不使用交叉驗證。雖然如此,本研究初步反映了川西亞高山生態(tài)系統(tǒng)植被原生演替過程中優(yōu)勢樹種單株及各器官異速生長過程,相關結果對于該區(qū)域植被原生演替過程及其生態(tài)功能研究具有促進作用。

    3.3 最優(yōu)模型與影響因素

    用于生物量預測的模型很多,最常用的不外乎以D為自變量的一元模型以及以D、H相結合的二元模型。大量研究證明胸徑D是生物量預測的最好因素[8,30-31],這也在本研究中得到了驗證。除個別枝葉的擬合外,以D為自變量的方程能解釋超過90 %的生物量變異。雖然以D2H作為自變量的方程能提高方程的擬合優(yōu)度,但由于H的測量相對麻煩,所以現(xiàn)實條件下較少使用。此外,曾偉生[10]和張繪芳[9]曾提出以D3/H為自變量的方程能較好擬合枝葉和根的生物量,這個觀點在本研究中也得到了證實。

    4 結 論

    (1)不同植株器官生物量分配存在一定的差異,本研究得出峨眉冷杉、川滇柳、沙棘3個樹種不同器官生物量分配比例均表現(xiàn)為:干>根>枝>葉,冬瓜楊各器官生物量分配比例表現(xiàn)為干>枝>根>葉。地下生物量約占14 %~21 %,地上生物量約占79 %~86 %,4個樹種得到相似的結論。

    (2)同一模型對地上和樹干生物量估計精度要優(yōu)于枝、葉和地下生物量的估計。究其原因可能是胸徑D和樹高H都與樹干生長密切相關,枝葉通常受光照和林分密度的影響較大,在原生演替過程中種間和種內(nèi)競爭劇烈,導致樹冠冠幅和根系分布發(fā)生相應變化。

    (3)引入樹高H的二元模型擬合效果要優(yōu)于一元模型,以D2H為自變量的方程對樹干和地上生物量的擬合效果更好,而以D3/H為自變量的方程更適合枝葉和根的生物量擬合。由于立木的生長受到多種因素的影響,增加變量能更好的擬合其生物量的變化,特別是不同器官的生物量最好選擇合適的變量來擬合。

    (4)盡管有很多異速生長模型可用于生物量的預測,但應該權衡精度、簡潔和實用性,考慮到野外測量的工作量較大,且樹高測量存在一定范圍的誤差,所以,在實際工作中,對精度要求不是很高的情況下,建議采用一元模型。

    猜你喜歡
    樣木沙棘樹干
    沙棘種植讓日子甜滋滋
    沙棘頌
    通化縣林業(yè)局天然蒙古櫟林生長變化情況研究
    為什么樹干不是方的?
    軍事文摘(2021年16期)2021-11-05 08:49:16
    沙棘在西藏的發(fā)展前景探討
    沙棘產(chǎn)業(yè)的直銷之路
    鎮(zhèn)賚縣人工楊樹5年間蓄積量變化情況分析
    鎮(zhèn)賚縣人工楊樹5 年間胸徑蓄積量變化情況探討
    為什么要在樹干上刷白漿
    為什么要在樹干上刷白漿
    日韩一区二区视频免费看| 国产 一区精品| xxx大片免费视频| 久久人人爽av亚洲精品天堂| 丝袜脚勾引网站| 欧美 日韩 精品 国产| 1024香蕉在线观看| 母亲3免费完整高清在线观看| 一本—道久久a久久精品蜜桃钙片| 啦啦啦中文免费视频观看日本| 国产片特级美女逼逼视频| 人人妻,人人澡人人爽秒播 | 80岁老熟妇乱子伦牲交| 母亲3免费完整高清在线观看| 中文字幕色久视频| 婷婷色综合www| 精品国产超薄肉色丝袜足j| 美女中出高潮动态图| av.在线天堂| 99久久精品国产亚洲精品| av国产精品久久久久影院| 国产精品国产三级国产专区5o| 精品一品国产午夜福利视频| 一级片免费观看大全| 大码成人一级视频| 亚洲国产精品成人久久小说| 国产精品久久久人人做人人爽| 成人毛片60女人毛片免费| 成人国产av品久久久| 黄片小视频在线播放| 国产免费福利视频在线观看| 18在线观看网站| 欧美国产精品一级二级三级| 在线免费观看不下载黄p国产| 丁香六月天网| 美女脱内裤让男人舔精品视频| 午夜老司机福利片| 日韩免费高清中文字幕av| av线在线观看网站| 免费观看av网站的网址| 高清av免费在线| 欧美人与性动交α欧美软件| av线在线观看网站| 欧美xxⅹ黑人| 一级毛片电影观看| 国产精品免费视频内射| 久久av网站| 国产极品天堂在线| 国产成人av激情在线播放| 青春草国产在线视频| 赤兔流量卡办理| 日韩av在线免费看完整版不卡| 夫妻午夜视频| 久久精品久久久久久久性| 欧美老熟妇乱子伦牲交| 69精品国产乱码久久久| 国产免费一区二区三区四区乱码| 久久国产精品男人的天堂亚洲| 在线天堂最新版资源| 日本vs欧美在线观看视频| 丝袜美足系列| 国产野战对白在线观看| 最近最新中文字幕大全免费视频 | 国产精品久久久久成人av| 精品亚洲成a人片在线观看| 亚洲人成77777在线视频| 亚洲欧洲国产日韩| 一二三四中文在线观看免费高清| 新久久久久国产一级毛片| 久久国产精品大桥未久av| av片东京热男人的天堂| 男人爽女人下面视频在线观看| 亚洲国产看品久久| 超色免费av| 9色porny在线观看| 伊人久久大香线蕉亚洲五| 久久精品人人爽人人爽视色| 视频区图区小说| 毛片一级片免费看久久久久| 国产又爽黄色视频| 亚洲七黄色美女视频| 亚洲综合色网址| a级毛片黄视频| 日韩一区二区视频免费看| 1024视频免费在线观看| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 女性生殖器流出的白浆| 自拍欧美九色日韩亚洲蝌蚪91| 在现免费观看毛片| 久久毛片免费看一区二区三区| 亚洲激情五月婷婷啪啪| 久久久久人妻精品一区果冻| 熟妇人妻不卡中文字幕| kizo精华| 欧美亚洲 丝袜 人妻 在线| 国产男人的电影天堂91| 少妇猛男粗大的猛烈进出视频| 天天躁夜夜躁狠狠躁躁| 久久亚洲国产成人精品v| 波野结衣二区三区在线| 亚洲精品,欧美精品| 啦啦啦在线观看免费高清www| 久久久久网色| 亚洲第一青青草原| 成人国产麻豆网| 天天躁日日躁夜夜躁夜夜| 一个人免费看片子| 汤姆久久久久久久影院中文字幕| 午夜福利视频精品| 日韩电影二区| 成人手机av| 免费观看性生交大片5| 日韩制服丝袜自拍偷拍| 成人毛片60女人毛片免费| 久久久精品94久久精品| 狠狠精品人妻久久久久久综合| 亚洲成色77777| 十分钟在线观看高清视频www| 女性被躁到高潮视频| 日韩一区二区视频免费看| 又大又爽又粗| 色播在线永久视频| 免费少妇av软件| 久久久久久久久久久久大奶| 天美传媒精品一区二区| 女人高潮潮喷娇喘18禁视频| 天天躁夜夜躁狠狠久久av| 伦理电影免费视频| 伊人久久大香线蕉亚洲五| 叶爱在线成人免费视频播放| 男女国产视频网站| 欧美成人精品欧美一级黄| 亚洲av电影在线进入| 日韩中文字幕视频在线看片| 最近的中文字幕免费完整| 国产伦理片在线播放av一区| e午夜精品久久久久久久| 亚洲一区中文字幕在线| 日韩,欧美,国产一区二区三区| 久久久久精品久久久久真实原创| 丝袜脚勾引网站| 热99久久久久精品小说推荐| 久久这里只有精品19| 免费在线观看视频国产中文字幕亚洲 | 日韩熟女老妇一区二区性免费视频| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久成人av| 两个人免费观看高清视频| 如日韩欧美国产精品一区二区三区| 久久 成人 亚洲| 91成人精品电影| 亚洲精品aⅴ在线观看| 国产日韩一区二区三区精品不卡| 国产亚洲精品第一综合不卡| 亚洲人成电影观看| 国产1区2区3区精品| 亚洲成人免费av在线播放| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 免费黄网站久久成人精品| 成年美女黄网站色视频大全免费| 日韩成人av中文字幕在线观看| 久久久久久人妻| 欧美日韩一区二区视频在线观看视频在线| 国产成人av激情在线播放| 十八禁高潮呻吟视频| 亚洲av欧美aⅴ国产| 人人妻人人添人人爽欧美一区卜| 国产日韩欧美在线精品| 日韩大码丰满熟妇| 中文字幕精品免费在线观看视频| 欧美变态另类bdsm刘玥| 少妇人妻久久综合中文| 久热爱精品视频在线9| 老汉色av国产亚洲站长工具| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲三区欧美一区| av.在线天堂| 搡老乐熟女国产| 狠狠精品人妻久久久久久综合| 久久国产精品大桥未久av| 一级片'在线观看视频| 99热全是精品| 男女之事视频高清在线观看 | 久久韩国三级中文字幕| 考比视频在线观看| av片东京热男人的天堂| 亚洲精品国产区一区二| 免费观看av网站的网址| 啦啦啦在线免费观看视频4| 亚洲欧美色中文字幕在线| 免费高清在线观看视频在线观看| 日韩免费高清中文字幕av| 亚洲视频免费观看视频| 最新在线观看一区二区三区 | 午夜精品国产一区二区电影| 亚洲精品视频女| 观看美女的网站| 精品亚洲成a人片在线观看| 亚洲精品久久久久久婷婷小说| 精品卡一卡二卡四卡免费| 亚洲国产看品久久| 亚洲av成人精品一二三区| 国产精品久久久久成人av| 男女床上黄色一级片免费看| 国产成人91sexporn| 亚洲欧美精品综合一区二区三区| 人体艺术视频欧美日本| 亚洲精品久久午夜乱码| 国产午夜精品一二区理论片| 成人国语在线视频| 在线亚洲精品国产二区图片欧美| 这个男人来自地球电影免费观看 | 日韩人妻精品一区2区三区| 一二三四中文在线观看免费高清| 尾随美女入室| 中文字幕色久视频| 黑人巨大精品欧美一区二区蜜桃| 大香蕉久久成人网| 欧美国产精品va在线观看不卡| 亚洲男人天堂网一区| 中国三级夫妇交换| 免费看av在线观看网站| 久久久亚洲精品成人影院| 卡戴珊不雅视频在线播放| 久久精品熟女亚洲av麻豆精品| av免费观看日本| av女优亚洲男人天堂| 成年人免费黄色播放视频| 亚洲综合色网址| 999久久久国产精品视频| 男女下面插进去视频免费观看| 久久青草综合色| 一区福利在线观看| 桃花免费在线播放| 亚洲美女搞黄在线观看| 亚洲av国产av综合av卡| 曰老女人黄片| 欧美日韩亚洲综合一区二区三区_| 久久午夜综合久久蜜桃| 秋霞伦理黄片| 操美女的视频在线观看| 在线免费观看不下载黄p国产| 午夜福利,免费看| 哪个播放器可以免费观看大片| 色播在线永久视频| 国产在线视频一区二区| 2018国产大陆天天弄谢| 丁香六月天网| 精品久久久久久电影网| 另类精品久久| 只有这里有精品99| 午夜激情av网站| 国产高清不卡午夜福利| 在线观看免费高清a一片| 国产又色又爽无遮挡免| 国产精品人妻久久久影院| 国产精品蜜桃在线观看| 日韩一本色道免费dvd| 人妻人人澡人人爽人人| 欧美最新免费一区二区三区| 亚洲 欧美一区二区三区| av视频免费观看在线观看| 亚洲免费av在线视频| 超碰97精品在线观看| 天天操日日干夜夜撸| 男男h啪啪无遮挡| 蜜桃在线观看..| 爱豆传媒免费全集在线观看| 成年av动漫网址| 一区二区三区精品91| 日本欧美国产在线视频| 一级毛片黄色毛片免费观看视频| 亚洲精品av麻豆狂野| 少妇被粗大的猛进出69影院| 国产男人的电影天堂91| 一级爰片在线观看| 又大又黄又爽视频免费| 黄网站色视频无遮挡免费观看| 成人毛片60女人毛片免费| av国产久精品久网站免费入址| 国产毛片在线视频| 搡老乐熟女国产| 亚洲国产看品久久| 777米奇影视久久| 日韩制服丝袜自拍偷拍| 久久精品久久久久久噜噜老黄| 99热全是精品| 亚洲av中文av极速乱| 日韩免费高清中文字幕av| 五月开心婷婷网| 精品第一国产精品| 丝袜脚勾引网站| 丝袜美腿诱惑在线| 王馨瑶露胸无遮挡在线观看| 波野结衣二区三区在线| 国产av国产精品国产| 国产不卡av网站在线观看| 国产成人91sexporn| 久久久久久久大尺度免费视频| 超色免费av| 下体分泌物呈黄色| 欧美精品一区二区大全| 老司机亚洲免费影院| 夫妻性生交免费视频一级片| 看免费av毛片| 老汉色av国产亚洲站长工具| 免费观看性生交大片5| 久久亚洲国产成人精品v| 欧美中文综合在线视频| 丁香六月天网| 成人国语在线视频| 97在线人人人人妻| 亚洲av男天堂| 久久久国产欧美日韩av| 操美女的视频在线观看| 男女午夜视频在线观看| 最近最新中文字幕免费大全7| 久久久久久人妻| 亚洲天堂av无毛| 黄网站色视频无遮挡免费观看| 国产精品 国内视频| 欧美 亚洲 国产 日韩一| 亚洲成色77777| 日日摸夜夜添夜夜爱| 国产毛片在线视频| 欧美精品av麻豆av| 午夜福利一区二区在线看| 色综合欧美亚洲国产小说| 少妇精品久久久久久久| 爱豆传媒免费全集在线观看| 亚洲国产毛片av蜜桃av| 国产成人欧美在线观看 | 国产熟女欧美一区二区| 女人爽到高潮嗷嗷叫在线视频| 色94色欧美一区二区| 啦啦啦中文免费视频观看日本| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 欧美日本中文国产一区发布| 精品少妇久久久久久888优播| 免费在线观看完整版高清| 欧美精品av麻豆av| 久久亚洲国产成人精品v| 伦理电影免费视频| 久久av网站| 中国国产av一级| 欧美变态另类bdsm刘玥| 久久影院123| 欧美变态另类bdsm刘玥| 成人18禁高潮啪啪吃奶动态图| 精品人妻在线不人妻| 一边亲一边摸免费视频| 欧美变态另类bdsm刘玥| 青春草亚洲视频在线观看| 国产av一区二区精品久久| 美女扒开内裤让男人捅视频| 欧美变态另类bdsm刘玥| 亚洲精品国产色婷婷电影| 精品人妻在线不人妻| 男女边吃奶边做爰视频| 国产一卡二卡三卡精品 | 国产免费又黄又爽又色| av网站在线播放免费| 国产色婷婷99| 国产激情久久老熟女| 精品一品国产午夜福利视频| 伦理电影免费视频| 人人妻人人爽人人添夜夜欢视频| xxxhd国产人妻xxx| 亚洲av日韩精品久久久久久密 | 国产片特级美女逼逼视频| 最近最新中文字幕大全免费视频 | 在线观看免费视频网站a站| 亚洲精品国产一区二区精华液| 久久国产精品男人的天堂亚洲| 精品午夜福利在线看| 波多野结衣av一区二区av| 成人国语在线视频| 男女无遮挡免费网站观看| 高清av免费在线| 一区福利在线观看| 国产高清国产精品国产三级| 一本一本久久a久久精品综合妖精| 久久热在线av| 国产欧美亚洲国产| 毛片一级片免费看久久久久| 欧美激情高清一区二区三区 | 亚洲av国产av综合av卡| 另类精品久久| 日韩不卡一区二区三区视频在线| 91成人精品电影| 成年女人毛片免费观看观看9 | 午夜免费男女啪啪视频观看| 黑人欧美特级aaaaaa片| 1024香蕉在线观看| 黄片播放在线免费| 国产欧美亚洲国产| 咕卡用的链子| 国产免费福利视频在线观看| 亚洲欧美中文字幕日韩二区| 日韩一本色道免费dvd| 少妇 在线观看| 一区二区三区精品91| 亚洲精品国产区一区二| 人人妻人人澡人人看| 嫩草影视91久久| 色婷婷av一区二区三区视频| 国产一区有黄有色的免费视频| av网站在线播放免费| 久热这里只有精品99| 黑人猛操日本美女一级片| 亚洲视频免费观看视频| 啦啦啦中文免费视频观看日本| 女性生殖器流出的白浆| 欧美人与性动交α欧美软件| 纵有疾风起免费观看全集完整版| 亚洲精品第二区| 免费高清在线观看视频在线观看| 女人被躁到高潮嗷嗷叫费观| 国产成人免费观看mmmm| 久久婷婷青草| 黄片播放在线免费| 90打野战视频偷拍视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲欧美色中文字幕在线| 日韩av免费高清视频| 国产不卡av网站在线观看| 观看美女的网站| av有码第一页| 超碰成人久久| 国产黄色免费在线视频| 中文字幕制服av| 日韩一卡2卡3卡4卡2021年| 一级爰片在线观看| 欧美精品av麻豆av| 国产在线免费精品| 免费在线观看视频国产中文字幕亚洲 | 青春草视频在线免费观看| 欧美少妇被猛烈插入视频| 亚洲,欧美精品.| 国产av国产精品国产| 色网站视频免费| 久久久久网色| 成人手机av| 亚洲美女搞黄在线观看| 国产男人的电影天堂91| 尾随美女入室| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人| 免费黄色在线免费观看| 欧美日韩视频精品一区| 搡老乐熟女国产| 亚洲成人国产一区在线观看 | 欧美成人精品欧美一级黄| 亚洲伊人久久精品综合| 欧美97在线视频| 国产高清国产精品国产三级| 女性生殖器流出的白浆| 色婷婷久久久亚洲欧美| 成人手机av| 久久精品国产a三级三级三级| xxx大片免费视频| kizo精华| 丝袜美足系列| 亚洲av欧美aⅴ国产| 久久精品久久久久久久性| 少妇人妻 视频| 一个人免费看片子| 制服人妻中文乱码| 90打野战视频偷拍视频| 丁香六月欧美| 大话2 男鬼变身卡| 天堂俺去俺来也www色官网| 日本一区二区免费在线视频| 国产日韩欧美在线精品| 9色porny在线观看| 色网站视频免费| 黄色视频在线播放观看不卡| 日本色播在线视频| 精品亚洲成a人片在线观看| 亚洲成国产人片在线观看| 午夜免费男女啪啪视频观看| 午夜福利视频精品| 国产精品久久久久久人妻精品电影 | 亚洲成人一二三区av| 午夜老司机福利片| 久久婷婷青草| 亚洲,一卡二卡三卡| 国产av码专区亚洲av| 777久久人妻少妇嫩草av网站| 日本av免费视频播放| 亚洲专区中文字幕在线 | 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 亚洲国产欧美一区二区综合| 蜜桃国产av成人99| 国产精品亚洲av一区麻豆 | 亚洲,一卡二卡三卡| 丝袜美腿诱惑在线| 狂野欧美激情性bbbbbb| 99热国产这里只有精品6| 一级爰片在线观看| 观看美女的网站| 黑丝袜美女国产一区| a 毛片基地| 99香蕉大伊视频| 国产乱人偷精品视频| 国产亚洲欧美精品永久| 捣出白浆h1v1| netflix在线观看网站| 欧美在线一区亚洲| 久久久久久人人人人人| 久久久国产一区二区| 99re6热这里在线精品视频| 亚洲精品国产av成人精品| 麻豆av在线久日| 自线自在国产av| 夫妻午夜视频| 午夜福利免费观看在线| 一级毛片我不卡| 国产成人精品无人区| 婷婷色av中文字幕| xxx大片免费视频| 国产成人一区二区在线| 日韩免费高清中文字幕av| 只有这里有精品99| 亚洲成人免费av在线播放| 国产人伦9x9x在线观看| 看十八女毛片水多多多| 一二三四在线观看免费中文在| 日韩不卡一区二区三区视频在线| 80岁老熟妇乱子伦牲交| 亚洲国产欧美网| 麻豆av在线久日| 亚洲精品一区蜜桃| 中国国产av一级| 在线天堂最新版资源| 久久久精品区二区三区| 可以免费在线观看a视频的电影网站 | 久久久久国产一级毛片高清牌| 尾随美女入室| 欧美人与性动交α欧美软件| 制服诱惑二区| 国产日韩欧美在线精品| 亚洲,一卡二卡三卡| 中文天堂在线官网| 亚洲美女搞黄在线观看| 欧美精品av麻豆av| 国产熟女欧美一区二区| 超碰97精品在线观看| 日日撸夜夜添| 国产免费又黄又爽又色| 男女国产视频网站| 午夜福利免费观看在线| 国产日韩欧美视频二区| 18禁裸乳无遮挡动漫免费视频| 国产精品二区激情视频| av不卡在线播放| 婷婷色综合大香蕉| 午夜福利乱码中文字幕| 国产免费视频播放在线视频| 国产精品 国内视频| 亚洲精品一区蜜桃| 欧美久久黑人一区二区| 一级毛片黄色毛片免费观看视频| 蜜桃国产av成人99| 免费久久久久久久精品成人欧美视频| 久久久国产精品麻豆| 国产亚洲av高清不卡| 最新在线观看一区二区三区 | 肉色欧美久久久久久久蜜桃| 蜜桃国产av成人99| 日本猛色少妇xxxxx猛交久久| 90打野战视频偷拍视频| 日韩一区二区三区影片| 可以免费在线观看a视频的电影网站 | 欧美另类一区| 欧美日韩综合久久久久久| 日韩中文字幕欧美一区二区 | 色婷婷久久久亚洲欧美| bbb黄色大片| 亚洲欧美成人精品一区二区| 99久国产av精品国产电影| 最黄视频免费看| 国产午夜精品一二区理论片| 婷婷色综合www| 中文字幕人妻熟女乱码| 一区二区日韩欧美中文字幕| 一级毛片我不卡| 亚洲美女黄色视频免费看| 晚上一个人看的免费电影| e午夜精品久久久久久久| 午夜福利乱码中文字幕| 国产99久久九九免费精品| 老鸭窝网址在线观看| 国产深夜福利视频在线观看| 搡老乐熟女国产| 亚洲欧美精品自产自拍| 成年av动漫网址| 成人漫画全彩无遮挡| 国产精品av久久久久免费| 中文乱码字字幕精品一区二区三区| 亚洲人成电影观看| 飞空精品影院首页| 国产精品.久久久| 麻豆乱淫一区二区| 桃花免费在线播放| 在线天堂中文资源库| 免费黄色在线免费观看| 中文字幕色久视频| 亚洲av日韩精品久久久久久密 | 女人高潮潮喷娇喘18禁视频| 精品少妇久久久久久888优播|