Next Article in Journal
Low Salt Treatment Results in Plant Growth Enhancement in Tomato Seedlings
Next Article in Special Issue
Modeling for Predicting the Potential Geographical Distribution of Three Ephedra Herbs in China
Previous Article in Journal
Plant Growth Regulators in Tree Rooting
Previous Article in Special Issue
Towards a Stochastic Model to Simulate Grapevine Architecture: A Case Study on Digitized Riesling Vines Considering Effects of Elevated CO2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crop Model Parameterisation of Three Important Pearl Millet Varieties for Improved Water Use and Yield Estimation

by
Petrus A. Ausiku
1,2,*,
John G. Annandale
1,
Joachim Martin Steyn
1 and
Andrew J. Sanewe
1
1
Department of Plant and Soil Sciences, University of Pretoria, Private Bag X20, Pretoria 0028, South Africa
2
Department of Crop Production and Agriculture Technologies, University of Namibia, Private Bag 13301, Windhoek 9000, Namibia
*
Author to whom correspondence should be addressed.
Plants 2022, 11(6), 806; https://doi.org/10.3390/plants11060806
Submission received: 18 January 2022 / Revised: 28 February 2022 / Accepted: 7 March 2022 / Published: 18 March 2022

Abstract

:
Pearl millet is an important crop for food security in Asia and Africa’s arid and semi-arid regions. It is widely grown as a staple cereal grain for human consumption and livestock fodder. Mechanistic crop growth and water balance models are useful to forecast crop production and water use. However, very few studies have been devoted to the development of the model parameters needed for such simulations for pearl millet. The objectives of the study were to determine crop-specific model parameters for each of three pearl millet varieties (landrace, hybrid, and improved), as well as to calibrate and validate the Soil Water Balance (SWB) model for predicting pearl millet production and water use based on weather data. The SWB was chosen because it is widely used in southern Africa; however, the developed parameters should benefit other models as well. The presented crop-specific parameter values were derived from field observations and literature. Varieties with different phenology, maturity dates and tillering habits were grown under well-watered and well-fertilised conditions for calibration purposes. The calibrated model was used to predict biomass production, grain yield and crop water use. The hybrid’s water use efficiency was higher than that of the landrace and improved variety.

1. Introduction

Pearl millet (Pennisetum glaucum (L.) R.Br.) is an important crop in subsistence agriculture in large areas of the semi-arid tropics of Africa and India, where it is widely grown for grain, fodder and fuel [1,2,3]. In terms of annual production, pearl millet is the sixth most important cereal crop in the world [4,5]. In southern Africa, pearl millet is mainly cultivated in Botswana, Namibia, Zambia, Zimbabwe and South Africa [6]. People in northern Namibia are almost entirely dependent on pearl millet for food [7]. Although pearl millet is widely cultivated by resource-poor farmers as a traditional staple grain crop, it is not as well-known as a major grain crop outside subsistence agriculture. Relatively little is grown commercially as it is mostly used as a fodder crop by commercial farmers. Pearl millets are known for their high levels of resilience to climate change effects, and good nutritional properties [8]. Pearl millets, therefore, play a critical role in improving food security for resource-poor farmers. Pennisetum glaucum, a C4 grass, has excellent photosynthetic efficiency and biomass production potential. Millet grain is highly nutritious, with 8–19% protein, low starch, high fibre (1.2%) [9], and higher micronutrient concentrations (iron and zinc) than rice, wheat, maize and sorghum [10]. Despite the clear importance of this crop, pearl millet is mainly grown under dryland conditions, characterised by marginal production environments, and with minimal use of commercial inputs, such as adequate irrigation and fertilisers. Pearl millet genotypes differ greatly in growth pattern and plant structure and generally require irrigation for optimum yield [11,12,13,14].
Recently, models have been extensively used to predict future events for different scenarios [15] and provide information on resource dynamics at individual sites, over regions and the globe [16,17]. According to Kasampalis et al. [18], models play a crucial role in the development of sustainable land management strategies for maximising the sustainability goals of land managers and policymakers. Models fall into two general categories: empirical and more mechanistic [19]. The empirical models attempt to portray observed phenomena without hypothesising how they occurred. Mechanistic models, on the other hand, seek to describe cause and effect. Mechanistic models, therefore, require understanding and mathematical descriptions of the underlying processes [20,21].
Crop models have become increasingly useful for different purposes, primarily in education and technology transfer [22,23], decision making, and as agronomic research tools to predict crop growth and development, reduce yield gaps, increase the understanding of water management and support interpretation of experimental results [24,25,26,27]. Crop growth models are important tools for scientific research and are widely applied in agriculture to make predictions about the agronomic, environmental, and economic consequences of the complex interactions between crop management, soil and weather conditions [28,29,30,31,32,33,34,35]. Crop models, along with short-term field experiments, are developed as tools to assess yield losses associated with the effects of climate change, changes in temperature, environmental stresses, poor crop management and rainfall on crop yield. Moreover, lengthy and expensive field experiments, particularly with a high number of treatments, can be pre-evaluated through a well-proven model to improve field tests and to lower their overall costs [36]. There are many models, and they all need either crop-specific or even variety-specific model parameters, which are often not available, but need to be determined from growth analysis and water balance studies. These parameters help to project future crop growth and yields for different regions, irrigation demand and to understand the effect of crop and soil type on food productivity and future soil fertility [37,38].
Crop growth simulation models have been used to investigate the effects of crop management options such as irrigation timing, amount and fertiliser applications, under different environmental conditions on long-term mean yield [39,40,41,42]. For the models to simulate yields across the full range of possible yields, it must be robust, thus it is important to evaluate model performance under high-yield conditions in which yields approach the yield potential ceiling, as well as in environments that produce lower yield levels under stress. The projection of crop yield potential, the attainable yield and the corresponding yield gap, is vital to assist in meeting the challenge of increasing food production demands of a growing world population [43]. This necessitates an in-depth understanding of crop development and growth, which is influenced by a variety of climatic, edaphic, hydrological, physiological and managerial factors. Agricultural yields are well below attainable levels in many parts of the world and closing yield gaps is widely viewed as an important strategy to secure a sufficient and reliable food supply [38,44]. Fundamentally, yield gaps are caused by deficiencies in biophysical crop growth environments that are not addressed by agricultural management practices [45]. Crop models can be used to estimate potential yields for a site based on weather conditions and soil water-holding characteristics of the site, and then a systems approach can be taken to determine causes of, and possible remedies for minimising the yield gap between the potential and the actual yields. A good understanding of crop water use is needed to optimise irrigation water management at the field level, for irrigation system design, and to influence water and energy savings [46,47]. Field-scale crop models have been subject to a broad range of uncertainty analyses, varying from grain crop yield response to biomass performance [48].
This paper specifically looked into the Soil Water Balance (SWB) model because it is widely used in southern Africa. The SWB model is a mechanistic crop growth model that is well established and widely used [46,47,49,50,51,52,53,54]. SWB is designed to simulate crop growth, water use and soil water balance components from weather, soil and crop data. Evapotranspiration (ET) is calculated from weather data using Penman–Monteith grass reference daily evapotranspiration (ETo) [55]. Soil water movement is simulated through a layered cascading model, where each layer can have its particular physical properties [22]. Precipitation interception, runoff and drainage are also accounted for in the model. Dry mass production is calculated daily from crop emergence to maturity as either transpiration or radiation-limited [56], whereas assimilates are partitioned to distinctive plant organs, depending on the phenological stage determined from thermal time and water stress level. Robust crop parameters are needed to successfully apply models and should be obtained from careful measurements under field conditions [57]. Some model parameters have already been published for several important cultivated crops [58,59,60,61], but there is a lack of crop-specific model parameters for pearl millet varieties. This is particularly true for pearl millet varieties that are grown under irrigation. Therefore, the objectives of this study were to measure growth and water use of improved, hybrid and landrace pearl millet and to derive crop-specific growth parameters for each variety to be used with mechanistic growth simulation models. Furthermore, the parameters are used to calibrate and validate the SWB generic crop model for widely grown varieties in Namibia (Kangara and Kantana) and South Africa (Agrigreen) to enable the establishment of accurate irrigation regimes under well-watered conditions. The model can then also be used to run simulations to predict crop biomass and grain yield response to various precipitation conditions in agricultural fields, which will give better insight into irrigation scheduling and production planning.

2. Results

2.1. Soil Water Balance

The components of the soil water balance for December 2017–April 2018 for the three pearl millet varieties are shown in Table 1 and Table 2. On 14 December 2017, there was 26 mm of rainfall before germination.
It was observed that volumetric soil water content in the root zone of the landrace variety increased by 21 mm from 328 to 350 mm. Similar trends were observed in both hybrid and improved varieties during the same period. During the seedling stage, there was an increase in soil water content due to irrigation and rainfall.
Seasonal ET varied between 598 and 670 mm (Table 2) and varied between 7 and 89 mm for weekly ET (Table 2). The total water input (rainfall and irrigation) in all plots during the crop seasons (December–April) of 2017/2018 were 947 mm for the landrace, 1035 mm for the improved variety and 1111 mm for the hybrid. Irrigation amount increased with decreasing precipitation and increased as canopy cover increased [62]. About 70 mm of rainfall was received in December, 78 mm in January, 91 mm in February, 230 mm in March and 73 mm in April. Generally, genotypic differences in ET were mostly related to the length of the crop growing period (Table 2).
Higher drainage was observed during the last week of March. In all varieties this was due to a high rainfall event. In March about 208 mm of rain fell during the late growth period, concentrated mainly in three days, when the soil was wet from previous irrigations, and while water consumption by the crop was low. The amount of drainage below the root zone (1 m soil depth) varied from 0 to 88 mm during growing seasons. The total drainage observed in the landrace was higher than improved and hybrid. Total drainage value represents 18.9% (landrace), 14.7% (improved) and 13.4% (hybrid) of the water applied by irrigation and rainfall.

2.1.1. Radiation Limited Dry Matter Production

The radiation use efficiency (RUE) is the slope of the relationship between the total dry matter produced and cumulative intercepted radiation. High R2 values (>0.97) were achieved for all three pearl millet varieties (Figure 1), highlighting the very robust positive relationship between these two variables. The calculated RUE was 2.6 × 10−3 kg MJ−1 for the landrace, 1.9 × 10−3 kg MJ−1 for the hybrid and 2 × 10−3 kg MJ−1 for the improved variety grown under optimal conditions.
The values for both the hybrid and improved variety compares well with those found in the literature for pearl millet (1.96 × 10−3 kg MJ−1 [63], 2.0 × 10−3 kg MJ−1 [64] and 2.1 ± 2.4 × 10−3 kg MJ−1 [65,66]). Similarly, the higher RUE calculated for the landrace is just above that reported by Ong and Monteith [67] (2.5 × 10−3 kg MJ−1) but well below the 4 × 10−3 kg MJ−1 found by Ram et al. [68]. The calculated radiation extinction coefficients (Ks) for the three pearl millet varieties vary from 0.32 to 0.53. These values are closely related to those for field measurements taken at noon reported in the literature (0.29 [69] and 0.5 [67]). The landrace had lower Ks due to its vertical and clumped leaves, which allowed more solar radiation to penetrate through the canopy, leading to more dry matter production [5].

2.1.2. Water Limited Dry Matter Production

Dry matter water ratio (DWR) values determined for the three varieties are presented in Table 3. The hybrid (Agrigreen) had substantially higher water use efficiency characteristics compared to the improved and landrace varieties for grain. The value of WUE (9.16 kgm−3) for landrace biomass was in the range of that reported by Payne [69].
Table 3. Measured crop-specific growth parameters for three pearl millet varieties: Hybrid (Agrigreen), landrace (Kantana) and improved (Kangara) and other values were sourced from literature *.
Table 3. Measured crop-specific growth parameters for three pearl millet varieties: Hybrid (Agrigreen), landrace (Kantana) and improved (Kangara) and other values were sourced from literature *.
Crop ParameterValueLiterature ValueReference
Hybrid (Agrigreen)Landrace (Kantana)Improved (Kangara)
Canopy extinction coefficient for photosynthetic active radiation KPAR0.430.400.420.64–0.42[70,71]
Canopy extinction coefficient for total solar radiation Ks0.310.280.300.49[72]
Radiation use efficiency, RUE (kg MJ−1)0.00190.00260.0020.0003–0.00261[64,65,73,74,75,76]
Dry matter/transpiration ratio corrected for vapour pressure deficit, DWR (Pa)11.315.811.2
Water use efficiencygrain, WUE (kg m−3)1.431.291.233.16–10.4[11,12]
Water use efficiencybiomass, WUE (kg m−3)4.839.164.698.56–39.6[11,12]
Specific leaf area, SLA (m2 kg−1)22.4919.9123.3411.98–33[77,78,79]
Leaf-stem partition parameter, PART (m2 kg−1)2.761.272.48
Maximum root depth (m)1.001.001.001.8[80,81,82]
Maximum crop height (m)2.874.222.804.9[6]
Maximum transpiration rate (mm d−1)9 *9 *9 *9.2[83]
Base temperature (°C)10 *10 *10 *10–12[67,84,85]
Optimum temperature (°C)33 *33 *33 *33–34[67,84]
Maximum temperature (°C)45 *45 *45 *45–47[67,84,85]
Emergence day degrees (°C d)60646060[67]
Flowering day degrees (°C d)9001058832954–1265[84]
Maturity day degrees (°C d)1686212414801552–1714[84]
Transition day degrees (°C d)670780655415–621[84]
Total dry matter yield at emergence (kg m−2)0.00190.00190.0019
Stress index0.30 *0.30 *0.30 *0.30–0.50[86]
Values with superscript star (*) are adopted from literature.
The grain yield WUE of the hybrid (1.43 kgm−3) was higher than that recorded for both the landrace (1.29 kgm−3) and improved (1.23 kgm−3) variety, but lower than the seasonal average determined by Payne [69].

2.1.3. Growing Degree Days

Growing degree days (GDD) are widely used for describing the developmental response of crops to temperature. GDDs recorded from emergence until senescence for the three varieties ranged between 1480–2124 day degrees (Table 3). All three of the pearl millet varieties had similar emergence GDD. The improved variety needed the lowest thermal time to flowering and maturity., The improved (Kangara) pearl millet variety reached flowering after 832 °C d, while for the landrace, the flowering stage was reached at 1058 °C d and the hybrid required 900 °C d to flower. Similar findings were also reported by Singh et al. [84]. The variation in GDDs is due to differences in developmental rates among the three varieties.

2.1.4. Radiation Interception

Maximum fractional interception of radiation of the three pearl millet varieties was almost identical and reached a maximum of 0.94–0.98, similar to what was reported by Squire et al. [65]. The highest KPAR for the hybrid pearl millet variety was 0.43, with an LAI of 6.79 which intercepted 96% of solar radiation. The landrace and improved variety recorded a lower KPAR value than for the hybrid. The extinction coefficient for solar radiation (Ks) was 0.30 for the improved pearl millet variety.
The radiation interception and canopy LAI varied from planting to harvesting and showed a similar trend to those observed by other authors [84]. The LAI (9.59) and FIPAR (98%) were generally higher for the landrace at the heading stage. Similarly, the improved variety recorded its highest LAI (6.99) and FIPAR (94%) also at heading. However, for the hybrid, the greatest LAI and FIPAR were observed at grain filling. Calculated Ks values for improved, landrace and hybrid varieties for well-watered treatments are reported in Table 3 and Figure 2.

2.1.5. Above-Ground Biomass Partitioning

Partitioning of total dry matter to different plant parts was observed to differ between varieties. For example, the leaf to stem (stem + leaf sheath) partition parameter was 2.76 m2 kg−1 for the hybrid, 2.48 m2 kg−1 for the improved variety and 1.27 m2 kg−1 for the landrace (Figure 3 and Table 3). The specific leaf area (SLA) and the leaf-stem dry matter partitioning parameter (p) have to be known to compute DM partitioning with SWB [46]. The dry matter produced is partitioned among roots, leaves, stems and grain. Partitioning depends on the phenological development stage of the crop, as well as water stress [87]. Figure 3 illustrates the correlation between canopy dry matter (CDM) and (SLA CDM)/LAI-1 for each of the three pearl millet varieties. The slope of the regression line, which is forced through the origin, represents p in m2 kg−1. Values of p for the three pearl millet varieties are summarised in Table 3.

2.1.6. Specific Leaf Area

The values of SLA presented in Table 3 and Figure 4 are averages obtained over the entire growing season. The average values recorded for SLA were 22.49 m2 kg−1 for the hybrid, 19.91 m2 kg−1 for the landrace and 23.34 m2 kg−1 for the improved variety. Specific leaf area (SLA) is vital in leaf and plant functioning and regulates the strategy of resource acquisition and plant growth [88]. SLA is presented in Figure 4 as a mean of all leaves measured over the growing season. The measured SLA (m2 kg–1) ranged between 19.91 and 23.34 m2 kg−1, which are in accordance with values of 18 and 33 m2 kg−1 reported by Singh [77], Van Heemst [78] and Monteith [79]. It is clear that SLA is higher at the beginning of the season and the use of a single value for this parameter in SWB is perhaps not ideal.

2.2. Model Calibration

The parameters derived from the calibration data set were used to calibrate the SWB model for each of the three pearl millet varieties (Table 3). Field measurements of LAI, FIPAR, biomass produced, grain yield and measured soil water deficits were used to calibrate the model. Root depth was not measured but estimated from depth of soil water extraction. The crop-specific growth parameters of the three varieties are shown in the list of values in Table 3 under well-watered conditions. Generally, most parameters correspond well to those reported by others [11,12,72] for pearl millet. The three varieties were grown in mean temperature ranges between 13 °C and 32 °C. Pearl millet is mostly cultivated at an optimum temperature of 25 °C [86] to grow and partition new dry matter into different plant parts. Values (1 m) of maximum root depth (RDmax) were generally in the range (0.9 m) of those reported by Gregory and Reddy [82]. The values of the base, optimal and maximum temperatures for pearl millet development were taken from Ong and Monteith [67], Garcia-Huidobro, Monteith, and Squire [85], and Singh, Joshi and Singh [84]. Maximum plant height (Hmax) for the landrace was 4.22 m, while for the hybrid and improved variety, measured Hmax values were 2.87 m and 2.80 m, respectively. On average, the hybrid and improved pearl millet variety were 32% shorter compared to the landrace.
The SWB model was calibrated for the Kantana, Kangara or Agrigreen varieties using the parameters obtained from the calibration data set (Table 3). Calibration was based on field measurements of measured soil water deficits. Soil water deficits were predicted with reasonable accuracy for the calibration data set, with R2 = 0.70, D = 0.85 and MAE = 47% for Kantana, R2 = 0.62, D = 0.83 and MAE = 44% for Kangara and R2 = 0.70, D = 0.78 and MAE = 93% for Agrigreen (Figure 5a–c).

2.3. Model Validation

SWB validation for leaf area index (LAI) of the varieties grown under different irrigation and weather conditions was obtained from three years of independent experimental data and are presented in Figure 6. LAI was predicted with reasonable to good accuracy for all three varieties grown under a wide range of conditions, giving confidence in the reliability of model output.
The validation of the above-ground dry matter (AGDM) and harvestable dry matter (HDM) yields were predicted with reasonable to good accuracy for pearl millet varieties that were irrigated only once in two weeks (Figure 7c,f,i). Throughout the 2016 to 2017 growing season, the simulated HDM (grain) and AGDM yields followed the same trend as the actual values but were frequently overestimated by the model (Figure 7). The overestimate in 2016 may be due to late planting of the experiment, while the 2017 experiment suffered from frost, resulting in lower actual yield then what is projected by model.
The soil water deficits were less accurately estimated, particularly from mid-season onwards, when the crop was rainfed or irrigated less regularly, using the validation data from the 2016–2018 growing season (Figure 8; R2 = 0.82, D = 0.72 and MAE = 75% for hybrid; R2 = 0.87, D = 0.81 and MAE = 83% for improved; and R2 = 0.77, D = 0.84 and MAE = 64% for landrace). These findings imply that the SWB model was able to accurately simulate crop growth, water use and yield for any of the varieties studied under various water supply conditions. As a result, the SWB model can be used with confidence to predict growth, water use and yield of these varieties under various scenarios for future planning.

3. Discussion

3.1. Soil Water Balance

A total of 542 mm of rain was captured between the date of sowing and harvest, 8% more than the optimum rainfall requirement of pearl millet, indicating that the season was good. The monthly distribution of rainfall in terms of rainy days per month was variable. The seasonal rainfall intensity and monthly distribution, along with cumulative seasonal rainfall variation, cause production uncertainties [89]. During the growing season, ET was different for each of the pearl millet varieties. Seasonal ET for each of the varieties was always less than cumulative precipitation, indicating that soil water supply was not being exhausted. However, at harvest, all three varieties’ seasonal ET was greater than total rainfall, suggesting some extraction of water stored at lower soil depths when no supplemental irrigation was supplied. Under semi-arid conditions, irrigation is one of the most important contributions to ET [90] and this contribution differed depending on the amount of water stored in the soil and crop variety. The water balance data, therefore, suggest an increased risk of exhausting available soil water content due to greater ET, especially for the local landrace and hybrid variety. The greater ET was associated with substantially increased grain yield and above-ground dry matter production. The higher seasonal crop ET recorded for the landrace (670 mm) than the hybrid and improved variety was in part due to late maturity and genotypic differences. In the first two months, no heavy rain occurred early in the growth period, so the crop ET was similar in all varieties.
ET of the landrace was higher due to the increased number of days to maturity. During the seedling stage, there was an increase in soil water content due to irrigation and rainfall. Irrigation amount increased with decreasing precipitation and increased with crop growth [62].
Seasonal irrigation amount was high for the hybrid variety and lower for the landrace. In the first two months, no heavy rain occurred early in the growth period, so the crop ET was similar in all varieties. Higher drainage was observed during the last week of March in all varieties due to a high rainfall event. In March about 208 mm of rain fell during the late growth period, concentrated mainly in three days, when the soil was wet from previous irrigations, and while water consumption by the crop was low. The higher seasonal crop ET recorded for the landrace (670 mm) than the hybrid and improved variety was in part due to late maturity.

3.2. Radiation Limited Dry Matter Production

The top dry matter accumulated was high in the landrace variety as it absorbed more radiation over time. This finding is supported by Monteith [66], Sinclair and Muchow [91], who argue that plant growth rate depends on the quantity of radiation intercepted by the canopy which is converted into dry matter accumulation. Further, it was observed among the varieties that the landrace required the highest growing degree days to attain physiological maturity. This is due to the longer growing duration required by the landrace to reach maturity. The RUEs were different for each of the three varieties. According to Monteith [92], Carretero et al. [93], Koester et al. [94] and Teixeira et al. [95], the RUE value depends on the cultivar, water supply, nutrient status, and disease presence. In addition, RUE is rooted in the fundamental relationship between radiant energy use and the accumulation of plant biomass and is influenced by plant phenology and environmental factors [92,96]. Jovanovic et al. [46] also reported that the projected values of RUE denote a lower limit for utilisation of radiation to produce dry matter because root dry matter was excluded in the computations of dry matter for the pearl millet varieties.
Dry matter production is very sensitive to the DWR and moderately sensitive to RDMAX and available water, both of which determine the amount of water used by the plant. Therefore, dry matter production can also be estimated by considering water use efficiency, another important model parameter. In SWB, a simple gas exchange model that acknowledges the tight link between transpiration corrected for dry matter production and vapour pressure deficit (VPD) is used [44]. It was not possible to accurately measure root dry matter, hence this was excluded from the calculation of DWR. Since it was difficult to determine transpiration and total dry matter, ET and above-ground dry matter were used to predict a lower limit value of DWR, after which DWR was adjusted upwards to account for root dry matter yield during model calibration.

3.3. Light Interception and Biomass Accumulation

The landrace had lower Ks due to its vertical and clumped leaves, which allowed more solar radiation to penetrate through the canopy, leading to more dry matter production [5]. The improved variety had higher Ks due to its erect leaves, which suggests that PAR was not widely distributed within the canopy, resulting in lower dry matter production. The results are consistent with numerous other studies, which indicate that Ks are correlated with increases in dry matter production as well as final grain yield [6,7,8]. The high RUE observed for the landrace variety could have been associated with the low Ks, as negative correlations between RUE and Ks have been reported for wheat [97] and pigeon pea [98]. The extinction coefficient for Kantana was low (Table 3), suggesting less interception than in both Agrigreen and Kangara. This relationship further suggested that the LAI required to intercept the same amount of PAR was greater in Kantana (landrace) than in Agrigreen (hybrid) and Kangara (improved). The landrace intercepted more PAR than both the hybrid and improved varieties because the landrace had a 20% greater LAI than the improved variety. Pearl millet, a C4 plant, is an efficient crop with a large leaf area index [99]. However, the higher efficiencies were observed due to lower atmospheric saturation vapour pressure deficits of 0.78, 0.83 and 0.82 kPa in Table 3, than other experimental conditions (3.7–5.2 kPa) reported by McIntyre [100]. According to Squire [74], RUE ranged from 1.4 to 2.5 g MJ−1 when daily maximum saturation vapour pressure deficits did not exceed 3 kPa. However, lower RUEs (0.3 and 0.8 g MJ−1) were reported where daily maximum saturation vapour pressure deficits were between 3 and 5 kPa [74], less than the finding in this study.
Crop genotype influenced fraction of intercepted photosynthetic active radiation (FIPAR) and LAI in this study. Photosynthetic active radiation (PAR) is the basic source of energy for dry matter production [101], and LAI followed the pattern of intercepted amount of incoming PAR in the canopy [102]. The landrace had a large amount of foliage which resulted in a higher interception of PAR. This finding is supported by Jonckheere et al. [103] and Huang et al. [4,104], who argued that PAR depends on the amount of foliage represented by LAI, and the orientation and distribution of foliage. More PAR reached the lower canopy in the improved (short-season) pearl millet than in both the landrace and hybrid. This finding is in line with Timlin et al. [105], who reported that short-season varieties produce low biomass and intercept less radiation.

3.4. Specific Leaf Area

Specific leaf area (SLA) for all varieties differed during the growing cycles and decreased with time due to the senescence of old leaves. However, a calculated seasonal average value of SLA was used in the SWB model, as supported by Rinaldi [106], that can be used in crop models, without a great source of error in the simulation. SLA has received great attention for accurately simulating the production of assimilates and their distribution to the various plant organs [107]. Both the landrace and hybrid varieties had low SLA and reached leaf senescence later than the improved variety. The SLA values also vary with genotype, modified by weather conditions and factors affecting the plant during the growing period [108,109,110]. Reich et al. [111] reported that SLA is negatively correlated with leaf life span and leaves with low SLA and long life span have lower assimilation rates. The lower SLA of both the hybrid and landrace varieties contributed to long leaf survival, which promotes nutrient retention, enhances long-term photosynthetic rate, nitrogen use efficiency, and protection from desiccation [112,113,114]. The long leaf life span leads to a greater assimilation period, providing the basis for the effective absorption of more solar radiation [115]. The landrace and hybrid also had high leaf weight per leaf area, resulting in low SLA, which is caused by thicker leaves [116]. The high SLA value of the improved variety reduces the amount of assimilation required to produce a given leaf area, which results in an earlier ground cover and consequently a greater light harvest that produces a higher assimilation rate early in the season.

4. Materials and Methods

4.1. Experiment Description

The field experiments of a detailed growth analysis and calibration of three pearl millet varieties were conducted in an open field at the Hatfield Experimental Farm (25°45′ S, 28°16′ E, 1327 m above sea level) of the University of Pretoria, South Africa. The trial was laid out as a randomised complete block design and treatments were randomly assigned to blocks in all years of the experiment. The field was divided into three blocks with nine plots of 45 m2 (5 × 9 m2) each in size, with 1 m paths between plots. Plots were demarcated with raised soil bunds between them to avoid the surface movement of water between adjacent plots. The soil was a Hutton sandy clay loam [11] (Soil Classification Working Group, 1991) (loamy, kaolinitic, mesic, Typic Eutrustox) with depths generally in excess of 1.2 m [72]. Nine treatment combinations were set up with three replications. The experiment consisted of three varieties and three irrigation regimes (Table 4). The variety treatments were a hybrid (V1: Agrigreen), open pollinated variety (OPV) landrace (V2: Kantana) and OPV improved (V3: Kangara). The irrigation regimes were well-watered (I1: irrigated every week to field capacity), an intermediate irrigation level (I2: irrigated every second week to field capacity) and a zero-irrigation control (I0: rainfed).
For calibration of the SWB model, the OPV landrace (Kantana), OPV improved (Kangara) and hybrid (Agrigreen) were planted in a well-managed plot. The well-watered plots were used for growth analysis. The pearl millet crops in the experiment were planted at a density of 18 plants per m2. For plant growth analysis data were measured from eight plants per plot per sampling event at fourteen days after planting (DAP), thereafter every week until the plants reached physiological maturity. Plant heights and destructive harvests were executed fortnightly for each replicate. The plant samples were separated into leaves, stems and panicles, whereafter the leaf area (LA) was measured using a belt-driven leaf area meter model LI-3100 (LI-COR, Lincoln, NE, USA) calibrated to 0.01 cm2. Leaf area at the plant level was calculated as the sum of the areas of each green leaf on a plant, and leaf area index (LAI) was calculated. An electronic balance was used for weighing samples that were oven-dried at 70 °C to constant mass to determine leaf dry matter (LDM) and stem dry matter (SDM). At maturity, grain yield was harvested from 1 m2 per plot.
For monitoring the soil water status as part of the soil water balance, a neutron water meter model 503 DR CPN Hydroprobe (Campbell Pacific Nuclear, Martinez, CA, USA) that was calibrated for the experimental site was used to measure soil water contents at 0.2 m increments to a depth of 1.0 m. Root-zone soil water deficit calculations for irrigation were made over an assumed rooting depth of 1.0 m. Weekly soil water content measurements were taken, and thereafter the field was watered to field capacity using a high-density drip system with drip lines spaced 0.45 m apart with an in-line dripper spacing of 0.30 m, and a delivery rate of 8.9 mm h−1. The soil water deficit to field capacity was calculated for each layer and the average for the four different layers determined. A profile pit was dug at the experimental site and soil samples were taken at 0.20 m intervals to a depth of 1.0 m to determine bulk density (ρb) soil texture and volumetric soil water content (θ) at field capacity (FC) and permanent wilting point (PWP). The ρb of the soils was determined by gently driving a cylinder of known volume horizontally into the side of the profile pit. The soil was removed and dried at 105 °C for 24 h. Volumetric water content at FC was determined, before planting, by saturating a portion of the field which was then left for 48 h to drain before sampling, while PWP was determined at the end of the season by withholding irrigation on a section of the experimental plot until the plants died. Soil samples of a known volume were taken and θ was calculated from the mass loss before and after the soil samples were dried at 105 °C for 48 h. Bulk density varied from 1152–1506 kg·m−3 for the 0–1 m soil profile [11].

4.2. Soil Water Balance

The water applied to the entire root zone by irrigation or rainfall was used to calculate the change in soil water storage at the end of every week in each plot. The soil water balance method [55,117,118,119] was used for the estimation of evapotranspiration (ET) from the experimental plots. The water balance method is based on the conservation of mass, which states that inputs minus outputs equal the change in soil water content (ΔS) of the root zone of a crop, and this is equal to the difference between the water content at end of measurement period Qo, and the water content at the beginning of measurement in the root zone, Qi [120], in a given time interval expressed as:
Δ S = Q o Q i
ET = P + I R D Δ S
where ΔS shows the change in soil water storage over the 1 m soil profile, P is precipitation, I is irrigation, R is runoff and D is deep percolation, all expressed in millimetres. The term R was assumed to be negligible because the experiment was conducted on flat terrain and no high-intensity precipitation was experienced. In addition, irrigation was applied using a drip irrigation system and the rate of application did not surpass that of the soil infiltrability. A negative value for ΔS shows a decrease in the amount of water stored in the soil.

4.3. Soil Water Content Estimation

Soil water was measured weekly at 0.20 m intervals from the soil surface down to 1.00 m. Aluminium access tubes that protruded 0.10 m above soil level were covered to prevent rain or soil from entering. Maximum root depth was estimated from weekly measurements of soil water extraction with the neutron water meter model 503 DR CPN Hydroprobe (Campbell Pacific Nuclear, Martinez, CA, USA) that was calibrated for the experimental site. Root-zone soil water deficit calculations for irrigation were made over an assumed rooting depth of 1.00 m.

4.4. Model Calibration

Weather data, including daily minimum and maximum air temperature, as well as humidity, wind speed, incoming solar radiation and precipitation, were recorded by a nearby automated weather station. The automatic weather station consisted of an LI-200SA pyranometer (Li-Cor, Lincoln, NE, USA) to measure solar radiation, a cup anemometer (R.M. Young, Traverse, MI, USA) to measure average wind speed, a tipping bucket rain gauge (Texas Electronics, Dallas, TX, USA), an HMP60 relative humidity and temperature sensor (Vaisala, Woburn, MA, USA), and a CR200X datalogger (Campbell Scientific, Inc., Logan, UT, USA). For calibration of the SWB model, weather, soil and crop data collected in the field trial were used to determine crop-specific growth parameters for well-irrigated landrace, improved and hybrid pearl millet varieties for the 2018 growing season.

4.4.1. VPD and SI Parameters

Vapour pressure deficit was calculated according to the following equation by Jovanovic, Annandale and Mhlauli [46]:
VPD = e s T m a x e s T m i n 2 e a
where esTmax is saturated vapour pressure at maximum air temperature (Pa), esTmin is saturated vapour pressure at minimum air temperature (Pa) and ea is actual vapour pressure (Pa). Saturated vapour pressure (es) at maximum (Tmax) and minimum air temperature (Tmin) was calculated by replacing Ta (actual air temperature) with Tmax and Tmin (°C) [119], as presented in Equation (4):
e s = 0.611 e x p 17.27 T a T a + 237.3
ea (actual vapour pressure) was calculated from the measured daily Tmax, Tmin, RHmax and RHmin, using the following equation from Allen, Pereira, Raes and Smith [119]:
e a = e s T m i n R H m a x 100 + e s T m a x R H m i n 100 2
As water becomes limiting, actual evapotranspiration will fall below the potential rate. A measure of the ratio of actual to reference evapotranspiration should therefore be an index of crop water status. Water stress occurs when available soil water does not meet potential transpiration demand, and it is calculated as the ratio of actual to potential transpiration. A water stress index can also be used to reduce the accumulation of growing day degrees, depending on phenology calculated with thermal time and modified by water stress. This dry matter is partitioned into roots, stems, leaves, grains, or fruits (Annandale, [49]). Such stress index (SI) is defined by Jackson et al. [121,122] as ranging from 0 (ample water) to 1 (maximum stress) and was computed from a soil water balance [86,123,124]:
SI = T Fi transp P E T
where SI is stress index, T is actual evapotranspiration (mm), and PET is potential evapotranspiration (mm).

4.4.2. Radiation Limited Dry Matter Production

Many crop growth models are radiation driven and require an estimate of radiation use efficiency (RUE). In SWB, the radiation conversion efficiency is a crop-specific parameter used to calculate dry matter production under conditions of radiation limited growth, from Monteith [125], as follows:
DM = RUE i = 1 n Fi RAD R s
where DM is total dry matter production in kg m−2, which includes root mass. RUE is the slope of the regression between intercepted solar radiation and DM production, forced through the origin (kg MJ−1). RUE is also known as the radiation conversion efficiency (Ec), FIRAD is the fraction of daily value incremented each day (i = 1 to n) of solar radiation intercepted by the green, transpiring canopy and Rs is the total fraction of daily value incremented each day (i = 1 to n) of incoming solar radiation in MJ m−2 day−1.

4.4.3. Water Limited Dry Matter Production

Transpired water vapour leaves the plant by diffusion through stomata, and carbon dioxide is taken up for photosynthesis following the same pathway but in the opposite direction. The CO2 taken up is converted to organic carbon, and this process which is accompanied by transpirational water loss can be described by a simple gas exchange dry matter production model [126]. Transpiration limited growth is calculated using the relationship between dry matter accumulation and transpiration, from Tanner and Sinclair [127].
DWR = DM × VPD ET
where, DWR (dry matter water ratio) is the water productivity [124] of the dry matter production per unit crop transpiration (T), adjusted to account for vapour pressure deficit (VPD) [127], and DM is total dry matter produced (kg m−2). In SWB, a simple gas exchange model that acknowledges the tight link between transpiration corrected for dry matter production and vapour pressure deficit (VPD) is used [44].
VPD is the seasonal average vapour pressure deficit. ET is the seasonal (from sowing to maturity) total evapotranspiration obtained using Equation (1). Both VPD and DWR are in Pa, while ET is in mm. Due to the difficulty in determining root dry matter, only top dry matter was measured, resulting in a smaller DWR than required in Equation (8). Similarly, transpiration (T) is difficult to measure, so ET was used, which also gives a lower DWR than the actual value. Evapotranspiration and top dry matter were therefore used to estimate a lower limit value of DWR. The DWR lower limit value was increased during model parameterisation until simulation data were observed in line with measured data. The estimation of lower limit DWR values was performed using measured top dry matter accumulation data at physiological maturity to avoid errors caused by plants losing leaves, as leaf senescence usually occurs after physiological maturity.

4.4.4. Growing Degree Days

Growing degree days for the transition period between planting and emergence, vegetative and reproductive growth and leaf senescence were estimated by calibration against field measurements of phenology for all pearl millet varieties. The growing day degrees (GDD) (d °C) were calculated using the daily average air temperature, from Monteith [125]:
GDD = i = 1 n T m a x + T m i n 2 T b
where GDD is growing day degrees, Tmax is the maximum air temperature in °C, Tmin is the minimum air temperature in °C, Tb is the base temperature in °C, n is the number of days and i is the daily increment until n days. The base temperature for calculating GDD is the minimum threshold temperature below which plant development ceases. Growing degree days were accumulated from the date of sowing to maturity of each pearl millet cultivar. In this research, a Tb value of 10 °C and a maximum temperature of 45 °C were chosen for pearl millet [67]. If the daily average temperature was less than the base temperature, then the GDD value was assumed equal to zero [84]. GDD was used to estimate the duration of the different crop development stages. Growing day degrees for emergence were computed as the difference between average daily temperature and base temperature, starting from sowing date to emergence.

4.4.5. Radiation Interception

A Decagon Sunfleck ceptometer (Decagon Devices, Pullman, WA, USA) was used to measured fractional interception of photosynthetically active radiation on a weekly basis. On clear-sky days, measurements were performed between 10:00 and 12:00. Five reference values above and five readings below the canopy were averaged in a series of measurements. FIPAR (fractional interception of photosynthetic active radiation) was determined using the following formula:
FI PAR = 1 PAR   below   canopy PAR   above   canopy
The Beer–Bouguer law was used to determine the canopy extinction coefficient for PAR [128]:
FI PAR = 1 e K PAR × LAI
The canopy extinction coefficient for PAR is represented by KPAR. Field measurements of LAI and FIPAR were used to calculate KPAR values. Dry matter production can be estimated using KPAR as a function of intercepted PAR. SWB requires the canopy extinction coefficient for solar radiation (Ks), as fractional interception of solar radiation (FIRAD) is used to predict radiation limited dry matter production [125], and for partitioning evapotranspiration into evaporation from the soil surface and crop transpiration [129]. The procedure recommended by Campbell and Van Evert [128] was used to convert KPAR into Ks:
K s = K bd a s
K bd = K PAR a p
a s = a p a n
where Kbd is the canopy radiation extinction coefficient for black leaves and diffuse radiation, ap is the absorptance of PAR, and as is leaf absorptance of near-infrared radiation (NIR, 0.7–3 m). The value of ap was set to 0.8, while the value of an was set to 0.2 [130]. The geometric mean of the absorptances in the PAR and NIR spectrums is referred to as as.

4.4.6. Above-Ground Biomass Partitioning

In SWB, the daily increment in DM is computed to be either limited by transpiration or limited by radiation [56], with the partitioning of assimilates among different plant organs influenced by plant water status. Specific leaf area (SLA) was obtained as the total leaf area (LA) divided by the dry mass of leaves per sampling area.
Specific   leaf   area   SLA = LA LDM
The SLA, LAI, and canopy dry matter (CDM) were used to calculate the leaf-stem partitioning parameter (p) [46]. The leaf-stem partitioning parameter is represented as the slope of the regression line between CDM and (SLA.CDM)/LAI—1 in m2 kg−1. LAI is the leaf area index (m2 m−2), and SLA is the specific leaf area (m2 kg−1). SDM is stem dry matter.
LDM = CDM / 1 + pCDM
SDM = CDM LDM
Grain yield, LDM, CDM and SDM are all measured in kg m−2.

4.5. Data Collection for Model Calibration

For SWB model calibration the following crop growth parameters for each of the crops were recorded: canopy extinction coefficient (Ks) for total solar radiation, dry matter: transpiration ratio, adjusted for vapour pressure deficit (DWR, kg·kg−1·Pa), radiation use efficiency (RUE, kg·MJ−1), leaf-stem dry matter partitioning (p, m2·kg−1) and specific leaf area (SLA, m2·kg−1).

4.6. Model Validation

A three-year dryland field experiment was established to investigate crop response to dry spells at the Hatfield Experimental Farm of the University of Pretoria in Pretoria, South Africa, with planting dates in February 2016, February 2017 and December 2018. Three pearl millet varieties with contrasting agronomic characteristics, namely Kantana, Kangara and Agrigreen, were grown. Kantana is a landrace variety with a longer duration to reach maturity, high plant height and a greater biomass. Landraces have high adaptation to prevalent abiotic stresses. Kangara is an improved early maturity variety released in Namibia in 1998 [6], characterised by high grain yield under rainfed conditions. Agrigreen is a high-yielding hybrid variety, predominantly grown as a fodder crop but is also valued for its grain. Hybrids are developed due to their yield superiority over open pollinated varieties; however, they are often outperformed by landrace and improved varieties under rainfed conditions [131].
These pearl millet varieties were subjected to three different irrigation regimes (rainfed, irrigated every week and every second week to field capacity). The soil was similar to that used to parameterise the model. The experiment was set up in a three-replication randomised complete block. Pearl millet was planted at a 0.45 m inter-row spacing with a 0.14 m intra-row spacing, in plots measuring 4.95 m by 9.5 m. A calibrated neutron probe was used to measure the amount of water in the soil on a weekly basis. Drip irrigation was used to irrigate the plots. The rainfed treatments and every second week irrigated treatments used to validate the model were as follows:
1.
Irrigated weekly to field capacity until the end of the growing season.
2.
Irrigated to field capacity on a fortnightly basis until the end of the growing season.
3.
Rainfed (dryland) until the end of the growing season.

5. Conclusions

All the varieties had almost similar water use efficiency (for grain yield), radiation use efficiency and degree days to emergence. The specific leaf area, canopy extinction coefficient for total solar radiation and flowering day degrees parameters were different among the varieties. The hybrid variety portrayed a determinate growth habit; therefore, if similar varieties grow in this way, hybrid parameters must be used. However, if the variety is late maturing and indeterminate, modellers are advised to use landrace parameters. If the modellers have an early maturity variety that is indeterminate, then the improved variety parameters should be used. The average values for canopy extinction coefficient for photosynthetic active radiation differ largely among the three varieties; therefore, using average values across the varieties is not advisable. The parameters need to be used with the model to evaluate results. Despite the lack of actual long-term yield data to compare to simulations, the authors are confident that the calibrated SWB model will be useful in estimating pearl millet yields under various production conditions. It is, however, recommended that crop parameters should also be determined for other pearl millet varieties to make the SWB model more useful for the management of irrigation, long-term yield projection and planning purposes.

Author Contributions

P.A.A. and J.G.A. conceptualised the presented idea, methodology, design and implementation of the research. J.G.A. and J.M.S. encouraged P.A.A. to investigate and supervised the findings of his work. J.M.S. provided critical feedback and helped shape the research, analysis and manuscript. P.A.A. wrote the manuscript with support from J.G.A., J.M.S. and A.J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Details regarding data availability will be provided during review.

Acknowledgments

The authors acknowledge support by the University of Namibia, University of Pretoria, and German Academic Services and Via Farm (ACIAR).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Andrews, D.; Kumar, K. Pearl millet for food, feed, and forage. Adv. Agron. 1992, 48, 89–139. [Google Scholar]
  2. Bouajila, A.; Lamine, M.; Rahali, F.; Melki, I.; Prakash, G.; Ghorbel, A. Pearl millet populations characterized by Fusarium prevalence, morphological traits, phenolic content, and antioxidant potential. J. Sci. Food Agric. 2020, 100, 4172–4181. [Google Scholar] [CrossRef] [PubMed]
  3. Varshney, R.K.; Shi, C.; Thudi, M.; Mariac, C.; Wallace, J.; Qi, P.; Zhang, H.; Zhao, Y.; Wang, X.; Rathore, A. Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments. Nat. Biotechnol. 2017, 35, 969–976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Brunken, J.; De Wet, J.; Harlan, J. The morphology and domestication of pearl millet. Econ. Bot. 1977, 31, 163–174. [Google Scholar] [CrossRef]
  5. Chandana, P.; Lata, A.M.; Khan, M.A.; Krishna, A. Influence of Nutrient Management Practices on Growth and Yield of Pearl Millet in Melia dubia Based Agri-Silvi System. Int. J. Curr. Microbiol. App. Sci. 2018, 7, 443–454. [Google Scholar] [CrossRef]
  6. Upadhyaya, H.; Reddy, K.; Gowda, C. Pearl millet germplasm at ICRISAT genebank-status and impact. J. SAT Agric. Res. 2007, 3, 5. [Google Scholar]
  7. Baltensperger, D.D. Progress with proso, pearl and other millets. In Trends in New Crops and New Uses; ASHS Press: Alexandria, VA, USA, 2002; pp. 100–103. [Google Scholar]
  8. Jukanti, A.; Gowda, C.L.; Rai, K.; Manga, V.; Bhatt, R. Crops that feed the world 11. Pearl Millet (Pennisetum glaucum L.): An important source of food security, nutrition and health in the arid and semi-arid tropics. Food Secur. 2016, 8, 307–329. [Google Scholar] [CrossRef]
  9. Nambiar, V.S.; Dhaduk, J.; Sareen, N.; Shahu, T.; Desai, R. Potential functional implications of pearl millet (Pennisetum glaucum) in health and disease. J. Appl. Pharm. Sci. 2011, 1, 62. [Google Scholar]
  10. Tako, E.; Reed, S.M.; Budiman, J.; Hart, J.J.; Glahn, R.P. Higher iron pearl millet (Pennisetum glaucum L.) provides more absorbable iron that is limited by increased polyphenolic content. J. Nutr. 2015, 14, 11. [Google Scholar] [CrossRef] [Green Version]
  11. Ausiku, A.P.; Annandale, J.G.; Steyn, J.M.; Sanewe, A.J. Improving Pearl Millet (Pennisetum glaucum) Productivity through Adaptive Management of Water and Nitrogen. Water 2020, 12, 422. [Google Scholar] [CrossRef] [Green Version]
  12. Maman, N.; Lyon, D.J.; Mason, S.C.; Galusha, T.D.; Higgins, R. Pearl millet and grain sorghum yield response to water supply in Nebraska. Agron. J. 2003, 95, 1618–1624. [Google Scholar] [CrossRef] [Green Version]
  13. Gulia, S.; Wilson, J.; Carter, J.; Singh, B. Progress in grain pearl millet research and market development. In Issues in New Crops and New Uses; ASHS Press: Alexandria, VA, USA, 2007; pp. 196–203. [Google Scholar]
  14. Basnayake, J.; Fukai, S.; Ouk, M. Contribution of potential yield, drought tolerance and escape to adaptation of 15 rice varieties in rainfed lowlands in Cambodia. In Proceedings of the Australian Agronomy Conference, Australian Society of Agronomy, Birsbane, Australia, 10–14 September 2006. [Google Scholar]
  15. Singels, A.; Annandale, J.G.; Jager, J.D.; Schulze, R.E.; Inman-Bamber, N.; Durand, W.; Rensburg, L.V.; Heerden, P.V.; Crosby, C.T.; Green, G.C. Modelling crop growth and crop water relations in South Africa: Past achievements and lessons for the future. S. Afr. J. Plant Soil 2010, 27, 49–65. [Google Scholar] [CrossRef]
  16. Bonan, G. Climate Change and Terrestrial Ecosystem Modeling; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  17. Ciais, P.; Sabine, C.; Bala, G.; Bopp, L.; Brovkin, V.; Canadell, J.; Chhabra, A.; DeFries, R.; Galloway, J.; Heimann, M. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P.M., Eds.; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  18. Kasampalis, D.A.; Alexandridis, T.K.; Deva, C.; Challinor, A.; Moshou, D.; Zalidis, G. Contribution of remote sensing on crop models: A review. J. Imaging 2018, 4, 52. [Google Scholar] [CrossRef] [Green Version]
  19. Rengel, Z. Mechanistic simulation models of nutrient uptake: A review. Plant Soil 1993, 152, 161–173. [Google Scholar] [CrossRef]
  20. Van Noordwijk, M.; De Willigen, P.; Ehlert, P. A simple model of P uptake by crops as a possible basis for P fertilizer recommendations. Neth. J. Agric. Sci. 1990, 38, 317–332. [Google Scholar] [CrossRef]
  21. Hoffland, E.; Bloemhof, H.; Leffelaar, P.; Findenegg, G.; Nelemans, J. Simulation of nutrient uptake by a growing root system considering increasing root density and inter-root competition. In Plant Nutrition—Physiology and Applications; Springer: Berlin/Heidelberg, Germany, 1990; pp. 9–15. [Google Scholar]
  22. Vahrmeijer, J.; Annandale, J.G.; Steyn, J.M.; Bristow, K.L. Model parameters of four important vegetable crops for improved water use and yield estimation. Water SA 2018, 44, 528–538. [Google Scholar] [CrossRef] [Green Version]
  23. Murthy, V.R.K. Crop growth modeling and its applications in agricultural meteorology. In Satellite Remote Sensing and GIS Applications in Agricultural Meteorology; World Meteorological Organisation: Geneva, Switzerland, 2004; Volume 235. [Google Scholar]
  24. Diekkrüger, B.; Söndgerath, D.; Kersebaum, K.; McVoy, C. Validity of agroecosystem models a comparison of results of different models applied to the same data set. Ecol. Model. 1995, 81, 3–29. [Google Scholar] [CrossRef]
  25. Eitzinger, J.; Trnka, M.; Hösch, J.; Žalud, Z.; Dubrovský, M. Comparison of CERES, WOFOST and SWAP models in simulating soil water content during growing season under different soil conditions. Ecol. Model. 2004, 171, 223–246. [Google Scholar] [CrossRef]
  26. Kersebaum, K.C. Modelling nitrogen dynamics in soil–crop systems with HERMES. In Modelling Water and Nutrient Dynamics in Soil–Crop Systems; Springer: Berlin/Heidelberg, Germany, 2007; pp. 147–160. [Google Scholar]
  27. Boote, K.J.; Jones, J.W.; Pickering, N.B. Potential uses and limitations of crop models. Agron. J. 1996, 88, 704–716. [Google Scholar] [CrossRef]
  28. Acutis, M.; Ducco, G.; Grignani, C. Stochastic use of the LEACHN model to forecast nitrate leaching in different maize cropping systems. Eur. J. Agron. 2000, 13, 191–206. [Google Scholar] [CrossRef]
  29. Donatelli, M.; Acutis, M.; Danuso, F.; Mazzetto, F.; Nasuelli, P.; Nelson, R.; Omicini, A.; Speroni, M.; Trevisan, M.; Tugnoli, V. Integrated procedures for evaluating technical, environmental and economical aspects in farms: The SIPEAA project. Proc. Seventh ESA 2002, 15–18. [Google Scholar]
  30. Bechini, L.; Bocchi, S.; Maggiore, T. Spatial interpolation of soil physical properties for irrigation planning. A simulation study in northern Italy. Eur. J. Agron. 2003, 19, 1–14. [Google Scholar] [CrossRef]
  31. Gömann, H.; Kreins, P.; Kunkel, R.; Wendland, F. Model based impact analysis of policy options aiming at reducing diffuse pollution by agriculture—A case study for the river Ems and a sub-catchment of the Rhine. Environ. Model. Softw. 2005, 20, 261–271. [Google Scholar] [CrossRef]
  32. Lewis, D.; McGechan, M.; McTaggart, I.P. Simulating field-scale nitrogen management scenarios involving fertiliser and slurry applications. Agric. Syst. 2003, 76, 159–180. [Google Scholar] [CrossRef]
  33. Johnsson, H.; Larsson, M.; Mårtensson, K.; Hoffmann, M. SOILNDB: A decision support tool for assessing nitrogen leaching losses from arable land. Environ. Model. Softw. 2002, 17, 505–517. [Google Scholar] [CrossRef]
  34. Rossing, W.; Meynard, J.M.; Van Ittersum, M. Model-based explorations to support development of sustainable farming systems: Case studies from France and the Netherlands. Eur. J. Agron. 1997, 7, 271–283. [Google Scholar] [CrossRef]
  35. Wolf, J.; Beusen, A.; Groenendijk, P.; Kroon, T.; Rötter, R.; Van Zeijts, H. The integrated modeling system STONE for calculating nutrient emissions from agriculture in the Netherlands. Environ. Model. Softw. 2003, 18, 597–617. [Google Scholar] [CrossRef]
  36. Whisler, F.; Acock, B.; Baker, D.; Fye, R.; Hodges, H.; Lambert, J.; Lemmon, H.; McKinion, J.; Reddy, V. Crop simulation models in agronomic systems. In Advances in Agronomy; Elsevier: Amsterdam, The Netherlands, 1986; Volume 40, pp. 141–208. [Google Scholar]
  37. Evans, L.T. Crop Evolution, Adaptation and Yield; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  38. De Vries, F.P.; Rabbinge, R.; Groot, J. Potential and attainable food production and food security in different regions. Philos. Trans. R. Soc. B Biol. Sci. 1997, 352, 917–928. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, E.; Robertson, M.; Hammer, G.; Carberry, P.S.; Holzworth, D.; Meinke, H.; Chapman, S.; Hargreaves, J.; Huth, N.; McLean, G. Development of a generic crop model template in the cropping system model APSIM. Eur. J. Agron. 2002, 18, 121–140. [Google Scholar] [CrossRef]
  40. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.J.; Holzworth, D.; Huth, N.I.; Hargreaves, J.N.; Meinke, H.; Hochman, Z. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef] [Green Version]
  41. Aggarwal, P.K.; Kalra, N.; Chander, S.; Pathak, H. InfoCrop: A dynamic simulation model for the assessment of crop yields, losses due to pests, and environmental impact of agro-ecosystems in tropical environments. I. Model description. Agric. Syst. 2006, 89, 1–25. [Google Scholar] [CrossRef]
  42. Ewert, F.; Rötter, R.P.; Bindi, M.; Webber, H.; Trnka, M.; Kersebaum, K.C.; Olesen, J.E.; van Ittersum, M.K.; Janssen, S.; Rivington, M. Crop modelling for integrated assessment of risk to food production from climate change. Environ. Model. Softw. 2015, 72, 287–303. [Google Scholar] [CrossRef]
  43. Asseng, S.; Ewert, F.; Martre, P.; Rötter, R.P.; Lobell, D.B.; Cammarano, D.; Kimball, B.A.; Ottman, M.J.; Wall, G.; White, J.W. Rising temperatures reduce global wheat production. Nat. Clim. Change 2015, 5, 143. [Google Scholar] [CrossRef]
  44. Phalan, B.; Green, R.; Balmford, A. Closing yield gaps: Perils and possibilities for biodiversity conservation. Philos. Trans. R. Soc. B Biol. Sci. 2014, 369, 20120285. [Google Scholar] [CrossRef] [Green Version]
  45. Mueller, N.D.; Gerber, J.S.; Johnston, M.; Ray, D.K.; Ramankutty, N.; Foley, J.A. Closing yield gaps through nutrient and water management. Nature 2012, 490, 254–257. [Google Scholar] [CrossRef]
  46. Jovanovic, N.; Annandale, J.; Mhlauli, N. Field water balance and SWB parameter determination of six winter vegetable species. Water SA 1999, 25, 191–196. [Google Scholar]
  47. Beletse, Y.; Annandale, J.; Steyn, J.; Hall, I.; Aken, M. Can crops be irrigated with sodium bicarbonate rich CBM deep aquifer water? Theoretical and field evaluation. Ecol. Eng. 2008, 33, 26–36. [Google Scholar] [CrossRef]
  48. Asseng, S.; Ewert, F.; Rosenzweig, C.; Jones, J.W.; Hatfield, J.L.; Ruane, A.C.; Boote, K.J.; Thorburn, P.J.; Rötter, R.P.; Cammarano, D. Uncertainty in simulating wheat yields under climate change. Nat. Clim. Change 2013, 3, 827–832. [Google Scholar] [CrossRef] [Green Version]
  49. Annandale, J.; Jovanovic, N.; Benade, N.; Tanner, P. Modelling the long-term effect of irrigation with gypsiferous water on soil and water resources. Agric. Ecosyst. Environ. 1999, 76, 109–119. [Google Scholar] [CrossRef]
  50. Jovanovic, N.; Annandale, J.; Nel, A. Calibration and validation of the SWB model for sunflower (Helianthus annuus L.). S. Afr. J. Plant Soil 2000, 17, 117–123. [Google Scholar] [CrossRef]
  51. Du Sautoy, N.; Jovanovic, N.; Annandale, J. Water balance simulation of a peach orchard using the SWB model. S. Afr. J. Plant Soil 2003, 20, 169–175. [Google Scholar] [CrossRef] [Green Version]
  52. Annandale, J.; Jovanovic, N.; Campbell, G.; Du Sautoy, N.; Lobit, P. Two-dimensional solar radiation interception model for hedgerow fruit trees. Agric. For. Meteorol. 2004, 121, 207–225. [Google Scholar] [CrossRef]
  53. Fessehazion, M.K.; Annandale, J.G.; Everson, C.S.; Stirzaker, R.J.; Tesfamariam, E.H. Evaluating of soil water balance (SWB-Sci) model for water and nitrogen interactions in pasture: Example using annual ryegrass. Agric. Water Manag. 2014, 146, 238–248. [Google Scholar] [CrossRef]
  54. Tesfamariam, E.H.; Annandale, J.G.; Steyn, J.M.; Stirzaker, R.J.; Mbakwe, I. Use of the SWB-Sci model for nitrogen management in sludge-amended land. Agric. Water Manag. 2015, 152, 262–276. [Google Scholar] [CrossRef] [Green Version]
  55. Pereira, L.S.; Allen, R.G.; Smith, M.; Raes, D. Crop evapotranspiration estimation with FAO56: Past and future. Agric. Water Manag. 2015, 147, 4–20. [Google Scholar] [CrossRef]
  56. Goudriaan, J.; Van Laar, H. Modelling Potential Crop Growth Processes: Textbook with Exercises; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 2. [Google Scholar]
  57. Bechini, L.; Bocchi, S.; Maggiore, T.; Confalonieri, R. Parameterization of a crop growth and development simulation model at sub-model components level. An example for winter wheat (Triticum aestivum L.). Environ. Model. Softw. 2006, 21, 1042–1054. [Google Scholar] [CrossRef]
  58. Fessehazion, M.K.; Stirzaker, R.J.; Annandale, J.G.; Everson, C.S. Improving nitrogen and irrigation water use efficiency through adaptive management: A case study using annual ryegrass. Agric. Ecosyst. Environ. 2011, 141, 350–358. [Google Scholar] [CrossRef] [Green Version]
  59. Donatelli, M.; Stöckle, C.; Ceotto, E.; Rinaldi, M. Evaluation of CropSyst for cropping systems at two locations of northern and southern Italy. Eur. J. Agron. 1997, 6, 35–45. [Google Scholar] [CrossRef]
  60. Giardini, L.; Berti, A.; Morari, F. Simulation of two cropping systems with EPIC and CropSyst models. Ital. J. Agron. 1998, 2, 29–38. [Google Scholar]
  61. Bellocchi, G.; Silvestri, N.; Mazzoncini, M.; Meini, S. Using the CropSyst Model in Continuous Rainfed Maize (Zea mais L.) under Alternative Management Options. Ital. J. Agron. 2002, 6, 43–56. [Google Scholar]
  62. Sun, H.; Shen, Y.; Yu, Q.; Flerchinger, G.N.; Zhang, Y.; Liu, C.; Zhang, X. Effect of precipitation change on water balance and WUE of the winter wheat–summer maize rotation in the North China Plain. Agric. Water Manag. 2010, 97, 1139–1145. [Google Scholar] [CrossRef]
  63. Squire, G.; Gregory, P.; Monteith, J.; Russell, M.; Singh, P. Control of water use by pearl millet (Pennisetum typhoides). Exp. Agric. 1984, 20, 135–149. [Google Scholar] [CrossRef] [Green Version]
  64. Reddy, M.; Willey, R. Growth and resource use studies in an intercrop of pearl millet/groundnut. Field Crop. Res. 1981, 4, 13–24. [Google Scholar] [CrossRef] [Green Version]
  65. Squire, G.; Marshall, B.; Terry, A.; Monteith, J. Response to temperature in a stand of pearl millet: VI. Light interception and dry matter production. J. Exp. Bot. 1984, 35, 599–610. [Google Scholar] [CrossRef]
  66. Craufurd, P.; Bidinger, F. Effect of the duration of the vegetative phase on crop growth, development and yield in two contrasting pearl millet hybrids. J. Agric. Sci. 1988, 110, 71–79. [Google Scholar] [CrossRef] [Green Version]
  67. Ong, C.; Monteith, J. Response of pearl millet to light and temperature. Field Crop. Res. 1985, 11, 141–160. [Google Scholar] [CrossRef] [Green Version]
  68. Ram, N.; Sheoran, K.; Sastry, C.S. Radiation efficiency and its efficiency in dry biomass production of pearl millet cultivars. Ann. Agric. Res. 1999, 20, 286–291. [Google Scholar]
  69. Payne, W.A. Managing yield and water use of pearl millet in the Sahel. Agron. J. 1997, 89, 481–490. [Google Scholar] [CrossRef] [Green Version]
  70. Jarwal, S.; Singh, P.; Virmani, S. Influence of planting geometry on photosynthetically active radiation interception and dry matter production relationships in pearl millet. Biomass 1990, 21, 273–284. [Google Scholar] [CrossRef] [Green Version]
  71. Jovanovic, N.; Annandale, J. Measurement of radiant interception of crop canopies with the LAI-2000 plant canopy analyzer. S. Afr. J. Plant Soil 1998, 15, 6–13. [Google Scholar] [CrossRef] [Green Version]
  72. Annandale, J.; Campbell, G.; Olivier, F.; Jovanovic, N. Predicting crop water uptake under full and deficit irrigation: An example using pea (Pisum sativum L. cv. Puget). Irrig. Sci. 2000, 19, 65–72. [Google Scholar] [CrossRef]
  73. Azam-Ali, S.; Gregory, P.; Monteith, J. Effects of planting density on water use and productivity of pearl millet (Pennisetum typhoides) grown on stored water. II. Water use, light interception and dry matter production. Exp. Agric. 1984, 20, 215–224. [Google Scholar] [CrossRef] [Green Version]
  74. Squire, G.R. The Physiology of Tropical Crop Production; CAB International: Wallingford, UK, 1990. [Google Scholar]
  75. Begue, A.; Desprat, J.; Imbernon, J.; Baret, F. Radiation use efficiency of pearl millet in the Sahelian zone. Agric. For. Meteorol. 1991, 56, 93–110. [Google Scholar] [CrossRef]
  76. Craufurd, P.; Bidinger, F. Potential and realized yield in pearl millet (Pennisetum americanum) as influenced by plant population density and life-cycle duration. Field Crop. Res. 1989, 22, 211–225. [Google Scholar] [CrossRef] [Green Version]
  77. Singh, A.; Goyal, V.; Mishra, A.; Parihar, S. Performance evaluation of Crop Syst simulation model for pearlmillet (Pennisetum glaucum)-chickpea (Cicer arietinum) cropping system. Indian J. Agron. 2015, 60, 516–523. [Google Scholar]
  78. Van Heemst, H. Plant Data Values Required for Simple Crop Growth Simulation Models: Review and Bibliography; Centre for Agrobiological Research and Agricultural University: Wageningen, The Netherlands, 1988. [Google Scholar]
  79. Monteith, J.; Huda, A.; Midya, D. Modeling Sorghum and Pearl Millet. Modeling the Growth and Development of Sorghum and Pearl Millet; International Crops Research Institute for the Semi-Arid Tropics: Andhra Pradesh, India, 1989. [Google Scholar]
  80. Wet, J.; Bidinger, F.; Peacock, J. Pearl Millet (Pennisetum glaucum)—A Cereal of the Sahel; Academic Press: Cambridge, MA, USA, 1992. [Google Scholar]
  81. Klaij, M.; Vachaud, G. Seasonal water balance of a sandy soil in Niger cropped with pearl millet, based on profile moisture measurements. Agric. Water Manag. 1992, 21, 313–330. [Google Scholar] [CrossRef] [Green Version]
  82. Gregory, P.; Reddy, M. Root growth in an intercrop of pearl millet/groundnut. Field Crop. Res. 1982, 5, 241–252. [Google Scholar] [CrossRef] [Green Version]
  83. McIntyre, B.; Riha, S.; Flower, D. Water uptake by pearl millet in a semiarid environment. Field Crop. Res. 1995, 43, 67–76. [Google Scholar] [CrossRef] [Green Version]
  84. Singh, R.; Joshi, N.; Singh, H. Pearl millet phenology and growth in relation to thermal time under arid environment. J. Agron. Crop Sci. 1998, 180, 83–91. [Google Scholar] [CrossRef]
  85. Garcia-Huidobro, J.; Monteith, J.; Squire, G. Time, temperature and germination of pearl millet (Pennisetum typhoides S. & H.) I. Constant temperature. J. Exp. Bot. 1982, 33, 288–296. [Google Scholar]
  86. Teowolde, H.; Voigt, R.L.; Osman, M.; Dobrenz, A.K. Water Stress Indices for Research and Irrigation Scheduling in Pearl Millet; College of Agriculture, University of Arizona: Tucson, AZ, USA, 1987. [Google Scholar]
  87. Spitters, C. Crop growth models: Their usefulness and limitations. In Proceedings of the VI Symposium on the Timing of Field Production of Vegetables 267, Wageningen, The Netherlands, 21–25 August 1989; pp. 349–368. [Google Scholar]
  88. Yulin, L.; Johnson, D.A.; Yongzhong, S.; Jianyuan, C.; Zhang, T. Specific leaf area and leaf dry matter content of plants growing in sand dunes. Bot. Bull. Acad. Sin. 2005, 46, 127–134. [Google Scholar]
  89. Silungwe, F.R.; Graef, F.; Bellingrath-Kimura, S.D.; Tumbo, S.D.; Kahimba, F.C.; Lana, M.A. The management strategies of pearl millet farmers to cope with seasonal rainfall variability in a semi-arid agroclimate. Agronomy 2019, 9, 400. [Google Scholar] [CrossRef] [Green Version]
  90. Agrawal, S.; Jain, A. Sustainable deployment of solar irrigation pumps: Key determinants and strategies. Wiley Interdiscip. Rev. Energy Environ. 2019, 8, e325. [Google Scholar] [CrossRef]
  91. Sinclair, T.R.; Muchow, R.C. Radiation use efficiency. In Advances in Agronomy; Elsevier: Amsterdam, The Netherlands, 1999; Volume 65, pp. 215–265. [Google Scholar]
  92. Monteith, J.L. Validity of the correlation between intercepted radiation and biomass. Agric. For. Meteorol. 1994, 68, 213–220. [Google Scholar] [CrossRef]
  93. Carretero, R.; Serrago, R.A.; Bancal, M.O.; Perelló, A.E.; Miralles, D.J. Absorbed radiation and radiation use efficiency as affected by foliar diseases in relation to their vertical position into the canopy in wheat. Field Crop. Res. 2010, 116, 184–195. [Google Scholar] [CrossRef]
  94. Koester, R.P.; Skoneczka, J.A.; Cary, T.R.; Diers, B.W.; Ainsworth, E.A. Historical gains in soybean (Glycine max Merr.) seed yield are driven by linear increases in light interception, energy conversion, and partitioning efficiencies. J. Exp. Bot. 2014, 65, 3311–3321. [Google Scholar] [CrossRef]
  95. Teixeira, E.I.; George, M.; Herreman, T.; Brown, H.; Fletcher, A.; Chakwizira, E.; de Ruiter, J.; Maley, S.; Noble, A. The impact of water and nitrogen limitation on maize biomass and resource-use efficiencies for radiation, water and nitrogen. Field Crop. Res. 2014, 168, 109–118. [Google Scholar] [CrossRef]
  96. Arkebauer, T.; Weiss, A.; Sinclair, T.; Blum, A. In defense of radiation use efficiency: A response to Demetriades-Shah et al. (1992). Agric. For. Meteorol. 1994, 68, 221–227. [Google Scholar] [CrossRef]
  97. Green, C. Genotypic differences in the growth of Triticum aestivum in relation to absorbed solar radiation. Field Crop. Res. 1989, 19, 285–295. [Google Scholar] [CrossRef]
  98. Robertson, M.; Silim, S.; Chauhan, Y.; Ranganathan, R. Predicting growth and development of pigeonpea: Biomass accumulation and partitioning. Field Crop. Res. 2001, 70, 89–100. [Google Scholar] [CrossRef] [Green Version]
  99. Ullah, A.; Ahmad, A.; Khaliq, T.; Akhtar, J. Recognizing production options for pearl millet in Pakistan under changing climate scenarios. J. Integr. Agric. 2017, 16, 762–773. [Google Scholar] [CrossRef] [Green Version]
  100. McIntyre, B.; Flower, D.; Riha, S. Temperature and soil water status effects on radiation use and growth of pearl millet in a semi-arid environment. Agric. For. Meteorol. 1993, 66, 211–227. [Google Scholar] [CrossRef]
  101. Verdoodt, A.; Van Ranst, E.; Ye, L. Daily simulation of potential dry matter production of annual field crops in tropical environments. Agron. J. 2004, 96, 1739–1753. [Google Scholar] [CrossRef] [Green Version]
  102. Liu, J.; Fan, Y.; Ma, Y.; Li, Q. Response of photosynthetic active radiation interception, dry matter accumulation, and grain yield to tillage in two winter wheat genotypes. Arch. Agron. Soil Sci. 2020, 66, 1103–1114. [Google Scholar] [CrossRef]
  103. Jonckheere, I.; Fleck, S.; Nackaerts, K.; Muys, B.; Coppin, P.; Weiss, M.; Baret, F. Methods for leaf area index determination. Part I: Theories, techniques and instruments. Agric. For. Meteorol. 2004, 121, 19–35. [Google Scholar] [CrossRef]
  104. Huang, S.; Gao, Y.; Li, Y.; Xu, L.; Tao, H.; Wang, P. Influence of plant architecture on maize physiology and yield in the Heilonggang River valley. Crop J. 2017, 5, 52–62. [Google Scholar] [CrossRef] [Green Version]
  105. Timlin, D.; Kim, S.-H.; Fleisher, D.; Wang, Z.; Reddy, V.R. Maize water use and yield in the solar corridor system: A simulation study. In The Solar Corridor Crop System; Elsevier: Amsterdam, The Netherlands, 2019; pp. 57–78. [Google Scholar]
  106. Rinaldi, M. Variation of specific leaf area for sugar beet depending on sowing date and irrigation. Ital. J. Agron. 2003, 7, 23–32. [Google Scholar]
  107. Danalatos, N.; Mitsios, I.; Archontoulis, S. Growth and Biomass Productivity of Kenaf as Biomass Crop in Central Greece. Available online: http://www.cres.gr/biokenaf/files/fs_inferior01_h_files/pdf/articles/Kenaf,%20Danalatos%20et%20al.,%202005.pdf (accessed on 15 January 2022).
  108. Amanullah, M.J.H.; Nawab, K.; Ali, A. Response of specific leaf area (SLA), leaf area index (LAI) and leaf area ratio (LAR) of maize (Zea mays L.) to plant density, rate and timing of nitrogen application. World Appl. Sci. J. 2007, 2, 235–243. [Google Scholar]
  109. Wadas, W.; Kosterna, E. Effect of perforated foil and polypropylene fibre covers on assimilation leaf area of early potato cultivars. Plant Soil Environ. 2007, 53, 299. [Google Scholar] [CrossRef] [Green Version]
  110. Milla, R.; Reich, P.B.; Niinemets, Ü.; Castro-Díez, P. Environmental and developmental controls on specific leaf area are little modified by leaf allometry. Funct. Ecol. 2008, 22, 565–576. [Google Scholar] [CrossRef]
  111. Reich, P.B.; Walters, M.B.; Ellsworth, D.S. From tropics to tundra: Global convergence in plant functioning. Proc. Natl. Acad. Sci. USA 1997, 94, 13730–13734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  112. Field, C.B. On the role of photosynthetic responses in constraining the habitat distribution of rainforest plants. Functl. Plant Biol. 1988, 15, 343–358. [Google Scholar]
  113. Chapin III, F.S.; Autumn, K.; Pugnaire, F. Evolution of suites of traits in response to environmental stress. Am. Nat. 1993, 142, S78–S92. [Google Scholar] [CrossRef]
  114. Mooney, H.; Dunn, E. Photosynthetic systems of Mediterranean-climate shrubs and trees of California and Chile. Am. Nat. 1970, 104, 447–453. [Google Scholar] [CrossRef]
  115. Szulc, P.; Piechota, T.; Jagła, M.; Kowalski, M. A comparative analysis of growth in maize (Zea mays L.) hybrids of different genetic profiles depending on type of nitrogen fertilizer and magnesium dose. Commun. Biometry Crop Sci. 2015, 10, 73–81. [Google Scholar]
  116. Castro-Díez, P.; Puyravaud, J.-P.; Cornelissen, J. Leaf structure and anatomy as related to leaf mass per area variation in seedlings of a wide range of woody plant species and types. Oecologia 2000, 124, 476–486. [Google Scholar] [CrossRef] [PubMed]
  117. Federer, C.A. Measuring Forest Evapotranspiration: Theory and Problems; Northeastern Forest Experiment Station, Forest Service, U.S. Department of Agriculture: Washington, DC, USA, 1970. [Google Scholar]
  118. Ward, R. Measuring evapotranspiration; a review. J. Hydrol. 1971, 13, 1–21. [Google Scholar] [CrossRef]
  119. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper 56. Fao Rome 1998, 300, D05109. [Google Scholar]
  120. Hillel, D. Environmental Soil Physics: Fundamentals, Applications, and Environmental Considerations; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  121. Jackson, R.D.; Idso, S.; Reginato, R.; Pinter, P., Jr. Canopy temperature as a crop water stress indicator. Water Resour. Res. 1981, 17, 1133–1138. [Google Scholar] [CrossRef]
  122. Jackson, R.D.; Kustas, W.P.; Choudhury, B.J. A reexamination of the crop water stress index. Irrig. Sci. 1988, 9, 309–317. [Google Scholar] [CrossRef]
  123. Kacira, M.; Ling, P.P.; Short, T.H. Establishing Crop Water Stress Index (CWSI) threshold values for early, non-contact detection of plant water stress. Trans. ASAE 2002, 45, 775. [Google Scholar] [CrossRef]
  124. Uçak, A.B. Identification of Water Usage Efficiency for Corn (Zea mays L.) Lines Irrigated with Drip Irrigation under Green House Conditions as Per Plant Water Stress Index Evaluations. Turk. Tarımsal Arastırmalar Derg. Turk. J. Agric. Res. 2017, 4, 1–9. [Google Scholar] [CrossRef] [Green Version]
  125. Monteith, J.L. Climate and the efficiency of crop production in Britain. Philos. Trans. R. Soc. B Biol. Sci. 1977, 281, 277–294. [Google Scholar]
  126. Campbell, G.S.; Norman, J.M. The light environment of plant canopies. In An Introduction to Environmental Biophysics; Springer: Berlin/Heidelberg, Germany, 1998; pp. 247–278. [Google Scholar]
  127. Tanner, C.; Sinclair, T. Efficient water use in crop production: Research or research? In Limitations to Efficient Water Use in Crop Production; Wiley: Hoboken, NJ, USA, 1983; pp. 1–27. [Google Scholar]
  128. Campbell, G.; Van Evert, F. Light interception by plant canopies: Efficiency and architecture. In Resource Capture by Crops; Nottingham University Press: Nottingham, UK, 1994; pp. 35–52. [Google Scholar]
  129. Ritchie, J.T. Model for predicting evaporation from a row crop with incomplete cover. Water Resour. Res. 1972, 8, 1204–1213. [Google Scholar] [CrossRef] [Green Version]
  130. Goudriaan, J. Crop Micrometeorology: A Simulation Study; Wageningen University and Research: Wageningen, The Netherlands, 1977. [Google Scholar]
  131. Yadav, O.; Bidinger, F. Dual-purpose landraces of pearl millet (Pennisetum glaucum) as sources of high stover and grain yield for arid zone environments. Plant Genet. Resourc. 2008, 6, 73. [Google Scholar] [CrossRef]
Figure 1. Cumulative top dry matter (TDM) production of the three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid as a function of the cumulative (∑) product of fractional interception and solar radiation (FI × Rs). Radiation use efficiency (RUE) and the coefficient of determination (R2) are shown. FI values were determined on the day of measurement.
Figure 1. Cumulative top dry matter (TDM) production of the three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid as a function of the cumulative (∑) product of fractional interception and solar radiation (FI × Rs). Radiation use efficiency (RUE) and the coefficient of determination (R2) are shown. FI values were determined on the day of measurement.
Plants 11 00806 g001
Figure 2. The relationship between leaf area index (LAI) and fractional interception (FI) of photosynthetically active radiation (PAR) for three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid. The coefficient of determination (R2), as well as the canopy extinction coefficient for photosynthetically active radiation (KPAR) and total solar radiation (Ks) are shown.
Figure 2. The relationship between leaf area index (LAI) and fractional interception (FI) of photosynthetically active radiation (PAR) for three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid. The coefficient of determination (R2), as well as the canopy extinction coefficient for photosynthetically active radiation (KPAR) and total solar radiation (Ks) are shown.
Plants 11 00806 g002
Figure 3. Determination of the stem-leaf dry matter partitioning parameter (p) as a function of canopy dry matter (CDM), specific leaf area (SLA) and leaf area index (LAI) for pearl millet varieties: (a) landrace, (b) improved and (c) hybrid, with the slope of the regression line (p) and the coefficient of determination (R2).
Figure 3. Determination of the stem-leaf dry matter partitioning parameter (p) as a function of canopy dry matter (CDM), specific leaf area (SLA) and leaf area index (LAI) for pearl millet varieties: (a) landrace, (b) improved and (c) hybrid, with the slope of the regression line (p) and the coefficient of determination (R2).
Plants 11 00806 g003
Figure 4. Measured values of specific leaf area (SLA) during the growing season of three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid, the value shown on horizontal line is the average SLA for each variety.
Figure 4. Measured values of specific leaf area (SLA) during the growing season of three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid, the value shown on horizontal line is the average SLA for each variety.
Plants 11 00806 g004
Figure 5. Simulated (solid lines) and measured values (points) of soil water deficits of three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid the value shown on horizontal line is the average total water deficit.
Figure 5. Simulated (solid lines) and measured values (points) of soil water deficits of three pearl millet varieties: (a) landrace, (b) improved and (c) hybrid the value shown on horizontal line is the average total water deficit.
Plants 11 00806 g005
Figure 6. Simulated (solid lines) and measured (points) leaf area index. Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly excluding (a).
Figure 6. Simulated (solid lines) and measured (points) leaf area index. Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly excluding (a).
Plants 11 00806 g006
Figure 7. Simulated (solid lines) and measured (points) harvestable dry matter (HDM) and above-ground dry matter (AGDM). Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly every season excluding (a).
Figure 7. Simulated (solid lines) and measured (points) harvestable dry matter (HDM) and above-ground dry matter (AGDM). Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly every season excluding (a).
Plants 11 00806 g007
Figure 8. Simulated (solid lines) and measured (points) soil profile deficit to field capacity. Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly every season excluding (a).
Figure 8. Simulated (solid lines) and measured (points) soil profile deficit to field capacity. Landrace (Kantana): (a) 2018 rainfed, (b) 2017, (c) 2018. Improved (Kangara): (d) 2016, (e) 2017, (f) 2018. Hybrid (Agrigreen): (g) 2016, (h) 2017, (i) 2018. All three varieties were irrigated fortnightly every season excluding (a).
Plants 11 00806 g008
Table 1. Weekly soil water balance components (mm) for three pearl millet varieties (December–April of 2017/2018). Precipitation (P) and irrigation (I), evapotranspiration (ET), drainage (D), runoff (R), amount of water contents at the beginning of a measurement period in the root zone (Qo), amount of water contents at the end of a measurement period in the root zone (Qi) and change in water storage (ΔS).
Table 1. Weekly soil water balance components (mm) for three pearl millet varieties (December–April of 2017/2018). Precipitation (P) and irrigation (I), evapotranspiration (ET), drainage (D), runoff (R), amount of water contents at the beginning of a measurement period in the root zone (Qo), amount of water contents at the end of a measurement period in the root zone (Qi) and change in water storage (ΔS).
Landrace Pearl Millet (Kantana)Improved Pearl Millet (Kangara)Hybrid (Agrigreen)
DatePIETDRQiQo∆SIETDRQiQoΔSIETDRQiQoΔS
14 December 20172601300315328−1401400264276−1301200292306−14
22 December 2017280700328350−2128800276324−480700306327−21
29 December 20171501300350352−241200324332−8018103273224
5 January 20187102900352340128230033232480381032229032
12 January 201811235000340323173147003243195405000290291−1
19 January 201826465500323339−17405400319332−12515400291314−23
31 January 2018343981003393318728100332356−25718180314330−16
5 February 201812374500331336−4503600356383−27603640330363−32
18 February 2018794687280336345−10087190383355274289280363366−3
5 March 20181560816034533412818170355363−86684603663579
10 March 20183313110334336−2633100363398−34523110357380−23
16 March 20182594010336357−21604010398420−22523910380395−14
24 March 2018209044880357434−77044640420521−101045510395508−113
2 April 2018838403404344062904153052143586244129050847038
8 April 2018116283040639114 0283047044030
14 April 201864025170391413−22
Table 2. Seasonal soil water balance (mm) for three pearl millet varieties (December–April of 2017/2018). Precipitation (P) and irrigation (I), evapotranspiration (ET), drainage (D), runoff (R), evapotranspiration (ETo) and average vapour pressure deficit is in Pascal (Pa) (Avg. VPD).
Table 2. Seasonal soil water balance (mm) for three pearl millet varieties (December–April of 2017/2018). Precipitation (P) and irrigation (I), evapotranspiration (ET), drainage (D), runoff (R), evapotranspiration (ETo) and average vapour pressure deficit is in Pascal (Pa) (Avg. VPD).
Landrace Pearl Millet (Kantana)Improved Pearl Millet (Kantana)Hybrid (Agrigreen)
PIETDREToAvgVPDIETDREToAvgVPDIETDREToAvgVPD
54240567017904360.7843759814404120.8345865313404250.81
Table 4. Treatments in 2016, 2017 and 2017/2018 varieties (V1, V2, V3) and water irrigation regimes (I0, I1, I2).
Table 4. Treatments in 2016, 2017 and 2017/2018 varieties (V1, V2, V3) and water irrigation regimes (I0, I1, I2).
Pearl Millet VarietiesIrrigation Regime
I0I1I2
V1012
V2012
V3012
Note: V1 = hybrid (Agrigreen), V2 = open pollinated landrace (Kantana), V3 = open pollinated improved (Kangara).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ausiku, P.A.; Annandale, J.G.; Steyn, J.M.; Sanewe, A.J. Crop Model Parameterisation of Three Important Pearl Millet Varieties for Improved Water Use and Yield Estimation. Plants 2022, 11, 806. https://doi.org/10.3390/plants11060806

AMA Style

Ausiku PA, Annandale JG, Steyn JM, Sanewe AJ. Crop Model Parameterisation of Three Important Pearl Millet Varieties for Improved Water Use and Yield Estimation. Plants. 2022; 11(6):806. https://doi.org/10.3390/plants11060806

Chicago/Turabian Style

Ausiku, Petrus A., John G. Annandale, Joachim Martin Steyn, and Andrew J. Sanewe. 2022. "Crop Model Parameterisation of Three Important Pearl Millet Varieties for Improved Water Use and Yield Estimation" Plants 11, no. 6: 806. https://doi.org/10.3390/plants11060806

APA Style

Ausiku, P. A., Annandale, J. G., Steyn, J. M., & Sanewe, A. J. (2022). Crop Model Parameterisation of Three Important Pearl Millet Varieties for Improved Water Use and Yield Estimation. Plants, 11(6), 806. https://doi.org/10.3390/plants11060806

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