基于WOFOST模型的蘇南地區春小麥種植適應性分析

江曉東 張濤 陳佳鑫 楊沈斌 李旭冉 吳可人 楚淏然

摘要：【目的】研究春小麥在蘇南地區的種植適應性，為增加該地區周年糧食生產的穩定性提供理論依據。【方法】基于2016—2017年的田間分期播種試驗和WOFOST作物生長模型（簡稱WOFOST模型），采用數值模擬方法，分析江蘇南部代表地區南京地區1—4月不同播期春小麥的生長發育動態和產量表現。以1980—2010年氣象數據驅動WOFOST模型，對春小麥產量進行動態模擬，分析最佳播期，并計算最佳播期的適宜播種量。【結果】在1—4月隨著播期的推遲，春小麥的生育期長度縮短，其中出苗—開花期階段最大縮短23 d，開花—成熟期最大縮短8 d，出苗—開花階段縮短的時間大于開花—成熟階段，導致春小麥葉片和莖稈的干物質積累量明顯減少，產量降低。在冬小麥無法播種的條件下，南京地區春小麥的適宜播種時間為1月1—20日，在該時間段內播種，合理的播種量為180 kg/ha，最高產量為4124.80 kg/ha，生育期長度為146 d，對下季水稻種植和生長無影響。【結論】在蘇南地區種植春小麥具有可行性，在冬前無法正常播種冬小麥的情況下，可將種植春小麥作為備選方案。

關鍵詞： WOFOST模型;春小麥;播期;產量;蘇南地區

中圖分類號： S314? ? ? ? ? ? ? ? ? ? ? 文獻標志碼： A 文章編號：2095-1191（2020）02-0335-07

Adaptability of spring wheat planting in the south of Jiangsu

based on WOFOST model

JIANG Xiao-dong1，ZHANG Tao1， CHEN Jia-xin2， YANG Shen-bin1， LI Xu-ran1，

WU Ke-ren1， CHU Hao-ran1

（1Jiangsu Key Laboratory of Agricultural Meteorology/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Nanjing University of Information Science & Technology， Nanjing? 210044， China; Cangzhou Meteorological Bureau， Cangzhou， Heibei? 061001）

Abstract：【Objective】The adaptability of spring wheat in southern Jiangsu was studied to provide reference for stabi-lity of crop production in this area. 【Method】Based on the field staging experimental data from 2016 to 2017， the growth and development and yield of spring wheat in Nanjing during January and April， a representative area of southern Jiangsu Province， were analyzed through numerical simulation. The WOFOST crop growth model（herein after referred to as WOFOST model） was employed in the simulation. It was also driven by meteorological data from 1980 to 2010 for the analysis of optimum sowing date and sowing rate. 【Result】The results showed that the growth period of spring wheat was shortened with the delay of sowing date from January to April， the maximum shortening time was 23 d from emergence to anthesis stage， and the maximum shortening time was 8 d from anthesis to maturity stage. The shortening of emergence-anthesis stage was greater than that of anthesis-maturity stage， which led to a reduction of dry matter accumulation in both leaves and stems as well as yield of spring wheat. If winter wheat could not be sown， the suitable sowing date of spring wheat in Nanjing was from January 1 to 20 and the reasonable seeding rate was 180 kg/ha， with the highest yield predic-ted to be 4124.80 kg/ha and growth period to be 146 d， which showed no evident impact on the rice planting of the next season. 【Conclusion】The results indicate that it is feasible to plant spring wheat in the south of Jiangsu， which can be adopted as an alternative in case winter wheat is not sown before winter.

Key words： WOFOST model; spring wheat; sowing date; yield; south of Jiangsu

Foundation item： General Program of National Natural Science Foundation of China（41875140）; Open Research Project of Key Laboratory of Agrometerological Ensuring and Applied Technique，China Meteorological Administration（AMF201602）

0 引言

【研究意義】小麥和水稻是江蘇南部的主要農作物，稻麥輪作是主要耕作方式。近年來，為提高水稻的品質和產量，蘇南地區晚（粳）稻推廣種植面積不斷擴大，水稻成熟期延遲至11月中旬，導致該地區冬小麥播種期推遲到11月中下旬甚至12月，冬小麥晚播現象普遍發生，稻麥生長爭季節的茬口矛盾突出（葛勝等，2009;高德榮等，2014）。隨著經濟的發展和城市化進程的加快，大量農村勞動力進入城市，農業勞動力供給不足，加劇了水稻收獲與小麥播種的茬口矛盾;此外，蘇南地區小麥播種期間經常發生連陰雨天氣，加劇了冬小麥晚播現象，如遇“雪壓稻”等災害性天氣發生，則會出現冬季撂荒現象（徐德利等，2016）。相對于冬小麥品種，春小麥品種在0～12 ℃下經過5～15 d即可通過感溫階段，在8～12 h光照條件下經過16 d以上即可通過光照階段而抽穗（楊文鈺和屠乃美，2011），具有生長周期短、播種安排靈活、適應性強的特點。因此，研究春小麥在蘇南地區的栽培適宜性，可增強該地區糧食生產的穩定性。【前人研究進展】春小麥或春性較強的冬小麥品種是小麥生產中常用的應急備選品種。馬香花等（2003）研究表明，春小麥可作為特殊年份的備用品種在鄭州地區春季種植。張耀蘭等（2012）指出，在淮北遲播小麥的品種布局上，可考慮選用一些播期彈性大、半冬性、豐產性好的小麥品種。蘇南地區地處亞熱帶季風氣候區，冬春季節氣溫較高，最冷月平均氣溫高于0 ℃，完全滿足春小麥的春化與種植要求，春小麥種植具有可行性。WOFOST作物生長模型（簡稱WOFOST模型）是由荷蘭瓦赫寧根大學和世界糧食研究中心共同研發的動態解釋型作物模型，該模型通過輸入特定的土壤和氣象條件，基于作物的生理生態過程，以天為步長，可對一年生作物的生長發育過程和產量進行模擬。在我國小麥的研究中，WOFOST模型得到廣泛應用。鄔定榮等（2003）、黃健熙等（2017）、張弘等（2019）通過WOFOST模型模擬了河南、華北及全國的冬小麥生長過程及產量，證實WOFOST模型可應用于不同面積小麥的生長發育和產量模擬;張鐵楠等（2015a）利用WOFOST模型模擬東北春小麥的生長及產量。WOFOST模型也可用于小麥的災害評估，張雪芬等（2006）利用WOFOST模型定量評估了商丘市冬小麥晚霜凍災害，結果表明該評估方法更具生物學意義且結果更準確;張建平等（2013）基于WOFOST模型對華北地區干旱災害進行評估，發現該模型對干旱評估具有較高的準確性。WOFOST模型能模擬小麥的水分平衡，張鐵楠等（2015b）使用WOFOST模型對東北春小麥的水分平衡進行了模擬與驗證，結果表明模型對水分平衡的模擬效果較好。【本研究切入點】春小麥在我國的種植區域主要分布在東北和西北等地區，在蘇南地區基本沒有種植（楊文鈺和屠乃美，2011;俄有浩等，2013）。盡管蘇南地區存在部分弱春性小麥品種，但小麥的適播期在11月左右，至今鮮見該地區春小麥種植適應性的研究報道。【擬解決的關鍵問題】蘇南各地區地理經緯度差異小，氣候相似，因此本研究以南京地區為代表，采用春小麥田間分期播種試驗與WOFOST模型模擬相結合，利用田間試驗獲取相關作物數據對WOFOST模型進行參數本地化，進而利用參數本地化的WOFOST模型對不同播種時期的春小麥生長進行模擬，開展春小麥在蘇南地區種植適應性研究，為增加蘇南地區周年糧食生產的穩定性及保障糧食安全提供理論基礎。

1 材料與方法

1. 1 試驗材料

供試春小麥品種為吉春24，由吉林省農業科學院提供。1980—2010年南京的歷史氣象數據來源于中國氣象數據網（http：//data.cma.cn/）。2016—2017年田間試驗期間的氣象數據及土壤數據，由南京信息工程大學農業氣象試驗站提供。

1. 2 試驗方法

試驗于2016和2017年在南京信息工程大學農業氣象試驗站（東經118°42′17″，北緯32°12′24″）進行，試驗地海拔約22 m，地處蘇南地區，屬亞熱帶季風氣候，多年平均降水量約1100 mm，多年平均溫度15.6 ℃。2016和2017年度小麥生長期間，總降水量分別為479和404 mm，日平均輻射分別為9.2和10.4 MJ/（m2·d），日平均氣溫為13.2和14.3 ℃，兩年氣象條件較相似。供試農田0～20 cm土壤質地為壤質黏土，黏粒含量26.1%，土壤pH 6.2，有機碳、全氮含量分別為19.4和11.5 g/kg。

試驗設6個播期處理，分別為Q1（2016年3月29日）、Q2（2016年4月5日）、Q3（2016年4月12日）、Q4（2016年12月16日）、Q5（2017年1月13日）和Q6（2017年2月19日）。春小麥采用條播種植，行距25 cm，種植密度300株/m2，小區面積3 m×3 m，采用隨機區組排列，3次重復。試驗田肥料為45%復合肥（N∶P∶K=15∶15∶15），用量800 kg/ha，基肥∶追肥=1∶1，基肥在土壤耕作時施用，追肥在春小麥拔節時施用，田間管理同當地高產田。

1. 3 測定項目及方法

1. 3. 1 春小麥發育期觀測 觀測記錄不同播期春小麥的出苗期、拔節期、抽穗期、開花期和成熟期。

1. 3. 2 春小麥干物質積累量和莖蘗動態測定 在每個試驗小區選擇長勢均勻的1 m2樣點做標記，在春小麥出苗期、拔節期、抽穗期、開花期和成熟期定點調查小麥莖蘗數，然后在小區其他地點隨機連續取樣（40個莖蘗），樣品葉片和莖稈分開，75 ℃烘干至恒重，測定干物質積累量。

1. 3. 3 春小麥產量及其構成因素測定 成熟期每小區選擇長勢均勻的2 m2地塊，先清點有效穗數，然后進行人工收割測產、晾曬和脫粒，折算成籽粒水含率為13%時的產量;取40個單莖進行考種，計數穗粒數和千粒重。

1. 4 WOFOST模型參數調整及驗證

為評估模型模擬值與試驗實測值間的擬合程度和模擬效果，采用RMSE（均方根誤差）和EF（模擬性能指數）2個指標進行評價。計算方法如下：

RMSE=[1ni=1nPi-Oi2] （1）

EF=[i=1nOi-O2-i=1nPi-Oi2i=1nOi-O2] （2）

式中，n為樣本容量，O為實測值的平均值，Oi、Pi分別為實測值和模擬值。RMSE越接近0、EF越靠近1，說明模擬效果越好。

在前人對WOFOST模型參數研究（張建平等，2012;張鐵楠等，2015a;黃健熙等，2017）的基礎上，利用2016年度Q1、Q2和Q3播期的試驗數據計算發育階段（DVS）及各時期的比葉面積（SLATB）、葉的分配系數（FLTB）、貯藏器官的分配系數（FOTB）和莖的分配系數（FSTB）。在近似確定部分作物參數后對WOFOST模型參數進行敏感性分析。本研究的敏感性分析主要采用初始變量擾動法，即將某一參數上調10%進行模擬，將其模擬結果與不上調參數的值進行比較和分析，計算其變化百分率，之后選取敏感性較大的參數進行本地化調參，結果見表1，最終模擬過程對AMAXTB、SPAN、Q10、TSUM和TDWI等9個參數的敏感性較強，SPAN甚至超過100.00%。因此，需根據資料對敏感性較強的參數進行調整。主要的參數調整見表2，然后利用2017年度Q4、Q5和Q6播期的試驗數據進行WOFOST模型驗證，使模型可模擬南京地區的春小麥生長發育。

1. 5 統計分析

采用Excel 2016進行數據處理和圖表繪制，采用DPS7.05進行數據方差分析，采用最小顯著性法（LSD）對不同處理間的數據進行5%水平上的統計檢驗。

2 結果與分析

2. 1 WOFOST模型驗證

作物參數調整和模型檢驗是模型應用的前提條件。圖1為參數調整后的WOFOST模型對2017年度Q4、Q5和Q6播期春小麥生育期、地上總干物質的模擬值與實測值對比圖。從圖1可看出，WOFOST模型對春小麥生育期誤差最大值為3 d，最小值為1 d，RMSE和EF分別為1.89 d和0.99，對生育期的模擬效果較好（圖1-A）;WOFOST模型對春小麥地上總干物質模擬中（圖1-B），RMSE為1639 kg/ha，EF為0.95，3個播期的誤差相對較小，可見，參數調整后的WOFOST模型對春小麥生長發育動態的模擬效果較好，各項模擬結果的誤差較小，即WOFOST模型參數本地化成功，可用于蘇南地區春小麥生長發育和產量的模擬。

2. 2 不同播期對春小麥生長發育的影響

利用2016年的氣象資料輸入WOFOST模型，模擬該年條件下不同播期春小麥的生育期、莖干重、葉干重及產量。模擬結果表明（圖2-A），春小麥播種—成熟持續時間隨播期的延遲而縮短，最長的第1播期（播種日序為1 d）為146 d，最短的第12播期（播種日序為111 d）只有78 d，比第1播期縮短68 d;春小麥出苗—開花期階段持續時間在42～65 d內變動，隨播期推遲而縮短，最長的第1播期比最短的第12播期長23 d;相對于出苗—開花期，開花—成熟期持續時間相對穩定，在24～32 d間變動，第1播期比第12播期長8 d。可見，隨播期的推遲，春小麥出苗—開花階段縮短的時間大于開花—成熟階段，春小麥生育期長度的縮短主要是由于出苗—開花期階段的持續天數縮短所造成。

發育期的變化影響了春小麥的干物質積累與分配。圖2-B顯示，春小麥的葉片干物質積累量在987～1849 kg/ha，莖稈干物質積累量在1568～3878 kg/ha，穗干物質積累量在2503～3878 kg/ha，皆隨播期的延遲而減小;葉片干物質積累量最大值（第1播期）比最小值（第12播期）高87.36%，莖稈干物質積累量最大值（第1播期）比最小值（第12播期）高147.32%，產量（穗干重）最大值（第1播期）是最小值（第12播期）的1.91倍，出苗—開花階段是小麥營養器官生長發育的主要階段，播期越早該階段持續的時間越長，越有利于葉片生長、提高葉面積指數、增加莖稈干物質的儲存量，有利于產量的提高，同時，開花—成熟階段持續時間播期早的春小麥高于播期晚的春小麥，該階段的延長有利于小麥籽粒的干物質積累，提高產量。

2. 3 不同播期對春小麥產量的影響

利用參數本地化后的WOFOST模型進行不同播期春小麥產量模擬。根據蘇南地區如12月中旬小麥無法播種將導致無法種植冬小麥的實際情況，將模型模擬春小麥的播種時間確定為：從每年1月1日開始每隔10 d設為一個播期，共設12個播期。圖3為利用1980—2010年歷史氣象資料輸入后WOFOST模型模擬的春小麥產量。從圖3可看出，盡管不同年份氣象因素的變化對春小麥產量有較明顯的影響，但不同播期小麥的產量變化趨勢一致，春小麥的產量隨播期的延遲而下降，最高產量為模擬的第1播期，為4124.80 kg/ha，最低產量為模擬的第12播期，為2379.07 kg/ha，第12播期產量僅為第1播期的57.68%，減產顯著（P<0.05）。從圖3還可看出，在第1～3播期，春小麥產量相對穩定，分別為4124.80、4106.66和3999.22 kg/ha，從第4播期（播種日序為31 d）開始，春小麥產量明顯下降，第4播期產量為3845.40 kg/ha，比第1播期減產6.77%。因此，可認為每年的1月1～20日為春小麥的最適播期。

2. 4 春小麥合理播種量分析

利用參數本地化后的WOFOST模型模擬在合理播期范圍內（播種日序分別為1、11和21 d）春小麥產量與播種量間的關系（圖4）。從圖4可看出，春小麥產量隨播種量的增加而增加，增加速率先快后慢。產量與播種量的關系可用分段函數來擬合，即前半部分用對數擬合，后半部分用增幅較小的線性擬合。以第1播期（播種日序為1 d）為例，以相關系數最大為依據，可計算出分段函數的分界點為播種量210 kg/ha處，分段函數見方程（3），通過計算，可得出最佳的播種量為177 kg/ha。將第2和第3播期的產量與播種量采用相同的方法進行對數擬合，可得到第2、第3播期的播種量分別為182和184 kg/ha。由此得出在合理播期內，最佳播種量為180 kg/ha左右。

Y=[1229lnx-1514.4，? x<2100.83x+4703， x≥210]? ? ? ? （3）

3 討論

播期改變了小麥不同發育階段的長度，進而對小麥的產量產生影響。營養生長和生殖生長的協調發展是小麥高產的保證，張凱等（2012）指出春小麥早播，春小麥播種—抽穗期日數增加，葉面積指數增大，產量提高;史曉芳等（2017）研究表明，小麥適期播種出苗早，營養生長期長，分蘗數和葉面積指數增大。本研究結果與以上結論一致，春小麥適期播種，明顯延長出苗—開花期持續時間，促進了春小麥葉片、莖稈、根系和分蘗的分化和生長，春小麥莖蘗干物質積累高、葉面積指數大、群體質量好，保證了春小麥高產。該階段持續時間的延長還有利于春小麥穗的分化，有利于提高春小麥的穗粒數，與韓金玲等（2011）研究認為正常播期延長了小麥穗分化始期—雌雄蕊分化期的持續時間，穗分化時間增加，增加了麥穗的小穗數的結論一致。同時，適期播種春小麥開花—成熟期持續時間亦有延長，增加了春小麥灌漿持續時間，也促進了產量的提高，與安霞等（2018）的研究結論一致。本研究結果還表明，在適宜的播期內（播種日序1～21 d），隨播期的推后，春小麥最適播種量由177 kg/ha增至184 kg/ha，產量由4124.80 kg/ha降至3999.22 kg/ha，基本保持穩定，說明適當增加播量可提高小麥群體數量，增大葉面積指數和光合生產量，從而彌補播期推遲對小麥產量的影響（屈會娟等，2009;張耀蘭等，2012;史曉芳等，2017;安霞等，2018）。

綜合本研究結果可知，在蘇南地區無法播種冬小麥的情況下，春小麥可作為備播作物進行種植，春小麥的最佳播期為1月1日，全生育期長度為146 d，在5月下旬收獲，不會影響后季水稻生產。在最適播期范圍內（1月上旬—中旬），春小麥播種量為117～180 kg/ha，產量為3999.22～4124.80 kg/ha，可有效保障該地區周年糧食生產的穩定性。本研究只選用了1個春小麥品種，對蘇南地區遲播春小麥的播期和播種量進行相應研究，關于春小麥在該地區的品種篩選、品質表現及適用栽培技術等均有待進一步研究。

4 結論

春小麥適宜在蘇南地區播種，年后春小麥的最佳播種時間為1月上旬～中旬，合理的播種量為180 kg/ha，最高產量為4124.80 kg/ha，生育期長度為146 d，對下季水稻正常栽插無影響。因此，蘇南地區在冬前無法正常播種冬小麥的情況下，可將種植春小麥作為備選方案。

參考文獻：

安霞，張海軍，蔣方山，呂連杰，陳軍. 2018. 播期播量對不同穗型冬小麥群體及子粒產量的影響[J]. 作物雜志，（5）：138-142. [An X，Zhang H J，Jiang F S，Lü L J，Chen J. 2018. Effects of different sowing dates and sowing densities on the population structure and yield of two spike type winter wheats[J]. Crops，（5）：138-142.]

俄有浩，霍治國，馬玉平，徐蕓皎. 2013. 中國北方春小麥生育期變化的區域差異性與氣候適應性[J]. 生態學報，33（19）：6295-6302. [E Y H，Huo Z G，Ma Y P，Xu Y J. 2013. The regional diversity of changes in growing duration of spring wheat and its correlation with climatic adaptation in Northern China[J]. Acta Ecologica Sinica，33（19）：6295-6302.]

高德榮，張曉，康建鵬，別同德，張伯橋，張曉祥，程順和. 2014. 長江中下游麥區小麥遲播的不利影響及育種對策[J]. 麥類作物學報，34（2）：279-283. [Gao D R，Zhang X，Kang J P，Bie T D，Zhang B Q，Zhang X X，Cheng S H. 2014. Negative effects of late sowing on wheat production in middle and lower reaches of Yangtze River Valley and breeding strategies[J]. Journal of Triticeae Crop，34（2）：279-283.]

葛勝，童小云，余媛. 2009. 晚茬小麥生育特性及高產栽培技術初探[J]. 現代農業科學，16（3）：89-90. [Ge S，Tong X Y，Yu Y. 2009. Late crop growth characteristics of wheat and of high-yielding cultivation techniques[J]. Modern Agricultural Sciences，16（3）：89-90.]

韓金玲，楊晴，王文頗，李彥生，周印富. 2011. 播期對冬小麥莖蘗幼穗分化及產量的影響[J]. 麥類作物學報，31（2）：115-119. [Han J L，Yang Q，Wang W P，Li Y S，Zhou Y F. 2011. Effects of sowing date on the caulis and tillers differentiation of young spike and yield in winter wheat[J]. Journal of Triticeae Crops，31（2）：115-119.]

黃健熙，賈世靈，馬鴻元，侯英雨，何亮. 2017. 基于WOFOST模型的中國主產區冬小麥生長過程動態模擬[J]. 農業工程學報，33（10）：222-228. [Huang J X，Jia S L，Ma H Y，Hou Y Y，He L. 2017. 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，33（10）： 222-228.]

馬香花，喬旭，周秋峰，王保林. 2003. 春小麥在鄭州地區春播適應性研究初報[J]. 河南農業科學，（10）：8-9. [Ma X H，Qiao X，Zhou Q F，Wang B L. 2003. Preliminary report on the adaptability of spring wheat in spring in Zhengzhou[J]. Journal of Henan Agricultural Sciences（10）：8-9.]

屈會娟，李金才，沈學善，魏鳳珍，王成雨，郅勝軍. 2009. 種植密度和播期對冬小麥品種蘭考矮早八干物質和氮素積累與轉運的影響[J]. 作物學報，35（1）：124-131. [Qu H J，Li J C，Shen X S，Wei F Z，Wang C Y，Zhi S J. 2009. Effects of plant density and seeding date on accumulation and translocation of dry matter and nitrogen in winter wheat cultivar Lankaoaizao 8[J]. Acta Agronomica Sinica，35（1）：124-131.]

史曉芳，仇松英，史忠良，謝福來，高煒，宋立紅. 2017. 播期和播量對冬小麥堯麥16群體性狀和產量的影響[J]. 麥類作物學報，37（3）：357-365. [Shi X F，Qiu S Y，Shi Z L，Xie F L，Gao W，Song L H. 2017. Effect of sowing date and sowing amount on population traits and yield of winter wheat cultivar Yaomai 16[J]. Journal of Triticeae Crop，37（3）：357-365.]

鄔定榮，歐陽竹，趙小敏，于強，羅毅. 2003. 作物生長模型WOFOST在華北平原的適用性研究[J]. 植物生態學報，27（5）： 594-602. [Wu D R，Ouyang Z，Zhao X M，Yu Q，Luo Y. 2003. The applicability research of WOFOST model in north China plain[J]. Acta Phytoecologica Sinica，27（5）： 594-602.]

徐德利，王艷，王兵，李國權，劉冬玲，許光輝. 2016. 基于“雪壓稻”災情對江蘇省淮河以北地區稻麥周年協調發展的思考[J]. 作物雜志，（3）：128-133. [Xu D L，Wang Y，Wang B，Li G Q，Liu D L，Xu G H. 2016. Ideas on coordinated development for annual production of rice and wheat in North of Huai River of Jiangsu region based on the disaster of “rice over whelmed by snow”[J]. Crops，（3）：128-133.]

楊文鈺，屠乃美. 2011. 作物栽培學各論（南方本）[M]. 第2版. 北京：中國農業出版社. [Yang W Y，Tu N M. 2011. Crop cultivation（South China edition）[M]. The 2nd Edition. Beijing： China Agriculture Press.]

張弘，劉偉昌，李樹巖. 2019. WOFOST模型對河南省冬小麥模擬的適用性分析[J]. 氣象與環境科學，42（1）：36-42. [Zhang H，Liu W C，Li S Y. 2019. Applicability analysis of WOFOST model for winter wheat in Henan[J]. Meteo-rological and Environmental Sciences，42（1）： 36-42.]

張建平，趙艷霞，王春乙，楊曉光，王靖. 2012. 不同發育期干旱對冬小麥灌漿和產量影響的模擬[J].中國生態農業學報，20（9）：1158-1165. [Zhang J P，Zhao Y X，Wang C Y，Yang X G，Wang J. 2012. Impact simulation of drought disaster at different developmental stages on winter wheat grain-filling and yield[J]. Chinese Journal of Eco-Agriculture，20（9）： 1158-1165.]

張建平，趙艷霞，王春乙，楊曉光，王靖. 2013. 基于WOFOST作物生長模型的冬小麥干旱影響評估技術[J]. 生態學報，33（6）：1762-1769. [Zhang J P，Zhao Y X，Wang C Y，Yang X G，Wang J. 2013. Evaluation technology on drought disaster to yields of winter wheat based on WOFOST crop growth model[J]. Acta Ecologica Sinica，33（6）：1762-1769.]

張凱，李巧珍，王潤元，郭林，雷俊，王艷斌，王鶴齡. 2012. 播期對春小麥生長發育及產量的影響[J]. 生態學雜志，31（2）：324-331. [Zhang K，Li Q Z，Wang R Y，Guo L，Lei J，Wang Y B，Wang H L. 2012. Effects of sowing data on the growth and yield of spring wheat[J]. Chinese Journal of Ecology，31（2）：324-331.]

張鐵楠，許為政，魏湜，顧萬榮，賀丹，蘆玉雙，錢誠，陳琦，李晶. 2015a. WOFOST模型對東北春麥區春小麥生長和產量的模擬效果[J]. 麥類作物學報，35（1）：120-128. [Zhang T N，Xu W Z，Wei S，Gu W R，He D，Lu Y S，Qian C，Chen Q，Li J. 2015a. Simulation and verification of growth and yield of spring wheat in Northeast Region based on WOFOST model[J]. Journal of Triticeae Crops，35（1）：120-128.]

張鐵楠，許為政，魏湜，顧萬榮，賀丹，蘆玉雙，錢誠，陳琦，李晶. 2015b. 基于WOFOST模型的東北地區春小麥水分平衡模擬與驗證[J]. 水土保持研究，22（4）：45-51. [Zhang T N，Xu W Z，Wei S，Gu W R，He D，Lu Y S，Qian C，Chen Q，Li J. 2015b. Simulation and verification of spring wheat water balance in Northeast Region based on WOFOST model[J]. Research of Soil and Water Conservation，22（4）：45-51.]

張雪芬，余衛東，王春乙，白凌霞. 2006. WOFOST模型在冬小麥晚霜凍害評估中的應用[J]. 自然災害學報，15（12）：337-341. [Zhang X F，Yu W D，Wang C Y，Bai L X. 2006. Application of WOFOST model to assessment of winter wheats chilling damage by late frost[J]. Journal of Natural Disasters，15（12）：337-341.]

張耀蘭，曹承富，喬玉強，杜世州，趙竹，李華偉. 2012. 超晚播條件下密度和追氮時期對淮北地區小麥產量及品質的影響[J]. 中國農學通報，28（36）：160-164. [Zhang Y L，Cao C F，Qiao Y Q，Du S Z，Zhao Z，Li H W. 2012. Effect of density and topdressing time of nitrogen on yield and quality of wheat in Huaibei area under super-late so-wing[J]. Chinese Agricultural Science Bulletin，28（36）：160-164.]