Extreme rainfall events alter the trophic structure in bromeliad tanks across the neotropics

ABSTRACT Changes in global and regional precipitation regimes are among the most pervasive components of climate change. Intensification of rainfall cycles, ranging from frequent downpours

to severe droughts, could cause widespread, but largely unknown, alterations to trophic structure and ecosystem function. We conducted multi-site coordinated experiments to show how

variation in the quantity and evenness of rainfall modulates trophic structure in 210 natural freshwater microcosms (tank bromeliads) across Central and South America (18°N to 29°S). The

biomass of smaller organisms (detritivores) was higher under more stable hydrological conditions. Conversely, the biomass of predators was highest when rainfall was uneven, resulting in

top-heavy biomass pyramids. These results illustrate how extremes of precipitation, resulting in localized droughts or flooding, can erode the base of freshwater food webs, with negative

implications for the stability of trophic dynamics. SIMILAR CONTENT BEING VIEWED BY OTHERS ASYNCHRONOUS RECOVERY OF PREDATORS AND PREY CONDITIONS RESILIENCE TO DROUGHT IN A NEOTROPICAL

ECOSYSTEM Article Open access 19 May 2022 WARMING INDIRECTLY SIMPLIFIES FOOD WEBS THROUGH EFFECTS ON APEX PREDATORS Article Open access 05 October 2023 RECOVERY OF FRESHWATER MICROBIAL

COMMUNITIES AFTER EXTREME RAIN EVENTS IS MEDIATED BY CYCLIC SUCCESSION Article 28 January 2021 INTRODUCTION Climate change is predicted to dramatically alter precipitation regimes and global

hydrological cycles1,2,3. Although changes in the spatial distribution of rainfall can both mitigate and amplify differences between dry and wet regions, there is a consensus that many

regions will suffer severe impacts of increased variability and magnitude of precipitation1,2,3,4. These climatic fluctuations can cause extreme hydrological events, such as flooding and

drought, which can lead to widespread, though largely unknown, shifts in ecosystem structure and function, particularly in freshwater ecosystems5,6,7,8. Such events can expose some

ecosystems to conditions with no recent historical precedent6. Whereas experimental research has focused on incremental changes in mean conditions, fluctuations or extreme events such as

floods and droughts may have more profound ecosystem consequences6. In addition, most studies simulating precipitation extremes at the multi-site scale have been limited to a single trophic

level, particularly producers5. Therefore, there is a clear and urgent need to identify which food web compartments (decomposers, primary consumers, predators) are most vulnerable to

rainfall fluctuations and extremes5,6. Whereas multitrophic research has focused on a single direction of extreme change, especially drought7,8,9,10,11,12,13,14,15,16, the ecological

consequences of drought and flooding have rarely been explored in concert17. Such studies suggest that drought can substantially alter aquatic7,8,9,10,11,12,13,14,15 and terrestrial food

webs16,17,18,19, with consequences for community structure and ecosystem function. Drought also weakens trophic cascades and the strength of biological interactions (e.g., competition,

predation), and disproportionately threatens top predators, often resulting in communities dominated by smaller organisms6,7,8,9,10,11,12,13. However, previous work has provided limited

mechanistic understanding of differences in ecosystem sensitivity to global change5, has been conducted at a local scale6, and has used dissimilar experimental approaches and methods5. This

hampers our ability to predict global impacts of drought and flooding on multiple taxa and trophic levels, including standing stock biomass, trophic biomass pyramids20,21,22, and biomass

fluxes through the food web23, across large geographic regions. Standing stock biomass is a common metric in food web research, with trophic structure represented by the distribution of

biomass across different trophic levels21,22. These biomass pyramids can summarize changes in complexity and biomass flux through food webs20,21,22. Pyramid size and shape exhibit highly

variable patterns across different types of ecosystems worldwide20, but it has been shown that climatic stability can change the shape of biomass pyramids24. Thus, quantifying biomass

pyramids improves mechanistic understanding of climate change effects on food web structure, resource partitioning, and energy use. Here we conducted a geographically coordinated experiment5

in seven sites across Central and South America (18°N to 29°S, Fig. 1a) to investigate the effects of rainfall changes on trophic structure24,25. We used natural, detritus-based microcosms

(bromeliad phytotelmata) as model systems due to their widespread distribution and ease of manipulation10,11,12. Bromeliad aquatic ecosystems are inhabited by a diverse fauna10,11,12,

comprising top predators, mesopredators, and detritivores15,24,25. We contrasted rainfall-mediated changes in hydrological stability of the study system with the effects of two main rainfall

components: (i) the mean daily amount of rainfall, _μ_; and (ii) distribution of rainfall events around this mean through time, _k_ (i.e., a measure of evenness in the frequency

distribution of rainfall; hereafter “frequency”; Fig. 1b). For instance, reductions in both mean daily rainfall (low _μ_) and the even frequency of rain (low _k_, hereafter “infrequent

rainfall”) increase the proportion of days that bromeliads are empty of water. Regular (current) variability of the rainfall components in each site were first determined using recent

meteorological data (see Methods). We applied a negative binomial distribution to these data to estimate the parameters μ and k. We applied ten levels of _μ_ (ranging from 0.1 to 3.0) and

three levels of k ranging from 0.5 to 2.0 in a fully factorial experimental design at each of our seven sites for a total of 10 × 3 × 7 = 210 food webs in individual bromeliads. This allowed

us to compare ambient, baseline conditions (_μ_ = 1, _k_ = 1) and extreme fluctuations of rainfall quantity (10–300%) and frequency (50–200%) to average historical levels of daily

variability for each site (“Methods”). We hypothesized that precipitation would affect trophic structure in bromeliads by altering hydrology, that is, the temporal dynamics of water within

the bromeliad. We therefore also quantified the underlying hydrological dynamics within bromeliads in the field (see “Methods”) as a potential proximate driver, and then projected it in

multivariate space using principal component analysis, after standardization between sites (Supplementary Table 2, Supplementary Fig. 1). We defined top predators as species without natural

predators within the aquatic food web. All bromeliads were open to colonization and extinction throughout the two-month experimental duration and thus communities could dynamically assemble

or disassemble. The multi-site approach allowed us to explore generalities and site contingencies in food web responses to climate change. We expected stronger impacts of rainfall events

when their components departed from current scenarios (_μ_ ≠ 1, _k_ ≠ 1). If bigger predators are more sensitive to drought (here measured as lower values of _μ_ and/or _k_) than smaller

organisms (e.g., mesopredators and detritivores)6,7,8,24, then drought could have stronger ecological impacts than heavy rainfall6. Consequently, under drought or extremely infrequent

rainfall events we expected communities dominated by smaller organisms7 (especially detritivores and filter feeders), which tend to be suppressed via top-down control under more favorable

conditions (Fig. 1c). Conversely, an excessive amount of rainfall combined with higher frequency of rainfall among consecutive days could impact lower trophic levels through hydrodynamic

perturbations, e.g., if important nutrients and basal resources (detrital organic matter and microorganisms) are lost to flooding26,27. Bigger predators could be more resistant to flooding,

but their biomass could still decrease if they have fewer resources to support them. Consequently, such changes in standing stock biomass could alter predator–prey biomass ratios (i.e.,

pyramid shape; Fig. 1d)20,21,24. Because drought and high rainfall frequency are predicted to favor lower and higher trophic levels, respectively, drought is expected to be associated with

bottom-heavy biomass pyramids, and high rainfall frequency is expected to be associated with top-heavy biomass pyramids (Fig. 1d). Both the climate and the regional species pool are

different among sites12,25,27, which may lead to site-specific contingencies, e.g., communities from arid regions and/or regions with large natural amplitudes in rainfall may be more

resistant to drought than those from wet regions and/or regions with narrow climatic amplitudes27. Extreme precipitation events and underlying hydrological dynamics influence each trophic

level in a different manner. Whereas the biomass of detritivores, often the smaller organisms in the study system, is higher under more stable hydrological conditions, the biomass of

predators is highest under uneven rainfall (drought conditions). Higher resource concentration under drought conditions fuels these higher trophic levels, resulting in top-heavy biomass

pyramids. Our results demonstrate that organisms from lower trophic levels may be the most susceptible to ongoing climate change. RESULTS STANDING STOCK Top predators, mesopredators, and

detritivores responded differently to rainfall components, and only detritivores were influenced by hydrological stability (Table 1). Predator standing stock biomass decreased with

increasing rainfall frequency across all sites in a remarkably consistent response (i.e., no site vs. rainfall interactions, Fig. 2a), despite large differences between sites in average

standing stock (Table 1) and taxonomic composition27. These emergent patterns indicate that basic properties and processes are recurring in different food webs and deserve deeper

understanding. First, top predators account for the majority of predator biomass due to their large body sizes. Second, contrary to initial predictions, total predator biomass (top predators

and mesopredators) was not negatively affected by infrequent rainfall. Instead, the standing stock of this trophic group increased as rainfall distribution became more infrequent, and

decreased under more even rainfall dispersion (frequent), relative to ambient conditions (_P_ = 0.045, Fig. 2a, Table 1). A similar pattern was observed if only top predators were evaluated,

though with marginally non-significant results (top predators; _P_ = 0.067, Table 1, Supplementary Fig. 2). In contrast to top predators, detritivore standing stock decreased under extremes

of rainfall frequency at some sites and increased in others, thus indicating strong site-specific responses (Table 1, Fig. 2b). In addition, detritivores were more sensitive to hydrological

instability, with biomass being higher under more stable hydrological conditions (Fig. 2c). In contrast, mesopredators were not affected by rainfall or hydrological stability (Table 1).

BIOMASS PYRAMIDS Differential responses among the trophic levels resulted in consistent rainfall-driven shifts in the shape of biomass pyramids across sites (i.e., no site vs. rainfall

interactions for biomass ratios; Fig. 2d, Table 1). Biomass pyramids comprised of all predators (meso and top predators) consistently became more top-heavy in many communities (i.e.,

increased total predator-detritivore mass ratio [PDMR]), and even inverted (PDMR > 1), under infrequent (low k) rainfall conditions. A similar pattern emerged for top predator-detritivore

ratios (Table 1, Supplementary Fig. 2). This pattern was driven primarily (i) by top predators, rather than by mesopredators (Table 1), and also (ii) by predator body size (Fig. 3) rather

than predator abundance. Predator abundance declined with increasing amount of rainfall (_μ_; Supplementary Fig. 3), but was not influenced by rainfall frequency (_P_ < 0.05, backward

selection). We also found a strong non-linear increase of PDMR with increase in whole system biomass, but the curves plateaued with accumulating biomass (Fig. 4a). This pattern was

determined exclusively by predator biomass (Fig. 4b), whereas an increase in detritivore biomass decreased PDMR (Fig. 4c). DISCUSSION We explored the simultaneous impacts of drought and

flooding extremes on lentic ecosystems. Our results contrast with findings from lotic ecosystems8,9,10,11, which found higher trophic levels to be more susceptible to drought or hydrological

disturbances than lower trophic levels9,15. Like many ponds and other wetlands, bromeliads are lentic water bodies that are naturally prone to partial desiccation and overflow. Larger

predators from these environments might be better adapted to drought, bouncing back quickly from perturbations. For such lentic ecosystems, flooding events seem to be at least as harmful as

drought, as flooding can leach out important nutrients, basal resources (organic matter, microorganisms)26,27 and even small macroinvertebrates (e.g., detritivores)27 from the ecosystem.

Indeed, we found turbidity, a measure that integrates organic and inorganic suspended matter in ecosystems, including free-living algae and particulate nutrients (C, N, P), resulting from

detritivore activity, to be lower under more frequent rainfall conditions in our experiments (Supplementary Fig. 4). This can make ecosystems less productive and thus less efficient at

sustaining higher trophic levels. Predators may still be able to persist in such drought conditions because of high biomass turnover by fast-growing organisms24,25. Moreover, larger

organisms living in intermittent ecosystems may already be selected to withstand drought12,13. Since top predators are often the largest aquatic organisms in bromeliad ecosystems, they

likely have the highest metabolic demands28, and thus must maximize foraging in a constrained space. Indeed, some bromeliad-living predators are known to display higher prey capture rates

under drought or infrequent rainfall11,29. When rainfall is more infrequent (_k_ < 1), water levels may drop, concentrating resources (as observed by turbidity results) and benefitting

top predators through higher prey encounter rates. This mechanism is not restricted to bromeliad food webs and can be observed in several other intermittent freshwater ecosystems29,30.

Whereas detritivore standing stock exhibited strong site-specific responses to extremes of rainfall frequency, mesopredators were not affected by rainfall or hydrological stability. The

patterns observed for detritivores may be explained by the different ambient rainfall patterns among sites12,25,27, with indications that communities from some regions are more resistant to

drought and flooding than others26,27. Conversely, mesopredators are known to be opportunistic species with generalized feeding and habitat requirements and are less sensitive to habitat

features and climate change24,31,32. Thus, it seems that rainfall extremes modulate different foraging strategies among consumers of different trophic levels. Further experiments should

determine if consumption rates and desiccation-related mortality are both negatively related to the rainfall frequency gradient. Top-heavy biomass pyramids and higher predator:prey body size

ratios consistently became more common under infrequent rainfall scenarios. Two major pathways likely generated top-heaviness in our system: (i) increased energy transfer across trophic

boundaries within the ecosystem (endogenous pathway) and (ii) increased energy transfer across ecosystem boundaries (exogenous pathway)21. Endogenous pathways were evidenced by increasing

ecosystem productivity (resource concentration, supplementary Fig. 4) and prey encounter rates11,29,30 under infrequent rainfall. Likewise, exogenous pathways may promote top-heaviness when

biomass turns over more slowly as trophic rank increases22,28. Top predators in these ecosystems have long lifecycles (>1 year) relative to their prey (often <1 month for _Culicidae_

and _Chironomidae_). Therefore, standing stock of prey may be low, but they normally replenish rapidly in the face of predation pressure through high rates of oviposition by terrestrial

adult flies (with such eggs representing terrestrial subsidies24,25,33,34). However, infrequent rains reduce the window for oviposition, as insect eggs in this system generally need to be

kept humid to survive10. This reduction in oviposition quickly affects the standing stock of the fast-turnover prey, leading to inverted pyramids, whereas the long-lived predators may

survive from an influx of terrestrial prey. Other experiments that reduce oviposition rates have shown similar inversion of biomass pyramids35. More generally, as inverted trophic pyramids

and higher predator:prey body size ratios are both associated with greater interaction strength and unstable predator–prey dynamics20,21,36, these findings indicate that extreme reductions

in rainfall frequency have the potential to destabilize food webs. We show how extremes of precipitation affected each trophic level in a different manner. There were remarkably consistent

changes in the standing stock of top predators and trophic biomass pyramids in response to altered frequency of rainfall across all sites, despite site-specific differences in detritivore

biomass. This implies that similar processes from higher trophic levels may buffer inherent variability within lower trophic levels, and drive consistent food web patterns across the

continental scale. In contrast, detritivores were confined to more stable hydrological regimes, greatly affected by rainfall frequency but with strong site-specific contingency. Detritivores

either benefited or were impaired by extreme rainfall frequency, whereas large predators mostly benefited from infrequent rainfall. This implies that organisms from lower trophic levels may

be more susceptible to rainfall changes in certain geographic regions, and may be restricted to the most stable ecosystems in the near future. As the manipulated changes in extreme rainfall

were within predictions for the next few decades1,2,3,4, we can predict strong changes in food web structure and dynamics in the near future5,6,7,8,15,37. In transient water bodies, such as

small streams, pools, ponds or phytotelmata, this may intensify trophic interaction and result in less stable ecosystems21. METHODS STUDY SITES AND EXPERIMENTAL COMMUNITIES We performed a

geographically coordinated experiment5, manipulating both the amount and temporal distribution of rainfall entering tank bromeliads. Our goal was to investigate the effects of variation in

rainfall on macroinvertebrate communities (standing stock at three adjacent trophic levels and the shape of trophic pyramids). We replicated the experiment at seven sites across Central and

South America (from 29°S to 18°N), including Las Gamas (Argentina), Cardoso and Macae (Brazil), Colombia, Pitilla (Costa Rica), French Guiana, and Puerto Rico. Detailed site descriptions,

the experimental manipulations, and complete list of macroinvertebrate families composing each functional group are provided in ref. 27 and Supplementary Table 3. While the taxonomic

compositions of macroinvertebrate communities are site specific, all taxa can be assigned to three trophic levels: detritivores, mesopredators, and top-predators24,25. The detritivores are

typically represented by larvae of Diptera (_Chironomidae_ (except _Tanypodinae_), _Culicidae_ (except _Toxorhynchites_), _Syrphidae_, _Tipulidae_), and _Coleoptera_ (_Scirtidae_). The

largest top predators are represented by larvae of damselfly (_Coenagrionidae_), horsefly (_Tabanidae_), adult Coleoptera (_Dytiscidae_), and leeches (_Hirudinea_). The mesopredators often

include larvae of _Ceratopogonidae_ (_Bezzia_ spp.), _Corethrellidae_, _Toxorhynchites_ (_Culicidae_) and _Tanypodinae_ (_Chironomidae_). However, the large predators are not present at all

sites25. Thus, we defined top predators on a site-specific basis, as the species without aquatic predators themselves. Consequently, mesopredators at some sites can act as top predators in

other sites (e.g., _Toxorhynchites_ spp., _Corethrella_ spp.). RAINFALL MANIPULATION We established the experimental rainfall treatments for each site using procedures described in ref. 38,

Supplementary Note 1 (Supplementary Table 1) and in ref. 27. Briefly, we used daily rainfall data from the last five years at each site to calculate the site-specific number of days on which

a given amount of rainfall was recorded and fit a negative binomial distribution described by the parameters _μ_ and _k_38. The parameter _μ_ represents the mean daily amount of rainfall

and the parameter k represents the frequency of rainfall events around this mean through time (a measure of evenness in the frequency distribution of rainfall). As climate change affects

individual sites relative to current conditions, our manipulations of precipitation are intentionally proportional to current site conditions and maintain the temporal autocorrelation

structure of each site. The treatment combinations spanned a 30-fold change in _μ_ and a 4-fold change in _k_ around the ambient conditions, using a response surface design with 30 unique

_μ_ by _k_ combinations (Fig. 1b). The “μ1k1” represents ambient treatment combination (control), while the other treatment combinations were derived by multiplying the control values of _μ_

by 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 2.5, and 3, and the control values of _k_ by 0.5, 1 and 2 (Fig. 1b). The range of experimentally imposed values of _μ_ and _k_ were generally larger

than recently observed values in the sites over the experimental months27. Thus, extreme values in these ranges represent extreme rainfall conditions. The experiment lasted for 60 days.

EXPERIMENTAL SETUP AND SAMPLING In each site, we selected thirty bromeliads of the most abundant species and with the most common size27. We used bromeliads that had more than 100 ml of tank

capacity and thus can be colonized by the large predators. We washed each bromeliad with spring water to remove detritus and organisms; we retained coarse (>850 μm) and fine (<850 μm)

detritus and sorted macroinvertebrates into species groups. To remove any residual invertebrates, we hung the bromeliads upside down, and let them dry for 7 days. Then, to initiate the

community assembly in the experimental ecosystems we evenly divided the fine and coarse detritus between the 30 bromeliads and stocked each bromeliad with the same community in terms of

invertebrate families and functional groups27. We employed individual transparent plastic shelters above each bromeliad to prevent natural rainfall into the plants. The rain shelters were

settled high enough to ensure that they did not alter macroinvertebrate colonization or temperature within the bromeliads. We randomly divided the 30 treatment combinations into three blocks

of 10 bromeliads, and initiated each block on one of three consecutive days. This procedure also allowed enough time to sample invertebrates at the end of the experiment. In order to

estimate the key hydrological parameters15, we measured water depth in the central and two lateral leaf wells of each bromeliad every two days, and used average values per bromeliad. The

hydrological parameters for each bromeliad included: (i) coefficient of variation of water depth across the entire experiment, (ii) proportion of overflow days, i.e., the number of days

water depth was ≥maximum depth recorded divided by the total number of measurements, (iii) proportion of dried-out days, i.e., the number of days water depth was <5 mm divided by the

total number of measurements. These hydrological parameters were used to create a metric of hydrological stability, using Principal Component Analysis (PCA)15. We used the scores of the

first axis of the PCA (Supplementary Table 2), for each site, to summarize these parameters into a single variable of hydrological stability. This axis quantified a gradient of habitat

permanence and stability, where increasing scores represent more stable ecosystems (i.e., ecosystems that dried out less often, and held more water throughout the experiment, Supplementary

Table 2). At the end of the experiment (60th day), we recorded with hand-held data loggers water turbidity, a measure that integrates organic and inorganic suspended matter, including

free-living algae and particulate nutrients (C, N, P), resulting from detritivore activity. Thus, turbidity represents a surrogate for total nutrient availability in freshwater

ecosystems39,40. Then, we dissected each bromeliad by removing and washing each leaf separately in running water and then filtered this water through 125 and 850 µm sieves. We recorded the

morphospecies and abundance of all aquatic macroinvertebrates (body size larger than 0.5 mm). We recorded 4–38 morphospecies per bromeliad (mean per bromeliad ± SD: 14.7 ± 9.4). We

determined the body size and trophic position of each individual organism surveyed. Trophic position was determined from our own feeding trials, gut contents, stable isotope analyses and

from the literature41,42,43. To calculate invertebrate body mass, we used allometric equations between the body length and dry mass, or mean of dry mass for very small insects. STATISTICAL

ANALYSES We tested the independent and interactive effects of site (categorical, seven levels), ecosystem size (maximum volume, continuous), rainfall evenness (k, categorical, three levels),

quantity (μ, continuous) and hydrology (PCA axis 1, continuous), on (i) the standing stock of the three trophic levels, (ii) the structure of trophic pyramids (predator-detritivore biomass

ratio), (iii) the predator-detritivore body size ratios, and (iv) ecosystem productivity (turbidity). We used a negative binomial response distribution, which best represents the variation

in these responses27. Site was treated as a fixed effect in all the models. To account for non-linear trends, we included quadratic functions for all the continuous predictors. Additive and

interactive models were analysed using type II and III sum of squares (SS), respectively. Probabilities were calculated using likelihood ratio tests (LRT, _χ_2). All the models were reduced

using backward selection, and only the final, simplified models are presented. Model simplification consisted of removing more complex non-significant interactions (third order), followed by

more simple non-significant interactions (second order), and non-significant quadratic functions44. Contingent and consistent responses across all sites were tested as an interaction

(linear and quadratic functions) with site, where non-significant interactions indicated consistency across the sites. In these cases, the models were fitted using generalized linear mixed

models (GLMM), with site as random factor, for graphical modeling45. The data are graphically displayed using the final, reduced (after backward selection) models. All the analyses were

performed using R 3.2.246. The established significance level was _α_ = 0.05. We checked variance heterogeneity, normality, and outliers through graphic inspections to assure the model

assumptions were met. We used the functions _rda_ (_vegan_ package) and _prcomp_ (_stats_ package) to perform the PCA (hydrological stability). This statistical method differed from a

previous analysis of functional feeding groups27 in that ecosystem size was allowed to leave the model, and backwards selection instead of AIC model selection was employed, so we do not

expect results to be completely comparable between the two studies. REPORTING SUMMARY Further information on research design is available in the Nature Research Reporting Summary linked to

this article. DATA AVAILABILITY The data that support the findings of this study are available at https://doi.org/10.5281/zenodo.1124951. This data was collated and hydrologic metrics

calculated by a custom-built R package, BWGTools, available at: https://doi.org/10.5281/zenodo.1120418. CODE AVAILABILITY The R code used to calculate the precipitation treatments is

publically archived at https://doi.org/10.5281/zenodo.18548. REFERENCES * Pendergrass, A. G. What precipitation is extreme? _Science_ 360, 1072–1073 (2018). Article  ADS  CAS  Google Scholar

  * Fischer, E. M. & Knutti, R. Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes. _Nat. Clim. Change_ 5, 560–564 (2015). Article  ADS

  Google Scholar  * Donat, M. G., Lowry, A. L., Alexander, L. V., O’Gorman, P. A. & Maher, N. More extreme precipitation in the world’s dry and wet regions. _Nat. Clim. Change_ 6,

508–513 (2016). Article  ADS  Google Scholar  * _IPCC Climate Change 2013: The Physical Science Basis_ (eds Stocker, T. F. et al) (Cambridge University Press, 2014). * Knapp, A. K. et al.

Pushing precipitation to the extremes in distributed experiments: recommendations for simulating wet and dry years. _Glob. Change Biol._ 23, 1774–1782 (2016). Article  ADS  Google Scholar  *

Woodward, G. et al. The effects of climatic fluctuations and extreme events on running water ecosystems. _Philos. Trans. R. Soc. B_ 371, 20150274 (2016). Article  Google Scholar  * Ledger,

M. E., Brown, L. E., Edwards, F. K., Milner, A. M. & Woodward, G. Drought alters the structure and functioning of complex food webs. _Nat. Clim. Change_ 3, 223–227 (2013). Article  ADS 

Google Scholar  * Lu, X. et al. Drought rewires the cores of food webs. _Nat. Clim. Change_ 6, 875–878 (2016). Article  ADS  Google Scholar  * McHugh, P. A., Thompson, R. M., Greig, H. S.,

Warburton, H. J. & McIntosh, A. R. Habitat size influences food web structure in drying streams. _Ecography_ 38, 700–712 (2015). Article  Google Scholar  * Dézerald, O., Céréghino, R.,

Corbara, B., Dejean, A. & Leroy, C. Functional trait responses of aquatic macroinvertebrates to simulated drought in a Neotropical bromeliad ecosystem. _Freshw. Biol._ 60, 1917–1929

(2015). Article  Google Scholar  * Amundrud, S. L. & Srivastava, D. S. Trophic interactions determine the effects of drought on an aquatic ecosystem. _Ecology_ 97, 1475–1483 (2016).

Article  Google Scholar  * Trzcinski, M. K. et al. The effects of food web structure on ecosystem function exceeds those of precipitation. _J. Anim. Ecol._ 85, 1147–1160 (2016). Article 

Google Scholar  * Dorn, N. J. & Cook, M. I. 2015. Hydrological disturbance diminishes predator control in wetlands. _Ecology_ 96, 2984–2993 (2015). Article  Google Scholar  *

González-Megías, A. & Menéndez, R. Climate change effects on above- and below-ground interactions in a dryland ecosystem. _Philos. Trans. R. Soc. B_ 367, 3115–3124 (2012). Article 

Google Scholar  * Marino, N. A. C. et al. Rainfall and hydrological stability alter the impact of top predators on food web structure and function. _Glob. Change Biol._ 23, 673–685 (2017).

Article  ADS  Google Scholar  * Stevnbak, K. et al. 2012. Interactions between above- and belowground organisms modified in climate change experiments. _Nat. Clim. Change_ 2, 805–808 (2012).

Article  ADS  CAS  Google Scholar  * Romo, C. M. & Tylianakis, J. M. Elevated temperature and drought interact to reduce parasitoid effectiveness in suppressing hosts. _PLoS ONE_ 8,

e58136 (2013). Article  ADS  CAS  Google Scholar  * Liu, S. et al. Spider foraging strategy affects trophic cascades under natural and drought conditions. _Sci. Rep._ 5, 12396 (2015).

Article  ADS  Google Scholar  * Suttle, K. B., Thomsen, M. A. & Power, M. E. Species interactions reverse grassland responses to changing climate. _Science_ 315, 640–642 (2007). Article

  ADS  CAS  Google Scholar  * Gilbert, B. et al. A bioenergetic framework for the temperature dependence of trophic interactions. _Ecol. Lett._ 17, 902–914 (2014). Article  Google Scholar  *

McCauley, D. J. et al. On the prevalence and dynamics of inverted trophic pyramids and otherwise top-heavy communities. _Ecol. Lett._ 21, 439–454 (2018). Article  Google Scholar  * Hatton,

I. A. et al. The predator-prey power law: biomass scaling across terrestrial and aquatic biomes. _Science_ 349, aac6284 (2015). Article  Google Scholar  * O’Gorman, E. J. et al. Unexpected

changes in community size structure in a natural warming experiment. _Nat. Clim. Change_ 7, 659–663 (2017). Article  ADS  Google Scholar  * Romero, G. Q., Piccoli, G. C. O., de Omena, P. M.

& Gonçalves-Souza, T. Food web structure shaped by habitat size and climate across a latitudinal gradient. _Ecology_ 97, 2705–2715 (2016). Article  Google Scholar  * Petermann, J. S. et

al. Dominant predators mediate the impact of habitat size on trophic structure in bromeliad invertebrate communities. _Ecology_ 96, 428–439 (2015). Article  Google Scholar  *

Benavides-Gordillo, S., Farjalla, V. F., González, A. L. & Romero, G. Q. Changes in rainfall level and litter stoichiometry affect aquatic community and ecosystem processes in bromeliad

phytotelmata. _Freshw. Biol._ 64, 1357–1368 (2019). Article  CAS  Google Scholar  * Srivastava, D. et al. Ecological response to altered rainfall differs across the Neotropics. _Ecology_.

101, e02984 (2020). Article  Google Scholar  * Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Towards a metabolic theory of ecology. _Ecology_ 85, 1771–1789

(2004). Article  Google Scholar  * Amundrud, S. L. et al. Drought alters the trophic role of an opportunistic generalist in an aquatic ecosystem. _Oecologia_ 189, 733–744 (2019). Article 

ADS  Google Scholar  * Dalal, A., Cuthbert, R. N., Dick, J. T. A. & Gupta, S. Water depth-dependent notonectid predatory impacts across larval mosquito ontogeny. _Pest Manag. Sci._ 75,

2610–2617 (2019). Article  CAS  Google Scholar  * Brose, U. et al. Climate change in size-structured ecosystems. _Philos. Trans. R. Soc. B_ 367, 2903–2912 (2012). Article  Google Scholar  *

Prugh, L. R. et al. The rise of the mesopredator. _BioScience_ 9, 779–790 (2009). Article  Google Scholar  * de Omena, P. M., Srivastava, D. S. & Romero, G. Q. Consumptive effects and

mismatch in predator–prey turnover rates cause inversion of biomass pyramids. _Oecologia_ 190, 159–168 (2019). Article  ADS  Google Scholar  * Srivastava, D. S. & Lawton, J. H. Why more

productive sites have more species: an experimental test of theory using tree‐hole communities. _Am. Naturalist_ 152, 510–529 (1998). Article  CAS  Google Scholar  * LeCraw, R. M.,

Srivastava, D. S. & Romero, G. Q. Metacommunity size influences aquatic community composition in a natural mesocosm landscape. _Oikos_ 123, 903–911 (2014). Article  Google Scholar  *

Emmerson, M. C. & Raffaelli, D. Predator–prey body size, interaction strength and the stability of a real food web. _J. Anim. Ecol._ 73, 399–409 (2004). Article  Google Scholar  *

Amundrud, S. L. & Srivastava., D. S. Disentangling how climate change can affect an aquatic food web by combining multiple experimental approaches. _Glob. Change Biol._ 25, 3528–3538

(2019). Article  ADS  Google Scholar  * MacDonald, A. A. M. & Srivastava D. S. Precipitation schedules for the Bromeliad Working Group rainfall experiment [Data set]. _Zenodo_.

https://doi.org/10.5281/zenodo.18548 (2015). * Kim, J. & Furumai, H. Improved calibration of a rainfall-pollutant-runoff model using turbidity and electrical conductivity as surrogate

parameters for total nitrogen. _Water Environ. J._ 27, 79–85 (2013). Article  CAS  Google Scholar  * Lannergård, E. E. et al. An evaluation of high frequency turbidity as a proxy for

riverine total phosphorus concentrations. _Sci. Total Environ._ 651, 103–113 (2019). Article  ADS  Google Scholar  * Kitching, R. L. _Food Webs and Container Habitats: The Natural History

and Ecology Of Phytotelmata_. (Cambridge University Press, 2000). * Merritt, R. W. & Cummins, K. W. _An Introduction to the Aquatic Insects of North America_. (Kendall Hunt Publishing

Company, Dubuque, 1996) . * Dézerald, O. et al. Functional traits and environmental conditions predict community isotopic niches and energy pathways across spatial scales. _Funct. Ecol_. 32,

2423–2434. * Crawley, M. J. _The R Book_, 2nd edn. (John Wiley & Sons, 2013). * Breheny, P. and Burchett, W. Visualizing regression models using visreg. http://myweb.

uiowa.edu/pbreheny/publications/visreg.pdf (2013) * R Core Team. 2019. _R: A language and environment for statistical computing_. (R Foundation for Statistical Computing, Vienna, 2019)

https://www.Rproject.org/ Download references ACKNOWLEDGEMENTS We acknowledge financial support for research provided by the São Paulo Research Foundation (FAPESP: grants 2012/51143-3,

2017/09052-4, and BPE Proc. 2016/01209-9) and CNPq through research grants to G.Q.R, by the Agence Nationale de la Recherche (ANR), Rainwebs project (grant ANR-12-BSV7-0022-01) to R.C.,

C.L., B.C., J.F.C. and D.S.S. and an Investissement d’Avenir grant (Labex CEBA, ref. ANR-10-LABX-25-01), by the Agencia Nacional de Promoción Científica y Tecnológica (grant PICT-2010-1614)

and Secretaría de Ciencia y Tecnología de la Universidad Nacional de Rosario (grant AGR-139) through grants to I.B. and G.M., by the Brazilian Council for Research, Development and

Innovation (CNPq) for research funds (Pesquisador Visitante Especial, PVE, Research Grant 400454/2014-9) and productivity grants to V.F.F., by the Natural Sciences and Engineering Research

Council of Canada (NSERC) through Discovery and Accelerator grant to D.S.S., and the Facultad de Ciencias, Universidad de los Andes, Colombia through a grant (2012-1) to F.O. This research

was further supported by scholarship and fellowship support from the ANR to M.K.T., from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) to R.F., from FAPESP to P.A.P.A.

(Proc. 2014/04603-4, 2017/26243-8), from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) to A.B.A.C, J.S.L., and N.A.C.M. (PNPD-CAPES 2013/0877) and P.M.O. (PNPD-CAPES

2014/04603-4), an Investissement d’Avenir grant (Labex CEBA, ref. ANR-10-LABX-25-01) and a PhD fellowship from the Fond Social Européen to O.D., from NSERC to A.A.M.M. (CGS-D) and D.S.S.

(EWR Steacie Memorial Fellowship), from the University of British Columbia to S.L.A., and from the Departamento Administrativo de Ciencia, Tecnología e Innovación (COLCIENCIAS) support of

F.O.B. (COLCIENCIAS grant 567). G.Q.R. and P.K. gratefully acknowledge funding from the Royal Society, Newton Advanced Fellowship (grant no. NAF/R2/180791). E.J.O’G. gratefully acknowledge

funding from NERC (NE/L011840/1). This paper is a contribution of the Brazilian Network on Global Climate Change Research funded by CNPq (grant #550022/2014-7) and FINEP (grant

#01.13.0353.00). AUTHOR INFORMATION Author notes * A. Andrew M. MacDonald Present address: Université de Montreal, Montreal, Québec, Canada * Paula M. de Omena Present address: Institute of

Biological Sciences, Universidade Federal do Pará, Belém, PA, Brazil AUTHORS AND AFFILIATIONS * Laboratory of Multitrophic Interactions and Biodiversity, Department of Animal Biology,

Institute of Biology, University of Campinas (UNICAMP), Campinas, SP, 13083-862, Brazil Gustavo Q. Romero, Pablo A. P. Antiqueira & Paula M. de Omena * Programa de Pós-Graduação em

Ecologia, Universidade Federal do Rio de Janeiro (UFRJ), CP 68020, Rio de Janeiro, RJ, Brazil Nicholas A. C. Marino, Juliana S. Leal & Alice B. A. Campos * Departamento de Ecologia,

Instituto de Biologia, Centro de Ciências da Saúde, Universidade Federal do Rio de Janeiro, PO Box 68020, Rio de Janeiro, RJ, Brazil Nicholas A. C. Marino & Vinicius F. Farjalla * Centre

for the Synthesis and Analysis of Biodiversity (CESAB-FRB), Aix-en-Provence, France A. Andrew M. MacDonald * Laboratoire Ecologie Fonctionnelle et Environnement, Université de Toulouse,

CNRS, Toulouse, France Régis Céréghino * Department of Forest and Conservation Sciences, University of British Columbia, Vancouver, BC, Canada M. Kurtis Trzcinski * Department of Biology,

University of Puerto Rico - Mayagüez Campus, Mayagüez, PR, 00681, USA Dimaris Acosta Mercado * AMAP, Univ. Montpellier, CIRAD, CNRS, INRAE, IRD, Montpellier, France Céline Leroy * UMR

ECOFOG, CIRAD, CNRS, INRAE, AgroParisTech, Université de Guyane, Université des Antilles, 97379, Kourou, France Céline Leroy * Université Clermont-Auvergne, CNRS, LMGE (Laboratoire

Microorganismes: Génome et Environnement), F-63000, Clermont-Ferrand, France Bruno Corbara & Jean-François Carrias * Facultad de Ciencias Agrarias, Instituto de Investigaciones en

Ciencias Agrarias de Rosario, IICAR-CONICET-UNR, Universidad Nacional de Rosario, S2125ZAA, Zavalla, Argentina Ignacio M. Barberis, Guillermo Montero & Rodrigo Freire * ESE, Ecology and

Ecosystem Health, INRAE, Agrocampus Ouest, 35042, Rennes, France Olivier Dézerald * Department of Watershed Sciences and the Ecology Center, Utah State University, Logan, 84322, USA Edd

Hammill & Trisha B. Atwood * Department of Zoology and Botany, University of São Paulo State (UNESP/IBILCE). 15054 - 000, São José do Rio Preto, SP, Brazil Gustavo C. O. Piccoli *

Departmento de Ciencias Biológicas, Universidad de los Andes, Bogotá, 111711, Colombia Fabiola Ospina Bautista & Emilio Realpe * Departmento de Ciencias Biológicas, Universidad de

Caldas, Manizales, 170004, Colombia Fabiola Ospina Bautista * Department of Zoology & Biodiversity Research Centre, University of British Columbia, Vancouver, BC, V6T 1Z4, Canada Sarah

L. Amundrud & Diane S. Srivastava * Queen Mary University of London, School of Biological and Chemical Sciences, London, UK Pavel Kratina * School of Life Sciences, University of Essex,

Colchester, UK Eoin J. O’Gorman Authors * Gustavo Q. Romero View author publications You can also search for this author inPubMed Google Scholar * Nicholas A. C. Marino View author

publications You can also search for this author inPubMed Google Scholar * A. Andrew M. MacDonald View author publications You can also search for this author inPubMed Google Scholar * Régis

Céréghino View author publications You can also search for this author inPubMed Google Scholar * M. Kurtis Trzcinski View author publications You can also search for this author inPubMed 

Google Scholar * Dimaris Acosta Mercado View author publications You can also search for this author inPubMed Google Scholar * Céline Leroy View author publications You can also search for

this author inPubMed Google Scholar * Bruno Corbara View author publications You can also search for this author inPubMed Google Scholar * Vinicius F. Farjalla View author publications You

can also search for this author inPubMed Google Scholar * Ignacio M. Barberis View author publications You can also search for this author inPubMed Google Scholar * Olivier Dézerald View

author publications You can also search for this author inPubMed Google Scholar * Edd Hammill View author publications You can also search for this author inPubMed Google Scholar * Trisha B.

Atwood View author publications You can also search for this author inPubMed Google Scholar * Gustavo C. O. Piccoli View author publications You can also search for this author inPubMed 

Google Scholar * Fabiola Ospina Bautista View author publications You can also search for this author inPubMed Google Scholar * Jean-François Carrias View author publications You can also

search for this author inPubMed Google Scholar * Juliana S. Leal View author publications You can also search for this author inPubMed Google Scholar * Guillermo Montero View author

publications You can also search for this author inPubMed Google Scholar * Pablo A. P. Antiqueira View author publications You can also search for this author inPubMed Google Scholar *

Rodrigo Freire View author publications You can also search for this author inPubMed Google Scholar * Emilio Realpe View author publications You can also search for this author inPubMed 

Google Scholar * Sarah L. Amundrud View author publications You can also search for this author inPubMed Google Scholar * Paula M. de Omena View author publications You can also search for

this author inPubMed Google Scholar * Alice B. A. Campos View author publications You can also search for this author inPubMed Google Scholar * Pavel Kratina View author publications You can

also search for this author inPubMed Google Scholar * Eoin J. O’Gorman View author publications You can also search for this author inPubMed Google Scholar * Diane S. Srivastava View author

publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS G.Q.R., D.S.S., R.C., V.F.F., F.O.B., I.M.B., J-F.C., A.A.M.M., B.C., D.A.M., C.L., N.A.C.M., and

G.C.O.P. conceived the idea, designed the global experiment, and developed it with M.K.T., O.D., E.H., T.B.A., J.S.L., G.M., P.A.P.A., R.F., S.L.A., P.M.O., and A.B.A.C. D.S.S., A.A.M.M.,

and N.A.C.M. collated the data and developed hydrological metrics. G.Q.R. drafted the paper with inputs from all authors. G.Q.R. and N.A.C.M. performed the statistical analyses and drafted

the figures. G.Q.R., P.K., D.S.S. and E.J.O’G. interpreted the results. A.A.M.M. and D.S.S. created rainfall schedules, G.Q.R., D.S.S., I.M.B., F.O.B., R.C., V.F.F., E.H., D.A.M., G.M., and

E.R. organized site-level experiments. All authors participated in discussing and editing the paper. CORRESPONDING AUTHOR Correspondence to Gustavo Q. Romero. ETHICS DECLARATIONS COMPETING

INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PEER REVIEW INFORMATION _Nature Communications_ thanks Alonso Ramirez and the other, anonymous, reviewers for

their contribution to the peer review of this work. Peer reviewer reports are available. PUBLISHER’S NOTE Springer Nature remains neutral with regard to jurisdictional claims in published

maps and institutional affiliations. SUPPLEMENTARY INFORMATION SUPPLEMENTARY INFORMATION PEER REVIEW FILE REPORTING SUMMARY RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under

a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate

credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article

are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and

your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this

license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Romero, G.Q., Marino, N.A.C., MacDonald, A.A.M. _et al._ Extreme

rainfall events alter the trophic structure in bromeliad tanks across the Neotropics. _Nat Commun_ 11, 3215 (2020). https://doi.org/10.1038/s41467-020-17036-4 Download citation * Received:

23 January 2020 * Accepted: 03 June 2020 * Published: 25 June 2020 * DOI: https://doi.org/10.1038/s41467-020-17036-4 SHARE THIS ARTICLE Anyone you share the following link with will be able

to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing

initiative