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

    水肥耦合驅(qū)動下的番茄植株形態(tài)模擬模型

    2015-02-21 18:46:32常毅博李建明尚曉梅張大龍潘銅華杜清潔
    關(guān)鍵詞:葉數(shù)節(jié)位株高

    常毅博,李建明,尚曉梅,張大龍,潘銅華,杜清潔

    (西北農(nóng)林科技大學(xué) a 園藝學(xué)院,b 經(jīng)濟管理學(xué)院,陜西 楊凌 712100)

    水肥耦合驅(qū)動下的番茄植株形態(tài)模擬模型

    常毅博a,李建明a,尚曉梅b,張大龍a,潘銅華a,杜清潔a

    (西北農(nóng)林科技大學(xué) a 園藝學(xué)院,b 經(jīng)濟管理學(xué)院,陜西 楊凌 712100)

    【目的】 建立土壤水肥影響下的番茄株高、單株展開葉數(shù)及單葉面積隨輻熱積的變化模型,為番茄施肥和灌溉提供依據(jù)?!痉椒ā?以灌水上限和施肥量為因素,采用二元二次正交旋轉(zhuǎn)組合設(shè)計,建立基于輻熱積的番茄株高及單葉面積的logistic模型和單株展開葉數(shù)的指數(shù)函數(shù)模型,并根據(jù)建立的單株展開葉數(shù)模型和單葉面積模型,模擬番茄葉面積指數(shù)?!窘Y(jié)果】 所建模型對不同水肥處理下的番茄株高、單株展開葉數(shù)、單葉面積和葉面積指數(shù)的預(yù)測效果較好,其預(yù)測值與實測值之間基于1∶1直線的決定系數(shù)R2分別為0.962 7,0.947 1,0.854 8和 0.926 3,相對誤差(RE)分別為9.41%,7.50%,17.80%和18.20%。【結(jié)論】 所建模型對番茄株高、單株展開葉數(shù)、單葉面積、葉面積指數(shù)4個指標(biāo)的預(yù)測精度均達80%以上,能較好地預(yù)測灌水上限和施肥量對番茄形態(tài)的動態(tài)影響。

    水肥耦合;株高;單株展開葉數(shù);單葉面積;葉面積指數(shù);輻熱積

    作物生長發(fā)育模型是將環(huán)境因子與作物生長發(fā)育相結(jié)合,實現(xiàn)作物生產(chǎn)管理經(jīng)驗數(shù)字化,從而指導(dǎo)溫室作物高質(zhì)高效生產(chǎn)的必要手段[1]。其中溫度、光輻射是影響作物生長發(fā)育的重要因素[2-3]。李建明等[4]研究了溫度、光輻射及水分對甜瓜幼苗生長的影響,建立了光輻射、日溫差和有效積溫驅(qū)動下的甜瓜幼苗干物質(zhì)積累與分配模型;李永秀等[5-6]用“輻熱積”法模擬黃瓜葉面積,構(gòu)建了光輻射和熱效應(yīng)影響下的不同整枝方式葉面積模擬模型,模擬精度較高;李青林等[7]基于“輻熱積”法模擬了黃瓜葉片形態(tài)特性、葉柄長度和直徑以及節(jié)間高度,建立了黃瓜的可視化生長模型。根系水肥環(huán)境對于作物的生長發(fā)育也有顯著影響[8]。李靈芝等[9]研究了不同氮濃度對溫室番茄生長發(fā)育的影響,結(jié)果表明,番茄的葉面積指數(shù)、株高、葉片數(shù)隨著營養(yǎng)液中氮濃度的增加而增大,但當(dāng)?shù)獫舛瘸^一定限度時,上述指標(biāo)會受到抑制;韋澤秀等[10]研究了水肥組合對日光溫室黃瓜葉片生長的影響,結(jié)果表明,隨土壤含水率和肥料施用量的增加,黃瓜葉面積和葉片擴展速率均會增加。然而,前人的研究主要集中在輻熱積對作物生長發(fā)育的動態(tài)變化的影響[4-7]上,或是研究不同水分、施肥單一因素處理對于作物形態(tài)的靜態(tài)影響[9-10],而未能綜合考慮水肥互作和輻熱積影響下作物生長發(fā)育的動態(tài)變化規(guī)律。為此,本試驗建立了水肥耦合下的番茄植株株高、單株展開葉數(shù)和單葉面積隨輻熱積的變化模型,實現(xiàn)了番茄植株形態(tài)隨光、溫和水肥環(huán)境變化的動態(tài)模擬,以期為溫室環(huán)境控制、番茄根系施肥和灌溉提供理論依據(jù)。

    1 材料與方法

    1.1 試驗設(shè)計

    試驗1于2013-03-07在西北農(nóng)林科技大學(xué)北校區(qū)園藝場3號塑料大棚內(nèi)進行。供試土壤理化性質(zhì)為:土壤體積質(zhì)量1.38 g/cm3,最大田間持水量24.3%,有機質(zhì)13.97 g/kg,堿解氮 85.76 mg/kg,速效磷19.42 mg/kg,速效鉀 100.17 mg/kg,pH值7.5。番茄植株于2013-03-09定植,07-15拉秧,行距為35 cm,株距為80 cm。試驗以土壤水分和施肥量為因素,采用二元(1/2 實施)二次正交旋轉(zhuǎn)組合設(shè)計,各處理因素及水平見表1。灌溉下限、灌溉上限均用最大田間持水量的百分比表示,各處理灌水下限均為田間最大持水量的50%,當(dāng)土壤水分含量降低至灌溉下限時,灌水將其補充到試驗設(shè)計的灌溉上限。根據(jù)顏冬云等[11]的研究結(jié)果,番茄栽培中m(N)∶m(P2O5)∶m(K2O)為1∶0.65∶0.93時肥效最好,且肥料利用率高,而本試驗番茄在生長期內(nèi)溫度光照充足,長勢快,因此施肥采用m(N)∶m(P2O5)∶m(K2O)=2∶1∶2的比例,提高了氮肥和鉀肥的比例,能夠保證番茄生長發(fā)育及果實產(chǎn)量品質(zhì)對氮鉀元素的需求。除磷肥全部作基肥一次性施入外,氮肥和鉀肥均分4次等量施入。為防止相鄰小區(qū)水分及養(yǎng)分互相影響,各小區(qū)之間埋50 cm深的薄膜隔離,并在兩頭設(shè)保護行。試驗1的數(shù)據(jù)用于建立模型。

    注:x1為灌溉上限對應(yīng)編碼,x2為施肥量對應(yīng)編碼。灌溉上限為基于田間持水量的百分比。

    Note:x1stands for the code of irrigation upper limit,x2stands for the code of fertilizer. Irrigation upper limit is the percentage of field moisture capacity.

    試驗2于2013-03-06在西北農(nóng)林科技大學(xué)北校區(qū)園藝場日光溫室內(nèi)進行。溫室東西長51 m,南北跨度8.3 m,脊高3.6 m。供試土壤體積質(zhì)量1.37 g/cm3,有機質(zhì)9.83 g/kg,堿解氮 107.36 mg/kg,速效磷28.15 mg/kg,速效鉀 99.84 mg/kg,土壤pH 7.3~7.5。試驗設(shè)置3個灌溉上限水平(最大田間持水量的70%,80%,90%)、3個施肥量水平(N、P2O5、K2O用量為200,100,200和375,187,375;550,275,550 kg/hm2),進行完全組合設(shè)計,共9個處理。試驗2栽培管理與試驗1相同。試驗2數(shù)據(jù)用于檢驗?zāi)P汀?/p>

    1.2 測定指標(biāo)及方法

    1.2.1 株 高 用卷尺測量從番茄莖基部至生長點的距離,記為株高(mm)。

    1.2.2 葉面積 2013-03-27起,從各處理中選取 4~6 株植株,每隔7~10 d剪取1次番茄葉片,采用剪紙法[6]測量供試葉片的面積,并測量剪取葉片的葉長(l)和葉寬(w),與所測實際葉面積回歸得到葉片面積(LA)模型:LA=0.495(l×w)0.911。

    (1)

    1.2.3 土壤相對質(zhì)量含水量 利用TDR水分測量儀,測量土壤水分相對體積含水量并轉(zhuǎn)化成相對質(zhì)量含水量。每次測量在小區(qū)內(nèi)隨機選取5個測試點,其平均值即為當(dāng)時的土壤相對質(zhì)量含水量。每次灌水量用公式(2)計算,在營養(yǎng)生長和開花坐果前期,計劃濕潤層深度為20 cm。

    V=rs×S×h×Q×(q1-q2) ×η。

    (2)

    式中:V為灌水量,L;rs為土壤體積質(zhì)量,g/cm3;S為小區(qū)面積,m2;h為計劃濕潤層深度,cm;Q為田間持水量;q1為土壤水分上限;q2為土壤水分下限;η為灌溉效率,取100%。

    1.2.4 環(huán)境監(jiān)測 采用物格環(huán)境記錄儀,每0.5 h記錄1次溫室溫度(℃)和光合有效輻射強度(W/m2)。

    2 番茄植株生長模擬模型的建立

    2.1 模型建立的基礎(chǔ)數(shù)據(jù)

    2.1.1 累積輻熱積 溫度和光輻射對于作物生長形態(tài)的影響可用輻熱積(Product of thermal effectiveness and PAR,簡稱TEP)來量化。倪紀恒等[12]將輻熱積定義為熱效應(yīng)和光合有效輻射的乘積,具體計算公式如下:

    RTE(T)=

    (3)

    PAR=0.5×Q′。

    (4)

    (5)

    TEP=∑(RTEP)。

    (6)

    式中:RTE(T)表示溫度為T時的相對熱效應(yīng);Tb為生長下限溫度,℃;Tab為生長的最適溫度下限,℃;Tou為生長的最適溫度上限,℃;Tm為生長上限溫度,℃;PAR為1 h 內(nèi)的總光合有效輻射,J/(m2·h);Q′為該時段內(nèi)的太陽總輻射,J/(m2·h);RTEP為每日相對輻熱積,MJ/(m2·d);RTE(i)、PAR(i)分別為1 d內(nèi)第i小時的相對熱效應(yīng)和總光合有效輻射,J/(m2·h);106是將J/(m2·h)換算成MJ/(m2·h)的單位換算系數(shù);TEP為累積輻熱積,MJ/m2。

    番茄生長發(fā)育的最適溫度白天為20~25 ℃,夜間為15~18 ℃,生長上限溫度為40 ℃,下限為10 ℃。

    2.1.2 傳統(tǒng)logistic曲線 對于作物生長的模擬,常用到傳統(tǒng)的logistic模型。傳統(tǒng)logistic模型[13]的解析表達式為:

    (7)

    式中:y表示植株的生長形態(tài)(y>0);k、a、b為方程的參數(shù),其中k為一定時期內(nèi)的生長上限(k>0),a為內(nèi)稟生長率(a>0),b為與曲線位置有關(guān)的參數(shù);t表示作物定植后的時間,d。

    2.2 株高的模擬

    傳統(tǒng)的logistic模型是解釋建立在生長時間基礎(chǔ)上的作物生長過程,忽略了每天的環(huán)境差異對于作物生長的影響。因此本研究用累積輻熱積(TEP,MJ/m2)模擬株高的生長,同時考慮水肥因子對于作物株高變化的影響,即:

    t=TEP。

    (8)

    k=kh(w,f)。

    (9)

    b=bh(w,f)。

    (10)

    a=ah(w,f)。

    (11)

    式中:kh(w,f)為基于水肥影響下的番茄株高上限;bh(w,f) 為基于水肥影響下的番茄株高位置參數(shù);ah(w,f)為基于水肥影響下的番茄內(nèi)稟生長率。

    由(7)~(11)式可得,基于水肥影響下的番茄株高模型為:

    (12)

    用不同水肥處理下實測的番茄營養(yǎng)生長期株高,擬合公式(12)得到結(jié)果如表2所示。以灌溉上限編碼值(x1)、施肥量編碼值(x2)為自變量,logistic函數(shù)的參數(shù)kh(w,f)、bh(w,f)、ah(w,f)為目標(biāo)函數(shù),利用DPS 7.5數(shù)據(jù)處理軟件回歸,配置水肥耦合效應(yīng)函數(shù),在α=0.10顯著水平下剔除不顯著項后的回歸方程為:

    (13)

    (14)

    式中:W為灌溉上限,%;M為土壤施肥量(m(N)∶m(P2O5)∶m(K2O)=2M∶M∶2M),kg/hm2。

    對(13)式進行F檢驗,經(jīng)計算得F1=0.28F0.01(6,5)=4.95。說明水肥耦合處理與番茄株高生長的logistic曲線參數(shù)回歸關(guān)系達到了顯著水平,因此用此水肥耦合回歸模型模擬番茄株高的變化比較可靠。

    注:R2為決定系數(shù),表4同。

    Note:R2stands for the coefficient of determination.The same for table 4.

    2.3 葉面積的模擬

    2.3.1 番茄單株展開葉數(shù)的模擬 試驗測得的番茄營養(yǎng)生長期單株展開葉數(shù)與累積輻熱積如圖1所示。

    圖1顯示,番茄單株展開葉數(shù)隨累積輻熱積的升高呈指數(shù)曲線增加。因此單株展開葉數(shù)可用下式表示:

    N=exp[c(w,f)×TEP]+N0-1。

    (15)

    式中:N為單株展開葉數(shù);c(w,f)為基于水肥影響的單株展開葉數(shù)模擬模型的參數(shù);N0為初始葉片數(shù),本試驗取為7。

    以灌溉上限編碼值(x1)、施肥量編碼值(x2)為自變量,c(w,f)為目標(biāo)函數(shù),利用DPS 7.5數(shù)據(jù)處理軟件回歸,配置水肥耦合效應(yīng)函數(shù),剔除α=0.10顯著水平上的不顯著項后,可得回歸方程為:

    (16)

    對回歸模型(16)進行F檢驗,經(jīng)計算得F1=0.95F0.05(6,5)=4.95,擬合較好。由公式(15)、(16)可以看出,在番茄生長之初,即TEP較小時,c(w,f)×TEP較小,即土壤水肥因子對番茄單株展開葉數(shù)的影響較??;隨著輻熱積的不斷積累,c(w,f)×TEP變大,則土壤水肥差異造成的番茄單株展開葉數(shù)差異也變大。

    2.3.2 番茄單葉面積的模擬 以T6處理為例,說明各個節(jié)位葉片面積隨累積輻熱積的變化規(guī)律,結(jié)果見圖2。從圖2可以看出,番茄單葉面積變化符合logistic函數(shù)規(guī)律,且葉片節(jié)位越高單葉面積增長速率越大。番茄各節(jié)位單葉面積隨累積輻熱積變化的logistic模型如下:

    (17)

    式中:i為葉片節(jié)位,指從植株底端到頂端的葉序,i=1時,表示植株底端第1片葉;ITEPi為第i片葉開始生長時的初始TEP;kla(i,w,f)、bla(i,w,f)、ala(i,w,f)為與葉片節(jié)位、灌溉上限和施肥量編碼值有關(guān)的logistic函數(shù)參數(shù)。

    對灌溉上限、施肥量和葉片節(jié)位與kla(i,w,f)、bla(i,w,f)、ala(i,w,f)進行相關(guān)性分析,結(jié)果見表3。表3顯示,葉片節(jié)位對kla(i,w,f)、bla(i,w,f)、ala(i,w,f)值的影響均達到極顯著水平(P<0.01)。因此,首先建立kla(i,w,f)、bla(i,w,f)、ala(i,w,f)與葉片節(jié)位之間的關(guān)系式。依據(jù)實測的番茄葉面積數(shù)據(jù),計算得kla(i,w,f)、bla(i,w,f)和ala(i,w,f) 的模擬值如圖3所示。

    注:“*”表示相關(guān)性顯著(P<0.05),“**”表示相關(guān)性極顯著(P<0.01)。

    Note:“*” Correlation is significant at theP=0.05 level,“**” Correlation is extremely significant at theP=0.01 level.

    圖3 水肥試驗中不同處理kla(i,w,f)、bla(i,w,f)及ala(i,w,f)值隨番茄葉片節(jié)位的變化
    Fig.3 Dynamic changes ofkla(n,w,f),bla(n,w,f),andala(n,w,f) values along with leaf order

    圖3顯示,kla(i,w,f)和bla(i,w,f)值隨葉片節(jié)位的增加呈logistic曲線增大,ala(i,w,f)隨葉片節(jié)位增加呈直線增大。建立各節(jié)位葉片葉面積參數(shù)kla(i,w,f)、bla(i,w,f)、ala(i,w,f)與葉片節(jié)位、灌水上限和施肥量編碼值的關(guān)系式為:

    (18)

    (19)

    (20)

    式中:kkla(w,f)、bkla(w,f)、akla(w,f)、kbla(w,f)、bbla(w,f)、abla(w,f)、pla(w,f)及qla(w,f)為與土壤灌溉上限和施肥量有關(guān)的參數(shù)。

    依據(jù)實測的不同水肥處理番茄各節(jié)位葉片面積,擬合公式(18)、(19)、(20)得到結(jié)果見表4。

    (21)

    對回歸模型(21)中各式進行F檢驗,經(jīng)計算得各式F1F0.05(6,5)=4.95,擬合較好。說明灌溉上限編碼值與施肥量編碼值能很好地模擬kkla(w,f)、kbla(w,f)、bkla(w,f)、bbla(w,f)、akla(w,f)、abla(w,f)、pla(w,f)和qla(w,f)。

    2.3.3 番茄葉面積指數(shù)的模擬 葉面積指數(shù)(LAI)是指單位土地面積上植物葉片總面積占土地面積的倍數(shù),可用下式計算:

    LAI=LAt×d/106,

    (22)

    (23)

    式中:LAt為番茄植株第1~n節(jié)位葉片總面積,mm2;LAi為第i個節(jié)位葉片面積,mm2;d為種植密度,株/m2;106為從mm2到m2的單位換算系數(shù)。

    3 番茄形態(tài)指標(biāo)模擬模型的驗證

    檢驗?zāi)P蜁r,常用回歸估計標(biāo)準(zhǔn)誤差(RMSE)和相對誤差(RE)對模擬值與實測值之間的符合度進行分析[14]。RMSE和RE可分別用式(24)和(25)計算:

    (24)

    (25)

    式中:OBSi為實測的番茄株高、單株展開葉數(shù)或單葉面積;SIMi為預(yù)測的番茄株高、單株展開葉數(shù)或單葉面積;n為樣本容量。RMSE值越小,表明模擬值與觀測值間的偏差越小,模型的預(yù)測精度越高。

    圖4結(jié)果表明,模型對于番茄株高、單株展開葉數(shù)、不同節(jié)位單葉面積和葉面積指數(shù)的預(yù)測值與實測值之間基于1∶1直線的決定系數(shù)(R2)分別為 0.962 7,0.947 1,0.854 8和0.926 3。由表5可以看出,本研究構(gòu)建的基于輻熱積、土壤水肥影響的番茄株高及單株展開葉數(shù)模型的模擬精度較高,其中株高模型標(biāo)準(zhǔn)誤差RMSE=5.22 mm,相對誤差RE=9.41%,小于10%;單株展開葉數(shù)模型標(biāo)準(zhǔn)誤差RMSE=0.92,相對誤差RE=7.5%,小于10%。而番茄單葉面積模型標(biāo)準(zhǔn)誤差RMSE=490.15 mm2,相對誤差RE=17.80%,大于15%;葉面積指數(shù)模型標(biāo)準(zhǔn)誤差RMSE=17 253.86 mm2,相對誤差RE=18.20%,大于15%。模型的預(yù)測精度均可以達到80%以上,能夠真實地反映水肥對番茄植株形態(tài)的影響。

    預(yù)測誤差除了來源于除模型自身因素外,可能還與其他未考慮在內(nèi)的因素有關(guān),如土壤理化性質(zhì)、番茄對土壤水肥的吸收利用等,因此需要通過進一步試驗來完善構(gòu)建的模型。

    4 結(jié)論與討論

    作物的形態(tài)生長受到其自身和外界環(huán)境的綜合影響,不同的作物種類、品種、栽培管理方法和外界環(huán)境均會造成作物生長過程的差異。本研究整合了包括溫度和光輻射氣象因子以及灌水上限和施肥情況等因子的綜合影響,對于番茄株高、單株展開葉數(shù)以及葉面積模型進行了初步的模擬,提出了基于輻熱積和水肥耦合的番茄株高、單株展開葉數(shù)以及單葉面積和葉面積指數(shù)模型。本研究通過二元二次正交旋轉(zhuǎn)組合設(shè)計,結(jié)合作物形態(tài)生長的曲線,認為水肥耦合對于番茄株高及單葉面積的影響可以通過logistic方程的參數(shù)變化來體現(xiàn),本研究建立的模型對單株展開葉數(shù)模擬精度較高,相對誤差小。本研究構(gòu)建的模型預(yù)測精度均可以達到80%以上,不僅能夠真實地反映番茄水肥對于番茄植株形態(tài)的影響,而且能夠清晰地闡釋水肥互作效應(yīng)對于植株形態(tài)的影響。

    將土壤灌溉上限和施肥量影響融入到logistic曲線中,直觀地表現(xiàn)了水肥對于番茄形態(tài)變化過程的影響,而未能深層次描述番茄根系對水肥的吸收[15-16]、同化及消耗利用[17],并轉(zhuǎn)化成植株形態(tài)的過程[18]。今后將進一步研究水肥對于番茄植株生理過程影響在形態(tài)上的外在表現(xiàn)。另外,本試驗著重模擬了水肥對番茄株高、單株展開葉數(shù)以及單葉面積、葉面積指數(shù)的模擬,忽略了番茄莖、花和果實等[19]其他形態(tài)指標(biāo)的變化。要建立能夠指導(dǎo)番茄栽培和管理的形態(tài)模型,精確模擬溫室作物動態(tài)空間生長過程,應(yīng)從多方位多角度動態(tài)模擬番茄整體形態(tài)變化過程[20-21],同時要考慮品種、人工管理和其他環(huán)境因子對于番茄生長的影響。

    [1] Marcelis L F M,Heuvelink E,Goudriaan J.Modelling biomass production and yield of horticultural crops:A review [J].Sci Hortic,1998,74(1):83-111.

    [2] Karlsson M G,Heins R D,Gerberick J O,et al.Temperature driven leaf unfolding rate in Hibiscus rosa-sinensis [J].Sci Hortic,1991,45(3/4):323-331.

    [3] Kinet J M.Effect of light conditions on the development of the inflorescence in tomato [J].Sci Hortic,1977,6(1):15-26.

    [4] 李建明,鄒志榮. 溫度、光輻射及水分對溫室甜瓜幼苗干物質(zhì)積累與分配的影響及其模擬模型 [J].應(yīng)用生態(tài)學(xué)報,2007,18(12):2715-2721.

    Li J M,Zou Z R.Effects of air temperature, solar radiation and soil water on dry matter accumulation and allocation of greenhouse muskmelon seedlings and related simulation models [J].Chinese Journal of Applied Ecology,2007,18(12):2715- 2721.(in Chinese)

    [5] 李永秀,羅衛(wèi)紅,倪紀恒,等.用輻熱積法模擬溫室黃瓜葉面積、光合速率與干物質(zhì)產(chǎn)量 [J].農(nóng)業(yè)工程學(xué)報,2005,21(12):131-136.

    Li Y X,Luo W H,Ni J H,et al.Simulation of leaf area,photosynthetic rate and dry matter production in greenhouse cucumber based on product of thermal effectiveness and photosynthetically active radiation [J].Transactions of the CSAE,2005,21(12):131-136.(in Chinese)

    [6] 李永秀,羅衛(wèi)紅,倪紀恒,等.基于輻射和溫度熱效應(yīng)的溫室水果黃瓜葉面積模型 [J].植物生態(tài)學(xué)報,2006,30(5):861-867.

    Li Y X,Luo W H,Ni J H,et al.Simulation of greenhouse cucumber leaf area based on radiation and thermal effectiveness [J].Journal of Plant Ecology,2006,30(5):861-867.(in Chinese)

    [7] 李青林,毛罕平,李萍萍.黃瓜地上部分形態(tài)-光溫響應(yīng)模擬模型 [J].農(nóng)業(yè)工程學(xué)報,2011,27(9):122-127.

    Li Q L,Mao H P,Li P P.Simulation of cucumber organ above-ground with relation to light and temperature [J].Transaction of the CSAE,2011,27(9):122-127.(in Chinese)

    [8] Jones H G,Tardieu F.Modeling water relations of horticultural crops: A review [J].Sci Hortic,1998,74(1):21-45.

    [9] 李靈芝,郭 榮,李海平,等.不同氮濃度對溫室番茄生長發(fā)育和葉片光譜特性的影響 [J].植物營養(yǎng)與肥料報,2010,16(4):965-969.

    Li L Z,Guo R,Li H P,et al.Effects of nitrogen concentration in hydroponics on growth and development of tomato and spectral characteristics of leaf in greenhouse [J].Plant Nutrition and Fertilizer Science,2010,16(4):965-969.(in Chinese)

    [10] 韋澤秀,梁銀麗,周茂娟,等.水肥組合對日光溫室黃瓜葉片生長和產(chǎn)量的影響 [J].農(nóng)業(yè)工程學(xué)報,2010,26(3):69-74.

    Wei Z X,Liang Y L,Zhou M J,et al.Physiological characteristics of leaf growth and yield of cucumber under different watering and fertilizer coupling treatments in greenhouse [J].Transactions of the CSAE,2010,26(3):69-74.(in Chinese)

    [11] 顏冬云,張 民.控釋復(fù)合肥對番茄生長效應(yīng)的影響研究 [J].植物營養(yǎng)與肥料學(xué)報,2005,11(1):110-115.

    Yan D Y,Zhang M.Effects of controlled-release and common compound fertilizers potted tomato [J].Plant Nutr Fert Sci,2005,11(1):110-115.(in Chinese)

    [12] 倪紀恒,羅衛(wèi)紅,李永秀,等. 溫室番茄葉面積與干物質(zhì)生產(chǎn)的模擬 [J].中國農(nóng)業(yè)科學(xué),2005,38(8):1629-1635.

    Ni J H,Luo W H,Li Y X,et al.Simulation of leaf area and dry matter production in greenhouse tomato [J].Scientia Agricultura Sinica,2005,38(8):1629-1635.(in Chinese)

    [13] 趙 龍,王振鳳,郭忠升,等.黃土丘陵半干旱區(qū)檸條林株高生長過程新模型 [J].生態(tài)學(xué)報,2013,33(7):2093-2103.

    Zhao L,Wang Z F,Guo Z S,et al.A new plant height growth process model ofCaraganaforest in semi-arid loess hilly region [J].Acta Ecologica Sinica,2013,33(7):2093-2103.(in Chinese)

    [14] 倪紀恒,毛罕平.電導(dǎo)率對溫室黃瓜葉面積和干物質(zhì)生產(chǎn)影響的動態(tài)模擬 [J].農(nóng)業(yè)工程學(xué)報,2011,27(12):105-109.

    Ni J H,Mao H P.Dynamic simulation of leaf area and dry matter production of greenhouse cucumber under different electrical conductivity [J].Transactions of the CSAE,2011,27(12):105-109.(in Chinese)

    [15] 楊啟良,張富倉,劉小剛,等.溝灌方式和水氮對玉米產(chǎn)量與水分傳導(dǎo)的影響 [J].農(nóng)業(yè)工程學(xué)報,2011,27(1):15-21.

    Yang Q L,Zhang F C,Liu X G,et al.Effects of different furrow irrigation patterns,water and nitrogen supply levels on hydraulic conductivity and yield of maize [J].Transactions of the CSAE,2011,27(1):15-21.(in Chinese)

    [16] 李淑文,文宏達,周彥珍,等.不同氮效率小麥品種氮素吸收和物質(zhì)生產(chǎn)特性 [J].中國農(nóng)業(yè)科學(xué),2006,39(10):1992-2000.

    Li S W,Wen H D,Zhou Y Z,et al.Characterization of nitrogen uptake and dry matter production in wheat varieties with different N efficiency [J].Scientia Agricultura Sinica,2006,39(10):1992-2000.(in Chinese)[17] 劉柯楠,楊啟良,郭 浩,等.基于不同滴灌方式的蘋果幼樹根系生長的三維可視化模擬 [J].中國生態(tài)農(nóng)業(yè)學(xué)報,2013,21(6):744-751.

    Liu K N,Yang Q L,Guo H,et al.Modeling and visualization of three-dimensional young apple tree root growth under different modes of drip irrigation [J].Chinese Journal of Eco-Agriculture,2013,21(6):744-751.(in Chinese)

    [18] 祁有玲,張富倉,李開峰.水分虧缺和施氮對冬小麥生長及氮素吸收的影響 [J].應(yīng)用生態(tài)學(xué)報,2009,20(10):2399-2405.

    Qi Y L,Zhang F C,Li K F.Effects of water deficit and nitrogen fertilization on winter wheat growth and nitrogen uptake [J].Chinese Journal of Applied Ecology,2009,20(10):2399-2405.(in Chinese)

    [19] 羅新蘭,李天來,仇家奇,等. 北方日光溫室長季節(jié)番茄莖節(jié)生長模擬模型 [J].農(nóng)業(yè)工程學(xué)報,2009,25(8):174-179.

    Luo X L,Li T L,Qiu J Q,et al.Simulation model of stem node growth of tomato in long season cultivation in northern solar greenhouse [J].Transactions of the CSAE,2009,25(8):174-179.(in Chinese)

    [20] 董喬雪,王一鳴,楊麗麗,等.番茄三維形態(tài)結(jié)構(gòu)的參數(shù)提取及模擬 [J].農(nóng)業(yè)工程學(xué)報,2010,26(S2):38-42.

    Dong Q X,Wang Y M,Yang L L,et al.Parameter identification of tomato 3D architectural model and simulation [J].Transactions of the CSAE,2010,26(S2):38-42.(in Chinese)

    [21] 王蕓蕓,溫維亮,郭新宇,等.煙草地上部植株三維重構(gòu)與可視化 [J].中國農(nóng)業(yè)科學(xué),2013,46(1):37-44.

    Wang Y Y,Wen W L,Guo X Y,et al.Research on three-dimensional reconstruction and visualization of above-ground tobacco plant [J].Scientia Agricultural Sinica,2013,46(1):37-44.(in Chinese)

    Simulation of tomato morphology growth with water-fertilizer coupling

    CHANG Yi-boa,LI Jian-minga,SHANG Xiao-meib,ZHANG Da-longa, PAN Tong-huaa,DU Qing-jiea

    (aCollegeofHorticulture,bCollegeofEconomicManagement,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

    【Objective】 This study established models for tomato stem length,number of unfolding leaves and single leaf area on product of thermal effectiveness and PAR (TEP) with water-fertilizer coupling to improve fertilization and irrigation of tomato roots.【Method】 Taking the upper limit of irrigation and fertilizer amount as factors,two quadratic general rotational combination design was used to build logistic regression model to assess the relationship between stem length,leaf area andTEPand exponential model between number of unfolding leaves andTEP.A leaf area index model was also established based on number of unfolding leaves and single leaf area.【Result】 Prediction accuracies (mean relative error) of tomato stem length model,number of unfolding leaves model,single leaf area model and leaf area index model were 9.41%,7.50%,17.80% and 18.20%,respectively.Their determination coefficients (R2) between simulated and measured values were 0.962 7,0.947 1,0.854 8 and 0.926 3,respectively.【Conclusion】 The prediction accuracies of the established models on the 4 indexes were all larger than 80%,indicating that the model can well predict the variations of tomato morphology characteristics.

    water and fertilizer coupling;stem length;number of unfolding leaves;single leaf area;leaf area index;thermal effectiveness and PAR

    2013-10-08

    國家“863”計劃項目(2011AA100504);國家“十二五”科技支撐計劃項目(2011BAD29B01)

    常毅博(1990-),男,河南靈寶人,碩士,主要從事設(shè)施作物生理生態(tài)研究。 E-mail:cyb2020@yeah.net

    李建明(1966-),男,陜西洛川人,教授,博士,博士生導(dǎo)師,主要從事設(shè)施園藝研究。E-mail:lijianming66@163.com

    時間:2015-01-05 08:59

    10.13207/j.cnki.jnwafu.2015.02.015

    Q945.12;S641.2

    A

    1671-9387(2015)02-0126-08

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20150105.0859.015.html

    猜你喜歡
    葉數(shù)節(jié)位株高
    皖南煙區(qū)烤煙新品種中川208的適宜留葉數(shù)研究
    玉米保護性耕作技術(shù)在遼陽地區(qū)的應(yīng)用效果研究
    ‘陽光玫瑰’葡萄不同節(jié)位冬芽激素含量差異分析
    煙草株高和葉數(shù)性狀QTL定位及候選基因預(yù)測
    摘葉對不同株型木薯品種產(chǎn)量和氮磷鉀素積累分配的影響
    陸地棉各節(jié)位種子油分含量的遺傳分析
    留果節(jié)位對哈密瓜果實品質(zhì)及產(chǎn)量的影響
    蔬菜(2020年1期)2020-01-17 03:27:50
    介紹四個優(yōu)良小麥品種
    不同栽培密度對柴胡生長的影響
    玉米骨干親本及其衍生系中基因的序列變異及與株高等性狀的關(guān)聯(lián)分析
    男人爽女人下面视频在线观看| 美女中出高潮动态图| 99久久精品国产亚洲精品| 午夜影院在线不卡| 我要看黄色一级片免费的| 成年女人毛片免费观看观看9 | 久久鲁丝午夜福利片| 亚洲av综合色区一区| 热re99久久国产66热| 99久久精品国产亚洲精品| 国产精品 国内视频| 国产探花极品一区二区| 18禁动态无遮挡网站| 91精品国产国语对白视频| 亚洲av成人精品一二三区| 精品免费久久久久久久清纯 | 免费观看a级毛片全部| 亚洲男人天堂网一区| svipshipincom国产片| 久久国产精品大桥未久av| 亚洲精品中文字幕在线视频| av国产精品久久久久影院| 成人国产麻豆网| 你懂的网址亚洲精品在线观看| 日韩电影二区| 18禁观看日本| 一本大道久久a久久精品| 激情五月婷婷亚洲| 中文字幕色久视频| 香蕉国产在线看| 永久免费av网站大全| 黄频高清免费视频| 两个人免费观看高清视频| 亚洲情色 制服丝袜| 国产精品秋霞免费鲁丝片| 中文欧美无线码| 亚洲av男天堂| 少妇人妻 视频| 精品国产一区二区久久| 一级,二级,三级黄色视频| 看十八女毛片水多多多| 国产在线一区二区三区精| 女性被躁到高潮视频| 国产一区亚洲一区在线观看| 不卡视频在线观看欧美| 午夜免费男女啪啪视频观看| 日韩欧美一区视频在线观看| 午夜福利一区二区在线看| 欧美av亚洲av综合av国产av | 亚洲国产日韩一区二区| a 毛片基地| 91精品国产国语对白视频| 国产亚洲午夜精品一区二区久久| 国产极品粉嫩免费观看在线| 一级片免费观看大全| 欧美在线黄色| 精品国产超薄肉色丝袜足j| 啦啦啦啦在线视频资源| 欧美日韩亚洲高清精品| 18禁动态无遮挡网站| 伊人久久国产一区二区| 自线自在国产av| 99热全是精品| 亚洲av综合色区一区| 国产片特级美女逼逼视频| 亚洲精品av麻豆狂野| 亚洲国产欧美一区二区综合| 熟女av电影| 韩国精品一区二区三区| 一边摸一边抽搐一进一出视频| 校园人妻丝袜中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 在线亚洲精品国产二区图片欧美| 自拍欧美九色日韩亚洲蝌蚪91| av.在线天堂| 午夜福利,免费看| 国产淫语在线视频| 欧美精品高潮呻吟av久久| 天天躁夜夜躁狠狠躁躁| 亚洲精品国产区一区二| 国产 精品1| 久久精品国产综合久久久| 男女高潮啪啪啪动态图| 最近最新中文字幕免费大全7| 亚洲美女黄色视频免费看| 亚洲精品视频女| 黄频高清免费视频| 午夜福利免费观看在线| 人妻一区二区av| 免费看av在线观看网站| 又黄又粗又硬又大视频| 9191精品国产免费久久| 黑人欧美特级aaaaaa片| 成人18禁高潮啪啪吃奶动态图| 欧美另类一区| e午夜精品久久久久久久| 美女主播在线视频| 伊人亚洲综合成人网| 成人免费观看视频高清| 老熟女久久久| 国产一级毛片在线| 久久久精品区二区三区| 亚洲,一卡二卡三卡| 国产精品亚洲av一区麻豆 | 亚洲一级一片aⅴ在线观看| 99精品久久久久人妻精品| 成人毛片60女人毛片免费| 亚洲一码二码三码区别大吗| 如日韩欧美国产精品一区二区三区| 免费高清在线观看视频在线观看| 婷婷色麻豆天堂久久| 中文字幕制服av| 在线免费观看不下载黄p国产| 一本久久精品| 国产av国产精品国产| 男女国产视频网站| 日韩一卡2卡3卡4卡2021年| 久久99精品国语久久久| 欧美精品一区二区大全| 色94色欧美一区二区| videosex国产| 国产毛片在线视频| 精品一品国产午夜福利视频| 国产精品久久久人人做人人爽| 1024香蕉在线观看| a级片在线免费高清观看视频| 人妻一区二区av| 久久久久久久久久久免费av| 国产成人精品久久二区二区91 | 伦理电影大哥的女人| 91精品三级在线观看| 久久av网站| 国产淫语在线视频| 无遮挡黄片免费观看| 国产日韩欧美亚洲二区| 欧美成人精品欧美一级黄| 亚洲精品在线美女| 国产一级毛片在线| 亚洲伊人色综图| 精品人妻一区二区三区麻豆| 午夜影院在线不卡| 亚洲精华国产精华液的使用体验| 精品视频人人做人人爽| 80岁老熟妇乱子伦牲交| 高清视频免费观看一区二区| 亚洲,欧美,日韩| 啦啦啦 在线观看视频| 久久女婷五月综合色啪小说| 精品国产一区二区三区久久久樱花| 99国产综合亚洲精品| 你懂的网址亚洲精品在线观看| 亚洲精品日韩在线中文字幕| 精品第一国产精品| a级毛片黄视频| 男女无遮挡免费网站观看| 亚洲成av片中文字幕在线观看| 午夜福利网站1000一区二区三区| 在线观看免费日韩欧美大片| 精品久久蜜臀av无| 91精品国产国语对白视频| 青春草国产在线视频| 女人高潮潮喷娇喘18禁视频| 乱人伦中国视频| 精品国产一区二区三区四区第35| 亚洲国产成人一精品久久久| 一区二区三区激情视频| 超碰97精品在线观看| 天天影视国产精品| 国产av精品麻豆| 一级片免费观看大全| 男女边摸边吃奶| 一级片免费观看大全| 日韩人妻精品一区2区三区| 久久久亚洲精品成人影院| 国产男女超爽视频在线观看| 午夜福利免费观看在线| 日韩熟女老妇一区二区性免费视频| 美女高潮到喷水免费观看| 精品酒店卫生间| 欧美日韩亚洲综合一区二区三区_| 久久久国产精品麻豆| 久久人妻熟女aⅴ| 成人手机av| 最近最新中文字幕免费大全7| 2018国产大陆天天弄谢| 国产亚洲午夜精品一区二区久久| 亚洲精品久久午夜乱码| 下体分泌物呈黄色| 欧美日韩亚洲高清精品| 国产野战对白在线观看| 精品视频人人做人人爽| 亚洲第一青青草原| 人体艺术视频欧美日本| 精品福利永久在线观看| 又大又黄又爽视频免费| 各种免费的搞黄视频| 欧美精品av麻豆av| 久久精品久久精品一区二区三区| 校园人妻丝袜中文字幕| 啦啦啦视频在线资源免费观看| 青青草视频在线视频观看| 亚洲综合色网址| 久久久久视频综合| 亚洲av福利一区| 国产亚洲欧美精品永久| 国产亚洲午夜精品一区二区久久| 搡老乐熟女国产| 91老司机精品| 亚洲 欧美一区二区三区| 日日摸夜夜添夜夜爱| 亚洲国产欧美网| 2021少妇久久久久久久久久久| 男女午夜视频在线观看| av在线老鸭窝| 黄色一级大片看看| 美女大奶头黄色视频| 久久这里只有精品19| 国产亚洲午夜精品一区二区久久| 久久国产亚洲av麻豆专区| 欧美日韩av久久| 91成人精品电影| 欧美变态另类bdsm刘玥| 欧美日本中文国产一区发布| 欧美日韩亚洲综合一区二区三区_| 91国产中文字幕| 国产一区亚洲一区在线观看| 久久影院123| av女优亚洲男人天堂| 国产野战对白在线观看| 亚洲欧美一区二区三区黑人| 别揉我奶头~嗯~啊~动态视频 | 久久 成人 亚洲| 亚洲综合精品二区| 老司机亚洲免费影院| 国产成人欧美| 99久久人妻综合| 国产淫语在线视频| 亚洲精华国产精华液的使用体验| 久久鲁丝午夜福利片| 中文字幕人妻丝袜制服| 久久久久久久久久久久大奶| 国产精品无大码| 一区二区av电影网| 午夜激情久久久久久久| 久久人人爽av亚洲精品天堂| 国产1区2区3区精品| av女优亚洲男人天堂| 老司机靠b影院| 亚洲成人国产一区在线观看 | 国产一区二区激情短视频 | 亚洲自偷自拍图片 自拍| 亚洲av成人不卡在线观看播放网 | 久久性视频一级片| √禁漫天堂资源中文www| netflix在线观看网站| 丰满饥渴人妻一区二区三| 蜜桃在线观看..| 国产成人精品久久二区二区91 | 黄片无遮挡物在线观看| 香蕉国产在线看| 十分钟在线观看高清视频www| 日韩av免费高清视频| 观看美女的网站| 一本久久精品| 国产精品久久久久久精品古装| 日日啪夜夜爽| 亚洲人成网站在线观看播放| 国产毛片在线视频| 免费观看a级毛片全部| www.精华液| 国产乱来视频区| av线在线观看网站| 80岁老熟妇乱子伦牲交| 国产福利在线免费观看视频| 青草久久国产| 国产精品三级大全| 婷婷色综合www| 自线自在国产av| 成人亚洲精品一区在线观看| 色精品久久人妻99蜜桃| 一本色道久久久久久精品综合| 中文字幕av电影在线播放| 亚洲精品国产区一区二| 久久亚洲国产成人精品v| 性高湖久久久久久久久免费观看| 日韩大码丰满熟妇| 亚洲综合色网址| 亚洲国产最新在线播放| 在线亚洲精品国产二区图片欧美| 人妻 亚洲 视频| 一级毛片黄色毛片免费观看视频| 日本色播在线视频| 欧美久久黑人一区二区| 亚洲四区av| 国产乱人偷精品视频| 久久鲁丝午夜福利片| 午夜日本视频在线| 熟女少妇亚洲综合色aaa.| 精品国产国语对白av| 丝袜在线中文字幕| 亚洲视频免费观看视频| 制服丝袜香蕉在线| 极品人妻少妇av视频| 亚洲婷婷狠狠爱综合网| 欧美日韩亚洲高清精品| 在线天堂中文资源库| 免费人妻精品一区二区三区视频| 亚洲av福利一区| 美女福利国产在线| 老司机靠b影院| 大香蕉久久网| 欧美精品一区二区免费开放| 亚洲精品久久午夜乱码| 999精品在线视频| 啦啦啦视频在线资源免费观看| 国产人伦9x9x在线观看| 国产片特级美女逼逼视频| 新久久久久国产一级毛片| 亚洲国产日韩一区二区| 亚洲欧美成人综合另类久久久| 校园人妻丝袜中文字幕| 亚洲精品第二区| 欧美黑人精品巨大| 国产精品av久久久久免费| 午夜av观看不卡| 一级毛片电影观看| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品在线美女| 国产又爽黄色视频| 天堂中文最新版在线下载| 久久精品久久久久久噜噜老黄| 啦啦啦中文免费视频观看日本| 亚洲av成人不卡在线观看播放网 | 色婷婷av一区二区三区视频| 亚洲成av片中文字幕在线观看| 中文乱码字字幕精品一区二区三区| 色婷婷久久久亚洲欧美| 国产日韩欧美亚洲二区| 国产成人午夜福利电影在线观看| a级片在线免费高清观看视频| 亚洲av日韩精品久久久久久密 | 王馨瑶露胸无遮挡在线观看| 天堂俺去俺来也www色官网| 免费少妇av软件| 亚洲成人av在线免费| 亚洲一卡2卡3卡4卡5卡精品中文| 国产片内射在线| 亚洲一区中文字幕在线| 伊人久久国产一区二区| 国产高清不卡午夜福利| 少妇人妻 视频| 国产片特级美女逼逼视频| 宅男免费午夜| 性少妇av在线| 精品久久久精品久久久| 伊人久久国产一区二区| av.在线天堂| 大片电影免费在线观看免费| 少妇人妻精品综合一区二区| 国产又色又爽无遮挡免| 亚洲美女黄色视频免费看| 天美传媒精品一区二区| 18禁观看日本| 亚洲av日韩在线播放| 又粗又硬又长又爽又黄的视频| 亚洲一区中文字幕在线| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看 | 日韩制服丝袜自拍偷拍| 精品国产乱码久久久久久男人| 最新在线观看一区二区三区 | 99精品久久久久人妻精品| 校园人妻丝袜中文字幕| 亚洲av成人不卡在线观看播放网 | 女人爽到高潮嗷嗷叫在线视频| 精品国产一区二区三区四区第35| 狠狠婷婷综合久久久久久88av| 亚洲av成人精品一二三区| 国产成人免费无遮挡视频| 如何舔出高潮| 岛国毛片在线播放| 制服诱惑二区| 久久久久人妻精品一区果冻| 亚洲精品成人av观看孕妇| 另类亚洲欧美激情| 男女边摸边吃奶| 丰满乱子伦码专区| 在线免费观看不下载黄p国产| 国产av一区二区精品久久| 男人操女人黄网站| 亚洲精品久久午夜乱码| av线在线观看网站| 欧美日韩一区二区视频在线观看视频在线| 2018国产大陆天天弄谢| 妹子高潮喷水视频| 成年人免费黄色播放视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美清纯卡通| 欧美黑人精品巨大| 日本午夜av视频| 午夜福利影视在线免费观看| 国产日韩欧美在线精品| 中文天堂在线官网| 男女国产视频网站| 亚洲成人一二三区av| 亚洲成色77777| 国产精品国产三级国产专区5o| 亚洲在久久综合| 国产熟女午夜一区二区三区| 99久久综合免费| 成人午夜精彩视频在线观看| 国产黄色免费在线视频| 久久人人爽av亚洲精品天堂| 欧美 日韩 精品 国产| 久久精品国产亚洲av高清一级| 久久ye,这里只有精品| 欧美成人午夜精品| 亚洲欧美成人精品一区二区| 久久性视频一级片| av在线播放精品| 亚洲精品国产av蜜桃| 亚洲人成77777在线视频| 七月丁香在线播放| 午夜福利影视在线免费观看| 国产日韩欧美在线精品| 美女大奶头黄色视频| 中国国产av一级| 美女视频免费永久观看网站| 成年人免费黄色播放视频| 日韩伦理黄色片| 久久久久网色| 热99久久久久精品小说推荐| 少妇人妻 视频| 久久久久网色| 黄色一级大片看看| 一本色道久久久久久精品综合| 国产精品99久久99久久久不卡 | 亚洲国产av影院在线观看| 国产精品久久久久久精品电影小说| 捣出白浆h1v1| 国产男女内射视频| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 欧美少妇被猛烈插入视频| 久久久久久久大尺度免费视频| 天堂中文最新版在线下载| 久久影院123| 国产日韩欧美亚洲二区| 免费高清在线观看视频在线观看| av国产久精品久网站免费入址| 午夜免费观看性视频| 一区在线观看完整版| 好男人视频免费观看在线| 精品视频人人做人人爽| 国产 一区精品| 欧美 亚洲 国产 日韩一| 欧美人与善性xxx| 国产有黄有色有爽视频| 大片电影免费在线观看免费| 亚洲三区欧美一区| 亚洲av日韩精品久久久久久密 | 久久精品aⅴ一区二区三区四区| 视频区图区小说| 色婷婷久久久亚洲欧美| 日本黄色日本黄色录像| 欧美国产精品一级二级三级| 丝瓜视频免费看黄片| 亚洲精品视频女| 王馨瑶露胸无遮挡在线观看| a级片在线免费高清观看视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人午夜精品| 叶爱在线成人免费视频播放| 成人漫画全彩无遮挡| 在线观看免费视频网站a站| 久久久久视频综合| 天天躁夜夜躁狠狠躁躁| 青草久久国产| 欧美精品亚洲一区二区| 男女之事视频高清在线观看 | 狠狠婷婷综合久久久久久88av| 青春草亚洲视频在线观看| 国产亚洲av片在线观看秒播厂| 国产黄频视频在线观看| av电影中文网址| 高清视频免费观看一区二区| 各种免费的搞黄视频| www.自偷自拍.com| 国产极品天堂在线| 下体分泌物呈黄色| 精品国产露脸久久av麻豆| 亚洲国产精品一区三区| 少妇精品久久久久久久| 国产精品成人在线| 老司机深夜福利视频在线观看 | av又黄又爽大尺度在线免费看| 亚洲欧美中文字幕日韩二区| 亚洲精品美女久久久久99蜜臀 | 日韩一卡2卡3卡4卡2021年| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 欧美黑人欧美精品刺激| 免费久久久久久久精品成人欧美视频| 丝瓜视频免费看黄片| 男的添女的下面高潮视频| 王馨瑶露胸无遮挡在线观看| 国产在线视频一区二区| 天美传媒精品一区二区| 不卡视频在线观看欧美| 国产一区二区 视频在线| 丁香六月天网| av福利片在线| 国产精品一二三区在线看| 亚洲一码二码三码区别大吗| 啦啦啦中文免费视频观看日本| 久久久精品免费免费高清| 国产精品秋霞免费鲁丝片| 女的被弄到高潮叫床怎么办| 大香蕉久久成人网| 久久久久国产精品人妻一区二区| 亚洲国产精品一区二区三区在线| 嫩草影院入口| 九草在线视频观看| 日韩成人av中文字幕在线观看| 美女午夜性视频免费| 亚洲国产av新网站| 丝袜在线中文字幕| 一区二区三区精品91| 国产精品久久久久久精品电影小说| 久久国产精品大桥未久av| 黄色怎么调成土黄色| 80岁老熟妇乱子伦牲交| 午夜激情久久久久久久| 亚洲国产精品国产精品| 国产熟女午夜一区二区三区| 99热全是精品| 亚洲精品aⅴ在线观看| 人妻人人澡人人爽人人| 色网站视频免费| 人妻人人澡人人爽人人| 天天添夜夜摸| 这个男人来自地球电影免费观看 | 尾随美女入室| 成人漫画全彩无遮挡| 国产不卡av网站在线观看| 丰满迷人的少妇在线观看| 国产精品国产av在线观看| 日本猛色少妇xxxxx猛交久久| 久久国产精品男人的天堂亚洲| 国产成人免费观看mmmm| 国产乱人偷精品视频| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 国产极品粉嫩免费观看在线| 免费黄网站久久成人精品| 观看av在线不卡| 亚洲精品中文字幕在线视频| 久久免费观看电影| 在线免费观看不下载黄p国产| 老鸭窝网址在线观看| 午夜激情av网站| 天天躁夜夜躁狠狠久久av| 水蜜桃什么品种好| 久久久久久久久久久久大奶| a 毛片基地| 大话2 男鬼变身卡| www日本在线高清视频| 99精国产麻豆久久婷婷| 大片免费播放器 马上看| 大香蕉久久网| 黄色怎么调成土黄色| 十八禁网站网址无遮挡| 天美传媒精品一区二区| 岛国毛片在线播放| 精品少妇内射三级| 亚洲中文av在线| 国产成人免费无遮挡视频| 人人妻人人爽人人添夜夜欢视频| 国产精品三级大全| tube8黄色片| 精品一区二区三区四区五区乱码 | 亚洲国产av影院在线观看| 成年av动漫网址| 亚洲一卡2卡3卡4卡5卡精品中文| 热re99久久精品国产66热6| 大香蕉久久成人网| 亚洲成人一二三区av| 日本av免费视频播放| 超碰97精品在线观看| 国产精品成人在线| 午夜激情av网站| 最近最新中文字幕大全免费视频 | 日日摸夜夜添夜夜爱| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| 黄色一级大片看看| av在线播放精品| 国产欧美亚洲国产| 色播在线永久视频| 欧美激情极品国产一区二区三区| 日日爽夜夜爽网站| 午夜福利免费观看在线| 精品少妇内射三级| 黑丝袜美女国产一区| 人妻人人澡人人爽人人| 免费观看人在逋| 国产精品嫩草影院av在线观看| 高清黄色对白视频在线免费看| 曰老女人黄片| 日韩一本色道免费dvd| 青春草国产在线视频| 国产 一区精品| 久久久久久久国产电影| 性少妇av在线|