Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Year-round breeding equatorial Larks from three climatically-distinct populations do not use rainfall, temperature or invertebrate biomass to time reproduction

  • Henry K. Ndithia ,

    hndithia@gmail.com

    Affiliations Ornithology Section, Department of Zoology, National Museums of Kenya, Nairobi, Kenya, Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, the Netherlands

  • Kevin D. Matson,

    Current address: Resource Ecology Group, Department of Environmental Sciences, Wageningen University, Droevendaalsesteeg 3a, Wageningen, The Netherlands

    Affiliation Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, the Netherlands

  • Maaike A. Versteegh,

    Affiliation Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, the Netherlands

  • Muchane Muchai,

    Current address: Department of Clinical Studies (Wildlife and Conservation), College of Agriculture and Veterinary Sciences, University of Nairobi, Nairobi, Kenya.

    Affiliation Ornithology Section, Department of Zoology, National Museums of Kenya, Nairobi, Kenya

  • B. Irene Tieleman

    Affiliation Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, the Netherlands

Abstract

Timing of reproduction in birds is important for reproductive success and is known to depend on environmental cues such as day length and food availability. However, in equatorial regions, where day length is nearly constant, other factors such as rainfall and temperature are thought to determine timing of reproduction. Rainfall can vary at small spatial and temporal scales, providing a highly fluctuating and unpredictable environmental cue. In this study we investigated the extent to which spatio-temporal variation in environmental conditions can explain the timing of breeding of Red-capped Lark, Calandrella cinerea, a species that is capable of reproducing during every month of the year in our equatorial east African study locations. For 39 months in three climatically-distinct locations, we monitored nesting activities, sampled ground and flying invertebrates, and quantified rainfall, maximum (Tmax) and minimum (Tmin) temperatures. Among locations we found that lower rainfall and higher temperatures did not coincide with lower invertebrate biomasses and decreased nesting activities, as predicted. Within locations, we found that rainfall, Tmax, and Tmin varied unpredictably among months and years. The only consistent annually recurring observations in all locations were that January and February had low rainfall, high Tmax, and low Tmin. Ground and flying invertebrate biomasses varied unpredictably among months and years, but invertebrates were captured in all months in all locations. Red-capped Larks bred in all calendar months overall but not in every month in every year in every location. Using model selection, we found no clear support for any relationship between the environmental variables and breeding in any of the three locations. Contrary to popular understanding, this study suggests that rainfall and invertebrate biomass as proxy for food do not influence breeding in equatorial Larks. Instead, we propose that factors such as nest predation, female protein reserves, and competition are more important in environments where weather and food meet minimum requirements for breeding during most of the year.

Introduction

Ecological and environmental factors, such as food availability and weather, shape reproductive decisions in many bird species. These factors can act alone or in combination, and they may fluctuate in predictable or unpredictable ways within and between years. For birds living at mid-latitude locations in the temperate zone, predictable seasonal changes in day length and other environmental conditions function as reproductive cues. These birds can potentially time reproduction with regularity [1, 2, 3, 4]. Increasing photoperiod, a characteristic of temperate zone spring, triggers reproduction in birds via neuroendocrine mechanisms [5, 4, 6]. Food availability and temperature serve as supplementary cues to fine-tune the timing of reproduction to local environmental conditions [7, 8, 5, 4]. The net result is synchronized spring breeding within and among species living at the same location [7, 5].

In contrast, birds living in equatorial locations experience little predictable intra-annual variation in day length [9, 10], but instead experience large, and frequently unpredictable, variation in rainfall and food availability [11, 12, 13, 10, 14, 15, 16, 17, 18]. By living in environments without predictable seasonal cues, equatorial birds are thought to time reproduction based on shorter-term and more irregular factors, such as rainfall and food availability [19, 15]. Additionally, these birds tend to have more flexible breeding schedules and may breed opportunistically [12, 10]. For example, initiation of breeding with the onset of rain [19, 20, 21] greatly promotes nesting success in low-latitude birds [22, 23, 24]. Likewise in some equatorial birds, timing of reproduction coincides with peaks in food supply [19, 21, 25, 2, 14, 17]. These observations that in unpredictable equatorial environments birds preferentially breed at times of the year with higher rainfall and food, match with the general pattern that environments that are more arid have lower primary productivity and select for reduced reproductive effort [26, 27]. Other factors, such as wind and mist, both of which effectively lower ambient temperature, may also be important [28, 29].

One possible consequence of commencing breeding in response to unpredictable localized conditions is that a single species living in distinct environments might show variation in breeding patterns on a small geographical scale [17]. Exploiting such small-scale variation in environmental aridity within the tropics, we intensively investigated year-round breeding in three equatorial populations of Red-capped Larks Calandrella cinerea (Gmelin 1789) in Kenya for 39 consecutive months. These resident populations, despite their close geographic proximity, experience different patterns of temperature and rainfall, with climates ranging from warm and dry to cool and wet, and representing an expected gradient of increasing primary productivity [30]. Thus, the study system allows for a comparative, intraspecific analysis of environmentally induced spatio-temporal variation in reproduction. Rarely have studies assessed year-round breeding activities of equatorial species and related breeding to biotic and abiotic characteristics of climatically-distinct locations.

The overall objective of our study was to compare and understand breeding in Red-capped Larks in relation to spatio-temporal variation in weather conditions and food resources. Specifically, we 1) compared spatial variation in rainfall, temperature, invertebrate biomass and breeding across our three study locations, 2) described within-location year-round patterns of rainfall, temperature and invertebrate biomass and how these variables co-vary with breeding and, 3) determined which, if any, biotic and abiotic factors are related to occurrence and intensity of breeding by Red-capped Larks in each location. We predicted that the drier and warmer the location, or the drier and warmer the time of the year, the lower the productivity of invertebrates and the lower the intensity of breeding by Larks.

Materials and methods

Study system

Red-capped Larks are small grassland birds that are widely distributed across Africa. They prefer habitats dominated by short grasses or almost bare ground, including fallow and cultivated agricultural fields. Red-capped Larks feed mostly on invertebrates (including beetles, wasps, caterpillars, butterflies and moths, earthworms, and grasshoppers) and occasionally on grass seeds (pers. obs.). Pairs build ground-level open-cup nests that are placed next to a scrub or grass tuft. They typically lay two eggs per clutch (mean 1.89 ± 0.33 (SD) eggs, n = 279, range 1–3 eggs; pers. obs.). During breeding, birds defend the area around the nest but neighboring nests can be as close as 10 m; outside breeding they occur in flocks (pers. obs.). Before our study, nothing had been documented about timing, number of breeding attempts and other breeding parameters at the individual or population level.

From January 2011 to March 2014, we worked simultaneously in multiple plots in South Kinangop, North Kinangop and Kedong (see Table 1 for details per plot), three locations in central Kenya with distinct climates. Distances between locations are 19 km (South Kinangop—North Kinangop), 29 km (South Kinangop—Kedong) and 34 km (North Kinangop–Kedong). Accessible plots within locations were chosen based on observations of Red-capped Larks made by local bird watchers and by us (H.K.N., B.I.T.). We set up a weather station (Alecto WS-3500, Den Bosch, Netherlands) at each location (Table 1) to measure daily rainfall (mm) and minimum (Tmin) and maximum (Tmax) temperatures (°C). Using these measurements from three full calendar years (March 2011 –February 2014), we calculated annual and monthly rainfall, and annual and monthly Tmin and Tmax.

thumbnail
Table 1. Coordinates, altitude (m ASL), surface area (km2) and distance to weather station for each plot in our three study locations South Kinangop, North Kinangop and Kedong.

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

South and North Kinangop lie on a plateau of montane grassland along the Aberdare mountain ranges. Study plots in South and North Kinangop flood periodically during rains and standing water remains after rains have stopped (pers. obs. 2010–2014). In South Kinangop, birds bred only in Seminis, despite initially observing them also in the other two plots (Table 1). Flooding made Seminis unavailable for breeding from April–December 2012 and April 2013. Flooding in North Kinangop affected nests located in parts of Joshua and Ndarashaini in October 2011 and October 2012; these plots also received heavy rainfall in April 2013 that affected nesting activities. Kedong, a privately owned ranch in the Rift Valley in Naivasha, consists of large grassland patches that did not flood (pers. obs. 2010–2014).

The study species involved is not and endangered or protected species. The National Museums of Kenya approved this research and owners of the land gave permission to conduct the study on their respective sites.

Estimating invertebrate availability

To assess invertebrate biomass as proxy for food availability in each location, we used pitfall traps and sweep nets to sample ground and flying invertebrates respectively once per month [31]. We assessed within-location variation in invertebrate biomass using data from plots within a location. For pitfall traps, we used plastic cups with a 26 cm circumference that contained ≈100 ml of 5% formaldehyde solution and that we buried so that rims of the cups were at ground level. We placed five pitfall traps in each plot and left them in place for five days each month. We placed the traps 70 m apart along a 280m-long transect in all plots except one plot in South Kinangop (Seminis), where instead we equally spaced 10 pitfall traps along a 630m-long transect. For sweep-netting (net diameter 0.4m), we established permanent 50m long transects in each plot, subjectively selected as representative for the plot. Per location, one field assistant collected invertebrates between 9:00–10:00 am. If it rained during this hour, we postponed sweep-netting to the same hour on a day without rain. All field assistants were trained to sample in the same manner. We standardized the analyses of pitfall and sweep net sampling data per location (see statistical analysis section below). To calculate annual average and monthly average biomasses, we used two complete calendar years (24 months, March 2011-February 2012 and March 2013-February 2014), excluding the year in which flooding caused multi-month gaps in the data (March 2012-February 2013). Our resulting 24-month data sets had three missing values for ground invertebrates (October 2013, North Kinangop; October and December 2011, Kedong) and three missing values for flying invertebrates (October 2011, October 2013, North Kinangop; February 2014, South Kinangop). For calculations of annual and monthly averages, we substituted each missing value with the average value of the preceding and subsequent months.

We preserved collected specimens in 70% alcohol, later identifying and sorting them based on morphology [32]. For biomass estimation, we classified invertebrates into 10 categories based on size and shape: ants; bees and wasps; beetles and bugs; butterflies and moths; caterpillars, caddisflies, and stoneflies; diplura, millipede, centipede, and earthworms; flies; grasshoppers, crickets, and mantis; spiders, ticks, and mites; and the rest (woodlice, cicadas, cockroaches and earwigs).

We estimated biomass of each of our invertebrate categories as a proxy for food availability. To do this, we first used a subsample of 2198 invertebrate specimens, representing all invertebrate categories from all locations, to develop a category-specific calibration curve relating dry mass as a function of length and width [31, 33, 34]. For every individual in the subset, we measured length (anterior-most part of the head to the tip of abdomen) and width (the widest point of abdomen) using vernier calipers, dried them in an oven for 48 hours at 65°C [35, 33, 34], and measured dry mass on an analytical balance (model KERN ACS 220-4N, KERN and Sohn of Belingen, Germany). We used a log-transformed power model to describe the length-width-mass relationship; the power model has been shown to give the highest adjusted r2 compared to length-mass and length-area relationships [35, 33, 34]:

biomass = a + b log(length) + c log (width), where a, b and c are coefficients of the model from each of the invertebrate categories whose biomass we estimated

We used calibration curves per invertebrate category to predict body mass from length and width (for details on the adjusted r2 and the range of length and width, see S1 Appendix).

Overall, we collected, measured, dried, and applied the biomass estimation protocol to 23,628 specimens from pitfall traps and 3260 captured by sweep-netting (including calibration subset).

Lark reproduction

To determine the year-round breeding activities of Larks, we spent on average 134 person-hours per month searching for nests in the three locations combined (Table 2 provides a breakdown for effort per location). Our nest search strategy included observing breeding behavior (e.g., transport of nest materials or food, breeding-related alarm calls, nervous parental behavior around nest sites) and routinely walking plots to flush parents incubating eggs or brooding young. We quantified nest-searching effort as person-hours, i.e., number of hours searching for nests multiplied by the number of persons searching. For each month we calculated a “nest index” (i.e., level of breeding) value, which we defined as the total number of nests found in a month per 10 person-hours of effort. To calculate annual average and monthly average nest indices we used the two complete calendar years (24 months) of March 2011-February 2012 and March 2013-February 2014, excluding the year in which flooding caused multi-month gaps in the data (March 2012-February 2013). Our resulting 24-month data set had one missing value (April 2013, South Kinangop), for which we substituted the average of March and May 2013.

thumbnail
Table 2. Search effort (in days (days had a minimum of 2 searching hours) and hours per month) for nests of Red-capped Larks Calandrella cinerea in our three study locations South Kinangop, North Kinangop and Kedong, from January 2011 to March 2014.

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

Statistical analyses

For all analyses, we tested and confirmed that the dependent variable and the final models observed the assumptions of normality and homoscedasticity of variance through graphical and statistical methods. We tested for among-location differences in rainfall, Tmin and Tmax, ground and flying invertebrates, and nest index (continuous variable) using mixed models (R-package lme) with location as fixed factor and month as random factor. To compare invertebrate biomasses among plots and locations, we standardized ground and flying invertebrate sampling by expressing biomass per five pitfall traps and one sweep net session per plot or location per month. For among-location comparisons, we log-transformed ground and flying invertebrate data because they were not normality distributed. We found no significant among-plot differences in ground invertebrate biomasses within South Kinangop (F2, 47 = 0.89, P = 0.42) or in Kedong (X2 = 3.98, P = 0.26). For North Kinangop, among-plot differences in ground invertebrate biomasses were significant (X2 = 6.49, P = 0.04), although post-hoc tests showed no significant differences among-plots. There were no significant among-plot differences in flying invertebrate biomasses for any of the three locations (all X2 < 5.74, P > 0.06). Therefore, we used the mean monthly biomass per location to test for among-location differences.

We investigated if and how environmental conditions in the month before breeding (“prior”) and in the month of breeding (“current”) were associated with the occurrence and intensity of breeding. We calculated pairwise correlation coefficients between the environmental factors per location (supplementary material 2), to identify potential collinearity. We used model selection based on the Akaike Information Criterion corrected for sample size (AICc) because this allows for exploration of multiple models simultaneously (Burnham and Anderson 2002). To investigate what determines occurrence of breeding, we transformed the continuous variable nest index into the new variable “occurrence of breeding” (binomial: presence/absence). We then used generalized linear mixed models with a binomial distribution to construct for each location a “full” model with the new dependent variable “occurrence of breeding”. These full models included ten explanatory variables (i.e., prior and current rain, Tmin, Tmax, ground invertebrate biomass, and flying invertebrate biomass) and four two-way interactions, i.e., prior and current rain and corresponding ground invertebrate biomass and prior and current rain and corresponding flying invertebrate biomass. We compared all the possible models and ranked them in order of their AICc, such that the lowest values were considered to have more statistical power [36]. The model with the highest weight and the lowest AICc value was considered the most parsimonious, although all models within 2 AICc of the best model were included in further analysis (Grueber et al. 2011). We explored the relative contribution of the various environmental parameters to breeding by applying model averaging and standardization based on all models with ΔAICc values < 2 (the “best model-set”), compared with the top model (Grueber et al. 2011). Although AICc values of models in the best model-set without month as random effect were higher than models with random effect, we added month as a random effect to the models to correct for potential seasonal effects. We analyzed all data using R statistical software (version 3.0.3) [37].

In the second part of our analysis we investigated how environmental conditions in the month before breeding and in the month of breeding were associated with the intensity of breeding. For this, we analyzed only the months in which breeding occurred (i.e. nest index > 0). We used linear models with a Gaussian distribution and constructed a “full” model with continuous variable “nest index”, for each location. We used the same explanatory variables and statistical approach as in the analysis above. Because month never improved the models in the analysis of occurrence of breeding (see above), and in order to maximize power for the tests of the effects of environmental factors we did not include month as a random effect in these models. In addition, because of low sample size in South Kinangop (n = 7 months), we performed the analyses of intensity of breeding only for North Kinangop and Kedong.

Additional analyses, in which we explored the effects of different time windows and time lags of the environmental variables (up to six months preceding breeding) on breeding occurrence and intensity, did not result in qualitatively different results (see supplementary material 1).

Results

Spatial differences in environmental factors, invertebrates and breeding of Larks

Mean annual and monthly rainfall were highest in South Kinangop, intermediate in North Kinangop and lowest in Kedong (Table 3; Figs 1A, 2A and 3A). South Kinangop received on average 123% more rain than Kedong, while North Kinangop received 40% more rain than Kedong. Tmin and Tmax were lowest in South Kinangop (annual average Tmin = 5.5°C, annual average Tmax = 24.7°C), intermediate in North Kinangop (annual average Tmin = 9.1°C, annual average Tmax = 25.4°C) and highest in Kedong (annual average Tmin = 10.5°C, annual average Tmax = 28.6°C) (Table 3; Figs 1B, 2B and 3B).

thumbnail
Fig 1.

Temporal variation during January 2011-March 2014 of A. rainfall (mm) and nest index (number of nests/10 search h) of Red-capped Larks Calandrella cinerea, B. average monthly maximum (Tmax) and minimum (Tmin) temperature (°C), C. biomass (g dry weight) of ground-dwelling and flying invertebrates in South Kinangop. Horizontal open rectangles represent periods of flooding (i.e. standing water in the study location). Data gaps other than from flooding in ground and flying invertebrates represent missing data due to e.g., vandalism (see Methods).

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

thumbnail
Fig 2.

Temporal variation during January 2011-March 2014 of A. rainfall (mm) and nest index (number of nests/10 search h) of Red-capped Larks Calandrella cinerea, B. average monthly maximum (Tmax) and minimum (Tmin) temperature (°C), C. biomass (g dry weight) of ground-dwelling and flying invertebrates in North Kinangop. Horizontal open rectangles represent periods of flooding (i.e. standing water in the study location). Data gaps other than from flooding in ground and flying invertebrates represent missing data due to e.g., vandalism (see Methods).

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

thumbnail
Fig 3.

Temporal variation during January 2011-March 2014 of A. rainfall (mm) and nest index (number of nests/10 search h) of Red-capped Larks Calandrella cinerea, B. average monthly maximum (Tmax) and minimum (Tmin) temperature (°C), C. biomass (g dry weight) of ground-dwelling and flying invertebrates in Kedong. Data gaps in ground and flying invertebrates represent missing data due to e.g., vandalism (see Methods).

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

thumbnail
Table 3. Annual (n = 3 years) and monthly (n = 36 months) rainfall (average ± SD, and range), and monthly minimum and maximum temperatures (n = 36 months, average ± SD, and range) as measured by our weather stations in South Kinangop, North Kinangop and Kedong, during March 2011 –February 2014.

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

Biomasses of ground invertebrates (log pitfall, mg) did not differ significantly among locations, but biomasses of flying invertebrates (log sweepnet, mg) were highest in Kedong, intermediate in South Kinangop, and lowest in North Kinangop (Table 4). Flying invertebrate biomasses were on average 42% lower in North Kinangop than in Kedong and 27% lower in South Kinangop than in Kedong; flying invertebrate biomasses did not differ significantly between South and North Kinangop (Table 4).

thumbnail
Table 4. Annual and monthly biomass of ground invertebrates (log(mg dry mass/5 pitfalls)) and flying invertebrates (log(mg dry mass/sweep-net session)) in South Kinangop, North Kinangop and Kedong.

Annual values (average ± SD) were based on two calendar years (March 2011-February 2012 and March 2013-February 2014). Data for the third year (March 2012-February 2013) were excluded because flooding caused incomplete data sets for South and North Kinangop (see Methods). Likewise, monthly values (average ± SD, range) were based on 24 months (March 2011- February 2012 and March 2013-February 2014).

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

Despite differences in climate and invertebrate biomass, Red-capped Larks bred in all three locations (Table 5). In the period January 2011 –March 2014, we found 74 nests in South Kinangop, 63 nests in North Kinangop and 153 nests in Kedong (Table 5). Calculating nest index corrected for search effort, we found the highest numbers in Kedong, followed by North Kinangop (63% lower than Kedong) and South Kinangop (84% lower than in Kedong).

thumbnail
Table 5. Total number of nests found, percentage successful nests, and annual and monthly nest index (number of nests/10 h search effort) in South Kinangop, North Kinangop and Kedong.

Total number of nests found is based on the period January 2011 –March 2014. Nests failed due to a variety of reasons, including predation, abandonment, flooding and human destruction. Average nest indices are based on two complete calendar years (March 2011—February 2012, and March 2013 –February 2014; 24 months) and excluded the year in which flooding occurred (see Methods). Superscripts indicate subsets of significant differences (P<0.05) among locations in post-hoc tests, after mixed model analyses.

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

Year-round variation in environmental conditions, invertebrates and breeding of Larks

In all three locations, rainfall occurred in all calendar months of the year, but the amount of rainfall in any given month was highly variable and unpredictable among years (Figs 1A, 2A and 3A). The only consistent annually recurring observation was that January and February were dry in all three locations in all four years (with the exception of Kedong in 2014). Outside of this annually recurring dry season, there was no month without rain in North Kinangop and only one month was without rain in South Kinangop (March 2014). However, Kedong received no rain at all during six months in 2013 (June-November), in contrast to 2011 and 2012 when this location received rain every month. Average monthly Tmin and average monthly Tmax varied unpredictably throughout the year in all locations and years, but generally, Tmin were lowest and Tmax were highest in January and February each year (Figs 1B, 2B and 3B). Average monthly Tmax varied between 21.2°C and 30.0°C in South Kinangop, between 22.1°C and 30.5°C in North Kinangop, and between 25.3°C and 34.9°C in Kedong (Table 3). Likewise, average monthly Tmin varied between 3.0°C and 8.2°C in South Kinangop, between 3.0°C and 13.7°C in North Kinangop, and between 6.2°C and 15.7°C in Kedong (Table 3).

Ground and flying invertebrates were present in all months in all locations, but biomasses varied among months and among years in an unpredictable manner (Figs 1C, 2C and 3C).

Overall, we observed Red-capped Larks breeding in all calendar months, but they did not breed in every month in every year in any of the three locations (Figs 1A, 2A and 3A). In South Kinangop, we found nests in the calendar months January-April and June-August, and in 10 out of 30 months total (33%). In North Kinangop, we found nests in all calendar months except June and July, and in 21 out of 39 months total (54%). Finally, in Kedong, we found nests in all calendar months and in 20 out of 39 months total (51%). In all locations, year-to-year variation in nest index was present, with highest nest indices in 2012 (Figs 1A, 2A and 3A).

Associations between nest index, environmental factors and invertebrate biomass

The abiotic (rainfall, Tmin and Tmax) and biotic (ground and flying invertebrate biomass) factors, in the month before or the month of breeding, explained little variation in occurrence or intensity of breeding of Larks, in any of the three locations (Table 6). For each location, and for both sets of analyses, no single best model emerged. For the analysis of occurrence of breeding there were subsets of six (South Kinangop), 13 (North Kinangop), and four (Kedong) models with ΔAICc < 2 (Table 6A). These models contained zero to two of the ten abiotic and biotic factors included in the analysis; each of the subsets contained four to seven of the ten factors (Table 6A). The explained variation (D2) of models was consistently low, and for South Kinangop and North Kinangop, the model sets included the intercept-only “null” model (i.e., no environmental factors; Table 6A).

thumbnail
Table 6. Model selection results of (A) occurrence of breeding of Red-capped Larks in South Kinangop, North Kinangop and Kedong, and (B) intensity of breeding in North Kinangop and Kedong as a function of biotic and abiotic factors in the month prior to the breeding observation and the month of the breeding observation (see Methods for details).

Parameters, degrees of freedom, AICc, ΔAICc, weights and the explained variation (D2 in (A) and R2 in (B)) of the top-model and models with a ΔAICc <2.

https://doi.org/10.1371/journal.pone.0175275.t006

Using model averaging and standardization to explore the relative contribution of the environmental parameters to breeding, we could not identify one or more significant environmental factors that explained breeding in Larks (Table 7A; Fig 4A–4D).

thumbnail
Fig 4. Average monthly (± SE) values during non-breeding and breeding months, for the biotic and abiotic factors that were selected in the model selection analysis (Table 6, Table 7) for South Kinangop (black symbols), North Kinangop (grey symbols), and Kedong (white symbols).

A. Minimum (Tmin) and maximum (Tmax) temperature (°C) in current month, B. Minimum (Tmin) and maximum (Tmax) temperature (°C) in prior month, C. Flying invertebrate biomass (log(g dry weight)) in prior month. D. Ground invertebrates biomass (log(g dry weight)) in current month.

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

thumbnail
Table 7. Model averaging and standardization results of occurrence of breeding of Red-capped Larks in South Kinangop, North Kinangop and Kedong, as a function of biotic and abiotic factors in the month prior to the breeding observation and the month of the breeding observation (see Methods for details).

Standardized estimates (± SE, 95% confidence intervals) for biotic and abiotic factors of the average model per location are based on the model subsets in Table 6.

https://doi.org/10.1371/journal.pone.0175275.t007

For the analysis of intensity of breeding there were subsets of six (North Kinangop) and three models (Kedong) with AICc < 2 (Table 6B). These models contained zero to two of the ten abiotic and biotic factors included in the analysis; each of the subsets contained two to seven of the ten factors (Table 6B). The explained variation (R2) of models was consistently low, and for North Kinangop, the model sets included the intercept-only “null” model (i.e., no environmental factors; Table 6B). Using model averaging and standardization to explore the relative contribution of the environmental parameters to breeding, we could only identify one significant environmental factor that explained breeding in Larks (i.e. Tmax in North Kinangop; Table 7B).

Discussion

This study showed that year-round breeding activities of Red-capped Larks in three climatically-distinct equatorial populations were not associated with rainfall, temperature and invertebrate biomass, across and within locations. Across locations that represent a gradient of rainfall and temperature, we found no support for the prediction that drier and warmer locations had lower invertebrate biomass and less breeding activity of Larks. In line with these results, within each location we also found no evidence that breeding was timed to co-occur with rainfall, higher or lower temperatures or invertebrate biomass. Instead, we observed highly unpredictable and irregular variation in environmental variables, invertebrate biomass and breeding of Larks, among months and among years. Red-capped Larks bred in all calendar months overall, but they did not breed in every month in every year in every location. Our findings raise the question of which factors do trigger equatorial Red-capped Larks to initiate breeding. They suggest that environmental conditions vary at a small spatio-temporal scale and often provide minimum requirements for breeding, at least for some individuals in each population.

Despite the small geographical distances (19–34 km) separating the three locations, our weather data allowed us to quantitatively confirm the distinct spatial patterns in rainfall and temperatures that we expected among locations due to orography (Table 1). However, against our predictions, the warmest and driest location, Kedong, had the highest flying invertebrate biomasses and the highest Lark nest index. In contrast, the cool and wet climatic extreme, South Kinangop, had intermediate ground and flying invertebrate biomasses and the lowest Lark nest index. Overall, the general global trend of lower primary productivity in more arid environments [38, 30, 27] is not reflected in our invertebrate and Lark breeding data.

Other factors that, in combination with rainfall and temperature, may affect Lark ecology and that differ among locations include soil type, land use, and the occurrence of climatic excesses. South Kinangop’s climate and fertile soil allow for intensive crop cultivation, but the combination of heavy rains and poor soil drainage also make the area prone to flooding. Seminis, a study site in South Kinangop with communally grazed land, experiences multiple months of standing water in some years, which might negatively affect productivity and diversity of plants and invertebrates [39]. In contrast, Kedong’s grasslands are extensively grazed by diverse wild and domestic herbivores, a process that might increase primary productivity and favor invertebrate and Lark reproduction [40, 41, 27, 42].

All three locations showed a lack of predictable intra-annual patterns in rainfall, temperature and invertebrate biomasses. Lark breeding was most unpredictable in Kedong where we observed breeding in 51% of the months but in all 12 calendar months, and most predictable in South Kinangop, where breeding occurred in 33% of the months restricted to seven calendar months. In addition, variation among years contributed to the unpredictability in all locations. Although we did not identify relationships among rainfall and breeding during the 39 months of our study, we did observe that breeding activities were affected by multi-month rainfall excesses: during a six-month drought in Kedong breeding was almost (but not fully) absent, while extended flooding prevented birds from breeding in South Kinangop and destroyed nests in both South and North Kinangop. These observations are in line with other studies in highly variable and unpredictable environments that also show that flooding from heavy rainfall disrupts breeding by destroying nests, by directly killing eggs and young [28, 43], and by reducing food availability or foraging efficiency [28, 29, 10, 39]. Moreover, drought can result in food and foraging limitations, and birds may avoid breeding during dry seasons [12, 20, 44, 10, 45, 46, 23].

The results of our model selection analysis, which did not strongly identify any of the environmental factors as relevant for Lark breeding, have two potential explanations. Either the spatial scale at which we sampled rainfall, temperature, and invertebrates was not the scale that Red-capped Larks used to time breeding, or Red-capped Larks used other factors to initiate breeding. We sampled rainfall and temperature with one weather station per location, placed centrally among the multiple plots within a location but as a result also at varying distance to each plot (Table 1). Likewise, we sampled invertebrates at only one transect in each plot. One might wonder if the small-scale spatio-temporal variation in environmental and ecological factors that we discovered calls for measuring rainfall, temperature and invertebrates at the level of Lark territories. While the addition of territory-level data might be interesting, we do not believe that sampling scale underlies our population-level findings about Lark breeding. Moving around in small flocks when not breeding, Red-capped Larks do not appear to stay within their breeding territory year-round (pers. obs. based on color-ringed individuals), a behavior that contrasts with some other Lark species (e.g., Hoopoe Larks, Alaemon alaudipes, in the Arabian Desert [45].

We therefore propose that Red-capped Larks use other factors to time reproduction, with prime candidates being nest predation [47, 48, 49, 50, 51], female protein reserves [52, 53, 54, 55], or social factors. Although a detailed study remains to be done, nest predation in our study sites is generally high, with only 53 nests that fledged from the total of 290 nests found in different stages of the nesting cycle (pers. obs.). In the face of such intense nest predation Red-capped Larks may breed opportunistically [11, 56, 57, 45, 42]. If environmental conditions are always permissive, birds may breed whenever they have resupplied their reserves after a failed nest attempt [58]. This would also be in line with earlier studies on equatorial passerines with opportunistic breeding schedules, in which the authors suggested that protein reserves of individual females may determine whether and when they breed, leading to asynchronous year-round breeding activities at the population level [52, 53, 54, 55, 56]. Social factors can affect breeding decisions in multiple ways. For example, some individuals may time their breeding to benefit from peak food availability [19, 25, 2, 14], while others may avoid competition for food or nesting space by conspecifics or other species by choosing “unpopular” times between food peaks [59]. In addition, timing of breeding may be influenced by predation-avoidance strategies (i.e., a pair avoids being the only one breeding at a particular time due to the high predation risk associated with that position [60, 61] and prolonged nestling dependence on parents [62]. A better understanding of these possible mechanisms will come from studies at the individual level that complement our population level findings.

Supporting information

S1 Appendix. Number of samples (n) and range of length (minimum L in mm, maximum L in mm) and width (minimum W in mm, maximum W in mm) of invertebrates used in generating calibration curves, and statistics per curve (coefficients a, b, c; adjusted r2, degrees of freedom, F and P-values), per invertebrate category we used to predict body mass from length and width.

Calibration curves were fit with the formula log(biomass in mg dry weight) = log(a) + b log(length in mm) + c log (width in mm), where a, b and c are coefficients of the model (see Methods).

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

(DOCX)

S1 File. Outcome of the investigation of the effect of maximum temperature, minimum temperature, rainfall, ground and flying invertebrates on nest index at time windows and time lags up until 6 month prior to the current month in Red-capped Larks living in three climatically-distinct Kenya locations, South Kinangop, North Kinangop and Kedong during March 2011 –February 2014.

Please note the inclusion of Table A and Figure A in S2 File.

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

(DOCX)

S1 Table. Correlation coefficients (below the diagonal) and P-values (above the diagonal) of pairwise correlations among the 3 weather variables and 2 invertebrate variables in three Kenya locations, South Kinangop, North Kinangop and Kedong during March 2011 –February 2014.

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

(DOCX)

S2 Table. Coefficients of variation of monthly (n = 36 months) rainfall, minimum temperature and maximum temperatures as measured by our weather stations in South Kinangop, North Kinangop and Kedong, during March 2011 –February 2014.

F-tests of equal variances indicated that there were no significant differences among locations (P<0.05; indicated by superscripts).

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

(DOCX)

Acknowledgments

We thank the landowners in South Kinangop (Maina Irungu, Nairobi Water and Sewage Company, and local authorities of Seminis) and North Kinangop (Kimani Mbae, Isaac Gathitu, Joshua Kimani and Francis Kagai) and the owners and manager (Amos Omondi) of Kedong Ranch for allowing us to carry out research on their land. We thank our field assistants Abraham Mwangi Kuria, Paul Maina Kimani, and Peter Kinyanjui Gachigi for field data collection and Dr. N.P.C. Horrocks for reading our manuscript and offering valuable input. Dr. A. van der Plas contributed ideas about the analysis of invertebrate biomass data. Sarah Higgins of Lake Naivasha Riparian Association provided a base to the research team during the years of fieldwork.

This study was funded by The Netherlands Fellowship Programme of Nuffic (grant No. CF6833/2010 to B.I.T and H.K.N.), the Netherlands Organization for Scientific Research (NWO-VIDI to B.I.T.), The Young Academy project grant (B.I.T.), the University of Groningen, the Schure-Beijerinck-Popping Fonds, and Dr. J.L. Dobberke foundation. The National Museums of Kenya provided paid study leave and organized permission letters for access to study areas. Capture and handling of birds was done by skilled bird ringers in adherence to the animal welfare protocols of National Museums of Kenya.

Author Contributions

  1. Conceptualization: BIT HKN.
  2. Data curation: HKN.
  3. Formal analysis: HKN BIT MAV KDM.
  4. Funding acquisition: HKN BIT.
  5. Investigation: HKN BIT MAV.
  6. Methodology: HKN BIT.
  7. Project administration: HKN.
  8. Resources: HKN BIT.
  9. Supervision: BIT MAV KDM MM.
  10. Writing – original draft: HKN.
  11. Writing – review & editing: HKN BIT MAV KDM MM.

References

  1. 1. Colwell RK. Predictability, constancy and contingency of periodic phenomena. Ecology. 1974; 55():1148–1153
  2. 2. Hau M, Wikelski M, Wingfield JC. Visual and nutritional food cues fine-tune timing of reproduction in a Neotropical rainforest bird. Journal of Experimental Zoology. 2000; 286: 494–504 pmid:10684573
  3. 3. Wikelski M, Hau M, Wingfield JC. Seasonality of reproduction in a Neotropical rain forest bird. Ecology 2000; 81: 2458–2472
  4. 4. Hau M. Timing of breeding in variable environments: tropical birds as model systems. Hormones and behavior. 2001; 40: 281–290 pmid:11534993
  5. 5. Lambrechts MM, Perret P, Blondel J. Adaptive differences in the timing of egg laying between different populations of birds results from variation in photoresponsiveness. Proceedings of the Royal Society of London. 1996; 263: 19–22
  6. 6. Gwinner E. Circannual rhythms in birds. Current Opinions in Neurobiology. 2003; 13: 770–778
  7. 7. Nager RG, van Noordwijk AJ. Proximate and ultimate aspects of phenotypic plasticity in timing of great Tit breeding in a heterogeneous environment. The American Naturalist. 1995; 146: 454–474
  8. 8. Van Noordwijk AJ, Mccleery RH, Perrins CM. Selection for the timing of great Tit breeding to caterpillar growth and temperature. Journal of Animal Ecology. 1995; 64: 451–458
  9. 9. Moreau RE. The breeding seasons of African birds–Landbirds. Ibis. 1949; 92: 223–267
  10. 10. Young BE. The effects of food, nest predation and weather on the timing of breeding in tropical house Wrens. The Condor. 1994; 96: 341–353
  11. 11. Grant PR, Boag PT. Rainfall on the Galápagos and the demography of Darwin’s Finches. Auk. 1980; 97: 227–244
  12. 12. Boag PT, Grant PR. Darwin’s finches (Geospiza) on Isla Daphne major, Galapagos: breeding and feeding ecology in a climatically variable environment. Ecological Monographs. 1984; 54: 463–489
  13. 13. Wrege PH, Emlen ST. Breeding seasonality and reproductive success in white-fronted Bee-eaters in Kenya. Auk. 1991; 108: 673–687
  14. 14. Scheuerlein A, Gwinner E. Is food availability a circannual zeitgeber in tropical birds? A field experiment on Stonechat in tropical Africa. Journal of Biological Rhythms. 2002; 17: 171–180 pmid:12002164
  15. 15. Stutchbury BJM, Morton ES. Behavioral ecology of tropical birds. San Diego: Academic Press; 2001
  16. 16. Conway D, Allison E, Felstead R, Goulden M. Rainfall variability in east Africa: implication for natural resources management and livelihoods. Philosophical Transactions of the Royal Society A. 2005; 363: 49–54
  17. 17. Moore IT, Bonier F, Wingfield JC. Reproductive asynchrony and population divergence between two tropical bird populations. Behavioral Ecology. 2005; 16:755–762
  18. 18. Ndithia H, Perrin MR, Waltert M. Breeding biology and nest site characteristics of the Rosy-faced Lovebird Agapornis roseicollis in Namibia. Ostrich. 2007; 78:13–20
  19. 19. Dittami JP, Gwinner E. Annual cycle in the African stonechat Saxicola torquata rubicola and their relationship to environmental factors. Journal of Zoological Society of London. 1985; 207: 357–370
  20. 20. Dittami JP. Seasonal reproduction, moult and their endocrine correlates in two tropical Ploceidae species. Journal of Comparative Physiology B. 1986; 156: 641–647
  21. 21. Chapman A. Breeding and moult of four bird species in tropical West Africa. Tropical Zoology. 1995; 8: 227–238
  22. 22. Lloyd P, Little RM, Crowe TM, Simmons RE. Rainfall and food availability as factors influencing the migration and breeding activity of Namaqua Sandgrouse Pterocles namaqua. Journal of African Ornithology. 2001; 72:50–62
  23. 23. Lepage D, Lloyd P. Avian clutch size in relation to rainfall seasonality and stochasticity along an aridity gradient across South Africa. Journal of African Ornithology. 2004; 75: 259–268
  24. 24. Monadjem A, Bamford AJ. Influence of rainfall and timing and success of reproduction in Marabous Stork Leptoptilos crumeniferus. Ibis. 2009; 151:344–351
  25. 25. Komdeur J. Seasonal timing of reproduction in a tropical bird, the Seychelles Warbler: a field experiment using translocation. Journal of Biological Rhythms. 1996; 11: 333–346 pmid:8946261
  26. 26. Tieleman BI, Williams JB, Bloomer P. Adaptation of metabolism and evaporative water loss along an aridity gradient. Proceedings of the Royal Society London B. 2003a; 270: 207–214
  27. 27. Tieleman BI, Williams JB, Visser GH. Energy and water budgets of Larks in a life-history perspective: parental effort varies with aridity. Ecology. 2004; 85: 1399–1420
  28. 28. Foster MS. Rain, feeding behavior and clutch size in tropical birds. Auk. 1974; 91: 722–726
  29. 29. Tye H. Reversal of breeding season by lowland birds at higher altitudes in western Cameroon. Ibis. 1991; 114:154–161
  30. 30. Peñuelas J, Prieto P, Beier C, Cesaraccio C, De Angelis P, De Dato G, et al. Response of plant species richness and primary productivity in shrublands along a north–south gradient in Europe to seven years of experimental warming and drought: reductions in primary productivity in the heat and drought year of 2003. Global Change Biology. 2007; 13:2563–2581
  31. 31. Ganihar SR. Biomass estimates of terrestrial invertebrates based on body length. Journal of Bioscience. 1997; 22: 219–224
  32. 32. Picker M, Griffiths C, Weaving A. Field guide to insects of S. Africa. Struik Publishers, South Africa; 2004
  33. 33. Benke AC, Huryn AD, Smock LA, Wallace JB. Length-mass relationships for freshwater macroinvertebrates in North America with particular reference to the Southeastern United States. Journal of North American Benthological Society. 1999; 18:308–343
  34. 34. Gruner DS. Regressions of length and width to predict invertebrate biomass in the Hawaiian Islands. Pacific Science. 2003; 57:325–336
  35. 35. Sample BE, Cooper RJ, Greer RD, Whitmore RC. Estimation of insect biomass by length and width. American Midland Naturalist. 1993; 120: 234–240
  36. 36. Burnham KP, Anderson DR. Model selection and multimodel inference: A practical information-theoretic approach, 2nd ed. Springer-Verlag, New York; 2002
  37. 37. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria; 2014
  38. 38. Pérez-Sánchez AJ, Lattke JE, Viloria AL. Patterns of Ant (Hymenoptera: Formicidae) richness and relative abundance along and aridity gradient in western Venezuela. Neotropical Entomology. 2013; 42: 128–136 pmid:23949745
  39. 39. Mcmullen LE, Lytle DA. Quantifying invertebrate resistance to floods: a global-scale meta-analysis. Ecological Applications. 2012; 22: 2164–2175 pmid:23387117
  40. 40. Veen GFC, Olff H, Duyts H, van der Putten WH. Vertebrate herbivores influences soil nematodes by modifying plant communities. Ecology. 2010; 91: 828–835 pmid:20426340
  41. 41. Van Klink R, van der Plas F, van Noordwijk CGET, WallisDeVries MF, Olff H. Effects of large herbivores on grassland arthropod diversity. Biological Reviews. 2015; 90: 347–366 pmid:24837856
  42. 42. Tieleman BI, Williams J. To breed or not to breed: that is the question, decisions facing hoopoe Larks in the Arabian desert. Seeking Nature’s Limits. KNNV publishing, Utrecht; 2005
  43. 43. Poiani A. Effects of floods on distribution and reproduction of aquatic birds. Advances in Ecological Research. 2006; 9: 63–83
  44. 44. Poulin B, Lefebvre G, Mcneil R. Tropical avian phenology in relation to abundance and exploitation of food resources. Ecology. 1992; 73: 2295–2309
  45. 45. Tieleman BI, Williams JB. Effects of food supplementation on behavioural decisions of hoopoe-larks in the Arabian Desert: balancing water, energy and thermoregulation. Animal Behaviour. 2002; 63: 519–529
  46. 46. Tieleman BI, Williams JB, Visser GH. Variation in allocation of time, water and energy in Hoopoe Larks from the Arabian Desert. Functional Ecology. 2003b; 17: 869–876.
  47. 47. Martin TE. Avian life history evolution in relation to nest sites, nest predation and food. Ecological Monographs. 1995; 65: 101–127
  48. 48. Ghalambor CK, Martin TE. Comparative manipulation of predation risk in incubating birds reveals variability in the plasticity of responses. Behavioral Ecology. 2002; 13: 101–108
  49. 49. Decker KL, Conway CJ, Fontaine JJ. Nest predation, food and female age explain seasonal declines in clutch size. Evolutionary Ecology. 2012; 26:983–699
  50. 50. Praus L, Hegemann A, Tieleman BI, Weidinger K. Predators and predation rates of Skylark Alauda arvensis and Woodlark Lullula arborea nests in a semi-natural area in The Netherlands. Ardea. 2014; 102: 87–94.
  51. 51. Toyama M, Kotaka N, Koizumi I. Breeding timing and nest predation rate of sympatric Scops owls with different dietary niche breadth. Canadian Journal of Zoology. 2015; 93: 841–847
  52. 52. Ward P. The annual cycle of yellow-vented Bulbul Pycnonotis goiavier in humid equatorial environment. J. Zool. Lond. 1969; 157: 25–45
  53. 53. Fogden MPL. The seasonality and population dynamics of equatorial forest birds in Sarawak. Ibis. 1972; 114: 307–343
  54. 54. Jones PJ, Ward P. The level of reserve protein as the proximate factor controlling the timing of breeding and clutch-size in Red-billed Quelea Quelea quelea. Ibis. 1976; 118: 547–574
  55. 55. Fogden MPL, Fogden PM. The role of fats and protein reserves in the annual cycle of Grey-backed Camaroptera in Uganda (Aves: Sylvidae). Journal of Zoology. 1979; 189: 233–258
  56. 56. Zann RA, Morton SR, Jones KR, Burley NT. The timing of breeding by Zebra Finches in relation to rainfall in Central Australia. Emu. 1996; 95: 208–222
  57. 57. Grant PR, Grant BR, Keller LF, Petren K. Effects of El Nino events on Darwin’s Finch productivity. Ecology. 2000; 81: 2442–2457
  58. 58. Scott DM, Lemon RE, Darley JA. Relaying interval after nest failure in Gray catbirds and Northern Cardinals. The Wilson Bulleti. 1987; 99: 708–712
  59. 59. Minot EO. Effects of interspecific competition for food in the breeding Blue and Great tits. Journal of Animal Ecology. 1981; 50: 375–385
  60. 60. Sansom A, Lind J, Cresswell W. Individual behavior and survival: the role of predator avoidance, foraging success and vigilance. Behavioral Ecology. 2009; 20: 1168–1174
  61. 61. Fontaine JJ, Martin TE. Parent birds assess nest predation risk and adjust their reproductive strategies. Ecology Letters. 2006; 9: 428–434 pmid:16623728
  62. 62. Langen TA. Prolonged offspring dependence and cooperative breeding in birds. Behavioral Ecology. 2000; 11: 367–377