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

    模擬氮沉降對(duì)東北地區(qū)興安落葉松樹(shù)干呼吸的影響

    2015-02-07 06:51:16劉瑞鵬李興歡毛子軍董利利
    生態(tài)學(xué)報(bào) 2015年11期
    關(guān)鍵詞:中氮落葉松樹(shù)干

    孫 濤,劉瑞鵬, 李興歡, 毛子軍, 黃 唯, 董利利

    東北林業(yè)大學(xué)森林植物生態(tài)學(xué)教育部重點(diǎn)實(shí)驗(yàn)室, 哈爾濱 150040

    模擬氮沉降對(duì)東北地區(qū)興安落葉松樹(shù)干呼吸的影響

    孫 濤,劉瑞鵬, 李興歡, 毛子軍*, 黃 唯, 董利利

    東北林業(yè)大學(xué)森林植物生態(tài)學(xué)教育部重點(diǎn)實(shí)驗(yàn)室, 哈爾濱 150040

    2011年的5月至2012年的12月,對(duì)興安落葉松 (Larixgmelinii)人工林進(jìn)行了模擬氮沉降實(shí)驗(yàn),設(shè)置了4個(gè)氮沉降水平處理,分別為:對(duì)照(0 g N m-2a-1)、低氮(5 g N m-2a-1)、中氮(5 g N m-2a-1)和高氮(15 g N m-2a-1),每個(gè)處理設(shè)置3個(gè)重復(fù)樣地。所施氮肥為NH4NO3, 以溶液方式在生長(zhǎng)季內(nèi)每月噴施1次。采用紅外氣體分析法于2012年4月底—10月份,每隔15d左右測(cè)定1次樹(shù)干呼吸,共測(cè)定12次。結(jié)果表明:各個(gè)氮處理下的樹(shù)干呼吸速率基本與樹(shù)干溫度的變化一致,均呈單峰型季節(jié)變化模式,其中7月底最高、10月底最低。氮處理均提高了落葉松的樹(shù)干呼吸速率,且隨著氮輸入水平的增加,平均樹(shù)干呼吸速率呈現(xiàn)逐漸增強(qiáng)的趨勢(shì)。生長(zhǎng)季內(nèi),對(duì)照、低氮、中氮、和高氮處理?xiàng)l件下的樹(shù)干表面釋放的CO2通量分別為67.3、72.5、78.9 g C/m2和86.5 g C/m2。樹(shù)干溫度與樹(shù)干呼吸速率存在顯著的指數(shù)函數(shù)關(guān)系,而且溫度敏感性(Q10)隨著氮輸入的逐漸升高亦隨之增強(qiáng),對(duì)照、低氮、中氮、和高氮處理下的Q10值分別為1.67、1.80、2.01和2.54。另外,伴隨氮輸入的逐漸增加,樹(shù)干的氮含量也逐漸升高,樹(shù)干氮含量與樹(shù)干呼吸速率之間也具有一定的相關(guān)性,能夠解釋樹(shù)干呼吸變化的38.3%。

    樹(shù)干呼吸; 氮沉降; 樹(shù)干溫度; 季節(jié)動(dòng)態(tài); 氮

    氮元素是影響森林生態(tài)系統(tǒng)生產(chǎn)力的主要限制因子[1- 2]。但自工業(yè)革命以來(lái),由于人類活動(dòng)向大氣中排放了大量的氮化物,導(dǎo)致這些氮化物在大氣中不斷累積并向陸地和水域生態(tài)系統(tǒng)沉降[3]。預(yù)計(jì)到2050年,全球氮沉降量將比1995年高出1倍,增加至200Tg N/a[4]。目前我國(guó)已成為繼歐洲和北美之后的第三大氮沉降區(qū),氮沉降現(xiàn)象在局部地區(qū)非常嚴(yán)重,而且伴隨著人類活動(dòng)的進(jìn)一步加強(qiáng),我國(guó)氮沉降正處于持續(xù)增加的狀態(tài)[5- 6]。大氣氮沉降的不斷增加,將會(huì)對(duì)陸地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能產(chǎn)生不同程度的影響。

    樹(shù)木的自養(yǎng)呼吸可消耗自身葉片光合作用固定碳的50%以上[7],對(duì)維持森林生態(tài)系統(tǒng)碳平衡具有非常重要的作用。其中,樹(shù)干呼吸消耗的碳占樹(shù)木自養(yǎng)呼吸的14%—48%[8- 10],其等于或有時(shí)已高于葉片呼吸釋放的碳通量[11]。另外,樹(shù)干呼吸也是樹(shù)木生長(zhǎng)活力及碳代謝狀況的重要物候?qū)W指標(biāo)之一[12- 13]。樹(shù)干釋放的CO2的來(lái)源和去向十分復(fù)雜,受多種因素的制約,主要包括生物因素和非生物因素,前者如樹(shù)干的生長(zhǎng)速率、液流密度、非結(jié)構(gòu)性碳水化合物含量、氮含量、直徑和邊材厚度等[9, 12, 14- 15],后者如樹(shù)干的溫度、大氣CO2濃度、樹(shù)干內(nèi)的氧氣濃度、土壤濕度及養(yǎng)分等[11, 16- 21]。其中氮元素含量與樹(shù)干呼吸密切相關(guān),這是因?yàn)橹参锝M織中90%的氮以蛋白質(zhì)的形式存在[22],樹(shù)木的生命活動(dòng)大約消耗總維持呼吸的20%[23]。因此,大氣氮沉降的持續(xù)增加會(huì)對(duì)樹(shù)干呼吸產(chǎn)生重要影響,然而關(guān)于這方面的研究還很缺乏,且研究結(jié)果不一致。Maier報(bào)道施肥后火炬松的樹(shù)干呼吸速率顯著增加, 與之伴隨的是樹(shù)干生長(zhǎng)速率和氮含量都顯著增加,而且發(fā)現(xiàn)樹(shù)干生長(zhǎng)速率和氮含量可解釋樹(shù)干基礎(chǔ)呼吸季節(jié)動(dòng)態(tài)的75%[14]。而別的研究結(jié)果與之相反,發(fā)現(xiàn)施肥后在沒(méi)有引起樹(shù)干生長(zhǎng)加快的情況下, 樹(shù)干呼吸速率也沒(méi)有顯著變化[24- 25]。

    落葉松(Larixgmelinii)是我國(guó)東北林區(qū)重要的建材樹(shù)種,落葉松人工林已構(gòu)成了東北地區(qū)主要人工林類型。目前Wang等[13]、石新立等[18]、許飛等[20]和王秀偉等[26]已對(duì)興安落葉松樹(shù)干呼吸的日變化和季節(jié)變化、種內(nèi)差異、溫度敏感性、與液流速度及胸徑的關(guān)系進(jìn)行了系統(tǒng)研究,但尚未開(kāi)展關(guān)于氮素添加對(duì)落葉松樹(shù)干呼吸影響方面的研究。我國(guó)東北地區(qū)的氮沉降量雖然低于較為嚴(yán)重的中南部地區(qū),但近年來(lái)隨著化肥施用量的逐漸增加等人類活動(dòng)的影響,氮沉降量也在持續(xù)增長(zhǎng),可能會(huì)對(duì)森林生態(tài)系統(tǒng)碳平衡產(chǎn)生重要影響。我們通過(guò)對(duì)黑龍江帽兒山2011年大氣年氮濕沉降進(jìn)行監(jiān)測(cè),發(fā)現(xiàn)該地區(qū)氮濕沉降為19.16 kg N hm-2a-1,其結(jié)果要明顯高于劉世榮[27]在1992年在該區(qū)測(cè)定的氮濕沉降值(12.9 kg N ha-1)。本研究以東北典型造林樹(shù)種興安落葉松為研究對(duì)象,在興安落葉松人工林內(nèi)進(jìn)行了不同水平的氮沉降實(shí)驗(yàn),旨在分析氮沉降對(duì)落葉松樹(shù)干呼吸通量的影響,并探討其影響機(jī)理,為全球變化尤其在氮沉降持續(xù)增加的背景下,理解溫帶森林碳循環(huán)過(guò)程及其模型構(gòu)建提供基礎(chǔ)數(shù)據(jù)和科學(xué)依據(jù)。

    1 研究方法

    1.1 研究地概況

    本試驗(yàn)地設(shè)立在黑龍江省尚志市帽兒山林場(chǎng)老山試驗(yàn)站(45°20′N, 127°34′E),地處長(zhǎng)白山系張廣才嶺西北部小嶺余脈,平均海拔340 m。該地區(qū)氣候?qū)儆诖箨懶约撅L(fēng)氣候,年降雨量600—800 mm,雨量主要集中于7、8月份,年蒸發(fā)量1093.9 mm, 年相對(duì)濕度70%;年平均溫度2.7 ℃,全年大于或等于10 ℃的積溫為2526 ℃,年平均總?cè)照諘r(shí)數(shù)2471.3 h。早霜一般在8月末出現(xiàn),晚霜截止到5月末,無(wú)霜期約120—140 d。地帶性土壤為典型暗棕壤[13]。

    該研究選取40年生興安落葉松人工林為研究對(duì)象,主要喬木層樹(shù)種為落葉松,伴有少量的水曲柳、白樺、胡桃楸和紫椴等,土壤為典型暗棕色森林土壤。土壤表層(0—10 cm)的有機(jī)質(zhì)、全氮和磷分別為(13.5±0.83)%、(0.96±0.12)%和(0.15±0.01)%;土壤表層(0—10 cm)的pH值為(6.28±0.09),土壤容重為(0.85±0.06) g/cm3。

    1.2 樣地設(shè)置和處理

    2011年4月在興安落葉松人工林內(nèi)設(shè)立12個(gè)10 m×20 m 的標(biāo)準(zhǔn)樣地。樣方之間留有15m左右的緩沖帶以防止相互之間干擾。實(shí)驗(yàn)分為4個(gè)處理組,分別為對(duì)照處理(CK)、低氮處理(LN)、中氮處理(MN)和高氮處理(HN),每個(gè)處理組分成3個(gè)重復(fù),即同一林分同一處理由隨機(jī)選擇的3個(gè)樣方組成。在2011年5月至10月份,2012年5月份至10月份間,每月底噴灑硝酸銨(NH4NO3),每個(gè)樣方每次將所需的氮肥溶解在50 L自來(lái)水中,然后采用背式噴霧器人工有序的均勻噴灑在林地,對(duì)照樣方噴施等量的自來(lái)水。連續(xù)2a內(nèi)共噴灑12次氮肥,生長(zhǎng)季內(nèi),每月噴灑氮量保持一致,分別按0、50、100和150 kg N hm-2a-1濃度進(jìn)行施氮處理。

    1.3 樣木選擇

    為減少樹(shù)干呼吸的種內(nèi)變異,在每塊樣地內(nèi)嚴(yán)格按照以下標(biāo)準(zhǔn)篩選出3株落葉松作為測(cè)定樣木:(1)樹(shù)干胸徑在(20 ± 1.5) cm范圍之內(nèi);(2)樹(shù)高和冠幅基本一致;(3)樹(shù)干通直,生長(zhǎng)狀態(tài)正常;(4)土壤條件排水較為良好。

    1.4 樹(shù)干呼吸和溫度的測(cè)定

    于2012年4月底至10月底測(cè)定樹(shù)干呼吸,半個(gè)月左右測(cè)定一次,生長(zhǎng)季內(nèi)共測(cè)定12次。采用石新立等[18]的測(cè)定方法,測(cè)定前,在選取的樣木正北面距離地面1.3 m處安裝1個(gè)內(nèi)徑為10.2 cm,高5 cm的PVC呼吸環(huán)。安裝前,首先除去松脫的表層樹(shù)皮、苔蘚等雜物,避免傷及樹(shù)皮活組織,影響樹(shù)干呼吸測(cè)定;然后,按照要安裝的樹(shù)干表面弧度將呼吸環(huán)的一端打磨成形狀與樹(shù)干表面相吻合的弧形;利用100%中性硅酮膠將制打磨好的PVC環(huán)固定在樹(shù)干表面上,并保持PVC呼吸環(huán)在整個(gè)測(cè)定期間位置不變且不漏氣。

    利用 LI- 6400 便攜式 CO2分析系統(tǒng)(LI-COR Inc. USA)進(jìn)行樹(shù)干呼吸的測(cè)定。考慮到樹(shù)干呼吸具有明顯的日變化格局,為更準(zhǔn)確的比較樹(shù)干呼吸對(duì)不同水平氮沉降的響應(yīng),將每次測(cè)定限定在晴天的10:00—12:00時(shí)間段內(nèi)進(jìn)行。測(cè)定時(shí),直接在LI- 6400-09同化室固定好的PVC環(huán)上進(jìn)行測(cè)定,每次測(cè)定3個(gè)循環(huán)。同步,在PVC環(huán)上方5 cm處放置接觸式溫度記錄儀測(cè)定樹(shù)干表皮層下1 cm深處的溫度。

    1.5 樹(shù)干氮素含量的測(cè)定

    在生長(zhǎng)季較為旺盛的7月末,利用內(nèi)徑為5 mm的生長(zhǎng)錐在每株樣木PVC呼吸環(huán)的下方5 cm處鉆取邊材最外側(cè)2 cm的樹(shù)干樣品,鉆取之前首先需將樹(shù)干表層的樹(shù)皮輕輕除掉,并在取樣之后立即用直徑像吻合的樹(shù)枝堵住樹(shù)干傷口。取回的樹(shù)干樣品65℃恒溫條件下至少烘48 h。采用半微量凱氏定氮法進(jìn)行樹(shù)干氮含量的測(cè)定,即稱磨細(xì)烘干的樹(shù)干樣品(0.25 mm 篩)0.1 g,置于100 mL的消化管中,用濃硫酸和硝酸混合樣消煮后,使用凱氏定氮儀進(jìn)行蒸餾并用標(biāo)準(zhǔn)酸溶液滴定[28- 30]。

    1.6 數(shù)據(jù)分析

    由于樹(shù)干表面的不規(guī)則性導(dǎo)致所測(cè)定的樹(shù)干表面積和呼吸環(huán)的體積可能會(huì)不一致,因此所測(cè)得的樹(shù)干呼吸速率(μmol CO2m-2s-1)需要進(jìn)行校正[17- 18, 31],校正采用以下公式:

    Rs=RM×FV×FA

    (1)

    其中:

    FA=S/S0

    (2)

    (3)

    式中,RS為校正后的樹(shù)干呼吸(μmol CO2m-2s-1),RM為校正前LI- 6400系統(tǒng)測(cè)定的樹(shù)干呼吸,S為系統(tǒng)默認(rèn)的同化室面積(81.7 cm2),S0為實(shí)際測(cè)定的樹(shù)干表面積(cm2),V1為同化室容積(991 cm3),V2為呼吸環(huán)的實(shí)際體積(cm3);D為呼吸環(huán)直徑(10.2 cm),H為呼吸環(huán)高度(5 cm)。

    采用公式(4)擬合樹(shù)干呼吸速率與樹(shù)干溫度間的關(guān)系[32]:

    Rs=R0eβT

    (4)

    溫度系數(shù)Q10值的計(jì)算公式如下:

    Q10= e10β

    (5)

    式中,RS為樹(shù)干呼吸速率(μmol m-2s-1);T為樹(shù)干溫度(℃);R0和β為常數(shù)。

    數(shù)據(jù)分析采用SPSS13.0完成。利用LSD多重檢驗(yàn)法檢驗(yàn)不同處理間的差異顯著性。RS生長(zhǎng)季內(nèi)的累積碳通量是基于實(shí)測(cè)數(shù)據(jù)建立的RS統(tǒng)計(jì)模型與連續(xù)樹(shù)干溫度監(jiān)測(cè)數(shù)據(jù)相結(jié)合的方法進(jìn)行估測(cè)。統(tǒng)計(jì)模型所需要的溫度參數(shù)來(lái)自于在研究樣地安裝的溫度自動(dòng)記錄儀Thermo Recorder Tr- 71U (T&D Corp.Japan),每隔30 min記錄1次樹(shù)干表皮層下1 cm深處的溫度。生長(zhǎng)季期間的碳通量根據(jù)上述模型以天為步長(zhǎng)累加而得。

    2 結(jié)果

    2.1 樹(shù)干呼吸的季節(jié)動(dòng)態(tài)

    不同氮處理?xiàng)l件下的樹(shù)干呼吸速率均存在明顯的季節(jié)變化,整個(gè)生長(zhǎng)季內(nèi)的呼吸速率呈單峰曲線,與樹(shù)干溫度的變化格局基本一致(圖1)。不同氮處理?xiàng)l件下的樹(shù)干溫度差異不限制(P> 0.05)。樹(shù)干呼吸速率從春季4月末開(kāi)始,隨著樹(shù)干溫度的升高也逐漸上升,在7月底均到達(dá)峰值,8月份以后隨著樹(shù)干溫度的下降,呼吸速率也開(kāi)始逐漸降低,10末達(dá)到最低值。生長(zhǎng)季內(nèi),對(duì)照、低氮、中氮、和高氮處理下的樹(shù)干呼吸速率最小值分別為1.33、1.60、1.85 μmol m-2s-1和1.73 μmol m-2s-1,最高值分別為4.71、5.01、5.38 μmol m-2s-1和5.96 μmol m-2s-1。相應(yīng)地,樹(shù)干呼吸的離散程度也呈現(xiàn)出類似的季節(jié)變化趨勢(shì),7月份的標(biāo)準(zhǔn)誤差是10月份的3.59倍。

    圖1 興安落葉松樹(shù)干溫度與呼吸速率的季節(jié)變化特征

    其中,CK、LN、MN和HN分別表示對(duì)照處理(0 kg N hm-2a-1)、低氮處理(50 kg N hm-2a-1)、中氮處理(100 kg N hm-2a-1)和高氮處理(150 kg N hm-2a-1)。

    2.2 不同氮處理對(duì)樹(shù)干呼吸的影響

    圖2 生長(zhǎng)季內(nèi)不同氮處理?xiàng)l件下的樹(shù)干表面CO2通量比較

    經(jīng)過(guò)連續(xù)2個(gè)生長(zhǎng)季的施氮處理,均顯著提高了落葉松的樹(shù)干呼吸速率(圖1b)。生長(zhǎng)季內(nèi),對(duì)照、低氮、中氮、和高氮處理?xiàng)l件下的平均樹(shù)干呼吸速率分別為(2.95±0.57)、(3.16±0.72)、(3.45±0.93) μmol m-2s-1和(3.82±0.82) μmol m-2s-1。而且隨氮處理水平增加, 生長(zhǎng)季內(nèi)樹(shù)干表面釋放的CO2通量也呈現(xiàn)逐漸上升的趨勢(shì),在高氮處理?xiàng)l件下達(dá)到峰值,為86.5g C/m2, 比對(duì)照高出了將近30%(圖2)。高氮處理下的樹(shù)干表面CO2通量與對(duì)照的差異達(dá)到顯著水平(P< 0.05),但與低氮和中氮處理下的差異不顯著。中氮處理?xiàng)l件下的碳通量也與對(duì)照的差異達(dá)到顯著水平(P< 0.05),但與其它處理?xiàng)l件下的差異均不顯著。低氮處理下的碳通量與各處理間的差異均未達(dá)到顯著水平。

    2.3 不同氮處理對(duì)樹(shù)干呼吸Q10值的影響

    對(duì)不同氮處理?xiàng)l件下的樹(shù)干呼吸速率與樹(shù)干溫度進(jìn)行的相關(guān)性分析表面,呼吸速率均與樹(shù)干溫度有顯著性的正相關(guān),并且樹(shù)干1 cm深處的溫度可以解釋樹(shù)干呼吸速率64%—82%的變化(表1)。

    研究發(fā)現(xiàn),溫度敏感性(Q10)對(duì)不同濃度的氮沉降的響應(yīng)也不同,隨著氮處理濃度增加,Q10值逐漸變大(表1)。對(duì)照、低氮、中氮、和高氮處理下的Q10值分別為1.67、1.80、2.01和2.54,高氮處理?xiàng)l件下的Q10與對(duì)照、低氮的差異達(dá)到顯著水平(P< 0.05),但與中氮處理下的差異不顯著。中氮處理?xiàng)l件下的Q10也與對(duì)照、低氮的差異達(dá)到顯著水平(P< 0.05)。低氮處理?xiàng)l件下的Q10與對(duì)照組間的差異未達(dá)到顯著水平,與其它各組間的差異均達(dá)到顯著水平(P< 0.05)。

    表1 生長(zhǎng)季內(nèi)不同氮處理?xiàng)l件下樹(shù)干呼吸速率與樹(shù)干溫度間的關(guān)系

    其中,CK、LN、MN和HN分別表示對(duì)照處理(0 kg N hm-2a-1)、低氮處理(50 kg N hm-2a-1)、中氮處理(100 kg N hm-2a-1)和高氮處理(150 kg N hm-2a-1)。

    2.4 樹(shù)干呼吸速率與樹(shù)干氮元素含量的關(guān)系

    圖3 樹(shù)干呼吸速率與樹(shù)干氮濃度之間的關(guān)系

    隨著氮處理濃度提高,樹(shù)干氮含量也隨之上升。其中對(duì)照、低氮、中氮、和高氮處理下氮元素的濃度分別為1.00、1.32、2.02和2.35 g/kg。高氮和中氮處理?xiàng)l件下的樹(shù)干氮含量與對(duì)照、低氮的差異均達(dá)到顯著水平(P< 0.05),但兩者之間的差異不顯著。低氮處理?xiàng)l件下的樹(shù)干氮含量也與對(duì)照組間的差異達(dá)到顯著水平(P< 0.05)。對(duì)樹(shù)干呼吸速率和樹(shù)干氮濃度的關(guān)系進(jìn)行線性回歸分析,結(jié)果表明樹(shù)干氮含量與樹(shù)干呼吸速率之間的相關(guān)性較弱,樹(shù)干氮濃度僅解釋了生長(zhǎng)季內(nèi)樹(shù)干呼吸速率38.3%的變化。

    3 討論與結(jié)論

    本研究結(jié)果表明,測(cè)定期間樹(shù)干呼吸速率表現(xiàn)為明顯的單峰型季節(jié)變化模式,生長(zhǎng)季的樹(shù)干呼吸顯著高于非生長(zhǎng)季。而且不同氮濃度處理均不同程度的促進(jìn)了落葉松樹(shù)干呼吸速率,其中以高氮濃度處理?xiàng)l件下樹(shù)干表面釋放的CO2通量最大,這與以往的研究結(jié)果相一致[14, 22, 24- 25]。樹(shù)干呼吸主要由提供能量合成新組織的生長(zhǎng)呼吸和維持活細(xì)胞正常生命活動(dòng)的維持呼吸兩部分組成[14]。本研究中氮沉降促進(jìn)樹(shù)干呼吸的可能原因主要有兩個(gè):一是氮元素增加促進(jìn)了植物組織的生長(zhǎng),從而引起植物呼吸尤其是生長(zhǎng)呼吸的增加[33]。目前大多數(shù)研究認(rèn)為溫帶森林生長(zhǎng)主要是氮元素的限制[34],大氣氮沉降通過(guò)提高土壤的有效氮水平,進(jìn)而促進(jìn)了植物的生長(zhǎng)[35],而植物的生長(zhǎng)速率與樹(shù)干的生長(zhǎng)呼吸存在密切的正相關(guān)關(guān)系[14, 24]。這很好的解釋了本研究中生長(zhǎng)季旺盛季節(jié),氮處理促進(jìn)了樹(shù)干呼吸速率的現(xiàn)象;二是外界氮輸入導(dǎo)致了植物細(xì)胞蛋白質(zhì)含量的增加。植物組織中90%的氮參與了蛋白質(zhì)的合成,外界氮輸入之后會(huì)使植物組織蛋白質(zhì)合成數(shù)量增加,而這些蛋白質(zhì)的修復(fù)和更新需要大量的能量支持,進(jìn)而促使樹(shù)干呼吸尤其是維持呼吸的上升[22]。這與本研究中生長(zhǎng)季末期,氮處理仍提高了樹(shù)干呼吸速率這一結(jié)果相吻合。

    本文研究還首次比較了不同氮處理?xiàng)l件下樹(shù)干呼吸的溫度敏感性(Q10)指數(shù),結(jié)果表明隨著氮輸入增加,Q10也隨之增加,這可能是由于樹(shù)干呼吸中維持呼吸和生長(zhǎng)呼吸對(duì)氮沉降的響應(yīng)模式不同引起的[36]。樹(shù)干組織氮含量對(duì)呼吸的影響(解釋38.3%的變化)可能也會(huì)間接影響樹(shù)干呼吸對(duì)溫度的敏感性。盡管不同氮處理的樹(shù)干呼吸對(duì)溫度的響應(yīng)程度不一,但總體上兩者均存在極顯著的指數(shù)函數(shù)關(guān)系(P< 0.01)。樹(shù)干溫度是影響樹(shù)干呼吸的重要環(huán)境因子,一方面是由于在一定范圍內(nèi)溫度升高增強(qiáng)了參與樹(shù)干呼吸作用的酶活性,促進(jìn)了樹(shù)干呼吸速率[24];另一方面,主要是溫度升高促進(jìn)了植物組織細(xì)胞呼吸、降低了CO2在水中的溶解度、提高了CO2的擴(kuò)散系數(shù),從而導(dǎo)致樹(shù)干表面釋放更多的CO2通量[37]。另外,本研究基于季節(jié)尺度計(jì)算而得到的樹(shù)干呼吸溫度敏感性與以往溫帶森林研究的結(jié)果范圍相符合[16, 19, 25, 37],但要明顯高于石新立等[18]在該地區(qū)測(cè)定的興安落葉松樹(shù)干呼吸溫度敏感性結(jié)果,盡管采用的研究方法相一致。導(dǎo)致不同研究中Q10變異較大的原因可能由于本研究中樹(shù)干呼吸的測(cè)定期間為4—10月份,測(cè)定間隔為15 d左右,而石新立等[18]的測(cè)定期間為5—9月份,測(cè)定時(shí)間尺度不同,直接會(huì)造成結(jié)果的很大差異;另外,本研究所選取樣木的胸徑在20 cm左右,明顯小于石新立等[18]測(cè)定樣木的胸徑(波動(dòng)在15.5—39.3 cm,平均為27.6 cm)。

    隨著氮輸入增加,樹(shù)干的氮含量也逐漸升高,這與以往的研究結(jié)果相一致[14, 22, 24]。例如Maier對(duì)火炬松林連續(xù)3a施肥后,發(fā)現(xiàn)氮輸入樣地的樹(shù)干氮含量比對(duì)照高出了1.5倍左右[14]。這可能是因?yàn)槊吭聡姙⒌腘H4NO3為可溶性氮,該部分氮很容易被植物組織吸收利用,進(jìn)而導(dǎo)致樹(shù)干氮含量升高。但是,目前很多研究,即使都采用同一方法,其結(jié)果仍然表明樹(shù)干氮含量與樹(shù)干呼吸之間的關(guān)系存在較大的差異。例如,本文與Lavigne等[38]很多研究的結(jié)果都表明樹(shù)干氮含量與樹(shù)干呼吸速率之間的相關(guān)較弱,而Vose等[9]的研究表明樹(shù)干氮含量與樹(shù)干呼吸速率之間存在顯著的正相關(guān)關(guān)系。導(dǎo)致樹(shù)干氮含量與樹(shù)干呼吸之間相關(guān)較弱的原因之一是由于樹(shù)干呼吸釋放的CO2的來(lái)源和去向十分復(fù)雜,實(shí)際上采用紅外線氣體分析儀LI- 6400同化室所測(cè)得的變量只是實(shí)際樹(shù)干呼吸產(chǎn)生的CO2的一部分,而剩余的部分CO2可能還溶解于樹(shù)干液流中,或被帶至冠層上方,或被樹(shù)皮光合作用所利用[15, 37];其二是進(jìn)入植物組織的氮素并沒(méi)有向蛋白質(zhì)分配,而是以游離氨基酸的形態(tài)貯存著,這部分氮將可能不會(huì)對(duì)呼吸速率造成影響[9]。

    通過(guò)模擬氮沉降對(duì)興安落葉松樹(shù)干呼吸的影響,發(fā)現(xiàn)施肥2a后,隨著氮輸入的增多,樹(shù)干呼吸速率、樹(shù)干氮含量和溫度敏感性都呈現(xiàn)逐漸增加的趨勢(shì)??梢?jiàn),氮沉降是影響溫帶森林樹(shù)干呼吸速率非常重要的因素之一。

    [1] LeBauer D S, Treseder K K. Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology, 2008, 89(2): 371- 379.

    [2] 李德軍, 莫江明, 方運(yùn)霆, 彭少麟, Gundersen P. 氮沉降對(duì)森林植物的影響. 生態(tài)學(xué)報(bào), 2003, 23(9): 1891- 1900.

    [3] Vitousek P M, Aber J D, Howarth R W, Likens G E, Matson P A, Schindler D W, Schlesinger W H, Tilman G D. Human alteration of the global nitrogen cycle: sources and consequences. Ecological Applications, 1997, 7(3): 737- 750.

    [4] Galloway J N, Townsend A R, Erisman J W, Bekunda M, Cai Z C, Freney J R, Martinelli L A, Seitzinger S P, Sutton M A. Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science, 2008, 320(5878): 889- 892.

    [5] Galloway J N, Dentener F J, Capone D G, Boyer E W, Howarth R W, Seitzinger S P, Asner G P, Cleveland C C, Green P A, Holland E A, Karl D M, Michaels A F, Porter J H, Townsend A R, V?r?smarty C J. Nitrogen cycles: past, present, and future. Biogeochemistry, 2004, 70(2): 153- 226.

    [6] Liu X J, Zhang Y, Han W X, Tang A H, Shen J L, Cui Z L, Vitousek P, Erisman J W, Goulding K, Christie P, Fangmeier A, Zhang F S. Enhanced nitrogen deposition over China. Nature, 2013, 494(7438): 459- 462.

    [7] Ryan M G. Effects of climate change on plant respiration. Ecological Applications, 1991, 1(2): 157- 167.

    [8] 曾小平, 彭少麟, 趙平.廣東南亞熱帶馬占相思林呼吸量的測(cè)定. 植物生態(tài)學(xué)報(bào), 2000, 24(4): 420- 424.

    [9] Vose J M, Ryan M G. Seasonal respiration of foliage, fine roots, and woody tissues in relation to growth, tissue N, and photosynthesis. Global Change Biology, 2002, 8(2): 182- 193.

    [10] 肖復(fù)明, 汪思龍, 杜天真, 陳龍池, 于小軍. 湖南會(huì)同林區(qū)杉木人工林呼吸量測(cè)定. 生態(tài)學(xué)報(bào), 2005, 25(10): 2514- 2519.

    [11] Edwards N T, Hanson P J. Stem respiration in a closed-canopy upland oak forest. Tree Physiology, 1996, 16(4): 433- 439.

    [12] Edwards N T, Tschaplinski T J, Norby R J. Stem respiration increases in CO2-enriched sweetgum trees. New Phytologist, 2002, 155(2): 239- 248.

    [13] Wang W J, Yang F J, Zu Y G, Wang H M, Kentaro T, Kaichiro S, Takayoshi K. Stem respiration of a larch (Larixgmelinii) plantation in Northeast China. Acta Botanica Sinica, 2003, 45(12): 1387- 1397.

    [14] Maier C A. Stem growth and respiration in loblolly pine plantations differing in soil resource availability. Tree Physiology, 2001, 21(16): 1183- 1193.

    [15] Teskey R O, Saveyn A, Steppe K, McGuire M A. Origin, fate and significance of CO2in tree stems. New Phytologist, 2008, 177(1): 17- 32.

    [16] 王淼, 武耀祥, 武靜蓮. 長(zhǎng)白山紅松針闊葉混交林主要樹(shù)種樹(shù)干呼吸速率. 應(yīng)用生態(tài)學(xué)報(bào), 2008, 19(5): 956- 960.

    [17] 王淼, 姬蘭柱, 李秋榮, 肖冬梅, 劉海良. 長(zhǎng)白山地區(qū)紅松樹(shù)干呼吸的研究. 應(yīng)用生態(tài)學(xué)報(bào), 2005, 16(1): 7- 13.

    [18] 石新立, 王傳寬, 許飛, 王興昌. 四個(gè)溫帶樹(shù)種樹(shù)干呼吸的時(shí)間動(dòng)態(tài)及其影響因子. 生態(tài)學(xué)報(bào), 2010, 30(15): 3994- 4003.

    [19] 王文杰, 王慧梅, 祖元?jiǎng)? 李雪瑩, 小池孝良. 林木非同化器官與土壤呼吸的溫度系數(shù)Q10值的特征分析. 植物生態(tài)學(xué)報(bào), 2005, 29(4): 680- 691.

    [20] 許飛, 王傳寬, 王興昌. 東北東部14個(gè)溫帶樹(shù)種樹(shù)干呼吸的種內(nèi)種間變異. 生態(tài)學(xué)報(bào), 2011, 31(13): 358- 3589.

    [21] Etzold S, Zweifel R, Ruehr N K, Eugster W, Buchmann N. Long-term stem CO2concentration measurements in Norway spruce in relation to biotic and abiotic factors. New Phytologist, 2013, 197(4): 1173- 1184.

    [22] Ryan M G. Foliar maintenance respiration of subalpine and boreal trees and shrubs in relation to nitrogen content. Plant, Cell and Environment, 1995, 18(4): 765- 772.

    [23] Wieser G, Bahn M. Seasonal and spatial variation of woody tissue respiration in aPinuscembratree at the alpine timberline in the central Austrian Alps. Trees, 2004, 18(5): 576- 580.

    [24] Stockfors J, Linder S. Effect of nitrogen on the seasonal course of growth and maintenance respiration in stems of Norway spruce trees. Tree Physiology, 1998, 18(3): 155- 166.

    [25] Ryan M G, Hubbard R M, Pongracic S, Raison R J, McMurtrie R E. Foliage, fine-root, woody-tissue and stand respiration in Pinus radiata in relation to nitrogen status. Tree Physiology, 1996, 16(3): 333- 343.

    [26] 王秀偉, 毛子軍, 孫濤, 吳海軍. 春、秋季節(jié)樹(shù)干溫度和液流速度對(duì)東北3樹(shù)種樹(shù)干表面CO2釋放通量的影響. 生態(tài)學(xué)報(bào), 2011, 31(12): 3358- 3367.

    [27] 劉世榮. 興安落葉松人工林生態(tài)系統(tǒng)營(yíng)養(yǎng)元素生物地球化學(xué)循環(huán)特征. 生態(tài)學(xué)雜志1992, 11(5): 1- 6.

    [28] Sun T, Mao Z J. Functional relationships between morphology and respiration of fine roots in two Chinese temperate tree species. Plant and Soil, 2011, 346(1/2): 375- 384.

    [29] Sun T, Mao Z J, Dong L L, Hou L L, Song Y, Wang X W. Further evidence for slow decomposition of very fine roots using two methods: litterbags and intact cores. Plant and Soil, 2013, 366(1/2): 633- 646.

    [30] Sun T, Mao Z J, Han Y Y. Slow decomposition of very fine roots and some factors controlling the process: a 4-year experiment in four temperate tree species. Plant and Soil, 2013, 372(1/2): 445- 458.

    [31] Zha T, Kellomǎki S, Wang K Y, Ryyppǒ A, Niinistǒ S. Seasonal and annual stem respiration of Scots pine trees under boreal conditions. Annals of Botany, 2004, 94(6): 889- 896.

    [32] Davidson E A, Belk E, Boone R D. Soil water content and temperature as independent or confounded factors controlling soil respiration in a temperate mixed hardwood forest. Global Change Biology, 1998, 4(2): 217- 227.

    [33] Vitousek P M, Aber J D, Howarth R W, Likens G E, Matson P A, Schindler D W, Schlesinger W H, Tilman D G, Human alteration of the global nitrogen cycle: sources and consequences. Ecological Applications, 1997, 7(3): 737- 750.

    [34] Vitousek P M, Howarth R W. Nitrogen limitation on land and sea: how can it occur? Biogeochemistry,1991, 13(2): 87- 115.

    [35] Frink C R, Waggoner P E, Ausubel J H. Nitrogen fertilizer: retrospect and prospect. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96(4): 1175- 1180.

    [36] Bredu S A, Yokota T, Hagihara A. Temperature effect on maintenance and growth respiration coefficients of young, field- grown Hinoki cypress (Chamaecyparisobtusa). Ecological Research, 1997, 12(3): 357- 362.

    [37] McGuire M A, Cerasoli S, Teskey R O. CO2fluxes and respiration of branch segments of sycamore (PlatanusoccidentalisL.) examined at different sap velocities, branch diameters and temperatures. Journal of Experimental Botany, 2007, 58(3): 2159- 2168.

    [38] Lavigne M B, Ryan M G. Growth and maintenance respiration rate of aspen, black spruce and jack pine stems at northern and southern BORE AS sites. Tree Physiology, 1997, 17(8/9): 543- 551.

    Effects of simulated nitrogen deposition on stem respiration ofLarixgmeliniiRupr in Northeastern China

    SUN Tao, LIU Ruipeng, LI Xinghuan, MAO Zijun*, HUANG Wei, DONG Lili

    KeyLaboratoryofForestPlantEcology,MinistryofEducation,NortheastForestryUniversity,Harbin150040,China

    Anthropogenic activities such as fossil fuel emissions and ammonium from production and application of fertilizers have elevated atmospheric nitrogen (N) deposition by three- to five-fold over the past century. Nitrogen deposition also has been accounted as one major component of global change and is predicted to increase in the future. Nitrogen additions in the atmosphere and soil have been demonstrated to alter regional and global environments and affect forest ecosystem carbon cycle, which is projected to alter the trajectory of future climate change through climate-carbon cycle feedback. Concern about the ecological effects of elevated nitrogen deposition on terrestrial ecosystems is currently increasing. The stem surface carbon flux is one major component of tree autotrophic respiration, and may represent tree growth activities and metabolism. Despite the importance of stem respiration, few studies have been conducted to examine the effects of simulated nitrogen deposition on stem respiration in forest ecosystem. In this study, response of simulated nitrogen deposition on stem respiration in a 40-year-oldLarixgmeliniiplantation in Northeastern China was studied. The objective was to test the hypothesis that N addition would increase stem respiration in N limited temperate forests. Nitrogen deposition experiments were initiated in April 2011. Four nitrogen deposition treatments (in three replicates) were established: Control (without N added), Low-N (5 g N m-2a-1), Medium-N (10 g N m-2a-1), High-N (15 g N m-2a-1). Twelve plots, each with 10 m×20 m dimensions were set up, surrounded by an approximately 15 m wide buffer strip. All plots and treatments were installed randomly. Ammonium nitrate solution was sprayed monthly as 6 equal applications from May 2011 to October 2011, and May 2012 to October 2012. During each application, fertilizer was weighed, mixed with 50 L of water, and applied to each plot below the canopy using a backpack sprayer. Two passes were made across each plot in order to ensure an even distribution of fertilizer. The Control plot received 50 L water without N. A polyvinyl chloride collar (10.2 cm inside diameter×5 cm height) was cut and polished to fit the stem shape of each sampled tree,and set up on the northward side with silicon adhesive at breast height. During the growing season between April 2012 and October 2012, stem respiration was measured once a week from 10:00 to 12:00 with a Li- 6400 portable CO2infrared gas analyzer equipped with a Li- 6400-09 chamber. Stem temperature at 1 cm depth beneath the bark was simultaneously measured with a digital thermometer. The measured stem respiration was calibrated with the polyvinyl chloride collar volume and stem surface area of each sampled tree. Our results showed that stem respiration exhibited a strong seasonal pattern, with the highest rates found in July and the lowest rates in October. The estimated stem respiration during the measuring period varied from 1.33 to 4.71 μmol m-2s-1in the Control treatment, 1.60 to 5.01 μmol m-2s-1in the Low-N treatment, 1.85 to 5.38 μmol m-2s-1in the Medium-N and 1.73 to 5.96 μmol m-2s-1in the High-N treatment, respectively. Stem respiration rates showed a significant positive exponential relationship with stem temperature within all treatments. We fitted site-based models and used continuous measurements of stem temperature to estimate cumulative stem respiration for the growing season of 2012. Cumulative stem respiration in the growing season was estimated to 67.3 g C/m2, 72.5 g C/m2, 78.9 g C/m2and 86.5 g C/m2in the Control, Low-N, Medium-N and High-N treatment, respectively. The measured apparent temperature sensitivity of stem respiration (Q10) was the lowest in the Control stand, and increased from the Control treatment to High-N treatment. TheQ10forLarixgmeliniiin the Control, Low-N, Medium-N and High-N treatment during the measuring period was 1.67, 1.80, 2.01 and 2.54, respectively. Stem nitrogen concentration also increased from Control treatment to High-N treatment, and peaked at High-N treatment. The estimated for stem nitrogen concentrationLarixgmeliniiin the Control, Low-N, Medium-N and High-N treatment during the measuring period was 1.00,1.32,2.02 and 2.35 g/kg, respectively. However, stem respiration was weakly correlated with stem nitrogen concentration. Stem nitrogen concentration explained 38.3% of variations during the measuring period of 2012. Our results suggest that simulated nitrogen deposition increased stem respiration inLarixgmeliniiplantation in Northeastern China, but it may vary depending on the rate of nitrogen deposition.

    stem respiration; nitrogen deposition; stem temperature; seasonal dynamics; nitrogen

    國(guó)家自然科學(xué)基金項(xiàng)目(31270494, 31070350); 國(guó)家“973”重點(diǎn)基礎(chǔ)研究項(xiàng)目發(fā)展計(jì)劃項(xiàng)目(2010CB951301)

    2013- 08- 04;

    2014- 06- 12

    10.5846/stxb201308042017

    *通訊作者Corresponding author.E-mail: zijunm@ nefu.edu.cn

    孫濤,劉瑞鵬, 李興歡, 毛子軍, 黃唯, 董利利.模擬氮沉降對(duì)東北地區(qū)興安落葉松樹(shù)干呼吸的影響.生態(tài)學(xué)報(bào),2015,35(11):3684- 3691.

    Sun T, Liu R P, Li X H, Mao Z J, Huang W, Dong L L.Effects of simulated nitrogen deposition on stem respiration ofLarixgmeliniiRupr in Northeastern China.Acta Ecologica Sinica,2015,35(11):3684- 3691.

    猜你喜歡
    中氮落葉松樹(shù)干
    山西落葉松雜交良種逾10萬(wàn)畝
    為什么樹(shù)干不是方的?
    軍事文摘(2021年16期)2021-11-05 08:49:16
    長(zhǎng)白落葉松離體再生體系的建立
    神奇的落葉松提取物
    為什么要在樹(shù)干上刷白漿
    為什么要在樹(shù)干上刷白漿
    為什么要在樹(shù)干上刷一層白漿
    清水江水體和沉積物中氮、磷的分布及變化趨勢(shì)
    蒸餾滴定法測(cè)定高溫合金中氮的含量
    西藏濕地生態(tài)系統(tǒng)中氮循環(huán)微生物數(shù)量和多樣性研究
    西藏科技(2015年4期)2015-09-26 12:12:59
    国产精品乱码一区二三区的特点| 男人舔女人的私密视频| 无限看片的www在线观看| 老司机在亚洲福利影院| 国产亚洲av高清不卡| 欧美三级亚洲精品| av天堂在线播放| 女人高潮潮喷娇喘18禁视频| 啦啦啦免费观看视频1| 久久精品国产亚洲av香蕉五月| av欧美777| 天堂动漫精品| 俄罗斯特黄特色一大片| 又大又爽又粗| 色播亚洲综合网| tocl精华| 99久久国产精品久久久| 午夜福利成人在线免费观看| 正在播放国产对白刺激| 波多野结衣巨乳人妻| 一本大道久久a久久精品| 国产亚洲精品久久久久久毛片| 欧美在线黄色| 久久精品影院6| 久久性视频一级片| 特大巨黑吊av在线直播| 高清在线国产一区| 欧美激情久久久久久爽电影| 一二三四在线观看免费中文在| 久久婷婷成人综合色麻豆| av天堂在线播放| 国产在线观看jvid| 日本免费a在线| 久久精品成人免费网站| 99久久国产精品久久久| 国产亚洲精品综合一区在线观看 | 琪琪午夜伦伦电影理论片6080| 久久久精品欧美日韩精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美成人免费av一区二区三区| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区| 一本精品99久久精品77| 欧美性长视频在线观看| 日韩免费av在线播放| 欧美在线黄色| or卡值多少钱| 两性夫妻黄色片| 国产视频内射| 精品久久久久久久久久免费视频| 91九色精品人成在线观看| 国产精品1区2区在线观看.| 亚洲中文字幕日韩| 国内揄拍国产精品人妻在线| 国产精品一区二区三区四区久久| 精品免费久久久久久久清纯| 午夜精品在线福利| 欧美中文综合在线视频| 亚洲在线自拍视频| 亚洲专区国产一区二区| 日韩成人在线观看一区二区三区| 欧美午夜高清在线| 51午夜福利影视在线观看| 亚洲五月婷婷丁香| 校园春色视频在线观看| 国产成人精品久久二区二区免费| 国产熟女xx| 波多野结衣巨乳人妻| 欧美乱码精品一区二区三区| 欧美日本视频| 色哟哟哟哟哟哟| 村上凉子中文字幕在线| 亚洲精品国产精品久久久不卡| 18禁美女被吸乳视频| 成人手机av| 欧美黄色淫秽网站| 精品欧美一区二区三区在线| 又黄又爽又免费观看的视频| 嫩草影视91久久| 亚洲av成人av| 69av精品久久久久久| 亚洲人成电影免费在线| xxx96com| 精品不卡国产一区二区三区| 91字幕亚洲| 在线视频色国产色| 成在线人永久免费视频| 高清毛片免费观看视频网站| aaaaa片日本免费| 日本一本二区三区精品| 久久香蕉激情| 高清在线国产一区| 99久久久亚洲精品蜜臀av| 天堂√8在线中文| 色播亚洲综合网| 日本三级黄在线观看| 岛国在线观看网站| 国产成人系列免费观看| 久久久精品欧美日韩精品| a级毛片在线看网站| 国产欧美日韩一区二区精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区二区三区不卡视频| 日本黄大片高清| 色在线成人网| 国产精品 国内视频| 一进一出抽搐动态| 国产亚洲av嫩草精品影院| 国产v大片淫在线免费观看| 脱女人内裤的视频| 亚洲欧美精品综合一区二区三区| 欧美午夜高清在线| 一边摸一边抽搐一进一小说| 国产亚洲欧美98| 正在播放国产对白刺激| 欧美+亚洲+日韩+国产| 香蕉av资源在线| 中文字幕熟女人妻在线| 国产久久久一区二区三区| 欧美高清成人免费视频www| 成人特级黄色片久久久久久久| 亚洲美女视频黄频| 亚洲五月天丁香| 黄色片一级片一级黄色片| 伦理电影免费视频| 无限看片的www在线观看| av在线播放免费不卡| 亚洲国产欧美人成| 午夜精品久久久久久毛片777| 男插女下体视频免费在线播放| 久久久久国产一级毛片高清牌| 麻豆一二三区av精品| 在线观看免费视频日本深夜| 亚洲精品美女久久久久99蜜臀| 亚洲国产精品成人综合色| 桃色一区二区三区在线观看| 亚洲人成网站高清观看| 亚洲美女黄片视频| 亚洲精品中文字幕一二三四区| 国产成人影院久久av| 免费在线观看日本一区| 亚洲一区二区三区色噜噜| 国产精品香港三级国产av潘金莲| 日韩大码丰满熟妇| 久久人妻av系列| 777久久人妻少妇嫩草av网站| 国产精品,欧美在线| a级毛片在线看网站| 精品久久久久久成人av| 亚洲精品一卡2卡三卡4卡5卡| 在线永久观看黄色视频| 中国美女看黄片| 婷婷丁香在线五月| 国产精品一区二区三区四区免费观看 | 国产成人啪精品午夜网站| 亚洲成av人片在线播放无| 一区二区三区高清视频在线| 午夜免费观看网址| 两人在一起打扑克的视频| 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕在线视频| 国产欧美日韩一区二区三| 亚洲成人久久爱视频| 久久久久久久久免费视频了| 亚洲 欧美 日韩 在线 免费| 午夜a级毛片| 后天国语完整版免费观看| 制服人妻中文乱码| 全区人妻精品视频| 久久99热这里只有精品18| 久久久久久久精品吃奶| 欧美成人免费av一区二区三区| 在线十欧美十亚洲十日本专区| 村上凉子中文字幕在线| 欧美性猛交╳xxx乱大交人| 国产成人av教育| 国产av麻豆久久久久久久| 国产一区二区三区视频了| 两个人免费观看高清视频| 精品免费久久久久久久清纯| 老司机午夜十八禁免费视频| 欧美成人免费av一区二区三区| 老司机深夜福利视频在线观看| 久久久久久九九精品二区国产 | 国内精品久久久久精免费| 一级黄色大片毛片| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看 | 精品福利观看| 一级a爱片免费观看的视频| 国产精品爽爽va在线观看网站| 后天国语完整版免费观看| 久久午夜综合久久蜜桃| 韩国av一区二区三区四区| 午夜视频精品福利| 中文字幕人成人乱码亚洲影| 男女之事视频高清在线观看| 国产区一区二久久| 天天一区二区日本电影三级| 亚洲欧美日韩无卡精品| 99久久无色码亚洲精品果冻| 美女扒开内裤让男人捅视频| 国产精品国产高清国产av| 亚洲国产精品sss在线观看| 国产精品自产拍在线观看55亚洲| 老司机午夜十八禁免费视频| 午夜日韩欧美国产| 日韩欧美国产一区二区入口| 久久久久久大精品| 欧美色欧美亚洲另类二区| 亚洲av片天天在线观看| 亚洲乱码一区二区免费版| 日本免费a在线| 首页视频小说图片口味搜索| 精品一区二区三区视频在线观看免费| 天堂动漫精品| 99精品久久久久人妻精品| 国内精品久久久久精免费| 岛国在线观看网站| 男女做爰动态图高潮gif福利片| 欧美在线一区亚洲| 婷婷丁香在线五月| 久久香蕉精品热| 国产欧美日韩一区二区三| 91国产中文字幕| 精品久久久久久久久久久久久| 国产精品1区2区在线观看.| 国产免费av片在线观看野外av| 亚洲国产精品成人综合色| 欧美中文综合在线视频| 欧美成人午夜精品| 在线观看www视频免费| 五月玫瑰六月丁香| 很黄的视频免费| 香蕉久久夜色| 少妇人妻一区二区三区视频| 亚洲欧美一区二区三区黑人| 久久久久精品国产欧美久久久| 国产亚洲av高清不卡| 两个人免费观看高清视频| 亚洲av第一区精品v没综合| 国产成人欧美在线观看| 久热爱精品视频在线9| 国产在线观看jvid| 久久性视频一级片| 久久 成人 亚洲| 他把我摸到了高潮在线观看| a级毛片在线看网站| 中文在线观看免费www的网站 | 亚洲男人天堂网一区| 国产av又大| 日日摸夜夜添夜夜添小说| 日韩欧美国产一区二区入口| 九色成人免费人妻av| av福利片在线| 国产伦在线观看视频一区| 女人高潮潮喷娇喘18禁视频| 久久精品国产亚洲av高清一级| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩瑟瑟在线播放| 亚洲精品久久成人aⅴ小说| 一级片免费观看大全| 一本精品99久久精品77| 日本免费一区二区三区高清不卡| 久久久久久大精品| 又爽又黄无遮挡网站| 国产精品一区二区三区四区免费观看 | 国产免费av片在线观看野外av| 丝袜人妻中文字幕| 成人欧美大片| 97人妻精品一区二区三区麻豆| 一个人观看的视频www高清免费观看 | 国产三级在线视频| 亚洲国产欧美一区二区综合| 欧美中文日本在线观看视频| 国产精品一区二区精品视频观看| 99热只有精品国产| 国产精品爽爽va在线观看网站| 日韩欧美在线二视频| 亚洲精品中文字幕在线视频| 丝袜美腿诱惑在线| 色综合欧美亚洲国产小说| 国内精品一区二区在线观看| 国产主播在线观看一区二区| 国产精品久久视频播放| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品久久久久久毛片| 精品乱码久久久久久99久播| 免费在线观看亚洲国产| 亚洲一区中文字幕在线| 亚洲精品粉嫩美女一区| 黄色a级毛片大全视频| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 男男h啪啪无遮挡| 日韩欧美精品v在线| 中国美女看黄片| 日日爽夜夜爽网站| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 精品福利观看| 国产高清有码在线观看视频 | 黑人巨大精品欧美一区二区mp4| 日韩三级视频一区二区三区| 麻豆成人午夜福利视频| 动漫黄色视频在线观看| 日韩欧美免费精品| 欧美一级a爱片免费观看看 | av片东京热男人的天堂| 岛国在线观看网站| 欧美zozozo另类| 成年女人毛片免费观看观看9| 天天躁夜夜躁狠狠躁躁| 一夜夜www| 99国产精品99久久久久| 亚洲精品久久成人aⅴ小说| 99在线视频只有这里精品首页| 无限看片的www在线观看| 99热只有精品国产| 精品欧美一区二区三区在线| 此物有八面人人有两片| 日本 欧美在线| 日日爽夜夜爽网站| 啦啦啦免费观看视频1| 欧美黄色淫秽网站| 久久国产精品影院| 国产精品av久久久久免费| 色综合欧美亚洲国产小说| 999精品在线视频| 精品国产乱码久久久久久男人| 日本在线视频免费播放| 草草在线视频免费看| 国产熟女午夜一区二区三区| 亚洲一区中文字幕在线| 少妇裸体淫交视频免费看高清 | 欧美另类亚洲清纯唯美| 久久精品人妻少妇| 国产成人影院久久av| 亚洲成av人片在线播放无| 久久精品夜夜夜夜夜久久蜜豆 | 国产又黄又爽又无遮挡在线| 国产av在哪里看| www.www免费av| 九色国产91popny在线| 婷婷精品国产亚洲av| 一二三四社区在线视频社区8| 桃色一区二区三区在线观看| 国产精品一区二区精品视频观看| 宅男免费午夜| 亚洲人成伊人成综合网2020| 女人爽到高潮嗷嗷叫在线视频| 最近最新中文字幕大全电影3| 亚洲人成网站在线播放欧美日韩| 成在线人永久免费视频| 精品欧美一区二区三区在线| 99久久久亚洲精品蜜臀av| 听说在线观看完整版免费高清| 精品欧美一区二区三区在线| 精品电影一区二区在线| 亚洲人与动物交配视频| 国产三级中文精品| 久久久久久亚洲精品国产蜜桃av| 午夜影院日韩av| 嫁个100分男人电影在线观看| 久久久久久久久中文| 久久精品国产亚洲av香蕉五月| 久久久精品国产亚洲av高清涩受| 亚洲欧美精品综合一区二区三区| 视频区欧美日本亚洲| 亚洲精品国产一区二区精华液| 亚洲成人久久性| 51午夜福利影视在线观看| 91老司机精品| 老汉色∧v一级毛片| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 亚洲专区中文字幕在线| 最近最新免费中文字幕在线| 久久精品国产综合久久久| 国产av一区二区精品久久| 久久久久国内视频| 99久久国产精品久久久| 国产1区2区3区精品| 又紧又爽又黄一区二区| 两个人免费观看高清视频| 高清在线国产一区| av国产免费在线观看| 99国产精品一区二区蜜桃av| 757午夜福利合集在线观看| 91成年电影在线观看| 制服人妻中文乱码| 国产爱豆传媒在线观看 | 无限看片的www在线观看| 欧美午夜高清在线| 国产精品一区二区免费欧美| 美女午夜性视频免费| 欧美日韩一级在线毛片| 又黄又爽又免费观看的视频| 母亲3免费完整高清在线观看| 黑人操中国人逼视频| 国产三级中文精品| 51午夜福利影视在线观看| 国产精品一区二区三区四区免费观看 | 波多野结衣高清作品| 精品日产1卡2卡| 亚洲全国av大片| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 五月玫瑰六月丁香| 级片在线观看| 午夜影院日韩av| 好男人在线观看高清免费视频| 婷婷亚洲欧美| 琪琪午夜伦伦电影理论片6080| 国内久久婷婷六月综合欲色啪| 日韩欧美国产一区二区入口| 午夜精品在线福利| 亚洲av五月六月丁香网| 一个人免费在线观看的高清视频| 50天的宝宝边吃奶边哭怎么回事| 精华霜和精华液先用哪个| 制服诱惑二区| www国产在线视频色| a在线观看视频网站| 亚洲av五月六月丁香网| 在线观看美女被高潮喷水网站 | 午夜老司机福利片| 亚洲精品在线观看二区| 免费在线观看黄色视频的| 午夜两性在线视频| 制服丝袜大香蕉在线| 国内精品一区二区在线观看| 精品国产乱码久久久久久男人| 丁香六月欧美| 在线永久观看黄色视频| 校园春色视频在线观看| 97人妻精品一区二区三区麻豆| 老鸭窝网址在线观看| 熟妇人妻久久中文字幕3abv| 欧美性猛交黑人性爽| 精品一区二区三区av网在线观看| 人人妻人人澡欧美一区二区| 亚洲精品av麻豆狂野| 婷婷六月久久综合丁香| 国产精品一区二区免费欧美| 91成年电影在线观看| 免费在线观看视频国产中文字幕亚洲| 18美女黄网站色大片免费观看| 日本黄色视频三级网站网址| 免费在线观看黄色视频的| 2021天堂中文幕一二区在线观| 国产精品亚洲av一区麻豆| 欧美精品亚洲一区二区| av有码第一页| 国产片内射在线| 国产av在哪里看| 国产欧美日韩一区二区精品| 国产视频一区二区在线看| 欧美性猛交黑人性爽| 国产探花在线观看一区二区| 欧美黑人欧美精品刺激| 国产高清videossex| 精华霜和精华液先用哪个| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 久久九九热精品免费| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 90打野战视频偷拍视频| 午夜成年电影在线免费观看| 午夜福利18| 欧美色视频一区免费| 国产视频一区二区在线看| 老鸭窝网址在线观看| 欧美日韩瑟瑟在线播放| 亚洲欧美一区二区三区黑人| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 中文字幕高清在线视频| 亚洲精品粉嫩美女一区| 老司机深夜福利视频在线观看| 精品国产乱码久久久久久男人| 亚洲av片天天在线观看| 国产欧美日韩一区二区三| 久久 成人 亚洲| 久久久水蜜桃国产精品网| 久久热在线av| 少妇被粗大的猛进出69影院| 怎么达到女性高潮| 国产欧美日韩精品亚洲av| 国内久久婷婷六月综合欲色啪| 久久人人精品亚洲av| 在线看三级毛片| 国内少妇人妻偷人精品xxx网站 | 18禁裸乳无遮挡免费网站照片| www日本在线高清视频| 亚洲性夜色夜夜综合| 国产成人啪精品午夜网站| 久久精品aⅴ一区二区三区四区| 精品欧美一区二区三区在线| 神马国产精品三级电影在线观看 | 最新美女视频免费是黄的| 亚洲美女视频黄频| 亚洲avbb在线观看| 五月玫瑰六月丁香| 一级毛片女人18水好多| 九色成人免费人妻av| 伦理电影免费视频| 中文字幕人妻丝袜一区二区| 精品久久久久久久久久免费视频| 久久99热这里只有精品18| e午夜精品久久久久久久| 日韩欧美 国产精品| 女警被强在线播放| 久久精品91蜜桃| 高清毛片免费观看视频网站| 精品不卡国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 国产精品免费视频内射| 亚洲五月天丁香| av福利片在线观看| 色精品久久人妻99蜜桃| 亚洲欧美精品综合一区二区三区| 久久精品91无色码中文字幕| 亚洲专区字幕在线| 亚洲aⅴ乱码一区二区在线播放 | 88av欧美| 国产精品亚洲美女久久久| 极品教师在线免费播放| 午夜日韩欧美国产| 1024视频免费在线观看| 午夜两性在线视频| 欧美成狂野欧美在线观看| 在线a可以看的网站| 18禁裸乳无遮挡免费网站照片| 香蕉久久夜色| 两性夫妻黄色片| 老司机靠b影院| 日韩欧美免费精品| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 婷婷亚洲欧美| 亚洲五月婷婷丁香| 麻豆久久精品国产亚洲av| a在线观看视频网站| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 日韩免费av在线播放| 亚洲国产精品久久男人天堂| 1024香蕉在线观看| 成人18禁高潮啪啪吃奶动态图| 欧美日韩亚洲综合一区二区三区_| 好看av亚洲va欧美ⅴa在| 国产黄色小视频在线观看| 亚洲 欧美一区二区三区| a级毛片a级免费在线| 黄色丝袜av网址大全| 日韩精品免费视频一区二区三区| 久久久国产成人精品二区| 波多野结衣高清无吗| 国产真实乱freesex| 亚洲自拍偷在线| 色综合欧美亚洲国产小说| 亚洲精品久久成人aⅴ小说| 可以免费在线观看a视频的电影网站| 久久99热这里只有精品18| 最新美女视频免费是黄的| 色综合亚洲欧美另类图片| 国产亚洲精品一区二区www| 我要搜黄色片| 在线国产一区二区在线| 又大又爽又粗| 麻豆久久精品国产亚洲av| 精品熟女少妇八av免费久了| 久久99热这里只有精品18| 黑人欧美特级aaaaaa片| 亚洲男人天堂网一区| www日本在线高清视频| 精品久久久久久久末码| 亚洲成人免费电影在线观看| 首页视频小说图片口味搜索| 在线观看一区二区三区| 人妻久久中文字幕网| 在线观看www视频免费| 69av精品久久久久久| 精品不卡国产一区二区三区| tocl精华| 国产男靠女视频免费网站| 国产单亲对白刺激| 色噜噜av男人的天堂激情| 一本精品99久久精品77| 变态另类丝袜制服| 国产高清videossex| 99久久精品热视频| 成人18禁高潮啪啪吃奶动态图| 亚洲无线在线观看| 国产av一区在线观看免费| 黑人巨大精品欧美一区二区mp4| 97超级碰碰碰精品色视频在线观看| 久久精品人妻少妇| 后天国语完整版免费观看| 国产精品一区二区精品视频观看| 真人一进一出gif抽搐免费| 亚洲全国av大片| 国产亚洲欧美98| 村上凉子中文字幕在线| 午夜福利成人在线免费观看| 日韩欧美三级三区| 国内精品久久久久精免费| 麻豆国产97在线/欧美 | 两性夫妻黄色片| 精品人妻1区二区| 欧美大码av| 搡老妇女老女人老熟妇|