The Relationship between Changes in Hydro-Climate Factors and Maize Crop Production in the Equatorial African Region from 1980 to 2021

: Agricultural production across the African continent is subjected to various effects of climate variability. One of the main staple foods in Sub-Saharan Africa is maize. However, limited scientific research has recently focused on understanding the possible effects of hydro-climatic variability on maize production. The aim of the present work was to contribute to policy and climate adaptation, thus reducing the vulnerability of maize production to climate change over Equatorial Africa. This study firstly examined long-term trends of precipitation (PRE), soil moisture (SM), actual evapotranspiration (E), and potential evapotranspiration (Ep), as well as surface air temperatures, including the minimum (TMIN) and maximum (TMAX). Secondly, the relationship between maize production and these climate variables was quantified for 18 Equatorial African countries (EQCs) over 1980 − 2021. To assess the linear trends, Mann–Kendall and Sen’s slope tests were used to quantify the magnitude of the hydro-climatic variable trends at the 5% significance level, and Pearson’s correlation coefficient was used to evaluate the relation of these climate parameters with the maize production. The annual mean PRE declined at 0.03 mm day − 1 10a − 1 . Other climate variables increased at different rates: SM at 0.02 mmday − 1 10a − 1 , E at 0.03 mm day − 1 10a − 1 , Ep at 0.02 mm day − 1 10a − 1 , TMIN and TMAX at 0.01 ◦ C day − 1 10a − 1 . A regional analysis revealed heterogeneous significant wet–dry and warm–cool trends over the EQCs. While, spatially, dry and warm climates were observed in the central to eastern areas, wet and warm conditions dominated the western regions. Generally, the correlations of maize production with the E, Ep, TMAX, and TMIN were strong (r > 0.7) and positive, while moderate (r > 0.45) correlations of maize production with PRE and SM were obvious. These country-wide analyses highlight the significance of climate change policies and offer a scientific basis for designing tailored adaptation strategies in rainfed agricultural regions.


Introduction
The rising global average temperature and changing precipitation regimes are expected to alter moisture conditions in various global land regions [1].However, there are currently numerous differences between such changes in different global regions-with either positive, negative, or stationary trends in the moisture conditions [1,2].Changes in moisture conditions widely affect the agricultural sector [3,4].Simultaneously, increasing food supplies to meet the demands of the growing population has been of great concern, as climate change poses significant threats to sustainable food production [1,2,5].Moreover, countries in sub-Saharan Africa (SSA) are vulnerable to climate change threats as 70% of their economy relies heavily on agriculture [6][7][8].Crop production in many SSA countries is mainly rainfed and carried out by small-holder farmers who lack adequate resources [6,9,10].In SSA, cereal crops are predominantly grown on a subsistence basis to meet nutritional needs [11].The maize crop represents the majority of the dietary or nutritional components of food in SSA, highlighting its high use [6][7][8][11][12][13].According to the Food and Agriculture Organization [FAO], maize (Zea mays L.) is among the most widely cultivated cereal crops and a major dietary component in food [12,14].However, in a warming climate, maize yields are affected by drought or flood events [4,[15][16][17], which has been documented across different global land areas [14,15,18] using crop models [19], regional climate scenarios based on simulation data [20,21], and statistical methods based on historical data [22].
From a statistical perspective, the relationships between crop yield or production and climate variables (mainly precipitation and temperature) have been documented.In SSA countries, due to data availability, precipitation deficits and temperature anomalies are often used to quantify the roles of climate variables in explaining changes in crop yields [22][23][24].Furthermore, many studies have focused on single national-or countrylevel analyses [22,25], but there are little transboundary or regional studies [7].Maize production is influenced by moisture conditions (as a limiting factor for plant growth under certain temperature conditions) [4,15].However, as gridded datasets are available from various sources (e.g., gauges, reanalyses, and satellites), we can assess how other climate variables such as atmospheric water demand (i.e., evapotranspiration) [26] and soil moisture availability to plants (soil moisture) [27] influence maize growth in certain regions.Assessing the trends and seasonality of several climate variables has gained attention [28,29].For example, previous studies have deduced the trends and seasonality of several climate parameters in transboundary research on SSA [29], West Africa [28], and East Africa [30], but there are still differences in the reported trend dynamics due to the different spans of study areas and data periods used [31].The impact of climate change on maize production or yields has also not been comprehensively studied across SSA [7,21].The few studies on the impacts of climate change on crop production or yield (i.e., focusing on precipitation and temperature variables) were limited to single-country analyses [22,[32][33][34].Since climate variables have different dynamics, their impact is still uncertain considering that, to date, not many national statistical analyses have been conducted on climate dynamics and crop production in the Equatorial Africa region, where food insecurity still persists [8].There are still gaps in the research on the impact of climate variability on maize production in the Equatorial Africa region.This study not only includes a trend analysis, as analyzed in previous studies, but we also analyzed the dynamics of six climate variables in 18 countries, which was previously unexplored in past studies.In addition, this study provides analyses on the responses to climate seasonality and trends in these six climate variables.It was necessary to conduct an in-depth investigation into the impacts of climate trends and dynamics on maize production for multiple countries in the Equatorial Africa region.To the best of our knowledge, this study presents the first multiple-country analysis of these six (6) climate variables based on the annual data and growing seasons of maize production in 18 countries in this region.Thus, the aim of this paper was to identify significant trends in multiple climate variables and their association with maize production in Equatorial African countries.The study objectives were as follows: (1) to investigate the annual and inter-annual variability and trends of precipitation, soil moisture, actual and potential evapotranspiration, and air surface temperatures (minimum and maximum) from 1980 to 2021; (2) to evaluate maize production and climate trends in 18 selected Equatorial African countries (EQCs); and (3) to determine the relationships between these climatic variables and maize production.
The remainder of this paper is structured as follows: Section 2 introduces the study area, data sources, and methods used in this study.Section 3 presents the study results.Section 4 discusses the results of the study.Section 5 concludes the research work.

Study Area
The Equatorial Africa region is located in the tropics and this study focused on regions at a longitude of    N-20 • N. The study area (Figure 1a) is divided into 8 major climatic zones, namely semi-arid, arid, and humid zones (see Figure 1b).The study area shows that lowlands are distributed in the western regions and mountains are located in the eastern region (Figure 1c).The region is dominated by oceanic (Atlantic and Indian Oceans) and remote (teleconnections) influences [35].The terrain of the study area ranges from −380 to 5885 m above sea level, and the highest altitude of 5800 m above sea level is in the Ethiopian Highlands (Figure 1c).The dominant vegetation types in the region are forests, grasses, shrubs, and crops (Figure 1d).The cereal crops grown in this tropical African region are mainly interspersed in dominant vegetations such as forests, grasses, and shrubs [36] (Figure 1d).All of these complex regional features, together with other climate features, define the region's distribution of rainfall and temperatures, which influences agricultural crop production [28,37].
Atmosphere 2024, 15, x FOR PEER REVIEW 3 of 24 the annual and inter-annual variability and trends of precipitation, soil moisture, actual and potential evapotranspiration, and air surface temperatures (minimum and maximum) from 1980 to 2021; (2) to evaluate maize production and climate trends in 18 selected Equatorial African countries (EQCs); and (3) to determine the relationships between these climatic variables and maize production.The remainder of this paper is structured as follows: Section 2 introduces the study area, data sources, and methods used in this study.Section 3 presents the study results.Section 4 discusses the results of the study.Section 5 concludes the research work.

Study Area
The Equatorial Africa region is located in the tropics and this study focused on regions at a longitude of 18 °W-55 °E and latitude of 2 °N-20 °N.The study area (Figure 1a) is divided into 8 major climatic zones, namely semi-arid, arid, and humid zones (see Figure 1b).The study area shows that lowlands are distributed in the western regions and mountains are located in the eastern region (Figure 1c).The region is dominated by oceanic (Atlantic and Indian Oceans) and remote (teleconnections) influences [35].The terrain of the study area ranges from -380 to 5885 m above sea level, and the highest altitude of 5800 m above sea level is in the Ethiopian Highlands (Figure 1c).The dominant vegetation types in the region are forests, grasses, shrubs, and crops (Figure 1d).The cereal crops grown in this tropical African region are mainly interspersed in dominant vegetations such as forests, grasses, and shrubs [36] (Figure 1d).All of these complex regional features, together with other climate features, define the region's distribution of rainfall and temperatures, which influences agricultural crop production [28,37].[38], (c) elevation [39], and (d) land use/land cover (LULC) [40].

Climate Datasets
We used the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center (CPC) Global Unified Gauge-Based Analysis of Daily Precipitation and Temperature for this study [41].The CPC Unified_v1.0/TD dataset is a high-resolution  [38], (c) elevation [39], and (d) land use/land cover (LULC) [40].

Data Sources 2.2.1. Climate Datasets
We used the National Oceanic and Atmospheric Administration (NOAA) Climate Prediction Center (CPC) Global Unified Gauge-Based Analysis of Daily Precipitation and Temperature for this study [41].The CPC Unified_v1.0/TD dataset is a high-resolution gauge-based gridded dataset developed by the NOAA CPC with a spatial resolution of 0.5 • by 0.5 • that includes monthly timestep data on climate variables [42].
This dataset is publicly available and was downloaded from the NOAA portal [43].We also obtained the following satellite-based datasets.The Global Land Evaporation Amsterdam Model (GLEAM) provides satellite-based data available from 1980 to 2021, with a high resolution of 0.25 • × 0.25 • .Daily scales for the surface soil moisture and actual and potential evapotranspiration were obtained from datasets on the GLEAM website [44].For detailed information on the GLEAM datasets, readers are directed to the following references: [45][46][47].

Maize Crop-Based Data
Country-level data on maize statistics were obtained from the US Department of Agriculture (USDA) website [48].We obtained data on the maize production and harvested area for countries in the EQA region (Table 1), with time-series data from 1950 to the present.The period of data used for statistical analysis was from 1980 to 2021.The climate data analyses were performed using Climate Data Operators (CDO) version 1.9cl [49].In addition, daily data were converted into annual values to match the maize crop data.

Maize Yield Computation
We computed the maize yields as an average ratio of production (1000 MT) per harvested area (1000 MHa).A limitation of this approach to computing yield is that it does not reflect the years in which maize production losses occurred.We computed the correlation at the country level and across the EQA region.

Index Anomalies
All datasets were standardized using the standardized anomalies formula (Equation (1)) based on the 1983-2021 climatology: where α std denotes the standardized α, µ is the average, and σ is the standard deviation of the variable α.

Linear Trends
We computed the linear trends using the Mann-Kendall trend test (Kendall's coefficient, Z) to determine the significance of the trends and the Sen's slope test to compute the magnitude of the trends [50,51].These non-parametric methods assume datasets to have independent and random distribution.An advantage of using non-parametric estimators is that they can handle outliers and these methods are widely used in many areas [35,52].All statistically significant trends were calculated using Student's t-test at a 95% confidence level.For more details on the original M-K and Sen's slope tests, refer to Mann [50] and Sen [51].

Correlation Analyses
We computed the correlation between country-wide maize production and climate anomalies from 1980 to 2021 based on Equation (3): where r denotes the correlation coefficient, x i denotes the climate variable, i denotes the time, and n denotes the sample size.

Results and Discussion
3.1.Annual Climatology and Trend Analysis 3.1.1.Climatology Figure 2 displays the spatial distributions of the average annual values of E, Ep, PRE, SM, TMAX, and TMIN during the period 1980-2021.E and Ep ranged from 0.8 to 5.6 mm and from 1.6 to 5.6 mm, respectively.The annual mean daily PRE was 4.5 mm and ranged from 1 to 8 mm (Figure 2c).The PRE showed a higher rainfall distribution.The annual mean daily surface SM was 0.25 m 3 m −3 and varied between 0.06 and 0.45 m 3 m −3 (Figure 2d).The annual mean daily TMAX showed a mean of 28.8 and ranged from 20 to 37.5 • C (Figure 2e).The TMIN showed an annual daily mean of 17.5 • C and ranged from 10 to 25 • C (Figure 2f).Generally, a gradient was observed with high values in humid regions (i.e., 2 • N-8 • N and 18 • W-18 • E) followed by semi-arid regions (i.e., 8 • N-14 • N) and low values in arid regions (14-20 • N) for E, Ep, PRE, and SM.On other hand, however, the temperature values TMAX and TMIN showed a latitudinal gradient: arid regions > semi-arid regions > humid regions (Figure 2e,d, respectively).

Distribution of Annual Linear Climate Trends
Figure 4 shows the spatial patterns of the annual linear trends of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.Overall, we observed a spatial heterogeneity of positive and negative trends across the study area.In arid regions, Ep (Figure 4b) displayed pronounced negative trends and E (Figure 4a) showed mixed trends.In semi-arid regions, the pattern of significant spatial trends for E and Ep were opposite (Figure 4a,c).For example, areas in the western Sahel (8 The spatial patterns of PRE and SM importantly showed spatial heterogeneity, which could be due to a possible interaction between these two variables.For example, the semi-arid regions showed that PRE (Figure 4c) displayed no changing trends, while the SM (Figure 4d) displayed mixed trends (positive trends in the western and eastern Sahel, and negative trends in the central Sahel).Within the humid regions, negative trends were observed in small areas encompassing the Guinean coast and in the northern Congo Basin area (Figure 4c,d).The spatial patterns of linear trends in TMAX (Figure 4e) and TMIN (Figure 4f) showed generally similar results in the western EQCs (18 • W-10 • E).Meanwhile, opposite trends were evidenced in the region located between 10 • and 52 • E.
Figure 5 displays the annual temporal trends of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.Generally, except for PRE, which exhibited a downward trend, the remaining variables displayed upward trends, albeit with differences in the magnitude of the trend values.E and Ep displayed upward trends of 0.03 mm (day −1 10a −1 ) and 0.02 mm (day −1 10a −1 ), respectively (Figure 5a,b).TMAX and TMIN showed a statistically significant increasing trend at 0.01 mm (day −1 10a −1 ) (Figure 5c,d).In contrast, PRE exhibited a downward trend of 0.03 mm (day −1 10a −1 ) (Figure 5c), while SM showed no change in trend (Figure 5d).Except for PRE and SM, all of the variables were statistically significant (p < 0.05).

Distribution of Annual Linear Climate Trends
Figure 4 shows the spatial patterns of the annual linear trends of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.Overall, we observed a spatial heterogeneity of positive and negative trends across the study area.In arid regions, Ep (Figure 4b) displayed pronounced negative trends and E (Figure 4a) showed mixed trends.In semiarid regions, the pattern of significant spatial trends for E and Ep were opposite (Figure negative trends were observed in small areas encompassing the Guinean coast and in the northern Congo Basin area (Figure 4c,d).The spatial patterns of linear trends in TMAX (Figure 4e) and TMIN (Figure 4f) showed generally similar results in the western EQCs (18° W-10° E).Meanwhile, opposite trends were evidenced in the region located between 10° and 52° E. Figure 5 displays the annual temporal trends of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.Generally, except for PRE, which exhibited a downward trend, the remaining variables displayed upward trends, albeit with differences in the magnitude of the trend values.E and Ep displayed upward trends of 0.03 mm (day −1 10a −1 ) and 0.02 mm (day −1 10a −1 ), respectively (Figure 5a,b).TMAX and TMIN showed a statistically significant increasing trend at 0.01 mm (day −1 10a −1 ) (Figure 5c,d).In contrast, PRE exhibited a downward trend of 0.03 mm (day −1 10a −1 ) (Figure 5c), while SM showed no change in trend (Figure 5d).Except for PRE and SM, all of the variables were statistically significant (p < 0.05).The spatial distribution of TMAX (Figure 6e) and TMIN (Figure 6f) showed similar patterns of warm-and-cool trends, except for the divergent trends observed in Ghana, Mali, Nigeria, Chad, Sudan, South Sudan, and Somalia.

Country-Level Maize Production and Yield Estimates
Figure 7 presents a time-series analysis of the maize production, yield, and harvested land area in selected countries located in the study area.High year-to-year variability is displayed in the maize production (Figure 7, line graph) and harvested area (Figure 7, bar graph).
The maize yield behavior, depicted in Figure 7, was highly variable between the countries and years, with most countries displaying relatively stable yields (Benin, Burkina Faso, Cameroon, Cote d'Ivoire, Chad, Ghana, Guinea, Mali, Ethiopia, Nigeria, Senegal, Somalia, and Togo) throughout time.Few countries displayed strong fluctuations (e.g., Guinea-Bissau, the CAR, Mauritania, Eritrea, Somalia, and Gambia).The con-

Country-Level Maize Production and Yield Estimates
Figure 7 presents a time-series analysis of the maize production, yield, and harvested land area in selected countries located in the study area.High year-to-year variability is displayed in the maize production (Figure 7, line graph) and harvested area (Figure 7, bar graph).

Descriptive Statistics of Annual Maize and Climate
To reveal the temporal differences in maize production and hydro-climate trends, the mean (standard deviation) values were computed for the selected countries from 1980 to 2021 (Table 2).Overall, Nigeria presented the largest average maize production of 6221.7 ± 3325.8MT.The annual daily climate variables for Nigeria were SM (0.27 ± 0.09 m −3 m −3 ), PRE (2.92 ± 2.9 mm), TMIN (22.8 °C), TMAX (33.19 ± 2.9 °C), E (2.06 ± 0.9 mm), and Ep (3.05 mm).The overall statistics across the EQCs ranged between 0.45 × 10 4 and 3.8 × 10 4 The maize yield behavior, depicted in Figure 7, was highly variable between the countries and years, with most countries displaying relatively stable yields (Benin, Burkina Faso, Cameroon, Cote d'Ivoire, Chad, Ghana, Guinea, Mali, Ethiopia, Nigeria, Senegal, Somalia, and Togo) throughout time.Few countries displayed strong fluctuations (e.g., Guinea-Bissau, the CAR, Mauritania, Eritrea, Somalia, and Gambia).The consistently higher maize production across the EQCs corroborates other studies that support maize as one of the vital food crops in this region [20,25] and other global land areas [53][54][55].

Country-Wide Trend Analysis of Climate Variables
Table 3 presents the trends in the hydro-climate variables, expressed according to the M-K test (Z) and Theil-Sen slope (β) values on an annual scale and a growing-season scale.Here, we present Z values that were statistically significant at α = 0.05.Most national annual climate parameters remained generally stable (trends without statistical significance).Our results focus on countries that experienced statistically significant positive/negative trends.
During the growing seasons, PRE increased significantly in Burkina Faso, Cote d'Ivoire, Ghana, Guinea, Guinea-Bissau, Mali, and Mauritania (p < 0.05).SM increased significantly in Burkina Faso, Eritrea, Gambia, Guinea-Bissau, Mali, Mauritania, and Senegal, while displaying negative trends in Cameroon, the CAR, and Ghana.E and Ep showed increasing trends in Benin, Burkina Faso, Cameroon, the CAR, Cote d'Ivoire, Ethiopia, Ghana, Guinea, Guinea-Bissau, Mali, Mauritania, Nigeria, Senegal, Somalia, and Togo.However, some countries showed increasing trends for E (Gambia) and Ep (Chad).TMAX and TMIN showed similar increasing trends in Benin, Burkina Faso, Cameroon, Cote d'Ivoire, Ghana, Guinea, Guinea-Bissau, and Mauritania.However, increasing trends were found for TMAX in Chad and Mali, and for TMIN in Gambia, Nigeria, Senegal, Somalia, and Togo.All trends were statistically significant (p < 0.05).8 shows the heat-map plots of the correlation coefficient (r) computed between the country-level maize production and climate anomalies.The results showed that the highest correlation was recorded for maize and TMAX (0.77) in Cote d'Ivoire, and the lowest was recorded for maize and PRE (−0.01) in Guinea-Bissau.SM and maize showed mixed correlation results, with a positive correlation in Burkina Faso (0.52), Gambia (0.67), Mali (0.45), Mauritania (0.59), and Senegal (0.43) and a negative correlation in Cameroon (−0.60) and the CAR (−0.64).The remaining correlation coefficient (r) values were insignificant at the 95% confidence level.On the other hand, PRE and maize displayed a moderately positive correlation coefficient (r) for Burkina Faso (0.47), Chad (0.43), Mali (0.58), and Senegal (0.48).
the country-level maize production and climate anomalies.The results showed that the highest correlation was recorded for maize and TMAX (0.77) in Cote d'Ivoire, and the lowest was recorded for maize and PRE (−0.01) in Guinea-Bissau.SM and maize showed mixed correlation results, with a positive correlation in Burkina Faso (0.52), Gambia (0.67) Mali (0.45), Mauritania (0.59), and Senegal (0.43) and a negative correlation in Cameroon (−0.60) and the CAR (−0.64).The remaining correlation coefficient (r) values were insignif icant at the 95% confidence level.On the other hand, PRE and maize displayed a moder ately positive correlation coefficient (r) for Burkina Faso (0.47), Chad (0.43), Mali (0.58) and Senegal (0.48).
Regarding E and Ep, our results showed a significant positive correlation with maize production in Benin, Burkina Faso, Cameroon, Cote d'Ivoire, Ghana, Guinea, Mali, Mau ritania, Nigeria, Senegal, and Togo.However, only a few countries showed a significan positive correlation with E (Gambia) or Ep (Ethiopia, Eritrea, Chad), and a negative one was observed for Somalia (Ep).Similarly, a significant positive correlation was found fo the temperature variables (TMIN and TMAX) in Benin, Burkina Faso, Cote d'Ivoire Ghana, Guinea, Mali, and Togo.TMIN also showed a positive trend for Cameroon, the CAR, Eritrea, Gambia, Mauritania, and Nigeria.All trends were statistically significant a the α = 0.05 level.Regarding E and Ep, our results showed a significant positive correlation with maize production in Benin, Burkina Faso, Cameroon, Cote d'Ivoire, Ghana, Guinea, Mali, Mauritania, Nigeria, Senegal, and Togo.However, only a few countries showed a significant positive correlation with E (Gambia) or Ep (Ethiopia, Eritrea, Chad), and a negative one was observed for Somalia (Ep).Similarly, a significant positive correlation was found for the temperature variables (TMIN and TMAX) in Benin, Burkina Faso, Cote d'Ivoire, Ghana, Guinea, Mali, and Togo.TMIN also showed a positive trend for Cameroon, the CAR, Eritrea, Gambia, Mauritania, and Nigeria.All trends were statistically significant at the α = 0.05 level.

Discussion
Climate change and natural variability threatens agricultural crop production, especially in water-and energy-limited regions, where food crops are rainfed and harvested through subsistence farming by small-holder farmers.
We investigated the climatology and linear trends of six hydro-climate variables and their impacts on maize production in selected countries located in the Equatorial Africa region.The spatial pattern analysis captured and revealed the observed climatology and observed trends of hydro-climate variables across the study area to be consistent with those reported in a past study [28].Some locations were captured with more pronounced rainfall patterns and less pronounced temperature patterns, or vice versa.The EQCs showed annual cycles of precipitation (with amounts ranging from 0 to 2.0 mm day −1 ) and temperatures (i.e., maximum air surface temperatures between 20 and 38 • C day −1 and minimum air surface temperatures in the range of 10 and 25 • C day −1 ) that are consistent with the ideal rainfall and temperatures for maize crop production [56].According to the FAO [9], the sowing seasons are in May-August and September to December; however, the intra-annual variability captured in this study did not show that pattern.This could be due to the area-averaged values we applied or the selected region not being similar to that used by the FAO.The annual monthly distributions of the hydro-climate variables suggest that the maize crop could have at least two cropping calendars [13,24].Maize is one of the most widely grown crops in EQA, which encompasses arid regions that transition through to rainforest regions [57].Spatially, the precipitation displayed significant positive changes mostly in the Sahel region, meaning that the Sahel experienced increases in precipitation.This confirms the earlier findings for this region.The southern part of the Sahel (Guinea Coast) experienced negative trends, thus culminating in an average negative trend of 0.03 mm/10a −1 , consistent with previous regional studies [30,58].Linear trends in the soil moisture show an augmentation in the region's precipitation over certain regions and a reduction in it in arid and humid regions, which is in line with past studies [59][60][61].The spatial pattern of evapotranspiration showed similar trends for actual and potential evapotranspiration.However, both types of evapotranspiration increased at rates of 0.03 mm/10a −1 and 0.02 mm10a −1 , respectively.In the context of a warming climate, increasing evapotranspiration may lead to an increase in surface water losses and increased water availability, as supported by previous studies [37,52,62].The spatial patterns of the minimum and maximum temperatures showed heterogeneity in their trends, but the overall values showed an identical trend rate of 0.01 • C 10a −1 .Our results are supported by previous regional studies [30,63,64] and studies on global climate patterns [1,5].TMAX and TMIN were reported with rates of 1.9 • C and 1.2 • C in East Africa by Gebrechorkos, Hülsmann and Bernhofer [30], and rates of 0.2 • C 10a −1 for the mean minimum temperature and 0.10 • C10a −1 for the mean maximum temperature were reported in the IPCC AR5 report [5].Overall, the values of trends and their significance depend on the region and time period under consideration [31].Although the results of the temporal pattern analysis (Figure 7) of annual maize production differed among the EQCs (Table 2), the mean (standard deviation) climatic variables suggest good ideal conditions for maize production, which largely explains the popularity of maize production in this region, as maize is a staple food in the region [20,25] and throughout the world [53][54][55].The quantity of maize produced by each country maybe explains the largely high inter-annual variations in the harvested area.Nigeria and Mauritania produced the highest and lowest amount of maize over the four decades (i.e., 1980-2021).Maize cultivation in many of these countries is rainfed and the yields of maize are maybe impacted by the harvested area.The harvested area mostly influences decisions on the space (i.e., where) and time (i.e., when) to grow crops [16].This in turn is partly due to a combination of different hydro-climate variables [55] or other limiting factors (e.g., the number of crops grown, the maize survival rate, the sowing date [16,17], general farming inputs (e.g., investments [65,66]) or profitability [6], technology [11,67], improved seed varieties [23], or fertilizer applications [68]).At the national level, the spatial pattern analysis showed significant (p < 0.05) wetting and warming over Burkina Faso, Mali, Mauritania, Senegal, Gambia, Guinea-Bissau, and Togo (Table 3).Also, the country-level analysis showed that the drying and warming over Benin, Ghana, Cameroon, and the CAR was significant at α = 0.05 (Table 3).Wet and warm conditions are essential in promoting agricultural development due to crops' water requirements [57,69].Our correlation analysis showed that the climate drivers' influence on maize production differed from country and country.Focusing on countries with significant values, Burkina Faso, Chad, Mali, and Senegal showed their maize production and precipitation to be positively correlated.Soil moisture was negatively correlated in Cameroon and the CAR, and positively correlated in Burkina Faso, Gambia, Mali, and Mauritania.Our analyses suggested that maize production is highly sensitive to different water conditions, which may influence the sowing dates [16,17].The temperature's correlation with maize showed a strongly positive correlation in nearly all of the countries (significant at p < 0.05), which is not surprising.These findings are supported by studies reporting that temperatures (i.e., minimum and maximum) control plant metabolism [32,70], where maximum temperatures enhance plant photosynthesis and minimum temperatures improve plant respiration and nutrient consumption to subsequently increase crop production (or yields) [71].However, countries that exhibit weakly negative and non-significant correlations between maize and precipitation and soil moisture may be prone to flash droughts.This is consistent with the regions that experienced high inter-annual variations in precipitation (Figures 4-6), with higher temperatures (Figure 4e,f) and water loss through ET (Figure 4a,b) further exacerbating the soil moisture conditions (Figure 4d) when there was a lack of precipitation (Figures 4c and 5c).Our findings confirmed that evapotranspiration, minimum and maximum temperatures, and soil moisture influence maize production across the EQCs.Overall, our correlation analysis showed that climate drivers impact maize production differently across the countries considered, thus suggesting the need for tailor-made responsive measures, consistent with the literature [22,72,73].The correlation analysis indicated that maize production in many of the studied regions is related to the trends of these climate drivers [22].Additional analyses showed that large-scale climatic phenomena such as the El Niño Southern Oscillation (ENSO) also influence maize production (Figures S1-S5).
Previous research on climate drivers' impacts on crop production has focused largely on mean temperatures and precipitation in different parts of SSA [7,22].However, there is limited research examining other factors such as evapotranspiration, temperatures (i.e., minimum and maximum), and soil moisture in the Equatorial Africa region.This study has analyzed the influence of climate variability on maize production.Also, we have explored the effects of other factors such as evapotranspiration, temperatures (i.e., minimum and maximum), and soil moisture on maize production.Future studies should quantify the relationship between climate change and crop production and examine the causality of the trends in climatic factors and global warming affecting maize production.
In summary, climate drivers mostly influence farmers' decisions on where (i.e., space) and when (i.e., time) to grow crops in rainfed regions.This study, consistent with previous studies, demonstrated that climate change influences maize production [20,74].However, we advise readers to treat these findings with caution as there are other unexplained limiting factors.Maize production is known to be influenced by factors such as the maize variety [23,25,68,75], technological investments [25,66,67], etc.Our study contributes to the ongoing discussion on the state of Africa's warming climate and its impact on different sectors of the economy.This study provides in-depth analyses of country-wide spatiotemporal changes in climate patterns in relation to maize production.These country-wide analyses are vital for effectively managing agricultural systems, and provide a scientific basis for adapting tailored policies and programs in rainfed agricultural regions.

Conclusions
The existing research on the effects of hydro-climate parameters on maize production over Equatorial African countries, especially considering country-wide analyses, is very limited.Hence, this study analyzed changes in hydro-climate variables (i.e., actual (E) and
Figure2displays the spatial distributions of the average annual values of E, Ep, PRE, SM, TMAX, and TMIN during the period 1980-2021.E and Ep ranged from 0.8 to 5.6 mm and from 1.6 to 5.6 mm, respectively.The annual mean daily PRE was 4.5 mm and ranged from 1 to 8 mm (Figure2c).The PRE showed a higher rainfall distribution.The annual mean daily surface SM was 0.25 m 3 m −3 and varied between 0.06 and 0.45 m 3 m −3 (Figure2d).The annual mean daily TMAX showed a mean of 28.8 and ranged from 20 to 37.5 • C (Figure2e).The TMIN showed an annual daily mean of 17.5 • C and ranged from 10 to 25 • C (Figure2f).Generally, a gradient was observed with high values in humid regions (i.e., 2 • N-8 • N and 18 • W-18 • E) followed by semi-arid regions (i.e., 8 • N-14 • N) and low values in arid regions (14-20 • N) for E, Ep, PRE, and SM.On other hand, however, the temperature values TMAX and TMIN showed a latitudinal gradient: arid regions > semi-arid regions > humid regions (Figure2e,d, respectively).Figure3shows box plots of the annual cycles of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.E displays a multi-year mean of 2.0 (range: 1-3.0 mm/day) (Figure3a).Ep shows a mean value of 3.25 (range: 1-3.0 mm day −1 ) (Figure3b).PRE records a multi-year mean of 1.0 (0-2.0 mm day −1 ) (Figure3c).SM displays a multi-year mean of 0.175 (range: 0.125-0.225m 3 m −3 day −1 ).TMAX and TMIN show mean daily values of 32 (range: 28-36 • C) and 24 (range: 22-26 • C day −1 ), respectively.The temporal pattern shows that E and Ep follow similar patterns, displaying high values in March and low values in July.PRE shows high values in November and low values in January.SM shows high values in November and low values in May.TMAX peaks from July to January, while TMIN peaks from August to January.

Figure 3
Figure3shows box plots of the annual cycles of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.E displays a multi-year mean of 2.0 (range: 1-3.0 mm/day) (Figure3a).Ep shows a mean value of 3.25 (range: 1-3.0 mm day −1 ) (Figure3b).PRE records a multi-year mean of 1.0 (0-2.0 mm day −1 ) (Figure3c).SM displays a multi-year mean of 0.175 (range: 0.125-0.225m 3 m −3 day −1 ).TMAX and TMIN show mean daily values of 32 (range: 28-36 °C) and 24 (range: 22-26 °C day −1 ), respectively.The temporal pattern shows that E and Ep follow similar patterns, displaying high values in March and low values in July.PRE shows high values in November and low values in January.SM shows high values in November and low values in May.TMAX peaks from July to January, while TMIN peaks from August to January.
Figure4shows the spatial patterns of the annual linear trends of E, Ep, PRE, SM, TMAX, and TMIN over the period 1980-2021.Overall, we observed a spatial heterogeneity of positive and negative trends across the study area.In arid regions, Ep (Figure4b) displayed pronounced negative trends and E (Figure4a) showed mixed trends.In semi-arid regions, the pattern of significant spatial trends for E and Ep were opposite (Figure4a,c).For example, areas in the western Sahel (8• -14 • N, 18 • W-15 • E) and central Sahel (8 • -14 • N, 15 • -35 • E) displayedpositive and negative trends for E and Ep, respectively.In humid regions, E and Ep showed similar spatial patterns with mixed trends along the Guinean coast (2 • -8 • N, 18 • W-8 • E) and the Horn of Africa (HOA) region (2 • -8 • N, 36-52 • E).We also observed opposite trends in E and Ep over the northern Congo Basin region (2 • -8 • N, 8 • -35 • E).The spatial patterns of PRE and SM importantly showed spatial heterogeneity, which could be due to a possible interaction between these two variables.For example, the semi-arid regions showed that PRE (Figure4c) displayed no changing trends, while the SM (Figure4d) displayed mixed trends (positive trends in the western and eastern Sahel, and negative trends in the central Sahel).Within the humid regions, negative trends were observed in small areas encompassing the Guinean coast and in the northern Congo Basin area (Figure4c,d).The spatial patterns of linear trends in TMAX (Figure4e) and TMIN (Figure4f) showed generally similar results in the western EQCs (18 • W-10 • E).Meanwhile, opposite trends were evidenced in the region located between 10 • and 52 • E.

24 Figure 3 .
Figure 3. Box plots of the monthly climatology of (a) actual evapotranspiration (E), (b) potential evapotranspiration (Ep), (c) precipitation (PRE), (d) surface soil moisture (SM), (e) maximum temperature (TMAX), and (f) minimum temperature (TMIN) across Equatorial Africa (EQA) from 1980 to 2021.Error bars represent one standard deviation of uncertainty, while the circles represent the outliers.The green dashed lines and orange lines indicate the median values and mean values, respectively.

Figure 3 .
Figure 3. Box plots of the monthly climatology of (a) actual evapotranspiration (E), (b) potential evapotranspiration (Ep), (c) precipitation (PRE), (d) surface soil moisture (SM), (e) maximum temperature (TMAX), and (f) minimum temperature (TMIN) across Equatorial Africa (EQA) from 1980 to 2021.Error bars represent one standard deviation of uncertainty, while the circles represent the outliers.The green dashed lines and orange lines indicate the median values and mean values, respectively.

Atmosphere 2024 , 24 Figure 6 .
Figure 6.Spatial patterns of the M-K trends of (a) actual evapotranspiration (E), (b) potential evapotranspiration (Ep), (c) precipitation (PRE), (d) surface soil moisture (SM), (e) maximum temperature (TMAX), and (f) minimum temperature (TMIN) across Equatorial Africa (EQA) from 1980 to 2021.The unit is mm d −1 .Analyzed based on the MK and Sen's slope tests.The red and blue colors denote decreasing and increasing trends, respectively.The dots denote the regions where the trends are statistically significant at a confidence level of 95%.

Figure 6 .
Figure 6.Spatial patterns of the M-K trends of (a) actual evapotranspiration (E), (b) potential evapotranspiration (Ep), (c) precipitation (PRE), (d) surface soil moisture (SM), (e) maximum temperature (TMAX), and (f) minimum temperature (TMIN) across Equatorial Africa (EQA) from 1980 to 2021.The unit is mm d −1 .Analyzed based on the MK and Sen's slope tests.The red and blue colors denote decreasing and increasing trends, respectively.The dots denote the regions where the trends are statistically significant at a confidence level of 95%.

Figure 7 .
Figure 7. Changes in maize production (MT, green color), harvested land area (HA, blue color), and yield (MT/ha, red color).The dotted and dashed lines in the plots denote the linear trends of the maize production and yields, respectively; the bar plots indicate the harvested area during 1980 to 2021.

Figure 7 .
Figure 7. Changes in maize production (MT), harvested land area (HA), and yield (MT/ha).The dotted and dashed lines in the plots denote the linear trends of the maize production and yields, respectively; the bar plots indicate the harvested area during 1980 to 2021.

Figure 8 .
Figure 8.Heat map of Pearson correlation coefficient (r) values between the annual maize produc tion and hydro-climate variables.* denote statistical significance at p < 0.05.

Figure 8 .
Figure 8.Heat map of Pearson correlation coefficient (r) values between the annual maize production and hydro-climate variables.* denote statistical significance at p < 0.05.

Figure A1 .
Figure A1.Changes in maize production (MT), harvested land area (ha), and yield (MT/ha) in Equatorial African countries (EQCs).The dotted and dashed lines in the plots denote the linear trends for maize production and yields, respectively; the bar plots indicate the harvested area during 1980 to 2021.

Table 1 .
List of selected countries in the Equatorial Africa region.

Table 2 .
Descriptive statistics of annual maize production and climate variables during 1980-2021.

Table 3 .
Country-level summary statistics according to the M-K test (Z) and Sen's slope estimator (β), annually and in the growing seasons, for the maize and climate variables over the period 1980-2021.* = significance level at 0.01 < α ≤ 0.05.The beta coefficient (β) indicates the changes per year or changes per growing season.