Effect of different soil moisture regimes on plant growth and water use efficiency of Sunflower: experimental study and modeling

Climate change and increasing demand in non-agricultural sectors profoundly affect the availability and quality of water resources for irrigated agriculture. The FAO AquaCrop simulation model provides a sound theoretical framework to investigate crop yield response to environmental stress. This model has successfully simulated crop growth and yield as influenced by varying soil moisture environments for crops. Integrating crop models that simulate the effects of water on crop yield with targeted experimentation can facilitate the development of irrigation strategies for high yield procurement and improving farm level water management and water use efficiency (WUE) under climatic condition of District Hyderabad, Sindh, Pakistan. This study was based on completely randomized block design with three treatments including T1 (30% soil moisture depletion), T2 (50% soil moisture depletion) and T3 (70% soil moisture depletion) with three replicates. In order to determine the crop water requirements under desired treatments, the gypsum blocks were used for computing the daily soil moisture depletion. The result shows that total volume of water applied to crop under T1, T2 and T3 was 9689, 5200 and 2045 m3 ha−1, respectively. As a result, the grain yield under T1, T2 and T3 was 13.2, 12.1 and 14.3 t ha−1, respectively. These results advocate that total yield of crop under T1 and T2 was less as compared to T3. The T3 gave higher yield and WUE compared than other treatments. On the other hand, results revealed that the simulated sunflower yields showed a good agreement with their measured under T3. The simulated grain yield was 15.5 t ha−1, while the measured yield varied from 12.1 to 14.3 t ha−1. This study suggested that WUE under T3 was more as compared to T1 and T2. The results showed that the T3 gave the highest crop yield in relation to WUE and optimize yield of sunflower crop under water scarcity. The Aquacrop model could very well predict crop yield and WUE at T3 under experiential region for sunflower production.


Background
Sunflower (Helianthus annuus L.) belongs to the family Asteraceae. The sunflower plant originated in the eastern North America. According to the literature, in the late 1800s the sunflower was introduced in the Russian Federation, where it become a food crop and Russian farmers made significant improvements in the way of the cultivated sunflower. Since 3000 B.C., a wide range of uses of sunflower have been reported throughout the world such as ornamental plant, medicinal, alimentary, feedstock, fodder, dyes for textile industry, body painting, decorations. So it is one of the 4th most important edible oil crops in the world that due to its high content of

Open Access
Bulletin of the National Research Centre *Correspondence: ashutus.iwm@sau.ac.bd 2 Chinese Academy of Agricultural Sciences, Beijing, China Full list of author information is available at the end of the article Page 2 of 8 Soothar et al. Bull Natl Res Cent (2021) 45:121 unsaturated fatty acids and low cholesterol level (Ismail and El-Nakhlawy 2018). In Pakistan, although it was introduced as an oilseed crop 40 years back, but its expansion in acreage and production is fluctuating due to various production and socio-economic constraints.
Research work on this crop has shown that there is a great potential of growing it under all the soil and climatic conditions in rain-fed as well as irrigated farming system in different agro-ecological zones. It has great potential to bridge the gaps between the production and consumption of edible oil. Population growth, land use change, climate change and increasing demand in non-agricultural sectors profoundly affect the availability and quality of water resources for irrigated agriculture. Amid increasing concerns that water scarcity and food security are among the main problems to be faced by many societies in the Twenty-first century, a global challenge for the agricultural sector is to produce more food with less water. Irrigation strategies focusing on increasing agricultural water efficiency, like as deficit irrigation practices couple with crop modeling to investigate multiple alternatives, have a pivotal role towards in the sustainable agricultural development by providing less than the exact crop water requirements, specifically during drought-tolerant at different crop growth stages, crop grain yields can be stabilized and maximum crop water productivity attained. Judicious planning is required as supplying crops with less than their water that, can significantly affect crop growth and development, inevitably affecting crop yield, if water stress occurs during the susceptible growth stage. In each case, various water amounts can be supplied. Owing to their cost and time effectiveness, crop-water simulation models are ideally suited for the evaluation of irrigation strategies during different stages.
According to FAO, AquaCrop simulation model provides a sound theoretical framework examine crop grain yield response to environmental stress. This model has successfully simulated crop growth and yield as influenced by varying soil moisture environments for crops, lie as sunflower. Geerts et al. (2009) suggested that the AquaCrop model maintains a good balance between robustness and accuracy, and a noteworthy feature of the model compared to other simulated models. It requires advanced skill for its calibration or operation and does not require a large numbers of inputs. The relatively small number of input data describes the soil-crop-atmosphere climatic in which the crop growths, most of which can be derived by simple ways. Heng et al. (2009) andHsiao et al.(2009) reported that the AquaCrop model simulated maize development, grain yield and water variables, such as the evapotranspiration and WUE, to name a few, reasonably well in cases of non-limiting conditions. Still, some studies reported that model performance declines in estimating some variables in severe water stress environments.
The availability of a model adapted to local conditions should have a strong impact on planned expansions and would aid in assessing competing management alternatives and possible constraints. Integrating crop models that simulate the effects of water on crop yield with targeted experimentation can facilitate the development of irrigation strategies for high yield procurement and improving farm level water management and water use efficiency (Soothar et al. 2021).

Methods
The field trial was carried out at the experimental station of Department Irrigation and Drainage, Faculty of Agricultural Engineering, Sindh Agriculture University Tandojam, Pakistan during Rabi season. The experimental site is located at latitude: 25° 25′ 28″ N; longitude: 68° 32′ 25″ E; altitude: 26 m. The design of the study was based on complete Randomized block design (CRBD) with including 3-Treatments, i.e., T1 = 30 percentage soil moisture depletion (SMD), T2 = 50 percentage soil moisture depletion and T3 = 70 percentage soil moisture depletion under basin irrigation method with three replications (R). Total nine plots with an average field size of 3 by 3 m were prepared. The experimental layout is shown in Fig. 1.

Soil properties
Before experiment, the soil physical properties such as dry bulk density, field capacity, soil porosity and infiltration rate were measured. Similarly, the collected soil samples were analyzed and basis soil proprieties are present in Table 1.

Crop variety and fertilizers application
Hysun 39 variety of Sunflower crop was selected in this research work. Hysun 39 sunflower variety is very much popular in local market. Fertilizers application was applied as the recommendations of On-Farm Water Management authorized by (MINFAL 2005).

Depth and frequency of irrigation
In order to determine the crop water requirements under different treatments (30, 50 and 70 percentage soil moisture depletion), the soil moisture meter with gypsum blocks was used to compute the daily soil moisture depletion level and irrigated plots at designed irrigation regimes.

Irrigation water measurements
In order to apply the required depth of water to experimental plot, a cut-throat flume was installed at the head of the ridge for the measurement of irrigation water to be applied. The following formula was used to determine the flow rate described by Skogerboe et al. (1997), and calculated discharge was verified through volumetric method at source of irrigation water. where f = subscript denoting free flow; u = subscript denoting upstream; Q f = free flow, discharge rate, L 3 /T; where S = subscript denoting submerged flow; d = subscript denoting downstream; Q s = submerged flow discharge rate, L 3 /T; n f = submerged flow exponent = 1.939 (for 8″ × 1.5′ size flume); C s = submerged flow coefficient = 1.606 (for 8″ × 1.5′ size flume); h d = downstream flow depth, L; n s = submerged flow exponent, dimensionless = 1.728 (for 8″ × 1.5′ size flume); S = submergence which is defined by the ratio of downstream head with upstream; head = h d /h u ; dimensionless.

Measurements of plant growth and yield
The experimental site was visited on daily basis and measurement of different crop growth parameters at maturity stage, like as plant height (cm), stem girth (cm), flower diameter, seed index and crop yield, selected and tagged. According to Abdou et al. (2011), the girth is the circumference of the plant around stem. Stem girth was measured with the help of vernier calliper at 10 cm height from soil surface. Measuring tape was used for measuring the height (cm) from the ground level to the tip of the disk.

Validation of Aquacrop model
Evaluation is an important step of model verification. It involves a comparison between independent field measurements (data) and output created by the model. Soil water content over the root depth during entire period of experiment was considered in this study for model evaluation. Different statistic indices including coefficient of determination of normalized root mean square error (RMSE) and agreement (D-index) were employed for comparison of simulated against observed data. The normalized RMSE expressed in percent, and it was calculated according procedure given by Loague and Green (1991). where P i and O i refer to simulated and observed values of the study variables, respectively, e.g., days from planting to anthesis, days from anthesis to physiological maturity. M is the mean of the observed variable. Normalized RMSE gives a measure (%) of the relative difference of simulated versus observed data. The simulation is considered excellent if a normalized RMSE is less than 10%, good if the normalized RMSE is greater than 10% and less than 20%, fair if normalized RMSE is greater than 20 and less than 30%, and poor if the normalized RMSE is greater than 30% described by Jamieson et al. (1991). The index of agreement (D-index) proposed by Willmott et al. (1985) was estimated through following formula. According to the d-statistic, the closer the index value is to one, the better the agreement between the two variables that are being compared and vice versa.
where n is the number of observations, P i the predicted observation, O i is a measured observation, P

Water use efficiency
The water use efficiency (WUE) for all treatments was calculated by using following formula (Soothar et al. 2019); where WUA = Water use efficiency (kg m −3 ), Y = Yield of crop (kg ha −1 ), WR = Total water consumed for crop production (m 3 ha −1 ).

Statistical analysis
The data collected were statistically analyzed using (analysis of variance) techniques following CRBD during field experiment. Tukey test was performed to see the significance of plant growth parameters and yield using the grand mean under different irrigation regimes. The software packages STATISTIX 8.1 was used for statistical analysis.

Irrigation water used
Total volume of water applied to crop under all treatments is shown in Fig. 2, and detailed date vise volumes of water applied are presented in Table 2. It is apparent from Table 2 that the total volume of water applied to crop under T1, T2 and T3 was 8.72, 4.68 and 1.84 m 3 , respectively. It was further calculated as 9689, 5200 and 2045 m 3 ha −1 . These results revealed that total volume of water used under T3 was less as compared to T1 and T2.

Plant growth parameters
The obtained data in Fig. 3 showed that the average plant height under T1, T2 and T3 was 134, 135 and 128 cm, respectively (Fig. 3a). However, the average stem girth plant was 2.52, 2.48 and 2.27 cm under T1, T2 and T3, respectively (Fig. 3b). Sunflower flower diameter was observed 15.3, 14.1 and 14.4 cm under T1, T2 and T3, respectively (Fig. 3c). These results imply that plant growth response was better and approximately same under all treatments (Fig. 3). The average value seed Index was non-significant difference among all treatments (Fig. 3d). The results of the experiment revealed that the crop irrigated by T1 under basin irrigation produced higher seed Index 7.2 g followed by T2 and T3 (6.93 and 6.92 g, respectively). These results implied that seed Index of dry grain was better T1 over T2 and T3. The results revealed that the grain yield of sunflower under T1, T2 and T3 was 13.2, 12.1 and 14.3 t ha −1 under different treatments, respectively. These

Aquacrop model validation for sunflower yield under different irrigation regimes
Validation is an important step in model verification and is done by comparing independent field measurements (data) with the outputs created by the model. Soil water content and sunflower crop yield were considered in this study for model validation. Table 3 showed that the measured and simulated results for the validated data sets for sunflower yield during 2016-2017. The simulated sunflower yields showed a good agreement with their measured under T3. The simulated crop yield was 15.6 t  Page 6 of 8 Soothar et al. Bull Natl Res Cent (2021) 45:121 ha −1 , while the measured yield varied from 12.1 to 14.3 t ha −1 under T1, T2 and T3 during experimental period. This could possibly be due to the fact that the senescence of the canopy accelerates under severe water stress, and the underground root system might be restricted and prevented from extracting more deeply stored soil water, thereby limiting its water uptake.

Water use efficiency
Results of water use efficiency of sunflower under all treatments are illustrated in Fig. 4. These result showed that water used efficiency in T1, T2 and T3 was 2.36, 3.60 and 10.24. Therefore, it is clear from the above results that T3 used less water as compare to T1 and T2 and gave higher yield and water used efficiency than that of T1 and T2. Similarly, water use efficiency simulated through Aquacrop model was found 4.2 kg m −3 of water evapo-transpired. However, it consistently tended to overestimate water use efficiency by T1 and T2, similarly, underestimate by T3.

Discussion
Water irrigation plays an important parameter for plant growth and performance of sunflower production. Therefore, optimum water irrigation scheduling became very necessary for high yield of sunflower (Abdou et al.

2011
). Lamm et al. (2013) found that the sunflower yield declines with deficit irrigation and yield increases with marginal irrigation system. (Ghani 2000) observed that 60% of the depletion of available water is optimum for better sunflower production. Bashir and Mohamed (2014);Buriro et al. (2015) reported that the maximum yield of sunflower grain was obtained with full irrigation.
Our study quantity of water applied to crop under T1, T2 and T3 was 8.72, 4.68 and 1.84 m 3 , respectively. Less amount of water was used in T3 treatment. Therefore, highest WUE was 10.24 obtained from T3 treatment and it was 65 to 77% higher than T1 and T2 treatment. Most recent study found that full irrigation produced the sunflower yield of 2049 kg ha −1 , where deficit irrigation produced 1710 kg ha −1 . However, deficit irrigation increased 22% more WUE than full irrigation (Eltarabily et al. (2020). The highest yield of 3119 kg ha −1 was obtained with six irrigations in contrast to 2200 kg ha −1 with two irrigations (Ghani 2000). Four times irrigations were an optimal irrigation system for attaining higher economical yields (Buriro et al. 2015). Similarly, approximately maximum yields were achieved despite water deficit only sufficient irrigation water provided during flowering stage. Applying deficit irrigation excluding flowering stage saves more irrigation water, which can be efficiently covered more land with the objective of stable yield (Ali et al. 2007). In our results revealed that the total grain yield of sunflower under T1, T2 and T3 was 13.2, 12.1 and 14.3 t ha −1 , respectively, under different treatments. These results indicated that the highest yield of crop produced by T3 in comparison with T1 and T2. However, our simulated Aquacrop model found 15.5 t ha −1 of sunflower yield. From above discussion found that T3 used less amount of water and gave higher yield and WUE than that of T1 and T2 treatment. Similarly, WUE simulated  Page 7 of 8 Soothar et al. Bull Natl Res Cent (2021) 45:121 through Aquacrop model was found 4.2 kg m −3 of water evapo-transpired. Several studies (Damdar et al. 2003;Prasad et al. 1999;Taha et al. 2001) reported that plant height was increased with increase in irrigation levels. These results are in line with Imam et al. (1999) reported that head diameter of sunflower was increased with increase in number of irrigations. Taha et al. (2001) reported that the seed Index was linearly related to the amount of irrigations. Heng et al. (2009) reported that the much greater deviations under severe water stress, as compared to wellwatered treatments for maize and canola crops simulated by AquaCrop. Table 3 shows the statistical assessment of the AquaCrop model for experimental time under all treatments. The results of NRMSE (9%), d-index (0.99) and percent deviation (8.8%) are comparable with those obtained by others authors. According to Paredes et al. (2014), they suggested that the AquaCrop model performance was obtained relative to biomass and yield predictions, with RMSE lower than 11 and 9% of the average observed biomass and yield, respectively. Iqbal et al. (2014) also suggested that the AquaCrop is a valid model and can be used with a reliable degree of accuracy for optimizing winter wheat grain yield production and water requirement on the NCP in terms of model validation.

Conclusion
These results revealed that total volume of water used under T3 was less as compared to T1 and T2. The simulated sunflower yields showed a good agreement with their measured under T3 and the results shown that the NRMSE is 9%, d (0.99) and percent deviation with 8.8%. The AquaCrop model successfully predicted cop yield of sunflower. The crop parameters are adjusted to simulate crop yield for sunflower under different irrigation regimes. These adjustments were made to obtain more stable and closer relationships between the simulated values and the measured values. Plant growth parameters were non-significantly affected. Seed Index under T1, T2 and T3 was 7.2, 6.93 and 6.92 g, respectively. Similarly, the yield of sunflower crop under T1, T2 and T3 was 13.2, 12.1 and 14.3 t ha −1 , respectively. The T3 gave higher yield and water used efficiency than other Treatments. Similarly, water use efficiency simulated through Aquacrop model was found 4.2 kg m −3 of water evapotranspired. However, it consistently tended to overestimate WUE under T1 and T2, similarly, underestimate by T3. The results of this study on T3 give more crop yield in relation to water use efficiency and optimize yield of sunflower crop under water scarcity. The AquaCrop model could very well predict crop yield and water use efficiency at 70% soil moisture depletion level under experiential region for sunflower production. From the comparison of results, it is concluded that the AquaCrop model is able to accurately simulate sunflower yield under climate conditions. Abbreviations SMD: Soil moisture depletion; CRBD: Completely randomized block design; RMSE: Root mean square error; WUE: Water use efficiency; ANOVA: Analysis of variation.