Next Article in Journal
Identification of Donors for Fresh Seed Dormancy and Marker Validation in a Diverse Groundnut Mini-Core Collection
Previous Article in Journal
Analysis and Research on Production Effects of Full-Film Double Ridge-Furrow Mulching with Polyethylene Film and Biodegradable Film Based on AquaCrop
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulation of Winter Wheat Growth Dynamics and Optimization of Water and Nitrogen Application Systems Based on the Aquacrop Model

1
School of Water Conservancy, North China University of Water Resources and Electric Power, Zhengzhou 450046, China
2
Collaborative Innovation Center for the Efficient Utilization of Water Resources, Zhengzhou 450046, China
*
Author to whom correspondence should be addressed.
Agronomy 2024, 14(1), 110; https://doi.org/10.3390/agronomy14010110
Submission received: 10 December 2023 / Revised: 27 December 2023 / Accepted: 28 December 2023 / Published: 1 January 2024
(This article belongs to the Section Innovative Cropping Systems)

Abstract

:
Winter wheat is the main grain crop in the Yellow River Basin, and optimizing water and nitrogen management can not only improve the yield, but also reduce water and fertilizer waste and environmental pollution. Using two years of winter wheat field trials, the Aquacrop model was calibrated and validated, and 36 different water and nitrogen application scenarios were simulated and analyzed. Among them, the field trials were divided into three N application levels, and three lower limits of water control. The results showed that: (1) the model has good applicability to simulate the growth process of winter wheat in the Yellow River Basin. (2) At a constant level of irrigation or nitrogen application, crop yield initially increased and subsequently decreased with the enhancement of either irrigation or nitrogen application. Finally, through response surface analysis of regression equations, the recommended irrigation volume was determined to be between 296 and 303 mm, and the fertilizer application rate was set between 229 and 261 kg/ha. At these levels, the PFP could achieve over 60% of its maximum. This research provides a theoretical basis for the water-saving, fertilizer-saving and high-yield water and fertilizer management system of winter wheat in the irrigated area of the Yellow River Basin.

1. Introduction

The irrigation area of the Yellow River Basin holds a pivotal role as a grain production base in China, with winter wheat being the primary food crop cultivated in this region. Ensuring high crop yields within this irrigation district is critically significant for safeguarding China’s food security [1]. Meng Zhaogang et al. systematically studied the coupling effect of water and fertilizer between the upper and lower crops in the one-year two-crop system of winter wheat and summer maize in the irrigation area of the Yellow River Basin, and put forward the water and fertilizer optimization scheme under different decision-making objectives [2]. Currently, a notable challenge facing China’s agricultural production is the scarcity of water resources, coupled with low water and fertilizer utilization efficiency [3]. Two principal factors contribute to this inefficiency: firstly, the prevalent irrigation methodology, with flood irrigation being the predominant choice in China’s irrigation zones. This method leads to a considerable waste of water resources and loss of soil fertility. Secondly, the excessive application of fertilizers causes the fertilizer waste and an increase in pest infestations and diseases, which ultimately hinder crop yield. Therefore, achieving both high and stable yields while ensuring rational utilization of water and fertilizers has continually been a direction of research for a number of researchers [4]. Within the irrigated area of the Yellow River Basin, the total availability of water resources is markedly limited. The water allocation from the Yellow River often falls short of meeting the irrigation requirements, especially during dry seasons. Concurrently, a gradual decline in soil fertility adversely affects both the enhancement of agricultural yield and the sustainable advancement of agriculture. Therefore, optimizing the water and nitrogen management regimes emerges as a paramount agenda for enhancing agricultural development in the irrigation area of the Yellow River Basin [5].
The exploration of optimal irrigation and fertilization strategies through field trials can be time-consuming and susceptible to external influences [6] hence, employing crop models to simulate crop growth has evolved as a vital tool in modern agriculture. Currently, prevalent crop models include DSSAT, WOFOST, and AquaCrop, among others. Notably, by focusing on water use efficiency (WUE), AquaCrop has the advantages of requiring fewer input parameters, and having a simple and intuitive operation interface [7,8]. Numerous researchers have utilized the AquaCrop model to simulate crop growth across various global regions, thereby evaluating its applicability. For instance, Farahani et al. [9] calibrated parameters suitable for simulating cotton yield in the Mediterranean region, establishing a foundation for future applications of this model in the area. Similarly, Araya et al. [10] demonstrated the model’s efficacy in simulating the aboveground biomass and yield of barley under varying planting dates. Toumi et al. [11] employed the model to simulate the canopy cover and soil moisture content of winter wheat, validating the AquaCrop model’s feasibility in capturing the growth dynamics of winter wheat throughout its lifecycle. Andarzian et al. [12] demonstrated, based on simulated data for winter wheat in southern Iran, that under normal precipitation conditions, the AquaCrop model can accurately predict the pattern of soil moisture change.
Furthermore, many studies have been conducted to explore the effects of distinct climates on crop yields using AquaCrop model. At present, in China, several scholars have applied the Aquacrop model to simulate crop growth in different regions and evaluated the applicability of the model: He et al. [13] examined the impacts on tobacco yield in Yunnan province, showcasing high simulation accuracy. Ni et al. [14] evaluated the suitability of summer maize production in the Loess Plateau area, which is greatly beneficial for simulating water optimization management of summer maize in this region. Pan et al. [15] applied the model to simulate maize yield in northwest China employing full-film double-furrow sowing technology and concluded the model’s suitability for simulating local dry farming. Additional studies by Meng et al. [16] and Du et al. [17] successfully simulated the growth and yield of winter wheat in Hebei province of China, achieving high accuracy in simulation results. Zhu et al. [18] explored alternating fresh and saline irrigation methods for winter wheat during different fertility periods, calibrating crop growth and yield accordingly. Cui et al. [19] utilized the one factor at a time (OTA) method to conduct a sensitivity analysis of model parameters, then optimized the sensitive parameters, concluding that the calibrated model could accurately simulate corn and soybean yields in northeast China.
Moreover, the AquaCrop model has been employed to simulate the growth and yield of various crops such as corn [20,21], soybean [22] peanut [23], and rapeseed [24] all achieving high accuracy. Some researchers have extended the model’s application to crop economic efficiency analysis and evaluation [25] and irrigation regime optimization [26,27,28,29], achieving favorable outcomes. Overall, the AquaCrop model, post-parameter calibration, exhibits substantial applicability to different crops and demonstrates higher simulation accuracy compared to other models [30].
In this study, we selected winter wheat in the irrigation area of the Yellow River Basin as the research subject to evaluate the applicability of the AquaCrop model. Experimental treatments were established with varying irrigation water and amounts of nitrogen applied, over a span of two years, to gather field test data and corresponding data analysis results. Subsequently, the AquaCrop model was calibrated and validated, employing the collected experimental data. The calibrated model was then utilized to perform scenario simulations of winter wheat irrigation during different fertility periods to determine the optimum irrigation period for enhancing winter wheat fertility. The experimental data were fitted to the water–fertilizer coupling regression equation [31], and through this equation, a response surface map illustrating the effects of water and fertilizer on yield, WUE, and fertilizer partial factor productivity (PFP) was generated [32]. This visual representation provides an intuitive and comprehensive understanding of the impacts of varying irrigation volumes and fertilizer application rate on crop growth, as well as the overall trend of water–fertilizer coupling. By doing so, we aimed to obtain the optimal irrigation volume and fertilizer application rate for winter wheat in the irrigation area of the Yellow River Basin, to ensure high yield while conserving water and fertilizer resources. Our investigation not only provides a reference for the applicability of the AquaCrop model in the irrigation area of the Yellow River Basin, but also highlights the importance of optimizing water and nitrogen use. Through the integration of scenario simulations, we can ensure a higher yield of winter wheat in the region, improve the effectiveness of water and nitrogen utilization, mitigate soil deterioration, and reduce wasteful water use.

2. Materials and Methods

2.1. Overview of the Experimental Area

The study was conducted at the Experimental Field for Efficient Water Use in Agriculture, situated within the Longzi Lake Campus of the North China University of Water Resources and Hydropower, Zhengzhou. This location lies in the heart of the Central China Plain with the geographic coordinates 34.78° N latitude and 113.76° E longitude, at an elevation of 110 m above sea level. The local climate is characterized by a warm temperate continental monsoon climate with an average annual temperature of approximately 14.5 °C. The region receives an average annual precipitation of 637.1 mm, enjoys an average sunshine duration of 6.57 h per day, and experiences a frost-free period extending over 220 days. The experimental site features flat terrain with clay loam being the predominant soil type.
The winter wheat cultivation cycle in this region spans from planting in October of one year to harvesting in June of the following year. Accordingly, the experimental periods for this study were specified as October 2021 to June 2022 and October 2022 to June 2023. The precise location of the experimental site is depicted in Figure 1.

2.2. Experimental Design

In this study, the amounts of irrigation and fertilizer application were determined by the climatic conditions of the Yellow River Irrigation District. In addition, a study on the optimal allocation of water and nitrogen for winter wheat in north China based on the quantification of crop water-fertilizer coupling types was conducted [33]. Three levels of nitrogen application were established: N1 (120 kg/ha), N2 (220 kg/ha), and N3 (320 kg/ha). Concurrently, three distinct lower water control thresholds were designated at each N application level: W1 (60%θf), W2 (70%θf), and W3 (80%θf), where θf represents the field water holding capacity. For each treatment, nitrogen application was divided into a base fertilizer and two top-dressing fertilizations (base fertilizer for the first basic application and two surface applications for the two supplementary applications). The base fertilizer application was consistent at 67.5 kg/ha for all treatments, utilizing a ternary compound fertilizer (comprising 60 kg/ha of N, 60 kg/ha of K2O, and 60 kg/ha of P2O5). The top-dressing fertilizer utilized was urea with a nitrogen content of 46.3%. The timings for top-dressing fertilizations were during the re-greening (returning-green stage) and tasseling (heading stage) periods of winter wheat, respectively [34].
Irrigation treatments were conducted according to the targeted soil moisture content within the wet layer during each fertility period, initiating irrigation whenever the soil moisture descended below the lower threshold of the targeted soil moisture content. The irrigation quota for winter wheat was set at a fixed 45 mm, and Jimai 22 was chosen as the seed variety. Conventional border irrigation was employed for the field experiments, with each plot measuring 10 m × 1.5 m. To ensure experimental accuracy, a buffer zone of 1.5 m width was established between adjacent plots. Table 1 outlines the specific water and fertilizer treatment configurations for the experimental plots:

2.3. Assessment Indicators and Methods

2.3.1. Aboveground Biomass and Yield

(1)
Aboveground biomass
The aboveground biomass was assessed every 10 days following the overwintering phase of winter wheat. In each treatment, five plants were selected for measurement. Plants were cut at ground level, and the aboveground parts were initially placed in a drying oven at 105 °C for 30 min to halt the greening process. The samples were then kept in the oven and dried at 75 °C until a constant mass was achieved.
(2)
Yield
The yield of winter wheat was calculated by multiplying the number of spikes per unit area, the number of grains per spike, and the mass of 1000 grains at harvest. The effective number of spikes within a 1 m2 plot was determined. In each treatment, five representative plants were chosen, and the number of grains per spike on each plant was counted. Additionally, a set of 1000 wheat seeds harvested randomly was weighed to obtain the thousand-grain mass value.

2.3.2. Canopy Cover

The canopy cover of winter wheat during its reproductive period can be derived from the leaf area index (LAI) values. These LAI values were acquired using the values of several measured parameters shown in Equation (1) [17]:
L A I = 0.75 × ρ × i = 1 m i = 1 n L i j × B i j m
where m is the number of winter wheat plants per unit area; ρ is the planting density; Lij and Bij are the length and the width (at the widest part) of the leaf blade (in mm), respectively; n is the number of leaves on the selected plant; 0.75 is the correction coefficient.
Canopy cover was then calculated from the derived LAI value, as shown in Equation (2):
C C = 1 exp 0.65 × L A I

2.3.3. Soil Water Content

The soil water content was monitored throughout the reproductive cycle of winter wheat using the time domain reflectometry with intelligent microelements (TRIME) method [35]. The TRIME method combines the proven principles of time domain reflectometry (TDR) with modern electronics and software algorithms to provide a robust tool for measuring soil water content. Measurements were conducted every 10 days following the overwintering phase. Prior to the commencement of the experiment, the TRIME probes were calibrated using the drying method. Measurements were obtained at a depth of 100 cm and at intervals of 20 cm.

2.3.4. WUE and Fertilizer PFP

The WUE is a critical metric that reflects the effectiveness of water utilization in crop production. It is calculated using the following formula [33]:
W U E = Y / E T
where Y stands for the yield (kg/ha); ET is the crop water consumption (mm), which can be calculated in the following equation:
E T = P + U + I R D Δ W
where P is the rainfall during the reproductive period (mm); U is the groundwater recharge (mm); I is the irrigation volume during the reproductive period (mm); D is the deep seepage (mm); R is the surface runoff (mm); and Δ W represents the change in soil moisture from the beginning to the end of the experiment (mm). According to the actual measurements, there was no significant change in the soil moisture content at a depth of 1 m in the experimental area. Additionally, given the deep groundwater level and relatively flat terrain in the experimental area, the terms U, D, and R in the formula can be neglected.
The expression for calculating fertilizer PFP is represented as follows:
P F P = Y / T
where Y is the yield (kg/ha); and T is the total amount of fertilizer applied during the reproductive period (kg/ha).

3. Results

3.1. Model Fundamentals

AquaCrop is a water-driven, daily timestep, crop model issued by the Food and Agriculture Organization of the United Nations (FAO) in 2009, with the aim of achieving an optimal balance between simplicity, accuracy, and robustness. To foster its utilization, AquaCrop is configured to require few parameters and intuitive input data that are either readily accessible or obtainable through straightforward measurements.
The operation principle of the AquaCrop model was elucidated comprehensively by Raes et al. [35] and Steduto et al. [36], encompassing the following principal simulation aspects: (1) simulation of crop canopy development; (2) simulation of crop transpiration; (3) simulation of crop aboveground biomass growth; (4) simulation of crop yields. The main equations are given as:
Y = B × H I
B = W P * × T r
where Y is the winter wheat yield (kg/ha); B is the aboveground biomass (kg/ha); HI is the winter wheat harvest index (%); WP* is the standardized water productivity (g/m2); and Tr is the crop transpiration (mm).
AquaCrop model adopts a semi-quantitative approach to illustrate the effects of fertilizer stress on crop growth without explicitly addressing nutrient cycling and balancing [30]. The model incorporates a scaling factor, Brel, relative to the aboveground biomass that can be obtained under conditions without fertilizer and water stress, as expressed in Equation (8).
B r e l = B s t r B r e f × 100 %
where Bstr is the aboveground biomass (t/ha) in the absence of both fertilizer and water stress; Bref is the aboveground biomass (t/ha) under fertilizer stress but devoid of water stress.
To investigate the impact of fertilizer stress on crop growth, four fertilizer stress coefficients were incorporated in the model: canopy expansion stress coefficient, maximum canopy cover, canopy decay stress coefficient, and biomass water productivity. This corresponded to the reduced maximum canopy cover, the decrease in canopy growth coefficient, the average canopy reduction, and the standardized crop water productivity reduction within the parameter calibration results as seen in Table 2. The calibrated fertilizer stress coefficients are tabulated in Table 2.

3.2. Creation of the AquaCrop Modeling Database

3.2.1. Meteorological Data

The meteorological data, encompassing minimum/maximum air temperatures (°C), daily precipitation (mm), and wind velocity (m/s) within the experimental area, were acquired from a high-precision automatic meteorological station (HM-HL08). Evapotranspiration values, derived from the ET0 calculator integrated within the AquaCrop version 6.1 software, were used as a reference. This dataset facilitated the constitution of a meteorological database within the software. Daily weather data spanning two years are shown in Figure 2 and Figure 3.

3.2.2. Crop Data

The crop parameters were obtained during the growth and development of winter wheat within the experimental area. They were then input into the crop module of the AquaCrop model to formulate a crop database. This database includes information such as initial canopy cover, planting density, canopy cover across the entire lifecycle, and key crop developmental phases including attainment of maximum canopy cover, flowering, and the timing of maturity and senescence. The AquaCrop model offers two modes to depict crop growth and development: calendar mode and growing degree-day (GDD) mode. To accommodate varying climatic conditions across different years and to provide a more precise depiction of crop growth, the GDD mode was employed for inputting crop data.
Subsequently, following the research approach of Xing et al. [37] which entailed a sensitivity analysis of winter wheat parameters based on the extended Fourier amplitude sensitivity test (EFAST) method, this study utilized the growth data of winter wheat from the experimental area for the years 2021 to 2022 to adjust and calibrate the model parameters. The “trial and error method” was employed to debug the model parameters, aligning them within the parameter ranges provided in the model operation manual. Some of the calibrated model crop parameters are presented in Table 3. The model parameters were then validated using the growth data of winter wheat from the experimental area for the years 2022 to 2023.

3.2.3. Field Management Data and Soil Parameter Data

The crop management practices in the experimental area were conducted in accordance with the local winter wheat field management guidelines, with a particular emphasis on pest and weed control. The respective field management parameters and their values were compiled into a database. Soil parameters, including soil depth, soil particle size and mass fraction corresponding to each soil depth, and physical properties such as field water holding capacity, saturated water content, and permanent wilting point, were incorporated into the AquaCrop model to establish a soil parameter database. The details pertaining to these soil parameters are presented in Table 4.

3.3. Testing Methods for Evaluating AquaCrop Model Accuracy

The AquaCrop model is capable of simulating variations in canopy cover, aboveground biomass, yield, and soil water content throughout the crop reproductive period. Upon the establishment of the crop, meteorological and soil databases within the model, the calibration of model crop parameters, and the growth data of winter wheat from the experimental area for the period 2022 to 2023 was compared against the data generated from the model run. The following indices were employed to assess the accuracy of model calibration and verification in this study: the coefficient of determination (R2), normalized root mean square error (NRMSE), and consistency index (d).
The simulation effect is considered better as the values of d and R2 approach 1, with the optimal value being 1. An NRMSE value below 10% indicates excellent model simulation; a value between 10% and 20% indicates good simulation; a value between 20% and 30% suggests the simulation is generally effective; and a value greater than 30% denotes poor simulation. The model evaluation indices are computed using the formulas:
R 2 = i = 1 n S i S ¯ O i O ¯ 2 i = 1 n S i S ¯ 2 i = 1 n O i O ¯ 2
N R M S E = 1 n i = 1 n S i O i 2 × 100 / O ¯
d = 1 i = 1 n S i O i 2 i = 1 n S i O ¯ + O i O ¯ 2
where n is the number of measured values; S i is the analog value; O i is the measured value; S ¯ is the mean of the analog values; O ¯ is the mean of the measured values.

3.4. Calibration and Validation of the Model and Crop Data Calibration

3.4.1. Calibration and Validation of Canopy Cover and Aboveground Biomass

Following the sequential calibration approach by Vanuytrecht et al. [38], this paper first calibrated the canopy cover, followed by aboveground biomass, and finally the crop yield. Measurements from 2021 to 2022 were used to calibrate the model, measurements from 2022 to 2023 were used to validate the model, and the growth data of winter wheat in the experimental area from 2022 to 2023 were employed for validating the model parameters. Figure 4 and Figure 5, along with Table 5, illustrate the results of both processes. During calibration, for each water and fertilizer treatment of winter wheat, the coefficient of determination, R2, ranged between 0.92 to 0.98 for canopy cover and 0.96 to 0.98 for aboveground biomass. The NRMSE varied from 7.2 to 15.5 for canopy cover and 3.2 to 8.6 for aboveground biomass, while the consistency index, d, ranged from 0.91 to 0.96 for canopy cover and 0.96 to 0.99 for aboveground biomass.
The growth data from the experimental area for the period 2022 to 2023 was used for model validation. For each water-fertilized treatment of winter wheat, the R2 values ranged from 0.92 to 0.98 for canopy cover and 0.98 to 0.99 for aboveground biomass. The NRMSE values varied between 5.1 to 11.7 for canopy cover and 2.8 to 6.5 for aboveground biomass, while the d values ranged from 0.90 to 0.97 for canopy cover and 0.91 to 0.99 for aboveground biomass. The above results showed that the Aquacrop model can simulate the canopy cover, biomass, and yield of winter wheat under different irrigation and fertilization regimes better where the coefficient of determination R2, the standard root mean square error NRMSE, and the consistency index d between the simulated and measured values were within acceptable limits. See Figure 6, Figure 7 and Figure 8 for specific comparison charts.

3.4.2. Calibration of Yield, WUE, and PFP

WUE serves as an indicator of water absorption and utilization by winter wheat, while fertilizer PFP serves as a comprehensive index to assess fertilizer application and the local soil base nutrient situation. The accuracy of the model simulation for yield, WUE, and PFP data of winter wheat is demonstrated through comparisons of measured and simulated values, as shown in Table 6.
For each treatment over two years, the R2 values for yield were 0.99 and 0.99; NRMSE values were 1.78 and 0.64; and d values were 0.98 and 0.99, respectively. Similarly, for WUE, the R2 values were 0.95 and 0.96; NRMSE values were 1.78 and 2.63; and d values were 0.98 and 0.98, respectively. For PFP, the R2 values were 0.99 and 0.99; NRMSE values were 1.60 and 0.62; and d values were 0.99 and 0.99, respectively.
These statistical values underscore the high accuracy achieved in the calibration of yield, WUE, and PFP data for winter wheat. Thus, the AquaCrop model proves a viable tool for simulating various irrigation and fertilization scenarios for winter wheat to determine optimal irrigation and fertilization systems, aiming to enhance yield while conserving water and fertilizer resources.

4. Discussion

4.1. Simulation Scenarios

Overwintering irrigation refers to irrigation before winter to stabilize the ground temperature, increase soil moisture, and enable winter wheat to overwinter safely. The growth trajectory of winter wheat is significantly influenced by moisture levels, especially during the period extending from the jointing to the filling stage [39]. To further understand the impact of irrigation on the yield and other indices of winter wheat across different fertility phases, a multi-scenario simulation was conducted. Firstly, we calibrated the model parameters based on two years of field experimental data, and then explored the optimal irrigation and nitrogen application scheme based on this model.
In the scenario design, this study used the 2022–2023 meteorological data as the base data, uniformly setting the crop sowing date as October 15, determined the start of the overwintering period W as 15 December, the start of the greening period G as 15 March, the start of the nodulation period J as 15 April, and the start of the filling period F as 15 May. Three nitrogen application levels were selected: N1 (120 kg/ha), N2 (220 kg/ha), and N3 (320 kg/ha), along with three irrigation levels: W1 (60%θf), W2 (70%θf), and W3 (80%θf) (where θf is the field water holding capacity). Then, altering a single irrigation event for the growth period with an irrigation volume of 45 mm, resulted in a total of 36 irrigation and nitrogen application combination schemes, as shown in Table 7.

4.2. Scenario Analysis Results and Determination of Optimal Water and Fertilizer Management System

4.2.1. Scenario Analysis Results

Table 8 shows that, with Programs T1–T12 for instance, under the same irrigation level, the crop yield initially increased and then decreased with the rise in nitrogen application. Similarly, under a constant nitrogen application level, the crop yield initially increased and then decreased with the increase in irrigation volume.
Under identical irrigation and fertilization conditions, irrigation at varied fertility stages also resulted in certain yield differences. For instance, with Programs T1–T4, the corresponding yields were 5025, 5045, 5043, and 5028 kg/ha, respectively, indicating higher yields during the jointing or re-greening stages compared to the overwintering or filling stages. This highlights the influence of irrigation timing on yield.
Moreover, contrasting Programs T17–T20 with T29–T32 reveals that with sufficient irrigation, increasing irrigation volume diminishes the yield variance caused by irrigation at different fertility stages, without enhancing yield, but reducing WUE instead.

4.2.2. Impact of Water–Fertilizer Coupling on Yield, WUE, and PFP

The investigation of water–fertilizer coupling effects on yield, WUE, and PFP was carried out by selecting simulated scenarios yielding optimal results under the same water–fertilizer treatments, with significance analyses performed using IBM SPSS.
Significant differences (p < 0.05) in yield, WUE, and PFP of winter wheat under the same and different fertilizer treatments are shown in Table 9 and Figure 7. Both yield and WUE showed a tendency to increase and then decrease under the same irrigation conditions or the same fertilization conditions. Furthermore, Table 9 reveals the highly significant effects (p < 0.01) of irrigation and fertilization treatments on winter wheat yield and WUE, while irrigation treatment significantly affected PFP (p < 0.05). The interaction between irrigation and fertilization treatments also displayed significant effects on yield, WUE, and PFP of winter wheat (p < 0.05).

4.2.3. Determination of Optimal Water and Fertilizer Management System Based on Yield and WUE

To analyze the overall trend of irrigation and fertilization effects on yield, WUE, and PFP, a synthesis of scenario analysis results was performed, followed by regression analyses. Yield, WUE, and PFP were treated as dependent variables, while irrigation and fertilization served as independent variables. Regression models were established and expressed using binary quadratic regression equations, as presented in Equation (12).
z = a 0 + a 1 × x + a 2 × y + a 12 × x × y + a 11 × x 2 + a 22 × y 2
where z represents yield, WUE, and PFP, respectively; x and y are dimensionless variables; a0 is the constant term; a1 and a2 are the coefficients of the primary term; a12 is the coefficient of the interaction term; and a12 and a12 are the coefficients of the quadratic term [33].
Utilizing Origin 2021, response surface relationship plots concerning irrigation, fertilization, and yield, WUE, and PFP were drawn and analyzed, with the outcomes displayed in Figure 8. Examination of Table 10 and Figure 8 reveals that winter wheat yield is influenced by both water and fertilizer factors, which either promote or inhibit each other. The two-dimensional projection of the response surface in Figure 6 vividly shows that, at certain irrigation or fertilizer levels, the yield of winter wheat follows a pattern of initial increase followed by a decline. The regression equation shows that the regression coefficients for irrigation volume and fertilizer application volume were positive (41 and 30 respectively), indicating that within a certain range, both factors contribute to an increase in yield, although the impact of irrigation volume (x) is more significant than that of fertilizer application (y). At the point of maximum yield (7189 kg/ha), the irrigation volume was 317 mm, and the fertilizer application amount was 242 kg/ha, aligning with the findings from treatments T17–T20.
The varying irrigation and fertilization treatments had notable effects on both the WUE and PFP of winter wheat. Increasing the irrigation amount or fertilizer application initially led to an increase in WUE, which peaked at 14.022 kg/m3 (when irrigation was 281.61 mm and fertilizer was 236.12 kg/ha, comparable to treatments T17–T20), before declining. Similarly, an increment in irrigation volume initially led to an increase in PFP, which then descended, while an increase in fertilizer application rate caused a consistent decrease in PFP. Merely pursuing the maximum PFP would significantly diminish both yield and WUE, which contradicts the practical goal of achieving high yield alongside water conservation. Hence, the confidence interval at which both yield and WUE attain their optimum was taken into account.
Based on calculations, a yield exceeding 99% of the maximum value can be achieved when the irrigation volume ranges from 296 to 339 mm and the fertilizer application rate is between 229 to 261 kg/ha. Concurrently, a WUE nearing 99% of the maximum value is attainable with an irrigation volume of between 267 to 303 mm and a fertilizer application rate ranging from 210 to 268 kg/ha. Considering the objectives of high yield, water, and fertilizer conservation, an optimal interval was selected with an irrigation volume of 296 to 303 mm and a fertilizer application rate of 229 to 261 kg/ha. During this interval, the PFP ranged from 27 to 31 kg/kg, which was about 60% of the maximum PFP value, aligning more with the actual scenario. This interval closely corresponded with the outcomes of the W2N2 treatment in the two-year field experiment, suggesting that the regression equation, based on scenario design results, fitted well. This interval can thus offer a reference for the irrigation and fertilization system of winter wheat in the irrigation area of the Yellow River Basin.

5. Discussion

5.1. Model Calibration and Validation

In this study, the AquaCrop model was calibrated and validated using measured values of canopy cover, yield, above-ground biomass, and soil moisture content from two years of field experiments. Following validation, this model was employed to simulate the changing patterns of yield, WUE, and PFP of winter wheat under various irrigation schemes across different fertility periods. Scenario simulations were conducted to optimize water and nitrogen regimes. The results demonstrated that the AquaCrop model could accurately simulate the canopy cover and aboveground biomass of winter wheat under different irrigation and fertilization treatments. The R2, NRMSE, and d were all found to be within an acceptable range. Specifically, the ranges of R2, NRMSE, and d for canopy cover across various winter wheat treatments were 0.92–0.98, 7.2–15.5, and 0.91–0.96, respectively. For aboveground biomass, these ranges were 0.96–0.98, 3.2–8.6, and 0.96–0.99, respectively. These findings align with the study of Teng et al. [39], which also reported high accuracy in predicting canopy cover and aboveground biomass of winter wheat.
Additionally, for all treatments, the R2 and NRMSE values for winter wheat yield were 0.99 and 1.78, respectively; for WUE, they were 0.95 and 1.78, respectively. It was observed that during the calibration stage, the measured values of winter wheat yield were greater than the simulated values. However, during the validation stage, some simulated values of winter wheat yields were higher than the measured values in treatments after increasing water and fertilizer, which aligned with the findings of Guo et al. [40]. This discrepancy might be due to the increased model error when soil moisture or fertilizer stress appeared as a result of raising water and fertilizer application in the treatment.
Regarding the canopy cover of winter wheat, the simulated value of canopy cover before the overwintering period was larger compared to the measured value, similar to the results of the applicability evaluation of winter wheat by Chen et al. [30] This discrepancy could be associated with the fact that in reality, winter wheat almost halted growth during the overwintering period, keeping the canopy cover nearly unchanged, whereas the model continued to indicate an increase in canopy cover. During the crop’s aboveground biomass accumulation process, the simulated values of the model and the measured values followed a nearly similar trend. However, the model’s simulation of aboveground biomass during the reproductive period was higher than the measured values overall. This was unlike the results of Andarzian et al. [12], where the model simulated aboveground biomass lower than the measured values at the late reproductive stage. This discrepancy might be tied to the fact that the distribution of irrigation dates in the present experiments was more evenly spread out, leading to minor changes in soil moisture stress, as opposed to the experiments of Andarzian et al. [12], where water stress in the soil decreased with the reduction of nitrogen application in late irrigation. The model might not have adequately considered the water and nitrogen stress in the soil. Nonetheless, overall, the model simulation was satisfactory, and the calibrated model can be employed for local winter wheat yield prediction.

5.2. Scenario Simulation and Treatment Analysis

The scenario simulation in this paper using the calibrated model illustrated that under different irrigation and fertilization conditions, the T17–T20 scenarios exhibited the highest yield and WUE, aligning with the experimental results from the two-year field experiment, which showcased the highest yield and WUE under the W2N2 treatment. This study also deduced that under identical irrigation and fertilization conditions, the crop yield would be larger when irrigated at the jointing or re-greening stage as opposed to the overwintering or filling stage. Li et al. [34] demonstrated that irrigation at the jointing stage could significantly enhance leaf WUE and seed yield of winter wheat, with the yield increase attributed to a notable rise in the number of spikes. Similarly, Han et al. [41] found that irrigation at the jointing stage was more conducive to improving WUE and dry matter accumulation in winter wheat, thereby boosting the yield of winter wheat. This aligns well with the simulation results of irrigation at different fertility stages under the same irrigation and fertilization conditions. Concurrently, under the same irrigation conditions, winter wheat yield exhibited a trend of increasing and then decreasing with the rise in fertilizer application. Under identical fertilization conditions, a similar trend was observed with increasing irrigation, which is in line with the experimental results of Jia et al. [42], that under the same amount of irrigation, the dry matter accumulation of wheat increased with the rise in nitrogen application.

5.3. Water–Fertilizer Coupling Effects and Optimal Regimes

In the process of exploring the impacts of water–fertilizer coupling on yield, WUE, and PFP, it was found that both irrigation and fertilization treatments, along with the interaction between these two factors, exhibited highly significant effects (p < 0.01) on the yield and WUE of winter wheat. Additionally, the irrigation treatment and the interaction between the two factors demonstrated significant effects (p < 0.05) on the PFP of winter wheat. These findings were consistent with the study by Albrizio et al. [43] who showed that both water and nitrogen application were pivotal factors influencing the growth of winter wheat, and that appropriate irrigation coupled with increased nitrogen application can increase the yield of winter wheat.
In this paper, the employment of fitted coupled water–fertilizer regression equations facilitated the construction of response surface plots, which were helpful in identifying optimal irrigation and fertilization intervals. It was concluded that excessive irrigation could significantly reduce WUE, whereas an appropriate amount of irrigation could enhance the water adoption efficiency of the crop. This claim aligns with the findings of Chen et al. [44], who suggested that excessive soil water could impair nutrient uptake of the crop due to its impact on the respiration of crop roots. They further suggested that appropriate water stress could be more beneficial to crop growth and is more conducive to WUE.
In this investigation, it was identified that within the irrigation amount range of 296.13 to 303.39 mm, and the fertilizer application rate range of 229.67 to 261.94 kg/ha, the yield and WUE could attain more than 99% of the maximum value. These results aligned with the data acquired from two years of field experiments, thereby providing a valuable reference for the determination of irrigation and fertilizer application amounts in the irrigation area of the Yellow River Basin. Furthermore, these findings can be employed as a reference for fertilizer application strategies within the irrigation area of the Yellow River Basin, thereby contributing to the broader objective of optimized agricultural water management in this region.

6. Conclusions

Based on the extensive use of the Aquacrop model globally and the more accurate simulation results of the model in different regions of the world, it can be seen that the Aquacrop model has strong applicability to different regions and crops, and attention should be paid to the accuracy of the data inputted into the model, including climatic data, crop data, field management data, etc., when using the model.
The Aquacrop model has proven to be effective at simulating the canopy cover, yield, and aboveground biomass among other parameters of winter wheat in the irrigation area of the Yellow River Basin. The model calibration assessment indices have all satisfied the requisite criteria, demonstrating its reliability. Moreover, utilizing the calibrated localized model to predict the yield of winter wheat, when irrigated at varying fertility periods, aligns well with the outcomes of the two-year field trials. The simulation reveals that both yield and WUE follow a pattern of initial increase, followed by a decrease, with the escalation of irrigation volume and fertilizer application.
Furthermore, through response surface regression analysis, it was deduced that to attain 99% of the maximum value for yield and WUE, the irrigation volume should range between 296.13 mm and 303.39 mm, while the fertilizer application rate should be within the span of 229.67 kg/ha to 261.94 kg/ha. These findings are in line with the results from the field trials and serve as valuable reference points for determining the irrigation and fertilizer application rates for winter wheat in the irrigation area of the Yellow River Basin. Such insights not only confirm the findings of the field experiment but also extend a practical framework for optimized agricultural water management in this particular region.

Author Contributions

Conceptualization, T.L.; methodology, M.Y. and S.L.; validation, Y.L. (Yulong Liu); formal analysis, M.L.; data curation, Y.L. (Yuan Li); writing—original draft preparation, D.W.; writing—review and editing, W.Z.; visualization, K.L.; project administration, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the General Project of National Natural Science Foundation of China, No. 52079051, supported by the Program for Innovative Research Team (in Science and Technology) at the University of Henan Province (24IRTSTHN012) Key Scientific Research Project of Henan Province Colleges and Universities, No.: 22A570004 & 23A570006.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. Due to confidentiality of experimental data, the data are not publicly available.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhao, Y.; Song, H.; Mao, Y. Study on Food Security and Water Resource Security Countermeasures in the Yellow River Basin. People’s Yellow River 2011, 33, 47–49. [Google Scholar]
  2. Meng, Z.; Liu, A.; Wang, H. Research on the mathematical model of coupled water and fertiliser for water conservation and high yield of winter wheat and summer maize in a two-crop system. China Rural Water Conserv. Hydropower 1999, 11, 13–15+46. [Google Scholar]
  3. Zhou, J.; Ma, Y.; Wu, M.; Peng, Z.; Wang, Y.; Li, H.; Sheng, L. Water and Nitrogen Utilization and Biological Effects of Winter Wheat under Different Water and Fertilizer Measures. J. Irrig. Drain. 2019, 38, 36–41. [Google Scholar]
  4. Shang, W.; Zhang, Z.; Zheng, E.; Liu, M. Nitrogen-water coupling affects nitrogen utilization and yield of film-mulched maize under drip irrigation. J. Irrig. Drain. 2019, 38, 49–55. [Google Scholar]
  5. Liu, J.; Qiu, H.; Zhang, W.; Cai, J.; Wang, X.; Lv, M. The Effects of Irrigation and Biochar Amendment on Yield and Water and Nitrogen Use Efficiency of Winter Wheat. J. Irrig. Drain. 2021, 40, 59–65. [Google Scholar]
  6. Dai, J.; Li, R.; Li, C.; Lu, Y.; Zou, C. Simlation of Effects of Different Water and Fertilizer Treatments on Maize Growth with Aquacrop in Hetao Irrigated Area. J. Soil Water Conserv. 2021, 35, 312–319. [Google Scholar]
  7. Hui, L.L.Y.; CaiJia-bing, M.X.M. The Applicability and Application of AquaCrop Model. J. Irrig. Drain. 2011, 30, 28–33. [Google Scholar]
  8. Lianxi, W.; Jiansheng, W.; Qi, L.; Jiayi, G.; Hongxi, X. A Review on the Research and Application of AquaCrop Model. Adv. Earth Sci. 2015, 30, 1100–1106. [Google Scholar]
  9. Farahani, H.J.; Izzi, G.; Oweis, T.Y. Parameterization and Evaluation of the Aqua Crop Model for Full and Deficit Irrigated Cotton. Agron. J. 2009, 101, 469–476. [Google Scholar] [CrossRef]
  10. Araya, A.; Habtu, S.; Hadgu, K.M.; Kebede, A.; Dejene, T. Test of AquaCrop model in simulating biomass and yield of water deficient and irrigated barley (Hordeum vulgare). Agric. Water Manag. 2010, 97, 1838–1846. [Google Scholar] [CrossRef]
  11. Toumi, J.; Er-Raki, S.; Ezzahar, J.; Khabba, S.; Jarlan, L.; Chehbouni, A. Performance assessment of AquaCrop model for estimating evapotranspiration, soil water content and grain yield of winter wheat in Tensift Al Haouz (Morocco): Application to irrigation management. Agric. Water Manag. 2016, 163, 219–235. [Google Scholar] [CrossRef]
  12. Andarzian, B.; Bannayan, M.; Steduto, P.; Mazraeh, H.; Barati, M.E.; Barati, M.A.; Rahnama, A. Validation and testing of the AquaCrop model under full and deficit irrigated wheat production in Iran. Agric. Water Manag. 2011, 100, 1–8. [Google Scholar] [CrossRef]
  13. He, J. Study of the Effect of Different Climatic Years on Tobacco Yield Based on AquaCrop Modeling. Master’s Thesis, Nanjing University of Information Science and Technology, Nanjing, China, 2022. [Google Scholar]
  14. Ni, L.; Feng, H.; Ren, X. Applicable evaluation of crop model AquaCrop for summer maizeproduction in Loess Plateau Region. Agric. Res. Arid Areas 2015, 33, 40–45. [Google Scholar]
  15. Pan, H.; Dai, F.; Sui, R. Soil Moisture Modeling of Corn Full-Film Double Ridge Ditch Sowing Based on AquaCrop Model. Trans. Chin. Soc. Agric. Mach. 2022, 53, 78–86. [Google Scholar]
  16. Meng, Y. Simulation of Winter Wheat Growth in Hebei Piedmont Plain Based on AquaCrop Model. Master’s Thesis, Hebei Agricultural University, Baoding, China, 2021. [Google Scholar]
  17. Du, W.; He, X.; Shamaila, Z. Yield and Biomass Prediction Testing of AquaCrop Model for Winter Wheat. Trans. Chin. Soc. Agric. Mach. 2011, 42, 174–178+183. [Google Scholar]
  18. Zhu, C.; Xu, Y.; Huang, M. Simulation and Evaluation of Cycle Irrigation with Brackish and Fresh Water for Winter Wheat Based on AquaCrop Model. Trans. Chin. Soc. Agric. Mach. 2022, 53, 330–342. [Google Scholar]
  19. Cui, Y.; Lin, H.-H.; Xie, Y. Application study of crop yield prediction based on AquaCrop model in black soil region of Northeast China. Acta Agron. Sin. 2021, 47, 159–168. [Google Scholar] [CrossRef]
  20. 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]
  21. Heng, L.K.; Hsiao, T.; Evett, S.; Howell, T.; Steduto, P. Validating the FAO AquaCrop model for irrigated and water deficient field maize. Agron. J. 2009, 101, 488–498. [Google Scholar] [CrossRef]
  22. Khoshravesh, M.; Mostafazadeh-Fard, B.; Heidarpour, M.; Kiani, A. AquaCrop model simulation under different irrigation water and nitrogen strategies. Water Sci. Technol. 2012, 67, 232–238. [Google Scholar] [CrossRef]
  23. Karunaratne, A.S.; Azam-Ali, S.N.; Izzi, G.; Steduto, P. Calibration and validation of FAO-AquaCrop model for irrigated and water deficient bambara groundnut. Exp. Agric. 2011, 47, 509–527. [Google Scholar] [CrossRef]
  24. Zeleke, K.T.; Luckett, D.; Cowley, R. Calibration and testing of the FAO AquaCrop model for canola. Agron. J. 2011, 103, 1610–1618. [Google Scholar] [CrossRef]
  25. Margarita, G.V.; Elías, F. Combining the simulation crop model AquaCrop with an economic model for the optimization of irrigation management at farm level. Eur. J. Agron. 2012, 36, 21–31. [Google Scholar]
  26. Nirman, S.; Dirk, R.; Shrawan, K.S. Strategies to improve coreal production in the Terai region (Nepal) during dry season: Simulations with Aquacrop. Procedia Environ. Sci. 2013, 19, 767–775. [Google Scholar]
  27. Shao, D.G.; Le, Z.H.; Xu, B.L.; Hu, N.J.; Tian, Y.N. Optimization of irrigation scheduling for organic rice based on AquaCrop. Trans. Chin. Soc. Agric. Eng. 2018, 34, 114–122. [Google Scholar]
  28. Wang, Q.J.; He, H.; Li, L. Research on Soybean Irrigation Schedule Based on AquaCrop Model. Sci. Agric. Sin. 2022, 55, 3365–3379. [Google Scholar]
  29. Song, J.; Li, J.; Yang, Q.H.; Mao, X.; Yang, J.; Wang, K. Multi-objective optimization and its applicationon irrigation scheduling based on AquaCrop and NSGA-II. J. Hydraul. Eng. 2018, 49, 1284–1295. [Google Scholar]
  30. Chen, C. Irrigation and Nitrogen Management for Winter Wheat-Summer Maize in Guanzhong Region Based on AquaCrop Model. Master’s Thesis, College of Water Resources and Architectural Engineering, Northwest A&F University, Xianyang, China, 2019. [Google Scholar]
  31. Ding, C. Water and Nitrogen Interaction and Optimization of Irrigation System under Degradable Film Mulching in Dripped Sandy Farmland. Master’s Thesis, Inner Mongolia Agricultural University, Hohhot, China, 2018. [Google Scholar]
  32. Zhang, F.C.; Yan, F.L.; Fan, X.K.; Li, G.; Liu, X.; Lu, J.S.; Wang, Y.; Ma, M. Effects of irrigation and fertilization levels on grain yield and water-fertilizer use efficiency of drip-fertigation spring maize in Ningxia. Trans. Chin. Soc. Agric. Eng. 2018, 34, 111–120. [Google Scholar]
  33. Zhang, J.; Lv, L.; Zhang, L. Application of a quantitative method for crop water-fertiliser coupling type in water and nitrogen allocation for winter wheat in North China. China Agric. Sci. 2019, 52, 2997–3007. [Google Scholar]
  34. Li, Q.; Shen, J.; Zhao, D. Effect of irrigation frequency on yield and leaf water use efficiency of winter wheat. Trans. CSAE 2011, 27, 33–36. [Google Scholar]
  35. Raes, D.; Steduto, P.; Hsiao, T.C.; Fereres, E. AquaCrop the FAO crop model to simulate yield response to water: II. Main Algorithms and Software Description. Agron. J. 2009, 101, 438–447. [Google Scholar] [CrossRef]
  36. 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]
  37. Xing, H.; Xiang, S.; Xu, X.; Feng, H.; Yang, G.; Chen, Z. Global Sensitivity Analysis of AquaCrop Crop Model Parameters Based on EFAST Method. Sci. Agric. Sin. 2017, 50, 64–76. [Google Scholar]
  38. Vanuytrecht, E.; Raes, D.; Steduto, P.; Hsiao, T.C.; Fereres, E.; Heng, L.K.; Vila, M.G.; Moreno, P.M. AquaCrop: FAO’s crop water productivity and yield response model. Environ. Model. Softw. 2014, 62, 351–360. [Google Scholar] [CrossRef]
  39. Teng, X.W.; Dong, Y.S.; Shen, J.X.; Meng, L.M.; Feng, H.K. Winter Wheat Irrigation Simulation in Arid Area Based on AquaCrop Model. Sci. Agric. Sin. 2015, 48, 4100–4110. [Google Scholar]
  40. Guo, D. Growth Simulation and Regional Irrigation Schedule Optimization for Winter Wheat and Summer Maize Based on AquaCrop Model in the Fenwei Plain. Ph.D. Thesis, College of Water Resources and Architectural Engineering, Northwest A&F University, Xianyang, China, 2021. [Google Scholar]
  41. Han, H.F.; Li, Q.Q.; Dong, B.D.; Liu, M.Y. Effects of irrigation frequency and stages on grain yield and quality characteristics of winter wheat. Acta Ecol. Sin. 2010, 30, 1548–1555. [Google Scholar]
  42. Jia, Z.-L.; Bai, Y.-G.; Xiao, J. Coupling Effect of Water and Fertilizer on Winter Wheat under Drip Irrigation in Northern Xinjiang. J. Irrig. Drain. 2014, 33, 12–15. [Google Scholar]
  43. Albrizio, R.; Todorovic, M.; Matic, T.; Stellacci, A.M. Comparing the interactive effects of water and nitrogen on durum wheat and barley grown in a Mediterranean environment. Field Crops Res. 2010, 115, 179–190. [Google Scholar] [CrossRef]
  44. Chen, K.; Zhao, J.; Huang, H. Effect of Different Irrigation Quota on Water Consumption Characteristics and Yeld of Wheat. J. Irrig. Drain. 2017, 36, 65–68, 84. [Google Scholar]
Figure 1. Location of the test area.
Figure 1. Location of the test area.
Agronomy 14 00110 g001
Figure 2. Maximum and minimum temperature data across the two-year.
Figure 2. Maximum and minimum temperature data across the two-year.
Agronomy 14 00110 g002
Figure 3. Rainfall and reference crop evapotranspiration data across the two-year.
Figure 3. Rainfall and reference crop evapotranspiration data across the two-year.
Agronomy 14 00110 g003
Figure 4. Comparison of analog and measured values of canopy cover for (a) calibration and (b) validation phases.
Figure 4. Comparison of analog and measured values of canopy cover for (a) calibration and (b) validation phases.
Agronomy 14 00110 g004
Figure 5. Comparison of analog and measured values of aboveground biomass for (a) calibration and (b) validation phases.
Figure 5. Comparison of analog and measured values of aboveground biomass for (a) calibration and (b) validation phases.
Agronomy 14 00110 g005
Figure 6. Comparison of analog and measured values of soil moisture content for (a) calibration and (b) validation phases.
Figure 6. Comparison of analog and measured values of soil moisture content for (a) calibration and (b) validation phases.
Agronomy 14 00110 g006
Figure 7. Effect of water–fertilizer coupling on (a) yield, (b) WUE, and (c) PFP. Note: Different lowercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying fertilization treatment levels within the same irrigation treatment level; different uppercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying irrigation treatment levels within the same fertilization treatment level.
Figure 7. Effect of water–fertilizer coupling on (a) yield, (b) WUE, and (c) PFP. Note: Different lowercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying fertilization treatment levels within the same irrigation treatment level; different uppercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying irrigation treatment levels within the same fertilization treatment level.
Agronomy 14 00110 g007
Figure 8. Response of irrigation and fertilization on the effect of (a) yield, (b) WUE, and (c) PFP.
Figure 8. Response of irrigation and fertilization on the effect of (a) yield, (b) WUE, and (c) PFP.
Agronomy 14 00110 g008
Table 1. Water and fertilizer treatment configurations.
Table 1. Water and fertilizer treatment configurations.
TreatmentWater Control LevelNumber of IrrigationsTotal Irrigation (mm)Irrigation Quota (mm)Nitrogen Application (kg/ha)Base Fertilizer (kg/ha)Top-Dressing Fertilizer (kg/ha)
W1N160%θf4454 × 4567.5 + 52.567.552.5
W1N260%θf67.5 + 152.5152.5
W1N360%θf67.5 + 252.5252.5
W2N170%θf66 × 4567.5 + 52.552.5
W2N270%θf67.5 + 152.5152.5
W2N370%θf67.5 + 252.5252.5
W3N180%θf99 × 4567.5 + 52.552.5
W3N280%θf67.5 + 152.5152.5
W3N380%θf67.5 + 252.5252.5
Table 2. Fertilizer stress parameters and their calibration values.
Table 2. Fertilizer stress parameters and their calibration values.
Nitrogen TreatmentParameter Calibration ValuesParameter Calibration Results
Relative Biomass
Brel
(%)
Maximum Canopy Cover under Fertilizer Stress
(%)
Degree of Canopy DecayReduced Maximum Canopy Cover (%)Decrease in Canopy Growth Coefficient (%)Average Canopy Reduction (%)Standardized Crop Water Productivity Reduction
(%)
N1 (120 kg/ha)7181.7Small1440.4034
N2 (220 kg/ha)8290.2Small530.2828
N3 (320 kg/ha)7584.5Small1130.3034
Table 3. Selected calibrated parameters of the Aquacrop model.
Table 3. Selected calibrated parameters of the Aquacrop model.
Model ParameterNotationCalibrated ValueDefault Parameter ValuesUnit
Reference harvest indexHI04845–50%
Base temperatureTbase00°C
Upper temperatureTupper2626°C
Canopy growth coefficientCGC3.601–5%/d
Canopy attenuation coefficientCDC0.3130.1–0.5%/d
Maximum effective root depthZX10.8–1.5m
Standardized water productivityWP*1711–22g/m3
Sowing to seedling emergence time/277150–280GDD
Seeding to maximum canopy cover time/1734Time to emergence +1000–2000GDD
Sowing to senescence time/1815Time to emergence +1000–2000GDD
Sowing to ripening time/2425Time to emergence +1500–2900GDD
Upper threshold for the effect of water stress on canopy growthPexp,upper0.200.14–0.26/
Lower threshold for the effect of water stress on canopy growthPexp,lower0.550.45–0.84/
Upper threshold for the effect of water stress on stomatal closurePsto,upper0.650.60–0.67/
Upper threshold for the effect of water stress on canopy senescencePsen,upper0.701.75–3.25/
Table 4. Physical properties of soil at 0–100 mm depth in the experimental area.
Table 4. Physical properties of soil at 0–100 mm depth in the experimental area.
Soil Depth (mm)Soil Particle Size Mass Fraction (%)Physical Parameter
SandSiltClayPermanent Wilting Point (%)Field Water Holding Capacity (%)Saturated Water Content (%)
0–200.170.640.1920.433.336.5
20–400.110.650.2420.935.538
40–600.090.650.2619.735.440
60–800.080.650.2721.135.638
80–1000.050.610.3421.435.738
Table 5. Indices of model calibration and validation accuracy.
Table 5. Indices of model calibration and validation accuracy.
Model Calibration and ValidationWater and Fertilizer TreatmentCanopy Cover (%)Aboveground Biomass (t/ha)Soil Moisture Content (%)
R2NRMSEdR2NRMSEdR2NRMSEd
CalibrationW1N10.9215.50.910.968.60.960.922.40.96
W1N20.987.90.970.967.70.960.964.90.90
W1N30.9215.30.910.983.40.970.882.90.94
W2N10.969.40.960.964.60.980.723.70.90
W2N20.987.20.960.9860.960.943.90.92
W2N30.9810.70.960.984.40.980.973.60.90
W3N10.9711.70.940.963.70.990.962.90.95
W3N20.9811.40.950.984.20.970.963.90.91
W3N30.9811.90.940.973.20.980.943.40.94
ValidationW1N10.9211.70.900.993.90.980.986.60.92
W1N20.9410.40.920.993.40.980.955.70.93
W1N30.948.50.930.994.60.970.975.10.96
W2N10.968.70.930.984.10.980.9730.94
W2N20.986.30.960.985.10.970.982.70.91
W2N30.986.10.960.984.60.910.9920.95
W3N10.987.50.940.992.80.990.983.20.94
W3N20.985.10.970.986.50.950.984.10.91
W3N30.9870.950.993.70.980.993.10.95
Table 6. Accuracy test indices for water and fertilizer indicators.
Table 6. Accuracy test indices for water and fertilizer indicators.
YearTreatmentYield (kg/ha)R2NRMSEdWUE (kg/m3)R2NRMSEdPFP (kg/kg)R2NRMSEd
MeasuredAnalogMeasuredAnalogMeasuredAnalog
2021–2022W1N1616860120.991.780.9813.2712.990.951.780.9825.7025.050.991.610.99
W1N26586656613.9013.9727.4427.36
W1N36371634313.1613.3726.5526.43
W2N17699760514.0314.4317.5017.28
W2N28138805314.9315.3518.5018.30
W2N37900768014.2714.3617.9617.45
W3N17538747712.1412.1111.7811.68
W3N27842773012.3412.1612.2512.08
W3N37670763312.5112.6111.9811.93
2022–2023W1N1500250300.990.640.999.549.520.990.750.9941.92 41.69 0.990.610.99
W1N25912590210.6910.6649.18 49.27
W1N35573554310.0610.0346.19 46.45
W2N16330630310.3710.2828.65 28.78
W2N27094700311.4611.3931.83 32.25
W2N36650662410.8810.7830.11 30.23
W3N1598159999.219.3318.75 18.69
W3N26545659110.1110.0620.60 20.46
W3N362326224.9.819.6919.45 19.48
Table 7. Simulation scenarios for nitrogen application by flooding.
Table 7. Simulation scenarios for nitrogen application by flooding.
Program NumberWater and Fertilizer TreatmentFlooding TimeIrrigation Quota (mm)Nitrogen Application (kg/hm2)Program NumberWater and Fertilizer TreatmentFlooding TimeIrrigation Quota (mm)Nitrogen Application (kg/hm2)
T1W1N1W180120T19W2N2J270220
T2W1N1G180120T20W2N2F270220
T3W1N1J180120T21W2N3W270320
T4W1N1F180120T22W2N3G270320
T5W1N2W180220T23W2N3J270320
T6W1N2G180220T24W2N3F270320
T7W1N2J180220T25W3N1W405120
T8W1N2F180220T26W3N1G405120
T9W1N3W180320T27W3N1J405120
T10W1N3G180320T28W3N1F405120
T11W1N3J180320T29W3N2W405220
T12W1N3F180320T30W3N2G405220
T13W2N1W270120T31W3N2J405220
T14W2N1G270120T32W3N2F405220
T15W2N1J270120T33W3N3W405320
T16W2N1F270120T34W3N3G405320
T17W2N2W270220T35W3N3J405320
T18W2N2G270220T36W3N3F405320
Table 8. Simulated values of winter wheat yield, WUE, and PFP under scenario analysis.
Table 8. Simulated values of winter wheat yield, WUE, and PFP under scenario analysis.
Program NumberYield
(kg/hm2)
WUE
(kg/m3)
PFP
(kg/kg)
Program NumberYield
(kg/hm2)
WUE
(kg/m3)
PFP
(kg/kg)
T1502510.7741.88T19699614.1331.80
T2504510.8942.04T20699114.1131.78
T3504310.8742.03T21652813.1720.40
T4502810.8041.90T22661413.3820.67
T5583811.9226.54T23661113.3720.66
T6590612.0826.85T24658213.3020.57
T7590112.0726.82T25579610.1548.31
T8586812.0026.67T26580410.2148.37
T9551411.2917.23T27579910.2048.33
T10554211.3617.32T28579810.1748.32
T11553911.3617.31T29658411.3329.93
T12552811.3317.28T30659011.3629.96
T13603412.0950.28T31658911.3629.95
T14606712.3050.56T32658811.3529.95
T15606612.2650.55T33622110.9419.44
T16605412.1950.45T34622310.9719.45
T17698514.0731.75T35622310.9719.45
T18699814.1531.81T36622210.9519.44
Table 9. Effect of water–fertilizer coupling on yield, WUE, and PFP.
Table 9. Effect of water–fertilizer coupling on yield, WUE, and PFP.
Program NumberFlooding TreatmentFertilizer TreatmentYield
(kg/hm2)
WUE
(kg/m3)
PFP
(kg/kg)
T2W1N15035 ± 10 cC10.83 ± 0.057 cB41.96 ± 0.086 aC
T6W1N25878 ± 31 aC12.02 ± 0.074 aB26.72 ± 0.144 bC
T10W1N35530 ± 12 bC11.34 ± 0.033 bB17.28 ± 0.04 cC
T14W2N16055 ± 15 cA12.21 ± 0.092 cA50.46 ± 0.128 aA
T18W2N26992 ± 5 aA14.12 ± 0.034 aA31.79 ± 0.027 bA
T22W2N36584 ± 39 bA13.31 ± 0.097 bA20.58 ± 0.125 cA
T26W3N15799 ± 3 cB10.18 ± 0.028 cC48.33 ± 0.027 aB
T30W3N26588 ± 2 aB11.35 ± 0.014 aC29.95 ± 0.012 bB
T34W3N36222 ± 1 bB10.96 ± 0.015 bC19.44 ± 0.005 cB
Flooding treatment*****
Fertilizer treatment******
Flooding treatment × Fertilizer treatment*****
Note: Different lowercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying fertilization treatment levels within the same irrigation treatment level; different uppercase letters indicate significant differences (p < 0.05) in winter wheat parameters at varying irrigation treatment levels within the same fertilization treatment level. Symbols * and ** denote significance levels (p < 0.05) and (p < 0.01), respectively.
Table 10. Regression relationship between water and fertilizer supply and yield, WUE, and PFP.
Table 10. Regression relationship between water and fertilizer supply and yield, WUE, and PFP.
Dependent Variable zRegression EquationR2Combination When Z Is the Maximum Value
X (mm)Y (kg/hm2)Zmax
Yield (kg/ha)Z = −2917.942 + 41.090x + 30.087y − 0.00196xy − 0.06412x2 − 0.06156y20.9959317.900242.5807189.340
WUE (kg/m3)Z = −5.267 + 0.096x + 0.047y + 5.100 × 10−6xy − 1.705 × 10−4x2 − 1.012 × 10−4y20.9858281.610236.12014.022
PFP (kg/kg)Z = 41.392 + 0.213x − 0.271y − 8.303× 10−5xy − 3.349 × 10−4x2 + 3.530 × 10−4y20.9904332.000120.00050.547
Note: x represents the amount of irrigation water, and y represents the amount of fertilizer applied.
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

Wang, S.; Wang, D.; Liu, T.; Liu, Y.; Luo, M.; Li, Y.; Zhou, W.; Yang, M.; Liang, S.; Li, K. Simulation of Winter Wheat Growth Dynamics and Optimization of Water and Nitrogen Application Systems Based on the Aquacrop Model. Agronomy 2024, 14, 110. https://doi.org/10.3390/agronomy14010110

AMA Style

Wang S, Wang D, Liu T, Liu Y, Luo M, Li Y, Zhou W, Yang M, Liang S, Li K. Simulation of Winter Wheat Growth Dynamics and Optimization of Water and Nitrogen Application Systems Based on the Aquacrop Model. Agronomy. 2024; 14(1):110. https://doi.org/10.3390/agronomy14010110

Chicago/Turabian Style

Wang, Shunsheng, Diru Wang, Tengfei Liu, Yulong Liu, Minpeng Luo, Yuan Li, Wang Zhou, Mingwei Yang, Shuaitao Liang, and Kaixuan Li. 2024. "Simulation of Winter Wheat Growth Dynamics and Optimization of Water and Nitrogen Application Systems Based on the Aquacrop Model" Agronomy 14, no. 1: 110. https://doi.org/10.3390/agronomy14010110

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