Search USGS⌕ Search

SEARCH · Search USGS

Results for “Earth Surfaces Processes and Landforms”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources
USGS2006

The effects of Hurricane Hugo on suspended-sediment loads, Lago Loiza Basin, Puerto Rico

In the two main tributaries that enter Lago Loiza, Rio Grande de Loiza and Rio Gurabo, 99 600 tonnes of suspended sediment was transported by 58.2??106 m3 of runoff in a 48 h period. The storm-average suspended-sediment concentration in the Rio Grande de Loiza for Hurricane Hugo was 2290 mgl-1, the second lowest for the 12 storms that have been monitored at this site. In Rio Gurabo the storm-average suspended-sediment concentration was 1420 mg l -1, the sixth lowest recorded out of 15 monitored storms. In Quebrada Salvatierra, a small tributary to Rio Grande de Loiza, suspended-sediment concentrations were as low as 33 mg l-1 during peak runoff of 20m3s-1. Normally the suspended-sediment concentrations at this discharge are 300 mg l-1. Hurricane force winds seem to be the most important factor contributing to the lower than expected suspended-sediment loads. High winds caused vegetation and debris to be dislodged and displaced. Debris accumulated on hillslopes and in small channels, blocked bridges and formed debris dams. These dams caused local backwater effects that reduced stream velocities and decreased suspended-sediment loads. -from Author

Earth Surface Processes and Landforms↗
USGS2012

Developing a new stream metric for comparing stream function using a bank-floodplain sediment budget: a case study of three Piedmont streams

A bank and floodplain sediment budget was created for three Piedmont streams tributary to the Chesapeake Bay. The watersheds of each stream varied in land use from urban (Difficult Run) to urbanizing (Little Conestoga Creek) to agricultural (Linganore Creek). The purpose of the study was to determine the relation between geomorphic parameters and sediment dynamics and to develop a floodplain trapping metric for comparing streams with variable characteristics. Net site sediment budgets were best explained by gradient at Difficult Run, floodplain width at Little Conestoga Creek, and the relation of channel cross-sectional area to floodplain width at Linganore Creek. A correlation for all streams indicated that net site sediment budget was best explained by relative floodplain width (ratio of channel width to floodplain width). A new geomorphic metric, the floodplain trapping factor, was used to compare sediment budgets between streams with differing suspended sediment yields. Site sediment budgets were normalized by floodplain area and divided by the stream's sediment yield to provide a unitless measure of floodplain sediment trapping. A floodplain trapping factor represents the amount of upland sediment that a particular floodplain site can trap (e.g. a factor of 5 would indicate that a particular floodplain site traps the equivalent of 5 times that area in upland erosional source area). Using this factor we determined that Linganore Creek had the highest gross and net (floodplain deposition minus bank erosion) floodplain trapping factor (107 and 46, respectively) that Difficult Run the lowest gross floodplain trapping factor (29) and Little Conestoga Creek had the lowest net floodplain trapping factor (–14, indicating that study sites were net contributors to the suspended sediment load). The trapping factor is a robust metric for comparing three streams of varied watershed and geomorphic character, it promises to be a useful tool for future stream assessments.

Piedmont;Difficult Run;Little Conestoga Creek;Ling↗
USGS2013

Velocity Mapping Toolbox (VMT): a processing and visualization suite for moving-vessel ADCP measurements

The use of acoustic Doppler current profilers (ADCP) for discharge measurements and three-dimensional flow mapping has increased rapidly in recent years and has been primarily driven by advances in acoustic technology and signal processing. Recent research has developed a variety of methods for processing data obtained from a range of ADCP deployments and this paper builds on this progress by describing new software for processing and visualizing ADCP data collected along transects in rivers or other bodies of water. The new utility, the Velocity Mapping Toolbox (VMT), allows rapid processing (vector rotation, projection, averaging and smoothing), visualization (planform and cross-section vector and contouring), and analysis of a range of ADCP-derived datasets. The paper documents the data processing routines in the toolbox and presents a set of diverse examples that demonstrate its capabilities. The toolbox is applicable to the analysis of ADCP data collected in a wide range of aquatic environments and is made available as open-source code along with this publication.

Earth Surface Processes and Landforms↗
USGS2012

Where fast weathering creates thin regolith and slow weathering creates thick regolith

Weathering disaggregates rock into regolith – the fractured or granular earth material that sustains life on the continental land surface. Here, we investigate what controls the depth of regolith formed on ridges of two rock compositions with similar initial porosities in Virginia (USA). A priori, we predicted that the regolith on diabase would be thicker than on granite because the dominant mineral (feldspar) in the diabase weathers faster than its granitic counterpart. However, weathering advanced 20 deeper into the granite than the diabase. The 20  -thicker regolith is attributed mainly to connected micron-sized pores, microfractures formed around oxidizing biotite at 20 m depth, and the lower iron (Fe) content in the felsic rock. Such porosity allows pervasive advection and deep oxidation in the granite. These observations may explain why regolith worldwide is thicker on felsic compared to mafic rock under similar conditions. To understand regolith formation will require better understanding of such deep oxidation reactions and how they impact fluid flow during weathering.

Virginia↗
USGS2013

Mass-balance modeling of mineral weathering rates and CO2 consumption in the forested, metabasaltic Hauver Branch watershed, Catoctin Mountain, Maryland, USA

Mineral weathering rates and a forest macronutrient uptake stoichiometry were determined for the forested, metabasaltic Hauver Branch watershed in north-central Maryland, USA. Previous studies of Hauver Branch have had an insufficient number of analytes to permit determination of rates of all the minerals involved in chemical weathering, including biomass. More equations in the mass-balance matrix were added using existing mineralogic information. The stoichiometry of a deciduous biomass term was determined using multi-year weekly to biweekly stream-water chemistry for a nearby watershed, which drains relatively unreactive quartzite bedrock. At Hauver Branch, calcite hosts ~38 mol% of the calcium ion (Ca 2+ ) contained in weathering minerals, but its weathering provides ~90% of the stream water Ca 2+ . This occurs in a landscape with a regolith residence time of more than several Ka (kiloannum). Previous studies indicate that such old regolith does not typically contain dissolving calcite that affects stream Ca 2+ /Na + ratios. The relatively high calcite dissolution rate likely reflects dissolution of calcite in fractures of the deep critical zone. Of the carbon dioxide (CO 2 ) consumed by mineral weathering, calcite is responsible for approximately 27%, with the silicate weathering consumption rate far exceeding that of the global average. The chemical weathering of mafic terrains in decaying orogens thus may be capable of influencing global geochemical cycles, and therefore, climate, on geological timescales. Based on carbon-balance calculations, atmospheric-derived sulfuric acid is responsible for approximately 22% of the mineral weathering occurring in the watershed. Our results suggest that rising air temperatures, driven by global warming and resulting in higher precipitation, will cause the rate of chemical weathering in the Hauver Branch watershed to increase until a threshold temperature is reached. Beyond the threshold temperature, increased recharge would produce a shallower groundwater table and reduced chemical weathering rates.

Maryland↗
USGS2013

Probing the deep critical zone beneath the Luquillo Experimental Forest, Puerto Rico

Recent work has suggested that weathering processes occurring in the subsurface produce the majority of silicate weathering products discharged to the world's oceans, thereby exerting a primary control on global temperature via the well-known positive feedback between silicate weathering and CO 2 . In addition, chemical and physical weathering processes deep within the critical zone create aquifers and control groundwater chemistry, watershed geometry and regolith formation rates. Despite this, most weathering studies are restricted to the shallow critical zone (e.g. soils, outcrops). Here we investigate the chemical weathering, fracturing and geomorphology of the deep critical zone in the Bisley watershed in the Luquillo Critical Zone Observatory, Puerto Rico, from two boreholes drilled to 37.2 and 27.0 m depth, from which continuous core samples were taken. Corestones exposed aboveground were also sampled. Weathered rinds developed on exposed corestones and along fracture surfaces on subsurface rocks slough off of exposed corestones once rinds attain a thickness up to ~1 cm, preventing the corestones from rounding due to diffusion limitation. Such corestones at the land surface are assumed to be what remains after exhumation of similar, fractured bedrock pieces that were observed in the drilled cores between thick layers of regolith. Some of these subsurface corestones are massive and others are highly fractured, whereas aboveground corestones are generally massive with little to no apparent fracturing. Subsurface corestones are larger and less fractured in the borehole drilled on a road where it crosses a ridge compared with the borehole drilled where the road crosses the stream channel. Both borehole profiles indicate that the weathering zone extends to well below the stream channel in this upland catchment; hence weathering depth is not controlled by the stream level within the catchment and not all of the water in the watershed is discharged to the stream

Earth Surface Processes and Landforms↗
USGS2013

Large wood budget and transport dynamics on a large river using radio telemetry

Despite the abundance of large wood (LW) river studies there is still a lack of understanding of LW transport dynamics on large low gradient rivers. This study used 290 radio frequency identification tagged (RFID) LW and 54 metal (aluminum) tagged LW, to quantify the percent of in-channel LW that moves per year and what variables play a role in LW transport dynamics. Aluminum tags were installed and monitored on LW in-transit during the rising limb of a flood, the mean distance traveled by those pieces during the week was 13.3 river kilometers (km) with a maximum distance of 72 km. RFID tagged LW moved a mean of 11.9 km/yr with a maximum observed at 101.1 km/yr. Approximately 41% of LW low on the bank moves per year. The high rate of transport and distance traveled is likely due to the lack of interaction between LW floating in the channel and the channel boundaries, caused primarily by the width of the channel relative to length of the LW. Approximately 80% of the RFID tags moved past a fixed reader during the highest 20% of river stage per year. LW transport and logjam dynamics are complicated at high flows as pieces form temporary jams that continually expand and contract. Unlike most other studies, key members that create a logjam were defined more by stability than jam size or channel/hydrologic conditions. Finally, using an existing geomorphic database for the river, and data from this study, we were able to develop a comprehensive LW budget showing that 5% of the in-channel LW population turns over each year (input from mass wasting and fluvial erosion equals burial, decomposition, and export out of system) and another 16% of the population moving within the system.

North Carolina;Virginia↗
USGS2015

Measurements of the initiation of post-wildfire runoff during rainstorms using in situ overland flow detectors

Overland flow detectors (OFDs) were deployed in 2012 on a hillslope burned by the 2010 Fourmile Canyon fire near Boulder, Colorado, USA. These detectors were simple, electrical resistor-type instruments that output a voltage (0–2·5 V) and were designed to measure and record the time of runoff initiation, a signal proportional to water depth, and the runoff hydrograph during natural convective rainstorms. Initiation of runoff was found to be spatially complex and began at different times in different locations on the hillslope. Runoff started first at upstream detectors 56% of the time, at the mid-stream detectors 6%, and at the downstream detectors 38% of the time. Initiation of post-wildfire runoff depended on the time-to-ponding, travel time between points, and the time to fill surface depression storage. These times ranged from 0·5–54, 0·4–1·1, and 0·2–14 minutes, respectively, indicating the importance of the ponding process in controlling the initiation of runoff at this site. Time-to-ponding was modeled as a function of the rainfall acceleration (i.e. the rate of change of rainfall intensity) and either the cumulative rainfall at the start of runoff or the soil–water deficit. Measurements made by the OFDs provided physical insight into the spatial and temporal initiation of post-wildfire runoff during unsteady flow in response to time varying natural rainfall. They also provided data that can be telemetered and used to determine critical input parameters for hydrologic rainfall–runoff models.

Earth Surface Processes and Landforms↗
USGS2015

Export of fine particulate organic carbon from redwood-dominated catchments

Recently, researchers have recognized the significant role of small mountainous river systems in the transport of carbon from terrestrial environments to the ocean, and the scale of such studies have ranged from channel bed units to continents. In temperate zones, these mountain river systems commonly drain catchments that are largely forested. However, the magnitude of carbon export from rivers draining old-growth redwood forests has not been evaluated to date. Old-growth redwood stands support some of the largest quantities of biomass in the world, up to 350 000 Mg of stem biomass km -2 and soil organic carbon can reach 46 800 Mg km -2 . In north coastal California, suspended sediment samples were collected at three gaging stations for two to four years on streams draining old-growth redwood forests. Carbon content, determined through loss-on-ignition tests, was strongly correlated with turbidity, and continuous turbidity records from the gaging stations were used to estimate annual carbon exports of 1 · 6 to 4 · 2 Mg km -2  yr -1 . These values, representing 13 to 33% of the suspended sediment load, are some of the highest percentages reported in the global literature. The fraction of organic carbon as part of the suspended sediment load decreased with discharge, but reached an asymptote of 5 to 10% at flows 10 to 20 times the mean annual flows. Although larger rivers in this region exhibit high sediment yields (up to 3600 Mg km -2  yr -1 ), mainly attributed to high rates of uplift, mass movement, and timber harvest, the small pristine streams in this study have sediment yields of only 8 to 100 Mg km -2  yr -1 . Because the current extent of old-growth redwood stands is less than 5% of its pre-European-settlement distribution, the present organic carbon signature in suspended sediment loads in this region is likely different from that in the early 20th century. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.

California↗
USGS2015

Soil carbon storage following road removal and timber harvesting in redwood forests

Soil carbon storage plays a key role in the global carbon cycle and is important for sustaining forest productivity. Removal of unpaved forest roads has the potential for increasing carbon storage in soils on forested terrain as treated sites revegetate and soil properties improve on the previously compacted road surfaces. We compared soil organic carbon (SOC) content at several depths on treated roads to SOC in adjacent second-growth forests and old-growth redwood forests in California, determined whether SOC in the upper 50 cm of soil varies with the type of road treatment, and assessed the relative importance of site-scale and landscape-scale variables in predicting SOC accumulation in treated road prisms and second-growth redwood forests. Soils were sampled at 5, 20, and 50 cm depths on roads treated by two methods (decommissioning and full recontouring), and in adjacent second-growth and old-growth forests in north coastal California. Road treatments spanned a period of 32 years, and covered a range of geomorphic and vegetative conditions. SOC decreased with depth at all sites. Treated roads on convex sites exhibited higher SOC than on concave sites, and north aspect sites had higher SOC than south aspect sites. SOC at 5, 20, and 50 cm depths did not differ significantly between decommissioned roads (treated 18–32 years previous) and fully recontoured roads (treated 2–12 years previous). Nevertheless, stepwise multiple regression models project higher SOC developing on fully recontoured roads in the next few decades. The best predictors for SOC on treated roads and in second-growth forest incorporated aspect, vegetation type, soil depth, lithology, distance from the ocean, years since road treatment (for the road model) and years since harvest (for the forest model). The road model explained 48% of the variation in SOC in the upper 50 cm of mineral soils and the forest model, 54%

California↗
USGS2015

Controls on valley spacing in landscapes subject to rapid base-level fall

What controls the architecture of drainage networks is a fundamental question in geomorphology. Recent work has elucidated the mechanisms of drainage network development in steadily uplifting landscapes, but the controls on drainage-network morphology in transient landscapes are relatively unknown. In this paper we exploit natural experiments in drainage network development in incised Plio-Quaternary alluvial fan surfaces in order to understand and quantify drainage network development in highly transient landscapes, i.e. initially unincised low-relief surfaces that experience a pulse of rapid base-level drop followed by relative base-level stasis. Parallel drainage networks formed on incised alluvial-fan surfaces tend to have a drainage spacing that is approximately proportional to the magnitude of the base-level drop. Numerical experiments suggest that this observed relationship between the magnitude of base-level drop and mean drainage spacing is the result of feedbacks among the depth of valley incision, mass wasting and nonlinear increases in the rate of colluvial sediment transport with slope gradient on steep valley side slopes that lead to increasingly wide valleys in cases of larger base-level drop. We identify a threshold magnitude of base-level drop above which side slopes lengthen sufficiently to promote increases in contributing area and fluvial incision rates that lead to branching and encourage drainage networks to transition from systems of first-order valleys to systems of higher-order, branching valleys. The headward growth of these branching tributaries prevents the development of adjacent, ephemeral drainages and promotes a higher mean valley spacing relative to cases in which tributaries do not form. Model results offer additional insights into the response of initially unincised landscapes to rapid base-level drop and provide a preliminary basis for understanding how varying amounts of base-level change influence valley network morphology.

Earth Surface Processes and Landforms↗
USGS2016

Relations between rainfall–runoff-induced erosion and aeolian deposition at archaeological sites in a semi-arid dam-controlled river corridor

Process dynamics in fluvial-based dryland environments are highly complex with fluvial, aeolian, and alluvial processes all contributing to landscape change. When anthropogenic activities such as dam-building affect fluvial processes, the complexity in local response can be further increased by flood- and sediment-limiting flows. Understanding these complexities is key to predicting landscape behavior in drylands and has important scientific and management implications, including for studies related to paleoclimatology, landscape ecology evolution, and archaeological site context and preservation. Here we use multi-temporal LiDAR surveys, local weather data, and geomorphological observations to identify trends in site change throughout the 446-km-long semi-arid Colorado River corridor in Grand Canyon, Arizona, USA, where archaeological site degradation related to the effects of upstream dam operation is a concern. Using several site case studies, we show the range of landscape responses that might be expected from concomitant occurrence of dam-controlled fluvial sand bar deposition, aeolian sand transport, and rainfall-induced erosion. Empirical rainfall-erosion threshold analyses coupled with a numerical rainfall–runoff–soil erosion model indicate that infiltration-excess overland flow and gullying govern large-scale (centimeter- to decimeter-scale) landscape changes, but that aeolian deposition can in some cases mitigate gully erosion. Whereas threshold analyses identify the normalized rainfall intensity (defined as the ratio of rainfall intensity to hydraulic conductivity) as the primary factor governing hydrologic-driven erosion, assessment of false positives and false negatives in the dataset highlight topographic slope as the next most important parameter governing site response. Analysis of 4+ years of high resolution (four-minute) weather data and 75+ years of low resolution (daily) climate records indicates that dryland erosion is dependent on short-term, storm-driven rainfall intensity rather than cumulative rainfall, and that erosion can occur outside of wet seasons and even wet years. These results can apply to other similar semi-arid landscapes where process complexity may not be fully understood.

Arizona↗
USGS2016

Three-dimensional flow structure and patterns of bed shear stress in an evolving compound meander bend

Compound meander bends with multiple lobes of maximum curvature are common in actively evolving lowland rivers. Interaction among spatial patterns of mean flow, turbulence, bed morphology, bank failures and channel migration in compound bends is poorly understood. In this paper, acoustic Doppler current profiler (ADCP) measurements of the three-dimensional (3D) flow velocities in a compound bend are examined to evaluate the influence of channel curvature and hydrologic variability on the structure of flow within the bend. Flow structure at various flow stages is related to changes in bed morphology over the study timeframe. Increases in local curvature within the upstream lobe of the bend reduce outer bank velocities at morphologically significant flows, creating a region that protects the bank from high momentum flow and high bed shear stresses. The dimensionless radius of curvature in the upstream lobe is one-third less than that of the downstream lobe, with average bank erosion rates less than half of the erosion rates for the downstream lobe. Higher bank erosion rates within the downstream lobe correspond to the shift in a core of high velocity and bed shear stresses toward the outer bank as flow moves through the two lobes. These erosion patterns provide a mechanism for continued migration of the downstream lobe in the near future. Bed material size distributions within the bend correspond to spatial patterns of bed shear stress magnitudes, indicating that bed material sorting within the bend is governed by bed shear stress. Results suggest that patterns of flow, sediment entrainment, and planform evolution in compound meander bends are more complex than in simple meander bends. Moreover, interactions among local influences on the flow, such as woody debris, local topographic steering, and locally high curvature, tend to cause compound bends to evolve toward increasing planform complexity over time rather than stable configurations.

Earth Surface Processes and Landforms↗
USGS2016

Episodic bedrock erosion by gully-head migration, Colorado High Plains, USA

This study explores the frequency of bedrock exposure in a soil-mantled low-relief (i.e. non-mountainous) landscape. In the High Plains of eastern Colorado, gully headcuts are among the few erosional features that will incise through the soil mantle to expose bedrock. We measured the last time of bedrock exposure using optically stimulated luminescence dating of alluvial sediment overlying bedrock in gully headcuts. Our dating suggests that headcuts in adjacent gullies expose bedrock asynchronously, and therefore, the headcuts are unlikely to have been triggered by a base-level drop in the trunk stream. This finding supports the hypothesis that headcuts can develop locally in gullies as a result of focused scour in locations where hydraulic stress during a flash flood is sufficiently high, and/or ground cover is sufficiently weak, to generate a scour hole that undermines vegetation. Alluvium dating also reveals that gullies have been a persistent part of this landscape since the early Holocene.

Colorado↗
USGS2016

Characterizing supraglacial meltwater channel hydraulics on the Greenland Ice Sheet from in situ observations

Supraglacial rivers on the Greenland ice sheet (GrIS) transport large volumes of surface meltwater toward the ocean, yet have received relatively little direct research. This study presents field observations of channel width, depth, velocity, and water surface slope for nine supraglacial channels on the southwestern GrIS collected between 23 July and 20 August, 2012. Field sites are located up to 74 km inland and span 494-1485 m elevation, and contain measured discharges larger than any previous in situ study: from 0.006 to 23.12 m 3 /s in channels 0.20 to 20.62 m wide. All channels were deeply incised with near vertical banks, and hydraulic geometry results indicate that supraglacial channels primarily accommodate greater discharges by increasing velocity. Smaller streams had steeper water surface slopes (0.74-8.83%) than typical in terrestrial settings, yielding correspondingly high velocities (0.40-2.60 m/s) and Froude numbers (0.45-3.11) with supercritical flow observed in 54% of measurements. Derived Manning's n values were larger and more variable than anticipated from channels of uniform substrate, ranging from 0.009 to 0.154 with a mean value of 0.035 +/- 0.027 despite the absence of sediment, debris, or other roughness elements. Ubiquitous micro-depressions in shallow sections of the channel bed may explain some of these roughness values. However, we find that other, unobserved sources of flow resistance likely contributed to these elevated n values: future work should explicitly consider additional sources of flow resistance beyond bed roughness in supraglacial channels. We conclude that hydraulic modelling for these channels must allow for both sub- and supercritical flow, and most importantly must refrain from assuming that all ice-substrate channels exhibit similar hydraulic behavior, especially for Froude numbers and Manning's n. Finally, this study highlights that further theoretical and empirical work on supraglacial channel hydraulics is necessary before broad scale understanding of ice sheet hydrology can be achieved. This article is protected by copyright. All rights reserved.

Earth Surface Processes and Landforms↗
USGS2016

Introduction to special issue on carbon and landscape dynamics

In October, 2013, at the Geological Society of America annual meeting, a theme session focused on carbon and landscape dynamics. That event led to interest in producing a special issue in ESPL compiling papers on this subject. The 13 papers collected for this special issue reflect the diversity of recent geomorphic research, across a range of climatic and geomorphic settings, addressing some aspect of carbon dynamics.

Earth Surface Processes and Landforms↗
USGS2016

Parsing anomalous versus normal diffusive behavior of bedload sediment particles

Bedload sediment transport is the basic physical ingredient of river evolution. Formulae exist for estimating transport rates, but the diffusive contribution to the sediment flux, and the associated spreading rate of tracer particles, are not clearly understood. The start-and-stop motions of sediment particles transported as bedload on a streambed mimic aspects of the Einstein–Smoluchowski description of the random-walk motions of Brownian particles. Using this touchstone description, recent work suggests the presence of anomalous diffusion, where the particle spreading rate differs from the linear dependence with time of Brownian behavior. We demonstrate that conventional measures of particle spreading reveal different attributes of bedload particle behavior depending on details of the calculation. When we view particle motions over start-and-stop timescales obtained from high-speed (250 Hz) imaging of coarse-sand particles, high-resolution measurements reveal ballistic-like behavior at the shortest (10 −2 s) timescale, followed by apparent anomalous behavior due to correlated random walks in transition to normal diffusion (>10 −1 s) – similar to Brownian particle behavior but involving distinctly different physics. However, when treated as a ‘virtual plume’ over this timescale range, particles exhibit inhomogeneous diffusive behavior because both the mean and the variance of particle travel distances increase nonlinearly with increasing travel times, a behavior that is unrelated to anomalous diffusion or to Brownian-like behavior. Our results indicate that care is needed in suggesting anomalous behavior when appealing to conventional measures of diffusion formulated for ideal particle systems.

Earth Surface Processes and Landforms↗
USGS2016

Channel-planform evolution in four rivers of Olympic National Park, Washington, U.S.A.: The roles of physical drivers and trophic cascades

Identifying the relative contributions of physical and ecological processes to channel evolution remains a substantial challenge in fluvial geomorphology. We use a 74-year aerial photographic record of the Hoh, Queets, Quinault, and Elwha Rivers, Olympic National Park, Washington, U.S.A., to investigate whether physical or trophic-cascade-driven ecological factors—excessive elk impacts after wolves were extirpated a century ago—are the dominant controls on channel planform of these gravel-bed rivers. We find that channel width and braiding show strong relationships with recent flood history. All four rivers have widened significantly in recent decades, consistent with increased flood activity since the 1970s. Channel planform also reflects sediment-supply changes, evident from landslide response on the Elwha River. We surmise that the Hoh River, which shows a multi-decadal trend toward greater braiding, is adjusting to increased sediment supply associated with rapid glacial retreat. In this sediment-routing system with high connectivity, such climate-driven signals appear to propagate downstream without being buffered substantially by sediment storage. Legacy effects of anthropogenic modification likely also affect the Quinault River planform. We infer no correspondence between channel geomorphic evolution and elk abundance, suggesting that trophic-cascade effects in this setting are subsidiary to physical controls on channel morphology. Our findings differ from previous interpretations of Olympic National Park fluvial dynamics and contrast with the classic example of Yellowstone National Park, where legacy effects of elk overuse are apparent in channel morphology; we attribute these differences to hydrologic regime and large-wood availability.

Washington↗