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

    施用生物炭對菜地土壤氮素?fù)p失影響的模擬

    2022-04-15 08:46:28李雙建李桂花薛冬梅張貴龍
    關(guān)鍵詞:施肥量作物生物

    李雙建,李 虎,李 潔,李桂花,薛冬梅,張貴龍

    (1.天津師范大學(xué) 天津市水資源與水環(huán)境重點(diǎn)實(shí)驗(yàn)室,天津 300387;2.天津師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院,天津 300387;3.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081;4.農(nóng)業(yè)農(nóng)村部環(huán)境保護(hù)科研監(jiān)測所,天津 300191)

    蔬菜種植過程中普遍存在過量施肥現(xiàn)象.多次過量施肥會導(dǎo)致肥料中的氮、磷養(yǎng)分不能完全被蔬菜吸收,殘留在土壤中,并隨降水、灌溉等過程產(chǎn)生徑流進(jìn)入地表水或以淋溶的方式進(jìn)入深層土壤,進(jìn)而增加了氮、磷養(yǎng)分進(jìn)入地下水的風(fēng)險[1].長期施肥也會造成土壤板結(jié)、通氣性變差等問題[2].為減緩菜地土壤的劣化進(jìn)程,在減少施肥量的同時需對已被破壞的土壤進(jìn)行修復(fù),此外還要兼顧作物產(chǎn)量及經(jīng)濟(jì)效益.生物炭是近幾年來較為常用的土壤改良劑,較高的含碳率、較大的比表面積、多微孔結(jié)構(gòu)及特殊的分子組成使其具有較強(qiáng)的吸附能力和抗氧化能力[3],這些特點(diǎn)是生物炭改善土壤理化性質(zhì)、固持土壤養(yǎng)分、減少氮損失的重要基礎(chǔ).

    目前針對農(nóng)業(yè)面源污染區(qū)氮損失的研究較多,為提高研究精度并降低工作量,不同用途的氮損失模擬模型作為科學(xué)工具應(yīng)運(yùn)而生,常用的包括Root Zone Water Quality Model(RZWQM)、水和氮質(zhì)量平衡模型、DeNitrification-DeComposition(DNDC)模型等.DNDC模型是運(yùn)用農(nóng)業(yè)系統(tǒng)中一系列控制碳和氮遷移轉(zhuǎn)化的生物化學(xué)及地球化學(xué)反應(yīng)機(jī)制的計(jì)算機(jī)模擬模型[4],在農(nóng)田氮損失研究中得到了廣泛應(yīng)用[5-7].本研究將施用生物炭后土壤的理化性質(zhì)變化輸入模型進(jìn)行菜地氮損失情況模擬,驗(yàn)證施用生物炭對菜地氮損失的影響.

    本研究借助DNDC模型:①確定影響作物產(chǎn)量及氮損失的重要因素;②定量生物炭與不同施肥量處理的綜合氮損失;③驗(yàn)證施用生物炭降低氮肥施用量的效率以及對氮損失的削減效率.

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)為位于黃淮海平原地區(qū)的山東省濟(jì)寧市金鄉(xiāng)縣雞黍鎮(zhèn)大蒜/辣椒套作菜地(34.96°N,116.23°E).該地區(qū)屬暖溫帶季風(fēng)氣候,具有明顯的冬夏季風(fēng)特點(diǎn),年平均氣溫為13.3~14.1℃,平均無霜期為199 d,年平均降水量在700 mm左右,其中夏季降水量占全年降水量的65%左右[8].研究區(qū)土壤類型為潮土,土壤黏粒質(zhì)量分?jǐn)?shù)約為37.8%~40.5%,砂礫質(zhì)量分?jǐn)?shù)約為7.3%~8.8%.供試土壤理化性質(zhì)如表1所示.

    表1 供試土壤理化性質(zhì)Tab.1 Physical and chemical properties of tested soil

    1.2 實(shí)驗(yàn)處理

    研究區(qū)內(nèi)大蒜種植周期一般為當(dāng)年10月至次年5月;辣椒種植周期為當(dāng)年2月至9月,其中2月至4月為育苗期,4月移栽套種.實(shí)驗(yàn)共設(shè)置常規(guī)種植(CK)、生物炭+推薦施肥(B-N)2種處理,每種處理3次重復(fù),共6個實(shí)驗(yàn)小區(qū),每個小區(qū)面積為44.65 m2.CK處理的種植管理措施與當(dāng)?shù)赝瓿R?guī)種植平均水平保持一致;B-N處理的種植管理措施為:播種前按100 g/m2的量添加生物炭作為土壤改良劑,并按推薦施肥量施肥,其他管理措施與CK處理相同.不同處理種植過程的施肥情況如表2所示.

    表2 不同處理施肥措施Tab.2 Fertilization measures of different treatments

    CK處理的施肥方案:大蒜種植季一次施入復(fù)合肥(尿素、P2O5、K2O),大蒜收獲后在辣椒生長期內(nèi)再適量追肥;B-N處理的施肥方案:在大蒜種植季將CK處理的單次施肥量按原施肥量的2/3和1/5分2次施用,并在大蒜種植前按100 g/m2的用量施用生物炭.

    1.3 研究方法

    DNDC模型在初次運(yùn)行前需要進(jìn)行參數(shù)本地化設(shè)置,將樣地的實(shí)際監(jiān)測參數(shù)輸入模型,并以此驅(qū)動模型運(yùn)行,根據(jù)輸出結(jié)果與實(shí)際監(jiān)測參數(shù)的誤差進(jìn)一步確定模擬參數(shù).利用模型中的敏感性分析模塊,借助獨(dú)立參數(shù)擾動法[9],以生物量、作物產(chǎn)量以及氮損失等為目標(biāo),對輸入模型的參數(shù)進(jìn)行敏感性分析,進(jìn)而確定影響生物量、作物產(chǎn)量以及氮損失等的關(guān)鍵敏感性參數(shù).根據(jù)敏感性分析結(jié)果,重點(diǎn)對影響目標(biāo)結(jié)果的關(guān)鍵參數(shù)進(jìn)行本地初始化設(shè)置并進(jìn)行模型模擬,根據(jù)模擬結(jié)果確定部分參數(shù),進(jìn)而提高模擬精度.本研究分別采用常規(guī)種植和生物炭+推薦施肥兩種處理方式,通過DNDC模型模擬作物產(chǎn)量以及土壤表層無機(jī)氮含量變化等來驗(yàn)證模型模擬精度,并推算此精度下不同種植方式過程中的氮損失情況.

    將陰干土按水土質(zhì)量比1∶2.5振蕩后靜置過濾,用便攜式水質(zhì)參數(shù)分析儀(美國Thermo Scientific Orion)測定土壤pH值和鹽度.用2 mol/L的KCl溶液從鮮土中浸提NH4+-N和NO3--N,土壤與浸提液的質(zhì)量比為1∶5,混合振蕩后靜置過濾,之后利用連續(xù)流動分析儀(Auto Seal 3,德國Seal公司)測定其含量.大蒜植株樣品的總碳(TC)和總氮(TN)含量使用元素分析儀(Vario EL III,德國Elementar公司)進(jìn)行測定.

    1.4 參數(shù)敏感性分析及模擬精度分析

    參數(shù)的敏感性表示模型輸入?yún)?shù)的變化對于輸出結(jié)果影響的強(qiáng)弱.依據(jù)Walker等[10]提出的公式計(jì)算敏感性指數(shù):

    式中:S為相對敏感性指數(shù);I1、I2分別為輸入?yún)?shù)的最小值和最大值;Iavg是I1和I2的平均值;O1、O2分別為對應(yīng)于I1和I2的輸出值;Oavg為O1和O2的平均值.S的絕對值大小表明影響程度強(qiáng)弱,絕對值越大,輸入值對輸出值的影響就越大;S值為負(fù),說明輸入值與輸出值呈負(fù)相關(guān).

    把菜地所在地的氣溫、降水量、土壤本底養(yǎng)分含量、播種與收獲日期、施肥量、灌溉量、耕地深度等田間監(jiān)測數(shù)據(jù)輸入到DNDC模型中,并進(jìn)行模型初始運(yùn)行,將模擬結(jié)果與田間實(shí)測結(jié)果進(jìn)行比較,進(jìn)而確定模型模擬的準(zhǔn)確性.根據(jù)Gjettermanna等[11]提出的歸一化均方根誤差(NRMSE)和Willmott等[12]提出的一致性指標(biāo)d驗(yàn)證模型模擬精度.這兩個評價指標(biāo)的具體公式為:

    式中:Pi、Oi、分別為模擬值、實(shí)測值、實(shí)測值平均值;n為實(shí)測次數(shù).NRMSE<10%時,模擬效果為優(yōu);10%≤NRMSE<20%時,模擬效果為良;20%≤NEMSE<30%時,模擬效果為中等;NRMSE≥30%時,模擬效果為差[11].d∈(-∞,1),越接近于1,表明模擬效果越好[12].

    本研究中所需的土壤養(yǎng)分?jǐn)?shù)據(jù)、施肥灌溉等種植管理數(shù)據(jù)以及作物植株C/N等數(shù)據(jù)均來自于研究區(qū)實(shí)地取樣監(jiān)測與記錄;土壤的容重、作物生長積溫、最適生長溫度和作物需水量等數(shù)據(jù)均通過查閱相關(guān)文獻(xiàn)獲得.模擬所使用的模型為DNDC95,采用SPSS 25對數(shù)據(jù)進(jìn)行處理與分析,采用Origin 2018繪制圖表.

    2 結(jié)果與分析

    2.1 模型參數(shù)確定

    利用2018—2019年常規(guī)種植處理下大蒜/辣椒套作模式的實(shí)際監(jiān)測數(shù)據(jù)進(jìn)行模型的參數(shù)校正,用作物產(chǎn)量、土壤無機(jī)氮含量模擬值與實(shí)測值的擬合度來校正作物參數(shù).由于模型作物數(shù)據(jù)庫中沒有大蒜數(shù)據(jù),因此需在大蒜及辣椒收獲后取樣測定相關(guān)參數(shù),并根據(jù)結(jié)果對模型內(nèi)置辣椒相關(guān)參數(shù)進(jìn)行調(diào)整.最終確定的作物相關(guān)參數(shù)如表3所示.其中,作物最大生物量根據(jù)當(dāng)?shù)貧v史記錄數(shù)據(jù)確定,通過測定植株鮮重和干重計(jì)算植株含水率,將作物最大產(chǎn)量轉(zhuǎn)換為干重,使用模型推薦生物量轉(zhuǎn)換參數(shù)(0.6)將其轉(zhuǎn)換為最大生物量.

    表3 作物相關(guān)參數(shù)Tab.3 Crop related parameters

    2.2 DNDC模型模擬效果驗(yàn)證

    2.2.1 作物產(chǎn)量驗(yàn)證

    借助DNDC模型的作物生長子模型模擬生長過程,將產(chǎn)量的模擬結(jié)果與田間實(shí)測結(jié)果進(jìn)行對比,用以驗(yàn)證模型對作物生長過程的模擬精度.將不同處理的大蒜和辣椒的模型模擬產(chǎn)量分別與實(shí)測產(chǎn)量進(jìn)行比較,檢驗(yàn)?zāi)P蛯ΨN植作物的產(chǎn)量模擬是否準(zhǔn)確,結(jié)果如表4所示.

    表4 各處理中作物產(chǎn)量的模擬結(jié)果Tab.4 Simulation results of crop yield of each treatment

    由表4模擬結(jié)果可以看出,CK處理的菜地中大蒜產(chǎn)量模擬值比實(shí)測值低1.84%,辣椒產(chǎn)量模擬值比實(shí)測值高16.85%.B-N處理的菜地中大蒜產(chǎn)量模擬值比實(shí)測值高1.55%,辣椒產(chǎn)量模擬值比實(shí)測值高13.91%.按作物種類來看,模型對大蒜產(chǎn)量模擬的誤差相對較小,在2%以內(nèi);對辣椒產(chǎn)量模擬的誤差相對較大,但總體誤差控制在17%以內(nèi).對模擬效果進(jìn)行評估,結(jié)果顯示,大蒜產(chǎn)量模擬結(jié)果的NRMSE為3.18%,d為0.97;辣椒產(chǎn)量模擬結(jié)果的NRMSE為14.32%,d為0.93,模擬效果均較好.綜合來看,模型對CK處理和B-N處理的作物產(chǎn)量模擬結(jié)果充分證明了模型對大蒜/辣椒套作模式的產(chǎn)量具有一定的模擬精度,為進(jìn)一步模擬2種處理的氮損失情況提供了可靠的支撐.

    2.2.2 土壤無機(jī)氮含量驗(yàn)證

    向DNDC模型輸入土壤理化性質(zhì)本底參數(shù),模擬作物生長的同時模擬土壤0~40 cm深度部分10 cm分層的無機(jī)氮含量(以N計(jì))隨施肥、灌溉及作物生長等變化的逐日變化.在作物生長周期內(nèi)采集不同小區(qū)內(nèi)土壤樣品并測定其無機(jī)氮含量,從而驗(yàn)證模型對土壤無機(jī)氮的模擬精度.由于所采樣品為20 cm分層,而模型模擬結(jié)果為10 cm分層,故使用0~10 cm及10~20 cm土壤無機(jī)氮模擬結(jié)果的平均值替代0~20 cm的數(shù)據(jù),20~30 cm和30~40 cm土壤無機(jī)氮模擬結(jié)果的平均值替代20~40 cm的數(shù)據(jù),結(jié)果如圖1和圖2所示.

    圖1 土壤NH4+-N含量模擬結(jié)果Fig.1 Simulation results of NH4+-N content in soil

    圖2 土壤NO3-N含量模擬結(jié)果Fig.2 Simulation results of NO3-N content in soil

    由圖1和圖2可以看出,模型對2個處理組無機(jī)氮含量模擬效果均是0~20 cm土層深度的較好,20~40 cm土層深度的較差.進(jìn)一步驗(yàn)證模擬精度發(fā)現(xiàn),模型對CK處理0~20 cm土壤中NH4+-N含量模擬結(jié)果的NRMSE為21.23%,一致性指標(biāo)d為0.67;對NO3--N含量模擬結(jié)果的NRMSE為10.86%,一致性指標(biāo)d為0.91.對B-N處理0~20 cm土壤中NH4+-N含量模擬結(jié)果的NRMSE為20.65%,一致性指標(biāo)d為0.63;對NO3--N含量模擬結(jié)果的NRMSE為25.53%,一致性指標(biāo)d為0.35.模型對CK處理20~40 cm土壤中NH4+-N含量模擬結(jié)果的NRMSE為27.40%,一致性指標(biāo)d為0.56;對NO3--N含量模擬結(jié)果的NRMSE為32.65%,一致性指標(biāo)d為0.17.對B-N處理20~40 cm土壤中NH4+-N含量模擬結(jié)果的NRMSE為24.81%,一致性指標(biāo)d為0.47;對NO3--N含量模擬結(jié)果的NRMSE為32.68%,一致性指標(biāo)d為0.37.綜合來看,模擬結(jié)果的誤差表現(xiàn)為中等,誤差總體在可允許范圍內(nèi).總的來看,模型對土壤0~20 cm深度無機(jī)氮含量的模擬效果較為理想,但對20~40 cm深度無機(jī)氮的模擬效果相對較差.由于0~20 cm深度為作物根系的主要活動區(qū)域,因此可以認(rèn)為模型對作物生長過程的氮收支過程具有較好的模擬效果.

    2.3 模型的參數(shù)敏感性分析

    本研究主要從氣象因素、土壤的理化性質(zhì)參數(shù)、作物參數(shù)以及種植管理措施4個方面對種植作物的產(chǎn)量及種植過程的NH3揮發(fā)損失、N2O排放、NO3--N淋溶影響的強(qiáng)弱進(jìn)行分析.氣象因素主要為溫度、降水量;土壤的理化性質(zhì)參數(shù)為土壤pH值及土壤總有機(jī)碳(SOC)、NO3--N和NH4+-N含量;種植管理措施參數(shù)主要為耕地深度、施肥量、施肥深度.除土壤NO3--N含量參數(shù)變異范圍為±50%外,其他參數(shù)變異范圍均為±20%,模型運(yùn)行次數(shù)默認(rèn)為4 000次.不同參數(shù)對作物產(chǎn)量、N2O釋放、NH3揮發(fā)、NO3--N淋溶具有不同的敏感性,分析結(jié)果如圖3所示.

    圖3 參數(shù)敏感性分析結(jié)果Fig.3 Results of parameter sensitivity analysis

    從圖3可以看出,對作物產(chǎn)量影響較大的參數(shù)為施肥量、土壤本底SOC、降水量、土壤pH值等,敏感性指數(shù)值分別為0.62、0.15、0.13、0.11、0.09.綜合來看,施肥、灌溉等種植管理措施對作物產(chǎn)量影響最大.對土壤N2O釋放影響較大的參數(shù)分別為施肥量、土壤pH值、溫度、土壤本底SOC、土壤本底NO3--N等,其敏感性指數(shù)值分別為2.3、2.27、0.53、0.52、0.41.施肥量、土壤本底SOC、土壤pH值、耕地深度等參數(shù)對土壤NH3揮發(fā)具有較強(qiáng)的影響,敏感性指數(shù)值分別為1.11、0.87、0.29、0.22.氣象因素、土壤理化性質(zhì)參數(shù)以及種植管理參數(shù)均對菜地NO3--N淋溶有較強(qiáng)的影響,其中影響較強(qiáng)的有施肥量、溫度、降水量,敏感性指數(shù)值分別為2.86、2.21、1.18.綜合來看,施肥量對作物產(chǎn)量及N2O釋放、NH3揮發(fā)、NO3--N淋溶等氮損失的影響均表現(xiàn)為最強(qiáng).除施肥量外,不同途徑的氮損失主要影響因素存在差異,土壤pH值對N2O釋放具有重要影響,而NH3揮發(fā)主要受土壤本底SOC的影響,NO3--N淋溶則主要受溫度及降水量的影響.

    2.4 施用生物炭效果的分析模擬

    2.4.1 施用生物炭對作物產(chǎn)量的影響

    作物收獲后的實(shí)際產(chǎn)量及氮、磷投入情況如表5所示.CK處理的大蒜產(chǎn)量比B-N處理的產(chǎn)量高2.2%,CK處理的辣椒產(chǎn)量比B-N處理的產(chǎn)量高2.6%,對不同處理間、不同作物間的產(chǎn)量進(jìn)行差異統(tǒng)計(jì)學(xué)意義分析,結(jié)果顯示,相同作物不同處理間產(chǎn)量的差異不具有統(tǒng)計(jì)學(xué)意義(大蒜產(chǎn)量P=0.914,辣椒產(chǎn)量P=0.888).從氮、磷投入量來看,B-N處理相較于CK處理,氮投入量減少了11.9%,磷投入量減少了12.7%,作物產(chǎn)量降低在3%以內(nèi),表明添加生物炭并改進(jìn)施肥模式能在有效減少氮、磷投入的同時保證作物的產(chǎn)量,即在保證經(jīng)濟(jì)效益的同時降低菜地土壤的氮、磷流失風(fēng)險.施用生物炭能降低土壤容重,提高土壤孔隙度,使作物根系獲得更大的生長空間,進(jìn)而增強(qiáng)生理功能,提高作物產(chǎn)量[14],這也是生物炭能提高作物產(chǎn)量的原因之一.

    表5 作物產(chǎn)量及氮磷投入量Tab.5 Crop yield and nitrogen and phosphorus input

    由于本研究中生物炭和施肥方式2個變量同時存在,不能明確B-N處理中的作物減產(chǎn)情況與施肥量及生物炭間的關(guān)系,故需要補(bǔ)充驗(yàn)證農(nóng)戶常規(guī)種植+推薦施肥(CK-N)處理及使用生物炭+常規(guī)種植施肥(B)處理下的作物產(chǎn)量情況.DNDC模型不僅可以用于模擬氮損失,也可用于田間實(shí)驗(yàn)的模擬研究[15-16].借助DNDC模型模擬CK-N及B兩種處理下的作物產(chǎn)量,模擬結(jié)果如圖4所示.

    圖4 菜地模擬實(shí)驗(yàn)作物產(chǎn)量Fig.4 Crop yield of vegetable field of simulation experiment

    從模擬實(shí)驗(yàn)作物產(chǎn)量來看,無論何種處理對大蒜產(chǎn)量的影響都較小,其中CK、CK-N、B處理的大蒜產(chǎn)量一致,而B-N處理的作物產(chǎn)量與其他處理相比提高了0.03%,幾乎沒有差距.不同處理間辣椒產(chǎn)量具有較大差異,其中CK與B處理的辣椒產(chǎn)量相同,且為最高水平,這表明生物炭的使用并不會導(dǎo)致作物的減產(chǎn);而CK-N處理及B-N處理相較于常規(guī)施肥模式下的CK處理與B處理均發(fā)生減產(chǎn),其中CK-N處理減產(chǎn)44.85%,B-N處理減產(chǎn)4.99%,這表明本研究中的作物減產(chǎn)主要受施肥量的影響.添加了生物炭處理的減產(chǎn)程度低于常規(guī)種植的減產(chǎn)程度,間接證明了生物炭的增產(chǎn)效果.推測導(dǎo)致減產(chǎn)的原因?yàn)椋菏褂猛扑]施肥方式后,投入的氮養(yǎng)分減少,導(dǎo)致土壤中累積的養(yǎng)分減少.5月份在大蒜收獲后,辣椒生長期長期高溫多雨的氣象條件導(dǎo)致土壤氮養(yǎng)分流失加劇,辣椒養(yǎng)分短缺;而在添加了生物炭的菜地,生物炭的吸附作用可將養(yǎng)分有效存儲在土壤中,并在后續(xù)的作物生長過程供作物吸收利用.

    2.4.2 施用生物炭對土壤表層無機(jī)氮的影響

    根據(jù)DNDC模型輸出的土壤無機(jī)氮含量變化(圖1、圖2)可以看出,施肥后CK與B-N處理的土壤表層(0~20 cm)中NH4+-N含量存在較大差異,而NO3--N的含量基本保持一致.受施肥量影響,土壤NH4+-N和NO3--N含量均表現(xiàn)出CK處理高于B-N處理,但在作物收獲后兩種處理方式的無機(jī)氮含量均降到較低且含量接近.CK處理的表層土壤NH4+-N含量在10月上旬因?yàn)楣喔鹊挠绊懓l(fā)生斷崖式下降,而同樣灌溉的B-N處理的表層土壤NH4+-N含量僅發(fā)生了較小的波動后便平緩下降,表明B-N處理中施用生物炭有效減緩了土壤表層NH4+-N向深層的運(yùn)移.

    2.4.3 施用生物炭對氮損失的影響

    農(nóng)田的氮損失主要以徑流、淋溶、氨揮發(fā)和硝化-反硝化4種方式進(jìn)行.其中:徑流損失主要損失可溶性速效氮及吸附于泥沙等顆粒的固態(tài)無機(jī)態(tài)氮和有機(jī)態(tài)氮;淋溶損失主要以NO3--N淋失為主;氨揮發(fā)和硝化-反硝化的損失主要以NH3、NO、N2O、N2等氣態(tài)損失為主[17].DNDC模型內(nèi)置的地表徑流模擬模塊為大尺度流域模擬,對小流域模擬數(shù)據(jù)要求高且參數(shù)校正困難[18],模擬結(jié)果誤差可能偏大.因此,本研究氮損失(以N計(jì))模擬主要以氣態(tài)氮損失(NO、N2O、N2和NH3)和淋溶氮損失(NO3--N)為主,模擬結(jié)果如表6所示.由表6可知,B-N處理的氮損失均低于常規(guī)種植CK處理,同種途徑的氮損失在不同處理間的差異具有統(tǒng)計(jì)學(xué)意義(P=0.002 76).模型模擬出的5種氮損失途徑整體顯示出NH3>NO3--N淋失>N2O>NO>N2的特征,其中兩種處理的NH3揮發(fā)損失占總損失的比重均達(dá)到97%以上.CK處理和B-N處理的氮損失模擬結(jié)果顯示,B-N處理的N2O損失通量較CK處理減少0.29 kg/hm2,NO損失通量減少0.22 kg/hm2,NH3揮發(fā)損失通量減少39.59 kg/hm2,NO3--N淋溶損失通量減少0.39 kg/hm2,除N2之外的其余4種途徑的損失通量降低率均大于25%,部分甚至大于40%.CK處理中氮損失模擬結(jié)果占總體氮投入量的24.8%,B-N處理中該數(shù)值為18.0%.B-N處理的氮損失量占總體氮投入量比重小于CK處理,表明B-N處理下的氮肥損失更少,因此B-N處理的氮肥利用率也相對更高.

    表6 菜地氮損失量模擬結(jié)果Fig.6 Simulation results of nitrogen loss in vegetable field

    3 討論

    本研究首次使用DNDC模型模擬施用生物炭后菜地的大蒜和辣椒產(chǎn)量及種植過程的氮損失.結(jié)果表明,DNDC模型能較準(zhǔn)確地模擬施用生物炭并降低施肥量后的作物產(chǎn)量,驗(yàn)證了施用生物炭能促進(jìn)作物對肥料養(yǎng)分利用的特點(diǎn),這與張偉明等[14]研究的生物炭能提高作物氮肥利用率的結(jié)論一致.本研究對研究區(qū)作物產(chǎn)量進(jìn)行模擬,大蒜產(chǎn)量誤差相對較小,辣椒產(chǎn)量誤差相對較大.可能是由于辣椒種植周期為4—9月,期間降水較多,地表養(yǎng)分隨地表徑流運(yùn)移出菜地,導(dǎo)致總體可利用養(yǎng)分減少,進(jìn)而使辣椒生長期內(nèi)受到養(yǎng)分制約.DNDC模型能較好地模擬出氮素在垂直方向上的動態(tài)運(yùn)移,卻不能模擬土壤水分在其他方向的運(yùn)移情況,因此該模型忽略了表層養(yǎng)分隨徑流的流失,產(chǎn)生模擬作物產(chǎn)量高于實(shí)測的情況.李虎等[19]也曾在研究中指出過模型的這一缺點(diǎn),并提出使用其他水文空間分布模型如SWAT、MIKESHE等輔助DNDC模型更準(zhǔn)確地模擬氮素的運(yùn)移及環(huán)境效應(yīng).

    有關(guān)施用生物炭減少氮損失的研究在國內(nèi)外均有報(bào)道.生物炭的獨(dú)特結(jié)構(gòu)使其對NH4+具有較強(qiáng)的吸附能力[20],降低了NH4+轉(zhuǎn)化為NH3的比例,一方面使得土壤中具有較高含量的NH4+供植物生長所需,另一方面也減少了土壤中NH3的含量,減弱了NH3在固-液-氣界面的遷移轉(zhuǎn)化[21],進(jìn)而降低了土壤的NH3揮發(fā)損失.此外生物炭對NH4+的吸附能力減緩了NH4+向NO3-的轉(zhuǎn)化,導(dǎo)致同樣土壤理化性質(zhì)及施氮水平下,施用生物炭菜地的NO3-含量稍低于未施用生物炭菜地的含量.而前文提到淋溶主要以NO3-為主,因此NO3-含量相對較低的施用生物炭的菜地氮流失風(fēng)險低于未施用生物炭的菜地,進(jìn)而表明施用生物炭有效降低了土壤的NO3-淋溶損失,這與高德才等[22]的研究結(jié)果一致.由于生物炭具有多孔結(jié)構(gòu),添加至土壤后能有效降低土壤容重,提高土壤通氣性,改善土壤厭氧狀況[23],進(jìn)而降低了土壤中反硝化細(xì)菌將NO3-還原為NO、N3O、N3等氣體的排放通量.Zhang等[24]的研究中也曾得出施用生物炭減少反硝化過程含氮?dú)怏w排放的結(jié)論.Ding等[25]在使用竹炭對土柱中氮保留和淋失評估的研究中發(fā)現(xiàn)添加竹炭能有效阻止NH4+-N向下運(yùn)移,在20 cm深度處的NH4+-N浸出損失減少了15.2%.Cox等[26]也在研究中發(fā)現(xiàn)生物炭對土壤中的NH4+-N和NH3有很強(qiáng)的吸附作用.由于作物的根系大多位于土壤的淺層,因此添加生物炭能將較多的NH4+-N吸附在土壤淺層,有利于作物對無機(jī)氮的吸收,進(jìn)而提高了肥料的利用率[14].鐘雪梅等[27]和Magrinibair等[28]的研究中均指出生物炭能提高氮肥利用率,這也與本研究中施用生物炭、減少施肥量而作物產(chǎn)量沒有明顯降低的結(jié)果對應(yīng),這表明施用生物炭能有效固持土壤中易被作物吸收的無機(jī)氮,進(jìn)而提高作物吸氮量,達(dá)到了一定的穩(wěn)產(chǎn)效果.

    4 結(jié)論

    (1)DNDC模型能較準(zhǔn)確地模擬出施用生物炭并減少施肥量后的作物產(chǎn)量,表明DNDC模型在模擬施用生物炭的方面具有一定的可行性,其中作物產(chǎn)量模擬誤差在17%以內(nèi).

    (2)DNDC模型的敏感性分析結(jié)果顯示,施肥量、土壤本底SOC、土壤pH值等是影響作物產(chǎn)量和種植過程中氮損失的重要因素.

    (3)施用生物炭并減少施肥量能在保證作物產(chǎn)量的前提下減少11.9%~12.7%的氮、磷投入,表明施用生物炭能有效減少施肥量并達(dá)到穩(wěn)產(chǎn)效果,降低了氮、磷流失風(fēng)險.

    (4)施用生物炭并減少施肥量有效減少了菜地的氮損失,其中除N2損失方式降低率為1.14%以外,NO、N2O、NH3、NO3--N淋失4種氮損失方式降低率達(dá)26.37%~41.26%.

    猜你喜歡
    施肥量作物生物
    生物多樣性
    蔬菜施肥量控制技術(shù)指南
    生物多樣性
    上上生物
    作物遭受霜凍該如何補(bǔ)救
    四種作物 北方種植有前景
    第12話 完美生物
    航空世界(2020年10期)2020-01-19 14:36:20
    內(nèi)生微生物和其在作物管理中的潛在應(yīng)用
    無人機(jī)遙感在作物監(jiān)測中的應(yīng)用與展望
    水果紅薯無土高產(chǎn)栽培的最佳施肥量研究
    精品亚洲成国产av| 成人亚洲欧美一区二区av| 久久99精品国语久久久| 一级a做视频免费观看| 热re99久久国产66热| 99热6这里只有精品| 男人添女人高潮全过程视频| 在线天堂最新版资源| 欧美精品一区二区免费开放| 免费久久久久久久精品成人欧美视频 | 亚洲人成网站在线观看播放| 高清黄色对白视频在线免费看| www.色视频.com| 免费观看无遮挡的男女| 999精品在线视频| 久久狼人影院| 满18在线观看网站| 日韩中文字幕视频在线看片| 国内精品宾馆在线| 精品一区二区三卡| 免费观看在线日韩| 国产伦精品一区二区三区视频9| 精品久久久久久久久av| 一级二级三级毛片免费看| 久久国产亚洲av麻豆专区| 青春草国产在线视频| 蜜桃在线观看..| 成人毛片a级毛片在线播放| 国产精品国产三级国产av玫瑰| 精品卡一卡二卡四卡免费| 亚洲av不卡在线观看| 伊人久久精品亚洲午夜| 26uuu在线亚洲综合色| 视频中文字幕在线观看| 18禁动态无遮挡网站| 麻豆成人av视频| 免费高清在线观看视频在线观看| 黄片播放在线免费| 欧美 日韩 精品 国产| 国产毛片在线视频| 五月伊人婷婷丁香| 亚洲成人手机| 精品亚洲成国产av| 91久久精品电影网| 男女高潮啪啪啪动态图| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 久久精品久久久久久噜噜老黄| 天美传媒精品一区二区| 超碰97精品在线观看| 亚洲欧美成人精品一区二区| 欧美日韩视频精品一区| 大片免费播放器 马上看| 久久久久久久国产电影| 欧美日韩综合久久久久久| 亚洲伊人久久精品综合| 热99国产精品久久久久久7| 丁香六月天网| 高清黄色对白视频在线免费看| 久久人人爽人人爽人人片va| 午夜激情福利司机影院| 一个人看视频在线观看www免费| 欧美亚洲 丝袜 人妻 在线| 伦理电影免费视频| 99热全是精品| 人妻制服诱惑在线中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 久久韩国三级中文字幕| 丝瓜视频免费看黄片| 亚洲中文av在线| 观看av在线不卡| 久久久久久久精品精品| 999精品在线视频| 久久鲁丝午夜福利片| 欧美日韩成人在线一区二区| 久久人人爽av亚洲精品天堂| 自线自在国产av| 久久久国产精品麻豆| 国产国语露脸激情在线看| 亚洲,一卡二卡三卡| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 人妻夜夜爽99麻豆av| 欧美日韩亚洲高清精品| 亚洲欧美清纯卡通| 亚洲不卡免费看| 又黄又爽又刺激的免费视频.| 熟女av电影| 亚洲中文av在线| 伦理电影大哥的女人| 黄片无遮挡物在线观看| 久久人人爽人人爽人人片va| 久久人人爽av亚洲精品天堂| 国产黄片视频在线免费观看| 一本久久精品| 大片免费播放器 马上看| 在线观看www视频免费| 99久久人妻综合| 午夜福利网站1000一区二区三区| 欧美三级亚洲精品| 男女啪啪激烈高潮av片| 视频在线观看一区二区三区| 一级毛片aaaaaa免费看小| 中文字幕人妻丝袜制服| 亚洲熟女精品中文字幕| 黄色怎么调成土黄色| 欧美日韩av久久| 免费黄网站久久成人精品| 97超碰精品成人国产| 国产女主播在线喷水免费视频网站| 国产女主播在线喷水免费视频网站| 亚洲av男天堂| 91久久精品国产一区二区成人| 超碰97精品在线观看| 熟女av电影| 纵有疾风起免费观看全集完整版| 成人国产av品久久久| 不卡视频在线观看欧美| av专区在线播放| 国产在线免费精品| 欧美日韩一区二区视频在线观看视频在线| 只有这里有精品99| 桃花免费在线播放| av卡一久久| 狂野欧美白嫩少妇大欣赏| 久久精品国产鲁丝片午夜精品| 国产不卡av网站在线观看| av.在线天堂| 久久99蜜桃精品久久| av又黄又爽大尺度在线免费看| 亚洲欧洲精品一区二区精品久久久 | 久久久久国产网址| 日本免费在线观看一区| 永久免费av网站大全| 人人妻人人澡人人看| 9色porny在线观看| 简卡轻食公司| 十八禁高潮呻吟视频| 伊人久久精品亚洲午夜| 日韩av在线免费看完整版不卡| 中国三级夫妇交换| 天美传媒精品一区二区| 最近的中文字幕免费完整| 高清在线视频一区二区三区| 99re6热这里在线精品视频| 亚洲精品自拍成人| 亚洲综合色网址| 春色校园在线视频观看| 午夜视频国产福利| 久久精品国产自在天天线| 国产黄片视频在线免费观看| 日韩电影二区| av福利片在线| 午夜激情av网站| 精品人妻熟女毛片av久久网站| 色婷婷av一区二区三区视频| 99热全是精品| 狠狠婷婷综合久久久久久88av| 免费久久久久久久精品成人欧美视频 | 国产成人免费无遮挡视频| 高清在线视频一区二区三区| 日本与韩国留学比较| 91精品国产国语对白视频| 成人综合一区亚洲| 精品久久久久久久久亚洲| videossex国产| 亚洲精品久久成人aⅴ小说 | 丰满饥渴人妻一区二区三| 波野结衣二区三区在线| 91久久精品国产一区二区三区| 中文字幕亚洲精品专区| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 天天影视国产精品| 日本午夜av视频| 一级毛片我不卡| 在线观看免费视频网站a站| 一区二区三区精品91| 91精品国产国语对白视频| 亚洲三级黄色毛片| 91国产中文字幕| 亚洲婷婷狠狠爱综合网| 久久久久久人妻| 人妻 亚洲 视频| 久久久久久久久久久免费av| 成人手机av| 国产精品久久久久久av不卡| 乱码一卡2卡4卡精品| 亚洲欧美一区二区三区黑人 | 国产成人精品久久久久久| 22中文网久久字幕| a级毛片免费高清观看在线播放| 卡戴珊不雅视频在线播放| 亚洲av在线观看美女高潮| 高清午夜精品一区二区三区| 啦啦啦在线观看免费高清www| 免费日韩欧美在线观看| 久久久a久久爽久久v久久| 高清午夜精品一区二区三区| 亚洲精品色激情综合| 国产亚洲精品久久久com| 日韩视频在线欧美| 边亲边吃奶的免费视频| 国产精品免费大片| 日韩 亚洲 欧美在线| 在线观看免费视频网站a站| 天天躁夜夜躁狠狠久久av| 欧美国产精品一级二级三级| 一级黄片播放器| 热99久久久久精品小说推荐| 一区二区三区精品91| 久久国产亚洲av麻豆专区| 熟妇人妻不卡中文字幕| 黑丝袜美女国产一区| 日日撸夜夜添| 成年人午夜在线观看视频| 国产伦精品一区二区三区视频9| 97超视频在线观看视频| 亚洲欧洲日产国产| 免费日韩欧美在线观看| 日韩欧美一区视频在线观看| 亚洲精品自拍成人| 免费av中文字幕在线| 欧美xxxx性猛交bbbb| 夜夜爽夜夜爽视频| 国产免费一级a男人的天堂| 一级a做视频免费观看| 欧美 亚洲 国产 日韩一| 天天影视国产精品| 永久网站在线| 免费不卡的大黄色大毛片视频在线观看| 久久午夜福利片| 欧美一级a爱片免费观看看| 五月开心婷婷网| 国产 一区精品| 亚洲欧洲精品一区二区精品久久久 | av免费观看日本| 九九久久精品国产亚洲av麻豆| 国产成人免费无遮挡视频| 香蕉精品网在线| 欧美激情极品国产一区二区三区 | 99久久中文字幕三级久久日本| 热re99久久国产66热| 亚洲av在线观看美女高潮| 亚洲人成77777在线视频| 亚洲美女视频黄频| 99久久综合免费| 黑人巨大精品欧美一区二区蜜桃 | 国产成人精品久久久久久| 国产精品免费大片| 午夜激情福利司机影院| 老女人水多毛片| 精品少妇黑人巨大在线播放| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 国产一级毛片在线| 午夜日本视频在线| 在现免费观看毛片| 色94色欧美一区二区| 精品一区在线观看国产| 国产日韩欧美亚洲二区| 国产极品粉嫩免费观看在线 | 免费观看在线日韩| videossex国产| 99国产精品免费福利视频| 一边亲一边摸免费视频| 午夜福利视频在线观看免费| 日韩不卡一区二区三区视频在线| 国产精品一二三区在线看| 亚洲av日韩在线播放| 天美传媒精品一区二区| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 熟女av电影| 超色免费av| 久久久久人妻精品一区果冻| 妹子高潮喷水视频| 男人添女人高潮全过程视频| 黑丝袜美女国产一区| 大码成人一级视频| 欧美日韩av久久| 两个人免费观看高清视频| 国产成人精品久久久久久| 国产精品一区二区在线观看99| 黑人高潮一二区| 亚洲人与动物交配视频| 久久国产精品大桥未久av| 精品一区在线观看国产| 中文字幕免费在线视频6| 国产男女超爽视频在线观看| 欧美另类一区| 女的被弄到高潮叫床怎么办| 我的老师免费观看完整版| 精品一区在线观看国产| 亚洲av不卡在线观看| 久久久国产精品麻豆| 一区二区av电影网| 毛片一级片免费看久久久久| 在线亚洲精品国产二区图片欧美 | xxx大片免费视频| 国模一区二区三区四区视频| 一区二区日韩欧美中文字幕 | videosex国产| 777米奇影视久久| 韩国高清视频一区二区三区| 老司机影院毛片| 久久ye,这里只有精品| 国产欧美日韩一区二区三区在线 | 午夜老司机福利剧场| 日韩视频在线欧美| 久久久久精品久久久久真实原创| 老司机影院成人| 亚洲av二区三区四区| 午夜久久久在线观看| 国产精品熟女久久久久浪| av专区在线播放| 亚洲精品,欧美精品| 人妻人人澡人人爽人人| 久久久a久久爽久久v久久| 国产黄频视频在线观看| 精品国产国语对白av| 欧美精品高潮呻吟av久久| 免费观看a级毛片全部| 热99国产精品久久久久久7| 亚洲精品一二三| 精品一区二区免费观看| 少妇熟女欧美另类| 国产白丝娇喘喷水9色精品| 日韩大片免费观看网站| 久久国内精品自在自线图片| 人妻 亚洲 视频| 欧美精品一区二区免费开放| 成人无遮挡网站| 国产色婷婷99| 欧美日韩视频精品一区| 一本大道久久a久久精品| 久久久欧美国产精品| 99热国产这里只有精品6| 国产亚洲精品久久久com| 午夜福利,免费看| 国产一区二区在线观看av| 日韩 亚洲 欧美在线| 成人免费观看视频高清| 国产国拍精品亚洲av在线观看| 亚洲av综合色区一区| 插逼视频在线观看| 精品一区在线观看国产| 国产在线一区二区三区精| 一级毛片黄色毛片免费观看视频| 3wmmmm亚洲av在线观看| 久久精品久久精品一区二区三区| av黄色大香蕉| 日韩av不卡免费在线播放| 美女福利国产在线| 青春草视频在线免费观看| 午夜av观看不卡| 春色校园在线视频观看| 欧美激情国产日韩精品一区| 各种免费的搞黄视频| 在线观看www视频免费| 久久久久久久久久人人人人人人| 国产精品欧美亚洲77777| 亚洲精品乱码久久久v下载方式| 中文天堂在线官网| 久久国产亚洲av麻豆专区| 一区二区日韩欧美中文字幕 | 我的老师免费观看完整版| 国产亚洲精品久久久com| 狠狠婷婷综合久久久久久88av| 婷婷色综合www| 大陆偷拍与自拍| 好男人视频免费观看在线| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| 久久久久国产网址| 好男人视频免费观看在线| 亚洲人成77777在线视频| 国产精品久久久久久av不卡| 插逼视频在线观看| 国产精品嫩草影院av在线观看| 日本-黄色视频高清免费观看| 最近手机中文字幕大全| 纵有疾风起免费观看全集完整版| 日韩 亚洲 欧美在线| 亚洲精品乱码久久久久久按摩| 午夜激情av网站| 亚洲激情五月婷婷啪啪| 啦啦啦啦在线视频资源| 美女主播在线视频| 亚洲国产精品国产精品| 亚洲一级一片aⅴ在线观看| 日本91视频免费播放| 爱豆传媒免费全集在线观看| h视频一区二区三区| 色哟哟·www| 日本91视频免费播放| 日韩av不卡免费在线播放| 国产男女超爽视频在线观看| 亚洲国产av新网站| 日韩 亚洲 欧美在线| 久久久国产精品麻豆| 国产国拍精品亚洲av在线观看| 亚洲国产成人一精品久久久| 亚洲一区二区三区欧美精品| 最近手机中文字幕大全| 欧美日本中文国产一区发布| 国产男女内射视频| 国产日韩欧美亚洲二区| 五月伊人婷婷丁香| 欧美日韩av久久| 久久鲁丝午夜福利片| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 久久99精品国语久久久| 美女国产高潮福利片在线看| 亚洲成色77777| 人人妻人人添人人爽欧美一区卜| 午夜免费观看性视频| 极品人妻少妇av视频| 三上悠亚av全集在线观看| 亚洲欧洲国产日韩| 国产成人精品一,二区| 精品久久国产蜜桃| 久久久久久久久久成人| 亚洲色图综合在线观看| 国产av码专区亚洲av| 我要看黄色一级片免费的| 男的添女的下面高潮视频| 久久精品国产亚洲网站| 久久精品国产自在天天线| 97在线人人人人妻| 午夜久久久在线观看| 久久亚洲国产成人精品v| 久久99热这里只频精品6学生| 亚洲av二区三区四区| 婷婷色综合www| 男的添女的下面高潮视频| 国产乱来视频区| 亚洲精品日韩av片在线观看| 99久久人妻综合| 99久久中文字幕三级久久日本| 免费人成在线观看视频色| 激情五月婷婷亚洲| 国产探花极品一区二区| 好男人视频免费观看在线| 国产免费福利视频在线观看| 三级国产精品片| 日韩伦理黄色片| 老女人水多毛片| 国产亚洲午夜精品一区二区久久| 一个人看视频在线观看www免费| 大片免费播放器 马上看| 制服丝袜香蕉在线| 麻豆精品久久久久久蜜桃| 国产日韩欧美视频二区| 寂寞人妻少妇视频99o| 亚洲av免费高清在线观看| 亚洲av欧美aⅴ国产| 99热全是精品| 99热国产这里只有精品6| 丝袜喷水一区| 黑人猛操日本美女一级片| 十八禁网站网址无遮挡| 最近手机中文字幕大全| 成人毛片a级毛片在线播放| 亚洲欧美一区二区三区国产| 免费高清在线观看视频在线观看| 亚洲美女搞黄在线观看| 亚洲一级一片aⅴ在线观看| 成人午夜精彩视频在线观看| 国精品久久久久久国模美| 日本免费在线观看一区| 一级黄片播放器| 日本午夜av视频| 老女人水多毛片| 男女国产视频网站| 91精品三级在线观看| 午夜免费鲁丝| 国产av精品麻豆| 99九九在线精品视频| 狠狠婷婷综合久久久久久88av| 最近的中文字幕免费完整| 国产熟女午夜一区二区三区 | 69精品国产乱码久久久| 天天影视国产精品| 黄片播放在线免费| 99视频精品全部免费 在线| 亚洲国产精品成人久久小说| 免费观看a级毛片全部| 国产亚洲av片在线观看秒播厂| 一级毛片电影观看| 成人免费观看视频高清| 老司机影院成人| 免费高清在线观看视频在线观看| 蜜桃久久精品国产亚洲av| 又粗又硬又长又爽又黄的视频| 国产精品国产av在线观看| 久久精品久久精品一区二区三区| 亚洲欧美精品自产自拍| 国产精品久久久久久精品古装| 亚洲国产av影院在线观看| 精品一区二区三卡| 18禁在线播放成人免费| 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 亚洲色图 男人天堂 中文字幕 | 伊人久久国产一区二区| 男女国产视频网站| 国产精品一区二区在线观看99| 人妻少妇偷人精品九色| 久久久久精品久久久久真实原创| 免费人成在线观看视频色| 最近中文字幕高清免费大全6| 久久人妻熟女aⅴ| 丰满饥渴人妻一区二区三| 99热网站在线观看| 久久这里有精品视频免费| av一本久久久久| 精品酒店卫生间| 在现免费观看毛片| 欧美日韩成人在线一区二区| 国产日韩欧美亚洲二区| 国产视频首页在线观看| 久久久久久久久久久免费av| 亚洲人与动物交配视频| 亚洲国产av影院在线观看| 少妇的逼水好多| 99久久综合免费| 欧美日韩视频高清一区二区三区二| 久久精品国产鲁丝片午夜精品| 视频区图区小说| 亚洲成人一二三区av| 超色免费av| 国产一级毛片在线| 日韩伦理黄色片| 久久精品国产a三级三级三级| 精品久久久久久电影网| 精品亚洲成a人片在线观看| 黑人高潮一二区| av专区在线播放| av天堂久久9| 一个人看视频在线观看www免费| 国精品久久久久久国模美| 久久久国产精品麻豆| 一区二区三区乱码不卡18| 你懂的网址亚洲精品在线观看| 春色校园在线视频观看| 边亲边吃奶的免费视频| 在线免费观看不下载黄p国产| 视频在线观看一区二区三区| 天堂俺去俺来也www色官网| 街头女战士在线观看网站| 久热久热在线精品观看| 久久久久久伊人网av| 免费看不卡的av| videossex国产| 一区二区三区精品91| 少妇精品久久久久久久| 国产黄片视频在线免费观看| 成人毛片a级毛片在线播放| 内地一区二区视频在线| 狂野欧美激情性xxxx在线观看| 国产男女超爽视频在线观看| 久久99热6这里只有精品| 女人精品久久久久毛片| 91精品三级在线观看| av女优亚洲男人天堂| 26uuu在线亚洲综合色| av视频免费观看在线观看| 丰满少妇做爰视频| 在线观看免费日韩欧美大片 | 最近最新中文字幕免费大全7| 日日摸夜夜添夜夜添av毛片| 亚洲av免费高清在线观看| 精品国产乱码久久久久久小说| 日韩精品免费视频一区二区三区 | 九色亚洲精品在线播放| 亚洲综合色惰| 在线精品无人区一区二区三| 国产精品嫩草影院av在线观看| 又大又黄又爽视频免费| 亚洲欧美成人精品一区二区| 插阴视频在线观看视频| 日本午夜av视频| 午夜视频国产福利| 午夜福利影视在线免费观看| 69精品国产乱码久久久| 亚洲人成网站在线观看播放| 中文字幕精品免费在线观看视频 | 精品久久蜜臀av无| 久久久久久久亚洲中文字幕| 免费看光身美女| 国产日韩一区二区三区精品不卡 | 亚洲国产精品一区二区三区在线| 国产亚洲一区二区精品| 一级片'在线观看视频| 国产精品一区二区在线观看99| 日本黄色片子视频| 少妇的逼好多水| 大陆偷拍与自拍| 91精品国产九色| 精品一区在线观看国产| 天堂中文最新版在线下载| 蜜桃久久精品国产亚洲av| 久久午夜综合久久蜜桃| 亚洲国产精品一区二区三区在线| 久久精品国产亚洲av天美| 青青草视频在线视频观看| 欧美精品人与动牲交sv欧美| videos熟女内射| 日日爽夜夜爽网站| 欧美亚洲 丝袜 人妻 在线| 青春草视频在线免费观看| 99九九在线精品视频| 五月伊人婷婷丁香|