Next Article in Journal
On-Farm Multi-Environment Evaluation of Selected Cassava (Manihot esculenta Crantz) Cultivars in South Africa
Next Article in Special Issue
Enhanced-Efficiency Fertilizers Impact on Nitrogen Use Efficiency and Nitrous Oxide Emissions from an Open-Field Vegetable System in North China
Previous Article in Journal
Growth, Biomass Partitioning, and Photosynthetic Performance of Chrysanthemum Cuttings in Response to Different Light Spectra
Previous Article in Special Issue
The Nitrogen Dynamics of Newly Developed Lignite-Based Controlled-Release Fertilisers in the Soil-Plant Cycle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simulating the Effects of Different Textural Soils and N Management on Maize Yield, N Fates, and Water and N Use Efficiencies in Northeast China

1
College of Resources and Environmental Sciences, Key Laboratory of Straw Comprehensive Utilization and Black Soil Conservation, Jilin Agricultural University, Changchun 130118, China
2
College of Land Science and Technology, Key Laboratory of Arable Land Conservation (North China), Ministry of Agriculture and Rural Affairs, China Agricultural University, Beijing 100193, China
*
Authors to whom correspondence should be addressed.
Plants 2022, 11(23), 3338; https://doi.org/10.3390/plants11233338
Submission received: 30 October 2022 / Revised: 23 November 2022 / Accepted: 28 November 2022 / Published: 1 December 2022
(This article belongs to the Special Issue Water and Nitrogen Management in Soil-Crop System)

Abstract

:
Determining the best management practices (BMPs) for farmland under different soil textures can provide technical support for improving maize yield, water- and nitrogen-use efficiencies (WUE and NUE), and reducing environmental N losses. In this study, a two-year (2013–2014) maize cultivation experiment was conducted on two pieces of farmland with different textural soils (loamy clay and sandy loam) in the Phaeozems zone of Northeast China. Three N fertilizer treatments were designed for each farmland: N168, N240, and N312, with N rates of 168, 240, and 312 kg ha−1, respectively. The WHCNS (soil Water Heat Carbon Nitrogen Simulator) model was calibrated and validated using the observed soil water content, soil nitrate concentration, and crop biological indicators. Then, the effects of soil texture combined with different N rates on maize yield, water consumption, and N fates were simulated. The integrated index considering the agronomic, economic, and environmental impacts was used to determine the BMPs for two textural soils. Results indicated that simulated soil water content and nitrate concentration at different soil depths, leaf area index, dry matter, and grain yield all agreed well with the measured values. Both soil texture and N rates significantly affected maize yield, N fates, WUE, and NUE. The annual average grain yield, WUE, and NUE under three N rates in sandy loam soil were 8257 kg ha−1, 1.9 kg m−3, and 41.2 kg kg−1, respectively, which were lower than those of loam clay, 11440 kg ha−1, 2.7 kg m−3, and 46.7 kg kg−1. The order of annual average yield and WUE under two textural soils was N240 > N312 > N168. The average evapotranspiration of sandy loam (447.3 mm) was higher than that of loamy clay (404.9 mm). The annual average N-leaching amount of different N treatments for sandy loam ranged from 5.1 to 13.2 kg ha−1, which was higher than that of loamy clay soil, with a range of 1.8–5.0 kg ha−1. The gaseous N loss in sandy loam soil accounted for 14.7% of the fertilizer N application rate, while it was 11.1%in loamy clay soil. The order of the NUEs of two textural soils was: N168 > N240 > N312. The recommended N fertilizer rates for sandy loam and loamy clay soils determined by the integrated index were 180 and 200 kg ha−1, respectively.

1. Introduction

Maize is one of the three major grain crops, with a yield ranking first in the world [1]. Maize production in China accounts for 20% of the total maize yield and 22% of the maize planting area in the world. The maize planting area in Northeast China is about 12.48 million hectares, accounting for 32.3% of the Chinese maize planting area and 30.2% of the total maize yield [2]. In recent years, due to excessive reclamation and unreasonable cultivation, the black soil region of Northeast China has experienced serious soil degradation. The soil fertility has decreased significantly, which has seriously hindered agricultural sustainable development and also caused a negative impact on food security in China [3].
Nitrogen (N) is the main nutrient element for maize growth, and also a major factor influencing crop yield. Since the 1990s, farmers have applied excessive N fertilizer in pursuit of high yield [4]. The amount of N applied is far higher than the N demand of crops, even twice the crop N requirement in some areas [5]. A large number of experimental studies have shown that excessive N application cannot improve crop yield, but reduces N fertilizer utilization and brings about a series of environmental problems. First of all, excess N accumulates in the soil profile. Once irrigation and heavy rainfall occur, soil residual N will easily leach into the groundwater, causing the nitrate content in groundwater to exceed the drinking water standard [6]. Secondly, a large amount of fertilizer N can be lost to the atmosphere in the forms of NH3, N2O, and other gases [7]. Therefore, it is important to determine the best management practices (BMPs) for improving crop yield and N use efficiency, reducing environmental pollution and maintaining sustainable agricultural development in the black soil region of Northeast China.
Most of the black soil region in Northeast China is a typical rain-fed agricultural region. Water and nutrients are the two main factors limiting crop productivity in this region. Soil quality has a significant influence on soil water holding capacity, water availability, hydraulic conductivity, nutrient transport, and transformation processes [8], which ultimately lead to different crop yields. Generally, there are many capillary micro-pores in clay or silty clay soil, with strong water- and nutrient-holding capacities, which can lead to a surplus of sufficient water and N for crop growth [9]. However, sandy or sandy loam soil has a low water holding capacity and larger hydraulic conductivity, which easily leads to water drainage and N leaching. Soil water and solute have the characteristics of preferentially passing through macropores. With the increase in soil clay, soil residual nitrate content shows an increasing trend. The abundant meso- and micro-pores in loamy clay lead to an increased retention of nitrate. Although the availability of water and nitrate is not better than that of sandy loam, it can be kept longer in soil [10]. In addition, soil texture directly affects soil’s available water content [11], further affecting soil N transformation, crop absorption, and N loss. Studies have shown that the N mineralization rate in sandy soil is usually higher than that in silty or clay soil due to the stronger aeration of active organic matter and less physical protection. Compared with sandy loam soil, clay soil has a stronger fixation capacity for NH4+-N. The N absorption efficiency and response to N fertilizer of maize in clay soil are also more sensitive [12]. Therefore, even if the fertilizer N applied rate is the same, soil N concentration and its distribution will also vary depending on the soil texture. This may affect the relationship between N supply and maize growth under specific field conditions. Many studies have shown that the effects of soil texture on soil N transformation, crop yield, N absorption, mineral N residue, and N loss are different [13]. Therefore, soil texture is a key factor affecting crop yield and water- and N-use efficiencies (WUE and NUE).
Optimizing fertilization is helpful to improve crop yield, WUE, and NUE and reduce N loss. However, it is difficult to accurately determine the N demand of crops at different growth stages. The soil–crop system model can simulate crop growth and the dynamic of soil water and N under different climates, soil properties, and field management practices. At present, some models are used to simulate crop growth and N leaching in the black soil region of Northeast China and optimize water and N management. Liu et al. (2012) used the DSSAT-CERES-Maize model to simulate the aboveground dry matter, maize yield, and N uptake in Gongzhuling, Jilin Province, in order to better predict the maize yield [14]. Yang et al. (2013) used the DSSAT model to simulate the effects of the combined application of organic and inorganic fertilizer on maize yield, soil carbon, and N dynamics [15]. The results showed that the combined application of inorganic and organic fertilizer with a total N rate of 165 kg ha−1 increased the content of soil mineral N, and nitrate leaching decreased by 10 kg ha−1 compared with the treatment without organic fertilizer. Li et al. (2016) used the RZWQM model to simulate the effect of different textural soils on the potential yield of maize. The results showed that, under the same weather conditions, the soil textures affected the potential yield through soil moisture content, and the average yield of black soil (8 Mg ha−1) was higher than that of eolian sandy soil (7 Mg ha−1) [16]. Jiang et al. (2019) used DNDC and DSSAT models to simulate maize yield under different fertilization scenarios. The results showed that 180–210 kg ha−1 was more suitable for maize production in Northeast China [17]. Peng et al. (2022) used an APSIM model to evaluate and optimize the long-term impact of different soil textures on WUE and economic benefits under alfalfa–maize rotation [18]. At present, these models are mainly used to simulate the spring maize in the black soil region of Northeast China, taking into account the factors such as rain-fed conditions, straw returning, and fertilization. However, relevant studies on crop growth, N fate, and the optimization of water- and N management practices in farmlands with different soil textures are rare.
Recently, according to the characteristics of climate, soil, and field management in China, Liang et al. (2016) developed a model for the integrated water and N management of a soil and crop system (WHCNS, soil Water Heat Carbon Nitrogen Simulator) [19], which has been successfully applied to the water use, N loss, WUE, and NUE of crops in different regions such as Northwest, North and South China [20,21,22]. However, the model has not been verified and applied in the black soil area of Northeast China. Therefore, this study aims to achieve the following: (i) evaluate the applicability of the WHCNS model in simulating the fate of water and N and crop growth processes in farmland with different textural soils in the black soil region of Northeast China; (ii) simulate and analyze the effects of different textural soils and N fertilizer application on maize yield, water consumption, N fate, WUE, and NUE; and (iii) determine the reasonable amount of N fertilizer for farmlands with different textural soils by comprehensively considering agronomic effects, environmental effects, and economic benefits.

2. Results

2.1. Model Calibration and Validation

In this study, the two-year field-measured data (soil water content, nitrate N concentration, LAI, dry matter weight, and crop yield) from the N312 treatment with two textural soils were used to calibrate the model. The calibration process was carried out in the following order: soil moisture, soil N, and crop growth process. The “trial and error method” was used to adjust the model parameters until the simulated values were as consistent as possible with the measured values [23]. After calibration, all parameters were fixed, and the model was verified with measured data from other treatments. Due to space limitation, this paper only shows the comparison results of the simulated and measured values of the N312 treatment of loamy clay from 2013 to 2014, and the rest of the results are shown in the appendix (Figures S1–S20).
The comparison results of the simulated and measured soil water content and nitrate concentration at different depths under N312 treatment in loamy clay soil are shown in Figure 1 and Figure 2, respectively. As seen in Figure 1, the peak value of soil water content was in good agreement with the rainfall events. Except for a few points, the simulated values of soil moisture in most periods agreed well with the measured values. However, the simulated soil water content of each layer in sandy loam soil fluctuates more drastically than that of loamy clay soil (Figure 1 and Figure S2), which was mainly related to the strong hydraulic conductivity and low water holding capacity of sandy loam.
Figure 2 shows that the nitrate content of loamy clay soil changed greatly when it was above 40 cm, while the change in the 40–100 cm soil layer was relatively slight, which was related to the effects of precipitation, fertilization, and other factors on the surface soil, followed by the mineralization of organic matter in the soil and the activities of microorganisms mainly in the 0–30 cm surface soil. The nitrate concentration fluctuated with the events of fertilization, which was similar to soil water content (Figure 1 and Figure 2). In the jointing stage, urea was applied as topdressing at the rate of 208 kg ha−1, and the nitrate concentration in the surface soil simulated by the model showed an obvious peak at the corresponding time, indicating that the simulation results matched well with the measured values. The comparison results between the simulated and measured values of nitrate in each soil layer of the two textural soils for validation treatments are shown in the appendix (Figures S6, S10, S14 and S18). In general, the simulated values are in good agreement with the measured values.
Figure 3 and Figure 4 show the simulated results of LAI and total dry matter, respectively. Figure 3 shows that the model performance of LAI under N312 treatment is perfectly good. The simulated LAI of N168 treatment was lower than that of the other two treatments, which might be caused by insufficient N fertilizer supply (Figures S15 and S19). It can be seen from Figure 4 that the simulated values of total dry matter under N312 for loamy clay soil were also in good agreement with the measured values.
The nRMSE, E, and d of yield, TDM, and crop N uptake of calibration treatment (N312) for sandy loam were 14.26–26.56%, 0.42–0.86, and 0.76–0.98, respectively. The loamy clay were 6.61–22.56%, 0.44–0.88, and 0.79–0.96, respectively; The results of the validation treatments (N168 and N240) for sandy loam were 16.52–29.62%, 0.36–0.46, and 0.71–0.96, respectively; The values of loamy clay were 10.56–28.64%, 0.38–0.98, and 0.73–0.98, respectively. The nRMSE values of yield, TDM, and crop N uptake of sandy loam and loamy clay were all less than 30% (Table 1). According to the thresholds suggested by van Liew and Garbrecht [24], an acceptable simulation should have a value of E>0.36 and d>0.70. The E values of yield, TDM, and crop N uptake were all greater than 0.36, and the d values were all greater than 0.70. These results indicated that the WHCNS model performed well in simulating the yield, TDM, and crop N uptake in this black soil.
Figure 5 shows the correlation between the simulated and measured soil water content, nitrate content, yield, and LAI of the soil clay validation treatment. The determination coefficients of the four indicators reached 0.99, 0.94, 0.87, and 0.99, respectively. Their determination coefficients for sandy loam soil are 0.80, 0.80, 0.99, and 0.97, respectively (Figure S21), and their slopes of regression equations are close to 1, which indicated that the model performance was good, and the WHCNS model can be used to simulate and analyze the water consumption, N fate, and crop growth of farmlands with different textural soils in the black soil region of Northeast China.

2.2. Crop Yield, Water Consumption, and Water Use Efficiency under Two Textural Soils

According to the measured yields (Table 2), the yields of the three treatments of sandy loam in 2014 decreased by 30% compared with that in 2013, while those of loamy clay increased by 12%. The annual average yields of sandy loam (8257 kg ha−1) decreased by 28% compared with that of loamy clay (11,440 kg ha−1). Under the same N fertilizer rate condition, the yield of loamy clay farmland was higher than that of sandy loam. The order of the annual average yields was N240 > N312 > N168. The annual average yields of N240 (10,341 kg ha−1) treatment in two years were 10% and 5% higher than those of N168 (9367 kg ha−1) and N312 (9838 kgha−1), respectively (Table S2). This shows that excessive N fertilizer application would not increase maize yield, but lead to yield reduction.
Table 2 shows the simulation results of the water balance of a 0–1 m soil profile for two textural soils and different N rates. The precipitations in 2013 and 2014 was 615 mm and 420 mm respectively. The water drainage and runoff in 2014 were much lower than that of 2013, mainly because there was no irrigation in this area, and the main water resource came from rainfall. With the increase in rainfall, the corresponding water drainage and runoff also increased. ET was the largest water consumption, accounting for 86.4% and 78.2% of the water input from sandy loam and loamy clay in farmland, respectively. The WUE of sandy loam was 1.8 kg m−3, and that of loamy clay was 2.8 kg m−3. The WUE of loamy clay was significantly higher than that of sandy loam. The order of WUE under two textural soils and different N fertilizer rates was N240 > N312 > N168.
Table 3 shows the simulated N balance in a 1 m soil profile under different N fertilizer rates. N output from farmland mainly includes crop N uptake, N leaching, and gaseous N loss. The amount of crop N uptake is the largest output of N, accounting for 54.1% and 70.6% of the N input in sandy loam and loamy clay, respectively. The amount of crop N uptake in loamy clay was significantly higher than that in sandy loam under the same N fertilizer rate (Table S2). N leaching was mainly affected by the N fertilizer rate and rainfall, but the increasing proportion for different textural soils was different. When the N fertilizer rate increased from 168 to 312 kg ha−1, the N-leaching amount of sandy loam increased by 162.2%, while it increased by 173.0% for loamy clay soil. Under the same N fertilizer rate, the N-leaching amount of sandy loam increased by 6.4–16.2 kg ha−1 compared with loamy clay in the abundant rainfall year (2013), with an increasing ratio of 165.9−177.0%. In the normal rainfall year (2014), under the same N fertilizer rate, the N-leaching amount of sandy loam was only 0.09–0.32 kg ha−1 more than that of loamy clay. The annual average N-leaching amount of sandy loam was 170.3% higher than that of loamy clay, indicating that determining a reasonable N fertilizer rate for different textural soils was the key to ensuring crop yield and environmental friendly.
Gaseous N emission is another main pathway of N loss, mainly including NH3 volatilization and denitrification (Table 3). The NH3 volatilization of both textural soils increased with the increase in the N fertilizer rate. The NH3 volatilization amounts of the N168, N240, and N312 treatments were 21.8, 32.3 kg, and 44.3 kg ha−1, respectively. Compared with the N168, the amounts of N240 and N312 increased by 32.7 and 103.5%, respectively. The average NH3 volatilization of sandy loam was 12.2 kg ha−1 higher than that of loamy clay, and the total gaseous N loss in sandy loam and loamy clay accounted for 14.6 and 11.1% of N fertilizer input, respectively.
From the results of NUE, the NUE for both textural soils decreased with the increase in the N fertilizer rate: N168 > N240 > N312. The NUE was negatively correlated with total N loss and positively correlated with maize yield (Table 3). The N loss of N312 was relatively high, resulting in low NUE. In addition, there was a large difference in NUE between different years. In 2013, the NUEs of sandy loam and loamy clay were 30.2 and 40.3 kg kg−1, respectively, while in 2014, they were 48.8 and 50.3 kg kg−1, respectively.

2.3. Optimizing N Fertilizer Management for Two Textural Soils under Rainfed Conditions

From the simulation results, soil texture, rainfall, and fertilization management have obvious effects on maize yield, water consumption, WUE, N leaching, gaseous N loss, and NUE. In order to assess the agronomic and environmental effects, the integrated index of each treatment was calculated (Table 3). It can be seen that the N168 treatment on loamy clay and sandy loam was the best among the three N fertilizer treatments. In order to further optimize the N fertilizer rate, the normal rainfall year 2014 was selected as a case study, and different N fertilizer management scenarios were set. Under the condition of the same field management practices, a total 20 N fertilizer management scenarios for each textural soil were simulated using the WHCNS model.
The simulation results showed that with the increase in the N fertilizer rate, the maize yield for two textural soils was also increasing at first (Figure 6a); when the N fertilizer rate exceeded 180 kg ha−1 for sandy loam soil, the yield reached the highest, which was 7180 kg ha−1, while that of loamy clay reached the highest when the N fertilizer rate reached 200 kg ha−1, which was 11164 kg ha−1. The ET of loamy clay was 7.0% higher than that of sandy loam, while water drainage was 88.9% lower than that of sandy loam. WUE also showed a similar trend. WUE increased with increases in the N fertilizer rate. When the N fertilizer rate of sandy loam exceeded 180 kg ha−1, and that of loamy clay exceeded 200 kg ha−1, WUE stayed at a constant level. (Figure 6b).
With the increase in N fertilizer rate, the N leaching of sandy loam slowly increased at first, and then increased sharply (Figure 6c), while it increased slightly for loamy clay soil. Ammonia volatilization in both textural soils increased exponentially, and it was higher for sandy loam than that for loamy clay (Figure 6d). The NUE of sandy loam and loamy clay increased linearly with the N fertilizer rate at first. When the N fertilizer rate of sandy loam exceeded about 180 kg ha−1, and the N fertilizer rate of loamy clay exceeded about 200 kg ha−1, the NUE decreased rapidly with the increase in the N rate (Figure 6e).
According to the simulation results in 2014, the integrated index of 20 fertilizer N management scenarios for each textural soil was calculated (Table 4). The integrated index of sandy loam and loamy clay soils was the lowest when the fertilization rate was 40 kg ha−1. Under this scenario, the WUEs of sandy loam and loamy clay soils were 0.3 kg and 0.1 kg m−3, respectively. Their total N loss was 6.2 and 0.46 kg ha−1, and their crop yields only reached 17.4 and 5.0% of the maximum. Their NUEs were 16.7 and 7.9 kg kg−1, respectively. When the N fertilizer rate was 180 kg ha−1 for sandy loam soil, and it was 200 kg ha−1 for loamy clay soil, their WUEs were 1.8 and 2.7 kg m−3, respectively, and their NUEs reached 32.0 and 40.9 kg kg−1, respectively. The crop yields reached the maximum, and the integrated index also reached the highest, which could be recommended as BMPs of the two textural soils for maize yield in this region.

3. Materials and Methods

3.1. Experimental Site Description

Field experiments were conducted from April 2013 to September 2014 in Sankeshu (43°20′17.4″ N,124°00′29.1″ E) and Fujiajie (43°21′48.1″ N,124°05′02″0 E) experimental sites with a distance of no more than 3km, which are located in Lishu County, Jilin Province of China. The region is characterized as a semi-humid continental monsoon climate in the north temperate zone. The annual mean temperature is 12 °C and the mean annual precipitation is 531 mm; 70–80% of the annual precipitation occurs from July to September. The soil textures of Sankeshu and Fujiajie are loamy clay and sandy loam, respectively. Soil physical–chemical properties in the soil profile for each site are listed in Table 5, and soil surface chemical properties are listed in Table S1.

3.2. Experiment Design

The spring maize variety Liangyu 99 was sown on 25 April 2013 and 1 May 2014 at a density of approximately 65,000 plants ha−1, and harvested on 5 October 2013 and 20 September 2014, respectively. On each soil texture, the field experiment was designed as a randomized block with three replicates of three fertilizer N rates: 168, 240, and 312 kg ha−1 by urea, which are respectively recorded as N168, N240, and N312. One-third of the fertilizer N was applied in furrows before ploughing as basal fertilizer, together with P (100 kg P2O5 ha−1 as Ca(H2PO4)2) and K (120 kg K2O ha−1 as K2SO4). The remaining N fertilizer was side-dressed in the topsoil layer at the jointing stage. Other farmland management practices were the same as local farmers.

3.3. Field Sampling and Laboratory Analysis

Soil profile pits were excavated to 100 cm soil depth. Soil samples for each 20 cm layer were collected before spring maize sowing. The bulk density and soil texture were measured. The soil hydraulic parameters were estimated by using the Pedo transfer function (PTF) method [25]. The 0–100 cm soil profile samples were collected at the key growth stages of spring maize (jointing stage, silking stage, filling stage, and maturity stage) for each plot, and the soil moisture content was determined by the drying method. Fresh soil samples were extracted with 1 mol L−1 KCl to determine the concentrations of NH4+-N and NO3-N by using a continuous flow analyzer (AA3, Bran and Luebbe, Norderstedt, Germany).
At the same time, at the key growth stage of spring maize, both the leaf length and width were measured for all the leaves from ten representative plants in each plot to calculate leaf area index (LAI), and they were dried at 70 °C to determine the aboveground biomass. The dried plant samples were then grounded into powder to measure the N content of various organs (stems, leaves, and grains) by Kjeldahl Apparatus (KDY-9820, KETUO, China). At harvest time, the total number of cobs and total cob weight of maize were recorded in the production area. A total of 10 representative plants were selected according to the average weight of a single cob to determine the maize yield (calculated by 14% water content). Meteorological data were obtained from the local small meteorological station, mainly recording daily rainfall, maximum temperature, minimum temperature, average temperature, relative humidity, daily average wind speed, etc.

3.4. Model Description and Inputs

The WHCNS model includes five main modules: soil water, soil temperature, soil carbon, soil N, and crop growth. The model is driven by daily meteorological data and crop biological parameters, and computes daily processes including soil evaporation, crop transpiration, soil water movement, runoff, soil temperature, soil net mineralization, nitrification, ammonia volatilization, denitrification, and crop growth. The potential evapotranspiration is calculated using the Penman–Monteith method [26]. The water infiltration and redistribution in the soil profile were described by the Green–Ampt method [27] and Richards equation, respectively. The soil heat-transfer simulation is taken from the HYDRUS-1D model, using the convection–dispersion equation [28]. The soil C-and N-cycling concepts are taken from the DAISY model [29]. The crop model is taken from the PS123 model, which simulates crop growth processes, including crop development stage, LAI, photosynthate accumulation and allocation, maintenance respiration, and growth respiration [30]. The water stress factor is calculated by the actual transpiration divided by the potential transpiration. The N stress factor is computed based on the simulation of the crop N demand, crop N uptake, and actual soil N supply.
This model can be used to simulate and analyze the effects of various field management practices (irrigation, fertilization, mulching, tillage, straw returning, etc.) on farmland water consumption, N fate, organic matter turnover, and crop growth [31,32]. The basic input data include the latitude and altitude of the experiment site, basic physical and chemical properties of the soil, C/N ratio, field management, initial soil conditions, and daily meteorological data. For more detailed processes, see the literature [33].

3.5. Evaluation of Model Simulation Effect

In this study, three statistical indices were used to assess the agreement between the predicted and observed values:
(i)
The normalized root-mean-squared error (nRMSE):
n R M S E = 1 n i = 1 n ( P i O i ) 2 × 100 x ¯        
(ii)
Nash–Sutcliffe efficiency (E):
E = 1 i = 1 n P i O i 2 i = 1 n O i O 2    
(iii)
Agreement index (d):
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 data pairs. Pi and Oi are the predicted and observed values, respectively. O is the mean of the observed data. nRMSE represents the percentage of the average deviation to the measured mean. nRMSE < 15%, nRMSE = 15–30%, and nRMSE > 30% are considered “good,” “moderate,” and “poor” agreement, respectively [31]. E ranges between −∞ and 1 and allows for a negative error. The d value is a descriptive measure that ranges from 0 to 1, wherein a value close to 1 indicates satisfactory model performance. Quantities of d ≥ 0.75 and E ≥ 0 are minimum threshold values for crop growth. Quantities of d ≥ 0.60 and E ≥ −1 are minimum threshold values for N-output evaluation [34].

3.6. Integrated Index Method

Quantitative analysis of model outputs (crop yield, water drainage, WUE, N leaching, gaseous N loss and NUE) is helpful to obtain BMPs. In this study, the integrated index method was used to obtain BMPs in the three treatments [31]. In order to maximize crop yield and economic benefits, and minimize environmental impact, the evaluation indicators include agronomy factors (AF), environment factors (EF), and value-to-cost ratio (VCR). Referring to the research of Hu et al. (2010), crop yield, WUE, and NUE were selected as input variables to calculate AF [10]. The weights of crop yield, WUE, and NUE are set to +0.6, +0.2, and +0.2, respectively, and the values of these variables are normalized to a range of 0–1. Subsequently, AF was calculated by summing the products of the three normalized variables and their corresponding weights. Gaseous N loss and nitrate leaching were used to calculate EF, and their weights were set as +0.5 and +0.5. EF was calculated in the same way as AF. The VCR was calculated using the following equation:
V C R = Y × Y P F × F P        
where F is the amount of fertilizer (kg), FP is the price of the fertilizer (yuan kg−1), Y is the crop yield (kg ha−1), and YP is the price of the grain (yuan kg−1). The prices of the fertilizer and maize are 4.5 yuan and 2.0 yuan kg−1, respectively. In order to calculate the integrated index, the weights of AF, EF and VCR are set to + 0.7, −0.2, and +0.5 respectively [31].

3.7. N Fertilizer Scenarios

In order to analyze the impact of different N fertilizer management practices on crop yield and N fate under two textural soils, the selected simulation years should be able to represent the rainfall situation in most of the local years. This study took 2014 as an example, where the rainfall year type was normal, and set up different N fertilizer management practices scenarios. As it is rainfed, irrigation scenarios are not considered. The specific N fertilizer input scenarios are as follows: the range of N fertilizer rate is 20–400 kg ha−1, and the interval of the N fertilizer rate is 20 kg ha−1, so the total number of N fertilizer rates is 20 for each textural soil. Then, the calibrated WHCNS model was used to simulate crop yield, water drainage, nitrate leaching, gaseous N loss, WUE, and NUE under each scenario

4. Discussion

4.1. Effects of Soil Texture and N Management on Crop Yield, ET, and WUE

Soil texture affects the water movement and N transport, thus affecting crop yield [35]. Related researchers have shown that physical properties of soil, such as soil texture and soil water retention, can lead to differences in crop yields, such as the yield of crops in loamy soil with a relatively higher water holding capacity [36]. In this study, the average maize yield of loamy clay was 11,440 kg ha−1, which was significantly higher than that of 8257 kg ha−1 for sandy loam (Table S2), which was consistent with other researches results [37], but lower than that of the results reported by Lu et al. [35].
Crop yields were most sensitive to precipitation variability in coarse-textured soils and less responsive to these weather parameters in medium- and fine-textured soils. Increasing SOC concentration (>2%) contributed to crop yields being less sensitive to precipitation [38]. For example, the rainfall in July 2014 was only 58.1 mm, which is a severe drought, resulting in a significant decline in the maize yield of sandy loam soil, with an annual average yield of 30% less than that of 2013. In contrast, the yield of loamy clay increased by 12% in 2014, compared with sandy loam, due to increased water retention, presumably due to increased buffering capacity against increased water loss through evapotranspiration in loamy clay. However, sandy loam has strong air permeability and lower SOC (Table S1), and the growth of crops on sandy loam is more vulnerable to drought stress, thus affecting the yield. Huang et al. (2020) found excessive N fertilization could not achieve a higher yield in Northeast China [39]. In this study, the annual average yields of N168, N240, and N312 were 9367 kg, 10341, and 9838 kg ha−1, respectively (Table S2), and their order was N240 > N312 > N168. The results showed that fertilization had an obvious effect on crop yield, but there exists a law of diminishing returns.
The hydraulic conductivity of sandy loam was high and its water holding capacity was poor, which led to higher ET and lower WUE [39]. Wang et al. (2020) found that ET and WUE under ridge furrow cultivation with medium or fine soil texture increase by 29 and 17%, respectively, compared with coarse soil [40]. Simulation results showed that water drainage was directly related to rainfall and soil texture. In 2013, the precipitation was abundant (615 mm), and ET and water drainage were the main water consumption. The average ETs of sandy loam and loamy clay were 503 and 394 mm, respectively, and their average water drainages were 157mm and 140mm, respectively. In 2014, due to the low rainfall, the difference between the water drainage and ET of the two textural soils was not obvious. The water consumption of sandy loam was higher than that of loamy clay. The main reason was that loamy clay has the stronger water holding capacity and buffer capacity to soil alternates between dry and wet, and the soil profile could be rapidly replenished after rainfall, resulting in lower ET. The loamy clay soil can improve WUE by reducing soil water drainage [12]. The difference in nutrient supply capacity caused by different soil textures also affects the WUE, and the nutrient supply capacity is positively correlated with WUE. The mineral nutrients in the loamy clay are sufficient, and its transpiration coefficient of plant is smaller, so its WUE is higher [9]. In this study, the average WUE of sandy loam is higher than that of loamy clay. The WUE of sandy loam is slightly lower than that of other research results at a nearby site, while that of loamy clay is higher (i.e., 1.9 and 2.7 kg m−3, respectively), mainly due to the influence of rainfall under different meteorological conditions [29].

4.2. Effects of Soil Texture and N Management on N Loss

Soil texture plays an important role in nitrate retention, transport, and leaching [41]. Su et al. (2014) observed that nitrate leaching was the lowest in loam soil sand highest in loamy sand [42]. Michalczyk et al. (2020) used the HERMES model to simulate the long-term effect of different textural soils on N leaching, and found that the amount of nitrate leaching was as follows: sandy loam > silty loam > clay loam, and the nitrate leaching decreased from 113.9 to 41.9 kg ha−1 [43]. Milroy et al. (2008) proposed that leaching varied markedly between the soil textures, and N leaching in loamy sand accounted for 20% that of acid loamy sand or sand [44]. Wang et al. (2022) pointed out that N leaching loss in the sandy-textured soils with high sand content was 1.65 times higher than that with low sand content [45]. In this study, the annual average nitrate leaching of each treatment of sandy loam is 9.3 kg ha−1, which is 3.4 kg ha−1 higher than that of loamy clay. This is due to soil texture combined with rainfall dynamics altering soil moisture dynamics, and consequently regulating soil N responses to precipitation changes. The clay loam soil more effectively buffered water stress during relatively long intervals between precipitation events, particularly after a large rainfall event. The large hydraulic conductivity of coarse texture soil and its poor ability to retain water and nutrients mean it is more likely to cause nitrate leaching [41,46]. In addition, optimizing the amount of N fertilizer for different soil textures can effectively reduce nitrate leaching. Azad et al. (2020) reported that the amount of nitrate leaching in sandy soil decreased from 122.18 to 44.98 kg ha−1, and loam and clay decreased from 22.1 to 13.9 kg ha−1 by optimized fertilization [47]. Reducing the N fertilizer rate can significantly reduce nitrate leaching. Zhang et al. (2015) showed that the 150–240 kg ha−1 fertilizer rate could maintain the high maize yield meanwhile reducing nitrate leaching [48]. These results in the literature are clearly supported by our findings. In our study, nitrate leaching was related to the soil texture and N application rate. The nitrate leaching of a 1 m soil profile of sandy loam and loamy clay increases with increases in the N fertilizer rate. When the N fertilizer rate of sandy loam exceeds 260 kg ha−1, the increase rate is fast (Figure 6). The risk of nitrate leaching is further increased under the condition of a high N fertilizer rate. The maximum amount of nitrate leaching for N312 is 25.9 kg ha−1. The N leaching is lower than that in North China, which may be related to the frequent irrigation in North China [49].
Gaseous N loss (NH3 volatilization and denitrification) is another main pathway of N loss in farmland. Ammonia volatilization in farmland is closely related to soil texture, fertilizer application rate, application method, etc. Soil texture and pH have different influences on NH3 volatilization. Awale et al. (2016) found that, across N amendments, cumulative NH3 losses from sandy loam soil ranged from 0.7 to 4.3% of applied N, and were higher than those in silty clay soil (0.1–0.4% of applied N). Cumulative N2O emissions did not differ between soil textures and ranged from 3.7 to 7.4% of applied N across N sources [50]. In this study, the ranges of NH3 volatilization for sandy loam and loamy clay were 21.8–1−56.5 and 14.5–42.2 kg ha−1, respectively, on average accounting for 13.4% and 8.9% of the N fertilizer rate, respectively. The main reason for that may be that loamy clay has a stronger adsorption capacity for NH4+. The results are higher than previous studies, mainly due to different types of N fertilizer [49]. We also found that NH3 volatilization loss in the northeast black soil region is lower than that in southern rice fields and calcareous soil in the North China Plain, which is closely related to different soil textures, soil pH, and climate conditions [51].
The effect of soil texture on denitrification is mainly through the change in soil water holding capacity, porosity, and aeration. Jamali et al. (2016) found that denitrification increased with soil texture becoming finer [52]. This is influenced by several factors such as the differences in soil physical and chemical properties, climate, and agricultural technology. Gaillard et al. (2016) found that the N2O emission by denitrification ranged from 74 to 98% of the total N2O emission, and the N2O emission from silty sandy loam was the highest, with an average emission raised from 80% to 158% that of loamy sand and sandy loam [53]. Pihlatieet al. (2004) found that, within a certain range, the N2O emission from soil was positively correlated to the soil moisture content in organic soil, clay, and loamy sand [54]. Soil moisture content is directly related to clay content, while soil aeration is negatively correlated to clay content. Soil water content plays an indirect role in regulating N2O emission, mainly by affecting the oxygen supply and demand in soil to control nitrification and denitrification processes. Our research confirmed this point. In this study, the denitrification emission of sandy loam is 3.65 kg ha−1, which is smaller than that of loam clay (6.81 kg ha−1). The results are consistent with the above studies. In addition, due to the difference in water content of sandy loam soil and loamy clay soil caused by rainfall (Figure 1 and Figure S1), the denitrification rate of sandy loam soil is 7.6 kg ha−1 lower, while in 2014, the denitrification rate decreased significantly, and the difference between the two textural soils is not obvious. Our results were also consistent with previous studies.

4.3. Effects of Soil Texture and N Management on N Uptake and N Use Efficiency

The soil texture and N management will result in the difference in crop N uptake. Previous studies have found that soil water content in sandy soil (0–60 cm root zone in 1–2 days after rainfall) for irrigated maize is higher than that in clay or loam soil [12,36]. However, soil water tends to infiltrate into the deeper layers rapidly, which is not conducive for maize root to absorb or utilize water and nutrients. Contrarily, a high water content in clay soil is more favorable for maize to adsorb nitrate. The accumulation of nitrate in the soil root zone is crucial to improve soil N-supplying capacity and crop N uptake [12]. In this study, the crop N uptake from 2013 to 2014 was 121.2–182.4 kg ha−1 in sandy loam and 174.5–249.4 kg ha−1 in loamy clay. The average crop N uptake of the three treatments in loamy clay was 55.04 kg N ha−1 higher than that in sandy loam (Table 3). In addition, except that the N240 treatment had the highest N uptake under sandy loam in 2013, the N uptake by crops increased with increases in the N fertilizer application rate, which was also consistent with the previous studies [55].
Improving NUE is a critical strategy to mitigate the impact of N loss on the environment [56]. Cao et al. (2017) optimized the N fertilizer management of maize on silt loam, reduced the N fertilizer application rate by 36%, and improved the NUE by 20–61% [57]. Zheng et al. (2016) showed that the NUE of summer maize in sandy loam soil was 33.8–47.3% [58]. Optimizing N fertilizer use is necessary to avoid the over-application of N fertilizer, high N loss, or decreased NUE. Li et al. (2007) proposed that reducing fertilizer N input can reduce N loss and improve NUE [59]. It is also possible to achieve greater NUE by adjusting the amount and effective period of fertilization for different soil textures to synchronize crop demand [60]. Sandy loam will usually lead to water drainage and nitrate leaching, thus decreasing NUE. In this study, crop N uptake and NUE were higher on loamy clay than on sandy loam (Figure 5). In addition, N loss increased with increases in N fertilizer application according to the simulation results. When the N fertilizer rate exceeds the critical value, crop yield tends to become stable. Continuing to increase N input will significantly increase N loss and reduce NUE (Figure 6).
The N fertilizer application presents inaccuracy on soil N supply capacity and the yield level of the different regions. Previous studies rarely considered the difference caused by soil textures and climate conditions in different regions An N fertilizer positioning experiment in Jilin Province showed that the N fertilizer application rate for maize should be 180–240 kg ha−1 under the yield level of 12,000–14,000 kg ha−1 maize [61]. Wu et al. (2015) recommended a 114–224 kg ha−1 N fertilizer application rate in the spring maize region of Northeast China [62]. However, in this study, the simulation results showed that the optimal N fertilizer rates for sandy loam and loamy clay soil were 180 and 200 kg ha−1, respectively. Compared with the conventional N fertilizer rate of 221 kg N ha−1, the N fertilizer rate decreased by 18.6 and 9.5%, respectively (Table 4), without reducing maize yield due to the reduction in the N fertilizer rate.

5. Conclusions

The WHCNS model was calibrated and evaluated based on field-measured data for spring maize under different N fertilizer management and two textural soils. The results indicated that the WHCNS model could be used to quantitatively analyze the maize yield, water consumption, N fates, and WUE and NUE, and recommended the BMPs for maize yield in the black soil region of Northeast China.
The simulation results showed that average maize yields, WUE, and NUE were 3984 kg ha−1, 0.9 kg m−3, and 16.1 kg kg−1 higher in loamy clay soil than those of sandy loam soil, respectively, while the total N loss in loamy clay soil was 2.4 kg ha−1 lower than that in sandy loam soil. Compared with sandy loam, maize planted in loamy clay had a higher yield, crop N uptake, WUE, and NUE, and lower N loss.
Under the same N fertilizer rate condition, the yield of loamy clay farmland was higher than that of sandy loam. The order of the annual average yields and WUE was N240 > N312 > N168, while the order of NUE of two textural soils was N168 > N240 > N312. This showed that excessive N fertilizer application would not increase maize yield, but lead to yield reduction. Therefore, the integrated index method was used to determine the BMPs for the two textural soils. The optimal N fertilizer rates in sandy loam and loamy clay soils were 180 and 200 kg ha−1, respectively, which not only ensured high crop yield, along with high WUE and NUE, but also reduced environmental N loss.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11233338/s1, Figure S1: Comparison of simulated and measured soil water content at different depths under N312 treatment in Sandy loam soil; Figure S2: Comparison of simulated and measured soil nitrate concentration at different depths under N312 treatment in Sandy loam soil; Figure S3: Comparison of measured and simulated leaf area index (LAI) under N312 treatment in Sandy loam soil; Figure S4: Comparison of measured and simulated total dry matter (TDM) under N312 treatment in Sandy loam soil; Figure S5: Comparison of simulated and measured soil water content at different depths under N240 treatment in loamy clay soil; Figure S6: Comparison of simulated and measured soil nitrate concentration at different depths under N240 treatment in loamy clay soil; Figure S7: Comparison of measured and simulated leaf area index (LAI) under N240 treatment in loamy clay soil; Figure S8: Comparison of measured and simulated total dry matter (TDM) under N240 treatment in loamy clay soil; Figure S9: Comparison of simulated and measured soil water content at different depths under N240 treatment in Sandy loam soil; Figure S10: Comparison of simulated and measured soil nitrate concentration at different depths under N240 treatment in Sandy loam soil; Figure S11: Comparison of measured and simulated leaf area index (LAI) under N240 treatment in Sandy loam soil; Figure S12: Comparison of measured and simulated total dry matter (TDM) under N240 treatment in Sandy loam soil; Figure S13: Comparison of simulated and measured soil water content at different depths under N168 treatment in loamy clay soil; Figure S14: Comparison of simulated and measured soil nitrate concentration at different depths under N168 treatment in loamy clay soil; Figure S15: Comparison of measured and simulated leaf area index (LAI) under N168 treatment in loamy clay soil; Figure S16: Comparison of measured and simulated total dry matter (TDM) under N168 treatment in loamy clay soil; Figure S17: Comparison of simulated and measured soil water content at different depths under N168 treatment in Sandy loam soil; Figure S18: Comparison of simulated and measured soil nitrate concentration at different depths under N168 treatment in Sandy loam soil; Figure S19: Comparison of measured and simulated leaf area index (LAI) under N168 treatment in Sandy loam soil; Figure S20: Comparison of measured and simulated total dry matter (TDM) under N168 treatment in Sandy loam soil; Figure S21: The relationship between measured and simulated soil water contents, nitrate concentrations, crop yields, and LAI under the validation treatments in Sandy loam soil; Table S1: Soil physical-chemical properties for soil profile in two experimental sites; Table S2: Analysis of significant different for yield and nitrogen uptake under different N rates for two textural soils in 2013 and 2014.

Author Contributions

Data curation, G.F.; methodology, validation, and writing original draft preparation, F.M.; writing—review and editing, K.H. and P.F.; project administration, Q.G. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the by the national natural science foundation of China (No. U19A2035) and the 2115 Talent Development Program of China Agricultural University (No. 1191-00109011).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO. 2020. Available online: http://www.fao.org/faostat/zh/#data (accessed on 19 December 2020).
  2. National Bureau of Statistics of China (NBSC). China Statistics Yearbook; China Statistics Press: Beijing, China, 2020. [Google Scholar]
  3. Liu, X.B.; Zhang, X.Y.; Wang, Y.X.; Sui, Y.Y.; Zhang, S.L.; Herbert, S.J.; Ding, G. Soil degradation: A problem threatening the sustainable development of agriculture in Northeast China. Plant. Soil. Environ. 2010, 56, 87–97. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, W.F.; Dou, Z.X.; He, P.; Ju, X.T.; Powlson, D.; Chadwick, D.; Norse, D.; Lu, Y.L.; Zhang, Y.; Wu, L.; et al. New technologies reduce greenhouse gas emissions from nitrogenous fertilizer in China. Proc. Natl. Acad. Sci. USA 2013, 110, 8375–8380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cui, Z.L.; Chen, X.P.; Zhang, F.S. Current nitrogen management status and measures to improve the intensive wheat-maize system in China. Ambio 2010, 39, 376–384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wick, K.; Heumesser, C.; Schmid, E. Ground water nitrate contamination: Factors and indicators. J. Environ. Manage. 2012, 111, 178–186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zhang, X.Y.; Chen, S.Y.; Sun, H.Y.; Shao, L.; Wang, Y. Changes in evapotranspiration over irrigated winter wheat and maize in North China Plain over three decades. Agr. Water. Manage. 2011, 98, 1097–1104. [Google Scholar] [CrossRef]
  8. Chen, L.F.; He, Z.B.; Zhao, W.Z. Soil structure and nutrient supply drive changes in soil microbial communities during conversion of virgin desert soil to irrigated cropland. Eur. J. Soil. Sci. 2020, 71, 768–781. [Google Scholar] [CrossRef]
  9. Zhang, X. A plan for efficient use of nitrogen fertilizers. Nature 2017, 543, 322–323. [Google Scholar] [CrossRef]
  10. Hu, K.L.; Li, Y.; Chen, W.P.; Chen, D.L.; Wei, Y.P.; Edis, R.; Li, B.G.; Huang, Y.F.; Zhang, Y.P. Modeling nitrate leaching and optimizing water and nitrogen management under irrigated maize in desert oases in Northwestern China. J. Environ. Qual. 2010, 39, 667–677. [Google Scholar] [CrossRef]
  11. Monzon, J.P.; Sadras, V.O.; Andrade, F.H. Fallow soil evaporation and water storage as affected by stubble in sub-humid (Argentina) and semi-arid (Australia) environments. Field Crops Res. 2006, 98, 83–90. [Google Scholar] [CrossRef]
  12. Cambouris, A.N.; Ziadi, N.; Perron, I.; Alotaibi, K.D.; Luce, M.S.; Tremblay, N. Corn yield components response to nitrogen fertilizer as a function of soil texture. Can. J. Soil. Sci. 2016, 96, 386–399. [Google Scholar] [CrossRef]
  13. Witt, C.; Ottesen, E.A. Lymphatic filariasis: An infection of childhood. Trop. Med. Int. Health. 2010, 6, 582–606. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, H.L.; Yang, J.Y.; He, P.; Bai, Y.L.; Jin, J.Y.; Drury, C.F.; Zhu, Y.P.; Yang, X.M.; Li, W.J.; Xie, J.G.; et al. Optimizing Parameters of CSM-CERES-Maize Model to Improve Simulation Performance of Maize Growth and Nitrogen Uptake in Northeast China. J. Integr. Agr. 2012, 11, 1898–1913. [Google Scholar] [CrossRef]
  15. Yang, J.M.; Yang, J.Y.; Dou, S.; Yang, X.M.; Hoogenboom, G. Simulating the effect of long-term fertilization on maize yield and soil C/N dynamics in northeastern China using DSSAT and CENTURY-based soil model. Nutr. Cycl. Agroecosys. 2013, 95, 287–303. [Google Scholar] [CrossRef]
  16. Li, Z.; Sun, Z. Optimized single irrigation can achieve high corn yield and water use efficiency in the Corn Belt of Northeast China. Eur. J. Agron. 2016, 75, 12–24. [Google Scholar] [CrossRef]
  17. Jiang, R.; He, W.; Zhou, W.; Hou, Y.; Yang, J.Y.; He, P. Exploring management strategies to improve maize yield and nitrogen use efficiency in northeast China using the DNDC and DSSAT models. Comput. Electro. Agr. 2019, 166, 104988. [Google Scholar] [CrossRef]
  18. Peng, Y.X.; Li, Z.Z.; Sun, T.; Zhang, F.; Wu, Q.; Du, M.; Sheng, T. Modeling long-term water use and economic returns to optimize alfalfa-corn rotation in the corn belt of Northeast China. Field Crops Res. 2022, 276, 108379. [Google Scholar] [CrossRef]
  19. Liang, H.; Hu, K.L.; Batchelor, W.D.; Qi, Z.M.; Li, B.G. An integrated soil-crop system model used for water and nitrogen management in North China. Sci. Rep. 2016, 6, 25755. [Google Scholar] [CrossRef] [Green Version]
  20. Liang, H.; Hu, K.L.; Qin, W.; Zuo, Q.; Guo, L.; Tao, Y.Y.; Lin, S. Ground cover rice production system reduces water consumption and nitrogen loss and increases water and nitrogen use efficiencies. Field Crops Res. 2019, 233, 70–79. [Google Scholar] [CrossRef]
  21. Shi, X.R.; Batchelor, W.D.; Liang, H.; Li, S.E.; Li, B.G.; Hu, K.L. Determining optimal water and nitrogen management under different initial soil mineral nitrogen levels in northwest China based on a model approach. Agr. Water. Manage. 2020, 234, 106110. [Google Scholar] [CrossRef]
  22. Xu, Q.; Hu, K.; Liang, H.; Leghari, S.J. Incorporating the WHCNS model to assess water and nitrogen footprint of alternative cropping systems for grain production in the North China Plain. J. Clean. Prod. 2020, 263, 121548. [Google Scholar] [CrossRef]
  23. Ma, L.W.; Ahuja, L.R.; Nolan, B.T.; Malone, R.W.; Trout, T.J.; Qi, Z. Root zone water quality model (RZWQM2): Model use, calibration, and validation. Trans. ASABE. 2012, 55, 1425–1446. [Google Scholar] [CrossRef]
  24. Van Liew, M.W.; Garbrecht, J. Hydrologic simulation of the Little Washita River experimental watershed using SWAT. J. Am. Water. Resour. As. 2003, 39, 413–426. [Google Scholar] [CrossRef]
  25. Schaap, M.G. Rosetta (1.0). 1999. Available online: http//.www.ussl.ars.usda.gov/MODELS/rosetta/rosetta.html (accessed on 19 December 2020).
  26. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration. Guidelines for computing crop water requirements. In FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; 300 D05109. [Google Scholar]
  27. Green, W.H.; Ampt, G.A. Studies on soil physics: Part 1. The flow of air and water through soils. J. Agr. Sci. 1911, 4, 1–24. [Google Scholar]
  28. Šimůnek, J.; Saito, H.; Sakai, M.; Genuchten, T.M. The HYDRUS-1D Software Package for Simulating the One-dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media; Department of Environmental Sciences University of California: Riverside, CA, USA, 2008. [Google Scholar]
  29. Hansen, S.; Abrahamsen, P.; Petersen, C.T.; Styczen, M.E. Daisy: Model use, calibration and validation. Trans. ASABE. 2012, 55, 1315–1333. [Google Scholar] [CrossRef]
  30. Driessen, P.M.; Konjin, N.T. Land Use System Analysis; Wageningen Agricultural University: Wageningen, The Netherlands, 1992. [Google Scholar]
  31. Li, Z.; Hu, K.; Li, B.; He, M.; Zhang, J. Evaluation of water and nitrogen use efficiencies in a double cropping system under different integrated management practices based on a model approach. Agr. Water. Manage. 2015, 159, 19–34. [Google Scholar] [CrossRef]
  32. Liang, H.; Hu, K.L.; Qin, W.; Zuo, Q.; Zhang, Y.A. Modelling the effect of mulching on soil heat transfer, water movement and crop growth for ground cover rice production system. Field Crop Res. 2017, 201, 97–107. [Google Scholar] [CrossRef]
  33. Liang, H.; Qi, Z.M.; Hu, K.L.; Prasher, S.O.; Zhang, Y. Can nitrate contaminated groundwater be remediated by optimizing flood irrigation rate with high nitrate water in a desert oasis using the WHCNS model. J. Environ. Manage. 2016, 181, 16–25. [Google Scholar] [CrossRef]
  34. Yang, J.M.; 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, 81e89. [Google Scholar] [CrossRef]
  35. Lu, J.; Bai, Z.H.; Velthof, G.L.; Wu, Z.G.; Chadwick, D.; Ma, L. Accumulation and leaching of nitrate in soils in wheat-maize production in China. Agric. Water. Manage. 2019, 212, 407–415. [Google Scholar] [CrossRef]
  36. Ziadi, N.; Cambouris, A.N.; Nyiraneza, J. Across a landscape, soil texture controls the optimum rate of N fertilizer for maize production. Field Crops Res. 2013, 148, 78–85. [Google Scholar] [CrossRef]
  37. Feng, G.Z.; Wang, Y.; Yan, L.; Mi, G.H.; Gao, Q. Effect of soil type and nitrogen application rate on yield of monocropping spring maize and nitrogen balance in crop field. Acta Pedol Sin 2017, 54, 444–455, (in Chinese with English abstract). [Google Scholar]
  38. Huang, J.; Hartemink, A.E.; Kucharik, C.J. Soil-dependent responses of US crop yields to climate variability and depth to groundwater. Agr. Syst. 2021, 190, 103085. [Google Scholar] [CrossRef]
  39. Huang, J.; Hartemink, A.E. Soil and environmental issues in sandy soils. Earth-Sci. Rev. 2020, 208, 103295–103316. [Google Scholar] [CrossRef]
  40. Wang, Y.; Guo, T.; Qi, L.; Zeng, H.; Liang, Y.; Wei, S.; Gao, F.; Wang, L.; Zhang, R.; Jia, Z. Meta-analysis of ridge-furrow cultivation effects on maize production and water use efficiency. Agric. Water. Manage. 2020, 234, 106144. [Google Scholar] [CrossRef]
  41. Gu, C.; Riley, W.J. Combined effects of short term rainfall patterns and soil texture on soil nitrogen cycling—A modeling analysis. J. Contam. Hydrol. 2010, 112, 141–154. [Google Scholar] [CrossRef] [Green Version]
  42. Su, Y.Z.; Yang, X.; Yang, R. Effect of soil texture in unsaturated zone on soil nitrate accumulation and groundwater nitrate contamination in a marginal oasis in the middle of Heihe River basin. Environ. Sci. 2014, 35, 3683–3690, (in Chinese with English abstract). [Google Scholar]
  43. Michalczyk, A.; Kersebaum, K.C.; Dauck, H.P.; Roelcke, M.; Yue, S.C.; Chen, X.-P.; Zhang, F.S. Quantifying nitrogen loss and water use via regionalization and multiple-year scenario simulations in the North China Plain. J. Plant. Nutr. Soil. Sci. 2020, 183, 718–733. [Google Scholar] [CrossRef]
  44. Milroy, S.P.; Asseng, S.; Poole, M.L. Systems analysis of wheat production on low water-holding soils in a Mediterranean-type environment II. Drainage and nitrate leaching. Field Crops Res. 2008, 107, 211–220. [Google Scholar] [CrossRef]
  45. Wang, L.S.; He, Z.B.; Zhao, W.Z.; Wang, C.; Ma, D.K. Fine Soil Texture Is Conducive to Crop Productivity and Nitrogen Retention in Irrigated Cropland in a Desert-Oasis Ecotone, Northwest China. Agronomy 2022, 12, 1509. [Google Scholar] [CrossRef]
  46. Li, J.S.; Yang, F.Y.; Li, Y.F. Water and nitrogen distribution under subsurface drip fertigation as affected by layered-textural soils. Trans. CSAE. 2009, 25, 25–31. [Google Scholar]
  47. Azad, N.; Behmanesh, J.; Rezaverdinejad, V.; Abbasi, F.; Navabian, M. An analysis of optimal fertigation implications in different soils on reducing environmental impacts of agricultural nitrate leaching. Sci. Rep. 2020, 10, 1–15. [Google Scholar] [CrossRef] [PubMed]
  48. Zhang, Y.T.; Wang, H.Y.; Liu, S.; Lei, Q.L.; Liu, J.Q.; He, J.Q.; Zhai, L.M.; Ren, T.Z.; Liu, H.B. Identifying critical nitrogen application rate for maize yield and nitrate leaching in a Haplic Luvisol soil using the DNDC model. Sci. Total. Environ. 2015, 514, 388–398. [Google Scholar] [CrossRef] [PubMed]
  49. Wang, Y.; Li, C.; Li, Y.; Zhu, L.; Liu, S.; Yan, L.; Feng, G.; Gao, Q. Agronomic and environmental benefits of nutrient expert on maize and rice in Northeast China. Environ. Sci. Pollut. R. 2020, 27, 28053–28065. [Google Scholar] [CrossRef] [PubMed]
  50. Awale, R.; Chatterjee, A. Enhanced Efficiency Nitrogen Products Influence Ammonia Volatilization and Nitrous Oxide Emission from Two Contrasting Soils. Agron. J. 2016, 109, 47–57. [Google Scholar] [CrossRef]
  51. Cabrera, M.L.; Kissel, D.E.; Rema, J.A.; Newsome, J.F.; Calvert, V.H., II. Ammonia volatilization from urea-based fertilizers applied to tall fescue pastures in Georgia, USA. Soil. Sci. Soc. Am. J. 2008, 72, 1665–1671. [Google Scholar]
  52. Jamali, H.; Quayle, W.; Scheer, C.; Rowlings, D.; Baldock, J. Effect of soil texture and wheat plants on N2O fluxes: A lysimeter study. Agr. Forest. Meteorol. 2016, 223, 17–29. [Google Scholar] [CrossRef] [Green Version]
  53. Gaillard, R.; Duval, B.D.; Osterholz, W.R.; Kucharik, C.J. Simulated effects of soil texture on nitrous oxide emission factors from corn and soybean agroecosystems in Wisconsin. J. Environ. Qual. 2016, 45, 1540–1548. [Google Scholar] [CrossRef]
  54. Pihlatie, M.; Syvasalo, E.; Simojoki, A.; Esala, M.; Regina, K. Contribution of nitrification and denitrification to N20 production in peat, clay and loamy sand soils under different soil moisture conditions. Nutr. Cycl. Agroecosys. 2004, 70, 135–141. [Google Scholar] [CrossRef]
  55. Abbasi, M.K.; Tahir, M.M.; Sadiq, A.; Iqbal, M.; Zafar, M. Yield and nitrogen use efficiency of rainfed maize response to splitting and nitrogen rates in Kashmir, Pakistan. Agron. J. 2012, 104, 448–457. [Google Scholar] [CrossRef]
  56. Cui, Z.L.; Wang, G.; Yue, S.; Wu, L.; Zhang, W.; Zhang, F.S.; Chen, X.P. Closing the N-use efficiency gap to achieve food and environmental security. Environ. Sci. Technol. 2014, 48, 5780–5787. [Google Scholar] [CrossRef]
  57. Cao, Q.; Miao, Y.X.; Feng, G.H.; Gao, X.; Liu, B.; Liu, Y.; Li, F.; Khosla, R.; Mulla, D.J.; Zhang, F. Improving nitrogen use efficiency with minimal environmental risks using an active canopy sensor in a wheat-maize cropping system. Field Crops Res. 2017, 214, 365–372. [Google Scholar] [CrossRef]
  58. Zheng, W.K.; Zhang, M.; Liu, Z.G.; Zhou, H.; Lu, H.; Zhang, W.; Yang, Y.; Li, C.; Chen, B. Combining controlled-release urea and normal urea to improve the nitrogen use efficiency and yield under wheat-maize double cropping system. Field Crops Res. 2016, 197, 52–62. [Google Scholar] [CrossRef]
  59. Li, Y.; White, R.; Chen, D.; Zhang, J.; Li, B.; Zhang, Y.; Huang, Y.; Edis, R. A spatially referenced water and nitrogen management model (WNMM) for (irrigated) intensive cropping systems in the North China Plain. Ecol. Model. 2007, 203, 395–423. [Google Scholar] [CrossRef]
  60. Chen, X.P.; Cui, Z.L.; Vitousek, P.M.; Cassman, K.G.; Matson, P.A.; Bai, J.-S.; Meng, Q.F.; Hou, P.; Yue, S.C.; Römheld, V.; et al. Integrated soil-crop system management for food security. Proc. Natl. Acad. Sci. USA 2011, 108, 6399–6404. [Google Scholar] [CrossRef] [Green Version]
  61. Chen, Y.L.; Zhang, J.; Li, Q.; He, X.L.; Su, X.P.; Chen, F.J.; Yuan, L.X.; Mi, G.H. Effects of nitrogen application on post-silking root senescence and yield of maize. Agron. J. 2015, 107, 835–842. [Google Scholar] [CrossRef]
  62. Wu, L.Q.; Wu, L.; Cui, Z.L.; Chen, X.P.; Zhang, F.S. Basic NPK fertilizer recommendation and fertilizer formula for maize production regions in China. Acta. Pedol. Sin. 2015, 52, 802–817, (in Chinese with English abstract). [Google Scholar]
Figure 1. Comparison of simulated and measured soil water content at different depths under N312 treatment in loamy clay soil.
Figure 1. Comparison of simulated and measured soil water content at different depths under N312 treatment in loamy clay soil.
Plants 11 03338 g001
Figure 2. Comparison of simulated and measured soil nitrate concentration at different depths under N312 treatment in loamy clay soil.
Figure 2. Comparison of simulated and measured soil nitrate concentration at different depths under N312 treatment in loamy clay soil.
Plants 11 03338 g002
Figure 3. Comparison of measured and simulated leaf area index (LAI) under N312 treatment in loamy clay soil.
Figure 3. Comparison of measured and simulated leaf area index (LAI) under N312 treatment in loamy clay soil.
Plants 11 03338 g003
Figure 4. Comparison of measured and simulated total dry matter (TDM) under N312 treatment in loamy clay soil.
Figure 4. Comparison of measured and simulated total dry matter (TDM) under N312 treatment in loamy clay soil.
Plants 11 03338 g004
Figure 5. The relationship between measured and simulated soil water contents, nitrate concentrations, crop yields, and LAI under the validation treatments in loamy clay soil.
Figure 5. The relationship between measured and simulated soil water contents, nitrate concentrations, crop yields, and LAI under the validation treatments in loamy clay soil.
Plants 11 03338 g005
Figure 6. Responses of simulated maize yield (a), WUE (b), nitrate leaching (c), ammonia volatilization (d), and NUE (e) to different N fertilizer rates for sandy loam and loamy clay soils in 2014.
Figure 6. Responses of simulated maize yield (a), WUE (b), nitrate leaching (c), ammonia volatilization (d), and NUE (e) to different N fertilizer rates for sandy loam and loamy clay soils in 2014.
Plants 11 03338 g006aPlants 11 03338 g006b
Table 1. Statistical criteria for simulation of total dry matter, yield, and crop N uptake.
Table 1. Statistical criteria for simulation of total dry matter, yield, and crop N uptake.
Soil TextureCalibrationYieldTotal Dry MatterCrop N Uptake
ValidationnRMSEEdnRMSEEdnRMSEEd
Sandy loamC14.260.420.7612.160.860.9826.560.560.86
V16.520.360.7119.620.920.9629.620.460.82
Loamy clayC12.180.440.796.610.880.9622.560.680.92
V10.560.380.7316.560.960.9828.640.580.86
Note: C, calibration; V, validation; nRMSE, normalized root mean square error; E, Nash–Sutcliffe modeling efficiency; d, index of agreement.
Table 2. Water balance and water use efficiency in 100 cm soil profiles under different N rates for two textural soils in 2013 and 2014.
Table 2. Water balance and water use efficiency in 100 cm soil profiles under different N rates for two textural soils in 2013 and 2014.
YearSoil TextureTreatmentPETDRWbalYieldWUE
(mm)(mm)(mm)(mm)(mm)(kg ha−1)(kg m−3)
2013Sandy loamN168615.1501.2151.422.5−60.089221.8
N240615.1503.7163.322.5−74.397671.9
N312615.1502.5155.122.5−65.098262.0
Loamy clayN168615.1400.5142.764.87.199472.5
N240615.1391.2138.364.820.811,4782.9
N312615.1391.2138.364.820.810,8962.8
2014Sandy loamN168420.1392.26.34.916.766771.7
N240420.1392.16.34.916.875411.9
N312420.1392.26.34.916.768081.7
Loamy clayN168420.1415.60.84.9−1.111,9232.9
N240420.1415.70.84.9−1.212,5763.0
N312420.1415.70.84.9−1.211,8222.8
MeanSandy loamN168517.6446.778.913.7−21.778001.8
N240517.6447.984.813.7−28.886541.9
N312517.6447.480.713.7−24.283171.9
MeanLoamy clayN168517.6408.171.834.93.010,9352.7
N240517.6403.569.634.99.812,0273.0
N312517.6403.569.634.99.811,3592.8
Note: P is precipitation; ET is evapotranspiration; D is drainage.; R is runoff; water balance, Wbal = P-ET-D-R; water use efficiency, WUE = measured yield/ET/10.3.3. N fate and N use efficiency under two textural soils.
Table 3. N budgets, N-use efficiencies, and the integrated index under different N rates for two textural soils from 2013 to 2014.
Table 3. N budgets, N-use efficiencies, and the integrated index under different N rates for two textural soils from 2013 to 2014.
SoilTreatNferNnetNvolNdenNupNleaBalanceWUENUEAFEFVCRII
(kg ha−1)(kg ha−1)(kg ha−1)(kg ha−1)(kg ha−1)(kg ha−1)(kg ha−1)(kg m−3)(kg kg−1)
2013SLN16816883.728.34.3140.29.968.91.848.80.90.923.61.0
N24024084.242.84.5182.419.175.41.939.31.01.218.10.9
N31231283.956.494.6180.525.9128.41.436.70.91.514.00.7
2013LCN168168109.614.512.4174.53.672.62.548.50.90.926.31.0
N240240108.321.911.9201.77.0105.82.747.31.01.221.30.9
N312312108.329.211.9204.79.7164.72.842.61.01.515.50.8
2014SLN16816816.221.81.7121.20.239.41.746.10.90.617.71.1
N24024016.235.32.7146.40.271.61.940.91.00.914.01.0
N31231216.249.24.1171.50.5103.02.530.20.91.59.70.7
2014LCN16816811.522.41.5213.30.1−57.72.650.31.00.928.91.1
N24024011.529.41.5228.80.1−8.43.048.41.01.220.70.9
N31231211.942.21.5249.40.230.62.840.30.91.515.40.7
Mean SLN16816850.025.13.0130.75.154.21.847.50.90.820.71.1
N24024050.239.13.6164.49.773.51.940.11.01.116.11.0
N31231250.152.84.4176.013.2115.72.033.50.91.511.90.7
Mean LCN16816860.618.57.0193.91.97.52.649.41.00.927.61.1
N24024059.925.76.7215.33.648.72.947.91.01.221.00.9
N31231260.135.76.7227.15.097.72.841.51.01.515.50.8
Note: SL, sandy loam; LC, loamy clay; Nfer, N fertilizer; Nnet, net mineralization; Nvol, ammonia volatilization; Nden, N denitrification; Nup, crop N uptake; Nlea, N-leaching; balance = Nfert + Nnet −Nvol – Nden – Nup − Nlea; NUE (N use efficiency) = grain yield/(Nvol+ Nden + Nup + Nlea); WUE, water use efficiency; AF, agronomy factor; EF, environment factor; VCR, value-to-cost ratio; II, Integrated index.
Table 4. Integrated index simulated by WHCNS model under different N rates for two soil textural soils in 2014.
Table 4. Integrated index simulated by WHCNS model under different N rates for two soil textural soils in 2014.
SoilNfert
(kg ha−1)
Nvol
(kg ha−1)
Nden
(kg ha−1)
Nlea
(kg ha−1)
Yield
(kg ha−1)
WUE
(kg m−3)
NUE
(kg kg−1)
VCRII
Sandy loam202.51.40.15200.310.311.60.3
16020.61.50.166801.831.218.31.0
18024.51.70.271801.832.017.31.1
20028.11.90.271801.830.016.01.0
40066.92.81.071801.824.78.40.7
Loamy clay201.71.80.03910.18.68.70.2
18019.31.50.093662.338.123.11.0
20022.41.50.0111642.639.924.21.1
22025.91.5 111642.735.822.61.0
40057.81.50.1111642.727.312.40.7
Note: Nfer, N fertilizer; Nvol, ammonia volatilization; Nden, N denitrification; Nlea, N leaching; VCR value-to-cost ratio; II, Integrated index.
Table 5. Soil physical and hydraulic properties for soil profiles in two experimental sites.
Table 5. Soil physical and hydraulic properties for soil profiles in two experimental sites.
SiteSoil Layer (cm)Particle Fraction (%)Texture (USDA)BD (g
cm−3)
θfc (cm3
cm−3)
θwp (cm3
cm−3)
θs (cm3
cm−3)
ks (cm
d−1)
SandSiltClay
Fujiajie0–2070.615.611.8Sandy loam1.500.220.070.27060.6
20–4082.29.78.5Loamy sand1.460.220.060.26370.4
40–6087.75.56.8Loamy sand1.490.230.070.27570.6
60–8083.911.76.4Loamy sand1.460.220.050.26980.0
80–10080.012.57.5Loamy sand1.460.260.040.30780.0
Sankeshu0–2032.824.143.1Loamy clay1.580.420.160.4632.24
20–4031.537.230.3Loamy clay1.580.420.180.4714.91
40–6033.030.936.1Loamy clay1.520.400.190.4574.52
60–8021.751.127.2Clay loam1.540.380.180.40622.7
80–10022.051.526.5Clay loam1.600.390.190.42620.9
Note: BD, bulk density; θfc, field capacity; θwp, soil wilting point; θs, soil saturated water content; ks, soil saturated hydraulic conductivity.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/plants11233338

AMA Style

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(23):3338. https://doi.org/10.3390/plants11233338

Chicago/Turabian Style

Meng, Fanchao, Kelin Hu, Puyu Feng, Guozhong Feng, and Qiang Gao. 2022. "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 11, no. 23: 3338. https://doi.org/10.3390/plants11233338

APA Style

Meng, F., Hu, K., Feng, P., Feng, G., & Gao, Q. (2022). 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, 11(23), 3338. https://doi.org/10.3390/plants11233338

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