Next Article in Journal
Effects of Elevated CO2 on Wheat Yield: Non-Linear Response and Relation to Site Productivity
Previous Article in Journal
Uncertainty of CERES-Maize Calibration under Different Irrigation Strategies Using PEST Optimization Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensitivity Analysis of Plant- and Cultivar-Specific Parameters of APSIM-Sugar Model: Variation between Climates and Management Conditions

1
United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24 Korimoto, Kagoshima-shi, Kagoshima 890-0065, Japan
2
Faculty of Agriculture, Rajarata University of Sri Lanka, Puliyankulama, Anuradhapura 50000, Sri Lanka
3
Faculty of Agriculture, University of the Ryukyus, 1 Senbaru, Nishihara-cho, Okinawa 903-0213, Japan
4
Faculty of Agriculture, Kagoshima University, 1-21-24 Korimoto, Kagoshima-shi, Kagoshima 890-8580, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 18 April 2019 / Revised: 7 May 2019 / Accepted: 11 May 2019 / Published: 14 May 2019

Abstract

:
With increasing demand for food and energy, there is a great need for improving sugarcane productivity. New cultivars and management strategies can be assessed using process-based crop models. Information on cultivars needs to be updated frequently, but it is still limited in most crop models. Therefore, it is important to identify possible candidates for varietal parameterization and calibration. Because sensitivity analysis is computationally expensive, we used a less expensive emulator-based approach to conduct a global sensitivity analysis using the apsimr package and GEM-SA software. We studied the sensitivity of four yield outputs of the APSIM-Sugar model to 13 parameters in rainfed and irrigated conditions in Japan and Sri Lanka. Unlike previous studies, our aim was to give comprehensive insights into the variation in sensitivity due to variation in climate. The results confirmed distinct variation of parameter influence between climates and between management conditions. We identify possible candidates for parameterization and calibration of new cultivars for APSIM-Sugar under different environments, and show the effect of variation in climate on variation in parameter influence under different management conditions. It was confirmed that both radiation use efficiency and transpiration efficiency were sensitive and have to be examined to use new cultivars, though these are not listed as cultivar parameters.

1. Introduction

The production of sugarcane is increasing in importance as both a food and a source of energy. Thus, productivity needs to improve continuously. The frequent introduction of new cultivars [1,2] and management strategies [3,4,5] around the world necessitates new studies to validate them.
Several process-based crop models have been developed to simulate sugarcane growth and yield. The widely cited Agricultural Production Systems Simulator (APSIM)-Sugar model [6,7] and the Decision Support System for Agrotechnology Transfer (DSSAT)-Canegro model [8] can simulate the growth and yield of sugarcane under different environmental (climatic and soil) and management (irrigation, fertilization, etc.) conditions with different cultivars [2]. Crop models have been used for decision support in sugarcane farming, including irrigation scheduling [9,10], harvest scheduling [1,11], and fertilizer management [12]. They are also used in projecting the influence of climate change [13,14,15] and guiding varietal improvement [16]. However, most crop models are limited to old cultivars, so studies of varietal differences are limited [17]. Further, model based assessments of currently grown popular cultivars are still rare, for instance, no currently grown commercial cultivars are listed in APSIM-Sugar 7.10 [2]. Process-based crop models should be well parameterized and calibrated to achieve high accuracy. However, the measurement of a wide range of parameters is practically difficult, and not enough data are available from field experiments (including breeding trials) to parameterize all the required varietal information in APSIM-Sugar. Therefore, it is important to identify the most influential parameters in the simulation of outputs through sensitivity analysis (SA). Some parameters are easily measurable and available, whereas some are not, but they can be estimated through varietal calibration. SA can guide crop modelers in parameterizing their cultivars using available data and identify important parameters that need to be estimated by varietal calibration.
SA techniques can be broadly categorized as local or global [18]; local SA considers a single parameter at a time, and global SA considers the combined effect of multiple parameters. Saltelli and Annoni [19] reported the advantages of global SA over local SA; several global SA methods are available to estimate the sensitivity of process-based crop models to parameters, but they are computationally expensive. To minimize the computational cost, Sexton et al. [2] used an emulator-based approach to study the sensitivity of APSIM-Sugar to cultivar parameters. An emulator is a simplified statistical approximation of a more complex model [20] used in place of computationally expensive models. An emulator with a high enough accuracy can replace an actual simulator to perform SA [21]. In this approach, the simulator initially runs for relatively few simulations to build the emulator, then the emulator is used for the SA.
Here, we aimed at assessing the sensitivity of four yield outputs—total aboveground biomass, fresh cane weight, the weight of plant sucrose, and commercial cane sugar—to variations in 13 parameters used in APSIM-Sugar model under different environment and management conditions using emulator-based global sensitivity analysis. By assessing variations in the influence of parameters, we aimed at identifying important candidates for the parameterization and calibration of APSIM-Sugar in different environments and management conditions. We also investigated the effect of radiation use efficiency (RUE) and transpiration efficiency (TE) on the parameterization and calibration of APSIM-Sugar. The relationship between the variations in influence and climate is examined.

2. Materials and Methods

2.1. Study Area

Two locations with different climates were selected for this study: Itoman city (26°7′58″ N, 127°40′52″ E), Okinawa prefecture, Japan, and Sevanagala town (6°22′13″ N, 80°54′47″ E), Monaragala district, Sri Lanka. In Japan, Okinawa is the foremost cane sugar producer, accounting for about 59% of the country’s production. In the Köppen climate classification, Okinawa is classified as Cfa (humid subtropical) [22]. In Sri Lanka, sugar is an important subsector in the economy, with great potential for employment and income generation and for the development of the country’s dry zone. Sugarcane is grown mainly in the southern dry zone, where most processing plants are located. This zone of Sri Lanka is classified as As (tropical with dry summer) [23]. In Okinawa, the 30-year average monthly maximum and minimum temperatures both differ by about 15°C between summer and winter (Figure 1a). In Sri Lanka, they differ by only 2 to 3 °C (Figure 1b). Solar radiation shows high seasonal variation in Okinawa (15 MJ/m2/day), but less variation (6 MJ/m2/day) in Sri Lanka. Average monthly rainfall shows a contrasting unimodal rainfall pattern in Okinawa but a bimodal pattern in Sri Lanka. The average monthly rainfall is 212 mm (± 72 mm SD) in Itoman and 169 mm (± 67) in Sri Lanka. APSIM uses the Priestley–Taylor method to estimate potential evapotranspiration (ET). The monthly potential ET varies slightly in Sri Lanka (2 mm/day) and moderately in Okinawa (4 mm/day; Figure 1).

2.2. Apsim Simulation

APSIM is a process-based dynamic crop model that combines biophysical and management modules within a central engine to simulate diverse cropping systems [6,24]. The model is driven by daily climate data and can simulate the growth, development, and yield of crops and their interactions with soil.
APSIM-Sugar model simulates sugarcane growth via dry weight accumulation due to intercepted radiation in a daily time step. Dry weight accumulation in APSIM-Sugar is determined by RUE [7]. The model partitions the daily accumulated biomass into leaf, immature stem top, structural stem, roots, and sucrose. Then it simulates the key outputs (fresh cane yield, sugar yield, and sucrose contents) [2,7]. This process is controlled by environmental (soil and climate), plant or ratoon, and cultivar-specific parameters [7,25,26].
Sugarcane growth was simulated from 1 January 2000 to 31 December 2010, using rainfed conditions and irrigated conditions (assuming 50% management-allowed deficit). Soil data for Sri Lanka were derived using pedotransfer functions developed by Gunarathna et al. [27], and other data were gathered from a report by the Soil Science Society of Sri Lanka [28]. Soil data for Okinawa were collected through comprehensive soil analysis. Meteorological data for Sri Lanka were extracted from the AgMERRA global gridded climate dataset [29] by the NetCDF-Extractor v. 2.0 tool of AgriMetSoft (https://www.agrimetsoft.com). Those for Okinawa were obtained from the Japan Meteorological Agency website (http://www.data.jma.go.jp/gmd/risk/obsdl/index.php). In both locations, the data on daily rainfall, maximum temperature, minimum temperature, and solar radiation covered 1980 to 2010. Figure 1 shows the variation of climatic variables used for these simulations. Table 1 summarizes the soil and management conditions of the two locations.

2.3. Gaussian Emulation Machine for Sensitivity Analysis (GEM-SA)

GEM-SA is an open-source software used to build software emulators from a set of inputs and outputs so as to perform predictions, uncertainty analysis, and SA using far fewer code runs than Monte Carlo–based methods [30]. It uses Bayesian analysis of computer code outputs [31]. The underlying mathematical procedures used in GEM-SA and analytical procedures used to conduct SA with GEM-SA are described in detail by Kennedy and O’Hagan [31], Kennedy et al. [32] and Kennedy and Petropoulos [30]. GEM-SA estimates two variance-based sensitivity indices—the main effect and the total effect—by partitioning the total output variance induced by variations in all input parameters [33]. Gunarathna et al. [34] reported the accuracy of emulators developed by GEM-SA for subtropical environments. Sexton et al. [2] used GEM-SA to assess the sensitivity of sugarcane biomass and yield to ten parameters in two regions in Australia; despite little variation in climate between the regions, sensitivities differed between the regions, as one region grows rainfed sugarcane and the other grows irrigated sugarcane. Gunarathna et al. [35] used GEM-SA to assess the sensitivity of outputs of APSIM-Oryza to soil parameters in different climatic conditions and reported different sensitivities among regions.
We used GEM-SA to assess the sensitivity of total (green + trash) aboveground biomass, fresh cane weight (canefw), the weight of plant sucrose (sucrose_wt), and commercial cane sugar (% ccs) to 13 selected parameters (Table 2). We assessed emulator accuracy, parameters that influence outputs, the variability of those parameters between years and between rainfed and irrigated conditions, and the effect of climate. Figure 2 shows an overview of the procedure we used for this study.

2.4. Global Sensitivity Analysis

SA indicates which input parameters have the most influence on model outputs. We conducted SA using daily climate data from 2000 to 2010 to assess the sensitivity of the 4 yield outputs to 13 parameters (11 cultivar-specific parameters, RUE, and TE; Table 2) under rainfed and irrigated conditions in two distinct environments (Figure 2). Initially, 300 test points (of each and every parameter and outputs) evenly distributed between lower and upper bounds (Table 2) and related outputs of APSIM were generated by the apsimr package [36] of R software [37]. We chose the upper and lower bounds in consideration of the range of parameter values of cultivars in APSIM 7.10 Sugar model. We used the Gaussian Process emulator in GEM-SA [32] to develop 160 emulators (10 years × 4 outputs × 2 environments × 2 management conditions). Variance based sensitivity indices (Main, Si and total effects, STi) were estimated by partitioning the total output variance induced by variations in all input parameters with the assumption that all input uncertainties are unknown but uniform. The main effect index (Si) is defined as:
S i = V a r { E ( f ( X | x i ) } V a r { f ( X ) }
where, Var{f(X)} is the total variance in the output given variations in all parameters; V a r { E ( f ( X | x i ) } is the variance in the expected output f(X) given xi is known. Hence, Si represents the expected reduction in output variance if parameter xi were known [2]. The relative importance of each parameter in terms of its effect on output uncertainty can be ranked using this Si values of selected parameters [38]. The total sensitivity index (STi) is defined as:
S T i = 1 V a r { E ( f ( X | x i ) } V a r { f ( X ) }
where, V a r { E ( f ( X | x i ) } is the variance in the expected output f(X) if all parameters except xi are known. Saltelli and Annoni [39] suggested to use STi when sensitivity analysis aims to set non-influential parameters to default values and removing them from potential calibrations.
The prior mean option for each input was set as linear. Models were assessed by using the leave-one-out cross-validation procedure of GEM-SA. In cross-validation procedure, a series of left-out points were estimated as output from the code which we actually know the true values. Therefore, the error values are readily available [40]. GEM-SA calculates the cross-validation root-mean-squared error (RMSE, Equation (3)) and root-mean-squared standardized error (RMSSE, Equation (4)) from the results of the cross-validation [30], and the sigma squared (σ2) value. We used these inbuilt diagnostics to assess the accuracy of the emulator approximations. We evaluated the variation in sensitivity of model outputs to changes in parameter values under two management conditions and two environments using the variances (as a percentage) of the main effect index (Si) and total effect index (STi) provided by GEM-SA.
R M S E = i = 1 n ( y i y ^ ) 2 n
R M S S E = i = 1 n ( ( y i y ^ ) / s i ) 2 n
where, yi is the true output for the ith training run, y ^ is the corresponding emulator approximation, si is the standard deviation calculated with the ith training point removed and n is the number of runs [30].

3. Results and Discussion

3.1. Emulator Accuracy

A series of internally calculated statistical measures (σ2, RMSE and RMSSE) of GEM-SA were used to quantify the uncertainty of the sensitivity analysis due to emulation of model simulation [30]. We used σ2 (variance of the emulator after standardization of output) to evaluate the linearity of the GEM-SA emulators, as σ2 ranges near to 0 when a model shows linearity and show higher values when a model shows moderate to high nonlinearity, albeit with no defined cutoff values. Petropoulos et al. [41] reported that their emulators showed linearity or moderate nonlinearity at σ2 values between 0.13 and 1.6. We got σ2 values of 0.15 to 1.43 for Okinawa and 0.10 to 0.59 for Sri Lanka (Figure 3), so the models showed good to moderate linearity in both environments, and higher linearity in Sri Lankan conditions. Hence emulators can successfully replace the simulators.
RMSSE values are near 1 when emulator results are close to simulator results; lower and higher values respectively indicate over- and underestimation. Our RMSSE values of cross-validation results were close to 1: 0.89 to 1.09 in Okinawa and 0.85 to 1.12 in Sri Lanka (Figure 3). These values are lower than previously reported values of emulators assessed as satisfactory [30,41,42].

3.2. Sensitivity of Crop Growth and Yield to Changes in Plant and Cultivar Specific Parameters

Both main effect and total effect sensitivity indices of parameters varied significantly across environments and management conditions (Figure 4, Figure 5, Figure 6 and Figure 7) and among years. There were no big differences between the main and total effect indices, and the main effect index was able to explain a major portion of the variability, so we neglected combined effects. Under rainfed conditions, soil moisture deficit is possible and may affect other processes. Thus, parameters showed less interannual variation under rainfed conditions than under irrigated conditions. Interannual variability was greater in Okinawa than in Sri Lanka as a result of wide variation in temperature and solar radiation during the growing season. We ranked the influence of parameters according to the median value of main effects over the 10 study years (Table 3). Significant contributions (over 1%) is highlighted (bold) in the Table 3.
RUE of growth stage 4 (from beginning of cane growth to flowering; RUE4) was the most influential on all the four tested model outputs despite locations and management (Table 3 and Figure 4, Figure 5, Figure 6 and Figure 7). In both environments, biomass and cane fresh weight showed a different pattern of parameter influence from sucrose weight and ccs for all other parameters. RUE of growth stage 3 (from emergence to beginning of cane growth; RUE3) was the second most influential parameter for biomass and cane fresh weight. It was ranked two to three for CCS and sucrose weight in Okinawa and two to six for CCS and sucrose weight in Sri Lanka. For CCS and sucrose weight, MSS has a higher effect than RUE3. Green leaf number (GLN) was the third most influential parameter for biomass, and cane fresh weight, however, it was ranked four to five in Okinawa and three to five in Sri Lanka for CCS and sucrose weight. Cane fraction (CF) and thermal time between emergence and beginning of cane (EB) also influenced biomass in both environments in both water regimes.
Transpiration efficiency of growth stage 4 (TE4) has a greater influence in biomass under rainfed conditions than under irrigation conditions [43]. This was more pronounced under Sri Lankan growing conditions, probably because higher water-holding capacity of the Sri Lankan soil (Table 1). In both environments, minimum structural stem sucrose content (MSS) and sucrose fraction under stress (SF) influenced CCS and sucrose weight under both water regimes. CF also influenced CCS in both climates and sucrose weight in Sri Lanka under rainfed conditions. MSS reduction (MSSR) influenced ccs under rainfed conditions in Okinawa but not in Sri Lanka. Sexton et al. [2] reported higher sensitivity of biomass production to RUE and GLN and of sucrose yield to RUE and SF under Australian conditions.
In almost all cases, RUE4, RUE3, TE4 (plant-specific parameters), GLN, CF, MSS, and SF (cultivar-specific parameters) explained >90% of variability (Table 3). RUE4 had the greatest influence, and was more influential in Sri Lanka than in Okinawa in both water regimes. In both climates, growth stage 4 fell between July and the following February or March, when solar radiation, potential ET, and maximum and minimum temperatures varied much more in Okinawa than in Sri Lanka.

3.3. Role of RUE4 on APSIM-Sugar Simulations

Daily dry matter production (DDMP; g/m2) of sugarcane under irrigated conditions was simulated using the default values for RUE (1.8 g/MJ) and TE (8.7 g/kPa/kg) in APSIM-Sugar. Maximum possible daily dry matter production (MDMP; g/m2) was calculated as the product of daily intercepted solar radiation and RUE. APSIM estimates intercepted solar radiation as a function of leaf area index and radiation extinction coefficient of 0.38 (Equation (5)).
I = I 0 × E X P ( k × L A I )
where, I is the Intercepted solar radiation, I0 is the total solar radiation at top of the canopy, k is extinction coefficient and LAI is the leaf area index.
During growth stage 4, DDMP had a close relationship with solar radiation in both environments (Figure 8), confirming the high dependency of DDMP on solar radiation. However, it had greater uncertainty throughout the range of intercepted solar radiation in Okinawa than in Sri Lanka, indicating limitation by other factors. This explains the higher sensitivity of DDMP to RUE4 in Sri Lanka. RUE is highly sensitive to nitrogen stress and to high and low temperatures [44]. No nitrogen stress was reported during this period in either environment (see the next paragraph). Therefore, temperature and solar radiation contributed most to the variation in the sensitivity of DDMP to RUE4 in both environments.
Lower DDMP values than MDMP values confirm that APSIM-Sugar did not use the maximum RUE values to simulate DDMP during growth stage 4 specially in days with higher intercepted solar radiation (Figure 9c). During the study period, the minimum and maximum temperatures were within the optimum range (15–45 °C, Figure 9a) [45], and no nitrogen deficit and only slight water stress were recorded (Figure 9b). There was no water stress on days with a radiation level of <15 MJ/m2 (Figure 9d), and APSIM-Sugar operates with maximum RUE. RUE can be maximized with favorable water, nitrogen, and temperature conditions [44], and some authors reported higher RUE values than 2 g/MJ [46,47,48] while APSIM-Sugar remains at 1.8 g/MJ. This might cause APSIM to underestimate yield, especially in simulation studies based on modern commercial sugarcane cultivars.
Since RUE is a standard parameter in plant and ratoon crops, it is usually unchanged in varietal parameterization. However, cultivars show a range of RUE values [48], so it needs to be parameterized to achieve good simulations. Sexton et al. [2,49] suggested to add RUE and TE as varietal parameters in upcoming APSIM versions; our results prompt us to agree. Our results under Okinawan conditions are similar to those reported by Sexton et al. [2], but there are no published results to compare with the Sri Lankan environment. Therefore, we suggest the need for studies to determine the influence of sugarcane cultivar parameters in different tropical environments to confirm our findings. Sexton et al. [49] reported the inability of APSIM-Sugar to differentiate the yields of four commercial cultivars. The lack of enough cultivar parameters that influence the growth and yield of sugarcane means that APSIM-Sugar may not be able to distinguish varietal differences if those parameters are not parameterized and calibrated properly.

3.4. Investigation of Interannual Variation in Parameter Influence

Sexton et al. [2] found high interannual variation of parameter influence and suggested the contribution of climatic variation. Several plant- and cultivar-specific parameters showed considerable interannual variation in influence (Figure 4, Figure 5, Figure 6 and Figure 7). To study the influence of climatic parameters on interannual variation in sensitivity in Sri Lanka, we investigated the sensitivity of canefw to highly influential parameters (Table 3). Considering climatic factors, we selected three growing seasons (Y2 2001–02, Y7 2006–07, Y10 2009–10) for comparison (Figure 10). We calculated cumulative growing degree-days, assuming a base temperature of 9 °C. Y10 had the highest cumulative growing degree-days (6663) and Y7 had the lowest (6425). Y7 had the highest cumulative rainfall (2232 mm) and Y2 had the lowest (1504 mm). Y2 had the highest cumulative solar radiation (6703 MJ/m2) and cumulative potential ET (1688 mm). Y10 had the lowest cumulative solar radiation (6179 MJ/m2) and cumulative potential ET (1566 mm).
As a tropical country, Sri Lanka receives high solar radiation even in high-rainfall years. Therefore, under irrigated conditions, canefw had high sensitivity to RUE4 in all three years. Since irrigation supplies all crop water requirements, RUE4 had greater influence in Y2, a low-rainfall year with high solar radiation, than in the other years. Under rainfed conditions, RUE4 had a stronger relationship to canefw in Y7 (highest rainfall) than in the other years. This result confirms the higher sensitivity of canefw to RUE4, but the sensitivity is directly linked to moisture availability for plants. Solar radiation was similar among years in the first 3 months of the growing season but then differed among years. Under irrigated conditions, the sensitivity of canefw to RUE3 was similar among years. Under rainfed conditions, however, it was less in Y2 than in the other years. In Y2, RUE3 became insensitive to increasing RUE as controlled by soil moisture deficit due to less rainfall and high potential ET. Under irrigated conditions, GLN, TE4, CF, and EB showed little or no variation in influence among years. Under rainfed conditions, however, GLN, TE4, and CF showed greater variation in influence, more so in Y2, owing to both higher solar radiation and moisture stress.
To study the influence of climatic parameters on interannual variation in sensitivity in Okinawa, we selected three growing seasons (Y2 2001–2002, Y4 2003–2004, and Y9 2008–2009) for comparison (Figure 11). Y4 had the highest cumulative growing degree-days (4847) and Y2 had the lowest (4833). Since Okinawa receives high rainfall due to typhoons, the cumulative rainfall does not reflect cropping conditions, so we calculated the effective rainfall by using a water balance approach. Y2 had the highest cumulative rainfall (2478 mm) and Y9 had the lowest (1215 mm), but Y4 had the highest cumulative effective rainfall and Y2 had the lowest. Y9 had the highest cumulative solar radiation (4999 MJ/m2) and cumulative potential ET (1247 mm). Y2 had the lowest cumulative solar radiation (4690 MJ/m2) and cumulative potential ET (1170 mm).
With lower solar radiation than in Sri Lanka, the sensitivity of canefw to RUE4 was lower, but still significant, in Okinawa in all three years. The sensitivity of canefw to RUE4 was lowest in Y2 owing to the lower solar radiation and effective rainfall than in the other two years. Under the lower solar radiation, the influence of RUE4 differed little between irrigated and rainfed conditions. This result confirms the higher sensitivity of canefw to RUE4, but the sensitivity is directly linked to solar radiation and moisture availability for plants. RUE3 also showed less influence on canefw in Y2 than in the other two years, on account of lower effective rainfall and solar radiation. Although the climatic conditions were similar, under rainfed conditions canefw was less sensitive to RUE3 in Y4 than in Y9. Canefw was more sensitive to GLN in Y4 and Y9 than in Y2 owing to higher effective rainfall and solar radiation. TE4, CF, and EB had little or no variation in influence among years under either irrigated or rainfed conditions.

3.5. Relationship between Statistical Dispersion and Climatological Parameters

We examined the statistical dispersion of emulator canefw outputs in response to climatic parameters (Figure 12). In Sri Lanka, under irrigated conditions, cumulative solar radiation (CSR) had the closest relationship with the average output, and cumulative growing degree-days had the closest relationship with SD. Under rainfed conditions, cumulative rainfall had the closest relationship with average and CSR had the closest relationship with SD. In Okinawa, under both water regimes, CSR had the closest relationship with average and SD.
Solar radiation is a key factor governing the sensitivity of canefw to parameters irrespective of climatic conditions. Under high solar radiation with abundant soil moisture, canefw depended mainly on solar radiation, and its variability in sensitivity depended on temperature. Under high solar radiation with water stress, canefw depended mainly on rainfall, and its variability in sensitivity was governed mainly by solar radiation. Under low solar radiation, irrespective of water availability, canefw and its sensitivity were governed mainly by solar radiation. Similarly, Grossi et al. [50] reported the higher sensitivity of sorghum yield to rainfall, solar radiation, and CO2 in DSSAT simulations. Hence, solar radiation, rainfall, and temperature have the greatest influence on canefw in crop models.

4. Conclusions

We used GEM-SA to assess the influence of 13 parameters (11 cultivar-specific parameters, RUE, and TE) of APSIM-Sugar on predicted biomass, fresh cane yield, sucrose weight, and commercial cane sugar yield (ccs) under rainfed and irrigated conditions in two distinctive environments. In both environments, all four outputs were highly sensitive to the RUE of crop growth stage 4 (from the beginning of cane growth to flowering) and growth stage 3 (from emergence to the beginning of cane growth) and to green leaf number, irrespective of water regime. Biomass was sensitive to cane fraction and thermal time from emergence to the beginning of cane (EB) in both environments. In Okinawa, biomass and fresh cane yield were sensitive to TE of growth stage 4 under rainfed conditions, but less sensitive under irrigated conditions. In Sri Lanka, they were sensitive under both water regimes. In both environments, ccs and sucrose weight were sensitive to minimum structural stem sucrose content (MSS) and sucrose fraction under stress condition (SF) under both water regimes. In Okinawa, ccs and sucrose weight were slightly sensitive to MSS reduction (MSSR) and cane fraction under rainfed but less sensitive under irrigated conditions. In Sri Lanka, biomass, fresh cane yield, and sucrose yield were sensitive to TE of growth stage 4 under rainfed but less sensitive under irrigated conditions. These results confirm distinct variations in parameter influence across climates, management conditions, and outputs. This shows why SA conducted in similar environments is vital to identifying parameters important for parameterization and calibration of sugarcane cultivars. In both environments, green leaf number and cane fraction were important candidates for parameterization of cultivars. We suggest that attention to calibration of EB, MSS, MSSR, and SF will improve the accuracy of simulations of sugarcane growth and yield in both environments. Although they are not listed as cultivar parameters in the APSIM-Sugar model, if reliable and ample data available, it is advisable to calibrate TE of growth stage 4 and RUE of growth stages 3 and 4 also. Interannual variations in solar radiation, rainfall, and temperature explained the variation of parameter influence. Therefore, variations in climatic parameters must be accounted for in the modeling of sugarcane growth and yield using APSIM.

Author Contributions

M.H.J.P.G., K.S. and M.K.N.K conceptualized, conceived and performed, and K.S., T.N. and K.M. supervised the study. M.H.J.P.G. and K.S. interpreted data and developed the manuscript. All authors read and approved the final manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lisson, S.N.; Inman-Bamber, N.G.; Robertson, M.J.; Keating, B.A. The historical and future contribution of crop physiology and modelling research to sugarcane production systems. Field Crop. Res. 2005, 92, 321–335. [Google Scholar] [CrossRef]
  2. Sexton, J.; Everingham, Y.L.; Inman-Bamber, G. A global sensitivity analysis of cultivar trait parameters in a sugarcane growth model for contrasting production environments in Queensland, Australia. Eur. J. Agron. 2017, 88, 96–105. [Google Scholar] [CrossRef]
  3. Gunarathna, M.H.J.P.; Sakai, K.; Nakandakari, T.; Kazuro, M.; Onodera, T.; Kaneshiro, H.; Uehara, H.; Wakasugi, K. Optimized subsurface irrigation system (OPSIS): Beyond traditional subsurface irrigation. Water 2017, 9, 599. [Google Scholar] [CrossRef]
  4. Gunarathna, M.H.J.P.; Sakai, K.; Nakandakari, T.; Momii, K.; Onodera, T.; Kaneshiro, H.; Uehara, H.; Wakasugi, K. Optimized subsurface irrigation system: The future of sugarcane irrigation. Water 2018, 10, 314. [Google Scholar] [CrossRef]
  5. Surendran, U.; Jayakumar, M.; Marimuthu, S. Science of the total environment low cost drip irrigation: Impact on sugarcane yield, water and energy saving in semiarid tropical agro ecosystem in India. Sci. Total Environ. 2016, 573, 1430–1440. [Google Scholar] [CrossRef]
  6. Holzworth, D.P.; Huth, N.I.; deVoil, P.G.; Zurcher, E.J.; Herrmann, N.I.; McLean, G.; Chenu, K.; van Oosterom, E.J.; Snow, V.; Murphy, C.; et al. APSIM—Evolution towards a new generation of agricultural systems simulation. Environ. Model. Softw. 2014, 62, 327–350. [Google Scholar] [CrossRef]
  7. Keating, B.A.; Robertson, M.J.; Muchow, R.C.; Huth, N.I. Modelling sugarcane production systems I. Development and performance of the sugarcane module. Field Crop. Res. 1999, 61, 253–271. [Google Scholar] [CrossRef]
  8. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  9. Everingham, Y.; Baillie, C.; Inman-Bamber, G.; Baillie, J. Forecasting water allocations for Bundaberg sugarcane farmers. Clim. Res. 2008, 36, 231–239. [Google Scholar] [CrossRef] [Green Version]
  10. Inman-Bamber, N.G.; McGlinchey, M.G. Crop coefficients and water-use estimates for sugarcane based on long-term bowen ratio energy balance measurements. Field Crop. Res. 2003, 83, 125–138. [Google Scholar] [CrossRef]
  11. Bocca, F.F.; Rodrigues, L.H.A.; Arraes, N.A.M. When do I want to know and why? Different demands on sugarcane yield predictions. Agric. Syst. 2015, 135, 48–56. [Google Scholar] [CrossRef]
  12. Gaydon, D.S.; Wang, E.; Poulton, P.L.; Ahmad, B.; Ahmed, F.; Akhter, S.; Ali, I.; Amarasingha, R.; Chaki, A.K.; Chen, C.; et al. Evaluation of the APSIM model in cropping systems of Asia. Field Crop. Res. 2017, 204, 52–75. [Google Scholar] [CrossRef]
  13. Singels, A.; Jones, M.; Marin, F.; Ruane, A.; Thorburn, P. Predicting climate change impacts on sugarcane production at sites in Australia, Brazil and South Africa using the canegro model. Sugar Tech. 2014, 16, 347–355. [Google Scholar] [CrossRef]
  14. Thorburn, P.; Biggs, J.; Jones, M.R.; Singels, A.; Marin, F.; Martine, J.-F.; Chinorumba, S.; Viator, R.; Nunez, O. Evaluation of the APSIM-Sugar model for simulation sugarcane yield at sites in seven countries: Initial results. Proc. S. Afr. Sugar Technol. Assoc. 2014, 318–322. [Google Scholar]
  15. Zubair, L.; Nissanka, S.P.; Weerakoon, W.M.W.; Herath, D.I.; Karunaratne, A.S.; Agalawatte, P.; Herath, R.M.; Yahiya, S.Z.; Punyawardhene, B.V.R.; Vishwanathan, J.; et al. Climate change impacts on rice farming systems in Northwestern Sri Lanka. Handb. Clim. Chang. Agroecosyst. 2015, 3, 315–352. [Google Scholar]
  16. Basnayake, J.; Jackson, P.A.; Inman-Bamber, N.G.; Lakshmanan, P. Sugarcane for water-limited environments. Genetic variation in cane yield and sugar content in response to water stress. J. Exp. Bot. 2012, 63, 6023–6033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Inman-Bamber, N.G.; Jackson, P.A.; Stokes, C.J.; Verrall, S.; Lakshmanan, P.; Basnayake, J. Sugarcane for water-limited environments: Enhanced capability of the APSIM sugarcane model for assessing traits for transpiration efficiency and root water supply. Field Crop. Res. 2016, 196, 112–123. [Google Scholar] [CrossRef]
  18. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley & Sons: Hoboken, NJ, USA, 2008; ISBN 9780470059975. [Google Scholar]
  19. Saltelli, A.; Annoni, P.; Azzini, I.; Campolongo, F.; Ratto, M.; Tarantola, S. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index. Comput. Phys. Commun. 2010, 181, 259–270. [Google Scholar] [CrossRef]
  20. O’Hagan, A. Bayesian analysis of computer code outputs: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 1290–1300. [Google Scholar] [CrossRef] [Green Version]
  21. Uusitalo, L.; Lehikoinen, A.; Helle, I.; Myrberg, K. An overview of methods to evaluate uncertainty of deterministic models in decision support. Environ. Model. Softw. 2015, 63, 24–31. [Google Scholar] [CrossRef] [Green Version]
  22. Gima, H.; Yoshitake, T. A comparative study of energy security in okinawa prefecture and the state of hawaii. Evergr. Jt. J. Nov. Carbon Resour. Sci. Green Asia Strateg. 2016, 3, 36–44. [Google Scholar] [CrossRef]
  23. Buysse, M. Irrigation water supply for non-irrigation purposes in Uda Walawe: Policy and effects. Int. Water Manag. Inst. Rep. 2002, 1–20. [Google Scholar]
  24. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.J.; Holzworth, D.; Huth, N.I.; Hargreaves, J.N.G.; Meinke, H.; Hochman, Z.; et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef]
  25. Dias, H.B.; Inman-Bamber, G.; Bermejo, R.; Sentelhas, P.C.; Christodoulou, D. New APSIM-Sugar features and parameters required to account for high sugarcane yields in tropical environments. Field Crop. Res. 2019, 235, 38–53. [Google Scholar] [CrossRef]
  26. Keating, B.A.; Robertson, M.J.; Muchow, R.C.; Huth, N.I. Modelling sugarcane production systems: II: Analysis of system performance and methodology issues. Field Crop. Res. 2000, 68, 1–20. [Google Scholar]
  27. Gunarathna, M.H.J.P.; Sakai, K.; Nakandakari, T.; Momii, K.; Kumari, M.K.N.; Amarasekara, M.G.T.S. Pedotransfer functions to estimate hydraulic properties of tropical Sri Lankan soils. Soil Tillage Res. 2019, 190, 109–119. [Google Scholar] [CrossRef]
  28. Dassanayake, A.R.; De Silva, G.G.R.; Mapa, R.B.; Kumaragamage, D. Benchmark Soils of the Dry Zone of Sri Lanka: Factsheets; Mapa, R.B., Somasiri, S., Dassanayake, A.R., Eds.; Soil Science Society of Sri Lanka: Peradeniya, Sri Lanka, 2010. [Google Scholar]
  29. Ruane, A.C.; Goldberg, R.; Chryssanthacopoulos, J. Climate forcing datasets for agricultural modeling: Merged products for gap-filling and historical climate series estimation. Agric. For. Meteorol. 2015, 200, 233–248. [Google Scholar] [CrossRef]
  30. Kennedy, M.C.; Petropoulos, G.P. GEM-SA: The Gaussian Emulation Machine for Sensitivity Analysis; Elsevier Inc.: Amsterdam, The Netherlands, 2017; ISBN 9780128030318. [Google Scholar]
  31. Kennedy, M.C.; O’Hagan, A. Bayesian calibration of computer models. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2001, 63, 425–464. [Google Scholar] [CrossRef] [Green Version]
  32. Kennedy, M.C.; Anderson, C.W.; Conti, S.; O’Hagan, A. Case studies in Gaussian process modelling of computer codes. Reliab. Eng. Syst. Saf. 2006, 91, 1301–1309. [Google Scholar] [CrossRef]
  33. Bakar, I.A.A.; Kamaruddin, M.K. Uncertain mesoscopic parameters in dry woven fabric with Gaussian Emulation Machine for Sensitivity Analysis (GEM-SA). J. Phys. Conf. Ser. 2018, 1130. [Google Scholar] [CrossRef]
  34. Gunarathna, M.H.J.P.; Sakai, K.; Kumari, M.K.N. Sensitivity Analysis of APSIM model using GEM-SA software: Emulator accuracy assessment. In Proceedings of the ICESD International Conference, Sydney, Australia, 5–6 March 2019; pp. 6–10. [Google Scholar]
  35. Gunarathna, M.H.J.P.; Sakai, K.; Kumari, M.K.N. Can crop modeling sucess with estimated soil hydraulic parameters? In Proceedings of the PAWEES-INWEPF International Conference Nara 2018, Nara, Japan, 20–22 November 2018; pp. 461–470. [Google Scholar]
  36. Stanfill, B. apsimr: Edit, Run and Evaluate APSIM Simulations Easily Using R. R Package Version 1.2. 2015. Available online: https://CRAN.R-project.org/package=apsimr (accessed on 5 July 2018).
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  38. Oakley, J.E.; Hagan, A.O. Probabilistic sensitivity analysis of complex models: A Bayesian approach. J. R. Stat. B 2004, 66, 751–769. [Google Scholar] [CrossRef]
  39. Saltelli, A.; Annoni, P. How to avoid a perfunctory sensitivity analysis. Environ. Model. Softw. 2010, 25, 1508–1517. [Google Scholar] [CrossRef]
  40. Petropoulos, G.P.; Ireland, G.; Griffiths, H.M.; Kennedy, M.C.; Ioannou-Katidis, P.; Kalivas, D.P. Extending the global sensitivity analysis of the SimSphere model in the context of its future exploitation by the scientific community. Water 2015, 7, 2101–2141. [Google Scholar] [CrossRef]
  41. Petropoulos, G.; Wooster, M.J.; Carlson, T.N.; Kennedy, M.C.; Scholze, M. A global Bayesian sensitivity analysis of the 1d SimSphere soil-vegetation-atmospheric transfer (SVAT) model using Gaussian model emulation. Ecol. Modell. 2009, 220, 2427–2440. [Google Scholar] [CrossRef]
  42. Lee, L.A.; Carslaw, K.S.; Pringle, K.J.; Mann, G.W.; Spracklen, D.V. Emulation of a complex global aerosol model to quantify sensitivity to uncertain parameters. Atmos. Chem. Phys. 2011, 11, 12253–12273. [Google Scholar] [CrossRef] [Green Version]
  43. Sexton, J.; Everingham, Y.L. Global sensitivity analysis of key parameters in a process-based sugarcane growth model—A Bayesian approach. In Proceedings of the 7th International Congress on Environmental Modelling and Software, San Diego, CA, USA, 15–19 June 2014. [Google Scholar]
  44. Zhao, G.; Bryan, B.A.; Song, X. Sensitivity and uncertainty analysis of the APSIM-wheat model: Interactions between cultivar, environmental, and management parameters. Ecol. Modell. 2014, 279, 1–11. [Google Scholar] [CrossRef]
  45. Martiné, J.-F.; Siband, P.; Bonhomme, R. Simulation of the maximum yield of sugar cane at different altitudes: Effect of temperature on the conversion of radiation into biomass Simulation of the maximum yield of sugar cane on the conversion of radiation into biomass. Agronomie 1999, 19, 3–12. [Google Scholar] [CrossRef]
  46. Ferreira, R.A., Jr.; De Souza, J.L.; Lyra, G.B.; Escobedo, J.F.; Santos, M.V.C. Energy conversion efficiency in sugarcane under two row spacings in northeast of Brazil. Rev. Bras. Eng. Agrícola e Ambient. 2015, 19, 741–747. [Google Scholar] [CrossRef] [Green Version]
  47. Martin, M.; Acreche, A. Nitrogen-, water- and radiation-use efficiencies affected by sugarcane breeding in Argentina. Plant Breed. 2017, 136, 174–181. [Google Scholar] [CrossRef]
  48. De Silva, A.L.C.; De Costa, W.A.J.M. Growth and radiation use efficiency of sugarcane under irrigated and rain-fed conditions in sri lanka. Sugar Tech. 2012, 14, 247–254. [Google Scholar] [CrossRef]
  49. Sexton, J.; Inman-bamber, N.G.; Everingham, Y.; Basnayake, J.; Lakshmanan, P.; Jackson, P. Detailed trait characterisation is needed for simulation of cultivar responses to water stress. Proc. Aust. Soc. Sugar Cane Technol. 2014, 36, 82–92. [Google Scholar]
  50. Grossi, M.C.; Justino, F.; Rodrigues, R.D.Á.; Andrade, C.L.T. Sensitivity of the sorghum yield to individual changes in climate parameters: Modelling based approach. Bragantia 2015, 74, 341–349. [Google Scholar] [CrossRef]
Figure 1. Average monthly climate data of (a) Itoman, Okinawa, Japan (1980–2010) and (b) Sevanagala, Monaragala, Sri Lanka (1980–2010); rain, mean monthly rainfall (mm); radn, mean daily solar radiation (MJ/m2); maxt, mean daily maximum temperature (°C); mint, mean daily minimum temperature (°C); eo, Potential evapotranspiration (mm/day).
Figure 1. Average monthly climate data of (a) Itoman, Okinawa, Japan (1980–2010) and (b) Sevanagala, Monaragala, Sri Lanka (1980–2010); rain, mean monthly rainfall (mm); radn, mean daily solar radiation (MJ/m2); maxt, mean daily maximum temperature (°C); mint, mean daily minimum temperature (°C); eo, Potential evapotranspiration (mm/day).
Agronomy 09 00242 g001
Figure 2. Flowchart of the analysis procedure used for the study.
Figure 2. Flowchart of the analysis procedure used for the study.
Agronomy 09 00242 g002
Figure 3. (a) Sigma squared value, (b) Cross validated Cross-validation root mean squared standardized error of results of selected outputs under two climatic conditions (Okinawa and Sri Lanka). IR, Irrigated; RF, rainfed.
Figure 3. (a) Sigma squared value, (b) Cross validated Cross-validation root mean squared standardized error of results of selected outputs under two climatic conditions (Okinawa and Sri Lanka). IR, Irrigated; RF, rainfed.
Agronomy 09 00242 g003
Figure 4. Sensitivity of APSIM-Sugar biomass (g/m2) to parameters in (a) Okinawa and (b) Sri Lanka; Si, main effect; STi, total effect; IR, irrigated; RF, rainfed; TLS1–5, tiller leaf size; TE1–6, transpiration efficiency; SF2, sucrose fraction under stress; SD, sucrose delay; RUE3–5, radiation use efficiency; MSS, minimum structural stem sucrose; MSSR, MSS reduction; LS1–3, leaf size; GLN, green leaf number; FC, thermal time from flowering to crop end; EB, thermal time from emergence to beginning of cane; CF, cane fraction; BF, thermal time from beginning of cane to flowering.
Figure 4. Sensitivity of APSIM-Sugar biomass (g/m2) to parameters in (a) Okinawa and (b) Sri Lanka; Si, main effect; STi, total effect; IR, irrigated; RF, rainfed; TLS1–5, tiller leaf size; TE1–6, transpiration efficiency; SF2, sucrose fraction under stress; SD, sucrose delay; RUE3–5, radiation use efficiency; MSS, minimum structural stem sucrose; MSSR, MSS reduction; LS1–3, leaf size; GLN, green leaf number; FC, thermal time from flowering to crop end; EB, thermal time from emergence to beginning of cane; CF, cane fraction; BF, thermal time from beginning of cane to flowering.
Agronomy 09 00242 g004
Figure 5. Sensitivity of APSIM-Sugar fresh cane yield (t/ha) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Figure 5. Sensitivity of APSIM-Sugar fresh cane yield (t/ha) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Agronomy 09 00242 g005
Figure 6. Sensitivity of APSIM-Sugar commercial cane sugar (%) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Figure 6. Sensitivity of APSIM-Sugar commercial cane sugar (%) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Agronomy 09 00242 g006
Figure 7. Sensitivity of APSIM-Sugar sucrose yield (g/m2) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Figure 7. Sensitivity of APSIM-Sugar sucrose yield (g/m2) to parameters in (a) Okinawa and (b) Sri Lanka. Si, main effect; STi, total effect; IR, irrigated; RF, rainfed. For parameter names, see Figure 4.
Agronomy 09 00242 g007
Figure 8. Relationship of intercepted solar radiation (MJ/m2) during growth stage 4 with simulated daily dry matter production (g/m2) of irrigated sugarcane and calculated maximum possible daily dry matter production (g/m2) of sugarcane in (a) Okinawa and (b) Sri Lanka.
Figure 8. Relationship of intercepted solar radiation (MJ/m2) during growth stage 4 with simulated daily dry matter production (g/m2) of irrigated sugarcane and calculated maximum possible daily dry matter production (g/m2) of sugarcane in (a) Okinawa and (b) Sri Lanka.
Agronomy 09 00242 g008
Figure 9. (a) Maximum and minimum temperature variation during crop stage 4 in Sri Lanka. (b) Soil water deficit (swdef_photo) and nitrogen deficit (nfact_photo) for photosynthesis of irrigated sugarcane (1 = no stress, 0 = full stress). (c) Relationship of intercepted solar radiation (MJ/m2) with simulated daily dry matter production (DDMP; g/m2) of irrigated sugarcane and calculated maximum possible daily dry matter production (MDMP; g/m2) of sugarcane. (d) Relationship between solar radiation (MJ/m2) and soil water deficit for photosynthesis of irrigated sugarcane.
Figure 9. (a) Maximum and minimum temperature variation during crop stage 4 in Sri Lanka. (b) Soil water deficit (swdef_photo) and nitrogen deficit (nfact_photo) for photosynthesis of irrigated sugarcane (1 = no stress, 0 = full stress). (c) Relationship of intercepted solar radiation (MJ/m2) with simulated daily dry matter production (DDMP; g/m2) of irrigated sugarcane and calculated maximum possible daily dry matter production (MDMP; g/m2) of sugarcane. (d) Relationship between solar radiation (MJ/m2) and soil water deficit for photosynthesis of irrigated sugarcane.
Agronomy 09 00242 g009
Figure 10. Depiction of interannual sensitivity variation of sugarcane in Sri Lanka: cumulative rainfall (mm), cumulative solar radiation (MJ/m2), cumulative growing degree-days, cumulative potential evapotranspiration (mm), and parameter sensitivity variation among 2001–2002, 2006–2007, and 2009–2010.
Figure 10. Depiction of interannual sensitivity variation of sugarcane in Sri Lanka: cumulative rainfall (mm), cumulative solar radiation (MJ/m2), cumulative growing degree-days, cumulative potential evapotranspiration (mm), and parameter sensitivity variation among 2001–2002, 2006–2007, and 2009–2010.
Agronomy 09 00242 g010
Figure 11. Depiction of interannual sensitivity variation of sugarcane in Okinawa: cumulative rainfall (mm), cumulative effective rainfall (mm), cumulative solar radiation (MJ/m2), cumulative growing degree-days, cumulative potential evapotranspiration (mm), and parameter sensitivity variation among 2001–2002, 2003–2004, and 2008–2009.
Figure 11. Depiction of interannual sensitivity variation of sugarcane in Okinawa: cumulative rainfall (mm), cumulative effective rainfall (mm), cumulative solar radiation (MJ/m2), cumulative growing degree-days, cumulative potential evapotranspiration (mm), and parameter sensitivity variation among 2001–2002, 2003–2004, and 2008–2009.
Agronomy 09 00242 g011
Figure 12. Relationships between average and SD of emulator canefw predictions with most closely related climatic conditions during the study period.
Figure 12. Relationships between average and SD of emulator canefw predictions with most closely related climatic conditions during the study period.
Agronomy 09 00242 g012
Table 1. Soil and management conditions of selected locations used for the simulations.
Table 1. Soil and management conditions of selected locations used for the simulations.
LocationItoman, Okinawa, JapanSevanagala, Monaragala, Sri Lanka
26°7′58″N 127°40′52″E6°22′13″N 80°54′47″E
SoilShimajiri Mahji
Depth: 110 cm
PAWC: 68.4 mm
Average bulk density 1.11 g/cm3
Solodized Solonetz
Depth: 100 cm
PAWC: 91.6 mm
Average bulk density 1.37 g/cm3
PlantingApril 01 (Spring planting)April 01 (Yala season planting)
Crop duration315 days360 days
Stalk density7 stalks/m28 stalks/m2
Fertilizer190 kg/ha as NH4-N200 kg/ha as Urea
Fertilizer application time31 and 62 days after planting45 and 90 days after planting
IrrigationAutomatic irrigation
Fraction of ASW below which irrigation is applied = 0.5
Efficiency of the irrigation = 0.5
PAWC, Plant available water content; ASW, Available soil water.
Table 2. Selected parameters used to assess the parameter sensitivity on total crop above-ground biomass, fresh cane weight, the weight of plant sucrose and commercial cane sugar.
Table 2. Selected parameters used to assess the parameter sensitivity on total crop above-ground biomass, fresh cane weight, the weight of plant sucrose and commercial cane sugar.
Parameter as Listed in APSIM-Sugar Model (Description)LevelCodeUnitLower and Upper Bound
leaf_size (Leaf area of the respective leaf)Leaf_size_no = 1LS1mm2500–2000
Leaf_size_no = 14LS2mm225,000–70,000
Leaf_size_no = 20LS3mm225,000–70,000
cane_fraction (Fraction of accumulated biomass partitioned to cane)CFgg−10.65–0.80
sucrose_fraction_stalk (Fraction of accumulated biomass partitioned to sucrose)Stress factor = 1SFgg−10.50–0.70
sucrose_delay (Sucrose accumulation delay)SDgm−20–600
min_sstem_sucrose (Minimum stem biomass before partitioning to sucrose commences)MSSgm−2450–1500
min_sstem_sucrose_redn (reduction to minimum stem sucrose under stress)MSSRgm−20–20
tt_emerg_to_begcane (Accumulated thermal time from emergence to beginning of cane)EB°C day1200–1900
tt_begcane_to_flowering (Accumulated thermal time from beginning of cane to flowering)BF°C day5500–6500
tt_flowering_to_crop_end (Accumulated thermal time from flowering to end of the crop)FC°C day1750–2250
green_leaf_no (Maximum number of fully expanded green leaves)GLNNo.9–14
tillerf_leaf_size (Tillering factors according to the leaf numbers)Tiller_leaf_size_no = 1TLS1mm2 mm−21–6
Tiller_leaf_size_no = 4TLS2mm2 mm−21–6
Tiller_leaf_size_no = 10TLS3mm2 mm−21–6
Tiller_leaf_size_no = 16TLS4mm2 mm−21–6
Tiller_leaf_size_no = 26TLS5mm2 mm−21–6
transp_eff (Transpiration efficiency)Stage_code = 1TE1kg kPa/kg0.008–0.014
Stage_code = 2TE2kg kPa/kg0.008–0.014
Stage_code = 3TE3kg kPa/kg0.008–0.014
Stage_code = 4TE4kg kPa/kg0.008–0.014
Stage_code = 5TE5kg kPa/kg0.008–0.014
Stage_code = 6TE6kg kPa/kg0.008–0.014
rue (Radiation use efficiency)Stage_code = 3RUE3g/MJ1.2–2.5
Stage_code = 4RUE4g/MJ1.2–2.5
Stage_code = 5RUE5g/MJ1.2–2.5
Table 3. Parameters most influential in simulation of biomass, canefw, ccs, and sucrose_wt in APSIM-Sugar under different environments (bold indicates parameters that contributed >1% of variation).
Table 3. Parameters most influential in simulation of biomass, canefw, ccs, and sucrose_wt in APSIM-Sugar under different environments (bold indicates parameters that contributed >1% of variation).
ParameterBiomassCane Fresh WeightCCSSucrose Weight
RFIRRFIRRFIRRFIR
SiRankSiRankSiRankSiRankSiRankSiRankSiRankSiRank
Okinawa, Japan
RUE449.0157.7158.3166.3141.2142.3158.5164.61
RUE318.0223.6212.6217.225.437.039.2311.22
GLN7.736.736.235.333.452.754.944.34
CF3.142.540.660.562.361.760.760.66
TE41.750.562.040.550.1100.0110.770.19
EB1.461.851.050.540.280.280.290.28
MSS0.0220.0180.0220.01628.5227.3210.128.83
MSSR0.0170.0250.0210.0251.270.570.380.37
SF0.0140.0160.0160.0195.044.842.251.65
Monaragala, Sri Lanka
RUE460.5171.1167.5176.4144.8154.2170.1179.31
RUE314.5213.0213.1211.725.142.569.427.42
GLN6.736.235.745.033.252.955.234.53
CF5.353.740.750.361.963.041.070.76
TE46.041.466.731.440.680.572.060.18
EB1.561.650.560.750.870.380.480.47
MSS0.0180.0110.0230.01319.5214.823.842.54
MSSR0.0160.0120.0210.0200.390.190.190.09
SF0.0230.080.0260.01810.4310.532.652.35
Biomass, total aboveground biomass (g/m2); Cane fresh weight (t/ha); CCS, commercial cane sugar (%); sucrose weight (g/m2); RF, rainfed; IR, irrigated; Si, main effect; RUE, radiation use efficiency; GLN, green leaf number; CF, cane fraction; TE, transpiration efficiency; EB, thermal time from emergence to beginning of cane; MSS, minimum structural stem sucrose; MSSR, MSS reduction; SF, sucrose fraction under stress.

Share and Cite

MDPI and ACS Style

Gunarathna, M.H.J.P.; Sakai, K.; Nakandakari, T.; Momii, K.; Kumari, M.K.N. Sensitivity Analysis of Plant- and Cultivar-Specific Parameters of APSIM-Sugar Model: Variation between Climates and Management Conditions. Agronomy 2019, 9, 242. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9050242

AMA Style

Gunarathna MHJP, Sakai K, Nakandakari T, Momii K, Kumari MKN. Sensitivity Analysis of Plant- and Cultivar-Specific Parameters of APSIM-Sugar Model: Variation between Climates and Management Conditions. Agronomy. 2019; 9(5):242. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9050242

Chicago/Turabian Style

Gunarathna, M.H.J.P., Kazuhito Sakai, Tamotsu Nakandakari, Kazuro Momii, and M.K.N. Kumari. 2019. "Sensitivity Analysis of Plant- and Cultivar-Specific Parameters of APSIM-Sugar Model: Variation between Climates and Management Conditions" Agronomy 9, no. 5: 242. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy9050242

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop