Next Article in Journal
Phytophthora Antagonism of Endophytic Bacteria Isolated from Roots of Black Pepper (Piper nigrum L.)
Previous Article in Journal
Growth, Photosynthesis, and Physiological Responses of Ornamental Plants to Complementation with Monochromic or Mixed Red-Blue LEDs for Use in Indoor Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Functional Analysis of Pedotransfer Functions Developed for Sri Lankan soils: Applicability for Process-Based Crop Models

by
M. H. J. P. Gunarathna
1,*,
Kazuhito Sakai
2,*,
M. K. N. Kumari
1,3,* and
Manjula Ranagalage
4,5
1
Faculty of Agriculture, Rajarata University of Sri Lanka, Puliyankulama, Anuradhapura 50000, Sri Lanka
2
Faculty of Agriculture, University of the Ryukyus, 1 Senbaru, Nishihara-cho, Okinawa 903-0213, Japan
3
United Graduate School of Agricultural Sciences, Kagoshima University, 1-21-24 Korimoto, Kagoshima-shi, Kagoshima 890-0065, Japan
4
Faculty of Social Science and Humanities, Rajarata University of Sri Lanka, Mihintale 50300, Sri Lanka
5
Faculty of Life and Environmental Sciences, University of Tsukuba, 1-1-1, Tennodai, Tsukuba, Ibaraki 305-8572, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 14 January 2020 / Revised: 14 February 2020 / Accepted: 14 February 2020 / Published: 17 February 2020

Abstract

:
As measurements are expensive and laborious, the estimation of soil hydraulic properties using pedotransfer functions (PTFs) has become popular worldwide. However, the estimation of soil hydraulic properties is not the final aim but an essential input value for other calculations and simulations, mostly in environmental and crop models. This modeling approach is a popular way to assess agricultural and environmental processes. However, it is rarely used in Sri Lanka because soil hydraulic data are rare. We evaluated the functionality of PTFs (developed to estimate field capacity (FC) and the permanent wilting point (PWP) of Sri Lankan soils) for process-based crop models. We used the Agricultural Production Systems sIMulator (APSIM) as the test model. Initially, we confirmed the importance of PWP (LL15) and FC (DUL) by assessing the sensitivity of the soil input parameters on the growth and yield of rice under rainfed conditions. We simulated the growth and yield of rice and the four selected outputs related to the APSIM soil module using the measured and estimated values of FC and PWP. These simulations were conducted for ten years in 16 locations of Sri Lanka, representing wet, intermediate, and dry zones. The simulated total aboveground dry matter and weight of the rough rice, using both input conditions (the measured and PTF-estimated soil hydraulic properties), showed good agreement, with no significant differences between each other. Outputs related to the soil module also showed good agreement, as no significant differences were found between the two input conditions (measured and PTF-estimated soil hydraulic properties). Although the DUL and LL15 are the most influential parameters for the selected outputs of APSIM–Oryza, the estimated FC and PWP values did not change the predictive ability of APSIM. In this way, the functionality of PTFs for APSIM crop modeling is confirmed.

1. Introduction

The direct measurement of some soil hydraulic properties is laborious and expensive. Hence, information about soil hydraulic properties at an appropriate scale is limited in most parts of the world, especially in tropical regions [1]. Pedotransfer functions (PTFs) are predictive functions used to estimate difficult-to-measure soil parameters with easily measurable soil parameters [2]. Point-based PTFs are used to determine soil parameters at specific, conventional values of matric potential [3]. The moisture content at −10 and −33 kPa (representing the field capacity (FC)) and moisture content at –1500 kPa (representing the permanent wilting point (PWP)) are the most common reference values for the soil moisture characteristic curve (SMCC) used for point PTFs [4].
Point PTFs were developed [5] for Sri Lankan soil and confirmed the statistical accuracy of these PTFs for Sri Lankan conditions. However, these PTFs cannot perfectly calculate soil hydraulic parameters. In other words, there are always some errors that could be over or underestimated according to soil type. These errors may be different according to the application of the PTFs. Hence, the final estimation of the model quality should be defined by the model’s applicability [6] because the estimation of hydraulic properties is not the final aim. Instead, these properties are essential inputs for other calculations and simulations [7]. Wosten et al. [7] proposed to use a functional criteria analysis, which evaluates the functional behavior of estimated parameters in different applications. Subsequently, a functional analysis (or functional criteria analysis) was carried out to determine the performance of the estimated parameters on different applications [7,8,9]. However, this type of analysis is rarely used for process-based crop models. A functional analysis, coupled with a stochastic model of the soil hydraulic parameters, can be effectively employed to assess the relative performance of different parameterization methods by evaluating how model parameter uncertainties propagate through the stochastic model [10,11]. Some authors have discussed the importance of functional analysis for the uses of PTFs and have reported the importance of functionally evaluating PTFs for real field-scale applications [7,9,12,13].
One of the most popular applications of PTFs is using model simulations for environmental and crop processes. Studies have been conducted to evaluate the functionality of PTFs for environmental simulations, such as water flow simulations [13,14,15,16,17], evapotranspiration simulations, water balance studies [10,11,18,19], and irrigation water estimations [5,20], using different models. However, no such comprehensive studies have been conducted for process-based crop models [5].
Process-based crop models derive their output based on soil, climate, plant, and management parameters; therefore, the accuracy model output depends on the accuracy of the input parameters. At present, the acquisition of historic climate information for specific locations is easy due to the availability of ground-level observations and grid-based generated datasets. Similar performances of gridded climatic data and observed data in yield predictions using process-based crop models have been reported [21,22,23]. Therefore, ground-level observations may be successfully replaced with grid-based generated data. However, knowledge of soil parameters is crucial, as these parameters are largely responsible for the variability of process-based crop model outputs [24,25,26] and are often taken from nearby soil surveys or expert knowledge. However, due to its inherent variations, this information may need to be acquired at a finer scale for better estimations using process-based crop models [25].
Crop modeling aims to describe and understand the relevant plant processes and their interactions with abiotic and biotic factors [27]. Agricultural Production Systems sIMulator (APSIM) [28], Decision Support Systems for Agrotechnology Transfer (DSSAT) [29], and Simulateur mulTIdisciplinaire pour les Cultures Standard (STICS) [30] are some of the most widely used dynamic crop models to analyze agricultural processes under different input parameters, such as abiotic and biotic factors, crop parameters, and management conditions. However, APSIM has been applied widely because it has a user-friendly interface with advanced flexibility in establishing simple and complex simulations [31].
APSIM is an open-source (for non-commercial uses) crop modeling software that can be used to model the growth and yield of many crops [28,32]. This model is driven by daily climate data and can simulate the growth, development, and yield of crops and their interactions with soil and management conditions. Further, it includes modeling functions that allow one to simulate the soil water, nutrients, and much more. Plant models in APSIM simulate major physiological processes, such as phenology, water, and nutrient uptake, the development of organs and responses to abiotic stresses, etc. The soil module in APSIM simulates water movements, such as infiltration, capillary rise, evaporation, surface runoff, and drainage [28,32]. APSIM requires soil hydraulic inputs, such as the drained upper limit (DUL), drained lower limit, and crop lower limit (LL15 and crop LL), saturation (SAT), hydraulic conductivity (KS), and drainage coefficient (SWCON), for its simulations. It also uses some soil physicochemical parameters, such as soil bulk density (BD), organic carbon percentage (OC), pH, and the initial nutrient levels for simulations. Since the level of influence of these parameters may be different, it is vital to study the influence of these parameters on the model output using sensitivity analysis [33,34]. Sensitivity analyses are used to evaluate the responses of model outputs to changes in model input parameters and to quantify the relative importance of those responses [35,36]. Many studies have been conducted to evaluate the sensitivity of plant parameters and the management practices of different crops under different models [24,37,38,39,40,41,42,43]. Although the sensitivity assessment of soil is essential [25], limited attention has been paid to evaluating the sensitivity of soil parameters [42]. Hence, this study aims to evaluate the sensitivity of soil parameters to the growth and yield of rice using the APSIM–Oryza model in APSIM 7.10.
Gunarathna et al. [5] developed point PTFs to estimate the FC and PWP of Sri Lankan soils. Although they reported statistical validity, no comprehensive studies were conducted to test the functionality of these PTFs.
In this study, we aimed to assess the functionality of these PTFs for process-based crop models. We compared the rice growth and yield with different soil input conditions (such as measured soil data and PTF-estimated soil data) to assess the functionality of PTFs for APSIM–Oryza simulations under tropical Sri Lankan conditions. At first, we conducted a sensitivity analysis to confirm whether the parameters calculated by PTFs are sensitive in a simulation using APSIM–Oryza. Then, we confirmed the statistical agreement between the simulations using the measured soil data and calculated soil data with the PTFs. Moreover, we examined the required improvements of the PTFs to simulate the yield more accurately.

2. Materials and Methods

As a staple food, rice is the most commonly grown agricultural crop and accounts for 34% of the total cultivated area of Sri Lanka. Considering the social, cultural, and economic importance of rice, many studies (including crop modeling) have focused on the rice production systems of Sri Lanka. Hence, we selected rice to study the functionality of the developed PTFs in process-based crop models, particularly using APSIM version 7.10. Rice cultivation requires a vast amount of water. Consequently, rice is extremely vulnerable to drought damages due to the high spatial and temporal variation of the rainfall in Sri Lanka. In Sri Lanka, rice is cultivated in two major seasons (known as Maha and Yala) based on rainfall distribution and water availability [44]. In the Maha season (October to March), primarily rainfed rice cultivation is practiced in the lowlands with or without supplementary irrigation. In the Yala season, rainfed rice cultivation is limited due to low rainfall and irrigated rice cultivation is only practiced where irrigation is possible. Therefore, we ignored the Yala season (April to September) and rainfed rice cultivation in the Maha season (without supplementary irrigation) for this simulation study. The elements of this study are twofold [6]:
I. Sensitivity analysis is used to confirm the influential ability of FC and PWP for the outputs of the APSIM–Oryza model;
II. A study of the functionality of PTFs is developed to estimate FC and PWP for APSIM crop modeling.
For the sensitivity analysis, we selected three locations (1–3) representing the three climatic zones of Sri Lanka (a wet zone, with a mean annual rainfall of > 2500 mm, a dry zone with a mean annual rainfall of < 1750 mm, and an intermediate zone, with a mean annual rainfall of 1750–2500 mm), while 16 locations (including the above-mentioned three locations) were selected for the functional analysis (Figure 1). For both parts of the analysis, the required APSIM files were prepared as follows.

2.1. Meteorological Data

The AgMERRA climate forcing dataset provides daily, high-resolution meteorological data for the period of 1980–2010. We extracted the daily rainfall, maximum temperature, minimum temperature, and solar radiation of 16 selected locations from the AgMERRA global gridded climate dataset [45,46] using the NetCDF-Extractor V2.0 tool of AgriMetSoft (https://www.agrimetsoft.com).

2.2. Soil Data

The factsheets of the Sri Lankan, Canadian Soil Resource (SRICANSOL) project of the Soil science society of Sri Lanka contain a comprehensive dataset of Sri Lankan soils, including information on land use, as well as taxonomic, physical, and chemical properties [5,47,48,49]. This dataset includes information at the soil horizon level; therefore, the depth of the layers varies from location to location. We used the SRICANSOL database to prepare soil profiles of the above-mentioned 16 locations. We considered the depth of different soil horizons as the depth of the layer of the soil profiles. We collected information on the soil pH, organic carbon (OC), bulk density (BD), and moisture content of −33 kPa (DUL) and −1500 kPa (LL15) saturated hydraulic conductivity (KS) to develop soil profiles for the pertinent locations. The crop LL was considered to be similar to LL15. Other than the aforesaid measured data, air drying (considered as half of LL15) and saturation (considered as equal to the porosity) were also used as input parameters. Porosity (n) was calculated using the following equation, assuming 2.65 g/cm3 as the true density (TD). BD denotes the bulk density in Equation (1).
n = 1 B D T D
Gunarathna et al. [5] used point PTFs to estimate volumetric water contents using the multiple linear regression method. They developed PTFs to estimate volumetric water contents at −33 and −1500 kPa, representing the field capacity and permanent wilting point using various levels of input parameters. Considering their practical importance, the PTFs developed to estimate −33 and −1500 kPa, using the sand percentage of soil (Equation (2) and (3)), were used to estimate the DUL and LL15 for the input uncertainty assessment due to the parameter estimation using PTFs. We duplicated the above-mentioned 16 soil profiles and modified them by changing their respective PTF-estimated DUL (VWC33) and LL15 (VWC1500) values. Air-dried and crop-LL values were also modified accordingly.
VWC33 = 0.4357 − 0.035 × Sand%
VWC1500 = 0.3426 − 0.003 × Sand%

2.3. 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 [28,32]. To model crop growth, we used specific modules linked to the engine of APSIM, such as the specific crop model (APSIM–Oryza) and the soil water module (SOILWAT). The APSIM-crop model simulates crop-growth processes using the information on water from the SOILWAT module. APSIM uses the Oryza model [50,51] to simulate the daily rate of dry matter production and the phenological development of rice plants under different environments, including the ponding conditions [50].
The simulation of soil and crop variables started on 1 April of the year 2000 and continued until 31 December 2010, as this period represents dry, wet, and normal years. We considered the rice direct-seeding for this simulation, as it is the most common practice of field establishment for rice crops in Sri Lanka [52]. The cultivation decisions of Sri Lankan rainfed rice farmers are dependent on the onset of rainfall [52]. In simulations, we sowed rice when 20 mm of rainfall was accumulated over three days between November 15 to 30 or on December 1 in each year from 2000 to 2009. Considering the regular practices of Sri Lankan rice farmers, the planting density was set as 180 plants per square meter. We selected the Sri Lankan local cultivar BG357 for this study. We modified the Oryza model in APSIM 7.10 by adding the new cultivar BG357 and parameterized the cultivar using the values reported by Zubair et al. [53] and Nissanka et al. [54]. We set the maximum ponding as 100 mm during the period of 1 December until the end of February. Considering the fertilizer recommendations of the Department of Agriculture-Sri Lanka for rainfed rice, we set the fertilizer application to 30, 45, and 20 kg for the Urea nitrogen per hectare at 14, 28, and 42 days after the emergence, respectively. We assumed all the other management practices, including weed control, pest, and disease control, would follow the standard management practices recommended by the Department of Agriculture. We selected the total above-ground dry matter (WAGT) (kg/ha) and dry weight of rough rice (WRR) (kg/ha) as the output variables for this simulation study.

2.4. Sensitivity Analysis

Sensitivity analysis is (SA) a type of informative analysis used to determine which input parameter/s had the most influence on the model output. However, the purpose of the SA in this study was to confirm the influence of the parameters estimated by PTFs using the simulation of APSIM–Oryza. The SA could be local or global, according to the strategy used to explore the parameter space. Local methods could be relatively easily implemented and also require low computational costs compared to global methods. However, these methods have some shortcomings, including the inability to detect interactions between parameters and a heavy dependence on the base values of the input parameters [35,55]. A global sensitivity analysis overcomes the demerits of local methods, as it explores the entire multi-dimensional parameter space simultaneously. Hence, for a specific output variable, the influence of single parameters, as well as the interactions between parameters can be evaluated [35]. However, his process is computationally expensive, especially for process-based complex models [56].
An emulator is a simplified statistical approximation of a more complex model [57] used as a substitute for computationally expensive models. An emulator with sufficient accuracy can successfully replace an actual simulator to perform a sensitivity analysis [58]. In this emulator-based approach, a simulator that runs for a relatively smaller number of simulations could be used to build the emulator; then, the emulator could be used for the sensitivity analysis. Many studies used the GEM-SA software for sensitivity analyses of different process-based complex models [56,59,60,61,62]. Kennedy and O’Hagan [63]; Kennedy et al. [64]; and Kennedy and Petropoulos [59] reported a full description of the underlying mathematical procedures used in the GEM-SA software.
In this study, eight soil parameters (Table 1) were selected to study the sensitivity of the total above-ground dry matter (WAGT) and dry weight of rough rice (WRR) to rice crops under rainfed conditions. Initially, we generated 300 test points representing the entire range of respective inputs and related outputs (WAGT and WRR) of the APSIM simulator using the APSIM-R package [65] of the R software [66]. When generating 300 test points, the ranges were set by considering the typical ranges of Sri Lankan soils (Table 1). In this study, we used the Gaussian Process (GP) emulator in the the GEM-SA software [64] to develop 60 emulators (10 years × 2 outputs × 3 environments). Kennedy and O’Hagan [63]; Kennedy et al. [64]; and Kennedy and Petropoulos [59] offer a more detailed discussion of the analytical procedures of GEM-SA. In this study, we estimated the main effects and sum effects, assuming that all input uncertainties are unknown but uniform [62]. Further, we set the prior mean option as linear for each input [62]. We assessed the models using the leave-one-out cross-validation procedure embedded in GEM-SA [62]. With leave-one-out cross-validation, the emulator was used to estimate the emulator variance of each of the training runs in turn [59]. A slightly modified emulator was used to estimate each training point using all other training points. The GEM-SA calculates several diagnostics, including the sigma square value, cross-validation root mean squared error (CV-RMSE), and cross-validation root mean squared standardized (CV-RMSSE) using the results of the cross-validation [59]. These diagnostics could be used to assess the accuracy of the emulator approximations. The sigma squared value also provides a better understanding of the emulator accuracy. In this study, we evaluated the emulator accuracy using the above built-in diagnostic options of GEM-SA [62].
C V R M S E = i = 1 n ( y i y ^ ) 2 n
C V R M S S E = i = 1 n ( ( O 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 [59,62].
The variance-based sensitivity index (main effect, Si) was estimated by partitioning the output variance induced by the variations in all input parameters under the assumption that all input uncertainties are unknown but uniform. Si represents the expected reduction in the output variance if parameter xi is known [61,62], and the relative importance of each parameter in terms of its effect on output uncertainty can be ranked using the Si values of the selected parameters [67]. 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, and Var{E(f((Xx_i)} is the variance in the expected output f(X) if xi is known.

2.5. Functional Analysis

We simulated the total above-ground dry matter (WAGT) (kg/ha) and the dry weight of rough rice (WRR) (kg/ha) for the BG357 cultivar of rice for the 16 selected locations. At first, we used the soil profiles developed using the measured soil data (hereafter known as the “measurement-based simulation”). Then, we modified the soil profiles by changing the DUL and LL15 using the PTF-derived values. Further, the air-drying (half of LL15) was also modified accordingly. Secondly, we simulated the WAGT and the WRR of BG357, the cultivar of rice for all 16 locations using the modified soil profiles (hereafter known as ‘PTF-based simulation’). We compared the results of both measurement-based and PTF-based simulations to study the effects of input uncertainty on the model outputs.
Four selected outputs of the APSIM soil module (daily infiltration (infiltration), daily runoff (runoff), the soil water content of the topsoil layer (SWS(1)), and the potential evapotranspiration after modification for green cover and residue (EOS)) were also simulated for both soil profiles and compared.
A simple linear regression was used to compare the paired data points obtained for the measurement-based and PTF-based simulations. Further, we checked the results of the variances and then evaluated them using a t-test assuming unequal variances. The root mean square error (RMSE) (Equation (6)), a widely used measurement to study the goodness of fit of observations and estimations [68], was used for this study. Nash–Sutcliffe efficiency (NSE) (Equation (7)) [69], the efficiency calculated as one minus the sum of the absolute squared differences between the predictions and the observations normalized by the variance of the observations [70], was also used to evaluate the goodness of fit of the two estimations. The index of agreement (d) (Equation (8)), a standardized measure of the degree of the model prediction errors, which varies between 0 (no agreement) and 1 (perfect fit) [71], was also used to study the goodness of fit of the observations and estimations.
R M S E = i = 1 n ( S i O i ) 2 n
N S E = 1 i = 1 n ( O i S i ) 2 i = 1 n O i O ¯ ) 2
d = 1 i = 1 n ( S i O i ) 2 i = 1 n ( | S i O ¯ | + | O i O ¯ | ) 2
where Si is the simulated value, Oi is the observed value, and O ¯ is the mean of the observations.
The Lin’s concordance correlation coefficient (CCC) integrates precision using Pearson’s correlation coefficient—which represents the proportion of the total variance in the observed data that can be explained by the model, and accuracy with bias—which indicates how far the regression line deviates from the concordance line [72]. The CCC ranges from −1 to 1, with a perfect agreement at 1. According to the CCC, this agreement could be classified as almost perfect (>0.99), substantial (0.95–0.99), moderate (0.90–0.95), or poor (<0.90) [73]. This method can legitimately calculate accuracy with few observations for agreement using a continuous measure obtained by the two methods [74]. We calculated the CCC using the epiR package [74] in the R software [66].

3. Results and Discussion

3.1. Sensitivity Analysis

First, we evaluated the accuracy of the emulators. The sigma squared value, a measure of the nonlinearity of the emulator, records values near zero when an emulator model shows linearity. It shows higher values when the emulator models show moderate or high nonlinearity; however, no cut off values were reported to explain the linearity or nonlinearity. Petropoulos et al. [56] reported that their emulators showed linearity or moderate nonlinearity when they recorded sigma squared values between 0.13 and 1.6. Gunarathna et al. [62] reported sigma squared values between 0.10 to 1.43, reported as good to moderate linearity. In this study, we observed sigma squared values between 0.02 and 0.48 (Appendix Table A1 and Table A2); hence, we confirmed the linearity of the emulator models, though no clear cutoff values were determined to define the linearity. The CV-RMSE values reported far lower values compared to the standard deviation of the simulator results. Thus, we confirm that the results of the emulator are close to the simulator results (Appendix Table A1 and Table A2). CV-RMSSE shows values near one when the emulator results are close to the simulator results, while higher and lower values indicate over and underestimations, respectively. The CV-RMSSE values were between 0.92 and 1.26 (Appendix Table A1 and Table A2). All of these diagnostics reported comparatively better values than the previously reported values of the emulator accuracy assessments, which were reported to be satisfactory emulators [56,59,62,75].
Figure 2 shows a boxplot diagram for the sensitivity index of the eight parameters selected for the model outputs (total above-ground dry matter (WAGT) (kg/ha) and the dry weight of rough rice (WRR) (kg/ha)). The results revealed that DUL is the most influential parameter for the simulation of WAGT and WRR for rice crops, while LL15 showed the second-highest sensitivity to the selected outputs. SWCON and SAT also showed slight sensitivity to the selected outputs. All other selected parameters recorded low sensitivity indices; hence, BD, OC, KS, and pH were identified as less-sensitive soil parameters for WAGT and WRR.
Similarly, Varella et al. [25] showed a higher sensitivity of field capacity when simulating the growth and yield of wheat and sugar beet using the STICS model. Chen et al. [76] also showed a higher sensitivity of wheat yield to the available water capacity.
The results confirm that DUL and LL15 are the most critical soil parameters for APSIM–Oryza model simulations. Hence, an accurate estimation of DUL and LL15 is of prime importance in APSIM–Oryza simulations compared to other selected soil parameters. SWCON and SAT also showed little sensitivity to the output; hence, a reliable method is required to estimate SWCON with reasonable accuracy. Further, a reasonable estimation of the SAT may help increase the accuracy of the simulations. Other selected soil parameters, such as BD, OC, pH, and KS, showed the lowest sensitivity for the selected model outputs. Hence, a rough estimation of these parameters may not contribute significant changes to the outputs. However, the estimation of most of those parameters is not laborious or expensive compared to DUL and LL15.

3.2. Functional Analysis

Few studies have reported on the development of PTFs for tropical Sri Lankan soils, such as an evaluation of existing PTFs [77] and developing PTFs to estimate PWP and FC [5,78]. In this study, we checked only the functionality of PTFs developed to estimate FC and PWP [5] because no empirical equations have been developed to estimate other soil parameters. However, these values may not significantly affect the outputs because the sensitivity of these parameters remains low compared to DUL and LL15. We used the measured parameters as the other input variables (except crop LL and air dry), neglecting the uncertainty of those parameters.
Measurement-based simulations showed a strong relationship with PTF-based simulations for rice crop growth and yield using the APSIM–Oryza model (Figure 3). Further, the selected outputs of the APSIM soil module (daily infiltration (infiltration), daily runoff (runoff), soil water content of the topmost soil layer (SWS(1)), and potential evapotranspiration after modification for green cover and residue (EOS)) also showed good agreement between measurement-based and PTF-based simulations (Figure 3).
Measurement-based and PTF-based simulations showed 0.89 and 0.90 R2 values for WAGT and WRR, respectively (Table 2). This ensures no significant differences between PTF-based simulations and measurement-based simulations. A comparison between the outputs of the rice simulations (WAGT and WRR) under both simulation scenarios also confirmed the acceptability of PTF-based simulations, as no significant difference (p > 0.05) between both types of simulations under a t-test was observed (Table 2). The reported Nash–Sutcliffe efficiency was 0.89 and 0.90 for WAGT and WRR, respectively. Hence, we confirmed the similarity of PTF-based and measurement-based simulations. Wilmott’s index of agreement also confirmed the agreement of the results of PTF-based simulations with measurement-based simulations for both WAGT (d = 0.97) and WRR (d = 0.99). The RMSEs for WAGT and WRR were 811 and 445 kg/ha, respectively, which also confirmed the similarity of both simulation scenarios, which were assumed to be within the bounds (half of the standard deviation (SD) of the measurement-based simulations) (WAGT, SD = 2485 kg/ha; WRR, SD = 1424 kg/ha).
The measurement-based and PTF-based simulations showed over 0.90 each for the R2 values for all selected soil module related daily outputs (SWS(1), Infiltration, Runoff, and EOS) (Table 2). This ensures that PTF-based simulations can perform as well as measurement-based simulations. Further, the comparison of soil module related outputs between both simulation scenarios also confirmed the acceptability of PTF-based simulations, as no significant difference (p > 0.05) was found between both simulation scenarios for infiltration and runoff (Table 2). The Nash–Sutcliffe efficiency showed a better agreement between the two simulation scenarios, as all outputs related to the soil module recorded an NSE over 0.85. Wilmott’s index of the agreement also confirmed the agreement of the results of the PTF-based simulation with the measurement-based simulation for all soil module related outputs, reporting over 0.95 for d. The RMSE for SWS(1), Infiltration, Runoff, and EOS were 0.05 mm/mm, 0.70 mm/day, 0.45 mm/day, and 0.35 mm/day, respectively, which also confirmed the similarity of both simulation scenarios, which were assumed to be within the bounds.
This result proved the ability of PTF-based soil profiles to simulate WAGT, WRR, SWS(1), infiltration, runoff, and EOS reasonably equally to the measured soil profiles. As both parameters derived using PTFs are the soil parameters most sensitive to the growth and yield of rice, the functionality of PTFs are remarkable for APSIM crop modeling. Therefore, PTF-based estimations of DUL and LL15 may provide a solution for the data unavailability for crop and environmental simulation purposes. As measuring LL15 and DUL is laborious and expensive, the application of PTF to develop soil profiles may help cut down unnecessary expenses and delays.
There are no studies on functionality assessments for crop models. Moeys et al. [13] conducted a functional assessment to check the reproducibility of PTF-estimated soil parameters to predict water flow and solute transport. They reported the capability of PTF to be reasonably good, with reported CCC values of 0.54–0.72. In this study, we observed CCC values from 0.91 to 1.00.
We confirmed the statistical agreement between the APSIM simulation with measured soil parameters and APSIM simulation with PTF-estimated soil parameters. However, WAGT, WRR, and SWS(1) showed a slight underestimation of PTF-based simulations compared to the measurement-based simulations. Hence, it is important to study the causes of underestimation and the possibility to minimize these underestimations for precise modeling assessments.
Figure 4 shows that the developed PTFs slightly underestimated the LL15 and DUL of the selected Sri Lankan soils. The underestimation of DUL and LL15 in the PTF estimations may be attributed to an underestimation of the yields.
We assessed the accuracy of PTFs using CCC. The CCC values for LL15 and DUL were 0.81 and 0.79, respectively, which were classified as poor [73]. Figure 5 shows the average (2000–2010) maha season WAGT and WRR of the 16 selected locations. The estimated WAGT and WRR showed substantially higher CCC values (0.95 and 0.97 for WAGT and WRR, respectively, which is classified as substantial [73]), thereby confirming a relationship close to 1:1 between the measurement-based simulations and the PTF-based simulations. Nasta et al. [11] conducted a functional assessment of PTFs to predict the soil water budget. They reported the poor performances of the water budget assessment (NSE = 75%) compared to the PTF-estimated variables (NSE = 85%). Therefore, we can conclude that the functionality of PTF is remarkable for APSIM simulations. Nemes et al. estimated the irrigation water requirement [20] and water transport modeling [79] using measured and PTF-estimated soil hydraulic parameters and reported the similarity of both input conditions. Gunarathna et al. [5] also reported the similar irrigation water estimations using PTF-based and measurement-based soil hydraulic data estimations. However, no details on the accuracy assessment have been reported.
Figure 6 shows the error variations of the different locations of PTFs and PTF-based simulations. The blue and red parts of the graph represent the overestimation and underestimation of AWC (DUL-LL15), respectively. From these graphs, we can confirm that the underestimation of AWC may cause underestimation of the yield (WAGT and WRR), and the overestimation of AWC may cause overestimation of the yield. The underestimation of AWC’s influence can also be confirmed. Since the DUL is most sensitive in the APSIM simulation, the improvement of the PTF with DUL may be important for more accurate APSIM simulations.
In most locations, both LL15 and DUL showed underestimations (Figure 4 and Figure 6). However, in PTF-based APSIM simulations, WAGT and WRR showed both overestimations and underestimations but without any patterns. Hence, the error distribution is random, confirming that the effect of the regular underestimation (of hydraulic parameters) by the PTFs developed by Gunarathna et al. [5] did not appear in the output of the APSIM yield simulations. In this way, the functionality of PTFs on APSIM simulations is further confirmed.

4. Conclusions

Among the soil input parameters, the APSIM–Oryza model is highly sensitive to the drained upper limit and drained lower limit and is moderately sensitive to the drainage coefficient and saturation. Therefore, the proper estimation of those parameters is vital to the growth and yield simulations of rice using the APSIM–Oryza model.
The characterization of the sensitivity of the model input parameters to the model outputs has essential implications for model design and interpretation, as well as for future data collection efforts. Therefore, the results of these sensitivity analyses can be used as a guide for future crop modeling research work.
Despite the higher sensitivity levels of the drained upper and lower limits, the simulated PTF-based soil profiles were comparable to the measurement-based soil profiles. As the variability in the simulated yields using PTF-based soil profiles is within the expected limits, the PTF can be considered a better model; we conclude that PTF-based estimations of the drained upper limit and lower limit can be successfully used to simulate crop growth and yield under similar environments. Hence, the functionality of PTFs for APSIM simulations has been confirmed, and the laborious and expensive measurements of drained upper and lower limits can be successfully replaced by PTF estimations.

5. Limitations

Due to data limitations, we used gridded climatic data for this study without any proper validation before using them. However, this may not significantly affect our results, as we compared two simulations using the same climatic data. Although the use of the gridded data has been validated in other environments, it is not yet validated in Sri Lanka; therefore, we suggest conducting studies to check the applicability of gridded climatic data for crop modeling studies.

Author Contributions

Conceptualization, M.H.J.P.G. and K.S.; data curation, M.H.J.P.G. and M.K.N.K.; formal analysis, M.H.J.P.G. and M.K.N.K.; methodology, M.H.J.P.G., M.K.N.K., and M.R.; supervision, K.S.; writing—original draft, M.H.J.P.G. and M.K.N.K.; Writing—review and editing, M.H.J.P.G., M.K.N.K., K.S., and M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The accuracy of the emulators developed to estimate the total above-ground dry matter (kg/ha) of rice (WAGT) for different environments in Sri Lanka.
Table A1. The accuracy of the emulators developed to estimate the total above-ground dry matter (kg/ha) of rice (WAGT) for different environments in Sri Lanka.
YearDry Zone SoilIntermediate Zone SoilWet Zone Soil
Sigma sqCV RMSECV RMSSESigma sqCV RMSECV RMSSESigma sqCV RMSECV RMSSE
Year 10.13615.71.120.171272.31.040.171272.31.04
Year 20.03233.11.110.14927.11.060.14927.11.06
Year 30.10732.61.010.281736.41.040.281736.41.04
Year 40.03150.90.920.04371.21.110.04371.21.11
Year 50.05450.21.100.241594.61.060.241594.61.06
Year 60.13998.21.110.371742.51.040.371742.51.04
Year 70.05284.11.180.261346.71.060.261346.71.06
Year 80.03373.41.010.181056.51.090.181056.51.09
Year 90.06438.51.080.11919.41.030.11919.41.03
Year 100.06268.31.140.211358.20.960.211358.20.96
Table A2. The accuracy of the emulators developed to estimate the weight of rough rice (kg/ha) of rice (WRR) for different environments in Sri Lanka.
Table A2. The accuracy of the emulators developed to estimate the weight of rough rice (kg/ha) of rice (WRR) for different environments in Sri Lanka.
YearDry Zone SoilIntermediate Zone SoilWet Zone Soil
Sigma sqCV RMSECV RMSSESigma sqCV RMSECV RMSSESigma sqCV RMSECV RMSSE
Year 10.14398.21.160.21867.91.040.14737.81.05
Year 20.07182.51.260.19631.41.030.11546.21.00
Year 30.04309.21.010.38777.01.050.22845.60.98
Year 40.12108.30.990.11303.11.040.09197.11.13
Year 50.02195.81.090.39942.71.000.21720.81.05
Year 60.04383.51.040.44752.71.020.4210091.05
Year 70.08169.01.100.38970.31.010.20584.31.07
Year 80.09197.81.190.22553.81.040.25524.61.06
Year 90.12180.51.130.26815.31.060.11315.11.10
Year 100.05176.41.160.48778.10.950.20616.91.01

References

  1. Rustanto, A.; Booij, M.J.; Wösten, H.; Hoekstra, A.Y. Application and recalibration of soil water retention pedotransfer functions in a tropical upstream catchment: Case study in Bengawan Solo, Indonesia. J. Hydrol. Hydromech. 2017, 65, 307–320. [Google Scholar] [CrossRef] [Green Version]
  2. Bouma, J. Using Soil Survey Data for Quantitative Land Evaluation. Adv. Soil Sci. 1989, 9, 177–213. [Google Scholar] [CrossRef]
  3. Patil, N.G.; Singh, S.K. Pedotransfer Functions for Estimating Soil Hydraulic Properties: A Review. Pedosphere 2016, 26, 417–430. [Google Scholar] [CrossRef]
  4. Botula, Y.-D.; Van Ranst, E.; Cornelis, W.M. Pedotransfer functions to predict water retention for soils of the humid tropics: A Review. Rev. Bras. Cienc. Solo 2014, 38, 679–698. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Shein, E.V.; Arkhangel’skaya, T.A. Pedotransfer functions: State of the art, problems, and outlooks. Eurasian Soil Sci. 2006, 39, 1089–1099. [Google Scholar] [CrossRef]
  7. Wosten, J.H.M.; Bannink, M.H.; Gruijter, J.J.D.E.; Bouma, J. A procedure to identify different groups of hydraulic-conductivity and moisture-retention curves for soil horizons. J. Hydrol. 1986, 86, 133–145. [Google Scholar] [CrossRef]
  8. Wosten, J.H.M.; van Genuchten, M.T. Using Texture and Other Soil Properties to Predict the Unsaturated Soil Hydraulic Functions. Soil Sci. Soc. Am. J. 1988, 52, 1762–1770. [Google Scholar] [CrossRef]
  9. Vereecken, H.; Diels, J.; Orshoven, J.; van Feyen, J.; Bouma, J. Functional evaluation of pedotransfer functions for the estimation of soil hydraulic properties. Soil Sci. Soc. Am. J. 1992, 56, 1371–1378. [Google Scholar] [CrossRef]
  10. Chirico, G.B.; Medina, H.; Romano, N. Functional evaluation of PTF prediction uncertainty: An application at hillslope scale. Geoderma 2010, 155, 193–202. [Google Scholar] [CrossRef]
  11. Nasta, P.; Romano, N.; Chirico, G.B. Functional evaluation of a simplified scaling method for assessing the spatial variability of soil hydraulic properties at the hillslope scale. Hydrol. Sci. J. 2013, 58, 1059–1071. [Google Scholar] [CrossRef] [Green Version]
  12. Basile, A.; Bonfante, A.; Coppola, A.; Mascellis, R.; De Bolognesi, S.F.; Terribile, F.; Manna, P. How does PTF interpret soil heterogeneity? A stochastic approach applied to a case study on maize in Northern Italy. Water 2019, 11. [Google Scholar] [CrossRef] [Green Version]
  13. Moeys, J.; Larsbo, M.; Bergström, L.; Brown, C.D.; Coquet, Y.; Jarvis, N.J. Functional test of pedotransfer functions to predict water flow and solute transport with the dual-permeability model MACRO. Hydrol. Earth Syst. Sci. 2012, 16, 2069–2083. [Google Scholar] [CrossRef] [Green Version]
  14. Varvaris, I.; Pittaki-Chrysodonta, Z.; Moldrup, P.; de Jonge, L.W.; Iversen, B.V. Combining Visible−Near-Infrared and Pedotransfer Functions for Parameterization of Tile Drain Flow Simulations. Vadose Zone J. 2019, 18. [Google Scholar] [CrossRef] [Green Version]
  15. Schneider, S.; Jacques, D.; Mallants, D. Inverse modelling with a genetic algorithm to derive hydraulic properties of a multi-layered forest soil. Soil Res. 2013, 51, 372. [Google Scholar] [CrossRef]
  16. Soet, M.; Stricker, J.N.M. Functional behaviour of pedotransfer functions in soil water flow simulation. Hydrol. Process. 2003, 17, 1659–1670. [Google Scholar] [CrossRef]
  17. Sun, W.; Yao, X.; Cao, N.; Xu, Z.; Yu, J. Integration of soil hydraulic characteristics derived from pedotransfer functions into hydrological models: Evaluation of its effects on simulation uncertainty. Hydrol. Res. 2016, 47, 964–978. [Google Scholar] [CrossRef] [Green Version]
  18. Garrigues, S.; Boone, A.; Decharme, B.; Olioso, A.; Albergel, C.; Calvet, J.-C.; Moulin, S.; Buis, S.; Martin, E. Impacts of the Soil Water Transfer Parameterization on the Simulation of Evapotranspiration over a 14-Year Mediterranean Crop Succession. J. Hydrometeorol. 2017, 19, 3–25. [Google Scholar] [CrossRef] [Green Version]
  19. Baroni, G.; Facchi, A.; Gandolfi, C.; Ortuani, B.; Horeschi, D.; Van Dam, J.C. Uncertainty in the determination of soil hydraulic parameters and its influence on the performance of two hydrological models of different complexity. Hydrol. Earth Syst. Sci. 2010, 14, 251–270. [Google Scholar] [CrossRef] [Green Version]
  20. Nemes, A.; Czinege, E.; Farkas, C. Use of simulation modeling and pedotransfer functions to evaluate different irrigation scheduling scenarios in a heterogeneous field. In Proceedings of the 19th World Congress of Soil Science, Soil Solutions for a Changing World, Brisbane, Australia, 1–6 August 2010; pp. 186–189. [Google Scholar]
  21. Lashkari, A.; Salehnia, N.; Asadi, S.; Paymard, P.; Zare, H.; Bannayan, M. Evaluation of different gridded rainfall datasets for rainfed wheat yield prediction in an arid environment. Int. J. Biometeorol. 2018, 62, 1543–1556. [Google Scholar] [CrossRef]
  22. Shelia, V.; Hansen, J.; Sharda, V.; Porter, C.; Aggarwal, P.; Wilkerson, C.J.; Hoogenboom, G. A multi-scale and multi-model gridded framework for forecasting crop production, risk analysis, and climate change impact studies. Environ. Model. Softw. 2019, 115, 144–154. [Google Scholar] [CrossRef]
  23. Matthew, O.J.; Abiodun, B.J.; Salami, A.T. Modelling the impacts of climate variability on crop yields in Nigeria: Performance evaluation of RegCM3-GLAM system. Meteorol. Appl. 2015, 22, 198–212. [Google Scholar] [CrossRef]
  24. Liu, M.; He, B.; Lü, A.; Zhou, L.; Wu, J. Parameters sensitivity analysis for a~crop growth model applied to winter wheat in the Huanghuaihai Plain in China. Geosci. Model Dev. Discuss. 2014, 7, 3867–3888. [Google Scholar] [CrossRef]
  25. Varella, H.; Buis, S.; Launay, M.; Guerif, M. Global sensitivity analysis for choosing the main soil parameters of a crop model to be determined. Agric. Sci. 2012, 3, 949–961. [Google Scholar] [CrossRef] [Green Version]
  26. Folberth, C.; Skalský, R.; Moltchanova, E.; Balkovič, J.; Azevedo, L.B.; Obersteiner, M.; Van Der Velde, M. Uncertainty in soil data can outweigh climate impact signals in global crop yield simulations. Nat. Commun. 2016, 7, 1–13. [Google Scholar] [CrossRef] [Green Version]
  27. Waha, K.; Huth, N.; Carberry, P.; Wang, E. How model and input uncertainty impact maize yield simulations in West Africa. Environ. Res. Lett. 2015, 10, 024017. [Google Scholar] [CrossRef]
  28. 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] [Green Version]
  29. 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. Eur. J. Agron. 2003, 18, 235–265.
  30. Brisson, N.; Gary, C.; Justes, E.; Roche, R.; Mary, B.; Ripoche, D.; Zimmer, D.; Sierra, J.; Bertuzzi, P.; Burger, P.; et al. An overview of the crop model STICS. Eur. J. Agron. 2003, 18, 309–332. [Google Scholar] [CrossRef]
  31. Archontoulis, S.V.; Miguez, F.E.; Moore, K.J. Evaluating APSIM maize, soil water, soil nitrogen, manure, and soil temperature modules in the Midwestern United States. Agron. J. 2014, 106, 1025–1040. [Google Scholar] [CrossRef]
  32. 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]
  33. Lawless, C.; Semenov, M.A.; Jamieson, P.D. Quantifying the effect of uncertainty in soil moisture characteristics on plant growth using a crop simulation model. Field Crops Res. 2008, 106, 138–147. [Google Scholar] [CrossRef]
  34. Gijsman, A.; Jagtap, S.; Jones, J. Wading through a swamp of complete confusion: How to choose a method for estimating soil water retention parameters for crop models. Eur. J. Agron. 2002, 18, 77–106. [Google Scholar] [CrossRef]
  35. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley & Sons: West Sussex, UK, 2008; ISBN 9780470059975. [Google Scholar]
  36. 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]
  37. Confalonieri, R.; Bellocchi, G.; Tarantola, S.; Acutis, M.; Donatelli, M.; Genovese, G. Sensitivity analysis of the rice model WARM in Europe: Exploring the effects of different locations, climates and methods of analysis on model sensitivity to crop parameters. Environ. Model. Softw. 2010, 25, 479–488. [Google Scholar] [CrossRef]
  38. Lamsal, A.; Amgai, L.; Giri, A. Modeling the sensitivity of CERES-Rice model: An experience of Nepal. Agron. J. Nepal 2013, 3, 11–22. [Google Scholar] [CrossRef] [Green Version]
  39. Xing, H.M.; Xu, X.G.; Li, Z.H.; Chen, Y.J.; Feng, H.K.; Yang, G.J.; Chen, Z.X. Global sensitivity analysis of the AquaCrop model for winter wheat under different water treatments based on the extended Fourier amplitude sensitivity test. J. Integr. Agric. 2017, 16, 2444–2458. [Google Scholar] [CrossRef] [Green Version]
  40. Corbeels, M.; Chirat, G.; Messad, S.; Thierfelder, C. Performance and sensitivity of the DSSAT crop growth model in simulating maize yield under conservation agriculture. Eur. J. Agron. 2016, 76, 41–53. [Google Scholar] [CrossRef]
  41. Eweys, O.A.; Elwan, A.A.; Borham, T.I. Integrating WOFOST and Noah LSM for modeling maize production and soil moisture with sensitivity analysis, in the east of The Netherlands. Field Crops Res. 2017, 210, 147–161. [Google Scholar] [CrossRef]
  42. Tan, J.; Cui, Y.; Luo, Y. Global sensitivity analysis of outputs over rice-growth process in ORYZA model. Environ. Model. Softw. 2016, 83, 36–46. [Google Scholar] [CrossRef]
  43. Zhao, G.; Bryan, B.A.; Song, X. Sensitivity and uncertainty analysis of the APSIM-wheat model: Interactions between cultivar, environmental, and management parameters. Ecol. Model. 2014, 279, 1–11. [Google Scholar] [CrossRef]
  44. Gunarathna, M.H.J.P.; Kumari, M.K.N. Rainfall trends in Anuradhapura: Rainfall analysis for agricultural planning. Rajarata Univ. J. 2013, 1, 38–44. [Google Scholar]
  45. 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] [Green Version]
  46. Ruane, A.C.; Cecil, L.D.W.; Horton, R.M.; Gordón, R.; McCollum, R.; Brown, D.; Killough, B.; Goldberg, R.; Greeley, A.P.; Rosenzweig, C. Climate change impact uncertainties for maize in Panama: Farm information, climate projections, and yield sensitivities. Agric. For. Meteorol. 2013, 170, 132–145. [Google Scholar] [CrossRef]
  47. Dassanayake, A.R.; Somasiri, L.L.W.; Mapa, R.B. Benchmark Soils of the Intermediate Zone: Factsheets; Mapa, R.B., Dassanayake, A.R., Nayakekorale, H.B., Eds.; Soil Science Society of Sri Lanka: Kandy, Sri Lanka, 2005. [Google Scholar]
  48. Dassanayake, A.R.; Silva, G.G.R.; De 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: Kandy, Sri Lanka, 2010. [Google Scholar]
  49. Senarath, A.; Dassanayake, A.R.; Mapa, R.B. Bench Mark Soils of the Wet Zone: Factsheets; Mapa, R.B., Somasiri, S., Nagarajah, S.L., Eds.; Soil Science Society of Sri Lanka: Kandy, Sri Lanka, 1999. [Google Scholar]
  50. Gaydon, D.S.; Probert, M.E.; Buresh, R.J.; Meinke, H.; Suriadi, A.; Dobermann, A.; Bouman, B.; Timsina, J. Rice in cropping systems-Modelling transitions between flooded and non-flooded soil environments. Eur. J. Agron. 2012, 39, 9–24. [Google Scholar] [CrossRef]
  51. Bouman, B.A.M.; Van Laar, H.H. Description and evaluation of the rice growth model ORYZA2000 under nitrogen-limited conditions. Agric. Syst. 2006, 87, 249–273. [Google Scholar] [CrossRef]
  52. Amarasingha, R.P.R.K.; Suriyagoda, L.D.B.; Marambe, B.; Gaydon, D.S.; Galagedara, L.W.; Punyawardena, R.; Silva, G.L.L.P.; Nidumolu, U.; Howden, M. Simulation of crop and water productivity for rice (Oryza sativa L.) using APSIM under diverse agro-climatic conditions and water management techniques in Sri Lanka. Agric. Water Manag. 2015, 160, 132–143. [Google Scholar] [CrossRef]
  53. 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. In Handbook of Climate Change and Agroecosystems; Imperial College Press: London, UK, 2015; pp. 315–352. [Google Scholar]
  54. Nissanka, S.P.; Karunaratne, A.S.; Perera, R.; Weerakoon, W.M.W.; Thorburn, P.J.; Wallach, D. Calibration of the phenology sub-model of APSIM-Oryza: Going beyond goodness of fit. Environ. Model. Softw. 2015, 70, 128–137. [Google Scholar] [CrossRef]
  55. Wallach, D.; Makowski, D.; Jones, J.W.; Brun, F. Uncertainty and Sensitivity Analysis. In Working with Dynamic Crop Models; Elsevier B.V.: Amsterdam, The Netherlands, 2014; pp. 161–204. ISBN 9780123970084. [Google Scholar]
  56. 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. Model. 2009, 220, 2427–2440. [Google Scholar] [CrossRef]
  57. O’Hagan, A. Bayesian analysis of computer code outputs: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 1290–1300. [Google Scholar] [CrossRef]
  58. 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]
  59. Kennedy, M.C.; Petropoulos, G.P. GEM-SA: The Gaussian Emulation Machine for Sensitivity Analysis; Elsevier: Amsterdam, The Netherlands, 2017; ISBN 9780128030318. [Google Scholar]
  60. 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]
  61. 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]
  62. 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. [Google Scholar] [CrossRef] [Green Version]
  63. 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]
  64. 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]
  65. 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 February 2019).
  66. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 5 February 2019).
  67. Oakley, J.E.; O’Hagan, A. Probabilistic sensitivity analysis of complex models: A Bayesian approach. J. R. Stat. Soc. Ser. B Stat. Methodol. 2004, 66, 751–769. [Google Scholar] [CrossRef] [Green Version]
  68. Wallach, D.; Makowski, D.; Jones, J.W.; Brun, F. Model Evaluation. In Working with Dynamic Crop Models; Elsevier B.V.: Amsterdam, The Netherlands, 2014; pp. 345–406. ISBN 9780123970084. [Google Scholar]
  69. Krause, P.; Boyle, D.P.; Bäse, F. Advances in Geosciences Comparison of different efficiency criteria for hydrological model assessment. Adv. Geosci. 2005, 5, 89–97. [Google Scholar] [CrossRef] [Green Version]
  70. Nash, J.E.; Sutcliffe, J.V. River Flow Forecasting Through Conceptual Models Part I—A Discussion of Principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  71. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
  72. Ojeda, J.J.; Volenec, J.J.; Brouder, S.M.; Caviglia, O.P.; Agnusdei, M.G. Evaluation of Agricultural Production Systems Simulator as yield predictor of Panicum virgatum and Miscanthus x giganteus in several US environments. GCB Bioenergy 2017, 9, 796–816. [Google Scholar] [CrossRef]
  73. McBride, G. A Proposal for Strength-of-Agreement Criteria for Lin’s Concordance Correlation Coefficient; National Institute of Water & Atmospheric Research: Hamilton, New Zealand, 2005. [Google Scholar]
  74. Stevenson, M.; Nunes, T.; Heuer, C.; Marshall, J.; Sanchez, J.; Thornton, R.; Reiczigel, J.; Robison-Cox, J.; Sebastiani, P.; Solymos, P.; et al. epiR: Tools for the Analysis of Epidemiological Data; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  75. 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]
  76. Chen, C.; Wu, R.; Fletcher, A.; Lawes, R.; Oliver, Y. Sensitivity of simulated yield of dryland wheat to uncertainty in estimated plant-available water capacity. In Proceedings of the 22nd International Congress on Modelling and Simulation (MODSIM-2017), Hobart, Australia, 3–8 December 2017; pp. 244–249. [Google Scholar]
  77. Gunarathna, M.H.J.P.; Sakai, K. Evaluation of pedotransfer functions for estimating soil moisture constants: A study on soils in dry zone of tropical Sri Lanka. Int. J. Adv. Sci. Eng. Technol. 2018, 6, 15–19. [Google Scholar]
  78. Gunarathna, M.H.J.P.; Sakai, K.; Nakandakari, T.; Momii, K.; Kumari, M.K.N. Machine learning approaches to develop pedotransfer functions for tropical Sri Lankan soils. Water 2019, 11, 1940. [Google Scholar] [CrossRef] [Green Version]
  79. Nemes, A.; Schaap, M.G.; Wösten, J.H.M. Functional Evaluation of Pedotransfer Functions Derived from Different Scales of Data Collection. Soil Sci. Soc. Am. J. 2003, 67, 1093. [Google Scholar] [CrossRef]
Figure 1. A map of Sri Lanka showing the study sites and climatic zones defined based on the mean annual rainfall.
Figure 1. A map of Sri Lanka showing the study sites and climatic zones defined based on the mean annual rainfall.
Agronomy 10 00285 g001
Figure 2. Variation of the the sensitivity index (Si) of selected soil parameters on the simulated total above-ground dry matter (WAGT) and simulated dry-weight of rough rice (WRR) for rice in Sri Lanka.
Figure 2. Variation of the the sensitivity index (Si) of selected soil parameters on the simulated total above-ground dry matter (WAGT) and simulated dry-weight of rough rice (WRR) for rice in Sri Lanka.
Agronomy 10 00285 g002
Figure 3. Relationship between use of measurement-based inputs and pedotrandfer functions derived (PTF-based) inputs on simulations of rice, (a) Total above-ground dry matter (WAGT) (kg/ha) of rice; (b) Dry weight of rough rice (WRR) (kg/ha); (c) Soil water content of the topsoil layer (SWS(1)) (mm/mm); (d) Daily infiltration (infiltration) (mm); (e) Daily runoff (runoff) (mm); and f) Potential evapotranspiration after modifications of green cover and residue (EOS) (mm).
Figure 3. Relationship between use of measurement-based inputs and pedotrandfer functions derived (PTF-based) inputs on simulations of rice, (a) Total above-ground dry matter (WAGT) (kg/ha) of rice; (b) Dry weight of rough rice (WRR) (kg/ha); (c) Soil water content of the topsoil layer (SWS(1)) (mm/mm); (d) Daily infiltration (infiltration) (mm); (e) Daily runoff (runoff) (mm); and f) Potential evapotranspiration after modifications of green cover and residue (EOS) (mm).
Agronomy 10 00285 g003
Figure 4. Relationship between the measured and PTF-estimated soil hydraulic parameters (whole soil profile) of the selected soil profiles; (a) lower limit (LL15) and (b) drainage upper limit (DUL).
Figure 4. Relationship between the measured and PTF-estimated soil hydraulic parameters (whole soil profile) of the selected soil profiles; (a) lower limit (LL15) and (b) drainage upper limit (DUL).
Agronomy 10 00285 g004
Figure 5. Relationship between measurement-based and PTF-based simulations of average yields (2000–2010); (a) The total above-ground dry matter of rice (kg/ha); WAGT, (b) dry weight of rough rice (kg/ha); WRR.
Figure 5. Relationship between measurement-based and PTF-based simulations of average yields (2000–2010); (a) The total above-ground dry matter of rice (kg/ha); WAGT, (b) dry weight of rough rice (kg/ha); WRR.
Agronomy 10 00285 g005
Figure 6. The error variation of PTF-based simulations of average yield (2000–2010) in different locations; (a) total above-ground dry matter (kg/ha); WAGT and (b) weight of rough rice kg/ha; WRR; LL15: Lower limit (mm/mm), DUL: drainage upper limit (mm/mm); The size of the bubble represents the magnitude of the error of the estimated yield.
Figure 6. The error variation of PTF-based simulations of average yield (2000–2010) in different locations; (a) total above-ground dry matter (kg/ha); WAGT and (b) weight of rough rice kg/ha; WRR; LL15: Lower limit (mm/mm), DUL: drainage upper limit (mm/mm); The size of the bubble represents the magnitude of the error of the estimated yield.
Agronomy 10 00285 g006
Table 1. Selected soil parameters for the sensitivity analysis.
Table 1. Selected soil parameters for the sensitivity analysis.
ParameterDescriptionUnitRange
MinMax
BDBulk densityg/cm30.82.0
LL15Drained lower limit (equals the permanent wilting point)mm/mm0.050.30
DULDrained upper limit (equals the field capacity)mm/mm0.0 *1.0 *
SATSaturationmm/mm0.0 **1.0 **
KSSaturated hydraulic conductivitymm/day0.110,000
OCOrganic carbon%0.11.5
pH 4.09.0
SWCONDrainage coefficient 01
* Change according to the LL15 and SAT; ** Change according to the BD, LL15, and DUL.
Table 2. Performance evaluation of the simulation accuracies of the measurement-based and PTF-based simulations.
Table 2. Performance evaluation of the simulation accuracies of the measurement-based and PTF-based simulations.
CriteriaWAGTWRRSWS(1)InfiltrationRunoffEOS
R20.900.900.910.990.980.95
RMSE 811 kg/ha445 kg/ha0.05 mm/mm0.74 mm0.45 mm0.35 mm
NSE0.890.900.850.990.980.95
d0.970.970.951.001.000.99
CCC0.940.950.911.000.980.99
SD2485 kg/ha1424 kg/ha0.12 mm7.84 mm3.47 mm1.63 mm
P(T <=t)0.420.600.060.830.740.02
WAGT: Total above-ground dry matter of rice (kg/ha), WRR: Dry weight of rough rice of rice (kg/ha), SWS (1): Soil water content of the topsoil layer (mm/mm), infiltration: Daily infiltration (mm), runoff: Daily runoff (mm), EOS: Potential evapotranspiration after modification for green cover and residue (mm); R2: Coefficient of determination, RMSE: Root mean square error, NSE: Nash–Sutcliffe efficiency, d: Index of agreement, CCC: Lin’s concordance correlation coefficient, SD: Standard deviation, P: Probability of t-test.

Share and Cite

MDPI and ACS Style

Gunarathna, M.H.J.P.; Sakai, K.; Kumari, M.K.N.; Ranagalage, M. A Functional Analysis of Pedotransfer Functions Developed for Sri Lankan soils: Applicability for Process-Based Crop Models. Agronomy 2020, 10, 285. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10020285

AMA Style

Gunarathna MHJP, Sakai K, Kumari MKN, Ranagalage M. A Functional Analysis of Pedotransfer Functions Developed for Sri Lankan soils: Applicability for Process-Based Crop Models. Agronomy. 2020; 10(2):285. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10020285

Chicago/Turabian Style

Gunarathna, M. H. J. P., Kazuhito Sakai, M. K. N. Kumari, and Manjula Ranagalage. 2020. "A Functional Analysis of Pedotransfer Functions Developed for Sri Lankan soils: Applicability for Process-Based Crop Models" Agronomy 10, no. 2: 285. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10020285

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