Next Article in Journal
Challenges for Protected Areas Management in China
Previous Article in Journal
Process Planning in Industry 4.0—Current State, Potential and Management of Transformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Viability of Agroecological Systems under Climatic Uncertainty

by
Aïchouche Oubraham
1,
Patrick Saint-Pierre
2 and
Georges Zaccour
1,3,*
1
GERAD, HEC Montréal, 3000 Chemin de la Côte-Sainte-Catherine, Montreal, QC H3T 2A7, Canada
2
LASTRE, 14 rue Domat, F-75005 Paris, France
3
Chair in Game Theory and Management, 3000 Chemin de la Côte-Sainte-Catherine, Montreal, QC H3T 2A7, Canada
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(15), 5880; https://doi.org/10.3390/su12155880
Submission received: 6 June 2020 / Revised: 17 July 2020 / Accepted: 20 July 2020 / Published: 22 July 2020
(This article belongs to the Section Sustainable Agriculture)

Abstract

:
To cope with ever-increasing demand and ensure food security, agronomic systems have shifted over time from traditional agriculture, based on the organic fertilization of soils, to intensive and specialized farming that use chemical fertilization. This resulted in increased soil productivity in the short term, but caused serious ecological drawbacks over time (degradation of soil quality, pollution of water and air, loss of biodiversity, erosion, etc.), and even reversed the trend of agricultural productivity. In this paper, we propose a viability theory–based model to study the sustainability of an agricultural system subject to climate uncertainty. Our objective is to determine what farming practices and activity sequences restore soil quality to a desired level while ensuring an acceptable level of productivity in the presence of the risk of major climatic disasters. The model is applied to Guadeloupe, an island in the West French Indies. We found that the results are highly sensitive to the direct effect of hurricanes on the soil’s quality, which, in turn, strongly affects the impact of the other parameters and that the export oriented sector is more vulnerable and less resilient to climatic uncertainties than the sector aimed at the local market.

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 τ T = { 0 , . . . , T } . The state of the farm is described by two variables: (i) the cash flow W R , 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 0 , 1 and is based on 54 variables measuring these characteristics (Velásquez et al. [19] and Camacho et al. [20]). Let I be the GISQ value of the parcel, with I I = { 0 , . . . , N I 1 } , which is, I can take a finite number of values.
The evolution of the soil quality depends on its current value I and on the following two control variables: (i) the crop σ grown on the parcel, with σ Σ = { σ 1 , . . . , σ N σ } ; and, (ii) the agricultural practice π Π = { π 1 , . . . , π N π } . 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:
ϕ : I × Σ × Π I , ( I , σ , π ) ϕ ( I , σ , π ) .
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:
ψ : Σ × Π 0 , δ ( σ , π ) , ( σ , π ) ψ ( σ , π ) .
Moreover, to each crop and practice σ , π , we associate its first harvest time f c σ , π and the time duration between two successive harvests p c σ , π . The duration of the first production cycle is f c σ , π , and the following cycles last p c σ , π . If a crop has a single harvest, then f c σ , π = δ ( σ , π ) and p c σ , π = 0 . 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
: I × Σ × Π R , ( I , σ , π ) ( I , σ , π ) .
Denote by s S = { 1 , . . . , 12 } 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 Σ ( s ) Σ 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 τ ( n ) to its timing, then the next event will happen at
τ ( n + 1 ) = τ ( n ) + ψ ( σ n , π n ) .
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:
  • e 1 : the remaining time before the next end/start of a production cycle,
  • e 0 : the remaining time before the next potential hurricane strikes, and
  • ϵ = min { e 0 , e 1 } : the remaining time before the next change in the system’s state,
Denote by e = ( e 0 , e 1 , ϵ ) the vector of events.
Denote by H ( t ) the intensity of a climatic event at time t. We suppose that H ( t ) takes its values from a finite and discrete set H = { 0 , 1 , . . . , N h } , where 0 represents the absence of a climatic event and N h is the highest possible intensity. In the case studies, we assume that H = { 0 , 1 , 2 } , 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 β S and ending at γ S . Denote by P ( s , h ) the probability that a hurricane of intensity h H happens during season s S .
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 θ ( n ) [ 0 , 1 ] be the level of damage, with 0 and 1 corresponding to no damage to crops and total destruction of crops, respectively. We specify θ ( n ) , as follows:
θ ( n ) = Θ ( b , σ ( n ) , a ( n ) , E ) ,
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; a ( n ) 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 a ( n ) 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 v ( n ) on the parcel at step n will then be defined by
v ( n ) = ( σ ( n ) , π ( n ) , a ( n ) , b , c ) Σ × Π × [ 0 , δ m a x ] × H × [ 0 , 1 ] .
where c [ 0 , 1 ] 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
I ( n + 1 ) = ζ H ( n ) , I ( n ) , v ( n ) , ϵ ( n ) , E .
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
T ( v ( n ) ) = f ( H ( n ) , σ ( n ) , a ( n ) , θ ( n ) ) .
The earnings generated by the emergency control v ( n ) are defined by
L ( I , v ( n ) , ϵ ( n ) ) = g ( H ( n ) , σ ( n ) , a ( n ) , θ ( n ) , ϵ ( n ) ) ,
that is, it depends on the hurricane intensity H ( n ) , the crop type σ ( n ) and age a ( n ) , the level of damage θ ( n ) , and the time elapsed since the last event ϵ ( n ) .

2.1. Dynamical System

We let the hurricane intensity H ( t ) 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
x = ( I , W , H , e , s , τ ) .
We recall that the control variable is
v = ( u , a , b , c ) V ,
where u = ( σ , π ) U = Σ × Π is the control applied to the parcel.
To account for crop seasonality, we define by
U ( s ) U , U ( s ) = Σ ( s ) × Π ,
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 G :
G = { (3a) τ ( n + 1 ) = τ ( n ) + ϵ ( n ) (3b) s ( n + 1 ) = ( s ( n ) + ϵ ( n ) ) m o d   12 , (3c) H ( n + 1 ) = { Generated following P , if ϵ ( n ) = e 0 ( n ) , 0 , otherwise , (3d) I ( n + 1 ) = ζ ( I ( n ) , v ( n ) , ϵ ( n ) , E ) , (3e) v ( n + 1 ) { { ( σ ( n ) , π ( n ) , a ( n ) + ϵ ( n ) , H ( n + 1 ) , c ( n ) ) } , if ϵ ( n ) = e 0 ( n ) , { ( σ ( n ) , π ( n ) , a ( n ) + ϵ ( n ) , 0 , c ( n ) ) } if ϵ ( n ) = e 1 ( n ) and θ ( n ) < θ m a x and a ( n ) < δ ( σ ( n ) , π ( n ) ) , { ( σ , π , 0 , 0 , I ( n + 1 ) ) | ( σ , π ) U ( s ( n + 1 ) ) } otherwise , (3f) W ( n + 1 ) = W ( n ) + L ( I ( n ) , v ( n ) , ϵ ( n ) ) , (3g) e 0 ( n + 1 ) = { 12 ( s ( n ) β ) If H ( n ) > 0 Or s ( n ) = γ , 1 If H ( n ) = 0 And β s ( n ) < γ , e 0 ( n ) ϵ ( n ) Otherwise (3h) e 1 ( n + 1 ) = { T ( v ( n ) ) , if ϵ ( n ) = e 0 ( n ) , p c ( σ ( n ) , π ( n ) ) , if ϵ ( n ) = e 1 ( n ) and θ ( n ) < θ m a x and a ( n ) < δ ( σ ( n ) , π ( n ) ) , f c ( σ ( n + 1 ) , π ( n + 1 ) ) otherwise , (3i) ϵ ( n + 1 ) = min { e 0 ( n + 1 ) , e 1 ( n + 1 ) }
The initial condition is x ( 0 ) = ( I ( 0 ) , W ( 0 ) , H ( 0 ) , e ( 0 ) , s ( 0 ) , τ ( 0 ) ) , with τ ( 0 ) = 0 , s ( 0 ) = 1 , e 0 ( 0 ) = β 1 , H ( 0 ) = 0 , e 1 ( 0 ) = 0 , ϵ ( 0 ) = min i = 0 , 1 ¯ e i ( 0 ) , δ ( 0 ) = 0 , I ( 0 ) and W ( 0 ) are given and c ( 0 ) = I ( 0 ) .
Equation (3i) serves to update the event vector at each step. Equation (3a) sets the clock on time at each step ( n + 1 ) by advancing it with ϵ ( n ) (the time elapsed between steps n and n + 1 ). 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 ( ϵ ( n ) = e 0 ( n ) ), meaning that the event in question (at step n + 1 ) is a hurricane, and then the hurricane intensity is generated following the probability distribution P (see (3c)). The emergency control { ( σ ( n ) , π ( n ) , a ( n ) + ϵ ( n ) , H ( n + 1 ) , c ( n ) ) } 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 ( a ( n ) < δ ( σ ( n ) , π ( n ) ), 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 p c ( σ ( n ) , π ( n ) ) (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 σ ( n ) ), then a new crop and agricultural practice v ( n + 1 ) has to be chosen from the set U ( s ( n + 1 ) ) for this parcel, as expressed in (3e). The new time until the next event on the parcel ( e 1 ( n + 1 ) ) is equal to the length of the first production cycle of the new crop, which is, f c ( σ ( n + 1 ) , π ( n + 1 ) ) ; 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:
  • Ensuring an acceptable income throughout the whole planning horizon.
  • Restoring the quality of the soil to an acceptable level I * by the end of the planning horizon.
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 W m i n , below which cash flow should not fall at any step. Consequently, we have the following admissible set:
K = { ( I , W , H , e , s , τ ) τ T , W W m i n 0 } .
The second objective is bio-ecological and it can be translated into a target to reach, i.e.,
C = { ( I , W , H , e , s , τ ) τ T , I I * } .
An alternative to (5) could be to restore the soil quality to at least a certain level d % of the initial value I ( 0 ) , which is,
C d = { ( I , W , H , e , s , τ ) τ T , I min { d 100 I ( 0 ) ; 1 } } .
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 G defined by (3a)–(3i), for which there exists at least one viable evolution in K (4) until it reaches the target C (5), or C d (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:
C a p t G , P ( K , C ) = { X 0 K x ( . ) s . t . x ( 0 ) = X 0 & t > 0 s . t . P ( x ( t ) C & n [ 0 , t ] , x ( n ) K ) p } .

3. The Emergency Control

The emergency control v = ( σ , π , a , b , c ) 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 ( θ ( n ) = 1 ). 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 t c ( H ( n ) ) and c c ( H ( n ) ) , respectively. The time and cost depend on the type of crop occupying the parcel σ ( n ) and its age a ( n ) . Denote the time and cost by t r ( σ ( n ) , a ( n ) ) and c r ( σ ( n ) , a ( n ) ) , respectively. Consequently, the duration of the emergency control v = ( σ , π , a , b , c ) is
T ( v ( n ) ) = t c ( b ) + t r ( σ ( n ) , a ) .
The associated earnings correspond to the total financial flow generated during the time elapsed since the last update ρ ( I ( n ) , v ( n ) , ϵ ( n ) ) minus the total amount incurred for the cleaning and rehabilitating of the parcel ( c c ( H ( n ) ) and c r ( σ ( n ) , a ( n ) ) ). The term ρ ( I ( n ) , v ( n ) , ϵ ( n ) ) includes the fixed costs and the costs for maintenance, labour, etc., and any other nonrecoverable cost. Therefore, we have
L ( I , v , ϵ ) = ρ ( I , v , ϵ ) c c ( b ) + c r ( σ , a ) .
     
Partial degradation of the crops. In this case, the climatic event only destroys a portion of the crops on the parcel ( θ ( n ) < 1 ). 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 θ m a x 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 θ m a x be a function of all other parameters.
If crop degradation caused by the hurricane is sufficiently low ( θ ( n ) < θ m a x ), 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 v = ( σ , π , a , b , c ) is given by
T ( v ( n ) ) = M a x { 1 { t c ( b ) δ ( σ ( n ) , π ( n ) ) a ( n ) } t c ( b ) , e 1 ( n ) ϵ ( n ) } .
The earnings that are associated with v = ( σ , π , a , b , c ) 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 ( θ = 0 ) generates earnings
( I ( n ) , σ ( n ) , π ( n ) ) = m ( I ( n ) , σ ( n ) , π ( n ) ) h c ( σ ( n ) , π ( n ) ) ,
where m ( I ( n ) , σ ( n ) , π ( n ) ) represents the income from selling the harvest and h c ( σ ( n ) , π ( n ) ) the harvesting cost, then the earnings from the surviving crop after a hurricane that caused a level θ of degradation are equal to
( 1 θ ) m ( I ( n ) , σ ( n ) , π ( n ) ) h c ( σ ( n ) , π ( n ) ) .
Consequently, we have the following expression for the earnings from v = ( σ , π , a , b ) :
L ( I , v , ϵ ) = ( 1 θ ) ( I , σ , π ) θ h c ( σ , π ) + ρ ( I , v , ϵ ) c c ( b ) .
Subsequently, the general form of the function T ( v ( n ) ) for v ( n ) = ( σ , π , a , b , c ) is given by
T ( v ( n ) ) = t c ( b ) + t r ( σ ( n ) , a ) , if θ ( n ) θ m a x M a x { 1 { t c ( b ) δ ( σ ( n ) , π ( n ) ) a ( n ) } t c ( b ) , e 1 ( n ) ϵ ( n ) } , otherwise
and for L ( I , v , ϵ ) by
L ( I , v , ϵ ) = ρ ( I , v , ϵ ) c c ( b ) + c r ( σ , a ) , if θ θ m a x , ( 1 θ ) ( I , σ , π ) θ h c ( σ , π ) + ρ ( I , v , ϵ ) c c ( b ) , otherwise .
Note that this expression of revenues is also valid for any other control than the emergency control. Indeed, in the absence of a hurricane, b = 0 so θ ( n ) = 0 and c c ( b ) = 0 , which leads to L ( I , v , ϵ ) = ( I , σ , π ) + ρ ( I , v , ϵ ) , which is the earnings from a normally completed production cycle.

4. The Climatic Events

This section presents all of the relevant information regarding the climatic events and their related uncertainties.

4.1. Types of Climatic Events

A tropical cyclone is a rotating weather system that forms over the ocean when specific water temperature and air humidity conditions arise simultaneously. A significant amount of ocean water then evaporates quickly, causing an atmospheric disturbance. This results in winds revolving around an axis, under the effect of the earth’s rotation, ultimately causing very strong winds, spiral thunderstorms, and heavy rain. When the maximum sustained surface winds of such formations are less than 39 miles per hour (mph), the cyclone is called a tropical depression, and between 39 and 74 mph, it is called a tropical storm. Beyond 74 mph, the event is referred to as a hurricane if it forms over the North Atlantic or Northeast Pacific Ocean, a typhoon if it forms over the Northwest Pacific Ocean, and simply a cyclone if it forms over the South Pacific or Indian Ocean (Source: NOAA website, Wikipedia, and others.)
The intensity of a hurricane is measured by the Saffir–Simpson Hurricane Wind Scale (SSHWS). It has a 1 to 5 categorization based on the hurricane’s intensity at an indicated time (see Table 1 Source: https://www.nhc.noaa.gov/aboutsshws.php.) We talk about a major hurricane for events of categories 3 to 5.

4.2. Climatic Uncertainty

As the Caribbean Basin is one of the most hurricane-prone regions in the world, many meteorological studies are available. Surprisingly, however, the French West Indies seem to be an exception, with relatively few studies.
One of the few dedicated studies of this region is Garnier et al. [22]. The authors obtained data from the French archives in France and the islands of Guadeloupe and Martinique to survey all f the hurricanes that hit the French West Indies between 1635 and 2007. They used the data to build a long-term chronology of hurricane severity and damage in this area and to draw some statistical conclusions from these events, in particular yearly probabilities of hurricanes by category, which are of direct interest to our study. We make use of the results of this work because of its geographical relevance and the reliability of its data sources. In addition, the authors only included hurricanes that actually hit the islands of Guadeloupe and Martinique, which makes the results extremely relevant geographically. Table 2 presents the yearly hurricane probabilities extracted from [22], which will be used in our numerical illustrations.

4.3. Hurricane Effects

The effects of a hurricane on soil and agriculture depend strongly on the type of agriculture and crops, and the soil’s type, relief, etc. Strong winds may affect the physical characteristics of the soil itself. Strong rainfall may induce flooding that can change the chemical composition of the soil and cause landslides in sloped areas. Ocean wind pressure may result in storm surges that can cause severe salt contamination of the soil in costal areas. All of these factors may disturb the microfauna and ultimately change the biological characteristics of the soil ([23]).
There are currently no studies estimating the effect of hurricanes on the GISQ. Therefore, in our numerical applications, we will look at three arbitrary situations: no effect, moderate effect, and significant effect of a hurricane on the GISQ (the precise meaning of these effects is given in Section 6.1.1).
Strong winds can cause significant structural damage to or uproot crops, and flooding may cause certain crops (especially those growing underground) to gorge with water or to rot, and changes to soil composition can compromise the normal growth of the crops. Spencer and Polachek [24] studied the effects of hurricanes on local crop production in Jamaica and found that crops grown above ground suffer greater damage than those grown below ground. In fact, the only below-ground crops that experience a drop of productivity are yams and potatoes. The vast majority of hurricanes that were considered in this study were minor (category 1 and 2) and their effects on crops are displayed in the first row of Table 3.
We consider, not unrealistically, that the losses are total after a major hurricane strike. Accordingly, our choice to consider a degradation of more than 60% ( θ m a x = 0.6 ) as a total loss will have no impact on the results as studies indicate that hurricanes of category 3 or more cause at least an 80% loss of production.
Table 3 displays the percentage loss in productivity caused by certain types of climatic events on three categories of crops: the above-ground category, which includes all the crops growing above ground; the water-damaged category, which includes above-ground crops that are vulnerable to flooding; and the below-ground category, which includes crops growing underground that are flood resistant.

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
X 0 = ( I ( 0 ) , W ( 0 ) , H ( 0 ) , e ( 0 ) , s ( 0 ) , τ ( 0 ) ) ,
satisfying the constraints stated in (7), which is, there exists an evolution starting from X 0 that remains in the set K during the entire planning horizon and reaches the target C at the end of that time horizon with a probability greater than or equal to the confidence level.
As τ ( 0 ) = 0 , s ( 0 ) = 1 , H ( 0 ) = 0 , and e ( 0 ) = 0 are given data, we can characterize any initial state X 0 only by its I ( 0 ) and W ( 0 ) values, and write
X 0 = ( I ( 0 ) , W ( 0 ) ) .
Proposition 1.
Let ( I , W ) be an initial state from the stochastic capture basin C a p t G , P ( K , C ) ; then any initial state ( I , W ) with W W belongs to C a p t G , P ( K , C ) .
( I , W ) C a p t G , P ( K , C ) ( I , W ) C a p t G , P ( K , C ) W W .
Proof. 
Let B ( x , I , p ) be the minimum budget needed for the evolution x ( . ) to be viable:
B ( x , I , p ) = min { W x ( 0 ) = ( I , W ) & t > 0 s . t . P ( x ( t ) C & n [ 0 , t ] , x ( n ) K ) p }
We have,
( I , W ) C a p t G , P ( K , C ) x ( . ) s . t . B ( x , I , p ) W .
Therefore,
W W , B ( x , I , p ) W ,
and consequently
( I , W ) C a p t G , P ( K , C ) .
 □
Corollary 1.
If it is not possible to reach the target from an initial state ( I , W ) , then it is not possible to reach it with less initial treasury.
( I , W ) C a p t G , P ( K , C ) ( I , W ) C a p t G , P ( K , C ) W W
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 W i n f ( I , p ) needed for the restoration, i.e.,
W i n f ( I , p ) = inf ( I , W ) C a p t G , P ( K , C ) W = min x B ( x , I , p ) .
Therefore, the problem amounts to finding a viable evolution of the system x ( · ) with the minimum possible restoration budget B ( x , I , p ) . 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 ( W m i n = 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 C a p t C i , T ( I * , p ) the capture basin of Case i for a time horizon T, where the GISQ target is I * with a confidence level p. Concretely, the capture basin is the set of all couples of initial GISQ I 0 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 Sustainability 12 05880 i001 represents the capture basin C a p t C 1 , 40 ( 0.8 , 0.9 ) , whereas Sustainability 12 05880 i002 represents the capture basin C a p t C 1 , 40 ( 0.8 , 0.2 ) .
In all cases, as expected, the capture basins get smaller when the confidence level is increased, which is,
C a p t C i , T ( I * , p 1 ) C a p t C i , T ( I * , p 2 ) , p 1 p 2 .
Indeed, if it is possible to reach the target starting from a certain GISQ and using a certain budget with a confidence level p 1 , then it is possible to reach it with a confidence level p 2 p 1 .
The higher is the initial GISQ I 0 , the lower is the minimum budget needed to reach the target. However, the budget is always strictly positive, even when the initial GISQ I 0 is larger than the target I * , 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 0.01 , while a major hurricane reduces the GISQ by 0.02 . 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 0.01 , but a major hurricane reduces the GISQ by 0.05 .
Figure 3 displays the capture basins for a series of cases, where the time horizon is 40 years and the GISQ target I * = 0.8 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 0.2 , with a confidence level of 0.6 , 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 0.2 to a GISQ of 0.8 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 0.6 to GISQ 0.8 in the case of a moderate hurricane effect on the GISQ with a probability of 0.8 (upper part) and in the case of a bigger hurricane impact on the GISQ with a probability of 0.5 (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 0.8 with probability 0.8 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,
C a p t C 1 , T ( I * , p ) C a p t C 2 , T ( I * , p ) and C a p t C 3 , T ( I * , p ) C a p t C 4 , T ( I * , p ) .
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 U i in case C i is included in the set U j of controls in case C j , then the capture basin of C i is contained in that of C j , i.e.,
U i U j C a p t C i , T ( I * , p ) C a p t C j , T ( I * , p ) .
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.,
U i U j B u d g e t C i , T ( I * , p ) B u d g e t C j , T ( I * , p )
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 0.2 within 20 years and with a probability of 0.8 , 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 ( T = 480 months), for a GISQ target I * = 0.8 with an initial GISQ of 0.6 , and a confidence level p = 0.8 . 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 I 0 at time 0 to reach I T I * 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 0.6 , 0.1 , 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.

7. Conclusions

The main qualitative takeaways of our study are as follows:
  • The results seem to be highly sensitive to the direct effect of hurricanes on the GISQ, which in turn strongly affects the impact of the other parameters.
  • Cases C1 and C2 are globally less efficient than the others, which suggests that the export-oriented sector is more vulnerable and less resilient to climatic uncertainties than the sector aimed at the local market. This result is quite worrisome, given that the export-oriented sector represents a fairly large portion of the economy of small tropical islands such as Guadeloupe and Martinique.
  • Because it conflicts with profitability and food-security objectives, preserving the soil is far from an easy task. This difficulty is accentuated by climatic uncertainties, which are the most damaging hazards to soil preservation, not only in terms of cost, but also due to their direct impact on the soil itself. Our results highlight that restoration possibilities and cost are very sensitive to the direct impact of hurricanes on soil quality. Consequently, it is crucial to have studies assessing the impact of hurricanes on the GISQ (and hence on soil quality).
Future developments include collecting additional data and refining the estimations of hurricane-related parameters (impact on the GISQ, recovery time, cost, and occurrence probabilities). A second interesting extension would be to consider a multi-parcel setting to assess the potential advantages of multi-cropping, both in terms of efficiency and risk management. Finally some cropping systems, such as monocultures, often lead to gradual soil degradation, but this can be decreased or even reversed to a certain extent with no or minimum tillage, cover cropping and organic amendment strategies. It would be then interesting to take into account, in future research, the effect of successive use of the same crop or type of crop *.
* We are grateful to one Reviewer for suggesting this extension.
** We would like to thank three anonymous Reviewers for their very helpful comments.

Author Contributions

Data curation, A.O.; Formal analysis, A.O. and G.Z.; Methodology, A.O., P.S.-P. and G.Z.; Resources, P.S.-P.; Supervision, P.S.-P. and G.Z.; Validation, P.S.-P.; Writing—original draft, A.O.; Writing—review & editing, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

NSERC Canada, grant RGPIN-2016-04975.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Lexicon

  
Capture basin: The set of all initial states of the system from which it exists at least one evolution that satisfies the viability constraints and reaches the target at the end of the time horizon with a certain probability.
  
Agricultural practice/type of agriculture: The type of practices adopted for the agricultural activities. We considered two types: conventional and agroecological.
  
Conventional practice: Agricultural practice based on modern practices/tools and on chemical fertilization.
  
Agroecological practice: Traditional practice based on organic fertilization of the soil or other practices for improved management of agroecosystems. Wezel et al. [21]
  
Fallow: Leaving land fallow consists in leaving it unsown for a certain period of time. We considered two types of fallow: the one coupled with the conventional practice (conventional fallow) and the one coupled with the agroecological practice (agroecological fallow).
  
Agroecological fallow: (short fallow and free long fallow) the land is left in simple fallow.
  
Conventional fallow: (improved short fallow and improved long fallow) the fallow is improved with chemical additives and fertilizers.

Appendix B. Assumptions

We made the following assumptions in our study:
A1: 
There can be at most one hurricane per year in September ( β = γ = 9 ).
A2: 
The probabilities of yearly hurricane occurrence are given and Markovian (independent of the hurricane occurrence history).
A3: 
Due to lack of information about hurricanes, we chose arbitrary values for some parameters, that is,
  • The rehabilitation time after a hurricane is 2 months after a minor hurricane and 6 months after a major one. (Since the resilience of crops is not at the center of our study, these parameters are chosen such that only the crops with multiple harvest and whose degradation is not total can recover and that the effects of a hurricane are only felt at the harvest following it)
  • The rehabilitation costs after a hurricane are set to 0, i.e., these costs are covered by insurance policies whose price is already included in the fixed costs.
  • The degradation caused by a hurricane on a crop with a single harvest is nonreversible, while a crop with multiple harvests recovers after the harvest that follows the hurricane strike.
A4: 
When a production cycle is interrupted by a hurricane, the effect of the crop on the GISQ is proportional to the crop’s age (to the time the crop spent on the parcel).
A5: 
The maximum degradation of all the crops is θ m a x = 60 % (the one from which the farmer decides to replace the crop). With the data considered in the numerical applications corresponds to the degradation caused by a major hurricane.
A6: 
The direct hurricane impact on the GISQ can have three levels: no effect, moderate effect, and large effect.

Appendix C. Parameter Values and Data

Appendix C.1. Agronomic Data and Parameters

For any control ( σ , π ) Σ × Π , we have listed the agronomic parameters and transition functions in this section. The values of all the parameters used in this section are displayed in the tables of Appendix C.3.

Appendix C.1.1. Crops Yield

  • R M = R M ( σ , π ) : Mean yield (tons/Ha). Corresponds to the mean yield of the crop when the soil is of average quality (GISQ = 0.5 ).
  • r i = r i ( σ , π ) : The control’s sensitivity to the quality of the soil.
  • R ( I , σ , π ) : The yield of the control in the absence of hurricanes.
If r i 0 . 5 the yield is given by
R ( I , σ , π ) = 2 R M I ( 1 2 r i ) + r i ( r i ) 2 + 2 I ( 1 2 r i ) 2 r i 1 If I 0 . 5 , R M 1 + 2 r i ( 1 2 r i ) ( I 0 . 5 ) ( 1 r i ) + ( 1 r i ) 2 2 ( I 0 . 5 ) ( 1 2 r i ) 2 r i 1 Otherwise
If r i = 0 . 5 the yield is given by
R ( I , σ , π ) = 2 I R M If I 0 . 5 , R M ( I + 0 . 5 ) Otherwise
  • R E ( I , σ , π ) : The effective yield of the control.
    R E ( σ , π ) = 0 If θ θ m a x , ( 1 θ ) R ( I , σ , π ) Otherwise

Appendix C.1.2. GISQ Transition Functions

  • r d : The sensitivity of the soil to the loss of organic matter by crop.
  • r p : The damage or improvement to the soil caused by the agricultural practice (effects that are dependent on the initial soil quality).
  • r a p : The damage or improvement to the soil caused by the agricultural practice (effects that are independent of the initial soil quality).
  • Δ I b ( I , σ , π ) : The change in GISQ induced by the crop σ
    Δ I b ( I , σ , π ) = R ( I , σ , π ) 2 R M ( σ , π ) r d ( σ , π ) If R M ( σ , π ) 0 , 0 Otherwise
  • Δ I p ( I , σ , π ) : The change in GISQ induced by the agricultural practice π
    Δ I p ( I , σ , π ) = I . r p ( σ , π ) + r a p ( σ , π )
  • ϕ ( I , σ , π ) : The transition function of the GISQ in the absence of hurricanes.
    ϕ ( I , σ , π ) = min { max { I Δ I b ( I , σ , π ) + Δ I p ( I , σ , π ) , 0 } , 1 }
  • δ ( σ , π ) : The cycle length of control ( σ , π ) (the normal cycle duration when not interrupted).
  • D ( h ) : The direct impact of a hurricane of intensity h on the GISQ.
  • ζ ( I , v , ϵ , E ) : The transition function of the GISQ if a hurricane of intensity b hits a parcel on which the control ( σ , π ) is applied when the time elapsed since the last change in the GISQ is ϵ ( v = ( σ , π , a , b , c ) ).
    ζ ( I , v , ϵ , E ) = min { max { I + ϵ δ ( σ , π ) ϕ ( c , σ , π ) D ( b ) , 0 } , 1 }

Appendix C.2. Economic Data and Parameters

The parameters values used in this section are displayed in the tables of Appendix C.3.
For any control ( σ , π ) Σ × Π , we have the following economic parameters and transition functions:
  • S a : Subsidy for the crop (€/month).
  • S p : Subsidy for the quantity produced (€/ton).
  • C t : Monthly labor cost (€/ton).
  • C i : Installation input cost (seeds, fertilizer, pesticides) (€/Ha).
  • C e : Other expenses (€/Ha).
  • C c : Harvesting cost (€/Ha).
  • C m : Maintenance cost (€/Ha/month).
  • C f : Fixed cost that are not related to a particular crop (€/Ha/month).
  • P: Selling price (€/ton).
  • f c : Date of the first sale.
  • p c : The time between successive sales.
  • ρ 0 : The monthly recurring payoff.
    ρ 0 = S a C f C t C m
  • ρ ( I , v , ϵ ) : The cash flow generated during the time that elapsed since the last event.
    ρ ( I , v , ϵ ) = C i , if a = 0 , ϵ ρ 0 , otherwise ,
  • ( I , σ , π ) : The earning generated by each harvest sale.
    ( I , v , ϵ ) = m ( I , σ , ϵ ) h c ( I , v , ϵ )
where,
  • m ( I , σ , ϵ ) : The income from selling the harvest.
    m ( I , σ , ϵ ) = R E ( I , σ , π ) ( P + S p ) if a = f c or ( a f c ) mod P c = 0 , 0 otherwise ,
  • h c ( I , v , ϵ ) : The harvesting cost.
    h c ( I , v , ϵ ) = C c R E ( I , σ , π ) R M ( σ , π ) if a = f c or ( a f c ) mod p c = 0 , C c R E ( I , σ , π ) R M ( σ , π ) + C e if a δ ( σ , π ) , 0 otherwise ,

Appendix C.3. Crop Data

All the parameters are normalized for one parcel of unit area (1 ha); production data are in tons per ha; economic data are provided by experts in euros per ha or per tons produced (fixed and variable input costs, depreciations, subsidies, labor, sale prices) or per month (labor). Data are extracted from Durand et al. [18] and displayed in the tables below.
Table A1 displays the planting/sowing seasons of the seasonal crops.
Table A2 displays the values of the agronomic parameters and Table A3 those of the economic parameters.
Table A1. Planting months of crops.
Table A1. Planting months of crops.
Crop010203040506070809101112
Plantain xxxxx
Export Banana xxxxx
Eggplant xx
Tomato xx
Table A2. Agronomic parameter values.
Table A2. Agronomic parameter values.
Crop and Practice δ ( σ , π ) r i r d r p r ap f c p c R M
Export banana Agro600.750.20.120.02121224
Export banana Conv600.30.5−0.080121230.6
Plantain Agro600.70.190.150635
Plantain Conv240.250.45−0.02−0.01637.85
Cabbage Agro30.50.210.1503013
Cabbage Conv30.250.3−0.0203018
Carrot Agro40.650.140.0704011
Carrot Conv40.050.15−0.03−0.0054014
Cassava Agro120.50.150.07012022
Cassava Conv120.10.2−0.01−0.0112025
Cucumber Agro20.80.280.1502013
Cucumber Conv20.050.25−0.05−0.0252018
Eggplant Agro70.80.20.1207012
Eggplant Conv70.150.2−0.02−0.0757020
Green beans Agro40.60.090.010.0024010
Green beans Conv30.050.1−0.0103013
Lettuce Agro20.60.180.120207
Lettuce Conv20.30.25−0.0102013
Melon Agro30.70.170.10309
Melon Conv20.20.35−0.02−0.0152016
Sugar cane Agro600.30.080.080.035121260.5
Sugar cane Conv600.10.05−0.01−0.01121263.75
Tomato Agro50.80.20.130.0025011
Tomato Conv40.10.4−0.05−0.064015
Turban squash Agro30.70.170.09903014
Turban squash Conv30.10.15−0.05−0.023022
Yam (yellow) Agro80.60.180.110.0028011.7
Yam (yellow) Conv80.150.4−0.0108015
Yam (Grosse Caille) Agro90.70.130.080.0049011
Yam (Grosse Caille) Conv80.140.25−0.02−0.028014
Short fallow10000100
Free long fallow120000.0251200
Improved short fallow30.100.010.01300
Improved long fallow60000.025600
Table A3. Economic parameter values.
Table A3. Economic parameter values.
Crop and Practice S a S p C t C i C m P C c C e
Export banana Agro040040015,1432506008680.71840
Export banana Conv0400392.917,143366.835508680.71840
Plantain Agro002705000465800452120
Plantain Conv002605500465600452160
Cabbage Agro0010002500013001200500
Cabbage Conv0085022005012001200500
Carrot Agro0055010002413002800503
Carrot Conv0050015054120028,500503
Cassava Agro002205301065033900
Cassava Conv00221.55281461033930
Cucumber Agro001800280020085018000
Cucumber Conv001300250065079018320
Eggplant Agro006003500013001900500
Eggplant Conv0030035005011001600500
Green beans Agro009003000020003000500
Green beans Conv0070025005018006000500
Lettuce Agro0017004000022001000500
Lettuce Conv00170032005020001000500
Melon Agro00600560020012003870
Melon Conv00450560075010003970
Sugar cane Agro013.2340192830601464150
Sugar cane Conv014.7640187250551362.5200
Tomato Agro00180053000160022000
Tomato Conv0011005300400130022000
Turban squash Agro0085060007010008760
Turban squash Conv0080065001309008760
Yam (yellow) Agro00150014,00010270025000
Yam (yellow) Conv0090013,200100250025000
Yam (Grosse Caille) Agro00140611,5500210031920
Yam (Grosse Caille) Conv001000781582200024000
Short fallow00000000
Free fallow00000000
Improved short fallow0015040050000
Improved long fallow0010080050000

References

  1. Blanco-Canqui, H.; Lal, R. Soil and water conservation. In Principles of Soil Conservation and Management; Springer: New York, NY, USA, 2008; pp. 1–19. [Google Scholar]
  2. Trautmann, N.M.; Porter, K.S.; Wagenet, R.J. Modern Agriculture: Its Effects on the Environment; Cornell Cooperative Extension: Plattsburgh, NY, USA, 1985. [Google Scholar]
  3. Conway, G. The Doubly Green Revolution: Food for All in the Twenty-First Century; Cornell University Press: Ithaca, NY, USA, 1998. [Google Scholar]
  4. Altieri, M.A. Ecological impacts of industrial agriculture and the possibilities for truly sustainable farming. Mon. Rev. 1998, 50, 60. [Google Scholar] [CrossRef]
  5. Kremen, C.; Iles, A.; Bacon, C. Diversified farming systems: An agroecological, systems-based alternative to modern industrial agriculture. Ecol. Soc. 2012, 17. [Google Scholar] [CrossRef]
  6. Angeon, V.; Bates, S. L’agriculture, facteur de vulnérabilité des petites économies insulaires. Région et Développement 2015, 42, 105–130. [Google Scholar]
  7. Angeon, V.; Bates, S.; Chia, E.; Diman, J.L.; Fanchone, A.; Ozier-Lafontaine, H.; Saint-Pierre, P. Détermination des contraintes de viabilité des exploitations agricoles: Application aux Antilles françaises. In Proceedings of the 50e Colloque de L’association de Science Regionale de Langue Française, Mons, Belgium, 8–11 July 2013. [Google Scholar]
  8. Aubin, J.-P. A survey of viability theory. SIAM J. Control Optim. 1990, 28, 749–788. [Google Scholar] [CrossRef]
  9. Oubraham, A.; Zaccour, G. A survey of applications of viability theory to the sustainable exploitation of renewable resources. Ecol. Econ. 2018, 145, 346–367. [Google Scholar] [CrossRef]
  10. Tichit, M.; Hubert, B.; Doyen, L.; Genin, D. A viability model to assess the sustainability of mixed herds under climatic uncertainty. Anim. Res. 2004, 53, 405–417. [Google Scholar] [CrossRef]
  11. Sabatier, R.; Doyen, L.; Tichit, M. Modelling trade-offs between livestock grazing and wader conservation in a grassland agroecosystem. Ecol. Model. 2010, 221, 1292–1300. [Google Scholar] [CrossRef]
  12. Sabatier, R.; Doyen, L.; Tichit, M. Action versus result-oriented schemes in a grassland agroecosystem: A dynamic modelling approach. PLoS ONE 2012, 7, e33257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Mouysset, L.; Doyen, L.; Jiguet, F. From population viability analysis to coviability of farmland biodiversity and agriculture. Conserv. Biol. 2014, 28, 187–201. [Google Scholar] [CrossRef] [PubMed]
  14. Tichit, M.; Doyen, L.; Lemel, J.-Y.; Renault, O.; Durant, D. A co-viability model of grazing and bird community management in farmland. Ecol. Model. 2007, 206, 277–293. [Google Scholar] [CrossRef]
  15. Baumgärtner, S.; Quaas, M.F. Ecological-economic viability as a criterion of strong sustainability under uncertainty. Ecol. Econ. 2009, 68, 2008–2020. [Google Scholar] [CrossRef] [Green Version]
  16. Martin, S.; Deffuant, G.; Calabrese, J.M. Defining resilience mathematically: From attractors to viability. In Viability and Resilience of Complex Systems; Springer: Berlin/Heidelberg, Geramny, 2011; pp. 15–36. [Google Scholar]
  17. Sabatier, R.; Oates, L.G.; Jackson, R.D. Management flexibility of a grassland agroecosystem: A modeling approach based on viability theory. Agric. Syst. 2015, 139, 76–81. [Google Scholar] [CrossRef]
  18. Durand, M.; Désilles, A.; Saint-Pierre, P.; Angeon, V.; Ozier-Lafontaine, H. Agroecological transition: A viability model to assess soil restoration. Natl. Resour. Model. 2017, 30, e12134. [Google Scholar] [CrossRef] [Green Version]
  19. Velásquez, E.; Lavelle, P.; Andrade, M. GISQ, a multifunctional indicator of soil quality. Soil Biol. Biochem. 2007, 39, 3066–3080. [Google Scholar] [CrossRef]
  20. Camacho, N.R.; Elena, V.; Pando, A.; Thibaud, D.; Florence, D.; Patrick, L. Indicateurs synthétiques de la qualité du sol. Etude et Gestion des Sols 2009, 16, 323–338. [Google Scholar]
  21. Wezel, A.; Casagrande, M.; Celette, F.; Van, J.F.; Ferrer, A.; Peigné, J. Agroecological practices for sustainable agriculture. A review. Agron. Sustain. Dev. 2014, 34. [Google Scholar] [CrossRef] [Green Version]
  22. Garnier, E.; Desarthe, J.; Moncoulon, D. The historic reality of the cyclonic variability in French Antilles, 1635–2007. Clim. Past Discuss. 2015, 11, 1519–1550. [Google Scholar] [CrossRef]
  23. Strobl, E. Impact of hurricane strikes on local cropland productivity: Evidence from the Caribbean. Natl. Hazards Rev. 2012, 13, 132–138. [Google Scholar] [CrossRef] [Green Version]
  24. Spencer, N.; Polachek, S. Hurricane watch: Battening down the effects of the storm on local crop production. Ecol. Econ. 2015, 120, 234–240. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Global assessment of human-induced soil degradation. Source: https://databasin.org/datasets/7254137cabb042298cae0b769cba589f.
Figure 1. Global assessment of human-induced soil degradation. Source: https://databasin.org/datasets/7254137cabb042298cae0b769cba589f.
Sustainability 12 05880 g001
Figure 2. Capture basins by confidence level ( T = 40 , I * = 0.8 ).
Figure 2. Capture basins by confidence level ( T = 40 , I * = 0.8 ).
Sustainability 12 05880 g002
Figure 3. Capture basins for different hurricane impacts on the General Indicator of Soil Quality (GISQ) ( T = 40 , I * = 0.8 ).
Figure 3. Capture basins for different hurricane impacts on the General Indicator of Soil Quality (GISQ) ( T = 40 , I * = 0.8 ).
Sustainability 12 05880 g003
Figure 4. Capture basins for different time horizons ( I * = 0.8 ).
Figure 4. Capture basins for different time horizons ( I * = 0.8 ).
Sustainability 12 05880 g004
Figure 5. Capture basins for different time horizons ( I * = 0.8 ).
Figure 5. Capture basins for different time horizons ( I * = 0.8 ).
Sustainability 12 05880 g005
Figure 6. Total restoration cost for different time horizons.
Figure 6. Total restoration cost for different time horizons.
Sustainability 12 05880 g006
Figure 7. Examples of GISQ and treasury evolution over restoration period ( T = 40 , I = 0.6 , I * = 0.8 ).
Figure 7. Examples of GISQ and treasury evolution over restoration period ( T = 40 , I = 0.6 , I * = 0.8 ).
Sustainability 12 05880 g007
Figure 8. Capture basins of C4 with different values for fixed cost ( T = 40 , I * = 0.8 ).
Figure 8. Capture basins of C4 with different values for fixed cost ( T = 40 , I * = 0.8 ).
Sustainability 12 05880 g008
Figure 9. Minimum budget for improvement of 0.2 ( T = 20 ).
Figure 9. Minimum budget for improvement of 0.2 ( T = 20 ).
Sustainability 12 05880 g009
Figure 10. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 0.6 , I * = 0.8 , p = 0.8 ).
Figure 10. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 0.6 , I * = 0.8 , p = 0.8 ).
Sustainability 12 05880 g010
Figure 11. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 0.1 , I * = 0.8 , p = 0.8 ).
Figure 11. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 0.1 , I * = 0.8 , p = 0.8 ).
Sustainability 12 05880 g011
Figure 12. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 1 , I * = 0.8 , p = 0.8 .)
Figure 12. Evolution of the system for a particular realization of the uncertainty (Case C4, T = 40 , I 0 = 1 , I * = 0.8 , p = 0.8 .)
Sustainability 12 05880 g012
Table 1. Saffir–Simpson Hurricane Wind Scale (SSHWS).
Table 1. Saffir–Simpson Hurricane Wind Scale (SSHWS).
CategoryWindsDamage
174–95 mphMinimal
119–153 Km/h
296–110 mphModerate
154–177 Km/h
3 111–129 mphIntense
Major178–208 Km/h
4 130–156 mphExtreme
Major209–251 Km/h
5 ≥157 mphCatastrophic
Major≥252 Km/h
Table 2. Yearly probability of a hurricane strike in the French West Indies, by category.
Table 2. Yearly probability of a hurricane strike in the French West Indies, by category.
SSHWS CategoryYearly Probability
50.02
40.05
30.08
20.12
10.16
Table 3. Damage caused by hurricanes to crops. (Based on [24] for damages and on [22] for the probabilities).
Table 3. Damage caused by hurricanes to crops. (Based on [24] for damages and on [22] for the probabilities).
Climatic EventYearly ProbabilityAbove GroundBelow GroundWater Damaged
Minor0.289.2 %0%60.6%
Categories 1–2
Major0.15100 %100%100%
Categories 3–5
Table 4. Description of the different cases.
Table 4. Description of the different cases.
CaseCrops
Case 1Plantain, Export banana
Banana
Case 2Plantain, Export banana, Sugar cane
Banana & sugar cane
Case 3Yam(Yellow), Yam(Grosse Caille), Tomato
Tomato and yam
Case 4Yam (Yellow), Yam(Grosse Caille), Tomato, Eggplant, Lettuce, Carrot, Green bean, Cabbage, Cassava, Melon, Cucumber, Turban squash
Multicrop
Case 5Plantain, Export banana, Sugar cane, Yam (Yellow), Yam(Grosse Caille), Tomato, Eggplant, Lettuce, Carrot, Green bean, Cabbage, Cassava, Melon, Cucumber, Turban squash
All crops

Share and Cite

MDPI and ACS Style

Oubraham, A.; Saint-Pierre, P.; Zaccour, G. Viability of Agroecological Systems under Climatic Uncertainty. Sustainability 2020, 12, 5880. https://doi.org/10.3390/su12155880

AMA Style

Oubraham A, Saint-Pierre P, Zaccour G. Viability of Agroecological Systems under Climatic Uncertainty. Sustainability. 2020; 12(15):5880. https://doi.org/10.3390/su12155880

Chicago/Turabian Style

Oubraham, Aïchouche, Patrick Saint-Pierre, and Georges Zaccour. 2020. "Viability of Agroecological Systems under Climatic Uncertainty" Sustainability 12, no. 15: 5880. https://doi.org/10.3390/su12155880

APA Style

Oubraham, A., Saint-Pierre, P., & Zaccour, G. (2020). Viability of Agroecological Systems under Climatic Uncertainty. Sustainability, 12(15), 5880. https://doi.org/10.3390/su12155880

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