999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

氮添加對中國草地生物量和土壤有機碳含量影響的Meta分析

2025-08-03 00:00:00劉暢陳積山朱瑞芬孫萬斌姚博董世魁
草地學報 2025年7期
關鍵詞:施氮土壤有機氮量

中圖分類號:S153 文獻標識碼:A 文章編號:1007-0435(2025)07-2078-1

doi:10.11733/j.issn.1007-0435.2025.07.004

引用格式:,等.氮添加對中國草地生物量和土壤有機碳含量影響的Meta分析[J].草地學報,2025,33(7):2078-2089LIU Chang,CHEN Ji-shan,ZHURui-fen,et al.Meta-analysis of the Impactsof Nitrogen Adition on Biomas andSoil Organic Carbon Content of China Grasslands[J].Acta Agrestia Sinica,2O25,33(7):2O78-2089

Meta-analysis of the Impacts of Nitrogen Addition on Biomass and Soil Organic Carbon Content of China Grasslands

LIU Chang1.2.3,CHEN Ji-shan23, ZHU Rui-fen2.3, SUN Wan-bin2,3, YAO Bo2,3,DONG Shui-ku1il4* (1.VeterinaryMedicineand Qinghai AcademyofAnimal Science,VeterinaryMedicineand AcademyofAnimalScience, QinghaiUniversityXining,QinghaiProvince8016,China;2.ChongqingAcademyofAnimalScience,Chongqing4246Cina; 3.PrataculturalEneigndtcoloesearchcnterofogig,ogg246ina;4oofrassandScieein ForestryUniversity,BeijinglOoo83,China)

Abstract:In order to explore the effects ofdifferent nitrogen addition states on grassland biomass,soil organic carbon(SOC),and microbial biomass carbon(MBC),a meta-analysis method was used to select 65 docu ments and 416 sets of data,and the effects of nitrogen addition conditions(Duration of Naddition,Rate of N addition and Duration of N addition),grassland types (Alpine grassland/meadow,desert grassand,typical grassland/meadow,meadow grassand),and climate factors (mean annual temperature (MAT),mean annual precipitation(MAP)) on the aboveground and belowground biomass,SOC and MBC of grassland ecosystems were investigated through subgroup analysis.The results showed that nitrogen application significantly affected grassland biomass and SOC;both urea and NH4NO3 nitrogen application significantly increased aboveground and underground biomass and SOC content ( .Plt;0.05 ;the durationofNaddition was more than5years,and the biomass,SOC and MBC contents increased. In grasslands with MAT gt;0°C ,aboveground and underground biomass and SOC content increased after nitrogen application. Therefore,when conducting nitrogen addition in grassland ecosystems,the effects of both biotic and abiotic factors on soil organic carbon should be considered.The research findings have important implications for the scientific formulation of grassland conservation and management strategies.

Key words:Nitrogen adition;Soil organic carbon; Grassland ecosystem; Biomass;Mean annual temperature; Mean annual precipitation

近年來,全球大氣氮沉降速率逐年增加,預計到2050年將比工業革命前增加近3倍1]。我國作為全球三大氮沉降區之一[2],大氣氮沉降對生物多樣性喪失、土壤酸化[3]和陸地生態系統中的碳和氮循環產生了嚴重影響[4-5]

在我國,草地面積約占國土面積的 41%[6] ,草地生態系統在全球碳氮循環和氣候變化中發揮重要作用,草地碳儲量約占我國陸地生態系統總碳儲量的 40%17 。研究表明,氮沉降(氮添加)會影響植物群落多樣性[8]、土壤理化性質[9-10]、土壤微生物群落結構、功能及作用過程[11-12]。氮沉降對土壤有機碳的影響呈現多元化,既能激發正效應,促進有機碳積累,也能表現出中性作用,對有機碳庫影響甚微,此外也會加速有機碳分解[13]。有研究表明,氮沉降(氮添加)抑制了微生物活性[14],增加土壤碳氮礦化速率[15],進而引起草地生態系統結構組成與功能的改變[12.16]。但長期氮添加會提高微生物活性,加速有機養分的分解,從而使有機碳含量下降[17-18]。氮輸入使得有機養分的輸人與有機養分的分解相平衡,從而對土壤有機碳含量無顯著影響[19]。氮沉降(氮添加)通過增加土壤有效氮,增加植物比葉面積、促進植物光合酶的合成,進而提高植物光合作用,提高草地地上生物量[20-21]和調落物量,從而增加土壤有機碳的含量[22]。這種促進作用在草地生態系統中尤為顯著,因為草地植物對氮素的需求較高。然而Fay等23]研究發現,氮添加對草地地上生物量的影響存在時間差異,短期內生物量增加,而長期觀察結果顯示生物量可能減少[24]。

氮添加對草地土壤有機碳的影響是一個復雜的多因素調控過程,包括施氮量、施氮時間、土壤性質、植物群落和微生物群落等多個方面。目前,不同草地類型對氮添加的響應是否存在一致的碳積累模式尚不清楚。特別是對氮添加狀態的不同(氮素處理持續時間、年施氮量和施氮類型)的響應尚不清楚。此外,不同地區之間的年均溫、年降雨量的差異會改變植物和土壤環境,并對土壤碳庫產生不同的影響[25-26]。王譽陶[27]研究發現,降雨改變了植物多樣性、土壤有機碳儲量和土壤微生物生物量碳。Chen等[28]研究發現,增溫可增加植物生物量,植物多樣性的提高促進了生態系統固碳。

我國草地類型多樣,不同的地理區域、氣候條件、土壤性質發生改變,都會使草地土壤有機碳含量發生變化[29-30]。因此,研究不同施氮條件下土壤性質、草地類型、溫度和降水等對土壤有機碳的影響對草地生態系統的適應性管理至關重要。本研究基于全國所展開的模擬氮沉降野外試驗,通過Meta分析草地土壤有機碳對氮沉降的響應特征,并探討了不同氮添加狀態(氮素處理持續時間、年施氮量和施氮類型)草地類型和氣候因子(年均溫、年降雨量)下土壤有機碳對氮添加的響應。為未來氮沉降背景下草地生態系統固碳潛力的評估提供科學依據,為不同草地類型固碳增匯能力的適應性管理提供理論基礎。

1 材料與方法

1. 1 數據收集

本研究利用中國知網、維普期刊和萬方數據庫對已發表的文章進行搜索,以“(草地 + 草場十草原 + 草甸)和(氮添加 + 氮肥 + 氮素添加 + 氮沉降 + 施氮十尿素 +NH4NO3) ”為主題詞進行期刊論文檢索,獲取2023年9月1日前發表的研究論文65篇,獲取416組數據。匯總數據的生物量指標包括地上生物量(Above-groundbiomass,AGB)地下生物量或根生物量(Below-groundbiomass,BGB),土壤有機碳指標包括微生物生物量碳(Microbialbiomasscarbon,MBC)、土壤有機碳(Soilorganiccarbon,SOC)。使用以下標準篩選合適的文獻(1)試驗數據基于中國天然草地,實驗方式為野外試驗,排除了室內培養試驗,且試驗持續時間確保不少于一年,以準確捕捉和評估草地生態系統中氮沉降的長期影響;(2)同一野外實驗中,實驗組和對照組必須在實驗地點、土壤質地以及氣候特征方面保持一致;(3)在相同的時間和空間尺度上,明確記錄了對照組和氮添加組的氮添加狀態(氮素處理持續時間、年施氮量和施氮類型),且至少有3個重復;

(4)文獻記錄了生物量或有機碳的指標;(5)對于多因素處理實驗,只選取對照和氮添加處理,避免其他交互實驗處理數據的影響;(6)實驗組和對照組中變量參數的平均值、標準差或標準誤以及樣本量可以直接從圖、表或者文字中進行提取。

此外,還獲取了環境數據和其他相關信息,包括試驗地點(經度、緯度和海拔)、年均溫(Meanannualtemperature,MAT)、年平均降水量(Meanannualprecipitation,MAP)、施氮類型(FormofNaddition)、年施氮量(Rate ofNaddition) (kg?ha-1?yr-1 N)和氮素處理持續時間(DurationofNaddition)、草地類型(Grasslandtype)。文獻中未給出試驗地點的年均溫和年平均降雨量信息時,通過全球氣候數據庫獲取(http://www.worldclim.org/)。通過GetDate(http://getdata-graph-digitizer.com/)獲取試驗數據中的平均值、標準差(或標準誤)。若文獻中提供數據為標準誤(SE),使用公式將其換算成標準差(SD):

式中: n 為試驗的重復次數。

圖1文獻研究的地理分布圖FiglGeographical distribution of the literature study areas.注:該圖基于國家基礎地理信息公共服務平臺下載的審圖號為GS (2024)0650號的標準地圖制作,底圖無修改 Note:Themap isbased onthe standard mapwith reviewnumberGS (2024)O65O downloaded from theNational Basic Geographic Informa tionPublic Service Platform,and the base map is not modified

1.2 數據分析

采用效應比值的自然對數(LRR)來計算施氮對生物量和土壤有機碳含量影響的效應[31]:

式中, Xt 為處理值的平均值, Xc 為對照組的平

均值。

每個個體效應大小的方差(u)計算如下:

式中, nt 和 nc 分別為實驗組和對照組的樣本量,St 和 Sc 分別為實驗組和對照組的標準差。

隨機效應模型應用限制最大似然法(REML)對加權響應比及置信區間進行計算:

wi*=1/(vi2

式中, τUi* 為單個研究的權重, vi 為研究內方差,τ2 為研究間方差。

累積效應值計算如下:

式中, 累計效應值, yi :單個研究的自然對數響應比。

總體標準誤計算如下:

式中,SE為總體標準誤。

累計效應值的 95% 置信區間計算如下:

當累計效應值 95% 置信區間(Confidenceinterval,CI均大于O,則認為施氮對土壤有機碳的影響顯著 ?Plt;0.05) ;若 95% 置信區間均小于0,則認為施氮對土壤有機碳的影響無顯著增加作用;若 95% 置信區間包含0,則認為施氮對土攘有機碳的影響無顯著影響 (Pgt;0.05)[31]

通過Q-test方法檢驗效應值整體異質性 Γ(Qt) 以及解釋變量對效應值的影響 (Qm[32] 。當效應值的整體異質性很強(且 Plt;0.05) 或者I系數 gt;60% 時,說明各施氮試驗的效應值受其他相關因素影響顯著[-34],引入解釋變量[35],使用混合效應(Mixed-effect)模型來計算加權效應大小并用R語言(version4.4.2)“metafor\"包中的“rma.mv\"函數進行亞組分析[36]。引人分類型變量氮素處理持續時間(DurationofNaddition)施氮量(RateofNaddition)施氮類型(FormofNaddition)、草地類型(Grasslandtype)、年均溫(MAT)年平均降水量(MAP)對數據進行亞組分析。草地類型包括高寒草原(草甸)(Alpinegrassland/meadow)荒漠草原(Desert grassland)、典型草原(草甸)(Typicalgrassland/meadow)、草甸草原(Meadowgrassland);施氮類型為尿素(Urea)、硝酸銨 (NH4NO3) ;施氮期限分為短、長,分別為 ?5 年、 .gt;5 年;施氮量分為3個水平: ? 60,60~120 和 gt;120kg-1?ha-1?yr-1N ;年均溫分為 ? 0,0~5 和 1gt;5°C ;年均降雨量分為 ?400 和 gt;400mm 。

利用多元回歸(rma.glmulti)分析各因子對效應值影響的大小。

2 結果與分析

發表性偏倚檢驗:本研究采用Egger's回歸分析偏倚檢驗。 Pgt;0.05 表示漏斗圖形狀對稱,未受到發表性偏倚的影響或受其影響很小,結果可行性很強, Plt;0.05 提示可能存在發表偏倚。如果失安全系數(fail-safenumber,Nfs)大于臨界值 (5n+10 n 為文獻中提取到的數量的組數),則表示無論是否存在潛在的發表偏倚,說明結果可靠,受偏倚性的影響較小[31,37-38]

用R語言(version4.4.2)采用分段結構方程模型\"piecewiseSEM\"包構建不同影響因子對地上、地下生物量、土壤有機碳、微生物生物量碳的結構方程模型[39]

2.1生物量和有機碳對響應變量發表性偏倚檢驗

經Egger檢驗(表1),施氮對地下生物量和土壤有機碳效應的發表性偏倚檢驗漏斗圖的整體對稱性顯著 (Plt;0.05) ,且失安全系數遠高于臨界值(5n+10,n 為文獻中提取到的數量的組數),綜合兩種檢驗結果,施氮對地下生物量和土壤有機碳的效應不存在發表偏移。地上生物量和微生物生物量碳的失安全系數遠高于臨界值,表明研究結果不存在顯著發表偏移。

表1發表性偏倚檢驗結果 Table1Results of publication bias test
注:n為文獻中提取到的數量的組數: ;P 為漏斗圖的非對稱性檢驗結果;Nfs為失安全系數 Note:n is the number of groups for the quantities extracted from the literature; P is the result of the asymmetry test for the funnel plot;Nfs is the Rosenthal's fail-safe number

2.2施氮對生物量和土壤有機碳的影響

施氮對草地地上、地下生物量、土壤有機碳含量、微生物生物量碳含量的效應值如圖2所示,施氮對地上生物量的平均效應值為 0.3149(95% 的置信區間為 0.2804~0.3493) ,地下生物量的平均效應值為 0.1516(95% 的置信區間為 0.1018~0.2014) ,土壤有機碳含量的平均效應值為 0.0801(95% 的置信區間為 0.0481~0.1121) ,上述的效應值均為正,且 95% 置信區間均大于0。表明施氮顯著增加了草地地上、地下生物量以及土壤有機碳含量( Plt; 0.05)。微生物生物量碳含量的平均效應值為0.024×95% 的置信區間為 -0.0803~0.1284) ,雖然效應值為正,但是置信區間包含0,說明施氮可能對土壤微生物生物量碳的含量無顯著影響。

表2顯示的異質性檢驗結果表明, I2 均大于60% ,異質性檢驗均達顯著水平( Plt;0.0001) ,說明施氮對生物量、土壤有機碳和微生物生物量碳還受其他相關因素影響,需要納入其他因素進行Meta亞組分析以分析其異質性來源。

Meta亞組分析結果表明,在各影響因素中(表3),施氮對地上生物量的影響受氮素處理持續時間、年施氮量、施氮類型、草地類型、年降雨量和年均溫的影響達到極顯著水平( Plt;0.001) ;施氮對地下生物量的效應受氮素處理持續時間、草地類型、年降雨量和年均溫的影響比較顯著 (Plt;0.01) ;施氮對土壤有機碳的效應受年施氮量的影響較顯著(Plt;0.01) ,受氮素處理持續時間、施氮類型、年降雨量和年均溫的影響較顯著 ?Plt;0.05) ;施氮對土壤微生物生物量碳的效應受施氮類型、草地類型的影響較顯著 (Plt;0.01) 。

2.3草地生物量對氮添加的響應

氮素處理持續時間5年以上,地上生物量和地下生物量分別增加了49. 52% 和 26.68% (圖3),處理時間 ?5 年時,地上生物量和地下生物量可增加29.39% 和 10.30% 。地上生物量和地下生物量隨著施氮量的增加隨之增加,施氮量 ?60kg?ha-1?a-1N 地上生物量和地下生物量分別增加了21. 64% 和

10.89% ;施氮量在 60~120kg?ha-1?a-1N 時地上生物量和地下生物量分別增加了 31.02%.14.11% 施氮量 gt;60kg?ha-1?yr-1N 地上生物量和地下生物量分別增加了 40.57% 和 18.55% 。在施加尿素時,地上生物量和地下生物量分別增加了 31.41% 和17.81% ;施用 NH4NO3 地上生物量和地下生物量分別增加了 32.04% 和 11.02% 。在典型草原/草甸中,施氮后地上生物量和地下生物量分別增加了16.80% 和 16.36% ;草甸草原中,施氮可增加地上生物量 50.95% ,但是地下生物量降低了 9.67% ;施氮增加了荒漠草原地上生物量 60% 和地下生物量 21.16% ;施氮增加了高寒草原(草甸)地上生物量 34.50% 和地下生物量 23.64% 。

年降雨量大于 400mm 時,施氮分別增加了地上生物量和地下生物量 26.02% 和 6.84% (圖4);年降雨量小于 400mm 時,地上生物量和地下生物量分別增加了 35.45% 和 19.27% 。年均溫小于 0°C 時,施氮分別增加了地上生物量和地下生物量47.26% 和 0.17% ;年均溫在 0°C~5°C 時,施氮分別增加了地上生物量和地下生物量 16.42% 和18.79% ;年均溫大于 5°C 時,地上生物量和地下生物量均達到最大分別是 48.80% 和 29.97% O

圖2氮添加影響草地地上生物量、地下生物量、土壤有機碳、微生物生物量碳含量效應值Fig.2Forest plots of the effect sizes of nitrogen addition affecting above-and below-ground biomass,soil organic carbon,andmicrobialbiomasscarboncontent

注:紅色實線為平均效應值Note:Red solid line,Effect value

表2異質性檢驗結果Table2 Results of heterogeneity test
注: Qt 為異質性檢驗的統計值: P 為異質性檢驗的顯著性;I為 研究間方差占總體方差的比例: ?τ2 為研究間方差 Note: Qt ,total heterogeneity; P ,Significance of heterogeneity test; I2 ,residual heterogeneity/unaccounted variability; τ2 ,estimated amount of residual heterogeneity

2.4土壤有機碳對氮添加的響應

施氮可以增加土壤有機碳 (8.01%) 和微生物生物量碳含量 (2.40%) 。氮素處理持續時間5年以上,土壤有機碳含量增加了 6.17% ,微生物生物量碳含量降低了 28.54% (圖5);處理時間 ? 5年時,可增加土壤有機碳和微生物生物量含量10.60% 和 7.59% 。施氮量 ?60kg?ha-1?yr-1N 土壤有機碳含量增加了 6.25% ,微生物生物量碳含量降低了 0.37% ;施氮量在 60~120kg?ha-1?yr-1N 時土壤有機碳和微生物生物量碳含量分別增加了 12.64% 和 0.34% ;施氮量 gt;60kg?ha-1?yr-1N 土壤有機碳含量增加了 11.77% ,微生物生物量碳含量降低了 9.65% 。在施加尿素時,土壤有機碳含量增加了 11.85% ,微生物生物量碳含量降低了 25.69% ;施用 NH4NO3 土壤有機碳和微生物生物量碳含量分別增加了 7.22% 和 27.08% 。在典型草原/草甸中,土壤有機碳含量增加了14.80% ,但是微生物生物量碳含量降低了24.32% ;草甸草原中,施氮可增加土壤有機碳含量 3.79% 、微生物生物量碳含量 67.24% ;施氮增加了荒漠草原土壤有機碳含量 8.21% ,但是微生物生物量碳含量降低了 7.35% ;施氮增加了高寒草原/草甸土壤有機碳含量 17.97% 、微生物生物量碳含量 8.97% 。

年降雨量大于 400mm 時,施氮分別增加了土壤有機碳和微生物生物量碳含量 11.49% 和 5.06% (圖6);年降雨量小于 400mm 時土壤有機碳增加了8.02% ,微生物生物量碳含量降低了 6.76% 。年均溫小于 0°C 時,施氮分別增加了土壤有機碳和微生物生物量碳含量 4.64% 和 4.64% ;年均溫在 0°C~5°C 時,施氮增加了土壤有機碳含量 14.04% ,微生物生物量碳含量降低了 6.71% ;年均溫大于 5°C 時,施氮分別增加了土壤有機碳和微生物生物量碳含量 10.09% 和 4.63% O

表3各影響因素的Meta亞組分析統計檢驗結果Table 3Results of meta subgroup analysis for each influencing factor
注:‘***”‘**”**’‘△\"NS'分別表示統計檢驗的 P 值 *lt;0.001mlt;0.01\"*lt;0.05\"*lt;0.1mgt;0.1\" 專Note: ************…Δ and ‘NS'indicate ≤0.001,Plt;0.01,Plt;0.05,Plt;0.1,Pgt;0.1, respectively
圖3不同草地類型和氮添加狀態(氮素處理持續時間、 年施氮量和施氮類型)下氮添加對草地生物量的影響 Fig.3Effectsof nitrogenaddition on grassland biomass underdifferent grassland types and nitrogen addition statuses (DurationofNaddition,RateofNadditionandDurationof Naddition)

多元回歸分析表明,影響地上生物量的主要因素有施氮量、草地類型、年均溫和年降雨量(圖7a)。施氮時間是影響地下生物量的主要因子(圖7b)。各因子對土壤有機碳的影響均不顯著,其中施氮量以及年均溫對土壤有機碳影響較大(圖7c)。施氮類型、草地類型和年均溫顯著影響微生物生物量碳(圖7d)。

圖4不同氣候因子下氮添加對草地生物量的影響Fig.4Effects of nitrogen addition on grassland biomassunderdifferentclimaticfactors

結構方程表明,MAT直接影響地上生物量 (Plt; 0.01),直接效果為0.211(圖8a);MAP對地上生物量有極顯著直接影響 (Plt;0.01) ,直接效果為0.308,因而MAP是地上生物量的最重要環境決定因子。施氮類型、施氮量、氮添加時間與地上生物量之間具有顯著的通徑關系,且均為顯著正相關關系。MAT對地下生物量有顯著負影響(圖8b),影響效果為一0.328。施氮類型對地下生物量有顯著的正效應,施氮量、氮添加時間對地下生物量有顯著的負影響。施氮中氮素類型是影響草地地上、地下生物量的主要決定因子。MAT和MAP通過影響草地類型影響微生物生物量碳含量(圖8c),施氮類型對微生物生物量碳之間有顯著的負效應,直接效果為一0.233,MAT和MAP與微生物生物量碳之間無顯著通徑關系。MAT對土壤有機碳含量有顯著負效應(圖8d),施氮類型和草地類型與土壤有機碳含量之間具有顯著的正效應,施氮時間與土壤有機碳含量之間無顯著通徑關系。

圖5不同草地類型和氮添加狀態(氮素處理持續時間、年施氮量和施氮類型)下氮添加對草地有機碳、微生物生物量碳含量的影響Fig.5Effects of nitrogen addition on organic carbon and microbial biomass carbon content ofgrasslandunderdifferent grassland typesand nitrogenadditionstatuses(Duration ofN addition,RateofNadditionand Duration ofNaddition)
圖6不同氣候因子下氮添加對草地有機碳、微生物生物量碳含量的影響Fig.6Effects of nitrogen addition on organic carbon and microbial biomass carbon content of grassland underdifferent climatic factors

圖7影響因子對地上、地下生物量、土壤有機碳、微生物生物量碳的相對重要性排序

Fig.7Ranking of relative importance of impact factors on above-and below-ground biomass,soil organic carbon,and microbial biomass carbon

注:Type為草地類型;Form為施氮類型;Rate為施氮量;Duration為 氮添加時間;MAT為年均溫;MAP為年平均降水量。紅線為相對 重要線,數值 gt;0.8 表示該因子對效應值影響顯著

Note:Type,Grassland type;Form,FormofNaddition;Rate,Rate of N addition;Duration,Duration ofNaddition;MAT,Annual MeanTemperature;MAP,MeanAnnualPrecipitation.Theredline istherelativeimportance line,thevalue gt;0.8 indicates that the factor has a significant effect on the effect value

3討論

3.1施氮對生物量、土壤有機碳的影響

氮添加刺激了植物的生長,提高植物地上地下生物量(圖2a和b),可能是由于氮添加后土壤中的氮有效性得到改善,使得植物能夠攝取到更為充足的氮養分,進而為其自身生長提供強勁動力[40-42]。進而有利于土壤有機碳的形成與累積過程[43-45]。然而,植物生物量的增加也可能是氮添加引發微生物群落結構與功能的改變,這些微生物群落層面的變動同樣可能作用于植物的生長過程,間接促使植物生物量上升[46-47]

施肥持續時間、施氮量以及氮肥添加種類的不同對草地土壤有機碳的影響不同[48],表明植物和微生物對氮的偏好存在差異。這與本研究結果一致,施氮量和時間是影響生物量的重要因素[49-50],短期氮添加量增加了微生物生物量碳,可能是由于減輕了微生物對氮的限制[51-52]。 NH4NO3 和尿素對生物量、有機碳、微生物生物量碳含量的響應模式不一致[53],可能是由于它們的化學性質、分解速率以及對土壤微生物活動的不同影響。 NH4NO3 在土壤中的分解速率較快, NH4+ 和 NO3- 可以直接參與植物和微生物的代謝過程。這種快速的氮素供應方式可能導致短期內地上生物量和微生物生物量碳含量的增加。本研究發現,施加尿素后土壤有機碳含量顯著增加,尿素在分解過程中會釋放碳元素,這些碳元素可以被土壤微生物利用,促進微生物的生長和活動,進而影響有機碳的分解和轉化。 NH4NO3 主要影響無機氮的供應,對有機碳的直接影響較小[52-54]

"

施氮量和氮素處理持續時間是調節生物量的重要因素[49],其中施氮量對地上生物量的影響更為重要(圖7a)。本研究結果顯示植物生物量、土壤有機碳和微生物生物量碳對氮添加持續時間和氮添加量的響應有差異,表明草地存在氮飽和效應。植物的地上和地下生物量以及土壤微生物的數量和活動都有不同的氮飽和閾值,植物和微生物的氮磷比例也不平衡[55-56],如Yao等[57]在對羊草草地和Zhang等[58]在黃土高原白羊草研究中發現的氮添加的閥值效應。本研究中也發現了添加氮超過120kg?hm-2?yr-1 后土壤微生物生物量碳降低。前人關于氮沉降對草地研究表明,當草地生態系統未達到氮飽和時,適量的氮素沉降能夠降低土壤微生物的營養壓力,過量的氮沉降將增加土壤環境變化對微生物群落的脅迫程度,從而使微生物的生長和活性受到抑制[59-60]

長期氮添加通過緩解植物的養分限制,提高了植物地上生物量及凋落物量[61]。Avolio等[62]的實驗表明長期氮添加對草地植物群落結構有非常顯著的影響,景明慧等63在進行長期的氮添加試驗發現地上生物量增加了。本研究發現施氮時間高于5年可增加草地生物量、有機碳和微生物生物量碳含量,與Bai和Lan等人[64-65]的研究結果一致。但長期氮添加導致土壤酸化進而導致土壤微生物生物量碳含量隨時間降低[66],長期氮添加累計的酸化效應也會進一步損害土壤微生物和微生物介導的土壤呼吸[67]。

3.2氣候因子調控氮添加對生物量、土壤有機碳的影響

水熱因子是影響草地生態系統生產力的主要因素[68]。本研究發現平均溫度和降水量對生物量和土壤有機碳有顯著影響 (Plt;0.05) (圖4、圖6)。氣候決定了草地類型,即植物群落組成和生物量,從而決定了微生物生物量碳。meta分析顯示,氮的添加對所有類型草地的植物生物量和有機碳含量均有持續的提高(圖2,圖5),荒漠草原的地上生物量對氮添加的響應最大,高寒草甸的地下生物量對氮添加的響應最大。施氮增加了高寒草甸植物地上和地下產量[69]和有機碳含量[70-71]。可能是由于荒漠草原地區氣候干燥、土壤缺乏營養物質,水分和養分是促進植物生長的主要因子[72],施氮解決了土壤氮素缺乏的狀態[73],增加了土壤營養從而地上生物量得到增加。而在降水和氮添加的共同作用下,一定程度上解決了缺水帶來的影響,水分更好的激發了氮素的肥力,使氮肥得到更有效的釋放,增加有機質含量[72]。在年均溫或年降水量較高的區域,土壤中的養分含量豐富,導致土壤趨向于“氮素充足\"甚至是“氮素過飽和”的狀態[74]。在這種情況下,追加氮肥反而會抑制微生物的代謝活動,進而導致土壤中有機碳的顯著累積。研究表明高寒草甸施氮降低微生物生物量碳[75],與我們得到的結論一致。

水分是影響草地生物量的關鍵因素[76],鐘澤坤[77的研究表明,降雨增加可以提高土壤有機碳含量,賀云龍等的研究表明,施氮下土壤微生物生物量碳隨著降水的豐富顯著增加。我們的研究發現,降水增加會導致土壤有機碳和微生物生物量碳含量的增加,這可能是降水增加使地表水文格局改變,調節植物生長及減緩枯落物分解過程,進而影響土壤碳庫儲量及碳固定過程[79-81]。降水增加了土壤中水分土壤水分有效性,從而提高了土壤微生物活性和植物生長[82],增加陽離子的可用性[83],促進了土壤微團聚體的形成和穩定[84],為土壤有機碳提供了物理保護,避免被微生物所降解,從而增加土壤有機碳含量[85-86]。

增溫可以刺激植物生長[87-88]。本研究Meta分析結果顯示隨著溫度的升高增加了草地生物量(圖4)。適宜的溫度有利于植物的生長和光合作用,從而促進生物量的積累和有機碳的固定。適宜的氣候條件有利于生物量的增加和土壤有機碳的積累,而極端的氣候條件可能會對生態系統造成負面影響。因此,在氮沉降背景下,了解氣候因子對生物量和土壤有機碳的影響,對于草地生態系統的保護修復和可持續發展具有重要意義。

4結論

氮添加可以增加不同草地類型的地上、地下生物量以及土壤有機碳含量。同時,隨著施氮量的增加,顯著增加了草地地上生物量的效應值。施氮年限的增加顯著增加了地上、地下生物量的效應值。年均溫和年平均降水量會顯著影響草地生物量和土壤有機碳對氮添加的響應,對土壤微生物生物量碳無顯著影響。因此,在草地生態系統進行氮添加時,應考慮生物和非生物因素對土壤有機碳的影響。本研究為未來氮沉降背景下草地生態系統固碳潛力的評估提供科學依據,也為不同草地類型固碳增匯能力的適應性管理提供理論基礎。

參考文獻

[1]KANAKIDOUM,MYRIOKEFALITAKISS,DASKALAKIS N,et al. Past present and future atmospheric nitrogen deposition[J]. Journal of the Atmospheric Sciences, 2016,73 (5):2039-2047

[2]LIUL,ZHANG T,GILLIAMF S,et al. Interactive Effcts of nitrogen and phosphorus on soil microbial communities in a tropical forest[J].Plos One,2013,8(4):1-10

[3]魏芙蓉.氮添加影響白羊草草地土壤碳氮轉化的機制[D]. 楊凌:西北農林科技大學,2023:3-11

[4]GALLOWAY JN,TOWNSEND AR,ERISMAN JW,et al.Transformation of the nitrogen cycle:recent trends,questions,and potential solutions[J].Science,20O8,32o(2878): 889-892

[5] CHENG Y,WANG J,GE Z,et al. Background nitrogen deposition controls the efects of experimental nitrogen addition on soil gross N transformations in forest ecosystems[J].Biogeochemistry,2020,151(1):335-341

[6]XUL,YUGR,HE NP,et al. Carbon storage in China’s terrestrial ecosystems:A synthesis[J].Scientific Reports,2018, 8(1):1-13

[7]TANG X L,ZHAO X,BAI Y F,et al. Carbon pools in China's terrestrial ecosystems:New estimates based on an intensive field survey[J].Proceedings of the National Academy of Sciences,2018,115(16):4021-4026

[8] STEVENS C J, DISE N B,MOUNTFORD J O,et al. Impact of nitrogen deposition on the species richness of grasslands[J]. Science,2004,303(5665):1876-1879

[9]LU XF,HOU EQ,GUO JY,et al.Nitrogen addition stimulates soil aggregation and enhances carbon storage in terrestrial ecosystems ofChina:A meta-analysis[J].Global Change Biology,2021,27(12):2780-2792

[10]張乃木,宋婭麗,王克勤,等.模擬氮沉降下滇中亞高山森林 凋落物養分元素釋放特征[J].生態環境學報,2021,30(5): 920-928

[11] LEFF JW,JONES S E,PROBER S M,et al. Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe[J].Proceeding of the National ofSciencesof theUnited StatesofAmerica,2015, 112(35):10967-10972

[12]劉春華,吳東梅,劉雨暉,等.氮沉降對米褚天然林土壤有機 碳及微生物群落結構的影響[J].林業科學研究,2021,34 (2):42-49

[13] XU C H,XU X,JU C H,et al.Long-term,amplified responses of soil organic carbon to nitrogen addition worldwide [J].Global Change Biology,2021,27(6):1170-1180

[14]秦潔.長期氮沉降對貝加爾針茅草原土壤微生物群落的影響 機制[D].北京:中國農業科學院,2021:53-60

[15]陳天,程瑞梅,沈雅飛,等.氮添加對三峽庫區馬尾松人工林 土壤團聚體有機氮組分和氮礦化的影響[J].應用生態學報, 2023,34(10):2601-2609

[16]ZAK D R,FREEDMAN Z B,UPCHURCH R,et al. Anthropogenic N deposition increases soil organic mattr accumulation without altering its biochemical composition [J]. Global Change Biology,2017,23(2):933-944

[17]MACKMC,SCHUUREA,BRETHARTEMS,etal. Ecosystem carbon storage in arctic tundra reduced by long-term nutrient fertilization[J].Nature,2004,431(7007):440-443

[18]李煥茹,朱瑩,田紀輝,等.碳氮添加對草地土壤有機碳氮磷 含量及相關酶活性的影響[J].應用生態學報,2018,29(8): 2470-2476

[19]WANG C H,ZHU F,ZHAO X,et al. The effects of N and P additions on microbial N transformations and biomass on salinealkaline grassland of loess plateau of northern china[J].Geoderma,2014,213(1):419-425

[20]CHEN J,LUO Y,VAN GROENIGEN KJ,et al.A keystone microbial enzyme for nitrogen control of soil carbon storage[J]. Science Advances,2018,4(8):1689-1695

[21]WAN D,MA M,PENG N,et al. Efects of long-term fertilization on calcium-associated soil organic carbon:Implications for C sequestration in agricultural soils[J]. Science of the Total Environment,2021,772(1):1-9

[22] SHI Y,SHENG L X,WANG ZQ,et al.Responses of soil enzyme activity and microbial community compositions to nitrogen addition in bulk and microaggregate soil in the temperate steppe of Inner Mongolia[J].Eurasian Soil Science,2O16,49 (10):1149-1160

[23]FAY P, PROBER S,HARPOLE W,et al. Grassland productivity limitedbymultiple nutrients[J].NaturePlants,2015,1 (7):1-5

[24]SEABLOOM E,ADLERP,ALBERTI J,et al. Increasing effects of chronic nutrient enrichment on plant diversity loss and ecosystem productivity over time[J]. Ecology,2021,102(2): 1-11

[25]GRIMM NB,CHAPINFS,BIERWAGENB,et al.The impacts of climate change on ecosystem structure and function [J].Frontiers in Ecology and the Environment,2Ol3,11(9): 474-482

[26]WALTHER G R. Community and ecosystem responses to recent climate change[J].Philosophical Transactions of the Royal Society B-Biological sciences,2010,365(1549):2019- 2024

[27]王譽陶.降水變化與氮沉降對黃土高原典型草原生態系統多 功能性的影響[D].銀川:寧夏大學,2023:104-106

[28]CHENJ,LUO YQ,XIAJY,et al.Warming effects on ecosystem carbon fluxes are modulated by plant functional types [J].Ecosystems,2017,20(3):515-526

[29]傅華,陳亞明,王彥榮,等.阿拉善主要草地類型土壤有機碳 特征及其影響因素[J].生態學報,2004,24(3):469-476

[30]JINVL,HANEY RL,FAY PA,et al. Soil type and moistureregime control microbial C and N mineralization in grassland soils more than atmospheric CO2-induced changes in liter quality[J].Soil Biology and Biochemistry,2013,58(1): 172-180

[31] HEDGES L V,GUREVITCH J,CURTIS P S.The metaanalysis of response ratios in experimental ecology[J].Ecology,1999,80(4):1150-1156

[32]VIECHTBAUER W. Conducting meta-analyses in R with the meta for package[J]. Journal of Statistical Software,2010,36 (3):1-48

[33]BENITEZ-LOPEZ A,ALKEMADE R,SCHIPPER A M, et al.The impact of hunting on tropical mammal and bird populations[J]. Science,2017,356(6334):180-183

[34] ZHAO J,YANG Y D, ZHANG K,et al. Does crop rotation yield more in China?A meta-analysis [J].Field Crops Research,2020,245(1):107659-107668

[35]HIGGINS JPT,THOMPSON SG,DEEKSJJ,et al. Measuring inconsistency in meta-analyses[J].British Medical Journal,2003,327(7414):557-560

[36]VIECHTBAUER W. Conducting meta-analyses in R with the metafor package[J]. Journal of Statistical Software,2O10,36 (3):1-48

[37]EGGER M,SMITHG D,SCHNEIDER M,et al.Bias in meta-analysis detected by a simple,graphical test-Reply[J]. BritishMedicalJournal,1997,315(7109):629-634

[38]ERISMANBE,BOLSERDG,ILICHA,et al.A metaanalytical review of the effects of environmental and ecological drivers on the abundance of red snapper (Lutjanus campechanus)n the U.S.Gulfof Mexico[J]. Reviews inFish Biolgy and Fisheries,2020,30(2):437-462

[39]LEFCHECKJS.piecewiseSEM:Piecewise structural equation modeling in R for ecology,evolution,and systematics[J]. Methods in Ecology and Evolution,2016,7:573-579

[40]PREGITZERKS,BURTONAJ,ZAKDR,etal.Simulated chronic nitrogen deposition increases carbon storage in NorthernTemperate forests[J].Global Change Biology, 2007,14(1):142-153

[41] SCHULTE-UEBBING L,DE VRIES W. Global - scale impacts of nitrogen deposition on tree carbon sequestration in tropical,temperate,and boreal forests:A meta-analysis[J]. Global Change Biology,2017,24(5878):416-431

[42] YUE K,FORNARA D A,LI W,et al. Nitrogen addition affects plant biomass allocation but not allometric relationships among different organs across the globe[J].Journal of Plant Ecology,2020,14(3):361-371

[43]KELLERA B,BORERET,COLLINS SL,et al.Soil carbon stocks in temperate grasslands differ strongly across sites but are insensitive to decade-long fertilization[J].Global Change Biology,2022,28(4):1659-1677

[44]FORNARA D A,TILMAN D. Soil carbon sequestration in prairie grasslands increased by chronic nitrogen addition[J]. EcologyLetters,2012,93(9):2030-2036

[45]張丹,鄧健,朱運,等.氮添加對黃土丘陵區草地土壤團聚體 有機碳庫的影響[J].草地學報,2023,31(7):2031-2040

[46]SEABLOOMEW,KINKEL L,BORER E T,et al.Food webs obscure the strength of plant diversity effcts on primary productivity[J]. Ecology Letters,2017,20(4):505-512

[47]EBELINGA,STRAUSSAT,ADLERPB,etal.Nutrient enrichment increases invertebrate herbivory and pathogen damage in grassands[J]. Journal of Ecology,2021,110(2): 327-339

[48]肖勝云,董云社,齊玉春,等.草地生態系統土壤有機碳庫對 人為干擾和全球變化的響應研究進展[J].地球科學進展, 2009,24(10):1138-1148

[49] ZHANG T A,CHEN HY,RUAN H.Global negative effcts of nitrogen deposition on soil microbes[J].The ISME Journal,2018,12(7):1817-1825

[50]YANGY,LIT,POKHARELP,etal.Global effects on soil respiration and its temperature sensitivity depend on nitrogen additionrate[J].Soil Biologyand Biochemistry,2O22,174 (5):1-11

[51]MANZONI S,TAYLOR P,RICHTER A,et al. Environmental and stoichiometric controls on microbial carbon-use efficiency in soils[J]. New Phytologist,2012,196(1):79-91

[52]SPOHN M,POTSCH E M,EICHORST SA,et al.Soil microbial carbon use efficiency and biomass turnover in a longterm fertilization experiment in a temperate grassand[J].Soil Biology and Biochemistry,2016,97(1):168-175

[53]張祚恩,刁華杰,郝杰,等.不同氮化合物添加對農牧交錯帶 草地土壤凈氮礦化速率的影響[J].草地學報,2023,31(5): 1322-1330

[54]HARRISON K A,BOL R,BARDGETTRD.Preferences fordifferent nitrogen forms bycoexisting plant species and soil microbes[J].Ecology,2007,88(4):989-999

[55]PENG YF,LIF,ZHOU GY,et al. Linkages of plant stoichiometry to ecosystem production and carbon fluxes with increasing nitrogen inputs in an alpine steppe[J].Global Change Biology,2017,23(12):5249-5259

[56]YANGG J,STEVENS C,ZHANG ZJ,et al.Diferent nitrogen saturation thresholds for above-,below-,and total net primary productivity in a temperate steppe[J]. Global Change Biology,2023,29(16):4586-4594

[57]YAO M,RUI J,LI J,et al. Rate-specific responses of prokaryotic diversity and structure to nitrogen deposition in the leymus chinensis steppe[J].Soil Biologyand Biochemisty, 2014,79(1):81-90

[58] ZHANG JY,AI Z M,LIANGC T,et al. Response of soil microbial communities and nitrogen thresholds of Bothriochloa ischaemum to short-term nitrogen addition on the Loess Plateau [J].Geoderma,2017,308(1):112-119

[59]劉蔚秋,劉濱揚,王江,等.不同環境條件下土壤微生物對模 擬大氣氮沉降的響應.生態學報,2010,30(7):1691-1698

[60]張成霞,南志標.土壤微生物生物量的研究進展[J].草業科 學,2010,27(6):50-57

[61]刁勵瑋,李平,劉衛星,等.草地生態系統生物量在不同氣候 及多時間尺度上對氮添加和增雨處理的響應[J].植物生態學 報,2018,42(08):818-830

[62]AVOLIO M L,KOERNER S E,PIERRE KJL,et al. Changes in plant community composition,not diversity,during a decade of nitrogen and phosphorus additions drive aboveground productivity ina talgrass prairie[J]. Journal of Ecology,2014,102(6):1540-1553

[63]景明慧,賈曉彤,張運龍,等.長期氮添加對內蒙古典型草原 植物地上、地下生物量及根冠比的影響[J].生態學雜志, 2020,39(10):3185-3193

[64]BAI YF,WUJG,CLARK CM,et al. Tradeoffs and thresholds in the effects of nitrogen addition on biodiversity and ecosystem functioning:evidence from innerMongolia Grasslands [J].Global Change Biology,2010,16(1):358-372

[65]LAN Z C,BAI Y F. Testing mechanisms of N-enrichmentinduced species loss in a semiarid Inner Mongolia grassland: critical thresholds and implications for long-term ecosystem responses[J]. Philosophical Transactions of the Royal Society B-Biological Sciences,2012,367(1606):3125-3134

[66] CHEN D,LAN Z,HU S,et al. Effects of nitrogen enrichment on belowground communities in grassland:relative role of soil nitrogenavailabilitys.oilacidificationJ]SoilBolog and Biochemistry,2015,89:99-108

[67]WILCOTS M E, SCHROEDER K M,DELANCEY LC,et al.Realistic rates of N addition increase carbon flux rates but do not change soil carbon stocks in a temperate grassand.Global ChangeBiology[J],2022,28,4819-4831

[68]白永飛.降水量季節分配對克氏針茅草原群落初級生產力的 影響[J].植物生態學報,1999,23(2):60-65

[69]LIU Y W,XU R,XU XL,et al. Plant and soil responses of an alpine steppe on the Tibetan Plateau to multi-level nitrogen addition[J].Plant and Soil,2013,373(1):515-529

[70]LIJH,HOU YL,ZHANG SX,et al.Fertilization with nitrogen and/or phosphorus lowers soil organic carbon sequestration in alpine meadows[J]. Land Degradation and Development,2017,29(6):1634-1641

[71]GUO Q,HU Z,LI S,et al. Nitrogen-induced changes in carbon fluxes are modulated by water availability in a temperate grassland[J]. Journal of Geophysical Research Biogeosciences,2021,126(12):1-14

[72]霍利霞,紅梅,趙巴音那木拉,等.氮素和水分添加對荒漠草 原土壤理化性質及微生物量碳氮的影響[J].北方農業學報, 2018,46(3):54-59

[73]何宇,盛茂銀,王軻,等.土地利用變化對西南喀斯特土壤團 聚體組成、穩定性以及C、N、P化學計量特征的影響[J].環境 科學,2022,43(7):3752-3762

[74]陳潔,駱土壽,周璋,等.氮沉降對熱帶亞熱帶森林土壤氮循 環微生物過程的影響研究進展[J].生態學報,2020,40(23): 8528-8538

[75]王麗娜,羅久富,楊梅香,等.氮添加對退化高寒草地土壤微 生物量碳氮的影響[J].草業學報,2019,28(7):38-48

[76]BAI Y,WUJ,PAN Q,et al.Positive linear relationship between productivity and diversity:Evidence from the Eurasian Steppe[J].Journal ofAppliedEcology,2007,44(5):1023- 1034

[77]鐘澤坤.增溫和降雨改變對黃土丘陵區擢荒草地土壤碳循環 關鍵過程的影響[D].楊凌:西北農林科技大學,2021: 108-109

[78]賀云龍,齊玉春,彭琴,等.外源碳和氮輸入對降水變化下土 壤呼吸的短期影響[J].環境科學,2018,39(4):1934-1942

[79]CRAFTCB,WASHBURNC,PARKERA.Latitudinal Trendsin Organic Carbon Accumulationin Temperate FreshwaterPeatlands[J].SpringerNetherlands,2OO8,23(3):23-31

[80]LOVELOCKCE,ADAMEMF,BENNIONV,etal.Contemporary rates of carbon sequestration through vertical accretionof sedimentsin mangrove forestsand saltmarshesof south east queensland,Australia[J].Estuariesand Coasts,2Ol4,37 (3):763-771

[81]SHANG LN,ZHANG ZS,SONG XL,etal.Carbon sequestration and nutrient accumulation(N and P)in two typical wetlands in SanjiangPlain,Northeast China[J].Fresenius EnvironmentalBulletin,2015,24(2):422-428

[82]DIAOH,CHENX,ZHAOX,et al.Effects of nitrogen addition and precipitation alteration on soil respiration and its componentsinasaline-alkalinegrassland[J].Geoderma,2O22,406 (1):1-11

[83]BAITS,WANGP,HALLSJ,etal.Interactive global change factorsmitigate soil aggregation and carbon change ina semi-arid grassland[J].Global Change Biology,2O2O,26(9): 5320-5332

[84]BLANCO-CANQUI H,LAL R.Mechanisms of carbon sequestration in soil aggregates[J].Critical Reviewsin Plant Sciences,2004,23(6):481-504

[85]ZHONGZK,HANXH,XUYD,etal.Effects of landuse change on organic carbon dynamicsassociated with soil aggregate fractions on the Loess Plateau,China[J]. Land Degradation andDevelopment,2019,30(9):1070-1082

[86]ZHONG ZK,WUSJ,LUXQ,etal.Organic carbon,nitro gen accumulation,and soil aggregate dynamics as affected by vegetation restoration patterns in the Loess Plateau of China [J].Catena,2021,196(1):1-9

[87]馬志良,趙文強.植物群落向土壤有機碳輸入及其對氣候變 暖的響應研究進展[J].生態學雜志,2020,39(1:270-281

[88]彭琴,齊玉春,董云社,等.干旱半干旱地區草地碳循環關鍵 過程對降雨變化的響應[J].地理科學進展,2012,31(11): 1510-1518

(責任編輯 彭露茜)

猜你喜歡
施氮土壤有機氮量
施氮對棉花生長發育、產量及棉田土壤養分的影響
施氮改變了三江平原濕地小葉章對 CO2 升高的生理響應
林業科技(2025年4期)2025-08-18 00:00:00
氮用量和基追肥氮比例對稻茬烤煙氮、磷、鉀積累和產量、產值的影響
秸稈直接還田及炭化還田對土壤活性有機碳組分的影響綜述
農田管理對微生物源土壤有機碳轉化與穩定性影響的研究進展
施氮水平對植絲式混合草坪生長及氮磷養分去向的影響
草地學報(2025年7期)2025-08-03 00:00:00
不同施氮量對大寧曬煙生長發育和品質的影響
主站蜘蛛池模板: 麻豆精品在线| 米奇精品一区二区三区| 国产亚洲精品yxsp| 女人一级毛片| 波多野结衣AV无码久久一区| 久久网欧美| 亚洲av无码专区久久蜜芽| 欧美中文字幕一区| 一区二区三区四区日韩| 国产女人在线观看| 国产精品国产三级国产专业不| 亚洲无码免费黄色网址| 亚洲丝袜第一页| 都市激情亚洲综合久久| 91青青草视频在线观看的| 国产自无码视频在线观看| 自拍偷拍欧美日韩| 青青操国产| 日韩精品亚洲一区中文字幕| 91麻豆精品国产高清在线| 国产精品精品视频| 在线欧美国产| 在线日本国产成人免费的| 欧美不卡二区| 国产粉嫩粉嫩的18在线播放91| 韩日免费小视频| 在线观看无码a∨| 亚洲一区二区视频在线观看| 婷婷五月在线视频| 亚洲精品爱草草视频在线| 国产v欧美v日韩v综合精品| 四虎AV麻豆| a毛片在线免费观看| 91视频国产高清| 久久综合色播五月男人的天堂| 一级做a爰片久久免费| www中文字幕在线观看| 凹凸精品免费精品视频| 91小视频在线观看免费版高清| 国产免费网址| 亚洲三级色| 全午夜免费一级毛片| 日韩精品欧美国产在线| 欧美在线导航| 国产91透明丝袜美腿在线| 国产精品不卡永久免费| 精品三级网站| 欧美日韩高清| 国产精品无码作爱| 亚洲伦理一区二区| 一级毛片在线免费视频| 看看一级毛片| 欧美精品影院| 日本高清免费一本在线观看 | 5555国产在线观看| 亚洲天堂首页| 色悠久久综合| a毛片在线免费观看| 波多野结衣爽到高潮漏水大喷| 婷婷色在线视频| 国产精品人成在线播放| 国产情精品嫩草影院88av| 亚洲AV一二三区无码AV蜜桃| 国产黄在线免费观看| 成年人视频一区二区| 国产全黄a一级毛片| 欧美一级高清片久久99| 日韩av无码DVD| 亚洲女同一区二区| 97国产在线播放| 国产午夜看片| 伊人网址在线| 亚洲Av综合日韩精品久久久| 亚洲永久精品ww47国产| 国内毛片视频| 亚洲美女视频一区| 国产欧美性爱网| 欧洲在线免费视频| 中日韩欧亚无码视频| 波多野结衣久久高清免费| 精品丝袜美腿国产一区| 中文字幕首页系列人妻|