Original Research ARTICLE
OliveCan: A Process-Based Model of Development, Growth and Yield of Olive Orchards
- 1Departamento de Agronomía, ETSIAM, Universidad de Córdoba, Córdoba, Spain
- 2Centre for Crop Systems Analysis, Wageningen University & Research, Wageningen, Netherlands
- 3Department of Agronomy, Institute for Sustainable Agriculture, Spanish National Research Council (CSIC), Córdoba, Spain
- 4Departamento de Ciencias e Engenharia de Biossistemas, Instituto Superior de Agronomia, Universidade de Lisboa, Lisbon, Portugal
Several simulation models of the olive crop have been formulated so far, but none of them is capable of analyzing the impact of environmental conditions and management practices on water relations, growth and productivity under both well-irrigated and water-limiting irrigation strategies. This paper presents and tests OliveCan, a process-oriented model conceived for those purposes. In short, OliveCan is composed of three main model components simulating the principal elements of the water and carbon balances of olive orchards and the impacts of some management operations. To assess its predictive power, OliveCan was tested against independent data collected in two 3-year field experiments conducted in Córdoba, Spain, each of them applying different irrigation treatments. An acceptable level of agreement was found between measured and simulated values of seasonal evapotranspiration (ET, range 393 to 1016 mm year-1; RMSE of 89 mm year-1), daily transpiration (Ep, range 0.14–3.63 mm d-1; RMSE of 0.32 mm d-1) and oil yield (Yoil, range 13–357 g m-2; RMSE of 63 g m-2). Finally, knowledge gaps identified during the formulation of the model and further testing needs are discussed, highlighting that there is additional room for improving its robustness. It is concluded that OliveCan has a strong potential as a simulation platform for a variety of research applications.
Olive orchards represent the main component of agricultural systems in many semiarid regions with Mediterranean climate, reaching 10.1 Mha worldwide in 2011 (FAOSTAT, 2014). In countries where the cultivation of this tree species is done in extensive areas, olive cropping systems have become of high relevance not only from an economic perspective, but also from an ecological one. Olive orchards have been traditionally cultivated at low planting densities under low-input rainfed conditions. However, the increase in the demand for oil of recognized and consistently high quality in recent years has triggered the development and adoption of farming techniques aimed to improve productivity, such as localized irrigation, fertigation and mechanical pruning and harvesting. As a result, traditional rainfed olive orchards (<200 trees ha-1) coexist nowadays with new intensive (250–850 trees ha-1) or super-intensive (1200–3000 trees ha-1) irrigated plantations. The rapid changes in olive farming have raised questions on the economic and environmental sustainability of the different olive cropping systems under present and future climate scenarios. Given that an olive orchard is a complex system, its quantitative study via modeling is a crucial step in understanding its behavior in response to climatic and management factors.
To our knowledge, Abdel-Razik (1989) was the first researcher to describe a model for estimating the productivity of olive orchards. The model describes the growth of different organs by simulating radiation interception, photosynthesis, respiration, and applying simple allocation rules, but it does not consider the effects of planting density, canopy structure or pruning, and many of its equations lack a consistent theoretical basis. Villalobos et al. (2006) proposed a simpler approach to estimate biomass production and yield in olive canopies, based on the concept of annual radiation use efficiency and partitioning coefficients, yet this approach does not give insight about the dynamics of the system, its response to climatic variables (besides solar radiation) or the effect of management. More recently, Morales et al. (2016) presented a mechanistic model of olive oil production in the absence of any biotic or abiotic stress, based on a three-dimensional model of canopy photosynthesis and respiration and dynamic distribution of assimilates among organs. However, water stress is the main limiting factor for biomass production in rainfed and deficit-irrigated olive orchards (Moriana et al., 2003; Iniesta et al., 2009).
Simulating the water balance of an irrigated olive orchard is a particularly challenging task as the trees are typically watered by point-source emitters that keep a small fraction of the surface frequently wet while the remaining area remains dry, unless it rains. This fact results in differences between these two soil areas in relation to soil water content, the water fluxes determining the water balance (i.e., runoff, drainage, redistribution along the soil profile, soil evaporation, and root water uptake) and root length density (Fernández et al., 1991). Therefore, traditional modeling approaches based on the use of the average soil water content can lead to large errors, besides giving a poor insight into the system. One alternative consists of using a two-compartment model that solves the water balance separately for each zone of the soil. In this regard, Testi et al. (2006) proposed a model capable of simulating potential transpiration, separately calculating runoff, drainage and soil evaporation from the wet and dry fractions of the soil surface under localized irrigation. The model was developed to determine the potential irrigation needs of olive orchards, so its use is unfortunately limited to unstressed conditions. Lately, García-Tejera et al. (2017a) have formulated a soil-plant-atmosphere-continuum (SPAC) model capable of calculating root water uptake from soils with spatially heterogeneous distributions of water content and root length densities. Such a model also discretizes the soil into different soil zones and layers and, for the canopy, it considers two leaf classes (i.e., sunlit and shaded). Furthermore, the model by García-Tejera et al. (2017a) provides estimates of gross assimilation (A), offering an opportunity to link the water and carbon balances of olive trees.
The goal of this study is to present and test a process-oriented model integrating existing knowledge on the growth and development processes of olive orchards and capable to account for the impacts of water stress, management and climate on their productivity, in the absence of nutrient deficiencies, diseases and pests. The model, hereafter named ‘OliveCan’ -which comes from ‘Olive Canopy-,’ was formulated using the models by Testi et al. (2006); Morales et al. (2016) and García-Tejera et al. (2017a) as starting point.
Materials and Methods
This section provides an overview of the main features and processes within OliveCan. An in-depth description of the model, along with its equations and scientific rationale is given as Supplementary Material. The code of OliveCan was written in Visual Basic 6.0.
OliveCan is subdivided into three main components (Supplementary Figure S1) that are devoted to the computation of the water and carbon balances of the olive orchard and to simulate the impacts of some management operations. The water and carbon balance components are interdependent (i.e., each one needs data provided by the other) and both of them require information on soil traits and weather data.
Although most processes in the model run at daily time steps, others (i.e., root water uptake, photosynthesis, maintenance respiration and chilling accumulation) are computed over the diurnal course and integrated to yield daily values. The number of sub-day periods per day to be considered is customizable through a user-defined parameter (N). The meteorological input data required consist of daily values of the following variables: maximum and minimum air temperatures (Tmax and Tmin, respectively), average vapor pressure (ea), solar radiation (IG,D) average wind speed (U) and precipitation (P). For those processes simulated at sub-day intervals, OliveCan incorporates routines that disaggregate the daily values of weather data into theoretical diurnal time series.
Water Balance Component
This modeling component simulates different physical and physiological processes relevant to olive water use (Figure 1). The model solves separately the water balance for two soil zones representing the dry and wetted (by localized irrigation) surface fractions. This approach enables the model to simulate the spatial heterogeneities in soil water content dynamics associated to the use of point-source emitters. Hence, the model considers that the water supplied by irrigation (Irr) is only applied to the wetted soil zone, whereas the dry soil zone is only watered by rainfall. On the other hand, runoff (Rf), soil evaporation (Es), root water uptake (RWU), and drainage (D) represent the water effluxes for both soil zones and are computed independently for each of those. Besides, each soil zone is subdivided into a user-defined number of soil layers (n) which are also customizable in thickness. Vertical soil water fluxes between adjacent soil layers are simulated, but no lateral flow between soil zones is considered.
FIGURE 1. Relational diagram of the water balance component, highlighting the processes simulated. Solid arrows represent fluxes of water. In OliveCan, the water balance is solved separately for two soil compartments (i.e., dry and wetted soil zones), as irrigation (Irr) is applied to only one of them. For each soil zone, soil water content (i.e., 𝜃dry and 𝜃wetted) are updated every day after calculating the fluxes of effective precipitation (Peff,dry and Peff,wetted), runoff (Rfdry and Rfwetted), drainage (Ddry and Dwetted), water redistribution between the different soil layers, soil evaporation (Es,dry and Es,wetted) and root water uptake (RWUdry and RWUwetted). Direct evaporation of rain water intercepted by the canopy (Ecan) is also considered.
Daily effective precipitation (Peff) is calculated by discounting rainfall interception by the canopy (Pint) from total daily precipitation (P). Pint is calculated using a simplified version of the model of Gómez et al. (2001) and the resulting Peff is distributed proportionally between the two soil zones as a function of the surface fractions that remain rainfed or are wetted by localized irrigation. With regard to Pint, the canopy is treated as a capacitor capable of storing rain water up to a certain limit determined by canopy dimensions and leaf area index (LAI), according to Gómez et al. (2001). The stored water is subsequently lost by direct evaporation, which is simulated based on the Penman–Monteith equation assuming a null canopy resistance. As in Testi et al. (2006), the aerodynamic resistance is deduced from the model proposed by Raupach (1994), parametrized and validated specifically for olive orchards following Verhoef et al. (1997). The direct evaporation from wet foliage prevents tree transpiration (Ep), until the intercepted water is totally lost.
Runoff and infiltration are calculated following a Soil Conservation Service curve number methodology that was specifically calibrated and validated for different typologies of olive orchards (Romero et al., 2007). The approach requires information on the canopy ground cover (GC) and the soil hydrological condition (SHC) -i.e., an indicative of the capacity of infiltration of the soil when it is wet. The water content at field capacity (𝜃UL), wilting point (𝜃LL) and saturation (𝜃sat) are also needed for the computation of infiltration and all the remaining simulated processes.
Drainage and soil water redistribution processes are simulated by CERES-type sub-models (Jones and Kiniry, 1986), while soil evaporation rates (Es) are estimated using the model of Bonachela et al. (2001). The latter, calculates Es with a modified Penman-FAO equation for stage-one evaporation and uses the model of Ritchie (1972) for the soil-limited evaporation stage. For the wetted soil zone, microadvection effects are considered (Bonachela et al., 2001).
The model by García-Tejera et al. (2017a) is used to compute root water uptake (RWU) from each layer in the two soil zones, canopy transpiration (Ep) and gross assimilation (A′). By analogy with the Ohm’s law for electric circuits, the model assumes that water transport through the SPAC is driven by differences in water potential and hydraulic resistances. In this regard, three hydraulic resistances are considered: from the soil to the root-soil-interface (Rs), from the soil-root interface to the root xylem (Rr) and from the root xylem to the canopy (Rx). Rs depends on soil texture, root length density (Lv), soil water content (𝜃) (Gardner, 1960). Rr is a function of Lv and root permeability, the latter being mediated by 𝜃 (Bristow et al., 1984) and temperature (García-Tejera et al., 2016). Finally, Rx is calculated from xylem anatomical traits and tree height. In the canopy, two leaf populations are considered (i.e., sunlit and shaded). For each one, gross assimilation (A′), stomatal conductance (gs), intercellular CO2 concentration (Ci) and leaf water potential (Ψl) are calculated iteratively, considering both the models by Farquhar et al. (1980) and Tuzet et al. (2003). In doing so, the environmental CO2 concentration (Ca) is explicitly taken into account for calculating both A′ and gs on the one hand. On the other, the model requires information on the intercepted photosynthetically active radiation (IPAR) as well as the sunlit and shaded fractions of the canopy. These inputs are provided by a simple geometric model of radiation interception which assumes a spheroidal shape for the crown and accounts for the shadowing from neighboring trees. Finally, Ep is estimated from the imposed evaporation equation assuming that the canopy is coupled to the atmosphere, whereas RWU is deduced in each layer of each soil zone from the corresponding water potential differences and hydraulic resistances.
Carbon Balance Component
This modeling component is aimed to simulate the growth and development of the trees and the carbon exchange of the orchard (Figure 2). First, the model calculates the daily pool of assimilates from the rate of A′ determined by the SPAC model. Based on experimental evidence, OliveCan also assumes that such pool can be increased to some extent by reserve remobilization (Bustan et al., 2011) or fruit photosynthesis (Proietti et al., 1999). The available assimilates are allocated to growth and respiration, the latter being segregated into maintenance and growth respiration (RESPM and RESPg, respectively), of the different organs. In this regard, OliveCan considers six organ types: leaves, shoots (i.e., stems of up to 3 years-old), branches (including the trunk), coarse roots (i.e., roots with secondary growth), fine roots (i.e., absorbing roots) and fruits.
FIGURE 2. Relational diagram of the carbon balance component, highlighting the processes simulated. Squares represent carbon stocks and solid arrows the fluxes of carbon through the system. Circles and thin dashed arrows indicate key factors determining the calculation of some fluxes. Pruning residues can either be incorporated into the soil or exported, at user’s choice (dotted arrows).
RESPM is calculated as a function of temperature and biomass, and it is subtracted directly from the pool of assimilates. Whenever maintenance respiration exceeds the pool of assimilates, the deficit is discounted from the reserve pool. The remaining assimilates are distributed among the different organs with partitioning rules being mediated by phenology. The loss of carbon during the synthesis of new biomass was included by calculating a production value (PV) (Penning de Vries et al., 1974) for each type of organ according to its biochemical composition.
Two phenological stages are considered for the vegetative organs: (i) a dormant stage characterized by an absence of growth that is induced by chilling accumulation during autumn and (ii) a phase of active growth that starts in late winter, by the time average temperature is above a threshold. In relation to the reproductive growth, the date of flowering is determined with the two-phase model by De Melo-Abreu et al. (2004). Fruit growth is assumed to start after a given amount of thermal time is accumulated from the date of flowering and ceases when either maturity or the harvest date is reached.
During the vegetative rest period and provided that fruits are not present, all the available assimilates after discounting maintenance respiration are allocated to a virtual pool of reserves. Such reserve pool is subsequently used for the growth of vegetative organs and fruits during the growth season. Fruit growth can either be source-limited or sink-limited. In the former case, the associated partitioning coefficient is fixed whereas in the latter, it is calculated as a function of the number of fruits (FN), which in turn is modeled as a function of the number of fruits and nodes produced in the previous year. In doing so, the model may be prone to errors in the estimates of productivity and vegetative growth for a given year when performing long runs, but such errors are to be compensated if those model outputs are averaged over biennia. With regard to the vegetative organs, fixed partitioning coefficients are adopted. Whenever fruits are present, the model considers that they become the prioritary sink of assimilates, thus the vegetative partitioning coefficients are applied after discounting the fruit demand from the daily pool of assimilates. Therefore, partitioning coefficients to vegetative organs are assumed to be independent of tree size, management factors and environmental conditions, as in the model of Morales et al. (2016). As a final remark, inspired by the CERES-type models (Jones and Kiniry, 1986), the growth of fine roots is distributed among the different layers in the two soil zones as a function of the size and water content of each soil compartment.
Senescence of leaves and fine roots are simulated using a similar approach to that in the model by Morales et al. (2016). OliveCan takes also into account the conversion of shoots into branches when they exceed 3 years-old. Besides that, the model considers some of the effects of frost events and heat stress. Frost damage is simulated by assuming that a fraction of the standing leaves is defoliated when minimum air temperature falls below a certain temperature threshold. A similar approach is used for simulating the effect of extremely high temperatures during flowering on fruit set: when maximum air temperature exceeds a given threshold, a reduction in the final FN is triggered.
Variables related to canopy characteristics such as leaf area index (LAI) or GC are updated from the estimates of biomass of leaves assuming that the crowns present an spheroidal shape with constant leaf area density (LAD) and ratio of vertical to horizontal canopy radiuses (Rzx). Similarly, the biomass of fine roots in each soil compartment is used to compute root length density (Lv) by adopting a constant specific root length (SRL).
Finally, the soil carbon balance and heterotrophic respiration (RESPH) are computed with an adaptation of the model proposed by Huang et al. (2009) and modified to take into account the effect of soil moisture on the rate of decomposition according to Verstraeten et al. (2006). Then, by considering the different computed fluxes of assimilation and respiration within the orchard, OliveCan provides estimates of the ecosystem respiration (RESPeco) and net ecosystem exchange (NEE).
Four management operations are considered in OliveCan: tillage, irrigation, harvest and pruning. In the model, tillage operations have an impact on CN whereas irrigation provides an additional water input for the wetted soil zone. Irrigation amounts and dates can either be defined explicitly by the users or implicitly calculated through a dedicated routine that, at customizable intervals, applies a fraction of the maximum ET lost since the last irrigation. Harvesting takes place on a user-defined day of the year and results in the removal of fruits. At harvest, the model provides an estimate of oil yield (Yoil) by multiplying the dry biomass of fruits and a fixed coefficient representing the ratio of oil content to dry matter. Finally, pruning is simulated by setting a customizable fraction of LAI to be removed (Fprune) and an interval between pruning operations. The model also reduces the biomasses of shoots and branches by the same fraction Fprune. The user should indicate whether pruning residues are incorporated into the soil or exported.
Apart from the weather dataset and some orchard (e.g., planting density, age, and latitude) and soil (e.g., depth, 𝜃UL, 𝜃LL) basic traits, the user is required to enter the initial values of GC and Lv to deduce the biomasses of the different organs following simple criteria (see Supplementary Material). For the computation of FN in the first season, an estimate of dry yield for the year preceding the start of the simulation is also needed. To initialize the state variables related to phenology, simulations must start at the beginning of a year and the temperature records of the preceding 3 months must be provided. Some simulation settings such as the number of years to simulate and N must also be provided. Finally, the user is to indicate the management operations to be implemented and provide values to their parameters.
When available, the values of the different parameters were taken from the literature. Supplementary Table S2 provides a complete list with the parameter values used for the simulations and the source from which they were taken. In short, the parameters of the SPAC model were taken from García-Tejera et al. (2017a,b), who, in turn, gathered most of the parameter values from different sources. Parameters related to phenology were obtained from reports by De Melo-Abreu et al. (2004) and López-Bernal et al. (2014, 2017). The studies by Mariscal et al. (2000) and Pérez-Priego et al. (2014) were used for setting the maintenance respiration and PV coefficients, respectively. Parameters related to the calculation of fruit number and yield were taken from several sources, including experimental data (see section “Number of Fruits and Alternate Bearing” in Supplementary Material). The coefficient of oil yield to dry fruit matter was taken from experimental data collected in a hedgerow cv. ‘Arbequina’ orchard (López-Bernal et al., 2015). Partitioning coefficients were based on findings by Mariscal et al. (2000); Villalobos et al. (2006) and Scariano et al. (2008). Reports from Barranco et al. (2005) and Koubouris et al. (2009) were used to parametrize the routines modeling the impacts of frost damage and heat stress, respectively. Coefficients modulating fine root growth distribution were directly taken from Jones and Kiniry (1986). Finally, parameters implied in the soil carbon balance were taken from Verstraeten et al. (2006); Huang et al. (2009) and, to a lesser extent, from other studies.
Experimental measurements conducted in two mature olive orchards located in the Alameda del Obispo Research Station, Córdoba, Spain (37.8°N, 4.8°W, 110 m) were used for assessing the reliability of OliveCan. The climate in the area is typically Mediterranean, with around 600 mm of average annual rainfall and 1390 mm and average annual ET0 of 1390 mm (Testi et al., 2004), respectively. The soil for both orchards is classified as a Typic Xerofluvent of sandy loam texture and exceeds 2 m in depth, with field capacity (𝜃UL) and permanent wilting point (𝜃LL) water contents of 0.23 m3 m-3 and 0.07 m3 m-3, respectively (Testi et al., 2004). Weather data were collected using a station placed 500 m away from the orchards. Within both orchards, irrigation experiments comprising several irrigation treatments were performed. Each irrigation treatment was simulated separately with OliveCan.
Extensive information on the orchard characteristics and dataset of Experiment I is provided by Iniesta et al. (2009). In short, the experiment was performed between 2004 and 2006 in a high density cv. ‘Arbequina’ olive orchard (tree spacing was 7 m × 3.5 m, i.e., 408 trees ha-1) planted in 1997. Irrigation was applied 5 days a week by drip, with seven emitters of 4 L h-1 per tree. A randomized complete-block design was used with three replications of 12 trees each, and the following irrigation treatments:
• Control irrigation (CON), which applied the required water to match the maximum ET, discounting rainfall. The maximum ET was estimated using the model of Orgaz et al. (2006).
• Continuous deficit irrigation (CDI), which applied 25% of the irrigation supplied to CON, distributed throughout the irrigation season.
• Regulated deficit irrigation (RDI), which applied the same seasonal water as CDI, with a midsummer (July 1st to September 10th–15th) deficit period without irrigation.
The measurements (only performed for the central trees of the replicates) used for testing the model were oil yield (Yoil, g m-2), seasonal ET and daily transpiration (Ep, mm d-1). With regard to the former, each tree was manually harvested and the fresh yield weighed in the field. Yoil was subsequently determined from sub-samples of 5 kg of fresh fruits. Cumulative ET was determined by water balance for the whole 2005 and 2006 seasons by measuring soil water content with a neutron probe (model 503, Campbell Pacific Nuclear Corp, Pacheco, CA, United States). Eight access tubes were installed between two trees per replicate, normal to tree rows. Measurements were taken at several depths (from 0.075 to 2.65 m deep). Finally, Ep was measured in 2006 with a sap-flow system device developed and assembled in the IAS-CSIC in Córdoba and described by Testi and Villalobos (2009). The system uses the Compensated Heat Pulse (CHP) method in combination with the Calibrated Average Gradient (CAG) procedure. The probes performed readings every 15 min at 4 depths in the xylem, spaced 10 mm. Six RDI, six CDI and four CON trees were instrumented with two probes per tree, at a height of 30 cm. The outputs of each probe were integrated first along the trunk radius and then around the azimuth angle. Average sap flow records for each treatment were calibrated against the estimates of Ep deduced from the difference between the measured ET and soil evaporation in a period of several weeks with no rainfall events during the summer. The model of Orgaz et al. (2006) was used to calculate soil evaporation. The calibrated sap flow data have not been published so far.
Values of GC, LAD, and Rzx required to initialize the model were taken from measurements of tree silhouettes. A record of Ydry of the year preceding simulations was also considered. Initial Lv values were taken from records measured by Moriana (2001) for the trees of Experiment II.
Extensive information on the orchard characteristics and dataset of Experiment II is provided by Moriana et al. (2003). In short, the experiment was performed between 1997 and 1999 in a high density cv. ‘Picual’ olive orchard (tree spacing was 6 m × 6 m, i.e., 278 trees ha-1) of 18 years of age. Irrigation was applied 5 days a week by drip, with four emitters of 4 L h-1 per tree. A randomized complete-block design was used with three replications of 16 trees each, and the following irrigation treatments:
• Control irrigation (CON), which applied the required water to match the maximum ET, based on the fully replenishing soil water extraction from April to October.
• Regulated deficit irrigation (RDI), which applied 75% of the water received by CON (i.e., rainfall plus irrigation) with a midsummer deficit period (15 July to 15 September) without irrigation.
• Continuous deficit irrigation (CDI), which also applied 75% of the water received by CON (i.e., rainfall plus irrigation), but for the whole irrigation season.
• Alternate year irrigation (AYI), which was rainfed in the year of low crop load (1998) and fully irrigated, as CON, during the heavy crop years (1997 and 1999).
• Rainfed (DRY), which received no irrigation during the whole experiment.
The measurements (only performed for the central trees of the replicates) used for the model were Yoil and seasonal ET. On the one hand, trees were harvested between December 15th and January 15th for the 3 years. Individual fruit weight of each tree was measured and a subsample of 150 fruits from each tree was used for determining oil content. On the other, cumulative ET was determined by water balance for each season by measuring soil water content with a neutron probe (model 503, Campbell Pacific Nuclear Corp, Pacheco, CA, United States). Eight access tubes were installed between two trees per replicate in the four irrigation treatments and six tubes were placed in the rainfed treatment. Measurements were taken were performed at several depths (from 0.075 to 2.4 m deep).
Values of GC, LAD, and Rzx required to initialize the model were taken from dedicated measurements. A record of Ydry of the year preceding simulations was also considered. Initial Lv values were taken from records measured by Moriana (2001).
Model performances in reproducing measured data were assessed using mean absolute error (MAE, from 0 to +∞, optimum 0), root mean square error (RMSE; from 0 to +∞, optimum 0), coefficient of residual mass (CRM, from -∞ to +∞, optimum 0) and modeling efficiency (EF, from -∞ to 1, optimum 1):
Where Mi is the ith measured variable, is the average value of all measurements, Si is the ith simulated variable and n is the number of measured values. In addition, the slope, intercept and coefficient of determination (r2) obtained by regressing the simulated and measured values were also used.
Measured and simulated values of Yoil for the two experiments are also presented in Table 1 in relation to the year and the irrigation treatment. Field measurements showed two common patterns, irrespective of the experiment. On the one hand, there were consistent differences in Yoil between treatments, with the fully irrigated treatments (i.e., CON) exhibiting higher values than the deficit or rainfed ones. On the other hand, measurements revealed “high” Yoil in the first and third experimental seasons and “low” Yoil in the intermediate one. The only exception to this rule occurred in the second biennia of Experiment II for the DRY and RDI treatments, probably as a consequence of the high level of water stress reached in 1999 (Moriana et al., 2003). OliveCan reproduced well both the differences between treatments and the general alternating trend in Yoil for both experiments. Pooling all the data together, MAE and CRM were close to zero, RMSE was 63.2 g m-2 and EF was 0.48. The regression analysis yielded an r2 of 0.51 and a slope of 0.65. Slightly better results in terms of regression parameters, RMSE and EF were found when analyzing the data grouped in biennia (Table 2).
TABLE 1. Observed (Obs.) and simulated (Sim.) values of seasonal evapotranspiration (ET) and oil yield (Yoil) in the model tests using the experimental data from Iniesta et al. (2009) (Experiment I) and Moriana et al. (2003) (Experiment II).
TABLE 2. Performances of OliveCan in reproducing the experimental data gathered for the model testing.
A good agreement was found between model estimates of seasonal ET and the measurements by Moriana et al. (2003) and Iniesta et al. (2009) (Table 1). Simulated values were close to experimental data and tracked reasonably well the differences between treatments. MAE and RMSE were 36.04 mm and 89.4 mm, respectively. The statistical analysis revealed no systematic bias, with CRM close to 0, and values of EF of 0.62 (Table 2). Nevertheless, the model was found to underestimate seasonal ET for the treatments receiving more irrigation, which led to a rather low slope for the linear regression fit (0.59) between simulated and observed values.
Beyond the seasonal scale, OliveCan also provided good estimates of daily Ep, as evidenced by comparing the model outputs with the calibrated sap flow measurements recorded in Experiment I (Figure 3). The model was found to slightly underestimate Ep in winter and late autumn, particularly for the CON treatment, which led to a CRM of 0.14. In addition, Ep was slightly overestimated during the RDI midsummer non-irrigated period. However, simulated and observed values were generally very close in the three treatments. Considering the three treatments together (n = 1095), the comparison between simulated and observed Ep values yielded a MAE of 0.17 mm d-1, a RMSE of 0.32 mm d-1, a EF of 0.84 and a satisfactory linear regression fit (Table 2).
FIGURE 3. Time course of simulated and observed daily transpiration (Ep) throughout the 2006 season for the three irrigation treatments of Experiment I: the fully irrigated control (CON) (A), the continuous deficit irrigation (CDI) (B), and the regulated deficit irrigation (RDI) (C) treatments.
Figure 4 illustrates the time course of measured and modeled Ep for the CON and CDI treatments over a typical sunny summer day of Experiment I (July 21st 2006, DOY 202). The four plotted Ep curves follow a bell shape, the CON ones exhibiting higher values than those of CDI. Apart from that, Figure 4 shows a remarkable lag between the observed and the simulated daily course of Ep, regardless of the treatment, with the diurnal trends of Ep being anticipated by the model. Hence, OliveCan tended to overestimate Ep in the hours following sunrise and to underestimate it around sunset.
FIGURE 4. Daily course of simulated and observed transpiration (Ep) in a typical sunny summer day (July 21st 2006, DOY 202) for the control (CON) and continuous deficit irrigation (CDI) treatments of Experiment I.
Model tests generally revealed a high level of agreement between simulations and experimental measurements. Given the variety of the simulated treatments and the many assumptions that a model like OliveCan must take, we found the results satisfactory. Notwithstanding that, there were situations in which model estimates departed from observations. For example, some discrepancies were found for some of the simulations of Yoil, ET (Table 1) and Ep (Figure 3), but, considering that the general trends and differences between treatments were captured by the model, we believe that the results are highly acceptable. Some of the divergences between measured and simulated Yoil might be attributed to the fact that the approach followed by OliveCan to simulate alternate bearing is limited, as far as the physiological bases of alternate bearing are not completely understood yet (Connor, 2005; Dag et al., 2010). However, biennial comparisons (Table 2) only improved slightly the results. Apart from that, the remarkable lag between the simulated and measured diurnal courses of Ep (Figure 4) was to be expected: measurements were performed in the trunk with sap flow sensors and OliveCan does not simulate the buffering effect of the water stored in aboveground organs (Cermák et al., 2007). Also, the model assumes that stomatal conductance responds instantaneously to changes in environmental conditions, but the slow dynamics of stomatal opening and closing can cause lags in diurnal transpiration (Vialet-Chabrand et al., 2013).
Considering all the simulations together, the maximum simulated oil yield was 358 g m-2 (Table 1), which is comparable to the maximum values estimated by the model of Morales et al. (2016) and to available experimental data (Villalobos et al., 2006; Pastor et al., 2007). Simulated values of radiation use efficiency for oil production (i.e., the amount of oil produced per unit of intercepted PAR) averaged over biennia ranged between 0.17 and 0.10 g MJ-1. These estimates are within the range of variation found by Villalobos et al. (2006) across a wide range of commercial orchards in Southern Spain.
Overall, the results of all the aforementioned comparisons suggest that model performance is fairly satisfactory. However, further testing against experimental data taken from different environmental conditions and orchard characteristics seems highly desirable. This would help to provide additional evidence on the predictive power of OliveCan, or else to identify situations for which model accuracy could be improved through either better calibrations or reformulation of some routines. Apart from that, it should be noted that the reliability of OliveCan for estimating certain output parameters (e.g., NEE, RESPH) has not been tested specifically in the present study, which should also be the focus of future research efforts.
Considering its mechanistic approach, the vast quantity of simulated processes and its potential uses, OliveCan represents a momentous step forward in relation to previous olive growth simulation models. In this regard, OliveCan enables one to assess the combined effects of management operations and weather over crop performance for different olive orchard and soil typologies both under unstressed and water deficit conditions. Thus, the model shows potential for a broad range of research applications. For instance, OliveCan seems particularly suitable for assessing the performance of olive orchards under future climatic scenarios, as the model explicitly accounts for the multiple effects of reduced rainfall and increased environmental CO2 and temperature for the water and carbon balances of the orchard and the development of trees.
Obviously, the comprehensive nature and the wide range of simulated processes come at the expense of both model complexity and high input requirements. The latter is likely to be its main limitation, as far as some of the inputs (e.g., soil depth, Lv distribution) are not easy to measure in the field. In any case, it is noteworthy to emphasize that OliveCan has not been primarily conceived as a decision support system for farmers, but as a research tool.
During the development of the model, it became apparent that our current understanding of some of the physiological processes to be simulated was limited. For example, timing of vegetative bud break, dynamics of leaf senescence, fruit photosynthesis and the use of reserves are among the phenomena that have received less attention in the literature. Also, OliveCan is missing a sub-model aimed to properly simulate the dynamics of oil accumulation during the fruit growth period. Further research on these and other topics (e.g., alternate bearing) are clearly needed and might result in model improvements through either a more consistent parametrization or the formulation of better equations for simulating such processes.
Further research regarding genetic variability in model parameters is also desirable. With the exception of those related to the simulation of flowering date (De Melo-Abreu et al., 2004) and frost damage (Barranco et al., 2005), all parameters have been taken from past experiments carried out either with only one cultivar each (‘Arbequina’ being the most frequent) or averaging the results obtained for a few of them. Although the scarce literature does not allow us to disentangle how many of these crop parameters are cultivar-specific, it is clear that exploring their genetic variability might be important for enhancing model reliability. Moreover, the quantification of such cultivar variability may be used for evaluating its impact on tree physiology and productivity under different management, weather or orchard characteristics using OliveCan, which may be useful for breeding purposes.
Finally, future improvements of OliveCan might include additional sub-models for simulating nutrient uptake and the impact of pests and diseases. Apart from that, the model shows potential for being adapted to other tree species, so its interest may not be only restricted to olive researchers.
The model presented here targets the simulation of the interactions between olive trees and their environment through a detailed characterization of the water and carbon balances of the orchard as affected by weather variables, soil attributes and management operations. The generally high level of agreement found between measured and simulated data evidence the suitability of OliveCan for estimating olive orchard dynamics. These results encourage the application of the model to simulate the growth, carbon exchange and water relations of olive orchards in a wide range of research contexts, including studies on the performance of olive trees under climate change scenarios. The development of OliveCan has also highlighted significant knowledge gaps in relation to some physiological processes and the cultivar specificity of some of the parameters. Further research on these aspects may contribute to improve the reliability of the model.
All authors played a significant role in the conception and development of the model. FV led out the coding, with contributions from ÁL-B, AM, OG-T, and LT. ÁL-B, LT, and FV gathered the datasets for testing the model. ÁL-B led out the writing with significant contributions from all co-authors.
The research leading to these results has received funding from Ministerio de Economía y Competitividad (Grant Nos. AGL-2010-20766 and AGL2015-69822), from Junta de Andalucía (Grant No. P08-AGR-04202), from the European Community’s Seven Framework Programme-FP7 (KBBE.2013.1.4-09) under Grant Agreement No. 613817. 2013–2016 “MODelling vegetation response to EXTREMe Events” (MODEXTREME, modextreme.org) and from ERA-NET FACCE SURPLUS (Grant No. 652615, project OLIVE-MIRACLE), the latter co-funded by INIA (PCIN-2015-259). Besides, ÁL-B was funded by a postdoctoral fellowship (‘Juan de la Cierva-Formación 2015’ Programme, FJCI-2015-24109) from Ministerio de Economía y Competitividad.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00632/full#supplementary-material
Bristow, K. L., Campbell, G. S., and Calissendorff, C. (1984). The effects of texture on the resistance to water-movement within the rhizosphere. Soil Sci. Soc. Am. J. 48, 266–270. doi: 10.2136/sssaj1984.03615995004800020007x
Bustan, A., Avni, A., Lavee, S., Zipori, I., Yeselson, Y., Schaffer, A., et al. (2011). Role of carbohydrate reserves in yield production of intensively cultivated oil olive (Olea europaea L.) trees. Tree Physiol. 31, 519–530. doi: 10.1093/treephys/tpr036
Cermák, J., Kucera, J., Bauerle, W. L., Phillips, N. G., and Hinckley, T. M. (2007). Tree water storage and its diurnal dynamics related to sap flow and changes in stem volume in old-growth Douglas-fir trees. Tree Physiol. 27, 181–198. doi: 10.1093/treephys/27.2.181
Dag, A., Bustan, A., Avni, A., Tzipori, I., Lavee, S., and Riov, J. (2010). Timing of fruit removal affects concurrent vegetative growth and subsequent return bloom and yield in olive (Olea europaea L.). Sci. Hort. 123, 469–472. doi: 10.1016/j.scienta.2009.11.014
De Melo-Abreu, J. P., Barranco, D., Cordeiro, A. M., Tous, J., Rogado, B. M., and Villalobos, F. J. (2004). Modelling olive flowering date using chilling for dormancy release and thermal time. Agric. For. Meteorol. 125, 117–127. doi: 10.1016/j.agrformet.2004.02.009
FAOSTAT (2014). Available at: http://faostat.fao.org
Fernández, J. E., Moreno, F., Cabrera, F., Arrue, J. L., and Martin-Aranda, J. (1991). Drip irrigation, soil characteristics and the root activity of olive trees. Plant Soil 133, 239–251. doi: 10.1007/BF00009196
García-Tejera, O., López-Bernal, A., Testi, L., and Villalobos, F. J. (2017a). A soil-atmosphere continuum (SPAC) model for simulating tree transpiration with a soil multi-compartment solution. Plant Soil 412, 215–233. doi: 10.1007/s11104-016-3049-0
García-Tejera, O., López-Bernal, A., Testi, L., and Villalobos, F. J. (2017b). Analysing the interactions between wetted area and irrigation volume in an olive orchard using a SPAC model with a multi-compartment soil solution. Irrig. Sci. 35, 409–423. doi: 10.1007/s00271-017-0549-5
García-Tejera, O., López-Bernal, A., Villalobos, F. J., Orgaz, F., and Testi, L. (2016). Effect of soil temperature on root resistance: implications for different trees under Mediterranean conditions. Tree Physiol. 36, 469–478. doi: 10.1093/treephys/tpv126
Huang, Y., Yu, Y., Zhang, W., Sun, W., Liu, S., Jiang, J., et al. (2009). Agro-C: a biogeophysical model for simulating the carbon budget of agroecosystems. Agric. For. Meteorol. 95, 203–223. doi: 10.1016/j.agrformet.2008.07.013
Iniesta, F., Testi, L., Orgaz, F., and Villalobos, F. J. (2009). The effects of regulated and continuous deficit irrigation on the water use, growth and yield of olive trees. Eur. J. Agron. 30, 258–265. doi: 10.1016/j.eja.2008.12.004
Koubouris, G. C., Metzidakis, I. T., and Vasilakakis, M. D. (2009). Impact of temperature on olive (Olea europea L.) pollen performance in relation to relative humidity and genotype. Environ. Exp. Bot. 67, 209–214. doi: 10.1016/j.envexpbot.2009.06.002
López-Bernal, A., Garcia-Tejera, O., Orgaz, F., Testi, L., and Villalobos, F. J. (2014). “Olive bud dormancy is induced by low temperatures” in Proceedings of the XIIIth Congress of the European Society for Agronomy, eds P. Pepó and J. Csajbók (Debrecen: University of Debrecen), 253–254.
López-Bernal, A., García-Tejera, O., Vega, V. A., Hidalgo, J. C., Testi, L., Orgaz, F., et al. (2015). Using sap flow measurements to estimate net assimilation in olive trees under different irrigation regimes. Irrig. Sci. 33, 357–366. doi: 10.1007/s00271-015-0471-7
López-Bernal, A., Villalobos, F. J., García-Tejera, O., Testi, L., and Orgaz, F. (2017). Do olive vegetative buds undergo a real dormant state in Winter? Acta Hortic. 1160, 227–230. doi: 10.17660/ActaHortic.2017.1160.33
Mariscal, M. J., Orgaz, F., and Villalobos, F. J. (2000). Radiation-use efficiency and dry matter partitioning of a young olive (Olea europaea) orchard. Tree Physiol. 20, 65–72. doi: 10.1093/treephys/20.1.65
Morales, A., Leffelaar, P. A., Testi, L., Orgaz, F., and Villalobos, F. J. (2016). A dynamic model of potential growth of olive (Olea europaea L.). Eur. J. Agron. 74, 93–102. doi: 10.1016/j.eja.2015.12.006
Orgaz, F., Testi, L., Villalobos, F. J., and Fereres, E. (2006). Water requirements of olive orchards-II: determination of crop coefficients for irrigation scheduling. Irrig. Sci. 24, 77–84. doi: 10.1007/s00271-005-0012-x
Pastor, M., García-Vila, M., Soriano, M. A., Vega, V., and Fereres, E. (2007). Productivity of olive orchards in response to tree density. J. Hortic. Sci. Biotechnol. 82, 555–562. doi: 10.3389/fpls.2017.01280
Penning de Vries, F. W. T., Brunsting, A. H. M., and van Laar, A. H. (1974). Products requirements and efficiency of biosynthesis, a quantitative approach. J. Theor. Biol. 45, 339–377. doi: 10.1016/0022-5193(74)90119-2
Pérez-Priego, O., Testi, L., Kowalski, A. S., Villalobos, F. J., and Orgaz, F. (2014). Aboveground respiratory CO2 effluxes from olive trees (Olea europaea L.). Agrofor. Syst. 88, 245–255. doi: 10.1007/s10457-014-9672-y
Raupach, M. R. (1994). Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index. Bound. Layer Meteorol. 71, 211–216. doi: 10.1007/BF00709229
Scariano, L., Lo Bianco, R., Di Marco, L., and Policarpo, M. (2008). “Dynamics of dry matter partitioning in young ‘Nocellara del Belice’ olive trees,” in Proceedings of the Fifth International Symposium on Olive Growing, Vols 1 and 2, ed. M. T. L. S. F. L. Ozkaya (Leuven: Acta Horticulturae), 397–401. doi: 10.17660/ActaHortic.2008.791.58
Testi, L., Villalobos, F. J., Orgaz, F., and Fereres, E. (2006). Water requirements of olive orchards: I simulation of daily evapotranspiration for scenario analysis. Irrig. Sci. 24, 69–76. doi: 10.1007/s00271-005-0011-y
Verhoef, A., McNaughton, K. G., and Jacobs, A. F. G. (1997). A parameterization of momentum roughness length and displacement height for a wide range of canopy densities. Hydrol. Earth Syst. Sci. 1, 81–91. doi: 10.5194/hess-1-81-1997
Verstraeten, W., Veroustraete, F., and Feyen, J. (2006). On the temperature and water limitation of net ecosystem productivity: implementation in the C-Fix model. Ecol. Modell. 199, 4–22. doi: 10.1016/j.ecolmodel.2006.06.008
Vialet-Chabrand, S., Dreyer, E., and Brendel, O. (2013). Performance of a new dynamic model for predicting diurnal time courses of stomatal conductance at the leaf level. Plant Cell Environ. 36, 1529–1546. doi: 10.1111/pce.12086
Keywords: carbon assimilation, crop model, Olea europaea L., SPAC model, water stress, water uptake
Citation: López-Bernal Á, Morales A, García-Tejera O, Testi L, Orgaz F, De Melo-Abreu JP and Villalobos FJ (2018) OliveCan: A Process-Based Model of Development, Growth and Yield of Olive Orchards. Front. Plant Sci. 9:632. doi: 10.3389/fpls.2018.00632
Received: 04 September 2017; Accepted: 23 April 2018;
Published: 09 May 2018.
Edited by:Mariela Torres, Instituto Nacional de Tecnología Agropecuaria (INTA), Argentina
Reviewed by:Mukhtar Ahmed, Pir Mehr Ali Shah Arid Agriculture University, Pakistan
Hartmut Stützel, Leibniz University of Hanover, Germany
Peter S. Searles, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Argentina
Copyright © 2018 López-Bernal, Morales, García-Tejera, Testi, Orgaz, De Melo-Abreu and Villalobos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luca Testi, firstname.lastname@example.org