Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment

Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd.     

Geometric constraints on composition of sediment derived from erosional landscapes

Although unroofing sequences are well known in the stratigraphic record, there is no general theory for estimating relevant basic quantities such as the time history of sediment production from a particular unit or the degree of mixing between successive units. Here we investigate the production of sediment from layered source rocks that are milled off by steady-state erosional topography. The shape of the sediment-production function for milling off a thin horizontal layer is given by the derivative of the hypsometric function, in the form of area contained within contours as a function of contour altitude. The time-scale for the production function, the ‘topographic mixing time’, is set by the topographic relief divided by the uplift rate. The production function for a sharp transition from one unit to another is given directly by the hypsometric function. The effects of stratal dip parallel to the mean slope of the erosional topography and finite layer thickness can be accounted for to a first approximation by simple geometric corrections to the mixing time. Finite layer thickness also has the effect of smoothing the production function although most natural hypsometric functions are smooth enough that this effect is relatively weak. The quality of an unroofing sequence can be measured in terms of the ‘sharpness’ of separation of successive peaks in sediment production produced by milling off a sequence of geometrically similar layers. This peak sharpness can be parameterized by a ratio of the interval between successive peaks in sediment production to topographic mixing time. By this measure, the quality of unroofing sequences is controlled by two parameters: the ratio of layer thickness to topographic relief, and the dip angle. The dip angle in concert with topographic mixing exerts a strong control on the degree of signal segregation; in particular, production of cleanly segregated signals for dip angles greater than about 15° requires very high ratios of layer thickness to relief. Hence identification of distinct unroofing sequences may place significant and useful constraints on the attitude and/or thickness of units in the eroding stratigraphy.     

Steering of experimental channels by lateral basin tilting

A major issue in tectonics and sedimentation is the role of cross‐stream tectonic tilting in steering channels. The general idea is that channels will be attracted to lateral maxima in subsidence rate. A physical experiment performed in 1999 at the St. Anthony Falls Laboratory, however, was in conflict with the idea and showed that fluvial channels and resulting stratigraphy can be insensitive to even relatively strong lateral variation in subsidence. Here, we present results from an experiment which uses a simplified relay‐ramp geometry with laterally variable uplift and subsidence to test a hypothesis developed from the earlier experiment: Tectonic tilting steers channels only when the ratio of the time scales describing lateral channel mobility to tectonic deformation is sufficiently large. Occupation time by experimental channels and sand fraction in the deposit (a proxy for channel deposition) both increase with subsidence rate indicating strong steering of channels by tectonic forcing. We also found that, due to local incision, uplift lengthened the time scale for lateral channel migration relative to subsidence. Comparing channel mobility at the beginning of the experiment, with no tectonic forcing, to later tectonic stages of the experiment indicates that active tectonics increased the channel time scale. The interplay of channel steering with uplift and subsidence led to cyclic appearance and disappearance of an autogenic lake in the hanging‐wall basin. This lake was associated with alternation between channels going around vs. across the adjoining upstream uplifted footwall region. This creation and filling of the lake under constant tectonic forcing (constant fault slip rate) in the hanging wall created subaerial fan‐delta parasequences separated by fluvial deposits.     

Reliability of the automatic procedures for locating earthquakes in southwestern Alps and northern Apennines (Italy)

Reliable automatic procedure for locating earthquake in quasi-real time is strongly needed for seismic warning system, earthquake preparedness, and producing shaking maps. The reliability of an automatic location algorithm is influenced by several factors such as errors in picking seismic phases, network geometry, and velocity model uncertainties. The main purpose of this work is to investigate the performances of different automatic procedures to choose the most suitable one to be applied for the quasi-real-time earthquake locations in northwestern Italy. The reliability of two automatic-picking algorithms (one based on the Characteristic Function (CF) analysis, CF picker, and the other one based on the Akaike’s information criterion (AIC), AIC picker) and two location methods (“Hypoellipse” and “NonLinLoc” codes) is analysed by comparing the automatically determined hypocentral coordinates with reference ones. Reference locations are computed by the “Hypoellipse” code considering manually revised data and tested using quarry blasts. The comparison is made on a dataset composed by 575 seismic events for the period 2000–2007 as recorded by the Regional Seismic network of Northwestern Italy. For P phases, similar results, in terms of both amount of detected picks and magnitude of travel time differences with respect to manual picks, are obtained applying the AIC and the CF picker; on the contrary, for S phases, the AIC picker seems to provide a significant greater number of readings than the CF picker. Furthermore, the “NonLinLoc” software (applied to a 3D velocity model) is proved to be more reliable than the “Hypoellipse” code (applied to layered 1D velocity models), leading to more reliable automatic locations also when outliers (wrong picks) are present.     

A model for the viscosity of rhyolite as a function of H2O-content and pressure: A calibration based on centrifuge piston cylinder experiments

The Newtonian viscosity of synthetic rhyolitic liquids with 0.15-5.24 wt% dissolved water was determined in the interval between 580 and 1640 °C and pressures of 1 atm and 5-25 kbar. Measurements were performed by combining static and accelerated (up to 1000g) falling sphere experiments on water-bearing samples, with high temperature concentric cylinder experiments on 0.15 wt% H₂O melts. These methods allowed viscosity determinations between 10² and 10⁷ Pa s, and cover the complete range of naturally occurring magmatic temperatures, pressures, and H₂O-contents for rhyolites.Our viscosity data, combined with those from previous studies, were modeled by an expression based on the empirical Vogel-Fulcher-Tammann equation, which describes viscosities and derivative properties (glass transition temperature T_g, fragility m, and activation volume of viscous flow Va) of silicic liquids as a function of P-T-X_(H2O). The fitted expressions do not account for composition-dependent parameters other than X_(H2O) and reproduce the entire viscosity database for silicic liquids to within 3.0% average relative error on log η (i.e. std. error of estimate of 0.26 log units).The results yield the expected strong decrease of viscosity with temperature and water content, but show variable pressure dependencies. Viscosity results to be strongly affected by pressure at low pressures; an effect amplified at low temperatures and water contents. Fragility, as a measure for the deviation from Arrhenian behavior, decreases with H₂O-content but is insensitive to pressure. Activation volumes are always largely negative (e.g., less than −10 cm³/mol) and increase strongly with H₂O-content. Variations in melt structure that may account for the observed property variations are discussed.     

Flexural deformation controls on Late Quaternary sediment dispersal in the Garo-Rajmahal Gap,NW Bengal Basin

Subsurface deformation is a driver for river path selection when deformation rates become comparable to the autogenic mobility rate of rivers. Here we combine geomorphology, soil and sediment facies analyses, and geophysical data of the Late Quaternary sediments of the central Garo-Rajmahal Gap in Northwest Bengal to link subsurface deformation with surface processes. We show variable sedimentation characteristics, from slow rates (<0.8 mm/year) in the Tista megafan at the foot of the Himalaya to nondeposition at the exposed surface of the Barind Tract to the south, enabling the development of mature soils. Combined subsidence in the Tista fan and uplift of the Barind Tract are consistent with a N-S flexural response of the Indian plate to loading of the Himalaya Mountains given a low value of elastic thickness (15–25 km). Provenance analysis based on bulk strontium concentration suggests a dispersal of sediment consistent with this flexural deformation—in particular the abandonment of the Barind Tract by a Pleistocene Brahmaputra River and the current extents of the Tista megafan lobes. Overall, these results highlight the control by deeply rooted deformation patterns on the routing of sediment by large rivers in foreland settings.     

The Vema Transverse Ridge (Central Atlantic)

Transverse ridges are elongate reliefs running parallel and adjacent to transform/fracture zones offsetting mid-ocean ridges. A major transverse ridge runs adjacent to the Vema transform (Central Atlantic), that offsets the Mid-Atlantic Ridge by 320 km. Multibeam morphobathymetric coverage of the entire Vema Transverse ridge shows it is an elongated (300 km), narrow (<30 km at the base) relief that constitutes a topographic anomaly rising up to 4 km above the predicted thermal contraction level. Morphology and lithology suggest that the Vema Transverse ridge is an uplifted sliver of oceanic lithosphere. Topographic and lithological asymmetry indicate that the transverse ridge was formed by flexure of a lithospheric sliver, uncoupled on its northern side by the transform fault. The transverse ridge can be subdivided in segments bound by topographic discontinuities that are probably fault-controlled, suggesting some differential uplift and/or tilting of the different segments. Two of the segments are capped by shallow water carbonate platforms, that formed about 3–4 m.y. ago, at which time the crust of the transverse ridge was close to sea level. Sampling by submersible and dredging indicates that a relatively undisturbed section of oceanic lithosphere is exposed on the northern slope of the transverse ridge. Preliminary studies of mantle-derived ultramafic rocks from this section suggest temporal variations in mantle composition. An inactive fracture zone scarp (Lema fracture zone) was mapped south of the Vema Transverse ridge. Based on morphology, a fossil RTI was identified about 80 km west of the presently active RTI, suggesting that a ridge jump might have occurred about 2.2 m.a. Most probable causes for the formation of the Vema Transverse ridge are vertical motions of lithospheric slivers due to small changes in the direction of spreading of the plates bordering the Vema Fracture Zone.     

Seismic response of L��Aquila downtown from comparison between 2D synthetics spectral ratios of SH, P-SV and Rayleigh waves and observations of the 2009 earthquake sequence

This is the first part of a study on the seismic response of the L’Aquila city using 2D simulation and experimental data. We have studied two velocity-depth models with the aim of outlining the behavior of a velocity reversal in the top layer, which is associated with the stiff Brecce de L’Aquila unit (BrA). In this setting, the SMTH model is topped by a layer with about 2:1 impedance contrast with the underlying layer while the NORV model has no velocity reversal. We have simulated the propagation of SH and P-SV wavefields in the range 0–10 Hz for incidence 0°–90°. Earthquake spectral ratios of the horizontal and vertical components at six sites in L’Aquila downtown are compared to corresponding synthetics spectral ratios. The vertical component of P-SV synthetics enables us to investigate a remarkable amplification effect seen in the vertical component of the recorded strong motion. Sites AQ04 and AQ05 are best matched by synthetics from the NORV model while FAQ5 and AQ06 have a better match with synthetics spectral ratios from the SMTH model. All simulations show this behavior systematically, with horizontal and near-horizontal incident waves predicting the overall pattern of matches more clearly than vertical and near-vertical incidence. The model inferences are in agreement with new geological data reporting lateral passages in the top layer from the stiff BrA to softer sediments. Matches are good in terms of frequency of the first amplification peak and of spectral amplitude: the horizontal components have spectral ratio peaks predominantly at 0.5 Hz in the simulations and at 0.7 Hz in the data, both with amplitudes of 4, while the vertical component spectral ratios reach values of 6 at frequencies of about 1 Hz in both data and simulations. The vertical component spectral ratios are very well matched using Rayleigh waves with incidence at 90°. The NORV model without the velocity reversal predicts spectral ratio peaks for the horizontal components at frequencies up to 6 Hz. The reversal of velocity acts as a low-pass frequency filter on the horizontal components reducing the amplification effect of the sediment filled valley.