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

    水氮耦合條件下番茄臨界氮濃度模型的建立及氮素營養(yǎng)診斷

    2015-01-27 22:33:35曹紅霞柳美玉劉世和
    關(guān)鍵詞:營養(yǎng)模型

    楊 慧, 曹紅霞, 柳美玉, 劉世和

    (西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室,陜西楊凌 712100)

    水氮耦合條件下番茄臨界氮濃度模型的建立及氮素營養(yǎng)診斷

    楊 慧, 曹紅霞*, 柳美玉, 劉世和

    (西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室,陜西楊凌 712100)

    【目的】臨界氮濃度是指在一定的生長時(shí)期內(nèi)獲得最大生物量時(shí)的最小氮濃度值,具有明確的生物學(xué)意義。探究不同水氮供應(yīng)對(duì)番茄地上部生物量、氮素累積的影響,構(gòu)建臨界氮濃度稀釋曲線模型,并基于氮素吸收和氮營養(yǎng)指數(shù)模型進(jìn)行番茄氮素營養(yǎng)診斷,可為番茄水肥一體化提供一定的理論依據(jù)?!痉椒ā坑?013年在日光溫室內(nèi)進(jìn)行了盆栽試驗(yàn),供試番茄品種為金鵬M6088。設(shè)置3個(gè)灌水量為低水 W1(60%~70%θf)、中水 W2(70%~80%θf)和高水 W3(80%~90%θf),θf為田間持水率;施氮量設(shè)置3個(gè)水平為低氮 N1(N 0.24 g/kg土)、中氮 N2(N 0.36 g/kg土)和高氮 N3(N 0.48 g/kg土),試驗(yàn)采用完全隨機(jī)區(qū)組設(shè)計(jì),共9個(gè)處理,每個(gè)處理重復(fù)15次,研究了不同水氮條件下番茄的地上部生物量、氮素累積及氮濃度的動(dòng)態(tài)變化,構(gòu)建了番茄不同水分條件下的臨界氮濃度稀釋曲線模型?!窘Y(jié)果】番茄地上部生物量、氮累積量隨移栽時(shí)間的動(dòng)態(tài)變化符合Logistic模型,不同水氮供應(yīng)對(duì)番茄地上部生物量理論最大值的影響不同,中水和高水條件下,番茄地上部生物量理論最大值隨著施氮量的增加呈先增加后減小的趨勢(shì);而在低水條件下呈遞增趨勢(shì),說明適量增施氮肥可以減輕干旱對(duì)干物質(zhì)量累積的抑制;番茄地上部生物量快速累積起始日較氮快速累積起始日晚8~17 d,且不同水氮處理番茄地上部生物量最大生長速率、氮累積量最大累積速率均出現(xiàn)在中水中氮(W2N2)處理;在相同的水分條件下,番茄地上部生物量氮濃度隨施氮量的增加而提高,隨生育進(jìn)程的推移呈下降趨勢(shì);氮濃度與地上部生物量之間符合冪指數(shù)關(guān)系,適當(dāng)增大灌水量可以提高植株對(duì)氮的容納能力,并且可以緩解氮濃度隨植株生物增長量下降,使植株穩(wěn)步有序地生長;不同的水氮供應(yīng)對(duì)番茄產(chǎn)量影響顯著,隨著灌水量和施氮量的增加,產(chǎn)量顯著提高,但當(dāng)灌水量和施氮量達(dá)到一定數(shù)量時(shí)產(chǎn)量不僅沒有提高反而隨其增加而降低?!窘Y(jié)論】基于臨界氮濃度構(gòu)建的氮營養(yǎng)指數(shù)、氮吸收模型對(duì)番茄的適宜施氮量診斷結(jié)果一致,均以中水中氮(W2N2)為最佳條件,即當(dāng)灌水量和施肥量分別為62.1 L/plant、15.1 g/plant時(shí),番茄單株產(chǎn)量達(dá)到最大(1602 g),構(gòu)建的模型合理可行。

    番茄; 水氮耦合; 干物質(zhì)累積; 臨界氮濃度; 氮營養(yǎng)指數(shù)

    確定作物不同生長發(fā)育階段的適宜施氮量對(duì)提高氮肥利用效率,解決環(huán)境污染問題具有重要意義[1]。目前,葉綠素含量和遙感技術(shù)被用來診斷作物氮素虧缺,然而這兩種方法在檢測(cè)過量的氮素吸收時(shí)尚有限制[2]。臨界氮濃度被定義為在一定的生長時(shí)期內(nèi)獲得最大生物量時(shí)的最小氮濃度值,確定作物的臨界氮濃度值是作物氮素營養(yǎng)診斷的基本方法之一[3]。國內(nèi)外學(xué)者已針對(duì)牧草[4]、向日葵[5]、包心菜[6]、水稻[7]、紅花[8]、馬鈴薯[9]、玉米[10-11]、棉花[12-13]等作物構(gòu)建了臨界氮濃度稀釋曲線模型,并基于氮素吸收和氮營養(yǎng)指數(shù)模型進(jìn)行了作物需氮量和氮素營養(yǎng)診斷。有學(xué)者對(duì)番茄的氮素營養(yǎng)診斷開展了研究,結(jié)果表明,番茄在生長過程中亦存在臨界、最低和最高氮稀釋曲線[3,14-15],并基于臨界氮濃度模型建立番茄氮素吸收模型、氮素營養(yǎng)指數(shù)模型來診斷番茄氮素營養(yǎng)狀況[3]。前人構(gòu)建的臨界氮濃度模型是在氮素單一因素下建立的,而水氮耦合對(duì)該模型的影響尚不清楚。本研究通過溫室盆栽試驗(yàn),探討了不同水氮條件下番茄地上部生物量及氮素動(dòng)態(tài)累積特征,并構(gòu)建了不同水分條件下番茄臨界氮濃度、氮素吸收、氮素營養(yǎng)指數(shù)模型,探究了上述模型的可靠性,并指導(dǎo)番茄的灌水施肥管理,為番茄水肥一體化提供理論依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)區(qū)概況

    試驗(yàn)于2013年4~7月底在西北農(nóng)林科技大學(xué)旱區(qū)農(nóng)業(yè)水土工程教育部重點(diǎn)實(shí)驗(yàn)室的灌溉試驗(yàn)站進(jìn)行,為溫室盆栽試驗(yàn)。該試驗(yàn)站位于北緯34°18′,東經(jīng)108°40′,海拔高度521 m,多年平均氣溫12.5℃,多年平均蒸發(fā)量1500 mm,年降水量 550~600 mm。供試番茄品種為金鵬M6088。試驗(yàn)用缽下底直徑25 cm,上口直徑30 cm,高度30 cm。每缽裝風(fēng)干土18 kg,控制裝土容重1.3 g/cm3,為防止滯水,缽底部裝河沙1 kg,且每缽垂向裝2根PVC管(直徑2.5 cm、長30 cm)用于灌水。灌水管縱向均勻鉆三排圓孔,外層用網(wǎng)孔直徑1 mm的紗網(wǎng)纏繞兩層,PVC管距缽底部5 cm。供試土壤為重壤土,取自西北農(nóng)林科技大學(xué)節(jié)水灌溉實(shí)驗(yàn)站大田0—20 cm耕層。土壤基本理化性狀為: 田間持水率25.5%、有機(jī)質(zhì)6.18 g/kg、全氮0.81 g/kg、全磷0.42 g/kg、全鉀13.8 g/kg、堿解氮10.93 mg/kg、速效磷4.18 mg/kg、速效鉀102.3 mg/kg。番茄于2013年4月1日移栽定植,2013年7月28日結(jié)束。

    1.2 試驗(yàn)設(shè)計(jì)

    試驗(yàn)處理設(shè)置水分(W)和氮素(N)2個(gè)因素,其中灌水量設(shè)置3個(gè)水平為低水(W1)、中水(W2)和高水(W3);施氮量設(shè)置3個(gè)水平為低氮(N1)、中氮(N2)和高氮(N3),具體灌水量和施氮量見表1。試驗(yàn)采用完全隨機(jī)區(qū)組設(shè)計(jì),共9個(gè)處理,即W1N1、W1N2、W1N3、W2N1、W2N2、W2N3、W3N1、W3N2、W3N3,每個(gè)處理15個(gè)重復(fù),每缽定植1株。當(dāng)番茄進(jìn)入果實(shí)膨大期開始水肥處理直至采摘結(jié)束。

    注(Note):θf為田間持水率Means field capacity.

    試驗(yàn)采用稱重法灌水,當(dāng)土壤含水率降至該處理水分控制下限時(shí),灌水至控制上限。試驗(yàn)用氮、磷、鉀肥分別為尿素(含N量為46%),過磷酸鈣(含P2O5為15%)和硫酸鉀(含K2O為50%);有機(jī)肥料為腐熟的雞糞。磷肥、鉀肥和有機(jī)肥各處理用量相同,分別為P2O50.198 g/kg、K2O 0.315 g/kg 和腐熟的雞糞30 g/kg。磷肥和有機(jī)肥作為基肥一次性施入,氮肥和鉀肥按照基追比1 ∶2施用,追肥在第一穗果膨大期和第二穗果膨大期分別等量隨水灌施。

    1.3 試驗(yàn)測(cè)定項(xiàng)目及方法

    1.3.1 器官干物質(zhì)量和產(chǎn)量的測(cè)定 定苗后分別于苗期、開花期、果實(shí)膨大期、成熟采摘期和拉秧期進(jìn)行破壞性取樣,每次取樣3株。取樣后,分別稱量地上部莖、葉、果鮮重,在105℃下殺青15 min,72℃下烘至恒重,計(jì)算干物質(zhì)量。在果實(shí)采摘期,每次采摘當(dāng)日用電子天平以單株為單位記錄番茄產(chǎn)量。

    1.3.2 植株各器官含氮量的測(cè)定 各處理的干植株樣分器官用小型粉碎機(jī)粉碎過篩,利用濃H2SO4-H2O2法消煮植物樣品,采用AA3型流動(dòng)分析儀測(cè)定消解液中全氮含量。

    各器官氮累積量(g/plant)= 器官氮濃度(g/g)× 器官干物質(zhì)量(g/plant);

    莖、葉、果中氮累積量相加得到植株地上部氮累積量;

    地上部氮濃度(%)= 地上部氮累積量(g/plant)/地上部干物質(zhì)量(g/plant)。

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

    數(shù)據(jù)采用EXCEL和DPS數(shù)據(jù)處理軟件進(jìn)行試驗(yàn)數(shù)據(jù)統(tǒng)計(jì)分析,方差分析使用最小顯著差異法(LSD)進(jìn)行,所有數(shù)據(jù)用OriginPro 8.5作圖。

    2 模型的描述

    2.1 生物量累積過程模擬

    本文采用Logistic模型[16]定量描述番茄地上部生物量和氮累積量的動(dòng)態(tài)變化,以定植后天數(shù)(t)為自變量,番茄各生育期地上部生物量或氮累積量為因變量(y),通過 Logistic方程 y = DMM/(1+ae-bt)[式中DMM為番茄地上部生物量或氮累積量的理論最大值(g/plant)],進(jìn)行模擬。

    將上述 Logistic方程進(jìn)行一階和二階求導(dǎo),得到相應(yīng)生長曲線快速累積期起始時(shí)間T1;快速累積期結(jié)束時(shí)間T2及最大累積速率Vmax.

    (1)

    (2)

    (3)

    2.2 臨界氮濃度稀釋曲線模型

    作物在生長過程中,若植株氮濃度值在臨界氮濃度以下,其生長將受到氮營養(yǎng)的制約;在臨界氮濃度以上,則說明施氮量已超過作物的需求量,作物不受氮素限制;只有植株氮濃度值等于臨界氮濃度時(shí),施氮量最為適宜[12]。按照J(rèn)utes等[17]提出的臨界氮濃度稀釋曲線計(jì)算方法,綜合薛曉萍等[12-13]、王新等[3]關(guān)于棉花、番茄臨界氮濃度稀釋曲線模型的建模思路,構(gòu)建了單株番茄臨界氮濃度稀釋曲線模型,計(jì)算公式如下:

    (4)

    式中,Nc(%)為臨界氮濃度值;a為當(dāng)番茄地上部生物量為 1 kg/plant 時(shí)植株的臨界氮濃度;DWmax為番茄地上部生物量的最大值(kg/plant);b 為決定臨界氮濃度稀釋曲線斜率的統(tǒng)計(jì)學(xué)參數(shù)。

    2.3 氮素吸收

    番茄植株的氮吸收量(Nupt,g/plant)與累積的地上部最大生物量(DW,kg/plant)之間的關(guān)系可用公式(5)表示:

    (5)

    將(4)式代入(5)式得到番茄臨界氮吸收模型:

    (6)

    2.4 氮素營養(yǎng)指數(shù)(NNI)

    為了進(jìn)一步明確作物的氮素營養(yǎng)狀況,采用氮素營養(yǎng)指數(shù)(nitrogen nutrition index,NNI)來反映植株體內(nèi)的氮素情況,可用公式(7)來表示:

    (7)

    式中,NNI為氮素營養(yǎng)指數(shù);Nt為地上部生物量氮濃度的實(shí)測(cè)值(g/100g);Nc為根據(jù)臨界氮濃度稀釋曲線模型求得的在相同的地上部生物量時(shí)的氮濃度值(g/100g)。NNI 可以直觀的反映作物體內(nèi)氮素的營養(yǎng)狀況,NNI=1,氮素營養(yǎng)狀況最為適宜;NNI>1,表現(xiàn)為氮素營養(yǎng)過剩;NNI<1,表現(xiàn)為氮素營養(yǎng)虧缺。

    3 結(jié)果與分析

    3.1 不同水氮條件下番茄地上部生物量及氮素動(dòng)態(tài)累積特征分析

    3.1.1 不同水氮條件下番茄地上部干物質(zhì)累積量 表2為各生育期不同水氮條件下番茄地上部干物質(zhì)的累積情況,可以看出,各水氮處理干物質(zhì)累積量隨著生育期的推進(jìn)而增加。不同水氮處理對(duì)各生育期干物質(zhì)累積量的影響不同,在苗期和花期階段,中氮和高氮處理干物質(zhì)累積量要顯著高于低氮處理,而花期階段中水和高水處理干物質(zhì)累積量也高于低水處理,且隨施氮量增加而增加,說明在營養(yǎng)生長初期,灌水施氮有助于干物質(zhì)快速累積,促進(jìn)植株生長發(fā)育。在果實(shí)膨大期、成熟采摘期和拉秧期,在中水和高水條件下,中氮和高氮處理干物質(zhì)累積量要顯著高于低氮處理,而中氮和高氮處理之間干物質(zhì)累積量無顯著差異;在各氮素條件下,干物質(zhì)累積量在中水條件達(dá)到最大,說明適量灌水施氮有利于植株干物質(zhì)的累積。

    注(Note): 表中數(shù)據(jù)為平均值±標(biāo)準(zhǔn)誤差Data of the table represents average value±standard error;同列數(shù)據(jù)具有相同字母的表示處理間沒有達(dá)到顯著性檢驗(yàn)(P<0.05)Values followed by same letters are not significantly different(P<0.05).

    3.1.2 不同水氮條件下番茄地上部生物量及氮素動(dòng)態(tài)累積特征分析 為了定量研究不同水氮條件下番茄地上部生物量隨生育進(jìn)程的動(dòng)態(tài)累積變化規(guī)律,采用Logistic方程進(jìn)行擬合,得到如下方程(見表3)??梢钥闯?,施氮量相同時(shí),番茄地上部生物量理論最大值隨著灌水量的增加呈先增加后減小的趨勢(shì),說明過量灌水不利于干物質(zhì)量的累積;中水和高水條件下,番茄地上部生物量理論最大值隨著施氮量的增加也呈先增加后減小的趨勢(shì);而在低水條件下呈遞增趨勢(shì),說明適量增施氮肥可以減輕干旱對(duì)干物質(zhì)量累積的抑制。從表中還可以看出,不同水氮條件下番茄地上部生物量快速累積期起始日在定植后的39~43 d(開花期末期)出現(xiàn),且隨著施氮量的增加,快速累積期起始日出現(xiàn)越早,說明較高的施氮量能促進(jìn)植株提早進(jìn)入旺盛的營養(yǎng)生長期。各處理快速累積期持續(xù)時(shí)間為41~54 d,高氮處理干物質(zhì)快速累積持續(xù)時(shí)間要長于中氮和低氮處理。各處理干物質(zhì)最大生長速率出現(xiàn)在W2N2處理,為4.4 g/(plant·d),且隨著灌水量和施氮量的增加呈先上升后下降趨勢(shì)。

    注(Note): T1—為快速累積期起始時(shí)間Starting time in rapid accumulation period;T2—為快速累積期結(jié)束時(shí)Terminating time in rapid accumulation period間;ΔT—為快速累積持續(xù)時(shí)間Duration time in rapid accumulation period;Vmax—為最大累積速率Maximum accumulation rate;**—相關(guān)性在α=0.01水平上差異顯著Significant difference at the 0.01 level.

    由表3中相關(guān)系數(shù)R2可知,不同水氮條件下番茄地上部氮累積量隨生育進(jìn)程的動(dòng)態(tài)變化也符合Logistic生長曲線,可以看出,相同施氮量下,隨著灌水量的增加,番茄地上部氮累積量理論最大值呈先增加后減小的趨勢(shì);相同灌水量下,番茄地上部氮累積量理論最大值隨著施氮量的增加也呈先增加后減小的趨勢(shì),這與不同水氮條件下地上部生物量的變化規(guī)律基本一致。從表中還可以看出,不同水氮條件下番茄地上部氮累積量最大累積速率出現(xiàn)在W2N2處理,為0.098 g/(plant·d),且隨著灌水量和施氮量的增加呈先上升后下降趨勢(shì)??焖倮鄯e期持續(xù)時(shí)間為35~50 d,快速累積期起始日在定植后的23~34 d,較生物量累積早8~17 d,這與宋海星等[18]、王新等[19]的研究結(jié)果一致,說明生物量的增長是以充足的養(yǎng)分吸收為基礎(chǔ)的,且隨著施氮量的增加,地上部氮素快速累積期起始日出現(xiàn)越早。

    3.2 不同水氮條件下番茄臨界氮濃度稀釋曲線模型的建立

    3.2.1 不同水氮條件下番茄氮濃度的動(dòng)態(tài)變化 圖1為不同水氮條件下番茄地上部氮濃度的動(dòng)態(tài)變化,可以看出,不同水氮條件下番茄地上部氮濃度值均表現(xiàn)為隨移栽天數(shù)的推移而降低,即其氮濃度值存在稀釋現(xiàn)象,并且在相同的水分條件下,氮濃度值隨施氮量的增加而增大,說明施氮可以增大植株對(duì)氮素養(yǎng)分的吸收。

    [注(Note): 同一采樣時(shí)期,不同字母表示氮肥處理間氮質(zhì)量分?jǐn)?shù)差異達(dá)0.05水平At the same sampling date, different letters indicate significant difference between N concentrations in different nitrogen treatments at 0.05 level.]

    3.2.2 不同水氮條件下番茄臨界氮濃度稀釋模型和氮素吸收模型的建立 根據(jù)公式(4)、(6)得到不同水氮條件下番茄地上部臨界氮濃度稀釋模型和氮素吸收模型,模型參數(shù)見表4,由相關(guān)系數(shù)可以看出,模型擬合結(jié)果較好。通過番茄地上部生物量與氮濃度值的分析結(jié)果表明,在各水分條件下,同樣的生物量其氮濃度值有很大的變異性,利用各取樣日氮濃度的最大、最小值(%Nmax、%Nmin),得到2個(gè)氮稀釋邊界模型,模型參數(shù)見表4。

    由表4中各水分條件下番茄地上部臨界、最高、最低氮濃度稀釋曲線參數(shù)可以看出,不同水分條件下參數(shù)值a不同,且隨著灌水量的增加呈先增大后減小的趨勢(shì),說明適度的灌水可以增大植株臨界氮濃度值,增加植株對(duì)氮的容納能力,使其對(duì)應(yīng)的需氮量也較高。原因可能是在低水條件下,番茄植株受水分的制約而長勢(shì)減弱,導(dǎo)致氮吸收能力降低,植株內(nèi)含氮量下降;而在高水條件下植株需氮量雖高,但植株含水量也增大,故臨界氮濃度值反而降低。其次,不同水分條件下各稀釋曲線的斜率(參數(shù)b)亦不相同,且中水處理要小于低水和高水處理,表明中水條件下番茄地上部氮濃度隨干物重的增加降低緩慢,說明適宜的灌水可以緩解氮濃度隨植株生物增長量下降,使植株穩(wěn)步有序地生長。

    3.3 基于臨界氮濃度的不同水氮條件下番茄氮素營養(yǎng)狀況分析

    3.3.1 基于氮素吸收模型的番茄適宜水氮條件分析 利用表4中建立的氮素吸收模型,可獲得不同水氮條件下番茄地上部生物量與氮吸收量的關(guān)系圖(圖2),可以看出,在相同的水分條件下,隨著施氮量的增加各氮素水平氮吸收量呈增加的趨勢(shì),而臨界氮吸收量曲線始終在氮素水平為N2~N3之間,即實(shí)際施氮量14.087~18.783 g/plant。采用標(biāo)準(zhǔn)誤差(RMSE)對(duì)圖2中各水氮條件下實(shí)測(cè)氮累積量和臨界氮累積量之間的偏離程度進(jìn)行統(tǒng)計(jì)分析,低氮、中氮、高氮水平下動(dòng)態(tài)氮累積曲線與臨界氮累積曲線的RMSE,低水條件下分別為1.087、0.316和0.154 g/plant;中水條件下分別為0.730、0.119和0.351 g/plant;高水條件下分別為0.694 g/plant、0.284 g/plant、0.163 g/plant,即在中水條件下中氮水平動(dòng)態(tài)氮累積曲線與臨界氮累積曲線的RMSE最小,說明在中水條件下中氮為較適宜的施氮量,而在低水和高水條件下高氮為較適宜的施氮量,原因可能是低水條件下番茄植株受水分制約,而適量增施氮肥能夠緩解干旱對(duì)植株生長的抑制;高水條件下可能由于水分促增長效應(yīng),使植株大幅度增長,需要吸收更多氮素來滿足作物生長對(duì)養(yǎng)分的需求。

    注(Note): **表示相關(guān)性在α=0.01水平上差異顯著 Significant difference at the 0.01 level.

    3.3.2 基于氮營養(yǎng)指數(shù)的番茄適宜水氮條件分析 根據(jù)公式(7)計(jì)算得到不同水氮條件下番茄氮素營養(yǎng)指數(shù)(NNI)隨移栽天數(shù)的動(dòng)態(tài)變化。由圖3可以看出,不同水肥條件下的NNI變化趨勢(shì)相同,均表現(xiàn)出一定的波動(dòng)性,且隨施氮量增加NNI值也變大;施氮量相同時(shí),不同水分處理下NNI變化趨勢(shì)也相同,且表現(xiàn)為W2>W1>W3。在低水條件下,低氮和中氮處理的NNI基本均小于1,表明氮素養(yǎng)分供應(yīng)不足,出現(xiàn)了氮虧缺;而高氮處理下NNI均大于1,表明氮素養(yǎng)分供應(yīng)充足,為氮素營養(yǎng)過剩。在高水條件下,各施氮處理NNI基本均小于1,說明高水條件下,植株出現(xiàn)缺氮的現(xiàn)象。在中水條件下,低氮處理NNI均小于1,表明氮素養(yǎng)分供應(yīng)不足;高氮處理NNI均大于1,說明氮素養(yǎng)分供應(yīng)充足,而中氮處理的NNI始終最接近1或在1附近變化,表明此時(shí)的水氮條件較為適宜。

    3.3.3 不同水氮供應(yīng)與番茄產(chǎn)量的關(guān)系 為了驗(yàn)證以上氮素吸收模型和氮營養(yǎng)指數(shù)模型對(duì)番茄水氮狀況診斷的合理性,對(duì)不同的水氮供應(yīng)與單株產(chǎn)量之間關(guān)系作回歸分析,得如下二元二次方程:

    0.5x1x2-4321.5

    (8)

    式中,Y為番茄單株產(chǎn)量(g/plant);x1為各處理的灌水量(L/plant);x2為各處理的施氮量(g/plant)。

    [注(Note): 同一采樣時(shí)期,不同字母表示氮肥處理間氮營養(yǎng)指數(shù)差異達(dá)0.05水平At the same sampling date, different letters indicate significant difference betweenNNIin different nitrogen treatments at 0.05 level.]

    F檢驗(yàn)結(jié)果表明,(8)式達(dá)到5%顯著水平(F=341.4>F0.05(5,3)=9.01),說明該式能較好地表達(dá)番茄水氮供應(yīng)與單株產(chǎn)量之間的關(guān)系。將(8)式繪圖表達(dá)如圖4。對(duì)(8)式求偏導(dǎo),當(dāng)灌水量和施肥量分別為62.1 L/plant、15.1 g/plant時(shí),番茄單株產(chǎn)量達(dá)到最大,為1602 g/plant,這與基于臨界氮濃度構(gòu)建的氮營養(yǎng)指數(shù)、氮吸收模型對(duì)番茄的水氮狀況診斷結(jié)果一致,均以中水中氮為最佳條件。從圖4可以看出,隨著灌水量和施氮量的增加,產(chǎn)量顯著提高,但是當(dāng)灌水量和施氮量達(dá)到一定數(shù)量時(shí)產(chǎn)量不僅沒有提高反而隨著其增加而降低,說明只有適宜的水分和氮肥才能有助于高產(chǎn)。

    4 結(jié)論

    本文通過對(duì)不同水氮條件下番茄的地上部生物量、氮素累積及氮濃度的動(dòng)態(tài)變化進(jìn)行分析,表明番茄地上部生物量增長、氮吸收累積均受灌水量和施氮量的影響,且其隨移栽時(shí)間的動(dòng)態(tài)變化符合Logistic模型,在中水中氮條件下可以獲得較高的地上部干物質(zhì)累積量,從而獲得較大的養(yǎng)分吸收速率,地上部生物量快速累積起始日較氮快速累積起始日晚8~17 d。在相同的水分條件下,番茄地上部氮濃度隨施氮量的增加而提高,隨生育進(jìn)程的推移呈下降趨勢(shì)。本文還在3種水分條件下分別構(gòu)建了番茄地上部生物量的臨界、最低和最高氮濃度稀釋曲線模型,結(jié)果表明氮濃度與地上部最大生物量之間符合冪指數(shù)關(guān)系,適當(dāng)增大灌水量可以提高植株對(duì)氮的容納能力。根據(jù)番茄單株產(chǎn)量與水氮供應(yīng)的關(guān)系,由擬合曲面得到當(dāng)灌水量和施肥量分別為62.1 L/plant、15.1 g/plant時(shí),番茄單株產(chǎn)量達(dá)到最大,為1602 g/plant?;谂R界氮濃度構(gòu)建的氮營養(yǎng)指數(shù)、氮吸收模型對(duì)番茄的氮素營養(yǎng)狀況診斷均以中水中氮最優(yōu),這與上述試驗(yàn)結(jié)論相符,說明構(gòu)建的模型合理可行。

    [1] Galloway J N, Townsend A R, Erisman J Wetal. Transformation of the nitrogen cycle: recent trends, questions, and potential solutions[J]. Science, 2008,(320): 889-892.

    [2] Zhao B. Determining of a critical dilution curve for plant nitrogen[J]. Field Crops Research, 2014,(160): 64-72.

    [3] 王新,馬富裕,刁明,等.滴灌番茄臨界氮濃度、氮素吸收和氮營養(yǎng)指數(shù)模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013, 29(18): 107-116. Wang X, Ma F Y, Diao Metal. Simulation of critical nitrogen concentration, nitrogen uptake and nitrogen nutrition index of processing tomato with drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering, 2013, 29(18): 107-116.

    [4] Lemaire G, Gastal F. N uptake and distribution in plant canopies[M]. Berlin Heidelberg Springer-Verlag, 1997. 3-43.

    [5] Debaeke P, Oosterom E J, Justes Eetal. A species-specific criti- cal nitrogen dilution curve for sunflower[J]. Field Crops Research, 2012,(136): 76-84.

    [6] Ekbladh G, Witter E. Determination of the critical nitrogen conce- ntration of white cabbage[J]. European Journal of Agronomy, 2010,(33): 276-284.

    [7] Ata-Ul-Karim S T, Yao X, Liu X Jetal. Development of critical nitrogen dilution curve of Japonica rice in Yangtze River Reaches[J]. Field Crops Research, 2013,(149): 149-158.

    [8] Dordas C A, Sioulas C. Dry matter and nitrogen accumulation, partitioning, and retranslocation in safflower(CarthamustinctoriusL.)as affected by nitrogen fertilization[J]. Field Crops Research, 2009, 110(1): 35-43.

    [9] Giletto C M, Echeverría H E. Critical nitrogen dilution curve for processing potato in Argentinean humid pampas[J]. American Journal of Potato Research, 2012, 89(2): 102-110.

    [10] Plenet D, Lemaire G. Relationships between dynamics of nitrog- en uptake and dry matter accumulation in maize crops. Determin- ation of critical N concentration[J]. Plant and Soil, 2000, 216(1-2): 65-82.

    [11] 梁效貴, 張經(jīng)廷, 周麗麗, 等.華北地區(qū)夏玉米臨界氮稀釋曲線和氮營養(yǎng)指數(shù)研究[J]. 作物學(xué)報(bào), 2013, 39(2): 292-299. Liang X G, Zhang J T, Zhou L Letal. Critical nitrogen dilution curve and nitrogen nutrition index for summer maize in north China plain[J]. Acta Agronomica Sinica, 2013, 39(2): 292-299.

    [12] 薛曉萍,沙奕卓,郭文琦,等.棉花蕾花鈴生物量、氮累積特征及臨界氮濃度稀釋模型[J]. 生態(tài)學(xué)報(bào), 2008, (12): 423-430. Xue X P, Sha Y Z, Guo W Qetal. Accumulation characteristics of biomass and nitrogen and critical nitrogen concentration dilution model of cotton reproductive organ[J]. Acta Ecologica Sinica, 2008, (12): 423-430.

    [13] 薛曉萍,王建國,郭文琦,等.棉花花后果枝葉生物量和氮累積特征及臨界氮濃度稀釋模型的研究[J]. 作物學(xué)報(bào), 2007, (4): 151-158. Xue X P, Wang J G, Guo W Qetal. Accumulation characters of biomass and nitrogen and critical nitrogen concentration dilution model of cotton fruit-branch leaf after flowering[J]. Acta Agronomica Sinca, 2007, (4): 151-158.

    [14] Tei F, Benincasa P, Guiducci M. Critical nitrogen concentration in processing tomato[J]. European Journal of Agronomy, 2002, 18(1): 45-55.

    [15] Hartz T K, Bottoms T G. Nitrogen requirements of drip-irrigated processing tomatoes[J]. HortScience, 2009, 44(7): 1988-1993.

    [16] 彭世彰.農(nóng)業(yè)高效節(jié)水灌溉理論與模式[M]. 北京: 科學(xué)出版社, 2009. Peng S Z.Theory and model on agricultural efficient water saving irrigation[M].Beijing: Science Press,2009.

    [17] Justes E, Mary B, Meynard J M. Determination of a critical nitrogen dilution curve for winter wheat crops[J]. Annals of Botany, 1994, 74(4): 397-407.

    [18] 宋海星,李生秀.玉米生長量、養(yǎng)分吸收量及氮肥利用率的動(dòng)態(tài)變化[J]. 中國農(nóng)業(yè)科學(xué), 2003, 3(61): 71-76. Song H X, Li S X. Dynamics of nutrient accumulation in maize plants under different water and N supply conditions[J]. Scientia Agricultura Sinica, 2003, 3(61): 71-76.

    [19] 王新, 馬富裕, 刁明, 等.不同施氮水平下加工番茄植株生長和氮素積累與利用率的動(dòng)態(tài)模擬[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(4): 1043-1050. Wang X, Ma F Y, Diao Metal. Dynamics simulation on plant growth, N accumulation and utilization of processing tomato at different N fertilization rates[J]. Chinese Journal of Applied Ecology, 2014, 25(4): 1043-1050.

    Simulation of critical nitrogen concentration and nitrogen nutrition index of tomato under different water and nitrogen conditions

    YANG Hui, CAO Hong-xia*, LIU Mei-yu, LIU Shi-he

    (KeyLaboratoryofAgriculturalSoilandWaterEngineeringinAridandSemiaridAreaofMinistryofEducation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

    【Objectives】 The critical nitrogen(N)concentration in plant aboveground biomass is defined as the minimum N concentration required for maximum plant growth. This study investigated the effects of different water and nitrogen supply on tomato aboveground biomass, nitrogen accumulation, and drew a critical N concentration dilution curve. The N status of tomato plant was analyzed based on a model of N uptake and nitrogen nutrition index(NNI), which provided a theoretical basis for optimal water and nitrogen management. 【Methods】A pot experiment was conducted in greenhouse of the Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Area of Ministry of Education in Northwest Agriculture and Forestry University in 2013. Cultivar of tomato(M6088)was used in this experiment. The treatment comprised three levels of irrigation(W1, 60%-70%θf; W2, 70%-80%θf; W3, 80%-90%θf),θfis the field capacity, and three levels of nitrogen(N1, N 0.24 g/kg; N2, N 0.36 g/kg; N3, N 0.48 g/kg). For determining the critical N concentration dilution curves under different water conditions, the treatments were replicated fifteen times in random complete block designs to examine the dynamic changes of tomato aboveground biomass and nitrogen accumulation under different water and nitrogen conditions. 【Results】 The aboveground biomass and N accumulations presented a Logistic curve over time. Different water and nitrogen supply had different effects on maximum theoretical value of tomato aboveground biomass: the maximum theoretical value of tomato aboveground biomass increased firstly and decreased with the increase of nitrogen rate under two levels of irrigation(W2, W3). It also increased with the increase of nitrogen rate under the level of irrigation(W1), which indicated that moderate nitrogen supply could enhance the inhibiting effect of drought on aboveground biomass accumulation of tomato. The beginning time of fast accumulation period for nitrogen was 8-17 days earlier than those for biomass, the maximum accumulation rates of tomato aboveground biomass and nitrogen were both found in W2N2 treatment. Under the same water supply condition, the nitrogen concentration of tomato aboveground biomass increased with the improving of applied N rates, and decreased in the growing process. The relationship between the aboveground biomass and N concentration could be described by the power equation, appropriate irrigation could improve the capacity of plant for nitrogen absorption and relieved the decline of nitrogen concentration with the aboveground biomass growth to ensure a steady and orderly growth of tomato. The yield was significantly affected by water and nitrogen supply, appropriate condition of water and nitrogen achieved maximum yield. 【Conclusions】 Based on the model of nitrogen nutrition(NNI)and the model of N uptake, the W2N2 treatment was the optimal option with irrigation amount of 62.1 L/plant, nitrogen rate of 15.1 g/plant, and the highest yield was 1602 g/plant. Thus, the models built in this study were reasonable and feasible for the research objectives.

    tomato; coupling of water and nitrogen; biomass accumulation; critical nitrogen concentration; nitrogen nutrition index

    2014-05-24 接受日期: 2014-10-16 網(wǎng)絡(luò)出版日期: 2015-05-21

    水利部公益性行業(yè)科研專項(xiàng)(201001061); 陜西省自然科學(xué)基金項(xiàng)目(2012JM3004); 中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(QN2011022)資助。

    楊慧(1989—),女,青海西寧人,碩士研究生,主要從事農(nóng)業(yè)節(jié)水理論研究。E-mail: yh2438130@163.com *通信作者E-mail: chx662002@163.com

    S641.2; S606

    A

    1008-505X(2015)05-1234-09

    猜你喜歡
    營養(yǎng)模型
    一半模型
    夏季給豬吃點(diǎn)“苦”營養(yǎng)又解暑
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    這么喝牛奶營養(yǎng)又健康
    營養(yǎng)Q&A
    幸福(2018年33期)2018-12-05 05:22:46
    眼睛需要營養(yǎng)嗎
    蔬菜與營養(yǎng)
    3D打印中的模型分割與打包
    春天最好的“營養(yǎng)菜”
    海峽姐妹(2016年5期)2016-02-27 15:20:20
    欧美精品一区二区大全| 久久精品成人免费网站| 久久这里只有精品19| 欧美日韩视频高清一区二区三区二| 免费日韩欧美在线观看| 极品少妇高潮喷水抽搐| 女性生殖器流出的白浆| 成人手机av| 多毛熟女@视频| 王馨瑶露胸无遮挡在线观看| 精品人妻1区二区| 欧美激情 高清一区二区三区| 天天操日日干夜夜撸| 欧美黄色淫秽网站| 欧美日韩国产mv在线观看视频| 中文字幕亚洲精品专区| 久久久久精品人妻al黑| 观看av在线不卡| 中文字幕最新亚洲高清| 777久久人妻少妇嫩草av网站| 日本91视频免费播放| 看十八女毛片水多多多| 欧美日韩福利视频一区二区| 亚洲成人手机| videos熟女内射| 人人妻人人澡人人爽人人夜夜| 国产熟女午夜一区二区三区| 日日夜夜操网爽| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区四区激情视频| 丝瓜视频免费看黄片| 久久鲁丝午夜福利片| 国产色视频综合| 久久ye,这里只有精品| 国产精品一区二区精品视频观看| 国产一卡二卡三卡精品| 母亲3免费完整高清在线观看| 国产精品人妻久久久影院| 一本久久精品| 少妇裸体淫交视频免费看高清 | 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩综合在线一区二区| 丁香六月欧美| 国产日韩欧美亚洲二区| 亚洲九九香蕉| 久久免费观看电影| 一本—道久久a久久精品蜜桃钙片| 啦啦啦在线观看免费高清www| 热re99久久国产66热| 悠悠久久av| 老汉色av国产亚洲站长工具| 少妇人妻 视频| 99国产精品99久久久久| 黑人巨大精品欧美一区二区蜜桃| 国产精品av久久久久免费| 老司机午夜十八禁免费视频| av不卡在线播放| 午夜福利视频在线观看免费| 国产精品久久久久久人妻精品电影 | 欧美黄色片欧美黄色片| 99热国产这里只有精品6| 看免费av毛片| 国产亚洲精品第一综合不卡| 免费观看人在逋| 男女边摸边吃奶| 电影成人av| 精品亚洲成a人片在线观看| av天堂在线播放| 欧美人与善性xxx| 日本av免费视频播放| 欧美久久黑人一区二区| 天天躁夜夜躁狠狠躁躁| 你懂的网址亚洲精品在线观看| 国产成人免费观看mmmm| 日韩免费高清中文字幕av| 免费黄频网站在线观看国产| 久久天躁狠狠躁夜夜2o2o | 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版| 久久午夜综合久久蜜桃| 午夜视频精品福利| 精品国产一区二区三区久久久樱花| 色婷婷av一区二区三区视频| 亚洲午夜精品一区,二区,三区| 男人添女人高潮全过程视频| 大香蕉久久成人网| 久久青草综合色| 999久久久国产精品视频| 亚洲成色77777| 精品久久久久久电影网| 婷婷色综合www| 欧美在线黄色| 国产成人影院久久av| 欧美国产精品一级二级三级| 久久青草综合色| 婷婷色综合www| 亚洲少妇的诱惑av| 老鸭窝网址在线观看| 欧美另类一区| 中文字幕制服av| 黑丝袜美女国产一区| 十分钟在线观看高清视频www| 99九九在线精品视频| 欧美精品高潮呻吟av久久| 国产免费福利视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 一边摸一边做爽爽视频免费| 美女高潮到喷水免费观看| 日韩制服骚丝袜av| 日本欧美国产在线视频| 亚洲成人国产一区在线观看 | 亚洲国产精品国产精品| 日韩人妻精品一区2区三区| av天堂久久9| 欧美亚洲日本最大视频资源| 午夜福利免费观看在线| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 午夜福利,免费看| 国产欧美日韩综合在线一区二区| 欧美人与性动交α欧美软件| 欧美黑人精品巨大| 中文字幕人妻熟女乱码| www日本在线高清视频| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 日日摸夜夜添夜夜爱| 亚洲图色成人| 国产成人精品无人区| 亚洲精品国产色婷婷电影| 性色av一级| 看免费成人av毛片| 久久久久久亚洲精品国产蜜桃av| 777久久人妻少妇嫩草av网站| 午夜福利,免费看| 国产欧美日韩综合在线一区二区| 9热在线视频观看99| 国产伦理片在线播放av一区| 桃花免费在线播放| 久久久久国产精品人妻一区二区| 日韩熟女老妇一区二区性免费视频| 日韩大片免费观看网站| 久久久精品区二区三区| 叶爱在线成人免费视频播放| 国产精品 欧美亚洲| 国产精品久久久久久精品古装| 成人免费观看视频高清| 精品欧美一区二区三区在线| 欧美 日韩 精品 国产| 夫妻性生交免费视频一级片| 久久久久精品国产欧美久久久 | tube8黄色片| 亚洲激情五月婷婷啪啪| 国产av国产精品国产| 91成人精品电影| 国产精品久久久久久人妻精品电影 | 国产精品久久久久久人妻精品电影 | 视频在线观看一区二区三区| 咕卡用的链子| 一二三四社区在线视频社区8| 精品一区二区三区四区五区乱码 | 午夜福利免费观看在线| 国产亚洲午夜精品一区二区久久| 亚洲专区中文字幕在线| 亚洲国产欧美一区二区综合| 精品国产国语对白av| 一区二区av电影网| 精品少妇内射三级| 国产精品国产av在线观看| 免费在线观看黄色视频的| 国产精品熟女久久久久浪| 久久久欧美国产精品| 嫁个100分男人电影在线观看 | 高清黄色对白视频在线免费看| 黑人猛操日本美女一级片| 欧美 日韩 精品 国产| 黄频高清免费视频| 一级毛片女人18水好多 | www.熟女人妻精品国产| 午夜免费成人在线视频| 欧美日韩成人在线一区二区| 老司机亚洲免费影院| 国产午夜精品一二区理论片| 99热网站在线观看| 精品福利观看| 久久国产精品大桥未久av| 亚洲国产中文字幕在线视频| 婷婷成人精品国产| 国产女主播在线喷水免费视频网站| 色视频在线一区二区三区| 涩涩av久久男人的天堂| 波多野结衣一区麻豆| 国产在视频线精品| 国产在线视频一区二区| 多毛熟女@视频| 成年av动漫网址| 国产成人精品无人区| 久久久久久久精品精品| 久热这里只有精品99| 观看av在线不卡| 亚洲精品日本国产第一区| 一本一本久久a久久精品综合妖精| 欧美在线一区亚洲| 亚洲国产欧美在线一区| 久久 成人 亚洲| 国产视频一区二区在线看| 丁香六月天网| 观看av在线不卡| 男人爽女人下面视频在线观看| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 性色av一级| 免费女性裸体啪啪无遮挡网站| 中国国产av一级| 97在线人人人人妻| 69精品国产乱码久久久| 9热在线视频观看99| 女人久久www免费人成看片| √禁漫天堂资源中文www| 纯流量卡能插随身wifi吗| 亚洲一码二码三码区别大吗| 99国产精品99久久久久| 高清视频免费观看一区二区| 可以免费在线观看a视频的电影网站| 建设人人有责人人尽责人人享有的| 波野结衣二区三区在线| 久久影院123| 极品少妇高潮喷水抽搐| 操出白浆在线播放| 秋霞在线观看毛片| 久久99一区二区三区| av在线播放精品| 久久狼人影院| 欧美黄色淫秽网站| 日韩制服丝袜自拍偷拍| 久久天躁狠狠躁夜夜2o2o | 欧美精品一区二区大全| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀 | 久久av网站| 爱豆传媒免费全集在线观看| 亚洲 国产 在线| 天天躁夜夜躁狠狠久久av| 美女福利国产在线| 欧美黑人欧美精品刺激| 丰满人妻熟妇乱又伦精品不卡| 极品少妇高潮喷水抽搐| 麻豆国产av国片精品| 欧美日韩精品网址| 国产成人精品久久久久久| 免费一级毛片在线播放高清视频 | 亚洲精品一区蜜桃| 中文字幕人妻熟女乱码| 啦啦啦 在线观看视频| 国产熟女午夜一区二区三区| 精品少妇一区二区三区视频日本电影| 男女下面插进去视频免费观看| 亚洲色图综合在线观看| 免费看不卡的av| 午夜老司机福利片| 午夜视频精品福利| 男人操女人黄网站| 午夜两性在线视频| 久久天躁狠狠躁夜夜2o2o | 欧美人与善性xxx| 国产男女内射视频| 免费黄频网站在线观看国产| 国产一区二区在线观看av| 欧美成人午夜精品| 国产成人精品无人区| 国产免费一区二区三区四区乱码| 国产精品.久久久| 久久影院123| 久久狼人影院| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久国产电影| 婷婷色综合大香蕉| 性色av一级| 亚洲黑人精品在线| 亚洲一区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 久久 成人 亚洲| 国产伦理片在线播放av一区| av天堂久久9| 精品国产一区二区三区四区第35| 18禁裸乳无遮挡动漫免费视频| 成年av动漫网址| 香蕉国产在线看| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 亚洲色图综合在线观看| 欧美变态另类bdsm刘玥| 青草久久国产| 午夜福利视频精品| 国产日韩欧美亚洲二区| 免费女性裸体啪啪无遮挡网站| av又黄又爽大尺度在线免费看| 老司机影院成人| 精品人妻熟女毛片av久久网站| 国产成人欧美在线观看 | 亚洲av综合色区一区| 女性生殖器流出的白浆| 免费久久久久久久精品成人欧美视频| 9191精品国产免费久久| 亚洲成人免费电影在线观看 | 国产成人91sexporn| 亚洲精品国产av成人精品| 久久久久精品人妻al黑| 制服人妻中文乱码| 热99久久久久精品小说推荐| 2018国产大陆天天弄谢| 天天影视国产精品| 啦啦啦在线观看免费高清www| 丁香六月天网| 免费不卡黄色视频| 国产日韩欧美亚洲二区| 99精品久久久久人妻精品| 亚洲三区欧美一区| 国产精品成人在线| 另类精品久久| 午夜av观看不卡| 亚洲欧美日韩另类电影网站| 日韩一本色道免费dvd| av电影中文网址| 国产精品.久久久| 亚洲 欧美一区二区三区| 亚洲成人国产一区在线观看 | www.av在线官网国产| 亚洲av电影在线观看一区二区三区| 两个人免费观看高清视频| av网站在线播放免费| 亚洲欧美色中文字幕在线| 国产一区二区在线观看av| 三上悠亚av全集在线观看| 久久九九热精品免费| 中文乱码字字幕精品一区二区三区| 亚洲成色77777| 高清不卡的av网站| 99国产精品免费福利视频| 久久国产精品人妻蜜桃| 久久av网站| 国产精品三级大全| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 啦啦啦在线免费观看视频4| 波多野结衣av一区二区av| 青草久久国产| 亚洲欧美清纯卡通| 另类亚洲欧美激情| 精品久久久久久电影网| 亚洲 国产 在线| 日韩伦理黄色片| 晚上一个人看的免费电影| 亚洲熟女精品中文字幕| 国产欧美日韩一区二区三区在线| 午夜老司机福利片| 女人高潮潮喷娇喘18禁视频| 最新在线观看一区二区三区 | 国产色视频综合| 亚洲三区欧美一区| 国产精品成人在线| 少妇粗大呻吟视频| 日本一区二区免费在线视频| 亚洲伊人色综图| 老司机午夜十八禁免费视频| 一个人免费看片子| 国产亚洲欧美精品永久| 精品一品国产午夜福利视频| 伦理电影免费视频| 欧美人与善性xxx| 色视频在线一区二区三区| 欧美人与善性xxx| 我的亚洲天堂| 亚洲国产日韩一区二区| 91精品国产国语对白视频| 亚洲中文av在线| 91精品伊人久久大香线蕉| 色94色欧美一区二区| 精品卡一卡二卡四卡免费| 老司机影院毛片| 亚洲国产最新在线播放| 天堂中文最新版在线下载| 一级毛片我不卡| 在线精品无人区一区二区三| 色婷婷av一区二区三区视频| 成人国语在线视频| 午夜激情av网站| 国产成人精品无人区| 国产激情久久老熟女| 成年av动漫网址| 亚洲美女黄色视频免费看| 免费在线观看视频国产中文字幕亚洲 | 午夜91福利影院| 亚洲成人国产一区在线观看 | 少妇精品久久久久久久| 一区在线观看完整版| 免费在线观看视频国产中文字幕亚洲 | 精品人妻1区二区| 汤姆久久久久久久影院中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲第一青青草原| 精品久久蜜臀av无| 婷婷成人精品国产| 高清黄色对白视频在线免费看| 亚洲,一卡二卡三卡| 一区福利在线观看| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 国产女主播在线喷水免费视频网站| 国产成人免费无遮挡视频| av欧美777| 日本一区二区免费在线视频| 亚洲中文字幕日韩| 51午夜福利影视在线观看| 亚洲 欧美一区二区三区| 久久久欧美国产精品| 亚洲成国产人片在线观看| 婷婷成人精品国产| 亚洲av欧美aⅴ国产| 国产成人免费无遮挡视频| 两个人免费观看高清视频| 成人国产av品久久久| 看十八女毛片水多多多| 久久ye,这里只有精品| 丝袜喷水一区| 一级,二级,三级黄色视频| 精品久久久精品久久久| 欧美老熟妇乱子伦牲交| 国产97色在线日韩免费| 国产成人欧美| 麻豆av在线久日| 在现免费观看毛片| 岛国毛片在线播放| 男女国产视频网站| 搡老岳熟女国产| 观看av在线不卡| 麻豆乱淫一区二区| 国产成人精品久久二区二区91| 99国产精品99久久久久| 宅男免费午夜| 色综合欧美亚洲国产小说| www.自偷自拍.com| 午夜av观看不卡| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人| 精品高清国产在线一区| 国产精品二区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美+亚洲+日韩+国产| 久久性视频一级片| 日韩 亚洲 欧美在线| 一级毛片黄色毛片免费观看视频| 亚洲精品久久午夜乱码| 国产av一区二区精品久久| 久久久久视频综合| 这个男人来自地球电影免费观看| 久久午夜综合久久蜜桃| 91成人精品电影| 99精品久久久久人妻精品| 国产精品国产三级国产专区5o| 精品免费久久久久久久清纯 | 黄网站色视频无遮挡免费观看| 大话2 男鬼变身卡| 精品卡一卡二卡四卡免费| 黄色一级大片看看| 人人澡人人妻人| 亚洲,欧美精品.| 国产亚洲午夜精品一区二区久久| 少妇精品久久久久久久| 日韩制服骚丝袜av| 中文字幕亚洲精品专区| 国产色视频综合| 国产亚洲av高清不卡| 男的添女的下面高潮视频| 日韩免费高清中文字幕av| 久久女婷五月综合色啪小说| 看免费av毛片| 91精品三级在线观看| 国产精品 欧美亚洲| 国产xxxxx性猛交| 七月丁香在线播放| 男人舔女人的私密视频| 免费观看a级毛片全部| 超碰97精品在线观看| 久久精品亚洲av国产电影网| 亚洲午夜精品一区,二区,三区| 久久人妻熟女aⅴ| av国产久精品久网站免费入址| 国产精品一区二区在线观看99| 久久精品久久精品一区二区三区| 一本一本久久a久久精品综合妖精| 中文字幕最新亚洲高清| 午夜福利一区二区在线看| 亚洲欧洲日产国产| 青春草视频在线免费观看| 成人三级做爰电影| av在线老鸭窝| 亚洲国产av新网站| 91字幕亚洲| 国产极品粉嫩免费观看在线| 欧美日韩黄片免| 日本av免费视频播放| 久久久精品区二区三区| 久久精品亚洲av国产电影网| 久久久精品区二区三区| 久久九九热精品免费| 欧美亚洲日本最大视频资源| 大香蕉久久成人网| 国产视频一区二区在线看| 超色免费av| 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 亚洲精品一卡2卡三卡4卡5卡 | 日韩欧美一区视频在线观看| 一级片'在线观看视频| 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 熟女少妇亚洲综合色aaa.| 精品视频人人做人人爽| 欧美日韩亚洲国产一区二区在线观看 | 男女边吃奶边做爰视频| 精品一品国产午夜福利视频| 免费久久久久久久精品成人欧美视频| 国产真人三级小视频在线观看| 777米奇影视久久| 日韩 欧美 亚洲 中文字幕| av又黄又爽大尺度在线免费看| 日本欧美视频一区| 午夜激情久久久久久久| 男女无遮挡免费网站观看| 啦啦啦在线观看免费高清www| 精品一区二区三卡| 亚洲精品国产一区二区精华液| avwww免费| 欧美国产精品一级二级三级| 赤兔流量卡办理| 在线观看国产h片| 高清视频免费观看一区二区| 制服人妻中文乱码| 99国产精品免费福利视频| 久久精品久久精品一区二区三区| 国产欧美亚洲国产| 午夜日韩欧美国产| 久久国产精品大桥未久av| 婷婷色综合www| 交换朋友夫妻互换小说| 黄色视频在线播放观看不卡| 日韩av不卡免费在线播放| 一级a爱视频在线免费观看| 汤姆久久久久久久影院中文字幕| 啦啦啦啦在线视频资源| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 国产成人免费观看mmmm| 欧美在线黄色| 多毛熟女@视频| 国产成人一区二区在线| 免费在线观看日本一区| 国产一区二区三区综合在线观看| 黑人欧美特级aaaaaa片| 精品人妻熟女毛片av久久网站| 国精品久久久久久国模美| 老司机影院毛片| 五月天丁香电影| 亚洲精品自拍成人| 日韩中文字幕欧美一区二区 | 国产精品久久久久久精品古装| 99热国产这里只有精品6| 国产野战对白在线观看| 777米奇影视久久| 热99久久久久精品小说推荐| 久久ye,这里只有精品| 精品人妻熟女毛片av久久网站| 午夜福利视频精品| 国产伦人伦偷精品视频| 欧美日韩精品网址| 只有这里有精品99| 三上悠亚av全集在线观看| 国产麻豆69| 欧美精品高潮呻吟av久久| 亚洲欧洲日产国产| 亚洲,欧美,日韩| 丝袜美足系列| 男女免费视频国产| 97精品久久久久久久久久精品| 91国产中文字幕| 国产爽快片一区二区三区| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 50天的宝宝边吃奶边哭怎么回事| 亚洲人成77777在线视频| www日本在线高清视频| 汤姆久久久久久久影院中文字幕| 老司机靠b影院| 99热网站在线观看| h视频一区二区三区| 在线观看国产h片| 久久鲁丝午夜福利片| 丰满少妇做爰视频| 国产精品 国内视频| 婷婷成人精品国产| 操美女的视频在线观看| 青春草视频在线免费观看| 一边亲一边摸免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲av片天天在线观看| 人成视频在线观看免费观看| 侵犯人妻中文字幕一二三四区| 99香蕉大伊视频| 久久久国产欧美日韩av| av一本久久久久| 美女高潮到喷水免费观看|