Original Research ARTICLE
Determining the Pace and Magnitude of Lake Level Changes in Southern Ethiopia Over the Last 20,000 Years Using Lake Balance Modeling and SEBAL
- 1Department of Geosciences, Eberhard Karls University Tübingen, Tübingen, Germany
- 2Senckenberg Centre for Human Evolution and Paleoenvironment (S-HEP), Tübingen, Germany
- 3Department of Climate Geochemistry, Max Planck Institute for Chemistry, Mainz, Germany
- 4Earth Observation Centre, German Aerospace Center (DLR), Wessling, Germany
- 5Institute of Geography Education, University of Cologne, Cologne, Germany
- 6School of Earth Sciences, Addis Ababa University, Addis Ababa, Ethiopia
- 7Institute of Geography, University of Leipzig, Leipzig, Germany
- 8Institute of Geosciences, University of Potsdam, Potsdam, Germany
The Ethiopian rift is known for its diverse landscape, ranging from arid and semi-arid savannahs to high and humid mountainous regions. Lacustrine sediments and paleo-shorelines indicate water availability fluctuated dramatically from deep fresh water lakes, to shallow highly alkaline lakes, to completely desiccated lakes. To investigate the role lakes have played through time as readily available water sources to humans, an enhanced knowledge of the pace, character and magnitude of these changes is essential. Hydro-balance models are used to calculate paleo-precipitation rates and the potential pace of lake level changes. However, previous models did not consider changes in hydrological connectivity during humid periods in the rift system, which may have led to an overestimation of paleo-precipitation rates. Here we present a comprehensive hydro-balance modeling approach that simulates multiple rift lakes from the southern Ethiopian Rift (lakes Abaya, Chamo, and paleo-lake Chew Bahir) simultaneously, considering their temporal hydrological connectivity during high stands of the African Humid Period (AHP, ~15–5 ka). We further used the Surface Energy Balance Algorithm for Land (SEBAL) to calculate the evaporation of paleo-lake Chew Bahir's catchment. We also considered the possibility of an additional rainy season during the AHP as previously suggested by numerous studies. The results suggest that an increase in precipitation of 20–30% throughout the southern Ethiopian Rift is necessary to fill paleo-lake Chew Bahir to its overflow level. Furthermore, it was demonstrated that paleo-lake Chew Bahir was highly dependent on the water supply from the upper lakes Abaya and Chamo and dries out within ~40 years if the hydrological connection is cut off and the precipitation amount decreases to present day conditions. Several of such rapid lake level fluctuations, from a freshwater to a saline lake, might have occurred during the termination of the AHP, when humid conditions were less stable. Fast changes in fresh water availability requires high adaptability for humans living in the area and might have exerted severe environmental stress on humans in a sub-generational timescale.
Eastern African precipitation is characterized by large interannual to centennial variability (Lamb et al., 1998; Nicholson, 2000; Junginger et al., 2014; Liu et al., 2017), controlling water resource availability in the region. This is accentuated by the uneven spatial distribution of precipitation, from the lowlands (<500 m) to the high plateaus (>4,000 m), due to the extreme topography of the East African Rift System (EARS) (Nicholson, 2017). In recent decades a decline in precipitation has led to an increase in the intensity and frequency of drought periods (Funk et al., 2008; Viste, 2012; Rowell et al., 2015).
Water availability has always played a fundamental role in eastern Africa, and fluctuations in eastern Africa's hydroclimate in the past have had a dramatic impact on human societies. Changes in orbital insolation during the late Pleistocene to Middle Holocene caused an “African Humid Period” (AHP; ~15–5 ka BP; deMenocal et al., 2000; Barker et al., 2004), when increased precipitation triggered changes in vegetation (Lamb et al., 2004), dust dynamics (Zielhofer et al., 2017), lake expansion and increased river flow (Barker et al., 2004; Costa et al., 2014; Revel et al., 2014; Bloszies et al., 2015). In the wider region, these environmental changes during the AHP have been associated with multiple migration dynamics, such as the reoccupation of the Sahara ~10.5 ka BP (Kuper and Kröpelin, 2006; Drake et al., 2013; Larrasoaña et al., 2013; Manning and Timpson, 2014), cultural innovation (Lario et al., 1997) and overlapped with a demographic transition to the Neolithic age (Manning and Timpson, 2014; Honegger and Williams, 2015). An increase in arid conditions has been documented during the termination of the AHP ~5 ka BP in eastern Africa (e.g., deMenocal et al., 2000; Barker et al., 2004; Foerster et al., 2012; Junginger et al., 2014; Bloszies et al., 2015). This aridification trend coincides with societal adaptations, such as the introduction of pastoralism to eastern Africa, and the emergence of highly organized and complex state-level societies (demenocal, 1995; Brooks, 2006; Kuper and Kröpelin, 2006; Garcin et al., 2012; Foerster et al., 2015). Environmental changes likely provide both opportunities and constraints, due to variability in water resources and food availability, necessitating adaptation. To better understand potential links between cultural changes and climate variability, paleoclimate records can be used in conjunction with archaeological records to reconstruct major influences on modern civilizations.
In eastern Africa, however, there are still spatial and chronological gaps in paleoclimate records, which are largely non-quantitative and often regionally specific. Discrepancies particularly occur between records from marine archives, and lacustrine sites and also whether lacustrine sites are located within or outside the EARS. Existing records covering the AHP, report conflicting accounts of past climate variability in the region, with some studies suggesting a rapid onset and termination of the AHP (deMenocal et al., 2000; Tierney et al., 2013; Collins et al., 2017), while others indicate a rapid onset but more gradual termination and climate evolution in response to insolation forcing (Renssen et al., 2006; Asrat et al., 2007; Foerster et al., 2012; Junginger et al., 2014; Trauth et al., 2018). While lakes with a constant groundwater/geothermal supply were mostly stable (e.g., Lake Tilo on the Ethiopian Plateau, Telford and Lamb, 1999), river-fed EARS lakes, such as Chew Bahir or Ziway-Shala in Ethiopia (Gillespie et al., 1983; Foerster et al., 2012) or Lakes Elementeita-Nakuru (Richardson and Dussinger, 1986), paleo-Lake Suguta in Kenya (Junginger et al., 2014) or Lake Manyara (Bachofer et al., 2015) demonstrated greater variability in lake levels in response to short-term precipitation-evaporation changes (Figure 1).
Figure 1. Overview of the study area, lake levels and lake-balance modeling results from previous an current studies: (A) Elevation map of Africa and the positions of major air boundaries and convergence zones during the NH summer. Indian Summer Monsoon (ISM), West African Monsoon (WAM), the Inter Tropical Convergence Zone (ITCZ), and Congo Air Boundary (CAB). (B) Elevation map of eastern Africa with the study area, major lakes, study area (white box) and the CAB position in July. (C,D) Known paleo-lake depths within the EARS during the African Humid Period (AHP) and published paleo-precipitation rates estimates using different lake-balance approaches and associated references. (E) Transect of the EARS from north-east to south, with modern and paleo level lake extents (light and dark blue, respectively), the study area (red box), and major drainage direction (dotted line). North-South transect of the eastern branch of the EARS, with modern and paleo-lake extents (light and dark blue, respectively), the study area (red box), and major drainage direction (dotted line). The figure was modified after Junginger and Trauth (2013).
The first steps toward quantifying precipitation rates during the AHP and lake responses to a changing precipitation-evaporation budget were conducted for several basins in the EARS by applying Lake Balance Modeling (LBM) (Figure 1; Hastenrath and Kutzbach, 1983; Bergner et al., 2003; Dühnforth et al., 2006; Kniess, 2006; Borchardt and Trauth, 2012; Junginger and Trauth, 2013). LBM uses field observations, such as 14C dated shorelines and overflow sills high above modern lakes (indicative of higher moisture availability), to calculate paleo-lake volumes, levels and paleo-precipitation rates. Such attempts of quantifying environmental parameters using field observations (e.g., shorelines) or proxy data (e.g., leaf waxes, stable isotopes, dust, pollen) through modeling approaches are defined as so-called Proxy System Modeling (PSM), as presented by Evans et al. (2013) or Dee et al. (2018).
For the AHP, the reconstruction of environmental parameters, such as paleo-temperatures, is based for example on GDGT and Tex86 proxies from lake sediments (e.g., Verschuren, 2004; Berke et al., 2012a; Damsté et al., 2012; Loomis et al., 2012, 2014, 2015). Precipitation reconstructions are usually created using pollen (e.g., Lamb et al., 2004; Vincens et al., 2005; Umer et al., 2007; Rucina et al., 2009; Marchant et al., 2018), stalagmites (e.g., Baker et al., 2010), leaf waxes (e.g., Berke et al., 2012b; Costa et al., 2014), and the aforementioned LBMs. Due to the lack of the aforementioned proxy data in our study region, Chew Bahir in southern Ethiopia, an LBM is chosen to calculate paleo-precipitation rate changes.
Previously developed LBMs for the EARS, however, only looked at individual basins and therefore did not consider the hydrological connection via surface overflow of many rift basins (Figure 1). Consequently, the hydrological inflow budget of these interconnected lakes may have been underestimated, affecting the calculated paleo-precipitation estimates from the LBMs. Another aspect that was not considered in previous LBMs is change in seasonality. For the AHP it is hypothesized that an additional rainy season developed due to a shift of the Congo Air Boundary (CAB) during the northern hemisphere summer months, therefore tremendously impacting seasonality, environment and associated feedback mechanisms (e.g., Junginger and Trauth, 2013; Costa et al., 2014; Bloszies et al., 2015; Beck et al., 2019).
Here we present a comprehensive hydro-balance modeling approach that simulates firstly multiple rift lakes from the southern Ethiopian Rift (Figure 1; Abaya, Chamo, Chew Bahir). For this purpose, we developed a LBM for the region, which is a further development of the former models (e.g., Blodgett et al., 1997; Lenters and Cook, 1999). We have applied a combined catchment approach including surface and subsurface flows. Two approaches were used to estimate the evapotranspiration; a bulk transfer equation and the Surface Energy Balance Algorithm for Land (SEBAL). The paleo-precipitation required during the AHP to fill lakes to their previous high-stands will be calculated under different scenarios including the incorporation of an additional rainy season as previously suggested by multiple studies (e.g., Junginger et al., 2014; Bloszies et al., 2015; Beck et al., 2019) and the theoretical feedback due to vegetation changes. We furthermore provide lake regression and transgression times under abrupt and gradual changes, in order to demonstrate the sensitivity of the investigated lake systems to changes in precipitation. Based on the outcome of the model, we compare lake response times to abrupt and gradual precipitation changes to aridity proxy data of an 11 m-long sediment core from the margin of paleo-lake Chew Bahir, which was published by Foerster et al. (2015) to reconstruct the paleo-lake level evolution of Chew Bahir during the past 20 ka.
The Investigated Catchments
There are three adjoining North-South trending, endorheic basins at the southern end of the Main Ethiopian Rift (MER) which are part of the greater (4,000 km) EARS. These are the basins of Lake Abaya, Lake Chamo, and paleo-lake Chew Bahir, from north to south, respectively (Figure 1). Lake Abaya (6.4°N, 37.95°E, 7 m deep, and lake surface area of 1,081 km2) is the northern-most and highest elevated (1,176 m a.s.l.) basin with a catchment size of 16,203 km2 (Figure 2). The major rivers draining into Lake Abaya include the Gelana, Bilate, Gidabo, Hare, Baso, and Amesa Rivers (Tiruneh, 2005). Lake Chamo (5.8°N,37.6°E, ~1,109 m a.s.l., 14 m deep, 310 km2 lake surface area) to the south is separated from Lake Abaya by the Quaternary Tosa Sucha Volcano (Ebinger et al., 1993; Figure 2J). The Kulfo River to the west is the main tributary, while Sile and Wezeka rivers are smaller tributaries into Lake Chamo (Tiruneh, 2005). The Kulfo River is also connected to Lake Abaya during wet periods when lake levels are higher (Tiruneh, 2005). The catchment area of Lake Chamo is ~1,793 km2. Both lakes, Chamo and Abaya are bounded by the Gamo-Gidole horst to the west and the Amaro horst to the east reaching elevations over 3,000 m a.s.l. and consisting mainly of Cenozoic rift volcanics unconformably overlaying Precambrian basement rocks (Davidson, 1983; Ebinger et al., 1993). Paleo-lake Chew Bahir (4.1–6.3°N, 36.5–38.1°E, ~498 m a.s.l, Figures 1, 2) to the south is most of the year a 210 km2 saline mudflat and the lowest in elevation of the herein investigated basins. The Chew Bahir surface comprises deltaic ephemeral swamps at its northerly and north easterly reaches, where the perennial Weyto and Segan rivers drain the north-western and north-eastern sides of the catchment, respectively (Figure 2). The catchment (20,650 km2) is bound by the MER rift shoulders, in the east by the Teltele Plateau (up to 1,600 m a.s.l.) and in the west by the Hammar Range (up to 1,300 m a.s.l.). The elevation of the rift shoulders decreases in a southerly direction. The northern side of the catchment is bound by uplifted Proterozoic metamorphic rocks and Cenozoic rift volcanics (Davidson, 1983; Figures 1, 2), hydrologically separating it from the Chamo and Abaya catchments. To the south, Cenozoic rift volcanics from the foothills of the Hammar Range prevent a direct hydrological link with the Turkana basin, today. The dry southern-most part of the catchment borders Kenya.
Figure 2. Geographical features of the studied lake basins: (A) Modern catchment boundaries, monthly temperature means in °C and precipitation in mm/month (IRI, last accessed 11/2019), overflow sill locations (white-blue circles), and paleo-lake dimensions with modern/paleo-depths in numbers. (B–E) Modeling outputs of paleo-lake Chew Bahir‘s dimension during different lake filling steps from modern to +5, +20, and +45 m when reaching the overflow level. (F–J) Pictures from study sites by A. Junginger.
Several paleo-shorelines, wave cut notches and sediment characteristics in the Ethiopian Rift and plateau regions, such as Lake Abhé (Gasse and Street, 1978), Lake Ziway-Shala-Abiyata (Gillespie et al., 1983; Chalié and Gasse, 2002), Lake Dendi (Wagner et al., 2018), or Lake Ashenge (Marshall et al., 2009) show evidence of deep lakes during the AHP. During those high stands, these lakes were overflowing into rivers that connected to other lake basins. Awulachew (2006) reported that lakes Abaya and Chamo were previously connected via surface streams. Here, Lake Abaya connected via an old channel through the Kulfo river with Lake Chamo (Figure 2). Awulachew (2006) and Kassa (2015) report of evidence that during high stands, Lake Chamo overflowed into the river Metenafesha and into the Sermale stream, both being tributaries of the Segan river, which in turn feeds paleo-lake Chew Bahir. During a field campaign in November 2018, paleo-shoreline deposits (shell beds) at Chew Bahir could be identified at the same elevation as the overflow sill (Figure 2). As calculated in section Digital Elevation Model Analysis, Lake Abaya's overflow sill to the Chamo catchment is at 1,194 m a.s.l. (6°0′40.38″N, 37°34′50.45″E) 18 m above its present lake level. Lake Chamo reaches its overflow sill to the Chew Bahir catchment at 1,123 m a.s.l. (5°51′49.66″N, 37°38′31.78″E), which is 14 m above its modern lake surface. The Chew Bahir basin has an overflow sill at 543 m a.s.l. (4°13′30.58″N, 36°39′35.19″E) to Lake Turkana, 45 m above the basin floor.
Table 1. Moderate-Resolution Imaging Spectroradiometer (MODIS) products and bands used for SEBAL calculation.
Modern vegetation in the Chew Bahir catchment is typically sparse and follows an east-west gradient (Friis et al., 2010). The eastern part of the Chew Bahir basin is classified as Acacia-Commiphora woodland and bushland, whereas the western part is described as Combretum-Terminalia woodland and wooded grassland (Friis et al., 2010). The higher elevated parts of the study area, such as close to Arba Minch (Figure 2), comprise evergreen Afro-Montane forest and grassland complex, which transitions to moist evergreen Afro-Montane forest in the north. The vegetation in the even higher mountainous areas are dominated by Afro-alpine vegetation with an Ericaceous Belt (Friis et al., 2010).
Seasonal rainfall in Chew Bahir is associated with the annual passage of a tropical rain belt (often referred to as the Inter Tropical Convergence Zone, ITCZ), which is predominantly fueled by the Indian Ocean (Figure 1, Levin et al., 2009; Nicholson, 2017). Generally, the tropical rain belt migrates between 10° North and South of the equator, resulting in a bimodal precipitation pattern close to the equator and an unimodal pattern at its limits (Nicholson, 1996). Consequently, the lower elevations of the Chew Bahir catchment precipitation is bimodal, with the “Belg” rains from March to May (long rains) and the “short rains” in October–November (Figure 2). The highlands northwest of Chew Bahir, including the Abaya and Chamo catchments, experience an unimodal rainfall pattern with only one longer wet season from March to November (Segele and Lamb, 2005; Williams and Funk, 2011). This wet season is, fueled by two different sources, with the “Belg” rains lasting from March to May and the “Kiremt” rains from June to September. The Kiremt rains are hypothesized to originate from the Atlantic Ocean-derived southwestern humid Congo air stream (Nicholson, 1996; Camberlin, 1997; Lamb et al., 2000). The eastern limit of this air stream is marked by the Congo Air Boundary, which separates the humid air masses from the west from dry air masses in the east (Nicholson, 1996). Annual to decadal fluctuations in the intensity and strength of precipitation are possibly related to a direct response to sea-surface temperature variations in the Indian and Atlantic Oceans (Nicholson, 2017). Compared to other areas at similar latitudes in Africa, Ethiopia receives considerably more rainfall due to its position in the MER and the resulting orographic effects (Nicholson, 2017).
Materials and Methods
We developed a LBM to quantify the present day and paleo-water balances of paleo-lake Chew Bahir, Lake Chamo and Lake Abaya. The model was coded in R and is available at GitHub (https://github.com/MLFischer/Lake-Balance-Model, 03/2019). In summary, a Digital Elevation Model (DEM) was used to delineate the catchment boundaries, lake bathometry, lake volumes and overflow sills (section Digital Elevation Model Analysis). The annual modern water balances were calculated for Lakes Abaya, Chamo and paleo-lake Chew Bahir (section Lake Balance of Lakes Abaya, Chamo, and Paleo-Lake Chew Bahir). Evaporation was estimated for the paleo-lake Chew Bahir catchment using SEBAL [section Surface Energy Balance Algorithm for Land (SEBAL)]. After parametrizing the model based on present-day conditions, we modeled the most recent lake high-stands of Lake Chamo, Abaya and paleo-lake Chew Bahir during the AHP (section Paleo-Lake Modeling). In section Calculation of Amplifier Lake Characteristics, we calculated the Hypsometric Integral (HI) and Aridity Index (AI) for the paleo-lake Chew Bahir, Lake Abaya and Chamo catchments which are used to define the Amplifier Lake characteristic of lake basins in the EARS following Olaka et al. (2010). In section Lake Level Reconstruction of Paleo-Lake Chew Bahir, we apply the modeled lake responses of paleo-lake Chew Bahir to various changes is precipitation and thus lake-level scenarios to changes in a K-proxy record from a drill-core taken from the paleo-lake Chew Bahir basin (Foerster et al., 2012, 2015; Trauth et al., 2018) and reconstruct the paleo-lake level of the past 20 ka.
Digital Elevation Model Analysis
The catchment boundary delineation was based on the DEM resulting from the Shuttle Radar Topography Mission (SRTM C-band; http://www2.jpl.nasa.gov/srtm/). The flow directions, flow accumulations, catchment boundaries and hydro-junctions (maximum possible lake extent) were calculated with the workflow for deranged terrains from the ArcHydro toolbox of ArcGIS 10.3 (Merwade, 2012). The DEM of each catchment was processed using the freq function in the R raster-package to produce an elevation frequency table, which then is used to compute the basin hypsometry and paleo-lake bathymetry (R Core Team, 2019).
Lake Balance Model
Three processing steps were applied for calculating the modern day annual water balance as later summarized in Equation (3). The parameters used in this model are climate parameter based on interpolated global gridded climate data (New et al., 2002), remote sensing data (Friedl et al., 2010; Platnick et al., 2015) and evaporation parametrization approaches (Schmugge and André, 1991). The input parameters are annual averages and they are summarized in Table 1.
Table 2. Moderate-Resolution Imaging Spectroradiometer (MODIS) products and bands used for SEBAL calculation.
The first was to calculate surface temperatures by a radiation-based surface energy flux (Equation 1). The energy balance was calculated using Newton's method to determine the surface temperature. The energy fluxes were assumed to be in an equilibrium state and include: the incoming short-wave radiation (Rsw), the incoming long-wave radiation (Rld) with the surface emissivity (ε), the long-wave radiation (Rlu) emitted by the Earth's surface, the sensible heat flux (H) and the latent heat flux (L· ETa), with L being the latent heat of vaporization, and ETa the rate of actual evaporation (Blodgett et al., 1997; Lenters and Cook, 1999).
In the second step (Equation 2), ETa was calculated using the bulk transfer formula described by Brutsaert (1982), with the surface-drag coefficient (CD), the wind speed (U), the soil moisture availability (f), the gas constant for dry air (R), the surface (Ts) and air (Ta) temperatures, the relative humidity (rh), and the saturation vapor pressure (es(Ts)) as a function of surface and air temperature.
The resulting ETa on land and water was calculated with Equations (1) and (2) from the environmental parameters summarized in Table 1. A sensitivity analysis of ETa for paleo-lake Chew Bahir to major environmental parameters, such as air temperature (Ta), cloud cover, relative humidity (rh), and wind speed (U) was performed.
The third step involved calculating the annual lake water balance of lakes Abaya, Chamo and paleo-lake Chew Bahir using ETa on land (El) and water (Ew), the precipitation amount from the Tropical Rainfall Measuring Mission (TRMM, Bookhagen, in review) (Pbas) and the basin-wide subsurface outflow/inflow (Sbas) following Equation (3). A positive and negative water balance suggest an increase and decrease in lake water volume, respectively. The modern-day conditions of Lake Abaya, Lake Chamo and paleo-lake Chew Bahir were assumed to be in equilibrium, where ΔV = 0 as suggested in previous studies, such as Bookhagen et al. (2001); Bergner et al. (2003), and Junginger and Trauth (2013). The water surplus was assumed to be groundwater recharge in contrast to previous LBM studies, which used a fixed precipitation to groundwater ratio (Blodgett et al., 1997), a fix value (Bergner et al., 2003) or that the groundwater recharge rate is zero (e.g., Bookhagen et al., 2001; Dühnforth et al., 2006; Junginger and Trauth, 2013).
The SEBAL Method
Southern Ethiopia is a data-sparse region, with few continuously monitored weather stations, making quantification of the evaporation difficult on a catchment scale. Hence, we used a second approach called Surface Energy Balance Algorithm for Land (SEBAL; Bastiaanssen et al., 1998a,b). SEBAL calculates the Earth's surface instantaneous energy flux (Equation 4) at the moment of satellite overpass, to determine the latent heat flux as the residual (Bastiaanssen et al., 1998a,b).
In order to populate the SEBAL input parameters, satellite remote sensing data products of the Moderate-Resolution Imaging Spectroradiometer (MODIS) were identified and resampled to an uniform spatial resolution of 500 m (Table 2). To quantify monthly evaporation (Zhang et al., 2011; Borchardt and Trauth, 2012), the available MODIS datasets for an average reference year (2010 CE, Kiptala et al., 2013) were processed. The selection criterion applied was a clear sky coverage of higher than 90 percent. The Land Cover MCD12Q1 was further utilized to estimate the surface roughness length and the surface drag coefficient, following Wieringa (1992).
To calculate ETa based on the surface energy fluxes we employed the net radiation (Rn) equation (Equation 4). ETa was calculated using the soil heat flux G, which accounts for the amount of radiant energy released or absorbed at the soil surface per unit time, the sensible heat flux H and the latent heat flux λE, which resolves as residual of Equation (4):
Rn was estimated following the clear sky approach described in Allen et al. (2007) as a function of the day of the year, longitude, latitude, altitude, α, ε and atmospheric emissivity.
G was calculated with Equation (5) as a function of Rn, α and the Normalized Difference Vegetation Index (NDVI), which is a standardized proxy for vegetation density and greenness (Bastiaanssen, 2000).
H was determined using Equation (6), with the specific heat of air at a constant pressure (Cp), density of moist air (p), the vertical temperature gradient (dT), and the aerodynamic heat transfer resistance (rah) (Bastiaanssen et al., 1998a).
To solve Equation (6) we applied an approach by Zhang et al. (2011) that utilizes the wind field data by New et al. (2002). The frictional velocity and the aerodynamic resistance to heat transport were calculated for each pixel according to the wind-field and the logarithmic wind-profile. Both variables converge through an iterative atmospheric stability correction depending on the Monin-Obukov-Length (Bastiaanssen et al., 1998a,b). The resulting calculated dT leads to H (Bastiaanssen et al., 1998a,b). The evaporation fraction (Δ) was calculated to solve for daily evaporation (ET24) from the instantaneous heat fluxes (Bastiaanssen, 2000; Bastiaanssen et al., 2005) with Equation (7).
ET24 is a function of daily net radiation (Rn24) and Δ (Bastiaanssen, 2000). Rn24 was calculated following De Bruin and Stricker (2000) through the integrated sinusoidal daily illumination and the extra-terrestrial solar radiation.
The daily evaporation (ET24) was averaged each month to calculate the monthly evaporation. Missing data for months with only cloudy days were interpolated from the adjacent months. To calculate ETa, we applied a slight modification of the water balance by Thornthwaite and Mather (1955). We defined the modeled monthly evaporation as ETp. Following that, we calculated the catchment average ETp per month and compared them to the precipitation average from the TRMM dataset (Bookhagen, in review). Two assumptions have been considered when calculating the monthly water balance: (1) The total monthly ETa is smaller than the total monthly precipitation (P) and the available soil water (SW) in the catchment. (2) 50% of the monthly water surplus, if occurring, is stored in the soil for the next month (SW) and 50% of the water surplus results in groundwater recharge (GW). In months where ETp < P, we set ETa = ETp. However, in months where ETp > P the monthly ETa was calculated by multiplying ETp by the ratio of available water and ETp (Thornthwaite and Mather, 1955; Borchardt and Trauth, 2012). The ETa spatial distribution of humid months does not need to be corrected, as it remains the ETp distribution. During dry months ETa was further modified using the normalized Topographic Wetness Index (TWI; Kirkby and Beven, 1979) as shown in Equation (9). This spatial correction factor using the TWI reflects the orographic influence on the soil water distribution.
In a last step, we compared the SEBAL and the bulk transfer LBM derived ETa estimates on land, as well as the subsurface recharge rates for the Chew Bahir catchment to assess the two different approaches. Lake evaporation estimates for Chamo and Abaya were compared to measured pan evaporation data over the years 1985–2005 from the Arba Minch weather station (6°2′N, 37°33′E), which is located between lakes Abaya and Chamo (Belete, 2009).
Paleo-Lake Modeling Method
In order to assess the sensitivity of lake levels to changes in precipitation and the impact of hydrological connectivity between basins, we combined the LBM's of Lake Abaya, Lake Chamo and paleo-lake Chew Bahir. The annual water balance of each catchment (Equation 3) and the consequent lake volumes were calculated to derive the potential change of lake surface areas. Annual changes in the total lake surface area were considered at the end of each annually modeled time-step and the ratio of land and water were adjusted to reflect the new surface evaporation fluxes for each modeled year. Furthermore, at each annual time step the overflow threshold was considered and an additional drainage function was activated if the lake water level exceeded the height of the overflow sill that was calculated in the DEM analysis (section Digital Elevation Model Analysis). For example, if the surface water level for Lake Abaya equals the overflow threshold, the water surplus is accounted for as inflow into the Chamo catchment and added to the water balance. The TRMM derived precipitation rates were multiplied by the total area of each catchment to calculate the volume of the total water influx. The computed ETa over water was multiplied by the lake area to derive the evaporated water volume above each lake and the ETa over land was multiplied by the remaining catchment areas to calculate the total ETa.
Modeling Increased Precipitation Scenarios
We modeled various changes in precipitation over the Abaya, Chamo, and Chew Bahir catchments to estimate the increase in precipitation that was required for paleo-lake Chew Bahir to reach its overflow sill, the point where paleo-lake Chew Bahir drained into Lake Turkana during lake high-stands in the AHP (Foerster et al., 2012). Covering both decadal and centennial transitions of paleo-lake Chew Bahir (Trauth et al., 2018), we ran 500-years simulations of the LBM. The precipitation is increasing by 0.1% in each simulation, providing a total of 300 simulations with modeled precipitation increases from 0 to 30% of present day. This allowed us to calculate both the transition times for paleo-lake Chew Bahir; to go from “no lake” conditions to a “flooded basin” and the amount of increased precipitation that would be required. In the simulations, the exchange rates between the basins were calculated and the amount of the water volume influx from Lake Abaya to Lake Chamo and Lake Chamo to paleo-lake Chew Bahir were quantified as a proportion of its total water budget. Due to a lack of paleo-vegetation information and its impact on the actual evaporation and resulting precipitation in the study region for the AHP, we used the estimates of a 7–15% increased precipitation during the AHP in the Kenyan Rift published by Bergner et al. (2003). Those estimates are presented as vegetation threshold within our results and are based on the same parametrization approach as our LBM and are coming from a comparable landscape composition.
Modeling Seasonality Changes in Precipitation
Previous studies have suggested a change in the seasonality of precipitation during the AHP, due to a shift in the Congo Air Boundary, resulting in additional precipitation during July to September in equatorial regions in eastern Africa (e.g., Junginger, 2011; Junginger and Trauth, 2013; Costa et al., 2014; Junginger et al., 2014; Bloszies et al., 2015; Beck et al., 2019). This would have resulted in a change of the evaporative potential on land, compared to present day, due to an increased potential moisture source over what were previously dry months (ETp > P). We calculate the difference between ETa and ETp using SEBAL from July to September and add this evaporation amount to Chew Bahir's catchment actual evaporation rate on land. We omit the catchment of the lakes Abaya and Chamo from this procedure as they don't have this modern day lack (ETp < P) of evaporable water for these months.
Transition Times for Onset and Termination of AHP
The timing of precipitation change was also modeled using different scenarios to estimate the temporal reaction of paleo-lake Chew Bahir's volume to abrupt and gradual precipitation changes during the onset and termination of the AHP. The first scenario (S1) was an abrupt transition, simulated by an instantaneous increase in precipitation to the threshold (TS) amount required to force paleo-lake Chew Bahir to the overflow sill. However, this only represents the minimum precipitation required to fill Paleo-lake Chew Bahir to the overflow level, therefore we additionally considered two scenarios, in which the threshold could have been exceeded (S2: TS+5%, S3: TS+10%). We simulated a gradual increase in precipitation from modern day conditions to TS over 50, 100, 250, and 500 years. The termination of the AHP was also simulated, where the abrupt scenario saw an immediate transition from increased precipitation rates (TS) to modern day precipitation. The gradual offset was simulated by decreasing the precipitation amount linearly from TS to the modern day amount within 50, 100, 250, and 500 years.
Amplifier Lake Characteristics
In order to compare the basins of Lake Abaya, Lake Chamo and paleo-lake Chew Bahir with other lakes within the EARS, we calculated the hypsometric integral (HI; Pike and Wilson, 1971) following Equation (10).
We calculated the AI of each, individual catchment as the quotient precipitation (mm a−1) and potential evaporation (mm a−1, Thornthwaite, 1948) within each catchment. The term amplifier lake characterizes lakes in the EARS that react very sensitively to even moderate climate changes (Street-Perrott, 1985; Olaka et al., 2010; Trauth et al., 2010). We compared the HI and AI of paleo-lake Chew Bahir, lakes Abaya, and Chamo with other lake basins and classified amplifier lakes in the EARS based on Olaka et al. (2010).
Chew Bahir Lake Level Reconstruction
To translate changes indicated by climate proxies in the Chew Bahir drill core sediments into actual lake level changes, we used the established aridity proxy potassium (K) (Foerster et al., 2012, 2015, 2018; Trauth et al., 2018). High K concentrations in the sedimentary record have been documented to be strongly controlled by the hydrochemistry of the paleo-lake and porewaters. Such changes in salinity and alkalinity are typically caused by a drier climate and associated lake water evaporation (Foerster et al., 2018). For our study, we are using the most complete and high resolution record from Chew Bahir for the last 20 ka BP, a composite of pilot core CB-01 and CB-03, as discussed in Foerster et al. (2015). In our study, we compare the timing of K-changes with modeled response times of paleo-lake Chew Bahir's lake level, and interpret whether long-term and short-term fluctuations can be explained by abrupt or gradual changes in precipitation based on their duration and also K content.
DEM Analysis Overflow Regime
As summarized in Table 3, an additional 25.2 km3 of water would be required to initiate drainage at the overflow sill from Lake Abaya toward Lake Chamo, resulting in an increased lake surface area of 3.6%. To initiate overflow from Lake Chamo toward Chew Bahir, only an additional 5.2 km3 is required, which would increase the lake area by 5%. Paleo-lake Chew Bahir would require 83.2 km3 for the lake level to reach the overflow sill and initiate drainage into Lake Turkana. In this scenario, paleo-lake Chew Bahir would cover 12% of the total catchment area.
Lake Balance Model Results
Sensitivity tests of the major environmental parameters in the LBM including temperature, cloud cover, humidity and wind speed were calculated for the paleo-lake Chew Bahir catchment (Figure 3). For example, a modeled mean temperature rise of 1°C with an error below 1°C (New et al., 2002) would result in an increase of 4.7% (42 mm a−1) in ETa for the average land evaporation in the paleo-lake Chew Bahir catchment (Figure 3A). Cloud cover measurements from MODIS typically have an error below 5% (Platnick et al., 2015). An increase in average cloud cover by 5% would decrease ETa by ~4% (367 mm a−1, Figure 3B). The dataset of the relative humidity shows errors up to ±15% (New et al., 2002). Hence, an increase in the average relative humidity of 15% would decrease ETa on land by 8.7% (78 mm a−1, Figure 3C). The windspeed dataset has the highest errors and according to New et al. (2002) the maximum deviation for specific measures sites in eastern Africa ranges between 40 and 60%. For example, an increase in wind speed of 1 m s−1 (corresponding to a 45% deviation) would enhance ETa by 19.5% (174 mm a−1, Figure 3C).
Figure 3. The effects of environmental parameters on lake and land evaporation for the Chew Bahir basin: (A) air temperature, (B) cloud cover, (C) relative humidity, and (D) windspeed. Dashed lines indicate used input parameter and resulting annual actual evaporation of the LBM. Solid thin lines indicate maximum/minimum errors with resulting maxima/minima of the annual actual evaporation.
Using 43 days with a clear sky of the reference year, the calculated mean Rn was 459 W m−2 (σ = 54 W m−2), the average G was 75 W m−2 (σ = 8.6 W m−2), and H was 130 W m−2 (σ = 27.3 W m−2). Those fluxes result in Δ of 0.64 (σ = 0.074) and λE of 0.38 mm h−1 (σ = 0.07). RN24 was 137 W m−2 (σ = 7.07 W m−2), resulting in ET24 of 3.18 mm d−1 (σ = 0.33 mm d−1).
The annual water balance shows a bimodal precipitation (916 mm a−1) and ETa (847 mm a−1) and an unimodal ETp (see Figure 4B). The spatial distribution (see Figure 4A) of ETa ranges from 52 to 1,360 mm a−1. The southern part of the catchment has the lowest values, whereas the northern and elevated rift shoulders, with streams and adjacent wetlands, has the highest ETa due to greater water availability.
Figure 4. SEBAL results for paleo-lake Chew Bahir catchment: (A) Map of actual evaporation distribution throughout the catchment. (B) Monthly precipitation and monthly water balance of the Chew Bahir catchment.
Evaporation Estimate Comparison
Using the SEBAL approach, ETa was 847 mm a−1 for the paleo-lake Chew Bahir catchment, compared to 892 mm a−1 using the bulk transfer method, which is a relative difference of 5% (45 mm a−1). This leads to a shift in the closed water balance of the paleo-lake Chew Bahir catchment and a difference in the estimated catchment-wide groundwater recharge rates between 69 mm a−1 (SEBAL) and 24 mm a−1 (LBM). We calculated annual lake evaporation for Lake Abaya of 1,513 mm and for Lake Chamo of 1,550 mm a−1 using the bulk transfer Lake Balance Model. At the Arba Minch metrological station (6°2'N, 37°33'E), class A pan evaporation was measured from 1985 to 2005 AD, with an average of 2,191 mm a−1. This results, with a correction factor of 0.85, in a 20-years average annual lake evaporation of 1,862 mm a−1 as calculated by Belete (2009). These values differ from our model results by 16.7 to 18.7%, which could be related to a high uncertainty of the individually determined correction factor for pan evaporation data.
Combined Lake Balance Model and Paleo-Water Balances Results
Modeling Increased Precipitation Scenarios
The transition of paleo-lake Chew Bahir from “no lake” to “overflow conditions” at a lake surface level of 543 m a.s.l., 45 m above basin floor, was simulated over a period of 500 years. The transition occurs at a minimum precipitation increase at a TS of 6.5% compared to modern day 12-years TRMM average and occurs within 302 years (Figure 5A). The transition time required to go from “no lake” to “overflow conditions” decreases with increasing precipitation to, for example, 39 years (TS+5% = 11.5%) and 21 years (TS+10% = 16.5%), as indicated by stars in Figure 5A.
Figure 5. The effects of precipitation increase on the investigated lake basins: (A) Transition times of paleo-lake Chew Bahir from “no lake” to “overflowing lake” into Lake Turkana under minimum threshold (TS) conditions where seasonality has not changed; under the consideration of an additional rainy season caused by a shift of the location of the Congo Air Boundary (TSCAB); and the incorporation of theoretically vegetation feedback (TSVEG) with values taken from Bergner et al. (2003), since no local paleo-data are available. The graph shows the increasing pace of lake level rise with increasing precipitation. (B) Water mass exchange between the three investigated basins for the scenarios, TS and TSCAB. (C) Precipitation and overflowing Lake Chamo-Abaya as relative water sources of paleo-lake Chew Bahir showing the importance of an enlarged catchment of Chew Bahir during overflow times of all lakes. Whereas, under TS conditions Chew Bahir inflow ratio is at 6.3%, the inflow doubles up to 12.7% when considering an additional rainy season.
The water flux from one basin to the next basin increases linearly as soon as the overflow threshold of each basin is reached (Figure 5B). For example, the water flux from Lake Chamo to paleo-lake Chew Bahir increases starting from zero when precipitation is enhanced by 1.1%, and rises up to 3.9 km3 a−1 when precipitation is 16.5% higher than today.
The calculated inflow ratio for paleo-lake Chew Bahir describes the ratio between the annual water derived by precipitation within the catchment and the inflow from Lake Chamo (Figure 5C). The inflow ratio quantifies the high importance of the extended catchment of Lake Abaya and Lake Chamo for the water balance of paleo-lake Chew Bahir. The inflow ratio of Chew Bahir increases starting from 0 to 1.1% enhanced precipitation and reaches 6.3% when precipitation enhances to the threshold precipitation (TS = 6.5%) for overflow conditions.
Seasonality Changes Results
The difference in annual evaporation on land of 167.6 mm for the paleo-lake Chew Bahir catchment that we assessed needs to be considered to compensate for the additional rainy months from June to September according to the SEBAL processing. This difference in annual evaporation on land results in an increased threshold to achieve the transition from “no lake” to “overflow conditions.” When considering a third rainy season to have been active during the AHP, the threshold precipitation to fill paleo-lake Chew Bahir up to its overflow level would have increased by 6.9% (TSCAB = TS+6.9% = 13.4%) compared to the modern day precipitation amount (Figures 5A,C). The consideration of a changing vegetation on the evaporation and the resulting necessity of an increased precipitation is 7–15%, as suggested by Bergner et al. (2003). This leads to a combined threshold TSVEG of 20.6 to 28.6% (TSVEG = TS+TSCAB+ 7–15% = 20.6–28.6%; 20–30% rounded) compared to the modern day annual precipitation (Figures 5A,C).
AHP Transition Times Results
For abrupt “lake fill” scenarios (Figure 6A), the lake filling curve using the minimum precipitation amount of 13.4% (TSCAB) required to fill paleo-lake Chew Bahir up its overflow sill, leads to a slow, around 300 years, asymptotic reaction of the lake surface level over time as soon as paleo-lake Chew Bahir derives a constant water flux from Lake Chamo. Whereas, an increase of the precipitation amount for example by 5% (TSCAB+5%) or 10% (TSCAB+10%) above the threshold (TSCAB) leads to a fast transition of paleo-lake Chew Bahir from “no lake” to “overflow conditions” within 30 and 22 years, respectively. Lake Abaya and Lake Chamo, however, are showing for all simulated precipitation signals decadal to sub-decadal transition times.
Figure 6. Modeled pace of lake transitions under abrupt and gradual precipitation changes: (A) Abrupt increase in precipitation from modern values to TSCAB values (13.4%) as a minimum possible amount of rainfall to fill up the investigated lake basins. TSCAB +5% and TSCAB +10% are provided as examples for even higher precipitation that was provided during the AHP, for example. (B) Lake fill scenario under gradual increase in precipitation from modern values to TSCAB values over 50, 100, 250, and 500 years. (C) Abrupt decrease in precipitation from TSCAB to modern values. (D) Lake regression scenario under gradual decrease in precipitation from TSCAB to modern values over 50, 100, 250, and 500 years.
For gradual “lake fill” scenarios (Figure 6B) within 50, 100, 250, and 500 years from the modern day precipitation amount to the threshold, the same asymptotic lake surface level reaction pattern in time is identified, but it is delayed through the gradual precipitation signal. This gradual delay factor leads to corresponding centennial transition (>350 years) times. Lake Abaya and Lake Chamo show for all simulated gradual “lake fill” scenarios filling times between 40–160 and 50–180 years, respectively.
We simulated abrupt “lake drain” scenarios (Figure 6C) starting with lakes at their specific overflow sill: Lake Abaya at +18 m, Lake Chamo at +14 m, and paleo-lake Chew Bahir at +45 m above the modern day lake surface levels. For this scenario paleo-lake Chew Bahir would have been dried out within 45 years, whereas Lake Chamo and Lake Abaya would need more than 100 years to reach their new equilibrium state.
The simulation of gradual “lake drain” scenarios within 50, 100, 250, or 500 years to modern day conditions leads to delayed transition times of paleo-lake Chew Bahir from 70 up to 370 years before drying up (Figure 6D). Under gradual precipitation decrease Lake Abaya dries up between 170 and >500 years, and Lake Chamo takes the longest with 370 to >500 years.
Amplifier Lake Characteristic Results
Following the main principles for amplifier lake characterization (HI 0.23–0.3; Olaka et al., 2010; Figure 7), Chew Bahir (HI: 0.23, AI: 0.83) and Abaya (HI: 0.25, AI: 1.72) can be characterized as amplifier lakes. Chew Bahir is with an AI <1 characterized as a lake with slightly higher potential evaporation than precipitation, whereas Abaya is the wettest of all investigated basins with AI>1 (Figure 7). Lake Chamo (HI: 0.15, AI: 1.33) cannot be considered as an amplifier lake.
Figure 7. The characterization of amplifier lakes: Hypsometric integral (HI) and Aridity Index (AI) of major lake basins within the EARS as published by Olaka et al. (2010) including new calculated values of paleo-lake Chew Bahir, Lake Abaya, and Lake Chamo (in yellow) from this study. Amplification of precipitation changes occurs in lakes with HI of 0.23–0.3 according to Olaka et al. (2010). We postulate that AI has changed in the course of climatic changes and thus amplification influences, such as groundwater and precipitation.
Advantages and Limitation of the LBM
The results from the hydro-balance modeling show that a surplus in precipitation of 20–30% compared to present day values is required to fill the presently dry paleo-lake Chew Bahir to its overflow sill at 545 m a.s.l. and a depth of ~45 m, resulting in a ~2,500 km2 large freshwater lake. These final estimates are the result of three different assumptions: (1) TS of +6.5% accounts for the minimum precipitation required to fill paleo-lake Chew Bahir to its overflow level (classical hydro-balance approach), (2) TSCAB equals the sum of TS+6.9% to include extra precipitation outside the current rainy season from the theoretical shift of the CAB over the rift region during the AHP as suggested by previous studies (e.g., Junginger, 2011; Costa et al., 2014; Junginger et al., 2014; Bloszies et al., 2015; Beck et al., 2019), and (3) TSVEG equals the sum TS+TSCAB +7–15% precipitation to include the impact of vegetation feedback (Figure 5). Vegetation feedback, which increases transpiration and thus precipitation was approximated in our model, based on biosphere-feedback modeling results (7–15%; Bergner et al., 2003) from the nearby Lake Naivasha basin in the central Kenya rift, due to a lack of direct paleo-vegetation data. Numerous studies from eastern Africa strongly suggest a pronounced vegetation change between the AHP and today (e.g., Lamb et al., 2004; Vincens et al., 2005; Umer et al., 2007; Rucina et al., 2009; Marchant et al., 2018) and that vegetation feedback played an important role in precipitation (e.g., Claussen et al., 2017). However, local and regional information during the AHP (e.g., pollen, non-pollen palynomorphs, charcoal or phytoliths from sediments) are currently lacking and would help refine estimates of paleo-precipitation changes due to vegetation feedback, allowing for improved estimations of past hydrological budgets in southern Ethiopia.
Our modeling results of precipitation increase (20–30%) for paleo-lake Chew Bahir are of the same order of magnitude as results for similar studies within the EARS during the AHP (Figure 1), such as Ziway–Shala (+28%, Gillespie et al., 1983), Lake Turkana (+20%, Hastenrath and Kutzbach, 1983), Suguta Valley (+26%, Junginger and Trauth, 2013), Lake Nakuru–Elmenteita (+23–45%%, Dühnforth et al., 2006; Kniess, 2006), and Lake Naivasha (+29–33%, Hastenrath and Kutzbach, 1983; Bergner et al., 2003). However, almost all results stem from different modeling approaches. None of the previous models have considered hydrological connectivity or the impact on evaporation (ETa) due to a third rainy season. The latter would have had especially affected the catchments of Lake Turkana, Lake Nakuru-Elementeita, Lake Naivasha and paleo-Lake Suguta, as they most likely experienced a precipitation increase due to the change in the position of the CAB during the AHP (Junginger et al., 2014). We therefore propose that lake catchments and LBM results in the EARS under the potential influence of the CAB during the AHP would require additional increased precipitation.
Another important aspect in hydrological modeling of eastern African rift lakes is the inclusion of basin surface connectivity in major humid periods (Figure 1). Paleo-lake Chew Bahir received 15–25% of its additional water during the AHP from overflowing lakes Abaya and Chamo via the Segan River. Such hydrological connectivity during the AHP almost doubled the catchment size of Chew Bahir from 20,650 to 38,647 km2 (Figure 2). However, hydrological basin connectivity was likely intermittent, only occurring during strong humid phases (Figure 8). Intermittent connectivity during high-stands could have been caused by short-term abrupt changes in the position of the CAB and associated temporal cut offs from the surplus of water. Such intermittent surficial hydrological connectivity is likely relevant for all EARS lakes in the mid and lower altitudes, such as lakes Abhé, Turkana, Suguta, Baringo-Bogoria, Magadi, and Manyara (Figure 1E).
Figure 8. Chew Bahir LBM-based paleo-lake level reconstruction and comparison to adjacent lakes: (A) CB01 K proxy-record by Foerster et al. (2015) with dry events discussed in the text and transition modes. (B) Reconstructed lake level curve of Chew Bahir based on the lake-balance model (LBM) applied to changes in the K concentration of CB01. (C) Human occupation history of the Ethiopian highlands implies enhanced activities when lakes in the lower elevated rift basins, such as paleo-lake Chew Bahir became to saline or even dried up. The Ethiopian highlands may thus have been a refuge area during unfavorable times as Foerster et al. (2015) suggested. Cow indicates cultural transition from fishing/hunting/gathering to herding in the Ethiopian Highlands. (D) Lake level reconstruction of Suguta Valley (Junginger et al., 2014) as the mirror basin of Chew Bahir. (E) Lake level reconstruction of Lake Turkana by Garcin et al. (2012) and Beck et al. (2019) (blue, shaded curve), and by Bloszies et al. (2015) (green). (F) Strontium isotopes measured on aquatic fossils from lake Turkana sediments (van der Lubbe et al., 2017).
The water balance has varied due to increased precipitation and changes in other environmental parameters (e.g., as mentioned in Table 1) which influence the inflows and outflows into the basin. It is difficult to quantify these parameters (e.g., cloud cover, wind speed, and relative humidity) during the AHP. Proxy record reconstructions over the AHP however, can provide some insights. The GDGT based temperature reconstruction for eastern Africa by Loomis et al. (2012) reveals a colder climate at the beginning of the AHP at around 15 ka BP and a climate which was up to 3°C between 12 and 5 ka BP. Based on our sensitivity analysis for environmental parameters (section Lake Balance Model), an increase of around three Kelvin could increase the actual evaporation by up to 14%. A decrease in temperature of three Kelvin could lower the evaporation by the same amount. In terms of the resulting precipitation estimates, the early AHP with a slightly colder climate would need a slightly lower precipitation amount to reach a positive water budget, whereas the late phase of the AHP would need a slight higher precipitation amount to reach a positive water budget.
In summary, our study reveals the same magnitude of precipitation increase (+20–30% of modern day) in eastern Africa for the AHP as most other studies, but separates this threshold into parts. Those parts would need to be added (precipitation during modern day dry months) or subtracted (paleo hydrological connectivity) for other lake basins within the EARS, but can only be quantified in a site specific manner.
Evaporation Estimates Discussion
We compensated the lack of ground measured climate data by using remote sensing MODIS products, modeled climate data and calculated the actual evaporation additionally through SEBAL. This algorithm as proposed by Bastiaanssen et al. (1998a,b) has been applied in numerous variations depending on the data, landscapes and research designs. Studies showed inaccuracy of SEBAL at a field-scale and suggested to use METRIC or SEBS instead (e.g., Losgedaragh and Rahimzadegan, 2018). Mkhwanazi et al. (2015) implemented advective conditions within the SEBAL for crop evaporation studies on a field-scale, because the algorithm underestimate actual evaporation with a mean bias Error of 17.1%. SEBAL for catchment-scale studies instead showed an accuracy above 95% (Bastiaanssen et al., 2005). A 3 years SEBAL application in Tanzania and Kenya with a comparable research design as our application achieved a difference between the SEBAL results and the water balance of the catchment of 12%. For a SEBAL application in the Suguta Valley the SEBAL ETa is 5% lower than the bulk transfer calculated ETa of the LBM (Borchardt and Trauth, 2012; Junginger and Trauth, 2013).
The evaporation estimates of the SEBAL, the bulk-transfer formula and the pan-evaporation measurement have been compared and vary between 5 and 18.7%. This means the input water fluxes for the LBM could vary at the same magnitude. But those fluxes are not independent, instead the modern day water budget for each year is in balance. An increase in ETaon land would necessarily imply, ETa on water or the groundwater recharge rate is overestimated. The same yields for the sensitivity analysis of the bulk-transfer derived ETa. The variation of the evaporation rate depending on the environmental parameters shows possible inter-annual variability and possible micro-climate caused differences within the catchment. This spatial variety is mapped by SEBAL. However, besides the sensitivity, spatial variety and inter-annual variability, the accuracy of each multi-annual average water flux rate is due to the multi-method verification and the closed water balance of each catchment with an definite error below 20%. The different subdivisions of precipitation within each catchment between ETa on water and land affects the paleo-precipitation threshold in approximately one-tenth of the maximum deviation between ETa on water and land. This leads to a robust error of the threshold of 2% and an overall paleo-precipitation threshold for paleo-lake Chew Bahir of plus 20–30% precipitation in a multi-annual average compared to the modern day conditions.
Subsurface Flows Discussion
Endorheic basins, and their subsequent lakes, can be classified into three groups; (1) flow-through lakes, (2) discharge lakes, and (3) recharge lakes (e.g., Olaka, 2011 and references therein). Flow-through lakes are balanced (Sbas = 0). Discharge lakes are positively affected by groundwater flowing into the lake (Sbas > 0) and recharge lakes lake water seeps into the groundwater (Sbas < 0). During the AHP, when the overall moisture availability in the atmosphere was much higher and lake levels rose, changes in the groundwater levels may have changed the classification and status of the basin. Changes in groundwater inflows and outflows, would have subsequent impacts on the overall lake water balance. For example, groundwater discharge into lake basins could reduce their sensitivity to reduced precipitation input in short-term dry periods and potentially delay lake desiccation in protracted dry phases. This is a possible mechanism to explain why extremely arid catchments, such as paleo-Lake Suguta or Lake Nakuru-Elementeita, developed and maintained paleo-lakes which were 10–300 m deeper than today (Olaka et al., 2010). For example, Olaka (2011) estimated response times of 2–2.7 ka for groundwater to drain a 41 km3 large aquifer underneath the Eburro volcano complex located between Lake Naivasha and Lake Nakuru in the central Kenya Rift (Figure 1). Unfortunately, the southern MER lack this kind of information, despite preliminary evidence of significant groundwater aquifers in the central part of the Ethiopian rift (Kebede and Travi, 2012).
In terms of the groundwater recharge rate and the subsequent effect on the paleo-precipitation threshold for the African Humid Period, we parameterised the groundwater recharge rate as the surplus of the water balance budget, which is in contrast to previous LBM studies (see section Lake Balance of Lakes Abaya, Chamo, and Paleo-Lake Chew Bahir). Based on Chernet (1993) recharge rates between 50 and 150 mm a−1 are typical, which would be ~5–15% of the annual precipitation within the Chew Bahir catchment. We calculated a groundwater recharge rate of 24 mm a−1 using the LBM and 69 mm a−1 using SEBAL. These comparisons show the possible errors of our modeling results due to this uncertainty, which are a quantity smaller than the annual fluxes in precipitation and evaporation (Dühnforth et al., 2006).
Lake Response Times Discussion
Amplifier Lake Characteristics
Olaka et al. (2010) postulated that amplifier lakes, where AI <1 (ETp > P), likely received greater groundwater or surface inputs during the AHP, in order to maintain expansive lake systems in semiarid to arid regions in the EARS. In contrast, lakes with AI > 1 (ETp < P), which are prevalent on the high plateaus of the rift margins and Ethiopian Domes (Figure 1), are and were during the AHP amplified through precipitation. However, these calculations are based on present day conditions and may be irrelevant in a paleo context where P and ETp would have been vastly different, and consequently the AI. We expect in a paleo context, that the position of the AI for the respective paleo lakes would likely shift vertically along the y-axes of Figure 7 over time periods of 10–10,000 years. Based on our LBM, for example, the rate of groundwater feeding paleo-lake Chew Bahir today remains smaller than the rate of groundwater discharge within the catchment and is furthermore a magnitude smaller than the water mass exchange rates due to precipitation, evaporation and surface basin connectivity. However, the groundwater system could act as a millennial-scale buffer, which is recharged during major humid periods may slow the transition to lake desiccation in drier periods (e.g., Garcin et al., 2009; Olaka, 2011).
Although Lake Chamo is not characterized as an amplifier lake, it can receive additional surface inflows from Lake Abaya, due to the very shallow overflow between the two lakes (Figures 1, 2). Due to the extremely shallow overflow sill of Lake Abaya toward Lake Chamo and toward paleo-lake Chew Bahir, we hypothesize that even sub-decadal precipitation shifts may have triggered overflow to the paleo-lake Chew Bahir basin, and thus would have caused abrupt changes in the water budget of paleo-lake Chew Bahir. It is therefore highly likely that paleo-lake Chew Bahir was able to change during the AHP from a fresh-water lake to a desert within only a couple of decades. Such sensitivity to even moderate climate changes makes the paleo-lake Chew Bahir basin responsive to even sub-decadal climate fluctuations in contrast to, for example the Suguta Valley, which has a greater buffering capacity to short-term variability in precipitation, due to its larger size and depth of the paleo-lake (paleo-lake depth of 295 m), and higher overflow sill (Figures 1, 8).
We understand paleo-lake Chew Bahir as a double amplifier lake, pronouncing the unique and high sensitivity of the landscape and lakes water-level reaction to precipitation signals. This is due to (1) the classification of its catchment as an amplifier lake (Figure 1) and (2) the importance of the catchment of lakes Chamo and Abaya for the paleo water balance of paleo-lake Chew Bahir. This expansion of paleo-lake Chew Bahirs catchment in case of over-spilling lakes Abaya and Chamo into the Chew Bahir basin amplifies the lakes reaction to precipitation changes and causes rapid filling if there is a significant water drain, or rapid draining if this hydrological connectivity gets cut off.
Pace and Magnitude of Lake Level Changes
In a previous study by Trauth et al. (2018) a change point analysis was applied to the aridity proxy K-record by Foerster et al. (2012, 2015) to investigate the timing and duration of climatic changes. Our study now is able to add the character and magnitude of lake level changes and thus allows to infer direct implications to biosphere impacts.
Trauth et al. (2018) stated that the most dramatic changes occurred during the onset of the AHP and the interruption by the Younger Dryas (YD), a global cold event, known as a pronounced dry episode in eastern Africa (e.g., Barker et al., 2004). The onset of the AHP in the Chew Bahir basin, so Trauth et al. (2018), occurred within ~240 years (15.7–15.46 ka BP). Within such period of time, paleo-lake Chew Bahir would have been able to reach the overflow sill toward Lake Turkana (Figures 1, 8). The extreme arid interval in Chew Bahir between 13.2 and 11.73 ka BP, largely coinciding with the YD chronozone, shows abrupt changes in the K- record. Those were presumably expressed in an abrupt lake level drop (within 45 years) and lake level rise (within 250 years) in Chew Bahir. After the return to full humid conditions, the AHP main phase was interrupted by four dry spells at 11–10.7, 9.8–9.0, 7.8–7.5 and 7.1–6.8 ka BP. Those dry events are also reflected in records from lakes Turkana and paleo-lake Suguta (Figure 8), Mount Kilimanjaro (Thompson et al., 2002), Arabia (Fleitmann et al., 2003), as well as marine records off the NE coast of Africa (e.g., Gupta et al., 2003). Most of these dry episodes have lasted around 300 years. We expect, from a modeling perspective, these dry spells to have been long enough to desiccate paleo-lake Chew Bahir completely.
Trauth et al. (2018) highlights fourteen 20–80 years lasting dry events within the K-proxy record during the overall gradual termination of the AHP. According to our LBM, decadal dry events of 40 years or longer were able to dry up the Chew Bahir basin completely. Dry periods of 40 years or shorter still caused a considerable reduction in lake level though not a complete desiccation (Figure 8). Over the entire length of the AHP termination, Chew Bahir appears to have been under severe desiccation pressure. Such rainfall variability was also recorded in stalagmites in the near-by south-eastern rift shoulders of Ethiopia during times of a progressively drier climate (Baker et al., 2010). However, short-term precipitation changes may have filled lakes Abaya and Chamo up to their shallow overflow sills and provided Chew Bahir every now and then with enough water to delay the drying up process. The short-term flickers at the end of the AHP suggest, that paleo-lake Chew Bahir had not enough time between those flickers to have reached the overflow sill again. This assumption is supported by the finding of van der Lubbe et al. (2017), who used strontium isotopes as a water provenance proxy over the AHP and found no indication of a contribution of Chew Bahir water after 6 ka into Lake Turkana (Figure 8F).
The recent dry period after the termination of the AHP, which has left the Chew Bahir basin dry for most of the time, is interrupted by brief excursions to wetter conditions at 3, 2.2, and 1.3 ka BP (Figure 8A). Such a lowering of K concentration in the sediments of paleo-lake Chew Bahir points toward a lake transgression and thus a positive water budget during this time. Based on the results of our LBM, we conclude that even though the duration of these events would have allowed a rise of Chew Bahir's paleolake level, it is unlikely that the lake has reached maximum lake depths during these short-term humid episodes. This conclusion is supported by two arguments: (1) We expect a different climatic forcing mechanism than during the AHP, where a precession minimum caused enhanced humidity and the shift of rain belts over the study region (e.g., Costa et al., 2014; Junginger et al., 2014; Bloszies et al., 2015; Beck et al., 2019). For the past ~5 ka, for example, short-term changes in sea surface temperatures and associated circulation systems could account for these short-term increases in precipitation (e.g., Nicholson, 2017; Bayon et al., 2019). (2) The second argument against maximum water levels in Chew Bahir is that K concentrations do not reach AHP minima. A clear signal whether our assumption, that lake Chew Bahir has not reached its overflow level, is true may only be proved by poxies that could reflect paelo-water provenance, such as strontium isotopes in aquatic fossils, as van der Lubbe et al. (2017) have shown. In summary, paleo-lake Chew Bahir was a large paleo-lake during the AHP that dried up abruptly at least 19 times within decades resulting in a desert-like rift floor similar to today. Such rapid environmental changes must have had a tremendous impact on the biosphere including humans.
Implications for Human Adaptations
A first attempt of comparing an available archaeological record from hypothesized refuge areas in the region with inferred phases of climatic stress from the Chew Bahir K-proxy record was made by Foerster et al. (2015, 2016) for the last 20 ka. The study used scarce but available radiocarbon frequencies of documented archaeological sites in the lush mountainous regions of SW Ethiopia as a possible indicator for changes in settlement activities in the highlands during dry periods in the lowlands, which are indicated by high K content in the drill core sediments of Chew Bahir (note inverse scale, Figure 8). Even though age model uncertainties, the indefinite incompleteness and natural biases in the archaeological record naturally constrain the possibility to directly correlate climatic and archaeological data sets (see Foerster et al., 2015, 2016) and references therein), the patterns found in the comparative study are a valuable starting point to indirectly infer shifts in human settlement activity, bearing the role of external factors in mind. The results tentatively suggest that both, long and short-term climatic change could have affected settlement patterns and cultural innovation differently, though the factor of human decision-making within environmental boundaries played an important but further incalculable role (Foerster et al., 2015).
On the one hand, short-term episodes of pronounced aridity in the lower elevated lake basins, such as Chew Bahir could have been a push-factor for a refugium-directed vertical movement of groups with highly mobile hunter-gatherers (Figure 8C). A now further specified variable is the climatic component in this comparison as shown in our study. The Southern Ethiopian Rift has responded sensitively to even shorter dry spells, so that living conditions would have deteriorated quickly when the rift floor became too dry. One of the reasons why hunter-gatherers might have returned to the South Ethiopian Rift region after such dry spells, could have been their dietary style at that time, that was mainly devoted to fishing (e.g., Owen et al., 1982; Hildebrand et al., 2018). The long-term transitions, on the other hand, driven by changes in orbital controlled insolation, could have fueled cultural adaptation and significant changes in the social organization within groups. Only after the gradual end of the AHP with the near-desiccation of almost all lakes in the Kenyan and Ethiopian rift over a longer period (>1,000 years), a cultural transition from hunter-gathering to pastoralism occurred (e.g., Marshall and Hildebrand, 2002; Garcin et al., 2012; Lesur et al., 2014). The change from fishing to herding seems to have been a dynamic process, since the introduction of cattle was catalyzed by the immigration of herders that were escaping the progressively drying Sahara region, where pastoralism had been introduced much earlier (e.g., Kuper and Kröpelin, 2006; Hildebrand and Grillo, 2012). The dried up lake beds and former lake margins could have provided new land and grazing grounds for this new life style, which also made people independent from fishing (e.g., Garcin et al., 2012).
Trauth et al. (2010) introduced the theoretical concept of allopatric speciation through precession forced long-term (>10,000 years) population separation and remixing within the EARS due to appearing and disappearing rift lakes. Based on our study results, the Chew Bahir, Abaya, and Chamo basins can be excluded from this concept, due to their extreme sensitive amplifier lake characteristics. We expect, from a modeling perspective, that these lakes reacted on much shorter time scales (<100 a), suggesting that intensive adaption was required from humans living at the lakes margins (cultural buffering; Galway-Witham et al., 2019). One of those adaptive strategies could have been short-term migration, either vertically to the nearby more humid higher elevated grounds (Foerster et al., 2015; Ossendorf et al., 2019) or longitudinally and latitudinally to lake or river refugia that were not affected by the dry spells. Both would have contributed to a periodic cultural and genetic exchange, which is thought to be one of the key drivers in cultural innovation and, on longer time scales, evolution (Lahr and Foley, 1998; Scerri et al., 2018; Galway-Witham et al., 2019), with innovations being generally favored by the exchange of cultural information (Ackermann et al., 2016).
We developed a comprehensive Lake Balance Model (LBM) for the southern Ethiopian rift with focus on paleo-lake Chew Bahir and its catchment as a first contribution to better understand the potential response of paleo-lake Chew Bahir to precipitation changes. We conclude that the following paleo-precipitation estimates would have been necessary to fill paleo-lake Chew Bahir until its overflow level: (1) the classic LBM approach resulted in a 6.5% (TS) precipitation increase, (2) the addition of a third rainy season in the region from July to September would cause a 6.9% (TSCAB) increase, leading to a total of TS+TSCAB = 13.4%. (3) The inclusion of vegetation feedback of +7–15%, resulted in an additional precipitation, which sums up to a robust estimate of ~20–30% increased precipitation compared to the modern day amount for the paleo-lake Chew Bahir during the African Humid Period. We furthermore determined the amplifier lake characteristics of Lakes Abaya (HI = 0.25, AI = 1.72), Chamo (HI = 0.15, AI = 1.33) and Chew Bahir (HI = 0.23, AI = 0.83). We found that the sensitivity of lake levels is greatly increased during periods of high humidity due to overflowing lakes along the EARS axis and enhanced surface flow. Such additional water resources should be taken into account in the analysis of the paleo-hydrology of other lakes in the EARS when abrupt or gradual changes in the proxies are observed. Accordingly, we implemented all additional factors in our new LBM and were able to calculate lake level response times to abrupt and gradual precipitation changes and characterize lake level changes from K-proxy changes in a drill core from Chew Bahir. Based on our LBM results, we can now support the proposed abruptness of the onset of the AHP in southern Ethiopia by Trauth et al. (2018), similar to the beginning and end of the millennial-scale dry episode during the Younger Dryas. The results of our LBM furthermore support the gradual termination of the AHP, with the model indicating that the reported 20–80 years lasting dry events that are punctuating the termination could have been able to dry up the paleo-lake completely, when caused by abrupt precipitation changes.
Data Availability Statement
The hydro-balance model coded in R is available at GitHub (https://github.com/MLFischer/Lake-Balance-Model, 03/2019).
This study was designed and financed by AJ and MF. MF developed and performed the final LBM. MM, FB, VF, MT, and AJ contributed the data and analysis tools. All authors contributed to the article and approved the submitted version.
The research has been funded by the Ministry of Culture and Science (MWK) of Baden Württemberg, Germany, the Stiftung der deutschen Wirtschaft and the Open Access Publishing Fund of the University of Tübingen.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This study was conducted as part of the project Wet Feet or Walking on Sunshine of AJ funded by the Ministry of Culture and Science (MWK) of Baden Württemberg, Germany and by Stiftung der deutschen Wirtschaft financing the doctoral project of MF. We thank the three reviewers for their helpful comments, which improved the quality of this manuscript. We also acknowledge support by the Open Access Publishing Fund of the University of Tübingen. We thank Bodo Bookhagen from the University of Potsdam for his advises in using MODIS data. We were also grateful for valuable insights from long-term field experience in southern Ethiopia given by Frank Schäbitz from the University of Cologne and Henry Lamb from the University of Aberystwyth. Walter Düsing (University of Potsdam) and Fabian Sittaro (University of Leipzig): thank you for your inspiring discussions about paleo-data analysis, hydro-balance models, vegetation models, and life.
Allen, R. G., Tasumi, M., Morse, A., Trezza, R., Wright, J. L., Bastiaanssen, W., et al. (2007). Satellite-based energy balance for mapping evapotranspiration with internalized calibration (METRIC)—applications. J. Irrig. Drainage Eng. 133, 395–406. doi: 10.1061/(ASCE)0733-9437(2007)133:4(395)
Asrat, A., Baker, A., Mohammed, M. U., Leng, M. J., Calsteren, P. V., and Smith, C. (2007). A high-resolution multi-proxy stalagmite record from Mechara, Southeastern Ethiopia: palaeohydrological implications for speleothem palaeoclimate reconstruction. J. Quat. Sci. 22, 53–63. doi: 10.1002/jqs.1013
Awulachew, S. B. (2006). Investigation of physical and bathymetric characteristics of Lakes Abaya and chamo, Ethiopia, and their management implications. Lakes Reserv. Res. Manag. 11, 133–140. doi: 10.1111/j.1440-1770.2006.00300.x
Bachofer, F., Quénéhervé, G., Maerker, M., and Hochschild, V. (2015). Comparison of SVM and boosted regression trees for the delineation of lacustrine sediments using multispectral ASTER data and topographic indices in the Lake Manyara basin. Photogrammetrie-Fernerkundung-Geoinformation 2015, 81–94. doi: 10.1127/pfg/2015/0251
Baker, A., Asrat, A., Fairchild, I. J., Leng, M. J., Thomas, L., Widmann, M., et al. (2010). Decadal-scale rainfall variability in Ethiopia recorded in an annually laminated, holocene-age, stalagmite. Holocene 20, 827–836. doi: 10.1177/0959683610365934
Barker, P. A., Talbot, M. R., Street-Perrott, F. A., Marret, F., Scourse, J., and Odada, E. O. (2004). “Late Quaternary climatic variability in intertropical Africa,” in Past Climate Variability through Europe and Africa, Developments in Paleoenvironmental Research, vol. 6. eds R. W. Battarbee, F. Gasse, C. E. Stickley (Dordrecht: Springer), 117–138. doi: 10.1007/978-1-4020-2121-3_7
Bastiaanssen, W. G., Menenti, M., Feddes, R. A., and Holtslag, A. A. M. (1998a). A remote sensing surface energy balance algorithm for land (SEBAL). 1. Formulation. J. Hydrol. 212, 198–212. doi: 10.1016/S0022-1694(98)00253-4
Bastiaanssen, W. G., Pelgrum, H., Wang, J., Ma, Y., Moreno, J. F., Roerink, G. J., et al. (1998b). A remote sensing surface energy balance algorithm for land (SEBAL).: Part 2: validation. J. Hydrol. 212, 213–229. doi: 10.1016/S0022-1694(98)00254-6
Bastiaanssen, W. G. M., Noordman, E. J. M., Pelgrum, H., Davids, G., Thoreson, B. P., and Allen, R. G. (2005). SEBAL model with remotely sensed data to improve water-resources management under actual field conditions. J. Irrig. Drainage Eng. 131, 85–93. doi: 10.1061/(ASCE)0733-9437(2005)131:1(85)
Bayon, G., Schefuß, E., Dupont, L., Borges, A. V., Dennielou, B., Lambert, T., et al. (2019). The roles of climate and human land-use in the late Holocene rainforest crisis of Central Africa. Earth Planet. Sci. Lett. 505, 30–41. doi: 10.1016/j.epsl.2018.10.016
Beck, C. C., Feibel, C. S., Wright, J. D., and Mortlock, R. A. (2019). Onset of the African humid period by 13.9 kyr BP at kabua gorge, turkana Basin, Kenya. Holocene 29, 1011–1019. doi: 10.1177/0959683619831415
Belete, A. (2009). Climate change impact on lake abaya water level (M.Sc. thesis), Addis Ababa University, Ethiopia. Available online at: http://etd.aau.edu.~et/dspace/bitstream/123456789/2635/1/17730248189668230838464212617211563806
Bergner, A. G., Trauth, M. H., and Bookhagen, B. (2003). Paleoprecipitation estimates for the Lake Naivasha basin (Kenya) during the last 175 ky using a lake-balance model. Global Planet Change. 36, 117–136. doi: 10.1016/S0921-8181(02)00178-9
Berke, M. A., Johnson, T. C., Werne, J. P., Grice, K., Schouten, S., and Damsté, J. S. S. (2012b). Molecular records of climate variability and vegetation response since the late pleistocene in the lake victoria basin, East Africa. Quat. Sci. Rev. 55, 59–74. doi: 10.1016/j.quascirev.2012.08.014
Berke, M. A., Johnson, T. C., Werne, J. P., Schouten, S., and Damsté, J. S. S. (2012a). A mid-Holocene thermal maximum at the end of the African Humid Period. Earth Planet. Sci. Lett. 351–352, 95–104. doi: 10.1016/j.epsl.2012.07.008
Blodgett, T. A., Isacks, B. L., and Lenters, J. D. (1997). Constraints on the origin of paleolake expansions in the central andes. Earth Interact. 1, 1–28.doi: 10.1175/1087-3562(1997)001<0001:COTOOP>2.3.CO;2
Bloszies, C., Forman, S. L., and Wright, D. K. (2015). Water level history for lake turkana, Kenya in the past 15,000 years and a variable transition from the African humid period to holocene aridity. Glob. Planet. Change 132, 64–76. doi: 10.1016/j.gloplacha.2015.06.006
Bookhagen, B. (in review). High resolution spatiotemporal distribution of rainfall seasonality and extreme events based on a 12-year TRMM time series, in review. [Data set] Available online at: http://www.geog.ucsb.edu/~bodo/TRMM/ (accessed Agust 1, 2018).
Bookhagen, B., Haselton, K., and Trauth, M. H. (2001). Hydrological modelling of a pleistocene landslide-dammed lake in the santa maria basin, NW Argentina. Palaeogeogr. Palaeoclimatol. Palaeoecol. 169, 113–127. doi: 10.1016/S0031-0182(01)00221-8
Borchardt, S., and Trauth, M. H. (2012). Remotely-sensed evapotranspiration estimates for an improved hydrological modeling of the early holocene mega-lake suguta, northern Kenya Rift. Palaeogeogr. Palaeoclimatol. Palaeoecol. 361, 14–20. doi: 10.1016/j.palaeo.2012.07.009
Camberlin, P. (1997). Rainfall anomalies in the source region of the Nile and their connection with the Indian summer monsoon. J. Clim. 10, 1380–1392. doi: 10.1175/1520-0442(1997)010<1380:RAITSR>2.0.CO
Chalié, F., and Gasse, F. (2002). Late glacial–holocene diatom record of water chemistry and lake level change from the tropical east african rift lake abiyata (Ethiopia). Palaeogeogr. Palaeoclimatol. Palaeoecol. 187, 259–283. doi: 10.1016/S0031-0182(02)00480-7
Claussen, M., Dallmeyer, A., and Bader, J. (2017). “Theory and modeling of the African humid period and the green Sahara,” in Oxford Research Encyclopedia of Climate Science. doi: 10.1093/acrefore/9780190228620.013.532
Collins, J. A., Prange, M., Caley, T., Gimeno, L., Beckmann, B., Mulitza, S., et al. (2017). Rapid termination of the african humid period triggered by northern high-latitude cooling. Nat. Commun. 8:1372. doi: 10.1038/s41467-017-01454-y
Costa, K., Russell, J., Konecky, B., and Lamb, H. (2014). Isotopic reconstruction of the african humid period and congo air boundary migration at lake tana, Ethiopia. Quat. Sci. Rev. 83, 58–67. doi: 10.1016/j.quascirev.2013.10.031
Damsté, J. S. S., Ossebaar, J., Schouten, S., and Verschuren, D. (2012). Distribution of tetraether lipids in the 25-ka sedimentary record of Lake Challa: extracting reliable TEX86 and MBT/CBT palaeotemperatures from an equatorial African lake. Quat. Sci. Rev. 50, 43–54. doi: 10.1016/j.quascirev.2012.07.001
Dee, S. G., Russell, J. M., Morrill, C., Chen, Z., and Neary, A. (2018). PRYSM v2.0: a proxy system model for lacustrine archives. Paleoceanogr. Paleoclimatol. 33, 1250–1269. doi: 10.1029/2018PA003413
deMenocal, P., Ortiz, J., Guilderson, T., Adkins, J., Sarnthein, M., Baker, L., et al. (2000). Abrupt onset and termination of the African Humid Period: rapid climate responses to gradual insolation forcing. Quat. Sci. Rev. 19, 347–361. doi: 10.1016/S0277-3791(99)00081-5
Didan, K. (2015). MOD13Q1 MODIS/Terra Vegetation Indices 16-Day L3 Global 250m SIN Grid V006 [Data set]. NASA EOSDIS Land Processes DAAC. Available online at: https://doi.org/10.5067/MODIS/MOD13Q1.006 (accessed August 1, 2018).
Drake, N. A., Breeze, P., and Parker, A. (2013). Palaeoclimate in the saharan and arabian deserts during the middle palaeolithic and the potential for hominin dispersals. Quat. Int. 300, 48–61. doi: 10.1016/j.quaint.2012.12.018
Ebinger, C. J., Yemane, T., Woldegabriel, G., Aronson, J. L., and Walter, R. C. (1993). Late eocene–recent volcanism and faulting in the southern main Ethiopian rift. J. Geol. Soc. 150, 99–108. doi: 10.1144/gsjgs.150.1.0099
Evans, M. N., Tolwinski-Ward, S. E., Thompson, D. M., and Anchukaitis, K. J. (2013). Applications of proxy system modeling in high resolution paleoclimatology. Quat. Sci. Rev. 76, 16–28. doi: 10.1016/j.quascirev.2013.05.024
Fleitmann, D., Burns, S. J., Mudelsee, M., Neff, U., Kramers, J., Mangini, A., et al. (2003). Holocene forcing of the Indian monsoon recorded in a stalagmite from southern Oman. science 300, 1737–1739. doi: 10.1126/science.1083130
Foerster, V., Deocampo, D. M., Asrat, A., Günter, C., Junginger, A., Krämer, K. H., et al. (2018). Towards an understanding of climate proxy formation in the chew bahir basin, southern Ethiopian Rift. Palaeogeogr. Palaeoclimatol. Palaeoecol. 501, 111–123. doi: 10.1016/j.palaeo.2018.04.009
Foerster, V., Junginger, A., Langkamp, O., Gebru, T., Asrat, A., Umer, M., et al. (2012). Climatic change recorded in the sediments of the chew bahir basin, southern Ethiopia, during the last 45,000 years. Quat. Int. 274, 25–37. doi: 10.1016/j.quaint.2012.06.028
Foerster, V., Vogelsang, R., Junginger, A., Asrat, A., Lamb, H. F., Schaebitz, F., et al. (2015). Environmental change and human occupation of southern Ethiopia and northern Kenya during the last 20,000 years. Quat. Sci. Rev. 129, 333–340. doi: 10.1016/j.quascirev.2015.10.026
Foerster, V., Vogelsang, R., Junginger, A., Asrat, A., Lamb, H. F., Schaebitz, F., et al. (2016). Reply to the comment on “Environmental change and human occupation of southern Ethiopia and northern Kenya during the last 20,000 years. Quaternary Science Reviews 129: 333-340”. Quat. Sci. Revi. 141, 130–133. doi: 10.1016/j.quascirev.2016.04.003
Friedl, M. A., Sulla-Menashe, D., Tan, B., Schneider, A., Ramankutty, N., Sibley, A., et al. (2010). MODIS Collection 5 global land cover: algorithm refinements and characterization of new datasets. Remote Sens. Environ. 114, 168–182. doi: 10.1016/j.rse.2009.08.016
Funk, C., Dettinger, M. D., Michaelsen, J. C., Verdin, J. P., Brown, M. E., Barlow, M., et al. (2008). Warming of the Indian ocean threatens eastern and southern african food security but could be mitigated by agricultural development. Proc. Natl. Acad. Sci. U.S.A. 105, 11081–11086. doi: 10.1073/pnas.0708196105
Garcin, Y., Junginger, A., Melnick, D., Olago, D. O., Strecker, M. R., and Trauth, M. H. (2009). Late pleistocene–holocene rise and collapse of lake suguta, northern Kenya Rift. Quat. Sci. Rev. 28, 911–925. doi: 10.1016/j.quascirev.2008.12.006
Garcin, Y., Melnick, D., Strecker, M. R., Olago, D., and Tiercelin, J. J. (2012). East African mid-holocene wet–dry transition recorded in palaeo-shorelines of lake turkana, northern Kenya Rift. Earth Planet. Sci. Lett. 331, 322–334. doi: 10.1016/j.epsl.2012.03.016
Gasse, E., and Street, F. A. (1978). Late Quaternary lake-level fluctuations and environments of the northern Rift Valley and Afar region (Ethiopia and Djibouti). Palaeogeogr. Palaeoclimatol. Palaeoecol. 24, 279–325. doi: 10.1016/0031-0182(78)90011-1
Gupta, A. K., Anderson, D. M., and Overpeck, J. T. (2003). Abrupt changes in the Asian southwest monsoon during the Holocene and their links to the North Atlantic Ocean. Nature 421:354. doi: 10.1038/nature01340
Hildebrand, E. A., Grillo, K. M., Sawchuk, E. A., Pfeiffer, S. K., Conyers, L. B., Goldstein, S. T., et al. (2018). A monumental cemetery built by eastern Africa's first herders near Lake Turkana, Kenya. Proc. Natl. Acad. Sci. U.S.A. 115, 8942–8947. doi: 10.1073/pnas.1721975115
Honegger, M., and Williams, M. (2015). Human occupations and environmental changes in the Nile valley during the holocene: the case of kerma in upper nubia (northern Sudan). Quat. Sci. Rev. 130, 141–154. doi: 10.1016/j.quascirev.2015.06.031
Junginger, A., Roller, S., Olaka, L. A., and Trauth, M. H. (2014). The effects of solar irradiation changes on the migration of the congo air boundary and water levels of paleo-lake suguta, northern Kenya Rift, during the African humid period (15–5 ka BP). Palaeogeogr. Palaeoclimatol. Palaeoecol. 396, 1–16. doi: 10.1016/j.palaeo.2013.12.007
Junginger, A., and Trauth, M. H. (2013). Hydrological constraints of paleo-lake suguta in the northern kenya rift during the African humid period (15–5 ka BP). Glob. Planet. Change 111, 174–188. doi: 10.1016/j.gloplacha.2013.09.005
Kiptala, J. K., Mohamed, Y., Mul, M. L., and Van der Zaag, P. (2013). Mapping evapotranspiration trends using MODIS and SEBAL model in a data scarce and heterogeneous landscape in Eastern Africa. Water Resour. Res. 49, 8495–8510. doi: 10.1002/2013WR014240
Lahr, M., and Foley, R. A. (1998). Towards a theory of modern human origins: geography, demography, and diversity in recent human evolution. Am. J. Phys. Anthropol. 107, 137–176. doi: 10.1002/(SICI)1096-8644(1998)107:27+<137::AID-AJPA6>3.0.CO;2-Q
Lamb, A. L., Leng, M. J., Lamb, H. F., and Mohammed, M. U. (2000). A 9000-year oxygen and carbon isotope record of hydrological change in a small Ethiopian crater lake. Holocene 10, 167–177. doi: 10.1191/095968300677444611
Lamb, A. L., Leng, M. J., Mohammed, M. U., and Lamb, H. F. (2004). Holocene climate and vegetation change in the main ethiopian rift valley, inferred from the composition (C/N and δ13C) of lacustrine organic matter. Quat. Sci. Rev. 23, 881–891. doi: 10.1016/j.quascirev.2003.06.010
Lamb, P. J., Bell, M. A., and Finch, J. D. (1998). “Varibility of Sahelian disturbance lines and rainfall during 1951-87,” in Water Resources Variability in Africa during XXth Century, eds E. Servat et al. (IAHS), 19–26.
Lario, J., Sanchez-Moral, S., Fernandez, V., Jimeno, A., and Menendez, M. (1997). Palaeoenvironmental evolution of the blue Nile (Central Sudan) during the early and mid-holocene (Mesolithic-Neolithic transition). Quat. Sci. Rev. 16, 583–588. doi: 10.1016/S0277-3791(96)00053-4
Lenters, J. D., and Cook, K. H. (1999). Summertime precipitation variability over South America: Role of the large-scale circulation. Monthly Weather Rev. 127, 409–431. doi: 10.1175/1520-0493(1999)127<0409:SPVOSA>2.0.CO;2
Lesur, J., Hildebrand, E. A., Abawa, G., and Gutherz, X. (2014). The advent of herding in the horn of africa: new data from Ethiopia, Djibouti and Somaliland. Quat. Int. 343, 148–158. doi: 10.1016/j.quaint.2013.11.024
Levin, N. E., Zipser, E. J., and Cerling, T. E. (2009). Isotopic composition of waters from Ethiopia and Kenya: insights into moisture sources for eastern Africa. J. Geophys. Res. Atmospheres. 114:D23. doi: 10.1029/2009JD012166
Liu, C., Ikeda, K., Rasmussen, R., Barlage, M., Newman, A. J., Prein, A. F., et al. (2017). Continental-scale convection-permitting modeling of the current and future climate of North America. Clim. Dynam. 49, 71–95. doi: 10.1007/s00382-016-3327-9
Loomis, S. E., Russell, J. M., Eggermont, H., Verschuren, D., and Damsté, J. S. S. (2014). Effects of temperature, pH and nutrient concentration on branched GDGT distributions in East African lakes: Implications for paleoenvironmental reconstruction. Org. Geochem. 66, 25–37. doi: 10.1016/j.orggeochem.2013.10.012
Loomis, S. E., Russell, J. M., Ladd, B., Street-Perrott, F. A., and Damsté, J. S. S. (2012). Calibration and application of the branched GDGT temperature proxy on East African lake sediments. Earth Planet. Sci. Lett. 357, 277–288. doi: 10.1016/j.epsl.2012.09.031
Loomis, S. E., Russell, J. M., and Lamb, H. F. (2015). Northeast African temperature variability since the Late Pleistocene. Palaeogeogr. Palaeoclimatol. Palaeoecol. 423, 80–90. doi: 10.1016/j.palaeo.2015.02.005
Losgedaragh, S. Z., and Rahimzadegan, M. (2018). Evaluation of SEBS, SEBAL, and METRIC models in estimation of the evaporation from the freshwater lakes (Case study: Amirkabir dam, Iran). J. Hydrol. 561, 523–531. doi: 10.1016/j.jhydrol.2018.04.025
Marchant, R., Richer, S., Boles, O., Capitani, C., Courtney-Mustaphi, C. J., Lane, P., et al. (2018). Drivers and trajectories of land cover change in East Africa: Human and environmental interactions from 6000 years ago to present. Earth Sci. Rev. 178, 322–378. doi: 10.1016/j.earscirev.2017.12.010
Marshall, M. H., Lamb, H. F., Davies, S. J., Leng, M. J., Kubsa, Z., Umer, M., et al. (2009). Climatic change in northern Ethiopia during the past 17,000 years: a diatom and stable isotope record from Lake Ashenge. Palaeogeogr. Palaeoclimatol. Palaeoecol. 279, 114–127. doi: 10.1016/j.palaeo.2009.05.003
Merwade, V. (2012). Watershed and Stream Network DelineationUsing Arc Hydro Tools. Available online at: https://web.ics.purdue.edu/~vmerwade/education/terrain_processing.pdf (accessed August 2018).
Mkhwanazi, M., Chavez, J. L., and Andales, A. A. (2015). SEBAL-A: A remote sensing ET algorithm that accounts for advection with limited data. Part I: Development and validation. Remote Sens. 7, 15046–15067. doi: 10.3390/rs71115046
Ossendorf, G., Groos, A. R., Bromm, T., Tekelemariam, M. G., Glaser, B., Lesur, J., et al. (2019). Middle stone age foragers resided in high elevations of the glaciated bale mountains, Ethiopia. Science 365, 583–587. doi: 10.1126/science.aaw8942
Pike, R. J., and Wilson, S. E. (1971). Elevation-relief ratio, hypsometric integral, and geomorphic area-altitude analysis. Geol. Soc. Am. Bull. 82, 1079–1084. doi: 10.1130/0016-7606(1971)821079:ERHIAG2.0.CO
Platnick, S., Ackerman, S. A., King, M. D., Meyer, K., Menzel, W. P., Holz, R. E., et al. (2015). MODIS atmosphere L2 cloud product (06_L2), NASA MODIS Adaptive Processing System. Goddard Space Flight Center. doi: 10.5067/MODIS/MYD06_L2.006
R Core Team (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/
Renssen, H., Brovkin, V., Fichefet, T., and Goosse, H. (2006). Simulation of the holocene climate evolution in northern Africa: the termination of the African humid period. Quat. Int. 150, 95–102. doi: 10.1016/j.quaint.2005.01.001
Revel, M., Colin, C., Bernasconi, S., Combourieu-Nebout, N., Ducassou, E., Grousset, F. E., et al. (2014). 21,000 Years of Ethiopian African monsoon variability recorded in sediments of the western nile deep-sea fan. Reg. Environ. Change 14, 1685–1696. doi: 10.1007/s10113-014-0588-x
Rucina, S. M., Muiruri, V. M., Kinyanjui, R. N., McGuiness, K., and Marchant, R. (2009). Late quaternary vegetation and fire dynamics on mount Kenya. Palaeogeogr. Palaeoclimatol. Palaeoecol. 283, 1–14. doi: 10.1016/j.palaeo.2009.08.008
Scerri, E. M., Thomas, M. G., Manica, A., Gunz, P., Stock, J. T., Stringer, C., et al. (2018). Did our species evolve in subdivided populations across Africa, and why does it matter? Trends Ecol. Evol. 33, 582–594. doi: 10.1016/j.tree.2018.05.005
Telford, R. J., and Lamb, H. F. (1999). Groundwater-mediated response to Holocene climatic change recorded by the diatom stratigraphy of an Ethiopian crater lake. Quat. Res. 52, 63–75. doi: 10.1006/qres.1999.2034
Thompson, L. G., Mosley-Thompson, E., Davis, M. E., Henderson, K. A., Brecher, H. H., Zagorodnov, V. S., et al. (2002). Kilimanjaro ice core records: evidence of holocene climate change in tropical Africa. Science 298, 589–593. doi: 10.1126/science.1073198
Tierney, J. E., Smerdon, J. E., Anchukaitis, K. J., and Seager, R. (2013). Multidecadal variability in East African hydroclimate controlled by the Indian Ocean. Nature 493:389. doi: 10.1038/nature11785
Trauth, M. H., Foerster, V., Junginger, A., Asrat, A., Lamb, H. F., and Schaebitz, F. (2018). Abrupt or gradual? Change point analysis of the late pleistocene–holocene climate record from chew bahir, southern Ethiopia. Quat. Res. 90, 321–330. doi: 10.1017/qua.2018.30
Trauth, M. H., Maslin, M. A., Deino, A. L., Junginger, A., Lesoloyia, M., Odada, E. O., et al. (2010). Human evolution in a variable environment: the amplifier lakes of Eastern Africa. Quat. Sci. Rev. 29, 2981–2988. doi: 10.1016/j.quascirev.2010.07.007
Umer, M., Lamb, H. F., Bonnefille, R., Lézine, A. M., Tiercelin, J. J., Gibert, E., et al. (2007). Late pleistocene and holocene vegetation history of the Bale mountains, Ethiopia. Quat. Sci. Rev. 26, 2229–2246. doi: 10.1016/j.quascirev.2007.05.004
van der Lubbe, H. J. L., Krause-Nehring, J., Junginger, A., Garcin, Y., Joordens, J. C. A., Davies, G. R., et al. (2017). Gradual or abrupt? Changes in water source of Lake Turkana (Kenya) during the African Humid Period inferred from Sr isotope ratios. Quat. Sci. Rev. 174, 1–12. doi: 10.1016/j.quascirev.2017.08.010
Vermote, E., and Wolfe, R. (2015). “MOD09GA MODIS/terra surface reflectance daily L2G global 1kmand 500m SIN grid V006 [Data set],” in NASA EOSDIS Land Processes DAAC. Available online at: https://doi.org/10.5067/MODIS/MOD09GA.006 (accessed January 8, 2020).
Verschuren, D. (2004). “Decadal and century-scale climate variability in tropical Africa during the past 2000 years,” in Past Climate Variability through Europe and Africa. Developments in Paleoenvironmental Research, vol 6, eds R. W. Battarbee, F. Gasse, C. E. Stickley (Dordrecht: Springer), 139–158. doi: 10.1007/978-1-4020-2121-3_8
Vincens, A., Buchet, G., Williamson, D., and Taieb, M. (2005). A 23,000 yr pollen record from Lake Rukwa (8 S, SW Tanzania): new data on vegetation dynamics and climate in Central Eastern Africa. Rev. Palaeobot. Palynol. 137, 147–162. doi: 10.1016/j.revpalbo.2005.06.001
Wagner, B., Wennrich, V., Viehberg, F., Junginger, A., Kolvenbach, A., Rethemeyer, J., et al. (2018). Holocene rainfall runoff in the central Ethiopian highlands and evolution of the River Nile drainage system as revealed from a sediment record from Lake Dendi. Glob. Planet. Change 163, 29–43. doi: 10.1016/j.gloplacha.2018.02.003
Wan, Z., Hook, S., and Hulley, G. (2015). “MOD11A1 MODIS/Terra land surface temperature/emissivity daily L3 global 1km SIN grid V006 [Data set],” in NASA EOSDIS Land Processes DAAC. Available online at: https://doi.org/10.5067/MODIS/MOD11A1.006 (accessed January 8, 2020).
Williams, A. P., and Funk, C. (2011). A westward extension of the warm pool leads to a westward extension of the walker circulation, drying eastern Africa. Clim. Dyn. 37, 2417–2435. doi: 10.1007/s00382-010-0984-y
Zhang, X. C., Wu, J. W., Wu, H. Y., and Li, Y. (2011). Simplified SEBAL method for estimating vast areal evapotranspiration with MODIS data. Water Sci. Eng.4, 24–35. doi: 10.3882/j.issn.1674-2370.2011.01.003
Zielhofer, C., von Suchodoletz, H., Fletcher, W. J., Schneider, B., Dietze, E., Schlegel, M., et al. (2017). Millennial-scale fluctuations in Saharan dust supply across the decline of the African humid period. Quat. Sci. Rev. 171, 119–135. doi: 10.1016/j.quascirev.2017.07.010
Keywords: African humid period, precipitation changes, abrupt and gradual changes, Chew Bahir, Lake Abaya, Lake Chamo, human-environment interaction
Citation: Fischer ML, Markowska M, Bachofer F, Foerster VE, Asrat A, Zielhofer C, Trauth MH and Junginger A (2020) Determining the Pace and Magnitude of Lake Level Changes in Southern Ethiopia Over the Last 20,000 Years Using Lake Balance Modeling and SEBAL. Front. Earth Sci. 8:197. doi: 10.3389/feart.2020.00197
Received: 12 March 2020; Accepted: 18 May 2020;
Published: 30 June 2020.
Edited by:David K. Wright, University of Oslo, Norway
Reviewed by:Nick Andrew Drake, King's College London, United Kingdom
Seifu Admassu Tilahun, Bahir Dar University, Ethiopia
Sylvia G. Dee, Rice University, United States
Copyright © 2020 Fischer, Markowska, Bachofer, Foerster, Asrat, Zielhofer, Trauth and Junginger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Markus L. Fischer, firstname.lastname@example.org