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

    基于作物模型與最佳季節(jié)法的錦州地區(qū)玉米最佳播種期分析*

    2016-04-06 12:01:03張玉書紀瑞鵬于文穎張淑杰
    中國農(nóng)業(yè)氣象 2016年1期

    米 娜,張玉書,紀瑞鵬,蔡 福,于文穎,張淑杰

    (中國氣象局沈陽大氣環(huán)境研究所,沈陽110166)

    ?

    基于作物模型與最佳季節(jié)法的錦州地區(qū)玉米最佳播種期分析*

    米 娜,張玉書**,紀瑞鵬,蔡 福,于文穎,張淑杰

    (中國氣象局沈陽大氣環(huán)境研究所,沈陽110166)

    摘要:為了探究錦州地區(qū)平均氣候狀態(tài)下玉米最佳播期,同時檢驗作物模型法和最佳季節(jié)法確定最佳播期的適用性,利用遼寧錦州農(nóng)田試驗站3a分期播種試驗數(shù)據(jù),在對作物生長模型CERES-Maize進行參數(shù)校正與模擬效果檢驗的基礎上,應用模型模擬不同播期下玉米30a(1981-2010年)的產(chǎn)量,同時應用最佳季節(jié)法分析該地區(qū)的玉米最佳播期,結合作物模型法的研究結果,提出對最佳季節(jié)法的改進辦法。結果表明:CERES-Maize模型能夠較好地模擬不同播期玉米的物候期和產(chǎn)量,其歸一化均方根誤差(NRMSE)小于10.3%,對不同播期下30a的產(chǎn)量模擬結果顯示,當播期從4月10日推遲至5月10日,玉米平均產(chǎn)量增加6%,當播期從5月10日推遲至5月30日,玉米產(chǎn)量中值從9112kg·hm-2降至8619kg·hm-2,4月25日和4 月30日播種玉米的平均產(chǎn)量與5月10日播種的玉米產(chǎn)量無顯著差異。結果顯示最佳季節(jié)法確定的錦州地區(qū)玉米最佳播期較為滯后,與作物模型法的研究結果及實際生產(chǎn)的播期有較大出入,因此,提出了對最佳季節(jié)法的改進辦法即將灌漿期間的不利氣溫條件考慮在內(nèi),改進后得到的最佳播期與作物模型法研究結果較一致。從30a平均氣候狀況看,該地區(qū)玉米的最佳播期在4月25日-5月10日,作物模型法有較好的適用性,最佳季節(jié)法經(jīng)過改進后也可實際應用。

    關鍵詞:CERES-Maize模型;最佳播期;最佳季節(jié)法

    米娜,張玉書,紀瑞鵬,等.基于作物模型與最佳季節(jié)法的錦州地區(qū)玉米最佳播種期分析[J].中國農(nóng)業(yè)氣象,2016,37(1):68-76

    玉米品質(zhì)及產(chǎn)量受生態(tài)環(huán)境(水、光照、溫度和土壤等)和栽培措施(密度、播種期、收獲期和肥料等)的影響較大[1],適宜的播期是實現(xiàn)作物高產(chǎn)的必要條件[2-3]。目前研究不同播期對作物產(chǎn)量影響的方法主要采用分期播種試驗法[4-6]、作物生長模型法[7-10]、最佳季節(jié)法[11]等。分期播種試驗具有直觀的優(yōu)點,但其結果只能反映一個地區(qū)試驗年份的狀況,不具有外推性。作物生長模型法通常在模型得到充分驗證之后用于模擬不同播期對作物產(chǎn)量的影響,因作物模型具有較強的機理性[12],因此,可以將短期的田間試驗結果外推至其它年份。最佳季節(jié)法是指依據(jù)玉米生長發(fā)育進程與最佳季節(jié)同步原理,使開花吐絲及灌漿成熟關鍵期處于最佳光溫水生態(tài)條件下的概率最大[11],據(jù)此確定玉米的最佳播期??傮w上,以往應用前兩種方法對作物最佳播期的研究較多,應用最佳季節(jié)法的研究報道并不多見,該方法是否具有較好的適用性,是否需要進一步完善,還有待深入。東北地區(qū)是中國春玉米的主產(chǎn)區(qū),玉米播種面積約占全區(qū)糧食作物面積的60%。東北地區(qū)玉米的主要生長季節(jié)為4-9月,期間限制玉米產(chǎn)量的主要農(nóng)業(yè)氣象災害有低溫冷害[13]、干旱[14]、霜凍[15]以及開花吐絲期的連陰雨等,因此,正確選擇播期是優(yōu)化作物產(chǎn)量的重要管理措施。本研究以遼寧錦州農(nóng)田試驗站為研究站點,利用3a分期播種試驗數(shù)據(jù)校正和檢驗作物生長模型CERES-Maize,利用錦州地區(qū)30a(1981-2010年)氣象數(shù)據(jù)驅(qū)動模型運行,模擬不同播期對玉米產(chǎn)量的影響,同時將應用最佳季節(jié)法確定的最佳播期與模型模擬結果進行比對,提出最佳季節(jié)法確定最佳播期方法的改進,旨在為農(nóng)業(yè)氣象業(yè)務中玉米播種期預報及最佳播期的確定提供方法參考與依據(jù)。

    1 材料與方法

    1.1 試驗設置及觀測

    2011-2013年在錦州農(nóng)田試驗站(41°08′N,121°07′E)設置玉米分期播種試驗,并觀測玉米生長狀況用以校正模型參數(shù)及檢驗模型。該研究區(qū)屬溫帶季風型大陸性氣候,四季分明,多年(1981-2010年)平均氣溫9.9℃,平均降水量568mm。2011和2012年各設置8個播期,播種日期分別為4月10、15、20、25、30日及5月10、20、30日,2013年設置7個播期,分別為4月16、20、25、30日及5 月10、20、30日。試驗品種為丹玉39,該品種為稀植大穗型晚熟品種,適宜在遼寧省南部及沿海地區(qū)種植。3a試驗中播種行距均為60cm,株距為40cm,播種時施底肥750kg·hm-2,底肥為復合肥料,N、P2O5、K2O的含量分別為12%、18%、15%,總養(yǎng)分≥45%。利用2011年數(shù)據(jù)校正模型參數(shù),2012年和2013年數(shù)據(jù)檢驗模型模擬效果。

    記錄所有播期的玉米發(fā)育期(出苗、三葉、七葉、拔節(jié)、抽雄、開花、乳熟、成熟),另外,針對4月10日(2013年為4月16日)、4月20日、4月30日和5月10日播種的處理,在三葉、七葉、拔節(jié)、抽雄和乳熟期各取3株代表性植株進行生物量測定,測定時,將植株齊地面割下,分離為莖、葉、果實部分,稱取鮮重,放入恒溫干燥箱內(nèi)殺青(100~105℃,1h),以后維持在70~80℃,烘干至恒重,稱取干重。乳熟期所測歷年(2011-2013年)植株密度分別為4.80、4.19和4.32株·m-2。玉米成熟后按照《農(nóng)業(yè)氣象觀測規(guī)范》測定籽粒產(chǎn)量,同時測定籽粒含水率,最終將籽粒產(chǎn)量校正至含水率0%,并與模型模擬產(chǎn)量進行比較。

    試驗期間的逐日氣象資料包括最高和最低氣溫、降雨時間和降雨量、日照時數(shù),均來自遼寧省氣象信息中心。

    1.2 作物模型校正和評價

    使用DSSAT(Decision Support System for Agrotechnology Transfer)中CERES-Maize(Version 4.5)[16-17]模擬玉米發(fā)育期和產(chǎn)量。利用2011年氣象數(shù)據(jù)和8個播期的玉米試驗數(shù)據(jù)進行參數(shù)調(diào)試,獲取該品種的遺傳參數(shù)(表1)。利用2012和2013年氣象資料驅(qū)動模型模擬15個播期的發(fā)育期、產(chǎn)量及生物量,將開花期、成熟日期及生物量、產(chǎn)量的模擬值與觀測值進行比對,利用誤差百分比、均方根誤差(RMSE)、歸一化均方根誤差(NRMSE)和決定系數(shù)(R2)等指標評價模型的模擬效果。其中歸一化均方根誤差可表征觀測值與模擬值的相對差異,NRMSE小于20%表示模擬效果好,20%~30%表示模擬效果較好,超過30%表示模擬效果較差[18]。

    式中,Pi為第i個模擬值,Oi為第i個觀測值,Oavg和Pavg分別為觀測平均值和模擬平均值,n為觀測值或模擬值的個數(shù)。

    表1 CERES-Maize模型中“丹玉39”遺傳參數(shù)取值Table 1 Genetic coefficients and its value for Danyu39 in the CERES-Maize model

    1.3 作物模型法

    使用錦州站30a(1981-2010年)逐日氣象觀測資料驅(qū)動作物模型CERES-Maize,模擬當?shù)?0a該品種不同播期的產(chǎn)量性狀,應用SPSS10軟件對不同播期的產(chǎn)量、單株籽粒數(shù)、百粒重進行統(tǒng)計,采用最小顯著性差異法(LSD)進行多重比較檢驗,分析各播期的產(chǎn)量特點。

    氣象資料包括日最高、最低氣溫、日照時數(shù)、降水時間和降水量,太陽輻射數(shù)據(jù)根據(jù)日照時數(shù)及每日的天文輻射量計算[19]。設置8個播期,分別為4月10、15、20、25、30日和5月10、20和30日。播種日施底肥,施肥量N、P、K分別為90、60、93kg·hm-2,施肥深度8cm,播種密度設為4.32株·m-2。為使玉米生長免受早霜影響,在模型中設置5月20日和5月30日兩個播期的收獲日期不晚于播種后134d和 124d(根據(jù)錦州地區(qū)最早初霜日在10月1日出現(xiàn)設置)。模擬產(chǎn)量結果用統(tǒng)計箱式圖表征,箱圖本身包括中間50%的數(shù)據(jù),箱圖上、下邊界分別代表75% 和25%的數(shù)據(jù),產(chǎn)量中值用箱圖中的一條水平線表示,線圖最上面和最下面的橫線代表90%和10%的產(chǎn)量數(shù)據(jù)。

    1.4 玉米生長發(fā)育過程與最佳季節(jié)同步法

    利用錦州站30a(1981-2010年)逐日平均氣象資料,采用玉米生長發(fā)育過程與最佳季節(jié)同步法[11]分析該品種生長發(fā)育過程與當?shù)貧夂蛸Y源的匹配情況,確定最適播種期。該方法確定玉米最佳播種期的原理是保證玉米生長發(fā)育進程與最佳季節(jié)同步,關鍵是使玉米的開花、結實期處于最佳的溫光水環(huán)境條件下,然后根據(jù)一定條件推算得到播種期,具體方法參考文獻[11]。

    首先,假定一個計算的起始日期,計算一定時間范圍內(nèi)光溫水組合情況。原則是使玉米在抽雄-吐絲期TED天內(nèi)、灌漿期GFD天內(nèi)處于最佳光溫水條件下的時間最長,即光溫的貢獻度(con1)與光溫水的貢獻度(con2)之和最大,則該日期為最佳開花日期。各指標的計算式為:

    式中,TED為玉米抽雄-吐絲期天數(shù),GFD為灌漿期天數(shù),SUMT_YPP為品種所需積溫,VM為品種早熟性,本研究取1,d1為TED內(nèi)日平均氣溫在25~28℃的天數(shù),d2為TED天內(nèi)日照時數(shù)介于6~10h的天數(shù),d3為GFD天內(nèi)日平均氣溫在22~24℃的天數(shù),d4為GFD天內(nèi)日照時數(shù)介于6~10h的天數(shù),d5為GFD天內(nèi)降雨量(mm),200為GFD天內(nèi)最適降雨量。

    具體做法為,利用遼寧錦州地區(qū)的常年(1981-2010年)氣象資料和玉米品種丹玉39的品種特征,從7月8日(開花日期)-8月5日每2d計算(Con1+Con2)的值,以確定玉米最佳開花日期。丹玉39的TED為7d,取開花日期前3d和后3d以及開花日當天,GFD為45d,起始日期為開花日期后10d。根據(jù)2011-2013年分期播種試驗得到的播種-開花日期所需積溫,從開花日期向前推算確定播種日期。

    2 結果與分析

    2.1 典型年作物模型模擬效果檢驗

    利用CERES-Maize模型及其參數(shù)調(diào)試結果對2012和2013年不同播期玉米生育期進行模擬,以播種后天數(shù)表示玉米開花和成熟日期,結果見圖1。由圖可見,實測不同播期玉米的開花日期和成熟日期分布在播種后67~90d和131~150d,而相應的模擬值范圍在63~87d和129~155d,觀測值與模擬值間相關系數(shù)分別為0.98和0.80,均通過0.01水平的顯著性檢驗。但與1:1線相比,15個播期開花期的模擬值普遍低于實測值,即模擬開花日期比實測提前,但總體誤差較小,NRMSE=3.4%,成熟期模擬值與實測值分布在1:1線兩測,總體誤差也不大,NRMSE=4.2%。

    CERES-Maize模型對2012和2013年指定播期(10個)玉米生物量的模擬如圖2所示。由圖可見,各播期玉米生物量累積過程均表現(xiàn)出“S”型變化特點,歸一化均方根誤差范圍在13.3%~31.2%,除2012年4月16日播種(圖2b)和2013年5月10日播種(圖2h)處理NRMSE較高外,其它處理均在10%~20%,模型對多數(shù)播期玉米生物量的模擬效果較好。

    CERES-Maize模型對2012和2013年共15個播期的產(chǎn)量模擬及檢驗結果如表2所示,2012和2013年產(chǎn)量模擬的誤差百分比分別為7.2%和9.4%,均方根誤差分別為752kg·hm-2和872kg·hm-2??傮w來看,該模型對玉米產(chǎn)量的模擬效果很好,2012年NRMSE 為8.1%;2013年為10.3%。

    圖1 2012和2013年玉米模擬發(fā)育期與實測發(fā)育期的比較Fig. 1 Comparison of simulated and observed anthesis date (a) and physiology maturity date (b) of maize in 2012 and 2013注:發(fā)育期以播種后的天數(shù)表示,實線為1:1線,虛線為擬合線Note: The growing date is expressed as days after sowing. Solid line and dotted line in this figure represent 1:1 and fitted line, respectively

    圖2 2012年和2013年不同播期玉米地上生物量的觀測值(點)與模擬值(線)的比較Fig. 2 Comparison of simulated (lines) and observed (dots) top biomass at different sowing dates in 2012 and 2013

    表2 2012年和2013年不同播期玉米產(chǎn)量模擬值與觀測值的比較Table 2 Comparison of observed and simulated average yield at different sowing dates in 2012 and 2013

    2.2 不同播期玉米產(chǎn)量模擬比較

    用統(tǒng)計箱式圖表征模擬產(chǎn)量結果,箱圖本身包括中間50%的數(shù)據(jù),箱圖上、下邊界分別代表75% 和25%的數(shù)據(jù),產(chǎn)量中值用箱圖中的一條水平線表示,線圖最上面和最下面的橫線代表90%和10%的產(chǎn)量數(shù)據(jù)。30a模擬結果顯示,從4月10日播種開始,隨著播期的延后玉米平均產(chǎn)量逐漸增加,直至播期為5月10日,玉米平均產(chǎn)量達到一個高值(8965kg·hm-2),產(chǎn)量中值也表現(xiàn)出明顯的躍升(從8754kg·hm-2升至9112kg·hm-2),之后隨著播期延后平均產(chǎn)量呈下降趨勢(圖3)。具體表現(xiàn)為,當播期從4月10日推遲至5月10日,玉米的開花日期從7月7日推遲到7月20日(表3),玉米平均產(chǎn)量增加6%;當播期從5月10日推遲至5月30日,開花日期從7月20日推遲至8月2日,玉米平均產(chǎn)量減少5%。播期為4月10日時,25%的年份玉米產(chǎn)量低于8179kg·hm-2,25%的年份玉米產(chǎn)量高于8960kg·hm-2;而播期為5月10日時,25%的年份玉米產(chǎn)量低于8529kg·hm-2,25%的年份玉米產(chǎn)量高于9257kg·hm-2。當播期從5月10日推遲至5 月30日,玉米產(chǎn)量中值從9112kg·hm-2降至8619kg·hm-2。因此,從過去30a不同播期對玉米產(chǎn)量影響的模擬結果來看,5月10日前后播種玉米更易獲得較高產(chǎn)量。

    圖3 不同播期玉米產(chǎn)量的模擬值Fig. 3 Simulated yield of maize for different sowing dates

    從不同播期玉米平均產(chǎn)量的變化來看(表4),以5月10日為播期的玉米產(chǎn)量與4月25日、4月30日和5月20日為播期的產(chǎn)量無顯著性差異。4月10、15、20日和5月30日為播期的玉米平均產(chǎn)量較5月10日為播期的產(chǎn)量約低5%(P<0.05)。各播期之間單株籽粒數(shù)差異不大,5月10日播期的單株籽粒數(shù)與其它播期均無顯著性差異,4月播種的玉米單株籽粒數(shù)較5月30日偏少9%左右,且差異均達到顯著水平(P<0.05)。百粒重則隨播期的變化呈拋物線型變化,5月20日以后播種的玉米其百粒重較之前播種的約偏輕6%,與4月25日、4月30日和5 月10日播種的差異均達到顯著水平(P<0.05)。

    表3 CERES-Maize模型對不同播期玉米物候期30a平均狀況的模擬Table 3 Simulated 30 years’ average penology conditions of maize at different sowing dates by CERES-Maize model

    表4 不同播期玉米平均產(chǎn)量、單株籽粒數(shù)和百粒重的變化Table 4 Changes of grain yield, kernels number per plant and kernel weight under different sowing dates

    2.3 基于最佳季節(jié)法的玉米最佳播期分析及方法改進

    應用玉米生長發(fā)育進程與最佳季節(jié)同步的方法對當?shù)赜衩咨L發(fā)育過程與氣候資源的匹配情況進行分析,分別計算30a(1981-2010)平均氣候狀態(tài)下玉米抽雄-吐絲期的光溫條件Con1和灌漿期的光溫水條件Con2(表5)。由表5可見,當開花期處于7月28日-8月1日時,(Con1+Con2)值最大,大于1.5。2011-2013年玉米分期播種試驗結果表明,播種-開花期所需≥10℃積溫為1562~1688℃·d,從開花期(7月28日-8月1日)向前推算,則玉米播種期為5月14-26日,該結果與作物模型法對最佳播期的研究結果相比有所差異,且與實際生產(chǎn)中的播期也有較大出入。對郭銀巧等[11]最佳季節(jié)法的進一步分析發(fā)現(xiàn),該方法沒有考慮玉米灌漿期內(nèi)低溫對產(chǎn)量的不利影響。研究表明[20],玉米灌漿期內(nèi)氣溫如果低于16℃,將影響籽粒中淀粉酶的活性,養(yǎng)分的運輸和積累不能正常進行,此期間日平均氣溫在21~22℃的年份籽粒灌漿速度快,成熟早,產(chǎn)量高,而日平均氣溫低于20℃的年份籽粒灌漿速度慢,成熟晚,收獲時籽粒水分含量高,產(chǎn)量相對較低。因此,本研究在應用最佳季節(jié)法時,考慮灌漿期內(nèi)日平均氣溫低于20℃的情況,定義d6為GFD內(nèi)日平均溫度低于20℃的天數(shù),則

    定義(con1+con2-con3)值最大所對應的日期為最佳開花日期。利用改進的最佳季節(jié)法對錦州玉米最佳播期的分析表明,7月20-22日為該地區(qū)最佳開花日期,對應的最佳播期為5月2-13日,該結果與作物模型法的研究趨于一致。

    當開花日期處于7月8-18日,對應的播期為4 月12日-5月5日,此期間(Con1+Con2-Con3)的值在0.8~1.1,低于開花日期為7月22日的(Con1+Con2-Con3)值,主要原因是玉米抽雄-吐絲期氣溫偏低(低于25℃),不利于玉米高產(chǎn)的形成。當播期推遲至5月15日以后時,(Con1+Con2-Con3)值小于1.1,也相對偏低,主要原因在于玉米灌漿期日均溫小于20℃的天數(shù)增多以及降水偏少,也不利于玉米高產(chǎn)的形成。

    表5 基于最佳季節(jié)法不同開花日期的播期設計結果Table 5 Analysis results of different anthesis dates based on optimum season method

    3 結論與討論

    3.1 討論

    Kwabiah[21]通過分期播種試驗研究播期對玉米產(chǎn)量的影響,其結果表明,當播期從5月7日推遲至6月4日,玉米產(chǎn)量的平均下降速率為1.5%·d-1,國外更多的播期影響研究是基于作物模型開展的,如Anapalli等[8]應用CERES-Maize模型和RZWQM (Root Zone Water Quality Model)模型的研究結果表明,播期推遲會使玉米產(chǎn)量下降,且生育期長的品種產(chǎn)量下降速度大于生育期短的品種,CERES-Mazie模型在巴西對作物的模擬顯示,當玉米播期從2月1日推遲至4 月15 日時,雨養(yǎng)地塊和灌溉地塊的平均產(chǎn)量分別下降55%和21%,雨養(yǎng)地塊玉米平均產(chǎn)量的變化呈現(xiàn)先增加后減少的拋物線型[7]。國內(nèi)的播期影響研究多集中于分期播種試驗,如對海城地區(qū)的研究表明,自4月29日開始,隨著玉米播期的推遲,生育進程的加快,各生育時期也相應縮短,產(chǎn)量逐漸降低[22]。CERES-Maize模型中單株籽粒數(shù)由開花吐絲期間單個植株的平均光合速率,以及作物模型參數(shù)單株最大穗粒數(shù)G2共同決定,本研究模擬結果表明,隨著播期的推遲,玉米單株籽粒數(shù)明顯增加,該結果在分期播種試驗中同樣得到證實,劉明等[23]研究表明,晚播(5月15日播種)春玉米穗粒數(shù)較早播(4月24日播種)顯著提高。本研究中8個播期玉米百粒重的變化均呈現(xiàn)先增加后降低的趨勢,CERES-Maize模型中粒重通過最優(yōu)生長速率函數(shù)與灌漿時間的乘積而得,由此可見,5月20日和5月30日為播期的玉米百粒重的降低主要由灌漿時間縮短引起,充足的籽粒灌漿時間有利于玉米獲得較高的產(chǎn)量[24]。從模擬結果看,當播期從4月10日推至5月10日時,玉米平均產(chǎn)量并不降低,這是因為隨著播期的推遲,玉米灌漿成熟期溫光水條件匹配得更好,百粒重增加,因此,玉米平均產(chǎn)量增加。當播期推至5月20日時,單株籽粒數(shù)雖有所增加,但百粒重下降明顯,且灌漿期內(nèi)遇低溫的概率增大,因此玉米的平均產(chǎn)量降低。綜合各個播期玉米的產(chǎn)量及產(chǎn)量組成來看,該地區(qū)4月25日-5月10日播種玉米較適宜,此期間播種玉米,產(chǎn)量無顯著性差異,且能獲得較適宜的單株籽粒數(shù)和百粒重,5月10日前后播種,利于在單株籽粒數(shù)與百粒重之間獲得較好的平衡,玉米平均產(chǎn)量最高。

    使用改進前的最佳季節(jié)法確定的玉米最佳播期(5月20-26日)較滯后,與實際生產(chǎn)不符,本研究在郭銀巧等[11]方法的基礎上,提出了將灌漿期內(nèi)不利溫度條件考慮在內(nèi)的改進方法,對玉米灌漿期的氣候資源進行了重新定量表達,因為溫度對籽粒形成和灌漿有重要影響,在籽粒灌漿成熟期間,不適宜的溫度會影響淀粉酶的活性,使養(yǎng)分的運轉(zhuǎn)和積累不能正常進行,造成結實不飽滿。研究表明,在玉米灌漿期進行平均19℃的低溫處理,玉米的干物質(zhì)積累速度減緩,灌漿速度明顯下降[20]。使用改進后的最佳季節(jié)法確定的最佳播期與作物模型法的研究結果趨于一致,可在今后的研究中進一步應用并檢驗效果。與作物模型法相比,最佳季節(jié)法確定玉米最佳播期僅用到溫度、日照、降水等氣象資料,可省去作物模型法的參數(shù)校準及模型檢驗等繁瑣環(huán)節(jié),易于推廣應用,同時其農(nóng)業(yè)氣象原理也較明確,易于理解。

    本研究中,通過作物模型法和改進的最佳季節(jié)法確定的適宜播期具有較好的一致性,即在4月25 日-5月10日前后播種玉米有利于獲得較高產(chǎn)量,該適宜播期僅代表基于多年(1981-2010年)氣候狀況的分析結果,然而影響玉米最終經(jīng)濟產(chǎn)量的因素較多,如人為管理、病蟲害、當年的氣候條件等,因此確定播期還應結合當年的實際情況因地制宜。由于田間試驗結果的限制,本研究僅考慮了錦州地區(qū)單一玉米品種(晚熟品種)的最佳播期設置,今后將針對其它熟性品種進一步驗證最佳季節(jié)法確定播期的可靠性。

    3.2 結論

    CERES-Maize模型能夠較好地模擬錦州地區(qū)丹玉39不同播期物候和產(chǎn)量狀況,其NRMSE小于10.3%,應用該模型分析不同播期對玉米產(chǎn)量的影響表明,該地區(qū)玉米適宜在4月25日-5月10日播種,此期間播種玉米,其產(chǎn)量無顯著性差異,玉米產(chǎn)量較其它播期平均高出5%。最佳季節(jié)法應用于錦州地區(qū)最佳播期分析時,應對方法加以改進,即需要考慮灌漿期內(nèi)不利溫度條件的影響,利用改進的最佳季節(jié)法得到的最佳播期與作物模型法的研究結果趨于一致??傊?,在得到充分驗證的條件下,作物生長模型可以作為玉米播期影響分析的有力工具;最佳季節(jié)法以保證玉米生長發(fā)育進程與最佳季節(jié)同步為原則,利用該方法確定的最佳播期也具有一定的科學性。

    參考文獻References

    [1]武艷芍,郝建平.不同播期對玉米(強盛49)出苗速度及生育期的影響[J].中國農(nóng)學通報,2009,25(4):119-121. Wu Y S,Hao J P.Effect on the speed of seeding and growth period of Qiangsheng49 in the different sowing date[J]. Chinese Agricultural Science Bulletin,2009,25(4):119-121.(in Chinese)

    [2]馬國勝,薛吉全,路海東,等.播種時期與密度對關中灌區(qū)夏玉米群體生理指標的影響[J].應用生態(tài)學報,2007,18(6): 1247-1253. Ma G S,Xue J Q,Lu H D,et al.Effects of planting date and density on population physiological indices of summer corn in central Shaanxi irrigation area[J].Chinese Journal of Applied Ecology,2007,18 (6):1247-1253.(in Chinese)

    [3]李向嶺,趙明,李從鋒,等.播期和密度對玉米干物質(zhì)積累動態(tài)的影響及其模型的建立[J].作物學報,2010,36(12): 2143-2153. Li X L,Zhao M,Li C F,et al.Effect of sowing-date and planting density on dry matter accumulation dynamic and establishment of its simulated model in Maize[J].Acta Agronomica Sinica,2010,36(12):2143-2153.(in Chinese)

    [4]Kamara A Y,Ekeleme F,Chikoye D.Planting date and cultivar effects on grain yield in dryland corn production[J].Agronomy Journal,2009,101(1):91-98.

    [5]王如芳,張吉旺,呂鵬.多分蘗玉米分蘗發(fā)生規(guī)律及密度和播期的影響[J].作物學報,2012,38(2):322-332 Wang R F,Zhang J W,Lv P.Tillering characteristics of muti-tiller maize and influence of plant density and sowing date[J].Acta Agronomica Sinica,2012,38(2):322-332.(in Chinese)

    [6]楊寧,潘學標,張立禎,等.適宜播期提高農(nóng)牧交錯帶春小麥產(chǎn)量和水氮利用效率[J].農(nóng)業(yè)工程學報,2014,30(8):81-90. Yang N,Pan X B,Zhang L Z,et al.Optimal sowing dates improving yield,water and nitrogen use efficiencies of spring wheat in agriculture and pasture ecotone[J].Transactions of the CSAE,2014,30 (8):81-90.(in Chinese)

    [7]Soler C M T,Sentelhas P C,Hoogenboom G.Application of the CSM-CERES-Maize model for planting date evaluation and yield forecasting for maize grown off-season in a subtropical environment[J].Europ.J.Agronomy,2007,27:165-177.

    [8]Anapalli S S,Ma L,Nielsen D C.Simulating planting date effects on corn production using RZWQM and CERES-Maize models[J].Agron.J.2005,97:58-71.

    [9]Andarzian B,Hoogenboom G,Bannayan M,et al.Determining optimum sowing date of wheat using CSM-CERES-Wheat model[J].Journal of Saudi Society of Agricultural Sciences,2015,(14):189-199.

    [10]薛昌穎,楊曉光,陳懷亮,等.基于ORYZA2000模型的北京地區(qū)旱稻適宜播種期分析[J].生態(tài)學報,2010,30(24): 6970-6979. Xue C Y,Yang X G,Chen H L,et al.Determining suitable sowing dates for aerobic rice in Beijing area using theORYZA2000 mode[J].Acta Ecologica Sinica,2010,30(24): 6970-6979.(in Chinese)

    [11]郭銀巧,郭新宇,趙春江,等.玉米適宜品種選擇和播期確定動態(tài)知識模型的設計與實現(xiàn)[J].中國農(nóng)業(yè)科學, 2006,39(2): 274-280. Guo Y Q,Guo X Y,Zhao C J,et al.Design and implementation of a dynamic knowledge model for suitable cultivars and sowing time of maize[J].Scientia Agricultura Sinica,2006,39(2):274-280.(in Chinese)

    [12]Boote K J,Jones J W,White J W,et al.Putting mechanisms into crop production models[J].Plant,Cell and Environment,2013,(36):1658-1672.

    [13]張建平,王春乙,趙艷霞,等.基于作物模型的低溫冷害對我國東北三省玉米產(chǎn)量影響評估[J].生態(tài)學報,2012,32(13): 4132-4138. Zhang J P,Wang C Y,Zhao Y X,et al.Impact evaluations of low temperature to yields of maize in Northeast China based on crop growth model[J].Acta Ecologica Sinica,2012,32(13): 4132-4138.(in Chinese)

    [14]康蕾,張紅旗.我國五大糧食主產(chǎn)區(qū)農(nóng)業(yè)干旱態(tài)勢綜合研究[J].中國生態(tài)農(nóng)業(yè)學報,2014,22(8):928-937. Kang L,Zhang H Q.Comprehensive research on the state of agricultural drought in five main grain producing areas in China[J].Chinese Journal of Eco-Agriculture,2014,22(8): 928-937.(in Chinese)

    [15]馮玉香,何維勛.我國玉米霜凍害的時空分布[J].中國農(nóng)業(yè)氣象,2000,21(3):6-10. Feng Y X,He W X.The temporal and spatial distribution of maize frost damages in China[J].Chinese Journal of Agrometeorology,2000,21(3):6-10.(in Chinese)

    [16]Hoogenboom G,Jones J W,Wilkens P W,et al.Decision Support System for Agrotechnology Transfer (DSSAT)version 4.5[M].Honolulu:University of Hawaii,2010.

    [17]Jones J W,Hoogenboom G,Porter C H,et al.DSSAT cropping system model[J].European Journal of Agronomy,2003,(18): 235-265.

    [18]Jamieson P D,Porter J R,Wilson D R.A test of the computer simulation model ARC-WHEAT on wheat crops grown in New Zealand[J].Field Crops Res,1991,27:337-350.

    [19]Allen R G,Pereira L S,Raes D,et al.Crop evapotranspiration: guidelines for computing crop requirements[R].Rome:FAO,1998:56.

    [20]李維岳,才卓.吉林玉米[M].長春:吉林科學技術出版社,2000:301. Li W Y,Cai Z.Maize of Jinlin[M].Changchun:Jinlin Scientific and Technical Publishers,2000:301.(in Chinese)

    [21]Kwabiah A B.Growth,maturity,and yield responses of silage maize to hybrid,planting date and plastic mulch[J].Journal of New Seeds,2005,(7):37-59.

    [22]陳曉威.不同播期對玉米生長發(fā)育及產(chǎn)量的影響[J].遼寧農(nóng)業(yè)職業(yè)技術學院學報,2009,11(1):10-11. Chen X W.Effects of different planting dates on maize growth and yield[J].Journal of Liaoning Agricultural College,2009,11(1):10-11.(in Chinese)

    [23]劉明,陶洪斌,王璞,等.播期對春玉米生長發(fā)育、產(chǎn)量及水分利用的影響[J].玉米科學,2009,17(2):108-111. Liu M,Tao Y B,Wang P,et al.Effects of sowing date on growth, yield formation and water utilization of spring maize[J].Journal of Maize Sciences,2009,17(2):108-111.(in Chinese)

    [24]李向嶺,李從鋒,侯玉虹.不同播期夏玉米產(chǎn)量性能動態(tài)指標及其生態(tài)效應[J].中國農(nóng)業(yè)科學,2012, 45(6):1074-1083. Li X L,Li C F,Hou Y H,et al.Dynamic characteristics of summer maize yield performance in different planting dates and its effect of ecological factors[J].Scientia Agicultura Sinica,2012,45(6):1074-1083.(in Chinese)

    Analysis on Optimum Sowing Date of Maize in Jinzhou Using Crop Growth Model and Optimum Season Method

    MI Na, ZHANG Yu-shu, JI Rui-peng, CAI Fu, YU Wen-ying, ZHANG Shu-jie
    (Institute of Atmospheric Environment, China meteorological Administration (Shenyang), Shenyang 110166, China)

    Abstract:To investigate the maize optimum sowing date under the average climate pattern in Jinzhou as well as evaluate the adaptability of crop growth model method and optimum season method, three years’ experimental data collected during 2011-2013 growing season conducted in Jinzhou crop site was employed in this study to validate the cultivar parameters and evaluate the model. Eight planting dates were simulated using 30 years of historical weather data from Jinzhou meteorological observational site to analyze the effect of different sowing dates on maize yield. Optimum season method was also used to investigate the optimum sowing date of maize in this site. An improved method was proposed to improve the ability of optimum season method. The evaluation of CERES-Maize showed that the model was able to simulated penology and grain yield accurately, with normalized RMSE less than 10.3%. The sowing date analysis showed that a delayed planting date from April 30 to May 10 caused an increase in average yield of 6%. Maize median yield decreased from 9112kg·ha-1to 8619 kg·ha-1when planting date delayed from May 10 to May 30. The yield of sowing date at 10 May had no significant difference with that of sowing at 20book=69,ebook=72April and 25 April. The optimum sowing date determined by optimum season analysis method was delayed,comparing to that determined by the crop model method and actual planting date. So an improved method was proposed, which took the unfavorable temperature during grain filling into account. Then the results were accordance with the real situation. In conclusion, the optimum sowing date of maize in this region was from 25 April to 10 May. The results could provide a reference and scientific basis for determining optimum sowing date of maize.

    Key words:CERES-Maize model; Optimum sowing date; Optimum season method

    作者簡介:米娜(1979-),女,副研究員,主要從事農(nóng)業(yè)氣象學研究。E-mail:mina7921@126.com

    基金項目:沈陽大氣環(huán)境研究所公益性科研院所基本科研業(yè)務費專項基金(2014IAE-CMA01);遼寧省農(nóng)業(yè)領域青年科技創(chuàng)新人才培養(yǎng)計劃項目(2014060);遼寧省科學技術廳農(nóng)業(yè)攻關及成果產(chǎn)業(yè)化項目(2014210003);公益性行業(yè)(氣象)科研專項(GYHY201206018)

    * 收稿日期:2015-06-10**通訊作者。E-mail:yushuzhang@126.com

    doi:10.3969/j.issn.1000-6362.2016.01.009

    日日摸夜夜添夜夜添小说| 成人特级av手机在线观看| 欧美性感艳星| av女优亚洲男人天堂| 亚洲精品色激情综合| 欧美不卡视频在线免费观看| 蜜桃久久精品国产亚洲av| 亚洲av免费高清在线观看| 九色成人免费人妻av| 精品一区二区三区视频在线观看免费| 人人妻人人看人人澡| 十八禁国产超污无遮挡网站| 中国美白少妇内射xxxbb| 亚洲av二区三区四区| 九九在线视频观看精品| 亚洲精品一卡2卡三卡4卡5卡| 一级毛片aaaaaa免费看小| 亚洲在线自拍视频| 精品国内亚洲2022精品成人| 午夜视频国产福利| 色综合色国产| 亚洲电影在线观看av| 久久精品91蜜桃| 人人妻人人澡人人爽人人夜夜 | 不卡一级毛片| 国产精品一区二区性色av| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 国产一级毛片七仙女欲春2| 成年版毛片免费区| 白带黄色成豆腐渣| 1000部很黄的大片| av天堂中文字幕网| 老熟妇乱子伦视频在线观看| 一级黄色大片毛片| 黄色一级大片看看| 又黄又爽又免费观看的视频| 亚洲无线观看免费| 噜噜噜噜噜久久久久久91| 99久久精品热视频| 长腿黑丝高跟| 中国美女看黄片| 国产美女午夜福利| 亚洲精品456在线播放app| 观看免费一级毛片| 午夜免费男女啪啪视频观看 | avwww免费| 国产一区二区亚洲精品在线观看| 一级毛片久久久久久久久女| 欧美zozozo另类| 亚洲欧美清纯卡通| 亚洲av电影不卡..在线观看| a级一级毛片免费在线观看| 小说图片视频综合网站| 成人av一区二区三区在线看| 日日撸夜夜添| 最近视频中文字幕2019在线8| 狂野欧美激情性xxxx在线观看| 国产一区二区三区在线臀色熟女| 一级毛片久久久久久久久女| 色播亚洲综合网| 久久午夜福利片| 亚洲四区av| 别揉我奶头 嗯啊视频| 久久草成人影院| 亚洲天堂国产精品一区在线| av专区在线播放| 久久久久国产网址| 亚洲欧美精品自产自拍| 无遮挡黄片免费观看| 国产精品一及| 亚洲av免费在线观看| 五月玫瑰六月丁香| 免费搜索国产男女视频| 欧美精品国产亚洲| 黑人高潮一二区| 久久久久久久午夜电影| 欧美性感艳星| 亚洲最大成人手机在线| 午夜老司机福利剧场| 国产麻豆成人av免费视频| 最近的中文字幕免费完整| 久久久国产成人免费| 给我免费播放毛片高清在线观看| 日本a在线网址| 91在线精品国自产拍蜜月| videossex国产| 69人妻影院| 免费看美女性在线毛片视频| 精品一区二区三区视频在线| 可以在线观看的亚洲视频| 在线看三级毛片| 国产精品美女特级片免费视频播放器| 久久午夜亚洲精品久久| 亚洲av免费高清在线观看| 18禁裸乳无遮挡免费网站照片| 免费人成视频x8x8入口观看| 亚洲国产色片| 亚洲人成网站在线播| 三级毛片av免费| 久99久视频精品免费| 99久久精品国产国产毛片| 97热精品久久久久久| 国产精品1区2区在线观看.| 伦理电影大哥的女人| 日韩国内少妇激情av| 女同久久另类99精品国产91| 国产精品亚洲美女久久久| 97超碰精品成人国产| 免费av观看视频| 午夜免费男女啪啪视频观看 | 婷婷色综合大香蕉| 欧美日韩在线观看h| 精品午夜福利在线看| 久久精品国产清高在天天线| av中文乱码字幕在线| 日韩精品青青久久久久久| 此物有八面人人有两片| 99国产极品粉嫩在线观看| 神马国产精品三级电影在线观看| 成人av一区二区三区在线看| 男人舔女人下体高潮全视频| 69人妻影院| 男人的好看免费观看在线视频| 99在线视频只有这里精品首页| 亚洲欧美中文字幕日韩二区| 亚洲成人久久爱视频| 国产精品美女特级片免费视频播放器| 午夜激情福利司机影院| 简卡轻食公司| 国产黄色视频一区二区在线观看 | 日日啪夜夜撸| 亚洲自拍偷在线| 国产蜜桃级精品一区二区三区| 麻豆国产97在线/欧美| 久久久久久久久大av| 国产精品精品国产色婷婷| 久久久久国内视频| 国产日本99.免费观看| 精品无人区乱码1区二区| 国产色婷婷99| 人人妻人人看人人澡| 亚洲精品一卡2卡三卡4卡5卡| 亚洲真实伦在线观看| 长腿黑丝高跟| 在线观看美女被高潮喷水网站| 男人和女人高潮做爰伦理| 亚洲精品456在线播放app| 尾随美女入室| 国产三级中文精品| 欧美bdsm另类| 久久精品国产自在天天线| 3wmmmm亚洲av在线观看| 日韩中字成人| 看非洲黑人一级黄片| 国产高清激情床上av| 国产一区亚洲一区在线观看| 男女下面进入的视频免费午夜| 男人舔女人下体高潮全视频| 亚洲美女视频黄频| 99久久精品热视频| 久久精品国产清高在天天线| 国产私拍福利视频在线观看| 亚洲婷婷狠狠爱综合网| 又爽又黄a免费视频| 精品福利观看| 男女啪啪激烈高潮av片| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 久久久久免费精品人妻一区二区| 亚洲乱码一区二区免费版| av免费在线看不卡| 国产黄a三级三级三级人| 久久久久久国产a免费观看| 女生性感内裤真人,穿戴方法视频| 乱人视频在线观看| 18禁在线无遮挡免费观看视频 | 国产白丝娇喘喷水9色精品| 黄片wwwwww| 两个人视频免费观看高清| 最近视频中文字幕2019在线8| 国国产精品蜜臀av免费| 亚洲内射少妇av| 国产黄a三级三级三级人| 床上黄色一级片| 99在线人妻在线中文字幕| 国产不卡一卡二| 欧美人与善性xxx| 免费无遮挡裸体视频| 搡老熟女国产l中国老女人| 99久国产av精品国产电影| 偷拍熟女少妇极品色| 亚洲欧美成人综合另类久久久 | 日韩av不卡免费在线播放| 不卡视频在线观看欧美| 直男gayav资源| 国产黄色视频一区二区在线观看 | 91久久精品国产一区二区成人| 欧美最新免费一区二区三区| 啦啦啦啦在线视频资源| 免费大片18禁| av国产免费在线观看| 在线a可以看的网站| 国产精品综合久久久久久久免费| 日韩av不卡免费在线播放| 久久这里只有精品中国| 日韩欧美一区二区三区在线观看| 亚洲精品粉嫩美女一区| 免费搜索国产男女视频| 国产成人aa在线观看| 校园春色视频在线观看| 欧美国产日韩亚洲一区| 国产 一区 欧美 日韩| 可以在线观看毛片的网站| 赤兔流量卡办理| 看黄色毛片网站| 国产精品亚洲美女久久久| 中国美白少妇内射xxxbb| 中文资源天堂在线| 白带黄色成豆腐渣| 插逼视频在线观看| 久久婷婷人人爽人人干人人爱| 亚洲欧美日韩卡通动漫| 亚洲精品乱码久久久v下载方式| 中文字幕av在线有码专区| 特级一级黄色大片| 欧美日韩综合久久久久久| 国内少妇人妻偷人精品xxx网站| 久久久午夜欧美精品| 精品国产三级普通话版| АⅤ资源中文在线天堂| 最近中文字幕高清免费大全6| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 免费观看在线日韩| 美女免费视频网站| 国产精品精品国产色婷婷| 国产麻豆成人av免费视频| 成人综合一区亚洲| 人人妻,人人澡人人爽秒播| 国产精品三级大全| 1024手机看黄色片| 国产aⅴ精品一区二区三区波| 日本成人三级电影网站| 男女啪啪激烈高潮av片| 波多野结衣高清作品| 国产欧美日韩一区二区精品| 国产aⅴ精品一区二区三区波| 黄色一级大片看看| 村上凉子中文字幕在线| 国产av不卡久久| 国产欧美日韩一区二区精品| 啦啦啦观看免费观看视频高清| av在线观看视频网站免费| a级毛片免费高清观看在线播放| 欧美中文日本在线观看视频| 看免费成人av毛片| 两个人的视频大全免费| 亚洲成人精品中文字幕电影| 久久这里只有精品中国| 国产黄a三级三级三级人| 国产黄色小视频在线观看| 国模一区二区三区四区视频| av女优亚洲男人天堂| 一级av片app| 老师上课跳d突然被开到最大视频| 天堂av国产一区二区熟女人妻| 看免费成人av毛片| 三级国产精品欧美在线观看| 欧美日本亚洲视频在线播放| 男女做爰动态图高潮gif福利片| 日韩av在线大香蕉| 亚洲三级黄色毛片| 亚洲国产精品合色在线| 成人特级黄色片久久久久久久| 亚洲精品亚洲一区二区| 日韩三级伦理在线观看| 狠狠狠狠99中文字幕| 国产69精品久久久久777片| 日韩国内少妇激情av| 2021天堂中文幕一二区在线观| 久久精品国产亚洲av香蕉五月| 亚洲av第一区精品v没综合| 在线看三级毛片| 国产成人freesex在线 | 亚洲欧美中文字幕日韩二区| 亚洲av中文字字幕乱码综合| 久久精品国产亚洲av天美| 一个人看视频在线观看www免费| 18禁黄网站禁片免费观看直播| 91久久精品国产一区二区成人| 国产亚洲欧美98| 九九爱精品视频在线观看| 少妇的逼水好多| 亚洲经典国产精华液单| 日本免费一区二区三区高清不卡| 看非洲黑人一级黄片| 99在线人妻在线中文字幕| 亚洲aⅴ乱码一区二区在线播放| 好男人在线观看高清免费视频| 啦啦啦啦在线视频资源| 日本一本二区三区精品| 中文字幕av在线有码专区| av天堂中文字幕网| 国产av在哪里看| 国产私拍福利视频在线观看| 俄罗斯特黄特色一大片| 最近手机中文字幕大全| 欧美丝袜亚洲另类| 欧美xxxx性猛交bbbb| 欧洲精品卡2卡3卡4卡5卡区| 看黄色毛片网站| 熟女电影av网| 亚洲精品日韩av片在线观看| 变态另类丝袜制服| 丝袜美腿在线中文| 午夜激情福利司机影院| 国产精品美女特级片免费视频播放器| 亚洲av不卡在线观看| 精品久久久久久久人妻蜜臀av| 中文亚洲av片在线观看爽| 91av网一区二区| 91久久精品国产一区二区三区| 亚洲图色成人| 人妻夜夜爽99麻豆av| a级毛片a级免费在线| 性欧美人与动物交配| 91在线观看av| 我的老师免费观看完整版| 精品一区二区免费观看| а√天堂www在线а√下载| 一进一出抽搐动态| 淫妇啪啪啪对白视频| 一个人看的www免费观看视频| 国产三级在线视频| 国国产精品蜜臀av免费| 国产亚洲av嫩草精品影院| 久久精品国产亚洲av天美| 日韩国内少妇激情av| av天堂在线播放| 色哟哟·www| 亚洲欧美清纯卡通| 欧洲精品卡2卡3卡4卡5卡区| 国产成人freesex在线 | 国产高清有码在线观看视频| 国产日本99.免费观看| 欧美潮喷喷水| av.在线天堂| 国产成人91sexporn| 91久久精品国产一区二区成人| 成年免费大片在线观看| 中文字幕av在线有码专区| 一进一出抽搐gif免费好疼| 亚洲av美国av| 亚洲精品久久国产高清桃花| 精品午夜福利视频在线观看一区| 中国美白少妇内射xxxbb| 亚洲精品亚洲一区二区| 色哟哟·www| 成人性生交大片免费视频hd| 亚洲在线自拍视频| 成人特级黄色片久久久久久久| 丰满的人妻完整版| 美女内射精品一级片tv| 精品久久久久久成人av| 人妻制服诱惑在线中文字幕| 日本色播在线视频| 国产黄色小视频在线观看| 麻豆av噜噜一区二区三区| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 女生性感内裤真人,穿戴方法视频| 亚洲av第一区精品v没综合| 一个人观看的视频www高清免费观看| 岛国在线免费视频观看| 精品久久久久久久末码| 最近手机中文字幕大全| 可以在线观看的亚洲视频| 成年免费大片在线观看| 午夜精品国产一区二区电影 | 亚洲欧美清纯卡通| 99在线视频只有这里精品首页| 亚洲人与动物交配视频| 国产69精品久久久久777片| 亚洲精品国产av成人精品 | 亚洲不卡免费看| 别揉我奶头 嗯啊视频| 亚洲av成人精品一区久久| 看片在线看免费视频| 国产伦一二天堂av在线观看| 自拍偷自拍亚洲精品老妇| av在线天堂中文字幕| 97超级碰碰碰精品色视频在线观看| 日韩在线高清观看一区二区三区| 久久久午夜欧美精品| 舔av片在线| 国产69精品久久久久777片| 亚洲无线观看免费| 日韩欧美精品v在线| 校园人妻丝袜中文字幕| 欧美成人精品欧美一级黄| 不卡视频在线观看欧美| 3wmmmm亚洲av在线观看| 欧美zozozo另类| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜 | 日韩精品中文字幕看吧| 久久精品久久久久久噜噜老黄 | 熟妇人妻久久中文字幕3abv| 深夜a级毛片| 两个人视频免费观看高清| 欧美极品一区二区三区四区| 国产精品乱码一区二三区的特点| eeuss影院久久| 国产精品久久视频播放| 精品久久久久久久久久久久久| 中文亚洲av片在线观看爽| 十八禁网站免费在线| 99久国产av精品国产电影| 可以在线观看毛片的网站| 亚洲精品日韩在线中文字幕 | 日韩三级伦理在线观看| 麻豆国产av国片精品| 天堂av国产一区二区熟女人妻| 国产av在哪里看| av在线播放精品| 久久精品国产亚洲av天美| 亚洲人成网站在线播放欧美日韩| 免费电影在线观看免费观看| 1024手机看黄色片| 亚洲中文日韩欧美视频| 99国产极品粉嫩在线观看| 久久九九热精品免费| 大香蕉久久网| 99热网站在线观看| 国产亚洲av嫩草精品影院| 日韩一区二区视频免费看| 亚洲真实伦在线观看| 日日干狠狠操夜夜爽| 日本在线视频免费播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品粉嫩美女一区| 久久久久性生活片| 看免费成人av毛片| 色吧在线观看| 国产精品一区www在线观看| 国产成人影院久久av| 久久草成人影院| 欧美日韩综合久久久久久| АⅤ资源中文在线天堂| 麻豆乱淫一区二区| 伦理电影大哥的女人| 亚洲av美国av| 久久精品夜色国产| 国内精品一区二区在线观看| 国产成人91sexporn| 成人av在线播放网站| 伦精品一区二区三区| 草草在线视频免费看| 久久九九热精品免费| 亚洲国产日韩欧美精品在线观看| 一a级毛片在线观看| 春色校园在线视频观看| 亚洲人成网站在线观看播放| 精品人妻视频免费看| 99riav亚洲国产免费| 亚洲人成网站在线观看播放| 看非洲黑人一级黄片| 亚洲最大成人手机在线| 18+在线观看网站| 精品国产三级普通话版| 亚洲欧美日韩高清专用| 欧美人与善性xxx| 国产精品免费一区二区三区在线| 美女 人体艺术 gogo| 精品福利观看| 国产一区二区亚洲精品在线观看| 极品教师在线视频| 日本黄色视频三级网站网址| 亚洲精品国产成人久久av| 欧美高清性xxxxhd video| 麻豆久久精品国产亚洲av| 免费观看人在逋| av在线老鸭窝| 精品人妻视频免费看| 毛片女人毛片| 欧美+日韩+精品| 国产精品,欧美在线| 精品无人区乱码1区二区| 国产高潮美女av| 免费观看的影片在线观看| 在线播放无遮挡| 色综合亚洲欧美另类图片| 狂野欧美白嫩少妇大欣赏| 白带黄色成豆腐渣| 欧美一区二区国产精品久久精品| 国产精品久久电影中文字幕| 国内少妇人妻偷人精品xxx网站| 久久6这里有精品| 99热精品在线国产| 欧美日韩精品成人综合77777| 国产亚洲精品久久久com| 亚洲,欧美,日韩| 国产高清视频在线观看网站| 91精品国产九色| 亚洲成人久久性| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 国产亚洲91精品色在线| 欧美xxxx性猛交bbbb| 久久久久国产精品人妻aⅴ院| 亚洲精品粉嫩美女一区| 欧洲精品卡2卡3卡4卡5卡区| 国产三级中文精品| 久久久久国内视频| 欧美人与善性xxx| 99在线人妻在线中文字幕| 免费黄网站久久成人精品| 此物有八面人人有两片| 午夜免费激情av| 久久久久国产网址| 久久久精品94久久精品| 中文在线观看免费www的网站| 亚洲美女黄片视频| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 舔av片在线| 人妻丰满熟妇av一区二区三区| 在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 在线天堂最新版资源| а√天堂www在线а√下载| 真人做人爱边吃奶动态| 亚洲综合色惰| 女的被弄到高潮叫床怎么办| 97碰自拍视频| 午夜免费激情av| 99久久九九国产精品国产免费| 噜噜噜噜噜久久久久久91| 精品一区二区三区av网在线观看| 五月伊人婷婷丁香| 麻豆乱淫一区二区| 悠悠久久av| 一本久久中文字幕| 老熟妇乱子伦视频在线观看| 听说在线观看完整版免费高清| 色av中文字幕| av天堂在线播放| 亚洲国产精品成人久久小说 | 搡老熟女国产l中国老女人| 哪里可以看免费的av片| 精品久久久久久久久久久久久| 欧美日韩在线观看h| 老司机影院成人| 九九在线视频观看精品| 最近中文字幕高清免费大全6| 日韩欧美在线乱码| 波多野结衣高清作品| 天堂影院成人在线观看| eeuss影院久久| 欧美3d第一页| 男人舔奶头视频| 人妻少妇偷人精品九色| 久久热精品热| 能在线免费观看的黄片| 日本-黄色视频高清免费观看| 午夜福利在线观看吧| 久久久久国产精品人妻aⅴ院| 欧美性感艳星| 国产精品人妻久久久久久| 床上黄色一级片| 亚洲欧美日韩高清在线视频| 99久久中文字幕三级久久日本| 九九爱精品视频在线观看| 美女免费视频网站| 最新在线观看一区二区三区| 嫩草影视91久久| 可以在线观看的亚洲视频| 国产精品一区二区性色av| 国产精品久久视频播放| 麻豆一二三区av精品| 免费av毛片视频| 精品久久久久久成人av| 少妇丰满av| 亚洲欧美精品综合久久99| 久久国产乱子免费精品| 又爽又黄无遮挡网站| 国产真实伦视频高清在线观看| 一进一出抽搐动态| 又粗又爽又猛毛片免费看| 国产精品美女特级片免费视频播放器| 成人av一区二区三区在线看| 最近在线观看免费完整版| 两个人的视频大全免费| 精品一区二区三区视频在线| 1000部很黄的大片| 国产片特级美女逼逼视频| 两个人的视频大全免费| 国产片特级美女逼逼视频| 国产高潮美女av| 最近2019中文字幕mv第一页| 美女大奶头视频| 亚洲av成人精品一区久久| 中国国产av一级| 少妇高潮的动态图| 国产亚洲av嫩草精品影院| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 午夜福利高清视频| 能在线免费观看的黄片| 嫩草影视91久久| 久久99热6这里只有精品| 欧美日韩在线观看h| 成熟少妇高潮喷水视频| 国产高清激情床上av|