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

    基于校正WOFOST模型的棗樹(shù)生長(zhǎng)模擬與水分利用評(píng)價(jià)

    2021-01-17 23:20:30白鐵成王濤張楠楠
    智慧農(nóng)業(yè)(中英文) 2021年2期
    關(guān)鍵詞:水分利用效率棗樹(shù)

    白鐵成 王濤 張楠楠

    摘要:為實(shí)現(xiàn)定量化分析溫、光和水資源對(duì)果樹(shù)生長(zhǎng)的影響,本研究以成齡駿棗樹(shù)為研究對(duì)象,提出了基于校正WOFOST模型的棗樹(shù)生長(zhǎng)和水分運(yùn)移模擬方法。利用2016和2017年的田間試驗(yàn)觀測(cè)數(shù)據(jù),重點(diǎn)校正WOFOST模型的物候?qū)W發(fā)育、初始化、綠葉、CO同化、干物質(zhì)分配、呼吸作用和水分利用參數(shù)。在田間尺度,完成總地上生物量(TAGP)、葉面積指數(shù)(LAI)和土壤水分含量的動(dòng)態(tài)模擬和精度驗(yàn)證;在縣域尺度,使用55個(gè)果園的最大LAI、單產(chǎn)、實(shí)際蒸散量(ET)和水分利用效率(WUE)數(shù)據(jù)評(píng)價(jià)模型區(qū)域尺度的模擬性能。結(jié)果表明,在田間尺度,校正模型模擬不同灌溉梯度TAGP的決定系數(shù)R范圍為0.92~0.98,歸一化均方根誤差NRMSE為8.7%~20.5%;模擬LAI的R范圍為0.79~0.97,NRMSE為8.3%~21.1%;模擬土壤水分含量的決定系數(shù)R范圍為0.29~0.75,NRMSE為4.1%~6.1%。在縣域尺度,兩年模擬最大LAI與實(shí)測(cè)LAI的R分別為0.64和0.78,NRMSE分別為13.3%和10.7%;模擬單產(chǎn)的R分別為0.48和0.60,NRMSE分別是12.1%和11.9%;模擬ET均方根誤差分別為36.1mm(7.9%)和30.8mm(7.4%);模型也表現(xiàn)了較高的WUE模擬精度(10%<NRMSE<20%),均方根誤差EMSE值分別為0.23和0.28kg/m。WOFOST模型在田間和縣域尺度都取得了較高的棗樹(shù)生長(zhǎng)和水分運(yùn)移模擬精度,可為土壤、氣象、灌溉管理和棗樹(shù)生長(zhǎng)耦合影響的定量化分析提供新思路。

    關(guān)鍵詞:作物生長(zhǎng)模型;參數(shù)校正;WOFOST;棗樹(shù);水分利用效率

    中圖分類(lèi)號(hào):S126文獻(xiàn)標(biāo)志碼:A文章編號(hào):202103-SA008

    引用格式:白鐵成,王濤,張楠楠.基于校正WOFOST模型的棗樹(shù)生長(zhǎng)模擬與水分利用評(píng)價(jià)[J].智慧農(nóng)業(yè)(中英文),2021,3(2):55-67.

    BAI Tiecheng, WANG Tao, ZHANG Nannan. Dynamic simulation of jujube tree growth and water use evaluation based on the calibrated WOFOST model[J]. Smart Agriculture, 2021, 3(2): 55-67. (in Chinese with English abstract)

    1引言

    紅棗富含蛋白質(zhì)、糖類(lèi)、維生素和礦物質(zhì)等營(yíng)養(yǎng)成分,具有重要的食用和醫(yī)學(xué)價(jià)值。中國(guó)紅棗種植區(qū)域主要有新疆、山西、陜西、河北、河南和山東等地區(qū)。新疆南部地區(qū)光熱資源豐富、日照時(shí)間長(zhǎng)、晝夜溫差大,是中國(guó)最重要的優(yōu)質(zhì)紅棗生產(chǎn)基地,年產(chǎn)量372.8萬(wàn)噸,占全國(guó)總產(chǎn)量的50%[1],已成為南疆農(nóng)民增收和人口聚集的重要支柱產(chǎn)業(yè)。然而,新疆地區(qū)夏季溫度高、蒸發(fā)蒸騰量大、干旱少雨,水資源緊缺,如何定量評(píng)價(jià)蒸散和水分利用效率是干旱區(qū)農(nóng)業(yè)節(jié)水的關(guān)鍵技術(shù)問(wèn)題?,F(xiàn)有關(guān)于棗樹(shù)生長(zhǎng)水效應(yīng)和灌溉管理研究主要集中在灌水處理對(duì)產(chǎn)量和品質(zhì)的影響[2]、管理制度對(duì)土壤水分的影響[3]、灌溉模式優(yōu)化水肥管理[4]、灌溉和施肥對(duì)產(chǎn)量的耦合影響[5-7]等。然而,這些研究大多基于田間試驗(yàn)數(shù)據(jù),灌溉方案通常只適用于特定土壤和氣候環(huán)境,無(wú)法實(shí)現(xiàn)定量化分析灌溉方案對(duì)紅棗產(chǎn)量和品質(zhì)的影響。

    作物生長(zhǎng)模型是基于生理和生態(tài)學(xué)原理,并使用數(shù)學(xué)模型量化受氣候、土壤和管理?xiàng)l件影響的作物生長(zhǎng)發(fā)育過(guò)程。在過(guò)去的幾十年中,已經(jīng)開(kāi)發(fā)的幾個(gè)比較成熟的作物模型包括WOFOST(WOrldFOod Studies)、DSSAT(Decision Support System for Agrotechnology Transfer)和EPIC(Environmental. Policy Integrated Climate)等[8]。這些模型各具特點(diǎn)和優(yōu)勢(shì),但通常都包含了物候?qū)W、CO同化(光合生產(chǎn))、呼吸作用、同化物分配和土壤水循環(huán)等模塊,被廣泛應(yīng)用于農(nóng)業(yè)生產(chǎn)管理以及與社會(huì)相關(guān)的應(yīng)用中,如氣候變化響應(yīng)[9],作物產(chǎn)量預(yù)測(cè)[10]以及環(huán)境、水分和作物生長(zhǎng)的耦合作用[11]等。已有的作物生長(zhǎng)模型主要聚焦一年生的大田作物,關(guān)于果樹(shù)作物模型研究方面,F(xiàn)ishman和Génard[12]基于水和干物質(zhì)運(yùn)輸機(jī)制開(kāi)發(fā)了桃樹(shù)果實(shí)生長(zhǎng)模型,計(jì)算干旱和新增質(zhì)量變化所表示的季節(jié)性和晝夜水果生長(zhǎng)情況,以分析作物負(fù)荷下的水分脅迫條件。Lopez等[13]基于形態(tài)學(xué)開(kāi)發(fā)了桃樹(shù)生長(zhǎng)的計(jì)算機(jī)模擬模型,實(shí)現(xiàn)桃樹(shù)生長(zhǎng)的三維描述,可用于樹(shù)木修剪管理。Lescourret等[14]用全局模型解釋了雌性和雄性藤的開(kāi)花、授粉和果實(shí)生長(zhǎng)之間的關(guān)系,并集成了生產(chǎn)管理輸入,仿真了果園管理對(duì)琳猴桃果實(shí)數(shù)量和尺寸的影響,該模型對(duì)氣候和技術(shù)操作的變化都很敏感。Green等[15]開(kāi)發(fā)了果樹(shù)冠層蒸騰和光能截獲模型,用于證明蒸騰作用如何響應(yīng)氣象環(huán)境的變化,敏感性分析表明,光攔截受葉面積和葉光學(xué)特性變化的影響最大,而蒸騰作用受葉面積和葉電導(dǎo)變化的影響最大。Costes等[16]使用混合隨機(jī)和生物力學(xué)模型模擬了果樹(shù)發(fā)育。

    然而,已報(bào)道的果樹(shù)生長(zhǎng)模型缺少生理生化過(guò)程的機(jī)理性描述。WOFOST是一種機(jī)理模型,可以根據(jù)物候?qū)W、CO同化、蒸騰作用、呼吸作用以及這些過(guò)程如何受環(huán)境條件影響來(lái)解釋作物生長(zhǎng)[8]。該模型的過(guò)程描述是通用的,通過(guò)參數(shù)校正可用于不同作物的生長(zhǎng)模擬[17]和水分限制產(chǎn)量評(píng)估[18],不僅被廣泛用于定量分析一年生大田作物的生長(zhǎng)過(guò)程[8],而且具有棗樹(shù)生長(zhǎng)模擬和產(chǎn)量評(píng)估的潛力[19]。基于WOFOST實(shí)現(xiàn)多年生果樹(shù)作物的水分限制生長(zhǎng)模擬,進(jìn)而實(shí)現(xiàn)土壤、氣象、灌溉管理和果樹(shù)作物生長(zhǎng)耦合影響的定量化分析,有望使果樹(shù)灌溉管理從試驗(yàn)統(tǒng)計(jì)分析邁入定量、機(jī)理性描述階段。因此,本研究以矮化密植駿棗為研究對(duì)象,在WOFOST模型作物參數(shù)校正基礎(chǔ)上,實(shí)現(xiàn)棗樹(shù)生長(zhǎng)和水分運(yùn)移模擬過(guò)程。并使用大田試驗(yàn)測(cè)量的總地上生物量(Total. Above-Ground Production,TAGP)、葉面積指數(shù)(Leaf Area Index, LAI)和土壤水分含量驗(yàn)證田間尺度的模擬精度;使用縣域尺度測(cè)量的55個(gè)果園的最大葉面積指數(shù)、單產(chǎn)、實(shí)際蒸散量和水分利用效率評(píng)價(jià)校正WOFOST模型在縣域尺度不同棗園的模擬性能。

    2材料與方法

    2.1研究區(qū)概況

    研究區(qū)位于新疆維吾爾自治區(qū)阿拉爾市(80°30′E—81°58′E,40°22′N(xiāo)—40°57′N(xiāo)),紅棗種植面積約4.5萬(wàn)公頃,約占全國(guó)紅棗產(chǎn)量的1/8。屬于暖溫帶極端大陸性干旱荒漠氣候,2016年和2017年平均降雨量分別為106和98mm,主要分布在夏季,太陽(yáng)輻射年均133.7~146.3 kcal/cm,年均日照2556.3~2991.8h,年平均氣溫在10.8~12.5℃之間,最大日溫差為20℃。豐富的光熱資源為優(yōu)質(zhì)紅棗種植提供了良好的自然條件。

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

    在田間試驗(yàn)區(qū)域(圖1藍(lán)色方塊定位區(qū)域)選擇成齡棗園在2016和2017年開(kāi)展灌溉試驗(yàn),采用滴灌方式在生育期設(shè)計(jì)5個(gè)灌溉梯度(不包括冬灌和春灌):灌溉定額CK(375 mm,D1)、90% CK(338 mm,D2)、80% CK(300 mm,D3)、70% CK (263 mm,D4)和60% CK(225 mm,D5)。灌溉定額參考了當(dāng)?shù)亟?jīng)驗(yàn)灌溉量,生長(zhǎng)季灌溉10次,各小區(qū)施肥量和其他農(nóng)藝措施相同。

    2.3測(cè)試項(xiàng)目及方法

    在田間試驗(yàn)區(qū)域測(cè)量的主要數(shù)據(jù)包括以下內(nèi)容。

    ①物候?qū)W發(fā)育時(shí)間:記錄萌芽,開(kāi)花(坐果)和成熟日期。

    ②干重測(cè)量:每10d取樣一次,生長(zhǎng)季取樣10次,每次選均勻的棗樹(shù)3棵,取全樹(shù)葉片、新枝和果實(shí)帶回實(shí)驗(yàn)室在80℃下烘干至恒重后測(cè)量各器官干重。

    ③光合有效輻射和葉面積指數(shù):每10d使用SunScan冠層分析儀(Delta公司,英國(guó))測(cè)定不同冠層空間的有效葉面積指數(shù)。

    ④光合作用:每10d測(cè)量一次,使用LI6400XT便攜式光合作用測(cè)試儀分層測(cè)試凈光合速率、氣孔導(dǎo)度、胞間CO濃度和蒸騰速率。

    ⑤土壤水分:灌溉前后在0~100cm土壤深度每20cm用根鉆取土樣,測(cè)定土壤容積含水量,每年測(cè)量14次。

    ⑥其他土壤參數(shù):土壤田持、容積密度、飽和土壤含水率、土壤水響應(yīng)曲線(xiàn)和滲透系數(shù)等直接取樣帶回實(shí)驗(yàn)室進(jìn)行測(cè)量,同步記錄灌溉日期和灌溉量。

    ⑦氣象數(shù)據(jù):由塔里木大學(xué)紅棗研究基地架設(shè)的小型氣象站提供,包括15min間隔的溫度、濕度、輻射、風(fēng)速、降雨和氣壓數(shù)據(jù)等。

    在縣域尺度觀測(cè)55個(gè)棗園(圖1紅點(diǎn)坐標(biāo))的初始棗樹(shù)干重、最大發(fā)育階段的葉面積指數(shù)、最終產(chǎn)量和實(shí)際蒸散。葉面積指數(shù)觀測(cè)時(shí)間分別為2016年7月24日和2017年7月27日,每年11月份觀測(cè)棗園產(chǎn)量。實(shí)際蒸散量ET。計(jì)算如公式(1)所示[20]。

    ET=P+I-D-?SM(1)

    其中,P、I、D和?SM分別為有效降雨量、灌溉量、深層土壤滲透、生長(zhǎng)季開(kāi)始和結(jié)束時(shí)土壤平均水分含量差值。新疆地區(qū)駿棗樹(shù)根系存在于0~100cm[21](主要分布在垂直60 cm內(nèi)土層),研究中測(cè)量100cm土壤層滲透量作為深層土壤滲透量。研究區(qū)地下水較深,地表徑流量極少,所以研究中未考慮地下水和地表徑流對(duì)ET的影響。

    水分利用效率(Water Use Efficiency,WUE)是節(jié)水農(nóng)業(yè)的重要指標(biāo),尤其在干旱和半干旱地區(qū),本研究中也將使用WUE進(jìn)一步驗(yàn)證校正模型的區(qū)域尺度模擬精度。WUE等于實(shí)測(cè)產(chǎn)量與ET的比值。

    2.4WOFOST模型參數(shù)校正

    WOFOST模型參數(shù)包括氣象、土壤和作物參數(shù),由于研究區(qū)域?yàn)榭h域尺度,氣象數(shù)據(jù)空間變異性較小,因此直接使用氣象站觀測(cè)數(shù)據(jù),土壤基本屬性參數(shù)直接使用測(cè)量值。本研究主要校正作物參數(shù),對(duì)于敏感性較低的參數(shù),如發(fā)芽的最低溫度TBASEM和葉片老化的最低溫度TBASE通過(guò)查閱文獻(xiàn)[22]確定,敏感性較高的參數(shù),如初始生物量TDWI和葉片死亡率SPAN通過(guò)2016年D1灌溉處理的田間觀測(cè)數(shù)據(jù)進(jìn)行校正。參考Wit等[23]的關(guān)于WOFOST模型的校準(zhǔn)方法,結(jié)合棗樹(shù)生理特點(diǎn),WOFOST模型主要參數(shù)的校正過(guò)程如下。

    (1)物候?qū)W參數(shù)校正。模型校準(zhǔn)從物候?qū)W參數(shù)開(kāi)始,通過(guò)有效積溫法確定[23]。棗樹(shù)開(kāi)始發(fā)育的最低閾值溫度(Lower Threshold Temperature for Emergence,TBASEM)被設(shè)置為10℃[22]。WOFOST給出的萌芽時(shí)最高有效溫度(Max Effective Temperature for Emergence,TEFFMX)范圍為18~32℃。由于2016和2017年棗樹(shù)萌芽期最高溫度低于30℃,萌芽期并無(wú)高溫脅迫,因此,將最高有效溫度設(shè)置為30℃。從開(kāi)始發(fā)育至萌芽的有效積溫和(Temperature Sum from Sowing to Emergence,TSUMEM)、萌芽至開(kāi)花的有效積溫和(Temperature Sum from Emergence to Anthesis,TSUM1)以及開(kāi)花至成熟的有效積溫和(Temperature Sum from Anthesis to Maturity,TSUM2)3個(gè)指標(biāo)通過(guò)觀測(cè)每日平均溫度和萌芽、開(kāi)花以及成熟日期進(jìn)行校準(zhǔn)。有效積溫的每日增加量(Daily Increase in Temperature Sum,DTSMTB)根據(jù)最高適宜溫度(35.5℃)[24]和最低發(fā)育溫度(10℃)[22]計(jì)算。

    (2)初始干重校正。棗樹(shù)與一年生作物不同,如果考慮上一年的莖,可能會(huì)導(dǎo)致初始作物干重值過(guò)高。本研究中初始棗樹(shù)干重被定義為初始新器官(初始芽和根)的重量,通過(guò)模擬和實(shí)測(cè)的總地上生物量和葉面積指數(shù)值進(jìn)行初始干重TDWI的校正。縣域尺度的55個(gè)果園的初始棗樹(shù)干重值根據(jù)實(shí)測(cè)的芽重和芽的分配系數(shù)計(jì)算獲得。其中芽的分配系數(shù)來(lái)自于田間試驗(yàn)的測(cè)量數(shù)據(jù)。

    (3)綠葉參數(shù)校正。在所有試驗(yàn)小區(qū)中,萌芽時(shí)(定義為芽上第五片葉子展開(kāi)時(shí))的葉面積指數(shù)(Leaf Area Index at Emergence,LAIEM)的測(cè)量值均小于0.004。另外,仿真結(jié)果表明,當(dāng)LAIEM從0.0007(WOFOST給出的最小值)增加到0.01時(shí),總地上生物量、活葉干重、活莖干重、果實(shí)干重和LAI的模擬結(jié)果差異極小,因此,LAIEM參數(shù)可以設(shè)置為0.0007至0.01之間的任何值。測(cè)量的LAI最大相對(duì)增加量(Maximum Relative Increase in LAI,RGRLAI)均小于0.05,而且當(dāng)RGRLAI從0.05增加到默認(rèn)最大值0.5時(shí),總地上生物量和葉面積指數(shù)模擬結(jié)果幾乎不變,因此RGRLAI被設(shè)置為最小值0.05。葉面積與干重比值(Specific Leaf Area, SLATB)根據(jù)測(cè)量的總地上生物量和葉面積指數(shù)進(jìn)行校準(zhǔn)。葉片老化速率(Life Span of Leaves Growing at 35 Celsius,SPAN)根據(jù)生長(zhǎng)季后期田間測(cè)量的LAI進(jìn)行了校正。

    (4)CO同化參數(shù)校準(zhǔn)。分別在棗樹(shù)最佳發(fā)育下限溫度19.5℃和上限溫度35.5℃擬合光響應(yīng)曲線(xiàn),擬合獲得最大的CO同化速率(Maximum Leaf CO Assimilation Rate,AMAXTB)和光能利用效率(Light-use Efficiency for Single Leaf,EFFTB)初始值,再通過(guò)實(shí)測(cè)總地上生物量和葉面積指數(shù),對(duì)消光系數(shù)(Extinction Coefficient for Diffuse Visible Light,KDIFTB)、AMAXTB和EFFTB進(jìn)行了校準(zhǔn)。

    (5)干物質(zhì)分配參數(shù)。FSTB(分配到莖的系數(shù))、FOTB(分配到果實(shí)的系數(shù))和FLTB(分配到葉的系數(shù))為發(fā)育階段的函數(shù),使用測(cè)量和模擬的地上各器官干物質(zhì)重量進(jìn)行校準(zhǔn)。

    (6)水分利用參數(shù)校正。使用測(cè)量的土壤水分含量校正蒸騰速率修正系數(shù)(Correction Factor Transpiration Rate,CFET)、水分脅迫敏感性校正系數(shù)(Correction Factor for Crop Water Stress Sensitivity,DEPNR)、初始根深(Initial. Rooting Depth,RDI)、每日根深最大增量(Maximum Daily Increase in Rooting Depth,RR1)和最大根深(Maximum Rooting Depth, RD- MCR)。

    2.5模型驗(yàn)證

    模型驗(yàn)證包括兩部分。在田間試驗(yàn)區(qū)域,使用2017年D1~D5處理測(cè)量的時(shí)間序列的總地上生物量、LAI和土壤水分含量驗(yàn)證棗樹(shù)生長(zhǎng)和土壤水分運(yùn)移的模擬精度??h域測(cè)量的55個(gè)果園最大發(fā)育階段的LAI、單產(chǎn)、實(shí)際蒸散量和水分利用效率用于驗(yàn)證模型區(qū)域模擬性能。使用決定系數(shù)(Coefficient of Determination,R)、均方根誤差(Root Mean Square Error, RMSE)和歸一化均方根誤差(Normalized Root Mean Square Error,NRMSE)評(píng)估田間試驗(yàn)處理的模擬性能,使用R,RMSE,NRMSE,相對(duì)偏差(Relative Bias Error,RBE)頻率分布和變異系數(shù)(Coefficient of Variation,CV)評(píng)價(jià)縣域尺度棗園的最大LAI和單產(chǎn)模擬性能。使用RBE頻率分布、RMSE和NRMSE評(píng)價(jià)實(shí)際蒸散和水分利用效率的模擬精度。其中,R表達(dá)實(shí)測(cè)值與模擬值之間的一致性,RMSE和NRMSE表達(dá)測(cè)量值和模擬值相對(duì)誤差和絕對(duì)誤差,NRMSE≤10%表示極高精度,10%<NRMSE≤20%表示高精度,20%<NRMSE≤30%表示中等精度,NRMSE>30%表示低精度[10]。CV反映了模擬結(jié)果空間變異性,CV<10%表示弱變異性,10%≤CV≤100%表示中等變異性,CV>100%表示強(qiáng)變異性[10]。RBE頻率分布表示樣本被高估和低估的比例。R、RMSE、NRMSE、RBE和CV計(jì)算見(jiàn)公式(2)~(6)。

    (2)

    (3)

    (4)

    (5)

    (6)

    其中,為第i個(gè)樣本模擬值;y表示第i個(gè)樣本測(cè)量值;表示測(cè)量樣本平均值;表示模擬樣本平均值;SD表示模擬結(jié)果標(biāo)準(zhǔn)差;n為樣本總數(shù),個(gè)。指標(biāo):總地上生物量,kg/ha;葉面積指數(shù),m/m;土壤體積含水量,cm/cm;單產(chǎn),kg/ha;最大葉面積指數(shù),m/m;實(shí)際蒸散量,mm;水分利用效率,kg/m。

    3結(jié)果與分析

    3.1模型參數(shù)校正

    根據(jù)上述模型參數(shù)校正方法,校正的棗樹(shù)作物參數(shù)結(jié)果如表1所示。參數(shù)分別來(lái)自測(cè)量值(m)、測(cè)量值基礎(chǔ)上的校正值(m-c)、估測(cè)值(e)、校正值(c)和文獻(xiàn)參考值。

    3.2田間驗(yàn)證

    在田間試驗(yàn)區(qū)域,通過(guò)2016年測(cè)量D1處理的田間數(shù)據(jù)校正模型參數(shù),2017年生長(zhǎng)季內(nèi)D1~D5處理的時(shí)間序列的總地上生物量、葉面積指數(shù)和土壤水分含量被用于驗(yàn)證模擬精度。

    3.2.1TAGP和LAI驗(yàn)證

    TAGP和LAI模擬精度是評(píng)價(jià)棗樹(shù)生長(zhǎng)模擬過(guò)程的重要指標(biāo)。表2顯示了不同處理的TAGP 和LAI模擬結(jié)果。校正模型模擬的D1~D3處理的TAGP與實(shí)測(cè)值一致性較好(R大于0.94),模擬精度較高(8.7%≤NRMSE≤16.9%),模擬D4和D5處理的總地上生物量的性能低于D1、D2和D3處理,R為0.92,NRMSE約等于20%。模擬D1~D4處理LAI的一致性和精度較高,R大于0.9,NRMSE小于15%,D5處理的模擬結(jié)果略差,R等于0.79,NRMSE大于20%。表明隨著灌水量的降低,水分脅迫強(qiáng)度將增加,校正模型對(duì)水分脅迫響應(yīng)的模擬能力可能下降??傊?,在70%~100%對(duì)照灌溉量時(shí),校正模型都表現(xiàn)了較好的總地上生物量和葉面積指數(shù)模擬性能。

    3.2.2土壤水分含量驗(yàn)證

    模擬和實(shí)測(cè)D1~D5處理的土壤水分含量的散點(diǎn)圖見(jiàn)圖2,測(cè)量值與模擬值總體趨勢(shì)保持一致,模擬與實(shí)測(cè)土壤水分含量的R在0.29~0.75之間,RMSE在0.010~0.013之間,NRMSE在4.1%~6.1%之間。D4和D5處理的模擬和實(shí)測(cè)土壤水分含量一致性略低于D1~D3處理,長(zhǎng)低于0.5。所有處理的土壤水分含量模擬精度較高,NRMSE均小于6.1%,絕對(duì)誤差較小。

    綜合總地上生物量、葉面積指數(shù)和土壤水分含量模擬的評(píng)價(jià)指標(biāo)R和NRMSE,當(dāng)灌溉量不低于經(jīng)驗(yàn)灌溉量的80%時(shí),模型具有較高的一致性和模擬精度。雖然在70%和60%的經(jīng)驗(yàn)灌溉量時(shí),模型模擬棗樹(shù)生長(zhǎng)過(guò)程總地上生物量和葉面積指數(shù)的性能略低,但最大絕對(duì)誤差NRMSE≤21.1%,仍可以實(shí)現(xiàn)中等精度的模擬。在所有灌溉處理中,校正模型都表現(xiàn)了極高的土壤水分含量模擬精度(NRMSE<10%)。

    3.3縣域尺度驗(yàn)證

    3.3.1模擬最大LAI驗(yàn)證

    分別輸入實(shí)測(cè)55個(gè)棗園的初始棗樹(shù)干重參數(shù)運(yùn)行模型,55個(gè)棗園的最大LAI模擬結(jié)果散點(diǎn)圖見(jiàn)圖3,2016和2017年模擬最大LAI與實(shí)測(cè)LA1的R分別為0.64和0.78,一致性較好,NRMSE分別為13.3%和10.7%,模擬精度較高。模擬LAI的變異系數(shù)CV值分別等于21.6%和21.2%,表達(dá)了顯著的空間變異性,也解釋了初始棗樹(shù)干重區(qū)域空間分布變異。

    3.3.2模擬單產(chǎn)驗(yàn)證

    2016和2017年模擬與實(shí)測(cè)單產(chǎn)的一致性指標(biāo)R分別為0.48和0.60(圖4),單產(chǎn)預(yù)測(cè)精度RMSE分別為0.83和0.94 kg/ha,NRMSE分別是12.1%和11.9%,具有較高的單產(chǎn)預(yù)測(cè)精度(NRMSE<20%)。圖4也表明高產(chǎn)區(qū)域明顯被低估,影響模擬產(chǎn)量的因素較多,包括CO同化參數(shù)、干物質(zhì)分配、葉面積與干重比值和葉片老化速率。本研究結(jié)果可能的原因是葉片老化速率參數(shù)設(shè)置值過(guò)低,田間試驗(yàn)測(cè)量的葉片老化速率值在40~60之間,縣域尺度部分果園的葉片老化速率值可能高于設(shè)置的固定值50,甚至可能高于60,所以導(dǎo)致高產(chǎn)區(qū)域被低估。

    單產(chǎn)模擬的相對(duì)誤差頻率分布如圖5所示,大部分樣本頻率分布在y=0附近,平均相對(duì)誤差分別為-3.23%和-2.94%,相對(duì)誤差絕對(duì)值的平均值分別為9.72%和9.68%,單產(chǎn)預(yù)測(cè)精度較高。CV分別是12.5%和11.9%,具有一定產(chǎn)量空間分布差異,不同樹(shù)齡和種植密度導(dǎo)致的初始棗樹(shù)干重的空間異質(zhì)性可能導(dǎo)致產(chǎn)量分布差異。然而,2016和2017年分別有61.8%和63.6%的樣本被低估,葉片老化速率空間變異性可能是導(dǎo)致模擬精度下降的原因。

    3.3.3模擬實(shí)際蒸散量和水分利用效率驗(yàn)證

    模擬的55個(gè)觀測(cè)果園的實(shí)際蒸散量和水分利用效率與測(cè)量值的相對(duì)百分比誤差頻率分布如圖6所示。所有果園的實(shí)際ET模擬誤差小于20%,其中49個(gè)果園的模擬誤差小于10%。2016和2017年模擬ET均方根誤差分別為36.1mm(7.9%)和30.8mm(7.4%)(表3),表現(xiàn)了極高的模擬精度(NRMSE<10%)。2016年模擬ET平均相對(duì)誤差略高于測(cè)量值,2017年低于測(cè)量值,平均相對(duì)誤差的絕對(duì)值分別為6.44%和5.98%,校正模型的整體ET模擬精度較高。模擬的WUE相對(duì)誤差略低于ET,在2016和2017年兩個(gè)生長(zhǎng)季,模擬WUE的相對(duì)百分比誤差小于10%的樣本分別占50.9%和47.3%,小于20%的分另占85.5%和80%,最大值分別為37.4%和38.3%。2016年模擬WUE的平均相對(duì)誤差低于測(cè)量值,2017年略高于測(cè)量值,平均相對(duì)誤差的絕對(duì)值分別為11.39%和12.73%。校正的模型表現(xiàn)了較高的WUE模擬精度(10%<NRMSE<20%),RMSE值分別為0.23和0.28 kg/m。

    4討論

    4.1田間尺度模擬精度分析

    田間尺度驗(yàn)證結(jié)果表明,基于時(shí)間序列的總地上生物量和葉面積指數(shù)的模擬結(jié)果均達(dá)到中等精度以上,土壤水分含量達(dá)到極高的模擬精度。然而,D1~D3的模擬性能略高于D4和D5處理,說(shuō)明校正的模型隨著灌溉量的降低、水分脅迫的增強(qiáng),模擬精度可能呈現(xiàn)下降的趨勢(shì)。原因可能是研究中僅僅使用2016年D1處理數(shù)據(jù)校正模型參數(shù),D1處理接近潛在水分供應(yīng),2017年D1~D3處理灌溉量大于80%CK,與2016年D1處理的灌溉量差異不大,棗樹(shù)面臨水分脅迫較小,對(duì)棗樹(shù)產(chǎn)量影響也較小。D4和D5處理可能導(dǎo)致更嚴(yán)重的干旱脅迫,實(shí)際生產(chǎn)中對(duì)棗樹(shù)生長(zhǎng)的影響可能高于模擬過(guò)程的干旱脅迫修正結(jié)果,進(jìn)而導(dǎo)致總地上生物量和LAI的模擬誤差。LAI的模擬偏差也將導(dǎo)致D4和D5處理的棗樹(shù)蒸騰和土壤蒸發(fā)的計(jì)算誤差,進(jìn)而導(dǎo)致土壤水分含量模擬精度略低于D1~D3處理。

    另外,在大田試驗(yàn)過(guò)程中,棗樹(shù)生長(zhǎng)通常會(huì)收到施肥、病蟲(chóng)害、風(fēng)、熱脅迫等因素影響。在本研究中,模型無(wú)法響應(yīng)這些因素對(duì)棗樹(shù)生長(zhǎng)的影響,在已校正的棗樹(shù)生長(zhǎng)模型中耦合經(jīng)典的氮素運(yùn)移模塊或病蟲(chóng)害響應(yīng)模型有望進(jìn)一步提高模擬精度和機(jī)理性。

    4.2模型區(qū)域尺度應(yīng)用的不確定性

    在運(yùn)行WOFOST作物生長(zhǎng)模型進(jìn)行生長(zhǎng)模擬和產(chǎn)量評(píng)估時(shí),輸入?yún)?shù)、氣象驅(qū)動(dòng)數(shù)據(jù)的不確定性以及模型模擬過(guò)程的簡(jiǎn)化都會(huì)影響單產(chǎn)估算的準(zhǔn)確性[25]。本研究預(yù)期模擬誤差主要來(lái)自輸入?yún)?shù)的不確定性,初始棗樹(shù)干重參數(shù)強(qiáng)烈影響農(nóng)作物的初始生長(zhǎng)速率,并顯示出高度的不確定性,這可能會(huì)影響初始LAI和最大LAI的生長(zhǎng)速率[26],進(jìn)而影響生長(zhǎng)模擬和單產(chǎn)評(píng)估精度。通過(guò)實(shí)測(cè)55個(gè)果園的初始棗樹(shù)干重值驅(qū)動(dòng)模型在縣域尺度模擬最大LAI和單產(chǎn)精度較高,說(shuō)明縣域尺度的棗園初始棗樹(shù)干重空間異質(zhì)性較大。在已報(bào)道的研究中,當(dāng)使用相同樹(shù)齡果園的初始棗樹(shù)干重平均值驅(qū)動(dòng)模型時(shí),181個(gè)區(qū)域棗園的單產(chǎn)驗(yàn)證NRMSE為16.3%(2016年)和17.2%(2017年)[19],明顯高于本研究中的12.1%和11.9%的驗(yàn)證誤差,這也說(shuō)明初始棗樹(shù)干重為棗樹(shù)生長(zhǎng)模擬過(guò)程的高敏感參數(shù)。另外,盡管決定葉片衰老速率和時(shí)間的葉片老化速率參數(shù)主要取決于作物遺傳學(xué)特性,但該參數(shù)也受干旱、營(yíng)養(yǎng)脅迫、病蟲(chóng)害以及作物管理等因素的影響[27]。但WOFOST模型無(wú)法模擬這些因素對(duì)葉片老化速率的影響[28]。研究中為葉片老化速率設(shè)置一個(gè)固定值(50)無(wú)法表達(dá)所有55個(gè)樣本的實(shí)際葉片老化或損傷情況,模擬最大LAI的誤差最可能來(lái)自于葉片老化速率的空間分布誤差。在以往研究中,通過(guò)遙感反演LAI同化潛在模式的WOFOST模型可以降低初始棗樹(shù)干重和葉片老化速率的不確定性,提高區(qū)域尺度的棗樹(shù)單產(chǎn)模擬精度,集合卡爾曼濾波EnKF和SUBPLEX算法同化遙感反演LAI后模擬單產(chǎn)的分別為9.2%和8.3%[29]。表明初始棗樹(shù)干重和葉片老化速率為區(qū)域尺度棗樹(shù)生長(zhǎng)模擬的不確定性參數(shù),遙感同化方法可以被考慮降低這一不確定性。

    另外,區(qū)域尺度的土壤屬性的空間異質(zhì)性也會(huì)導(dǎo)致最大LAI和水分運(yùn)移過(guò)程的模擬精度。同化遙感反演的LAI和ET狀態(tài)變量可以降低土壤水分含量模擬的不確定性,進(jìn)而提高產(chǎn)量評(píng)估精度。因此,縣域尺度單產(chǎn)和蒸散量模擬誤差可能來(lái)自灌溉、施肥和修剪管理導(dǎo)致的初始棗樹(shù)干重、葉片老化速率和土壤屬性參數(shù)的不確定性。在今后的研究中,可通過(guò)遙感同化土壤水分、蒸發(fā)蒸騰和葉面積指數(shù)等狀態(tài)變量對(duì)區(qū)域尺度的初始棗樹(shù)干重、葉片老化速率和土壤屬性參數(shù)進(jìn)行校正,提高區(qū)域尺度模擬精度,為在區(qū)域尺度分析水分脅迫對(duì)棗樹(shù)生長(zhǎng)和產(chǎn)量的影響提供機(jī)理性模擬模型。

    5結(jié)論

    本研究以WOFOST模型為基礎(chǔ),以縣域尺度的棗樹(shù)生長(zhǎng)模擬和水分利用評(píng)價(jià)為目標(biāo),使用田間實(shí)測(cè)棗樹(shù)生理生化指標(biāo)和土壤理化數(shù)據(jù)校正模型參數(shù)。在田間試驗(yàn)果園獲得了較高的TAGP 和LAI模擬精度(8.3%≤NRMSE≤21.1%),以及極高的土壤水分含量模擬精度(4.1%≤NRMSE≤6.1%)。使用校正的模型和55個(gè)果園測(cè)量的初始棗樹(shù)干重驅(qū)動(dòng)模型進(jìn)行了縣域尺度的棗樹(shù)生長(zhǎng)模擬和水分使用評(píng)價(jià)。結(jié)果表明,2016和2017年模擬最大LAI和單產(chǎn)的NRMSE均小于20%,實(shí)現(xiàn)了較高的模擬精度。校正模型也表現(xiàn)了極高的ET模擬精度(NRME<10%)和較高的WUE模擬精度(NRMSE<20%)??傊谔镩g和區(qū)域尺度的棗園,WOFOST模型都表現(xiàn)了較為理想的棗樹(shù)生長(zhǎng)模擬精度和水分利用評(píng)價(jià)性能,可為分析氣象、土壤、灌溉管理和棗樹(shù)生長(zhǎng)耦合影響提供一種新思路。

    參考文獻(xiàn):

    [1]國(guó)家統(tǒng)計(jì)局.數(shù)據(jù)查詢(xún)[EB/OL].[2021-06-07]. http://www.stats.gov.cn/.

    [2]胡家?guī)?,鄭旭榮,王振華,等.不同滴灌灌水處理對(duì)南疆沙區(qū)紅棗產(chǎn)量及品質(zhì)的影響[J].節(jié)水灌溉,2016(9): 40-43.

    HU J, ZHENG X, WANG Z, et al. Effects of different drip irrigation treatment on yield and quality in southern Xinjiang sand jujube[J]. Water Saving Irrigation,2016 (9):40-43.

    [3] WANG J, JUN H, WU P, et al. Effects of soil managements on surface runoff and soil water content in jujube orchard under simulated rainfalls[J]. Catena, 2015, 135: 193-201.

    [4] WANG Z, BIAN Q, ZHANG J, et al. Optimized water and fertilizer management of mature jujube in Xinjiang arid area using drip irrigation[J]. Water (Switzerland), 2018, 10(10): ID 1467.

    [5] DAI Z, FEI L, HUANG D, et al. Coupling effects of irrigation and nitrogen levels on yield, water and nitrogen use efficiency of surge-root irrigated jujube in a semiarid region[J]. Agricultural. Water Management, 2019, 213: 146-154.

    [6]王振華,扁青永,李文昊,等.南疆沙區(qū)成齡紅棗水肥一體化滴灌的水肥適宜用量[J].農(nóng)業(yè)工程學(xué)報(bào),2018,34(11): 96-104.

    WANG Z, BIAN Q, LI W. Suitable water and fertilizer amount for mature jujube with drip-irrigation under fertigation in southern Xinjiang sandy area[J]. Transactions of the CSAE, 2018,34(11): 96-104.

    [7]扁青永,王振華,胡家?guī)洠?水肥供應(yīng)對(duì)南疆沙區(qū)滴灌紅棗生理、生長(zhǎng)及產(chǎn)量的影響[J].干旱地區(qū)農(nóng)業(yè)研究,2018, 36(4): 165-171.

    BIAN Q, WANG Z, HU J. Effects of water and fertilizer supply on physiology, growth and yield of drip-irrigation jujube in the southern Xinjiang sandy area[J]. Agricultural. Research in the Arid Areas, 2018, 36(4): 165-171.

    [8] DE WIT A, BOOGAARD H, FUMAGALLI D, et al. 25 years of the WOFOST cropping systems model[J]. Agricultural. Systems, 2019, 168: 154-167.

    [9] EWERT F, R?TTER R.P, BINDI M, et al. Crop modelling for integrated assessment of risk to food production from climate change[J]. Environmental. Modelling and Software, 2015, 72: 287-303.

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

    HUANG J, JIA S, MA H. Dynamic simulation of growth process of winter wheat in main production areas of China based on WOFOST model[J]. Transactions of the CSAE, 2017, 33(10): 222-228.

    [11]閆鈴.基于CERES-Wheat模型的小麥生長(zhǎng)發(fā)育過(guò)程模擬研究[D].楊凌:西北農(nóng)林科技大學(xué),2020.

    YAN L. Simulation of wheat growth and development process based on the CERES-Wheat model[D]. Yang ling: Northwest A&F University, 2020.

    [12] FISHMAN S, GENARD M. A biophysical. model of fruit growth: Simulation of seasonal. and diurnal. dynamics of mass[J]. Plant, Cell & Environment, 1998, 21:739-752.

    [13] LOPEZ G, FAVREAU R R, SMITH C, et al. L-PEACH: A computer-based model to understand how peach trees grow [J]. Horttechnology, 2010, 20: 983-990.

    [14] LESCOURRET F, BLECHER N, HABIB R, et al. Development of a simulation model for studying kiwi fruit orchard management[J]. Agricultural. Systems, 1999, 59:215-239.

    [15] GREEN S, MCNAUGHTON K, WUNSCHE J N, et al. Modeling light interception and transpiration of apple tree canopies[J]. Agronomy Journal, 2003, 95(6): 1380-1387.

    [16] COSTES E, SMITH C, RENTON M. MAppleT: Simulation of apple tree development using mixed stochastic and biomechanical. models[J]. Functional. Plant Biology, 2008, 35: 936-950.

    [17]董智強(qiáng),王萌萌,李鴻怡,等.WOFOST模型對(duì)山東省夏玉米發(fā)育期與產(chǎn)量模擬的適用性評(píng)價(jià)[J].作物雜志,2019(5): 159-165.

    DONG Z, WANG M, LI H. Applicability assessment of WOFOST model of growth and yield of summer maize in Shandong province[J]. Crops, 2019(5): 159-165.

    [18] KULIG B, BARBARA S, AGNIESZKA K, et al. The use of the WOFOST model to simulate water-limited yield of early potato cultivars[J]. Agronomy, 2020, 10 (1): ID 81.

    [19] BAI T, WANG T, ZHANG N, et al. Growth simulation and yield prediction for perennial. jujube fruit tree by integrating age into the WOFOST model[J]. Journal. of Integrative Agriculture, 2020, 19(3): 721-734.

    [20] DAI Z, FEI L, HUANG D, et al. Coupling effects of irrigation and nitrogen levels on yield, water and nitrogen use efficiency of surge-root irrigated jujube in a semiarid region[J]. Agricultural. Water Management,2019,213: 146-154.

    [21]李朝陽(yáng),王興鵬,楊玉輝.滴灌帶鋪設(shè)模式對(duì)成齡棗樹(shù)根系再分布及產(chǎn)量的影響[J].灌溉排水學(xué)報(bào). 2021,40(5): 1-7.

    LI Z, WANG X, YANG Y. Location of drip irrigation tapes affects root redistribution and yield of mature jujube[J]. Journal. of Irrigation and Drainage, 2021, 40 (5): 1-7.

    [22]孫愛(ài)良.氣溫變化對(duì)黃驊市冬棗生產(chǎn)的影響[J].貴州農(nóng)業(yè)科學(xué),2019, 47(6): 88-91.

    SUN A. Effect of temperature change on winter jujube production in Huanghua city [J]. Guizhou Agricultural. Sciences, 2019,47(6): 88-91.

    [23] DE WIT A, BOOGAARD H, SUPIT I. System description of the WOFOST 7.2. cropping system model[EB/OL]. [2021-05-15]. https://www.wur.nl/en/Research-Results/Research-Institutes/Environmental-Res- earch/Facilities-Tools/Software-models-and-databases/WOFOST/Documentation-WOFOST.htm

    [24]王慶江,溫陟良,賈彥麗.贊皇大棗幼樹(shù)葉片光合特性的研究[J].河北農(nóng)業(yè)大學(xué)學(xué)報(bào),2002(81): 120-121.

    WANG Q, WEN S, JIA Y. Study on the photosynthetic characteristics of young Jujube[J]. Journal. of Hebei Agricultural. University, 2002(Sl): 120-121.

    [25] HUANG J, JOSE L G, HUANG H, et al. Assimilation of remote sensing into crop growth models: Currentstatus and perspectives[J]. Agricultural. and Forest Meteorology, 2019, 276-277: ID 107609.

    [26] BAI T, WANG S, BENOIT MERCATORIS, et al. Assimilation of remotely-sensed LAI into WOFOST model with the SUBPLEX algorithm for improving the field-scale jujube yield forecasts[J]. Remote Sensing,2019, 11(16): ID 1945.

    [27] HUANG J, TIAN L, LIANG S, 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.

    [28] CURNEL Y, DE WIT A J W, DUVEILLER G, et al. Potential. performances of remotely sensed LAI assimilation in WOFOST model based on an OSS experiment[J]. Agricultural. and Forest Meteorology,2011, 151: 1843-1855.

    [29] BAI T, WANG S, MERCATORIS B, et al. Assimilation of remotely-sensed LAI into WOFOST model with the SUBPLEX algorithm for improving the field-scale jujube yield forecasts[J]. Remote Sensing, 2019,11(16): ID 1945.

    Dynamic Simulation of Jujube Tree Growth and Water Use Evaluation Based on the Calibrated WOFOST Model

    BAI Tiecheng1*, WANG Tao1, ZHANG Nannan1,2

    (1. College of Information Engineering/Southern Xinjiang Research Center Information Technology in Agriculture. Tarim University. Alaer 843300, China; 2. College of Information and Electrical. Engineering. China Agricultural. University, Beijing 100083, China)

    Abstract: Irrigation schemes determined based on statistical. analysis of field trials are usually only applicable to specific soils and meteorological. environments. It is difficult to quantitatively analyze the impact of irrigation strategies on the growth of jujube trees. In order to realize the quantitative analysis of the influence of temperature, light and water resources on the growth of fruit trees, WOrldFOod Studies (WOFOST) model parameters were calibrated to simulate the jujube tree growth and water migration process. Firstly, the observed data obtained from field trials in 2016 and 2017 were used to calibrate the phenology development, initialization, green leaf, CO assimilation, dry matter partitioning, respiration, and water use parameters of the WOFOST model. Secondly, the time series of total. above-ground biomass, leaf area index (LAI) and soil moisture content in field trials were dynamically simulated, and accuracy verification and analysis were also performed. Finally, the maximum LAI,yield, actual. evapotranspiration (ET) and water use efficiency (WUE) data of 55 orchards were employed to evaluate the performance of the calibrated model at the county scale. The results showed that the coefficient of determination R of TAGP simulated in the field test area was between 0.92 and 0.98, and the normalized root mean square error (NRMSE) was between 8.7% and 20.5%, the R of simulated LAI ranged from 0.79 to 0.97, and the NRMSE ranged from 8.3% to 21.1%. The R of the simulated soil moisture content was between 0.29 and 0.75, and the NRMSE ranged from 4.1% and 6.1%, The model could well simulate the time series of jujube tree growth dynamics and soil moisture content changes. At the county scale, the R between the simulated and measured maximum LAI were 0.64 and 0.78, and the NRMSE were 13.3% and 10.7% in 2016 and 2017, respectively. The simulated yield showed R value of 0.48 and 0.60, and NRMSE of 12.1% and 11.9%, respectively. RMSE of the simulated versus measured ET were 36.1 mm (7.9%) and 30.8 mm (7.4%), respectively. The model also showed high WUE simulation accuracy (10%<NRMSE<20%) with RMSE values of 0.23 and 0.28 kg/m in 2016 and 2017, respectively. In short, WOFOST model achieved accurate simulation of jujube tree growth and water transport at the field and county scales, which may provide new ideas for the quantitative and mechanism analysis of the coupled effects of soil, weather, irrigation management and jujube tree growth.

    Key words: crop growth model; parameter calibration; WOFOST; jujube tree; water use efficiency

    (登陸www.smartag.net.cn免費(fèi)獲取電子版全文)

    *通訊作者:白鐵成(1982—),男,博士,教授,研究方向?yàn)樽魑锷L(zhǎng)模擬和遙感同化方法。電話(huà):15299409008。E-mail:baitiechengl983@163.com。

    猜你喜歡
    水分利用效率棗樹(shù)
    大棗樹(shù)
    七八月棗樹(shù)管理莫放松
    棗樹(shù)綠盲蝽象的發(fā)生與防治
    棗樹(shù)的故事
    作品(2018年10期)2018-10-19 11:30:56
    葡萄大棚滴灌耗水規(guī)律研究
    旱澇交替對(duì)水稻生理生長(zhǎng)及水分利用效率的影響
    玉米根—冠及葉片水分利用效率對(duì)土壤水分的響應(yīng)
    三個(gè)引進(jìn)歐洲李品種光合特性對(duì)比研究
    微生物菌肥對(duì)草莓水分利用效率和產(chǎn)量的影響
    濱海設(shè)施黃瓜的減蒸技術(shù)研究
    精品午夜福利在线看| 亚洲国产欧美人成| 精品少妇久久久久久888优播| 日本一二三区视频观看| 亚洲人成网站在线播| 久久精品国产亚洲av天美| 秋霞伦理黄片| 欧美精品一区二区大全| 亚洲国产av新网站| 嫩草影院新地址| 一级毛片我不卡| 欧美zozozo另类| 2021少妇久久久久久久久久久| 免费看日本二区| 黄色欧美视频在线观看| 久久久久久久久久久免费av| 午夜日本视频在线| 亚洲精品乱久久久久久| 王馨瑶露胸无遮挡在线观看| 中文字幕制服av| 在线播放无遮挡| 亚洲欧美一区二区三区黑人 | 国产毛片a区久久久久| 亚洲一区二区三区欧美精品 | 亚洲丝袜综合中文字幕| 日本-黄色视频高清免费观看| 日韩伦理黄色片| 精品久久国产蜜桃| 97精品久久久久久久久久精品| 极品少妇高潮喷水抽搐| 最近中文字幕2019免费版| 精华霜和精华液先用哪个| 精品一区在线观看国产| 麻豆精品久久久久久蜜桃| 别揉我奶头 嗯啊视频| 亚洲精品一二三| 白带黄色成豆腐渣| 免费看光身美女| 亚洲最大成人av| 国产综合精华液| 天天躁日日操中文字幕| 能在线免费看毛片的网站| 亚洲人成网站在线播| 夜夜看夜夜爽夜夜摸| 国产一区二区三区综合在线观看 | 在线 av 中文字幕| 国产精品99久久99久久久不卡 | 国产黄色免费在线视频| 九色成人免费人妻av| 亚洲四区av| 新久久久久国产一级毛片| 国产白丝娇喘喷水9色精品| 精品国产乱码久久久久久小说| 欧美日韩视频高清一区二区三区二| 免费在线观看成人毛片| 欧美人与善性xxx| 永久免费av网站大全| 日韩人妻高清精品专区| 久久人人爽人人片av| 丰满少妇做爰视频| 欧美激情在线99| 在线观看免费高清a一片| 欧美亚洲 丝袜 人妻 在线| 最近最新中文字幕大全电影3| 免费看av在线观看网站| 国产成人免费无遮挡视频| 欧美日韩国产mv在线观看视频 | 综合色av麻豆| 国产精品女同一区二区软件| 91精品伊人久久大香线蕉| 国语对白做爰xxxⅹ性视频网站| 国产 精品1| 干丝袜人妻中文字幕| 蜜桃久久精品国产亚洲av| 在线观看人妻少妇| 国模一区二区三区四区视频| 91精品伊人久久大香线蕉| 国精品久久久久久国模美| 国产黄a三级三级三级人| 国产成年人精品一区二区| 亚洲美女视频黄频| 亚洲国产欧美人成| 亚洲av免费高清在线观看| 久久久久九九精品影院| 精品人妻一区二区三区麻豆| 人体艺术视频欧美日本| eeuss影院久久| 久久6这里有精品| 人妻少妇偷人精品九色| 秋霞伦理黄片| 91久久精品电影网| 一级二级三级毛片免费看| 国产片特级美女逼逼视频| 男女无遮挡免费网站观看| 在线a可以看的网站| 国产一区二区在线观看日韩| 搡女人真爽免费视频火全软件| 国内少妇人妻偷人精品xxx网站| 黑人高潮一二区| 天天躁日日操中文字幕| 男人添女人高潮全过程视频| 国产精品不卡视频一区二区| 国产又色又爽无遮挡免| 日本-黄色视频高清免费观看| 日产精品乱码卡一卡2卡三| 成人特级av手机在线观看| 亚洲最大成人手机在线| 成人亚洲精品av一区二区| 亚洲国产成人一精品久久久| 伦精品一区二区三区| 免费高清在线观看视频在线观看| 日日摸夜夜添夜夜爱| 久久精品久久久久久久性| 久久久a久久爽久久v久久| 国产成人freesex在线| 可以在线观看毛片的网站| 日日啪夜夜撸| 久久热精品热| 亚洲欧洲国产日韩| 又大又黄又爽视频免费| 亚洲综合色惰| 一级毛片 在线播放| 欧美 日韩 精品 国产| 国产在视频线精品| 国产一区有黄有色的免费视频| 精品酒店卫生间| 成人毛片a级毛片在线播放| 只有这里有精品99| 69av精品久久久久久| 国产视频内射| .国产精品久久| 日本午夜av视频| 搡老乐熟女国产| 日韩 亚洲 欧美在线| 国产一区有黄有色的免费视频| 国产精品一区二区在线观看99| 我的女老师完整版在线观看| 美女cb高潮喷水在线观看| 婷婷色av中文字幕| 国产探花极品一区二区| 只有这里有精品99| 国产精品久久久久久av不卡| 熟妇人妻不卡中文字幕| 亚洲av日韩在线播放| 国产亚洲午夜精品一区二区久久 | 成年版毛片免费区| 免费av毛片视频| 日韩亚洲欧美综合| 精品久久久久久久人妻蜜臀av| 三级国产精品欧美在线观看| 国产黄片视频在线免费观看| 国产在线一区二区三区精| 日本熟妇午夜| 国产一区亚洲一区在线观看| 久久久久国产精品人妻一区二区| h日本视频在线播放| av免费观看日本| 午夜老司机福利剧场| 欧美成人午夜免费资源| 成年女人看的毛片在线观看| 久久人人爽人人爽人人片va| 久热久热在线精品观看| 观看免费一级毛片| 边亲边吃奶的免费视频| 在线观看国产h片| 女的被弄到高潮叫床怎么办| 亚洲欧美一区二区三区黑人 | 久久精品国产亚洲av涩爱| 汤姆久久久久久久影院中文字幕| 2018国产大陆天天弄谢| 日本黄色片子视频| 久久久久久久久久久免费av| 在线观看三级黄色| 韩国av在线不卡| 熟女av电影| 免费高清在线观看视频在线观看| 美女cb高潮喷水在线观看| 国产日韩欧美在线精品| 少妇丰满av| 一个人观看的视频www高清免费观看| 免费在线观看成人毛片| 久久久久国产网址| 精品久久久久久久久亚洲| 亚洲欧美日韩卡通动漫| 国产一区有黄有色的免费视频| 午夜日本视频在线| 国产成人免费无遮挡视频| 国产免费又黄又爽又色| 男女啪啪激烈高潮av片| 卡戴珊不雅视频在线播放| 97热精品久久久久久| 人妻系列 视频| 国产白丝娇喘喷水9色精品| 一区二区av电影网| 男人爽女人下面视频在线观看| 久久精品国产a三级三级三级| h日本视频在线播放| 精品久久久久久久久av| 国产一级毛片在线| 成人毛片60女人毛片免费| av一本久久久久| 久久久久久久国产电影| 欧美xxⅹ黑人| 嫩草影院新地址| 最近2019中文字幕mv第一页| 成人综合一区亚洲| 日韩一本色道免费dvd| 精品午夜福利在线看| 女人被狂操c到高潮| 五月开心婷婷网| 免费观看a级毛片全部| 99久久精品国产国产毛片| 国产淫语在线视频| 国产精品成人在线| 亚洲精品国产色婷婷电影| 午夜激情福利司机影院| 日韩一区二区视频免费看| 又爽又黄无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 国产精品99久久久久久久久| 免费av不卡在线播放| 麻豆久久精品国产亚洲av| 三级国产精品欧美在线观看| av.在线天堂| 一本一本综合久久| 男人和女人高潮做爰伦理| 国产精品蜜桃在线观看| 一个人看视频在线观看www免费| 内射极品少妇av片p| 可以在线观看毛片的网站| 色综合色国产| 人妻制服诱惑在线中文字幕| 久久久午夜欧美精品| 激情 狠狠 欧美| 最后的刺客免费高清国语| 热re99久久精品国产66热6| 夜夜爽夜夜爽视频| 你懂的网址亚洲精品在线观看| 免费播放大片免费观看视频在线观看| xxx大片免费视频| 久久影院123| 久久久久网色| 亚洲精品乱码久久久v下载方式| 一个人看的www免费观看视频| 欧美日韩在线观看h| 99九九线精品视频在线观看视频| 麻豆成人午夜福利视频| 成年av动漫网址| 97在线人人人人妻| 男人添女人高潮全过程视频| 网址你懂的国产日韩在线| 国产成人精品久久久久久| 婷婷色综合大香蕉| 在线观看国产h片| 秋霞在线观看毛片| 乱码一卡2卡4卡精品| 久久精品国产亚洲av涩爱| 99精国产麻豆久久婷婷| 在线a可以看的网站| 少妇人妻久久综合中文| 成人漫画全彩无遮挡| 婷婷色综合www| 夜夜爽夜夜爽视频| 22中文网久久字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 在线观看一区二区三区激情| 国产伦理片在线播放av一区| 大话2 男鬼变身卡| 26uuu在线亚洲综合色| 97超视频在线观看视频| 91在线精品国自产拍蜜月| 亚洲精品第二区| 免费人成在线观看视频色| 久久久久国产网址| 国产午夜福利久久久久久| 久久精品久久久久久噜噜老黄| 女人十人毛片免费观看3o分钟| 久久久久网色| 中文字幕人妻熟人妻熟丝袜美| av在线亚洲专区| 久久久久久久午夜电影| 精品久久久精品久久久| 亚洲av不卡在线观看| 十八禁网站网址无遮挡 | 国产成人福利小说| 91精品伊人久久大香线蕉| 国产片特级美女逼逼视频| 亚洲美女搞黄在线观看| 亚洲在线观看片| 十八禁网站网址无遮挡 | 综合色丁香网| 国产亚洲一区二区精品| 新久久久久国产一级毛片| 人妻制服诱惑在线中文字幕| 1000部很黄的大片| 亚洲怡红院男人天堂| 亚洲精品影视一区二区三区av| 婷婷色av中文字幕| 麻豆成人av视频| 26uuu在线亚洲综合色| 免费黄色在线免费观看| 欧美激情久久久久久爽电影| 2018国产大陆天天弄谢| 中文在线观看免费www的网站| 国产女主播在线喷水免费视频网站| 亚洲精品色激情综合| 亚洲成人一二三区av| 国产有黄有色有爽视频| 欧美国产精品一级二级三级 | 菩萨蛮人人尽说江南好唐韦庄| 老司机影院成人| 色5月婷婷丁香| 国产一区亚洲一区在线观看| 亚洲国产成人一精品久久久| 国产 一区精品| 精品久久久精品久久久| 18+在线观看网站| 精品久久久久久电影网| 热99国产精品久久久久久7| 亚洲国产精品999| 中国三级夫妇交换| 麻豆精品久久久久久蜜桃| 插逼视频在线观看| 一级毛片久久久久久久久女| 成人欧美大片| 高清日韩中文字幕在线| 极品教师在线视频| 在线观看av片永久免费下载| 久久久色成人| 国产精品一区二区三区四区免费观看| 欧美性猛交╳xxx乱大交人| 97在线人人人人妻| 亚洲欧洲国产日韩| 亚洲精品成人av观看孕妇| 看免费成人av毛片| 日本爱情动作片www.在线观看| 老女人水多毛片| 久久热精品热| 国产欧美另类精品又又久久亚洲欧美| 久久人人爽av亚洲精品天堂 | 国产一区二区三区综合在线观看 | 亚洲精品国产成人久久av| 好男人在线观看高清免费视频| 极品教师在线视频| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 精华霜和精华液先用哪个| 少妇猛男粗大的猛烈进出视频 | 亚洲内射少妇av| 特大巨黑吊av在线直播| 97人妻精品一区二区三区麻豆| 国产亚洲精品久久久com| 欧美激情久久久久久爽电影| 欧美另类一区| 国产乱人偷精品视频| av国产免费在线观看| 精品一区二区三区视频在线| 国产免费福利视频在线观看| 嫩草影院入口| av天堂中文字幕网| 久久99热这里只频精品6学生| 欧美日本视频| 国产亚洲91精品色在线| 亚洲三级黄色毛片| 精品一区二区免费观看| 亚洲国产精品国产精品| av国产免费在线观看| 91精品国产九色| 99re6热这里在线精品视频| 偷拍熟女少妇极品色| 成年女人在线观看亚洲视频 | 久久精品国产鲁丝片午夜精品| av播播在线观看一区| 久久久久精品性色| 欧美xxⅹ黑人| 国产精品一区www在线观看| 波野结衣二区三区在线| av卡一久久| 色5月婷婷丁香| 舔av片在线| 欧美一级a爱片免费观看看| 最近最新中文字幕免费大全7| 久久久久久久大尺度免费视频| 国产精品偷伦视频观看了| 黄色视频在线播放观看不卡| 欧美日韩一区二区视频在线观看视频在线 | 美女内射精品一级片tv| 搞女人的毛片| 91精品伊人久久大香线蕉| 国产日韩欧美在线精品| 精品久久久久久久久亚洲| 国产午夜精品久久久久久一区二区三区| 七月丁香在线播放| 卡戴珊不雅视频在线播放| 晚上一个人看的免费电影| 久久女婷五月综合色啪小说 | 亚洲国产av新网站| 99久久精品国产国产毛片| 成年女人在线观看亚洲视频 | 欧美少妇被猛烈插入视频| 国产乱来视频区| 国产欧美亚洲国产| 亚洲av.av天堂| 色视频www国产| 久久鲁丝午夜福利片| xxx大片免费视频| 国产亚洲一区二区精品| 全区人妻精品视频| 欧美另类一区| 嘟嘟电影网在线观看| 精品人妻视频免费看| 性色avwww在线观看| 国产精品人妻久久久久久| 国产精品成人在线| 男插女下体视频免费在线播放| 午夜视频国产福利| 久久久成人免费电影| 亚洲精品第二区| 边亲边吃奶的免费视频| 三级国产精品片| 在线看a的网站| 国产美女午夜福利| 国产女主播在线喷水免费视频网站| 熟女人妻精品中文字幕| 国产爱豆传媒在线观看| 婷婷色综合www| 99久久精品一区二区三区| 永久免费av网站大全| 国产免费福利视频在线观看| 亚洲怡红院男人天堂| 国产91av在线免费观看| 天天躁夜夜躁狠狠久久av| 久久午夜福利片| 99久久九九国产精品国产免费| 97超视频在线观看视频| 欧美成人午夜免费资源| 91aial.com中文字幕在线观看| 男男h啪啪无遮挡| 尤物成人国产欧美一区二区三区| 边亲边吃奶的免费视频| 最后的刺客免费高清国语| 中文资源天堂在线| 久久精品综合一区二区三区| 国产视频内射| 国产午夜福利久久久久久| av国产久精品久网站免费入址| 精品人妻一区二区三区麻豆| 免费观看a级毛片全部| 久久久久精品性色| 一二三四中文在线观看免费高清| 日韩欧美一区视频在线观看 | 男男h啪啪无遮挡| 成人综合一区亚洲| 自拍偷自拍亚洲精品老妇| av女优亚洲男人天堂| 欧美日本视频| 99久久中文字幕三级久久日本| 男人狂女人下面高潮的视频| 日本免费在线观看一区| 久久久精品94久久精品| 亚洲av一区综合| 亚洲综合精品二区| 小蜜桃在线观看免费完整版高清| 美女主播在线视频| 狂野欧美激情性bbbbbb| 网址你懂的国产日韩在线| 一级毛片电影观看| 真实男女啪啪啪动态图| 亚洲欧美精品专区久久| 国产免费视频播放在线视频| 欧美97在线视频| 少妇人妻 视频| 99热这里只有是精品在线观看| 欧美精品人与动牲交sv欧美| 久久99热这里只频精品6学生| 秋霞伦理黄片| 亚洲欧美日韩无卡精品| 2021少妇久久久久久久久久久| 最近中文字幕2019免费版| 久久人人爽人人爽人人片va| 禁无遮挡网站| 久久久久久久久久成人| 高清av免费在线| 久久久久久久久久久丰满| 久久久久精品久久久久真实原创| 国产黄色视频一区二区在线观看| 夫妻性生交免费视频一级片| 亚洲精品久久午夜乱码| 日本一本二区三区精品| 哪个播放器可以免费观看大片| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载 | 日韩成人av中文字幕在线观看| 另类亚洲欧美激情| 日韩伦理黄色片| 水蜜桃什么品种好| 欧美三级亚洲精品| 亚洲av电影在线观看一区二区三区 | a级毛色黄片| 亚洲av欧美aⅴ国产| 亚洲欧美清纯卡通| 成人亚洲欧美一区二区av| 国产精品女同一区二区软件| 国产成人午夜福利电影在线观看| 丰满人妻一区二区三区视频av| 新久久久久国产一级毛片| 国产精品成人在线| 亚洲精品日本国产第一区| 天天一区二区日本电影三级| 亚洲自偷自拍三级| 国产中年淑女户外野战色| 两个人的视频大全免费| 天天躁日日操中文字幕| 欧美高清成人免费视频www| 国产成人福利小说| 亚洲不卡免费看| 亚洲成人精品中文字幕电影| 国产视频内射| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 精品少妇黑人巨大在线播放| 啦啦啦在线观看免费高清www| 精品亚洲乱码少妇综合久久| 久久精品国产自在天天线| 建设人人有责人人尽责人人享有的 | 高清午夜精品一区二区三区| 国产在视频线精品| 免费高清在线观看视频在线观看| 精品人妻一区二区三区麻豆| 大片电影免费在线观看免费| 搞女人的毛片| 好男人视频免费观看在线| 亚洲精品乱码久久久久久按摩| 少妇熟女欧美另类| 看非洲黑人一级黄片| 亚洲人成网站在线播| 一边亲一边摸免费视频| 久久鲁丝午夜福利片| 国语对白做爰xxxⅹ性视频网站| 日产精品乱码卡一卡2卡三| 又爽又黄无遮挡网站| 亚洲av免费在线观看| 日韩一区二区视频免费看| 我要看日韩黄色一级片| 久久鲁丝午夜福利片| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区性色av| 日韩一区二区三区影片| 色吧在线观看| 欧美日本视频| 好男人视频免费观看在线| 国产爱豆传媒在线观看| 国产片特级美女逼逼视频| 精品熟女少妇av免费看| 成年免费大片在线观看| 人体艺术视频欧美日本| 天堂俺去俺来也www色官网| 中文天堂在线官网| 午夜爱爱视频在线播放| 少妇熟女欧美另类| 国产av码专区亚洲av| 91在线精品国自产拍蜜月| 熟女av电影| 我要看日韩黄色一级片| 国产精品99久久久久久久久| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 超碰av人人做人人爽久久| 精品视频人人做人人爽| 国产精品无大码| 亚洲av免费高清在线观看| 亚洲成人av在线免费| 男女边吃奶边做爰视频| 国产成人精品一,二区| 日本与韩国留学比较| 综合色丁香网| 久久久久久久久久成人| 嫩草影院入口| 夫妻性生交免费视频一级片| 大香蕉97超碰在线| 欧美区成人在线视频| 国语对白做爰xxxⅹ性视频网站| 交换朋友夫妻互换小说| 亚洲四区av| 天天一区二区日本电影三级| 久久综合国产亚洲精品| 尾随美女入室| 久久韩国三级中文字幕| 伦理电影大哥的女人| 女人十人毛片免费观看3o分钟| 三级经典国产精品| 欧美性感艳星| 最近最新中文字幕大全电影3| 国国产精品蜜臀av免费| 80岁老熟妇乱子伦牲交| 啦啦啦中文免费视频观看日本| 国产一区二区三区综合在线观看 | 最后的刺客免费高清国语| 午夜福利在线观看免费完整高清在| 亚洲va在线va天堂va国产| 精品少妇黑人巨大在线播放| 国产精品久久久久久av不卡| av天堂中文字幕网| 好男人在线观看高清免费视频| 免费观看av网站的网址| 18禁在线无遮挡免费观看视频| 国产av码专区亚洲av| 大片免费播放器 马上看| 成人高潮视频无遮挡免费网站| 欧美成人一区二区免费高清观看| 伦理电影大哥的女人| 亚洲av免费高清在线观看| 寂寞人妻少妇视频99o| 亚洲欧美一区二区三区国产| 国产精品成人在线|