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

    馬尾松林葉面積指數(shù)季節(jié)動(dòng)態(tài)的測(cè)定

    2024-09-13 00:00:00冉佳璇羅光浪何世強(qiáng)丁貴杰譚偉戚玉嬌
    關(guān)鍵詞:生長(zhǎng)研究

    摘 要:【目的】評(píng)價(jià)光學(xué)儀器法測(cè)定馬尾松林葉面積指數(shù)(LAI)的精度,為準(zhǔn)確監(jiān)測(cè)馬尾松林LAI季節(jié)動(dòng)態(tài)提供技術(shù)支持?!痉椒ā恳缘蚵湮锓y(cè)定值代表馬尾松林真實(shí)的LAI,對(duì)展葉調(diào)查、異速生長(zhǎng)方程法和光學(xué)儀器法(半球攝影法和LAI-2200植物冠層分析儀)獲取的馬尾松林LAI季節(jié)動(dòng)態(tài)進(jìn)行評(píng)估,并量化不同因素包括集聚指數(shù)(ΩE)、木質(zhì)比例(α)和針簇比(γE)對(duì)光學(xué)儀器測(cè)定值的影響。【結(jié)果】馬尾松林內(nèi)主要樹種的平均葉壽命為1.53~1.94,比葉面積(SLA)差異較大(99.15 ~278.33 cm2/g)。不同方法測(cè)定馬尾松林LAI的季節(jié)動(dòng)態(tài)均呈單峰形變化模式,且在7月達(dá)到峰值。在各調(diào)查時(shí)期,凋落物法的測(cè)定值與異速生長(zhǎng)方程法的測(cè)定值差異不顯著。半球攝影法測(cè)定的有效葉面積指數(shù)(Le)比凋落物法、異速生長(zhǎng)方程法的測(cè)定值分別低估了49%~64%和47%~61%;LAI-2200植物冠層分析儀法測(cè)定的Le低估了25%~44%和20%~42%。半球攝影法測(cè)定的Le經(jīng)自動(dòng)曝光、α、ΩE及γE校正后,精度明顯提高,比凋落物法、異速生長(zhǎng)方程法平均低估了8%和5%;LAI-2200植物冠層分析儀法測(cè)定的Le經(jīng)α、ΩE及γE校正后,其精度大幅度提高,比凋落物法、異速生長(zhǎng)方程法平均低估了6%和2%?!窘Y(jié)論】考慮影響光學(xué)儀器法測(cè)定LAI的因素后,半球攝影法和LAI-2200植物冠層分析儀法測(cè)定的LAI精度得到明顯提高,表明兩種光學(xué)儀器法均能相對(duì)準(zhǔn)確地測(cè)定馬尾松林LAI的季節(jié)動(dòng)態(tài),可為不同區(qū)域馬尾松林LAI季節(jié)動(dòng)態(tài)的研究提供方法參考。

    關(guān)鍵詞:馬尾松人工林;葉面積指數(shù);季節(jié)動(dòng)態(tài);凋落物法;半球攝影法;LAI-2200植物冠層分析儀

    中圖分類號(hào):S791.248 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1673-923X(2024)05-0056-11

    基金項(xiàng)目:國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0600302);貴州省科技計(jì)劃項(xiàng)目(黔科合支撐[2019]2874);貴州省一流學(xué)科建設(shè)項(xiàng)目(GNYL[2017]007)。

    Determination of seasonal dynamics of leaf area index in the Pinus massoniana forest

    RAN Jiaxuan1a, LUO Guanglang2, HE Shiqiang3, DING Guijie1a,b, TAN Wei1a, QI Yujiao1a

    (1.a College of Forestry; b. Guizhou Provincial Forest Resources and Environment Research Center, Guizhou University, Guiyang 550025, Guizhou, China; 2. Guizhou Singularity Forest Survey and Design Co., Anshun 561000, Guizhou, China; 3. Guizhou Academy Forestry, Guiyang 550005, Guizhou, China)

    Abstract:【Objective】Evaluate the accuracy of the leaf area index (LAI) of Pinus massoniana by optical instrumentation to provide technical support for accurately monitoring the LAI seasonal dynamics of P. massoniana forest.【Method】We evaluated the seasonal dynamics of LAI estimated by leaf-spreading surveys, the anisotropic growth equation method, the apomictic method, and optical instrumental methods (hemispherical photography and the LAI-2200 plant canopy analyser) in the P. massoniana forest by using the apomictic method measured values as representative of the true LAI values of P. massoniana forest, and we considered the sources of error of the LAI determined by the optical instrumental methods (clumping index (ΩE), woody-to-total area ratio (α), and needle-cluster ratio (γE)).【Result】The mean leaf lifespan of the major species within the P. massoniana forest ranged from 1.53 to 1.94, and SLA varied widely (99.15 cm2/g-278.33 cm2/g). The seasonal dynamics of LAI in P. massoniana forests determined by different methods showed a single-peak change pattern, and reached a peak in July. In each period, there was no significant difference between the value determined by litterfall method and allometry method. During the entire investigation period, the effective leaf area index (Le) measured by DHP was underestimated by 49%-64% and 47%-61%, respectively, compared with the values measured by litterfall method and allometry method. Le measured by LAI-2200 underestimates 25%-44% and 20%-42%, respectively. The accuracy of Le measured by DHP was obviously improved after corrected by automatic exposure, α, ΩE and γE. Compared with litterfall method and allometry method, the accuracy was average underestimated by 8% and 5%. Generally speaking, after corrected by α, ΩE and γE, the accuracy of LAI-2200 was greatly improved, compared with litterfall method and allometry method, it was underestimated by an average of 6% and 2%, respectively.【Conclusion】After considering the factors that affect the determination of LAI by optical instrument method, the accuracy of LAI measured by DHP and LAI-2200 methods had been significantly improved, indicating that both optical instrument methods can accurately determine the seasonal dynamics of LAI in P. massoniana forest and can provide a methodological reference for the study of seasonal dynamics of LAI in P. massoniana forests in different regions.

    Keywords: Pinus massoniana plantation; leaf area index; seasonal dynamics; litter collection; digital hemispherical photography; LAI-2200 plant canopy analyzer method

    葉面積指數(shù)(leaf area index,LAI)通常被定義為單位面積上總?cè)~片表面積的一半[1],它控制著植被的生理生態(tài)過程,如光合、呼吸、蒸騰、碳循環(huán)、降雨截留及物質(zhì)能量交換等[2-4],是生態(tài)和生物化學(xué)模型中的關(guān)鍵參數(shù)。作為植物群體和群落生長(zhǎng)分析的重要指標(biāo),LAI在農(nóng)業(yè)[5]、生態(tài)學(xué)[6]和環(huán)境科學(xué)[7]等領(lǐng)域得到廣泛研究。因此,準(zhǔn)確測(cè)定LAI是開展地球系統(tǒng)碳、氮、水循環(huán)模擬等研究的關(guān)鍵[8],也是評(píng)估森林生態(tài)系統(tǒng)對(duì)全球氣候變化響應(yīng)的基礎(chǔ)[9]。

    測(cè)定LAI的方法分為直接法和間接法[10-11]。直接法主要包括破壞性取樣法、異速生長(zhǎng)方程法、斜點(diǎn)樣方法和凋落物法等,這些方法比較經(jīng)典且相對(duì)成熟,所測(cè)結(jié)果準(zhǔn)確性高[12-13],常代表真實(shí)值校驗(yàn)間接法的測(cè)定值。破壞性取樣法適合于對(duì)農(nóng)田和草原生態(tài)系統(tǒng)進(jìn)行LAI測(cè)量[14-15],由于其破壞力大,很難獲得冠層結(jié)構(gòu)復(fù)雜的森林LAI[10,16]。異速生長(zhǎng)方程法雖然具有很低的破壞性,但是不能用于LAI季節(jié)動(dòng)態(tài)測(cè)定,且極易受林分特性的影響[13,17]。斜點(diǎn)樣方法不具有破壞性,但近年來的應(yīng)用卻不多,因?yàn)樵摲椒ㄍǔP枰銐虼蟮臉颖緮?shù)量才能保證測(cè)量精度,實(shí)施難度大[18]。凋落物法則是在落葉林中應(yīng)用較多的一種方法,即定期對(duì)凋落物進(jìn)行收集,獲取落葉季節(jié)LAI動(dòng)態(tài)以及葉面積最大時(shí)期的LAI值[10,19]。目前,學(xué)者們將凋落物法與常綠樹種的平均葉壽命相結(jié)合得到常綠林葉面積最大時(shí)期的LAI[20-22],這為常綠落葉混交林LAI季節(jié)變化的精準(zhǔn)測(cè)定提供了新的思路[18,22]。間接法(地面測(cè)定)主要指光學(xué)儀器法,其可以測(cè)定不同森林類型的LAI,具有方便快捷、無破壞性、不受時(shí)間限制、獲取數(shù)據(jù)量大等特點(diǎn)[19]。常用的光學(xué)儀器設(shè)備主要有:TRAC、LAI2200植物冠層分析儀(LAI-2200)、半球攝影圖像(digital hemispherical photography,DHP)和Sunfleck ceptometer等[11-12,23]。但光學(xué)儀器法也有其不足之處:1)將木質(zhì)部劃分為葉片而高估 LAI;2)忽視了冠層和針簇內(nèi)水平的集聚從而過低估計(jì)LAI[11,13,24]?;诖?,研究者們總結(jié)了影響光學(xué)儀器法測(cè)定精度的3種誤差源,分別是木質(zhì)部、冠層水平集聚及簇內(nèi)集聚[25-26],其中,簇內(nèi)集聚效應(yīng)僅存在針葉林內(nèi)。這3種誤差源分別用集聚指數(shù)(clumping index, ΩE)、木質(zhì)比例(woody-to total area ratio, α)及針簇比(needle-to shoot area ratio, γE)3個(gè)參數(shù)來表示[25-27]。此外,DHP測(cè)定精度還受曝光設(shè)置的影響[27-28]。因此,為獲得相對(duì)準(zhǔn)確的LAI,需要對(duì)光學(xué)儀器測(cè)定值進(jìn)行合理校正以提高其測(cè)定精度。

    近年來,學(xué)者們逐漸對(duì)不同森林類型LAI開展研究,如千煙洲針葉林[29]、廣西桉樹Eucalyptus robusta人工林[30]和青海云杉Picea crassifolia林[31]等,涉及的方法多是光學(xué)儀器法和異速生長(zhǎng)方程法。隨著研究的深入,凋落物法與光學(xué)儀器法相結(jié)合測(cè)定不同森林類型LAI的季節(jié)動(dòng)態(tài)得到廣泛應(yīng)用,如闊葉紅松Pinus koraiensis林[32]、帽兒山落葉闊葉混交林[33]以及黑龍江涼水國家級(jí)自然保護(hù)區(qū)紅松和興安落葉松Larix gmelinii人工林[27]等。馬尾松Pinus massoniana林是我國南方的主要森林類型之一,目前對(duì)該森林的研究主要集中在生物量空間特征[34]、病蟲害研究及預(yù)測(cè)模型構(gòu)建[35-36]、凋落物及土壤理化性質(zhì)[37-38]和生長(zhǎng)模型[39]等方面,而關(guān)于LAI季節(jié)動(dòng)態(tài)方面鮮有研究,尤其是將凋落物法和光學(xué)儀器法相結(jié)合測(cè)定LAI。因此,本研究以馬尾松林為研究對(duì)象,以展葉調(diào)查和凋落物法測(cè)定值為基礎(chǔ),對(duì)光學(xué)儀器法測(cè)定LAI進(jìn)行精度評(píng)估,分析LAI季節(jié)動(dòng)態(tài)變化特征,以期為科學(xué)經(jīng)營(yíng)和管理馬尾松林提供理論指導(dǎo)。

    1 材料與方法

    1.1 研究區(qū)概況與樣地設(shè)置

    研究區(qū)位于貴州省貴陽市花溪區(qū)(26°11′~26°34′N、106°27′~106°52′E),該區(qū)位于貴陽市中部,地勢(shì)呈東西高、中間低,海拔在1 400 m以上。氣候區(qū)劃屬于亞熱帶季風(fēng)氣候區(qū),冬無嚴(yán)寒,夏無酷熱,年均氣溫15.3 ℃,年降水量1 197~1 248 mm。研究區(qū)內(nèi)主要喬木樹種有馬尾松、油茶Camellia oleifera、猴樟Cinnamomum bodinieri、杉木Cunninghamia lanceolata和楓香樹Liquidambar formosana,同時(shí)還存在少量的化香樹Platycarya strobilacea、栲Castanopsis fargesii、樸樹Celtis sinensis和野漆Toxicodendron succedaneum。

    在該研究區(qū)馬尾松林設(shè)置1個(gè)40 m×50 m的固定樣地,并劃分為20個(gè)10 m×10 m的小樣方,對(duì)每個(gè)樣方內(nèi)胸徑≥1 cm的樹種進(jìn)行每木檢尺,記錄樹種名稱、胸徑、樹高和坐標(biāo)等因子,樣地描述性統(tǒng)計(jì)見表1。在小樣方中心位置布設(shè)1個(gè)凋落物收集器,距離地面高度約為0.5 m。該收集裝置由鐵絲(直徑8 mm)和尼龍網(wǎng)(1 mm的孔徑,0.5~0.6 m的深度)編織而成,網(wǎng)口為正方形,面積為1.0 m2。同時(shí),在每個(gè)收集器旁布設(shè)3根PVC管(長(zhǎng)約30 cm)以固定光學(xué)儀器的位置。本研究數(shù)據(jù)收集分為3個(gè)時(shí)期:1)2019年3—12月,利用光學(xué)儀器(DHP和LAI-2200)測(cè)定Le;2)2019年3—7月,進(jìn)行展葉調(diào)查;3)2018年9月—2020年8月,進(jìn)行凋落物的采集。

    1.2 用直接法測(cè)定的葉面積指數(shù)

    1.2.1 最大時(shí)期葉面積指數(shù)(LAIm)的測(cè)定

    本研究將各樹種不同時(shí)期凋落葉數(shù)據(jù)與比葉面積(specific leaf area,SLA)結(jié)合得到最大時(shí)期葉面積指數(shù)(LAIm)。SLA指葉片單位干質(zhì)量的葉面積[40]。落葉闊葉樹種的LAIm通過累加落葉季節(jié)凋落葉產(chǎn)生的LAI得出,針葉樹種的LAIm則通過測(cè)定1年內(nèi)(本研究為2018年9月—2019年8月)凋落葉產(chǎn)生的LAI再乘以各樹種的平均葉壽命得到[22]。

    闊葉樹種和針葉樹種使用不同的方法測(cè)定SLA。闊葉樹種:隨機(jī)獲取50個(gè)成熟時(shí)期的平整健康的葉片,隨后采用掃描法得到LAI,即葉片所占像素和像素大小的乘積。針葉樹種:采集大量針葉帶回室內(nèi)后使用體積替換法[11]測(cè)定其半表面積。馬尾松SLA的具體測(cè)定方法及公式見前期研究結(jié)果[41]。將處理過的闊葉和針葉放置于65 ℃烘箱烘干至恒質(zhì)量,隨后使用電子天平對(duì)樣本稱質(zhì)量(精確到0.001 g),根據(jù)SLA的定義計(jì)算出各主要樹種的SLA[22]。

    馬尾松平均葉壽命根據(jù)不同年齡凋落針葉產(chǎn)生的LAI進(jìn)行加權(quán)得到[11],本研究采用實(shí)地采集樣枝法獲取馬尾松針葉的存活率。具體步驟為:1)根據(jù)樣地馬尾松胸徑的統(tǒng)計(jì)結(jié)果,將馬尾松分為3個(gè)林層:主林層(DBH≥40 cm)、次林層(40>DBH≥20 cm)和被壓層(DBH<20 cm)。2)在每個(gè)林層中選擇一棵樣樹,將樣樹林冠分為上、中、下3層,每層取6個(gè)樣枝。樣枝上各小枝的年齡是根據(jù)其生物學(xué)特性來決定的,而小枝的年齡則是該小枝上馬尾松針葉的壽命[21]。3)采用最大針葉數(shù)量進(jìn)行規(guī)范化處理(最大針葉樹為1.0),對(duì)馬尾松不同年齡針葉存活率進(jìn)行統(tǒng)計(jì)。4)計(jì)算3個(gè)層次針葉的平均壽命,并根據(jù)各層次胸高斷面積比例加權(quán)得到整個(gè)林分內(nèi)馬尾松針葉的平均葉壽命。常綠闊葉樹種的平均葉壽命根據(jù)生活型得到[42]。

    1.2.2 生長(zhǎng)季節(jié)葉面積指數(shù)的測(cè)定

    利用展葉調(diào)查確定馬尾松林主要樹種葉片物候變化規(guī)律。本研究于2019年3—7月對(duì)林內(nèi)主要針闊葉樹種開展展葉調(diào)查。

    針葉樹種:選擇馬尾松樣樹3株,從每株樣樹選擇20個(gè)樣枝,在各調(diào)查時(shí)期需測(cè)定樣枝長(zhǎng)度、單位長(zhǎng)度樣枝上針葉數(shù)量和每個(gè)樣枝上5~10個(gè)針葉的平均長(zhǎng)度。馬尾松針葉形狀近似橢圓柱體,針葉先端尖或成漸尖的長(zhǎng)尖頭,因此可忽略針葉頂端的葉面積。馬尾松t時(shí)期各樣枝上的總?cè)~面積如公式(2)所示:

    闊葉樹種:在每個(gè)調(diào)查期內(nèi)測(cè)定樣枝上所有新生葉片(長(zhǎng)度大于0.5 cm的葉片)的葉長(zhǎng)(最長(zhǎng)處的葉長(zhǎng))和葉寬(最寬處的葉寬),再結(jié)合各樹種的葉面積調(diào)整系數(shù)(m)得到各樹種的葉面積[22]。本研究選取馬尾松林中主要闊葉樹種,每個(gè)樹種選3株樣樹,每株樣樹選擇1個(gè)健康的樣枝,并對(duì)健康樣枝進(jìn)行測(cè)定。本研究闊葉樹種葉面積計(jì)算公式參考劉志理[11]的研究。

    1.2.3 落葉季節(jié)葉面積指數(shù)的測(cè)定

    對(duì)于針葉和闊葉樹種,不同落葉階段的LAI均可按每一樹種的LAIm依次減去落葉季各個(gè)階段凋落葉產(chǎn)生的LAI計(jì)算得到,本研究參考劉志理等[22]的研究方法。

    1.2.4 異速生長(zhǎng)方程法測(cè)定的葉面積指數(shù)

    本研究中的異速生長(zhǎng)方程采用汪珍川等[43]的研究結(jié)果,構(gòu)建以胸徑為自變量、葉生物量為因變量的回歸方程,葉生物量通過SLA被轉(zhuǎn)化為各樹種LAIm,結(jié)合不同生長(zhǎng)期總LAI的增加比例及落葉季節(jié)的凋落物數(shù)據(jù),得出異速生長(zhǎng)方程法測(cè)定的LAI動(dòng)態(tài)。而馬尾松則采用以下方法:在所劃分的3層林層中,各林層分別選取3棵馬尾松,即優(yōu)勢(shì)木、標(biāo)準(zhǔn)木及被壓木,全部伐倒后測(cè)定葉生物量,并建立胸徑與葉生物量的回歸關(guān)系方程。

    1.3 光學(xué)儀器法測(cè)定的葉面積指數(shù)

    DHP通過Winscanopy2006冠層分析儀采集圖像,該裝置主要由數(shù)碼相機(jī)(Canon EOS 6D)和180°魚眼鏡頭組成,使用三腳架將數(shù)碼相機(jī)固定在距離水平地面1.3 m處。為避免陽光直射對(duì)圖像采集的影響,選擇在陰天或近日出和日落前后采集圖像,并將相機(jī)設(shè)置為自動(dòng)曝光模式,圖像存為JPEG圖像格式。采用LAI-2200與DHP獲取數(shù)據(jù)時(shí)應(yīng)保持時(shí)間(2019年3—12月)、地點(diǎn)和高度一致,并保持光學(xué)感應(yīng)探頭處于水平狀態(tài)。

    1.4 數(shù)據(jù)處理

    采用DHP軟件處理半球攝影圖像,處理時(shí)提取30°~60°天頂角范圍內(nèi)的Le[44]。對(duì)DHP測(cè)定的Le采用自動(dòng)曝光校正公式進(jìn)行校正[28]。使用Photoshop7.0軟件計(jì)算木質(zhì)比例α [11],用DHPTRACWin軟件獲取ΩE[45],計(jì)算時(shí)選取的天頂角為40°~45°。采用FV2200軟件分析LAI-2200數(shù)據(jù),處理時(shí)提取3~4環(huán)內(nèi)的Le。對(duì)馬尾松林內(nèi)展葉調(diào)查法和凋落物法測(cè)定的LAI使用配對(duì)t檢驗(yàn)進(jìn)行差異性顯著檢驗(yàn)(α=0.05)。本研究所有數(shù)據(jù)處理在SPSS22.0軟件中完成,作圖在Origin 2018軟件中完成。

    2 結(jié)果與分析

    2.1 主要樹種的平均葉壽命和比葉面積

    馬尾松林內(nèi)不同樹種的SLA差異較大(表2)。闊葉樹種中楓香樹的SLA最大,為278.33 cm2/g;猴樟的最小,為134.83 cm2/g。馬尾松的SLA為99.15 cm2/g,馬尾松的平均葉壽命為1.53 a。

    2.2 主要闊葉樹種的葉面積調(diào)整系數(shù)

    葉形不同導(dǎo)致不同樹種的葉面積調(diào)整系數(shù)存在顯著差異(表3)。本研究中3種主要闊葉樹種的葉面積調(diào)整系數(shù)為0.44~0.68,平均值為0.60(標(biāo)準(zhǔn)差為0.04)。葉面積調(diào)整系數(shù)由葉形決定,心形或橢圓形葉片的長(zhǎng)寬乘積更接近真實(shí)葉面積,即其葉形的葉面積調(diào)整系數(shù)值高于其他形狀葉片,如楓香樹(掌狀形葉片)的葉面積調(diào)整系數(shù)為0.44,而油茶(橢圓形葉片)的葉面積調(diào)整系數(shù)為0.68。

    2.3 針簇比(γE)

    3個(gè)林層馬尾松針簇比為(1.41±0.31)~(1.94±0.63)(平均值±標(biāo)準(zhǔn)差),各林層γE平均值從大到小排序?yàn)椋褐髁謱樱?.84±0.53)>次林層(1.75±0.52)>被壓層(1.52±0.43)(表4)。在各林層內(nèi),樹冠位置越高,γE值越大。

    2.4 集聚指數(shù)(ΩE)及木質(zhì)部分所占比例(α)的季節(jié)動(dòng)態(tài)

    馬尾松林冠層水平上的集聚效應(yīng)表現(xiàn)出較小的季節(jié)波動(dòng)(圖1)。ΩE最大值出現(xiàn)在3月(ΩE=0.93),最小值出現(xiàn)在11月(ΩE=0.89),各月份ΩE均值為0.91。木質(zhì)比例α為0.05~0.09,平均值為0.07。α值呈先減小后增加的季節(jié)變化模式,隨著葉片生長(zhǎng)逐漸減小,隨葉片凋落逐漸增加。

    2.5 馬尾松林主要樹種LAI的季節(jié)變化

    總體來看,針葉樹種和闊葉樹種葉面積增加比例不同(圖2)。馬尾松在3月份出現(xiàn)展葉高峰,猴樟和油茶在5月份出現(xiàn)展葉高峰,但3種樹種在6月初均已完成全年總?cè)~面積90%以上的生長(zhǎng)量。各樹種的LAI隨季節(jié)變化呈先升高后降低的趨勢(shì),其峰值均出現(xiàn)在7月份(圖3),表現(xiàn)為馬尾松(7.33)>猴樟(0.63)>油茶(0.61)。相對(duì)闊葉樹種,整個(gè)調(diào)查期內(nèi),馬尾松的LAI隨季節(jié)變化差異顯著,展葉期3、4和5月份LAI隨葉片的生長(zhǎng)逐漸增大,6和7月份展葉基本完成,其LAI波動(dòng)較小,進(jìn)入落葉季節(jié),葉片凋落,其LAI逐漸減小。而猴樟和油茶LAI在整個(gè)季節(jié)隨季節(jié)變化波動(dòng)不存在顯著差異。

    2.6 用直接法和間接法測(cè)定的LAI的比較

    由圖4可知,使用凋落物法和異速生長(zhǎng)方程法測(cè)定的馬尾松LAI在各季節(jié)內(nèi)的差異不顯著,且都表現(xiàn)出隨著季節(jié)變化呈現(xiàn)出單峰形的變化模式,兩種方法所測(cè)LAI值均在7月份達(dá)到最高值,分別為8.60和8.27。其中,在整個(gè)調(diào)查期內(nèi)凋落物法測(cè)定值均高于異速生長(zhǎng)方程法測(cè)量值。

    不同方法測(cè)定常綠針葉林的LAI隨季節(jié)變化差異明顯,且均為單峰形變化,在7月達(dá)到峰值(圖5)。生長(zhǎng)季節(jié)的LAI隨新展出葉片的生長(zhǎng)而增大,落葉季節(jié)的LAI隨葉片凋落逐漸減小。總體來看,光學(xué)儀器法測(cè)定的Le在各時(shí)期均小于凋落物法和異速生長(zhǎng)方程法的測(cè)定值。真實(shí)的LAI(凋落物法)為5.13~8.60,異速生長(zhǎng)方程法的LAI為5.05~8.27。DHP測(cè)定的Le為2.03~3.58,LAI-2200測(cè)定的Le為3.78~5.51。整個(gè)調(diào)查期內(nèi),DHP測(cè)定的Le比真實(shí)的LAI低估了49%~64%,平均低估57%;LAI-2200則低估了25%~44%,平均低估32%(圖6)。相對(duì)而言,DHP測(cè)定的Le比異速生長(zhǎng)方程法的LAI低估47%~61%,平均低估56%;LAI-2200測(cè)定的Le低估20%~42%,平均低估30%(圖6)。

    整個(gè)季節(jié)內(nèi),校正后的兩種光學(xué)儀器測(cè)定值精度均明顯提高。其中,校正后的DHP測(cè)定值比真實(shí)LAI低估4%~16%,平均低估8%,其測(cè)定精度高達(dá)84%;校正后的LAI-2200測(cè)定值低估2%~6%,平均低估4%。相對(duì)異速生長(zhǎng)方程法,校正后的DHP測(cè)定值在整個(gè)調(diào)查期內(nèi)低估2%~12%,表明DHP經(jīng)誤差源校正后,其精度高于88%;而LAI-2200測(cè)定的Le經(jīng)校正后在3、5、7和10月份比異速生長(zhǎng)方程法測(cè)定的LAI分別高估了1%、2%、4%及6%(圖6),其他時(shí)期比異速生長(zhǎng)方程法測(cè)定的LAI平均低估2%。

    3 討 論

    凋落物法和異速生長(zhǎng)方程法是測(cè)定不同森林類型LAI的常用直接法,但利用兩種直接法測(cè)定并對(duì)比馬尾松林LAI季節(jié)動(dòng)態(tài)的研究較少。本研究對(duì)比分析兩種直接法測(cè)定馬尾松林LAI的季節(jié)動(dòng)態(tài),結(jié)果顯示,各調(diào)查期內(nèi),兩種方法的測(cè)定值均不存在顯著差異,表明異速生長(zhǎng)方程法的構(gòu)建雖受到林分特性如地域、樹種及林齡等的限制[46],但只要合理運(yùn)用,仍有較廣泛的適用性。

    各樹種的SLA和常綠樹種的葉壽命是精準(zhǔn)測(cè)定森林LAI的關(guān)鍵因子[47]。有研究發(fā)現(xiàn),闊葉樹種和針葉樹種之間的SLA差異較大[48]。本研究中,闊葉樹種的平均SLA是針葉樹種的1.9倍。有學(xué)者在對(duì)黑龍江省涼水國家級(jí)自然保護(hù)區(qū)內(nèi)針闊葉樹種LAI研究時(shí)同樣發(fā)現(xiàn),闊葉樹種的平均SLA是針葉樹種的3.6倍,針闊葉樹種SLA差異顯著[22]。此外,同一樹種的SLA隨季節(jié)變化也存在一定差異[49],Nouvellon等[50]研究發(fā)現(xiàn)桉樹Eucalyptus林平均SLA在1年內(nèi)的變異約為20%。本研究只考慮了各樹種葉面積最大時(shí)期的SLA,為更準(zhǔn)確地利用凋落物法測(cè)定馬尾林LAI的季節(jié)變化,后續(xù)研究中將進(jìn)一步探討。本研究通過采集樣枝法測(cè)定了馬尾松的葉壽命,其平均葉壽命為1.53,劉志理等[22]的研究結(jié)果顯示,紅松、臭冷杉Abies nephrolepis及云杉Picea asperata的平均葉壽命分別為3.08、3.68及3.91 a。周自宗等[42]的研究表明,葉片壽命與葉片形狀和結(jié)構(gòu)、葉片出生率和死亡率及生境異質(zhì)性等相關(guān),這可能是導(dǎo)致馬尾松平均葉壽命與紅松、臭冷杉及云杉差異顯著的原因。

    各時(shí)期內(nèi),DHP測(cè)定的有效葉面積指數(shù)均明顯低于直接法(凋落物法和異速生長(zhǎng)方程法)。這與其他學(xué)者得出的結(jié)果類似,如王寶琦等[51]分別使用光學(xué)儀器法和凋落物法對(duì)紅松人工林內(nèi)有效葉面積指數(shù)進(jìn)行測(cè)定,結(jié)果發(fā)現(xiàn),DBH方法比凋落物法最大時(shí)期低估77.06%,平均低估75.5%。馬澤清等[52]利用DHP測(cè)定了千煙洲人工林中馬尾松林、杉木林和濕地松Pinus elliottii林內(nèi)的有效葉面積指數(shù),結(jié)果顯示有效葉面積指數(shù)的降低幅度分別為61%、57%和68%;謝賢勝等[30]對(duì)廣西桉樹人工林葉面積指數(shù)進(jìn)行估測(cè)時(shí)也得到相似結(jié)論。與DHP法相比,LAI-2200測(cè)定值的低估程度要小得多。如郝佳等[53]對(duì)華北落葉松Larix principis-rupprechtii林進(jìn)行測(cè)定,發(fā)現(xiàn)有效葉面積指數(shù)平均值低于凋落物法的29%;劉志理等[27]測(cè)定常綠針葉林內(nèi)的有效葉面積指數(shù)比凋落物法低估19%~35%;周明等[26]測(cè)定興安落葉松人工林的有效葉面積指數(shù)比凋落物法低估19%~32%,比異速生長(zhǎng)方程法低估18%~29%。

    光學(xué)儀器法能省時(shí)、方便且省力,在測(cè)定不同類型森林的LAI季節(jié)動(dòng)態(tài)變化方面被廣泛應(yīng)用。與直接法相比,光學(xué)儀器法存在低估的現(xiàn)象,因此需要綜合考慮影響光學(xué)儀器法測(cè)定LAI的誤差源,進(jìn)而獲得更加準(zhǔn)確的LAI[47]。本研究采用PS法量化馬尾松林α值的季節(jié)變化,結(jié)果表明α值呈單峰形變化模式。在生長(zhǎng)期,各樹種葉片生長(zhǎng)葉面積增大,木質(zhì)部被新生葉片遮擋進(jìn)而減小了其對(duì)LAI的貢獻(xiàn)率,而隨著葉片的逐漸凋落,葉面積減小,被遮擋的木質(zhì)部重新暴露出來從而加大了對(duì)LAI的貢獻(xiàn)率。對(duì)冠層水平上的集聚效應(yīng)進(jìn)行季節(jié)性量化,結(jié)果顯示集聚指數(shù)為0.89~0.93,與其他結(jié)果相符合,如Chen等[44]量化不同森林類型的集聚指數(shù)為0.88~0.96。對(duì)于馬尾松林,還存在葉簇內(nèi)的集聚現(xiàn)象,本研究采用實(shí)地樣枝法測(cè)定馬尾松的針簇比平均值為1.70,與其他學(xué)者得到的結(jié)論存在異同。如王寶琦等[51]對(duì)黑龍江省涼水國家級(jí)自然保護(hù)區(qū)內(nèi)紅松人工林的針簇比測(cè)定結(jié)果為1.73,而劉志理等[27]測(cè)定的興安落葉松人工林內(nèi)針簇比結(jié)果為1.00~0.30,可能與樹種生物學(xué)特性、生境等有關(guān)。

    對(duì)DHP和LAI-2200測(cè)定的有效葉面積指數(shù)進(jìn)行校正后,其精度顯著提高。與凋落物法測(cè)定值相比,在整個(gè)調(diào)查期內(nèi),DHP低估了4%~16%,LAI-2200測(cè)定的有效葉面積指數(shù)經(jīng)校正后比凋落物法低估2%~6%。與異速生長(zhǎng)方程法測(cè)定值相比,DHP測(cè)定精度高達(dá)88%,LAI-2200測(cè)定有效葉面積指數(shù)經(jīng)校正后,其精度也明顯提高,比異速生長(zhǎng)方程法平均低估2%。其他研究結(jié)果也表明光學(xué)儀器法測(cè)定的有效葉面積指經(jīng)過上述誤差源校正后,其精度大幅度提高[22,26]。

    4 結(jié) 論

    本研究采用展葉調(diào)查、異速生長(zhǎng)方程法、凋落物法、光學(xué)儀器法(半球攝影法和LAI-2200植物冠層分析儀)探討了馬尾松林LAI的季節(jié)動(dòng)態(tài),并評(píng)估了2種光學(xué)儀器法測(cè)定值的精度。結(jié)果發(fā)現(xiàn),不同方法測(cè)定的LAI存在差異,但均為單峰形變化,且在7月達(dá)到峰值。在整個(gè)調(diào)查期內(nèi),光學(xué)儀器法測(cè)定的Le在各時(shí)期均小于凋落物法和異速生長(zhǎng)方程法的測(cè)定值,平均低估30%~57%。校正后的兩種光學(xué)儀器測(cè)定值精度均明顯提高,校正后的DHP測(cè)定值比真實(shí)LAI低估4%~16%,平均低估8%,其測(cè)定精度高達(dá)84%;校正后的LAI-2200測(cè)定值低估2%~6%,平均低估4%。相對(duì)異速生長(zhǎng)方程法,校正后的DHP測(cè)定值在整個(gè)調(diào)查期內(nèi)低估2%~12%,表明DHP經(jīng)以上誤差源校正后,其精度高于88%。研究結(jié)果為光學(xué)儀器法在未來相對(duì)準(zhǔn)確測(cè)定相似林型LAI季節(jié)動(dòng)態(tài)提供科學(xué)的理論依據(jù)。本研究沒有考慮不同森林類型條件下樹種LAI的季節(jié)動(dòng)態(tài)變化,后期可結(jié)合實(shí)際情況對(duì)不同森林類型進(jìn)行研究。

    參考文獻(xiàn):

    [1] CHEN J M, BLACK T A. Defining leaf area index for non-flat leaves[J]. Plant, Cell and Environment,1992,15(4):421-429.

    [2] PARKER G G. Tamm review: Leaf area index (LAI) is both a determinant and a consequence of important processes in vegetation canopies[J]. Forest Ecology and Management,2020,477:118496.

    [3] MAJASALMI T, RAUTIAINEN M, STENBERG P. Modeled and measured fPAR in a boreal forest: Validation and application of a new model[J]. Agricultural and Forest Meteorology, 2014,189:118-124.

    [4] WU J, ALBERT L P, LOPES A P, et al. Leaf development and demography explain photosynthetic seasonality in Amazon evergreen forests[J]. Science,2016,351:972-976.

    [5] 肖強(qiáng)智,徐維新,代娜,等.高寒牧草葉面積的冬季衰減過程與估算[J].草地學(xué)報(bào),2023,31(8):2436-2445. XIAO Q Z, XU W X, DAI N, et al. The temporal pattern of leaf area index in withered procedure and its estimation in winter in Qinghai-Tibet Plateau[J]. Acta Agrestia Sinica,2023,31(8): 2436-2445.

    [6] 陳卓,范文義,于穎,等.運(yùn)用機(jī)載激光雷達(dá)數(shù)據(jù)估算不同森林類型葉面積指數(shù)[J].東北林業(yè)大學(xué)學(xué)報(bào),2020,48(11):39-45. CHEN Z, FAN W Y, YU Y, et al. Estimating leaf area index of different forest types using airborne laser scanner[J]. Journal of Northeast Forestry University,2020,48(11):39-45.

    [7] 孫越,顧祝軍,李棟梁.無人機(jī)與衛(wèi)星影像的葉面積指數(shù)遙感反演研究[J].測(cè)繪科學(xué),2021,46(2):106-112,145. SUN Y, GU Z J, LI D L. Study on remote sensing retrieval of leaf area index based on unmanned aerial vehicle and satellite imagery[J]. Science of Surveying and Mapping,2021,46(2): 106-112,145.

    [8] 段利民,李瑋,羅艷云,等.錫林河流域植被葉面積指數(shù)時(shí)空變化特征及其對(duì)地形因子的響應(yīng)[J].水土保持研究, 2019,26(3):224-231. DUAN L M, LI W, LUO Y Y, et al. Spatiotemporal variations of vegetation leaf area index and its response to topographical factors in Xilin river basin[J]. Research of Soil and Water Conservation, 2019,26(3):224-231.

    [9] BEQUET R, KINT V, CAMPIOLI M, et al. Influence of stand, site and meteorological variables on the maximum leaf area index of beech, oak and Scots pine[J]. European Journal of Forest Research,2012,131(2):283-295.

    [10] BRéDA N J J. Ground-based measurements of leaf area index: a review of methods, instruments and current controversies[J]. Journal of Experimental Botany,2003,54(392):2403-2417.

    [11] 劉志理.東北典型森林葉面積指數(shù)的時(shí)空動(dòng)態(tài)[D].哈爾濱:東北林業(yè)大學(xué),2015. LIU Z L. Spatial and temporal dynamics of leaf area index in typical forests in northeastern China[D]. Harbin: Northeast Forestry University,2015.

    [12] 鄒杰,閻廣建.森林冠層地面葉面積指數(shù)光學(xué)測(cè)量方法研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2010,21(11):2971-2979. ZOU J, YAN G J. Optical methods for in situ measuring leaf area index of forest canopy: a review[J]. Chinese Journal of Applied Ecology,2010,21(11):2971-2979.

    [13] CHEN J M, RICH P M, GOWER S T, et al. Leaf area index of boreal forests: Theory, techniques, and measurements[J]. Journal of Geophysical Research,1997,102:29429-29443.

    [14] LIU J G, PATTEY E, ADMIRAL S. Assessment of in situ crop LAI measurement using unidirectional view digital photography[J]. Agricultural and Forest Meteorology,2013,169: 25-34.

    [15] 王猛,李貴才,王軍邦.典型草原通量塔通量貢獻(xiàn)區(qū)地上生物量和葉面積指數(shù)的時(shí)空變異[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(3): 637-643. WANG M, LI G C, WANG J B. Spatiotemporal variations of aboveground biomass and leaf area index of typical grassland in tower flux footprint[J]. Chinese Journal of Applied Ecology,2011,22(3):637-643.

    [16] 王曉紅,申方圓,周鑫,等.撫遠(yuǎn)地區(qū)三種林型葉面積指數(shù)的估測(cè)[J].溫帶林業(yè)研究,2018,1(3):36-41. WANG X H, SHEN F Y, ZHOU X, et al. Estimation of leaf area index of three forest types in Fuyuan study area[J]. Journal of Temperate Forestry Research,2018,1(3):36-41.

    [17] GOWER S T, REICH P B, SON Y. Canopy dynamics and aboveground production of five tree species with different leaf longevities[J]. Tree Physiology,1993,12:327-345.

    [18] 向洪波,郭志華,趙占輕,等.不同空間尺度森林葉面積指數(shù)的估算方法[J].林業(yè)科學(xué),2009,45(6):139-144. XIANG H B, GUO Z H, ZHAO Z Q, et al. Estimating method of forest leaf area index on different space scales[J]. Scientia Silvae Sinicae,2009,45(6):139-144.

    [19] 劉志理,金光澤.小興安嶺三種林型葉面積指數(shù)的估測(cè)[J].應(yīng)用生態(tài)學(xué)報(bào),2012,23(9):2437-2444. LIU Z L, JIN G Z. Estimation of leaf area index of three forest types in Xiaoxing’an Mountains of Northeast China[J]. Chinese Journal of Applied Ecology,2012,23(9):2437-2444.

    [20] 蘇宏新,白帆,李廣起.3類典型溫帶山地森林的葉面積指數(shù)的季節(jié)動(dòng)態(tài):多種監(jiān)測(cè)方法比較[J].植物生態(tài)學(xué)報(bào),2012,36(3): 231-242. SU H X, BAI F, LI G Q. Seasonal dynamics in leaf area index in three typical temperate montane forests of China: a comparison of multi-observation methods[J]. Chinese Journal of Plant Ecology, 2012,36(3):231-242.

    [21] REICH P B, FRELICH L E, VOLDSETH R A, et al. Understorey diversity in southern boreal forests is regulated by productivity and its indirect impacts on resource availability and heterogeneity[J]. Journal of Ecology,2012,100(2):539-545.

    [22] 劉志理,金光澤,周明.利用直接法和間接法測(cè)定針闊混交林葉面積指數(shù)的季節(jié)動(dòng)態(tài)[J].植物生態(tài)學(xué)報(bào),2014,38(8): 843-856. LIU Z L, JIN G Z, ZHOU M. Measuring seasonal dynamics of leaf area index in a mixed conifer-broadleaved forest with direct and indirect methods[J]. Chinese Journal of Plant Ecology,2014, 38(8):843-856.

    [23] FANG H L, YE Y C, LIU W W, et al. Continuous estimation of canopy leaf area index (LAI) and clumping index over broadleaf crop fields: an investigation of the PASTIS-57 instrument and smartphone applications[J]. Agricultural and Forest Meteorology, 2018,253-254: 48-61.

    [24] JONCKHEERE I, FLECK S, NACKAERTS K, et al. Review of methods for in situ leaf area index determination: Part I. Theories, sensors and hemispherical photography[J]. Agricultural and Forest Meteorology,2004,121(1-2):19-35.

    [25] CHEN J M. Optically-based methods for measuring seasonal variation of leaf area index in boreal conifer stands[J]. Agricultural and Forest Meteorology,1996,80:135-163.

    [26] 周明,劉志理,金光澤.利用3種校正方案提高間接法測(cè)定興安落葉松人工林葉面積指數(shù)的精度[J].植物生態(tài)學(xué)報(bào), 2016,40(6):574-584. ZHOU M, LIU Z L, JIN G Z. Improving the accuracy of indirect methods in estimating leaf area index using three correction schemes in a Larix gmelinii plantation[J]. Chinese Journal of Plant Ecology,2016,40(6):574-584.

    [27] 劉志理,金光澤.光學(xué)儀器法測(cè)定葉面積指數(shù)季節(jié)變化的誤差分析[J].林業(yè)科學(xué),2016,52(9):11-21. LIU Z L, JIN G Z. Bias analysis of seasonal changes of leaf area index derived from optical methods[J]. Scientia Silvae Sinicae, 2016,52(9):11-21.

    [28] ZHANG Y Q, CHEN J M, MILLER J R. Determining digital hemispherical photograph exposure for leaf area index estimation[J]. Agricultural and Forest Meteorology,2005,133: 166-181.

    [29] 李軒然,劉琪璟,蔡哲,等.千煙洲針葉林的比葉面積及葉面積指數(shù)[J].植物生態(tài)學(xué)報(bào),2007,31(1):93-101. LI X R, LIU Q J, CAI Z, et al. Specific leaf area and leaf area index of conifer plantations in Qianyanzhou station of subtropical China[J]. Journal of Plant Ecology,2007,31(1):93-101.

    [30] 謝賢勝,盧峰,楊元征,等.廣西桉樹人工林葉面積指數(shù)的估測(cè)及其校正[J].生態(tài)學(xué)報(bào),2022,42(11):4451-4462. XIE X S, LU F, YANG Y Z, et al. Leaf area index estimation and correction in Eucalyptus plantations in Guangxi Autonomous Region[J]. Acta Ecologica Sinica,2022,42(11):4451-4462.

    [31] 趙傳燕,沈衛(wèi)華,彭煥華.祁連山區(qū)青海云杉林冠層葉面積指數(shù)的反演方法[J].植物生態(tài)學(xué)報(bào),2009,33(5):860-869. ZHAO C Y, SHEN W H, PENG H H. Methods for determining canopy leaf area index of Picea crassifolia forest in Qilian mountains, China[J]. Chinese Journal of Plant Ecology,2009,33(5):860-869.

    [32] QI Y J, JIN G Z, LIU Z L. Optical and litter collection methods for measuring leaf area index in an old-growth temperate forest in northeastern China[J]. Journal of Forest Research,2013,18: 430-439.

    [33] LIU Z L, WANG C K, CHEN J M, et al. Empirical models for tracing seasonal changes in leaf area index in deciduous broadleaf forests by digital hemispherical photography[J]. Forest Ecology and Management,2015,351:67-77.

    [34] 丘麗萍,張帆,林嘉源,等.南方典型紅壤區(qū)馬尾松生物量空間特征及驅(qū)動(dòng)力研究[J].西南林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)), 2024,44(1):106-115. QIU L P, ZHANG F, LIN J Y, Et al. Spatial characteristics and driving forces of biomass of Pinus massoniana forest in a typical red soil region of southern China [J]. Journal of Southwest Forestry University (Natural Science),2024,44(1):106-115.

    [35] 蘆偉,談家金.松材線蟲侵染對(duì)馬尾松莖干內(nèi)生細(xì)菌群落結(jié)構(gòu)的影響[J].東北林業(yè)大學(xué)學(xué)報(bào),2022,50(1):105-110,122. LU W, TAN J J. Effects of pine wood nematode infection on the endophytic bacterial community structure in the stems of Pinus massoniana[J]. Journal of Northeast Forestry University,2022,50(1): 105-110,122.

    [36] 田開慧,陳怡帆,周宏威.湖南湘西州馬尾松毛蟲和松材線蟲病發(fā)生非線性建模與預(yù)測(cè)[J].森林工程,2022,38(6):38-44. TIAN K H, CHEN Y F, ZHOU H W. Prediction of occurrence trend of Dendrolimus punctatus and pine wilt disease in Xiangxi prefecture[J]. Forest Engineering,2022,38(6):38-44.

    [37] 李鵬,陳璇,楊章旗,等.不同密度馬尾松人工林枯落物輸入對(duì)土壤理化性質(zhì)的影響[J].水土保持學(xué)報(bào),2022,36(2): 368-377. LI P, CHEN X, YANG Z Q, et al. Effect of litter inputs on soil physical and chemical properties of Pinus massoniana plantation with different densities[J]. Journal of Soil and Water Conservation, 2022,36(2):368-377.

    [38] 劉平,李鵬,楊章旗,等.廣西不同林齡和區(qū)域馬尾松人工林的土壤C、N、P化學(xué)計(jì)量特征[J].廣西科學(xué),2022,29(1):192-200. LIU P, LI P, YANG Z Q, et al. Stoichiometric characteristics of soil C, N and P in Pinus massoniana plantations of different ages and regions in Guangxi[J]. Guangxi Science,2022,29(1):192-200.

    [39] 陳哲夫,肖化順,龍時(shí)勝.基于混合效應(yīng)的湖南馬尾松次生林單木生長(zhǎng)模型[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2021,41(1):100-108.CHEN Z F, XIAO H S, LONG S S. Growth model for individual tree of secondary Pinus massoniana forest in Hunan province based on mixed effect[J]. Journal of Central South University of Forestry Technology,2021,41(1):100-108.

    [40] MACFARLANE C, HOFFMAN M, EAMUS D, et al. Estimation of leaf area index in Eucalypt forest using digital photography[J]. Agricultural and Forest Meteorology,2007,143(3-4):176-188.

    [41] 羅光浪,何世強(qiáng),譚偉,等.不同林齡馬尾松人工林葉面積指數(shù)的測(cè)定[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2022,42(2):55-64. LUO G L, HE S Q, TAN W, et al. Determination of leaf area index in Pinus massoniana plantations of different ages[J]. Journal of Central South University of Forestry Technology,2022,42(2):55-64.

    [42] 周自宗,袁莉,王震洪.貴陽市常綠樹種葉子壽命的研究[J].生態(tài)科學(xué),2008,27(3):148-153. ZHOU Z Z, YUAN L, WANG Z H. Leaf longevity of main evergreen tree species of Guiyang city[J]. Ecological Science, 2008,27(3):148-153.

    [43] 汪珍川,杜虎,宋同清,等.廣西主要樹種(組)異速生長(zhǎng)模型及森林生物量特征[J].生態(tài)學(xué)報(bào),2015,35(13):4462-4472. WANG Z C, DU H, SONG T Q, et al. Allometric models of major tree species and forest biomass in Guangxi[J]. Acta Ecologica Sinica,2015,35(13):4462-4472.

    [44] CHEN J M, GOVIND A, SONNENTAG O, et al. Leaf area index measurements at Fluxnet-Canada forest sites[J]. Agricultural and Forest Meteorology,2006,140(1/4):257-268.

    [45] LEBLANC S G, CHEN J M, FERNANDES R, et al. Methodology comparison for canopy structure parameters extraction from digital hemispherical photography in boreal forests[J]. Agricultural and Forest Meteorology,2005,129:187-207.

    [46] CHEN J M, CIHLAR J. Plant canopy gap-size analysis theory for improving optical measurements of leaf-area index[J]. Applied Optics,1995,34(27):6211-6222.

    [47] 羅光浪.貴州省三種典型森林類型葉面積指數(shù)的測(cè)定研究[D].貴陽:貴州大學(xué),2022. LUO G L. Study on the determination of leaf area index of three typical forest types in Guizhou province[D]. Guiyang: Guizhou University,2022.

    [48] LONG W X, ZANG R G, SCHAMP B S, et al. Within-and among-species variation in specific leaf area drive community assembly in a tropical cloud forest[J]. Oecologia,2011,167(4): 1103-1113.

    [49] ISHIHARA M I, HIURA T. Modeling leaf area index from litter collection and tree data in a deciduous broadleaf forest[J]. Agricultural and Forest Meteorology,2011,151:1016-1022.

    [50] NOUVELLON Y, LACLAU J P, EPRON D, et al. Withinstand and seasonal variations of specific leaf area in a clonal Eucalyptus plantation in the Republic of Congo[J]. Forest Ecology and Management,2010,259:1796-1807.

    [51] 王寶琦,劉志理,戚玉嬌,等.利用不同方法測(cè)定紅松人工林葉面積指數(shù)的季節(jié)動(dòng)態(tài)[J].生態(tài)學(xué)報(bào),2014,34(8):1956-1964. WANG B Q, LIU Z L, QI Y J, et al. Seasonal dynamics of leaf area index using different methods in the Korean pine plantation[J]. Acta Ecologica Sinica,2014,34(8):1956-1964.

    [52] 馬澤清,劉琪璟,曾慧卿,等.南方人工林葉面積指數(shù)的攝影測(cè)量[J].生態(tài)學(xué)報(bào),2008,28(5):1971-1980. MA Z Q, LIU Q J, ZENG H Q, et al. Estimation of leaf area index of planted forests in subtropical China by photogrammetry[J]. Acta Ecologica Sinica,2008,28(5):1971-1980.

    [53] 郝佳,熊偉,王彥輝,等.華北落葉松人工林葉面積指數(shù)實(shí)測(cè)值與冠層分析儀讀數(shù)值的比較和動(dòng)態(tài)校正[J].林業(yè)科學(xué)研究, 2012,25(2):231-235. HAO J, XIONG W, WANG Y H, et al. The comparison and dynamic calibration between the LAI values of a Larix principisrupprechtii plantation determined by canopy scanner and litterfall collection[J]. Forest Research,2012,25(2):231-235.

    [本文編校:謝榮秀]

    猜你喜歡
    生長(zhǎng)研究
    FMS與YBT相關(guān)性的實(shí)證研究
    2020年國內(nèi)翻譯研究述評(píng)
    遼代千人邑研究述論
    碗蓮生長(zhǎng)記
    小讀者(2021年2期)2021-03-29 05:03:48
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    共享出行不再“野蠻生長(zhǎng)”
    生長(zhǎng)在哪里的啟示
    EMA伺服控制系統(tǒng)研究
    野蠻生長(zhǎng)
    NBA特刊(2018年21期)2018-11-24 02:48:04
    生長(zhǎng)
    文苑(2018年22期)2018-11-19 02:54:14
    国产xxxxx性猛交| 新久久久久国产一级毛片| 巨乳人妻的诱惑在线观看| 人人妻人人澡人人爽人人夜夜| 久久久久国内视频| svipshipincom国产片| 久久精品国产综合久久久| 国产老妇伦熟女老妇高清| 999久久久国产精品视频| 亚洲精品一二三| 99国产精品一区二区蜜桃av | 精品国内亚洲2022精品成人 | 亚洲精品在线美女| 中文字幕人妻丝袜一区二区| 十八禁人妻一区二区| 少妇的丰满在线观看| 亚洲欧美一区二区三区久久| 亚洲成人免费av在线播放| 一二三四社区在线视频社区8| 1024视频免费在线观看| 国产在线精品亚洲第一网站| 久久婷婷成人综合色麻豆| 久久av网站| 母亲3免费完整高清在线观看| 日本一区二区免费在线视频| 啦啦啦视频在线资源免费观看| 国产视频一区二区在线看| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 国产男女超爽视频在线观看| 美女主播在线视频| 亚洲av成人不卡在线观看播放网| 亚洲精品粉嫩美女一区| 成人黄色视频免费在线看| 久久国产精品人妻蜜桃| 999精品在线视频| 美女高潮喷水抽搐中文字幕| 欧美日韩福利视频一区二区| 黑人操中国人逼视频| 男人舔女人的私密视频| 久久国产精品大桥未久av| 99久久国产精品久久久| 亚洲专区字幕在线| 波多野结衣av一区二区av| 成在线人永久免费视频| 亚洲第一欧美日韩一区二区三区 | 人妻一区二区av| 久久人妻福利社区极品人妻图片| 99久久精品国产亚洲精品| 午夜福利视频精品| 黑人巨大精品欧美一区二区蜜桃| 91av网站免费观看| 操出白浆在线播放| 免费久久久久久久精品成人欧美视频| 中文字幕av电影在线播放| 免费在线观看影片大全网站| 久久精品人人爽人人爽视色| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| 一二三四在线观看免费中文在| 新久久久久国产一级毛片| 成在线人永久免费视频| 免费久久久久久久精品成人欧美视频| 亚洲成人免费电影在线观看| 国产精品.久久久| 18禁黄网站禁片午夜丰满| 日日夜夜操网爽| 少妇粗大呻吟视频| 成人免费观看视频高清| 狂野欧美激情性xxxx| 老司机亚洲免费影院| 日本五十路高清| 国产免费av片在线观看野外av| 91精品三级在线观看| 男女边摸边吃奶| 18禁裸乳无遮挡动漫免费视频| 在线观看舔阴道视频| 国产野战对白在线观看| 亚洲国产精品一区二区三区在线| 最近最新中文字幕大全电影3 | 在线观看人妻少妇| 国产精品 国内视频| √禁漫天堂资源中文www| 曰老女人黄片| 热99久久久久精品小说推荐| 久久午夜亚洲精品久久| 99九九在线精品视频| 亚洲少妇的诱惑av| 夜夜爽天天搞| 日本av手机在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 少妇裸体淫交视频免费看高清 | 另类亚洲欧美激情| 天堂动漫精品| 大陆偷拍与自拍| 在线看a的网站| 1024视频免费在线观看| 久久久精品国产亚洲av高清涩受| 夫妻午夜视频| 啦啦啦在线免费观看视频4| 757午夜福利合集在线观看| 老司机靠b影院| 一区福利在线观看| 成人亚洲精品一区在线观看| 国产黄色免费在线视频| av片东京热男人的天堂| 久久久久久久久免费视频了| 激情在线观看视频在线高清 | 国产精品麻豆人妻色哟哟久久| 国产精品.久久久| 男人操女人黄网站| 热re99久久国产66热| 久久久国产欧美日韩av| 亚洲黑人精品在线| 脱女人内裤的视频| 国产又爽黄色视频| 久久精品亚洲av国产电影网| 脱女人内裤的视频| 国产真人三级小视频在线观看| 亚洲黑人精品在线| 精品国内亚洲2022精品成人 | av视频免费观看在线观看| 99riav亚洲国产免费| 精品国产乱子伦一区二区三区| 亚洲全国av大片| 俄罗斯特黄特色一大片| 成人18禁在线播放| 国产成人欧美| 18禁美女被吸乳视频| 国产亚洲午夜精品一区二区久久| 久久久久久免费高清国产稀缺| 蜜桃在线观看..| 母亲3免费完整高清在线观看| 十分钟在线观看高清视频www| 50天的宝宝边吃奶边哭怎么回事| 亚洲av成人不卡在线观看播放网| 一级,二级,三级黄色视频| 最新美女视频免费是黄的| 每晚都被弄得嗷嗷叫到高潮| 日本a在线网址| 久久ye,这里只有精品| 久久婷婷成人综合色麻豆| 露出奶头的视频| 日韩一卡2卡3卡4卡2021年| 日日摸夜夜添夜夜添小说| 国产野战对白在线观看| a级毛片在线看网站| 91大片在线观看| 激情在线观看视频在线高清 | 不卡一级毛片| 极品少妇高潮喷水抽搐| 国产一区二区三区综合在线观看| 国产三级黄色录像| 精品久久久久久久毛片微露脸| 国产精品99久久99久久久不卡| 男女午夜视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 亚洲国产av新网站| 久久久久精品国产欧美久久久| 色老头精品视频在线观看| 日本欧美视频一区| 午夜福利在线免费观看网站| 精品午夜福利视频在线观看一区 | 国产午夜精品久久久久久| 亚洲国产av新网站| 天堂俺去俺来也www色官网| av在线播放免费不卡| 狠狠婷婷综合久久久久久88av| 法律面前人人平等表现在哪些方面| 超碰97精品在线观看| 国产精品自产拍在线观看55亚洲 | 老司机在亚洲福利影院| 少妇精品久久久久久久| 露出奶头的视频| 欧美黄色淫秽网站| 丝袜在线中文字幕| a级毛片在线看网站| 99精品在免费线老司机午夜| h视频一区二区三区| 69精品国产乱码久久久| 国产老妇伦熟女老妇高清| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 精品国产亚洲在线| 亚洲国产精品一区二区三区在线| 在线观看免费高清a一片| 丁香六月天网| 后天国语完整版免费观看| 91九色精品人成在线观看| 欧美变态另类bdsm刘玥| 一区二区av电影网| 成人特级黄色片久久久久久久 | xxxhd国产人妻xxx| 欧美变态另类bdsm刘玥| 狂野欧美激情性xxxx| 天堂俺去俺来也www色官网| 天天添夜夜摸| 精品一区二区三区av网在线观看 | 欧美人与性动交α欧美精品济南到| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇熟女久久| 男女高潮啪啪啪动态图| 免费黄频网站在线观看国产| 69精品国产乱码久久久| 老司机在亚洲福利影院| 免费高清在线观看日韩| 国产区一区二久久| 在线观看www视频免费| 亚洲性夜色夜夜综合| 国产日韩一区二区三区精品不卡| 国产精品麻豆人妻色哟哟久久| 日本av免费视频播放| 中文字幕高清在线视频| 午夜免费成人在线视频| 国产精品久久久av美女十八| 精品国产一区二区三区四区第35| 90打野战视频偷拍视频| 一本—道久久a久久精品蜜桃钙片| 99久久99久久久精品蜜桃| 亚洲av片天天在线观看| 国产日韩欧美在线精品| 亚洲av欧美aⅴ国产| 搡老乐熟女国产| 激情视频va一区二区三区| 成年人免费黄色播放视频| 在线观看免费日韩欧美大片| 啦啦啦中文免费视频观看日本| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲综合一区二区三区_| 亚洲综合色网址| 一二三四在线观看免费中文在| 日韩欧美三级三区| 欧美日韩福利视频一区二区| 精品国产国语对白av| 国产一区有黄有色的免费视频| 交换朋友夫妻互换小说| 色尼玛亚洲综合影院| 亚洲三区欧美一区| 欧美激情极品国产一区二区三区| 热99re8久久精品国产| 一本综合久久免费| 日韩视频在线欧美| 亚洲精华国产精华精| 久久热在线av| 久久久久久久国产电影| 天天添夜夜摸| 2018国产大陆天天弄谢| 亚洲三区欧美一区| 国产一区二区激情短视频| 99久久精品国产亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 欧美黑人精品巨大| 高清毛片免费观看视频网站 | 国产xxxxx性猛交| 久久免费观看电影| 黄片播放在线免费| 在线亚洲精品国产二区图片欧美| 五月开心婷婷网| 成人国产av品久久久| 十八禁网站网址无遮挡| 黄色片一级片一级黄色片| 国产在视频线精品| 成年动漫av网址| tocl精华| 欧美乱码精品一区二区三区| 男女无遮挡免费网站观看| 男女边摸边吃奶| 高潮久久久久久久久久久不卡| 久久人人97超碰香蕉20202| 国产成人精品在线电影| 少妇的丰满在线观看| 欧美人与性动交α欧美软件| 国产xxxxx性猛交| 亚洲精品中文字幕一二三四区 | avwww免费| tocl精华| 成人精品一区二区免费| 国产在线观看jvid| 天天操日日干夜夜撸| 成人18禁在线播放| 妹子高潮喷水视频| 黑丝袜美女国产一区| 色综合欧美亚洲国产小说| 国产精品99久久99久久久不卡| 亚洲精品美女久久av网站| 国产精品美女特级片免费视频播放器 | 久久精品亚洲av国产电影网| 天堂8中文在线网| 女人被躁到高潮嗷嗷叫费观| 老司机影院毛片| 成人三级做爰电影| 黄网站色视频无遮挡免费观看| 女人久久www免费人成看片| 桃花免费在线播放| 又紧又爽又黄一区二区| 亚洲成av片中文字幕在线观看| 国产欧美日韩综合在线一区二区| av国产精品久久久久影院| 日本五十路高清| 交换朋友夫妻互换小说| 国产精品 欧美亚洲| 久久久久久免费高清国产稀缺| 午夜福利在线免费观看网站| 一本—道久久a久久精品蜜桃钙片| 欧美日韩国产mv在线观看视频| tube8黄色片| 大型av网站在线播放| 免费av中文字幕在线| 久热爱精品视频在线9| 亚洲,欧美精品.| 99香蕉大伊视频| 精品久久久精品久久久| 久久久精品免费免费高清| 这个男人来自地球电影免费观看| 午夜免费鲁丝| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| 精品国产乱码久久久久久小说| 无限看片的www在线观看| 在线观看免费视频网站a站| 嫁个100分男人电影在线观看| 日日夜夜操网爽| 国产男靠女视频免费网站| 最近最新中文字幕大全免费视频| 免费久久久久久久精品成人欧美视频| 一区二区三区激情视频| 国产三级黄色录像| 国产单亲对白刺激| 国产精品99久久99久久久不卡| www.999成人在线观看| 精品久久久久久电影网| www日本在线高清视频| 少妇被粗大的猛进出69影院| 成年人黄色毛片网站| 中文字幕av电影在线播放| 少妇精品久久久久久久| 午夜福利在线观看吧| 窝窝影院91人妻| 丁香欧美五月| 黄片大片在线免费观看| 男男h啪啪无遮挡| 新久久久久国产一级毛片| 国产精品 国内视频| www.熟女人妻精品国产| 婷婷成人精品国产| 大陆偷拍与自拍| 一级毛片电影观看| 日韩免费av在线播放| 亚洲av第一区精品v没综合| 欧美精品av麻豆av| 女人被躁到高潮嗷嗷叫费观| 国产亚洲精品久久久久5区| 我的亚洲天堂| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 热re99久久国产66热| 一区福利在线观看| 亚洲天堂av无毛| 国产精品久久久久成人av| 久久天堂一区二区三区四区| 久久香蕉激情| 色婷婷久久久亚洲欧美| 欧美黑人欧美精品刺激| 免费av中文字幕在线| 日韩中文字幕欧美一区二区| 午夜福利免费观看在线| 亚洲国产欧美在线一区| 亚洲成人免费电影在线观看| 国产精品一区二区在线观看99| 免费av中文字幕在线| www.自偷自拍.com| 免费看十八禁软件| 又黄又粗又硬又大视频| 国产午夜精品久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品亚洲一级av第二区| 成在线人永久免费视频| 男女高潮啪啪啪动态图| 好男人电影高清在线观看| 久久亚洲真实| 黄色毛片三级朝国网站| 国产日韩欧美亚洲二区| 最近最新免费中文字幕在线| 99国产精品99久久久久| 日本欧美视频一区| 国产精品成人在线| 午夜福利影视在线免费观看| 国产不卡av网站在线观看| 国产极品粉嫩免费观看在线| 夜夜骑夜夜射夜夜干| 日本黄色日本黄色录像| 狂野欧美激情性xxxx| 纯流量卡能插随身wifi吗| 亚洲精品自拍成人| 午夜福利影视在线免费观看| 99精品欧美一区二区三区四区| 大码成人一级视频| 老熟妇仑乱视频hdxx| 大型av网站在线播放| 脱女人内裤的视频| av欧美777| 国产精品美女特级片免费视频播放器 | 欧美精品人与动牲交sv欧美| 黄色毛片三级朝国网站| 考比视频在线观看| 美女主播在线视频| 一二三四在线观看免费中文在| 久久狼人影院| 后天国语完整版免费观看| av不卡在线播放| 久久久久精品人妻al黑| e午夜精品久久久久久久| 精品少妇一区二区三区视频日本电影| 狠狠婷婷综合久久久久久88av| 国精品久久久久久国模美| 国产又色又爽无遮挡免费看| 中文字幕人妻熟女乱码| 欧美日韩av久久| 黄色视频,在线免费观看| 9色porny在线观看| 国产xxxxx性猛交| 国产精品秋霞免费鲁丝片| 国产精品一区二区精品视频观看| 亚洲中文av在线| 制服人妻中文乱码| 久久精品熟女亚洲av麻豆精品| 男人舔女人的私密视频| 大型av网站在线播放| 亚洲欧美一区二区三区久久| 欧美老熟妇乱子伦牲交| 欧美日韩一级在线毛片| 精品第一国产精品| 国产高清视频在线播放一区| 久久99一区二区三区| 中文字幕高清在线视频| 亚洲男人天堂网一区| 自拍欧美九色日韩亚洲蝌蚪91| 狠狠精品人妻久久久久久综合| 亚洲精品自拍成人| 久久精品国产综合久久久| www.999成人在线观看| 一级片'在线观看视频| 国产精品影院久久| 欧美日韩视频精品一区| 激情在线观看视频在线高清 | 性色av乱码一区二区三区2| 啦啦啦在线免费观看视频4| 在线永久观看黄色视频| 国产精品亚洲av一区麻豆| 黑人巨大精品欧美一区二区蜜桃| 久久性视频一级片| 亚洲一区中文字幕在线| 国产一区二区三区在线臀色熟女 | 香蕉丝袜av| 日日夜夜操网爽| 高清视频免费观看一区二区| 美女扒开内裤让男人捅视频| 国产亚洲精品第一综合不卡| 日本一区二区免费在线视频| 蜜桃国产av成人99| 亚洲五月婷婷丁香| 少妇 在线观看| 亚洲av欧美aⅴ国产| 女性被躁到高潮视频| 精品免费久久久久久久清纯 | 国产成+人综合+亚洲专区| 成年人免费黄色播放视频| 制服人妻中文乱码| 久久国产精品男人的天堂亚洲| 午夜老司机福利片| 丝袜在线中文字幕| 一区二区三区精品91| a级毛片黄视频| 欧美激情极品国产一区二区三区| 制服人妻中文乱码| 久热这里只有精品99| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 亚洲av日韩在线播放| 午夜福利视频精品| 男女床上黄色一级片免费看| 亚洲欧洲精品一区二区精品久久久| 国产在线观看jvid| av片东京热男人的天堂| 日韩制服丝袜自拍偷拍| 日韩有码中文字幕| 欧美日韩黄片免| 91国产中文字幕| 亚洲五月色婷婷综合| 亚洲黑人精品在线| 不卡av一区二区三区| 日韩视频一区二区在线观看| 天堂8中文在线网| 欧美午夜高清在线| 丰满人妻熟妇乱又伦精品不卡| 日韩 欧美 亚洲 中文字幕| 久久久久视频综合| 女人高潮潮喷娇喘18禁视频| 国产成人av教育| 女人被躁到高潮嗷嗷叫费观| 如日韩欧美国产精品一区二区三区| 女同久久另类99精品国产91| 久久精品成人免费网站| 97人妻天天添夜夜摸| 高清黄色对白视频在线免费看| 我要看黄色一级片免费的| 香蕉丝袜av| 欧美人与性动交α欧美软件| 国产国语露脸激情在线看| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说 | 成人特级黄色片久久久久久久 | 怎么达到女性高潮| 亚洲精品美女久久久久99蜜臀| 50天的宝宝边吃奶边哭怎么回事| 亚洲男人天堂网一区| 最新美女视频免费是黄的| 宅男免费午夜| 我的亚洲天堂| 99在线人妻在线中文字幕 | 在线观看舔阴道视频| 叶爱在线成人免费视频播放| 色播在线永久视频| 黑丝袜美女国产一区| 午夜日韩欧美国产| 黑丝袜美女国产一区| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 亚洲九九香蕉| 亚洲专区字幕在线| 免费女性裸体啪啪无遮挡网站| 欧美国产精品va在线观看不卡| 香蕉丝袜av| 757午夜福利合集在线观看| 亚洲七黄色美女视频| 手机成人av网站| 老鸭窝网址在线观看| 在线观看66精品国产| 99riav亚洲国产免费| 婷婷成人精品国产| 精品国产乱码久久久久久小说| 日韩制服丝袜自拍偷拍| 日本vs欧美在线观看视频| 乱人伦中国视频| 欧美黄色淫秽网站| 日本欧美视频一区| 少妇的丰满在线观看| 久久精品熟女亚洲av麻豆精品| 一级毛片电影观看| 国产一卡二卡三卡精品| 久久久久久久久免费视频了| 成年人免费黄色播放视频| 99国产精品99久久久久| 精品国产一区二区三区久久久樱花| 天天躁日日躁夜夜躁夜夜| 国产av精品麻豆| 国产在视频线精品| 肉色欧美久久久久久久蜜桃| 757午夜福利合集在线观看| 亚洲国产欧美一区二区综合| 五月天丁香电影| 国产成人精品久久二区二区91| 激情在线观看视频在线高清 | 成人手机av| 狠狠精品人妻久久久久久综合| 热99re8久久精品国产| 丁香六月天网| 亚洲avbb在线观看| 日韩欧美一区二区三区在线观看 | 色婷婷av一区二区三区视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人啪精品午夜网站| 欧美老熟妇乱子伦牲交| xxxhd国产人妻xxx| 午夜福利免费观看在线| 99国产精品一区二区三区| 久久精品国产亚洲av高清一级| 国产精品久久久久久精品电影小说| 亚洲一区二区三区欧美精品| 欧美亚洲日本最大视频资源| 三级毛片av免费| 日韩一区二区三区影片| 中文字幕人妻熟女乱码| 日本五十路高清| 欧美另类亚洲清纯唯美| 精品乱码久久久久久99久播| 免费一级毛片在线播放高清视频 | 国产成人免费观看mmmm| 99re在线观看精品视频| 搡老乐熟女国产| 国产欧美亚洲国产| 欧美乱码精品一区二区三区| 欧美精品一区二区大全| 亚洲专区国产一区二区| 男人操女人黄网站| 午夜免费成人在线视频| www.精华液| 十八禁网站网址无遮挡| 精品国产国语对白av| 男女下面插进去视频免费观看| 在线观看66精品国产| 久久香蕉激情| tube8黄色片| 国产人伦9x9x在线观看| 男女床上黄色一级片免费看| 亚洲国产欧美网| 人成视频在线观看免费观看| 妹子高潮喷水视频| 免费不卡黄色视频| 国产不卡av网站在线观看| 亚洲av成人一区二区三| 纵有疾风起免费观看全集完整版| 丰满饥渴人妻一区二区三|