• <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ù)值模擬
    国产aⅴ精品一区二区三区波| 真人一进一出gif抽搐免费| 亚洲欧美日韩卡通动漫| АⅤ资源中文在线天堂| 免费电影在线观看免费观看| 国产精品一及| 日本精品一区二区三区蜜桃| 亚洲内射少妇av| 亚洲国产欧美人成| 精品久久久久久久毛片微露脸| 国产伦人伦偷精品视频| 久久久久久久精品吃奶| 欧美xxxx黑人xx丫x性爽| 久久久久九九精品影院| 亚洲黑人精品在线| 国模一区二区三区四区视频| 成人无遮挡网站| 精品欧美国产一区二区三| 1024手机看黄色片| 国产69精品久久久久777片| 特大巨黑吊av在线直播| 日韩欧美在线二视频| xxx96com| 黑人欧美特级aaaaaa片| 99热精品在线国产| 18+在线观看网站| 久久久久国产精品人妻aⅴ院| av福利片在线观看| 国产精品综合久久久久久久免费| 国内精品美女久久久久久| 99国产极品粉嫩在线观看| 国产不卡一卡二| 久久人人精品亚洲av| 叶爱在线成人免费视频播放| 欧美区成人在线视频| 国产精品久久久久久久电影 | 国产av不卡久久| 亚洲欧美日韩卡通动漫| 国产日本99.免费观看| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| 国产av在哪里看| 久久精品影院6| 欧美xxxx黑人xx丫x性爽| 国内精品久久久久久久电影| 国产在视频线在精品| 欧美激情久久久久久爽电影| 老熟妇仑乱视频hdxx| 久久精品国产亚洲av香蕉五月| 国产精品自产拍在线观看55亚洲| tocl精华| 亚洲人成网站在线播| 中文字幕高清在线视频| 中文亚洲av片在线观看爽| 国产成人系列免费观看| 国产成年人精品一区二区| a级毛片a级免费在线| 日韩欧美精品免费久久 | 国产视频内射| 亚洲avbb在线观看| 香蕉久久夜色| 欧美日韩黄片免| 国产成人系列免费观看| 成人特级av手机在线观看| 欧美成人a在线观看| 97超级碰碰碰精品色视频在线观看| 精品一区二区三区av网在线观看| 特级一级黄色大片| 精品电影一区二区在线| 丰满的人妻完整版| 精品久久久久久久毛片微露脸| 亚洲av熟女| 国产淫片久久久久久久久 | 超碰av人人做人人爽久久 | 特级一级黄色大片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 怎么达到女性高潮| 少妇的丰满在线观看| 99热这里只有精品一区| 高潮久久久久久久久久久不卡| 人妻久久中文字幕网| 一级毛片女人18水好多| 国产黄片美女视频| www日本黄色视频网| 国产精品一区二区三区四区久久| 日日干狠狠操夜夜爽| e午夜精品久久久久久久| 熟女电影av网| 又黄又爽又免费观看的视频| 日本黄色视频三级网站网址| 国产男靠女视频免费网站| 国产蜜桃级精品一区二区三区| 嫩草影院精品99| 高清毛片免费观看视频网站| 熟女人妻精品中文字幕| 欧美中文日本在线观看视频| 黑人欧美特级aaaaaa片| a级一级毛片免费在线观看| 一区二区三区激情视频| 亚洲美女黄片视频| 天堂动漫精品| 精品久久久久久久末码| 蜜桃亚洲精品一区二区三区| 免费观看精品视频网站| 99热只有精品国产| 久久伊人香网站| 一级作爱视频免费观看| 欧美一级毛片孕妇| 成人18禁在线播放| 欧美一级a爱片免费观看看| 国产探花极品一区二区| 男人舔奶头视频| 岛国视频午夜一区免费看| 亚洲性夜色夜夜综合| 国产成人av教育| 一级毛片女人18水好多| 国产精品一区二区三区四区免费观看 | 69人妻影院| 亚洲精品影视一区二区三区av| 国产精品嫩草影院av在线观看 | 欧美日韩亚洲国产一区二区在线观看| 亚洲第一电影网av| 尤物成人国产欧美一区二区三区| 亚洲欧美激情综合另类| 真人一进一出gif抽搐免费| 内地一区二区视频在线| 国内久久婷婷六月综合欲色啪| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美三级三区| 免费观看人在逋| 18禁裸乳无遮挡免费网站照片| 美女 人体艺术 gogo| 国产真实乱freesex| 最新美女视频免费是黄的| 精品久久久久久久久久免费视频| 欧美绝顶高潮抽搐喷水| 1000部很黄的大片| 女人高潮潮喷娇喘18禁视频| 一区福利在线观看| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| 亚洲av一区综合| 制服丝袜大香蕉在线| 一夜夜www| 天天躁日日操中文字幕| 亚洲一区高清亚洲精品| 欧美一级a爱片免费观看看| 国产在视频线在精品| 精品人妻一区二区三区麻豆 | 制服人妻中文乱码| 日韩大尺度精品在线看网址| 欧美一级a爱片免费观看看| 免费人成在线观看视频色| 韩国av一区二区三区四区| 色尼玛亚洲综合影院| 久久伊人香网站| av黄色大香蕉| 中文字幕高清在线视频| 精品福利观看| 日韩欧美国产一区二区入口| 一进一出抽搐gif免费好疼| 欧美一区二区精品小视频在线| 国产私拍福利视频在线观看| 一个人免费在线观看的高清视频| 国产av不卡久久| 久久久成人免费电影| 此物有八面人人有两片| 99久久精品热视频| 90打野战视频偷拍视频| 亚洲一区二区三区不卡视频| 天堂影院成人在线观看| 搡老熟女国产l中国老女人| 国产精品亚洲一级av第二区| 午夜福利18| 国产精品三级大全| 亚洲五月婷婷丁香| 日本在线视频免费播放| 国产熟女xx| 99国产极品粉嫩在线观看| 色综合站精品国产| 免费在线观看亚洲国产| 一区二区三区国产精品乱码| 色哟哟哟哟哟哟| 一个人看的www免费观看视频| 国产精品久久久久久久电影 | 老司机午夜十八禁免费视频| 欧美乱码精品一区二区三区| 女人被狂操c到高潮| 欧美极品一区二区三区四区| 国产伦一二天堂av在线观看| h日本视频在线播放| 最近最新中文字幕大全电影3| 中文字幕高清在线视频| h日本视频在线播放| 无遮挡黄片免费观看| 热99在线观看视频| 叶爱在线成人免费视频播放| 无人区码免费观看不卡| 日日摸夜夜添夜夜添小说| 国产极品精品免费视频能看的| 日韩欧美一区二区三区在线观看| av视频在线观看入口| 亚洲成人久久爱视频| 日本一本二区三区精品| 丁香六月欧美| 亚洲人成网站在线播| 99国产精品一区二区三区| 色在线成人网| 欧美乱妇无乱码| 91九色精品人成在线观看| 精品国产美女av久久久久小说| 热99re8久久精品国产| 成年版毛片免费区| 丰满乱子伦码专区| 99热这里只有是精品50| 神马国产精品三级电影在线观看| 免费看美女性在线毛片视频| 成年免费大片在线观看| 免费av观看视频| 日本免费一区二区三区高清不卡| 国产精品日韩av在线免费观看| 日韩亚洲欧美综合| 国产伦精品一区二区三区四那| 国产一区二区三区视频了| 高清日韩中文字幕在线| 日韩亚洲欧美综合| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| 18美女黄网站色大片免费观看| 狂野欧美激情性xxxx| 国产亚洲精品综合一区在线观看| 亚洲乱码一区二区免费版| 九九热线精品视视频播放| 久9热在线精品视频| 最近最新中文字幕大全电影3| 亚洲精品一卡2卡三卡4卡5卡| 精华霜和精华液先用哪个| 免费大片18禁| 女警被强在线播放| 好看av亚洲va欧美ⅴa在| 大型黄色视频在线免费观看| 欧美性猛交╳xxx乱大交人| 国产美女午夜福利| 亚洲一区二区三区不卡视频| 亚洲精品成人久久久久久| 床上黄色一级片| 老司机午夜十八禁免费视频| 老汉色∧v一级毛片| 欧美成人性av电影在线观看| 亚洲午夜理论影院| 国产黄a三级三级三级人| 亚洲人成网站在线播| 亚洲 国产 在线| 亚洲狠狠婷婷综合久久图片| 麻豆国产av国片精品| 成熟少妇高潮喷水视频| 久久午夜亚洲精品久久| 久久久精品大字幕| 无遮挡黄片免费观看| 老鸭窝网址在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 性色av乱码一区二区三区2| 女人被狂操c到高潮| 亚洲欧美日韩东京热| 亚洲国产精品999在线| 亚洲无线在线观看| а√天堂www在线а√下载| 亚洲七黄色美女视频| 国产 一区 欧美 日韩| 人妻夜夜爽99麻豆av| 国产精品三级大全| 女警被强在线播放| 国产欧美日韩精品一区二区| 国产高清视频在线播放一区| 欧美一区二区国产精品久久精品| 很黄的视频免费| 好看av亚洲va欧美ⅴa在| av国产免费在线观看| 老鸭窝网址在线观看| 制服人妻中文乱码| 亚洲成人久久性| 最新在线观看一区二区三区| 少妇的逼水好多| 精品国产三级普通话版| 欧美成人一区二区免费高清观看| 国产久久久一区二区三区| 中文字幕高清在线视频| 在线观看免费午夜福利视频| 观看免费一级毛片| 变态另类成人亚洲欧美熟女| 精品乱码久久久久久99久播| 欧美色视频一区免费| 午夜视频国产福利| 亚洲真实伦在线观看| 国产高清有码在线观看视频| 亚洲国产精品久久男人天堂| 熟女电影av网| 麻豆国产97在线/欧美| 51国产日韩欧美| 一个人看视频在线观看www免费 | 禁无遮挡网站| a级一级毛片免费在线观看| 日韩成人在线观看一区二区三区| 99久久无色码亚洲精品果冻| 两个人看的免费小视频| 国产男靠女视频免费网站| 久久久国产成人免费| 白带黄色成豆腐渣| 18禁美女被吸乳视频| 亚洲第一电影网av| 欧美+亚洲+日韩+国产| 超碰av人人做人人爽久久 | 欧美日韩综合久久久久久 | 午夜福利在线在线| 男女之事视频高清在线观看| 久久久久精品国产欧美久久久| 亚洲人成网站在线播放欧美日韩| 久久久久久国产a免费观看| 美女高潮的动态| 一个人看的www免费观看视频| 丰满人妻一区二区三区视频av | 在线播放无遮挡| 日韩成人在线观看一区二区三区| 国产精华一区二区三区| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 老司机午夜十八禁免费视频| 午夜两性在线视频| 欧美性感艳星| 在线观看舔阴道视频| 国产精品日韩av在线免费观看| 久久久精品大字幕| 国产欧美日韩精品亚洲av| 国产色爽女视频免费观看| 国产高清激情床上av| 中出人妻视频一区二区| 日本五十路高清| 亚洲成人久久爱视频| 三级国产精品欧美在线观看| 亚洲av免费高清在线观看| 久久国产精品影院| 一本久久中文字幕| 精品一区二区三区视频在线观看免费| 国产精华一区二区三区| 在线天堂最新版资源| 啪啪无遮挡十八禁网站| 少妇高潮的动态图| 亚洲国产精品成人综合色| 国产成+人综合+亚洲专区| 天天添夜夜摸| 又黄又粗又硬又大视频| 狠狠狠狠99中文字幕| 国产成+人综合+亚洲专区| 欧美性猛交╳xxx乱大交人| 看黄色毛片网站| 欧美日韩中文字幕国产精品一区二区三区| av专区在线播放| 少妇的逼好多水| 99热这里只有是精品50| 99在线人妻在线中文字幕| 国产色爽女视频免费观看| 国产精品免费一区二区三区在线| 精品久久久久久久久久免费视频| АⅤ资源中文在线天堂| 麻豆一二三区av精品| 床上黄色一级片| 免费人成在线观看视频色| 久久精品影院6| 一级毛片女人18水好多| 国产精品日韩av在线免费观看| 可以在线观看毛片的网站| 国产黄a三级三级三级人| 脱女人内裤的视频| 国产精品三级大全| 久久精品国产综合久久久| 成人特级av手机在线观看| 日日摸夜夜添夜夜添小说| 亚洲精品国产精品久久久不卡| 淫秽高清视频在线观看| 两个人看的免费小视频| АⅤ资源中文在线天堂| 中文字幕av在线有码专区| 有码 亚洲区| 看片在线看免费视频| 国产一区二区亚洲精品在线观看| 国产亚洲精品综合一区在线观看| 制服丝袜大香蕉在线| 欧美乱色亚洲激情| 成人午夜高清在线视频| 一级黄色大片毛片| 亚洲最大成人手机在线| 悠悠久久av| 熟女人妻精品中文字幕| 国产精品久久久久久亚洲av鲁大| 欧美黄色淫秽网站| 午夜a级毛片| 亚洲精品乱码久久久v下载方式 | 国产欧美日韩精品一区二区| 国产一区二区三区在线臀色熟女| 一个人看视频在线观看www免费 | 久久久久久人人人人人| 狂野欧美白嫩少妇大欣赏| 国产精品影院久久| 久久久国产成人精品二区| 亚洲欧美日韩无卡精品| 国产成人啪精品午夜网站| 久久99热这里只有精品18| 国内久久婷婷六月综合欲色啪| 国产精品av视频在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 美女黄网站色视频| 国产欧美日韩精品一区二区| 99国产综合亚洲精品| 一本综合久久免费| 99久久无色码亚洲精品果冻| 在线天堂最新版资源| 日韩亚洲欧美综合| 神马国产精品三级电影在线观看| 久久草成人影院| avwww免费| 在线播放国产精品三级| 欧美黄色片欧美黄色片| 成人欧美大片| 欧美性感艳星| 成年版毛片免费区| 国产亚洲精品久久久com| 身体一侧抽搐| 老司机午夜十八禁免费视频| 3wmmmm亚洲av在线观看| 午夜久久久久精精品| 中文字幕人妻丝袜一区二区| 午夜日韩欧美国产| 久久6这里有精品| 欧美最新免费一区二区三区 | 老司机午夜十八禁免费视频| 婷婷精品国产亚洲av在线| 亚洲一区高清亚洲精品| 婷婷六月久久综合丁香| 在线a可以看的网站| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 精品国产超薄肉色丝袜足j| 欧美zozozo另类| 一本综合久久免费| 午夜免费观看网址| 亚洲狠狠婷婷综合久久图片| 蜜桃亚洲精品一区二区三区| 日本与韩国留学比较| 国产精品国产高清国产av| 99在线人妻在线中文字幕| 亚洲男人的天堂狠狠| 波多野结衣高清无吗| 男人和女人高潮做爰伦理| 国产精品国产高清国产av| 国产亚洲欧美在线一区二区| 亚洲男人的天堂狠狠| 日本五十路高清| 操出白浆在线播放| 久久久久久人人人人人| 天天躁日日操中文字幕| 在线播放无遮挡| 国产高清视频在线观看网站| 又爽又黄无遮挡网站| 成人av一区二区三区在线看| 色精品久久人妻99蜜桃| 伊人久久精品亚洲午夜| 90打野战视频偷拍视频| 亚洲av成人精品一区久久| 一进一出好大好爽视频| 非洲黑人性xxxx精品又粗又长| 一个人看视频在线观看www免费 | 99久久99久久久精品蜜桃| 久久久成人免费电影| 久久久久性生活片| 又紧又爽又黄一区二区| 99久久无色码亚洲精品果冻| 免费人成在线观看视频色| 午夜福利在线在线| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 青草久久国产| 国语自产精品视频在线第100页| 丝袜美腿在线中文| 一进一出好大好爽视频| 久久久久性生活片| 免费看a级黄色片| 国产97色在线日韩免费| av天堂在线播放| 久99久视频精品免费| 动漫黄色视频在线观看| 欧美xxxx黑人xx丫x性爽| 伊人久久精品亚洲午夜| 欧美性猛交黑人性爽| 国产又黄又爽又无遮挡在线| 久久亚洲精品不卡| 色综合欧美亚洲国产小说| 成人三级黄色视频| 成人鲁丝片一二三区免费| 午夜免费成人在线视频| 狂野欧美白嫩少妇大欣赏| 每晚都被弄得嗷嗷叫到高潮| 欧美成狂野欧美在线观看| 在线观看66精品国产| 亚洲aⅴ乱码一区二区在线播放| 免费观看人在逋| netflix在线观看网站| 99热这里只有是精品50| 亚洲精品美女久久久久99蜜臀| 亚洲无线观看免费| 成人无遮挡网站| 午夜福利高清视频| 熟女少妇亚洲综合色aaa.| 五月玫瑰六月丁香| 特级一级黄色大片| 亚洲精品在线美女| 成年免费大片在线观看| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 黄片小视频在线播放| 9191精品国产免费久久| 国产精品久久视频播放| 国产精品久久久久久亚洲av鲁大| 国产精品永久免费网站| 久久久久久大精品| 黄色成人免费大全| 黄色片一级片一级黄色片| 色噜噜av男人的天堂激情| 免费在线观看日本一区| 搡老妇女老女人老熟妇| 色哟哟哟哟哟哟| 免费看日本二区| 久久久久久久久久黄片| 国产精品亚洲美女久久久| 日韩高清综合在线| 日韩精品青青久久久久久| 国产成人啪精品午夜网站| 成年女人永久免费观看视频| 国产乱人伦免费视频| 一进一出抽搐动态| 搡女人真爽免费视频火全软件 | 国产精品影院久久| 尤物成人国产欧美一区二区三区| or卡值多少钱| 三级男女做爰猛烈吃奶摸视频| 午夜激情福利司机影院| 他把我摸到了高潮在线观看| 在线播放无遮挡| 日韩高清综合在线| 午夜福利18| 久久久久久久亚洲中文字幕 | 免费看日本二区| 亚洲18禁久久av| 日本熟妇午夜| 欧美色视频一区免费| 国产真人三级小视频在线观看| 夜夜躁狠狠躁天天躁| 91在线精品国自产拍蜜月 | 欧美3d第一页| 亚洲av免费在线观看| 欧美绝顶高潮抽搐喷水| 99在线视频只有这里精品首页| 3wmmmm亚洲av在线观看| 美女被艹到高潮喷水动态| 老鸭窝网址在线观看| 97碰自拍视频| 国产综合懂色| 欧美一区二区亚洲| 精品国产亚洲在线| 亚洲欧美精品综合久久99| 国产成人福利小说| av福利片在线观看| 国产亚洲欧美98| 欧美黑人欧美精品刺激| 欧美激情久久久久久爽电影| 国产真实乱freesex| 久久亚洲精品不卡| 亚洲国产日韩欧美精品在线观看 | 午夜福利在线在线| 国产高清三级在线| 狠狠狠狠99中文字幕| 久久久久久九九精品二区国产| 欧美另类亚洲清纯唯美| 一个人免费在线观看的高清视频| 熟女电影av网| 性欧美人与动物交配| 亚洲18禁久久av| 婷婷亚洲欧美| 免费在线观看亚洲国产| 国产精品免费一区二区三区在线| 一个人看视频在线观看www免费 | 18禁黄网站禁片免费观看直播| 窝窝影院91人妻| 欧美成人一区二区免费高清观看| 村上凉子中文字幕在线| 夜夜躁狠狠躁天天躁| 啪啪无遮挡十八禁网站| 日本五十路高清| 亚洲在线自拍视频| 久久6这里有精品| 久久欧美精品欧美久久欧美| 国产又黄又爽又无遮挡在线| 狠狠狠狠99中文字幕| 欧美+亚洲+日韩+国产| 国产色爽女视频免费观看| 亚洲av熟女| 亚洲精品在线观看二区| 亚洲av成人精品一区久久| 亚洲国产精品成人综合色| av天堂在线播放| 国内精品久久久久精免费| 亚洲精品影视一区二区三区av| 午夜影院日韩av| 国产欧美日韩精品一区二区| 狂野欧美白嫩少妇大欣赏| 日本黄大片高清|