Longitudinal pattern of resource utilization by aquatic consumers along a disturbed subtropical urban river: Estimating the relative contribution of resources with stable isotope analysis

Abstract The utilization of food resources by aquatic consumers reflects the structure and functioning of river food webs. In lotic water systems, where food availability and predator–prey relationships vary with gradient changes in physical conditions, understanding diet assimilation by local communities is important for ecosystem conservation. In the subtropical Liuxi River, southern China, the relative contribution of basal resources to the diet assimilation of functional feeding groups (FFGs) was determined by stable carbon (13C) and nitrogen (15N) isotope analyses. The output of Bayesian mixing models showed that diatom‐dominated periphyton (epilithic biofilm), aquatic C3 plants (submerged hydrophytes), and suspended particulate organic matter (SPOM) associated with terrestrial C3 plants contributed the most to the diet assimilation of FFGs in the upper, middle, and lower reaches, respectively. The relative contribution of consumer diet assimilation was weighted by the biomass (wet weight, g/m2) of each FFG to reflect resource utilization at the assemblage level. From the upper to the lower reaches, the spatial variation in the diet assimilation of fish and invertebrate assemblages could be summarized as a longitudinal decrease in periphyton (from 57%–76% to <3%) and an increase in SPOM (from <7% to 51%–68%) with a notable midstream increase in aquatic C3 plants (23%–48%). These results indicate that instream consumers in the Liuxi River rely more on autochthonous production (e.g., periphyton and submerged hydrophytes) than on terrestrially derived allochthonous matter (e.g., terrestrial plants). The pattern of resource utilization by consumers in the mid‐upper Liuxi River is consistent with findings from other open subtropical and neotropical rivers and provides evidence for the riverine productivity model. Our study indicates that protecting inherent producers in rivers (e.g., periphyton and submerged hydrophytes) and restoring their associated habitats (e.g., riffles with cobble substrate) are conducive to aquatic ecosystem management.


| INTRODUC TI ON
To manage the integrity and functioning of river ecosystems, it is important to recognize the region-specific characteristics of resource utilization by aquatic consumers (Wang et al., 2021;Welcomme et al., 2006;Zeni & Casatti, 2014). Several different theories have been developed to explain resource utilization in large river food webs and to provide a theoretical basis for the ecological functioning of large rivers, such as the river continuum concept (Vannote et al., 1980), the flood pulse concept (Junk, 1999), the riverine productivity model (Thorp & Delong, 1994), and riverine ecosystem synthesis . Each of these models focuses on a different dimension of the riverscape. The river continuum concept has a longitudinal perspective and describes ecosystem processes from the upstream to downstream reaches, while the flood pulse concept highlights the importance of energy transfer from lateral floodplains.
In contrast, the riverine productivity model suggests that carbon sources supporting food webs in large river systems derive from within the river channel itself. The riverine ecosystem synthesis emphasizes that although these resources are predominantly autochthonous in nature, they can be influenced by the ecological nature of river hydrogeomorphic units. Therefore, resource subsidies to river food webs should have local particularity and be explored over a range of climates, hydrologic regimes, watershed geochemistry conditions, and anthropogenic pressures. Ecologists examined the origin of potential resources that fuel riverine ecosystem processes, and the results emphasized the relative contribution of allochthonous versus autochthonous organic matter to consumers (Pingram et al., 2012a). For example, in forested headwater streams, allochthonous inputs such as leaves may form the base of the food web; however, in wider rivers with low elevations, autochthonous algae may be more important to consumers. Beyond the acknowledged longitudinal pattern in increased reliance on algae in lotic systems, tropical headwater streams have been shown to have greater reliance on algal resources than their temperate counterparts (Lewis et al., 2001;March & Pringle, 2003). Collectively, these studies demonstrated specific resource utilization in different habitats, regions, and climate zones.
Studies examining the spatial pattern of resource utilization in temperate streams, characterized by fishes and insects, support different concepts (Hoeinghaus et al., 2007). However, previous studies based on behavioral-morphologic feeding habits and stomach/gut content analysis (Aarts & Nienhuis, 2003) leave questions regarding which and how resources are assimilated by consumers and their relative contribution to consumers' diet composition unanswered (Buchheister & Latour, 2015). Particularly, along a longitudinal gradient, the changing physical structure of a river determines the variable distribution of resources; thus, it is hypothesized that covariation exists between resource utilization by consumers and environmental changes across space (Chang et al., 2012;Hoeinghaus et al., 2007). As riverine organisms display different food sources across habitats, especially migratory assemblages and generalist feeders (e.g., omnivores), it is challenging to precisely quantify resource utilization by region-or site-specific consumers.
Stable isotopes can be used to identify basal food resources used by consumers and determine broad-scale linkages between consumers and organic matter sources. Terrestrial and aquatic plants often differ in their 13 C signatures (Fry, 1991). Because the 13 C signatures of primary consumers reflect those of the plants they eat, analyses of 13 C signatures can be used to estimate the relative contributions of different basal resources in stream food webs (Finlay, 2001;March & Pringle, 2003;Rounick & Hicks, 1985).
Analysis of 13 C signatures provides an advantage over stomach/gut content analysis because the former measures the amount of carbon assimilated from each food source as opposed to the amount ingested. This differentiation is important because food sources can differ greatly in quality. Furthermore, stable isotopes are also useful in determining trophic levels because the 15 N signature is generally enriched by 2.5-3.5‰ with every trophic transfer (Neres-Lima et al., 2016;Peterson & Fry, 1987).
Due to spatially unbalanced development along disturbed river systems (e.g., upstream pristine forests, midstream agricultural areas, and downstream urbanized regions; Chang et al., 2012;Kaymak et al., 2018;Wang et al., 2018a), we hypothesized that the food resources used by local instream communities would be region-or site-specific. We examined the δ 13 C and δ 15 N values of basal resources and aquatic consumers in heterogeneous habitats along a subtropical river in southern China. A Bayesian mixing model was used to analyze the relative contributions of basal resources, including periphyton, aquatic C 3 plants (i.e., hydrophytes), terrestrial C 3 and C 4 plants (i.e., pasture grasses), and suspended particulate organic matter (SPOM), to the diet assimilation of fish and aquatic invertebrates. Differentiation between terrestrial C 3 and C 4 plants often indicates a catchment land-use transition from low intensity/fellow land cover (C 3 ) to high intensity/improved agriculture (C 4 ).
We aimed to test the hypothesis by addressing the following questions: (1) What are the relative contributions of various resources to aquatic consumers and assemblages? (2) Do the contributions of these resources vary along rivers whose natural continuum is disrupted by human activities? (3) Which concept best explains the spatial pattern of resource utilization along the studied subtropical river? By answering these questions, we can gain insight for evaluating the effects of ongoing management activities and potential restoration measures on stream or river ecosystems.

| Study area and sampling sites
The Liuxi River is regarded as the mother river (of high cultural significance) of Guangzhou, the capital city of Guangdong Province in southern China. The study area has a typical subtropical monsoon climate. The mean annual precipitation of the watershed is 1800 mm, mostly occurring in April-September. The Liuxi River watershed (with an area of 2300 km 2 ) is situated in the northeastern corner of the Pearl River Delta, which has experienced rapid development in the last two decades. The watershed spans four countylevel districts (Huadu, Luogang, Baiyun, and Conghua), which occupy 70% of the watershed area ( Figure 1). The water of the Liuxi River is used for a wide range of purposes, such as drinking water supply, agriculture and industry water use, and recreation. However, the source-water intake near the river outlet (site 8) was abandoned in 2010 due to pollution, and only intakes located in the middle river (sites 5-6) are still in use.
Eight sampling sites, extending from headwaters to the estuary, were chosen along the main channel of the Liuxi River ( Figure 1).
Environmental data, including habitat characteristics and physicochemical parameters of water quality (see Tables S1 and S2 of the Supporting Information), were provided by a nationally accredited (China Metrology Accreditation) third-party testing agency. Forests are the dominant land use type in the watershed, accounting for 55% of the total area. Approximately 33% of the watershed is used for agriculture, including orchards (20%), paddy rice fields (9%), and vegetable lands (4%). The rest of the watershed is occupied by built-up areas (9%) and water (3%). The development in the watershed is spatially unbalanced. The upstream areas (sites 1-2) are covered by dense forests, while agricultural activities are concentrated in the midstream areas (sites 3-5). The area close to the watershed outlet (sites 6-8) is highly urbanized.

| Sample collection and functional feeding group classification
Consumer and basal resource samples were collected concurrently from headwater site 1 to downstream site 8 in September 2020, which is the transition period between the rainy (June to July) and dry (November to December) seasons. The methods for collecting primary producers (e.g., periphyton, aquatic C 3 plants, and terrestrial C 3 and C 4 plants), aquatic consumers (fish, shrimp, crabs, insects, bivalves, gastropods, annelids, and zooplankton), and suspended particulate organic matter (SPOM) are provided in the supplementary methods of the Supporting Information. Macroinvertebrate and fish specimens were sampled and quantified in wet weight (g/ m 2 ) according to the methods for building mass-balance Ecopath models by Wang et al. (2018a). All consumer species were identified to the lowest possible taxonomic level and assigned to FFGs F I G U R E 1 Locations of the eight sampling sites (#1-#8) along the main channel of the Liuxi River following the classification of Morse et al. (1994) and Wang et al. (2019). Zooplankton specimens included cyclopoid copepods and malacostracan nauplii that could be selected while alive under stereomicroscopy. The dominant species/genera/families and FFGs of fish and invertebrates are provided in Table S3 of the Supporting Information.

| Preparation and determination for stable isotope analysis
According to previous basin-scale faunal surveys along the Liuxi River (data provided by Guangzhou Environmental Monitoring Center Station), the fish and invertebrate species in each FFG can represent the dominant populations of the local aquatic community (Wang et al., 2020c). Aquatic insects, tadpoles, snails, amphipods, and polychaetes were kept for 1 day in clean water to allow for gut evacuation. For small insect larvae, zooplankton, and amphipods, multiple individuals from the same taxon were pooled as replicates, and the entire composite was used. For large shrimp, crabs, snails, odonates, and fish, individual specimens were treated as replicates, and only the muscle tissues were used. All samples were freezedried, ground into homogeneous powder, and filtered through a 150-mesh sieve (pore size 100 µm). The sample size for stable isotope measurement of each FFG at each site is provided in Table S3 of the Supporting Information.
Samples of SPOM, zooplankton, bivalves, gastropods, and crustaceans were tested for the presence of carbonates by adding a few drops of 10% HCl to each subsample. Given that acidification can affect δ 15 N signals, a nonacidified subsample was always kept for separate nitrogen isotope analysis (Carabel et al., 2006). C:N ratios were used to normalize δ 13 C for lipid content. For the fish and inver- where δ 13 C normalized is the normalized estimate of δ 13 C for the effects of lipid concentration and δ 13 C untreated is the δ 13 C of samples determined without lipid extraction.
Stable isotope data were reported as the difference between the ratios of a sample and standard, using the standard notation of δ 13 C or δ 15 N (‰) = [(X sample /X standard ) − 1] × 10 3 , where X = 13 C/ 12 C or 15 N/ 14 N. δ 13 C or δ 15 N is the per-mil (‰) deviation of that sample from the recognized isotope standard: Vienna Pee Dee Belemnite for δ 13 C and atmospheric N 2 for δ 15 N. The δ 13 C and δ 15 N values of the samples were determined with a continuous-flow isotope ratio mass spectrometer (Finnigan Delta S) coupled with an elemental analyzer (NA 1500; Fison, Milan, Italy) in the Stable Isotope Laboratory, Third Institute of Oceanography, China. The precision of repeated measurements of the standard reference materials was ±0.1‰ for both the carbon and nitrogen stable isotopes. The standard error of the mean for replicates of the same tissue was 0.13‰ for δ 13 C and 0.15‰ for δ 15 N.

| Statistical analysis
A Bayesian mixing model (R package MixSIAR, R Core Team; Parnell et al., 2010) was used to estimate the relative contribution (%) of basal resources to the diet assimilation of invertebrate and fish FFGs based on stable isotope values. Prior to analysis, the δ 13 Cδ 15 N space of consumer FFGs and five basal resources at sampling sites 1-8 were plotted to diagnose the potential solutions (see the convex hull in Figure S1 of the Supporting Information). The δ 13 C and δ 15 N values of each consumer specimen were input as raw data, whereas those of basal resources were input as the means and standard deviations. Although the stable isotope data were correlated, MixSIAR assumes multivariate normality, which accounts for the fact that tracer values can covary. The model is fit via Markov chain Monte Carlo, which estimates entire posterior distributions for each variable. Aquatic C 3 plants were not used as a source in the MixSIAR model at site 1 because submerged hydrophytes did not appear in the headwaters. To back-calculate the underlying sources, we used a δ 13 C trophic fractionation of 0.96 ± 0.26‰ for all consumers and δ 15 N trophic fractionation of 1.07 ± 0.32‰ for invertebrates and 2.38 ± 0.37‰ for fish (Wang et al., 2020a).
To reflect the overall trend of resource utilization by local assemblages, the relative contributions of basal resources to the diet composition of fish and aquatic invertebrate assemblages were calculated as: where RC m is the relative contribution of resource m (e.g., periphyton) to the diet of a consumer assemblage (i.e., invertebrate or fish); B j is the biomass (fresh weight, g/m 2 ) of functional feeding group j where n is the highest number of local functional feeding groups; and R m is the relative contribution of resource m to the diet assimilation of functional feeding group j.
Nonmetric multidimensional scaling (NMDS) was used to visualize the relationships between consumers and their food resources based on the diet assimilation data matrix output by the MixSIAR model (i.e., relative contributions [%] of five resources to the diet assimilation of fish and invertebrate FFGs). NMDS is a nonparametric ordination technique that relies on the rank order of pairwise dissimilarities (Euclidean distance in this study) and does not make any underlying distributional assumptions about the data (Borcard et al., 2011). NMDS was chosen over parametric ordination approaches because the diet assimilations of consumer FFGs at sites 1-8 were skewed and not normally distributed. The sampling sites were plotted in the ordination space with the distance among points positively related to the dissimilarity of the consumer diet assimilation; thus,

| Relative contributions of basal resources to the diet assimilation of invertebrate FFGs
The diet assimilation of invertebrate FFGs showed two spatial patterns ( Figure 2), one along the longitudinal gradient (i.e., from the upper to the lower reaches) and the other across the lateral gradient (i.e., from the riparian or revetment zones to the central channel).
The shrimp scrapers (Atyidae: Neocaridina denticulata), which live near the littoral zone that has lush water grass, showed less longitudinal variation in diet assimilation. Thus, even if shrimp scrapers were sampled at downstream sites 7-8, they were grouped with the insect groups (e.g., insect c-g, c-f, scraper, and c-g/scraper) that were typi- assimilation of Cluster 4. Cluster 5 was an exception because it was composed of insect shredders that preferred to eat terrestrial C 3 and C 4 plant debris. In addition, insect c-gs at sites 5-6 could be sampled only in the littoral zone due to the high water depth, and decaying terrestrial plant debris might provide the dominant food sources for the local insect c-gs.

| Relative contributions of basal resources to the diet assimilation of fish FFGs
Similar to the pattern observed for invertebrate FFGs, fish FFGs were also divided into five clusters according to the relative contributions of five basal resources in the diet assimilation ( Figure 3).
Clusters 1, 2, and 4 were identified by the greatest contribution (in terms of the mean) of periphyton, aquatic C 3 plants, and SPOM to diet assimilation, respectively; Cluster 3 was identified by the mixed contribution of periphyton and aquatic C 3 plants to diet assimilation; and Cluster 5 was identified by the high contribution of terrestrial C 4 plants to diet assimilation. Compared with that of invertebrates, the longitudinal variation in resource utilization of fish showed a more obvious changing trend from periphyton to submerged hydrophytes and to SPOM along the river, suggesting that the diet assimilation of fish might reflect the long-term status of the local resource supply.
Resource utilization by fish assemblages showed significant differences among FFGs and exhibited three spatial patterns from the headwaters to the lower reaches: decreased periphyton with increased SPOM in the diet assimilation of periphytivorous, omnivorous, and molluscivorous fish; decreased periphyton with increased aquatic C 3 plants in the diet assimilation of insectivorous and invertivorous fish; and decreased aquatic and terrestrial C 3 plants with increased SPOM in the diet assimilation of crustaceavorous, piscivorous, detritivorous, phytoplanktivorous, and zooplanktivorous fish.
These results indicated the dependence of fish FFGs on local food resources. Considering the longitudinal shift in the resources assimilated by invertebrate FFGs, in addition to local basal resources, the diet assimilation of fish may also be influenced by their consumption of invertebrates with different isotope signatures.

| Comparison of diet assimilation among the site-specific consumer assemblages
The NMDS ordination ranked consumer assemblages composed of fish and invertebrate FFGs based on the relative contributions (%) of five resources to diet assimilation (Figure 4). The NMDS stress value of 0.06 shows a good fit. From site 1 to site 8, there was a longitudinal shift in dominant resource utilization by site-specific consumer assemblages, which was reflected by the clockwise distribution of five resources from the 3rd to the 2nd quadrant in the NMDS biplot.
In the upper stream (sites 1-2), periphyton, located near the center left of the biplot composed of the NMDS1 and NMDS2 axes, was an important resource used by local consumers. In the mid-upper reaches (sites 3-5), aquatic C 3 plants could be regarded as an important resource because they were close to the center of the 1st quadrant in the NMDS biplot. In the mid-lower reaches (sites 6-8), terrestrial C 3 plants, located near the center right of the NMDS biplot, were dominant resources used by local consumers. In addition to terrestrial C 3 plants, SPOM was an important resource for the consumers at estuarial site 8. Notably, the increased overlap area of a convex hull and Bayesian ellipses between sampling sites indicated a high degree of commonality in resource utilization by the consumers in the lower reaches.

| Evaluating the effects of environmental variables on system attributes
The accepted SEM explained a moderate amount of the variance in richness (60%) and fit the data well, with an RMSEA of 0.061 (p = .056) and a goodness-of-fit index (GFI) value of 0.94 ( Figure 6).
These index values showed that there was a strong correlation be-

| Region-specific resource utilization in the upper stream
Terrestrial plants are rarely limited by the CO 2 supply, and thus have δ 13 C fractionation consistent with that of atmospheric CO 2 , resulting in relatively constant δ 13 C values of −28.0 ± 0.9‰ for C 3 plants and −13.5 ± 1.5‰ for C 4 plants (Finlay, 2001). This result contrasts with the periphyton δ 13 C values of −18.5 ± 2.7‰ (commonly −30‰ to −15‰), which vary as a function of productivity and flow velocity (Finlay, 2004). In addition, periphyton δ 13 C values might also be influenced by periphyton composition, that is, seasonally dependent taxonomic composition of bacteria and algae, which can store large amounts of carbon and nutrients (Wu, 2016). Aquatic (e.g., periphyton) and terrestrial vegetation (e.g., vascular C 3 plants) form the diets of invertebrates or fish in streams with different climates, producing local coupling between consumers and local resources (Jardine & Bunn, 2012;March & Pringle, 2003;Pingram et al., 2012a). Diatomdominated periphyton (epilithic biofilm) appeared to be the main autochthonous carbon source driving secondary production in the upper wadable stream (sites 1-3) of this subtropical urban river in the Northern Hemisphere. Previous investigations in tropical systems also reported that periphyton, despite forming an inconspicuous component of the total carbon biomass, is largely responsible for metazoan biomass (Bunn et al., 2003;Thorp & Delong, 2002;Wang et al., 2018b), whereas detritus from terrestrial vascular C 3 and C 4 plants plays only a minor role (Lau et al., 2009;Lewis et al., 2001).
Our finding is consistent with those of studies on large rivers in Australia (Bunn et al., 2003;Medeiros & Arthington, 2011), the subtropics (e.g., the East River in southern China, Wang et al., 2020b), and the neotropics (e.g., the Parana´ and Orinoco Rivers in South America, Hoeinghaus et al., 2007;Jepsen & Winemiller, 2007), where benthic algae provide important carbon sources for food webs in lower order streams. In contrast, some large continental rivers (e.g., F I G U R E 3 Cluster analysis based on the relative contributions (%) of five basal resources to the diet assimilation of fish functional feeding groups (FFGs). SPOM, suspended particulate organic matter. Per, periphytivore; Ins, insectivore; Omn, omnivore; Det, detritivore; Inv, invertivore; Mol, molluscivore; Cru, crustaceavore; Pis, piscivore; Phy, phytoplanktivore; Zoo, zooplanktivore; Her, herbivore; Ann, annelidivore the Mississippi and Ohio Rivers in North America, Herwig et al., 2004;Thorp et al., 1998) in temperate and dry climates rely more on the terrestrial carbon supply that provides inputs as coarse organic particulates. In temperate streams, the low temperature, seasonal dryness, and dense canopy limit autochthonous production; thus, most primary consumers use allochthonous matter (e.g., riparian leaf litter, see Chang et al., 2012;Kaymak et al., 2018). In the tropics and subtropics, however, constant sunlight, an open canopy, perennial precipitation, and a wet climate stimulated the productivity of epilithic diatoms (e.g., Melosira and Aulacoseira, Wang et al., 2018a). Notably, the canopy immediately above the headwater streams was not covered by dense forests, because the riparian vegetation was composed of straight-growing eucalyptus and bamboo. As a result, diatom-dominated periphyton was the major carbon source for consumers in the upper Liuxi River, which explained the high contribution of periphyton to the diet assimilation of consumer FFGs at sites 1-3.

| Longitudinal shifts in resource utilization by consumer FFGs along the stream
An interesting finding was that submerged hydrophytes (e.g., H. verticillata and M. verticillatum) played an important role as a resource supply in the transition zone connecting the upper wadable and middle nonwadable stream sections (Figure 7a), which enhanced the effects of autochthonous resources on the energy supply of the local food web. Similarly, Hoeinghaus et al. (2007) examined carbon flow in 10 food webs of the Paraná River, Brazil, and found that aquatic C 3 plants and phytoplankton were the dominant (>75%) carbon sources in low-gradient rivers and reservoirs, respectively. Moreover, the longitudinal pattern of resource utilization by consumer FFGs in the Liuxi River was almost the same as that reported in the Waikato River of New Zealand (Pingram et al., 2012b), where the morphological features of hydrophytes and feeding habitats of fish/invertebrates in the two regions were similar. By comparing the habitat F I G U R E 4 Nonmetric multidimensional scaling (NMDS) based on the relative contributions (%) of five resources to the diet assimilation of sitespecific consumer assemblages composed of fish and invertebrate functional feeding groups. AC 3 , aquatic C 3 plants; TC 3 , terrestrial C 3 plants; TC 4 , terrestrial C 4 plants; SPOM, suspended particulate organic matter F I G U R E 5 Relative contributions (%) of five basal resources to the diet assimilation of invertebrate and fish assemblages. S1-S8, sampling sites 1-8; SPOM, suspended particulate organic matter were the determinants influencing the region-specific carbon flows ( Figure 6, also see Wang et al., 2018b;2020c). Given that aquatic vegetation serves as a key substrate for biofilms as well as a habitat and refugium for aquatic organisms, riverine invertebrates prefer to dwell and hide in submerged hydrophytes, and thus, they often scrape the periphyton attached to the blade surface (Wu, 2016). This may also explain why periphyton dominates the consumers' diet assimilation in the upper river. The availability and accessibility to such essential habitats, which support the high productivity of benthic periphyton and biofilms, might be critical to carbon flow transmission along periphyton-based food chains.
From headwaters to the lower reaches, a downstream decrease in periphyton and aquatic C 3 plants with an increase in SPOM and terrestrial C 3 plants, which indicated the longitudinal replacement of autochthonous by allochthonous organic matter along the river, was the basic spatial pattern for the diet assimilation of widespread invertebrate and fish FFGs (Figures 4 and 5).
In contrast to the features of the upstream riffle habitats, the in- The numeric values are the standardized path coefficients and variance estimates. Riffle habitat, % of habitat area; SF, seasonal floodplain area (m 2 ); Rev, revetment (% of riverbank length); cobble substrate, % of substrate area; urban land, % of land use area. DO, dissolved oxygen (mg/L); flow velocity (m/s); TP, total phosphorus (mg/L); COD Mn , chemical oxygen demand determined using the permanganate method (mg/L); NH 3 -N, ammonia nitrogen (mg/L). Crab (g/m 2 ); insect scraper (g/m 2 ); bivalve (g/m 2 ); Ins F, insectivorous fish (g/m 2 ); Det F, detritivorous fish (g/m 2 ). AC 3 and SPOM to I, relative contribution (%) of aquatic C 3 plants and suspended particulate organic matter to the diet assimilation of invertebrate assemblages. AC 3 , TC 3 , and Per to F, relative contribution (%) of aquatic C 3 plants, terrestrial C 3 plants, and periphyton to the diet assimilation of fish assemblages particularly with respect to vegetation degradation (Figure 7b and hydrophyte productivity, such as riffles, cobble/pebble substrates, and floodplains (see Figure 6).

| A suitable concept for explaining the pattern of resource utilization
Carbon flows in large river food webs are dependent on variations in habitat conditions, both temporally in terms of flow variability and spatially in relation to habitat characteristics. Based on recent work using isotope analysis to quantify carbon flows in large river food webs (Pingram et al., 2012a), it appears that autochthonous carbon of benthic algae and, to a lesser extent, hydrophytes and phytoplankton provide the dominant carbon sources fuelling the riverine food webs. This pattern was also observed within the present study for the diet assimilation of invertebrate assemblages at sites 1-3 and of fish assemblages at sites 1-5 ( Figure 5)

| CON CLUS IONS
This study assessed the spatial pattern of resource utilization by aquatic consumers along a subtropical river. Autochthonous resources, including cobble-attached periphyton and aquatic C 3 plants Our results demonstrated that aspects of the riverine productivity model are most suitable for revealing the longitudinal pattern of resource utilization along this large, disturbed river system, which is more similar to other examples in the neotropics and tropics than to those in temperate areas. The findings in the Liuxi River have important implications for revealing energy supply and transmission processes in river food webs, which may have pronounced effects on the functioning of these aquatic ecosystems.

ACK N OWLED G M ENTS
This research was financially supported by the National Natural Science Foundation of China (51909107)

CO N FLI C T O F I NTE R E S T
The authors have no conflicts of interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
All data supporting this study are provided as supporting information accompanying this manuscript.