Next Article in Journal
Photo-Fenton Catalyzed by Cu2O/Al2O3: Bisphenol (BPA) Mineralization Driven by UV and Visible Light
Previous Article in Journal
Performance Analysis of an Eductor-Based Membrane Distillation Unit
Previous Article in Special Issue
Performance Management of Natural Resources: A Systematic Review and Conceptual Framework for China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Linear Programming Model for Operational Optimization of Agricultural Activity Considering a Hydroclimatic Forecast—Case Studies for Western Bahia, Brazil

by
Igor Boninsenha
1,*,
Everardo Chartuni Mantovani
1,
Marcos Heil Costa
1 and
Aziz Galvão da Silva Júnior
2
1
Department of Agricultural Engineering, The Federal University of Viçosa, Viçosa 36570-900, MG, Brazil
2
Department of Rural Economy, The Federal University of Viçosa, Viçosa 36570-900, MG, Brazil
*
Author to whom correspondence should be addressed.
Water 2022, 14(22), 3625; https://doi.org/10.3390/w14223625
Submission received: 12 May 2022 / Revised: 28 October 2022 / Accepted: 31 October 2022 / Published: 10 November 2022
(This article belongs to the Special Issue Optimal Utilization and Management of Natural Resources)

Abstract

:
The water crisis is a reality in Western Bahia. In this region, a hydroclimatic monitoring system capable of providing water availability information in advance for water users was implemented objectively to avoid water conflicts. In this study, we proposed the integration between the hydroclimatic monitoring system and a linear programming method to optimize the agricultural net benefit considering the scenarios of rainfall delay or reduction. Case studies were demonstrated in five farms and three municipalities of Western Bahia. The results show that in irrigated areas, the model optimizes the net economic benefit by the possibility of a continuous or double-cropping system, even in rainfall delay or reduction, where irrigation can supply the water demand of crops. In rainfed areas, it is noticeable that in rainfall delay or reduction scenarios, the model response is not to crop, due to the accentuated water deficit in crops, which may cause a significant yield reduction. It was found in a farm-level analysis, where the model response was not to crop, and farmers had a yield reduction of 61% in rainfed areas. This study opens the integration between the optimization methodologies and the hydroclimatic monitoring system with new insights into how this integration can guide water governance actions in regions where the water crisis is a reality.

1. Introduction

One of the main risks for society in the following years is the water crisis [1]. The projections of population increase [2], climate change [3], and the attribution to agriculture are mainly responsible for global water withdrawal [4], supporting this concern. Thus, global food production must increase in uncertain scenarios of water resources, and without local actions, the water crisis scenarios may become a reality rather than a risk.
Recent estimates by the Pacific Institute [5] show that global water conflicts have increased. These conflicts are reported mainly in developing countries [6,7,8,9] with a solid agricultural dependency on national economies. Gilbert [10] declared that to avert a global water crisis, we need to have better data related to water demand and availability and use these data for assertive management decisions.
The relationship between agriculture and the water crisis in Brazil is evident in a specific mesoregion, Western Bahia. Western Bahia is located in MATOPIBA (the acronym stands for the first two letters of the Brazilian states of Maranhão, Piauí, Tocantins, and Bahia), which is considered the last agricultural frontier in the country. The region is characterized by rapid changes in land cover and land use, predominantly over the natural vegetation of the Cerrado biome, through significant agricultural intensification [11]. The agrarian model practiced in Western Bahia has an impressive performance in the regional economy. The 2.4 million hectares destined for crop production concentrate 34.2% of the state agricultural gross domestic product [12]. This agricultural model is well known for the high concentration of irrigated areas. Despite covering only 5.81% of the cropland, they are responsible for approximately 30% of its agricultural production value [13]. However, despite the economic benefits of irrigated areas, these are the leading cause of water stress and conflicts in Western Bahia [14]. On the other hand, most of the croplands in the region are rainfed, with the adoption of intensive double-cropping systems in some areas where the rainy season is long enough. However, the sizeable interannual variability on the onset of the rainy season and an observed shortening of its duration significantly increases the risk of such intensive practices [15].
To avoid water conflicts in Western Bahia, Pousa et al. [14] and Mantovani et al. [16] suggested a hydroclimatic monitoring system capable of predicting water availability in the months of water scarcity and providing information about the surface and groundwater availability, land use, hydrological features, policy, and governance. A prototype of this system is available on the OBahia—Territorial and Water Intelligence Platform for Western Bahia (http://obahia.dea.ufv.br, accessed on 7 July 2020) [17]. For the information available on the OBahia platform to become a guide for assertive management decisions, we need to integrate this system with data analysis tools. For example, the “rainy season onset forecast” tool provides information that can be useful for planning the crop season, helping farmers to decide the best sowing dates for each crop, for resource allocation during crop season, and predicting harvest losses. However, without an effective tool for analyzing these data, this type of forecast would only be another piece of information whose usefulness would not be used by the system’s users.
Several studies have been developed to optimize agricultural output. Linear programming (LP) has superb usability in these studies, with the problem formulation adapted to different objectives such as (1) maximizing net economic benefits [18]; (2) reducing crop production costs [19,20]; (3) indicating resources allocation [21,22,23,24,25,26,27]; and (4) improving economic and environmental indicators [28,29]. The LP method was also demonstrated in different management scales, with multiple objectives, under parameter uncertainty, small data series, and various other conditions [30,31,32,33,34]. Despite the significant contribution of LP models, the integration of models with hydroclimatic forecast models is non-existent. This combination can favor agricultural water management in a region, as the forecast results can be analyzed automatically, allowing for in-depth analysis and decision making at different management levels.
Considering the hydroclimatic forecast provided for Western Bahia, it becomes necessary to study how the LP methods can be applied to improve the water management in the region, guiding farmers to better crop decisions and regional water authorities to effective management strategies. This study proposes an LP model for operational optimization in Western Bahia that is fully integrable with the hydroclimatic monitoring system information. The model optimizes agricultural output considering the in-farm restrictions and water availability that depend on the rainy season onset forecast. Specifically, this work formulates an LP model and solves it for applications at distinct water availability and use levels. It also indicates strategies to make the model application effective. The work structure is divided to describe the region, present how the dataset was acquired and processed, how the model was formulated, present the case studies, and discuss how the application can be helpful for farmers and regional water authorities.

2. Materials and Methods

This section describes the acquisition and processing of the dataset and the linear programming model (LPM) formulation. In Figure 1, the complete work flowchart also shows the steps adopted in the case studies, results, and discussion sections.

2.1. Region Description

Western Bahia (43.23° to 46.61° W, 10.10° to 15.26° S), represented in Figure 2, borders the Brazilian states of Minas Gerais, Goiás, Tocantins e Piauí. The region’s representative climate is tropical with a dry season, which favors the development of irrigated agriculture using the center pivot system combined with the flat topography. From 2000 to 2020, the planted area of soybean, maize, cotton, and bean represented 94% ± 2% of the total cultivated area in Western Bahia [35]. According to Pimenta et al. [36], in 2020, the total agricultural area was equal to 5.14 × 106 ha, and irrigated cropland was equal to 2.18 × 105 ha. The water grant in this region follows Federal Law n° 9.433/97 and State Law n° 11.612/09, conceding the use right of water resources by limited time to the applicant. This process in Western Bahia is coordinated by INEMA (Environment and Water Resources Institute of Bahia), which considers, for superficial withdrawal of 80% of Q90, the minimum amount of water flow in a river that is present 90% of the time as the maximum of water grant for each river basin.

2.2. Dataset Acquisition and Processing

The data acquisition process depends on the level of water governance considered in the case studies section. At the farm governance level, through interviews and local consultation, obtained data are related to irrigated and rainfed area availability, water grant, weather data, unit production cost, technical coefficients (unit need for machinery and manual labor), and adopted crops. At the municipal governance level, weather data were obtained from “Banco de Dados Meteorológicos do INMET” [37], selecting at least 30 years of daily data for each weather station. Unit production cost, technical coefficients, adopted crops, and revenue value were obtained from “Série Histórica de Dados” [38]. Irrigated and rainfed area availability and water grant were obtained from the OBahia platform [17]. Each crop’s unit revenue value was obtained from “Preços agrícolas, da sociobio e da pesca” [39]. Weather data refer to both cases’ maximum and minimum temperature, insolation, average relative humidity, and average wind speed.
At the farm level, the LPM was applied considering the original weather condition and in four scenarios of rainfall delay (D1, D2, D3, and D4, representing 1, 2, 3, and 4 fortnights in the onset of the rainy season, respectively) and rainfall reduction (R5, R10, R15, R20 representing 5%, 10%, 15%, and 20% of rainfall reduction, respectively). These scenarios were adopted to represent the possible responses of the hydroclimatic forecast. At the municipal level, the weather data were analyzed to determine the average and quartiles (Q25 and Q75 representing the 25% and 75% quantiles, respectively) for each fortnight to represent each municipality’s rainfall condition; average values of weather data represented the original climatic condition. In addition to analyzing rainfall delays and reductions, the LPM was also applied for quartile conditions. When available, the actual area occupation was compared with the results of the LPM application. In both cases, the resource availability and needs were converted to represent the fortnightly period, making the model application compatible with the hydroclimatic forecast available for the region on the OBahia platform [17]. All characteristics of the LPM application are fully described in the case studies section.

2.3. LPM Formulation

The LPM formulation is fully described by Vanderbei [40] and Luenberger and Ye [41]. An LPM model is applied to optimize a process considering a linear correspondence of variables. In other words, it is a method to optimize linear processes subject to actual decision variables and constraints. This subsection highlights the definition of decision variables, the objective function, constraints, and the solution method.

2.3.1. Decision Variables

This study’s decision variables (DV) are defined as the planted area of a crop, over irrigated or rainfed condition, with a specific crop cycle, produced in a particular fortnight. Mathematically, DV can be demonstrated, as X i j k l . X (in hectares) is the optimum value for the planted area; i is the crop type, which may be soybeans, maize, beans, wheat, or cotton; j is the planting system, which may be irrigated or rainfed; k is the crop cycle type, which may be short, average, or long; and l is the sowing period, limited by agroclimatic zoning made by the Brazilian National Supply Company (CONAB) [42]. The crop cycle type refers to the period between the sowing and harvest of the crop, defined by the genetics of the adopted variety.

2.3.2. Objective Function

The objective function (OF) of this optimization model is to maximize the net benefit (NB) (BRL ha−1) of agricultural activity, which has been mathematically expressed as:
Maximize   NB = R i j k l   ×   X i j k l     C i j k l   ×   X i j k l ,
where R (BRL ha−1) is the unitary revenue value, and C (BRL ha−1) is the unitary production cost.

2.3.3. Constraints and Solution Method

The constraints are based on the agricultural use of resources (Res) divided into land, work, machines, and capital. Land constraints refer to area (irrigated or rainfed) occupation; work constraints refer to labor and administrative activities; machines constraints refer to agricultural machinery’s availability as a planter or harvester; capital constraints mean the working capital of activity. In all cases, the resource demand was computed but was used as constraints when that information was substantial for the management level; all constraints’ setups are defined in the case study section. In all constraint conditions, the use of the resources must be equal to or less than availability. In general, the resource demand is described by Equation (2) as the product between the optimum DV and the unitary resource amount needed.
Res   demand = X i j k l   ×   Y ¯ i j k l ,
where Ȳ (dimensionless) is the adopted technical coefficient for each resource, representing the unitary amount needed.
Crop water demand (CWD) was the only constraint parameter calculated differently than described by Equation (2). CWD is calculated in steps, first, the reference evapotranspiration ( E T o , mm day−1) and the crop evapotranspiration ( E T c , mm day−1) accumulated for the period. E T c is the product between E T o and the crop coefficient ( K c , dimensionless), as indicated by Allen et al. [43] and described in Equation (3).
E T c = E T o   ×   K c i j k l
In Cerrado’s condition (Cerrado’s region includes Western Bahia), the reductions of 20% in E T c for cotton and soybean, 12% for maize, and 13% for beans do not cause a significant reduction in crop yields [44]. This work defined CWD as the product between the related E T c and a potential crop evapotranspiration reduction that does not cause a significant yield reduction ( K r , adimensional). In turn, the adopted K r for each crop was 0.80 for cotton and soybean, 0.88 for maize, 0.87 for beans, and 0.85 for wheat. Therefore, CWD for each DV is mathematically described by Equation (4).
C W D X i j k l = X i j k l   ×   E T c   ×   K r
The Solver tool of Microsoft Excel was used to set all the LPM applications. The “LP Simplex” was the solution method applied on all occasions.

3. Case Studies

3.1. Farms and Municipalities Characterization

Due to data availability, five farms served as case studies at the farm governance level: (1) Sama, where the LPM was applied to optimize the Soybean crop season with different cycle durations and sowing dates in rainfed and irrigated areas; (2) Floryl, for multiple crops with multiple cycle duration, (3) DRB MOR (“Decisão Rio Branco e Morena Farm”) and (4) DECSD (“Decisão São Desidério Farm”) for multiple crops during six seasons and finally (5) Busato, for multiple crops in a season. Sama, DRB MOR, DECSD, and Busato farms are included in the regions of intense irrigation systems growth presented by Pousa et al. [14]. Three municipalities served as case studies at the municipal governance level: (1) Barreiras, with the largest irrigated area; (2) Correntina, a region with the largest rainfed agricultural area; and (3) Santa Rita de Cássia, a region with the lowest technology agriculture. In all cases, cotton, beans, soybeans, and maize were evaluated during six crop seasons. For the study sites, the resources’ availability, crops, and season are summarized in Table 1.
The minimum and maximum values of ETo and rainfall observed for each farm are presented in Table 2.

3.2. Crop Characterization

Table 3 shows the adopted characteristics for each crop. These data are usually adopted by companies that provide irrigation consulting services, such as Valley Scheduling® (called in Brazil Irriger Connect®). The adopted technical coefficients, as labor and machinery needs for each condition, and the adopted revenue values are demonstrated in Supplementary Materials S1.

3.3. Model’s Constraints Setup

For each application and LPM analysis, a different constraint setup was adopted. These constraints may be the irrigated (Ia) or rainfed area (Ra), CWD irrigated (CwdI) or rainfed areas (DwdR), labor (L), or machinery (M). This setup difference was necessary due to the solution tool’s limitations and can be thoroughly checked in Table 4.

4. Results

This section presents the results of applying the LPM for the Western Bahia farms and municipalities.

4.1. LPM Application at Farms

At Sama farm, LPM results differ from real decisions made by farmers. In irrigated areas, the LPM’s sowing recommendation was later compared to the actual conditions, while the harvest indication was similar to the real condition. In rainfed areas, without considering crop water demand as a constraint, only a tiny difference in the harvest date is noticeable in comparing the LPM and actual condition due to changes in the crop cycle. These results are demonstrated in Figure 3.
In irrigated areas, considering rainfall delay conditions, the LPM’s responses considered the best price for crop commercialization and distribution of labor and machinery. The LPM’s response considering rainfall reduction scenarios did not differ between them because irrigation can supply the crop water demand during the crop season. When considering the crop water demand in rainfed areas, the LPM’s response was not to plant. This recommendation occurred because the crop’s water stress, in simulated conditions, was higher than the one proposed by Santos et al. [44], which could cause yield losses. According to Sama farm managers, the crops in this harvest suffered an average yield reduction of 12% and up to 61% (considering the expected and actual) in rainfed areas. The LPM’s application results for irrigated areas are demonstrated in Figure 4. Results for rainfed areas were omitted due to the null LPM response.
At Floryl farm, the LPM optimizes (in all conditions of rainfall delay or reduction) the area occupied in irrigated areas recommending the sowing of soybeans (from October to February) followed by maize first crop (from February to June) and then by maize second crop (from June to October). Double-cropping systems are standard in Brazil, where almost 58% of maize is produced under these conditions [15], but the adoption of continuous monoculture (maize first and second crop) may not be recommended due to the increase in pathogen population in the area, which can cause yield reduction of the second crop [45]. In rainfed areas, the LPM’s response was to plant maize first crop on all occasions and adjust the sowing date in rainfall delay conditions. Sowing at the D4 scenario was not recommended due to the CONAB’s sowing calendar [42]. This result is also consistent with the assessment made by Abrahão and Costa [15], where rainfall delay is one of the threats to double-cropping systems. All LPM results of area occupation can be verified in Figure 5.
The results for the DRB MOR farm were similar to the Floryl farm. In irrigated areas, the LPM’s responses favor higher area occupation by sowing cotton and bean second crops during November 2017 and September 2018, followed by sowing bean first crop, cotton, and wheat during November 2018 and January 2020, in all scenarios of rainfall delay or reduction. In rainfed areas, the recommendation was cotton after September 2019 and December 2020, adjusting the sowing dates according to rain delay conditions. It is perceptible that from 2017 to 2019, there were no recommended crops in rainfed areas due to the accentuated water deficit. This response shows that crop yield in rainfed areas can be significantly reduced, especially in a drier year. At the DECSD farm, rainfed areas are not cultivated, and the LPM responses are similar to those obtained for the DRB MOR farm. LPM results for the DRB MOR and DECS farms are shown in Figure 6.
At the Busato farm, in irrigated areas, the LPM recommends a double-cropping system composed of cotton and maize second crops for all occasions of rainfall delay or reduction. In rainfed areas, only cotton is recommended, adjusting the sowing date in rainfall delay scenarios. Again, the double-cropping system is threatened by accentuated water deficit, as pointed out by Abrahão and Costa [15]. All results are shown in Figure 7.

4.2. LPM Application at Municipalities

Results of the LPM application at Barreiras show an excellent optimization of area occupation. The predominant recommendation is the double-cropping system, with beans first crop and cotton over irrigated and rainfed areas. Irrigated areas maintain the same results for all evaluated conditions. Over rainfall delay conditions, rainfed lands cannot be capable of a double-cropping system. In extremely dry year scenarios (Q25), the LPM recommendation is not to plant in rainfed areas, indicating a possibility of a decline in crop yield. Correntina and Santa Rita de Cássia results resembled the Barreiras results, with cotton and beans in a double-cropping system for irrigated areas. The model indicated an accentuated water deficit in rainfed areas that may cause a significant yield reduction, threatening the double-cropping systems.
Among those evaluated, Barreiras is the municipality with the most irrigated area to total cropland (12.78%). On the other hand, Correntina and Santa Rita de Cássia have smaller proportions (3.54% and 0.81%, respectively). Despite not being in regions of intense irrigation growth [14], Correntina has been the scene of conflicts over water resources [46]. Santa Rita de Cássia and other regions not evaluated in the study have a vast potential for irrigated agriculture development. Still, suitability and technical limitations must be considered (as energy availability, topography, climatic conditions, and logistical aspects), and beyond these, as well as the development’s environmental, economic, and social implications. The complete results of the LPM application to Barreiras, Correntina, and Santa Rita de Cássia are shown in Figure 8.

5. Discussion

This work presents an integrated approach between the hydroclimatic monitoring system and a linear programming model to guide actions of water governance in Western Bahia. Given the proposed model formulation, optimizing agricultural activity’s net economic benefit in rainfall delay and reduction scenarios was possible. This optimization was made by the recommendation of crops, cycle type, and sowing dates, considering the limitations of water, land, capital, and labor of each farm or municipality. In this way, it becomes possible to plan a crop season based on information on water availability in the region and avoid significant water stress in crops. This information is substantial to proactive water management decisions in Western Bahia.
The presented results in this work are consistent with those found in other studies. There was success in applying the model at different levels of water use, such as farms and municipalities, in the same way, demonstrated by several authors [21,22,23,24,25,26,27]. At the farm level, it becomes perceptible that the model can recommend or not crops, given the water availability scenarios. In the irrigated area, the model response is generally directed to select the crop with the maximum net benefit. In rainfed areas, the model recommendation was not to crop occasionally. It occurred due to the expected water deficit in a crop season or to the scenarios of rainfall delay when the beginning of the rainy season is out of the sowing calendar. This occurrence was noticed in the Sama farm simulation, where considering the constraint of water demand, the model response was not to crop in rainfed areas. These areas were cropped, but the yield reduction stood between 12% and 61%. The linear programming response at the municipal farm level reinforces the possibility of a continuous cropping system in irrigated areas and a double or single cropping system over rainfed areas. Analyzing the model responses in more depth, we should declare the limitations and potentials of the proposed integration for the study region.
The proposed LPM has limitations. First, being a uni-objective model, with this formulation, it is not possible to optimize other essential points for agricultural water management, such as increasing crop yields or efficient water allocation. Nonetheless, this limitation opens the possibility of developing and comparing other optimization methodologies, such as multi-objective linear programming, making integrating the systems even more helpful. Second, agricultural systems have a linear response to some variables, such as crop yield response to water stress and variations in general cost. The nonlinearity of the variables implies that we cannot guarantee that the forecasts of net benefit from the activity will be the same under actual conditions. We must also declare that the biweekly period considered in all case studies impacts the model’s assertiveness, given that, in reality, crop management decisions such as irrigation are taken daily. Third, the LPM application does not consider the spatial variation of climatic parameters, representing a limitation for governance in the highest and most extensive water levels as in the entire mesoregion. Furthermore, the Microsoft Excel® Solver tool, used to solve the linear programming model, has a severe limitation of 200 decision variables and 100 restrictions for each evaluation. Therefore, it was necessary to divide the evaluation into shorter periods to make the model’s size compatible with the tool.
As for the potential of this integration, we can cite the possibility of evaluating the area’s productive potential over rainfall delay and reduction scenarios occasioned by climate change. The effects of rainfall delay significantly affect the double-cropping system practiced over rainfed areas, as demonstrated by Abrahão and Costa [15]. With the advanced information on the hydroclimatic system, we could apply the proposed model to simulate crops and variety with different water stress tolerance over the expected climate scenarios. Another consideration of the integration is to simulate the economic impacts of increasing or decreasing the water grants and, consequently, the irrigated area over a region. Analyzing the model results, especially at the municipal level, it becomes possible to estimate the labor, machinery, and capital needs, guiding job creation, industry strengthening, and funding actions for farmers. Finally, being simple and requiring few input variables, the model can be easily understood and quickly disseminated to farmers, institutions, and decision makers that contribute to the water management process.

Strategy for LPM Effectiveness

The following actions are recommended to ensure that the model presented is effectively used as a water resources management tool in Western Bahia: (1) Develop and make available an online tool with the presented linear programming model, which can be made available on the OBahia platform [17] integrating water availability data and the rainy season onset forecast; (2) adjust the technical crop coefficients (need of resources, yield, and K c for example) for regional conditions and make them available by default on the platform; (3) promote training for farmers and other decisions makers aiming at the dissemination and correct use of the online tool, demonstrating its potential for use and the possibilities of analysis.

6. Conclusions

Concerns about the water crisis and its environmental, social, and economic impacts on agricultural activity have increased in recent years, especially in regions where irrigated agriculture has extensive development, such as Western Bahia, Brazil. In this region, a hydroclimatic monitoring system was implemented, which, through data availability, aims to avoid water conflicts over water use, especially in periods of greater water scarcity. Considering this monitoring system, it became necessary to develop integrated tools that combine the water availability data and optimization tools, generating a valuable guide for water management actions. This study proposes a linear programming model that can be used for this purpose. The model maximizes the net benefit of agricultural activity considering land, labor, capital, machinery, and water restrictions. Its application was demonstrated at farm and municipality levels of water use and with different simulated rainfall delay and reduction conditions. The results show that the model response can guide decision making, especially when the adoption of double or continuous cropping systems is expected. For farmers, the model can provide information about which crop would be the best option for sowing and when to do so, considering the rainfall delay or reduction scenario given by the hydroclimatic forecast. The application can also estimate the economic impacts of new water grants, which is helpful for regional water authorities. Furthermore, the presented integration between a linear programming model and a hydroclimatic monitoring system is a scientific novelty. Integrating a linear programming model and a hydroclimatic monitoring system can better guide water management decisions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w14223625/s1, Spreadsheets S1: adopted coefficients.

Author Contributions

Conceptualization, project administration and funding acquisition, M.H.C.; field data collection, I.B., E.C.M. and A.G.d.S.J.; formal analysis, investigation and data curation, I.B.; writing—original draft preparation, I.B.; writing—review and editing, I.B., E.C.M., M.H.C. and A.G.d.S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Coordination for the Improvement of Higher Education Personnel (CAPES), Finance Code: 001, by PRODEAGRO, grant numbers 011/2016 and 045/2019, and by the National Council for Scientific and Technological Development—Brazil (CNPq), process 132424/2019-3 and 441.210/2017-1.

Data Availability Statement

Example of LPM application and instructions are available on the OBahia Platform (http://obahia.dea.ufv.br/#/ruralprofit, accessed on 7 July 2020).

Acknowledgments

Gustavo B. Braga has read and commented on previous versions of this manuscript. We are thankful for his time and expertise. We also thank the farmers from Western Bahia and AIBA for providing data and support during data acquirement.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study.

References

  1. World Economic Forum. The Global Risks Report 2020; World Economic Forum: Geneva, Swizetzerland, 2020; ISBN 9782940631247. [Google Scholar]
  2. Department of Economic and Social Affairs of United Nations. World Population Prospects 2019; United Nations, Ed.; United Nations: New York, NY, USA, 2019; ISBN 9789211483161. [Google Scholar]
  3. IPCC. Summary for Policymakers. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  4. WWAP. The United Nations World Water Development Report 2016: Water and Jobs; UNESCO: Paris, France, 2016; ISBN 9789231001550. [Google Scholar]
  5. Pacific Institute Water Conflict Chronology. Available online: https://www.worldwater.org/conflict/map/ (accessed on 10 June 2022).
  6. Rawlins, J. Political economy of water reallocation in South Africa: Insights from the Western Cape water crisis. Water Secur. 2019, 6, 100029. [Google Scholar] [CrossRef]
  7. Muñoz, A.A.; Klock-Barría, K.; Alvarez-Garreton, C.; Aguilera-Betti, I.; González-Reyes, Á.; Lastra, J.A.; Chávez, R.O.; Barría, P.; Christie, D.; Rojas-Badilla, M.; et al. Water crisis in petorca basin, Chile: The combined effects of a mega-drought and water management. Water 2020, 12, 648. [Google Scholar] [CrossRef]
  8. Varisco, D. Pumping Yemen Dry: A History of Yemen’s Water Crisis. Hum. Ecol. 2019, 47, 317–329. [Google Scholar] [CrossRef]
  9. Shan, V.; Singh, S.K.; Haritash, A.K. Water Crisis in the Asian Countries: Status and Future Trends. In Resilience, Response, and Risk in Water Systems; Springer: Singapore, 2020; pp. 173–194. [Google Scholar] [CrossRef]
  10. Gilbert, N. How to avert a global water crisis. Nature 2010. [Google Scholar] [CrossRef]
  11. Dias, L.C.P.; Pimenta, F.M.; Santos, A.B.; Costa, M.H.; Ladle, R.J. Patterns of land use, extensification, and intensification of Brazilian agriculture. Glob. Chang. Biol. 2016, 22, 2887–2903. [Google Scholar] [CrossRef]
  12. AIBA. Dados e Pesquisa: Região Oeste. Available online: https://aiba.org.br/regiao-oeste/ (accessed on 25 January 2021).
  13. AIBA. Atlas do Oeste da Bahia; AIBA: Barreiras, Brazil, 2019. [Google Scholar]
  14. Pousa, R.; Costa, M.H.; Pimenta, F.M.; Fontes, V.C.; Castro, M. Climate change and intense irrigation growth in Western Bahia, Brazil: The urgent need for hydroclimatic monitoring. Water 2019, 11, 933. [Google Scholar] [CrossRef]
  15. Abrahão, G.M.; Costa, M.H. Evolution of rain and photoperiod limitations on the soybean growing season in Brazil: The rise (and possible fall) of double-cropping systems. Agric. For. Meteorol. 2018, 256, 32–45. [Google Scholar] [CrossRef]
  16. Mantovani, E.C.; da Silva Júnior, A.G.; Costa, M.H.; Marques, E.A.G.; da Silva Júnior, G.C.; Pruski, F.F. Relatório Técnico Final: Estudo do Potencial Hídrico da Região Oeste da Bahia: Quantificação e Monitoramento da Disponibilidade dos Recursos do Aquífero Urucuia e Superficiais nas Bacias dos rios Grande, Corrente e Carinhanha; UFV: Viçosa, Brazil, 2019. [Google Scholar]
  17. Pimenta, F.M.; Silva, M.L.; Castro, M.; Fontes, V.C.; Costa, M.H. OBahia—Western Bahia Map Server. 2020. Available online: https://www.escavador.com/sobre/7173300/fernando-martins-pimenta (accessed on 12 May 2022).
  18. Zhang, C.; Guo, P. An inexact CVaR two-stage mixed-integer linear programming approach for agricultural water management under uncertainty considering ecological water requirement. Ecol. Indic. 2018, 92, 342–353. [Google Scholar] [CrossRef]
  19. Ma, T.; Wang, J.; Liu, Y.; Sun, H.; Gui, D.; Xue, J. A mixed integer linear programming method for optimizing layout of irrigated pumping well in Oasis. Water 2019, 11, 1185. [Google Scholar] [CrossRef]
  20. Sheibani, H.; Shourian, M. Determining optimum reliability for supplying agricultural demand downstream of a reservoir using an explicit method with an economic objective function. Water Resour. Econ. 2019, 26, 100131. [Google Scholar] [CrossRef]
  21. Aljanabi, A.A.; Mays, L.W.; Fox, P. Optimization model for agricultural reclaimed water allocation using mixed-integer nonlinear programming. Water 2018, 10, 1291. [Google Scholar] [CrossRef]
  22. Chen, S.; Xu, J.; Li, Q.; Tan, X.; Nong, X. A copula-based interval-bistochastic programming method for regional water allocation under uncertainty. Agric. Water Manag. 2019, 217, 154–164. [Google Scholar] [CrossRef]
  23. Freire-González, J.; Decker, C.A.; Hall, J.W. A linear programming approach towater allocation during a drought. Water 2018, 10, 363. [Google Scholar] [CrossRef]
  24. Fu, J.; Zhong, P.A.; Zhu, F.; Chen, J.; Wu, Y.N.; Xu, B. Water resources allocation in transboundary river based on asymmetric Nash-Harsanyi Leader-Follower game model. Water 2018, 10, 270. [Google Scholar] [CrossRef]
  25. Fu, Q.; Li, L.; Li, M.; Li, T.; Liu, D.; Cui, S. A simulation-based linear fractional programming model for adaptable water allocation planning in the main stream of the Songhua River basin, China. Water 2018, 10, 627. [Google Scholar] [CrossRef]
  26. Jin, L.; Fu, H.; Kim, Y.; Long, J.; Huang, G. A Bi-objective pseudo-interval T2 linear programming approach and its application towater resources management under uncertainty. Water 2018, 10, 1545. [Google Scholar] [CrossRef]
  27. Veintimilla-Reyes, J.; De Meyer, A.; Cattrysse, D.; Tacuri, E.; Vanegas, P.; Cisneros, F.; Van Orshoven, J. MILP for optimizing water allocation and reservoir location: A case study for the Machángara river basin, Ecuador. Water 2019, 11, 1011. [Google Scholar] [CrossRef]
  28. Bekchanov, M.; Ringler, C.; Bhaduri, A.; Jeuland, M. Optimizing irrigation efficiency improvements in the Aral Sea Basin. Water Resour. Econ. 2016, 13, 30–45. [Google Scholar] [CrossRef]
  29. Zhao, J.; Li, M.; Guo, P.; Zhang, C.; Tan, Q. Agricultural water productivity oriented water resources allocation based on the coordination of multiple factors. Water 2017, 9, 490. [Google Scholar] [CrossRef]
  30. Suo, M.; Wu, P.; Zhou, B. An integrated method for interval multi-objective planning of a water resource system in the eastern part of Handan. Water 2017, 9, 528. [Google Scholar] [CrossRef]
  31. Li, M.; Fu, Q.; Singh, V.P.; Liu, D. An interval multi-objective programming model for irrigation water allocation under uncertainty. Agric. Water Manag. 2018, 196, 24–36. [Google Scholar] [CrossRef]
  32. Li, M.; Fu, Q.; Singh, V.P.; Liu, D.; Li, T.; Zhou, Y. Managing agricultural water and land resources with tradeoff between economic, environmental, and social considerations: A multi-objective non-linear optimization model under uncertainty. Agric. Syst. 2020, 178, 102685. [Google Scholar] [CrossRef]
  33. Zhang, F.; Guo, S.; Liu, X.; Wang, Y.; Engel, B.A.; Guo, P. Towards sustainable water management in an arid agricultural region: A multi-level multi-objective stochastic approach. Agric. Syst. 2020, 182, 102848. [Google Scholar] [CrossRef]
  34. Niu, G.; Zheng, Y.; Han, F.; Qin, H. The nexus of water, ecosystems and agriculture in arid areas: A multiobjective optimization study on system efficiencies. Agric. Water Manag. 2019, 223, 105697. [Google Scholar] [CrossRef]
  35. IBGE. Produção Agrícola Municipal. Available online: https://sidra.ibge.gov.br/tabela/5457 (accessed on 9 June 2022).
  36. Pimenta, F.M.; Speroto, A.T.; Costa, M.H.; Dionizio, E.A. Historical changes in land use and suitability for future agriculture expansion in Western Bahia, Brazil. Remote Sens. 2021, 13, 1088. [Google Scholar] [CrossRef]
  37. MAPA INMET. Available online: https://portal.inmet.gov.br/ (accessed on 19 August 2020).
  38. CONAB Série Histórica das Safras. Available online: https://www.conab.gov.br/info-agro/safras/serie-historica-das-safras (accessed on 25 January 2021).
  39. CONAB Preços Agrícolas, da Sociobio e da Pesca. Available online: http://sisdep.conab.gov.br/precosiagroweb/ (accessed on 20 January 2021).
  40. Vanderbei, R.J. Linear Programming. In International Series in Operations Research & Management Science, 4th ed.; Hillier, F.S., Ed.; Springer: Berlin/Heidelberg, Germany, 2014; Volume 196, ISBN 9781461476290. [Google Scholar]
  41. Luenberger, D.G.; Ye, Y. Linear and Nonlinear Programming, 4th ed.; Price, C.C., Zhu, J., Hillier, F.S., Eds.; Springer: Cham, Switzerland, 2016; Volume 228, ISBN 978-3-319-18841-6. [Google Scholar]
  42. CONAB. Calendário de Plantio e Colheita de Grãos no Brasil 2019; CONAB: Brasilia, Brazil, 2019. [Google Scholar]
  43. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration (Guidelines for Computing Crop Water Requirements); FAO: Rome, Italy, 1998; Volume 300, p. D05109. [Google Scholar]
  44. Santos, I.S.; Mantovani, E.C.; Venancio, L.P.; da Cunha, F.F.; Aleman, C.C. Controlled water stress in agricultural crops in Brazilian Cerrado. Biosci. J. 2020, 36, 886–895. [Google Scholar] [CrossRef]
  45. Bullock, D.G. Crop rotation. CRC Crit. Rev. Plant Sci. 1992, 11, 309–326. [Google Scholar] [CrossRef]
  46. G1 Grupo Invade Fazendas e Incendeia Galpão em Protesto no Oeste da Bahia. 2017. Available online: https://g1.globo.com/bahia/noticia/grupo-invade-fazendas-e-incendeia-galpao-em-protesto-no-oeste-da-bahia.ghtml (accessed on 7 July 2020).
Figure 1. Work flowchart.
Figure 1. Work flowchart.
Water 14 03625 g001
Figure 2. General view of Western Bahia region. (a) This region is defined by the union of basins of the Grande, Corrente, and Carinhanha Rivers. Part of the Carinhanha River basin is over the territory of Minas Gerais State. (b) Location of farms and municipalities where the linear programming model was applied in case studies (characterized in Section 3.1). Red dots represent the farms’ locations. Green areas correspond to municipal areas. Red lines pattern fill represent areas more susceptible to the water crisis in Western Bahia, evaluated by Pousa et al. [14].
Figure 2. General view of Western Bahia region. (a) This region is defined by the union of basins of the Grande, Corrente, and Carinhanha Rivers. Part of the Carinhanha River basin is over the territory of Minas Gerais State. (b) Location of farms and municipalities where the linear programming model was applied in case studies (characterized in Section 3.1). Red dots represent the farms’ locations. Green areas correspond to municipal areas. Red lines pattern fill represent areas more susceptible to the water crisis in Western Bahia, evaluated by Pousa et al. [14].
Water 14 03625 g002
Figure 3. Planted area, in hectares, recommended by LPM application for irrigated and rainfed areas of Sama farm.
Figure 3. Planted area, in hectares, recommended by LPM application for irrigated and rainfed areas of Sama farm.
Water 14 03625 g003
Figure 4. Planted area, in hectares, recommended by LPM application for irrigated areas of Sama farm considering rainfall delay and reduction scenarios.
Figure 4. Planted area, in hectares, recommended by LPM application for irrigated areas of Sama farm considering rainfall delay and reduction scenarios.
Water 14 03625 g004
Figure 5. Planted area, in hectares, recommended by LPM application at Floryl farm. (a) Shows the results for irrigated areas. (b) Shows the results for rainfed areas.
Figure 5. Planted area, in hectares, recommended by LPM application at Floryl farm. (a) Shows the results for irrigated areas. (b) Shows the results for rainfed areas.
Water 14 03625 g005
Figure 6. Planted area, in hectares, recommended by LPM application at DRB MOR and DECSD farms. (a) Shows the results for irrigated areas of DRBMOR farm. (b) Shows the results for rainfed areas of the DRB MOR farm. (c) Shows the results for irrigated areas of DECSD farm.
Figure 6. Planted area, in hectares, recommended by LPM application at DRB MOR and DECSD farms. (a) Shows the results for irrigated areas of DRBMOR farm. (b) Shows the results for rainfed areas of the DRB MOR farm. (c) Shows the results for irrigated areas of DECSD farm.
Water 14 03625 g006
Figure 7. Planted area, in hectares, recommended by LPM application at Busato farm. (a) Shows the results for irrigated areas. (b) Shows the results for rainfed areas.
Figure 7. Planted area, in hectares, recommended by LPM application at Busato farm. (a) Shows the results for irrigated areas. (b) Shows the results for rainfed areas.
Water 14 03625 g007
Figure 8. Planted area, in hectares, recommended by LPM application for irrigated and rainfed conditions of Barreiras, Correntina, and Santa Rita de Cássia. (a) Shows the results for irrigated areas of Barreiras. (b) Shows the results for rainfed areas of Barreiras. (c) Shows the results for irrigated areas of Correntina. (d) Shows the results for rainfed areas of Correntina. (e) Shows the results for irrigated areas of Santa Rita de Cássia. (f) Shows the results for rainfed areas of Santa Rita de Cássia.
Figure 8. Planted area, in hectares, recommended by LPM application for irrigated and rainfed conditions of Barreiras, Correntina, and Santa Rita de Cássia. (a) Shows the results for irrigated areas of Barreiras. (b) Shows the results for rainfed areas of Barreiras. (c) Shows the results for irrigated areas of Correntina. (d) Shows the results for rainfed areas of Correntina. (e) Shows the results for irrigated areas of Santa Rita de Cássia. (f) Shows the results for rainfed areas of Santa Rita de Cássia.
Water 14 03625 g008
Table 1. Description of resource availability, crops, and crop season adopted in LPM’s analysis in municipalities and farms.
Table 1. Description of resource availability, crops, and crop season adopted in LPM’s analysis in municipalities and farms.
Farm or MunicipalityIrrigated AreaRainfed AreaWater GrantEvaluated CropSeasonLabor
Availability
Tractor
Availability
Truck
Availability
Sprayer
Availability
Harvester
Availability
(Hectares)(Hectares)(m³ day−1)(h day−1)(h day−1)(h day−1)(h day−1)(h day−1)
Sama f1221.10878.5039,297.50Soybean2018/1916016163216
Floryl f950.00120.0068,437.00Maize 1−2; Soybean2019/20Not applicable
DRB MOR f4267.00646.00101,760.00Cotton; Bean 1−2; Maize 1−2; Soybean; Wheat2017 aw; 2017 ss; 2018 aw; 2018 ss; 2019 aw; 2019 ss
DECDS f3299.000.00156,552.00
Busato f4246.80676.0023,240.00Cotton; Maize 1−2; Soybean2019/20
Barreiras m42,760.00291,948.004,586,256.52Cotton; Bean 1−2; Maize 1−2; Soybean2018 ss; 2019 aw; 2019 ss; 2020 aw; 2020 ss; 2021 aw
Correntina m15,008.00408,022.002,444,307.23
Santa Rita de Cássia m60.007293.0044,000.00
f Represents the farms, m represents the municipality, aw means the autumn/winter season, ss illustrates the spring/summer season, and 1–2 represents the crop cultivated in 1st or 2nd season in a double-cropping system. At Sama farm, the labor, tractor, truck, sprayer, and harvester availability, in h day−1 represents the total resource availability in one day, given the product between the amount of machinery or employees and the work journey in a day, equal to eight hours.
Table 2. Minimum and maximum rainfall values observed for each farm or municipality.
Table 2. Minimum and maximum rainfall values observed for each farm or municipality.
Farm or MunicipalityMinimum EToMaximum EToMinimum RainfallMaximum Rainfall
(mm fortnight−1)(mm fortnight−1)(mm fortnight−1)(mm fortnight−1)
Sama f45.296.90.0172.8
Floryl f42.786.70.0142.8
DRB MOR f46.2107.50.0142.8
DECDS f43.189.90.0232.2
Busato f41.194.70.0168.8
Barreiras m54.788.60.2103.7
Correntina m54.584.80.2108.7
Santa Rita de Cássia m59.091.60.0101.5
f Represents the farms, m represents the municipality.
Table 3. Crop, cycle type, duration Kc, Kr, and ranging sowing time.
Table 3. Crop, cycle type, duration Kc, Kr, and ranging sowing time.
CropCycle TypeCycle DurationInitial KcAverage KcFinal KcKr [44]Range Sowing Time
(Fortnight, or 15 Days)(Fortnight–Month) [42]
Soybeanshort80.600.700.800.801–10 to 2–02
average9
long10
Maize 1st seasonaverage90.651.000.600.881–10 to 2–02
long120.600.50
Maize 2nd Seasonaverage90.651.000.600.881–05 to 2–06
long120.600.50
Cottonlong140.500.900.380.801–11 to 2–02
Bean 1st seasonaverage70.701.200.600.871–10 to 2–02
Bean 2nd seasonaverage70.701.200.600.871–04 to 2–06
Wheataverage80.701.200.400.851–08 to 2–09
Table 4. Constraint setup for each LPM application and analysis.
Table 4. Constraint setup for each LPM application and analysis.
AnalysisFarms Municipalities
SamaFlorylDRB MORDECSDBusatoBarreirasCorrentinaSanta Rita de Cássia
OriginalIa; Ra; CwdI----Ia; Ra; CwdI; CwdRIa; Ra; CwdI; CwdRIa; Ra; CwdI; CwdR
Irrigated originalIa; CwdI; L; MIa; CwdIIa; CwdIIa; CwdIIa; CwdIIa; CwdIIa; CwdIIa; CwdI
Irrigated rain delay
Irrigated rain reduction
Rainfed originalRa; CwdR; L; MRa; CwdRRa; CwdR-Ra; CwdRRa; CwdRRa; CwdRRa; CwdR
Rainfed rain delay
Rainfed rain reduction
Ia means irrigated area demand, Ra means rainfed area demand, CwdI means crop water demand in irrigated areas, CwdR means crop water demand in rainfed areas, L means labor demand, and M means machinery demand.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boninsenha, I.; Mantovani, E.C.; Costa, M.H.; da Silva Júnior, A.G. A Linear Programming Model for Operational Optimization of Agricultural Activity Considering a Hydroclimatic Forecast—Case Studies for Western Bahia, Brazil. Water 2022, 14, 3625. https://doi.org/10.3390/w14223625

AMA Style

Boninsenha I, Mantovani EC, Costa MH, da Silva Júnior AG. A Linear Programming Model for Operational Optimization of Agricultural Activity Considering a Hydroclimatic Forecast—Case Studies for Western Bahia, Brazil. Water. 2022; 14(22):3625. https://doi.org/10.3390/w14223625

Chicago/Turabian Style

Boninsenha, Igor, Everardo Chartuni Mantovani, Marcos Heil Costa, and Aziz Galvão da Silva Júnior. 2022. "A Linear Programming Model for Operational Optimization of Agricultural Activity Considering a Hydroclimatic Forecast—Case Studies for Western Bahia, Brazil" Water 14, no. 22: 3625. https://doi.org/10.3390/w14223625

APA Style

Boninsenha, I., Mantovani, E. C., Costa, M. H., & da Silva Júnior, A. G. (2022). A Linear Programming Model for Operational Optimization of Agricultural Activity Considering a Hydroclimatic Forecast—Case Studies for Western Bahia, Brazil. Water, 14(22), 3625. https://doi.org/10.3390/w14223625

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop