Next Article in Journal
Homogeneity Assessment and Correction Methodology for the 1980–2022 Daily Temperature Series in Padua, Italy
Next Article in Special Issue
Safe Sowing Windows for Smallholder Farmers in West Africa in the Context of Climate Variability
Previous Article in Journal
Multi-Hazard Extreme Scenario Quantification Using Intensity, Duration, and Return Period Characteristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linking Climate Change Information with Crop Growing Seasons in the Northwest Ethiopian Highlands

by
Gashaw Bimrew Tarekegn
1,
Addis A. Alaminie
2 and
Sisay E. Debele
3,*
1
Institute of Disaster Risk Management and Food security Studies, Bahir Dar University, Bahir Dar P.O. Box 5501, Ethiopia
2
Faculty of Civil and Water Resources Engineering, Bahir Dar University, Bahir Dar P.O. Box 400, Ethiopia
3
Global Centre for Clean Air Research (GCARE), School of Sustainability, Civil and Environmental Engineering, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7XH, UK
*
Author to whom correspondence should be addressed.
Climate 2023, 11(12), 243; https://doi.org/10.3390/cli11120243
Submission received: 16 October 2023 / Revised: 30 November 2023 / Accepted: 6 December 2023 / Published: 15 December 2023

Abstract

:
In Ethiopia, the impacts of climate change are expected to have significant consequences for agriculture and food security. This study investigates both past (1981–2010) and future (2041–2070) climate trends and their influence on the length of the growing season (LGS) in the North-Western Ethiopian highlands. Climate observations were obtained from the National Meteorological Agency of Ethiopia, while the best performing and highest resolution models from the CMIP5 experiment and RCP6 (Coupled Models Intercomparison Project and representative concentration pathway 6) were used for the analysis. Standard statistical methods were applied to compute soil water content, evaluate climate variability and trends, and assess their impact on the length of the growing season. Maximum temperature (tasmax) and minimum temperature (tasmin) inter-annual variability anomalies show that the region has experienced cooler years than hotter years in the past. However, in the future, the coolest years are expected to decrease by −1.2 °C, while the hottest years will increase by +1.3 °C. During the major rainfall season (JJAS), the area has received an adequate amount of rainfall in the past and is very likely to receive similar rainfall in the future. On the other hand, the rainfall amount in the season February to May (FMAM) is expected to assist only with early planting. For the season October to January (ONDJ), the rainfall amount may help lengthen the growing season of JJAS if properly utilized; otherwise, the season has the potential to destroy crops before and during the harvesting time. The soil water content changes in the future remain close to those of the past period. The length of growing seasons has less variable onset and cessation dates, while in the future, the length of the growing period (LGP) of 174 to 177 days will be suitable for short- and long-cycle crops, as well as double cropping, benefiting crop production yield in the North-Western Ethiopian highlands in the future.

1. Introduction

There is a consensus that the impacts of climate change on agriculture will significantly contribute to the development challenges of ensuring food security and reducing poverty, particularly in Africa. Climate change plays a crucial role in agricultural production, directly impacting processes from the start of land preparation to the final harvest [1,2,3,4]. In the past half-century, Ethiopia has experienced a late onset and early cessation of seasonal rains, leading to the associated failure of the crop growing season [5,6]. To optimally utilize seasonal rainfall for agricultural production, additional knowledge and understanding are needed regarding how shifts in rainfall seasons may affect crop yields. Addressing this issue requires tailored predictive information on the start and end of the growing season, which is timely and important and can assist in adopting appropriate adaptation and coping measures.
Climate change has various impacts on the growth and development processes of crops. For instance, an increase in carbon dioxide (CO2) can stimulate photosynthesis rates, sometimes resulting in higher yields [7,8]. Changes in temperature and precipitation affect crop photosynthesis, plant development rates, as well as water and nutrient budgets in the field [9,10]. One of the critical processes affected by climate change is photosynthesis. The direct effects of CO2 enrichment on plants include an increase in the rate of photosynthesis and water-use efficiency (WUE). As CO2 concentrations increase, the transpiration intensity of plants reduces by partially closing the stomata, leading to improved water use efficiency and a lower probability of water stress occurrence. These physiological responses are known as the CO2 fertilization effect or the direct effect of increased CO2.
The increase in temperature due to climate change has both positive and negative impacts on crop production. For example, in the middle and higher latitudes and at high elevations in the tropics, global warming will extend the potential growing season’s length, allowing for earlier planting, maturation, and harvesting. Increased CO2 may enable the completion of two or more cropping cycles during the same season and the expansion of crop-producing areas to higher latitudes and elevations, provided soil fertility is maintained [7,11,12]. In warmer, lower latitude regions, increased temperatures may accelerate the release of CO2 through plant respiration, creating less-than-optimal conditions for net growth. High temperatures reduce yield by hastening maturation, not allowing crops to progress slowly through the season for maximum resource capture and assimilated partitioning to reproductive structures [9,13].
Therefore, under conditions that are physiologically ‘too warm’ for certain crops, yields may decrease. An increase in average temperature can have several effects, including lengthening the growing season in regions with a relatively cool spring and fall, and adversely affecting crops in regions where summer heat already limits production. Soil evaporation rates may increase, and the chances of severe droughts may also rise [14,15]. Warmer temperatures in high-altitude areas may allow more insects to overwinter in these areas. Crop damage from plant diseases is likely to increase in temperate regions because many fungal and bacterial diseases have a greater potential to reach severe levels when temperatures are warmer, or precipitation increases [16].
Changes in rainfall can affect soil erosion rates and soil moisture, both of which are crucial for growth and yields. In the case of the NW Ethiopian highlands, the impact of climate change is likely to be both positive and negative due to the area’s different topography with high-, middle-, and low-altitude elevations.
Moreover, improved forecasts of seasonal precipitation with adequate spatial resolution could potentially increase agricultural production and reduce production risks [17]. Therefore, predictions of local determinants of climate change, climate variability, and the growing season are crucial for livelihoods dependent on rain-fed agriculture. Hence, this study uses fine-resolution climate change models to simulate and project climatic conditions and associated growing seasons in the NW Ethiopian highlands.
According to the World Bank’s [18] categorization, the NW Ethiopian highlands consist of majority areas with substantially higher rainfall and more availability of arable land, where the average crop yield is double that of the whole country. Despite being potential zones for crops and livestock, the NW Ethiopian highlands experience erratic rainfall in terms of both spatial and temporal distribution, with dry spells significantly reducing crop yields and sometimes leading to total crop failure. This has caused the highlands to face recurrent drought [18].
The Ethiopian highlands, major contributors to crop production for the country, were affected by inter-annual and seasonal variability of rainfall, causing fluctuations in cereal production in the region [19,20]. Hence, research on climate change’s impact on crop production and agricultural water management is timely and important to enhance food security in the country. Although some studies have been conducted to identify the spatial and temporal climate changes of the past and future, these have not specifically been tied to the growing season at the local level [3,4,21,22].
Moreover, the impact analysis of climate change on soil water content, and the onset and cessation date and length of the growing period derived from soil water content, provide more reliable results than other methods used to determine crop growing seasons. Hence, this research analyzes climate change and its impact on the length of the growing season (based on soil water content) in the Ethiopian highlands. The research results can be used to design suitable adaptation and mitigation strategies, as well as coping mechanisms. Thus, in this study, we address the past (1981–2010) and future (2041–2070) CMIP5 RCP6 model projections of trends and changes in specific climatic parameters (rainfall, maximum and minimum temperatures, PET, and soil water content) and their impact on the length of crop growing seasons in the NW Ethiopian highlands.

2. Materials and Method

2.1. Study Area

Figure 1 shows the location of NW Ethiopia and the distribution of the 35 meteorological stations utilized in the study. The NW Ethiopian highlands are situated in the north-western part of Ethiopia, located within 8–14° N and 35–39° E (Figure 1).

Agriculture of the Study Area

Agriculture has been practiced for centuries in Ethiopia, with the overall farming system strongly oriented towards grain production as a source of livelihood and a way of life. Ethiopian peasant agriculture is characterized by polyculture, involving the traditional rotation of crops and livestock. This agricultural approach has persisted for centuries and continues to be beneficial for both the land and the people, meeting complete livelihood and nutritional needs, as well as ensuring community security. Despite remaining subsistent, peasant agriculture allows individuals to produce everything required for themselves and their families, with a variety of crops being cultivated.
The proximity of areas with different altitudes has contributed to the diversity of crops grown. Farmers primarily focus on food crops, with only a few cash crops intended for the local market or barter system. However, Ethiopia has experienced droughts in the recent half-century [6,23], mainly due to the late onset and early cessation of seasonal rains, leading to the failure of the crop growing season. This is particularly challenging as the agricultural practices of most farmers in the country rely on rain-fed agriculture. To optimally utilize seasonal rainfall for agricultural production, additional knowledge is needed on how shifts in rainfall seasons may affect crop yields.
On the other hand, the soil type in the Ethiopian highlands is predominantly Vertisols, accounting for 12.7 million hectares in Ethiopia, with 7.6 million hectares in the NW Ethiopian highlands. Waterlogging of Vertisols is more severe in the Ethiopian highlands, where rainfall is higher and evaporative demands are lower. The technical efficiency of traditionally applied surface drainage techniques is not sufficient to fully utilize the potential of these soils.

2.2. Data Sets and Research Design

The data sources utilized for analysis included gauge data from weather stations, gridded observations (reanalysis) data, and CMIP5 rcp6 (representative concentration path 6) model simulation datasets extracted from the KNMI-Climate Explorer website (https://climexp.knmi.nl/start.cgi) (accessed on 23 May 2023). The Royal Netherlands Meteorological Institute (KNMI) hosts the Climate Explorer website, providing access to CMIP5 data. KNMI collaborates closely with international meteorological networks, such as the European Centre for Medium-Range Weather Forecasts (ECMWF) and the European meteorological satellite network EUMETSAT. The KNMI Climate Explorer (CE) web page, organized by [24], is a valuable resource. For the list of CMIP5 model groups open to the public on KNMI CE, see Supplementary Information (SI) Table S1.
This study focused on using CMIP5 rcp6 model data sources developed by the AIM modeling team at the National Institute for Environmental Studies (NIES) in Japan. The decision to use only CMIP5 rcp6 model groups is explained as follows. (1) CMIP5 rcp6 is a stabilization scenario where total radiative forcing stabilizes shortly after 2100, without overshot, through the application of various technologies and strategies to reduce greenhouse gas emissions [25,26]. This aligns with the global consensus on reducing carbon emissions and mitigating climate change. (2) CMIP5 rcp6 models consist of 23 to 25 models serving as data sources for 21 climate variables. All these models’ outputs were considered with equal weight given to each model or version with different parameterizations. Describing each CMIP5 rcp6 model is beyond the scope of this study; readers are referred to [27] and the Program for Climate Model Diagnosis & Intercomparison (https://shorturl.at/zHLPV) (accessed on 10 April 2023). (3) Furthermore, studies conducted in the Ethiopian highlands by [28,29] confirm minimal bias using CMIP5 rcp6 models. The average CMIP5 rcp6 model data annual cycle correlation results with GPCC V6 rainfall (r ≥ 0.6) and with CRU TS v3.22 maximum temperature (r ≥ 0.6) in the Ethiopian highlands indicate the skill of rcp6 in capturing rainfall and temperature (Figure 2a,b). Due to the above justifications and limited capacity, only CMIP5 rcp6 model data sources were used to project the climatic conditions of the NW Ethiopian highlands. It should be noted that the results shown in this study are from the CMIP5 rcp6 scenario only. Nevertheless, when assessing impacts, changes in climate should be examined in the remaining three scenarios as well.
Each of the variables—rainfall, sensible heat flux (SHF), and temperature—was area-averaged for the NW Ethiopian highlands to facilitate temporal and spatial analysis. The procedures and methods devised to calculate soil water content, incorporating Runoff (RO) (not included in the Thornthwaite soil water balance method) and utilizing the Thornthwaite and Mather soil moisture retention table, play a significant role in studying soil moisture in the NW Ethiopian highlands [30]. The soil moisture changes obtained from the computed soil water content were applied to determine the onset date, cessation date, and length of growing period (LGP) for the study area.

2.3. Determining the Best Reference for Validation

Before assessing the models, the optimal reference values for temperature, rainfall, and potential evapotranspiration (PET) variables were determined. For tasmax, the Climate Research Unit Time Series (CRU TS v3.22) interpolated observations compared with the average of 35 stations over the study area from 1981 to 2000. The correlation was found to be r = 0.95, with a minor cold bias of −0.09 °C. This suggests that CRU TS v3.22 is a reliable data source for temperature, suitable for validating and evaluating the CMIP5 rcp6 models’ temperature data sources, and for selecting the most appropriate model.
Comparing Global Precipitation Climatology Center version 6 (GPCC v6) rainfall with the area-average of 35 station observations from 1981 to 2000, a correlation of r = 0.99 was observed, with a minimal dry bias of GPCC −0.6 mm/day based on the scatter plot. Therefore, GPCC v6 at 50 km resolution effectively accounts for Ethiopia’s orographic effect on rainfall. It can be used for validating and evaluating the CMIP5 rcp6 models’ data sources for precipitation and selecting the most suitable model. The selected climate model for rainfall simulation was then applied to derive the soil water content and the length of the growing season in the NW Ethiopian highlands for both the past (1981–2010) and the future (2041–2070). In line with this, GPCC rainfall products have demonstrated good agreement with reference gauge data in studies conducted by [31,32] on the Ethiopian highlands.
PET was calculated using 18 stations from [33,34] in the NW Ethiopian highlands, with Climate Research Unit–Potential Evapotranspiration (CRU-PET) used as a reference for comparison with other datasets and proxies, such as maximum temperature and surface heat flux. The relation between station PET and CRU maximum temperature was significant but considering the known lag of temperature behind solar radiation (Ojo, 1977), it was excluded. In comparison with ECMWF reanalysis-adjusted sensible heat flux (adjshf), a robust relation was indicated by a regression of R2 = 0.95. Weaker correlations were found between station PET and Modern-Era Retrospective Analysis (MERRA) reanalysis-adjusted sensible heat flux. Thus, the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim adjshf served as the validation reference for CMIP5 rcp6-sensible heat flux, following adjustment to match ECMWF.

2.4. Model Validation and Evaluation

The evaluation involved comparing monthly gridded observations, reanalysis-gridded data from GPCC v6 (for rainfall), CRU TS v3.22 (tasmax), and ERA-Interim (adjshf) with station observation data. The comparison revealed a correlation coefficient (r) greater than 0.90, which was used to validate the CMIP5 rcp6 model’s area-averaged data over the NW Ethiopian highlands. Stepwise criteria were applied for comparison, and the models that met most criteria were CCSM4 (for tasmax) and HadGEM2-ES (for precipitation and PET/adjshf). Consequently, these two CMIP5 models were utilized for projecting climate change and calculating the soil water content. Additionally, they were employed to determine the length of the growing season for both the past and future in the study area. Following the stepwise criteria mentioned above, the evaluation and justification of CMIP5 rcp6 model simulations against GPCC V6 (rainfall), CRU TS3.22 (tasmax)-gridded observations, and ERA-Interim (shf) reanalysis-gridded observations were proposed in the context of the Ethiopian highlands’ climate. This approach is consistent with the work of [29,32]. Details on how the reference models were selected for precipitation, temperature, and sensible heat flux are presented in Table 1, Table 2, and Table 3, respectively.

2.5. Climate Change Analysis

Among the various approaches to comprehend climate change, the technique of subtracting past data from future projections is frequently employed. While this method is straightforward, it is plagued by uncertainties stemming from natural climate variability, contingent upon the chosen variable and time periods. In this study, the climatological reference period of 30 years [35] from 1981 to 2010 is utilized. This timeframe strikes a balance between the extended duration necessary for precipitation assessments and a shorter one sufficient for temperature analyses. Another consideration is that trends observed in the most recent 30-year period hold greater relevance for predicting the subsequent 30-year period, particularly considering the acknowledged acceleration of greenhouse warming. To evaluate the signal-to-noise ratio, linear trends are employed instead of comparing differences between two periods. These trends can be expressed as proportional to CO2 concentration or global mean temperature. Despite the accelerating trends in CO2 concentration (Figure 2), the temporal response of climate is assumed to be linear during the study period. This assumption helps eliminate ambiguities arising from second-order responses. Examining trends in rainfall and temperature can reveal ‘drift’ attributed to climate change, whether natural or anthropogenic. Detecting such trends in the data is valuable for decision making at the policy level in agriculture and food security. Trend analysis for temperature, rainfall, PET, and anomalies involved fitting a linear equation through regression. The InStat + v3.36 software [36] applied t-tests to determine the statistical significance of the changes. Bias was assessed by comparing the simulated or projected annual cycle against reference data for all the variables mentioned. Negative bias occurs when the model minus reference is less than 0, representing dry bias for rainfall and cold bias for temperature. Statistical measures of association, indicating the strength or degree of the relationship between variables, were determined spatially and temporally. The Royal Netherlands Meteorological Institute-Climate Explorer (KNMI CE) climexp.knmi.nl webpage and Excel 2007, along with InStat +v3.36 software, were used for this purpose. The 95% significance level served as the threshold for validity.

2.6. Potential Evapotranspiration

Among the various methods employed for calculating potential evapotranspiration (PET), the FAO Penman–Monteith method is recommended by the FAO as the standard [37]. However, in developing countries such as Ethiopia, obtaining the necessary climate parameters for calculating PET is challenging. Due to this difficulty, an alternative method was employed, involving the comparison of other climate variables expected to serve as substitutes for PET. To implement this surrogate method, PET values calculated from station observations by [33,34] in the NW Ethiopian highlands, along with PET estimates from the CRU data source, were used as references. These values were then compared with CRU (tasmax), ECMWF, and MERRA reanalysis data for sensible heat flux (shf), latent heat flux (lhf), and the sum of shf and lhf. The comparison results, which showed good agreement with PET station observations and CRU-PET estimates, were used to produce a surrogate PET. Subsequently, this surrogate PET was utilized to justify the data from CMIP5 rcp6 models for the purpose of evaluation and analysis.

2.7. Soil Water Content Proxy

The spatial distribution and temporal variation of soil water content in the top meter of the Earth’s surface are crucial for numerous environmental studies. Soil water content can be determined through (i) point measurements; (ii) soil water content models; and (iii) remote sensing. In this study, the soil water content calculated from CMIP5 rcp6 models’ data was utilized to estimate soil water balance and changes in soil water content for both the past and the future. Subsequently, using precipitation and a surrogate for PET from CMIP5 rcp6 models, the soil water balance estimate (P-Es-RO) was calculated for both historical and future periods. These estimates were employed for the analysis of soil water balance in the region as shown in Equation (1). RO (runoff) will be calculated as follows:
R O = C r o × P
where the runoff coefficient (Cor) calculate by estimating Cro for each site, taking into account the soil type, steepness, and the slope range, as per [38], Table S1.

2.8. Growing Season Analysis

The LGP methodology draws on insights from [39,40,41,42]. The initiation of the growing season necessitates an increase in soil water content, contingent on both soil and crop type. When the change in soil water first becomes positive (ΔST > 0), the season is considered to have commenced with germination. Conversely, the cessation time is determined when the change in soil water content becomes zero or negative (ΔST ≤ 0). This methodology is applied consistently from past simulations to future projections. Despite Ethiopia being a tropical country, the highlands experience minimum temperatures below 9 °C, a critical threshold beyond which crops can no longer extract available soil water [39,40,41,42].
Upon identifying the CMIP5 model closest to the observations, the study area’s growing season length is estimated using a proxy of water available for storage (P-Es-RO) and RO calculated from (P) and the runoff coefficient (Cor) Equation (1). The calculation considers Cro between 0.15 and 0.32, representing heavy soil, with slope ranging from 2% to 7%, following [29] for details. The computation continues until the last month of the wet season, when water available for storage equals or exceeds the soil’s water-holding capacity. The length of the growing season is then estimated from the CMIP5 rcp6.0 projection of soil water change (ΔST) obtained through the calculation of (P-Es-RO), as outlined below, drawing on [40,43,44,45].
The initial step involves determining the soil, vegetation, crop type, root depth, and water-holding capacity of the study area to obtain corresponding values from the soil moisture retention tables of [30]. Given that the major crops in the NW Ethiopian highlands are corn and cereal grains, the area is classified as moderately deep-rooted crops with soil type silt loam, water holding capacity 200 mm/m, root zone 100 cm, and an applicable soil moisture retention table of 200 mm, used to estimate soil water content.
Accumulated potential water loss represents the potential deficiency of soil moisture associated with moisture contents below the water-holding capacity of the soil. This accumulation increases during dry seasons, decreases during wet seasons due to soil moisture recharge, and equals zero when the soil moisture storage matches the water-holding capacity. The next step involves computing the accumulated potential water loss (Acc. Pot. WL), starting from an estimated value of potential water deficiency in the first month after the rainy season when (PEsRo) turns negative (as per [30]). Subsequent monthly values are obtained by progressively adding (PEsRo) with negative values (Equation (2)).
A c c . P o t . W L i + 1 = A c c . P o t . W L i + 1 + ( P E s R o ) i + 1
where i is the current month and i + 1 the next month from the current month. If the soil moisture exceeds the water-holding capacity, the soil moisture is set equal to the water-holding capacity:
S m j = S m j 1 + ( P E s R o ) ( j )   f o r   S m j 1 + ( P E s R o ) ( j ) < W W H C W W H C                                                         f o r     S m j 1 + ( P E s R o ) ( j ) W W H C  
where Sm(j) is the soil moisture storage (in millimeters) for the j-th month and WWHC is the water-holding capacity of the soil.
The moisture retained in the soil is determined by converting each value of accumulated potential water loss into soil water content storage (ST) with the help of a soil moisture retention table. The change in soil moisture storage is calculated as the difference between the current and the previous month’s soil moisture storage, starting from December of each year and moving into the next growing season. The change in soil moisture is considered positive if the soil moisture has increased and negative if it has decreased. Once ΔST is obtained, the onset date, cessation, and subsequent length of growing period (LGP) are estimated from monthly data. To better resolve the time scale, ΔST is interpolated to daily data using INSTAT v3.36 with a 10-day running mean. The daily onset date, cessation date, and LGP are then estimated. Finally, a minimum temperature threshold is applied, which may induce an early cessation of LGP. The flow chart illustrating how the research steps interlink can be found in Supplementary Information (SI), Figure S1, and a summary of the methods used in the research is presented in Table S2.

3. Results and Discussion

3.1. Changes in Temperature

In the area of study, according to the CCSM4 model projection, the average tasmax is expected to rise from 26.5 °C in the past (1981–2010) to 27.8 °C in the future (2041–2070), indicating a projected increase of 1.3 °C when comparing the past and future tasmax. For average tasmin, the model projection shows a change from 15.6 °C in the past to 16.8 °C in the future, with an expected rise of 1.2 °C in the future. The temperature trend test for both tasmax and tasmin across all time periods (past, future, and forward projection) revealed statistically significant and increasing trends.
The increasing trend for tasmax (Figure 3b) anticipates tasmax rises of about +0.023 °C/yr, consistent with [46,47] using the RCP6 scenario (+0.025 °C/yr). Similar results have been reported in IPCC assessments of tasmax over Ethiopian/NW Ethiopian highlands by [29,48]. The [49] projection of average annual temperature using RCP4.5 at the end of this century is +2.0 °C, almost in the middle of the projection difference (+2.8 °C) obtained from the RCP6 scenario used in this study. It should be noted that RCP4.5 (650 ppm CO2 eq) and RCP6 (850 ppm CO2 eq) are stabilization scenarios without an overshot pathway after 2100. The trend analysis for tasmin (Figure 3c) and the intra-seasonal trend analysis of the hot season (February to May) and the cool season (November to January) show significant changes with temperature increases of +0.0271 °C/yr and +0.0226 °C/yr, respectively. These indicate that the hot season will heat up faster, while the cold season and tasmin will experience slower temperature increases. The annual cycle comparison of past and future tasmax and tasmin shows a high correlation (R2 = 0.99) with higher temperature changes in the future for all months (≈+1.3 °C). Regarding the hot and cool seasons, the hot season in the past (February-May) remains consistent in the future, with an average temperature rise from 29.6 °C to 31.1 °C (+1.5 °C). The cool season, from November to January, also remains consistent in the past and future, with an average temperature rise from 13 °C to 14.3 °C (+1.3 °C) in the future (Figure 3a).
The evaluation of spatial trend maps calculated for the time period 1981–2100 using the CCSM4 model projection indicates fast warming in the NW Ethiopian highlands, with faster warming in the central highlands (>+0.026 °C/yr) slowing down (+0.020 to +0.025 °C/yr) toward the edges of the highlands, as shown in Figure 4a. Furthermore, tasmin shows slow warming in the western part (<+0.021 °C/yr), while in the remaining part of the study area, there is modest warming (+0.021 to 0.024 °C/yr), and the fastest warming is projected to occur in the northern extreme part (+0.024 to +0.025 °C/yr), as depicted in Figure 4b. These trends align with [6,48] and may not be entirely harmful in the region [46]. Nevertheless, projections from [50] of mean annual temperature increase by the 2050s may result in teff, barley, and sorghum yields declining by 8–17%, and wheat and maize increasing by 1–3%.
The spatial distribution of tasmax (Figure 5a) depicts most of the NW Ethiopian highlands tasmax ranging between 23 and 26 °C, with the hottest part from 28 to 31 °C along the western border and the extreme north of the region, while the coolest (<22 °C) part is the south-central part of the study area. The spatial temperature change shows an increase in tasmax in the future by >+1.3 °C, but in the south-central part of the highlands, it will most likely increase faster by +1.45 °C (Figure 5b). The spatial pattern of mean surface air minimum temperature (tasmin) simulated by CCSM4 over Ethiopia clearly identifies cool highlands (<15 °C) surrounded by warmer lowlands (<19 °C). There are also areas in the southern central part that were the coolest (9 to 11 °C) (Figure 6a). The change between past and future CCSM4 minimum air temperatures over the NW Ethiopian highlands indicates warming (+1.13 to +1.37 °C), except for the central western part which warms slowly (+0.97 to +1.13 °C), while the northern extreme part warms faster (>+1.45 °C) (Figure 6b). The coolest zone (9 to 11 °C) in the past is anticipated to warm by +1.2 to 1.4 °C, implying the minimum temperature of this coolest zone will range between 10.2 and 12.4 °C in the future. The rising tasmin may improve the physiological development of crops, especially at the end of the growing season.
There is high confidence in future warming, so there is a need for more information on the impacts of higher temperatures on evaporation (and soil water availability), agriculture, health, and vector-borne diseases [50,51]. The increase in temperature due to climate change has positive impacts on crop production at high elevations in the tropics by extending the length of the growing season, allowing earlier planting of crops, later maturation and harvesting, and even providing the opportunity for double cropping in a single season.

3.2. Changes in Rainfall

The results from the HadGEM2-ES simulation indicate that the mean average annual rainfall in the NW Ethiopian highlands was 1210 mm (1981–2010) and 1267 mm (2041–2070), reflecting an increase of +4.7% in the future. This finding aligns with [29,48,49,50,52,53]. The slope of the annual rainfall in the HadGEM2-ES forward projection (1981–2100) is upward by +0.0023 mm/day, demonstrating a statistically significant trend (Figure 7a), consistent with [29].
In the forward projection (1981–2100) of the seasons, the major rainfall season (JJAS) exhibits a trend of +0.0023 mm/day, and the dry season (ONDJ) shows a trend of +0.0006 mm/day, both indicating a significant increasing trend. The short rainfall season (FMAM) demonstrates a trend of +0.00031 mm/day but is not statistically significant. Anticipating wetter conditions, the livelihood of NW Ethiopian farmers may benefit from the increased rainfall, assisting in improving crop production (Figure 7a). The rainfall annual cycle of the past (1981–2010) and future (2041–2070) shows a unimodal nature, consistent with results reported by [5,54,55]. The main change in the future is a seasonal broadening with more rain in May-June and September-October (Figure 7b). The seasons’ contribution to the annual rainfall during the forward projection (1981–2100) shows JJAS (67.1%), FMAM (22.6%), and ONDJ (10.3%), indicating that JJAS will be the major rainfall season, followed by FMAM, and ONDJ as the dry season of the study area. Commonly, JJAS and FMAM seasons are merged seasons between May and June.
The rainfall spatial trend analysis in the forward projection shows an increasing trend (>+0.0005 mm/day), except for a small area around the eastern Blue Nile basin, which increases by +0.0005 to +0.0010 mm/day. The majority part of NW Ethiopian highlands (western region) exhibits a fast-increasing trend (>+0.0025 mm/day), consistent with the time series analysis of the forward projection, which depicts an increasing trend of rainfall by +0.0023 mm/day (Figure 7c).
In Figure 8a, the spatial distribution of the past depicts a decreasing pattern from the south-central part (4.25 mm/day) of the region to the driest point in the northern extreme (1.75 to 2.25 mm/day), with modest rainfall along the western and eastern boundaries (2.25 to 3.25 mm/day) of the study area. However, most of the region has received sufficient rainfall between 3.25 and 4.75 mm/day, and the Shewa highlands (8–9° N and 37–38° E) around Wolkitie, with elevations between 2000 and 3500 m, are the wettest part (>5.25 mm/day). The mean annual rainfall spatial distribution (middle) of the past and change in the future (all units in mm/day)/right. Figure 8b shows the spatial distribution trend of wetter conditions (>+0.05 mm/day) in the study area, except for the Wollo highlands (east-central part) with no change. The greatest increase in wetter conditions (>+0.2 mm/day) is expected in the west-southern part of the region (including the Wollega and Benshagul highlands in the western Blue Nile basin). Additionally, in the south and west-south part of the highland regions, rainfall amounts between 4.5 and 5.6 mm/day are projected, indicating that this area will continue to have the highest rainfall potential, as confirmed by [54,55,56,57]. Thus, both the time series and spatial analysis of annual rainfall reveal an increase in rainfall in the NW Ethiopian highlands, consistent with the study by [58].

3.3. Changes in Potential Evapotranspiration

PET is valuable for assessing moisture deficits and identifying periods with a high demand for irrigation and water harvesting. PET is derived from the surrogate sensible heat flux (SHF), which is obtained from the CMIP5 RCP6 HadGEM2-ES model’s SHF output, adjusted with adjSHF (+2 mm/day). The annual PET cycle for both the past and the future indicates a mean of 3.64 mm/day, with a minimum of 3.62 mm/day and a maximum of 3.81 mm/day. Conversely, the mean for the future is 3.55 mm/day, with a minimum of 3.51 mm/day and a maximum of 3.70 mm/day. These values align with [37] suggestion of an average PET ranging from 3 mm/day to 5 mm/day for humid and sub-humid regions, with a mean daily temperature of around 20 °C.
The annual PET cycle pattern remains relatively stable from June to October, followed by an increase leading to the maximum during the hot season (February to April) (see Figure 9a). PET correlates with temperature, consistent with the findings of [59] in their study of Ethiopia’s climate. To assess the change in PET from the past to the future, a trend test was conducted. The results indicate a decreasing trend of −0.0003 mm/day from 1981 to 2070 and an increasing trend of +0.0012 mm/day from 2041 to 2070. However, neither change is statistically significant, with p-values greater than 0.1517. Notably, PET exhibits a significant decreasing trend of −0.0036 mm/day in the past (1981–2010). Subsequently, all spatial maps in this chapter are presented in SHF (W/m2/month). To maintain consistency, the presentation is adjusted by converting SHF in W/m2 to PET/SHF in mm. The forward projection (1981–2070) regression trend spatial map for PET/SHF indicates minimal change (−0.0010 to +0.0010 mm/day) in most areas and a decreasing trend (−0.0010 to −0.0031 mm/day) in the north-western extreme (see Figure 9b). In general, the forward projection of the spatial trend aligns with the time series trend.

3.4. Changes in Soil Water Balance

The water balance technique, according to [60], involves quantitatively estimating the water circulation within the hydrosphere, lithosphere, and atmosphere. The soil water balance method is employed to calculate soil water balance by extracting precipitation (P) and PET (Es) data from the HadGEM2-ES CMIP5 rcp6 model for both the past (2041–2070) and the future. Before evaluating the water balance (P-Es) for the past and future, a trend analysis is conducted to determine the statistical significance of the change. The forward projection (2041–2070) and the past reveal a significantly increasing trend of +0.035 mm/day and +0.139 mm/day, respectively. However, the trend test for the future indicates a decreasing trend (−0.017 mm/day), but the change is not statistically significant. The computational procedures used to estimate changes in soil water content (ΔST > 0) are likely to yield reliable results regarding soil water content and its implications in the study area. The results show excess forcing (ΔST > 0) during the months of June to September (4 months) and deficit forcing (ΔST ≤ 0) from October to May (8 months), both in the past and in the future (Table 4).

3.5. Variability in Length of Growing Season

The annual cycle of soil water change (ΔST) was computed for both past and future years. In the NW Ethiopian highlands, estimates of the average reserve of soil water storage were 82 mm/month in the past and a projected 84 mm/month in the future. Both are positively affected by ΔST > 0, starting from the onset day of the year (DOY). However, soil water content changes from the first day (offset DOY) (ΔST < 0) influence the driving soil water content to become deficient. The trend analysis test revealed a decreasing trend in onset date (−0.021 day/yr) and cessation date (−0.300 day/yr) in the past. Still, LGP showed an increasing trend (+0.033 day/yr). In the future, there is a decreasing trend in cessation date (−0.039 day/yr) and LGP (−0.335 day/yr), while the onset date shows an increasing trend (+0.295 day/yr). However, these changes are not statistically significant for any of the parameters at a 0.05 significance level. The statistical median result for the onset date was 139 DOY, for the cessation date it was 263 DOY, and for LGP it was 130 days in the past. In the future, the onset date, cessation date, and LGP are likely to be 137 DOY, 270 DOY, and 136 days, respectively.
Our results align with [61], who emphasized the significant variability in seasonal rainfall, especially during the cropping season, potentially resulting in moisture stress for maize plants and impacting the LGP. A slightly similar result for LGP, ranging between 137 and 205 days for the NW Ethiopian highlands, has been reported in the mapping of climate vulnerability and poverty in Africa (current condition of 2000) by [62,63] also suggested a planting date in April and an LGP of 180 days for maize (grain) in East Africa’s high-altitude areas. The LGP estimation obtained in both time periods has a minimum LGP ≥ 118 days, which is above the threshold LGP (90 days). Ref. [39], most crops grown in Ethiopia, with the exception of some pulses and very low-yielding varieties of Teff (a common cereal food in Ethiopia) and wheat, require a growing period of at least 90 days. Additionally, ref. [64] recommended a 90-day LGP as the optimum LGP in the case of Ethiopia. Figure 9 illustrates the evaluation of the variability of the length of growing seasons, showing highly variable onset dates, followed by less variable LGP, and offset dates in the past. In the future, offset dates are likely to be highly variable, while LGP and onset dates are less variable. In general, the evaluation of the onset date, end date, and LGP of crop growing seasons under a climate change scenario emphasizes the importance of early planting in May and harvesting at the end of September, with LGP suitable for single and, to some extent, double cropping both in the past and future. There is a slightly longer growing season in the future, which is likely to benefit agricultural production in the NW Ethiopian highlands.

4. Conclusions and Summary

This study provides in-depth analyses of climate change and its impact on soil water balance and the length of growing seasons in the NW Ethiopian highlands. Utilizing data sets from station observations, gridded observations, reanalysis, and simulations of CMIP5 rcp6 models, various analysis tools, including trends, spatial maps, and measures of bias and dispersion, were employed for a comprehensive result analysis. The validation criteria identified CCSM4 (for tasmax) and HadGEM2-ES (for precipitation and PET/adjshf) as the models meeting most criteria, serving as data sources for climate change, variability, extremes, and soil water balance analysis.
Temperature exhibits a significant increasing trend across all periods, with tasmax increasing faster than tasmin, and the hot season surpassing the cold season in both historical and future contexts. The spatial maps for tasmax and tasmin projections indicate a faster warming trend in the region. The rainfall trend analysis in the NW Ethiopian highlands is mixed, with significant increases in JJAS season across all periods and non-significant trends in other seasons and annual rainfall. The spatial distribution aligns with the time series trend, suggesting a wetter future, particularly in the southwest.
The extreme decrease in minimum rainfall necessitates substantial adjustments in agriculture and water management. Conversely, the maximum rainfall changes benefits agriculture, except in flood-prone areas. The extreme probabilities for the FMAM season suggest that it may not be viable for rain-fed crops alone but could support early planting. The ONDJ season, anticipated to bring substantial rainfall, could serve as an additional growing season, or it could extend JJAS, but early warning is essential to avoid disruptions during harvest.
PET trends are mixed, with significant changes only in the past. Comparisons with rainfall and temperature indicate low PET rates during low temperatures and vice versa. The water balance estimation reveals a significant increasing trend in the past but a non-significant decreasing trend in the future. The annual total water balance conditions show PET exceeding rainfall in both periods.
This study underscores the importance of early planting in April/May and harvesting in September, suitable for single and, to some extent, double cropping in the past and future. The growing season variability indicates highly variable LGP dates, with onset and offset dates less variable in the past. In the future, onset dates are expected to be highly variable, while LGPs and offset dates remain less variable. Overall, this study suggests crucial implications for farming in the NW Ethiopian highlands. Farmers should adapt by adjusting planting and harvesting times, benefiting from the projected increase in rainfall. Efficient water management, particularly in irrigation, can capitalize on the observed rise in soil water balance, enhancing agricultural resilience. Diversifying crops and introducing new varieties aligns with the anticipated lengthening of the growing season, optimizing productivity. Integrating these findings into climate-resilient planning is essential for regional authorities, with a focus on sustainable water resource management and infrastructure development. Ongoing research remains vital for refining adaptation strategies and staying abreast of evolving climate dynamics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cli11120243/s1, Figure S1: Flow chart of how the research steps interlink is illustrated. Table S1: Runoff coefficient: empirical coefficient representing fraction of rainfall that runs off (adapted from [42]). Table S2: Summary for the methods how the research is performed.

Author Contributions

G.B.T.: Conceptualization, Methods, Supervision, Analysis, Writing—Original Draft, Writing—reviewing & editing. A.A.A.: Conceptualization, Methods, Data extraction, Data analysis (Figures, Tables), Writing—Original Draft, Writing—reviewing & editing. S.E.D.: Overall supervision, conceptualization, Writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this article are provided based on the inquiries related to the data may be directed to Gashaw Bimrew Tarekegn, and Addis A. Alaminie at [email protected] or [email protected].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Akinseye, F.M.; Ajayi, V.O.; Oladitan, T.O. Assessing the impacts of climate variability on crop yield over Sudano-Sahelian zone in Nigeria. Int. J. Agric. Sci. 2013, 1, 91–98. [Google Scholar]
  2. Mesike, C.S.; Esekhade, T.U. Rainfall variability and rubber production in Nigeria. Afr. J. Environ. Sci. Technol. 2014, 8, 54–57. [Google Scholar] [CrossRef]
  3. Alemayehu, A.; Bewket, W. Local climate variability and crop production in the central highlands of Ethiopia. Environ. Dev. 2016, 19, 36–48. [Google Scholar] [CrossRef]
  4. Wubie, A.; Tadesse, Z.; Tesfaye, G.; Lakew, A. Fish Distribution and Abundance in the Upper Awash River, West Shewa, Ethiopia. J. Sci. Sustain. Dev. (JSSD) 2017, 5, 109–118. [Google Scholar]
  5. Cheung, W.H.; Senay, G.B.; Singh, A. Trends and spatial distribution of annual and seasonal rainfall in Ethiopia. Int. J. Climatol. A J. R. Meteorol. Soc. 2008, 28, 1723–1734. [Google Scholar] [CrossRef]
  6. Alebachew, A.; Woldeamlak, B. A Climate Change Country Assessment Report for Ethiopia. In Submitted to Forum for Environment; Epsilon International R&D: Addis Ababa, Ethiopia, 2011. [Google Scholar]
  7. Benedicta, Y.F.-M. Modelling Maize (Zea mays L.) Productivity and Impact of Climate Change on yield and Nutrient Utilization in Sub-humid Ghana; Ecology and Development Series; Zentrum für Entwicklungsforschung: Bonn, Germany, 2012. [Google Scholar]
  8. Hultine, K.R.; Hernández-Hernández, T.; Williams, D.G.; Albeke, S.E.; Tran, N.; Puente, R.; Larios, E. Global change impacts on cacti (Cactaceae): Current threats, challenges and conservation solutions. Ann. Bot. 2023, mcad040. [Google Scholar] [CrossRef]
  9. Lobell, D.B.; Gourdji, S.M. The influence of climate change on global crop productivity. Plant Physiol. 2012, 160, 1686–1697. [Google Scholar] [CrossRef] [PubMed]
  10. Gao, Y.; Wang, M.; Jiang, L.; Zhao, F.; Gao, F.; Zhao, H. Dynamics of carbon budget and meteorological factors of a typical maize ecosystem in Songnen Plain, China. Span. J. Agric. Res. 2023, 21, e0301. [Google Scholar] [CrossRef]
  11. Rosenzweig, C.; Strzepek, K.M.; Major, D.C.; Iglesias, A.; Yates, D.N.; McCluskey, A.; Hillel, D. Water resources for agriculture in a changing climate: International case studies. Glob. Environ. Chang. A 2004, 14, 345–360. [Google Scholar] [CrossRef]
  12. Horgan, F.G. Potential for an impact of global climate change on insect herbivory in cereal crops. In Crop Protection under Changing Climate; Springer: Cham, Switzerland, 2020; pp. 101–144. [Google Scholar]
  13. Lesk, C.; Anderson, W.; Rigden, A.; Coast, O.; Jägermeyr, J.; McDermid, S.; Davis, K.F.; Konar, M. Compound heat and moisture extreme impacts on global crop yields under climate change. Nat. Rev. Earth Environ. 2022, 3, 872–889. [Google Scholar] [CrossRef]
  14. Hatfield, J.L.; Prueger, J.H. Temperature Extremes: Effect on Plant Growth and Development. Weather. Clim. Extrem. 2015, 10, 4–10. [Google Scholar] [CrossRef]
  15. Kukal, M.S.; Irmak, S. US agro-climate in 20th century: Growing degree days, first and last frost, growing season length, and impacts on crop yields. Sci. Rep. 2018, 8, 6977. [Google Scholar] [CrossRef] [PubMed]
  16. Gautam, H.R.; Bhardwaj, M.L.; Kumar, R. Climate change and its impact on plant diseases. Curr. Sci. 2013, 105, 1685–1691. [Google Scholar]
  17. Robertson, A.W.; Ines, A.V.; Hansen, J.W. Downscaling of seasonal precipitation for crop simulation. J. Appl. Meteorol. Climatol. 2007, 46, 677–693. [Google Scholar] [CrossRef]
  18. World Bank. Ethiopian: Country Economic Memorandum, Background Report: Fourth Ethiopia: A Regional Characterization Assessing Ethiopia Growth Potential and Development Obstacles; World Bank: Washington, DC, USA, 2004. [Google Scholar]
  19. World Bank. Ethiopia: Climate Risk Factsheet (Draft Summary Report); World Bank: Addis Ababa, Ethiopia, 2009. [Google Scholar]
  20. Thornton, P.K.; Jones, P.G.; Alagarswamy, G.; Andresen, J. Spatial variation of crop yield response to climate change in East Africa. Glob. Environ. Chang. 2009, 19, 54–65. [Google Scholar] [CrossRef]
  21. Gebrehiwot, T.; Van Der Veen, A. Farm level adaptation to climate change: The case of farmer’s in the Ethiopian Highlands. Environ. Manag. 2013, 52, 29–44. [Google Scholar] [CrossRef]
  22. Asrat, P.; Simane, B. Farmers’ perception of climate change and adaptation strategies in the Dabus watershed, North-West Ethiopia. Ecol. Process. 2018, 7, 7. [Google Scholar] [CrossRef]
  23. Senbeta, A.F.; Worku, W. Ethiopia’s wheat production pathways to self-sufficiency through land area expansion, irrigation advance, and yield gap closure. Heliyon 2023, 9, e20720. [Google Scholar] [CrossRef]
  24. Trouet, V.; Van Oldenborgh, G.J. KNMI Climate Explorer: A Web-Based Research Tool for KNMI Climate Explorer (CE) Instruction Manual. 2013. Available online: http://climexp.knmi.nl/help.cgi (accessed on 1 December 2012).
  25. Fujino, J.; Nair, R.; Kainuma, M.; Masui, T.; Matsuoka, Y. Multigas mitigation analysis on stabilization scenarios using aim global model. Energy J. Spec. Issue 2006, 3, 343–354. [Google Scholar]
  26. Hijioka, Y.; Matsuoka, Y.; Nishimoto, H.; Masui, T.; Kainuma, M. Global GHG emission scenarios under GHG concentration stabilization targets. J. Glob. Environ. Eng. 2008, 13, 97–108. [Google Scholar]
  27. Taylor, K.E.; Stouffer, R.J.; Meehl, G.A. An overview of CMIP5 and the experiment design. Bull Amer. Met. Soc. 2012, 93, 485–498. [Google Scholar] [CrossRef]
  28. Becker, A.; Finger, P.; Meyer-Christoffer, A.; Rudolf, B.; Schamm, K.; Schneider, U.; Ziese, M. A description of the global land-surface precipitation data products of the Global Precipitation Climatology Centre with sample applications including centennial (trend) analysis from 1901–present. Earth Syst. Sci. Data (ESSD) 2013, 5, 71–99. [Google Scholar] [CrossRef]
  29. Jury, M.R. Statistical evaluation of CMIP5 climate change model simulations for the Ethiopian highlands. Int. J. Climatol. 2015, 35, 37–44. [Google Scholar] [CrossRef]
  30. Thornthwaite, C.W.; Mather, J.R. Instructions and Tables for Computing Evapotranspiration and the Water Balance, 5th ed.; Climatology: Centerton, NJ, USA, 1957. [Google Scholar]
  31. Fuchs, T.; Schneider, U.; Rudolf, B. Global Precipitation Analysis Products of the GPCC. 2007. Available online: http://www.dwd.de/en/FundE/Klima/KLIS/int/GPCC/ReportsPublications/QR/GPCC-intro-products-2007.pdf (accessed on 1 December 2012).
  32. Dinku, T.; Connor, S.J.; Ceccato, P.; Ropelewski, C.F. Comparison of global gridded precipitation products over a mountainous region of Africa. Int. J. Climatol. A J. R. Meteorol. Soc. 2008, 28, 1627–1638. [Google Scholar] [CrossRef]
  33. Muluneh, B. Spatial Variability of Impact of Climate Change on Water Availability in BlueNile Basin. Master’s Thesis, Arbaminch University, Arbaminch, Ethiopian, 2008. [Google Scholar]
  34. Banchiamlak, K. Analysis of the Variability of Rainfall Distribution and Evapotranispiration over Eastern Amahara Region. Master’s Thesis, Arbaminch University, Arbaminch, Ethiopia, 2009. [Google Scholar]
  35. World Meteorology Organization. Scientists Urge More Frequent Updates of 30-Year Climate Baseline to Keep Pace with the Rapid Climate Change; Press Release No. 997; World Meteorology Organization: Heidelberg, Germany, 2014. [Google Scholar]
  36. Stern, R.; Rijks, D.; Dale, I.; Knock, J. INSTAT Climatic Guide; University of Reading: Reading, UK, 2006; 330p, Available online: http://www.reading.ac.uk/ssc/software/instat/climatic.pdf/oct.27/2010 (accessed on 1 December 2012).
  37. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration—Guidelines for Computing Crop Water Requirements—FAO Irrigation and Drainage Paper 56; Food and Agricultural Organization: Rome, Italy, 1998. [Google Scholar]
  38. Henricksen, B.L.; Durkin, J.W. Moisture Availability, Cropping Period and the Prospects for Early Warning of Famine in Ethiopia; Computer Unit, ILCA: Addis Ababa, Ethiopia, 1987. [Google Scholar]
  39. Masresha, E. Growing season belts of Ethiopia. In Proceedings of the National Sensitization Workshop on Agrometeorology and GIS, Addis Ababa, Ethiopia, 17–18 December 2001; EARO: Addis Ababa, Ethiopia, 2003; pp. 35–36. [Google Scholar]
  40. Belay, S. Agroclimatic Analysis in Relation to Crop Production; Ethiopian Agricultural Research Organization of Ethiopia: Bahir Dar, Ethiopia, 2001. [Google Scholar]
  41. Araya, S.D.; Keesstra, L.; Stroosnijder, L. A new agro-climatic classification for crop suitability zoning in northern semi-arid Ethiopia. Agric. For. Meteorol. 2010, 150, 1057–1064. [Google Scholar] [CrossRef]
  42. McBean, E.A.; Rovers, F.A.; Farquhar, G.J. Solid Waste Landfill Engineering and Design; Prentice Hall PTR: Englewood Cliffs, NJ, USA, 1995. [Google Scholar]
  43. Food and Agricultural Organization. Agro-Ecological Zoning Guidelines, FAO Soils Bulletin 76, Soil Resources, Management and Conservation Service; FAO Land and Water Development Division: Rome, Italy, 1996. [Google Scholar]
  44. Thornton, P.K.; Jones, P.G.; Owiyo, T.; Kruska, R.L.; Herrero, M.T.; Kristjanson, P.M.; Notenbaert, A.M.O.; Bekele, N.; Omolo, A. Mapping Climate Vulnerability and Poverty in Africa; Department for International Development: Nairobi, Kenya, 2006; p. 171. [Google Scholar]
  45. Thornthwaite, C.W. An approach toward a rational classification of climate. Geographical Review. Am. Geogr. Soc. 1948, 38, 55–94. [Google Scholar]
  46. Nata, T. Surface Waters Potential of the Hantebet Basin, Tigray, Northern Ethiopia. Agric. Eng. Int. CIGRE J. 2006, VIII, LW0501. [Google Scholar]
  47. Tesfaye, M. Water Balance Applications in the Identification and Analysis of Agricultural Droughts in South Wollo Zone, Ethiopia. Ph.D. Thesis, The University the Andhra University, Visakhapatnam, India, 2014. [Google Scholar]
  48. Knutson, T.R.; Zeng, F.; Wittenberg, A.T. Multi-Model Assessment of Regional Surface Temperature Trends: CMIP3 and CMIP5 20th Century Simulations. J. Clim. 2013. [Google Scholar] [CrossRef]
  49. Federal Democratic Republic of Ethiopia. Climate Resilient Green Economy: Mission Statement; FDRE: Addis Ababa, Ethiopia, 2011. [Google Scholar]
  50. National Meteorological Services Agency. Climate Change National Adaptation Programme of Action (NAPA) of Ethiopia; NMSA: Addis Ababa, Ethiopia, 2007. [Google Scholar]
  51. Intergovernmental Panel on Climate Change. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Pachauri, R.K., Meyer, L.A., Eds.; IPCC: Geneva, Switzerland, 2014; 151p. [Google Scholar]
  52. Conway, D.; Schipper, L.; Mahmud, Y.; Menale, K.; Aurélie, P.; Bereket, K. Reducing Vulnerability in Ethiopia: Addressing the Implications of Climate Change Integration of Results from Phase I; Report prepared for the UK Department for International Development (DFID) and the Canadian International Development Agency (CIDA); Ethiopian Development Research Institute (EDRI): Addis Ababa, Ethiopia, 2007. [Google Scholar]
  53. Cramer, W.; Guiot, J.; Marini, K.; Azzopardi, B.; Balzan, M.V.; Cherif, S.; Doblas-Miranda, E.; Dos Santos, M.J.P.L.; Drobinski, P.; Fader, M.; et al. MedECC 2020 summary for policymakers. In Climate and Environmental Change in the Mediterranean Basin–Current Situation and Risks for the Future. First Mediterranean Assessment Report; Union for the Mediterranean, Plan Bleu, UNEP/MAP: Marseille, France, 2020; pp. 11–40. [Google Scholar]
  54. Conway, D.; Schipper, E.L.F. Adaptation to Climate Change in Africa: Challenges and Opportunities; (Draft) Final Revised Version, June, 2012; Elsevier: Addis Ababa, Ethiopia, 2011. [Google Scholar]
  55. Jury, M.R. Ethiopian highlands crop-climate prediction 1979–2009. J. Appl. Meteorol. Climatol. 2013, 52, 1116–1126. [Google Scholar] [CrossRef]
  56. National Meteorological Services Agency. Climatic and Agro-Climatic Resources of Ethiopia; Meteorological Research Report Series; National Meteorological Services Agency: Addis Ababa, Ethiopia, 1996; Volume 1, pp. 1–37. [Google Scholar]
  57. Yilma, S.; Ulrich, Z. Recent changes in rainfall and rainy days in Ethiopia. Int. J. Climatol. 2004, 24, 973–983. [Google Scholar] [CrossRef]
  58. Segele, Z.T.; Lamb, P.J. Characterization and variability of JJAS rainy season over Ethiopia. Meteorol. Atmos. Phys. 2005, 89, 153–180. [Google Scholar] [CrossRef]
  59. Gizaw, M.T. High-Resolution Monthly Rainfall Database for Ethiopia: Homogenization, Reconstruction, and Gridding. J. Climatol. 2012, 25, 8422–8443. [Google Scholar]
  60. McSweeney, C.; New, M.; Lizcano, G. UNDP Climate Change Country Profiles—Ethiopia; United Nations Development Programme: New York, NY, USA, 2008. [Google Scholar]
  61. Fazzini, M.; Carlo, B.; Billi, P. (Eds.) The Climate of Ethiopia: Landscapes and Landforms of Ethiopia, World Geomorphological Landscapes; Springer: Dordrecht, The Netherlands, 2015. [Google Scholar] [CrossRef]
  62. Subrahmanyam, V.P. Water Balance and Its Applications; Andhra University Press: Waltair Visakhapatnam, India, 1980. [Google Scholar]
  63. Markos, D.; Worku, W.; Mamo, G. Spatio-temporal variability and rainfall trend affects seasonal calendar of maize production in southern central Rift Valley of Ethiopia. PLoS Clim. 2023, 2, e0000218. [Google Scholar] [CrossRef]
  64. Ayele, G.; Jabbar, M.A.; Teklewold, H.; Mulugeta, E.; Kebede, G. Seasonal and Inter-Market Differences in Prices of Small Ruminants in Ethiopia; Research Reports 182171; International Livestock Research Institute: Nairobi, Kenya, 2006. [Google Scholar]
Figure 1. Location of NW Ethiopian highlands and distribution of the 35 meteorological stations utilized in this study.
Figure 1. Location of NW Ethiopian highlands and distribution of the 35 meteorological stations utilized in this study.
Climate 11 00243 g001
Figure 2. Average CMIP5 rcp6 correlation (a) with GPCC V6 for rainfall and (b) with CRU TS v3.22 for tasmax.
Figure 2. Average CMIP5 rcp6 correlation (a) with GPCC V6 for rainfall and (b) with CRU TS v3.22 for tasmax.
Climate 11 00243 g002
Figure 3. Annual cycle comparison of CCSM4 projection of past (1981–2010) (a) and the future (2041–2070) and regression trend test of CCSM4 forward projection (b) for tasmax and (c) for tasmin from 1981 to 2100.
Figure 3. Annual cycle comparison of CCSM4 projection of past (1981–2010) (a) and the future (2041–2070) and regression trend test of CCSM4 forward projection (b) for tasmax and (c) for tasmin from 1981 to 2100.
Climate 11 00243 g003
Figure 4. Spatial trend projection of CCSM4 for (a) tasmax and (b) tasmin from 1981 to 2100.
Figure 4. Spatial trend projection of CCSM4 for (a) tasmax and (b) tasmin from 1981 to 2100.
Climate 11 00243 g004
Figure 5. Average maximum temperature (Tasmax): (a) spatial distribution of the past and (b) spatial change in the future.
Figure 5. Average maximum temperature (Tasmax): (a) spatial distribution of the past and (b) spatial change in the future.
Climate 11 00243 g005
Figure 6. Average minimum temperature (tasmin): (a) spatial distribution of the past and (b) spatial change in the future.
Figure 6. Average minimum temperature (tasmin): (a) spatial distribution of the past and (b) spatial change in the future.
Climate 11 00243 g006
Figure 7. (a) Regression trend of rainfall forward projection from 1981 to 2100; (b) annual cycle precipitation monthly relative (%) curve from the HadGEM2-ES projection of past (1981–2010) and future (2041–2070); and (c) spatial trend of annual mean rainfall using Hadgem2-ES forward projection from 1981 to 2100.
Figure 7. (a) Regression trend of rainfall forward projection from 1981 to 2100; (b) annual cycle precipitation monthly relative (%) curve from the HadGEM2-ES projection of past (1981–2010) and future (2041–2070); and (c) spatial trend of annual mean rainfall using Hadgem2-ES forward projection from 1981 to 2100.
Climate 11 00243 g007
Figure 8. Spatial trend of annual mean rainfall using Hadgem2-ES for the past (a) and comparison of the difference between the past and the future (b).
Figure 8. Spatial trend of annual mean rainfall using Hadgem2-ES for the past (a) and comparison of the difference between the past and the future (b).
Climate 11 00243 g008
Figure 9. Annual cycle of sensible heat flux of the past (1981–2010) and the future (2041–2070) (a), and spatial trend of sensible heat flux using Hadgem2-ES for the past (b).
Figure 9. Annual cycle of sensible heat flux of the past (1981–2010) and the future (2041–2070) (a), and spatial trend of sensible heat flux using Hadgem2-ES for the past (b).
Climate 11 00243 g009
Table 1. NW Ethiopian highlands CMIP5 rcp6 model’s rainfall stepwise evaluation taking GPCC v6 as reference from 1981–2010. CMIP5 rcp6 rainfall models validation taking GPCC V6 as reference from 1981–2010.
Table 1. NW Ethiopian highlands CMIP5 rcp6 model’s rainfall stepwise evaluation taking GPCC v6 as reference from 1981–2010. CMIP5 rcp6 rainfall models validation taking GPCC V6 as reference from 1981–2010.
NoModelGPCC v6(r)JJAS- Diff. (GPCC v6 vs. rcp6)Season Peak- JA/Day)Annual Pattern FitnessRemark
1bcc-csm1-10.65−3.23.8poorScreened out
2bcc-csm1-1-m0.90−1.46.3poorScreened out
3CCSM40.76−2.04.3poorScreened out
4CESM1-CAM50.75−0.95.7poorScreened out
5CSIRO-Mk3-6-00.91−1.18.8moderateScreened out
6FIO-ESM0.80−2.14.6poorScreened out
7GFDL-CM30.87−1.56.2poorScreened out
8GFDL-ESM2G0.86−0.36.7moderateScreened out
9GFDL-ESM2M0.87−0.56.4moderateScreened out
10GISS-E2-H_p10.96−5.33.2poorScreened out
11GISS-E2-H_p20.96−5.52.8poorScreened out
12GISS-E2-H_p30.96−4.74.0poorScreened out
13GISS-E2-R_p10.97−5.62.5poorScreened out
14GISS-E2-R_p20.97−5.72.6poorScreened out
15GISS-E2-R_p30.96−5.32.7poorScreened out
16HadGEM2-AO0.96−0.48.1high2nd selected
17HadGEM2-ES0.97−0.58.0high1st selected
18IPSL-CM5A-LR0.89−2.76.3poorScreened out
19IPSL-CM5A-MR0.90−3.35.3poorScreened out
20MIROC50.9811.521.1poorScreened out
21MIROC-ESM0.86−0.57.3highScreened out
22MIROC-ESM-CHEM0.88−0.57.4highScreened out
23MRI-CGCM30.93−2.56.0moderateScreened out
24NorESM1-M0.66−2.73.4poorScreened out
25NorESM1-ME0.62−2.33.3poorScreened out
Table 2. CMIP5 rcp6 models’ maximum temperature (tasmax) stepwise evaluation taking CRU TS3.22 as reference from 1981–2010.
Table 2. CMIP5 rcp6 models’ maximum temperature (tasmax) stepwise evaluation taking CRU TS3.22 as reference from 1981–2010.
CMIP5 rcp6 Tasmax Model Stepwise Evaluation Using CRUTS 3.22 as Reference from 1981–2010
NoModel TypeCRUTS 3.22 Tasmax (r)Seasonal Peak/MAM Diff. (CRUTS 3.22 vs. rcp6)Annual Pattern FitnessRemark
1bcc-csm1-10.902.0PoorScreened out
2bcc-csm1-1-m0.35−1.4PoorScreened out
3CCSM40.98−0.5High1st selected
4CESM1-CAM50.90−2.3PoorScreened out
5CSIRO-Mk3-6-00.43−1.8PoorScreened out
6FIO-ESM0.811.6PoorScreened out
7GFDL-CM30.903.6PoorScreened out
8GFDL-ESM2G0.91−1.0ModerateScreened out
9GFDL-ESM2M0.92−0.7Moderate4th selected
10GISS-E2-H_p10.944.0PoorScreened out
11GISS-E2-H_p20.904.0PoorScreened out
12GISS-E2-H_p30.953.9PoorScreened out
13GISS-E2-R_p10.974.7PoorScreened out
14GISS-E2-R_p20.934.8PoorScreened out
15GISS-E2-R_p30.954.8PoorScreened out
16HadGEM2-AO0.910.2Moderate3rd selected
17IPSL-CM5A-LR0.20−1.1PoorScreened out
18IPSL-CM5A-MR0.221.1PoorScreened out
19MIROC50.93−3.3PoorScreened out
20MIROC-ESM0.933.2PoorScreened out
21MIROC-ESM-CHEM0.942.8PoorScreened out
22MRI-CGCM30.850.9PoorScreened out
23NorESM1-M0.940.3High2nd selected
Table 3. Stepwise evaluation results comparing CMIP5 rcp6 model sensible heat flux (SHF) with adjshf ECMWF-ERA-int. to simulate adjshf for the NW Ethiopian highlands based on the period 1981–2010.
Table 3. Stepwise evaluation results comparing CMIP5 rcp6 model sensible heat flux (SHF) with adjshf ECMWF-ERA-int. to simulate adjshf for the NW Ethiopian highlands based on the period 1981–2010.
NoCMIP5 rcp6 ModelsERA Vs rcp6
PET/Adjshf (r)
Seasonal
Diff. FMA (ERA vs. rcp6)
Annual Cycle
Fitness
Remark
1bcc-csm1-1-m0.790.46moderateFiltered out
2bcc-csm1-10.730.72poorFiltered out
3CCSM40.960.63high2nd selected
4CESM1-CAM50.780.06moderateFiltered out
5CSIRO-Mk3-6-00.791.04poorFiltered out
6FIO-ESM0.950.81high4th selected
7GFDL-CM30.921.26poorFiltered out
8GFDL-ESM2G0.90.92moderateFiltered out
9GFDL-ESM2M0.920.88moderateFiltered out
10GISS-E2-H_p10.951.48poorFiltered out
11GISS-E2-H_p20.961.32poorFiltered out
12GISS-E2-H_p30.971.21poorFiltered out
13GISS-E2-R_p10.951.29poorFiltered out
14GISS-E2-R_p20.971.11poorFiltered out
15GISS-E2-R_p30.971.04moderateFiltered out
16HadGEM2-AO0.930.34highFiltered out
17HadGEM2-ES0.960.24very high1st selected
18IPSL-CM5A-LR0.871.52poorFiltered out
19IPSL-CM5A-MR0.741.34poorFiltered out
20MIROC50.95−0.95moderateFiltered out
21MIROC-ESM0.870.11moderateFiltered out
22MIROC-ESM-CHEM0.850.05moderateFiltered out
23MRI-CGCM30.80.31moderateFiltered out
24NorESM1-M0.950.58high3rd selected
25NorESM1-ME0.910.44highFiltered out
Table 4. Computation of monthly Acc. Pot. WL, soil water storage, soil water storage change, and actual evapotranspiration (mm/month).
Table 4. Computation of monthly Acc. Pot. WL, soil water storage, soil water storage change, and actual evapotranspiration (mm/month).
Time PeriodParameterJanuaryFebruaryMarchAprilMayJuneJulyAugustSeptemberOctoberNovemberDecember
PastRainfall7112989147197217231171742314
RO2.23.59.328.547.063.069.473.954.723.77.44.5
(1981–2010)PET/Es110121126123112928181869399103
P-Es-RO−105−114−106−62−1242677630−43−83−93
ACC.P.WL−325−439−545−607−619 −43−126−220
ST382213985011719320016110566
∆ST−28−16−9−4−14267767−39−56−39
AE (P + ∆ST)3527389311292818186937953
FuturePrecipitation892692161211216227188862615
RO2.62.98.329.451.567.569.172.660.227.58.34.8
(2041–2070)PET/Es123140152139111948482848897108
P-Es-RO−118−134−134−76−249637244−30−79−98
ACC.P.WL−325−459−593−669−671 −30−109−207
ST392010775611919220017211570
∆ST−31−19−10−304963728−28−57−45
AE (P + ∆ST)3928369511194848284888360
Ro = runoff, PET = potential evapotranspiration, Es = actual evapotranspiration, ACC.P.WL = accumulated potential water loss, ST = soil moisture storage, ∆ST = soil moisture change.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tarekegn, G.B.; Alaminie, A.A.; Debele, S.E. Linking Climate Change Information with Crop Growing Seasons in the Northwest Ethiopian Highlands. Climate 2023, 11, 243. https://doi.org/10.3390/cli11120243

AMA Style

Tarekegn GB, Alaminie AA, Debele SE. Linking Climate Change Information with Crop Growing Seasons in the Northwest Ethiopian Highlands. Climate. 2023; 11(12):243. https://doi.org/10.3390/cli11120243

Chicago/Turabian Style

Tarekegn, Gashaw Bimrew, Addis A. Alaminie, and Sisay E. Debele. 2023. "Linking Climate Change Information with Crop Growing Seasons in the Northwest Ethiopian Highlands" Climate 11, no. 12: 243. https://doi.org/10.3390/cli11120243

APA Style

Tarekegn, G. B., Alaminie, A. A., & Debele, S. E. (2023). Linking Climate Change Information with Crop Growing Seasons in the Northwest Ethiopian Highlands. Climate, 11(12), 243. https://doi.org/10.3390/cli11120243

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop