• <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ù)值模擬
    一二三四社区在线视频社区8| 亚州av有码| av国产免费在线观看| 深夜a级毛片| 久久国产乱子免费精品| 国产精品久久久久久精品电影| 亚洲内射少妇av| 欧美成人性av电影在线观看| 欧美日韩综合久久久久久 | 国模一区二区三区四区视频| 毛片一级片免费看久久久久 | 我的老师免费观看完整版| 国产白丝娇喘喷水9色精品| 少妇高潮的动态图| 国产精品98久久久久久宅男小说| 国产亚洲精品久久久久久毛片| 国产视频一区二区在线看| 身体一侧抽搐| 此物有八面人人有两片| 国产男靠女视频免费网站| 久久精品国产自在天天线| 一进一出好大好爽视频| 哪里可以看免费的av片| bbb黄色大片| 午夜精品一区二区三区免费看| 小说图片视频综合网站| 亚洲自偷自拍三级| 亚洲欧美清纯卡通| 美女被艹到高潮喷水动态| 精品国产亚洲在线| 国产黄a三级三级三级人| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| 亚洲 欧美 日韩 在线 免费| 五月玫瑰六月丁香| 一本一本综合久久| 99久久精品一区二区三区| 精品一区二区三区人妻视频| 亚洲国产欧美人成| 久久久国产成人精品二区| 在线播放无遮挡| 嫩草影视91久久| 亚洲内射少妇av| 麻豆国产av国片精品| 久久国产精品人妻蜜桃| 国产高潮美女av| 国产国拍精品亚洲av在线观看| 亚洲av免费高清在线观看| 老女人水多毛片| 午夜两性在线视频| 九色成人免费人妻av| 午夜精品久久久久久毛片777| 国产成人a区在线观看| 精品国产亚洲在线| 午夜福利免费观看在线| 桃色一区二区三区在线观看| 欧美精品啪啪一区二区三区| 亚洲综合色惰| 简卡轻食公司| 99久久99久久久精品蜜桃| 最近在线观看免费完整版| 亚洲第一电影网av| 国产一级毛片七仙女欲春2| 亚洲av美国av| 国产成人aa在线观看| 直男gayav资源| 一个人看的www免费观看视频| 久久国产乱子免费精品| 久久亚洲真实| 99热这里只有是精品50| 亚洲av成人不卡在线观看播放网| 性色av乱码一区二区三区2| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 午夜福利18| 男女之事视频高清在线观看| 国产精品永久免费网站| 身体一侧抽搐| 国产精品女同一区二区软件 | 搞女人的毛片| av欧美777| 久久亚洲真实| 欧美潮喷喷水| 村上凉子中文字幕在线| 欧美黑人巨大hd| 欧美国产日韩亚洲一区| 国产一区二区三区视频了| 亚洲熟妇中文字幕五十中出| 日本撒尿小便嘘嘘汇集6| 一本综合久久免费| 99久久精品热视频| av女优亚洲男人天堂| 美女被艹到高潮喷水动态| 波多野结衣高清无吗| 黄色视频,在线免费观看| 国产单亲对白刺激| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 老司机午夜十八禁免费视频| 白带黄色成豆腐渣| 91麻豆精品激情在线观看国产| 亚洲成av人片免费观看| 中文字幕久久专区| 日韩免费av在线播放| or卡值多少钱| 欧美3d第一页| 日本成人三级电影网站| 久久精品综合一区二区三区| 国产精品久久久久久久久免 | 久久中文看片网| 免费电影在线观看免费观看| 成人高潮视频无遮挡免费网站| 国产精品一区二区三区四区免费观看 | 老司机午夜福利在线观看视频| 亚洲在线自拍视频| 久久久久国产精品人妻aⅴ院| 悠悠久久av| 欧美在线黄色| 淫秽高清视频在线观看| 久久久久国内视频| aaaaa片日本免费| 看免费av毛片| 亚洲成人免费电影在线观看| 亚洲av免费在线观看| 久99久视频精品免费| 国产探花在线观看一区二区| 精品国产亚洲在线| 国产亚洲精品久久久久久毛片| 成年女人毛片免费观看观看9| 啦啦啦韩国在线观看视频| 国产在线男女| 日日摸夜夜添夜夜添av毛片 | 亚洲男人的天堂狠狠| 国产精品一区二区免费欧美| 欧美+日韩+精品| 99久久久亚洲精品蜜臀av| 3wmmmm亚洲av在线观看| 亚洲内射少妇av| 18美女黄网站色大片免费观看| 久久久久性生活片| 天堂网av新在线| 一个人看的www免费观看视频| 午夜影院日韩av| 哪里可以看免费的av片| 亚洲欧美清纯卡通| 国内久久婷婷六月综合欲色啪| 中文字幕高清在线视频| 99久久99久久久精品蜜桃| 日韩亚洲欧美综合| 老司机午夜福利在线观看视频| 欧美中文日本在线观看视频| 内射极品少妇av片p| 身体一侧抽搐| 99在线人妻在线中文字幕| 在线观看午夜福利视频| 色尼玛亚洲综合影院| 成人高潮视频无遮挡免费网站| 在线观看一区二区三区| 国产高清有码在线观看视频| 亚洲成a人片在线一区二区| 高潮久久久久久久久久久不卡| 国产精品1区2区在线观看.| 午夜影院日韩av| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 午夜久久久久精精品| 国产黄片美女视频| 国产三级中文精品| 丰满的人妻完整版| 久久6这里有精品| 在线看三级毛片| 丁香六月欧美| 午夜两性在线视频| 久久精品久久久久久噜噜老黄 | 日本黄色视频三级网站网址| 日本在线视频免费播放| 国内精品久久久久精免费| 免费无遮挡裸体视频| 一进一出好大好爽视频| 三级毛片av免费| 伊人久久精品亚洲午夜| 哪里可以看免费的av片| 国内久久婷婷六月综合欲色啪| 国产高清视频在线播放一区| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| 久久亚洲精品不卡| 美女 人体艺术 gogo| 白带黄色成豆腐渣| 黄色视频,在线免费观看| 国产av一区在线观看免费| 国产乱人视频| 最后的刺客免费高清国语| 午夜视频国产福利| 色吧在线观看| 午夜福利成人在线免费观看| 国产高清视频在线播放一区| 一边摸一边抽搐一进一小说| 中文资源天堂在线| 波多野结衣高清作品| www.熟女人妻精品国产| 欧美潮喷喷水| 床上黄色一级片| 精品无人区乱码1区二区| 国产欧美日韩精品一区二区| 国产伦一二天堂av在线观看| 美女大奶头视频| 国产精品女同一区二区软件 | 在线看三级毛片| 18美女黄网站色大片免费观看| 成人美女网站在线观看视频| 亚洲精品在线观看二区| 欧美成人性av电影在线观看| 亚洲一区高清亚洲精品| 天天一区二区日本电影三级| 亚洲成人精品中文字幕电影| 噜噜噜噜噜久久久久久91| 少妇人妻精品综合一区二区 | 99riav亚洲国产免费| 国产精品久久电影中文字幕| 成人高潮视频无遮挡免费网站| 深夜a级毛片| or卡值多少钱| 免费在线观看亚洲国产| 91狼人影院| 91久久精品电影网| 九色成人免费人妻av| 在现免费观看毛片| 久久久精品大字幕| 99久久久亚洲精品蜜臀av| 天堂√8在线中文| 俺也久久电影网| 色av中文字幕| 给我免费播放毛片高清在线观看| 蜜桃亚洲精品一区二区三区| 日韩大尺度精品在线看网址| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 国模一区二区三区四区视频| 尤物成人国产欧美一区二区三区| 中出人妻视频一区二区| 亚洲欧美日韩高清在线视频| 日韩欧美 国产精品| 床上黄色一级片| 非洲黑人性xxxx精品又粗又长| 成人永久免费在线观看视频| 精品人妻视频免费看| 69av精品久久久久久| 国产久久久一区二区三区| 亚洲av二区三区四区| 搡女人真爽免费视频火全软件 | 亚洲成人精品中文字幕电影| 美女 人体艺术 gogo| 九九在线视频观看精品| 免费搜索国产男女视频| 国产熟女xx| 久久人人爽人人爽人人片va | av在线观看视频网站免费| 给我免费播放毛片高清在线观看| 麻豆成人av在线观看| 一级毛片久久久久久久久女| 成人美女网站在线观看视频| 日韩人妻高清精品专区| 伊人久久精品亚洲午夜| 精品日产1卡2卡| 欧美日韩国产亚洲二区| 丁香欧美五月| 一边摸一边抽搐一进一小说| 精品一区二区三区av网在线观看| 两个人的视频大全免费| 国产一区二区激情短视频| 午夜激情福利司机影院| 久久久久久久久久成人| 国产不卡一卡二| www.www免费av| 九色国产91popny在线| 两个人的视频大全免费| 亚洲狠狠婷婷综合久久图片| 最新在线观看一区二区三区| 国产毛片a区久久久久| 美女黄网站色视频| 在线免费观看不下载黄p国产 | 国产免费男女视频| 国产精品爽爽va在线观看网站| 国产综合懂色| 久久精品久久久久久噜噜老黄 | 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 赤兔流量卡办理| 黄片小视频在线播放| 特大巨黑吊av在线直播| 午夜影院日韩av| 18禁黄网站禁片午夜丰满| 国产一区二区三区视频了| 欧美xxxx黑人xx丫x性爽| 国产免费男女视频| 亚洲av一区综合| 国产激情偷乱视频一区二区| 精品久久久久久久久av| 午夜激情福利司机影院| 精品人妻偷拍中文字幕| 熟女人妻精品中文字幕| 12—13女人毛片做爰片一| 色精品久久人妻99蜜桃| 亚洲欧美日韩东京热| 欧美性猛交╳xxx乱大交人| 丁香六月欧美| 国产精品三级大全| 99精品在免费线老司机午夜| 床上黄色一级片| 亚洲在线自拍视频| 天堂动漫精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 很黄的视频免费| 我的老师免费观看完整版| 男人的好看免费观看在线视频| 国产精品国产高清国产av| 观看美女的网站| 在线播放无遮挡| 久久精品综合一区二区三区| 看黄色毛片网站| 精品久久国产蜜桃| 可以在线观看毛片的网站| 亚洲,欧美精品.| 精品人妻一区二区三区麻豆 | 黄色日韩在线| 免费搜索国产男女视频| 99久久九九国产精品国产免费| 偷拍熟女少妇极品色| 少妇的逼水好多| 九色成人免费人妻av| 欧美一级a爱片免费观看看| 久久久久久久久大av| 国产久久久一区二区三区| 少妇丰满av| 特大巨黑吊av在线直播| 亚洲国产精品sss在线观看| 97超级碰碰碰精品色视频在线观看| 久久热精品热| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在 | 99热这里只有是精品50| 波多野结衣高清作品| 国产在线精品亚洲第一网站| 亚洲18禁久久av| 观看免费一级毛片| 国产精品,欧美在线| 久久伊人香网站| 午夜福利欧美成人| 精品午夜福利视频在线观看一区| 国产午夜精品久久久久久一区二区三区 | 久久国产乱子伦精品免费另类| 日本撒尿小便嘘嘘汇集6| 久久欧美精品欧美久久欧美| 国产成人av教育| 一区二区三区高清视频在线| 国产成人aa在线观看| 午夜福利高清视频| 美女 人体艺术 gogo| 97碰自拍视频| 美女高潮的动态| 久久人人精品亚洲av| 国产亚洲欧美98| 国产精品一区二区性色av| 欧美成人性av电影在线观看| 久久久国产成人免费| 真人一进一出gif抽搐免费| 看黄色毛片网站| 国产不卡一卡二| 99精品久久久久人妻精品| 成人高潮视频无遮挡免费网站| 亚洲av成人av| 成熟少妇高潮喷水视频| 色哟哟哟哟哟哟| 91久久精品国产一区二区成人| 国产精品野战在线观看| 久久久成人免费电影| 很黄的视频免费| 亚洲成人免费电影在线观看| 亚洲av第一区精品v没综合| 免费观看精品视频网站| 男女下面进入的视频免费午夜| 99热6这里只有精品| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 亚洲在线观看片| 一区福利在线观看| 99热6这里只有精品| 中文亚洲av片在线观看爽| 性色av乱码一区二区三区2| 蜜桃亚洲精品一区二区三区| 国产欧美日韩精品一区二区| 老熟妇乱子伦视频在线观看| 免费大片18禁| www.熟女人妻精品国产| 人人妻,人人澡人人爽秒播| 久久久久久久精品吃奶| 亚洲欧美日韩高清在线视频| 午夜福利欧美成人| 国产中年淑女户外野战色| 国产乱人伦免费视频| 乱人视频在线观看| 亚洲三级黄色毛片| 国产伦精品一区二区三区视频9| 在线看三级毛片| 身体一侧抽搐| 久久午夜福利片| 老司机深夜福利视频在线观看| 国产三级黄色录像| 日本三级黄在线观看| 村上凉子中文字幕在线| 国产人妻一区二区三区在| 国产老妇女一区| 久久久久久久精品吃奶| 99久久99久久久精品蜜桃| 长腿黑丝高跟| 老女人水多毛片| 国产精品av视频在线免费观看| 热99re8久久精品国产| 久久香蕉精品热| 亚州av有码| 搡女人真爽免费视频火全软件 | 久久伊人香网站| 日韩欧美精品v在线| 成人av在线播放网站| av福利片在线观看| 嫩草影院新地址| 真人做人爱边吃奶动态| 在线免费观看的www视频| 欧美日韩亚洲国产一区二区在线观看| 男人舔奶头视频| 精品一区二区三区视频在线| 国产精品亚洲一级av第二区| 亚洲av第一区精品v没综合| 国产亚洲精品综合一区在线观看| 亚洲自偷自拍三级| 又爽又黄a免费视频| 少妇丰满av| 又爽又黄无遮挡网站| 日本五十路高清| av女优亚洲男人天堂| 成年女人看的毛片在线观看| 在线观看午夜福利视频| 看十八女毛片水多多多| 成人特级av手机在线观看| 九色国产91popny在线| 精品久久久久久久久av| 久久久久久久久久成人| 成人无遮挡网站| 999久久久精品免费观看国产| 51国产日韩欧美| 国产午夜精品久久久久久一区二区三区 | 少妇被粗大猛烈的视频| 真实男女啪啪啪动态图| 国产精品日韩av在线免费观看| 久久久久久久久久成人| 午夜老司机福利剧场| 精品一区二区三区人妻视频| 成年免费大片在线观看| 热99在线观看视频| 高潮久久久久久久久久久不卡| 噜噜噜噜噜久久久久久91| 亚洲国产精品合色在线| 午夜精品在线福利| 亚洲美女黄片视频| 男人狂女人下面高潮的视频| 91午夜精品亚洲一区二区三区 | 久久久久免费精品人妻一区二区| 啪啪无遮挡十八禁网站| 午夜福利在线观看免费完整高清在 | 一级黄色大片毛片| 可以在线观看的亚洲视频| 国产野战对白在线观看| 亚洲国产精品久久男人天堂| 精品久久久久久久久av| 成人高潮视频无遮挡免费网站| a在线观看视频网站| 91字幕亚洲| 中出人妻视频一区二区| 久久久久久久久久黄片| 91九色精品人成在线观看| 久久亚洲真实| ponron亚洲| 久久久国产成人精品二区| 他把我摸到了高潮在线观看| 成人精品一区二区免费| 亚洲国产精品成人综合色| 国产一级毛片七仙女欲春2| 日本撒尿小便嘘嘘汇集6| 久久亚洲精品不卡| 丁香六月欧美| 欧美zozozo另类| 国产一区二区三区视频了| 桃色一区二区三区在线观看| 亚洲激情在线av| 99久久精品一区二区三区| 日本a在线网址| 日本黄色视频三级网站网址| 99热这里只有是精品在线观看 | 真人做人爱边吃奶动态| 日韩av在线大香蕉| 国产高清视频在线播放一区| 婷婷丁香在线五月| 欧美精品国产亚洲| 一个人看视频在线观看www免费| 嫁个100分男人电影在线观看| 神马国产精品三级电影在线观看| 免费黄网站久久成人精品 | 成人永久免费在线观看视频| 一本久久中文字幕| 2021天堂中文幕一二区在线观| 亚洲国产精品久久男人天堂| 欧美日韩中文字幕国产精品一区二区三区| 免费无遮挡裸体视频| 99精品久久久久人妻精品| av国产免费在线观看| 夜夜躁狠狠躁天天躁| 俄罗斯特黄特色一大片| 中文亚洲av片在线观看爽| aaaaa片日本免费| 他把我摸到了高潮在线观看| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 欧美又色又爽又黄视频| www.熟女人妻精品国产| 亚洲狠狠婷婷综合久久图片| 91九色精品人成在线观看| 少妇裸体淫交视频免费看高清| 免费观看精品视频网站| 欧美激情久久久久久爽电影| 色综合站精品国产| av专区在线播放| 精品久久久久久久久av| 成人欧美大片| 麻豆国产97在线/欧美| 看免费av毛片| 成年女人永久免费观看视频| 看黄色毛片网站| 99国产综合亚洲精品| 一个人看视频在线观看www免费| 岛国在线免费视频观看| 婷婷亚洲欧美| 一区二区三区免费毛片| 全区人妻精品视频| 好男人在线观看高清免费视频| 成人毛片a级毛片在线播放| 热99re8久久精品国产| 2021天堂中文幕一二区在线观| 国产v大片淫在线免费观看| 草草在线视频免费看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产国拍精品亚洲av在线观看| 精品久久久久久久人妻蜜臀av| or卡值多少钱| 伊人久久精品亚洲午夜| 老鸭窝网址在线观看| 日本黄大片高清| 成年免费大片在线观看| 国产主播在线观看一区二区| 国产av一区在线观看免费| 床上黄色一级片| 国产黄片美女视频| 欧美色欧美亚洲另类二区| 亚洲,欧美精品.| 日韩国内少妇激情av| 国产人妻一区二区三区在| 婷婷六月久久综合丁香| netflix在线观看网站| 国产精华一区二区三区| 无遮挡黄片免费观看| 精品久久久久久久久久免费视频| 成人午夜高清在线视频| 国产精品综合久久久久久久免费| 直男gayav资源| 深爱激情五月婷婷| 人妻夜夜爽99麻豆av| 内地一区二区视频在线| 黄色视频,在线免费观看| 999久久久精品免费观看国产| 色av中文字幕| 露出奶头的视频| 久久精品91蜜桃| 亚洲精品456在线播放app | 一级作爱视频免费观看| 婷婷六月久久综合丁香| 一本综合久久免费| 午夜精品久久久久久毛片777| 国产精品av视频在线免费观看| 中文资源天堂在线| 真人一进一出gif抽搐免费| 日本 av在线| 亚洲无线观看免费| 久9热在线精品视频| 九九久久精品国产亚洲av麻豆| 欧美日韩国产亚洲二区| 99精品久久久久人妻精品| 他把我摸到了高潮在线观看| 欧美绝顶高潮抽搐喷水| 欧美乱妇无乱码| 亚洲真实伦在线观看| 国产久久久一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲精品久久久久久毛片| 日本黄色片子视频| 欧美午夜高清在线| 窝窝影院91人妻| 偷拍熟女少妇极品色| 黄色日韩在线| 别揉我奶头~嗯~啊~动态视频| 午夜免费男女啪啪视频观看 | 床上黄色一级片| 免费在线观看亚洲国产| 三级国产精品欧美在线观看| 国产中年淑女户外野战色|