• <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久久久久小说| 狠狠狠狠99中文字幕| 老司机午夜十八禁免费视频| 午夜免费鲁丝| 亚洲va日本ⅴa欧美va伊人久久| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产欧洲综合997久久, | 老司机福利观看| 97超级碰碰碰精品色视频在线观看| 成人免费观看视频高清| 欧美黑人欧美精品刺激| 亚洲精品美女久久av网站| 欧美日韩中文字幕国产精品一区二区三区| 黄色视频不卡| av在线天堂中文字幕| cao死你这个sao货| 中国美女看黄片| 久99久视频精品免费| 国产区一区二久久| svipshipincom国产片| 一本大道久久a久久精品| 国产乱人伦免费视频| 国产精品免费视频内射| 亚洲色图av天堂| 一边摸一边做爽爽视频免费| 亚洲国产欧美日韩在线播放| 精品不卡国产一区二区三区| 亚洲国产欧美日韩在线播放| 999久久久国产精品视频| 亚洲全国av大片| 国产精品久久久av美女十八| 免费在线观看完整版高清| 久久精品91无色码中文字幕| 亚洲avbb在线观看| 成人国产综合亚洲| 日韩欧美一区视频在线观看| 久久天躁狠狠躁夜夜2o2o| 99久久综合精品五月天人人| 脱女人内裤的视频| 欧美乱码精品一区二区三区| 国产免费男女视频| 国产成+人综合+亚洲专区| 亚洲一区中文字幕在线| 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 免费av毛片视频| 视频区欧美日本亚洲| 亚洲国产精品成人综合色| 真人一进一出gif抽搐免费| 成人国语在线视频| 国产成人啪精品午夜网站| 久久青草综合色| 国产精品98久久久久久宅男小说| 少妇被粗大的猛进出69影院| 久久久久国内视频| 日韩大尺度精品在线看网址| 天天躁夜夜躁狠狠躁躁| 国产精品自产拍在线观看55亚洲| 欧美黄色淫秽网站| 欧美最黄视频在线播放免费| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利18| 久久亚洲精品不卡| 日韩高清综合在线| 国产精品电影一区二区三区| 日韩欧美国产一区二区入口| 国产亚洲精品久久久久5区| 少妇的丰满在线观看| 一本大道久久a久久精品| 亚洲第一av免费看| 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产在线观看| 亚洲精品在线美女| 19禁男女啪啪无遮挡网站| 在线观看免费午夜福利视频| 中文资源天堂在线| 美女免费视频网站| 精品久久久久久久末码| 校园春色视频在线观看| 欧美日韩精品网址| 国产精品二区激情视频| 免费在线观看黄色视频的| 看片在线看免费视频| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲| 国产成人欧美在线观看| 亚洲欧美精品综合一区二区三区| 亚洲欧美日韩无卡精品| 亚洲中文av在线| 久久久精品欧美日韩精品| 在线观看午夜福利视频| 国产成人一区二区三区免费视频网站| 精品午夜福利视频在线观看一区| 欧美最黄视频在线播放免费| 国产精品国产高清国产av| 老汉色∧v一级毛片| 精华霜和精华液先用哪个| 午夜福利免费观看在线| 国产成人影院久久av| 中出人妻视频一区二区| 一二三四在线观看免费中文在| 免费av毛片视频| 精品欧美国产一区二区三| 亚洲精品av麻豆狂野| 好男人电影高清在线观看| 国产精品自产拍在线观看55亚洲| 日韩欧美在线二视频| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 人人妻人人澡人人看| 韩国av一区二区三区四区| 亚洲国产毛片av蜜桃av| 男女视频在线观看网站免费 | 欧美日韩精品网址| 亚洲精品美女久久av网站| 亚洲人成伊人成综合网2020| 日本熟妇午夜| 男女床上黄色一级片免费看| 国内精品久久久久久久电影| 中文字幕精品免费在线观看视频| 国产免费av片在线观看野外av| a在线观看视频网站| 免费观看精品视频网站| 亚洲国产精品sss在线观看| av超薄肉色丝袜交足视频| 两个人看的免费小视频| 久久久久亚洲av毛片大全| 国产高清有码在线观看视频 | 精品久久久久久久久久免费视频| 亚洲专区国产一区二区| 美女大奶头视频| 国产蜜桃级精品一区二区三区| 国产精品av久久久久免费| 欧美日韩一级在线毛片| www.熟女人妻精品国产| 亚洲午夜理论影院| 色哟哟哟哟哟哟| 午夜福利一区二区在线看| 免费看美女性在线毛片视频| 中文字幕人妻熟女乱码| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区字幕在线| 97碰自拍视频| www.999成人在线观看| 女警被强在线播放| 国产aⅴ精品一区二区三区波| 久久精品国产综合久久久| 身体一侧抽搐| 国产精品,欧美在线| 久久久久久久午夜电影| 18美女黄网站色大片免费观看| 性欧美人与动物交配| 亚洲av成人一区二区三| 老汉色∧v一级毛片| 在线永久观看黄色视频| 国产黄色小视频在线观看| 美女 人体艺术 gogo| 看免费av毛片| 欧美日韩黄片免| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 桃红色精品国产亚洲av| 9191精品国产免费久久| 免费av毛片视频| 操出白浆在线播放| 亚洲午夜精品一区,二区,三区| 欧美黄色片欧美黄色片| 狂野欧美激情性xxxx| 国产亚洲av高清不卡| 亚洲av日韩精品久久久久久密| 精品日产1卡2卡| 亚洲精品国产区一区二| а√天堂www在线а√下载| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 国产精品自产拍在线观看55亚洲| 很黄的视频免费| 欧美激情久久久久久爽电影| 日韩精品免费视频一区二区三区| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 免费av毛片视频| 久久久精品国产亚洲av高清涩受| 中文亚洲av片在线观看爽| 免费高清在线观看日韩| 国产激情久久老熟女| 国产精品99久久99久久久不卡| 日韩欧美 国产精品| 成人精品一区二区免费| 999久久久精品免费观看国产| 亚洲第一av免费看| 美女高潮到喷水免费观看| 欧美人与性动交α欧美精品济南到| 91大片在线观看| 国产亚洲精品久久久久久毛片| 久久婷婷人人爽人人干人人爱| 满18在线观看网站| 特大巨黑吊av在线直播 | 久久中文看片网| 国产亚洲av嫩草精品影院| 观看免费一级毛片| 一级黄色大片毛片| 国产av一区二区精品久久| 国产黄色小视频在线观看| 级片在线观看| 精华霜和精华液先用哪个| 男人操女人黄网站| 一级a爱视频在线免费观看| 亚洲av电影在线进入| 欧美国产日韩亚洲一区| 久久久国产欧美日韩av| 最近在线观看免费完整版| 欧美绝顶高潮抽搐喷水| 视频区欧美日本亚洲| 黄色丝袜av网址大全| 久久久久亚洲av毛片大全| 18禁裸乳无遮挡免费网站照片 | 亚洲久久久国产精品| 性欧美人与动物交配| 身体一侧抽搐| 在线观看午夜福利视频| 搞女人的毛片| 淫秽高清视频在线观看| 91老司机精品| 久久狼人影院| 久久久久久亚洲精品国产蜜桃av| 12—13女人毛片做爰片一| 欧美成人性av电影在线观看| 国产精品永久免费网站| 欧美精品亚洲一区二区| 欧美乱色亚洲激情| 两个人看的免费小视频| 不卡一级毛片| 国产精品久久久av美女十八| 国产精品精品国产色婷婷| 亚洲色图 男人天堂 中文字幕| 日日夜夜操网爽| 亚洲中文日韩欧美视频| 久久久久亚洲av毛片大全| 国产欧美日韩一区二区三| 午夜精品久久久久久毛片777| 久久国产亚洲av麻豆专区| 国产乱人伦免费视频| 欧美午夜高清在线| 午夜精品在线福利| 日本一本二区三区精品| 欧美亚洲日本最大视频资源| 国产一区二区在线av高清观看| 香蕉国产在线看| 成人18禁在线播放| 国产伦一二天堂av在线观看| 国产av在哪里看| 久久中文字幕一级| 女人被狂操c到高潮| 国产精品乱码一区二三区的特点| 欧美在线黄色| 啪啪无遮挡十八禁网站| 午夜久久久久精精品| 国产精品98久久久久久宅男小说| 国产成年人精品一区二区| 午夜精品久久久久久毛片777| 看黄色毛片网站| 国产男靠女视频免费网站| 欧美在线一区亚洲| 欧美日本视频| 亚洲狠狠婷婷综合久久图片| 亚洲 欧美一区二区三区| 国产av一区二区精品久久| 国产97色在线日韩免费| 99国产极品粉嫩在线观看| 俺也久久电影网| 757午夜福利合集在线观看| 一个人免费在线观看的高清视频| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 夜夜夜夜夜久久久久| 欧美亚洲日本最大视频资源| 成人手机av| 亚洲欧美精品综合一区二区三区| 亚洲国产精品成人综合色| 一本综合久久免费| 亚洲五月婷婷丁香| 国产精品电影一区二区三区| 久久精品影院6| 久久午夜综合久久蜜桃| 日韩欧美一区视频在线观看| 国产激情久久老熟女| 亚洲自拍偷在线| 男人操女人黄网站| 欧美亚洲日本最大视频资源| 精品欧美国产一区二区三| 亚洲五月色婷婷综合| 大型av网站在线播放| 亚洲午夜理论影院| 男女床上黄色一级片免费看| 午夜久久久久精精品| 久久精品夜夜夜夜夜久久蜜豆 | 免费在线观看黄色视频的| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久人妻精品电影| 9191精品国产免费久久| 亚洲狠狠婷婷综合久久图片| 亚洲午夜理论影院| 视频区欧美日本亚洲| 中文字幕久久专区| 亚洲久久久国产精品| 欧美大码av| 99re在线观看精品视频| 中出人妻视频一区二区| 在线观看www视频免费| 国产成人精品久久二区二区91| 亚洲av成人一区二区三| 淫秽高清视频在线观看| 91在线观看av| 国产一区二区在线av高清观看| 国产成人av激情在线播放| 亚洲激情在线av| 亚洲国产中文字幕在线视频| 日韩精品免费视频一区二区三区| 午夜久久久在线观看| 午夜影院日韩av| 国产日本99.免费观看| 久久久久久人人人人人| 人人妻人人澡人人看| 老熟妇仑乱视频hdxx| 久久午夜亚洲精品久久| 亚洲成人久久性| 在线观看免费日韩欧美大片| 成人特级黄色片久久久久久久| 老鸭窝网址在线观看| 99国产精品一区二区三区| 亚洲性夜色夜夜综合| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| 观看免费一级毛片| 亚洲精品美女久久久久99蜜臀| 日韩欧美国产在线观看| 午夜福利18| 国产色视频综合| 免费在线观看黄色视频的| 母亲3免费完整高清在线观看| 特大巨黑吊av在线直播 | 久久精品影院6| 国产精品亚洲av一区麻豆| 亚洲电影在线观看av| 中文字幕另类日韩欧美亚洲嫩草| 国内久久婷婷六月综合欲色啪| 伊人久久大香线蕉亚洲五| 大型av网站在线播放| 久久午夜综合久久蜜桃| 国产爱豆传媒在线观看 | 亚洲熟妇熟女久久| 国产亚洲欧美精品永久| 白带黄色成豆腐渣| 一级毛片精品| 亚洲av成人不卡在线观看播放网| 久久中文字幕人妻熟女| 日韩欧美在线二视频| 亚洲 国产 在线| 日韩av在线大香蕉| 淫妇啪啪啪对白视频| 美女大奶头视频| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| АⅤ资源中文在线天堂| 精品久久蜜臀av无| 亚洲一区二区三区色噜噜| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 99re在线观看精品视频| 精品一区二区三区四区五区乱码| 精品不卡国产一区二区三区| 听说在线观看完整版免费高清| 丁香六月欧美| 午夜福利免费观看在线| 一本大道久久a久久精品| 黄频高清免费视频| 777久久人妻少妇嫩草av网站| 色精品久久人妻99蜜桃| 黄网站色视频无遮挡免费观看| 嫁个100分男人电影在线观看| 精品国产亚洲在线| 18禁黄网站禁片午夜丰满| 国内少妇人妻偷人精品xxx网站 | 亚洲七黄色美女视频| 欧美日韩瑟瑟在线播放| 欧美成人一区二区免费高清观看 | 欧美激情久久久久久爽电影| 亚洲精华国产精华精| 亚洲专区国产一区二区| 青草久久国产| 国产一区二区三区视频了| 日本免费a在线| 欧美乱码精品一区二区三区| 亚洲七黄色美女视频| 悠悠久久av| 高清毛片免费观看视频网站| 精品国产亚洲在线| 人人妻人人澡人人看| 日本五十路高清| 国产亚洲欧美98| 啦啦啦免费观看视频1| 国产成人啪精品午夜网站| 欧美午夜高清在线| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 国产熟女xx| 51午夜福利影视在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 久热爱精品视频在线9| 久久久久国内视频| 亚洲成人久久性| 精品无人区乱码1区二区| 在线视频色国产色| 精品午夜福利视频在线观看一区| 国产欧美日韩一区二区三| 欧美一级a爱片免费观看看 | 国产日本99.免费观看| 美女大奶头视频| 国产久久久一区二区三区| 日韩欧美三级三区| 国产激情偷乱视频一区二区| 久久中文字幕一级| 国产精品免费视频内射| 国产欧美日韩一区二区精品| av中文乱码字幕在线| 精品日产1卡2卡| 国产欧美日韩一区二区三| 日韩大尺度精品在线看网址| av在线天堂中文字幕| 色老头精品视频在线观看| 中国美女看黄片| 麻豆久久精品国产亚洲av| x7x7x7水蜜桃| 久久国产精品男人的天堂亚洲| 日韩欧美在线二视频| 久久久久久国产a免费观看| 91老司机精品| 天天添夜夜摸| 欧美日韩精品网址| 真人一进一出gif抽搐免费| 国产精品久久久久久亚洲av鲁大| 一卡2卡三卡四卡精品乱码亚洲| 国产一卡二卡三卡精品| 亚洲精品av麻豆狂野| 成在线人永久免费视频| 国产伦一二天堂av在线观看| 中国美女看黄片| 成人av一区二区三区在线看| 在线看三级毛片| 无遮挡黄片免费观看| www国产在线视频色| 日本 av在线| 我的亚洲天堂| 亚洲人成77777在线视频| 亚洲最大成人中文| 久久九九热精品免费| 日韩欧美 国产精品| 黄色片一级片一级黄色片| 午夜福利18| √禁漫天堂资源中文www| 国产精品一区二区三区四区久久 | 国产男靠女视频免费网站| 一边摸一边抽搐一进一小说| 亚洲成av人片免费观看| 国产伦一二天堂av在线观看| 久久久久久国产a免费观看| 精品第一国产精品| 亚洲精品美女久久久久99蜜臀| 成人一区二区视频在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲中文日韩欧美视频| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 精品久久久久久,| 在线看三级毛片| 韩国精品一区二区三区| 欧美人与性动交α欧美精品济南到| 999久久久精品免费观看国产| 国产亚洲精品久久久久5区| 久久精品亚洲精品国产色婷小说| 精品乱码久久久久久99久播| 国产亚洲欧美精品永久| 国产精品久久久久久亚洲av鲁大| 禁无遮挡网站| 在线观看免费视频日本深夜| 狠狠狠狠99中文字幕| 好看av亚洲va欧美ⅴa在| or卡值多少钱| 在线国产一区二区在线| 露出奶头的视频| 国产又黄又爽又无遮挡在线| 国产三级在线视频| 精品一区二区三区四区五区乱码| 亚洲国产高清在线一区二区三 | 后天国语完整版免费观看| 欧美不卡视频在线免费观看 | 视频在线观看一区二区三区| 韩国av一区二区三区四区| 亚洲av电影在线进入| 久久久久久久久中文| 黄色女人牲交| 国产人伦9x9x在线观看| 91麻豆精品激情在线观看国产| 国产高清videossex| 午夜福利在线在线| av福利片在线| 国产视频一区二区在线看| 给我免费播放毛片高清在线观看| 老汉色∧v一级毛片| 久久精品成人免费网站| 精品少妇一区二区三区视频日本电影| 亚洲国产精品成人综合色| 少妇熟女aⅴ在线视频| 性色av乱码一区二区三区2| 不卡av一区二区三区| 黑丝袜美女国产一区| 久久精品国产综合久久久| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 成人免费观看视频高清| 国产激情偷乱视频一区二区| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 色综合站精品国产| 看片在线看免费视频| 在线天堂中文资源库| 精品卡一卡二卡四卡免费| 免费观看人在逋| 久久久久国产精品人妻aⅴ院| 亚洲七黄色美女视频| 777久久人妻少妇嫩草av网站| 18禁裸乳无遮挡免费网站照片 | 亚洲精品美女久久久久99蜜臀| 夜夜躁狠狠躁天天躁| 国产精品二区激情视频| 可以在线观看毛片的网站| 亚洲中文av在线| 波多野结衣巨乳人妻| 国产av不卡久久| 国产亚洲精品综合一区在线观看 | 久久久久九九精品影院| 亚洲五月色婷婷综合| 国产色视频综合| 久久久久久免费高清国产稀缺| 99在线视频只有这里精品首页| 91字幕亚洲| aaaaa片日本免费| 天天躁狠狠躁夜夜躁狠狠躁| 色av中文字幕| 精品国产一区二区三区四区第35| 黄色丝袜av网址大全| 天堂√8在线中文| 亚洲 欧美一区二区三区| 一个人免费在线观看的高清视频| 久久精品亚洲精品国产色婷小说| 日本a在线网址| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 国产午夜福利久久久久久| 99精品久久久久人妻精品| 神马国产精品三级电影在线观看 | 久久草成人影院| 天天躁狠狠躁夜夜躁狠狠躁| 国产精华一区二区三区| 久久香蕉精品热| 一本久久中文字幕| 少妇粗大呻吟视频| 日本撒尿小便嘘嘘汇集6| www.自偷自拍.com| 大香蕉久久成人网| 天天躁夜夜躁狠狠躁躁| 国产视频内射| 亚洲精品久久成人aⅴ小说| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久av美女十八| 国产精品 欧美亚洲| 黄片小视频在线播放| 中文字幕久久专区| 黑丝袜美女国产一区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久亚洲av毛片大全| 久久精品国产99精品国产亚洲性色| 亚洲片人在线观看| 久久久久久久久免费视频了| 亚洲中文av在线| 天天一区二区日本电影三级| 日韩高清综合在线| 999久久久国产精品视频| 成年免费大片在线观看| 亚洲欧美日韩无卡精品| 欧美另类亚洲清纯唯美| 美女免费视频网站| 精品久久久久久久久久久久久 | 后天国语完整版免费观看| 波多野结衣av一区二区av| 99精品在免费线老司机午夜| av超薄肉色丝袜交足视频| 此物有八面人人有两片| 日韩欧美国产在线观看| 亚洲一码二码三码区别大吗| www.自偷自拍.com|