The importance of allochthonous organic matter quality when investigating pulse disturbance events in freshwater lakes: a mesocosm experiment

Extreme precipitation is occurring with greater frequency and intensity as a result of climate change. Such events boost the transport of allochthonous organic matter (allo-OM) to freshwater ecosystems, yet little is known about the impacts on dissolved organic matter (DOM) quality and seston elemental stoichiometry, especially for lakes in warm climates. A mesocosm experiment located in a Turkish freshwater lake was designed to simulate a pulse event leading to increased inputs of allo-OM by examining the individual effects of increasing water colour (HuminFeed®, HF), the direct effects of the extra energetic inputs (alder tree leaf leachate, L), and the interactions of the single treatment effects (combination of both sources, HFL), along with a comparison with unmanipulated controls. Changes in the DOM quality and nutrient stoichiometry of the allo-OM treatment additions was examined over the course of the experiments. Results indicated that there was an increase of high recalcitrant DOM components in the HF treatment, in contrast to an increase in less aromatic microbially derived molecules for the L treatment. Unexpectedly, seston C:P ratios remained below a severe P-limiting threshold for plankton growth and showed the same temporal pattern in all mesocosms. In contrast, seston N:P ratios differed significantly between treatments, with the L treatment reducing P-limiting conditions, whilst the HF treatment increased them. The effects of the combined HFL treatment indicated an additive type of interaction and chlorophyll-a was highest in the HFL treatment. Our results demonstrate that accounting for the optical and stoichiometric properties of experimental allo-OM treatments is crucial to improve the capacity to explain extrapolated conclusions regarding the effects of climate driven flooding on freshwater ecosystems in response to global climate change.

elemental stoichiometry, especially for lakes in warm climates. A mesocosm experiment located in a Turkish freshwater lake was designed to simulate a pulse event leading to increased inputs of allo-OM by examining the individual effects of increasing water colour (HuminFeedÒ, HF), the direct effects of the extra energetic inputs (alder tree leaf leachate, L), and the interactions of the single treatment effects (combination of both sources, HFL), along with a comparison with unmanipulated controls. Changes in the DOM quality and nutrient stoichiometry of the allo-OM treatment additions was examined over the course of the experiments. Results indicated that there was an increase of high recalcitrant DOM components in the HF treatment, in contrast to an increase in less aromatic microbially derived molecules for the L treatment. Unexpectedly, seston C:P ratios remained below a severe P-limiting threshold for plankton growth and showed the same temporal pattern in all mesocosms. In contrast, seston N:P ratios differed significantly between treatments, with the L treatment reducing P-limiting conditions, whilst the HF treatment increased them. The effects of the combined HFL treatment indicated an additive type of interaction and chlorophyll-a was highest in the HFL treatment. Our results demonstrate that accounting for the optical and stoichiometric properties of experimental allo-OM treatments is crucial to improve the capacity to explain extrapolated conclusions regarding the effects of climate driven flooding on freshwater ecosystems in response to global climate change.

Introduction
Global climate change is altering precipitation patterns worldwide with reported widespread increases in both the frequency and intensity of daily rainfall (Toreti et al., 2013). Extreme precipitation events have become more intense in shorter time periods (Alpert et al., 2002). Consequently, there is a higher risk of flooding events but also more chances of summer droughts, which can directly affect soil moisture and stream flow (Cramer et al., 2018). Such conditions (i.e. warmer temperatures and intensified precipitation patterns) may boost the transport of terrestrial organic matter to freshwater ecosystems, varying in both quantity and quality, which are capable of producing cascade effects to aquatic food webs (Quante and Colijn, 2016).
Autochthonous organic matter (auto-OM) consists of carbon-based components directly derived from inlake primary producers, mainly comprising biopolymers of non-humic substances (Findlay & Sinsabaugh, 2003). Alternatively, allochthonous organic matter (allo-OM) is the fraction derived from the surrounding terrestrial ecosystems, typically originating from plant tissues, and predominantly based on structural cellulose and lignin-like compounds (Thurman, 1985). Allo-OM loading in freshwater ecosystems has direct and indirect effects due to an increase in both water colour (a phenomenon referred to as 'browning') and the extra energetic input to the base of trophic food webs supporting catabolic (i.e. respired) and anabolic (i.e. biosynthesis) metabolisms (Solomon et al., 2015). Browning of freshwater bodies reduces euphotic depth, which can strongly limit phytoplankton growth and affect community composition depending on nutrient availability . In contrast, dissolved allo-OM can enhance heterotrophic bacteria productivity, selectively allocating this terrestrial dissolved organic carbon subsidy to biosynthesis (Guillemette et al., 2016). Furthermore, particulate (seston) allo-OM can be directly or indirectly (through microbial loop components) available to higher trophic levels, increasing or maintaining primary consumers biomass (Cole et al., 2011;Tanentzap et al., 2017).
Studies investigating the consequences of increased allo-OM loading in lakes tend to only focus on the quantity of dissolved organic carbon (DOC) fraction (Solomon et al., 2015;de Wit et al., 2016). Dissolved organic matter (DOM) quality, however, also plays a key role (Kothawala et al., 2014;Stadler et al., 2020). Intense rainfall events can drive shifts towards lower quality DOM molecular composition in freshwater bodies by relatively increasing concentrations of terrestrially derived high recalcitrant molecules (Zhou et al., 2020(Zhou et al., , 2021b, directly impacting microbe-mediated transformation of DOM and consequently altering energy flow through the base of aquatic ecosystems (Zhou et al., 2021a). The optical properties of DOM are related to its biochemical characteristics and, therefore, supply information on its origin and quality for planktonic communities (McKnight et al., 2001). These properties distinguish between coloured dissolved organic matter (CDOM), organic compounds that absorb light (chromophores), and fluorescent dissolved organic matter (FDOM), chromophores that re-emit light (fluorophores). Despite the complexity of DOM, specific patterns of FDOM have been associated with ecologically meaningful characteristics (Jaffé et al., 2014). Thus, FDOM constituents can range from humic-like fluorescence peaks strongly associated with high aromaticity, colour and highly persistent (i.e. difficult to degrade) poor-quality molecules mainly derived from terrestrial sources; to protein-like fluorescence peaks more associated with non-coloured, non-aromatic, and very labile (i.e. easy to degrade) high quality molecules mainly originating from in-lake primary producers (Fellman et al., 2010). For these reasons, efforts at understanding the effects of allo-OM loads to aquatic ecosystems would benefit from considering DOM quality in addition to its quantity.
Increases in extreme precipitation events leading to higher river inflows also mobilise particulate carbon (C) and export macro-nutrients such as nitrogen (N) and phosphorus (P) to lakes (Sebestyen et al., 2008;Dhillon & Inamdar, 2013). Factors controlling the type of allo-OM entering freshwater ecosystems are broad and catchment dependent (Qualls & Richardson, 2003). Phosphorus and N mobilisation mechanisms are complex, with a range of molecular weight distributions and consisting of both inorganic and organic fractions (Asam, 2012;Ged & Boyer, 2013;Kelly et al., 2019). Following high-intensity storms, particulate organic carbon (POC) export can be 6 to 8 times higher than DOC, dramatically altering the proportion of DOC and POC inputs to aquatic ecosystems (Dhillon & Inamdar, 2013;Hitchcock & Mitrovic, 2015). In addition, organic C exported from agricultural catchments after stormflows can produce a significant load (i.e. 50%) of relatively photoreactive and aromatic DOC into freshwater bodies (Caverly et al., 2013). In general, observations clearly indicate that POC and DOC can be highly mobilised during extreme rainfall events, with differences between relative contributions of dissolved and particulate fractions depending on the catchment type and land use (Hope et al., 1994;Jeong et al., 2012). In turn, N and P exports are also often higher during peak streamflow (Hitchcock & Mitrovic, 2015;Kelly et al., 2019). Indeed, dissolved nutrient subsides can have stimulatory (higher nutrient availability) and inhibitory (higher nutrient limitation) effects for primary producers and bacteria (Hitchcock & Mitrovic, 2015). In addition, examining elemental stoichiometry of seston is crucial for interpreting the effects of resource quality on zooplankton primary consumers (Sterner & Elser, 2002).
To date, laboratory experiments have been designed to better understand the complex effects of allo-OM inputs on planktonic communities (Brett et al., 2009;Mcmeans et al., 2015;Vad et al., 2020). Furthermore, a commercially available humic substance called HuminFeedÒ (Humintech, Germany) has been widely used in mesocosm experiments investigating allo-OM inputs to lakes (Rasconi et al., 2015;Urrutia-Cordero et al., 2017;Lebret et al., 2018;Minguez et al., 2020). Considering the complexity of catchment-dependent allo-OM characteristics, other mesocosm studies have mimicked allo-OM using a range of different types of natural mixtures such as i) humic soil and inorganic N and P additions (Lefébure et al., 2013), ii) leaf leachate and inorganic N and P additions (Hitchcock et al., 2016;Brighenti et al., 2018;Tonetta et al., 2018), iii) filtered and concentrated humic lake water (Nicolle et al., 2012;Cooke et al., 2015;Hamdan et al., 2021), and iv) organic matter-rich top soils from inflowing rivers produced by preserving natural fine soil particles (Liess et al., 2015). In addition, in order to separate the browning effect from allo-OM as energy source effects, some authors simulated light reduction by adding a black net to reduce solar irradiance (Urabe et al., 2002;Brighenti et al., 2018).
Overall, experimental studies have shown that despite a general pattern of increasing heterotrophic activity following allo-OM disturbances in aquatic ecosystems, there are still a wide range of possible responses from primary producers and consumers. However, characterisation of DOM quality (Lennon et al., 2013) and changes in seston nutrient stoichiometry (Liess et al., 2015;Minguez et al., 2020) have received little attention, despite their potential importance in explaining the divergent findings regarding the effects of allo-OM on aquatic ecosystems (Creed et al., 2018;Senar et al., 2021). None of the published mesocosm experiments, to date, have attempted to characterise how quality features of different allo-OM sources would impact freshwater planktonic communities (Lennon & Cottingham, 2008;Geddes, 2015;Hitchcock et al., 2016;Lebret et al., 2018;Fonseca et al., 2022), despite the known variability of allo-OM sources in nature (Kellerman et al., 2014). Finally, studies of allo-OM pulses following extreme precipitation are scarce for lakes in warmer climates compared to those carried out in north temperate and boreal aquatic systems (Monteith et al., 2007;de Wit et al., 2016).
Experiments which address past experimental limitations are essential for gaining greater consensus on allo-OM effects to the base of aquatic pelagic food webs. Therefore, a mesocosm experiment located in a Turkish freshwater lake was carried out to explore the system responses to allo-OM pulse disturbance events. The experimental treatments were designed to examine the individual effects of increasing water colour (HuminFeedÒ, ''HF''), the extra energetic inputs (alder tree leaf leachate, ''L''), and the interactions of the single treatment effects (combination of both sources, ''HFL''), along with a comparison with unmanipulated controls. Specifically, the work presented here aimed to first characterise the quality of different allo-OM source additions during the experiment by examining DOM fluorescence indicators and seston food quality for primary consumers (i.e. seston C:P and N:P ratios); to compare the overall treatment effects on water quality parameters with the log response ratio (LRR) detailed in Hillebrand & Gurevitch (2016); and subsequently, investigate the response in primary producers by examining chlorophyll-a (Chl-a) trends over the course of the experiment. We expected the following: i. The HF treatment would increase water colour and provide a low quantity of high recalcitrant DOM , compared to the L treatment that would add higher quantities of more labile DOM, without significantly increasing water colour (Hitchcock et al., 2016), relative to controls. ii. The HF treatment would directly increase seston C:P ratios, without affecting seston N:P composition (Minguez et al., 2020), whilst the L treatment would add both seston C and P (Navarro et al., 2019), increasing seston C:P ratios, but less than in the HF treatment, and reduce seston N:P, relative to controls. iii. Synergistic interactions would prevail in the combined HFL treatment in that DOM fluorescence indicators and seston elemental ratios LRR effects would be significantly greater than the addition of the respective single treatments due to biotic and abiotic processes occurring when both sources were mixed. iv. Chl-a would, relative to the controls, decrease in the HF treatment due to light and nutrient constraints (Lefébure et al., 2013;Rasconi et al., 2015), and increase in the L treatment due to nutrient additions arising from the leaf leachate and low light limitation compared to the HF treatment (Hitchcock et al., 2016;Brighenti et al., 2018). In the combined HFL treatment we expected Chl-a to remain higher than the HF but lower than the L treatments due to single treatment counteracting effects (i.e. increase of light limitation but addition of nutrients).

Study area and experimental design
At the beginning of June 2019, 16 cylindrical fibreglass enclosures (1.2 m diameter 9 2.2 m height 9 4 mm thick) were fixed on a floating platform called Middle East Technical University (METU) Mesocosm System (39°52 0 13.18 00 N, 32°46 0 31.92 00 E), in Ankara, Turkey. The experimental set-up consisted of a control and three different treatments in a 4 9 4 factorial design (i.e. 4 treatments and 4 replicates per treatment). Water from the experimental lake, sieved through a 500 lm mesh, was used to fill the enclosures (2480 L initial volume) and ensure a natural chemistry of the water (Landkildehus et al., 2014). A water pump (3 W, Sobo wp-50M) was placed in the middle of each enclosure (150 cm below surface water) to ensure sufficient mixing of the whole water volume (i.e. avoiding water column stratification and sedimentation processes). A bird protecting net (3 9 3 cm mesh size) was placed on top to prevent the entrance of any bigger entity. The mesocosms did not include sediments, macrophytes, or fish to target effects on pelagic planktonic ecosystems. Mesocosms were inoculated on the 11th of June with a mixed sample of phytoplankton and zooplankton communities from five other nearby lakes whose trophic states varied from oligotrophic to eutrophic. To collect inoculated phytoplankton samples, 10 l of whole water column from the five other nearby lakes using a Ruttner sampler were sieved through 55 lm mesh to remove zooplankton and prevent grazing of phytoplankton before inoculation of mesocosms. For inoculated zooplankton samples, five vertical hauls covering the entire water column and single horizontal haul over the longest horizontal transect of the lake were taken by using a 55 lm mesh plankton net. Concentrated zooplankton samples were stored in 10 l of surface lake water before inoculation of mesocosms. Until the inoculation day (maximum 3 days after collection), the plankton samples were oxygenated (via gentle bubbling) and kept in a temperature controlled (22 ± 1°C) room under low, indirect light. For the inoculation, planktonic samples were pooled together and well mixed. Two litres of the planktonic sample were added into each mesocosms. After inoculation, the mesocosms were incubated for a period of nine days before the pulse disturbance event to allow potential development of a diverse planktonic community with similar starting conditions. The additions of the three different allo-OM sources were carried out on the 20th of June (day 0). After the pulse disturbance the experiment ran until day 36, the 26th of July.

Allo-OM treatments
HuminFeedÒ is a commercially available lignite derived product and its detailed chemical composition has been fully characterised by Meinelt et al. (2007). In our experiment, the HuminFeedÒ stock solution was freshly prepared a few hours before the pulse by manually mixing for approximately 30 min at room temperature 99.2 g of HuminFeedÒ powder (stored in the fridge at 4°C) in 40 l of distilled water (at room temperature). Final HuminFeedÒ stock solution concentration was 2480 mgHF l -1 . To simulate local terrestrial carbon sources, approximately 12 kg of fresh alder leaves were collected and spread out evenly in a ventilated greenhouse to allow them to dry quickly for 48 h five days prior to the mesocosms additions. After the drying period, leaves were placed in a dark cold room (4°C) with distilled water (30 g of dried leaves per each litre of distilled water added) for 48 h. Darkness and cold environment was chosen following methodology in previous studies (Hitchcock et al., 2016;Fonseca et al., 2022) to reduce biotic and abiotic degradation of the prepared leaf leachate. After sieving the leaves from the water through 45 lm mesh size to exclude larger leaf material, a final stock solution was produced with a final DOC concentration of 2506 mgC l -1 . Alder tree leaf leachate was chosen because it had been used in similar studies, could be found in a nearby forested area, and were logistically viable to obtain enough DOC concentration of the stock solution in a relative short period of time (Mutschlecner et al., 2018). Therefore, both stock solutions were added to the respective mesocosms immediately after preparation without any further sterilisation (i.e. autoclave) either filtration (i.e. through 0.22-0.45 lm) processes, such as carried out in other studies (Nicolle et al., 2012;Geddes, 2015;Brighenti et al., 2018). The experimental setup consisted of four replicates for each of the following treatments applying concentrations used in other published work: (1) unmanipulated controls (Cntl); (2) the HF treatment which received 4 l of Humin-FeedÒ stock solution resulting in a final HuminFeedÒ addition of 4 mgHF l -1 (equivalent to * 1.5 mgC l -1 as DOC) to increase water colour (Urrutia-Cordero et al., 2017;Lebret et al., 2018); (3) the L treatment by adding 7.9 l of alder tree leaf extract stock solution resulting in a corresponding added DOC concentration of * 8 mgC l -1 to test for an extra energetic input (Hitchcock et al., 2016); and (4) the combined HFL treatment composed of the same volume of stock solutions as added in the individual treatments, resulting in a pulse addition of 4 l of HuminFeedÒ stock solution (addition of * 1.5 mgC l -1 as DOC) and 7.9 l of alder tree leaf extract stock solution (addition of * 8 mgC l -1 as DOC) to compare the interactions of the single treatment effects (Fonseca et al., 2022). Owing to the focus on examining qualitative characteristics of the allo-OM, DOC quantities in this study are not directly comparable.

Mesocosm sampling strategy
Dissolved, seston, total C, N, and P, and Chl-a were measured every 4 days. In addition, parameters were sampled before (i.e. day 0) and after (i.e. day 1) the pulse treatment additions. The enclosures were sampled from the surface to the bottom using integrated tube samplers (one tube per treatment to avoid contamination). For each enclosure, three depthintegrated samples were taken across the diameter of each mesocosm and were pooled in 15 l plastic buckets. A subsample of the mixed water was transferred to acid washed 500 ml plastic bottles for water chemistry analyses. Three litres of 45 lm presieved water was transferred to pre-washed plastic bottles to collect the seston fraction. Collected samples were stored in dark in a cooler until further analyses or preservation step in the laboratory (within 3-4 h).

Laboratory analyses
Dissolved samples were obtained by filtering water through Whatman GF/C glass fibre filters. Dissolved organic carbon concentrations were analysed with a Shimadzu TOC-L/CPN analyser (Schimadzu Scientific Instruments, Japan). Total nitrogen (TN; unfiltered), ammonium (NH4 ? ; filtered), and nitrite-nitrate (NO 2,3 filtered) were determined using an automated wet chemistry analyser (San??, Skalar Analytical, The Netherlands) (Baird & Bridgewater, 2017). Ammonium, nitrite, and nitrate concentrations were summed and reported as dissolved inorganic nitrogen (DIN). Total phosphorus (TP; unfiltered) and soluble reactive phosphorous (SRP; filtered) concentrations were determined spectrophotometrically using the molybdenum blue method (Mackerth, 1978). Between 200 and 600 ml of pre-sieved water through a 45 lm mesh were filtered onto pre-combusted and preweighted Whatman GF/F glass fibre filters for seston C, N, and P content analyses. In total 2 filters per enclosure were collected (one for C and N analyses and the other for P analyses). Filters were immediately dried at 60°C for 24 h and then weighted for dry weight calculations. Filters were preserved in dark at room temperature inside zip lock plastic bags until further analyses. Prior to C and N analyses, seston filters were left for 6 h in a desiccator containing a beaker of hydrochloric acid (HCL, 37%) to remove inorganic C (Harris et al., 2001). Filters were then measured for C and N content using an elemental analyser (vario EL cube, Elementar, Germany). Standard reagent used for this analysis was Acetanilide (C 8 H 9 NO) with a composition of 10.36% of N and 71.09% of C. Seston P content was determined as specified for TP following combustion at 550°C for 2 h and digestion step with potassium persulfate under high pressure and temperature (Autoclave for 30 min at 15 psi and 121°C). Three blank filters per batch were considered for all C, N, and P analyses (McCarthy et al., 2006). Chlorophyll-a samples were filtered onto Whatman GF/C glass fibre filters and then extracted with ethanol before reading the absorbances in a UV-Vis spectrophotometer (Perkin Elmer Lamb-da35) (limit of detection [LoD] of 0.04 lgChl-a l -1 ) (Jespersen & Christoffersen, 1987).

Optical properties of DOM
For logistical reasons, only water samples from day 4 to day 24 were filtered through Whatman 0.45 lm filters to carry out full UV-Vis spectrophotometer scans and fluorescence excitation-emission matrices (EEMs). DOM absorbance was measured using an UV-Vis spectrophotometer (UV-1800, Shimadzu, Germany) carrying out full scans across wavelengths from 200 to 900 nm using 1 cm glass cuvette and Milli-Q water as blank. Absorbance units measured at 420 nm were used to characterise water colour as a proxy of brownification . Fluorescence EEMs were characterised using a fluorescence spectrophotometer (Varian Cary Eclipse, Agilent technologies, USA) in 1 cm quartz cuvette. Fluorescence intensities were measured at excitation wavelengths ranging from 220 to 450 nm (5 nm increments) and emission wavelengths ranging from 240 to 600 nm (2 nm increments). Fluorescence EEMs were analysed using the drEEM toolbox (Murphy et al., 2013) in MATLAB version R2020b (The Math Works, Inc., 2020) and available code can be found in supplementary material. EEMs were corrected for instrument-specific biases using manufacturer supplied instrument correction factors. Spectra were also blank subtracted and corresponding absorbance measurements were used to correct for inner filter effects. Arbitrary fluorescence intensities were then normalised to a standard scale of Raman Units (RU). Any residuals of Raman and Rayleigh non-trilinear scatter were removed, without further interpolating the missing values, before multivariate parallel factor analysis (PARAFAC) (Murphy et al., 2013). Specific fluorescence components were identified by carrying out PARAFAC models. These models assume that each EEM sample set can be decomposed depending on a predefined number of components (Murphy et al., 2013). Steps to design the optimal PARAFAC model include the following: (1) the identification and removal of outliers, (2) development of preliminary models to explore the appropriate number of components, (3) determination of the correct number of components by finding chemically meaningful spectral loadings, and (4) internal model validation using split-half analysis by comparing models with different subsets of the original data (Murphy et al., 2013). Accordingly, the best model identified four fluorescence components, labelled C1-C4. Components identified ranged from high recalcitrant humic-like substances of terrestrial origin to labile protein-like freshly produced components, C3, C1, C4, and C2, respectively (Table 1). Finally, the optimum validated PARAFAC model was exported to OpenFluor, an online repository of published organic fluorescence spectra, to link the mesocosm fluorescence DOM patterns with trends observed in other aquatic systems (Murphy et al., 2013). PARAFAC model will be publicly available in OpenFluor once it can be linked with a published study.

Statistical methods
Pairwise comparisons of the four different treatments (four-factor ANOVA), with the Bonferroni method to adjust P-values by the number of comparisons (Heckler, 2005), were applied to different nutrient fraction variables every sampling day of the experiment (days 0, 1, 4, 8, 12, 16, 20, 24, 28, 32, and 36) to determine statistical differences amongst treatments. Pairwise t-tests were performed using the t_test function in the rstatix package in R (Kassambara, 2020).
A non-metric multidimensional analysis (NMDS) followed by a pairwise adonis test was used to describe and compare changes in the PARAFAC components through the mesocosm experiment (from day 4 to day 24). A NMDS is a robust ordination technique based on a distance matrix by representing the pairwise (dis)similarity between samples in a low-dimensional space (Clarke, 1993). The Bray Curtis dissimilarity matrix was applied to relative intensity means of replicates of PARAFAC components which were fourth root transformed to down weigh the effect of the most dominant component. Samples that are ordinated closer to one another are likely to be more similar than those further apart. However, the scale of the axes of a NMDS is arbitrary as is the orientation of the plot. In addition, solutions with high stress values (usually [ 0.2) should be interpreted with caution, whilst stress values \ 0.05 indicate that the ordination is an excellent representation of the high-dimensional assemblage structure and there is low probability of misinterpretation (Buttigieg & Ramette, 2014). A pairwise adonis post-analysis was run to test for statistical differences between treatments (Anderson, 2017). The NMDS was computed with the function metaMDS and the pairwise adonis test with the function pairwise.adonis (Martinez Arbizu, 2017) of the vegan package in R (Oksanen, 2011).
We fitted generalised additive models (GAMs) with a Gaussian distribution (Wood, 2017) to seston nutrient ratios and Chl-a concentrations (all response variables were ln-transformed as recommended when dealing with nutrient ratios (Isles, 2020) and to avoid removal of Chl-a outliers). We analysed the effect of the treatments (factor variable of 4 levels) as well as temporal trends (i.e. days of the experiment: DOE) considering ordered-factor-smooth interaction and using the control treatment as the reference factor (Rose et al., 2012). After checking for the violation of the assumptions, if violated, we added a simple temporal correlation structure with the days of the experiment (correlation = corAR1(form = * 1|DOE)) to account for remaining temporal dependence of the observations (i.e. generalised additive mixed models, GAMMs). All GAM/Ms were carried out considering restricted maximum likelihood method (REML) as the parameter estimation procedure and optimum cubic regression splines (bs = 'cr') for the explanatory variables (Zuur et al., 2009). Specific mesocosm replicates' variability was checked by adding individual enclosures (n = 16) as a model smoother with random effect splines (bs = 're') (Pedersen et al., 2019), but these fixed random effects were dropped from all models as enclosures variability effects were insignificant. Models were compared by means of the Akaike Information Criterion (AIC), with lower AIC values indicating better models (Wood, 2017) (Table S2, S3, and S4). GAMs and GAMMs were performed using the gam and gamm functions, respectively, in the mgcv package in R (Wood & Wood, 2015).
A principal component analysis (PCA) was run to explore differences in mesocosm treatments based on colour, PARAFAC components, dissolved and seston nutrient concentrations and ratios as descriptors (15 variables), including sampling points from day 4 to day 24 and all treatment replicates (n = 96). The PCA summarises the general pattern between treatment observations and the available descriptive variables in a low-dimensional space illustrated as principal components and representing a high percentage variability of the data (Buttigieg & Ramette, 2014). The PCA was computed with the function prcomp in base R (R Core Team., 2019).
The magnitude of the effect of the addition of different allo-OM sources on carbon and nutrient quantity and quality parameters was calculated at each sampling time using the log response ratio (Wyatt & Rober, 2020). The LRR quantifies the proportionate change resulting from experimental manipulation and is calculated as the natural log (ln) of the treatment mean divided by the control mean (Hillebrand & Gurevitch, 2016). The direction of the effect size indicates whether the treatment had a positive or negative effect on the correspondent tested variable. In this study, calculated effect sizes at different days of the experiment have been grouped by treatments in boxplots of between 8 and 10 observations (n), depending on the variable studied and based on the sampling days, from day 1 to day 36. In turn, the respective LRR measures at day 0 (before allo-OM additions) were used as benchmarks and when day 0 was not available, the benchmark value was 0 (Hillebrand et al., 2018). With these effect sizes we wanted to test if, overall, (1) the treatment effects were significantly different from benchmark of day 0 and (2) the treatment effects were significantly different between treatments.
All R analyses were performed in version 4.0.3 (R Core Team., 2019) and available code can be found in supplementary material.

Allo-OM additions
The initial (day 0, before treatment addition) DOC concentrations in all enclosures were the same, ranging from 3.4 to 2.2 mgC l -1 . Immediately following the allo-OM pulse (day 1), the mean (± SD) DOC concentration significantly increased by 6.2 ± 1.4 mgC l -1 in the L (P-value \ 0.05) and 6.6 ± 0.2 mgC l -1 in the HFL (P-value \ 0.01) treatments, contrasting with insignificant changes in the HF treatment, relative to the control mesocosms ( Fig. 1a; Table S1). The control treatment had a maximum DOC concentration of 5.4 ± 1.5 mgC l -1 at day 20 whilst HF enclosures had maximum values of 6.3 ± 1.4 mgC l -1 at day 28. In contrast, the DOC concentrations in the L and HFL treatments reached minimum values of 7.0 ± 0.1 mgC l -1 at day 36 and 6.3 ± 0.6 mgC l -1 at day 32, respectively (Fig. 1a).
The initial concentrations of seston C at day 0 were not significantly different across the treatments and varied between 0.6 and 0.9 mgC l -1 . The allo-OM additions then resulted in a significant increase of 0.3 ± 0.0 mgC l -1 , 0.4 ± 0.1 mgC l -1 , and 1.4 ± 0.2 mgC l -1 , respectively for the HF, L, and HFL treatments (P-values \ 0.01; 0.01; 0.001) ( Fig. 1b; Table S1). The seston C concentrations in the control mesocosms gradually decreased from 0.7 ± 0.1 to 0.3 ± 0.1 mgC l -1 during the experiment. In the HF, L, and HFL treatments, seston C also gradually decreased over the experiment to concentrations ranging from 0.3 to 0.5 mgC l -1 at day 36. None of the treatments were significantly different from the controls at the end of the experiment (Fig. 1b; Table S1).
HuminFeedÒ added significantly less water colour than the leaf leachate source (P-values \ 0.01; 0.001 for HF and L, respectively), and colour in the HF treatment changed little throughout the experiment having an overall value of 0.007 ± 0.002 AU cm -1 (Fig. 1c; Table S1). In contrast, both the L and HFL treatments had the highest colour values immediately after the addition (Fig. 1c). Water colour of the L and HFL treatments decreased as the experiment progressed, from 0.018 ± 0.001 and 0.015 ± 0.001 AU cm -1 to 0.007 ± 0.003 and 0.013 ± 0.001 AU cm -1 , respectively, from day 4 to day 24. Water colour in the control treatment was significantly lower than all treatments throughout the experiment having an overall stable value of 0.001 ± 0.001 AU cm -1 ( Fig. 1c; Table S1).
Initial DIN concentrations (day 0, before treatment addition) in all enclosures were highly variable and ranged from 15 to 117 lgN l -1 . Overall, the DIN pool did not increase after the allo-OM pulse disturbance (Fig. 1d). The HFL treatment had significant lower DIN concentrations at day 8 (P-value \ 0.05), whilst the HF treatment had significant higher values at day 12 (P-value \ 0.05), respective to controls ( Fig. 1d; Table S1). For all treatments, the highest mean DIN values were observed between days 24 and 28. Final DIN concentrations in all mesocosms ranged from 37 to 83 lgN l -1 , which was within day 0 range (Fig. 1d). The initial seston N content of all replicates ranged from 325 to 649 lgN l -1 . The control treatment had a seston N minimum concentration of 219 ± 125 lgN l -1 at day 8 of the experiment, increasing again to the initial concentrations of 497 ± 65 lgN l -1 at the end of the experiment. After the allo-OM additions, the HF and L treatments experienced a similar increase in N concentration of 248 ± 208 and 356 ± 154 lgN l -1 , respectively. In the HFL treatment, seston N concentration increased to 676 ± 362 lgN l -1 . Only the L treatment increase was significant (P-value \ 0.05) because the HF and HFL treatments experienced high variability amongst replicates ( Fig. 1e; Table S1). At day 12, the seston N concentrations in the different treatments were already reduced to pre-disturbance conditions and remained stable until the end of the experiment (Fig. 1e). Total nitrogen values were highly variable amongst treatments. There was an overall increase in TN values in all mesocosms as the experiment progressed, with concentration ranging between 92 and 693 lgN l -1 from day 0 to day 36 ( Fig. 1f; Table S1).
Unlike N, the HF treatment did not add more P to the system. Therefore, P fractions in the control and HF treatments remained constantly low and were not significantly different throughout the experiment (Fig. 1g, h, i; Table S1). The overall average SRP concentrations was 3.6 ± 3.2 lgP l -1 whilst the seston P content was 6.1 ± 4.4 lgP l -1 . In contrast, leaf leachate additions added P in both the L and HFL treatments, and increased the SRP and seston P concentrations of these two treatments by 18.4 ± 8.2 and 42.2 ± 11.7 lgP l -1 , respectively, from day 0 to day 1 (Fig. 1g, h; Table S1). Overall concentrations of SRP in the HFL treatment after the pulse remained high throughout the 36 days of the experiment (i.e. 19.4 ± 7.4 lgP l -1 ), whilst SRP in the L treatment gradually decreased from 25.0 ± 5.5 lgP l -1 at day 1 to pre-pulse conditions (i.e. 2.8 ± 1.0 lgP l -1 ) from day 20 onwards (Fig. 1g). The concentrations of seston P in the L and HFL treatments followed the same pattern, decreasing from maximum values at day 1 (i.e. 51.9 ± 10.5 lgP l -1 ) to 7.0 ± 0.8 lgP l -1 at day 36 (Fig. 1h). In turn, TP values clearly replicated the pattern of the sum of the dissolved (SRP) and particulate (seston) fractions. Overall, the added P pool seemed to be consumed as the experiment progressed (Fig. 1i).
DOM quality PARAFAC analysis of the fluorescent components of the DOM pool highlighted the dominance of more labile molecules (C2) in control mesocosms compared with the allo-OM treatments (HF, L, and HFL) (Fig. 2). The HF treatment contained a high proportion of more recalcitrant humic-like molecules (C3 and C1). In contrast, the L treatment had a more even distribution of each of the four components, with a slightly greater prevalence of the less recalcitrant microbial derived humic-like C4 component (Fig. 2). In turn, the HFL treatment showed evidence indicating the presences of molecules from both the HF and L sources with an overall increase in humic-like components (C3, C1, and C4) (Fig. 2).
PARAFAC component composition revealed significant differences amongst treatments (pairwise adonis test, P-values \ 0.05) (Fig. 3 and Fig. S1). The more labile protein-like C2 component found in all treatments (Fig. 2) dominated in the control enclosures during the whole study period (Fig. 3 and  Fig. S1). The HFL treatment clustered with the C3 component whilst the HF treatment clustered closer to the C1 component, indicating a higher overall proportion of terrestrial humic-like and high recalcitrant poor-quality DOM in the HFL treatment (Fig. 3). In contrast, the L treatment grouped closer to the C4 component, which correspond to microbial humiclike, more labile DOM (Fig. 3 and Fig. S1). Temporal patterns from different treatments showed that, towards the end of the experiment, relative proportions of DOM components barely changed for each respective treatment obtaining different DOE observations grouped together to specific different components ( Fig. 3 and Fig. S1).

Seston quality
Before the allo-OM pulse (day 0), all mesocosms experienced moderate P limitation based on published Component EEM peaks and PARAFAC contours of emission and excitation matrices are given seston C:P nutrient limitation thresholds for zooplankton primary consumers (ln-transformed seston C:P molar ratios [ 4.86) (Hecky et al., 1993) (Fig. 4). Using GAMM models to assess temporal trends, all treatments were found to have the same significant temporal pattern in seston C:P (P-value \ 0.01), with decreases occurring directly after the pulse additions, but gradually returning to pre-disturbance conditions during the experiment ( Fig. 4; Table S2). In turn, the overall values in HF did not differ from controls, whilst being significantly lower in the L and HFL treatments (P-values \ 0.001). At the end of the experiment, seston C:P in the controls and HF treatments approached the severe P-limitation threshold (ln-transformed seston C:P molar ratios [ 5.56) for zooplankton herbivores (Urabe et al., 2002) ( Fig. 4). The model explained 55.9% of the deviance in seston C:P ratios over the experimental days (Table S2). Despite the allo-OM additions, seston N:P values of all enclosures remained above the theoretical severe seston P-limitation threshold, even though the treatments receiving the leaf leachate (i.e. L and HFL) had significant lower values of N:P relative to the control (P-values \ 0.001), and the HF treatment significantly higher (P-value \ 0.05) ( Fig. 5; Table S3). In contrast to the seston C:P temporal dynamics, the trends in seston N:P for the L and HFL treatments were significantly different from the control (Fig. 5; Table S3). Seston N:P ratios gradually increased in the control and HF treatments during the experiment, specially between day 12 and 20. In the L treatment, seston N:P rapidly decreased after the pulse, but quickly increased to initial conditions by day 12 and remained high until the end of the experiment (Pvalue \ 0.01). In contrast, the HFL treatment changes were weaker, with dynamics similar to the control and HF treatments but with absolute values resembling the L mesocosms (P-value \ 0.05) (Fig. 5). Overall, this model explained 68.9% of the deviance in seston N:P ratios during the mesocosm experiment (Table S3).

Overall treatment effects
There was a clear separation between treatments with regards carbon and nutrient quality and quantity parameters (Fig. 6). The separation by the two principal components, PC1 with a weight of 49.7% and PC2 with a weight of 14.8%, was mediated mainly by a differentiation in the DOM quality in both axis and by the additions of C and P fractions highly represented by PC1. Therefore, the control treatment observations clustered together and were characterised by relatively high concentrations of N compared to low P concentrations and a more labile DOM pool (C2). In turn, the HF treatment seemed to be differentiated by high seston C:P and N:P ratios and the recalcitrant fulvic-like DOM pool (C1) and, therefore, aggregated separately from the control samples. In contrast, the L treatment samples appeared aggregated in the opposite direction to the HF treatment and strongly associated with high colour, DOC and P concentrations. Finally, the HFL treatment observations grouped in between HF and L treatments and in the opposite direction to the control cluster, directly associated with high seston C, N, and P content and the high recalcitrant humic-like C3 DOM component (Fig. 6).
For all the included variables, which were those relating to carbon quantity as well as those related to nutrient seston quality and quantity, the HF treatment always had the lowest effect compared to the L and HFL treatments (Fig. 7a, b, d), with the exception of DOM quality high recalcitrant humic-like C3 and fulvic-like C1 components (Fig. 7c). In turn, the LRR effects of the HFL treatment indicated an additive type interaction (i.e. the sum of the effects of the treatments involved) with this treatment having more similar overall effects with the L than the HF treatments (i.e. LRR were not significantly different from the L treatment) (Fig. 7a, b, c, d). The LRR effect sizes of the seston elemental ratios showed that the HF and the treatments receiving the leaf leachate source evolved in opposite directions over the course of the experiment (i.e. positive vs negative LRR) (Fig. 7d).

Response of Chl-a
Chlorophyll-a temporal dynamics in the different treatments did not significantly differ from the control treatment. However, the L and HFL treatments had significantly higher Chl-a concentrations relative to the control (P-values \ 0.05; 0.001, respectively). The GAM model explained 28.9% of the Chla deviance (Fig. 8; Table S7). Chlorophyll-a in all the enclosures during the experiment ranged from \ LoD to 5.6 lgChl-a l -1 (Fig. 8). Before the pulse (day 0), overall Chl-a concentrations were low with mean ± SD of 1.8 ± 1.0 lgChl-a l -1 . At day 12 all mesocosms reached minimum values \ LoD. From day 16 to the end of the experiment Chl-a in the control and HF treatments remained low (i.e. 0.4 ± 0.5 lgChl-a l -1 ). In contrast, from day 16 onwards, the treatments receiving the leaf leachate source (L and HFL) reached higher Chl-a concentrations of 1.0 ± 1.0 lgChl-a l -1 . The highest value of 5.6 lgChl-a l -1 in the HFL treatment was recorded at day 32, in only one replicate for this treatment (Fig. 8).

Discussion
As concerns relating to ongoing changes in allo-OM export to lakes and reservoirs have increased, there has been a parallel increase in the use of mesocosm studies to investigate these effects. Here we characterised the quality of different allo-OM additions to a mesocosm experiment carried out in a warm climate freshwater lake. To our knowledge, this is the first study to examine changes in both the constituent DOM components and seston nutrient ratios of different allo-OM source additions over the course of a mesocosm experiment. We found that the different allo-OM sources resulted in significantly different treatment LRR effect sizes of the evaluated water quality aspects. These distinct treatment effect sizes should be key for understanding associated responses in the aquatic food web and elucidate divergent responses to these inputs within the system.
Importance of DOM quality along quantity As expected, the addition of HuminFeedÒ increased the proportion of high recalcitrant DOM in the system, mainly the terrestrial fulvic-like C1 component and to a lesser degree the terrestrial humic-like C3  Table S1. P-values significant codes differing from controls (Cntl): P \ 0.001 = '***'; P \ 0.01 = '**'; P \ 0.05 = '*'; non-significant (Table S1). The treatment with significant difference is indicated by the respective colour (grey HF; yellow L; red HFL) component. The humic-like C3 component has also been found in agricultural-forested catchments (Derrien et al., 2018), but in very low amounts compared to other components. On the other hand, the fulvic-like C1 component has been recurrently found in soil organic matter (Chai et al., 2019) and in a wide range of freshwater lakes and streams (Kida et al., 2019), being particularly dominant in forested and agricultural catchment streams (Stedmon and Markager, 2005). In our mesocosm experiment, C3 and C1 components that were associated with HF additions did not seem to be consumed by planktonic communities or to be abiotically degraded over the experiment. The leaf leachate source, on the other hand, enhanced the microbially derived humic-like C4 component, a fraction also found in yellow poplar leaf litter leachate (Wheeler et al., 2017), and has been observed in other studies using leaf leachate as a DOM source (Cuss & Guéguen, 2015). High amounts of this component have been found across a diverse gradient of lakes in Greenland and it is particularly dominant in soil DOM extracts (Osburn et al., 2017). Other authors have related the C4 component to macrophyte DOM exudates, which are highly reactive to both sunlight and bacterial activity (Cory et al., 2013). To date, only one other mesocosm study has examined specific fluorescence properties (i.e. PARAFAC components) of different allo-OM treatments (Mutschlecner et al., 2018). Those authors identified PARAFAC components in leaf leachates derived from three different types of trees (i.e. alder, poplar, and spruce leaves). Similar to our study, they also found relatively high proportions of more labile DOM molecules, particularly for alder tree leaf leachates, the same type of leaves used in our experiment. Finally, in our study, the fact that all treatments maintained high proportions  Table 1. From high recalcitrant humic-like component (C3), less recalcitrant terrestrial fulvic-like component (C1), microbial derived humic-like component (C4), to more labile proteinlike component (C2). Box = 25th and 75th percentiles, whiskers = 1.5*inter-quartile range. Black line = median of the protein-like C2 component as the experiment progressed is supported by the observed increase in DOC concentrations in all treatments including the control treatment as the experiment progressed. The observed increase in DOC concentration in control enclosures that did not receive any allo-OM additions suggests that internal food web processes such as DOC exudation from algae or excretion from higher trophic levels were important contributors to DOC concentrations, as also reported in another similar experiment (Geddes, 2015). In our study, quantities of Humin-FeedÒ and leaf leachate were added based on previously published experiments utilising these allo-OM sources with similar designs (Geddes, 2015;Rasconi et al., 2015;Hitchcock et al., 2016;Urrutia-Cordero et al., 2017). The HF treatment added only small quantities of DOC, as was expected based on previous studies Minguez et al., 2020), compared with the higher quantities of DOC added from the leaf leachate source as intended based on the experimental design. Therefore, monitoring DOM quality in this study allowed individual DOM constituents in the different treatments to be elucidated, each of which will have differing outcomes and responses for planktonic communities. To this end, considering DOM quality and not just quantity allows for a better understanding of the consequences of increased allo-OM inputs to aquatic ecosystems.

Browning phenomenon
There are a number of approaches which can be used in order to examine the browning effect in aquatic systems which are thought to reduce light conditions for phytoplankton growth. These include the use of covers across the top of the treatment enclosures in order to reduce solar irradiation (Urabe et al., 2002;Villar-Argaiz et al., 2018). However, these methods may not properly mimic potential phytoplankton light limitation in aquatic ecosystems as predicted with increases in water colour. In addition, by using covers, algae growth could be unexpectedly enhanced due to a reduction of photo inhibitory effects, as found in a mesocosm study carried out in a tropical lake (Brighenti et al., 2018). HuminFeedÒ has previously been used to examine the browning effect in mesocosm experiments (Urrutia-Cordero et al., 2017;Lebret et al., 2018) and was, therefore, chosen for this purpose in our experiment. The C1 fulvic-like component clearly increased in our HF treatment and it has been previously reported to have highly photorefractory components because of the rapid attenuation of UVB light in water which increases as bulk DOM photodegrades (Kida et al., 2019). This would reinforce our initial assumption that HF mimics well freshwater browning scenarios. However, in contrast to what we expected, and mainly due to higher DOC pulse additions, the L treatment added more water colour than the HF treatment. Consequently, the L treatment behaved in a similar way to what was expected from the combined HFL treatment (i.e. an increase in water colour and DOC quantity). In turn, the HFL treatment experienced more than double the shading effects compared with the HF treatment. It is clear, therefore, that despite the initial intention of showing the temporal trends in ln-transformed seston C:P ratios during the days of the experiment (DOE). Points are observations for each respective treatment (Cntl blue points, HF grey triangles; L yellow diamonds; HFL red squares). Green dashed horizontal line indicates the Redfield ratio (lnC:P = 4.66) and pink one the threshold elemental ratio for seston P nutrient deficiency, lnC:P [ 5.56 (Healey & Hendzel, 1980). Black dashed vertical line indicates when the allo-OM pulse event occurred. Model summary output is in Table S2. P-values significant codes: \ 0.001 = '***'; \ 0.01 = '**'; \ 0.05 = '*'; non-significant = ns Fig. 5 GAM smoothers (solid lines) showing the temporal trends in ln-transformed seston N:P ratios during the days of the experiment (DOE). Points are observations for each respective treatment (Cntl blue points; HF grey triangles; L yellow diamonds; HFL red squares). Green dashed horizontal line indicates the Redfield ratio (lnN:P = 2.77) and pink one the threshold elemental ratio for seston P nutrient deficiency, lnN:P [ 3.09 (Healey & Hendzel, 1980). Black dashed vertical line indicates when the allo-OM pulse event occurred. Model summary output is in Table S3. P-values significant codes: \ 0.001 = '***'; \ 0.01 = '**'; \ 0.05 = '*'; non-significant = ns the experimental design to increase DOC in the absence of significant increases in colour via the L treatment as per the methodology employed by Hitchcock et al. (2016), colour was, nevertheless, added in substantial concentrations. Colour is not always measured following additions of allo-OM leaf leachate sources in mesocosm experiments (Hitchcock et al., 2016). This may lead to erroneous conclusions regarding planktonic responses, as the effect of increased colour is not separated from effects of increased DOC in these cases. Colour associated with DOC additions from allo-OM sources such as leaf leachates should be determined prior to experimental additions to ensure that the effect of colour following allo-OM pulses can be appropriately assessed. In other studies, tested separate shading effects of leaf leachates as an added DOC source is not generally well achieved due to the use of covers as an extra treatment (Brighenti et al., 2018), which is different from adding colour to the water, or by using a completely different uncoloured labile DOC source such as sucrose as a comparable treatment (Geddes, 2015). On the other hand, there are also studies which only focus on the browning effects, without considering the quality and quantity of potential DOC added, such as when commercially available humic substances are used (Rasconi et al., 2015;Urrutia-Cordero et al., 2017;Lebret et al., 2018). In our study, the aim to disentangle the independent and combined physical and chemical effects of browning (i.e. light limitation) and external energy exports (i.e. carbon and nutrient quantity) was also hardly achieved because of the unexpected high colour added with the leaf leachate source. Future studies investigating how different allo-OM sources impact aquatic ecosystems should more carefully consider DOC-colour relationships of these extracts and compare individual effects by adding the same amount of DOC or colour, aspect not considered in this work.

Seston elemental content and stoichiometry
Despite the large number of studies quantifying seston nutrient content and ratios in lakes (Prater et al., 2017;Bergström et al., 2018), direct analyses of this fraction in mesocosm studies are not that common (Schulhof et al., 2019;Minguez et al., 2020). In our study the HF treatment added seston N, also found in a recently published study (Minguez et al., 2020), whereas the alder tree leaf leachate source increased concentrations of P, as reported elsewhere (Mutschlecner et al., 2018;Navarro et al., 2019). As expected, the HF HF grey triangles, L yellow diamonds, HFL red squares treatment enclosures had higher seston C:P ratios compared to the L and HFL treatments, similar to those found in other experiments using HuminFeedÒ as an allo-OM source analogue (Minguez et al., 2020). High seston C:P values were observed in the HF treatment, but also in the control treatments, implying potential nutrient constraints for primary consumers based on previously published seston nutrient ratio thresholds (Urabe & Watanabe, 1992). Moreover, high seston N:P values throughout the experimental period indicated that the HF and control mesocosms were under severe P limited conditions relative to N, even before starting the allo-OM additions and contrasting with results reported by Minguez et al. (2020) who observed no P limitation in the HF treatment or in the unmanipulated control treatments. In our study, seston C:P and N:P ratios in treatments which received the leaf leachate source (L and HFL enclosures) were lower than the HF treatment. The nutrient concentrations and relative proportions of leaf leachates are highly variable and different from those found elsewhere in leaf litter from the same tree species (Osborne et al., 2007;Schreeg et al., 2013). Due to the higher solubility of P compared to N and C, it is common to observe greater influence of P in leaf Fig. 7 Effect sizes (LRR) of the three different treatments to different water chemistry quantity and quality parameters. a Carbon quantity: DOC, seston C and colour; b nutrient quantity: DIN, seston N, SRP and seston P; c DOM quality: PARAFAC components; and d seston quality: C:P and N:P ratios. Dashed horizontal black lines at 0 indicate the general benchmark (notice than when available we used benchmark of day 0, close to 0 but not represented visually. See Table S4). P-values significant codes differing from benchmark values: P \ 0.001 = '***'; P \ 0.01 = '**'; P \ 0.05 = '*'; P [ 0.05 = 'ns', non-significant (Table S5). Different letters between treatments indicate significant differences (P-value \ 0.05) amongst them (Table S6). Box = 25th and 75th percentiles, whiskers = 1.5*inter-quartile range. Black line = median leachates directly affecting nutrient ratios (Schreeg et al., 2013).
Significant non-linear trends were observed for both ratios (C:P and N:P) during the 36 days of experiment. The pulse disturbance event clearly affected these nutrient ratio values. Seston C:P reached minimum values at day 4 but increased and stabilised again towards the end of the experiment. Similarly, seston N:P values were lowest directly following the allo-OM additions. In natural freshwater ecosystems, seston C:P and N:P ratios tend to decrease during mixing and the rainy season, or in autumn and winter, when there is likely a greater nutrient supply to the lake surface mixed layer as a consequence of increased fluvial inputs, whilst higher values are often associated with spring and summer (Prater et al., 2017;Ngochera & Bootsma, 2018). Most of these studies relate seasonal patterns in lake seston quality with changes in particulate P supply and subsequent P-uptake during the growing season, despite sitespecific physical and biological drivers (Elser et al., 1995;Hessen et al., 2005). Overall, extensive temporal variations of these ratios suggest the need to better understand mechanistic changes in the biochemistry of lake seston (Kreeger et al., 1997;Calderó-Pascual et al., 2020). In our study, the unexpected decrease of seston C:P ratio in all treatments the first week after the pulse could be explained by 1) relatively high P additions in the treatments receiving the P rich leaf leachate source and 2) the collapse of primary producers indicated by Chl-a in the control and HF mesocosms, which accounted for rapid declines in seston C relative to seston P. Seston C concentrations in the control enclosures gradually decreased during the experiment to below 0.5 mgC l -1 , the threshold at which Daphnia growth is considered food-limited (Lampert, 1987;Urabe et al., 2002). This suggests that regardless of food quality in terms of carbon constituents and elemental ratios, food quantity for primary consumers was likely to have declined in the control treatment as the experiment progressed to levels generally considered insufficient to support higher trophic levels. As the experiment progressed, increases in seston elemental ratios to pre-disturbance conditions, and even higher values, would suggest the uptake of the available dissolved P by primary consumers, which were likely already highly P limited even before the pulse disturbance event took place.

Mimic natural allo-OM sources
In our experiment, in an effort to mimic natural additions of nutrients following inputs of allo-OM and to optimise the number of treatments versus the number of replicates per treatment, inorganic nutrients were not added to match those introduced with the  Table S7. P-values significant codes: \ 0.001 = '***'; \ 0.01 = '**'; \ 0.05 = '*'; non-significant = ns HuminFeedÒ and leaf leachate sources, which contrasts with other authors who added inorganic nutrients to the enclosures (Lefébure et al., 2013;Geddes, 2015;Hitchcock et al., 2016;Urrutia-Cordero et al., 2017) or created an extra treatment controlling the single effect of nutrient additions (Brighenti et al., 2018). Consequently, in this study, the control and HF treatments experienced extreme P limitation over the course of the experiment, potentially masking the directly tested browning effects, as found in other studies (Tonetta et al., 2018). Prior to the described experiment, we aimed to prepare a local (i.e. top bare soil) allo-OM stock solution. However, it was logistically complex to obtain the designed high DOC concentrations for the targeted treatment. Therefore, for practical reasons, we collected alder tree leaves from a national park (40°39 0 22.1'' N, 31°37 0 47.89'' E), in Bolu, Turkey. This source is unlikely, therefore, to be directly representative of natural allo-OM additions to the lake in which the mesocosms were located. Efforts to mimic natural allo-OM additions would have been improved by preparing and characterising treatment stock solutions similar to the surroundings catchment sources. In turn, new studies considering other types of allo-OM such as agricultural soil (Liess et al., 2015), different types of leaf leachates (Lennon & Cottingham, 2008), and burned soil leachates (Cunillera-Montcusí et al., 2019) should be further explored as common catchment type systems affecting aquatic environments. However, more attention should be paid to mimicking natural inputs of DOM. Realistically, an increase of 8 mgC l -1 as DOC in pelagic systems is not a predicted scenario for lakes owing to higher buffering capacity and lower resident time (Leech et al., 2018) compared to concentration in streams, which often show high peaks in POC and DOC concentrations after flooding events (Dhillon & Inamdar, 2013;Hitchcock & Mitrovic, 2015). Other authors mimicking potential predicted scenarios ranged from ? 1.2 to ? 4.8 mgC l -1 (Sanders et al., 2015;Brighenti et al., 2018;Lebret et al., 2018), which places their maximum DOC additions at half the value of the additions with the leaf leachate source in this study which were used to follow intermediate concentrations used by Hitchcock et al, (2016).

Overall treatment effects
In this study, the LRR effect sizes approach recommended by Hillebrand & Gurevitch (2016) was utilised in order to contrast impacts of different allo-OM sources (i.e. examination of changes with reference to the control treatment). The HuminFeedÒ effects were mainly lower, and even non-statistically significant for some of the variables studied, compared to the stronger effects arising from the additions of the leaf leachate source. The only variables to have the same strong effects in the HF treatment as observed in the L treatment were DOM quality parameters (i.e. PARAFAC components). In turn, the HFL treatment experienced additive effects, which were different to our expectations of finding multiplicative influences on the outcomes. The HF treatment showed lower seston nutrient ratio effects relative to the L and HFL treatments. In addition, in this case, the HF treatment showed positive effects (i.e. increased overall seston C:P and N:P) whilst the L and HFL treatments' effects were negative (i.e. decreased overall seston C:P and N:P). For all these reasons, the overall effect sizes of the combined HFL treatment were mainly influenced by the presence of the leaf leachate source, since HuminFeedÒ acted more like the control and different from the expected synergistic effects of the combined sources.

Implications for primary productivity
Different from what we hypothesised, we found maximum Chl-a concentrations in the HFL treatment followed by the L treatment, whilst Chl-a concentrations and dynamics in the HF treatment did not differ from controls. Previous studies have simulated increases in water colour using HuminFeedÒ additions and have shown that, similar to what we found but different to what we expected, the respective percentage increase in water colour did not have any effect on total phytoplankton biomass (Ratcovich, 2014;Urrutia-Cordero et al., 2017;Lebret et al., 2018). Our results suggest that, as expected, primary producers were mainly stimulated by P additions associated with the leaf leachate source. Therefore, it would appear that the addition of nutrients exerted a greater influence than light reduction due to browning, as was also reported by other authors (Kankaala et al., 2010;Sanders et al., 2015;Hitchcock et al., 2016;Degerman et al., 2018;Tonetta et al., 2018), and especially in already P limited ecosystems (Villar-Argaiz et al., 2018), as in our investigated system. In addition, although there was no phytoplankton community composition data available for our experiment, other authors found greater abundance of mixotrophic, e.g. Cryptophytes and nanoplanktonic flagellates in the enclosures receiving leaf leachate on top of HuminFeedÒ in contrast to clear controls and brown HF individual treatments (Fonseca et al., 2022). This would suggest that, in our study, a potential increase in mixotrophic algae in the treatments receiving the leaf leachate source could explain the unexpected maximum increase in Chl-a response in the HFL treatment. Regarding reported overall temporal dynamics, minimum Chl-a concentrations observed around day 12 in all mesocosms could be explained by high grazing pressure since the beginning of the experiment as this decrease does not seem to be related with differences in treatment water chemistry, but more influenced by biological interactions (Yıldız et al., in prep).

Conclusion
Alterations in the supply of allo-OM to freshwater lakes due to global change is evident. However, there is still considerable uncertainty in consolidating projected scenarios for in-lake processes. Our study highlights that it is important to characterise the resources used to mimic allo-OM inputs in mesocosm experiments in order to understand and predict the impacts to aquatic ecosystems. In this study, the fact that leaf leachate treatments added P to a likely P limited system suggests that allo-OM disturbance events could increase nutrient availability to nutrient limited planktonic communities, as observed in other temperate lakes (Senar et al., 2021). On the other hand, the same events could enhance eutrophication in water bodies which are already nutrient sufficient, as reported elsewhere (Gao et al., 2021). In contrast, other types of allo-OM sources, which are more similar to HuminFeedÒ characteristics, such as those associated with humic waters (Hamdan et al., 2021), could lead to higher nutrient constraints (i.e. high C relative to N and P) at all trophic levels. Our results demonstrated that, the HF treatment behaved as expected compared to the unexpected increases in water colour with the L treatment. In addition, the HF treatment overall LRR effect sizes were weaker than expected, but stronger for the L treatment which added high amounts of P, adding complexity to the capacity of the experiment to disentangle browning effects from those of direct external C inputs effects. However, by characterising the nutrient components within the treatment, some of this complexity can be directly addressed. The indirect effect of nutrient stoichiometry associated with dissolved and particulate allo-OM inputs as well as DOM quality, in addition to browning and carbon quantity repercussions, are all key factors regulating the ecological responses to projected increases in allo-OM in freshwater systems and should be characterised in future mesocosm experiments.
Availability of data and material The datasets generated and analysed during the current study are available within this article and its supplementary information files through the presented figures and tables. Raw experimental metadata will be publicly available after three years of the experiment (i.e. 2023) through the AQUACOSM mesocosm data repository (https://www. aquacosm.eu/data/) but can be available from the corresponding authors on reasonable request.
Code availability All code available can be found in electronic supplementary material.

Declarations
Conflict of interest The authors declare no conflict of interest.
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.