Due to climate change, extreme weather events are becoming more frequent and severe. These extreme events have been documented to affect avian predators in stream ecosystems. To better understand the mechanisms behind this effect, we used a decade-long dataset from a mountain stream in Taiwan to assess the effects of extreme flooding caused by typhoons on invertebrate abundance in different periods of the year and the resulting effects on reproductive output of an avian predator of these invertebrates, the Brown Dipper (Cinclus pallasii). In this study stream, all extreme floods occurred between June and October, and these floods negatively affected invertebrate density. Consequently, average invertebrate density was lowest in October at the end of the typhoon season, and highest 4 mo later. Because invertebrate density increases over time after a flood, the length of the recovery period between floods that occurred between June and October was more important than each flood's magnitude in determining invertebrate density in October. October invertebrate density then positively correlated with invertebrate density, the number of dipper breeding pairs, and the proportion of breeding individuals in the following breeding season, which lasted from January to early April. The effects of lower invertebrate densities in October and then February percolated all the way through the system, affecting laying date, fledgling population, and the next winter's population. Given our results, an increase in the frequency of typhoons, especially late-season typhoons, will have a negative effect on Brown Dipper reproductive output through bottom-up effects in stream ecosystems.
Climate change is altering the phenology, distribution, and demographic rates of many species (e.g., Parmesan et al. 2000, Hussell and Brittingham 2003, Lof et al. 2012, Dybala et al. 2013, Jenouvrier 2013). One of the consequences of climate change is that, in addition to an overall warming trend, extreme weather events have become more frequent and severe, and these trends will increase (IPCC 2013). Weather-related disturbances, such as heat waves, droughts, storms, and floods, can change or limit a species' reproduction, distribution, and abundance (Parmesan et al. 2000, McCarty 2001, Elkins 2010). Extreme weather events can thus directly affect the survival and reproduction of individuals, or indirectly affect them by limiting the availability of habitat or other resources (Jenouvrier 2013). Since extreme weather events are rare and difficult to predict, their effects on species' survival and reproduction are still poorly understood (Reed et al. 2003, Jenouvrier 2013).
The 5 species of dipper (Cinclus) occur on 5 continents, where they inhabit fast-flowing rivers in which they feed mainly on aquatic macroinvertebrates (Tyler and Ormerod 1994, Ormerod and Tyler 2005). Dippers are sensitive to acidification and pollution, and are therefore considered bioindicators of stream habitats (Ormerod et al. 1985, Tyler and Ormerod 1992, Ormerod and Tyler 1993, Sorace et al. 2002, Morrissey et al. 2004). They are also sensitive to the average air and water temperature, which influences population dynamics and laying dates of White-throated Dippers (Cinclus cinclus) (Sæther et al. 2000, Hegelbach 2001, D'Amico et al. 2003). Additionally, dippers are sensitive to extreme weather events: floods can cause their reproduction to fail (Shaw 1978, Price and Bock 1983, Morrissey 2004, Royan et al. 2014), change their foraging behavior and diet (Da Prato 1981, Taylor and O'Halloran 2001), and decrease their survival (Lebreton et al. 1992, Marzolin 2002, Chiu et al. 2013), abundance (Chiu et al. 2008), and distribution (Royan et al. 2015). Changes in flow regime are also associated with changes in the presence and absence of White-throated Dippers (Royan et al. 2013).
A river's flow regime is also one of the key factors affecting macroinvertebrate composition and abundance in stream ecosystems. Based on previous work, we know that extreme floods triggered by typhoons in Taiwan result in a large reduction in macroinvertebrate abundance, which, in turn, causes the numbers (Chiu et al. 2008, Chiu and Kuo 2012) as well as the survival rate (Chiu et al. 2013) of one of their avian predators, the Brown Dipper (Cinclus pallasii), to decrease. After such an extreme flooding event, the total abundance of invertebrates generally recovers again over a matter of months (Greenwood and Booker 2014), and the length of this recovery period affects macroinvertebrate abundance and composition (Chiu and Kuo 2012). Thus, the time period between a flood and the Brown Dipper's breeding season may be a crucial factor that determines prey availability, and therefore breeding onset and output.
In this study, we used long-term data on ambient temperatures and flow regimes of streams in Taiwan to investigate the effect of extreme weather events precipitated by typhoons on invertebrate abundance and the abundance and reproduction of the Brown Dipper. In Taiwan, typhoons usually cause extreme floods from June to October (Tu et al. 2009), while the breeding season of the Brown Dipper lasts from January to early April (Hong et al. 2011). We therefore hypothesized that extreme floods would affect the dipper's breeding phenology and reproductive output via bottom-up impacts caused by decreased prey abundance prior to the breeding season, or even lasting into the breeding season. We used information-theoretic model selection methods (Burnham and Anderson 2002) to determine whether ambient temperature, the magnitude of the flood, or the length of the recovery period after the flood influenced invertebrate abundance and, in turn, the abundance and reproductive output of dippers.
Study Area and Environmental Variables
Our study was conducted from 2003 to 2014 in the upstream drainage of the Tachia River in central Taiwan, which includes the Cijiawan (also spelled Chichiawan), the Gaoshan, and the Yousheng streams (Figure 1). The watershed covers 77 km2, and elevation ranges from 1,700 to 2,000 m. Air temperature data is recorded at Siyuan weather station (24.3900°N, 121.3482°E), which is located near the Yousheng Stream, ∼5 km from our study area. During the study period, daily mean air temperature ranged from 4.9 ± 4.2°C in January to 16.7 ± 3.1°C in July (mean ± SD). To calculate the winter temperature prior to the Brown Dipper's breeding season, we averaged daily air temperature from December 1 to January 31.
Daily mean water flow, measured in cubic meters per second (m3 s−1), was recorded by the Taiwan Power Company at a site downstream of the confluence of the Cijiawan and Yousheng streams (Figure 1). From all the daily means, we took the maximum value for each of 2 periods, namely June 1–October 31 and November 1–January 31, and used these maximum values as estimates of maximum summer flow and maximum winter flow, respectively. In other words, we used the water flow from the day with the largest water discharge as a measure of the flood with the largest magnitude, because we assumed a priori that the largest flood would have the strongest impact on the invertebrate fauna. The summer period was chosen because the typical typhoon season in Taiwan usually lasts from June to October (Tu et al. 2009), and the winter period was selected because it precedes the breeding season.
Brown Dipper Population, Breeding, and Food Surveys
Dippers are habitat specialists that occupy a linear territory. In our case, this linear habitat was always a rocky and narrow (5–20 m wide) riverbed with very little vegetation that was always bounded by steep, rocky, vegetated cliffs (Figure 2). The probability of observers visually detecting dippers by walking along a stream is very high. For example, D'Amico and Hemery (2003) showed that 3 visits (i.e. surveys) were sufficient to detect the entire breeding population of White-throated Dippers in mountainous stream habitats.
We conducted monthly population surveys from November to January to census the potential breeding population for the following breeding season. During each survey, we counted individual Brown Dippers by slowly walking along the stream's edge through the entire defined study area (8.5 km), which was one sampling unit. Our study area extended 7 km up the Cijiawan Stream from its confluence with the Yousheng Stream and 1.5 km up the Gaoshan stream from its confluence with the Cijiawan stream (Figure 1).
Individual dippers almost invariably move away from the observer before they “double back” and fly past the observer (and never through the high vegetation adjacent to the riverbed). Therefore, we ignored individuals that flew ahead of the observer to avoid re-counting birds. Each survey lasted from 08:00 to 16:00, with the lower part of the study area covered on the first day, and the upper part on the second day. Further details of and justifications for our survey methods, including a detailed description of the dipper's territorial behavior, are given in Chiu et al. (2008) and Chen and Wang (2010).
We are confident that our surveys resulted in a complete (or almost complete) census of the dipper population present at the time. Furthermore, we calculated the repeatability (Krebs 1999) of the 3 monthly censuses, which yielded a relatively high repeatability of 0.859 (95% CI: 0.814–0.926, n = 27 mo). Most of the remaining variation was due to some floater individuals that moved rather than remaining on territories during the prebreeding season. We wanted to obtain the potential breeding population during the prebreeding season, so we included floaters. Therefore, we used the maximum (instead of the mean) number of dippers counted during the 3 monthly surveys as our estimate of the winter population.
During the breeding season, dippers typically nest on cliff ledges, in natural or artificial cavities, and under bridges (Severinghaus et al. 2010, Hong et al. 2011). We conducted at least 2 breeding surveys per month for dipper nests from January to March by slowly walking along the streams and following birds with binoculars when they flew to their nests. Dippers are territorial, with a typical distance between adjacent nests of 500–1,000 m in Taiwan, with variation explained by food abundance (Chen and Wang 2010). Because of our prior winter population censuses (see above) and our constant observations during the breeding season, it was highly unlikely that we missed any nests of pairs that continuously occupied a territory. We thus determined the number of breeding pairs as the total number of nests into which a female placed eggs. In some years, a few pairs built a nest but then abandoned it before laying eggs; these “attempted nests” were not included in our number of breeding pairs. In a few cases, we observed a second brood within the same breeding season, but we never observed a pair building a second nest in another location. Since most breeding adults (and some of the floaters) were color-banded, we were able to ascertain if a pair moved to another location. We then used the number of breeding pairs to calculate the breeding percentage as (2 × breeding pairs) / winter population.
Once a nest was found, we observed the nest every 3–5 days to determine the breeding status. We determined the laying date by direct observation or by backdating from the hatching date or based on nestling weight (Hong et al. 2011). For each nest, we counted the number of days from January 1 until a nest contained the first egg. We then averaged the laying dates of the first 2 nests for our variable of first laying date, because we believed that an average from 2 nests would be a more reliable estimate of the onset of laying for the entire population.
For accessible nests, we counted the number of almost-fledged nestlings when we banded them at 18–20 days of age. We carefully observed inaccessible nests to count the fledglings when they left the nest. Over our study years, the percentage of such inaccessible nests ranged from 0% to ∼12%, averaging ∼6%. We thus determined the fledgling population as the total number of fledglings produced by all the nests.
To examine the main food source for dippers, we sampled benthic macroinvertebrates using a Surber sampler (area = 30.48 cm × 30.48 cm, mesh size = 250 μm) at 4 sites (Figure 1) every 2 mo from February 2003 to February 2014. All intervals were sampled in all years, except December, which was only sampled in 2003, 2004, and 2006. In each sampling period, we took 6 samples from each site at randomly selected locations within runs and riffles. Samples were preserved in 75% ethanol. We identified most individuals to genus or species using published keys (Kang 1993, Merritt and Cummins 1996, Kawai and Tanida 2005). According to Chiu et al. (2009), dippers in Taiwan feed mainly on large invertebrates, primarily Diptera, Ephemeroptera, Plecoptera, and Trichoptera. To estimate food abundance per square meter (m−2), we multiplied the mean number of invertebrates of these 4 taxa caught in the 24 samples (4 sites × 6 samples) by a factor of 10.764 (i.e. 10,000 cm2 / [30.48 cm × 30.48 cm]) in order to derive an estimate of invertebrate density. In our model selection analysis, we then used only June, August, October, and February invertebrate density (Table 1).
Environmental variables and biotic response variables used in our model selection analysis of the influence of extreme floods on Brown Dippers and their macroinvertebrate prey. Sample size (n) is the number of years for which a variable was measured; the sampling period was 2003–2014, thus n =12 unless data are missing.
Because typhoon-triggered floods are known to affect the dipper's invertebrate prey, we also determined the recovery period as the number of days between the date of the most recent flood during the typhoon season from June to October and the date when invertebrates were sampled in October for the first time at the end of the typhoon season. A flood was defined as a maximum water flow of >100 m3 s−1, which was a convenient threshold because the 3 yr without a typhoon, and thus without a flood, were 2003, 2010, and 2011, with maximum flows of 15, 58, and 41 m3 s−1, respectively, while all other years had at least 1 typhoon, with the smallest resulting maximum flow being 141 m3 s−1 in 2005. The largest maximum flow of 610 m3 s−1 occurred in 2004. For years without a flood, we used the beginning of the typhoon season (June 1) as the starting date for the recovery period.
We tested for various reasonable relationships among the 4 environmental variables and 8 biotic responses (Table 1) using multiple linear regression. To achieve a normal distribution for each variable, all independent variables and response variables (i.e. all variables in Table 1) were log-transformed prior to all analyses (see Chiu et al. 2008), except for August invertebrate density and fledgling population, which required a Box-Cox transformation (Krebs 1999) to achieve a normal distribution (one-sample Kolmogorov-Smirnov test or Shapiro-Wilk P > 0.05).
For each biotic response variable (i.e. the dependent variable) listed in Table 2 by capital letters A–H, we first used all of the other 12 variables listed in Table 1 as independent variables (the correlation matrix for all variables is given in Appendix Table 4). However, in each case, only a few variables made an important contribution according to AICc value and Akaike weight (wi; see below). Therefore, in the second step, we selected the 3 independent variables from among the 12 variables presented in Table 1 which (1) were the most important according to ΔAICc, but also (2) were most likely to influence the response variable using prior knowledge (see Introduction and Methods). We thus followed Burnham and Anderson's (2002) suggestion to only test for relationships which make biological sense and restricted our analysis to testing only the most ecologically reasonable models. For each response variable (Table 2), we then tested the 8 candidate models resulting from all of the possible combinations of the 3 independent variables plus the null (intercept only) model.
Ranking of all possible models for 8 biotic response variables (see definitions in Table 1) of Brown Dippers and their macroinvertebrate prey to factors influenced by extreme flood events in Taiwan, 2003–2014. Models are ranked by the difference from the top model in Akaike's Information Criterion corrected for small ample size (ΔAICc) and Akaike weight (wi). The parameter count (K) and model deviance (−2logL) is also shown. Models with substantial support (ΔAICc < 2) are in bold font. Response variables with positive relationships are in roman font; italic font indicates negative relationships.
The AICc values of the top models are as follows: A. Aug. invert. = −17.6; B. Oct. invert. = 4.3; C. Feb. invert. = 4.9; D. Winter pop. = −14.1; E. Breed. Pair = −20.0; F. Breed. perc. = −19.6; G. First laying = −22.6; H. Fledg. pop. = 72.1.
For this analysis, summer flow and recovery period were defined for the period from June 1 to the invertebrate sampling date in August.
The winter population of the previous year.
We ranked models based on Akaike's Information Criterion adjusted for small sample size (AICc). We accepted only those models with a ΔAICc < 2 as having “substantial” support (Burnham and Anderson 2002, Grabowska-Zhang et al. 2011, Majić et al. 2011). If there was model selection uncertainty (i.e. >1 model with ΔAICc < 2), we averaged the accepted high-ranking models to obtain parameter estimates. To assess the relative importance of each predictor, we also calculated the relative variable importance (RVI) by summing the Akaike weights (wi) for all of the accepted models with each selected variable (Burnham and Anderson 2002), thereby ranking variables by importance according to their contribution to the entire set of accepted models. We conducted all model ranking and model averaging procedures with the MuMIn package (Bartoń 2012) applied in program R 3.1.0 (R Core Team 2012).
During the study period, there were 15 flooding events with a water flow of >100 m3 s−1, each lasting from 1 to 3 days (Figure 3A, Appendix Figure 6). These floods all occurred between June and October, while the water flow during the remaining months was comparatively low and stable. The higher magnitude of water flow was mirrored by a decrease in invertebrate density; the period from June to October had the lowest mean invertebrate densities (Figure 3B).
Each of the 8 response variables was effectively predicted by models with only 1 or 2 predictors (Table 2). For 5 of the 8 response variables, only 1 independent variable had substantial model support as a correlated variable, while the breeding percentage was correlated with 2 independent variables. For the 2 remaining response variables, August invertebrate density and breeding pairs, 2 candidate models had substantial support. Table 3 presents the results from averaging these 2 models. Using the relative variable importance as a measure, summer flow was more strongly correlated with August invertebrate density than was June invertebrate density, and October invertebrate density had a higher correlation with breeding pairs than winter population. All of the correlations with substantial support in Table 2 are depicted graphically in Figure 4.
Results of model-averaging the top models (ΔAICc < 2) with >1 variable selected to determine the influence of extreme floods on Brown Dippers and their macroinvertebrate prey in Taiwan, 2003–2014 (see Tables 1 and 2 for corresponding capital letter labels and variable definitions). Coeff. = coefficient of explanatory variable; SE = unconditional standard error; CL = 95% confidence limit; RVI = relative variable importance.
For this analysis only, summer flow was defined for the period from June 1 to the invertebrate sampling date in August.
To summarize the information contained in Tables 2 and 3 and Figure 4, we constructed a flow diagram depicting the directions of the relationships (Figure 5). Higher summer flow and lower June invertebrate density were associated with lower August invertebrate density (Table 2, Figure 4A), but summer flow was of greater importance (Table 3). Only the length of the recovery period positively influenced October invertebrate density (Table 2, Figure 4B). October invertebrate density positively influenced February invertebrate density (Table 2, Figure 4C), and higher February invertebrate density was related to an earlier laying date (Table 2, Figure 4G), which then led to a higher fledgling population (Table 2, Figure 4H). Therefore, the influence of the recovery period filtered all the way down to the fledgling population in the subsequent spring breeding season. The fledgling population then had a positive effect on the population in the following winter (Table 2, Figure 4D), which in turn had a positive influence on the number of breeding pairs in the subsequent breeding season (Table 2, Figure 4E-2). However, the number of breeding pairs was influenced more by the October invertebrate density before the breeding season (Tables 2 and 3, Figure 4E-1) than by the winter population just before the breeding season (Tables 2 and 3, Figure 4E-2). Finally, more individuals of the entire population bred (i.e. the breeding percentage was higher) if the October invertebrate density was higher (Table 2, Figure 4F-1), but fewer individuals of the entire population bred (i.e. the breeding percentage was lower) if the winter population was higher (Table 2, Figure 4F-2).
Our long-term study monitored the population and breeding parameters of the Brown Dipper, as well as its invertebrate prey and some weather-related variables which might have influenced this study system, for more than a decade. The strengths of this particular study system are that (1) population counts of dippers are highly accurate because of the visual openness and constrained linear shape of the habitat; (2) dippers feed mainly on aquatic macroinvertebrates, which can be reliably sampled; and (3) this study and other studies have demonstrated the strong influence of water flow on this macroinvertebrate community. Therefore, the correlations summarized in Figure 5 are most likely causal relationships, although correlative evidence is naturally always open to alternative explanations.
Previous studies have shown negative correlations between water flow, invertebrate density, and dipper populations (Chiu et al. 2008), and the negative influence of extreme flooding on dipper survival (Chiu et al. 2013). In addition to using a temporally longer dataset, our study also provides a comprehensive and more mechanistic understanding of how extreme flooding affects the reproductive output and subsequent abundance of dippers.
First, if the initial population of macroinvertebrates sampled in June was higher, the population sampled in August was also higher. Second, logically, invertebrate density in August was negatively influenced by summer flow, especially if the summer flow was a major flood. This influence of water flow (and especially extreme water flow) on macroinvertebrate communities has also been demonstrated in previous studies (Chiu et al. 2008, Chiu and Kuo 2012). Furthermore, Chiu and Kuo (2012) showed that, after a typhoon-triggered flood, invertebrate density increased over time. Consequently, the length of the recovery period was the only important variable in our study that determined invertebrate density in October, at the end of the typhoon season, while a flood's magnitude had no influence on invertebrate density in October. This result meant that October invertebrate abundance was lower when extreme flooding was caused by a typhoon later in the typhoon season than if the flooding occurred earlier in the typhoon season. We also demonstrated that invertebrate abundance in October influenced all of the subsequent breeding parameters and outputs of dippers (Figure 5). Therefore, flood timing was more influential than flood magnitude on Brown Dipper breeding.
Interestingly, a higher proportion of individuals bred if October invertebrate density was higher; presumably, the higher food supply improved individual condition (Wendeln and Becker 1999, Cubaynes et al. 2011, Lehikoinen et al. 2011, Reichert et al. 2012, Millon et al. 2014). Alternatively, the higher food supply may have decreased territory size, thus allowing more individuals to breed (Chen and Wang 2010). However, a lower proportion of individuals bred if the winter population was higher, suggesting that the total number of territories was limiting. Dippers also could have been limited by the number of suitable nest sites (Loegering and Anthony 2006, Severinghaus et al. 2010). Dippers in our study usually began building their nests in December, but started as early as mid-November in a few years. Thus, competition for territories and nest sites began long before breeding. Similar observations were made by Price and Bock (1983), who reported that the poorest average reproduction occurred in the year with the largest study population of American Dippers (Cinclus mexicanus).
In our study population, a sufficiently high number of successful fledglings managed to survive during the summer and autumn to supplement the next winter's population and the next spring's number of breeding pairs (Figure 5). Similar results have been documented for other species (e.g., Sillett et al. 2000, Sæther et al. 2002). Demographic responses to interannual variations in weather may be a typical pattern revealed with large datasets (e.g., Dugger et al. 2016).
A species' laying period usually coincides with favorable weather conditions or seasonal peaks in food abundance. In many regions, these favorable periods are shifting with higher spring temperatures, and species are adjusting their phenology accordingly (e.g., Van Noordwijk et al. 1995, Parmesan and Yohe 2003). Brown Dippers begin to breed in January, which is the coldest month in our study area, but is also when stream invertebrates have recovered from floods caused by typhoons, and this timing further allows fledglings sufficient time to mature before the next typhoon season. Unlike the findings of several other studies on dippers (Sæther et al. 2000, Hegelbach 2001, D'Amico et al. 2003), winter temperature had no influence on invertebrate density or laying date in our study (Table 2). Perhaps the effects of temperature were masked by the significant effects of flooding, although we found no statistical effect of temperature. The only other variable that we examined which had no discernable influence on our study system was maximum winter flow (Table 2). The most likely explanation for this result is that the period of winter flow does not include the typhoon season and therefore has a much lower and more stable flow regime, with little influence on invertebrate density (Figure 3).
Many studies of dippers and other bird species have demonstrated the positive influence of higher food abundance on laying dates and fledgling populations, and of earlier laying dates on clutch sizes and fledgling populations (e.g., Price and Bock 1983, Van Noordwijk et al. 1995, Morrissey 2004, Visser et al. 2006, Lehikoinen et al. 2011; but see Wilson 1996). Surprisingly, the fledgling population was correlated only with the first laying date and not with the number of breeding pairs within the same breeding season (Table 2H). Since other studies have found a positive relationship between the number of breeding pairs and the total number of fledglings that they produce (e.g., Holland and Yalden 2002, Millon et al. 2014), our result is somewhat surprising. The positive correlation between the number of breeding pairs and the number of fledglings produced explained only 34% of the variation in our model. Thus, there was some effect, but the effect was small compared with the effect of first laying date, especially given the fact that adding the number of breeding pairs into the model with first laying date added only another 6% of explained variation.
In some years, stochastic factors, such as bad weather or predators, cause breeding failures (Hong et al. 2011, Y.-H. Sun personal observation) and thus diminish the influence of the number of breeding pairs on the production of fledglings. On the other hand, during years when breeding started early (usually years with no floods before the breeding season), a few pairs produced 2 clutches, thus increasing the total number of fledglings. These double broods may explain why the fledgling population was positively correlated with first laying date, but not with the number of breeding pairs.
As always, correlational evidence is open to alternative explanations. For example, February invertebrate density could conceivably be correlated with a third variable that also influences first laying date (cf. Table 2); perhaps higher temperatures in February cause both higher invertebrate densities and earlier laying dates. In lieu of another explanation, we conclude that bottom-up effects percolate all the way through our study system, from floods to invertebrate density to dipper abundance. A similar bottom-up effect of benthic prey distribution has also been implicated in the distributional shifts of an Arctic duck (Sexson et al. 2016).
Given our results, it is highly likely that an increase in the frequency and severity of extreme climatic events, such as typhoons, will have a negative effect on both the stream invertebrate communities and their main avian predator, the Brown Dipper. These conclusions, naturally, assume no evolutionary adaptation of the invertebrates or the Brown Dipper to these future environmental changes (Visser 2008, Hendry et al. 2010).
The Brown Dipper is the only dipper species which lives in a region that is frequently affected by typhoons. Because typhoons hit Taiwan annually, severe river floods are common (Chen and Chen 2003), and the number and severity of these typhoons is predicted to increase due to climate change (Allan and Soden 2008, Elsner et al. 2008). Typhoons in the western North Pacific are becoming more intense and longer-lasting (Webster et al. 2005, Chiang and Chang 2011). In Taiwan, records show that rainfall from typhoons of all intensities has increased significantly, while other rainfall, particularly lighter rain, has tended to decrease (Tsai and Huang 2011, Tu and Chou 2013, Chu et al. 2014). The intense rainfall from typhoons usually causes extreme flooding of streams, and these floods have increased in magnitude (Appendix Figure 6; Chang et al. 2011, Chiang and Chang 2011, Chiu and Kuo 2012). Given our results, we can thus make a tentative prediction that climate change will negatively impact the Brown Dipper's reproductive output in Taiwan due to increased frequency and severity of floods. Interestingly, if an increase in typhoon frequency were also to be associated with a shift toward earlier or later typhoon occurrence, then the effects on the Brown Dipper would correspondingly change. Naturally, other effects of climate change, such as higher summer temperatures, heat waves, or droughts, may also affect the dipper's food and breeding output, but through different mechanisms.
Streams and rivers not only satisfy human need for freshwater, but also harbor one of the earth's most biodiverse ecosystems. Climate change already affects them, and will increasingly do so (Cumming et al. 2012, Royan et al. 2013). Our results add to the mounting evidence that, at least in the short term, many biological communities will be negatively affected by the impacts of climate change and, in particular, extreme weather events resulting from climate change.
We thank C.-C. Chen, T.-S. Ding, R.-S. Lin, and 2 anonymous reviewers who gave useful comments on the manuscript.
Funding statement: Our research was supported by grants from the Shei-Pa National Park, Taiwan, and partly supported by the National Science Council of Taiwan. B.A.W. acknowledges financial support from Taipei Medical University. None of the funders had any input into the content of the manuscript. None of the funders required their approval of the manuscript before submission or publication.
Ethics statement: This project was conducted under permits from Shei-Pa National Park and the Wildlife Conservation Act of Taiwan.
Author contributions: Y.-H.S. conceived the idea, design, and experiment (supervised research, formulated question or hypothesis); S.-Y.H. and M.-C.C. performed the experiments (collected data, conducted the research); S.-Y.H. and B.A.W. wrote the paper (or substantially edited the paper); Y.-H.S., M.-C.C., and M.-H.K. developed or designed the methods; and S.-Y.H. and B.A.W. analyzed the data.
APPENDIX TABLE 4
Pearson correlation coefficient matrix for the 13 variables used in our analyses of the effects of extreme flood events on Brown Dippers and their macroinvertebrate prey in Taiwan, 2003–2014 (see Table 1 for variable definitions).