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

    基于WOFOST模型的中國(guó)主產(chǎn)區(qū)冬小麥生長(zhǎng)過程動(dòng)態(tài)模擬

    2017-07-07 00:43:38黃健熙賈世靈馬鴻元侯英雨
    關(guān)鍵詞:主產(chǎn)區(qū)單產(chǎn)冬小麥

    黃健熙,賈世靈,馬鴻元,侯英雨,何 亮

    (1. 中國(guó)農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院,北京 100083;2. 國(guó)家氣象中心,北京 100081)

    基于WOFOST模型的中國(guó)主產(chǎn)區(qū)冬小麥生長(zhǎng)過程動(dòng)態(tài)模擬

    黃健熙1,賈世靈1,馬鴻元1,侯英雨2,何 亮2

    (1. 中國(guó)農(nóng)業(yè)大學(xué)信息與電氣工程學(xué)院,北京 100083;2. 國(guó)家氣象中心,北京 100081)

    大區(qū)域尺度WOFOST(world food studies)模型的動(dòng)態(tài)模擬是作物模型區(qū)域應(yīng)用的重要基礎(chǔ)。該文以中國(guó)冬小麥主產(chǎn)區(qū)為研究對(duì)象,利用中國(guó)冬小麥主產(chǎn)區(qū)內(nèi) 174個(gè)農(nóng)業(yè)氣象站多年觀測(cè)數(shù)據(jù)以及氣象站點(diǎn)觀測(cè)數(shù)據(jù),重點(diǎn)優(yōu)化WOFOST模型中與品種相關(guān)的積溫參數(shù),即出苗至開花有效積溫與開花至成熟有效積溫。在冬小麥主產(chǎn)區(qū)分區(qū)的基礎(chǔ)上,以2012—2015年氣象數(shù)據(jù)驅(qū)動(dòng)WOFOST模型,在站點(diǎn)尺度進(jìn)行冬小麥的物候期、葉面積指數(shù)(leaf area index,LAI)和單產(chǎn)動(dòng)態(tài)模擬和精度分析。結(jié)果表明:WOFOST模型模擬出苗至開花天數(shù)的決定系數(shù)R2為0.89~0.94,均方根誤差RMSE為7.87~11.52 d,模型模擬開花至成熟天數(shù)的R2為0.63~0.77,RMSE為2.99~4.65 d; 模型模擬LAI的R2為0.70~0.83,RMSE為0.89~1.46 m2/m2;灌溉區(qū)WOFOST模擬的單產(chǎn)精度R2為0.45~0.59,RMSE為734~1 421 kg/hm2;雨養(yǎng)區(qū)WOFOST模擬的單產(chǎn)精度R2為0.48~0.61,RMSE為1 046~1 329 kg/hm2。結(jié)果表明,WOFOST模型在全國(guó)尺度取得了較高模擬精度,為區(qū)域尺度作物模型的農(nóng)業(yè)應(yīng)用提供了堅(jiān)實(shí)的過程模型基礎(chǔ)。

    模型;優(yōu)化;溫度;WOFOST;冬小麥;參數(shù)標(biāo)定;物候期;動(dòng)態(tài)模擬

    黃健熙,賈世靈,馬鴻元,侯英雨,何 亮. 基于 WOFOST模型的中國(guó)主產(chǎn)區(qū)冬小麥生長(zhǎng)過程動(dòng)態(tài)模擬[J]. 農(nóng)業(yè)工程學(xué)報(bào),2017,33(10):222-228. doi:10.11975/j.issn.1002-6819.2017.10.029 http://www.tcsae.org

    Huang Jianxi, Jia Shiling, Ma Hongyuan, Hou Yingyu, He Liang. Dynamic simulation of growth process of winter wheat in main production areas of China based on WOFOST model[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2017, 33(10): 222-228. (in Chinese with English abstract)

    doi:10.11975/j.issn.1002-6819.2017.10.029 http://www.tcsae.org

    0 引 言

    冬小麥?zhǔn)侵袊?guó)的 3大糧食作物之一,主要分布在長(zhǎng)城以南,長(zhǎng)江以北,其種植面積占全國(guó)耕地總面積約18%~27%,糧食作物總面積的18%~24%?;谧魑锷L(zhǎng)模型的方法是開展長(zhǎng)勢(shì)監(jiān)測(cè)與產(chǎn)量估測(cè)的重要技術(shù)手段之一。其區(qū)域尺度的應(yīng)用極大依賴于作物模型的標(biāo)定精度。

    本文選擇WOFOST(world food studies)模型作為作物生長(zhǎng)動(dòng)態(tài)過程模型。WOFOST作物模型由世界糧食研究中心和瓦赫寧根農(nóng)業(yè)大學(xué)共同研發(fā),能夠以天為步長(zhǎng)定量模擬氣象和其他環(huán)境因子影響下的作物生長(zhǎng)過程[1]。在過去的幾十年里,WOFOST模型已經(jīng)在諸多國(guó)家和地區(qū)的多個(gè)領(lǐng)域得到廣泛應(yīng)用。例如產(chǎn)量風(fēng)險(xiǎn)分析、年際間產(chǎn)量變化分析、土壤狀況對(duì)產(chǎn)量的影響、氣象條件對(duì)產(chǎn)量的影響、不同作物品種與耕作制度對(duì)產(chǎn)量的影響、氣象條件對(duì)產(chǎn)量的影響等。且模型對(duì)過程的描述是通用的,可以通過改變參數(shù)模擬不同的地理位置上的不同作物,因此對(duì)作物模型的參數(shù)進(jìn)行標(biāo)定校準(zhǔn),使其適應(yīng)于當(dāng)?shù)氐闹付ㄗ魑?,是進(jìn)行模型區(qū)域應(yīng)用的重要前提。WOFOST能夠模擬潛在生長(zhǎng)、水分脅迫和養(yǎng)分脅迫三種水平[1]。

    國(guó)內(nèi)外學(xué)者在作物模型區(qū)域標(biāo)定和模型應(yīng)用做了許多研究與探索。Wit等[2]利用集合卡爾曼濾波方法將遙感的土壤水分估測(cè)值同化到WOFOST模型中,糾正該模型土壤水平衡誤差,對(duì)西班牙、法國(guó)、意大利和德國(guó)1992?2000年間的冬小麥和玉米進(jìn)行產(chǎn)量估測(cè),結(jié)果表明,數(shù)據(jù)同化明顯改善了 66%地區(qū)的冬小麥產(chǎn)量模擬和56%地區(qū)的玉米產(chǎn)量模擬。Ma等[3]在華北平原冬小麥標(biāo)定與區(qū)域化WOFOST模型,通過葉面積指(leaf area index,LAI)耦合SAIL-PROSPECT模型來模擬土壤調(diào)節(jié)植被指數(shù)(soil adjust vegetation index, SAVI),最小化模擬與合成的SAVI之間的差異,重新初始化出苗日期,結(jié)果表明,該方法在將模擬應(yīng)用到區(qū)域尺度方面具有應(yīng)用潛力。Boogaard等[4]采用WOFOST模型評(píng)估歐盟25個(gè)國(guó)家的秋播小麥產(chǎn)量差距的優(yōu)勢(shì)與限制,結(jié)果表明,WOFOST水分脅迫模式下估算秋播小麥的產(chǎn)量精度較高。張建平等[5-10]利用WOFOST模型分析了氣候變化與低溫冷害對(duì)東北地區(qū)春玉米產(chǎn)量的影響。張素青等[11]在河南省夏玉米主產(chǎn)區(qū)對(duì)WOFOST模型進(jìn)行了校準(zhǔn)與驗(yàn)證。孫琳麗等

    [12]在內(nèi)蒙古河套灌區(qū)玉米種植區(qū)對(duì)WOFOST模型進(jìn)行了適應(yīng)性分析。陳思寧等[13]分析了PyWOFOST模型在東北玉米種植區(qū)的適應(yīng)性。張建平等[14]選擇華北平原冬小麥為研究對(duì)象,在WOFOST模型區(qū)域適應(yīng)性分析與檢驗(yàn)的基礎(chǔ)上,利用區(qū)域化的WOFOST模型模擬結(jié)果實(shí)現(xiàn)干旱災(zāi)害對(duì)作物影響的定量分析與動(dòng)態(tài)評(píng)估。綜上所述,目前作物模型標(biāo)定與驗(yàn)證方面的工作主要集中于單個(gè)站點(diǎn)和若干站點(diǎn)的尺度上。尚未見中國(guó)主產(chǎn)區(qū)尺度WOFOST作物模型標(biāo)定和適應(yīng)性研究報(bào)道,其主要挑戰(zhàn)在于中國(guó)冬小麥主產(chǎn)區(qū)的WOFOST模型輸入?yún)?shù)和初始狀態(tài)的空間變異性。

    本文基于中國(guó)冬小麥主產(chǎn)區(qū)內(nèi)的 174個(gè)農(nóng)業(yè)氣象站點(diǎn)觀測(cè)數(shù)據(jù),在站點(diǎn)尺度,評(píng)估WOFOST模型生育期、葉面積指數(shù)和單產(chǎn)動(dòng)態(tài)模擬與精度。評(píng)價(jià)WOFOST模型在全國(guó)冬小麥主產(chǎn)區(qū)的動(dòng)態(tài)模擬的適應(yīng)性。

    1 材料與方法

    1.1 研究區(qū)與數(shù)據(jù)

    冬小麥主要分布在長(zhǎng)城以南,長(zhǎng)江以北,本文的研究區(qū)選擇中國(guó)主要冬小麥種植區(qū),主要包括河北、山西、江蘇、安徽、山東、河南、湖北、重慶、四川、貴州、云南、陜西、甘肅、寧夏等省區(qū)。

    研究區(qū)內(nèi)共包括174個(gè)農(nóng)業(yè)氣象站點(diǎn),其中包括15個(gè)農(nóng)業(yè)氣象試驗(yàn)站點(diǎn),簡(jiǎn)稱試驗(yàn)站,位置分布如圖 1所示。本研究收集了2011至2015年各站點(diǎn)的冬小麥生育期及單產(chǎn)的觀測(cè)數(shù)據(jù)。其中關(guān)鍵生育期用于積溫參數(shù)的計(jì)算,并對(duì)模擬生育期的驗(yàn)證,單產(chǎn)數(shù)據(jù)用于對(duì)模擬結(jié)果精度的檢驗(yàn)。此外,試驗(yàn)站還提供不同生育期冬小麥根、莖、葉、貯存器官的干物質(zhì)質(zhì)量和實(shí)測(cè)LAI等生長(zhǎng)率參數(shù)。

    圖1 研究區(qū)冬小麥分區(qū)及農(nóng)業(yè)氣象站點(diǎn)分布Fig.1 Winter wheat zone and spatial distribution of agro-meteorological stations in study area

    WOFOST模型的氣象要素采用中科院青藏高原研究所生產(chǎn)的中國(guó)區(qū)域地面氣象數(shù)據(jù)集[15-16],主要包括 7個(gè)要素,近地面氣溫,地表氣壓,近地面空氣比濕,近地面全風(fēng)速,向下短波輻射,向下長(zhǎng)波輻射和降水率??臻g分辨率為 0.1°,數(shù)據(jù)獲取網(wǎng)址為:“http://www. Tpedatabase.cn/portal/MetaDataInfo.jsp?MetaDataId=249369”。選擇2011—2015年的氣象數(shù)據(jù),進(jìn)行要素計(jì)算與格式轉(zhuǎn)換,獲得WOFOST模型所需6個(gè)氣象要素,包括逐日的輻射量、平均水汽壓、日最高溫、日最低溫、風(fēng)速及降水。

    1.2 WOFOST模型參數(shù)的校準(zhǔn)方法

    由于主產(chǎn)區(qū)內(nèi)的冬小麥品種和種植模式存在差異,因此有必要將整個(gè)冬小麥主產(chǎn)區(qū)劃分為相對(duì)均質(zhì)的子區(qū)域進(jìn)行作物模型標(biāo)定。選取了單產(chǎn)水平、土壤類型、氣象條件和種植結(jié)構(gòu)為指標(biāo),采用空間聚類的方法,獲得冬小麥分區(qū)。

    WOFOST模型中輸入?yún)?shù)包括氣象、作物、土壤和田間管理,參數(shù)較多,難以實(shí)現(xiàn)對(duì)每個(gè)參數(shù)的標(biāo)定與校準(zhǔn)。因此,需要根據(jù)WOFOST輸入?yún)?shù)的敏感性和物理含義進(jìn)行分類標(biāo)定與校準(zhǔn),對(duì)于不敏感的參數(shù)或低敏感性的參數(shù),采用WOFOST模型默認(rèn)值或通過文獻(xiàn)查閱確定;對(duì)于與品種有關(guān),敏感性較高且空間變異性較大的參數(shù),先通過觀測(cè)數(shù)據(jù)計(jì)算其取值范圍,再通過優(yōu)化算法確定。

    對(duì)于每個(gè)分區(qū),通常還包含若干農(nóng)業(yè)氣象站點(diǎn)。在站點(diǎn)尺度上,通過每個(gè)農(nóng)氣站點(diǎn)記錄的生育期和鄰近氣象站點(diǎn)觀測(cè)的日平均溫度,計(jì)算出與品種相關(guān)的積溫參數(shù),即出苗至開花有效積溫TSUM1與開花至成熟有效積溫TSUM2。同時(shí),假設(shè)鄰近的2 a冬小麥種植的品種不發(fā)生變化,因此,把 2011—2014冬小麥生育期標(biāo)定的TSUM1和TSUM2參數(shù)值,分別賦予到WOFOST模擬的2012—2015生育期。此外,WOFOST的出苗日期和重要土壤參數(shù)(田間持水量、凋萎系數(shù)和初始可利用水含量)是通過每個(gè)農(nóng)業(yè)氣象站點(diǎn)觀測(cè)給定。對(duì)于每個(gè)試驗(yàn)站點(diǎn),比葉面積(specific leaf area,SLA)根據(jù)試驗(yàn)站不同生育期的干物質(zhì)質(zhì)量和LAI計(jì)算。不同生育期的根、莖、葉、貯存器官的干物質(zhì)量分配系數(shù),也是通過觀測(cè)數(shù)據(jù)計(jì)算獲得。對(duì)于分區(qū)中缺失試驗(yàn)站的情況,則采用最近鄰站點(diǎn)賦值的方法確定。其他參數(shù)值通過文獻(xiàn)查閱[17-27]確定或采用WOFOST默認(rèn)值。表1為WOFOST模型中主要作物參數(shù)校準(zhǔn)值。

    1.3 模型參數(shù)的檢驗(yàn)方法

    模型模擬檢驗(yàn)包括了 2個(gè)部分:散點(diǎn)圖比較以及選擇統(tǒng)計(jì)評(píng)價(jià)指標(biāo)對(duì)模擬值和實(shí)測(cè)值進(jìn)行定量評(píng)價(jià)。散點(diǎn)圖為出苗至開花天數(shù)、開花至成熟天數(shù)、單產(chǎn)及LAI模擬值與實(shí)測(cè)值的回歸分析圖。統(tǒng)計(jì)評(píng)價(jià)指標(biāo)選擇了決定系數(shù)(R2)、一致性系數(shù)(D)、變異系數(shù)(coefficient of variation,CV)、均方根誤差(root mean square error,RMSE)、歸一化均方根誤差(normalized root mean square error,NRMSE)。其中R2和D反映了實(shí)測(cè)值與模擬值之間的一致性,越接近1表示模擬效果越好。CV反映了數(shù)值離散程度,值越大越能體現(xiàn)數(shù)據(jù)的空間變異性[28],可將其進(jìn)行粗略分級(jí):CV<10%,為弱變異性;10%≤CV≤100%;為中等變異性,CV>100%,為強(qiáng)變異性[29]。RMSE和 NRMSE反映了模擬值與實(shí)測(cè)值之間的相對(duì)誤差和絕對(duì)誤差[30],值越小表示模擬效果越好,其中,NRMSE≤10%,為極高精度;10%30%,為低精度[11]。D、CV、RMSE和NRMSE的計(jì)算公式如式(1)~(4)。

    式中i表示第i個(gè)樣本;Yi和Xi分別為第i個(gè)樣本模擬值和實(shí)測(cè)值;為全部樣本實(shí)測(cè)數(shù)據(jù)平均值;n為樣本數(shù);SD為模擬結(jié)果的標(biāo)準(zhǔn)差,為全部樣本模擬結(jié)果的平均值。

    表1 WOFOST模型中主要作物參數(shù)校準(zhǔn)值范圍Table 1 Range of calibrated values of main crop parameters of WOFOST model

    2 結(jié)果與分析

    2.1 模型參數(shù)的校準(zhǔn)結(jié)果

    根據(jù)上述模型參數(shù)校準(zhǔn)方法,進(jìn)行模型的校準(zhǔn)。表1為所有冬小麥分區(qū)的部分關(guān)鍵參數(shù)校準(zhǔn)值范圍。

    2.2 WOFOST模型檢驗(yàn)

    為驗(yàn)證WOFOST模型在中國(guó)冬小麥主產(chǎn)區(qū)的動(dòng)態(tài)模擬精度,在站點(diǎn)尺度,以2012—2015年當(dāng)年的站點(diǎn)觀測(cè)出苗日期為模擬初始日期,以氣象、土壤、作物等參數(shù)驅(qū)動(dòng)WOFOST模型進(jìn)行冬小麥生長(zhǎng)模擬,并對(duì)模型模擬的出苗至開花天數(shù)、開花至成熟的天數(shù)、LAI和單產(chǎn)進(jìn)行模擬結(jié)果精度分析與驗(yàn)證。

    2.2.1 生育期驗(yàn)證

    開花期和成熟期分別是冬小麥營(yíng)養(yǎng)生長(zhǎng)和生殖生長(zhǎng)階段的結(jié)束日期,是評(píng)價(jià)WOFOST模型模擬的重要生育期。該文分別選擇出苗期至開花期的天數(shù)和開花期至成熟期的天數(shù)來進(jìn)行關(guān)鍵生育期的驗(yàn)證。

    2012—2015年,模型對(duì)生育期天數(shù)的模擬,具有較為相似的模擬精度。由表2可知,出苗至開花的R2在0.89以上,D在0.96以上,說明模擬值與實(shí)測(cè)值具有較好的一致性,NRMSE在7%以下,模擬誤差在7.87~11.52 d之間,表明WOFOST模型能準(zhǔn)確模擬冬小麥開花期。開花至成熟天數(shù)的R2位于0.63與0.77之間,D在0.87~0.93之間,NRMSE在8%~12%之間,模擬誤差在2.99到4.65 d之間。不同熱量條件的地區(qū),開花到成熟期的天數(shù)差異較大。模擬誤差主要依賴于開花期的誤差和開花到成熟期的積溫精度。以2012年為例(圖2),模型模擬出苗至開花天數(shù)、開花至成熟天數(shù)分別與實(shí)測(cè)值之間具有較好的相關(guān)性,各點(diǎn)均勻的分布在回歸線兩側(cè)。同時(shí),出苗至開花天數(shù)與開花至成熟天數(shù)的 CV均在 10%以上,具有顯著的空間變異性,能充分解釋模擬冬小麥生育期的區(qū)域空間變異。

    表2 不同年份WOFOST模擬生育期的驗(yàn)證結(jié)果精度對(duì)比(2012—2015)Table 2 Comparison of simulated growth stages accuracies in different years (2012—2015)

    圖2 WOFOST模型模擬出苗到開花期天數(shù)和開花到成熟期天數(shù)對(duì)比(2012)Fig.2 Comparison of simulated and measured days from emergence to anthesis and anthesis to maturity (2012)

    2.2.2 LAI的驗(yàn)證

    由圖3可知,2012—2015年模擬LAI與實(shí)測(cè)值之間的R2在0.70~0.83之間,D在0.88~0.96之間,WOFOST模擬LAI值與實(shí)測(cè)值之間的一致性較好,RMSE在0.89~1.46 m2/m2之間,NRMSE在50%~63%之間。由敏感性分析結(jié)果可知,對(duì)LAI最大值敏感的參數(shù)主要有葉片最大CO2同化速率、SLATB、初始生物量和葉片衰老系數(shù)[19]。本研究中,TDWI和SPAN都采用了默認(rèn)值,AMAXTB雖然根據(jù)參考文獻(xiàn)確定,但有一些冬小麥種植區(qū)缺少觀測(cè)數(shù)據(jù),采取模型默認(rèn)值,導(dǎo)致某些站點(diǎn)的WOFOST模擬誤差較大。

    圖3 WOFOST模型模擬LAI值與實(shí)測(cè)值對(duì)比(2012—2015)Fig.3 Comparison of WOFOST simulated and field-measured LAI (2012—2015)

    2.2.3 單產(chǎn)的驗(yàn)證

    考慮到冬小麥光熱和降雨條件的差異,將主產(chǎn)區(qū)劃分為灌溉區(qū)和雨養(yǎng)區(qū)分別進(jìn)行WOFOST模擬,其中黃淮海設(shè)定灌溉區(qū),采用WOFOST潛在模式,西北和西南地區(qū)設(shè)定為雨養(yǎng)區(qū),采用WOFOST的水分脅迫模式。

    灌溉區(qū)WOFOST模擬的單產(chǎn)精度R2為0.45~0.59,RMSE為734~1 421 kg/hm2;雨養(yǎng)區(qū)WOFOST模擬的單產(chǎn)精度R2為0.48~0.61,RMSE為1 046~1 329 kg/hm2。相比較而言,WOFOST模型模擬的灌溉區(qū)單產(chǎn)精度總體要高于雨養(yǎng)區(qū),具有更低的RMSE值(表3)。

    表3 灌溉區(qū)不同年份WOFOST模型模擬單產(chǎn)的驗(yàn)證結(jié)果精度對(duì)比(2012—2015)Table 3 Comparison of simulated yield accuracies in multiple years in irrigation area (2012—2015)

    對(duì)于某些年份單產(chǎn)偏低的可能原因是設(shè)定的品種相關(guān)的參數(shù) TSUM1和 TSUM2等具有較大誤差,導(dǎo)致WOFOST模擬的生育期和產(chǎn)量誤差較大。由表3可知,灌溉區(qū)與雨養(yǎng)區(qū)的D位于0.73~0.99之間,說明模擬值與實(shí)測(cè)值具有很好的一致性??傮w而言,而單產(chǎn)模擬的誤差主要在于,對(duì)單產(chǎn)敏感的參數(shù)難以獲得準(zhǔn)確的空間分布值。

    圖4 2012和2015年WOFOST模型模擬單產(chǎn)和實(shí)測(cè)單產(chǎn)相對(duì)誤差的空間分布圖Fig.4 Spatial distribution of relative error of WOFOST simulated and field-measured yield per unit in 2012 and 2015

    從單產(chǎn)的空間分布差異來看(圖 4),模擬單產(chǎn)精度較高的站點(diǎn)主要分布于黃淮海灌溉區(qū)。而模擬單產(chǎn)相對(duì)誤差較大的點(diǎn),主要集中在雨養(yǎng)區(qū)。可能原因是農(nóng)業(yè)氣象站點(diǎn)分布稀疏和關(guān)鍵土壤參數(shù)難以準(zhǔn)確標(biāo)定。灌溉區(qū)模擬單產(chǎn)的 CV在 14%~22%之間,雨養(yǎng)區(qū)模擬單產(chǎn)的CV在25%~40%之間,能解釋模擬冬小麥單產(chǎn)的空間變異性,能解釋模擬冬小麥單產(chǎn)的空間變異性。

    3 討 論

    由于中國(guó)冬小麥主產(chǎn)區(qū)光熱條件和種植品種的差異,WOFOST模型模擬的生育期空間差異性較大,國(guó)家農(nóng)氣站點(diǎn)記錄的觀測(cè)數(shù)據(jù)表明,主產(chǎn)區(qū)內(nèi)冬小麥出苗到開花期歷時(shí)天數(shù)在89~200 d之間,開花到成熟期時(shí)天數(shù)在20~70 d之間,參數(shù)校準(zhǔn)后的WOFOST模型能準(zhǔn)確撲捉這一差異;WOFOST模型模擬LAI的誤差主要來源于對(duì)LAI敏感的WOFOST模型輸入?yún)?shù)(例如:TDWI和SPAN)未考慮參數(shù)的空間變異性;WOFOST單產(chǎn)模擬方面,我們考慮了雨養(yǎng)區(qū)和灌溉區(qū)的差異,其他影響因素未予考慮,例如:營(yíng)養(yǎng)和病蟲害脅迫等。研究表明,合理的標(biāo)定WOFOST模型的順序?yàn)樯?,LAI和單產(chǎn)。今后的研究將通過衛(wèi)星遙感數(shù)據(jù)和作物模型數(shù)據(jù)同化獲得WOFOST的關(guān)鍵輸入?yún)?shù)的空間分布的優(yōu)化值,從而進(jìn)一步提高大區(qū)域作物模型的模擬能力。標(biāo)定后的WOFOST模型將為區(qū)域尺度的溫度脅迫或水分脅迫對(duì)產(chǎn)量的影響提供動(dòng)態(tài)過程模型。

    4 結(jié) 論

    本文以WOFOST為動(dòng)態(tài)生長(zhǎng)模型,中國(guó)冬小麥主產(chǎn)區(qū)為研究對(duì)象,在分區(qū)的基礎(chǔ)上,基于農(nóng)業(yè)氣象站點(diǎn)觀測(cè)數(shù)據(jù)標(biāo)定WOFOST模型的敏感參數(shù),在站點(diǎn)尺度,動(dòng)態(tài)模擬生育期、LAI和單產(chǎn)。驗(yàn)證結(jié)果表明,模型模擬出苗-開花天數(shù)的NRMSE在4%~7%之間,模型模擬開花-成熟天數(shù)的NRMSE在8%~12%之間,具有較高的模擬精度,CV在14%~20%之間,具有空間變異性。模型模擬的LAI的NRMSE在50%~63%之間。模型模擬單產(chǎn)的NRMSE在11%~28%之間,CV在14%~40%之間,能較好地體現(xiàn)單產(chǎn)的空間差異性??傮w來說,WOFOST模型取得了較為理想的模擬精度,具有較好的適應(yīng)性。

    [1]Diepen C A, Wolf J, Keulen H, et al. WOFOST: A simulation model of crop production[J]. Soil Use &Management, 1989, 5(1): 16-24.

    [2]Wit A J W De, Diepen C A Van. Crop model data assimilation with the Ensemble Kalman filter for improving regional crop yield forecasts[J]. Agricultural & Forest Meteorology, 2007, 146(1-2): 38-56.

    [3]Ma Yuping, Wang Shili, Zhang Li, et al. Monitoring winter wheat growth in North China by combining a crop model and remote sensing data[J]. International Journal of Applied Earth Observation & Geoinformation, 2008, 10(4): 426-437.

    [4]Boogaard Hendrik, Wolf Joost, Supit Iwan, et al. A regional implementation of WOFOST for calculating yield gaps of autumn-sown wheat across the European Union[J]. Field Crops Research, 2013, 143(143): 130-142.

    [5]張建平,王春乙,楊曉光,等. 溫度導(dǎo)致的我國(guó)東北三省玉米產(chǎn)量波動(dòng)模擬[J]. 生態(tài)學(xué)報(bào),2009,29(10):5516-5522.Zhang Jianping, Wang Chunyi, Yang Xiaoguang, et al.Simulation of yield fluctustion caused by the temperature in Northeast China[J]. Acta Ecologica Sinica, 2009, 29 (10):5516-5522. (in Chinese with English abstract)

    [6]張建平,王春乙,趙艷霞,等. 基于作物模型的低溫冷害對(duì)我國(guó)東北三省玉米產(chǎn)量影響評(píng)估[J]. 生態(tài)學(xué)報(bào),2012,32 (13):4132-4138.Zhang Jianping, Wang Chunyi, Zhao Yanxia, et al. Impact evaluation 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 with English abstract)

    [7]陳振林,張建平,王春乙,等. 應(yīng)用 WOFOST模型模擬低溫與干旱對(duì)玉米產(chǎn)量的綜合影響[J]. 中國(guó)農(nóng)業(yè)氣象,2007,28(4):440-442.Chen Zhenlin, Zhang Jianping, Wang Chunyi, et al.Application of WOFOST model in simulation of integrated impacts of low temperature and drought on maize yield[J].Chinese Journal of Agrometeorology, 2007, 28(4): 440-442.(in Chinese with English abstract)

    [8]張建平,趙艷霞,王春乙,等. 氣候變化情景下東北地區(qū)玉米產(chǎn)量變化模擬[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2008,16(6):1448-1452.Zhang Jianping, Zhao Yanxia, Wang Chunyi, et al.Simulation of maize production under climate change scenario in Northeast China[J]. Chinese Journal of Eco-Agriculture, 2008, 16(6): 1448-1452. (in Chinese with English abstract)

    [9]張建平,趙艷霞,王春乙,等. 不同時(shí)段低溫冷害對(duì)玉米灌漿和產(chǎn)量的影響模擬[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2012, 40(9):115-121.Zhang Jianping, Zhao Yanxia, Wang Chunyi, et al. Modeling the impact of low temperature disaster during different development stages on grain filling processes and yields of maize[J]. Journal of Northwest A&F University (Nat.Sci.Ed.),2012, 40(9): 115-121. (in Chinese with English abstract)

    [10]張建平,趙艷霞,王春乙,等. 不同發(fā)育期低溫冷害對(duì)玉米灌漿和產(chǎn)量影響模擬[J]. 中國(guó)農(nóng)學(xué)通報(bào),2012, 28(36):176-182.Zhang Jianping, Zhao Yanxia, Wang Chunyi, et al. The simulation of the effects of low temperature disasters on maize growth and yield during different development stages[J]. Chinese Agricultural Science Bulletin, 2012,28(36): 176-182. (in Chinese with English abstract)

    [11]張素青,張建濤,李繼蕊,等. WOFOST模型在河南省夏玉米主產(chǎn)區(qū)的校準(zhǔn)與驗(yàn)證[J]. 河南農(nóng)業(yè)科學(xué),2014, 43(8):152-156.Zhang Suqing, Zhang Jiantao, Li Jirui, et al. Calibration andvalidation of WOFOST in main maize-producing regions in Henan[J]. Journal of Henan Agricultural Sciences, 2014,43(8): 152-156. (in Chinese with English abstract)

    [12]孫琳麗,侯瓊,馬玉平,等. WOFOST模型在內(nèi)蒙古河套灌區(qū)模擬玉米生長(zhǎng)全程的適應(yīng)性[J]. 生態(tài)學(xué)雜志,2016,35(3):1-10.Sun Linli, Hou Qiong, Ma Yuping, et al. Adaptability of WOFOST model to simulate the whole growth period of maize in Hetao irrigation region of Inner Mongolia[J].Chinese Journal of Ecology, 2016, 35(3): 1-10. (in Chinese with English abstract)

    [13]陳思寧,趙艷霞,申雙和,等. 基于 PyWOFOST作物模型的東北玉米估產(chǎn)及精度評(píng)估[J]. 中國(guó)農(nóng)業(yè)科學(xué),2013,46(14):2880-2893.Chen Sining, Zhao Yanxia, Shen Shuanghe, et al. Study on maize yield estimation and accuracy assessment based on PyWOFOST crop model in Northeast China[J]. Scientia Agricultura Sinica, 2013, 46(14): 2880-2893. (in Chinese with English abstract)

    [14]張建平,趙艷霞,王春乙,等. 基于 WOFOST作物生長(zhǎng)模型的冬小麥干旱影響評(píng)估技術(shù)[J]. 生態(tài)學(xué)報(bào),2013,33(6):1762-1769.Zhang Jianping, Zhao Yanxia, Wang Chunyi, et al.Evaluation technology on drought disaster to yields of winter wheat based on WOFOST crop growth model[J]. Acta Ecologica Sinica, 2013, 33(6): 1762-1769. (in Chinese with English abstract)

    [15]Chen Yingying, Yang Kun, He Jie, et al. Improving land surface temperature modeling for dry land of China[J].Journal of Geophysical Research Atmospheres, 2011,116(116): 999-1010.

    [16]Yang Kun, He Jie, Tang Wenjun, et al. On downward shortwave and longwave radiations over high altitude regions:Observation and modeling in the Tibetan Plateau[J].Agricultural and Forest Meteorology, 2010, 150(1): 38-46.

    [17]Huang Jianxi, Sedano Fernando, Huang Yanbo, et al.Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to improve regional winter wheat yield estimation[J]. Agricultural and Forest Meteorology,2016(216): 188-202.

    [18]Huang Jianxi, Tian Liyan, Liang Shunlin, et al. Improving winter wheat yield estimation by assimilation of the leaf area index from Landsat TM and MODIS data into the WOFOST model[J]. Agricultural and Forest Meteorology, 2015(204):106-121.

    [19]Ma Guannan, Huang Jianxi, Wu Wenbin, et al. Assimilation of MODIS-LAI into the WOFOST model for forecasting regional winter wheat yield[J]. Mathematical & Computer Modelling, 2013, 58(3–4): 634-643.

    [20]Ma Hongyuan, Huang Jianxi, Zhu Dehai, et al. Estimating regional winter wheat yield by assimilation of time series of HJ-1 CCD NDVI into WOFOST-ACRM model with Ensemble Kalman filter[J]. Mathematical & Computer Modelling, 2013, 58(3–4): 759-770.

    [21]Huang Jianxi, Ma Hongyuan, Su Wei, et al. Jointly assimilating MODIS LAI and ET products into the SWAP model to estimate winter wheat yield[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing. 20158(8): 4060-4071.

    [22]Tian Liyan, Li Zhongxia, Huang Jianxi, et al. Comparison of two optimization algorithms for estimating regional winter wheat yield by integrating MODIS leaf area index and world food studies model[J]. Sensor Letters. 2013, 11(6–7): 1261-1268.

    [23]何亮,侯英雨,趙剛,等. 基于全局敏感性分析和貝葉斯方法的 WOFOST作物模型參數(shù)優(yōu)化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(2):169-179.He Liang, Hou Yingyu, Zhao Gang, et al. Parameters optimization of WOFOST model by integration of global sensitivity analysis and Bayesian calibration method[J].Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2016, 32(2): 169-179. (in Chinese with English abstract)

    [24]黃健熙,馬鴻元,田麗燕,等. 基于時(shí)間序列 LAI和 ET同化的冬小麥遙感估產(chǎn)方法比較[J]. 農(nóng)業(yè)工程學(xué)報(bào),2015,31(4):197-203.Huang Jianxi, Ma Hongyuan, Tian Liyan, et al. Comparison of remote sensing yield estimation methods for winter wheat based on assimilating time-sequence LAI and ET[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(4): 197-203.(in Chinese with English abstract)

    [25]黃健熙,武思杰,劉興權(quán),等. 基于遙感信息與作物模型集合卡爾曼濾波同化的區(qū)域冬小麥產(chǎn)量預(yù)測(cè)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012(04):142-148.Huang Jianxi, Wu Sijie, Liu Xingquan, et al. Regional winter wheat yield forecasting based on assimilation of remote sensing data and crop growth model with Ensemble Kalman method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE),2012(04): 142-148. (in Chinese with English abstract)

    [26]邱美娟,宋迎波,王建林,等. 山東省冬小麥產(chǎn)量動(dòng)態(tài)集成預(yù)報(bào)方法[J]. 應(yīng)用氣象學(xué)報(bào),2016, 27(2):191-200.Qiu Meijuan, Song Yingbo, Wang Jianlin, et al. Integrated technology of yield dynamic prediction of winter wheat in Shandong province[J]. Journal of applied meteorological science, 2016, 27(2): 191-200.(in Chinese with English abstract)

    [27]Wu Dingrong, Yu Qiang, Lu Changhe, et al. Quantifying production potentials of winter wheat in the North China Plain[J]. European Journal of Agronomy, 2006, 24(3): 226-235.

    [28]Cheng Zhiqiang, Meng Jihua, Wang Yiming. Improving spring maize yield estimation at field scale by assimilating time-series HJ-1 CCD data into the WOFOST model using a New Method with Fast Algorithms[J]. Remote Sensing, 2016,8(4): 303.

    [29]Mallants D, Mohanty P B, Jacques D, et al. Spatial Variability Of Hydraulic Properties In A Multi-Layered Soil Profile[J]. Soil Science, 1996, 161(3): 167-181.

    [30]戴彤,王靖,赫迪,等. APSIM模型在西南地區(qū)的適應(yīng)性評(píng)價(jià):以重慶冬小麥為例[J]. 應(yīng)用生態(tài)學(xué)報(bào),2015, 26(4):1237-1243.Dai Tong, Wang Jing, He Di, et al. Adaptability of APSIM model in Southwestern China: A case study of winter wheat in Chongqing City. Chinese Journal of Applied Ecology,2015, 26(4): 1237-1243. (in Chinese with English abstract)

    Dynamic simulation of growth process of winter wheat in main production areas of China based on WOFOST model

    Huang Jianxi1, Jia Shiling1, Ma Hongyuan1, Hou Yingyu2, He Liang2
    (1.College of Information and Electrical Engineering, China Agricultural University, Beijing100083,China;2.National Meteorological Center, Beijing100081, China)

    Crop model calibration and parameterization are essential for model evaluation and agricultural application. It is important for model application to accurately estimate the values of crop model parameters and further improve the performance of model prediction. WOFOST (world food studies) is a well-known, widely applied simulation model to analyze quantitatively the growth and production of field crops, which was originally developed for crops in European countries. It is the base model for Monitoring Agricultural Resources (MARS) Crop Growth Monitoring System (CGMS) in operational use for yield estimation in European Union. Dynamic simulation of WOFOST model in large regional scale is an important basis for regional crop modeling. In this study, we selected the main winter wheat production areas of China as the study area, and the data from 174 agricultural meteorological stations from 2011 to 2014 were used to calibrate several key WOFOST input parameters, especially 2 parameters related with variety, namely the effective accumulated temperature from emergence to flowering (TSUM1) and the effective accumulated temperature from flowering to maturity (TSUM2). On the basis of the zoning of the main winter wheat production areas, we used the meteorological data from 2012 to 2015 to drive the WOFOST model at a single-point scale, to simulate the winter wheat growth and dynamic development. The simulated phenology, LAI(leaf area index) and yield at the station level were evaluated with the field measured data. Results showed that the NRMSE(normalized root mean square error) of LAI ranged from 50% to 63%. The NRMSE of simulated days was 4%-7% from emergence to anthesis period and 8%-12% from anthesis to maturity period, and then CV (coefficient of variation) of the phenology was between 14% and 20%, which meant significant spatial variability. We simulated the yield respectively in irrigated area and rainfed area. And the NRMSE of simulated yield in irrigated area ranged from 11% to 23%, while the NRMSE of simulated yield in rain-fed area ranged from 22% and 28%, and the CV ranged from 14% to 22% for irrigated areas and from 25% to 40% for rain-fed areas, which exhibited significant spatial variability. The NRMSE of simulated LAI was between 50% and 63%, which could be explained that the LAI during different growth stages was all included into the accuracy analysis. Several important input parameters (such as TDWI (initial biomass) and SPAN (leaf senescence coefficient))could be optimized through assimilating remote sensing data into crop model, which could greatly improve the performance of crop model at the regional scale. Our results showed that the WOFOST model is of great potential for simulating the dynamic growth process of winter wheat in China. The calibrated WOFOST provides the dynamic model basis for regional applications,such as assimilating remote sensing data into crop model for crop yield estimation and climate change prediction with crop model.

    models; optimization; temperature; WOFOST; winter wheat; parameter calibration; phendogy; dynamic simulation

    10.11975/j.issn.1002-6819.2017.10.029

    S127

    A

    1002-6819(2017)-10-0222-07

    2016-10-07

    2017-05-05

    國(guó)家自然科學(xué)基金(41671418,41471342,41371326)

    黃健熙,博士,博士生導(dǎo)師,主要從事農(nóng)業(yè)定量遙感研究。北京中國(guó)農(nóng)業(yè)大學(xué),100083。Email:jxhuang@cau.edu.cn

    猜你喜歡
    主產(chǎn)區(qū)單產(chǎn)冬小麥
    主產(chǎn)區(qū)小麥?zhǔn)召忂M(jìn)度過七成
    農(nóng)大農(nóng)企聯(lián)手創(chuàng)山西小麥最高單產(chǎn)新紀(jì)錄
    油菜“不務(wù)正業(yè)”,單產(chǎn)3.4噸
    單產(chǎn)948.48千克!“金種子”迸發(fā)大能量
    我國(guó)玉米單產(chǎn)紀(jì)錄第七次被刷新
    全國(guó)麥?zhǔn)者M(jìn)度過半 主產(chǎn)區(qū)機(jī)收率超98%
    關(guān)于支持主產(chǎn)區(qū)發(fā)展專用糧食生產(chǎn)的建議
    甘肅冬小麥田
    冬小麥和春小麥
    中學(xué)生(2015年4期)2015-08-31 02:53:50
    冬小麥——新冬18號(hào)
    制服诱惑二区| 精品国产亚洲在线| 久久久国产成人免费| 18禁裸乳无遮挡免费网站照片| 久久婷婷成人综合色麻豆| 国产精品野战在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲五月天丁香| 天堂动漫精品| 一边摸一边做爽爽视频免费| 看免费av毛片| 在线观看午夜福利视频| 香蕉av资源在线| 一边摸一边做爽爽视频免费| 精品第一国产精品| 男人的好看免费观看在线视频 | 51午夜福利影视在线观看| 日本 av在线| 亚洲中文字幕日韩| 中文字幕高清在线视频| 很黄的视频免费| 色播亚洲综合网| 日韩欧美一区二区三区在线观看| 色哟哟哟哟哟哟| 亚洲精品av麻豆狂野| 色综合欧美亚洲国产小说| 国产精华一区二区三区| 九色成人免费人妻av| 12—13女人毛片做爰片一| 国产午夜精品论理片| 亚洲成人国产一区在线观看| 不卡一级毛片| 老司机深夜福利视频在线观看| 欧美久久黑人一区二区| www国产在线视频色| av有码第一页| 最近最新中文字幕大全免费视频| 美女大奶头视频| 在线观看www视频免费| 黑人操中国人逼视频| 亚洲国产精品久久男人天堂| 三级毛片av免费| 中亚洲国语对白在线视频| 十八禁网站免费在线| 国产真人三级小视频在线观看| 999精品在线视频| www日本在线高清视频| 国产av麻豆久久久久久久| √禁漫天堂资源中文www| 18禁黄网站禁片午夜丰满| 亚洲五月天丁香| 亚洲性夜色夜夜综合| 在线免费观看的www视频| 村上凉子中文字幕在线| 动漫黄色视频在线观看| 好看av亚洲va欧美ⅴa在| 91老司机精品| 亚洲avbb在线观看| 狂野欧美白嫩少妇大欣赏| 午夜福利视频1000在线观看| 97碰自拍视频| 国产精品一区二区三区四区免费观看 | 日本 av在线| 久久久久亚洲av毛片大全| 国产精品亚洲一级av第二区| 午夜日韩欧美国产| 国产高清videossex| 亚洲一区高清亚洲精品| 高清毛片免费观看视频网站| 日韩欧美免费精品| 十八禁网站免费在线| 国产午夜精品久久久久久| 国产亚洲av高清不卡| netflix在线观看网站| 一本精品99久久精品77| 国产亚洲欧美98| 国产精品一及| 日韩精品青青久久久久久| 小说图片视频综合网站| 欧美乱码精品一区二区三区| 国产在线精品亚洲第一网站| 国产精品 欧美亚洲| 亚洲免费av在线视频| 日韩高清综合在线| 18禁美女被吸乳视频| 久久精品影院6| 少妇被粗大的猛进出69影院| 又粗又爽又猛毛片免费看| 熟女电影av网| 不卡av一区二区三区| 日日爽夜夜爽网站| 校园春色视频在线观看| 女警被强在线播放| 一级片免费观看大全| 听说在线观看完整版免费高清| 国产私拍福利视频在线观看| 亚洲欧美日韩无卡精品| 老司机靠b影院| 亚洲专区国产一区二区| 国产亚洲精品一区二区www| 99热这里只有是精品50| 熟女少妇亚洲综合色aaa.| 亚洲免费av在线视频| 日韩有码中文字幕| 日韩大码丰满熟妇| 黄频高清免费视频| 琪琪午夜伦伦电影理论片6080| 亚洲男人天堂网一区| 久久这里只有精品中国| 色在线成人网| 最好的美女福利视频网| 国产三级中文精品| 成人av一区二区三区在线看| 亚洲成人免费电影在线观看| 99热6这里只有精品| 国产av一区二区精品久久| 国产三级中文精品| 12—13女人毛片做爰片一| 身体一侧抽搐| 一级黄色大片毛片| 天天一区二区日本电影三级| 日韩欧美在线乱码| 一区二区三区国产精品乱码| 精品国产乱子伦一区二区三区| svipshipincom国产片| 欧美成人免费av一区二区三区| 欧美3d第一页| 美女 人体艺术 gogo| 操出白浆在线播放| 亚洲精品久久国产高清桃花| 午夜亚洲福利在线播放| 久久中文字幕人妻熟女| 国产亚洲欧美在线一区二区| 男女床上黄色一级片免费看| 91麻豆精品激情在线观看国产| 国产一区二区在线观看日韩 | 久久精品国产亚洲av香蕉五月| 美女免费视频网站| 国产三级在线视频| 好男人电影高清在线观看| 麻豆久久精品国产亚洲av| 亚洲国产高清在线一区二区三| 欧美不卡视频在线免费观看 | 久久香蕉国产精品| 亚洲熟妇熟女久久| 久久久久久久久久黄片| 久久精品aⅴ一区二区三区四区| 免费看a级黄色片| 免费在线观看亚洲国产| 在线观看美女被高潮喷水网站 | 国产在线观看jvid| 国产成人系列免费观看| 亚洲美女黄片视频| 2021天堂中文幕一二区在线观| 国产男靠女视频免费网站| 国产男靠女视频免费网站| 久久九九热精品免费| 男插女下体视频免费在线播放| 午夜福利在线观看吧| 国模一区二区三区四区视频 | 欧美在线一区亚洲| 啪啪无遮挡十八禁网站| 免费人成视频x8x8入口观看| 久久亚洲真实| 欧美成人一区二区免费高清观看 | 欧美zozozo另类| 夜夜爽天天搞| a在线观看视频网站| 禁无遮挡网站| 成年人黄色毛片网站| 欧美乱妇无乱码| 久久香蕉激情| 午夜激情福利司机影院| 免费在线观看亚洲国产| 成年免费大片在线观看| 国产97色在线日韩免费| 国产97色在线日韩免费| 久久久久久亚洲精品国产蜜桃av| 国产69精品久久久久777片 | 狂野欧美白嫩少妇大欣赏| 久久这里只有精品19| 国产97色在线日韩免费| 欧美成人免费av一区二区三区| 亚洲国产中文字幕在线视频| 91大片在线观看| 亚洲天堂国产精品一区在线| 黄色毛片三级朝国网站| 成熟少妇高潮喷水视频| 亚洲国产中文字幕在线视频| 丰满人妻熟妇乱又伦精品不卡| 成人18禁高潮啪啪吃奶动态图| 精品国产乱子伦一区二区三区| 久久 成人 亚洲| 中文字幕高清在线视频| 无限看片的www在线观看| 婷婷亚洲欧美| 亚洲专区国产一区二区| 999久久久精品免费观看国产| 12—13女人毛片做爰片一| www.自偷自拍.com| 日本免费一区二区三区高清不卡| 国产精品影院久久| 亚洲国产精品999在线| 一a级毛片在线观看| 精品国产美女av久久久久小说| 日本精品一区二区三区蜜桃| 哪里可以看免费的av片| 欧美性猛交╳xxx乱大交人| 69av精品久久久久久| 日韩大码丰满熟妇| 麻豆国产av国片精品| 老司机午夜十八禁免费视频| 少妇的丰满在线观看| 高清在线国产一区| 精品久久久久久久久久免费视频| 国产亚洲精品综合一区在线观看 | 国产精品,欧美在线| 亚洲精品国产精品久久久不卡| 国产精品精品国产色婷婷| 精品免费久久久久久久清纯| 国产精品久久视频播放| 日韩有码中文字幕| 成人三级做爰电影| 精品一区二区三区av网在线观看| 国产亚洲精品av在线| 国产精品一区二区精品视频观看| 国产一区二区三区视频了| 午夜影院日韩av| 1024香蕉在线观看| 国产一区二区三区在线臀色熟女| 免费观看人在逋| 久久久久久人人人人人| 国产午夜精品论理片| 在线a可以看的网站| 久久精品人妻少妇| 国产一区二区三区视频了| 欧美三级亚洲精品| 国产三级中文精品| 亚洲欧美日韩无卡精品| 在线观看www视频免费| 国产亚洲av高清不卡| 欧美性长视频在线观看| 亚洲熟女毛片儿| 男人舔奶头视频| 国产精品久久久人人做人人爽| 哪里可以看免费的av片| 亚洲七黄色美女视频| 亚洲欧美日韩高清专用| www国产在线视频色| 最近最新中文字幕大全电影3| 真人做人爱边吃奶动态| 少妇熟女aⅴ在线视频| 日本黄色视频三级网站网址| 人人妻人人澡欧美一区二区| 一二三四社区在线视频社区8| 久久久久久久午夜电影| 在线观看美女被高潮喷水网站 | 国产真人三级小视频在线观看| 五月玫瑰六月丁香| 国产精品自产拍在线观看55亚洲| 精品久久久久久久毛片微露脸| 十八禁人妻一区二区| 久久久国产成人精品二区| 99热这里只有是精品50| 国产一区二区三区在线臀色熟女| 国产精品av视频在线免费观看| 国产午夜精品久久久久久| 久久精品国产亚洲av高清一级| 男插女下体视频免费在线播放| 桃色一区二区三区在线观看| 观看免费一级毛片| 国产免费男女视频| 欧美三级亚洲精品| 婷婷六月久久综合丁香| 国产高清视频在线观看网站| 亚洲九九香蕉| 久久久久九九精品影院| aaaaa片日本免费| 久久 成人 亚洲| 亚洲国产精品合色在线| 日本黄大片高清| 国产精品av视频在线免费观看| 国产欧美日韩一区二区精品| 久久久久久亚洲精品国产蜜桃av| 久久热在线av| 99国产精品一区二区三区| 最好的美女福利视频网| 丰满人妻一区二区三区视频av | 麻豆成人午夜福利视频| 亚洲精品国产一区二区精华液| 国产成人精品久久二区二区免费| 国产精品免费视频内射| 国产亚洲精品第一综合不卡| 国产成年人精品一区二区| 夜夜爽天天搞| 桃色一区二区三区在线观看| 精品一区二区三区av网在线观看| 色老头精品视频在线观看| 亚洲人与动物交配视频| 又大又爽又粗| 国产亚洲欧美98| 久久热在线av| 国产精品国产高清国产av| 男人舔女人的私密视频| 精品国产美女av久久久久小说| 免费看美女性在线毛片视频| 在线观看日韩欧美| www.熟女人妻精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩精品青青久久久久久| 18禁国产床啪视频网站| 欧美zozozo另类| 桃色一区二区三区在线观看| 国内精品久久久久精免费| 精品福利观看| 亚洲国产精品sss在线观看| 90打野战视频偷拍视频| a在线观看视频网站| 亚洲欧美精品综合久久99| 高潮久久久久久久久久久不卡| 久久久久国产精品人妻aⅴ院| 妹子高潮喷水视频| 免费在线观看日本一区| 亚洲自偷自拍图片 自拍| 亚洲一区中文字幕在线| 久久香蕉精品热| 久久欧美精品欧美久久欧美| 国产成年人精品一区二区| 亚洲av电影在线进入| 久9热在线精品视频| 久久精品国产99精品国产亚洲性色| 欧美国产日韩亚洲一区| 看片在线看免费视频| 老司机深夜福利视频在线观看| av免费在线观看网站| 亚洲成人久久性| 亚洲成人精品中文字幕电影| 美女黄网站色视频| 琪琪午夜伦伦电影理论片6080| 搡老妇女老女人老熟妇| 国产成人精品无人区| 亚洲一卡2卡3卡4卡5卡精品中文| 黄片小视频在线播放| 亚洲午夜理论影院| 搞女人的毛片| 老司机靠b影院| 好男人在线观看高清免费视频| 麻豆国产av国片精品| 最新美女视频免费是黄的| 国产精品1区2区在线观看.| 色哟哟哟哟哟哟| 久久精品aⅴ一区二区三区四区| 两个人视频免费观看高清| 国产欧美日韩一区二区精品| 久久久久久亚洲精品国产蜜桃av| 怎么达到女性高潮| 一本综合久久免费| 日韩精品青青久久久久久| 中文字幕人成人乱码亚洲影| 免费在线观看完整版高清| 国产伦人伦偷精品视频| 免费人成视频x8x8入口观看| 午夜精品在线福利| 久久精品国产99精品国产亚洲性色| 天堂影院成人在线观看| 国产成人av激情在线播放| 亚洲av中文字字幕乱码综合| 亚洲av五月六月丁香网| 一进一出抽搐gif免费好疼| 丁香欧美五月| 制服丝袜大香蕉在线| 中国美女看黄片| 成人国语在线视频| 午夜精品一区二区三区免费看| 黑人操中国人逼视频| 亚洲电影在线观看av| 国产一区二区在线观看日韩 | 久久久久久久午夜电影| 两个人视频免费观看高清| 精品熟女少妇八av免费久了| 精品国内亚洲2022精品成人| 在线a可以看的网站| 99热只有精品国产| 亚洲aⅴ乱码一区二区在线播放 | 精品国产超薄肉色丝袜足j| 免费在线观看完整版高清| 久久香蕉国产精品| 18禁黄网站禁片午夜丰满| 校园春色视频在线观看| 黑人操中国人逼视频| 又黄又爽又免费观看的视频| www日本黄色视频网| 淫秽高清视频在线观看| 中文字幕精品亚洲无线码一区| 亚洲人成电影免费在线| 美女扒开内裤让男人捅视频| 亚洲av中文字字幕乱码综合| 国产精品99久久99久久久不卡| 99国产精品99久久久久| 变态另类成人亚洲欧美熟女| 午夜福利高清视频| 亚洲,欧美精品.| 啪啪无遮挡十八禁网站| av片东京热男人的天堂| 中文字幕最新亚洲高清| 国产精品一区二区三区四区免费观看 | 久久久久亚洲av毛片大全| 久久久久久国产a免费观看| 亚洲精品中文字幕一二三四区| 女警被强在线播放| 欧美乱码精品一区二区三区| 欧美日韩黄片免| 亚洲18禁久久av| 欧美zozozo另类| 免费看日本二区| 麻豆一二三区av精品| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品第一综合不卡| 视频区欧美日本亚洲| 别揉我奶头~嗯~啊~动态视频| 最新美女视频免费是黄的| 丁香六月欧美| 久久人妻av系列| 窝窝影院91人妻| 亚洲国产精品合色在线| 国产一区二区三区视频了| 欧美日本亚洲视频在线播放| 欧美黑人欧美精品刺激| 亚洲一码二码三码区别大吗| 性欧美人与动物交配| 国产久久久一区二区三区| 国产高清videossex| 国产人伦9x9x在线观看| 中文字幕av在线有码专区| 亚洲一区中文字幕在线| 搡老岳熟女国产| 在线免费观看的www视频| 国产欧美日韩一区二区三| 听说在线观看完整版免费高清| 免费在线观看黄色视频的| 国产又黄又爽又无遮挡在线| 日韩精品中文字幕看吧| 老熟妇仑乱视频hdxx| 变态另类丝袜制服| 18美女黄网站色大片免费观看| 成人av在线播放网站| 正在播放国产对白刺激| 在线十欧美十亚洲十日本专区| 身体一侧抽搐| 看片在线看免费视频| 两人在一起打扑克的视频| 18禁黄网站禁片午夜丰满| 欧美三级亚洲精品| 久久精品国产综合久久久| 又黄又粗又硬又大视频| 亚洲人成网站在线播放欧美日韩| 18禁裸乳无遮挡免费网站照片| 久久久水蜜桃国产精品网| 一进一出好大好爽视频| 国产亚洲欧美在线一区二区| 国产精品久久久久久精品电影| 亚洲 国产 在线| 男女做爰动态图高潮gif福利片| 哪里可以看免费的av片| 在线观看免费午夜福利视频| 精品高清国产在线一区| 999久久久国产精品视频| 变态另类成人亚洲欧美熟女| 午夜两性在线视频| 国产免费av片在线观看野外av| 欧美日韩乱码在线| 91麻豆av在线| 一本大道久久a久久精品| 国产成人啪精品午夜网站| cao死你这个sao货| 欧美大码av| 亚洲 欧美 日韩 在线 免费| 18禁国产床啪视频网站| 欧美日韩一级在线毛片| 国产亚洲av高清不卡| 老汉色av国产亚洲站长工具| 狂野欧美白嫩少妇大欣赏| 非洲黑人性xxxx精品又粗又长| 久久伊人香网站| 黑人欧美特级aaaaaa片| 成人av一区二区三区在线看| 波多野结衣高清作品| 午夜精品久久久久久毛片777| 毛片女人毛片| 久久热在线av| 人妻夜夜爽99麻豆av| 国产精品香港三级国产av潘金莲| 亚洲精品色激情综合| 国产精品 国内视频| 99久久国产精品久久久| 久久99热这里只有精品18| 99国产极品粉嫩在线观看| 午夜a级毛片| 88av欧美| 久久精品国产亚洲av高清一级| 可以在线观看的亚洲视频| 精品久久久久久久人妻蜜臀av| 欧美丝袜亚洲另类 | 中文字幕久久专区| 麻豆一二三区av精品| 级片在线观看| 老司机午夜十八禁免费视频| 免费在线观看黄色视频的| 国产成人影院久久av| 老司机福利观看| 久久久久国产一级毛片高清牌| 久久午夜综合久久蜜桃| 波多野结衣高清无吗| 日韩欧美三级三区| 狂野欧美白嫩少妇大欣赏| 又黄又爽又免费观看的视频| 宅男免费午夜| 亚洲欧美日韩无卡精品| 岛国在线免费视频观看| 国产成人啪精品午夜网站| 久久午夜综合久久蜜桃| 国产男靠女视频免费网站| 91在线观看av| 丰满人妻一区二区三区视频av | 国产一级毛片七仙女欲春2| 国产探花在线观看一区二区| 欧美国产日韩亚洲一区| 欧美国产日韩亚洲一区| 国产视频内射| 午夜久久久久精精品| 午夜免费成人在线视频| 欧美在线黄色| 亚洲av熟女| 级片在线观看| 丁香六月欧美| 亚洲中文日韩欧美视频| 国产三级中文精品| av福利片在线观看| 每晚都被弄得嗷嗷叫到高潮| 一二三四社区在线视频社区8| 俺也久久电影网| 丁香欧美五月| 成在线人永久免费视频| 欧美中文日本在线观看视频| 亚洲精品在线美女| 嫩草影院精品99| 成人亚洲精品av一区二区| 亚洲精品在线美女| 国产成年人精品一区二区| 一本精品99久久精品77| 欧美+亚洲+日韩+国产| 色播亚洲综合网| 国产在线精品亚洲第一网站| √禁漫天堂资源中文www| xxx96com| 亚洲精品在线美女| 亚洲国产看品久久| av超薄肉色丝袜交足视频| 我要搜黄色片| 久久国产精品人妻蜜桃| 天天一区二区日本电影三级| 十八禁网站免费在线| 亚洲国产精品999在线| 久久精品91无色码中文字幕| 久久久久久久久免费视频了| 日本在线视频免费播放| 久久香蕉精品热| 听说在线观看完整版免费高清| 精品欧美一区二区三区在线| 国语自产精品视频在线第100页| 久久久国产精品麻豆| 国产真实乱freesex| 精品少妇一区二区三区视频日本电影| 精品国产乱码久久久久久男人| 俺也久久电影网| 欧美乱码精品一区二区三区| 国产av不卡久久| 国产成人啪精品午夜网站| 亚洲最大成人中文| 亚洲国产精品999在线| av中文乱码字幕在线| tocl精华| 国产精品电影一区二区三区| 国产精品久久久人人做人人爽| 国产精品日韩av在线免费观看| 亚洲avbb在线观看| 久久人人精品亚洲av| 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 国产成人av教育| 午夜激情福利司机影院| 亚洲专区国产一区二区| 天堂√8在线中文| 午夜成年电影在线免费观看| 成人av一区二区三区在线看| 他把我摸到了高潮在线观看| 亚洲av成人不卡在线观看播放网| 国产精品av久久久久免费| 亚洲精品久久国产高清桃花| 99久久99久久久精品蜜桃| av国产免费在线观看| 91国产中文字幕| 日本熟妇午夜| 国产野战对白在线观看| 中亚洲国语对白在线视频| www.熟女人妻精品国产| 国产不卡一卡二| 69av精品久久久久久| 无限看片的www在线观看| 99re在线观看精品视频| 久久精品成人免费网站| 毛片女人毛片| 免费在线观看亚洲国产| 国产午夜福利久久久久久| 香蕉av资源在线|