Alluvial channel response to environmental perturbations: Fill-terrace formation and sediment-signal disruption

. The sensitivity of fluvial systems to tectonic and climatic boundary conditions allows us to use the 10 geomorphic and stratigraphic records as quantitative archives of past climatic and tectonic conditions. Thus, fluvial terraces that form on alluvial fans and floodplains as well as the rate of sediment export to oceanic and continental basins are commonly used to reconstruct paleo-environments. However, we currently lack a systematic and quantitative understanding of the transient evolution of fluvial systems and their associated sediment storage and release in response to changes in base level, water input, and sediment input. Such knowledge is necessary to quantify past environmental change from terrace records or 15 sedimentary deposits, and to disentangle the multiple possible causes for terrace formation and sediment deposition. Here, we use a set of seven physical experiments to explore terrace formation and sediment export from a single, braided channel that is perturbed by changes in upstream water discharge or sediment supply, or through downstream base-level fall. Each perturbation differently affects (1) the geometry of terraces and channels, (2) the timing of terrace cutting, and (3) the transient response of sediment export from the basin. In general, an increase in water discharge leads to near-instantaneous channel 20 incision across the entire fluvial system and consequent local terrace cutting, thus preserving the initial channel slope on terrace surfaces, and it also produces a transient increase in sediment export from the system. In contrast, a decreased upstream sediment-supply


Introduction
Sediment-routing systems are commonly subdivided into three zones: a sediment-production zone, typically a mountainous region; a transfer zone of alluvial and fluvial systems that transport and/or temporarily store sediment; and a sedimentation (or deposition) zone, comprising continental or oceanic basins (Fig. 1; Allen, 2017;Castelltort and Van Den Driessche, 2003).Because climate and tectonics can affect sediment production rates, any changes in those conditions may lead to the formation of fluvial terraces in the transfer zone or changes in sedimentation rates in the deposition zone (Alloway et al., 2007;Bull, 1990;Scherler et al., 2015;Zhang et al., 2001).Many past studies have used such records to reconstruct paleoenvironmental conditions (fluvial terraces: Litty et al., 2016;Poisson and Avouac, 2004;Schaller et al., 2004; sedimentation rates: Hay et al., 1988;Zhang et al., 2001).Quantitative interpretations of either record, however, require a clear understanding of how terraces are formed or how sedimentary signals are altered in the transfer zone (Romans et al., 2016 and references therein).In addition, both records suffer from ambiguity, because variability in different environmental parameters can produce similar sedimentary responses.For example, changes in either sediment or water inputs can create fill terraces (Scherler et al., 2015) and affect sediment deposition rates (e.g., Armitage et al., 2011;Simpson and Castelltort, 2012).
Alluvial rivers adjust their slope and width with respect to the local base level such that in a graded (steady) state, the incoming water discharge (Qw) can transport the incoming sediment supply (Qs,in) downstream (Buffington, 2012;Gilbert, 1877;Lane, 1955;Mackin, 1948).When graded, the slope (S) scales nearly linearly with the ratio of Qs,in and Qw (e.g., Blom et al., 2017;Malatesta and Lamb, 2018;Parker, 1979;Wickert and Schildgen, 2019): Changes in boundary conditions (Qw, Qs,in, and base level) therefore cause alluvial rivers to adjust their geometries through sediment deposition (aggradation) or incision, until a new graded profile is reached.Incision or aggradation result from the dependence of bedload-transport capacity on slope and water discharge (Meyer-Peter and Müller, 1948).For example, if Qw increases while Qs,in is held constant, the transport capacity exceeds Qs,in, which leads to the entrainment of additional sediment from the channel bed, thus incision.As incision proceeds, the channel slope decreases until the transport capacity drops to match Qs,in.Conversely, if Qs,in exceeds the transport capacity of the channel, sediment will be deposited to steepen the channel, thus increasing the transport capacity until it matches Qs,in.These adjustments can be recorded through (1) fillterrace formation in the transfer zone (e.g., Bridgland and Westaway, 2008;Bull, 1990;Merritts et al., 1994) and (2) changes in sediment export to basins (e.g., Allen, 2008;Castelltort and Van Den Driessche, 2003;Romans et al., 2016).
When studying terraces in the field, it can be difficult to distinguish between terraces that mark a sudden switch from aggradation or stable conditions to incision ("fill-top" terraces of Howard, 1959) and those that preserve surfaces that were cut by a river moving laterally during a period of overall incision ("fill-cut" terraces of Bull, 1990 andPazzaglia, 2013).In the latter case, there can be a substantial lag between the onset of the environmental perturbation and the abandonment of the terrace surface (e.g., Steffen et al., 2010Steffen et al., , 2009)).Consequently, from fill terraces alone, both the formation mechanism (change in Qw, Qs,in or base level) and the timing of the perturbation can be ambiguous.
Numerical and experimental work has demonstrated that the geometrical adjustment of alluvial rivers to external perturbations not only creates fluvial terraces, but also affects sediment discharge at the outlet (Qs,out; Allen and Densmore, 2000;Armitage et al., 2013Armitage et al., , 2011;;Bonnet and Crave, 2003;Simpson and Castelltort, 2012;Tucker and Slingerland, 1997;van den Berg van Saparoea and Postma, 2008;Wickert and Schildgen, 2019), which may be recorded by changes in sedimentation rates in the deposition zone.For example, increases in either Qw or Qs,in increase Qs,out (Allen and Densmore, 2000;Armitage et al., 2013Armitage et al., , 2011;;Bonnet and Crave, 2003;Simpson and Castelltort, 2012), but each has a characteristic signature.Whereas a change in Qs,in triggers a permanent change in Qs,out, a change in Qw leads to a transient change in Qs,out (Armitage et al., 2011;Bonnet and Crave, 2003;Wickert and Schildgen, 2019).However, because environmental forcings can be cyclic rather than step changes, it can be difficult to relate variability in sedimentation rates to a distinct forcing.Moreover, changes in Qs,out in response to changes in Qs,in or Qw may be buffered, amplified, or directly transmitted through sediment routing systems (Armitage et al., 2013;Godard et al., 2013;Romans et al., 2016 and references therein;Simpson and Castelltort, 2012).We propose that to correctly interpret changes in sedimentation rates, the modifications of Qs,in within the transfer zone (then referred to as Qs) due to sediment deposition (channel aggradation) and remobilization (channel incision) must be understood.
The long temporal and broad spatial scales of fill-terrace formation and sediment deposition preclude direct observations of their potential links in nature.Numerical models provide an inroad to understand the evolution of river long profiles and/or Qs,out after perturbations (Blom et al., 2017(Blom et al., , 2016;;Malatesta et al., 2017;Simpson and Castelltort, 2012;Slingerland and Snow, 1988;Wickert and Schildgen, 2019), but most simulate river-profile evolution without taking into account modifications of channel width or terrace formation (Blom et al., 2017(Blom et al., , 2016;;Simpson and Castelltort, 2012;Slingerland and Snow, 1988;Wickert and Schildgen, 2019).In addition, most numerical models for river-profile evolution rely on equations derived for the steady-state case.As such, they may not accurately simulate transient responses, which are important for capturing terrace formation and modifications of Qs,in in the transfer zone.Physical experiments provide an alternative approach to studying the dynamics of the transfer zone, including terrace formation (Baynes et al., 2018;Frankel et al., 2007;Gardner, 1983;Lewis, 1944;Mizutani, 1998;Schumm and Parker, 1973;Wohl and Ikeda, 1997) and the evolution of Qs,out (Bonnet and Crave, 2003;Van den Berg van Saparoea and Postma 2008).Most experimental studies have tested the cutting of terraces due to base-level fall (Frankel et al., 2007;Gardner, 1983;Schumm and Parker, 1973), or explained their cutting through autogenic processes (Lewis, 1944;Mizutani, 1998).Only one experimental study by Baynes et al. (2018) investigated terrace formation as a response to changes in Qs,in or Qw, but this study focused on vertical incision into bedrock and strath-terrace cutting.Van den Berg van Saparoea and Postma (2008) and Bonnet and Crave (2003) investigated the effects of variability in Qw and Qs,in on topograpic evolution and Qs,out, but neither considered how these processes may be linked to terrace formation.To our knowledge, there are no experimental studies that consider the combined evolution of two records of landscape evolutionfill terraces in the transfer zone and sediment discharge to the deposition zonein response to environmental perturbations.
In this study, we present results from seven physical experiments of the transfer zone, represented by a single braided channel in non-cohesive sediment, in which we perturb Qw, Qs,in, and base level.We investigate the timing and geometrical response (slope, width) of the alluvial channel in the transfer zone (with a particular focus on fill-terrace cutting) and patterns and response rates of Qs,out, with a particular focus on how the records may be linked and if a combination of both records can be diagnostic of specific changes in boundary conditions.

Methods
To test the impact of different external forcings on fill-terrace formation in the transfer zone and sediment export to the deposition zone, we performed seven experiments at the Saint Anthony Falls Laboratory in Minneapolis, USA, in 2015 (Table 1).The experimental setup consisted of a wooden box with dimensions of 4 m x 2.5 m x 0.4 m (Fig. 2A) that was filled with quartz sand with a mean grain size of 144 µm.At the inlet, we supplied sand and water through a cylindrical wire-mesh diffuser filled with gravel to ensure sufficient mixing of sand and water.Water discharge (Qw) and sediment supply (Qs,in) could be regulated independently of one another.At the downstream end, water and sand (Qs,out) exited through a 20 cm-wide gap that opened onto the basin floor below.This downstream sink was required to avoid deltaic sediment deposition that would, if allowed to grow, eventually raise the base level of the upstream fluvial system.At the beginning of each experiment, we shaped an initial channel by hand (Fig. 2A) and ran the experiment under reference conditions (Qw,ref = 95 ml/s, Qs,ref = 1.3 ml/s) for 240 minutes.This runtime was sufficient to reach a quasi-steady state in which the average Qs,out approximately equaled Qs,in.After this 'spin-up' phase, the channel had a uniform equilibrium slope of approximately 7%.
Every 30 min, we stopped the experiments to measure topography using a laser scanner mounted on the railing of the basin that surrounded the wooden box.Digital elevation models (DEMs) created from the scans have a horizontal and vertical resolution of 1 mm (Fig. 2B).Using those DEMs, we measured the evolution of channel cross-sectional profiles, longitudinal channel profiles, and surface slopes.Long profiles were calculated by extracting the lowest elevation point in each crosssection at 1 mm increments.By plotting elevation against the distance down the long axis of the box rather than against channel length, resulting slopes are slightly overestimated due to the minor sinuosity of the channels.To directly compare terrace and channel slopes, we extracted 5 cm wide swath profiles along the terrace surfaces and the equivalent stretch of the modern channel.Where terraces were narrower than 5 cm, we reduced this swath width.Slopes were calculated based on a linear fit through the mean elevation profiles.To assess uncertainties, the root mean square error (RMSE) was calculated between the linear model and the observed data.
Overhead photos were taken every 20 s with a fish-eye lens (Fig. 2C).Distortions of the photos were ortho-rectified in Adobe Photoshop and photos were resampled at 1 mm horizontal resolution to directly overlap with the laser scans.Photos were turned into binary images with values of 1 for wet pixels and 0 for dry pixels.This binarization was performed by transforming the rgb (red, green, blue) images into hsv (hue, saturation, value) images and then manually defining a hue cutoff for each experiment that best separates wet and dry pixels in the image (Fig. 2D).To distinguish wet and dry pixels by color, the supplied water was dyed blue (Fig. 2C).From the binary images, the number of wet pixels in each cross-section (perpendicular to the basin margin and therefore to the average flow direction) were counted.Analyses were restricted to the areas within the orange box (Fig. 2C, D), because terraces mainly developed in this part of the channel and because we considered this sector at the upstream side of the basin to be least affected by the fixed location of the outlet.To calculate average channel width, the average number of wet pixels in 1200 cross sections were counted and are reported with one standard deviation.No overhead photos were taken for the Ctrl_1 experiment because of an error in the camera installation.
We manually measured Qs,out at 10-minute intervals by collecting the discharged sediment in a container over a 10second period and measuring its volume.This approach allowed us to estimate whether the system had returned to steady state (Qs,in ≈ Qs,out) during the runs.At the same 10-minute interval, we measured bed elevation at the inlet and at the outlet to estimate the spatially-averaged channel slope.We interpreted a constant slope for over more than 30 minutes as additional evidence for a graded (steady state) channel.The data can be found in the supplementary material.
We ran seven experiments to monitor how changes in Qs,in, Qw, and base level affect channel adjustment, the evolution of fill terraces along the main-stem and sediment discharge at the outlet (Qs,out) through time.The experiments are summarized in Table 1.To investigate the effect of Qw, we ran two separate experiments: in one experiment we doubled Qw (IQw = increase discharge) to 190 mL/s at 240 min (end of the 'spin-up' time) and in the other experiment we first halved Qw to 48 ml/s at 240 min and then returned to the initial 95 mL/s at 480 min (DQw_IQw = decrease discharge, increase discharge).To test the effect of Qs,in, we ran one experiment in which we reduced the Qs,in by 83% to 0.22 ml/s (DQs,in = decrease sediment supply) at 240 min and another one in which we first doubled Qs,in to 2.6 ml/s at 240 min and then halved Qs,in again to the initial 1.3 ml/s at 480 min (IQs,in_DQs,in = increase sediment supply, decrease sediment supply).All Qs,in and Qw changes were imposed instantaneously, resulting in a step function in the forcing (Table 1).Immediately before imposing these changes, we covered the near-channel surface with a thin layer of red sand to optically identify the area that was reworked after the change.This method allowed us to distinguish visually between fill-top (covered in red sand) and fill-cut terraces (red sand removed due to continuous overwash).We ran one experiment in which we dropped the base level by 10 cm gradually over 20 min starting at 240 min, resulting in a base-level lowering rate of 0.5 cm/min (BLF).For this experiment, we started with a base level higher than in the initial setting by flooding the basin surrounding the wooden box (Fig. 2A).The final base level equaled that of the other experiments.In this experiment, the red sand was applied immediately before the onset of base-level lowering.
Additionally, we performed two control experiments in which we made no changes to the initial conditions in to investigate whether terraces would form in our experiment without any change in external forcing (Ctrl_1, Ctrl_2).

Results
Fluvial terraces were cut in the experimental runs IQw, DQw_IQw (in the IQw phase), DQs,in, IQs,in_DQs,in (in the DQs,in phase) and BLF (Fig. 3, 4).No terraces were formed after the 'spin-up' time of Ctrl_1 and Ctrl_2.The terraces visible in the cross-section of Ctrl_2 formed in response to incision during the 'spin-up' phase and did not substantially develop after 240 min (Fig. 4B, red line).We named the terraces to the left of the channel (in downstream direction) TL and the terraces to the right TR.In all terrace-forming experiments, both fill-top (red sand) and fill-cut terraces (red sand removed) formed (Fig. 3).
Only in the IQw, DQw_IQw and IQs,in_DQs,in experiments were the fill-top terraces preserved as the most extensive terrace surface, at least on one side of the channel (Fig. 3A, B, D).In the DQs,in experiment, only a fraction of the fill-top terrace (TL) survived the transient channel adjustment phase (Fig. 3C).In all experiments that experienced upstream perturbation, fill-top and fill-cut terraces formed only in the upstream half of the sandbox.In contrast, in the BLF experiment, terraces formed in the downstream channel reach immediately after the onset of base-level drop, but were mostly destroyed within 30 min (Fig. 3E, G).Later during the BLF experiment, terraces formed in the upstream portion of the sandbox (Fig. 3F).
Fill-cut terrace cutting lagged minutes to hours behind the onset of the imposed perturbation (Fig. 3).We determined lag times from overhead photos, defined as the time interval between the onset of the perturbation (at minute 240 or 480) and the final time that the future terrace surface was occupied by water.In the two experiments during which we changed Qw and in the IQs,in_DQs,in experiment, the cutting of fill-cut terraces began within 6 minutes after the change in boundary conditions (Fig. 3A, B, D).In the IQw experiment, for example, the majority of the TL terrace is a fill-top terrace (0 min lag time) and only a small part at the downstream end was occupied until 6 minutes after perturbation (Fig. 3A, H).In the DQs,in experiment, however, several fill-cut terraces formed successively with lag times between ~14 min and 289 min (Fig. 3C).This experiment was the only one in which a sequence of terraces, instead of a single major surface, developed.In the BLF experiment, terrace cutting in the upstream part of the basin began 112 and 117 min after the onset of base-level lowering (Fig. 3F).
Fill-terrace formation requires changes in the channel-bed elevation and width of the active floodplain.In our experiments, sediment deposition (aggradation) or erosion (incision) altered the channel-bed elevation (Fig. 5).However, these changes were not uniform along the channel.In the runs Ctrl_1 and Ctrl_2, the longitudinal profiles were stable over time and experienced only minor lowering in bed elevation (max.4 cm) at their upstream ends (Fig. 5A, B).A sudden increase in Qw (IQw, and the IQw phase of DQw_IQw) or a decrease in Qs,in (DQs,in, and the DQs,in phase of IQs,in_DQs,in) both led to channel incision (Fig. 5C, D, G, H).This incision was most pronounced at the upstream end, near the changing boundary condition, but not recognizable at the downstream end (Fig. 5D, G, H), where the channel-bed elevation was fixed due to the steady base level.Sediment deposition in the channels followed a decrease in Qw (DQw phase of DQw_IQw) or an increase in Qs,in (IQs,in phase of IQs,in_DQs,in), which was, again, most recognizable at the upstream end of the channel (Fig. 5E, F).The drop in base level, however, caused maximum incision at the downstream end, and the incision wave migrated upstream as a diffuse knickzone (Fig. 5I).
Channel slope and width changes were observed in the absence of external perturbations.Channel slopes in the Ctrl_1 and Ctrl_2 marginally decreased after the 240 min 'spin-up' time from ~0.074 and 0.071, respectively, to around 0.070 and 0.067 (~6 % reduction; Fig. 6A).As such, we consider any change in slope after the 'spin-up' time that is on the same order as those observed in Ctrl_1 and Ctrl_2 as ongoing adjustment to the reference condition as opposed to the result of an external perturbation.Channel width in the control experiments varied slowly between ca.20 cm and 35 cm.
External perturbations in water and sediment inputs forced the channel width and slope to evolve.An instant doubling of Qw (IQw; Fig. 6B) resulted in a rapid decrease in channel slope that decayed exponentially as the channel approached a new graded state.After approximately 480 min, the slope was reduced from ~0.072 to ~0.043 (40% reduction), and new stable conditions were reached.The doubling of Qw also triggered an instant narrowing of the channel from ~35 cm to ~15 cm (~57 % decrease), followed by subsequent slow widening.
In contrast, suddenly reducing Qw to half its initial value (DQw_IQw; Fig. 6C) increased channel slope from ~0.072 to ~0.085 (18% increase) between 240 and 480 min runtime, and caused the channel to widen from about 25 cm to about 45 cm (~80% increase) during the same time period.The subsequent doubling in Qw back to its initial value triggered a rapid reduction in slope back to the initial ~0.072 (~15% reduction, again following an exponential decay) and an instantaneous narrowing of the channel (~45% reduction) followed by slow widening.
Reducing Qs,in by 83% (DQs,in; Fig. 6D) triggered channel slope to decrease at a slower rate than in the IQw run; and the new slope stabilized around 0.060 (24% reduction from the initial 0.079).An instantaneous decrease in channel width also occurred, but this change was again less pronounced than what we observed in the IQw experiment (~33% reduction).We detected no subsequent widening of the channel.
Finally, increasing Qs,in (IQs,in _DQs,in; Fig. 6E) led to channel steepening from a slope of 0.070 to about 0.078 (11% increase) and increased channel width from about 30 cm to about 55 cm (~83% increase).The subsequent reduction in Qs,in decreased the channel slope and caused an instantaneous channel narrowing to < 30 cm, followed by subsequent widening back to the initial width of ~30 cm.
During the base-level fall experiment (BLF; Fig. 6F), mean channel slope instantly and rapidly increased after the onset of base-level fall from about 0.047 to 0.073 (55% increase), and continued to increase at a slower rate to about 0.08, before decreasing back to 0.072.These mean values, however, average over any spatial variability in incision, meaning that they do not resolve the details of the diffusive propagation of the knickzone.Beyond impacts on slope, the drop in base level resulted in a sudden decrease in channel width, followed by three cycles of channel widening and narrowing.In summary, we observed that an increase in Qw and a decrease in Qs,in, resulted in an immediate decrease in channel slope (through upstream incision) and an instant reduction in channel width, whereas a drop in base level caused an increase in channel slope (through downstream incision) and a reduction in channel width (Fig. 6).
To analyze how well terrace surfaces represent channel slopes immediately preceding the time of perturbation, we compared the elevation profiles of the terraces on each side of the channel (yellow and orange lines) with the channel that existed at the onset of the perturbation (red line) (Fig. 7).We sampled elevations across the most extensively preserved terrace surface, regardless of its lag time, in a way that is similar to terrace mapping in the field.In experiments with increasing Qw (IQw, IQw phase of DQw_IQw) or base-level changes (BLF), the elevation profiles of the terraces are similar to the initial floodplain profile (Fig. 7A, B and E).In cases of changes in Qs,in (DQs,in, DQs,in phase of IQs,in_DQs,in), the terraces were cut at lower elevations than the former channel (Fig. 7C, D).In the DQs,in experiment, fill-cut terraces on either side of the channel formed at different elevations, with one surface about 3 cm below the other (Fig. 7C, 4E).In contrast, terrace surfaces in the other four experiments are at approximately the same elevation (Fig. 4, 7).Despite similar elevations, the slope differences between TL and TR range from about 5% (IQs,in _DQs,in) to 33% (IQw).When comparing terrace slopes to the active channel slopes at the end of each run (blue lines), terrace slopes are steeper (by 20-122%) in all experiments in which upstream conditions (Qw, Qs,in) changed (Fig. 7 A-D).In contrast, the slopes of the terraces and the active channel in the BLF experiment are similar to one another (within 11%) (Fig. 7E).
Changes in boundary conditions also affected Qs,out (Fig. 6, lowest panels).An instantaneous doubling of Qw (IQw; Fig. 6B) resulted in an instant increase in Qs,out to more than 20 times Qs,in.This rapid increase was followed by an exponential decay down to the initial Qs,out value.A sudden reduction in Qw to half its initial value (DQw_IQw; Fig. 6C) resulted in a decrease in Qs,out.The subsequent doubling in Qw back to its initial value triggered a rapid increase in Qs,out that decayed over time.In contrast, neither the instantaneous reduction in Qs,in by 83% (DQs,in; Fig. 6D) nor the doubling in Qs,in (IQs,in_DQs,in; Fig. 6E) triggered a measurable change in Qs,out.For the base-level fall experiment (BLF; Fig. 6F), Qs,out could not be measured before and during the base level drop because the basin surrounding the wooden box was flooded for this experiment.Qs,out was only measured from minute 280 onwards, which corresponds to minute 40 after the 'spin-up' of the base level fall.At that time, Qs,out was still about 10 times higher than Qs,in, and Qs,out decreased approximately linearly from that time onwards.

Discussion
When attempting to use geomorphic or depositional records to reconstruct paleo-environmental conditions, we face a range of challenges.One challenge is to understand how the information on environmental boundary conditions is translated into the sedimentary record, considering the potential modification of sediment signals during fluvial transport.A second, related challenge is that depositional records and fluvial terraces often cannot unambiguously be associated with a particular forcing mechanism.In the following, we will discuss these two challenges, and we will focus on both records that we monitored in our experimentsfill terraces in the transfer zone and sediment discharge to the deposition zone (Qs,out).Further, we will discuss the use of an integrated set of observations to address these challenges, the limitations that arise when comparing the experimental work to natural settings, as well as potential implications of our observations for future field studies.The cutting of fluvial-fill terraces requires vertical incision and a simultaneous reduction of the active floodplain width.Whether this occurs depends on the response of alluvial channels to changing boundary conditions, which can include adjustments of slope, wetted perimeter (width and depth), and/or bed-surface texture (grain-size distribution) (e.g., Blom et al., 2017;Buffington, 2012 and references therein;Wickert and Schildgen, 2019).Because the grain-size distribution in our experiments remained constant, we focus our discussion on the externally forced adjustments of channel slope (S) and width (w) during terrace formation.
In our experiments, river incision (with terrace cutting) was driven by an increase in Qw, a decrease in Qs,in, or a fall in base level (Figs. 3 -6).In the case of base-level fall, incision began at the downstream boundary and diffused upstream, producing a transient steepening.Enhanced Qw or reduced Qs,in , on the other hand, decreased channel slope.The evolution of longitudinal channel profiles in our experiments is in agreement with earlier flume studies that investigated channel response to upstream (van den Berg van Saparoea and Postma, 2008) and downstream (Begin et al., 1981;Frankel et al., 2007) perturbations, as well as with numerical models that predict the evolution of longitudinal profiles following variations in Qs,in, Qw or base level (Blom et al., 2017;Simpson and Castelltort, 2012;Wickert and Schildgen, 2019).In all experiments, incision and terrace cutting coincided with an instantaneous decrease in channel width, while aggradation corresponded to an increase in channel width (Fig. 7).
A common application of fluvial-terrace mapping is to reconstruct paleo-longitudinal channel profiles from terrace remnants (e.g., Faulkner et al., 2016;Hanson et al., 2006;Pederson et al., 2006;Poisson and Avouac, 2004).These profiles are thought to be representative of the former channel profiles, ideally reflecting their geometries immediately prior to a perturbation.However, morphological adjustments of a channel to external perturbations require time, such that the geomorphological response can lag behind the changes in environmental parameters (e.g., Blum and Tornqvist, 2000;Tebbens et al., 2000;Vandenberghe, 2003Vandenberghe, , 1995)).The lag time between an external perturbation and the onset of terrace cutting determines how much time the fluvial system has to modify the terrace sediments before their abandonment.In the following, we first discuss the relationship between lag times and the preserved terrace profiles related to upstream perturbations (Qs,in, Qw), followed by those related to downstream perturbations (BLF).
In our experiments, the upstream-perturbation-related terrace surfaces (fill-top and fill-cut terraces) following an increase in Qw had short lag times (≤ 6 min; Fig. 3A, B and 6B, C) and preserved the channel elevation profiles prior to perturbation well (Fig. 7A, B).Similarly, terrace cutting in the IQs,in_DQs,in experiment was characterized by short (TR) or no (TL) lag times (Fig. 6E).The small discrepancy between terrace slopes and initial channel slopes in this experiment (Fig. 7D) is a result of slope variations between the center of the channel belt (where initial and final channel profiles were measured), and the sides of the channel belt, where the terrace slopes were measured.In contrast, terrace cutting in the DQs,in experiment occurred with a delay of several hours and the terraces were also cut successively (Fig. 3C, 6D).The difference in lag times between the TL and TR terrace of about two and a half hours resulted in different terrace elevations on both sides of the channel, with elevation profiles several cm below the channel profile prior to perturbation (Fig. 7C).These results illustrate how short lag times are critical to enable accurate reconstructions of the pre-perturbation channel-profile, which can potentially be used to reconstruct paleo-environmental conditions.But what determines the duration of the lag time?
The length of the lag time between the perturbation and the abandonment of a terrace surface is expected to depend on the ratio of vertical incision versus lateral erosion.Bufe et al. (2018) and Malatesta et al. (2017) demonstrated that the rate of lateral channel migration scales inversely with the height of valley walls (elevation difference between a terrace surface and the active channel).As such, higher incision rates after a perturbation lead to faster wall-height growth and greater reductions in lateral mobility.Accordingly, fast incision should result in short lag times between the onset of the perturbation and terrace cutting, guaranteeing good preservation of the channel profile that existed prior to the perturbation.In contrast, slow river incision and enhanced lateral channel movement can lead to long lag times, with terrace profiles that reflect a channel profile at some (unknown) phase of adjustment.The incision rate, on the other hand, is thought to be a function of the excess sediment transport capacity, and sediment transport capacity should be directly proportional to Qw (Wickert and Schildgen, 2019): doubling Qw should correspondingly double the excess sediment transport capacity, whereas halving Qs,in should increase the excess sediment transport capacity by a factor of 0.5.Therefore, increases in Qw should, in theory, cause more rapid incision, shorter lag times, and a higher preservation potential for the pre-perturbation channel profile than a proportionately equal reduction in Qs,in,.However, while one of the two experiments with a reduction in Qs,in (DQs,in) is consistent with this theory (Fig. 6D), in the other one (IQs,in_DQs,in), we observed relatively short lag times (Fig. 6E).These unexpectedly short lag times might be related to how the incision phase was preceded by an aggradation phase (due to an increase in Qs,in).Possibly, the system rapidly settled back to the initial conditions because it had not completely adjusted to the preceding increase in Qs,in.
The lag time between the onset of base-level fall and the cutting of terraces in the upstream reach of the channel is about 115 min (Fig. 6I), which was the time required for the knickzone to propagate upstream.As such, for terraces related to base-level fall, the temporal lag between the onset of the perturbation and terrace cutting increases with upstream distance.
Hence, terrace surfaces created through upstream knickpoint migration are diachronous, become progressively younger upstream despite being a physically continuous unit.Similar results have been reported from field studies (Faulkner et al., 2016;Pazzaglia, 2013).In comparison, incision initiated near-synchronously along the entire experimental channel when triggered by a change in upstream boundary conditions (IQw, DQs,in; Fig. 5C, D).In summary, lag times between the onset of a perturbation and terrace cutting depend on the combination of local incision rates after the perturbation and the trigger for incision (base-level fall vs. a change in upstream conditions).

Terrace geometry as an indicator of perturbation type
Because fluvial-fill terraces result from changes in Qw, changes in Qs,in, or a drop in base level, their presence alone does not indicate which of the parameters changed over time.However, our experimental results revealed differences in terrace geometry between changes in upstream (Qw, Qs,in) versus downstream (BLF) conditions.For terraces related to changes in Qs,in or Qw, the slopes of terrace surfaces are always steeper than the active channel (the new steady state channel after the perturbation), whereas the slope of terraces formed due to downstream perturbations is very similar to that of the active channel (Fig. 7).These observations concur with predictions from theoretical work that suggest a positive scaling of slope and Qs,in and a negative scaling of slope and Qw, while a drop in base level should, after the signal has propagated upstream, result in a slope similar to the channel before the perturbation because of a constant Qs,in/Qw ratio (Lane, 1955;Mackin, 1948;Malatesta and Lamb, 2018;Meyer-Peter and Müller, 1948;Wickert and Schildgen, 2019;Wobus et al., 2010).Similar observations have been made in the field.In the Tien Shan, Poisson and Avouac (2004)  Our findings support earlier observations that slope comparisons between the terrace surfaces and the active channel could indicate whether an upstream-or a downstream-sourced perturbation caused the cutting of the terraces (Faulkner et al., 2016;Pepin et al., 2013;Poisson and Avouac, 2004;Wobus et al., 2010).However, such comparisons are only informative if the active channel is still graded to the boundary conditions that initiated incision and terrace cutting.In addition, this approach to identifying the terrace-formation mechanism requires negligible or quantifiable tectonic tilting of the terraces after cutting.

Qs-signal modification in the transfer zone
In steady-state, the transfer zone experiences no net sediment deposition or removal of sediment.Hence, when averaged over a certain time, Qs,in equals Qs,out and the signal can be considered as faithfully transmitted to the deposition zone (e.g.Romans et al., 2016).During geometrical channel adjustments, however, when sediment is either deposited or eroded to adjust the channel slope, Qs,out differs from Qs,in (Fig. 6).This is schematically shown in figure 8 by the divergence between the solid line (Qs,in) and circles (Qs,out).Thus, Qs-signals can be considered as modified as long as the transfer zone is in a transient state.The total time of Qs-signal modification depends on two variables: (1) the time a certain transfer zone requires to reach graded conditions again, i.e. the channel response or equilibrium time (Castelltort and Van Den Driessche, 2003;Howard, 1982;Métivier and Gaudemer, 1999;Paola et al., 1992) and ( 2) the frequency at which boundary conditions (Qw, Qs,in base level) change.Consequently, if the period of the forcing is shorter than the required response time of the channel reach, the Qs signal will never be faithfully transmitted (Paola et al., 1992).

Observable changes in sediment export to the deposition zone (Qs,out)
Regardless of whether the Qs-signal is faithfully transmitted or modified, we observed changes in Qs,out in our experiments, which would likely be reflected by changes in sedimentation rates within the deposition zone.Enhanced Qs,out, for example, was generated both by an increase in Qw (Fig. 6B) and by a drop in base level (Fig. 6F).Hence, from the Qs,out record alone, the driving mechanism cannot be identified.The temporary Qs,out peak in the IQw experiment (Fig. 6B and schematically in Fig. 8C) resembles the observations of earlier numerical (Armitage et al., 2013(Armitage et al., , 2011;;Tucker and Slingerland, 1997) and experimental work (Bonnet and Crave, 2003;van den Berg van Saparoea and Postma, 2008).In both this earlier work and ours, the peak in Qs,out was generated during the transient phase of slope adjustment.According to equation 1, an increase in Qw will decrease channel slope and, therefore, trigger river incision.Because Qs,in was held constant during the experiment, the additional sediment that reached the outlet was remobilized from within the channel, in particular from the upstream part (Fig. 5C, G; Castelltort and Van Den Driessche, 2003;van den Berg van Saparoea and Postma, 2008;Wickert and Schildgen, 2019).In contrast, a decrease in Qw requires a steeper channel slope, which is achieved through sediment deposition within the channel (Fig. 5E).In our experiments, this adjustment appears as a reduction in Qs,out relative to the upstream sediment supply during the transient slope-adjustment phase (Fig. 6C and 8D).
A decrease in Qs,in should, following the achievement of a graded channel profile, reduce Qs,out, whereas an increase in Qs,in should increase Qs,out (Allen and Densmore, 2000;Armitage et al., 2011;Bonnet and Crave, 2003).According to equation 1, reducing Qs,in will trigger temporary incision because a lower slope is required to transport less sediment with the same Qw.Conversely, increasing Qs,in without changing Qw will require a steeper transport slope and thus trigger aggradation.
Channel incision and slope reduction occurred in the DQs,in experiments (Fig. 5D, H and 6D, E), whereas aggradation and slope increase followed an increase in Qs,in (Fig. 5F and 6E).However, in none of the experiments with variable Qs,in was a clear change in Qs,out recognizable during the transient phase of slope adjustment (Fig. 6D, E and 8E, F).We consider the negative feedback between Qs,in and the bed-elevation change during the transient channel-adjustment phase as the main reason for this lack of response (Simpson and Castelltort, 2012;van den Berg van Saparoea and Postma, 2008).The additional sediment supplied upstream is deposited within the channel, resulting in aggradation, and is therefore not detectable at the outlet.When less sediment is supplied upstream, the channel incises and complements the supplied upstream sediment with remobilized sediment from within the channel, such that once again, no change in Qs,out is detectable at the outlet during the adjustment phase.We did not run the experiments long enough to analyze the adjusted steady-state phase, but we would expect that once the channel has adjusted to new equilibrium conditions, Qs,out will eventually equal Qs,in (Fig. 8E, F; Allen and Densmore, 2000;Armitage et al., 2011;Bonnet and Crave, 2003;Wickert and Schildgen, 2019).
Internal channel dynamics can lead to variability in Qs,out even without external forcing.In the Ctrl_1 and Ctrl_2 experiments, scatter in the Qs,out signal was up to 5 times the value of Qs,in (Fig. 6a).This variability is due to continuous lateral movement of the channel and subsequent bank collapse, which results in stochastic contributions of additional sediment.
Lateral channel mobility of a stream varies with water and sediment discharge (Bufe et al., 2018;Wickert et al., 2013).
However, if the volume of sediment mobilized from valley walls due to lateral migration is much larger than the change in Qs,in, then no clear signal in Qs,out might be recognizable, even after channel adjustment.The channel instead will continually adjust to the stochastic lateral input of sediment.
Regarding Qs,out signals, we conclude that terraces, floodplains, and the channel itself act as a temporary storage space where sediment can be deposited or remobilized when boundary conditions change (Coulthard et al., 2005;Simpson and Castelltort, 2012;van den Berg van Saparoea and Postma, 2008).The consequence of sediment deposition or remobilization during transient response times is that Qs,in differs from Qs,out, such that the Qs-signal can be considered as modified during transient phases of channel adjustment.Our data also support earlier findings by Simpson and Castelltort (2012) and van den Berg van Saparoea and Postma ( 2008), who concluded from their respective numerical model and physical experiments that Qw variability creates an amplified, substantial response in Qs,out, whereas changes in Qs,in create a dampened response in Qs,out due to the a negative feedback between Qs,in and channel slope.Our experiments, illustrated schematically in Fig. 8, also suggest that Qw-driven Qs,out changes are temporary, and that as the channel slope adjusts to the new input Qw, Qs,out evolves back to its initial steady-state value.In contrast, Qs,out signals driven by changes in Qs,in may not be observable during transient channel adjustment, but will occur and persist once the channel has adjusted to new steady-state conditions.

Combining the records
Both fill terraces in the transfer zone and changes in sedimentation rates in the deposition zone record changes in boundary conditions.Our experiments provide an opportunity to investigate the links between terrace formation and sediment export to the deposition zone, as well as how both of these records, if available, can disambiguate changes in past Qs,in, Qw, or base level.
Sediment discharge at the outlet (Qs,out) is a function of (1) upstream sediment supply (Qs,in) and ( 2) sediment deposition or mobilization within the transfer zone.Doubling Qw while holding Qs,in constant triggers river incision to archive a lower channel slope.This incision leaves behind terraces that are steeper than the modern, graded channel (Fig. 9).The sediment mobilized in the transfer zone during the transient incision phase produces a transient peak in Qs,out.Reducing Qs,in also reduces the equilibrium transport slope, causing channel incision and terrace abandonment.However, no peak in Qs,out during the transient phase is visible, as the extra sediment remobilized from the transfer zone is compensated by the preceding reduction in Qs,in (Fig. 9).Finally, incision due to a fall in base level also generates a temporary peak in Qs,out due to the additional sediment remobilized within the transfer zone.The terraces left behind following this base-level-driven incision, however, have slopes parallel to that of the modern, graded channel.In summary, Qs,out reflects a combination of Qs,in and the geometrical adjustment of the transfer zone, which in turn is recorded by fill terraces.
Coupling the two records, if available, provides the opportunity to unambiguously identify the forcing mechanism, which is not possible using either the fill terraces or the deposits alone (Fig. 9).For example, the presence of terraces whose slopes are steeper than the present-day channel, combined with a simultaneous but transient peak in Qs,out, points towards a change in Qw as the main driver.In turn, terraces whose slopes are steeper than the main channel in combination with a lagged reduction in Qs,out point towards a change in Qs,in as the main driver.Finally, a temporary increase in Qs,out in combination with channel-parallel terraces that young in the upstream direction indicates past base-level fall.Complications may arise when the forcing includes a combination of changes in Qs,in and Qw.Nevertheless, our results point to the potential of combining terrace records with sediment-export data for the reconstructions of paleo-environmental conditions.

Limitations of experiments
Physical experiments allow for investigations of the isolated influence of individual key parameters on landscape evolution.However, a number of limitations arise when attempting to compare the experimental results to natural settings.
First of all, in natural sediment-routing systems, the three distinct zones of erosion, transfer and deposition (Fig. 1) are coupled to one another (Allen, 2017).Erosion processes on the hillslopes, for example, determine the amount of sediment provided to the transfer zone, i.e.Qs,in (e.g., Dixon et al., 2009;Tofelde et al., 2018).In turn, changes in channel-bed elevation in the transfer zone can affect hillslope-erosion processes (e.g., Hurst et al., 2012;Roering et al., 2007).In our experimental setup, however, we investigate the response of the transfer zone as an isolated feature and can thus not account for any hillslopechannel feedbacks that might lead to additional variations in sediment supply to the channel.
Second, we varied Qw and Qs,in separately, forcing them to remain independent of one another.In natural systems, however, they are commonly coupled.For example, changes in precipitation can alter both Qw and Qs,indirectly though changes in rainfall-driven sediment-transport rates from hillslopes to the channel (Bookhagen et al., 2006;Dey et al., 2016;Steffen et al., 2010Steffen et al., , 2009) ) and indirectly through long-term changes in hillslope-stabilizing vegetation types (Garcin et al., 2017;Langbein and Schumm, 1958;Schmid et al., 2018;Torres Acosta et al., 2015;Werner et al., 2018).Those feedback mechanisms between different key parameters (Qw, Qs,in) and between sub-zones of sediment-routing systems will likely complicate the forcing-response behavior of natural systems.
Third, we have only investigated a single, braided channel.Therefore, our experimental set-up does not allow us to investigate channel-geometry adjustments related to feedbacks between the main stem and adjacent tributaries (Schumm, 1979(Schumm, , 1973)), and related terraces forming at channel junctions (Faulkner et al., 2016;Larson et al., 2015;Schildgen et al., 2016).
Furthermore, we can draw no conclusions on terraces forming in meandering rivers, such as those related to meander-bend cut-off (e.g., Erkens et al., 2009;Gonzalez, 2001;Limaye and Lamb, 2016).As such, the lack of terrace formation in the two control experiments after the 'spin-up' time does not imply that autogenic terraces do not exist in natural systems, because several potential mechanisms of autogenic or complex-response terrace formation like meander-bend cut-off (Erkens et al., 2009;Gonzalez, 2001;Limaye and Lamb, 2016;Womack and Schumm, 1977) or internal feedbacks between the main-stem and tributaries (Schumm, 1979, 1973, Gardener 1983, Schumm and Parker 1973, Slingerland and Snow 1988) could not be tested with our experimental set-up.
Fourth, apart from the step changes in input parameters, Qs,in and Qw were held constant through time.As the experiments exhibit geomorphically effective flow conditions at all times (intermittency equals 1), we assume that the experiments integrate over a number of large floods in natural channels.Natural rivers in turn, experience a wide range of intermittencies.This variability in natural systems complicates any attempts to scale channel response times and lag times from experiments to real systems, but also complicates the comparison of real systems with each other.
Finally, we performed a limited number of experimental runs, with only the control experiments being repeated.
Although we did not repeat the experiments that included external perturbations, we consider the last phase of the two experiments during which we performed two changes (DQw_IQw and IQs,in_DQs,in) as repetition of the experiments with only one perturbation (IQw and DQs,in), albeit with different absolute values of Qw and Qs,in.Comparing those experiments reveals similar trajectories of channel evolution (longitudinal profiles, slope, width; Fig. , 5 and 6).In addition, the same boundary conditions (Qs,in, Qw) persisted during the 'spin-up' phases as well as at the end of the two experiments during which we performed two changes.During those conditions, the channel slopes always evolved to a value of ~0.07.Although not being exact repetitions of the same experiments, the evolution to the same equilibrium conditions indicates that the results are reproducible.But, we acknowledge that further repetitions would improve our ability to quantify variability that is internal to each system.
Despite the above limitations, the performance of physical experiments under controlled conditions allows us to directly link causes and their effects in fluvial systems.As the key parameters (Qs,in, Qw, base level) can be varied independently, physical experiments provide the opportunity to isolate the influence of different environmental parameters on the evolution of landscapes and to test theoretical models (Paola et al., 2001).As many processes acting in landscapes are scale-independent, experimental observations can help to decipher process behavior in natural systems (e.g., Cantelli et al., 2004).

Implications for field studies
Despite the restrictions when comparing the experimental work to natural systems, the general patterns that we observed do have implications for field studies.First, lag times between the perturbation and the onset of terrace cutting can be important when dating the surfaces of fluvial-fill terraces in the field.Common methods to date the onset of river incision include the dating of terrace surface material with cosmogenic exposure dating (e.g., Schildgen et al., 2016;Tofelde et al., 2017), dating sand or silt lenses with optically stimulated luminescence close to the terrace surface (OSL; e.g., Fuller et al., 1998;Schildgen et al., 2016;Steffen et al., 2009) or dating embedded organic material with 14 C (Farabaugh and Rigsby, 2005;Scherler et al., 2015).When transferring our observations to a field scenario, the ~2h or more of channel material reworking before terraces were cut within the upstream part of the channel in the BLF and the DQs,in experiment would result in terrace ages that are younger than the time of perturbation.The best temporal correlations between the perturbation and the terrace surface ages are achieved by those formed by changes in Qw due to the fast onset of vertical incision and minimal reworking of terrace surface material.To assess the significance of this time lag in natural systems requires more work on how to scale the experiment to larger channels.
Second, in tectonically active regions, both strath and fill terraces have been used to calculate river incision rates to infer tectonic uplift rates (e.g., Litchfield and Berryman, 2006;Maddy et al., 2001;Schildgen et al., 2012;Wegmann and Pazzaglia, 2009).Because the slope changed in our experiments of upstream perturbations, incision rates were not uniform along the channel (Fig. 5).Litchfield and Berryman (2006) also measured variable fluvial incision rates based on terrace heights at several locations along 10 major rivers located along the Hikurangi Margin, New Zealand.Accordingly, incision rates reconstructed from terraces should be interpreted in the context of potential changes to the shape of the longitudinal profile in addition to tectonic changes.
Third, despite information stored in variable sediment deposition rates in the transfer or deposition zone, information on landscape evolution can also be preserved in the chemical composition of the deposited sediment.Our findings have implications for geochemical signatures of sediment, for example the concentration of cosmogenic 10 Be, which is commonly measured to infer catchment mean denudation rates (Bierman and Steig, 1996;Brown et al., 1995;Granger et al., 1996).In cases of channel aggradation, Qs,out is reduced compared to Qs,in due to deposition within the channel (Fig. 8B, D, F).The exported sediment could be sourced from incoming sediment that is not deposited (grey circles) and/or mixes with remobilized sediment within the channel (yellow circles).In general, net deposition along the channel leads to the majority of the grains at the outlet being freshly delivered from local hillslopes, thus carrying the contemporaneous catchment chemical composition at the time of transport.In contrast, during incision, older material stored within the channel, floodplain, and/or terraces is remobilized and contributes to the temporary peak in Qs,out (Fig. 8A, C, E).Shortly after the perturbation, most of the remobilized sediment will be stratigraphically high and relatively young (yellow circles), but older material from deeper layers (orange and red circles) will be progressively remobilized and mixed with young material from upstream.Cosmogenic nuclide analyses along the eastern Altiplano margin (Hippe et al., 2012) and in the Amazon basin (Wittmann et al., 2011) indicate that sediment can be stored within the fluvial system over thousands to millions of years.Remobilization of formerly deposited material and subsequent mixing with fresh hillslope material (incoming sediment) can temporally modify signals stored in the geochemical composition of detrital river sediments (e.g., Tofelde et al., 2018;Wittmann et al., 2016Wittmann et al., , 2011)).The degree of modification is thought to be a function of the ratio between fresh and remobilized material exported at a certain time as well as of the difference in geochemical composition between the fresh and remobilized material.We conclude that modern chemical signals are more likely to be transmitted through the system during aggradation phases, whereas local sediment that has been transiently stored may strongly overprint the signal of modern sediments during times of incision.

Summary and Conclusion
We performed seven physical experiments to investigate the effects of changing boundary conditions (Qs,in, Qw, base level) on channel geometry and related fill-terrace cutting as well as on sediment discharge (Qs,out).To reliably reconstruct paleo-environmental conditions from terraces in the transfer zone or sedimentary deposits in the sedimentation zone, it is important to understand (1) how information on environmental conditions may be modified and eventually transferred into the geologic record, and (2) whether the geomorphic characteristics of terraces or the patterns of sedimentation rates are specific to the forcing mechanism.
Our experiments provide important insights into the coupling between sediment discharge to the deposition zone (Qs,out) and terrace formation in the transfer zone.Sediment discharge to the deposition zone (Qs,out) is a combination of sediment supply to the transfer zone (Qs,in) and its modification within the transfer zone through deposition or remobilization.
Deposition or remobilization of sediment within the transfer zone occurs mainly during the transient response phase after a perturbation.Hence, during the transient response, Qs,out does not equal Qs,in.One consequence of this finding for field studies is that the geochemical composition of sediment sampled during transient river incision does not represent direct hillslope inputs alone, but rather includes a significant component of remobilized sediment from some point in the geologic past.This remobilization can modify the 10 Be concentration in detrital sediment, thus impacting estimates of catchment-wide denudation rates.
The same modifications (sediment deposition and remobilization) that alter the Qs signal during transient times also form fill terraces in the transfer zone.Increases in Qw trigger channel incision to archive a lower equilibrium slope.The resulting temporary peak in Qs,out coincides with the cutting of terraces whose slopes are steeper than the main channel.
Reducing Qs,in also reduces the equilibrium transport slope, causing channel incision and terrace abandonment.However, no substantial increase in Qs,out occurs during the transient phase, as the extra sediment remobilized from the transfer zone is compensated by the preceding reduction in Qs,in.. Finally, a drop in base level causes a temporary peak in Qs,out and the formation of terraces parallel to the modern channel.
If both terrace and depositional records are available, the combination of the two can unambiguously identify the main mechanism driving channel adjustment.Identifying this mechanism can be important, for example, when using the height of terraces to infer channel incision rates.As upstream perturbations cause greater incision at the upstream end than at the downstream end, incision rates inferred from terrace heights are expected to vary along the profile.
The cutting of terraces following an upstream perturbation (Qs,in, Qw) requires a period of time (lag time) that is a function of incision rate, which in turn is expected to be a function of the excess transport capacity of a channel (Wickert and Schildgen, 2019).Indeed, our experiments showed that greater excess transport capacity leads to faster incision and shorter lag times, thus ensuring better preservation of the channel profile that existed prior to perturbation.These lag times can also be critical for field studies that attempt to link the ages of terraces surfaces to the timing of perturbations, as long lag times may lead to substantial temporal mismatches between forcing and fluvial-system response.

4. 1
Terrace formation in the transfer zone 4.1.1Conditions of terrace formation, lag times, and the preservation of pre-perturbation channel profiles found a successive reduction in slope within a terrace sequence, which they related to changes in Qw.In the Central Andes,Pepin et al. (2013) explained downstream-converging terraces (and thus a reduction in terrace-surface slopes) on a piedmont through variability in climatic drivers.In contrast, along the Chippewa River in the USA (a tributary to the Mississippi River), where terrace cutting is linked to base-level fall,Faulkner et al. (2016) found no substantial slope change between the longitudinal profile reconstructed from the terraces and the modern channel.

Figure 1 .
Figure1.Schematic summary of a sediment-routing system and records of landscape evolution.Sediment-routing systems are typically subdivided into three zones: sediment production, sediment transfer, and sediment deposition(Allen, 2017;Castelltort and Van Den Driessche, 2003).As sediment production is thought to vary with environmental conditions, changes in those conditions might be preserved in sedimentary records in the transfer zone (e.g., fill terraces) or deposition zone (e.g., sedimentation rates).Complications arise, however, as alluvial rivers within the transfer zone continuously adjust their channel geometry to the incoming water discharge (Qw) and 820

Figure 2 .
Figure 2. Experimental setup, data collection, and analysis.(A) Overview of experimental setup.Sediment supply (Qs,in) and water discharge (Qw) can be regulated separately.For all but the base-level fall (BLF) experiment, the base level was fixed.Water and sediment fell off of an edge at the outlet.For the BLF experiment (shown in the picture), the base level was controlled through the water level in the surrounding basin.(B) Digital elevation model (DEM) derived from laser scans showing the final topography of the increased water (IQw) experiment.(C) Overhead photograph of the IQw experiment taken directly before the scan shown in B. The surface was covered with a thin layer of red sand before the instant increase in Qw was performed.The remnants of red sand on the terraces indicate no further reworking after the onset of increased discharge.(D) Overhead photographs were turned into binary (wet, dry) images from which the average channel width within the analyzed area (orange frame) can be calculated.

Figure 3 .
Figure 3. Fill terraces formed during experimental runs.(A-D) Top-view of terraces that formed due to upstream perturbations.Remnants of red sand on the terrace surfaces indicate areas that have not been flooded after the change in boundary condition was performed (i.e.fill-top terrace).The other terrace surfaces were cut with the indicated lag times (fill-cut terraces).(E-F) During the base-level fall (BLF) experiment, terraces at the downstream end were abandoned instantly after the onset of base-level fall (250 min = 10 min after onset of BLF).840

Figure 4 .
Figure 4. Evolution of cross sections in the upper part of the channel (left panel).In each cross section, the lowest point is set

Figure 5 .Figure 6 .
Figure 5. Evolution of longitudinal river profiles from minute 240 (end of 'spin-up' phase) onwards.River profiles were extracted from the laser scans.Laser scans were recorded every 30 min, and an additional two scans at 10 and 20 minutes after the initiation of the base-level fall were conducted during the BLF experiment.Dashed arrows indicate down-basin distance along which terraces formed.Solid arrows indicate modes of aggradation or incision.Note that the DQw_IQw and IQs,in_DQs,in were split into two panels each, with one panel representing each phase.855

Figure 7 .Figure 8 .Figure 9 .
Figure 7. Elevation profile and slope comparison of terrace surfaces and active channels.Elevation profiles are given as mean (solid lines) and minimum and maximum values (dashed lines), extracted along a 5 cm wide swaths as indicated on the right panel.Swath width was reduced in two cases of too narrow terraces to 1 cm (DQw_IQw TL terrace) and 2 cm (DQs,in TL terrace).TL and TR indicate terraces 875