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

    東北黑土不同開墾年限稻田土壤有機(jī)氮礦化特征

    2022-05-16 05:17:30高佳蕊方勝志張玉玲安晶虞娜鄒洪濤
    關(guān)鍵詞:氮礦化黑土土壤有機(jī)

    高佳蕊,方勝志,張玉玲,安晶,虞娜,鄒洪濤

    東北黑土不同開墾年限稻田土壤有機(jī)氮礦化特征

    高佳蕊,方勝志,張玉玲*,安晶,虞娜,鄒洪濤

    沈陽農(nóng)業(yè)大學(xué)土地與環(huán)境學(xué)院/農(nóng)業(yè)農(nóng)村部東北耕地保育重點(diǎn)實(shí)驗(yàn)室,沈陽 110866

    【目的】分析東北黑土自然荒地開墾種稻后土壤礦化氮含量、氮凈礦化速率和氮凈礦化率,探討土壤供氮能力及其特點(diǎn),揭示土壤氮素的演變規(guī)律,為東北黑土的合理利用和培肥提供理論依據(jù)?!痉椒ā恳詵|北黑土自然荒地(0 a,為對(duì)照土壤,原始自然草甸植被)和不同開墾年限(12、35、62和85 a)的稻田(地形、種植制度、施肥和水分管理等大致相同)土壤為研究對(duì)象,采用長(zhǎng)期淹水密閉-間歇淋洗培養(yǎng)方法,研究黑土自然荒地開墾種稻后土壤有機(jī)氮礦化的特征?!窘Y(jié)果】在培養(yǎng)初期(約1個(gè)月),各年限土壤累積礦化氮量迅速增加,之后呈緩慢增加趨勢(shì);在淹水培養(yǎng)結(jié)束(297 d)時(shí),土壤累積礦化氮量為212.43—388.11 mg·kg-1,各開墾年限土壤累積礦化氮量大小順序?yàn)?、12、35、85和62 a。土壤有機(jī)氮礦化過程可用混合模型(Special 模型)很好地進(jìn)行描述,并可將土壤有機(jī)氮庫分為增量氮庫和其他組分氮庫,與對(duì)照土壤(0 a)相比,各年限稻田土壤增量氮庫的氮礦化勢(shì)(N)均呈下降趨勢(shì),其中,開墾12 a與開墾35、62和85 a稻田土壤的N之間均無顯著差異(>0.05),但開墾12和35 a稻田土壤的N顯著高于開墾62和85 a稻田土壤(<0.05),而礦化速率常數(shù)(k)均呈上升趨勢(shì),但各年限田土壤的k之間均無顯著差異(>0.05);開墾12和35 a稻田土壤其他組分氮庫的礦化速率常數(shù)(0)無顯著變化(>0.05),但開墾62和85 a稻田土壤的0則顯著下降(<0.05)。各年限土壤氮凈礦化速率在培養(yǎng)4 d時(shí)為最大,之后逐漸下降,在淹水培養(yǎng)結(jié)束(297 d)時(shí),土壤氮凈礦化速率大小順序與其累積礦化氮量的大小順序相一致;各年限土壤氮凈礦化率在淹水培養(yǎng)初期較高,之后緩慢增加,在培養(yǎng)結(jié)束(297 d)時(shí),土壤氮凈礦化率為78.60—101.82 mg·g-1,各開墾年限土壤氮凈礦化率的大小順序?yàn)? a、35 a、12 a、85和62 a;土壤全氮和C/N是影響土壤礦化氮量和氮凈礦化速率的重要因素(<0.05)。土壤初始礦質(zhì)氮與N之和可用來表征當(dāng)季稻田土壤供氮能力大小,與對(duì)照土壤(0 a)相比,各年限稻田土壤的供氮能力顯著下降(<0.05),開墾12和35 a稻田土壤的供氮能力顯著高于開墾62和85 a的稻田土壤(<0.05)?!窘Y(jié)論】黑土自然荒地開墾種稻85 a間,稻田土壤的供氮能力均有所下降,種稻大于35 a時(shí)下降顯著,因此在稻田土壤地力培育中應(yīng)注重土壤有機(jī)質(zhì)含量的提高。

    黑土區(qū);稻田土壤;種稻年限;礦化氮;土壤供氮能力

    0 引言

    【研究意義】東北黑土是我國(guó)重要的商品糧基地,約占全國(guó)糧食總產(chǎn)量的21%[1]。東北黑土的結(jié)構(gòu)良好,有機(jī)質(zhì)含量高,具有豐富的土壤養(yǎng)分及良好的物理性質(zhì)[2]。土壤碳氮是衡量土壤肥力的重要指標(biāo)[3],是維系土壤質(zhì)量的關(guān)鍵因素,對(duì)維持養(yǎng)分循環(huán)至關(guān)重要[4-5]。水稻在生長(zhǎng)發(fā)育過程中所吸收的氮素有50%—80%來源于土壤[6];土壤中90%以上氮素是以有機(jī)態(tài)氮的形式存在[7-8],有機(jī)態(tài)氮需礦化轉(zhuǎn)變成無機(jī)態(tài)氮(銨態(tài)氮或硝態(tài)氮)才能被作物吸收利用,土壤有機(jī)氮礦化過程生成的礦化氮(無機(jī)氮)是水稻生長(zhǎng)發(fā)育過程中最主要的氮素來源[9],因此土壤有機(jī)氮礦化能力可表征土壤的供氮潛力[10]。土壤含水量、溫度、土壤pH[11-13]等因素影響土壤微生物的生活環(huán)境,進(jìn)而影響土壤有機(jī)氮礦化過程及其礦化能力[14]。不同開墾年限稻田土壤的性質(zhì)及環(huán)境條件會(huì)發(fā)生明顯的變化,這可能會(huì)影響土壤有機(jī)氮的礦化能力,因此開展東北黑土不同開墾年限稻田土壤有機(jī)氮的礦化特征研究,對(duì)于明確黑土稻田土壤的有機(jī)氮轉(zhuǎn)化具有重要意義。【前人研究進(jìn)展】黑土開墾后土壤性質(zhì)發(fā)生了很大的變化,隨著種稻年限的延長(zhǎng),土壤有機(jī)碳、全氮含量和土壤質(zhì)地[15-18]發(fā)生了很大的變化。黑土開墾種稻85年間,隨著種稻年限的延長(zhǎng),土壤有機(jī)碳含量呈先增加后下降趨勢(shì),而全氮含量大體上呈現(xiàn)下降趨勢(shì)[19]。東北黑土區(qū)不同有機(jī)碳水平水稻土氮礦化勢(shì)及礦化速率常數(shù)具有明顯差異,土壤有機(jī)碳(氮)、C/N和pH是影響土壤有機(jī)氮礦化的重要因素[20-21],并認(rèn)為基于一級(jí)反應(yīng)理論的混合模型(Special模型)較一級(jí)模型(One-pool模型)和雙組分模型(Two-pool模型)能夠更好地描述黑土區(qū)水稻土有機(jī)氮的礦化過程[21]。土壤的供氮能力不僅包括土壤初始礦質(zhì)氮(無機(jī)氮,主要為NH4+-N 和NO3--N),而且也包括當(dāng)季作物生長(zhǎng)發(fā)育過程中土壤有機(jī)氮礦化生成的礦化氮(無機(jī)氮,主要為NH4+-N和NO3--N)[22]。土壤有機(jī)氮礦化過程生成的礦化氮是土壤礦質(zhì)氮的主要來源,土壤氮凈礦化速率是氮初級(jí)礦化速率和初級(jí)固定速率的綜合結(jié)果,可以用來表征土壤中無機(jī)氮供應(yīng)能力的大小[23-24]。李平等[25]研究表明,開墾2年土壤氮凈礦化速率顯著低于開墾30年土壤。顧春朝[26]研究發(fā)現(xiàn),80余年稻田土壤氮礦化勢(shì)和氮凈礦化速率顯著高于120余年土壤。由此可見,不同開墾年限對(duì)稻田土壤全氮及有機(jī)氮的礦化能力影響的研究結(jié)果不盡相同?!颈狙芯壳腥朦c(diǎn)】目前關(guān)于東北黑土不同有機(jī)碳水平水稻土有機(jī)碳礦化特征的研究有報(bào)道,但關(guān)于東北黑土不同開墾年限稻田土壤有機(jī)氮的礦化特征的研究比較薄弱,需進(jìn)一步深入探討?!緮M解決的關(guān)鍵問題】本文采用長(zhǎng)期淹水密閉培養(yǎng)-間歇淋洗法室內(nèi)培養(yǎng)試驗(yàn),研究東北黑土區(qū)不同開墾年限(0、12、35、62和85 a)稻田土壤有機(jī)氮的礦化特征,以期明確黑土區(qū)開墾種稻后土壤的供氮能力及供氮特點(diǎn),為該地區(qū)黑土的合理利用及稻田土壤的合理培肥提供理論依據(jù)。

    1 材料與方法

    1.1 供試土壤

    東北黑土區(qū)(122—132° E,43—50° N)主要分布在松嫩平原山前波狀起伏臺(tái)地和漫崗丘陵區(qū)。本研究土壤樣本采自黑龍江省綏化市慶安縣,該研究區(qū)域地處中國(guó)黑龍江省中部,小興安嶺南麓向松嫩平原的過渡地帶,地理位置為127°30'—128°35' E,46°30'—47°35' N,氣候?qū)俸疁貛Т箨懶约撅L(fēng)氣候,溫?zé)釢駶?rùn);年平均日照時(shí)數(shù)為2 599 h,年平均氣溫1.69 ℃,無霜期128 d,年平均降水量577 mm,研究區(qū)黑土均為黃土性沉積物發(fā)育母質(zhì)[21]。項(xiàng)目組于2015年10月通過詳細(xì)的實(shí)地調(diào)查,確定慶安縣勤勞鎮(zhèn)張家爐村為供試土壤的采集區(qū)域,在研究區(qū)內(nèi)選取典型黑土未開墾利用的荒地(0年,作為對(duì)照土壤,以自然草甸植被為主,坡度小于3°)和不同開墾種稻年限(12、35、62和85a)的種稻田塊;其中,不同年限稻田田塊的氣候和地形大致相同,在開墾種稻前植被與未開墾利用的自然草甸植被狀況相近,各年限種稻田塊每年均未施有機(jī)肥,每年施用氮、磷、鉀化肥用量(農(nóng)戶常規(guī)施用N 90—120 kg·hm-2、P2O545—60 kg·hm-2、K2O 45—75 kg·hm-2)和水分管理(農(nóng)戶常規(guī)灌溉的灌溉定額5 500—6 500 m3·hm-2)措施基本相同,但年限間存在差異;各年限種稻田塊每年秋季水稻秸稈移除,水稻殘茬隨深翻入田。土壤樣本采集時(shí),將每個(gè)相同年限的田塊作為一個(gè)采樣區(qū)域,并確定面積大致相同的 3 個(gè)采樣單元(0.13—0.2 hm2),每個(gè)采樣單元“S”形布點(diǎn),每點(diǎn)采樣深度為0—20 cm,采集5—7點(diǎn)均勻混合后作為一個(gè)土壤樣本,即每個(gè)年限的土壤樣本均為3次重復(fù)樣本。供試土壤的基本性質(zhì)見表1。

    表1 供試土壤基本性質(zhì)

    數(shù)據(jù)后不同小寫字母表示不同年限間差異達(dá)0.05顯著水平

    Different lowercase letters indicate the significance difference at 0.05 level between different years

    1.2 土壤有機(jī)氮礦化培養(yǎng)試驗(yàn)

    土壤有機(jī)氮礦化培養(yǎng)試驗(yàn)采用長(zhǎng)期淹水密閉培養(yǎng)-間歇淋洗法進(jìn)行[21,27]。稱取過2 mm篩的風(fēng)干土15.0 g,置于100 mL離心管中,每管加過 2 mm篩的石英砂 7.5 g后均勻混合,每一年限樣本 3 次重復(fù)。首先,初始礦質(zhì)氮(無機(jī)氮,主要為NH4+-N 和 NO3--N)淋洗:向管中加入 30 mL蒸餾水,蓋緊離心管蓋,淋洗土樣 1 次、再向管中加入 1 mol·L-1KCl 30 mL淋洗土樣 2 次、最后向管中加 30 mL蒸餾水淋洗土樣兩次;每次淋洗時(shí),加入蒸餾水或淋洗液后攪拌均勻,以 6 000 r/min轉(zhuǎn)速離心 10 min,每次淋洗后將上清液傾入 200 mL容量瓶,最后將淋洗液定容至 200 mL。其次,礦化氮淋洗:將初始淋洗后樣品加入 30 mL蒸餾水后攪拌均勻,置于 30 ℃恒溫箱中培養(yǎng)(培養(yǎng)溫度以土壤樣本供試區(qū)域水稻生長(zhǎng)旺盛時(shí)期的平均最高氣溫進(jìn)行設(shè)置),分別于培養(yǎng)的第 4、7、14、21、28、42、56、70、84、122、241、297天全部取出淋洗,淋洗方法與初始礦質(zhì)氮淋洗相同,并定容至 200 mL。初始淋洗液中測(cè)定礦質(zhì)氮(主要為NH4+-N 和 NO3--N),礦化淋洗液中測(cè)定銨態(tài)氮(NH4+-N)。

    采用非線性擬合建立實(shí)測(cè)累積礦化氮量與培養(yǎng)時(shí)間關(guān)系的混合模型(Special 模型)[21,28],Special 模型為:

    N=N[1 – exp(-kt)] +0

    式中,N為培養(yǎng)時(shí)間時(shí)的累積礦化氮量(mg·kg-1);Nk分別為培養(yǎng)時(shí)間趨于無限長(zhǎng)時(shí)增量氮庫(即土樣風(fēng)干過程中死亡微生物體氮)的氮礦化勢(shì)(mg·kg-1)和礦化速率常數(shù)(d-1);0為其他組分氮庫的礦化速率常數(shù)(mg·kg-1·d-1);為礦化培養(yǎng)時(shí)間(d)[21,29]。

    根據(jù)時(shí)間實(shí)測(cè)的累積礦化氮量,進(jìn)一步求得土壤氮凈礦化速率和氮凈礦化率;土壤氮凈礦化速率(mg·kg-1·d-1)表示單位培養(yǎng)時(shí)間(d)內(nèi)土壤礦化生成的礦化氮(NH4+-N)數(shù)量(mg·kg-1),土壤氮凈礦化率(mg·g-1)表示一段培養(yǎng)時(shí)間內(nèi)1 g土壤全氮礦化生成的礦化氮(NH4+-N)數(shù)量(mg),即一段培養(yǎng)時(shí)間內(nèi)土壤凈礦化氮量(mg·kg-1)與土壤全氮(g·kg-1)的比值。

    1.3 測(cè)定方法

    土壤有機(jī)碳、全氮采用元素分析儀測(cè)定(德國(guó),Elementar 公司),淋洗液中 NH4+-N、NO3--N采用 AA3 自動(dòng)分析儀(Seal Analytical USA)測(cè)定。

    1.4 統(tǒng)計(jì)分析

    利用Origin 2017軟件進(jìn)行作圖,Excel 2016 和 SPSS window version 19.0 軟件進(jìn)行數(shù)據(jù)處理和統(tǒng)計(jì)分析,LSD法進(jìn)行多重比較。文中數(shù)據(jù)為3次重復(fù)的平均值。

    2 結(jié)果

    2.1 土壤累積礦化氮量—時(shí)間曲線

    土壤有機(jī)氮礦化過程是有機(jī)態(tài)氮經(jīng)微生物作用轉(zhuǎn)化生成為銨態(tài)氮和硝態(tài)氮的過程;淹水培養(yǎng)過程中,淋洗液中僅檢測(cè)出銨態(tài)氮。由圖1可見,各年限土壤累積礦化氮量隨時(shí)間變化的趨勢(shì)大致相同,在淹水培養(yǎng)初期(約1個(gè)月)土壤有機(jī)氮礦化速率較快,累積礦化氮量迅速增加,之后隨著培養(yǎng)時(shí)間的延長(zhǎng)增加緩慢,淹水培養(yǎng)297 d時(shí)各年限土壤累積礦化氮量為212.43—388.11 mg·kg-1;在淹水培養(yǎng)的42 d內(nèi),開墾12和35 a的稻田土壤累積礦化氮量要高于對(duì)照土壤(0 a),在淹水培養(yǎng)42 d后,則明顯低于對(duì)照土壤(0 a),但在整個(gè)培養(yǎng)期間,開墾12和35 a的稻田土壤與對(duì)照土壤相比,土壤累積礦化氮量無顯著差異(>0.05);在培養(yǎng)1個(gè)月后,開墾62和85 a的稻田土壤累積礦化氮量均明顯低于對(duì)照土壤(0 a)和開墾12和35 a的稻田土壤??傮w上來看,在淹水培養(yǎng)42 d 后,各開墾年限土壤累積礦化氮量大小順序?yàn)?、12、35、85和62 a,其中,開墾12和35 a的稻田土壤累積礦化氮量下降不顯著(>0.05),而開墾62和85 a稻田土壤累積礦化氮量則顯著下降(<0.05),但二者之間無顯著差異(>0.05)。相關(guān)分析顯示,淹水培養(yǎng)297 d時(shí)土壤累積礦化氮量與其土壤全氮之間具有顯著線性正相關(guān)(= 0.955,<0.05),而與土壤C/N之間則呈顯著曲線負(fù)相關(guān)(= -0.927,<0.05),說明土壤全氮含量及C/N是影響土壤有機(jī)氮礦化的重要原因。

    圖中曲線為利用Special模型擬合的曲線;不同小寫字母表示相同取樣時(shí)間、不同年限差異達(dá) 0.05 顯著水平

    2.2 Special模型對(duì)土壤有機(jī)氮礦化過程的擬合參數(shù)

    利用Special模型對(duì)實(shí)測(cè)累積礦化氮量進(jìn)行擬合,其擬合參數(shù)見表2。從擬合方程的決定系數(shù)(2)和估計(jì)標(biāo)準(zhǔn)誤差()及擬合度(圖1)綜合來看,Special模型均能很好地描述各年限土壤有機(jī)氮的礦化過程。各年限土壤增量氮庫的氮礦化勢(shì)(N)為146.78— 254.50 mg·kg-1,礦化速率常數(shù)(k)為0.042—0.119 d-1;與對(duì)照土壤(0 a)相比,各年限稻田土壤的N均呈現(xiàn)下降趨勢(shì),其中,開墾12和35 a、62和85 a稻田土壤的N之間均無顯著差異(>0.05),但開墾12、35 a稻田土壤的N顯著高于開墾62 a,開墾35 a稻田土壤的N顯著高于開墾85 a(<0.05);各開墾年限稻田土壤的k均顯著高于對(duì)照土壤(0 a),但各年限稻田間土壤的k均無顯著差異(>0.05)。NN(297 d)的累積礦化氮量之間具有顯著的正相關(guān)(=0.974,<0.01),N與土壤C/N之間呈顯著曲線負(fù)相關(guān)(=-0.957,<0.05),k則與土壤C/N之間呈顯著曲線正相關(guān)(=0.925,<0.05),說明土壤C/N是影響土壤增量氮庫礦化的重要因素。

    表2 Special模型的擬合參數(shù)

    k分別為土樣增量氮庫的氮礦化勢(shì)和氮礦化速率常數(shù),0為其他組分氮庫的氮礦化速率常數(shù);2和分別為方程擬合的決定系數(shù)和估計(jì)標(biāo)準(zhǔn)誤差;同列不同小寫字母表示差異達(dá) 0.05 顯著水平;**表示方程擬合達(dá) 0.01 顯著水平

    andkare the potentially mineralizable nitrogen (N) and rate constant of mineralization of the increment N pool, respectively,kis the rate constant of mineralization of the resistant N pool;2andare the determination coefficient and the standard error of estimation of equation fitting, respectively; Different lowercase letters indicate the significantly difference at 0.05 levelamong different years; ** indicates that the equation fitting reaches a significant at 0.01 level

    各年限土壤其他組分氮庫的礦化速率常數(shù)(0)為0.21—0.48 mg·kg-1·d-1,稻田土壤的0隨種稻年限延長(zhǎng)呈下降趨勢(shì);其中,開墾12和35 a稻田土壤的0與對(duì)照土壤(0 a)之間無顯著差異(>0.05),但顯著大于開墾62和85 a的稻田土壤(<0.05),而開墾62與85 a稻田土壤的0之間則無顯著差異(>0.05)。0與土壤全氮、有機(jī)碳含量之間均具有線性正相關(guān)關(guān)系(=0.979,<0.01;=0.899,<0.05),但與土壤全氮的相關(guān)性更大,說明土壤有機(jī)碳、全氮含量影響著其他組分氮庫,但土壤全氮對(duì)其他組分氮庫的礦化影響更為明顯。

    2.3 土壤氮凈礦化速率和凈礦化率

    各年限土壤氮凈礦化速率曲線的變化趨勢(shì)基本相同(圖2-A)。在培養(yǎng)最初4 d時(shí)為最大(13.52—23.18 mg·kg-1·d-1),各開墾年限土壤氮凈礦化速率大小順序?yàn)?5、12、85、62和0 a;在培養(yǎng)約1個(gè)月時(shí),土壤氮凈礦化速率快速下降,隨后呈緩慢下降趨勢(shì),在培養(yǎng)結(jié)束的297 d時(shí)數(shù)值較小且趨于穩(wěn)定(0.72—1.31 mg·kg-1·d-1),各開墾年限土壤氮凈礦化速率的大小順序則表現(xiàn)為0、12、35、85和62 a,但各年限間土壤氮凈礦化速率差異較小。在培養(yǎng)結(jié)束(297 d)時(shí),各年限土壤氮凈礦化速率與土壤全氮呈顯著曲線正相關(guān)關(guān)系(=0.957,<0.05),與土壤C/N呈顯著曲線負(fù)相關(guān)(=-0.931,<0.05),這也說明土壤氮凈礦化速率大小也受土壤全氮和C/N的影響。

    土壤氮凈礦化率為一段時(shí)間內(nèi)土壤中的礦化氮量與土壤全氮的比值,其數(shù)值大小表征單位數(shù)量土壤全氮中含有可礦化氮的數(shù)量或單位數(shù)量土壤全氮礦化生成的礦化氮的數(shù)量,其值越大意味著土壤中可礦化氮數(shù)量越多或礦化生成的礦化氮數(shù)量越多。各開墾年限土壤氮凈礦化率隨時(shí)間變化的趨勢(shì)也大致相同(圖2-B)。在培養(yǎng)42 d前,各年限稻田土壤氮凈礦化率均要高于對(duì)照土壤(0 a),表現(xiàn)為開墾85 a>開墾35 a>開墾12 a>開墾62 a>自然荒地(0 a);在培養(yǎng)42 d后,各年限土壤氮凈礦化率逐漸增加,在培養(yǎng)至120 d左右,稻田土壤氮凈礦化率均低于對(duì)照土壤,在培養(yǎng)結(jié)束(297 d)時(shí),各年限土壤氮凈礦化率為78.60—101.82 mg·g-1,表現(xiàn)為自然荒地(0 a)>開墾35 a>開墾12 a>開墾85 a>開墾62 a,這說明在淹水培養(yǎng)前期各年限稻田土壤可礦化有機(jī)氮數(shù)量較多,這些可礦化有機(jī)氮經(jīng)過礦化生成銨態(tài)氮成為作物吸收利用的氮源。

    圖2 土壤氮凈礦化速率和凈礦化率

    2.4 土壤的供氮能力

    各年限土壤初始礦質(zhì)氮含量為42.82—137.92 mg·kg-1,除開墾35 a外,其他年限稻田土壤初始礦質(zhì)氮均顯著低于對(duì)照土壤(<0.05)(表1)。東北黑土區(qū)水稻生長(zhǎng)期約為130 d,本試驗(yàn)培養(yǎng)時(shí)間為297 d,遠(yuǎn)超過了水稻當(dāng)季生長(zhǎng)周期,因此本試驗(yàn)以Special模型擬合曲線計(jì)算求得各年限土壤培養(yǎng)130 d時(shí)的礦化氮量,即開墾0、12、35、62和85 a土壤在當(dāng)季水稻生長(zhǎng)過程中的礦化氮量分別301.47、261.38、262.86、172.19和197.00 mg·kg-1,其數(shù)值與各年限土壤增量氮庫的氮礦化勢(shì)相接近,且二者之間具有顯著線性正相關(guān)關(guān)系(=0.988,<0.01)。因此可以土壤初始礦質(zhì)氮與增量氮庫的氮礦化勢(shì)(N)之和來表征土壤的供氮能力(圖3)。與對(duì)照土壤(0 a)相比,各年限稻田土壤的供氮能力均顯著下降(<0.05),降低幅度分別為21.22%、12.61%、47.19%和44.71%,但開墾12與35 a、62與85 a稻田土壤的供氮能力之間均無顯著差異(>0.05)。另外,各年限土壤初始礦質(zhì)氮占土壤供氮能力的比例分別34.66%、31.98%、36.76%、29.05%和19.81%。由此可見,黑土自然荒地開墾種稻年限小于35 a時(shí),土壤的供氮能力較大,初始礦質(zhì)氮占比較大,而在種稻年限大于35 a時(shí),稻田土壤的供氮能力顯著下降,土壤初始礦質(zhì)氮占比也相對(duì)較小。

    柱上不同小寫字母表示差異達(dá)0.05顯著水平;NF為增量氮庫的氮礦化勢(shì);Nmin為初始礦質(zhì)氮

    3 討論

    3.1 開墾年限對(duì)黑土稻田土壤有機(jī)氮礦化特征的影響

    土壤有機(jī)氮礦化過程受土壤理化性質(zhì)的影響,土地利用方式、年限的不同,其理化性質(zhì)也會(huì)存在明顯的差異[15-19]。已研究表明,淹水培養(yǎng)161 d時(shí),土壤礦化氮數(shù)量與有機(jī)碳、全氮含量高度相關(guān)[20-21],土壤氮礦化勢(shì)與土壤C/N呈顯著正相關(guān),土壤氮礦化速率常數(shù)則與土壤C/N呈顯著負(fù)相關(guān)[20]。在本研究中,各年限土壤有機(jī)碳、全氮含量、C/N、初始礦質(zhì)氮等存在不同程度的差異(表1),在培養(yǎng)297 d時(shí)累積礦化氮量與土壤全氮含量呈顯著正相關(guān),與土壤C/N呈顯著負(fù)相關(guān),這與前人的研究結(jié)果相一致[20-21]。土壤增量氮庫的氮礦化勢(shì)(N)與土壤C/N呈顯著負(fù)相關(guān),而氮礦化速率常數(shù)(k)與土壤C/N呈顯著正相關(guān),其他組分氮庫的礦化速率常數(shù)(k)與土壤全氮、有機(jī)碳含量均呈顯著正相關(guān),說明有機(jī)碳、全氮含量及其C/N在影響土壤有機(jī)氮礦化數(shù)量或氮礦化潛勢(shì)方面具有非常重要的作用。

    氮凈礦化速率是氮初級(jí)礦化速率和初級(jí)固定速率的綜合結(jié)果[23-24],當(dāng)土壤中能源物質(zhì)缺乏或C/N較小時(shí),土壤氮初級(jí)礦化速率強(qiáng)于土壤生物的初級(jí)固定速率,此時(shí)土壤氮素轉(zhuǎn)化表現(xiàn)為土壤無機(jī)氮素的累積過程,即土壤氮凈礦化速率為正值;反之,當(dāng)土壤能源物質(zhì)充足或C/N較高時(shí),土壤氮素轉(zhuǎn)化主要表現(xiàn)為生物的固定,土壤氮素以消耗為主[30]。一般情況下土壤C/N越低,土壤氮素礦化速率一般較大[31]。自然黑土經(jīng)開墾種稻后,種稻62和85 a土壤全氮含量顯著降低,各年限稻田土壤C/N顯著增加(表1),說明不同種稻年限土壤有機(jī)碳、全氮含量發(fā)生了明顯變化,同時(shí)土壤有機(jī)質(zhì)的組成結(jié)構(gòu)、復(fù)雜程度也隨之發(fā)生了變化。本研究中,氮凈礦化速率與土壤C/N呈顯著負(fù)相關(guān)關(guān)系,與土壤全氮含量呈顯著正相關(guān)關(guān)系,因此,不同年限土壤氮凈礦化速率的差異可能與土壤全氮含量及C/N有關(guān)。開墾年限顯著影響稻田土壤氮凈礦化速率和NH4+-N含量[25-26]。已有研究表明,開墾年限對(duì)旱地土壤氮凈礦化速率影響顯著,好氧培養(yǎng)7 d時(shí)開墾30 a土壤氮凈礦化速率顯著高于開墾2 a土壤[25],開墾種稻年限越長(zhǎng),越不利于土壤NH4+-N的積累[26]。在淹水培養(yǎng)297 d時(shí),各年限土壤累積礦化氮(NH4+-N)量也具有顯著差異,種稻時(shí)間為62和85 a土壤累積礦化氮(NH4+-N)量顯著低于對(duì)照土壤和其他兩個(gè)開墾年限稻田土壤(圖1),說明自然黑土經(jīng)開墾種稻時(shí)間較長(zhǎng)時(shí)土壤中NH4+-N的累積數(shù)量有所下降。本研究中,在淹水培養(yǎng)4 d時(shí),各年限稻田土壤氮凈礦化速率均大于對(duì)照土壤,而在淹水培養(yǎng)120 d后,各年限稻田土壤氮凈礦化速率均低于對(duì)照土壤(圖2-a),這說明自然黑土開墾種稻后,短時(shí)間淹水可使稻田土壤氮凈礦化速率增大,礦化氮數(shù)量增加,但隨著淹水時(shí)間的延長(zhǎng),稻田土壤氮凈礦化速率減小,礦化氮數(shù)量降低,淹水時(shí)間長(zhǎng)短也是影響稻田土壤氮凈礦化速率大小的重要因素。

    土壤氮庫可分為增量氮庫和其他組分氮庫,其中,增量氮庫是指土樣風(fēng)干過程中死亡微生物體氮庫[21,29],增量氮庫的氮礦化勢(shì)(N)是在一定條件下土壤中這一部分氮庫能夠礦化成無機(jī)態(tài)氮的數(shù)量的最大值,可以用來表征土壤中這一部分氮庫的供氮容量指標(biāo),增量氮庫的礦化速率常數(shù)(k)可用來表征這一部分氮庫礦化速率快慢的供氮強(qiáng)度指標(biāo),而其他組分氮庫的礦化速率常數(shù)(0)則可用來表示土壤中這一部分氮庫的礦化速率快慢的供氮強(qiáng)度指標(biāo)。研究表明,土壤有機(jī)碳是影響氮礦化勢(shì)的主要因素[32-33],土壤有機(jī)碳可以為微生物的生存提供養(yǎng)分和能量,有機(jī)碳含量越高,微生物活性越大,土壤風(fēng)干過程中死亡微生物體氮庫越多,氮礦化勢(shì)也會(huì)越大[34-35]。本研究中,與對(duì)照土壤(0 a)相比,各年限稻田土壤增量氮庫的氮礦化勢(shì)(N)均呈現(xiàn)逐漸降低的趨勢(shì),開墾后種稻12和35 a的稻田土壤的氮礦化勢(shì)(N)顯著高于開墾62與85 a的稻田土壤(表2),這與土壤有機(jī)碳的變化趨勢(shì)不完全一致(表 1),各年限稻田土壤增量氮庫的礦化速率常數(shù)(k)顯著高于對(duì)照土壤,但各年限之間k則無顯著差異(表2),Nk兩個(gè)參數(shù)與土壤C/N關(guān)系密切。從土壤其他組分氮庫的礦化速度常數(shù)(k)來看,開墾后種稻62和85 a稻田土壤的0則顯著低于其他年限土壤(表2),0與土壤有機(jī)碳、全氮含量關(guān)系密切。這一結(jié)果說明自然黑土開墾種稻后土壤有機(jī)碳、全氮含量及其結(jié)構(gòu)等發(fā)生變化,一方面影響了稻田土壤增量氮庫,種稻年限越長(zhǎng),土壤增量氮庫的供氮容量減少、供氮強(qiáng)度增大;另一方面也影響了其他組分氮庫,種稻年限越長(zhǎng),土壤其他組分氮庫的供氮容量和強(qiáng)度均有所下降。但土壤增量氮庫的供氮潛力要明顯大于其他組分氮庫的供氮潛力。另外,土壤有機(jī)氮的礦化數(shù)量及速率與土壤溫度有關(guān),本研究中培養(yǎng)試驗(yàn)溫度(30 ℃)設(shè)置以供試區(qū)域水稻生長(zhǎng)旺盛時(shí)期平均最高氣溫進(jìn)行設(shè)置,這一溫度設(shè)置在一定程度上增加了水稻生長(zhǎng)非旺盛時(shí)期土壤有機(jī)氮的礦化數(shù)量和氮凈礦化速率,進(jìn)而也會(huì)高估各年限土壤的供氮潛力。

    3.2 開墾年限對(duì)黑土稻田土壤供氮能力的影響

    土壤無機(jī)態(tài)氮(以NH4+-N和NO3--N為主)可直接被作物吸收,是表征土壤供氮能力的重要指標(biāo)[36],也是國(guó)內(nèi)外研究土壤氮素有效性的熱點(diǎn)之一[37-38]。土壤氮凈礦化率是衡量氮素有效性的重要指標(biāo),氮凈礦化率是指在一定時(shí)間內(nèi)土壤礦化氮數(shù)量與全氮的比例,其比值的意義為1 g全氮可礦化的氮的毫克數(shù),其數(shù)值大小在一定程度上可以反映土壤有機(jī)氮的品質(zhì),其值越大說明可礦化的有機(jī)氮數(shù)量越多[39]。本研究中,各年限土壤初始礦質(zhì)氮含量(主要以NH4+-N含量為主)具有顯著差異,稻田土壤初始礦質(zhì)氮及其占全氮的比例均低于對(duì)照土壤(0 a)(表1),說明自然黑土開墾種稻后土壤中可供給當(dāng)季水稻吸收利用的無機(jī)態(tài)氮數(shù)量有所減少,種稻年限較長(zhǎng)時(shí)減少更為明顯。各年限土壤增量氮庫的氮礦化勢(shì)(N)占全氮的比例各不相同,開墾0、12、35、62和85 a土壤的N占全氮的比例分別為66.83、56.88、65.47、53.60和69.11 mg·g-1,分別占培養(yǎng)297 d時(shí)氮凈礦化率的65.64%、61.52%、65.84%、68.19%和74.89%,這一結(jié)果說明各年限土壤中增量氮庫和其他組分氮庫在土壤有機(jī)氮礦化過程中都發(fā)揮著作用,共同決定著一定時(shí)間內(nèi)土壤礦化氮數(shù)量的多少,但土壤增量氮庫的供氮潛力要明顯大于其他組分氮庫的供氮潛力。

    東北黑土區(qū)水稻種植制度為一年一熟制,土壤的供氮能力大小主要體現(xiàn)在水稻生長(zhǎng)季土壤能夠提供無機(jī)氮的數(shù)量,既包括土壤初始質(zhì)氮,也包括水稻生長(zhǎng)季土壤有機(jī)氮礦化生成的礦化氮[22]。本研究結(jié)果顯示,除了各年限土壤初始礦質(zhì)氮外,土壤增量氮庫的供氮潛力是衡量各年限土壤供氮能力大小的重要指標(biāo),各年限稻田土壤的供氮能力顯著低于對(duì)照土壤,且隨種稻年限延長(zhǎng)呈下降趨勢(shì),且種稻年限大于35 a稻田土壤的供氮能力顯著低于種稻年限小于35 a的稻田土壤(圖3)。這一結(jié)果顯示,黑土自然荒地開墾種稻一定年限(大于35 a)后,不僅土壤有機(jī)碳(氮)水平有所下降,而且土壤的供氮能力也會(huì)隨之下降,因此在東北黑土區(qū)稻田土壤地力培育中應(yīng)注意土壤有機(jī)質(zhì)含量的提高。

    4 結(jié)論

    開墾年限對(duì)東北黑土區(qū)稻田土壤有機(jī)氮的礦化量、氮凈礦化速率和氮凈礦化率以及供氮能力具有顯著的影響。淹水培養(yǎng)297 d,經(jīng)歷不同種植年限稻田土壤的礦化氮量和氮凈礦化速率均顯著低于自然荒地土壤。土壤初始礦質(zhì)氮和土壤氮增量氮庫的氮礦化勢(shì)之和可表征土壤的供氮能力,且不同年限稻田土壤的供氮能力均顯著低于自然荒地土壤,種稻大于35 a稻田土壤的供氮能力顯著低于種稻小于35 a的稻田土壤。

    [1] 國(guó)家統(tǒng)計(jì)局.2018 中國(guó)統(tǒng)計(jì)年鑒.北京: 中國(guó)統(tǒng)計(jì)出版社, 2018

    National Bureau of Statistics.China statistical yearbook.Beijing: China Statistics Press, 2018.(in Chinese)

    [2] 康日峰, 任意, 吳會(huì)軍, 張淑香.26年來東北黑土區(qū)土壤養(yǎng)分演變特征.中國(guó)農(nóng)業(yè)科學(xué), 2016, 49(11): 2113-2125.doi:10.3864/j.issn.0578-1752.2016.11.008.

    KANG R F, REN Y, WU H J, ZHANG S X.Changes in the nutrients and fertility of black soil over 26 years in northeast China.Scientia Agricultura Sinica, 2016, 49(11): 2113-2125.doi:10.3864/j.issn.0578- 1752.2016.11.008.(in Chinese)

    [3] HUANG Y, SUN W J, ZHANG W, YU Y Q.Changes in soil organic carbon of terrestrial ecosystems in China: a mini-review.Science China Life Sciences, 2010, 53(7): 766-775.doi:10.1007/s11427- 010-4022-4.

    [4] MASEK J G, HAYES D J, JOSEPH HUGHES M, HEALEY S P, TURNER D P.The role of remote sensing in process-scaling studies of managed forest ecosystems.Forest Ecology and Management, 2015, 355: 109-123.doi:10.1016/j.foreco.2015.05.032.

    [5] DENG L, WANG G L, LIU G B, SHANGGUAN Z P.Effects of age and land-use changes on soil carbon and nitrogen sequestrations following cropland abandonment on the Loess Plateau, China.Ecological Engineering, 2016, 90: 105-112.doi:10.1016/j.ecoleng.2016.01.086.

    [6] FAN H Z, ZHANG H, FENG W Q, ZHANG J, WANG C T.Effect of transplanting density on rice yield, nitrogen uptake and15N-fertilizer fate.Agricultural Science & Technology, 2012, 13(5): 1037-1039, 1054.doi:10.16175/j.cnki.1009-4229.2012.05.036.

    [7] LI D J, LIU J, CHEN H, ZHENG L, WANG K L.Soil gross nitrogen transformations in responses to land use conversion in a subtropical Karst region.Journal of Environmental Management, 2018, 212: 1-7.doi:10.1016/j.jenvman.2018.01.084.

    [8] LUCE M S, WHALEN J K, ZIADI N, ZEBARTH B J.Nitrogen dynamics and indices to predict soil nitrogen supply in humid temperate soils.Advances in Agronomy.Amsterdam: Elsevier, 2011: 55-102.doi:10.1016/b978-0-12-385538-1.00002-0.

    [9] SHARIFI M, ZEBARTH B J, BURTON D L, GRANT C A, PORTER G A.Organic amendment history and crop rotation effects on soil nitrogen mineralization potential and soil nitrogen supply in a potato cropping system.Agronomy Journal, 2008, 100(6): 1562-1572.doi:10.2134/agronj2008.0053.

    [10] KAUR J, CIHACEK L J, CHATTERJEE A.Estimation of nitrogen and sulfur mineralization in soils amended with crop residues contributing to nitrogen and sulfur nutrition of crops in the North Central US.Communications in Soil Science and Plant Analysis, 2018, 49(18): 2256-2266.doi:10.1080/00103624.2018.1499761.

    [11] 唐海龍, 王景燕, 黃帥, 龔偉, 周于波.華西雨屏區(qū)常綠闊葉林土壤氮礦化對(duì)溫度和濕度變化的響應(yīng).甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào), 2019, 54(2): 124-131.doi:10.13432/j.cnki.jgsau.2019.02.017.

    TANG H L, WANG J Y, HUANG S, GONG W, ZHOU Y B.Responses of soil nitrogen mineralization of evergreen broad-leaved forest in rainy area of Western China to moisture and temperature.Journal of Gansu Agricultural University, 2019, 54(2): 124-131.doi:10.13432/j.cnki.jgsau.2019.02.017.(in Chinese)

    [12] SIERRA J.Nitrogen mineralization and nitrification in a tropical soil: effects of fluctuating temperature conditions.Soil Biology and Biochemistry, 2002, 34(9): 1219-1226.doi:10.1016/s0038-0717(02) 00058-5.

    [13] HANAN E J, SCHIMEL J P, DOWDY K, D'ANTONIO C M.Effects of substrate supply, pH, and char on net nitrogen mineralization and nitrification along a wildfire-structured age gradient in chaparral.Soil Biology and Biochemistry, 2016, 95: 87-99.doi:10.1016/j.soilbio.2015.12.017.

    [14] KUYPERS M M M, MARCHANT H K, KARTAL B.The microbial nitrogen-cycling network.Nature Reviews Microbiology, 2018, 16(5): 263-276.doi:10.1038/nrmicro.2018.9.

    [15] 于君寶, 劉景雙, 王金達(dá), 劉淑霞, 齊小寧, 王洋, 王國(guó)平.不同開墾年限黑土有機(jī)碳變化規(guī)律.水土保持學(xué)報(bào), 2004, 18(1): 27-30.doi:10.13870/j.cnki.stbcxb.2004.01.007.

    YU J B, LIU J S, WANG J D, LIU S X, QI X N, WANG Y, WANG G P.Organic carbon variation law of black soil during different tillage period.Journal of Soil Water Conservation, 2004, 18(1): 27-30.doi:10.13870/j.cnki.stbcxb.2004.01.007.(in Chinese)

    [16] 賈樹海, 張佳楠, 張玉玲, 黨秀麗, 范慶鋒, 王展, 虞娜, 鄒洪濤, 張玉龍.東北黑土區(qū)旱田改稻田后土壤有機(jī)碳、全氮的變化特征.中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(7): 1252-1262.doi:10.3864/j.issn.0578-1752.2017.07.008.

    JIA S H, ZHANG J N, ZHANG Y L, DANG X L, FAN Q F, WANG Z, YU N, ZOU H T, ZHANG Y L.Changes of the characteristics of soil organic carbon and total nitrogen after conversation from upland to paddy field in black soil region of northeast China.Scientia Agricultura Sinica, 2017, 50(7): 1252-1262.doi:10.3864/j.issn.0578- 1752.2017.07.008.(in Chinese)

    [17] 馬原, 遲美靜, 張玉玲, 范慶峰, 虞娜, 鄒洪濤.黑土旱地改稻田土壤水穩(wěn)性團(tuán)聚體有機(jī)碳和全氮的變化特征.中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(8): 1594-1605.doi:10.3864/j.issn.0578-1752.2020.08.009.

    MA Y, CHI M J, ZHANG Y L, FAN Q F, YU N, ZOU H T.Change characteristics of organic carbon and total nitrogen in water-stable aggregate after conversion from upland to paddy field in black soil.Scientia Agricultura Sinica, 2020, 53(8): 1594-1605.doi:10.3864/j.issn.0578-1752.2020.08.009.(in Chinese)

    [18] 谷思玉, 胡洋, 聶艷龍, 蔡海森, 宋秀麗, 何鑫, 崔凱歌.農(nóng)墾九三管理局不同開墾年限農(nóng)田黑土團(tuán)聚體變化.東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2015, 46(11): 56-62.doi:10.19720/j.cnki.issn.1005-9369.2015.11.009.

    GU S Y, HU Y, NIE Y L, CAI H S, SONG X L, HE X, CUI K G.Change of aggregate in Agricultural Reclamation Jiusan Administration Bureau black soil with different reclamation years.Journal of Northeast Agricultural University, 2015, 46(11): 56-62.doi:10.19720/j.cnki.issn.1005-9369.2015.11.009.(in Chinese)

    [19] 遲美靜, 侯瑋, 孫瑩, 范慶鋒, 虞娜, 鄒洪濤, 安晶, 張玉玲.東北黑土區(qū)荒地開墾種稻后土壤養(yǎng)分及pH值的變化特征.土壤通報(bào), 2018, 49(3): 546-551.doi:10.19336/j.cnki.trtb.2018.03.07.

    CHI M J, HOU W, SUN Y, FAN Q F, YU N, ZOU H T, AN J, ZHANG Y L.Characteristics of soil nutrients and pH value of paddy fields with different planted years in black soil region of northeast China.Chinese Journal of Soil Science, 2018, 49(3): 546-551.doi:10.19336/j.cnki.trtb.2018.03.07.(in Chinese)

    [20] 叢耀輝, 張玉玲, 張玉龍, 虞娜, 鄒洪濤, 范慶鋒, 王展.黑土區(qū)水稻土有機(jī)氮組分及其對(duì)可礦化氮的貢獻(xiàn).土壤學(xué)報(bào), 2016, 53(2): 457-467.doi:10.11766/trxb201508220362.

    CONG Y H, ZHANG Y L, ZHANG Y L, YU N, ZOU H T, FAN Q F, WANG Z.Soil organic nitrogen components and their contributions to mineralizable nitrogen in paddy soil of the black soil region.Acta Pedologica Sinica, 2016, 53(2): 457-467.doi:10.11766/trxb201508220362.(in Chinese)

    [21] ZHANG Y L, XU W J, DUAN P P, CONG Y H, AN T T, YU N, ZOU H T, DANG X L, AN J, FAN Q F, ZHANG Y L.Evaluation and simulation of nitrogen mineralization of paddy soils in Mollisols area of Northeast China under waterlogged incubation.PLoS One, 2017, 12(2): e0171022.doi:10.1371/journal.pone.0171022.

    [22] 張玉玲, 陳溫福, 虞娜, 付時(shí)豐, 張玉龍, 鄒洪濤.不同利用方式下土壤有機(jī)氮素礦化特征的研究.土壤通報(bào), 2013, 44(1): 52-56.doi:10.19336/j.cnki.trtb.2013.01.008.

    ZHANG Y L, CHEN W F, YU N, FU S F, ZHANG Y L, ZOU H T.Long-term effects of different land use patterns on mineralizing characteristic of soil organic nitrogen.Chinese Journal of Soil Science, 2013, 44(1): 52-56.doi:10.19336/j.cnki.trtb.2013.01.008.(in Chinese)

    [23] DAEBELER A, BODELIER P L E, HEFTING M M, RüTTING T, JIA Z J, LAANBROEK H J.Soil warming and fertilization altered rates of nitrogen transformation processes and selected for adapted ammonia-oxidizing Archaea in sub-Arctic grassland soil.Soil Biology and Biochemistry, 2017, 107: 114-124.doi:10.1016/j.soilbio.2016.12.013.

    [24] LANG M, LI P, TI C P, ZHU S X, YAN X Y, CHANG S X.Soil gross nitrogen transformations are related to land-uses in two agroforestry systems.Ecological Engineering, 2019, 127: 431-439.doi:10.1016/ j.ecoleng.2018.12.022.

    [25] 李平, 郎漫.開墾年限對(duì)黑土氮初級(jí)轉(zhuǎn)化速率和凈轉(zhuǎn)化速率的影響.土壤學(xué)報(bào), 2020, 57(1): 165-173.doi:10.11766/trxb201902180022.

    LI P, LANG M.Effect of cultivation on gross and net N transformation rates in black soil relative to duration.Acta Pedologica Sinica, 2020, 57(1): 165-173.doi:10.11766/trxb201902180022.(in Chinese)

    [26] 顧春朝.施肥類型和種植年限對(duì)稻田土壤氮礦化及可溶性有機(jī)碳、氮的影響[D].延吉: 延邊大學(xué), 2016.

    GU C C.Effects of different fertilization types and cultivation years on nitrogen mineralization and DOC, SON in paddy soil[D].Yanji: Yanbian University, 2016.(in Chinese)

    [27] KEENEY D R, BREMNER J M.Comparison and evaluation of laboratory methods of obtaining an index of soil nitrogen availability 1.Agronomy Journal, 1966, 58(5): 498-503.doi:10.2134/agronj1966.00021962005800050013x.

    [28] RICHTER J, NUSKE A, BOEHMER M, WEHRMANN J.Simulation of nitrogen mineralization and transport in loess-parabrownearthes: plot experiments.Plant and Soil, 1980, 54(3): 329-337.doi:10.1007/ bf02181829.

    [29] CABRERA M L.Modeling the flush of nitrogen mineralization caused by drying and rewetting soils.Soil Science Society of America Journal, 1993, 57(1): 63-66.doi:10.2136/sssaj1993.03615995005700010012x.

    [30] 姚槐應(yīng), 黃昌勇.土壤微生物生態(tài)學(xué)及其實(shí)驗(yàn)技術(shù).北京: 科學(xué)出版社, 2006.

    YAO H Y, HUANG C Y.Soil Microbial Ecology and Its Experimental Techniques.Beijing: Science Press, 2006.(in Chinese)

    [31] BOOTH M S, STARK J M, RASTETTER E.Controls on nitrogen cycling in terrestrial ecosystems: a synthetic analysis of literature data.Ecological Monographs, 2005, 75(2): 139-157.doi:10.1890/04-0988.

    [32] UEDA M U, KACHINA P, MAROD D, NAKASHIZUKA T, KUROKAWA H.Soil properties and gross nitrogen dynamics in old growth and secondary forest in four types of tropical forest in Thailand.Forest Ecology and Management, 2017, 398: 130-139.doi:10.1016/j.foreco.2017.05.010.

    [33] BARRETT J E, BURKE I C.Potential nitrogen immobilization in grassland soils across a soil organic matter gradient.Soil Biology and Biochemistry, 2000, 32(11/12): 1707-1716.doi:10.1016/s0038-0717 (00)00089-4.

    [34] ACCOE F, BOECKX P, BUSSCHAERT J, HOFMAN G, VAN CLEEMPUT O.Gross N transformation rates and net N mineralisation rates related to the C and N contents of soil organic matter fractions in grassland soils of different age.Soil Biology and Biochemistry, 2004, 36(12): 2075-2087.doi:10.1016/j.soilbio.2004.06.006.

    [35] CONRAD K A, DALAL R C, DALZELL S A, ALLEN D E, MENZIES N W.The sequestration and turnover of soil organic carbon in subtropical-grass pastures.Agriculture, Ecosystems & Environment, 2017, 248: 38-47.doi:10.1016/j.agee.2017.07.020.

    [36] 周碧青, 陳成榕, 楊文浩, 張黎明, 邢世和.茶樹對(duì)可溶性有機(jī)和無機(jī)態(tài)氮的吸收與運(yùn)轉(zhuǎn)特性.植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2017, 23(1): 189-195.doi:10.11674/zwyf.16067.

    ZHOU B Q, CHEN C R, YANG W H, ZHANG L M, XING S H.Uptake and transport characteristics of soluble organic and inorganic nitrogen by tea plant.Plant Nutrition and Fertilizer Science, 2017, 23(1): 189-195.doi:10.11674/zwyf.16067.(in Chinese)

    [37] MORRIS K A.Nitrogen cycling in the rhizosphere of cheatgrass and crested wheatgrass: contributions of root exudates and senescence.Longan: MS Thesis of Utah State University, 2014.

    [38] WANG F M, LI J, WANG X L, ZHANG W, ZOU B, NEHER D A, LI Z A.Nitrogen and phosphorus addition impact soil N?O emission in a secondary tropical forest of South China.Scientific Reports, 2014, 4: 5615.doi:10.1038/srep05615.

    [39] 張玉玲.東北地區(qū)水稻土供氮能力研究.北京: 中國(guó)農(nóng)業(yè)出版社, 2020.

    ZHANG Y L.Study of Nitrogen Supplying Capacity of Paddy Soil in Northeast China.Beijing: Chinese Agriculture Press, 2020.(in Chinese)

    Characteristics of Organic Nitrogen Mineralization in Paddy Soil with Different Reclamation Years in Black Soil of Northeast China

    GAO JiaRui, FANG ShengZhi, ZHANG YuLing*, AN Jing, YU Na, ZOU HongTao

    College of Land and Environment, Shenyang Agricultural University/Key Laboratory of Northeast Arable Land Conservation, Ministry of Agriculture and Rural Affairs, Shenyang 110866

    【Objective】 The aim of this study was to analyze the mineralized nitrogen (N) content, net N mineralized rate and net N mineralized ratio (the ratio of mineralized N to total N) in paddy soil with different reclamation years in black soil, and to explore the soil N supply capacity and its characteristics, and to reveal the soil N evolution law, so as to provide the theoretical basis for rational utilization and fertilization of black soil in Northeast China.【Method】 The natural wasteland (0 years, as the control soil, original natural meadow vegetation) and paddy soils with different reclamation years (12, 35, 62 and 85 a) (topography and cropping system, fertilization, and water management, roughly the same) in black soil region were selected as the research object, and the characteristics of soil organic N mineralization after cultivated rice form natural wasteland in black soil were studied by the water-logged incubation method.【Result】 During the early stages of incubation (about 1 month) , the cumulative mineralized N increased rapidly in each year, then showed a slow increase trend.At the end of incubation (297 d), the cumulative mineralized N ranged from 212.43 to 388.11 mg·kg-1, and the order of cumulative mineralized N was 0, 12, 35, 85 and 62 a.The mineralization process of soil organic N could be well described by a hybrid model (Special model), and the soil organic N pools could be divided into the increment N pool (the N pool made available after a drying and rewetting event) and the resistant N pool.Compared with the control soil (0 a), the potentially mineralisable N (N) of the increment N pool in all paddy soils showed a decreasing trend in each year.There was no significant difference between paddy soils of 12 and 35 years, as well as 62 and 85 years, but theNin paddy soils of 12 and 35 years were significantly higher than that of 62 and 85 years (<0.05).The rate constant (k) of mineralization of the increment N pool in paddy soils all showed an upward trend, but there was no significant difference betweenkof all paddy soils in each year (>0.05).Compared with the control soil (0 a), the rate constant (0) of mineralization of the resistant N pool in paddy soils of 12 and 35 years did not change significantly (>0.05), but0in 62 and 85 years decreased significantly (<0.05).The net N mineralization rate of the soils in each year were the largest at 4 days of incubation, and then decreased gradually.At the end of water-logged incubation (297 d), the order of the soil net N mineralization rate was consistent with that of the cumulative mineralization N.The net N mineralized ratio was relatively high at the beginning of incubation, and then increased slowly.At the end of incubation (297 d), the net N mineralized ratio ranged from 78.60 to 101.82 mg·g-1, and the order was 0, 35, 12, 85 and 62 a.Soil total N and C/N were important factors affecting the amount of mineralization N and the net N mineralization rate in paddy soils with different reclamation years (<0.05).The sum of initial mineral N andNcould be used to characterize the N supply capacity of paddy soil in rice growing season;compared with the control soil (0 a), the N supply capacity of paddy soil in each year decreased significantly (<0.05), and the soils of 12 and 35 years were significantly higher than that of 62 and 85 years (<0.05).【Conclusion】 During 85 years of rice cultivation from natural wasteland in black soil, the N supply capacity in paddy soil have declined, and the decline was significant after 35 years rice cultivation.Therefore, the improvement of soil organic matter content should be paid attention in the soil fertility cultivation in paddy fields.

    black soil region; paddy soil; rice cultivation year; mineralizated nitrogen; soil nitrogen suppying capacity

    2021-03-12;

    2021-11-11

    國(guó)家自然科學(xué)基金面上項(xiàng)目(41571280)

    高佳蕊,E-mail:1203997865@qq.com。通信作者張玉玲,E-mail:zhangyuling@syau.edu.cn

    (責(zé)任編輯 李云霞)

    猜你喜歡
    氮礦化黑土土壤有機(jī)
    自然資源部:加強(qiáng)黑土耕地保護(hù)
    氮磷添加對(duì)溫帶和亞熱帶森林土壤碳氮礦化的影響
    基于Sentinel-2遙感影像的黑土區(qū)土壤有效磷反演
    溫度變化對(duì)尕海濕地不同退化梯度土壤氮礦化的影響
    生物炭對(duì)杉木人工林土壤碳氮礦化的影響
    西雙版納橡膠林土壤有機(jī)碳分布特征研究
    秸稈還田的土壤有機(jī)碳周轉(zhuǎn)特征
    土壤有機(jī)碳轉(zhuǎn)化研究及其進(jìn)展
    AMDIS在土壤有機(jī)污染物鑒別中表征性統(tǒng)計(jì)量的探究
    典型黑土區(qū)不同尺度觀測(cè)場(chǎng)地融雪徑流
    国产成人免费观看mmmm| 91国产中文字幕| 80岁老熟妇乱子伦牲交| 一区二区三区激情视频| 在线av久久热| 黑丝袜美女国产一区| 在线观看一区二区三区激情| 久久精品久久久久久久性| av电影中文网址| 国产成人av教育| 亚洲九九香蕉| 午夜福利一区二区在线看| 亚洲成人免费av在线播放| 免费看av在线观看网站| 日韩免费高清中文字幕av| 亚洲欧美日韩高清在线视频 | 亚洲精品国产av蜜桃| 午夜激情av网站| e午夜精品久久久久久久| 人成视频在线观看免费观看| 国产成人精品在线电影| 精品少妇黑人巨大在线播放| 亚洲,欧美精品.| 久久精品久久久久久久性| 丰满少妇做爰视频| 国产男女内射视频| 久久精品亚洲av国产电影网| 老司机影院毛片| 国产亚洲欧美在线一区二区| 亚洲成av片中文字幕在线观看| 美女中出高潮动态图| 国产在线一区二区三区精| 国产精品一二三区在线看| 亚洲国产最新在线播放| 亚洲中文日韩欧美视频| 美女福利国产在线| 看免费av毛片| 成年人黄色毛片网站| 99re6热这里在线精品视频| 久久国产精品影院| 咕卡用的链子| 这个男人来自地球电影免费观看| 老汉色av国产亚洲站长工具| 波多野结衣av一区二区av| 晚上一个人看的免费电影| 亚洲欧美一区二区三区国产| 丝袜美足系列| 黄色a级毛片大全视频| 最新在线观看一区二区三区 | 国产欧美日韩综合在线一区二区| 中文字幕另类日韩欧美亚洲嫩草| 青春草视频在线免费观看| 国产精品久久久久久精品古装| 日韩免费高清中文字幕av| 欧美久久黑人一区二区| 少妇人妻久久综合中文| 汤姆久久久久久久影院中文字幕| 精品久久久久久久毛片微露脸 | 天天躁日日躁夜夜躁夜夜| 久久ye,这里只有精品| 国产一区二区在线观看av| 永久免费av网站大全| 国产成人精品无人区| 欧美在线一区亚洲| 波野结衣二区三区在线| 免费少妇av软件| 亚洲精品中文字幕在线视频| 妹子高潮喷水视频| 精品少妇一区二区三区视频日本电影| 丝袜美足系列| 婷婷色综合www| 男人爽女人下面视频在线观看| 啦啦啦在线免费观看视频4| 黄色a级毛片大全视频| 国产成人a∨麻豆精品| 人妻人人澡人人爽人人| 一区二区三区乱码不卡18| 午夜免费鲁丝| 色婷婷av一区二区三区视频| 久9热在线精品视频| 汤姆久久久久久久影院中文字幕| 久久精品成人免费网站| 一级片免费观看大全| 日本黄色日本黄色录像| 久久精品aⅴ一区二区三区四区| 欧美av亚洲av综合av国产av| 久久人人97超碰香蕉20202| 丝袜在线中文字幕| 在线观看人妻少妇| 日韩精品免费视频一区二区三区| av在线老鸭窝| 国产精品国产三级国产专区5o| 婷婷色综合www| 午夜福利影视在线免费观看| 一本综合久久免费| 午夜激情av网站| 国产精品久久久久久精品古装| 亚洲av综合色区一区| 国产精品熟女久久久久浪| 一区二区三区乱码不卡18| 又黄又粗又硬又大视频| 久久久久网色| 亚洲精品在线美女| 超碰成人久久| 91精品国产国语对白视频| 两个人看的免费小视频| av天堂久久9| 伦理电影免费视频| www.熟女人妻精品国产| 国产片内射在线| 黑人猛操日本美女一级片| 大码成人一级视频| 国产男人的电影天堂91| 中文字幕最新亚洲高清| 久久久久精品国产欧美久久久 | 亚洲国产精品国产精品| 欧美日韩亚洲高清精品| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一出视频| 久久精品国产综合久久久| 国产免费一区二区三区四区乱码| 人人澡人人妻人| 午夜免费鲁丝| 丰满少妇做爰视频| 水蜜桃什么品种好| netflix在线观看网站| 看免费av毛片| 亚洲欧洲精品一区二区精品久久久| 精品一区二区三区av网在线观看 | 欧美精品啪啪一区二区三区 | 国产成人精品在线电影| 久久人妻福利社区极品人妻图片 | av一本久久久久| 国产男人的电影天堂91| 国产激情久久老熟女| 丰满人妻熟妇乱又伦精品不卡| 校园人妻丝袜中文字幕| 一级毛片我不卡| 曰老女人黄片| 黄色毛片三级朝国网站| 真人做人爱边吃奶动态| 黄色a级毛片大全视频| 看免费av毛片| 午夜日韩欧美国产| 97人妻天天添夜夜摸| 午夜精品国产一区二区电影| 女人高潮潮喷娇喘18禁视频| 日本91视频免费播放| 考比视频在线观看| 欧美黄色片欧美黄色片| 久久av网站| www.999成人在线观看| 人人妻人人澡人人爽人人夜夜| 成年美女黄网站色视频大全免费| 美女高潮到喷水免费观看| 久久亚洲精品不卡| 另类亚洲欧美激情| 大香蕉久久网| 肉色欧美久久久久久久蜜桃| 久久鲁丝午夜福利片| 亚洲精品美女久久av网站| 亚洲精品美女久久久久99蜜臀 | 亚洲精品国产一区二区精华液| 欧美精品av麻豆av| xxxhd国产人妻xxx| 国产免费一区二区三区四区乱码| 成年人午夜在线观看视频| 亚洲欧美一区二区三区国产| 亚洲欧美日韩另类电影网站| 国产精品免费视频内射| 亚洲三区欧美一区| 超碰97精品在线观看| 亚洲国产成人一精品久久久| 91字幕亚洲| 精品人妻一区二区三区麻豆| 无遮挡黄片免费观看| 美女午夜性视频免费| 国产av一区二区精品久久| 18禁国产床啪视频网站| 精品人妻熟女毛片av久久网站| 亚洲精品美女久久av网站| 色94色欧美一区二区| 色网站视频免费| 男女免费视频国产| 日韩av在线免费看完整版不卡| 国产精品偷伦视频观看了| 精品国产乱码久久久久久小说| 成年女人毛片免费观看观看9 | 免费女性裸体啪啪无遮挡网站| av网站免费在线观看视频| 日韩人妻精品一区2区三区| 18禁黄网站禁片午夜丰满| 亚洲欧洲精品一区二区精品久久久| 咕卡用的链子| 一级黄片播放器| 亚洲欧美日韩另类电影网站| 欧美日韩亚洲高清精品| 波多野结衣av一区二区av| 欧美黄色片欧美黄色片| 满18在线观看网站| 国产成人一区二区三区免费视频网站 | 亚洲欧美一区二区三区黑人| av福利片在线| 欧美少妇被猛烈插入视频| 2021少妇久久久久久久久久久| 亚洲五月色婷婷综合| 亚洲精品一卡2卡三卡4卡5卡 | 看免费av毛片| 亚洲一区中文字幕在线| 啦啦啦啦在线视频资源| 性色av一级| 一个人免费看片子| 80岁老熟妇乱子伦牲交| 中文字幕另类日韩欧美亚洲嫩草| 汤姆久久久久久久影院中文字幕| 日本wwww免费看| 色视频在线一区二区三区| 欧美人与性动交α欧美软件| 一级a爱视频在线免费观看| 999久久久国产精品视频| 欧美变态另类bdsm刘玥| 久久久国产精品麻豆| 免费看十八禁软件| 两人在一起打扑克的视频| 啦啦啦在线免费观看视频4| 天堂俺去俺来也www色官网| 一个人免费看片子| 亚洲精品美女久久av网站| 欧美人与善性xxx| 你懂的网址亚洲精品在线观看| 最黄视频免费看| 欧美另类一区| 99国产精品一区二区三区| 久久国产精品影院| 日本91视频免费播放| 高清视频免费观看一区二区| 99国产综合亚洲精品| 国产精品成人在线| 日韩大片免费观看网站| 我的亚洲天堂| 亚洲午夜精品一区,二区,三区| 免费在线观看影片大全网站 | 免费在线观看完整版高清| 亚洲一码二码三码区别大吗| 亚洲国产av新网站| 搡老乐熟女国产| 99九九在线精品视频| 久久人人爽人人片av| 国产深夜福利视频在线观看| 尾随美女入室| 国产成人精品在线电影| 一区在线观看完整版| 成人18禁高潮啪啪吃奶动态图| 老司机在亚洲福利影院| 天堂俺去俺来也www色官网| 亚洲国产成人一精品久久久| 在线观看免费视频网站a站| 夫妻性生交免费视频一级片| 国产一区二区三区av在线| 亚洲精品国产色婷婷电影| 丝袜脚勾引网站| 国产深夜福利视频在线观看| 久久av网站| 亚洲激情五月婷婷啪啪| 亚洲精品国产一区二区精华液| 婷婷色综合www| 精品人妻1区二区| 在线天堂中文资源库| 久久天躁狠狠躁夜夜2o2o | 男人爽女人下面视频在线观看| 国产精品一二三区在线看| 国产精品久久久av美女十八| 岛国毛片在线播放| 波多野结衣一区麻豆| 人人澡人人妻人| 一区二区三区激情视频| 午夜免费男女啪啪视频观看| 男女之事视频高清在线观看 | xxx大片免费视频| 丁香六月欧美| 亚洲欧美精品综合一区二区三区| 国产亚洲av片在线观看秒播厂| 90打野战视频偷拍视频| 午夜老司机福利片| 午夜福利视频精品| 啦啦啦啦在线视频资源| 熟女av电影| 成年人黄色毛片网站| xxxhd国产人妻xxx| 一边摸一边做爽爽视频免费| 日韩视频在线欧美| 国产成人精品久久二区二区91| av又黄又爽大尺度在线免费看| 亚洲精品国产av成人精品| 亚洲第一av免费看| 中文字幕人妻丝袜制服| 黄色毛片三级朝国网站| a级毛片黄视频| 午夜影院在线不卡| 久久精品国产亚洲av高清一级| 女人精品久久久久毛片| 精品久久久精品久久久| 久久久精品区二区三区| kizo精华| 男的添女的下面高潮视频| 一级毛片 在线播放| 99国产精品99久久久久| 少妇粗大呻吟视频| 中文字幕人妻丝袜一区二区| 丝袜人妻中文字幕| 男女边吃奶边做爰视频| 色婷婷久久久亚洲欧美| videos熟女内射| 捣出白浆h1v1| 成人手机av| 久久99热这里只频精品6学生| 国产亚洲一区二区精品| 精品视频人人做人人爽| 精品久久久久久久毛片微露脸 | 日韩 欧美 亚洲 中文字幕| 天天添夜夜摸| 日本色播在线视频| 亚洲国产精品一区三区| 秋霞在线观看毛片| 亚洲精品久久成人aⅴ小说| cao死你这个sao货| 欧美精品av麻豆av| 亚洲五月色婷婷综合| 亚洲久久久国产精品| 1024视频免费在线观看| 日韩大片免费观看网站| 国产成人a∨麻豆精品| 亚洲伊人色综图| videosex国产| 日韩一卡2卡3卡4卡2021年| 亚洲欧美清纯卡通| 国产精品久久久av美女十八| 久久女婷五月综合色啪小说| 久久这里只有精品19| 国产在视频线精品| 男人添女人高潮全过程视频| 1024视频免费在线观看| 人人澡人人妻人| 亚洲精品国产区一区二| 久久久欧美国产精品| 天堂中文最新版在线下载| 18禁黄网站禁片午夜丰满| 午夜福利视频精品| av在线老鸭窝| 国产精品久久久久久人妻精品电影 | 高清av免费在线| 91精品伊人久久大香线蕉| 每晚都被弄得嗷嗷叫到高潮| tube8黄色片| 亚洲七黄色美女视频| 欧美中文综合在线视频| 黄网站色视频无遮挡免费观看| 成年人免费黄色播放视频| 性少妇av在线| 欧美日韩精品网址| 国产淫语在线视频| 伊人亚洲综合成人网| 黑人欧美特级aaaaaa片| 久久精品国产a三级三级三级| 亚洲av电影在线进入| 两人在一起打扑克的视频| 99国产精品一区二区三区| 国产精品久久久人人做人人爽| 在线观看免费日韩欧美大片| 黄色毛片三级朝国网站| 在线观看一区二区三区激情| 最近最新中文字幕大全免费视频 | 97精品久久久久久久久久精品| 亚洲欧美色中文字幕在线| 永久免费av网站大全| 在线观看国产h片| 999久久久国产精品视频| 午夜福利免费观看在线| 精品国产一区二区久久| 真人做人爱边吃奶动态| 乱人伦中国视频| 成在线人永久免费视频| 精品欧美一区二区三区在线| 精品福利观看| 欧美黄色片欧美黄色片| 99精品久久久久人妻精品| 9191精品国产免费久久| 男女午夜视频在线观看| 亚洲一码二码三码区别大吗| 夫妻午夜视频| 亚洲精品国产区一区二| 视频区欧美日本亚洲| 欧美另类一区| 蜜桃国产av成人99| 国产av精品麻豆| 免费av中文字幕在线| 亚洲熟女毛片儿| 久久久久精品人妻al黑| 精品国产乱码久久久久久小说| 丝袜美足系列| av在线老鸭窝| 男的添女的下面高潮视频| 一级毛片电影观看| 黄网站色视频无遮挡免费观看| 欧美人与善性xxx| 操出白浆在线播放| 亚洲欧美一区二区三区久久| 欧美精品高潮呻吟av久久| 国产成人精品久久二区二区91| 欧美另类一区| 热re99久久精品国产66热6| 五月开心婷婷网| 最新在线观看一区二区三区 | 国产成人欧美| 久久久国产精品麻豆| 成人国产av品久久久| 欧美黑人精品巨大| 欧美日韩成人在线一区二区| 脱女人内裤的视频| 精品人妻1区二区| 国产成人精品无人区| 搡老乐熟女国产| 在线观看免费午夜福利视频| 大片免费播放器 马上看| 国产又爽黄色视频| 又大又爽又粗| 精品国产国语对白av| 午夜福利乱码中文字幕| 欧美日韩福利视频一区二区| 国产老妇伦熟女老妇高清| 国产欧美日韩精品亚洲av| 一本—道久久a久久精品蜜桃钙片| 美国免费a级毛片| 亚洲精品久久成人aⅴ小说| 欧美成狂野欧美在线观看| 久久国产亚洲av麻豆专区| 国产又爽黄色视频| 亚洲人成电影免费在线| 丝袜脚勾引网站| 性高湖久久久久久久久免费观看| 国产精品人妻久久久影院| 亚洲午夜精品一区,二区,三区| 国产欧美日韩精品亚洲av| 大片电影免费在线观看免费| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 国精品久久久久久国模美| 国产免费又黄又爽又色| 麻豆国产av国片精品| 一区二区三区四区激情视频| 日本五十路高清| 日韩一本色道免费dvd| 精品高清国产在线一区| 国产不卡av网站在线观看| 日韩大片免费观看网站| 亚洲精品中文字幕在线视频| 母亲3免费完整高清在线观看| 黄色毛片三级朝国网站| 日韩电影二区| 国产色视频综合| 2021少妇久久久久久久久久久| 在线精品无人区一区二区三| 欧美黄色淫秽网站| 一边摸一边做爽爽视频免费| 欧美激情高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 精品卡一卡二卡四卡免费| 国产精品久久久av美女十八| 免费在线观看黄色视频的| 一级a爱视频在线免费观看| 99精国产麻豆久久婷婷| 久久久久久久久免费视频了| 丝袜喷水一区| 中文字幕人妻熟女乱码| 狂野欧美激情性bbbbbb| 国产三级黄色录像| 久久影院123| 久久久国产欧美日韩av| 99热国产这里只有精品6| 99热全是精品| 高清黄色对白视频在线免费看| 男女高潮啪啪啪动态图| 曰老女人黄片| 午夜福利免费观看在线| 精品久久久久久电影网| 国产精品久久久久成人av| 在线精品无人区一区二区三| av网站在线播放免费| 嫩草影视91久久| 久久精品亚洲av国产电影网| 女性被躁到高潮视频| 美女高潮到喷水免费观看| 丰满人妻熟妇乱又伦精品不卡| 国产野战对白在线观看| 一区在线观看完整版| 狂野欧美激情性bbbbbb| 叶爱在线成人免费视频播放| 久久精品久久精品一区二区三区| 久久精品aⅴ一区二区三区四区| 女人精品久久久久毛片| 中文字幕亚洲精品专区| 日韩av在线免费看完整版不卡| 国产在视频线精品| 久久久久精品人妻al黑| 亚洲国产中文字幕在线视频| av线在线观看网站| 大香蕉久久成人网| www日本在线高清视频| 捣出白浆h1v1| 精品一区在线观看国产| 丁香六月天网| 两性夫妻黄色片| 9热在线视频观看99| 夫妻性生交免费视频一级片| 久久综合国产亚洲精品| 美女中出高潮动态图| 亚洲国产av影院在线观看| 精品国产超薄肉色丝袜足j| 日本a在线网址| 日日摸夜夜添夜夜爱| 1024香蕉在线观看| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 欧美日韩成人在线一区二区| 男人爽女人下面视频在线观看| 女性生殖器流出的白浆| 欧美日韩视频高清一区二区三区二| 91麻豆精品激情在线观看国产 | 久久久国产欧美日韩av| 亚洲天堂av无毛| 麻豆乱淫一区二区| 纵有疾风起免费观看全集完整版| 99久久99久久久精品蜜桃| 老司机影院毛片| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 日日夜夜操网爽| 国产av精品麻豆| 飞空精品影院首页| 手机成人av网站| a级毛片黄视频| 久久久久精品人妻al黑| 欧美黑人欧美精品刺激| 人人妻人人澡人人看| 亚洲专区中文字幕在线| 亚洲欧洲国产日韩| 国产成人精品久久二区二区免费| 久久毛片免费看一区二区三区| 波多野结衣一区麻豆| 午夜日韩欧美国产| 亚洲av成人不卡在线观看播放网 | 狠狠婷婷综合久久久久久88av| 夫妻午夜视频| 国产淫语在线视频| 日韩av免费高清视频| 成人午夜精彩视频在线观看| 91麻豆av在线| av天堂在线播放| 日韩av不卡免费在线播放| 啦啦啦中文免费视频观看日本| 成年动漫av网址| 人人妻人人澡人人看| 啦啦啦啦在线视频资源| 国产有黄有色有爽视频| 91精品伊人久久大香线蕉| 国产三级黄色录像| 少妇猛男粗大的猛烈进出视频| 91麻豆精品激情在线观看国产 | 天天操日日干夜夜撸| 中文字幕精品免费在线观看视频| 欧美日韩综合久久久久久| netflix在线观看网站| www.av在线官网国产| 久久久久久久精品精品| 欧美黑人精品巨大| 性少妇av在线| 午夜激情av网站| 国产精品av久久久久免费| 精品一区二区三区av网在线观看 | 日韩熟女老妇一区二区性免费视频| 手机成人av网站| 一区二区av电影网| 国产欧美日韩一区二区三 | 手机成人av网站| 亚洲国产av影院在线观看| 国产野战对白在线观看| 99久久精品国产亚洲精品| 成人国产一区最新在线观看 | 热99国产精品久久久久久7| av有码第一页| 两个人看的免费小视频| 久久精品熟女亚洲av麻豆精品| 麻豆国产av国片精品| 黄网站色视频无遮挡免费观看| 亚洲一码二码三码区别大吗| 久久精品亚洲熟妇少妇任你| 中国美女看黄片| av在线app专区| 精品一区二区三区av网在线观看 | 久久人人97超碰香蕉20202| 日韩一卡2卡3卡4卡2021年| 晚上一个人看的免费电影| 午夜免费鲁丝| 蜜桃在线观看..| 我的亚洲天堂| 黄色视频不卡| 久久热在线av| 欧美变态另类bdsm刘玥| 夫妻午夜视频| 国产欧美日韩精品亚洲av| 日日夜夜操网爽| 大香蕉久久网| 日日摸夜夜添夜夜爱| 国产一区二区激情短视频 | 啦啦啦在线观看免费高清www|