Next Article in Journal
Greater Nitrogen Availability, Nitrous Oxide Emissions, and Vegetable Yields with Fall-Applied Chicken Relative to Horse Manure
Next Article in Special Issue
Assessing the Impact of Best Management Practices in a Highly Anthropogenic and Ungauged Watershed Using the SWAT Model: A Case Study in the El Beal Watershed (Southeast Spain)
Previous Article in Journal
Distributions of Alien Invasive Weeds under Climate Change Scenarios in Mountainous Bhutan
Previous Article in Special Issue
AquaCropR: Crop Growth Model for R
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 9
Issue 8
10.3390/agronomy9080443
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain

by
Jie Zhao
1,2,†,
Qingquan Chu
1,†,
Mengjie Shang
1,
Manyowa N. Meki
2,
Nicole Norelli
2,
Yao Jiang
1,
Yadong Yang
1,
Huadong Zang
1,
Zhaohai Zeng
1,* and
Jaehak Jeong
2,*
1
College of Agronomy and Biotechnology, China Agricultural University, No. 2 Yuanmingyuan West Road, Beijing 100193, China
2
Blackland Research Center, Texas A&M AgriLife Research, Texas A&M University, 720 East Blackland Road, Temple, TX 76502, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Agronomy 2019, 9(8), 443; https://doi.org/10.3390/agronomy9080443
Submission received: 18 June 2019 / Revised: 26 July 2019 / Accepted: 6 August 2019 / Published: 10 August 2019
(This article belongs to the Special Issue Model Application for Sustainable Agricultural Water)
Browse Figures
Review Reports Versions Notes

Abstract

:
Spring peanut is a valuable alternative crop to mitigate water scarcity caused by excessive water use in conventional cropping systems in the North China Plain (NCP). In the present study, we evaluated the capability of the Agricultural Policy Environmental eXtender (APEX) model to predict spring peanut response to sowing dates and seeding rates in order to optimize sowing dates, seeding rates, and irrigation regimes. Data used for calibration and validation of the model included leaf area index (LAI), aboveground biomass (ABIOM), and pod yield data collected from a field experiment of nine sowing dates and seeding rate combinations conducted from 2017 to 2018. The calibrated model was then used to simulate peanut yield responses to extended sowing dates (5 April to 4 June with a 5-day interval) and seeding rates (15 plants m−2 to 50 plants m−2 with a 5 plants m−2 interval) using 38 years of weather data as well as yield, evapotranspiration (ET), and water stress days under different irrigation regimes (rainfed, one irrigation before planting (60 mm) or at flowering (60 mm), and two irrigation with one time before planting and one time at flowering (60 mm each time) or at pod set (60 mm each time)). Results show that the model satisfactorily simulates pod yield of peanut based on R2 = 0.70, index of agreement (d value) being 0.80 and percent bias (PBIAS) values ≤4%. Moreover, the model performed reasonably well in predicting the emergence, LAI and ABIOM, with a R2 = 0.86, d = 0.95 and PBIAS = 8% for LAI and R2 = 0.90, d = 0.97 and PBIAS = 1% for ABIOM, respectively. Simulation results indicate that the best combination of sowing dates and seeding rates is a density of 35–40 plants m−2 and dates during early-May to mid-May due to the influence of local climate and canopy structure to the growth and yield of peanut. Under the optimal sowing date and plant density, an irrigation depth of 60 mm during flowering gave a pod yield (5.6 t ha−1) and ET (464 mm), which resulted in the highest water use efficiency (12.1 kg ha−1 mm−1). The APEX model is capable of assessing the effects of management practices on the growth and yield of peanut. Sowing 35–40 plants m−2 during early-May to mid-May with 60 mm irrigation depth is the recommended agronomic practice for peanut production in the water-constrained NCP.

1. Introduction

The increasing need for high-quality food with the minimum environmental impact has led to renewed interest in grain legumes to maintain sustainable crop production [1,2,3]. Peanut (Arachis hypogaea L.) plays a critical role in oil security as well as being an important cash crop in China. China is the largest producer of peanut in the world [4]. The North China Plain (NCP), which is often referred to as the breadbasket of China, provides more than 72% of the nation’s wheat and 33% of its maize production [5]. Such achievements have heavily depended on continuous overexploitation of groundwater for irrigation and excessive mineral nitrogen (N) fertilizer application to meet the water and N requirement of intensive wheat–maize cropping systems [6,7]. Adjustment and optimization of cropping systems could be an effective strategy to mitigate groundwater depletion and reduce N input to ensure sustainable food production [8,9,10,11]. Alternative cropping systems include spring maize monoculture [12,13], a 2-year system of winter wheat/summer maize–spring maize [11] or double maize cropping systems [14] to allow replenishment of soil water reserves and fertility. Peanut has significant beneficial effects on subsequent crops in rotation with its inherent capacity for symbiotic atmospheric nitrogen fixation and considerable economic return, making it a valuable alternative crop in cropping structure adjustment [15,16].
Under the policy of “fallow in winter wheat season and rainfed in rainy season”, which aimed to reduce the winter wheat planting area [12,17], the cultivated area of spring peanut (or early sowing peanut, planting during early April to late May) increased steadily. Contrast to the fixed sowing dates (5 June to 20 June, planting after winter wheat) of summer peanut, spring peanut can be planted over a considerable part of the season under suitable conditions in the NCP. However, the choosing of a sowing date must have consideration of alleviating drought or high-temperature stress during critical stages of peanut growth [18]. Sowing dates for peanut have been studied extensively in most peanut-growing countries [19,20,21,22]. These studies indicate that matching the phenology of the peanut crop to the duration of favorable conditions by selecting the most suitable sowing dates to avoid periods of stress is crucial for maximizing peanut yield. Moreover, appropriate planting density is needed to fully explore the peanut’s genetic potential for better yield. Planting density not only determines competition for light, water, and nutrients, but also regulates vegetative growth and reproductive development and ultimately controls the distribution of dry matter between the organs and pod yield [23,24,25]. Peanut seed is also more expensive than the seed of most row crops, making it important to reduce seed requirement and costs. In addition, supplementary irrigation during the critical period of water requirement is crucial to increase the productivity of peanut [26,27]. Despite these resources, research on optimum sowing dates, seeding rates, and irrigation regimes of spring peanut is limited particularly in the NCP.
Experimental approaches aimed at optimizing management practices are costly and time-consuming because of the necessity in capturing the potential impact of long-term climate. Opportunities exist to examine the effect of management practices through process-based crop models. The Agricultural Policy/Environmental eXtender (APEX) is a process-based model for managing landscapes, whole farms, and small watersheds to achieve sustainable production and maintain environmental quality [28]. APEX is capable of simulating crop growth, hydrology, soil erosion, soil carbon sequestration, nutrient cycling, and losses [29]. Due to its wide applicability and strong capability, the APEX model has been extensively used to simulate grain yield, irrigation regimes, carbon, and nitrogen management under a variety of management conditions across the world. Wang et al. (2008) successfully predicted maize grain yield and soil organic carbon in an erosion-prone area in Iowa, US [30]. Cavero et al. (2012) assessed irrigation and N management strategies using APEX in three Mediterranean watersheds [31]. Meki et al. (2011; 2013) simulated maize residue removal effects on nutrient, soil organic carbon losses, and subsequent soybean yield and nitrogen dynamics in the Upper Mississippi River Basin [32,33]. Soybean yield potential, yield gap under rainfed conditions, and crop-water production functions were identified in a humid region of Mississippi, USA [34,35,36]. Luo and Wang (2019) evaluated corn, rice, and wheat grain yields and sediments under various irrigation scenarios and tillage methods in southwest China [37]. Plotkin et al. (2013) used the APEX model to simulate peanut yield, hydrology, and pesticide losses for a peanut/cotton rotation under the conditions of the Southern Atlantic Coastal Plain [38]. Van Liew et al. (2017) evaluated crop yields and sediments under peanuts and cotton dominated cropping systems on the Little River watersheds in Georgia [39]. Overall, these studies indicate that the APEX model is capable of successfully modeling the impacts of climate, soil, cropping system, tillage, irrigation, and management practices on crop growth and productivity.
Although simulation models have the potential for optimizing crop management, they must be calibrated and validated before they can be applied in the NCP. We conducted a field experiment to compare peanut growth and yield response to seasons, sowing dates, and seeding rates. Data from this experiment was used to calibrate and validate the APEX model. The objectives of this study are to (1) evaluate the performance of APEX on estimating peanut growth and crop yield, (2) assess peanut yield under extended scenarios of sowing dates and seeding rates, and (3) identify the best practice for maximizing water productivity under different irrigation regime scenarios. It is of great significance to improve peanut water productivity and further the application of the APEX model under water constrained areas in the NCP.

2. Materials and Methods

2.1. Experimental Site Description

The field experiment was conducted at the Wuqiao Experimental Station of the China Agricultural University in Wuqiao County (37°37′ N, 116°26′ E, altitude 18–21 m) of Hebei Province, China (Figure 1). The climate of the study area is a warm temperate and continental monsoon climate, typically with hot, rainy summers and dry, cold winters. Monthly trends in rainfall and average air temperature during the last 38 years including 2017 and 2018, the experiment years, are shown in Figure 2. The annual mean air temperature over the last 38 years was 12.9 °C, the annual accumulated temperature (≥0 °C) was 4826 °C, the frost-free period was 201 days, and the annual number of sunshine hours was 2724 h. The average annual precipitation over the last 38 years was 554 mm, 70% of which occurred from June to August (Figure 2). The predominant soil is classified as Calcaric Fluvisol with a silt loam texture. Table 1 lists the soil texture and hydraulic properties for each layer of the 1.2 m profile. The station is representative of the overall agricultural production and climate conditions of the NCP. The winter wheat–summer maize rotation is the main cropping system and irrigated cropland accounts for about 70% of the total land in this area. Groundwater is the major source for irrigation as almost no surface water is available. The groundwater table is approximately 17 m below ground surface.

2.2. Experimental Design and Data Collection

Measured data for evaluating the APEX model’s performance in simulating peanut growth and yield were obtained from an experiment conducted during the growing seasons of 2017 and 2018. A split-plot design with three replicates was used each year for data collection. The main plot had three sowing dates, i.e., late April (25 April, T1), early May (5 May, T2), and mid-May (15 May, T3); and sub-plots had three levels of seeding rates, 24 (D1), 30 (D2), and 36 (D3) plants m−2. Individual sub-plots were built with 7.5 cm high, 40 cm wide peanut ridging beds with two rows, the inner rows were spaced 30 cm apart and outer rows were spaced 81 cm for D1, 59 cm for D2, and 44 cm for D3, respectively (Figure 3). Plant spacing was 15 cm for all treatments and each plot was 5 m by 7 m in size. Each plot had 6, 7, and 9 raised beds for D1, D2, and D3, respectively. The rows were oriented north–south. Yuhua 25, a high yield median-large size peanut cultivar, planted with two seeds per hill, was selected from among 9 peanut varieties used in a previous experiment [40]. Flood irrigation (60 mm, recorded with totalizing water meters) was applied to all treatments to facilitate seed germination. All plots were fertilized prior to planting with 171 kg N ha−1 in the form of urea (225 kg ha−1, 46% N) and diammonium phosphate (375 kg ha−1, 18% N), 173 kg P2O5 ha−1 in the form of diammonium phosphate (375 kg ha−1, 46% P2O5), and 150 kg K2O ha−1 in the form of potassium sulfate (300 kg ha−1, 50% K2O). The herbicide Pendimethalin was applied at a rate of 0.99 kg ha−1 after sowing. The insecticides Chlorantraniliprole and Bacillus thuringiensis were applied as needed at a rate of 0.03 kg ha−1 at pegging stage and 600 mL ha−1 (8 IU mL−1) at early podding stage to control grubs (Phyllophaga spp.) and aphids (Aphidoidea), respectively. The fungicides Carbendazim and Mancozeb were applied both at an early podding stage at a rate of 0.68 kg ha−1 and 1.13 kg ha−1 to control peanut leafspot (Cercospora arachidicola), respectively. Although measures were implemented in both experimental years to avoid diseases and pests damage, there was still some plant damage due to the inefficiency of pesticides caused by frequent and heavy rainfall events during the 2018 growing season.
Data were collected on plant growth and development. The vegetative and reproductive growth stages including emergence (VE, cotyledons near the soil surface with some part of the plant visible in half of the seedlings.) and harvest maturity (R8, plants with 70%–80% mature pods) were determined using the system developed for peanut by Boote (1982) [41]. These stages were assumed to have occurred when at least 50% of the plants in a plot had reached that stage. Maximum plant height was measured at the R8 growth stage [41]. Growth analysis data were collected four times during the peanut growing season (Table 2). Measured plant variables included dry weight and leaf area index. Five consecutive plants in the central rows of each plot were sampled each time and washed to remove soils attached to the pods. Plant components were separated to determine leaf area and to partition dry matter into plant components such as leaves, stems, pegs, pods, and seeds. After recording the leaf area of each sample with a leaf area meter (Yaxin-1241, Beijing Yaxin Science Instrument Technology Co., Ltd., Beijing, China). The samples were oven-dried at 65 °C for 48 h and weighed to determine dry weight. At each harvest date, the peanut pods were dug up by hand with conventional tools. The commercial yield was assessed by harvesting 3-m lengths of row by hand from the central four rows, and drying until pods reached about 12% moisture content. Then, 500 g of pods as representative subsamples were hand-shelled to determine kernel yields and shelling percentages.

2.3. APEX Model Description and Input Parameters

APEX 1501 (last updated on 15 January 2019) was used in this study to develop APEX models for the peanut field trials. It has a user-friendly Windows interface, which is designed and maintained by the Blackland Research and Extension Center Modeling Team [28]. APEX was developed to expand EPIC’s capabilities to simulate multiple subareas simultaneously [42], which means it can be run for single fields similar to EPIC [43]. APEX uses EPIC’s general crop growth model that is capable of simulating various types of over 100 annual and perennial crops [44]. Each crop is defined by 56 crop parameters with unique values. APEX crop parameters and their values used for peanut simulation in this study are listed in Table 3. The APEX model is a continuous model which operates on a daily time step and can simulate up to hundreds of years [30,45]. Crop growth is simulated based on daily heat unit accumulation [44]. Potential increase in biomass is calculated daily based on intercepted photosynthetic active radiation and radiation use efficiency, and is then adjusted to actual growth depending on the level of daily stresses caused by water, nutrients, temperature, and aeration. If the crop is under the full stress condition on a day, the crop does not gain biomass for the day. In APEX, crop yield is estimated by using the harvest index, which is defined as the ratio of economic yield to aboveground biomass at maturity [46]. Water or temperature stress may cause a reduction in harvest index if they occur during sensitive stages such as flowering or grain-filling periods.
All the agricultural production processes are accounted for in the APEX model in a very detailed fashion [29]. APEX inputs required for model construction include weather data, soil properties, crop data, and management information. APEX simulation is mainly driven by weather data includes daily solar radiation, precipitation, maximum air temperature, and minimum air temperature. In addition, measured daily wind speed and relative humidity were added to the daily weather input to select the Penman or Penman–Monteith method for estimating potential evapotranspiration [46,47]. Weather data were added into the existing APEX weather database by WeatherImport, a stand-alone weather generator. Soil properties include layer depth, bulk density, field capacity, permanent wilting point, percentage sand, percentage silt, percentage organic carbon, and saturated hydraulic conductivity. Crop data includes growth, yield, leaf area index (LAI), harvest index, root depth, crop height, biomass, and plant population. Crop management inputs include the rate and/or date of plowing, sowing, fertilizer, irrigation, pest control, cultivar, and harvest. We used the harvesting plant density, around 90% of designed density in both years, as the input density in the operation file. As the current version cannot simulate twin rows in one raised bed, one row in one bed was adopted in the model with the average row interval of 0.56, 0.44, and 0.37m for D1, D2, and D3, respectively. The simulated peanut seed yield was adjusted to pod yield based on 5% of water content in seeds and 71% of shelling percentage for uniformity of comparisons.

2.4. APEX Model Calibration and Validation

Before a model is applied to simulate various scenarios, model calibration is often needed to obtain reliable results [48]. In this study, APEX models for the peanut trials were calibrated and validated using a field measured pod yield, above-ground plant biomass (ABIOM), and LAI. APEX was calibrated for the 2017 growing season and then validated using data collected in 2018, with a warm-up simulation of 6 years (2011–2016) to stabilize the soil pools (e.g., carbon and nitrogen). Calibration was carried out manually using a trial-and-error method in which sensitive influential parameters and inputs were adjusted within their recommended ranges. Model outputs were analyzed graphically by comparing simulated versus measured LAI, ABIOM, and yield data, in addition to simultaneously assessing model performance statistics [49]. Parameter sensitivity assessed during manual adjustments indicated that peanut biomass and yield were responsive to several input parameters, including crop biomass energy ratio (WA), harvest index (HI), the maximum LAI (DMLA), and the fraction of growing season when LAI starts to decline (DLAI) and, as expected, soil properties, weather, and management practices.

2.5. Simulation Scenarios

To optimize seeding rates and sowing dates, APEX scenarios that combine 8 seeding rates (15–50 plants m−2, 5 plants m−2 interval) and 13 sowing dates (5 April to 4 June, 5 days interval) were developed. In total, 104 scenarios under rainfed or full irrigation conditions were constructed and assessed for 38 years (1981 to 2018). These scenarios were made sufficiently broad to evaluate the effects of irrigation regimes on peanut pod yield and water use efficiency under the best seeding rates (40 plants m−2) and sowing dates (10 May). Five irrigation regime scenarios were designed, namely W0 (rainfed), W11 (60 mm depth before planting), W12 (60 mm depth at flowering), W21 (60 mm depth before planting and 60 mm depth at flowering), and W22 (60 mm depth before planting and 60 mm depth at pod set). In addition, in order to explore crop responses to irrigation in different rainfall category years, simulation years (1981 to 2018) were grouped into three precipitation regimes according to the rainfall amount during each growing seasons using an empirical frequency analysis [50,51]: Wet years (25% guaranteed rate of precipitation), normal years (50% guaranteed rate of precipitation), and dry years (75% guaranteed rate of precipitation). Model simulation results were evaluated for yield responses as well as water use efficiency (WUE). WUE (kg ha−1 mm−1) for the pod yield (kg ha−1) was calculated as follows [52]:
WUE = P Y E T
where PY is the pod yield (kg ha−1) and ET (mm) was simulated by the APEX model using Penman–Monteith methods. Water stress (WS), defined as actual plant uptake/potential plant water use rate, was used to identify the timing and duration of water deficit under the three precipitation regimes.

2.6. Data Analysis and Model Performance Criteria

Data analysis of the field experiments was conducted using the R statistical package [53,54,55]. A linear mixed-effects model with sowing dates (date) and seeding rates (density) as a fixed effect, and block and year as a random effect was fit by a reduced mixed model (no significant interaction between date and density):
Yield   ~   date   +   density   +   ( 1 | block )   +   ( 1 | year ) .
Analysis of variance (ANOVA) and multiple comparisons were conducted to examine the treatment differences.
To evaluate model performance and accuracy in prediction, statistical indicators that included the coefficient of determination (R2), root mean square error (RMSE), normalized root mean square error (NRMSE), and the Willmott (1981) index of agreement (d) were computed from observed and simulated variables (LAI, ABIOM, and pod yield) [56]. The RMSE is the standard deviation of residuals (i.e., the differences between observed and predicted values), showing how concentrated the predicted values are around the best fit. RMSE is calculated using Equation (3).
RMSE = ∑ i = 1 n ( P i − O i ) 2 n
where n is the number of observations, Pi and Oi refer to predicted and observed values, respectively. NRMSE is calculated using Equation (4).
NRMSE = R M S E M
where M is the mean of the observed variable. The index of agreement (d) is computed using Equation (5).
d = 1 − [ ∑ i = 1 n ( P i − O i ) 2 ∑ i = 1 n ( | P ′ i | + | O ′ i | ) 2 ]
where P’i = Pi − M and O’i = Oi − M. The d value is a better indicator of model performance, particularly relative to 1:1 line and values closer to 1 indicate better agreement between the two variables while a d value of zero indicates no predictability. Percent bias (PBIAS) measures the average tendency of the simulated data to be greater or less than the measured data. PBIAS is calculated by Equation (6).
PBIAS = [ ∑ i = 1 n ( O i − P i ) ∑ i = 1 n O i ] × 100
The optimal value of PBIAS is zero. A value of low magnitude indicates accurate model simulation. Positive values indicate model underestimation bias, and negative values indicate model overestimation bias. In addition, the model performance was assessed using the means and standard deviations of measured and predicted values.

3. Results

3.1. Sowing Date and Plant Density Effects on Yield

ANOVA results show that sowing dates had no significant effect (p = 0.08) on pod yield, but seeding rates had a significant effect (p < 0.001). Multiple comparison results show that D2 and D3 significantly increased the pod yield by 11.0% and 13.3% compared to D1, but no difference was observed between D2 and D3, although the pod yield under D3 is 2.1% higher than that under D2.

3.2. Model Calibration

Table 3 presents calibrated crop parameters in APEX, which were used to conduct further validation. Default values were kept for WA, HI, and DMLA because of the good fit between simulated and measured biomass, seed yield, and LAI. To increase the response of LAI to sowing density, plant population for crops and grass-1st point on curve (PPLP1) and plant population for crops and grass-2nd point on curve (PPLP2), the first and second of two points on the plant population curve were adjusted to 10.55 and 40.95, respectively. Then, the fraction of growing season when leaf area declines (DLAI), first and second point on optimal leaf area development curve (DLAP1 and DLAP2) were adjusted to fit the LAI curve. Leaf area index decline rate parameter (RLAD) and biomass–energy ratio decline rate parameter (RBMD) were adjusted to fit the ABIOM curve. Potential heat units (PHUs) for peanut were adjusted until the heat unit index (HUI, the fraction of total base-zero heat units at which operation takes place) at harvest closed to 1.0. The PHUs were set at 2000 to 2050 degree-days for different sowing dates. These above steps were repeated iteratively until reasonably accurate results were achieved.
APEX uses soil moisture content and soil temperature to estimate the timing of crop emergence. The simulated emergence was 14, 11, and 9 days for T1, T2, and T3 after sowing, as compared to the measured values of 13, 10, and 9 days (Table 2), indicating a good predictive capability. Simulated values for pod yield, LAI, and ABIOM were in good accordance with those measured in the calibration dataset (Table 4). Simulated pod yield was close to measured pod yield for calibration (4.96 ± 0.13 t ha−1 versus 4.96 ± 0.32 t ha−1, Figure 4a). The performance statistics achieved in this study on pod yield (R2 = 0.71, PBIAS = 0.04%) compare well with the criterion for judging satisfactory model performance (R2 ≥ 0.60, PBIAS < 25%) suggested by Wang et al. (2012) [49]. The model performed very well in LAI projection with a high value of R2 (0.95) and d (0.98) and low PBIAS (−1.67%) and NRMSE (15.33%). The trend line of simulated LAI conformed well to peanut growth, suggesting the calibration result for LAI growth curve was acceptable (Figure 5). APEX performed well with ABIOM simulation during early growth stages though errors increased after establishing full canopy (Figure 6). This may be caused by the out-of-sync of leaf senescence and defoliation. Performance statistics for ABIOM indicated the model performance was satisfactory (R2 = 0.83, d = 0.95, PBIAS = 8.53%, and NRMSE = 25.66%).

3.3. Model Validation

The APEX model for peanut was validated using data from the same treatments for the following year (2018) after calibration. Unlike the 2017 crop year, pest damage was reported during the 2018 crop year at these field trials. Thus, the APEX setting was modified to allow for simulating pest damage. Pest (insects and disease) factor (PST) was set to 0.85 based on field investigation results. The calibrated parameters along with the pest damage simulation yielded a good agreement between measured and simulated pod yields (4.66 ± 0.4 t ha−1 versus 4.82 ± 0.18 t ha−1, Figure 4b) with R2 of 0.72, PBIAS of −3.47% and d of 0.82 (Table 4). APEX performed better on LAI during the calibration period than the validation (Figure 7). The model performed well in ABIOM projection with a high value of R2 (0.90) and d (0.97) and a relatively low PBIAS (−1.21%) and NRMSE (18.11%) (Table 4). However, the model overestimated ABIOM during the late portion of the growing season in all T3 (mid-May sowing) scenarios (Figure 8g–i).

3.4. Optimum Seeding Rate and Sowing Rate

Pod yield and the maximum LAI under different sowing dates with rainfed irrigation were, on average, 22.8% and 15.4% lower than that with the full irrigation (Figure 9). The low pod yield (3.64 t ha−1) and maximum LAI (2.73) at 15 plants m−2 sharply increased in scenarios which have higher seeding rates up to 40 plants m−2 were evaluated. Thereafter, the pod yield and maximum LAI stagnated under rainfed conditions (Figure 9a). A stronger response of pod yield and maximum LAI to seeding rates were found under full irrigation scenarios. Pod yield of peanut increased rapidly from 4.67 t ha−1 at 15 plants m−2 to 5.90 t ha−1 at 40 plants m−2. Similarly, maximum LAI increased rapidly from 3.17 at 15 plants m−2 to 4.3 at 40 plants m−2. Thereafter, no increase in pod yield or maximum LAI was observed under an increased seeding rate with full irrigation (Figure 9b). A greater variability was observed in pod yield and maximum LAI of all sowing dates under rainfed irrigation than full irrigation scenarios (Figure 10). The smallest pod yield (3.73 t ha−1) of peanut sown under rainfed on 5 April sharply increased with the postponement of planting time to attain a peak (4.67 t ha−1) on 20 May (Figure 10a). From 25 May onward, pod yields started to decrease. In contrast, with the delay of planting time, maximum LAI of peanut under no irrigation firstly decreased and then increased, and minimized LAI (3.32) was observed on 5 May (Figure 10a). Pod yields of peanut under full irrigation increased firstly and then decreased as sowing date was delayed further into late May. The highest pod yield was 5.87 t ha−1 when peanut was sown on 10 May (Figure 10b). From 15 May onward, pod yields declined sharply as sowing dates were further delayed. Maximum LAI of peanut under full irrigation gradually declined as the planting date was delayed, ranging from 4.16 on 5 April to 3.83 on 4 June (Figure 10b).

3.5. Irrigation Scenarios, Simulation, and Analysis

Total precipitation for dry, normal, and wet years was 355 mm, 425 mm, and 510 mm, respectively. Thirty-eight simulation years were identified and categorized as 9 dry years, 19 normal years, and 10 wet years. In general, peanut pod yield, water consumption, and water use efficiency were enhanced and WS was decreased with the increase of irrigation volume and the optimum irrigation time (Figure 11). Effects of five irrigation regimes on pod yield under no other stress in different precipitation years are shown in Figure 11a. The simulated pod yields broadly ranged from 2.06 t ha−1 in a dry year (250 mm) under rainfed condition to 7.03 t ha−1 in a normal year (443 mm) under W21 treatment. The simulated peanut pod yield with rainfed irrigation was on average 4.80 t ha−1 across the 38 years. The irrigation volume of 60 mm applied at sowing or flowering increased the pod yield by 14%, compared with the rainfed peanut. However, the yield increases declined to 7% when the irrigation volume increased from 60 mm to 120 mm. The simulated pod yield for peanut with W12 was on average 9% higher compared with that with W11 in dry years, but no yield increases in normal and wet years. The same trend was found in peanut yield with W12 and W22. Greater yield increases were observed when supplemental irrigation was applied in dry years (32%) than that in normal (16%) and wet years (14%) when compared to that with no irrigation.
The simulated water consumption of peanut with different irrigation regimes during 1981–2018 ranged from 310 mm in a dry year (332 mm) under rainfed conditions to 593 mm in a wet year (622 mm) under W21 treatment (Figure 11b). The water consumption of peanut was positively proportional to the irrigation rate. Similar to the case of pod yield, increased water consumption for peanut with full irrigation was observed in dry years (19%) than in normal (11%) or wet years (10%) when compared with no irrigation scenarios. The simulated water use efficiency of peanut with different irrigation regimes ranged from 6.4 kg ha−1 mm−1 in a dry year (250 mm) under rainfed to 15.4 kg ha−1 mm−1 in a normal year (425 mm) under W12 (Figure 11c). The water use efficiency of peanut increased with an increase in irrigation rate in dry years, but not in normal or wet years. Specifically, the simulated WUE was on average 5% higher with W11, and 16% higher with W12, than rainfed peanut in dry years. In that case, any further increase in the irrigation water supply will have no effect. As expected, the simulated pod yield of peanut was responsive to the WS during the growing season. The simulated WS of peanut with different irrigation regimes ranging from 65 d in a dry year (332 mm) under rainfed conditions to 0 d in wet years with irrigation (Figure 11d). The mean WS among five irrigation regimes was 24 d in dry category years, which accounted for about 19% of total peanut growing days. The normal and the wet category years had WS values of 11 d and 7 d, respectively. The WS of peanut decreased from 23 d to 7 d with an increase in irrigation rate in all precipitation categories. Impact of irrigation on WS was greater for the dry year category (24 d) than that for the normal (11 d) or the wet years (7 d) when compared with no irrigation. Overall, under the optimal sowing dates and plant density, irrigation of 60 mm during flowering (W21) led to acceptable pod yield (5.6 t ha−1) and ET (464 mm), which resulted in the highest water use efficiency (12.1 kg ha−1 mm−1).

4. Discussion

4.1. Model Calibration and Validation

The days after planting (DAP) to the emergence of peanut was positively associated with the mean air temperature [57]. We found that the DAP to emergence reduced from 13 d for an early sowing date (T1) to 7 d for a late date (T3) with the rise of air temperature. The APEX model reasonably replicated the trend of emergence date in all three sowing dates. The model successfully replicated pod yields in both calibration and validation periods as evidenced by R2 values being greater than 0.7, though the pod yield was slightly underestimated during the calibration phase and slightly overestimated during the validation phase (Table 4). The low slope of yield during calibration and validation periods could be due to the moderate response of yield to seeding rate in APEX. Crop parameters PPLP1 and PPLP2 determine the response of maximum LAI to plant population and further influence dry matter production and yield formation. However, the yield did not change much when we adjusted these two parameters to fit the LAI curve. These results indicated that the yield increases slowly as the seeding rate increases. This means that the response of simulated yield is rather moderate. In contrast, increases in seeding rates can lead to significant increases in yield under field conditions. Consequently, the model ended up overestimating low yield values and underestimating high values. The differences between predicted and observed pod yield may also be driven by potential incidences of diseases and pests. This is especially true for the 2018 growing season, which had more precipitation and higher air temperature (Figure 3). Treatments with late sowing date and high plant density had greater incidences of disease and pests under a hotter and more humid environment, resulting in more yield losses. Therefore, the model tends to overestimate treatments with late sowing date and high plant density and underestimate the rest of the treatments, corresponding to the high and low yield points, respectively (Figure 4b). The performance of APEX model on yield prediction was very good when compared with good simulation (PBIAS = 12.6%) of peanut rotated with cotton under conventionally tilled or strip-tilled conditions [38] and satisfactory yield prediction (underestimate peanut yields by 20.4%) on the Little River watersheds in Georgia [39]. APEX simulations of LAI and ABIOM under 9 treatments were fairly accurate for the two growing seasons. The model predicted the ABIOM with good accuracy until 85 and 95 DAP during 2017 and 2018, respectively. After that, however, the model consistently overpredicted the ABIOM accumulation (Figure 6 and Figure 8). This occurred as the leaf dry-weight began to decline due to disease effects and defoliation around that time for which the model did not simulate. APEX would perform better if the simulation of disease-induced defoliation is added to the crop simulation processes. A more sophisticated pest damage function than the current one in APEX can improve the simulation accuracy of peanut growth and yield as well [58,59,60,61]. The late leaf spot function incorporated into CSM-CROPGRO-Peanut model by Singh et al. (2013) is a good example for improving the accuracy of peanut growth and development simulation [61].

4.2. Optimum Plant Density and Sowing Date

Increased seeding rate had a positive effect on pod yield. A satisfactory yield level can be achieved with 35–40 plants m−2 across 13 sowing dates. This may be explained by the increased LAI found in the optimum seeding rate inducing the better distribution of light and greater light interception in the canopy as well as contributing to higher absorption of nutrients from soil [20,24]. The yield was unaffected by further increasing planting density and with less than 35 plants m−2 it was significantly reduced. These negative effects on pod yield could be due to greater competition for light and water among plants at high plant densities or the lower efficiency of resource use at extremely low plant density [23,62]. Consistent with previous observations [25,63,64], these results indicate that optimum plant density greatly improves the canopy structure and light interception, which enhance nutrient and water use efficiency, and finally help produce the highest pod yield.
Mid-May sowings under rainfed irrigation gave the highest pod yield across 38 years and 8 planting densities, despite the lowest maximum LAI, and yield declined progressively as sowing occurred later in the year until early June despite increasing maximum LAI. Irrigation improved the pod yield and maximum LAI across 13 sowing dates, and impacted yield and LAI response to sowing dates. The optimum sowing date under full irrigation appeared 10 days earlier than that under rainfed due to the mitigation of water deficit by irrigation. Our results suggest that the most suitable sowing period falls between early to mid-May. It would be a disadvantage to the yields if sowed earlier to later. Previous studies indicate that the optimum air temperature for vegetative growth of peanut is between 30 °C and 35 °C and the ideal temperature for reproductive growth ranges between 25 °C and 28 °C. High temperatures above 35 °C during the reproductive stage reduce dry matter production, yield component, and pod yield [65,66,67]. Very early sowing before 30 April did not generate any advantage for yield in NCP due to water stress for vegetative growth under rainfed conditions and high-temperature stress during the pod-filling stage. In addition, a lower temperature during the early growth stage delayed flower initiation and extended maturity [21]. Therefore, it is inappropriate to sow early under the local climate and irrigation water availability. In relation to late sowing, a reduction in solar radiation and temperature during pod-filling led to a decrease in 100-seed weight and pod yield. These findings are well supported by Haro et al. (2007) who reported that delayed sowing causes a reduction in radiation use efficiency and pod yield of peanut [68]. This negative effect can be reduced by adjusting plant populations [19]. Overall, an optimum sowing date allows peanut to be exposed to appropriate temperature regimes and receive more solar radiation as well as to receive natural rainfall required for growth during the vegetative and the reproductive growth stages, enabling the crop to fulfill its yield potential.

4.3. Optimum Irrigation Regime

The yield reduction under rainfed ranged from 30% in dry years to 16% in wet years due to serious water stress compared to the full-irrigation (W21). There are two potential reasons for this result. In the first instance, water deficit reduced flowering and the formation of pods due to increased soil strength promoted by dry soil surface, resulting in a significant decline in pod numbers [26,69]. Secondly, water deficit inhibited canopy expansion and reduced radiation use efficiency, leading to reduced biomass production [26]. Although the ET reduction under rainfed ranged from 20% in dry years to 11% in wet years compared to the full-irrigated condition (W21), the water use efficiency of peanut did not increase. Instead, the WUE of peanut under rainfed was 14% lower than that under W21. These results indicate that rainfed peanut helped save water by not irrigating, but resulted in lower yield and WUE. Therefore, it is not recommended to grow peanut without irrigation during whole growing seasons in the NCP.
Pod yields of peanut under W12, which only irrigated 60 mm during flowering, were 29%, 15%, and 12% higher than that under rainfed in dry, normal, and wet years, respectively. Yet, pod yields of peanut under W12 were close to that under the W21 and W22 in all precipitation years. These results are well supported by previous research that the flowering and pod filling stages are the most sensitive phenological stages to soil water deficit [70,71]. Our findings also show that supplementary irrigation during the flowering stage is more efficient than that during pod filling stage. Moreover, ET reduction under W12 ranged from 10% in dry years to 6% in wet years compared to that under well-irrigated condition (W21) and resulted in higher WUE. Based on the combined effects of pod yield, ET, and WUE, peanuts can be irrigated 60 mm during the flowering stage (W12) under water shortage conditions.

4.4. Limitations

The current version of APEX is not able to simulate a twin-row planting pattern in one raised bed. Thus, we simplified it to a single row in the model simulation. These two kinds of planting patterns may have the same [62] or different [72] impacts on pod yield of peanut due to differences in canopy structure and light interception. Such a simplified representation of the planting pattern may influence the simulation accuracy of the APEX model.

4.5. Future Research

Accounting for pest damage is crucial for models to provide useful predictions of crop production under real conditions. Although yield losses due to pest and disease damage were simulated by the model in 2018 growing season given that a fraction of the yield remained after damage, the simulated aboveground biomass was, in general, less accurate in the late-season than early to mid-season. APEX was limited in simulating the reduction of live biomass by leaf defoliation. Therefore, an enhancement is needed in the pest damage module to better represent biomass losses associated with leaf defoliation. In addition, the peanut is a legume that has considerable economic return and N-fixation as well as the forage value. Further research is needed to investigate the possibility of replacing the conventional system of winter wheat and summer maize, with a 2-year system of winter wheat/summer maize-spring peanut at the cropping system level.

5. Conclusions

Spring peanut is a promising alternative crop to mitigate water scarcity and environmental pollution caused by excessive water use and N inputs in conventional cropping systems in the NCP. In the present study, we evaluated the capability of the APEX model in predicting the responses of peanut to sowing dates and seeding rate in order to optimize sowing dates, seeding rate, and irrigation regimes. The results indicate that the model performed satisfactorily in simulating pod yield of peanut. The model performed reasonably well in predicting the emergence, LAI, and aboveground biomass during the majority of the growing season, with the exception of the late portion of the growing season. In the model applications for optimizing the sowing dates and seeding rates, the plant density of 35–40 plants m−2 and sowing dates between early May (5 May) to mid-May (15 May) are the best management practices for achieving maximum yield due to the optimized climatic environment and canopy structure in respect to the growth and yield of peanut. Under optimal sowing date and plant density, a 60 mm irrigation during the flowering period led to acceptable pod yield (5.6 t ha−1) and ET (464 mm), which resulted in the highest water use efficiency (12.1 kg ha−1 mm−1). Overall, the sowing density of 35–40 plants m−2 during early May (5 May) to mid-May (15 May) with 60 mm of irrigation water can be recommended for optimum peanut performance in water constrained areas of the NCP.

Author Contributions

Conceptualization, J.Z., Z.Z., Q.C., and Y.Y.; investigation, Y.J., and M.S.; writing—original draft preparation, J.Z.; writing—review and editing, J.J., M.N.M., N.N., and H.Z.; visualization, J.Z.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2016YFD0300205-01 and the National Natural Science Foundation of China, grant number 31671640.

Acknowledgments

We are grateful for the constructive comments on this manuscript from the anonymous reviewers and editors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Calles, T. The International Year of Pulses: What are they and why are they important. Agric. Dev. 2016, 26, 40–42. [Google Scholar]
  2. Foyer, C.H.; Lam, H.-M.; Nguyen, H.T.; Siddique, K.H.M.; Varshney, R.K.; Colmer, T.D.; Cowling, W.; Bramley, H.; Mori, T.A.; Hodgson, J.M.; et al. Neglecting legumes has compromised human health and sustainable food production. Nat. Plants 2016, 2, 16112. [Google Scholar] [CrossRef] [PubMed]
  3. Zander, P.; Amjath-Babu, T.S.; Preissel, S.; Reckling, M.; Bues, A.; Schläfke, N.; Kuhlman, T.; Bachinger, J.; Uthes, S.; Stoddard, F.; et al. Grain legume decline and potential recovery in European agriculture: A review. Agron. Sustain. Dev. 2016, 36, 26. [Google Scholar] [CrossRef]
  4. FAOSTAT: Crop. Available online: http://www.fao.org/faostat/zh/#data/QC (accessed on 10 April 2019).
  5. NBS. China’s National Bureau of Statistics. 2018. Available online: http://www.stats.gov.cn/tjsj/ndsj/2018/indexeh.htm (accessed on 10 April 2019).
  6. Meng, Q.F.; Sun, Q.P.; Chen, X.P.; Cui, Z.L.; Yue, S.C.; Zhang, F.S.; Roemheld, V. Alternative cropping systems for sustainable water and nitrogen use in the North China Plain. Agric. Ecosyst. Environ. 2012, 146, 93–102. [Google Scholar] [CrossRef]
  7. Van Oort, P.A.J.; Wang, G.; Vos, J.; Meinke, H.; Li, B.G.; Huang, J.K.; van der Werf, W. Towards groundwater neutral cropping systems in the Alluvial Fans of the North China Plain. Agric. Water Manag. 2016, 165, 131–140. [Google Scholar] [CrossRef]
  8. An, P.L.; Ren, W.; Liu, X.L.; Song, M.M.; Li, X.M. Adjustment and optimization of the cropping systems under water constraint. Sustainability 2016, 8, 1207. [Google Scholar] [CrossRef]
  9. Davis, K.F.; Rulli, M.C.; Seveso, A.; D’Odorico, P. Increased food production and reduced water use through optimized crop distribution. Nat. Geosci. 2017, 10, 919–924. [Google Scholar] [CrossRef]
  10. Davis, K.F.; Chiarelli, D.D.; Rulli, M.C.; Chhatre, A.; Richter, B.; Singh, D.; DeFries, R. Alternative cereals can improve water use and nutrient supply in India. Sci. Adv. 2018, 4, eaao1108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Zhao, Z.G.; Qin, X.; Wang, Z.M.; Wang, E.L. Performance of different cropping systems across precipitation gradient in North China Plain. Agric. For. Meteorol. 2018, 259, 162–172. [Google Scholar] [CrossRef]
  12. Liang, H.; Qin, W.; Hu, K.L.; Tao, H.B.; Li, B.G. Modelling groundwater level dynamics under different cropping systems and developing groundwater neutral systems in the North China Plain. Agric. Water Manag. 2019, 213, 732–741. [Google Scholar] [CrossRef]
  13. Yan, P.; Tao, Z.Q.; Chen, Y.Q.; Zhang, X.P.; Sui, P. Spring maize kernel number and assimilate supply responses to high-temperature stress under field conditions. Agron. J. 2017, 109, 1433–1442. [Google Scholar] [CrossRef]
  14. Meng, Q.F.; Wang, H.F.; Yan, P.; Pan, J.X.; Lu, D.J.; Cui, Z.L.; Zhang, F.S.; Chen, X.P. Designing a new cropping system for high productivity and sustainable water usage under climate change. Sci. Rep. 2017, 7, 41587. [Google Scholar] [CrossRef] [PubMed]
  15. Zeng, Z.H.; Lu, Z.Y.; Jiang, Y.; Zhang, K.; Yang, Y.D.; Zhao, P.Y. Legume-cereal crop rotation systems in China. In Crop Rotations: Farming Practices, Monitoring and Environmental Benefits, 1st ed.; Ma, B., Ed.; Nova Science Publishers: New York, NY, USA, 2016; pp. 51–69. [Google Scholar]
  16. Zhang, K.; Zhao, J.; Wang, X.Q.; Zang, H.D.; Liu, J.N.; Hu, Y.G.; Zeng, Z.H. Estimates on nitrogen uptake in the subsequent wheat by above-ground and root residue and rhizodeposition of using peanut labeled with 15N isotope on the North China Plain. J. Integr. Agric. 2019, 18, 571–579. [Google Scholar] [CrossRef]
  17. Lu, H.; Xie, H.L.; Lv, T.G.; Yao, G.R. Determinants of cultivated land recuperation in ecologically damaged areas in China. Land Use Policy 2019, 81, 160–166. [Google Scholar] [CrossRef]
  18. Dangthaisong, P.; Banterng, P.; Jogloy, S.; Vorasoot, N.; Patanothai, A.; Hoogenboom, G. Evaluation of the CSM-CROPGRO-Peanut model in simulating responses of two peanut cultivars to different moisture regimes. Asian J. Plant Sci. 2006, 5, 913–922. [Google Scholar]
  19. Bell, M. Effect of sowing date on growth and development of irrigated peanuts, Arachis hypogaea L. cv. Early Bunch, in a monsoonal tropical environment. Aust. J. Agric. Res. 1986, 37, 361–373. [Google Scholar] [CrossRef]
  20. Gardner, F.; Auma, E. Canopy structure, light interception, and yield and market quality of peanut genotypes as influenced by planting pattern and planting date. Field Crop Res. 1989, 20, 13–29. [Google Scholar] [CrossRef]
  21. Banterng, P.; Patanothai, A.; Pannangpetch, K.; Jogloy, S.; Hoogenboom, G. Seasonal variation in the dynamic growth and development traits of peanut lines. J. Agric. Sci.-Camb. 2003, 141, 51–62. [Google Scholar] [CrossRef]
  22. Caliskan, S.; Caliskan, M.; Arslan, M.; Arioglu, H. Effects of sowing date and growth duration on growth and yield of groundnut in a Mediterranean-type environment in Turkey. Field Crop. Res. 2008, 105, 131–140. [Google Scholar] [CrossRef]
  23. Singh, P.; Boote, K.; Virmani, S. Evaluation of the groundnut model PNUTGRO for crop response to plant population and row spacing. Field Crop. Res. 1994, 39, 163–170. [Google Scholar] [CrossRef] [Green Version]
  24. Rasekh, H.; Asghari, J.; Massoumi, S.; Zakerinejad, R. Effect of planting pattern and plant density on physiological characteristics and yield of peanut (Arachis hypogaea L.) in Iran. Res. J. Biol. Sci. 2010, 5, 542–547. [Google Scholar] [CrossRef]
  25. Zhao, C.X.; Shao, C.L.; Yang, Z.L.; Wang, Y.F.; Zhang, X.J.; Wang, M.L.; McGiffen, M.E., Jr. Effects of planting density on pod development and yield of peanuts under the pattern of precision planted peanuts. Legume Res. 2017, 40, 901–905. [Google Scholar] [CrossRef]
  26. Haro, R.J.; Dardanelli, J.L.; Otegui, M.E.; Collino, D.J. Seed yield determination of peanut crops under water deficit: Soil strength effects on pod set, the source–sink ratio and radiation use efficiency. Field Crop. Res. 2008, 109, 24–33. [Google Scholar] [CrossRef]
  27. Amiri, E.; Abdzad Gohari, A.; Mianabadi, A. Evaluation of water schemes for peanut, using CSM-CROPGRO-Peanut model. Arch. Agron. Soil Sci. 2015, 61, 1439–1453. [Google Scholar] [CrossRef]
  28. Williams, J.R.; Arnold, J.G.; Kiniry, J.R.; Gassman, P.W.; Green, C.H. History of model development at Temple, Texas. Hydrol. Sci. J. 2008, 53, 948–960. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, X.; Tuppad, P.; Williams, J.R. Modeling agricultural management systems with APEX. In Soil Hydrology, Land Use, and Agriculture: Measurement and Modelling, 1st ed.; Shukla, M.K., Ed.; CABI Press: Wallingford, UK, 2011; pp. 117–136. [Google Scholar]
  30. Wang, X.Y.; Gassman, P.W.; Williams, J.R.; Potter, S.; Kemanian, A.R. Modeling the impacts of soil management practices on runoff, sediment yield, maize productivity, and soil organic carbon using APEX. Soil Till. Res. 2008, 101, 78–88. [Google Scholar] [CrossRef]
  31. Cavero, J.; Barros, R.; Sellam, F.; Topcu, S.; Isidoro, D.; Hartani, T.; Lounis, A.; Ibrikci, H.; Cetin, M.; Williams, J. APEX simulation of best irrigation and N management strategies for off-site N pollution control in three Mediterranean irrigated watersheds. Agric. Water Manag. 2012, 103, 88–99. [Google Scholar] [CrossRef] [Green Version]
  32. Meki, M.; Marcos, J.; Atwood, J.; Norfleet, L.; Steglich, E.; Williams, J.; Gerik, T. Effects of site-specific factors on corn stover removal thresholds and subsequent environmental impacts in the Upper Mississippi River Basin. J. Soil Water Conserv. 2011, 66, 386–399. [Google Scholar] [CrossRef] [Green Version]
  33. Meki, M.N.; Atwood, J.D.; Norfleet, L.M.; Williams, J.R.; Gerik, T.J.; Kiniry, J.R. Corn residue removal effects on soybean yield and nitrogen dynamics in the Upper Mississippi River Basin. Agroecol. Sustain. Food 2013, 37, 379–400. [Google Scholar] [CrossRef]
  34. Zhang, B.B.; Feng, G.; Kong, X.B.; Lal, R.; Ouyang, Y.; Adeli, A.; Jenkins, J.N. Simulating yield potential by irrigation and yield gap of rainfed soybean using APEX model in a humid region. Agric. Water Manag. 2016, 177, 440–453. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, B.B.; Feng, G.; Read, J.J.; Kong, X.B.; Ouyang, Y.; Adeli, A.; Jenkins, J.N. Simulating soybean productivity under rainfed conditions for major soil types using APEX model in East Central Mississippi. Agric. Water Manag. 2016, 177, 379–391. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, B.B.; Feng, G.; Ahuja, L.R.; Kong, X.B.; Ouyang, Y.; Adeli, A.; Jenkins, J.N. Soybean crop-water production functions in a humid region across years and soils determined with APEX model. Agric. Water Manag. 2018, 204, 180–191. [Google Scholar] [CrossRef]
  37. Luo, Y.; Wang, H.Y. Modeling the impacts of agricultural management strategies on crop yields and sediment yields using APEX in Guizhou Plateau, southwest China. Agric. Water Manag. 2019, 216, 325–338. [Google Scholar] [CrossRef]
  38. Plotkin, S.; Wang, X.; Potter, T.L.; Bosch, D.D.; Williams, J.R.; Hesketh, E.S.; Bagdon, J.K. APEX calibration and validation of water and herbicide transport under U.S. southern atlantic coastal plain conditions. Trans. ASABE 2013, 56, 43–60. [Google Scholar] [CrossRef]
  39. Van Liew, M.W.; Wortmann, C.S.; Moriasi, D.N.; King, K.W.; Flanagan, D.C.; Veith, T.L.; McCarty, G.W.; Bosch, D.D.; Tomer, M.D. Evaluating the APEX Model for simulating streamflow and water quality on ten agricultural watersheds in the US. Trans. ASABE 2017, 60, 123–146. [Google Scholar]
  40. Jiang, Y.; Wang, X.Q.; Zhao, J.; Jiang, M.; Cao, L.S.; Li, C.; Zeng, Z.H. High-yield variety screening on direct-seeding spring peanut in Heilonggang low plain area. J. Hebei Agric. Univ. 2017, 40, 26–31. [Google Scholar]
  41. Boote, K. Growth stages of peanut (Arachis hypogaea L.). Peanut Sci. 1982, 9, 35–40. [Google Scholar] [CrossRef]
  42. Williams, J.R.; Izaurralde, R.C. The APEX model. In Watershed Models; Singh, V.P., Frevert, D.K., Eds.; CRC Press, Taylor & Francis: Boca Raton, FL, USA, 2006; pp. 437–482. [Google Scholar]
  43. Gassman, P.W.; Williams, J.R.; Wang, X.; Saleh, A.; Osei, E.; Hauck, L.M.; Izaurralde, R.C.; Flowers, J.D. The Agricultural Policy Environmental eXtender (APEX) Model: An Emerging tool for Landscape and Watershed Environmental Analyses; Center for Agricultural and Rural Development, Iowa State University: Ames, IA, USA, 2009. [Google Scholar]
  44. Williams, J.R.; Jones, C.A.; Kiniry, J.R.; Spanel, D.A. The EPIC crop growth model. Trans. ASAE 1989, 32, 497–511. [Google Scholar] [CrossRef]
  45. Powers, S.; Ascough, J., II; Nelson, R.; Larocque, G. Modeling water and soil quality environmental impacts associated with bioenergy crop production and biomass removal in the Midwest USA. Ecol. Model. 2011, 222, 2430–2447. [Google Scholar] [CrossRef]
  46. Williams, J.; Izaurralde, R. Agricultural Policy/Environmental eXtender Model Theoretical Documentation; Blackland Research and Extension Center: Temple, TX, USA, 2008. [Google Scholar]
  47. Penman, H. Natural evaporation from open water, bare soil and grass. Proc. R. Soc. A 1947, 193, 120–145. [Google Scholar]
  48. Gassman, P.W.; Williams, J.R.; Benson, V.W.; Izaurralde, R.C.; Hauck, L.; Jones, C.A.; Atwood, J.D.; Kiniry, J.; Flowers, J.D. Historical Development and Applications of the EPIC and APEX Models; Center for Agricultural and Rural Development, Iowa State University: Ames, IA, USA, 2005. [Google Scholar]
  49. Wang, X.Y.; Williams, J.R.; Gassman, P.W.; Baffaut, C.; Izaurralde, R.C.; Jeong, J.; Kiniry, J.R. EPIC and APEX: Model use, calibration, and validation. Trans. ASABE 2012, 55, 1447–1462. [Google Scholar] [CrossRef]
  50. Zhang, Q.P.; Yang, X.G.; Xue, C.Y.; Yan, W.X.; Yang, J.; Zhang, T.Y.; Bauman, B.A.M.; Wang, H.Q. Analysis of coupling degree between crop w ater requirement of aerobic rice and rainfall in Beijing areas. Transf. CSAE 2007, 23, 51–56. [Google Scholar]
  51. Zhao, Y.Q.; Xiao, D.P.; Qi, Y.Q.; Bai, H.Z. Crop yield and water consumption of different cropping patterns under different precipitation years in North China Plain. Transf. CSAE 2018, 34, 108–116. [Google Scholar]
  52. Hussain, G.; Aljaloud, A.A.; Alshammary, S.F.; Karimulla, S. Effect of saline irrigation on the biomass yield and the protein, nitrogen, phosphorus and potassium composition of alfalfa in a pot experiment. J. Plant Nutr. 1995, 18, 2389–2408. [Google Scholar] [CrossRef]
  53. The Analysis of Split-Plot Experiments. Available online: http://www2.compute.dtu.dk/courses/02429/enotepdfs/eNote-7.pdf (accessed on 21 July 2019).
  54. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  55. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  56. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
  57. Bell, M.; Shorter, R.; Mayer, R. Cultivar and environmental effects on growth and development of peanuts (Arachis hypogaea L.). I. Emergence and flowering. Field Crop. Res. 1991, 27, 17–33. [Google Scholar] [CrossRef]
  58. Batchelor, W.; Jones, J.; Boote, K.; Pinnschmidt, H. Extending the use of crop models to study pest damage. Trans. ASAE 1993, 36, 551–558. [Google Scholar] [CrossRef]
  59. Naab, J.; Singh, P.; Boote, K.; Jones, J.; Marfo, K. Using the CROPGRO-peanut model to quantify yield gaps of peanut in the Guinean Savanna Zone of Ghana. Agron. J. 2004, 96, 1231–1242. [Google Scholar] [CrossRef]
  60. Adomou, M.; Prasad, P.; Boote, K.; Detongnon, J. Disease assessment methods and their use in simulating growth and yield of peanut crops affected by leafspot disease. Ann. Appl. Biol. 2005, 146, 469–479. [Google Scholar] [CrossRef]
  61. Singh, M.P.; Erickson, J.E.; Boote, K.J.; Jones, J.W.; Tillman, B.L.; van Bruggen, A.H. Using the CSM-CROPGRO-peanut model to simulate late leaf spot effects on peanut cultivars of differing resistance. Agron. J. 2013, 105, 1307–1316. [Google Scholar] [CrossRef]
  62. Lanier, J.E.; Jordan, D.L.; Spears, J.F.; Wells, R.; Johnson, P.D.; Barnes, J.S.; Hurt, C.A.; Brandenburg, R.L.; Bailey, J.E. Peanut response to planting pattern, row spacing, and irrigation. Agron. J. 2004, 96, 1066–1072. [Google Scholar] [CrossRef]
  63. Papastylianou, I. Spacing of peanut plants (Arachis hypogaea L.) under irrigation. Eur. J. Agron. 1995, 4, 101–107. [Google Scholar] [CrossRef]
  64. Vadez, V.; Halilou, O.; Hissene, H.M.; Sibiry-Traore, P.; Sinclair, T.R.; Soltani, A. Mapping water stress incidence and intensity, optimal plant populations, and cultivar duration for African groundnut productivity enhancement. Front Plant Sic. 2017, 8, 432. [Google Scholar] [CrossRef]
  65. Ketring, D. Temperature effects on vegetative and reproductive development of peanut. Crop Sci. 1984, 24, 877–882. [Google Scholar] [CrossRef]
  66. Prasad, P.V.; Craufurd, P.; Summerfield, R. Effect of high air and soil temperature on dry matter production, pod yield and yield components of groundnut. Plant Soil 2000, 222, 231–239. [Google Scholar] [CrossRef]
  67. Craufurd, P.; Prasad, P.; Summerfield, R. Dry matter production and rate of change of harvest index at high temperature in peanut. Crop Sci. 2002, 42, 146–151. [Google Scholar] [CrossRef]
  68. Haro, R.J.; Otegui, M.E.; Collino, D.J.; Dardanelli, J.L. Environmental effects on seed yield determination of irrigated peanut crops: Links with radiation use efficiency and crop growth rate. Field Crop. Res. 2007, 103, 217–228. [Google Scholar] [CrossRef]
  69. Collino, D.; Dardanelli, J.; Sereno, R.; Racca, R. Physiological responses of argentine peanut varieties to water stress. Water uptake and water use efficiency. Field Crop. Res. 2000, 68, 133–142. [Google Scholar] [CrossRef]
  70. Stirling, C.; Ong, C.; Black, C. The response of groundnut (Arachis hypogaea L.) to timing of irrigation. J. Exp. Bot. 1989, 40, 1145–1153. [Google Scholar] [CrossRef]
  71. Reddy, T.; Reddy, V.; Anbumozhi, V. Physiological responses of groundnut (Arachis hypogaea L.) to drought stress and its amelioration: A critical review. Plant Growth Regul. 2003, 41, 75–88. [Google Scholar] [CrossRef]
  72. Tubbs, R.S.; Beasley, J.P., Jr.; Culbreath, A.K.; Kemerait, R.C.; Smith, N.B.; Smith, A.R. Row pattern and seeding rate effects on agronomic, disease, and economic factors in large-seeded runner peanut. Peanut Sci. 2011, 38, 93–100. [Google Scholar] [CrossRef]
Agronomy 09 00443 g001 550
Figure 1. Location of Wuqiao County and the North China Plain.
Figure 1. Location of Wuqiao County and the North China Plain.
Agronomy 09 00443 g001
Agronomy 09 00443 g002 550
Figure 2. Monthly precipitation and air temperature (MP = mean precipitation, MT = mean air temperature) during 2017, 2018, and the last 38 years (1981–2018) at Wuqiao County, Hebei Province, China.
Figure 2. Monthly precipitation and air temperature (MP = mean precipitation, MT = mean air temperature) during 2017, 2018, and the last 38 years (1981–2018) at Wuqiao County, Hebei Province, China.
Agronomy 09 00443 g002
Agronomy 09 00443 g003 550
Figure 3. Schematic diagram of the arrangement of planting rows and raised beds under the three seeding rates.
Figure 3. Schematic diagram of the arrangement of planting rows and raised beds under the three seeding rates.
Agronomy 09 00443 g003
Agronomy 09 00443 g004 550
Figure 4. Simulated and measured pod yield of peanut under three sowing dates * three planting densities during calibration ((a), 2017) and validation ((b), 2018) stage. Values are means of three replicates. Broken line, regression line; solid line, 1:1 line. ○ represents T1, □ represents T2, △ represents T3; red represents D1, pink represents D2, blue represents D3.
Figure 4. Simulated and measured pod yield of peanut under three sowing dates * three planting densities during calibration ((a), 2017) and validation ((b), 2018) stage. Values are means of three replicates. Broken line, regression line; solid line, 1:1 line. ○ represents T1, □ represents T2, △ represents T3; red represents D1, pink represents D2, blue represents D3.
Agronomy 09 00443 g004
Agronomy 09 00443 g005 550
Figure 5. Simulated and measured leaf area index (LAI) of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2017. Values are means of three replicates. DAP is the number of days after planting.
Figure 5. Simulated and measured leaf area index (LAI) of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2017. Values are means of three replicates. DAP is the number of days after planting.
Agronomy 09 00443 g005
Agronomy 09 00443 g006 550
Figure 6. Simulated and measured aboveground biomass of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2017. Values are means of three replicates. DAP is the number of days after planting.
Figure 6. Simulated and measured aboveground biomass of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2017. Values are means of three replicates. DAP is the number of days after planting.
Agronomy 09 00443 g006
Agronomy 09 00443 g007 550
Figure 7. Simulated and measured LAI of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2018. Values are means of three replicates. DAP is the number of days after planting.
Figure 7. Simulated and measured LAI of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2018. Values are means of three replicates. DAP is the number of days after planting.
Agronomy 09 00443 g007
Agronomy 09 00443 g008 550
Figure 8. Simulated and measured aboveground biomass of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2018. Values are means of three replicates. DAP is the number of days after planting.
Figure 8. Simulated and measured aboveground biomass of peanut under (a) T1D1, sowing on 25 April with 24 plants m−2, (b) T1D2, sowing on 25 April with 30 plants m−2, (c) T1D3, sowing on 25 April with 36 plants m−2, (d) T2D1, sowing on 5 May with 24 plants m−2, (e) T1D1, sowing on 5 May with 30 plants m−2, (f) T1D1, sowing on 5 May with 36 plants m−2, (g) T1D1, sowing on 15 May with 24 plants m−2, (h) T1D1, sowing on 15 May with 30 plants m−2, (i) T1D1, sowing on 15 May with 36 plants m−2 during 2018. Values are means of three replicates. DAP is the number of days after planting.
Agronomy 09 00443 g008
Agronomy 09 00443 g009 550
Figure 9. Simulated pod yield and maximum LAI of peanut of different seeding rates under rainfed (a) and full irrigation (b). Values are means of 494 data points.
Figure 9. Simulated pod yield and maximum LAI of peanut of different seeding rates under rainfed (a) and full irrigation (b). Values are means of 494 data points.
Agronomy 09 00443 g009
Agronomy 09 00443 g010 550
Figure 10. Simulated pod yield and maximum LAI of peanut of different sowing dates under rainfed (a) and full irrigation (b). Values are means of 304 data points.
Figure 10. Simulated pod yield and maximum LAI of peanut of different sowing dates under rainfed (a) and full irrigation (b). Values are means of 304 data points.
Agronomy 09 00443 g010
Agronomy 09 00443 g011 550
Figure 11. Simulated ranges of peanut (a) pod yield, (b) evapotranspiration (ET), (c) water use efficiency (WUE), and (d) water stress (WS) days under different irrigation scenarios, i.e., W0 (rainfed), W11 (60 mm depth before planting), W12 (60 mm depth at flowering), W21 (60 mm depth before planting and 60 mm depth at flowering), and W22 (60 mm depth before planting and 60 mm depth at pod set). The numbers displayed in parentheses are categorized for precipitation years. The box plots show the 25%, 50%, and 75% percentiles and the caps indicate the 5th and 95th percentiles, with the average shown by a red point. The black points are outliers.
Figure 11. Simulated ranges of peanut (a) pod yield, (b) evapotranspiration (ET), (c) water use efficiency (WUE), and (d) water stress (WS) days under different irrigation scenarios, i.e., W0 (rainfed), W11 (60 mm depth before planting), W12 (60 mm depth at flowering), W21 (60 mm depth before planting and 60 mm depth at flowering), and W22 (60 mm depth before planting and 60 mm depth at pod set). The numbers displayed in parentheses are categorized for precipitation years. The box plots show the 25%, 50%, and 75% percentiles and the caps indicate the 5th and 95th percentiles, with the average shown by a red point. The black points are outliers.
Agronomy 09 00443 g011
Table
Table 1. Soil texture and hydraulic properties of the soil profile at the Wuqiao experimental site.
Table 1. Soil texture and hydraulic properties of the soil profile at the Wuqiao experimental site.
Soil Layer (m)BD
(g cm−3)		Particle Fraction (%)TextureSATC
(mm h−1)		UW
(m m−1)		FC
(m m−1)		pH
SandSiltClay
0–0.21.3316.269.214.6Silt Loam4.40.100.308.1
0.2–0.41.5218.067.714.4Silt Loam1.10.120.278.3
0.4–0.61.4612.182.15.8Silt Loam5.40.120.278.3
0.6–0.81.458.975.915.2Silt Loam7.90.120.298.2
0.8–1.01.464.781.014.3Silt Loam7.10.130.318.3
1.0–1.21.412.686.311.1Silt8.30.180.328.6
Note: BD is bulk density, SATC is soil saturated hydraulic conductivity, UW is the wilting point, FC is field capacity.
Table
Table 2. Agronomic and potential heat unit (PHU) data of peanut in 2017 and 2018.
Table 2. Agronomic and potential heat unit (PHU) data of peanut in 2017 and 2018.
SeasonSowing DateVEoVEsSampling DateHarvesting DateGrowth DurationP
(mm)		PHU (°C)
20174/255/75/86/36/237/269/19/21313622000
5/55/145/156/106/308/69/99/101293552050
5/155/235/236/217/78/109/179/181263552030
20184/255/75/86/37/178/149/19/21316342040
5/55/155/156/67/198/209/039/51246212040
5/155/225/236/137/238/209/89/91186192040
Note: VEo is observed emergence date, VEs is simulated emergence date, P is growing season precipitation.
Table
Table 3. Peanut plant parameters used to calibrate the Agricultural Policy Environmental eXtender (APEX) model. The adjusted parameters are shown in boldface.
Table 3. Peanut plant parameters used to calibrate the Agricultural Policy Environmental eXtender (APEX) model. The adjusted parameters are shown in boldface.
ParametersDescriptionDefaultAdjusted
WABiomass-Energy Ratio30.0030.00
HIHarvest index0.400.40
DMLAMaximum potential leaf area index5.005.00
DLAIFraction of growing season when leaf area declines0.850.75
DLAP1First point on optimal leaf area development curve15.0112.05
DLAP2Second point on optimal leaf area development curve50.9550.70
RLADLeaf area index decline rate parameter1.000.20
RBMDBiomass–energy ratio decline rate parameter0.501.00
PPLP1Plant Population for Crops and Grass-1st Point on curve3.1010.55
PPLP2Plant Population for Crops and Grass-2nd Point on curve10.9040.95
Table
Table 4. Summary statistics of calibration and validation results for pod yield, leaf area index (LAI), and above-ground plant material (ABIOM).
Table 4. Summary statistics of calibration and validation results for pod yield, leaf area index (LAI), and above-ground plant material (ABIOM).
StageIndexnR2PBIASdRMSENRMSE (%)
CalibrationYield90.710.040.790.203.97
LAI360.95−1.670.980.2915.33
ABIOM360.838.530.950.9325.66
ValidationYield90.72−3.470.820.36.37
LAI360.867.950.950.5121.69
ABIOM360.90−1.210.970.8218.11
Note: n is the number of data-points. Unit of RMSE for Yield and ABIOM is t ha−1.

Share and Cite

MDPI and ACS Style

Zhao, J.; Chu, Q.; Shang, M.; Meki, M.N.; Norelli, N.; Jiang, Y.; Yang, Y.; Zang, H.; Zeng, Z.; Jeong, J. Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain. Agronomy 2019, 9, 443. https://doi.org/10.3390/agronomy9080443

AMA Style

Zhao J, Chu Q, Shang M, Meki MN, Norelli N, Jiang Y, Yang Y, Zang H, Zeng Z, Jeong J. Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain. Agronomy. 2019; 9(8):443. https://doi.org/10.3390/agronomy9080443

Chicago/Turabian Style

Zhao, Jie, Qingquan Chu, Mengjie Shang, Manyowa N. Meki, Nicole Norelli, Yao Jiang, Yadong Yang, Huadong Zang, Zhaohai Zeng, and Jaehak Jeong. 2019. "Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain" Agronomy 9, no. 8: 443. https://doi.org/10.3390/agronomy9080443

APA Style

Zhao, J., Chu, Q., Shang, M., Meki, M. N., Norelli, N., Jiang, Y., Yang, Y., Zang, H., Zeng, Z., & Jeong, J. (2019). Agricultural Policy Environmental eXtender (APEX) Simulation of Spring Peanut Management in the North China Plain. Agronomy, 9(8), 443. https://doi.org/10.3390/agronomy9080443

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