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

    古爾班通古特沙漠固沙灌木的冠幅預(yù)測(cè)模型

    2022-08-04 09:10:08路麗寧張定海李功麟
    生態(tài)科學(xué) 2022年5期
    關(guān)鍵詞:固沙梭梭冠幅

    路麗寧, 張定海, 李功麟

    古爾班通古特沙漠固沙灌木的冠幅預(yù)測(cè)模型

    路麗寧, 張定海*, 李功麟

    甘肅農(nóng)業(yè)大學(xué)理學(xué)院數(shù)量生物研究所, 蘭州 730070

    以古爾班通古特沙漠中三種類型沙丘(固定沙丘、半固定沙丘和流動(dòng)沙丘)上主要固沙灌木梭梭()、白梭梭()、沙拐棗()、蛇麻黃()和油蒿()為研究對(duì)象, 利用12個(gè)基礎(chǔ)模型建立了不同沙丘類型上所有灌木和不同類型灌木的冠幅模型。模型擬合結(jié)果表明: 不論是不同沙丘類型上所有灌木的擬合效果還是不同類型灌木的擬合效果, 固定沙丘和半固定沙丘均優(yōu)于流動(dòng)沙丘。同時(shí), 兩個(gè)近“S”形模型M7(Gompertz Model)和M12(HossfeldⅠ Model)對(duì)三種類型沙丘上的固沙灌木均有較好擬合效果; M1(Linear Model)、M2(Quadratic Model)、M3(Power Model)、M6(Logistic Model)、M7、M9(Gauss Model)、M10(Mitscherlinch Model)和M12八個(gè)模型適用于研究區(qū)絕大部分固沙灌木冠幅模型的模擬。但不同沙丘類型上所有灌木和不同類型的灌木都有不同的最優(yōu)冠幅模型。研究可為今后該地區(qū)相關(guān)生態(tài)模型建立和防風(fēng)固沙效益的評(píng)價(jià)提供科學(xué)依據(jù)。

    古爾班通古特; 固沙灌木; 沙丘類型; 冠幅預(yù)測(cè)模型

    0 前言

    植物冠層是植物進(jìn)行光合作用和有機(jī)物積累的重要場(chǎng)所[1-3], 植物冠層的大小與植物的生物量和生長(zhǎng)狀況密切相關(guān), 反映了植物的生存和競(jìng)爭(zhēng)能力[4], 是描述植物活力與健康的重要指標(biāo)[5-7]。衡量植冠大小的冠幅是預(yù)測(cè)植冠表面積、生物量和形狀的重要變量[8], 冠幅與葉表面積和體積顯著相關(guān)[9], 植冠的表面積和體積又與植物光合的規(guī)模相關(guān)。冠幅尺寸的測(cè)量和分析對(duì)于量化和鑒定植物生長(zhǎng)期、穩(wěn)定性和植物的生產(chǎn)效率具有重要意義[10]。目前, 針對(duì)森林生態(tài)系統(tǒng)中樹(shù)冠長(zhǎng)度、樹(shù)冠率和冠幅, 已經(jīng)開(kāi)展了大量的模擬研究[10-13]。冠幅作為模型預(yù)測(cè)的重要變量在樹(shù)木個(gè)體的死亡模型、地上生物量模型、立木材積和削度方程等模型研究中發(fā)揮了重要的作用[14-16]。

    對(duì)于冠幅模型的研究, 通常的做法是通過(guò)多元回歸的方法建立冠幅和林木變量及林分變量的各種形式的預(yù)測(cè)模型[17-21]。Gill等人建立了區(qū)域性的主要樹(shù)種冠幅的線性預(yù)測(cè)模型, 同時(shí)指出樹(shù)木的胸徑和樹(shù)高是影響冠幅的顯著因子[17]。周志翔等人選擇線性函數(shù)和幾種可以轉(zhuǎn)換為線性函數(shù)的非線性函數(shù)研究了不同密度濕地松幼齡林冠幅和栽植年齡之間的關(guān)系[8, 17]; 韓艷剛等人收集了702 棵樟子松立木數(shù)據(jù), 構(gòu)建了胸徑與相關(guān)林木因子和林分變量的冠幅—胸徑模型[22]。目前, 常見(jiàn)的冠幅—胸徑模型的形式有十幾種, 世界各地研究人員利用這些模型開(kāi)展了大量的有關(guān)冠幅—胸徑模型的研究[21, 23-26]。

    當(dāng)今世界面臨最大的環(huán)境、經(jīng)濟(jì)問(wèn)題之一就是土地沙漠化。土地沙漠化面積的迅速擴(kuò)展, 造成環(huán)境退化和巨大的經(jīng)濟(jì)損失, 使之成為全球廣泛關(guān)注的熱點(diǎn)[27]。古爾班通古特沙漠是中國(guó)土地沙漠化較嚴(yán)重的地區(qū)之一, 固沙灌木對(duì)干旱半干旱地區(qū)生態(tài)系統(tǒng)的恢復(fù)有著十分重要的作用。冠層是固沙灌木進(jìn)行光合作用的重要場(chǎng)所, 它決定灌木的活力和生產(chǎn)力, 也是灌木在長(zhǎng)期的水分脅迫下水分競(jìng)爭(zhēng)的重要指標(biāo)。冠幅是冠層結(jié)構(gòu)的重要特征因子之一, 它對(duì)于建立沙區(qū)灌木的生長(zhǎng)模型、計(jì)算灌木的競(jìng)爭(zhēng)指數(shù)、沙區(qū)植被的生物量反演以及評(píng)價(jià)固沙效益的優(yōu)劣等方面具有重要的意義。對(duì)固沙灌木冠幅預(yù)測(cè)模型的探究有利于分析固沙植被的活力、穩(wěn)定性和固沙植被生物量模型等許多方面的研究, 而這些研究都對(duì)干旱半干旱生態(tài)系統(tǒng)防風(fēng)固沙有十分重要的作用。但是, 對(duì)干旱半干旱生態(tài)系統(tǒng)健康發(fā)揮重要作用的固沙灌木冠幅預(yù)測(cè)模型的相關(guān)研究幾乎未見(jiàn)報(bào)道。這主要是由于目前專門針對(duì)干旱生態(tài)系統(tǒng)而開(kāi)展的相關(guān)模型研究比較缺乏。因此, 開(kāi)發(fā)準(zhǔn)確的冠幅預(yù)測(cè)模型對(duì)今后專門針對(duì)干旱半干旱生態(tài)系統(tǒng)模型開(kāi)發(fā)和研究是至關(guān)重要的。本研究以古爾班通古特沙漠三種沙丘類型(包括固定沙丘、半固定沙丘和流動(dòng)沙丘)上的主要固沙灌木為研究對(duì)象, 建立了該區(qū)域不同沙丘類型上固沙灌木的冠幅—株高模型, 研究結(jié)果可為今后該區(qū)域相關(guān)生態(tài)模型研究和防風(fēng)固沙效益的評(píng)價(jià)提供科學(xué)的依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    古爾班通古特沙漠研究區(qū)位于新疆維吾爾族自治區(qū)準(zhǔn)噶爾盆地中南部(44°11′—146°20′N, 84°31′—190°00′E), 面積5.11萬(wàn)km2, 海拔高度為300—1600 m, 是我國(guó)第二大沙漠, 也是我國(guó)最大的固定和半固定沙漠。該沙漠沙丘形態(tài)主要呈壟狀,也可見(jiàn)梁窩狀和蜂窩狀, 高度10—50 m。沙漠區(qū)為溫帶干旱荒漠, 是典型的內(nèi)陸干旱氣候, 年平均氣溫7.19 ℃,年降水量100—150 mm, 沙漠腹地僅有70—100 mm, 主要集中于春季, 年蒸發(fā)量在2000 mm以上, 空氣相對(duì)濕度平均為50%—60%[28-32]。該地區(qū)絕大部分為固定和半固定沙丘, 最有代表性的沙丘類型是沙垅, 占沙漠面積的50%以上。該地區(qū)不同類型沙丘上優(yōu)勢(shì)固沙灌木的種類和數(shù)量有顯著的差異[32]: 固定沙丘上主要分布的灌木為梭梭()和白梭梭()等, 半固定沙丘和流動(dòng)沙丘上主要分布的固沙灌木為白梭梭、沙拐棗()、蛇麻黃()和油蒿()等。

    1.2 數(shù)據(jù)來(lái)源

    2013年8月在古爾班通古特沙漠沿地形分別在固定沙丘北沙窩(87o55′12″E, 44o22′24″N)、半固定沙丘石油管理處(88o16′35″E, 44o33′34″N)和流動(dòng)沙丘彩南(88o32′28″E, 44o56′49″N)設(shè)置調(diào)查樣地。三塊樣地均包含迎風(fēng)坡、背風(fēng)坡、丘頂和丘底四種地貌類型, 樣地的大小分別為40 m×148 m、40 m×180 m和40 m×148 m, 將三塊樣地均劃分為4 m×4 m的調(diào)查樣方, 對(duì)樣方內(nèi)的每一株固沙灌木的冠幅和株高進(jìn)行了調(diào)查。分別沿東西和南北兩個(gè)垂直方向測(cè)量每株灌木在兩個(gè)方向上冠幅[10–11, 33], 每株灌木冠幅(CW)的大小等于兩個(gè)方向上冠幅的均值。需要注意的是不論兩個(gè)方向上冠幅的大小是否相等, 灌木垂直投影到地面上的面積均可通過(guò)以冠幅均值為直徑圓的面積計(jì)算。因此, 灌木冠層的大小可以通過(guò)灌木的冠幅計(jì)算。三塊樣地分別收集到2344個(gè)、1651個(gè)和1187個(gè)固沙灌木冠幅和株高的數(shù)據(jù), 三種沙丘上固沙灌木的基本信息見(jiàn)表1。

    1.3 基礎(chǔ)模型

    選擇文獻(xiàn)中12個(gè)常用的冠幅—株高模型作為本研究的基礎(chǔ)模型(表2)[1, 2, 7–8, 10, 22, 34–35], 利用得到的數(shù)據(jù)對(duì)模型進(jìn)行擬合。最優(yōu)模型的評(píng)價(jià)標(biāo)準(zhǔn)為: 決定系數(shù)(2)較大, 均方根誤差(RMSE)、平均絕對(duì)誤差(MAE) 和平均絕對(duì)百分比誤差(MAPE)較小的模型作為最優(yōu)模型。模型評(píng)價(jià)指標(biāo)的計(jì)算公式分別為:

    表1 不同沙丘類型上固沙灌木的基本情況信息表

    注: 表中冠幅、株高數(shù)據(jù)表示均值±標(biāo)準(zhǔn)差, 變異系數(shù)數(shù)據(jù)表示冠幅變異系數(shù)±株高變異系數(shù)。

    表2 冠幅—株高基礎(chǔ)模型

    注: 模型中表示灌木的冠幅,表示每株灌木的株高。

    2 結(jié)果與分析

    利用12個(gè)基礎(chǔ)模型對(duì)不同沙丘類型上的固沙灌木和不同沙丘類型上不同灌木的冠幅進(jìn)行模型擬合。觀察所有擬合結(jié)果(表3—表8)和表1中變異系數(shù)數(shù)據(jù)表明: 模型擬合結(jié)果并沒(méi)有因?yàn)楣诜椭旮邤?shù)據(jù)的離散程度較大而擬合效果差或離散程度較小而擬合效果好, 數(shù)據(jù)的離散程度和模型擬合優(yōu)度之間沒(méi)有明顯的聯(lián)系。由于模擬結(jié)果較多, 模擬結(jié)果僅保留了模型評(píng)價(jià)較好的模型, 具體結(jié)果如下:

    2.1 不同沙丘類型上固沙灌木冠幅—株高預(yù)測(cè)模型

    雖然不同類型沙丘上固沙灌木的種類不盡相同, 但同一類型沙丘上固沙灌木生長(zhǎng)的土壤環(huán)境、氣候特點(diǎn)和光照等都十分相似, 因此可以考慮是否存在一種較優(yōu)模型可以對(duì)同一沙丘類型上所有固沙灌木進(jìn)行擬合。利用選取的12個(gè)基礎(chǔ)模型分別對(duì)研究區(qū)固定沙丘上的主要固沙灌木梭梭和白梭梭以及半固定沙丘和流動(dòng)沙丘上的主要固沙灌木白梭梭、沙拐棗、油蒿和蛇麻黃進(jìn)行模擬研究。結(jié)果表明, 篩選出的基本模型對(duì)固定沙丘和半固定沙丘上的固沙灌木的擬合效果(決策系數(shù)均大于0.7)要優(yōu)于流動(dòng)沙丘上的固沙灌木(決策系數(shù)在0.3左右)(見(jiàn)表3)。值得注意的是, 12種基本模型對(duì)流動(dòng)沙丘上的固沙灌木均沒(méi)有很好的擬合效果??傮w來(lái)看, 12種基本模型中M6、M7、M10和M12模型對(duì)三種沙丘上的主要固沙灌木均有相對(duì)較好的模擬效果。對(duì)固定沙丘上的主要固沙灌木而言, M1、M2、M3、M6、M7、M9、M10和M12八個(gè)模型模擬效果較好。其中, M7對(duì)固定沙丘上的灌木擬合效果最優(yōu)。M7模型常用于描述事物的發(fā)展由萌芽、成長(zhǎng)到飽和的周期過(guò)程, 該曲線的特點(diǎn)是: 初期增長(zhǎng)緩慢, 以后逐漸加快, 當(dāng)?shù)竭_(dá)一定程度后, 增長(zhǎng)率又逐漸下降, 最后接近一條水平線。對(duì)半固定沙丘上的主要固沙灌木而言, M2、M3、M6、M7、M10和M12六個(gè)模型模擬效果較好。其中, M12對(duì)半固定沙丘上的灌木擬合效果最優(yōu)。對(duì)流動(dòng)沙丘上的主要固沙灌木而言, M5(Logarithmic Model)、 M6、M7、M10和M12五個(gè)基本模型模擬效果較好。其中, M10對(duì)流動(dòng)沙丘上的灌木擬合效果最優(yōu), M10曲線比較簡(jiǎn)單, 它無(wú)拐點(diǎn), 相當(dāng)于理想的生長(zhǎng)曲線, 生長(zhǎng)曲線類似于“肩形”, 是一種近似的“S”形, 它適用于描述一開(kāi)始生長(zhǎng)較快, 無(wú)拐點(diǎn)的灌木。三種沙丘上主要固沙灌木的冠幅—株高觀測(cè)值和幾個(gè)較好的擬合曲線如圖1所示。

    表3 三種沙丘上主要固沙灌木的基本模型擬合結(jié)果

    圖1 三種沙丘上主要固沙灌木的冠幅—株高觀測(cè)值和相應(yīng)的擬合曲線

    Figure 1 Observation values of crown width-plant height and fitting curves of sand-fixing shrubs on three types of sand dunes

    2.2 不同固沙灌木在不同沙丘上的冠幅—株高預(yù)測(cè)模型

    2.2.1 梭梭冠幅—株高預(yù)測(cè)模型

    梭梭是研究區(qū)固定沙丘特有的灌木類型, 其冠幅—株高模型的擬合結(jié)果見(jiàn)表4。綜合來(lái)看, M1、M2、M3、M6、M7、M9、M10和M12八個(gè)模型對(duì)梭梭的擬合效果較好(決策系數(shù)均大于0.76), 且誤差較小。其中, M2模型較其他模型擬合效果最優(yōu)。M2為二次函數(shù)曲線且參數(shù)擬合結(jié)果均大于零, 即梭梭最佳的模型曲線為一個(gè)開(kāi)口向上的拋物線形式, 梭梭的冠幅與株高的關(guān)系為: 冠幅隨株高的增大而減低, 降低到一定數(shù)值后又隨株高的增長(zhǎng)而增長(zhǎng)。

    2.2.2 白梭梭冠幅—株高預(yù)測(cè)模型

    白梭梭在研究區(qū)三種沙丘上均有分布, 其冠幅—株高模型的擬合結(jié)果見(jiàn)表5。結(jié)果表明, M1、M2、M3、M6、M7、M9、M10和M12八個(gè)模型對(duì)白梭梭均有擬合效果較好(決策系數(shù)均大于0.7)。其中, M12對(duì)固定沙丘上的白梭梭固沙植被表現(xiàn)最優(yōu), M10對(duì)半固定沙丘和流動(dòng)沙丘上的白梭梭固沙植被表現(xiàn)最優(yōu), 說(shuō)明半固定沙丘和流動(dòng)沙丘的白梭梭較符合一開(kāi)始生長(zhǎng)較快, 無(wú)拐點(diǎn)的曲線模型。

    2.2.3 沙拐棗冠幅—株高預(yù)測(cè)模型

    沙拐棗主要分布在研究區(qū)的半固定沙丘和流動(dòng)沙丘上, 其冠幅—株高模型的擬合結(jié)果見(jiàn)表6。結(jié)果表明, M1、M2、M3、M6、M7、M9、M10和M12八個(gè)模型對(duì)固沙灌木沙拐棗均有較好擬合效果(決策系數(shù)均大于0.4)。其中, M3對(duì)半固定沙丘上的沙拐棗固沙植被表現(xiàn)最優(yōu), 由擬合參數(shù)結(jié)果可知, 最佳模型曲線為上凸型的冪函數(shù)曲線, 前期冠幅隨株高增長(zhǎng)而快速增長(zhǎng), 后期冠幅隨株高增長(zhǎng)而較緩慢增長(zhǎng)。 M1模型對(duì)流動(dòng)沙丘上的沙拐棗表現(xiàn)較好, 即簡(jiǎn)單的一元線性函數(shù)就可以很好的擬合流動(dòng)沙丘上的沙拐棗。

    表4 固沙灌木梭梭的基本模型擬合結(jié)果

    表5 固沙灌木白梭梭的基本模型擬合結(jié)果

    表6 固沙灌木沙拐棗的基本模型擬合結(jié)果

    2.2.4 蛇麻黃冠幅—株高預(yù)測(cè)模型

    固沙植被蛇麻黃的冠幅—株高模型的擬合結(jié)果見(jiàn)表7。結(jié)果表明, M2、M4 (Exponential Model)、M7、M8 (Growth Model)、M11(Hossfeld Model)五個(gè)模型對(duì)半固定沙丘上的蛇麻黃均擬合效果較好(決策系數(shù)均大于0.3)。其中, M11對(duì)半固定沙丘上的蛇麻黃固沙植被表現(xiàn)最優(yōu)。M2、M3、M5、M7、M10和M12六個(gè)模型對(duì)流動(dòng)沙丘上的蛇麻黃均有較好擬合效果較好(決策系數(shù)均在0.5左右)。其中, M5對(duì)流動(dòng)沙丘上的蛇麻黃固沙植被表現(xiàn)最優(yōu), 由參數(shù)擬合結(jié)果可知M5為上凸的遞增曲線。

    2.2.5 油蒿冠幅—株高預(yù)測(cè)模型

    固沙植被油蒿的冠幅—株高模型的擬合結(jié)果見(jiàn)表8。結(jié)果表明, M1、M2、M3、M6、M7、M9、M10和M12八個(gè)模型對(duì)半固定沙丘和流動(dòng)沙丘上的油蒿均擬合效果較好(決策系數(shù)均在0.7左右)。其中, M1即簡(jiǎn)單的一元線性模型對(duì)半固定沙丘上的油蒿擬合最優(yōu), M10對(duì)流動(dòng)沙丘上的油蒿表現(xiàn)最優(yōu)。

    表7 固沙灌木蛇麻黃的基本模型擬合結(jié)果

    表8 固沙灌木油蒿的基本模型擬合結(jié)果

    3 討論

    本研究以古爾班通古特沙漠三種沙丘上主要的固沙灌木為例, 利用12個(gè)常用的冠幅—株高模型選出了對(duì)該地區(qū)三種沙丘上的所有灌木和不同類型的灌木擬合效果較好的基礎(chǔ)模型。結(jié)果表明: 篩選出的基本模型對(duì)固定沙丘和半固定沙丘上固沙灌木的擬合效果要顯著的優(yōu)于流動(dòng)沙丘上的固沙灌木。這一方面可能是由于本研究在固定沙丘和半固定沙丘樣地中調(diào)查的固沙灌木的數(shù)量較流動(dòng)沙丘上固沙灌木的數(shù)量多, 另一方面可能是因?yàn)榱鲃?dòng)沙丘上的固沙灌木遭受了更強(qiáng)烈的干旱、沙蝕、風(fēng)蝕等氣候和環(huán)境因素的影響。因此, 要準(zhǔn)確的建立流動(dòng)沙丘上固沙灌木的冠幅預(yù)測(cè)模型, 今后可能需要在更多采樣的基礎(chǔ)上在預(yù)測(cè)模型中加入氣候和環(huán)境等相關(guān)變量。同時(shí), 兩個(gè)近“S”型的模型M10(Gompertz Model)和M12(HossfeldⅠ Model)對(duì)三種類型沙丘上的固沙灌木均有較好的模擬效果。兩個(gè)模型對(duì)應(yīng)的生長(zhǎng)曲線較Logistic曲線在初始階段更加陡峭, 這表明固沙灌木的生長(zhǎng)過(guò)程均遵循先快后慢的生長(zhǎng)趨勢(shì), 即固沙灌木在生長(zhǎng)初期的生長(zhǎng)速度較快, 到生長(zhǎng)后期固沙灌木的生長(zhǎng)速度會(huì)逐漸放緩, 這也是固沙灌木為適應(yīng)沙區(qū)的氣候特點(diǎn)而演變出來(lái)的生長(zhǎng)機(jī)制。

    對(duì)于不同類型沙丘而言, M6、M7、M10 和M12 四個(gè)模型可以較好的擬合三種沙丘類型上所有固沙灌木。對(duì)于不同類型固沙灌木而言, M1、M2、M3、M6、M7、M9、M10 和M12八個(gè)模型均可以較好擬合所有沙丘類型上的梭梭、白梭梭、沙拐棗和油蒿, M2和M7可以較好擬合所有沙丘類型上的蛇麻黃。梭梭、白梭梭和油蒿模型的擬合優(yōu)度較高均大于0.7, 而沙拐棗和蛇麻黃的擬合優(yōu)度卻在0.4左右。這主要與固沙灌木的類型有關(guān), 由于沙拐棗和蛇麻黃均屬于典型的衰退型部落[32], 測(cè)得數(shù)據(jù)較少且數(shù)據(jù)之間差異性較大, 導(dǎo)致其擬合優(yōu)度較差, 誤差較大。整體來(lái)看, 對(duì)干旱半干旱生態(tài)系統(tǒng)而言, 固沙灌木的冠幅不僅與其所處地貌類型有關(guān), 而且還與固沙灌木的類型有關(guān), 僅用單個(gè)變量株高作為其冠幅的預(yù)測(cè)因子還遠(yuǎn)遠(yuǎn)不夠。今后, 針對(duì)不同干旱區(qū)不同類型的固沙灌木, 需要進(jìn)一步開(kāi)發(fā)適合不同干旱生態(tài)系統(tǒng)中固沙灌木的冠幅預(yù)測(cè)模型并開(kāi)展綜合沙區(qū)生物因子和非生物因子的冠幅預(yù)測(cè)模型的研究。

    4 結(jié)論

    本文研究結(jié)果顯示: 對(duì)于固定沙丘和半固定沙丘而言, 利用株高預(yù)測(cè)冠幅寬度針對(duì)本文選取的大部分模型而言具有較高的擬合優(yōu)度(R2在0.7左右); 流動(dòng)沙丘模型的擬合優(yōu)度均在0.3左右。除固沙灌木蛇麻黃和沙拐棗外其他灌木的模型擬合優(yōu)度也均大于0.7。其中, M7、M12和M10分別為固定沙丘、半固定沙丘和流動(dòng)沙丘上灌木的最優(yōu)冠幅預(yù)測(cè)模型; M2和M12分別為固定沙丘上梭梭和白梭梭的最優(yōu)冠幅預(yù)測(cè)模型; M10、M3、M11和M1

    分別為半固定沙丘上白梭梭、沙拐棗、蛇麻黃和油蒿的最優(yōu)冠幅預(yù)測(cè)模型; M10、M1、M5和M10分別為流動(dòng)沙丘上白梭梭、沙拐棗、蛇麻黃和油蒿的最優(yōu)冠幅預(yù)測(cè)模型。這些模型可用于對(duì)沙漠地區(qū)生態(tài)系統(tǒng)、防風(fēng)固沙效益、固沙灌木的生物量模型和生長(zhǎng)模型的研究。之后的研究需要更多的考慮地理位置、環(huán)境及其他林分因子這些影響因素, 建立更為精確的冠幅預(yù)測(cè)模型。

    [1] IBRAHIM E M, OSMAN E H. Diameter at breast height-crown width prediction models for(DC.) Guill &Perr and Combretum Hartmannianum Schweinf[J]. Research Article, 2014, 3(4): 191–197.

    [2] YANG Yuqing, HUANG Shongming. Allometric modelling of crown width for white spruce by fixed-and mixed-effects models[J]. The Forestry Chronicle, 2017, 93(2): 138–147.

    [3] BUCHACHER R, LEDERMANN T. Interregional crown width models for individual trees growing in pure and mixed stands in Austria[J]. Forests, 2020, doi: 10.3390/f11010114.

    [4] BRCHTOLD W A. Crown-diameter prediction models for 87 species of stand-grown trees in the eastern United States[J]. Southern Journal of Applied Forestry, 2003, 27(4): 269–278.

    [5] 張冬燕, 王冬至, 范冬冬, 等. 不同立地類型華北落葉松人工林冠幅與胸徑關(guān)系研究[J]. 林業(yè)資源管理, 2019, 4: 69–73.

    [6] LEI Yakai, FU Liyong, AFFLECK D, et al. Additivity of nonlinear tree crown width models: Aggregated and disaggregated model structures using nonlinear simultaneous equations[J]. Forest Ecology and Manage-ment, 2018, 427: 372–382.

    [7] XU Hao, SUN Yujun, WANG Xinjie, et al. Linear mixed-effects models to describe individual tree crown width for China-fir in Fujian Province, southeast China[J]. Plos One, 2015, 10(4): e0122257.

    [8] 呂樂(lè), 董利虎, 李鳳日. 黑龍江省東部地區(qū)天然椴樹(shù)單木冠幅預(yù)測(cè)模型[J]. 東北林業(yè)大學(xué)學(xué)報(bào), 2019, 47(7): 37–42.

    [9] 吳云霄, 雷忻. 樹(shù)木冠層消光度與三維綠量的關(guān)系[J]. 生態(tài)科學(xué), 2016, 35(1): 167–171.

    [10] FU Liyong, SHARMA R P, Hao Kaijie, et al. A generalized interregional nonlinear mixed-effects crown width model for Prince Rupprecht larch in northern China[J]. Forest Ecology and Management, 2017, 389: 364–373.

    [11] MARSHALL D D, JOHNSON G P, Hann D W. Crown profile equations for stand-grown western hemlock trees in northwestern Oregon[J]. Canadian Journal of Forest Research, 2003, 33(11): 2059–2066.

    [12] TAHVANAINEN T, FORSSS E. Individual tree models for the crown biomass distribution of Scots pine, Norway spruce and birch in Finland[J]. Forest Ecology and Management, 2008, 255(3/4): 455–467.

    [13] SHARMA R P, VACEK Z, VACEK S. Individual tree crown width models for Norway spruce and European beech in Czech Republic[J]. Forest Ecology and Manage-ment, 2016, 366: 208–220.

    [14] MONSERUD R A, STERBA H. A basal area increment model for individual trees growing in even-and uneven- aged forest stands in Austria[J]. Forest Ecology and Management, 1996, 80(1/3): 57–80.

    [15] CARVALHO J P, PARRESOL B R. Additivity in tree biomass components of Pyrenean oak (Willd)[J]. Forest Ecology and Management, 2003, 179(1/3): 269–276.

    [16] JIANG Lichun, LIU Ruilong. Segmented taper equations with crown ratio and stand density for Dahurian Larch (Larix gmelinii) in Northeastern China[J]. Journal of Forestry Research, 2011, 22(3): 347.

    [17] 周志翔, 徐永榮, 王鵬程, 等. 不同密度濕地松紙漿原料試驗(yàn)林早期冠幅生長(zhǎng)模型研究[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 1998, 17(3): 289–293.

    [18] GILL S J, BIGING G S, Murphy E C. Modeling conifer tree crown radius and estimating canopy cover[J]. Forest Ecology and Management, 2000, 126(3): 405–416.

    [19] BECHTOLD W A. Largest-crown-width prediction models for 53 species in the western United States[J]. Western Journal of Applied Forestry, 2004, 19(4): 245–251.

    [20] 雷相東, 張則路, 陳曉光. 長(zhǎng)白落葉松等幾個(gè)樹(shù)種冠幅預(yù)測(cè)模型的研究[J]. 北京林業(yè)大學(xué)學(xué)報(bào), 2006, 28(6): 75–77.

    [21] CONDES S, STERBA H. Derivation of compatible crown width equations for some important tree species of Spain[J]. Forest Ecology and Management, 2005, 217(2/3): 203– 218.

    [22] 韓艷剛, 雷澤勇, 趙國(guó)軍, 等. 樟子松人工固沙林冠幅——胸徑模型[J]. 干旱區(qū)研究, 2018, 35(5): 1129–1137.

    [23] 張樹(shù)森, 賈煒瑋, 王玉霞, 等. 大興安嶺地區(qū)天然興安落葉松疏開(kāi)木冠幅預(yù)測(cè)模型[J]. 森林工程, 2017, 33(3): 33–38.

    [24] RUSSELL M B, WEISKITTEL A R. Maximum and largest crown width equations for 15 tree species in Maine[J]. Northern Journal of Applied Forestry, 2011, 28(2): 84–91.

    [25] SATTLER D F, LEMAY V. A system of nonlinear simultaneous equations for crown length and crown radius for the forest dynamics model SORTIE-ND[J]. Canadian Journal of Forest Research, 2011, 41(8): 1567–1576.

    [26] 彭娓, 李鳳日, 董利虎. 黑龍江省長(zhǎng)白落葉松人工林單木生長(zhǎng)模型[J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 42 (3): 19–27.

    [27] 李靖宇, 張琇. 騰格里沙漠不同生物土壤結(jié)皮微生物多樣性分析[J]. 生態(tài)科學(xué), 2017, 36(3): 36–42.

    [28] ZHANG Liyun, CHEN Changdu. On the general charac-teristics of plant diversity of Gurbantunggut sandy desert[J]. Acta Ecologica Sinica, 2002, 22(11): 1923–1932.

    [29] QIAN Yibing, WU Zhaoning, ZHANG Liyun, et al. Impact of habitat heterogeneity on plant community pattern in Gurbantunggut Desert[J]. Journal of Geographical Sciences, 2004, 14(4): 447–455.

    [30] CHEN Yanfen, ZHANG Lingwei, SHI Xiang, et al. Life history responses of spring-and autumn-germinated ephemeral plants to increased nitrogen and precipitation in the Gurbantunggut Desert[J]. Science of The Total Environ-ment, 2019, 659: 756–763.

    [31] 傅思華, 胡順軍, 李浩, 等. 古爾班通古特沙漠南緣梭梭()群落優(yōu)勢(shì)植物水分來(lái)源[J]. 中國(guó)沙漠, 2018, 38(5): 1024–1032.

    [32] 石亞飛, 張志山, 黃磊, 等. 古爾班通古特沙漠半固定沙丘植物群落物種組成和種群結(jié)構(gòu)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2016,27(4): 1024–1030.

    [33] BRAGG D C. A local basal area adjustment for crown width prediction[J]. Northern Journal of Applied Forestry, 2001, 18(1): 22–28.

    [34] RAPTIS D, KAZANA V, KAZAKLS A, et al. A crown width-diameter model for natural even-aged black pine forest management[J]. Forests, 2018, 9(10): 610.

    [35] 符利勇, 孫華. 基于混合效應(yīng)模型的杉木單木冠幅預(yù)測(cè)模型[J]. 林業(yè)科學(xué), 2013, 49(8): 65–74.

    Crown prediction models of the sand-fixing shrubs in Gurbantunggut Desert

    LU Lining, ZHANG Dinghai*, LI Gonglin

    Center for Quantitative Biology, College of Science, Gansu Agricultural University, Lanzhou 730070, China

    Taking the main sand-fixing shrubs,,,andon the three types of dunes (the fixed sand dunes, semi-fixed dunes and mobile dunes) in the Gurbantunggut Desert fortheresearchsubjects, 12 basic models were used to establish crown models of shrubs on different types of sand dunes and shrubs. The simulation results show that both fixed and semi-fixed dunes have better fitting goodness than mobile dunes. Meanwhile, two nearly "S" shaped models M7 (Gompertz Model) and M12 (Hossfeld Ⅰ Model) have good simulation effects for sand-fixing shrubs on three types of sand dunes. Eight models, M1 (Linear Model), M2 (Quadratic Model), M3 (Power Model), M6 (Logistic Model), M7, M9 (Gauss Model), M10 (Mitscherlinch Model) and M12 are suitable for simulating the crown models of most sand-fixing shrubs in the study area. However, shrubs on different sand dunes and different types of shrubs have different optimal crown models. The study provides a scientific basis for the future research of ecological models and the evaluation of wind-break and sand-fixation benefits in this region.

    Gurbantunggut Desert; sand-fixing shrubs;sand dune type; crown prediction models

    路麗寧, 張定海, 李功麟. 古爾班通古特沙漠固沙灌木的冠幅預(yù)測(cè)模型[J]. 生態(tài)科學(xué), 2022, 41(5): 81–89.

    LU Lining, ZHANG Dinghai, LI Gonglin. Crown prediction models of the sand-fixing shrubs in Gurbantunggut Desert[J]. Ecological Science, 2022, 41(5): 81–89.

    10.14108/j.cnki.1008-8873.2022.05.010

    Q939.5

    A

    1008-8873(2022)05-081-10

    2020-08-19;

    2020-09-10

    甘肅省科技計(jì)劃項(xiàng)目(21JR7RA831); 國(guó)家自然基金項(xiàng)目(41661022); 甘肅農(nóng)業(yè)大學(xué)盛彤笙創(chuàng)新基金項(xiàng)目(GAU-CX1121); 甘肅農(nóng)業(yè)大學(xué)學(xué)科建設(shè)專項(xiàng)基金項(xiàng)目(GAU-XKJS-2018-144)

    路麗寧(1995—), 女, 河北邢臺(tái)人, 碩士, 主要從事應(yīng)用統(tǒng)計(jì)學(xué)研究, E-mail: 5084231@qq.com

    張定海, 男, 博士, 副教授, 主要從事生態(tài)水文模型研究, E-mail: zhangdh@gsau.edu.cn

    猜你喜歡
    固沙梭梭冠幅
    無(wú)人機(jī)遙感影像提取的單木冠幅數(shù)據(jù)在桉樹(shù)林分蓄積量估測(cè)中的應(yīng)用1)
    梭梭的建筑課
    城市綠地微環(huán)境對(duì)土壤動(dòng)物群落多樣性的影響
    長(zhǎng)青葉
    化學(xué)固沙材料固沙作用機(jī)理研究的必要性
    施肥對(duì)三江平原丘陵區(qū)長(zhǎng)白落葉松人工林中齡林單木樹(shù)冠圓滿度影響
    與生命賽跑的“沙漠植被之王”——梭梭
    疏水締合聚丙烯酰胺固沙劑的制備及性能研究
    沙漠梭梭的守望者
    沙漠鐵路固沙磚設(shè)計(jì)參數(shù)的數(shù)值模擬
    亚洲国产精品一区二区三区在线| 国产亚洲欧美精品永久| 国产99久久九九免费精品| 老司机深夜福利视频在线观看 | 亚洲欧美色中文字幕在线| 国产日韩欧美亚洲二区| 免费黄色在线免费观看| 久久久久网色| 在线 av 中文字幕| 欧美亚洲日本最大视频资源| 丁香六月天网| 精品人妻一区二区三区麻豆| 国产精品麻豆人妻色哟哟久久| 另类亚洲欧美激情| 啦啦啦 在线观看视频| 亚洲国产最新在线播放| 免费av中文字幕在线| 欧美少妇被猛烈插入视频| 国产精品人妻久久久影院| 亚洲av男天堂| 99热网站在线观看| 99久久精品国产亚洲精品| 最新在线观看一区二区三区 | 国产成人精品在线电影| 亚洲一码二码三码区别大吗| 亚洲一码二码三码区别大吗| 欧美中文综合在线视频| 侵犯人妻中文字幕一二三四区| 久久精品熟女亚洲av麻豆精品| 精品免费久久久久久久清纯 | 亚洲婷婷狠狠爱综合网| 免费黄频网站在线观看国产| 国产免费福利视频在线观看| 午夜免费观看性视频| 叶爱在线成人免费视频播放| 日韩不卡一区二区三区视频在线| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 99精国产麻豆久久婷婷| 久久久久国产一级毛片高清牌| 午夜福利影视在线免费观看| 日韩视频在线欧美| 大话2 男鬼变身卡| 国产亚洲午夜精品一区二区久久| 美女国产高潮福利片在线看| 久久狼人影院| 一本一本久久a久久精品综合妖精| 亚洲国产日韩一区二区| 在线精品无人区一区二区三| 久久精品熟女亚洲av麻豆精品| 久久亚洲国产成人精品v| 九草在线视频观看| 男女下面插进去视频免费观看| 午夜日本视频在线| 十八禁人妻一区二区| 亚洲精品自拍成人| 大码成人一级视频| 五月天丁香电影| 巨乳人妻的诱惑在线观看| 免费女性裸体啪啪无遮挡网站| 一级毛片 在线播放| 久久精品国产综合久久久| av线在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| av又黄又爽大尺度在线免费看| 一本大道久久a久久精品| 在线观看人妻少妇| 免费观看av网站的网址| 一本色道久久久久久精品综合| av电影中文网址| 波多野结衣av一区二区av| 丝袜美腿诱惑在线| 最近最新中文字幕免费大全7| 亚洲欧美清纯卡通| 国产av一区二区精品久久| 日韩人妻精品一区2区三区| 激情五月婷婷亚洲| 欧美成人精品欧美一级黄| 在线天堂中文资源库| 看非洲黑人一级黄片| 啦啦啦在线免费观看视频4| 午夜精品国产一区二区电影| 中文字幕另类日韩欧美亚洲嫩草| 午夜激情久久久久久久| 黄网站色视频无遮挡免费观看| 国产毛片在线视频| 久久久精品国产亚洲av高清涩受| 国产精品久久久久成人av| 高清在线视频一区二区三区| 亚洲欧美精品自产自拍| 一边亲一边摸免费视频| 秋霞在线观看毛片| 久久久久精品人妻al黑| 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜爱| 亚洲精品国产av蜜桃| 欧美激情极品国产一区二区三区| 中文字幕色久视频| 欧美日韩精品网址| 欧美黑人欧美精品刺激| bbb黄色大片| netflix在线观看网站| 国产野战对白在线观看| 欧美精品人与动牲交sv欧美| 我的亚洲天堂| 在线观看人妻少妇| 亚洲精品,欧美精品| 亚洲 欧美一区二区三区| 久久精品亚洲熟妇少妇任你| 国产又爽黄色视频| 丝袜美腿诱惑在线| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 精品卡一卡二卡四卡免费| 99久久综合免费| 国精品久久久久久国模美| 啦啦啦啦在线视频资源| 久久久精品免费免费高清| 国产精品 国内视频| www.熟女人妻精品国产| 一二三四在线观看免费中文在| 欧美激情 高清一区二区三区| 如日韩欧美国产精品一区二区三区| 一区二区三区四区激情视频| 亚洲欧美色中文字幕在线| 交换朋友夫妻互换小说| 精品少妇久久久久久888优播| 久久久久久人人人人人| 国产福利在线免费观看视频| 久久精品久久精品一区二区三区| 91成人精品电影| 男女之事视频高清在线观看 | 欧美国产精品一级二级三级| 美女视频免费永久观看网站| 国产精品av久久久久免费| 极品人妻少妇av视频| 高清欧美精品videossex| 久久婷婷青草| 中文字幕人妻熟女乱码| 在线观看免费高清a一片| 国产亚洲精品第一综合不卡| 国产极品天堂在线| 丰满少妇做爰视频| 亚洲精品一二三| 99热国产这里只有精品6| 美国免费a级毛片| 日韩精品免费视频一区二区三区| 又大又黄又爽视频免费| 美女国产高潮福利片在线看| 美女高潮到喷水免费观看| 最新在线观看一区二区三区 | 亚洲精品在线美女| 你懂的网址亚洲精品在线观看| 免费久久久久久久精品成人欧美视频| 国产欧美日韩一区二区三区在线| 亚洲精品aⅴ在线观看| 在线观看国产h片| 久久久久久久国产电影| 考比视频在线观看| 国产成人精品在线电影| 男女免费视频国产| 亚洲熟女精品中文字幕| 国产一区亚洲一区在线观看| 在线观看人妻少妇| 捣出白浆h1v1| 国产精品av久久久久免费| 日本av免费视频播放| 熟女少妇亚洲综合色aaa.| 国产99久久九九免费精品| 亚洲国产精品成人久久小说| 乱人伦中国视频| 九草在线视频观看| 免费高清在线观看视频在线观看| 多毛熟女@视频| 日韩欧美精品免费久久| 亚洲欧美日韩另类电影网站| videos熟女内射| 亚洲中文av在线| 免费人妻精品一区二区三区视频| 在线观看三级黄色| 国产成人免费无遮挡视频| 国产乱来视频区| a级片在线免费高清观看视频| av网站免费在线观看视频| 在线观看免费午夜福利视频| 国产人伦9x9x在线观看| 高清欧美精品videossex| 欧美日韩亚洲高清精品| 日韩 欧美 亚洲 中文字幕| 99国产综合亚洲精品| av片东京热男人的天堂| 国精品久久久久久国模美| 99热全是精品| av视频免费观看在线观看| 最近最新中文字幕大全免费视频 | 高清不卡的av网站| 人妻 亚洲 视频| 看非洲黑人一级黄片| 热99国产精品久久久久久7| 亚洲成人国产一区在线观看 | 亚洲综合精品二区| 尾随美女入室| 777米奇影视久久| 免费看av在线观看网站| 亚洲av欧美aⅴ国产| 国精品久久久久久国模美| 美女中出高潮动态图| 欧美日韩视频精品一区| www日本在线高清视频| 欧美黄色片欧美黄色片| 高清av免费在线| 亚洲国产成人一精品久久久| 婷婷色麻豆天堂久久| 精品少妇内射三级| 国产精品国产av在线观看| 日韩一卡2卡3卡4卡2021年| 久久久久国产精品人妻一区二区| 黄频高清免费视频| 大片免费播放器 马上看| av视频免费观看在线观看| 亚洲国产av影院在线观看| 秋霞在线观看毛片| 人成视频在线观看免费观看| 国语对白做爰xxxⅹ性视频网站| 免费在线观看完整版高清| 亚洲一级一片aⅴ在线观看| 日韩中文字幕欧美一区二区 | 大码成人一级视频| 久久久久久免费高清国产稀缺| 久久久久精品人妻al黑| 亚洲美女视频黄频| 看免费av毛片| 精品少妇一区二区三区视频日本电影 | 国产福利在线免费观看视频| 各种免费的搞黄视频| 欧美日韩国产mv在线观看视频| 国产亚洲一区二区精品| 精品亚洲成国产av| 制服丝袜香蕉在线| 中文字幕色久视频| 国产精品熟女久久久久浪| 国产色婷婷99| 大话2 男鬼变身卡| 性少妇av在线| 热99国产精品久久久久久7| 午夜福利视频在线观看免费| 亚洲人成77777在线视频| 日韩精品有码人妻一区| 老司机深夜福利视频在线观看 | 精品一品国产午夜福利视频| 少妇的丰满在线观看| 如何舔出高潮| 超色免费av| av不卡在线播放| 日本黄色日本黄色录像| 久久热在线av| 成人午夜精彩视频在线观看| av网站在线播放免费| 日韩制服丝袜自拍偷拍| 欧美日韩视频精品一区| 亚洲三区欧美一区| 亚洲一区二区三区欧美精品| 丝袜脚勾引网站| 90打野战视频偷拍视频| 日日摸夜夜添夜夜爱| 蜜桃在线观看..| 各种免费的搞黄视频| 女的被弄到高潮叫床怎么办| 日韩大码丰满熟妇| 亚洲精品日本国产第一区| 看十八女毛片水多多多| xxxhd国产人妻xxx| 国产精品二区激情视频| 成人漫画全彩无遮挡| 国产淫语在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲专区中文字幕在线 | 国产爽快片一区二区三区| 男女边摸边吃奶| 日韩免费高清中文字幕av| 国产欧美亚洲国产| 看免费av毛片| 精品少妇黑人巨大在线播放| 欧美久久黑人一区二区| kizo精华| 激情五月婷婷亚洲| 成人午夜精彩视频在线观看| 免费久久久久久久精品成人欧美视频| 国产又爽黄色视频| 国产一区亚洲一区在线观看| 免费黄频网站在线观看国产| 最新在线观看一区二区三区 | 精品一区二区三卡| 十分钟在线观看高清视频www| 亚洲精品中文字幕在线视频| 街头女战士在线观看网站| 黄色视频在线播放观看不卡| 欧美精品一区二区大全| 一边摸一边抽搐一进一出视频| 黄色一级大片看看| 一本色道久久久久久精品综合| 日本欧美国产在线视频| 涩涩av久久男人的天堂| 宅男免费午夜| 街头女战士在线观看网站| 成年人免费黄色播放视频| 亚洲精品视频女| 久久久久久久精品精品| 亚洲av在线观看美女高潮| 午夜福利免费观看在线| 久久这里只有精品19| 久久99热这里只频精品6学生| 妹子高潮喷水视频| 亚洲中文av在线| 各种免费的搞黄视频| 韩国av在线不卡| 大香蕉久久成人网| 国产精品女同一区二区软件| 久久精品国产a三级三级三级| 免费观看人在逋| 亚洲欧洲精品一区二区精品久久久 | 国产高清国产精品国产三级| 欧美在线黄色| 在线观看三级黄色| 欧美 日韩 精品 国产| 欧美xxⅹ黑人| 丰满迷人的少妇在线观看| 丝袜在线中文字幕| 欧美日韩一区二区视频在线观看视频在线| 丰满迷人的少妇在线观看| 人妻 亚洲 视频| 午夜激情久久久久久久| 桃花免费在线播放| 久久久久网色| av.在线天堂| 好男人视频免费观看在线| 亚洲精品美女久久av网站| 国产日韩一区二区三区精品不卡| 亚洲男人天堂网一区| 国产人伦9x9x在线观看| 免费在线观看完整版高清| av一本久久久久| 欧美黑人精品巨大| 国产视频首页在线观看| kizo精华| 青青草视频在线视频观看| 久久天堂一区二区三区四区| 亚洲国产精品一区三区| 久久亚洲国产成人精品v| 男女高潮啪啪啪动态图| 熟女av电影| 搡老乐熟女国产| 人妻 亚洲 视频| 国产亚洲午夜精品一区二区久久| 久久毛片免费看一区二区三区| 国产精品 欧美亚洲| 亚洲精品久久成人aⅴ小说| 亚洲综合色网址| 欧美精品高潮呻吟av久久| 亚洲一级一片aⅴ在线观看| 国产亚洲精品第一综合不卡| av又黄又爽大尺度在线免费看| tube8黄色片| 一本色道久久久久久精品综合| 国产野战对白在线观看| av又黄又爽大尺度在线免费看| 少妇精品久久久久久久| 亚洲一级一片aⅴ在线观看| 亚洲成人免费av在线播放| 99热网站在线观看| 亚洲欧美色中文字幕在线| 可以免费在线观看a视频的电影网站 | 99精国产麻豆久久婷婷| 婷婷色综合www| 波野结衣二区三区在线| 黄片播放在线免费| 国产精品香港三级国产av潘金莲 | 黄片无遮挡物在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲av在线观看美女高潮| 国产免费福利视频在线观看| √禁漫天堂资源中文www| 日韩制服丝袜自拍偷拍| 久久久久久久久久久久大奶| 日韩成人av中文字幕在线观看| 新久久久久国产一级毛片| a级毛片在线看网站| 精品国产国语对白av| 久久久精品免费免费高清| 亚洲人成电影观看| 久久久久精品性色| 亚洲欧美一区二区三区黑人| 日本欧美国产在线视频| 久久久亚洲精品成人影院| 九色亚洲精品在线播放| 一区二区三区精品91| 性色av一级| 欧美亚洲 丝袜 人妻 在线| 麻豆乱淫一区二区| 欧美日韩视频精品一区| www.自偷自拍.com| 亚洲成国产人片在线观看| 国产xxxxx性猛交| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品熟女亚洲av麻豆精品| 男女边吃奶边做爰视频| 国产成人一区二区在线| 卡戴珊不雅视频在线播放| 老司机影院毛片| 国产在视频线精品| 伊人亚洲综合成人网| 一边摸一边抽搐一进一出视频| a级毛片在线看网站| 亚洲欧美激情在线| 在线观看免费午夜福利视频| videos熟女内射| 精品午夜福利在线看| 男女无遮挡免费网站观看| 熟女少妇亚洲综合色aaa.| av视频免费观看在线观看| 国产一区二区 视频在线| 欧美变态另类bdsm刘玥| 亚洲综合色网址| 新久久久久国产一级毛片| 精品久久蜜臀av无| 精品人妻一区二区三区麻豆| 制服丝袜香蕉在线| 老司机在亚洲福利影院| 一区二区三区精品91| 欧美 日韩 精品 国产| 国产麻豆69| 少妇精品久久久久久久| av一本久久久久| 亚洲,欧美,日韩| 国产av国产精品国产| 欧美最新免费一区二区三区| 成年人午夜在线观看视频| 亚洲激情五月婷婷啪啪| 亚洲精品国产区一区二| 久久久精品国产亚洲av高清涩受| 少妇被粗大猛烈的视频| av线在线观看网站| 一级片'在线观看视频| 国产成人午夜福利电影在线观看| 亚洲色图综合在线观看| 欧美激情高清一区二区三区 | 亚洲人成电影观看| 天天躁夜夜躁狠狠躁躁| 亚洲av福利一区| 自线自在国产av| 精品一区二区免费观看| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 99九九在线精品视频| 天天添夜夜摸| 亚洲欧美中文字幕日韩二区| 亚洲熟女精品中文字幕| 精品亚洲成a人片在线观看| 十分钟在线观看高清视频www| 亚洲精品一区蜜桃| 极品人妻少妇av视频| 久久久久久人人人人人| 欧美日韩视频精品一区| 久久久精品区二区三区| 日韩人妻精品一区2区三区| 久久久久网色| 亚洲国产欧美日韩在线播放| av视频免费观看在线观看| av有码第一页| 搡老乐熟女国产| 久久久精品免费免费高清| 熟女少妇亚洲综合色aaa.| 在线观看一区二区三区激情| 国产成人免费观看mmmm| 韩国精品一区二区三区| 哪个播放器可以免费观看大片| 国产乱来视频区| 亚洲激情五月婷婷啪啪| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人| 我要看黄色一级片免费的| 男女之事视频高清在线观看 | 久久综合国产亚洲精品| 女人久久www免费人成看片| 国产免费视频播放在线视频| 大香蕉久久网| 天堂8中文在线网| 一本一本久久a久久精品综合妖精| 亚洲情色 制服丝袜| 亚洲第一青青草原| 欧美日韩视频高清一区二区三区二| 亚洲国产精品999| 亚洲精品国产一区二区精华液| 国产精品香港三级国产av潘金莲 | 色视频在线一区二区三区| 精品人妻在线不人妻| 国产亚洲午夜精品一区二区久久| 一区二区日韩欧美中文字幕| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 2018国产大陆天天弄谢| 成年动漫av网址| www.精华液| 欧美在线黄色| www.精华液| 免费在线观看完整版高清| 亚洲欧洲精品一区二区精品久久久 | 热99国产精品久久久久久7| 久久久国产精品麻豆| xxxhd国产人妻xxx| 高清视频免费观看一区二区| 亚洲国产精品一区三区| 欧美激情 高清一区二区三区| 国产片内射在线| 欧美激情 高清一区二区三区| 9色porny在线观看| 在线精品无人区一区二区三| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区| 国产亚洲最大av| 免费在线观看完整版高清| 91老司机精品| 99国产精品免费福利视频| 高清不卡的av网站| 久久女婷五月综合色啪小说| 国语对白做爰xxxⅹ性视频网站| 欧美日韩一区二区视频在线观看视频在线| 91精品伊人久久大香线蕉| 国产高清不卡午夜福利| av卡一久久| 精品国产露脸久久av麻豆| 咕卡用的链子| 美女脱内裤让男人舔精品视频| 欧美日韩国产mv在线观看视频| 色网站视频免费| 久久精品亚洲熟妇少妇任你| 一级片免费观看大全| 狠狠精品人妻久久久久久综合| 国产精品免费大片| 国产在线免费精品| a 毛片基地| av福利片在线| 成人午夜精彩视频在线观看| 精品卡一卡二卡四卡免费| 亚洲激情五月婷婷啪啪| 亚洲精品一区蜜桃| 国产日韩一区二区三区精品不卡| 免费观看人在逋| 亚洲av综合色区一区| xxx大片免费视频| 黄片小视频在线播放| 日韩欧美精品免费久久| 成人手机av| 国产成人系列免费观看| 亚洲 欧美一区二区三区| 欧美精品一区二区大全| 丝瓜视频免费看黄片| 国产高清国产精品国产三级| 香蕉丝袜av| 999精品在线视频| 高清av免费在线| 久久久久久久久久久免费av| 国产精品免费视频内射| 少妇人妻久久综合中文| 欧美激情 高清一区二区三区| 又黄又粗又硬又大视频| 制服诱惑二区| 精品一区二区三卡| 91精品伊人久久大香线蕉| 国产黄色免费在线视频| 无遮挡黄片免费观看| 久久青草综合色| 国产女主播在线喷水免费视频网站| 国产爽快片一区二区三区| 国产片内射在线| 免费观看人在逋| 晚上一个人看的免费电影| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 天堂俺去俺来也www色官网| 777米奇影视久久| 伦理电影大哥的女人| 麻豆乱淫一区二区| 日韩中文字幕视频在线看片| 精品酒店卫生间| 侵犯人妻中文字幕一二三四区| 亚洲免费av在线视频| 国产在线视频一区二区| 亚洲欧美一区二区三区国产| 最新在线观看一区二区三区 | 久久女婷五月综合色啪小说| 在线观看www视频免费| 久久久久国产一级毛片高清牌| 伦理电影大哥的女人| 国产精品国产三级国产专区5o| 日韩中文字幕视频在线看片| 女人久久www免费人成看片| 国产精品国产三级国产专区5o| 777米奇影视久久| 国产成人精品在线电影| 超碰97精品在线观看| √禁漫天堂资源中文www| 国产有黄有色有爽视频| 亚洲国产成人一精品久久久| 国产精品一区二区在线不卡| 男女免费视频国产| 亚洲国产精品国产精品| 青青草视频在线视频观看| 电影成人av| 亚洲国产毛片av蜜桃av| 午夜福利视频精品|