Sediment supplies are a fundamental component of alluvial river systems, but the importance of sustained supplies of externally derived sediments for the evolution of meandering planforms remains unclear. Here we demonstrate the importance of sediment supply in enhancing the growth of point bars that influence the rate of sinuosity increase through flow deflections in meander bends. We use an archive of Landsat images of 16 meandering reaches from across the Amazon Basin to show that rivers transporting larger sediment loads increase their sinuosity more rapidly than those carrying smaller loads. Sediment-rich rivers are dominated by downstream-rotating meanders that increase their sinuosity more rapidly than both extensional and upstream-rotating meanders. Downstream-rotating meanders appear to establish larger point bars that expand throughout the meander, in contrast to extensional meanders, which have smaller bars, and upstream rotating meanders, which are characterized by deposition over the bar head. These observations demonstrate that the size and position of point bars within meander bends influences flow routing and thus controls the dominant direction of meander growth. Rivers with low sediment supplies build smaller point bars, which reduces their capacity to increase meander curvature and the resulting sinuosity.

Meandering rivers are characterized by their sinuous planforms and tendency to migrate. A range of studies have highlighted sedimentological (e.g., Schumm and Khan, 1972), biotic (e.g., Tal and Paola, 2007), valley-slope (e.g., Ouchi, 1985), and planform-curvature controls (e.g., Hickin and Nanson, 1975) on the evolution of the meandering planform, yet there remains no quantitative statement on how channel sinuosity should adjust in response to meander migration. Hooke (1984) defined patterns by which the meandering planform can adjust, but the multitude of patterns suggested that the planform response to migration is controlled by a number of variables that might prevent generalization. Schwendel et al. (2015) highlighted the sensitivity of river migration to the heterogeneity of floodplain deposits along the Beni River (Bolivia), documenting that fine-grained sediment matured in distal parts of the floodplain could affect the migration rates and forms of meanders in reaches with generally high channel mobility. Still, a formal evaluation of the nature of the relationship between meander migration and channel sinuosity has been prevented by the absence of detailed observations across a range of meandering river settings.

Meander migration is principally accomplished through curvature-driven flow divergence initiated at the bend entrance. The transverse water-surface slope facilitates sediment excavation from the bank toe, while obliquely oriented near-bed currents transport sediment inwards to the point bar (Dietrich and Smith, 1984). Point bar growth increases the length and curvature of bends, increasing the cross-stream centrifugal force (Hickin and Nanson, 1975) and enhancing cross-stream velocity by topographic steering (Dietrich and Smith, 1983; Legleiter et al., 2011). Constantine et al. (2014) observed that sediment fluxes correlated positively with rates of riverbank erosion and oxbow lake production, implying that the construction of point bars is an influence on meander migration and sinuosity adjustment. Empirical observations have documented the importance of sediment supplies in planform adjustment through bedform development (Braudrick et al., 2009; Bufe et al., 2016). In the River Ystwyth (Wales, UK), sediment supplies led to the accumulation of point-bar complexes, which intensified the development of a sinuous planform despite repeated attempts to artificially straighten the channel (Lewin, 1976). Physically based numerical models have produced similar results, with enhanced bar migration and riverbank erosion predicted along reaches that experience increased bar construction (e.g., Dunne et al., 2010; Schuurman et al., 2016). Because point bars exert a strong control on bank erosion—which determines how sinuosity changes—and rely upon sediment supply for their growth, we assess the relationship between meander migration and channel lengthening, and quantify how they vary with sediment supply.

The Amazon Basin provides an ideal opportunity to evaluate the relationship between sediment supply, meander migration, and sinuosity adjustment over a large range of scales in the absence of engineering influences on channel mobility. The physiography of the basin leads to systematically varying sediment supplies throughout the drainage network, at least as measured in total suspended sediment fluxes (Guyot et al., 1989, 1994; Latrubesse et al., 2005; Filizola and Guyot, 2009). We selected 16 meandering reaches for study, sampling across the physiographic provinces of the Amazon (Fig. 1A; Table DR1 in the GSA Data Repository1). We ensured that the reaches had wide valley margins and were free from visible evidence of confining elements (e.g., bedrock exposures) as discerned from satellite imagery, although we have no way of assessing subtler erodibility differences in alluvium of the kind referred to by Schwendel et al. (2015), which are incorporated into our average effects. Total suspended sediment (TSS) flux data compiled from various sources by Constantine et al. (2014) were normalized by bankfull channel width, because wider channels can convey larger sediment loads, and used as a proxy for the bed material supply to each reach, which was necessary in the absence of direct measurements of bed material flux (Table DR2).

Multispectral Landsat images from 1984 to 2014 were used to delineate bankfull channel margins. Reach imagery for three rivers (Mamoré, Beni, and Madre de Díos) was obtained at a near-annual resolution (Fig. 1A), while the remaining reaches had a temporal spacing of ∼10 yr. The average bankfull channel width of each reach was measured along straight sections of river (n > 20 per reach). The bankfull channel boundary was delineated by selecting a representative pixel threshold at the channel-bank interface using the normalized difference vegetation index from composite images. The boundary was vectorized and used to generate centerlines for each year on record. After Micheli et al. (2004), overlapping centerlines between sequential time steps were then used to calculate reach-averaged channel migration rates (MR), reported in units of channel widths per year.

The near-annual resolution of Landsat imagery available for the Mamoré, Beni, and Madre de Díos rivers allowed us to directly measure the fractional change in sinuosity resulting from channel migration (S*) as the following:
where S1 and S2 are channel sinuosity at the first and second time steps, measured as the ratio of reach length to valley length, and t is time. For the decadal data set where cutoffs were unavoidable, we determined S* using a modified relation derived by Constantine and Dunne (2008, their equation 4), in which:
where M* is the fractional change in channel length measured directly from imagery between time steps, f is the fractional change in channel length resulting from meander cutoff, L is the characteristic length of a cutoff event in channel-width units, M1 is the reach length at the first time step in channel-width units, and n is the number of cutoff events between time steps. Equation 2 was solved using information derived from our Landsat imagery (Table DR3) and from Constantine and Dunne (2008).

We developed a meander symmetry index (σ) to ensure a quantitative and repeatable method for describing the direction of planform geometry change resulting from meander migration (Fig. 1B). A symmetry index was calculated for each meander as defined by the ratio of the eroded areas created in the upstream versus downstream portions of the meander (see the methods in the Data Repository). Large σ values are indicative of downstream rotation, while small values describe more extensional or upstream bend deformation. In total, 281 measurements were made of reach-averaged channel migration rate (MR) and S*, while 16 measurements of σ were formed from reach-median values of individual σ measurements. We selected the reach-median metric to capture broad-scale meandering behavior that would otherwise be obscured by inter-meander variability.

Fractional change in sinuosity resulting from channel migration (S*) increases for rivers with greater rates of sediment supply, a pattern that is consistent for both the annual and decadal data sets (Fig. 2A). Rivers with TSS fluxes >0.05 Mt m−1 yr−1 increase their sinuosity significantly more rapidly than those with lower sediment loads. This relationship suggests that rivers characterized by high TSS fluxes undergo more rapid meander migration and sinuosity increase (Fig. 2B).

Reach-median σ indices positively correlate with increasing sediment load (Fig. 3A). Median σ values range from 0.95 to 3.25 (with variability described by the interquartile range [IQR] 0.50–2.46) for the decadal data, and from 1.39 to 3.55 (IQR 0.43–3.81) for the annual data (Table DR4). S* is positively correlated with differences in meander deformation style (Fig. 3B), with rivers characterized by σ >1.50 experiencing, on average, more than four times greater sinuosity change than rivers where σ <1.50. Frequency distributions of σ values for each reach (Fig. DR1 in the Data Repository) illustrate that most meanders (36%–76%) rotated downstream (1.05 < σ < 5.0). The percentage of meanders with symmetry indices >5.0 increased for sediment-rich rivers in the basin, with up to 14% of their meanders having indices between 10.0 and 15.0. Although median values are all >1.0, a fraction of meanders within each reach (14%–47%) exhibit σ values <0.90, reflecting upvalley meander migration. Most upstream-rotating meanders are in reaches with moderate (<0.035 Mt m−1 yr−1) sediment loads (Fig. DR1; Table DR2). Extensional meanders are consistently rare, comprising only 5%–24% of meanders across all reaches. In three sediment-rich reaches (M1, Hu, and Uc; Fig. 1A), 18%–24% of the meanders are extensional, which may be explained by the large number of temporally clustered cutoff events (2–13), which reduces the population of mature asymmetrical meanders on the reach.

Our results suggest a link between sediment supply and meander deformation, driven by sediment sequestration on point bars. Bar growth augments the curvature-driven components of the force balance driving sediment transport and bank erosion (Dietrich and Smith, 1983). Rivers with greater sediment loads can deposit this material as water shoals over the point bar and momentum diminishes in the lee, or as material is transported across the channel by transverse near-bed currents (Dietrich and Smith, 1984; Frothingham and Rhoads, 2003). Experimental evidence suggests that point bars grow preferentially in the downstream direction as sediment is deposited at the tail or on the margins of the bar (Dietrich et al., 1979; Pyrce and Ashmore, 2005; Braudrick et al., 2009). Frequent deposition associated with greater sediment fluxes should result in faster outward and downstream bar growth, enhancing the displacement and retention of high-velocity flows near the outer bank downstream of the meander apex. Downstream-rotational meanders experience more-sustained curvature-driven transfer of momentum to the outer bank and have well developed secondary circulation cells that facilitate the growth of longer, more developed point bars that begin at the upstream inflection (Abad and Garcia, 2009a, 2009b).

Upstream-rotating meanders were observed in all reaches and were especially common in reaches with low sediment loads. Although reported less frequently in the literature, upstream-rotating meanders have been previously documented (e.g., Seminara et al., 2001; Duan and Julien, 2010). Upstream rotation has been attributed to spatial variations in the erodibility of bank materials (Posner and Duan, 2012; Schwendel et al., 2015), to backwater effects caused by meander cutoff (Schwenk and Foufoula-Georgiou, 2017), and to the upstream valley orientation of tightly curved meanders (Abad and Garcia, 2009b). The role of low sediment supply and potentially slow bar growth in the initiation of upstream rotation requires further examination. However, our observations suggest that upstream-rotating meanders are associated with compound bend development or with sediment deposition at the bar head (Fig. DR2). Protruding sediment lobes are then generated that promote localized bank erosion upstream of the meander apex, similar to the observations of Nanson (1980) and Hagstrom et al. (2018).

We assessed the low-flow subaerial extent of point bars for three meanders on the Beni River (Fig. 4) to discern whether the direction of planform geometry change was related to systematic patterns of bar growth. In meanders with high σ values, bars tended to be longer, more continuous, and wider between inflection points. Conversely, extensional meanders had narrower, sometimes discontinuous, bars with smaller surface areas (Table DR5). Our observations are similar to those made in the laboratory experiments of Abad and Garcia (2009b), who observed larger, better-developed bar forms in downstream-rotational meanders. They attributed this development to the rapid cross-stream transition of the high-velocity core in response to the change in local curvature, which facilitated the development of stronger helical circulation cells and quicker establishment of the inner-bank bar.

Channel stability analysis suggests that in low-sinuosity rivers with poorly developed meanders, point bars are suppressed in favor of migrating alternate bars (Tubino and Seminara, 1990). In the Amazon Basin, this has been observed on rivers with low sediment supplies (e.g., Xingu River), while more sinuous channels (e.g., Beni River) are characterized by fixed point bars and high sediment loads (Monegaglia et al., 2017). Rivers enriched with sediment can construct point bars more rapidly, increasing bend length and channel curvature while facilitating channel migration. A positive feedback is initiated wherein greater sediment loads encourage bar growth and channel curvature, ultimately exceeding the curvature threshold required for the persistence of fixed point bars over migrating bars. Rivers in which meanders are unable to grow as quickly fail to attain the necessary curvature required to establish permanent fixed bars, therefore limiting sinuosity growth.

The direction of meander deformation determines how channel sinuosity evolves in response to meander migration. The fractional rate of sinuosity change is observed to vary as a function of sediment supply for a range of rivers across the Amazon Basin. This relationship between sinuosity change and sediment supply is manifested through the rate of alluvial sequestration onto point bars. Alluvium is typically sequestered downstream of the bar apex, facilitating the down-channel growth of the bar. For upstream-rotating meanders, sediment accumulates near the head of the bar, creating lobes that protrude into the channel, locally increasing curvature and inducing bank erosion upstream of the meander apex. Extensional meanders are rare, reflecting an early stage of meander development. The meander symmetry index, which quantifies the geometric change in meander growth, is also correlated with sediment supply and rates of channel growth. Our results suggest that the evolution of meandering rivers is strongly tied to upstream supplies of sediment. These are responsible for the turnover of floodplain sediments, nutrients, and habitats and the creation and preservation of alluvial architecture important for recording the river’s evolutionary history through time.

We would like to thank Bodo Bookhagen for assistance with the development of channel extraction methodologies, Alex Horton for assistance with curvature data generation, Joan Teng for assistance with Landsat imagery collection, and Rolf Aalto for useful discussions relating to discharge variability. We also thank T.C. Hales and the reviewers for constructive comments that helped improve the clarity and quality of the manuscript .

1GSA Data Repository item 2019289, extended methods, Figures DR1–DR3, and Tables DR1–DR6, is available online at http://www.geosociety.org/datarepository/2019/, or on request from [email protected].