Heterogeneity influences on stream water–groundwater interactions in a gravel-dominated floodplain

ABSTRACT

Floodplains are composed of complex depositional patterns of ancient and recent stream sediments, and research is needed to address the manner in which coarse floodplain materials affect stream–groundwater exchange patterns. Efforts to understand the heterogeneity of aquifers have utilized numerous techniques typically focused on point-scale measurements; however, in highly heterogeneous settings, the ability to model heterogeneity is dependent on the data density and spatial distribution. The objective of this research was to investigate the correlation between broad-scale methodologies for detecting heterogeneity and the observed spatial variability in stream/groundwater interactions of gravel-dominated alluvial floodplains. More specifically, this study examined the correlation between electrical resistivity (ER) and alluvial groundwater patterns during a flood event at a site on Barren Fork Creek, in the Ozark ecoregion of Oklahoma, USA, where chert gravels were common both as streambed and as floodplain material. Water table elevations from groundwater monitoring wells for a flood event on 1–5 May 2009 were compared to ER maps at various elevations. Areas with high ER matched areas with lower water table slope at the same elevation. This research demonstrated that ER approaches were capable of indicating heterogeneity in surface water–groundwater interactions, and that these heterogeneities were present even in an aquifer matrix characterized as highly conductive. Portions of gravel-dominated floodplain vadose zones characterized by high hydraulic conductivity features can result in heterogeneous flow patterns when the vadose zone of alluvial floodplains activates during storm events.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR X. Chen     

Approximating constant‐Q seismic propagation in the time domain

In this study, we investigate the accuracy of approximating constant‐Q wave propagation by series of Zener or standard linear solid (SLS) mechanisms. Modelling in viscoacoustic and viscoelastic media is implemented in the time domain using the finite‐difference (FD) method. The accuracy of numerical solutions is evaluated by comparison with the analytical solution in homogeneous media. We found that the FD solutions using three SLS relaxation mechanisms as well as a single SLS mechanism, with properly chosen relaxation times, are quite accurate for both weak and strong attenuation. Although the RMS errors of FD simulations using a single relaxation mechanism increase with increasing offset, especially for strong attenuation (Q = 20), the results are still acceptable for practical applications. The synthetic data of the Marmousi‐II model further illustrate that the single SLS mechanism, to model constant Q, is efficient and sufficiently accurate. Moreover, it benefits from less computational costs in computer time and memory.     

Classification of Harmonic Functions in the Exterior of a Unit Sphere and Its Application to the Description of the External Magnetic Field of the Earth

Geomagnetism and Aeronomy - The magnetic fields of a subterranean source observed on the Earth’s surface or from the orbit of a satellite are studied. We analyze the space of all harmonic...     

Fast and direct nonparametric procedures in the L-moment homogeneity test

Regional frequency analysis is an important tool to properly estimate hydrological characteristics at ungauged or partially gauged sites in order to prevent hydrological disasters. The delineation of homogeneous groups of sites is an important first step in order to transfer information and obtain accurate quantile estimates at the target site. The Hosking–Wallis homogeneity test is usually used to test the homogeneity of the selected sites. Despite its usefulness and good power, it presents some drawbacks including the subjective choice of a parametric distribution for the data and a poorly justified rejection threshold. The present paper addresses these drawbacks by integrating nonparametric procedures in the L-moment homogeneity test. To assess the rejection threshold, three resampling methods (permutation, bootstrap and Pólya resampling) are considered. Results indicate that permutation and bootstrap methods perform better than the parametric Hosking–Wallis test in terms of power as well as in time and procedure simplicity. A real-world case study shows that the nonparametric tests agree with the HW test concerning the homogeneity of the volume and the bivariate case while they disagree for the peak case, but that the assumptions of the HW test are not well respected.     

Condensed hyperelements method of non-vertical consistent boundaries for wave propagation analysis in irregular media

The study of wave propagation in finite/infinite media has many applications in geotechnical and structural earthquake engineering and has been a focus of research for the past few decades. This paper presents an analysis of 2D anti- plane problems （Love waves） and 2D in-plane problems （Rayleigh waves） in the frequency domain in media consisting of a near-field irregular and a far-field regular part. The near field part may contain structures and its boundaries with the far-field can be of any shape. In this study, the irregular boundaries of the near-field are treated as consistent boundaries, extending the concept of Lysmer＇s vertical consistent boundaries. The presented technique is called the Condensed Hyperelements Method （CHM）. In this method, the irregular boundary is limited to a vertical boundary at each end that is a consistent boundary at the far-field side. Between the two ends, the medium is discretized with hyperelements. Using static condensation, the stiffness matrix of the far-field is derived for the nodes on the irregular boundary. Examples of the application of the CHM illustrate its excellent accuracy and efficiency.     

Estimate of fine sediment deposit dynamics over a gravel bar using photography analysis

Three different methods to analyse fine sediment deposits on a gravel bar using pictures are presented in this paper. A manual digitization and deposits zone delineation are performed as well as two different automated procedures. The three methods are applied on aerial pictures taken in 2006 by a drone from a height around 150 m above the study site. Two other sets of pictures taken in 2010 are also studied: the first set was obtained from the left side bank of the river at approximately 15m above the gravel bar whereas the second one was taken from a helicopter flying 600～m above the ground. These methods were used to estimate the surface of fine sediment deposits before and after flushing events. They yield similar results even if the first automated procedure is able to capture smaller patches of fine sediments. The total surface of fine sediment deposits seems to be similar before and after a flushing event, but the distribution appears quite different. Before a flushing event, a significant amount of fine sediment deposits are mixed with coarser sediments. After the flushing event, one can observe more large fine sediment deposits located on the downstream part of the secondary channel and at the channel margin. Most of the small fine sediment deposit patches were washed out. A short discussion is provided on the possible dynamics of fine sediment deposits over the gravel bar.     

The Aterno valley strong-motion array: seismic characterization and determination of subsoil model

The paper focuses on the strong motion array deployed in the upper Aterno River Valley, in the immediate outskirts north-west of the town of L’Aquila, which is part of the Italian Strong Motion Network operated by the Department of Civil Protection. The array is composed of six accelerometric stations located along a cross section of the valley. The importance of this array relies on the fact that a large amount of high-quality records were obtained during the 2009 L’Aquila seismic sequence, from both the mainshock and several aftershocks. These data are especially important to investigate site effects in sediment-filled valleys during moderate earthquakes in epicentral area because well-documented observational studies are very limited in the literature. However, the main drawback for the study of site effects in the Aterno valley is the lack of a detailed knowledge of the geometry of the valley, soil layering and dynamic properties of materials. The main motivation for this study stems from the need to provide a reliable subsoil model of the valley coupled with high-quality strong motion data. Based on the above, in the framework of S4 project, a major effort was undertaken to get a trustworthy cross section of the valley by an ad hoc investigation, comprising geological and geotechnical surveys as well as an extensive geophysical campaign, characterized by both active and passive measurements. These results were complemented by additional geological and geotechnical data available in the literature. By merging all the information acquired, a 2D subsoil model of the transversal section of the upper Aterno valley has been produced. The valley is characterised by an asymmetric shape with a shallower rock basement at the western edge of the valley that deepens at the valley centre. Moreover, based on the results of geophysical tests, representative V_s values were assigned to the different lithologic units forming the alluvial deposits filling the valley. Shear wave velocity is a fundamental parameter for ground response studies and it is also effective in classifying the accelerometric station from a seismic point of view. The 2D model may be therefore, considered a benchmark model for future studies of site effects. It will offer the possibility to examine site effects with a complex underlying geology and to validate the results of numerical simulations of site response analyses with the numerous observations from earthquake recordings, both for weak and strong ground motion conditions.     

Determining the location of buried plastic water pipes from measurements of ground surface vibration

‘Mapping the Underworld’ is a UK-based project, which aims to create a multi-sensor device that combines complementary technologies for remote buried utility service detection and location. One of the technologies to be incorporated in the device is low-frequency vibro-acoustics, and techniques for detecting buried infrastructure, in particular plastic water pipes, are being investigated. One of the proposed techniques involves excitation of the pipe at some known location with concurrent vibrational mapping of the ground surface in order to infer the location of the remainder of the pipe. In this paper, measurements made on a dedicated pipe rig are reported. Frequency response measurements relating vibrational velocity on the ground to the input excitation were acquired. Contour plots of the unwrapped phase revealed the location of the pipe to within 0.1-0.2 m. Magnitude contour plots revealed the excitation point and also the location of the pipe end. By examining the unwrapped phase gradients along a line above the pipe, it was possible to identify the wave-type within the pipe responsible for the ground surface vibration. Furthermore, changes in the ground surface phase speed computed using this method enabled the location of the end of the pipe to be confirmed.