# A Prediction Model of Forest Preliminary Precision Fertilization Based on Improved GRA-PSO-BP Neural Network

Authors: Chen Zuxing; Wang Dian
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The optimal amount of fertilizer application which was needed by the trees and the factors that influence the fertilization have an intricated nonlinear relationship. According to the problems that the traditional fertilization prediction model has, such as lacking of the scalability and practicality, this paper initiates an accurate fertilization prediction model that was based on the GRA-PSO-BP neural network which can make the accurate fertilization come true and improve the economic benefits of forest industry. This paper uses the GRA method to determine the input of the neural network as the site index and make the forest age, nutrient content of the advantage trees, biomass of the advantage trees, biomass of average trees, and target yield as the output numbers of the Actual amount of fertilizer applied. During the calculation process, the global particle swarm optimization algorithm is used to optimize the initial numbers and threshold numbers of BP neural network which build a phased GRA-PSO-BP accurate fertilization model. Compared with the prediction algorithm of full input variate that is based on the single BP neural network and the prediction algorithm of full input variate that is based on PSO-BP Neural Network, the GRA method can determine the key factors that influence the amount of fertilizer applied in different forest areas and modify the prediction model to improve the scalability and accuracy of the prediction and finally achieve the precision fertilization as the data of different forests updated, so we can see that the prediction result of this paper is more accurate. The result demonstrates that the GRA-PSO-BP neural network Segment fertilization model is more accurate than the traditional BP neural network and BP Neural Network that was optimized by the PSO algorithm, and specifically, the error of the predicted amount of fertilizer application and the actual amount of fertilizer application is less than 5%, which can effectively guide the fertilization in stages.


---

## Body

## 1. Introduction

In the process of forestry production, it is fundamentally required to calculate the amount of fertilizer, such as nitrogen, phosphorus, and potassium fertilizer, for different forest species and different soil conditions [1]. Cunninghamia lanceolata enjoys the characteristics of fast growth, wide application, high economic value, and easy reproduction [2]. Cunninghamia lanceolata is dependent on fertilizer and humidity. In the condition of 25°C–35°C, Cunninghamia lanceolata can grow rapidly. In this regard, fertilization is the key link for high yield of Chinese fir forest. Fertilization generally applies from young forest to middle-aged forest, which plays an important role in forestry economy. Forest fertilization functions as a critical technique to enhance soil fertility, improve tree nutrition, and facilitate rapid growth and high yield [3, 4]. At present, forestry production is increasingly transformed to precision operations. Cunninghamia lanceolata production is no exception. The premise of precision production of Chinese fir forest is precise fertilization since such an agrotechnique can help balance the nutrients in soil and increase the nutrients required by trees, so as to achieve the target yield. In addition, the implementation of precision fertilization can save the fertilizer and phase down environmental pollution.Fertilization prediction model, as the key of precision fertilization inCunninghoya lanceolata forest, has become a hot spot in forestry research. In response to this issue under discussion, Hu Yueli et al. proposed a site-nutrient effect fertilization model [5–9]. However, sometimes this model can hardly help to get satisfactory accuracy and stability. The traditional fertilization model can only describe the static relationship between the growth of forest and the required fertilization nutrients, which can hardly reflect the different growth periods of trees. In view of the differences in nutrient absorption capacity and fertilizer application in different growth stages, this model cannot be applied to guide the stage fertilization of trees. In the past, it was difficult to collect data due to the long growth period of forest. Therefore, there were few models of forest fertilization. However, in recent years, sensor network can be used to collect real-time data in forest areas with the burgeoning development of artificial intelligence. It is time to construct a predictive model of accurate fertilization of forest through the collected data. Similarly, as the traditional one, the model can only exert real-time and accurate small-scale prediction in forest areas considering the complexity of different forest areas and site indices, the diversified factors affecting forest fertilization. At the same time, since there is a high degree of nonlinearity between site index and dominant wood-related nutrient content, forest age, and other influencing factors, the traditional fertilization model often has inaccurate results due to manifold parameters or low goodness of fitting. As a result, the large amount of accumulated data can hardly serve to guide production practice. In the later stage, expert experience is still needed, resulting in a waste of manpower, material, and financial resources.Neural network is considered to be a powerful tool to solve nonlinear problems [10–15]. In order to ensure the integrity of information and improve the precision of fertilization prediction, the prediction model of multifactor precision fertilization based on neural network can be used to solve these problems. It can reflect different absorption of nutrients in different growth periods of young and middle ages of trees and apply fertilizer efficiently in the optimal fertilization period of trees. The BP algorithm shows great potential to introduce the prediction model of forest fertilization. With the updating of the data of forest fertilization experiment, the BP algorithm modifies the existing prediction model of fertilization, which makes the forecasts of the model more precise and get closer to the actual situation. However, the accuracy of predicted results and actual values is not fairly reliable, on account of some big defects, such as local optimization, slow convergence speed, and large dependence on training data, when a single BP neural network algorithm is used to predict the amount of fertilizer applied to trees. In this paper, the analyzing method of grey relativity [16–19] is introduced to extract the original variable data of multiple indexes in advance, and the factors with high correlation degree are taken as the network input layer, and then the global searching ability of PSO [20–25] is used to optimize the BP neural network, which greatly avoids the defect of the model falling into the local optimum and further improves the accuracy. Although neural network prediction has been applied in medical engineering, electric power, system, and other fields at home and abroad, and some gratifying results have been achieved, there are only few studies on the application of neural network prediction model to forest fertilizer application. Therefore, an improved GRA-PSO-BP algorithm was used to establish the accurate fertilization model of trees in stages.


## 2. Methodology

### 2.1. Test Area Overview and Plot Settings

The experimental data are from the state-run forest farm in Renhua County, Guangdong Province, and Liujiashan Forest Farm in Shixing County, Guangdong Province, which is the main production area of fir trees, with an altitude of less than 500 m, and the soil type is the red and yellow soil developed on slates. Large-scale fertilization experiments were carried out for young and mature timber of fir forests, and the experimental forests with different site indexes were classified during the experiment according to the utilization of the last round of forest land and the Nanling Mountains fir trees in the different site conditions; the experimental forest can be divided into four types: (1) triple-tilled woodlands ofCunninghamia lanceolata, with a site index of 5; (2) cultivated Chinese fir double-tilled woodlands of Cunninghamia lanceolata, with a site index of 18; (3) first ploughing woodlands after shrub felling, with a site index of 21; and (4) first ploughing woodlands after logged broad-leaf forests, with a site index of 23. Under the condition of the same site index, the slope direction of the test site is the same, the terrain difference is small, and the trees are arranged at the same height. The young forest and the middle-aged forest test area are, respectively, set in the forest stands with appropriate site index, and the area of each test site is 1 hm2. The initial planting density of young trees in the experimental forests was 3,600 plants/1 hm2. In the first 6 to 8 years, the first thinning was carried out, and the forests in the middle-aged forests were retained at 2,700 plants/hm2. The size, distribution, and density of the forests in each experimental area were basically the same. There are fixed stakes around each test area, and the dominant wood and average wood in the test area have specific numbers. The growth and nutrient dynamics of dominant wood and average wood in the test area are measured regularly every year, and nitrogen, phosphorus, and potash fertilizers are applied according to the fertilization method in the test area.16 plots with age of 5a, 19a, 22a, and 25a were selected in the above four different site types. The selected average and dominant trees in different areas were logged and then divided by leaves, branches, stems, barks, and roots to weigh each organ of them. The age of leaves and branches were weighed, respectively, and the root system was dug to a depth of 60 cm in the range of the canopy. Samples of all organs were evenly sampled and dried at 80°C to determine water and nutrient content. Meanwhile, four soil sampling points of different stand ages were set up in four different site index types, 16 soil sampling points in total. Take soil samples in 0∼60 cm soil layer of each soil sampling point. The physicochemical analysis of plant and soil samples was carried out in accordance with conventional methods, and the determination of the physicochemical properties of the soil in the experimental forest land is shown in Table1.Table 1
Test results of physical and chemical properties of the soil of experimental forest land.
Stand age (t/a)Site indexSTextureOrganic matter (g/kg−1)Total nitrogen (g/kg−1)Fast-acting nitrogen (mg/kg−1)Fast-acting phosphorus (g/kg−1)Fast-acting potassium (g/kg−1)PH5a5Middle soil9.00.3162.1574.618Middle soil17.40.8927.73.0804.921Middle soil22.01.12324.0925.023Middle soil27.01.35485.01005.219a5Middle soil8.90.3216.42.3544.518Middle soil17.20.9283.1884.821Middle soil22.11.14324.1905.123Middle soil27.31.3847.65.0965.322a5Middle soil9.00.3116.12.2524.418Middle soil17.10.8728.13.2914.721Middle soil22.01.11324.2875.423Middle soil27.11.3348.15.2945.125a5Middle soil9.10.316.22.0554.818Middle soil17.30.9128.23.0834.921Middle soil22.21.10324.3915.323Middle soil27.01.3647.95.11005.1The discrepancies of site level are due to the different locations of the samples, which results in different climatic zones, slopes, slope positions, and soil types. The differentiated use of forest land at different site levels leads to significant differences in the organic matter content, nutrient content, and fertilizer utilization rate of the humus layer in the woodlands. Consequently, it can be concluded that the availability of soil nutrients is one of the essential factors for different experimental forests.In the selection of the test plots, the age of the forest stands has been replaced by space instead of time. The accuracy of analysis is, to some extent, affected by the varied environmental conditions between the dominant wood and the average wood in the same site type. However, the effects of random errors on the basic trend of tree growth and nutrient absorption changing with forest age were relatively eliminated due to the significant difference between the growth of dominant and average stands.


### 2.2. Variation of Nutrient Content in Organs ofCunninghamia lanceolata

Table2 shows the measured values of nitrogen (N), phosphorus (P), and potassium (K) in different organs of dominant Chinese fir trees. The nutrient content of different Chinese firs varies greatly. The order of nutrient content in the various parts of the tree is as leaf, branch, bark, root, and stem. At the same time, it was also shown that the nutrient contents of different organs of Cunninghamia lanceolata with different site indexes were different at different ages. The content of nitrogen in the leaves increased with the increase of the site index and the age of the forest; the content of potassium in the leaves, nitrogen, phosphorus, potassium in the trunk, branches, and bark and the content of the phosphorus and potassium in the roots decreased with the increase of the age.Table 2
Nutrient content of dominant Chinese fir of different forest ages under different site conditions, g/kg.
OrganStand age (t/a)NitrogenPhosphorusPotassiumSite indexSSite indexSSite indexS518212351821235182123Leaf56.810.712.613.80.751.021.091.123.24.74.95.0198.612.112.412.70.841.161.281.063.34.34.74.9229.413.213.513.70.911.121.250.933.44.24.64.82510.113.914.114.20.921.041.070.973.24.14.54.7Branch53.96.77.17.30.590.710.740.773.14.04.14.2192.84.55.75.90.490.550.570.592.73.63.73.9222.74.25.55.70.470.520.550.582.63.83.93.6252.54.15.35.60.450.540.560.542.53.74.03.5Stem51.11.92.12.30.110.200.230.240.81.71.51.4191.21.41.92.00.080.090.120.140.40.91.11.3221.01.31.61.80.070.080.100.090.50.81.01.1250.91.21.41.30.040.060.080.070.40.70.90.8Bark52.84.64.54.70.210.320.340.354.73.63.33.1192.41.93.53.30.100.150.160.184.51.91.81.7222.33.63.43.10.090.110.130.154.11.71.61.5252.23.43.22.90.070.090.110.143.81.51.31.4Root51.73.53.83.90.100.270.290.313.42.12.01.8191.41.92.22.50.060.110.140.162.41.61.51.4221.61.72.12.40.050.120.120.172.11.71.61.5251.21.52.02.20.080.100.140.152.01.41.31.2Table3 shows the biomass range from average tree to dominant tree and the percentage of each organ in Chinese fir forest. Based on the data in Tables 2 and 3, the average nutrient content C of the Chinese fir forest obtained by the weighted average method is calculated as C=∑Ciwi/∑wi, where Ci and wi represent the nutrient content and biomass of leaves, branches, stems, bark, and roots of a single tree, respectively. Its product sum ∑Ciwi is the total amount of nutrients absorbed by a single tree, the total organ biomass ∑wi is the total biomass of a single tree, and the relationship between the total nutrient absorption X per unit area of Chinese fir forest and the average nutrient content C is X=nwC/1000=CW, where, n is the stand density, w is the average biomass per tree (kg/plant), and W is the biomass per unit area (t/hm2).Table 3
The biomass range (w) of a single Chinese fir and the percentage of its organs.
Stand age (t/a)OrganSite indexS = 5Site indexS = 18Site indexS = 21Site indexS = 23W (kg)Percentage (%)W (kg)Percentage (%)W (kg)Percentage (%)W (kg)Percentage (%)5aLeaf1.67 ± 0.5821.443.48 ± 0.9524.343.98 ± 1.0422.044.12 ± 1.3918.7Branch1.25 ± 0.4416.052.36 ± 0.8416.502.78 ± 1.2115.43.47 ± 1.3215.7Stem2.60 ± 0.833.404.40 ± 1.0830.775.64 ± 1.9431.236.37 ± 2.5133Bark0.53 ± 0.126.8040.96 ± 0.336.711.46 ± 0.688.081.98 ± 0.839Root1.74 ± 0.6622.333.10 ± 0.9221.684.20 ± 1.4623.255.24 ± 2.0723.6Tree7.79 ± 2.6100.0014.30 ± 4.12100.0018.06 ± 6.33100.0022.04 ± 8.12100.0019aLeaf2.64 ± 0.924.115.27 ± 1.475.045.98 ± 2.095.146.34 ± 2.874.83Branch4.43 ± 2.056.907.61 ± 2.447.288.22 ± 3.087.069.12 ± 4.437.00Stem40.36 ± 15.763.0064.13 ± 20.661.3472.33 ± 24.362.1281.2 ± 29.762Bark6.78 ± 2.1110.6011.30 ± 4.8710.8413.11 ± 5.1111.315.10 ± 7.0512Root9.97 ± 3.011616.20 ± 5.5615.5016.80 ± 6.7614.3819.40 ± 8.1114.8Tree64.18 ± 23.79100.00104.51 ± 34.94100.00116.44 ± 41.34100.00131.16 ± 44.86100.0022aLeaf4.22 ± 1.274.666.37 ± 2.024.896.89 ± 2.944.767.26 ± 3.164.53Branch5.11 ± 2.315.648.78 ± 3.116.749.89 ± 3.846.8410.43 ± 4.386.50Stem56.7 ± 20.462.677.32 ± 26.2159.3687.34 ± 32.3060.4598.53 ± 41.161.4Bark9.12 ± 4.110.114.40 ± 5.9811.115.21 ± 6.1210.5316.32 ± 7.1410.2Root15.5 ± 6.5217.123.40 ± 9.7017.9125.16 ± 11.617.4127.89 ± 13.117.4Tree90.65 ± 34.6100.00130.27 ± 47.02100.00144.49 ± 56.8100.00160.43 ± 68.8100.0025aLeaf5.64 ± 1.625.38.65 ± 2.615.709.16 ± 3.165.6010.38 ± 4.215.73Branch6.10 ± 2.655.79.76 ± 3.466.4310.78 ± 3.916.6011.70 ± 4.616.50Stem63.2 ± 24.76088.70 ± 28.7058.596.80 ± 31.359.10108.3 ± 35.959.8Bark11.7 ± 5.221115.90 ± 6.2110.516.80 ± 7.4110.2617.70 ± 8.019.8Root20.4 ± 7.5919.128.70 ± 9.5318.8730.24 ± 11.218.533.12 ± 13.4118.3Tree107.04 ± 41.78100.00151.71 ± 50.51100.00163.78 ± 56.98100.00181.2 ± 66.14100.00


### 2.3. Analysis of Influencing Factors

There are a series of important indexes to calculate the amount of fertilizer applied to trees, such as site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood biomass, fertilizer utilization rate, and target yield. The site index is a natural environmental factor required for forest production and is a collection of conditions such as moisture, temperature, light intensity, and soil fertility. In the forest area with high site index or low site index, the fertilization effect is difficult to show. In the forest area with high site index, the soil fertility, water, temperature, light, and other conditions are very good. After fertilization, the effect on the increase of forest yield is small. On the contrary, the soil fertility, water, temperature, light, and other conditions are very poor in the forest land with low site index. After fertilization, due to natural factors, the fertilization effect is also very poor. In the effective site index interval, the theoretical limit value of forest production is the dominant wood biomass of forest land in practice, so the dominant wood nutrient content is the limit value that can be reached by the average wood nutrient content after fertilization. If the difference between the dominant wood biomass and the average wood biomass is greater, the fertilization effect will be better. In addition, the target yield is also related to the amount of fertilization. The target yield is determined by the input-output ratio of forest land in previous years, and the target yield determines the target increase yield, that is, the amount of fertilization. The relative intensities of the absorbed nutrients of different forest ages are different. The best fertilization period, the more fertilizer is absorbed. At the same time, the utilization rate of the fertilizer is related to the physical and chemical properties of the soil. The physical and chemical properties of the soil in different forest land are different, so the utilization rate of the fertilizer is also different. The higher the utilization rate, the better the fertilizer effect. Generally, three kinds of fixed nitrogen, phosphorus, and potassium fertilizers are applied in forest land. The utilization rate of fertilizer is observed according to the different physical and chemical properties of soil in different forest land. Therefore, seven factors including site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood biomass, fertilizer utilization rate, and target yield were selected as the influencing factors of forest fertilizer application.


### 2.4. Data Normalization

Three hundred tree growth and nutrient measurement records were obtained. The data of each group wereX1 site index, X2 forest age, X3 dominant wood-related nutrient concentration, X4 dominant wood biomass, X5 average wood biomass, X6 target production, and X7 fixed fertilizer utilization rate. Since the input parameters have different dimensions and the orders of magnitude of units are also greatly different, these will have a great impact on the training of neural network. Therefore, the input parameters need to be normalized. In order to improve the training efficiency, an accurate fertilization model is established. The minimum maximum transformation method is used to normalize different data and deal with the dimensional influence between data. The parameters are in the range of [−1, 1], and then the parameters of neural network are trained. The following formula is adopted for normalization:(1)y=x−xminxmax−xmin.Map the data to (−1, +1) and replace the formula with(2)y=y∗2−1,where y is the normalized value, x is the original value of a parameter, xmax is the maximum value in the set of values, and xmin is the minimum value in the set of values.


### 2.5. Basic BPNN Model

The structure of the basic BPNN consists of the input layer, the hidden layer, and the output layer, as shown in Figure1. The BPNN algorithm consists of two parts: the forward transfer of the input data and the back propagation of the error between the output data and the expected data [26, 27]. One part is the input data through the input layer to the hidden layer, the hidden layer to the output layer. In the process, the calculation will be based on the given initial weight and threshold, and finally the output data will be obtained. The other part is to calculate the error change value δj (k) between the actual output and the expected value and then turn to the back propagation, and the error signal δj (k) is propagated back through the original connection path through the network to modify the weight of each layer of neurons. The value ωij (jk) and the threshold θj (k) until the target accuracy are reached.Figure 1
BP neural network structure.The hidden layer of BPNN can be one or more layers. It is proved theoretically that the neural network of a single hidden layer can approximate the nonlinear function with arbitrary precision, so that the model can realize the nonlinear mapping from input to output. With the increase of the number of hidden layers, the output error of the network will decrease [28]. Therefore, the increase of hidden layers will improve the accuracy of the network but will make the network structure become complex, the running time become longer, and even lead to the overfitting phenomenon. Therefore, after the number of neurons in the hidden layer is determined by empirical formula, the structure of BPNN tree fertilization model is 3-layer structure of 7-12-1.The BPNN has strong nonlinear mapping ability, self-adaptive and good self-learning ability, and strong fault tolerance. In the absence of a mathematical model, the nonlinear mapping between the amount of fertilizer applied to trees and the factors affecting the effect of fertilization can be realized. By learning and training the relationship between input and output data, when the actual output value and the expected value have a large error, the weights and thresholds of each layer of the neural network are constantly updated until the accuracy of network error is reached, so as to better improve the model.In the figure,i, j, and k are the number of neurons in each layer; xi is any input signal of each group of data in the input p in the input layer; ωij, ωjk and θj, θk, respectively, represent the weight and threshold of each layer; φ (·) and τ (·) represent the hidden layer and output layer activation functions, respectively. In the BP neural network algorithm, the hidden layer neurons input the signal netj when the sample p acts.(3)netj=∑ωijxi+θj.The hidden layer neuron node outputs a signaloj when the sample p acts:(4)oj=φnetj=φ∑ωijxi+θj.The output layer neuron node inputs the signal netk when the sample p acts:(5)netk=∑ωjkOj+θk=∑ωjkφ∑ωijxi+θj+θk.The output layer neurons output the signalyk when the sample p acts:(6)yk=τnetk=τ∑ωjkφ∑ωjxi+θj+θk.In the formula,η is the error back propagation learning rate (0 < η < 1), δj (k) is the error signal of each layer, and the weight and threshold are updated and iterated according to(7)Δωijjk=ηδjkOij,(8)Δθjk=ηδjk.The error variation value of each layer gradually corrects the fertilization model during the back propagation process.


### 2.6. Particle Swarm Optimization BPNN Model

When a single BPNN model is used to predict the precise fertilization amount of forest trees, the error between the actual output and the expected value is still large, because in the actual use, the convergence speed of the network is relatively slow, and it is easy to fall into the local optimization, even when there are few training samples, and there may be a fitting problem. Particle swarm optimization algorithm has better global optimization ability. Therefore, the particle swarm optimization algorithm was proposed to optimize the BPNN model. The particle swarm optimization (PSO) algorithm is a swarm intelligence optimization algorithm [29, 30]. The velocity characteristics, fitness value, and position were used to represent the motion characteristics of the particle. After initializing the potential optimal solution particles, the speed of the particles determines the search direction and distance of the particles, the fitness value determines the quality of the particles. When the particles move in the preset space, the position is changed according to the individual optimal solution and the global optimal solution. By updating the fitness value of the particle, the fitness value is the smallest. The position corresponding to the particle is the optimal solution, as shown in Figure 2, the motion of particles in the particle swarm leaving its current position in the search area, the velocity vi+1 that this particle is about to update is a combination of various factors, It includes vi, vp, and vG, which are the current particle speed, the past best experience of individual particles, and the best experience of population particles.Figure 2
Particle movement in a swarm.The weights and thresholds of each layer in the BPNN can be represented by individual particle swarms. The individual swarms are coded, represented by the position vector of the particle swarm, and the optimal population particles are output by the iterative algorithm. After decoding, the BPNN global optimal is obtained. Weight and threshold establish a PSO-BPNN algorithm model. After decoding, the global optimal weight and threshold of BPNN are obtained, and the pso-bpnn algorithm model with 3-layer structure of 7-12-1 is established.


### 2.7. GRA-PSO-BPNN Model Construction Process

The degree of correlation can delicate the correlation in mutual matters and factors, and the Grey correlation degree analysis method can express the changeable situation in mutual matters and factors. Taking the timeliness and regionality of the fertilization model into consideration, different forests have different factors that affect the accuracy of fertilization, resulting that there is no certain model to implement the fertilization model to achieve the accurate prediction so that the accurate rate of the prediction is not so well. This research adopts the analysis method of Grey Relational degree which can analyze the key factors that influence the fertilizer application amount in different areas and on the basis of the annual forest data updated by artificial forest, which can make the accurate fertilization come true. The grey relation analysis method is used to calculate the grey relation between input variables and output results. The specific steps are as follows:Step 1: Using the formula (9), establish a data matrix xi(9)xi=xi1,xi2,xi3,…,xi300,whereXi (1), Xi (2), Xi (3),…,Xi (300) represents the number of sequences. Among them, i represents site index, forest age, nutrient concentration related to dominant wood, biomass of dominant wood, average biomass of wood, fixed fertilizer utilization rate, and target yield, which are different factors affecting the fertilization effect.Step 2: Establish an initialization change matrixxi′ using(10)xi′=xi1xi1,xi2xi1,xi3xi1,…,xi300xi1=xi′1,xi′2,xi′3,…,xi′300.Step 3: Use formula (11) to calculate the difference sequence, where k represents the sequence number:(11)Δoik=absxik−xi′=Δoi1,Δoi2,…,Δoi300.Step 4: Calculate the correlation coefficient and the grey correlation degree by using equations (12) and (13).Step 4: Use formulas (12) and (13) to calculate correlation coefficient and grey correlation degree. φ is the resolution coefficient, 0 < φ < 1. In this study, φ = 0.5:(12)ξoik=miniminkΔoik+φmaximaxkΔoikΔoik+φminiminkΔoik,(13)γoi=1300∑1300ξoik.According to the calculation formula of grey correlation coefficient, the minimum differenceminiminkΔoik and the maximum difference maximaxkΔoik of N, P, and K fertilizers were 8.822e − 06, 3.483; 7.1719e − 05, 7.307; 2.388e − 04, 1.201, respectively; Table 4 shows the grey correlation coefficient and grey correlation degree between the input variables and the output of the prediction model of forest precise fertilization of nitrogen, phosphorus, and potassium.Table 4
Grey correlation degree and grey correlation coefficient of each influencing factor and output fertilization amount.
XNiεoiroiXPiεoiroiXKiεoiroiX1εo1 = (0.9178,…,0.8813)0.9105X1εo1 = (0.9512,…,0.8731)0.9334X1εo1 = (0.7483,…,0.4639)0.7512X2εo2 = (0.9259,…,0.8947)0.9122X2εo2 =(0.9554,…,0.8793)0.9357X2εo2 =(0.7644,…,0.4744)0.7559X3εo3 = (0.8972,…,0.9323)0.9115X3εo3 =(0.9372,…,0.9191)0.9441X3εo3 = (0.7312,…,0.6561)0.7606X4εo4 =(0.9240,…,0.9390)0.9232X4εo4 =(0.9544,…,0.8992)0.9398X4εo4 =(0.7605,…,0.5119)0.7715X5εo5 =(0.9338,…,0.9554)0.9009X5εo5 =(0.9594,…,0.9063)0.9321X5εo5 =(0.7800,…,0.5262)0.7372X6εo6 =(0.9194,…,0.9516)0.9119X6εo6 =(0.9520,…,0.9046)0.9368X6εo6 =(0.7514,…,0.5229)0.7561X7εo7 =(0.6660,…,0.6646)0.6142X7εo7 =(0.6016,…,0.6144)0.6736X7εo7 =(0.3952,…,0.3762)0.5868According to the above analysis, the grey correlation degree ofX1 site index, X2 forest age, X3 dominant wood-related nutrient concentration, X4 dominant wood biomass, X5 average wood biomass, and X6 target yield is higher, and X7 fixed fertilizer utilization rate is smaller; this study uses data reduction, grey correlation analysis, to reduce the influencing factors of different forest areas for predicting and calculating forest fertilization, identifying key factors affecting forest fertilization and providing more effective input for the prediction model.Step 1 defines the input and output of the GRA-PSO-BPNN model.Compared with the PSO-BPNN model, using the grey correlation analysis, the site index, the forest age, the dominant wood-related nutrient content, the dominant wood biomass, the average wood biomass, and the target yield are used as the network input, and the actual fertilization amount is used as the output. Therefore, there are 300 records, 200 for training and 100 for validation. In the formula, the number of neurons in the input layerx is 6, and the number of neurons in the output layer y is 1. The empirical formula is used to determine the number m of hidden neurons:(14)m=x+y+a,a∈1,10.The reduction of network error can not only adjust the number of hidden layer neurons, but also increase the number of hidden layer. However, increasing the number of hidden layers will complicate the network, increase the training time, and even overfitting. It is biased towards selecting the single hidden layer BPNN [31]. As mentioned above, the structure of GRA-PSO-BPNN model is determined to be a 3-layer structure of 6-12-1.Step 2: Particle Swarm Initialization.The weights and thresholds in the BPNN model are the parameters to be optimized by the PSO algorithm. At the beginning, a certain number of particle swarm individuals can be randomly generated to represent, and the global optimal position gbest and the individual optimal position pbest of the particle are initialized. The sum of ownership value and threshold value in network structure is the dimensiond of particle swarm individual search space, that is,(15)d=xm+my+m+y.Thed-dimensional vector xit=xi1t,xi2t,…,xidt represents the position vector of the ith particle in the tth generation population.Step 3: Particle Velocity Location Update.According to the following formula, the particle velocity is updated according to the individual optimal solution and the global optimal solution.(16)vt+1=w⋅vt+c1⋅r1⋅pbestt−xt+c2⋅r2⋅gbestt−xt.Among them, pbest (t) and gbest (t) are the individual optimal solution of the ith particle in the tth generation of the d-dimensional space, and the global optimal solution in the tth generation, r1 and r2 are in [0, 1], c1, c2 are learning factors, w is the inertia weight, and the particle position is updated according to the following formula:(17)xt+1=xt+vt+1.The particle group can be evaluated by calculating the fitness value. The smaller the fitness value of the particle group is, the higher the fitness is. The speed and position of each particle are adjusted based on the fitness value. The particle fitness is calculated according to the following formula. ValueF is(18)F=∑i=1Nabsyi−ti,whereN is the number of samples in the test area; abs is the absolute value function; yi is the actual value of the sample i; and ti is the predicted value of the sample i.Step 4: Optimal Population Particles:When the number of population evolution reaches the upper limitT or the iteration error reaches the set precision e, when the algorithm stops, the global optimal solution can be obtained and mapped to the weight and threshold of the BP network.Step 5: GRA-PSO-BPNN Model Training.After determining the GRA-PSO-BPNN weights and thresholds, the training data are input to train the GRA-PSO-BPNN model. The entire model flow is shown in Figure3.Figure 3
GRA-PSO-BPNN model flow chart.


## 2.1. Test Area Overview and Plot Settings

The experimental data are from the state-run forest farm in Renhua County, Guangdong Province, and Liujiashan Forest Farm in Shixing County, Guangdong Province, which is the main production area of fir trees, with an altitude of less than 500 m, and the soil type is the red and yellow soil developed on slates. Large-scale fertilization experiments were carried out for young and mature timber of fir forests, and the experimental forests with different site indexes were classified during the experiment according to the utilization of the last round of forest land and the Nanling Mountains fir trees in the different site conditions; the experimental forest can be divided into four types: (1) triple-tilled woodlands ofCunninghamia lanceolata, with a site index of 5; (2) cultivated Chinese fir double-tilled woodlands of Cunninghamia lanceolata, with a site index of 18; (3) first ploughing woodlands after shrub felling, with a site index of 21; and (4) first ploughing woodlands after logged broad-leaf forests, with a site index of 23. Under the condition of the same site index, the slope direction of the test site is the same, the terrain difference is small, and the trees are arranged at the same height. The young forest and the middle-aged forest test area are, respectively, set in the forest stands with appropriate site index, and the area of each test site is 1 hm2. The initial planting density of young trees in the experimental forests was 3,600 plants/1 hm2. In the first 6 to 8 years, the first thinning was carried out, and the forests in the middle-aged forests were retained at 2,700 plants/hm2. The size, distribution, and density of the forests in each experimental area were basically the same. There are fixed stakes around each test area, and the dominant wood and average wood in the test area have specific numbers. The growth and nutrient dynamics of dominant wood and average wood in the test area are measured regularly every year, and nitrogen, phosphorus, and potash fertilizers are applied according to the fertilization method in the test area.16 plots with age of 5a, 19a, 22a, and 25a were selected in the above four different site types. The selected average and dominant trees in different areas were logged and then divided by leaves, branches, stems, barks, and roots to weigh each organ of them. The age of leaves and branches were weighed, respectively, and the root system was dug to a depth of 60 cm in the range of the canopy. Samples of all organs were evenly sampled and dried at 80°C to determine water and nutrient content. Meanwhile, four soil sampling points of different stand ages were set up in four different site index types, 16 soil sampling points in total. Take soil samples in 0∼60 cm soil layer of each soil sampling point. The physicochemical analysis of plant and soil samples was carried out in accordance with conventional methods, and the determination of the physicochemical properties of the soil in the experimental forest land is shown in Table1.Table 1
Test results of physical and chemical properties of the soil of experimental forest land.
Stand age (t/a)Site indexSTextureOrganic matter (g/kg−1)Total nitrogen (g/kg−1)Fast-acting nitrogen (mg/kg−1)Fast-acting phosphorus (g/kg−1)Fast-acting potassium (g/kg−1)PH5a5Middle soil9.00.3162.1574.618Middle soil17.40.8927.73.0804.921Middle soil22.01.12324.0925.023Middle soil27.01.35485.01005.219a5Middle soil8.90.3216.42.3544.518Middle soil17.20.9283.1884.821Middle soil22.11.14324.1905.123Middle soil27.31.3847.65.0965.322a5Middle soil9.00.3116.12.2524.418Middle soil17.10.8728.13.2914.721Middle soil22.01.11324.2875.423Middle soil27.11.3348.15.2945.125a5Middle soil9.10.316.22.0554.818Middle soil17.30.9128.23.0834.921Middle soil22.21.10324.3915.323Middle soil27.01.3647.95.11005.1The discrepancies of site level are due to the different locations of the samples, which results in different climatic zones, slopes, slope positions, and soil types. The differentiated use of forest land at different site levels leads to significant differences in the organic matter content, nutrient content, and fertilizer utilization rate of the humus layer in the woodlands. Consequently, it can be concluded that the availability of soil nutrients is one of the essential factors for different experimental forests.In the selection of the test plots, the age of the forest stands has been replaced by space instead of time. The accuracy of analysis is, to some extent, affected by the varied environmental conditions between the dominant wood and the average wood in the same site type. However, the effects of random errors on the basic trend of tree growth and nutrient absorption changing with forest age were relatively eliminated due to the significant difference between the growth of dominant and average stands.


## 2.2. Variation of Nutrient Content in Organs ofCunninghamia lanceolata

Table2 shows the measured values of nitrogen (N), phosphorus (P), and potassium (K) in different organs of dominant Chinese fir trees. The nutrient content of different Chinese firs varies greatly. The order of nutrient content in the various parts of the tree is as leaf, branch, bark, root, and stem. At the same time, it was also shown that the nutrient contents of different organs of Cunninghamia lanceolata with different site indexes were different at different ages. The content of nitrogen in the leaves increased with the increase of the site index and the age of the forest; the content of potassium in the leaves, nitrogen, phosphorus, potassium in the trunk, branches, and bark and the content of the phosphorus and potassium in the roots decreased with the increase of the age.Table 2
Nutrient content of dominant Chinese fir of different forest ages under different site conditions, g/kg.
OrganStand age (t/a)NitrogenPhosphorusPotassiumSite indexSSite indexSSite indexS518212351821235182123Leaf56.810.712.613.80.751.021.091.123.24.74.95.0198.612.112.412.70.841.161.281.063.34.34.74.9229.413.213.513.70.911.121.250.933.44.24.64.82510.113.914.114.20.921.041.070.973.24.14.54.7Branch53.96.77.17.30.590.710.740.773.14.04.14.2192.84.55.75.90.490.550.570.592.73.63.73.9222.74.25.55.70.470.520.550.582.63.83.93.6252.54.15.35.60.450.540.560.542.53.74.03.5Stem51.11.92.12.30.110.200.230.240.81.71.51.4191.21.41.92.00.080.090.120.140.40.91.11.3221.01.31.61.80.070.080.100.090.50.81.01.1250.91.21.41.30.040.060.080.070.40.70.90.8Bark52.84.64.54.70.210.320.340.354.73.63.33.1192.41.93.53.30.100.150.160.184.51.91.81.7222.33.63.43.10.090.110.130.154.11.71.61.5252.23.43.22.90.070.090.110.143.81.51.31.4Root51.73.53.83.90.100.270.290.313.42.12.01.8191.41.92.22.50.060.110.140.162.41.61.51.4221.61.72.12.40.050.120.120.172.11.71.61.5251.21.52.02.20.080.100.140.152.01.41.31.2Table3 shows the biomass range from average tree to dominant tree and the percentage of each organ in Chinese fir forest. Based on the data in Tables 2 and 3, the average nutrient content C of the Chinese fir forest obtained by the weighted average method is calculated as C=∑Ciwi/∑wi, where Ci and wi represent the nutrient content and biomass of leaves, branches, stems, bark, and roots of a single tree, respectively. Its product sum ∑Ciwi is the total amount of nutrients absorbed by a single tree, the total organ biomass ∑wi is the total biomass of a single tree, and the relationship between the total nutrient absorption X per unit area of Chinese fir forest and the average nutrient content C is X=nwC/1000=CW, where, n is the stand density, w is the average biomass per tree (kg/plant), and W is the biomass per unit area (t/hm2).Table 3
The biomass range (w) of a single Chinese fir and the percentage of its organs.
Stand age (t/a)OrganSite indexS = 5Site indexS = 18Site indexS = 21Site indexS = 23W (kg)Percentage (%)W (kg)Percentage (%)W (kg)Percentage (%)W (kg)Percentage (%)5aLeaf1.67 ± 0.5821.443.48 ± 0.9524.343.98 ± 1.0422.044.12 ± 1.3918.7Branch1.25 ± 0.4416.052.36 ± 0.8416.502.78 ± 1.2115.43.47 ± 1.3215.7Stem2.60 ± 0.833.404.40 ± 1.0830.775.64 ± 1.9431.236.37 ± 2.5133Bark0.53 ± 0.126.8040.96 ± 0.336.711.46 ± 0.688.081.98 ± 0.839Root1.74 ± 0.6622.333.10 ± 0.9221.684.20 ± 1.4623.255.24 ± 2.0723.6Tree7.79 ± 2.6100.0014.30 ± 4.12100.0018.06 ± 6.33100.0022.04 ± 8.12100.0019aLeaf2.64 ± 0.924.115.27 ± 1.475.045.98 ± 2.095.146.34 ± 2.874.83Branch4.43 ± 2.056.907.61 ± 2.447.288.22 ± 3.087.069.12 ± 4.437.00Stem40.36 ± 15.763.0064.13 ± 20.661.3472.33 ± 24.362.1281.2 ± 29.762Bark6.78 ± 2.1110.6011.30 ± 4.8710.8413.11 ± 5.1111.315.10 ± 7.0512Root9.97 ± 3.011616.20 ± 5.5615.5016.80 ± 6.7614.3819.40 ± 8.1114.8Tree64.18 ± 23.79100.00104.51 ± 34.94100.00116.44 ± 41.34100.00131.16 ± 44.86100.0022aLeaf4.22 ± 1.274.666.37 ± 2.024.896.89 ± 2.944.767.26 ± 3.164.53Branch5.11 ± 2.315.648.78 ± 3.116.749.89 ± 3.846.8410.43 ± 4.386.50Stem56.7 ± 20.462.677.32 ± 26.2159.3687.34 ± 32.3060.4598.53 ± 41.161.4Bark9.12 ± 4.110.114.40 ± 5.9811.115.21 ± 6.1210.5316.32 ± 7.1410.2Root15.5 ± 6.5217.123.40 ± 9.7017.9125.16 ± 11.617.4127.89 ± 13.117.4Tree90.65 ± 34.6100.00130.27 ± 47.02100.00144.49 ± 56.8100.00160.43 ± 68.8100.0025aLeaf5.64 ± 1.625.38.65 ± 2.615.709.16 ± 3.165.6010.38 ± 4.215.73Branch6.10 ± 2.655.79.76 ± 3.466.4310.78 ± 3.916.6011.70 ± 4.616.50Stem63.2 ± 24.76088.70 ± 28.7058.596.80 ± 31.359.10108.3 ± 35.959.8Bark11.7 ± 5.221115.90 ± 6.2110.516.80 ± 7.4110.2617.70 ± 8.019.8Root20.4 ± 7.5919.128.70 ± 9.5318.8730.24 ± 11.218.533.12 ± 13.4118.3Tree107.04 ± 41.78100.00151.71 ± 50.51100.00163.78 ± 56.98100.00181.2 ± 66.14100.00


## 2.3. Analysis of Influencing Factors

There are a series of important indexes to calculate the amount of fertilizer applied to trees, such as site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood biomass, fertilizer utilization rate, and target yield. The site index is a natural environmental factor required for forest production and is a collection of conditions such as moisture, temperature, light intensity, and soil fertility. In the forest area with high site index or low site index, the fertilization effect is difficult to show. In the forest area with high site index, the soil fertility, water, temperature, light, and other conditions are very good. After fertilization, the effect on the increase of forest yield is small. On the contrary, the soil fertility, water, temperature, light, and other conditions are very poor in the forest land with low site index. After fertilization, due to natural factors, the fertilization effect is also very poor. In the effective site index interval, the theoretical limit value of forest production is the dominant wood biomass of forest land in practice, so the dominant wood nutrient content is the limit value that can be reached by the average wood nutrient content after fertilization. If the difference between the dominant wood biomass and the average wood biomass is greater, the fertilization effect will be better. In addition, the target yield is also related to the amount of fertilization. The target yield is determined by the input-output ratio of forest land in previous years, and the target yield determines the target increase yield, that is, the amount of fertilization. The relative intensities of the absorbed nutrients of different forest ages are different. The best fertilization period, the more fertilizer is absorbed. At the same time, the utilization rate of the fertilizer is related to the physical and chemical properties of the soil. The physical and chemical properties of the soil in different forest land are different, so the utilization rate of the fertilizer is also different. The higher the utilization rate, the better the fertilizer effect. Generally, three kinds of fixed nitrogen, phosphorus, and potassium fertilizers are applied in forest land. The utilization rate of fertilizer is observed according to the different physical and chemical properties of soil in different forest land. Therefore, seven factors including site index, forest age, nutrient concentration of dominant trees, dominant wood biomass, average wood biomass, fertilizer utilization rate, and target yield were selected as the influencing factors of forest fertilizer application.


## 2.4. Data Normalization

Three hundred tree growth and nutrient measurement records were obtained. The data of each group wereX1 site index, X2 forest age, X3 dominant wood-related nutrient concentration, X4 dominant wood biomass, X5 average wood biomass, X6 target production, and X7 fixed fertilizer utilization rate. Since the input parameters have different dimensions and the orders of magnitude of units are also greatly different, these will have a great impact on the training of neural network. Therefore, the input parameters need to be normalized. In order to improve the training efficiency, an accurate fertilization model is established. The minimum maximum transformation method is used to normalize different data and deal with the dimensional influence between data. The parameters are in the range of [−1, 1], and then the parameters of neural network are trained. The following formula is adopted for normalization:(1)y=x−xminxmax−xmin.Map the data to (−1, +1) and replace the formula with(2)y=y∗2−1,where y is the normalized value, x is the original value of a parameter, xmax is the maximum value in the set of values, and xmin is the minimum value in the set of values.


## 2.5. Basic BPNN Model

The structure of the basic BPNN consists of the input layer, the hidden layer, and the output layer, as shown in Figure1. The BPNN algorithm consists of two parts: the forward transfer of the input data and the back propagation of the error between the output data and the expected data [26, 27]. One part is the input data through the input layer to the hidden layer, the hidden layer to the output layer. In the process, the calculation will be based on the given initial weight and threshold, and finally the output data will be obtained. The other part is to calculate the error change value δj (k) between the actual output and the expected value and then turn to the back propagation, and the error signal δj (k) is propagated back through the original connection path through the network to modify the weight of each layer of neurons. The value ωij (jk) and the threshold θj (k) until the target accuracy are reached.Figure 1
BP neural network structure.The hidden layer of BPNN can be one or more layers. It is proved theoretically that the neural network of a single hidden layer can approximate the nonlinear function with arbitrary precision, so that the model can realize the nonlinear mapping from input to output. With the increase of the number of hidden layers, the output error of the network will decrease [28]. Therefore, the increase of hidden layers will improve the accuracy of the network but will make the network structure become complex, the running time become longer, and even lead to the overfitting phenomenon. Therefore, after the number of neurons in the hidden layer is determined by empirical formula, the structure of BPNN tree fertilization model is 3-layer structure of 7-12-1.The BPNN has strong nonlinear mapping ability, self-adaptive and good self-learning ability, and strong fault tolerance. In the absence of a mathematical model, the nonlinear mapping between the amount of fertilizer applied to trees and the factors affecting the effect of fertilization can be realized. By learning and training the relationship between input and output data, when the actual output value and the expected value have a large error, the weights and thresholds of each layer of the neural network are constantly updated until the accuracy of network error is reached, so as to better improve the model.In the figure,i, j, and k are the number of neurons in each layer; xi is any input signal of each group of data in the input p in the input layer; ωij, ωjk and θj, θk, respectively, represent the weight and threshold of each layer; φ (·) and τ (·) represent the hidden layer and output layer activation functions, respectively. In the BP neural network algorithm, the hidden layer neurons input the signal netj when the sample p acts.(3)netj=∑ωijxi+θj.The hidden layer neuron node outputs a signaloj when the sample p acts:(4)oj=φnetj=φ∑ωijxi+θj.The output layer neuron node inputs the signal netk when the sample p acts:(5)netk=∑ωjkOj+θk=∑ωjkφ∑ωijxi+θj+θk.The output layer neurons output the signalyk when the sample p acts:(6)yk=τnetk=τ∑ωjkφ∑ωjxi+θj+θk.In the formula,η is the error back propagation learning rate (0 < η < 1), δj (k) is the error signal of each layer, and the weight and threshold are updated and iterated according to(7)Δωijjk=ηδjkOij,(8)Δθjk=ηδjk.The error variation value of each layer gradually corrects the fertilization model during the back propagation process.


## 2.6. Particle Swarm Optimization BPNN Model

When a single BPNN model is used to predict the precise fertilization amount of forest trees, the error between the actual output and the expected value is still large, because in the actual use, the convergence speed of the network is relatively slow, and it is easy to fall into the local optimization, even when there are few training samples, and there may be a fitting problem. Particle swarm optimization algorithm has better global optimization ability. Therefore, the particle swarm optimization algorithm was proposed to optimize the BPNN model. The particle swarm optimization (PSO) algorithm is a swarm intelligence optimization algorithm [29, 30]. The velocity characteristics, fitness value, and position were used to represent the motion characteristics of the particle. After initializing the potential optimal solution particles, the speed of the particles determines the search direction and distance of the particles, the fitness value determines the quality of the particles. When the particles move in the preset space, the position is changed according to the individual optimal solution and the global optimal solution. By updating the fitness value of the particle, the fitness value is the smallest. The position corresponding to the particle is the optimal solution, as shown in Figure 2, the motion of particles in the particle swarm leaving its current position in the search area, the velocity vi+1 that this particle is about to update is a combination of various factors, It includes vi, vp, and vG, which are the current particle speed, the past best experience of individual particles, and the best experience of population particles.Figure 2
Particle movement in a swarm.The weights and thresholds of each layer in the BPNN can be represented by individual particle swarms. The individual swarms are coded, represented by the position vector of the particle swarm, and the optimal population particles are output by the iterative algorithm. After decoding, the BPNN global optimal is obtained. Weight and threshold establish a PSO-BPNN algorithm model. After decoding, the global optimal weight and threshold of BPNN are obtained, and the pso-bpnn algorithm model with 3-layer structure of 7-12-1 is established.


## 2.7. GRA-PSO-BPNN Model Construction Process

The degree of correlation can delicate the correlation in mutual matters and factors, and the Grey correlation degree analysis method can express the changeable situation in mutual matters and factors. Taking the timeliness and regionality of the fertilization model into consideration, different forests have different factors that affect the accuracy of fertilization, resulting that there is no certain model to implement the fertilization model to achieve the accurate prediction so that the accurate rate of the prediction is not so well. This research adopts the analysis method of Grey Relational degree which can analyze the key factors that influence the fertilizer application amount in different areas and on the basis of the annual forest data updated by artificial forest, which can make the accurate fertilization come true. The grey relation analysis method is used to calculate the grey relation between input variables and output results. The specific steps are as follows:Step 1: Using the formula (9), establish a data matrix xi(9)xi=xi1,xi2,xi3,…,xi300,whereXi (1), Xi (2), Xi (3),…,Xi (300) represents the number of sequences. Among them, i represents site index, forest age, nutrient concentration related to dominant wood, biomass of dominant wood, average biomass of wood, fixed fertilizer utilization rate, and target yield, which are different factors affecting the fertilization effect.Step 2: Establish an initialization change matrixxi′ using(10)xi′=xi1xi1,xi2xi1,xi3xi1,…,xi300xi1=xi′1,xi′2,xi′3,…,xi′300.Step 3: Use formula (11) to calculate the difference sequence, where k represents the sequence number:(11)Δoik=absxik−xi′=Δoi1,Δoi2,…,Δoi300.Step 4: Calculate the correlation coefficient and the grey correlation degree by using equations (12) and (13).Step 4: Use formulas (12) and (13) to calculate correlation coefficient and grey correlation degree. φ is the resolution coefficient, 0 < φ < 1. In this study, φ = 0.5:(12)ξoik=miniminkΔoik+φmaximaxkΔoikΔoik+φminiminkΔoik,(13)γoi=1300∑1300ξoik.According to the calculation formula of grey correlation coefficient, the minimum differenceminiminkΔoik and the maximum difference maximaxkΔoik of N, P, and K fertilizers were 8.822e − 06, 3.483; 7.1719e − 05, 7.307; 2.388e − 04, 1.201, respectively; Table 4 shows the grey correlation coefficient and grey correlation degree between the input variables and the output of the prediction model of forest precise fertilization of nitrogen, phosphorus, and potassium.Table 4
Grey correlation degree and grey correlation coefficient of each influencing factor and output fertilization amount.
XNiεoiroiXPiεoiroiXKiεoiroiX1εo1 = (0.9178,…,0.8813)0.9105X1εo1 = (0.9512,…,0.8731)0.9334X1εo1 = (0.7483,…,0.4639)0.7512X2εo2 = (0.9259,…,0.8947)0.9122X2εo2 =(0.9554,…,0.8793)0.9357X2εo2 =(0.7644,…,0.4744)0.7559X3εo3 = (0.8972,…,0.9323)0.9115X3εo3 =(0.9372,…,0.9191)0.9441X3εo3 = (0.7312,…,0.6561)0.7606X4εo4 =(0.9240,…,0.9390)0.9232X4εo4 =(0.9544,…,0.8992)0.9398X4εo4 =(0.7605,…,0.5119)0.7715X5εo5 =(0.9338,…,0.9554)0.9009X5εo5 =(0.9594,…,0.9063)0.9321X5εo5 =(0.7800,…,0.5262)0.7372X6εo6 =(0.9194,…,0.9516)0.9119X6εo6 =(0.9520,…,0.9046)0.9368X6εo6 =(0.7514,…,0.5229)0.7561X7εo7 =(0.6660,…,0.6646)0.6142X7εo7 =(0.6016,…,0.6144)0.6736X7εo7 =(0.3952,…,0.3762)0.5868According to the above analysis, the grey correlation degree ofX1 site index, X2 forest age, X3 dominant wood-related nutrient concentration, X4 dominant wood biomass, X5 average wood biomass, and X6 target yield is higher, and X7 fixed fertilizer utilization rate is smaller; this study uses data reduction, grey correlation analysis, to reduce the influencing factors of different forest areas for predicting and calculating forest fertilization, identifying key factors affecting forest fertilization and providing more effective input for the prediction model.Step 1 defines the input and output of the GRA-PSO-BPNN model.Compared with the PSO-BPNN model, using the grey correlation analysis, the site index, the forest age, the dominant wood-related nutrient content, the dominant wood biomass, the average wood biomass, and the target yield are used as the network input, and the actual fertilization amount is used as the output. Therefore, there are 300 records, 200 for training and 100 for validation. In the formula, the number of neurons in the input layerx is 6, and the number of neurons in the output layer y is 1. The empirical formula is used to determine the number m of hidden neurons:(14)m=x+y+a,a∈1,10.The reduction of network error can not only adjust the number of hidden layer neurons, but also increase the number of hidden layer. However, increasing the number of hidden layers will complicate the network, increase the training time, and even overfitting. It is biased towards selecting the single hidden layer BPNN [31]. As mentioned above, the structure of GRA-PSO-BPNN model is determined to be a 3-layer structure of 6-12-1.Step 2: Particle Swarm Initialization.The weights and thresholds in the BPNN model are the parameters to be optimized by the PSO algorithm. At the beginning, a certain number of particle swarm individuals can be randomly generated to represent, and the global optimal position gbest and the individual optimal position pbest of the particle are initialized. The sum of ownership value and threshold value in network structure is the dimensiond of particle swarm individual search space, that is,(15)d=xm+my+m+y.Thed-dimensional vector xit=xi1t,xi2t,…,xidt represents the position vector of the ith particle in the tth generation population.Step 3: Particle Velocity Location Update.According to the following formula, the particle velocity is updated according to the individual optimal solution and the global optimal solution.(16)vt+1=w⋅vt+c1⋅r1⋅pbestt−xt+c2⋅r2⋅gbestt−xt.Among them, pbest (t) and gbest (t) are the individual optimal solution of the ith particle in the tth generation of the d-dimensional space, and the global optimal solution in the tth generation, r1 and r2 are in [0, 1], c1, c2 are learning factors, w is the inertia weight, and the particle position is updated according to the following formula:(17)xt+1=xt+vt+1.The particle group can be evaluated by calculating the fitness value. The smaller the fitness value of the particle group is, the higher the fitness is. The speed and position of each particle are adjusted based on the fitness value. The particle fitness is calculated according to the following formula. ValueF is(18)F=∑i=1Nabsyi−ti,whereN is the number of samples in the test area; abs is the absolute value function; yi is the actual value of the sample i; and ti is the predicted value of the sample i.Step 4: Optimal Population Particles:When the number of population evolution reaches the upper limitT or the iteration error reaches the set precision e, when the algorithm stops, the global optimal solution can be obtained and mapped to the weight and threshold of the BP network.Step 5: GRA-PSO-BPNN Model Training.After determining the GRA-PSO-BPNN weights and thresholds, the training data are input to train the GRA-PSO-BPNN model. The entire model flow is shown in Figure3.Figure 3
GRA-PSO-BPNN model flow chart.


## 3. Case Study

300 groups of tree growth and nutrient dynamic test data were measured in the test area, and the test data were divided into two parts. Among them, 200 groups of experimental data were used as training data to train the stage accurate fertilization model of forest based on neural network, and 100 groups of experimental data were used as the verification set of test fertilization model.BPNN models optimized by different algorithms are compared and verified in the experiment. Tansig function and purelin function are, respectively, assigned as the activation functions of hidden layer and output layer. In the process of adjusting parameters of BPNN model, the accuracy is found higher when the maximum training times of BPNN model reaches about 200 times, while big error and inconsistencies occur when the training times are more than or less than 200 times. Therefore, the max of training times of the model is set to 200, the target error is 0.00001, and the learning rate is 0.1. Even so, the target accuracy is still unachievable, and the training time is about 5∼8 s. Based on the BPNN model, the PSO algorithm was used for optimization. When adjusting the parameters of the PSO-BPNN model, it is found that accuracy was higher when the population size was about 200 and the training time was about 8∼14 s. When the population size is significantly lower than 200, the accuracy becomes lower; when the population size is significantly higher than 200, the accuracy remains the same, but the training time becomes longer. At the same time, the accuracy is higher when the max of training times of the PSO-BPNN model is in the range of about 200 times; the accuracy however becomes lower when it is significantly lower than 200 times. When the training times greatly exceed 200, the accuracy remains the same and the training time becomes longer. Therefore, the maximum training times of the model is set as 200. The population size of the particles in the PSO algorithm is 200, the inertia weight is 1, and the learning factorc1 = c2 = 1.5. The parameters remain unchanged and the training time is about 4∼10 s when the PSO-BP model is improved to GRA-PSO-BPNN model for training. With regard to the basic BPNN, PSO-BPNN, and GRA-PSO-BPNN training models, Figures 4–6 show the prediction situation and error percentage of tree fertilization.Figure 4
Prediction results and error analysis of nitrogen, phosphorus, and potassium fertilizer application rates in BP model.
(a)(b)(c)(d)(e)(f)Figure 5
Prediction results and error analysis of nitrogen, phosphorus, and potassium fertilizer application rates in the PSO-BP model.
(a)(b)(c)(d)(e)(f)Figure 6
Prediction results and error analysis of nitrogen, phosphorus, and potassium fertilizer application rates in GRA-PSO-BP model.
(a)(b)(c)(d)(e)(f)It can be seen from Figures4–6 that the basic BPNN model can only predict the change level of forest fertilization amount and not accurately predict the actual fertilization amount, the accuracy is poor, the error range is about 20%, and the PSO-BPNN model predicts the amount of forest fertilization. The accuracy is obviously improved, the error between predicted fertilization amount and actual fertilization amount is reduced, and the error range is about 10%, which indicates that the optimization of PSO algorithm has a great influence on accurately predicting the amount of forest fertilization; the main influencing factors of determining the amount of forest fertilization by using grey correlation analysis method. After that, the GRA-PSO-BP prediction model of this paper further enhances the ability to predict the amount of forest fertilization. The predicted value of the GRA-PSO-BP neural network model has small fluctuations near the actual value. The model has a good prediction effect. The error between the predicted fertilization amount and the actual fertilization amount is within 5%, which can reflect the different forests. The change in nutrient demand during the growth phase can well guide the staged precision fertilization.The three models were used to predict the nitrogen, phosphorus, and potassium fertilizer application rates of 8 experimental sites and compared with the actual fertilization use. The results are shown in Table5.Table 5
Comparison of measured results with predicted results (kg/hm2).
NumberStand age (t/a)Site indexSActual amount of nitrogen fertilizerBPPercentage of error (%)PSO-BPPercentage of error (%)GRA-PSO-BPPercentage of error (%)15a517.3816.316.55117.671.6617.350.1752218a1984.1992.228.70686.843.14284.320.155321a22218.12009.082217.30.3702218.30.0666423a25202.32197.628203.30.47842010.668455a514.9416.328.41615.080.907414.950.0589618a19119.9130.58.12124.63.8761200.0650721a22216.5199.38.605214.60.8908216.40.0367823a25207222.16.79209.91.402209.11.0215a51.6091.79410.321.6492.451.6110.1245218a199.81110.859.5589.6591.5529.8180.0637321a2228.5530.717.02728.911.27828.40.705423a2539.5843.128.20639.750.425539.880.747155a51.4041.5016.4421.3811.6341.3940.7023618a1914.9113.788.23614.840.515114.920.0377721a2224.9427.459.12524.90.174724.740.8051823a2539.7342.937.44439.840.261939.450.699315a57.416.7449.8747.3970.16937.3880.2956218a19111104.16.643110.60.3362110.20.664321a22155.91625.541155.80.07406155.50.3252423a25237.9253.86.257241.61.5572351.21255a58.0457.4837.5148.0390.078328.0350.1311618a1987.8893.756.26791.213.79187.90.0249721a22172.8183.45.767172.60.081861730.1017823a25281.2261.87.395269.14.289283.20.7244It can be seen from Table5 that the error percentage between the predicted fertilization amount and the actual fertilization use amount of the GRA-PSO-BP prediction model is the smallest and the error is within 5%. The prediction accuracy of the GRA-PSO-BP prediction model is better than other fertilization models. The prediction accuracy is high, the error between the predicted fertilization amount and the actual fertilization usage is small, and the generalization performance is good. The model can fully describe the mapping relationship between input and output. According to the target production requirements of Chinese fir forest biomass, the practical amount of nitrogen, phosphorus, and potassium fertilizers can be predicted, which can guide the staged precision fertilization.On the basis of the above data of actual fertilizer application of nitrogen, phosphorus, and potassium and predicted fertilizer application of nitrogen, phosphorus, and potassium, the paper calculated the mean square error between fertilization prediction results and expected values of different neural network algorithms, different absolute error values, and correlations one by one. The formula for calculating the mean square error is shown in (19). Besides the paper also compared the error and correlation to analyze the traditional BP Neural Network, BP Neural Network that was optimized by PSO algorithm, and the BP Neural Network Prediction Model that was optimized by GRA-PSO algorithm:(19)σ=ε12+ε22+⋯+εn2n.In the formulation,n represents the number of calculated mean square error data, ε1, ε2, …, εn indicates the difference between the actual fertilizer application amount of nitrogen, phosphorus, and potassium and the predicted fertilizer application amount of nitrogen, phosphorus, and potassium and σ means the square error.Tables6–8 show the prediction models which are the traditional BP Neural Network, BP Neural Network that was optimized by PSO algorithm, and the BP Neural Network that was optimized by GRA-PSO algorithm, which predict the correlation and error value between the predicted fertilizer application of nitrogen, phosphorus, and potassium and the actual fertilizer application rate of nitrogen, phosphorus, and potassium.Table 6
Correlation and error value of three models for predicting nitrogen fertilizer.
Fertilization modelMean square errorMaximum absolute errorMinimum absolute errorMean absolute errorCorrelationTraditional BP neural network12.79718.11.0711.020.9881BP neural network optimized by PSO algorithm2.31514.70.141.7970.9997BP neural network optimized by GRA-PSO algorithm0.87872.10.010.4960.9999Table 7
Correlation and error value of three models for predicting phosphorus fertilizer.
Fertilization modelMean square errorMaximum absolute errorMinimum absolute errorMean absolute errorCorrelationTraditional BP neural network2.12533.540.0971.7320.9986BP neural network optimized by PSO algorithm0.15900.360.0230.1210.9999BP neural network optimized by GRA-PSO algorithm0.16990.30.0020.1190.9999Table 8
Correlation and error value of three models for predicting potassium fertilizer.
Fertilization modelMean square errorMaximum absolute errorMinimum absolute errorMean absolute errorCorrelationTraditional BP neural network10.37719.40.5618.2430.9939BP neural network optimized by PSO algorithm4.628712.10.0062.4810.9989BP neural network optimized by GRA-PSO algorithm1.28702.90.010.7940.9999According to the data in above table, the mean square error of nitrogen, phosphorus, and potassium fertilizer of traditional BP neural network are 12.797; 2.1253; and 10.37. In these data, the maximum absolute errors are 18.1; 3.54; and 19.4. The minimum absolute errors are 1.07; 0.097; and 0.561, The average absolute errors are 11.02; 1.732; and 8.243, The correlation numbers are 0.9881; 0.9986; and 0.9939. The mean square error of nitrogen, phosphorus, and potassium fertilizer of BP Neural Network that was optimized by PSO algorithm are 2.3151; 0.1590; and 4.6287. In these data, the maximum absolute errors are 4.7; 0.36; and 12.1. The minimum absolute errors are 0.14; 0.023; and 0.006. The average absolute errors are 1.797; 0.121; and 2.481. The correlation numbers are 0.9997; 0.9999; and 0.9989. The mean square error of nitrogen, phosphorus, and potassium fertilizer of the BP Neural Network that was optimized by GRA-PSO algorithm are 0.8787; 0.1699; and 1.2870. In these data, the maximum absolute errors are 2.1; 0.3; and 2.9. The minimum absolute errors are 0.01; 0.002; and 0.01. The average absolute errors are 0.496; 0.119; and 0.794. The correlation numbers are 0.9999; 0.9999; and 0.9999.Through the error value and correlation between the predicted results of N, P, and K fertilizer amount and the actual fertilization amount, it can be seen that the prediction accuracy of predication mode of BP neural network optimized by GRA-PSO algorithm is the best. However, the 8 groups of experimental data are only a small part of the test set. And the poor data sampling may achieve the best predication accuracy of the BP neural network model optimized by the GRA-PSO algorithm. To assure the best predication accuracy of the the BP neural network model optimized by the GRA-PSO algorithm under the maximum probability, the probability of significant difference between BP neural network model and PSO-BP neural network model and that between the PSO-BP neural network model and the BP neural network prediction model optimized by the GRA-BP algorithm must be calculated. DM test provides a calculation method from the perspective of statistics.The following are the specific methods of DM test.Because three different algorithms predict models in the span ofT = 8, the error sequence of predicted and actual values of BP network prediction model, PSO-BP network prediction model, and GRA-PSO optimized BP network prediction model are calculated on the span of T = 8, respectively:(20)Ea=a1,a2,…,aT,(21)Eb=b1,b2,…,bT,(22)Ec=c1,c2,…,cT.Next, calculate the difference sequenceDab = [dab1, dab2, …,dabT], Dbc = [dbc1, dbc2, …, dbcT], where dabi = ai − bi, dbci = bi − ci.Finally, the mean and standard deviation ofDab and Dab are calculated, and DM statistics are calculated. The formulas are as follows:(23)dab,bcmean=∑i=1Tdab,bciT,(24)dab,bcstd=∑i=1Tdab,bci−dab,bcmean2T−1,(25)DM=dab,bcmeandab,bcstd.The DM test theory holds that the distribution of DM is in accordance with the requirements of the standard normal distribution. Therefore, by querying the confidence values corresponding to DM in the standard normal distribution table, the confidence values of the significant differences between BP network model and PSO-BP network model and between PSO-BP network model and BP network prediction model optimized by GRA-PSO algorithm can be obtained while the predication of the fertilization amount of nitrogen, phosphorus, and potassium is advancing. It is shown in the results that when predicting the fertilization amount of nitrogen, phosphorus, and potassium, the confidence values with significant differences between BP network model and PSO-BP network model and between PSO-BP network model and BP network prediction model optimized by GRA-PSO algorithm are 53%, 68%; 81%, 64%; 61%, and 53%, while due to the smallest prediction errors of BP network prediction model optimized by the GRA-PSO algorithm, it enjoys the highest accuracy and the best effect.The prediction data of the network training are limited, so the different trained fertilization models are timeliness and regionality, while as for this research, because it is on the basis of the annual forest data that was updated by artificial forest and adopted the Analysis method of Grey Relational degree to analyze the key factors of fertilizer application amount in different areas, this can make the precision fertilization come true. After comparing these three prediction algorithms, the result came out that the square error, maximum absolute error, minimum absolute error, and average absolute error numbers of the BP Neural Network that were optimized by GRA-PSO algorithm are small, while the correlation number of it is large, so that the accuracy rate of the BP Neural Network that was optimized by GRA-PSO algorithm is high, and the square error, maximum absolute error, minimum absolute error, and average absolute error numbers of the traditional BP Neural Network and BP Neural Network that were optimized by PSO algorithm are relatively large, but the correlation number is relatively small, which draws the conclusion that the accuracy rate of these two prediction methods is relatively low.


## 4. Conclusions

It is inefficient to determine the amount of fertilization by traditional experience or basic theoretical model, which may also cause environmental pollution, so it is difficult to adapt to the new situation of precision forestry advocated by the contemporary world. As the current application of sensor network is getting more diversified, data acquisition in forest area is easier than before. Through real-time data update, real-time prediction can be carried out using high-tech to build fertilizer model. But different fertilization models will have timeliness and regionality. In different forest ecological areas, it seems to be impossible to accurately make a prediction at present using the definite fertilization model with an eye to different factors that affect the precision fertilization of trees.The traditional methods to determine the amount of fertilizer application are empirical method, fertilizer effect function method, and nutrient balance method. The empirical method means that forest farmers fertilize trees on their willing which is according to their experience and without basis. This method is simple and sometimes feasible, but in most cases, it is not accurate enough, which often results in insufficient or excessive fertilizer application and make it not fit the needs of high yield and efficiency of forestry construction. The fertilizer effect function method needs a lot of experimental data, but the successful rate of the fitting is not high, and a lot of data resources are just wasted, so that a lot of financial wealth, material, and human labors are wasted too. The nutrient balance method often needs to come up with expert experience, because it requires too many parameters and also there is a relative error in the calculation of fertilization. It is important to find an effective and accurate fertilization model. Artificial neural network has a strong ability to solve nonlinear problems. In order to solve the problem of nonlinear precision fertilization, this paper introduces the neural network modeling method and uses BPNN, PSO-BPNN, and GRA-PSO-BPNN to construct three kinds of forest precise fertilization models for comparative verification and result analysis. Through comparison, the accuracy of the BPNN model is low and the error range is about 20%, and it is unstable. Although the PSO-BPNN model has high accuracy, the error range is about 10%, considering that the fertilization model mentioned above is time-sensitive and regional, the factors that may affect the fertilization are different in different time and different regions. In some geographically complex areas, more factors need to be considered to ensure accurate fertilization. Since it is difficult for the PSO-BPNN model to identify the key influencing factors of forest fertilizer application in different regions, the PSO-BPNN model has poor scalability. However, the key factors influencing forest fertilizer application in different regions can be determined by grey relational analysis. The constructed GRA-PSO-BPNN forest tree precision fertilization model not only has high accuracy, the error range is about 5%, and it has a faster calculation speed, but also applicable to the determination of forest fertilizer application in different regions and can constantly update the forest fertilizer model with good expansibility.


## 5. Discussion

With the in-depth development of modern forestry, precision forestry is receiving unprecedented attention. The emergence of forest precise fertilization model will greatly promote the forestry production. Forest economy plays an important role in the society improvement, and to improve the forest economy, only the development of artificial fast-growing and high-yielding forest can be depended, and the effective way to improve the productivity of artificial forest is just precision fertilization. The prediction of precision fertilization is a scientific and innovative technology which can achieve five goals: ensuring high quality, increasing the production, improving soil quality and ecology, and high efficiency; taking the factors which affect the amount of fertilizer application as the basis and making the fertilizer application rate and fertilization target clear before implementing it; linking the fertilization with economic benefit; achieving the quantitative fertilization of forest trees under the realization of target benefit, making a prediction of the benefit that earned by the fertilization, and according to the investigation of basic data and soil chemical analysis in forest area, so that the ratio of tree fertilization under the target benefit can be calculated; making forest fertilization become an important safeguard measure to improve the quality and yield of commercial forest. At the same time, as a direction of precision forestry technology, precision fertilization model will attract more and more attention of forestry science and technology personnel and forestry production managers. In this paper, the GRA-PSO-BPNN prediction model of forest periodic precise fertilization is of far-reaching significance to solve the problem of forest fertilization that has been troubling people for a long time, which can be applied by forest-related workers in forestry ecological zones where conditions are basically the same. If the generalization is carried out in a larger and more complex area to improve the generalization ability of the model, it is necessary to consider other factors affecting the amount of forest fertilization, such as soil type, organic matter, pH, and other factors. In some cases, the GRA-PSO-BPNN precise fertilization model for trees may need to be updated based on the experience of experts, which is what needs to be further studied.



---

*Source: 1356096-2020-08-27.xml*