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

    山東蒙山森林冠層綠度與樹干徑向生長(zhǎng)的關(guān)系

    2018-01-03 08:20:22方歐婭
    生態(tài)學(xué)報(bào) 2017年22期
    關(guān)鍵詞:綠度植被指數(shù)年輪

    汪 舟, 方歐婭

    1 中國(guó)科學(xué)院植物研究所植被與環(huán)境變化國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100093 2 中國(guó)科學(xué)院大學(xué),北京 100049

    山東蒙山森林冠層綠度與樹干徑向生長(zhǎng)的關(guān)系

    汪 舟1,2, 方歐婭1,*

    1 中國(guó)科學(xué)院植物研究所植被與環(huán)境變化國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100093 2 中國(guó)科學(xué)院大學(xué),北京 100049

    森林冠層綠度和樹木年輪寬度是描述森林生長(zhǎng)過程的重要指標(biāo),它們之間存在怎樣的關(guān)系以及這種關(guān)系的穩(wěn)定性如何目前還沒有清晰的回答。森林冠層綠度通過遙感影像計(jì)算,在空間上連續(xù),而樹木年輪寬度是樹木健康的綜合指標(biāo),樣點(diǎn)上具有代表性。森林冠層綠度和樹木年輪寬度的關(guān)系的研究能增進(jìn)對(duì)森林生長(zhǎng)的多角度理解和森林生長(zhǎng)狀況的尺度轉(zhuǎn)換。在山東蒙山地區(qū)采集了4個(gè)赤松 (Pinusdensiflora) 林樣點(diǎn)的樹木年輪樣本,獲得了樹木年輪寬度數(shù)據(jù),分析了增強(qiáng)型植被指數(shù) (Enhanced Vegetation Index, EVI) 與樹木年輪寬度的關(guān)系。結(jié)果顯示:1)對(duì)于健康森林,4月和6月的冠層綠度與樹木年輪寬度存在因果關(guān)系;森林不健康時(shí),兩者關(guān)系較為復(fù)雜;2)其他月份冠層綠度與樹干徑向生長(zhǎng)不存在因果關(guān)系,而是共同受其他環(huán)境因子,如氣候因子的驅(qū)動(dòng);3)弱冠層綠度降低后5年內(nèi)有顯著的徑向生長(zhǎng)恢復(fù),但是恢復(fù)年份少;強(qiáng)冠層綠度降低之前,樹干徑向生長(zhǎng)已經(jīng)開始降低,之后的5年內(nèi)有著持續(xù)的徑向生長(zhǎng)降低。這些結(jié)果表明森林冠層綠度的降低并不能反映樹干徑向生長(zhǎng)降低的開始,只有健康的森林冠層綠度和年輪寬度有相關(guān)關(guān)系。冠層綠度的降低對(duì)森林健康有強(qiáng)烈的影響,冠層綠度降低導(dǎo)致的徑向生長(zhǎng)的降低很難恢復(fù)。

    年輪寬度;增強(qiáng)型植被指數(shù);格蘭杰因果分析;疊加事件分析;森林衰退;赤松林

    森林的健康生長(zhǎng)關(guān)系到其生態(tài)服務(wù)功能的發(fā)揮和維護(hù),對(duì)森林健康的評(píng)估是保護(hù)生態(tài)學(xué)研究的重要內(nèi)容之一。對(duì)于森林的健康狀況,從外觀上可以通過森林冠層綠度得到直接反映,從樹木內(nèi)部也可通過樹木年輪寬度來反映。過去的研究表明森林冠層綠度與樹木年輪寬度之間存在一定的相關(guān)關(guān)系[1- 3],但另一些研究發(fā)現(xiàn)冠層綠度的短暫變化并不能監(jiān)測(cè)到病蟲害導(dǎo)致的樹干徑向生長(zhǎng)的變化[4]。目前對(duì)于森林冠層綠度與樹干徑向生長(zhǎng)的關(guān)系還缺乏清晰的認(rèn)識(shí)。

    森林冠層綠度是森林光合作用活力和冠層結(jié)構(gòu)的代用指標(biāo)[5],可以表示森林外在的健康狀況,用增強(qiáng)型植被指數(shù)(Enhanced Vegetation Index, EVI)度量。衰退的森林在外觀上主要表現(xiàn)為新葉、葉量、綠葉數(shù)的減少,不正常落葉,樹體結(jié)構(gòu)不完整,小枝和枝干的死亡甚至有些情況下整株樹的死亡[6- 7]。EVI能較好的反映冠層結(jié)構(gòu)包括葉面積指數(shù)和冠層結(jié)構(gòu)的變化[5]。同類的植被指數(shù)在森林健康狀況的研究中也有著廣泛的應(yīng)用[5,7- 11]。植被指數(shù)可以通過遙感數(shù)據(jù)獲得,在時(shí)間和空間上連續(xù),有利于評(píng)估大范圍植物生長(zhǎng)的狀況。

    樹木年輪寬度可以表示樹木徑向生長(zhǎng),其時(shí)間尺度較長(zhǎng)。樹木年輪寬度指數(shù)(Tree Ring Index, TRI)是測(cè)量樹木生理健康的指標(biāo)[12- 13],廣泛應(yīng)用于森林衰退的研究中[14- 20]。外觀上冠層綠度降低而表現(xiàn)出來的森林衰退是容易觀察到的,但是徑向生長(zhǎng)量的變化需要通過每年監(jiān)測(cè)或者鉆取樹芯精確檢測(cè)。樹木年輪寬度與植被指數(shù)在反映森林健康的不同方面,樹木年輪相比于植被指數(shù)能反應(yīng)較長(zhǎng)時(shí)間尺度上的趨勢(shì),能比較準(zhǔn)確的反應(yīng)森林健康程度[12- 13];植被指數(shù)有年周期, 能反應(yīng)短時(shí)間上冠層綠度的變化[7]。對(duì)于不同健康狀況的森林,植被指數(shù)與樹木年輪寬度指數(shù)的關(guān)系并不清楚,冠層綠度對(duì)樹木年輪寬度的影響也不清楚。

    2014年在山東蒙山地區(qū)進(jìn)行森林生態(tài)調(diào)查時(shí)發(fā)現(xiàn)該地的赤松林存在不同程度的斑塊狀衰退現(xiàn)象,外觀上表現(xiàn)為冠層綠度降低、落葉、小枝和枝干的枯死,甚至個(gè)體的死亡。森林冠層綠度的降低是否一定同時(shí)伴隨著樹木徑向生長(zhǎng)的降低?冠層綠度與樹木徑向生長(zhǎng)的關(guān)系是否在時(shí)間上穩(wěn)定?針對(duì)這些問題,選擇了不同健康程度的赤松林為研究對(duì)象,采集了樹木年輪樣本,并分析冠層綠度與樹木年輪的關(guān)系。

    1 研究區(qū)概況

    蒙山位于山東省中部,位于117°35′—118°20′E和35°10′—38°00′N之間,蒙山山脈呈西北-東南走向,是泰沂山脈的支脈,綿延75 km,面積1125 km2,跨費(fèi)縣、沂南、平邑和蒙陰4縣。蒙山主峰龜蒙頂海拔1156 m,為山東省第二高峰。研究區(qū)內(nèi)地質(zhì)構(gòu)造復(fù)雜,山體主要為片麻巖、花崗片麻巖等,山腳有石灰?guī)r覆蓋。土壤類型以棕壤為主,中性至微酸[21]。蒙山屬暖溫帶半濕潤(rùn)季風(fēng)氣候,年平均氣溫約13.2℃,無霜期約196 d,年降雨量約750 mm,最高溫度和最高降雨量均出現(xiàn)于7月份。蒙山森林覆蓋率達(dá)85%以上,因較長(zhǎng)時(shí)間的封山育林,蒙山植被發(fā)育良好,主要有人工落葉闊葉林、人工針葉林、山地灌草叢及小面積的次生落葉闊葉雜木林[22]。

    本研究的樹種為赤松,赤松是我國(guó)暖溫帶沿海地區(qū)溫性針葉林的主要建群種之一,赤松林是山東地帶性的針葉林之一,也是山東省面積最大、分布最廣、資源最豐富的天然森林群落。赤松的年輪清晰,尤其適合樹木年輪生態(tài)學(xué)研究[23]。采樣點(diǎn)所在的塔山林場(chǎng)和大青山林場(chǎng)主要樹種為赤松,造林時(shí)間為20世紀(jì)50年代末60年代初。林分密度,覆蓋度較為一致,種源和林齡相同,后期人為干擾較少。

    圖1 采樣點(diǎn)概況圖Fig.1 Forests of the tree-ring sampling sites(a)采樣點(diǎn)位置和地勢(shì)示意圖;(b)大青山林場(chǎng)五彩山(FXWC)森林衰退景觀,山頂灰褐色為衰退的林分,山下近處可見衰退樹;(c)大青山林場(chǎng)小旺(FXXW)道路旁森林衰退景觀

    2 數(shù)據(jù)與方法

    2.1 樹木年輪數(shù)據(jù)

    選擇4個(gè)樹木年輪采樣點(diǎn),分別在蒙山塔山林場(chǎng)內(nèi)水庫(kù)上的低山山頂(FXTS),塔山林場(chǎng)葫蘆崖山腰(FXHL),大青山林場(chǎng)內(nèi)的五彩山(FXWC)和小旺(FXXW),見圖1。在這4個(gè)采樣點(diǎn)中,FXTS樣點(diǎn)的樹木外觀比較健康,樹體結(jié)構(gòu)完整,針葉深綠,FXHL樣點(diǎn)的樹木有落葉、褐色枯葉的衰退現(xiàn)象,而費(fèi)縣大青山林場(chǎng)的兩個(gè)采樣點(diǎn)FXWC和FXXW的樹木在外觀上表現(xiàn)出高比例的落葉、褐色枯葉、死枝等衰退現(xiàn)象。FXXW和FXWC的樹木在1976年爆發(fā)病蟲害。

    采樣時(shí),用生長(zhǎng)錐在樹木胸高處與山坡走向平行方向鉆取樹芯。將取到的樹芯放置在塑料管內(nèi),并在管上標(biāo)注編碼。在實(shí)驗(yàn)室內(nèi),將樹芯置于陰涼處風(fēng)干,用白乳膠將樣條粘到木槽上,使樹芯的木纖維方向與槽面垂直,用細(xì)繩將樹芯纏繞固定在木槽上。待白乳膠干后,拆下繩子,用280—600目的砂紙進(jìn)行打磨,直到樹輪在顯微鏡下清晰可見。各點(diǎn)采樣的樣本量見表1。

    在實(shí)驗(yàn)室內(nèi),使用年輪寬度測(cè)量?jī)xLinTab5[24]測(cè)量每個(gè)樹芯樣本的年輪寬度,測(cè)量精度為0.001 mm。根據(jù)樣芯最后一年的解剖學(xué)特征確定靠近樹皮最外一個(gè)年輪的形成年份。用樹木年輪學(xué)通用的標(biāo)記符號(hào)標(biāo)記偽輪,窄輪,缺年等。之后用骨架示意圖進(jìn)行交叉定年[25]。用COFECHA軟件對(duì)定年和測(cè)量結(jié)果進(jìn)行了檢驗(yàn),保證定年與測(cè)量的準(zhǔn)確性。從樹木年輪上看,大青山林場(chǎng)五彩山和小旺的部分樹木在2014年采樣時(shí)有極窄的年輪或已經(jīng)不再有徑向生長(zhǎng)。大青山林場(chǎng)五彩山有3個(gè)樹芯在1985年和1992年存在缺輪,大青山林場(chǎng)小旺有4個(gè)樹芯在1991年,1992年和1993年存在缺輪。

    平均采樣點(diǎn)各樣芯年輪寬度指數(shù)可以得到理論平均值為1的無量綱樹輪寬度標(biāo)準(zhǔn)年表[26]。因?yàn)橐芯抗趯泳G度與樹木徑向生長(zhǎng)的關(guān)系,本文主要采用了水平去趨勢(shì)的方法,即樹芯每年的年輪寬度除以整個(gè)時(shí)間序列的平均年輪寬度,這樣可以保留長(zhǎng)時(shí)間樹木個(gè)體的生長(zhǎng)趨勢(shì)[17,27]。對(duì)有明顯年齡趨勢(shì)的樹芯仍用負(fù)指數(shù)函數(shù)去趨勢(shì)。

    表1 采樣點(diǎn)信息表

    FXTS: 費(fèi)縣塔山林場(chǎng)水庫(kù);FXHL: 費(fèi)縣塔山林場(chǎng)葫蘆崖; FXXW: 費(fèi)縣大青山林場(chǎng)小旺; FXWC: 費(fèi)縣大青山林場(chǎng)五彩山

    2.2 植被指數(shù)數(shù)據(jù)及冠層綠度降低

    植被指數(shù)產(chǎn)品選用 MOD13Q1-EVI[28]數(shù)據(jù),引用自美國(guó)航空航天局(NASA)數(shù)據(jù)中心(https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod13q1_v006),其空間分辨率為250 m×250 m,時(shí)間分辨率為16 d。MOD13Q1在每個(gè)像元上提供了6種植被指數(shù),本文采用的是其中的增強(qiáng)植被指數(shù)(EVI)。該指數(shù)已做太陽(yáng)輻射矯正,大氣矯正和氣溶膠的影響的低,最大限度地減少了冠層背景變化,保持高植被密度的敏感性,更適用于高植被覆蓋區(qū)。

    EVI應(yīng)用于植物群落時(shí)代表該地區(qū)森林綠度,即反映該地區(qū)植被數(shù)量和植被健康程度或生長(zhǎng)狀態(tài)。EVI的單年變化顯示了植被冠層由春季開始綠度上升至夏季結(jié)束綠度下降的過程,反映了植被的物候特征。異常的氣候條件或突發(fā)的蟲害會(huì)顯著影響植被的頻譜特征[7]。本研究中,年際尺度的EVI變化用于監(jiān)測(cè)和發(fā)現(xiàn)森林冠層綠度的變化。

    盡管16 d最大值合成的MODIS地表反射率數(shù)據(jù)在一定程度上去除了云和大氣的影響,但仍然可能受太陽(yáng)高度角、觀測(cè)角度、氣溶膠等一些隨機(jī)因素的干擾,使得時(shí)間序列數(shù)據(jù)的變化呈不規(guī)則狀態(tài),通常表現(xiàn)為EVI曲線短期起伏很大[29],這是不符合植被正常生長(zhǎng)狀況的,因而需要對(duì)合成的MODIS EVI數(shù)據(jù)產(chǎn)品進(jìn)一步濾波,去除噪聲干擾。本文選用Savitzky-Golay濾波函數(shù)平滑時(shí)間序列EVI[30](圖2)。

    部分研究者[31]用生長(zhǎng)季植被指數(shù)的變化量表征冠層綠度降低,由于EVI短期起伏大,部分年份存在“高尖”現(xiàn)象,會(huì)產(chǎn)生誤差。部分研究者[32]利用生長(zhǎng)季EVI的累積EVIgs(growth season EVI)表示單年植被生長(zhǎng)季綠度的積累,這種綠度的積累一定程度上表示了光合產(chǎn)物的積累[7,33]。

    圖2 EVI原始時(shí)間序列及Savitzky-Golay 濾波后的時(shí)間序列(以塔山林場(chǎng)葫蘆崖為例)Fig.2 Raw time series of 16-day EVI and time series through Savitzky-Golay smoothing (FXHL as an example)

    綜上,本文預(yù)選兩種冠層綠度降低指標(biāo),一種是生長(zhǎng)季4—10月份的逐月EVI值。另一種是生長(zhǎng)季EVI的累積EVIgs,表示年植被生長(zhǎng)季冠層綠度的積累。EVI值大于0.25的時(shí)間段作為生長(zhǎng)季(表2),繼而對(duì)生長(zhǎng)季的濾波后的EVI進(jìn)行積分,積分后的值作為植被健康指標(biāo),以年為時(shí)間分辨率。這些備選的冠層綠度指標(biāo)用于探究冠層綠度與年輪寬度的穩(wěn)定關(guān)系。

    表2 蒙山4個(gè)采樣點(diǎn)森林從2000—2014年從16d EVI提取的生長(zhǎng)季長(zhǎng)度

    2.3 增強(qiáng)型植被指數(shù)與年輪寬度指數(shù)的因果關(guān)系分析

    格蘭杰因果關(guān)系檢驗(yàn)以隨機(jī)過程的線性回歸模型來度量時(shí)間序列之間的因果性[34]。如果一個(gè)時(shí)間序列的過去值對(duì)另外一個(gè)時(shí)間序列的將來值的預(yù)測(cè)有促進(jìn)作用,則這兩個(gè)時(shí)間序列就存在著一定的因果關(guān)系。年輪寬度指數(shù)序列(TRI)與增強(qiáng)型植被指數(shù)序列(EVI)為兩個(gè)在生理上相關(guān)的時(shí)間序列,本文采用格蘭杰因果關(guān)系檢驗(yàn)探索年輪寬度與冠層綠度的因果關(guān)系。

    在使用格蘭杰檢驗(yàn)方法進(jìn)行檢驗(yàn)時(shí),欲判斷可能的原因要素是否引起了一個(gè)結(jié)果要素,則是在考察被觀察的結(jié)果要素自身過去值對(duì)當(dāng)前值的解釋度的基礎(chǔ)上,加入原因要素的滯后值,看新的回歸方程是否能夠顯著增加對(duì)結(jié)果要素的解釋程度。因果關(guān)系的成立需要滿足以下兩個(gè)條件:一是原因要素的加入能夠顯著增加回歸方程對(duì)結(jié)果要素的解釋度;二是結(jié)果要素不能夠有助于對(duì)原因要素的預(yù)測(cè)。

    在格蘭杰因果檢驗(yàn)前,需要檢驗(yàn)時(shí)間序列的平穩(wěn)性,通過自相關(guān)和偏相關(guān)檢驗(yàn)、單位根檢驗(yàn),單月的EVI和濾波后EVI生長(zhǎng)期累積EVIgs為平穩(wěn)時(shí)間序列。一般來說樹木年輪也是平穩(wěn)時(shí)間序列。

    構(gòu)建格蘭杰因果關(guān)系檢驗(yàn)的隨機(jī)過程的線性回歸模型:

    EVIijt=α+βEVIijt-1+μit

    (1)

    EVIijt=α+β1EVIijt-1+β2TRIit-1+μit

    (2)

    TRIit=α+βTRIit-1+μit

    (3)

    TRIit=α+β1TRIit-1+β2EVIijt-1+μit

    (4)

    (5)

    其中公式(1)和(3)是限制條件的回歸,即公式(5)中的R,公式(2)和(4)是無限制條件的回歸,即公式(5)中的UR。其中EVIijt表示i采樣點(diǎn),t年j月份,TRIit是i采樣點(diǎn)t年的年輪寬度指數(shù)。α,β,β1,β2是回歸系數(shù),μ是回歸殘差項(xiàng)。假設(shè)檢驗(yàn)為F檢驗(yàn)(公式(5)),其中q為自回歸滯后,t為時(shí)間序列長(zhǎng)度。

    2.4 冠層綠度降低與徑向生長(zhǎng)的疊加事件分析

    疊加事件分析用來檢驗(yàn)平均樹木生長(zhǎng)對(duì)某事件響應(yīng)的顯著性[35- 37],該方法將年輪寬度根據(jù)“關(guān)鍵事件”分成幾類,徑向生長(zhǎng)對(duì)關(guān)鍵事件的響應(yīng)會(huì)從均值中顯現(xiàn)出來,而其他的噪音會(huì)消失,其置信度是利用分塊自舉法重新取樣10000次計(jì)算得到。疊加事件分析的關(guān)鍵事件是每個(gè)采樣點(diǎn)冠層綠度降低發(fā)生的年份。

    3 結(jié)果

    3.1 增強(qiáng)型植被指數(shù)與年輪寬度指數(shù)的因果關(guān)系

    通過格蘭杰因果關(guān)系分析,本研究對(duì)蒙山4個(gè)采樣點(diǎn)的增強(qiáng)型植被指數(shù)與年輪寬度指數(shù)的關(guān)系進(jìn)行了探究。

    圖3 蒙山4個(gè)采樣點(diǎn)的年輪寬度指數(shù)(TRI)和增強(qiáng)型植被指數(shù)(EVI)的格蘭杰因果關(guān)系檢驗(yàn)(P =0.1)Fig.3 The Granger Causality Test between Ring Width Index and Enhanced Vegetation Index of four sample sites in Mengshan Mountains (P =0.1)

    塔山林場(chǎng)水庫(kù)上的低山山頂(FXTS)采樣點(diǎn)的樹木外觀比較健康,格蘭杰因果檢驗(yàn)結(jié)果表明4月份,6月份和7月份,EVI引起了TRI的變化,而TRI沒有引起EVI的變化。EVI和TRI是第1種關(guān)系:冠層綠度的變化導(dǎo)致了年輪寬度的變化。

    塔山林場(chǎng)葫蘆崖(FXHL)采樣點(diǎn)的樹木有落葉、褐色枯葉的衰退現(xiàn)象,格蘭杰因果檢驗(yàn)結(jié)果表明4月份和6月份EVI引起了TRI的變化,而TRI沒有引起EVI的變化。屬于第1種關(guān)系:冠層綠度的變化導(dǎo)致了年輪寬度的變化。5月份,7—10月份以及EVIgs格蘭杰檢驗(yàn)得EVI和TRI是第4種關(guān)系:某種變量導(dǎo)致了冠層綠度和年輪寬度兩者的變化。大青山林場(chǎng)小旺(FXXW)采樣點(diǎn)的樹木外觀表現(xiàn)出嚴(yán)重的衰退,5—10月份以及EVIgs格蘭杰因果檢驗(yàn)結(jié)果得EVI和TRI是第4種關(guān)系:某種變量導(dǎo)致了冠層綠度和年輪寬度兩者的變化。大青山林場(chǎng)五彩山(FXWC)采樣點(diǎn)的樹木外觀表現(xiàn)出嚴(yán)重的衰退,5—10月份以及EVIgs格蘭杰因果檢驗(yàn)結(jié)果得EVI和TRI是第4種關(guān)系:某種變量導(dǎo)致了冠層綠度和年輪寬度兩者的變化。綜上,對(duì)于健康的森林,4月份和6月份EVI是TRI變化的格蘭杰原因(P<0.1),對(duì)于衰退的森林,EVI和TRI互相不是變化的原因,可能由于某種變量導(dǎo)致了EVI和TRI的變化。

    圖4 FXTS的年輪寬度和逐月溫度、降水的相關(guān)性分析Fig.4 The correlation of ring width index of FXTS and monthly temperature and precipitation in meteorological stations

    圖5 FXHL的年輪寬度和逐月溫度、降水的相關(guān)性分析Fig.5 The correlation of ring width index of FXHL and monthly temperature and precipitation in meteorological stations

    圖6 FXXW的年輪寬度和逐月溫度、降水的相關(guān)性分析 Fig.6 The correlation of ring width index of FXXW and monthly temperature and precipitation in meteorological stations

    氣候因子與年輪寬度指數(shù)生長(zhǎng)季4—10月份的滑動(dòng)相關(guān)系數(shù)用自舉法計(jì)算,窗口長(zhǎng)度為20,滑動(dòng)步長(zhǎng)為2年。如圖4—圖7,從逐月的臨沂氣象站(1951—1997年)和費(fèi)縣氣象站(1998—2014年)實(shí)測(cè)溫度和降水與年輪寬度指數(shù)的相關(guān)關(guān)系上來看,4月份溫度,6月份的溫度和6月份的降水與年輪寬度指數(shù)沒有在時(shí)間上穩(wěn)定的相關(guān)性,有些年份正相關(guān),有些年份負(fù)相關(guān)。4月份的降水和年輪寬度指數(shù)時(shí)間上穩(wěn)定的相關(guān)性。

    3.2 冠層綠度降低年份

    根據(jù)格蘭杰因果檢驗(yàn)的結(jié)果,4月份和6月份的森林冠層綠度是徑向生長(zhǎng)的原因,故用4月份和6月份的EVI值來確定冠層綠度降低事件發(fā)生的年份。4月份和6月份的EVI極低值大多發(fā)生在相同的年份,或者是相差1年。若6月份的EVI極低值發(fā)生在當(dāng)年,后一年4月發(fā)生EVI極低值時(shí),EVI極低年份確定為當(dāng)年,因?yàn)楹芸赡茉诠趯泳G度快速增加的4月受前一年綠度達(dá)到最高的月份6—8月份的生物量積累影響,如果前一年的物質(zhì)積累不夠,當(dāng)年生長(zhǎng)季開始綠度增加也會(huì)不夠。若極低值在4月份和6月份發(fā)生的年份不同,兩者的極低值年都作為極端年。根據(jù)極低值的排序,確定弱的冠層綠度降低事件和強(qiáng)的冠層綠度降低事件。強(qiáng)的冠層綠度降低事件設(shè)定為4月份EVI的降低和6月份EVI的降低接連發(fā)生或4, 6兩月份中有極低的EVI。在閾值的設(shè)定上,參考了Tovar對(duì)衰退森林植被綠度的研究[38]。弱的冠層綠度降低事件設(shè)定為4月份EVI小于0.23且6月份EVI小于0.37或4月份和6月份兩月都有較低的EVI。根據(jù)冠層綠度降低年份的確定條件,塔山林場(chǎng)水庫(kù)(FXTS),塔山林場(chǎng)葫蘆崖(FXHL),大青山林場(chǎng)小旺(FXXW),大青山林場(chǎng)五彩山(FXWC)的強(qiáng)的冠層綠度降低事件分別發(fā)生在2000年,2009年,2012年;2009年,2013年; 2007年,2013年;2003年,2008年,2012年。弱的冠層綠度降低事件分別發(fā)生在2003年,2006年;2001年,2014年;2000年,2003年,2010年;2001年,2005年。

    3.3 冠層綠度降低與徑向生長(zhǎng)的疊加事件分析

    確定了冠層綠度降低事件發(fā)生的時(shí)間后,通過疊加事件分析來分析冠層綠度降低事件對(duì)徑向生長(zhǎng)的影響。在冠層綠度降低事件的前5年內(nèi),徑向生長(zhǎng)可能已經(jīng)有降低的預(yù)兆,在冠層綠度降低事件的后5年內(nèi),徑向生長(zhǎng)可能在一些年份恢復(fù),也可能很難恢復(fù)。冠層綠度降低的強(qiáng)弱程度以及森林健康狀況對(duì)樹木徑向生長(zhǎng)的恢復(fù)有顯著的影響。

    在弱的冠層綠度降低事件之后,4個(gè)采樣點(diǎn)的森林徑向生長(zhǎng)都有某些年份的顯著恢復(fù),兩個(gè)健康的采樣點(diǎn)在顯著的恢復(fù)年恢復(fù)力較強(qiáng),FXTS恢復(fù)年的年輪寬度指數(shù)相比前一年增長(zhǎng)0.88,FXHL年輪寬度指數(shù)增長(zhǎng)達(dá)0.91。兩個(gè)衰退的采樣點(diǎn)在顯著的恢復(fù)年恢復(fù)力弱,年輪寬度指數(shù)增長(zhǎng)未達(dá)0.50(圖9)。

    圖7 FXWC的年輪寬度和逐月溫度、降水的相關(guān)性分析Fig.7 The correlation of ring width index of FXWC and monthly temperature and precipitation in meteorological stations

    圖8 4個(gè)采樣點(diǎn)的年輪寬度指數(shù)(TRI),4月份,6月份增強(qiáng)型植被指數(shù)(EVI)的時(shí)間序列Fig.8 Average ring width index and enhanced vegetation index time series of April and June of four sample sites

    圖9 4個(gè)采樣點(diǎn)年輪寬度指數(shù)與弱冠層綠度降低的疊加事件分析Fig.9 Superposed Epoch Analysis between ring width index and moderate canopy green loss events of four sample sites in Mengshan Mountains黑色表示顯著的年份,灰色表示不顯著的年份,P =0.05

    圖10 4個(gè)采樣點(diǎn)年輪寬度指數(shù)與強(qiáng)冠層綠度降低的疊加事件分析 Fig.10 Superposed Epoch Analysis between ring width index and severe canopy green loss events of four sample sites in Mengshan Mountains黑色表示顯著的年份,灰色表示不顯著的年份,P =0.05

    對(duì)于健康的森林,年輪寬度指數(shù)與冠層綠度降低事件有相關(guān)關(guān)系,冠層綠度降低當(dāng)年年輪寬度指數(shù)顯著降低。因?yàn)楣趯泳G度降低前后都有年輪寬度指數(shù)顯著降低,故冠層綠度降低與樹干徑向生長(zhǎng)降低不存在因果關(guān)系。在強(qiáng)的冠層綠度降低事件之后,4個(gè)采樣點(diǎn)年輪寬度都有顯著的降低,并且是連續(xù)的持續(xù)時(shí)間長(zhǎng)的顯著降低。4個(gè)采樣點(diǎn)中強(qiáng)的冠層綠度降低事件之前的1年徑向生長(zhǎng)都降低了。在強(qiáng)的冠層綠度降低事件之后徑向生長(zhǎng)在近5年幾乎不存在恢復(fù)。2000—2014年中 2—3個(gè)強(qiáng)的冠層綠度降低事件影響的時(shí)間長(zhǎng),強(qiáng)度大(圖10)。

    4 討論4.1 森林健康與EVI

    4個(gè)采樣點(diǎn)之間EVI閾值法定義的生長(zhǎng)季長(zhǎng)度有顯著的差異(P<0.01),塔山林場(chǎng)兩個(gè)較健康采樣點(diǎn)的森林比大青山林場(chǎng)兩個(gè)衰退采樣點(diǎn)的森林的生長(zhǎng)季平均長(zhǎng)了38 d,38 d是通過MODIS EVI數(shù)據(jù)計(jì)算出來的,并不是準(zhǔn)確的天數(shù),因?yàn)镸ODIS EVI數(shù)據(jù)的時(shí)間分辨率是16天。4個(gè)采樣點(diǎn)之間4月份和6月份的EVI也都有顯著差異(4月份P=0.0180,6月份P=0.0033)。4個(gè)采樣點(diǎn)每年的植被指數(shù)從最低到最高的變化值差異不顯著(P=0.6304),若作為森林健康程度的指標(biāo)可能并不能真正表示森林的健康程度,冠層綠度降低通常是短暫的,可能出現(xiàn)在個(gè)別月份但是對(duì)整年植被指數(shù)從低到高的變化沒有太大影響,這與Debeurs等的結(jié)果一致,Debeurs認(rèn)為粗的時(shí)間分辨率很難動(dòng)態(tài)監(jiān)測(cè)病蟲害導(dǎo)致的短暫的冠層綠度降低事件,短暫冠層綠度降低往往不影響植被指數(shù)的最低最高值[4]。

    4.2 EVI與TRI的關(guān)系

    4月份和6月份的EVI是TRI變化的格蘭杰原因在健康森林中明顯有所體現(xiàn)。這與Kaufmann等[1- 3]對(duì)北美和亞歐北方48個(gè)點(diǎn)年輪寬度和NDVI格蘭杰因果關(guān)系計(jì)算的結(jié)果較為一致。Kaufmann等認(rèn)為針葉林為主的北方森林6月份和7月份的NDVI與TRI之間相關(guān)關(guān)系最強(qiáng),3—5月份及8—10月份的相關(guān)關(guān)系弱,然而生長(zhǎng)季TRI與NDVI序列之間不存在顯著的因果關(guān)系[2]。氣候驅(qū)動(dòng)TRI和NDVI的關(guān)系,這種氣候的控制區(qū)域與緯度更相關(guān)而不是樹種[3]。格蘭杰因果關(guān)系檢驗(yàn)結(jié)果的差異主要來自3個(gè)方面:樹木年輪去趨勢(shì)方法,回歸方程以及植被指數(shù)的數(shù)據(jù)源[39]。在做相關(guān)性分析時(shí),去趨勢(shì)方法對(duì)結(jié)果的影響需要考慮[40]。Kaufmann 等的年輪寬度序列為標(biāo)準(zhǔn)化年表,一般標(biāo)準(zhǔn)化年表去掉了單株樹木年輪寬度的年齡趨勢(shì)和內(nèi)源性的干擾影響[41-42],去趨勢(shì)的過程中很可能去掉了與冠層綠度相關(guān)的一些趨勢(shì)。本文主要用水平去趨勢(shì),有明顯年齡趨勢(shì)的用負(fù)指數(shù)函數(shù)去趨勢(shì),這樣可以保留長(zhǎng)時(shí)間樹木個(gè)體的生長(zhǎng)趨勢(shì)[17,27]。在考慮植被指數(shù)和年輪寬度的回歸關(guān)系時(shí),太陽(yáng)高度角和氣溶膠光學(xué)厚度也作為回歸因子,以降低其對(duì)NDVI的影響[1- 2]。研究植被指數(shù)和年輪寬度序列的關(guān)系時(shí)大多使用8 km分辨率的GIMMS NDVI數(shù)據(jù)集[1- 2,43]和PAL NDVI數(shù)據(jù)集[39],8 km × 8 km像元中的地物可能不僅是植被,也可能不僅是單一物種的植被,同時(shí)樹木年輪采樣能代表的空間觀測(cè)范圍并沒有這么大,粗的分辨率大大降低了植被指數(shù)和樹木年輪寬度關(guān)系的準(zhǔn)確性,影響植被指數(shù)和年輪寬度關(guān)系的結(jié)果。本文使用的是MOD13Q1時(shí)間分辨率16天,空間分辨率250 m×250 m 的EVI三級(jí)遙感產(chǎn)品,太陽(yáng)輻射和氣溶膠影響的較低。

    對(duì)于健康的森林,4月份和6月份EVI是TRI變化的格蘭杰原因(P<0.1)可能是因?yàn)樯L(zhǎng)季初期森林冠層的快速變綠(green-up)與徑向生長(zhǎng)關(guān)系密切,森林冠層快速變綠時(shí)期的冠層綠度反映了生長(zhǎng)季初期植被變化的活力[5,33]。Rossi等的研究顯示新的管胞的形成在夏季早期之后就會(huì)降低[44],因?yàn)槟贻唽挾绕毡榕c總的管胞數(shù)量強(qiáng)相關(guān)[45],TRI也和夏季開始階段的樹木狀態(tài)有強(qiáng)相關(guān)。對(duì)于衰退的森林,可能是某種因子驅(qū)動(dòng)EVI和TRI的變化。在衰退中徑向生長(zhǎng)體現(xiàn)的是樹木生長(zhǎng)的趨勢(shì),EVI體現(xiàn)的是年內(nèi)冠層的周期變化,冠層產(chǎn)生的營(yíng)養(yǎng)物質(zhì)儲(chǔ)存在莖中,莖為生長(zhǎng)季的冠層生長(zhǎng)提供營(yíng)養(yǎng)物質(zhì),冠層綠度和徑向生長(zhǎng)的因果關(guān)系和營(yíng)養(yǎng)物質(zhì)分配有關(guān)[46-47]。

    Beck等采用了重復(fù)很好的4個(gè)阿拉斯加和加拿大西部沿著東西梯度的靠近北極地區(qū)的林線采樣的樹木年輪數(shù)據(jù)做研究,結(jié)果顯示TRI和NDVI以及溫度都沒有很好的相關(guān)性,最大晚材密度MXD和NDVI在生長(zhǎng)季早期有顯著正相關(guān)(P<0.05),可能是因?yàn)橛行┎蓸狱c(diǎn)脫離了溫度限制的區(qū)域[43]。Bunn等研究了19對(duì)西伯利亞針葉林TRI和MXD數(shù)據(jù)集和NDVI的關(guān)系,結(jié)果顯示TRI和春季及夏季的NDVI有中等程度和中等時(shí)間連續(xù)性的相關(guān),MXD與NDVI有更明確的相關(guān)[40]。

    4.3 冠層綠度降低事件對(duì)樹木徑向生長(zhǎng)的影響

    雖然4個(gè)采樣點(diǎn)森林在外觀上的健康水平不同,但是通過4月份和6月份EVI值篩選出來2000—2014年這15年間強(qiáng)冠層綠度降低的年份個(gè)數(shù)都為2—3個(gè),強(qiáng)的冠層綠度降低事件影響的時(shí)間長(zhǎng),強(qiáng)度大。年輪寬度對(duì)強(qiáng)的冠層綠度降低事件的響應(yīng)表現(xiàn)為事件當(dāng)年徑向生長(zhǎng)顯著降低,前一年徑向生長(zhǎng)降低,在前2—5年中徑向生長(zhǎng)有顯著或者不顯著的生長(zhǎng)增加。在強(qiáng)的冠層綠度降低事件之后,4個(gè)采樣點(diǎn)年輪寬度降低顯著,持續(xù)時(shí)間長(zhǎng),5年內(nèi)徑向生長(zhǎng)幾乎沒有恢復(fù)的年份。強(qiáng)的冠層綠度降低事件對(duì)于森林的健康是非常不利的,持續(xù)的徑向生長(zhǎng)降低的結(jié)果就是死亡[16]。費(fèi)縣大青山小旺(FXXW)和五彩山(FXWC)兩個(gè)采樣點(diǎn)在2013年和2014年就發(fā)生了樹木死亡。

    在弱的冠層綠度降低之后,兩個(gè)健康的采樣點(diǎn)在顯著的恢復(fù)年恢復(fù)力較強(qiáng),兩個(gè)衰退的采樣點(diǎn)在顯著的恢復(fù)年恢復(fù)力弱。弱冠層綠度降低之后的5年4個(gè)采樣點(diǎn)的森林徑向生長(zhǎng)都有恢復(fù)且某一年份恢復(fù)顯著。弱的冠層綠度降低后雖然有恢復(fù),但是恢復(fù)力和森林的健康狀況有關(guān)。短暫的干擾如蟲害和氣候變化對(duì)冠層的傷害可以通過保持氮限制顯著影響生態(tài)系統(tǒng)的功能,氮限制最終會(huì)導(dǎo)致森林受到干旱,蟲害,氣候變化等后續(xù)短暫干擾的脅迫從而易于死亡[48]。有研究指出,蟲害干擾,如松葉蜂危害各個(gè)齡級(jí)的森林,造成歐洲赤松和黑松嚴(yán)重的冠層落葉,胸徑斷面積降低40%—70%[49]。Poljan?ek等追蹤歐洲黑松和赤松的針葉動(dòng)態(tài)以及測(cè)量樹木年輪寬度,結(jié)果顯示針葉簇的減少導(dǎo)致了當(dāng)年徑向生長(zhǎng)和高生長(zhǎng)的衰退。氣候脅迫事件導(dǎo)致的針葉存留比例在葉齡上有區(qū)別[50]??釤岷头磸?fù)的干旱事件同時(shí)導(dǎo)致了生物量的降低,缺少糖類及營(yíng)養(yǎng),森林的恢復(fù)至少需要4年[51]。還會(huì)導(dǎo)致弱的枝條長(zhǎng)度增加,不正常的短針葉,每簇針葉數(shù)量少,少的分枝和結(jié)果率[51]。目前的研究仍無法準(zhǔn)確獲知冠層綠度降低以及森林衰退的原因,但個(gè)體樹木的衰退通常是長(zhǎng)時(shí)間尺度上的脅迫和短時(shí)間尺度上的脅迫的聯(lián)合效應(yīng),森林衰退是一系列環(huán)境脅迫,易感因子,誘導(dǎo)因子以某種發(fā)生順序引起的森林外貌的病態(tài)[6],這種衰退很難恢復(fù)。

    5 結(jié)論

    蒙山地區(qū)健康森林中,冠層綠度與樹木年輪寬度在4月份和6月份存在因果關(guān)系;森林不健康時(shí),兩者無因果關(guān)系,關(guān)系較為復(fù)雜。

    冠層綠度降低與樹干徑向生長(zhǎng)降低不存在因果關(guān)系,但是冠層綠度降低對(duì)樹干徑向生長(zhǎng)有影響:在強(qiáng)的冠層綠度降低事件造成赤松徑向生長(zhǎng)的持續(xù)降低且難恢復(fù)。強(qiáng)冠層綠度降低事件可能造成的徑向生長(zhǎng)的持續(xù)下降似乎是不可逆轉(zhuǎn)的,甚至可能導(dǎo)致樹木死亡。弱冠層綠度降低不一定伴隨著樹干徑向生長(zhǎng)的降低,之后某些年份恢復(fù)顯著,恢復(fù)力與森林的健康狀況有關(guān)。

    致謝:張齊兵老師對(duì)研究工作給予幫助,邱紅巖老師幫助測(cè)定樹木年輪,山東省林業(yè)廳、費(fèi)縣塔山林場(chǎng)和大青山林場(chǎng)給予工作支持,聶昌遠(yuǎn)博士和梁寒雪博士對(duì)野外工作給予幫助,特此致謝。

    [1] Kaufmann R K, Zhou L, Myneni R B, Tucker C J, Slayback D, Shabanov N V, Pinzon J. The effect of vegetation on surface temperature: a statistical analysis of NDVI and climate data. Geophysical Research Letters, 2003, 30(22): 2147.

    [2] Kaufmann R K, D′Arrigo R D, Laskowski C, Myneni R B, Zhou L, Davi N K. The effect of growing season and summer greenness on northern forests. Geophysical Research Letters, 2004, 31(9): L09205.

    [3] Kaufmann R K, D′Arrigo R D, Paletta L F, Tian H Q, Jolly W M, Myneni R B. Identifying climatic controls on ring width: the timing of correlations between tree rings and NDVI. Earth Interactions, 2008, 12(14): 1- 14.

    [4] De Beurs K M, Townsend P A. Estimating the effect of gypsy moth defoliation using MODIS. Remote Sensing of Environment, 2008, 112(10): 3983- 3990.

    [5] Pettorelli N, Vik J O, Mysterud A, Gaillard J M, Tucker C J, Stenseth N C. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology & Evolution, 2005, 20(9): 503- 510.

    [6] Ciesla W M, Donaubauer E. Decline and dieback of trees and forests: a global overview. Fao Forestry Paper, 1994, 120: 4- 17.

    [7] Lambert J, Drenou C, Denux J P, Balent G, Cheret V. Monitoring forest decline through remote sensing time series analysis. Giscience & Remote Sensing, 2013, 50(4): 437- 457.

    [8] Linderholm H W. Growing season changes in the last century. Agricultural and Forest Meteorology, 2006, 137(1/2): 1- 14.

    [9] Carnicer J, Coll M, Ninyerola M, Pons X, Sánchez G, Peuelas J. Widespread crown condition decline, food web disruption, and amplified tree mortality with increased climate change-type drought. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108(4): 1474- 1478.

    [10] 高廣磊, 信忠保, 丁國(guó)棟, 李叢叢, 張佳音, 梁文俊, 安云, 賀宇, 肖萌, 李文葉. 基于遙感技術(shù)的森林健康研究綜述. 生態(tài)學(xué)報(bào), 2013, 33(6): 1675- 1689.

    [11] Zoran M A, Dida A I, Zoran L F V. Remote sensing of climate changes effects on forest biophysical variables//Neale C M U, Maltese A, eds. Remote Sensing for Agriculture, Ecosystems, and Hydrology XVI. Edinburgh, United Kingdom: SPE, 2014.

    [12] Fritts H C. Tree rings and climate. Scientific American, 1972, 226(5): 92- 100.

    [13] Cook E R, Kairiukstis L A. Methods of Dendrochronology. Netherlands: Springer, 1990: 224- 384.

    [15] Ogle K, Whitham T G, Cobb N S. Tree-ring variation in pinyon predicts likelihood of death following severe drought. Ecology, 2000, 81(11): 3237- 3243.

    [16] Hogg E H, Brandt J P, Michaelian M. Impacts of a regional drought on the productivity, dieback, and biomass of western Canadian aspen forests. Canadian Journal of Forest Research, 2008, 38(6): 1373- 1384.

    [17] Mundo I A, El Mujtar V A, Perdomo M H, Gallo L A, Villalba R, Barrera M D. Austrocedrus chilensis growth decline in relation to drought events in northern Patagonia, Argentina. Trees, 2010, 24(3): 561- 570.

    [18] Amoroso M M, Daniels L D, Larson B C. Temporal patterns of radial growth in decliningAustrocedruschilensisforests in Northern Patagonia: the use of tree-rings as an indicator of forest decline. Forest Ecology and Management, 2012, 265: 62- 70.

    [19] Linares J C, Camarero J J. From pattern to process: linking intrinsic water-use efficiency to drought-induced forest decline. Global Change Biology, 2012, 18(3): 1000- 1015.

    [20] Williams A P, Allen C D, Macalady A K, Griffin D, Woodhouse C A, Meko D M, Swetnam T W, Rauscher S A, Seager R, Grissino-Mayer H D, Dean J S, Cook E R, Gangodagamage C, Cai M, McDowell N G. Temperature as a potent driver of regional forest drought stress and tree mortality. Nature Climate Change, 2013, 3(3): 292- 297.

    [21] 趙遵田, 王錫華, 李京東, 肖素榮. 山東省蒙山種子植物區(qū)系研究. 山東科學(xué), 2005, 18(4): 42- 47, 51- 51.

    [22] 高遠(yuǎn), 慈海鑫, 邱振魯, 陳玉峰. 山東蒙山植物多樣性及其海拔梯度格局. 生態(tài)學(xué)報(bào), 2009, 29(12): 6377- 6384.

    [23] Schweingruber F H. Trees and Wood in Dendrochronology: Morphological, Anatomical, and Tree-Ring Analytical Characteristics of Trees Frequently Used in Dendrochronology. Johnson S, trans. Berlin Heidelberg: Springer-Verlag, 1993: 117- 120.

    [24] Rinn F. TSAP-Win Time Series Analysis and Presentation for Dendrochronology and Related Applications, Version 0.53 for Microsoft Windows. Heidelberg, Germany: Rinn Tech, 2003

    [25] Stokes M A, Smiley T L. An Introduction to Tree-Ring Dating. Tucson, AZ, USA: The University of Arizona Press, 1996: 402- 404.

    [26] Cook E R. A time-series analysis approach to tree-ring standardization[D]. Tucson, AZ, USA: University of Arizona, 1985.

    [27] Veblen T T, Hadley K S, Reid M S, Rebertus A J. Methods of detecting past spruce beetle outbreaks in Rocky Mountain subalpine forests. Canadian Journal of Forest Research, 1991, 21(21): 242- 254.

    [28] Didan K. MOD13Q1: MODIS/Terra vegetation indices 16-Day L3 global 250m grid SIN V006. NASA EOSDIS Land Processes DAAC. 2015. https://doi.org/10.5067/MODIS/MOD13Q1.006.

    [29] 王乾坤, 于信芳, 舒清態(tài), 尚珂, 文可戈. MODIS EVI時(shí)序數(shù)據(jù)重建方法及擬合分析. 地球信息科學(xué)學(xué)報(bào), 2015, 17(6): 732- 741.

    [30] Savitzky A, Golay M J E. Smoothing and differentiation of data by simplified least squares procedures. Analytical Chemistry, 1964, 36(8): 1627- 1639.

    [31] John A G, Christopher B F, Michael L G, Kevin L G, Anne E H, Geeske J, Josep P, Riccardo V, Relationships between Ndvi, Canopy Structure, and Photosynthesis in 3 Californian Vegetation Types. Ecological Applications, 1995(5), 28- 41.

    [32] White M A, Nemani R R. Canopy duration has little influence on annual carbon storage in the deciduous broad leaf forest. Global Change Biology, 2003, 9(7): 967- 972.

    [33] Reed B C. Trend analysis of time-series phenology of north america derived from satellite data. Giscience & Remote Sensing, 2006, 43(1): 24- 38.

    [34] Granger C W J. Investigating causal relations by econometric models and cross-spectral methods. Econometrica, 1969, 37(3): 424- 438.

    [35] Haurwitz M W, Brier G W. A critique of the superposed epoch analysis method: its application to solar-weather relations. Monthly Weather Review, 1981, 109(10): 2074- 2079.

    [36] Swetnam T W, Betancourt J L. Temporal patterns of El Nio/Southern Oscillation-wildfire and teleconnections in the southwestern United States//Diaz H F, Markgraf V, eds. El Nio: Historical and Paleoclimatic Aspects of the Southern Oscillation. Cambridge: Cambridge University Press, 1992: 259- 279.

    [37] Swetnam T W, Baisan C H. Historical fire regime patterns in the southwestern United States since AD 1700//Allen C D, ed. Fire Effects in Southwestern Forests. Proceedings of the Second La Mesa Fire Symposium. Los Alamos, New Mexico, March 29- 31, 1994. Fort Collins, CO, USDA Forest Service, Rocky Mountain Forest and Range Experiment Station, General Technical Report RM-GTR- 286. Fort Collins, CO: USDA Forest Service, 1981: 11- 32.

    [38] Tovar C L M. Case Studies on Measuring and Assessing Forest Degradation, Analysis of the Normalized Differential Vegetation Index (NDVI) for the Detection of Degradation of Forest Cover in Mexico 2008- 2009. Forest Resources Assessment Working Paper 173, Rome, Italy: Food and Agriculture Organization of the United Nations, 2009.

    [39] Lopatin E, Kolstr?m T, Spiecker H. Determination of forest growth trends in Komi Republic (northwestern Russia): combination of tree-ring analysis and remote sensing data. Boreal Environment Research, 2006, 11(5): 341- 353.

    [40] Bunn A G, Hughes M K, Kirdyanov A V, Losleben M, Shishov V V, Berner L T, Oltchev A, Vaganov E A. Comparing forest measurements from tree rings and a space-based index of vegetation activity in Siberia. Environmental Research Letters, 2013, 8(3): 035034.

    [41] Cook E R, Johnson A H, Blasing T J. Forest decline: modeling the effect of climate in tree rings. Tree Physiology, 1987, 3(1): 27- 40.

    [42] Johnson A H, Cook E R, Siccama T G. Climate and red spruce growth and decline in the northern appalachians. Proceedings of the National Academy of Sciences of the United States of America, 1988, 85(15): 5369- 5373.

    [43] Beck P S A, Juday G P, Alix C, Barber V A, Winslow S E, Sousa E E, Heiser P, Herriges J D, Goetz S J. Changes in forest productivity across Alaska consistent with biome shift. Ecology Letters, 2011, 14(4): 373- 379.

    [44] Rossi S, Deslauriers A, Anfodillo T, Morin H, Saracino A, Motta R, Borghetti M. Conifers in cold environments synchronize maximum growth rate of tree-ring formation with day length. New Phytologist, 2006, 170(2): 301- 310.

    [45] Vaganov E A, Hughes M K, Shashkin A V. Growth Dynamics of Conifer Tree Rings: Images of Past and Future Environments. Berlin Heidelberg: Springer, 2006: 112.

    [46] Hunt E R Jr, Martin F C, Running S W. Simulating the effects of climatic variation on stem carbon accumulation of a ponderosa pine stand: comparison with annual growth increment data. Tree Physiology, 1991, 9(1/2): 161- 171.

    [47] D′Arrigo R D, Malmstrom C M, Jacoby G C, Los S O, Bunker D E. Correlation between maximum latewood density of annual tree rings and NDVI based estimates of forest productivity. International Journal of Remote Sensing, 2000, 21(11): 2329- 2336.

    [48] McNeil B E, de Beurs K M, Eshleman K N, Foster J R, Townsend P A. Maintenance of ecosystem nitrogen limitation by ephemeral forest disturbance: an assessment using MODIS, Hyperion, and Landsat ETM+. Geophysical Research Letters, 2007, 34(19): L19406.

    [49] Langstr?m B, Annila E, Hellqvist C, Varama M, Niemel? P. Tree mortality, needle biomass recovery and growth losses in scots pine following defoliation byDiprionpini(L.) and subsequent attack byTomicuspiniperda(L.). Scandinavian Journal of Forest Research, 2001, 16(4): 342- 353.

    [51] Vennetier M, Girard F, Taugourdeau O, Cailleret M, Caraglio Y, Sabatier S A, Ouarmim S, Didier C, Thabeet A. Climate change impact on tree architectural development and leaf Area//Singh B R, ed. Climate Change-Realities, Impacts Over Ice Cap, Sea Level and Risks. Rijeka, Croatia: InTech, 2013: 103- 126.

    RelationshipsbetweenforestcanopygreennessandtreeradialgrowthintheMengshanMountainsofShandongProvince

    WANG Zhou1, 2, FANG Ouya1,*

    1LaboratoryofVegetationandEnvironmentalChange,InstituteofBotany,ChineseAcademyofSciences,Beijing100093,China2UniversityofChineseAcademyofSciences,Beijing100049,China

    Forest canopy greenness and tree-ring widths are important indices describing the processes of forest growth. Questions about the relationships between the two indices and their temporal stability are still not completely answered. The canopy greenness index calculated from remote sensing maps is continuous in space, whereas tree-ring width is a comprehensive index that represents the health of the sample sites. Research concerning the relationship between canopy greenness and tree-ring width could enhance the understanding of growth conditions from multi-perspectives and help transform growth conditions of forests at different spatial scales. We collected increment cores from aPinusdensifloraforest at four sample sites with different health conditions in the Mengshan Mountains of Shandong Province and analyzed the relationships between enhanced vegetation indices and tree-ring widths. We found that 1) the variation in tree-ring width was significantly caused by the variation in forest canopy greenness in April and June for healthy forests as calculated using the Granger causality analysis, whereas the relationship was complicated for unhealthy forests. April and June represent the green-up period of the growth season, and according to the results, this green-up period was closely related to radical growth; 2) canopy greenness in other months and radial growth did not have a causality relationship, rather they were both driven by other environmental factors, especially climatic factors; and 3) we used superposed epoch analysis to determine the influence of canopy greenness decrease events on radical growth, and radical growth recovery following the canopy greenness decrease events. Moderate canopy greenness decreases were not always accompanied by a decrease in radial growth; however, if radial growth had already started to decline before the severe greenness loss, then the decrease would continue for five years after canopy greenness decrease event. Radial growth experienced difficult recovery after moderate canopy greenness decrease events. We identified two or three severe canopy greenness decrease events during 2000 to 2014. After these events, radial growth significantly decreased and continued to decrease for years. Radial growth was not completely recovered five years after severe canopy greenness decrease events. The results showed that changes in forest canopy greenness could not represent the beginning of forest stress or the process of forest restoration. Correlations between canopy greenness and tree radial growth existed only in healthy forests. Declines in canopy greenness strongly influenced radical growth and was difficult to recover after losses in greenness.

    tree-ring widths (TRI); enhanced vegetation index (EVI); Granger causality analysis; superposed epoch analysis; forest decline;Pinusdensifloraforest

    國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(31330015)

    2016- 09- 27;

    2017- 04- 11

    *通訊作者Corresponding author.E-mail: oyfang@ibcas.ac.cn

    10.5846/stxb201609271955

    汪舟, 方歐婭.山東蒙山森林冠層綠度與樹干徑向生長(zhǎng)的關(guān)系.生態(tài)學(xué)報(bào),2017,37(22):7514- 7527.

    Wang Z, Fang O Y.Relationships between forest canopy greenness and tree radial growth in the Mengshan Mountains of Shandong Province.Acta Ecologica Sinica,2017,37(22):7514- 7527.

    猜你喜歡
    綠度植被指數(shù)年輪
    寧夏草地綠度時(shí)空變化圖譜及驅(qū)動(dòng)因素分析
    全球森林綠度異常對(duì)干旱事件的響應(yīng)
    2003~2018年干旱對(duì)中國(guó)森林綠度異常的影響
    年輪
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    年輪
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    淺談唐卡中綠度母的造型特點(diǎn)
    心事
    滇池(2017年5期)2017-05-19 03:22:18
    為什么樹有年輪
    a在线观看视频网站| 国产精品久久久av美女十八| 女警被强在线播放| 天天添夜夜摸| 一级a爱片免费观看的视频| 欧美av亚洲av综合av国产av| 91老司机精品| 国内精品久久久久久久电影| 免费大片18禁| 婷婷六月久久综合丁香| 999精品在线视频| 香蕉久久夜色| 成熟少妇高潮喷水视频| 亚洲五月婷婷丁香| 亚洲熟妇熟女久久| 色在线成人网| 黄色成人免费大全| 窝窝影院91人妻| av黄色大香蕉| 黄片小视频在线播放| 日本成人三级电影网站| 曰老女人黄片| 亚洲国产日韩欧美精品在线观看 | 欧美午夜高清在线| 亚洲第一欧美日韩一区二区三区| 日韩有码中文字幕| 久久精品91无色码中文字幕| 成年女人看的毛片在线观看| 麻豆一二三区av精品| 午夜激情福利司机影院| 97人妻精品一区二区三区麻豆| 曰老女人黄片| 手机成人av网站| 欧美国产日韩亚洲一区| 无遮挡黄片免费观看| 国产精品1区2区在线观看.| 一级毛片高清免费大全| 亚洲精品456在线播放app | 噜噜噜噜噜久久久久久91| 久久国产乱子伦精品免费另类| 免费无遮挡裸体视频| 国产熟女xx| 老司机午夜十八禁免费视频| 男女之事视频高清在线观看| 亚洲美女视频黄频| 亚洲熟女毛片儿| 中文资源天堂在线| 日韩欧美国产在线观看| 操出白浆在线播放| 18禁美女被吸乳视频| 宅男免费午夜| 在线国产一区二区在线| 久久久久国产一级毛片高清牌| 别揉我奶头~嗯~啊~动态视频| 欧美成狂野欧美在线观看| 免费看十八禁软件| 一夜夜www| 国产高清视频在线观看网站| 欧美一级a爱片免费观看看| 午夜福利高清视频| 女人被狂操c到高潮| 日本免费a在线| 舔av片在线| 国产不卡一卡二| 激情在线观看视频在线高清| 一区二区三区国产精品乱码| 国产精品久久电影中文字幕| 非洲黑人性xxxx精品又粗又长| 国产成人精品无人区| 天天添夜夜摸| 精品久久久久久,| 又黄又爽又免费观看的视频| 一区二区三区国产精品乱码| 国产又黄又爽又无遮挡在线| 国产69精品久久久久777片 | 在线国产一区二区在线| 欧美大码av| 少妇的丰满在线观看| 成人av在线播放网站| 免费人成视频x8x8入口观看| 成年女人永久免费观看视频| 国产91精品成人一区二区三区| 色在线成人网| 伊人久久大香线蕉亚洲五| 国语自产精品视频在线第100页| 国产欧美日韩精品一区二区| 欧美中文日本在线观看视频| 成人av在线播放网站| 日本一二三区视频观看| av视频在线观看入口| 国产午夜福利久久久久久| 校园春色视频在线观看| 很黄的视频免费| 国产精品美女特级片免费视频播放器 | 亚洲色图av天堂| 欧美精品啪啪一区二区三区| av欧美777| 伊人久久大香线蕉亚洲五| 男人舔女人的私密视频| 久久久久久人人人人人| 我的老师免费观看完整版| 日本一二三区视频观看| 国产精品爽爽va在线观看网站| 九九热线精品视视频播放| 噜噜噜噜噜久久久久久91| 黄片小视频在线播放| 欧美av亚洲av综合av国产av| 久久久久久人人人人人| 国产精品99久久99久久久不卡| 不卡av一区二区三区| 亚洲av成人一区二区三| 国产精品电影一区二区三区| 在线播放国产精品三级| 国产精品亚洲av一区麻豆| 热99在线观看视频| 婷婷精品国产亚洲av| 一本综合久久免费| 黑人操中国人逼视频| 窝窝影院91人妻| 天堂网av新在线| xxx96com| 久久久久久久久免费视频了| 俺也久久电影网| 天堂√8在线中文| av视频在线观看入口| 日日干狠狠操夜夜爽| 亚洲av电影在线进入| 叶爱在线成人免费视频播放| 日韩欧美三级三区| АⅤ资源中文在线天堂| 精品欧美国产一区二区三| 国产成人系列免费观看| 国产精品乱码一区二三区的特点| 香蕉久久夜色| 精品乱码久久久久久99久播| 国产私拍福利视频在线观看| 在线视频色国产色| 国产成人系列免费观看| 搡老岳熟女国产| 成人三级黄色视频| 国产高潮美女av| 一级a爱片免费观看的视频| 国产又黄又爽又无遮挡在线| 首页视频小说图片口味搜索| 高清毛片免费观看视频网站| 欧美zozozo另类| 日本a在线网址| 一个人免费在线观看电影 | 欧美中文日本在线观看视频| 欧美3d第一页| 日本撒尿小便嘘嘘汇集6| 日本撒尿小便嘘嘘汇集6| 色视频www国产| 少妇的丰满在线观看| 99久久综合精品五月天人人| 最近在线观看免费完整版| 在线观看日韩欧美| 成人av在线播放网站| 亚洲人与动物交配视频| 欧美黑人巨大hd| 日韩高清综合在线| 国产精华一区二区三区| 欧美日韩一级在线毛片| 色在线成人网| 九九久久精品国产亚洲av麻豆 | 日本 欧美在线| 在线观看日韩欧美| 国产69精品久久久久777片 | 在线永久观看黄色视频| 中文字幕人成人乱码亚洲影| 亚洲国产精品sss在线观看| 丁香欧美五月| АⅤ资源中文在线天堂| 国产精品日韩av在线免费观看| 黄频高清免费视频| 精品久久久久久久毛片微露脸| 亚洲中文av在线| 亚洲狠狠婷婷综合久久图片| 深夜精品福利| 岛国视频午夜一区免费看| 亚洲精品色激情综合| 99久久精品热视频| 日本免费一区二区三区高清不卡| 国产午夜福利久久久久久| 一二三四在线观看免费中文在| 欧美日韩瑟瑟在线播放| 欧美日韩瑟瑟在线播放| 欧美绝顶高潮抽搐喷水| 亚洲无线观看免费| 99久久99久久久精品蜜桃| 欧美乱码精品一区二区三区| 好看av亚洲va欧美ⅴa在| 国产熟女xx| 一个人看视频在线观看www免费 | 夜夜夜夜夜久久久久| 久久这里只有精品19| 欧美日韩一级在线毛片| 非洲黑人性xxxx精品又粗又长| 夜夜夜夜夜久久久久| 老司机福利观看| 一本一本综合久久| 久久九九热精品免费| 一本久久中文字幕| 99久久精品国产亚洲精品| 桃色一区二区三区在线观看| 老汉色av国产亚洲站长工具| 亚洲av成人不卡在线观看播放网| 日本熟妇午夜| 国产精品久久电影中文字幕| 1024手机看黄色片| 非洲黑人性xxxx精品又粗又长| 日本撒尿小便嘘嘘汇集6| 久久人妻av系列| 欧美极品一区二区三区四区| 久久精品国产亚洲av香蕉五月| 国产免费av片在线观看野外av| 国产蜜桃级精品一区二区三区| 十八禁网站免费在线| 舔av片在线| 母亲3免费完整高清在线观看| 午夜精品在线福利| 99久久精品国产亚洲精品| 看黄色毛片网站| 国语自产精品视频在线第100页| 看黄色毛片网站| 超碰成人久久| 欧美日韩精品网址| 久久人人精品亚洲av| 老汉色av国产亚洲站长工具| 久久99热这里只有精品18| 欧美黄色片欧美黄色片| 99热只有精品国产| 天天躁日日操中文字幕| 最近在线观看免费完整版| 欧美国产日韩亚洲一区| 天堂av国产一区二区熟女人妻| 日本成人三级电影网站| 久久亚洲精品不卡| 免费av毛片视频| 久久这里只有精品中国| 日韩人妻高清精品专区| 日韩欧美一区二区三区在线观看| 九九热线精品视视频播放| 国产亚洲av嫩草精品影院| 女人被狂操c到高潮| 国产亚洲欧美98| 成人三级做爰电影| 黑人欧美特级aaaaaa片| 99热精品在线国产| 一本综合久久免费| 精品电影一区二区在线| 亚洲激情在线av| 欧美激情在线99| 一夜夜www| 日本一二三区视频观看| 在线免费观看不下载黄p国产 | a级毛片a级免费在线| 色精品久久人妻99蜜桃| 美女 人体艺术 gogo| 国产单亲对白刺激| 美女扒开内裤让男人捅视频| 亚洲av成人精品一区久久| 国产精品99久久久久久久久| 精品国产亚洲在线| 国产野战对白在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区激情短视频| 夜夜爽天天搞| x7x7x7水蜜桃| 亚洲人成网站高清观看| 国产精品乱码一区二三区的特点| 国产又色又爽无遮挡免费看| 久久精品aⅴ一区二区三区四区| 51午夜福利影视在线观看| 国产高清视频在线观看网站| 久久这里只有精品中国| 国产精品日韩av在线免费观看| 一级作爱视频免费观看| 精品久久久久久久久久久久久| 免费看光身美女| 亚洲乱码一区二区免费版| 国产精品免费一区二区三区在线| 亚洲精华国产精华精| 亚洲第一电影网av| 啦啦啦观看免费观看视频高清| 日本三级黄在线观看| 97超级碰碰碰精品色视频在线观看| 日韩有码中文字幕| 99久国产av精品| 熟女少妇亚洲综合色aaa.| 精品福利观看| 欧美午夜高清在线| 亚洲无线在线观看| 欧美成狂野欧美在线观看| 色av中文字幕| 天堂影院成人在线观看| 亚洲成av人片在线播放无| 国产熟女xx| 老司机在亚洲福利影院| 久久久久亚洲av毛片大全| 亚洲狠狠婷婷综合久久图片| 后天国语完整版免费观看| 国产成人欧美在线观看| 法律面前人人平等表现在哪些方面| xxx96com| 免费无遮挡裸体视频| 欧美日本亚洲视频在线播放| 老汉色av国产亚洲站长工具| 中亚洲国语对白在线视频| 中文字幕精品亚洲无线码一区| 亚洲人成网站高清观看| 精品福利观看| 国产精品久久久av美女十八| 久久久国产精品麻豆| 午夜精品在线福利| 精品久久久久久久久久久久久| 日韩国内少妇激情av| 婷婷精品国产亚洲av在线| 精品乱码久久久久久99久播| 欧美zozozo另类| 男插女下体视频免费在线播放| 欧美xxxx黑人xx丫x性爽| 欧美在线一区亚洲| 免费在线观看亚洲国产| 日日摸夜夜添夜夜添小说| 老司机在亚洲福利影院| 欧美在线黄色| 51午夜福利影视在线观看| 两人在一起打扑克的视频| 亚洲av成人一区二区三| 日本一本二区三区精品| 一卡2卡三卡四卡精品乱码亚洲| 99久国产av精品| 99久久99久久久精品蜜桃| 国产乱人视频| 免费电影在线观看免费观看| 露出奶头的视频| 国产精华一区二区三区| 熟女少妇亚洲综合色aaa.| 国内久久婷婷六月综合欲色啪| 色老头精品视频在线观看| 午夜福利18| 19禁男女啪啪无遮挡网站| 麻豆av在线久日| 精品一区二区三区视频在线 | www.999成人在线观看| 久久精品aⅴ一区二区三区四区| 最近视频中文字幕2019在线8| 亚洲色图 男人天堂 中文字幕| 国产久久久一区二区三区| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 亚洲无线在线观看| 国产麻豆成人av免费视频| 可以在线观看毛片的网站| 成人av一区二区三区在线看| 99精品欧美一区二区三区四区| 久久天躁狠狠躁夜夜2o2o| 国产久久久一区二区三区| 白带黄色成豆腐渣| 精华霜和精华液先用哪个| 国产精品女同一区二区软件 | 嫁个100分男人电影在线观看| 桃色一区二区三区在线观看| netflix在线观看网站| 又粗又爽又猛毛片免费看| 亚洲狠狠婷婷综合久久图片| 国内精品久久久久久久电影| 看黄色毛片网站| 国产精品影院久久| 国产精品自产拍在线观看55亚洲| 国产日本99.免费观看| 国产精品99久久久久久久久| www日本黄色视频网| 久久久久久大精品| 无人区码免费观看不卡| 国产1区2区3区精品| 少妇人妻一区二区三区视频| 亚洲五月天丁香| 怎么达到女性高潮| 2021天堂中文幕一二区在线观| 久久久久久久精品吃奶| www.熟女人妻精品国产| 老汉色∧v一级毛片| or卡值多少钱| 18禁裸乳无遮挡免费网站照片| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品久久男人天堂| 又黄又爽又免费观看的视频| 91av网一区二区| 亚洲精品久久国产高清桃花| 免费在线观看亚洲国产| 国产精品99久久久久久久久| 巨乳人妻的诱惑在线观看| 99久久无色码亚洲精品果冻| 精品午夜福利视频在线观看一区| 国产午夜精品论理片| 欧美一级a爱片免费观看看| 亚洲精品一卡2卡三卡4卡5卡| 在线观看一区二区三区| 欧美日韩综合久久久久久 | 日韩大尺度精品在线看网址| 夜夜看夜夜爽夜夜摸| av国产免费在线观看| 亚洲国产欧美一区二区综合| 亚洲天堂国产精品一区在线| 亚洲中文字幕一区二区三区有码在线看 | 久久久久久久午夜电影| 久99久视频精品免费| 免费一级毛片在线播放高清视频| 精品人妻1区二区| 国产伦人伦偷精品视频| 婷婷丁香在线五月| 又紧又爽又黄一区二区| 亚洲熟妇中文字幕五十中出| 18禁黄网站禁片午夜丰满| 国产伦在线观看视频一区| 99re在线观看精品视频| 久久久久国内视频| 啦啦啦观看免费观看视频高清| 国产精品久久久久久亚洲av鲁大| 国产黄色小视频在线观看| 色尼玛亚洲综合影院| 97碰自拍视频| 露出奶头的视频| 精品久久蜜臀av无| 免费无遮挡裸体视频| 久久久久久久午夜电影| 国产aⅴ精品一区二区三区波| 神马国产精品三级电影在线观看| 麻豆av在线久日| 人人妻人人看人人澡| 黄色女人牲交| 久99久视频精品免费| 国产黄片美女视频| 亚洲性夜色夜夜综合| 法律面前人人平等表现在哪些方面| 亚洲精品一卡2卡三卡4卡5卡| 免费高清视频大片| av女优亚洲男人天堂 | 美女黄网站色视频| 日本黄色片子视频| 国产极品精品免费视频能看的| 国产精品精品国产色婷婷| 国产精品亚洲av一区麻豆| 中文字幕久久专区| 男人舔女人下体高潮全视频| 俺也久久电影网| 一区二区三区国产精品乱码| 丁香欧美五月| 国产亚洲av嫩草精品影院| 91九色精品人成在线观看| 97超视频在线观看视频| 特大巨黑吊av在线直播| 老熟妇仑乱视频hdxx| 淫秽高清视频在线观看| 久久久久国产一级毛片高清牌| 亚洲国产高清在线一区二区三| 亚洲狠狠婷婷综合久久图片| 午夜成年电影在线免费观看| 国产高潮美女av| 欧美成人性av电影在线观看| 国产在线精品亚洲第一网站| 日韩精品青青久久久久久| 亚洲黑人精品在线| 美女被艹到高潮喷水动态| 一级a爱片免费观看的视频| 在线免费观看不下载黄p国产 | 亚洲精品456在线播放app | 国产精品久久久av美女十八| 操出白浆在线播放| 熟女电影av网| 国产成人精品久久二区二区免费| 久久这里只有精品中国| 桃色一区二区三区在线观看| 亚洲自偷自拍图片 自拍| 成人三级黄色视频| 99精品久久久久人妻精品| 黑人操中国人逼视频| av天堂在线播放| 亚洲色图 男人天堂 中文字幕| 激情在线观看视频在线高清| 淫秽高清视频在线观看| 最近最新中文字幕大全电影3| 巨乳人妻的诱惑在线观看| 无人区码免费观看不卡| 在线a可以看的网站| 最近视频中文字幕2019在线8| 国内精品久久久久精免费| 久久久成人免费电影| 欧美一级毛片孕妇| 神马国产精品三级电影在线观看| 两性夫妻黄色片| 在线观看66精品国产| 欧美性猛交╳xxx乱大交人| 一本精品99久久精品77| 午夜福利18| 男人和女人高潮做爰伦理| 超碰成人久久| 亚洲中文av在线| 在线播放国产精品三级| 亚洲精品一区av在线观看| 男插女下体视频免费在线播放| 国产免费男女视频| 国产成年人精品一区二区| 欧美3d第一页| 国产高清视频在线播放一区| 久久亚洲真实| 欧美精品啪啪一区二区三区| 国产一区二区在线av高清观看| 色尼玛亚洲综合影院| 又粗又爽又猛毛片免费看| 午夜精品在线福利| 精品不卡国产一区二区三区| 亚洲av第一区精品v没综合| 国产熟女xx| 久久精品综合一区二区三区| 脱女人内裤的视频| 日本精品一区二区三区蜜桃| 日韩欧美国产在线观看| 国产高清激情床上av| 午夜福利成人在线免费观看| 亚洲国产精品sss在线观看| 久久久久性生活片| 1000部很黄的大片| 母亲3免费完整高清在线观看| 国产探花在线观看一区二区| 在线免费观看的www视频| 一卡2卡三卡四卡精品乱码亚洲| 搡老岳熟女国产| 欧美丝袜亚洲另类 | 免费一级毛片在线播放高清视频| 精品国产三级普通话版| www国产在线视频色| 国产三级在线视频| 亚洲欧美激情综合另类| av在线蜜桃| 性色avwww在线观看| 免费无遮挡裸体视频| 亚洲av熟女| 久久久久国产精品人妻aⅴ院| 午夜亚洲福利在线播放| 中文资源天堂在线| 国产精品自产拍在线观看55亚洲| 中文字幕av在线有码专区| 中文字幕人妻丝袜一区二区| 黄色丝袜av网址大全| 国产真人三级小视频在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲成人久久性| 免费观看的影片在线观看| 欧美不卡视频在线免费观看| 精品国产乱码久久久久久男人| 亚洲av中文字字幕乱码综合| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 日韩 欧美 亚洲 中文字幕| 成熟少妇高潮喷水视频| avwww免费| 法律面前人人平等表现在哪些方面| 国产不卡一卡二| 亚洲av电影不卡..在线观看| 国产不卡一卡二| 校园春色视频在线观看| 久久亚洲精品不卡| 久久伊人香网站| 色综合亚洲欧美另类图片| 国内久久婷婷六月综合欲色啪| 久久婷婷人人爽人人干人人爱| 亚洲激情在线av| 一卡2卡三卡四卡精品乱码亚洲| 18禁黄网站禁片免费观看直播| 国产成人欧美在线观看| 久久久国产欧美日韩av| 又爽又黄无遮挡网站| 日本与韩国留学比较| 国内精品久久久久久久电影| 好男人电影高清在线观看| 熟女人妻精品中文字幕| 亚洲 国产 在线| 很黄的视频免费| 欧美日本亚洲视频在线播放| 欧美黄色淫秽网站| 欧美av亚洲av综合av国产av| 亚洲欧美精品综合久久99| 亚洲乱码一区二区免费版| 欧美色欧美亚洲另类二区| 成人国产一区最新在线观看| 国产成人aa在线观看| 大型黄色视频在线免费观看| 在线观看免费午夜福利视频| 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| av在线蜜桃| 午夜精品久久久久久毛片777| 12—13女人毛片做爰片一| 国产淫片久久久久久久久 | 日韩欧美一区二区三区在线观看| 高清在线国产一区| 亚洲激情在线av| 最近最新中文字幕大全电影3| 在线免费观看的www视频| 久久久久久久久免费视频了| 久久天堂一区二区三区四区| xxx96com| 哪里可以看免费的av片| 看片在线看免费视频| 国产精品爽爽va在线观看网站| 久久精品亚洲精品国产色婷小说| 岛国在线免费视频观看| 在线观看66精品国产| 免费观看的影片在线观看| 国产高清视频在线播放一区| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频|