versão On-line ISSN 0718-9516
J. Soil Sci. Plant Nutr. vol.12 no.1 Temuco 2012
Journal of Soil Science and Plant Nutrition, 2012, 12 (1), 87-97
Use of nonlinear regression and nonlinear mathematical programming in the formulation of substrate mixtures for soil-less culture - a review -
J.A. Díaz-García1* and R. Requejo-López2
1Universidad Autónoma Agraria Antonio Narro, Department of Statistics and Computation, 25315 Buenavista, Saltillo, Coahuila, México. 'Corresponding author: email@example.com
2Universidad Autónoma Agraria Antonio Narro, Department of Soil Science, 25315 Buenavista, Saltillo, Coahuila, México.
With respect to four inorganic and organic substrates - tezontle (volcanic rock), perlite, peat and coir - we studied the relevant physical and chemical properties of mixtures of two of these substrates. A nonlinear regression model was used to analyse the physical properties of each possible mixture. Finally, a methodology is proposed for formulating mixtures of substrates as a deterministic nonlinear mathematical programming problem.
Key Words: Coir; Peat; Nonlinear regression, Nonlinear optimisation; Substrate mixtures; Tezontle.
Soil serves as a support for plants and is the medium from which they obtain water, mineral substances and oxygen, which are essential for growth and vegetative development. An ideal soil has a porosity and particle disposition that enables root penetration while retaining appropriate amounts of water and air [Foth and Ellis (1988), Hablin et al. (1999), and Russell (2008)].
However, such soils may not be available, and then artificial soils may be used. This is the case of most flowerpot and containerised crops. Such plants have limited root growth, but they still need nutrients, air and water. For this reason, artificial substrates are sometimes required. These are often obtained by mixtures of various products, and the combination of the characteristics of each component provides the overall (and hopefully, optimal) characteristics of the substrate Resh (2002a), [Resh (2002b) and Russell (2008)].
The apparent density of a substrate should be low, so that the roots can penetrate it easily, and at the same time, to minimise the weight of the containerised crops. Organic matter has high porosity, ionic interchange and water retention capabilities, together with low density. All these properties are recommended for a substrate, which explains why most artificial substrates are essentially composed of organic matter; the remaining components are usually constituted of natural minerals or of artificial minerals such as tezontle, sand, perlite or vermiculite. These mineral products are very porous, with a high real density and a very low apparent density. In general, the grain size of an artificial substrate is higher than that of the soil, which facilitates air infiltration; at the same time, however, this property reduces water retention. This consideration explains why mixtures are usually composed of both organic and mineral substances, to achieve an equilibrium between air infiltration and water retention Resh (2002a), [Resh (2002b) and Russell (2008)].
The chemical properties of a substrate are important because they are the main factors responsible for the availability of nutrients. However, new crops are increasingly being constituted in such a way that the physical properties of a substrate prevail as its desired characteristics, with all the additional plant necessities being provided by nutritive substances, which are permanently available, being recycled to the plant by the irrigation system [Hablin et al. (1999), Resh (2002b) and Russell (2008)].
In many countries, substrates are still obtained by traditional trial and error methods; mixtures are formulated by arbitrarily establishing the proportions of the respective components, and the resulting mixture is then accepted if it improves the overall crop yield. The resulting mixture is not necessarily the optimal one, as the method does not explore all the possible combinations of the available components. The time employed in evaluating the mixtures, the cost of the materials used and the search for optimal resources are among the aspects contributing to the complexity of substrate formulation, given the number of variables involved in the problem Zamora Morales et al. 2005.
Linear programming is a commonly-used empirical technique for formulating substrate mixtures used in a soli-less culture, and generally provides good approximations (Burés et al. (1988), Burés (1997) and Zamora Morales et al. (2005), among many others). However, most authors accept that the constraints imposed in the linear programming are not completely suitable, because they are established over design variables which do not present a linear behaviour.
This paper proposes a general methodology, based on nonlinear regression models and nonlinear mathematical programming, for obtaining substrat mixtures for soil-less culture comprised of two inorganic or organic substrates, to be used for protected and containerised crops. For this purpose, we study the physical and chemical properties of four inorganic and organic substrates - tezontle (volcanic rock), perlite, peat and coir - together with those of the mixtures comprising different proportions of two of these four substrates, see Section 2.
Subsequently, the physical properties of the mixtures are modelled by nonlinear regression models, according to the respective substrate proportions in each mixture, see Section 3. These models are then used to replace the linear restrictions corresponding to the physical properties expressed in the linear mathematical programming problem; thus, we obtain a nonlinear mathematical programming solution. The corresponding linear and nonlinear problems are solved, compared and analysed using the commercial software LINGO.
2. Physical and chemical properties of substrates and mixtures
Four substrates were used in the mixture formulation tests: perlite, tezontle, peat and coir. Following the methodology proposed by Pire and Pereira (2003), Gabriel et al. (1993) and North Carolina State University (2007), we determined, in triple analysis, the following physical and chemical properties: total porous space (TPS), air porosity (AP), humidity retention (HR), apparent density (AD), real density (RD), ashes(A), organic matter (OM), pH and electric conductivity (EC), which are shown in Table 1.
We also formulated six mixtures (denoted by a M1 - M6 and defined in Tables 2 and 3) composed of two substrates with eleven different proportions. Using the above-mentioned procedure, their physical and chemical properties to be used in the formulation of mixtures of substrates were determined (Burés et al., 1988); see Tables 2 and 3.
3. Modelling the physical and chemical properties of the mixtures
As remarked above, in general the physical and chemical properties of substrates (and their mixtures) do not present a linear behaviour. In this section, we study the functional behaviour of the physical properties TPS, HR and OM in the mixtures M1 - M6 1.
The results were obtained by the implementation of a program in S-PLUS (Lam (2011)) to fit the different nonlinear regression models. The least squares method was applied to the following fitting models: linear, polynomial, reciprocal, logarithmic and exponential; for which the particular explicit forms are:
â¢ Linear regression model
â¢Â Polynomial regression model
â¢Â Reciprocal regression model
â¢Â Logarithmic regression model
â¢Â Exponential regression model
The best model was selected by means of the minimum mean square error criterion (Graybill and Iyer (1994)).
In general, the models fitted take the form
where Y is the mixture variable to be modelled (TPS, HR or OM);ƒ (X1,X2) is the functional form of the model to be fitted;X, X2 denote the proportions of two of the substrates(selected from perlite, tezontle, peat and coir) making up the mixture, where O≤ Xk ≤1, k = 1,2, and ξ is the error. Note also that X1, X2 satisfy X1+ X2 =1 which means that X1 and X2 are linearly dependent, moreover, X1 = 1 - Χ2 or X2 = 1 - X1 . A priori, all that is necessay is to fit a model in terms of X1 (or X2). However, in this study we wish to obtain a mathematical model for computing the physical property Y of the mixture, when the proportions and the individual physical properties of the corresponding component of the substrate are known. This model is then used in the nonlinear optimisation program formulation. Strictly speaking, we do not estimate the model for any possible value of X1 and X2, but rather, in the subspace defined by the line X1 + X2 = 1. A possible drawback is the fact that few experimental points (eleven) are used for fitting the TPS, HR and OM characteristic models. Fortunately, the experimental data present a smooth behaviour, and, most importantly, they provide information over the whole domain of X1 and X2 (specifically, the domain is the line X1 + X2 = 1 in the first octant of R3). Moreover, it is important to take into account that models of the form Y = p0 + p1z1 + p2z2 +ξ (with p0 ≠ 0), or similar should not be considered, because when χ1 = X2 = 0 we obtain Y = p0 , but the correct result is Y = 0. Finally, the most complicated aspect of this study is to propose the same functional form for each characteristic (TPS, HR and OM) and for all the possible mixtures obtained.
As shown in Tables 2-3, the variable OM presents a linear behaviour. Then, we just need to determine the behaviour of TPS and HR. In both cases, a well fitted model for TPS and HR for the six mixtures has the form
where Zk= dxk, k = 1,2, and d denotes the individual TPS or HR mixture component of the k -th substrate, respectively.
Table 4. Estimators of the model parameters and the mean square of the error for Y = TPS
Table 5. Estimators of the model parameters and the mean square of the error for Y = HR
As an example, we include in Figure 1 the graphs of the estimated model and the experimental points for the TPS characteristic of mixture M
Figure 1.Graphs of the estimated model and the observed TPS data for the mixture M1.
4. Formulation of mixtures by nonlinear mathematical programming
We consider the hypothesis proposed by Zamora Morales et al. (2005) as the basis for our formulation of the nonlinear mathematical programming problem. Since the mixture would essentially be used as a support, where the plant receives the nutrients via nutritious solutions, we only consider the physical properties as the design variables, except the variable OM, which is an important representative of the substrate's chemical properties and is therefore incorporated into the general model. In addition, we wish to achieve optimum control of the chemical properties of the nutritious solutions.
Taking these considerations into account, we study the following design variables: the balance of volumes, OM, TPS and HR. The minimising objective function is the cost of the mixture, which is studied via linear functions of the substrate proportions in the mixtures and their individual costs. A basic problem type takes the form (2), where Ck is the cost in dollars per litre of the k -th substrate (perlite: 0.073, tezontle: 0.218, peat: 0.21 and coir: 0.106); Xk denotes the proportion in litres of the k -th substrate; mk is the OM percentage per litre of the k -th substrate; ek is the TPS percentage per litre of the k -th substrate and Ck is the HR percentage per litre of the k -th substrate. Now, OM and OM. are limits or levels of the upper and lower restrictions for OM, respectively. In this study, we consider the limits (0-45) and (46-90); similarly, we define the upper and lower limits for TPS and HR, establishing the following two levels: (55-70) and (71-86) for TPS, and (45-60) and (61-76) for HR. These intervals or levels of restriction, for each of the variables OM, TPS and HR, were established following the recommendations of Abad et al. (1988), for substrates used in a soil-less culture and considering labeling specifications for commercial sustrates. Finally, the respective were obtained from Table 5.
Potentially, if we combined the four substrates in the two-substrate mixtures 4C2 with all the possible combinations of the two upper and lower levels for OM, TPS and HR (23) we would have 48=6x8 possible mixtures. However, Table 6 shows only the mixtures selected by linear or nonlinear mathematical programming. In other words, mixtures for which neither a linear nor a nonlinear optimisation program present a feasible region (i.e. the corresponding optimisation program has no solution) have been removed from Table 6, leaving only the substrate mixtures for soil-less culture with a feasible solution, as either a linear or a nonlinear optimisation program. In cases where the mixture of substrates has no feasible solution in the linear or nonlinear optimisation program, this is noted as NF (there is no feasible solution to the optimisation problem). In addition, the restrictions to the coded variables OM, TPS and HR are coded by assigning "1" to the range with the smaller values and "2" to the higher range.
The results were computed using the commercial software Hyper LINGO/PC, release 6.0 (Winston ,1995. The default optimisation method used by LINGO to solve the nonlinear mathematical programming is that of Generalised Reduced Gradient (GRG) (Bazaraa et al. (2006). As can be seen, the linear programming reduces the total of mixtures by 66.66%, while the nonlinear mathematical programming reduces them by 68.75%; thus, both methods reduce the number of mixtures by very similar amounts. However, note the following two aspects:
Table 6. Mixtures of substrates via nonlinear mathematical programming
a)Â The mixture substrate proportions assigned by the two techniques can vary significantly, as can be seen, for example, in peat-tezontle with restriction levels 2 1 2.
b)Â Linear mathematical programming is a useful and widely used tool in mixture problems, however, in general the nonlinear programming model achieves a better mixture formulation. But, if a user insists in using a linear solution, he/she could face two situations: One, is that the optimal linear solution might not be feasible under the nonlinear model, and second, that the optimal linear solution might be suboptimal with the respect to the nonlinear solution. Thus, from a practical point of view, the interpretation of these results would be that on occasions non-feasible substrate mixtures used in a soil-less culture are tested or produced, which results in a waste of research time and/or money. Worse still, when linear programming leads us to reject a mixture that would have been accepted by the nonlinear method, we might be rejecting the best possible mixture.
The aim of this study is to propose a general methodology which can be applied to various situations, with only simple modifications. For example, a mixture formulation in the absence of nutritive substances, or one comprised of substrates that are rich in nutrients and artificial nutrition, could be studied by modifying nonlinear mathematical programming, through restrictions on chemical properties and/or via alternative objective functions.
Alternatively, the presence in the mixture of certain undesirable (desirable) ions or cations could be minimised (maximised). This methodology could be extended to the simultaneous optimisation of several objective functions; for example, to minimise the cost and to maximise the amount of nitrogen in a mixture, which implies the use of nonlinear multiobjective optimisation methods (Steuer, 1986).
Note that the procedure described in Section 3 can be used for modelling various kinds of physical and chemical properties of the mixtures. In addition, the nonlinear mathematical programming approach can be adapted by adding new restrictions and/or objective functions for a particular characteristic.
Finally, note that there are few experimental points and a large number of parameters. Nevertheless, on the one hand, we have a sample of equally-spaced experimental points throughout the domain of the variables X1 and X2 (i.e. the line X1 + X2 = 1), and the mathematical behaviour considered is a smoother of these variables. Furthermore, our aim is to obtain the best mathematical approximation possible to these points and thus obtain the best results.
A potential user of this methodology (e.g. the commercial producer of a protected crop) could consider using a greater number of experimental points, taking into account the higher initial investment that would be required.
The authors wish to express their gratitude to the editor and the referees for their helpful comments and suggestions. This work was partially supported by the research project CONACYT-45974-F, Mexico.
1 The interest in this particular modelling is detailed in Section 4.
Abad Berjón, M. and Noguera Murray, P. 1998. Sustrato para el cultivo sin suelo y feriirrigación, Chapter 8. In, Fertirrigación, Cultivos hortícolas y ornamentales, Ed. Cadahia López, C. Ediciones Mundi-Orensa, Bacelona. [ Links ]
Arthanari, T. S., Dodge Y. 1981. Mathematical Programming in Statistics. John Wiley & Sons, New York. [ Links ]
Bazaraa, M. S., Sherali, H. D., Shetty, C. M. 2006. Nonlinear Programming: Theory and Algorithms, 3rd Edition. Wiley-Interscience. [ Links ]
Burés, S. (1997). Sustratos. Agrotecnia. Madrid, España. 339 pp. [ Links ]
Burés S., Martínez, F. X., Llorca, M. 1988. Preliminary study of the application of parametric linear programming in formulation of substrate mixtures. Acta Horticulturae 221, 141-152. [ Links ]
Foth, H.D., Ellis, B.G. 1988. Soil fertility. John Wiley & S. USA. [ Links ]
Hablin, J. L., Beaton, J. D., Tisdale, S. L., Nelson, W. 1999. Soil fertility and fertilizers. An introduction to nutrient management. Prentice Hall, New Jersey. USA [ Links ]
Gabriel, R., Keirsbulk, W., Engels, H. 1993. A rapid method for the determination of physical properties of growing media. Acta Horticulturae 342, 243-248. [ Links ]
Graybill, F. A., Iyer, H. K. 1994. Regression Analysis. Duxbury Pr. USA. [ Links ]
Pire, R., Pereira, A. 2003. Propiedades físicas de componentes de sustratos de uso comÃºn en la horticultura del estado de Lara, Venezuela. Propuesta metodológica. Bioagro 15(1),1316-1361. [ Links ]
Rao, S. S. 1979. Optimization Theory and Applications. Wiley Eastern Limited. [ Links ]
Resh, H.M. 2002a. Hydroponic food production. Newconcept Press Inc. New Jersey. USA. [ Links ]
Resh, H.M. 2002b. Hydroponics. Questions and answers for successful growing. Newconcept Press Inc. New Jersey. USA. [ Links ]
Russell, E.J. 2008. Soil conditions and plant growth. Biotech books. USA. [ Links ]
Steuer, R. E. 1986. Multiple criteria optimization: Theory, computation and applications. John Wiley, New York. [ Links ]
Winston, W. L. 1995. Introduction to mathematical programming: Applications and algorithms. Dux-bury Press. [ Links ]
Zamora Morales, B. P., Galvis Spínola, A., Volke Haller, V. H., Sánchez García, P., Espinosa Victoria, D. 2005. Formulación de mezclas de sustratos mediante programación lineal. Interciencia: Revista de ciencia y tecnología de América, 30(6), 365-369. [ Links ]