Benthic community structure and ecosystem functions in above- and below-waterfall pools in Borneo

Waterfalls are geomorphic features that often partition streams into discrete zones. Our study examined aquatic communities, litter decomposition and periphyton growth rates for above- and below-waterfall pools in Ulu Temburong National Park, Brunei. We observed higher fish densities in below-waterfall pools (0.24 fish m−2 vs. 0.02 fish m−2 in above-waterfall pools) and higher shrimp abundance in above-waterfall pools (eight shrimp/pool vs. less than one shrimp/pool in below-waterfall pools). However, macroinvertebrate densities (excluding shrimp) were similar among both pool types. Ambient periphyton was higher in below-waterfall pools in 2013 (4.3 vs. 2.8 g m−2 in above-waterfall pools) and 2014 (4.8 vs. 3.4 g m−2 in above-waterfall pools), while periphyton growth rates varied from 0.05 to 0.26 g m−2 days−1 and were significantly higher in below-waterfall pools in 2014. Leaf litter decomposition rates (0.001 to 0.024 days−1) did not differ between pool types, suggesting that neither shrimp nor fish densities had consistent impacts on this ecosystem function. Regardless, this research demonstrates the varied effects of biotic and abiotic factors on community structure and ecosystem function. Our results have highlighted the importance of discontinuities, such as waterfalls, in tropical streams.

While availability of food resources is important for regulating both community structure and ecosystem functions, physical factors are also important (Bond & Downes, 2000). Waterfalls are fast flowing, rocky biotopes with distinct geomorphic structure (e.g. channel slopes and stream channels typically erode down to parent material). These biotopes flow vertically without obstruction and are generally more than 1 m in height (Newson & Newson, 2000). In spite of their potential importance in tropical streams, waterfalls have not received as much attention as other biotopes (Rackemann et al., 2013;Clayton & Pearson, 2016). In high elevation localities, waterfalls become common because the stream gradient steepens and the river channels form discrete, sequential pools and waterfalls (or cascades and riffles). The discontinuities in channel form caused by waterfalls are barriers to the upstream dispersal of fish and thus create discrete communities (Creed, 2006;Covich et al., 2009;Hein & Crowl, 2010;Karssing et al., 2012;El-Sabaawi et al., 2015). In tropical streams, below-waterfall pools support a range of predacious and herbivorous fish (Choy, 1996); in contrast, above-waterfall pools tend to have few fish and often become a refuge for shrimp (Covich et al., 2009;Hein & Crowl, 2010). However, Bass (2007) notes that different aquatic communities in above-and below-waterfall pools may be the result of differences in microhabitats or other environmental conditions.
The different community structures of fish and shrimp above and below waterfalls can influence a range of ecosystem functions including organic matter decomposition and periphyton growth (Greathouse & Pringle, 2006;Moulton et al., 2010;Ho & Dudgeon, 2016). Below waterfalls, high numbers of fish can play an important role in controlling periphyton growth by clearing rock surfaces after frequent tropical flood events (Pringle & Hamazaki, 1998). Above waterfalls, where shrimp can be abundant, organic matter breakdown rates may be increased by consumption and shredding of leaf litter by resident shrimp (Crowl et al., 2001;March et al., 2001;Coughlan et al., 2010). Fish and shrimp can also reduce macroinvertebrate densities, as many are insectivorous (Pringle & Hamazaki, 1998). Macroinvertebrates, a term which refers to invertebrates excluding shrimp in this study, are also important in restricting periphyton growth (Feminella & Hawkins, 1995;Steinman, 1996;Moulton et al., 2015) and breaking down organic matter (Cheshire et al., 2005;Masese et al., 2014). Therefore, fish and shrimp may have an indirect impact on ecosystem functions by reducing macroinvertebrate densities. The separation of communities and its impact on ecosystem structure and function have been researched to some extent, with mixed results, in the Neotropics (Pringle & Hamazaki, 1998;March et al., 2002;Ramirez & Hernandez-Cruz, 2004;Covich et al., 2009). However, there have been few similar studies conducted in tropical Asia despite the similarity of the region's biotic and abiotic conditions to those of the Neotropics (Mantel et al., 2004;Ho & Dudgeon, 2016).
The natural fragmentation of tropical streams by waterfalls provides a ''natural experiment'' for examining these biotic and physical factors. This study investigated four streams in Ulu Temburong National Park, Brunei Darussalam, all with waterfalls that reduce the upstream dispersal of fish, creating a refuge for shrimp in above-waterfall pools. The aim of this study was to investigate the extent to which waterfalls influence fish community structure and how expected declines in fish densities in above-waterfall pools affect the density and biomass of both shrimp and nonshrimp macroinvertebrates. In addition, we tested the hypothesis that ecosystem functions (i.e. periphyton growth and leaf litter breakdown rates) differed between above-and below-waterfall pools. We hypothesized that both shrimp and macroinvertebrates occur in higher densities in the above-waterfall pools lacking fish and that periphyton growth rates are higher in above-waterfalls pools owing to the lower number of herbivorous fish. Leaf litter breakdown rates were also hypothesized to be higher in abovewaterfall pools owing to the higher numbers of shrimp and macroinvertebrates.

Study sites
Small tributaries of Sungai Belalong and Sungai Temburong in the Temburong district, Brunei Darussalam (in northern Borneo), were the focus of this study (Fig. 1). -The four streams (Sungai Lower Apan, Sungai Baki, Sungai Esu and Sungai Apan Threelan; Fig. 1) were in the vicinity of the Kuala Belalong Field Studies Centre (KBFSC). They were chosen because they have perennial flows and catchments covered in primary rainforest (Cranbrook & Edwards, 1994). The climate in Brunei Darussalam is tropical and weakly influenced by the Southeast Asia monsoon (Dykes, 1996). The daily weather pattern in the rainforest of Ulu Temburong is very erratic, with most rain originating from convection cells formed over the ocean. As the cells move inland, they are forced to rise over the mountains of Bukit Belalong and Bukit Pagon, thereby condensing and producing daily heavy rain (Cranbrook & Edwards, 1994). Topography ranges from 30 m asl at the Kuala Belalong Field Study Centre but rises to 913 m asl at Bukit Belalong and to 1,850 m asl at Bukit Pagon (Dykes, 1994). The streams drain V-shaped valleys with steep bank slopes.
Water quality can have a strong impact on benthic communities and ecosystem functions (Giller & Malmqvist, 1998;Everaert et al., 2014). This study was therefore conducted in streams that are not under any anthropogenic impacts in the catchment (Cranbrook & Edwards, 1994). Further, the similar geology and catchment conditions suggest water quality should be relatively uniform among all the study streams. In a pilot study, nitrate and phosphate were found to occur at concentrations less than 0.55 and 0.08 mg l -1 , respectively. Any variation in ecosystem structure and function in the streams of Ulu Temburong National Park are thus assumed to be attributable to biotic interactions among resident taxa (e.g. competition for resources and predation) and to natural abiotic conditions, including physical disturbances cause by variable river discharges.
Physical disturbances are assumed to be similar across all pool biotopes owing to the consistent gradient of the streams. Light levels were also assumed to be similar across our study pools, with gaps in canopies (of approximately 70-90%) created by the waterfalls extending to both the above-and below-waterfall pools. Lower Apan was the longest reach (90 m) of the four reaches, followed by Apan Threelan (75 m), Esu (70 m) and Baki (55 m). The waterfalls at Threelan and Baki were furthest from their respective confluences with Temburong compared to the other study streams (separated by over 300 m). This may affect the numbers of fish present in the tributaries, with fish dispersing into the tributaries as refuges during high discharges. The waterfall at Esu was approximately 150 m from its confluence with Belalong, and Lower Apan was approximately 60 m from its confluence with Temburong (Online Appendix 1).
Background community structure A total of 44 species of fish from 10 families and 30 genera have been identified in the streams and rivers of Ulu Temburong National Park (Choy & Chin, 1994). Species of Cyprinoidea were predominant, with families Balitoridae (river loaches) and Cyprinidae being the most commonly represented (Choy & Chin, 1994). In general, cyprinids have a mixed diet that changes with their life stage and includes insects, algae, diatoms and higher plants (McConnell & Lowe-McConnell, 1987), while the endemic genus Gastromyzon (Balitoridae) comprises herbivorous fish that feed on periphyton exclusively (Tan, 2006).
Three genera of freshwater shrimp live in the streams of Ulu Temburong National Park with the most common being Macrobrachium (Decapoda: Palaemonidae; Wowor, personal communication). Macrobrachium are omnivorous, but predominantly predaceous (Crowl & Covich, 1994). The other shrimp in Ulu Temburong, which were not present in this study, include Caridina (Decapoda: Atyidae), which are found in headwater streams, and Atyopsis (Decapoda: Atyidae), which are present in streams with strong currents (Wowor, personal communication). A previous study of non-shrimp macroinvertebrates in these streams (Baker et al., in press) identified 14 orders, with Diptera being the most abundant in number of individuals (approximately 40%), followed by Coleoptera (approximately 20%), Ephemeroptera (approximately 20%) and Trichoptera (9%).

Physical habitat measurements
Wetted and bankfull widths and channel depths were measured in three pools above the waterfall and three pools below the waterfall, except at Esu and Lower Apan, where there were only two pools below the waterfall. At these locations, the benthic substrate was assessed visually and characterized by percentage cover of gravel, cobble, boulder and bedrock. A formal cross-section of velocity and depth was taken at the widest point of each pool every 0.5 m following common methods described by Gordon et al. (2004). Stream velocity was measured using an electromagnetic flowmeter (Valeport Ò model 801; Valeport, Totnes, UK). A two-way analysis of variance (ANOVA) was used to assess differences in the width, wetted width, depth and velocity of the study pools among streams and between the pools above and below the waterfalls. This analysis was carried out in the statistical computing environment R (R Core Team 2013).

Macroinvertebrates, shrimp and fish sampling
Macroinvertebrates and shrimp were sampled in three pools above each waterfall and in all the pools below the waterfall to the confluence. Surber samples (0.1 m 2 ; 250-lm mesh) were used to collect macroinvertebrates, excluding shrimp. Three Surber samples were randomly selected in each pool where it was shallow enough to sample and composited together to create one sample. A combination of three Surber samples and a single 3-min kick sample with a standard D-frame net (500-lm mesh) was used to sample freshwater shrimp in each pool. Shrimp and macroinvertebrates were considered separately, as shrimp are significantly larger than the other macroinvertebrates and are more vagile. Samples were preserved in 70% ethanol in the field and processed under (910) magnification in the laboratory. Individuals were identified to the lowest practical taxonomic level (most often to genus or morphotype to an assumed similar level) and measured to the nearest 0.5 mm. Macroinvertebrate taxonomic knowledge in northern Borneo is poor; therefore, identifications were conducted using the few keys available (e.g. Dudgeon, 1999;Yule & Yong, 2004;Sangpradub & Boonsoong, 2007), alongside open source identification methods. Specifically, taxa were photographed and uploaded onto the Flickr website (http://flickr.com/photos/tropical-streams/ sets/), where interested experts commented or requested to see the actual specimens. Ash-free dry mass (AFDM) of macroinvertebrates was estimated using length-mass regressions (Benke et al., 1999;Sabo et al., 2002;McNeely et al., 2007). In cases in which no taxon-specific equations were available, estimates were made using equations from taxa with similar body shapes . Where only dry mass (DM) estimates were available, values were converted to AFDM following the calculation described by Waters (1977). All macroinvertebrates were considered to be capable of feeding on periphyton and organic matter except known predators, i.e. Odonata, Hemiptera and Hydracarina (Jinggut, 2015). Unfortunately macroinvertebrate and shrimp specimens from Baki were severely damaged in transit, so the analysis for macroinvertebrates and shrimp includes only those from Lower Apan, Esu and Threelan.
In contrast to the macroinvertebrate sampling procedures, fish were only sampled in pools directly above and below the first waterfall upstream from each confluence. Fish were sampled with a cast net (mesh size, 0.64 cm 2 ) thrown into each of the sample pools three times. Fish were photographed and identified following the methods described by Choy & Chin (1994). Fish were categorized as herbivorous or predacious using appropriate scholarly literature, and percentages of fish belonging to these categories were calculated (McConnell & Lowe-McConnell, 1987;Tan, 2006;Kottelat, 2013). Once processed, fish were immediately returned to streams.
ANOVAs were used to assess differences in biodiversity of macroinvertebrates and shrimp among study sites and between pools above and below the waterfalls, followed by Tukey's tests to identify significant factors. As fish were only sampled in pools directly above and below the waterfall, it was not possible to assess differences among study sites. For the analysis, each pool was considered as a replicate, rather than the samples within each pool. To satisfy the assumptions of normality and homoscedasticity, fish and shrimp data were transformed [using a log 10 (n ? 1) transformation]. Macroinvertebrate assemblage structure in above-and below-waterfall pools were examined via a non-metric multidimensional scaling (NMDS) analysis using abundance data followed by an analysis of similarities (ANOSIM; Clarke, 1993) to test for differences in abundance and composition of macroinvertebrates among the above-and below-waterfall pools. Surber samples were combined to give a single estimate per pool. The statistical computing environment R (R Core Team 2013) and the R packages mass and vegan were used to conduct the ANOVA, NMDS and ANOSIM procedures. The Bray-Curtis dissimilarity index, a dissimilarity index for ecological data, was used for the NMDS procedure (Borcard et al., 2012).

Ecosystem function estimates
Leaf litter decomposition rates and periphyton growth estimates were conducted in pools directly above and below the waterfall at each study stream.

Leaf litter decomposition estimates
Leaves from the genus Campnosperma (Family: Anacardiaceae) were used to estimate decomposition rates in pools above and below waterfalls. Although there are approximately 250 species of trees in Ulu Temburong National Park (Cranbrook & Edwards, 1994), Campnosperma was most abundant along the study streams and therefore chosen for this experiment. Leaves from one large branch that had recently fallen into one of the study streams were used in this experiment. Decomposition rates were estimated using standard litter bag techniques (Irons et al., 1994;Benstead, 1996). Prior to the experiment, leaves were dried at 70°C for 24 h, cooled and then weighed. Three grams of leaves were placed into nylon litter bags (mesh size, approximately 8 mm). Leaf litterbags were placed at Lower Apan, Esu, Apan Threelan and Baki. Litter bags were collected from the study sites after 7, 23, 26 or 71 days in the field and processed immediately in the KBFSC laboratory. The leaves were picked out by hand and the remaining material was processed through a 250-lm sieve, with macroinvertebrates picked out and preserved in 70% ethanol. Macroinvertebrates were then identified and measured using the methods described above. Leaf material remaining was dried at 70°C for 24 h, cooled and weighed. Breakdown rates (k d ) were calculated using the equation: Litter bag experiments were conducted in both 2013 and 2014. In the first year, three replicates per study pool were left in study pools for 7-and 23-day exposures. In the second year, the same methods were repeated except four replicates were used, and leaf litter bags were deployed for longer periods of time, 26 and 71 days, in order to take account of the scouring effects of high discharge.

Periphyton growth
Periphyton growth rates were determined as the standing stocks through time and estimated from rocks located in the same pools used for the leaf decomposition experiment. Random rocks from each stream were selected, placed in boiling water for 10 min and scrubbed with a wire brush to remove all periphyton. Rocks were then placed into the study pools with plastic wire used to secure them to the riverbed. After exposure, periphyton abundance was estimated from rock scrubs. A fixed area (0.002 m 2 )of each rock was sampled with material removed and collected onto GFF filters. Filters were dried at KBFSC and transported back to the UK for AFDM analysis. AFDM was calculated by taking the difference between the dry biofilm weight and ash weight (obtained from heating samples to 500°C for 1 h). AFDM was then used as an estimate of periphyton abundance. Periphyton growth rates were calculated as the difference between final and initial standing stocks (AFDM) divided by the total exposure time (days). It should be noted that our estimates of periphyton cannot distinguish between living and dead algal cells (Wellnitz & Poff, 2006). This method was selected because of the remote location and limited facilities at KBFSC. Periphyton experiments were conducted in both 2013 and 2014. In the first year, four rocks per study pool were left for 23-day exposures. In the second year, the same methods were repeated, but with five replicate rocks, and each was deployed for 71-day exposures. In each year, ambient levels of periphyton were estimated from four random rocks in each study pool in order to ensure that the periphyton levels from the experiment were consistent with the background standing stock.
The statistical computing environment R (R Core Team 2013) was used to conduct a one-way ANOVA to evaluate differences in leaf litter breakdown rates and periphyton growth rates between the above-and below-waterfall pools. Assumptions of normality and homoscedasticity were tested with graphical methods, and in order to satisfy the statistical assumptions, many of the datasets collected were transformed (Thomas et al., 2013). A log 10 (n ? 1) transformation was applied to the 7-day exposure leaf litter bag macroinvertebrate data and the 23-day exposure leaf litter bag macroinvertebrate biomass and abundance data from 2013. A square root transformation was applied to the 2013 and 2014 ambient periphyton rock scrub data as well as the 2014 litter bag data on macroinvertebrates, including the 28-day exposure macroinvertebrate biomass and richness data, the 71-day exposure macroinvertebrate biomass data, the 28-day exposure non-predacious macroinvertebrate abundance data and the 71-day exposure macroinvertebrate biomass data.

Physical habitat
There were no significant statistical differences in pool widths, depths or velocities either among streams or between above-and below-waterfall pools (Online Appendix 1). The substrates of above-waterfall pools were 36% bedrock compared to 7% bedrock in belowwaterfall pools. In contrast, there was a higher substrate percentage of cobbles in below-waterfall pools (31%) compared to above-waterfall pools (17%).
It was not possible to identify all fish to the species level, particularly the smallest individuals. However, six distinct taxa were sampled across the four streams with a total of 75 individuals collected during the survey (Online Appendix 6). The most abundant fish were cyprinids, comprising 97% of observed fish, with the most abundant species being Nematabramis steindachneri Popta, 1905 (32 fish) and Rasbora Macroinvertebrate density was estimated using Surber samples and shrimp data from Surber and kick samples for each pool both above and below the waterfalls. Error bars represent standard errors (all below-waterfall pools, n = 8; all above-waterfall pools, n = 9; Sungai Apan Threelan, n = 6; Sungai Esu, n = 5; Sungai Lower Apan n = 6). No error bars are presented for fish density because all samples were combined in the field. There was a significant difference between above-and belowwaterfall densities of fish (P = 0.04) and shrimp (P B 0.01), but no difference in macroinvertebrate densities (P = 0.89). There was also a significant difference in macroinvertebrate densities among the study streams (P = 0.01) argyrotaenia (Bleeker, 1850) (24 fish). The other 3% of fish were Gastromyzon. Average fish densities were significantly lower in above-waterfall pools (0.02 ± 0.02 fish/m 2 vs. 0.24 ± 0.27 fish/m 2 ; F 1,6 = 6.00, P = 0.04). There were no fish in above-waterfalls pools at Baki and Threelan, but fish were present in above-waterfalls pools at Esu (Gastromyzon) and at Lower Apan (Cyprinidae).
In 2014, there were significantly higher leaf litter breakdown rates in above-waterfall pools over 26 days (0.005 vs. 0.004 days -1 ; F 1,29 = 24.80, P B 0.001), but no difference over 71 days (0.008 vs. 0.005 days -1 ; F 1,28 = 3.40, P = 0.07; Fig. 3). Over the 71-day exposure, there was a higher abundance of macroinvertebrates accumulated in the leaf litter bags of below-waterfall pools (F 1,28 = 12.78, P = 0.002). However, there were no significant correlations found between leaf litter breakdown rate and macroinvertebrate richness (R 2 = 0.17), abundance (R 2 = 0.13) or biomass (R 2 \ 0.01). In both 2013 and 2014, some of the leaf litter bags were scoured out during storm events. Overall in both experiments, leaf litter bag recovery was [83%. Leaf litter breakdown rates significantly differed among streams in 2013 over 23-day exposures (P B 0.01). Higher leaf litter decay rates were also found in 2013 over 23-day exposures below the waterfalls (P B 0.01). There was no significant difference in macroinvertebrate biodiversity or shredder diversity between pools above and below the waterfalls

Discussion
Waterfalls can affect the distribution of aquatic animals by partitioning streams into discrete zones (Covich et al., 2009;Hein & Crowl, 2010). Tropical headwater streams in Ulu Temburong National Park contain a large number of shrimp and fish, with waterfalls acting as semi-permeable filters that create different aquatic communities in reaches above and below waterfalls. We found below-waterfall pools supported a higher density of fish and significantly higher ambient periphyton levels in 2013 and 2014 (Table 4). Higher growth rates of periphyton were also found in below-waterfall pools in 2014. Above-waterfalls pools supported a higher abundance and biomass of shrimp, but there was no difference in macroinvertebrate biodiversity between above-and below-waterfall pools and no consistent statistical difference in leaf litter decomposition rates in 2013 or 2014. These findings illustrate how geomorphology can isolate biotopes, creating discontinuities in streams that affect community structure, specifically with respect to fish and shrimp. However, the corresponding impact on ecosystem functions is less straightforward, suggesting differing degrees of influence of biotic and abiotic factors on ecosystem functions. There was a significant difference in leaf litter breakdown rates between pools above and below waterfalls over 71-day exposures (P B 0.01). Over 71-day exposures, a higher abundance of macroinvertebrates accumulated in below-waterfall pools (P = 0.002). A higher abundance of non-predacious macroinvertebrates also occurred in pools below waterfalls (P = 0.01), whilst non-predacious invertebrate biomass was higher in above-waterfall pools (P = 0.001) Biodiversity in above-and below-waterfall pools As expected, the highest densities of fish were found in below-waterfall pools. We found waterfalls ([6 m in height) acted as a barrier to fish at Esu, Threelan and Baki, likely creating a refuge for shrimp in the pools above waterfalls (Covich et al., 2009;Hein & Crowl, 2010). Choy (1996) studied the distribution ecology of fish in the Ulu Temburong National Park and reported similar findings, with no Cyprinidae taxa recorded above any waterfall greater than 5 m in height.
Gastromyzon fish were an exception; these herbivorous fish have a ventral sucker along their bodies and are therefore able to ascend vertical bedrock (Tan, 2006). However, despite this trait, surprisingly few Gastromyzon were observed in pools above waterfalls. This may be the result of the sampling design, as many Gastromyzon species inhabit faster flowing biotopes (Sheldon, 2011), which were not sampled in our study.
Macrobrachium, which were collected in Lower Apan, Esu and Threelan, are important members of aquatic communities. However, their impact on benthic environments has been mainly studied in Neotropical streams, with only a few studies conducted in Asia (e.g. Mantel et al., 2004). In the Neotropics, Macrobrachium have been found to exhibit aggressive behaviour, which can cause other shrimp to take shelter (Crowl & Covich, 1994;March et al., 2001). This behaviour may explain the absence of other shrimp taxa sampled in our study sites. There were no significant differences in shrimp biomass and density among the streams. Among all pools above the waterfalls, there were significantly higher shrimp densities (Fig. 2). These results concur with similar findings in Puerto Rico, where shrimp occur in high abundances in above-waterfall pools (Pringle et al., 1993;March et al., 2002;Covich et al., 2009).
However, the high abundance of shrimp in abovewaterfall pools did not appear to have a significant impact on macroinvertebrate communities, with the NMDS and ANOSIM procedures finding no difference in community patterns between macroinvertebrates in pools above and below the waterfalls. In the headwater streams of Puerto Rico, Ramirez & Hernandez-Cruz (2004) conducted a shrimp exclusion experiment, which revealed Macrobrachium had mixed impacts on aquatic insects (i.e. the effects of shrimp were negligible in some streams). In another study in Puerto Rico, March et al. (2001) found that the impact of shrimp exclusion on insect biomass varied and was dependent on the genus of shrimp. In the highaltitude pools, exclusion of shrimp, which included Atya spp. and Xiphocaris elongata (Guérin-Méneville, 1855), resulted in significantly lower insect biomass. In contrast, exclusion of Macrobrachium spp. in the low altitude pools promoted higher insect biomass. In our study, the low numbers of Macrobrachium in belowwaterfall pools had no apparent effect on Growth rates were higher below the waterfalls in 2014 (P B 0.01) but not in 2013 (P = 0.3). In 2014, there was also a difference in growth rates among streams (P = 0.02). *Only one experimental rock in the pool above the waterfall on Sungai Threelan due to the rocks being scoured out during flood events macroinvertebrate biomass, while in Costa Rica, fish and shrimp significantly reduced aquatic insects with combined effects greater than the summed effect of both groups individually (Pringle & Hamazaki, 1998). These studies conducted in the Neotropics suggest that freshwater shrimp have an important role in structuring benthic community composition, but their effects vary and may depend on the presence or absence of other biota (March et al., 2002). In Hong Kong, Mantel et al. (2004) found no effect of the removal of Macrobrachium on aquatic insects; however, the authors were unable to determine the reason there was no effect. One explanation is that the macroinvertebrates consistently evaded fish and shrimp predators; this has been shown in studies that artificially doubled the natural densities of fish and failed to observe a decrease in some invertebrate taxa (Gilliam et al., 1989;Dudgeon, 1991). Our study also did not reveal a consistent pattern in macroinvertebrate diversity of above-and below-waterfall pools, suggesting that communities of shrimp and fish do not influence these biotic patterns; however, it is possible that fish and shrimp predators are substitutable and thus functionally redundant (sensu Ho & Dudgeon, 2016). In Northern Australia, Garcia et al. (2015) found that top-down effects on macroinvertebrates were context dependent with respect to factors such as benthic habitat, stream velocity and community structure, each of which influenced invertebrate diversity. Ulu Temburong has frequent flood events; for example, after a storm, it is not uncommon for the level of Belalong to rise by 1 m in a 30-min period (Cranbrook, 1993). These 'flashy flows' scour out habitats, physically removing taxa, and therefore may be more important in structuring macroinvertebrates than the presence of shrimp or fish (Bond & Downes, 2000). The hypothesis that predation and competition are less important in highly disturbed streams and rivers has been suggested for systems in other parts of the world (Bishop, 1973;Peckarsky et al., 1990;Yang & Dudgeon, 2010). Flash floods can create relatively unpredictable systems (Boulton et al., 2008) characterized by a continuum of conditions that may ultimately result in two distinct environments occurring at various times within a given portion of a stream. For example, during low flows, tropical streams are complex systems exhibiting a mix of flow biotopes (i.e. pools, riffles and cascades) and functional habitats (i.e. wood debris, leaf litter, cobbles and gravel; sensu Harper et al., 1995, Harvey et al., 2008; however, during a flood event, these streams become homogeneous as water rises to form a uniform flood biotope.

Ecosystem functions
It was assumed that differences in fish and shrimp densities would affect macroinvertebrate patterns along with leaf litter decomposition rates, ambient periphyton levels and periphyton growth rates. Specifically, past studies in the tropics have suggested that leaf litter decomposition can be increased by higher shrimp densities (March et al., 2001;Wright & Covich, 2005), and herbivorous fish can increase periphyton levels (Pringle & Hamazaki, 1997). We found ambient periphyton standing stocks and associated growth rates to be linked with higher fish densities, while leaf litter decomposition rates were unaffected by varying fish or shrimp densities. The observed leaf litter decomposition rates were in the range of those from other tropical studies reported by Wantzen et al. (2008) with processing rates ranging from 0.001 days -1 (Mathooko, 1998;Rueda-Delgado et al., 2006) to 0.789 days -1 (Mathooko & Kariuki, 2000). However, the lack of a correlation between fish, shrimp and macroinvertebrate occurrence parameters and breakdown rates suggests that leaf litter loss in our study sites was caused by other factors. Microbial processing has been suggested to be more important in the tropics compared to temperate regions (Irons et al., 1994). However, the importance of microbes to the decomposition of leaf litter is poorly understood, especially in the tropics (Boyero et al., 2011a, b), and thus, more research is required to understand the role of bacteria and fungi in breakdown rates. Mechanical abrasion from frequent floods in catchments can affect leaf litter breakdown, as has been demonstrated in other tropical studies (Pearson et al., 1989). Fast flows, which occur frequently in streams within Ulu Temburong, have been reported to cause physical abrasion of leaf litter, breaking debris down into smaller pieces and therefore influencing decomposition rates (Wantzen et al., 2008).
We found no significant differences in leaf litter decay between above-and below-waterfall pools. Many streams that flow through dipterocarp forests, such as our study sites, have very high tree diversity and leaves that are high in lignin and low in protein (Yule et al., 2009); such leaves possess an unappealing combination of traits for shredders. Landeiro et al. (2008) excluded fish and shrimp (mainly Macrobrachium) from leaf litter bags and found greater leaf litter breakdown in the control relative to the exclusion area. This could be linked with higher insect shredder abundances; however, there was no difference in shredder abundance between treatments, and the results were inconclusive. In spite of more recent research suggesting insect shredders do exist in tropical streams (Boyero et al., 2011a, b), our regression analysis showed no correlation between leaf litter breakdown rates with macroinvertebrate diversity. In our study, some leaf litter bags lacked any macroinvertebrates, for example, those in Lower Apan (e.g. the 2013, 23-day exposure bags in the below-waterfall pools; Table 3), suggesting that leaf litter is an unappealing food source for local macroinvertebrates.
Periphyton growth rates appear to be low compared to ambient periphyton levels, which is unsurprising because of the abundant herbivorous grazers and frequent scouring flows. However, some tropical fish have been found to be vitally important in increasing periphyton growth rates (Pringle & Hamazaki, 1997;Moulton et al., 2010), which may explain high levels of ambient periphyton in below-waterfall pools and low periphyton growth rates in the above-waterfall pool of Apan Threelan (0.04 g m -2 days -1 ) where no fish were present. Overall abundance of ambient periphyton was correlated with the occurrence of fish. Pringle & Hamazaki (1997) made similar findings; they observed fish influencing algal community composition, whereas shrimp had no significant effect. Studies have found herbivorous fish to be efficient in cleaning rock surfaces of deposited sediment after storm events, therefore helping increase periphyton levels (Pringle & Hamazaki, 1997). This is important in tropical streams, such as our study streams, where fast flows are frequent and can both cause physical scouring and cover rocks in fine sediment (Wantzen et al., 2008). However, rates of accrual including grazing and scouring losses of periphyton among our study streams are not known.
Although many of the fish sampled in our study streams were not herbivorous, the literature contains studies in which even predacious fish have been shown to create trophic cascades. Moulton et al. (2010) found that predacious fish inhibit shrimp and baetid mayflies from grazing in two Neotropical streams, thus increasing both abundance and quality of periphyton in pools containing fish. However, other experiments have found that top-down effects of fish on algae have less impact relative to bottom-up effects such as increases in nutrient levels (Garcia et al., 2015;Ho & Dudgeon, 2016). Interestingly, Ho & Dudgeon (2016) found no impact of high numbers of fish and shrimp on algal biomass or periphyton accumulation in three Hong Kong streams. This was suggested to be caused by the removal of macroconsumers, which may have reduced nutrient levels and therefore contributed to the apparently weak top-down effects on fish and shrimp (Ho & Dudgeon, 2016). Other factors, such as variation in sunlight across sites, may also affect periphyton levels; this requires further investigation across our study sites. However, as our experiments were conducted in pools directly above and below individual waterfalls, the pairs of sites experienced similar light levels.

Conclusion
As hypothesized, this study clearly illustrates that waterfalls affect habitat patterns of fish and shrimp, with waterfalls (generally [5 m in height) acting as a barrier to the upstream dispersal of fish. This natural habitat fragmentation influenced algal biomass, but not litter breakdown. Our results show that even in a highly disturbed system, biotic interactions are still important. Further studies must be conducted to understand the effects of other factors, such as frequent scouring flows, which can have greater impacts than biotic factors. These findings suggest that streams divided into naturally discrete units by waterfalls have considerable ecological and conservation significance. Recent land-use changes caused by the widespread growth of the palm oil industry in Southeast Asia have increased the urgency of identifying and studying the remaining pristine rivers and streams. Understanding how natural discontinuities, such as waterfalls, can affect habitat patterns and ecosystem functions is vitally important to the successful management and conservation of these systems. helpful field assistant. This work was funded by a Natural Environment Research Council PhD studentship. Special thanks are also extended to Roy Whiles, Paul Clark, David Reid, Stephen Butler, Albert Orr, Manfred Jäch, Rüdiger Wagner, Virginija Podeniene, Marija Ivković, Herbert Zettel, Ignac Sivec, Ben Price, Michel Sartori and Jean Luc Gattolliat for their time in assisting with taxa identification. We thank Prof. Naselli-Flores, Dr. Checo Colón-Gaud and the two anonymous reviewers for their insightful comments and help to improve this text.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.