Rainfall trends and variation in the Maasai Mara ecosystem and their implications for animal population and biodiversity dynamics

Abstract

Rainfall exerts a controlling influence on the availability and quality of vegetation and surface water for herbivores in African terrestrial ecosystems. We analyse temporal trends and variation in rainfall in the Maasai Mara ecosystem of East Africa and infer their implications for animal population and biodiversity dynamics. The data originated from 15 rain gauges in the Mara region (1965–2015) and one station in Narok Town (1913–2015), in Kenya’s Narok County. This is the first comprehensive and most detailed analysis of changes in rainfall in the region of its kind. Our results do not support the current predictions of the International Panel of Climate Change (IPCC) of very likely increases of rainfall over parts of Eastern Africa. The dry season rainfall component increased during 1935–2015 but annual rainfall decreased during 1962–2015 in Narok Town. Monthly rainfall was more stable and higher in the Mara than in Narok Town, likely because the Mara lies closer to the high-precipitation areas along the shores of Lake Victoria. Predominantly deterministic and persistent inter-annual cycles and extremely stable seasonal rainfall oscillations characterize rainfall in the Mara and Narok regions. The frequency of severe droughts increased and floods intensified in the Mara but droughts became less frequent and less severe in Narok Town. The timings of extreme droughts and floods coincided with significant periodicity in rainfall oscillations, implicating strong influences of global atmospheric and oceanic circulation patterns on regional rainfall variability. These changing rainfall patterns have implications for animal population dynamics. The increase in dry season rainfall during 1935–2015 possibly counterbalanced the impacts of resource scarcity generated by the declining annual rainfall during 1965–2015 in Narok Town. However, the increasing rainfall extremes in the Mara can be expected to create conditions conducive to outbreaks of infectious animal diseases and reduced vegetation quality for herbivores, particularly when droughts and floods persist over multiple years. The more extreme wet season rainfall may also alter herbivore space use, including migration patterns.

About the Authors:

Gundula S. Bartzke

Contributed equally to this work with: Gundula S. Bartzke, Joseph O. Ogutu

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Biostatistics Unit, Institute of Crop Science, University of Hohenheim, Stuttgart, Germany

ORCID logo http://orcid.org/0000-0002-1249-4998

Joseph O. Ogutu

Contributed equally to this work with: Gundula S. Bartzke, Joseph O. Ogutu

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Software, Supervision, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Biostatistics Unit, Institute of Crop Science, University of Hohenheim, Stuttgart, Germany

Sabyasachi Mukhopadhyay

Roles Data curation, Methodology, Writing – review & editing

¶‡ These authors also contributed equally to this work.

Affiliation: Biostatistics Unit, Institute of Crop Science, University of Hohenheim, Stuttgart, Germany

Devolent Mtui

Roles Data curation, Writing – review & editing

¶‡ These authors also contributed equally to this work.

Affiliation: Directorate of Research, Tanzania Wildlife Research Institute, Arusha, Tanzania

Holly T. Dublin

Roles Data curation, Funding acquisition, Writing – review & editing

Affiliation: Wasaa Conservation Centre, IUCN Eastern and Southern Africa Regional Office, Nairobi, Kenya

Hans-Peter Piepho

Roles Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

Affiliation: Biostatistics Unit, Institute of Crop Science, University of Hohenheim, Stuttgart, Germany

Introduction

A better understanding of rainfall dynamics is indispensable for developing biodiversity conservation measures likely to be effective under climate change [1]. Such understanding requires carefully verified observational data to ensure accuracy and reliability. This is especially pertinent for Africa where high-quality observational rainfall datasets with sufficiently high spatial and temporal resolutions are rare [2] and noteworthy discrepancies often exist between digital datasets and original weather records [3] even for the same weather stations [4].

Here, we use carefully verified station rainfall data for the Maasai Mara ecosystem to answer the following questions. (1) Are there temporal trends in the monthly, annual and seasonal rainfall components? (2) Are there shifts in rainfall seasonality? (3) What are the dominant cycle periods of oscillations in the rainfall components and are the periods changing? (4) Are severe droughts and floods becoming more frequent and severe and do they persist over multiple years? (5) How might the changing rainfall patterns affect animal population and biodiversity dynamics based on known responses of animal abundance, reproduction, survival, disease susceptibility and migration to rainfall?

Rainfall is the principal driver of the population dynamics of savanna herbivores [5,6] because it controls plant biomass production [7,8] and plant nutrient concentration [9], which affect herbivore birth [6] and survival [10] rates, susceptibility to predation [11] and, ultimately, biomass [12,13]. Not surprisingly, oscillatory dynamics in ungulate population size [5] and ungulate fecundity [14] are coupled with inter-annual and seasonal rainfall oscillations in African savannas, respectively.

Droughts can cause substantial herbivore mortality and often regulate population size [15]. For example, 75% of wildebeest (Connochaetes taurinus mearnsi (Burchell)) deaths in Serengeti were caused by undernutrition and rainfall was the most important factor determining food supply [16]. Concentrations of herbivores around water points during droughts [17] can elevate vegetation damage [18], and result in increased competition and predation [19].

Droughts and floods also facilitate infestations by parasites [20] and diseases such as the canine distemper virus outbreak among Serengeti lions following a severe drought in 1993 [21]. Excessive rainfall can adversely affect small herbivores that require high-quality forage through diluting plant nutrient concentration [22]. High rainfall also promotes fires because it increases fuel loads in grasslands [23]. But low rainfall can lead to more destructive and extensive fires in savanna woodlands and forests [24].

Rainfall distribution and seasonality principally drives animal migration [25] and dispersal [26] in savannas. Consequently, animals may alter both their migratory [25] and short-term movements rapidly in response to localised rainfall patterns. During low rainfall years, animals are forced to travel longer distances between water and foraging grounds, making their offspring more vulnerable to predation [27].

Temperatures have risen in recent decades in most parts of the eastern African region [2,28–30] but the contemporaneous changes in rainfall seem subtle and largely unpredictable [2]. Rainfall increased in parts of East Africa during 1951–2001 [31] and during 1979–2010 [2]. Thus, rainfall trends around Lake Victoria were predominantly positive over the 20th century [32]. Likewise, General Circulation Models project increasing rainfall [33–36], more intense wet seasons and less severe droughts for most of East Africa [37]. Climate models also project strengthening of the El Niño Southern Oscillation (ENSO) and more frequent occurrences of the positive phase of the Indian Ocean Dipole [38,39] as temperatures rise, although evidence for the strengthening of the ENSO phenomenon remains controversial [40,41]. Such conditions facilitate moisture export from the Indian Ocean towards East Africa by weakening westerly winds [42,43] and lead to more intense wet seasons and floods.

The above patterns are in contrast to the findings of other studies of temporal trends in East African rainfall, such as decreasing annual [30,33], wet season [44–46] and dry season rainfall [36] in recent decades. Concurrently, droughts became more severe during 1970–2006 in Eastern and Southern Africa [44]. La Niña events, which often follow extreme El Niño events [47], typically lead to severe drying in East Africa [48]. Rainfall also declines during negative phases of the Indian Ocean Dipole [49]. For example, the 2005–2006 East African drought was associated with both a strong negative Indian Ocean Dipole and La Niña-like conditions [50,51].

Contrasting rainfall trends have also been documented for the Mara-Serengeti ecosystem. Ritchie et al. [52] reported a decrease in the total annual and wet season rainfall during 1960–2001 but an increase in the dry season rainfall in the Serengeti during 1913–2001. In contrast, Ogutu et al. [53] reported a decline in the dry season rainfall in the Maasai Mara Reserve (Mara Reserve) during 1975–2003. However, dry season rainfall in Narok, a Kenyan town located 75 km north-east of the Mara Reserve, increased during 1940–2004 after a protracted drought during 1930–1939 [53]. These contrasting findings demonstrate considerable uncertainty inherent in trends and variation in the past and anticipated future rainfall scenarios [35].

Climate warming can change rainfall seasonality and cycle periods by modulating ocean-atmosphere circulations. Climate warming alters ocean temperatures, cloud and ice cover, leading to shifts in the movement of the Inter-Tropical Convergence Zone [54]. This belt of rising and convecting air masses driven by solar radiation is the causal agent for rainfall seasonality as it moves southwards from East Africa during the transition between the dry and the wet season (July-January) and northwards during the transition between the wet and the dry season (January-July) [55].

The dominant East African rainfall cycle periods ranging between about 2 to 12.5 years appear very variable in space and time [56–58]. Besides the ENSO phenomenon and the Indian Ocean Dipole, oscillations in Atlantic ocean temperatures can influence East African rainfall considerably [59] through its teleconnections (an influence occurring over large distances, typically thousands of kilometres) via the Indian Ocean [60] and the West African monsoon [61]. The highly variable local topography of the African Rift Valley may also contribute to temporal and spatial variability in rainfall [62].

The aim of this study was to quantify trends and variation in rainfall in the Maasai Mara ecosystem in East Africa as a background for understanding their past and possible future implications for animal population and biodiversity dynamics.

Materials and methods

Study area

The Maasai Mara ecosystem is situated along the international border between Kenya and Tanzania in equatorial East Africa (34.7° to 35.4° E, 1.2° to 1.7° S, Fig 1). The elevation in the Mara Reserve ranges from about 1,450 m to about 2,100 m above sea level (Fig 1). The Mara Reserve was established in 1961 and the adjoining Serengeti National Park in 1951 to protect the rainfall-driven migration of the numerous wildebeest, zebra (Equus quagga burchellii) and Thomson’s gazelle (Gazella thomsonii) [63]. Multiple buffer zones with various degrees of protection and land use types now surround the protected areas [9] (Fig 1). The Maasai Mara ecosystem of Kenya consists mostly of grasslands with the cover of shrubs and thorny bushes increasing towards its northern and eastern extremes [64]. The dominant vegetation type in the Serengeti is savanna with a mixture of grasses, shrubs and trees [63]. Riverine forests fringe various streams and drainage lines in the area [63] (Fig 1) but much of this forest has been lost.

[Figure omitted. See PDF.]

Fig 1. The Mara-Serengeti ecosystem (orange borders) straddling the international border (black line) between Tanzania and Kenya.

https://doi.org/10.1371/journal.pone.0202814.g001

The entire Mara-Serengeti ecosystem supports an extremely diverse and abundant community of herbivores [65] and carnivores [66]. Populations of many wildlife species are declining, concurrent with changes in climate, growing pressures from cultivation, livestock grazing and other anthropocentric influences that are driven by accelerating human population growth [30,63,67]. The major land use type in the Maasai Mara ecosystem is pastoralism [30]. Higher rainfall areas in the buffer zone more distant from the Mara Reserve are now largely converted into agriculture [64]. Settlements, extensive agriculture and sedentary livestock holdings characterize the western side of the Serengeti National Park today [63].

The climatic year in the Mara starts in November and ends in October of the following year. The seasonal rainfall distribution is strongly bimodal, with the wet season spanning November-June and the dry season spanning July-October. The wet season consists of the short (November-December) and the long (January-June) rains. January-February trends to be dry and hence is sometimes called the short dry season in contrast to the long dry season from July-October [68]. The annual rainfall in the Mara reserve follows a spatial continuum from about 650 mm in the southeast to about 1300 mm in the northwest [68]. The dry season rainfall in particular is higher and more stable close to Lake Victoria ([55], Fig 1), the largest lake in Africa. Climate models suggest that the difference in air temperature between the land and the lake water creates a local convergence zone [69] that interferes with large-scale atmospheric and oceanic circulation patterns [70] and brings rain to the north-western extensions of the Mara-Serengeti ecosystem [55]. Consequently, rainfall decreases away from Lake Victoria, reaching about 750 mm in Narok Town in Kenya and 350–450 mm on the south-eastern Serengeti plains in Tanzania [63].

During the dry season, the seasonal watercourses dry out and the remaining isolated pools, ponds, springs and the Mara River, the only permanent river in the Mara-Serengeti ecosystem, serve as the only sources of water for wildlife and livestock across the area [71]. As water levels drop, the water quality of these remaining sources becomes very poor, in part because wildlife concentrate near water as the dry season progresses [71]. Forage quality and quantity progressively decline with time after the wet season [72]. Consequently, wildebeest migrate back from the south-eastern Serengeti Plains to the Mara in Kenya where the dry season rainfall is higher [55]. Poor rainfall can force wildebeest to leave these nutrient-rich grass plains in Serengeti earlier than usual but increased rainfall can cause them to return earlier [25].

Direct correlations between rainfall and ENSO appear weak in the Mara-Serengeti ecosystem [53,73] but major El Niño events can cause substantial mortality of herbivores [74]. Similarly, La Niña events [50,51] are sometimes followed by marked reductions in herbivore biomass [75]. The dominant rainfall cycles in the Mara-Serengeti ecosystem have periods of about 5 to 10 years but cycle lengths can vary widely [76].

In addition to rainfall variability, water flow in the Mara River has been declining as a consequence of upstream deforestation of the Mau forest and excessive water abstraction for irrigation in Kenya [77].

Data sources and processing

We obtained total monthly rainfall data for 15 gauges in the Mara spanning 1965–2014 (S1 Data) and for one gauge for Narok Town spanning 1913–2015 (S2 Data) from the sources listed in Table 1. Thirteen of the 15 gauges in the Mara were operated by the World Wide Fund for Nature (WWF) and Friends of Conservation (FOC) as part of the Maasai Mara Ecological Monitoring Programme from 1989 to 2003. The records for 5 of the 15 gauges were taken daily and then summed to obtain monthly totals. Rainfall for the other gauges was measured at monthly intervals.

[Figure omitted. See PDF.]

Table 1. The sources of rainfall records for 15 rain gauges in the Mara and a rain gauge in Narok Town in Kenya.

https://doi.org/10.1371/journal.pone.0202814.t001

Due to some obvious inconsistencies and discrepancies in the rainfall data sets supplied by government institutions, administrative bodies, research institutions and individual researchers [4], the Kenya Meteorological Department verified 93% of the monthly values for Narok Town Meteorological Station and 97% of the monthly rainfall values for Keekorok Hydromet Station availed to us against the original records that were handwritten either on data cards or in ledger books. When data cards were not available, we used the values recorded in the ledger books or other digital files. Thus we were able to compile, in our estimation, the most accurate historical rainfall dataset available for the study region for 1913–2015. For periods when no records were available for any of the gauges in the Mara (e.g. March 1968 to April 1969) we imputed total monthly rainfall values using the state space method of Piepho and Ogutu [78]. Only one missing monthly record (January 2014) was imputed for Narok Town.

We derived a single time series of rainfall across all available gauges in the Mara after adjusting monthly rainfall by the mean monthly rainfall at individual gauges (Eq 3 in S1 Text). This method performed best among four different standardisation approaches considered for accounting for spatial variation in rainfall in the Mara region (S1 Text). The standardisation methods (S1 Fig) were assessed in terms of how well they accounted for the level shift in rainfall due to the closure of several recording stations after 2003 and whether they reproduced the well-documented incidences of extreme events in the Mara-Serengeti ecosystem such as the extreme floods of 1998, the extreme droughts of 1982 and 2006, and the severe droughts of 1984, 1993, 1999, 2000 and 2009 [29,75,79–82].

We summed up the dry season rainfall spanning the months July-October in the Mara (Fig 2A) and June-October in Narok Town (Fig 2B) based on the distribution of monthly rainfall. In June, the rainfall in Narok Town is already much lower than rainfall for the wet season months unlike for the Mara, where migratory wildebeest normally arrive in July to occupy their dry season range [25]. The wet season rainfall spans two consecutive calendar years; November of the preceding year up to June of the current year in the Mara but up to May in Narok Town.

[Figure omitted. See PDF.]

Fig 2. The distribution of total monthly rainfall across months.

(A) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1965–2015. (B) Rainfall in Narok Town in Kenya was recorded during 1913–2015.

https://doi.org/10.1371/journal.pone.0202814.g002

Statistical analyses

Temporal trend and seasonality in monthly rainfall.

We used the unobserved components model (UCM, [83]), a special case of the linear Gaussian state space model, that is well suited to simultaneously analyse changes in trends () and seasonality (γn) in monthly rainfall levels (xn) while taking account of autocorrelation () by decomposing the time series for the n-th months as follows:(1)

This model assumes statistical independence of the different model components.

We first assume a random walk model for the time trend, or equivalently, that the trend remains approximately constant through time. The random walk trend model can be specified for the n-th month as(2)in which are independent and identically distributed (i.i.d.) normal errors or disturbances having zero mean and variance , assuming that is a Gaussian white noise process. It is noteworthy that = 0 implies that is a constant. A significant disturbance (error) variance implies that the level component is stochastic.

Besides the random walk model (2), we modelled the trend component using a locally linear time trend incorporating level and slope components and specified by(3)(4)where the disturbances and are assumed to be independent. We dropped the locally linear time trend when it was non-significant. The seasonal oscillations were modelled with trigonometric components (S2 Text).

Temporal trend and variation in the annual and seasonal rainfall components.

Analogous to the time series of the total monthly rainfall records we used the unobserved components model to estimate the overall annual, as well as the wet and dry season rainfall trends (μt) and cycles (φt) simultaneously by decomposing the time series (rt):(5)in which ϵt are i.i.d. normally distributed errors or disturbances having zero mean and variance σ2ϵ so that ϵt is a Gaussian white noise process. The trend component μt is defined analogously to Eq 2, where ηt are now the i.i.d. normally distributed stochastic disturbances.

Oscillations in annual and seasonal rainfall components.

In addition to the trend components we estimated the periods (p > 2) and damping factors (ρ) of the stochastic cycle components (φt) with a time-varying amplitude and phase given by(6)where 0 < ρ ≤ 1, λ = 2×π/p is the angular frequency of the cycle with 0 < λ < π, υt and υ*t are independent Gaussian disturbances with zero mean and variance σ2υ. Values of ρ, p and σ2υ are estimated from the data alongside the other model parameters. Significant stochastic cycle disturbance variances (σ2υ) imply stochastic and transient cycles.

The damping factor ρ governs the stationarity properties of the random sequence φt such that φt has a stationary distribution with mean zero and variance σ2υ/(1 ‒ ρ2) if ρ < 1 but is nonstationary if ρ = 1. A damping factor close to one indicates stable and persistent cycles. We specified and tested for significance of up to two cycles in the annual, wet season and dry season rainfall components.

The UCM models (1) and (6) were fitted by the diffuse Kalman filtering and smoothing algorithm [84] in the SAS UCM procedure [83].

Analysis of changing periodicity in monthly rainfall oscillations.

To reveal temporal changes in the periodicity of monthly rainfall oscillations that are hard to expose using UCM models, we applied wavelet analysis techniques using the R-package biwavelet Version 0.20.10 [85]. For the wavelet analysis, we de-seasonalized the monthly rainfall values by subtracting from each month the average of the total monthly rainfall values for that month. The de-seasonalized values were subsequently transformed with the Morlet wavelet function [86]:(7)using Eq 1 in S3 Text where κ is a nondimensional “time” parameter and ω0 is the nondimensional frequency. Significant oscillations in the time series of monthly rainfall were derived using normalized wavelet power spectra (Eq 2 in S3 Text) and plotted.

Frequency, severity and timing of severe droughts and floods.

We classified drought and flood years or seasons as severe or extreme. We define a year as a severe drought or flood year if the annual or seasonal rainfall component for the year (season) does not reach the estimated 10-year return level, implying that such low rainfall occurs with a probability of less than 0.1 Rainfall does not reach the 20-year return level estimated from the low rainfall values during extreme droughts. We use the term floods to denote very wet years or seasons. A severe flood year is one in which rainfall exceeds the 10-year extreme return level for the annual or seasonal component. However, rainfall exceeds the 20-year return level estimated from the high rainfall values during extreme flood years or seasons. This classification characterizes the extent of rainfall deficit or surfeit and portrays the associated broad transitions in rainfall influences on vegetation production and quality for herbivores.

To find out if droughts and floods are becoming more frequent, we first estimated return levels of droughts and floods for the annual, wet season and dry season rainfall components by applying extreme value analysis using the threshold excess approach in the R-package extRemes [87]. The return level is the estimated amount of rainfall that is not reached or exceeded once within a specified return period. Return levels are comparable to quantiles but are estimated statistically rather than empirically. The generalized Pareto distribution was used to approximate the distribution of extremely low annual and seasonal rainfall components as detailed in S4 Text.

Finally, we assessed if droughts and floods are becoming more frequent by counting their occurrences before and after the midpoints of the annual, wet season and dry season rainfall time series.

Multiannual persistence of severe droughts and floods.

To detect multiannual persistence of severe droughts and floods, we derived the extremal index based on the runs estimator [88] in the R-package extRemes (Eq 1 in S5 Text). An extremal index smaller than 1 indicates that droughts and floods persist over multiple years. To determine if severe droughts and floods are clustered in time without necessarily persisting over consecutive years, we additionally obtained the extremal index based on the intervals estimator [88] in the R-package extRemes (Eq 2 in S5 Text). Additionally, we tested for non-randomness in the sequence of severe drought and flood years for the annual and seasonal rainfall components using exact runs tests [89] in the R-package randtests version 1.0 [90]. Significant runs tests and values smaller than 1 for the extremal indices indicate that droughts and floods cluster in time.

Trends in the severity of droughts and floods.

To establish if droughts and floods are becoming more severe, we analysed temporal trends in extremely low (0.05 and 0.10 quantiles), high (0.90 and 0.95 quantiles) and intermediate (0.50 quantile) annual, wet season and dry season rainfall components with quantile regression (S6 Text) using the R-package quantreg version 5.26 [91]. Quantile regression is similar to ordinary linear regression in principle but can use other quantiles than the median, the only quantile handled by ordinary regression assuming normally distributed errors.

We started with an intercept-only null model and sequentially added polynomial terms of increasing order in time (yt). We retained the higher-order terms in the model if they reduced AICc (corrected Akaike Information Criterion) values compared to corresponding models with lower-order term(s) by at least two units [92]. Standard errors of the coefficients were estimated by bootstrapping using the xy-pair method and 10,000 bootstrap replications. The xy-pair method randomly resamples with replacement, pairs of the explanatory (year yt) and response (rainfall rt) variable [93].

The significance threshold was set at 0.05. All analyses (S1 File) were done in R version 3.3.1 [94] except for the UCM models (S2 File) that were implemented in the SAS UCM Procedure, SAS version 9.4 [95], SAS/ETS version 14.1 [83].

Results

Temporal trend and variation in rainfall

Temporal trend and seasonality in monthly rainfall.

The average (mean ± 1sd) total monthly rainfall was higher in the Mara (82 ± 54 mm; Fig 2A) than in Narok Town (61 ± 63 mm; Fig 2B). The level of monthly rainfall remained nearly constant in the Mara during 1965–2015 (S2A Fig) and in Narok Town during 1913–2015 (S2B Fig). There was a slight peak in the level of monthly rainfall around 1963 in Narok Town according to UCM analysis (S2B Fig).

Rainfall was markedly bimodal throughout the study period (Fig 2). There was a minor peak in December during the short rains (November-December) and a major peak in April during the long rains spanning January-June in the Mara (Fig 2A) and January-May in Narok Town (Fig 2B). The dry season spanning July-October in the Mara (Fig 2A) and June-October in Narok Town (Fig 2B) separates the long rains from the subsequent short rains. The average monthly rainfall for the dry season was much lower and more variable for Narok Town (24 ± 25 mm, CV = 104%) than for the Mara (53 ± 34 mm, CV = 64%), even though June was part of the dry season for Narok Town but not for the Mara. The average wet season rainfall was also more variable but only marginally lower for Narok Town (88 ± 69 mm, CV = 78%) than for the Mara (96 ± 56 mm, CV = 58%), despite June being part of the wet season for the Mara but not for Narok Town.

There was no significant shift in rainfall seasonality over the recording period for both the Mara and Narok Town. Rainfall seasonality was strongly deterministic (non-random) and persistent both in the Mara and in Narok Town as indicated by the highly significant seasonal components (S1 Table) that were characterized by highly insignificant disturbance variances (S2 Table). Rainfall seasonality was also remarkably stable based on cumulative averages of the total monthly rainfall for both the Mara (S3A Fig) and Narok Town in Kenya (S3B Fig).

The random sequence of standardized rainfall shows that periods with below-average dry season rainfall frequently received above-average wet season rainfall both in the Mara and in Narok Town (S4 Fig). Conversely, several periods of below-average wet season rainfall were followed by periods of above-average dry season rainfall (S4 Fig).

Temporal trend and variation in the annual and seasonal rainfall components.

There was considerable inter-annual variation in the wet season and dry season rainfall components both in the Mara (S4A Fig) and in Narok Town (S4B Fig). However, only the level components (Eq 5) of the annual (Fig 3B) and the dry season rainfall (Fig 3F) in Narok Town changed notably over time. From 1917 onwards there was a decreasing trend in the annual rainfall component that switched to an increasing trend in 1935, after which rainfall increased from a low of about 600 mm to a high of about 850 mm in 1962 (Fig 3B). Thereafter rainfall declined persistently to about 700 mm in 2015 (Fig 3B). The dry season component initially decreased from about 125 mm in 1913 to about 105 mm in 1940 and then increased to about 135 mm in 2015 (Fig 3F).

[Figure omitted. See PDF.]

Fig 3. Smoothed level component based on the structural time series analysis of the annual, wet season and dry season rainfall.

(A, C, D) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B, D, F) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B) annual, (C, D) wet season and (E, F) dry season rainfall components were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.g003

Oscillations in annual and seasonal rainfall components

The rainfall oscillations in the Mara had approximately 3-year cycles for the annual, wet season and dry season rainfall components based on the structural (UCM) time series analyses (Table 2, S5A, S5C and S5E Fig). For Narok Town, the oscillations had approximate 5.2-year cycles for the annual (S5B Fig) and wet season (S5D Fig) rainfall components and 2.3-year cycles for the dry season (S5F Fig) rainfall component (Table 2). The amplitude of the oscillations in the primary wet season cycles became wider during 1995–2015 than during the earlier period 1965–1994 in the Mara (S5C Fig) and was widest during 1950–1965 in Narok Town (S5D Fig). Structural time series analysis also identified high frequency secondary 2.2-year cycles in the annual rainfall component for the Mara (Table 2, S6A Fig) and 2.3-year cycles for Narok Town (Table 2, S6B Fig). Secondary 2.2-year cycles for the Mara (Table 2, S6C Fig) and 2.5-year cycles for Narok Town (Table 2, S6D Fig) were identified for the wet season rainfall component.

[Figure omitted. See PDF.]

Table 2. Estimated disturbance variances, damping factors and periods of the cycles in the annual, wet and dry season rainfall components.

https://doi.org/10.1371/journal.pone.0202814.t002

The stochastic cycle disturbance variances were insignificant except for the secondary 2.3-year annual rainfall cycle in Narok Town (Table 2). The stochastic disturbance variance of the primary annual cycle did not quite reach significance (Table 2). The 3-year cycles for the primary annual and wet season rainfall and the secondary 2.2-year cycles for the annual and wet season were significant for the Mara in the final state model at the end of the estimation span (Table 3). For Narok Town, the cyclic rainfall components were significant for the 2.5-year secondary wet season rainfall cycle and for the 2.3-year dry season rainfall cycle (Table 3).

[Figure omitted. See PDF.]

Table 3. Significance analysis of components (based on the final state) of annual, wet and dry season rainfall components.

https://doi.org/10.1371/journal.pone.0202814.t003

The damping factors for the rainfall cycles were approximately 1 except for the primary (0.76) and secondary (0.53) annual rainfall cycles in Narok Town (Table 2). These properties collectively imply that the annual and wet season rainfall cycles were deterministic and persistent in the Mara. For Narok Town, the primary dry season cycle and the secondary wet season cycle were also non-random and persistent. The annual rainfall cycles and the primary wet season cycle in Narok Town contained stochastic and transient elements.

Changing periodicity in monthly rainfall oscillations

Wavelet analysis established the existence of oscillations in rainfall with significantly time-varying cycle periods (periodicity). The statistically significant cycles identified for rainfall in the Mara by wavelet analysis ranged from about 0.75 to 1.5, 1.5 to 2.5 and 2.5 to 3.5 years (Fig 4A). For Narok Town, statistically significant cycles of about 0.75 to 1.5, 1 to 2.5, 2 to 3, 4 to 6 and 4.5 to 8 years were identified (Fig 4B). Seasonal rainfall oscillations with significant cycles of less than one year were distributed throughout the full span of the time series as expected (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Periodicity in monthly rainfall derived from wavelet analysis.

(A) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during the period 1965–2015. (B) Rainfall in Narok Town in Kenya was recorded during 1913–2015. The wavelet power spectrum (bias-corrected and normalized by the variance) for varying rainfall cycle periods is given for each month and year. Areas of high power are indicated in warm colours (red), whereas areas of low power are indicated in cold colours (blue). The semi-transparent area represents the cone of influence, where edge effects become important due to padding with zeros at the end of the time series. The power in this region is reduced. Significantly higher wavelet power spectra than expected under a red-noise process AR(1) (lag ‒ 1 = 0.05 for the Mara and lag ‒ 1 = 0.16 for Narok Town) are encircled by black lines.

https://doi.org/10.1371/journal.pone.0202814.g004

During 1970–1984 oscillations with a periodicity of about 0.75 to 1.5 years were apparent during three distinct episodes in the Mara: 1970–1976, 1979–1981 and 1982–1984 (Fig 3A). Significant but longer 2.5 to 3.5-year cycles characterised the period around the turn of the second millennium (1996–2002; Fig 4A). Furthermore, significant rainfall cycles of about 1.5 to 2.6 years were evident in the Mara during 2005–2007 (Fig 4A).

For Narok Town, significant 2 to 3-year rainfall cycles were apparent during the periods 1916–1919 and 1941–1946. The latter period was followed by significant oscillations of about 1.5 to 2 years during 1946–1949 (Fig 4B). Significant but longer 4 to 6-year cycles occurred during 1929–1935 and during 1997–2005 (Fig 4B). Additionally, long and significant 4.5 to 8-year cycles characterized the period 1951–1971 (Fig 4B). The episode partially overlapped with a significant, about 1 to 2.5-year cycle during 1959–1968 (Fig 4B). An oscillation with a short, approximately 0.75 to 1.5-year cycle was apparent during 2010–2012 and another during 1977–1979 (Fig 4B).

Frequency, severity and timing of droughts and floods

The temporal sequences of all severe droughts (below the 10-year return level) and floods (above the 10-year return level) of the annual and seasonal rainfall components for the Mara during 1966–2014 (dry season: 1965–2014) and for Narok Town during 1914–2015 (dry season: 1913–2015) are given in Fig 5. The return levels of extreme annual, wet season and dry season rainfall components for return periods up to 50 years for the Mara and Narok Town estimated using extreme value analysis techniques are shown in S7 Fig.

[Figure omitted. See PDF.]

Fig 5. Extreme values (circles) for the annual, wet season and dry season rainfall components.

The thin lines are the 50-year (dotted), 20-year (dashed) and 10-year (continuous) return levels. (A, C, E) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B, D, F) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B) annual, (C, D) wet season and (E, F) dry season rainfall components (thick blue lines) were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.g005

The number of severe annual (Fig 5A) and wet season (Fig 5C) droughts in the Mara rainfall time series doubled from 2 during 1965–1989 to 4 during 1990–2014. The worst annual and wet season drought did not reach the 50-year return level in 1982 (Fig 5A), followed by 2006, in which rainfall was below the 20-year return level (Fig 5C). The extreme drought of 2006 was coincident with rainfall oscillations with significant 1.5 and 2.5-year cycles (Fig 4A).

All the annual and wet season floods occurred during 1998–2012 in the Mara except for one severe annual flood in 1974 (Fig 5A and 5C). Of those, the flood in 1998 was the most extreme, exceeding the 50-year return levels (Fig 5A and 5C). This flood occurred during a time of significant rainfall periodicity of 2.5 to 3.5 years in the Mara (Fig 4A) and 4 to 6 years in Narok Town (Fig 4B).

Severe dry season droughts or floods did not increase in frequency (Fig 5E). The worst dry season drought was in 1969 with rainfall falling below the 50-year return level (Fig 5E). Another extreme drought, in which rainfall did not reach the 20-year return level, occurred during 1993 (Fig 5E). The highest dry season rainfall exceeded the 50-year return level in 2014 in the Mara (Fig 5E).

In contrast to the Mara, the frequency of severe droughts or floods did not increase in Narok Town (Fig 5B, 5D and 5F). Annual (Fig 5B), wet season (Fig 5D) and dry season (Fig 5F) droughts were more frequent up to 1953, the year with the lowest annual and wet season rainfall. The most severe drought during the dry season occurred in 1933 in Narok (Fig 5F). The droughts of 1933 and 1953 were contemporaneous with rainfall oscillations with significant, approximately 6-year cycles (Fig 4B). The majority of the severe annual (Fig 5B) and wet season (Fig 5D) floods were distributed around the most extreme flood year in 1962, in which rainfall was well above the 50-year return levels in Narok. This flood was associated with rainfall oscillations with long and significant 4.5 to 8-year cycles lasting from 1951 to 1971 (Fig 4B). The most extreme dry season rainfall exceeded the 50-year return level in 1912 and in 2011 (Fig 5F).

These results demonstrate that the severe annual and wet season droughts and floods became more frequent during 1991–2014 than during the earlier years (1965–1990) in the Mara. The severe annual and wet season floods in Narok Town were clustered around the most extreme flood year of 1962.

Multiannual persistence of severe droughts and floods

The extremal runs index established multiannual persistence in the sequence of severe droughts (below the 10-year return level) and floods (above the 10-year return level) except for the dry season droughts in the Mara, the wet season droughts in Narok Town and the annual and wet season floods in the Mara (Table 4). The extremal interval index established clustering in the sequence of severe wet season droughts in the Mara, and in the sequence of severe dry season droughts in Narok Town (Table 4). The confidence intervals for the extremal interval index for the severe dry season droughts in Narok Town did not include unity (value 1), indicating significant clustering (Table 4). All the other confidence intervals for the estimated extremal indices included unity (indicating no clustering; Table 4). The runs index and its confidence limits are unity if there is no single instance of multiannual persistence of severe drought or flood years in the annual or seasonal rainfall components. The interval index and its confidence limits are unity if the time series contains only 3 or fewer severe drought and flood years separated by 2 or less inter-exceedance times as apparent from Eq 2 in S5 Text.

[Figure omitted. See PDF.]

Table 4. Extremal runs and interval indices and significance of runs tests for clustering of severe drought and flood years.

https://doi.org/10.1371/journal.pone.0202814.t004

The p-value for the runs test was marginally significant (Table 4), suggesting that the sequence of extreme flood years in which the wet season rainfall exceeded the 20-year return level in Narok Town deviated from complete randomness. Otherwise, runs tests did not establish any other significant deviation from randomness in the occurrence of severe droughts and floods in the time series of rainfall for the Mara or Narok Town (Table 4).

Trends in the severity of droughts and floods

There was no evidence for a significant trend in the severity of droughts in the Mara from quantile regression during 1965–2014 (Table 5, Fig 6A, 6C and 6E). In Narok Town, the 0.10 quantile for the annual rainfall component increased significantly during 1914–2015 (Table 5, Fig 6B), indicating decreasing severity of droughts. The apparent decrease in the severity of the extreme (0.05 quantile) dry season droughts in the Mara (Fig 6E) and the decrease in extreme annual (Fig 6B) and wet season (Fig 6D) droughts in Narok Town were not significant (Table 5).

[Figure omitted. See PDF.]

Fig 6. Temporal trends in the severity of droughts and floods.

Blue lines are the 0.95/0.05 (dotted), 0.90/0.10 (dashed) and 0.50 (continuous) quantiles. (A, C, E) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B, D, F) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B) annual, (C, D) wet season and (E, F) dry season rainfall components (points) were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.g006

[Figure omitted. See PDF.]

Table 5. Temporal trends in rainfall extremes (0.05, 0.10, 0.50, 0.90 and 0.95 quantiles).

https://doi.org/10.1371/journal.pone.0202814.t005

Extreme wet season floods (0.95 quantile) increased significantly in the Mara during 1966–2014 (Table 5, Fig 6C). The increase in severe annual floods (0.90 quantile) was close to significance in the Mara (Table 5, Fig 6A). The increase in the extreme (0.95 quantile) annual rainfall component (Fig 6A) and increase in the 0.90 quantile of the wet season rainfall component (Fig 6C) in the Mara were not significant (Table 5).

The median (0.50 quantile) dry season rainfall component increased substantially though non-significantly from about 90 mm in 1913 to about 125 mm in 2015 in Narok Town (Table 5, Fig 6F). Otherwise, we found no evidence for significant trend over time in the median annual or seasonal rainfall components for the Mara or Narok Town (Table 5, Fig 6).

Discussion

We analysed temporal trends and variation in rainfall in the Maasai Mara ecosystem in Kenya as a background for understanding animal population and biodiversity dynamics in African terrestrial ecosystems. In contrast to IPCC’s predictions [2] and projections of the General Circulation Models [33–35], which have difficulties capturing small-scale orographic rainfall variation [96], we found only minor empirical support for increases in rainfall over two areas in Eastern Africa.

Our results show that only the dry season rainfall component for Narok Town increased during 1940–2015 based on UCM models. Likewise, the dry season rainfall increased during 1913–2001 in the Serengeti [52]. But our finding of decreasing annual rainfall during 1962–2015 in Narok Town reinforces results of earlier studies [36,44–46] that attributed decreasing East African rainfall to increasing temperatures of the northern hemisphere [36] and of the central Indian or west Pacific Oceans [44–46]. Similarly, Ritchie et al. [52] documented evidence of decreasing annual and wet season rainfall during 1960–2001 in the Serengeti. Rainfall also declined during 1960–2014 in 14 Kenyan counties but in another 6 Kenyan counties it declined initially and then switched to an upward trend according to Ogutu et al. [30], demonstrating spatial distinctions in rainfall trends. There was, however, no systematic change in the annual and seasonal rainfall components during 1965–2014 in the Mara according to our UCM analysis. The apparent discrepancies in the contrasting findings may partly reflect the strong spatial variation evident in East-African rainfall [55,97] and calls for considerable caution in applying results of climate modelling at large spatial scales to particular localities.

Our result demonstrating virtually constant rainfall seasonality in the Mara or Narok Town suggests that regular movements of the Inter-Tropical Convergence Zone provide an extremely stable modulation of the rainfall seasons. The higher and less variable rainfall in the Mara than in Narok Town, in particular the higher and less variable dry season rainfall, is likely due to the closer proximity of the Mara to the high-precipitation areas near the eastern shores of Lake Victoria [55] and the influence of the Lake Victoria System. The modulation of Mara’s rainfall by this system and the resulting higher dry season rainfall likely sustains the high abundance and diversity of wildlife and livestock there. The seasonal variation in rainfall both in the Mara and Narok Town agree well with the monthly discharge patterns for the Mara River, with peaks in December and April-May [98], implying that rainfall governs much of the Mara River discharge.

The compensatory pattern we detected for the amounts of the wet and dry season rainfall, where periods with below-average dry season rainfall frequently received above-average wet season rainfall and vice versa, may possibly emerge from the interactions between the local Lake Victoria circulation and hemispheric level climate drivers. This compensatory dampening of rainfall seasonality has likely generated the stability and resilience of the Mara and allowed the high abundance and diversity of wildlife it supports. The positive Indian Ocean Dipole Mode and the El Niño–Southern Oscillation can increase East African rainfall primarily during the wet season by weakening westerly winds and enabling moisture transport from the Indian Ocean towards the land [42,43,99]. But weakened westerly winds may, in turn, reduce the dry season rainfall that typically originates from Lake Victoria and precipitates east of the lake in the Mara [55]. Correspondingly, La Niña-like conditions and negative Indian Ocean Dipole phases may enhance westerly winds from Lake Victoria that lead to above-average dry season rainfall but concurrently block moisture export from the Indian Ocean towards East Africa resulting in less wet season rainfall.

The predominantly persistent rainfall cycle periods we estimated of between 2.1 and 3 years for the Mara and 2.3 to 5.2 years for Narok Town agree with the dominant cycle periods reported for East Africa ranging from about 2 to 11 years [56–58,76]. They are also within the range of the periods characteristic of the Indian Ocean Dipole (about 2 years), the El Niño-Southern Oscillation (3–6 years; [100]) and the quasi-biennial oscillations in the lower equatorial stratospheric zonal winds [101].

By means of wavelet analysis we were able to establish that periodicity in the monthly rainfall oscillations is not constant but varies over time and space in the Mara-Serengeti ecosystem. The cycle periods varied from just 0.75 to 1.5 years for several distinct episodes during 1970–1984 for the Mara to periods ranging from 4.5 to 8 years during 1951–1971 for Narok Town. Likewise, the amplitude and frequency of the Indian Ocean Dipole and the El Niño-Southern Oscillation can be highly variable [100,102].

In general, the wavelet analysis showed that extreme droughts and floods coincided with the times when the regional rainfall oscillations had statistically significant cycle periods. These results strongly suggest a possible association between oscillations in hemispheric atmospheric and oceanic circulations and regional rainfall oscillations in the Mara-Serengeti ecosystem that are similarly evident in East African rainfall [103–105]. For example, the first half of the wet season of 2006 was associated with both a strong negative Indian Ocean Dipole and La Niña-like conditions [50,51] with consequent severe drying in East Africa [48,75,82]. The flood of 1962 coincided with strong easterly winds [106] and an extreme dipole reversal of the Indian Ocean sea surface temperatures [42]. But there was no simultaneous El Niño effect during that time [42]. Similarly, the flood of 1998 was associated with the second strongest El Niño event on instrumental record during 1950–2016 [50], and coincided with a strong and positive Indian Ocean Dipole [42,50,51]. These circulations may also explain the simultaneous occurrence of overlapping rainfall cycles that were apparent both in UCM models and wavelet analysis.

But the extreme dry season droughts of 1969 and 1993, and the extreme annual and wet season drought of 1982 in the Mara were not associated with any significant periodicity in rainfall oscillation. There was a moderate El Niño event during 1968–1969 and 1993, and the longest recorded El Niño episode during 1990–1995 [50]. The interactions between the El Niño events and Lake Victoria circulation patterns that probably underlie the compensatory pattern of the wet and dry season rainfall could also have contributed to the dry season droughts of 1969 and 1993.

Our results indicating increasing frequency and severity of floods in the Mara during 1965–2014 and reduction in drought frequency and severity in Narok Town during 1913–2015 agree with the projections of General Circulation Models for most of East Africa of more intense wet seasons [107] and less severe droughts [74]. Precipitation extremes are expected as temperatures rise [107]. Such patterns may emerge from the projected strengthening of the El Niño Southern Oscillation [108,109] and more frequent occurrences of the positive phases of the Indian Ocean Dipole [38,39].

In contrast to Narok, severe droughts are apparently becoming more frequent in the Mara. Although some of the inconsistencies in our results might have arisen from using only one gauge in Narok Town versus 15 in the Mara, they show that the trends in extreme rainfall events can vary both locally and over time. Simulation models indicate that, although droughts may increase with increasing temperatures on a global scale, the changes expected for East Africa appear uncertain and comparatively small [47]. One possible cause of the recent severe droughts in the Mara could be due to the increasing frequency of extreme La Niña events that often follow extreme El Niño events [47]. These hemispheric circulation patterns may also underlie the multiannual persistence of severe droughts and floods that we detected in the annual and seasonal rainfall components in the Mara and in Narok Town.

Implications for animal population and biodiversity dynamics

The decreasing annual rainfall in Narok Town may have negative long-term effects on the reproductive performance of herbivores, including livestock, in that area through limitation of food and surface water availability [5,6,10,110]. Wildlife and livestock populations are also likely to suffer more from water and food scarcity in the Mara despite the stable rainfall levels during 1965–2014 due to the rising temperatures and increasing human impacts [30]. In addition, the combination of decreasing rainfall, rising temperatures and hence increased rates of evapo-transpiration can make wildfires more destructive and severe [24]. Such conditions can have potentially adverse effects on animal populations and biodiversity [30]. By contrast, the increase in dry season rainfall in Narok Town during 1940–2015 can improve the survival prospects of ungulates when resources are most limiting [10,16].

The predominantly deterministic and persistent primary rainfall cycles in the Mara likely have important implications for the dynamics and management of animal populations and their vegetation resources [5,14]. The changing periodicity in the rainfall oscillations and the increasing amplitude in the wet season rainfall oscillations in the Mara may compound impacts of the recurrent seasonal and cyclic variations in water and forage availability that animals already have to cope with by adaptively adjusting their reproduction, foraging and migration patterns.

The increase in excessive wet season rainfall in the Mara can displace wildlife due to flooding [29,111], reduce recruitment [112] and forage quality due to excessive plant growth and dilution of plant nutrients and cause population declines [28,74]. The increasing severity of floods creates favourable conditions for the transmission of several diseases, including anthrax [113], Rift Valley Fever [114] and African horse sickness [115] and promotes infestation with parasites [20]. In the Serengeti, abnormally high wet season rainfall may enable wildebeest to use the nutrient-rich southern short grass plains for extended periods towards the dry season [25]. As a consequence, the excessive rainfall may reduce wildebeest’s occupancy of their dry season range in the Mara in Kenya [25]. Rainfall-mediated migrations and local concentrations of large herbivores can have cascading and difficult-to-predict effects on the ecosystem by modifying nutrient cycling through grazing, urine and dung deposition [116] and through prey availability for predators and scavengers [117,118].

If the pattern of increasing drought frequency continues in the Mara it will likely have adverse and immediate impacts on both wildlife and livestock populations in the ecosystem [16,29,75,80,82,119], although some uncertainty remains owing to the effect of the opposing trend in drought frequency in Narok Town. Droughts also tend to be associated with outbreaks of infectious diseases among large herbivores [120,121], such as anthrax that spreads when herbivores graze short grass close to the ground. In particular, multi-year droughts and floods can have much stronger adverse impacts on herbivores than single-year droughts [15,75]. The changing drought and flood frequencies can also affect biodiversity in the region through their impact on the Mara River flow levels.

Conclusions

For important and unique ecosystems such as the Maasai Mara in East Africa it is crucial to keep reliable historical and comprehensive contemporary climate records to provide a sound scientific basis for informing decisions on the nature and consequences of climate change to large herbivores and their predators, biodiversity and human livelihoods. Our analysis of verified station data for the Mara region and Narok Town in Kenya indicates that local trends and variation in rainfall can differ substantially from the patterns predicted for regional or continental scales. In particular, we found only minor support for IPCC’s large-scale predictions of very likely increases in East African rainfall. Although the rainfall cycle periods in our study area correspond to the dominant cycle periods evident in East African rainfall, we found strong spatial and temporal variation in rainfall periodicity in the Maasai Mara ecosystem. While droughts are apparently becoming more frequent in the Mara, we detected the opposite pattern for Narok Town located only 75 km away. Similarly, wet season floods became more severe in the Mara but not in Narok Town. Significant changes in rainfall amounts and periodicity can profoundly affect animal population and biodiversity dynamics by lowering the availability and increasing the variability of food and water resources, increasing the risk of outbreaks of infectious diseases and altering ungulate migration and dispersal patterns.

Supporting information

[Figure omitted. See PDF.]

S1 Data. Mara monthly rainfall 1965–2015.

Rainfall data comprising recordings from 15 stations in the Mara region of Kenya. KMD, Kenya Meteorological Department; MMEMP, Maasai Mara Ecological Monitoring Programme; NA, not available; Rainfall (mm), records verified for analysis; Source of verification; Rain_imp, rainfall records in mm including imputed values used in analyses.

https://doi.org/10.1371/journal.pone.0202814.s001

(XLSX)

S2 Data. Narok monthly rainfall 1913–2015.

Rainfall data comprising recordings from the gauge in Narok Town operated by the Kenya Meteorological Department.

https://doi.org/10.1371/journal.pone.0202814.s002

(XLSX)

S1 Text. Rainfall standardisation methods.

https://doi.org/10.1371/journal.pone.0202814.s003

(DOCX)

S2 Text. Modelling seasonal oscillations in rainfall.

https://doi.org/10.1371/journal.pone.0202814.s004

(DOCX)

S3 Text. Using wavelet analysis to detect changing periodicity in monthly rainfall oscillations.

https://doi.org/10.1371/journal.pone.0202814.s005

(DOCX)

S4 Text. Threshold selection for estimating return levels of droughts and floods.

https://doi.org/10.1371/journal.pone.0202814.s006

(DOCX)

S5 Text. Estimation of the runs and intervals estimator to detect multiannual persistence of droughts and floods.

https://doi.org/10.1371/journal.pone.0202814.s007

(DOCX)

S6 Text. Quantile regression to analyse trends in the severity of droughts and floods.

https://doi.org/10.1371/journal.pone.0202814.s008

(DOCX)

S1 File. R code used to analyse the time series of rainfall derived from 15 gauges in the Mara (1966–2014) and a gauge in Narok Town (1913–2015).

The code is provided for wavelet analysis for the total monthly rainfall, and extreme value analysis, extremal indices, runs tests and quantile regression for the annual and seasonal rainfall components. The standardisation methods for the Mara data are included.

https://doi.org/10.1371/journal.pone.0202814.s009

(TXT)

S2 File. SAS code used to fit UCM models for the time series of rainfall derived from 15 gauges in the Mara (1966–2014) and a gauge in Narok Town (1913–2015).

The analysis pertains to the time series of the total monthly rainfall, and to the annual and seasonal rainfall components.

https://doi.org/10.1371/journal.pone.0202814.s010

(TXT)

S1 Fig. Standardisation methods for the time series of annual rainfall derived from 15 gauges in the Mara (1966–2014).

Four different standardisation methods were applied: arithmetic mean (black line), mean-weighted average (orange line), gauge-adjusted mean (blue line) and Best Linear Unbiased Predictions (BLUPS) from a generalized linear mixed model with Tweedie error distribution and a log link function (red line).

https://doi.org/10.1371/journal.pone.0202814.s011

(TIFF)

S2 Fig. Smoothed level component based on the structural time series analysis of monthly rainfall.

(A) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) recording in the Mara during 1965–2014 and in (B) Narok Town in Kenya during 1913–2015. Significant components are marked in bold-faced font.

https://doi.org/10.1371/journal.pone.0202814.s012

(TIF)

S3 Fig. Cumulative plot of decadal averages of the total monthly rainfall.

(A) Rainfall in the Mara during 1965–2014. (B) Rainfall in Narok Town in Kenya during 1913–2015.

https://doi.org/10.1371/journal.pone.0202814.s013

(TIF)

S4 Fig. The standardized (divided by the mean) moving averages of the wet season (red lines) and dry season (blue lines) rainfall components.

The vertical needles are the standardized deviates and the solid curves are the 3-year (Mara), 5-year (wet season of Narok Town) and 2-year (dry season of Narok Town) moving averages. (A) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The wet season and dry season rainfall components were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.s014

(TIF)

S5 Fig. Smoothed primary cycles for standardized rainfall based on the structural time series analysis.

(A, C, D) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B, D, F) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B) annual, (C, D) wet season and (E, F) dry season rainfall components were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.s015

(TIF)

S6 Fig. Smoothed secondary cycles for standardized rainfall based on the structural time series analysis.

(A, C) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (B, D) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B) annual and (C, D) wet season rainfall components were summed from the monthly rainfall records.

https://doi.org/10.1371/journal.pone.0202814.s016

(TIF)

S7 Fig. Return periods and the corresponding return levels (blue lines) for extreme rainfall.

(A, B, E, F, I, J) Rainfall recordings in the Mara were derived from 15 gauges (Eq 3 in S1 Text) during 1966–2014 (dry season: 1965–2014). (C, D, G, H, K, L) Rainfall in Narok Town in Kenya was recorded during 1914–2015 (dry season: 1913–2015). The (A, B, C, D) annual, (E, F, G, H) wet season and (I, J, K, L) dry season rainfall components were summed from the monthly rainfall records. Blue polygons are the 95% normal approximate confidence bands.

https://doi.org/10.1371/journal.pone.0202814.s017

(TIFF)

S1 Table. Significance analysis of components (based on the final state) of monthly rainfall.

https://doi.org/10.1371/journal.pone.0202814.s018

(XLSX)

S2 Table. The estimated variances of the disturbance terms and the damping factor of the autoregressive component for the monthly rainfall series.

https://doi.org/10.1371/journal.pone.0202814.s019

(XLSX)

Acknowledgments

We thank Professor Kay E. Holekamp (Michigan State University) and the Kenya Meteorological Department (KMD) for providing the rainfall data as well as helpful advice. In particular, we are greatly indebted to Ms Christine Mahonga for her patience and excellent support in making the data available to us and coordinating the verification of the rainfall data at KMD. The Maasai Mara Ecological Monitoring Programme (MMEMP) was designed and supervised by Dr Holly T. Dublin, supported by Professor A.R.E. Sinclair (University of British Columbia), and executed by Messrs. Paul Chara (July 1989–1992), John Naiyoma (1989–1993), Alex Obara (1995–1997) and Charles Matankory (1991–2003). It was funded by the Worldwide Fund for Nature-East Africa Programme (WWF-EARPO) and Friends of Conservation (FOC). The programme also received financial, material or logistical support from WWF-Sweden, the Darwin Initiative through DICE, Cottar’s Camp, Kichwa Tembo and Ker and Downey Safaris. Many other individuals and organizations, too many to list individually, also helped in various ways over the last 40 years. We are grateful to Mr. Gordon Ojwang (Directorate of Resource Surveys and Remote Sensing of Kenya) for his assistance in acquiring the rainfall data for Narok for 1970–2014 from KMD. We are also grateful to Professor Andrew Dobson (Princeton University) for suggesting the use of wavelet analysis and to Professor Norman Owen-Smith (University of the Witwatersrand) for providing valuable insights. Three anonymous reviewers and Dr. Shimelis Dessu provided detailed constructive criticisms and suggestions that helped considerably improve an earlier draft of this paper. JOO was supported by a grant from the German Research Foundation (Grant No. OG 83/1-1). This project has received funding from the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No. 641918.

Citation: Bartzke GS, Ogutu JO, Mukhopadhyay S, Mtui D, Dublin HT, Piepho H-P (2018) Rainfall trends and variation in the Maasai Mara ecosystem and their implications for animal population and biodiversity dynamics. PLoS ONE 13(9): e0202814. https://doi.org/10.1371/journal.pone.0202814

References

1. Maclean IMD, Wilson RJ. Recent ecological responses to climate change support predictions of high extinction risk. Proc Natl Acad Sci. 2011;108: 12337–12342. pmid:21746924

2. Niang I, Ruppel OC, Abdrabo MA, Essel A, Lennard C, Padgham J, et al. 2014: Africa. In: Barros VR, Field CB, Dokken DJ, Mastrandrea MD, Mach KJ, Bilir TE, et al., editors. Climate change 2014: impacts, adaptation, and vulnerability. Part B: regional aspects. Contribution of working Group II to the fifth assessment report of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, NY, USA: Cambridge University Press 2014. pp. 1199–1265.

3. Sylla MB, Giorgi F, Coppola E, Mariotti L. Uncertainties in daily rainfall over Africa: assessment of gridded observation products and evaluation of a regional climate model simulation. Int J Climatol. 2013;33: 1805–1817.

4. Bartzke GS, Ogutu JO, Mtui D, Mukhopadhyay S, Piepho H-P. Modelled trends in climate and vegetation. 2016. Report No.: 3.1 Available from: https://www.researchgate.net/profile/Joseph_Ogutu2/publication/308802687_Modelling_trends_and_variation_in_climate_of_the_Serengeti-Mara_Ecosystem/links/57f3740508ae8da3ce51b363.pdf.

5. Ogutu JO, Owen-Smith N. Oscillations in large mammal populations: are they related to predation or rainfall? Afr J Ecol. 2005;43: 332–339.

6. Ogutu JO, Piepho H-P, Dublin HT, Bhola N, Reid RS. Rainfall influences on ungulate population abundance in the Mara-Serengeti ecosystem. J Anim Ecol. 2008;77: 814–829. pmid:18422558

7. Boutton TW, Tieszen LL, Imamaba SK. Biomass dynamics of grassland vegetation in Kenya. Afr J Ecol. 1988;26: 89–101.

8. Sankaran M, Hanan NP, Scholes RJ, Ratnam J, Augustine DJ, Cade BS, et al. Determinants of woody cover in African savannas. Nature. 2005;438: 846–849. pmid:16341012

9. Boutton TW, Tieszen LL, Imbamba SK. Seasonal changes in the nutrient content of East African grassland vegetation. Afr J Ecol. 1988;26: 103–115.

10. Owen-Smith N, Mason DR, Ogutu JO. Correlates of survival rates for 10 African ungulate populations: density, rainfall and predation. J Anim Ecol. 2005;74: 774–788.

11. Mills MGL, Biggs HC, Whyte IJ. The relationship between rainfall, lion predation and population trends in African herbivores. Wildl Res. 1995;22: 75–87.

12. Coe MJ, Cumming DH, Phillipson J. Biomass and production of large African herbivores in relation to rainfall and primary production. Oecologia. 1976;22: 341–354. pmid:28308896

13. Fritz H, Duncan P. On the carrying capacity for large ungulates of African savanna ecosystems. Proc Biol Sci. 1994;256: 77–82. pmid:8008761

14. Ogutu JO, Piepho H-P, Dublin HT. Reproductive seasonality in African ungulates in relation to rainfall. Wildl Res. 2014;41: 323–342.

15. Hillmann JC, Hillman AKK. Mortality of wildlife in Nairobi National Park, during the drought of 1973–1974. Afr J Ecol. 1977;15: 1–18.

16. Mduma SAR, Sinclair ARE, Hilborn R. Food regulates the Serengeti wildebeest: a 40-year record. J Anim Ecol. 1999;68: 1101–1122.

17. Chamaillé-Jammes S, Fritz H, Valeix M, Murindagomo F, Clobert J. Resource variability, aggregation and direct density dependence in an open context: the local regulation of an African elephant population. J Anim Ecol. 2008;77: 135–144. pmid:17986249

18. Landman M, Schoeman DS, Hall-Martin AJ, Kerley GIH. Understanding long-term variations in an elephant piosphere effect to manage impacts. PLoS One. 2012;7: e45334. pmid:23028942

19. Smuts GL. Interrelations between predators, prey, and their environment. Bioscience. 1978;28: 316–320.

20. Fosbrooke H. The stomoxys plague in Ngorongoro, 1962. Afr J Ecol. 1963;1: 124–126.

21. Roelke-Parker ME, Munson L, Packer C, Kock R, Cleaveland S, Carpenter M, et al. A canine distemper virus epidemic in Serengeti lions (Panthera leo). Nature. 1996;379: 441–445. pmid:8559247

22. Olff H, Ritchie ME, Prins HHT. Global environmental controls of diversity in large herbivores. Nature. 2002;415: 901–904. pmid:11859367

23. Van Wilgen BW, Govender N, Biggs HC, Ntsala D, Funda XN. Response of savanna fire regimes to changing fire-management policies in a large African national park. Conserv Biol. 2004;18: 1533–1540.

24. Hemp A. Climate change-driven forest fires marginalize the impact of ice cap wasting on Kilimanjaro. Global Change Biol. 2005;11: 1013–1023.

25. Pennycuick L. Movements of the migratory wildebeest population in the Serengeti area between 1960 and 1973. Afr J Ecol. 1975;13: 65–87.

26. Young K, Van Aarde R. Density as an explanatory variable of movements and calf survival in savanna elephants across southern Africa. J Anim Ecol. 2010;79: 662–673. pmid:20180876

27. Loveridge AJ, Hunt JE, Murindagomo F, Macdonald DW. Influence of drought on predation of elephant (Loxodonta africana) calves by lions (Panthera leo) in an African wooded savannah. J Zool. 2006;270: 523–530.

28. Ogutu JO, Owen-Smith N, Piepho H-P, Kuloba B, Edebe J. Dynamics of ungulates in relation to climatic and land use changes in an insularized African savanna ecosystem. Biodivers Conserv. 2012;21: 1033–1053.

29. Ogutu JO, Owen-Smith N, Piepho H-P, Said MY, Kifugo S, Reid RS, et al. Changing wildlife populations in Nairobi National Park and adjoining Athi-Kaputiei Plains: collapse of the migratory wildebeest. Open Conserv Biol J. 2013;7: 11–26.

30. Ogutu JO, Piepho H-P, Said MY, Ojwang GO, Njino LW, Kifugo SC, et al. Extreme wildlife declines and concurrent increase in livestock numbers in Kenya: What are the causes? PLoS ONE. 2016;11: e0163249. pmid:27676077

31. Schreck CJ, Semazzi FHM. Variability of the recent climate of eastern Africa. Int J Climatol. 2004;24: 681–701.

32. Kizza M, Rodhe A, Xu C-Y, Ntale HK, Halldin S. Temporal rainfall variability in the Lake Victoria basin in East Africa during the twentieth century. Theor Appl Climatol. 2009;98: 119–135.

33. Giannini A, Biasutti M, Held IM, Sobel AH. A global perspective on African climate. Clim Change. 2008;90: 359–383.

34. McHugh MJ. Multi-model trends in East African rainfall associated with increased CO2. Geophys Res Lett. 2005;32: L01707.

35. McSweeney CF, Jones RG. No consensus on consensus: the challenge of finding a universal approach to measuring and mapping ensemble consistency in GCM projections. Clim Change. 2013;119: 617–629.

36. Tierney JE, Ummenhofer CC, deMoneocal PB. Past and future rainfall in the Horn of Africa. Sci Adv. 2015;1: e1500682. pmid:26601306

37. Shongwe ME, van Oldenborgh GJ, van den Hurk B, van Aalst M. Projected changes in mean and extreme precipitation in Africa under plobal warming. Part II: East Africa. J Climate. 2011;24: 3718–3733.

38. Cai W, Cowan T, Sullivan A. Recent unprecedented skewness towards positive Indian Ocean Dipole occurrences and its impact on Australian rainfall. Geophys Res Lett. 2009;36: 1–4.

39. Dore MHI. Climate change and changes in global precipitation patterns: What do we know? Environ Int. 2005;31: 1167–1181. pmid:15922449

40. Collins M, An S-I, Cai W, Ganachaud A, Guilyardi E, Jin F-F, et al. The impact of global warming on the tropical Pacific Ocean and El Niño. Nat Geosci. 2010;3: 391–397.

41. Fedorov AV, Philander SG. Is El Niño changing? Science. 2000;288: 1997–2002. pmid:10856205

42. Saji NH, Goswami BN, Vinayachandran PN, Yamagata T. A dipole mode in the tropical Indian Ocean. Nature. 1999;401: 360–363. pmid:16862108

43. Webster PJ, Moore AM, Loschnigg JP, Leben RR. Coupled ocean–atmosphere dynamics in the Indian Ocean during 1997–98. Nature. 1999;401: 356–360. pmid:16862107

44. Funk C, Dettinger MD, Michaelsen JC, Verdin JP, Brown ME, Barlow M, et al. Warming of the Indian Ocean threatens eastern and southern African food security but could be mitigated by agricultural development. Proc Natl Acad Sci. 2008;105: 11081–11086. pmid:18685101

45. Lyon B, DeWitt DG. A recent and abrupt decline in the East African long rains. Geophys Res Lett. 2012;39: L02702.

46. Williams AP, Funk C. A westward extension of the warm pool leads to a westward extension of the Walker circulation, drying eastern Africa. Climate Dyn. 2011;37: 2417–2435.

47. Prudhomme C, Giuntoli I, Robinson EL, Clark DB, Arnell NW, Dankers R, et al. Hydrological droughts in the 21st century, hotspots and uncertainties from a global multimodel ensemble experiment. Proc Natl Acad Sci. 2014;111: 3262–3267. pmid:24344266

48. Hastenrath S, Polzin D, Mutai C. Diagnosing the 2005 drought in equatorial East Africa. J Climate. 2007;20: 4628–4637.

49. Owiti Z, Ogallo LA, Mutemi J. Linkages between the Indian Ocean Dipole and East African seasonal rainfall anomalies. J Kenya Meteorol Soc. 2008;2: 3–17.

50. Barnard PL, Hoover D, Hubbard DM, Snyder A, Ludka BC, Allan J, et al. Extreme oceanographic forcing and coastal response due to the 2015–2016 El Niño. Nat Commun. 2017;8: 1–8.

51. Ham Y-G, Choi J-Y, Kug J-S. The weakening of the ENSO–Indian Ocean Dipole (IOD) coupling strength in recent decades. Climate Dyn. 2016; 1–13.

52. Ritchie ME. Global environmental changes and their impact on the Serengeti. In: Sinclair ARE, Packer C, Mduma SAR, Fryxell JM, editors. Serengeti III: human impacts on ecosystem dynamics. Chicago, Illinois, USA and London, UK: University of Chicago Press 2008. pp. 183–208.

53. Ogutu JO, Piepho H-P, Dublin HT, Bhola N, Reid RS. El Niño-Southern Oscillation, rainfall, temperature and Normalized Difference Vegetation Index fluctuations in the Mara-Serengeti ecosystem. Afr J Ecol. 2008;46: 132–143.

54. Green DS, Roloff GJ, Heath BR, Holekamp KE. Temporal dynamics of the reponses by African mammals to prescribed fire. J Wildl Manage. 2015;79: 235–242.

55. Norton-Griffiths M, Herlocker D, Pennycuick L. The patterns of rainfall in the Serengeti ecosystem, Tanzania. Afr J Ecol. 1975;13: 347–374.

56. Ogallo L. Temporal fluctuations of seasonal rainfall patterns in East Africa. Mausam. 1984;35: 175–180.

57. Ogutu JO, Kuloba B, Piepho H-P, Kanga E. Wildlife population dynamics in human-dominated landscapes under community-based conservation: the example of Nakuru Wildlife Conservancy, Kenya. PLoS One. 2017;12: e0169730. pmid:28103269

58. Rodhe H, Virji H. Trends and periodicities in East African rainfall data. Mon Weather Rev. 1976;104: 307–315.

59. Mutai CC, Ward MN. East African rainfall and the tropical circulation/convection on intraseasonal to interannual timescales. J Climate. 2000;13: 3915–3939.

60. Mwale D, Gan TY. Wavelet Analysis of Variability, Teleconnectivity, and Predictability of the September-November East African Rainfall. J Appl Meteorol. 2005;44: 256–269.

61. Rodríguez-Fonseca B, Janicot S, Mohino E, Losada T, Bader J, Caminade C, et al. Interannual and decadal SST-forced responses of the West African monsoon. Atmos Sci Lett. 2011;12: 67–74.

62. Nicholson SE. A review of climate dynamics and climate variability in eastern Africa. In: Johnson TC, Odada EO, editors. The limnology, climatology and paleoclimatology of the East African lakes. Amsterdam, Netherlands: Gordon and Breach 1996. pp. 25–56.

63. Sinclair ARE, Hopcraft JGC, Olff H, Mduma SAR, Galvin KA, Sharam GJ. Historical and future changes to the Serengeti ecosystem. In: Sinclair ARE, Packer C, Mduma SAR, Fryxell JM, editors. Serengeti III: human impacts on ecosystem dynamics. Chicago, Illinois, USA and London, UK: University of Chicago Press 2008. pp. 7–46.

64. Serneels S, Said MY, Lambin EF. Land cover changes around a major East African wildlife reserve: the Mara ecosystem (Kenya). Int J Remote Sens. 2001;22: 3397–3420.

65. Mduma SAR, Hopcraft JGC. The main herbivorous mammals and crocodiles in the greater Serengeti ecosystem. In: Sinclair ARE, Packer C, Mduma SAR, Fryxell JM, editors. Serengeti III: human impacts on ecosystem dynamics. Chicago, Illinois, USA and London, UK: University of Chicago Press 2008. pp. 497–505.

66. Sinclair ARE, Mduma S, Brashares JS. Patterns of predation in a diverse predator–prey system. Nature. 2003;425: 288–290. pmid:13679915

67. Ogutu JO, Piepho H-P, Dublin HT, Bhola N, Reid RS. Dynamics of Mara–Serengeti ungulates in relation to land use changes. J Zool. 2009;278: 1–14.

68. Ogutu J, Dublin H. Spatial dynamics of lions and their prey along an environmental gradient. Afr J Ecol. 2004;42: 8–22.

69. Fraedrich K. A simple climatological model of the dynamics and energetics of the nocturnal circulation at Lake Victoria. Quart J Roy Meteor Soc. 1972;98: 322–335.

70. Anyah RO, Semazzi FHM. Variability of East African rainfall based on multiyear RegCM3 simulations. Int J Climatol. 2007;27: 357–371.

71. Wolanski E, Gereta E. Water quantity and quality as the factors driving the Serengeti ecosystem, Tanzania. Hydrobiologia. 2001;458: 169–180.

72. Murray MG. Specific Nutrient Requirements and Migration of Wildebeest. In: Sinclair ARE, editor. Serengeti II: dynamics, management and conservation of an ecosystem. Chicago, Illinois, USA and London, UK: The University of Chicago Press 1995. pp. 231–256.

73. Wolanski E, Gereta E, Borner M, Mduma S. Water, migration and the Serengeti ecosystem. Amer Sci. 1999;87: 526.

74. Ogutu JO, Piepho H-P, Kanga E. Dynamics of an insularized and compressed impala population: rainfall, temperature and density influences. Open Ecol J. 2012;5: 1–17.

75. Ogutu JO, Piepho H-P, Said MY, Kifugo SC. Herbivore dynamics and range contraction in Kajiado county Kenya: climate and land use changes, population pressures, governance, policy and human-wildlife conflicts. Open Ecol J. 2014;7: 9–31.

76. Pennycuick L, Norton-Griffiths M. Fluctuations in the rainfall of the Serengeti Ecosystem, Tanzania. J Biogeogr. 1976;3: 125–140.

77. Gereta E, Mwangomo E, Wolanski E. Ecohydrology as a tool for the survival of the threatened Serengeti ecosystem. Ecohydrol Hydrobiol. 2009;9: 115–124.

78. Piepho H-P, Ogutu JO. Simple state-space models in a mixed model framework. Am Stat. 2007;61: 224–232.

79. Angassa A, Oba G. Relating long-term rainfall variability to cattle population dynamics in communal rangelands and a government ranch in southern Ethiopia. Agric Syst. 2007;94: 715–725.

80. Dublin HT, Ogutu JO. Population regulation of African buffalo in the Mara–Serengeti ecosystem. Wildl Res. 2015;42: 382–393.

81. Kimwele CN, Graves JA. A molecular genetic analysis of the communal nesting of the ostrich (Struthio camelus). Mol Ecol. 2003;12: 229–236. pmid:12492891

82. Nkedianye D, de Leeuw J, Ogutu JO, Said MY, Saidimu TL, Kifugo SC, et al. Mobility and livestock mortality in communally used pastoral areas: the impact of the 2005–2006 drought on livestock mortality in Maasailand. Pastoralism. 2011;1: 1–17.

83. SAS Institute Inc. SAS/ETS® 14.1 user’s guide: the UCM procedure [software]. Cary, North Carolina, USA: SAS Institute Inc. 2015. Available from: https://support.sas.com/documentation/onlinedoc/ets/141/ucm.pdf

84. Durbin J, Koopman SJ. Time series analysis by state space methods. Oxford, UK and New York, USA: Oxford University Press 2012.

85. Gouhier TC, Grinsted A, Simko V. Biwavelet: conduct univariate and bivariate wavelet analyses [software]. 2016. Available from: http://github.com/tgouhier/biwavelet

86. Torrence C, Compo GP. A practical guide to wavelet analysis. Bull Amer Meteor Soc. 1998;79: 61–78.

87. Gilleland E, Katz RW. New software to analyze how extremes change over time. Eos. 2011;92: 13–14.

88. Ferro CAT, Segers J. Inference for clusters of extreme values. J R Stat Soc Series B Stat Methodol. 2003;65: 545–556.

89. Gibbons JD, Chakraborti S. Tests of randomness. In: Balakrishnan N, Schucany WR, editors. Nonparametric statistical inference. Boca Raton, Florida, USA: Chapman & Hall/CRC 2011. pp. 75–99.

90. Caeiro F, Mateus A. Randtests: testing randomness in R [software]. 2014. Available from: https://CRAN.R-project.org/package=randtests

91. Koenker R. Quantreg: quantile Regression [software]. 2016. Available from: https://cran.r-project.org/src/contrib/Archive/quantreg/

92. Burnham KP, Anderson D. Model selection and multi-model inference: a practical information-theoretic approach. New York, New York, USA: Springer, Inc. 2002.

93. Koenker R. Confidence intervals for regression quantiles. In: Mandl P, Hušková M, editors. Asymptotic statistics: proceedings of the fifth prague symposium, held from September 4–9, 1993. Heidelberg, Baden-Württemberg, Germany: Physica-Verlag HD 1994. pp. 349–359.

94. R Core Team. R: a language and environment for statistical computing [software]. Vienna, Austria: R Foundation for Statistical Computing 2016. Available from: https://www.R-project.org/

95. SAS Institute Inc. SAS system for windows [software]. Cary, North Carolina, USA: 2016.

96. Dessu SB, Melesse AM. Evaluation and Comparison of Satellite and GCM Rainfall Estimates for the Mara River Basin, Kenya/Tanzania. In: Younos T, Grady CA, editors. Climate Change and Water Resources. Berlin, Heidelberg: Springer Berlin Heidelberg 2013. pp. 29–45.

97. Omondi P, Awange JL, Ogallo LA, Okoola RA, Forootan E. Decadal rainfall variability modes in observed rainfall records over East Africa and their relations to historical sea surface temperature changes. J Hydrol. 2012;464–465: 140–156.

98. McClain ME, Subalusky AL, Anderson EP, Dessu SB, Melesse AM, Ndomba PM, et al. Comparing flow regime, channel hydraulics, and biological communities to infer flow–ecology relationships in the Mara River of Kenya and Tanzania. Hydrol Sci J. 2014;59: 801–819.

99. Tierney JE, Smerdon JE, Anchukaitis KJ, Seager R. Multidecadal variability in East African hydroclimate controlled by the Indian Ocean. Nature. 2013;493: 389–392. pmid:23325220

100. Behera SK, Yamagata T. Influence of the Indian Ocean dipole on the Southern Oscillation. J Meteor Soc Japan Ser II. 2003;81: 169–177.

101. Indeje M, Semazzi FHM. Relationships between QBO in the lower equatorial stratospheric zonal winds and East African seasonal rainfall. Meteorol Atmos Phys. 2000;73: 227–244.

102. Cobb KM, Charles CD, Cheng H, Edwards RL. El Niño/Southern Oscillation and tropical Pacific climate during the last millennium. Nature. 2003;424: 271–276. pmid:12867972

103. Behera SK, Luo J-J, Masson S, Delecluse P, Gualdi S, Navarra A, et al. Paramount impact of the Indian Ocean dipole on the East African short rains: A CGCM study. J Climate. 2005;18: 4514–4530.

104. Mistry VV, Conway D. Remote forcing of East African rainfall and relationships with fluctuations in levels of Lake Victoria. Int J Climatol. 2003;23: 67–89.

105. Mutai CC, Ward MN, Colman AW. Towards the Prediction of the East Africa short rains based on sea-surface temperature-atmosphere coupling. Int J Climatol. 1998;18: 975–997.

106. Reverdin G, Cadet DL, Gutzler D. Interannual displacements of convection and surface circulation over the equatorial Indian Ocean. Quart J Roy Meteor Soc. 1986;112: 43–67.

107. Fischer EMA, Knutti R. Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes. Nat Clim Chang. 2015;5: 560–564.

108. Lu J, Vecchi GA, Reichler T. Expansion of the Hadley cell under global warming. Geophys Res Lett. 2007;34: L06805.

109. Vecchi GA, Soden BJ, Wittenberg AT, Held IM, Leetmaa A, Harrison MJ. Weakening of tropical Pacific atmospheric circulation due to anthropogenic forcing. Nature. 2006;441: 73–76. pmid:16672967

110. Owen-Smith N. Demography of a large herbivore, the Greater Kudu Tragelaphus Strepsiceros, in relation to rainfall. J Anim Ecol. 1990;59: 893–913.

111. Watson R, Turner M. A count of the large mammals of the Lake Manyara National Park: results and discussion. Afr J Ecol. 1965;3: 95–98.

112. Ogutu JO, Piepho H-P, Dublin HT. Ostrich recruitment dynamics in relation to rainfall in the Mara–Serengeti ecosystem. Ostrich. 2012;83: 119–136.

113. Prins HHT, Weyerhaeuser FJ. Epidemics in populations of wild ruminants: anthrax and impala, rinderpest and buffalo in Lake Manyara National Park, Tanzania. Oikos. 1987;49: 28–38.

114. Davies FG, Linthicum KJ, James AD. Rainfall and epizootic Rift Valley fever. Bull World Health Organ. 1985;63: 941–943. pmid:3879206

115. Baylis M, Mellor PS, Meiswinkel R. Horse sickness and ENSO in South Africa. Nature. 1999;397: 574–574. pmid:10050850

116. Waal C van der, Kool A, Meijer SS, Kohi E, Heitkönig IMA, de Boer WF, et al. Large herbivores may alter vegetation structure of semi-arid savannas through soil nutrient mediation. Oecologia. 2011;165: 1095–1107. pmid:21225433

117. Cooper SM, Holekamp KE, Smale L. A seasonal feast: long-term analysis of feeding behaviour in the spotted hyaena (Crocuta crocuta). Afr J Ecol. 1999;37: 149–160.

118. Fryxell JM. Aggregation and migration by grazing ungulates in relation to resources and predators. In: Sinclair ARE, editor. Serengeti II: dynamics, management and conservation of an ecosystem. Chicago, Illinois, USA and London, UK: The University of Chicago Press 1995. pp. 257–273.

119. Wangai P, Muriithi JK, Koenig A. Drought related impacts on local people’s socioeconomic life and biodiversity conservation at Kuku Group Ranch, southern Kenya. Int J Ecosyst. 2013;3: 1–6.

120. Muoria PK, Muruthi P, Kariuki WK, Hassan BA, Mijele D, Oguge NO. Anthrax outbreak among Grevy’s zebra (Equus grevyi) in Samburu, Kenya. Afr J Ecol. 2007;45: 483–489.

121. Bett B, Kiunga P, Gachohi J, Sindato C, Mbotha D, Robinson T, et al. Effects of climate change on the occurrence and distribution of livestock diseases. Prev Vet Med. 2017;137: 119–129. pmid:28040271