• <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)不確定性的初步分析

    日本黄色视频三级网站网址 | 中亚洲国语对白在线视频| 欧美亚洲 丝袜 人妻 在线| 在线av久久热| 欧美乱码精品一区二区三区| videos熟女内射| 亚洲男人天堂网一区| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 国精品久久久久久国模美| 成年人午夜在线观看视频| 国产日韩欧美亚洲二区| 国产精品av久久久久免费| 免费在线观看黄色视频的| av一本久久久久| 一边摸一边抽搐一进一出视频| 国产深夜福利视频在线观看| 亚洲性夜色夜夜综合| 国产色视频综合| 国产人伦9x9x在线观看| 中文亚洲av片在线观看爽 | 人成视频在线观看免费观看| 成年女人毛片免费观看观看9 | 国产精品99久久99久久久不卡| 国产在线观看jvid| 亚洲国产精品一区二区三区在线| 国产成人欧美在线观看 | 在线观看人妻少妇| 男女免费视频国产| 法律面前人人平等表现在哪些方面| 老司机深夜福利视频在线观看| 蜜桃国产av成人99| 色在线成人网| 另类亚洲欧美激情| 日本精品一区二区三区蜜桃| a级片在线免费高清观看视频| 国产在线免费精品| 男女之事视频高清在线观看| 欧美亚洲 丝袜 人妻 在线| 无人区码免费观看不卡 | 国产精品久久久av美女十八| 免费在线观看完整版高清| 视频区欧美日本亚洲| cao死你这个sao货| 免费观看a级毛片全部| 亚洲第一欧美日韩一区二区三区 | 亚洲全国av大片| 亚洲精品美女久久av网站| 欧美乱妇无乱码| 一本色道久久久久久精品综合| 国产精品一区二区免费欧美| 中文字幕高清在线视频| 精品福利观看| 18禁观看日本| avwww免费| 国产午夜精品久久久久久| 成人黄色视频免费在线看| 两个人看的免费小视频| 欧美日韩亚洲高清精品| 午夜福利在线免费观看网站| 中文字幕色久视频| 精品福利观看| 国产精品二区激情视频| 国产欧美日韩综合在线一区二区| 黑人操中国人逼视频| 国产深夜福利视频在线观看| 99精国产麻豆久久婷婷| 老司机福利观看| 国产精品熟女久久久久浪| 国产成人系列免费观看| 夜夜爽天天搞| 黄色a级毛片大全视频| 妹子高潮喷水视频| av有码第一页| 成人国语在线视频| 欧美日韩精品网址| 亚洲久久久国产精品| 久久九九热精品免费| 欧美久久黑人一区二区| 午夜视频精品福利| 亚洲九九香蕉| 亚洲国产中文字幕在线视频| 热re99久久精品国产66热6| 深夜精品福利| 久久久精品国产亚洲av高清涩受| 欧美 亚洲 国产 日韩一| 国产成人系列免费观看| 汤姆久久久久久久影院中文字幕| 国精品久久久久久国模美| www.自偷自拍.com| 亚洲国产精品一区二区三区在线| 久久人妻av系列| 999久久久精品免费观看国产| 国产xxxxx性猛交| 新久久久久国产一级毛片| 亚洲中文日韩欧美视频| 国产精品久久久久久精品电影小说| 国产老妇伦熟女老妇高清| 国产真人三级小视频在线观看| 最近最新中文字幕大全免费视频| 91大片在线观看| 成人特级黄色片久久久久久久 | 精品国产乱码久久久久久小说| 中文字幕高清在线视频| 国产不卡一卡二| 久久精品成人免费网站| 午夜免费鲁丝| 男女午夜视频在线观看| 一级片免费观看大全| 国产片内射在线| 999久久久精品免费观看国产| 老汉色av国产亚洲站长工具| 黄色片一级片一级黄色片| 亚洲成人手机| 在线观看人妻少妇| 最近最新免费中文字幕在线| 精品一区二区三区av网在线观看 | 午夜福利在线观看吧| av天堂久久9| 12—13女人毛片做爰片一| 女性生殖器流出的白浆| 久久久国产一区二区| 少妇 在线观看| 久久av网站| 麻豆av在线久日| 最近最新中文字幕大全电影3 | 在线观看一区二区三区激情| 久久精品亚洲av国产电影网| 9191精品国产免费久久| 80岁老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 中文字幕av电影在线播放| 国产在线一区二区三区精| 免费女性裸体啪啪无遮挡网站| 午夜成年电影在线免费观看| 午夜视频精品福利| 欧美在线一区亚洲| 国产99久久九九免费精品| 国产精品影院久久| 午夜福利在线观看吧| 精品少妇久久久久久888优播| 欧美+亚洲+日韩+国产| 亚洲av第一区精品v没综合| 久久久精品94久久精品| 亚洲精品在线观看二区| 一边摸一边抽搐一进一小说 | 激情视频va一区二区三区| 午夜福利欧美成人| 正在播放国产对白刺激| 久久人妻av系列| 丁香六月欧美| 黑人巨大精品欧美一区二区mp4| 丁香欧美五月| 国产午夜精品久久久久久| 免费看十八禁软件| 国产成人啪精品午夜网站| 王馨瑶露胸无遮挡在线观看| 在线 av 中文字幕| 少妇 在线观看| 免费看a级黄色片| 少妇精品久久久久久久| 91成年电影在线观看| 久久久国产成人免费| www.熟女人妻精品国产| 激情在线观看视频在线高清 | 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 日韩大片免费观看网站| 中文字幕高清在线视频| 精品亚洲成国产av| 亚洲精品中文字幕一二三四区 | 在线观看一区二区三区激情| 久久亚洲真实| 多毛熟女@视频| 激情在线观看视频在线高清 | 免费高清在线观看日韩| 老司机午夜十八禁免费视频| 免费一级毛片在线播放高清视频 | 飞空精品影院首页| 中文字幕制服av| 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 久久国产精品人妻蜜桃| 每晚都被弄得嗷嗷叫到高潮| 宅男免费午夜| 国产精品自产拍在线观看55亚洲 | 脱女人内裤的视频| 日韩免费av在线播放| 19禁男女啪啪无遮挡网站| 日韩有码中文字幕| 日本av手机在线免费观看| 男人舔女人的私密视频| 黄色视频,在线免费观看| 国产精品成人在线| 免费av中文字幕在线| 美国免费a级毛片| 成年动漫av网址| 久久九九热精品免费| 亚洲欧洲精品一区二区精品久久久| 在线观看免费高清a一片| 激情视频va一区二区三区| 亚洲第一青青草原| 亚洲天堂av无毛| 一区二区三区国产精品乱码| 亚洲三区欧美一区| 欧美午夜高清在线| 最新的欧美精品一区二区| 伦理电影免费视频| 少妇 在线观看| 国产精品免费一区二区三区在线 | 国产色视频综合| 色播在线永久视频| 丰满饥渴人妻一区二区三| 亚洲视频免费观看视频| 成人黄色视频免费在线看| 精品卡一卡二卡四卡免费| 午夜两性在线视频| 国产在线免费精品| 91麻豆av在线| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 后天国语完整版免费观看| 91麻豆av在线| 亚洲av片天天在线观看| 91字幕亚洲| 一个人免费在线观看的高清视频| 欧美日韩一级在线毛片| 一级毛片精品| 日韩免费av在线播放| 久久久国产精品麻豆| 精品国产一区二区久久| 一区二区三区国产精品乱码| videosex国产| 女警被强在线播放| 女性被躁到高潮视频| 国产伦理片在线播放av一区| 亚洲人成伊人成综合网2020| 在线十欧美十亚洲十日本专区| 男女午夜视频在线观看| 大陆偷拍与自拍| 国产一区二区在线观看av| 久久国产精品影院| 日韩欧美免费精品| 在线观看免费高清a一片| 国产精品国产av在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲成人免费电影在线观看| 天天操日日干夜夜撸| 日日爽夜夜爽网站| 免费在线观看完整版高清| 中文字幕色久视频| 男女午夜视频在线观看| 99精品久久久久人妻精品| 老司机在亚洲福利影院| 久久精品国产亚洲av高清一级| 黄色毛片三级朝国网站| 考比视频在线观看| 一级片'在线观看视频| 老熟妇仑乱视频hdxx| 亚洲精华国产精华精| 亚洲av成人一区二区三| 国产一区二区三区综合在线观看| 最近最新中文字幕大全电影3 | 宅男免费午夜| 欧美国产精品一级二级三级| 亚洲国产看品久久| 亚洲精品乱久久久久久| 纯流量卡能插随身wifi吗| 亚洲精品自拍成人| 成年人午夜在线观看视频| 久久国产精品大桥未久av| 50天的宝宝边吃奶边哭怎么回事| 宅男免费午夜| 高清在线国产一区| av欧美777| 女警被强在线播放| 久久久精品区二区三区| 1024视频免费在线观看| 美女国产高潮福利片在线看| 久久精品国产亚洲av高清一级| 精品高清国产在线一区| 天堂动漫精品| 一级毛片电影观看| 两个人免费观看高清视频| 99精品久久久久人妻精品| 免费在线观看视频国产中文字幕亚洲| 欧美+亚洲+日韩+国产| 99国产精品99久久久久| 男女午夜视频在线观看| 男女之事视频高清在线观看| 国产成人啪精品午夜网站| 亚洲熟女毛片儿| 18禁国产床啪视频网站| 久久香蕉激情| 久久国产亚洲av麻豆专区| 丝瓜视频免费看黄片| 狠狠婷婷综合久久久久久88av| 午夜老司机福利片| 在线观看舔阴道视频| 纯流量卡能插随身wifi吗| 少妇被粗大的猛进出69影院| 欧美精品一区二区免费开放| 久久精品国产99精品国产亚洲性色 | 亚洲午夜理论影院| 美女国产高潮福利片在线看| 久久精品国产a三级三级三级| 99国产综合亚洲精品| av网站在线播放免费| 一边摸一边抽搐一进一出视频| 另类亚洲欧美激情| 国产精品偷伦视频观看了| 亚洲精品乱久久久久久| 欧美精品亚洲一区二区| 人妻 亚洲 视频| 一区二区三区精品91| 99国产精品99久久久久| 亚洲第一av免费看| 久久精品成人免费网站| 国产成人精品久久二区二区免费| 亚洲第一av免费看| 满18在线观看网站| 久久99热这里只频精品6学生| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 亚洲精品国产色婷婷电影| 精品少妇一区二区三区视频日本电影| 国产免费现黄频在线看| 亚洲精品乱久久久久久| 麻豆成人av在线观看| 欧美日本中文国产一区发布| 欧美日韩一级在线毛片| 香蕉久久夜色| 黑人巨大精品欧美一区二区蜜桃| 人成视频在线观看免费观看| 久久午夜综合久久蜜桃| 建设人人有责人人尽责人人享有的| 十八禁网站网址无遮挡| bbb黄色大片| 咕卡用的链子| 国产不卡av网站在线观看| 免费人妻精品一区二区三区视频| 啦啦啦在线免费观看视频4| 女人高潮潮喷娇喘18禁视频| 亚洲七黄色美女视频| xxxhd国产人妻xxx| 国产淫语在线视频| 这个男人来自地球电影免费观看| 久久午夜综合久久蜜桃| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 一进一出抽搐动态| 热99国产精品久久久久久7| 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| 91精品国产国语对白视频| 国产精品影院久久| 亚洲av第一区精品v没综合| 精品一区二区三区av网在线观看 | av视频免费观看在线观看| 热99re8久久精品国产| 搡老熟女国产l中国老女人| 亚洲视频免费观看视频| 国产极品粉嫩免费观看在线| 欧美激情久久久久久爽电影 | 一二三四社区在线视频社区8| 国产不卡av网站在线观看| 青草久久国产| 久久精品国产综合久久久| 精品一区二区三卡| 午夜精品久久久久久毛片777| 手机成人av网站| 国产亚洲欧美在线一区二区| 国产精品欧美亚洲77777| 亚洲专区中文字幕在线| 国产精品 国内视频| 亚洲av片天天在线观看| 18禁国产床啪视频网站| 一本—道久久a久久精品蜜桃钙片| 动漫黄色视频在线观看| 亚洲精品中文字幕一二三四区 | 亚洲欧美色中文字幕在线| 国产成人欧美在线观看 | 亚洲一码二码三码区别大吗| 国产精品免费视频内射| 成人手机av| 久久婷婷成人综合色麻豆| 国产欧美日韩一区二区精品| 丝瓜视频免费看黄片| 日韩一卡2卡3卡4卡2021年| 亚洲伊人色综图| 国产97色在线日韩免费| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 国产91精品成人一区二区三区 | 免费高清在线观看日韩| 欧美 亚洲 国产 日韩一| 国产福利在线免费观看视频| 精品少妇一区二区三区视频日本电影| 欧美乱妇无乱码| 成人av一区二区三区在线看| 一区二区av电影网| 国产野战对白在线观看| 亚洲综合色网址| 一区二区日韩欧美中文字幕| 黑人欧美特级aaaaaa片| 99在线人妻在线中文字幕 | 国产男靠女视频免费网站| 首页视频小说图片口味搜索| 久久精品熟女亚洲av麻豆精品| 99热网站在线观看| 亚洲精品中文字幕在线视频| 免费高清在线观看日韩| 视频区欧美日本亚洲| 国产97色在线日韩免费| 天堂俺去俺来也www色官网| 日本a在线网址| 精品卡一卡二卡四卡免费| 国产男女超爽视频在线观看| av天堂久久9| 80岁老熟妇乱子伦牲交| 精品国产一区二区久久| 亚洲成人免费电影在线观看| 日韩视频一区二区在线观看| 看免费av毛片| 欧美老熟妇乱子伦牲交| 一级毛片女人18水好多| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情久久久久久爽电影 | 久久久久久久国产电影| 9191精品国产免费久久| 91av网站免费观看| 97人妻天天添夜夜摸| 亚洲九九香蕉| 少妇猛男粗大的猛烈进出视频| 精品卡一卡二卡四卡免费| 男女午夜视频在线观看| 精品国产乱码久久久久久小说| 下体分泌物呈黄色| 香蕉丝袜av| 啦啦啦视频在线资源免费观看| 考比视频在线观看| 日韩大片免费观看网站| 飞空精品影院首页| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 精品久久久久久久毛片微露脸| 丰满迷人的少妇在线观看| 精品一区二区三卡| 性少妇av在线| 飞空精品影院首页| 汤姆久久久久久久影院中文字幕| 美女国产高潮福利片在线看| 久久中文字幕一级| 看免费av毛片| 欧美黄色片欧美黄色片| www.精华液| 天天影视国产精品| 国产片内射在线| 亚洲熟妇熟女久久| 国产亚洲精品一区二区www | 亚洲伊人久久精品综合| 在线观看免费视频网站a站| 中文字幕最新亚洲高清| 女性生殖器流出的白浆| 亚洲一区中文字幕在线| 欧美日韩视频精品一区| 午夜免费成人在线视频| 黑人猛操日本美女一级片| 一二三四社区在线视频社区8| 亚洲av成人一区二区三| 国产精品 国内视频| 久久婷婷成人综合色麻豆| 99热网站在线观看| 久久毛片免费看一区二区三区| 精品国内亚洲2022精品成人 | 十八禁人妻一区二区| 午夜精品国产一区二区电影| 久久精品国产99精品国产亚洲性色 | av网站免费在线观看视频| 精品国产超薄肉色丝袜足j| 亚洲国产av新网站| 精品人妻在线不人妻| 黄色a级毛片大全视频| 久久国产精品大桥未久av| 欧美日韩av久久| 国产欧美日韩一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕最新亚洲高清| 亚洲人成77777在线视频| 下体分泌物呈黄色| av福利片在线| 亚洲国产看品久久| av国产精品久久久久影院| 黄频高清免费视频| 精品一区二区三区视频在线观看免费 | 精品一区二区三卡| 日本av手机在线免费观看| 久久人妻福利社区极品人妻图片| 一边摸一边做爽爽视频免费| 老司机午夜十八禁免费视频| 少妇被粗大的猛进出69影院| 亚洲成人手机| 日韩人妻精品一区2区三区| 久久婷婷成人综合色麻豆| 18禁裸乳无遮挡动漫免费视频| 国产一卡二卡三卡精品| 午夜91福利影院| 麻豆av在线久日| 久久人人97超碰香蕉20202| 2018国产大陆天天弄谢| 19禁男女啪啪无遮挡网站| 欧美国产精品va在线观看不卡| 青草久久国产| 国产精品麻豆人妻色哟哟久久| 最近最新中文字幕大全电影3 | 99riav亚洲国产免费| 国产深夜福利视频在线观看| 免费av中文字幕在线| 国产亚洲欧美在线一区二区| 午夜激情av网站| 国产精品美女特级片免费视频播放器 | 久久久久久亚洲精品国产蜜桃av| 宅男免费午夜| 香蕉久久夜色| 国产色视频综合| 无遮挡黄片免费观看| 久久人妻av系列| 俄罗斯特黄特色一大片| 天堂俺去俺来也www色官网| 桃花免费在线播放| 午夜福利,免费看| 国产亚洲精品一区二区www | av国产精品久久久久影院| 亚洲七黄色美女视频| 国产成人欧美| 美女视频免费永久观看网站| 日韩中文字幕视频在线看片| 一二三四社区在线视频社区8| 亚洲va日本ⅴa欧美va伊人久久| 欧美精品一区二区免费开放| 在线观看一区二区三区激情| 精品国内亚洲2022精品成人 | 成人黄色视频免费在线看| 免费少妇av软件| av国产精品久久久久影院| cao死你这个sao货| 久久久久久久大尺度免费视频| 日韩大片免费观看网站| 一二三四社区在线视频社区8| 欧美精品啪啪一区二区三区| 免费观看a级毛片全部| 久久狼人影院| 在线观看免费高清a一片| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区蜜桃| 正在播放国产对白刺激| 欧美日韩亚洲综合一区二区三区_| 丰满饥渴人妻一区二区三| 亚洲第一欧美日韩一区二区三区 | 欧美大码av| 香蕉丝袜av| 欧美亚洲日本最大视频资源| 亚洲精品国产一区二区精华液| 亚洲av日韩精品久久久久久密| 亚洲欧美日韩另类电影网站| 成人国产av品久久久| 亚洲情色 制服丝袜| 久久中文字幕一级| 999精品在线视频| av网站在线播放免费| 日韩熟女老妇一区二区性免费视频| 亚洲精品一卡2卡三卡4卡5卡| 日本vs欧美在线观看视频| 午夜福利,免费看| 女性被躁到高潮视频| 黄色片一级片一级黄色片| 蜜桃在线观看..| 老司机深夜福利视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品香港三级国产av潘金莲| 亚洲国产欧美网| 成人手机av| 国产淫语在线视频| 9191精品国产免费久久| 日本av免费视频播放| 国产精品香港三级国产av潘金莲| 亚洲成人手机| 国产精品成人在线| 国产淫语在线视频| 久久精品国产综合久久久| 国产av又大| avwww免费| 丰满少妇做爰视频| 一区在线观看完整版| 久久精品国产亚洲av香蕉五月 | 午夜福利视频在线观看免费| 久久青草综合色| 亚洲精品自拍成人| 国产精品久久久av美女十八| 真人做人爱边吃奶动态| 欧美日韩成人在线一区二区| 美女午夜性视频免费| 99久久99久久久精品蜜桃| 久久九九热精品免费| 亚洲精品久久成人aⅴ小说| 一区二区三区乱码不卡18| 亚洲一卡2卡3卡4卡5卡精品中文| 另类亚洲欧美激情| 久久精品亚洲熟妇少妇任你| 午夜老司机福利片| 亚洲成av片中文字幕在线观看| 精品福利观看| 一本—道久久a久久精品蜜桃钙片| 国产精品九九99| 日本黄色视频三级网站网址 |