Litterfall Chemistry Is Modulated by Wet-Dry Seasonality and Leaf Phenology of Dominant Species in the Tropics

Litterfall has a large influence on carbon and nutrient cycling of ecosystems, particularly in light-limited forested streams, as most nutrients return in the form of litter. Although recent evidence points to the prevalence of seasonal litterfall in species-rich and evergreen tropical riparian forests, there is a limited understanding of how riparian plant diversity intersects with stream and riparian ecosystem functions. To explore this question, we investigate litterfall chemistry across wet and dry seasons and the congruence between litter traits and plant species composition of litterfall in the wet-dry tropics. Using generalized additive models, we observed consistent seasonal patterns of litterfall chemistry over 2 years, mostly influenced by dominant species in litterfall. While drier seasons showed litter lower in nutrients and structural compounds and higher in polyphenols, litter from wetter seasons were nutrient rich but lower in polyphenols. We also found contrasting seasonal patterns in litterfall chemistry, one showing that litterfall nutrient, structural compounds, and secondary metabolite concentrations declined in drier seasons while the other showed that mass-based litterfall inputs increased markedly in drier seasons. Our findings suggest that litterfall chemistry may be altered by shifts in the identity of dominant plant species and seasonality, possibly leading to changes in carbon and nutrient fluxes in tropical riparian ecosystems.


INTRODUCTION
Riparian zones are interfaces between forest and freshwater ecosystems tightly linked by the flux of energy and nutrients (Naiman et al., 2005). This connection is especially important to forested headwater streams that are light-limited ecosystems, in which the supply of terrestrial litter is the main energetic source for complex food webs (Wallace et al., 1997;Neres-Lima et al., 2017). Streams play a particularly important role in receiving, transporting and processing the litter, contributing to nutrient and carbon fluxes at local, regional and global scales (Raymond et al., 2013;Hotchkiss et al., 2015). Most litter enters streams in the form of senescent leaves-although green leaves may also be found-and through direct fall (hereafter litterfall) or laterally from the forest soil, however, litterfall is often the main route (>70% of total litter inputs; Tonin et al., 2017). Consequently, the amount of leaf litter entering streams depend on the timing of leaf senescence of multiple plant species present in the riparian forest (Reich, 1995;Gonçalves and Callisto, 2013).
Plants change the timing of leaf senescence according to environmental cues (e.g., temperature, precipitation, photoperiod), which vary in importance depending on species, biomes and environmental constraints (Tang et al., 2016). For example, in temperate deciduous and coniferous forests-where the majority of studies were done-most of leaf litterfall occurs in autumn as day length shortens and temperature cools (Abelho, 2001). In contrast, tropical forests-which encompass rain, dry and peat swamp forests, and savannas-are often assumed to have lower (or even unpredictable) seasonality of litterfall because of lower temperature and photoperiod seasonality. However, a recent study showed that seasonal litterfall in the tropics is consistently (and negatively) associated with precipitation, even in rainforests, which experiences no dry season ; but see also Zhang et al., 2014;Nakagawa et al., 2019). In the wet-dry tropics-which include the Cerrado biome in Brazil as well as southern Amazon, central Africa and northern Australia-, most of annual precipitation occurs in certain months of the year (Alvares et al., 2013). Consequently, massive litterfall is produced in the driest periods as a response of plants to reduce water stress (Reich, 1995). Although litterfall has a large influence on carbon and nutrient cycling in forests and streams because most nutrients return in the form of leaf litter (Prescott, 2002), there still are a number of knowledge gaps-especially in the understudied tropical region-including (i) the chemistry of leaf litterfall, and (ii) how it varies year-round and (iii) whether the diversity of leaf phenological patterns in plant communities correlates with the chemistry of litterfall.
Litterfall chemistry is expected to vary year-round because of several mechanisms related to species physiology and environmental controls acting on leaf senescence and/or directly on freshly fallen litter. First, plants are able to resorb nutrients from senescing leaves back into live plant tissues allowing the conservation of nutrients that otherwise would be lost to soils via litterfall (Aerts and Chapin, 1999). Nutrient resorption is then responsible for reducing the concentrations of nitrogen and phosphorus in freshly fallen litter and thus, drive changes in C:N:P balance, which has proven to be a powerful tool for understanding biogeochemical cycling in ecosystems (Elser et al., 2000). Lower soil fertility and drought conditions, both occurring in the Brazilian Cerrado, increase nutrient resorption (Brant and Chen, 2015). Therefore, most leaf litter from riparian forests of Cerrado are nutrient poor with high C:N and C:P mass ratios. Second, polyphenols, that represents the most abundant class of plant secondary compounds, have limited resorption during senescence and its production tends to be stimulated by climatic stress resulting from either drought or increased temperature (Tharayil et al., 2011). Finally, litter chemistry may also be the result of different plant species losing their senescent leaves at different periods of the year. For example, evergreen, semi-deciduous and deciduous species, all occurring in tropical riparian forests, may present distinct patterns of leaf senescence and of annual litter production. While litter of evergreen species may fall year-round or episodically, litter from deciduous species is mostly found in specific seasons (Reich, 1995). Thus, interspecific variation in the timing of litterfall may have important implications for the amounts and composition of carbon and nutrients input to riparian ecosystems, which will ultimately determine the rate of several key ecosystem processes, e.g., litter decomposition, primary production, nutrient cycling.
Our study explored the variation of litter chemical traits across wet and dry seasons, and the congruence between chemical traits and species composition of litter entering streams of the Cerrado biome. We quantified chemical traits related to litter structural compounds, nutrients, secondary metabolites as well as C:N:P stoichiometric mass ratios and the contribution of individual species to litterfall by sampling streams of the Cerrado biome over 2 years. We tested the following hypotheses: (1) Litter is nutrient (N and P) poor in drier periods because nutrient resorption occurs before leaf senescence and is a common mechanism of plants to avoid losing important and scarce nutrients. Litter has a greater concentration of secondary metabolites in drier periods because of reduced resorption and drought stress before leaf abscission. We have no a priori expectation about the variability of structural compounds of litter across drier and wetter periods.
(2) Despite the enormous diversity of species that contribute to litterfall in the tropics, litter chemical traits are mostly congruent with the composition of dominant plant species because they comprised the majority of litter biomass.

Study Sites
The study was conducted in three riparian ecosystems encompassing streams and their surrounding vegetation ( (Eiten, 1972). Streams are similar in size (mean ± SE: wetted channel width, 2.5 ± 0.3 m; water depth: 38 ± 12 cm), with a channel slope of 3.3 ± 0.9 • and containing a mixture of gravel and silt as streambed substrates.
Streams flow through dense (84-92% forest cover) strips of evergreen riparian forests (known as gallery forests) within the savanna (i.e., Cerrado sensu stricto). Riparian forests are evergreen, have trees 20-30 m tall and a dense understory of low stature (Oliveira-Filho et al., 2002). There is also the presence of semi-deciduous and deciduous species (33 and 13% of species, respectively, Tonin et al., 2019). Plant species diversity of phanerogams is often high in riparian forests of Cerrado as shown by a previous floristic study in the same area (68 species at CAP, 71 at CVE and 29 at RON in an area of approx. 0.12 ha, Bambi et al., 2017a). Common species in our sites were Protium spruceanum, Pseudomenia laevigata, and Tapirira obtusa at CAP; P. spruceanum, Matayba guianensis, and Cyathea villosa at CVE; and Xilopia emarginata, Richeria grandis, and Clusia criuva at RON. Complementary information on morphological and physical characteristics of these riparian ecosystems can be found in Bambi et al. (2017a) and Tonin et al. (2019). The studied ecosystems experienced a seasonal tropical climate (Aw, Köppen's climate classification) with a dry season that coincides with the coldest months of the year (June-August) and a wet season (November-March). Additionally, two transitional seasons are evident: the dry-wet (September-October) and the wet-dry (April-May). Annual rainfall is around 1,477 mm and mean air temperature is 21 • C (National Institute of Meteorology (1981-2010); Figure 1). During the 2 years of this study, annual rainfall was 1,418 mm and mean air temperature was 21.1 • C.

Field and Laboratory Work
We analyzed the chemical composition (hereafter chemical traits) and plant species composition of leaf litter entering the three streams. Our study focused on leaf litter because it represented the majority of inputs  and it is a biologically active component in forest streams that is replaced throughout the year (Webster et al., 1999). Litterfall was collected in each stream using suspended traps distributed at 5 equally distanced sampling sites along a 50-100 m reach. Sampling sites were selected to be representative of reaches, mostly avoiding large gaps in forest canopies or large trees. Litter traps were plastic buckets with small holes in the bottom for drainage (26 cm in diameter) that were placed 2m above the streambed. Traps were installed using ropes tied to trees along the streambank and in a perpendicular direction to the stream channel. At each site, a total of 18 litter traps were installed on three ropes, each with six traps, which is equivalent to a sampling effort of approximately 1 m 2 . Ropes were located several centimeters from one another and along each rope traps were equally distributed to cover the stream channel and therefore maximize the capture of litterfall.
Leaf litter (hereafter litter) was collected once a month for 2 years, comprising a total of 24 sampling events (from September 2010 to September 2012) (see Bambi et al., 2017b;Tonin et al., 2019). All litter from a site were pooled together (18 traps) per sampling event. In the laboratory, samples were first air-dried, identified to the species level, oven dried (60 • C for 72 h) and weighed separately per species. Species identification was done only for a subset of samples (i.e., samples of the first year of study collected between September 2010 and August 2011) due to logistical reasons, using the databases of the Tropics projects 1 , the Flora do Brazil species list 2 , the Field Museum 3 and the Link species network 4 , based on the APG III classification system. Secondly, all species from a sample were mixed, oven-dried (60 • C for 72 h), ground in a ball mill to powder and chemically analyzed for the concentration of structural compounds [carbon (C), fibers, lignin and cellulose], nutrients [nitrogen (N) and phosphorus (P)] and secondary metabolites (polyphenols). The concentration of C and N was determined in an elemental analyzer (Leco Corporation-TruSpec Micro CHN628) by total combustion of litter sample at 950 • C and detection of C and N in the form of CO 2 and N 2 , respectively. Fibers, lignin (as proximate lignin) and cellulose were obtained gravimetrically using the acid-detergent fiber procedures proposed by Goering and Van Soest (1970) with slight modifications (Gessner, 2005). Briefly, this method consists in the determination of residual weight of samples after successive extraction of tissue constituents. Phosphorus was quantified as phosphate (PO 4 -P) spectrophotometrically (measuring absorbance at 882 nm) using the ascorbic acid method after combustion (550 • C for 4 h) and digestion (with hydrochloric acid) of litter samples (Flindt and Lillebø, 2005). Polyphenols were determined by the Folin-Ciocalteu method (Folin and Ciocalteu, 1927) which consists in extracting phenol hydroxyl groups from litter samples with 70% acetone for 1 h at 4 • C, reading absorbance at 760 nm in spectrophotometer and determine phenolic concentration in terms of tannic acid equivalents (Bärlocher and Graça, 2005).

Calculation of Litter Traits
Litter traits were expressed as concentrations (in percentage) of each chemical compound or mass ratios of C:N, C:P, and N:P per month at each sampling site. The amount of litter constituents in terms of masses was estimated by multiplying the proportion of each constituent by litterfall mass. Litterfall was calculated as litter dry mass (in grams) divided by the area of traps in a site (i.e., 0.9558 m 2 ) and the elapsed time (in days), multiplied by 30 days to estimate monthly litterfall. The contribution of each species to litterfall was calculated similarly, but to the scale of a riparian ecosystem, as litter dry mass of a species (in grams) divided by the area of traps in an ecosystem (i.e., 4.779 m 2 ) and the elapsed time (in days) multiplied by 30 days.

Data Analyses
All analyses were conducted in R v. 4.0.1 (R Core Team, 2020) using the packages mgcv (Wood, 2011), nlme (Pinheiro et al., 2020), tidyverse (Wickham et al., 2019), and vegan (Oksanen et al., 2020) and figures were prepared using ggplot2 (Wickham, 2016). Initial data exploration revealed potential non-linear patterns of litter chemical traits over time and remarkably similarities of these data across riparian ecosystems,-supporting its inclusion as a random term in the generalized additive mixed models (GAMMs). Further data exploration also indicated (i) different variance over time, which was taken into account while structuring and validating the models and (ii) high variability and the presence of true extreme observations in N:P mass ratio data, which required a log-transformation prior to analysis.
We thus explored temporal dynamics of litter chemical traits in terms of concentrations and mass with GAMMs with time (number of each sampling event) as a continuous predictor and riparian ecosystems as a random term for each chemical compound separately. Models were fitted with time as a smooth term, with a normal distribution of response data, the identity-link function and restricted maximum likelihood method. This type of model (also known as smoothing models) was used because they allow for non-linear relationships of the response variable in relation to the predictor (Zuur et al., 2009). The optimal amount of smoothing was estimated with cross-validation and expressed as effective degrees of freedom (edf). The higher the edf, the more non-linear is the smoother. Variance was allowed to differ among sampling events using a constant variance function structure (VarIdent) to meet the homogeneity assumption (Zuur and Ieno, 2015). The need for this term was identified in initial data exploration and confirmed by comparison of the Akaike Information Criterion (AIC) of models with and without this component. Neither spatial autocorrelation within and among streams nor temporal autocorrelation between subsequent sampling events were detected. We checked for spatial and temporal autocorrelation using variograms and the autocorrelation function of normalized model residuals, respectively. Graphical diagnostics were used as validation tools for the models: QQ-plots and histograms of residuals to assess normality and residuals versus fitted values for homogeneity.
We investigated the congruence between plant species composition of litterfall and litter chemical traits over time with redundancy analysis (RDA). Species data was structured to represent one observation of each riparian ecosystem at each sampling event per species. Therefore, litter chemical traits were averaged across sites per riparian ecosystem per sampling event. We first fitted an RDA model using all plant species data as a response matrix (Hellinger-transformed) and litter traits as a matrix of predictors (z-score transformed). We then fitted a second RDA model using only the dominant species to test whether the litter trait variability is more congruent with the phenology of the dominant species. Four non-colinear litter traits (i.e., carbon, C:N mass ratio, cellulose and lignin) were retained and used in the matrix of predictors after evaluation of pairwise correlations between traits (Supplementary Table 2) and adequacy of variance inflation factors of RDA. Statistical significance of the global RDA models and of its axes were obtained by permutation tests using 999 iterations.

Variation of Litter Chemistry Across Seasons
All litter traits related to structure, nutrients (except N:P mass ratio) and secondary metabolism showed a consistent non-linear variation across dry and wet seasons (Table 1 and Figure 2). The concentration of structural compounds (carbon, fibers, lignin and cellulose) and nutrients (N and P) were highest in wet and wet-dry transition seasons and decreased toward the dry season or until the beginning of the dry-wet transition. In contrast, secondary metabolites (polyphenols), and C:N and C:P mass ratios revealed the opposite pattern with highest values in dry and dry-wet seasons and a decrease toward the wet seasons. Our models also showed higher seasonal variance of nutrients than of carbon. 1 | Results of generalized additive mixed models (GAMMs) exploring non-linear temporal patterns of litter traits in terms of concentrations (%) or mass (g m −2 month −1 ) of each compound or C:N, C:P, and N:P ratios over 2 years. The amount of smoothing of models is presented as the effective degrees of freedom (edf) in relation to time (i.e., sampling events). We show estimated intercepts for each variable, values of F and P for the smoother and adjusted R 2 for the GAMMs. *Log-transformed variable.
When the amount of litter constituents in terms of mass that enter streams is taken into account a different pattern was found. The mass of all litter constituents (structural, nutrients and secondary metabolites) was two to threefold higher in drywet (when litterfall is greatest) than in wet or wet-dry seasons (Table 1 and Figure 3).

Congruence Between Litter Traits and Species Composition of Litterfall
Litterfall entering the three streams and across the four seasons was comprised of 142 plant species (85 species at CVE, 84 at CAP and 39 at RON stream). However, 64-74% of species (or 25-63 species) contributed less than 1% of annual litterfall biomass. As a result, 9-12 species dominated litter phenological patterns and were responsible for more than 70% of annual litter biomass (Figure 4 and Supplementary  Figure 1). Although some of the dominant species contributed at certain periods of the year, others were present year-round (e.g., Vochysia pyramidalis, Protium spruceanum, Calophyllum brasiliense, Xylopia emarginata) (Figure 4). Dominant species or genus that were found in the three riparian ecosystems showed similar litterfall patterns as supported by Calophyllum brasiliense and Myrcia (Figure 4).
The redundancy analysis (RDA) using data from all species showed that 46% of total variation (i.e., constrained variance) of species composition data was related to litter trait variation across streams and seasons ( Table 2). The RDA ordination evidenced that species with lower C:N mass ratios and lower cellulose were more frequent at CVE (e.g., Ormosia arborea and Myrcia tenuifolia), species with higher C but lower lignin at CAP (e.g., Vochysia pyramidalis) and species with higher C:N mass ratios, higher cellulose and higher lignin at RON (e.g., Xylopia emarginata and Clusia criuva) (Figure 5A). Patterns across seasons were apparent within streams. For example, the dry season was mostly related to higher C:N ratios and lower lignin, while wet season showed the opposite pattern with lower C:N ratios and higher lignin content. Also, the dry-wet season was more similar (i.e., closest in the RDA ordination) to the wet season and the wet-dry to the dry season ( Figure 5A).
A second RDA using species composition of only dominant species indicated that dominant species were driving the variability in litter traits across streams and seasons. This second RDA explained a higher percentage of total variation (46% in the first vs. 49% in this second RDA) while maintaining the same patterns of the first RDA ordination, supporting the predominant role of dominant species in litterfall chemistry ( Figure 5B and Table 2).

DISCUSSION
Our study showed clear seasonal patterns of litterfall chemistry, which were largely driven by dominant species in riparian forests of Brazilian Cerrado. Importantly, we found contrasting seasonal patterns in the concentrations and mass of litter traits. Despite the importance of understanding leaf phenological patterns in forests and riparian ecosystems to carbon and nutrient cycling, most available literature is still restricted to temperate and boreal regions (e.g., Niinemets and Tamm, 2005;Montgomery et al., 2020) or focused on specific sets of plant species (e.g., habits, functional type; Panchen et al., 2015). This is surprising considering that the tropics cover 40% of the Earth's land surface, harbor a high diversity of plants, and are responsible for processing a considerable portion of organic matter globally (Wright, 1996;Tiegs et al., 2019). Within this context, our findings contribute to the understanding of basic aspects of carbon and nutrient dynamics and help to disentangle how riparian diversity influence stream and riparian ecosystem functions.
The seasonal patterns of litterfall chemistry observed here could due to (i) species seasonal variation in litter chemistry (e.g., species with different litter chemistry in drier and wetter seasons because of different rates of nutrient resorption or leaf development), (ii) different species composition of litterfall (e.g., higher litterfall mass of a nutrient or lignin-rich species), (iii) intra-species litter chemistry variability (e.g., individuals of a species that produce litter with different chemistry) or (iv) the occurrence of all the above-mentioned mechanisms. Our finding that the species composition of litterfall is strongly coordinated with seasonal variation in litterfall chemistry indicates that interspecific variation in leaf senescence plays a determinant role in riparian ecosystem functions. Although intra-species variability in litter chemistry is expected (and was not tested in this study due to the lack of data), it has been shown to be lower than inter-species variability (Anderegg et al., 2018).
FIGURE 2 | Temporal patterns of litter traits as concentrations of each compound and, ratios of C:N, C:P, and N:P across the four seasons (dry-wet, wet, wet-dry, and dry seasons) over 2 years. Blue lines represent the non-linear temporal trend from generalized additive mixed models and dashed lines the 95% confidence intervals of the models. Circles are means and vertical lines denote upper and lower limits of standard errors from data of multiple sampling sites within streams.
Frontiers in Forests and Global Change | www.frontiersin.org FIGURE 3 | Temporal patterns of litter traits as the mass of each compound (g m -2 month -1 ) across the four seasons (dry-wet, wet, wet-dry, and dry seasons) over 2 years. Solid lines represent the non-linear temporal trend from generalized additive mixed models and dashed lines the 95% confidence intervals of the models. Circles are means and vertical lines denote upper and lower limits of standard errors from data of multiple sampling sites within streams.
Taken together, these results suggest that litterfall chemistry (on an annual or seasonal basis) is primary driven by the chemical composition of dominant species occurring in riparian forests.
Seasonally distinct litter chemistry of a species may reflect the influence of intrinsic (e.g., plant development) and extrinsic environmental controls (e.g., climate). In our study, nutrient resorption may have reduced nitrogen and phosphorus of litter from wetter to drier periods, when their concentrations were lowest and carbon to nutrient ratios were highest, as shown elsewhere (Anaya et al., 2007 in tropical dry forests;Heartsill-Scalley et al., 2012 in tropical rainforests). Nitrogen and phosphorus resorption plays a key role in nutrient conservation of plants (Brant and Chen, 2015) and its importance was previously reported for plant species of Cerrado biome (Kozovits et al., 2007). It is also noteworthy that the seasonal patterns of carbon is likely to be associated with those of N and P, as carbon assimilation by plants through photosynthesis is maximized only if sufficient nutrients are allocated to leaves under optimum irradiance (Hikosaka, 2005). A plausible explanation for the higher seasonal variance of nutrients besides carbon is that the former are limiting elements in most plants-and are reabsorbed before leaf senescence and recycled multiple time (Han et al., 2013;Brant and Chen, 2015),-while the turnover of the carbon is faster given its assimilation from the atmosphere (Chambers et al., 2004).
In contrast, seasonal changes in structural compounds may be associated with the growing season and leaf development of most (or dominant) plant species that coincide with first rains (i.e., the transition between dry and wet season) (Reich and Borchert, 1984;Reich, 1995). During leaf development, there is a continuous production of cell walls consisting of strong fiberrich composite materials such as lignin and cellulose, explaining their increase from drier to wetter seasons, which promote mechanical resistance to protect leaves from herbivores, wind  and other risks of injury (Boyero et al., 2017). Increases in polyphenols content on litter from wetter to drier seasons might be due a stimulation of polyphenols production by climatic stressors such as drought or increased temperature (Tharayil et al., 2011). An alternative explanation underlying decreases in polyphenols would be leaching of senescing leaves, before or even after falling in traps, during periods of heavy rainfall, as demonstrated previously (Serrano, 1992). However, it is unlikely that leaching was responsible for such striking reductions of polyphenols given that thick cuticles of most leaves may attenuate leaching losses. We show that, despite litter having lower concentrations of nitrogen and structural compounds in drier seasons, the amount of all litter constituents added to soils and streams was remarkably (2-3 times) higher in drier seasons as a consequence of higher litterfall mass. This suggests the existence of two contrasting patterns providing complementary information to the understanding of ecosystem functions. First, exploring the seasonality of the concentrations of each litter constituent allows for the prediction of rates for key ecosystem processes in both forests and streams, such as decomposition and nutrient cycling that depend on physical and biological agents (Marks, 2019). For example, lower in-stream litter decomposition reported in the dry season (e.g., Tonin et al., 2019) may be related to the higher C:N or C:P ratios observed in this study, and thus a lower availability of nutrients to decomposers and detritivores. Second, quantifying the mass of litter constituents provides a more holistic view of the ecosystem in terms of fluxes of elements and organic molecules from litterfall to the forest soils and stream. Thus, periods of large accumulation of elements and organic compounds such as nutrients or fiber-rich materials may drive the chemistry of soils and water (especially in low-flow conditions) stimulating or retarding the biological activity (Lutz et al., 2012). These contrasting patterns also indicate that there is a trade-off between lower amounts of higher-quality litter and higher amounts of lower-quality litter over the year. Such a trade-off may FIGURE 5 | Biplot ordinations of redundancy analyses (RDA) exploring the congruence between plant species composition of litterfall and litter chemical traits. Four non-colinear traits summarizing the chemical variation of litter were used in both RDA analyses (i.e., carbon, C:N, cellulose and lignin). Species data are shown with non-filled red circles. Species acronyms are shown only for the dominant species. Acronyms are the three initial letters of each species (or four letters if a genus represent more than one species) described in Figure 4. enable litter-consuming detritivores and microbial decomposers to meet their nutritional demands, suggesting that abundant low-quality litter may provide a long-term resource once the more nutritious and scarce litter is gone. However, experimental evidence suggests that compensatory feeding on poorer resources may be necessary to alleviate nutritional imbalances of organisms (Ott et al., 2012;Flores et al., 2014).
We provide novel evidence of seasonal patterns of leaf litterfall chemistry, mostly driven by the composition of dominant plant species in litterfall, in wet-dry tropical riparian ecosystems. We also show that seasonal patterns in litterfall chemistry revealed two contrasting perspectives, if presented in terms of concentrations or mass of compounds. Thus, we can no longer assume that tropical riparian forests and streams have less clear seasonal patterns than temperate or boreal ecosystems. While our study emphasizes the contributions of dominant species in determining the concentrations and amounts of C, N, and P in litter, it is possible that rare species contribute disproportionately to ecosystem functions or in the provisioning of essential elements (such as micronutrients) not considered in this study (cf. Mouillot et al., 2013). The applicability of our findings to other tropical riparian ecosystems should interpreted with caution given the limitations of the spatial scale of this study. Our results provide a basis for understanding how plant litter resources influence the biogeochemistry of riparian ecosystems and for predicting the consequences of potential future shifts of rainfall seasonality for ecosystem processes.

AUTHOR CONTRIBUTIONS
AT, JG, and RR conceived the study. PB and AT collected field data. LL and MF performed the chemical analyses. AT managed and analyzed the data. AT wrote the manuscript with feedback from JG and other authors. All authors contributed to the article and approved the submitted version.