Investigating the Pearl Millet ( Pennisetum glaucum ) as a Climate-Smart Drought-Tolerant Crop under Jordanian Arid Environments

: To investigate drought tolerance under arid conditions, eleven pearl millet breeds (HHVBC tall B6; IP13150; IP19586; IP19612; IP22269; IP6110; IP7704; MC94C2; P. millet icms7709; Sudan-pop I; Sudan-pop III) were tested under arid water-scarce climatic conditions. A ﬁeld randomized complete block design experiment with three replicates per year was conducted at the Deir-Alla Regional Agriculture Research center in the middle Ghor within the Jordan Valley from 2010 to 2020. The plant-deﬁcit irrigation was maintained at 80% based on the crop water requirements using a time-domain reﬂectometer. The plant morphological characteristics, forage production, seed formation, and water-use efﬁciency (WUE) were monitored for ten years for two case scenarios: seed and forage production. The individual and combined drought indices of the precipitation, temperature, and vegetation were calculated and correlated with the millet morphological and yield parameters. Climate change analyses show signiﬁcant impacts, reaching a 1 mm/year reduction in precipitation and a 0.04 mm/year increase in air temperature, which causes the study area to be more prone to drought events. Along with the proven increase in the drought intensity over time, the millet breeds showed signiﬁcant drought-tolerance capacities under arid, drought-prone conditions by adjusting their system to tolerate salt, heat, and water stresses. For the seed production scenario, the WUE ranged from 27 to 57.3 kg/ha · mm, and from 7.1 to 14.9 kg/ha · mm for fresh and dry conditions, respectively. The IP13150 millet breed showed the highest capacity to tolerate the drought of Jordan’s environment, and it is thus recommended as a good substitute under water-scarcity situations, with an average production of 17.7 ton/ha. For the vegetative production scenario, the WUE ranged from 32.03 to 64.82 kg/ha · mm for the fresh biomass and from 10.8 to 24.6 kg/ha · mm for the dry biomass. Based on the WUEs and vegetative production results, the IP19586, IP22269, IP19612, IP7704, and HHVBC tall B6 millet breeds are recommended as forage support due to their phenological characteristics, which tolerate drought and heat conditions. In contrast to the vegetation drought index, both the precipitation and temperature drought indices show strong correlations (above r > 0.6) with the plant growth factors and a moderate correlation (0.3 < r < 0.6) with the yield factors. Both precipitation and temperature indices are capable of explaining the variations among millet breeds, especially as related to millets’ morpho-physiological characteristics.


Introduction
Along with rapid population growth, the increasing energy-food-water demand, and sudden refugee migrations, drought remains the most complex natural phenomenon threatening the sustainable development of all sectors, including the economy, agriculture, water, ecosystems, health, tourism, and even society, at the global, regional, and local levels [1,2]. Over the last two decades, many researchers have pointed to the increase in a tendency to tiller copiously and can compensate for uneven stand establishment. To support the developing plants, prop roots emerge from the lower nodes [32][33][34]. The Pearl Millet kind of millet has the biggest kernels overall. The color of the kernel can range from white to grey to pale yellow to slate blue to purple to brown. The plant has a four-meter height limit [35,36]. The amount of wheat, millet, canary seed, and other grains that Jordan imported from Argentina in 2018 was 15.13 thousand US dollars, according to the UN Comtrade Database on International Trade [37]. Due to the lack of improved variations, the introduction of millet into the plant species in marginal areas, and its usage as a rainfed crop, because it is a hardy plant, millet yields were seen to be lower in many nations.
Many studies are currently attempting to introduce new species with high productivity through the use of modified kinds that can multiply green fodder and, as a result, might feed twice as many animals per unit area as conventional fodder crops [38]. International research has been undertaken through a number of studies to examine the water-use efficiency of pearl millet in water-limited conditions [39,40]. According to certain research [30,41,42], a plant's rooting behavior impacts its capacity to use the soil's moisture reserves and how much salinity it can endure. Others have looked into the production of pearl millet in dry conditions and scenarios with changing climatic variables [43][44][45][46][47]. Pearl millet has been proposed as a Climate-Resilient Nutri-cereal for reducing covert hunger and ensuring nutritional security [44,[48][49][50][51][52]. Assessing the impact of drought on pearl millet's morphological and yield responses in dry environments would help guide the selection of prospective drought-tolerant pearl millet cultivars for drought adaptation choices. This research is justified by the actual evaluation of the diversity and effectiveness of many pearl millet types utilizing lengthy historical drought indices. The primary goal of this study was to link several drought indicators with the growth and yield characteristics of various types of pearl millet. The key claim is that drought indicators can identify and track Pearl Millet's ability for drought tolerance by inferring temporal plant phonological alterations under heat and drought loads.

Site Description
An experiment was carried out at the Deir-Alla Regional Agriculture Research Centre in the Jordan Valley's middle Ghor, 50 km west of Amman. The Centre is located at 32 • 13 E, 35 • 37 N, and 224 feet below sea level ( Figure 1). The climate in the area is semi-arid, with warm winters and hot summers and an average annual rainfall of 280 mm. The mean maximum annual temperature is 30 degrees Celsius centigrade, and the mean minimum annual temperature is 15 degrees Celsius centigrade, which sometimes goes to a few degrees in harsh winters. The mean maximum annual temperature is 30 degrees Celsius centigrade, and the mean minimum annual temperature is 15 degrees Celsius centigrade, which sometimes goes to a few degrees in harsh winters. The mean relative humidity ranges from 30% in summer to 70% in winter. This region is highly prone to drought with unevenly distributed rainfall. The long-term annual potential evaporation at the site is about 2200 mm.
The soil is located within the Ghor map unit representing a fine, mixed, hyper-thermic, deep family of Typic Ustochrepts. It is characterized by dark brown and dark yellowish brown (7.5 YR-10YR 4\4) deep (>80 cm) clay loam and light clay soils with weakly cracked surfaces and a strong medium sub-angular block structure. The soil is highly calcareous and non-saline, with slope gradients of <2%. More information on the soil is presented in Table 1. The available water resources at the site are the convey pressurized pipes from the King Talal Reservoir (KTR) at Zerqa River with moderate salinity ranging from 1.4 dS/m in winter to about 3.0 dS/m in summer, which may negatively impact the salt accumulation in soil profile and deterioration of land productivity. Due to climate change impacts, salinity is increasing over time, making potential deterioration more likely. The soil is located within the Ghor map unit representing a fine, mixed, hyper mic, deep family of Typic Ustochrepts. It is characterized by dark brown and dar lowish brown (7.5 YR-10YR 4\4) deep (>80 cm) clay loam and light clay soils with w cracked surfaces and a strong medium sub-angular block structure. The soil is high careous and non-saline, with slope gradients of <2%. More information on the soil sented in Table 1. The available water resources at the site are the convey pressu pipes from the King Talal Reservoir (KTR) at Zerqa River with moderate salinity ra from 1.4 dS/m in winter to about 3.0 dS/m in summer, which may negatively impa salt accumulation in soil profile and deterioration of land productivity. Due to c change impacts, salinity is increasing over time, making potential deterioration likely.

Drought Monitoring
Thirty years of daily climate data (1980-2020) were obtained from in-situ weather station. DrinC software version 1.5.73 was used to calculate the Drought Standardized Precipitation Index (SPI) [53]. The long-term monthly and seasonal rainfall records were initially normalized using the gamma distribution function. The temporal SPI magnitudes were calculated by dividing the difference between the normalized seasonal precipitation and its long-term seasonal precipitation mean by the standard deviation (Equation (1)).
where Xi is the seasonal precipitation at the rain gauge station, Xm is the long-term seasonal mean (i.e., 30 years of records), and where σ is its standard deviation. According to McKee et al., a drought event occurs any time the SPI is continuously negative and reaches an intensity of −1.0 or less [54]. The event ended as SPI became positive. In this study, drought severity was divided into seven classes: Extremely wet (SPI > 2), very wet (1.5 to 1.99), moderately wet (1.0 to 1.49), near normal (−0.99 to 0.99), moderate drought (−1.49 to −1), severe drought (−1.99 to −1.5), and extreme drought (SPI < −2). In addition, the Combined Drought Index (CDI) was computed as a linear weighted average from the Precipitation Drought Index (PDI), Temperature Drought Index, and Vegetation Drought Index (PDI) (Equation (2), Equation (3), and Equation (4), respectively). Based on Balint et al. recommendations, the weights were assigned as 50% for PDI, 25% for each TDI and VDI [55].
where P * is the modified annual precipitation, T * is the modified annual temperature, NDV I * is the modified annual average Normalized Difference Vegetation index, IP is the interest period, RRL m,i (run-length) as the maximum number of successive years below long-term average rainfall in the intended period, RL maximum number of successive years below long-term average NDVI in the IP, which indicates the number of years with relevant data, j is a summation running parameter covering the IP, and k is the summation parameter covering the years for which relevant data are available. The modified temperature, modified NDVI, and modified rainfall data were obtained using (Equations (5)-(8)) to avoid dividing by zero in certain cases, as rainfall in Jordan is mostly characterized by a distinct long dry season; this also helped to unify the ranges of the drought index values: NDV I * = NDV I − (NDV I min − 0.01) (7) where P, T, and NDVI are the original precipitation, temperature, and NDVI values and RL is the original run-length. The Normalized Difference Vegetation Index (NDVI) values were computed according to Kogan [56,57] using the ratio of responses in the near-infrared (NIR) and visible red portion of the spectrum (R) bands of the Advanced Very High-Resolution Radiometer (AVHRR) at the National Oceanic and Atmospheric Administration (NOAA) (Equation (9)).
The extracted monthly NDVI data from the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) Vegetation Indices (MOD13A3) Version 6 data at a 1 km spatial resolution [58]). A CDI of 1.0 thus represents average weather conditions; if the CDI is greater than 1.0, this represents wetter than average conditions, and if it is below 1.0, this represents dryer than average conditions. According to Balint et al. [55], five drought categories were implemented in this study, as presented in Table 2.

Crop Experiment
Eleven pearl millet breeds, including HHVBC tall B6, IP13150, IP19586, IP19612, IP22269, IP6110, IP7704, MC94C2, P. millet icms7709, Sudan pop I, and Sudan pop III, were cultivated at Deir-Alla Regional Agriculture Research Center from 2010 to 2020. Twelve hectares were used to set the field trials with a Complete Randomized Design (CRD) experiment, including three replications. The land seedbed preparation imitated the common practices adopted in the region and included chisel plowing followed by disk smothering for the top 20 cm. The seeds were planted in October and harvested at maturity by the end of February. The seeds were spaced 40 cm within rows and 40 cm between rows. Chemical fertilizers were applied prior to the plantation by 250 kg/ha triple superphosphate (46% P 2 O 5 ) and 200 kg ha −1 granular urea (46%N). At 30 to 40 cm plant height, another 200 kg ha −1 of urea was applied through fertigation practice.

Crop Quality, Productivity, and Water-Use Efficiency
Two scenarios were used to evaluate the crop productivity and quality of millet breeds: Seed Production and Forage Production Scenarios. For the seed production scenario, the plants were solely tested under rainfed conditions. After maturity, the fruit length (cm), fruit diameter (cm), root depth (cm), mature fresh yield (ton/ha), mature dry yield (ton/ha), and seed yield (ton/ha) were measured per breed. A destructive 1.5 m × 1.5 m Quadrat sampling method was used to collect the plant samples after maturing. Random selection was used for plant aboveground sampling from each plot. The plant samples were chopped, fresh weighed, oven dried at 80 • C for 24 h, and dry weighed. Furthermore, the millet seeds were gathered by hand, fresh weighed, dried using a continuous flow drier for 48 h, and dry weighed. To investigate plants' reproductive efficiency, especially under heat and drought stresses, the Harvest Index (HI) was calculated as the ratio of grain yield to total aboveground biomass [59].
For the vegetative production scenario, all of the breeds underwent three cuts to investigate the progressive vegetative growth capacity (i.e., fresh forage at different cutting). This experiment was conducted under supplemental irrigations, in which the plants were irrigated using a drip system fitted by inline emitters of 4 l ph that are spaced 40 cm meters and 40 cm between lateral lines. The irrigation schedule was achieved once every ten days based on crop water demand.
The crop water requirements were estimated using the water balance equation (Equation (10)). Revised Penman-Monteith equation used to calculate the potential evapotranspiration (Equation (11)).
where S is the soil water stored in the root zone, P is the precipitation, I is the irrigation amount, R is the runoff, D is the drainage, and ET o is the potential evapotranspiration [mm/day], Rn is the net radiation at the crop surface (MJ/m 2 ·day), G is the soil heat flux density (MJ/m 2 ·day), T is the mean daily air temperature at 2-m height ( • C), U 2 is the wind speed at 2-m height (m/s), e s is the saturation vapor pressure (kPa), e a is the actual vapor pressure (kPa), e s − e a is the saturation vapor pressure deficit (kPa), ∆ is the slope vapor pressure curve (kPa/ • C), and γ is the psychrometric constant (kPa/ • C). The actual crop evapotranspiration (ET c ) is derived by multiplying the potential evapotranspiration (ET o ) by the crop factor (K c ) based on both growth stages and soil water measurements. Daily soil water content for the five soil increments (0-20, 20-40, 40-80, 80-120, and 120-160 cm) was measured at each site using the time domain reflectometer "Trime-FM3 TDR" (Imko GmbH, Ettlingen, Germany). Supplemental irrigation events were achieved regarding soil water deficit and crop water requirements. When the soil reaches the critical stage, irrigation was attained until reaching 80% of the available water. Throughout the temporal study, the amount applied varies from plot to plot and thus from replicate to replicate. The irrigation depth is the sum of the effective rainfall and irrigation minus the drainage water.
Crop measurements for the vegetative production scenario included the number of branches, plant height (cm), plant diameter (cm), fresh weight (ton/ha), and dry weight (ton/ha) per cut. The water-use efficiencies (WUE) were estimated by dividing the biomass weights (fresh biomasses and dry biomasses) per amount of water received (m 3 ), including the effective rainfall and supplemental water applied.

Investigating the Impact of Drought on Plant Growth
A Pearson Correlation Matrix was generated between the drought indices and millet production to investigate the drought indices' capacities and to represent the temporal variations among millet breeds (Equation (12)).
where r xy is the correlation coefficient, x i is the drought indicator at each specific year, x is the mean of the drought indicator across all the years of the study, y i is the millet production (ton/ha) for each year, and y x is the millet production mean across all the years of the study. A comparison between means for all crop parameters and among millet breeds was achieved using Tukey-Kramer test within JMP statistical software version 11.0.0 [60]. The Tukey-Kramer test represents an honest test at a 95% confidence level.

Climate Temporal Variability
Recorded climate variables during the experiment period show temporal dispersity as indicated by the associated high standard deviations and the coefficient of variation, especially for precipitation ( Table 3). The annual precipitation varies significantly from 99.8 mm during the dry period to 500 mm during the wet season ( Figure 2).
Precipitation decreases insignificantly over time at a rate of 1 mm per year (p-value of 0.4732), thus suggesting a future reduction of about 25% by the end of the 21st century. Similarly, the annual and inter-annual air temperature varies significantly, ranging from 1.8 to a maximum of 50.1 • C ( Figure 3). The trend of temporal change indicates a significant increase (p-value > 0.0001) by a rate of 0.04 • C per year, thus increasing up to 3.2 • C by the end of the 21st century.

Drought Intensity and Severity
All of the drought indicators confirm that the study area is highly prone to drought events. As indicated in Table 4, the precipitation drought indicator is the highest dispersed variable compared to temperature and vegetation drought indicators. This illustrates how the study area is subject to significant climate changes, especially in terms of precipitation followed by temperature.
The SPI trend in Figure 4 shows a frequent alternation from extremely wet (SPI > 2) to severe drought (−1.99 to −1.5) with the exception of lack of extreme drought (SPI < −2). Based on the SPI results, 14.6% of the period between the years 1980 to 2020 is classified as drought (SPI < −1), while 70.7% of the records are classified as near normal, and 14.6% are classified as wet. Local drought events' occurrence frequency tends to reach once every seven years, which persists from one to two consecutive years.

Drought Intensity and Severity
All of the drought indicators confirm that the study area is highly prone to d events. As indicated in Table 4, the precipitation drought indicator is the high persed variable compared to temperature and vegetation drought indicators. Th trates how the study area is subject to significant climate changes, especially in t precipitation followed by temperature. The SPI trend in Figure 4 shows a frequent alternation from extremely wet ( to severe drought (−1.99 to −1.5) with the exception of lack of extreme drought (SP Based on the SPI results, 14.6% of the period between the years 1980 to 2020 is cl as drought (SPI < −1), while 70.7% of the records are classified as near normal, and are classified as wet. Local drought events' occurrence frequency tends to reac every seven years, which persists from one to two consecutive years.  In terms of drought timing, the PDI always agrees with the SPI based on the temporal variation in the drought indicators ( Figure 5). The highly severe precipitation droughts occurred mainly in 1995, 1999, and 2017, while severe precipitation droughts occurred in 1993, 2008, and 2010. On the other hand, the TDI temporal trend indicates that 2001 and 2010 were categorized as extreme temperature droughts, while 1999, 2007, and 2016 were categorized as severe temperature droughts. Finally, the VDI temporal trend shows that the year 2013 can be categorized as an extreme drought, while the years 1985 and 2017 can be categorized as severe drought events (Table 5). Although the PDI, TDI, and VDI show different magnitudes of temporal changes, they all have similar linear trends of increasing intensity and frequency over time. Figure 6 shows the increasing CDI rate is about 0.009 magnitude per year, thus shifting the drought mean by one level towards the severe category by 2050 and towards extreme droughts by 2100. This illustrates the severe threats of future climate impacts in the study area.  In terms of drought timing, the PDI always agrees with the SPI based on the temporal variation in the drought indicators ( Figure 5). The highly severe precipitation droughts occurred mainly in 1995, 1999, and 2017, while severe precipitation droughts occurred in 1993, 2008, and 2010. On the other hand, the TDI temporal trend indicates that 2001 and 2010 were categorized as extreme temperature droughts, while 1999, 2007, and 2016 were categorized as severe temperature droughts. Finally, the VDI temporal trend shows that the year 2013 can be categorized as an extreme drought, while the years 1985 and 2017 can be categorized as severe drought events (Table 5). Although the PDI, TDI, and VDI show different magnitudes of temporal changes, they all have similar linear trends of increasing intensity and frequency over time. Figure 6 shows the increasing CDI rate is about 0.009 magnitude per year, thus shifting the drought mean by one level towards the severe category by 2050 and towards extreme droughts by 2100. This illustrates the severe threats of future climate impacts in the study area. Sustainability 2022, 14, x FOR PEER REVIEW 12 of 23

Irrigation Water Variability
Based on the soil water budget determined using both TDR and crop needs, the amount of supplemental water applied is significantly different across the millet breeds ( Figure 6) and across the study period (Figure 7), ranging from a minimum of 275 mm for a wet year to 365 mm in the driest year. On the other hand, some millet breeds display high water demands, as IP6110 and Sudan pop I regardless of the season (Table 6). Although the statistical analyses show no significant differences between water applied between breeds; however, the IP19586 and IP13150 millet breeds show low potential for water consumption and temporal variability as compared to other breeds.

Irrigation Water Variability
Based on the soil water budget determined using both TDR and crop needs, the amount of supplemental water applied is significantly different across the millet breeds ( Figure 6) and across the study period (Figure 7), ranging from a minimum of 275 mm for a wet year to 365 mm in the driest year. On the other hand, some millet breeds display high water demands, as IP6110 and Sudan pop I regardless of the season (Table 6). Although the statistical analyses show no significant differences between water applied between breeds; however, the IP19586 and IP13150 millet breeds show low potential for water consumption and temporal variability as compared to other breeds.

Comparison of General Parameters among Millet Breeds
At a 95% confidence level, the comparison of the means for the number of branches, plant height, fruit length, fruit diameter, and root depth reveal significant differences. (p < 0.001) regarding millet breeds and seasons. The number of branches ranges from 16 for the IP22269 breed to only seven branches for the IP19612 millet breed (Table 6). Similarly, the recorded millet height ranges from 113 cm for the IP22269 breed to 271 cm for the IP7704 breed.
In terms of fruit length and diameter, the IP13150 and HHVBC tall B6 millet breeds have the best records, regardless of seasonal variations, with an average of 26 cm and 3 cm, respectively. Seasonal variations are caused by drought-related conditions, while breed variations are caused by phenology and physiology. Some pearl millet breeds, such as IP13150 and HHVBC tall B6, showed high records and low dispersion in terms of seasonal variations, indicating that they can withstand drought conditions. These breeds can adjust their phenology to rainfall patterns and are thus unaffected by drought conditions [61][62][63].
The root depth records show significant differences between breeds and seasons. The root depth varies from 120 cm for the IP22269 breed to 183 cm for the IP7704 breed. The inter-seasonal variation among pearl millet breeds, on the other hand, is very low, implying that pearl millet breeds can adjust their root system under drought conditions. This coincides with Ajithkumar et al., in which pearl millet breeds can adjust their root system to become longer under drought conditions [41]. Pearl millet produces deep and profuse root systems that can withstand drought effects by extracting water from the lower soil profile. On the other hand, Zegada-Lizarazu and Iijima indicated that higher WUE could explain the drought resistance of pearl millet [64].

Comparison of Millet Production and WUE among Millet Breeds
According to the first scenario, millet production varies greatly and significantly by breed and season. The average maturity fresh yield ranges from 8.3 to 17.7 ton/ha, with the maximum associated with the IP13150 millet breed (Table 7). Furthermore, the average dry yield ranges from 2.2 to 4.6 ton/ha, with the maximum associated with the IP13150 millet breed. The WUE ranges from 27 to 57.3 kg/ha·mm and from 7.1 to 14.9 kg/ha·mm for fresh and dry conditions, respectively. These values coincide with many researchers' findings for millet production under drip irrigation practices [14,40,65,66]. Despite temporal variations in millet production, as shown in Figure 8, the promising millet breed IP13150 could potentially be a good substitute in water-stress situations. ICBA reached similar conclusions, reporting that pearl millet salt-tolerant hybrids produced 11% more grain and 38% more fodder than the most widely cultivated dual-purpose commercial hybrids. [67]. On the other hand, P. millet icms7709 and Sudan pop I breeds are unable to produce enough yield per unit of limited water, and thus cannot tolerate the arid environments of Jordan. According to the second scenario, the total fresh and dry biomasses have increased as fodder production has increased due to cutting practices. This suggests that millet breed variability toward cutting is due to the physiological behavior of vegetative growth rather than the surrounding environment. The total yields of millet breeds ranged from 9.9 to 20.2 tons per acre for fresh yield and from 3.3 to 7.7 tons per acre for dry yield. The average WUE of millet breeds in the second scenario ranges from 32.03 to 64.82 kg/ha·mm for fresh biomass and 10.8 to 24.6 kg/ha·mm for dry biomass.
In contrast to the first scenario, millet breeds can show good forage production (i.e., animal feed) as evidenced by insignificant variation in fresh and dry bio-mass and WUE, particularly for the IP19586, IP22269, IP19612, IP7704, and HHVBC tall B6 millet breeds. However, there are differences in the ratio of dry to fresh yield, which ranges from 0.31 at the MC94C2 breed to 0.53 at the HVBC tall B6 millet breed, despite the fact that the season varies greatly (for example, dry conditions). Figures 9 and 10 show that the millet breeds IP19586, IP22269, and IP19612 are significant good producers that can be used for forage production in Jordan's arid environment.  According to the second scenario, the total fresh and dry biomasses have increased as fodder production has increased due to cutting practices. This suggests that millet breed variability toward cutting is due to the physiological behavior of vegetative growth rather than the surrounding environment. The total yields of millet breeds ranged from 9.9 to 20.2 tons per acre for fresh yield and from 3.3 to 7.7 tons per acre for dry yield. The average WUE of millet breeds in the second scenario ranges from 32.03 to 64.82 kg/ha·mm for fresh biomass and 10.8 to 24.6 kg/ha·mm for dry biomass.
In contrast to the first scenario, millet breeds can show good forage production (i.e., animal feed) as evidenced by insignificant variation in fresh and dry bio-mass and WUE, particularly for the IP19586, IP22269, IP19612, IP7704, and HHVBC tall B6 millet breeds. However, there are differences in the ratio of dry to fresh yield, which ranges from 0.31 at the MC94C2 breed to 0.53 at the HVBC tall B6 millet breed, despite the fact that the season varies greatly (for example, dry conditions). Figures 9 and 10 show that the millet breeds IP19586, IP22269, and IP19612 are significant good producers that can be used for forage production in Jordan's arid environment. According to the second scenario, the total fresh and dry biomasses have increased as fodder production has increased due to cutting practices. This suggests that millet breed variability toward cutting is due to the physiological behavior of vegetative growth rather than the surrounding environment. The total yields of millet breeds ranged from 9.9 to 20.2 tons per acre for fresh yield and from 3.3 to 7.7 tons per acre for dry yield. The average WUE of millet breeds in the second scenario ranges from 32.03 to 64.82 kg/ha·mm for fresh biomass and 10.8 to 24.6 kg/ha·mm for dry biomass.
In contrast to the first scenario, millet breeds can show good forage production (i.e., animal feed) as evidenced by insignificant variation in fresh and dry bio-mass and WUE, particularly for the IP19586, IP22269, IP19612, IP7704, and HHVBC tall B6 millet breeds. However, there are differences in the ratio of dry to fresh yield, which ranges from 0.31 at the MC94C2 breed to 0.53 at the HVBC tall B6 millet breed, despite the fact that the season varies greatly (for example, dry conditions). Figures 9 and 10 show that the millet breeds IP19586, IP22269, and IP19612 are significant good producers that can be used for forage production in Jordan's arid environment. Figure 9. Total fresh and dry biomass variabilities regarding millet breeds. Figure 9. Total fresh and dry biomass variabilities regarding millet breeds. Levels not connected by the same letter within the same column are significantly different at 95% confidence level. Figure 10. WUEs of total fresh and dry biomass regarding millet breeds.

Comparison of Seed Parameter among Millet Breeds
Millet breeds differ in terms of seed production and arid climate tolerance (Table 8). Under the Jordanian climate conditions and regardless of seasonal variations, the HHVBC tall B6 millet breed has the highest seed yield of 0.12 ton/ha and WUE of 0.39 07 kg/ha·mm. As a result, it is able to tolerate drought conditions and reproduce properly. In contrast to the preceding, the IP22269 breed, with the lowest seed yield of 0.02 ton/ha and seed WUE of 0.05 kg/ha·mm, fails to tolerate Jordan's drought conditions and propagate properly. Levels not connected by the same letter within the same column are significantly different at 95% confidence level.
In this study, the millet breeds' propagation capability can be categorized into three classes based on HI; (1) the high capability class (HI of above 0.04) as HHVBC tall B6 breed, (2) the moderate capability class (0.02 < HI < 0.04) as IP6110, MC94C2, Sudan pop

Comparison of Seed Parameter among Millet Breeds
Millet breeds differ in terms of seed production and arid climate tolerance (Table 8). Under the Jordanian climate conditions and regardless of seasonal variations, the HHVBC tall B6 millet breed has the highest seed yield of 0.12 ton/ha and WUE of 0.39 07 kg/ha·mm. As a result, it is able to tolerate drought conditions and reproduce properly. In contrast to the preceding, the IP22269 breed, with the lowest seed yield of 0.02 ton/ha and seed WUE of 0.05 kg/ha·mm, fails to tolerate Jordan's drought conditions and propagate properly. In this study, the millet breeds' propagation capability can be categorized into three classes based on HI; (1) the high capability class (HI of above 0.04) as HHVBC tall B6 breed, (2) the moderate capability class (0.02 < HI < 0.04) as IP6110, MC94C2, Sudan pop I, IP22269, IP19612, and P. millet icms7709 breeds, and (3) the poor capability class (HI < 0.02) as IP7704, Sudan pop III, IP13150, and IP19586 breeds. Kumari explained that the high WUE for biomass yield purposes obtained under drought stress is mainly due to the improved HI [68].

Correlation of Drought Indicators with Millet Growth and Production Parameters
The correlation matrix between the millet growth parameters and production with drought indices shows that the SPI and PDI precipitation drought indices are the most significant indices that can represent the smallest variations between millet breeds, followed by TDI. VDI, on the other hand, showed little correlation with millet growth parameters and production. This implies that the CDI is preferred for combining only the PDI and TDI indices (Table 9). According to the correlation matrix, there is a strong relationship between the amount of water used and the precipitation and drought indices. The drier the environmental conditions, the more water required to compensate for the soil-water-plant deficit. However, most of the applied water under drought conditions does not necessary uptake by the plant rather than being evaporated from the soil system.
As a result, the water use term is not recommended for correlating plant growth and yield factors because this flawed expression may lead to the crop being less water efficient when it could be the most promising climate-resilient crop. As a result, the water use term is not recommended for correlating plant growth and yield factors because this flaw expression may lead to the crop being less water efficient when it could be the most promising climate resilient crop. Furthermore, the SPI and PDI indices show a moderate (0.3 ≥ r ≥ 0.6) correlation with all other millet growth parameters and yield, indicating that millet is susceptible to drought conditions. According to the matrix, pearl millet tolerates drought conditions significantly by increasing root elongation and density to increase soilwater uptake from deeper soil layers. Furthermore, its phenological characteristics improve its drought tolerance by reducing leaf area and adjusting osmotic potential. A similar explanation was found by Cattivelli et al. pearl millets holding morpho-physiological traits can tolerate drought through stomatal conductance, photosynthetic capacity, the timing of the phenological phases, stem reserve mobilization in drought stress, reduced leaf area, rooting depth and density, cuticular resistance and surface roughness, osmotic adjustment, membrane composition, antioxidative defense, and accumulation of stressrelated proteins [69].
The TDI indicator was able to moderately correlate with water use, plant diameter and height, root depth, and fruit length and diameter, which could be attributed to crop physiology influences during heatwaves or extreme temperatures. The presence of higher-than-average temperatures simulates stress factors in millet, increasing its ability to withstand stress. According to others, Pearl millet is a hardy, climate-smart grain crop that is ideal for environments prone to drought and heat stress [43,48]. According to Dai, the millet phenological character compensates for its potential in case of any stress [4]. Depending on the critical stage occurrence during the plant growth stages, the millet breeds impact differently. According to Aparna et al., daytime maximum temperatures above 42 • C, as well as the associated vapor pressure deficit (VPD) during flowering time, directly reduce seed setting in pearl millet. [62]. Increased air temperature, on the other hand, may result in decreased reproduction due to heat damage to the pollen grain and increased sterility [70][71][72].
The precipitation and heat-stress drought indices have a moderate correlation with the WUEs, indicating that millets can withstand climatic changes to varying degrees depending on millet variety. The WUE is more influenced by edaphic than climatic conditions [64]. Furthermore, there are several categories of millets: (1) Agronomy-related traits, that is known as yield and yield components; (2) Morphology-related traits that are drought tolerance through changing its morphological and biochemical properties as flag leaf and leaf tensile strength; (3) Physiology-related traits that tolerate drought through osmotic adjustment; and (4) Biochemical-related traits that produce antioxidants to cope with unfavorable reactive oxygen species during stress conditions [49]. Finally, according to Sultan et al., photoperiod-sensitive pearl millet traditional cultivars appear to be more resistant to future climate conditions than improved cultivars with high genetic yield potential [73]. As a result, it is recommended that pearl millet cultivation in Jordan be strategically managed and adapted by recognizing production options in light of anticipated climate change scenarios.

Conclusions
Some pearl millet breeds demonstrated moderate to high resilience to extreme environmental conditions, particularly insufficient moisture and heat stress, and thus can withstand climate change. Although there are differences in millets breeds' responses to drought and heat stresses, as evidenced by differences in growth and plant production indicators, they have proven to adapt to the arid Jordanian environment and thus sustain their productivity in drought-prone conditions. The IP19586, IP22269, and IP19612 millet breeds are preferred for vegetative production because they provide good forage support. Their total yields were approximately 19 ton/ha fresh yield and approximately 7 ton/ha dry yield, which is twice as much as other millet breeds. Furthermore, the breeds mentioned above were bred to withstand drought conditions in Jordan, where they provided the highest WUEs (above 60 kg/ha·mm) for fresh biomass and (above 20 kg/ha·mm) for dry biomass. In terms of propagation characteristics, the HHVBC tall B6 breed exhibits good seed propagation with an average production of 0.12 ton/ha and a WUE of 0.39 07 kg/ha·mm, whereas the IP22269 breed demonstrated significantly lower seed yield, indicating a failure to propagate properly in arid-drought conditions. The correlation of millet growth and drought indices (VDI) revealed a very weak relationship (r < 0.3) and thus cannot be used for drought monitoring in Jordan. Temperature and precipitation drought indices both had strong correlations with millet growth factors (r > 0.6) and a moderate correlation with yield factor (0.3 < r < 0.6). As a result, the findings of this study suggest that a combined precipitation-temperature drought index be used to investigate, monitor, and assess millet growth parameters and production under climatic stress. Furthermore, based on these findings, the authors rec-ommend that other drought indices be investigated in order to develop a pearl millet monitoring plan.
Despite the fact that this study demonstrated that pearl millet breeds were able to tolerate drought conditions significantly due to their phenological characteristics, and thus set as a climate smart grain crop in Jordan, more research and replications are needed to investigate spatial and temporal variability under various Jordan climate zones.  Data Availability Statement: Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.