Solar Irradiance and Temperature Variability and Projected Trends Analysis in Burundi

: This paper assessed the variability and projected trends of solar irradiance and temperature in the East of Burundi. Observed temperature from meteorological stations and the MERRA-2 data set provided by NASA / Goddard Space Flight Center are used over the historical period 1976–2005. In addition, solar irradiance data provided by SoDa database were considered. Furthermore, projection data from eight Regional Climate Models were used over the periods 2026–2045 and 2066–2085. The variability analysis was performed using a standardized index. Projected trends and changes in the future climate were respectively detected through Mann-Kendall and t-tests. The ﬁndings over the historical period revealed increase temperature and decrease in solar irradiance over the last decades of the 20th century. At a monthly scale, the variability analysis showed that excesses in solar irradiance coincide with the dry season, which led to the conclusion that it may be a period of high production for solar energy. In the future climate, upward trends in temperature are expected over the two future periods, while no signiﬁcant trends are forecasted in solar irradiance over the entire studied region. However, slight decreases and signiﬁcant changes in solar irradiance have been detected over all regions. both mean at EAP, ELL and NLL. downward


Introduction
All over the world, energy is indispensable for economic and socio-development activities. Many African countries, especially those in the process of development, are not able to produce enough energy, less expensive and which respects the standards of environment protection despite the potential of the continent. This insufficiency of energy causes the poverty and under development activities for smallholder peoples. In developing countries, populations in generally use woods energy that contributes to the destruction of the environment [1]. On the other hand, due to low adaptive capacity and high sensitivity of socio-economic systems, Africa is one of the most vulnerable regions highly affected and to be affected by the impacts of climate changes [2]. The renewable energies in general, considered as sustainable and cleaned energies, must be one of the solutions of this environmental and socio-economic challenges. Furthermore, solar energy is an important renewable energy source, both in the generation of PV electricity and as heat [3]. Indeed, the photovoltaic power (PV power) in particularly is the fastest growing electricity generation technology in the world [4]. Solar radiation is the principal source of energy for the climate system. Variability in the amount of energy received at the surface of the Earth, therefore, has implications for climate changes at global, regional and local scales, as well as for water resource, agriculture, architectural design, solar thermal devices [5].
The most important variables in solar energy are solar irradiance and surface air temperature compared with other climatic and environmental factors [6]. Then changes of solar radiation intensity are proportionally related to PV power output, whereas increasing air ambient temperature negatively affects PV energy output [7]. A good understanding of climate changes and projected trends of these two variables is very important for policymakers and/or industries managers to undertake planning projects of installing solar cells in a considered locality. Many studies have been carried out for highlighting the changes, trends and variability analysis of the solar irradiance and temperature.
At a global scale, some authors [8] assessed 20th century changes in surface solar irradiance in simulations and observations. They analyzed 20th Century simulations using nine state-of-the-art climate models and showed that all models estimated a global annual mean reduction in downward surface solar radiation of 1-4 W/m 2 while the globe warms by 0.4-0.7 • C. Other available studies documented a widespread reduction (3-9 W/m 2 ) of surface solar radiation (SSR) from the 1950s-1980s called "global dimming" [9] and a subsequent increase (1-4 W/m 2 ) since the 1980s called "brightening" [10,11]. These variations are mostly due to changes in the transparency of the atmosphere due to variations in cloudiness and/or changes in anthropogenic aerosol emissions [12,13]. Indeed Dhomse et al. [14] conducted a study on Estimates of ozone return dates from Chemistry-Climate Model Initiative simulations which showed the impacts of greenhouse gases emission (anthropogenic aerosols emission) on stratospheric ozone layer depletion and concluded that different RCP scenarios may affect surface UV radiation and therefore the population.
A study conducted by J.A. Crook [15] assessed the climate changes impacts on future photovoltaic (PV) and concentrated solar power (CSP) energy output. The results indicated that PV output from 2010-2080 is likely to increase by a few percent in Europe and China, see little change in Algeria and Australia, and decrease by a few percent in Western USA and Saudi Arabia. CSP energy output is likely to increase by more than 10% in Europe, increase by several percent in China, and a few percent in Algeria and Australia, and decrease by a few percent in Western USA and Saudi Arabia. Other studies conducted in West Africa, References [6,16], analyzed the energy output considering the impact of climate changes, variability and trends of solar radiation intensity and surface air temperature. They found that the trends of solar irradiation are negative and photovoltaic power output is estimated to decrease in all the parts of West Africa. A study conducted in Australia [17] assessed the Interannual Variability of the Solar Resource across Australia and their results indicated that the Southeastern and central coastal regions of Queensland experience the highest levels of interannual variability, ranging between 3.5% and 6.5% for global horizontal irradiance.
At regional and local scales, there are no many studies conducted in East Africa related to solar radiation features. However, spatial and temporal variability in global, diffuse, and horizontal direct irradiance and sunshine duration has been evaluated at eight stations in South Africa and two stations in Namibia, where the time series range between 21 and 41 years. Global and direct irradiance and sunshine duration decrease from Northwest to Southeast; diffuse irradiance increases toward the East [5]. In Kenya, a study conducted in 2013 investigated the potential of solar energy as a local source of clean and renewable resource for Nakuru city. Results revealed that Nakuru is a moderate to high solar energy potential region. The study concluded that Nakuru is endowed with abundant solar energy resources, favorable for tapping at both small and medium scale levels [18].
Concerning the surface air temperature, due to the global warming of the earth in these last decades, many researchers have been interested in studying this climatic parameter across the world and/or at the local scale. All results agreed that temperature have increased by 0.6 • C-0.8 • C during 20th century [19]. Studies conducted in Ethiopia [20,21] assessed trends and variability of rainfall and temperature in some regions and revealed an increasing trend of annual maximum, minimum and mean temperature and an increase of 1.08 • C was observed. Trend analysis of mean annual temperature conducted in Rwanda for fifty-two years revealed for all observatories a significant warming trend over the period after 1977-1979 where Kigali the capital of Rwanda presented the highest values of the slope [22].
In Burundi, some studies have been conducted on surface air temperature, rainfall, and wind speed in some regions of the country. Thus, recent studies carried out in Western highlands and lowlands of Burundi presented an upward trend of ambient temperature in last half of 20th century and projected increasing trends in temperature by 2050 using Regional Climate Models [23]. Burundi has also gone through many catastrophic events and climatic disasters which caused many problems to the local population [24]. Indeed, Burundi and especially the North region have been hit by the socio-political conflicts since 1988 which led to the creation of refugees' sites. The livelihoods of this peoples' grouping and demographic increasing induced the increase need forwoods, crops, and then led to the deforestation in addition to water scarcity in this region. In addition to the global warming and according to the study conducted by Calabrò & Magazù [25], the local climate change may have been connected to change of these environments and demographic factors which can lead to desertification. In fact, the drought burst in the Northeast of Burundi where the smallholder peoples were severely hit by famine following this dryness. Furthermore, another main consequence of climate changes in Burundi is a chronic shortage and deficit of electrical energy. Indeed, energy resources are not enough to deal with increased demand for energy, and these resources are mainly of hydroelectric origins. However, these hydropower plants do no longer produce the amount of energy expected during the installation due to the significant decrease in the water level in the reservoir [23]. This reduction of water in the retained lakes is mainly due to the climatic changes observed in recent years, for example, the decrease in rainfall and the extension of dry season. In the East and Northeast of Burundi there are no important hydropower plants installed which can produce sufficient energy regarding the need of local population. These challenges lead to the need of developing others forms of renewable energy resources like solar energy especially in the Eastern region of Burundi. In this goal, researches conducted in this area are needed to highlight climate changes and variability of the main factors of solar energy to predict the eventual potential of solar energy production. Our study aims to investigate two climatic parameters, surface air temperature and solar irradiance over three Eastern regions of Burundi. The main objective of our research is to analyze the variability, changes and projected trends of solar irradiance, and surface air temperature to predict the eventual potential availability of solar energy in this part of the country. The study is conducted in Eastern lowlands, Northern lowlands and Eastern arid plateaus of Burundi. Specifically, the study analyzes the variability of these two parameters in the historical period and detects their trends and changes considering two future periods, namely the near future (2026-2045) and far future (2066-2085), referring to the baseline period (1986-2005).

Study Area
The study area is the East and Northern East of Burundi presented in Figure 1. Burundi is a small country located in East Africa between longitudes 28.8 • -30.9 • East and latitudes 2.3 • -4.45 • South [26]. Burundi is bounded in the North by the Republic of Rwanda, in the West by the Democratic Republic of Congo, and in the East and South by the Republic of Tanzania. The area of Burundi is 27,834 km 2 and it belongs to two major watersheds: The Nile basin with an area of 13,800 km 2 and Congo River basin with 14,034 km 2 of area [27]. The climate in Burundi, as an East African country, is mainly influenced by the North-South movement of the Intertropical Convergence Zone (ITCZ), the topography of the country and El-Nino Southern oscillation (ENSO) [23]. Then the annual mean of climate variables such as precipitation, temperature and wind speeds, depends mainly to the climate zone location. The interest study area covers the East and Northeast of Burundi and is divided into three sub-regions considering their landforms: The Eastern lowlands of Kumoso (ELL), the Eastern arid plateaus (EAP) and the Northern lowlands of Bugesera (NLL). The ELL is a region located to the extension of the border with Republic of Tanzania whereas NLL is covering two provinces (Kirundo, Muyinga) located Northeast of the country. The main different between the three regions is elevation as shown in Figure 1. This elevation defines the climate of each region and then temperature is high in lowlands regions. The annual mean temperature varies from 19 • C for EAP to 22 • C for ELL and NLL, and their altitudes also vary from 1400-2000 m for EAP to 1100-1400 m for ELL and NLL [27]. The mean yearly solar irradiance for these regions is closely the same and varies from 232 W/m 2 -248 W/m 2 .
Climate 2019, 7, x FOR PEER REVIEW 4 of 21 covering two provinces (Kirundo, Muyinga) located Northeast of the country. The main different between the three regions is elevation as shown in Figure 1. This elevation defines the climate of each region and then temperature is high in lowlands regions. The annual mean temperature varies from 19 °C for EAP to 22 °C for ELL and NLL, and their altitudes also vary from 1400-2000 m for EAP to 1100-1400 m for ELL and NLL [27]. The mean yearly solar irradiance for these regions is closely the same and varies from 232 W/m 2 -248 W/m 2 .

Data Used
Data of two climate parameters are considered in this study: Surface air temperature and solar irradiance as observed and projected data. The observed data used are grouped into three sets of data. The first set of data is a collection of observed data from the synoptic stations of Burundi belongs to the Geography Institute of Burundi (IGEBU). Monthly and annual mean temperature data for seven stations located in the study area are used over the historical period of thirty years . These climate stations are presented in Figure 1 and their geographic information and description are given in Table 1.
However, some stations had uncompleted data and shorter duration than the 30 years recommended. To complete observed data, the second set of data used is MERRA-2data set available at http://gmao.gsfc.nasa.gov/reanalysis/MERRA-2 [28] and provided by NASA/Goddard Space Flight Center.

Data Used
Data of two climate parameters are considered in this study: Surface air temperature and solar irradiance as observed and projected data. The observed data used are grouped into three sets of data. The first set of data is a collection of observed data from the synoptic stations of Burundi belongs to the Geography Institute of Burundi (IGEBU). Monthly and annual mean temperature data for seven stations located in the study area are used over the historical period of thirty years . These climate stations are presented in Figure 1 and their geographic information and description are given in Table 1. However, some stations had uncompleted data and shorter duration than the 30 years recommended. To complete observed data, the second set of data used is MERRA-2data set available at http://gmao.gsfc.nasa.gov/reanalysis/MERRA-2 [28] and provided by NASA/Goddard Space Flight Center.
Due to the lack of solar irradiance data measured in the stations of our study area, the third set of observed data is provided by the SoDa database available at http://www.soda-is.com [29]. Furthermore, due to the lack of long-term solar irradiance measurement, data from 1985-2005 are used as an observational solar irradiance dataset with spatial resolution of 20 km. The SoDa server provides hourly, daily, monthly and annually irradiation over Europe, Africa, and Atlantic Ocean which is accessible on free basis [30]. The SoDa irradiance is satellite-derived data of HelioClim-1 Daily solar irradiance v4.0. HelioClim-1 databases use the Heliosat-2 method to process the Meteosat images. The Heliosat-2 method converts images acquired by meteorological geostationary satellites such as Meteosat (Europe), GOES (USA) or GMS (Japan), into data and maps of solar radiation received at ground level [30]. Table 2   Climate projections data used are daily surface air temperature and solar irradiance taken from eight regional climate models available in the context of Coordinated Regional Climate Downscaling Experiment (CORDEX) over Africa at 0.44 • resolution for the period 1950-2100 [31] and accessed online https://www.cordex.org [32]. Table 3 presents used climatic models, home institute, and model short name adopted in this paper. Climatic models are complex programs based on atmospheric circulation including its chemistry and radiation, oceanic circulation including its biochemistry, land-surface, river routing and sea ice modeling [6]. The differences between them are mainly related to the physical parameterization of each component of the model structure. IPCC has established four Representative Concentration Pathway (RCP) scenarios which are linked to the concentration of greenhouse gas emission during the 21st century. According to the socio-economic technological and policy development activities that disturb the concentration of carbon dioxide in the atmosphere, there are RCP2.6, RCP4.5, RCP6.0 and RCP8.5 where the associated number corresponds to the radiative forcing reached at the end of 21st century [33].
In this study, experiments were performed following RCP8.5, which is the highest forcing radiative scenario, are considered for two future periods (named near future and far future) in CORDEX database. Due to the short baseline period of solar irradiance the period 2026-2045 for near future and 2066-2085 for far future have been considered for these two parameters. The RCP 8.5 is based on the A2r scenario which combines assumptions about high population and relatively low income growth with modest rates of technological change and energy intensity improvements, leading in the long term to high energy demand and greenhouse gas emissions in the absence of climate change policies [34]. This set of data is the up-to-date ensemble of high-resolution Regional Climate Model (RCM) projections, and data have been bias corrected using historical observed data.

Bias Correction
RCM temperature and solar irradiance were corrected for their biases in mean and standard deviation for each calendar month, following the methodology presented in Haerter et al. [35]. The bias in the mean is corrected by subtracting the differences found between observed and modeled values and a correction to the model data is performed to conform to the variability of the historical data. This procedure takes the sequence of anomalies and scales them consistently with the observed historical variability. As used by Ioanna S. Panagea et al. [30], for data following a normal distribution, the transfer function is linear and is of the form shown in the following equation: where µ cor sc is the final adjusted time series, µ sc mod is the raw model predictions for the scenario period µ con obs and µ con mod are the mean of observed and modeled data for the control period, respectively, and σ con obs and σ con mod are the standard deviations of observed and modeled data for the control period, respectively. The final adjusted model time series exhibits the appropriate baseline mean and standard deviation with respect to the observed data [30].

Variability Analysis
The discrimination of high or low values of each climate parameter considered was done by computing the standardized variable index (I) created by Reference [36] and which can be found in many articles [37,38]. Then, the standardized variable index (I) is defined by Equation (2): where χ i , χ m and σ are respectively the value for the considered year (or month i), the average and the standard deviation of the time series. Thus, in this work we considered that a year is normal when its index is included between −0.5 and +0.5. It will be in surplus if its index is greater than +0.5, and will be considered as deficit when its index is below than −0.5. This interval is relatively weak and is criticized but it allows distinguishing of the deficit years from surplus years. The standardized variable index makes also possible to analyze the interannual variability of the variables at considered time scale.

Future Climate Analysis
The temporal assessment of the future climate has been done using Regional climate models (RCM) over two periods named near future and far future. The period 2026-2045 and 2066-2085 for surface air temperature and surface solar irradiance have been considered with respect of the reference period 1986-2005. Following RCP8.5 scenario, eight RCM (CCCma, CNRM, ICHEC, IPSL, MIROC, MPI-M, NCC and NOAA) are used to plot projected trends at annual scale and a linear trend has been plotted to show the multi model's trend, after bias correction according to observational data set. Furthermore, a multi-model mean has been computed and adopted at a monthly scale.
The analysis of climate changes was done in two steps according to two future periods considered in this study with reference to the baseline period. We compute and add delta change factors to monthly time series of observation records. The calculations of delta change factors were based on a selected projection interval (2026−2045) with a baseline period (1986-2005) for both temperature and solar irradiance using the Equations (3). The same computation was repeated with the selected projection interval (2066−2085). Then the change factors for temperature and solar irradiance were determined by: where ∆µ is the change factor ( • Cor w/m 2 ), µ is the monthly mean ( • Cor w/m 2 ), i is the RCM, j is the projection period, ref is the reference period and k is the month. Two statistics tests, Mann Kendall and t-test have been used respectively to clarify the trends detected of a time series and to test the mean equality of two times series for each variable. The Mann Kendall's test (MK) [39,40], which is used by many authors, was used for detecting trends in time series. MK checks increasing, decreasing, or no trends of the projected models mean time series of the three regions of our study area. If there is a linear trend in the time series, the true slope has been estimated using Sen's slope method [21] which is a simple non-parametric found in many articles [41,42]. Furthermore, the t-test [43] has been used to prove whether the difference between two ensemble models means over the baseline and future period is statistically significant for a considered region. Before using the t-test, which is a parametric test, normal distribution of the used data must be assessed.  The analysis of Figure 2a shows that June and July present low mean values of temperature, whereas September and October show high mean values but the pattern does not show remarkable variations. This figure reveals also that ELL has high values and EAP has low values of surface air temperature over the year.

Variability Analysis in the Historical Period
The means monthly temperature for NLL, EAP and ELL used to calculate the standardized variable index are respectively 20.56 °C, 19.38 °C and 21.75 °C. In the same way, the standard deviations are 0.54, 0.61, and 0.92 respectively for NLL, EAP and ELL. The analysis of Figure 2b shows that six months (April, May, June, July, November and December) are cooler at NLL while three months (May, June and July) are cooler at EAP and ELL. The coolest month with the lowest index value equal to −1.66 is July for ELL while the coolest month is June for NLL and EAP with the lowest index values equal to −1.27 and −1.51, respectively. The analysis shows also that four months (February, August, September and October) are hotter at NLL, while three months (August, September and October) are the warmest at EAP and only two months are the hottest (September and October) at ELL. Normal months in the year are January, February, March, November and December for all regions except NLL where normal months are January and March only. The hottest month is September for both regions with the highest index values equal to 1.87, 1.89, and 1.74 respectively for NLL, EAP, and ELL. Figure 3a presents the interannual temperature patterns for the period 1976-2005 at NLL, EAP and ELL. This figure shows that the ELL is the hottest region with about mean value of 21.7 °C whereas EAP is the coolest region of the study area with the 19.4 °C of annual mean value. This high temperature over low land is due to the low altitude and as we can remark that there is an upward trend with the increase changes. The analysis of Figure 2a shows that June and July present low mean values of temperature, whereas September and October show high mean values but the pattern does not show remarkable variations. This figure reveals also that ELL has high values and EAP has low values of surface air temperature over the year.

Interannual Temperature
The means monthly temperature for NLL, EAP and ELL used to calculate the standardized variable index are respectively 20.56 • C, 19.38 • C and 21.75 • C. In the same way, the standard deviations are 0.54, 0.61, and 0.92 respectively for NLL, EAP and ELL. The analysis of Figure 2b shows that six months (April, May, June, July, November and December) are cooler at NLL while three months (May, June and July) are cooler at EAP and ELL. The coolest month with the lowest index value equal to −1.66 is July for ELL while the coolest month is June for NLL and EAP with the lowest index values equal to −1.27 and −1.51, respectively. The analysis shows also that four months (February, August, September and October) are hotter at NLL, while three months (August, September and October) are the warmest at EAP and only two months are the hottest (September and October) at ELL. Normal months in the year are January, February, March, November and December for all regions except NLL where normal months are January and March only. The hottest month is September for both regions with the highest index values equal to 1.87, 1.89, and 1.74 respectively for NLL, EAP, and ELL. Figure 3a presents the interannual temperature patterns for the period 1976-2005 at NLL, EAP and ELL. This figure shows that the ELL is the hottest region with about mean value of 21.7 • C whereas EAP is the coolest region of the study area with the 19.4 • C of annual mean value. This high temperature over low land is due to the low altitude and as we can remark that there is an upward trend with the increase changes.   The analysis of the figure shows also 11 years of high temperature for NLL and EAP and 9 years of high temperature for ELL while 10 years of low temperature are observed at NLL and, 9 years of low temperature at EAP and ELL.

Interannual Temperature
These results reveal an upward trend of average temperature for the last sub-period1992-2005 and pointed out the hottest years occurred in 1999 for NLL and ELL and 1998 for EAP whereas the coolest years occurred in 1986 for EAP and ELL with 1978 for NLL. Furthermore, these findings agree with the results of many studies which show increases in temperature over the end of the last century, See References [19,22,44]. More recently, the study of Birara et al. [21] analyzed the trend and variability of Rainfall and temperature in the Tana basin region, Ethiopia, for the period covers 1980-2015 and found the significant increases of maximum, minimum and mean temperature for most of the stations considered and an increase of 1.08 °C was observed. These increases in temperature may affect negatively the production of solar energy and increase evaporation which reduces water availability.   The analysis of the figure shows also 11 years of high temperature for NLL and EAP and 9 years of high temperature for ELL while 10 years of low temperature are observed at NLL and, 9 years of low temperature at EAP and ELL.

Seasonal Solar Irradiance
These results reveal an upward trend of average temperature for the last sub-period1992-2005 and pointed out the hottest years occurred in 1999 for NLL and ELL and 1998 for EAP whereas the coolest years occurred in 1986 for EAP and ELL with 1978 for NLL. Furthermore, these findings agree with the results of many studies which show increases in temperature over the end of the last century, See References [19,22,44]. More recently, the study of Birara et al. [21] analyzed the trend and variability of Rainfall and temperature in the Tana basin region, Ethiopia, for the period covers 1980-2015 and found the significant increases of maximum, minimum and mean temperature for most of the stations considered and an increase of 1.08 • C was observed. These increases in temperature may affect negatively the production of solar energy and increase evaporation which reduces water availability.   Figure 4b presents the monthly standardized solar irradiance at NLL, EAP and ELL over the historical period. The analysis of this figure shows that four months (April, May, November and December)are in deficit of solar irradiance at NLL and four months also (January, May, November and December) are in deficit of solar irradiance at EAP and ELL. The month with the lowest index value equal to −1.4 is May for NLL while the month with low value of solar irradiance is December for ELL and EAP with the lowest index value equal to −1.26 and −1.13, respectively. The analysis shows also that four months (February, July, August and September) are in excess of solar irradiance at NLL while four months also(June, July, August and September) for ELL with three months (July, August, September) at EAP are in excess of solar irradiance. Normal months are February, March, April and October for EAP and ELL whereas January, March, June and October are normal for the region of NLL. The Month with highest index value of solar irradiance is July for all regions with 1.89, 1.95, and 1.65 for NLL, EAP and ELL, respectively.
The monthly means solar irradiance for NLL, EAP and ELL used to calculate the standardized variable index are respectively 248 W/m², 236 W/m², and 231 W/m², whereas standard deviationsare 10.9, 13.2, and 19.5 respectively for NLL, EAP, and ELL. This analysis shows that the months with highest values of solar irradiance coincide with the dry season, which leads to the conclusion that there is an excess of solar irradiance in the dry season and, therefore, it must be a period of high production of solar energy. Indeed, a solar power plant installed in this region could be in complementary with a hydroelectric plant with low production in the dry season because of the reduction of the water level in the reservoir.

Interannual solar irradiance
The Figure 5a gives the interannual solar irradiance patterns for the period 1985-2005 at NLL, EAP and ELL. This figure shows that NLL indicates high values of solar irradiance than other regions which have close values over the full period. In the other hand, this pattern does not show long variability of solar irradiance values through the historical period. We can also observe that EAP and ELL have likely the same values and there are not many differences, and some months show a mixture of solar irradiance values.  The monthly means solar irradiance for NLL, EAP and ELL used to calculate the standardized variable index are respectively 248 W/m 2 , 236 W/m 2 , and 231 W/m 2 , whereas standard deviationsare 10.9, 13.2, and 19.5 respectively for NLL, EAP, and ELL. This analysis shows that the months with highest values of solar irradiance coincide with the dry season, which leads to the conclusion that there is an excess of solar irradiance in the dry season and, therefore, it must be a period of high production of solar energy. Indeed, a solar power plant installed in this region could be in complementary with a hydroelectric plant with low production in the dry season because of the reduction of the water level in the reservoir.

Interannual Solar Irradiance
The Figure 5a    This analysis shows also a general upward of average solar irradiance for the first sub-period. Furthermore, these findings agree with the results of many studies which revealed a subsequent increase (1-4 W/m 2 ) since the 1980s called "brightening" [9] and decreasing of irradiation in some regions of the world. Solar irradiation is the principal resource of solar power output (PV); these decreases of this parameter affect negatively the production of solar energy.  This analysis shows also a general upward of average solar irradiance for the first sub-period. Furthermore, these findings agree with the results of many studies which revealed a subsequent increase (1-4 W/m 2 ) since the 1980s called "brightening" [9] and decreasing of irradiation in some regions of the world. Solar irradiation is the principal resource of solar power output (PV); these decreases of this parameter affect negatively the production of solar energy.

Monthly Scale Changes
Temperature Figure 6a presents the projected monthly surface air temperature over the period 2026-2045 and 2066-2085 according to reference period 1986-2005. The data used for future projected temperature are the mean values of eight RCMs after bias correction and the reference period were presented by observed data from stations and mean values of raw RCMs. For the near future, the analysis of Figure 6a projects an increase in surface air temperature from January to September but also projects a small decrease from October to December for all regions according to the reference period. The highest increase is projected to be in July which reaches the maximum of 3.88 • C at ELL. On the other hand, the highest decrease is projected in November and reaches the minimum of −0.93 • C at ELL also. The hottest month is projected to be September for NLL and ELL, except for EAP, which is projected to be August, and the highest value will reach 24.79 • C at ELL. Then, the coolest month is forecasted to be November for all regions, with the lowest value of 18.84 • C at EAP. The far future period 2066-2085 projects greatest increases of surface air temperature compared to the near future as observed in Figure 6a. The highest increase is projected to be in July also with a maximum change value of 4.85 • C and the lowest increase is forecasted to be in November with a change value of 0.06 • C. For this period, the hottest month is projected to be September for all regions where the highest value is projected to be 26.00 • C at ELL. Otherwise the coolest month is forecasted to be November for all regions and lowest value is forecasted to be 19.54 • C at EAP. The findings of this study are closely in accordance with many studies carried out before locally or at regional scale. Reference [23] assessed the trends and changes detection in rainfall, temperature and wind speed in Burundi using seven RCMs up the period 2021-2050 and found that all models predict increases in temperature and wind speeds compared to the baseline period, 1981-2010. Solar Irradiance Figure 6b indicates the projected monthly solar irradiance over the period 2026-2045 and 2066-2085 according to the baseline period 1986-2005 at NLL, EAP, and ELL. The data used for future projected solar irradiance are the mean values of eight RCMs used after bias correction, whereas observed data from Soda and raw CORDEX data are used to present the historical period. For the first period (2026-2045), the analysis of this figure distinguishes two parts of the pattern. The first one shows an increase in solar irradiance from January to June with the maximum change occurred in March for all regions. The highest percentage changes values are 7.97%, 8.76%, and 13.53%, respectively, at NLL, EAP, and ELL. The second part shows a decrease in solar irradiance from July to December with a maximum change occurred in July for NLL and EAP, and November for ELL. The lowest percentage changes values are −7.96% and −7.71% at NLL and EAP respectively and −13.85% at ELL. Moreover, the highest value of solar irradiance is projected to be in February for all regions with 270.63 w/m 2 , 258.27 w/m 2 and 262.97 w/m 2 at NLL, EAP and ELL, respectively. The lowest value of solar irradiance is projected to be in November with 221.37 w/m 2 , 203.58 w/m 2 and 183.45 w/m 2 at NLL, EAP and ELL, respectively. We can note that the highest monthly changes (increases or decreases) in solar irradiance are observed in ELL.
In the same way, for the far future period the figure gives the projected monthly solar irradiance over the period 2066-2085 according to the baseline period 1986-2005 at NLL, EAP, and ELL. The analysis of this figure also shows a mixture of increasing and decreasing over the pattern: The increases in solar irradiance occurred from January to June and September, with the maximum change found in March for all regions.
The highest percentage changes values are 8.27%, 9.53%, and 15.03% respectively at NLL, EAP and ELL. The decreases in solar irradiance appeared in July and August and from October to December with a maximum decrease change occurred in November for NLL, EAP, and ELL. The highest percentages of decreasing changes are −12.05%, −12.25%, and −17.29% at NLL, EAP and at ELL respectively. Moreover, the highest value of solar irradiance is projected to be in March for all regions with 264.8 w/m 2 , 254.1 w/m 2 and 258.4 w/m 2 at NLL, EAP, and ELL respectively. The lowest value of solar irradiance is projected to be in November with 211.5 w/m 2 , 194.75w/m 2 and 176.15 w/m 2 at NLL, EAP, and ELL respectively.

Temperature
The spatial distribution of observed, models raw mean values of temperature for the reference period and projected RCMs annual mean values of temperature after bias correction over the period 2026-2045 and 2066-2085 is given by Figure 7 following rcp 8.5 scenario. The annual mean temperature is expected to increase in the future all over the regions studied. The analysis of this figure shows that this increasing change is very marked in the far future rather than the near future. Furthermore, the region ELL which is experiencing high temperature is expected to have higher increase changes than other regions, whereas the low increase changes are expected at EAP for both future periods. The projected increase in temperature found in this paper is consistent with the findings by Brohan et al. [19] and may be attributed to global warming. These features will affect the lifestyle of the population in these regions and may present a negative effect on solar energy production because an increase of temperature affects negatively the PV power output. Compared with the findings shown in Table 7, these results predict low decrease of solar irradiance than far future as said before. In fact, the results obtained in period 2066-2085 presented in the Table 7 show projected changes of −7.6%, −6.0%, and −7.4% by the ensemble models mean at EAP, ELL, and NLL. Indeed, the highest change is forecasted by CNRM at EAP and the value is −10.9%. Even if NLL presents high values of solar irradiance, ELL and NLL are expected to have high deltas changes for almost all models and for ensemble models mean. The analysis of these two tables reveals that surface solar irradiance will experience small annual changes in the future which determine an interest of installing solar power plan in the regions. Despite the slight decreasing of solar irradiance predicted over all regions, the ensemble models mean differences between baseline and near future or far future are significant over all regions regarding t-test's results with a significant level α = 0.05. In fact, the p-value is very lower than α and the observed value |t| is greater than the critical value. Before using the parametric t-test, a normal distribution of solar irradiance data is evaluated for all periods and over all regions.
For trend detection over the period 2026-2045, the MK's test indicates no significant trend in forecasted surface solar irradiance at significance level α = 0.05 following rcp 8.5 scenario for all the regions. The same conclusion is observed over the far future for all the regions. In fact, the p-values found of all time series are greater than α = 0.05 for the far future and near future,so the null hypothesis which assumes a trend can be rejected.Therefore, the Sen's slope showed a downward trend in the ensemble model mean with small values of −0.012, −0.03 and −0.114 respectively at EAP, ELL and NLL over the near future period. In the same way, a downward trend is also detected by the Sen's slope with small values of −0.089, −0.077 and −0.167 respectively at EAP, ELL, and NLL over far future. Even if the MK's test is not significant, the trends of solar irradiance predict a small decreasing feature for a long time. For projected changes in interannual average temperature, Tables 4 and 5 show changes in projected average temperature at annual scale for the two periods considered. Used data are each model raw mean value for the reference period and projected RCMs annual mean values of temperature. The analysis of the Table 4 shows that at EAP, ELL, and NLL will experience a change of 0.4 • C, 3.6 • C, and 2.5 • C by 2045 in the ensemble models-mean according to the baseline period 1986-2005. In the same way, the analysis of the Table 5 shows that a change of 2.1 • C and 6.9 • C and 5.9 • C are respectively expected in the ensemble models mean at EAP, ELL and NLL by 2085 according to the reference period 1986-2005.The analysis of the two tables shows that higher changes are normally expected in the far future rather than the near future. We can note that there are few decrease changes in the near future, especially at EAP region predicted by CCCma and IPSL. Indeed, the highest changes occurred at ELL and are predicted by ICHEC in both future periods. In general, ELL presents high deltas for almost all models over all regions and for ensemble models mean. Therefore, the lowest change occurred at EAP with IPSL in the near future and with CCCma by 2085. Furthermore, with a significant level of 0.05, the ensemble models mean differences between the baseline period and each of two future periods are statistically significant over all regions according to the t-test with a p-value less than 0.0001. However, an ensemble models mean difference between 1986-2005 and 2026-2045 at EAP presents a slight significance compared to others regions, as well as observed value |t| is very close to critical value. Given that t-test is a parametric test, a normal distribution of all temperature data used has been controlled.  For trend detection in projected average temperature at significance level α = 0.05, the Mann Kendall's test(MK) revealed upward trends all over the regions for both periods 2026-2045 and 2066-2085 with a very small p-value = 0.0001. Indeed, the Kendall's tau for both periods is likely the same. In fact, over the period 2026-2045 the Kendall's tau equals to 0.81 for EAP, 0.75 for ELL and it equals 0.82 for NLL, while in period 2066-2085 the Kendall's tau is 0.77 for EAP and 0.80 for ELL and NLL. Combined with the Sen's slope method, these statistical tests revealed upward trends of surface air temperature in the considered periods all over the regions of the study area. Solar Irradiance Figure 8 presents a spatial distribution of annual mean of solar irradiance data from Soda database (observed), annual mean from original CORDEX in the reference period 1986-2005 and projected models mean of solar irradiance after bias correction over two future periods namely near future 2026-2045 and far future 2066-2085. Following rcp 8.5, the entire region does not show the highest change either in the near future or far future. However, the analysis of the figures shows that NLL is experiencing high values of solar irradiance whereas ELL presents low values. In addition, a slight downward trend of solar irradiation is also observed generally in far future over all regions. Furthermore, the analysis of this figure shows also that the highest changes are observed at EAP for both future periods.  The projected changes in solar irradiance for the period 2026-2045 are summarized in Table 6, while Table 7 presents the projected changes for far future 2066-2085 both using each raw model mean value over the reference period and projected RCMs annual mean values of solar irradiance. The two tables show that the projections are very close; therefore, the last table presents a higher decrease of solar irradiance in far future than near future.  The analysis of Table 6 shows that the multi models mean project a decrease change of −6.3%, −5.4%, and −5.8% at EAP, ELL, and NLL, respectively. Indeed, the highest decreasing change is projected by NCC at EAP.
Compared with the findings shown in Table 7, these results predict low decrease of solar irradiance than far future as said before. In fact, the results obtained in period 2066-2085 presented in the Table 7 show projected changes of −7.6%, −6.0%, and −7.4% by the ensemble models mean at EAP, ELL, and NLL. Indeed, the highest change is forecasted by CNRM at EAP and the value is −10.9%. Even if NLL presents high values of solar irradiance, ELL and NLL are expected to have high deltas changes for almost all models and for ensemble models mean. The analysis of these two tables reveals that surface solar irradiance will experience small annual changes in the future which determine an interest of installing solar power plan in the regions. Despite the slight decreasing of solar irradiance predicted over all regions, the ensemble models mean differences between baseline and near future or far future are significant over all regions regarding t-test's results with a significant level α = 0.05. In fact, the p-value is very lower than α and the observed value |t| is greater than the critical value. Before using the parametric t-test, a normal distribution of solar irradiance data is evaluated for all periods and over all regions.
For trend detection over the period 2026-2045, the MK's test indicates no significant trend in forecasted surface solar irradiance at significance level α = 0.05 following rcp 8.5 scenario for all the regions. The same conclusion is observed over the far future for all the regions. In fact, the p-values found of all time series are greater than α = 0.05 for the far future and near future, so the null hypothesis which assumes a trend can be rejected. Therefore, the Sen's slope showed a downward trend in the ensemble model mean with small values of −0.012, −0.03 and −0.114 respectively at EAP, ELL and NLL over the near future period. In the same way, a downward trend is also detected by the Sen's slope with small values of −0.089, −0.077 and −0.167 respectively at EAP, ELL, and NLL over far future. Even if the MK's test is not significant, the trends of solar irradiance predict a small decreasing feature for a long time.

Conclusions
The analysis of variability and projected trends of temperature and solar irradiance has been performed considering historical period and two future periods. Observed temperature data from meteorological stations and solar irradiance data from the Soda database were used over the historical period. In addition, projection data from eight Regional Climate Models were used over the periods 2026-2045 and 2066-2085 after bias correction for their mean and standard deviation in order to produce correctly the future climate and conform to the observed data. The characterization of the historical period was done using the standardized variable index for each climatic parameter. The interannual results revealed an upward trend in average temperature for the last sub-period 1992-2005. At the monthly scale, the analysis showed seasonal variability where an upward trend in average temperature appears from August-October, and the hottest month is September for all regions. A downward trend occurred from May-July, and the coolest month is July. For solar irradiance, the monthly analysis showed that the surplus months coincide with the dry season, which leads to the conclusion that there is an excess of solar irradiance in the dry season with the highest value appearing in July. The interannual results revealed a downward trend in average solar irradiance for the last sub-period 1999-2005 and showed also a general upward trend in average solar irradiance for the first sub-period.
For future analysis, at a monthly level, an increase in surface air temperature from January-September and a small decrease from October-December are projected for all regions according to the reference period. Therefore, the far future presents the same features as the near future with a high amplitude of increase over the far future. The highest increase is projected to reach a maximum change value of 4.85 • C. The hottest month is projected to be September for all regions where the highest value is projected to be 26.00 • C at ELL. For solar irradiance, the seasonal analysis distinguishes two parts of the pattern for both periods. The first one shows an increase in solar irradiance from January-June with the maximum change occurring in March for all regions. The second part shows a decrease in solar irradiance from July-December with a maximum change occurring in July for NLL and EAP, and November for ELL. The analysis of the results revealed that the far future presents the same features as the near future, where the changes are very marked over the far future.
At the interannual scale, considering air ambient temperature, the results show that lowlands are expected to be hotter than arid plateaus. For trend detection in projected average temperature at significance level α = 0.05, the MK's test revealed upward trends all over the regions for both the future periods. The ensemble models mean differences between the baseline period and each of the two future periods are statistically significant over all regions according to t-test. Considering projected changes, the analysis shows that at EAP, ELL, and NLL will experience a change of 0.4 • C, 3.6 • C, and 2.6 • C by 2045 in the ensemble models mean according to the baseline period 1986-2005, whereas changes of 2.1 • C and 6.9 • C and 5.9 • C are respectively expected in the ensemble models mean at EAP, ELL, and NLL by 2085. The analyses show also that higher changes are normally expected in the far future than the near future.
Interannual analysis of solar irradiance shows clearly a slight downward trend of solar irradiance for both periods. The results predict a downward average of solar irradiance over all regions studied and for all models used. For trend detection over both periods, the MK's test indicates no significant trend in forecasted solar irradiance at a significance level α = 0.05 over allregions. Therefore, the Sen's slope showed a downward trend in the ensemble models mean with small values over the near future and far future periods. The ensemble models mean differences between baseline and near future or far future are significant over all regions regarding t-test's results. The projected changes in solar irradiance show that the projections are very close over two periods. The multi models mean projected the decreasing changes of −6.3%, 5.4%, and 5.8% at EAP, ELL and NLL respectively over 2026-2045. The results obtained in the period 2066-2085 show the projected changes of −7.3%, 6.0%, and −7.4% by ensemble models mean at EAP, ELL and NLL. The findings show that projected downward trends in solar irradiance over the whole studied region are not significant at a considered threshold of 5%, while upward trends in air temperature are very significant.