# Simulation of Winter Wheat Yield with WOFOST in County Scale

Shangjie Ma, Zhiyuan Pei, Yajuan He, Lianlin Wang, Zhiping Ma

► To cite this version:

Shangjie Ma, Zhiyuan Pei, Yajuan He, Lianlin Wang, Zhiping Ma. Simulation of Winter Wheat Yield with WOFOST in County Scale. 9th International Conference on Computer and Computing Technologies in Agriculture (CCTA), Sep 2015, Beijing, China. pp.161-172, 10.1007/978-3-319-48357-3_16. hal-01557810

HAL Id: hal-01557810

https://inria.hal.science/hal-01557810

Submitted on 6 Jul 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Simulation of Winter Wheat Yield with WOFOST in County Scale

Shangjie Ma¹,ª, Zhiyuan Pei¹,ᵇ,*, Yajuan He¹, Lianlin Wang², Zhiping Ma²

¹ Chinese Academy of Agricultural Engineering, Beijing 100125, China;

² Yutian Agriculture Planning Office, in Hebei Province, Yutian 064100, China.

ªmashangjie@agri.gov.cn, ᵇpeizhiyuan@agri.gov.cn

**Abstract.** Winter wheat is mainly planted in water shortage area, such as North China and Northwest China. As a key field management measure, irrigation plays an important role in the production of winter wheat. This paper focuses on the improvement of regional winter wheat yield estimation technique in county scale by adjusting the irrigation management measure in crop growth model. The WOFOST (World Food Study) model was used by dividing the whole county into a number of EMUs element mapping units) and then running the model in each unit in sequence. While running, the measured soil moisture and LAI were used to rate the irrigation parameters. Finally, the calibrated irrigation parameters were used to run the model again. The results showed that the simulated winter wheat growth process was normal. During the whole growing period of winter wheat, the change trends of the time series of soil moisture and LAI were basically consistent with that of the measured. The precision of simulated yield was between 87.26% and 98.68% among the 5 units, and the average of the precision was 94.56%. The precision of simulated winter wheat yield was well, and could meet the needs of winter wheat yield estimation in county-wide. This study may provide basis for estimating crop yield in regional area by using the crop growth model.

**Keywords:** WOFOST, Irrigation Parameter, Winter Wheat Yield

## 1 Introduction

Winter wheat, one of the main crops in northern China, is mainly planted in water shortage area, such as North China and Northwest China. As a key management measure, irrigation plays an important role in the production of winter wheat. Driving by the light, temperature, water, soil and other environmental variables, the crop growth dynamic model can daily simulate the photosynthesis, respiration, transpiration and evaporation process in crop growth period, and their relationship with meteorological, soil, and other environmental conditions, and can quantitatively describe the dynamic process of crop growth, development, and yield information. WOFOST model originated in the framework of interdisciplinary studies on world food security and on the potential world food production by the Center for World Food Studies in cooperation with the Wageningen Agricultural University[1]. WOFOST is a mechanism model, and had been successfully used in daily business of agro meteorological monitoring and yield forecasting in European Union[2].
Meanwhile, WOFOST was also widely used in board as follows. Wu Dingrong[3] evaluated the applicability of WOFOST in the North China Plain. Xie Wenxia[4] simulated and validated the rice potential growth process in Zhejiang Province by utilizing WOFOST. According to climate characteristic and ecologic type of winter wheat variety in North China, Ma Yuping[1] adjusted the parameters of WOFOST, and simulated phenological phase of overwinter in the region. Chen Zhenlin[5] simulated the integrated impacts of low temperature and drought on maize yield. The abovementioned applications were in field scale. Combining the GIS technology, the application of WOFOST can be extended to regional scale. Ma Yuping[6] simulated the winter wheat growth in regional scale in North China Plain by scaling-up WOFOST model. Li Xiufen[7] evaluated the suitability of meteorological conditions during growing period of Maize in northeast of China. Li Xiufen[8], Gao Yonggang[9-10], Du Chunying[11] respectively simulated the growth process of rice, wheat, maize, soybean, beet, potato and so on using WOFOST in Heilongjiang Province. Huang Jianxi[12] estimated the winter wheat yield in region by the method of optimizing the parameters of WOFOST based on assimilating simulated LAI and MODIS LAI at the re-greening stage of winter wheat. However, there were few reports about the method of simultaneously using the measured LAI and soil moisture to rate the irrigation parameters of model, and to simulate crop growth processing in regional scale.

In this paper, Yutian County, a water shortage region in Hebei Province, was selected as the research area. Both the measured LAI and soil moisture were used to rate the irrigation parameters. The adjusted irrigation parameters were used to run the model to simulate the winter wheat growth processing and to estimate the yield of winter wheat in county scale. All these would provide the basis for regional crop yield estimation by using crop growth model.

# 2 MATERIALS AND METHODS

## 2.1 STUDY AREA

Yutian county locates between 117° 31′ E～117° 56′ E and 39° 30′ N～39° 58′ N, at the south of Yan Mountains, the center of Eastern Hebei Plain. The terrain of the county slopes gently from northeast to southwest, the altitude in north is about 50～400 meters, in center and south under 25 meters, forming the views of the northern hills, central plains, the southern low-lying land pattern. According to the second national soil survey, there are 4 soil types, 9 subtypes in the county. The total area of the main 5 subtypes, leached cinnamon soil, meadow cinnamon soil, swamp alluvial soil, alluvial soil, salted alluvial soil, accounts for nearly 95%. The county belongs to warm temperate semi humid region in the eastern monsoon region, with the average annual rainfall 655.4 mm, the average temperature 11.6 degrees, and the annual frost free period 196 days. The main crops include winter wheat, spring maize,
summer maize. The main varieties of winter wheat are Jingdong NO.8 and Lunxuan NO. 987, with average yield level 3750~6000 kg/hm².

<figure><img src="image_2.png" /></figure>

Legend
<caption>Fig. 1. Spatial distribution of experiment sites and soil types in study area</caption>

## 2.2 MEASURED DATA

According to the environmental condition and the spatial distribution of winter wheat, 36 sites were chosen to measure ground parameters in 2012. In order to ensure the purity of the remote sensing image pixels, the area of winter wheat must be not less than 200 m×200 m in each site. The distribution of the sites, shown in Fig.1, should be representative and homogeneous. Based on the phenological phase of winter wheat, some key stages from re-greening to mature were selected to measure soil moisture, LAI and other crop parameters. The three elements of yield, including the number of heads per hectare, grain number per head, and 1000 grain weight, were collected to estimate the yield from heading stage to mature stage. Dates were respectively April 9-12 re-greening stage, May 2-6 jointing stage, May 16-19 heading stage, June 4-6 grain filling stage, and June 13-17 mature stage. Only soil moisture was measured in grain filling stage, and only crop parameters were measured mature stage.

Method of measuring LAI is as follows. Three subareas about 20 m² for each one were homogeneously chosen in each site. Five values of LAI from different locations in each subarea were measured by LAI2000 which was produced by U.S. LI-COR Company. The average of the five values is read as the LAI value of the subarea, and the average of three subareas is written as that of the site. The measurement is done between 09:00 am and 17:00 pm.

Method of measuring soil moisture is as follows. Fifteen representative sites whose area of winter wheat beyond 500m×500m were selected as soil moisture measuring sites from the LAI measuring sites. The soil volume moisture of winter wheat field in the four angles of the site was measured respectively at the depth of 0~10、10~20、20~40、40~60cm by the TSC-I soil moisture measuring instrument which was produced by Chinese Agricultural University.

The winter wheat yield is calculated by the method called three agricultural elements estimating yield algorithm, which is gained from the number of heads per hectare, grain number per head, and 1000 grain weight multiplication[13-14]. The detailed steps are as follows. In each site, three 1m² rectangular areas were selected randomly. The number of heads in the three rectangular areas during heading stage and grain filling stage are counted respectively. The average of heads in the three rectangles at the same stage is used to calculate the number of heads per hectare at that stage. The average of the number of heads per hectare of the two stages is the number of heads per hectare of the site. At mature stage, total 50~80 heads wheat were randomly selected as the sample from the above three rectangular areas in each site. After drying, threshing, counting, weighing, grain number per head, and 1000 grain weight were calculated.

## 2.3 STUDY IDEAS

Based on the spatial discretization and irrigation parameters rating, the purpose of estimating winter wheat yield in county scale by using the WOFOST model was
implemented. The steps were shown as Fig. 2 in follows. First, overlaying the county
boundary and the map of soil type spatial distribution, dividing the whole county into
several EMUs in which the combination of meteorology and soil type was
homogeneous. Then, the parameters of weather, soil, crop, field management
measures including initial irrigation management measures extracted from local
management habits were prepared in each unit. After that, the model was run. The
ground sites in each unit were selected, and the measured LAI and soil moisture in
these sites were averaged as that of the unit. With the help of measured LAI and soil
moisture, the initial irrigation operations were calibrated. Finally, the calibrated
irrigation parameters were used to run the model again.

<figure><img src="image_3.png" /><figcaption>Fig. 2. Flow chart of study</figcaption></figure>

2.4 WOFOST

MODEL LOCALIZATION WOFOST needs parameters of meteorology, soil, crop, management and so on. Both the daily and monthly meteorology parameters can be accepted by the model[15]. The format of daily meteorological parameter is CABO, which contains the maximum and minimum temperature(°C), precipitation(mm), average wind speed(m/s), water vapor pressure at 08:00 am(kPa), solar radiation (MJ/m2·d). The daily measured meteorological data from January 1st, 2010 to June 30th, 2012, including the maximum and minimum temperature, precipitation and
sunshine duration, was got from Yutian Meteorology Office(117° 44' E、39° 53' N, Altitude 14.4 m). The period covered two whole growth stages of winter wheat, 2010/2011, and 2011/2012. The measured data from neighboring Tangshan Meteorology Observation Station (118° 9' E、39° 40' N) from January 1st, 1957 to June 30th, 2012, including daily average atmosphere pressure, maximum atmosphere pressure, minimum atmosphere pressure, average temperature, maximum temperature, minimum temperature, average relative humidity, minimum relative humidity, precipitation, average wind speed, sunshine duration and so on, was also downloaded on the internet as additional information[16]. The wind speed parameter needed by the model could be replaced by the observed data from Tangshan. The other needed parameters, water vapor pressure and solar radiation, which had no observed data, could be calculated. The water vapor was calculated by relative humidity multiplying the saturation vapor pressure. The saturation vapor pressure was calculated by the formula recommended by China Meteorological Administration[17]. The solar radiation was calculated by sunshine duration, the maximum and minimum temperature[18-19]. Harmonized World Soil Database from Food and Agricultural Organization of the United Nations[20] and the data from the Second National Soil Survey were also collected. The map of soil spatial distribution was extracted from the soil database of FAO. The cinnamon soil distributed in non cultivated land area, the northern mountainous area, was eliminated firstly. The smaller proportion soil types in cultivated land area were combined to the adjacent bigger types. The main five types, leached cinnamon soil, meadow cinnamon soil, swamp alluvial soil, alluvial soil, salted alluvial soil, were finally extracted. The most physical and chemical properties referred to that in model, and the important parameters, such as field capacity, saturated water content, wilting coefficient, were calculated by SPAM (soil-plant-atmosphere-water) using the texture, organic matter content from the soil database of FAO[21]. Some crop parameters such as the temperature sum from sowing to emergence(TSUMEM), the temperature sum from emergence to anthesis(TSUM1), and the temperature sum from anthesis to maturity(TSUM2), was calculated by the long time-series observed temperature data in Tangshan and the surveyed long-time phenology of winter wheat in local. The other parameters for crop were confirmed by looking up papers and documents[1,3,15,22-23] or debugging the model. The detailed crop parameters were listed in table1 as below. The needed filed management parameters were gained by surveying the planting and managing habit of winter wheat in local. In general, the winter wheat is sowing in early October, harvesting in the second or third ten-days in June. So, when running, the variable sowing date was selected, with the earliest Oct 1st and the latest Oct 10th, the harvesting date was calculated automatically by the model and the longest duration was less than 300 days.

<table><caption>Table 1. Crop parameters settings for WOFOST model</caption><thead><tr><td>Parameters</td><td>Description</td><td>Value</td></tr></thead><tbody><tr><td>TBASEM</td><td>lower threshold temperature for emergence</td><td>0</td></tr><tr><td>TEFFMX</td><td>maximum effective temperature for emergence</td><td>30.0</td></tr><tr><td>TSUMEM</td><td>temperature sum from sowing to emergence</td><td>113.0</td></tr><tr><td>TSUM1</td><td>temperature sum from emergence to anthesis</td><td>1150.4</td></tr></tbody></table>
<table><tr><td>TSUM2</td><td>temperature sum from anthesis to maturity</td><td>900.4</td></tr><tr><td>LAIEM</td><td>leaf area at emergence</td><td>0.1365</td></tr><tr><td>RGRLAI</td><td>maximum relative increase in LAI</td><td>0.00917</td></tr><tr><td>SPAN</td><td>life span of leaves growing at 35°C</td><td>31.3</td></tr><tr><td>TBASE</td><td>lower threshold temperature for aging of leaves</td><td>0</td></tr><tr><td>EFFTB</td><td>light-use efficiency of single leaf</td><td>0.45</td></tr><tr><td>CVL</td><td>efficiency of conversion into leaves</td><td>0.685</td></tr><tr><td>CVO</td><td>efficiency of conversion into storage organs</td><td>0.709</td></tr><tr><td>CVR</td><td>efficiency of conversion into roots</td><td>0.694</td></tr><tr><td>CVS</td><td>efficiency of conversion into stems</td><td>0.662</td></tr><tr><td>Q10</td><td>relative increase in respiration rate per 10°C temperature increase</td><td>2</td></tr><tr><td>RML</td><td>relative maintenance respiration rate of leaves</td><td>0.03</td></tr><tr><td>RMO</td><td>relative maintenance respiration rate of storage organs</td><td>0.01</td></tr><tr><td>RMR</td><td>relative maintenance respiration rate of roots</td><td>0.015</td></tr><tr><td>RMS</td><td>relative maintenance respiration rate of stems</td><td>0.015</td></tr><tr><td>TBASEM</td><td>maximum relative death rate of leaves due to water stress</td><td>0.03</td></tr></table>

**RATING IRRIGATION PARAMETERS** The irrigation parameters, including irrigation date and irrigation volume, were initially set on the base of surveying local irrigation condition. In general, there is 4 times irrigation during the winter wheat growth period. The volume of irrigation for each time is about 50~100 mm. The first time of irrigation often occurs before overwintering period in early December. The second time of irrigation will be implemented after re-greening period in middle ten-days in March next year. The third irrigation may occur before jointing stage, in middle or late April. The fourth time irrigation will be done in heading-anthesis stage in middle May. The times of irrigation may be reduced to three according to the actual weather condition. In the paper, the irrigation parameters were initiated four times, 70mm for each time, and the dates of which were respectively in November 30th, March 14th the next year, April 14th, and May 14th.

The assumed irrigation parameters were used to run the model firstly. While running, the parameters were rated with the measured LAI and soil moisture by using the trial-and-error method. The initial irrigation volume as the center, 10mm as the step size, irrigation volume was adjusted between 50~100 mm. The initial irrigation date was determined as the center, and 1 day as the step size. The irrigation date was changed from -5 to 5 days, and avoided the rainy day. The parameters which make the error of simulated and measured LAI and soil moisture minimum and the irrigation amount lowest were chosen as the best value. Overlaying the map of soil type distribution and boundary, the whole county was divided into 5 EMUs, which were respectively leached cinnamon soil unit, meadow cinnamon soil unit, swamp alluvial soil unit, alluvial soil unit, and salted alluvial soil unit. The irrigation parameters for each unit were shown as table2.

<table><caption>**Table 2. Calibrated irrigation parameters**</caption><thead><tr><td>Soil Type</td><td>First Irrigation</td><td>Second Irrigation</td><td>Third Irrigation</td><td>The 4th Irrigation</td></tr></thead></table>
<table><thead><tr><th></th><th>Date</th><th>Volume (mm)</th><th>Date</th><th>Volume (mm)</th><th>Date</th><th>Volume (mm)</th><th>Date</th><th>Volume (mm)</th></tr></thead><tbody><tr><td>Leached</td><td>11-</td><td>70</td><td>12-3-</td><td>80</td><td>12-4-</td><td>80</td><td>12-5-</td><td>70</td></tr><tr><td>Cinnamon Soil</td><td>11-30</td><td></td><td>19</td><td></td><td>18</td><td></td><td>17</td><td></td></tr><tr><td>Meadow</td><td>11-</td><td>70</td><td>12-3-</td><td>70</td><td>12-4-</td><td>80</td><td>12-5-</td><td>60</td></tr><tr><td>Cinnamon Soil</td><td>11-30</td><td></td><td>19</td><td></td><td>18</td><td></td><td>17</td><td></td></tr><tr><td>Swamp</td><td>11-</td><td>70</td><td>12-3-</td><td>70</td><td>12-4-</td><td>70</td><td>12-5-</td><td>70</td></tr><tr><td>Alluvial Soil</td><td>11-30</td><td></td><td>14</td><td></td><td>10</td><td></td><td>13</td><td></td></tr><tr><td>Alluvial Soil</td><td>11-30</td><td>70</td><td>12-3-</td><td>70</td><td>12-4-</td><td>70</td><td>12-5-</td><td>70</td></tr><tr><td></td><td></td><td></td><td>14</td><td></td><td>9</td><td></td><td>13</td><td></td></tr><tr><td>Salted</td><td>11-</td><td>70</td><td>12-3-</td><td>70</td><td>12-4-</td><td>50</td><td></td><td></td></tr><tr><td>Alluvial Soil</td><td>11-30</td><td></td><td>19</td><td></td><td>18</td><td></td><td></td><td></td></tr></tbody></table>

# 3 RESULTS AND ANALYSIS

The WOFOST model was run in each unit in sequence. Both the no irrigation and rated irrigation were run respectively to simulate soil moisture, LAI, and winter wheat yield.

## 3.1 SOIL MOISTURE SIMULATION

The simulation results shown as in Fig 3 showed that the precipitation cannot meet the needs of winter wheat. The soil moisture kept at low level during the whole growth stage. After irrigation, the simulated soil moisture was similar with that of measured. The peak of time-series curve could reflect the irrigation. In general, the average percentage error of the simulated and measured among twenty values, which were measured during the four field experiments for all five soil types, was -17.56%. For some values, the error of absolute value between the simulated and measured was large. But during the whole growth stage, the trends of simulated and measured soil moisture were basically similar.
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Comparison of simulated and measured soil moisture under different irrigation methods</figcaption></figure>

## 3.2 LAI SIMULATION

Under no irrigation, the results of the LAI simulation (shown as Fig4) showed that it is basically normal for the simulated LAI before heading stage of winter wheat. But, after heading stage, the LAI decreased quickly, leading the winter wheat died too early. After irrigation, the change trends of simulated and measured LAI were same. But in the early re-greening stage and lately maturity stage, there were relatively larger deviations between them. In the re-greening stage, the value of simulated was higher than that of measured, which might be caused by ignoring the influence of 'winter loss', a phenomenon that some leaves of winter wheat may be died in overwintering period. The simulated leaf might grow more quickly than the actual leaf in re-greening stage which affected by 'winter loss'. In the lately maturity stage, the simulated LAI was lower than that of measured. The measured LAI might be wrongly added by the dead leaves whose influence were difficulty to be eliminated while measuring. The maximum LAI values in the five units were respectively 4.64, 4.64, 4.46, 4.81, 3.69, and the corresponding simulated value were respectively 4.57, 4.57, 4.59, 4.57, 4.57. The RMSE(root-mean-square error) of the simulated and measured maximum LAI was 0.41, occurring 1.86% of the average maximum LAI. The date of the simulated maximum LAI was earlier about 10 days than that of the measured. This difference between them might be influenced by the different time step size of them. LAI would be measured only during the key stage, and the average time interval is about 20 days. But the model was simulated day by day. The different
temporal resolution between the simulated and measured LAI might lead the
difference of date.

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Comparison of simulated and measured winter wheat LAI under different irrigation methods</figcaption></figure>

3.3 YIELD SIMULATION

While no irrigation, the winter wheat died too early, and the yield was too little.
Added irrigation in the model, the simulated growth process was normal, and the
yield was also reasonable. The precision equals 1- abs(percentage error), which was
used to express the effect of simulation. The measured yields were respectively
4614.29, 4852.65, 5537.62, 5170.72, 4302.24 kg/hm² in the five EMUs, and that of
simulated were respectively 4943, 4857, 4832, 4970, 4359 kg/hm². The precisions for
the five EMUs were respectively 92.88%, 97.85%, 87.26%, 96.12%, 98.68%.
Averaging all the units, the measured yield was 4895.5 kg/hm², the simulated yield
was 4812.20 kg/hm², and the precision was 94.56%. The difference of simulated
yields among the EMUs was obviously lower than that of observed. For example, the
minimum measured yield was 4302.24 kg/hm², the maximum measured yield was
5537.62 kg/hm², the difference between them was 1235.38 kg/hm². Simultaneously,
the simulated value were respectively 4359 kg/hm², 4970 kg/hm², 611 kg/hm². It had
a lot to do with the considered parameters in model and the method of model
localization. The considered parameters in model included meteorology parameters,
soil parameters, crop parameters and irrigation parameters. But in the farm, the factors
that influencing yield were more complex, not only containing the above mentioned factors, but also soil fertility condition, local climate, field terrain, plant diseases and insect pests and so on. During the localization process, being limited by the data acquisition conditions, only soil and irrigation parameters were considered in different EMUs. But the crop and meteorology parameters were same in the whole county. It was not consistent with the actual situation. For example, the precipitation may be different in the whole county, in the region of 55km from north to south, 36km from east to west. But there is only one weather observing station in the whole county. Using the same meteorology parameters in all EMUs reduce the difference of simulated yield among the EMUs. But totally, the precision of simulation can meet the estimating yield in county scale.

<table><caption>Table 3. Comparison of simulated and measured winter wheat yield</caption><thead><tr><th>Soil Type</th><th>Measured Yield (kg/ha)</th><th>Simulated Yield No Irrigation (kg/ha)</th><th>Precision</th><th>Simulated Yield Calibrated Irrigation (kg/ha)</th><th>Precision</th></tr></thead><tbody><tr><td>Leached Cinnamon Soil</td><td>4614.29</td><td>1954</td><td>42.35%</td><td>4943</td><td>92.88%</td></tr><tr><td>Meadow Cinnamon Soil</td><td>4852.65</td><td>2352</td><td>48.47%</td><td>4957</td><td>97.85%</td></tr><tr><td>Swamp Alluvial Soil</td><td>5537.62</td><td>2390</td><td>43.16%</td><td>4832</td><td>87.26%</td></tr><tr><td>Alluvial Soil</td><td>5170.72</td><td>3413</td><td>66.01%</td><td>4970</td><td>96.12%</td></tr><tr><td>Salted Alluvial Soil</td><td>4302.24</td><td>1954</td><td>45.42%</td><td>4359</td><td>98.68%</td></tr><tr><td>Average</td><td>4895.5</td><td>2412.60</td><td>49.08%</td><td>4812.20</td><td>94.56%</td></tr></tbody></table>

## 4 CONCLUSION

As an example, Yutian County in Hebei Province was selected to divide into five EMUs, in which the combination of soil type and meteorology was same. The assumed irrigation parameters from surveying local irrigation condition were used to run the model firstly. Then, run the model, rated the irrigation parameters with the help of measured LAI and soil moisture by using the trial-and-error method. Finally, the rated irrigation parameters were used to run the model again, and completing the estimating yield in county scale. The conclusions were as follows.

1) The simulated growth process of winter wheat was basically normal. The change trends of the time-series of simulated soil moisture and LAI were basically consistent with the measured in the whole growth period of winter wheat. In the numerical, the difference between the simulated LAI and the measured LAI in the early re-greening stage and lately maturity stage was greater than that of in jointing stage and heading-
anthesis stage. The simulated maximum LAI values in the whole stage in all EMUs
were close to that of measured. But the date had an error about 10 days which might
be affected by the different time scale between simulation and measure.

2) The simulated results could meet the needs of winter wheat yield estimation in county-wide. The precision of simulated yield was between 87.26% and 98.68% among the 5 units, and the average of the precision is 94.56%.

Because being limited by the difficulty in obtaining parameters in large range, and ignoring some affected factors such as soil fertility condition, local climate, field terrain, plant diseases and insect pests and so on, the simulated difference among the EMUs was obviously lower than that of observed. In the more fine application, in order to improve the precision, the more parameters should be considered in the model and the more detailed input parameters should be used in the simulation.

Acknowledgment

Funds for this research was provided by the Research on the earth observation
system by high resolution remote sensing, and the cooperative application
technology of high resolution and multi sensors in agriculture(GF13/15-311-003),
and was supported by the Innovation Team of Crop Monitoring by Remote
Sensing(CMIT), in Chinese Academy of Agricultural Engineering(CAAE).

References

1. Yuping Ma, Shili Wang, Li Zhang. Study on improvement of WOFOST against overwinter of wheat in North China[J]. Chinese Journal of Agrometeorology, 2005,26(3):145-149.

2. Supit I,Hooijper A A,van Diepen C A, et al. System Description of WOFOST6.0 Crop Simulation Model Implemented in CGMS[M], Theory and algorithms. Wagenningen, Netherlands: The Winand Starting Center for Intergrated Land, Soil and Water Research(SC-DLO), 1994:1-144.

3. Dingrong Wu, Yangzhu Ou, Xiaomin Zhao et al. The applicability research of WOFOST model in North China plain[J]. Acta Phytoecologica Sinica, 2003, 27(5):594-602.

4. Xie Wenxia, Yan Lijiao, Wang Guanghuo. Simulation and validation of rice potential growth process in Zhejiang by utilizing WOFOST model[J]. Chinese Journal of Rice Science, 2006,20(3):319-323.

5. Zhenlin Chen, Jianping Zhang, Chunyi Wang 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

6. Yuping Ma, Shili Wang, Li Zhang et al. A preliminary study on a regional growth simulation model of winter wheat in North China based on scaling-up approach I potential production level[J]. Acta Agronomica Sinica, 2005,31(6):697-705.

7. Xiufen Li, Shuqing Ma, Lijuan Gong et al. Evaluation of meteorological suitability degree during maize growth period based on WOFOST in Northeast China[J]. Chinese Journal of Agrometeorology, 2013, 34(1): 43-49
8. Xiufen Li, Yuguang Wang, Shengtai Ji et al. Validation of crop growth monitoring system(CGMS) in Heilongjiang province[J]. Chinese Journal of Agrometeorology, 2005, 26(3):155-157.

9. Yonggang Gao, Jihai Na, Hong Gu et al. Characteristics and compartment of potential climatic productivity of potato(Solanum tuberosum L.) in Heilongjiang province[J]. Chinese Journal of Agrometeorology, 2007, 28(3):275-280.

10. Yonggang Gao, Rui Nan, Hong Gu et al. Simulation of beta vulgaris climatic productivity and its planting climatic zoning in Heilongjiang Province[J]. Chinese Journal of Ecology, 2009, 28(1):27-31.

11. Chunying Du, Jifen Li, Chenyi Wang et al. Study on dynamic yield forecasting of rice based on WOFOST model in Heilongjiang province[J]. Journal of Anhui Agri Sci 2011, 39 (24): 15093-15095, 15122

12. Jianxi Huang, Sijie Wu, Xingquan Liu 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 CASE), 2012, 28(4): 142-148.(in Chinese with English abstract)

13. Weicheng Zhang, Yu'e Wu, Yongguang Yang et al. Analysis of interactive relationship among yield combinational factors of winter wheat[J]. Journal of Henan Vocation-Technical Teachers College, 1992,20 (3): 1-4

14. Xiaohong Zhu, Kunqing Xie, Xiru Xu et al. The structure analysis of winter wheat yield and the principal of remote sensing estimation of winter wheat yield[J]. Remote Sensing of Environment China, 1989,4 (2):116-127

15. Boogaard H L, De Wit A J W, te Roller J A, et al. User's Guide for the WOFOST Control Center 1.8 and WOFOST 7.1.3 Crop Growth Simulation Model[M]. Wageningen: Wageningen University & Alterra Research Center, 2011. http://www.wofost.wur.nl

16. China Meteorological Data Sharing Service System[Z]. http://cdc.cma.gov.cn/home.do

17. Li Luo, Xiaolei Wang, Peng Yu. The compare and research of the calculate formula of the saturation water steam pressure[J]. Meteorological, Hydrological and Marine Instruments, 2003 (4): 24-27.

18. Fei Cao, Shuanghe Shen. Estimation daily solar radiation in China[J]. Journal of Nanjing Institute of Meteorology, 2008,31(4):587-591.

19. Fei Cao. Study on Earth Surface Total Solar Radiation in China in Recent 40 Years and Daily Solar Radiation Model[D]. Nanjing: Nanjing Institute of Meteorology, 2008.

20. Harmonized World Soil Database[Z]. http://www.fao.org/nr/land/soils/harmonized-world-soil-database/download-data-only/en/

21. Huaibin Wei, Zhanpang Zhang, Jinpeng Yang. Establishing method for soil database of SWAT model[J]. 2008,38(6):15-18.

22. Hu Zhao, Zhiyuan Pei, Shangjie Ma et al. Retrieving LAI by assimilating time series HJ CCD with WOFOST[J]. Transactions of the Chinese Society of Agricultural Engineering (Transaction of CASE), 2012,28(11):158-163.(in Chinese with English abstract)

23. Sijie Wu. Study on Winter Wheat Yield Prediction Based on Assimilating Remote Sensing Data and Crop Growth Model[D]. Changsha: Central South University, 2012.