Simulation of Crop Growth and Water-Saving Irrigation Scenarios for Lettuce : A Monsoon-Climate Case Study in Kampong Chhnang , Cambodia

Setting up water-saving irrigation strategies is a major challenge farmers face, in order to adapt to climate change and to improve water-use efficiency in crop productions. Currently, the production of vegetables, such as lettuce, poses a greater challenge in managing effective water irrigation, due to their sensitivity to water shortage. Crop growth models, such as AquaCrop, play an important role in exploring and providing effective irrigation strategies under various environmental conditions. The objectives of this study were (i) to parameterise the AquaCrop model for lettuce (Lactuca sativa var. crispa L.) using data from farmers’ fields in Cambodia, and (ii) to assess the impact of two distinct full and deficit irrigation scenarios in silico, using AquaCrop, under two contrasting soil types in the Cambodian climate. Field observations of biomass and canopy cover during the growing season of 2017 were used to adjust the crop growth parameters of the model. The results confirmed the ability of AquaCrop to correctly simulate lettuce growth. The irrigation scenario analysis suggested that deficit irrigation is a “silver bullet” water saving strategy that can save 20–60% of water compared to full irrigation scenarios in the conditions of this study.


Introduction
Humanity's environmental footprint is unsustainable within the Earth's limited natural resources and assimilative capacity [1].Climate change and growth in the global population are increasing pressure on these scarce environmental resources, notably water [2][3][4].Particularly, increasing relative evapotranspiration from flow regulation and irrigation over the past century raises the global human water consumption and footprint [5].Improving food production with less water and benchmarking efficiency of resource use is therefore a great challenge of our time, and urgently needed to ensure food security [1,6,7].
Cambodia is considered to be the country most vulnerable to climate change in Southeast Asia [8].In recent decades, extreme events, such as floods and droughts, have negatively affected the livelihoods of farmers, especially in terms of the loss of crop production [9].Cambodian farmers are generally conscious of these changes and challenges [9].Guidelines for agricultural adaptation to improve crop productivity and the sustainability of the farming system and to minimise vulnerability to climate change, are therefore crucial [8,10].Currently, the production of vegetables, like lettuce, poses more challenges in term of managing irrigation water efficiently, due to the crop's sensitivity to water shortage [11][12][13].Lettuce, the most widely consumed leaf vegetable, is also one of the most widely cultivated vegetables in the world [14].It is also an important to local vegetable production in Cambodia [15,16].Improving strategies for vegetable farming productivity, including lettuce, for Cambodian farmers, is being increasingly considered [17].
Many irrigation strategies have been investigated for improving irrigation water productivity (IWP) during recent decades, with IWP defined as the ratio of agricultural output to the amount of irrigation water use [18].Full irrigation via water application with the crop evapotranspiration requirements (ETc) method is an effective irrigation practice for crop production [19][20][21][22].In traditional irrigation scheduling, a technique to meet full irrigation, as well, the soil moisture in the root zone is allowed to fluctuate between an upper limit approximating "field capacity" and the lower limit of the readily accessible water (RAW), referred to as "the threshold", somewhat above where a crop begins to experience water stress [23,24].These methods have been applied to improve crop water productivity in various regions of the world, including Asian regions [25][26][27][28][29][30].Nevertheless, deficit irrigation, as an adaptation strategy for regions with limited water resources or prone to drought, has been proven to be worth considering [31,32].
Deficit irrigation is an irrigation practice whereby a crop is irrigated with an amount of water below the full requirement for optimal plant growth, thereby saving water and minimising the economic impact on the harvest [18,19].By limiting water applications to drought-sensitive growth stages such as, the vegetative stages and the late ripening period, the aims of this approach is to maximise water productivity and to stabilise, rather than maximise yields [33].Water deficit can be defined at five levels: severe deficit (with soil moisture (SM) less than 50% of field capacity (FC)), moderate deficit (SM < 50-60% of FC), mild deficit (SM < 60-70% of FC), no deficit or full irrigation (SM > 70% of FC), and overirrigation (application above water requirements) [34].Crops under deficit irrigation will experience some level of water stress, and often have lower yields than fully irrigated plants [35].Deficit irrigation can allow irrigation water savings of up to 20-40% at yield reductions below 10% [36], and has been widely investigated in dry regions [36].Deficit irrigation can be based on applying irrigation water under crop evapotranspiration.Patanè et al. [37] found that deficit irrigation at 50% of ETc for tomato plants resulted in no biomass (B) loss and high irrigation water-use efficiency.Experimental results obtained by Abd El-Wahed et al. [38] suggested that deficit irrigation at 85% of ETc is favourable to save 15% of water provided, with no reduction in the bean crop.The study results of Samperio et al. [39] offered deficit irrigation at 20% and 60% of ETc during stage II and postharvest, respectively, to "Angeleno" Japanese plum as a water-saving strategy, without negatively affecting crop yield.Results from Yang et al. [40] confirmed that the yield loss for cotton was less than 10% under deficit irrigation of 70% of ETc and 85% of ETc.Meanwhile, crop sensitivity to water deficit can be affected by many factors, including climatic conditions, crop species and cultivars, and agronomic management practices, amongst others [34].Payero et al. [41] suggested that deficit irrigation is not a good strategy for improving the crop water productivity of maize in a semi-arid climate.A study on deficit irrigation treatment on lettuce showed that water stress caused by deficit irrigation at 20% and 40% of ETc significantly reduced leaf number, leaf area index, and dry matter accumulation [42].Final fresh weight was reduced by 20% to 30% when compared with full irrigation.Kuslu et al. [43] concluded that for lettuce grown in semi-arid regions, full irrigation should be used under no water shortage, and deficit irrigation by 60% of ETc could be used for 40% water saving with a 35.8% yield loss where irrigation water supplies are limited.
Elaborating irrigation strategies merely on the basis of field research is difficult and time consuming [44].Crop models are effective decision-support tools to investigate irrigation scenarios and to develop improved irrigation strategies [7,45,46].They can provide a rapid and reasonable accurate prediction of the response of agriculture over a range of environmental conditions [47].The model AquaCrop, developed by the Food and Agricultural Organisation of the United Nations (FAO), is a water-driven crop model that simulates daily crop growth (e.g., canopy cover and biomass production) and final crop yield, with a balance between accuracy, simplicity, and robustness in incorporating various agronomy practices [48,49].It is considered as a valuable tool for improving irrigation water productivity in crop production planning [6,50].AquaCrop has been calibrated and parameterised to various crops under various environmental and irrigation conditions, including barley [51], soybean [52], sunflower [53], cotton [54,55], corn [56], sugar beet [57], wheat [58,59], potato [60,61], cabbage [62], and rice [63].However, this has not yet been done in the case of lettuce.Most of these studies proved that the model is capable of accurately simulating crop growth and yield.However, some case studies still report some flaws in simulation of crop evolution and yield, especially under severe deficit irrigation and heat stress conditions.Adeboye et al. [64] found that biomass of soybean simulated by AquaCrop was overestimated under deficit irrigation conditions.Zeleke et al. [65] found that AquaCrop simulated the canopy cover and biomass growth of canola well, but the model was less satisfactory under severe water stress conditions in a semi-arid region.Similarly, a reduction in model reliability in biomass and canopy cover prediction for maize under the severe stress conditions of deficit irrigation in a tropical environment was indicated in a study of Greaves et al. [66].AquaCrop performed well in biomass simulation of potato in the experiment under deficit irrigation at 120, 100, 80, and 60% of ETc [67].However, the potato yield simulation was overestimated due to the heat stress, with the authors suggesting the incorporation of a temperature stress coefficient into AquaCrop when a crop is affected by high temperatures.Further research is therefore required to improve the performance of AquaCrop.Furthermore, its performance simulating lettuce growth in Cambodian conditions has not yet been tested.The main objective of this study is to improve the water productivity of lettuce under limited irrigations in the Cambodian climate.More specific objectives are (i) to parameterise the crop model AquaCrop using data from farmer fields, since lettuce is not yet available in the AquaCrop catalogue; and (ii) to assess the impact of water-saving scenarios in full and deficit irrigation in silico using this calibrated model.

Experimental Sites
The field experiments were conducted with lettuce plants (Lactuca sativa var.crispa L.) which are widely used in the study area, during a period from August to September 2017 in two experimental sites located in the villages of Chea Rov (site S1) (104 • 38 54.442 E 12 • 9 15.482 N) and Ou Roung (site S2) (104 • 37 16.24 E 12 • 11 52.518N) in the province of Kampong Chhnang, Cambodia (Figure 1).The total land area of the plots was 400 m 2 .Lettuce seeds were sown in standard trays (with 123 holes).After 15 days, seedlings were transplanted into raised bed rows (0.30 m in height and with bed tops 0.50 cm wide) and covered with plastic mulch with a planting density of 12 plants m −2 .The compost was basally applied at the rate of 20 ton ha −1 before transplantation.
Irrigation was carried out using a drip system, with emitters of constructor maximum discharge of 3 L h −1 spaced 0.10 m apart.A plastic cover was used to protect the crops from heavy rainfall.Nevertheless, due to the intense rain which flowed between the crop rows, water ponding at 20 cm below the top bed row level was observed between the lettuce rows at both sites during almost the entire growing period.This ponding kept the soil wet during the growing period, and had to be factored into the calibration of the lettuce growing curve.At site S2, irrigation was not applied after a week after planting, due to the benefit of water ponding.At site S1, even though there was also water ponding in the field, the irrigation was applied every other day.The irrigation was determined by checking soil moisture (SM) using the feel and appearance method of Klocke et al. [68].The irrigation was done when the SM was depleted below field capacity in the root zone at 5 cm, as lettuce have a root depth between 5-10 cm.

Climate Data
Weather data for the experimental sites were collected from a local meteorological station (104 • 40 21.767 E; 12 • 10 45.965 N) (Figure 1).Daily maximum and minimum temperature, relative humidity, wind speed, rainfall, and solar radiation were recorded automatically at a five minute time step.The daily reference evapotranspiration (ETo) for the growing season, used as input data in AquaCrop, was calculated using the ETo calculator based on the FAO's Penman-Monteith method [69] (Figure 2).

Soil Data
The physical and chemical soil characteristics which were measured are listed in Tables 1 and 2. The soil texture was measured using the pipette method [70].The bulk density was measured by the core method [71].Field capacity, wilting point, and saturated hydraulic conductivity were derived from tension and soil moisture monitoring, using MPS2 and 10HS probes and using inverse modelling as presented in [72].

Crop Data
Canopy cover was measured at three-day intervals during the growing stage.Four pictures of 1 m 2 were taken randomly using a digital compact camera (Nikon Coolpix p600, Tokyo, Japan) at a fixed height of 1 m above ground level.The canopy cover was analysed using image processing with ImageJ ® software (https://imagej.nih.gov).Aboveground dry biomass was determined by harvesting 10 heads at the surface level of each site, oven-drying plant samples at 70 • C for 48 h, and weighing them [73].

AquaCrop Model
The AquaCrop model is a crop water-driven productivity model developed by the FAO in 2009.A detailed description is presented in [49].Water is the key limiting factor for crop production in this model [74].Inputs for the AquaCrop model consist of weather data, crop, and soil characteristics (soil profile and groundwater), and field management practice or irrigation management practices [49].
Canopy cover is a crucial feature of AquaCrop [49].Under unstressed condition, the exponential growth equation to simulate canopy development for the vegetative stage is where CC is the canopy cover at time t and is expressed as fraction of ground covered, CC o is initial canopy cover size (at t = 0) as a fraction (%), and CGC is the canopy growth coefficient in fraction per growing degree day (GDD), a constant for a crop under optimal conditions, but modulated by stresses.
In the condition of water stress, the CGC is multiplied by a water stress coefficient of expansive growth (Ks exp ) (Equation ( 2)).
CGC adj = Ks exp .CGC where Ks exp ranges from 1 to 0, canopy growth begins to slow down below the maximum rate when soil water depletion reaches the upper threshold, and stops completely when the depletion reaches the lower threshold.Crop transpiration is proportional to the canopy cover and given by Ks sto is the stress coefficient for stomatal closure.Kc Tr is the crop transpiration coefficient (determined by canopy cover and Kc Tr,x ), Kc Tr,x is the coefficient for maximum crop transpiration, and ETo is reference evapotranspiration (mm).
Biomass production is computed from crop transpiration and crop water productivity normalised for ETo and CO 2 (Equation ( 4)).The extreme effect of low temperature on crop phenology, biomass accumulation, and harvest index, is considered with adjustment factors [67,75].
where B is biomass, Tr is crop transpiration (mm day −1 ), ETo is reference evapotranspiration (mm day −1 ), and K s b is the stress coefficient for low-temperature effects on biomass production.f WP is the adjustment factor to account for differences, if any exist, in the chemical composition of the vegetative biomass and harvestable organs.WP * is normalised crop water productivity, defined as the ratio of biomass produced to water transpired, normalised for the evaporative demand and CO 2 concentration of the atmosphere.
The AquaCrop stress indicators include water storage (not enough water), waterlogging (too much water), air temperature (too high or too low), and soil salinity stress (too high).

Model Parameterisation
The process of parameterisation is illustrated in Figure 3.The vegetative stage of lettuce refers to the growing period of lettuce growth after germination until harvest.A growing period during the vegetative stage of 59 days after transplanting was simulated in this study.As lettuce is a crop which is not yet parameterised in AquaCrop, calibration of the model involved adjusting the model parameters to make them match the observed data [54,77].
The primary variables of lettuce growth, e.g., canopy cover and aboveground biomass were parameterised.For the calibration of the curves, the measured data in two experimental fields at the Chearov site (S1) (having sand soil) and Ourong site (S2) (with loam soil) were used, during the growing season in 2017.The AquaCrop model does not allow the use of observed data to build the canopy cover and biomass curves, but allows the data to be used to calibrate the canopy cover and biomass curves [78].
Canopy cover curves are a plot of the development of leaf expansion response to growing time per day, based on Equation (1).Biomass curves are a relationship plot of the growth of lettuce biomass response to growing time per day, based on Equation ( 4).The calibration of simulated canopy and biomass curves is based on one-at-a-time (OAT) methods (i.e., changing one parameter at a time while holding others constant) [79] and adjusting the parameters by trial and error, by comparing simulated and observed field data, and minimising the function of root mean square error.
We parameterised the canopy cover curve, which is important to the model for transpiration and evaporation [78].The main parameters of Equation ( 1), e.g., CCo and CGC for canopy cover curve determination, were adjusted to match the observed canopy cover data.In addition, adjusting the maximum canopy cover (CCx), time to reach maximum canopy cover, and time to recover, is crucial in order to obtain correct simulations of canopy cover growth.Subsequently, the focus was on adjusting the biomass curve of Equation ( 4).WP* and Kc Tr,x (coefficient for maximum crop transpiration) are the main parameters for regulating biomass curves in AquaCrop [74].As lettuce is a C3 crop type [80], the recommended values for WP* lie between 15 and 20 g m −2 .All calibrated crop parameters are shown in Table 3.The model performance for canopy cover and biomass simulation was evaluated using statistic indicators, including root mean square error (RMSE), Nash-Sutcliffe coefficient (N), and coefficient of determination (R 2 ), defined as below.
where O and S are the observed and simulated values at time i, respectively, and n is the total amount of the data.When N and R 2 are close to 1, it is considered to be satisfactory [81].RMSE should be close to 0.
AquaCrop requires the selection of inputs related to the irrigation method, such as sprinkler, drip, or surface.These methods determine the fraction of the soil surface made wet by irrigation [82] and the impact on irrigation efficiency [83].
Default AquaCrop settings for field management include mulching, and use an adjusted factor for the effect of mulches on soil evaporation.It varied between 0.5 for mulches derived from plant material, and 1.0 for plastic mulch [75].
The drip irrigation method with plastic mulch was applied as the input for field management in the model during the parameterisation, as this is the actual practice of the experiment in this study.
The soil water balance calculation, including soil moisture simulation in AquaCrop, is based on the storage capacity of the soil layers, described in Raes et al. [84], and previously in the BUDGET model [85].
During the experimental period, water ponding at 15 cm and 20 cm below the bed soil at site S1 and S2 respectively, which was observed during the experiment, was taken into account as a boundary condition during the parameterisation of the model.This water ponding resulted in wet soil during the growing period.The values of physical soil available data in the Section 2.2.2 were adopted to simulate soil moisture in this study.
It was noted that the plantation experiment was during the rainy season when irrigation was not needed.The crop parameters obtained after parameterisation are important for the investigation of the irrigation scenarios for water saving when irrigation is necessary, especially during the dry season.

Irrigation Scenarios
In the current study, AquaCrop was used to simulate the full and deficit irrigation scenarios described below (and in Table 4), in order to identify the optimal water use efficiency for lettuce.These irrigation scenarios applied irrigation scheduling based on soil moisture depletion [86] by applying readily available water depletion in the default option in AquaCrop.The time and irrigation dose were calculated with the criteria below: 1.
Soil water content depleted until a fixed lower threshold (RAW) and refill to field capacity (time criteria).
where ID is irrigation depth (mm), RAW = p TAW = p 1000(FC − PWP)Z r , p is soil water depletion threshold, set to 0.3 for lettuce recommended by [69], and Z r is root depth (m).TAW is the amount of water that a crop can extract from its root zone [88].FC is field capacity, that is, the amount of water well-drained soil should hold against gravitational forces (m 3 m −3 ) [88].PWP is permanent wilting point, referring to soil water content when a plant fails to recover its turgidity on watering (m 3 m −3 ) [88].RAW is readily available soil water, referring to the fraction of TAW that a crop can extract from the root zone without suffering water stress [88].AD is allowable depletion, defined as the percentage of RAW that can be depleted before irrigation water has to be applied.
Full irrigation scenarios with varied RAW thresholds were simulated by selecting allowable depletion levels at 0, 50, 80, 100% in AquaCrop, that avoid drought stress during the growth stage [41].The irrigation schedule is generated by selecting a so-called "time" and "depth" criterion, with "back to field capacity" and "allowable depletion", respectively.In other words, the different full irrigation scenarios result in decreasing irrigation frequency.
Deficit irrigation scenarios with varied RAW thresholds were similar to the full irrigation scenario criteria, but applied allowable depletion levels at 120, 130, 150, 180, and 200%.These levels result in drought stress during the growing stage, since soil moisture can decrease to a level below RAW before an irrigation event is triggered [41].The full irrigation scenario, based on a fixed irrigation frequency maintained the soil moisture in the root zone at field capacity on a daily basis, since the literature claims this is the optimal status to maximise lettuce yield [90].The irrigation schedule was generated with a fixed time interval (daily) (time criteria) and refill to field capacity (depth criteria).

Varied Field Capacity Threshold Irrigation Scenarios
Deficit irrigation scenarios with varied field capacity threshold reduce the irrigation dose below the dose at field capacity but keeping the same irrigation frequency, as in full irrigation scenario.Daily generated irrigation doses obtained in full irrigation scenario were reduced by 70, 60, 50, and 40%.
Irrigation water productivity (IWP) was used to evaluate the irrigation scenarios for efficient irrigation water use [31,91].IWP is the ratio between the yield and the irrigation water use [31].
where IWP is irrigation water productivity (kg m −3 ), Y is simulated yield (kg ha −1 ) and interest yield in this study is biomass, and I is irrigation water use (mm).The adjusted crop parameters obtained from the parameterisation process were used in the scenario simulation under the same weather conditions, using no soil surface cover in model field management, and no ground water at bottom soil profile boundary condition.

Plant Growth and Soil Moisture Status
Figure 6 shows both the lettuce growth measurement and simulation by AquaCrop.Biomass accumulated at a very low rate during the first two weeks of the growing season, and increased sharply in the final week.This trend accords with results obtained by Gallardo et al. [73].
The measured canopy cover and biomass yields were 34% and 0.11 ton ha −1 , respectively, at site S1 with sand soil, and 18.5% and 0.11 ton ha −1 , respectively, at site S2, which has loam soil.The measured results are comparable with Fazilah et al. [92], who found observed canopy cover of 33% and biomass yields of 0.22 ton ha −1 for lettuce under similar tropical conditions.Zhang et al. [93] found higher measured biomass for lettuce with a range of 0.33 to 0.63 ton ha −1 under lower temperatures of 20-25 • C. Thus, high day temperatures above 23 • C often limit lettuce production [94].Optimum growth for lettuce occurs between 15-20 • C [12].Unfavourable weather conditions, of high average temperature 33/25 • C (day/night) during the experiment, can be the reason of the low measured biomass yields for this study.

Model Parameterisation and Evaluation
The primary crop variables calibrated for daily lettuce growth were canopy cover and biomass, with the daily soil moisture simulated by AquaCrop, by adapting available physical soil data.
Table 5 presents the adjusted model parameters for canopy cover and biomass curve simulation of lettuce growth.The time to recovery of transplant, the time to reach the maximum canopy cover, the initial canopy cover (CCo), the maximum canopy cover growth coefficient (CCx), the coefficient for maximum crop transpiration (Kc Tr,x ), and the normalised biomass water productivity (WP*) were mainly calibrated.Note: C = conservative, NC = non-conservative.
WP* was adjusted at 16 gm −2 for both sites, within the recommended range.Kc Tr,x was adjusted at 0.65 and 0.5 for site S1 and S2, respectively.These adjusted Kc Tr,x are lower than crop coefficient for the mid-season (K cb,mid = 1) proposed by  The difference between the values proposed by FAO-56 and the adjusted Kc Tr,x values is due to the fact that the FAO crop coefficients were obtained for specific agroclimatic conditions, which are different from the conditions of this study [78].
In addition, Kc Tr,x is a major requisite for estimating crop transpiration and biomass.The low adjusted value of this parameter resulted in low simulated biomass yields to fit to measured values.
High temperature stress observed during the experiment could be the reason for the low observed lettuce biomass production [12].This observation leads to a recommendation for further development of a heat stress factor in relation to canopy cover and biomass simulations for lettuce.
The minimum root depth cannot be adjusted under 0.1 m, while the root development of lettuce was under this limit.Thus, root development in the model requires further modification [91].
The crop growth simulation of canopy cover and biomass fitted the observed data well (Figure 6).The statistical values for model evaluation in Table 6 were satisfactory, resulting in R 2 = 0.99, RMSE < 0.8%, N < 4.6 for canopy cover, and R 2 > 0.98, RMSE < 0.01 ton ha −1 , N < −0.07 for biomass.Thus, the model has ability to simulate well the growth of lettuce in both soil types at the two experimental sites.The measured and simulated soil moisture, at both soil depths of 5 and 15 cm in both sites, also matched well (Figure 7).The soil moisture simulation resulted in good accuracy with low RMSE of 0.18 and 0.14 m 3 m −3 at depths of 5 and 15 cm, respectively, at site S1, and 0.05 and 0.06 m 3 m −3 at depths of 5 and 15 cm, respectively, at site S2.

Irrigation and Soil Moisture Response
The cumulative irrigation in Figure 8, and the fluctuation of the soil moisture depletion in Figure 9, reflect the interaction between irrigation frequency and amount of water applied.In both varied RAW and field capacity threshold irrigation scenarios, the irrigation frequency decreased together with decreasing the amount of water applied per irrigation event.
In varied RAW threshold irrigation scenarios, the simulation of irrigation resulted in irrigation depths which ranged from 57 to 104 mm in site S1 (sand soil) and 46-82 mm in site S2 (loam soil) (Figure 8a,c).In varied field capacity threshold irrigation scenarios, irrigation depths ranged from 81-201 mm in site S1 and 83-209 mm in site S2 (Figure 8b,d).In varied RAW threshold irrigation scenarios, total simulated ETc ranged from 60 to 100 mm in site S1, and from 53 to 85 mm in site S2 (Figure 10a,c).The main reason for the higher ETc yield in site S1 is the higher adjusted transpiration characteristic of lettuce in sand soil as compared to loam soil.The simulated values of ETc fall within the range reported by Abdullah et al. [97] for lettuce, which varied from 43 mm to 285 mm in response to their different irrigation applications between 0 and 267 mm for open surface soil.

Crop Evapotranspiration and Biomass Growth Response
In varied field capacity threshold irrigation scenarios, simulated total crop evapotranspiration ranged from 77 to 205 mm in site S1, and from 83 to 211 mm in site S2 (Figure 10b,d).In both irrigation scenario classes, it was noted that while reducing irrigation events, crop evapotranspiration decreased simultaneously.
Figure 11 shows the response of biomass to the different irrigation scenarios.The varied RAW threshold irrigation scenarios (Figure 11a,c) resulted in biomass yield range from 0.88-1.77ton ha −1 at site S1, and 0.44-0.91ton ha −1 at site S2.By definition, biomass growth is closely related to crop evapotranspiration.Thus, the difference between biomass yields in the two experimental sites is due to the difference in the Kc Tr,x (coefficient for maximum crop transpiration) and CCx (maximum canopy cover) parameters between both sites.As expected, in varied RAW threshold irrigation scenarios, the simulations maintained biomass yield at 1.77 ton ha −1 at site S1 and 0.90 ton ha −1 at site S2 in the full irrigation scenarios with allowable depletion from 0-100% of RAW (e.g., S0RAW to S100RAW for site S1 and L0RAW to L100RAW for site S2), that is due to no-water stress condition.As the water stress started below the RAW line [41], with available depletion from 120-200% of RAW thresholds, the biomass yields decreased up to 50% in the S200RAW (200% of RAW threshold) scenario at site S1 and 52% in L200RAW scenario at site S2.
In varied field capacity threshold irrigation scenarios (Figure 11b,d), biomass yields ranged from 0.85 to 1.77 ton ha −1 at site S1, and 0.89 to 0.90 ton ha −1 at site S2.At site S1, reducing deficit irrigation at 50% of field capacity (S50FC scenario), the biomass yield started to decrease with 22% and deficit irrigation at 40% of field capacity (S40FC scenario), biomass yields decreased up to 51% compared to full irrigation scenario (S100FC).For site 2, deficit irrigation up to 40% of field capacity (L40FC) did not affect biomass yield.

Relationship between Water Productivity and Irrigation Scenarios
The responses of biomass yield and irrigation water productivity to irrigation depths in various scenarios are presented in Figure 12.Simulated water productivity of varied RAW threshold irrigation scenarios ranged from 1.5 to 2.1 kg m −3 for site S1 and 0.9 to 1.4 kg m −3 for site S2.In varied field capacity irrigation scenarios, simulated irrigation water productivity (IWP) ranged from 0.8 to 1.36 kg m −3 for site S1 and 0.43-1.08kg m −3 for site S2.The simulated irrigation water productivity results are comparable with other studies found in the literature.For instance, Gallardo et al. [98] found a measured IWP for lettuce dry matter of 1.86 kg m −3 .Figure 13 shows the relationship curves of biomass yield and irrigation water productivity response to irrigation scenarios.As expected, irrigation water productivity curve response to irrigation depths had parabolic relationships for both soil types in varied RAW threshold irrigation scenarios.Increasing water use efficiency can be enhanced by decreasing the irrigation to an optimum point.The optimum point, which resulted in 22% water saving for site S1, was found at the scenario with depletion of 150% of RAW (S150RAW), resulting in the irrigation water productivity = 2.07 kg m −3 , irrigation depth = 81 mm, and biomass yield = 1.68 ton ha −1 .For site S2, the optimum irrigation water productivity was at 130% of RAW scenario (L130RAW), resulting in irrigation water productivity = 1.42 kg m −3 , irrigation depth = 60 mm, and biomass yield = 0.85 ton ha −1 .In varied field capacity threshold irrigation scenarios, for site S1, the optimum irrigation water productivity with 39% water saving was found at deficit irrigation at 60% of field capacity (S60FC) with irrigation water productivity = 1.36 kg m −3 , irrigation depth = 130 mm, and biomass yield = 1.77 ton ha −1 .For site S2, the optimum water productivity resulted in 60% water saving, which was found at deficit irrigation at 40% of field capacity (L40FC scenario) irrigation water productivity = 1.08 kg m −3 , irrigation depth = 83 mm, and biomass yield = 0.89 ton ha −1 .
The varied RAW threshold irrigation scenarios resulted in higher simulated higher irrigation water productivity than the varied field capacity threshold scenarios in this study.Overall, deficit irrigation simulation scenarios in both irrigation scenario classes can provide a remarkable improvement in irrigation water productivity for water saving strategies.

Limitation
Crop models, like AquaCrop, are potentially valuable tools for answering questions primarily relating to research understanding, assessing crop management, and policy decision-making [49,99].However, it is essential to test the models in diverse field environments, such as those with varied temperatures, elevation transects, or amidst latitudinal variations [99].Particularly, AquaCrop has some limitations in terms of predicting crop yields only at the single growth cycle, single field scale, and only factoring in vertical water balance.The results of this study, obtained using climate data and field observation data relating to lettuce from a single growth cycle experiment at farm scale, allowed important information to be obtained in terms of calibrating lettuce crop parameters for sand and loam soil, and assessing limited water irrigation scenarios in the Cambodian context.However, it remains limited and the uncertainty on parameters has to be kept in mind.This study should be repeated in a contrasting range of diverse environments.Climate conditions and different cultural practices are the variables that differentiate the scenarios between different sites [99,100].It has been emphasised that uncertainty model simulation results are themselves uncertain, due to known inadequacies of the model (residual errors in measurement) and due to unknown inadequacies of the model (by inputting new cultivars or different types of management, the model may be wrong in unsuspected ways) [101].Despite such limitations, AquaCrop has already proven its usefulness in practical applications, and should still be tested widely in broader crop management applications, in diverse field environments [99,100].

Conclusions
An AquaCrop model was parameterised to simulate the canopy cover and aboveground biomass growth of lettuce under drip irrigation and plastic mulching for both sand and loam soil in the tropical monsoon climate of Cambodia.The model simulated canopy cover (RMSE < 0.8%) and aboveground biomass (RMSE < 0.01 ton ha −1 ) in a satisfactory way after adjusting several key parameters, as mentioned in Farahani et al. [54].
Additionally, the results suggested that the incorporation of a heat stress factor affecting canopy cover and biomass growth is necessary to meet the conditions encountered in a tropical climate context.
Shortage of water in Cambodian agriculture has increased due to climate change, and this is a significant challenge facing farmers in their crop production.In this study, the AquaCrop model has helped to develop the simulation process for limited irrigation management strategies to maximise irrigation water productivity.To test the impact of different irrigation scheduling and water saving strategies, two scenario classes were explored: (i) varied readily available water (RAW) threshold irrigation and (ii) varied field capacity threshold irrigation scenarios.The irrigation scenario analysis proposed optimal irrigation strategies for lettuce.
For varied RAW threshold irrigation scenarios, the analysis proposed optimal simulated irrigation water productivity at scenarios of 150% of RAW (irrigation water productivity = 2.1 kg m −3 ) for sand and 130% of RAW (irrigation water productivity = 1.4 kg m −3 ) for loam soil.This can save 22% of water, and resulted in a biomass yield reduction of 5 and 2%, respectively, for sand and loam soil.For varied field capacity threshold irrigation scenarios, the optimal deficit irrigation depth was found at 60% of field capacity (irrigation water productivity of 1.4 kg m −3 ) for sand soil, and at 40% of field capacity (irrigation water productivity of 1.0 kg m −3 ) for loam soil.It can save water up to 39% and 60%, for sand and loam soil, respectively, maintaining biomass yields compared to full irrigation.These results suggest that deficit irrigation is worth considering as a water saving strategy for lettuce in the monsoon climate of Cambodia.
Overall, AquaCrop is a valuable tool to predict lettuce growth and to investigate different scenarios for providing irrigation scheduling strategies for water saving in Cambodia.However, further research is necessary to standardise the model parameters for lettuce in various irrigation management, environmental, and climatic conditions.

Figure 3 .
Figure 3. Flow chart of parameterisation of AquaCrop in this study (adjusted from[76]).T is temperature, ETo is potential evapotranspiration, gs is stomatal conductance, WP is water productivity coefficient, Ks is stress coefficient, Es is soil temperature, Tr is crop transpiration.

Figure 4 .
Figure 4. Schematic representation of the crop response to varied RAW threshold irrigation scenarios simulated by AquaCrop (adjusted from[77]).RAW is readily available water content, TAW is total available water content, CC is the simulated canopy cover, CC o is initial canopy cover size, CGC is canopy growth coefficient in fraction per growing degree day (GDD), Ks sto is the water stress for stomatal closure, Kc Tr is the crop transpiration coefficient (determined by CC and Kc Tr,x at maximum canopy cover), ETo is the reference evapotranspiration, Ks b is the stress coefficient for low-temperature effects on biomass production, f WP is the adjustment factor to account for differences in chemical composition of the vegetative biomass and harvestable organs, WP* is the normalised water productivity.

Figure 5
Figure 5 illustrated concept of the varied field capacity threshold irrigation scenarios.

Figure 5 .
Figure 5. Schematic illustration of the soil water reservoir concepts of varied irrigation depth under field capacity irrigation scenarios (adjusted from [89]).FC is field capacity, full ID is full irrigation depth.

Figure 6 .
Figure 6.Observed (Obs) data and simulation (Sim) of lettuce growth of AquaCrop: (a) canopy cover at site S1 (CC-S) (sand soil) and site S2 (CC-L) (loam soil); (b) aboveground biomass at site S1 (B-S) and site S2 (B-S).The error bars were based on 10 biomass samples, except the last observed, which was based on 60 samples at harvest time.Sim CC-S is simulated canopy cover at site S1, Sim CC-L is simulated canopy cover at site S2, Obs CC-S is observed canopy cover at site S1, Obs CC-L is observed canopy cover at site S2, Sim B-S is simulated biomass at site S1, Sim B-L is simulated biomass at site S2, Obs B-S is observed biomass at site S1, Obs B-L is observed biomass at site S2.

Figure 7 .
Figure 7. Simulated soil moisture and observed soil moisture data measured at depths of 5 cm (H1) and 15 cm (H2) using soil moisture sensor 10HS and soil potential MPS-2: (a) soil moisture at site S1; (b) soil potential at site S1; (c) soil moisture at site S2; (d) soil potential at site S2.DAP is day after planting, Sim SM is simulated soil moisture, Obs SM is observed soil moisture, IRRI is irrigation, h is soil potential.

Figure 8 .
Figure 8. Irrigation accumulation response to different scenarios: (a) varied RAW threshold irrigation scenarios at site S1 (sand soil); (b) varied field capacity threshold irrigation scenarios at site S1; (c) varied RAW threshold irrigation scenarios at site S2 (loam soil); (d) varied field capacity threshold irrigation scenarios at site S2.RAW is readily available water content, S0RAW-S200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for sand soil.L0RAW-L200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for loam soil.S40FC-S100FC refers to deficit irrigation at 40-100% of field capacity for sand soil.L40FC-L100FC refers to deficit irrigation at 40-100% of field capacity for loam soil.

Figure 9 .
Figure 9. Daily soil moisture (VWC) response to different scenarios: (a) varied RAW threshold irrigation scenarios at site S1 (sand soil); (b) varied field capacity threshold irrigation scenarios at site S1; (c) varied RAW threshold irrigation scenarios at site S2 (loam soil); (d) varied field capacity threshold irrigation scenarios at site S2.RAW is readily available water content, S0RAW-S200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for sand soil.L0RAW-L200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for loam soil.S40FC-S100FC refers to deficit irrigation at 40-100% of field capacity for sand soil.L40FC-L100FC refers to deficit irrigation at 40-100% of field capacity for loam soil.

Figures 10 and 11
Figures 10 and 11 illustrate the cumulative crop evapotranspiration (ETc) and cumulative biomass of lettuce, respectively, under various irrigation scenarios simulated with AquaCrop calibrated for lettuce.In varied RAW threshold irrigation scenarios, total simulated ETc ranged from 60 to 100 mm in site S1, and from 53 to 85 mm in site S2 (Figure10a,c).The main reason for the higher ETc yield in site S1 is the higher adjusted transpiration characteristic of lettuce in sand soil as compared to loam soil.The simulated values of ETc fall within the range reported by Abdullah et al.[97] for lettuce, which varied from 43 mm to 285 mm in response to their different irrigation applications between 0 and 267 mm for open surface soil.In varied field capacity threshold irrigation scenarios, simulated total crop evapotranspiration ranged from 77 to 205 mm in site S1, and from 83 to 211 mm in site S2 (Figure10b,d).In both irrigation scenario classes, it was noted that while reducing irrigation events, crop evapotranspiration decreased simultaneously.Figure11shows the response of biomass to the different irrigation scenarios.The varied RAW threshold irrigation scenarios (Figure11a,c) resulted in biomass yield range from 0.88-1.77ton ha −1 at site S1, and 0.44-0.91ton ha −1 at site S2.By definition, biomass growth is closely related to crop evapotranspiration.Thus, the difference between biomass yields in the two experimental sites is due to the difference in the Kc Tr,x (coefficient for maximum crop transpiration) and CCx (maximum canopy cover) parameters between both sites.

Figure 10 .
Figure 10.Crop evapotranspiration accumulation responses to different scenarios: (a) varied RAW threshold irrigation scenarios at site S1 (sand soil); (b) varied field capacity threshold irrigation scenarios at site S1; (c) varied RAW threshold irrigation scenarios at site S2 (loam soil); (d) varied field capacity threshold irrigation scenarios at site S2.RAW is readily available water content, S0RAW-S200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for sand soil.L0RAW-L200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for loam soil.S40FC-S100FC refers to deficit irrigation at 40-100% of field capacity for sand soil.L40FC-L100FC refers to deficit irrigation at 40-100% of field capacity for loam soil.

Figure 11 .
Figure 11.Biomass accumulation responses to different scenarios: (a) varied RAW threshold irrigation scenarios at site S1 (sand soil); (b) varied field capacity threshold irrigation scenarios at site S1; (c) varied RAW threshold irrigation scenarios at site S2 (loam soil); (d) varied field capacity threshold irrigation scenarios at site S2.RAW is readily available water content, S0RAW-S200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for sand soil.L0RAW-L200RAW refers to irrigation scenarios with irrigation at 0-200% of RAW for loam soil.S40FC-S100FC refers to deficit irrigation at 40-100% of field capacity for sand soil.L40FC-L100FC refers to deficit irrigation at 40-100% of field capacity for loam soil.

Figure 12 .
Figure 12.Comparison of biomass and water productivity response (IWP) to different irrigation scenarios.RAW is readily available water content.S0RAW-S200RAW refers to irrigation at 0-200% of RAW threshold irrigation scenarios for sand soil.L0RAW-L200RAW refers to refers to irrigation at 0-200% of RAW threshold irrigation scenarios for loam soil.S40FC-S100FC refers to deficit irrigation at 40-100% of field capacity for sand soil.L40FC-L100FC refers to deficit irrigation at 40-100% of field capacity for loam soil.

Figure 13 .
Figure 13.Relationship between biomass and irrigation water productivity responses to different scenarios: (a) at site S1 (sand soil) and (b) at site S2 (loam soil).I is irrigation, B is biomass, IWP is irrigation water productivity, RAW-IS is varied readily available water content threshold irrigation scenarios, FC-IS is varied field capacity threshold irrigation scenarios.

Table 1 .
Measured physical soil characteristics.

Table 2 .
Measured chemical soil characteristics.
Note: EC is electrical conductivity; OM is organic matter content; N is total nitrogen; P is available phosphorous; K is exchangeable potassium; CEC is cation exchange capacity.

Table 3 .
Calibrated parameters of lettuce growth.
Figure4presents the calculation process of varied RAW threshold irrigation scenarios.

Table 6 .
Statistical evaluation of model simulation.