Introduction

Apples, which are widely planted on the Loess Plateau, are the main income-generating fruit tree, and apple orchards cover a planted area of more than 1.2 million ha across the Weibei Dryland Highlands of Shaanxi1. A total of 80% of these apple orchards are rainfed (dryland), and their intense water consumption due to apple transpiration has caused excessive soil water depletion and consequent deep soil desiccation. Furthermore, the excessive and widespread applications of nitrogen and phosphorus in apple orchards improve apple yield and quality but intensify deep soil desiccation and restrict the stable and healthy development of the apple industry in semi-humid regions2,3,4. Excess applications of nitrogen and phosphorus was observed in 74% and 54% of the apple orchards, respectively, in the Weibei Dryland Highland5. Management of N and P fertilization is pivotal for sustainable fruit production in North China6. Wang et al.7 found that in the semi-arid region of the Loess Plateau, excessive application of P could result in high P-fixation capacity and consequently, severely damaged environments. In terms of economic returns, Qinguan apples have a P-use efficiency that is 50% higher than that of corn8. Researchers in Zagora Pelion (Central Greece) found that in apple orchards, the N inputs were five times higher than the nitrogen outputs; thus, reductions in N fertilization were considered9.

Some researchers have carried out short-term field experiments to determine the effects of different fertilization treatments on fruit yield and the quality of apple. However, few long-term experiments have investigated apple yield, soil water use and soil desiccation under different fertilization treatments4, 10, 11. Several models have been proposed to monitor apple tree growth12,13,14, e.g., models for root water uptake and soil water recovery after apple tree harvesting. The Environmental Policy Integrated Climate (EPIC) model has been used as a simulation model of crop productivity systems for many crops and is a universal model for water and soil management as well as crop productivity evaluation15,16,17. Ko et al.18 found EPIC to be a useful method for managing irrigated cotton and maize. Michele19 used EPIC to simulate irrigation scheduling in sunflower and found that the model effectively compared management strategies. The soil water and yield dynamics of grain crop rotation systems in the semi-humid region of the Loess Plateau have been well described using the EPIC20 model. Previous model applications involving soil water and yield of apple orchards on the Loess Plateau of China indicate that the EPIC model accurately simulates apple production21. However, further studies are needed to simulate soil water and crop productivity in response to different fertilization regimes applied to apple orchards on the Loess Plateau.

Using the EPIC model, this study sought to quantitatively simulate yield fluctuations and soil desiccation under different fertilization treatments applied to apple orchards over a long period of time (45 years). The effects of fertilization on apple yield and deep soil desiccation were analyzed to determine the best fertilization regime relative to local precipitation and to promote sustainable soil water utilization. Hence, the main objective of this study was to guide fertilization management decisions for sustainable apple production in the Weibei Dryland Highland on the Loess Plateau and in similar regions of the world.

Results

Simulated yields in dryland apple orchards

The annual precipitation in Luochuan between 1965 and 2009 ranged from 343.50 mm to 899.60 mm and averaged 596.64 mm, with a standard error of 117.17 mm and a coefficient of variation of 19.64%. The precipitation decreased in a fluctuating manner and was 44.53 mm lower in the last 15 years compared to the first 15 years.

During the fruit-bearing stage of the apple trees (apple trees start to bear fruit at 4 years old, as simulated), the yields of 4- to 45-year-old rainfed apple orchards in the F0-F5 treatments tended to first increase to their maxima and then decrease in a fluctuating manner from 1968 to 2009 (Fig. 1). The yield curves of the different fertilization treatments generally peaked in the initial simulated years and in the rainy years and then exhibited a trough during the latter simulated years and in the dry years. The yields of orchards subjected to the respective F0-F5 treatments averaged 24.42, 27.27, 28.69, 29.63, 30.49 and 29.43 t/ha, with standard deviations of 6.14, 6.11, 6.83, 7.68, 8.12 and 8.36 t/ha. Compared to yields under F0, the yields under F1-F5 increased by 11.70%, 17.49%, 21.35%, 24.85% and 20.54%, respectively, which indicated that F3 and F4 performed the best. After 42 fruiting years, the yields decreased considerably under all fertilization treatments. The greatest decrease occurred under F0, and greater decreases occurred under higher compared with lower fertilization applications. The lowest yields occurred under the 6 fertilization treatments in 1995 when the precipitation was only 343.5 mm, and the highest yields under F0 and F1 occurred in 1971; under the other fertilization treatments, the highest yields occurred in 1975 when the precipitation reached 899.6 mm.

Figure 1
figure 1

Simulated annual fruit yields under different fertilization treatments and precipitation in the apple orchards of Luochuan.

Simulated drought stresses in rainfed apple orchards

Drought stress is referred to as a restriction of crop growth due to water deficit, and in EPIC, its duration indicate its severity during the growing season. Over the 45 study years, the number of drought stress days under the respective F0-F5 treatments ranged from 0 to 57, 0 to 67, 0 to 69, 0 to 72, 0 to 78 and 0 to 83 days and averaged 16.71, 20.95, 24.01, 27.41, 33.06 and 35.12 days annually, with standard deviations of 14.84, 16.92, 17.41 18.24, 19.45 and 20.62 days. The number of drought stress days increased as the fertilization amount increased, and drought stress occurred earlier under higher fertilization. In the apple orchards, the earliest year of drought stress arising under the 6 respective fertilization treatments occurred in 1977 (in the 13-year-old orchards), 1974 (in the 10-year-old orchards), 1972 (in the 8-year-old orchards), 1972 (in the 8-year-old orchards), 1971 (in the 7-year-old orchards) and 1971 (in the 8-year-old orchards). The number of drought stress days under the different fertilization treatments tended to fluctuate in a basically identical manner, opposite from the manner in which the local precipitation varied. The peak number of drought stress days appeared in the dry years of 1986, 1995, 1997 and 2004. The number of drought stress days in the six different fertilization treatments peaked in 1995, reaching 56.39, 66.50, 68.83, 72.42, 78.40 and 83.14 days, in the F0, F1, F2, F3, F4 and F5, respectively, because precipitation in 1995 was only 343.5 mm (Fig. 2).

Figure 2
figure 2

Simulated number of water stress days under different fertilization treatments in the apple orchards of Luochuan.

Simulated monthly available soil moisture in the 0–15 m soil layer in the apple orchards

From 1965 to 2009, the simulated monthly available soil moisture in the 0–15 m soil layer in the respective F0-F5 treatments ranged from 723.22 to 2573.23, 699.88 to 2557.88, 681.28 to 2533.53, 669.60 to 2508.09, 663.72 to 2486.96 and 661.60 to 2479.27 mm and averaged 1267.40, 1223.44, 1191.33, 1153.71, 1116.45 and 1102.67 mm, with standard deviations of 556.9, 563.1, 555.2, 544.7, 525.3 and 521.1 mm and coefficients of variation of 43.9%, 46.0%, 46.6%, 47.2%, 47.0% and 47.2%. No obvious differences existed in available soil moisture between the high fertilization treatments (i.e., F4 and F5), but significant differences in available soil moisture occurred in the other treatments. As fertilization increased, the available soil moisture decreased and presented annual increases and seasonal fluctuations.

During the 45 simulated years, the monthly available soil moisture in the 0–15 m soil layer under the different fertilization treatments tended to vary similarly (Fig. 3). In the early stage of the 45-year simulation, the simulated monthly available soil moisture in the respective F0-F5 treatments tended to decrease clearly in a fluctuating manner before 1990, 1985, 1984, 1982, 1980 and 1980, with decreases at 1282, 1130, 1315, 1283, 1315 and 1345 mm and soil desiccation rates of 51.68, 58.57, 65.75, 75.47, 87.67 and 89.67 mm/a. Later, the monthly available soil moisture in the 0–15 m soil layer in the F0-F5 treatments varied at a lower level with seasonal precipitation. No obvious difference occurred in the soil desiccation rate between F4 and F5. In the Weibei Dryland Highlands, it frequently rains from July to September, resulting in a recovery of available soil moisture in most years.

Figure 3
figure 3

Simulated monthly available soil moisture in the 0–15 m soil level under different fertilization treatments in the apple orchards of Luochuan.

Simulated soil moisture distribution in the 0–15 m soil layer in apple orchards

The soil moisture distributions in apple orchards during the initial (1965–1971), middle (1980–1986) and final (2003–2009) simulated years are shown in Fig. 4. The figure shows that there were approximately the same simulated annual soil moisture distributions in the 0–15 m soil layer under F0-F5. Compared with the soil moisture in the apple rhizospheres when the simulation began, in 1965, the soil moisture in the apple rhizospheres decreased yearly as the apple trees grew and as the dry soil layers appeared and then deepened and thickened each year until they stabilized. The dry soil layers in the respective F0-F5 treatments occurred in 1976, 1974, 1973, 1973, 1972 and 1972 and took 25, 23, 21, 18, 15 and 15 years to reach a depth of 11 m. As fertilization increased, the dry soil layers occurred earlier, accompanied by accelerated formation, and the depth of 11 m was reached earlier. Soon after the dry soil layers became stable, the effects of fertilization on the soil moisture in the apple orchards disappeared; the soil water at maximum soil depth was utilized by the apple trees; therefore, the field water consumption of the apple trees mainly relied on seasonal precipitation. There was considerable variation in soil moisture in the 0–2 m soil layer, which resulted from the influence of rainfall infiltration.

Figure 4
figure 4

Simulated annual soil moisture profiles in the 0–15 m soil layer under different fertilization treatments in the apple orchards of Luochuan.

Simulated water consumption in apple orchards

During the period from 1965 to 2009, the simulated average annual water consumption in the respective F0-F5 treatments in the dryland apple orchards in Luochuan ranged from 335.8 to 831.0, 334.8 to 846.1, 332.6 to 864.7, 332.7 to 876.3, 332.9 to 899.3 and 332.9 to 900.0 mm and averaged 623.2, 625.6, 627.8, 628.0, 628.5 and 629.7 mm, with standard deviations of 110.4, 115.3, 118.5, 120.0, 122.5 and 123.4 mm. In apple orchards with 45-year-old trees, the total increased water consumption in the respective F1-F5 treatments was approximately 108, 207, 216, 239 and 293 mm, i.e., annual increases of 2.4, 4.6, 4.8, 5.3 and 6.5 mm compared with the average water consumption under F0. The average annual water consumption in the orchards increased as fertilization increased.

The water consumption differed significantly among the different fertilization treatments before the apple trees were 25 years old and was almost identical after the age of 26, although the values fluctuated with precipitation (Fig. 5). During the 45 years, the water use efficiencies in the F0-F5 treatments in the apple orchards of Luochuan were 39.19, 43.59, 45.70, 47.18, 48.51 and 46.74 t/(mm·hm2), respectively. These results indicate lower water use efficiency under F0 compared with the other treatments.

Figure 5
figure 5

Simulated annual water consumption in the 0–15 m soil layer under different fertilization treatments in the apple orchards of Luochuan.

In the orchards, the annual soil water over-consumption was the difference between the annual water consumption and precipitation, a measurement of the soil water balance between soil water consumption intensity and soil water replenishment capacity from seasonal rainfall. The soil water over-consumption fluctuated around approximately 0 in the later stage of the 45-year simulation (Fig. 6). The soil water over-consumption in F1-F5 was 2.5, 4.6, 4.8, 5.3 and 6.5 mm higher, respectively, than that in F0. It follows from the standpoint of sustainable soil water utilization that F2 and F3 are the best choices, although under F4 and F5, the yields and water use efficiencies were higher, but soil desiccation occurred faster and more severely.

Figure 6
figure 6

Simulated water over-consumption in the 0–15 m soil layer under different fertilization treatments in the apple orchards of Luochuan.

Economic returns from apple orchards

Depending on the simulated yields and corresponding fertilizer rates under the different fertilization treatments, a regression equation was established between the yields and fertilizer rates as follows:

$${\rm{y}}=-2E-05{x}^{2}+0.0218x+24.692\quad {R}^{2}=0.98$$
(1)

In this equation, x is the application rate of N, and y is the yield of the apple orchard in question. The prices of nitrogen, phosphorus and apples were calculated depending on the prices and composition of the fertilizers. Because other management practices that accompanied fertilization were the same, their input costs were not taken into account, such that the overall economic returns of one apple orchard were calculated using the following equation:

$$Y=4y-4.5{X}_{N}-7.5{X}_{P}$$
(2)

In this equation, Y is the total economic return from the apple orchard, y is the yield of the apple orchard, XN is the application rate of nitrogen, and XP is the application rate of phosphorus. The overall economic returns from the apple orchard under the different fertilization treatments were calculated using this equation (Table 1). For comparison, we obtained the maximum economic returns and corresponding fertilization for the apple orchards. In Luochuan, the maximum economic returns were 117,985 yuan/hm2 orchard, and the corresponding orchard yield and nitrogen and phosphorus application rates were 30.49 t/hm2, 480 kg/hm2 and P 240 kg/hm2, respectively (Table 1).

Table 1 Economic returns from apple orchards under different fertilization treatments in Luochuan.

A suitable fertilization rate based on different statistical indicators of dryland apple orchards in Luochuan can be sourced from a comprehensive analysis of orchard yields, water use efficiency, soil water over-consumption and economic benefits (Table 2). In terms of orchard yield increases, sustainable soil water utilization and economic returns, the most suitable fertilization rates in dryland apple orchards in Luochuan were observed under F3 and F4, i.e., N 360–480 kg/ha and P 180–240 kg/ha.

Table 2 Recommended fertilizer amount for apple orchards in Luochuan based on different statistical indictors.

Discussion

Dryland apple orchard yield and soil desiccation

The Weibei Dryland Highlands, a primary apple-growing region, has various meteorological conditions that are advantageous for apple growth and production. In this region, more than 80% of the orchards are rainfed—these orchards are not irrigated and are subjected to severe soil water over-consumption, marked by deep soil desiccation, soil impoverishment and low fertility. As a result, proper fertilization is very important for soil fertility maintenance, yield increases and sustainable soil water utilization22. Holb et al.23 reported that better mobility of artificial fertilizers gives rise to higher macronutrient uptake in apple orchards compared to organic orchards in Debrecen-Pallag, Hungary. Jun et al.24 found that nitrogen and phosphorus applied at rates of 250–300 kg/ha and 150–200 kg/ha, respectively, achieved an orchard target yield of 25–30 t/ha. In the current study, the best fertilization treatments were N 360–480 kg/ha and P 180–240 kg/ha, i.e., the nitrogen rate is at a level between its optimum and a practical level, and the phosphorus rate is equal to the optimum. This result likely occurred because in this study, the orchards were dryland orchards that were not irrigated and thus had low fertilizer utilization. Another reason for the large amounts of applied nitrogen was that no apple orchards were ever measured. Nitrogen can increase apple yields, but its excessive application will result in a series of adverse results, such as a deterioration in apple quality, extensive growth of tree components, decreased soil pH, soil nitrate accumulation and groundwater pollution25, 26.

Our results add to previous findings, i.e., that water and nutrients are the most important factors that influence the yields of dryland orchards, and fertilization can effectively improve water use efficiency and yields of dryland orchards but can easily increase water consumption through transpiration, resulting in increased yields accompanied by aggravated soil water consumption27. In our study, under high fertilization levels, the available soil moisture stabilized earlier and soil desiccation occurred faster; furthermore, no obvious difference in the soil desiccation rate was observed between F4 and F5. The available soil moisture markedly decreased in the middle stage of the 45-year simulation, and the dry years somewhat recovered and increased during years with high precipitation. In Luochuan, the soil water storage obviously decreased. The precipitation led to fluctuating decreases in the average annual water consumption because of the decreases and fluctuations in annual precipitation, in addition to the influence of increasing soil desiccation; and this was one of the main reasons for the fluctuating decrease in the apple yields. We found that fertilization increases can increase water use efficiency but may result in soil desiccation and soil water over-consumption in deep soil layers. Proper fertilization can improve fruit quality, reduce production costs, maintain soil fertility, eliminate nutrient deficiencies, and maintain tree vigor28. In addition to fertilization, the growth and yield of “Cortland” apples also depended on the weather conditions and tree age29.

This study demonstrated that the yields under F1-F5 increased by 11.70%, 17.49%, 21.35%, 24.85% and 20.54%, respectively, compared with the average yield under F0 in the 45 years of study. The yields increased obviously as fertilization increased, but F5, which had excessive fertilization, resulted in the wasting of resources and had an adverse influence on apple yield. Furthermore, soil desiccation occurred faster, and soil drought stress occurred earlier and was more severe with higher levels of fertilization. The thickness of the dry soil layers differed obviously among the different fertilization treatments before the apple trees reached 23 years of age. The dry soil exceeded a depth of 11 m in all treatments except the control, and soil moisture was close to the wilting point over a long-term period. In the apple orchards, the dry soil layers formed earlier and deepened faster with high fertilization compared to low fertilization. During the period of time when soil water decreased obviously, the soil desiccation rates in the fertilization treatments of F0, F1, F2, F3, F4 and F5 were 51.68, 58.57, 65.75, 75.47, 87.67 and 89.67 mm/a. Higher fertilization readily resulted in deep soil desiccation, which led to a fluctuation in soil water storage under lower level precipitation over a long-term period.

EPIC simulation in this study

Recently, Li et al. developed a simulation verification and application of EPIC for the Loess Plateau, revealing dynamic water productivities of artificial forestlands and grain crops in dryland regions, providing new ideas for studying water productivity as well as the effects of soil water and fertilizer in apple orchards30, 31. Research on water productivity and soil water has mainly concentrated on short-term field experiments. Consequently, there is a lack of experiments that are continuously monitored over a long period of time. By combining a short-term field experiment with quantitative and dynamic simulations over a long-term period using WinEPIC, this study generated continuous dynamic data that could not be obtained using conventional research methods. In addition, the evolution of annual water productivities and deep soil water dynamics under different fertilization treatments in apple orchards was examined, thus providing a scientific basis for determining suitable fertilization treatments and growth years relative to local precipitation. This was a simple and effective quantitative research method. However, simulation differences occurred depending on the factors in different planting areas, the density of the apple trees, and the fertilization treatment. There were some differences between the simulated and practical growth conditions, which resulted from the omission of the effects of different site conditions, e.g., slope and aspect, in the simulation. In addition, we were unable to simulate the effects of orchard management measures on apple yields and soil water conservation, flower and fruit thinning, tree pruning and thinning, orchard grass planting and film mulching, etc. These defects might lead to differences between the simulated and measured results, although the simulation may still reflect water productivity and soil water variations under different fertilization treatments and may provide a scientific basis for soil water management and suitable fertilization. The simulation parameters must be further modified to more accurately simulate water productivity and the effects of deep soil desiccation.

Conclusions

During the simulated period from 1965 to 2009, the apple yields under F0-F5 tended to first increase and then decrease in a fluctuating manner; as the fertilization level increased, drought stress occurred earlier and was more intensive. The monthly available soil moisture in all of the fertilization treatments presented an obvious fluctuating decrease and then varied at a lower level with seasonal precipitation. As fertilization increased, drying of the soil layers occurred earlier and deepened more rapidly; and the the soil dry layers reached their maximum depths earlier, and the average annual water consumption, soil water over-consumption and water use efficiencies increased. In terms of yield increases, sustainable soil water utilization and orchard economic returns, F3 and F4 were the best choices. Therefore, a suitable fertilization rate for dryland apple orchards in Luochuan is N 360–480 kg/ha and P 180–240 kg/ha.

Materials and Methods

Study location

The simulation study was conducted in apple orchards in Luochuan, a county in Shaanxi, which is located in the central area of the apple production region on the Loess Plateau and in the northeastern part of the Weibei Dryland Highland (109°13′14″—109°45′47″E, 35°26′29″—36°04′12″N). With a warm temperate and semi-humid continental monsoon climate, the county has an average altitude of 1072 m, an average annual precipitation of 600 mm, an average annual temperature of 9.2 °C, and an average frost-free period of 167 d; the rainy season coincides with high temperatures. Luochuan has a broken topography typical of the high plain and gully loess regions of the Loess Plateau, which is characterized by uniformly textured Heilou clay soil and very thick loessial soil that can exceed a depth of 100 m. In Luochuan, most apple orchards are rainfed, i.e., they are not irrigated.

WinEPIC Profile

The Environmental Policy Integrated Climate (EPIC) model includes weather simulation, hydrology, eroded sediments, nutrient cycling, pesticide rate, crop growth, soil temperature, soil tillage, economic benefit and crop environmental control. Based on integrated growth parameters of more than 120 types of field crops, pastures, and forests, as well as apple trees, the EPIC model can be employed to simulate long-term dynamic changes in soil water and nutrient utilization as well as crop productivity on a daily basis. Furthermore, EPIC has been used to evaluate agro-ecological system management and its effect on soil and water resources32, 33. WinEPIC v. 3060, which was adopted in this study, is a new EPIC model that can be run on Windows with a user-friendly interface. Using WinEPIC, soil water dynamics can be described in more detail; the daily soil moisture that is thus obtained for different soil layers can be used to study crop production systems and, in particular, can be used to simulate ecological environmental effects related to soil moisture in drylands34. Our research team has simulated the biomass and soil desiccation effects of artificial Robinia pseudoacacia and Chinese pine forestlands on the Loess Plateau and has preliminarily investigated the yield and soil water dynamics in apple orchards on the Weibei Dryland Highlands35,36,37. The details of mathematical methods used for crop productivity as well as soil water and fertilizer movement have been reviewed in related literature38,39,40,41,42.

Dataset construction for WinEPIC

WinEPIC was run with datasets comprised of local daily weather variables, soil physical and chemical properties, crop growth parameters and the management measures of apple orchards. Adapted from real-time weather data observations (1965~2009) from the Luochuan weather station, the daily weather variables mainly included daily solar radiation (MJ/m2), maximum temperature (°C), minimum temperature (°C), precipitation (mm), relative humidity (%), and wind speed (m/s). These daily weather data were input into WinEPIC in the EPIC format; then, a daily weather database for Luochuan was constructed43. The soil physical and chemical properties, which consisted of more than 40 terms relating to 0–15 m deep Heilu clay soil in Luochuan (Table 3), were adapted from soil census data such as those contained in the Records of the Chinese Soil Survey, the Soil of Shaanxi44, 45. These data were transformed into a typical soil parameter database of the research region in EPIC. On the Loess Plateau, soil water more than 10 m deep is typically used by rainfed orchards, and in some rainfed orchards, water that is nearly 15 m deep can be used. The maximum soil depth in this study was set to 15 m to examine the deep soil water utilization of the apple trees in the orchards. There was a relatively uniform soil texture at a depth of 0 to 15 m, unlike the other soil layers. To facilitate the analysis of soil moisture dynamics, the soil profile in question was divided into 17 soil layers: 0–0.01 m, 0.01–0.5, 0.5–1 m and 1-m-thick layers that were distributed between 1 m and 15 m deep. Some important physiological and ecological parameters of apples were modified in EPIC (Table 4) in reference to related research and observed data from the Loess Plateau46.

Table 3 Some important physical and chemical properties of Heilu clay soil in Luochuan.
Table 4 Some important modified parameters of apple growth used in WinEPIC.

Simulation methods

In this study, 6 different fertilization treatments were applied to the local apple orchards on the same dates: F0 (N 0 kg/ha, P 0 kg/ha), F1 (N 120 kg/ha, P 60 kg/ha), F2 (N 240 kg/ha, P 120 kg/ha), F3 (N 360 kg/ha, P 180 kg/ha), F4 (N 480 kg/ha, P 240 kg/ha) and F5 (N 600 kg/ha, P 300 kg/ha). The soil water simulations commenced in the years when the apple trees were transplanted. In this simulation, the fertilizers were applied each year at a given time according to the fertilization treatment protocol. After the soil water parameters were entered in the model, the initial soil water condition was automatically set at 75% of the field water-holding capacity. Using WinEPIC, the daily dynamics of growth, soil water and fertilizer utilization from 1965 to 2009 were quantitatively simulated, and the daily biomass growth and soil water and nutrient balances were obtained. To fully reflect soil water consumption and the soil water balance conditions resulting from precipitation, the following was done: within the daily data sequence of simulated soil moisture in the 0–15 m layer, the data from the 15th of each month were considered as characteristic of the available soil moisture for the month, and the soil moisture profile distribution data from the 1st of November of each year were considered as characteristic for the year.

WinEPIC model verification

To verify the EPIC performance in simulating the yield and soil moisture of the dryland apple orchards, these values were simulated for apple orchards in Luochuan from 1980 to 2009. The differences between the simulated and measured yield and soil moisture (2–10 m) were compared to verify the performance of the EPIC model. The measured yields were adapted from orchard surveys conducted at different ages by the experimental station in Luochuan. A significant correlation was observed between the simulated and measured fruit yield in the apple orchards of Luochuan, with a relative deviation of 0.76%, an RMSE of 1.09 t/ha, and a correlation index of 0.94. With respect to soil water, the relative error was 13.1%, the RMSE was 0.019 m/m, and the correlation index was 0.97 (Table 5). In addition, the slope of the linear regression between the simulated and observed yield and soil water was close to 1, with an R2 value that showed a significant correlation (Fig. 7). EPIC revised with these parameters could be used to accurately simulate the yield response and soil moisture characteristics of dryland apple orchards in Weibei.

Table 5 Simulated and measured yield and soil moisture of apple orchards in Luochun.
Figure 7
figure 7

Simulated and measured yield (a) and soil moisture (b) in the rainfed apple orchards of Luochuan.