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

    華北落葉松人工林樹(shù)高和胸徑的坡位差異與坡面尺度效應(yīng)
    ——以六盤(pán)山香水河小流域?yàn)槔?/h1>
    2016-12-20 09:08:19鄧秀秀王云霓王彥輝王忠誠(chéng)于澎濤
    關(guān)鍵詞:坡長(zhǎng)坡位坡頂

    鄧秀秀 ,王云霓 ,王彥輝 ,王忠誠(chéng) ,熊 偉 ,于澎濤 ,張 桐 ,

    (1. 中南林業(yè)科技大學(xué) 林學(xué)院,湖南 長(zhǎng)沙 410004;2. 中國(guó)林業(yè)科學(xué)研究院 森林生態(tài)環(huán)境與保護(hù)研究所 國(guó)家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100091;3. 北京林業(yè)大學(xué) 水土保持學(xué)院,北京 100083)

    華北落葉松人工林樹(shù)高和胸徑的坡位差異與坡面尺度效應(yīng)
    ——以六盤(pán)山香水河小流域?yàn)槔?/p>

    鄧秀秀1,2,王云霓2,王彥輝2,王忠誠(chéng)1,熊 偉2,于澎濤2,張 桐2,3

    (1. 中南林業(yè)科技大學(xué) 林學(xué)院,湖南 長(zhǎng)沙 410004;2. 中國(guó)林業(yè)科學(xué)研究院 森林生態(tài)環(huán)境與保護(hù)研究所 國(guó)家林業(yè)局森林生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100091;3. 北京林業(yè)大學(xué) 水土保持學(xué)院,北京 100083)

    在六盤(pán)山香水河小流域選擇了一個(gè)水平長(zhǎng)398.2 m和生長(zhǎng)33年生華北落葉松人工林的典型坡面,調(diào)查分析了華北落葉松胸徑與樹(shù)高的坡位差異及隨坡面水平長(zhǎng)度變化的空間尺度效應(yīng)。結(jié)果表明:(1)林分的胸徑和樹(shù)高都存在明顯的坡位差異,從坡頂向下均表現(xiàn)為先增加后減小的變化趨勢(shì)。(2)林分平均胸徑、平均樹(shù)高、優(yōu)勢(shì)木平均高均存在坡面尺度效應(yīng),即平均胸徑的坡面滑動(dòng)平均值在離開(kāi)坡頂?shù)?~200 m和200~398 m的水平距離范圍內(nèi)每增加100 m時(shí)增大0.33 cm和減少0.21 cm;平均樹(shù)高的滑動(dòng)平均值在0~67 、67~305、305~398 m的水平距離范圍內(nèi)每增加100 m時(shí)分別減少1.12 m、增大0.31 m、減少0.10 m;優(yōu)勢(shì)木平均高的坡面滑動(dòng)平均值在0~250 m、250~398 m的水平距離范圍內(nèi)每增加100 m時(shí)增大0.42 m和0.07 m。(3)各樣地的平均胸徑、平均樹(shù)高、優(yōu)勢(shì)木平均高與坡面加權(quán)平均值的比值(Y1、Y2、Y3,小數(shù))隨離開(kāi)坡頂相對(duì)水平坡長(zhǎng)(X,小數(shù))變化的數(shù)量關(guān)系為:Y1= -0.312 8X2+ 0.239 1X+0.986 4(R2= 0.56)、Y2=-0.483 4X3+0.304 3X2+0.125 0X + 0.959 2 (R2= 0.43)、Y3= -0.177 5X2+ 0.230 8X+0.941 6(R2= 0.26),藉此可由特定坡位的樣地調(diào)查值推算整個(gè)坡面平均值,實(shí)現(xiàn)樹(shù)木生長(zhǎng)指標(biāo)從樣地到坡面的尺度轉(zhuǎn)換。

    華北落葉松;胸徑;樹(shù)高;坡面變化;尺度效應(yīng);六盤(pán)山

    林分結(jié)構(gòu)指標(biāo)是描述森林質(zhì)量和制定經(jīng)營(yíng)決策的重要基礎(chǔ),是對(duì)更新、生長(zhǎng)、競(jìng)爭(zhēng)、自然稀疏等林分發(fā)育過(guò)程和人類經(jīng)營(yíng)干擾活動(dòng)等的綜合反映[1],同時(shí)又是評(píng)價(jià)森林生態(tài)系統(tǒng)的生物多樣性保護(hù)[2]、水土保持[3]、水文調(diào)節(jié)[4]、土壤改良[5-6]等服務(wù)功能的依據(jù)。不同林分結(jié)構(gòu)會(huì)造成林下灌草生長(zhǎng)[7]、樹(shù)木病蟲(chóng)害抵抗能力[8]、土壤有機(jī)碳含量[9]、水源涵養(yǎng)能力[10]等的顯著差異,并造成森林的生物量、碳密度、木材蓄積量及價(jià)值等的顯著不同。

    胸徑和樹(shù)高是最主要和最基本的林分結(jié)構(gòu)指標(biāo),相關(guān)研究集中在林分胸徑和樹(shù)高的組成分布上[11-15],但這些均是在林分樣地尺度上開(kāi)展的。林木的生長(zhǎng)與立地及氣候條件密切相關(guān)[16]。在西北地區(qū),降水及溫度等氣象條件往往是限制樹(shù)木生長(zhǎng)的主要因子[17],尤其土壤水分具有直接限制作用,但其限制作用具有地形差異[18],如有關(guān)學(xué)者[19]研究了生物多樣性、葉面積指數(shù)、胸高斷面積等林分結(jié)構(gòu)指標(biāo)隨海拔的變化,以及胸徑和樹(shù)高的組成結(jié)構(gòu)[20]與分布范圍[21]的海拔差異。樹(shù)高和胸徑的生長(zhǎng)也存在著明顯的坡位差異,如福建壽寧縣的巨尾桉表現(xiàn)為坡下>坡中>坡上[22]。但是,當(dāng)前對(duì)胸徑和樹(shù)高等樹(shù)木生長(zhǎng)指標(biāo)空間變化的研究多集中在坡位比較和海拔差異方面,缺乏介于樣地尺度和山體(海拔)尺度之間的山坡尺度范圍內(nèi)的變化規(guī)律研究。事實(shí)上,在坡面尺度范圍內(nèi),存在著土壤厚度、土壤水分、土壤養(yǎng)分、微氣候等方面的明顯差異,并可能導(dǎo)致樹(shù)木生長(zhǎng)指標(biāo)的坡面變化,然而目前還未見(jiàn)到樹(shù)木胸徑和樹(shù)高的坡面變化規(guī)律及空間尺度效應(yīng)的研究,這限制著對(duì)坡面森林結(jié)構(gòu)的準(zhǔn)確描述和精細(xì)經(jīng)營(yíng)。

    華北落葉松Larix principis-rupprechtii是我國(guó)北方地區(qū)的主要造林樹(shù)種。在寧夏六盤(pán)山地區(qū),自20世紀(jì)70年代引入后開(kāi)展了大面積造林,已成為當(dāng)?shù)刈钪匾娜斯ち謽?shù)種,發(fā)揮著水源涵養(yǎng)、水土保持、木材生產(chǎn)和風(fēng)景游憩等方面的重要作用。本研究在具半濕潤(rùn)氣候的六盤(pán)山香水河小流域的一個(gè)華北落葉松同齡人工林典型坡面上,對(duì)不同坡位設(shè)立的空間連續(xù)樣地進(jìn)行林分結(jié)構(gòu)調(diào)查,分析胸徑和樹(shù)高的坡位變化及其坡面空間尺度效應(yīng),以期為森林資源清查、森林經(jīng)營(yíng)方案制定、森林生長(zhǎng)估計(jì)和森林生態(tài)水文功能評(píng)價(jià)等的精細(xì)化管理提供科學(xué)基礎(chǔ),并促進(jìn)森林生態(tài)水文研究,理解坡面尺度效應(yīng)和發(fā)展尺度轉(zhuǎn)換理論。

    1 研究區(qū)概況

    香水河小流域地處寧夏固原市涇源縣(E106°09′~ 106°30′,N35°15′~ 35°41′),呈西北—東南走向,海拔范圍為2 060~2 931 m,面積43.73 km2;屬溫帶半濕潤(rùn)氣候,有春寒無(wú)夏、秋短冬長(zhǎng)的特點(diǎn),年均氣溫6.8 ℃,月均最高和最低氣溫分別為7月份的17.4 ℃和1月份的-7.0 ℃。年均降水量636.8 mm,但集中于7~9月份。土壤類型以灰褐土為主,部分為高山草甸土;土壤厚度一般在60 cm以上,枯落物層厚度在3.6 cm以上。

    該小流域內(nèi)森林植被類型豐富,天然林主要樹(shù)種有遼東櫟Quercus liaotungensis、白樺Betula platyphylla、少脈椴Tilia paucicostata;人工林主要為華北落葉松林,占人工林總面積的90%及小流域面積的24%。灌叢主要有甘肅山楂Crataegus kansuensis、華西箭竹Fargesia nitida、蒙古莢迷Viburnum mongolicum、刺薔薇Rosa acicularis等。草本種類主要有鐵桿蒿Artemisa vestita、蕨Pteridium aquilum、東方草莓Fragaria orientalis、苔草Carex gmelinii、艾蒿Artemisa argyi、糙蘇Phlomis umbrosa等。

    2 研究方法

    香水河小流域內(nèi)72%的華北落葉松人工林均分布在半陽(yáng)坡。因此選擇了一個(gè)代表性的東南坡向的半陽(yáng)坡坡面,在2014年生長(zhǎng)季中期進(jìn)行樣帶布設(shè)及林分結(jié)構(gòu)調(diào)查。該坡面斜坡全長(zhǎng)480.6 m,水平坡長(zhǎng)398.2 m,平均坡度33.7°,海拔變化在2 258.5~2 524.8 m之間。從坡頂至坡腳,布設(shè)了一條寬30 m的樣帶,將其劃分為空間連續(xù)的16塊寬30 m和斜坡長(zhǎng)30 m的樣地。

    在坡頂和坡腳各設(shè)一林外自動(dòng)氣象站,將所測(cè)氣象要素根據(jù)各樣地海拔線性插值,求得各樣地的林外氣象要素。

    在各樣地中心用手持GPS測(cè)量并記錄經(jīng)緯度和海拔,用羅盤(pán)儀測(cè)量坡度和坡向。各樣地進(jìn)行每木調(diào)查:用胸徑卷尺測(cè)胸徑,精確到0.01 cm;用測(cè)高儀測(cè)量樹(shù)高和枝下高,精確到0.1 m;從樹(shù)干向東、南、西、北四個(gè)方向拉伸皮尺測(cè)定樹(shù)冠半徑,利用其平方平均值計(jì)算冠幅直徑,精確到0.01 m。用目測(cè)法估測(cè)林冠郁閉度、林下灌木覆蓋度和草本覆蓋度。

    根據(jù)樣地內(nèi)各株樹(shù)木測(cè)定數(shù)據(jù),計(jì)算樣地的胸徑、樹(shù)高、優(yōu)勢(shì)木高的平均值。

    樣地平均胸徑(D,cm)指各樣地中所有樹(shù)木胸徑的平方平均值(式1):

    樣地平均樹(shù)高(H,m)指各樣地內(nèi)所有樹(shù)木樹(shù)高的算術(shù)平均值(式2):

    式中:Di指樣地內(nèi)第i株樹(shù)木的胸徑(cm);n為各樣地的樹(shù)木株數(shù);Hi為樣地內(nèi)第i株樹(shù)木的樹(shù)高(m)。

    選各樣地最高的5株樹(shù)木作為優(yōu)勢(shì)木,用式(2)計(jì)算各樣地的優(yōu)勢(shì)木平均高。

    在得到各樣地平均胸徑和樹(shù)高后,從坡頂樣地到指定坡面水平長(zhǎng)度的樣地,進(jìn)行樣地胸徑和樹(shù)高的各樣地代表坡段水平長(zhǎng)度的加權(quán)平均,作為從坡頂?shù)街付拥氐钠旅婊瑒?dòng)平均胸徑(Ds)和坡面滑動(dòng)平均樹(shù)高(Hs),其中樹(shù)高進(jìn)行算術(shù)平均 (式3),胸徑進(jìn)行平方平均(式4)。當(dāng)所有樣地參與計(jì)算后,即為整個(gè)坡面的平均胸徑和平均樹(shù)高。

    式中:Hj和Dj指不同坡位第j個(gè)樣地的平均樹(shù)高(m)和胸徑(cm);Sj指第j個(gè)樣地代表的水平坡長(zhǎng)(m)。

    用于評(píng)價(jià)樹(shù)高和胸徑的空間尺度效應(yīng)的指標(biāo)是其沿坡滑動(dòng)平均值隨水平坡長(zhǎng)變化的速率,即水平坡長(zhǎng)每增加100 m對(duì)應(yīng)的胸徑和樹(shù)高變化量的絕對(duì)值。

    3 結(jié)果與分析

    3.1 坡面樹(shù)木生長(zhǎng)整體特征及坡位差異

    研究坡面上的華北落葉松人工林均為33年生,各樣地的林下灌木發(fā)育不明顯,覆蓋度在15%左右,主要種類為刺薔薇、華西箭竹、甘肅山楂等。林下草本層發(fā)育明顯,覆蓋度約40%;各樣地草本種類組成相近,以東方草莓和苔草為優(yōu)勢(shì)種。

    坡面上各樣地的基本信息見(jiàn)表1。在整個(gè)坡面上,華北落葉松林的林齡均為33 a,平均胸徑為20.42 cm,平均樹(shù)高為16.9 m;優(yōu)勢(shì)木平均高為20.6 m,郁閉度變化在0.62~0.81之間,平均為0.74;除個(gè)別樣地密度較低外,絕大多數(shù)樣地的林分密度為800~950株·hm-2,平均為891株·hm-2。

    表1 典型坡面上華北落葉松人工林樣地基本信息Table 1 The basic information of sample plots of L. principis-rupprechtii plantation on representative slope

    由于優(yōu)勢(shì)木高較少受密度影響,比平均樹(shù)高更能反映立地質(zhì)量,因此進(jìn)行了不同坡位樣地的胸徑和優(yōu)勢(shì)木高的單因素方差分析,得到結(jié)果見(jiàn)表2。胸徑有如下坡位變化規(guī)律:坡中上>坡中>坡中下>坡上>坡下,坡下樣地除了與坡上樣地?zé)o顯著差異外,與其它樣地均存在顯著差異, 而其他樣地間差異性不顯著。優(yōu)勢(shì)木平均高在不同坡位的變化為坡中>坡中下>坡下>坡中上>坡上,坡上樣地的優(yōu)勢(shì)木平均高與坡中樣地出現(xiàn)顯著的差異。

    表2 不同坡位處的胸徑、優(yōu)勢(shì)木高生長(zhǎng)情況?Table 2 Growing states of DBH and mean dominant tree height at different slope positions

    3.2 樣地樹(shù)木胸徑與樹(shù)高的坡面變化

    坡面上各樣地的平均胸徑和平均樹(shù)高都存在著明顯的差異,其中平均胸徑變化在18.47~21.57 cm之間,變幅為3.10 cm;平均樹(shù)高變化在15.0~18.6 m之間,變幅為3.6 m。各樣地優(yōu)勢(shì)木平均高變化為19.4~22.0 m,變幅為2.6 m。

    由圖1可知,從坡頂?shù)狡履_,隨坡面水平坡長(zhǎng)增加,樣地平均胸徑、平均樹(shù)高和優(yōu)勢(shì)木平均高都呈現(xiàn)先增加后減小的變化趨勢(shì),其中平均胸徑和平均樹(shù)高的最大值分別出現(xiàn)在水平坡長(zhǎng)150、220 m處;優(yōu)勢(shì)木平均高最大值出現(xiàn)在240 m處,且坡位差異小于平均樹(shù)高。

    為定量評(píng)價(jià)林分生長(zhǎng)指標(biāo)的坡位差異,進(jìn)一步分析了各樣地平均胸徑、樹(shù)高、優(yōu)勢(shì)木高與其坡面加權(quán)平均值的比值隨樣地離開(kāi)坡頂?shù)南鄬?duì)水平坡長(zhǎng)的變化(見(jiàn)圖2)。在離開(kāi)坡頂?shù)南鄬?duì)水平坡長(zhǎng)0~0.71范圍內(nèi),平均胸徑的比值先增加后降低,且其值除在相對(duì)水平坡長(zhǎng)0.04之前均大于1;之后從0.71到坡腳,比值持續(xù)下降且都小于1。平均樹(shù)高比值的坡面變化較胸徑一致,從坡頂?shù)诫x開(kāi)坡頂?shù)南鄬?duì)水平坡長(zhǎng)0.22范圍內(nèi),樹(shù)高比值持續(xù)增加但小于1,在0.22~0.80 范圍內(nèi)比值都大于1,呈現(xiàn)先增加后降低的變化;從0.80到坡腳,其值小于1,且持續(xù)下降。樣地優(yōu)勢(shì)木高比值的坡面變化與平均樹(shù)高有明顯差異,雖然整體變化趨勢(shì)一致,但各坡段具體情況不同,在從坡頂?shù)诫x開(kāi)坡頂相對(duì)水平坡長(zhǎng)的0.36范圍內(nèi),其比值都小于1,隨坡位下降而升高;在0.36~0.93范圍內(nèi),比值基本上都大于1,雖然有先增加后降低的變化,但坡位差異較??;在0.93至坡腳的小段范圍內(nèi),其比值又開(kāi)始小于1。

    3.3 樹(shù)木胸徑和樹(shù)高的坡面空間尺度效應(yīng)

    為分析樣地的樹(shù)木胸徑和樹(shù)高隨坡長(zhǎng)增加而變化的尺度效應(yīng),計(jì)算了它們從坡頂向坡腳沿水平坡長(zhǎng)增加的滑動(dòng)平均值,表明存在明顯的坡面尺度效應(yīng)(見(jiàn)圖3)。

    圖1 華北落葉松人工林樣地的平均胸徑、樹(shù)高和優(yōu)勢(shì)木高隨離坡頂水平距離的變化Fig. 1 Variation of mean DBH, tree height and height of dominant trees of Larix principis-rupprechtii plantation plots along horizontal distance from slope top

    隨著坡面水平坡長(zhǎng)增加,胸徑滑動(dòng)平均值先增加后降低。具體表現(xiàn)為:在離開(kāi)坡頂?shù)乃狡麻L(zhǎng)為0~200 m的范圍內(nèi),胸徑滑動(dòng)平均值持續(xù)增加。樣地平均胸徑的坡面滑動(dòng)平均值(y1,cm)與水平坡長(zhǎng)(x,m)的統(tǒng)計(jì)關(guān)系為:

    樣地樹(shù)木平均樹(shù)高的坡面滑動(dòng)平均值隨水平坡長(zhǎng)的變化與胸徑有所不同,由于坡頂樣地樹(shù)高明顯大于其低處幾個(gè)樣地,故而呈現(xiàn)出先降低后增加再降低的變化趨勢(shì)。在離開(kāi)坡頂?shù)乃狡麻L(zhǎng)0~67 m的比值基本上都大于1,雖然有先增加后降低的變化,但坡位差異較??;在0.93至坡腳的小段范圍內(nèi),其比值又開(kāi)始小于1。

    圖2 華北落葉松人工林樣地平均胸徑、樹(shù)高和優(yōu)勢(shì)木高與坡面平均值的比值隨相對(duì)水平坡長(zhǎng)的變化Fig. 2 Variation of ratios of DBH, tree height and dominant tree height of Larix principis-rupprechtii plantation plots to their whole slope average along relative horizontal slope length

    圖3 華北落葉松林樣地胸徑、樹(shù)高和優(yōu)勢(shì)木高的坡面滑動(dòng)平均值隨離開(kāi)坡頂水平距離的變化Fig.3 Variation of moving averages of DBH, tree height and dominant tree height of Larix principis-rupprechtii plantation plots along distance from slope top

    隨著坡面水平坡長(zhǎng)的增加,胸徑滑動(dòng)平均值先增加后降低。具體表現(xiàn)為:在離開(kāi)坡頂?shù)乃狡麻L(zhǎng)為0~200 m的范圍內(nèi),胸徑滑動(dòng)平均值持續(xù)增加。樣地平均胸徑的坡面滑動(dòng)平均值(y1,cm)與水平坡長(zhǎng)(x,m)的統(tǒng)計(jì)關(guān)系為:

    樣地樹(shù)木平均樹(shù)高的坡面滑動(dòng)平均值隨水平坡長(zhǎng)的變化與胸徑有所不同,由于坡頂樣地樹(shù)高明顯大于其低處幾個(gè)樣地,故而呈現(xiàn)出先降低后增加再降低的變化趨勢(shì)。在離開(kāi)坡頂?shù)乃狡麻L(zhǎng)0~67 m的范圍內(nèi),樹(shù)高滑動(dòng)平均值降低,之后在水平坡長(zhǎng)67~300 m范圍內(nèi)持續(xù)增加,增到最大值后又開(kāi)始急劇降低。樣地平均樹(shù)高的坡面滑動(dòng)平均值(y2,cm)與水平坡長(zhǎng)(x,m)的統(tǒng)計(jì)關(guān)系為:

    由于樣地優(yōu)勢(shì)木樹(shù)高較少受密度影響,其坡面滑動(dòng)平均值隨離開(kāi)坡頂?shù)乃狡麻L(zhǎng)增加的變化規(guī)律性明顯強(qiáng)于平均樹(shù)高,表現(xiàn)為隨水平坡長(zhǎng)增加的單調(diào)上升,但在水平坡長(zhǎng)0~250 m范圍內(nèi)為快速上升,在大于250 m的范圍內(nèi)轉(zhuǎn)為緩慢上升或漸趨穩(wěn)定。樣地優(yōu)勢(shì)木樹(shù)高坡面滑動(dòng)平均值(y3,cm)與水平坡長(zhǎng)(x,m)的統(tǒng)計(jì)關(guān)系為:

    4 結(jié)論和討論

    4.1 胸徑和樹(shù)高坡位變化的主要影響因子

    不少學(xué)者研究了胸徑和樹(shù)高這兩個(gè)重要生長(zhǎng)指標(biāo)的坡位差異,如吉靈波等[23]在云南騰沖的禿杉人工林研究表明,胸徑和樹(shù)高的生長(zhǎng)量均隨坡位下降呈遞增趨勢(shì),這主要與地勢(shì)較低的山坡下部匯集了坡上流入的水分和養(yǎng)分以及坡下部的日照相對(duì)短、風(fēng)力較弱導(dǎo)致蒸散量小和濕度大有關(guān)。甘國(guó)勇[24]對(duì)福建多個(gè)林場(chǎng)不同立地的紅豆樹(shù)人工林研究表明,林木的胸徑與樹(shù)高年生長(zhǎng)量從下坡、中坡和上坡依次遞減。這兩個(gè)研究結(jié)果與本文的優(yōu)勢(shì)木高坡位變化較為一致,但沒(méi)有樣地平均胸徑和平均樹(shù)高隨坡位下降而先增加后降低的復(fù)雜變化,這可能與本文研究地區(qū)的環(huán)境限制更突出有關(guān)。目前對(duì)于華北落葉松雖然沒(méi)有直接研究胸徑和樹(shù)高的坡位差異,但已有學(xué)者研究了其生長(zhǎng)量的坡位差異,如袁明龍等[25]在冀北山地的研究表明其材積連年及平均生長(zhǎng)速率均隨坡位上升依次減小,這主要是因在高坡位立地的土層薄和水分養(yǎng)分條件差,導(dǎo)致林木生長(zhǎng)較慢。

    林分的樹(shù)高和胸徑生長(zhǎng)受到土壤養(yǎng)分[26]和水分、林木間競(jìng)爭(zhēng)[27]、氣象(降水、溫度、輻射等)[28]、地形(海拔(坡位)、坡度、坡向)[29-30]等多因素的共同影響。已有研究[31]表明,六盤(pán)山區(qū)土層厚度、坡向、坡度、海拔高度等立地條件對(duì)華北落葉松生長(zhǎng)都有顯著影響,而在本文研究的坡面上,各樣地土壤厚度均在100 cm以上(見(jiàn)表1),能滿足樹(shù)木根系發(fā)育和營(yíng)養(yǎng)吸收的需求;且該坡面上各坡位樣地的土壤物理性質(zhì)變化不大[32],故土壤養(yǎng)分和物理性質(zhì)的坡位差異不會(huì)是影響樹(shù)木生長(zhǎng)特征坡位差異的主要因素,因而土壤水分的坡位差異影響可能比養(yǎng)分更重要,尤其在干旱年份和干旱季節(jié)。在該研究坡面上,各坡位樣地的坡向相同,所以不存在坡向的影響。樣地的海拔高度有較大差異(2276.6~2524.5 m),并由此帶來(lái)氣溫和降水的一定差異;樣地的坡度有明顯差異(26.5~37.2°),加之存在水分和養(yǎng)分的坡面再分配,所以會(huì)存在坡位(海拔)和坡度的樹(shù)木生長(zhǎng)影響。本研究坡面上生長(zhǎng)的是人工同齡純林,不存在樹(shù)種和林齡影響,但密度變化范圍較大(704~948株·hm-2),會(huì)對(duì)樹(shù)木胸徑及樹(shù)高產(chǎn)生一定影響。

    對(duì)地形和氣象因子與林分平均胸徑、平均樹(shù)高、優(yōu)勢(shì)木平均高的相關(guān)分析(見(jiàn)表3)表明,平均胸徑與坡度呈極顯著負(fù)相關(guān),與坡位(海拔)引起的降水和氣溫變化相關(guān)性較大,其中與海拔和降水正相關(guān),與氣溫負(fù)相關(guān)。林分平均樹(shù)高和優(yōu)勢(shì)木平均高與各氣象因子及地形因子均無(wú)顯著相關(guān),只是胸徑與樹(shù)高、樹(shù)高與優(yōu)勢(shì)木高間存在極顯著和顯著正相關(guān)。

    表3 坡面林分樣地的胸徑和樹(shù)高與地形、氣象、植被因子的相關(guān)分析?Table 3 Correlation of DBH and tree height of stand plots with factors of landform, climate and vegetation at the slope studied

    綜合來(lái)看,各樣地的坡度和坡位是造成樹(shù)木生長(zhǎng)坡位差異的主要因素,其直接影響可能在于它們會(huì)引起水分和養(yǎng)分的坡面再分配,以及坡位(海拔)會(huì)引起氣溫和降水的差異。然而,由于沒(méi)有進(jìn)行土壤水分、養(yǎng)分和各樣地氣象條件的長(zhǎng)期監(jiān)測(cè),不能進(jìn)行定量分析,需在未來(lái)結(jié)合森林生態(tài)定位研究揭示其數(shù)量關(guān)系和作用機(jī)理。

    4.2 坡面樣地樹(shù)高和胸徑變化的尺度效應(yīng)

    在分析坡面樣地的平均胸徑、樹(shù)高和優(yōu)勢(shì)木高的空間尺度效應(yīng)時(shí),將其順坡滑動(dòng)平均值隨水平坡長(zhǎng)增加的變化速率作為評(píng)價(jià)指標(biāo)。圖3表明,坡面樣地平均胸徑的坡面滑動(dòng)平均值隨水平坡長(zhǎng)增加呈現(xiàn)先增加后降低的非單調(diào)變化;平均樹(shù)高的坡面滑動(dòng)平均值表現(xiàn)出比胸徑更復(fù)雜的變化曲線,兼有一個(gè)低谷和一個(gè)高峰;優(yōu)勢(shì)木平均高的坡面滑動(dòng)平均值隨水平坡長(zhǎng)增加呈非線性單調(diào)增加,先是升高較快,后趨于緩慢增加或穩(wěn)定在最大值附近。表明在不同水平坡長(zhǎng)范圍內(nèi)的坡面尺度效應(yīng)存在明顯差異,具體表現(xiàn)為:平均胸徑的坡面滑動(dòng)平均值在離開(kāi)坡頂?shù)?~200 m和200~398 m(至坡腳)的水平距離范圍內(nèi)每增加100 m時(shí)增大0.33 cm和減少0.21 cm;平均樹(shù)高的滑動(dòng)平均值在0~67、67~305、305~398 m(至坡腳)的水平距離范圍內(nèi)每增加100 m時(shí)分別減少1.12 m、增大0.31 m、減少0.10 m;優(yōu)勢(shì)木平均高的坡面滑動(dòng)平均值在0~250、250~398 m(至坡腳)的水平距離范圍內(nèi)每增加100 m時(shí)增大0.42 m和0.07 m。這主要是由于坡度、坡位等的差異,造成土壤含水量的坡面再分配的結(jié)果[34]。在本研究中,坡中位置的坡度分別小于坡中上和坡中下的位置,因此更利于土壤水分的聚集,從而使樹(shù)木生長(zhǎng)較好,使滑動(dòng)平均值在坡中位置相應(yīng)增大。

    4.3 樣地值推算坡面均值的尺度轉(zhuǎn)換途徑

    胸徑和樹(shù)高是森林經(jīng)營(yíng)管理的最常用生長(zhǎng)指標(biāo),也是測(cè)算森林生物量、蓄積量的基礎(chǔ)指標(biāo)。由于存在坡面森林生長(zhǎng)的明顯差異,需對(duì)坡面森林的胸徑和樹(shù)高進(jìn)行精確測(cè)算與估計(jì)。由于實(shí)際森林調(diào)查中不可能在坡面上設(shè)立很多樣地進(jìn)行全面調(diào)查,特別是坡面很長(zhǎng)時(shí),因此往往將典型樣地的調(diào)查值視為整個(gè)坡面的平均值。本研究表明,華北落葉松林坡面樣地的胸徑、樹(shù)高、優(yōu)勢(shì)木高與其坡面加權(quán)平均值的比值存在非線性的沿坡變化。與坡面平均值的比值為1的代表性點(diǎn)位對(duì)胸徑為離開(kāi)坡頂 15 m(相對(duì)坡長(zhǎng)0.04)和274 m(相對(duì)坡長(zhǎng)0.71)處,對(duì)樹(shù)高為85 m(相對(duì)坡長(zhǎng)0.22)和309 m(相對(duì)坡長(zhǎng)0.80)處,對(duì)優(yōu)勢(shì)木高為139 m(相對(duì)坡長(zhǎng)0.36)和359 m(相對(duì)水平坡長(zhǎng)0.93)處,這些特殊坡位樣地的調(diào)查值可代表整個(gè)坡面。與常規(guī)想象的不同,代表性坡位不是在坡面中部,而是隨生長(zhǎng)指標(biāo)不同出現(xiàn)在坡面中下部或中上部。

    研究坡面上各樣地的平均胸徑、平均樹(shù)高和優(yōu)勢(shì)木平均高與其坡面平均值的比值(Y1、Y2、Y3,小數(shù))隨相對(duì)水平坡長(zhǎng)(樣地離開(kāi)坡頂水平距離/整個(gè)坡面水平坡長(zhǎng))(X,小數(shù))變化的回歸關(guān)系式分別為:Y1= -0.312 8X2+ 0.239 1X+0.98 64(R2= 0.56)、Y2=-0.483 4X3+ 0.304 3X2+0.125 0X +0.959 2(R2= 0.43)、Y3= -0.177 5X2+ 0.230 8X +0.941 6(R2= 0.26)。與同處六盤(pán)山區(qū)但具有半干旱氣候的疊疊溝小流域的坡面華北落葉松林的生物量及葉面積指數(shù)相比[33],香水河小流域坡面華北落葉松林的樹(shù)高和胸徑與坡面平均值的比值的坡位變化回歸關(guān)系式的R2值較低,是由于香水河小流域處于半濕潤(rùn)區(qū),土壤水分的坡位差異較小,導(dǎo)致樹(shù)木生長(zhǎng)的坡位差異也較小。即便如此,仍可利用這些關(guān)系將不同長(zhǎng)度坡面上不同坡位的樣地調(diào)查值推算出整個(gè)坡面的平均值。當(dāng)然,如果應(yīng)用地區(qū)的氣候條件、坡面長(zhǎng)度等與本文調(diào)查坡面相差較大,可能會(huì)因坡面水分再分配格局不同而產(chǎn)生較大估計(jì)偏差,因此一方面需在推廣應(yīng)用時(shí)有所限制,一方面需未來(lái)在更多地域和更多坡面開(kāi)展研究,以尋找樹(shù)木生長(zhǎng)指標(biāo)的坡位差異和尺度效應(yīng)與尺度轉(zhuǎn)換的普遍規(guī)律。

    4.4 結(jié) 論

    (1)生長(zhǎng)指標(biāo)存在坡位差異

    在整個(gè)調(diào)查坡面上,華北落葉松人工林的胸徑、樹(shù)高、優(yōu)勢(shì)木高平均值分別為20.42 cm、16.9 m、20.6 m,其變化范圍分別為18.47~21.57 cm、15.0~18.6 m、19.4~22.0 m。由于坡位和坡度等地形因子導(dǎo)致了土壤水分和養(yǎng)分的坡面再分配及降水和溫度等氣象條件的坡位差異,使得林分生長(zhǎng)指標(biāo)存在明顯坡面變化,表現(xiàn)為隨坡位下降先增加后減小的整體變化趨勢(shì),但具體指標(biāo)有所差異,平均胸徑、平均樹(shù)高和優(yōu)勢(shì)木平均高的最大值分別出現(xiàn)在水平坡長(zhǎng)150、220、240 m處。

    (2)生長(zhǎng)指標(biāo)具有坡面尺度效應(yīng)

    林分的胸徑、樹(shù)高以及優(yōu)勢(shì)木高的平均值存在明顯的坡面尺度效應(yīng),且在不同坡段表現(xiàn)不同。平均胸徑在0~200、200~398 m(至坡腳)范圍內(nèi),離開(kāi)坡頂水平距離每增加100 m,滑動(dòng)平均值分別增加0.33 cm和減少0.21 cm;平均樹(shù)高在0~67、67~305、305~398 m(至坡腳)范圍內(nèi),水平距離每增加100 m,滑動(dòng)平均值分別減少1.12 m、增加0.31 m、減少0.10 m;優(yōu)勢(shì)木平均高在0~250 m、250~398 m范圍內(nèi),水平距離每增加100 m,滑動(dòng)平均值分別增加0.42 m和0.07 m。

    (3)可實(shí)現(xiàn)樣地到坡面的尺度轉(zhuǎn)換

    建立了不同坡位各樣地的胸徑、樹(shù)高、優(yōu)勢(shì)木高的平均值與坡面加權(quán)平均值的比值隨離開(kāi)坡頂相對(duì)水平坡長(zhǎng)的關(guān)系式,藉此可由特定樣地的測(cè)定值推算整個(gè)坡面的平均值及其它樣地的數(shù)值,從而實(shí)現(xiàn)從樣地到坡面的尺度轉(zhuǎn)換。

    [1]寧 波. 樟子松人工林結(jié)構(gòu)動(dòng)態(tài)及生物量的研究 [D]. 哈爾濱:東北林業(yè)大學(xué), 2007.

    [2]Nagaike T, Kamitani T, Nakashizuka T. Plant species diversity in abandoned coppice forests in a temperate deciduous forest area of central Japan [J]. Plant Ecology, 2003,166: 45-156

    [3]楊 紅, 周 剛, 傅紹春, 等. 湘中丘陵不同林分結(jié)構(gòu)水土流失特征分析 [J]. 湖南林業(yè)科技, 1996, 23(4): 17-21.

    [4]高 崗. 以水源涵養(yǎng)為目標(biāo)的低功能人工林更新技術(shù)研究 [D].呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學(xué), 2009.

    [5]黨宏忠, 周澤福, 趙雨森, 等. 祁連山水源涵養(yǎng)林土壤水文特征研究 [J]. 林業(yè)科學(xué)研究, 2006, 19(1): 39-44.

    [6]吳家勝, 曾燕如, 黎章矩. 油茶豐產(chǎn)林土壤肥力與林分結(jié)構(gòu)調(diào)查 [J]. 北京林業(yè)大學(xué)學(xué)報(bào), 2009, 31(6): 203-208.

    [7]林 霞, 劉洪見(jiàn), 趙 雪, 等. 不同針闊葉林分結(jié)構(gòu)對(duì)方枝野海棠生長(zhǎng)的影響 [J]. 浙江林業(yè)科技, 2009, 29(3): 30-34.

    [8]李順平, 張 力. 不同林分結(jié)構(gòu)對(duì)天??剐缘恼{(diào)查研究 [J]. 寧夏農(nóng)林科技, 2011, 52(9): 41-42.

    [9]劉夢(mèng)云, 常慶瑞, 齊雁冰, 等. 黃土臺(tái)塬不同林分結(jié)構(gòu)土壤有機(jī)碳質(zhì)量分?jǐn)?shù)特征 [J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2012, 18(6):1418-1427.

    [10]魯紹偉, 劉鳳芹, 余新曉. 北京山區(qū)不同密度油松結(jié)構(gòu)與功能研究 [J]. 水土保持研究, 2007, 14(4): 46-51.

    [11]Bailey R L. Individual tree growth derived from diameter distribution models [J]. Forest Science, 1980, 26: 626-632.

    [12]Kilkki P, Maltamo M, Mykkanen R, et al. Use of the Weibull function in estimating the basal area DBH-distribution [J]. Silva Fennica, 1989, 23: 311-318.

    [13]Borders B E, Patterson W D. Projecting stand tables: A comparison of the Weibull diameter distribution method, a percentile-based projection method, and a basal area growth projection method [J]. Forest Science, 1990, 36(12): 413-424.

    [14]楊艷鋒, 鄭小賢, 梁 雨, 等. 北京八達(dá)嶺林場(chǎng)元寶楓人工林林分結(jié)構(gòu)研究 [J]. 林業(yè)資源管理, 2008, (2): 57-60.

    [15]王春香, 張秋良, 春 蘭, 等. 大青山落葉松人工林直徑分布規(guī)律的研究 [J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版, 2011, 42(3):349-355.

    [16]李書(shū)靖, 黨宏忠, 何虎林, 等. 華北落葉松生長(zhǎng)量與立地氣候因子關(guān)系的研究 [J]. 甘肅林業(yè)科技, 1999, 24(1): 2-6.

    [17]管 偉, 熊 偉, 王彥輝, 等. 六盤(pán)山北側(cè)華北落葉松樹(shù)干直徑生長(zhǎng)變化及其對(duì)環(huán)境因子的響應(yīng) [J]. 林業(yè)科學(xué), 2007,43(9): 1-6.

    [18]Cavelier J. Environmental factors and ecophysiological processes along altitudinal gradients in wet tropical mountains [M].Tropical Forest Plant Ecophysiology, Publisher: Chapman and Hall, 1996: 399-439.

    [19]Homeier J, Breckle S W, Gu¨nter S, et al. Tree diversity, forest structure and productivity along altitudinal and topographical gradients in a Species-Rich Ecuadorian Montane Rain Forest [J].Biotropica, 2010, 42(2): 140-148.

    [20]鄧?yán)っ? 邵 彬, 李 飛. 長(zhǎng)白山北坡云冷杉林胸徑、樹(shù)高結(jié)構(gòu)及其生長(zhǎng)規(guī)律的分析 [J]. 資源科學(xué), 1999, 21(1): 77-84.

    [21]馮愿楠. 晉西黃土區(qū)林草植被特征數(shù)量化研究 [D]. 北京:北京林業(yè)大學(xué), 2008.

    [22]韋 華. 巨尾桉生長(zhǎng)量與不同坡向、坡位關(guān)系的研究 [J]. 亞熱帶植物科學(xué), 2011, 40(4): 50-51.

    [23]吉靈波, 許彥紅, 李 驕, 等. 騰沖縣禿杉人工林立地條件與林分生長(zhǎng)關(guān)系分析 [J]. 林業(yè)調(diào)查規(guī)劃, 2014, 39(2): 147-154.

    [24]甘國(guó)勇. 不同立地質(zhì)量紅豆樹(shù)人工林造林效果分析 [J]. 福建熱作科技, 2011, 36(1): 8-11.

    [25]袁明龍, 姜 韜, 劉 琳, 等. 冀北山地華北落葉松人工林與天然次生油松林不同坡位生長(zhǎng)狀況研究 [J]. 河北林果研究,2011, 26(4): 327-333.

    [26]歐光龍, 王俊峰, 胥 輝, 等. 思茅松天然林胸徑與樹(shù)高結(jié)構(gòu)的變化 [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2014, 34(1): 37-41.

    [27]諶紅輝, 方升佐, 丁貴杰, 等. 馬尾松間伐的密度效應(yīng) [J]. 林業(yè)科學(xué), 2010, 46(5): 84-91.

    [28]姜榮春,王國(guó)蘭. 落葉松胸徑生長(zhǎng)量與氣象因子相關(guān)的初步研究 [J]. 吉林林業(yè)科技, 1985, (1):30-36.

    [29]王向榮, 孫海龍, 余 鑫, 等. 坡向和坡位對(duì)水曲柳中齡林生長(zhǎng)的影響 [J]. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版, 2011, 31(1):30-34.

    [30]王秀花, 馬麗珍, 馬雪紅, 等. 木荷人工林生長(zhǎng)和木材基本密度 [J]. 林業(yè)科學(xué), 2011, 47(7): 138-144.

    [31]姜華偉, 丁 保, 于培彥. 六盤(pán)山地區(qū)華北落葉松生長(zhǎng)與立地條件的關(guān)系分析 [J]. 寧夏農(nóng)林科技, 1988, (6): 16-20.

    [32]曹恭祥. 六盤(pán)山香水河小流域植被結(jié)構(gòu)水文影響及其坡面尺度效應(yīng) [D]. 北京:中國(guó)林業(yè)科學(xué)研究院, 2014.

    [33]李振華. 六盤(pán)山疊疊溝典型植被蒸散及水文要素的坡面尺度效應(yīng) [D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2014.

    [34]韓新生, 鄧?yán)蛱m, 王彥輝, 等. 六盤(pán)山疊疊溝華北落葉松人工林地上生物量的坡面變化[J]. 林業(yè)科學(xué),2015,51(3):132-139.

    Slope variation and scale effect of tree height and DBH of Larix principisrupprechtii plantations along a slope: A case study of Xiangshui watershed of Liupan Mountains

    DENG Xiu-xiu1,2, WANG Yun-ni2, WANG Yan-hui2, WANG Zhong-cheng1, XIONG Wei2, YU Peng-tao2, ZHANG Tong2,3
    (1. Colleague of Forestry, Central South University of Forestry and Technology, Changsha 410004, Hunan, China; 2. Institute of Forest Ecology, Environment and Protection, Key Laboratory of Forestry Ecology and Environment of State Forestry Administration, Chinese Academy of Forestry, Beijing 100091, China; 3. School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China)

    This paper selected a representative slope with a horizontal length of 398.2 m and covered by 33 years-old Larix principisrupprechtii plantation in the small watershed of Xiangshuihe of Liupan Mountains, the variation of tree height and DBH along the slope length and their scale effect were investigated and analyzed. Results showed that: (1) Both the DBH and tree height presented a changing tendency of fi rstly increase and then decrease along the slope from top to foot. (2) The regression relations showing the variation of the moving averages of mean DBH, mean tree height and mean height of dominant trees with increasing horizontal distance from slope top were fi tted. With these relations, the slope scale effect of the three growth parameters was calculated. With a 100 m increase of the horizontal distance downwards, the moving average of the stand mean DBH increased by 0.33 m in the horizontal distance range of 0~200 m, and decreased by 0.21 m in the range of 200~398 m; the moving average of tree height decreased by 1.12 m in the range of 0~67 m, increased by 0.31 m in the range of 67~305 m, and decreased by 0.10 m in the range of 305~398 m; the moving average of the height of dominant trees increased by 0.42 m in the range of 0~250 m, and increased by 0.07 m in the range of 250~398 m. (3)The relations were fi tted to express the slope variation of the ratios of the plots means of DBH, tree height and the height of dominant trees(Y1, Y2, Y3, decimal) to their slope average with the relative horizontal distance from slope top (X, decimal): Y1= -0.312 8X2+ 0.239 1X+0.986 4 (R2= 0.56), Y2=-0.483 4X3+0.304 3X2+0.125 0X+ 0.959 2 (R2= 0.43), Y3= -0.177 5X2+0.230 8X +0.941 6 (R2= 0.26). Using these relations, the whole slope averages of stand growth parameters can be calculated based on the measured values of a plot with a certain slope position, and the up-scaling of stand growth parameters can be realized.

    Larix principis-rupprechtii; DBH; tree height; slope variation; scale effect; Liupan Mountains

    S718.51

    A

    1673-923X(2016)05-0121-08

    10.14067/j.cnki.1673-923x.2016.05.022

    2015-07-10

    科技部“十二五”農(nóng)村領(lǐng)域國(guó)家科技計(jì)劃(2012BAD22B03);國(guó)家自然科學(xué)基金項(xiàng)目(41230852、41390461、41471029);長(zhǎng)沙市典型森林生態(tài)系統(tǒng)固土保肥價(jià)值研究;國(guó)家林業(yè)局寧夏六盤(pán)山森林生態(tài)站資助目

    鄧秀秀,碩士研究生 通訊作者:王忠誠(chéng),副教授,博士;E-mail:wzc366@163.com

    鄧秀秀,王云霓,王彥輝,等. 華北落葉松人工林樹(shù)高和胸徑的坡位差異與坡面尺度效應(yīng)——以六盤(pán)山香水河小流域?yàn)槔齕J].中南林業(yè)科技大學(xué)學(xué)報(bào),2016, 36(5): 121-128.

    [本文編校:謝榮秀]

    猜你喜歡
    坡長(zhǎng)坡位坡頂
    地邊截水地物對(duì)黑土區(qū)小流域坡長(zhǎng)因子計(jì)算的影響
    礦車路線迷宮
    礦車路線迷宮
    坡長(zhǎng)對(duì)貴州喀斯特區(qū)黃壤坡耕地土壤侵蝕的影響
    坡向坡位及郁閉度對(duì)森林公園內(nèi)林下南方紅豆杉生長(zhǎng)的影響
    綠色科技(2017年1期)2017-03-01 10:14:33
    不同坡位對(duì)毛竹林小氣候及地上生長(zhǎng)的影響
    種植地坡向和坡位對(duì)楸樹(shù)引種早期生長(zhǎng)的影響
    橋梁調(diào)坡頂升關(guān)鍵技術(shù)之頂升支撐技術(shù)探討
    南通通滬及東方大道連接線大跨度連續(xù)箱梁調(diào)坡頂升技術(shù)
    流域分布式坡長(zhǎng)不確定性的初步分析

    狂野欧美激情性xxxx| 美女被艹到高潮喷水动态| 国产亚洲精品综合一区在线观看| 日韩免费av在线播放| 一区福利在线观看| 亚洲七黄色美女视频| 午夜激情福利司机影院| 午夜久久久久精精品| 两性午夜刺激爽爽歪歪视频在线观看| 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 热99re8久久精品国产| 国产色婷婷99| 老汉色∧v一级毛片| 高清日韩中文字幕在线| 亚洲精品影视一区二区三区av| 国产成人欧美在线观看| 一本综合久久免费| 波野结衣二区三区在线 | 欧美成人一区二区免费高清观看| 亚洲avbb在线观看| 久久精品国产亚洲av涩爱 | 欧美性猛交╳xxx乱大交人| 欧美成狂野欧美在线观看| 天堂影院成人在线观看| 在线观看66精品国产| 欧美bdsm另类| 成人午夜高清在线视频| 欧美成人免费av一区二区三区| 国产亚洲精品综合一区在线观看| 最好的美女福利视频网| 亚洲18禁久久av| 午夜免费激情av| 国内久久婷婷六月综合欲色啪| 亚洲国产精品合色在线| 亚洲精品成人久久久久久| 超碰av人人做人人爽久久 | 日本黄大片高清| 中文字幕人成人乱码亚洲影| 女人高潮潮喷娇喘18禁视频| 免费看日本二区| 狂野欧美白嫩少妇大欣赏| 亚洲成a人片在线一区二区| 国产在视频线在精品| 免费人成在线观看视频色| 国产激情偷乱视频一区二区| 丁香六月欧美| 在线播放无遮挡| 亚洲国产精品合色在线| 日本黄色视频三级网站网址| 国产乱人视频| 一a级毛片在线观看| 亚洲熟妇中文字幕五十中出| 成年女人永久免费观看视频| 免费看美女性在线毛片视频| 亚洲av第一区精品v没综合| 午夜福利在线观看吧| 1024手机看黄色片| xxxwww97欧美| 日本在线视频免费播放| 免费观看精品视频网站| 亚洲中文日韩欧美视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲狠狠婷婷综合久久图片| 狂野欧美激情性xxxx| 欧美绝顶高潮抽搐喷水| 少妇人妻一区二区三区视频| 欧美成人a在线观看| 99热这里只有精品一区| 天天一区二区日本电影三级| 久久久久久久久大av| 51午夜福利影视在线观看| 午夜老司机福利剧场| 黄色片一级片一级黄色片| 成熟少妇高潮喷水视频| 午夜久久久久精精品| 18+在线观看网站| 老司机在亚洲福利影院| 日韩中文字幕欧美一区二区| 国产乱人视频| 在线免费观看的www视频| 男女床上黄色一级片免费看| 淫秽高清视频在线观看| 我要搜黄色片| 免费观看的影片在线观看| 热99re8久久精品国产| 午夜福利成人在线免费观看| 99热这里只有是精品50| 亚洲专区国产一区二区| 国产主播在线观看一区二区| 亚洲国产色片| 成人精品一区二区免费| 亚洲精品一区av在线观看| 国内毛片毛片毛片毛片毛片| 亚洲男人的天堂狠狠| 国产精品99久久久久久久久| av专区在线播放| 欧美zozozo另类| 51国产日韩欧美| 校园春色视频在线观看| 亚洲成av人片在线播放无| 国产99白浆流出| 国产主播在线观看一区二区| 综合色av麻豆| 国产亚洲欧美在线一区二区| 久久亚洲真实| 午夜精品久久久久久毛片777| 亚洲av电影不卡..在线观看| 亚洲在线自拍视频| 三级男女做爰猛烈吃奶摸视频| 久久久成人免费电影| 99riav亚洲国产免费| 夜夜躁狠狠躁天天躁| 久久精品人妻少妇| 香蕉丝袜av| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| www.999成人在线观看| 亚洲人成电影免费在线| 好看av亚洲va欧美ⅴa在| 欧美bdsm另类| 好男人在线观看高清免费视频| 成人欧美大片| 亚洲真实伦在线观看| 免费高清视频大片| 99riav亚洲国产免费| 免费在线观看日本一区| 女人高潮潮喷娇喘18禁视频| 一区二区三区免费毛片| 国产探花在线观看一区二区| 男人舔奶头视频| 制服丝袜大香蕉在线| 内地一区二区视频在线| 亚洲中文字幕一区二区三区有码在线看| 少妇的逼好多水| 丁香欧美五月| 搡老妇女老女人老熟妇| 国产久久久一区二区三区| 精品国产超薄肉色丝袜足j| 男女床上黄色一级片免费看| 中出人妻视频一区二区| 99视频精品全部免费 在线| 一二三四社区在线视频社区8| 亚洲人与动物交配视频| 伊人久久精品亚洲午夜| 国产91精品成人一区二区三区| 麻豆成人午夜福利视频| 搞女人的毛片| www.999成人在线观看| 精品无人区乱码1区二区| 日本熟妇午夜| 国产黄a三级三级三级人| 成年免费大片在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国产成年人精品一区二区| 热99re8久久精品国产| 欧美日韩瑟瑟在线播放| 久久精品国产清高在天天线| av片东京热男人的天堂| 精品福利观看| a级毛片a级免费在线| 特大巨黑吊av在线直播| 国产精品久久久久久久电影 | 久久久精品大字幕| 国产午夜精品久久久久久一区二区三区 | 中国美女看黄片| 国内精品久久久久精免费| 一进一出抽搐gif免费好疼| 精品久久久久久,| 美女被艹到高潮喷水动态| 国产老妇女一区| 悠悠久久av| 久久精品影院6| 禁无遮挡网站| 蜜桃久久精品国产亚洲av| 国产在线精品亚洲第一网站| 成人精品一区二区免费| 精品国产亚洲在线| 亚洲精品成人久久久久久| 床上黄色一级片| www.www免费av| 老司机福利观看| 亚洲一区二区三区不卡视频| 精品电影一区二区在线| 美女高潮喷水抽搐中文字幕| 亚洲精品久久国产高清桃花| 成年女人毛片免费观看观看9| 久久久久久久久中文| 亚洲一区二区三区色噜噜| 99久久精品热视频| 欧美成人a在线观看| 欧洲精品卡2卡3卡4卡5卡区| av天堂在线播放| 亚洲精华国产精华精| 精品国内亚洲2022精品成人| www.色视频.com| 黄色片一级片一级黄色片| 欧美日韩精品网址| 国产精品日韩av在线免费观看| 91麻豆av在线| 又粗又爽又猛毛片免费看| 丝袜美腿在线中文| 丁香六月欧美| 久久久久久久久久黄片| 亚洲一区高清亚洲精品| 在线天堂最新版资源| 最近视频中文字幕2019在线8| 欧洲精品卡2卡3卡4卡5卡区| 无遮挡黄片免费观看| 久久久久九九精品影院| 内射极品少妇av片p| 久久久久久久久久黄片| 蜜桃亚洲精品一区二区三区| 1000部很黄的大片| 亚洲国产精品成人综合色| 看免费av毛片| 一进一出好大好爽视频| 久久久国产成人精品二区| 动漫黄色视频在线观看| 日本a在线网址| 一级毛片女人18水好多| 男女下面进入的视频免费午夜| 精品熟女少妇八av免费久了| 午夜两性在线视频| 搡老妇女老女人老熟妇| 丁香六月欧美| 国产高清视频在线观看网站| 日本在线视频免费播放| 久久人人精品亚洲av| 日韩高清综合在线| 成人永久免费在线观看视频| 丰满乱子伦码专区| 一本精品99久久精品77| 每晚都被弄得嗷嗷叫到高潮| 搡女人真爽免费视频火全软件 | 精品久久久久久成人av| 91九色精品人成在线观看| 男人舔奶头视频| 九九热线精品视视频播放| 在线免费观看不下载黄p国产 | 最近最新免费中文字幕在线| 亚洲久久久久久中文字幕| 又黄又爽又免费观看的视频| 搡老岳熟女国产| 18禁美女被吸乳视频| 蜜桃亚洲精品一区二区三区| 成人高潮视频无遮挡免费网站| 精品国产美女av久久久久小说| 免费在线观看亚洲国产| 999久久久精品免费观看国产| 岛国在线免费视频观看| 小蜜桃在线观看免费完整版高清| 国产欧美日韩精品亚洲av| 亚洲av电影不卡..在线观看| 日韩欧美精品v在线| 一级作爱视频免费观看| 三级男女做爰猛烈吃奶摸视频| 国内毛片毛片毛片毛片毛片| 18+在线观看网站| 欧美黄色片欧美黄色片| 亚洲精品乱码久久久v下载方式 | 真实男女啪啪啪动态图| 国产精品香港三级国产av潘金莲| 欧美黄色淫秽网站| 一进一出抽搐gif免费好疼| 亚洲精品影视一区二区三区av| 在线观看舔阴道视频| 五月玫瑰六月丁香| 黄色丝袜av网址大全| 日韩成人在线观看一区二区三区| x7x7x7水蜜桃| 日韩免费av在线播放| 在线观看午夜福利视频| 成人国产综合亚洲| 日韩欧美免费精品| 麻豆一二三区av精品| 日韩高清综合在线| 婷婷亚洲欧美| 久久6这里有精品| 在线播放无遮挡| АⅤ资源中文在线天堂| 九色国产91popny在线| 午夜a级毛片| 桃色一区二区三区在线观看| av天堂在线播放| 一级作爱视频免费观看| 日本黄色视频三级网站网址| 一个人观看的视频www高清免费观看| 香蕉久久夜色| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 99久久九九国产精品国产免费| 国产色婷婷99| 村上凉子中文字幕在线| 18禁裸乳无遮挡免费网站照片| 免费一级毛片在线播放高清视频| 18禁黄网站禁片免费观看直播| 精品欧美国产一区二区三| 在线观看免费午夜福利视频| 成人国产综合亚洲| 久久亚洲精品不卡| 日韩有码中文字幕| 亚洲av中文字字幕乱码综合| 国产一区在线观看成人免费| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 国产成年人精品一区二区| 国产精品久久久久久久久免 | 色综合欧美亚洲国产小说| 少妇熟女aⅴ在线视频| av国产免费在线观看| 欧美成人a在线观看| 精品久久久久久久久久免费视频| 一区二区三区高清视频在线| 国产精品精品国产色婷婷| 99久国产av精品| 日本成人三级电影网站| 99精品欧美一区二区三区四区| 日韩有码中文字幕| 午夜日韩欧美国产| 波多野结衣高清无吗| 久久久久久久久大av| 国产亚洲av嫩草精品影院| 国产综合懂色| 国产精品影院久久| 女警被强在线播放| 久久久久久九九精品二区国产| 男插女下体视频免费在线播放| 亚洲av二区三区四区| 岛国在线免费视频观看| 欧美日韩瑟瑟在线播放| 成人午夜高清在线视频| 亚洲国产精品合色在线| 欧美zozozo另类| 伊人久久精品亚洲午夜| 亚洲av第一区精品v没综合| 国产熟女xx| 综合色av麻豆| 久久久久免费精品人妻一区二区| 免费观看人在逋| 国产熟女xx| 听说在线观看完整版免费高清| 色综合婷婷激情| 亚洲av日韩精品久久久久久密| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 精品久久久久久久久久免费视频| 天堂√8在线中文| 亚洲在线观看片| 国产精品,欧美在线| 国产av一区在线观看免费| 婷婷六月久久综合丁香| 法律面前人人平等表现在哪些方面| 久久精品国产自在天天线| 色综合欧美亚洲国产小说| www日本黄色视频网| 免费看日本二区| 成人三级黄色视频| 老司机午夜福利在线观看视频| 午夜影院日韩av| aaaaa片日本免费| 日韩精品青青久久久久久| 热99re8久久精品国产| 久久精品人妻少妇| 久久久久性生活片| 国产av在哪里看| 日本撒尿小便嘘嘘汇集6| 波多野结衣高清无吗| 精品无人区乱码1区二区| 神马国产精品三级电影在线观看| 热99re8久久精品国产| 国产三级中文精品| 国产97色在线日韩免费| 美女免费视频网站| 亚洲avbb在线观看| 国产精品久久久久久精品电影| 日本三级黄在线观看| 九九久久精品国产亚洲av麻豆| 国产久久久一区二区三区| 国产伦精品一区二区三区视频9 | 一a级毛片在线观看| 老司机福利观看| av欧美777| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 国产精品永久免费网站| 亚洲,欧美精品.| 99热这里只有是精品50| 亚洲精品粉嫩美女一区| 人人妻人人看人人澡| 亚洲av电影在线进入| 91在线观看av| 亚洲av第一区精品v没综合| 国内精品一区二区在线观看| 欧美xxxx黑人xx丫x性爽| 两性午夜刺激爽爽歪歪视频在线观看| 床上黄色一级片| 99国产综合亚洲精品| 最近在线观看免费完整版| 免费一级毛片在线播放高清视频| 美女大奶头视频| 一级作爱视频免费观看| 在线观看66精品国产| 免费大片18禁| 窝窝影院91人妻| 久久99热这里只有精品18| 首页视频小说图片口味搜索| 三级国产精品欧美在线观看| 免费人成视频x8x8入口观看| 色综合婷婷激情| 一二三四社区在线视频社区8| 亚洲久久久久久中文字幕| 人人妻人人澡欧美一区二区| 久久久国产成人免费| 欧美日韩综合久久久久久 | 日韩欧美在线二视频| 无遮挡黄片免费观看| 国产成人av教育| 亚洲乱码一区二区免费版| 99久久综合精品五月天人人| 精品国产美女av久久久久小说| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| netflix在线观看网站| 亚洲av一区综合| 亚洲国产精品合色在线| 久久久久九九精品影院| 国产乱人视频| 国产精品久久久久久久久免 | 成年女人永久免费观看视频| 无限看片的www在线观看| 亚洲国产精品成人综合色| 久久香蕉精品热| 麻豆成人午夜福利视频| 无人区码免费观看不卡| 亚洲成人免费电影在线观看| 熟女人妻精品中文字幕| 午夜福利18| 特大巨黑吊av在线直播| 色吧在线观看| 免费一级毛片在线播放高清视频| 一个人免费在线观看电影| 亚洲色图av天堂| 毛片女人毛片| 亚洲精品亚洲一区二区| 可以在线观看的亚洲视频| 国产精品98久久久久久宅男小说| 精品欧美国产一区二区三| 波多野结衣高清无吗| 国产精品久久久久久久久免 | 丰满乱子伦码专区| 欧美精品啪啪一区二区三区| 久久久久九九精品影院| 97碰自拍视频| 亚洲av不卡在线观看| 精品国产美女av久久久久小说| 怎么达到女性高潮| 在线观看美女被高潮喷水网站 | 一进一出抽搐gif免费好疼| 亚洲成a人片在线一区二区| 男人舔女人下体高潮全视频| 夜夜躁狠狠躁天天躁| 噜噜噜噜噜久久久久久91| 日韩精品中文字幕看吧| 亚洲五月婷婷丁香| 国产三级中文精品| 亚洲国产中文字幕在线视频| 亚洲熟妇中文字幕五十中出| 一本综合久久免费| 18禁裸乳无遮挡免费网站照片| ponron亚洲| 狂野欧美激情性xxxx| 亚洲激情在线av| 成人欧美大片| 成人午夜高清在线视频| 亚洲一区二区三区色噜噜| 亚洲成a人片在线一区二区| 国产主播在线观看一区二区| 免费人成在线观看视频色| 偷拍熟女少妇极品色| 欧美不卡视频在线免费观看| 麻豆国产av国片精品| 人人妻人人澡欧美一区二区| 丰满人妻熟妇乱又伦精品不卡| 变态另类丝袜制服| 少妇熟女aⅴ在线视频| 日本一二三区视频观看| 国产男靠女视频免费网站| 狂野欧美激情性xxxx| 日韩欧美精品v在线| 国产单亲对白刺激| 亚洲国产精品久久男人天堂| 成人特级黄色片久久久久久久| 久久久久九九精品影院| av福利片在线观看| 在线观看美女被高潮喷水网站 | 欧美日韩一级在线毛片| 在线观看av片永久免费下载| tocl精华| 亚洲一区二区三区色噜噜| 别揉我奶头~嗯~啊~动态视频| 日韩人妻高清精品专区| 国产美女午夜福利| 日韩人妻高清精品专区| 久久午夜亚洲精品久久| 欧美日韩综合久久久久久 | 内射极品少妇av片p| 日韩欧美免费精品| 天堂av国产一区二区熟女人妻| 亚洲精品美女久久久久99蜜臀| 久久久精品大字幕| 波多野结衣高清作品| 久久久久久久亚洲中文字幕 | 亚洲无线在线观看| 窝窝影院91人妻| 久久香蕉国产精品| 国产激情欧美一区二区| 有码 亚洲区| 日日摸夜夜添夜夜添小说| 国产av不卡久久| 人人妻人人看人人澡| 亚洲精华国产精华精| 亚洲 国产 在线| 精品99又大又爽又粗少妇毛片 | 丰满人妻一区二区三区视频av | 少妇的逼好多水| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 97人妻精品一区二区三区麻豆| 精品一区二区三区av网在线观看| 一夜夜www| 午夜福利在线观看免费完整高清在 | 少妇的丰满在线观看| 日韩欧美三级三区| 久久精品国产亚洲av涩爱 | 久久久久久久亚洲中文字幕 | 亚洲人成网站在线播| 国产成人福利小说| h日本视频在线播放| 非洲黑人性xxxx精品又粗又长| www.色视频.com| 人人妻人人看人人澡| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 久久久国产成人精品二区| 国产高清视频在线观看网站| 日本熟妇午夜| 欧美日韩福利视频一区二区| 婷婷六月久久综合丁香| 91在线观看av| 久久九九热精品免费| 五月伊人婷婷丁香| 国产伦精品一区二区三区四那| 天天躁日日操中文字幕| av天堂在线播放| xxx96com| 俄罗斯特黄特色一大片| 亚洲一区二区三区色噜噜| 色av中文字幕| 欧美性猛交黑人性爽| 在线观看免费视频日本深夜| 成人鲁丝片一二三区免费| 一二三四社区在线视频社区8| 色老头精品视频在线观看| 一级作爱视频免费观看| 国产伦精品一区二区三区视频9 | 日韩人妻高清精品专区| 搡老妇女老女人老熟妇| 精品99又大又爽又粗少妇毛片 | 欧美一区二区亚洲| 亚洲av熟女| 免费无遮挡裸体视频| 国模一区二区三区四区视频| 动漫黄色视频在线观看| 欧美色视频一区免费| 在线观看av片永久免费下载| 两人在一起打扑克的视频| 91av网一区二区| 高清在线国产一区| 亚洲18禁久久av| 日日摸夜夜添夜夜添小说| 天天添夜夜摸| 亚洲最大成人中文| 夜夜爽天天搞| 午夜福利在线观看免费完整高清在 | 国产高清视频在线播放一区| 美女 人体艺术 gogo| 丁香欧美五月| 精品人妻一区二区三区麻豆 | 听说在线观看完整版免费高清| 999久久久精品免费观看国产| 熟女少妇亚洲综合色aaa.| 欧美日韩中文字幕国产精品一区二区三区| 国产精品亚洲一级av第二区| 蜜桃亚洲精品一区二区三区| 久久精品综合一区二区三区| 日韩中文字幕欧美一区二区| 色av中文字幕| 国产精品自产拍在线观看55亚洲| 欧美日韩黄片免| 日韩欧美在线二视频| 深爱激情五月婷婷| 精品一区二区三区人妻视频| 国内久久婷婷六月综合欲色啪| 欧美日韩乱码在线| 好男人电影高清在线观看| 国产高清三级在线| 成人午夜高清在线视频| 99国产极品粉嫩在线观看| 人妻丰满熟妇av一区二区三区| 老鸭窝网址在线观看| 亚洲欧美日韩无卡精品| 窝窝影院91人妻| 日本一二三区视频观看| 日本 av在线|