Login   /   Register
Baldini, Danuso, Rocca, Bulfoni, Monti, and De Mastro: Jerusalem Artichoke (Helianthus tuberosus L.) Productivity in Different Italian Growing Areas: a Modelling Approach

Jerusalem Artichoke (Helianthus tuberosus L.) Productivity in Different Italian Growing Areas: a Modelling Approach


Jerusalem artichoke is considered one of the most interesting crops for inulin production. It has been used to date for the production of low calorie sweetening syrups, dietetic food enriched with fibre, medicines and cosmetics, while more recently, interesting prospects have been opening up for energy uses. The main aspects influencing its adaptability to different pedo-climatic conditions and cropping systems were analysed by implementing CSS (Cropping System Simulator, a simulation model describing crop biomass production) for this specific crop. Growth analysis experimental data of plant dry matter accumulation, obtained over two years of trials (1999-2000) in different Italian growing areas (Udine, Bologna, Bari) under irrigated and rain-fed conditions, were used for the parameterisation and calibration of the model. The biomass accumulation observed and simulated under rain-fed and irrigated conditions in the different growing areas is reported for the different plant organs, with good correspondence shown between simulated and measured values as reported by the statistical indices for the model calibration, particularly for biomass of tubers and leaves. The model studied, despite a simplified description of some processes, proves to represent the maximum biomass yield of Jerusalem artichoke satisfactorily, with an adequate response to the main environmental factors causing yield and biomass production variation among the years and locations. However, further model improvements are necessary in order to better represent the relationship between phenology and translocation of the assimilates between stalk and tuber during the development phases of the plant, suggesting shorter survey intervals over this growing phase.


Jerusalem artichoke can be considered as both a sugar and dietary fibre crop, as it accumulates linear polymers of fructose (fructans, also known as inulin) in its roots, tubers and stalks, with a highly variable degree of polymerisation (4-150 DP), which affects their end-use. Fructans with low DP and fructose, obtained by hydrolysis of the fructans, are often used directly as dietary fibres or as low-cal orie sweetening syrups, while inulin with high DP (10-30) is used entirely for industrial and non-food uses such as pharmaceuticals and cosmetics (Fuchs, 1993; Danuso, 2001). One of the main problems of today’s sugar processing industry is the need to extend the harvest season, in particular by anticipating it, with compensation given to farmers for the lower yield of earlier harvests, as in Belgium where the industry pays a higher price to compensate for this lower production (Danuso, 2001).

Jerusalem artichoke is relegated today to small areas mainly for the production of tubers for human or livestock consumption, so there could be interesting prospects for high earning potential with an early whole-plant harvest (stalks plus tubers); however, there are still problems connected to the perfecting of harvest machinery and processing methods (Baldini et al., 2006). The whole-plant harvesting method has always produced higher yields than the traditional harvesting of tubers alone, with the agronomic advantage of freeing the land earlier (Paolini et al., 1996; Baldini et al., 2004). The crop can be cultivated as: i) an annual crop, with the harvest of stalks and tubers together, or ii) a multi-year crop, with only the aerial biomass harvested each year (D’Egidio et al., 1998; Baldini et al., 2006). These considerations and the current need to find new renewable energy sources open new prospects for Jerusalem artichoke as a biomass crop for energy uses, particularly for liquid biofuel production (bioethanol; Curt et al., 2006), methane from anaerobic digestion (Lehtomaki et al., 2008) and gas from pyrolysis (Encinar et al., 2009).

A more detailed understanding of the distribution of assimilates during plant growth among different organs becomes necessary in order to predict, through scenario analysis and simulation experiments, growth and production of stems and tubers depending on pedo-climatic conditions, cropping techniques and genetic material. This may be particularly useful in analysing the feasibility of such possible crop uses in different environments. For this reason, the aim of our study was the improvement and calibration of the Cropping System Simulator (CSS) model (Danuso et al., 1999) to enable it to simulate Jerusalem artichoke productivity in different conditions (different Italian growing areas). CSS, in fact, has already been parameterised and calibrated with encouraging results for the most common annual crop (Danuso et al., 2009). CSS is a generic crop simulation model implemented with the Simple Easy to use Modelling Language (SEMoLa) software, an application for the development of simulation models and agro-ecological knowledge integration (Danuso, 2003). Both SEMoLa and CSS have been developed and are maintained at the Department of Agricultural and Environmental Sciences of the University of Udine, Italy. The calibration performed in our study has been done using data obtained from experimental trials conducted in different Italian areas within the framework of the PRIN Project (MIUR), entitled Colture per la produzione di inulina: modelli di risposta ambientale e strategie colturali.

Materials and Methods

Experimental trials

Trials on Jerusalem artichoke (Helianthus tuberosus L., variety Violet de Rennes) were conducted over the two years 1999 to 2000, by the Working Groups taking part in the PRIN project, in three areas of Italy: Udine, Policoro–Bari and Cadriano–Bologna (Table 1). In each environment, the effects of two different production factors were evaluated: time of harvest and irrigation regime with two treatments: i) replacement of total evapo-transpiration and ii) dry regime (rainfed crop with the aid of irrigation).

Tubers (average weight 50-70 g each) where manually planted at a depth of 3-4 cm, in rows 0.7 m apart and with a distance between plants in the row of 0.25 m (0.20 m in Bologna), giving a planting density of 5.7 tubers/m2 in Bari and Udine and 7 tubers/m2 in Bologna. The trial design was a split-plot with four replications. The experimental unit was 5 rows of 12 m in length in Bologna, giving an area of 42 m2, and 8 rows, 8 m long in Udine and Bari for an area of 44.8 m2. Dry matter accumulation was evaluated after collecting, by hand, five plants at different plant phenological phases for each plot, covering the total growing season until the standing crop was dried naturally. In particular, the number of samplings varied from 6 to 11, depending on the different locations, with about 10-15 days between two consecutive harvests. At each sampling, the fresh and dry weights of leaves, stems (with ramifications) and tubers (without stolons) were measured.

Model implementation

CSS is a generic simulator of the cropping system, with a daily time step, already used for forecasting and risk evaluation in different cropping areas (Danuso et al., 2009). The model has been implemented by the SEMoLa language (version 5.8; Danuso, 2003), which can be compiled in different ways in order to produce stand-alone executable programmes. CSS is freely available as source code written in the SEMoLa language or as an executable program. An important characteristic of this model is that, because it is created with this language, it can be modified and updated very easily.

CSS is formed by interconnected modules that simulate soil water balance (using a simple two-layer cascade method), nitrogen dynamics, crop phenology and growth as affected by environmental, pedological and agricultural factors (Figure 1). It is a generic model because the same equation structure allows the simulation of different crops (e.g. wheat, maize, soybean) in relation to crop parameters that have to be specifically calibrated. It consists of different sub-models (modules), among which the most relevant to crop production is CSS_CropYield, which simulates crop phenology (growing phases depending on growing degree, days accumulation and day length), aerial and root biomass (Wroot), crop yield and leaf area dynamics [based on the Simple and Universal CROp growth Simulator (SUCROS) model; van Laar et al., 1997]. The aerial biomass (or crop biomass, Wcrop) is partitioned in leaves (Wleav), stalks (Wstem) and storage organs (Wstor). CSS_CropYield describes the growth from photosynthesis and respiration, and allocates the daily dry matter increments to the different organs according to the partitioning factors introduced as a function of the development stage of the crop.

Soil is considered as divided into two layers (Figure 2). The depth of the upper layer (Dr, root layer) changes, according to the crop root growth, from the sowing depth to the maximum depth exploitable by the roots during the crop cycle (Ds). The soil layer between the rooting depth and maximum exploitable soil depth is the deep layer (Dd), maximum at crop sowing and decreasing as the crop roots continue growing. The maximum depth that can be reached by the root system depends on the crop characteristics, the presence of compact layers, gravel, rocks, or shallow water table.

CSS simulates drainage to the water table and capillary rise from it, according to Driessen and Konijn (1992). If the depth of the phreatic water table bed is known, it also gives a simple simulation of the dynamics of water table depth.

Crop residue decay produces new humus and mineral compounds; they are divided into easily decomposable residues and those resistant to decomposition, with different mineralisation and humus synthesis coefficients. Soil microbial biomass is also considered. The balance of soil organic matter and residues is simulated, for the working layer, by an implementation of the RothC model (Coleman and Jenkinson, 2008).

The CSS model has been modified in order to improve its capacity to represent crops like Jerusalem artichoke or potato, in which the translocation of assimilates from stalks to tubers, at the end of the cycle, is relevant. This is a specific physiological process not existing in grain crops like maize or wheat. The decreasing of leaves and stem weight at the late growth stages owing to senescence and translocation has been described in CSS through two interrelated aspects: i) no more biomass is allocated to leaves and stems after the physiological maturity stage and ii) part of their biomass is lost through senescence while the rest is translocated to tubers.

Model parameterisation and calibration

In our study, CSS has been parameterised and calibrated for the Jerusalem artichoke by comparing the simulation results for biomass accumulation (t/ha dry matter) in the different plant organs (leaf, stem and storage organ) with the experimental data obtained in the different locations. The first step of the parameterisation procedure has been done using information from previous research, common knowledge and experience; the second step with a trial-and-error procedure. The last step, for a more accurate calibration, has been performed using the automatic calibration routine of the SEMoLa software that uses an iterative procedure (Gauss-Newton linearisation method; Beck and Arnold, 1977; Draper and Smith, 1981) for the minimisation of the residual sum of squares between simulated and observed data.

This procedure was carried out by setting several multiple simulations for each scenario, following the real situation (pedo-climatic conditions and cropping practices) in the experimental trials for each location and thesis. To avoid the influence of crop parameters on calibration, they have been changed or fixed together for all scenarios.

Input data for simulations

To perform simulations, CSS requires meteorological data, soil parameters, crop parameters and parameters for the cropping scenario to be represented (Figure 1). These data have to be organised into four types of input data files: exofile containing exogenous variables (meteorological data), parfile with soil and crop parameters, evtfile declaring time, type and modality of every cropping practice (events) and actfiles (one for each event type), containing the actions or parameters of each application of cropping techniques. At run-time and before the beginning of the simulation, parameter values of the parfile are set and their values usually remain the same throughout the simulation; however, their values can be changed during the crop cycle by the event instances of the evtfile; these instances make reference to specific sets of parameter values in the corresponding actfile and then immediately modify the current parameter values. The modifications of the parameter values by events are considered the actions of the event; in this case, param eters act in the system like switches. For example, the simulation can start with the crop parameters for fallow; when the event Planting occurs, all the crop parameters (for phenology, light interception, growth, etc.) are changed, depending on the sown/planted crop. This allows the simulation of crop rotations with different types and amounts of crop inputs. Every event type of the evtfile has a related actfile, which contains various sets of parameters for the different application modes of the event. The user (and not just the modeller) can also change and increase the number of parameter sets in actfiles. In the case of organic fertilisation, the event OrgFert has the actfile OrgFert.act, which contains parameter sets specifying the characteristics of different organic fertilisers like manure, straw, slurry, etc.

Meteorological data

The daily meteorological data required in CSS to run simulations are minimum and maximum air temperatures (°C), precipitation (mm/d), reference evapo-transpiration (mm/d) and global radiation (MJ/m2/d), which are saved in exofile. These data have to be complete and without outliers; to avoid this problem, meterological data have been checked using Climatica software (Danuso and Sandra, 2006), which assists in the creation of the correct meteorological data file by automatic correction of the names of variables, changes the date format, and rebuilds missing data. It also advises when files are not correct for missing days, wrong time order and data out of range. The model calibration has been performed using the historical meteorological data corresponding to the years of the trials (1999-2000) in the three localities. The trends of maximum and minimum temperature and rainfall reported in Figure 3 show the different meteorological conditions.

Soil and crop parameters

The main soil characteristics measured for each site and used for the simulations are given in Table 2. The soil in Udine was very shallow (about 50 cm) and without groundwater useful for the crops, while in both Bologna and Bari the contribution of groundwater was important. According to the specific crop selected, crop parameter values were adjusted considering a reasonable range of variation, as dictated by previous research, knowledge or experience. A list of the crop parameters resulting from the parameterisation and used for the final simulations is reported in Tables 3 and 4. This information about soil and crops was used to customise the parameter files (parfile).

Crop management simulation

In CSS the cropping techniques are considered as events; that is, as phenomena that happen and instantaneously modify parameters and states of the system. At present, the following events can be selected to build cropping scenarios: planting, organic and mineral fertilisation, irrigation, harvesting, residue chopping, harrowing, hoeing, extirpation, chiselling and ploughing. The sequences of cropping techniques are saved in event files (evtfiles). The agronomic techniques adopted during experimental trials are inserted in the evtfiles to perform the simulation (Tables 5 and 6).


The agronomical results for the Jerusalem artichoke obtained in the two years of trials, about yield and biomass accumulation for the different crop organs, emphasise that in Bologna, thanks to the good pedo-climatic characteristics, irrigation appeared superfluous. The highest production of the entire research project was obtained corresponding to a very late harvest time; for the final yield of Jerusalem artichoke tubers, the Bari environment was the most productive. This crop adapts well to a hot arid climate, such as that in the south of Italy, on condition that the limited water availability is compensated for by good soil depth and a water table within reach of the exploratory capacities of its root apparatus (Baldini et al., 2006). The integration of the specific algorithm in CSS in order to represent translocation of assimilates from stalks to tubers was necessary because this specific physiological process is typical of this crop. Results demonstrate that this model adjustment does improve the simulation results especially in the last phase of the crop growth. Despite this model improvement to estimate a better distribution of biomass in the storage organs, further ameliorations are desirable especially considering the physiological factors.

The most relevant results of this study concern the model calibration, which led to good simulation results. In Figures 4 and 5, the biomass accumulation observed and simulated under rain-fed and irrigated conditions is reported for the different plant organs, where good correspondence between simulated and measured values is noted, particularly for biomass of tubers and leaves. In stalks, the simulated accumulation trend turns out to be slightly more advanced in comparison with that obtained experimentally. An inexact representation of the timing of storage of assimilates by the stems and the further reallocation of the storage reserves to the tuber when growing rapidly can be observed.

Differences observed between measured and simulated data for biomass accumulation in tubers during the last growing phase in Udine and Bologna in 1999 could be explained with the quick physiological changes owing to climatic variations. In fact, in October and November of 1999, a sudden temperature decrease corresponding to the last harvest times can be observed (Figure 3). This fact could have triggered a translocation rate increase, accelerating the remobilisation of reserves from stalks to tubers. This becomes more evident when the harvest was delayed in Bologna during the first year. These results emphasise the difficulties found in survey methodology and model implementation to represent these particular physiological adaptations, suggesting shorter survey intervals during this growing phase, when the translocation rate is increasing. On the other hand, comparing the total biomass accumulation for Jerusalem artichoke in the different growing areas, during the same crop season and under rain-fed conditions (Figure 6), it is evident that the model is able to represent the maximum yields in biomass satisfactorily. In some cases, the irrigation in the trials had a very limited effect on the biomass yield; this has also been well reproduced by the simulation (Figure 6). In these environments, in fact, the effect of the irrigation regime was very limited owing to the favourable climatic conditions, good soil characteristics or shallow water table.

Table 7 shows some statistical indices (root mean square error RMSE, modelling efficiency EF, coefficient of determination CD, coefficient residual mass CRM and maximum absolute error MaxAE) for the model calibration (Janssen and Heuberger, 1995) for the different growing areas, thesis and crop seasons, obtained by the SEMoLa framework. The lowest limit of RMSE is 0, that means full adherence between model estimates and measures. The use of EF is an extension to any class of models of the common R2 statistics of regression lines. EF can have either positive or negative values, 1 being the upper limit, while negative infinity is the theoretical lower bound. CD is the proportion of total variance of measurements explained by the estimates but it is not the same as R2, being possible values of CD greater than 1 (CD=1 is the best; that is, the deviation from the mean of measurements is the same for estimates and measurements). MaxAE is most sensitive to outliers, and in fact is the worst case measure. The indices MaxAE and CRM can be either positive or negative, zero being the optimal value.

Discussion and Conclusions

The capability of the CSS model to represent well the experimental results on Jerusalem artichoke obtained in very different environmental and agronomic conditions (from northern to southern Italy) suggests its usefulness as a tool for the planning and risk evaluation of introducing Jerusalem artichoke as a crop for energy or inulin production. The model, despite a simplified description of some processes (e.g. soil water dynamics), seems to be adequately responsive to the main environmental factors causing yield and biomass production variation among years and locations.

Indeed, the simulations confirm, at least for the yield in total biomass, the possibility of using the model to evaluate the different crop management techniques and Italian pedo-climatic conditions, even if the model has been developed primarily to simulate annual herbaceous crops and yielding seeds but not tubers and stalks as in the case of Jerusalem artichoke. However, further model improvements are necessary in order to better represent phenology, partitioning and translocation of the assimilates among the plant organs. In particular, the source-sink relationship between stalk and tuber during the development phases of the plant and in different environments will have to be clarified better and modelled. In fact, the tubers are grown by both current photosynthates and remobilisation of reserves from other plant parts, mainly from the stalk. The transfer from stalk to tuber is known to be up to 50% of final tuber dry weight, depending on many factors: temperatures, flowering time of the cultivar, aerial structural growth, tuber sink capacity, with complicated competing sinks, at plant level, and with changing hierarchical relations during the crop cycle (Denoroy, 1996). Efforts have to be made to understand better the mechanism of distribution of assimilates among aerial structures and tubers, in order to identify the optimal period to obtain the maximum biomass yield, adopting an integral or an aerial harvest, within the perspective of considering Jerusalem artichoke as a crop for producing raw material for energy use at competitive prices.

In conclusion, CSS allows several aspects related to the cropping system to be simulated. However, it is necessary to calibrate the model for each specific crop. Moreover, the development of a more detailed crop parameter database is desirable for further evaluation and improvement of the model for crops with high translocation of assimilates. CSS is also a fundamental part of the farming system simulation X-Farm model (Danuso et al., 2010) where, using new concepts, it can be run automatically for any number of fields on a farm.


Field trials were conducted within the framework and with the financial support of the miur prin 1999 Project Colture per la produzione di inulina: modelli di risposta ambientale e strategie colturali, while crop model development has been possible thanks to the Projects Filiere agroenergetiche in Friuli Venezia Giulia: valutazione della sostenibilità economica, energetica e ambientale a livello aziendale e territoriale, which were funded by the Friuli Venezia Giulia Region, ex art 17, LR n.26, 10/11/2005.



R.G. Allen, L.S. Pereira, D. Raes, M. Smith 1998. Crop evapotranspiration - guidelines for computing crop water requirements. FAO Irrigation and Drainage paper 56. Food and Agriculture Organization of the United Nations, Rome, Italy.


M. Baldini, F. Danuso, A. Monti, M.T. Amaducci, P. Stevanato, G. De Mastro 2006. Chicory and Jerusalem artichoke productivity in different areas in Italy, in relation to water availability and time of harvest. Ital. J. Agron. 1:291-309.


M. Baldini, F. Danuso, M. Turi, G.P. Vannozzi 2004. Evaluation of new clones of Jerusalem artichoke (Helianthus tuberosus L.) for inulin and sugar yield from stalks and tubers. Ind. Crops Prod. 19:25-40.


J.V. Beck, K.J. Arnold 1977. Parameter Estimation in Engineering and Science. Wiley & Sons, New York, NY, USA.


K. Coleman, D.S. Jenkinson 2008. ROTHC-26.3. A model for the turnover of carbon in soil. Model description and windows users’ guide. Rothamsted Research, Harpenden, Herts, UK.


M.D. Curt, P. Aguado, M. Sanz, G. Sànchez, J. Fernàndez 2006. Clone precocity and the use of Helianthus tuberosus L. stems for bioethanol. Ind. Crops Prod. 24:314-320.


F. Danuso 2001. Le colture per la produzione di inulina: attualità e prospettive. Inform. Agr. 3:176-187.


F. Danuso 2003. SEMoLa: uno strumento per la modellizzazione degli agroecosistemi. Pages 283-284 in Atti XXXV Congresso Società Italiana di Agronomia, Napoli, Italy.


F. Danuso, D. Franz, L. Bigot, G. Budoi 1999. CSS: a modular software for cropping system simulation. pp 287-288 in Proc. Int. Symp. Modelling cropping systems, ESA, Lleida, Spain.


F. Danuso, A. Rocca, V. Andreoni, E. Bulfoni 2010. Simulation of the agroenergy farm with the X-farm model: calibration of the crop module for sorghum yield. Ital. J. Agron. 5:275-284.


F. Danuso, A. Rocca, E. Bulfoni, F. Zuliani, L. Barbanti, A. Monti, E. Bonari, C. Tozzini, S. Cosentino, V. Copani 2009. Studio del Comportamento del Sorgo da Biomassa in Diversi Scenari Colturali Italiani con il Modello CSS. Pages 231-232 in Atti XXXVIII Congresso Società Italiana di Agronomia, Firenze, Italia.


F. Danuso, M. Sandra 2006. Climatica: un software per la gestione ed elaborazione delle informazioni climatiche. Riv. Ital. Agrometeorol. 2:26-33.


M.G. D’Egidio, C. Cecchini, T. Cervigni, B. Donini, V. Pignatelli 1998. Production of fructose from cereal stems and polyannual cultures of Jerusalem artichoke. Ind. Crops Prod. 7:113-119.


P. Denoroy 1996. The crop physiology of Helianthus tuberosus L.: a model orientated view. Biomass Bioenerg. 11:11-32.


N.R. Draper, H. Smith 1981. Applied Regression Analysis. 2nd ed. Wiley & Sons, New York, NY, USA.


P.M. Driessen, N.T. Konijn 1992. Land-use systems analysis. Department of Soil Science & Geology, Wageningen Agricultural University Ed., Waningen, The Netherlands.


J.M. Encinar, J.F. Gonzalez, G. Martinez, S. Roman 2009. Jerusalem artichoke pyrolysis: energetic evaluation. J. Anal. Appl. Pyrolysis 85:294-300.


A. Fuchs 1993. Inulin and inulin containing crops. Studies in Plant Science, Vol. 3. Elsevier, Amsterdam, The Netherlands.


P.H.M. Janssen, P.S.C. Heuberger 1995. Calibration of process oriented models. Ecol. Modelling 83:55-66.


A. Lehtomäki, T. Viinikainen, J. Rintala 2008. Screening boreal energy crops and crop residues for methane biofuel production. Biomass Bioenerg. 32:541-550.


R. Paolini, S. Del Puglia, V. Abbate, V. Copani, F. Danuso, G. De Mastro, N. Losavio, V. Marzi, P. Molfetta, V. Pignatelli, G. Venezia, A.V. Vonella 1996. Produttività del topinambur (Helianthus tuberosus L.) in relazione a fattori agronomici diversi. Agricoltura Ricerca 163:126-144.


H.H. van Laar, J. Goudriaan, H. van Keulen 1997. SUCROS97: Simulation of crop growth for potential and water-limited production situations. Quantitative approaches in systems analysis No. 14. AB-DLO Ed., Wageningen, The Netherlands.

Figure 1.

General structure of the Cropping System Simulator model.

Figure 2.

Soil layers as described by the Cropping System Simulator (CSS) model. Eta, evapo-transpiration; Dr, depth of the root layer; Ds, soil depth by roots; Dd, depth of the deep layer; WTdepth, water table depth; WTdepthIni, initial water table depth; Wtab/MWC, depth of water table/maximum water content; WTbed, depth of water table bed.

Figure 3.

Ten-day total rainfall, and maximum and minimum temperatures at the experimental sites during 1999 and 2000.

Figure 4.

Leaves, stalks and storage organs biomass accumulation, as simulated by the Cropping System Simulator (CSS), in comparison with the experimental growth data for the three locations and the two irrigation regimes (year 1999). __ stems d.m. weight simulated; ---- leaves d.m. weight simulated; _ _ _ storage organs d.m. weight simulated; * stems d.m. weight measured; Δ leaves d.m. weight measured; • storage organs d.m. weight measured.

Figure 5.

Leaves, stalks and storage organs biomass accumulation, as simulated by the Cropping System Simulator, in comparison with the experimental growth data for the three locations and the two irrigation regimes (year 2000). __ stems d.m. weight simulated; ---- leaves d.m. weight simulated; _ _ _ storage organs d.m. weight simulated; * stems d.m. weight measured; Δ leaves d.m. weight measured; • storage organs d.m. weight measured.

Figure 6.

Total biomass accumulation obtained by simulations, in comparison with data measured for irrigated and rain-fed Jerusalem artichoke, at Udine, Bologna and Bari, for the two years.

Table 1.

Working groups and locations of the experimental sites.

Working groups Locations Altitude
m a.s.l.
Latitude Longitude
University of Udine Udine (UD) 110 46°03’N 13°13’E
University of Bologna Cadriano (BO) 33 44°30’N 11°20’E
University of Bari Policoro (BA) 31 40°20’N 16°70’E

[i] m a.s.l., meters above sea level.

Table 2.

Values of the soil parameters measured in the different sites and adopted to run the Cropping system simulator model.

Name Soil parameter Cadriano
Gravel Gravel volumetric content, % 0 0 15
Sand Sand content, % 37 40 43
Clay Clay content, % 18 23 17
OM Organic matter, % 1.3 3.6 2.9
CaCO3 Total carbonates, % 1 6 3
Ds Max exploitable soil depth by roots, mm 3000 2500 500
MWC Max water content fine fraction, mm/mm 0.5 0.58 0.56
FC Field capacity, mm/mm 0.26 0.38 0.31
WP Wilting point, mm/mm 0.12 0.14 0.15
Dw Soil working depth, mm 400 400 400
Wtbed Depth of the water table bed, mm 20,000 20,000 60,000
Wtdepth Initial depth of the water table, mm 2500 2000 40,000
Table 3.

General crop parameters of the Cropping system simulator.

Parameter Description Value Estimation
AEffCoef Potential assimilation efficiency coefficient, gCO2/MJ 8 L°
AlloWS Allowed water stress 0.2 L°
CNcrit Critical Nitrogen concentration in plant, kgN/t 0.65 L°
CoefEst Light extinction coefficient 0.8 C
ConvLeCo Conversion of CH2O to leaf biomass, t/tCH2O 0.59 L°
ConvRoCo Conversion of CH2O to root biomass, t/tCH2O 0.71 L°
ConvSmCo Conversion of CH2O to stem biomass, t/tCH2O 0.75 L°
ConvSrCo Conversion of CH2O to storage organs biomass, t/tCH2O 0.71 L°
CRnc Nitrogen content in crop residue, kgN/kg 0.004 L°
Dsow Sowing depth, mm 90 A
FracResDPM Fraction of DPM in crop residues 0.6 L#
LAIc Critical leaf area index for stress 3.5 C
MDRWat Maximum leaf death rate for water stress 0.03 L/C
Qsow Seed sowing amount, t ha/d 0.7 A
Rmax Maximum root deepening rate, mm/°C/d 1.3 L#
Sls Specific leaf surface, haLeaf/t 2.5 L°
Taef1 Minimum temperature for CO2 assimilation, °C 5 L/C
Taef2 Temperature for 25% opt. CO2 assimilation, °C 10 L/C
Taef3 Minimum optimal temperature for CO2 assimilation, °C 15 L/C
Taef4 Maximum optimal temperature for CO2 assimilation, °C 38 L/C
Taef5 Maximum temperature for CO2 assimilation, °C 45 L/C
Tbase Base temperature for development, °C 9 L/C
Ync Nitrogen content in yield, kgN/kg 0.005 A

[i] The values adopted and calibrated for the simulation of the Jerusalem artichoke crop are reported; the estimation method is indicated: from literature (L), calibrated against experimental results of crop biomass (C) and from common agricultural knowledge (A);

[ii] °from SUCROS (Simple and Universal CROp growth Simulator);

[iii] #from a previous study of CSS calibration (Danuso et al., 1999); DPM, decomposable plant material.

Table 4.

Crop parameters of the Cropping System Simulator related to the development stage of the crop.

Parameter Symbol Emergence Complete soil covering End vegetative phase Start accumulation Physiological maturity Harvest maturity Estimation method
Growing degree days GDD 90 490 900 1200 2400 3000 C
Crop coefficients for water use Kc 0.35 0.6 1 1 0.35 - L°
Optimal Nitrogen concentration No 20 20 12 10 8 - L/C
Growth fraction leaves/shoots LF 0.6 0.4 0.3 0.15 0.1 - C
Growth fraction stem/shoots SF 0.4 0.6 0.7 0.8 0.4 - C
Total growth fraction roots RF 0.05 0.05 0.01 0.01 0 - C

[i] For each parameter, values have to be inserted, corresponding to the following stages: plant emergence from soil, complete soil covering, end of vegetative phase, start of accumulation in storage organs, physio logical maturity and harvest maturity. During simulation, the current values for these parameters are obtained by linear interpolation among these values and as a function of current value of accumulated GDD; the estimation method is indicated: from the literature (L), calibrated against experimental results of crop biomass (C);

[ii] °from FAO (Allen et al., 1998).

Table 5.

Main cropping management techniques adopted for Jerusalem artichoke.

Bari, 1999 Bari, 2000 Bologna, 1999 Bologna, 2000 Udine, 1999 Udine, 2000
Soil tillage Ploughing at 40 cm, harrowing (2) Ploughing at 40 cm, harrowing (2) Ploughing at 30 cm, harrowing, pre-planting puckering Ploughing at 30 cm, harrowing, pre-planting puckering Ploughing at 40 cm, harrowing (2) Ploughing at 40 cm, harrowing (2)
kg/ha N,
P2O5, K2O
Pre-sowing + post-emergence 150-120-0 Pre-sowing + post-emergence 150-120-0 Pre-sowing + post-emergence 100-100-0 Pre-sowing + post-emergence 100-100-0 Pre-sowing 80-200-200 Pre-sowing 80-0-125
Sowing time 12/4 27/03 12/04 14/04 25/03 21/03
Table 6.

Natural and artificial water supply during the crop cycle for the different cropping scenarios.

Environment Year Thesis Number of irrigations Irrigation water, mm Rainfall, mm° Total water received by crops, mm Groundwater availability
Udine 1999 Rain-fed 6 52 1145 1197 No
Irrigated 12 110 1145 1255 No
2000 Rain-fed 2 70 1104 1174 No
Irrigated 6 210 1104 1314 No
Bologna 1999 Rain-fed 0 0 516 516 Yes
Irrigated 12 318 516 834 Yes
2000 Rain-fed 0 0 410 410 Yes
Irrigated 11 361 410 771 Yes
Bari 1999 Rain-fed 14 143 237 380 Yes
Irrigated 14 570 237 807 Yes
2000 Rain-fed 11 115 403 518 Yes
Irrigated 11 451 403 854 Yes

[i] °Rainfall period: March-October.

Table 7.

Model calibration statistics of comparison of simulated and observed data for irrigated and rain-fed Jerusalem artichoke at Udine, Bologna and Bari, for the two years (total biomass).

Environment Year Thesis RMSE,
Udine 1999 Irrigated 4.88 0.52 0.25 0.06 10.87
Rain-fed 3.43 0.30 1.02 -0.22 5.66
2000 Irrigated 3.40 0.79 0.35 0.00 6.30
Rain-fed 1.54 0.87 0.69 -0.07 2.97
Bologna 1999 Irrigated 3.62 0.87 0.76 -0.13 6.63
Rain-fed 4.18 0.74 0.71 -0.18 7.81
2000 Irrigated 3.43 0.90 0.60 -0.05 5.97
Rain-fed 4.46 0.76 0.35 0.01 6.69
Bari 1999 Irrigated 9.35 0.48 0.19 0.11 7.49
Rain-fed 3.76 0.75 0.38 0.00 7.17
2000 Irrigated 5.25 0.45 0.66 -0.24 7.45
Rain-fed 4.92 0.13 0.55 -0.22 7.71

[i] RMSE, root mean square error; EF, modelling efficiency; CD, coefficient of determination; CRM, coefficient residual mass; MaxAE, maximum absolute error.

Abstract views:


Article Metrics

Metrics Loading ...

Metrics powered by PLOS ALM

Copyright (c) 2011 Mario Baldini, Francesco Danuso, Alvaro Rocca, Elena Bulfoni, Maria Teresa Amaducci, Andrea Monti, Giuseppe De Mastro

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
© PAGEPress 2008-2018     -     PAGEPress is a registered trademark property of PAGEPress srl, Italy.     -     VAT: IT02125780185