Abstract
Reconstruction of active channel geometry from fluvial strata is critical to constrain the water and sediment fluxes in ancient terrestrial landscapes. Robust methods—grounded in extensive field observations, numerical simulations, and physical experiments—exist for estimating the bankfull flow depth and channel-bed slope from preserved deposits; however, we lack similar tools to quantify bankfull channel widths. We combined high-resolution lidar data from 134 meander bends across 11 rivers that span over two orders of magnitude in size to develop a robust, empirical relation between the bankfull channel width and channel-bar clinoform width (relict stratigraphic surfaces of bank-attached channel bars). We parameterized the bar cross-sectional shape using a two-parameter sigmoid, defining bar width as the cross-stream distance between 95% of the asymptotes of the fit sigmoid. We combined this objective definition of the bar width with Bayesian linear regression analysis to show that the measured bankfull flow width is 2.34 ± 0.13 times the channel-bar width. We validated our model using field measurements of channel-bar and bankfull flow widths of meandering rivers that span all climate zones (R2 = 0.79) and concurrent measurements of channel-bar clinoform width and mud-plug width in fluvial strata (R2 = 0.80). We also show that the transverse bed slopes of bars are inversely correlated with bend curvature, consistent with theory. Results provide a simple, usable metric to derive paleochannel width from preserved bar clinoforms.
INTRODUCTION
Reconstruction of formative channel geometry from fluvial strata is critical to constrain the ancient hydrology and terrestrial mass fluxes on Earth and other planets (e.g., Bhattacharya et al., 2016), unravel fluvial responses to past climate change (Foreman et al., 2012), and aid hydrocarbon exploration (e.g., Miall and Tyler, 1991). Alluvial river geometry is characterized by streamwise channel-bed slope (S) and bankfull flow depth (hbf) and width (Bbf). Robust methods—tested with extensive modern observations, numerical simulations, and physical experiments—exist for estimating hbf and S from fluvial strata. For example, hbf has been estimated from the geometry of preserved river dune deposits (Paola and Borgman, 1991; Leclair and Bridge, 2001; Jerolmack and Mohrig, 2005) or the preserved bar-clinoform surfaces (relict bank-attached and free channel-bar surfaces; Mohrig et al., 2000; Hajek and Heller, 2012; Alexander et al., 2020). Channel-bed slope has been estimated from the elevations of correlative downstream channel architectures (Bhattacharya et al., 2016) or using empirical bankfull Shields stress criteria observed in modern rivers (Paola and Mohrig, 1996; Trampush et al., 2014). However, we lack similar tools to estimate Bbf from fluvial strata.
Current methods to reconstruct Bbf from fluvial strata leverage fully preserved channel architectures (Ielpi and Ghinassi, 2014), channel-belt dimensions (Gibling, 2006; Ielpi et al., 2017), and the scaling of Bbf and hbf in extant rivers (Leeder, 1973; Hayden et al., 2019). Fully preserved abandoned channel fills offer a direct estimate of Bbf, but they are rarely preserved due to rechannelization or limited outcrop exposure (Bridge, 2003; Bhattacharya et al., 2016; Toonen et al., 2012). Abandoned channel dimensions can also be measured in seismic data sets (Bhattacharya et al., 2016) and on inverted topography where planform channel architecture is preserved at the terrestrial surface (Ielpi and Ghinassi, 2014). In the absence of these observations, Bbf is estimated from channel-body widths, or in conjunction with estimated hbf based on the width-to-depth scaling observed in modern alluvial rivers (Bridge, 2003; Hayden et al., 2019). These estimates are hampered by substantial uncertainty because channel-body dimensions can differ significantly from formative channel dimensions (Hayden et al., 2019), and the scaling between Bbf and hbf depends on channel sinuosity (Bridge, 2003), varying by an order of magnitude in extant rivers (Trampush et al., 2014). While a more mechanistic method to estimate Bbf of single-threaded rivers exists (Lapôtre et al., 2019), it requires detailed measurements of bed and bank-material grain size and is yet to be expanded to include the cohesive effects of grain-size mixtures in the bank material and floodplain vegetation.
Following previous work (Allen, 1965; Ethridge and Schumm, 1978; Bridge, 2003; Bhattacharya et al., 2016), we propose that the geometry of meandering river channels is encoded in the size of bank-attached forced bars (point bars) and can be reconstructed from their deposits. Point bars are macroforms for which size scales with channel dimensions (Allen, 1965; Mohrig et al., 2000; Hajek and Heller, 2012). Point bar sediments can be readily identified in fluvial strata at outcrop scales (Fig. 1) and in seismic data (Jackson, 1976; Durkin et al., 2017), and their internal accretionary surfaces record lateral channel migration (Allen, 1965). Point bar deposits are also underrecognized in sandy deposits, suggesting a greater abundance of bar clinoform surfaces than is often considered (Hartley et al., 2015; Chamberlin and Hajek, 2019).
The ability to use preserved bar clinoforms as proxies for Bbf would unlock the potential for detailed morphodynamic reconstructions of ancient landscapes. We combined high-resolution topographic data from meandering rivers with Bayesian linear regression analysis to develop a robust relation between Bbf and the point bar surface width, Wbar, which facilitates Bbf reconstruction from fluvial strata.
ANALYSIS OF MODERN RIVERS
We used lidar elevation data from 11 meandering rivers across the United States, available through OpenTopography (https://opentopography.org) (sampled at 1 m) and the U.S. Geological Survey (USGS) 3D Elevation Program (https://www.usgs.gov/core-science-systems/ngp/3dep) (sampled at 1 m and 1/3 arc-seconds; Fig. 2A). These rivers exhibit active channel migration and natural levee development and are not incised into their floodplains. Our data set spans bankfull depths of 2 m to 26 m (field estimates were available for nine reaches; the empirical relation of Ielpi and Lapôtre [2020] was used for hbf estimation for two reaches) and reach lengths of 19 km to 183 km, and it covers six meander bends per study reach on average. The data set covers the temperate, arid, and cold Köppen climate classifications (Beck et al., 2018), bed slopes of 1 × 10−4 to 1 × 10−3, and median discharges of 15 to 18,000 m3/s.
For each study reach, we mapped channel centerlines from Landsat 8 satellite imagery using RivaMap (Isikdogan et al., 2017). This method produced discontinuous centerlines when applied to rivers with small channel widths (<60 m) relative to the image resolution (30 m). We manually corrected all programmatically generated centerlines to fill gaps and remove stems in the centerline from adjoining tributaries. We georeferenced the channel centerlines to the elevation products and systemically sampled each bend by generating at least three thalweg-perpendicular cross sections per bar at even spacing across the bend inlet, apex, and outlet (Fig. 2B). We computed Bbf as the cross-stream distance between manually picked channel levees (Fig. 2C).

RESULTS
Our model (Equation 3) predicted Bbf for the compiled natural and numerical rivers (R2 = 0.79; Fig. 3D) and performed significantly better than the existing power-law (R2 = 0.55) and linear (R2 = 0.23) width-to-depth scaling relations. Furthermore, the order-of-magnitude uncertainty range in the width-to-depth scaling relations is significantly greater than both the high probability density interval and the posterior predictive interval (Equations 3 and 4) used in our method.
APPLICATION TO FLUVIAL STRATA
To test the model's applicability to ancient strata, we compiled paired measurements of fully preserved bar clinoform widths and mud-plug widths from four published outcrop panels and two published seismic cross sections (details are provided in the Supplemental Material1). From these six images, we independently estimated Bbf from the lateral extent of fully preserved abandoned channel fills. We measured Wbar from fully preserved clinoforms that were correlative to the measured channel fills such that we sampled the same channel body (e.g., within the sixth-order bounding surfaces of Miall [1988]). We digitized the clinoform surfaces, fit Equation 1, and followed the same methodology applied to the modern bar surfaces to estimate Wbar.
Sample size, cross-section position relative to the bar apex, and the cross-section obliquity relative to the paleocurrent direction are factors that likely influence the application of our model to fluvial strata. To test model sensitivity to these effects, we assumed that the bar clino-form preservation in fluvial strata is randomly distributed with respect to cross-section obliquity and bend position, and we generated channel cross sections at varying angles, θ, to the centerline direction across Trinity River (Texas, USA) bars, yielding a data set of channel cross sections at θ ∈ [0°, 90°]. From this data set, we systematically sampled 1–15 cross sections at random angles and bend positions to assess sensitivity of α to sample size. We also sampled three cross sections at each θ to assess the sensitivity of α to θ.
Model predictions of Bbf using Equation 3 were consistent with the correlative measurements of abandoned channel widths (R2 = 0.80; Fig. 3D), validating our model application to ancient strata. Our sensitivity tests revealed that a robust representation of Wbar and inferred Bbf from fluvial strata can be achieved by sampling ≥3 bar clinoforms from the same channel body (e.g., sixth-order bounding surfaces of Miall, 1988; Fig. 4B). For cross sections with θ ≲ 70°, α remained consistent (Fig. 4A) but never fully converged to Equation 3 (Fig. 4B). Multiple measurements from the same lateral accretion set (e.g., third-order bounding surfaces of Miall [1988], or multiple clinoform surfaces like those pictured in Fig. 1A) could be problematic because this violates the assumption of random bend positions. In the absence of multiple samples of preserved bar clinoforms across a channel complex, Equation 4 can be used to estimate Bbf from a single Wbar measurement that has been geometrically corrected to paleoflow perpendicular direction. While our model is directly applicable to fully preserved clinoforms, the symmetry of Equation 1 may enable model application by fitting the sigmoid to a partially preserved bar clinoform, so long as at least half of the original bar surface is preserved.
DISCUSSION AND CONCLUSIONS
We present a robust method to estimate Bbf of meandering rivers from preserved point bar surfaces in fluvial strata. Using high-resolution lidar data and field observations of Wbar and Bbf from 19 rivers, and concurrent field measurements of preserved bar clinoform width and mud-plug widths from fluvial strata, we showed that the formative channel width is 2.34 ± 0.13× the measured channel-bar width. Our model is minimally sensitive to sample size and the position and angle of the preserved bar clinoforms with respect to the paleocurrent direction (Figs. 4A and 4B). These results enable the robust inversion of Bbf from fluvial strata and confirm previously proposed heuristic relations between Wbar and Bbf (Ethridge and Schumm, 1978; Bhattacharya et al., 2016). The consistency between the predicted channel width and measured mud-plug widths (Fig. 3D) indicates that Equation 3 can substantially reduce uncertainty in ancient Bbf estimates compared to existing methods. While Equation 4 provides a lower-precision estimate of Bbf when limited samples are available from fluvial strata, the credible interval is of higher precision than the indirect methods based on hbf. Further, our model can also be used to estimate Bbf from bar clinoforms documented in Martian fluvial strata (Goudge et al., 2018).
We also explored the controls on the variability of bar shapes between rivers. In Equation 1, L is set by hbf, and we found that k is related to the transverse bed slope of bars, ∂z/∂n (Fig. 4C). Theoretical, experimental, and field studies have explored the controls on ∂z/∂n (Struiksma et al., 1985; Talmon et al., 1995), and our results revealed that ∂z/∂n averaged across the bars of the individual study reaches is inversely proportional to their radius of curvature relative to channel depth (Fig. 4D), consistent with theory. While cross-section position relative to the channel apex will directly influence ∂z/∂n (Kleinhans et al., 2012), this observation suggests that quantitative characterization of bar clinoform shape (Eq. 1) may enable the comparison of average bend curvatures across paleochannels.
Finally, quantification of preserved bar deposit and paleochannel widths can enable reconstructions of paleoriver mobility. Robust Bbf estimates can inform empirical relations between channel width and the lateral migration rates of meandering rivers (Ielpi and Lapôtre, 2020). Equation 1 can also be used to quantify the extent of vertical preservation of formative bar topography, which encodes the relative time scales of channel migration and avulsion (Chamberlin and Hajek, 2019; Ganti et al., 2020). Thus, detailed measurements of the size and shape of preserved bars in fluvial strata can constrain the geometry and rates of ancient river evolution, which are central to unraveling fluvial responses to boundary condition perturbations on Earth and river mobility on Mars.
ACKNOWLEDGMENTS
We thank A. Ielpi and two anonymous reviewers for constructive feedback. This work was supported by National Science Foundation grants to Ganti (EAR 1935669) and Hajek (EAR 1935513).