River Styles and stream power analysis reveal the diversity of fluvial morphology in a Philippine tropical catchment

Characterisation of hydromorphological attributes is crucial for effective river management. Such information is often overlooked in tropical regions such as the Philippines where river management strategies mainly focus on issues around water quality and quantity. We address this knowledge gap using the River Styles Framework as a template to identify the diversity of river morphodynamics. We identify eight distinct River Styles (river types) in the Bislak catchment (586 km2) in the Philippines, showing considerable geomorphic diversity within a relatively small catchment area. Three River Styles in a Confined valley setting occupy 57% of the catchment area, another three in a partly confined valley setting occupy 37%, and two in the remaining 6% are found in a laterally unconfined valley setting. Five characteristic downstream patterns of River Styles were identified across the catchment. We observe that variation in channel slope for a given catchment area (i.e., total stream power) is insufficient to differentiate between river types. Hence, topographic analyses should be complemented with broader framed, catchment-specific approaches to river characterisation. The outputs and understandings from the geomorphic analysis of rivers undertaken in this study can support river management applications by explicitly incorporating understandings of river diversity and dynamics. This has the potential to reshape how river management is undertaken, to shift from reactive, engineering-based approaches that dominate in the Philippines, to more sustainable, ecosystem-based approaches to management.


Introduction
Recognition of morphological diversity and understanding of river processes (dynamics) are essential for effective river management (Brierley and Fryirs 2005;Gurnell et al. 2015;Rinaldi et al. 2016;Hohensinner et al. 2018). Differences between rivers are often contextualised at the reach scale (10 -1 -10 1 km; Belletti et al. 2017) where geomorphic structure and function are approximately uniform, determined by a set of boundary conditions within which the river operates (Frissell et al. 1986;Brierley and Fryirs 2005;Wyrick and Pasternack 2014).
Imposed boundary conditions (i.e., geological and terrain controls) dictate the valley setting and topography of the landscape within which rivers adjust, while flux boundary conditions such as the interaction of water discharge and sediment transport induce reach-scale variability in morphodynamics (Brierley and Fryirs 2005). Alterations to flux boundary conditions, including human-induced pressures and disturbances to water and sediment flows, may irreversibly impact fluvial systems (Rhoads 2020). Integrating principles from hydrology, geomorphology, and ecology strengthen the potential of development of sustainable river management programs (Brierley 2008;Brierley et al. 2019). This requires detailed information of river morphological diversity across multiple spatiotemporal scales (Gurnell et al. 2015).
A range of classification schemes has been developed to assess river morphological diversity. Generally, these schemes seek to categorise reaches by grouping similar process and form characteristics. Classification schemes vary in their approach, the environment for which they were developed, and the spatiotemporal scales over which they are applied (Kondolf et al. 2003;Buffington and Montgomery 2013). However, Kasprak et al. (2016) demonstrate that the underlying principles and premises of such geomorphological analyses are inherently consistent, with different approaches generating similar outputs (i.e., maps of river type). Classification based on hydrology and river geomorphology (hydromorphology) often provides a first step in the analysis of river systems (Fuller et al. 2013), and is a fundamental starting-point when integrating interdisciplinary components of analysis (e.g., Sear et al. 1995;Gilvear 1999;Kondolf et al. 2003;Downs and Gregory 2004;Brierley and Fryirs 2005;Brierley 2008;Meitzen et al. 2013;Tadaki et al. 2014;Rinaldi et al. 2016;Dallaire et al. 2019).
Spatially hierarchical classification frameworks support river management strategies, including those to maintain ecosystem functions (Dollar et al. 2007;Beechie et al. 2010), to mitigate the effects of flood hazards (Rinaldi et al. 2013, and to restore degraded rivers (Beechie et al. 2010;Hawley 2018). These frameworks use a nested spatial hierarchy to organise and structure complex river systems, wherein large-scale features (i.e., a river) are subdivided into sequentially smaller features (e.g., segment to river reach to geomorphic unit/s to microhabitat subsystems; Frissell et al. 1986). The River Styles Framework developed by Brierley and Fryirs (2005) incorporates spatially-hierarchical geomorphic analyses within a catchment-based approach to river management. Knowledge of changes in the catchment (historical and contemporary, from upstream to downstream) is essential to contextualise river adjustment (e.g., lateral, vertical and wholesale) (Brierley and Fryirs 2009), especially as disturbances caused by natural phenomenon and anthropogenic activities may occur anytime and at any position in a catchment (Gurnell et al. 2015). The River Styles Framework provides a set of consistent and generic procedures and guidelines for river assessment that can be applied at local to regional scales (e.g., Kuo and Brierley 2013;Marcal et al. 2017;Khan et al. 2021;Nardini et al. 2020). The framework has four stages: Stage One involves identifying and characterising River Styles; Stage Two uses geomorphic principles to assess evolution, and river condition; Stage Three assesses river recovery potential; Stage Four sets target conditions and priorities to realise effective river management (Brierley and Fryirs 2005).
While the River Styles Framework has been used to support river management in a variety of locations, there has been limited application in tropical regions (Fryirs and Brierley 2021). Tropical rivers exhibit a variety of morphological forms (Latrubesse et al. 2005) and significant geomorphic diversity (Sinha and Latrubesse 2020). Taking examples from across the Philippines, remote sensing analyses have indicated the dynamic planform morphologies of selected rivers. Reach-scale aerial imagery analysis in the Bucao and Santo Tomas Rivers (Zambales) showed spatially variable changes in river width between 2004 and 2013; a legacy effect of sediment supply from the Mt Pinatubo volcanic eruption in 1991 (Paz-Alberto et al. 2016). A national-scale assessment of river migration in the vicinity of 74 large river bridges using Landsat satellite imagery showed that some river reaches have remained approximately stable over the 30-year monitoring period whilst others have considerably adjusted through migration, contraction and expansion (Boothroyd et al. 2021b). River migration rates > 300 m per decade have been quantified along a ~ 85 km segment of the Cagayan River; these migration rates are typically higher than those observed in temperate rivers (Dingle et al. 2019) and can present a challenge for river management.
Rivers in the Philippines play a central role in the community (Fig. 1a), providing economic benefits such as agriculture on floodplains (Fig. 1b), fishing pens in channels (Fig. 1c, and gravel mining on riverbeds and bars (Fig. 1d). While some anthropogenic activities are temporary and have a limited impact on river form and function, the construction of large anthropogenic structures including dams, bridges and embankments ( Fig. 1e-g) imposes a particular width and/or alignment on the channel. These structures restrict the capacity for adjustment, effectively 'fighting' against the prevailing river behaviour (Brierley and Fryirs 2009).
Despite the importance of river morphological characteristics and their associated processes, this information is not routinely considered in the classification of water bodies in the Philippines. Instead, the classification is primarily based on water quality standards set by the Department of Environment and Natural Resources (DENR Administrative Order 2016-08) (DENR-EMB 2016). Furthermore, river basin management plans exist for the major catchments only (drainage area > 3000 km 2 ). Although some data on physiography, climate, geology, and land use exist, analyses on fluvial geomorphological detail (i.e., stream network characteristics, geomorphic units, bed material, sediment, and flow regime) are often lacking. As the result of the gaps in hydromorphological understanding in the catchment management plans, interventions tend to be limited to local and reactive approach such as misplaced dikes and river training measures that are expensive to build and maintain. Considering the increasing magnitude and variability of river flows and growing pressures on water supply from climate change and floodplain land use for agricultural and urban development (Tolentino et al. 2016;Eccles et al. 2019), such management responses could pose serious problems for the future.
Hence, geomorphologically informed approaches to river management in the Philippines are critical.
In this study, we use the River Styles Framework to appraise topographic controls on river morphological diversity and its subsequent implications on river management. First, we identify distinct River Styles through Fig. 1 a The Bislak River, like many rivers in the country, serves multiple purposes in Philippine communities, including: b agriculture on floodplains; c aquaculture/fish pens in channels; and, d the provision of construction materials from gravel mining on riverbeds and bars. Semi-permanent/ permanent structures such as bridges, embankments, dikes, and gabions (e-g) are commonly constructed along Philippine rivers. Examples shown are from the Bislak Catchment; locations are indicated by the corresponding labels (a-g) topographic analysis of a recent digital elevation model (DEM), interpretation of satellite and aerial imagery, and ground-truthing in the Bislak Catchment, Philippines. Then using stream power as a measure of the capacity of rivers to erode and transport sediment, we evaluate whether a simple stream power model can be used to predict patterns in River Styles. Based on the patterns of morphological river diversity in the Bislak Catchment, we recommend a range of transferrable, geomorphologically informed strategies that can contribute towards improved place-based understanding of river character and behaviour in the Philippines, essential for effective river management.

Bislak Catchment and Bislak River overview
The Bislak Catchment (586 km 2 ) is located in the Ilocos Region, Luzon Island, Philippines (Fig. 2a). Its headwaters are bounded by the Luzon Central Cordillera (LCC) Mountain Range with a maximum elevation of 1857 m (Fig. 2b). The catchment shares similar topographic characteristics with other small to medium-sized Philippine catchments with mountainous headwaters having steep terrain and high relief. The average catchment slope is 21.9° (Fig. 2c) and 10% of the catchment area is distributed above 1000 m. The river network flows west and discharges into the South China Sea (average channel slope = 0.044 m/m).
The Bislak River is a valuable resource for the local community; as a water source for domestic and agricultural purposes, a habitat for freshwater species, an economic resource for aggregate (i.e., sand and gravel) extraction, a food source, and a recreational and cultural amenity through the annual River Ritual. Almost half of the population of the municipalities of Bacarra (15,937) and Vintar (15,753) reside on the Bislak River floodplains (Philippines Statistics Authority 2015), in locations highly susceptible to flooding (Paringit and Pascua 2017). In attempts to mitigate the flood and erosion hazards, hard engineering structures have been installed along the riverbanks (e.g., Fig. 1f ), and the riverbed is frequently dredged to reduce the effects of aggradation and to increase channel conveyance. Similar hard engineering structures have been installed in other catchments where communities live adjacent to rivers (e.g., Talisay River; Paz-Alberto et al. 2017).

Climate
There are four climate types in the Philippines based on the distribution of rainfall. The Bislak Catchment has a Type I climate classified as having two distinct seasons, dry from November to April and wet from May to October (Coronas 1920). Annual rainfall totals are high (

Geology
Geology in the Bislak Catchment is heterogenous, with thick sequences of sedimentary rock units and local exposures of intrusive and volcanic igneous rocks (Fig. 2d). The headwaters of the catchment are underlain by conglomerates, breccias, sandstones, and intrusive rocks of the Bangui (Late Eocene-Late Oligocene) and Bojeador (Early Miocene) Formations. Less resistant sedimentary rocks belonging to the Pasuquin Limestone (Late Miocene) and Laoag (Late Early Pliocene-Pleistocene) Formation are exposed in the Rugged Hills (middle reaches), whilst recent alluvial deposits are generally found on the Lowland Plains. Northeast-trending lineaments are apparent along ridges in the headwaters and represent segments of the West Ilocos Fault System (WIFS) (i.e., Vigan-Aggao Fault) ( Fig. 2d; Philippine Institute of Volcanology and Seismology 2015). The Vigan-Aggao Fault is an active oblique sinistral-strike-slip fault which has contributed to landscape deformation in northern Luzon (Rimando and Rimando 2020). The heterogenous geology and presence of active tectonic structures are important controls on the morphostructural characteristics of the Bislak Catchment.

Landscape units and land cover
We identified and mapped three morphostructural regions (landscape units) in the Bislak Catchment: Steep Uplands, Rugged Hills and Lowland Plains based on physiographic setting, landscape position, and relief (Fig. 2e). The landscape units represent portions of the catchment with similar morphostructural characteristics. Steep Uplands which comprise headwater areas account for 59% of the catchment area (mean elevation 705 m). Rugged Hills, 36% of the catchment (mean elevation 170 m), are characterised by partly confined valleys ranging from 25 to 1300 m in width. Only 5% of the catchment area is classified as Lowland Plains (mean elevation 20 m), where the river occupies an alluvial plain with continuous floodplain. These physical attributes impose boundary conditions that control the River Styles and their position in the Bislak Catchment.
Wooded grasslands and grasslands occupy approximately 30% (171.95 km 2 ) and 26% (153.09 km 2 ) of the total land area (Fig. 2f ), respectively. Grasslands are dominant from the Steep Upland to Rugged Hills, while wooded grasslands only become dominant at the base of the Rugged Hills. Open forests occupy 23% of the catchment area (137.74 km 2 ) and tend to be positioned at the eastern edge of the catchment in steep headwaters. Shrublands occupy approximately 10% (59.81 km 2 ) of the catchment area and tend to be positioned on valley bottoms. About 7% (43.67 km 2 ) of the catchment area is occupied by agricultural activities (including annual croplands) positioned along gentler hillslopes and relatively flatter floodplain areas adjacent to the inland water. Barren areas cover 2% (12.38 km 2 ) of the catchment area. Only a small fraction (0.5%) of the catchment is classified as built-up, with scattered urban settlements throughout the rugged hill landscape and municipality centers at the Lowland Plains.

Topographic analysis
Topographic analysis used a nationwide digital elevation model (DEM) acquired in 2013 and generated through airborne IfSAR technology, with 5 m spatial resolution (1 m root-mean-square error vertical accuracy; Grafil and Castro 2014). TopoToolbox V2 was used to extract the stream network and calculate catchment areas using standard flow-routing algorithms (Schwanghart and Scherler 2014). An upstream area threshold value (1 km 2 ) delineated the transition from debris flow-dominated channels to alluvial channels (Montgomery and Foufoula-Georgiou 1993); only alluvial channels (> 1 km 2 ) were considered for analysis. Non-parametric quantile regression was applied to remove data artefacts and errors in longitudinal profiles (i.e., hydrological correction). Quantile carving along the central tendency of the stream network elevation data (τ = 0.5) ensured downstream decreasing elevations (Schwanghart and Scherler 2017). Because channel slope was variable over short distances, slope values were averaged over 0.2 km segment lengths (n = 1129). The mean elevation [m], mean slope [m/m], and median catchment area [km 2 ] were extracted for each segment length and exported to a stream network shapefile. The position of the stream network was validated against recent Google Earth imagery. Thematic maps for land cover, geology and faults provided by the National Mapping and Resource Information Authority (NAMRIA), Mines and Geosciences Bureau (MGB), and Philippine Institute of Volcanology and Seismology (PHIVOLCS) were used in the analysis of the catchment and regional setting.

River Style identification
The procedural tree that underpins Stage One of the River Styles Framework was used to identify different river types across the Bislak Catchment (Brierley and Fryirs 2005;Fryirs and Brierley 2018). The naming convention followed Fryirs and Brierley (2018). Segments of the river with similar confinement (position of the channel on the valley bottom), channel planform (continuity, number of channels, sinuosity), geomorphic unit assemblage, and bed material texture were classified into discrete reaches termed River Styles (Brierley and Fryirs 2005;Khan and Fryirs 2020). Each River Style represents a spatial unit with distinctive hydromorphological attributes. Interpretation of Google Earth imagery, viewed at the reach (~ 1:75,000 to 1:10,000) and unit scales (~ 1:5000 to 1:500) assisted River Style identification. Geomorphic attributes were ground-truthed and verified during field visits to accessible sites. The spatial distribution of River Styles was mapped and the longitudinal profiles for the major tributaries interpreted alongside boundaries from the catchment and regional setting analyses (e.g., Fig. 2bf; including landscape units, valley settings, geology, contemporary process zones, and sediment transport regimes). The approach enabled the identification of controls on the position and type of river across the catchment (Rhoads 2020). Given the multitude of criteria being considered in the River Styles Framework, the classification ultimately depends on the expert judgement and decision of the researchers (Tadaki et al. 2014).

Stream power
Stream power is widely used as an indicator of the capacity of rivers to erode and transport sediment (Jain et al. 2006;Bizzi and Lerner 2015). Total stream power is calculated as: where ). Q is a discharge that controls channel form, for which the bankfull discharge is often used (Petit et al. 2005). In landscape-scale studies, the relationship between catchment area, A, and Q often uses a catchment area-discharge relationship derived for the region. Using data from 14 gauges in NW Luzon (Irrigation Division: Philippines 1924) collected in the early twentieth Century, the 2-year discharge Q 2 was estimated using the R package fasstr (https:// github. com/ bcgov/ fasstr), and related to catchment area as Q 2 = 0.44A 1.03 (n = 14; R 2 = 0.93; se = 0.25) [Eq. 2]. The exponent is close to 1 (95% confidence interval 0.85-1.20), justifying substitution of A for Q 2 in the stream power equation: For constant total stream power, catchment area and slope are thus inversely related, A α S −1 . Thus, Eq. 2 can be solved for 0.2 km river segment lengths using the values for A and S derived from the topographic analysis ("Topographic analysis").

River diversity
Eight River Styles were identified along 246 km of mapped stream length (Table 1 and Fig. 3). River Styles in a Confined valley setting, where > 85% of either channel margin abuts the valley bottom margin, account for 57.1% of the stream length. The confined reaches are located in upstream tributaries of the Steep Upland landscape unit (Fig. 2d), have a single-threaded, low sinuosity channel planform and bedrock or coarse bed material texture (boulder to cobble sized) with sculpted highenergy erosional geomorphic units. In confined reaches, there is limited capacity for the river to adjust laterally or vertically, except locally where isolated pockets of floodplain are present. River Styles in a partly confined valley setting, where 10-85% of either channel margin abuts the valley bottom margin, account for 36.5% of the stream length. These partly confined reaches are mainly located on moderate slopes of tributaries in the Rugged (2) � = 0.44γ AS. Hills landscape unit, and have one to two channels and a mixed bed material texture (boulder to sand-sized).
Mid-channel and bank-attached depositional units occur more frequently in these reaches than in confined reaches upstream. The capacity for adjustment is moderate, particularly where discontinuous floodplain pockets are present, meaning that floodplain stripping and reworking can occur (enabling both lateral and vertical adjustment; Nanson 1986). River Styles in a laterally unconfined valley setting account for only 6.3% of the stream length. These occur within the Lowland Plains landscape unit (Gob et al. 2016). The channel planform is predominantly multi-threaded with a mixed bed material texture (gravel to clay-sized). Geomorphic units in these reaches are products of short and long-term sediment accumulation such as compound islands and bars. Here, the capacity for lateral adjustment is greatest for these rivers, where the floodplain is continuous on both banks of the channel (Fig. 3). Several anthropogenic modifications (Fig. 1) designed to restrict lateral channel adjustment are present in the laterally unconfined river types.

Downstream patterns in River Styles
The distribution of downstream patterns (Fig. 4a) and processes along the river are presented by plotting the stream power against longitudinal profile-catchment area and annotating boundary conditions along the Bislak River (Fig. 5b). Note that stream power is spatially variable (Fig. 5b); calculated at 0.2 km intervals. The sources of all tributaries in the Bislak Catchment are in the Steep Upland landscape unit, and these tributaries flow through confined valleys (Fig. 4a and c) and wider partly confined valleys before joining the trunk stream (Bislak River). Five downstream patterns or sequences of River Styles occur in the catchment (Fig. 4b)

Differentiating river diversity with stream power
There is a general inverse relationship between slope and catchment area (Fig. 5a). Different River Styles occupy different regions of the plot, with confined River Styles occurring on steeper slopes than partly confined River Styles. However, at any given catchment area, more than one River Style is encountered in both  confined and partly confined valley settings occurring at all catchment areas between 5 and 120 km 2 (Fig. 6a). This overlap between River Styles becomes more apparent when all data are plotted (Fig. 6b). The regression lines in Fig. 6b have gradients significantly less than the contours of constant stream power lines, and the gradient of all the data regression is − 0.70 [± 0.03]. The significance of stream power is further emphasised when considering the slope-area data for individual river patterns ( Fig. 6c; see Fig. 4 for pattern definitions). Overall, the slope-area relationships reflect the general shape of the topography (Fig. 6c) which also determines stream power. However, stream power is not the sole control over river morphology (Fig. 6a, b) suggesting that valley topography, stream network structure and hence lithology have significant influences on morphological diversity. Note that the gradient of the constant stream power lines will be steeper if river morphology is controlled by a less frequent flow than bankfull, as the gradient of relationships between discharge and area reduces as flow frequency decreases (Burgers et al. 2014).

River diversity and downstream patterns
Stage One of the River Styles Framework revealed the diversity of fluvial morphology in the Bislak Catchment. The identified River Styles operate under a combination of imposed boundary conditions (including landscape units, lithology, and valley-setting) which influences process zone distributions, interplaying with flux boundary conditions that affect stream power (as a function of slope and contributing catchment area).
In the Bislak Catchment, confined River Styles are characterised by small drainage areas, and steepest slopes (Fig. 6a). Initiating with a Steep Headwater River Style (Fig. 5), a change in lithology from metavolcanics to sedimentary bedrock results in transition to a Gorge River Style. With increasing catchment area and the channel slope remaining high (in the range 0.008-0.24 m/m; Table 1), the stream power gradually increases to a peak of 8340 W/m through the confined valley setting. Even though the stream power remains high (maximum 4305 W/m; Table 1), the storage potential is limited in Occasional Floodplain Pockets given the limited accommodation space in the valley bottom (Figs. 3 and 6c). Upon entering the Rugged Hills landscape unit, the valley confinement transitions to partly confined and the channel slope is reduced (in the range 0.002-0.03 m/m; Table 1). With discontinuous floodplain, there is a switch to Bedrock Margin-controlled and Planform-controlled, Wandering, Discontinuous Floodplain River Styles. Stream power falls to a maximum of 2640 W/m and remains approximately consistent (as increases in discharge from tributary inputs are offset by downstream reductions in slope). At the downstream end of the Lowland Plain landscape unit, the stream power is reduced to < 2270 W/m because of the decrease in channel slope (in the range < 0.001-0.002 m/m; Table 1). Such catchment-based fluvial geomorphic understanding is essential for explaining local hydromorphological patterns and processes.
The relationship between slope and catchment area shows the potential to determine the downstream pattern of River Styles in the Bislak Catchment (Fig. 6). Stream power analysis indicates the general topographic controls on the distribution of River Styles within the catchment. However, the overlap (variation in slope) between River Styles for a given catchment area (e.g., between confined and partly confined River Styles, Fig. 6) suggests that stream power alone is insufficient to differentiate and predict river character. To fully appraise the morphological diversity, a hierarchical approach such as the River Styles Framework can complement topographic analyses to understand river character and behaviour. In the Bislak Catchment, patterns of River Styles reflect geologic controls upon valley setting alongside determinants of flow energy expressed as stream power, including analysis of discharge sediment supply, and vegetation controls (Buffington 2012).
Diversity in River Styles and their downstream patterns have been observed in other tropical river catchments. In the similarly sized (616 km 2 ) but relatively steeper Liwu Basin in the eastern Central Range of Taiwan, confined River Styles dominate the channel network (82% of the stream length; Kuo and Brierley 2013). The downstream patterns of River Styles in the Liwu Basin are similar to those identified in the Bislak Catchment, with the transition from partly confined to laterally unconfined valley setting located towards the lower reaches of the catchment. The downstream patterns of River Styles along the Liwu River have been used as a guide to sediment storage; more than 95% of sediment storage occurred within the laterally unconfined valley setting (Kuo and Brierley 2013). Although comparable measurements of sediment storage have not been made for the Bislak Catchment, laterally unconfined valley settings represent a similarly small proportion of the total stream length between the catchments (6.3% in the Bislak, 13% in the Liwu). Consequently, we might expect sediment storage to be limited to a disproportionately small portion of the Bislak Catchment. Fluctuating sediment supply from monsoon and typhoon related landslides, earthquakes, and volcanic eruptions are reported to be key geomorphic drivers in other Philippine catchments (Gran et al. 2011;Catane et al. 2012;Dingle et al. 2019).
In the relatively larger (1800 km 2 ) Macaé Catchment in northern Rio de Janeiro State, Brazil, diversity in River Styles was most pronounced in the Mountain Escarpment landscape unit where there was variability in slope, valley confinement and valley sinuosity (Marcal et al. 2017). For the Bislak Catchment, six out of the eight identified River Styles were located in the Steep Upland and Rugged Hills landscape units, where the topography and geology were most variable (e.g., Table 1 and Fig. 2). Parallels between the lower reaches of the Macaé and Bislak Rivers can also be drawn. Lowland reaches of the Macaé River were channelized from the 1940s to 1980s for flood protection and the development of agricultural land, resulting in enhanced longitudinal connectivity (flow and sediment flux) and reduced lateral connectivity (Marcal et al. 2017). Due to these structural interventions, sections of the lower Macaé River have been irreversibly changed from a laterally unconfined, meandering sand bed River Style to a laterally unconfined, low sinuosity channelized River Style (Brierley et al. 2019). Anthropogenic structures (e.g., concrete dikes; Fig. 1d, f ) along the laterally unconfined reaches of the Bislak River could modify patterns of longitudinal and lateral connectivity.

Translation of River Styles principles to guide river management
The River Styles Framework, contextualised using topographic and stream power analysis, showed the diversity of River Styles in the Bislak Catchment and the differences in their character and behaviour. Bespoke management strategies will be required to work with the local river character, behaviour, and capacity for adjustment. It is contended that the principles introduced in this paper are required to support sustainable river management across the Philippines. Table 2 outlines several key principles of the River Styles Framework (column 1) that are translated to a simplified understanding (column 2) and can support sustainable river management practices (column 3). Indeed, we argue that effective river management cannot be reliably conducted independent from geomorphic insights into catchment-scale patterns of river character and behaviour.
The Bislak River serves a range of functions for the community (Fig. 1b-c). Neglecting the principles of fluvial geomorphology can lead to reactive management, exemplified by expensive and unsustainable attempts to control the river in addressing hydrometeorological risks (Fig. 1g). In the Bislak Catchment, river management has currently involved construction of multiple structures including flood defences, embankments, gabion walls, and flow deflectors ( Fig. 1e-g). The structures are frequently damaged and require regular repair following high-magnitude flows. The effects of the structures on the geomorphology of the river are often left unaccounted for but may negatively affect flow patterns and sediment flux (Brierley et al. 2011). Disregarding the principles of fluvial geomorphology often leads to mismanagement of a river, posing threats to the ecosystem the river is supporting (Brierley 2008).
Translating the principles of Stage One of the River Styles Framework to the management of the Bislak Catchment (and Philippine rivers more widely) revealed the diversity of River Styles in one small catchment. Although physical processes and principles which control the river form are universal, the context of the river, particularly catchment properties and anthropogenic history mean that the outcomes of these principles will not always be the same and recognising and understanding such differences are critical for effective management (Brierley et al. 2011). The new insights provided by the River Styles Framework (as we demonstrate in this study) can be used to develop pro-active and effective strategies for managing river systems that works with the river. These principles provide an entry point for hydromorphological information to be incorporated into river basin management and catchment action plans. Subsequent stages of the River Styles Framework, focusing on river evolution and the assessment of geomorphic condition as a basis to predict future trajectories and recovery potential Brierley and Fryirs 2016). Such insights can be used to inform policies for various river projects including but not limited to sustainable river gravel quarrying plans, upstream intervention for flood mitigation, and river use management. Lastly, since the River Styles Framework is an involved process and relies on interpretation (i.e., field-assessment, interpreting imagery, etc.), assessment of river character and behaviour over large areas will continue to remain a challenge. That said, the predictions based on topographic and stream power analyses that we demonstrate here could potentially contribute to further studies that automate and upscale these analyses for larger scale roll-out (e.g., Fryirs et al. 2019;Boothroyd et al. 2021a;Khan et al. 2021).

Conclusions
This study utilises the River Styles Framework to show that sustainable river management should account for geomorphic insight into both the character and behaviour of rivers. In the Bislak catchment (586 km 2 ), we identified eight distinct River Styles with Confined (57%) and Partly confined (37%) River Styles dominating the stream network by length. Five downstream patterns influenced by imposed and flux boundary conditions exist in the catchment. These findings indicate that a representative Philippine catchment such as Bislak (ranked 55 in the Philippines by catchment area) exhibits different River Styles and downstream patterns requiring different management strategies. Furthermore, variation in channel slope for a given catchment area (i.e., total stream power) is shown to be insufficient to differentiate the type of river in a given reach, hence the need to understand the differences between channel reaches beyond such metrics. The guiding principles of the River Styles Framework of reading the landscape, working with nature, and knowing your catchment by respecting the diversity and the different patterns of rivers at the catchment scale (Brierley and Fryirs 2005) are starting points for developing catchment-scale visions and moving towards geomorphologically informed, sustainable river management.