1. Introduction
It is self-evident to state that agricultural land is essential to life and a valuable resource for most, if not all, countries. Still, its importance seems to only be noticed when productivity declines and our food security is at stake. Soil has been poorly protected and overexploited for decades, which has resulted in its deterioration worldwide (see
Figure 1).
Modern agriculture and farming are important human-induced factors in this degradation. Whereas traditional practices are based on organic fertilization and diversified crop production, modern agriculture relies on intensive single-crop production and chemical fertilizers, which together progressively degrade soil quality and reduce its productivity by changing its physical, chemical, and biological composition. Indeed, fertilization and low-quality irrigation water alters the soil’s chemical makeup. Further, plowing, tillage, removal of vegetative cover, and overgrazing make soil more vulnerable to wind and water erosion, and intensive and specialized cultivation exhausts some minerals and water from the soil and damages its microfauna (Blanco-Canqui and Lal [
1]).
The shift from traditional practices to highly specialized agricultural systems was, in part, a result of population growth, which increased demand. This shift increased soil productivity in the short term, but generated serious ecological drawbacks in the long term (degradation of soil quality, pollution of water and air, loss of biodiversity, erosion, etc.) and even reversed the trend of agricultural productivity. Modern agricultural systems have since fallen in a vicious cycle, where increasing chemical fertilizers are used to compensate for the loss of productivity, causing more damage to the soil and the environment (pollution of water and air) (Trautmann et al. [
2], Conway [
3]).
Implementing eco-responsible agricultural practices, based on crop rotations and mixed crop-livestock associations that are less harmful for the soil, is one way to achieve resource (soil) sustainability (Altieri [
4], Kremen et al. [
5]). At a macro level, this transition is a long-term process that must take populations’ food needs and farming profitability into account. Such a transition is particularly needed in island regions, given the importance of agriculture to their economies (Angeon et al. [
6]). A survey of farmers in the French West Indies revealed that the population is aware of the problem that farmers are making soil quality central to their concerns, and they are willing to make efforts and even to sacrifice some of their financial benefits to restore the quality of their land (Angeon et al. [
7]).
In this paper, we adopt a micro point of view and consider the problem of the long-term sustainable operation of a single farm, from both the physical (soil quality) and economic (farmer’s revenues) perspectives. Our empirical terrain is Guadeloupe, an island that is part of the French West Indies. In this region of the world, the restoration of soil quality depends (as everywhere else) on the soil’s inherent properties, the type of agricultural practice adopted by the farmer and on recurrent climatic events such as hurricanes.
Our research questions are as follows:
Given some economic and soil quality constraints, and taking into account the possible occurrence of climatic events, what are the viable crop rotations to grow over a predefined planning horizon?
How sensitive are these solutions with respect to parameter values?
Which farming systems or sectors are the most vulnerable to climatic uncertainties?
Given our problem statement, we adopt viability theory (Aubin [
8]) as our methodological framework. There is now a long tradition of applying this theory to address problems that are related to the exploitation of renewable resources, e.g., fisheries and forests (see Oubraham and Zaccour [
9] for a literature review). A few have applied viability theory to farming; see, e.g., Tichit et al. [
10], Sabatier et al. [
11], Sabatier et al. [
12], Mouysset et al. [
13], Tichit et al. [
14], Baumgärtner and Quaas [
15], Martin et al. [
16], Sabatier et al. [
17], and Durand et al. [
18]; most of which have focused on herd and grassland management problems.
To the best of our knowledge, Durand et al. [
18] is the only study that looked at a soil-quality management problem. When considering a single plot of agricultural land, the authors proposed a viability theory-based deterministic model, and looked for agricultural strategies and crop rotations that would restore soil quality to an acceptable level while preserving the farm’s economic profitability. In this paper, we extend their model to a stochastic setting to account for the impact of uncertainty from major climatic events, such as hurricanes, on the system’s evolution. During one season of the year, it is almost a certainty that the West Indies region will be affected by climatic events. The uncertainty relates to their force and the damage they will cause. Given this state of affairs, we believe our model adds a level of realism to what was done in Durand et al. [
18].
The rest of the paper is organized, as follows: in
Section 2, we extend the model in Durand et al. [
18] by introducing climatic uncertainty. In
Section 3, we define so-called emergency control, that is, what needs to be done after a climatic event.
Section 4 is dedicated to the formal definition of all the elements related to climatic events. The solution method is introduced in
Section 5, and the empirical application is presented in
Section 6. Finally, we briefly conclude in
Section 7.
2. A Stochastic Bio-Economic Model
In this section, we extend the model in Durand et al. [
18] to take into account the uncertainty related to major climatic events, e.g., cyclones and hurricanes. As in Durand et al. [
18], we consider a one-hectare parcel farm managed over time by a single agent. It is a commercial farm whose production is intended for sale. The planning horizon is
T and the current time is denoted by
. The state of the farm is described by two variables: (i) the cash flow
, i.e., a continuous variable characterizing the financial situation; and, (ii) the soil quality of the parcel, which is an agronomic metric. Evaluating the quality of a soil is a complex operation, as it involves a series of physical, chemical, and biological characteristics. Here, we adopt the General Indicator of Soil Quality (GISQ), which has values in the interval
and is based on 54 variables measuring these characteristics (Velásquez et al. [
19] and Camacho et al. [
20]). Let
be the GISQ value of the parcel, with
, which is,
can take a finite number of values.
The evolution of the soil quality depends on its current value
and on the following two control variables: (i) the crop
grown on the parcel, with
; and, (ii) the agricultural practice
To illustrate, in the case studies to follow, the set of crops includes either all or some of the following: plantain, export banana, sugar cane, yam (yellow), yam (Grosse Caille), tomato, eggplant, lettuce, carrot, green bean, cabbage, cassava, melon, cucumber, and turban squash. The set of agricultural practices only includes two elements: (i) conventional practice, which is based on modern practices/tools and chemical fertilization; and, (ii) agroecological practice, i.e., traditional practice based on organic fertilization of the soil or other practices for improved management of agroecosystems. Wezel et al. [
21] Let the evolution of the soil quality be described by the following function:
Denote by
the duration of the whole agricultural cycle of a crop
using agricultural practice
. Subsequently, to each pair
, we can associate production cycles
, defined as follows:
Moreover, to each crop and practice , we associate its first harvest time and the time duration between two successive harvests . The duration of the first production cycle is , and the following cycles last . If a crop has a single harvest, then and . Additionally, we suppose that the produce is sold right away after it is harvested.
The revenue from each harvest depends on the crop and the agricultural practice, as well as on the soil quality at the beginning of the cycle. Denote by
ℓ the revenue function defined by
Denote by
the current month. This discrete variable is needed to deal with the seasonality of some crops, which can be planted or sown only at specific times of the year. Let
be the set of all crops that can be planted or sown during season
s. In the deterministic model presented in Durand et al. [
18], the only event that must be accounted for over time is the beginning/end of a production cycle. If we use
n to refer to an event, and
to its timing, then the next event will happen at
The above equation assumes that once an agricultural cycle starts, it will reach its end. Here, we suppose that a production cycle may be randomly interrupted by a climatic event, which can occur during a specific period of the year, i.e., hurricane season. Consequently, (
1) needs to be modified to account for such an event. We introduce the following additional variables:
: the remaining time before the next end/start of a production cycle,
: the remaining time before the next potential hurricane strikes, and
: the remaining time before the next change in the system’s state,
Denote by the vector of events.
Denote by the intensity of a climatic event at time t. We suppose that takes its values from a finite and discrete set where 0 represents the absence of a climatic event and is the highest possible intensity. In the case studies, we assume that where 0 corresponds to no hurricane, 1 to a minor hurricane, which is, a hurricane of intensity 1 or 2 on the Saffir-Simpson hurricane wind scale, and 2 to a major hurricane, that is, a hurricane of intensity 3, 4, or 5 on the Saffir-Simpson scale. In line with what is typically observed in the West Indies, which is the area targeted by our case studies, we assume that a climatic event can happen at most once a year during the hurricanes season, starting at and ending at . Denote by the probability that a hurricane of intensity happens during season .
Depending on its intensity, on the type and state of the planted crops, and on the geological characteristics of a parcel, a hurricane happening at step
n can interrupt a production cycle and totally or partially destroy the crop growing on that parcel. Let
be the level of damage, with 0 and 1 corresponding to no damage to crops and total destruction of crops, respectively. We specify
, as follows:
where
E is a vector of fixed and known parameters describing the geological characteristics of the soil parcel, e.g., its slope and wind exposure;
is the age of the crop occupying the parcel at step
n; and,
b is the intensity of the hurricane it was exposed to.
If a hurricane strikes during a production cycle, then the revenues and the GISQ will be affected. The revenue loss will depend on the age of the crop at the time the hurricane hits and on the level of damage. Let
be the age of the crop, and
b the intensity of the hurricane, respectively. If no hurricane happens during the production cycle, then the parameter
b is set equal to zero. The emergency control
on the parcel at step
n will then be defined by
where
is the GISQ of the soil at the beginning of the agronomic cycle of
(when planted or sown). The new parcel’s GISQ value is given by
Further, cleaning and rehabilitation work have to be carried out on the farm before normal agricultural activities resume, which takes time and involves costs. The duration of the rehabilitation work after a hurricane that happened at step
n is defined by
The earnings generated by the emergency control
are defined by
that is, it depends on the hurricane intensity
, the crop type
and age
, the level of damage
, and the time elapsed since the last event
.
2.1. Dynamical System
We let the hurricane intensity
be a state variable, which adds a new event, e.g., a potential hurricane strike, to the only one considered in the deterministic model, which is, the end/start of a production cycle. The state of the farm (cash flow and soil quality) is determined by the vector
We recall that the control variable is
where
is the control applied to the parcel.
To account for crop seasonality, we define by
the set of possible controls on the parcel during a particular season
s.
The evolution of the farm is then governed by the following discrete-time dynamical system
:
The initial condition is , with , , , , , , , and are given and .
Equation (
3i) serves to update the event vector at each step. Equation (
3a) sets the clock on time at each step (
) by advancing it with
(the time elapsed between steps
n and
). Equation (
3b) does the same thing with the season, while taking into account its cyclicity. The other variables are updated according to the type of current event. At any step
n, if the remaining time until the next event is equal to the remaining time until the next hurricane strike (
), meaning that the event in question (at step
) is a hurricane, and then the hurricane intensity is generated following the probability distribution
(see (
3c)). The emergency control
has to be applied using (
3e). The evolution of the GISQ is given in (
3d). The remaining time until the next hurricane is set to the hurricane period of next year, since only one hurricane can occur per year (
3g). Otherwise, a production cycle ends on the parcel, leading to the modification of the soil quality following (
3d). If the production cycle that just ended is not the last one (
), then the crop and practice remain the same. We merely update the age of the crop and reset the value of the parameter
b to 0. We supposed that a crop with multiple production cycles (harvests) entirely recover after one production cycle. The effect of the hurricane is only felt on the harvest following the hurricane (see (
3e)), and a new production cycle starts right away, with the remaining time until the next event on the parcel being
(see (
3h)). Otherwise, if the whole agricultural cycle ends (interruption by total distruction from the hurricane or the end of the last production cycle of
), then a new crop and agricultural practice
has to be chosen from the set
for this parcel, as expressed in (
3e). The new time until the next event on the parcel (
) is equal to the length of the first production cycle of the new crop, which is,
see (
3h).
Finally, (
3f) updates the economic state of the farm by accumulating at each step the earnings generated by the parcel when a production cycle ends and any expenses resulting from a hurricane strike.
2.2. Viability Constraints and Capture Problem
As in Durand et al. [
18], we assume that the farmer aims at achieving the following two main objectives:
The first objective is economic. Clearly, it can take different forms depending on how one defines an “acceptable” income. One practical option is to retain a minimum threshold
, below which cash flow should not fall at any step. Consequently, we have the following admissible set:
The second objective is bio-ecological and it can be translated into a target to reach, i.e.,
An alternative to (5) could be to restore the soil quality to at least a certain level
of the initial value
, which is,
With the uncertainty that is related to climatic events, we end up with a stochastic capture problem, where we look for initial states of system
defined by (
3a)–(
3i), for which there exists at least one viable evolution in
(
4) until it reaches the target
(
5), or
(
6), with a probability exceeding a certain threshold
p (the confidence level). This set is called the stochastic capture basin and it is defined as follows:
3. The Emergency Control
The emergency control includes everything that must be done after a climatic event. The duration, cost, and earnings that result from the emergency control depend on the hurricane’s intensity and the damage it caused. In principle, a large variety of cases could occur, but, for tractability, we focus on two of them:
-
Total destruction of the crops. In this case, the climatic event destroys the totality of the crops on the parcel (
). If this happens, then the farmer has to carry out cleaning and rehabilitation work on the farm, and remove the damaged crops in order to be able to farm the parcel again. We suppose that the cleaning and rehabilitation work directly depends on the hurricane’s intensity. The corresponding time and cost are denoted
and
, respectively. The time and cost depend on the type of crop occupying the parcel
and its age
. Denote the time and cost by
and
, respectively. Consequently, the duration of the emergency control
is
The associated earnings correspond to the total financial flow generated during the time elapsed since the last update
minus the total amount incurred for the cleaning and rehabilitating of the parcel (
and
). The term
includes the fixed costs and the costs for maintenance, labour, etc., and any other nonrecoverable cost. Therefore, we have
-
Partial degradation of the crops. In this case, the climatic event only destroys a portion of the crops on the parcel (). Here, the farmer must decide, based on a given criterion, whether to bring the crop to maturity, or abandon it and start a new agricultural cycle instead. In practice, this decision will depend on the type of crop, its age, and the earnings it will potentially generate. To keep the model parsimonious, we consider a threshold above which degradation is considered to be total and the farmer decides to drop the crop. To be more realistic, one could let the threshold be a function of all other parameters.
If crop degradation caused by the hurricane is sufficiently low (
), then the emergency control must last as long as necessary for the surviving crops to reach maturity and for the cleaning work to be completed. Because these two processes evolve simultaneously, the duration of the emergency control
is given by
The earnings that are associated with
correspond to the gain from the surviving crops minus the total cost (cleaning cost, fixed and set-up costs). We suppose that the earnings from a production cycle are inversely proportional to its level of degradation
. Knowing that a totally healthy production cycle (
) generates earnings
where
represents the income from selling the harvest and
the harvesting cost, then the earnings from the surviving crop after a hurricane that caused a level
of degradation are equal to
Consequently, we have the following expression for the earnings from
:
Subsequently, the general form of the function
for
is given by
and for
by
Note that this expression of revenues is also valid for any other control than the emergency control. Indeed, in the absence of a hurricane, so and , which leads to , which is the earnings from a normally completed production cycle.
5. Solution Method
To compute the stochastic capture basin defined in (
7), we adopt a dynamic programming approach. We exploit the following property established in Proposition 1 and Corollary 1: if it is possible (not possible) to reach the target with a certain initial treasury, then it is possible (not possible) to reach it with a higher (lower) treasury.
The stochastic capture basin defined in Equation (
7) is the set of initial states
satisfying the constraints stated in (
7), which is, there exists an evolution starting from
that remains in the set
during the entire planning horizon and reaches the target
at the end of that time horizon with a probability greater than or equal to the confidence level.
As
,
,
, and
are given data, we can characterize any initial state
only by its
and
values, and write
Proposition 1. Let be an initial state from the stochastic capture basin ; then any initial state with belongs to . Proof. Let
be the minimum budget needed for the evolution
to be viable:
Therefore,
and consequently
□
Corollary 1. If it is not possible to reach the target from an initial state then it is not possible to reach it with less initial treasury. The results presented in Proposition 1 and Corollary 1 imply that, in order to completely characterize the capture basin, it is sufficient to identify its border or, in other words, to find for each initial GISQ the minimum initial treasury
needed for the restoration, i.e.,
Therefore, the problem amounts to finding a viable evolution of the system with the minimum possible restoration budget This problem can be formulated as a stochastic shortest-path problem, whose objective is to minimize the expected restoration budget, and which can be solved by dynamic programming.
A feedback solution takes the form of a decision tree that gives a control to apply for each state of the system from where it is possible to reach the target. Knowing that each sequence of controls applied in the past, combined with the realization of the uncertainties (the climatic events that occurred) during that period, leads the system to a specific state, the feedback solution gives a sequence of controls to apply during the next periods. Each branch of the decision tree representing the solution describes the evolution of the system for one possible realization of the uncertainties during the exploitation period.
6. Numerical Results
In this section, we consider some farming systems, to which we will refer as cases in a single-parcel setting.
Table 4 lists the crops considered in these cases. Each crop can be used either with a conventional or agroecological practice. Additionally, farmers can choose to leave their land fallow. When coupled with the agroecological practice (short fallow and long fallow), the land is left in simple fallow. When coupled with the conventional practice (improved short fallow and improved long fallow), the fallow is improved with the use of chemical additives or fertilizers.
Each one of the cases listed above represents a type of farming system practiced in the French West Indies: the export sector, specialized in banana and sugar cane (cases C1 and C2); the local market sector, based on diversified vegetable farming (C3 and C4); and, a theoretical case (C5) that combines all of the crops.
Computations were made for various hurricane impacts on the GISQ, initial GISQ levels, time horizons, confidence levels, and fixed costs in order to be able to analyze the effect of each parameter. We always require that the farm remain self-sufficient on average, i.e., we set the economic constraint to 0 () and impose that the retained solution(s) allow for a positive mean treasury.
All of the case studies concern the Guadeloupe archipelago. For the common elements, we use the same data and parameter values, as in Durand et al. [
18]. The other required data and parameter values are estimated in
Section 4 and/or displayed in
Appendix B and
Appendix C.
The rest of the section is divided into two parts. In the first part (
Section 6.1), we focus on the capture basins. In particular, we analyze the impact of hurricanes on the GISQ (in
Section 6.1.1), the impact of the time horizon (in
Section 6.1.2), the impact of crop diversification (in
Section 6.1.3), the impact of fixed charges (in
Section 6.1.4), and the impact of the initial GISQ (in
Section 6.1.5).
In the second part (
Section 6.2), we look at viable evolutions. More specifically, we discuss the main characteristics of such evolutions and how they are affected by such features as direct hurricane impact, time horizon, fixed charges, and GISQ improvement level.
6.1. Analysis of Capture Basins
Denote by the capture basin of Case i for a time horizon T, where the GISQ target is with a confidence level p. Concretely, the capture basin is the set of all couples of initial GISQ and budget (initial treasury) that make it possible to reach the target while respecting the imposed constraints.
Each panel in
Figure 2 represents a superposition of capture basins of the same case for different confidence levels. For example, in the upper left panel of
Figure 2, the surface in
represents the capture basin
whereas
represents the capture basin
.
In all cases, as expected, the capture basins get smaller when the confidence level is increased, which is,
Indeed, if it is possible to reach the target starting from a certain GISQ and using a certain budget with a confidence level then it is possible to reach it with a confidence level .
The higher is the initial GISQ , the lower is the minimum budget needed to reach the target. However, the budget is always strictly positive, even when the initial GISQ is larger than the target , due to the fixed costs and the money needed for planting the first crop. In fact, even if nothing is done (i.e., leaving the parcel in free fallow), the farmer needs to cover the fixed cost that must be paid, regardless of activity.
6.1.1. Direct Effect of Hurricanes on the GISQ
The GISQ is a quite recent index describing the soil quality and, thus, there are no studies establishing how its value is affected by hurricanes. To work around this lack of information, our numerical simulations retained the following three arbitrary contexts:
- 1—No effect:
regardless of its intensity, a hurricane’s impact on the GISQ is set equal to 0.
- 2—Moderate effect:
a minor hurricane decreases the GISQ value by , while a major hurricane reduces the GISQ by . This case is considered to be the default one, which is, when the effect is not specified in a figure or in the discussion, then it is a moderate effect.
- 3–Large effect:
as above, a minor hurricane decreases the GISQ by , but a major hurricane reduces the GISQ by
Figure 3 displays the capture basins for a series of cases, where the time horizon is 40 years and the GISQ target
for the three different hurricane effects introduced above. A first observation is that when a hurricane has no direct impact on the GISQ (see the first column of
Figure 3), then it is always possible to reach the desired level of soil restoration with probability 1, provided that the time horizon is long enough (see
Section 6.1.2 ) and a budget is available. Note that the confidence level only affects the minimum budget needed for restoration, due to the assumption that the climatic event does not deteriorate the GISQ.
Now, if a hurricane affects the GISQ, then we notice that the greater the impact, the more difficult it is to restore the soil to the desired level. In particular, for some realizations of the uncertain event, the cost of restoring the soil to the desired level is high, and restoration itself may become infeasible for certain confidence levels. To illustrate, for all displayed cases in
Figure 3, it is possible to restore the soil with an initial GISQ of
, with a confidence level of
, when we suppose that the hurricane have no direct impact on the GISQ. Restoration is also possible when we assume a moderate effect, but it would take a higher minimum budget. However, restoration is no longer possible if the hurricane has a large impact on the GISQ.
Finally, we note that the capture basin for a large hurricane impact on the GISQ is included in the capture basin with a moderate effect, which is, in turn, included in the no-effect capture basin.
6.1.2. Impact of the Time Horizon
Figure 4 and
Figure 5 display the capture basins for cases with different planning horizons (40, 20, 10, and five years) and levels of hurricane impact on the GISQ.
A first takeaway from these figures is that the capture basins shrink and slide to the right as the time horizon gets shorter in all cases, regardless of the level of hurricane impact on the GISQ. This indicates that the less time we have, the more difficult it is to restore the poorest soils, which is quite intuitive. For example, it is possible to bring a soil with a GISQ of to a GISQ of in 40 or 20 years, but not if we only have 10 or 5 years.
When the improvement is feasible, the minimum restoration budget, for different time horizons, varies across scenarios. To better visualize this behavior, we display in
Figure 6, a sample of the total restoration costs as a function of the time horizon for a soil restoration from GISQ
to GISQ
in the case of a moderate hurricane effect on the GISQ with a probability of
(upper part) and in the case of a bigger hurricane impact on the GISQ with a probability of
(lower part).
As we would expect, for the same conditions (same initial GISQ and same confidence level), the total restoration budget decreases with the time horizon; see case C4 and its variants (with fixed cost 400 and 800) as well as C3 and C5. Indeed, the less lime is available to restore the soil, the quicker and more efficient the work must be, and the higher is the cost. Moreover, there is a threshold time horizon above which the minimum budget stabilizes and becomes constant. In particular, this happens when the farm is profitable enough: when it just needs a certain amount to start the first production cycle and then it becomes self-sufficient. The time horizon at which the minimum budget stabilizes is the minimum time it takes for farming to become totally self-sufficient.
The minimum time horizon for self-sufficiency depends on the farm’s efficiency level. That is, when the farm faces fewer charges and/or achieves higher incomes, then the timespan for reaching self-sufficiency is shorter. For example, the minimum budget in C4 with a fixed cost of 100 stabilizes before those with fixed costs of 400 and 800. However, for other cases, e.g., C1 and C2, we observe the opposite, which is, globally, the longer the planning horizon, the more expensive it becomes to restore the soil, which seems to be pretty counterintuitive. In fact, this happens when the farm is not profitable enough, i.e., it does not generate sufficient cash to cover all of the operational expenses. Therefore, tautologically, the longer is the time horizon, the higher are the costs to be covered; as a result, the minimum restoration budget increases. This is confirmed by the treasury evolution in the different cases. In those cases with a decreasing restoration budget over time (the profitable cases), the treasury trend is increasing and ends up with a surplus. On the other hand, the cases with an increasing restoration budget over time exhibit a decreasing treasury trend and end up with an almost zero final treasury. To illustrate, we show in
Figure 7 the GISQ and treasury over the restoration time for two cases: a profitable one (C4) in the upper part of the figure and an unprofitable one (C1) in the lower part. Note that the two chosen evolutions are computed with the same realization of the uncertainty.
In terms of confidence level, one would expect that, if it is possible to reach the target with a given probability in a certain time horizon, then it would be possible to reach this target with that same probability in a longer time horizon. This intuition is confirmed for no or moderate hurricane impact on the GISQ (see the relevant parts in
Figure 4 and
Figure 5). However, this result does not hold when the impact of hurricanes on the GISQ is large. Indeed, the lower parts of
Figure 4 and
Figure 5 show that it is possible to reach the target from an initial GISQ of
with probability
when we have 5, 10, or 20 years available, but not when the planning horizon is 40 years. The reason is that when hurricanes’ impact on the GISQ is large, then the damage is also large; having a longer planning horizon does not help, as the costs accumulate and eventually exceed the restoration capability of the feasible controls in this case. Consequently, it becomes impossible to reach the target with the desired level of confidence.
6.1.3. Impact of Crop Diversification
The comparison of cases C1 and C2 to C3 and C4 in
Figure 4 and
Figure 5 reveals that the capture basins of scenarios with more crops are always at least as large as those with fewer crops, which is,
In fact, adding a control (a degree of freedom) can only help in meeting the objectives. The comparison of C5, a case with all crops, to the other cases confirms this observation. We can then conclude that, if the set of possible controls
in case
is included in the set
of controls in case
, then the capture basin of
is contained in that of
, i.e.,
Figure 6 shows that the restoration budget in a case with more crops (a larger set of possible controls) is at least as low as the restoration budget of a case with fewer crop choices, i.e.,
The above inclusion holds for all planning horizons and hurricane impacts on the GISQ.
It is true that crop diversification makes it easier to achieve the soil restoration objective, but adding new crops to the set of possible controls of a scenario does not necessarily decrease the restoration budget. If the initial set already contains the right mix of crops, adding new crops will certainly not deteriorate the solution, but neither will it improve it. For example, the performance of C5 is, in almost all cases, identical to the performance of C4, even if it has more crop choices.
6.1.4. Impact of Fixed Charges
Figure 8 displays the capture basins of C4 with different fixed costs (100, 400, and 800) and different hurricane impacts on the GISQ. First, we note that these results are consistent with the conclusions reached in
Section 6.1.1. Second, the overall shape of the capture basins remains the same when the fixed cost changes. This means that, if it is possible to reach the target from a certain initial GISQ with a certain probability, then it remains possible to reach it if the fixed cost is increased. However, the minimum budget for restoration significantly changes, i.e., it gets higher for cases with higher fixed costs. Further, the difference in budget when considering two fixed-cost values is more pronounced when the initial GISQ is low. Indeed, if we compare the minimum restoration budgets for the same level of improvement within the same time horizon and with the same probability (
Figure 9), we see that the budgets for cases with higher fixed costs are higher than those with lower fixed costs, and also that the budget increment is much larger when the initial GISQ is low. This difference can be explained by the fact that the crops planted in a better quality soil yield more and, therefore, generate more income, which can be used to cover a larger portion of the fixed costs.
Comparing the minimum restoration budget of these three cases for different time horizons (
Figure 6), we confirm that the minimum budget increases with the fixed cost. Additionally, we notice that the difference in budgets becomes smaller as the time horizon gets longer. The budget difference gradually decreases and stabilizes after a certain value of time horizon.
6.1.5. Impact of the Initial GISQ
Figure 9 displays the evolution of the minimum budget needed to raise the GISQ by
within 20 years and with a probability of
, as a function of the initial GISQ, for different hurricane impact levels on the GISQ. In a nutshell, a soil quality improvement of the same level does not necessarily require the same budget when the initial GISQ is low or high.
The first two plots in
Figure 9, corresponding to no and moderate hurricane impact on the GISQ, show that the higher the initial GISQ, the lower is the restoration cost. The reason is that a crop planted in better soil generates a higher income than a crop grown in poor-quality soil. Consequently, a lower minimum budget is needed, as the income covers a larger portion of the exploitation costs. Note that the budget may be roughly constant regardless of the initial GISQ value, as in C2. This occurs when the crop’s yield is not very sensitive to the GISQ value.
When the hurricane impact on the GISQ is large (lower part of
Figure 9), the above result remains valid for mid-interval initial GISQ levels. However, the opposite behavior is observed for values that are either close to 0 or 1. For instance, in C4 and C2, the minimum budget curves are increasing near initial GISQ values of 0 and 1. One explanation for this is that, when the GISQ is too low, the degradation caused by a hurricane, regardless of the storm’s severity, cannot be very high. On the other hand, starting with a soil of medium quality leaves room for improvements that exceed the target, in anticipation of future hurricanes and compensating for their damage. However, if the initial GISQ is too high, then such an opportunity is not available, which increases the farm’s sensitivity to hurricane strikes.
6.2. Viable Evolutions
In this section, we exhibit (viable) solutions with the minimum initial budget that lies on the boundary of capture basins. These solutions are not the result of an optimization of, e.g., income or another criterion, they are simply viable. In the case of multiple solutions, we chose the one with the highest final treasury. Further, we recall that the problem of determining viable evolutions is represented by a decision tree, with each branch corresponding to the evolution of the state for one possible realization of the uncertainties during the exploitation period.
Figure 10 displays one branch of the solution for case C4 over 40 years (
months), for a GISQ target
with an initial GISQ of
, and a confidence level
. The upper part of the figure shows the evolution of the GISQ and the treasury over the time horizon, while the lower part exhibits the sequence of controls to apply. Finally, the red vertical lines indicate the instants of the hurricane’s strikes and their intensity (“m” stands for a minor hurricane and “M” for a major one).
A first observation is that each increase in the treasury comes with a decrease of the GISQ level and vice versa. This illustrates the conflict between the criteria of farm profitability and resource preservation, which is at the root of our problem. The GISQ level goes from its initial value at time 0 to reach at the terminal date. The treasury starts at the minimum restoration budget, goes up and down over time, and remains positive throughout the planning horizon.
We make the following remarks on the viable evolutions and the effect of some parameters:
-
Crop recommendations: in all viable evolutions, there is a prevalence of fallow as a recommended control. This is due to the fact that the solutions lie on the boundary of the capture basin. Even if it were possible to obtain higher revenues by investing a bit more, the algorithm returns the solution with the minimum budget. As the cheapest option is fallow, this is most likely to be chosen.
-
Direct hurricane impact on the GISQ: when the GISQ is more sensitive to hurricane strikes, we once again notice that fallow controls are very prevalent. This is intuitive as they offer the best GISQ improvement to price ratio. When the impact of hurricanes on the GISQ is low, the degradation over time becomes less significant and thus we can afford the use of more profitable crops even if they are less performant on the GISQ improvement side.
-
Time horizon: the impact of the planning horizon depends on the hurricane impact on the GISQ. When this impact is low, fallows are less recommended when the horizon is longer. However, when the direct effect of hurricanes on the GISQ is high, then fallow controls are recommended more often, due to the large accumulation of damage on the GISQ over longer periods of time, as alluded to in
Section 6.1.2.
-
Fixed charges: higher fixed costs call for more profitable crops, that is, other than fallow. This is expected, as a higher income helps to cover this cost. In fact, it is more appealing to increase revenues than to keep budget to cover fixed costs. However, this result no longer holds when the hurricane impact on the GISQ is too high. Indeed, in this case, achieving the GISQ target becomes harder, if even possible, using productive crops. Consequently, there may be no other choice than fallow controls to improve the GISQ, which adds the fixed cost to the initial budget.
-
Level of improvement: the comparison of
Figure 10,
Figure 11 and
Figure 12, which respectively show a branch of the feedback solution tree for an initial GISQ of
,
, and 1, shows that a higher level of improvement in the GISQ involves a fallow control. The intuition is crystal clear: a small level of improvement in the GISQ requires less effort and time to achieve, which leaves more time to improve the treasury.