Ecology

SettingThis study was conducted in approximately 40-year-old naturally regenerated P. densiflora stands in Wola National Experimental Forest in Gyeongnam province in South Korea (35°12′ N, 128°10′ E; Table 1). The productivity of this forest is low, with a dominant tree height of 10 m at 20 years of age. Over the last 10 years, the mean annual precipitation was 1490 mm, of which one third fell during summer (July–August), and the mean temperature was 13.1 °C. The vegetation growing season generally lasts for approximately 200 days, extending from early April to October. The soil texture is a silt loam originating from sandstone and shale (clay 13.0 ± 0.8%, silt 44.1 ± 1.3%, sand 42.9 ± 1.6%; n = 18). The given texture results in volumetric water contents at 13.4 ± 0.7% (m3 m−3) at permanent wilting point (1500 kPa) and 40.7 ± 1.2% at field capacity (10 kPa)55. The understory is covered with Lespedeza spp., Quercus variabilis, Q. serrata, Smilax china, and Lindera glauca.In 2010, we selected two adjacent P. densiflora stands approximately 100 m apart from each other (180 m and 195 m above sea level, on slopes of 15° and 33°, both stands face south). Following a completely randomized design, we established nine plots (10 × 10 m2 with a 5 m untreated buffer) within each stand, of which three were randomly assigned to annual NPK fertilization, three to PK fertilization, and the rest to a control treatment without fertilization. The fertilizer, composed of urea, fused superphosphate and potassium chloride (N3P4K1) or P4K1 was added manually by deposition on the forest floor for 3 years in April 2011, April 2012, and March 2013. Over these 3 years, the NPK plots received 33.9 g N, 45 g P, and 11.1 g K m−2, while the PK plots received 45 g P and 11.1 g K m−2.Tree and stand measurementsThe standing biomass of trees was estimated using a combination of site-specific allometric equations based on destructive harvesting56 and repeated measurements of the dimensions of all trees in each plot (5–18 trees plot−1). The stem diameter at 1.2 m (D) was measured for all trees in each plot for which D was ≥ 6 cm. Selecting a representative tree in size for each plot within the 4 × 4 m2 center of the plot, we measured the tree height (H) and crown base for the representative trees. Measurements were performed in April and September 2011, September 2012–2014, and October 2021. We observed no effect of fertilization on the relationship between D and H or between D and crown base, so we assumed no effect on the allometric functions for foliage or branch biomass. A dendrometer band (Series 5 Manual Band, Forestry Suppliers Inc., Jackson, MS, USA) was installed on 18 representative trees (one per plot) to monitor radial growth monthly.Three 0.25 m2 circular litter traps were installed 60 cm above the forest floor in each plot in April 2011. Litter was collected at 3-month intervals between June 2011 and March 2015. The litter from each trap was transported to the laboratory and then oven-dried at 65 °C for 48 h. All dried samples were separated into needles, bark, cones, branches, and miscellaneous components, and weighed separately.In September 2014, we estimated the biomass of understory vegetation, separately for woody plants and herbaceous plants. All woody plants  More

Description of the phenological, vegetative and yield traits of the accessions per habitatThe process of data management included the computation of mean squares for the assessed phenological, vegetative and yield traits of the accessions with the sampling habitats considered as the treatment factor. The error mean squares served in the multiple comparison of means reported in Table 1.Table 1 Means of the measured phonological, vegetative and flowering and yield traits of Cucurbita moschata genotypes sampled from seven habitats.Full size tableRegarding the phenological traits, the accessions from the habitat of Zh have the longest period from seeding to first male (102.39 d) and first female (108.14 d) flower appearances, and the longest period from seeding to physiological maturity (153.95 d). For those traits, the accessions from Tiassale and Soubre are not significantly different from those of Zh. And, accessions from Tiassale and Zh have the longest periods from seeding to 50% flowering. On the other hand, accessions from Korho, Ferke, Bondu and Burki develop their first male and female flowers and attain 50% flowering in a very short period. They also reach physiological maturity faster. Accessions from Korho, however, have the longest period from seeding to 50% emergence (6.07 d) and accessions from Bondu have the longest period from first female flower appearance to physiological maturity (53.04 d).
For the vegetative traits, accessions from Tiassale and Soubre have the largest girth size (4.43 cm and 4.63 cm, respectively). Accessions from Tiassale have the longest (24.98 cm) and widest (19.94 cm) leaves, the longest male (16.2 cm) and female (4.03 cm) peduncles and the longest petioles (34.94 cm). The measures for those organs on accessions from Soubre rank second to those of Tiassale. On the other hand, accessions from Korho, Ferke, Bondu and Burki are characterized by smaller girth size, smaller leaves, smaller petioles and smaller peduncles of male and female flowers. But the accessions from Bondu are the tallest (586.91 cm) followed by the accessions from Ferke (489.20 cm). And the accessions from Zh are the shortest (417.38 cm).For the flowering and yield traits, accessions from Tiassale and Soubre show the largest numbers of male (27.33 units and 22.58 units, respectively) and female (5.22 units and 6.05 units, respectively) flowers per plant, largest numbers of fruits per plant (2.78 units and 2.53 units, respectively) and largest measures of all fruit-related traits. Their seeds are very large, but in small numbers. In contrast, accessions from Korho, Ferke, Bondu and Burki have the smallest numbers of male and female flowers per plant, the smallest numbers of fruits per plant and the smallest measures of fruit-related traits. They have large numbers of seeds, but their seeds are smaller, except the seeds of the accessions from Burki. Refer to Table 1 for more detailed information.
Variability of the phenological, vegetative and yield traitsTable 2 shows the spread of the phenological and morphological traits of the assessed accessions of C. moschata. All the evaluated traits showed very wide ranges of distribution of the observations. Some conspicuously wide ranges of traits include number of days to 50% flowering (DTF) that goes from 52 to 152 d, plant height with a minimum of 48 cm and a maximum of 1510 cm, diameter of the fruit that is between 5.8 cm and 35 cm, weight of the fruit that varies between 150 g and 10,930 g and number of seeds per fruit that spreads in the interval from 32 units per fruit to 729 units per fruit. Excluding the number of days to 50% emergence (DTE), all the other assessed traits have remarkably wide ranges of phenotypic expressions (Table 2). All the traits but DTE, DTF, days from first female flower appearance to fruit maturity, fruit length and length of the dry seed, had outliers. The number of outliers ranged from 1 to 67. Except the outliers observed with the width of the dry seed, all the outliers were above 1.5*IQR + Q3 where IQR is the inter-quartile range and Q3 is the third quartile. The presence of outliers is indicative of the richness and large variability of the population of accessions. The outliers are exceptional performances that fall outside the normal distribution of the observations. They are a stock of unusual traits that can be used in a crop improvement program when beneficial. For example, the observed outliers for diameter of the fruit, weight of the fruit or thickness of the pulp can be used in a breeding program for the improvement of fruit yield. Similarly, outliers for beneficial traits related to the seed can be used to improve C. moschata crop for seed yield. Besides, the computed mean squares (data not reported) showed highly significant variations between accessions for the assessed traits. They all yielded p-values less than 0.01, providing additional support to the evidence of large variability among the accessions of C. moschata of Cote d’Ivoire. The computed standard deviation, and median absolute deviation for each trait are additional evidence. We should note that in most cases, the mean squares associated to year (data not reported) were not significant, indicating the relative stability of the assessed traits.Table 2 Minimum (Min), first quartile (Q1), median, third quartile (Q3), maximum (Max), standard deviation (SD), median absolute deviation (MAD) and outliers obtained from the phenological, vegetative and flowering and yield traits of 663 accessions of C. moschata.Full size tableThe components of variance, the quantitative genetic differentiation, the overall mean, and the coefficients of variation are reported in Table 3. The lme4 package37 used in the determination of the components of variance, does not provide p-values in the analysis of mixed or random models. The reported quantities in Table 3 are not accompanied with tests of significance. It is worth mentioning that the respective units of measure of the assessed traits are squared for the variances and the evaluated estimates will be reported without the units of measure. The phenotypic variance ((sigma_{p}^{2})) is partitioned into variance between morphotypes or genotypic variance ((sigma_{g}^{2})), and within morphotypes or residual variance ((sigma_{e}^{2})). For the class of phenological traits, considerable genotypic variances were observed with days to 50% flowering (266.21) and days to first male flower appearance (254.40), compared with their respective residual variances (148.13 and 199.50). Regarding the class of vegetative traits, only the peduncle length of male flowers had a genotypic variance (9.22) greater than its residual variance (8.86). In the class of flowering and yield traits, 8 of the 15 traits assessed showed large genotypic variances in comparison with their respective residual variances. They are number of female flowers per plant ((sigma_{g}^{2}) = 3.02 versus (sigma_{e}^{2}) = 2.36), length of the fruit ((sigma_{g}^{2}) = 53.96 versus (sigma_{e}^{2}) = 48.97), diameter of the fruit ((sigma_{g}^{2}) = 37.17 versus (sigma_{e}^{2}) = 16.76), volume of the fruit ((sigma_{g}^{2}) = 10,713,468 versus (sigma_{e}^{2}) = 3,904,590), weight of the fruit ((sigma_{g}^{2}) = 5,413,819 versus (sigma_{e}^{2}) = 1,420,187), diameter of the cavity enclosing the seed ((sigma_{g}^{2}) = 19.12 versus (sigma_{e}^{2}) = 7.75), thickness of the fruit pulp ((sigma_{g}^{2}) = 1.11 versus (sigma_{e}^{2}) = 0.94) and weight of the fruit pulp ((sigma_{g}^{2}) = 5,979,212 versus (sigma_{e}^{2}) = 1,088,750). For a trait to have a lager genotypic variance than the residual variance is synonymous to a relative ease of improvement of the crop for that trait through a breeding program.Table 3 Components of variances ((sigma_{p}^{2}), (sigma_{g}^{2}), (sigma_{e}^{2}), (sigma_{a}^{2})), quantitative genetic differentiation ((Q_{ST})), overall mean ((mu)), and coefficients of variation (%) ((CV_{p}),(CV_{g}),(CV_{e})), of the measured phenological, vegetative and yield traits of the accessions of C. moschata of Cote d’Ivoire.Full size tableThe coefficient of variation (CV) is another statistic that measures variation. It is actually the dispersion of a trait per unit measure of its mean, which can be used to compare variations of traits with different measurement units or different scales. As a rule-of-thumb, a coefficient of variation greater than 20% is indicative of large variation for the trait. The phenotypic coefficient of variation is considerably high for 25 of the 28 assessed traits. Only the number of days from seeding to physiological maturity, the first and second longest axes of the dry seed show coefficients of variation less than 20%. Traits with very large phenotypic coefficients of variation include the peduncle length of female flowers ((CV_{p}) = 93.98%), weight of the pulp ((CV_{p}) = 92.96%), volume of the fruit ((CV_{p}) = 89.17%), weight of the fruit ((CV_{p}) = 78.30%) and number of female flowers per plant ((CV_{p}) = 65.81%). With respect to the residual coefficients of variation, only the number of days from seeding to 50% emergence and number of days from first female flower appearance to physiological maturity have residual coefficients of variation greater than 20%, among the phenological traits. All the vegetative traits have residual coefficients of variation greater than 20%, and show a near-perfect linear relation (r = 0.98; p  More

Preparation of the pao cai soupFirst, 35 kg of white radish (Raphanus sativus), 35 kg of cabbage (Brassica oleracea), 2 kg of chili pepper (Capsicum frutescens), 1 kg of ginger (Zingiber officinale), 1 kg of peppercorns (Zanthoxylum bungeanum), 2.5 kg of rock sugar, and 210 kg of cold boiled water (containing 6% salt) were divided into six ceramic jars. After 7 days of natural fermentation at room temperature, the pao cai was filtered out with sterile gauze to obtain 200 kg of pao cai soup. To ensure an even distribution of microorganisms in the soup, the soup was mixed well and then left to rest for 12 h, the supernatant was taken, and the soup was left to rest for 12 h again.The plants used in this study were cultivated vegetables which purchased from the vegetable market at the study site. All local, national or international guidelines and legislation were adhered to in the production of this study.Establishment of the microcosm systemSeventy-eight for each size of 10 ml, 20 ml, 50 ml, 100 ml, 250 ml, 500 ml, and 1000 ml sterile glass culture flasks were filled with pao cai soup, the bottle mouth was sealed with sterile sealing film, and the bottle was capped without leaving any air (Fig. 1). Each flask became a microcosm and was cultured in a 25 °C incubator.Figure 1Schematic diagram of the establishment of the microcosmic system.Full size imageSample collectionBefore the microcosm system was established, a sample of well-mixed pao cai soup was taken as a reference to establish background biodiversity. The microbial community dynamics should change the fastest at the beginning of the microcosm system establishment and gradually become slower over time. Considering the workload and cost, this study collected samples daily for 1–10 day after the establishment of the microcosm and then collected every 2 days for 10–30 day and every 5 days for 30–60 day. Three different microcosms of the same volume were established. Monitoring was carried out for 60 days, and a total of 546 samples of 7 volumetric gradients were obtained at 26 time points. At the time of sampling, the pao cai soup in the microcosm was mixed, and 50 mL of sample (10 mL of sample was collected for microcosm systems with a volume of less than 50 mL) was collected. The sample was centrifuged at 8000 rpm for 10 min, the supernatant was collected for pH determination, and the pellet was stored in a − 80 °C freezer.Microbial analysesMicrobial DNA was extracted from pao cai samples using the E.Z.N.A.® Soil DNA Kit (Omega Biotek, Norcross, GA, U.S.) according to the manufacturer’s protocols. For bacteria, we targeted the V3-V4 region of the 16S ribosomal RNA (rRNA) gene, using the 338F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806R (5′-GGACTACHVGGGTWTCTAAT-3′) primer pairs. For fungi, we targeted the ITS1-1F region of the nuclear ribosomal internal transcribed spacer region (ITS rDNA) gene, using ITS1-1F-F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS-1F-R (5′-GCTGCGTTCTTCATCGATGC-3′). PCRs were performed in triplicate in a 20 μL mixture containing 4 μL of 5 × FastPfu Buffer, 2 μL of 2.5 mM dNTPs, 0.8 μL of each primer (5 μM), 0.4 μL of FastPfu Polymerase and 10 ng of template DNA. The PCR program for the 16S rRNA gene was as follows: 3 min of denaturation at 95 °C; 27 cycles of 30 s at 95 °C, 30 s of annealing at 55 °C, and 45 s of elongation at 72 °C; and a final extension at 72 °C for 10 min. For the ITS1-1F region, the PCR program was as follows: samples were initially denatured at 98 °C for 1 min, followed by 30 cycles of denaturation at 98 °C for 10 s, primer annealing at 50 °C for 30 s, and extension at 72 °C for 30 s. A final extension step of 5 min at 72 °C was added to ensure complete amplification of the target region. The resulting PCR products were extracted from a 2% agarose gel, further purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™-ST (Promega, Madison, WI, USA).Purified amplicons were pooled in equimolar amounts and paired-end sequenced (2 × 300) on an Illumina NovaSeq platform (Illumina, San Diego, CA, USA) according to standard protocols. The analysis was conducted by following the “Atacama soil microbiome tutorial” of QIIME2 docs along with customized program scripts (https://docs.qiime2.org/2019.1/). Briefly, raw data FASTQ files were imported in the QIIME2 system using the qiime tools import program. Demultiplexed sequences from each sample were quality filtered, trimmed, denoised, and merged, and then the chimeric sequences were identified and removed using the QIIME2 DADA2 plugin to obtain the feature table of amplicon sequence variants (ASVs)24. Compared with traditional OTU that clusters at 97% similarity, ASV has higher accuracy, equivalent to 99% similarity clustering. The QIIME2 feature-classifier plugin was then used to align ASV sequences to the pretrained GREENGENES 13_8 99% database (trimmed to the V3-V4 region bound by the 338F/806R primer pair for bacteria) and UNITE database (for fungi) to generate the taxonomy table25. Any contaminating mitochondrial and chloroplast sequences were filtered using the QIIME2 feature-table plugin. Based on the sequence number of the lowest sample, perform the resampling to make the sequence number equal for each sample. Due to the random nature of sequencing, ASVs specific to each sample in this study were present. To reduce the uncertainty introduced by the sequencing process, we filtered out rare ASVs with less than 0.001% of the total sequence volume.Data analysisIn this study, the data of fungi and bacteria were integrated and analyzed, and all microbial diversity appearing in the text represent the sum of all fungi and bacteria. Species richness is equal to the number of taxa, which is equal to the total number of all bacterial and fungal ASVs. The vegan package in R 4.2.1 was used to calculate the species richness of each sample based on the ASV feature table26. Using flask volume instead of area, SAR fitting was performed using a semi-logarithmic model, and its significance was tested. The semi-logarithmic model is the function S = c + b*logA, where S is species richness, A is area (in this case, volume is used instead), and b and c are fit parameters27.The microcosmic system in this study is hermetically sealed, and all microorganisms originate from a single portion of well-mixed paocai soup (ie species pool). The speciation process in the 60-day experimental system should be negligible due to the short experimental period. The extinction rate of a microcosm system is equal to the number of ASVs lost in the microcosm system compared to the species pool divided by the total number of ASVs in the species pool. The extinction rate is the number of extinct ASVs in each system compared to the species pool. Pearson correlation analysis was performed with volume as the independent variable and extinction rate as the dependent variable to determine the correlation between volume and extinction rate at each time point. When microorganisms of a microcosmic system disappear entirely or cannot be detected, the microcosm is recorded as an annihilated microcosm. The annihilation rate at a time point is equal to the number of microcosms annihilated at that time, divided by the total number of microcosms. The difference between the extinction rate and annihilation rate defined in this paper is that the extinction rate is for ASVs within each sample, and the annihilation rate is for microcosmic system at each sampling time point. The two indicators jointly characterize the local extinction of microorganisms from different perspectives. Non-linear regression with a bell-shaped form was performed with time as an independent variable and pH and annihilation rate as dependent variables, and regression lines were plotted based on R 4.2.1.According to the taxonomy table, bacterial ASVs were divided into acid-producing and non-acid-producing categories, and their extinction rates were calculated separately. The agricola, ggplot2, vegan and ggpubr packages were used to draw alpha diversity box plots and perform the Wilcoxon rank sum test for differences between groups26,28,29,30. Non-metric multidimensional scaling (NMDS) analysis was performed with the vegan package based on Bray–Curtis dissimilarity. In addition, the potential Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologue (KO) functional profiles of microbial communities were predicted with PICRUSt31. Resistance-related genes were screened using the gene function predictions. The relationship between the relative abundance of resistance-related genes and the volume of the microcosm was analysed by Pearson correlation, and a forest map was plotted to present the results. More

We investigated the complex relationships between the environmental characteristics of blue spaces and visit-related mental well-being in a multi-country study including 17 bluespace types and four facets of subjective well-being. Our aim was to operationalise, and consider the utility of, the Bratman et al.9 conceptual model that links ecosystem services (ESS) with mental health.Consistent with the proposed conceptual model, mental well-being outcomes relied on a complex interplay of individual, environmental, and visit characteristics.Summary of findingsOverall, bluespace visits were associated with better subjective mental well-being outcomes if the visits took place in nearby coastal areas or rural rivers, were perceived as safe and to have good water quality, and had a long duration. They could involve a range of activities such as playing with children, socialising, or walking. The degree to which the perceived presence of wildlife predicted visit satisfaction varied depending on the bluespace type, suggesting that the importance of ecosystem features such as biodiversity may vary by the setting.We can also identify the combination of environmental and visit characteristics associated with particularly high levels of well-being for a specific outcome. For example, an optimal visit in terms of happiness might be to sandy beaches where there are high levels of perceived safety and excellent water quality; with a visit lasting at least three hours; and possibly involving playing with children, socialising, sunbathing/paddling and/or walking with a dog; and has short travel times that do not involve public transport.RQ1—natural and environmental featuresResearch question 1a—Which bluespace type(s) were associated with the highest levels of recalled visit mental well-being?Four of the five bluespace types associated with the highest levels of visit satisfaction were coastal (sea cliffs, rocky shore, sandy beaches, rural river and seaside promenade), indicating that these environments may be particularly beneficial for well-being. Visits to these environments were also associated with the lowest levels of visit anxiety, with the exception of seaside promenade and sea cliffs, which were not significantly different to the grand mean. Seaside promenade was the only urban environment in the top five.In addition, only coastal sites were associated with significantly higher levels of visit happiness (compared to the grand mean), further highlighting the potential importance of these environments. Although not explored here, coastal scenes tend to be associated with particularly high aesthetic and scenic value25,26 which may also be positively related to subjective well-being.These findings are broadly consistent with other studies from the UK17,27, but are extended here to our international sample. White et al.28 also used data from the BlueHealth International Survey (BIS) and explored visit frequency to different environments and associations with general mental health and well-being outcomes, including the World Health Organisation five-item Well-being index referring to the two weeks prior to the survey. Consistent with the results here, they found that visit frequency to “coastal blue” environments was more strongly associated with psychological well-being in general than visit frequency to “inland blue” environments. Our study adds to these more general findings by showing that these associations may come as a direct result of the recalled well-being experienced on specific visits to these locations.Confidence in our results was strengthened as we included general mental well-being in our analysis to adjust for whether happier people tend to visit sandy beaches, for example. The results for visit anxiety were not always the inverse of the trends observed in the positive measures of well-being, supporting the need to look at multiple aspects of mental well-being when considering the effects of nature contact.Research question 1b—Which bluespace qualities were associated with the highest levels of recalled visit mental well-being?Of the range of qualities that we investigated as predictors, perceived safety and ‘excellent’ water quality (vs. ‘sufficient’) consistently exhibited the strongest relationships with subjective mental well-being. Perceived safety has been found to be important when visiting blue spaces in several qualitative studies29,30,31, as well as a quantitative study with older adults in Hong Kong14. Blue spaces have particular safety issues with respect to drowning32,33, but fear of crime29,30,33 or pedestrian safety34 may also be relevant.Water quality has also been found to be important in previous economic valuation studies of recreational use and enjoyment of lakes and estuaries in the USA and Australia35,36 as well as a contingent behaviour experiment carried out as part of the BlueHealth International Survey (in European countries only)37. We recognise that here we used a metric of perceived water quality, rather than measures based on biological or toxicological sampling. Nevertheless, perceptions have been reported to positively correlate with sampled water quality parameters38, although assessments can vary systematically such as by bluespace type39. Highly visible harmful algal blooms, for instance, have also been found to affect experiences of blue spaces40.Further, and again consistent with earlier work15,41,42, the presence of facilities and wildlife, and absence of litter, were generally associated with better subjective mental well-being. Both perceived presence of wildlife and facilities were also associated with higher levels of anxiety however, indicating complexities between environmental qualities and well-being. Some wildlife may be deemed unpleasant or an ecosystem disservice, for example. The presence of good facilities may indicate the presence of more people; and visitor density in natural environments can be related to preference43. These results highlight the importance of environmental quality and not just type, consistent with other frameworks12,37.Research question 2—How is exposure, as operationalised by visit duration, related to recalled visit mental well-being?Broadly consistent with research in the green and bluespace literature14,17,44, we found that mental well-being outcomes were generally higher with greater exposure as indicated by visit duration. For decreasing visit anxiety, this was only significant when visits were longer than an hour and a half. As we did not measure pre-visit anxiety levels, we are cautious about identifying this as a potential temporal threshold for reducing anxiety at this stage.Similarly, also using the BlueHealth International Survey, White et al.28 found that well-being outcomes were higher with greater visit exposure to green and blue spaces using a metric of visit frequency. However, in contrast to this and other research which looked at overall weekly aggregated time in nature (e.g.28,45), we have no evidence of diminishing marginal returns as the effect sizes associated with specific visit duration continued to increase with increasing duration.Research question 3—What experiences in blue spaces, in terms of activities (3a) and companions (3b), are associated with the most positive recalled visit mental well-being outcomes?Although walking was the most popular activity, the activity with the highest mental well-being ratings was playing with children, especially in certain locations such as beaches (Fig. 4). However, we also find that anxiety tended to be higher when children were present. We speculate that the purpose of the visit may be important. For example, many who go to the beach with children do so in order to play. However, if children are present on more adult-oriented activities such as hiking, this may increase adult anxiety during the visit. From a representative sample of English adults, White et al.17 found that recent nature visits with children were associated with the lowest levels of well-being. Therefore, visits with children may be associated with a more complex set of emotions, being both slightly more stressful, but also potentially more rewarding and ‘meaningful’46. Ecosystem features of beaches may be particularly supportive of high well-being activities. A qualitative study in the UK, for instance, highlighted the particular opportunities for adults and children to play together at the beach, including rock-pooling and making sandcastles as well as water-based activities47.Visits with other adults were associated with higher levels of both visit satisfaction and worthwhile-ness, and socialising as an activity was associated with better visit well-being for all outcomes compared to the grand mean. This is consistent with studies using the day reconstruction method, which link activities with experiential well-being, in the USA48 and Germany49 where socialising was associated with the highest, or second highest, levels of well-being for all the activities assessed. Further, social interactions have been recognised as an important benefit by many of those visiting freshwater blue spaces in a previous study18.Research question 4—Does the relationship between wildlife presence and recalled visit well-being vary by bluespace settings?The relationship between the presence of wildlife and visit satisfaction varied with bluespace type. The strongest positive association was found for fen, marsh and bog areas, which may also be related to the purpose of visit. For instance, those who visit places such as fens, marshes and bogs, may do so for the explicit purpose of observing wildlife (often birds) and the presence of wildlife would therefore be important for satisfaction with the visit.Perceptions towards wildlife have been found to vary by location in other studies. For example, in Sweden, greater prior experience with geese at beaches was associated with a negative attitude towards geese50. Further, the species present are likely to vary across different environments. In three urban areas in the UK, green spaces correlated with the abundance and species richness of birds considered to provide cultural services (songbirds and woodpeckers), while an abundance of birds considered to provide disservices (e.g. some gull species, feral pigeons) was independent of green spaces51. Preferences for some species over others may explain some of the negative or null relationships between the presence of wildlife at different blue spaces. These examples from the literature, alongside our own results, indicate the potential for benefits from the management of wildlife for psychological ecosystem services differentially across environments, although these should be considered alongside other conservation and ESS goals.MechanismsSeveral mechanisms potentially explain the beneficial effects of visiting blue spaces on mental health and well-being12, including the provision of opportunities for physical activity52,53; social interaction18; cognitive restoration and stress reduction17,54; emotion regulation55 and connecting with nature12. Consistent with these mechanisms, we found that respondents were using blue spaces for both physical activity and social interaction; and that playing with children and socialising were associated with particularly high levels of well-being.In addition to the positive association we find between some ESS and well-being, including presence of wildlife and water quality, additional bluespace ESS not considered here, may also affect mental health and well-being12. For example, the provision of a cooling effect56 and air pollution mitigation57.Strengths and limitationsA key strength is our operationalisation of the Bratman et al.9 conceptual model for mental health using data from a large, 18 country survey that included 17 different bluespace types, five quality metrics and four subjective mental well-being outcomes. The relatively high explanatory power of our models suggests all the variables we explored were important for subjective well-being.Despite the strengths, however, there were also several limitations. The survey was cross-sectional and causality cannot be inferred. For example, happier people may choose to visit a beach rather than another location, although we also controlled for general levels of subjective mental well-being in an attempt to control for this possibility (See Supplemental Materials). Further, although the majority of respondents (53%) recalled a visit within the last 7 days, some were recalling visits up to a month ago, with potential memory biases increasing in line with length of recall.Although our data were collected by an international market research company to be representative of age, gender and region within country, our online sample may not be fully representative across more characteristics and any country-level conclusions need to be treated with caution. We also acknowledge that there were no results from Africa, the Middle East or South America; and Hong Kong was the only representative from Asia. This suggests far more research is needed in other regions to better understand how bluespace ecosystems interact with subjective well-being globally.There may also be socioeconomic confounds that we did not include in our models which may account for some of the effects. Not everyone visits nature for recreation58, including about 4000 people here who did not visit a bluespace in the four weeks prior to responding to the survey. Some groups may therefore have been under-represented; and we should be careful in assuming that our findings generalise to all sub-population groups.Nevertheless, our visit sub-sample distributions were generally similar to that of the weighted percentages in the full sample, with the exception of age where those aged over 60 were under-represented (Table S2); therefore, we suspect these issues were not too influential for the overall results, although care needs to be extended to inferences with respect to older adults.A further limitation was that we only considered the qualities of places where people reported making recreational visits, with respondents presumably less likely to visit places where they feel really unsafe or lacking in facilities29. Further research may want to study responses to a broader range of bluespace settings, including those that are less visited, to determine the generalisability of the generally positive results found here. Such studies could use pre-existing tools to objectively assess the quality of blue spaces59.ImplicationsOur finding that coastal environments are particularly beneficial adds to the body of evidence linking coastal environments with health and well-being and suggests this is consistent across many countries. Previous research has found that greater proximity to blue spaces, especially coastal settings, predicts visit frequency14,60,61 as well as other health outcomes—e.g. reduced risk of mortality and better general health, well-being and physical activity53,62. Here, we found that shorter travel times also predict visit well-being, highlighting the importance of having equitable access to good quality natural environments near to people’s homes.We also identified that different types of coastal and inland blue spaces (e.g. seaside promande, rural rivers), with different qualities (e.g. wildlife present), involving particular types of activities in specific social configurations (e.g. playing with children), were especially good at promoting well-being. This moves beyond a simple location-based assessment of benefit to one that recognises the complex interplay between location, behavioural and social processes. Numerous commentators63 (including Bratman et al.9 on which we have based this paper) have argued that we need to go beyond the determinate effects of green and blue spaces and develop a far richer, more nuanced understanding. The approach we have taken here is intended as a step in this direction.In terms of policy applications, these results provide support for the potential health benefits of efforts to improve equitable access to high quality environments, such as the English Coast Path (https://englandcoastpath.co.uk/) and the creation of beaches in Barcelona with the Olympic project in 199264. Our results also hint at the importance of high-level legislation, such as the EU’s Bathing Waters Directive65 for mental well-being37. If conducted on a more fine-grained geographical level, results could have the potential to leverage public support for more localised conservation initiatives. Furthermore, such results could be used as a basis for integration into more systematic conservation planning66.Further researchAlthough we incorporate a range of variables in our analysis, and our pseudo-R2 values are relatively high for a social research context, considerable variation remains unexplained. Although other individual characteristics may be important, such as nature connectedness67 and memories68, further research could explore the specific ecosystem features and social contexts associated with the particular positive results from coastal spaces, which would be of interest to policy makers and environmental managers. We also speculated that purpose of visit may explain some of our findings. Further research could explore the interactions between motivations and location, experience, and well-being outcomes.The presence of wildlife was differentially important across bluespace types and further research could unpack this. Exploring similar possibilities for the other quality metrics, as well as considering additional ecosystem characteristics, would also be informative. For example, identifying which factors are important in perceptions of safety in blue spaces. Bratman et al.9 also considered effect modification by visitor characteristics and further research could include interactions, or sub-group analysis, by socio-demographic factors.Further research could also explore longer-term benefits of these features over repeated visits; the potential for ecosystem disservices, such as the relationships we find between an interaction of wildlife and ice rinks and well-being; the potential for negative outcomes associated with ecosystem degradation69; and the potential for positive mental health outcomes from ecological restoration70.We have demonstrated some of the complexities involved in the human-nature relationship and that many factors are related to the outcome from a visit. The conceptual model applied allows the investigation of a wide range of variables including natural features and other environmental qualities, and characteristics of the exposure and experience, as well as individual parameters. We suggest that other researchers can apply this conceptual model and design data collection accordingly to target specific research questions and hypotheses (as opposed to where we have fitted already collected data). More

Jayathilake PG, Jana S, Rushton S, Swailes D, Bridgens B, Curtis T, et al. Extracellular polymeric substance production and aggregated bacteria colonization influence the competition of microbes in biofilms. Front Microbiol. 2017;8:1865.Article 
PubMed 
PubMed Central 

Google Scholar 
Flemming H-C, Neu TR, Wingender J (eds). The Perfect Slime: Microbial Extracellular Polymeric Substances (EPS). IWA Publishing, London, 2016.Morales-García AL, Bailey RG, Jana S, Burgess JG. The role of polymers in cross-kingdom bioadhesion. Philos Trans R Soc Lond B Biol Sci. 2019;374:20190192.Article 
PubMed 
PubMed Central 

Google Scholar 
Davey ME, O’Toole GA. Microbial biofilms: From ecology to molecular genetics. Microbiol Mol. 2000;64:847–67.Article 
CAS 

Google Scholar 
Peters BM, Jabra-Rizk MA, O’May GA, Costerton JW, Shirtliff ME. Polymicrobial interactions: Impact on pathogenesis and human disease. Clin Microbiol Rev. 2012;25:193–213.Article 
PubMed 
PubMed Central 

Google Scholar 
Karygianni L, Ren Z, Koo H, Thurnheer T. Biofilm matrixome: Extracellular components in structured microbial communities. Trends Microbiol. 2020;28:668–81.Article 
CAS 
PubMed 

Google Scholar 
Morris BEL, Henneberger R, Huber H, Moissl-Eichinger C. Microbial syntrophy: Interaction for the common good. FEMS Microbiol Rev. 2013;37:384–406.Article 
CAS 
PubMed 

Google Scholar 
Decho AW, Gutierrez T. Microbial extracellular polymeric substances (EPSs) in ocean systems. Front Microbiol. 2017;8:922.Article 
PubMed 
PubMed Central 

Google Scholar 
Boleij M, Seviour T, Wong LL, van Loosdrecht MCM, Lin Y. Solubilization and characterization of extracellular proteins from anammox granular sludge. Water Res. 2019;164:114952.Article 
CAS 
PubMed 

Google Scholar 
Kim D, Barraza JP, Arthur RA, Hara A, Lewis K, Liu Y, et al. Spatial mapping of polymicrobial communities reveals a precise biogeography associated with human dental caries. Proc Natl Acad Sci USA. 2020;117:12375–86.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Sadiq FA, Burmølle M, Heyndrickx M, Flint S, Lu W, Chen W, et al. Community-wide changes reflecting bacterial interspecific interactions in multispecies biofilms. Crit Rev Microbiol. 2021;47:338–58.Article 
PubMed 

Google Scholar 
Liu W, Jacquiod S, Brejnrod A, Russel J, Burmølle M, Sørensen SJ. Deciphering links between bacterial interactions and spatial organization in multispecies biofilms. ISME J. 2019;13:3054–66.Article 
PubMed 
PubMed Central 

Google Scholar 
Bridier A, Le Coq D, Dubois-Brissonnet F, Thomas V, Aymerich S, Briandet R. The spatial architecture of Bacillus subtilis biofilms deciphered using a surface-associated model and in situ imaging. PLOS One. 2011;6:e16177.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Lee KWK, Periasamy S, Mukherjee M, Xie C, Kjelleberg S, Rice SA. Biofilm development and enhanced stress resistance of a model, mixed-species community biofilm. ISME J. 2014;8:894–907.Article 
CAS 
PubMed 

Google Scholar 
Myszka K, Czaczyk K. Characterization of adhesive exopolysaccharide (EPS) produced by Pseudomonas aeruginosa under starvation conditions. Curr Microbiol. 2009;58:541–6.Article 
CAS 
PubMed 

Google Scholar 
Harimawan A, Ting YP. Investigation of extracellular polymeric substances (EPS) properties of P. aeruginosa and B. subtilis and their role in bacterial adhesion. Colloids Surf B. 2016;146:459–67.Article 
CAS 

Google Scholar 
Yang X-R, Li H, Nie S-A, Su J-Q, Weng B-S, Zhu G-B, et al. Potential contribution of anammox to nitrogen loss from paddy soils in Southern China. Appl Environ Microbiol. 2015;81:938–47.Article 
PubMed 
PubMed Central 

Google Scholar 
Kartal B, van Niftrik L, Keltjens JT, Op den Camp HJM, Jetten MSM Chapter 3 – Anammox—Growth physiology, cell biology, and metabolism. In: Poole RK, editor. Adv Microb Physiol. 60: Academic Press; 2012. p. 211–62.Lu Y, Natarajan G, Nguyen TQN, Thi SS, Arumugam K, Seviour TW, et al. Species level enrichment of AnAOB and associated growth morphology under the effect of key metabolites. bioRxiv. 2020. 2020.02.04.934877Gonzalez-Gil G, Sougrat R, Behzad AR, Lens PN, Saikaly PE. Microbial community composition and ultrastructure of granules from a full-scale anammox reactor. Micro Ecol. 2015;70:118–31.Article 
CAS 

Google Scholar 
Kindaichi T, Yuri S, Ozaki N, Ohashi A. Ecophysiological role and function of uncultured Chloroflexi in an anammox reactor. Water Sci Technol. 2012;66:2556–61.Article 
CAS 
PubMed 

Google Scholar 
Qin Y, Han B, Cao Y, Wang T. Impact of substrate concentration on anammox-UBF reactors start-up. Bioresour Technol. 2017;239:422–9.Article 
CAS 
PubMed 

Google Scholar 
Chen Z, Meng Y, Sheng B, Zhou Z, Jin C, Meng F. Linking exoproteome function and structure to anammox biofilm development. Environ Sci Technol. 2019;53:1490–500.Article 
CAS 
PubMed 

Google Scholar 
Ali M, Shaw DR, Albertsen M, Saikaly PE. Comparative genome-centric analysis of freshwater and marine ANAMMOX cultures suggests functional redundancy in nitrogen removal processes. Front Microbiol. 2020;11:1637.Article 
PubMed 
PubMed Central 

Google Scholar 
Jia F, Yang Q, Liu X, Li X, Li B, Zhang L, et al. Stratification of extracellular polymeric substances (EPS) for aggregated anammox microorganisms. Environ Sci Technol. 2017;51:3260–8.Article 
CAS 
PubMed 

Google Scholar 
Hou X, Liu S, Zhang Z. Role of extracellular polymeric substance in determining the high aggregation ability of anammox sludge. Water Res. 2015;75:51–62.Article 
CAS 
PubMed 

Google Scholar 
Feng C, Lotti T, Lin Y, Malpei F. Extracellular polymeric substances extraction and recovery from anammox granules: Evaluation of methods and protocol development. Chem Eng J. 2019;374:112–22.Article 
CAS 

Google Scholar 
Lotti T, Carretti E, Berti D, Montis C, Del Buffa S, Lubello C, et al. Hydrogels formed by anammox extracellular polymeric substances: Structural and mechanical insights. Sci Rep. 2019;9:11633.Article 
PubMed 
PubMed Central 

Google Scholar 
Jakubovics NS, Goodman SD, Mashburn-Warren L, Stafford GP, Cieplik F. The dental plaque biofilm matrix. Periodontology 2000. 2021;86:32–56.Article 
PubMed 
PubMed Central 

Google Scholar 
Ðapa T, Leuzzi R, Ng YK, Baban ST, Adamo R, Kuehne SA, et al. Multiple factors modulate biofilm formation by the anaerobic pathogen Clostridium difficile. J Bacteriol. 2013;195:545–55.Article 
PubMed 
PubMed Central 

Google Scholar 
Honma K, Inagaki S, Okuda K, Kuramitsu HK, Sharma A. Role of a Tannerella forsythia exopolysaccharide synthesis operon in biofilm development. Micro Pathog. 2007;42:156–66.Article 
CAS 

Google Scholar 
Li X-R, Du B, Fu H-X, Wang R-F, Shi J-H, Wang Y, et al. The bacterial diversity in an anaerobic ammonium-oxidizing (anammox) reactor community. Syst Appl Microbiol. 2009;32:278–89.Article 
CAS 
PubMed 

Google Scholar 
Cho S, Takahashi Y, Fujii N, Yamada Y, Satoh H, Okabe S. Nitrogen removal performance and microbial community analysis of an anaerobic up-flow granular bed anammox reactor. Chemosphere 2010;78:1129–35.Article 
CAS 
PubMed 

Google Scholar 
Morgenroth E, Sherden T, Van Loosdrecht MCM, Heijnen JJ, Wilderer PA. Aerobic granular sludge in a sequencing batch reactor. Water Res. 1997;31:3191–4.Article 
CAS 

Google Scholar 
Wong LL, Natarajan G, Boleij M, Thi SS, Winnerdy FR, Mugunthan S, et al. Extracellular protein isolation from the matrix of anammox biofilm using ionic liquid extraction. Appl Microbiol Biotechnol. 2020;104:3643–54.Article 
CAS 
PubMed 

Google Scholar 
Law Y, Kirkegaard RH, Cokro AA, Liu X, Arumugam K, Xie C, et al. Integrative microbial community analysis reveals full-scale enhanced biological phosphorus removal under tropical conditions. Sci Rep. 2016;6:25719.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ondov BD, Bergman NH, Phillippy AM. Interactive metagenomic visualization in a Web browser. BMC Bioinform. 2011;12:385.Article 

Google Scholar 
Liu X, Arumugam K, Natarajan G, Seviour TW, Drautz-Moses DI, Wuertz S, et al. Draft genome sequence of a Candidatus brocadia bacterium enriched from activated sludge collected in a tropical climate. Genome Announc. 2018;6:e00406–18.Article 
PubMed 
PubMed Central 

Google Scholar 
Price MN, Dehal PS, Arkin AP. FastTree 2-approximately maximum-likelihood trees for large alignments. PloS One. 2010;5:e9490–e.Article 
PubMed 
PubMed Central 

Google Scholar 
Arkin AP, Cottingham RW, Henry CS, Harris NL, Stevens RL, Maslov S, et al. KBase: The United States Department of Energy Systems Biology Knowledgebase. Nat Biotechnol. 2018;36:566–9.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Buchfink B, Reuter K, Drost H-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat Methods. 2021;18:366–8.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.Article 
CAS 
PubMed 

Google Scholar 
Daims H, Nielsen JL, Nielsen PH, Schleifer K-H, Wagner M. In situ characterization of Nitrospira-like nitrite-oxidizing bacteria active in wastewater treatment plants. Appl Environ Microbiol. 2001;67:5273–84.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Seviour T, Wong LL, Lu Y, Mugunthan S, Yang Q, Shankari UDOCS, et al. Phase transitions by an abundant protein in the anammox extracellular matrix mediate cell-to-cell aggregation and biofilm formation. mBio 2020;11:e02052–20.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Protter DSW, Rao BS, Van Treeck B, Lin Y, Mizoue L, Rosen MK, et al. Intrinsically disordered regions can contribute promiscuous interactions to RNP granule assembly. Cell Rep. 2018;22:1401–12.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Fulton KM, Smith JC, Twine SM. Clinical applications of bacterial glycoproteins. Expert Rev Proteom. 2016;13:345–53.Article 
CAS 

Google Scholar 
Upreti RK, Kumar M, Shankar V. Bacterial glycoproteins: Functions, biosynthesis and applications. Proteomics 2003;3:363–79.Article 
CAS 
PubMed 

Google Scholar 
van Teeseling MCF, Maresch D, Rath CB, Figl R, Altmann F, Jetten MSM, et al. The S-layer protein of the anammox bacterium Kuenenia stuttgartiensiss is heavily O-glycosylated. Front Microbiol. 2016;7:1721.PubMed 
PubMed Central 

Google Scholar 
McGonigle JM, Lang SQ, Brazelton WJ, Parales RE. Genomic evidence for formate metabolism by Chloroflexi as the key to unlocking deep carbon in lost city microbial ecosystems. Appl Environ Microbiol. 2020;86:e02583–19.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Vuillemin A, Kerrigan Z, D’Hondt S, Orsi WD. Exploring the abundance, metabolic potential and gene expression of subseafloor Chloroflexi in million-year-old oxic and anoxic abyssal clay. FEMS Microbiol Ecol. 2020;96:fiaa223.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kartal B, de Almeida NM, Maalcke WJ, Op den Camp HJ, Jetten MS, Keltjens JT. How to make a living from anaerobic ammonium oxidation. FEMS Microbiol Rev. 2013;37:428–61.Article 
CAS 
PubMed 

Google Scholar 
Loera-Muro A, Guerrero-Barrera A, Tremblay DNY, Hathroubi S, Angulo C. Bacterial biofilm-derived antigens: A new strategy for vaccine development against infectious diseases. Expert Rev Vaccines. 2021;20:385–96.Article 
CAS 
PubMed 

Google Scholar 
Hobley L, Harkins C, MacPhee CE, Stanley-Wall NR. Giving structure to the biofilm matrix: An overview of individual strategies and emerging common themes. FEMS Microbiol Rev. 2015;39:649–69.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Elias S, Banin E. Multi-species biofilms: Living with friendly neighbors. FEMS Microbiol Rev. 2012;36:990–1004.Article 
CAS 
PubMed 

Google Scholar 
Teeseling MCFV, Almeida NMD, Klingl A, Speth DR, Camp HJMOD, Rachel R, et al. A new addition to the cell plan of anammox bacteria: Candidatus Kuenenia stuttgartiensis has a protein surface layer as the outermost layer of the cell. J Bacteriol. 2014;196:80–9.Article 
PubMed 
PubMed Central 

Google Scholar 
Paula AJ, Hwang G, Koo H. Dynamics of bacterial population growth in biofilms resemble spatial and structural aspects of urbanization. Nat Commun. 2020;11:1354.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kragelund C, Caterina L, Borger A, Thelen K, Eikelboom D, Tandoi V, et al. Identity, abundance and ecophysiology of filamentous Chloroflexi species present in activated sludge treatment plants. FEMS Microbiol Ecol. 2007;59:671–82.Article 
CAS 
PubMed 

Google Scholar 
Nierychlo M, Miłobędzka A, Petriglieri F, McIlroy B, Nielsen PH, McIlroy SJ. The morphology and metabolic potential of the Chloroflexi in full-scale activated sludge wastewater treatment plants. FEMS Microbiol Ecol. 2018;95.Kragelund C, Thomsen TR, Mielczarek AT, Nielsen PH. Eikelboom’s morphotype 0803 in activated sludge belongs to the genus Caldilinea in the phylum Chloroflexi. FEMS Microbiol Ecol. 2011;76:451–62.Article 
CAS 
PubMed 

Google Scholar 
Zhang J, Miao Y, Zhang Q, Sun Y, Wu L, Peng Y. Mechanism of stable sewage nitrogen removal in a partial nitrification-anammox biofilm system at low temperatures: Microbial community and EPS analysis. Bioresour Technol. 2020;297:122459.Article 
CAS 
PubMed 

Google Scholar 
Björnsson L, Hugenholtz P, Tyson GW, Blackall LL. Filamentous Chloroflexi (green non-sulfur bacteria) are abundant in wastewater treatment processes with biological nutrient removal. Microbiology 2002;148:2309–18.Article 
PubMed 

Google Scholar 
Boleij M, Pabst M, Neu TR, van Loosdrecht MCM, Lin Y. Identification of glycoproteins isolated from extracellular polymeric substances of full-scale anammox granular sludge. Environ Sci Technol. 2018;52:13127–35.Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pabst M, Grouzdev DS, Lawson CE, Kleikamp HBC, de Ram C, Louwen R, et al. A general approach to explore prokaryotic protein glycosylation reveals the unique surface layer modulation of an anammox bacterium. ISME J. 2022;16:346–57.Article 
CAS 
PubMed 

Google Scholar 
Berlanga M, Guerrero R. Living together in biofilms: The microbial cell factory and its biotechnological implications. Micro Cell Fact. 2016;15:165.Article 

Google Scholar 
Liu T, Tian R, Li Q, Wu N, Quan X. Strengthened attachment of anammox bacteria on iron-based modified carrier and its effects on anammox performance in integrated floating-film activated sludge (IFFAS) process. Sci Total Environ. 2021;787:147679.Article 
CAS 
PubMed 

Google Scholar  More

Zhe, W. et al. Probing the nature of soil organic matter. Crit. Rev. Environ. Sci. Technol. 52, 4072–4093 (2022).Article 

Google Scholar 
Blanco-Canqui, H. & Lal, R. Mechanisms of carbon sequestration in soil aggregates. Crit. Rev. Plant Sci. 23, 481–504 (2004).Article 
CAS 

Google Scholar 
Six, J., Paustian, K., Elliott, E. T. & Combrink, C. Soil structure and organic matter: I. Distribution of aggregate–size classes and aggregate–associated carbon. Soil Sci. Soc. Am. J. 64, 681–689 (2000).Article 
ADS 
CAS 

Google Scholar 
Lehmann, J. & Kleber, M. The contentious nature of soil organic matter. Nature 528, 60–68 (2015).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Six, J., Bossuyt, H., Degryze, S. & Denef, K. A history of research on the link between (micro)aggregates, soil biota, and soil organic matter dynamics. Soil Tillage Res. 79, 7–31 (2004).Article 

Google Scholar 
Tisdall, J. M. & Oades, J. M. Organic matter and water-stable aggregates in soils. Eur. J. Soil Sci. 33, 141–163 (1982).Article 
CAS 

Google Scholar 
Luo, Y. et al. Rice rhizodeposition promotes the build-up of organic carbon in soil via fungal necromass. Soil Biol. Biochem. 160, 108345 (2021).Article 
CAS 

Google Scholar 
Wang, X. et al. Organic amendments drive shifts in microbial community structure and keystone taxa which increase C mineralization across aggregate size classes. Soil Biol. Biochem. 153, 108062 (2021).Article 
CAS 

Google Scholar 
Duan, Y. et al. Long–term fertilisation reveals close associations between soil organic carbon composition and microbial traits at aggregate scales. Agric. Ecosyst. Environ. 306, 107169 (2021).Article 
CAS 

Google Scholar 
Christensen, B. T. Physical fractionation of soil and structural and functional complexity in organic matter turnover. Eur. J. Soil Sci. 52, 345–353 (2001).Article 
CAS 

Google Scholar 
Olk, D. C. & Gregorich, E. G. Overview of the symposium proceedings, “meaningful pools in determining soil carbon and nitrogen dynamics”. Soil Sci. Soc. Am. J. 70, 967–974 (2006).Article 
ADS 
CAS 

Google Scholar 
Courtier-Murias, D. et al. Unraveling the long–term stabilization mechanisms of organic materials in soils by physical fractionation and NMR spectroscopy. Agric. Ecosyst. Environ. 171, 9–18 (2013).Article 
CAS 

Google Scholar 
Rodrigues, L. A. T. et al. Short– and long–term effects of animal manures and mineral fertilizer on carbon stocks in subtropical soil under no–tillage. Geoderma 386, 114913 (2021).Article 
ADS 
CAS 

Google Scholar 
Mao, J., Dan, C. O., Fang, X., He, Z. & Schmidt-Rohr, K. Influence of animal manure application on the chemical structures of soil organic matter as investigated by advanced solid–state NMR and FT–IR spectroscopy. Geoderma 146, 353–362 (2008).Article 
ADS 
CAS 

Google Scholar 
Simonetti, G. et al. Characterization of humic carbon in soil aggregates in a long–term experiment with manure and mineral fertilization. Soil Sci. Soc. Am. J. 25, 880–890 (2012).Article 

Google Scholar 
Cambardella, C. A. & Elliott, E. T. Particulate soil organic-matter changes across a grassland cultivation sequence. Soil Sci. Soc. Am. J. 56, 777–783 (1992).Article 
ADS 

Google Scholar 
Conant, R. T., Six, J. & Paustian, K. Land use effects on soil carbon fractions in the southeastern United States. I. Management-intensive versus extensive grazing. Biol. Fertil. Soils 38, 386–392 (2003).Article 
CAS 

Google Scholar 
Blanco-Moure, N., Gracia, R., Bielsa, A. C. & López, M. V. Soil organic matter fractions as affected by tillage and soil texture under semiarid Mediterranean conditions. Soil Tillage Res. 155, 381–389 (2016).Article 

Google Scholar 
Yu, H. et al. Accumulation of organic C components in soil and aggregates. Sci. Rep. 5, 13804 (2015).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Schöning, I., Morgenroth, G. & Kögel-Knabner, I. O/N–alkyl and alkyl C are stabilised in fine particle size fractions of forest soils. Biogeochemistry 73, 475–497 (2005).Article 

Google Scholar 
Solomon, D., Lehmann, J., Kinyangi, J., Liang, B. & Schäfer, T. Carbon K-edge NEXAFS and FTIR–ATR spectroscopic investigation of organic carbon speciation in soils. Soil Sci. Soc. Am. J. 13, 107–119 (2005).Article 

Google Scholar 
Yan, H., Chen, C., Xu, Z., Williams, D. & Xu, J. Assessing management impacts on soil organic matter quality in subtropical Australian forests using physical and chemical fractionation as well as 13C NMR spectroscopy. Soil Biol. Biochem. 41, 640–650 (2009).Article 

Google Scholar 
Masoom, H. et al. Soil organic matter in its native state: Unravelling the most complex biomaterial on earth. Environ. Sci. Technol. 50, 1670–1680 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Vogel, C. et al. Clay mineral composition modifies decomposition and sequestration of organic carbon and nitrogen in fine soil fractions. Biol. Fertil. Soils 51, 427–442 (2015).Article 
CAS 

Google Scholar 
Sharma, S., Singh, P., Angmo, P. & Satpute, S. Total and labile pools of organic carbon in relation to soil biological properties under contrasting land-use systems in a dry mountainous region. Carbon Manage. 13, 352–371 (2022).Article 
CAS 

Google Scholar 
Six, J., Elliott, E., Paustian, K. & Doran, J. Aggregation and soil organic matter accumulation in cultivated and native grassland soils. Soil Sci. Soc. Am. J. 62, 1367–1377 (1998).Article 
ADS 
CAS 

Google Scholar 
Elliott, E. T. Aggregate structure and carbon, nitrogen, and phosphorus in native and cultivated soils. Soil Sci. Soc. Am. J. 50, 627–633 (1986).Article 
ADS 

Google Scholar 
Yu, H., Ding, W., Luo, J., Geng, R. & Cai, Z. Long-term application of organic manure and mineral fertilizers on aggregation and aggregate-associated carbon in a sandy loam soil. Soil Tillage Res. 124, 170–177 (2012).Article 

Google Scholar 
Carter, M. R. & Gregorich, E. G. (eds) Soil Sampling and Methods of Analysis 2nd edn, 230–233 (Taylor & Francis Group, CRC, 2007).
Google Scholar 
Lu, R. (ed.) Soil and Agro-chemistry Analytical Methods 146–149 (China Agricultural Science and Technology Press, 1999).
Google Scholar 
Wu, J., Joergensen, R. G., Pommerening, B., Chaussod, R. & Brookes, P. C. Measurement of soil microbial biomass C by fumigation extraction: An automated procedure. Soil Biol. Biochem. 22, 1167–1169 (1990).Article 
CAS 

Google Scholar 
Zhang, X., Zhu, A., Yang, W. & Zhang, J. Accumulation of organic components and its association with macroaggregation in a sandy loam soil following conservation tillage. Plant Soil. 416, 1–15 (2017).Article 
CAS 

Google Scholar 
Skjemstad, J. O., Clarke, P., Taylor, J. A., Oades, J. M. & Newman, R. H. The removal of magnetic materials from surface soils—a solid state 13C CP/MAS NMR study. Soil Res. 32, 1215–1229 (1994).Article 
CAS 

Google Scholar 
Ringle, C. M., Wende, S. & Becker, J. M. SmartPLS 3.” Boenningstedt: SmartPLS GmbH. Preprint at http://www.smartpls.com (2015).Jerbi, M., Labidi, S., Lounès-Hadj Sahraoui, A., Chaar, H. & Ben Jeddi, F. Higher temperatures and lower annual rainfall do not restrict, directly or indirectly, the mycorrhizal colonization of barley (Hordeum vulgare L.) under rainfed conditions. PLoS ONE 15, e0241794 (2020).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cohen, J. Statistical power analysis for the behavioral sciences 2nd edn, 407–530 (Erlbaum Associates, Berlin, 1988).MATH 

Google Scholar 
Singh, P. & Benbi, D. K. Physical and chemical stabilization of soil organic matter in cropland ecosystems under rice–wheat, maize–wheat and cotton–wheat cropping systems in northwestern India. Carbon Manag. 12, 603–621 (2021).Article 
CAS 

Google Scholar 
Kiem, R. & Kögel-Knabner, I. Contribution of lignin and polysaccharides to the refractory carbon pool in C–depleted arable soils. Soil Biol. Biochem. 35, 101–118 (2003).Article 
CAS 

Google Scholar 
Lutzow, M. V. et al. Stabilization of organic matter in temperate soils: mechanisms and their relevance under different soil conditions: a review. Eur. J. Soil Sci. 57, 426–445 (2006).Article 

Google Scholar 
Yudina, A. V., Klyueva, V. V., Romanenko, K. A. & Fomin, D. S. Micro- within macro: How micro-aggregation shapes the soil pore space and water-stability. Geoderma 415, 115771 (2022).Article 
ADS 
CAS 

Google Scholar 
Tisdall, J. M., Smith, S. E. & Rengasamy, P. Aggregation of soil by fungal hyphae. Soil Res. 35, 55–60 (1997).Article 

Google Scholar 
Li, T. et al. Contrasting impacts of manure and inorganic fertilizer applications for nine years on soil organic carbon and its labile fractions in bulk soil and soil aggregates. CATENA 194, 104739 (2020).Article 
CAS 

Google Scholar 
Liang, Y. et al. Effect of chemical fertilizer and straw-derived organic amendments on continuous maize yield, soil carbon sequestration and soil quality in a Chinese Mollisol. Agric. Ecosyst. Environ. 314, 107403 (2021).Article 
CAS 

Google Scholar 
Liang, C., Kästner, M. & Joergensen, R. G. Microbial necromass on the rise: The growing focus on its role in soil organic matter development. Soil Biol. Biochem. 150, 108000 (2020).Article 
CAS 

Google Scholar 
Sharma, S., Singh, P. & Kumar, S. Responses of soil carbon pools, enzymatic activity, and crop yields to nitrogen and straw incorporation in a rice-wheat cropping system in North-Western India. Front. Sustain. Food Syst. 4, 532704 (2020).Article 

Google Scholar 
Puget, P., Chenu, C. & Balesdent, J. Dynamics of soil organic matter associated with particle–size fractions of water–stable aggregates. Eur. J. Soil Sci. 51, 595–605 (2000).Article 

Google Scholar  More

Finkelman, S., Navarro, S., Rindner, M. & Dias, R. Effect of low pressure on the survival of Trogoderma granarium Everts, Lasioderma serricorne (F.) and Oryzaephilus surinamensis (L.) at 30°C. J. Stored. Prod. Res. 42, 23–30 (2006).Article 

Google Scholar 
Hosseininaveh, V. A., Bandani, A. P., Azmayeshfard, P. S., Hosseinkhani, S. & Kazzazi, M. Digestive proteolytic and amylolytic activities in Trogoderma granarium Everts (Dermestidae: Coleoptera). J. Stored. Prod. Res. 43, 515–522 (2007).Article 
CAS 

Google Scholar 
Burges, H. D. Development of the khapra beetle, Trogoderma granarium, in the lower part of its temperature range. J. Stored. Prod. Res. 44, 32–35 (2008).Article 

Google Scholar 
Hagstrum D. W & Subramanyam, B. Stored-Product Insect Resource (AACC International, 2009).Beal, R. S. Synopsis of the economic species of Trogoderma occurring in the United States with description of a new species (Coleoptera: Dermestidae). Ann. Entomol. Soc. Am. 49, 559–566 (1956).Article 

Google Scholar 
Kerr, J. A. Khapra beetle returns. Pest Control 49(12), 24–25 (1984).
Google Scholar 
Sinha, R. N. & Utida, S. Climatic areas potentially vulnerable to stored product insects in Japan. Appl. Entomol. Zool. 2, 124–132 (1967).Article 

Google Scholar 
Banks, H. J. Distribution and establishment of Trogoderma granarium Everts (Coleoptera: Dermestidae): Climatic and other influences. J. Stored. Prod. Res. 13, 183–202 (1977).Article 

Google Scholar 
Kavallieratos, N. G., Athanassiou, C. G., Guedes, R. N. C., Drempela, J. D. & Boukouvala, M. C. Invader competition with local competitors: Displacement or coexistence among the invasive khapra beetle, Trogoderma granarium Everts (Coleoptera: Dermestidae), and two other major stored-grain beetles?. Front. Plant. Sci. 8, 1837 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Lampiri, E., Baliota, G. V., Morrison, W. M., Domingue, M. J. & Athanassiou, C. Comparative population growth of the khapra beetle (Coleoptera: Dermestidae) and the warehouse beetle (Coleoptera: Dermestidae) on wheat and rice. J. Econ. Entomol. 115, 344–352 (2021).Article 

Google Scholar 
Athanassiou, C. G., Phillips, T. W. & Wakil, W. Biology and control of the khapra beetle, Trogoderma granarium, a major quarantine threat to global food security. Ann. Rev. Entomol. 64, 131–148 (2019).Article 
CAS 

Google Scholar 
Stibick, J. New pest response guidelines: khapra beetle. APHIS– PPQ–Emergency and Domestic Programs. (U.S Department of Agriculture, 2009).Myers, S. W. & Hagstrum, D. W. Quarantine, In Stored stored product protection, (ed. Hagstrum D.W. Phillips T.W. & Cuperus G.) 297–304 (Kansas State University, 2012).Day, C. & White, B. Khapra beetle, Trogoderma granarium interceptions and eradications in Australia and around the world. In SARE working papers 1609. (Crawley: School of Agricul. Res. Econ. 2016).Burges, H. D. Diapause, pest status and control of the Khapra beetle. Trogoderma Granar. Everts Ann. Appl. Biol. 50, 614–617 (1962).Article 

Google Scholar 
Nair, K. & Desai, A. The termination of diapause in Trogoderma granarium Everts (Coleoptera, Dermestidae). J. Stored. Prod. Res. 8, 275–290 (1973).Article 

Google Scholar 
Burges, H. D. Studies on the Dermestid beetle Trogoderma granarium Everts—IV. Feeding, growth, and respiration with particular reference to diapause larvae. J. Insect. Physiol. 5, 317–334 (1960).Article 
CAS 

Google Scholar 
Wilches, D., Laird, R. A., Floate, K. & Fields, P. G. A review of diapause and tolerance to extreme temperatures in dermestids (Coleoptera). J. Stored Prod. Res. 68, 50–62 (2016).Article 

Google Scholar 
Vick, K. W., Drummond, P. C. & Coffelt, J. A. Trogoderma inclusum and T. glabrum: Effects of time of day on production of female pheromone, male responsiveness and mating. Ann. Entomol. Soc. Am. 66, 1001–1004 (1973).Article 

Google Scholar 
Partida, G. J. & Strong, R. G. Distribution and relative abundance of Trogoderma spp. in relation to climate zones of California. J. Econ. Entomol. 63, 1553–1560 (1970).Article 

Google Scholar 
Hagstrum, D. W. Seasonal variation of stored wheat environment and insect populations. J. Econ. Entomol. 16, 77–83 (1987).
Google Scholar 
Mullen, M. A. & Arbogast, R. T. Insect succession in a stored-corn ecosystem in southeast Georgia. J. Econ. Entomol. 81, 899–912 (1988).
Google Scholar 
Partida, G. J. & Strong, R. G. Comparative studies on the biologies of six species of Trogoderma: T. inclusum. Ann. Entomol. Soc. Am. 68, 91–103 (1975).Article 

Google Scholar 
Beal, R. S. Biology and taxonomy of the nearctic species of Trogoderma. Univ. Calif. Misc. Publ. Entomol. 10, 35–102 (1954).
Google Scholar 
Castañé, C., Agustí, N., del Estal, P. & Riudavets, J. Survey of Trogoderma spp in Spanish mills and warehouses. J. Stored. Prod. Res. 88, 1061 (2020).Article 

Google Scholar 
Levinson, H. Z. & Mori, K. The pheromone activity of chiral isomers of trogodermal for male khapra beetles. Naturwissenschaften 67, 148–149 (1980).Article 
CAS 

Google Scholar 
Silverstein, R. M. et al. Perception by Trogoderma species of chirality and methyl branching at a site far removed from a functional group in a pheromone component. J. Chem. Ecol. 6, 911–917 (1980).Article 
CAS 

Google Scholar 
Vick, K. W. Effects of interspecific matings of Trogoderma glabrum and T. inclusum on oviposition and re-mating. Ann. Entomol. Soc. Am. 66, 237–239 (1973).Article 
MathSciNet 

Google Scholar 
Drijfhout, S. et al. Catalogue of abrupt shifts in intergovernmental panel on climate change climate models. Proc. Natl. Acad. Sci. USA 112, E5777–E5786 (2015).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Phillips, T. W., Pfannenstiel, L. & Hagstrum, D. Survey of Trogoderma species (Coleoptera: Dermestidae) associated with international trade of dried distiller’s grains and solubles in the USA. Julius-Kühn-Archiv 1, 233–238 (2018).
Google Scholar 
Hadaway, A. The biology of the beetles, Trogoderma granarium Everts and Trogoderma versicolor (Creutz). Bull. Entomol. Res. 46, 781–796 (1956).Article 
CAS 

Google Scholar 
Gorham, J. R. Insect and Mite Pests in Food: An Illustrated Key. Vols. 1 and 2, (U.S Department of Agriculture, 1991).Furui, S., Miyanoshita, A., Imamura, T., Minegishi, Y. & Kokutani, R. Qualitative real-time PCR identification of the khapra beetle, Trogoderma granarium (Coleoptera: Dermestidae). Appl. Entomol. Zool. 54, 101–107 (2019).Article 
CAS 

Google Scholar 
Olson, R. L., Farris, R. E., Barr, N. B. & Cognato, A. I. Molecular identification of Trogoderma granarium (Coleoptera: Dermestidae) using the 16s gene. J Pest Sci 87, 701–710 (2014).Article 

Google Scholar 
Wu, Y. et al. Development of an array of molecular tools for the identification of khapra beetle (Trogoderma granarium), a destructive beetle of stored food products. Sci. Rep. 13, 3327 (2023).Article 
PubMed 
PubMed Central 

Google Scholar 
Lampiri, E., Athanassiou, C. & Arthur, F. H. Population growth and development of the khapra beetle (Coleoptera: Dermestidae), on different sorghum fractions. J. Econ. Entomol. 114, 424–429 (2021).Article 
CAS 
PubMed 

Google Scholar 
Athanassiou, C. G., Kavallieratos, N. G. & Boukouvala, M. C. Population growth of the khapra beetle, Trogoderma granarium Everts (Coleoptera: Dermestidae) on different commodities. J. Stored. Prod. Res. 69, 72–77 (2016).Article 

Google Scholar 
Karnavar, G. K. Mating behaviour and fecundity in Trogoderma granarium (Coleoptera: Dermestidae). J. Stored. Prod. Res. 8, 65–69 (1972).Article 

Google Scholar 
Pray, L. A. & Goodnight, C. J. Genetic variation in inbreeding depression in the red flour beetle Tribolium castaneum. Evolution 49, 176–188 (1995).Article 
PubMed 

Google Scholar 
Barzin, S., Naseri, B., Fathi, S. A. A., Razmjou, J. & Aeinehchi, P. Feeding efficiency and digestive physiology of Trogoderma granarium Everts (Coleoptera: Dermestidae) on different rice cultivars. J. Stored. Prod. Res. 84, 101511 (2019).Article 

Google Scholar 
Naseri, B., Aeinehchi, P. & Ashjerdi, A. R. Nutritional responses and digestive enzymatic profile of Trogoderma granarium Everts (Coleoptera: Dermestidae) on 10 commercial rice cultivars. J. Stored. Prod. Res. 87, 101591 (2020).Article 

Google Scholar 
Sarwar, M. & Sattar, M. Varietals assessment of different wheat varieties for their resistance response to Khapra beetle Trogoderma granarium. Pak. J. Seed. Technol. 1(10), 1–7 (2007).
Google Scholar 
Wilches, D., Laird, R., Floate, K. & Fields, P. Effects of acclimation and diapause on the cold tolerance of Trogoderma granarium. Entomol. Exp. Appl. 165, 169–178 (2017).Article 
CAS 

Google Scholar 
Paini, D. R. & Yemshanov, D. Modelling the arrival of invasive organisms via the international marine shipping network: a Khapra beetle study. PLoS ONE 7(9), e44589 (2012).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Morrison, W. R., Grosdidier, R. F., Arthur, F. H., Myers, S. W. & Domingue, M. J. Attraction, arrestment, and preference by immature Trogoderma variabile and Trogoderma granarium to food and pheromonal stimuli. J. Pest Sci. 93, 135–147 (2020).Article 

Google Scholar 
Arthur, F. H. & Morrison, W. M. Methodology for assessing progeny production and grain damage on commodities treated with insecticides. Agronomy 10(6), 804 (2020).Article 
CAS 

Google Scholar  More

Post, D. M. Using stable isotopes to estimate trophic position: Models, methods, and assumptions. Ecology 83, 703–718 (2002).Article 

Google Scholar 
Fry, B. Stable Isotope Ecology (Springer, 2007).
Google Scholar 
Boon, P. I. & Bunn, S. E. Variations in the stable isotope composition of aquatic plants and their implications for food web analysis. Aquat. Bot. 48, 99–108 (1994).Article 

Google Scholar 
Kling, G. W., Fry, B. & O’Brien, W. J. Stable isotopes and planktonic trophic structure in arctic lakes. Ecology 73, 561–566 (1992).Article 

Google Scholar 
Nielsen, J. M., Clare, E. L., Hayden, B., Brett, M. T. & Kratina, P. Diet tracing in ecology: Method comparison and selection. Methods Ecol. Evol. 9, 278–291 (2018).Article 

Google Scholar 
Coulter, A. A., Swanson, H. K. & Goforth, R. R. Seasonal variation in resource overlap of invasive and native fishes revealed by stable isotopes. Biol. Invasions 21, 315–321 (2019).Article 

Google Scholar 
Jung, A. S., Van Der Veer, H. W., Van Der Meer, M. T. & Philippart, C. J. Seasonal variation in the diet of estuarine bivalves. PLoS One 14, e0217003 (2019).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Devlin, S. P., Vander Zanden, M. J. & Vadeboncoeur, Y. Depth-specific variation in carbon isotopes demonstrates resource partitioning among the littoral zoobenthos. Freshw. Biol. 58, 2389–2400 (2013).CAS 

Google Scholar 
Possamai, B., Vieira, J. P., Grimm, A. M. & Garcia, A. M. Temporal variability (1997–2015) of trophic fish guilds and its relationships with El Niño events in a subtropical estuary. Estuar. Coast. Shelf Sci. 202, 145–154 (2018).Article 
ADS 

Google Scholar 
Syvaranta, J., Hamalainen, H. & Jones, R. I. Within-lake variability in carbon and nitrogen stable isotope signatures. Freshw. Biol. 51, 1090–1102 (2006).Article 
CAS 

Google Scholar 
Janbu, A. D., Paasche, Ø. & Talbot, M. R. Paleoclimate changes inferred from stable isotopes and magnetic properties of organic-rich lake sediments in Arctic Norway. J. Paleolimnol. 46, 29 (2011).Article 
ADS 

Google Scholar 
Leng, M. et al. Late quaternary palaeoenvironmental reconstruction from Lakes Ohrid and Prespa (Macedonia/Albania border) using stable isotopes. Biogeosciences 7, 3109–3122 (2010).Article 
ADS 
CAS 

Google Scholar 
Jiang, Q., Shen, J., Liu, X., Zhang, E. & Xiao, X. A high-resolution climatic change since holocene inferred from multi-proxy of lake sediment in westerly area of China. Chin. Sci. Bull. 52, 1970–1979 (2007).Article 

Google Scholar 
Finlay, J. C. & Kendall, C. Stable isotope tracing of temporal and spatial variability in organic matter sources to freshwater ecosystems. Stable Isot. Ecol. Environ. Sci. 2, 283–333 (2007).Article 

Google Scholar 
Harvey, C. J. & Kitchell, J. F. A stable isotope evaluation of the structure and spatial heterogeneity of a Lake Superior food web. Can. J. Fish. Aquat. Sci. 57, 1395–1403 (2000).Article 
CAS 

Google Scholar 
Xu, D. et al. Spatial heterogeneity of food web structure in a large shallow eutrophic lake (Lake Taihu, China): Implications for eutrophication process and management. J. Freshw. Ecol. 34, 229–245 (2019).Article 
CAS 

Google Scholar 
Ruokonen, T., Kiljunen, M., Karjalainen, J. & Hämäläinen, H. Invasive crayfish increase habitat connectivity: A case study in a large boreal lake. Knowl. Manag. Aquat. Ecosyst. https://doi.org/10.1051/kmae/2013034 (2012).Article 

Google Scholar 
Veselý, L. et al. The crayfish distribution, feeding plasticity, seasonal isotopic variation and trophic role across ontogeny and habitat in a canyon-shaped reservoir. Aquat. Ecol. 54, 1169–1183 (2020).Article 

Google Scholar 
Kalff, J. Limnology: Inland Water Ecosystems Vol. 592 (Prentice Hall, 2002).
Google Scholar 
Polačik, M., Harrod, C., Blažek, R. & Reichard, M. Trophic niche partitioning in communities of African annual fish: Evidence from stable isotopes. Hydrobiologia 721, 99–106 (2014).Article 

Google Scholar 
Costalago, D., Navarro, J., Álvarez-Calleja, I. & Palomera, I. Ontogenetic and seasonal changes in the feeding habits and trophic levels of two small pelagic fish species. Mar. Ecol. Prog. Ser. 460, 169–181 (2012).Article 
ADS 

Google Scholar 
Matthews, B. & Mazumder, A. Consequences of large temporal variability of zooplankton δ15N for modeling fish trophic position and variation. Limnol. Oceanogr. 50, 1404–1414 (2005).Article 
ADS 
CAS 

Google Scholar 
Taipale, S., Kankaala, P., Tiirola, M. & Jones, R. I. Whole-lake dissolved inorganic 13C additions reveal seasonal shifts in zooplankton diet. Ecology 89, 463–474 (2008).Article 
PubMed 

Google Scholar 
Zohary, T., Erez, J., Gophen, M., Berman-Frank, I. & Stiller, M. Seasonality of stable carbon isotopes within the pelagic food web of Lake Kinneret. Limnol. Oceanogr. 39, 1030–1043 (1994).Article 
ADS 
CAS 

Google Scholar 
Stenroth, P. et al. Stable isotopes as an indicator of diet in omnivorous crayfish (Pacifastacus leniusculus): The influence of tissue, sample treatment, and season. Can. J. Fish. Aquat. Sci. 63, 821–831 (2006).Article 
CAS 

Google Scholar 
R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/ (2021).Moore, J. W. & Semmens, B. X. Incorporating uncertainty and prior information into stable isotope mixing models. Ecol. Lett. 11, 470–480 (2008).Article 
PubMed 

Google Scholar 
Stock, B. C. & Semmens, B. X. Unifying error structures in commonly used biotracer mixing models. Ecology 97, 2562–2569 (2016).Article 
PubMed 

Google Scholar 
Irz, P., Laurent, A., Messad, S., Pronier, O. & Argillier, C. Influence of site characteristics on fish community patterns in French reservoirs. Ecol. Freshw. Fish 11, 123–136 (2002).Article 

Google Scholar 
Sutela, T., Aroviita, J. & Keto, A. Assessing ecological status of regulated lakes with littoral macrophyte, macroinvertebrate and fish assemblages. Ecol. Indic. 24, 185–192 (2013).Article 

Google Scholar 
Hunt, P. & Jones, J. The effect of water level fluctuations on a littoral fauna. J. Fish Biol. 4, 385–394 (1972).Article 

Google Scholar 
Kaster, J. & Jacobi, G. Benthic macroinvertebrates of a fluctuating reservoir. Freshw. Biol. 8, 283–290 (1978).Article 

Google Scholar 
Kraft, K. The effect of unnatural water level fluctuations on benthic invertebrates in Voyageurs National Park. Research⁄Resources Management Report MWR-12. US Department of the Interior, National Park Service. International Falls, Minnesota (1988).Glon, M., Larson, E. R. & Pangle, K. Comparison of 13C and 15N discrimination factors and turnover rates between congeneric crayfish Orconectes rusticus and O. virilis (Decapoda, Cambaridae). Hydrobiologia 768, 51–61 (2016).Article 
CAS 

Google Scholar 
Hesslein, R. H., Hallard, K. & Ramlal, P. Replacement of sulfur, carbon, and nitrogen in tissue of growing broad whitefish (Coregonus nasus) in response to a change in diet traced by δ34S, δ13C, and δ15N. Can. J. Fish. Aquat. Sci. 50, 2071–2076 (1993).Article 
CAS 

Google Scholar  More