Spatio-Temporal Trends of Precipitation and Temperature Extremes across the North-East Region of C ô te d’Ivoire over the Period 1981–2020

: The northeast region of C ô te d’Ivoire, where agriculture is the main economic activity, is potentially vulnerable to extreme climatic conditions. This study aims to make a comprehensive spatio-temporal analysis of trends in extreme indices related to precipitation and temperature for the Zanzan region of C ô te d’Ivoire over the period of 1981–2020. The statistical signiﬁcance of the calculated trends was assessed using the non-parametric Mann–Kendall test, while Sen’s slope estimation was used to deﬁne the amount of change. For extreme precipitations, the results showed a decreasing trend in annual total precipitations estimated at 112.37 mm and in daily precipitations intensity indices. Furthermore, the consecutive dry days’ index showed an increasing trend estimated at 18.67 days. Unlike the trends in precipitation extremes, which showed statistically non-signiﬁcant trends, the trends in temperature extremes were mostly signiﬁcant over the entire study area. The cold spells indices all show decreasing trends, while the warm spells show increasing trends. Drawing inferences from the results, it becomes clear that the study area may be threatened by food insecurity and water scarcity. The results are aimed to support climate adaptation efforts and policy intervention in the region.


Introduction
The achievement of sustainable development goals (SDGs), especially those related to ensuring environmental sustainability, reduction of poverty and hunger, and access to clean water among others are under severe threat from climate change [1,2]. This is important, particularly for Sub-Saharan Africa (SSA), where agriculture is the primary source of income and livelihood for many people. With nearly 19% of the region's gross output produced in the agricultural sector, SSA is the most agriculturally dependent region in the world after South Asia [3]. However, there is a high dependence of crops on weather events, with 97% of agricultural land being rainfed [4], and only 1% of cultivated land being irrigated [5,6]. For several decades climate change has strongly affected this part of Africa [7]. This manifestation of climate change in SSA is perceived through the increase in temperature, the presence of storms and extreme precipitation [8], the increase in land salinity, the worsening of weather conditions, and tornadoes that result in the ruin of ecosystems [9]. In its third assessment report, the Intergovernmental Panel on Climate Change (IPCC) stated that the number of hydro-meteorological disasters has increased Guinean climate. The main climate has four seasons, two of which are rainy seasons from March to June and from September to October. These rainy seasons are interspersed with two dry seasons which run from November to February and from July to August. The annual rainfall varies between 800 and 1200 mm with an average temperature of about 28 • C [54]. The Zanzan region has a surface area of approximately 38,000 km 2 . Based on the 2014 national population and housing census, its population is estimated at 934,352 inhabitants [55]. The region's economy is mainly based on agriculture, livestock production, trade (food crops, handicrafts, etc.), and services. Cashew nuts are the main cash crop in the region with a contribution of 22% of national production [49]. The cashew industry constitutes 8.8% of the GDP of Côte d'Ivoire with a production of 700,000 tons [56]. In addition to this crop, cocoa and coffee are produced in small quantities in the southern part of the Zanzan region.

Study Area
Côte d'Ivoire is a West African country currently divided into 31 regions including the Zanzan region. The Zanzan region is in the northeast of Côte d'Ivoire between latitudes 7° and 10° North, and longitudes 4°30′ and 2°30′ West (Figure 1). This region, whose regional capital is the city of Bondoukou, is dominated by a humid tropical Sudanese-Guinean climate. The main climate has four seasons, two of which are rainy seasons from March to June and from September to October. These rainy seasons are interspersed with two dry seasons which run from November to February and from July to August. The annual rainfall varies between 800 and 1200 mm with an average temperature of about 28 °C [54]. The Zanzan region has a surface area of approximately 38,000 km 2 . Based on the 2014 national population and housing census, its population is estimated at 934,352 inhabitants [55]. The region's economy is mainly based on agriculture, livestock production, trade (food crops, handicrafts, etc.), and services. Cashew nuts are the main cash crop in the region with a contribution of 22% of national production [49]. The cashew industry constitutes 8.8% of the GDP of Côte d'Ivoire with a production of 700,000 tons [56]. In addition to this crop, cocoa and coffee are produced in small quantities in the southern part of the Zanzan region.

Data
The present research is based on a climatic database, spatially covering the entire Zanzan region. The data are a 40-year (1981 to 2020) time series of daily precipitation and temperature data from 12 measuring stations ( Figure 1) taken from the Côte d'Ivoire meteorological station (SODEXAM). The percentage of missing data was used to verify data

Data
The present research is based on a climatic database, spatially covering the entire Zanzan region. The data are a 40-year (1981 to 2020) time series of daily precipitation and temperature data from 12 measuring stations ( Figure 1) taken from the Côte d'Ivoire meteorological station (SODEXAM). The percentage of missing data was used to verify data quality. The percentages of missing data for minimum temperature (Tmin), maximum temperature (Tmax), and precipitation ranged from 0 to 6%, 0 to 7%, and 0 to 11%, respectively (Table 1). To fill in the missing data, the method of using a gridded meteorological database, based on reanalysis and satellite data was used in this study, especially the data from NASA Prediction of Worldwide Energy Resource (POWER). NASA POWER data are a suitable alternative when ground-based data are missing or of doubtful quality [57]. These data based on reanalysis models developed by NASA's Global Modeling and Assimilation Office (GMAO) [58] can be accessed and downloaded via the POWER website homepage at http://power.larc.nasa.gov (accessed on 11 March 2021). NASA POWER data have long been used and validated in the US in several studies [59]. Before using the NASA POWER data to fill in the missing data from the ground stations (SODEXAM), correlation tests were conducted between the different data sets using Pearson's correlation coefficient (r) for the twelve stations ( Table 2). The positive correlations closer to 1 illustrated by the results in Table 2 show a good correlation between the different data sets. It is at the end of the different tests that the POWER data were used to fill our missing data. Other studies have found a strong link between ground data and NASA POWER data, such as Bai et al. [60] in China, Bender and Sentelhas [61] in Brazil, and Aboelkhair et al. [62] in Egypt, or have used NASA POWER data to fill the gaps in their ground data, such Berhane et al. [63] in Ethiopia and West Africa [64]. Before calculating indices related to climate extremes, it is necessary to check the homogeneity of the daily input data, look for points of change in the time series, assess whether they are natural or not, and then rectify them. Homogeneity refers to a series' consistency through time, and it is an evident necessity for reliable climatic time series analysis [65]. There is presently no mechanism to homogenize the data in ClimPACT2, the software used to calculate the indices related to climatic extremes, except for simple quality control. The software application RHtest of the "Canadian Meteorological Service" and adopted by the ETCCDMI [32,66], was used in the current research work to correct the homogeneity of the daily data. This software is run in the R statistical and programming environment. A double regression model is used in this process [67].

Analysis of Climate Extremes Indices
For the characterization of the intensity, frequency, and duration of extreme climate events, indices have been developed and recommended by the Expert Team on Sectorspecific Climate Indices (ET-SCI) and the Expert Team on Climate Change Detection Monitoring and Indices (ETCCDMI), both from the WMO [65]. The ETCCDMI indices are integrated into the RClimDex software [68]. Those of ET-SCI, which already contain the ETCCDMI indices with additional indices to be used in sectoral applications based on specialists from the health, water, and agriculture sectors, are integrated into ClimPACT2 software [65]. ClimPACT2 software is an update of ClimPACT, which was based on ETCCDMI's RClimDex software. RClimDex and ClimPACT2 are written in R, statistical computing, and graphic language and environment. A series of 27 indices including 11 precipitation and 16 temperature-related indices were defined by ETCCDMI [65,69]. Climate indices are generated values that can be used to reflect the state of variations in the climate system. They allow statistical study and comparison of time series, extremes, trends, and means [70]. The use of indices to detect climate change has the advantage of allowing for a simple comparison of trends across different climate zones and areas [71]. In this study, the indices were calculated using the ClimPACT2 Master program after a quality control and homogeneity test of the daily input data, but only the indices deemed relevant and applicable to Côte d'Ivoire were selected for a detailed analysis ( Table 3). The quality control was performed by the software before the calculation of the climate indices in order to check for fundamental input mistakes, such as intentionally exaggerated results or rounding bias, and missing data. The indices listed in Table 2 were chosen for this study because they have already been applied several times in other studies in West Africa [27,[72][73][74] and Côte d'Ivoire [71,75,76] for the detection of change in climate extremes. These indices are therefore considered relevant and adaptable to the climatic context of the Zanzan region for the detection of changes in precipitation and temperature extremes.

Trend Analysis
After verification of data quality and coherency, the next step was to calculate the climate indices by ClimPACT2 software. Then a statistical analysis of the trends was performed. There are several methods of statistical trend analysis that other similar studies have already used, where we can mention among others the linear trend analysis [77], the linear trend analysis associated with the t-test [75,78], the non-parametric Kendall's tau test with linear least squares trend [33,79], the Sen's slope estimator with non-parametric Mann-Kendall test [63,74,80,81]. The non-parametric Mann-Kendall test has been utilized in various research works since it makes no assumptions regarding data distribution or trend linearity [82]. Unlike parametric trend tests, which demand that the data be both normally distributed and independent, non-parametric trend tests only require independent data [28]. The Sen's slope estimator and non-parametric Mann-Kendall (MK) test have been used in this study to assess the statistical significance of the trends and to accurately determine the magnitude of the trends in climate extremes indices due to their robustness in addition to their insensitivity to outliers in the times series. The MK test was performed using the MAKESENS software. The null hypothesis (H 0 ) in the calculation of this test asserts that there is no trend in climate extremes over time, whereas the alternative hypothesis (H 1 ) indicates that there is a monotonic decreasing or increasing trend. The mathematical equations used to calculate the Mann-Kendall test are as follows: In a data set of length n, X j and X k are consecutive data values.
The behavior of the S statistic is approximately the same for a normal distribution when n ≥ 8, so the test is performed with a normal distribution with E(S) = 0 with the variance as follows: Here t p is the number of input values inside the p-th affiliated group and q is the number of affiliated groups in the data set. The test statistic Z mk is calculated using the values of S and Var (S mk ) as follows: , i f S < 0 (4) Z mk is a statistically significant trend detection tool. A positive (negative) value of Z mk indicates the upward (downward) trend. To check for a monotonic upward or downward trend (two-tailed test) at the α level of significance, H 0 is rejected if the absolute value of is acquired from the usual normal cumulative distribution tables. The trend outcomes in this study were evaluated at the 5% level of significance, with the null hypothesis of no trend being rejected if |Z mk | > 1. 96. The estimator of Sen [83] which is defined by Equations (5) and (6), was used to calculate the real slope of an existing trend.
There will be N = n(n − 1)/2 slope Q i estimations if the time series has n x j values. The median of these N values of Q i is Sen's slope estimator.

Spatial Interpolation of Climate Indices
To produce a continuous surface, spatial distribution or interpolation improves the representation of a surface and forecasts values for other uncharted territories [84]. It is a mathematical approach or function for estimating values in situations when there are no measurable values [85]. Interpolation methods (e.g., spline, inverse distance weighting (IDW), Kriging, etc.) are numerous and vary greatly in complexity and effectiveness. The Kriging approach was employed in this study to interpolate the extreme climate indices from the ArcGIS software for spatial analysis. Kriging is a geostatistical method that provides unbiased estimates with minimum variance taking into account the spatial relationship between the data points [86]. The advantage of using geostatistical interpolation approaches over conventional techniques for spatial estimation of weather data has been reported by several authors [87][88][89]. One of the main advantages of Kriging over simpler techniques, such as IDW, besides providing a measure of estimation uncertainty (Kriging variance), is that it can use correlated dense secondary variables to improve the prediction of a sparsely sampled primary variable [86]. Considering the topography as a secondary variable, the prediction of meteorological data can be improved by using Kriging. In this study, among the interpolation methods, Kriging gives the best representativeness of our results.  Table 4 shows both the Sen slope and Mann-Kendall statistical test for all of the extreme rainfall indices studied for the 12 stations. In the Zanzan region, an examination of cumulative annual precipitations PRCPTOT shows a general decreasing trend. This is observed in 10 of the study stations, excluding the Doropo and Téhini sites, which show increasing trends. This downward trend is however statistically non-significant at the 5% level of significance. Over a decade, the largest decreases are estimated at 45.1 mm/decade at Kouassi-datékro, 41.5 mm/decade at Sandegué, and 38.2 mm/decade at Bondoukou. Over the period of 1981-2020, the annual total precipitations decreased by an average of 112.35 mm over the entire study area. The total annual precipitation index is shown in Figure 2 as a spatial distribution per decade. From 1981 to 2020, the cumulative annual rainfall ranged from 837.69 to 1404.3 mm. In the first decade (i.e., 1981 to 1990), the cumulative precipitation was between 947.1 and 1404.3 mm and from 1991 to 2000, the value was between 837.69 and 1276.9 mm, representing a decrease of 127.4 mm. For the 2001-2010 decade, the total annual precipitation varied between 852.14 and 1145.3 mm, for a decrease of 131.6 mm compared to the previous decade. From 2011 to 2020, the total precipitation was between 895.56 and 1327.2 mm, corresponding to an increase of 181.9 mm. Compared to the last two decades, the 2011-2020 decade shows a slight return of precipitation, but less important than the precipitation of the 1981-1990 decade. The  Figure 2 indicates that the largest overall precipitation was found in the Zanzan region's south, while the lowest was found in the study area's northern part.

Daily Precipitation Intensity Index (SDII)
The Mann-Kendall test of trends in the daily precipitation intensity index (Table 4) reveals general declining trends at most of the research sites, but are not significant at the 5% level of significance. The daily precipitation intensity in the studied area decreased by 0.47 mm/d on average across the study period. The largest drops were observed in Kouassi-datékro and Sandegué stations, with estimates of 1.17 mm/d and 0.78 mm/d, respectively. Figure 3 displays the spatial distribution of the index per decade. From 1981 to 1990, the daily rainfall intensity in the Zanzan region varied between 5.34 and 7.22 mm/d, and over the 1991-2000 decade this variation was between 4.82 and 6.59 mm/d, corresponding to a decrease of 0.63 mm/d. From 2001 to 2010, the values of the index were between 4.81 and 6.15 mm/d, which implies a decrease of 0.09 mm/d, and from 2011 to 2020 the variation was between 5.26 and 6.73 mm/d, which is an increase of 0.58 mm/d. From 1981 to 2020, the southwestern part of the Zanzan region experienced the highest daily precipitation intensity, while the north-eastern half experienced the lowest. 1990, the daily rainfall intensity in the Zanzan region varied between 5.34 and 7.22 mm/d, and over the 1991-2000 decade this variation was between 4.82 and 6.59 mm/d, corresponding to a decrease of 0.63 mm/d. From 2001 to 2010, the values of the index were between 4.81 and 6.15 mm/d, which implies a decrease of 0.09 mm/d, and from 2011 to 2020 the variation was between 5.26 and 6.73 mm/d, which is an increase of 0.58 mm/d. From 1981 to 2020, the southwestern part of the Zanzan region experienced the highest daily precipitation intensity, while the north-eastern half experienced the lowest. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo). Figure 4 and Table 4 show the outputs of the Mann-Kendall statistical test of the data series and the geographical extent of the consecutive wet day's index in the Zanzan region, respectively. All stations in the study area have a general dropping trend in consecutive wet days, and three of these stations have a significant declining trend at the 5% significance level, according to the descriptive statistics. Over the study period of 1981-2020, these significant decreases were estimated at 30.03 days at Bondoukou, 26.52 days at Sapli, and 22.62 days at Kouassi-datékro. Over the whole study area, the total average decline was 16.45 days. The spatial distribution of consecutive wet days per decade shows a variation from 37.   Table 4 show the outputs of the Mann-Kendall statistical test of the data series and the geographical extent of the consecutive wet day's index in the Zanzan region, respectively. All stations in the study area have a general dropping trend in consecutive wet days, and three of these stations have a significant declining trend at the 5% significance level, according to the descriptive statistics. Over the study period of 1981-2020, these significant decreases were estimated at 30.03 days at Bondoukou, 26.52 days at Sapli, and 22.62 days at Kouassi-datékro. Over the whole study area, the total average decline was 16.45 days. The spatial distribution of consecutive wet days per decade shows a variation from 37.  Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Index of Consecutive Dry Days (CDD)
In the Zanzan region, the trends in the time series of the consecutive dry days' index assessed using the Mann-Kendall test (Table 4) showed an increase in 10 of the stations tested, except for Koun Fao and Tanda, which showed non-significant declining trends. The almost general upward trend in consecutive dry days in the study area was however non-significant at the 5% significance level for most stations except for the Doropo station, which presented a significant increasing trend. In Doropo the increase was 0.89 days/year, 8.9 days per decade, and 34.71 days over the entire study period. On average, the increase was estimated at 18.67 days over the study period. Figure 5 presents the index of the con- Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Index of Consecutive Dry Days (CDD)
In the Zanzan region, the trends in the time series of the consecutive dry days' index assessed using the Mann-Kendall test (Table 4) showed an increase in 10 of the stations tested, except for Koun Fao and Tanda, which showed non-significant declining trends. The almost general upward trend in consecutive dry days in the study area was however nonsignificant at the 5% significance level for most stations except for the Doropo station, which presented a significant increasing trend. In Doropo the increase was 0.89 days/year, 8.9 days per decade, and 34.71 days over the entire study period. On average, the increase was estimated at 18.67 days over the study period. Figure Table 4 shows the findings of the Mann-Kendall statistical test, which demonstrated a general trend in the Zanzan region of increasing total yearly precipitation from very heavy rain days. At any of the sites studied, the general rise was non-significance at the 5% threshold of significance. Over the study period, the strong increases were estimated at 23.79 mm in Téhini, 10.14 mm in Bondoukou, and 8.97 mm in Tanda. Over the entire research area, the overall average was estimated at 4.61 mm. Figure 6 illustrates the spatial distribution of the total annual precipitations from very heavy rain days. Over the 1981-1990 decade, the evolution of extreme rainy days ranged from 37.4 to 98.5 mm, and from 1991 to 2000, it ranged from 34.1 to 66.9 mm, which represents a decrease of 31.6 mm. study witnessed a return to a significant quantity of extreme rainy days. These extremely wet days were more pronounced in the south of the Zanzan region across the study decades. The highest daily water intakes relative to the 99th percentile threshold are total annual precipitations from very heavy rain days (R99p).  Table 4 shows the findings of the Mann-Kendall statistical test, which demonstrated a general trend in the Zanzan region of increasing total yearly precipitation from very heavy rain days. At any of the sites studied, the general rise was non-significance at the 5% threshold of significance. Over the study period, the strong increases were estimated at 23.79 mm in Téhini, 10.14 mm in Bondoukou, and 8.97 mm in Tanda. Over the entire research area, the overall average was estimated at 4.61 mm. Figure 6 illustrates the spatial distribution of the total annual precipitations from very heavy rain days. Over the 1981-1990 decade, the evolution of extreme rainy days ranged from 37.4 to 98.5 mm, and from 1991 to 2000, it ranged from 34.1 to 66.9 mm, which represents a decrease of 31.6 mm.

Contribution of Extremely Wet Days (R99pTOT)
The contribution of extremely wet days analyzed using the Mann-Kendall test displays a general increasing trend that was non-significant at the 5% level of significance over the entire study area. The most significant increases were in the locality of Téhini estimated at 3.12%, Bondoukou estimated at 1.95% and Tanda estimated at 1.95% of days. On average, this increase was estimated at 1.41% of days from 1981 to 2020 for the entire region. The spatial distribution per decade of the contribution of extremely wet days illustrated in Figure 7 indicates a variation between 3.87% and 6.18% from 1981 to 1990, and from 1991 to 2000 the contribution was between 3.25% and 4.65%, hence a decrease of 1.53%. During the period 2001-2010, these values were between 2.82% and 3.09%, corresponding to a decrease of 1.56%, and from 2011 to 2020, the contribution of extremely wet days was between 7.58% and 10.08%, representing an increase of 6.99%. The last decade thus witnessed an important contribution of extremely wet days located in the east of the study area. Figure 6. Spatial distribution of the total annual precipitation from very heavy rain days (R99p) (mm) per decade in the Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Contribution of Extremely Wet Days (R99pTOT)
The contribution of extremely wet days analyzed using the Mann-Kendall test displays a general increasing trend that was non-significant at the 5% level of significance over the entire study area. The most significant increases were in the locality of Téhini estimated at 3.12%, Bondoukou estimated at 1.95% and Tanda estimated at 1.95% of days. On average, this increase was estimated at 1.41% of days from 1981 to 2020 for the entire region. The spatial distribution per decade of the contribution of extremely wet days illustrated in Figure 7 indicates a variation between 3.87% and 6.18% from 1981 to 1990, and from 1991 to 2000 the contribution was between 3.25% and 4.65%, hence a decrease of 1.53%. During the period 2001-2010, these values were between 2.82% and 3.09%, corresponding to a decrease of 1.56%, and from 2011 to 2020, the contribution of extremely wet days was between 7.58% and 10.08%, representing an increase of 6.99%. The last decade thus witnessed an important contribution of extremely wet days located in the east of the study area.

Contribution of Extremely Wet Days (R99pTOT)
The contribution of extremely wet days analyzed using the Mann-Kendall test displays a general increasing trend that was non-significant at the 5% level of significance over the entire study area. The most significant increases were in the locality of Téhini estimated at 3.12%, Bondoukou estimated at 1.95% and Tanda estimated at 1.95% of days. On average, this increase was estimated at 1.41% of days from 1981 to 2020 for the entire region. The spatial distribution per decade of the contribution of extremely wet days illustrated in Figure 7 indicates a variation between 3.87% and 6.18% from 1981 to 1990, and from 1991 to 2000 the contribution was between 3.25% and 4.65%, hence a decrease of 1.53%. During the period 2001-2010, these values were between 2.82% and 3.09%, corresponding to a decrease of 1.56%, and from 2011 to 2020, the contribution of extremely wet days was between 7.58% and 10.08%, representing an increase of 6.99%. The last decade thus witnessed an important contribution of extremely wet days located in the east of the study area.

Evolution of the Index of the Number of Cool Days (TX10p)
The Mann-Kendall statistical test and the Sen slope of all the indices of temperature extremes assessed for the 12 study sites are shown in Table 5. The analysis of the number of cool days indicates a general declining trend in the Zanzan region for all the stations studied. Within the general declining trend, eight of the studied stations showed significant trends at the level of 5%. The values of these declines over the 1981-2020 study period for the localities were estimated at 14.54% for Sapli, 9.98% for Bondoukou, 9.08% for Sandegué, 8.85% for Nassian, and 5.96% of days for Tanda sites. Generally, the average decline in the study area was 6.58% of days. To better appreciate this decreasing trend, Figure 8 illustrates the spatial distribution of the number of cool days per decade. From 1981 to 2020, the number of cool days was between 5.61% and 18.53% of days. In the first decade (i.e., 1981 to 1990), the number was between 15.57% and 18.53% and from 1991 to 2000 the values were between 6.92% and 10.51%, which means a decrease of 8.02%. Over the decade 2001 to 2010, the number of cool days ranged between 6.68% and 8.25%, representing a decrease of 2.26% compared to the previous decade. From 2011 to 2020 the number of days varied between 5.61% and 10.35%, corresponding to an increase of 2.1%. Compared to the last two decades, the decade 2011-2020 showed a slight return of cool days, but less than the number of cool days in the decade 1981-1990. The analysis of Figure 8 revealed that the high percentages of cool days in the decades 1981-1990 and 1991-2000 were concentrated in the center of the Zanzan region, but in the decades 2001-2010 and 2011-2020, they were concentrated in the study area's northern and southern parts.

Changes in the Index of the Number of Cold Nights (TN10p)
The Mann-Kendall test of trends in the index of the number of cold nights (Table 5) showed a general declining trend across all stations in the research area, which was significant at the 5% statistical significance. Over the research period from 1981 to 2020, the decline in the number of cold nights in the Zanzan region was estimated at an average of 11.83% of nights. The strongest decreases of cold nights were observed in Bondoukou (14.23%), Nassian (12.36%), Sandegué (12.48%), and Tanda (12.55%). Figure 9 shows the spatial distribution of the index per decade. From 1981 to 1990, the number of cold nights in the Zanzan region varied between 16.69% and 18.73%, and over the decade 1991-2000, the variation was between 8.79% and 11.52%, which represents a decrease of 7.21% of cold nights. From 2001 to 2010 the values of the index were between 6.18% and 7.18%, meaning a decrease of 4.34% and from 2011 to 2020, the variation was between 6.04% and 8.18%, indicating an insignificant increase of 0.99% of cold nights. From 1981 to 1990 and 1991 to 2000, the highest percentages of cold nights in the Zanzan region were recorded in the southeast and southwest of the region, respectively. From 2001 to 2010 and from 2011 to 2020, high percentages were seen in the northern part of Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Changes in the Index of the Number of Cold Nights (TN10p)
The Mann-Kendall test of trends in the index of the number of cold nights (Table 5) showed a general declining trend across all stations in the research area, which was significant at the 5% statistical significance. Over the research period from 1981 to 2020, the decline in the number of cold nights in the Zanzan region was estimated at an average of 11.83% of nights. The strongest decreases of cold nights were observed in Bondoukou (14.23%), Nassian (12.36%), Sandegué (12.48%), and Tanda (12.55%). Figure 9 shows the spatial distribution of the index per decade. From 1981 to 1990, the number of cold nights in the Zanzan region varied between 16.69% and 18.73%, and over the decade 1991-2000, the variation was between 8.79% and 11.52%, which represents a decrease of 7.21% of cold nights. From 2001 to 2010 the values of the index were between 6.18% and 7.18%, meaning a decrease of 4.34% and from 2011 to 2020, the variation was between 6.04% and 8.18%, indicating an insignificant increase of 0.99% of cold nights. From 1981 to 1990 and 1991 to 2000, the highest percentages of cold nights in the Zanzan region were recorded in the southeast and southwest of the region, respectively. From 2001 to 2010 and from 2011 to 2020, high percentages were seen in the northern part of Zanzan region.  Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Changes in the Index of the Number of Cold Nights (TN10p)
The Mann-Kendall test of trends in the index of the number of cold nights (Table 5) showed a general declining trend across all stations in the research area, which was significant at the 5% statistical significance. Over the research period from 1981 to 2020, the decline in the number of cold nights in the Zanzan region was estimated at an average of 11.83% of nights. The strongest decreases of cold nights were observed in Bondoukou (14.23%), Nassian (12.36%), Sandegué (12.48%), and Tanda (12.55%). Figure 9 shows the spatial distribution of the index per decade. From 1981 to 1990, the number of cold nights in the Zanzan region varied between 16.69% and 18.73%, and over the decade 1991-2000, the variation was between 8.79% and 11.52%, which represents a decrease of 7.21% of cold nights. From 2001 to 2010 the values of the index were between 6.18% and 7.18%, meaning a decrease of 4.34% and from 2011 to 2020, the variation was between 6.04% and 8.18%, indicating an insignificant increase of 0.99% of cold nights. From 1981 to 1990 and 1991 to 2000, the highest percentages of cold nights in the Zanzan region were recorded in the southeast and southwest of the region, respectively. From 2001 to 2010 and from 2011 to 2020, high percentages were seen in the northern part of Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Evolution of the Index of Cold Spell Duration Indicator (CSDI)
The Mann-Kendall test demonstrated a general declining trend in the index of cold spell duration indicator (Table 5) for the entire research area. At the 5% significance level, the general decreasing trends were significant at all study sites. Over the research period of 1981-2020, the largest significant declines were in Bondoukou estimated at 39 days, Kouassi-datékro estimated at 30.57 days, Nassian estimated at 29.99 days, and Sandegué estimated at 29.95 days. On average, this decline was estimated at 26.15 days from 1981 to 2020 for the entire region. The spatial distribution per decade of the cold spell duration indicator shown in Figure 10 indicates that from 1981 to 1990, the variation was between 36.65 and 46.97 days, and from 1991 to 2000 the variation was between 14.43 and 25.8 days, representing a decrease of 21.17 days. During the period 2001 to 2010, the values were between 10.21 and 14.1 days, which was a decrease of 11.7 days, and from 2011 to 2020, the variations were between 11.8 and 19.59 days, implying an increase of 5.49 days. In comparison with the preceding two decades, the cold spell duration indicator in the north of the region under consideration has made a slight return in the last decade.
The Mann-Kendall test demonstrated a general declining trend in the index of cold spell duration indicator (Table 5) for the entire research area. At the 5% significance level, the general decreasing trends were significant at all study sites. Over the research period of 1981-2020, the largest significant declines were in Bondoukou estimated at 39 days, Kouassi-datékro estimated at 30.57 days, Nassian estimated at 29.99 days, and Sandegué estimated at 29.95 days. On average, this decline was estimated at 26.15 days from 1981 to 2020 for the entire region. The spatial distribution per decade of the cold spell duration indicator shown in Figure 10 indicates that from 1981 to 1990, the variation was between 36.65 and 46.97 days, and from 1991 to 2000 the variation was between 14.43 and 25.8 days, representing a decrease of 21.17 days. During the period 2001 to 2010, the values were between 10.21 and 14.1 days, which was a decrease of 11.7 days, and from 2011 to 2020, the variations were between 11.8 and 19.59 days, implying an increase of 5.49 days. In comparison with the preceding two decades, the cold spell duration indicator in the north of the region under consideration has made a slight return in the last decade. Figure 10. Spatial distribution of cold spell duration indicator (CSDI) (%) per decade in Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo). Table 5 and Figure 11 give the outputs of the Mann-Kendall statistical test and the spatial distribution of the number of hot days' index in the Zanzan region, respectively. In all of the stations within the research area, the statistical analysis provides an overall increase in the number of hot days, and two of these stations show significant increasing trends at the 5% significance level. Over the study period of 1981-2020, the significant increases were estimated at 13.18% of days at Sapli, and 7.44% of days at Sandegué. The overall average increasing trend was 5.19% of days across the entire research area. The spatial distribution of the number of hot days per decade showed a variation from 5.74% to 9.11% over the period of 1981-1990, and from 1991 to 2000, the index varied from 8.07% to 12.3%, which represents an increase of 3.19% of days. The highest index values over the last two decades have been observed in the northern and southern parts of the study area. From 2001 to 2010 the variation of the index was estimated between 12.2% and 15.4%, Figure 10. Spatial distribution of cold spell duration indicator (CSDI) (%) per decade in Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo). (TX90p)   Table 5 and Figure 11 give the outputs of the Mann-Kendall statistical test and the spatial distribution of the number of hot days' index in the Zanzan region, respectively. In all of the stations within the research area, the statistical analysis provides an overall increase in the number of hot days, and two of these stations show significant increasing trends at the 5% significance level. Over the study period of 1981-2020, the significant increases were estimated at 13.18% of days at Sapli, and 7.44% of days at Sandegué. The overall average increasing trend was 5.19% of days across the entire research area. The spatial distribution of the number of hot days per decade showed a variation from 5.74% to 9.11% over the period of 1981-1990, and from 1991 to 2000, the index varied from 8.07% to 12.3%, which represents an increase of 3.19% of days. The highest index values over the last two decades have been observed in the northern and southern parts of the study area. From 2001 to 2010 the variation of the index was estimated between 12.2% and 15.4%, meaning an increase of 3.1% and from 2011 to 2020, the variation was between 6.72% and 14.5%, meaning a small decrease of 0.9% of days. The decade 2001-2010 appears to be the hottest in the spatial analysis of the number of relatively hot days, and this significant heating was more pronounced in the south-western part of the Zanzan region.

Evolution of the Index of the Number of Warm Nights (TN90p)
The trends in the temporal series of the index of the number of warm nights analyzed through the Mann-Kendall test (Table 5) in the Zanzan region revealed a general increase over all the stations studied. Except for the Sapli station, which indicated a non-significant rising trend, the general rising trend in the number of warm nights in the study region was significant at the 5% level of statistical significance for 11 of the study stations. Geographically over the study period, these increases were estimated at 14.89% in Bondoukou, 14.27% in Sandegué, 13.68% in Nassian, and 13.80% of nights in Tanda. On average in the study area, this increase was estimated at 11.50% of nights. The spatial distribution of the warm nights' index per decade is presented in Figure 12. Over the decade 1981-1990, the number of warm nights was between 4.84% and 7.77% and from 1991 to 2000 it varied between 8.1% and 10.25%, which is an increase of 2.48%. From 2001 to 2010, the variation in the number of warm nights was estimated between 11.19% and 14.73%, representing an increase of 4.48%, and from 2011 to 2020, the variation was between 10.44% and 14.56%, indicating a small decrease of 0.17% in warm nights. The comparison between the first decade of this study (i.e., [1981][1982][1983][1984][1985][1986][1987][1988][1989][1990]  meaning an increase of 3.1% and from 2011 to 2020, the variation was between 6.72% and 14.5%, meaning a small decrease of 0.9% of days. The decade 2001-2010 appears to be the hottest in the spatial analysis of the number of relatively hot days, and this significant heating was more pronounced in the south-western part of the Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Evolution of the Index of the Number of Warm Nights (TN90p)
The trends in the temporal series of the index of the number of warm nights analyzed through the Mann-Kendall test (Table 5) in the Zanzan region revealed a general increase over all the stations studied. Except for the Sapli station, which indicated a non-significant rising trend, the general rising trend in the number of warm nights in the study region was significant at the 5% level of statistical significance for 11 of the study stations. Geographically over the study period, these increases were estimated at 14.89% in Bondoukou, 14.27% in Sandegué, 13.68% in Nassian, and 13.80% of nights in Tanda. On average in the study area, this increase was estimated at 11.50% of nights. The spatial distribution of the warm nights' index per decade is presented in Figure 12. Over the decade 1981-1990, the number of warm nights was between 4.84% and 7.77% and from 1991 to 2000 it varied between 8.1% and 10.25%, which is an increase of 2.48%. From 2001 to 2010, the variation in the number of warm nights was estimated between 11.19% and 14.73%, representing an increase of 4.48%, and from 2011 to 2020, the variation was between 10.44% and 14.56%, indicating a small decrease of 0.17% in warm nights. The comparison between the first decade of this study (i.e., [1981][1982][1983][1984][1985][1986][1987][1988][1989][1990]   (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Evolution of the Warm Spell Duration Indicator (WSDI)
The data in Table 5 on the outputs of the Mann-Kendall statistical test showed a general increasing trend of the warm spell duration indicator in the Zanzan region. The general rising trend was nevertheless non-significant at the 5% significance level, except for the Sapli and Sandegué stations, which displayed significant rising trends. Over the study period, these significant increases were estimated at 24.29 days in Sapli, 14.39 days in Sandegué, and the overall average increase was estimated at 9.71 days over the entire study area. Figure 13 illustrates the spatial pattern of the warm spell duration indicator. Over the decade 1981-1990, the evolution of the warm spell duration indicator was between 13.5 and 26.7 days, and from 1991 to 2000, the evolution was between 16.8 and 31.3 days, meaning an increase of 4.6 days. Over the first decade, the strong warm spells were more intense in the southern and northern part of the research area, and over the second decade 1991-2000, the strong warm spells were more noticeable in the southern part of the region.

Extreme Precipitations
This study carried out in the Zanzan region showed a decrease in total annual precipitation (PRCPTOT) at most of the stations studied (10 of the stations), with an average decrease of 112 mm over the research period. Furthermore, a generalized decrease in daily precipitation intensity SDII was seen across the entire study area, with the example of the localities in the southeast of the region, namely, Kouassi-datékro and Sandegué with respective decreases estimated at 1.17 mm/d and 0.78 mm/d. The decrease of the PRCPTOT index was also observed in the spatial distribution following the decades even if in the last decade we notice a slight return of precipitations. The primary influence of the decrease in precipitation was observed through the low intensity of rainfall received by the study area, as the daily precipitation intensity was also decreasing even though the decrease seemed insignificant. The increasing trend in daily precipitation intensity over the last decade had also led to a slight overall increase in rainfall. The general decrease in rainfall in the study area followed the general trend of decreasing rainfall that is observed throughout Côte d'Ivoire and varies according to the regions of the country. Indeed, for the PRCPTOT index, Allechy et al. [76] in the central-western part of Côte d'Ivoire noted a decrease of 788 mm between 1984 and 2013, Atcheremi et al. [71] in the south-western Figure 13. Spatial distribution of warm spell duration indicator (WSDI) (days) per decade in Zanzan region. Stations (S1 = Koun Fao, S2 = Kouassi-datékro, S3 = Tanda, S4 = Sandegué, S5 = Bondoukou, S6 = Nassian, S7 = Sapli, S8 = PN Comoé S1, S9 = Bouna, S10 = Téhini, S11 = Doropo).

Extreme Precipitations
This study carried out in the Zanzan region showed a decrease in total annual precipitation (PRCPTOT) at most of the stations studied (10 of the stations), with an average decrease of 112 mm over the research period. Furthermore, a generalized decrease in daily precipitation intensity SDII was seen across the entire study area, with the example of the localities in the southeast of the region, namely, Kouassi-datékro and Sandegué with respective decreases estimated at 1.17 mm/d and 0.78 mm/d. The decrease of the PRCPTOT index was also observed in the spatial distribution following the decades even if in the last decade we notice a slight return of precipitations. The primary influence of the decrease in precipitation was observed through the low intensity of rainfall received by the study area, as the daily precipitation intensity was also decreasing even though the decrease seemed insignificant. The increasing trend in daily precipitation intensity over the last decade had also led to a slight overall increase in rainfall. The general decrease in rainfall in the study area followed the general trend of decreasing rainfall that is observed throughout Côte d'Ivoire and varies according to the regions of the country. Indeed, for the PRCPTOT index, Allechy et al. [76] in the central-western part of Côte d'Ivoire noted a decrease of 788 mm between 1984 and 2013, Atcheremi et al. [71] in the south-western part of Côte d'Ivoire revealed a decrease of 85 mm in Daloa and 493 mm in Sassandra as well as a decrease of 159 mm according to Balliet et al. [75] in Gagnoa. In the forested south of Côte d'Ivoire, N'guessan et al. [90] determined a decrease in PRCPTOT of 20 mm followed by a decrease of 0.5 mm/d in rainfall intensities over two decades. The generalized decline in rainfall, which had even created a major drought wave starting in the 1970s, as revealed in the study conducted by Paturel et al. [91], continued over the following decades in Côte d'Ivoire, particularly in the Zanzan region, even though a slight return of rainfall was observed over the last decade. However, this decline in total annual precipitation and daily precipitation intensity is more than problematic for agricultural development, especially for the main cash crop in the Zanzan region, cashew nuts. For acceptable cashew nut production, rainfall between 800 and 1800 mm is required to be considered sufficient [92], whereas in the Zanzan region this rainfall, which varies between 800 and 1400 mm, is decreasing. A future study of the impact of rainfall variations on agricultural products, especially cashew nuts, will provide information on the real economic threat in the Zanzan region. The observation of rainfall disruption is also shown elsewhere in other studies. In the West African sub-region, studies such as that of New et al. [33] showed that the regions as a whole experienced a decrease in total annual precipitation and the number of rainy days. Thus, Hountondji et al. [72] in Benin found a decrease in total annual precipitation PRCPTOT for 95% of the stations studied and stability of the daily precipitation intensity SDII. This decrease in rainfall extends even across Africa as Vondou et al. [74], reported that the cumulative annual rainfall and the daily precipitation intensity SDII are decreasing in northern Cameroon. Unlike the decrease in the daily precipitation intensity SDII found in our study area, Gebrechorkos et al. [81] in their work found an increase in the SDII index in Kenya.
The CWD index for consecutive wet days was declining by an average of 16.45 days across all of the stations analyzed in the Zanzan region, while the CDD index for consecutive dry days was rising by an average of 18.67 days throughout most of the stations studied, according to the findings of this study. The CDD index is a marker of the presence of drought in a region. Its increase in a locality should be a concern, especially in an area such as the Zanzan region where the source of income of the population is based on agriculture. It is also a factor in the threat to water resources. These indicators have been studied in other regions of Côte d'Ivoire. In the southwest of Côte d'Ivoire, contrary to our result, Atcheremi et al. [71] observed an increase in consecutive wet days in Daloa, but in Gagnoa in the same area, Balliet et al. [75] indicated that the CWD was decreasing while the CDD was increasing. Furthermore, in the west-central area of Côte d'Ivoire, Allechi et al. [76] demonstrated that the CWD index has decreased while the CDD index has increased by 0.812 days/year. The occurrence of dry sequences within the rainy season thus seems to be widespread in Côte d'Ivoire, which is a cause for concern in an essentially agricultural country. For a proper practice of rainfed agriculture, when during the rainy season there are sequences of dry days, this can cause the drying of the young plant at the primary stage of its growth, especially for crops such as maize and yam, which are cultivated in the Zanzan region. This impact may be minimal on cashew trees because they have more or less short-drought resistant characteristics. In the West African sub-region, this observation of the evolution of these indices seems to be the same. The study in Ghana by Larbi et al. [93] showed rising trends in CDD and a declining trend in the CWD index, and Ozer et al. [94] identified the same results in Niger. These rainfall pattern parallels would be worldwide in scope, as they have been confirmed in several parts of the Pacific [31], and Donat et al. [95] detected that the CDD index was trending toward arid conditions across South Africa, East Asia, and South America in their analytical research of precipitation and temperature extremes since the turn of the 20th century. The R99p and R99pTOT indices in this study were slightly increased in the Zanzan region. The average increase in total annual precipitation from very heavy rain days was 1.41% of days. The spatial distribution clearly showed this evolution over the decades. The last decade, 2011-2020, showed a significant increase for both indices. A rise in the R99p and R99pTOT indices is a very important water contribution. For rainfed agriculture, these high water inputs can contribute to the development of the plant, but only if the plant has already exceeded the primary growth period, at the risk of being destroyed by these heavy rains, as could be the case for maize or cashew. These indices are sometimes the cause of flooding and landslides. Balliet et al. [75] identified a rise in R99p index at Gagnoa in Côte d'Ivoire. These patterns of change in these indices are similar and sometimes different in other parts of Africa. Hountondji et al. [72] found more stable conditions of the R99p and R99pTOT index in Benin. In the Central Sahel, Panthou et al. [73] reported that the occurrence of extreme rainfall was increasing as well as in the West African Sahel [96]. In contrast to the increase in R99p index found in this study, Berhane et al. [63] found in their study a significant decrease of R99p in Ethiopia in the Tigray region. With the advent of climate change, some areas of the world will receive more rainfall, and sometimes the scarcity of rainfall in some areas will increase, hence the similarity and differences in the results of this study compared to others. In the Zanzan region, these high rainfall amounts have certainly influenced the slight return of rainfall observed in the decade 2011-2020.

Temperature Extremes
A considerable drop in the trends of the cold spells, such as the index of the number of cold nights TN10p, the index of the number of cool days TX10p, and the cold spell duration indicator CSDI, was noted at all research stations in the Zanzan region in the north-east of Côte d'Ivoire. On the other hand, the warm spells, such as the index of the number of hot days TX90p, the index of warm spell duration indicator WSDI, and the index of the number of warm nights TN90p increased in all the stations studied and were more significant in the stations such as Bondoukou, Sapli, and Sandegué. Regarding the spatial distribution per decade, the cold and warm spells have evolved in contrast. The first decade (i.e., [1981][1982][1983][1984][1985][1986][1987][1988][1989][1990] showed a large number of cool days and this has evolved in a decreasing way until the last decade (i.e., 2011-2020) which showed a very low number of cool days. For the warm spells, the evolution started with a low number of hot days before displaying a high number of hot days over the last two decades. However, the decade 2001-2010 has shown to be the warmest in terms of days and the last two decades in terms of nights in the study area. The humid climate of the Zanzan region has been transformed into a hot climate. In the climatic context of Côte d'Ivoire, Atcheremi et al. [71] found similar results concerning the increase in the WSDI index and the decrease in a CSDI index in the southwest of Côte d'Ivoire, precisely in Daloa, Gagnoa, and Sassandra. The various trends observed related to the warming of the study area are under the influence of climate change. For the Zanzan region, this situation is a concern. The shortening of the rainy seasons and the warming of the study area will have an influence on water resources in the region and will greatly impact agriculture, which is essentially rainfed in this area. For cashew farming, for example, the presence of high temperatures during the flowering period can be very problematic. The threshold temperature should not be above 32 • C in order not to compromise the flowering and fruiting of the cashew tree. Indeed, beyond this threshold value, the temperature can provoke scald and flower abortion phenomena and significantly reduce the fruiting of the cashew tree [97]. The lack of cold nights in favor of the warming could impact the health of the population as shown in other studies, namely, Ambouta et al. [98] in Niger and Abatan et al. [26] in Nigeria. The climate context related to global warming is a sub-regional observation. Barry et al. [27] found in West Africa that warm nights and hot days are becoming more common, whereas cold nights and cool days are becoming less common. Furthermore, Chaney et al. [99] concluded that in Sub-Saharan Africa, there is a more rapid rise in warm climatic occurrences and a decline in cold spells. This warming is even spreading across Africa and around the world. In a study conducted in the "Greater Horn of Africa" in countries such as Burundi, Sudan, Rwanda, Djibouti, Tanzania, Uganda, and Ethiopia, Omondi et al. [100] stated that there was a rise in extreme heat (TX90P, WSDI) while cold spells decreased. In South Africa [101] and Morocco [102,103], the scholars concluded that there was an increase in warm extreme events TX90p and WSDI while cold extreme conditions TX10p and CSDI declined. Popov et al. [104] evaluated the extreme temperature indices in Bosnia and Herzegovina and concluded that there was a significant rise in warm temperature (TN90p, TX90p, WSDI) and similar observations were made in Israel and Palestine by Salameh et al. [105], as well as in China by Zhou et al. [106]. The decade 2001-2010 was found to be warmer than other decades in this study. This increase in warmer conditions during the 2001-2010 decade according to the WMO [107] and Coumou and Rahmstorf [108] studies was felt globally in several regions of the world. According to the IPCC [31], the temperature change follows natural and anthropogenic factors, especially in the causal relationships between the emissions of greenhouse gas (GHG). For the IPCC, continued GHG emissions at the current rate or higher are expected to increase warming and profoundly alter the climate system in the 21st century. It should be noted that the decrease in rainfall and the temperature rise are also locally linked to the regression of forests following the ongoing deforestation in Côte d'Ivoire. Indeed, according to Yao [109], the Ivorian forest has declined from 13 million hectares in colonial times to less than two million hectares currently.

Conclusions
The current study, whose major objective was to spatiotemporally assess the climate change in precipitation and temperature extremes in the Zanzan region using ETCCDMI indices, has produced a summary of the region's extreme climatic trends. Daily observations of minimum and maximum temperatures, as well as daily rainfall, were used to calculate the indices. For the extreme precipitation trends, the stations studied showed indices that have mostly non-significant statistical trends. Indeed, the results showed a decline in annual total precipitation PRCPTOT, daily precipitation intensity SDII, and consecutive wet days CWD, indicating a decrease in rainfall in the region. These decreases were estimated to be an average over the study period at 112.35 mm for the PRCPTOT index and 16.45 days for the CWD index. Following that, the CDD (consecutive dry days) was increased, as well as the indices R99p and R99pTOT (total annual precipitation from very heavy rain days and contribution from extremely wet days, respectively). The latter indices can lead to flooding and sometimes landslides. The increase, for example, in the CDD index, was estimated at 18.67 days over the study period. The indices related to extreme temperatures were much more statistically significant than the indices related to extreme precipitation. For these temperature extremes, we observed a decline in the cold spells TX10p, TN10p, CSDI and a rise in the warm spells TX90p, TN90p, and WSDI. We can state that the decreases in the cold spells were estimated on average over the study period at 6.58% of the days for the TX10p index and 26.15 days for the CSDI index. As for the warm spells, the increases were estimated on average at 5.19% of days for the TX90p index and 9.71 days for the WSDI index. Decreased rainfall and increased temperatures in this agricultural region of Côte d'Ivoire may influence agricultural production through low yields and could also threaten water resources. Furthermore, the socio-economic situation of this region would be even more alarming if there is an intensification of the currently observed trends in the future regarding the practice of crops, such as cashew, which would provoke a vulnerable health situation for the population. The knowledge of these different trends in temperature and rainfall extremes for this region is of high importance in the measure of sustainable agricultural and economic planning by the government authorities. Measures to encourage the population to adopt crops that are more resistant to the current climatic conditions, to think about better planning of the agricultural calendar according to the rainfall conditions, and consider forest management can immediately be set up by the authorities of the Zanzan region.