Use of linear programming models in experimentation with plant nutrients

Uso de modelos de programação linear aplicado à experimentação nutricional de plantas

Use of linear programming models in experimentation with plant nutrients

Uso de modelos de programação linear aplicado à experimentação nutricional de plantas

Mauro Brino Garcia^I; Lucas Rezende Gomide^II

^IForest Engineer, MSc candidate in Soil Fertility and Plant Nutrition – Universidade Federal de Lavras/UFLA – Departamento de Ciência do Solo – Cx.P. 3037 – 37200-00 – Lavras, MG, Brasil – [e-mail]

^IIForest Engineer, Professor PhD in Forest Engineering – Universidade Federal de Lavras/UFLA – Departamento de Ciências Florestais – Cx.P. 3037 – 37200-00 – Lavras, MG, Brasil – [e-mail]

ABSTRACT

Nutrition is an important issue of plant cultivation and experimentation with plant nutrients is a supporting tool for agriculture. However, use of high purity grade reagents as nutrient sources can be expensive and increases the cost of an experiment. The objective of this study was to minimize the acquisition cost of high purity grade reagents in experiments on plant nutrient deficiency by using the missing element technique through linear programming models, and to generate recommendation tables for preparation of culture solutions, as well as to quantify gains through a simulated experiment. Two linear programming models were formulated containing concentration constraints for each nutrient in the culture solution. Model A was based on 16 reagents for preparation of the culture solution, while model B was based on 27 reagents, looking to increase choice options. Results showed that both models minimized the acquisition cost of reagents, allowing a 9.03% reduction in model A and a 25.98% reduction in model B. The missing sulfur treatment proved the most costly for reagent acquisition while the missing nitrogen treatment proved the least costly. It was concluded that the formulated models were capable of reducing acquisition costs of reagents, yet the recommendations generated by them should be tested and checked for practical viability.

Key words: Operating research, optimization, missing element.

RESUMO

A nutrição vegetal é um importante aspecto no cultivo de espécies, sendo a experimentação nutricional uma ferramenta de suporte à agricultura. Entretanto, o uso de reagentes p.a. como fontes de nutrientes é caro e aumenta os custos do experimento. Assim, objetivou-se, com este trabalho, minimizar o custo de aquisição de reagentes p.a. em experimentos de deficiência nutricional de plantas, considerando a técnica do elemento faltante, por meio de modelos de programação linear, gerar tabelas de recomendação para a elaboração de soluções de cultivo e quantificar os ganhos por um experimento simulado. Dois modelos de programação linear foram formulados, contendo restrições de concentração para cada nutriente na solução de cultivo. O modelo A baseou-se na utilização de 16 reagentes na elaboração da solução de cultivo, já, o modelo B empregou 27, com o intuito de aumentar as opções de escolha. Os resultados mostraram que os dois modelos minimizaram o custo na aquisição de reagentes, com redução de 9,03% (modelo A) e 25,98% (modelo B). O tratamento de omissão de Enxofre foi o mais oneroso para a aquisição de reagentes e o de omissão de Nitrogênio o menos oneroso. Conclui-se que os modelos formulados foram capazes de reduzir os custos na aquisição dos reagentes; porém deve-se testar as recomendações geradas pelos mesmos e verificar sua viabilidade prática.

Palavras-chave: Pesquisa operacional, otimização, elemento faltante.

1 INTRODUCTION

Soil is a biologically active layer of material resulting from complex transformations which in turn involve rock and mineral weathering, nutrient cycling, as well as biomass production and decomposition. It is thus the primary medium for plant growth (NOVAIS et al., 2007). A way of maintaining nutrient levels and soil stability is to use fertilization. However, suitable management of fertilization involves defining not only the correct dosages and sources of nutrients but also the best seasons and methods to apply liming and fertilizers to the soil (OLIVEIRA et al., 1991).

The learning process of soil-plant relationships calls for experiments with which to identify nutrient deficiencies in the target species (OLIVEIRA et al., 1991). The most widely used method in experiments on nutrient deficiency is the missing element technique, based on Liebig's law of the minimum (1840). For this experiment, high purity grade reagents (pure for analysis) were used as nutrient sources, that way ensuring maximum purity and precision although their prohibitive cost obstructs research.

Studies on mineral nutrition of plants are widely explored in literature (MAFFEIS et al., 2000; MARQUES et al., 2004; SILVA et al., 2005; WALLAU et al., 2008). However, defining optimal nutrient formulations in connection with the acquisition cost of reagents is a subject not sufficiently explored by researchers, despite being critical to the reduction of such costs without affecting the development and conduction of an experiment.

Use of tools such as linear programming (LP) to assist the decision-making can contribute toward finding optimal responses (BERGER et al., 2003; RODRIGUEZ et al., 1985; SILVA et al., 2003), whenever possible, helping researchers solve complex problems. A line of research in the field of mathematical programming, LP is intended to maximize or minimize a function subject to linear constraints that delimit a feasible region (GOLDBARG; LUNA, 2005).

The objective of this study was to apply LP models to experimentation with plant nutrients in an attempt to minimize the acquisition cost of high purity grade reagents in experiments involving nutrient deficiency by using the missing element technique, and also to generate recommendation tables for preparation of culture solutions as well as to quantify gains through a simulated experiment.

2 MATERIAL AND METHODS

The base nutrient formulation to be optimized in this work was directed at experiments on nutrient deficiency (missing element technique), in which one nutrient is eliminated per treatment. The method generates a large number of combinations between reagents until the desired concentrations are attained. The LP mathematical model was thus based on the classic mixture problem (GOLDBARG; LUNA, 2005). The model developed for the relevant problem is provided as follows, where (1) is objective function, (2) are constraints on nutrient concentration and (3) is non-negativity condition:

subject to constraints

where,

n - number of reagents used for providing nutrients tested in the experiment with the missing element technique.

K - nutrients tested in the experiment with the missing element technique, including N, P, K, Ca, Mg, S, Fe, B, Mn, Zn, Cu and Mo;

P_i - price of each reagent, expressed as Reais per gram (R$.g^-1);

x_i - decision variable representing the quantity of reagent i in gram (g) to be used in formulating one liter of culture solution;

R_ki - quantity of nutrient k available and present in reagent i, expressed as mass ratio of atomic weight of atoms of the nutrient present in the reagent to molecular mass of reagent;

C_k - desired concentration of nutrient k expressed as grams per liter of culture solution (g.l^-1) (PRADO, 2010).

Table 1 provides the price (R$.g^-1) of every reagent being analyzed, relevant molecular formulations and which reagents are being tested in models A and B. Each reagent was correlated to a decision variable (x_i) representing the quantity of reagent k to be used in formulating one liter of culture solution. Reagent prices came from a quote obtained from a company selling laboratory products, on November 18, 2010.

The adopted reference table was obtained from Prado (2010), which provides mg.l^-1 concentrations of each nutrient (C_k) in the culture solution, where: N (210.1011), P (30.9737), K (234.5898), Ca (200.3900), Mg (48.6100), S (64.1320), Fe (5.0103), B (0.5000), Mn (0.5025), Zn (0.0480), Cu (0.0189) and Mo (0.0107). With this information at hand, two LP models were developed, the first being termed model A and comprising 16 reagents, the second being termed model B and comprising 27 reagents in an attempt to increase model choice options.

Model constraints were formulated to define the required concentration of each nutrient in the culture solution, following the order: N; P; K; Ca; Mg; S; Fe; B; Mn; Zn; Cu and Mo. To eliminate one nutrient and form a new nutrient-omitting treatment (missing element), it was necessary to reset parameter C_k and thus form a new treatment scenario. The scenarios created in this study were analogous with experimentation treatments on plant nutrients. For solving the linear programming models, the commercial solver ILOG CPLEX 11 was used.

An economic analysis of acquisition cost of reagents was carried out, comparing the reagent quantity recommended by Prado (2010) with the quantity found by models A and B in this work. For that, an experiment was simulated which lasted four months and whose experimental design consisted of 13 treatments (as used in the experiment on nutrient deficiency with the missing element technique), with 8 replicates and 4 vases per replicate. Each vase had a capacity of 3 liters and the culture solution was replaced every 7 days, according to variations in concentration due to plant absorption (TAIZ; ZEIGER, 2004). The strategy was to progress with the experiment, computing its cost with the culture solution.

3 RESULTS AND DISCUSSION

Model A selected virtually every reagent suggested by Prado (2010), except potassium nitrate (x₂) and calcium nitrate (x₃), using larger quantities of calcium chloride (x₆), ammonium nitrate (x₈) and magnesium sulfate (x₄). It was noted that potassium nitrate (x₂) and calcium nitrate (x₃) were replaced with ammonium nitrate (x₈) which, despite its higher unit cost, has a greater mass ratio of nitrogen and thus meets the nitrogen requirement at a lower cost. The quantities used of micronutrient-supplying reagents remained constant, as there are no alternative sources of these for replacement.

With a larger number of reagent options, model B selected the variables potassium phosphate (x₁), potassium chloride (x₅), calcium chloride (x₆), monoammonium phosphate (x₇), ammonium sulfate (x₉), boric acid (x₁₁), manganese chloride (x₁₂), zinc chloride (x₁₃), copper chloride (x₁₄), molybdic acid (x₁₅) iron chloride (x₁₆), sodium nitrate (x₁₈), potassium sulfate (x₁₉), sodium sulfate (x₂₀), magnesium chloride (x₂₁), copper sulfate (x₂₃), zinc sulfate (x₂₄), manganese sulfate (x₂₅), iron sulfate (x₂₆), in which the most generally used were sodium nitrate (x₁₈), calcium chloride (x₆) and magnesium chloride (x₂₁). Following results, the potassium nitrate (x₂), calcium nitrate (x₃) and ammonium nitrate (x₈) used by Prado (2010) and in model A were replaced here with sodium nitrate (x₁₈) due to low cost.

In both models A and B, the main source of Ca was calcium chloride (x₆). In model B, magnesium chloride (x₂₁) fully replaced magnesium sulfate (x₄) and magnesium nitrate (x₁₀) due to having only magnesium as nutrient and due to low cost.

Tables 2 and 3 provide recommendations for treatment composition (missing element), and values of stock solutions for each reagent selected in models A and B respectively.

Based on scaled quantities of each reagent, required costs per treatment were then generated for preparing one liter of culture solution, following the experiment with the missing element technique, as illustrated in Table 4.

By using the LP model clearly helped reduce the acquisition cost of reagents, with model A allowing savings of 9.03%, using the same range of reagents suggested by Prado (2010), only with varying quantities of each reagent. In the treatment where N was omitted, no reduction was noted due to the poor choice of reagents without N in their composition. The treatment where K was omitted had a 12.55% cost reduction, mainly due to replacement of the Ca source, in which calcium chloride (x₆) replaced calcium nitrate (x₃). In addition, nine treatments attained a near to 10.50% reduction.

Model B promoted a 25.98% cost reduction against what is suggested by Prado (2010) and a 18.63% reduction in relation to model A, and this logic is explained by the fact that model B works with 11 more choices of reagents than model A. The treatment where S was omitted promoted the greatest cost reduction (37.19%) due to replacement of magnesium nitrate (x₁₀) with magnesium chloride (x₂₁) and to use of sodium nitrate (x₁₈) as source of N. Around 8 treatments attained savings of 25.90%, and the smallest reduction was 14.95%. However, the 14.95% reduction in model B was greater than the maximum attained in model A (12.55%).

A similar behavior was noted in all three cases regarding the costs of treatments omitting N and S, the first generating lower costs and the second generating higher costs among treatments. N-supplying reagents entail high cost for formulation of nutrient solutions, which explains the low cost of the treatment omitting N, when such reagents are not present. The high cost of the treatment omitting S is explained by the large number of reagents with this nutrient in their composition, 8 in total, thus forcing use of a limited number of reagents which, in turn, constitute a high-costing group.

Table 5 provides costs of reagent acquisition for conducting the simulated experiment, showing that values could be as high as around R$2,000.00, a significant sum considering the intended purpose of the research.

The cost reduction promoted by the LP mathematical models proposed in this work reached over R$ 500.00 (model B) per experiment. These savings could be the determining factor for the viability or not of a research project, and help shift the focus of research onto other priorities such as sampling intensity, increasing the number of vases per replicate (plot) or the number of replicates, or else extend the experiment duration (weeks) (Table 6).

The savings of R$ 332.75 in model A could allow addition of an extra replicate in the experimental design or else prolong the duration of the experiment for another three weeks. As for model B, the generated savings of R$ 510.08 could allow addition of an extra vase per replicate (plot), or two replicates, or else 5 weeks conducting the experiment. Here, savings generated by models A or B and applied to the experiment could bring benefits, for instance, improve experimental accuracy with increased sample intensity.

Considering the savings generated by model B only, and based on 30 experiments per year, a reduction of R$15,302.40/year in expenditure would be achieved, something real and practical in most laboratories of research institutions and universities. This statement is intended to serve as a reference for planning of laboratories exploring mineral nutrition of plants, where the number of experiments being conducted in a short time span is expressive.

That said, although linear programming has been known for a while, it has not been widely used for assisting with decision-making in laboratory planning, yet it could potentially make a difference where costs are concerned, bringing benefits to research.

4 CONCLUSIONS

The models formulated via LP proved to be effective in minimizing the acquisition cost of reagents for generating nutrient solutions in experimentation with plant nutrients, observing the desired concentrations of each nutrient. However, real trials should be conducted to check the recommendation suggested by the mathematical models for viability. The cost reduction was 9.03% in model A and 25.98% in model B, based on the reference table provided by Prado (2010).

The recommendation tables for formulating experiments on nutrient deficiency with the missing element technique can be easily generated based on results obtained in models A and B

5 REFERENCES

(received: February 9, 2011; accepted: December 21, 2012)

• BERGER, R.; CARNIERI, C.; LACOWICZ, P. G.; BRASIL, A. A. Minimização de custos de transporte florestal com a utilização da programação linear. Floresta, Curitiba, v. 33, n. 1, p. 53-62, 2003.
• GOLDBARG, M.; LUNA, H. Otimização combinatória e programação linear: modelos e algoritmos. 2. ed. Rio de Janeiro: Elsevier, 2005. 518 p.
• LIEBIG, J. von. Chemistry and its application to agriculture and physiology London: Tauylor and Walton, 1840.
• MAFFEIS, A. R.; SILVEIRA, R. L. V. de A.; BRITO, J. O. Reflexos das deficiências de macronutrientes e boro no crescimento de plantas, produção e qualidade de óleo essencial em Eucalyptus citriodora Scientia Forestalis, Piracicaba, n. 57, p. 87-98, jun. 2000.
• MARQUES, T. C. L. L. de S. e M.; CARVALHO, J. G. de; LACERDA, M. P. C.; MOTA, P. E. F. da. Crescimento inicial do Paricá (Schizolobium amazonicum) sob omissão de nutrientes e de sódio em solução nutritiva. Cerne, Lavras, v. 10, n. 2, p. 184-195, jul./dez. 2004.
• NOVAIS, R. F.; ALVAREZ, V. H. V.; BARROS, N. F.; FONTES, R. L. F.; CANTARUTTI, R. B.; NEVES, J. C. L. Fertilidade do solo Viçosa, MG: SBCS, 2007. 1017 p.
• OLIVEIRA, A. J.; GARRIDO, W. E.; ARAÚJO, J. D.; LOURENÇO, S. (Coord.). Métodos de pesquisa em fertilidade do solo Brasília: EMBRAPA-SEA, 1991. 392 p.
• PRADO, R. de M. Diagnose de deficiência nutricional em culturas Jaboticabal: UNESP, 2010. Disponível em: <http://www.nutricaodeplantas.agr.br/site/downloads/unesp_jaboticabal/trabalho_pratico.pdf>. Acesso em: 25 abr. 2011.
• RODRIGUEZ, L. C. E.; LIMA, A. B. H. P. M. de. A utilização da programação linear na determinação de uma estratégia ótima de reforma de um talhão florestal. IPEF, Piracicaba, n. 31, p. 47-53, dez. 1985.
• SILVA, E. de B.; GONÇALVES, N. P.; PINHO, P. J. de. Limitações nutricionais para crescimento de mudas de umbuzeiro em Latossolo Vermelho distrófico no Norte de Minas. Acta Scientiarum, Maringá, v. 27, n. 1, p. 55-59, jan./mar. 2005.
• SILVA, G. F. da; LEITE, H. G.; SILVA, M. L. da; RODRIGUES, F. L.; SANTOS, H. do N. Problemas com o uso de programação linear com posterior arredondamento da solução ótima, em regulação florestal. Revista Árvore, Viçosa, v. 27, n. 5, p. 677-688, set./out. 2003.
• TAIZ, L.; ZEIGER, E. Fisiologia vegetal 3. ed. Porto Alegre: Artmed, 2004. 719 p.
• WALLAU, R. L. R. de; BORGES, A. R.; ALMEIDA, D. R. de; CAMARGOS, S. L. Sintomas de deficiências nutricionais em mudas de mogno cultivadas em solução nutritiva. Cerne, Lavras, v. 14, n. 4, p. 304-310, out./dez. 2008.

Publication Dates

• Publication in this collection
  16 July 2013
• Date of issue
  June 2013

History

• Received
  09 Feb 2011
• Accepted
  21 Dec 2012