• <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ù)及主要病蟲害的防治對策
    人妻一区二区av| 国产精品一区二区在线观看99| 午夜免费鲁丝| 一级爰片在线观看| 日韩制服丝袜自拍偷拍| 欧美成人精品欧美一级黄| 日韩一区二区三区影片| 老司机亚洲免费影院| 成人国产av品久久久| 美女国产高潮福利片在线看| 午夜福利视频精品| 欧美另类一区| 狂野欧美激情性xxxx在线观看| 国产精品一区二区在线不卡| 王馨瑶露胸无遮挡在线观看| 麻豆精品久久久久久蜜桃| 国产亚洲av片在线观看秒播厂| 巨乳人妻的诱惑在线观看| 国产精品不卡视频一区二区| 久久av网站| 91精品三级在线观看| 伦精品一区二区三区| 自线自在国产av| 99久久中文字幕三级久久日本| 一级毛片 在线播放| 一级毛片电影观看| 午夜免费鲁丝| 国产1区2区3区精品| av卡一久久| 桃花免费在线播放| 两个人免费观看高清视频| 中文字幕免费在线视频6| 狠狠精品人妻久久久久久综合| 国产精品蜜桃在线观看| 国产一级毛片在线| 久久久久久久国产电影| 亚洲国产成人一精品久久久| 亚洲欧美日韩卡通动漫| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区黑人 | 国产在线免费精品| 母亲3免费完整高清在线观看 | 免费女性裸体啪啪无遮挡网站| 国产黄色免费在线视频| 日产精品乱码卡一卡2卡三| 男女边摸边吃奶| 日韩制服骚丝袜av| 搡老乐熟女国产| 又黄又粗又硬又大视频| 亚洲国产欧美日韩在线播放| 欧美亚洲日本最大视频资源| 99国产综合亚洲精品| 国产伦理片在线播放av一区| 99精国产麻豆久久婷婷| av网站免费在线观看视频| 大香蕉久久成人网| 菩萨蛮人人尽说江南好唐韦庄| 我的女老师完整版在线观看| 国产成人一区二区在线| av视频免费观看在线观看| 我的女老师完整版在线观看| 成人国语在线视频| 国产日韩一区二区三区精品不卡| 黄片播放在线免费| 国产精品久久久久久av不卡| 亚洲国产色片| 激情五月婷婷亚洲| 成人毛片60女人毛片免费| 男女边摸边吃奶| 欧美bdsm另类| 一边亲一边摸免费视频| 国产成人a∨麻豆精品| 日本欧美国产在线视频| 丝袜脚勾引网站| 黄色 视频免费看| 十八禁高潮呻吟视频| 午夜福利乱码中文字幕| 欧美日本中文国产一区发布| 高清视频免费观看一区二区| 久久久久久久久久成人| 国产成人精品福利久久| 精品午夜福利在线看| 精品久久久久久电影网| 各种免费的搞黄视频| 午夜久久久在线观看| 97在线视频观看| 欧美日韩国产mv在线观看视频| 亚洲精品日本国产第一区| 蜜桃国产av成人99| 精品一区二区三卡| 国精品久久久久久国模美| 日本猛色少妇xxxxx猛交久久| 啦啦啦在线观看免费高清www| 国产精品欧美亚洲77777| 美女内射精品一级片tv| 欧美国产精品va在线观看不卡| 啦啦啦啦在线视频资源| 中文字幕制服av| 极品少妇高潮喷水抽搐| 国产精品麻豆人妻色哟哟久久| a级片在线免费高清观看视频| 免费看光身美女| av网站免费在线观看视频| 极品人妻少妇av视频| av国产精品久久久久影院| 国产精品99久久99久久久不卡 | 少妇人妻久久综合中文| 欧美日韩视频高清一区二区三区二| 国产老妇伦熟女老妇高清| 在线观看人妻少妇| 黄色一级大片看看| 国产黄色免费在线视频| 国产精品欧美亚洲77777| 国产老妇伦熟女老妇高清| 91成人精品电影| 国产精品人妻久久久久久| 超碰97精品在线观看| 日本欧美视频一区| 亚洲av日韩在线播放| 国产一区亚洲一区在线观看| 大香蕉97超碰在线| 夜夜爽夜夜爽视频| 大话2 男鬼变身卡| 五月玫瑰六月丁香| 亚洲综合精品二区| 999精品在线视频| 99re6热这里在线精品视频| 大香蕉久久成人网| 日本爱情动作片www.在线观看| 欧美 日韩 精品 国产| 日韩熟女老妇一区二区性免费视频| 一级毛片黄色毛片免费观看视频| 最近最新中文字幕免费大全7| 99久久精品国产国产毛片| 国产一区有黄有色的免费视频| www.av在线官网国产| 日韩伦理黄色片| 亚洲欧洲精品一区二区精品久久久 | av一本久久久久| 亚洲婷婷狠狠爱综合网| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 日韩制服骚丝袜av| 午夜福利视频精品| 高清在线视频一区二区三区| 看十八女毛片水多多多| 观看av在线不卡| 日日摸夜夜添夜夜爱| 中文字幕制服av| 18禁在线无遮挡免费观看视频| 免费日韩欧美在线观看| 国产又爽黄色视频| av.在线天堂| 国产成人免费观看mmmm| 欧美性感艳星| 国产成人aa在线观看| 成年女人在线观看亚洲视频| 久久99一区二区三区| 精品一区二区免费观看| 精品少妇内射三级| 国产乱来视频区| 国产爽快片一区二区三区| 精品人妻熟女毛片av久久网站| 美女中出高潮动态图| 最近中文字幕高清免费大全6| 多毛熟女@视频| av女优亚洲男人天堂| 国产色爽女视频免费观看| 免费人成在线观看视频色| 国产精品久久久久久久电影| 少妇的丰满在线观看| 少妇 在线观看| 日韩在线高清观看一区二区三区| 国产成人精品无人区| 大片电影免费在线观看免费| 亚洲欧美中文字幕日韩二区| 精品一品国产午夜福利视频| 男的添女的下面高潮视频| 午夜激情久久久久久久| 午夜免费鲁丝| 晚上一个人看的免费电影| 性高湖久久久久久久久免费观看| 欧美bdsm另类| 一本大道久久a久久精品| 亚洲三级黄色毛片| 日韩av在线免费看完整版不卡| 久久午夜综合久久蜜桃| 亚洲av综合色区一区| 亚洲中文av在线| 日本与韩国留学比较| 有码 亚洲区| 在线观看国产h片| 一区二区三区精品91| 日韩精品有码人妻一区| 国产免费视频播放在线视频| 黄色 视频免费看| 亚洲国产日韩一区二区| 九色成人免费人妻av| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 国产在线免费精品| www.色视频.com| 99久久综合免费| 天堂俺去俺来也www色官网| 赤兔流量卡办理| 在线天堂最新版资源| 婷婷色综合大香蕉| 又粗又硬又长又爽又黄的视频| 日本av手机在线免费观看| 国产精品99久久99久久久不卡 | 成人免费观看视频高清| 视频在线观看一区二区三区| 爱豆传媒免费全集在线观看| 人妻人人澡人人爽人人| 男女午夜视频在线观看 | 美国免费a级毛片| 新久久久久国产一级毛片| 国产有黄有色有爽视频| 亚洲一级一片aⅴ在线观看| 国产精品久久久久成人av| 丁香六月天网| 国产免费福利视频在线观看| 日韩av免费高清视频| 狂野欧美激情性bbbbbb| a级片在线免费高清观看视频| 亚洲熟女精品中文字幕| 午夜福利视频精品| 99re6热这里在线精品视频| 亚洲综合色网址| 国产一级毛片在线| freevideosex欧美| 国产片特级美女逼逼视频| 亚洲中文av在线| 国产亚洲午夜精品一区二区久久| 久久久国产欧美日韩av| 成人国产麻豆网| 亚洲五月色婷婷综合| 大陆偷拍与自拍| 成人漫画全彩无遮挡| 韩国精品一区二区三区 | 亚洲伊人色综图| 国产国拍精品亚洲av在线观看| 18禁在线无遮挡免费观看视频| 高清毛片免费看| 久久久久久久久久久免费av| 观看av在线不卡| 久久青草综合色| 日日啪夜夜爽| 欧美97在线视频| 久久久久精品人妻al黑| 成人毛片a级毛片在线播放| 丝袜人妻中文字幕| 色哟哟·www| 男人爽女人下面视频在线观看| 精品亚洲乱码少妇综合久久| 欧美激情极品国产一区二区三区 | 少妇人妻精品综合一区二区| 亚洲中文av在线| 有码 亚洲区| 三上悠亚av全集在线观看| 高清黄色对白视频在线免费看| 久久热在线av| 九色亚洲精品在线播放| 亚洲久久久国产精品| 我的女老师完整版在线观看| 99热全是精品| av有码第一页| 亚洲熟女精品中文字幕| 亚洲综合精品二区| 80岁老熟妇乱子伦牲交| 国产极品天堂在线| 精品卡一卡二卡四卡免费| 国产又爽黄色视频| 国产 一区精品| 最后的刺客免费高清国语| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站| 成人影院久久| 久久精品久久久久久久性| www.av在线官网国产| 美女国产视频在线观看| 在线 av 中文字幕| 国产色婷婷99| 免费观看在线日韩| 多毛熟女@视频| 黄色 视频免费看| 99视频精品全部免费 在线| 中国三级夫妇交换| 美女国产高潮福利片在线看| 久久亚洲国产成人精品v| 99香蕉大伊视频| 国产亚洲精品久久久com| 久久狼人影院| av播播在线观看一区| av国产精品久久久久影院| 一个人免费看片子| 亚洲av国产av综合av卡| 欧美性感艳星| 成人二区视频| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 国产精品久久久久成人av| 热99久久久久精品小说推荐| 国产精品国产三级专区第一集| 啦啦啦中文免费视频观看日本| av在线app专区| 伊人亚洲综合成人网| 国产综合精华液| 国产精品久久久久久久久免| 九九在线视频观看精品| 激情五月婷婷亚洲| 久久久久久久国产电影| 狠狠婷婷综合久久久久久88av| 多毛熟女@视频| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| 免费av中文字幕在线| 香蕉精品网在线| 久久国产精品大桥未久av| av一本久久久久| 九九爱精品视频在线观看| 精品熟女少妇av免费看| 精品第一国产精品| 久久久a久久爽久久v久久| 只有这里有精品99| 午夜福利影视在线免费观看| 亚洲人成网站在线观看播放| 99久久中文字幕三级久久日本| 中文字幕制服av| 免费久久久久久久精品成人欧美视频 | 有码 亚洲区| 在线亚洲精品国产二区图片欧美| 寂寞人妻少妇视频99o| 黑人欧美特级aaaaaa片| 丝袜喷水一区| 婷婷色麻豆天堂久久| 啦啦啦视频在线资源免费观看| 在线观看人妻少妇| a级毛片在线看网站| 一二三四中文在线观看免费高清| 国产xxxxx性猛交| av免费观看日本| 在线 av 中文字幕| 免费大片黄手机在线观看| 久久人人爽人人爽人人片va| 香蕉精品网在线| 亚洲精品aⅴ在线观看| 一区二区三区四区激情视频| 一级爰片在线观看| 久久99蜜桃精品久久| 中文字幕精品免费在线观看视频 | 久久影院123| 久久精品国产a三级三级三级| 久久免费观看电影| 美女内射精品一级片tv| 国产不卡av网站在线观看| 一级毛片 在线播放| 天堂俺去俺来也www色官网| 国产精品 国内视频| 国产 精品1| 久久久久久久久久久久大奶| 最近2019中文字幕mv第一页| av在线观看视频网站免费| 韩国av在线不卡| 国产黄色免费在线视频| 亚洲成人一二三区av| 成人手机av| 成人二区视频| 99香蕉大伊视频| 女的被弄到高潮叫床怎么办| 亚洲国产av影院在线观看| 高清毛片免费看| 国产一区二区在线观看av| 国产爽快片一区二区三区| 国产极品粉嫩免费观看在线| av黄色大香蕉| 国产成人精品婷婷| 国产又色又爽无遮挡免| 天美传媒精品一区二区| videosex国产| av电影中文网址| 日韩熟女老妇一区二区性免费视频| 伦理电影免费视频| 久久国产精品男人的天堂亚洲 | 男女午夜视频在线观看 | 97精品久久久久久久久久精品| 黑丝袜美女国产一区| 夫妻午夜视频| 天堂俺去俺来也www色官网| 在线天堂最新版资源| 又粗又硬又长又爽又黄的视频| 你懂的网址亚洲精品在线观看| 日韩不卡一区二区三区视频在线| 国产精品国产三级国产专区5o| 狠狠婷婷综合久久久久久88av| 中国国产av一级| 欧美精品高潮呻吟av久久| 亚洲欧洲国产日韩| 亚洲av综合色区一区| 免费人妻精品一区二区三区视频| 成年人午夜在线观看视频| 黄网站色视频无遮挡免费观看| 男人舔女人的私密视频| 亚洲一码二码三码区别大吗| 久久久久视频综合| 天天躁夜夜躁狠狠躁躁| 男女免费视频国产| 天堂8中文在线网| 精品视频人人做人人爽| 两性夫妻黄色片 | 国语对白做爰xxxⅹ性视频网站| 国产成人精品在线电影| 色视频在线一区二区三区| 国产免费又黄又爽又色| 精品久久久久久电影网| 久久人人97超碰香蕉20202| 成人手机av| 久久久久网色| 久久久久久久久久久免费av| 黄片播放在线免费| 人妻一区二区av| 亚洲情色 制服丝袜| 免费看不卡的av| www.熟女人妻精品国产 | 国产成人91sexporn| 建设人人有责人人尽责人人享有的| 美女大奶头黄色视频| 啦啦啦视频在线资源免费观看| av视频免费观看在线观看| 晚上一个人看的免费电影| 成年美女黄网站色视频大全免费| 成人无遮挡网站| 在线天堂中文资源库| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 99热网站在线观看| 亚洲av综合色区一区| 国产黄频视频在线观看| 国产 精品1| 亚洲av电影在线进入| 青青草视频在线视频观看| 青春草视频在线免费观看| 久久这里只有精品19| 婷婷成人精品国产| 2022亚洲国产成人精品| 久久久久网色| 全区人妻精品视频| 国产免费一区二区三区四区乱码| 成年人免费黄色播放视频| 成人漫画全彩无遮挡| 成人午夜精彩视频在线观看| 亚洲国产欧美在线一区| 男女午夜视频在线观看 | 国产 精品1| 午夜福利在线观看免费完整高清在| 精品卡一卡二卡四卡免费| 黄色配什么色好看| 黄片播放在线免费| 亚洲精华国产精华液的使用体验| 一级片'在线观看视频| tube8黄色片| 亚洲国产精品专区欧美| 女人精品久久久久毛片| 亚洲av免费高清在线观看| 秋霞伦理黄片| av一本久久久久| 宅男免费午夜| 夫妻午夜视频| 亚洲第一av免费看| av片东京热男人的天堂| 亚洲精品中文字幕在线视频| 国产永久视频网站| 国产精品久久久久成人av| 国产免费现黄频在线看| 黄色毛片三级朝国网站| 日本黄色日本黄色录像| 欧美性感艳星| 久久亚洲国产成人精品v| 成人影院久久| 一二三四在线观看免费中文在 | 99香蕉大伊视频| 黄色 视频免费看| 欧美日韩视频精品一区| 妹子高潮喷水视频| 成人手机av| 婷婷色av中文字幕| 国产精品成人在线| 青春草国产在线视频| 亚洲精品日韩在线中文字幕| 日韩成人av中文字幕在线观看| 亚洲少妇的诱惑av| 日韩av免费高清视频| 亚洲成av片中文字幕在线观看 | 下体分泌物呈黄色| 免费观看a级毛片全部| 国产高清三级在线| 天堂8中文在线网| 99国产综合亚洲精品| 免费黄色在线免费观看| 日韩制服骚丝袜av| 久久久欧美国产精品| 亚洲成人手机| 亚洲精品自拍成人| 亚洲图色成人| 国产片内射在线| 久久精品人人爽人人爽视色| av女优亚洲男人天堂| 五月伊人婷婷丁香| 免费人成在线观看视频色| 久久99热这里只频精品6学生| 亚洲第一区二区三区不卡| 丰满饥渴人妻一区二区三| 亚洲av男天堂| 亚洲国产精品999| 亚洲婷婷狠狠爱综合网| 婷婷色麻豆天堂久久| 亚洲av在线观看美女高潮| 久久久精品94久久精品| 久久久国产精品麻豆| 欧美精品亚洲一区二区| 最新中文字幕久久久久| 精品一品国产午夜福利视频| 在线观看免费高清a一片| 永久免费av网站大全| 国产成人欧美| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻人人爽人人添夜夜欢视频| 午夜福利影视在线免费观看| 两性夫妻黄色片 | 精品一区二区免费观看| 国产不卡av网站在线观看| 久久精品国产综合久久久 | 老熟女久久久| 成人影院久久| 在线天堂最新版资源| 国产亚洲一区二区精品| 黑人高潮一二区| 高清视频免费观看一区二区| 亚洲av免费高清在线观看| 亚洲天堂av无毛| 美女中出高潮动态图| 日日摸夜夜添夜夜爱| 国产精品 国内视频| 亚洲图色成人| 成人漫画全彩无遮挡| 肉色欧美久久久久久久蜜桃| 国产一级毛片在线| 国产男女内射视频| 亚洲av电影在线进入| 国产乱来视频区| 免费不卡的大黄色大毛片视频在线观看| 色5月婷婷丁香| 欧美 亚洲 国产 日韩一| 香蕉国产在线看| 国产xxxxx性猛交| 中文字幕av电影在线播放| 国产精品 国内视频| 精品亚洲乱码少妇综合久久| 国产精品女同一区二区软件| 久久99一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产视频首页在线观看| 你懂的网址亚洲精品在线观看| 蜜桃国产av成人99| 一级a做视频免费观看| 大香蕉久久网| 少妇熟女欧美另类| 日韩三级伦理在线观看| 九草在线视频观看| 久久久久国产精品人妻一区二区| 女人精品久久久久毛片| 捣出白浆h1v1| a 毛片基地| 极品人妻少妇av视频| 亚洲性久久影院| 国产片特级美女逼逼视频| 一区二区av电影网| 久久av网站| 成年人午夜在线观看视频| 亚洲美女视频黄频| 国产毛片在线视频| 亚洲国产精品成人久久小说| 极品人妻少妇av视频| 欧美+日韩+精品| 久久久久精品性色| 亚洲国产精品专区欧美| 91精品三级在线观看| 黄色视频在线播放观看不卡| 99热网站在线观看| 十八禁网站网址无遮挡| 狂野欧美激情性bbbbbb| 亚洲成国产人片在线观看| 亚洲国产日韩一区二区| 97精品久久久久久久久久精品| 亚洲成国产人片在线观看| 国产在线免费精品| a级毛色黄片| 国产亚洲欧美精品永久| 色94色欧美一区二区| 激情视频va一区二区三区| 日韩三级伦理在线观看| videosex国产| 超碰97精品在线观看| 国产精品欧美亚洲77777| 九九在线视频观看精品| 久久久久久人人人人人| 国产高清国产精品国产三级| 满18在线观看网站| 黑人巨大精品欧美一区二区蜜桃 | 在线天堂最新版资源| 久久久久久人人人人人| 成人18禁高潮啪啪吃奶动态图| 国产成人精品无人区| av黄色大香蕉| 日日爽夜夜爽网站|