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

    福建地區(qū)馬尾松生物量轉(zhuǎn)換和擴(kuò)展因子的影響因素

    2017-10-13 12:32:58歐強(qiáng)新李???/span>
    生態(tài)學(xué)報(bào) 2017年17期
    關(guān)鍵詞:坡向馬尾松林分

    歐強(qiáng)新,李???*,楊 英

    1 中國(guó)林業(yè)科學(xué)研究院資源信息研究所, 北京 100091 2 國(guó)家林業(yè)局調(diào)查規(guī)劃設(shè)計(jì)院, 北京 100714

    福建地區(qū)馬尾松生物量轉(zhuǎn)換和擴(kuò)展因子的影響因素

    歐強(qiáng)新1,李???,*,楊 英2

    1 中國(guó)林業(yè)科學(xué)研究院資源信息研究所, 北京 100091 2 國(guó)家林業(yè)局調(diào)查規(guī)劃設(shè)計(jì)院, 北京 100714

    基于第8次國(guó)家森林資源清查福建省331塊馬尾松的固定樣地調(diào)查數(shù)據(jù),利用增強(qiáng)回歸樹(shù)法(BRT)研究地上生物量轉(zhuǎn)換和擴(kuò)展因子(BCEF)和地下BCEF的影響因素。研究結(jié)果表明: 林分特征因子和地形因子是影響地上BCEF以及地下BCEF的主導(dǎo)因素,二者對(duì)地上BCEF影響的相對(duì)貢獻(xiàn)率之和為87.20%、地下BCEF為86.59%。其中,齡組和坡向分別是林分特征因子和地形因子中影響地上BCEF的最大因素(41.13%和14.52%)。地上BCEF隨齡組的增大而逐漸減??;在東南坡最大、西坡最小。此外,齡組和坡向分別是林分特征因子和地形因子中影響地下BCEF的最大因素(41.54%和15.16%)。地下BCEF隨齡組的增大而逐漸減增大;在東南坡最小、西坡最大。土壤因子對(duì)地上BCEF以及地下BCEF的影響都較小(12.80%和13.41%),腐殖層厚度是土壤因子中影響地上BCEF以及地下BCEF的最大因素(9.02%和9.13%)。在所有的影響因素中,齡組對(duì)地上BCEF以及地下BCEF的影響均最大,依據(jù)齡組計(jì)算相應(yīng)的BCEF或者建立林齡普適的BCEF模型,可以有效地提高生物量的估算精度。

    馬尾松;生物量轉(zhuǎn)換和擴(kuò)展因子(BCEF);國(guó)家森林資源清查;增強(qiáng)回歸樹(shù)(BRT)

    Abstract: In this study, the boosted regression trees (BRT) method was used to analyze the main factors that influenced the above- and underground biomass conversion and expansion factor (BCEF) of Masson pine. Data from the permanent sample plot of the 8thChinese National Forest Inventory for Fujian Province was used. The results showed that the cumulative relative influences of stand characteristic and topographic factors on above- and underground BCEF were 87.20 and 86.59%, respectively. Stand characteristic and topographic factors were the most important factors for both above- and underground BCEF. Among the stand characteristic factors, the age group was dominant, contributing to 41.13% of the aboveground BCEF differences that decreased with an increase in the age group, and 41.54% of the underground BCEF differences that increased with an increase in the age group. Among the topographic factors, slope was the key variable, contributing to 14.52% of aboveground BCEF differences, which were the highest at the southeast slope and the lowest at the west slope, and 15.16% of the underground BCEF differences, which were the lowest at the southeast slope and the highest at the west slope. In addition, soil factors had a lower effect on both above and underground BCEF. The relative influences of soil factors on above- and underground BCEF were 12.80 and 13.41%, respectively. Among the soil factors, the thickness of the humus layer was the key variable, contributing to 9.02% of the aboveground BCEF differences and 9.13% of the underground BCEF differences. Among all influencing factors, age group had the highest effect on both above- and underground BCEF. Consequently, the ways, such as calculating BCEF based on age group or building universal BCEF models of all age, could achieve a better precision of biomass.

    KeyWords: Masson pine; biomass conversion and expansion factor (BCEF); national forest inventory; boosted regression trees (BRT)

    森林生物量是評(píng)估森林碳收支的重要參數(shù)[1- 2],準(zhǔn)確的生物量估算對(duì)量化碳固定率、評(píng)估氣候變化的潛在影響等內(nèi)容至關(guān)重要[3- 4]。生物量轉(zhuǎn)換和擴(kuò)展因子(BCEF)是估算森林生物量普遍使用的估算參數(shù)[5],BCEF可以直接將蓄積轉(zhuǎn)換并擴(kuò)展為生物量[6-13]。對(duì)林分BCEF的影響因素進(jìn)行研究,確定影響B(tài)CEF的主要因素,可提高林分生物量估算的準(zhǔn)確性[14]。已有的研究表明,BCEF并非是一種固定不變的常數(shù)[7,15-17],而與立地、林齡等因子密切相關(guān)[15,18-21]。如Lehtonen等[18]發(fā)現(xiàn)歐洲赤松(Pinussylvestris)和歐洲云杉(Piceaabies)總BCEF均隨著林分年齡的增加而變小,而地下BCEF均隨著林分年齡的增加而變大;Soares and Tomé[15]發(fā)現(xiàn)藍(lán)桉樹(shù)(Eucalyptusglobulus)總BCEF隨立地指數(shù)的增加而減?。籊onzález-García等[20]發(fā)現(xiàn)亮果桉(Eucalyptusnitens)地上BCEF與林分年齡呈現(xiàn)出十分顯著的負(fù)相關(guān)關(guān)系?,F(xiàn)有學(xué)者主要利用回歸模型、方差分析以及相關(guān)性分析等統(tǒng)計(jì)學(xué)方法,對(duì)BCEF的影響因素進(jìn)行了研究[5, 18- 20,22]。如左舒翟等[5]利用回歸模型方法分析了林分特征因子對(duì)杉木(Cunninghamialanceolata)地上BCEF的影響;羅云建等[19]利用回歸模型、相關(guān)性分析方法研究了林分特征因子對(duì)落葉松(Larixgmelinii)地上BCEF的影響;吳小山[22]利用方差分析方法研究了地形、土壤和人為因子對(duì)楊樹(shù)(Populus)地上BCEF的影響。然而,關(guān)于林分特征因子與環(huán)境因子如何協(xié)同影響B(tài)CEF的以及對(duì)BCEF影響的相對(duì)貢獻(xiàn)率是怎樣的,目前尚未見(jiàn)報(bào)道。

    馬尾松(Pinusmassoniana)在造林、用材等方面都有著較為廣泛的應(yīng)用[23]。由第八次國(guó)家森林資源清查可知,馬尾松總面積和總蓄積高達(dá)1.0×107hm2和5.9×108m3,分別占全國(guó)喬木林總面積和總蓄積的6.1%和4.0%。福建省森林覆蓋率為65.95%,其中福建省馬尾松面積占福建省喬木林總面積的比例為13.79%,福建省馬尾松蓄積占福建省喬木林總蓄積的比例為11.88%。在我國(guó),森林資源調(diào)查的主要內(nèi)容由立地與土壤、森林特征、森林功能、土地利用與覆蓋以及其他因素等5大方面,其中,林分特征因子、土壤因子以及地形因子這3類因子不僅與林木的生長(zhǎng)息息相關(guān),而且相對(duì)容易獲取。增強(qiáng)回歸樹(shù)(BRT)是一種利用隨機(jī)抽樣和自學(xué)習(xí)方法而建立的多重回歸樹(shù)[24],BRT方法提高了計(jì)算結(jié)果的穩(wěn)定性和精度,可獲得自變量對(duì)因變量影響的相對(duì)貢獻(xiàn)率,已經(jīng)在分析土地利用的影響因素[25]、區(qū)域生物量的影響因素[26]以及火災(zāi)的影響因素[27]等研究中得到成功應(yīng)用。本文以第8次國(guó)家森林資源清查福建省優(yōu)勢(shì)樹(shù)種為馬尾松的固定樣地?cái)?shù)據(jù)作為研究材料,將BRT方法和森林資源調(diào)查數(shù)據(jù)相結(jié)合,分別確定地上BCEF以及地下BCEF的影響因素及其相對(duì)貢獻(xiàn)率,以便提高林分生物量估算的準(zhǔn)確性。

    1 研究區(qū)概況與數(shù)據(jù)

    福建省位于23°33′—28°20′N,115°50′—120°40′E之間,屬亞熱帶海洋性季風(fēng)氣候,以丘陵和山地為主,廣泛分布著黃壤、紅壤以及磚紅壤,最高海拔達(dá)2158米。 森林資源主要分布于龍巖、三明以及南平3個(gè)區(qū)市。馬尾松、杉木以及桉樹(shù)(Eucalyptus)是福建省分布最為廣泛的3個(gè)喬木樹(shù)種,其面積占福建喬木林總面積的比例分別為13.79%、22.36%以及4.47%;天然林中,馬尾松、針葉混交林以及闊葉混交林占天然林面積比重相對(duì)較大,分別為9.44%、12.81%和57.27%;人工林中,桉樹(shù)、馬尾松林以及杉木占人工林面積比重相對(duì)較大,分別為9.77%、18.95%和40.35%。

    福建省采用系統(tǒng)抽樣方法、按4km×6km公里網(wǎng)格布設(shè)方形固定樣地,共5059塊。樣地面積為0.0667hm2。從5059塊樣地中篩選出優(yōu)勢(shì)樹(shù)種為馬尾松的樣地331塊,然后對(duì)篩選的樣地做進(jìn)一步的處理:去除其它伴生樹(shù)種以及檢尺類型分別為多測(cè)木、枯倒木、采伐木以及枯立木的樣木,只保留馬尾松活立木。樣地的調(diào)查因子包括:起源、齡組、地貌、海拔、坡向、坡位、坡度、土壤種類、土壤厚度和腐殖層厚度等。上述各因子的確定與劃分標(biāo)準(zhǔn)參照國(guó)家森林資源清查技術(shù)規(guī)定。

    2 方法

    2.1 單木生物量和材積計(jì)算

    本研究采用生物量模型和材積模型計(jì)算相應(yīng)的生物量和材積,模型形式以及參數(shù)選用中華人民共和國(guó)林業(yè)行業(yè)標(biāo)準(zhǔn)《LY/T 2263—2014:立木生物量模型及碳計(jì)量參數(shù)——馬尾松》中總體一的相關(guān)模型公式,該標(biāo)準(zhǔn)中依據(jù)總體一所建立的相關(guān)模型適用于福建省。上述標(biāo)準(zhǔn)的建模數(shù)據(jù)在選取過(guò)程中,充分考慮了在立地條件、起源以及齡組等方面具有廣泛代表性的樣本,建模數(shù)據(jù)調(diào)查取樣的方法詳見(jiàn)曾偉生[28]一文。由于國(guó)家森林資源清查只對(duì)胸徑大于等于5cm的樣木進(jìn)行了每木檢尺,而未進(jìn)行每木測(cè)高[29- 31],故本研究選用以胸徑(胸徑≥5cm)為唯一自變量的一元模型計(jì)算各單株木的生物量和材積。表1為具體的模型形式。

    表1 生物量模型和材積模型的模型形式以及評(píng)價(jià)指標(biāo)

    D:胸徑(cm且≥5cm)

    2.2 樣地水平BCEF計(jì)算

    BCEF可以將蓄積轉(zhuǎn)換并擴(kuò)展為生物量,單位為t/m3。在本研究中,依據(jù)下述公式計(jì)算樣地水平各項(xiàng)BCEF[18]:

    (1)

    式中,i表示林木第i個(gè)分項(xiàng)(如:地上部分、地下部分);k表示第k個(gè)樣地;j表示第j株樹(shù);Nk表示第k個(gè)樣地中樹(shù)木的總株樹(shù);BCEFik表示第k個(gè)樣地、分項(xiàng)i的樣地水平生物量轉(zhuǎn)換和擴(kuò)展因子;Bij表示第k個(gè)樣地、第j株樹(shù)、分項(xiàng)i的生物量;Vj表示第k個(gè)樣地、第j株樹(shù)的材積。

    各項(xiàng)BCEF的計(jì)算結(jié)果以及部分土壤、地形因子的統(tǒng)計(jì)見(jiàn)表2;部分土壤、地形和林分特征因子各級(jí)(類)別樣地?cái)?shù)的統(tǒng)計(jì)見(jiàn)表3。

    表2 各項(xiàng)BCEF以及部分土壤、地形因子的統(tǒng)計(jì)

    表3部分土壤、地形、林分特征因子樣地?cái)?shù)的統(tǒng)計(jì)

    Table3Statisticsofplotnumberforsomesoilfactors,topographyfactorsandstandcharacteristicfactors

    土壤因子Soilfactor級(jí)(類)別Grade(Sort)樣地?cái)?shù)Numberofplot土壤種類Soiltype磚紅壤Latosol8紅壤Redsoil276黃壤Yellowsoil47坡位Slopeposition山脊Mountainridge21上坡Upslope139中坡Mesoslope97下坡Downslope74坡向Slopeaspect北坡Northslope26東北坡Northeastslope42東坡Eastslope38東南坡Southeastslope70南坡Southslope39西南坡Southwestslope54西坡Westslope31西北坡Northwestslope31地貌Landform中山Middlemountain49低山Lowmountain127丘陵Hilly155齡組Agegroup幼齡Young42中齡Middleage136近熟NearMature105成熟Mature38過(guò)熟Overmature10起源Origin天然Natural120人工Plantation211

    2.3 BRT法分析主要調(diào)查因子對(duì)BCEF的影響

    BRT法是一種基于分類和回歸樹(shù)的、高級(jí)形式的機(jī)器學(xué)習(xí)方法[32]。BRT法利用訓(xùn)練數(shù)據(jù)分析各個(gè)自變量對(duì)因變量的影響、利用測(cè)試數(shù)據(jù)檢驗(yàn)擬合結(jié)果,最終獲得多重回歸結(jié)果的平均值[27,33]。BRT法不僅可獲得自變量對(duì)因變量影響的相對(duì)貢獻(xiàn)率,而且還可以獲得在其它自變量不變或取均值的前提下,因變量隨某一自變量變化而變化的邊際效應(yīng)[27]。

    BRT模型有4個(gè)重要的輸入?yún)?shù)需要設(shè)置:損失函數(shù)的形式(distribution),學(xué)習(xí)速率(learning.rate),再抽樣比率(bag.fraction),交叉驗(yàn)證折數(shù)(cv.folds)[24,32,34]。損失函數(shù)的形式容易設(shè)定,分類問(wèn)題選擇bernoulli分布,回歸問(wèn)題選擇gaussian分布[24]。學(xué)習(xí)速率決定模型達(dá)到最優(yōu)時(shí)所需的訓(xùn)練時(shí)間,學(xué)習(xí)速率值過(guò)小,則收斂速度慢、 訓(xùn)練時(shí)間長(zhǎng);學(xué)習(xí)速率值過(guò)大,則容易在抽樣時(shí)產(chǎn)生過(guò)多噪音以及使得響應(yīng)函數(shù)平滑性降低、穩(wěn)定性差[34]。交叉驗(yàn)證折數(shù)是指將全部數(shù)據(jù)隨機(jī)k等分,k-1部分作為訓(xùn)練數(shù)據(jù),剩下的部分作為測(cè)試數(shù)據(jù)。再抽樣比率是實(shí)際用于BRT模型建模的數(shù)據(jù)占交叉驗(yàn)證中訓(xùn)練數(shù)據(jù)的比例,從訓(xùn)練數(shù)據(jù)中再抽樣不僅減少了運(yùn)算量,而且在一定程度上避免了過(guò)度學(xué)習(xí)。依據(jù)De′Ath[32]和Elith等[24]推薦的經(jīng)驗(yàn)法則,并結(jié)合本研究的情況,BRT模型參數(shù)設(shè)置為:①損失函數(shù)的形式為“gaussian”;②學(xué)習(xí)速率為0.005;③交叉驗(yàn)證折數(shù)為8;④再抽樣比率為0.5。本文在統(tǒng)計(jì)軟件R 3.3.0下調(diào)用“gbm”包進(jìn)行BRT統(tǒng)計(jì)分析,此外,使用到的軟件還有Excel 2013。

    通過(guò)使用BRT法,可以量化林分特征因子(起源和齡組)、地形因子(地貌、海拔、坡向、坡位和坡度)以及土壤因子(土壤種類、土壤厚度和腐殖層厚度)分別對(duì)地上BCEF以及地下BCEF影響的相對(duì)貢獻(xiàn)率,并獲取相應(yīng)的邊際效應(yīng)。

    3 結(jié)果分析

    3.1 多重因素分別對(duì)地上BCEF和地下BCEF影響的相對(duì)貢獻(xiàn)率

    由圖1A可以看出,齡組對(duì)地上BCEF的影響最大,相對(duì)貢獻(xiàn)率為41.13%,其次為坡向(14.52%)和海拔(10.57%),土類最小(0.18%);由圖1B可以看出,齡組對(duì)地下BCEF的影響最大(41.54%),坡向(15.16%)和腐殖層厚度(9.13%)次之,土類最小(0.37%)。值得一提的是,海拔對(duì)地上BCEF的影響力(10.57%)明顯大于地下BCEF(8.87%),這可能是由于馬尾松地上部分和地下部分適應(yīng)生境變化的能力不同引起的。本文中馬尾松海拔差異較大(55—1412 m,表2),海拔不同,導(dǎo)致熱、水、光和氣等的綜合效應(yīng)產(chǎn)生差異[35- 36],植物的生境發(fā)生變化。

    圖1 各自變量分別對(duì)地上BCEF和地下BCEF影響的相對(duì)貢獻(xiàn)率Fig.1 Relative influences of predictor variables affecting aboveground BCEF and underground BCEF

    3.2 多重因素對(duì)地上BCEF影響的邊際效應(yīng)

    林分特征因子對(duì)地上BCEF影響的相對(duì)貢獻(xiàn)率為47.98%,其中,齡組41.13%,起源6.85%(圖1)。地上BCEF隨齡組由幼齡向過(guò)熟變化而減小,幼齡林最大、過(guò)熟林最小(圖2)。此外,從圖2可以看出,人工林地上BCEF顯著大于天然林,與之相似的結(jié)果也出現(xiàn)在其它森林類型的研究中,羅云建等[37]發(fā)現(xiàn)落葉松人工林地上BCEF顯著大于天然林,羅云建等[37]認(rèn)為光照是引起差異的主要原因,一般而言,天然林通常為異林齡,人工林通常為同齡林[38],兩者的林分結(jié)構(gòu)存在差異,受光不同,光照調(diào)控植物不同器官的生長(zhǎng),進(jìn)而影響著樹(shù)木的材積和生物量。

    地形因子對(duì)地上BCEF影響的相對(duì)貢獻(xiàn)率為39.23%,其中,坡向最大(14.52%),其次為海拔(10.57%)、坡度(8.93%)、坡位(4.53%),地貌最小(0.68%)(圖1)。對(duì)于坡向,地上BCEF在東南坡最大,在西坡最小(圖2)。對(duì)于海拔,地上BCEF隨海拔的升高而升高,并在海拔大于1000m后趨向平穩(wěn)(圖2)。一般海拔每升高100m氣溫下降0.5—0.6 ℃[39],隨海拔上升,環(huán)境條件逐漸趨于惡劣[35],植物的地上部分通過(guò)加強(qiáng)同化能力、積累有機(jī)物來(lái)適應(yīng)生境變化,故地上BCEF隨著海拔的升高而升高;一般而言,海拔1000—1200m是馬尾松分布的低溫極限,然而,少數(shù)福建省馬尾松林可以分布在海拔1100m以上[39],極端的生境條件下植物生長(zhǎng)受限可能是地上BCEF在海拔大于1000m后趨向平穩(wěn)的原因。此外,對(duì)于坡度,當(dāng)坡度小于30°時(shí),地上BCEF相對(duì)穩(wěn)定,當(dāng)處于30°—35°之間時(shí),地上BCEF先升高后下降,大于35°后,地上BCEF不再變化;對(duì)于坡位,地上BCEF在上坡處最大、中坡處最小;對(duì)于地貌,丘陵處的地上BCEF較中山和低山的大(圖2)。

    土壤因子對(duì)地上BCEF影響的相對(duì)貢獻(xiàn)率為12.80%,其中,腐殖層厚度最大(9.02%),土壤厚度次之(3.60%),而土壤種類最小,僅為0.68%(圖1)。從圖2可知,當(dāng)腐殖層厚度小于10cm時(shí),地上BCEF隨腐殖層厚度的增加而減小,超過(guò)10cm后,地上BCEF不再變化并趨于平穩(wěn);當(dāng)土壤厚度小于50cm時(shí),地上BCEF相對(duì)穩(wěn)定,當(dāng)介于50—110cm之間時(shí),地上BCEF先升高后下降,大于110cm后,地上BCEF不再變化;土壤種類對(duì)地上BCEF無(wú)影響。

    圖2 各自變量對(duì)地上BCEF影響的邊際效應(yīng)Fig.2 Marginal effect of predictor variables on aboveground BCEF坡向:ENs(東北坡),Es(東坡),ESs(東南坡),Ns(北坡),Ss(南坡),WNs(西北坡),Ws(西坡),WSs(西南坡)

    3.3 多重因素對(duì)地下BCEF影響的邊際效應(yīng)

    林分特征因子對(duì)地下BCEF影響的相對(duì)貢獻(xiàn)率為48.63%,其中,齡組和起源分別為41.54%和7.09%(圖1)。地下BCEF隨齡組由幼齡向過(guò)熟變化而增大,幼齡林最小、過(guò)熟林最大(圖3)。對(duì)于起源,人工林地下BCEF明顯小于天然林地下BCEF(圖3),地下BCEF隨起源而變化的趨勢(shì)與地上BCEF相反,這可能是因?yàn)椴煌鹪吹闹参锸芄獠煌?進(jìn)而調(diào)控有機(jī)物向地上部分和地下部分分配的策略不同。

    地形因子對(duì)地下BCEF影響的相對(duì)貢獻(xiàn)率為37.97%,其中,坡向最大(15.16%),坡度(8.96%)、海拔(8.87%)、坡位(4.30%)次之,地貌最小(0.68%)(圖1)。從圖3可知,對(duì)于坡向,地下BCEF在東南坡最小,在西坡最大;對(duì)于海拔,地下BCEF隨著海拔的升高而減小,在海拔大于1000m后逐漸趨向平穩(wěn);對(duì)于坡度,當(dāng)坡度小于30°時(shí),地下BCEF相對(duì)穩(wěn)定,處于30°—35°之間時(shí),地下BCEF先下降后升高,大于35°后,地下BCEF不再變化;對(duì)于坡位,地下BCEF在上坡處最小、中坡處最大;對(duì)于地貌,丘陵處的地下BCEF較中山和地山的小。

    土壤因子對(duì)地下BCEF影響的相對(duì)貢獻(xiàn)率為13.41%,其中,腐殖層厚度最大(9.13%),土壤厚度次之(3.91%),土壤種類最小,僅為0.37%,(圖1)。從圖3可知,對(duì)于腐殖層,當(dāng)腐殖層厚度小于10cm時(shí),地下BCEF隨腐殖層厚度的增加而增大,超過(guò)10cm后,地下BCEF不再變化并趨于平穩(wěn);對(duì)于土壤厚度,當(dāng)土壤厚度小于50cm時(shí),地下BCEF一直處于穩(wěn)定狀態(tài),介于50—100cm之間時(shí),地下BCEF先下降后升高,大于100cm后,地下BCEF不再變化;對(duì)于土壤種類,紅壤處的地下BCEF最大,黃壤處的地下BCEF最小。

    圖3 各自變量對(duì)地下BCEF影響的邊際效應(yīng)Fig.3 Marginal effect of predictor variables on underground BCEF

    4 結(jié)論與討論

    利用森林資源調(diào)查數(shù)據(jù),探究林分特征因子、地形因子以及土壤因子如何協(xié)同影響林分BCEF以及多重因素對(duì)林分BCEF影響的相對(duì)貢獻(xiàn)率,可為森林生物量和碳儲(chǔ)量的準(zhǔn)確估算提供科學(xué)根據(jù)[14,26]。本文研究結(jié)果表明,林分特征因子和地形因子是影響地上BCEF以及地下BCEF的主導(dǎo)因素,二者對(duì)地上BCEF影響的相對(duì)貢獻(xiàn)率之和為87.20%、地下BCEF為86.59%。其中,齡組是林分特征因子中影響地上BCEF以及地下BCEF的最大因素,也是所有影響因素中對(duì)地上BCEF以及地下BCEF影響相對(duì)貢獻(xiàn)率最大的因素(41.13%和41.54%)(圖1);坡向是地形因子中影響地上BCEF以及地下BCEF的最大因素(14.52%和15.16%)(圖1)。土壤因子對(duì)地上BCEF以及地下BCEF的影響都較小(12.80%和13.41%),腐殖層厚度是土壤因子中影響地上BCEF以及地下BCEF的最大因素(9.02%和9.13%)。

    植物的生長(zhǎng)與環(huán)境密切相關(guān),面對(duì)復(fù)雜、多變的環(huán)境,植物必須不斷地調(diào)節(jié)自身內(nèi)部機(jī)能,以實(shí)現(xiàn)對(duì)生境的適應(yīng)[40]。地下部分為地上部分提供水、礦質(zhì)元素,而地上部分為地下部分提供有機(jī)物[41],當(dāng)植物所處的生境發(fā)生變化時(shí),調(diào)控有機(jī)物向地上部分和地下部分分配的策略也會(huì)不同。本文發(fā)現(xiàn),當(dāng)林分的齡組依次從幼齡向過(guò)熟變化時(shí),地上BCEF逐漸減小,而地下BCEF的變化趨勢(shì)恰好相反。地上BCEF隨著林齡的增加而減小,這與羅云建等[19]、Pajtík等[42]、González-García等[20]的研究結(jié)果一致;地下BCEF隨著林齡的增加而增大,這與Lehtonen等[18]的研究結(jié)果一致。林木各器官生物量的累積隨著林齡的變化而發(fā)生變化[26],當(dāng)林分處于中幼齡階段時(shí),林木間的相互競(jìng)爭(zhēng)以及林分郁閉,需要更多的光合作用產(chǎn)生的有機(jī)物以促進(jìn)林分的快速生長(zhǎng),此時(shí)林分地上部分生長(zhǎng)迅速,地下部分的生長(zhǎng)受到一定的限制[43]。當(dāng)林分處于成過(guò)熟階段時(shí),林分地上部分的細(xì)胞分化成熟、不再具有分裂能力[41],地上部分的生長(zhǎng)就生長(zhǎng)放緩或停止;而植物的地下部分——根端具有特殊的頂端分生組織[41],可以無(wú)限分裂,故地下部分可以繼續(xù)生長(zhǎng)。地上BCEF在東南坡最大、西坡最小,這不同于羅云建[19]、吳小山[22]的研究結(jié)果,羅云建[19]、吳小山[22]發(fā)現(xiàn)不同坡向間的地上BCEF無(wú)顯著差異;地下BCEF在東南坡最小、西坡最大。坡向的差異會(huì)導(dǎo)致光照強(qiáng)度、水分、溫度、土壤濕度等生態(tài)環(huán)境的變化,不同坡向間日照時(shí)長(zhǎng)、生長(zhǎng)季長(zhǎng)度等因素也存在差異,從而對(duì)植物個(gè)體生長(zhǎng)產(chǎn)生影響[44- 45]。

    齡組對(duì)地上BCEF和地下BCEF的影響均極大,且當(dāng)林分分別處于幼齡和成過(guò)熟階段時(shí),幼齡林的地上BCEF會(huì)遠(yuǎn)遠(yuǎn)大于成過(guò)熟林的地上BCEF,幼齡林的地下BCEF會(huì)小于成過(guò)熟林的地下BCEF。本文建議,在進(jìn)行外業(yè)調(diào)查估算森林生物量和碳儲(chǔ)量時(shí),應(yīng)按照各齡組的林分比例進(jìn)行采樣,依據(jù)齡組計(jì)算相應(yīng)的BCEF或者建立林齡普適的BCEF模型。此外,在采樣時(shí),應(yīng)綜合考慮生物因素和非生物因素對(duì)林分的影響,盡可能選擇林分條件、地形、土壤等方面具有廣泛代表性以及多樣性的樣地,盡量避免只選生長(zhǎng)良好、環(huán)境條件單一的林分[26]。本文僅研究了土壤因子、地形因子以及林分特征因子對(duì)各項(xiàng)BCEF的影響,然而林木生長(zhǎng)還會(huì)受到諸如氣候因子、火因子、生物因子以及人為因子的影響[26],有待深入研究。

    [1] Keith H, Mackey B, Berry S, Lindenmayer D, Gibbons P. Estimating carbon carrying capacity in natural forest ecosystems across heterogeneous landscapes: addressing sources of error. Global Change Biology, 2010, 16(11): 2971- 2989.

    [2] 雷相東, 張會(huì)儒, 牟惠生. 東北過(guò)伐林區(qū)蒙古櫟林分相容性生物量模型研究. 第四紀(jì)研究, 2010, 30(3): 559- 565.

    [3] Temesgen H, Affleck D, Poudel K, Gray A, Sessions J. A review of the challenges and opportunities in estimating above ground forest biomass using tree-level models. Scandinavian Journal of Forest Research, 2015, 30(4): 326- 335.

    [4] Timmons D S, Buchholz T, Veeneman C H. Forest biomass energy: Assessing atmospheric carbon impacts by discounting future carbon flows. GCB Bioenergy, 2016, 8(3): 631- 643.

    [5] 左舒翟, 任引, 王效科, 張小全, 羅云建. 中國(guó)杉木林生物量估算參數(shù)及其影響因素. 林業(yè)科學(xué), 2014, 50(11): 1- 12.

    [6] 羅云建, 張小全, 王效科, 朱建華, 侯振宏, 張治軍. 森林生物量的估算方法及其研究進(jìn)展. 林業(yè)科學(xué), 2009, 45(8): 129- 134.

    [7] Schroeder P, Brown S, Mo J M, Birdsey R, Cieszewski C. Biomass estimation for temperate broadleaf forests of the united states using inventory data. Forest Science, 1997, 43(3): 424- 434.

    [8] Somogyi Z, Cienciala E, M?kip?? R, Muukkonen P, Lehtonen A, Weiss P. Indirect methods of large-scale forest biomass estimation. European Journal of Forest Research, 2007, 126(2): 197- 207.

    [9] Tobin B, Nieuwenhuis M. Biomass expansion factors for Sitka spruce (Piceasitchensis(Bong.) Carr.) in Ireland. European Journal of Forest Research, 2007, 126(2): 189- 196.

    [11] Magalh?es T M, Seifert T. Tree component biomass expansion factors and root-to-shoot ratio of Lebombo ironwood: measurement uncertainty. Carbon Balance and Management, 2015, 10(1): 1- 14.

    [12] Tolunay D. Coefficients that can be used to calculate biomass carbon storage of Turkish forests. Journal of the Faculty of Forestry Istanbul University, 2017, 67(2), doi: 10.17099/jffiu.xxxxx.

    [13] Di Cosmo L, Gasparini P, Tabacchi G. A national-scale, stand-level model to predict total above-ground tree biomass from growing stock volume. Forest Ecology and Management, 2016, 361: 269- 276.

    [14] Lam T Y, Li X D, Kim R H, Lee K H, Son Y M. Bayesian meta-analysis of regional biomass factors forQuercusmongolicaforests in South Korea. Journal of Forestry Research, 2015, 26(4): 875- 885.

    [15] Soares P, Tomé M. Analysis of the effectiveness of biomass expansion factors to estimate stand biomass//Proceedings of the International Conference on Modeling Forest Production. Austria, 2004.

    [16] 方精云, 劉國(guó)華, 徐嵩齡. 我國(guó)森林植被的生物量和凈生產(chǎn)量. 生態(tài)學(xué)報(bào), 1996, 16(5): 497- 508.

    [17] 劉國(guó)華, 傅伯杰, 方精云. 中國(guó)森林碳動(dòng)態(tài)及其對(duì)全球碳平衡的貢獻(xiàn). 生態(tài)學(xué)報(bào), 2000, 20(5): 733- 740.

    [18] Lehtonen A, M?kip?? R, Heikkinen J, Siev?nen R, Liski J. Biomass expansion factors (BEFs) for Scots pine, Norway spruce and birch according to stand age for boreal forests. Forest Ecology and Management, 2004, 188(1/3): 211- 224.

    [19] 羅云建. 華北落葉松人工林生物量碳計(jì)量參數(shù)研究[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2007.

    [20] González-García M, Hevia A, Majada J, Barrio-Anta M. Above-ground biomass estimation at tree and stand level for short rotation plantations ofEucalyptusnitens(Deane & Maiden) Maiden in Northwest Spain. Biomass and Bioenergy, 2013, 54: 147- 157.

    [21] Teobaldelli M, Somogyi Z, Migliavacca M, Usoltsev V A. Generalized functions of biomass expansion factors for conifers and broadleaved by stand age, growing stock and site index. Forest Ecology and Management, 2009, 257(3): 1004- 1013.

    [22] 吳小山. 楊樹(shù)人工林生物量碳計(jì)量參數(shù)研究[D]. 雅安: 四川農(nóng)業(yè)大學(xué), 2008.

    [23] 黃云奉, 劉屹, 黃世友, 何邦亮, 陳麗潔, 李燕瓊. 不同林齡馬尾松生長(zhǎng)及生物量分配研究. 四川林業(yè)科技, 2015, 36(4): 72- 75.

    [24] Elith J, Leathwick J R, Hastie T. A working guide to boosted regression trees. Journal of Animal Ecology, 2008, 77(4): 802- 813.

    [25] Weisberg P J, Shandra O, Becker M E. Landscape influences on recent timberline shifts in the Carpathian Mountains: Abiotic influences modulate effects of land-use change. Arctic, Antarctic, and Alpine Research, 2013, 45(3): 404- 414.

    [26] 閆晶, 羅云建, 鄭德福, 王水城. 桉樹(shù)生物量估算差異的源解析. 林業(yè)科學(xué), 2014, 50(2): 92- 98.

    [27] 焦琳琳, 常禹, 申丹, 胡遠(yuǎn)滿, 李春林, 馬俊. 利用增強(qiáng)回歸樹(shù)分析中國(guó)野火空間分布格局的影響因素. 生態(tài)學(xué)雜志, 2015, 34(8): 2288- 2296.

    [28] 曾偉生. 全國(guó)立木生物量方程建模方法研究[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2011.

    [29] Li H K, Zhao P X. Improving the accuracy of tree-level aboveground biomass equations with height classification at a large regional scale. Forest Ecology and Management, 2013, 289: 153- 163.

    [30] 傅煜, 雷淵才, 曾偉生. 區(qū)域尺度杉木生物量估計(jì)的不確定性度量. 林業(yè)科學(xué), 2014, 50(12): 79- 86.

    [31] 傅煜, 雷淵才, 曾偉生. 單木生物量模型估計(jì)區(qū)域尺度生物量的不確定性. 生態(tài)學(xué)報(bào), 2015, 35(23): 7738- 7747.

    [32] De′Ath G. Boosted trees for ecological modeling and prediction. Ecology, 2007, 88(1): 243- 251.

    [33] 李春林, 劉淼, 胡遠(yuǎn)滿, 徐巖巖, 孫鳳云. 基于增強(qiáng)回歸樹(shù)和Logistic回歸的城市擴(kuò)展驅(qū)動(dòng)力分析. 生態(tài)學(xué)報(bào), 2014, 34(3): 727- 737.

    [34] Zhang Y, Chen H Y H, Reich P B. Forest productivity increases with evenness, species richness and trait variation: a global meta-analysis. Journal of Ecology, 2012, 100(3): 742- 749.

    [35] 劉彥春, 張遠(yuǎn)東, 劉世榮, 張笑鶴. 川西亞高山針闊混交林喬木層生物量、生產(chǎn)力隨海拔梯度的變化. 生態(tài)學(xué)報(bào), 2010, 30(21): 5810- 5820.

    [36] 吳棟棟, 周永斌, 于大炮, 戴冠華. 不同海拔長(zhǎng)白山岳樺的生理變化. 生態(tài)學(xué)報(bào), 2009, 29(5): 2279- 2285.

    [37] 羅云建, 張小全, 侯振宏, 于彭濤, 朱建華. 我國(guó)落葉松林生物量碳計(jì)量參數(shù)的初步研究. 植物生態(tài)學(xué)報(bào), 2007, 31(6): 1111- 1118.

    [38] 羅云建, 王效科, 張小全, 朱建華, 張治軍, 侯振宏. 華北落葉松人工林的生物量估算參數(shù). 林業(yè)科學(xué), 2010, 46(2): 6- 11.

    [39] 張大鵬. 不同海拔高度對(duì)馬尾松莖木質(zhì)部結(jié)構(gòu)的影響. 生態(tài)學(xué)雜志, 1994, 13(2): 12- 18.

    [40] 姜漢僑, 段昌群, 楊樹(shù)華, 王崇云, 蘇文華. 植物生態(tài)學(xué). 北京: 高等教育出版社, 2010.

    [41] 胡金良. 植物學(xué). 北京: 中國(guó)農(nóng)業(yè)大學(xué)出版社, 2012.

    [42] Pajtík J, Konpka B, Lukac M. Individual biomass factors for beech, oak and pine in Slovakia: a comparative study in young naturally regenerated stands. Trees, 2011, 25(2): 277- 288.

    [43] Vanninen P, Ylitalo H, Siev?nen R, M?kel? A. Effects of age and site quality on the distribution of biomass in Scots pine (PinussylvestrisL.). Trees, 1996, 10(4): 231- 238.

    [44] 郭濱德, 張遠(yuǎn)東, 王曉春. 川西高原不同坡向云、 冷杉樹(shù)輪對(duì)快速升溫的響應(yīng)差異. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(2): 354- 364.

    [45] 王富剛, 徐偉洲, 亢福仁, 尚愛(ài)軍, 張靜. 不同坡向條件下小葉楊光合光響應(yīng)曲線研究. 水土保持研究, 2015, 22(6): 177- 182.

    [46] 李俊清. 森林生態(tài)學(xué)(第二版). 北京: 高等教育出版社, 2010.

    FactorsaffectingthebiomassconversionandexpansionfactorofmassonpineinFujianProvince

    OU Qiangxin1, LI Haikui1, *, YANG Ying2

    1ResearchInstituteofForestResourceInformationTechniques,ChineseAcademyofForestry,Beijing100091,China2AcademyofForestryInventoryandPlanning,StateForestryAdministration,Beijing100714,China

    國(guó)家自然科學(xué)基金 (31370634);中國(guó)清潔發(fā)展機(jī)制基金贈(zèng)款項(xiàng)目:2020年后林業(yè)增匯減排的行動(dòng)目標(biāo)研究(2013014)

    2016- 06- 13; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 24

    10.5846/stxb201606131140

    *通訊作者Corresponding author.E-mail: lihk@ifrit.ac.cn

    歐強(qiáng)新,李海奎,楊英.福建地區(qū)馬尾松生物量轉(zhuǎn)換和擴(kuò)展因子的影響因素.生態(tài)學(xué)報(bào),2017,37(17):5756- 5764.

    Ou Q X, Li H K, Yang Y.Factors affecting the biomass conversion and expansion factor of masson pine in Fujian Province.Acta Ecologica Sinica,2017,37(17):5756- 5764.

    猜你喜歡
    坡向馬尾松林分
    撫育間伐對(duì)油松林下灌木多樣性的影響
    馬尾松栽培技術(shù)及撫育管理
    綠色科技(2019年5期)2019-11-29 13:17:37
    DEM地表坡向變率的向量幾何計(jì)算法
    4種闊葉混交林的持水特性研究
    綠色科技(2019年6期)2019-04-12 05:38:42
    4種人工林的土壤化學(xué)性質(zhì)和酶活性特征研究
    綠色科技(2019年6期)2019-04-12 05:38:42
    馬尾松果糖-1,6-二磷酸酶基因克隆及表達(dá)模式分析
    馬尾松初級(jí)種子園復(fù)壯技術(shù)
    24年生馬尾松種子園自由授粉子代測(cè)定及家系選擇
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    亞熱帶地區(qū)典型林分氮保留能力的差異及δ15N空間垂直分異特征
    女人十人毛片免费观看3o分钟| 少妇人妻一区二区三区视频| 亚洲内射少妇av| 日本av手机在线免费观看| 在线 av 中文字幕| 国产综合懂色| 大香蕉久久网| 亚洲精品,欧美精品| 有码 亚洲区| 毛片一级片免费看久久久久| 能在线免费看毛片的网站| 在现免费观看毛片| 亚洲精品日本国产第一区| 麻豆成人午夜福利视频| 亚洲欧美一区二区三区国产| 毛片一级片免费看久久久久| 少妇裸体淫交视频免费看高清| 视频中文字幕在线观看| 亚洲精品影视一区二区三区av| 亚洲av一区综合| 国产综合精华液| 69av精品久久久久久| 久久99热这里只频精品6学生| 亚洲色图综合在线观看| 亚洲aⅴ乱码一区二区在线播放| 91久久精品国产一区二区成人| 99久久精品国产国产毛片| 国产v大片淫在线免费观看| av在线亚洲专区| 亚洲av国产av综合av卡| 男女国产视频网站| 99热全是精品| 国产探花在线观看一区二区| 男人舔奶头视频| 精品久久久久久久久av| 51国产日韩欧美| 精品久久久噜噜| 高清午夜精品一区二区三区| 五月玫瑰六月丁香| 久久久久网色| 免费人成在线观看视频色| 天堂中文最新版在线下载 | 视频中文字幕在线观看| 欧美丝袜亚洲另类| 国产午夜精品一二区理论片| 自拍偷自拍亚洲精品老妇| 综合色av麻豆| 免费看av在线观看网站| 日韩,欧美,国产一区二区三区| 黄片无遮挡物在线观看| 男男h啪啪无遮挡| 精品99又大又爽又粗少妇毛片| 日韩,欧美,国产一区二区三区| 亚洲精品日韩av片在线观看| 嫩草影院新地址| 蜜桃亚洲精品一区二区三区| 亚洲色图av天堂| 免费播放大片免费观看视频在线观看| 国产免费视频播放在线视频| av在线app专区| av天堂中文字幕网| 黑人高潮一二区| 日本一二三区视频观看| 80岁老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 日本av手机在线免费观看| 国产精品久久久久久久久免| 国产毛片a区久久久久| 国产一区二区在线观看日韩| 性色avwww在线观看| 老师上课跳d突然被开到最大视频| 在线a可以看的网站| 亚洲真实伦在线观看| 久久99热这里只有精品18| 国产精品女同一区二区软件| 亚洲四区av| 国产成人a区在线观看| 日本-黄色视频高清免费观看| 五月玫瑰六月丁香| 麻豆精品久久久久久蜜桃| 国产在线男女| 国产午夜精品一二区理论片| 欧美精品人与动牲交sv欧美| 久久综合国产亚洲精品| 亚洲图色成人| 三级国产精品欧美在线观看| 一级a做视频免费观看| 亚洲人与动物交配视频| 亚洲精华国产精华液的使用体验| 自拍偷自拍亚洲精品老妇| 午夜亚洲福利在线播放| 精品人妻一区二区三区麻豆| 久久久久久九九精品二区国产| 国产精品三级大全| 嘟嘟电影网在线观看| 伊人久久精品亚洲午夜| 久久久久久伊人网av| 亚洲最大成人手机在线| 熟女av电影| 国产真实伦视频高清在线观看| 国产黄片视频在线免费观看| 嘟嘟电影网在线观看| 99久久九九国产精品国产免费| 亚洲精品日韩在线中文字幕| 一本色道久久久久久精品综合| 听说在线观看完整版免费高清| 一个人观看的视频www高清免费观看| 日本一二三区视频观看| 在线免费观看不下载黄p国产| 新久久久久国产一级毛片| av专区在线播放| 三级男女做爰猛烈吃奶摸视频| 欧美精品一区二区大全| 日韩中字成人| 各种免费的搞黄视频| 一区二区三区四区激情视频| 毛片女人毛片| 热re99久久精品国产66热6| 亚洲人成网站高清观看| 亚洲最大成人中文| 久久精品人妻少妇| 亚洲欧洲日产国产| 又爽又黄无遮挡网站| 国产精品一区二区三区四区免费观看| 爱豆传媒免费全集在线观看| 日韩一区二区三区影片| 国产精品.久久久| 高清日韩中文字幕在线| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 高清日韩中文字幕在线| 99热这里只有是精品在线观看| 亚洲在久久综合| 日韩精品有码人妻一区| 亚洲国产精品成人久久小说| kizo精华| av免费在线看不卡| 少妇熟女欧美另类| 欧美最新免费一区二区三区| 国产亚洲91精品色在线| 成人二区视频| 亚洲av日韩在线播放| 国国产精品蜜臀av免费| 亚洲人成网站在线播| 99视频精品全部免费 在线| 国产黄a三级三级三级人| 免费在线观看成人毛片| 国产熟女欧美一区二区| 久久精品综合一区二区三区| 亚洲不卡免费看| 欧美日韩一区二区视频在线观看视频在线 | 国产黄频视频在线观看| 日本一二三区视频观看| 日韩国内少妇激情av| 免费av毛片视频| 国产淫片久久久久久久久| 亚洲国产精品999| 成人国产麻豆网| 亚洲精品国产av成人精品| 男女无遮挡免费网站观看| 国产亚洲5aaaaa淫片| 久久97久久精品| 尤物成人国产欧美一区二区三区| 久久久午夜欧美精品| av在线播放精品| av在线app专区| 天天躁夜夜躁狠狠久久av| 嫩草影院新地址| 国产大屁股一区二区在线视频| 22中文网久久字幕| 好男人视频免费观看在线| 国产淫语在线视频| 日韩精品有码人妻一区| 免费在线观看成人毛片| 在线天堂最新版资源| 国国产精品蜜臀av免费| 新久久久久国产一级毛片| 久久久久久九九精品二区国产| 精品少妇久久久久久888优播| 亚洲天堂国产精品一区在线| 亚洲成人精品中文字幕电影| 热99国产精品久久久久久7| 国产又色又爽无遮挡免| 国产精品一及| 久热久热在线精品观看| 中文天堂在线官网| 王馨瑶露胸无遮挡在线观看| 国产在线男女| 日韩三级伦理在线观看| 久久久久久久大尺度免费视频| 午夜福利在线观看免费完整高清在| 日韩强制内射视频| 老女人水多毛片| 国产精品国产三级专区第一集| 色播亚洲综合网| 精品一区二区三卡| 国产精品秋霞免费鲁丝片| 午夜爱爱视频在线播放| 欧美区成人在线视频| 精品久久国产蜜桃| 观看免费一级毛片| 五月伊人婷婷丁香| 51国产日韩欧美| 成人免费观看视频高清| 夫妻午夜视频| 国产综合懂色| 26uuu在线亚洲综合色| 免费观看的影片在线观看| 亚洲欧美成人精品一区二区| 又粗又硬又长又爽又黄的视频| 性色av一级| 99视频精品全部免费 在线| 久久久精品欧美日韩精品| 99久久人妻综合| 国产乱人视频| 丰满少妇做爰视频| 久久鲁丝午夜福利片| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 欧美精品一区二区大全| 丝袜喷水一区| 三级男女做爰猛烈吃奶摸视频| 少妇丰满av| 成年人午夜在线观看视频| 大又大粗又爽又黄少妇毛片口| 成人黄色视频免费在线看| 婷婷色综合www| 熟女av电影| 99久久中文字幕三级久久日本| 三级国产精品片| 国产有黄有色有爽视频| av在线天堂中文字幕| 国产高清三级在线| 亚洲怡红院男人天堂| 在线免费观看不下载黄p国产| 老司机影院毛片| 久久影院123| 免费黄频网站在线观看国产| 美女主播在线视频| 性色av一级| 亚洲欧美精品自产自拍| 欧美成人a在线观看| 午夜视频国产福利| 亚州av有码| 国产精品久久久久久精品古装| 三级国产精品欧美在线观看| 一区二区三区乱码不卡18| 欧美人与善性xxx| 大片电影免费在线观看免费| 亚洲久久久久久中文字幕| 22中文网久久字幕| 国产精品嫩草影院av在线观看| 欧美激情久久久久久爽电影| 18禁在线播放成人免费| 欧美变态另类bdsm刘玥| 自拍欧美九色日韩亚洲蝌蚪91 | 久久综合国产亚洲精品| 国产午夜精品一二区理论片| 国产国拍精品亚洲av在线观看| 免费黄网站久久成人精品| 午夜日本视频在线| 亚洲精品国产av成人精品| 国产成人freesex在线| 国国产精品蜜臀av免费| 久久久a久久爽久久v久久| 亚洲精品日韩av片在线观看| 成年免费大片在线观看| 1000部很黄的大片| 99久国产av精品国产电影| 亚洲最大成人中文| 欧美日韩视频高清一区二区三区二| 男人添女人高潮全过程视频| 日韩三级伦理在线观看| videos熟女内射| 3wmmmm亚洲av在线观看| 干丝袜人妻中文字幕| 久久热精品热| 国产片特级美女逼逼视频| 黄色日韩在线| 亚洲色图综合在线观看| 久久精品国产亚洲av涩爱| 97在线视频观看| 国产色爽女视频免费观看| 又大又黄又爽视频免费| 成人亚洲精品av一区二区| 蜜臀久久99精品久久宅男| 99久国产av精品国产电影| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 中文资源天堂在线| 91久久精品国产一区二区三区| av在线老鸭窝| 99热网站在线观看| 免费黄网站久久成人精品| 国产黄频视频在线观看| 久久精品国产亚洲网站| 91午夜精品亚洲一区二区三区| 日本一二三区视频观看| 在线天堂最新版资源| 国产免费视频播放在线视频| 九草在线视频观看| 观看美女的网站| 精品久久久久久久人妻蜜臀av| av免费在线看不卡| 涩涩av久久男人的天堂| 亚洲国产av新网站| 99热这里只有精品一区| 在线观看一区二区三区| 国产精品蜜桃在线观看| 97人妻精品一区二区三区麻豆| 少妇熟女欧美另类| 亚洲精品日韩在线中文字幕| 禁无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 啦啦啦中文免费视频观看日本| 亚洲欧美成人精品一区二区| 国产高清国产精品国产三级 | 中文天堂在线官网| 麻豆成人午夜福利视频| 欧美精品一区二区大全| av女优亚洲男人天堂| 最近中文字幕2019免费版| 91在线精品国自产拍蜜月| 欧美老熟妇乱子伦牲交| 国产精品爽爽va在线观看网站| 精品视频人人做人人爽| 国产成人精品婷婷| 免费人成在线观看视频色| 亚洲精品色激情综合| 亚洲图色成人| 美女国产视频在线观看| 久久久久久久久久久丰满| 黄片wwwwww| 婷婷色av中文字幕| 免费观看a级毛片全部| 色视频www国产| 日韩av不卡免费在线播放| 国内精品宾馆在线| 不卡视频在线观看欧美| 黑人高潮一二区| 久久精品夜色国产| 久热这里只有精品99| 22中文网久久字幕| 久久久久久国产a免费观看| 免费黄色在线免费观看| 小蜜桃在线观看免费完整版高清| 26uuu在线亚洲综合色| 欧美日韩国产mv在线观看视频 | 欧美97在线视频| 日韩人妻高清精品专区| 久久人人爽av亚洲精品天堂 | 一级av片app| 黄色日韩在线| 免费av毛片视频| 国产黄色视频一区二区在线观看| 在线播放无遮挡| 亚洲av福利一区| 老女人水多毛片| 国产成人免费观看mmmm| 亚洲人成网站高清观看| 亚洲国产欧美在线一区| 久久精品熟女亚洲av麻豆精品| 亚洲精品乱久久久久久| 中文字幕久久专区| 久久精品综合一区二区三区| 亚洲欧美成人精品一区二区| 超碰av人人做人人爽久久| 久久久欧美国产精品| 91久久精品国产一区二区三区| 国产在线男女| 欧美性感艳星| 国产精品99久久久久久久久| 亚洲,一卡二卡三卡| 欧美人与善性xxx| 少妇裸体淫交视频免费看高清| 久久久久久久久久成人| 免费观看性生交大片5| 亚洲在久久综合| 国产黄片视频在线免费观看| 精品酒店卫生间| 久久韩国三级中文字幕| 夫妻性生交免费视频一级片| 丰满少妇做爰视频| 国产免费福利视频在线观看| 亚洲成人久久爱视频| 91久久精品电影网| 一级黄片播放器| 亚洲av中文av极速乱| 最近的中文字幕免费完整| 亚洲一区二区三区欧美精品 | 欧美bdsm另类| 国产乱人偷精品视频| 久久久成人免费电影| 亚洲天堂国产精品一区在线| 中文精品一卡2卡3卡4更新| 成人综合一区亚洲| 日韩av免费高清视频| 欧美潮喷喷水| 中国国产av一级| 人妻一区二区av| 日本爱情动作片www.在线观看| 成人高潮视频无遮挡免费网站| 69av精品久久久久久| 亚洲av二区三区四区| 丰满少妇做爰视频| 麻豆久久精品国产亚洲av| 五月天丁香电影| 1000部很黄的大片| 久久精品国产亚洲网站| 亚洲精品成人久久久久久| 亚洲人成网站在线播| 精品一区二区免费观看| 黄色欧美视频在线观看| 日韩强制内射视频| 人妻制服诱惑在线中文字幕| 97热精品久久久久久| 精品熟女少妇av免费看| 欧美xxxx性猛交bbbb| 只有这里有精品99| 如何舔出高潮| 国产精品一及| 午夜精品一区二区三区免费看| 精品久久国产蜜桃| 国产成人一区二区在线| 少妇猛男粗大的猛烈进出视频 | 97精品久久久久久久久久精品| 午夜日本视频在线| 极品教师在线视频| 菩萨蛮人人尽说江南好唐韦庄| av黄色大香蕉| 午夜激情久久久久久久| 午夜免费男女啪啪视频观看| 免费观看性生交大片5| h日本视频在线播放| 国语对白做爰xxxⅹ性视频网站| 肉色欧美久久久久久久蜜桃 | 久久国产乱子免费精品| 精品亚洲乱码少妇综合久久| 国产成人精品婷婷| 亚洲无线观看免费| 男男h啪啪无遮挡| 高清午夜精品一区二区三区| 国产黄色视频一区二区在线观看| 欧美日韩视频精品一区| 精品人妻一区二区三区麻豆| 麻豆国产97在线/欧美| 中国美白少妇内射xxxbb| 亚洲精品久久午夜乱码| 人体艺术视频欧美日本| 成年人午夜在线观看视频| 神马国产精品三级电影在线观看| 中文字幕制服av| 亚洲伊人久久精品综合| 国产淫语在线视频| 欧美日韩视频高清一区二区三区二| 成人免费观看视频高清| 国产精品嫩草影院av在线观看| 亚洲精品乱码久久久久久按摩| 亚洲成人一二三区av| 日本wwww免费看| 国产免费又黄又爽又色| 久久久精品欧美日韩精品| 黄色配什么色好看| 高清欧美精品videossex| 日韩av不卡免费在线播放| 国产精品久久久久久精品古装| 国产伦精品一区二区三区视频9| 亚洲av电影在线观看一区二区三区 | 少妇丰满av| 男女那种视频在线观看| 国产精品爽爽va在线观看网站| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 亚洲一区二区三区欧美精品 | av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 亚洲怡红院男人天堂| 免费黄频网站在线观看国产| 国产视频首页在线观看| 18+在线观看网站| 欧美国产精品一级二级三级 | 日韩电影二区| 黄色怎么调成土黄色| 欧美最新免费一区二区三区| 水蜜桃什么品种好| 亚洲精品日韩av片在线观看| 免费黄频网站在线观看国产| 亚洲内射少妇av| 可以在线观看毛片的网站| 一级二级三级毛片免费看| 亚洲av二区三区四区| 国产片特级美女逼逼视频| 亚洲欧美精品专区久久| 97在线人人人人妻| 亚洲精品国产成人久久av| 大香蕉久久网| 国产淫片久久久久久久久| 亚洲自拍偷在线| 中文字幕人妻熟人妻熟丝袜美| 2018国产大陆天天弄谢| 国产探花极品一区二区| 高清欧美精品videossex| 汤姆久久久久久久影院中文字幕| 丝瓜视频免费看黄片| 精品久久久久久久久av| 日韩,欧美,国产一区二区三区| 午夜福利视频1000在线观看| 亚洲精品亚洲一区二区| 狂野欧美白嫩少妇大欣赏| 男女边摸边吃奶| 午夜激情福利司机影院| 午夜日本视频在线| 最近的中文字幕免费完整| 交换朋友夫妻互换小说| av女优亚洲男人天堂| 精品人妻视频免费看| 乱码一卡2卡4卡精品| 国产高清国产精品国产三级 | 草草在线视频免费看| 美女被艹到高潮喷水动态| 国产av国产精品国产| av在线蜜桃| 欧美3d第一页| 亚洲自偷自拍三级| 欧美人与善性xxx| 春色校园在线视频观看| 黄片无遮挡物在线观看| 男女那种视频在线观看| 午夜免费鲁丝| 三级国产精品欧美在线观看| 如何舔出高潮| 亚洲内射少妇av| 国产成人精品福利久久| 下体分泌物呈黄色| 久久久久精品久久久久真实原创| 男人添女人高潮全过程视频| 亚洲,一卡二卡三卡| 纵有疾风起免费观看全集完整版| 在线亚洲精品国产二区图片欧美 | 在线精品无人区一区二区三 | 日韩电影二区| 69av精品久久久久久| 男的添女的下面高潮视频| 99热这里只有精品一区| 久久久精品94久久精品| 一个人看的www免费观看视频| 狂野欧美激情性bbbbbb| 精品国产三级普通话版| 久久久欧美国产精品| 欧美日韩国产mv在线观看视频 | 青青草视频在线视频观看| 国产老妇女一区| 韩国高清视频一区二区三区| 黑人高潮一二区| 日本熟妇午夜| 国产老妇伦熟女老妇高清| 久久国产乱子免费精品| 国产淫语在线视频| 91精品伊人久久大香线蕉| 亚洲精品乱码久久久久久按摩| 亚洲av在线观看美女高潮| 91久久精品国产一区二区三区| 亚洲成人av在线免费| 人妻系列 视频| 色吧在线观看| av在线亚洲专区| 亚洲av免费高清在线观看| 五月伊人婷婷丁香| 人妻 亚洲 视频| 好男人视频免费观看在线| 蜜臀久久99精品久久宅男| 免费黄频网站在线观看国产| 久久久久久久精品精品| 中文天堂在线官网| 亚洲av在线观看美女高潮| 亚洲av电影在线观看一区二区三区 | 黄色怎么调成土黄色| 麻豆久久精品国产亚洲av| 成年av动漫网址| 国产一区二区亚洲精品在线观看| 久久久久久久久久人人人人人人| 欧美丝袜亚洲另类| 国产成年人精品一区二区| 九草在线视频观看| 少妇人妻精品综合一区二区| 亚洲综合色惰| 国产探花极品一区二区| 黄色视频在线播放观看不卡| 日韩制服骚丝袜av| 久久亚洲国产成人精品v| 热99国产精品久久久久久7| 99热国产这里只有精品6| 国产成人午夜福利电影在线观看| 午夜福利视频1000在线观看| 3wmmmm亚洲av在线观看| 午夜视频国产福利| 亚洲精品国产av成人精品| 欧美日韩精品成人综合77777| 可以在线观看毛片的网站| 欧美日韩国产mv在线观看视频 | 国产av国产精品国产| 亚洲激情五月婷婷啪啪| 国产色爽女视频免费观看| 久久6这里有精品| 日韩av不卡免费在线播放| 久久久a久久爽久久v久久| 亚洲欧美日韩卡通动漫| 男女边摸边吃奶| 又黄又爽又刺激的免费视频.| 亚洲色图综合在线观看| 婷婷色麻豆天堂久久| 午夜精品一区二区三区免费看| 天天躁日日操中文字幕| 只有这里有精品99| 联通29元200g的流量卡|