Next Article in Journal
Antioxidative and Metabolic Responses in Canola: Strategies with Wood Distillate and Sugarcane Bagasse Ash for Improved Growth under Abiotic Stress
Previous Article in Journal
Going or Returning to Nature? Wild Vegetable Uses in the Foraging-Centered Restaurants of Lombardy, Northern Italy
Previous Article in Special Issue
Modeling Nitrogen Fate and Water and Nitrogen Use Efficiencies under Different Greenhouse Vegetable Production Systems Using the WHCNS-Veg Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing the Fates of Water and Nitrogen on an Open-Field Intensive Vegetable System under an Expert-N System with EU-Rotate_N Model in North China Plain

1
College of Land Science and Technology, China Agricultural University, Key Laboratory of Arable Land Conservation (North China), Ministry of Agriculture and Rural Affairs, Beijing 100193, China
2
Agricultural Information Institute of Chinese Academy of Agricultural Science, Beijing 100081, China
3
College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2024, 13(15), 2150; https://doi.org/10.3390/plants13152150
Submission received: 20 April 2024 / Revised: 6 July 2024 / Accepted: 1 August 2024 / Published: 3 August 2024
(This article belongs to the Special Issue Water and Nitrogen Management in Soil-Crop System II)

Abstract

:
Nitrate leaching, greenhouse gas emissions, and water loss are caused by conventional water and fertilizer management in vegetable fields. The Expert-N system is a useful tool for recommending the optimal nitrogen (N) fertilizer for vegetable cultivation. To clarify the fates of water and N in vegetable fields, an open-field vegetable cultivation experiment was conducted in Dongbeiwang, Beijing. This experiment tested two irrigation treatments (W1: conventional and W2: optimal) and three fertilizer treatments (N1: conventional, N2: optimal N rate by Expert-N system, and N3: 80% optimal N rate) on cauliflower (Brassica oleracea L.), amaranth (Amaranthus tricolor L.), and spinach (Spinacia oleracea L.). The EU-Rotate_N model was used to simulate the fates of water and N in the soil. The results indicated that the yields of amaranth and spinach showed no significant differences among all the treatments in 2000 and 2001. However, cauliflower yield under the W1N2 and W1N3 treatments obviously reduced in 2001. Compared with the W1 treatment, W2 reduced irrigation amount by 27.9–29.8%, water drainage by over 76%, increased water use efficiency by 5–17%, and irrigation water use efficiency by 29–45%. Nitrate leaching was one of the main pathways in this study, accounting for 8.4% of the total N input; compared to N1, the input of fertilizer N under the N2 and N3 treatments decreased by over 66.5%, consequently reducing gaseous N by 48–72% and increasing nitrogen use efficiency (NUE) by 17–37%. Additionally, compared with the W1 treatments, gaseous N loss under the W2 treatments was reduced by 18–26% and annual average NUEs increased by 22–29%. The highest annual average NUEs were under W2N3 (169.6 kg kg−1) in 2000 and W2N2 (188.0 kg kg−1) in 2001, respectively. We found that optimizing fertilizer management allowed subsequent crops to utilize residual N in the soil. Therefore, we suggest that the W2N3 management should be recommended to farmers to reduce water and N loss in vegetable production systems.

1. Introduction

Vegetable crops are high-priced cash crops with short growing periods and high returns. They require large amounts of water and nitrogen (N) to grow rapidly and attain high yields [1]. Therefore, vegetable fields are particularly prone to water loss and N pollution [2,3]. The North China Plain (NCP) serves as a significant center for vegetable cultivation in China, where the shortage of water resources becomes the primary constraint. Previous research have proposed that conventional irrigation and N fertilizer practices lead to lower nitrogen use efficiency (NUE) and higher residual N in the soil profile in this area [4,5,6]. Furthermore, with the increasing areas of vegetable cultivation, the less-developed root systems of vegetable crops have caused large amounts of nitrate leaching, and resulted in the contamination of groundwater [7,8].
To obtain the best management practices for intensive vegetable systems, many water and N fertilizer practices were applied to improve water use efficiency (WUE) and NUE. Compared with conventional management, optimized N application rates significantly decreased the nitrate leaching loss by 47% without negative impacts on N uptake and tomato yields [8]. A meta-analysis revealed that reducing the input rates of N fertilizers or irrigation water brought advantages for maintaining vegetable yields, increasing NUE, decreasing N2O emissions, and reducing reactive N released into the environment [9]. Another vegetable crop rotation experiment indicated that a combination of drip irrigation and high irrigation water input could achieve a balance of high yields and low environmental risks [10]. N expert systems, such as the KNS system (Kulturbegleitende Nmin Software V4.6) are useful tools for managing N fertilizer applications in vegetable fields and have been widely used in Europe [11,12]. Later, a modified N-Expert system for vegetable cultivation was proposed by Chen et al. [13], and results indicated that this N recommendation system could reduce N fertilizer input and N residual in soil in the NCP. However, the fate of N under the N-Expert system remains unclear.
After the popularization of computers, human knowledge of agricultural systems was captured in computer models to guide farmers in making decisions on the different aspects of crop management [14]. In most cases, mechanistic models perform well in simulating the process of nitrate leaching, crop N uptake, and crop production quality [15,16,17]. For example, the NLEAP model was developed to estimate the potential of nitrate leaching in crop fields [18,19,20]. Alternatively, based on the N_ABLE model, the EU-Rotate_N model is a more powerful tool for simulating water and N dynamics in soil/plant systems [21]. The EU-Rotate_N model has been successfully used to simulate water and N dynamics and vegetable growth in some studies in Europe and China [22,23,24,25]. Their results suggested that the EU-Rotate_N model worked well in tracking water and N fates under different irrigation and fertilizer management. However, few studies have focused on the application of the EU-Rotate_N model for various continuous vegetable cultivation systems.
Open-field vegetable cultivation is one of the main vegetable-planting patterns in the NCP. The issues of water drainage, nitrate leaching, and gaseous N release in greenhouse vegetable planting systems have already been extensively explored in previous works [8,23,26,27], but those processes are more complex and less certain under open-field vegetable cultivation systems. Moreover, in the NCP, three varieties of vegetables are typically cultivated consecutively within a single year, leading to an increased frequency of water and N applications to the soil. Based on the above studies, we proposed the following hypotheses: (i) the improved N expert systems could save more water and fertilizer resources and (ii) optimized water and N management practices will significantly reduce water drainage and N loss, and improve WUE and NUE compared to the conventional practices in the NCP.
Therefore, our aims are (i) to calibrate and validate the EU-Rotate_N model with water and soil mineral N contents measured under the conventional and N-Expert systems on an open-field vegetable-growing system and (ii) to assess the fates of water and N, and WUE and NUE under conventional and optimized water and N management in an intensive vegetable system.

2. Results

2.1. Model Calibration and Validation

The soil water content, soil nitrate concentration, dry weight, N uptake, and yields of vegetables collected from the W1N1 treatment were used to calibrate the EU-Rotate_N model. The irrigation and N fertilizer applications, as detailed in Table S1, were incorporated into the model. Calibration was achieved through a ‘trial and error’ approach, adjusting parameters to align the simulated outcomes with the measured data [8,27,28,29]. The calibration was processed in the following order: soil water content, nitrate concentration, crop N uptake, vegetable yield, and crop dry matter [23,30]. Crop coefficients and luxury N consumption coefficients for the three vegetables were adjusted based on measured vegetable yield and N uptake (Table S2). Adjustments to N transformation parameters were guided by the comparison of the simulated and actual soil nitrate concentrations. The modified parameter values are listed in Table S2. Additionally, Figure 1 and Figure S2 compare the simulated and actual measurements of the soil water content, soil nitrate concentration, and dry matter weight for the W1N1 treatment, demonstrating a commendable correlation between the simulated and observed values.
After the calibration, the model was validated by the datasets collected from the other treatments (W1N2, W1N3, W2N1, W2N2, and W2N3). The comparisons of the simulated and measured soil water content and nitrate concentrations at different soil layers for all the validation treatments are presented in Figure 1, Figure 2, Figure 3 and Figure 4. Moreover, statistical indices for the simulated soil water content and nitrate concentration for these treatments are illustrated in Table 1.
As shown in Figure 1, Figure 2 and Figure 3, the surface soil layer (0–15 cm) exhibited frequent fluctuations in the water content, leading to elevated root mean square error (RMSE) values (Table 1). However, the RMSE values for the subsequent soil layers ranged between 0.01 and 0.02. Across all the soil layers, the values of Nash/Sutcliffe efficiency coefficient (E) were more than 0.36, and the values of d ranged from 0.83 to 0.98. These indices suggest a strong concordance between the simulated soil water content and the actual observations [31,32]. As shown in Figure 4, most of the simulated values of the soil nitrate concentration were in good agreement with the observed values. Furthermore, the values of E and agreement index (d) were both above 0.87, indicating a high level of accuracy in the simulation results.
The correlation coefficient between the simulated and measured crop N uptake was 0.965, and the p-value was significant (p < 0.001) (Figure 5). The results shown in Figures S1 and S2 indicate that the simulated crop dry matter weight agreed well with the measured values.
In conclusion, the EU-Rotate_N model performed well in predicting the soil water content, soil nitrate concentration, N uptake, and vegetable growth in this study. Therefore, the model could be used to predict the fate of water and N in vegetable production systems [25,33].

2.2. Vegetable Yield, Water Balance and Water Use Efficiency

Different irrigation and fertilizer practices did not result in significant differences in vegetable yield in 2000 and 2001, except for cauliflower (p < 0.05, Table 2 and Table 3). The average yield for three vegetables in 2000 and 2001 is as follows: cauliflower, 20,450 and 17,800 kg ha−1; amaranth, 13,933 and 19,700 kg ha−1; and spinach, 33,933 and 29,933 kg ha−1. Over the two years, the W1N1 treatment produced the highest vegetable yield.
The water balance within the 0.9 m soil profile across all the treatments for the years 2000 and 2001 is shown in Table 2 and Table 3. Rainfall and irrigation were the major water resources in this experiment. Notably, the precipitation in 2000 was higher than that in 2001, with the majority occurring during the amaranth-growing season. The total water inputs were as follows: for treatment W1, 769.5 mm in 2000 and 665.2 mm in 2001; for treatment W2, 627.8 mm in 2000 and 509.9 mm in 2001. Evapotranspiration (ET) was the predominant form of water consumption for all the treatments, accounting for 70.5% to 96.2% of the total water input over the two-year period. The ET difference between the different irrigation treatments was so obvious, but under the same irrigation treatments, the difference between different fertilization treatments was very small. Optimizing water input was found to significantly reduce ET within identical fertilizer treatments, although varying fertilizer regimes induced minimal alterations in ET under consistent irrigation management. For instance, in 2000, the total ET under the W1 treatment ranged from 542.7 to 562.0 mm, which decreased to a range of 489.0 mm to 501.3 mm under treatment W2. Among the three vegetables, the order of ET in 2000 was as follows: cauliflower > spinach > amaranth, while the ranking in 2001 was cauliflower > amaranth > spinach. The annual ET was primarily influenced by the ET during the cauliflower growth period.
The simulation results revealed that water drainage accounted for a minor fraction of the total water output in this study. The amount of water drainage was more substantial in 2000 than in 2001, with the highest annual drainage of 168.3 mm and 56.4 mm occurring under the W1N3 treatment in the respective years. Furthermore, the peaks of water drainage mainly occurred after irrigation or heavy rainfall during different vegetable growth seasons in 2000 and 2001, respectively (Figure 6). The study further demonstrated that optimizing irrigation schedules markedly reduced water drainage. The water drainage under the W2 treatment decreased by 74–78% in 2000 and 82–98% in 2001 compared to the W1 treatment, respectively. Moreover, the simulation results indicated that a small quantity of water remained within the 0–0.9 m soil profile after spinach harvest in both years under all the treatments. Conversely, the cauliflower exhibited signs of mild dehydration during the two experimental years.
Regardless of the irrigation and fertilizer management employed, water use efficiencies (WUEs) and irrigation water use efficiencies (IWUEs) in the spinach-growing season were the highest in both two years. However, the WUEs and IWUEs exhibited variable responses among the three vegetable cultivations under differing water and fertilizer regimes. For instance, the highest WUEs of cauliflower were observed under the W2N1 treatment in 2000 and 2001, reaching 7.6 and 7.5 kg m−3, respectively. In the case of amaranth, the W2N3 (in 2000) and W2N1 (in 2001) treatments produced optimal WUEs, reaching 16.3 and 19.9 kg m−3, respectively. Similarly, spinach achieved its maximum WUEs under the same treatments as those that benefited amaranth, reaching 29.4 and 26.9 kg m−3, respectively.
Additionally, the maximum annual WUE was recorded under the W2N2 treatment, while the W2N3 treatment produced the second-highest IWUE in 2000. In 2001, the W2N1 treatment observed the highest annual WUE and IWUE. When considering the same N fertilizer input treatments, there was a negative correlation with irrigation input. For example, the amount of irrigation under the W1N1 treatment increased by 27.9–29.7% compared with W2N1. Consequently, the annual average WUE and IWUE increased by 5–17% and 29–45%, respectively. Under the same irrigation treatments, the patterns of the annual average WUEs and IWUEs in response to the reduced N fertilizer input were ambiguous in 2000. However, the annual average WUEs and IWUEs slightly declined with the reductions in N fertilizer input in 2001. Therefore, alterations in the annual average WUEs and IWUEs were sensitive to irrigation management.

2.3. Nitrogen Balance and Nitrogen Use Efficiency

The N balances in the 0.9 m soil profile for all the treatments in both two years are summarized in Table 4 and Table 5. Besides manure, fertilizer, and irrigation water, the total net mineralization was another source of N, which ranged from 78.0 to 86.5 kg N ha−1 in 2000 and 97.0 to 103.1 kg N ha−1 in 2001. For all the treatments, most of the N was taken up by the crop, which accounted for 31.7–93.2% and 29.3–81.8% of the total N input in 2000 and 2001, respectively.
In our research, the results indicated that gaseous N loss (volatilization and denitrification) was the main pathway of N loss in the open-field system under most treatments. The amount of total gaseous N was 53.5–212.9 kg N ha−1 in 2000 and 37.5–179.2 kg N ha−1 in 2001, accounting for 14.8–23.8% and 9.5–16.3% of the total N input, respectively. The maximum gaseous N release occurred under the W1N1 treatment, reaching up to 212.9 and 179.2 kg N ha−1 over the two years. The gaseous N loss varied with the different fertilization and irrigation practices. Under the same irrigation management, optimizing N fertilizer application could effectively reduce gaseous N loss. For example, under the W1 treatment, the N2 and N3 treatments reduced total gaseous N loss by 56–59% and 48–65% by saving 68–73% and 67–80% of the input of N fertilizer compared to the N1 treatment.
The total amount of nitrate leaching ranged from 12.2 to 195.7 kg N ha−1 in 2000 and 0.5 to 57.0 kg N ha−1 in 2001, respectively. In both years, the W1N1 treatment resulted in the highest total nitrate leaching, accounting for 18.9% and 5.2% of the total N input. In our study, the interaction of both irrigation and fertilizer management had a comparable impact on the total nitrate leaching. For the same fertilizer treatment, reducing irrigation amount significantly decreased the total nitrate leaching in 2000 and 2001. For example, the amount of irrigation under the W2N1 treatment was reduced by 85–97% compared with the W1N1 treatment, which resulted in the amount of total nitrate leaching under the W2N1 treatment being reduced by 166.3 and 55.3 kg N ha−1 in 2000 and 2001, respectively. For the same irrigation treatment, total nitrate leaching under W1N2 was reduced by 41–61% compared with W1N1.
The effects of different irrigation and fertilizer management on the amount of soil N residual (N balance) in the 0.9 m soil profile were varied (Table 4 and Table 5). Compared with the conventional fertilizer treatment, optimizing fertilizer treatment could significantly reduce N accumulation in the 0.9 m soil profile under all the irrigation treatments. Especially in 2000, due to the amount and intensity of rainfall, the total soil N residual under N2 and N3 was negative. For the conventional fertilizer treatment (N1), optimizing irrigation increased N residual in the 0.9 m soil profile in 2000 and 2001. The amount of N residual in each vegetable-growing season varied with different irrigation and fertilizer management. Except for the W2N2 treatment in 2000, the N residual in cauliflower season was positive, and was the highest among the three vegetables. This was related to higher N application in this region. During the amaranth-growing season, the N residual was negative under most conditions, especially under the optimizing N fertilizer treatments.
Our simulation results indicated that the application of N does not exhibit a linear relationship with the N uptake by the vegetables. For example, despite the total N fertilizer application under the W2N2 and W2N3 treatments being lower than the W2N1 treatment in 2000, the N uptake displayed the opposite trend. Moreover, the N uptake of amaranth ranged from 33.2 to 64.4 kg N ha−1 under the optimized fertilizer management (Table 4 and Table 5), which was even higher than the application of N in these treatments. These results also suggested that amaranth growth possibly consumed a part of the residual N from the cauliflower season under the optimized fertilizer treatments. For the spinach season, the residual N patterns were similar to those observed in the amaranth season.
The greatest annual average NUEs reached up to 169.6 kg kg−1 in 2000 and 188.0 kg kg−1 in 2001, which occurred under the W2N3 and W2N2 treatments, respectively (Table 4 and Table 5). The simulation results indicated that for the same irrigation management, the optimizing fertilizer treatments increased annual average NUEs. For example, the annual average NUE under the W1N2 treatment increased by 28% in 2000 and 22% in 2001 compared with the W1N1 treatment. Similarly, the annual average NUE under the W2N2 treatment increased by 29% in 2000 and 20% in 2001 compared with the W2N1 treatment. On the other hand, optimizing irrigation also effectively increased annual average NUEs. For example, the annual average NUE under the W2N1 treatment increased by 29% in 2000 and 22% in 2001, compared with the W1N1 treatment. Moreover, the NUE of each vegetable varied in different years and treatments. In 2000, spinach had the highest NUE among the three vegetables, while in 2001, amaranth did. Over the two years, cauliflower obtained the highest NUE under the W2N2 treatment. For spinach, the W2N3 treatment yielded the maximum NUE. For amaranth, the highest NUE was observed under the W2N3 and W2N2 treatments in 2000 and 2001, respectively.

3. Discussion

3.1. Water Consumption Responses to Different Water and N Management Practices

Unlike greenhouse vegetable fields, open-field vegetable production relies on both rainfall and irrigation as the main sources of water input, while ET is the main process of water loss. Water drainage only occupies a small fraction of the total water consumption [34,35,36]. As Table 2 and Table 3 show, most of the water was drained during the amaranth-growing season in 2000. However, there was almost no water leakage during this time in 2001. This may be because of the heavy rainfall in 2000 (Figure 6).
In our research, the EU-Rotate_N model used the Penman/Monteith method of the FAO to simulate ET and the Ritchie water balance equation to calculate water drainage [24,37,38]. ET in the Penman/Monteith method is mainly determined by climatic factors and vegetable cultivar characteristics. As shown in Table 2 and Table 3, the ET values under different N treatments were similar and showed no significant differences, implying that the water uptake capacity of the vegetable crops was not significantly affected by N input. Additionally, ET was strongly correlated with water input in both years of our study, which is consistent with previous studies [39].
Water drainage is a major process of water consumption in agricultural systems [40]. It is mainly influenced by soil properties and irrigation practices, and it varies significantly among different water treatments. In particular, under W2 in 2001, the water drainage under each N treatment was significantly lower than the other treatments, and the yield reduction in the vegetable crops was minimal. Our results are consistent with those of Cheng et al. [41], who found that water-saving irrigation can effectively reduce water drainage in vegetable fields without significantly impacting yield. Besides, there are some differences between the different N treatments. In 2001, under W2, the water drainage under N3 was slightly higher than the other N treatments, possibly due to the reduced root development caused by N deficiency, which in turn decreased the water uptake capacity of the vegetables [42].
Climate change will intensify water scarcity, the structural and environmental conditions of open-field vegetable production are more complex than those of greenhouse vegetable production [43,44], and therefore, vegetable management must be adapted to save water [45]. It is necessary to apply greenhouse cultivation techniques for vegetable crops that require sufficient water supply, especially in regions with low rainfall. In addition to optimizing irrigation management and applying greenhouse cultivation techniques, farmers can also adopt drought-tolerant vegetable crops such as amaranth or spinach instead of cauliflower especially in arid areas to reduce irrigation costs and the drought risks.

3.2. The Fates of N Responses to Different Water and N Management Practices

In this study, gaseous N loss and nitrate leaching were identified as the main N loss pathways within the open-field vegetable system, with N fertilizer and manure constituting the predominant sources of N input. For all the treatments, crop N uptake was the major N output, which accounted for 31.7–93.2% and 29.3–81.8% of the total N input in 2000 and 2001, respectively. Except for the W2 treatment in 2000, crop N uptake was related positively to N fertilizer rates, which may be caused by greater precipitation in 2000 [7,22].
Apart from manure, fertilizer, and irrigation water, total net mineralization was another source of N. Despite the net mineralization of soil N being significantly lower than the other two sources, it may play a critical role for amaranth under the N2 and N3 treatments due to reduced N input. As indicated in Table 4 and Table 5, the total net mineralization observed in our investigation was slightly lower than the values reported by other researchers [25,46]. The variation in the results of our research could be explained by the heterogeneity of soil properties, manure application, and climatic conditions. The simulation results demonstrated that different irrigation practices influenced net mineralization rates, corroborating earlier research that soil moisture and temperature affect mineralization velocity [47]. However, our simulation results indicated that the soil net mineralization rate remained stable across varying fertilization conditions. Notably, different N applications are known to impact net mineralization [48], suggesting that the EU-Rotate_N model requires refinement in its approach to N net mineralization estimation.
Unlike greenhouse vegetable fields, open-field vegetable production had higher gaseous N losses than nitrate leaching, which was the main form of N loss under all the treatments in this study [23,49]. The amount of gaseous N release was the highest under the W1N1 treatment, up to 212.9 and 179.2 kg N ha−1 in two years, which was significantly higher than that observed in greenhouse vegetable fields [23,50,51]. Previous studies reported that the amount of gaseous N loss had a positive relationship with N fertilizer rates [23,49,50], which was also proven by our results. On the other hand, optimizing irrigation had a minor effect on reducing gaseous N release. Compared with the W1N1 treatment, the amount of irrigation under the W2N1 treatment was reduced by about 28–30%, and the gaseous N loss declined by 18–26%. This may be because reducing irrigation could decrease gaseous N loss and increase N uptake, as reported by Wang et al. [52]. Under optimized irrigation practices, the concentration of air in the soil remained at a higher level, which inhibited denitrifying bacteria activity. Furthermore, more than 50% of the total N fertilizer was applied in the cauliflower-growing season, which resulted in a higher potential for gaseous N loss.
Many previous studies reported that optimizing irrigation and N fertilizer input is effective in reducing nitrate leaching [8,9,10,22], which is consistent with our results. The contribution of nitrate leaching during each vegetable-growing season to the total nitrate leaching varied, more than half of the total nitrate leaching occurred in the amaranth-growing season in 2000. However, most of the nitrate leaching was observed in cauliflower and spinach seasons in 2001. Furthermore, the nitrate leaching in 2000 was significantly higher than that in 2001. Previous research proposed that the amount of nitrate leaching is mainly determined by the amount and intensity of rainfall [7,22]. Compared with 2001, rainfall was greater and more intense in 2000, and most of it occurred during the amaranth season in this study (Figure 6).
Previous research reported that less irrigation and drip irrigation could significantly increase N residual by saving water input [52,53,54]. In this study, although drip irrigation was not employed, water input declined by about 30% under the optimized irrigation treatment (W2). The amount of N residual in each vegetable-growing season varied with different irrigation and fertilizer management. Except for the W2N2 treatment in 2000, the N residual in the cauliflower season was positive, and was the highest among the three vegetables. This was related to the higher N application in this region [55,56]. During the amaranth-growing season, the N residual was negative under most conditions, especially under the optimized fertilizer treatments. Our results indicated that amaranth N uptake ranged from 33.2 to 64.4 kg N ha−1 under the optimized fertilizer management (Table 4 and Table 5) and could surpass the N input in the amaranth season. This means amaranth growth possibly consumed a part of residual N from the cauliflower season under the optimized fertilizer treatments.

3.3. Vegetable Yield, WUE, and NUE Responses to Different Water and N Management Practices

As is known, reduced irrigation and nutrient input mostly decreases the yield in most of agricultural systems; however, it is not applicable to regions with excessive water and nitrogen input [57,58]. Similarly, as shown in Table 2 and Table 3, the total yield of three vegetables generally follows the trend that the yield decreased with the decrease in water and N inputs, but there is no significant variation among the different irrigation and fertilizer treatments in both years. However, there was a large variation in the WUEs and NUEs of different treatments.
Previous studies confirmed that WUE had a negative relationship with irrigation input [10,59,60], which is consistent with our results. Conversely, investigations into field crops indicated that increasing nitrogen application can enhance the WUE in wheat due to enhanced root development [61]. However, this effect was not pronounced in our current study, which may be attributed to the less-developed root systems typically associated with vegetable crops.
The simulation results indicated that for the same irrigation management, optimizing the fertilizer treatments increased the annual average NUEs. Furthermore, the optimization of irrigation practices alone also significantly enhanced the annual average NUEs, corroborating the findings from prior research [52]. Our results were significantly higher than those reported by Wang et al. [7], but similar to the results of other studies [49,51]. In this study, the following vegetables grown could deplete N left by pre-season vegetables under optimized fertilizer management. Therefore, the effects of this fertilizer management on improving NUE were better than the results obtained by only reducing N fertilizer input.
The results of our study indicate that excess N or water does not significantly improve the yield of the vegetables, as it is mainly limited by the genetic potential of vegetables, and the less-developed root system of vegetable crops allows a large amount of water and N to be wasted [7,8]. However, the judicious timing of water and nitrogen applications can markedly improve the efficiency of resource utilization in the vegetable production system [62,63].

4. Materials and Methods

4.1. Site Description

The field experiment was conducted at the China Agricultural University experiment station, located in Dongbeiwang Town, Beijing, China (39.9° N, 116.4° E). The site belongs to a temperate semi-humid continental monsoon climate. The annual average air temperature and precipitation are 11.5 °C and 630 mm. The soil type is Calcaric Cambisol [64]. In the topsoil, soil organic matter is 15 g kg−1, total N is 1.2 g kg−1, and the available phosphorus and available potassium are 21.9 mg kg−1 and 109.6 mg kg−1, respectively.

4.2. Experimental Design and Field Management

Vegetable systems: The field experiment lasted for two years, from April 2000 to October 2001. Based on the climate and physiological characteristics of vegetables, three types of vegetables were cultivated each year, including cauliflower (Brassica oleracea L. var. botrytis cv. Xuefeng), amaranth (Amaranthus tricolor L. Local variety) and spinach (Spinacia oleracea L. cv. Boza 18) (Figure S3). The detailed planting and harvest dates of the three vegetables in both years are shown in Table 6. Cauliflower was planted with 50 cm between plants and 42.5 cm between rows. Amaranth and spinach were both sown in strips with 24.0 cm between rows.
Water and fertilization practices: There were two water treatments (W1, traditional irrigation; W2, optimal irrigation) and three N fertilizer treatments (N1, traditional N fertilizer; N2, optimal N fertilizer; N3, 80% N2) with three replications in this study. A total of six treatments (W1N1, W1N2, W1N3, W2N1, W2N2, and W2N3) were arranged in a randomized block design, with the area of each plot being 12 m × 12 m.
The traditional water (W1) treatment was based on local farmers’ practices. For the W2 treatment, the soil water content was kept at 50%-80% of field capacity. Time Domain Reflectometry (TDR) probes were installed to monitor the soil water content and calculate the irrigation amount needed in W2. The three N treatments were implemented as follows: (1) N1: same as farmers’ practice; (2) N2: the rate of optimal N was calculated by the Expert-N system (see Equation (1)) [65]; (3) based on the N2 treatment, the rate for N fertilizer of spinach was reduced by about 20%.
Nopt = Nupt + Nres + Nloss − Nini − Nhum − Nroot
where Nopt is the optimal N fertilizer dosage; Nupt is the expected amount of N taken up by the crop during the whole growing season; Nres is the amount of mineral N in the root zone after harvest; Nloss is the amount of N loss calculated with the following equation: Nloss = Ninput × 0.019 × growing season (week), Ninput refers to the quantity of nitrogen utilized in farmers’ practices; Nini is the amount of mineral N in the root zone at the beginning of each crop growing season; Nhum is the amount of N from soil organic N mineralization; and Nroot is the amount of N from crop residue mineralization. In this experiment, most of the crop residues were removed from the field, so this part was zero.
In all the treatments, manure was applied as base fertilizer and mixed with tillage into the topsoil 20 cm before cauliflower cultivation each year. The amount of manure was 85 and 120 kg N ha−1 in 2000 and 2001, respectively. Under the W1 and W2 treatments, the methods of irrigation were furrow and micro-sprinkler irrigation systems, respectively. For cauliflower, N fertilizer was topdressed at the end of April and mid-May each year. Amaranth did not receive any topdressing. For spinach, N fertilizer was topdressed in mid-September and early October each year. For each vegetable season, except for amaranth, the rate of fertilizer N application each time was one-third of the total amount. Urea was applied by broadcasting, followed by immediate irrigation. The amount of phosphorus and potassium were the same for all the treatments. Conventional tillage was applied before each vegetable sowing. The details of water and N management practices are shown in Table S1.

4.3. Soil and Crop Sampling and Analyzing

We excavated a soil profile pit to 120 cm, and soil samples were collected from four soil layers: 0–30, 30–60, 60–90, and 90–120 cm before the experiment. Basic soil physical/chemical properties, including soil bulk density, soil saturated water content, field capacity, wilting point, and pH, were measured. The soil bulk density and water content were measured by oven-drying for 24 h at 105 °C [66]. The soil water characteristic curves for estimating the saturated water content, θs; field capacity, θfc; and wilting water content, θwp were determined with a pressure plate apparatus at 0, −33 and −1500 kPa, respectively [67]. Soil pH was determined using a soil/water (1:2.5) slurry. The results are shown in Table 7.
In each plot, the soil volumetric water content from 0 to 120 cm at 30 cm intervals in the soil profile was measured by TDR (MP-917) every two days. The experiment initially aimed to investigate nitrate leaching at a 90 cm depth under various water and fertilizer regimes. Therefore, the suction cups were installed at only 90 cm soil depth. Soil solutions were periodically collected at intervals ranging from 7 to 15 days, and were analyzed for mineral N to ascertain soil nitrate concentrations with the continuous flow analyzer (TRAACS 2000). Plant biomass, fresh yield, and nitrogen content were quantified on a weekly basis through the collection of plant samples. The yields of vegetables were determined within a 1 m2 quadrat, replicated four times within each experimental plot.

4.4. Model Description and Input

In this study, the EU-Rotate_N model was utilized to simulate soil water movement, nitrogen transport, and vegetable growth. It contains modules for simulating various soil/plant processes including root development, N mineralization from soil organic matter and crop residues, water and N uptake by roots, and water movement and N transport in soil. The organization of the main model modules is shown in Figure S4. In the EU-Rotate_N model, the vegetable growth algorithm was derived from the N_ABLE model [68,69]. Water demand for crop growth was calculated by the Penman/Monteith method of the FAO [37], while soil water movement was simulated following Ritchie’s soil water balance approach [70]. Surface runoff estimation adhered to the methodology advocated by the U.S. National Resource Conservation Service [71]. The N mineralization simulation was based on the DAISY model’s algorithm [72]. Ammonia emissions from manure were quantified using the empirical formula from the ALFAM model [73]. The simulations of urea hydrolysis and gaseous N emissions were based on the AMOVOL model [74]. A detailed model description is available in the literature [24,68].
The model inputs include site location (latitude and altitude), basic soil physico-chemical properties (e.g., bulk density and hydraulic properties), crop information (species, row spacing, planting density, and sowing and harvesting schedules), field management practices, soil initial water and mineral nitrogen condition, and local meteorological data [22,24,28].

4.5. Model Evaluation Statistics and Data Analysis

Three statistical indices were used to evaluate the model performance: root mean square error (RMSE), Nash/Sutcliffe efficiency coefficient (E), and agreement index (d) [32,75].
R M S E = i = 1 n ( P i O i ) 2 n
E = 1 i = 1 n ( O i P i ) 2 i = 1 n ( O i O ) 2
d = 1 i = 1 n ( O i P i ) 2 i = 1 n ( | P i O | + | O i O | ) 2
where n is the number of samples, Pi and Oi are the simulated and observed values, and O is the mean of the measured data. The closer the value of root mean square error (RMSE) is to 0, the more accurate the model. Modeling efficiency (E) ranges from −∞ to 1. The agreement index (d) represents the ratio of the mean square error and the potential error. When d = 1, it indicates a perfect simulation of results. Based on van Liew and Garbrecht [31], an acceptable simulation should have E > 0.36 and d > 0.7.
The analysis of variance (ANOVA)was used to compare the effects of the different water and N management on vegetable yields using SPSS 23 (International Business Machines Corporation, Armonk, NY, USA). The graphs were generated using Origin 2018 (OriginLab Corporation, Northampton, MA, USA).

5. Conclusions

In this study, we used the soil water content, nitrate concentration, crop dry matter, and other data collected from an open-field vegetable field with different irrigation and N fertilizer managements to validate the EU-Rotate_N model. Statistical analysis indicated that the EU-Rotate_N model could simulate the dynamics of water and N in an open-field intensive vegetable system in the North China Plain.
The results showed the yields of amaranth and spinach had no significant differences among all the treatments in 2000 and 2001. However, the cauliflower yield under the W1N2 and W1N3 treatments obviously reduced in 2001, indicating that excess fertilizer N or water does not significantly improve the yield of vegetables.
Unlike greenhouse vegetable fields, open-field vegetable production had higher gaseous N losses than nitrate leaching, which was the main form of N loss under all the treatments in this study. Gaseous N loss accounted for 16.5% of the total N input, while the annual amount of nitrate leaching was less and accounted for only 8.4% of the total N input. Optimizing irrigation and N fertilizer practices could effectively reduce gaseous N release and nitrate leaching, and increase NUE. Compared with W1, gaseous N loss under W2 was reduced by 18–26% and the annual average NUEs increased by 22–29%. Compared to the N1 treatments, gaseous N under the N2 and N3 treatments reduced by 48–72% and increased NUE by 17–37%. Based on two years of field experiments, the simulation results showed that the W2N3 treatment was the best water and N fertilizer management, and should be recommended to local farmers.
Open-field vegetable production is significantly influenced by climatic conditions. In 2000, the rainfall reached 151 mm during the period of amaranth, which led to the highest drainage and nitrate leaching. In contrast, it was only 69.1 mm in 2001, resulting in no drainage or nitrate leaching. Climate change will intensify water scarcity, making it necessary to apply greenhouse cultivation techniques for vegetable crops that require sufficient water supply, especially in regions with low rainfall. Additionally, farmers can adopt drought-tolerant vegetable crops such as amaranth or spinach instead of cauliflower especially in arid areas to reduce irrigation costs and drought risks. Furthermore, rationally adjusting the cropping structure can allow subsequent crops to fully utilize the residual N in the soil, thereby reducing the risk of environmental pollution.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13152150/s1, Table S1. Application of irrigation and N fertilizer for field vegetable cultivation in 2000 and 2001. Table S2. Crop and N transformation parameters used in EU-Rotate_N model. Figure S1. Comparisons of simulated and measured crop dry matter weight for three vegetables under different treatments in 2000. Figure S2. Comparisons of simulated and measured crop dry matter weight for three vegetables under different treatments in 2001. Figure S3. Some photos of the experiment site. Figure S4. The organization of the main modules of EU-Rotate_N model.

Author Contributions

Formal analysis, investigation, and writing—original draft, Y.S.; formal analysis and writing—review and editing, S.C.; writing—review and editing, P.F.; data curation, resources, Q.C.; conceptualization, methodology, writing—review and editing, supervision, and funding acquisition, K.H. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the National Key Research and Development program of China (No. 2021YFD1700900) and the 2115 Talent Development Program of China Agricultural University (No. 1191–00109011).

Data Availability Statement

The data is contained within the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wang, X.; Dou, Z.; Shi, X.; Zou, C.; Liu, D.; Wang, Z.; Guan, X.; Sun, Y.; Wu, G.; Zhang, B. Innovative management programme reduces environmental impacts in Chinese vegetable production. Nat. Food 2021, 2, 47–53. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, X.; Davidson, E.A.; Mauzerall, D.L.; Searchinger, T.D.; Dumas, P.; Shen, Y. Managing nitrogen for sustainable development. Nature 2015, 528, 51–59. [Google Scholar] [CrossRef]
  3. Granados, M.; Thompson, R.; Fernández, M.; Martínez-Gaitán, C.; Gallardo, M. Prescriptive–corrective nitrogen and irrigation management of fertigated and drip-irrigated vegetable crops using modeling and monitoring approaches. Agric. Water Manag. 2013, 119, 121–134. [Google Scholar] [CrossRef]
  4. Cui, B.J.; Niu, W.Q.; Du, Y.D.; Zhang, Q. Response of yield and nitrogen use efficiency to aerated irrigation and N application rate in greenhouse cucumber. Sci. Hortic. 2020, 265, 109220. [Google Scholar] [CrossRef]
  5. Li, Y.; Xue, X.; Guo, W.; Wang, L.; Duan, M.; Chen, H.; Chen, F. Soil moisture and nitrate-nitrogen dynamics and economic yield in the greenhouse cultivation of tomato and cucumber under negative pressure irrigation in the North China Plain. Sci. Rep. 2019, 9, 4439. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, H.; Yuan, B.; Hu, X.; Yin, C.; Tang, X. Cucumber production and the economic revenues under various nitrogen applications in an unheated solar greenhouse on the North China Plain. Agron. J. 2021, 113, 3444–3459. [Google Scholar] [CrossRef]
  7. Wang, D.; Guo, L.; Zheng, L.; Zhang, Y.; Yang, R.; Li, M.; Ma, F.; Zhang, X.; Li, Y. Effects of nitrogen fertilizer and water management practices on nitrogen leaching from a typical open field used for vegetable planting in northern China. Agric. Water Manag. 2019, 213, 913–921. [Google Scholar] [CrossRef]
  8. Zhou, W.; Lv, H.; Chen, F.; Wang, Q.; Li, J.; Chen, Q.; Liang, B. Optimizing nitrogen management reduces mineral nitrogen leaching loss mainly by decreasing water leakage in vegetable fields under plastic-shed greenhouse. Environ. Pollut. 2022, 308, 119616. [Google Scholar] [CrossRef] [PubMed]
  9. Gu, J.; Wu, Y.; Tian, Z.; Xu, H. Nitrogen use efficiency, crop water productivity and nitrous oxide emissions from Chinese greenhouse vegetables: A meta-analysis. Sci. Total Environ. 2020, 743, 140696. [Google Scholar] [CrossRef]
  10. Fan, Y.; Hao, X.; Ding, R.; Kang, S. Soil water and nitrogen dynamics from interaction of irrigation and fertilization management practices in a greenhouse vegetable rotation. Soil Sci. Soc. Am. J. 2020, 84, 901–913. [Google Scholar] [CrossRef]
  11. Thompson, R.B.; Tremblay, N.; Fink, M.; Gallardo, M.; Padilla, F.M. Tools and strategies for sustainable nitrogen fertilisation of vegetable crops. In Advances in Research on Fertilization Management of Vegetable Crops; Springer: Cham, Switzerland, 2017; pp. 11–63. [Google Scholar]
  12. Jelle, V. Model-Based Approaches for Fertilizer Recommendations and Their Environmental Impact Due to Nutrient Leaching under Variable Weather Conditions in Open Field Horticultural Production. Ph.D. Thesis, KU Leuven, Leuven, Belgium, 2018. [Google Scholar]
  13. Chen, Q.; Zhang, H.; Li, X.; Christie, P.; Horlacher, D.; Liebig, H.P. Use of a modified N-expert system for vegetable production in the Beijing region. J. Plant Nutr. 2005, 28, 475–487. [Google Scholar] [CrossRef]
  14. Zhai, Z.; Martínez, J.F.; Beltran, V.; Martínez, N.L. Decision support systems for agriculture 4.0: Survey and challenges. Comput. Electron. Agric. 2020, 170, 105256. [Google Scholar] [CrossRef]
  15. Sha, S.; Zhao, X.; Li, Y.; Li, C.; Zhu, L.; Wang, Y.; Gao, Q. Nutrient expert system optimizes fertilizer management to improve potato productivity and tuber quality. J. Sci. Food Agric. 2022, 102, 1233–1244. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, X.; He, P.; Pampolino, M.F.; Li, Y.; Liu, S.; Xie, J.; Hou, Y.; Zhou, W. Narrowing yield gaps and increasing nutrient use efficiencies using the Nutrient Expert system for maize in Northeast China. Field Crops Res. 2016, 194, 75–82. [Google Scholar] [CrossRef]
  17. Zhang, J.; Ding, W.; He, P.; Xu, X.; Abbas, T.; Ullah, S.; Ai, C.; Li, M.; Cui, R.; Jin, C. Establishment and validation of nutrient expert system for radish fertilization management in China. Agron. J. 2019, 111, 2435–2444. [Google Scholar] [CrossRef]
  18. Delgado, J.A.; Ascough II, J.C.; Lighthart, N.; Neer, D. Potential use of a new nitrogen trading tool to assess nitrogen management practices to protect groundwater quality. Comput. Electron. Agric. 2020, 169, 105195. [Google Scholar] [CrossRef]
  19. Li, Z.; Wen, X.; Hu, C.; Li, X.; Li, S.; Zhang, X.; Hu, B. Regional simulation of nitrate leaching potential from winter wheat-summer maize rotation croplands on the North China Plain using the NLEAP-GIS model. Agric. Ecosyst. Environ. 2020, 294, 106861. [Google Scholar] [CrossRef]
  20. Shaffer, M.; Halvorson, A.; Pierce, F. Nitrate leaching and economic analysis package (NLEAP): Model description and application. In Managing Nitrogen for Groundwater Quality and Farm Profitability; Soil Science Society of America: Madison, WI, USA, 1991; pp. 285–322. [Google Scholar]
  21. Nendel, C. Evaluation of Best Management Practices for N fertilisation in regional field vegetable production with a small-scale simulation model. Eur. J. Agron. 2009, 30, 110–118. [Google Scholar] [CrossRef]
  22. Hua, B.; Cao, Z.; Zhang, K.; Xu, X.; Zhang, Y.; Dai, H.; Zhang, Z.; Jiang, J.; Miao, M. Simulation of greenhouse cucumber growth, water and nitrogen dynamics in areas with high groundwater (HG) levels using the HG EU-Rotate_N model. Veg. Res. 2022, 2, 1–11. [Google Scholar] [CrossRef]
  23. Sun, Y.; Zhang, J.; Wang, H.; Wang, L.; Li, H. Identifying optimal water and nitrogen inputs for high efficiency and low environment impacts of a greenhouse summer cucumber with a model method. Agric. Water Manag. 2019, 212, 23–34. [Google Scholar] [CrossRef]
  24. Ullah, I.; Cao, Z.; Bing, H.; Xu, X.; Miao, M. Simulation of Soil Water and Nitrogen Dynamics for Tomato Crop Using EU-Rotate_N Model at Different Nitrogen Levels in the Greenhouse. Agronomy 2023, 13, 2006. [Google Scholar] [CrossRef]
  25. Xu, X.; Wang, C.; Wang, H.; Zhang, Y.; Cao, Z.; Zhang, Z.; Dai, H.; Miao, M. Development and performance evaluation of an APP for vegetable fertilization and irrigation management originated from EU-Rotate_N. Agric. Water Manag. 2023, 289, 108520. [Google Scholar] [CrossRef]
  26. Ji, L.; Ting, X.; Han, H.; Guoying, D.; Yufei, L.; Xi, W.; Hefa, Y.; Shengnan, L. Ecological Mechanisms and Regulation of Greenhouse Organic Vegetable Production in North China. In Agroecology in China: Science, Practice, and Sustainable Management; CRC Press: Boca Raton, FL, USA, 2017; pp. 259–288. [Google Scholar]
  27. Li, J.; Liu, H.; Wang, H.; Luo, J.; Zhang, X.; Liu, Z.; Zhang, Y.; Zhai, L.; Lei, Q.; Ren, T. Managing irrigation and fertilization for the sustainable cultivation of greenhouse vegetables. Agric. Water Manag. 2018, 210, 354–363. [Google Scholar] [CrossRef]
  28. Øvsthus, I.; Thorup-Kristensen, K.; Seljåsen, R.; Riley, H.; Dörsch, P.; Breland, T.A. Calibration of the EU-Rotate_N model with measured C and N mineralization from potential fertilizers and evaluation of its prediction of crop and soil data from a vegetable field trial. Eur. J. Agron. 2021, 129, 126336. [Google Scholar] [CrossRef]
  29. Suárez-Rey, E.; Gallardo, M.; Romero-Gámez, M.; Giménez, C.; Rueda, F.J. Sensitivity and uncertainty analysis in agro-hydrological modelling of drip fertigated lettuce crops under Mediterranean conditions. Comput. Electron. Agric. 2019, 162, 630–650. [Google Scholar] [CrossRef]
  30. Meng, F.; Hu, K.; Feng, P.; Feng, G.; Gao, Q. Simulating the Effects of Different Textural Soils and N Management on Maize Yield, N Fates, and Water and N Use Efficiencies in Northeast China. Plants 2022, 11, 3338. [Google Scholar] [CrossRef]
  31. Han, L.; Li, Y.; Hou, Y.; Liang, H.; Feng, P.; Hu, K. Improvement of water and nitrogen use efficiencies by alternative cropping systems based on a model approach. Plants 2023, 12, 597. [Google Scholar] [CrossRef] [PubMed]
  32. Yang, J.; Yang, J.-Y.; Liu, S.; Hoogenboom, G. An evaluation of the statistical methods for testing the performance of crop models with observed data. Agric. Syst. 2014, 127, 81–89. [Google Scholar] [CrossRef]
  33. Suárez-Rey, E.; Romero-Gámez, M.; Giménez, C.; Thompson, R.; Gallardo, M. Use of EU-Rotate_N and CropSyst models to predict yield, growth and water and N dynamics of fertigated leafy vegetables in a Mediterranean climate and to determine N fertilizer requirements. Agric. Syst. 2016, 149, 150–164. [Google Scholar] [CrossRef]
  34. Arunadevi, K.; Ashok, A.; Singh, M. Estimation of evapotranspiration under polyhouse and open field conditions in capsicum. Indian J. Hortic. 2022, 79, 208–214. [Google Scholar] [CrossRef]
  35. Gong, X.; Qiu, R.; Ge, J.; Bo, G.; Ping, Y.; Xin, Q.; Wang, S. Evapotranspiration partitioning of greenhouse grown tomato using a modified Priestley–Taylor model. Agric. Water Manag. 2021, 247, 106709. [Google Scholar] [CrossRef]
  36. Wen, J.; Ning, S.; Wei, X.; Guo, W.; Sun, W.; Zhang, T.; Wang, L. The Impact of Insect-Proof Screen on Microclimate, Reference Evapotranspiration and Growth of Chinese Flowering Cabbage in Arid and Semi-Arid Region. Horticulturae 2022, 8, 704. [Google Scholar] [CrossRef]
  37. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; Fao: Rome, Italy, 1998; 300p. [Google Scholar]
  38. Flores-Velazquez, J.; Akrami, M.; Villagrán, E. The Role of Radiation in the Modelling of Crop Evapotranspiration from Open Field to Indoor Crops. Agronomy 2022, 12, 2593. [Google Scholar] [CrossRef]
  39. Shah, H.L.; Zhou, T.; Sun, N.; Huang, M.; Mishra, V. Roles of irrigation and reservoir operations in modulating terrestrial water and energy budgets in the Indian Subcontinental River basins. J. Geophys. Res. Atmos. 2019, 124, 12915–12936. [Google Scholar] [CrossRef]
  40. Liang, H.; Lv, H.; Batchelor, W.D.; Lian, X.; Wang, Z.; Lin, S.; Hu, K. Simulating nitrate and DON leaching to optimize water and N management practices for greenhouse vegetable production systems. Agric. Water Manag. 2020, 241, 106377. [Google Scholar] [CrossRef]
  41. Cheng, M.; Wang, H.; Fan, J.; Zhang, S.; Liao, Z.; Zhang, F.; Wang, Y. A global meta-analysis of yield and water use efficiency of crops, vegetables and fruits under full, deficit and alternate partial root-zone irrigation. Agric. Water Manag. 2021, 248, 106771. [Google Scholar] [CrossRef]
  42. Nemeskéri, E.; Helyes, L. Physiological responses of selected vegetable crop species to water stress. Agronomy 2019, 9, 447. [Google Scholar] [CrossRef]
  43. Rahimikhoob, H.; Sohrabi, T.; Delshad, M. Assessment of reference evapotranspiration estimation methods in controlled greenhouse conditions. Irrig. Sci. 2020, 38, 389–400. [Google Scholar] [CrossRef]
  44. Yan, H.; Zhao, S.; Zhang, C.; Zhang, J.; Wang, G.; Li, M.; Deng, S.; Liang, S.; Jiang, J. Calibration and assessment of evapotranspiration methods for cucumber plants in a Venlo-type greenhouse. Irrig. Drain. 2024, 73, 119–135. [Google Scholar] [CrossRef]
  45. Mirás-Avalos, J.M.; Rubio-Asensio, J.S.; Ramírez-Cuesta, J.M.; Maestre-Valero, J.F.; Intrigliolo, D.S. Irrigation-Advisor—A decision support system for irrigation of vegetable crops. Water 2019, 11, 2245. [Google Scholar] [CrossRef]
  46. Tei, F.; De Neve, S.; de Haan, J.; Kristensen, H.L. Nitrogen management of vegetable crops. Agric. Water Manag. 2020, 240, 106316. [Google Scholar] [CrossRef]
  47. Dettmann, U.; Kraft, N.N.; Rech, R.; Heidkamp, A.; Tiemeyer, B. Analysis of peat soil organic carbon, total nitrogen, soil water content and basal respiration: Is there a ‘best’drying temperature? Geoderma 2021, 403, 115231. [Google Scholar] [CrossRef]
  48. Zhang, X.; Zhu, B.; Yu, F.; Cheng, W. Plant inputs mediate the linkage between soil carbon and net nitrogen mineralization. Sci. Total Environ. 2021, 790, 148208. [Google Scholar] [CrossRef] [PubMed]
  49. Fan, D.; Wang, X.; Song, D.; Shi, Y.; Chen, Y.; Wang, J.; Cao, B.; Zou, G.; He, W. Optimizing nitrogen management to mitigate gaseous losses and improve net benefits of an open-field Chinese cabbage system. J. Environ. Manag. 2022, 318, 115583. [Google Scholar] [CrossRef] [PubMed]
  50. Qasim, W.; Xia, L.; Lin, S.; Wan, L.; Zhao, Y.; Butterbach-Bahl, K. Global greenhouse vegetable production systems are hotspots of soil N2O emissions and nitrogen leaching: A meta-analysis. Environ. Pollut. 2021, 272, 116372. [Google Scholar] [CrossRef] [PubMed]
  51. Ti, C.; Luo, Y.; Yan, X. Characteristics of nitrogen balance in open-air and greenhouse vegetable cropping systems of China. Environ. Sci. Pollut. Res. 2015, 22, 18508–18518. [Google Scholar] [CrossRef]
  52. Wang, S.; Xia, G.; Zheng, J.; Wang, Y.; Chen, T.; Chi, D.; Bolan, N.S.; Chang, S.X.; Wang, T.; Ok, Y.S. Mulched drip irrigation and biochar application reduce gaseous nitrogen emissions, but increase nitrogen uptake and peanut yield. Sci. Total Environ. 2022, 830, 154753. [Google Scholar] [CrossRef] [PubMed]
  53. Lv, H.; Lin, S.; Wang, Y.; Lian, X.; Zhao, Y.; Li, Y.; Du, J.; Wang, Z.; Wang, J.; Butterbach-Bahl, K. Drip fertigation significantly reduces nitrogen leaching in solar greenhouse vegetable production system. Environ. Pollut. 2019, 245, 694–701. [Google Scholar] [CrossRef]
  54. Ma, G.; Liu, W.; Li, S.; Zhang, P.; Wang, C.; Lu, H. Determining the optimal N input to improve grain yield and quality in winter wheat with reduced apparent N loss in the North China Plain. Front. Plant Sci. 2019, 10, 371354. [Google Scholar] [CrossRef]
  55. Zhang, Y.; Wang, H.; Lei, Q.; Luo, J.; Lindsey, S.; Zhang, J.; Zhai, L.; Wu, S.; Zhang, J.; Liu, X. Optimizing the nitrogen application rate for maize and wheat based on yield and environment on the Northern China Plain. Sci. Total Environ. 2018, 618, 1173–1183. [Google Scholar] [CrossRef]
  56. Qiu, H.; Yang, S.; Jiang, Z.; Xu, Y.; Jiao, X. Effect of irrigation and fertilizer management on rice yield and nitrogen loss: A meta-analysis. Plants 2022, 11, 1690. [Google Scholar] [CrossRef] [PubMed]
  57. Tolessa, E.S. A review on water and nitrogen use efficiency of potato (Solanum tuberosum L.) in relation to its yield and yield components. Arch. Agric. Environ. Sci. 2019, 4, 119–132. [Google Scholar] [CrossRef]
  58. Pabuayon, I.L.B.; Singh, S.; Lewis, K.L.; Ritchie, G.L. Water extraction and productivity of cotton, sorghum, and sesame under deficit irrigation. Crop Sci. 2019, 59, 1692–1700. [Google Scholar] [CrossRef]
  59. Zain, M.; Si, Z.; Li, S.; Gao, Y.; Mehmood, F.; Rahman, S.-U.; Mounkaila Hamani, A.K.; Duan, A. The coupled effects of irrigation scheduling and nitrogen fertilization mode on growth, yield and water use efficiency in drip-irrigated winter wheat. Sustainability 2021, 13, 2742. [Google Scholar] [CrossRef]
  60. Li, Y.; Huang, G.; Chen, Z.; Xiong, Y.; Huang, Q.; Xu, X.; Huo, Z. Effects of irrigation and fertilization on grain yield, water and nitrogen dynamics and their use efficiency of spring wheat farmland in an arid agricultural watershed of Northwest China. Agric. Water Manag. 2022, 260, 107277. [Google Scholar] [CrossRef]
  61. Chang, Y.-S.; Lin, K.-H.; Chen, Y.-J.; Wu, C.-W.; Chang, Y.-J. Saving water used for vegetable production by applying regulated deficit irrigation practices. Acta Sci. Polonorum. Hortorum Cultus 2021, 20, 27–36. [Google Scholar] [CrossRef]
  62. Sarkar, S.; Saha, R.; Patra, S.K. Irrigation and fertilizer management on growth, yield, water and fertilizer use efficiencies on cabbage in a sandy loam soil. Int. J. Curr. Microbiol. Appl. Sci. 2020, 9, 64–77. [Google Scholar] [CrossRef]
  63. Ullah, H.; Santiago-Arenas, R.; Ferdous, Z.; Attia, A.; Datta, A. Improving water use efficiency, nitrogen use efficiency, and radiation use efficiency in field crops under drought stress: A review. Adv. Agron. 2019, 156, 109–157. [Google Scholar]
  64. FAO, F. Unesco Soil Map of the World, Revised Legend, with corrections and updates. World Soil Resour. Rep. 1988, 60, 140. [Google Scholar]
  65. Fink, M.; Scharpf, H. N-Expert-A decision support system for vegetable fertilization in the field. In Proceedings of the Workshop on Ecological Aspects of Vegetable Fertilization in Integrated Crop Production in the Field, Einsiedeln, Switzerland, 7–11 September 1992; pp. 67–74. [Google Scholar]
  66. Grossman, R.; Reinsch, T. The solid phase, bulk density and linear extensibility. In Methods of Soil Analysis. Part 4; Soil Science Society of America: Madison, WI, USA, 2002; pp. 201–228. [Google Scholar]
  67. Givi, J.; Prasher, S.; Patel, R. Evaluation of pedotransfer functions in predicting the soil water contents at field capacity and wilting point. Agric. Water Manag. 2004, 70, 83–96. [Google Scholar] [CrossRef]
  68. Rahn, C.; Zhang, K.; Lillywhite, R.; Ramos, C.; Doltra, J.; De Paz, J.; Riley, H.; Fink, M.; Nendel, C.; Thorup Kristensen, K. EU-Rotate_N–a decision support system–to predict environmental and economic consequences of the management of nitrogen fertiliser in crop rotations. Eur. J. Hortic. Sci. 2010, 75, 20–32. [Google Scholar]
  69. Greenwood, D.; Rahn, C.; Draycott, A.; Vaidyanathan, L.; Paterson, C. Modelling and measurement of the effects of fertilizer-N and crop residue incorporation on N-dynamics in vegetable cropping. Soil Use Manag. 1996, 12, 13–24. [Google Scholar] [CrossRef]
  70. Ritchie, J.T. Soil water balance and plant water stress. In Understanding Options for Agricultural Production; Springer: Dordrecht, The Netherlands, 1998; pp. 41–54. [Google Scholar]
  71. NRCS; USDA. Estimation of direct runoff from storm rainfall. In National Engineering Handbook; U.S. Department of Agriculture (USDA): Washington, DC, USA, 2004; Part 630. [Google Scholar]
  72. Hansen, S.; Jensen, H.; Nielsen, N.; Svendsen, H. Daisy, a Soil Plant System Model. Danish Simulation Model for Transformation and Transport of Energy and Matter in the Soil Plant Atmosphere System; The National Agency for Environmental Protection: Copenhagen, Denmark, 1990; 369p. [Google Scholar]
  73. Bussink, D.; Nicholson, F. Ammonia volatilization from field-applied animal slurry—The ALFAM model. Atmos. Environ. 2002, 36, 3309–3319. [Google Scholar]
  74. Sadeghi, A.; Kissel, D.; Cabrera, M. Estimating molecular diffusion coefficients of urea in unsaturated soil. Soil Sci. Soc. Am. J. 1989, 53, 15–18. [Google Scholar] [CrossRef]
  75. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
Figure 1. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W1N1 and W1N2 treatments in 2000 and 2001.
Figure 1. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W1N1 and W1N2 treatments in 2000 and 2001.
Plants 13 02150 g001
Figure 2. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W1N3 and W2N1 treatments in 2000 and 2001.
Figure 2. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W1N3 and W2N1 treatments in 2000 and 2001.
Plants 13 02150 g002
Figure 3. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W2N2 and W2N3 treatments in 2000 and 2001.
Figure 3. Comparisons of simulated and measured soil water content in 0–120 cm soil profile for W2N2 and W2N3 treatments in 2000 and 2001.
Plants 13 02150 g003
Figure 4. Comparisons of simulated and measured soil nitrate concentration at 60–90 cm soil layer under different treatments in 2000 and 2001.
Figure 4. Comparisons of simulated and measured soil nitrate concentration at 60–90 cm soil layer under different treatments in 2000 and 2001.
Plants 13 02150 g004
Figure 5. Relationships of simulated and measured crop N uptake for three vegetables (cauliflower, amaranth, and spinach) under all treatments in 2000 and 2001.
Figure 5. Relationships of simulated and measured crop N uptake for three vegetables (cauliflower, amaranth, and spinach) under all treatments in 2000 and 2001.
Plants 13 02150 g005
Figure 6. Dynamics of water drainage and nitrate leaching at the 90 cm depth under all of the water and fertilizer treatments in 2000 and 2001.
Figure 6. Dynamics of water drainage and nitrate leaching at the 90 cm depth under all of the water and fertilizer treatments in 2000 and 2001.
Plants 13 02150 g006
Table 1. Statistical indices for the predicted soil water content and nitrate concentration at different depths for the validation treatments (W1N2, W1N3, W2N1, W2N2, and W2N3).
Table 1. Statistical indices for the predicted soil water content and nitrate concentration at different depths for the validation treatments (W1N2, W1N3, W2N1, W2N2, and W2N3).
Soil Layer (cm)RMSEEd
Soil water content0–150.040.380.83
15–300.020.470.87
30–600.010.460.88
60–900.010.860.97
90–1200.010.920.98
Soil NO3-N concentration60–9013.410.870.97
Note: RMSE, root mean square error; E, Nash/Sutcliffe efficiency coefficient; d, agreement index.
Table 2. Yield and soil water balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2000.
Table 2. Yield and soil water balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2000.
TreatmentsVegetableRF
mm
I
mm
ET
mm
D
mm
Y
kg ha−1
WUE
kg m−3
IWUE
kg m−3
Wbalance
mm
W1N1cauliflower45.1294.1327.820.423,400 a7.18.0−9.0
amaranth151.085.5104.0111.514,300 a13.816.721.0
spinach65.6128.2130.234.634,800 a26.727.129.0
total261.7507.8562.0166.572,50012.914.341.0
W1N2cauliflower45.1294.1323.120.719,600 ab6.16.7−4.6
amaranth151.085.599.4108.713,700 a13.816.028.4
spinach65.6128.2128.537.132,900 a25.625.728.2
total261.7507.8551.0166.566,20012.013.052.0
W1N3cauliflower45.1294.1319.219.520,500 ab6.47.00.5
amaranth151.085.595.8108.610,800 a11.312.632.1
spinach65.6128.2127.740.236,400 a28.528.425.9
total261.7507.8542.7168.367,70012.513.358.5
W2N1cauliflower45.1208.5281.22.921,400 ab7.610.3−30.5
amaranth151.068.3100.935.813,100 a13.019.282.6
spinach65.689.3114.84.033,100 a28.837.136.1
total261.7366.1496.942.767,60013.618.588.2
W2N2cauliflower45.1208.5273.31.219,700 ab7.29.4−20.9
amaranth151.068.3100.438.015,200 a15.122.380.9
spinach65.689.3115.34.132,500 a28.236.435.5
annual261.7366.1489.043.367,40013.818.495.5
W2N3cauliflower45.1208.5284.80.118,100 bc6.48.9−31.3
amaranth151.068.3101.130.716,500 a16.324.287.5
spinach65.689.3115.46.833,900 a29.438.032.7
total261.7366.1501.337.668,50013.718.788.9
Note: RF, rainfall. I, irrigation. ET, evapotranspiration. D, water drainage. Y, fresh yield. WUE (water use efficiency) = Y/ET. IWUE (irrigation water use efficiency) = Y/I. Wbalance = RF + I − ET − D. Different letters indicate significant differences in yield for same vegetable under different treatments at p < 0.05.
Table 3. Yield and soil water balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2001.
Table 3. Yield and soil water balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2001.
TreatmentsVegetableRF
mm
I
mm
ET
mm
D
mm
Y
kg ha−1
WUE
kg m−3
IWUE
kg m−3
Wbalance
mm
W1N1cauliflower21.9290.7310.530.519,200 a6.26.6−28.4
amaranth69.1107.4127.70.022,500 a17.620.948.8
spinach52.7123.4121.324.430,100 a24.824.430.4
total143.7521.5559.554.971,80012.813.850.8
W1N2cauliflower21.9290.7304.331.016,200 b5.35.6−22.7
amaranth69.1107.4129.60.018,600 a14.417.346.9
spinach52.7123.4118.016.830,200 a25.624.541.3
total143.7521.5551.947.865,00011.812.565.5
W1N3cauliflower21.9290.7297.129.514,500 c4.95.0−14.0
amaranth69.1107.4128.10.018,300 a14.317.048.4
spinach52.7123.4118.026.931,400 a26.615.431.2
total143.7521.5543.256.464,20011.812.365.6
W2N1cauliflower21.9207.6258.00.919,300 a7.59.3−29.4
amaranth69.174.0112.70.022,400 a19.930.330.4
spinach52.784.6117.30.431,600 a26.937.419.6
total143.7366.2488.01.373,30015.020.020.6
W2N2cauliflower21.9207.6257.60.219,000 a7.49.2−28.3
amaranth69.174.0114.00.021,200 a18.628.629.1
spinach52.784.6118.91.725,900 a21.830.616.7
annual143.7366.2490.51.966,10013.518.117.5
W2N3cauliflower21.9207.6254.90.618,600 a7.39.0−26.0
amaranth69.174.0113.70.015,200 a13.420.529.4
spinach52.784.6117.69.530,400 a25.935.910.2
total143.7366.2486.210.164,20013.217.513.6
Note: RF, rainfall. I, irrigation. ET, evapotranspiration. D, water drainage. Y, fresh yield. WUE (water use efficiency) = Y/ET. IWUE (irrigation water use efficiency) = Y/I. Wbalance = RF + I − ET − D. Different letters indicate significant differences in yield for same vegetable under different treatments at p < 0.05.
Table 4. N balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2000.
Table 4. N balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2000.
Treatments ManuFertNirrNnetNupNgasNleaNUE
kg kg−1
Nbalance
kg N ha−1
kg N ha−1
W1N1cauliflower75450.011.525.0162.5141.37.875.1249.9
amaranth0100.03.323.546.924.9129.271.1−74.2
spinach0309.05.032.2120.046.758.7154.4120.8
total75859.019.880.7329.4212.9195.798.2296.5
W1N2cauliflower75165.911.524.5158.665.86.984.745.6
amaranth026.03.321.644.112.779.0100.9−84.9
spinach081.45.032.5113.716.129.4206.7−40.3
total75273.319.878.6316.494.6115.3125.8−79.6
W1N3cauliflower75189.811.524.6163.679.26.382.351.8
amaranth026.23.321.433.214.688.979.0−85.8
spinach072.05.032.0117.716.032.4219.1−57.1
total75288.019.878.0314.5109.8127.6122.7−91.1
W2N1cauliflower75450.08.124.0159.8116.40.877.3280.1
amaranth0100.02.728.744.920.122.0150.644.4
spinach0309.03.532.7122.938.56.6197.0177.2
total75859.014.385.4327.6175.029.4127.1501.7
W2N2cauliflower7583.68.124.1165.332.00.399.7−6.8
amaranth025.72.728.349.615.215.5189.3−23.6
spinach083.93.534.1118.312.23.3242.9−12.3
total75193.214.386.5333.259.419.1163.7−42.7
W2N3cauliflower75101.78.123.8164.736.40.090.07.5
amaranth031.52.728.257.19.08.2222.1−11.9
spinach055.43.532.7116.38.14.0264.0−36.8
total75188.614.384.7338.153.512.2169.6−41.2
Note: Manu, manure; Fert, fertilizer; Nnet, net mineralization; Nirr, N input by irrigation; Nup, crop N uptake. Nlea, nitrate leaching. Ngas, gaseous N emission. NUE (nitrogen use efficiency) = Fresh yield/(Nup + Nlea + Ngas). Nbalance = N input − N output.
Table 5. N balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2001.
Table 5. N balance in 0.9 m soil profile under three vegetable continuous cropping systems with different water and fertilizer treatments in 2001.
Treatments ManuFertNirrNnetNupNgasNleaNUE
kg kg−1
Nbalance
kg N ha−1
kg N ha−1
W1N1cauliflower120450.011.332.8156.8113.326.664.7317.4
amaranth0100.04.233.961.616.00.0289.960.5
spinach0309.04.836.0104.249.930.4163.1165.3
total120859.020.3102.7322.6179.257.0128.5543.2
W1N2cauliflower120109.511.333.5135.844.013.883.780.7
amaranth00.04.234.655.67.00.0297.1−23.8
spinach0126.34.835.0102.922.98.4225.031.9
total120235.820.3103.1294.373.922.2166.588.8
W1N3cauliflower12071.611.329.4122.634.48.287.867.1
amaranth00.04.234.053.08.50.0297.6−23.3
spinach099.64.835.4102.519.315.2229.22.8
total120171.220.398.8278.162.223.4176.546.6
W2N1cauliflower120450.08.130.5156.681.91.180.6369.0
amaranth0100.02.933.762.419.50.0273.554.7
spinach0309.03.334.0112.932.10.6217.0200.7
total120859.014.398.2331.9133.51.7156.9624.4
W2N2cauliflower12062.08.130.5141.122.20.0116.457.3
amaranth00.02.933.264.43.20.0313.6−31.5
spinach098.93.335.0108.112.10.5214.616.5
total120160.914.398.7313.637.50.5188.042.3
W2N3cauliflower12068.58.129.4139.825.70.2112.360.3
amaranth00.02.933.151.84.50.0270.0−20.3
spinach071.83.334.5112.310.63.5240.5−16.8
total120140.314.397.0303.940.83.7184.323.2
Note: Manu, manure; Fert, fertilizer; Nnet, net mineralization; Nirr, N input by irrigation; Nup, crop N uptake. Nlea, nitrate leaching. Ngas, gaseous N emission. NUE (nitrogen use efficiency) = Fresh yield/(Nup + Nlea + Ngas). Nbalance = N input − N output.
Table 6. Growth periods of different vegetables in 2000 and 2001.
Table 6. Growth periods of different vegetables in 2000 and 2001.
CauliflowerAmaranthSpinach
YearSowing DateHarvest DateSowing DateHarvest DateSowing DateHarvest Date
20003 April6 June1 July30 July4 September14 October
200113 April6 June24 June22 July6 September29 October
Table 7. Soil physical and chemical properties in 0-120 cm soil profile.
Table 7. Soil physical and chemical properties in 0-120 cm soil profile.
Soil Layer (cm)Particle Fraction (%)Soil
Texture
BD
(g cm−3)
θfc
(cm3 cm−3)
θwp
(cm3 cm−3)
θs
(cm3 cm−3)
pH
SandSiltClay
0–3032.750.217.1Loam1.330.330.140.507.9
30–6016.456.826.8Silt loam1.520.350.150.437.9
60–9030.747.621.7Loam1.430.350.140.448.0
90–12035.046.418.6Loam1.620.330.140.438.0
Note: BD, bulk density; θs, saturated water content; θfc, field capacity; θwp, wilting water content.
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

Sun, Y.; Chen, S.; Feng, P.; Chen, Q.; Hu, K. Assessing the Fates of Water and Nitrogen on an Open-Field Intensive Vegetable System under an Expert-N System with EU-Rotate_N Model in North China Plain. Plants 2024, 13, 2150. https://doi.org/10.3390/plants13152150

AMA Style

Sun Y, Chen S, Feng P, Chen Q, Hu K. Assessing the Fates of Water and Nitrogen on an Open-Field Intensive Vegetable System under an Expert-N System with EU-Rotate_N Model in North China Plain. Plants. 2024; 13(15):2150. https://doi.org/10.3390/plants13152150

Chicago/Turabian Style

Sun, Yuan, Shaoqing Chen, Puyu Feng, Qing Chen, and Kelin Hu. 2024. "Assessing the Fates of Water and Nitrogen on an Open-Field Intensive Vegetable System under an Expert-N System with EU-Rotate_N Model in North China Plain" Plants 13, no. 15: 2150. https://doi.org/10.3390/plants13152150

APA Style

Sun, Y., Chen, S., Feng, P., Chen, Q., & Hu, K. (2024). Assessing the Fates of Water and Nitrogen on an Open-Field Intensive Vegetable System under an Expert-N System with EU-Rotate_N Model in North China Plain. Plants, 13(15), 2150. https://doi.org/10.3390/plants13152150

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