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

    山西省長治市過去150年森林的生長抑制和釋放歷史

    2017-06-22 14:04:16梁寒雪
    生態(tài)學(xué)報 2017年9期
    關(guān)鍵詞:樹輪年輪樹木

    張 啟,閆 明,梁寒雪

    1 山西師范大學(xué)生命科學(xué)學(xué)院,臨汾 041004 2 中國科學(xué)院植物研究所植被與環(huán)境變化國家重點實驗室,北京 100093 3 中國科學(xué)院大學(xué),北京 100049

    ?

    山西省長治市過去150年森林的生長抑制和釋放歷史

    張 啟1,閆 明1,梁寒雪2,3,*

    1 山西師范大學(xué)生命科學(xué)學(xué)院,臨汾 041004 2 中國科學(xué)院植物研究所植被與環(huán)境變化國家重點實驗室,北京 100093 3 中國科學(xué)院大學(xué),北京 100049

    森林生長與氣候變化有著緊密的關(guān)系,在全球變暖情形下了解樹木的干擾歷史對準確預(yù)測森林生長的變化具有促進作用。本文選擇山西黃土高原東南部長治地區(qū)保存較好的一個油松(Pinustabuliformis)和兩個白皮松(Pinusbungeana)森林為研究對象,利用樹木年輪學(xué)方法分析了其干擾歷史。結(jié)果表明:黎城縣、平順縣和屯留縣研究地點中年齡最老的樹木分別有227、185、102a;通過計算樹木徑向生長的變化幅度,發(fā)現(xiàn)該地區(qū)在過去150年中發(fā)生了3次大的生長抑制事件(分別發(fā)生在1873—1877、1925—1930和1994—1997年期間)和5次大的生長釋放事件(分別發(fā)生在1867—1871、1878—1884、1930—1935、1980—1985和1999—2004年期間)。樹木徑向生長與氣象觀測資料的相關(guān)分析顯示,該地區(qū)森林生長在年際尺度上主要與6月份溫度呈負相關(guān)而與4—5月份水分條件呈正相關(guān),揭示了由降水減少或高溫下水分蒸散所帶來的極端干旱事件是導(dǎo)致森林抑制現(xiàn)象的主要原因。這些森林歷史研究結(jié)果對區(qū)域林業(yè)管理具有實踐意義,營林建設(shè)更宜選擇在森林生長釋放時期;經(jīng)歷多次干擾而存留下來的老齡樹具有較強的抗干擾能力和豐富的歷史信息,加強對老齡樹的鑒定與針對性保護可有效維護區(qū)域森林生態(tài)服務(wù)功能。

    黃土高原;長治;樹木年輪;生長抑制;生長釋放

    干擾事件會對樹木生長產(chǎn)生影響,使樹木表現(xiàn)出生長抑制和釋放兩類不同的生長變化。徑向生長率顯著低于平均生長率的現(xiàn)象稱作生長抑制,而徑向生長率出現(xiàn)快速而持續(xù)增加的現(xiàn)象則被稱作生長釋放[1]。研究樹木生長抑制和釋放的規(guī)律和發(fā)生機制,可以推測過去的干擾歷史,為預(yù)測未來的干擾提供重要的依據(jù),對理解森林生長動態(tài)和森林保護至關(guān)重要。樹木年輪由于能準確記錄年際間的生長變化,具有重建歷史時間長、定年準確、材料廣泛容易獲得的優(yōu)點,所以可以有效地用于研究森林的釋放和抑制歷史[2]。目前國內(nèi)對森林的生長抑制和釋放已有一些相關(guān)研究,在不同的研究中確定生長抑制和釋放的方法存在差異,生長釋放和抑制事件的成因也因地區(qū)、樹種的不同存在不同的解釋[3- 6]。

    黃土高原東南部是中國東部季風區(qū)向西部干旱區(qū)的過渡地帶,既是氣候變化敏感區(qū),又是生態(tài)環(huán)境脆弱帶,還是黃河中上游水土保持重點區(qū)域[7]。20世紀以來,這一地區(qū)暖干化趨勢十分明顯,氣象災(zāi)害發(fā)生頻繁,對樹木的生長構(gòu)成極大威脅[8]。森林對于改善環(huán)境,涵養(yǎng)水源,防止水土流失的作用非常重要,保護森林植被是黃土高原治理的治本之策。對該地區(qū)森林生長干擾歷史的研究,有助于了解樹木的生長方式以及森林植被動態(tài)規(guī)律,對森林的保護也有參考意義。然而目前山西黃土高原上的樹木年輪學(xué)研究主要集中在重建過去的氣候變化與樹輪氣候響應(yīng)上[9- 12],對于森林干擾歷史的相關(guān)研究則很少。

    山西長治地區(qū)位于黃土高原東南部,是東亞季風影響的邊緣地帶。由于自然氣候變化和人類活動的影響,該地區(qū)的森林受到了嚴重的干擾,森林資源在逐年減少。本文在長治地區(qū)保存較好的一個油松(Pinustabuliformis)和兩個白皮松(Pinusbungeana)森林中采集了樹木年輪樣本,旨在研究以下兩個問題:1)該地區(qū)森林的年齡有多大,在其歷史上經(jīng)歷了哪些抑制和釋放干擾事件?2)造成該地區(qū)森林生長釋放和抑制的因素可能是什么?研究結(jié)果可為林業(yè)部門制定合理的經(jīng)營管理方案,提供森林生長歷史數(shù)據(jù),對該地區(qū)森林的保護也有參考意義。

    1 研究區(qū)概況和研究方法

    1.1 研究區(qū)概況

    長治地區(qū)位于中國山西省東南部,海拔最高處為沁源太岳山2453m,最低處為平順漳河谷380m[13]。長治地區(qū)植被類型多樣,屬于華北暖溫帶闊葉林地帶,植被垂直分布明顯,從下到上依次為:灌叢(或灌草叢)帶、低中山針葉林(或針闊混交林)帶、落葉闊葉林帶、亞高山草甸帶。土壤類型主為褐土、褐土性土、山地草甸土、棕壤[13]。該地區(qū)屬于典型中溫帶半濕潤大陸性季風氣候,全年冬無嚴寒,夏無酷暑,雨熱同季。年平均氣溫在4.9—10.4℃之間。本研究中的氣候數(shù)據(jù)來自長治市氣象站(36.13°N,113.31°E),多年的統(tǒng)計結(jié)果表明該地區(qū)1月份最冷,平均最低氣溫為-6.9℃;7月份最熱,平均最高氣溫為22.5℃。年日照時數(shù)2418—2616h,一般年降水量在537.4—656.7mm,7月最多,為132.2mm,1月最少,為5.5mm,年平均無霜期在156.8—181.9h,年平均風速為1.5—3.0m/s之間[14]。

    1.2 樹輪數(shù)據(jù)

    圖1 采樣點示意圖Fig.1 The location map of the sampling sites

    本研究選取長治地區(qū)保存較好,樹齡相對較長的三個采樣點,分別位于屯留縣紫金山(36.21°N,112.47°E)、平順縣南垴山(36.20°N,113.29°E)、黎城縣廣志山(36.31°N,113.15°E)(圖1)。其中屯留縣位于長治市上黨盆地西緣,屬于太行山系,總土地面積11.9萬hm2,森林覆蓋率28.30%;黎城縣位于長治市東北部,太行山中南段,有林地面積7.68萬hm2,森林覆蓋率達24.05%;平順縣長治市的東部,太行山南端,有林面積94.74萬畝,森林覆蓋率41.6。三個縣城都是山西省園林縣城[15]。采樣于2012年10月底進行,按照國際樹木年輪取樣的標準,我們在3個采樣點選取年齡較老的樹木,在其胸高處(1.3m)沿平行坡向取一根樣芯,共采集到了27棵樹的樣芯。其中黎城縣廣志山油松(Pinustabuliformis)10個樹芯,平順縣南垴山白皮松(Pinusbungeana)11個樹芯,屯留縣紫金山白皮松(Pinusbungeana)6個樹芯。

    樣品帶回實驗室后,按照國際樹木年輪分析的基本程序[16],將樣本進行自然干燥、棉線固定、砂紙打磨,然后利用LINTAB樹輪測量儀進行年輪寬度的測量,精度為0.001mm。通過比較不同樹木個體間的年輪序列,鑒別出缺輪、偽輪和測量誤差,使每棵樹上的每個年輪都被確定了其準確的生長年份。然后利用COFECHA[17]程序?qū)徊娑杲Y(jié)果進行了檢查。

    樹木年輪寬度序列包含了與年齡有關(guān)的生長趨勢以及樹木間的相互競爭等微環(huán)境所導(dǎo)致的其他信號[18- 19]。為了剔除這些信號,用ARSTAN[20]程序?qū)唽挃?shù)據(jù)去趨勢??紤]到樣本年齡較短、樣本量較少,而且本研究更加關(guān)注高頻的信息,故選用了步長為16a的樣條函數(shù)來擬合生長趨勢,去掉了周期較長的信號,再用實際年輪寬度值除以對應(yīng)的生長趨勢線的值得到每個樹輪序列的輪寬指數(shù)。

    樹木生長變化百分率被認為能夠準確判定生長抑制和釋放事件[21- 22]。采用Nowacki和Abrams的方法計算每個樣本輪寬指數(shù)序列的生長變化百分率[23],本文計算的是5a滑動平均值,即GC%=[(M2-M1)/M1]×100%,其中GC%是生長變化百分率,M1和M2分別是前5a(包括當年)和隨后5a(不包括當年)輪寬指數(shù)的平均值。由于生長釋放和抑制事件是一個連續(xù)多年的樹木生長狀態(tài),為了區(qū)別于個別年份的偶然因素所造成的生長變化對生長釋放和抑制事件檢測的影響。以5a的時間為窗口,對得到的GC%序列中的數(shù)值進行篩選。若5年內(nèi)生長變化率平均值大于0.25,則定義該時期發(fā)生了生長釋放事件;若5年內(nèi)生長變化率平均值小于-0.25,定義該時期為生長抑制事件。為了研究整個區(qū)域森林的徑向生長變化情況,又計算出各個樣點所有樹木每年生長變化百分率的中值,得到各樣點的生長變化百分率序列,并分析了樣點序列間的相關(guān)性來說明不同采樣點樹木徑向生長變化的一致性。同時統(tǒng)計了各點樹木發(fā)生生長抑制和釋放的時間段及樹木個數(shù)來分析其歷史狀況。

    本研究使用DendroClim 2002程序[24]分析了各樣點標準年表和氣候因子一階差序列的相關(guān)關(guān)系,以此來進而探討該地區(qū)森林生長在年際尺度上對當?shù)貧夂蜃兓捻憫?yīng)特征。氣候因子選取月平均溫度和月降水量,其中月平均溫度是基于長治氣象觀測站觀測的當月每天的平均氣溫的平均值,月降水量是是指當月每天降水量的總和。氣象數(shù)據(jù)選取公共時間段(1958—2012年),選取前一年10月到當年9月的月平均溫度和月降水量。

    2 研究結(jié)果

    2.1 樹輪樣本及交叉定年結(jié)果

    經(jīng)過交叉定年,獲得了可用于分析的樹輪數(shù)據(jù),COFECHA的檢驗結(jié)果表明定年準確可靠(表1)。其中黎城樣本總體年齡偏大,最大年齡為227a;平順樣本最大年齡為185a;屯留樣本總體年齡相對偏小,最大年齡為102a。各采樣點序列間平均相關(guān)系數(shù)都達到了0.5以上,表明該區(qū)域內(nèi)各點的樹木包含著大量的共同信號。

    表1 3個采樣點樹木年輪樣本的統(tǒng)計特征

    2.2 樹輪生長變化百分率

    經(jīng)計算得到的樹輪生長變化百分率序列如圖2所示,由于黎城1859年和平順1894年以前樣本量都不足3,所以計算結(jié)果分別從1860年和1895作為起始年。結(jié)果顯示3個樣點所有的樹木生長變化百分率的變化趨勢大致是一致的,高生長變化和低生長變化出現(xiàn)的時間段也基本一致。為了研究整體區(qū)域內(nèi)樹輪生長變化的特點,對各采樣點(1915—2007年)生長變化序列在公共時間段內(nèi)進行了相關(guān)分析,相關(guān)系數(shù)均達到了顯著正相關(guān),其中黎城和平順之間的相關(guān)系數(shù)為0.660(P<0.001),黎城和屯留之間的相關(guān)系數(shù)為0.560(P<0.001),平順和屯留之間的相關(guān)系數(shù)為0.564(P<0.001)。

    圖2 3個采樣點樹木生長變化百分率情況Fig.2 Percent growth change of the three sampling sites虛線為單個樹木的生長變化百分率,實線為各點所有樹木的生長變化百分率的中值

    2.3 生長抑制和釋放事件

    從各點樹輪生長變化百分率圖(圖2)及生長抑制和釋放樣本量來看(圖3),黎城所表現(xiàn)出的生長釋放事件年代分布最廣,最早的一次釋放事件出現(xiàn)在1867年,最后一次釋放事件出現(xiàn)在2000年;每10年期間都會有個體發(fā)生釋放,其中個體出現(xiàn)釋放的高峰期主要是5個,分別是1867—1871、1878—1883、1892—1896、1931—1935和1981—1985。平順樣點的樹木生長釋放事件比較集中,幾乎每年都會有個體發(fā)生釋放,但主要有4個高峰期,分別集中在1929—1935、1943—1947、1981—1985和1999—2004年期間;屯留樣點的樹木生長釋放事件主要集中在在1931—1935、1943—1947、1959—1963、1980—1984和2002—2004年間。

    黎城樣點的樹木生長抑制事件多集中在1864—1865、1873—1877、1889—1890、1925—1927和1977年期間。平順所表現(xiàn)的抑制大體上呈現(xiàn)出3個高峰:1923—1927、1938—1942和1994—1997年期間。而屯留生長抑制主要在1925—1928、1938—1941和1964—1966年期間。

    綜合統(tǒng)計了各點樹木發(fā)生生長抑制和釋放的時間段及樹木個數(shù),發(fā)現(xiàn)該地區(qū)森林在過去一百多年期間出現(xiàn)生長抑制的主要時期有3個,分別是1873—1877、1925—1930和1994—1997年期間;生長釋放事件共有5次,分別發(fā)生在1867—1871、1878—1884、1930—1935、1980—1985和1999—2004年期間。

    圖3 各點生長釋放和生長抑制樹木個數(shù)隨時間的變化Fig.3 Number of trees with release and suppression over time in the three sampling sites

    2.4 樹輪與氣候要素的相關(guān)分析

    3個樣點的標準年表與氣候因子(月平均溫度和月總降水量)的一階差相關(guān)結(jié)果如圖4所示。從圖中可以看出,研究區(qū)域樹輪寬度指數(shù)與月平均溫度的相關(guān)結(jié)果一致性較差。3個樣點只與當年6月份平均溫度都顯著負相關(guān),相關(guān)系數(shù)分別為黎城-0.47(n=10,p<0.05)、平順-0.48(n=11,P<0.05)和屯留-0.31(n=6,P<0.05)。其他月份的相關(guān)系數(shù)既有正值也有負值,相關(guān)系數(shù)也比較低。從樹輪與月總降水量的相關(guān)性來看,3個樣點的樹木生長均與當年4—5月份總降水量都達到顯著正相關(guān)。其中與4月份總降水量相關(guān)系數(shù)分別為黎城0.28(n=10,P<0.05)、平順0.36(n=11,P<0.05)、屯留0.43(n=6,P<0.05);與5月份總降水量分別為黎城0.35(n=10,P<0.05),平順0.32(n=11,P<0.05),屯留0.39(n=6,P<0.05)。

    圖4 樹輪徑向生長和氣候因子一階差的相關(guān)關(guān)系Fig.4 Correlation coefficients between the 1st differences of tree ring width chronologies and climate factorsP為上一年月份,C為當年月份,*表示該月份相關(guān)達到95%顯著水平

    3 討論

    3.1 樹木徑向生長對氣候變化的響應(yīng)

    該地區(qū)各點森林生長在年際尺度上與氣候響應(yīng)類似,都對生長季5—6月的平均溫度負相關(guān)和4—5月的降水正相關(guān),其它時期生長對氣候條件的響應(yīng)關(guān)系較弱。在我國黃土高原地區(qū),4—6月是油松開始和快速生長的時期[25],生長季前期過高的溫度加速了樹木的蒸騰量,使樹木體內(nèi)的水分散失過快,限制了樹木的生理代謝活動,不利于樹木的徑向生長,從而導(dǎo)致了較窄的輪寬。在黃土高原的崆峒山[26]、山西蘆芽山[27]、陜西黃龍山[28]等地區(qū)的相關(guān)研究也發(fā)現(xiàn)了類似的結(jié)果。生長季前期(4—5月)充足的降水能夠使樹木儲存足夠的水分供給生長以及養(yǎng)分的運輸,從而導(dǎo)致了樹木較高的生長量。而到了7—8月份是樹木生長旺盛時期,盡管溫度較高,但長治地區(qū)充足的降水能滿足樹木的生長,此時的溫度和降水都不成為它們生長的限制因子[29]。盡管3個樣點樹木生長都和6月的平均溫度顯著負相關(guān)(n=27,r=-0.31—-0.48,P<0.05),與4—5月的降水顯著正相關(guān)(n=27,r=0.28—0.43,P<0.05),但是對樹輪生長的方差解釋量和其他樹輪氣候?qū)W的研究相比系數(shù)較低,說明局部小環(huán)境對樣地樹木生長也產(chǎn)生了一定的影響。

    3.2 森林生長抑制與釋放的特征及解釋

    本研究發(fā)現(xiàn)長治地區(qū)森林的釋放與抑制現(xiàn)象在有些時間段上表現(xiàn)出很強的一致性。例如,3個點的樹木在1923—1928年期間都表現(xiàn)較強的生長抑制,之后直到1936年期間,都存在生長釋放的情況。同樣的在1970—1985年間也存在一致的生長變化。這種區(qū)域森林的抑制和釋放現(xiàn)象的一致性,極有可能是由于氣候事件導(dǎo)致的[6]。但是在個別時期也存在微小的差異,平順和屯留的白皮松在1938—1942期間都存在生長抑制,隨后在1943—1947期間都進入快速生長階段,而黎城的油松在此期間少數(shù)樹木表現(xiàn)出類似的生長變化。屯留的白皮松在1956—1957年間生長降低,之后1959—1963發(fā)生了生長釋放事件,而其它兩個地方?jīng)]有明顯的釋放表現(xiàn),這種情況可能是由于樣點微環(huán)境或者樹木自身生長特性的不同所導(dǎo)致的。

    將結(jié)果與史料記載和近代有氣象記錄以來氣候相關(guān)的研究進行了對比,以此來分析該地區(qū)生長釋放和抑制產(chǎn)生的原因。本研究發(fā)現(xiàn)黎城樹木在1860—1865以及1873—1877年間均發(fā)生了明顯的生長抑制。根據(jù)史料記載長治地區(qū)1862年發(fā)生了嚴重的干旱事件,同樣的在1875—1878年間的降雨量也非常小[30]。此時,根據(jù)樹輪重建的降水序列顯示山西寧武地區(qū)也處于持續(xù)的干旱時期[31]。3個點樹輪年表在1928—1931年形成極窄年輪及缺失輪,且在1925—1930期間都有生長抑制現(xiàn)象。20世紀20年代末期的旱災(zāi),是中國近代十大災(zāi)荒之一[32],而且多次被樹輪的生長異常檢測到[33,34]。根據(jù)李強等對山西寧武地區(qū)樹輪重建的降水序列顯示該地區(qū)在1901—1930年降水量平均值為399mm,降水量特別少導(dǎo)致持續(xù)的干旱時期[31]。長治地區(qū)20世紀70—90年代除春季降水處于增多趨勢外,夏、秋、冬季降水均處于不同程度的減少趨勢,處于長期干旱時期[35]。黎城樣本在70年代后生長緩慢,表現(xiàn)出較弱的生長抑制。而重建的降水序列也顯示1972年降水量只有324mm,遠遠低于多年平均值的413mm[31]。根據(jù)長治市氣象站的記錄,在1995—2000期間長治地區(qū)的溫度呈持續(xù)上升趨勢,2000年后溫度下降趨于平穩(wěn)。蔡秋芳等對山西呂梁山氣溫變化的研究中發(fā)現(xiàn),1994—2002年的溫度是過去近170年中溫度最高的時段[10]。平順和屯留的白皮松在這個時間段內(nèi)表現(xiàn)出生長抑制現(xiàn)象。持續(xù)異常的高溫會加速地表水分蒸發(fā)與樹木體內(nèi)水分的蒸騰,當溫度恢復(fù)正常水平后,樹木出現(xiàn)釋放事件。3個點生長抑制和釋放與極端氣候事件在時間上的一致性,表明該地區(qū)降水減少或高溫下水分蒸散所帶來的極端干旱是該區(qū)域森林生長抑制的主要原因。

    本研究區(qū)域內(nèi)樹齡超過150a的古樹有7棵,其中最大年齡為227a。樹木自身的生理功能也就是光合作用、水分傳導(dǎo)以及耐旱能力等隨著年齡而發(fā)生變化[36]。不同齡級的樹木的抗旱,抗高溫能力也是不同的。對于老齡樹,夏季的水分需要量更多,在水分條件一定的條件下,夏季較高的溫度,加強了蒸騰速率,造成樹木生長的干旱脅迫,因此老齡樹木對夏季氣溫的變化更敏感[37]。能夠經(jīng)歷多次各種干擾而存留下來的老齡樹具有較強的抗干擾能力。老齡樹復(fù)雜的年齡結(jié)構(gòu)常常能反映過去氣候的變化情況,也記錄了該樹種從生到死的總規(guī)律[38]。因此我們要加強對老齡樹的針對性保護以此來保存珍貴的歷史信息,維護區(qū)域森林生態(tài)服務(wù)功能。此外,森林生長釋放時期良好的環(huán)境適宜樹木的生長,營林建設(shè)更宜選擇在此時。

    4 結(jié)論與展望

    本文利用長治地區(qū)3個縣采樣點的樹輪數(shù)據(jù)計算了樹木生長的生長變化百分率,以此重建了過去150年間森林的生長抑制和釋放歷史,結(jié)果表明該地區(qū)森林在過去150年期間發(fā)生生長釋放和抑制的時間基本上是同步的,但是也存在由于樹種以及微環(huán)境的差異所產(chǎn)生的微小不同,主要有3次大的生長抑制事件和5次大的生長釋放事件。而通過樹木生長與氣候響應(yīng)分析結(jié)果以及其他相關(guān)歷史資料的佐證,表明該地區(qū)森林抑制現(xiàn)象主要是由于極端干旱事件引起的干擾導(dǎo)致的。主要的生長抑制和釋放時間與史料以及其他相關(guān)研究結(jié)果一致,表明本文的研究準確可靠,說明該地區(qū)利用樹輪來研究森林的干擾歷史是可行的。研究結(jié)果還為林業(yè)部門提供了合理的森林經(jīng)營建議,為今后進一步研究整個黃土高原的森林生長干擾歷史,以及群落的更新和保護提供了可供參考的方法。

    致謝:感謝中國科學(xué)院植物研究所樹木年輪實驗室提供的樹輪樣本、張齊兵研究員在數(shù)據(jù)分析和論文寫作上的悉心指導(dǎo)、邱紅巖老師對實驗室工作的熱情幫助。

    [1] 王樹芝, 趙秀海. 樹輪生態(tài)學(xué)研究進展. 世界林業(yè)研究, 2010, 23(2): 17- 21.

    [2] 封曉輝, 程瑞梅, 肖文發(fā), 王曉榮, 王瑞麗. 樹木年輪在干擾歷史重建中的應(yīng)用. 生態(tài)學(xué)報, 2011, 31(11): 3215- 3222.

    [3] 侯向陽, 韓進軒, 謝海生. 長白山紅松林干擾節(jié)律研究. 生態(tài)學(xué)報, 2000, 20(3): 409- 411.

    [4] 夏冰, 賀善安, 蘭濤, 鄧飛, 姚淦. 亞高山云冷杉混交林樹木生長釋放與干擾分析. 植物資源與環(huán)境, 1997, 6(1): 1- 8.

    [5] 劉錄三, 邵雪梅, 梁爾源, 王麗麗. 祁連山中部祁連圓柏生長與更新方式的樹輪記錄. 地理研究, 2006, 25(1): 53- 61.

    [6] 王曉春, 趙玉芳. 黑河勝山國家自然保護區(qū)紅松和紅皮云杉生長釋放判定及解釋. 生態(tài)學(xué)報, 2011, 31(5): 1230- 1239.

    [7] 姚玉璧, 王潤元, 楊金虎, 岳平, 陸登榮, 肖國舉, 王洋, 劉林春. 黃土高原陸地表層最大可能蒸散量的變化特征. 生態(tài)環(huán)境學(xué)報, 2011, 20(8- 9): 1189- 1195.

    [8] 李晉昌, 劉勇, 張彩霞. 山西省春季氣溫、降水及其極端事件的變化. 干旱區(qū)資源與環(huán)境, 2010, 24(10): 55- 60.

    [9] 張文濤, 江源, 王明昌, 張凌楠, 董滿宇, 郭媛媛. 蘆芽山陽坡不同海拔白杄徑向生長對氣候變暖的響應(yīng). 植物生態(tài)學(xué)報, 2013, 37(12): 1142- 1152.

    [10] 蔡秋芳, 劉禹, 包光, 雷鶯, 孫鉑. 樹輪記錄的呂梁山地區(qū)公元1836年以來5—7月平均氣溫變化. 科學(xué)通報, 2010, 55(20): 2033- 2039.

    [11] 易亮, 劉禹, 宋惠明, 李強, 蔡秋芳, 楊銀科, 孫軍艷. 山西蘆芽山地區(qū)樹木年輪記錄的1676 AD以來5—7月溫度變化. 冰川凍土, 2006, 28(3): 330- 336.

    [12] 李穎俊, 王尚義, ??〗? 方克艷, 李曉嵐, 栗燕, 布文麗, 李玉晗. 蘆芽山華北落葉松(Larixprincipis-rupprechtii)樹輪寬度年表對氣候因子的響應(yīng). 生態(tài)學(xué)報, 2016, 36(6): 1608- 1618.

    [13] 張桂萍, 栗建華. 長治市生態(tài)環(huán)境退化及成因分析. 山西大學(xué)學(xué)報: 自然科學(xué)版, 2010, 33(4): 620- 624.

    [14] 長治市地方志編纂委員會. 長治市地方志. 北京: 海潮出版社, 1995.

    [15] 劉鑫焱. 長治: 弘揚太行精神, 打造綠色家園. 人民日報, 2009-08- 20 [2015- 12- 25]. http://env.people.com.cn/GB/9893665.html

    [16] Cook E R, Kairiukstis L A. Methods of Dendrochronology: Applications in the Environmental Sciences. Dordrecht: Kluwer Academic Publishers, 1990: 51- 55.

    [17] Holmes R L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bulletin, 1983, 43: 69- 78.

    [18] Fritts H C. Tree Ring and Climate. London: Academic Press, 1976.

    [19] Schweingruber F H. Tree Rings and Environment Dendroecology. Berne: Paul Haupt, 1996.

    [20] Cook E R, Holmes R L. Users Manual for ARSTAN. Tucson: Laboratory of Tree-ring Research, University of Arizona, 1986.

    [21] Lorimer C G, Frelich L E. A methodology for estimating canopy disturbance frequency and intensity in dense temperate forests. Canadian Journal of Forest Research, 1989, 19(5): 651- 663.

    [22] Payette S, Filion L, Delwaide A. Disturbance regime of a cold temperate forest as deduced from tree-ring patterns: the Tantaré Ecological Reserve, Quebec. Canadian Journal of Forest Research, 1990, 20(8): 1228- 1241.

    [23] Nowacki G J, Abrams M D. Radial-growth averaging criteria for reconstructing disturbance histories from presettlement-origin oaks. Ecological Monographs, 1997, 67(2): 225- 249.

    [24] Biondi F, Waikul K. DENDROCLIM2002: a C++ program for statistical calibration of climate signals in tree-ring chronologies. Computers & Geosciences, 2004, 30(3): 303- 311.

    [25] 楊建偉, 梁宗鎖, 韓蕊蓮, 王培榛. 土壤干旱對油松生長及水分利用的影響. 西北農(nóng)林科技大學(xué)學(xué)報: 自然科學(xué)版, 2004, 32(4): 88- 92.

    [26] 侯迎, 王乃昂, 李鋼, 鄭芳. 利用樹輪資料重建1751- 2005年崆峒山地區(qū)夏季溫度變化. 氣候變化研究進展, 2007, 3(3): 172- 176.

    [27] 張文濤, 江源, 王明昌, 張凌楠, 董滿宇. 蘆芽山陽坡不同海拔華北落葉松徑向生長對氣候變化的響應(yīng). 生態(tài)學(xué)報, 2015, 35(19): 6481- 6488.

    [28] 蔡秋芳, 劉禹, 宋慧明, 孫軍艷. 樹輪記錄的陜西中-北部地區(qū)1826年以來4—9月溫度變化. 中國科學(xué) D輯: 地球科學(xué), 2008, 38(8): 971- 977.

    [29] 王婷, 李聰, 張弘, 任思遠, 李鹿鑫, 潘娜, 袁志良, 葉永忠. 寶天曼自然保護區(qū)不同針葉樹徑向生長對氣候的響應(yīng). 生態(tài)學(xué)報, 2016, 36(17): 1- 9.

    [30] 錢林清, 鄭炎謀, 胡慧敏. 山西氣候. 北京: 氣象出版社, 1991: 1- 458.

    [31] 李強, 劉禹, 蔡秋芳, 孫軍艷, 易亮, 宋惠明, 王雷. 山西寧武地區(qū)1686年以來年降水重建. 第四紀研究, 2006, 26(6): 999- 1006.

    [32] 陜西歷史自然災(zāi)害簡要紀實編委會. 陜西歷史自然災(zāi)害簡要紀實. 北京: 氣象出版社, 2002.

    [33] 梁爾源, 邵雪梅, 黃磊, 王麗麗. 中國中西部地區(qū)樹木年輪對20世紀20年代干旱災(zāi)害的指示. 自然科學(xué)進展, 2004, 14(4): 469- 474.

    [34] Liang E Y, Liu X H, Yuan Y J, Qin N S, Fang X Q, Huang L, Zhu H F, Wang L L, Shao X M. The 1920s drought recorded by tree rings and historical documents in the semi-arid and arid areas of northern China. Climatic Change, 2006, 79(3/4): 403- 432.

    [35] 趙雙巧, 李晶晶, 張卉, 張紅英, 楊梅. 長治市近54年降水時空變化特征分析. 科技與創(chuàng)新, 2004, (14): 142- 144.

    [36] Vieira J, Campelo F, Nabais C. Age-dependent responses of tree-ring growth and intra-annual density fluctuations ofPinuspinasterto Mediterranean climate. Trees, 2009, 23(2): 257- 265.

    [37] 宋來萍. 大興安嶺樟子松年輪氣候響應(yīng)及年齡影響探討[D]. 哈爾濱: 東北林業(yè)大學(xué), 2011.

    [38] 胡堅強, 夏有根, 梅艷, 王學(xué)勤. 古樹名木研究概述. 福建林業(yè)科技, 2004, 31(3): 151- 154.

    History of growth suppression and release events in forests in Changzhi Prefecture, Shanxi Province, China

    ZHANG Qi1, YAN Ming1, LIANG Hanxue2, 3,*

    1CollegeofLifeScience,ShanxiNormalUniversity,Linfen041004,China2StateKeyLaboratoryofVegetationandEnvironmentalChange,InstituteofBotany,ChineseAcademyofSciences,Beijing100093,China3UniversityofChineseAcademyofSciences,Beijing100049,China

    Forest growth is closely associated with climate change, and understanding the history of growth disturbances can help predict potential changes in forest growth caused by global warming. In this study, we utilized dendrochronological methods to study growth release and suppression events of trees in order to reconstruct the disturbance history of forests. A total of 27 samples were collected from onePinustabuliformisforest and twoPinusbungeanaforests in Changzhi Prefecture (southeast of Loess Plateau). After measuring and cross-dating the tree ring width, we found that the oldest trees in the Licheng, Pingshun, and Tunliu sampling sites were 227, 185, and 102 years old, respectively. Percentage of growth change (GC%) was calculated to identify the disturbance events. Growth release was defined as an averageGC% (five years) of more than 0.25, while growth suppression was defined as an averageGC% (five years) of less than-0.25. In this study, the results indicated that in the past 150 years, growth suppression occurred three times (1873—1877, 1925—1930, and 1994—1997), while growth release occurred five times (1867—1871, 1878—1884, 1930—1935, 1980—1985, and 1999—2004). Although the three forests had similarities in their disturbance histories, there were also differences in their spatial and temporal scales. These differences may be caused by the genetic characteristics of the different species and local environmental variations. To further understand the causes of the disturbances, a correlation analysis was performed for the first order difference of the three tree ring width chronologies and the meteorological records. The correlation coefficients showed that forest growth was negatively correlated with June temperature and positively correlated with April—May precipitation. This suggested that the historical growth suppression events might be related to extreme drought events, caused by precipitation reduction or excessive water evapotranspiration, whereas the growth release events might be related to an abundant supply of water. To validate these results, we compared our study with historical records and local paleoclimate studies. The drought events recorded in the historical materials, and reconstructed using other paleoclimate studies reasonably fit our results, which indicates that our study successfully rebuilt the historical growth history and that dendrochronological methods could be used in forest disturbance research in the east edge of the Loess Plateau. This study provided useful information for forestry management. For example, forest plantations should be carried out in moist conditions during the growing season. The old growth trees that survived multiple disturbance events have strong resilience and contain valuable information about historical environmental changes; identifying and protecting old growth trees is important for maintaining efficient and healthy forest ecosystem services.

    the Loess Plateau; Changzhi Prefecture; tree ring; growth suppression; growth release

    國家自然科學(xué)基金項目(31330015)

    2016- 01- 27; 網(wǎng)絡(luò)出版日期:2016- 12- 19

    10.5846/stxb201601270193

    *通訊作者Corresponding author.E-mail: dtxlhx@163.com

    張啟,閆明,梁寒雪.山西省長治市過去150年森林的生長抑制和釋放歷史.生態(tài)學(xué)報,2017,37(9):3115- 3123.

    Zhang Q, Yan M, Liang H X.History of growth suppression and release events in forests in Changzhi Prefecture, Shanxi Province, China.Acta Ecologica Sinica,2017,37(9):3115- 3123.

    猜你喜歡
    樹輪年輪樹木
    年輪
    樹輪研究方法與進展
    天山南北坡樹輪穩(wěn)定碳同位素對氣候的響應(yīng)差異
    樹木之最
    年輪
    辨認樹木
    樹木也睡覺
    心事
    滇池(2017年5期)2017-05-19 03:22:18
    樹木年輪對生態(tài)環(huán)境變化的響應(yīng)及應(yīng)用研究
    為什么樹有年輪
    欧美另类亚洲清纯唯美| 国产精品99久久久久久久久| 欧美又色又爽又黄视频| 久久精品综合一区二区三区| 99久久国产精品久久久| 亚洲中文日韩欧美视频| 每晚都被弄得嗷嗷叫到高潮| 国产三级黄色录像| 欧美日韩一级在线毛片| 国产伦精品一区二区三区四那| 精品国产超薄肉色丝袜足j| 亚洲美女视频黄频| 色播亚洲综合网| 亚洲成a人片在线一区二区| 国产 一区 欧美 日韩| 精品福利观看| 久久精品aⅴ一区二区三区四区| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 亚洲va日本ⅴa欧美va伊人久久| 女人高潮潮喷娇喘18禁视频| 久久久国产精品麻豆| 最近在线观看免费完整版| 日韩欧美国产在线观看| 午夜福利在线观看免费完整高清在 | 亚洲国产色片| 日日摸夜夜添夜夜添小说| 国产一级毛片七仙女欲春2| 精品国产美女av久久久久小说| 亚洲av第一区精品v没综合| 夜夜夜夜夜久久久久| 一级毛片高清免费大全| 国产三级在线视频| 岛国在线免费视频观看| 99热这里只有精品一区 | 九色成人免费人妻av| 国产亚洲精品久久久com| 制服人妻中文乱码| 999精品在线视频| 97超视频在线观看视频| 午夜精品久久久久久毛片777| 国产成人av激情在线播放| 亚洲成人免费电影在线观看| 国产成年人精品一区二区| 久久久久九九精品影院| 免费在线观看日本一区| 好男人电影高清在线观看| 久久久水蜜桃国产精品网| 真人一进一出gif抽搐免费| 亚洲在线自拍视频| 男人舔女人的私密视频| 日本免费a在线| 97超级碰碰碰精品色视频在线观看| 久久久久久久久中文| www.自偷自拍.com| 黄频高清免费视频| 成人18禁在线播放| 黄色 视频免费看| 久久精品影院6| 亚洲国产色片| 丁香六月欧美| 在线视频色国产色| 51午夜福利影视在线观看| 欧美中文日本在线观看视频| 欧美日韩精品网址| 99热只有精品国产| 麻豆国产av国片精品| 国产成人av激情在线播放| 国产一区二区三区视频了| 天堂√8在线中文| 一级作爱视频免费观看| 免费一级毛片在线播放高清视频| 久久国产乱子伦精品免费另类| 在线观看美女被高潮喷水网站 | 搡老妇女老女人老熟妇| 亚洲第一欧美日韩一区二区三区| 麻豆成人av在线观看| 国产91精品成人一区二区三区| 一本久久中文字幕| 国产精品野战在线观看| 搡老岳熟女国产| 国产成人精品久久二区二区91| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人免费电影在线观看| 亚洲精品美女久久av网站| 午夜精品久久久久久毛片777| 免费一级毛片在线播放高清视频| 丁香欧美五月| 男人的好看免费观看在线视频| 观看免费一级毛片| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 美女高潮的动态| 免费看光身美女| 美女午夜性视频免费| 日日干狠狠操夜夜爽| 成年女人毛片免费观看观看9| 久久久久国产精品人妻aⅴ院| 女同久久另类99精品国产91| 黄片小视频在线播放| 欧美黑人巨大hd| www.熟女人妻精品国产| 久久久久久人人人人人| 男女之事视频高清在线观看| 国产高清videossex| 免费av不卡在线播放| 久久精品国产清高在天天线| 99精品久久久久人妻精品| 看片在线看免费视频| 婷婷精品国产亚洲av在线| 香蕉丝袜av| 欧美成人免费av一区二区三区| 老司机深夜福利视频在线观看| 国产不卡一卡二| 色精品久久人妻99蜜桃| 成人精品一区二区免费| 最新中文字幕久久久久 | 搞女人的毛片| 亚洲av免费在线观看| 久9热在线精品视频| 天堂av国产一区二区熟女人妻| 精品人妻1区二区| 亚洲欧美一区二区三区黑人| xxx96com| 国产精品综合久久久久久久免费| 91在线精品国自产拍蜜月 | 国产精品1区2区在线观看.| 国产91精品成人一区二区三区| 麻豆国产av国片精品| 91在线观看av| 亚洲,欧美精品.| tocl精华| 国产精品乱码一区二三区的特点| 97超视频在线观看视频| 免费看日本二区| 欧美日韩精品网址| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 成人高潮视频无遮挡免费网站| 在线永久观看黄色视频| 99久久国产精品久久久| 免费av不卡在线播放| 亚洲av成人av| 亚洲av电影不卡..在线观看| 免费在线观看影片大全网站| 69av精品久久久久久| 国产精品一区二区免费欧美| 亚洲av成人不卡在线观看播放网| 男女视频在线观看网站免费| 九色国产91popny在线| 又粗又爽又猛毛片免费看| 亚洲精品中文字幕一二三四区| 校园春色视频在线观看| ponron亚洲| 黄频高清免费视频| 国产探花在线观看一区二区| 午夜a级毛片| 一边摸一边抽搐一进一小说| 国产精品 欧美亚洲| 亚洲av片天天在线观看| 免费无遮挡裸体视频| 十八禁人妻一区二区| 真人一进一出gif抽搐免费| 级片在线观看| 国产乱人伦免费视频| 美女免费视频网站| 亚洲狠狠婷婷综合久久图片| 可以在线观看的亚洲视频| 99久久精品一区二区三区| 大型黄色视频在线免费观看| 久久这里只有精品中国| 国产三级黄色录像| 一本久久中文字幕| 蜜桃久久精品国产亚洲av| 色综合婷婷激情| 亚洲电影在线观看av| 九九久久精品国产亚洲av麻豆 | 日韩有码中文字幕| 久久久久久久久免费视频了| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 精品乱码久久久久久99久播| 国产成人精品无人区| 国产精品亚洲av一区麻豆| 人人妻人人看人人澡| 久久久久久九九精品二区国产| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩卡通动漫| 国产探花在线观看一区二区| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 岛国在线观看网站| 久久久精品欧美日韩精品| 国产成人啪精品午夜网站| 黄色 视频免费看| 搡老妇女老女人老熟妇| 免费一级毛片在线播放高清视频| 国产欧美日韩一区二区三| 久久这里只有精品19| 老司机在亚洲福利影院| 综合色av麻豆| 亚洲人成网站高清观看| 久久精品aⅴ一区二区三区四区| 操出白浆在线播放| 黑人巨大精品欧美一区二区mp4| 亚洲av日韩精品久久久久久密| 精品99又大又爽又粗少妇毛片 | 校园春色视频在线观看| 日本五十路高清| 嫁个100分男人电影在线观看| 黑人操中国人逼视频| 99久久精品热视频| 又爽又黄无遮挡网站| 国产一区二区三区视频了| 性色av乱码一区二区三区2| 亚洲在线自拍视频| 精品国产美女av久久久久小说| 亚洲人成伊人成综合网2020| 在线观看免费视频日本深夜| 亚洲aⅴ乱码一区二区在线播放| 美女被艹到高潮喷水动态| 热99在线观看视频| 中出人妻视频一区二区| 日本五十路高清| 日本一二三区视频观看| 国产乱人伦免费视频| 久久久久久大精品| 怎么达到女性高潮| 国产精品,欧美在线| 欧美黄色淫秽网站| 亚洲专区中文字幕在线| 18禁裸乳无遮挡免费网站照片| 看片在线看免费视频| 免费看a级黄色片| 久久中文字幕一级| 在线观看美女被高潮喷水网站 | 久久久久久人人人人人| 色在线成人网| 九色成人免费人妻av| 在线观看免费午夜福利视频| netflix在线观看网站| 美女免费视频网站| 岛国在线免费视频观看| 精品久久久久久久久久久久久| 国产亚洲av高清不卡| av欧美777| 美女cb高潮喷水在线观看 | 少妇丰满av| 午夜免费激情av| 人妻丰满熟妇av一区二区三区| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| 成人av在线播放网站| 免费在线观看视频国产中文字幕亚洲| 丰满的人妻完整版| cao死你这个sao货| 国产视频内射| 国产免费男女视频| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩卡通动漫| 看免费av毛片| 国产男靠女视频免费网站| 亚洲国产欧美人成| 99久久精品国产亚洲精品| 成人18禁在线播放| 一二三四在线观看免费中文在| 国内揄拍国产精品人妻在线| 亚洲成av人片免费观看| 国产伦精品一区二区三区四那| 欧美日本亚洲视频在线播放| 免费av毛片视频| 毛片女人毛片| 韩国av一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| 两个人视频免费观看高清| 欧美中文日本在线观看视频| 91在线观看av| 一a级毛片在线观看| 精品久久久久久久末码| 蜜桃久久精品国产亚洲av| 久久精品91蜜桃| 亚洲无线在线观看| 中亚洲国语对白在线视频| 国产成人av激情在线播放| 国产一级毛片七仙女欲春2| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 成年女人看的毛片在线观看| 日本精品一区二区三区蜜桃| 欧美丝袜亚洲另类 | 亚洲午夜精品一区,二区,三区| 宅男免费午夜| 午夜精品久久久久久毛片777| 91在线观看av| 国产精品香港三级国产av潘金莲| 亚洲专区国产一区二区| 久久精品aⅴ一区二区三区四区| 久9热在线精品视频| 国产高清三级在线| 日韩人妻高清精品专区| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 欧美不卡视频在线免费观看| or卡值多少钱| 国产精品久久电影中文字幕| 在线观看免费视频日本深夜| 日本免费一区二区三区高清不卡| 在线观看美女被高潮喷水网站 | 国产三级中文精品| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 欧美日韩黄片免| 久久久久久久午夜电影| 99热6这里只有精品| 在线国产一区二区在线| 国产私拍福利视频在线观看| 国产成人一区二区三区免费视频网站| 国产精品乱码一区二三区的特点| 草草在线视频免费看| 成人av在线播放网站| 韩国av一区二区三区四区| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 手机成人av网站| 亚洲18禁久久av| e午夜精品久久久久久久| 精品福利观看| 国产欧美日韩精品亚洲av| 少妇的丰满在线观看| 黄色 视频免费看| 免费看十八禁软件| 久久久久免费精品人妻一区二区| 国产av不卡久久| 免费在线观看亚洲国产| 又黄又爽又免费观看的视频| 精品国产三级普通话版| 中亚洲国语对白在线视频| 性欧美人与动物交配| 免费在线观看视频国产中文字幕亚洲| 国产一区在线观看成人免费| 国产成人精品久久二区二区91| 老司机福利观看| 午夜两性在线视频| 深夜精品福利| 99久国产av精品| 亚洲va日本ⅴa欧美va伊人久久| 伊人久久大香线蕉亚洲五| 19禁男女啪啪无遮挡网站| 国产亚洲精品久久久com| 看黄色毛片网站| 97人妻精品一区二区三区麻豆| 欧美成人性av电影在线观看| 亚洲国产精品久久男人天堂| 人妻夜夜爽99麻豆av| 国内毛片毛片毛片毛片毛片| 99久久成人亚洲精品观看| 91字幕亚洲| 少妇裸体淫交视频免费看高清| 少妇的丰满在线观看| 亚洲精品一卡2卡三卡4卡5卡| 天堂√8在线中文| 99久久国产精品久久久| 亚洲成av人片在线播放无| 国产精品av视频在线免费观看| 久99久视频精品免费| 国产黄片美女视频| 97超视频在线观看视频| 91在线观看av| 熟女电影av网| 精品乱码久久久久久99久播| 欧美日韩一级在线毛片| 国产免费男女视频| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器 | 亚洲精品在线美女| 免费观看人在逋| 亚洲欧美日韩高清在线视频| 精品人妻1区二区| 欧美日韩乱码在线| 欧美大码av| 国产av在哪里看| 国产成人精品无人区| 男女那种视频在线观看| 国产野战对白在线观看| 国内精品久久久久久久电影| 亚洲电影在线观看av| 又大又爽又粗| 老司机在亚洲福利影院| 国产一区二区在线观看日韩 | 天堂动漫精品| 日本 欧美在线| svipshipincom国产片| 午夜福利在线观看免费完整高清在 | 一本久久中文字幕| 人妻丰满熟妇av一区二区三区| 大型黄色视频在线免费观看| 偷拍熟女少妇极品色| 国产高清三级在线| 欧美在线黄色| 好男人在线观看高清免费视频| 成人永久免费在线观看视频| 一进一出好大好爽视频| 毛片女人毛片| 国产精品99久久久久久久久| 波多野结衣高清作品| www.精华液| 我的老师免费观看完整版| 国内毛片毛片毛片毛片毛片| 91老司机精品| 国产精品久久久av美女十八| 性色av乱码一区二区三区2| 久久久久国产一级毛片高清牌| 国产毛片a区久久久久| 国产精品免费一区二区三区在线| 免费观看的影片在线观看| 亚洲人成网站高清观看| 99久久久亚洲精品蜜臀av| 久久中文字幕一级| 国产高清videossex| 成人特级黄色片久久久久久久| 成人特级av手机在线观看| 成人亚洲精品av一区二区| 999精品在线视频| 男插女下体视频免费在线播放| 搞女人的毛片| avwww免费| 久久精品91无色码中文字幕| 热99re8久久精品国产| 精品久久久久久久久久免费视频| 白带黄色成豆腐渣| 九色成人免费人妻av| 少妇的丰满在线观看| h日本视频在线播放| 国产av不卡久久| 午夜福利成人在线免费观看| 99久久久亚洲精品蜜臀av| 亚洲五月天丁香| 欧美高清成人免费视频www| 男女床上黄色一级片免费看| 成人18禁在线播放| 精品电影一区二区在线| 国产精品九九99| 亚洲欧洲精品一区二区精品久久久| 国产爱豆传媒在线观看| 国产精品久久视频播放| 亚洲国产欧美网| 长腿黑丝高跟| 日韩欧美国产在线观看| av国产免费在线观看| 很黄的视频免费| 老汉色∧v一级毛片| 高清在线国产一区| 两个人视频免费观看高清| 91在线观看av| 久久香蕉国产精品| 九色国产91popny在线| 亚洲激情在线av| 国产又黄又爽又无遮挡在线| 国产私拍福利视频在线观看| 黑人欧美特级aaaaaa片| 国产精品一及| 亚洲男人的天堂狠狠| 人人妻人人澡欧美一区二区| 精品无人区乱码1区二区| 国产精品影院久久| 99久久精品国产亚洲精品| 亚洲乱码一区二区免费版| 操出白浆在线播放| 女人高潮潮喷娇喘18禁视频| 成人鲁丝片一二三区免费| 91av网一区二区| 成人三级做爰电影| 一a级毛片在线观看| www.精华液| 丰满的人妻完整版| 一本精品99久久精品77| 老熟妇仑乱视频hdxx| 国产精品美女特级片免费视频播放器 | 精品电影一区二区在线| 欧美日韩福利视频一区二区| 男人舔奶头视频| 两人在一起打扑克的视频| 成人三级黄色视频| 99久久综合精品五月天人人| 亚洲美女视频黄频| 亚洲天堂国产精品一区在线| 黄片小视频在线播放| 精品久久久久久久久久免费视频| 亚洲18禁久久av| 脱女人内裤的视频| 琪琪午夜伦伦电影理论片6080| 亚洲av成人一区二区三| 久久久久久国产a免费观看| 日本在线视频免费播放| 午夜福利在线在线| 成人鲁丝片一二三区免费| 亚洲国产精品999在线| 久久久国产成人精品二区| 欧美一区二区精品小视频在线| 在线观看午夜福利视频| 亚洲欧美日韩东京热| 综合色av麻豆| 免费观看精品视频网站| 亚洲国产欧洲综合997久久,| 99精品在免费线老司机午夜| 色综合站精品国产| 国产精品亚洲美女久久久| 舔av片在线| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久精品吃奶| 久久精品夜夜夜夜夜久久蜜豆| aaaaa片日本免费| 亚洲美女黄片视频| 91麻豆av在线| 久久精品国产亚洲av香蕉五月| 精品久久久久久成人av| 热99在线观看视频| 观看免费一级毛片| 99在线人妻在线中文字幕| 中文字幕高清在线视频| 听说在线观看完整版免费高清| 亚洲最大成人中文| 男女视频在线观看网站免费| 久久香蕉精品热| 1024手机看黄色片| 嫩草影视91久久| 欧美黄色淫秽网站| 99re在线观看精品视频| 少妇的丰满在线观看| 久久亚洲精品不卡| 国产av一区在线观看免费| 国产91精品成人一区二区三区| 国内精品久久久久精免费| 欧美性猛交黑人性爽| 99国产精品一区二区三区| 成人午夜高清在线视频| 免费av不卡在线播放| 香蕉国产在线看| 亚洲专区中文字幕在线| 国产单亲对白刺激| 国产主播在线观看一区二区| 久久久久久久久免费视频了| 深夜精品福利| 国产精品电影一区二区三区| 久久久久久大精品| 国产精品99久久久久久久久| 搡老妇女老女人老熟妇| h日本视频在线播放| 日本一二三区视频观看| 男人舔女人的私密视频| 桃色一区二区三区在线观看| www日本黄色视频网| www.精华液| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 在线观看美女被高潮喷水网站 | 亚洲激情在线av| 精品国产乱子伦一区二区三区| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| 精品久久久久久久末码| 欧美日韩黄片免| 久久这里只有精品19| 美女大奶头视频| 国产精品久久电影中文字幕| 黄色女人牲交| 欧美一级a爱片免费观看看| 天堂√8在线中文| 99久久精品一区二区三区| 天天添夜夜摸| 亚洲午夜精品一区,二区,三区| 18禁裸乳无遮挡免费网站照片| 999久久久国产精品视频| 最新在线观看一区二区三区| www日本黄色视频网| 中出人妻视频一区二区| 草草在线视频免费看| 国内精品久久久久精免费| 国产精品自产拍在线观看55亚洲| 超碰成人久久| www.www免费av| 午夜视频精品福利| av片东京热男人的天堂| 99久久精品热视频| 99精品在免费线老司机午夜| 偷拍熟女少妇极品色| 国产精品98久久久久久宅男小说| 欧美zozozo另类| 亚洲午夜精品一区,二区,三区| 午夜福利视频1000在线观看| 不卡一级毛片| 非洲黑人性xxxx精品又粗又长| 此物有八面人人有两片| 欧美黄色淫秽网站| 最新中文字幕久久久久 | 国产午夜精品论理片| 日韩欧美免费精品| 色老头精品视频在线观看| 天天躁日日操中文字幕| 一区二区三区高清视频在线| 在线免费观看不下载黄p国产 | 精品一区二区三区视频在线 | 国产黄a三级三级三级人| 国产伦一二天堂av在线观看| 男插女下体视频免费在线播放| 国产男靠女视频免费网站| 中文字幕最新亚洲高清| 噜噜噜噜噜久久久久久91| 变态另类成人亚洲欧美熟女| 无人区码免费观看不卡| 久久精品国产亚洲av香蕉五月| 看片在线看免费视频| 国产精品 国内视频| 欧美+亚洲+日韩+国产| 欧美国产日韩亚洲一区|