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

    黃土高原西部3個降水量梯度近成熟油松人工林碳庫特征

    2017-06-26 11:28:26孫美美關(guān)晉宏吳春榮岳軍偉李國慶
    生態(tài)學報 2017年8期
    關(guān)鍵詞:碳庫油松人工林

    孫美美,關(guān)晉宏,,吳春榮,岳軍偉,李國慶,,杜 盛,,*

    1 西北農(nóng)林科技大學黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 楊凌 712100 2 中國科學院水利部水土保持研究所, 楊凌 712100 3 甘肅省治沙研究所, 蘭州 730070

    ?

    黃土高原西部3個降水量梯度近成熟油松人工林碳庫特征

    孫美美1,關(guān)晉宏1,2,吳春榮3,岳軍偉2,李國慶1,2,杜 盛1,2,*

    1 西北農(nóng)林科技大學黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室, 楊凌 712100 2 中國科學院水利部水土保持研究所, 楊凌 712100 3 甘肅省治沙研究所, 蘭州 730070

    森林在陸地生態(tài)系統(tǒng)吸收碳素方面起著主要作用,了解其固碳特征對研究地區(qū)之間的碳循環(huán)至關(guān)重要。油松人工林是黃土高原地區(qū)一種典型的退耕還林樹種,研究其固碳特征有利于綜合分析評價油松人工林的生態(tài)效益。為了研究黃土高原西部地區(qū)油松人工林碳儲量及碳密度空間分布特征因降水量不同引起的差異,以黃土高原西部地區(qū)3個典型栽培區(qū)域的近成熟油松人工林為對象,研究了群落內(nèi)各組成部分的生物量和碳庫特征。喬木層生物量的估算采用以胸徑和樹高為基礎(chǔ)變量的生物量方程,灌木、草本、凋落物采用樣方收獲法,土壤碳庫依據(jù)土壤剖面(1 m)和土鉆取樣相結(jié)合的方法測算。結(jié)果表明:在蘭州官蘑灘地區(qū)(372 mm)、太子山(519 mm)和小隴山(632 mm)3個不同降水量區(qū)域,油松人工林生物量碳密度分別為(49.08±2.86)t/hm2、(73.90±9.36)t/hm2和(82.55±7.36)t/hm2。小隴山地區(qū)的生態(tài)系統(tǒng)總碳密度和生物量碳密度與蘭州地區(qū)存在顯著性差異。在3個不同降水量區(qū)域,土壤有機碳密度僅在表層土壤(0—10 cm)差異達到顯著水平(P<0.05),而土壤總碳密度間差異未達到顯著水平(P>0.05)。黃土高原半干旱區(qū)近成熟油松人工林的生物量碳密度與年均降水量間呈現(xiàn)出顯著正相關(guān)關(guān)系。在半干旱地區(qū),降水量可能成為影響油松人工林生產(chǎn)力和碳固存的關(guān)鍵因素。

    油松人工林;生物量碳密度;降水量;黃土高原西部

    森林作為陸地生態(tài)系統(tǒng)的重要組成部分,儲存著陸地生態(tài)系統(tǒng)2/3以上的有機碳[1],并在減緩大氣CO2濃度上升、調(diào)節(jié)碳平衡等方面發(fā)揮著重要作用[2]。隨著退耕還林工程和三北防護林工程等重大生態(tài)項目的實施,我國人工林在森林總面積中的比重大幅增加[3],人工林在CO2吸收和固定方面的作用也引起越來越多的關(guān)注[4-5]。油松(Pinustabulaeformis)廣泛分布于中國北方地區(qū),耐旱耐貧瘠,適應性強,是黃土高原人工造林的主要樹種之一,具有良好的保持水土和生態(tài)防護等功能[6]。油松人工林是黃土高原一種典型植被,分析研究不同降水量區(qū)域其近成熟林的碳庫特征,有利于綜合評價該種人工林類型的碳匯效益。

    關(guān)于油松人工林生態(tài)系統(tǒng)碳密度的研究已有一些報道,例如不同林齡階段碳儲量比較和特定地區(qū)固碳潛力的估算[7],與其他樹種間碳匯情況的比較[8],碳密度與立地因子之間的關(guān)系[9],以及生態(tài)系統(tǒng)中碳氮儲量和分配格局[10]等方面。程小琴等[11]研究報道了山西近成熟油松人工林的生物量與碳庫空間分布特征,估算了其凈生產(chǎn)力與年均凈固碳量。目前研究較多的是不同林齡之間的比較,盡管涉及近成熟林,但特意針對近成熟林齡段油松人工林的研究卻甚少。對近成熟林的研究有助于了解該種森林類型在特定地區(qū)的最大生產(chǎn)力和固碳潛力。迄今為止,關(guān)于較大空間尺度上油松人工林的碳庫特征差異及其對降水量響應的研究仍未見報道。本研究選取黃土高原西部油松人工林的3個重要栽培區(qū)域,分析油松人工林生態(tài)系統(tǒng)各層次間碳密度特征,揭示降水量對油松人工林碳固存的影響,為估算不同地區(qū)油松人工林碳匯潛力和科學指導人工林經(jīng)營提供理論依據(jù)。

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

    1.1 研究區(qū)概況

    甘肅省蘭州市榆中縣官蘑灘,地處內(nèi)陸,大陸性季風氣候明顯,特點是降水量少,日照多,光能潛力大,氣候干燥,晝夜溫差大,年日照時數(shù)約2600 h,年平均降水量400 mm左右,并集中分布在6—9月,年平均氣溫7℃左右。調(diào)查區(qū)海拔為2360 m,森林土壤為山地灰褐土。樣地內(nèi)喬木層油松林樹種單一,林下幾乎沒有灌木,草本以長芒草(Stipabungeana)和細裂葉蓮蒿(Artemisiagmelinii)為主。

    甘肅省太子山林區(qū)(臨夏自治州康樂縣),氣候寒冷,生長期短,年均降水量在500—600 mm左右,土壤以山地褐土、棕壤為主[12]。調(diào)查區(qū)位于藥水保護站,海拔在2450 m左右,喬木層油松人工林樹種單一為絕對優(yōu)勢種,林下灌木主要有灰栒子(Cotoneasteracutifolius)和虎榛子(Ostryopsisdavidiana)等,草本主要以苔草(Carextristachya)、早熟禾(Poaannua)和三褶脈紫菀(Asterageratoides)為主。

    甘肅省小隴山林區(qū),地處秦嶺西段,是暖溫帶向亞熱帶的過渡地帶,年平均氣溫7—12℃,年均降水量600—900 mm,雨季集中在7、8、9月份,屬濕潤和半濕潤類型區(qū),土壤多屬壤土、森林褐土[13]。調(diào)查區(qū)位于隴南市兩當縣境內(nèi),海拔在1400 m左右,樣地內(nèi)油松為絕對優(yōu)勢種,林下灌木主要有胡枝子(Lespedezabicolor)、虎榛子(Ostryopsisdavidiana)和繡線菊(Spiraeasalicifolia)等,草本以苔草(Carextristachya)和唐松草(Thalictrumaquilegifolium)為主。

    1.2 研究方法

    1.2.1 調(diào)查樣地設(shè)置

    本研究遵循“代表性、均一性、連續(xù)性”原則設(shè)置調(diào)查標準地,原則上圍取標準地面積為1000 m2(50 m×20 m),但由于實際地形的限制,實際樣地中有部分設(shè)置為600 m2(30 m×20 m)。在甘肅省內(nèi)3個油松人工林重要栽培區(qū)(太子山、小隴山、蘭州)選取近成熟林(林齡40—60 a)的9塊樣地進行碳儲量和碳密度的研究。各標準地基本信息見表1。

    表1 不同降水量區(qū)域油松人工林標準地基本信息

    年降水量和年均氣溫的數(shù)據(jù)來自中國氣象數(shù)據(jù)網(wǎng)30 a(1981—2010)平均值[14]

    1.2.2 野外調(diào)查及樣品采集

    喬木層調(diào)查 對每個標準地內(nèi)胸徑大于2 cm的喬木進行每木檢尺,記錄樹高和胸徑,統(tǒng)計株數(shù)。按大、中、小徑級分別選擇3—5株樣木,分葉、枝、干、皮、根5個器官取樣,將相同器官的樣品混合,再取各器官不少于200 g重量的混合樣帶回實驗室待處理分析。

    灌木、草本和凋落物的調(diào)查 灌木調(diào)查是在每個喬木標準地內(nèi)沿對角線設(shè)置3個灌木樣方(2 m×2 m),調(diào)查樣方內(nèi)灌木種類,全部收獲后分葉、枝、根稱鮮重,將3個樣方內(nèi)枝、葉、根分別混合均勻后取混合樣,每個樣品重量不少于300 g;在每個灌木樣方內(nèi)設(shè)置一個1 m×1 m草本樣方,分地上與地下部分收獲后稱鮮重并取樣;凋落物的調(diào)查是將草本樣方內(nèi)的凋落物全部收獲后稱重并取重量不少于300 g的樣品。將灌木、草本和凋落物的樣品稱鮮重后帶回實驗室分析測定。

    土壤 在每個喬木標準地中選擇一個未受人為干擾、植被結(jié)構(gòu)和土壤具有代表性的地段,挖取土壤剖面深至100 cm,不足100 cm的挖至基巖為止,沿剖面按0—10,10—20,20—30,30—50,50—100 cm分層,用環(huán)刀取各個土層的原狀土,用于測定土壤容重。同時按照上述土壤分層標準,使用內(nèi)徑為4 cm的土鉆,分3個取樣點鉆取土鉆土,將等層樣品混合均勻后取200 g帶回實驗室分析測定。

    1.2.3 室內(nèi)樣品處理與碳、氮含量的分析測定

    將野外采集的喬木器官、灌木器官、草本及凋落物樣品置于85℃烘箱烘至恒重求得含水量。用于植物有機碳含量測定的烘干樣品需磨碎過0.25 mm篩。每個樣地的待測樣品均包括以下10部分:喬木的葉、枝、干、根,灌木的葉、枝、根,草本的地上、地下部分以及凋落物。待測土壤樣品在室內(nèi)陰干,壓碎大塊土,揀出植物根系等雜物,粉碎并過2 mm篩,將大于2 mm石礫集中稱重。再用四分法取部分樣品研磨直至全部土樣通過0.25 mm篩。植物樣品有機碳和土壤有機碳含量測定均采用重鉻酸鉀-硫酸氧化外加熱法,氮含量的測定采用凱氏定氮法[15]。

    1.2.4 生物量的測算

    油松生物量估算采用程堂仁[16]在這一地區(qū)建立的生物量方程(表2),構(gòu)建該方程的區(qū)域與本研究基本屬于同一地理區(qū),油松樹形無明顯變化。由胸徑和樹高求得單株油松各器官生物量,進而求得整株油松生物量和全樣地喬木生物量。

    表2 油松人工林生物量估算方程[16]

    B:生物量 Biomass;D:胸徑DBH (Diameter at Breast Height);H:樹高 Height

    灌木、草本及凋落物生物量的估算:將所采樣品帶回實驗室后置于85℃烘箱中,烘干至恒重,稱重計算含水率,進而推算出樣地內(nèi)單位面積(hm2)灌木各器官、草本的生物量及凋落物量。

    1.2.5 碳密度的估算

    各組分生物量碳密度采用各組分的生物量與其碳含量的乘積進行計算。土壤層(0- 100cm)碳密度為各土壤層碳密度之和,單獨某一土層的有機碳密度(Si,g/cm2)計算公式為:

    Si=Ci×Di×Ei×(1-Gi)

    式中,Ci為土壤碳含量(g/kg),Di為土壤容重(g/cm3),Ei為該層土壤的厚度(cm),Gi為直徑>2mm的石礫所占的體積百分比,i代表某一土層。本研究中全部土壤剖面無直徑>2mm的石礫,所以式中Gi取0。

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

    野外調(diào)查所得數(shù)據(jù)經(jīng)Excel 2010簡單整理后,用SPSS 20.0進行數(shù)據(jù)分析,選擇單因素方差分析(one-way ANOVA)比較參數(shù)之間的差異,用簡單相關(guān)分析和偏相關(guān)分析對林齡和年均降水量對生物量碳密度的關(guān)系進行分析,采用Sigmaplot 12.5軟件繪圖。

    2 結(jié)果與分析

    2.1 不同降水量地區(qū)油松人工林生物量、生物量碳密度差異

    蘭州(372 mm)和小隴山(632 mm)兩個降水量區(qū)域的油松人工林生態(tài)系統(tǒng)總生物量差異達到顯著性水平,而太子山(519 mm)與其他兩個地區(qū)的總生物量差異均不顯著(表3)。各調(diào)查區(qū)的生態(tài)系統(tǒng)生物量組分中喬木占比例最大,其中蘭州地區(qū)的喬木生物量與太子山、小隴山地區(qū)之間差異顯著。各降水量區(qū)域的凋落物生物量之間差異顯著性與總生物量之間趨勢一致,而草本生物量之間差異均未達到顯著水平(表3)。

    表3 不同降水量地區(qū)植被層生物量及其組成

    不同小寫字母表示同一組分在不同降水量區(qū)域間的差異性顯著(P<0.05)

    蘭州(372 mm)與小隴山(632 mm)地區(qū)的生物量碳密度之間呈顯著差異,但蘭州與太子山(519 mm)之間和太子山與小隴山之間差異均未達到顯著水平(圖1)。在蘭州地區(qū)喬木碳密度的貢獻率最大達到98.91%,主要是因為該地區(qū)的樣地內(nèi)幾乎沒有林下灌木。3個不同降水量區(qū)域,喬木碳密度的貢獻率總體上在95.21%—98.91%之間,灌木、草本以及凋落物碳密度貢獻率很小(0.43%—2.48%)。

    圖1 不同降水量區(qū)域生物量碳密度及各組分所占比例Fig.1 Biomass carbon density and component proportions in the three sites across a precipitation gradient不同字母表示同一生態(tài)系統(tǒng)組分不同降水量區(qū)域存在顯著性差異(P<0.05)

    2.2 不同降水量地區(qū)油松人工林土壤碳密度差異

    3個不同降水量地區(qū),樣地土壤表層(0—50 cm)各土層之間的碳密度存在較大差異,1 m深的整個土層內(nèi)總碳密度水平比較接近(圖2),差異沒有達到顯著性水平。降水量最多的小隴山林區(qū)(632 mm)表層土壤(0—10 cm)碳密度高于其他兩個地區(qū),但其50—100 cm土層的土壤碳密度與其他兩個地區(qū)差異不大,可能與該區(qū)域林地較高的凋落物生產(chǎn)量和因水分較充足的腐殖質(zhì)分解速率較快有關(guān)。

    2.3 不同降水量地區(qū)生態(tài)系統(tǒng)碳庫差異

    生態(tài)系統(tǒng)總碳庫中土壤層是最大的碳庫,土壤碳密度在184.90—207.07 t/hm2之間,是生物量碳密度的2—3倍(表4),且3個不同降水量地區(qū)之間的差異沒有達到顯著性水平。不同降水量地區(qū)僅在生物量碳密度之間存在一定顯著性差異,而土壤層碳庫和總碳庫之間差異不顯著。

    2.4 不同降水量地區(qū)近成熟油松人工林生物量碳密度與降水量之間的關(guān)系

    圖2 不同降水量區(qū)域各層土壤碳密度 Fig.2 Soil carbon density in the three sites across a precipitation gradient

    3個調(diào)查區(qū)中小隴山的年均降水量最大,分別是太子山和蘭州地區(qū)降水量的1.38和1.57倍,同時該區(qū)生物量碳密度是太子山和蘭州地區(qū)的1.12和1.68倍。簡單相關(guān)分析結(jié)果表明:近成熟油松人工林生物量碳密度與林齡之間不存在顯著相關(guān)性,而與年均降水量之間存在極顯著正相關(guān)關(guān)系,Pearson相關(guān)系數(shù)為0.820。由表5的偏相關(guān)分析結(jié)果可知,林齡與年均降水量影響近成熟油松人工林的生物量碳密的偏相關(guān)系數(shù)分別是-0.264和0.76,雙側(cè)檢驗結(jié)果顯示年均降水量對生物量碳密度的影響顯著,而林齡不顯著。3個調(diào)查區(qū)的近成熟油松人工林生物量碳密度隨降水量的增加而呈現(xiàn)出增加的趨勢,且降水量較少的蘭州地區(qū)的生物量碳密度與降水量較多的小隴山地區(qū)之間存在顯著性差異,表明降水量是決定區(qū)域間油松近成熟林生物量累積和碳固存差異的重要因素之一。

    表4 不同降水量地區(qū)生態(tài)系統(tǒng)碳庫變異

    不同小寫字母表示同一組分在不同降水量區(qū)域的差異性顯著(P<0.05)

    表5 近成熟油松人工林生物量碳密度與年均降水量和林齡的Pearson相關(guān)系數(shù)和偏相關(guān)系數(shù)

    *表示P<0.05, ** 表示P<0.01

    3 討論

    3.1 油松人工林生態(tài)系統(tǒng)碳庫特征

    本研究表明太子山和小隴山調(diào)查區(qū)的生物量碳密度分別是73.90 t/hm2和82.55 t/hm2,而蘭州地區(qū)的生物量碳密度則只有49.08 t/hm2,其中喬木層碳密度分別是71.01,78.60,48.54 t/hm2。這均比李??黐17]等人估算的中國油松林平均碳密度28.98 t/hm2分別高2.5,2.7倍和1.7倍,且只有蘭州地區(qū)的喬木層碳密度略低于我國溫帶針葉林和森林植被平均碳密度(分別為54.60 t/hm2和65.28 t/hm2)[18]。本研究調(diào)查的油松均是40—60 a的近成熟人工林,基本反映了相應區(qū)域該類人工林的碳庫特征??傮w而言,地處半干旱區(qū)的蘭州調(diào)查區(qū)油松人工林的碳固存現(xiàn)狀較低,而且,林下灌草植物也較少,應該已經(jīng)受到了水資源短缺的影響。

    土壤碳作為陸地生態(tài)系統(tǒng)最大的碳庫,在地球碳循環(huán)中起著至關(guān)重要的作用,稍有變動就會對全球碳儲量產(chǎn)生較大影響[19]。以往對中國土壤有機碳密度的研究結(jié)果差異較大,李克讓等[20]對中國常綠針葉林土壤有機碳密度的估算結(jié)果為179.80 t/hm2,孟蕾等[21]估算的子午嶺油松人工林土壤有機碳密度為93.83 t/hm2,本研究結(jié)果顯示油松人工林的土壤碳密度在200 t/hm2左右處于偏高水平。這可能是由于土壤類型、植被結(jié)構(gòu)、林齡的不同以及氣候和立地條件的差異造成了相關(guān)研究結(jié)果的不一致,說明今后針對不同地區(qū)、林分土壤有機碳庫的研究有著重要的意義。本研究表明3個降水量區(qū)域,土壤碳庫均是森林生態(tài)系統(tǒng)中最大的碳庫,占森林生態(tài)系統(tǒng)總碳庫的63.66%—78.57%。蘭州地區(qū)表層(0—10 cm)的碳密度較小,可能是因為蘭州地區(qū)的林下植被較少,凋落物少且降水量少,從而影響了土壤表層有機碳積累。隨著土層的加深土壤碳密度基本穩(wěn)定,這與生態(tài)系統(tǒng)中土壤有機碳的蓄積規(guī)律相吻合。

    本研究中3個降水量地區(qū)油松人工林生態(tài)系統(tǒng)碳庫之間差異不顯著,生物量碳密度和土壤層碳密度的不同差異使得總碳庫之間的差異微小。其中蘭州地區(qū)油松人工林生態(tài)系統(tǒng)碳庫最小251.70 t/hm2,太子山地區(qū)最大280.97 t/hm2,均比王寧等[22]研究報道的山西油松林生態(tài)系統(tǒng)平均碳密度127.53 t/hm2高,并且比馬欽彥等[23]研究的全國油松林平均碳密度138.00 t/hm2高。這可能是因為本研究中調(diào)查的油松人工林林齡均處于40—60 a之間,而山西省油松林90%屬于中幼齡林且森林質(zhì)量不高。本研究的結(jié)果與周玉榮等[17]研究的全國森林生態(tài)系統(tǒng)平均碳密度258.83 t/hm2基本一致。

    3個降水量區(qū)域的油松人工林生態(tài)系統(tǒng)碳密度均表現(xiàn)為土壤層碳密度>生物量碳密度,這與之前的研究結(jié)果一致[24-25]。3個降水量區(qū)域的土壤層碳密度分別是油松人工林生物量碳密度的4.13,2.80倍和2.24倍。蘭州地區(qū)(372 mm)生物量碳密度最小,所以土壤碳密度比生物量碳密度高出較多,隨著降水量增加,生物量碳密度增加。

    3.2 油松人工林生物量碳密度與降水量的關(guān)系

    Liu等[26]研究了全球尺度上成熟林的地上生物量碳密度與氣溫和降水量緯度地區(qū)之間的關(guān)系,結(jié)果表明在中緯度地區(qū)降水量1000—2500 mm之間的地區(qū)地上生物量達到最大值。Mehta等[27]研究指出年均降水量是空間尺度上均衡植被覆蓋、人工林及入侵物種等因素引起森林碳儲量變異的主要因素,且在研究4個不同降水量地區(qū)的熱帶森林時年均降水量與地上生物量呈現(xiàn)出正相關(guān)關(guān)系。本研究結(jié)果說明油松人工林的生長與年均降水量之間也呈現(xiàn)較強的線性正相關(guān)關(guān)系。張首軍[28]報道指出年均降水量是影響油松碳固存的主要因子之一,在油松適生的范圍內(nèi)較高的降水量更有利于其積累碳素,這與本研究結(jié)果一致。Chen[29]等研究結(jié)果表明在加拿大西部森林中,年均可利用降水每減少1 cm,演替初期針葉林和演替后期針葉林的地上生物量分別減少0.03 Mg hm-2a-1和0.17 Mg hm-2a-1。Wang[30]等研究報道了青藏高原冷杉林的活生物量碳在海拔梯度上的特征,結(jié)果表明杉木林活生物量碳與年均降水量呈現(xiàn)出正相關(guān)關(guān)系,這與我們的研究結(jié)果一致。

    Kimball等[31]研究指出溫度和降水量的變化可能會引起森林生長季長度的變化。特別是在干旱半干旱區(qū)土壤條件和水分條件是制約人工林生長的重要因素,降水量的增加可能誘導生長季延長,從而提高碳素的固定量。另外,降水量與人工林的凈初級生產(chǎn)力密切相關(guān),特別是在半干旱地區(qū),降水量可能成為影響油松人工林生產(chǎn)力和碳固存的關(guān)鍵因素,所以在降雨量不充足的半干旱地區(qū)應注意采取一些適當?shù)谋K胧﹣肀WC油松林生長所需要的水分,適當?shù)慕?jīng)營管理措施也可提高人工林生態(tài)系統(tǒng)的固碳潛力和綜合服務功能。

    [1] 王效科, 馮宗煒, 歐陽志云. 中國森林生態(tài)系統(tǒng)的植物碳儲量和碳密度研究. 應用生態(tài)學報. 2001, 12(1): 13- 16.

    [2] Moroni M T. Aspects of forest carbon management in Australia: A discussion paper. Forest Ecology and Management, 2012, 275: 111- 116.

    [3] 文仕知, 田大倫, 楊麗麗, 方晰. 榿木人工林的碳密度、碳庫及碳吸存特征. 林業(yè)科學, 2010, 46(6): 15- 21.

    [4] 楊玉姣, 陳云明, 曹揚. 黃土丘陵區(qū)油松人工林生態(tài)系統(tǒng)碳密度及其分配. 生態(tài)學報, 2014, 34(8): 2128- 2136.

    [5] 馮瑞芳, 楊萬勤, 張健. 人工林經(jīng)營與全球變化減緩. 生態(tài)學報, 2006, 26(11): 3870- 3877.

    [6] Li G Q, Xu G H, Guo K, Du S. Geographical boundary and climatic analysis ofPinustabulaeformisin China: Insights on its afforestation. Ecological Engineering, 2016, 86: 75- 84.

    [7] 劉迎春, 王秋鳳, 于貴瑞, 朱先進, 展小云, 郭群, 楊浩, 李勝功, 胡中民. 黃土丘陵區(qū)兩種主要退耕還林樹種生態(tài)系統(tǒng)碳儲量和固碳潛力. 生態(tài)學報, 2011, 31(15): 4277- 4286.

    [8] 張景群, 許喜明, 王曉芳, 王蕾, 張婷婷, 趙曉春. 黃土高原刺槐、油松人工幼林生態(tài)系統(tǒng)碳匯研究. 干旱區(qū)地理, 2011, 34(2): 201- 207.

    [9] 沈彪, 黨坤良, 武朋輝, 朱成功. 秦嶺中段南坡油松林生態(tài)系統(tǒng)碳密度研究. 生態(tài)學報, 2015, 35(6):1- 14.

    [10] 劉冰燕, 陳云明, 曹揚, 吳旭. 秦嶺南坡東段油松人工林生態(tài)系統(tǒng)碳、氮儲量及其分配格局. 應用生態(tài)學報, 2015, 26(3): 643- 652.

    [11] 程小琴, 韓海榮, 康峰峰. 山西油松人工林生態(tài)系統(tǒng)生物量、碳積累及其分布. 生態(tài)學雜志, 2012, 31(10): 2455- 2460.

    [12] 王弋博, 王春燕, 李勃, 楊東. 甘肅省太子山林區(qū)植被恢復措施與土壤酶活性和肥力的相關(guān)性研究. 蘭州大學學報:自然科學版, 2004, 40(2): 92- 96.

    [13] 程堂仁, 馮菁, 馬欽彥, 王玉濤, 康峰峰, 馮忠科, 張彥林, 鄧向瑞. 甘肅小隴山森林植被碳庫及其分配特征. 生態(tài)學報, 2008, 28(1): 33- 44.

    [14] 中國氣象數(shù)據(jù)網(wǎng). 中國地面累年值月值數(shù)據(jù)集(1981- 2010). (2012-08- 28) [2015- 12-04]. http://data.cma.cn/data/detail/dataCode/A.0029.0004.html.

    [15] 中國科學院南京土壤研究所. 土壤理化分析. 上海: 上??萍汲霭嫔? 1978: 132- 140.

    [16] 程堂仁, 馬欽彥, 馮仲科, 羅旭. 甘肅小隴山森林生物量研究. 北京林業(yè)大學學報, 2007, 29(1): 31- 36.

    [17] 李??? 雷淵才, 曾偉生. 基于森林清查資料的中國森林植被碳儲量. 林業(yè)科學, 2011, 47(7): 7- 12.

    [18] 周玉榮, 于振良, 趙士洞. 我國主要森林生態(tài)系統(tǒng)碳貯量和碳平衡. 植物生態(tài)學報, 2000, 24(50): 518- 522.

    [19] Post W M, Kwon K C. Soil carbon sequestration and land-use change: processes and potential. Global Change Biology, 2000, 6(3): 317- 327.

    [20] 李克讓, 王紹強, 曹明奎. 中國植被和土壤碳貯量. 中國科學(D輯), 2003, 33(1): 72- 80.

    [21] 孟蕾, 程積民, 楊曉梅, 韓娟娟, 范文娟, 胡秀娟. 黃土高原子午嶺人工油松林碳儲量與碳密度研究. 水土保持通報, 2010, 30(2): 133- 137.

    [22] 王寧, 王百田, 王瑞君, 曹曉陽, 王文靜, 遲璐. 山西省油松林生態(tài)系統(tǒng)碳密度與分配格局. 應用基礎(chǔ)與工程科學學報, 2014, 22(1): 58- 68.

    [23] 馬欽彥, 謝征鳴. 中國油松林儲碳量基本估計. 北京林業(yè)大學學報, 1996, 18(3): 31- 34.

    [24] Finér L, Mannerkoski H, Piirainen S, Starr M. Carbon and nitrogen pools in an old-growth, Norway spruce mixed forest in eastern Finland and changes associated with clear-cutting. Forest Ecology and Management, 2003, 174(1/3): 51- 63.

    [25] Tandel M B, Kukadia M U, Kolambe B N. Influence of tree cover on physical proper-ties of soil. Indian Forester, 2009, 135(3): 420- 424.

    [26] Liu Y C, Yu G R, Wang Q F, Zhang Y. How temperature, precipitation and stand age control the biomass carbon density of global mature forests. Global Ecology and Biogeography, 2014, 23: 323- 333.

    [27] Mehta N, Pandrya N R, Thomas V O,Krishnayya N S. Impact of rainfall gradient on aboveground biomass and soil organic carbon dynamics of forest covers in Gujarat, India. Ecological research, 2014, 29(6): 1053- 1063.

    [28] 張首軍. 油松純林喬木層碳密度影響因子分析. 山西林業(yè)科技, 2012, 41(3): 14- 17.

    [29] Chen H Y H, Luo Y. Net aboveground biomass declines of four major forest types with forest ageing and climate change in western Canada′s boreal forests. Global Change Biology, 2015, 21(10): 3675- 3684.

    [30] Wang G X, Ran F, Chang R Y, Luo J, Fan J R. Variations in the live biomass and carbon pools ofAbiesgeorgeialong an elevation gradient on the Tibetan Plaateau, China. Forest ecology and Management, 2014, 329: 255- 263.

    [31] Kimball J S, McDonald K C, Running S W, Frolking S E. Satellite radar remote sensing of seasonal growing seasons for boreal and subalpine evergreen forests. Remote Sensing of Environment, 2004, 90(2): 243- 258.

    Carbon storage features ofPinustabulaeformisplantations that are approaching maturity at three sites across a precipitation gradient in western Loess Plateau

    SUN Meimei1, GUAN Jinhong1, 2, WU Chunrong3, YUE Junwei2, LI Guoqing1, 2, DU Sheng1, 2,*

    1StateKeyLaboratoryofErosionandDrylandFarmingonLoessPlateau,NorthwestA&FUniversity,Yangling712100,China2InstituteofSoilandWaterConservation,ChineseAcademyofSciencesandMinistryofWaterResource,Yangling712100,China3GansuDesertControlResearchInstitute,Lanzhou730070,China

    Forest is an important vegetation type of terrestrial ecosystems that maintains the dynamic balance of the biosphere and geosphere, and it plays a key role in terrestrial ecosystem carbon sequestration. With the implementation of Grain for Green Project, the Three-North Shelterbelt Project, and other major ecological projects, the total area of plantations has dramatically increased; thus, the role of plantations in absorbing and fixing CO2has attracted more attention. Investigating carbon storage of major forest types is important for understanding regional carbon cycles.Pinus.tabulaeformisis a common plantation species in sub-humid and semi-arid regions for the restoration of forest ecosystems. In western Loess Plateau, artificial forests ofP.tabulaeformisprovide a substantial proportion of the terrestrial ecosystem carbon sink. To understand the carbon storage features of these plantations in response to the precipitation gradient, three typical distribution regions on western Loess Plateau were selected for plot surveys. Tree biomass was estimated separately for stems, branches, leaves, and roots, using previously constructed allometric biomass equations based on the diameter at breast height (DBH) and tree height. Carbon contents in the different parts and layers of the tree were each measured by collecting corresponding samples. Carbon density was calculated by multiplying the biomass of each plant part (from a specific forest area) and the corresponding carbon content. The results showed that in the three sites (Lanzhou, Taizishan, Xiaolongshan) with precipitation of 372 mm, 519 mm and 632 mm, respectively, biomass carbon densities ofP.tabulaeformisplantations were (49.08±2.86) t/hm2, (73.90±9.36) t/hm2, and (82.55±7.36) t/hm2, respectively. Both the carbon density of the ecosystem and biomass carbon density were significantly different (P<0.05) between Xiaolongshan and Lanzhou. Because the shrubs were sparse, the contribution of trees was greatest, reaching 95.21%—98.91% of the biomass carbon density across the precipitation gradient. Overall density of organic carbon in soil did not differ between the three sites. Soil organic carbon in the upper soil layers showed significant differences among the sites, but were not statistically different (P>0.05) in other soil layers. Correlative analysis showed that the biomass carbon density was positively correlated with annual precipitation, with a Pearson coefficient of 0.820. However, there were no significant positive correlations between the biomass carbon density and the age of stand. Similarly, partial correlative analysis showed that there was a significant positive correlation between the carbon density ofP.tabulaeformisplantations and precipitation in this sub-humid and semi-arid region of the Loess Plateau. These results indicate that precipitation is a key factor that affects biomass production and carbon fixation in semi-arid forests. Furthermore, it implies that appropriate water conservation measures are necessary, particularly in semi-arid areas, to satisfy the hydric demands ofP.tabulaeformisplantations and to maintain normal growth. Such measures may also be helpful for improving carbon sequestration potential and comprehensive services of plantation ecosystems. Therefore, this study provides valuable insight into how carbon storage characteristics of plantations respond to precipitation gradients, and provides useful information for regional forest management with respect to productivity and carbon storage.

    Pinustabulaeformis;biomass carbon density; precipitation; western Loess Plateau

    中國科學院戰(zhàn)略性先導科技專項資助項目(XDA05050202)

    2015- 12- 23; 網(wǎng)絡出版日期:2016- 10- 29

    10.5846/stxb201512232562

    *通訊作者Corresponding author.E-mail: shengdu@ms.iswc.ac.cn

    孫美美,關(guān)晉宏,吳春榮,岳軍偉,李國慶,杜盛.黃土高原西部3個降水量梯度近成熟油松人工林碳庫特征.生態(tài)學報,2017,37(8):2665- 2672.

    Sun M M, Guan J H, Wu C R, Yue J W, Li G Q, Du S.Carbon storage features ofPinustabulaeformisplantations that are approaching maturity at three sites across a precipitation gradient in western Loess Plateau.Acta Ecologica Sinica,2017,37(8):2665- 2672.

    猜你喜歡
    碳庫油松人工林
    油松容器育苗技術(shù)
    長期定位試驗下砒砂巖與沙復配土的碳庫管理指數(shù)
    綠色科技(2020年20期)2020-11-20 01:56:34
    北方油松育苗技術(shù)
    我國人工林生態(tài)輪伐期探討
    人工林生產(chǎn)力提升對策探析
    秸稈還田對農(nóng)田土壤碳庫和溫室氣體排放的影響研究進展
    大氣氮沉降對森林土壤碳庫的影響
    山東林業(yè)科技(2016年5期)2016-07-05 00:43:04
    濕地松人工林生長規(guī)律
    油松種植技術(shù)及主要病蟲害的防治對策
    欧美日本亚洲视频在线播放| 女人十人毛片免费观看3o分钟| 亚洲人成网站在线播| 国产综合懂色| 天天一区二区日本电影三级| 不卡一级毛片| 亚洲美女视频黄频| 最近中文字幕高清免费大全6| 最近的中文字幕免费完整| 亚洲人成网站高清观看| 亚洲av二区三区四区| 久久99精品国语久久久| 有码 亚洲区| 69av精品久久久久久| 麻豆乱淫一区二区| 久久久久久久亚洲中文字幕| 婷婷精品国产亚洲av| 2022亚洲国产成人精品| 又粗又硬又长又爽又黄的视频 | 直男gayav资源| 国产成人a区在线观看| 免费无遮挡裸体视频| 国产精品久久久久久亚洲av鲁大| 日韩欧美三级三区| 天天躁日日操中文字幕| 国产精品久久久久久久久免| 在线国产一区二区在线| 国产精品综合久久久久久久免费| 国产精品一及| 欧美在线一区亚洲| 亚洲一区高清亚洲精品| 赤兔流量卡办理| 亚洲中文字幕一区二区三区有码在线看| 成人国产麻豆网| 少妇裸体淫交视频免费看高清| 精品人妻熟女av久视频| 国产成人精品一,二区 | 免费av毛片视频| 联通29元200g的流量卡| 久久精品国产99精品国产亚洲性色| 麻豆av噜噜一区二区三区| 免费看日本二区| 精品日产1卡2卡| 晚上一个人看的免费电影| 干丝袜人妻中文字幕| 亚洲人成网站在线播| 一级毛片aaaaaa免费看小| 国产精品99久久久久久久久| 夫妻性生交免费视频一级片| 91久久精品电影网| 精品人妻视频免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品人妻偷拍中文字幕| 日本色播在线视频| 变态另类丝袜制服| 久久韩国三级中文字幕| 亚洲精品国产av成人精品| 日韩 亚洲 欧美在线| 又爽又黄无遮挡网站| 亚洲国产精品久久男人天堂| 精品久久久久久久人妻蜜臀av| 欧美变态另类bdsm刘玥| 欧美变态另类bdsm刘玥| 老司机影院成人| 性色avwww在线观看| 国产片特级美女逼逼视频| 午夜激情欧美在线| 欧美成人精品欧美一级黄| 国产成人a区在线观看| 伦理电影大哥的女人| 国产av在哪里看| 欧美性猛交黑人性爽| 国产片特级美女逼逼视频| 中文字幕制服av| 99久久成人亚洲精品观看| 免费人成在线观看视频色| 国产精品.久久久| 大型黄色视频在线免费观看| 国产蜜桃级精品一区二区三区| 狂野欧美激情性xxxx在线观看| 一级毛片我不卡| 久久久国产成人精品二区| 午夜视频国产福利| 99久久九九国产精品国产免费| 国产三级中文精品| 国产男人的电影天堂91| 久久午夜福利片| 精品人妻熟女av久视频| 亚洲天堂国产精品一区在线| 亚洲美女搞黄在线观看| av在线亚洲专区| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 热99re8久久精品国产| 看黄色毛片网站| 日本成人三级电影网站| 日本黄色视频三级网站网址| 国产精品无大码| 中文欧美无线码| 日本五十路高清| 亚洲成人久久性| 色5月婷婷丁香| 亚洲最大成人中文| 亚洲国产日韩欧美精品在线观看| 亚洲图色成人| 1024手机看黄色片| 一本久久中文字幕| 18禁裸乳无遮挡免费网站照片| 午夜亚洲福利在线播放| 99久久无色码亚洲精品果冻| 久久韩国三级中文字幕| 国产成人影院久久av| 成年免费大片在线观看| 国产精品一二三区在线看| 麻豆精品久久久久久蜜桃| 高清在线视频一区二区三区 | h日本视频在线播放| 国产伦精品一区二区三区四那| 成人综合一区亚洲| 精品久久久噜噜| 国内精品久久久久精免费| 老女人水多毛片| 日韩欧美国产在线观看| 国产亚洲av片在线观看秒播厂 | 亚洲va在线va天堂va国产| 久久6这里有精品| 卡戴珊不雅视频在线播放| 国产精品久久久久久av不卡| 亚洲欧美日韩卡通动漫| 国产精品美女特级片免费视频播放器| 国产麻豆成人av免费视频| 亚洲国产精品成人综合色| 精品免费久久久久久久清纯| 亚洲av第一区精品v没综合| 亚洲国产精品成人综合色| 夜夜爽天天搞| 中文精品一卡2卡3卡4更新| 波多野结衣高清作品| 五月玫瑰六月丁香| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲av嫩草精品影院| 久久精品国产99精品国产亚洲性色| 免费无遮挡裸体视频| 国产一区二区三区av在线 | 亚洲av中文字字幕乱码综合| 九九在线视频观看精品| 2021天堂中文幕一二区在线观| 成人永久免费在线观看视频| 成人漫画全彩无遮挡| 我的女老师完整版在线观看| 欧美精品国产亚洲| 久久久精品欧美日韩精品| 成人午夜高清在线视频| a级毛片免费高清观看在线播放| 一边亲一边摸免费视频| 五月玫瑰六月丁香| 亚洲婷婷狠狠爱综合网| 能在线免费看毛片的网站| 高清毛片免费观看视频网站| 麻豆国产av国片精品| 成人性生交大片免费视频hd| 99久久人妻综合| 国产av麻豆久久久久久久| 亚洲中文字幕日韩| av免费观看日本| 一本久久精品| 亚洲国产精品国产精品| 色综合色国产| 国产精品99久久久久久久久| 国产精品.久久久| 亚洲成人精品中文字幕电影| 久久人人精品亚洲av| 亚洲欧美成人综合另类久久久 | 成人亚洲欧美一区二区av| 久久精品91蜜桃| 国产高潮美女av| 久久国产乱子免费精品| 国产伦精品一区二区三区视频9| 深夜精品福利| 2022亚洲国产成人精品| av视频在线观看入口| 又爽又黄无遮挡网站| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 91av网一区二区| 亚洲最大成人av| 久久久精品94久久精品| 69人妻影院| 在线免费观看不下载黄p国产| 午夜福利成人在线免费观看| av黄色大香蕉| 成人午夜精彩视频在线观看| av在线亚洲专区| 日韩国内少妇激情av| 亚洲丝袜综合中文字幕| 直男gayav资源| 免费人成在线观看视频色| 免费av观看视频| 国产高清不卡午夜福利| 搡女人真爽免费视频火全软件| 久久久久久久久中文| 国产日本99.免费观看| 国产日本99.免费观看| 嘟嘟电影网在线观看| 国产av一区在线观看免费| 亚洲人成网站高清观看| 一区二区三区高清视频在线| 国产精品免费一区二区三区在线| 日韩中字成人| 超碰av人人做人人爽久久| av女优亚洲男人天堂| 联通29元200g的流量卡| 成人特级黄色片久久久久久久| 国产私拍福利视频在线观看| 欧美成人一区二区免费高清观看| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜精品论理片| 男的添女的下面高潮视频| 久久这里只有精品中国| 搞女人的毛片| 免费无遮挡裸体视频| 亚洲精品日韩av片在线观看| 亚洲自偷自拍三级| 91在线精品国自产拍蜜月| 麻豆乱淫一区二区| 日韩成人伦理影院| 91久久精品电影网| 久久亚洲精品不卡| 亚洲精品影视一区二区三区av| 男人狂女人下面高潮的视频| 国产在视频线在精品| 别揉我奶头 嗯啊视频| 最近最新中文字幕大全电影3| 成人综合一区亚洲| 白带黄色成豆腐渣| 成人毛片60女人毛片免费| 色噜噜av男人的天堂激情| 国产精品久久久久久亚洲av鲁大| 日日啪夜夜撸| 日本与韩国留学比较| 国产高清视频在线观看网站| 日韩欧美精品免费久久| 日本撒尿小便嘘嘘汇集6| 国产av一区在线观看免费| 亚洲欧美日韩高清在线视频| 久久久久免费精品人妻一区二区| 最近的中文字幕免费完整| 午夜爱爱视频在线播放| 性欧美人与动物交配| 校园春色视频在线观看| 啦啦啦观看免费观看视频高清| 中文字幕人妻熟人妻熟丝袜美| 国产综合懂色| 波野结衣二区三区在线| 中文精品一卡2卡3卡4更新| 人妻少妇偷人精品九色| 精品久久久噜噜| 少妇丰满av| 久久久久久国产a免费观看| 麻豆精品久久久久久蜜桃| 午夜免费男女啪啪视频观看| 色吧在线观看| 国产精品一区www在线观看| .国产精品久久| 国产亚洲91精品色在线| 天天躁日日操中文字幕| 国产伦一二天堂av在线观看| 一进一出抽搐gif免费好疼| 听说在线观看完整版免费高清| 日本av手机在线免费观看| 久久午夜福利片| 国产精品久久视频播放| 禁无遮挡网站| 国产一级毛片在线| 中文字幕av在线有码专区| 日本熟妇午夜| 天天躁夜夜躁狠狠久久av| 美女大奶头视频| 国产成人91sexporn| av国产免费在线观看| 精品国产三级普通话版| 国产男人的电影天堂91| 亚洲av成人av| 免费观看人在逋| 最近中文字幕高清免费大全6| 最近手机中文字幕大全| 日本免费a在线| 国产成人影院久久av| 久久草成人影院| 亚洲丝袜综合中文字幕| 国产成人精品久久久久久| 国产成年人精品一区二区| 国产片特级美女逼逼视频| 天堂av国产一区二区熟女人妻| 麻豆久久精品国产亚洲av| 在线免费观看的www视频| 插阴视频在线观看视频| 十八禁国产超污无遮挡网站| www日本黄色视频网| 亚洲五月天丁香| 久久99热6这里只有精品| 十八禁国产超污无遮挡网站| 亚洲不卡免费看| 国产乱人偷精品视频| 成人欧美大片| 亚洲国产高清在线一区二区三| 精品不卡国产一区二区三区| 婷婷色综合大香蕉| 亚洲精品粉嫩美女一区| 国产精品女同一区二区软件| 在线免费十八禁| 天堂av国产一区二区熟女人妻| 美女黄网站色视频| 一本一本综合久久| 蜜桃久久精品国产亚洲av| av在线老鸭窝| 日韩视频在线欧美| 熟女电影av网| 九色成人免费人妻av| 热99re8久久精品国产| 老女人水多毛片| 国产高清三级在线| 校园春色视频在线观看| 精品人妻熟女av久视频| 国产又黄又爽又无遮挡在线| av福利片在线观看| 床上黄色一级片| 少妇裸体淫交视频免费看高清| а√天堂www在线а√下载| 国产成人aa在线观看| 久久久久久久久久成人| 国产精品1区2区在线观看.| 亚洲激情五月婷婷啪啪| 哪里可以看免费的av片| 性色avwww在线观看| 国产白丝娇喘喷水9色精品| 久久久久久久久久黄片| 免费黄网站久久成人精品| 国产毛片a区久久久久| 久久精品人妻少妇| 久久久久久国产a免费观看| 亚洲五月天丁香| 成年av动漫网址| 国产精品99久久久久久久久| 99视频精品全部免费 在线| 男女下面进入的视频免费午夜| 午夜福利视频1000在线观看| 亚洲人成网站高清观看| 国产精品一区二区在线观看99 | 国产精品久久久久久av不卡| av在线老鸭窝| 只有这里有精品99| 97热精品久久久久久| 午夜精品国产一区二区电影 | 欧美日韩综合久久久久久| www.av在线官网国产| 五月伊人婷婷丁香| 欧美区成人在线视频| 久久热精品热| 欧美成人一区二区免费高清观看| 免费观看精品视频网站| 亚洲欧美日韩无卡精品| 亚洲人与动物交配视频| 一级黄色大片毛片| 国产精品久久久久久av不卡| 欧美高清成人免费视频www| 亚州av有码| 欧美色欧美亚洲另类二区| 三级男女做爰猛烈吃奶摸视频| 久久久精品94久久精品| 成人鲁丝片一二三区免费| 三级男女做爰猛烈吃奶摸视频| 一区二区三区四区激情视频 | 久久久久久久午夜电影| 人妻久久中文字幕网| 亚洲色图av天堂| 免费人成视频x8x8入口观看| 久久久欧美国产精品| 91久久精品国产一区二区三区| 天堂影院成人在线观看| 3wmmmm亚洲av在线观看| 少妇人妻精品综合一区二区 | 免费电影在线观看免费观看| 卡戴珊不雅视频在线播放| 亚洲欧美日韩高清在线视频| 日本五十路高清| 久久精品影院6| 啦啦啦韩国在线观看视频| 国产一区二区在线观看日韩| 精品人妻熟女av久视频| 欧美日本视频| 一边亲一边摸免费视频| 91精品国产九色| 成人特级黄色片久久久久久久| 爱豆传媒免费全集在线观看| 桃色一区二区三区在线观看| 国产色爽女视频免费观看| 国产成人a区在线观看| 99久久人妻综合| 午夜精品一区二区三区免费看| 99精品在免费线老司机午夜| 看非洲黑人一级黄片| 免费搜索国产男女视频| 欧美xxxx性猛交bbbb| 午夜福利成人在线免费观看| 亚洲在线自拍视频| 免费黄网站久久成人精品| 熟妇人妻久久中文字幕3abv| 国产av在哪里看| 一个人观看的视频www高清免费观看| 国产欧美日韩精品一区二区| 最新中文字幕久久久久| 一级毛片久久久久久久久女| 国产午夜精品论理片| 亚州av有码| 欧美丝袜亚洲另类| 在线观看66精品国产| 91久久精品国产一区二区三区| 校园春色视频在线观看| 久久婷婷人人爽人人干人人爱| 青春草国产在线视频 | 亚洲精品日韩av片在线观看| 国产老妇伦熟女老妇高清| 欧美日韩国产亚洲二区| 日韩av在线大香蕉| 日韩大尺度精品在线看网址| 国产极品天堂在线| 亚洲精品久久国产高清桃花| 在线播放国产精品三级| 国产真实伦视频高清在线观看| 一边摸一边抽搐一进一小说| 国产午夜福利久久久久久| 日韩高清综合在线| 麻豆成人av视频| 久久久久久久久久黄片| 亚洲av电影不卡..在线观看| 美女脱内裤让男人舔精品视频 | 亚洲精品乱码久久久v下载方式| 男人舔奶头视频| 久久久国产成人精品二区| 伦精品一区二区三区| АⅤ资源中文在线天堂| 老女人水多毛片| 国产一区二区在线观看日韩| 乱码一卡2卡4卡精品| 男女啪啪激烈高潮av片| 一级av片app| 国产一级毛片七仙女欲春2| 亚洲国产色片| 永久网站在线| 国产真实乱freesex| 在线观看av片永久免费下载| 国产老妇伦熟女老妇高清| 少妇人妻一区二区三区视频| 亚洲欧美日韩卡通动漫| 麻豆一二三区av精品| 免费不卡的大黄色大毛片视频在线观看 | 久久综合国产亚洲精品| 久久草成人影院| 国产精品久久久久久久电影| 一个人观看的视频www高清免费观看| 免费黄网站久久成人精品| 国产黄片美女视频| av国产免费在线观看| 日韩av在线大香蕉| 亚洲图色成人| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 久久人人爽人人片av| 国产午夜精品论理片| 综合色av麻豆| 成熟少妇高潮喷水视频| 少妇高潮的动态图| 一级毛片久久久久久久久女| 成年av动漫网址| 日本黄色视频三级网站网址| 国国产精品蜜臀av免费| 国产亚洲av片在线观看秒播厂 | 久久热精品热| 中文字幕熟女人妻在线| 99热这里只有是精品在线观看| www.av在线官网国产| 高清日韩中文字幕在线| 少妇熟女欧美另类| 国产一级毛片七仙女欲春2| 99热全是精品| 国产免费一级a男人的天堂| 夜夜夜夜夜久久久久| 午夜福利高清视频| 三级毛片av免费| 亚洲一区二区三区色噜噜| 在线免费观看不下载黄p国产| 午夜福利在线观看免费完整高清在 | 精品人妻一区二区三区麻豆| 国产高清激情床上av| 国产人妻一区二区三区在| 12—13女人毛片做爰片一| 亚洲av熟女| 久久久精品大字幕| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩卡通动漫| 又爽又黄无遮挡网站| 亚洲无线观看免费| 国产极品天堂在线| 色综合站精品国产| 综合色丁香网| 在现免费观看毛片| 极品教师在线视频| 久久鲁丝午夜福利片| 激情 狠狠 欧美| 欧美成人精品欧美一级黄| av.在线天堂| 国产精品爽爽va在线观看网站| 久久精品久久久久久久性| 日本黄色视频三级网站网址| 一级毛片久久久久久久久女| 亚洲在线自拍视频| 亚洲18禁久久av| 国产单亲对白刺激| 特级一级黄色大片| 色噜噜av男人的天堂激情| 69人妻影院| 日本撒尿小便嘘嘘汇集6| 国产av麻豆久久久久久久| 久久精品人妻少妇| 97在线视频观看| 久久精品夜色国产| 可以在线观看毛片的网站| 1000部很黄的大片| 丰满人妻一区二区三区视频av| 午夜精品在线福利| 在线观看66精品国产| 床上黄色一级片| 欧美日韩国产亚洲二区| 国产一区二区三区在线臀色熟女| 国产精品人妻久久久久久| 亚洲精品影视一区二区三区av| 亚洲国产精品合色在线| 69av精品久久久久久| 在线国产一区二区在线| 国产精品爽爽va在线观看网站| 中文字幕久久专区| 草草在线视频免费看| 99riav亚洲国产免费| 国产伦理片在线播放av一区 | 男插女下体视频免费在线播放| 老女人水多毛片| 91av网一区二区| 久久综合国产亚洲精品| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 1000部很黄的大片| 色播亚洲综合网| 国内精品久久久久精免费| 久久99精品国语久久久| 午夜福利视频1000在线观看| 中国美白少妇内射xxxbb| 99热只有精品国产| 少妇人妻精品综合一区二区 | 一个人看视频在线观看www免费| 在线免费观看的www视频| 精品无人区乱码1区二区| 三级经典国产精品| 在线观看66精品国产| 久久99热这里只有精品18| 99热这里只有是精品在线观看| 亚洲中文字幕日韩| 午夜福利视频1000在线观看| 国产女主播在线喷水免费视频网站 | 亚洲国产欧美在线一区| 欧美极品一区二区三区四区| 日韩制服骚丝袜av| 在线免费十八禁| 亚洲精品亚洲一区二区| 婷婷色av中文字幕| 99视频精品全部免费 在线| 超碰av人人做人人爽久久| 嫩草影院入口| 国产探花在线观看一区二区| 国内揄拍国产精品人妻在线| 综合色丁香网| 成人永久免费在线观看视频| 成年女人永久免费观看视频| 91狼人影院| 毛片一级片免费看久久久久| 在线播放无遮挡| 日日摸夜夜添夜夜爱| 全区人妻精品视频| 国产成人aa在线观看| 国产精品无大码| 蜜桃久久精品国产亚洲av| 波多野结衣高清作品| 中文字幕久久专区| 综合色av麻豆| 真实男女啪啪啪动态图| 欧美精品国产亚洲| 国产亚洲精品av在线| 99久久中文字幕三级久久日本| 国产v大片淫在线免费观看| 男的添女的下面高潮视频| 热99re8久久精品国产| 高清午夜精品一区二区三区 | 高清毛片免费观看视频网站| 自拍偷自拍亚洲精品老妇| 极品教师在线视频| 一级毛片电影观看 | 亚洲欧美清纯卡通| 神马国产精品三级电影在线观看| 国产麻豆成人av免费视频| 一区福利在线观看| 青春草视频在线免费观看| 综合色av麻豆| 中国美白少妇内射xxxbb| 久久午夜福利片| 亚洲国产精品sss在线观看|