ISSN 2280-6180 (print) © Firenze University Press ISSN 2280-6172 (online) www.fupress.com/bae Bio-based and Applied Economics 2(3): 301-321, 2013 An integrated PMP model to assess the development of agro-energy crops and the effect on water requirements Michele Donati1,*, Diego BoDini2, Filippo arFini2, annalisa ZeZZa3 1 Department of Biosciences, University of Parma, Italy 2 Department of Economics, University of Parma, Italy 3 National Institute of Agricultural Economics, Italy Abstract. This paper presents an integrated model for the economic and environmen- tal assessment of the use of natural resources when new activities (i.e. biomass crops for energy production) are introduced into the farm production plan. The methodol- ogy is based on the integration of positive mathematical programming (PMP) with the AquaCrop model developed by FAO. PMP represents farmer decision processes and evaluates how farms react to the biomass-sorghum activity option at different price lev- els. AquaCrop evaluates the relationship between water needs and biomass production and assesses the effect of the land allocation on water requirements at regional level. The integration of these two models assists global policy evaluation at regional level as it makes it possible to identify the economic threshold for biomass crops, the change in land allocation and total water requirement. The model can help policy makers to evalu- ate the impacts of variations in crop profitability and market innovations on farm profit- ability, land use and water consumption and the sustainability of the market scenario. Keywords. Agroenergy economics, Positive Mathematical Programming, water demand, AquaCrop JEL codes. C61, Q42, Q24, Q25, Q51 1. Introduction Farmer interest in agro-energy production is mainly due to the high level of financial subsidies for renewable energy production. European strategy against climate change indi- cates agriculture as one of the main sectors that should contribute through specific actions to achieve the objectives established for 2020 (European Commission, 2010). The 2014- 2020 CAP reform also states that European agriculture can play a major role in mitigating global warming through new entrepreneurial vision and the substitution of fossil energy with renewable energy (European Commission, 2009). Agro-energies are thus considered as an instrument to enhance energy security, reduce greenhouse gas emissions and raise farm income (Petersen, 2008). Agricultural biomass and agronomic management are the main resources to generate environmentally virtuous processes on farms. * Corresponding author: michele.donati@unipr.it. 302 M. Donati, D. Bodini, F. Arfini, A. Zezza Despite the theoretical advantages for European society, the threats as well as the benefits of agro-energy need to be carefully evaluated. On one hand, agricultural biomass used for producing fuel, heat and electricity has favourable effects on the environment, but on the other hand, it can determine end-consumption competition and excessive pressure on natural resources. Biofuel gives rise to concern about the sustainability of the non-food supply chains. Rising demand for biofuel has led to widespread exploita- tion of arable land for growing first and second generation biofuel crops. According to the internal biofuel strategy, 10% of total transport fuel should be biofuels (biodiesel or bioethanol) in European Union (EU) by 2020 (European Commission, 2006). If this is considered in terms of available land surface, the risk of agricultural land overexploita- tion is real (Zezza, 2008). Another example of how renewable energies can contribute to a sustainable society is biogas. Biogas in the Po Valley, in Italy, is very profitable as a result of large public sub- sidies and very short payback time, and is expanding rapidly. Most of the biogas plants maximize economic performance using biomass that in part originates from waste such as slurry and agroindustrial wastes, and often from specialized crops, in particular maize. About 200 ha of dedicated crops are used in a biogas plant of 1MW. This type of agro- energy farms present a low level of crop diversification, which brings the risk of creating monoculture requiring intensive use of pesticide. They also use large quantities of water for irrigation for the entire plant growing cycle. Overexploiting resources means reducing agricultural commodities for food and feeding purposes but also jeopardizing the equilib- rium of the ecosystem. Soil depletion, biodiversity reduction, aquifer pollution and water scarcity are important concerns which need to be taken into consideration alongside glob- al warming. Paradoxically, it sometimes appears that it is necessary to sacrifice the envi- ronment in order to save it (Doornbosch and Steenblick, 2007). To produce food and fibre, EU agriculture takes about 50% of the available fresh water, while at global level this percentage is 70% (OECD, 2010). The rapid growth of agro-energy poses the need for management of irrigation water in order to avoid excessive pressure on the resource and water crisis in drought periods (OECD, 2010). Increasing demand for water from other economic sectors and society worsens the problem of water scarcity. In line with the EU Water Framework Directive (WFD) (Directive n. 2000/60/ EC), Member States will have to apply administrative and economic tools for saving irri- gated water and controlling the use of water for agriculture. The objective of the WFD is threefold: protect water resources, optimize water use and sustain agricultural productiv- ity. As a limited resource, water should be managed through appropriate economic instru- ments to drive user behaviour to improve efficiency in water distribution and application. WFD encourages the use of a set of economic instruments, including the tariff system and the market for water entitlements to rationalize water allocation, which might assign an economic value to water. Without an explicit value, users pursuing private interests will use water as a free access resource compromising its sustainable exploitation (Har- din, 1968). But if water is perceived as a scarce resource and an economic good, farmers have to compare alternative crops to identify those maximising farm income and mini- mizing the cost of the water consumption (García-Vila and Fereres, 2012; Bazzani et al., 2007; Bazzani et al., 2004; Ward and Michelsen, 2002). The problem becomes particularly important in Mediterranean regions and in drought periods (Howell, 2001). 303PMP to assess agroenergy impact on water consumption Agricultural biomass production for agro-energy supply chain is an alternative that farmers are increasingly considering in their decision making, but according to the type of crop it has a big effect on water availability. In the case of sorghum for second generation bioethanol production, water consumption for irrigation is not high, but in other cases, like maize for biogas, where maize is intensively cultivated, it can be much higher. Evaluating the impact of agro-energy production on land and water use can support policy makers in planning intervention to optimize the use of scarce resources. Local water management authorities in particular need to be made aware of likely farm allo- cation decisions and the potential water consumption of agro-energy crop development. This information supports water distribution planning and implementation of tools to regulate water for irrigation. To support policy and management decisions, quantitative models based on mathematical programming (MP) are usually implemented. The literature shows that MP can be used for farm management and policy assess- ment. MP models generate an optimization process of an objective function subject to a set of constraints and can be implemented following normative and positive approaches. The normative approach can be considered the “classical” MP tool for farm management. The characteristics of normative MP models are related to the level of knowledge about the farmer technological set, prices and costs and farmer assumptions in a suboptimal condition. The model requires a large amount of technical and economic data. A nor- mative MP model will typically identify optimal production level and the optimal use of inputs to maximize revenue or minimize costs. These models have a prescriptive charac- ter, they indicate what the decision maker ought to do in order to optimise his objective and they do not reproduce what he is actually doing. Normative models are useful for their capacity to predict the use of inputs. In the case of water, demand can be calculated for single farms or for groups of farms. For groups of farms, models can reproduce water supply nodes in order to reorganize the network and water allocation more efficiently (Harou et al., 2009; Bartolini et al., 2007; Rosengrant et al., 2000; García-Villa and Fer- eres, 2012). Positive MP models, or positive mathematical programming (PMP) (Howitt, 1995; Paris and Howitt, 1998; Heckelei et al., 2012) are used on the other hand for policy assess- ment when the size of inputs and variable costs are not precisely known. In a context of poor information, they can predict farm use under the hypothesis that the observed pro- duction level is considered optimal by the entrepreneur. Variation in output market price, or in specific coupled (or decoupled) payments, leads to a new optimal production lev- el. The main feature of positive MP models is to calibrate, for a given farm, the observed production level and to estimate a non linear cost function that reproduces the cost of the inputs. The technological matrix defines the productivity level for the observed crops and activities. PMP typically identifies the cost that economic agents are willing to pay in order to be optimal at the observed level. These models can be developed for a single farm or a group of homogeneous farms belonging to a large region (Arfini and Donati, 2011). While normative MP models can describe in great detail the technological level and predict the use of inputs, they meet difficulties assessing the impacts of new market and policy scenarios when many farms in a region are considered. This is because of the dif- ficulty in collecting and differentiating technological information among farms. By con- trast, positive MP models can easily estimate the cost of technology and thus assess the 304 M. Donati, D. Bodini, F. Arfini, A. Zezza impact of new market and policy scenarios for a large group of farms. Such farms usually belong to large samples such as FADN which do not collect micro information related to the input use of each farm activity. For PMP models, the drawback is the lack of infor- mation related to the physical use of inputs. There are two possible strategies in order to provide information on input use: i) consider in PMP models the use of specific inputs (Helming and Peerlings, 2003), although there is the problem of considering different technological sets and input-output relationships for many farms; ii) integrate PMP mod- els with other methodologies that consider more specifically input-output relationships for some specific inputs. Water as an input can be introduced into MP models (including PMP) in vari- ous ways: i) water is considered as a constraint with fixed water requirement coefficients (Graveline and Mérel, 2012; Cortignani and Severini, 2009; Medellin-Azuara et al., 2010); ii) water is included within the production function linking crop yields to water applica- tion (García-Villa and Fereres, 2012; Graveline and Mérel, 2012) through the introduction of biophysical information including local environment characteristics (García-Villa and Fereres, 2012; Cortignani and Severini, 2009). In particular, a model estimating the rela- tionship between yields and water use has recently been developed by FAO: AquaCrop. AquaCrop (Steduto et al., 2009) permits evaluation and simulation of yield responses with regard to water application for a group of crops. It takes account of specific information on climate conditions, soil characteristics and irrigation management. PMP has shown to be very efficient for policy analysis purposes (Heckelei et al., 2012) and particularly for assessing the introduction of a new crop in the production pattern of a group of farms (Arfini and Donati, 2013). The present analysis was conducted consider- ing new activities as “latent” information that the entrepreneur can use for maximizing the objective function. PMP provides results using small amounts of information on the latent crop (yield and market price) without a detailed description of the technological set for all the inputs. It was used to evaluate the farmer’s ability to choose agricultural activi- ties not observed in the base year (Röhm and Dabbert, 2003; Blanco et al., 2008; Arfini and Donati, 2013) This paper presents an integrated framework of models answering research questions on assessment of impact of a new activity (e.g. agro-energy crops) on the agricultural pro- duction plan of a region. The assessment considers key issues for policy makers: land use, supply variations, territorial specialization, economic impact and environmental implica- tions of the change in water use by farmers. The framework integrates micro-based PMP and the AquaCrop model. More specifically the goal is to assess the effects of the intro- duction of sorghum for biomass production on land and water allocation in the province of Parma, in Emilia Romagna region in northern Italy. The model simulates different mar- ket scenarios for sorghum in order to evaluate the level of sustainability and policy impli- cations at different market prices. The paper consists of five sections. The first discusses PMP methodology using latent crop information. The second presents the AquaCrop model, adapted to the character- istics of the area under investigation and the purpose of this research. The third section describes the integration of the models and their structure. The fourth presents results of simulations using AquaCrop. The fifth section presents the main results and their implica- tions for policy. In the conclusion, opportunities for future research are discussed. 305PMP to assess agroenergy impact on water consumption 2. PMP model with latent information PMP appeared in the world of MP quite recently thanks to the pioneering studies of Howitt (1995) and Paris and Arfini (1995). Its impact on agricultural economists was important in generating a wide field of literature that has led over time to more sophisti- cated modelling. Various elements have influenced the development of PMP models over time: research objectives, characteristics of available information, number of farms in the sample, number of farms that are represented by the models, level of representativeness of the data at regional level, method of calibration, method of estimating the non linear cost function as well as theoretical assumptions underpinning the models. Two distinct strands can be seen in the literature: the first considers PMP as a calibration method (Howitt, 1995; Heckelei and Wolff, 2003), while the second considers PMP as a method of esti- mating variable costs (Arfini and Paris, 1995; Paris and Howitt, 1998; Paris and Arfini, 2000; Paris, 2012). This paper follows the approach proposed by Paris (2012), where PMP assesses the impacts of market scenarios and agricultural policy by estimating the variable cost function associated with the use of inputs. As regards the PMP models that consider latent information, the background hypoth- esis is based on the assumption that farmers have knowledge of a set of information regarding production activities, which is larger than the set of information that can be observed from the production plan. Some of this farmer information comes from their previous experience, some from the experience of their neighbours or from advice on new crops given by experts (agronomists). The economic and technological information regarding activities that are perceived as more costly or more risky will not be used in the production plan, but will remain latent until it becomes economically useful. The decision to change the production plan or technological choices is usually moti- vated by variables such as risk aversion, level of technical knowledge, availability of capi- tal, family structure, age of the farmer, the presence of support agencies in the territory, etc. All these variables affect the farmer’s decision process and lead to the selection of a certain combination of crops. Why does a farmer produce soft wheat and alfalfa and not, for instance, tomato and sugarbeet, which are produced on other farms in the area? The farmer could potentially insert tomato and sugar beet into his production plan, but he does not currently produce them because they are not profitable for his farm. Until they become profitable, information is “latent” in the sense that it is known by the farmer but not used. Economic information related to farm crops identified as “latent” can however be introduced into a PMP model in two different ways: 1. as “latent technology”, when a given crop is adopted only by a group of farms in the sample: each farm belonging to the same sample is cross-linked with the others through a common set of cost function parameters and shares the same technology. Because there is self-selection, a given crop existing in the production pattern of farm “A” can be adopted by farm “B” if there is economic advantage; 2. as “latent crop”, when the information is related to a given crop that does not exist in the farm production plan for any farms belonging to the sample. Of course, in the real world farmers are not alone. They are in an environment char- acterized by different production decisions deriving from different production possibili- 306 M. Donati, D. Bodini, F. Arfini, A. Zezza ties, among which the farmer selects what he assumes to be the best solution. The main driving force that leads farmers to select one production plan and not others is the cost function associated with each activity. The total cost function is the economic measure of the available technology and all the other factors. The total cost includes all variable costs perceived by the farmer for each activity. Some of these costs are clearly registered by bookkeeping and identified as explicit costs, while others are only perceived by entrepre- neur and are implicit costs of the decision. Farmer behaviour can be formally reproduced by a production plan where the observed activities are accompanied by latent activities. We present the PMP model, articulated in three phases (Paris and Howitt, 1998), includ- ing non-observed activities, and consider, for sake of simplicity, only one farm with two sets of products: activities realized r (for r=1,2,…,R) and latent activities l (for l=1,2,…,L). We assume that xr is the vector of the realized output quantities and x l the vector of latent output quantities. The observed quantity levels for realized and latent activity are known and correspond to xr and x l respectively. Farm activity is subject to limiting factors i (for i=1,2,…,I) and the upper bound values for each factor are included in vector b farm tech- nology is provided by the coefficients of the I by R matrix Ar for the realized activities, and of the I by L matrix A l for the latent activities. Given this information, we can devel- op the first PMP phase through a problem that maximizes the gross margin (GM) as fol- lows: = − + − ≥ ≥ GM p c x p c xmax ( ) ' ( ) 'r r r l l lx x0, 0r l (1) where pr and pl represent the vectors of realized output prices and latent output prices respectively, while, cr and cl the vectors of the exogenous specific costs for the realized and latent activities. The objective function (1) is maximized with respect to the non-negative variables xr and x l and is subject to the following constraints: + ≤A x A x b y( )r r l l (2) ε λ≤ +x x ( )r r r (3) ε λ≤ +x x ( )l l l (4) Constraint (2) identifies the relationship between the total demand of input to pro- duce xr and x l (left hand side) and the total input supply (right hand side). The shadow prices of the binding farm resources b are represented by the vector y. Constraints (3) and (4) are the PMP calibrating constraints, while λr λl are the dual values of the realized and latent activities. λr represents the hidden costs, i.e. the implicit marginal costs of the realized activities, and λl the implicit marginal cost associated with the latent activities. 307PMP to assess agroenergy impact on water consumption In the definition of the problem (1)-(4), the latent crop is added from the first phase, assigning it a very low production level x l close to zero, while the data related to the pric- es and specific costs are taken from market and must guarantee a condition of positive marginal profit. Yields are assumed by experts and by literature. This first PMP setting provides dual information for realized and latent crops to be used in the second phase, where a non-linear function is estimated. We choose the follow- ing quadratic cost function             x x Q x x 1 2 r l rl r l (5) where matrix Q is symmetric positive semidefinite and includes parameters to be estimat- ed by properly methods. In this work, the parameter estimation is carried out adopting the maximum entropy approach (Paris and Howitt, 1998) considering the following rela- tionship: λ λ         +         =         c c x x Qr l r l r l rl (6) where the explicit and implicit costs recovered in the previous phase should be equal to the marginal cost derived from the quadratic total cost function (5). In analytical terms: λ λ λ λ λ λ                             +                             =                                                         � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � � c c c c c c x x x x x x q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q r r rR l l lL r r rR l l lL r r rR l l lL r r r r r rR r l r l r lL r r r r r rR r l r l r lL rRr rRr rRrR r l rRl rRlL l r l r l rR l l l l l lL l r l r l rR l l l l l lL lLr lLr lLrR lLl lLl lLlL 1 2 1 2 1 2 1 2 1 2 1 2 1 1 1 2 1 1 1 1 2 1 2 1 2 2 2 2 1 2 2 2 1 2 1 1 2 1 1 1 2 1 1 1 1 2 1 2 1 2 2 2 2 1 2 2 2 1 2 1 2 (7) As shown in Equation (7), the parameters of the Q matrix provide the information about the substitution and complementarity relationships among activities and, thus, between realized and latent activities (Arfini and Donati, 2013; Paris and Howitt, 1998). The new non linear cost function estimated by using the maximum entropy technique is used in the third phase of PMP to calibrate the observed situation without the calibrat- ing constraints: 308 M. Donati, D. Bodini, F. Arfini, A. Zezza + −             + ≤ ≥ ≥ S t p x p x x x Q x x A x A x b max 1 2 . . r r l l r l rl l l r r l l x x ' 0, 0 ' r l (8) At this stage, all the information about latent activities is incorporated in a model that can be applied to evaluate policy and market scenarios. The model output shows the change in resource allocation (e.g. land), the marginal value of resources (dual values) and the dynamics in output levels as well as other important economic information on rev- enue, subsidy, total variable cost and gross margin variations. 3. An overview of AquaCrop In an economic system characterized by a limited availability of resources, models that can indicate a better allocation of inputs and resources are a key to making optimal decisions. Evaluating the use of sensitive inputs for the environment, like water, becomes more complex if carried out at regional level. Specific environmental characteristic of an area (e.g. rainfall, soil characteristics, temperature, etc.) have to be associated with the technical capability of farmers and with crop profitability. Moreover, yields are the sim- plest expression of crop productivity, but a direct relation between yield and input use is hard to identify especially if the information is not directly registered in a bookkeeping system. Researchers have thus developed models able to create a link between yields and input use on the basis of information collected over time or by experimental methods. One of the models which can best simulate water-limited attainable yield is AquaCrop developed by FAO (Steduto et al., 2009). The model provides and predicts the crop yield according to water needs and different irrigation methods. FAO methodology consid- ers the use of empirical production functions to evaluate crop yield response to water (Doorenbos and Kassam, 1979). The central feature is the following equation, relating yield to water used in the irrigation process: −    = −    Y Y Y k ET ET ET X a X y X a X (9) YX Ya are the maximum and actual yields, ETX and ETa are the maximum and actual evapotranspiration, and ky the proportionality factor between relative yield loss and rela- tive reduction in evapotranspiration. The problem of water scarcity led to a different theorization of equation (9), which now separates field crops from tree crops. In particular, for field crops, FAO research- ers proposed re-elaborating the equation in order to plan, manage and simulate different water management scenarios. So on the basis of work by Doorenbos and Kassam (1979), AquaCrop evolved by: i) dividing the ET into crop transpiration Tr and soil evaporation (E); ii) developing a model relating to canopy growth and senescence in order to estimate 309PMP to assess agroenergy impact on water consumption Tr and its separation from E; iii) considering the final yield Y as a function of final bio- mass B and harvest index HI and iv) considering and separating the effects of water stress into four components: canopy growth, canopy senescence, Tr and HI. The central equation of the new AquaCrop is thus: ∑= ×B WP Tr (10) where Tr is the crop transpiration (in mm) and WP is the water productivity parameter (kg of biomass per m2 and per mm of cumulated water transpired over the time period in which the biomass is produced). Equation (10) implies a shift from the seasonal or long- term evaluation used by Doorenbos and Kassam (1979) to a daily time relationship that is closer to the time scale of crop responses to water deficits. AquaCrop considers three important aspects: soil, crop and atmosphere. Soil is con- sidered by determining the level of fertility that can affect crop development as well as water balance. Crops and plants are measured and simulated using data relating to growth, development, and yield processes. Atmosphere is reproduced considering thermal regime, rainfall, evaporative demand, and carbon dioxide concentration. In order to better simulate reality, a wide range of parameters such as the irrigation system, water productiv- ity, crop adjustments to stress can be modified and reflect on the final yield. 4. Model architecture The model is divided into three parts or modules. Each module is devoted to a specif- ic task and interfaced with the others providing the input information. Figure 1 presents a scheme where the different phases of the model are specified and linked. The first module provides information about the characteristics of the farms belong- ing to the sample and information related to the water use. In a diversified territorial con- text where different agronomic areas are embedded, it becomes important to cover differ- ent farm types and to truthfully represent land use and productivity levels. There is thus a compromise between minimizing the amount of data entered into the model and the need to provide enough details of individual farmer behavior and technologies and pro- duction decisions at farm level. The use of two databases providing complementary farm information is very useful: IACS (Integrated Administration and Control System) database provides information on land use of each farm, while FADN (Farm Accountancy Data Network) database provides economic2 and technical information regarding farm type. The integration of IACS with FADN makes it possible to measure the exact dimension of agricultural production systems, gross marketable output, subsidies distributed and the amount of variable costs attributable to each process in areas smaller than NUTS3 (Arfini et al., 2005). The integration of FADN and IACS databases with AquaCrop into a single database gives a complete dataset of land use, technical and economic parameters for production processes and water use of all the farms included in the area considered by the model. The 2 FADN Farm Accountancy Data Network, provides information for Italy at farm level on the explicit variable cost (accounting cost) per activity. 310 M. Donati, D. Bodini, F. Arfini, A. Zezza aggregation is performed at macro-farm level, i.e. farms grouped by size and farm special- ization and by each agricultural area (a homogeneous altitude area belonging to the same province). More precisely, in each province three altitude levels, seven classes of size (0-10 ha, 10-20 ha, 20-30 ha, 30-50 ha, 50-100 ha, 100-300 ha, > 300 ha) and three economic sectors (arable crops, fruit and vegetables, and animal production) were considered. Figure 1. PMP model architecture. AquaCrop simulations IACS dataFADN data PMP Model Constraints Land constraints Agronomic constraints Land allocation Farm income dynamics Water allocation Input PMP model Output Water constraints CAP constraints Scenarios CAP and market prices Agroenergy crops The second module consists of the PMP optimization model, which estimates variable costs for all activities and assesses the impact of policy and market scenarios. There are two calibration phases: the first is obtained by n linear programming models (one for each macro-farm) adopting the calibration constraints, and the second calibration is achieved by using the non-linear cost function estimated in the second PMP phase. The second calibration allows to verify that the model reproduces the observed land allocation with- out calibrating constraints. Once verified, the model is ready to be implemented for the simulation phase of all sets of constraints: land, agronomic, water and policy constraints (CAP). The simulation phase considers at baseline level (year 2012) the main first pillar CAP measures, like decoupling and payment modulation. This part of the model can sup- port modification in CAP mechanisms, such as the transition from historical to regional- ized decoupling, as well as market price variation. The third module addresses the results of the simulation. The results are stored in spe- cific output files readable by statistical and spreadsheet software by adopting GDX rou- tines (GAMS, 2012). These routines generate an organized output comprising calibrating checks, land and water allocation, and farm economic variable dynamics. 311PMP to assess agroenergy impact on water consumption 5. AquaCrop simulations With regard to the research question of this paper, AquaCrop makes it possible to calculate the amount of water that each crop requires in relation to the level of observed yield. One of its key features concerns the climate characteristics of the area. The model requires information on daily maximum and minimum air temperatures, daily rainfall, daily evaporative demand of the atmosphere expressed as reference evapotranspiration, and the mean annual carbon dioxide concentration in the bulk atmosphere. The climatic parameters adopted were collected from the meteorological measurements recorded in the Emilia-Romagna region, from 1 January 2002 to 31 December 2009. For the annual car- bon dioxide value, we adopted the CO2 concentration in the Mauna Loa Observatory in Hawaii, a monitoring centre used by AquaCrop. The reference value of evapotranspiration was calculated using FAO criteria (Allen et al., 1998) and parameters were determined using the FAO’s model ETo-Calculator (Raes et al., 2009). In AquaCrop, the crop system is represented by five different linked and interrelated modules concerning: phenology, aerial canopy, rooting depth, biomass production and harvestable yield. Each different crop grows and develops over its cycle by expanding its canopy and its rooting system at the same time. Table 1 shows crop growing characteristics used in the evaluation of water requirement. The six crops are the arable crops for which AquaCrop provides information for simulations on response of yield to irrigation water use. Table 1. Crop growing characteristics adopted in the simulation process. Crop Sowing Date Harvest Date Plant Density (plant ha-1) Maize 01-apr 10 August 95 x 103 Silage Maize 01-apr 10 August 85 x 103 Sorghum 13 May 11 September 200 x 103 Sugar Beet 10-apr 29 August 100 x 103 Soybean 1 May 7 September 350 x 103 Tomato 01-apr 19 July 35 x 103 The characteristics of soil (hydraulic conductivity, water content, field capacity, per- manent wilting point) adopted in the study follow the standard profile suggested by the FAO model. This means that soil characteristics within the area are not differentiated. AquaCrop also has an irrigation management component which makes it possible to operate on soil fertility and on the type of irrigation system. Standard model values were used for soil fertility, while the period and the method of irrigation were differentiated according to crop for the irrigation system (Table 2). As noted above, the study focused on the interdependence between water require- ment and biomass in order to analyze how an increase in applied irrigation water (AIW) influences the production in terms of biomass (fresh yields). The AquaCrop model was used to generate non-linear yield response functions to AIW, able to highlight the production trend of the yield level in relation to AIW. The sim- ulations were carried out starting with a rain-field crop, and increasing AIW by season, 312 M. Donati, D. Bodini, F. Arfini, A. Zezza applying constant increments of irrigation water (different for each type of crop), in order to reach the maximum level of biomass. Table 3 reports an example of how results were organized to generate non-linear yield response functions for sugarbeet. Table 2. Irrigation management hypothesized for each crop considered in simulation procedure. Crop Sowing Date Harvest Date Irrigation treatments Method Maize 01-apr 10 August 4 Irr.Schedule - Sprinkler SilageMaize 01-apr 10 August 4 Irr.Schedule - Sprinkler Sorghum 13 May 11 September 2 Irr.Schedule - Sprinkler SugarBeet 10-apr 29 August 8 Irr.Schedule - Sprinkler Soybean 1 May 7 September 3 Irr.Schedule - Sprinkler Tomato 01-apr 19 July 16 Irr.Schedule - Drip Table 3. Production level according to different AIW - Sugarbeet (t/ha). Year Applied irrigation water (mm) 0 50 100 150 200 250 300 350 400 450 500 550 2002 35.1 39.7 43.6 44.6 44.5 44.0 43.5 42.7 42.1 41.6 41.2 41.0 2003 4.2 5.2 9.7 17.1 24.8 32.4 36.5 40.1 41.8 42.8 43.0 42.2 2004 4.6 5.3 13.3 27.4 33.1 38.2 42.5 43.6 44.2 42.7 41.6 41.4 2005 14.5 18.4 27.2 33.1 38.0 41.7 44.2 45.0 44.4 42.2 41.6 41.5 2006 3.6 4.8 9.0 16.7 28.3 36.1 41.4 44.4 45.1 44.8 43.5 42.3 2007 18.7 20.9 24.2 29.4 34.1 38.0 40.6 42.3 43.4 44.2 44.2 43.4 2008 26.2 27.7 30.8 34.4 36.9 39.2 40.7 41.5 42.0 42.2 42.0 41.7 2009 4.9 7.1 13.0 20.7 28.4 34.6 39.2 41.0 42.0 42.0 42.2 42.2 From the data in Table 3 and according to the methodology proposed by García-Vila and Fereres (2012), three non-linear crop-water production functions related to the 20th, the 50th and the 80th percentile were identified. The percentile values made it possible to obtain a non-linear function (polynomial), that reports the contribution of irrigation to the production of fresh yield. Table 4 shows the values of fresh yield obtained for the dif- ferent levels of AIW according to a distribution over three percentiles for sugarbeet. Table 4. Percentiles values from simulation for sugarbeet (t/ha). Percentile Applied irrigation water (mm) 0 50 100 150 200 250 300 350 400 450 500 550 20th 4.0 5.2 9.5 17.0 27.6 34.2 38.7 40.8 41.9 41.9 41.5 41.3 50th 9.7 12.8 18.8 28.4 33.6 38.1 41.0 42.5 42.7 42.5 42.1 41.9 80th 27.9 30.1 33.4 36.4 39.3 42.1 43.6 44.5 44.6 44.3 43.7 42.5 313PMP to assess agroenergy impact on water consumption The results obtained for each crop were represented by polynomial functions in order to better describe the evolution of fresh yield production in relation to the applied irriga- tion water, i.e. the yield-AIW functions. In Figure 2, the vertical axis of each chart shows the level of fresh yield produced (t/ha), and the horizontal axis shows the different levels of applied irrigation water. The results demonstrate that the polynomial functions fit the simulations responses provided by AquaCrop with a high level of R2. Figure 2. Simulated biomass production for 8 years in response to different levels of applied irrigation water (AIW): (a) Sugar beet; (b) Silage Maize; (c) Maize; (d) Sorghum; (e) Soybean; (f ) Tomato. 1   (a) (b) (c) (d) (e) (f) 2    3   314 M. Donati, D. Bodini, F. Arfini, A. Zezza 6. Impact assessment 6.1 Land allocation The integrated model was developed in the agricultural context of a lowland area of Emilia Romagna (Northern Italy). For the purposes of clarity and brevity, the impact assessment presents the results of a homogenous area, the province of Parma. The latent crop introduced in the simulation phase was sorghum for bioethanol production. We assume that this variety of sorghum is currently not present in the regional production plan. The information concerning average yields, price and specific production cost for sorghum was taken from an experimental study promoted by the Emilia Romagna Region aiming to evaluate the possibility of building a regional supply chain for bioethanol. The results presented were developed in the framework of the “Health Check” CAP reform. The market simulations consist of increasing the price of sorghum by 1 €/t in 200 steps, to reach the maximum price of 200 €/t. The results identify the economic threshold for sorghum, that is the starting price from which sorghum for biomass can be inserted into the production plan of farms in the province of Parma. Figure 3 shows the total area that would be grown with sorghum in the region at different level of sorghum price. The graph presents a curve that starts to increase from a level of 58 €/t, the profitable thresh- old for the crop. The price dynamics causes a big increase in sorghum acreage. The rota- tional constraints and the complementary and substitution relationships within the cost matrix prevent farm surfaces from becoming specialized in a single crop. The graphical representation of the simulation shows the different production levels with regard to dif- ferent price levels, so that it is possible to identify the price to pay to producers in order to obtain a certain quantity of raw material. So, for example, if the supply chain needs 20.000 ha of sorghum, the price that should be paid to farms is more or less 108 €/t. The simulation also shows variation in the relative incidence of sorghum compared to other crops in the regional production plan. Figure 4 highlights that the increase in the incidence of sorghum on the regional production plan is due to a big fall in the incidence of fodder crops and wheat. It is important to note that in this study, the prices of other crops are assumed to be constant throughout, but likely market price modifications would produce a change in profitability across the model, with different impacts on sorghum production responses. 6.2 Water consumption The simulations carried out by the integrated model capture the relationships between yields and the applied irrigation water. For the six crops in Table 1, this relationship was estimated by using AquaCrop as previously described, while for the other crops informa- tion has been taken from the literature. For sugarbeet, maize, silage cereals, tomato, soya and sorghum, yields are non-linearly related to AIW, while for the other irrigated crops (e.g. alfalfa) an exogenous fixed relationship was assumed. The polynomial regressions estimated for the six crops are used in the AquaCrop model to identify the quantity of irrigated water each farm needs to obtain the observed yield. Polynomial functions from the 50th percentile were used. The difference in terms of AIW among farms is reflected in 315PMP to assess agroenergy impact on water consumption the observed yield. To adapt the yield-AIW functions to the observed information, a cor- rection factor was applied as follows: = ⋅yield F w v( )n j n j n j j, , , (11) Figure 3. Evolution of the sorghum hectares due to price variation (Parma province). 0 20 40 60 80 100 120 140 160 180 200 0 5,000 10,000 15,000 20,000 25,000 30,000 35,000 40,000 45,000 50,000 (Ha) (€ /t ) Figure 4. Evolution of the specific incidence of each crop on the regional production plan (Parma province). 0 10 20 30 40 50 60 70 80 90 1 10 19 28 37 46 55 64 73 82 91 100 109 118 127 136 145 154 163 172 181 190 199 €/ton    -­‐  SORGHUM In c.  % Wheat Maize Other  cereals Fodder  crops Sorghum Others 316 M. Donati, D. Bodini, F. Arfini, A. Zezza where the yield for each irrigated crop j (for j=1,2,…,J) and each farm n (for n=1,2,…,N) is derived from the non linear function obtained by AquaCrop simulations multiplied by a weighting value v j calculated as the ratio between the average observed yield and the average yield provided by AquaCrop simulations. wn j, the variable related to AIW to be estimated for each farm and crop. This correction makes it possible to retain the shape of the yield-AIW function and model feasibility. In fact, the difference between observed and experimental yield data can be unavoidable (Dillon and Hardtacker, 1993; Cortignani and Severini, 2009). Table 5 presents the AIW estimated levels for the six crops simulated in AquaCrop. From experimental tests, sorghum yield is assumed to be 23 t/ha correspond- ing to 121.1 mm of AIW. Table 5. Estimated AIW per crop in Parma province. Area Altitude Farm Type Class of Size (ha) No. of farms Sugarbeet Silage cereals Maize Tomato Soya Sorghum AIW (mm) Parma Plain Arable crops 0-10 1233 158.8 28.3 159.8 86.9 121.6 121.1 Parma Plain Arable crops 10-20 452 185.4 234.5 169.4 131.4 121.1 Parma Plain Arable crops 20-30 175 118.8 8.9 171.7 240.7 151.6 121.1 Parma Plain Arable crops 30-50 129 220.4 32.1 21.2 211.2 168.0 121.1 Parma Plain Arable crops 50-100 97 271.9 418.9 63.7 151.3 407.8 121.1 Parma Plain Arable crops 100-300 20 51.4 145.2 234.5 240.7 121.1 Parma Plain Arable crops > 300 1 30.1 121.1 Parma Plain Horticulture 0-10 10 158.6 86.9 121.1 Parma Plain Horticulture 10-20 6 185.4 169.4 121.1 Parma Plain Horticulture 20-30 8 118.8 172.1 240.7 121.1 Parma Plain Horticulture 30-50 14 220.4 21.2 211.1 121.1 Parma Plain Horticulture 50-100 10 271.9 63.7 151.3 121.1 Parma Plain Dairy 0-10 90 160.0 211.1 121.1 Parma Plain Dairy 10-20 127 185.5 630.5 234.5 160.1 121.1 Parma Plain Dairy 20-30 121 118.8 8.9 171.7 240.7 121.1 Parma Plain Dairy 30-50 180 220.4 32.1 21.2 196.3 121.1 Parma Plain Dairy 50-100 97 271.9 418.9 63.7 164.9 121.1 Parma Plain Dairy 100-300 28 51.4 145.2 234.5 211.1 121.1 Parma Plain Dairy > 300 6 13.8 418.9 234.5     121.1 Parma Plain All All 2804 174.6 169.3 234.5 196.4 154.7 121.1 For five crops out of six, the AIW changes according to the farm, revealing a direct relationship with the specific crop yield. But for sorghum the AIW level is the same for all farms in relation to the supposed uniform yield in the area. The relationship yield-AIW was used to evaluate the impact of the potential introduc- tion of sorghum for biomass production in terms of water demand. Price rises may per- suade farmers to convert cereal land and grassland to sorghum with a negative effect on 317PMP to assess agroenergy impact on water consumption the water consumption. In particular, the process of substitution of a non-irrigated crop, like wheat, with an irrigated crop, like sorghum, generates an increase in total irrigated water for the province of Parma. More specifically, when a price of 70 €/t for sorghum is applied, about 30% of the wheat disappears in favour of sorghum. This shift produces an increase of the total irrigated water of 2.5%, which corresponds to an increase of 1.5 mil- lion m3 compared to before sorghum was introduced (see Figure 5). Figure 6 shows the trend in consumption of water for irrigation. Up to a price of 61 €/t irrigated water consumption declines by 0.5%, due to a process of substitution between sorghum and grassland with high AIW. This substitution is interrupted by the demand for forage from dairy farms, represented by a specific constraint in the model. Figure 6. Dynamics in water consumption with sorghum cultivation (province of Parma). 52,000 54,000 56,000 58,000 60,000 62,000 64,000 66,000 68,000 70,000 72,000 57 61 65 69 73 77 81 85 89 93 97 10 1 10 5 10 9 11 3 11 7 12 1 12 5 12 9 13 3 13 7 14 1 14 5 14 9 15 3 15 7 16 1 16 5 16 9 17 3 17 7 18 1 18 5 18 9 19 3 19 7 Price  Scenarios  -­‐  €/t. .0 00  m 3 Figure 5. Irrigated water allocation. 0 10,000 20,000 30,000 40,000 50,000 60,000 70,000 57 62 67 72 77 82 87 92 97 10 2 10 7 11 2 11 7 12 2 12 7 13 2 13 7 14 2 14 7 15 2 15 7 16 2 16 7 17 2 17 7 18 2 18 7 19 2 19 7 Price  Scenarios  -­‐  €/t. .0 00  m 3 Sorghum Grassland Soya Silage  maize Tomato Maize Sugarbeet 318 M. Donati, D. Bodini, F. Arfini, A. Zezza The model captures the impact of farm production decisions on the allocation of water, but it does not capture efficiency in the use of irrigated water. An implicit assump- tion of the model is that farmers in the province take the production decisions consider- ing that water is available, and water allocation is a consequence of their decision. From this point of view, the calibrated AIW represents the willingness to use water for each crop. This implies that there are no direct constraints on the total available water per farm and, therefore that trade of water entitlements cannot be evaluated in this model. 7. Conclusions This paper uses positive mathematical programming model to evaluate the effects of agro-energy crop cultivation on land and water allocation. The agro-energy crop is sorghum for the production of second generation bioethanol. The model simulates the introduction of sorghum as a new crop in the farm production plan in the lowlands of Parma, an agrar- ian region in northern Italy, and evaluates crop economic threshold and the change in farm crop set and irrigated water use. The model considers latent crop information in each farm production plan, so that in the simulation phase, new crops, such as sorghum for biomass production, are considered a new option with regard to the crops already activated in the observed situation (Arfini and Donati, 2013). The second PMP phase estimates non-line- ar cost functions where information about the latent crop is included. The main constraint of the model is related to the total farm land. Water appears in the third phase through the yield-AIW function which defines the relationship between crop yields and the quantity of AIW. The yield-AIW function is estimated by simulating the yield response of the crop in relation to different AIW quantities; simulations were carried out by implementing AquaCrop (Steduto et al., 2009) for six crops: sugarbeet, maize, silage cereals, tomato, soya and sorghum. The simulation results revealed the capacity of the model to evaluate the productive potential of sorghum for biomass in the region and the consequences on the regional pro- duction plan. The scenarios developed to estimate the sensitivity of the sorghum in rela- tion to a variation in its market price made it possible to identify the price level necessary for providing a sufficient quantity of raw material to feed a second generation bioetha- nol supply chain. The increase in sorghum price leads to a reduction in wheat and fodder crops; maize is also affected by a process of substitution. The PMP model makes it possible to evaluate the impact of sorghum cultivation on water consumption for irrigation. The results show that increasing land use for bio- mass sorghum in the lowlands of Parma province, a dairy and horticultural area, entails an increase in the total irrigated water quantity due to a decline in non-irrigated crops, like wheat. The information about water consumption will be useful for policy makers in preventing excessive demand for water in areas where the risk of drought might gener- ate non-efficient water allocation. Furthermore, the model can evaluate the impact of the CAP on irrigation water needs as well as land use. The 2014-2020 CAP reform, which emphasized environmental objectives in agricultural policies, could also be evaluated in order to predict effects on water resources using a similar approach. The model does however present certain limitations. The AquaCrop simulation was run on the crops included in the FAO model dataset. This is because yield-AIW functions for fodder crops, such as alfalfa and permanent meadows, have not been estimated. 319PMP to assess agroenergy impact on water consumption Furthermore, the PMP model presented in this study does not consider an explicit water constraint, and it is thus not possible to simulate reductions in water availability and identify the dual value associated with this resource. Optimizing the model to include this information would provide a value for levels of water scarcity in certain groups of farms or areas, and would be useful in identifying an appropriate price. Finally, a water con- straint in the model would make it possible to simulate the water entitlement market and put in place economic tools to make the use of water for irrigation more efficient. Acknowledgments An earlier version of this paper was presented at the 2nd AIEAA Conference “Between Crisis and Development: which Role for the Bio-Economy”, 6-7 June, 2013, Parma, Italy. The present work was carried out in the framework of the project BIOSEA “Optimization for bioenergy supply chains for an economic and an environmental sus- tainability” financed by the Italian Ministry of Agriculture and Forestry (www.biosea.dista. unibo.it). We thank the Azienda Agraria Sperimentale “Stuard” of Parma, Italy for provid- ing experimental data on sorghum H133 for biomass. The authors wish to thank the two anonymous referees and the Editor for their helpful suggestions and constructing com- ments on an earlier version of the paper. References Allen, R.G., Pereira, L.S., Raes, D., and Smith M. (1988). Crop evapotraspiration. Guide- lines for computing crop water requirements. Irrigation and Drainage Paper no. 56. Rome: FAO. Arfini F., Donati M., and Zuppiroli M. (2005). Un modele quantitatif pour l’evaluation des effets de la riforme de la PAC en Italie. Economie Rurale 285: 70-87. Arfini, F., and M., Donati (2011). The impact of the Health Check on structural change and farm efficiency: A comparative assessment of three European agricultural regions. Disaggregated impacts of CAP reforms. Paris: OECD, 75–90. Arfini, F., and Donati, M. (2013). Organic Production and the Capacity to Respond to Market Signals and Policies: An Empirical Analysis of a Sample of FADN Farms. Agroecology and Sustainable Food Systems 37: 149-171. Arfini, F., Donati, M., and Zuppiroli, M. (2005). Un modèle quantitatif pour l’évaluation des effets de la réforme de la Pac en Italie. Economie Rurale, 285: 70–87. Arfini, F., and Paris, Q. (1995). A positive mathematical programming model for regional analysis of agricultural policies. In Sotte, F. (ed.). The Regional Dimension in Agri- cultural Economics and Policies, EAAE, Proceedings of the 40th Seminar, June 26-28, Ancona: 17-35. Bartolini F., Bazzani G.M., Gallerani V., Raggi M., and Viaggi D. (2007).The impact of water and agriculture policy scenarios on irrigated farming systems in Italy: An analysis based on farm level multi-attribute linear programming models. Agricultur- al Systems, 93: 90-114. Bazzani G.M., Di Pasquale S., Gallerani V., Morganti S., Raggi M., and Viaggi D. (2005). The sustainability of irrigated agricultural systems under the Water Framework Directive: first results. Environmental Modelling & Software 20: 165-175. 320 M. Donati, D. Bodini, F. Arfini, A. Zezza Bazzani G.M., Di Pasquale S., Gallerani V., and Viaggi D. (2004). Irrigated agriculture in Italy and water regulation under the European Union water framework directive. Water Resources Research 40/7:1-8. Blanco, M., Cortignani, R., and Severini, S. (2008). Evaluating Changes in Cropping Patterns due to the 2003 CAP Reform. An Ex-post Analysis of Different PMP Approaches Considering New Activities. Paper prepared for presentation at the 107th EAAE Seminar “Modelling of Agricultural and Rural Development Policies”, Seville. Brooke, A., D. Kendrick, A. Meeraus, R. Raman, and R. E. Rosenthal (2011). GAMS: A User’s Guide. Washington, DC: GAMS Development Corporation. Cortignani, R., and Severini, S. (2009). Modeling farm-level adoption of deficit irriga- tion using Positive Mathematical Programming. Agricultural Water Management 96: 1785-1791. Dillon, J.L., and Hardtacker J.B. (1998). Farm management research for small farmer development. Rome: FAO. Doornbosch, R., and Steenblik, R. (2007). Biofuels: is the cure worse than the disease?. Round Table on Sustainable Development, Paris: OECD Publishing. Dorenbos, O.T., and Kassam, A.H. (1979). Yield response to water. Irrigation and Drain- age Paper n., 33. Rome: FAO. European Commission (2006). An EU Strategy for Biofuels, Brussels. European Commission (2009). Why do we need a Common Agricultural Policy?, Brussels. European Commission (2010). EUROPE 2020. A strategy for smart, sustainable and inclu- sive growth, Brussels. European Parliament and European Council (2010). Directive 2000/60/EC of the Europe- an Parliament and of the Council of 23 October 2000 establishing a framework for Community action in the field of water policy, Luxembourg. GAMS (2012). GAMS GDX facilities and tools. Washington, DC: GAMS Development Corporation. García-Vila M., and Fereres E. (2012), Combining the simulation crop model AquaCrop with an economic model for the optimization of irrigation management at farm lev- el. European Journal of Agronomy 36: 21-31. Graveline, N., and Mérel, P. (2012). How do farmers adapt to water scarcity? Intensive Margin Adjustments in Beauces’ agriculture. Paper presented at the EcoProd Semi- nar of INRA, Montpellier. Hardin, G. (1968). The tragedy of the commons. Science 162: 1243-1248. Harou, J.J, Pulido-Velazquez, M., Rosenberg, D.E., Medellín-Azuara, J., Lund, J.R., and Howitt, R.E (2009). Hydro-economic models: Concepts, design, applications, and future prospects. Journal of Hydrology 375: 627-643. Heckelei, T., Britz, W., and Zhang, Y. (2012). Positive Mathematical Programming Approaches – Recent Developments in Literature and Applied Modelling. Bio-based and Applied Economics 1: 109-124. Heckelei, T., and Wolff H. (2003). Estimation of constrained optimisation models for agri- cultural supply analysis based on generalised maximum entropy. European Review of Agricultural Economics 30(1): 27-50. Helming, J., and Peerlings, J. (2003), Effects of EU dairy policy reform for Dutch agriculture and economy; applying an agricultural programming / mixed input-output model, Paper presented at 25th International Conference of Agricultural Economists, Durban. 321PMP to assess agroenergy impact on water consumption Howell, T.A. (2001). Enhancing water use efficiency in irrigated agriculture. Agronomy Journal 93: 281-289. Howitt, R.E. (1995). Positive mathematical programming. American Journal of Agriculture Economics 77: 329–342. Medellín-Azuara, J., Harou, J.J., and Howitt R.E. (2010). Estimating economic value of agricultural water under changing conditions and the effects of spatial aggregation. Science of the Total Environment 408: 5639-5648. OECD (2010). Water resources in agriculture: outlook and policy issues. Paris: OECD Publishing. Paris, Q. (2012). Economic Foundations of Symmetric Programming. Cambridge Univer- sity Press, New-York. Paris, Q., and Arfini, F. (2000). Funzioni di costo di frontiera, auto-selezione, rischio di prezzo, PMP e Agenda 2000. Rivista di Economia Agraria 2: 211-242. Paris, Q., and Howitt, R.E. (1998). An analysis of ill-posed production problems using maximum entropy. American Journal of Agricultural Economics 80: 124-138. Petersen J.-E. (2008). Energy production with agricultural biomass: environmental impli- cations and analytical challenges. European Review of Agricultural Economics 35: 385-408. Raes. D., Steduto, P., Hsiao, T.C., and Fereres E. (2009). AquaCrop – The FAO crop model for predicting yield response to water: II. Main algorithms and software description. Agronomy Journal 101: 438-447. Röhm, O., and Dabbert, S. (2003). Integrating Agri-environmental programs into regional production models: An extension of positive mathematical programming. American Journal of Agricultural Economics 85(1): 254-265. Rosegrant, M.W., Ringler, C., McKinney, D.C., Caia, X., Keller, A., and Donosod, G. (2000). Integrated economic–hydrologic water modeling at the basin scale: the Maipo river basin. Agricultural Economics 24: 33-46. Steduto, P., Hsiao, T.c., RAes, D., and Fereres, E. (2009). AquaCrop – the FAO crop model to simulate yield response to water: I. Concepts and underlying principles. Agrono- my Journal 101: 426-437. Ward, F. A., and A.M., Michelsen (2002). The Economic Value of Water in Agriculture: Concepts and. Policy Applications. Water Policy 4: 423-446. Zezza, A. (2008). Bioenergie:quali opportunità per l’agricoltura italiana, Collana Studi e Ricerche INEA, Edizioni Scientifiche Italiane, Napoli.