Calibration and Evaluation of the FAO AquaCrop Model for Canola (Brassica napus) under Varied Moistube Irrigation Regimes

The AquaCrop model was calibrated and validated for canola (Brassica napus) under Moistube irrigation (MTI) and various water regimes [(i) 100%, (ii) 75%, and (iii) 55% of crop water requirement (ETc)] over two seasons, 2019 and 2020. The normalised root mean square (nRMSE), Model Efficiency (EF), R2, and the Willmot’s index of agreement (d) statistics were used to evaluate the model’s efficiency in simulating biomass (B), canopy cover (CC), yield (Y), and harvest index (HI). The calibration results indicated the model simulated with accuracy the CC (under 100% ETc R2 = 0.99, EF = 0.92, nRMSE = 6.4%, d = 0.98) and 75% ETc (R2 = 0.99, EF = 0.92, nRMSE = 10.3%, d = 0.98). The model simulated CC well for validation for 100% ETc (R2 = 0.97, EF = 0.93, nRMSE = 22.5%, d = 0.98) and 75% ETc (R2 = 0.84, EF = 0.45, nRMSE = 59.2%, d = 0.86) irrigation regimes. Final biomass simulations were reasonably good under 100% ETc, 75% ETc, and 55% ETc irrigation regimes (R2 > 0.90, d > 0.65). The study showed the usefulness of AquaCrop for assessing yield response of canola to full and deficit irrigation scenarios under MTI.


Introduction
The global agricultural water consumption utilises 70% of the world's freshwater. Water is a finite resource, and climate variability and change have exacerbated the natural resource's depletion. Burgeoning populations have also increased per capita water use, thus compounding the global freshwater water scarcity situation [1]. Modern agricultural practices in sub-Saharan Africa (SSA) are a dualistic exercise that meets the poor's food security needs and is also a primary economic driver [2]. As such, using advanced irrigation techniques will maximise water productivity (WP), and subsequently, increase yields. Climate variability and change threaten food security, and industrial crops are not spared. An increase in global temperatures will lead to high carbon concentrations and warmer temperature; this consequently impacts cool climate C3 (canola, flax, wheat, and soybean) industrial crops [3]. Expanding irrigation land under current irrigation technologies and strategies can accelerate water scarcity, i.e., an increase in demand will lead to water scarcity under the present climate change variability scenarios [4,5]. Henceforth, adopting efficient irrigation techniques and strategies can ameliorate the accelerated demand on the finite water resource [6].
Canola is a C3 crop of economic importance. The crop produces oilseed that is processed into oil products for human consumption [7], and it is also used for forage production and phytoremediation [8]. Canola is considered a "healthy" trade oil. It contains no cholesterol, thus reducing the risk of cardiovascular diseases [9]. This has subsequently increased AquaCrop is a water-driven model that simulates yield, biomass production, and beneficial water use [28][29][30]. Water-driven models are an attractive option compared to their counterparts because of the ease of use. They facilitate easy normalisation of WP parameter under different climatic conditions (evaporative demand and atmospheric carbon dioxide). AquaCrop use has been applied to different crops such as cowpea [31], groundnuts [32], wheat [33][34][35][36][37], maize [38][39][40], and much recently on leafy vegetables [41]. Zeleke et al. [42] used AquaCrop to simulate canola yields under rainfed and irrigated conditions. The percentage relative difference (%D) between the observed and simulated yield was 2.2%. This signified AquaCrop's capability in simulating canola grain yield. This study investigated the capability of MTI for canola production under full and deficit irrigation scenarios. Identifying optimal DI strategies can potentially save water without imposing yield penalties on the canola grower. To extend the study's applicability beyond location-specific results, the experiment adopted AquaCrop modelling software [28]. The model has been used in numerous studies [31][32][33][34]37,38,40,41,[43][44][45][46][47] to assess yield response to water stress; however, there is a need to calibrate and test the AquaCrop model for industrial crops such as canola under MTI. As mentioned prior, there is a gap in how canola performs under MTI water stress conditions. The study was premised on the hypothesis that AquaCrop cannot effectively simulate canola crop performance under varying MTI water regimes. The specific objectives for this study were to (i) calibrate AquaCrop for canola under MTI water stress conditions and (ii) evaluate its ability to simulate CC, biomass, yield, and evapotranspiration (ET) under local South African conditions.

Model Description
AquaCrop is a water-driven modelling software used to simulate plant growth processes such as canopy cover (CC), biomass accumulation, and yield [28]. The model simulates yield response to water, i.e., water productivity [41]. Water productivity (WP*) is one of the crucial variables together with simulated transpiration (T r ), and reference evapotranspiration (ET o ) required to compute daily biomass (B) production (Equation (1)) [42]. For this study, canola biomass referred only to the above-ground component. AquaCrop's calculation scheme includes simulating the water stored in the root zone. The water stress coefficient is instrumental in determining the harvest index (HI). Once B is determined, the crop yield (Y) is then computed as per Equation (2): where WP* = water productivity (g·m −3 ), B = biomass (g·m −2 ), T r = transpiration (mm), ET o = reference evapotranspiration (mm), HI = harvest index, and Y = yield (kg·ha −1 ). AquaCrop simulates water use as a function of four stress factors, namely (1) canopy expansion, (2) stomatal closure, (3) early canopy senescence, and (4) aeration stress [41,46,48]. The model is relatively easy to use, as it requires few explicit parameters and largely intuitive input variables [29]. AquaCrop is underpinned by two sets of parameters: conservative parameters and non-conservative parameters. The former does not change with time management and are applicable on a large spatial variation scale, whereas the latter change with time, management, and location [41,49]. (September-November), in a tunnel, which was not temperature-controlled but was designed to exclude rainfall. The 2019 growing season was used to calibrate the AquaCrop model, and the 2020 season was used for model validation. The experiment was a split-plot design that consisted of three MTI regimes, namely 100%, 75% ET c , and 55% crop water requirement (ET c ), under tunnel conditions measuring 30 m by 10 m. The ET c was computed according to Equation (3):

Experimental Design
where ET c = crop water requirement (mm·day −1 ), K c = crop coefficient, and ET o = evapotranspiration (mm·day −1 ). Each MTI regime comprised of four experimental plots measuring 2 m by 1 m. The MTI AquaCrop deficit irrigation schedules followed the procedure by Geerts et al. [50]. The varied irrigation scheduling is summarised in Table 1. The experiment was done under a tunnel to facilitate better control of water fluxes and exclude of rainfall. A 1-m buffer hydrologically separated each experimental plot; a 250-micron plastic sheeting was vertically inserted to a depth of 1 m in each buffer space. PR2/6 profile probe access tubes were installed in each plot for soil water measurement at depths of 10, 20, 30, 40, 60, and 100 cm. Soil water content (SWC) measurements were done weekly using a PR2/6 profile probe connected to an HH2 handheld moisture meter (Delta-T, Cambridge, UK). Kanda et al. [31] performed weekly SWC measurements for cowpea production under MTI and showed that there was minimal temporal and spatial SWC variability. The canola was nursed at the University of KwaZulu-Natal-Pietermaritzburg (UKZN-PMB) (29 • 37 34.0 S 30 • 24 11.9 E) Controlled Environment Facilities (CEF) for two months before transplanting to the Ukulinga farm. Soil water measurements commenced two weeks before transplanting. Each plot accommodated 18 plants resulting in 9 plants·m −2 . Heng et al. [40] adopted plant densities of 6-8 plants·m −2 to prevent canola lodging.

Model Parameters and Input Data
The following data were collected during the July 2019-October 2019 and the October 2020-December 2020 growing season.

Weather Data
HOBO temperature and relative humidity (RH) sensors (Onset Computer Corporation, USA) were installed in the greenhouse for additional data collection (Figures 1 and 2 and Table 2). The ET o for the local conditions (within), the greenhouse were calculated using the evapotranspiration (ET o ) calculator [51]. Some variables required for calculating ET o were obtained from the automatic weather station (AWS) situated 100 m away from the greenhouse. The AWS uses the CS-500 Vaisala probe (Campbell Scientific, Logan, UT, USA) to measure temperature and relative humidity (converted into vapour pressure deficit), L1-200 pyranometer (Campbell Scientific, Logan, UT, USA) to measure solar radiation, and the Penman-Monteith equation to calculate reference evapotranspiration. The signal was transmitted wirelessly, and downloadable files are available from the South African Sugarcane Research Institute (SASRI) weather data portal. USA) to measure temperature and relative humidity (converted into vapour pressure deficit), L1-200 pyranometer (Campbell Scientific, Logan, UT, USA) to measure solar radiation, and the Penman-Monteith equation to calculate reference evapotranspiration. The signal was transmitted wirelessly, and downloadable files are available from the South African Sugarcane Research Institute (SASRI) weather data portal.  The weather data were used to create the climate file (.CLI) in AquaCrop consisting of ETo (.ETO) and daily minimum and maximum temperature (.TMP). Solar radiation data were input into the ETo calculator [51] for computing ETo. There was no daily rainfall file (.PLU). The leaf area index (LAI) was measured every two weeks using the LAI 2200 Canopy Analyser (Li-Cor, Lincoln, NE, USA). Since AquaCrop uses canopy cover (CC), Equation USA) to measure temperature and relative humidity (converted into vapour pressure deficit), L1-200 pyranometer (Campbell Scientific, Logan, UT, USA) to measure solar radiation, and the Penman-Monteith equation to calculate reference evapotranspiration. The signal was transmitted wirelessly, and downloadable files are available from the South African Sugarcane Research Institute (SASRI) weather data portal.  The weather data were used to create the climate file (.CLI) in AquaCrop consisting of ETo (.ETO) and daily minimum and maximum temperature (.TMP). Solar radiation data were input into the ETo calculator [51] for computing ETo. There was no daily rainfall file (.PLU). The leaf area index (LAI) was measured every two weeks using the LAI 2200 Canopy Analyser (Li-Cor, Lincoln, NE, USA). Since AquaCrop uses canopy cover (CC), Equation  The weather data were used to create the climate file (.CLI) in AquaCrop consisting of ET o (.ETO) and daily minimum and maximum temperature (.TMP). Solar radiation data were input into the ET o calculator [51] for computing ET o . There was no daily rainfall file (.PLU).

Canopy Cover (CC)
The leaf area index (LAI) was measured every two weeks using the LAI 2200 Canopy Analyser (Li-Cor, Lincoln, NE, USA). Since AquaCrop uses canopy cover (CC), Equation (4) was used to convert LAI to CC. Mabhaudhi et al. [46] argued that diffuse non-interceptance (DIFN) (Equation (5)), which is an output of the LAI 2200, can be used to compute CC. The DIFN utilises gap fractions to estimate the sections not "fully" obscured by the growing canopy [46,52]. The DIFN value ranges from 0 (no sky visible to the sensor) to 1 (no canopy obscuring the sun): where cg f = canopy gap fraction at zenith angle θ (averaged over azimuth angle and horizontal area) [52]. Since the seedlings were transplanted, the initial canopy cover (CC o ) was calculated by Equation (6). The computed CC o was 4.5%: Soil samples were subjected to soil textural analyses using the hydrometer method. The experiment sampled five depths for textural analysis, and the resultant textural data were fed into the SPAW model (Saxton and Willey, 2005) to determine the saturated hydraulic conductivity (K s ) and the bulk density (BD) ( Table 3). Other soil hydraulic parameters total porosity (θ s ) and residual soil water content (θ r ) were laboratory determined using the soil-water retention pressure method [53][54][55]. The soil data were used to create the soil file in AquaCrop (.SOL). Table 3. Soil textural and soil hydraulic parameters. Notes: θ r = residual soil water content (SWC), θ s = total porosity, k s = saturated hydraulic conductivity, and BD = Bulk density.

Field and Water Management Practices
The experiment was done under greenhouse conditions; hence, no rainfed systems were considered. The canola was subjected to optimal and deficit irrigation (DI) regimes. The optimal conditions consisted of irrigating at 100% of the crop water requirements (100% ET c ), whereas the DI irrigation regimes consisted of 75% ET c and 55% ET c . The irrigation intervals were used to create the irrigation file (.IRR). The SWC data were one of the parameters used to create the observation file (.OBS). The .OBS file was used for calibration and validation, respectively. The experiment assumed zero fertility stress. Other field management practices considered were (i) no weeds, (ii) no mulch, and (iii) zero runoff.

Biomass (B)
The above-ground biomass (AGB) (g·m −2 ) was harvested three times during each growing season. To avoid border effects, the samples were collected from the middle row. The freshly collected leaves and stems were weighed and then oven-dried at 85 • C for four days until there was consistent mass. One of each irrigation regime's plot was dedicated to destructive sampling. The harvest index (HI) was calculated as per Equation (7): where HI = Harvest index (no units), Y = yield (g·m −2 ), and B = above ground biomass (g·m −2 ). Other crop parameters recorded were transplanting date, amount of irrigation water, agronomic practices, time to flowering, time to yield formation, time to senescence, and harvesting dates.
The experiments were carried out in a greenhouse; hence, rainfall was zero. MTI is a subsurface irrigation method; therefore, surface runoff assumed a zero value. The impermeable layer at Ukulinga farm lies at a depth of 60 cm; thus, it prevented drainage and capillary rise [31]. ET a was converted from mm to m 3 ·ha −1 by multiplying Equation (8) by (10) [51].

Model Calibration
The calibration involved fine-tuning the non-conservative parameters for the canola crop. Table 4 presents summarised conservative and non-conservative values derived from the experiment. The parameters were adopted by Zeleke et al. [42] for calibrating and testing the FAO AquaCrop model for canola in Wagga Wagga, Australia. The study adopted the canola crop files calibrated by researchers from Lethbridge University Alberta, Canada [57]. The crop file was calibrated for warmer and drier climates in Swift Current, Saskatchewan, Canada.
The study destructively measured the seedling leaf area (4.50 cm 2 ) of the canola shoots at 90% emergence. Other input parameters were minimum rooting depth at 90% emergence (5 cm) and maximum rooting depth at harvesting. The average maximum rooting depths were 15.69 cm, 16.24 cm, and 20.41 cm for the 100% ET c , 75% ET c , and 55% ET c irrigation regimes, respectively. The reference harvest index (HI o ) was computed using Equation (2). Fine-tuning the HI o resulted in adopting a value of 25% for good simulations.
The calibration involved adjusting the non-conservative parameters HI o , initial canopy cover (CC o ), and canopy growth coefficient (CGC) until the simulated CC, B, and Y closely matched the observed data. The time to flowering was measured from the day of transplanting, and it was defined as the time when 50% of the plants had visible yellow flowering. Length of the flowering stage was the date after 50% flowering to the date when 50% of the plants had formed pods [31,58]. The maximum rooting depth was measured from the fully matured harvestable plants.

Model Validation and Evaluation Statistics
The model was validated using an independent dataset from the 2020 growing season. The data consisted of optimum irrigation (100% ET c ) and two DI regimes of 75% ET c and 55% ET c . It was validated, similar to the calibration, for SWC, CC, final B, Y, and WP ET .
Statistical analyses were employed to assess the model's ability to simulate canola crop growth and yield under MTI. The study applied the following criteria; normalised root mean square error (nRMSE), Wilmott's index of agreement (d), Model efficiency (EF), and the R 2 to assess the model's performance. The selected criteria are defined by Equations (10)- (13) [54,59]: where O i and P i = observed and predicted value(s), respectively, O i = mean observed data, and m = number of observations. The error index nRMSE showed the model's performance but did not clearly indicate the degree of over or under-estimation, hence using the EF statistical tool in the analysis. The EF statistic measured the residual variance vs. the measured data variance, and it ranges from −∞ to 1. EF values between 0.0 and 1.0 are considered acceptable (Table 5); however, Yang et al. [60] asserted that there exists a positive and scattered correlation between EF and d. Thus, when estimating soil water content, a satisfactory agreement can be considered when EF is greater than or equal to −1 and when d is greater than or equal to 0.60. R 2 represents the goodness of fit between the observed and simulated values [31]. For R 2 , a range of 0.5-1.0 represents good collinearity between observed and simulated values [61]. Above ground biomass, yield, and ET differences were computed as percentage relative differences (%D) using Equation (14). Relative differences of ±10% were considered accurate, whilst differences of ±20% were deemed acceptable [29,32,45]:

The Effects of Water Regimes on Growth, Yield, and Water Productivity of Canola
The leaf area index, represented by CC x , was significantly high for the 100% ET c treatment ( Figure 3). CC x under 100% ET c was reached after approximately nine weeks after transplanting. Pavlista et al. [17] reached CC x after 10 weeks of planting under optimal irrigation conditions. The 75% ET c treatment recorded a 91% CC x and 85.7% CC x for season 1 (S1) and season 2 (S2), respectively. The 55% ET c treatment recorded a low 74% CC x during S1 and CC x of 86% for S2. The 55% ET c S2 observation contradicted the norm since severe deficit irrigation is reported to yield a reduced canopy cover. The CC x was reached at week eight and week seven after transplanting under 75% ET c and 55% ET c , respectively. Deficit irrigation allows early crop maturity and small canopy cover as a form of drought avoidance mechanism. Small canopy development occurs to minimise water losses through transpiration [31,62].
Agriculture 2021, 11, x FOR PEER REVIEW 10 of 18 1 (S1) and season 2 (S2), respectively. The 55% ETc treatment recorded a low 74% CCx during S1 and CCx of 86% for S2. The 55% ETc S2 observation contradicted the norm since severe deficit irrigation is reported to yield a reduced canopy cover. The CCx was reached at week eight and week seven after transplanting under 75% ETc and 55% ETc, respectively. Deficit irrigation allows early crop maturity and small canopy cover as a form of drought avoidance mechanism. Small canopy development occurs to minimise water losses through transpiration [31,62]. Soil water content (SWC) varied across the irrigation regimes (Table 6). Soil water content between the 100% ETc and 75% ETc did not differ significantly (p > 0.05). There was no significant difference between the 75% ETc and the 55% ETc irrigation regime. Under the 100% ETc irrigation regimes, the recorded yields were 1.32 ton·ha −1 , whilst under 75% ETc and 55% ETc, the yield was 0.73 ton·ha −1 and 0.56 ton·ha −1 , respectively, during S1. The recorded yields during S2 for the 100% ETc, 75% ETc, and the 55% ETc irrigation regimes were 1.48 ton·ha −1 , 1.15 ton·ha −1 , and 0.75 ton·ha −1 , respectively (Table 7). Soil water content (SWC) varied across the irrigation regimes (Table 6). Soil water content between the 100% ET c and 75% ET c did not differ significantly (p > 0.05). There was no significant difference between the 75% ET c and the 55% ET c irrigation regime. Under the 100% ET c irrigation regimes, the recorded yields were 1.32 ton·ha −1 , whilst under 75% ET c and 55% ET c , the yield was 0.73 ton·ha −1 and 0.56 ton·ha −1 , respectively, during S1. The recorded yields during S2 for the 100% ET c , 75% ET c , and the 55% ET c irrigation regimes were 1.48 ton·ha −1 , 1.15 ton·ha −1 , and 0.75 ton·ha −1 , respectively ( Table 7).
The recorded yields under 100% ET c were approximately consistent with Zeleke et al. [42], who obtained canola yields of about 1.75 ton·ha −1 using the Bln3343-Co0401 cultivar. Additionally, Zeleke et al. [63]  because limited irrigation water supply inhibits canopy growth. Small canopy size results in low biomass, which consequently affects pod formation and grain yield. Biomass accumulation was also influenced by deficit irrigation. Extreme deficit irrigation strategies are not suitable for canola crop growth and yield development. It is worth mentioning that the yields attained were under tunnel conditions and the referenced literature performed the experiments under field conditions. Thus, this study reveals that there is no significant effect in growing canola under tunnel conditions compared to field conditions under full and optimal irrigation.

Soil Water Content
Since AquaCrop is a water-driven model, the model was firstly calibrated for soil water content (SWC). SWC simulations, if done accurately, will improve the accuracy of the simulated biomass and yield [31]. The model satisfactorily simulated the SWC under the 100% ET c irrigation regime (R 2 = 0.99, nRMSE = 16.3%, and d = 0.44); thus, the model was successfully calibrated for SWC (Figure 4d). The EF was significantly low, considering that the model successfully simulated CC under the 100% ET c regime. The low EF can potentially be attributed to inherent errors experienced during the calibration for CC. Under 100% ET c , it is evident that the model over-estimated the SWC. This could be potentially attributed to discrepancies in initiating drainage under the continuous irrigation regime. Furthermore, MTI is a slow-release irrigation technology hence the delay in wetting the soil to field capacity. Zeleke et al. [42] attributed the same phenomenon to the lag in AquaCrop to initiate drainage. The model simulated the SWC under the 55% ET c regime well (R 2 = 0.98, EF = 0.93, nRMSE = 4.5%, and d = 0.98) (Figure 4f). Despite having simulated the CC under the 75% ET c irrigation regimes well, the model yielded average simulation statistics for SWC (R 2 = 0.30, EF = 0, nRMSE = 15.1%, and d = 0.53). Inherent modelling errors in simulating CC can be attributed to the poor R 2 value under the 75% ET c irrigation regime (Figure 4e).

Canopy Cover
The model successfully simulated the canopy cover for the 100% ET c treatment (R 2 = 0.99, EF = 0.92, nRMSE = 6.4%, and d = 0.98) (Figure 4a). Under the 75% ET c deficit irrigation regime ( Figure 4b The R 2 and d-index were within the acceptable range; however, the EF was very low, and the nRMSE was significantly high. This resulted from the model capturing poor plant establishment and poor crop development after the transplanting exercise. Zeleke et al. [42] noted a poor CC whenever AquaCrop picked poor crop establishment and development. A careful calibration for C3 crops under extreme water deficit is required to produce a smooth and fitting CC curve. AquaCrop has simulation inaccuracies when predicting CC under water stress conditions [39,42]. 0.50). The R and d-index were within the acceptable range; however, the EF was very low, and the nRMSE was significantly high. This resulted from the model capturing poor plant establishment and poor crop development after the transplanting exercise. Zeleke et al. [42] noted a poor CC whenever AquaCrop picked poor crop establishment and development. A careful calibration for C3 crops under extreme water deficit is required to produce a smooth and fitting CC curve. AquaCrop has simulation inaccuracies when predicting CC under water stress conditions [39,42].

Biomass (B) and Yield (Y)
The model generally simulated the biomass well ( Table 7). The model satisfactorily simulated the biomass accumulation for the 100% ETc and 75% ETc irrigation regimes (R 2 > 0.90, EF > 0.50, and d > 0.89). The nRMSE under 100 ETc was 6.1%, and under 75% ETc was 37.3%, signifying a high residual variance in estimating the biomass. The model simulated the biomass under the 55% ETc irrigation regime well (R 2 = 0.90, EF = 0.30, nRMSE = 26.9%, and d = 0.75) ( Figure 5). The nRMSE = 26.9% was seemingly high; however, Ahmadi et al. [39] asserted that an nRMSE < 30% could be acceptable for crop simulation models. Thus, AquaCrop was successfully calibrated for biomass accumulation. It is worth noting that under the 55% ETc irrigation regime, the model under-estimated biomass by 25.50% (Figure 5c). This is a common phenomenon with AquaCrop under deficit irrigation scenarios [31]. On the contrary, Zeleke et al. [42] showed that AquaCrop overestimated the canola biomass because of heat stress. This study, however, was carried out during the winter (cool) season. The model simulated biomass with deviations of −27.48%, −2.30%, and 20.31%. The deviations fell within the acceptable ranges; thus, further asserting that the model was successfully calibrated for biomass under MTI.

Biomass (B) and Yield (Y)
The model generally simulated the biomass well ( Table 7). The model satisfactorily simulated the biomass accumulation for the 100% ET c and 75% ET c irrigation regimes (R 2 > 0.90, EF > 0.50, and d > 0.89). The nRMSE under 100 ET c was 6.1%, and under 75% ET c was 37.3%, signifying a high residual variance in estimating the biomass. The model simulated the biomass under the 55% ET c irrigation regime well (R 2 = 0.90, EF = 0.30, nRMSE = 26.9%, and d = 0.75) ( Figure 5). The nRMSE = 26.9% was seemingly high; however, Ahmadi et al. [39] asserted that an nRMSE < 30% could be acceptable for crop simulation models. Thus, AquaCrop was successfully calibrated for biomass accumulation. It is worth noting that under the 55% ET c irrigation regime, the model under-estimated biomass by 25.50% (Figure 5c). This is a common phenomenon with AquaCrop under deficit irrigation scenarios [31]. On the contrary, Zeleke et al. [42] showed that AquaCrop over-estimated the canola biomass because of heat stress. This study, however, was carried out during the winter (cool) season. The model simulated biomass with deviations of −27.48%, −2.30%, and 20.31%. The deviations fell within the acceptable ranges; thus, further asserting that the model was successfully calibrated for biomass under MTI.
The model over-estimated yield simulations despite having simulated CC well. Yield simulations were in the over-estimation range of ±34-97% and an under-estimation under the 100% ET c water regime, all of which are deemed unacceptable ( Table 8). The inability of AquaCrop to simulate yield can be attributed to the low heat units available during the winter season in which the experiment was run. Spring canola cultivar is a cool-season crop that requires a substantial amount of heat units for optimal growth [17]. Hergert et al. [12] also attributed low grain yield for canola to frost. The model over-estimated yield simulations despite having simulated CC well. Yield simulations were in the over-estimation range of ±34-97% and an under-estimation under the 100% ETc water regime, all of which are deemed unacceptable ( Table 8). The inability of AquaCrop to simulate yield can be attributed to the low heat units available during the winter season in which the experiment was run. Spring canola cultivar is a cool-season crop that requires a substantial amount of heat units for optimal growth [17]. Hergert et al. [12] also attributed low grain yield for canola to frost. The model successfully predicted the grain WPET under 100% ETc and 75% ETc irrigation regimes and, more interestingly, under 55% ETc, since the yield was low. Under the 100% ETc, the observed WPET was 0.42 kg·m −3 , whilst the simulated WPET was 0.36 kg·m −3 (D = 14.29%). Under the 75% ETc, the observed and simulated WPET were 0.48 and 0.49 kg·m −3, respectively, whilst under the 55% ETc, the observed and simulated WPET was 0.26 kg·m −3 . The WPET under the 55% ETc represented an optimal calibration scenario, whilst under the 100% ETc irrigation regime, the model under-estimated the WPET by 13% and over-estimated WPET by 2.1% under the 75% ETc irrigation regime. The observed WPET under the 100% ETc and 75% ETc irrigation regimes slightly matched those reported by Kumar et al. [36] for potatoes (C3 crop) grown in saline soils. The reported WPET were in the range of 0.63-0.98 kg·m −3 , although the model exhibited a low EF of 0.27.

Model Validation
Model validation was done after the calibration exercise. An independent dataset from the 2020 growing season was used to validate the model. Canola was transplanted on 27/10/2020 and harvested on 4 Janyary 2021.

Water Productivity (WP ET )
The model successfully predicted the grain WP ET under 100% ET c and 75% ET c irrigation regimes and, more interestingly, under 55% ET c , since the yield was low. Under the 100% ET c , the observed WP ET was 0.42 kg·m −3 , whilst the simulated WP ET was 0.36 kg·m −3 (D = 14.29%). Under the 75% ET c , the observed and simulated WP ET were 0.48 and 0.49 kg·m −3, respectively, whilst under the 55% ET c , the observed and simulated WP ET was 0.26 kg·m −3 . The WP ET under the 55% ET c represented an optimal calibration scenario, whilst under the 100% ET c irrigation regime, the model under-estimated the WP ET by 13% and over-estimated WP ET by 2.1% under the 75% ET c irrigation regime. The observed WP ET under the 100% ET c and 75% ET c irrigation regimes slightly matched those reported by Kumar et al. [36] for potatoes (C3 crop) grown in saline soils. The reported WP ET were in the range of 0.63-0.98 kg·m −3 , although the model exhibited a low EF of 0.27.

Model Validation
Model validation was done after the calibration exercise. An independent dataset from the 2020 growing season was used to validate the model. Canola was transplanted on 27 October 2020 and harvested on 4 January 2021.

Soil Water Content
The model successfully simulated the SWC under the 100% ET c irrigation regime (R 2 = 0.90, EF = 0.37, nRMSE = 8.7%, and d = 0.83) and 75% ET c irrigation regime (R 2 = 0.91, EF = 0.17, nRMSE = 4.1%, and d = 0.79) (Figure 6d,e). The EF was relatively low, but it signified a generally good model performance for the crop models [60]. The model reasonably simulated the SWC under the 55% ET c irrigation regime (R 2 = 0.55, EF = 0.05, nRMSE = 9.6%, and d = 0.63). The observation was attributed to the poor CC simulations' errors in which the model under-estimated the canopy growth (Figure 6c). The model successfully simulated SWC during the flowering and yield formation stages across all three irrigation regimes. The evidence revealed the capability of AquaCrop to simulate soil water content with reasonable accuracy for canola grown under MTI.
The model successfully simulated the SWC under the 100% ETc irrigation regime (R = 0.90, EF = 0.37, nRMSE = 8.7%, and d = 0.83) and 75% ETc irrigation regime (R 2 = 0.91, EF = 0.17, nRMSE = 4.1%, and d = 0.79) (Figures 6d,e). The EF was relatively low, but it signified a generally good model performance for the crop models [60]. The model reasonably simulated the SWC under the 55% ETc irrigation regime (R 2 = 0.55, EF = 0.05, nRMSE = 9.6%, and d = 0.63). The observation was attributed to the poor CC simulations' errors in which the model under-estimated the canopy growth (Figure 6c). The model successfully simulated SWC during the flowering and yield formation stages across all three irrigation regimes. The evidence revealed the capability of AquaCrop to simulate soil water content with reasonable accuracy for canola grown under MTI.

Biomass and Yield
The model simulated the yield with accuracy. The observed deviations (%D) were 7.43%, −25.22%, and 12.0% for the 100% ETc, 75% ETc, and 55% ETc irrigation regimes, respectively ( Table 9). The findings concur with Zeleke et al. [42], who found the D = −2.1-12% for spring canola cultivars grown in Wagga Wagga, Australia. The results obtained during validation were relatively accurate than those obtained during calibration.

Biomass and Yield
The model simulated the yield with accuracy. The observed deviations (%D) were 7.43%, −25.22%, and 12.0% for the 100% ET c , 75% ET c , and 55% ET c irrigation regimes, respectively ( Table 9). The findings concur with Zeleke et al. [42], who found the D = −2.1-12% for spring canola cultivars grown in Wagga Wagga, Australia. The results obtained during validation were relatively accurate than those obtained during calibration. The model over-estimated the final biomass under 100% ET c and the 75% ET c irrigation regimes giving D ≥ ±54.47% deviations. The model reasonably simulated the biomass data under the 55% ET c irrigation regime (D = 20.12%). The presented evidence demonstrates that AquaCrop can confidently simulate crop yields and biomass for canola under various MTI regimes with necessary adjustments.

Water Productivity (WP ET )
The simulated WP ET across the three irrigation regimes matched those obtained during calibration; thus, AquaCrop was successfully calibrated and validated for simulating WP ET . Under the 100% ET c , the observed WP ET was 0.42 kg·m −3 , whilst the simulated was 0.36 kg·m −3 (D = 14.29%). Under the 75% ET c , the observed and simulated WP ET were 0.48 and 0.49 kg·m −3 , respectively, whilst under the 55% ET c , the observed and simulated WP ET was 0.26 kg·m −3 . The D values were within the good range (D ≤ ±15%). WP ET was high under the 75% ET c , signifying that optimal yields and WP ET can be achieved with optimal deficit irrigation management practices.

Conclusions and Recommendations
This study sought to calibrate and validate the FAO AquaCrop model for canola grown under MTI and local South Africa conditions. The study was premised on the hypothesis that AquaCrop cannot effectively simulate the yield response of canola under varying irrigation regimes. The study, thus, failed to reject the null hypothesis for the 100% ET c and 75% ET c and rejected the hypothesis for the 55% ET c irrigation regime. The model was successfully calibrated and validated for soil water content, canopy cover, biomass accumulation, final biomass, yield, and water productivity under 100% ET c and 75% ET c irrigation regimes. AquaCrop poorly simulated the canopy cover and the SWC under the extreme deficit irrigation regime (55% ET c ). The poor simulation results can potentially be attributed to canola's sensitivity to extreme deficit irrigation scenarios. The study revealed that good deficit irrigation regimes could achieve optimum canola growth. The 75% ET c irrigation regime had an optimal grain yield and relatively high water productivity (WP ET ) compared to the 100% ET c irrigation regime. Thus, appropriate deficit management practices can produce high biomass and lower yield penalties. The study revealed the capability of the AquaCrop model to simulate canola response to various irrigation regimes. It is recommended that the study be done in open field conditions and assess the reliability of the reported results in this study. Additionally, the authors recommend the study be carried over several DI regimes and investigate the WP ET and yield correlation. In addition, the field experiments will generate an independent dataset that will be used to further test AquaCrop.