Next Article in Journal
The Effect of Cooking and Simulated Digestion on the Antioxidants and Minerals in Rice Grains: A Predictor for Nutritional Efficiency
Previous Article in Journal
Germination and Vigor of Maize Seeds: Pilot-Scale Comparison of Low-Oxygen and Traditional Storage Methods
Previous Article in Special Issue
An Overview of Polymeric Hydrogel Applications for Sustainable Agriculture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Comprehensive Deficit Irrigation Strategies for Drip-Irrigated Cotton Using AquaCrop

by
Yalong Du
,
Qiuping Fu
,
Pengrui Ai
,
Yingjie Ma
* and
Yang Pan
College of Hydraulic and Civil Engineering, Xinjiang Agricultural University, Urumqi 830052, China
*
Author to whom correspondence should be addressed.
Agriculture 2024, 14(8), 1269; https://doi.org/10.3390/agriculture14081269
Submission received: 4 July 2024 / Revised: 26 July 2024 / Accepted: 30 July 2024 / Published: 2 August 2024

Abstract

:
The development of a crop production strategy through the use of a crop model represents a crucial method for the assurance of a stable agricultural yield and the subsequent enhancement thereof. There are currently no studies evaluating the suitability of the AquaCrop model for the drip irrigation of Gossypium barbadense in Southern Xinjiang, which is the primary planting region for Gossypium barbadense in China. In order to investigate the performance of the AquaCrop model in simulating the growth of cotton under mulched drip irrigation, the model was locally calibrated and validated according to different irrigation thresholds during a key growth period of two years. The results of the simulation for total soil water (TSW), crop evapotranspiration (ETc), canopy coverage (CC), aboveground biomass (Bio), and seed cotton yield demonstrated a high degree of correlation with the observed data, with a root mean square error (RMSE) of <11.58%. The Bio and yield simulations demonstrated a high degree of concordance with the corresponding measured values, with root mean square error (RMSE) values of 1.23 t ha−1 and 0.15 t ha−1, respectively. However, the predicted yield declined in the verification year, though the prediction error remained below 15%. Furthermore, the estimated evapotranspiration (ETc) value demonstrated a slight degree of overestimation. Generally, the middle and late stages of cotton growth led to an overestimation of the TSW content. However, the prediction error was less than 13.99%. Through the calculation of each performance index of the AquaCrop model, it is found that they are in the acceptable range. In conclusion, the AquaCrop model can be employed as a viable tool for predicting the water response of cotton to drip irrigation under mulched film in Southern Xinjiang. Based on 64 years of historical meteorological data, three years were selected as scenarios for simulation. Principal component analysis (PCA) showed that, in a local wet year in Southern Xinjiang, the irrigation quota was 520 mm, and the irrigation cycle was 6 days/time. In normal years, the irrigation quota was 520 mm, with an irrigation cycle of 6 days/time. In dry years, the irrigation quota was 595 mm, with an irrigation cycle of 10 days/time. This allowed for higher seed cotton yields and irrigation water productivity, as well as the maximization of cotton yields and net revenue in the arid oasis area of Southern Xinjiang.

1. Introduction

Global warming has led to a rise in the frequency of meteorological disasters, such as abnormally high temperatures, droughts, and heavy rainfall. These events create instability in agricultural production, which, in turn, affects crop yields [1,2]. Among the various abiotic stresses, drought is one of the most significant threats to agricultural production. Anomalous temperature increases and precipitation uncertainty directly affect soil moisture evaporation and plant transpiration in agricultural fields, which, in turn, affect crop water consumption. Furthermore, fluctuations in regional precipitation have implications for both local irrigated agriculture and water management [3,4]. Agriculture is a significant contributor to the economy of Xinjiang, and ensuring the continued production of food is crucial for maintaining social stability in the region. As a typical irrigated agricultural area in the arid oasis zone of inland Northwest China, Xinjiang tends to apply excessive amounts of water to ensure a suitable soil microenvironment for crop growth in arid or semi-arid regions due to high temperatures, and the high evaporative demand during the crop-growing season [5] further exacerbates water shortages in Xinjiang. Consequently, the efficient utilization of water resources in Xinjiang is crucial for the sustainable development of agriculture.
Gossypium barbadense, which is also known as long-staple cotton, is a high-end raw material for textiles. Due to its susceptibility to light and heat conditions, only the United States, Egypt, and a few other countries are able to cultivate it. With the advent of molecular genetics and biotechnology, new varieties of Gossypium barbadense with enhanced tolerance to the environment and superior fiber length have been successfully developed [6,7], paving the way for cultivation in arid regions. The demand for Gossypium barbadense has been on the rise in recent years due to its exceptional fiber quality [8]. It is easy to grow Gossypium barbadense in Southern Xinjiang because the weather is right for it. This plant has very good fibers and can be spun easily, so Southern Xinjiang has become China’s main source of high-quality Gossypium barbadense (cotton). The cultivation of cotton is susceptible to environmental fluctuations, with climatic conditions and agronomic practices exerting a significant influence on the crop’s yield and quality [9]. However, in the context of climate change and water scarcity [10], there is a dearth of research on irrigation strategies for cotton with low water consumption and high water use efficiency. Therefore, the effective reduction of irrigation water consumption for cotton represents a crucial approach to alleviating the shortage of water resources in the Southern Xinjiang region, improving the efficiency of water resources, optimizing the structure of the cotton industry in Xinjiang, and facilitating the upgrading of the cotton industry there [11].
When employed in conjunction with experimental studies, crop simulation models have proven to be invaluable tools for integrating a vast array of environmental and management factors that influence crop production. The most common crop models include WOFOST [12,13], APSIM [14,15], DSSAT [16,17], and AquaCrop, a water-driven model based on a conserved relationship between biomass and transpiration. This model is widely used for crop yield simulation due to its fewer input parameters, simplicity of operation, and accuracy of simulation [18,19,20]. To date, the model has been tested in a variety of soils, crops, and climates around the world. The objective of these tests was to determine crops’ response to soil moisture and to facilitate on-farm irrigation management under climate change [17,21,22]. Researchers from around the world have simulated a wide range of crops and have determined the growth parameters for them.
Research has shown that the AquaCrop model can reasonably replicate the growth processes of popular crops, such as cotton [23,24], potato [25], wheat [22], and maize [18], as well as those of vegetables, such as saffron [26], moss [27], and oilseed rape [28]. Furthermore, the model can simulate grain yields at the end of a growing season while accounting for biomass buildup and salt and fertility impacts. In addition, models may be used to determine the best fertility and irrigation schedules for a range of climatic circumstances by analyzing long-term meteorological data series, frequencies, and a variety of scenarios and field management techniques [29,30,31]. Usually, AquaCrop has to be calibrated and validated against regional crop types and environments. The application of AquaCrop to the modeling of deficit irrigation optimization for Gossypium hirsutum L. under long-term climate conditions has been investigated in the dry oasis area of Southern Xinjiang and the humid region of the North Plain in China (NPC) [23,24]. It has been discovered that variations in the model parameters for the same crop can also have an impact on the model’s simulated mistakes in yield and biomass [32]. Because of differences in plant height, leaf shape, and other geographical and temporal variables of plant structure and morphology, cotton displays several genotypes. These traits set it apart from Gossypium hirsutum L., which influences the traits of the canopy of cotton. As a result, the yield of cotton may not be precisely predicted by past studies on Gossypium hirsutum L. variants. Thus, it is essential to look into the characteristics of the cotton model in the dry oasis area near the southern border. This is necessary to close the gap in the knowledge of the AquaCrop simulation parameters for cotton in the vicinity of the desert oasis.
This study set out to achieve three goals: (1) determining irrigation strategies for years with varying precipitation categories; (2) assessing the total soil water, canopy cover, biomass, and crop yields in AquaCrop with different simulated water stresses; (3) using field measurements to calibrate and validate the AquaCrop model for cotton under various irrigation regimes.

2. Materials and Methods

2.1. Site Description

The research was conducted on an experimental farm (76°56′51″ E, 39°12′ 28″ N) in Yuepuhu County, Kashgar, Xinjiang, China (2021–2022). A map of the study area is shown in Figure 1. The soil in the cotton field was classified as silt loam with 0.40% organic matter, 0.30 g kg−1 available nitrogen, 11.5 mg kg−1 available phosphorus, 86 mg kg−1 available potassium, 17.3 g kg−1 total nitrogen, and 10.1 g kg−1 nitrogen. Other relevant soil characteristics (0–100 cm) before seeding are shown in Table S1. In order to ensure uniformity between test plots, soil and fertilizer management was not carried out at the experimental sites prior to the study. The climate of the region is characterized by a temperate, arid continental climate with hot summers and cold winters. The annual sunshine time of the test site is 2825 h, the annual average precipitation is 66.4 mm, and the evaporation is 1600 mm. Detailed information on the temperature, precipitation, and reference evapotranspiration (ET0) in the different growth stages of cotton can be found in Figure S1.

2.2. Design of the Test Pit Experiment and Crop Management

The cotton (Gossypium barbadense. Xin78) was sown on 15 April 2021 and 15 April 2022, and it was harvested on 15 October 2021 and 28 September 2022. A local machine-picked cotton planting mode is adopted. The irrigation method in this experiment was mulched drip irrigation using a single-wing labyrinth irrigator produced by Xinjiang Tianye Water Saving Company (Shihezi, China). The distance between adjacent emitters of drip irrigation tape was 20 cm, the rated working pressure of the emitters was 0.1 MPa, and the design flow rate of the emitters was 3.2 L h−1. Three drip irrigation tapes were arranged under a 2.05 m plastic film, and each tape was used to irrigate two rows of cotton (Figure S1c). The irrigation water source was well water; the buried depth of the well was below 15 m, the salinity of the water was about 0.954 g L−1, and the irrigation water was controlled with a water meter. The irrigation water quotas at the squaring stage (SQS) and the blooming stage (BLS) were based on different proportions of FC (field capacity) and were used as the upper limit of irrigation water and the lower limit of irrigation water. Each irrigation quota was based on formula 1. The crop water consumption (ET) was estimated with the following water balance formula (Formula (2)). Based on the different water requirements at the SQS and the BLS, the drip irrigation test schemes in 2021 and 2022 are shown in Figure S1a,b. The irrigation test plan and actual irrigation conditions are shown in Table 1. Fertilizer dripped into the soil at the root of the cotton along with water. The amount of fertilizer applied was determined according to the local planting method; the fertilizer applied at the SQS was applied at 30% of the total amount, the fertilizer applied at the BLS was applied at 70% of the total amount, and the total fertilizer was applied at 510 kg ha−1 N, 454 kg ha−1 P2O5, and 163 kg ha−1 K2O.
The quota for each drip irrigation event was determined with the following formula:
1000 M = h p (θmaxθmin)
where M is the individual irrigation quota (mm) under drip irrigation; h is the designed wetting depth (cm), which was 40 cm at the cotton budding stage and 60 cm from the flowering stage to the boll opening stage (BOS); p is the moisture ratio, which was 0.7 under drip irrigation; θmax and θmin are the upper and lower limits of soil moisture (volume moisture content) in the designed wetting layers, respectively.
E T   = P   + U   + I W   R     D   10   i = 1 n [ H i   (   θ i 1     θ i 2 ) ]
where ET, P, and U represent the amounts of crop water consumption, precipitation, and groundwater recharge, respectively (mm); IW, R, and D represent the amounts of irrigation, surface runoff, and deep leakage, respectively (mm); i is the number of soil layers; n is the total number of soil layers; Hi is the thickness of soil layer i (cm); θi1 and θi2 are the volumes of soil moisture content at the beginning and end of the calculation interval, respectively (%). Due to the fact that there was no surface runoff and deep leakage in the experimental area because of the limited drip irrigation and precipitation, R = 0 and D = 0.

2.3. Measurement Indicators and Methods

According to the characteristics of cotton growth and canopy development, in the seedling stage (SS), squaring stage (SQS), flowering stage (FS), boll stage (BS), and boll opening stage (BOS), the soil moisture content, above-ground biomass, and leaf area were monitored. The soil moisture content was measured with a Diviner 2000 soil moisture meter the day before irrigation and the day after irrigation; the instrument was calibrated with the drying method. The calibration formula is shown in Formula (3). The depth was measured at 0–100 cm, the 0–40 cm plowing layer was measured in 10 cm intervals, and the 40–100 cm range was measured in 20 cm intervals. The layout of the probe is shown in the planting model diagram in Figure S1b. The soil moisture content was determined by using the average soil moisture content at the point where a common water probe was placed [33]. On each measurement date, five representative plants were taken from each plot to obtain the seasonal aboveground biomass. The samples were killed at 105 °C for 30 min and then dried at 75 °C to a constant weight. In addition, 10 cotton plants were randomly selected in each plot. The leaf area was measured using the punching method [34], where the leaves of individual plants were punched into 40 pieces with a 15 mm punch. The LAI data were then converted into canopy cover (CC) using the LAI and CC relationship formula established by Hsiao et al. [21].
SF = 0.06623 × θ 0.1093
CC = 1.005 (1 − exp (−0.6 LAI))1.2 × 100%
LAI = total   area   of   the   leaf occupied   land   area
Here, SF is the capacitance–frequency reading ratio, and θ is the moisture content in the soil volume (cm3 cm−3).
HI = Ya/B
At the time of harvest, a 6.90 m2 (2.3 m × 3 m) sample of cotton was collected from the center of each plot. The samples were then dried and weighed in order to determine the seed cotton yield. Here, Ya is the seed cotton yield of cotton (t ha−1), B is the aboveground biomass of cotton (t ha−1), and HI (harvest index) is the ratio of the seed cotton yield to aboveground biomass (%).
The yields, amount of irrigation water applied, and consumptive water use were used to compute and analyze WPIrrig and WP. The current prices and production costs were collected under the conditions of a free market, excluding any crop subsidies received, based on prevailing market rates and actual expenses [35].
WPIrrig = Ya/IW
WP = Ya/ET
Here, WPIrrig is the water productivity per amount of irrigation water applied (kg m−3), IW is the amount of irrigation water applied (m3 ha−1), ET is the actual crop evapotranspiration (crop consumptive water use) calculated for the plots (m3 ha−1), WP is the water productivity with respect to the crop consumptive water use (kg m−3), NR (net revenue) was calculated by subtracting the total costs from the total revenue (CNY ha−1), and TC (total cost) included all costs of growing cotton, including labor costs, material costs, water costs, and electricity costs (CNY ha−1).

2.4. Description of the Model (AquaCrop 6.1)

This model establishes a functional connection between CC and ET by disaggregating ETc into components related to transpiration (Tr) and evapotranspiration (E). Transpiration or evapotranspiration coefficients are employed to estimate ET0. The aboveground biomass (B) is estimated using the computed Tr and the standard crop water productivity (WP*). The harvest index value is then utilized to convert B into the final yield (Ya) using the following formulas:
Tr = KcTr,x × CC × ET0
E = Kr (1 − CC) Kex × ET0
B = WP *   ×   ( Tr ET 0 )
Ya = fHI HI0 B
The standard crop transpiration coefficient (KcTr,x) is calculated using the following variables: canopy cover in CC (%), the maximum standard crop transpiration coefficient, and the evapotranspiration reduction coefficient. The Kr coefficient is employed to regulate the impact of surface water deficits. The maximum soil evapotranspiration coefficient is represented by Kex, while the water stress regulating factor is represented by fHI. Finally, the reference harvest index is represented by HI0.
Crop leaf and canopy development, stomatal opening and transpiration, canopy senescence, and the yield index are the main factors affected by water stress. The susceptibility of a crop to water stress is indicated by the water stress coefficient (Ks). The AquaCrop model employs four processes to explain how water stress affects crops: the yield index, stomatal closure, fast canopy senescence, and restricted canopy development. The implications of these processes depend on the sensitivity of the root zone to the total decrease in available soil water. The user may select their desired watering schedule, irrigation strategies, and irrigation schedule using the AquaCrop model. A film module within the site management documentation of the model allows the user to choose various film materials and coverage to replicate film circumstances. In addition to the effects of field surface features and weed control on soil moisture, the model also considers the impact of temperature, salinity stress, and soil fertility on crops. A comprehensive account of the underlying principles, assumptions, and hypotheses of these simulations is provided [19,21,36].

2.5. Model Input, Calibration, and Verification

The experimental field data from the 2021 and 2022 seasons were employed for calibration and verification purposes. A total of 14 groups of treatments were utilized for calibration in 2021, while 10 treatments were employed for verification in 2022. The three principal irrigation growth stages (SQS, FS, and BS) were each associated with 14 distinct irrigation thresholds. The meteorological data required for the model, comprising solar radiation (MJ m−2), maximum temperature (°C), minimum temperature (°C), and rainfall (mm), were provided by a HOBO weather station installed at the test site with an RX3000 data logger (Onset, Bourne, MA, USA). The ET0 value is estimated using the Penman–Monteith formula via AquaCrop’s ET0 calculator. Carbon dioxide concentrations were calculated using the Mauna Loa Observatory’s annual average. Table 2 presents the remaining parameters of the model.

2.6. Model Evaluation Indicators

To assess the model’s performance during calibration and validation, five statistical metrics were utilized, as listed in Formulas (13) to (17). These were the prediction error (Pe), coefficient of determination (R2), root mean square error (RMSE), Nash coefficient (EF), and Willmott’s consistency index (d).
Pe = ( S i     M i ) M i
R 2 = ( i = 1 n   ( S i     S )   ( O i O ) i = 1 n   ( S i S ) 2 i = 1 n   ( O i O ) 2 ) 2
RMSE = 1 n   i = 1 n   ( S i     O i ) 2
EF = i = 1 n   ( O i O ) 2 i = 1 n   ( S i     O i ) 2 i = 1 n   ( O i O ) 2
d = 1 i = 1 n   ( S i     O i ) 2 i = 1 n   (   S i O + O i O ) 2
In this context, n represents the number of observations, Si denotes the simulated value, Oi signifies the measured value, S stands for the simulated mean, and O represents the measured mean. A superior fit is indicated by lower Pe and RMSE values, as well as higher R2 and d values. If EF ≥ 0.80, the model performs well; if EF ranges from 0.60–0.79, the model performs well; if EF ranges from 0.40–0.59, the model performs acceptably; if EF ≤ 0.39, the model performs poorly.

2.7. Irrigation Scenario Design

Following calibration and validation, the model was employed to evaluate cotton growth on Lower Island in diverse field management and historical climate scenarios (Figure S1). The irrigation amount at the seedling stage was 70 mm, and the irrigation quotas at the key growth stage were 360, 405, 450, and 525 mm, respectively. The comprehensive implementation of historical climate scenarios in the AquaCrop model is shown in Table 3. A total of 16 irrigation schemes based on different stages were established for irrigation cycles of varying lengths, including 6, 8, 10, and 12 days. Based on the meteorological data from 1958 to 2022 collected by the Meteorological Data Center of the China Meteorological Administration (http://data.cma.cn (accessed on 20 March 2023)) in Yuephu County, Kashgar, the cumulative rainfall and accumulated ET0 during the cotton growing period were calculated (Figure S2a,b). The Pearson III distribution method was employed for the classification of three precipitation year types [37]. Tsakiris, G. et al. [38] calculated the RDIst value based on cumulative rainfall and cumulative ET0 and used the Survey drought index (RDIst) from 1958 to 2022 to categorize the severity of drought (Figure S2c). The analysis of precipitation categories and drought severity from 1958 to 2022 indicated that 2021 and 2022 were normal and wet years, respectively. Based on the Pearson III precipitation year type and RDIst value, 2004, 2013, 2015, 2016, and 2005 were selected as the wet years for the simulation. The normal years were 1962, 1992, 1963, 1984, and 1988. The simulated dry years were 1961, 1979, 1975, 1970, and 1969.
a k i = j = 1 k   P ij j = 1 k   P ET ij
RDI st i = y i y ¯ σ y
Here, P and PET are the precipitation and potential evapotranspiration of the j-th month of the i-th year, yi is lnak, y ¯ i is its arithmetic mean, and σy is its standard deviation; k is the number of growing seasons: 64.

2.8. Data Analysis

The data obtained in this experiment were processed using Excel 2021 (Microsoft Corp, Redmond, WA, USA). A planting pattern diagram was created using Microsoft PowerPoint 2021 (Microsoft Corp, Redmond, WA, USA). The results of each treatment were analyzed using SPSS statistics 25.0 (SPSS Inc. IBM Corp., Armonk, NY, USA). Graphs were generated using Origin 9.0 (Origin Lab Corp., Northampton, MA, USA).

3. Results

3.1. Total Soil Water And ETc

The simulated and measured values of TSW in 2021 were employed for parameter calibration (Figure 2), and validation was carried out using the simulated and measured values of TSW for different treatments in 2022 (Figure 3). Overall, the simulated and observed TSW values for all treatments in 2021 demonstrated a strong correlation with precipitation and irrigation (Pe = −12.37%–−6.82%, R2 = 0.27–0.72, RMSE = 7.88–27.75, EF = −0.34–−0.64, d = 0.60–0.91). In the model validation using the data from 2022, the model overestimated the TSW content in the middle and later stages of growth, and the indicators of goodness of fit (i.e., Pe, R2, RMSE, EF, d) were −22.75–13.99, 0.27–0.80, 3.98–21.95 mm, 0.11–0.75, and 0.55–0.94, respectively. During the early reproductive phase, there was an overall trend of the simulated values of TSW being underestimated. The discrepancy between the simulated and actual values of TSW gradually increased as the lower irrigation limit was approached when the top irrigation limit was 90% FC. For the extremely water-stressed WD1 and overwatered CK treatments, the discrepancies between the simulated and measured TSW values were even more pronounced, with WD1 exhibiting Pe, R2, RMSE, EF, and d values of −5.14–−7.97, 46.49%, 10.26 mm, 0.32, and 82.7, respectively. The Pe, R2, RMSE, EF, and d values for CK were −22.75–−12.43, 27.03%, and 21.95%, respectively. The Pe, R2, RMSE, EF, and d values for CK were 27.03%, 21.95 mm, 0.24, and 0.55, respectively. A comparison of the observed and simulated values of total ETc for the two years of the growing season of cotton is presented in Table 4. The simulated values of ET were greater than the observed values, and Pe was greater for all treatments, with a range of −10.06–2.47%. The simulated values of total ETc for these two years displayed a range of −10.06–2.47%, with an R2 value of 0.54, RMSE of 25.30 mm, EF of −0.22, and d of 0.77. These results indicate that the simulation of ETc was not satisfactory.

3.2. Canopy Cover

Figure 4 and Figure 5 illustrate the growth process of cotton. The AquaCrop model accurately simulated the development of the calibrated CC for each treatment. As indicated by the high R2 (≥0.97), the high d value (≥0.99), the high EF value (≥0.95), and the low estimation error (−6.86% ≤ Pe ≤ 26.17%, 3.78% ≤ RMSE ≤ 7.08%), the model demonstrated a high degree of accuracy. The maximum canopy cover (CC) that could be achieved with the 14 irrigation treatments for cotton decreased with the increase in water stress applied. Across all treatments, the impact of water deficiency stress on cotton growth was greater in 2022 than in 2021. The validation phase in 2022 revealed higher estimation errors than those observed during calibration, particularly in treatments with significant water deficits. In the 2022 validation, the ranges for Pe, R2, RMSE, EF, and d for canopy cover across treatments were 15.90%–56.80%, 0.93–1.00, 2.50–11.98%, 0.82–0.99, and 0.96–1.00, respectively. The experimental region was located in an arid inland zone at the southern border, where the WD1, WD4, and WD7 treatments demonstrated an inadequate lower limit of irrigation at the BLS, resulting in an acceleration of leaf senescence. In the validation phase, the model demonstrated a tendency to overestimate the development and senescence of CC. The larger error observed in the WD1 treatment was primarily attributed to a low lower limit of soil moisture, resulting in the model’s overestimation of canopy cover in the mid-to-late growth stage of cotton. Conversely, the irrigation floors of 65% FC and 75% FC were found to be highly accurate. Excessive water stress at the BLS was identified as a factor contributing to the decreased accuracy of the simulation of the canopy because the temperatures observed in the southern region of the border were higher than the upper limit of temperature for the growth of cotton (Tupper).

3.3. Aboveground Biomass

Figure 6 illustrates the results of the calibration of the biomass accumulation process for cotton. Figure 7 presents the results of the validation of the biomass accumulation process for cotton. The AquaCrop model estimated the progress of aboveground biomass with satisfactory accuracy under the different irrigation treatments. The difference between the simulated and measured values decreased with the lower irrigation limit in the treatment with severe water stress. In the calibration phase, the values of Pe, R2, RMSE, EF, and d for all treatments were within the ranges of −39.74–21.41%, 0.97–1.00, 1.53 t ha−1–3.34 t ha−1, and 0.94–1.00, respectively. The measured values were found to be closer to the simulated values in the validation phase, with the simulation effect becoming more credible. The corresponding value of the goodness-of-fit index was found to be −24.86%–−24.41%, with a value of 0.97. The model displayed a tendency to underestimate the observed values, particularly after BS, with the following values: −0.98, 1.91 t ha−1–4.36 t ha−1, 0.87–0.98, and 0.97–1.00, respectively. This overestimation was particularly evident in the water-stressed treatments.

3.4. Biomass and Final Seed Yield

Figure 8 presents the outcomes of the model calibration and validation for the biomass and seed cotton yield. During the calibration process, the Pe, R2, RMSE, EF, and d values for the simulated and measured aboveground biomass in all treatments were found to be −39.74–24.41%, 0.98, 1.23 t ha−1, 0.99, and 0.99, respectively. During verification, the Pe, R2, RMSE, EF, and d values were found to be −24.86%–24.41%, 0.96, 1.33 t ha−1, 0.95, and 0.99, respectively. The simulation error for the biomass in the WD1, WD4, and WD7 treatments was considerable. The Pe, R2, RMSE, EF, and d values between the simulated and measured values were −0.28–3.77%, 0.93, 0.15 t ha−1, −0.91, and 0.75, respectively. The Pe, R2, RMSE, EF, and d values were found to be in the ranges of −0.09–4.45%, 0.81–0.14 t ha−1, 0.54, and 0.84. At the lower limit of 75% FC irrigation, the model was capable of accurately simulating biomass and yield. The yield of cotton displayed a positive correlation with the increase in the irrigation amount, reaching a peak at the full irrigation level. However, when the irrigation amount exceeded this threshold, the yield displayed a decline, which was potentially due to the adverse effects of over-irrigation.

3.5. Model Application

Following the calibration and validation of the model, 16 irrigation combinations were selected for simulation. These were imported from meteorological data from 1958 to 2022, with the basic parameters of the model remaining unchanged (Table 2). The three scenario years were based on five typical years and 16 irrigation treatments, resulting in a total of 240 scenarios. The simulations of seed cotton yield, final aboveground biomass, ET, WP, and WPIrrig in different irrigation scenarios are presented in Table 5. A box diagram of the specific trends in the three scenario years is shown in Figure S3.
The average seed cotton yield in wet, normal, and dry years was 6.35 t ha−1, 5.34 t ha−1, and 4.98 t ha−1, respectively. When the irrigation cycle was certain, the seed cotton yields in the three scenario years increased with the increase in irrigation water. When the amount of irrigation water was certain, an irrigation period of six days resulted in a yield increase of 0.01–5.02% compared with the other treatments. When the irrigation period was constant with 525 mm of irrigation water, the yield increased by 0.72–20.18% compared with the other treatments. In years with high precipitation, for instance, when the volume of irrigation water was fixed, the biomass, evapotranspiration (ET), water use efficiency (WP), and irrigation water use efficiency (WPIrrig) of cotton displayed a tendency to decrease with the increase in the irrigation cycle. In the case of a six-day irrigation cycle, the biomass, ET, WP, and WPIrrig were observed to be higher than those of the other treatments by 0.73–5.61%. The biomass reached a maximum of 17.95 t ha−1 in the case of an irrigation water volume of 525 mm and an irrigation period of 12 d (T13). The maximum ET was 562 mm at 525 mm, with an irrigation period of 12 d (T13). The maximum WP was 1.27 kg m−3 at 525 mm, with an irrigation period of 12 d (T4). The maximum WP was 1.27 kg. At 525 mm, 6 d (T4), and 6 d (T4), the maximum WP was 0.73–5.61%, 0.86–6.47%, 0.32–3.28%, and 0.17–3.32%, respectively. The maximum WPIrrig was observed at 525 mm and 6 d (T1), with a value of 1.46 kg m−3.
The total costs and net revenue corresponding to the simulated treatment were calculated based on the cotton prices and local prices of agricultural production materials in 2021 and 2022, as shown in Figure 9. The blue bar chart shows the total costs in wet years; the green bar chart shows the total costs in normal years; the red bar chart shows the total costs in dry years. Each subgraph’s colored bar chart depicts the net revenue of each simulation process. It was determined that the dry years had a discernible impact on the net yield of cotton. In dry years, the cotton yield was demonstrably lower than that in wet years or normal years. Furthermore, the difference in yield among all treatments in dry years was statistically significant. When the irrigation cycle was held constant, the net revenue of cotton gradually increased with the increase in the irrigation amount in dry years. When the amount of irrigation was held constant, the net revenue of cotton displayed a decreasing trend with increasing irrigation cycle length, reaching a maximum of 2.48 × 104 CNY ha−1 in the T15 treatment.

3.6. A Comprehensive Evaluation of the AquaCrop Simulation Results Based on Different Irrigation Scenarios

In order to elucidate the underlying drivers of changes in seed cotton yield, correlation analyses were conducted on various agronomic and economic variables, including WPIrrig, WP, net revenue (NR), HI, total cost (TC), total irrigation (TI), irrigation frequency (IF), final aboveground biomass (Bio), and ET. The results of these analyses are presented in Figure 10. A significant correlation was found among Ya, Bio, TI, WPIrrig, WP, ET, and TC for cotton in the three scenario years, as well as between wet and normal years. Additionally, a significant correlation was found between seed cotton yield and Ya, Bio, TI, WPIrrig, WP, HI, ET, and TC. Furthermore, a strong correlation was observed among Ya, Bio, TI, ET, and TC (R2 > 0.88) for cotton in wet and normal years. However, efficiency indicators such as WPIrrig, WP, and HI demonstrated a negative correlation (R2 > −0.76). Figure 10 illustrates the relative importance of each factor for the yield of cotton, with TC being identified as the most significant factor, followed by Bio, ET, TI, and IF. It is noteworthy that the TI in dry years was greater than that in normal and wet years, which was likely due to the higher evaporation rates in dry years, which are also characterized by higher temperatures. The correlation analysis revealed that TI and TC displayed a highly significant positive correlation with seed cotton yield during normal and dry years. Additionally, TI and Ya demonstrated a highly significant positive correlation, while NR displayed a highly significant positive correlation with TC. It can be concluded that cotton is more dependent on TI, and ensuring the amount of water allotted to cotton during normal and dry years is crucial for stabilizing yield. In addition, the primary drivers of seed cotton yield were ET, TI, and TC. This may be attributed to the fact that seed cotton yield is the result of a process that couples both irrigation management and anthropogenic management and is influenced by a combination of physiological and environmental factors.
Principal component analysis (PCA) was employed to conduct a comprehensive evaluation of various indexes pertaining to cotton [39]. The first two components were deemed valid for the three scenarios, as their sum exceeded two-thirds of the total variance. In wet years, for instance, the first two components accounted for 97.95% of the total variation. PC1 showed a positive correlation with yield, Bio, ET, and TI, as well as a negative correlation with WPIrrig, WP, and HI. In PC2, which accounted for 12.12% of the total variance, IF, WPIrrig, and NR displayed a robust positive correlation, whereas TI, WPIrrig, and TC demonstrated a negative correlation (Table 6).
It was reported that PCA was used to obtain an important relationship between the yield and water amount [40,41]. This study used the first two components (PC1 and PC2) to create biplots of wet, normal, and dry years (Figure 11), explaining the influence of different scenario years on cotton yield, efficiency, moisture, and investment indexes. Upon evaluating the positive aspects of PC1 and PC2, a significant positive correlation was identified between Ya and Bio (Figure 11), particularly in years with low precipitation; here, maintaining a sufficient amount of final aboveground biomass played a pivotal role in the formation of seed cotton yield due to the high soil evaporation. This phenomenon may be explained by the strong correlation of TC, NR, and TI in dry seasons with seed cotton yield. Due to the high evapotranspiration of the crop itself, increasing the irrigation quota and decreasing the watering cycle can minimize the negative impacts on cotton yield due to climatic extremes. In contrast, Ya, TC, and NR were strongly correlated in wet and normal years. As long as the TI meets the normal growth and development of the crop and mild stress does not occur, the higher the investment, the higher the seed cotton yield. The results of the principal component analysis were consistent with those of the correlation analysis.
Aboveground biomass is indicative of the growth of cotton. Factors influencing agricultural irrigation management, including evapotranspiration (ET), total irrigation, and the number of irrigations, directly impact the seed cotton yield. Higher seed cotton yield necessitates human inputs, prompting concern regarding the investment and net revenue of economic factors for agricultural producers. Southern Xinjiang, a region where cotton is cultivated, is water-scarce. Therefore, improving water use efficiency and the harvest index is a crucial aspect of efficient water use in arid agriculture. In this study, a raw matrix was constructed based on simulated data from the 16 treatments in the three scenario years. The matrix was constructed using 10 indexes, namely, Bio (X1), Ya (X2), ET (X3), TI (X4), IF (X5), TC (X6), NR (X7), WPIrrig (X8), WP (X9), and HI (X10). The raw data were initially standardized to eliminate the effect of magnitude (Figure 12). Subsequent to this, the comprehensive score and ranking of each treatment were obtained by substituting the standardized data into a linear function, as shown in Table 7.
In the wet year, the three treatments with the highest comprehensive evaluation scores were T13, T14, and T15. In the normal year, the three treatments with the highest comprehensive evaluation scores were T13, T14, and T15. In the dry year, the three treatments with the highest comprehensive evaluation scores were T15, T13, and T14. However, the implementation of water resource regulation measures in Southern Xinjiang reached the maximum allowable precipitation of 450 mm during the critical growth period. In light of the aforementioned considerations, the cotton irrigation quota of 520 mm was treated according to T9 in both the wet and dry years in order to satisfy the water regulations and ensure stable yield. The critical period of growth was 450 mm, the irrigation cycle was 6 days/time, and the number of irrigations was 16 as per the final plan. In consideration of the stabilization of the crop yield in dry years, the final plan was to treat the cotton in accordance with T15, with an irrigation quota of 595 mm, a critical growth period of 525 mm, an irrigation cycle of 10 days/time, and 10 irrigations.
Wet year:
The first principal component was w1 = 0.115 X1 + 0.114 X2 + 0.116 X3 + 0.11 X4 + 0.034 X5 + 0.116 X6 + 0.111 X7 − 0.107 X8 − 0.114 X9 − 0.116 X10.
The second principal component was w2 = −0.019 X1 + 0.059 X2 + 0.099 X3 − 0.252 X4 + 0.783 X5 − 0.09 X6 + 0.122 X7 + 0.306 X8 − 0.136 X9 + 0.005 X10.
Normal year:
The first principal component was n1 = 0.118 X1 + 0.116 X2 + 0.117 X3 + 0.113 X4 + 0.032 X5 + 0.116 X6 + 0.116 X7 − 0.107 X8 −0.106 X9 − 0.115 X10.
The second principal component was n2 = 0.003 X1 − 0.058 X2 + 0.057 X3 − 0.1974 X4 + 0.703 X5 − 0.104 X6 − 0.043 X7 + 0.277 X8 − 0.293 X9 − 0.132 X10.
Dry year:
The first principal component was d1 = 0.132 X1 + 0.129 X2 + 0.131 X3 + 0.126 X4 + 0.039 X5 + 0.128 X6 + 0.128 X7 − 0.116 X8 − 0.074 X9 − 0.108 X10.
The second principal component was d2 = 0.012 X1 + 0.087 X2 − 0.025 X3 + 0.131 X4 − 0.423 X5 + 0.083 X6 + 0.099 X7 − 0.181 X8 + 0.368 X9 + 0.247 X10.

4. Discussion

4.1. Applicability of the AquaCrop Model

In this study, AquaCrop was used to simulate the soil water dynamics in the cotton growing season in Southern Xinjiang. Similarly to this study, other studies have also reported that AquaCrop has satisfactory performance in estimating soil moisture content in Gossypium hirsutum L. with different types of irrigation management [20]. AquaCrop made good predictions of soil water dynamics and crop response in all treatments of cotton; for example, the TSW was accurately simulated for the treatment with a lower irrigation limit of 65% FC. In addition, some studies have shown that AquaCrop has limitations in simulating severe water stress [20,42,43], but it also underestimated the TSW (−12.37% ≤ Pe ≤ 6.82%, 0.27 ≤ R2 ≤ 0.72, 7.88 mm ≤ RMSE ≤ 27.75 mm) for the treatment with a lower irrigation limit of 75% FC (−0.34 ≤ EF ≤ 0.64, 0.60 ≤ d ≤ 0.91). Because there is only so much water in a single irrigation cycle [36], the model may have trouble simulating a wet soil surface at a lower irrigation limit of 75% FC. This is because the model has to copy the complex water movement in the surface soil under drip irrigation. Secondly, the variations in the canopy of cotton may be different from those of Gossypium hirsutum L. due to differences in the variety, which further affected the variation rule for soil water evaporation in cotton [44,45]. Ultimately, the oversimplification of the model’s internal components, including soil water balance, root distribution pattern, canopy development and senescence, and biomass production and distribution, could be the cause [20,25,44].
The model poorly simulated ET (−10.06% ≤ Pe ≤ 2.43%), and it underestimated most of the ET in most of the treatments, which was related to the use of unadjusted or low evapotranspiration and evaporation coefficients [20]. On the one hand, a possible reason is that frequent irrigation with a small amount of water causes the evaporation of water from holes in membranes to increase, which increases the error of ETc and causes it to be underestimated. Conversely, stress only affects the canopy cover curve during the vegetative stage [46], leading to a replacement of the canopy cover curve with the classical crop coefficient curve for calculating crop transpiration and soil evaporation, thereby reducing the simulation accuracy of ET [20,45,46].
For the CC simulation results, R2 (≥0.97), EF (≥0.95), and d (≥0.99) were relatively high, while RMSE (3.78%–7.08%) and Pe (−6.86%–26.17%) were relatively low. In general, the AquaCrop model was able to more accurately simulate cotton’s CC growth process. However, for the treatment with a lower limit of 55% FC at the BLS, the CC simulation was poor, indicating that the model could not fully simulate the growth of CC under severe water stress. The main reason may be that canopy senescence is accelerated under water stress; senescence occurred in the simulated CC, but the growth rate of the measured CC was slowed down when water stress occurred in a short period of time. Zhang et al. [24] found that the model could effectively simulate the CC growth of cotton under brine irrigation (0.69 ≤ R2 ≤ 0.99) (4.83% ≤ RMSE ≤ 23.42%). Both Feng et al. and Ran et al. reported a similar trend. The reason may be that the CC curve in AquaCrop is not sensitive to water stress throughout the growing season but is more sensitive to water stress in the vegetative stage, resulting in lower simulated CC values in the middle and late stages [42,45]. In 2012, Tan et al. [20] highlighted the model’s limitations in simulating CC under severe water stress in a study of Gossypium hirsutum L. The simulation of CC under the condition of a 55% irrigation limit in the T4 treatment was found to be underestimated, while the CC was overestimated in this experiment. This discrepancy may be attributed to the differences in plant types. The development process of the canopy in the early and middle periods is slower than that of Gossypium hirsutum L. Additionally, the time required to reach CCmax differed from that of Gossypium hirsutum L., which affected the evapotranspiration process of soil water.
In the biomass simulation, crop transpiration was converted into biomass using a normalized water production efficiency (WP*) of 24 g m−2, which differed from the WP* values found in prior research using AquaCrop in terrestrial cotton. This discrepancy was likely due to the potential influence of varietal differences in WP*. Zhang et al. and Wang et al. [24,47] simulated reductions in biomass production of 20 g m−2 and 19 g m−2 due to the WP* of salinity stress and a filmless mulching strategy, respectively. The model was used to simulate the final biomass of Gossypium hirsutum L. at the end of the fertility period with different salinity gradients and moisture stresses in filmless mulching. The results showed that Zhang overestimated the biomass, and Wang underestimated it. It is evident that different treatments can cause differences in model calibration.
In our study, the aboveground biomass in the treatment with a lower limit of irrigation of 65% FC was underestimated in the middle of the simulation and overestimated in the BS. This discrepancy may be attributed to the utilization of a constant WP* value throughout the growth period, though it should gradually decline as the plant matures. Another possible explanation is that the model underestimated the early leaf senescence due to water stress, leading to an overestimation of the canopy and transpiration. The leaves of cotton subjected to severe water stress were removed through the application of defoliants at the maturity stage. However, the treatment with a lower irrigation limit of 65% FC at the BLS extended the growth period of cotton to a certain extent, thereby mitigating the impact of defoliants. This finding is consistent with the results of Feng et al. and Steduto et al. [18,45]. According to previous studies, due to the physiological regulation and compensatory growth of cotton, mild water stress is beneficial for the growth of cotton, and the BLS is more sensitive to soil water stress than other stages [48,49]. Overall, AquaCrop is capable of adequately simulating the differences in yield and aboveground biomass quality due to varietal differences for cotton treatments under mild moisture stress.

4.2. Application of the Model to Different Irrigation Scenarios

Crop modeling and effective field irrigation decisions may address the impact of climate change on the cotton industry in Xinjiang [50,51]. In the scenario analysis, all of the common scenarios with different irrigation quotas were simulated and analyzed for different scenario years. In order to achieve optimal crop yield and irrigation water productivity, it is essential to utilize appropriate irrigation quotas in normal years to attain the highest possible yield. As evidenced by the field experiment conducted in the mulched cotton cultivation of Aksu Oasis in Southern Xinjiang, the 75% FC field irrigation method demonstrated the highest yield and water use efficiency, a result that aligned with the findings of our modeling analysis [52]. The application of AquaCrop also demonstrated that, under mulched drip irrigation, the irrigation water quota was 45 mm during the SQS, FS, and BS. When there were 11 irrigations, the yield was 6186.17 kg ha−1, and the economic benefit was the largest [23].
In our scenario analysis, the irrigation system was optimized based on the critical period of cotton growth, and the simulation of different precipitation categories was realized. These appropriate water stress treatments may be effective measures to stabilize yield and improve water use efficiency. Cotton yield is related not only to varieties [32] but also to climatic factors [53,54]. Only two years of field trials served as the basis for the calibration and validation of AquaCrop due to data limitations. The model gave satisfactory estimates of the soil water dynamics, growth, and yield of cotton in normal and wet years. Nevertheless, the specific weather conditions in the single validation year may have affected the model performance and scenario analysis to some extent. In AquaCrop, too low of a soil moisture limit could seriously affect the simulation. In many previous studies, significant biases in the AquaCrop model were observed when simulating severe and extreme water stress [17,20,28,42,44]. In the scenario analysis, with the exception of 1978, 1990, and 1994, which were characterized by extreme drought, and 1959, 1973, and 1997, which were characterized by severe drought, the remaining years ranged from wet years to years with moderate drought conditions. Therefore, only under conditions of moderate aridity or humidity can we use the calibrated AquaCrop model with reliable accuracy. It is our contention that the AquaCrop model can be employed as an effective tool for providing information on the cultivation of island cotton under deficit irrigation, with the additional benefit of a supplementary study to verify the measured data of island cotton in drought years in Southern Xinjiang. Huang et al. [22] demonstrated that biochar has the potential to enhance soil texture, thereby augmenting cotton’s capacity to retain soil moisture. It is, therefore, necessary to adjust the soil texture by controlling the dosage and timing of biochar application in the future in order to improve the water use efficiency and fertilizer use efficiency of cotton in northwestern arid areas. In light of the water control and water restriction policies currently in place in Southern Xinjiang, as well as the associated economic costs, there is potential for improvement in the deficit irrigation system through the integration of an economic model.

4.3. Selection and Evaluation of Optimal Moisture Regulation for Machine-Picked Cotton in Southern Xinjiang

PCA is an unsupervised multivariate analysis method [55]. Many researchers have used PCA for data analysis [40,56,57]. For instance, Jarwar et al. [56] employed PCA to identify the Chinese Gossypium hirsutum L. varieties with the highest scores in terms of fiber and yield. The results of the principal component analysis (PCA) indicated that when combining various factors, including growth factors, yield factors, agricultural irrigation management factors, economic factors, efficiency factors, and local water resources in Southern Xinjiang, the irrigation quota for cotton in both wet and dry years was found to be 520 mm. Additionally, the critical growth period was determined to be 450 mm, the irrigation cycle was once every six days, and irrigation occurred 16 times as per the final plan. In consideration of the stabilization of crop yield in dry years, the final plan was to treat island cotton according to T15, with an irrigation quota of 595 mm, a critical growth period of 525 mm, an irrigation cycle of 10 days, and an irrigation frequency of 10 times. In addition, a strict method was used to choose a wide range of meteorological data based on the AquaCrop model in order to find the best irrigation plans for island cotton in several different scenario years. This study also underscored the significance of water regulation on the yield, benefit, and water use efficiency of island cotton in desert oasis areas of Xinjiang.

5. Conclusions

The AquaCrop model was used to simulate the growth performance of cotton (Gossypium barbadense) under mulched drip irrigation in Southern Xinjiang, China. The model was calibrated and validated at the local level based on two years of field data. When different water stress conditions were tested and compared, it was found that the model had a significant simulation effect on TSW, ETc, CC, Bio, and Ya, with an RMSE of less than 11.58%. The TSW content was found to be overestimated in the middle and late periods of cotton growth, with the RMSE value ranging from 3.98 mm to 27.75 mm. The simulated evapotranspiration (ETc) values were generally high, but the prediction error for each treatment was less than 11%. In the calibration year (2021), the model demonstrated satisfactory performance in simulating aboveground biomass and yield, with root mean square error (RMSE) values of 1.23 t ha−1 and 0.15 t ha−1, respectively. The simulation performance decreased in the verification year (2022), but the prediction error was less than 15%. AquaCrop can be a reliable tool for guiding the field management for irrigated cotton in Southern Xinjiang.
Following the verification of the dataset, which comprised 64 years of meteorological data, principal component analysis was employed to recommend cotton irrigation schemes in Southern Xinjiang, and three scenarios were presented. In wet years, the irrigation quota was 520 mm, and the critical period of growth irrigation was once every 6 days; in normal years, the irrigation quota was 520 mm, and the critical period of growth irrigation was once every 6 days. In dry years, an irrigation quota of 595 mm and irrigation once every 10 days in the critical growth period can provide higher seed cotton yield and irrigation water productivity so as to maximize cotton yield and net income in Southern Xinjiang. In addition, the model was only calibrated and validated on the basis of a two-year field trial under the weather conditions of normal and wet years. To confirm these conclusions, future confirmatory studies using long-term field data with large climatic variability will be necessary.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/agriculture14081269/s1. Figure S1. This experiment had irrigation treatment schedules for 2021 and 2022 (S1-a, S1-b). The most widely used method of cotton cultivation in Xinjiang is drip irrigation under film (S1-c). The temperature, ET0, and rainfall (d) for two years during the trial are provided. Figure S2. Changes in precipitation (S2-a), reference crop evapotranspiration (ET0) (S2-b), the drought index, and drought severity (S2-c) in the growing season of cotton in Yuepuhu County, Kashgar from 1958 to 2022. Figure S3. Comparison of the final yield, biomass, ET, WP, and WPIrrig of cotton in different comprehensive scenarios in wet years, normal years, and dry years. Box plots of simulated trends for each irrigation scenario in the three scenario years (S3-a–S3-o). Table S1. Main physical properties of soil at a depth of 0–100 cm.

Author Contributions

Conceptualization, Y.D.; methodology, Y.D.; software, Y.D. and Q.F.; validation, Y.D. and Q.F.; formal analysis, Y.D. and Y.P.; investigation, Y.D.; resources, Q.F. and Y.M.; data curation, Y.D.; writing—original draft preparation, Y.D.; writing—review and editing, Y.D., Q.F. and Y.M.; visualization, Y.D.; supervision, Q.F., Y.M. and P.A.; project administration, Q.F. and Y.M.; funding acquisition, Q.F. and Y.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key R&D Project of Xinjiang Uygur Autonomous Region (No. 2022B02009−3), the Major Science and Technology Project of the Xinjiang Uygur Autonomous Region (No. 2022A02003−6), the Major Science and Technology Project of the Xinjiang Uygur Autonomous Region (No. 2022A02011−2), the Earmarked Fund for the China Agriculture Research System (CARS−15−13), and the Key Talent Training Project “Three Rural Areas” in Xinjiang (No. 2022SNGGGCC016).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

This paper contains the original contributions made during the study. For more information, contact the corresponding author.

Acknowledgments

The authors would like to thank the corresponding editor and the anonymous reviewers for their insightful criticism and suggestions on how to improve the paper.

Conflicts of Interest

The authors declare that none of the work presented in this study may have been influenced by any known conflicting financial interests or personal relationships.

References

  1. Asseng, S.; Ewert, F.; Martre, P.; Rötter, R.P.; Lobell, D.B.; Cammarano, D.; Kimball, B.A.; Ottman, M.J.; Wall, G.W.; White, J.W.; et al. Rising temperatures reduce global wheat production. Nat. Clim. Chang. 2015, 5, 143–147. [Google Scholar] [CrossRef]
  2. Lesk, C.; Coffel, E.; Horton, R. Net benefits to US soy and maize yields from intensifying hourly rainfall. Nat. Clim. Chang. 2020, 10, 819–822. [Google Scholar] [CrossRef]
  3. Chavez, E.; Conway, G.; Ghil, M.; Sadler, M. An end-to-end assessment of extreme weather impacts on food security. Nat. Clim. Chang. 2015, 5, 997–1001. [Google Scholar] [CrossRef]
  4. Li, Y.; Guan, K.; Schnitkey, G.D.; DeLucia, E.; Peng, B. Excessive rainfall leads to maize yield loss of a comparable magnitude to extreme drought in the United States. Glob. Chang. Biol. 2019, 25, 2325–2337. [Google Scholar] [CrossRef]
  5. Li, Z.; Zou, H.; Lai, Z.; Zhang, F.; Fan, J. Optimal drip fertigation regimes improved soil micro-environment, root growth and grain yield of spring maize in arid northwest china. Agronomy 2023, 13, 227. [Google Scholar] [CrossRef]
  6. Hu, Y.; Chen, J.; Fang, L.; Zhang, Z.; Ma, W.; Niu, Y.; Ju, L.; Deng, J.; Zhao, T.; Lian, J.; et al. Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton. Nat. Genet. 2019, 51, 739–748. [Google Scholar] [CrossRef]
  7. Yang, Z.; Qanmber, G.; Wang, Z.; Yang, Z.; Li, F. Gossypium genomics: Trends, scope, and utilization for cotton improvement. Trends Plant Sci. 2020, 25, 488–500. [Google Scholar] [CrossRef]
  8. Allanov, K.; Shamsiev, A.; Durdiev, N.; Avliyakulov, M.; Karimov, A.; Khaitov, B. Improving nutrition and water use efficiencies of pima cotton (Gossypium barbadense L.) varieties under arid conditions of Uzbekistan. J. Plant. Nutr. 2020, 43, 2590–2600. [Google Scholar] [CrossRef]
  9. Zhu, Y.; Sun, L.; Luo, Q.; Chen, H.; Yang, Y. Spatial optimization of cotton cultivation in Xinjiang: A climate change perspective. Int. J. Appl. Earth Obs. Geoinf. 2023, 124, 103523. [Google Scholar] [CrossRef]
  10. Piao, S.; Ciais, P.; Huang, Y.; Shen, Z.; Peng, S.; Li, J.; Zhou, L.; Liu, H.; Ma, Y.; Ding, Y.; et al. The impacts of climate change on water resources and agriculture in China. Nature 2010, 467, 43–51. [Google Scholar] [CrossRef]
  11. Yang, Z.; Tang, J.; Yu, M.; Zhang, Y.; Abbas, A.; Wang, S.; Bagadeem, S. Sustainable Cotton Production through Increased Competitiveness: Analysis of Comparative Advantage and Influencing Factors of Cotton Production in Xinjiang, China. Agronomy 2022, 12, 2239. [Google Scholar] [CrossRef]
  12. Akbari Variani, H.; Afshar, A.; Vahabzadeh, M.; Molajou, A. A review on food subsystem simulation models for the water-food-energy nexus: Development perspective. Environ. Sci. Pollut. Res. 2023, 30, 95197–95214. [Google Scholar] [CrossRef] [PubMed]
  13. Berghuijs, H.N.; Silva, J.V.; Reidsma, P.; de Wit, A.J. Expanding the WOFOST crop model to explore options for sustainable nitrogen management: A study for winter wheat in the Netherlands. Eur. J. Agron. 2024, 154, 127099. [Google Scholar] [CrossRef]
  14. Lagerquist, E.; Vogeler, I.; Kumar, U.; Bergkvist, G.; Lana, M.; Watson, C.A.; Parsons, D. Assessing the effect of intercropped leguminous service crops on main crops and soil processes using APSIM NG. Agric. Syst. 2024, 216, 103884. [Google Scholar] [CrossRef]
  15. Shahid, M.R.; Wakeel, A.; Ullah, M.S.; Gaydon, D.S. Identifying changes to key APSIM-wheat constants to sensibly simulate high temperature crop response in Pakistan. Field Crop. Res. 2024, 307, 109265. [Google Scholar] [CrossRef]
  16. Guhan, V.; Annadurai, K.; Easwaran, S.; Marimuthu, M.; Balu, D.; Vigneswaran, S.; Navinkumar, C. Assessing the impact of climate change on water requirement and yield of sugarcane over different agro-climatic zones of Tamil Nadu. Sci. Rep. 2024, 14, 8239. [Google Scholar] [CrossRef]
  17. Wang, H.; Cheng, M.; Liao, Z.; Guo, J.; Zhang, F.; Fan, J.; Feng, H.; Yang, Q.; Wu, L.; Wang, X. Performance evaluation of AquaCrop and DSSAT-SUBSTOR-Potato models in simulating potato growth, yield and water productivity under various drip fertigation regimes. Agric. Water Manag. 2023, 276, 108076. [Google Scholar] [CrossRef]
  18. Ran, H.; Kang, S.; Hu, X.; Li, F.; Du, T.; Tong, L.; Li, S.; Ding, R.; Zhou, Z.; Parsons, D. Newly developed water productivity and harvest index models for maize in an arid region. Field Crop. Res. 2019, 234, 73–86. [Google Scholar] [CrossRef]
  19. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO crop model to simulate yield response to water: I. Concepts and underlying principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef]
  20. Tan, S.; Wang, Q.; Zhang, J.; Chen, Y.; Shan, Y.; Xu, D. Performance of AquaCrop model for cotton growth simulation under film-mulched drip irrigation in southern Xinjiang, China. Agric. Water Manag. 2018, 196, 99–113. [Google Scholar] [CrossRef]
  21. Hsiao, T.C.; Heng, L.; Steduto, P.; Rojas Lara, B.; Raes, D.; Fereres, E. AquaCrop—The FAO crop model to simulate yield response to water: III. Parameterization and testing for maize. Agron. J. 2009, 101, 448–459. [Google Scholar] [CrossRef]
  22. Huang, M.; Wang, C.; Qi, W.; Zhang, Z.; Xu, H. Modelling the integrated strategies of deficit irrigation, nitrogen fertilization, and biochar addition for winter wheat by AquaCrop based on a two-year field study. Field Crop. Res. 2022, 282, 108510. [Google Scholar] [CrossRef]
  23. Wang, H.; Li, Z.; Si, Z.; Hamani, A.K.M.; Huang, W.; Fan, K.; Wang, X.; Gao, Y. Alternative planting patterns of film-mulching cotton for alleviating plastic residue pollution in Aksu oasis, southern Xinjiang. Ind. Crop. Prod. 2023, 203, 117205. [Google Scholar] [CrossRef]
  24. Zhang, J.; Li, K.; Gao, Y.; Feng, D.; Zheng, C.; Cao, C.; Sun, J.; Dang, H.; Hamani, A.K.M. Evaluation of saline water irrigation on cotton growth and yield using the AquaCrop crop simulation model. Agric. Water Manag. 2022, 261, 107355. [Google Scholar] [CrossRef]
  25. Yin, J.; Yang, Y.; Eeswaran, R.; Yang, Z.; Ma, Z.; Sun, F. Irrigation scheduling for potatoes (Solanum tuberosum L.) under drip irrigation in an arid region using AquaCrop model. Front. Plant Sci. 2023, 14, 1242074. [Google Scholar] [CrossRef] [PubMed]
  26. Mirsafi, Z.; Sepaskhah, A.R.; Ahmadi, S.H.; Kamgar-Haghighi, A.A. Assessment of AquaCrop model for simulating growth and yield of saffron (Crocus sativus L.). Sci. Hortic. 2016, 211, 343–351. [Google Scholar] [CrossRef]
  27. Araya, A.; Keesstra, S.D.; Stroosnijder, L. Simulating yield response to water of Teff (Eragrostis tef) with FAO’s AquaCrop model. Field Crop. Res. 2010, 116, 196–204. [Google Scholar] [CrossRef]
  28. Razzaghi, F.; Babolhakami, A.; Sepaskhah, A.R. Is AquaCrop a useful tool for rapeseed growth and yield prediction in semi-arid regions: Model evaluation under different water-saving using long-term weather data. Theor. Appl. Climatol. 2023, 155, 1–21. [Google Scholar]
  29. Jiang, T.; Wang, B.; Xu, X.; Cao, Y.; Li Liu, D.; He, L.; Jin, N.; Ma, H.; Chen, S.; Zhao, K.; et al. Identifying sources of uncertainty in wheat production projections with consideration of crop climatic suitability under future climate. Agric. For. Meteorol. 2022, 319, 108933. [Google Scholar] [CrossRef]
  30. Solgi, S.; Ahmadi, S.H.; Sepaskhah, A.R.; Edalat, M. Wheat yield modeling under water-saving irrigation and climatic scenarios in transition from surface to sprinkler irrigation systems. J. Hydrol. 2022, 612, 128053. [Google Scholar] [CrossRef]
  31. Ul-Allah, S.; Rehman, A.; Hussain, M.; Farooq, M. Fiber yield and quality in cotton under drought: Effects and management. Agric. Water Manag. 2021, 255, 106994. [Google Scholar] [CrossRef]
  32. Wellens, J.; Raes, D.; Fereres, E.; Diels, J.; Coppye, C.; Adiele, J.G.; Ezui, K.S.G.; Becerra, L.; Selvaraj, M.G.; Dercon, G.; et al. Calibration and validation of the FAO AquaCrop water productivity model for cassava (Manihot esculenta Crantz). Agric. Water Manag. 2022, 263, 107491. [Google Scholar] [CrossRef]
  33. Shen, X.J.; Sun, J.S.; Zhang, J.Y.; Wang, J.L.; Li, M.S. Study on the placement of sensors for moisture content in soil profile for cotton under mulched drip irrigation condition. Agric. Res. Arid. Areas 2012, 30, 90–95. [Google Scholar]
  34. Hou, X.; Fan, J.; Zhang, F.; Hu, W.; Yan, F.; Xiao, C.; Li, Y.; Cheng, H. Determining water use and crop coefficients of drip-irrigated cotton in south Xinjiang of China under various irrigation amounts. Ind. Crop. Prod. 2022, 176, 114376. [Google Scholar] [CrossRef]
  35. Çetin, O.; Kara, A. Assesment of wa35ter productivity using different drip irrigation systems for cotton. Agric. Water Manag. 2019, 223, 105693. [Google Scholar] [CrossRef]
  36. Raes, D.; Steduto, P.; Hsiao, T.C.; Fereres, E. AquaCrop—The FAO crop model to simulate yield response to waicter: II. Main algorithms and software description. Agron. J. 2009, 101, 438–447. [Google Scholar] [CrossRef]
  37. Chow, V.T.; Maidment, D.R.; Mays, L.W. Applied hydrology. J. Eng. Educ. 1962, 308, 1959. [Google Scholar]
  38. Tsakiris, G.; Vangelis, H. Establishing a drought index incorporating evapotranspiration. Eur. Water 2005, 9, 3–11. [Google Scholar]
  39. Pierce, K.M.; Hope, J.L.; Hoggard, J.C.; Synovec, R.E. A principal component analysis based method to discover chemical differences in comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometry (GC×GC-TOFMS) separations of metabolites in plant samples. Talanta 2006, 70, 797–804. [Google Scholar] [CrossRef]
  40. Yavuz, D.; Seymen, M.; Süheri, S.; Yavuz, N.; Türkmen, Ö.; Kurtar, E.S. How do rootstocks of citron watermelon (Citrullus lanatus var. citroides) affect the yield and quality of watermelon under deficit irrigation? Agric. Water Manag. 2020, 241, 106351. [Google Scholar] [CrossRef]
  41. Yavuz, N. Can grafting affect yield and water use efficiency of melon under different irrigation depths in a semi-arid zone? Arab. J. Geosci. 2021, 14, 1–14. [Google Scholar] [CrossRef]
  42. Ran, H.; Kang, S.; Li, F.; Du, T.; Tong, L.; Li, S.; Ding, R.; Zhang, X. Parameterization of the AquaCrop model for full and deficit irrigated maize for seed production in arid Northwest China. Agric. Water Manage. 2018, 203, 438–450. [Google Scholar] [CrossRef]
  43. Tsakmakis, I.D.; Kokkos, N.P.; Gikas, G.D.; Pisinaras, V.; Hatzigiannakis, E.; Arampatzis, G.; Sylaios, G.K. Evaluation of AquaCrop model simulations of cotton growth under deficit irrigation with an emphasis on root growth and water extraction patterns. Agric. Water Manage. 2019, 213, 419–432. [Google Scholar] [CrossRef]
  44. Paredes, P.; de Melo-Abreu, J.P.; Alves, I.; Pereira, L.S. Assessing the performance of the FAO AquaCrop model to estimate maize yields and water use under full and deficit irrigation with focus on model parameterization. Agric. Water Manag. 2014, 144, 81–97. [Google Scholar] [CrossRef]
  45. Feng, D.; Li, G.; Wang, D.; Wulazibieke, M.; Cai, M.; Kang, J.; Yuan, Z.; Xu, H. Evaluation of AquaCrop model performance under mulched drip irrigation for maize in Northeast China. Agric. Water Manag. 2022, 261, 107372. [Google Scholar] [CrossRef]
  46. Sandhu, R.; Irmak, S. Assessment of AquaCrop model in simulating maize canopy cover, soil-water, evapotranspiration, yield, and water productivity for different planting dates and densities under irrigated and rainfed conditions. Agric. Water Manag. 2019, 224, 105753. [Google Scholar] [CrossRef]
  47. Wang, X.; Jiang, F.; Wang, H.; Cao, H.; Yang, Y.; Gao, Y. Irrigation Scheduling Optimization of Drip-irrigated without Plastic Film Cotton in South Xinjiang Based on AquaCrop model. Trans. Chin. Soc. Agric. Mach. 2021, 52, 293–301. [Google Scholar]
  48. Du, T.; Kang, S.; Zhang, J.; Li, F. Water use and yield responses of cotton to alternate partial root-zone drip irrigation in the arid area of north-west China. Irrig. Sci. 2008, 26, 147–159. [Google Scholar] [CrossRef]
  49. Xiao, C.; Ji, Q.; Zhang, F.; Li, Y.; Fan, J.; Hou, X.; Yan, F.; Liu, X.; Gong, K. Effects of various soil water potential thresholds for drip irrigation on soil salinity, seed cotton yield and water productivity of cotton in northwest China. Agric. Water Manag. 2023, 279, 108172. [Google Scholar] [CrossRef]
  50. Jiang, T.; Sun, S.; Li, Z.; Li, Q.; Lu, Y.; Li, C.; Wang, Y.; Wu, P. Vulnerability of crop water footprint in rain-fed and irrigation agricultural production system under future climate scenarios. Agric. For. Meteorol. 2022, 326, 109164. [Google Scholar] [CrossRef]
  51. Li, N.; Li, Y.; Yang, Q.; Biswas, A.; Dong, H. Simulating climate change impacts on cotton using AquaCrop model in China. Agr. Syst. 2024, 216, 103897. [Google Scholar] [CrossRef]
  52. He, P.; Yu, S.E.; Zhang, F.; Ma, T.; Ding, J.; Chen, K.; Chen, X.; Dai, Y. Effects of Soil Water Regulation on the Cotton Yield, Fiber Quality and Soil Salt Accumulation under Mulched Drip Irrigation in Southern Xinjiang, China. Agronomy 2022, 12, 1246. [Google Scholar] [CrossRef]
  53. An, J.; Zhang, Z.; Li, X.; Xing, F.; Lei, Y.; Yang, B.; Wang, Z.; Han, Y.; Chen, H.; Wang, G.; et al. Loose and tower-type canopy structure can improve cotton yield in the Yellow River basin of China by increasing light interception. Arch. Agron. Soil Sci. 2023, 69, 920–933(14). [Google Scholar] [CrossRef]
  54. Wang, Q.; Chen, H.; Han, Y.; Xing, F.; Wang, Z.; Feng, L.; Wang, G.; Yang, B.; Lei, Y.; Xiong, S.; et al. Effect of Spatial-Temporal Light Competition on Cotton Yield and Yield Distribution. Agronomy 2021, 11, 2346. [Google Scholar] [CrossRef]
  55. De Mello, C.S.; Van Dijk, J.P.; Voorhuijzen, M.; Kok, E.J.; Arisi, A.C.M. Tuber proteome comparison of five potato varieties by principal component analysis. J. Sci. Food Agric. 2016, 96, 3928–3936. [Google Scholar] [CrossRef] [PubMed]
  56. Jarwar, A.H.; Wang, X.; Iqbal, M.S.; Sarfraz, Z.; Wang, L.; Ma, Q.; Shuli, F. Genetic divergence on the basis of principal component, correlation and cluster analysis of yield and quality traits in cotton cultivars. Pak. J. Bot. 2019, 51, 1143–1148. [Google Scholar] [CrossRef]
  57. Yu, H.; Li, S.; Ding, J.; Yang, T.; Wang, Y. Water use efficiency and its drivers of two typical cash crops in an arid area of Northwest China. Agric. Water Manag. 2023, 287, 108433. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Agriculture 14 01269 g001
Figure 2. Calibration (2021) of the simulated and observed total soil water (TSW) in the 0–60 cm soil profile. The W1-W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50–100% FC, W2: 60–80% FC, W3: 60–90% FC, W4: 60–100% FC.
Figure 2. Calibration (2021) of the simulated and observed total soil water (TSW) in the 0–60 cm soil profile. The W1-W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50–100% FC, W2: 60–80% FC, W3: 60–90% FC, W4: 60–100% FC.
Agriculture 14 01269 g002
Figure 3. Validation (2022) of the simulated and observed values of total soil water (TSW) in the 0–60 cm soil profile. The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Figure 3. Validation (2022) of the simulated and observed values of total soil water (TSW) in the 0–60 cm soil profile. The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Agriculture 14 01269 g003
Figure 4. Simulation and measurement of the canopy cover (CC) curves for cotton under all treatments during the calibration period (2021). The W1–W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50%–100% FC, W2: 60%–80% FC, W3: 60%–90% FC, W4: 60%–100% FC.
Figure 4. Simulation and measurement of the canopy cover (CC) curves for cotton under all treatments during the calibration period (2021). The W1–W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50%–100% FC, W2: 60%–80% FC, W3: 60%–90% FC, W4: 60%–100% FC.
Agriculture 14 01269 g004
Figure 5. Simulation and measurement of canopy cover (CC) curves for cotton under all treatments during the validation period (2022). The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Figure 5. Simulation and measurement of canopy cover (CC) curves for cotton under all treatments during the validation period (2022). The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Agriculture 14 01269 g005
Figure 6. The development of the aboveground biomass of cotton was simulated and measured in all treatments during the calibration period (2021). The W1–W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50–100% FC, W2: 60–80% FC, W3: 60–90% FC, W4: 60–100% FC.
Figure 6. The development of the aboveground biomass of cotton was simulated and measured in all treatments during the calibration period (2021). The W1–W4 treatments had the same upper and lower limits of irrigation at the SQS and BLS. W1: 50–100% FC, W2: 60–80% FC, W3: 60–90% FC, W4: 60–100% FC.
Agriculture 14 01269 g006
Figure 7. The development of the aboveground biomass of cotton was simulated and measured in all treatments (aj) during the validation period (2022). The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Figure 7. The development of the aboveground biomass of cotton was simulated and measured in all treatments (aj) during the validation period (2022). The dots represent the measured data, and the line represents the simulated values. Three lower irrigation limits of 55%, 65%, and 75% FC were set at the SQS, and three lower irrigation limits of 60%, 70%, and 80% FC were set at the BLS. Nine treatments (WD1–WD9, i.e., (ai)) were completely combined, and CK (j) was the control treatment.
Agriculture 14 01269 g007
Figure 8. Comparison between the measured and simulated values of aboveground biomass (a) and seed cotton yield (b) in 2021 (calibration) and 2022 (validation). The dots represent the relationship between the observations and simulations. In the figure, * means p <= 0.05, *** means p <= 0.001.
Figure 8. Comparison between the measured and simulated values of aboveground biomass (a) and seed cotton yield (b) in 2021 (calibration) and 2022 (validation). The dots represent the relationship between the observations and simulations. In the figure, * means p <= 0.05, *** means p <= 0.001.
Agriculture 14 01269 g008
Figure 9. Comparison of the total costs and net revenue of cotton in wet years, normal years, and dry years. The blue bar chart shows the TC in wet years, the green bar chart shows the TC in normal years, and the red bar chart shows the TC in dry years. The colored bar chart in each subgraph shows the NR of each simulation process. Bars are the means + one standard error of the mean (n = 5). Different letters above the error bars indicate a significant difference at p < 0.05 according to the Duncan test.
Figure 9. Comparison of the total costs and net revenue of cotton in wet years, normal years, and dry years. The blue bar chart shows the TC in wet years, the green bar chart shows the TC in normal years, and the red bar chart shows the TC in dry years. The colored bar chart in each subgraph shows the NR of each simulation process. Bars are the means + one standard error of the mean (n = 5). Different letters above the error bars indicate a significant difference at p < 0.05 according to the Duncan test.
Agriculture 14 01269 g009
Figure 10. Correlation analysis of Ya, Bio, ET, TI, IF, WPIrrig, WP, HI, TC, and NR in different scenarios. (Ya: seed cotton yield, Bio: biomass, ET: evapotranspiration, TI: total irrigation, IF: irrigation frequency, WPIrrig: water productivity with respect to the amount of irrigation water applied, WP: water productivity with respect to crop consumptive water use, HI: harvest index, TC: total cost, NR: net revenue).
Figure 10. Correlation analysis of Ya, Bio, ET, TI, IF, WPIrrig, WP, HI, TC, and NR in different scenarios. (Ya: seed cotton yield, Bio: biomass, ET: evapotranspiration, TI: total irrigation, IF: irrigation frequency, WPIrrig: water productivity with respect to the amount of irrigation water applied, WP: water productivity with respect to crop consumptive water use, HI: harvest index, TC: total cost, NR: net revenue).
Agriculture 14 01269 g010
Figure 11. Fractional loading biplots created using PC1 and PC2 in three scenario years (Ya: seed cotton yield, Bio: biomass, ET: evapotranspiration, TI: total irrigation, IF: irrigation frequency, WPIrrig: water productivity with respect to the amount of irrigation water applied, WP: water productivity with respect to crop consumptive water use, HI: harvest index, TC: total cost, NR: net revenue).
Figure 11. Fractional loading biplots created using PC1 and PC2 in three scenario years (Ya: seed cotton yield, Bio: biomass, ET: evapotranspiration, TI: total irrigation, IF: irrigation frequency, WPIrrig: water productivity with respect to the amount of irrigation water applied, WP: water productivity with respect to crop consumptive water use, HI: harvest index, TC: total cost, NR: net revenue).
Agriculture 14 01269 g011
Figure 12. Schematic diagram of the comprehensive evaluation of irrigation options for machine-picked cotton in the dry oasis zone in different scenario years.
Figure 12. Schematic diagram of the comprehensive evaluation of irrigation options for machine-picked cotton in the dry oasis zone in different scenario years.
Agriculture 14 01269 g012
Table 1. Irrigation test scheme and actual irrigation situation.
Table 1. Irrigation test scheme and actual irrigation situation.
YearTreatmentIrrigation ThresholdIrrigation Quotas and Frequency of IrrigationIrrigation during Critical Fertility Period (mm)Whole Fertility Cycle
Squaring StageBlooming StageSeeding Stage
(mm)
Squaring Stage (mm)Blooming Stage (mm)Total Irrigation
(mm)
2021W150–100%50–100%7042 × 266 × 4346416
W260–80%60–80%7017 × 426 × 11356426
W360–90%60–90%7025 × 339 × 8390460
W460–100%60–100%7034 × 352 × 6415485
2022WD155–90%60–90%7029 × 239 × 6295365
WD270–100%70–100%7029 × 226 × 10321391
WD355–90%60–90%7029 × 213 × 21334404
WD455–90%70–90%7021 × 339 × 6299369
WD555–90%80–90%7021 × 326 × 10325395
WD665–90%80–90%7021 × 313 × 21338408
WD765–90%70–90%7013 × 539 × 7338408
WD865–90%80–90%7013 × 526 × 11351421
WD975–90%60–90%7013 × 713 × 20350420
CK75–90%70–90%70--575645
Note: The irrigation schemes for 2021 and 2022 were not the same; there were upper and lower irrigation restrictions during the square period and flowering period. The treatments selected in 2021 were only for model calibration. The irrigation amount and irrigation cycle refer to the squaring stage and blooming stage; the squaring stage and blooming stage are the key periods of fertility, and the total irrigation is the irrigation amount for the whole fertility cycle of cotton.
Table 2. Crop parameters of the AquaCrop model.
Table 2. Crop parameters of the AquaCrop model.
Crop ParameterValueRemarks
Base temperature (Tbase)/(°C)12Measured
Upper temperature (Tupper)/(°C)35Measured
Crop transpiration coefficient KcTR,x1.30Calibrated
Initial canopy cover (CC0)/(%)1.22Calibrated
Canopy growth coefficient (CGC)/(% d−1)6.9Calibrated
Maximum canopy cover (CCX)/(%)88Calibrated
Canopy decline coefficient (CDC)/(% d−1)8.8Calibrated
Normalized crop water productivity (WP)/(g m−2)24Calibrated
Reduction coefficient of WP* during yield formation (f yield)74Calibrated
Maximum effective rooting depth/(m)0.70Recommended
Reference harvest index (HI0)/(%)35Calibrated
Upper threshold for canopy expansion (Pexp,upper)0.2Calibrated
Lower threshold for canopy expansion (Pexp,lower)0.64Calibrated
Upper threshold for stomatal closure (Pclo,upper)0.32Calibrated
Upper threshold for canopy senescence (Psen,upper)0.6Calibrated
Lower threshold of the impact of salt on crop growth (Ece,lower)/(dS m−1)8Calibrated
Upper limit of impact threshold of salt on crop growth (Ece,upper)/(dS m−1)27Calibrated
Table 3. Implementation of comprehensive historical climate scenarios in the AquaCrop model.
Table 3. Implementation of comprehensive historical climate scenarios in the AquaCrop model.
Simulation ScenarioIrrigation Amount (mm)Irrigation Cycle (d)Irrigation FrequencySimulation ScenarioIrrigation Amount (mm)Irrigation Cycle (d)Irrigation Frequency
T1360616T9450616
T2360812T10450812
T33601010T114501010
T4360128T12450128
T5405616T13525616
T6405812T14525812
T74051010T155251010
T8405128T16525128
Note: The irrigation amount in the table refers to cotton irrigation during the critical growth period (mm); the irrigation cycle refers to the critical growth irrigation cycle (d). Irrigation frequency refers to the number of times the water was filled during the critical growth period. The key growth periods were the SQS and the BLS. In addition, the seedling stage was irrigated four times according to the irrigation volume of 70 mm, and the boll opening stage (BOS) was not irrigated according to the local irrigation habits.
Table 4. Simulated and measured crop evapotranspiration (ETc).
Table 4. Simulated and measured crop evapotranspiration (ETc).
YearTreatmentETc (mm)
SimulatedMeasuredPe (%)Dif
2021W1441445−0.70−3
W24904811.728
W34924802.4312
W4499512−2.70−14
2022WD1421436−3.53−15
WD2427448−4.80−22
WD3429456−6.11−28
WD4419423−1.17−5
WD5433447−2.97−13
WD6432462−6.47−30
WD7432468−7.67−36
WD8434474−8.50−40
WD9436485−10.06−49
Note: Dif: difference between the measured and simulated ETc; Pe: prediction error of the simulated ETc with respect to the measured values.
Table 5. Comparison of the final yield, biomass, Bio, ET, WP, and WPIrrig of cotton in different comprehensive scenarios in wet years, normal years, and dry years.
Table 5. Comparison of the final yield, biomass, Bio, ET, WP, and WPIrrig of cotton in different comprehensive scenarios in wet years, normal years, and dry years.
TreatmentScenario Simulation
Wet YearsNormal YearsDry Years
YaBioETWPWPIrrigHIYaBioETWPWPIrrigHIYaBioETWPWPIrrigHI
T16.25 e16.10 f509.52 fg1.23 cd1.46 a38.86 b5.16 fg12.93 h506.78 i1.02 ef1.20 a39.96 cde4.66 g11.44 h507.30 h0.92 de1.08 a40.78 def
T26.18 f15.70 g496.24 hi1.25 b1.44 b39.42 a4.99 hi12.37 i488.16 j1.02 de1.16 bc40.36 bcd4.51 h10.94 i486.80 j0.93 cd1.05 b41.28 cd
T36.20 ef15.65 g491.40 i1.26 a1.44 ab39.62 a5.09 gh12.52 i488.32 j1.04 ab1.18 ab40.72 ab4.62 g11.02 i486.74 j0.95 b1.07 a41.90 ab
T46.05 g15.24 h478.56 j1.27 a1.41 c39.74 a4.91 i11.96 j467.64 k1.05 a1.14 cd41.12 a4.46 h10.54 j463.52 k0.96 a1.04 bc42.38 a
T56.37 cd16.67 e525.48 de1.21 ef1.34 d38.20 d5.30 e13.47 f527.90 fg1.01 fg1.12 de39.38 ef4.85 e12.02 f532.24 f0.91 ef1.02 d40.34 fg
T66.35 cd16.55 e521.00 de1.22 de1.34 d38.42 cd5.23 ef13.15 g515.06 hi1.02 ef1.10 e39.86 de4.84 e11.73 g517.54 g0.94 c1.02 d41.25 cd
T76.32 d16.30 f510.36 f1.24 bc1.33 d38.76 bc5.28 ef13.11 g508.22 hi1.04 abcd1.11 e40.32 bcd4.87 e11.63 g510.24 h0.96 ab1.03 cd41.90 ab
T86.23 ef16.07 f502.26 gh1.24 b1.31 e38.82 b5.15 fg12.77 h495.44 j1.04 abc1.09 e40.48 bc4.76 f11.30 h494.32 i0.96 a1.00 e42.16 ab
T96.47 ab17.33 c543.80 bc1.19 g1.24 f37.34 g5.47 bc14.14 d549.96 cd1.00 gh1.05 f38.78 gh5.12 c12.74 d563.36 c0.91 ef0.98 f40.18 g
T106.46 ab17.23 c539.28 bc1.20 fg1.24 f37.46 fg5.45 cd13.92 e542.70 de1.01 fg1.05 f39.18 fg5.11 c12.47 e546.82 d0.94 c0.98 f40.98 de
T116.41 bc16.95 d528.88 d1.21 ef1.23 fg37.80 ef5.48 bc13.85 e535.00 ef1.03 cde1.05 f39.62 ef5.15 c12.36 e539.44 e0.96 ab0.99 ef41.72 bc
T126.35 cd16.65 e518.46 e1.23 cde1.22 g38.14 de5.33 de13.40 f518.30 gh1.03 bcde1.03 f39.84 de4.97 d11.94 f519.66 g0.96 ab0.96 g41.68 bc
T136.52 a17.95 a561.70 a1.16 h1.10 h36.32 i5.71 a15.17 a589.40 a0.97 i0.96 g37.70 i5.47 a13.83 a606.04 a0.90 f0.92 h39.58 h
T146.51 a17.82 a556.76 a1.17 h1.09 h36.54 hi5.68 a14.92 b577.68 b0.99 hi0.96 g38.36 h5.42 a13.52 b588.96 b0.92 de0.91 h40.18 g
T156.48 a17.58 b546.86 b1.19 g1.09 h36.88 h5.68 a14.79 b570.08 b1.00 gh0.95 g38.50 h5.47 a13.45 b583.94 b0.94 c0.92 h40.70 efg
T166.45 ab17.32 c538.38 c1.20 fg1.09 h37.24 g5.59 ab14.49 c556.92 c1.01 fg0.94 g38.66 gh5.30 b13.10 c565.74 c0.94 c0.89 i40.54 efg
Note: Ya is the seed cotton yield (t ha−1), Bio is the final aboveground biomass of cotton (t ha−1), ET is the evapotranspiration of cotton throughout its growing period (mm), WP is the water productivity with respect to the crop consumptive water use (kg m−3), WPIrrig is the water productivity with respect to the amount of irrigation water applied (kg m−3), and HI is the cotton harvest index (%). Letters after numbers indicate a significant difference at p < 0.05, according to the Duncan test.
Table 6. Factor load and variance contribution rate according to principal component analysis.
Table 6. Factor load and variance contribution rate according to principal component analysis.
EigenvectorScenario
Wet YearsNormal YearsDry Years
PC1PC2PC1PC2PC1PC2
Ya0.98 0.07 0.99 −0.08 0.97 0.19
Bio0.99 −0.02 1.00 0.00 1.00 0.03
ET0.99 0.12 1.00 0.08 1.00 −0.06
TI0.94 −0.31 0.96 −0.27 0.95 0.29
IF0.29 0.95 0.27 0.96 0.30 −0.94
WPIrrig−0.92 0.37 −0.90 0.38 −0.88 −0.40
WP−0.98 −0.17 −0.90 −0.40 −0.56 0.82
HI−1.00 0.01 −0.98 −0.18 −0.82 0.55
TC0.99 −0.11 0.99 −0.14 0.97 0.22
NR0.96 0.15 0.99 −0.06 0.97 0.19
Eigenvalue8.58 1.21 8.48 1.37 7.58 2.22
Percentage of variance85.83 12.12 84.78 13.65 75.76 22.22
Cumulative variance85.83 97.95 84.78 98.44 75.76 97.98
Note: Ya: seed cotton yield, Bio: biomass, ET: evapotranspiration, TI: total irrigation, IF: irrigation frequency, WPIrrig: water productivity with respect to the amount of irrigation water applied, WP: water productivity with respect to crop consumptive water use, HI: harvest index, TC: total cost, NR: net revenue.
Table 7. The comprehensive score of each treatment based on Ya, Bio, ET, WPIrrig, WP, HI, TI, IF, TC, and NR.
Table 7. The comprehensive score of each treatment based on Ya, Bio, ET, WPIrrig, WP, HI, TI, IF, TC, and NR.
TreatmentScenariosRank
Wet YearNormal YearDry YearWet YearNormal YearDry Year
ScorePC1PC2ScorePC1PC2ScorePC1PC2
T1−0.44−0.751.76−0.39−0.741.74−0.94−0.71−1.73111114
T2−1.00−1.210.51−0.89−1.140.69−1.03−1.08−0.88141415
T3−1.14−1.28−0.16−1.04−1.19−0.16−0.91−1.200.04151513
T4−1.65−1.76−0.86−1.48−1.60−0.73−1.06−1.570.67161616
T50.16−0.031.520.14−0.081.50−0.38−0.03−1.588711
T6−0.14−0.220.44−0.28−0.380.35−0.36−0.38−0.3191010
T7−0.46−0.48−0.34−0.55−0.55−0.55−0.32−0.590.6012129
T8−0.79−0.77−0.92−0.83−0.80−0.98−0.41−0.841.04131312
T90.800.721.330.690.601.240.230.67−1.25447
T100.520.560.260.360.390.160.250.320.03666
T110.180.28−0.530.090.22−0.720.300.110.96785
T12−0.17−0.04−1.11−0.26−0.11−1.180.14−0.161.181098
T131.471.560.841.611.701.061.141.74−0.91112
T141.191.39−0.231.201.41−0.141.121.410.14223
T150.891.15−0.990.961.26−0.901.191.260.94331
T160.590.89−1.530.681.01−1.371.041.031.07554
Linear functionw = 0.876 w1 +0.124 w2n = 0.862 n1 + 0.139 n2d = 0.773 d1 + 0.227 d2T13T13T15
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Du, Y.; Fu, Q.; Ai, P.; Ma, Y.; Pan, Y. Modeling Comprehensive Deficit Irrigation Strategies for Drip-Irrigated Cotton Using AquaCrop. Agriculture 2024, 14, 1269. https://doi.org/10.3390/agriculture14081269

AMA Style

Du Y, Fu Q, Ai P, Ma Y, Pan Y. Modeling Comprehensive Deficit Irrigation Strategies for Drip-Irrigated Cotton Using AquaCrop. Agriculture. 2024; 14(8):1269. https://doi.org/10.3390/agriculture14081269

Chicago/Turabian Style

Du, Yalong, Qiuping Fu, Pengrui Ai, Yingjie Ma, and Yang Pan. 2024. "Modeling Comprehensive Deficit Irrigation Strategies for Drip-Irrigated Cotton Using AquaCrop" Agriculture 14, no. 8: 1269. https://doi.org/10.3390/agriculture14081269

APA Style

Du, Y., Fu, Q., Ai, P., Ma, Y., & Pan, Y. (2024). Modeling Comprehensive Deficit Irrigation Strategies for Drip-Irrigated Cotton Using AquaCrop. Agriculture, 14(8), 1269. https://doi.org/10.3390/agriculture14081269

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

Article Metrics

Back to TopTop