Forecast of Maize Production in Henan Province

Considering the influence of meteorological factors on maize production, in order to improve the yield of maize in Henan Province, a grey combination model is constructed to predict the yield of maize in Henan Province. Firstly, the yield of maize in 2017 is obtained by GM (1, 1) model; secondly, the trend yield of maize is obtained by HP filter method, then the meteorological yield of maize is obtained, and the yield of maize reduction is determined according to the meteorological yield. Combined with Markov model, the maize yield reduction in various cities in Henan Province is forecasted. Finally, based on the reduction of production, policy recommendations are made for maize production in Henan Province.


Introduction
Henan Province is located in the middle and lower reaches of the Yellow River and the upper reaches of the Huaihe River.The land area is 167,000 square kilometers, accounting for 17.3 percent of the country's total area.Henan is a major grain producing Province, the stability of food production is related to people's basic needs.There are many factors affecting grain production: agricultural production technology, meteorological factors and policy price factors.Among them, as one of the main factors affecting grain production, the amount of meteorological factors investment is not controlled by humans.For meteorological disasters, all that people can do is to predict and prevent them in order to reduce their impact on grain production.The per unit area yield of maize in Henan Province was 5638 kg/hm 2 in 2012 and fell to 5285 kg/hm 2 in 2016 [1].
With the development of science and technology, the effect of agricultural production and policy price on maize yield is positive.Therefore, meteorological factors are the main factors that affect grain yield reduction.
The following scholars used different methods to study maize production, including: Ma Shuqing et al. [2] analyzed the correlation between maize seedling condition and maize yield, and provided a variety of evaluation models to evaluate the yield reduction of maize in spring drought.Sun Yulian et al. [3] explored the key meteorological factors affecting the growth of Linxia maize in Gansu Province and constructed a spatial-temporal dynamic climate prediction model for forecasting maize in different growth stages.Luo Zhengyuan et al. [4] used support vector machine (SVM) algorithm to predict maize yield in Heilongjiang Province.Zhang Wenying et al. [5] Li Guoqiang et al. [6] used grey GM (1, 1) model to predict maize yield in Hebei and Huanghuaihai.Jiang Hewen et al. [7] used a method of periodic correction of residual values to construct a grey prediction model to predict grain yield in Liaoning Province.Fang Huimin et al. [8] used BP neural network to predict the regional test yield.Su Hengqiang et al. [9] used the entropy weight method to determine the weight of a single model, and used combined forecasting model to predict maize yield.
Hou Yingyu, et al. [10] established a dynamic forecasting model to predict the yield of maize and winter wheat in the United States.Wei Jie et al. [11] used exponential smoothing method to predict maize yield in Shanxi Province.Ma Jianrui et al. [12] used water-fertilizer coupling model to predict maize yield.
In addition, some scholars have studied the relationship between meteorological change and grain production.Wang Futang [13] used three atmospheric circulation models to predict possible changes in rice, wheat and maize production in China's major crops in future climate scenarios.Using the WOFOST crop model and the future climate scenario data from the climate model BCC-T63, Zhang Jianping et al. [14] simulated and analyzed the effect of future climate change on maize yield in Northeast China.Jin Zhiqing et al. [15] and Xiong Wei [16] used CERES-Maize to simulate the impact of global climate change on maize production and agricultural production in ecologically fragile areas in China.Chen Shang et al. [17] combined historical meteorological data, K-NN algorithm and CERES-Maize model to construct a dynamic prediction method for maize growth period.Shang Zongbo [18] used the maize growth physiological and ecological simulation model (MPESM) to analyze the growth trend and yield variation of spring maize in Shenyang under the background of global climate change.Jiqi [19] combined multiple regression forecasting model with disaster reduction rate to predict maize yield and make the forecast of yield tend to be quantitative.
The above scholars predict maize yield single, or establish crop model to predict maize yield.Less prediction was made on maize yield reduction.The grey GM (1, 1) model works well for short-term prediction.The HP filter method is reasonable for separating crop trend yield.In addition, the Markov model is suitable for the prediction problem of data series with large random fluctuations.
Therefore, after using the GM (1, 1) model to obtain the corn yield in Henan Province in 2017, the HP filter method is used to obtain the trend yield of maize, and then the meteorological yield of maize is obtained.Finally, markov model is used to predict the yield of maize in Henan Province.And we provided policy recommendations for future maize production in Henan Province.

Basis of Theory and Method
In 1982, Professor Deng Julong published "The Control Problems of Grey systems" [20], which marked the birth of the emerging discipline of the grey system theory [21].After more than 30 years of development, the grey system theory has been applied to the fields of agriculture, economy, society and industry, solved many problems in production and life.

Grey GM (1, 1) Model
Grey GM (1, 1) model is mainly used in the prediction of grey system theory in prediction.Grey GM (1, 1) model requires less data, can be tested, and has a high precision in short-term prediction.It has been widely used in many scientific fields such as industry, agriculture, energy, transportation, geology, meteorology, hydrology, ecology, environment, medicine, military, economy, society and so on, and has successfully solved a large number of practical problems in production, life and scientific research.The modeling mechanism is as follows: Let an original sequence be: Of which: ( ) ( ) Perform an accumulation (1-AGO) generation on the sequence to get the sequence: ( ) 1 X adjacent mean generating sequence is: the mean form of the GM (1, 1) model.
Let the model's parameter list be the whitening differential equation of ( ) ( ) ( ) ( ) Then the solution of Equation ( 6) is Then the simulated prediction difference form of the original sequence is In this model, a is the development coefficient and b is the grey effect.

Markov Model
The research object of Markov Probability Matrix Prediction model is a stochastic dynamic system.Markov model is effective in predicting stochastic volatility problems.It predicts the development of the system according to the transition probability between states.The transition probability reflects the degree of influence of each random factor and the inherent regularity of each state.Markov Probability Matrix is suitable for predicting data columns with large random fluctuations [22].
The Markov forecasting model [23] is as follows: ∈ ,which denote the probability that the system is in the state j at time m k + under the condition that the system in the state i of at the moment m.That is, the transition probability of the k step is experienced; by sequentially sorting ( ) k ij P , the following matrix can be obtained.
Then the matrix is the k-step transition probability matrix of Markov chain.
Of which: ( ) the state transition probability.Since the final state of the data sequence is not clear, the last k data in the data sequence should be removed when calculating i M .

HP Filter Method
The HP filtering method is a decomposition method of time series in state space.It assumes that the time series consists of two parts, namely: long-term trend component and short-term fluctuation component.The HP filter method can be considered as a high-pass filter that can separate high-frequency components with periods below 8a [24].
Maize yield sequence { } t y is divided into trend yield t g and fluctuating yield t c by HP filter, in which t g is a stable high frequency trend signal and t c is an unstable low frequency disturbance signal.The principle is to remove the low frequency signal t c by filtering, thereby detecting the high frequency signal t g from the production sequence.HP filter decomposes t y into: The detection method of high frequency signal t g is often defined as the so- lution of the following minimization formula, namely: In the formula, λ is the HP filter parameter.When

Calculation of Meteorological Production
Crop yield is usually decomposed into three parts: trend yield, meteorological yield and random error.The trend output reflects the long-period production component of the level of productivity development in the historical period, also known as technical output.Meteorological production is a volatility that is affected by short-period changes in meteorological factors [25].Trend yields often use various linear or nonlinear simulations with "time" as the independent variable.The general formula can be expressed as follows: where Y represents the actual yield of maize, t Y is the trend yield, which represents the contribution of social and economic conditions and technological level of productivity to grain production.w Y is meteorological production, which characterizes the contribution of climate fluctuations to actual production.
e is the component of yield affected by some random factors, which accounts for a small proportion and is often neglected in practical calculation.In general, the level of agricultural technology is gradually improved, so trend yield is usually a function of time, and its influence on yield is a more smooth process in time series, especially in a larger geographical area [26].In this paper, the HP filter me-thod is used to obtain the trend yield of each year's maize, and then the meteorological yield of each year is calculated by w t Y Y Y = − .

Data Collection and Collation
The  2.
For the year when the meteorological yield is less than 0, we set the year of production reduction.For the year when the meteorological yield is greater than 0, we set the year of production increase.Therefore, according to the Markov prediction model, the reduction of yield in various places can be predicted.Take Zhengzhou City as an example, set a reduction of production to 1 M , and in- crease production to 2 M , therefore ( ) ( ) ( ) ( ) The state transition probability is The state transition probability matrix is ( ) As Zhengzhou maize is in the state of reduced production in 2017, according to the maximum probability criterion, it can be obtained Therefore, it can be found that Zhengzhou 2018 maize will increase production.
The same principle is available that the state transition probability matrix of Kaifeng City, Luoyang City, Pingdingshan City, Anyang City, Hebi City, Xinxiang City, Jiaozuo City, Fuyang City, Xuchang City, Luohe City, Sanmenxia City, Nanyang City, Shangqiu City, Xinyang City, Zhoukou City, Zhumadian City and Jiyuan City are as follows:
The geographical location of 18 cities in Henan Province is shown in Figure 1.By observing the map and combining the results of Markov forecast, the following conclusions can be obtained: Around Anyang are Xinxiang, Puyang, Hebi, Anyang including Hebi.According to Markov forecast, Xinxiang will increase production in 2018, while Puyang will reduce production and Hebi is uncertain, but Anyang obviously has more border areas with Xinxiang.As a result, Anyang may increase its maize production in 2018.Hebi is surrounded by Anyang and Xinxiang, Anyang and Xinxiang are both increasing production, so it is predicted that Hebi' maize will also increase production in 2018.Jiaozuo is surrounded by Xinxiang, Luoyang, Zhengzhou and Jiyuan, and mainly Xinxiang, Luoyang and Zhengzhou.In these three cities, maize in 2018 is predicted to increase production, so it is predicted that Jiaozuo maize may increase production.Around Jiyuan is Luoyang and Jiaozuo.As the two cities' maize production predictions may increase production, it is predicted that Jiyuan' maize may increase production in 2018.Around Sanmenxia, it is Luoyang and Nanyang.Both cities have increased production.Therefore, it is estimated that Sanmenxia' maize may increase production in 2018.Around Zhoukou, there are Shangqiu, Kaifeng, Xuchang, Luohe and Zhumadian.Except Zhumadian, all other cities are expected to increase production.From the perspective of contact area, it is predicted that Zhoukou's 2018 maize production may also increase production.According to the forecast, it is found that most areas of Henan Province increased production in 2018, which is more in line with the actual situation.Although the meteorological anomalies will make maize yield decrease to a certain extent, with the improvement of scientific and technological level, People take precautions against meteorological anomalies in advance and compensate the yield losses caused by meteorological anomalies to a certain extent, thus increasing maize production in Henan Province on the whole.

Policy Recommendations
Through improving the varieties of maize, selecting and cultivating new varieties with strong resistance, and adjusting the production structure of grain, it can improve its ability to defend against meteorological anomalies, reduce the losses caused by agricultural meteorological anomalies, and ensure stable production and increase of grain production.
We should vigorously strengthen the construction of farmland and water conservancy and do a good job in the construction of agricultural infrastructure projects such as well irrigation and water-saving irrigation.It is necessary to promote the comprehensive utilization of straw returning to the field, protect the ecological environment, avoid over-exploitation and utilization of land resources, and improve the quality of farmland.
Strengthen research on the prevention of agrometeorological disasters.It is necessary to increase the systematic analysis of disasters and explore its complex background to develop detailed supporting defense measures.Through the research on the disaster early warning system, the ability to predict and forecast disasters will be improved.Conduct in-depth research on meteorological anomalies and other disciplines to improve the response of maize to meteorological anomalies and reduce the impact of meteorological anomalies on maize yield.

Conclusion
In order to improve grain production and ensure national grain security, considering the influence of meteorological factors on grain production, the grain 0 λ = , the minimized solution of the function is { } t y sequence.With the increase of λ value, the trend of minimum solution estimation is smoother.The general experience of λ is as follows: Table 1.Yield per unit area of maize in 18 cities of Henan Province in 2006-2017.

Figure 1 .
Figure 1.Geographic map of 18 cities in Henan Province.
is the number of original data samples transferred from state i ⊗ yield data set of Henan maize in 2006-2016 comes from Henan Statistics Bureau [1].The research objects selected in this paper are 18 cities in Henan Province.In addition, the yield data of Henan Province in 2017 are obtained by GM (1, 1) model based on the data of the previous 10 years.The detailed data are shown in Table1.The trend yield of maize in 18 cities of Henan Province from 2006 to 2016 is obtained by using HP filter method.The calculated results are shown in Table