Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

Topographic forcing of flow partition and flow structures at river bifurcations

This paper presents the predicted flow dynamics from the application of a Reynolds-averaged NavierStokes model to a series of bifurcation geometries with morphologies measured during previous flume experiments. The topography of the bifurcations consists of either plane or bedform-dominated beds which may or may not possess discordance between the two bifurcation distributaries. Numerical predictions are compared with experimental results to assess the ability of the numerical model to reproduce the division of flow into the bifurcation distributaries. The hydrodynamic model predicts: (1) diverting fluxes in the upstream channel which direct water into the distributaries; (2) super-elevation of the free surface induced at the bifurcation edge by pressure differences; and (3) counter-rotating secondary circulation cells which develop upstream of the apex of the bifurcation and move into the downstream channels, with water converging at the surface and diverging at the bed. When bedforms are not present, weak transversal fluxes characterize the upstream channel for almost its entire length, associated with clearly distinguishable secondary circulation cells, although these may be under-estimated by the turbulence model used in the solution. In the bedform dominated case, the same hydrodynamic conditions were not observed, with the bifurcation influence restricted and depth scale secondary circulation cells not forming. The results also demonstrate the dominant effect bed discordance has upon flow division between the two distributaries. Finally, results indicate that in bedform dominated rivers. Consequently, we suggest that sand-bed river bifurcations are more likely to have an influence that extends much further upstream and have a greater impact upon water distribution. This may contribute to observed morphological differences between sand-bedded and gravel-bedded braided river networks. Copyright (C) 2012 John Wiley & Sons, Ltd.

Reservoir compensation releases : impact on the macroinvertebrate community of the Derwent River, Northumberland, UK-A longitudinal

River regulation for the purposes of public water supply causes the flow regime downstream of a dam to change. Traditionally, in the UK, such regulation was accompanied by requirements for reservoir releases to compensate downstream water users (e.g. industry) for the loss of natural flow (compensation flows). In this article, we compare a unique pre-impoundment macroinvertebrate data set for a regulated upland river with survey data post-impoundment. This allows a longitudinal assessment of the response of the system to regulation. The Derwent River, Northumberland, was impounded in 1966. Impacts on the hydrological regime were quantified by comparing long-term hydrographs, flow duration curves, flow ranges and flashiness indices for the pre-impoundment and post-impoundment periods. The comparison of changes in macroinvertebrate richness and diversity post-impoundment showed that the change in flow regime has had limited effect on the ecological community structure. The flow regime of the Derwent River has become less flashy with fewer extreme events, and the richness and the diversity of macroinvertebrates have, in some cases, increased and at worst have not deteriorated. We suggest that this reflects the strict compensation regime, which has guaranteed minimum flows at all times. Copyright (c) 2012 John Wiley & Sons, Ltd.

The relationship between Lamb weather types and long-term changes in flood frequency, River Eden, UK

Research has found that both flood magnitude and frequency in the UK may have increased over the last five decades. However, evaluating whether or not this is a systematic trend is difficult because of the lack of longer records. Here we compile and consider an extreme flood record that extends back to 1770. Since 1770, there have been 137 recorded extreme floods. However, over this period, there is not a unidirectional trend of rising extreme flood risk over time. Instead, there are clear flood-rich and flood-poor periods. Three main flood-rich periods were identified: 18731904, 19231933, and 1994 onwards. To provide a first analysis of what is driving these periods, and given the paucity of more sophisticated datasets that extend back to the 18th century, objective Lamb weather types were used. Of the 27 objective Lamb weather types, only 11 could be associated with the extreme floods during the gauged period, and only 5 of these accounted for > 80% of recorded extreme floods The importance of these five weather types over a longer timescale for flood risk in Carlisle was assessed, through calculating the proportion of each hydrological year classified as being associated with these flood-generating weather types. Two periods clearly had more than the average proportions of the year classified as one of the flood causing weather types; 19001940 and 19832007; and these two periods both contained flood-rich hydrological records. Thus, the analysis suggests that systematic organisation of the North Atlantic climate system may be manifest as periods of elevated and reduced flood risk, an observation that has major implications for analyses that assume that climatic drivers of flood risk can be either statistically stationary or are following a simple trend. Copyright (c) 2011 Royal Meteorological Society

Characterization of thirteen microsatellite markers in river and brook lampreys (Lampetra fluviatilis and L-planeri)

We describe the development based on 454 pyrosequencing technology of thirteen microsatellite markers for two closely related species of lamprey: Lampetra fluviatilis and L. planeri. The number of alleles per locus ranged from 2 to 5 in L. fluviatilis and from 2 to 6 in L. planeri. Gene diversity ranged from 0.062 to 0.718 in L. fluviatilis and from 0.322 to 0.677 in L. planeri. These markers will be helpful to study population genetic structure of both species and resolve their taxonomic status as separate species or ecotypes of a single species.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.

Mixing of Rhône River water in Lake Geneva (Switzerland-France) inferred from stable hydrogen and oxygen isotope profiles

Depth profiles were sampled at different locations throughout Lake Geneva on a monthly and seasonal basis over the course of 2 years and analysed for their stable hydrogen and oxygen isotope compositions. The isotopic compositions indicate an isotopic stratification in the metalimnion during summer and fall. This is related to mixing of Rhône River water, which in summer is dominated by snow and glacier melt waters, and lake water, with the latter having a homogenous isotopic composition. The observed interflow layer is 7-15 m thick and can be traced by the distinct stable isotope composition of the water for about 55 km throughout the lake as well as into shallow bay regions. Depth of the interflow layer close to the Rhône River mouth is similar to those previously described based on echo-soundings and turbidity profiles of sediment dispersion. In contrast to previous descriptions of the interflow within Lake Geneva, the stable isotope compositions allow for direct, natural tracing of the Rhône River water even in cases where the turbidity and conductivity measurements do not indicate such an interflow. In addition, the method allows for a quantification of the Rhône River and lake water in the interflow with the fraction of Rhône River water within the interflow estimated to be up to 37% in summer. The isotopic composition further indicates different vertical mixing processes within the two lake basins of Lake Geneva, related to the density gradients and local stability within the water column. The method may be applicable to other lakes in catchments with large differences in the topography as water that originates from high altitudes or glaciers has a distinct oxygen and hydrogen isotope composition compared to other sources of water originating at lower altitudes and/or from direct precipitation over the lake. Stable isotope measurements thus improve the understanding of the circulation of water within the lake, which is fundamental for an evaluation of the water residence times, dissolved pollutant and nutrient transport as well as oxygenation.

Diverging Cave- and River-Dwelling Newts Exert the Same Mate Preference in their Native Light Conditions

When colonizing a new habitat, populations must adapt their sexual behaviour to new ecological constraints. Because caves display drastically different conditions from surface habitats and cave animals are deprived from visual information, hypogean populations are expected to have modified their mate preference and signalling behaviour after cave colonization. Here, we experimentally examined the female preference and the sexual behaviour of brook newts Calotriton asper from different cave and river populations, either in light or in darkness. Our results suggest that females prefer large individuals in both hypogean and epigean populations, but that this preference is only expressed in the light conditions of their native habitat. Hence, some mate choice criteria would be maintained across genetically divergent populations and throughout dissimilar habitats. However, this sexual behaviour is likely to be expressed via a different sensory pathway in the different habitats, suggesting that a sensory shift has occurred in cave populations, enabling animals to communicate through a non-visual channel.

Organic matter processing and soil evolution in a braided river system

Traditionally, braided river research has considered flow, sediment transport processes and, recently, vegetation dynamics in relation to river morphodynamics. However, if considering the development of woody vegetated patches over a time scale of decades, we must consider the extent to which soil forming processes, particularly related to soil organic matter, impact the alluvial geomorphic-vegetation system. Here we quantify the soil organic matter processing (humification) that occurs on young alluvial landforms. We sampled different geomorphic units, ranging from the active river channel to established river terraces in a braided river system. For each geomorphic unit, soil pits were used to sample sediment/soil layers that were analysed in terms of grain size (<2mm) and organic matter quantity and quality (RockEval method). A principal components analysis was used to identify patterns in the dataset. Results suggest that during the succession from bare river gravels to a terrace soil, there is a transition from small amounts of external organic matter supply provided by sedimentation processes (e.g. organic matter transported in suspension and deposited on bars), to large amounts of autogenic in situ organic matter production due to plant colonisation. This appears to change the time scale and pathways of alluvial succession (bio-geomorphic succession). However, this process is complicated by: the ongoing possibility of local sedimentation, which can serve to isolate surface layers via aggradation from the exogenic supply; and erosion which tends to create fresh deposits upon which organic matter processing must re-start. The result is a complex pattern of organic matter states as well as a general lack of any clear chronosequence within the active river corridor. This state reflects the continual battle between deposition events that can isolate organic matter from the surface, erosion events that can destroy accumulating organic matter and the early ecosystem processes necessary to assist the co-evolution of soil and vegetation. A key question emerges over the extent to which the fresh organic matter deposited in the active zone is capable of significantly transforming the local geochemical environment sufficiently to accelerate soil development.

Estimating traveltimes and groundwater flow patterns using 3D time-lapse crosshole ERT imaging of electrical resistivity fluctuations induced by infiltrating river water

The infiltration of river water into aquifers is of high relevance to drinking-water production and is a key driver of biogeochemical processes in the hyporheic and riparian zone, but the distribution and quantification of the infiltrating water are difficult to determine using conventional hydrological methods (e.g., borehole logging and tracer tests). By time-lapse inverting crosshole ERT (electrical resistivity tomography) monitoring data, we imaged groundwater flow patterns driven by river water infiltrating a perialpine gravel aquifer in northeastern Switzerland. This was possible because the electrical resistivity of the infiltrating water changed during rainfall-runoff events. Our time-lapse resistivity models indicated rather complex flow patterns as a result of spatially heterogeneous bank filtration and aquifer heterogeneity. The upper part of the aquifer was most affected by the river infiltrate, and the highest groundwater velocities and possible preferential flow occurred at shallow to intermediate depths. Time series of the reconstructed resistivity models matched groundwater electrical resistivity data recorded on borehole loggers in the upper and middle parts of the aquifer, whereas the resistivity models displayed smaller variations and delayed responses with respect to the logging data. in the lower part. This study demonstrated that crosshole ERT monitoring of natural electrical resistivity variations of river infiltrate could be used to image and quantify 3D bank filtration and aquifer dynamics at a high spatial resolution.

Erosion and channel change as factors of landslides and valley formation in Champlain sea clays: The Chacoura river, Quebec, Canada

The Champlain Sea clays of Eastern Canada are incised by numerous rivers. Their slopes have been modified by landslides: on the Chacoura River near Trois-Rivières (Quebec), several large landslide scars, more or less recent, are visible. The role of erosion (channel incision, lateral channel migration and erosion of slopes due to agricultural drainage) as a trigger of these landslides is important. The aim of this study is to understand how erosion and landslides are related to valley development. From a detailed analysis of aerial photographs and DEMs, a map of the phenomena has been drawn by identifying various elements such as landslides, limits of the slope, position of the channel, and the area covered by forest. It is shown that channel change and erosion are strongly linked to landslides by the fact that they change the bank morphology in an unstable way. A slide in itself is a natural way for the slope to achieve stability. But when it occurs in a stream, it creates a disturbance to the stream flow enhancing local erosion which may change the river path and generate more erosion downstream or upstream resulting in more slides. Cross-valley sections and a longitudinal profile show that landslides are a major factor of valley formation. It appears that the upper part of the Chacoura River valley is still unaffected by landslides and has V-shaped sections. The lower part has been subject to intense erosion and many landslide scars can be seen. This shows that the valley morphology is transient, and that future activity is more likely to occur in the upper part of the river. Therefore the identification of areas prone to erosion will help determine the possible location of future large landslides just like the ones that occurred in the lower part.

The Dutung-Thaktote extensional fault zone and nappe structures documented by illite crystallinity and clay-mineral paragenesis in the Tethys Himalaya between Spiti river and Tso Morari, NW India

Numerous measurements by XRD of the Scherrer width at half-peak height (001 reflection of illite), coupled with analyses of clay-size assemblages, provide evidence for strong variations in the conditions of low temperature metamorphism in the Tethyan Himalaya metasediments between the Spiti river and the Tso Morari. Three sectors can be distinguished along the Spiti river-Tso Morari transect. In the SW, the Takling and Parang La area is characterised by a metamorphism around anchizone-epizone boundary conditions. Further north, in the Dutung area, the metamorphic grade abruptly decreases to weak diagenesis, with the presence of mixed-layered clay phases. At the end of the profile towards the NE, a progressive metamorphic increase up to greenschist facies is recorded, marked by the appearance of biotite and chloritoid. The combination of these data with the structural. observations permits to propose that a nappe stack has been crosscut by the younger Dutung-Thaktote extensional fault zone (DTFZ). The change in metamorphism across this zone helps to assess the displacements which occurred during synorogenic extension. In the SW and NE parts of the studied transect, a burial of 12 km has been estimated, assuming a geothermal gradient of 25 degrees C/km. In the SW part, this burial is due to the juxtaposition of the Shikar Beh and Mata nappes and in the NE part, solely to burial beneath the Mata nappe. In the central part of the profile, the effect of the DTFZ is to bring down diagenetic sediments in-between the two aforesaid metamorphic zones. The offset along the Dutung-Thaktote normal faults is estimated at 16 km.

3D crosshole ERT for aquifer characterization and monitoring of infiltrating river water

The hydrogeological properties and responses of a productive aquifer in northeastern Switzerland are investigated. For this purpose, 3D crosshole electrical resistivity tomography (ERT) is used to define the main lithological structures within the aquifer (through static inversion) and to monitor the water infiltration from an adjacent river. During precipitation events and subsequent river flooding, the river water resistivity increases. As a consequence, the electrical characteristics of the infiltrating water can be used as a natural tracer to delineate preferential flow paths and flow velocities. The focus is primarily on the experiment installation, data collection strategy, and the structural characterization of the site and a brief overview of the ERT monitoring results. The monitoring system comprises 18 boreholes each equipped with 10 electrodes straddling the entire thickness of the gravel aquifer. A multi-channel resistivity system programmed to cycle through various four-point electrode configurations of the 180 electrodes in a rolling sequence allows for the measurement of approximately 15,500 apparent resistivity values every 7 h on a continuous basis. The 3D static ERT inversion of data acquired under stable hydrological conditions provides a base model for future time-lapse inversion studies and the means to investigate the resolving capability of our acquisition scheme. In particular, it enables definition of the main lithological structures within the aquifer. The final ERT static model delineates a relatively high-resistivity, low-porosity, intermediate-depth layer throughout the investigated aquifer volume that is consistent with results from well logging and seismic and radar tomography models. The next step will be to define and implement an appropriate time-lapse ERT inversion scheme using the river water as a natural tracer. The main challenge will be to separate the superposed time-varying effects of water table height, temperature, and salinity variations associated with the infiltrating water.