Groundwater–surface‐water interactions in a braided river: a tracer‐based assessment

Natural tracers (alkalinity and silica) were used to infer groundwater–surface‐water exchanges in the main braided reach of the River Feshie, Cairngorms, Scotland. Stream‐water samples were collected upstream and downstream of the braided section at fortnightly intervals throughout the 2001–2002 hydrological year and subsequently at finer resolution over two rainfall events. The braided reach was found to exert a significant downstream buffering effect on the alkalinity of these waters, particularly at moderate flows (4–8 m³ s^?1/?Q_30–70). Extensive hydrochemical surveys were undertaken to characterize the different source waters feeding the braids. Shallow groundwater flow systems at the edge of the braided floodplain, recharged by effluent streams and hillslope drainage, appeared to be of particular significance. Deeper groundwater was identified closer to the main channel, upwelling through the hyporheic zone. Both sources contributed to the significant groundwater–surface‐water interactions that promote the buffering effect observed through the braided reach. Their impact was less significant at higher flows (>15 m³ s^?1/>Q₁₀) when acidic storm runoff from the peat‐covered catchment headwaters dominated, as well as under baseflow conditions (<4 m³ s^?1/<Q₇₀), when upstream alkalinity was already buffered owing to headwater groundwater sources assuming dominance. The significant temporally and spatially dynamic influence of these groundwater–surface‐water interactions was therefore seen to have important implications for both catchment functioning and instream ecology. Copyright © 2004 John Wiley & Sons, Ltd.     

Occurrence of “monalbite” in nature: a TEM study

Natural “monalbite” (MA) has been observed on a submicron-scale in jadeite along with high albite (HA) and low albite (LA) in jadeite quartzite in the Dabie ultrahigh-pressure (UHP) metamorphic terrane, China. Transmission electron microscope (TEM) observation of inclusions in jadeite crystals has revealed mineral phases with C2/m and C1? structures. The cell parameters of the C2/m crystal structure correspond to those of MA, whereas the cell parameters of C1? structure are identical to those of HA and LA. The existence of the C2/m structure shows that the jadeite underwent partial change to MA on a submicron-scale during retrograde metamorphism. Albite transformation twinning of HA suggests the transformation of MA into HA. MA occurrence provides powerful evidence of high-temperature metamorphism during exhumation, probably above 930 °C. The preservation of MA may be due to the armouring by the jadeite and implies a very rapid cooling (quenching) during retrogression.     

Least‐squares reverse‐time migration in a matrix‐based formulation

This paper describes least‐squares reverse‐time migration. The method provides the exact adjoint operator pair for solving the linear inverse problem, thereby enhancing the convergence of gradient‐based iterative linear inversion methods. In this formulation, modified source wavelets are used to correct the source signature imprint in the predicted data. Moreover, a roughness constraint is applied to stabilise the inversion and reduce high‐wavenumber artefacts. It is also shown that least‐squares migration implicitly applies a deconvolution imaging condition. Three numerical experiments illustrate that this method is able to produce seismic reflectivity images with higher resolution, more accurate amplitudes, and fewer artefacts than conventional reverse‐time migration. The methodology is currently feasible in 2‐D and can naturally be extended to 3‐D when computational resources become more powerful.     

Application of a multi‐temporal,LiDAR‐derived,digital terrain model in a landslide‐volume estimation

Sediments produced by landslides are crucial in the sediment yield of a catchment, debris flow forecasting, and related hazard assessment. On a regional scale, however, it is difficult and time consuming to measure the volumes of such sediment. This paper uses a LiDAR‐derived digital terrain model (DTM) taken in 2005 and 2010 (at 2 m resolution) to accurately obtain landslide‐induced sediment volumes that resulted from a single catastrophic typhoon event in a heavily forested mountainous area of Taiwan. The landslides induced by Typhoon Morakot are mapped by comparison of 25 cm resolution aerial photographs taken before and after the typhoon in an 83.6 km² study area. Each landslide volume is calculated by subtraction of the 2005 DTM from the 2010 DTM, and the scaling relationship between landslide area and its volume are further regressed. The relationship between volume and area are also determined for all the disturbed areas (V_L = 0.452A_L^1.242) and for the crown areas of the landslides (V_L = 2.510A_L^1.206). The uncertainty in estimated volume caused by use of the LiDAR DTMs is discussed, and the error in absolute volume estimation for landslides with an area >10⁵ m² is within 20%. The volume–area relationship obtained in this study is also validated in 11 small to medium‐sized catchments located outside the study area, and there is good agreement between the calculation from DTMs and the regression formula. By comparison of debris volumes estimated in this study with previous work, it is found that a wider volume variation exists that is directly proportional to the landslide area, especially under a higher scaling exponent. Copyright © 2013 John Wiley & Sons, Ltd.     

Damping in tall buildings—a mechanism and a predictor

Only when the data base is of sufficient quality can a theoretical model of any process be justified. Such a data base has recently been assembled in the field of damping in tall buildings. It is now possible to specify a damping value appropriate to a particular building vibrating at a specified amplitude which is well correlated with values which are measured in practice. The mechanism on which this predictor is based suggests that the predictor will be applicable to other types of structure as well. As yet the data base for these other structures does not exist, but some circumstantial evidence is encouraging.     

Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

Bedload tracing in a high‐sediment‐load mountain stream

This paper reports a radiofrequency identification (RFID) tracing experiment implemented in a high‐sediment‐load mountain stream typical of alpine gravel‐bed torrents. The study site is the Bouinenc Torrent, a tributary to the Bléone River in southeast France that drains a 38·9‐km² degraded catchment. In spring 2008, we deployed 451 tracers with b‐axis ranging from 23 to 520 mm. Tracers were seeded along eight cross‐sections located in the upstream part of the lowest 2·3 km of the stream. Three tracer inventories were implemented in July 2008, 2009 and 2010. Recovery rates calculated for mobile tracers declined from 78% in 2008 to 45% in 2009 and 25% in 2010. Observations of tracer displacement revealed very high sediment dispersion, with frontrunners having travelled more than 2 km only three months after their deployment. The declining recovery rate over time was interpreted as resulting from rapid dispersion rather than deep burial. We evaluated that 64% of the tracers deployed in the active channel were exported from the 2·3‐km study reach three years after the onset of the tracing experiment. Travel distances were characterized by right‐skewed and heavy‐tailed distributions, correctly fitted by a power‐law function. This supports the idea that in gravel‐bed rivers with abundant sediment supply relative to transport capacity, bedload transport can be viewed as a superdiffusive sediment dispersion process. It is also shown that tracers initially deployed in the low‐flow channel were characterized by a 15‐ to 30‐fold increase of mobility compared to tracers deployed in gravel bars. Copyright © 2011 John Wiley & Sons, Ltd.     

Pressure‐driven,shoreline currents in a perennially ice‐covered,pro‐glacial lake in Antarctica,identified from a LiCl tracer injected into a pro‐glacial stream

The distribution of streamwater within ice‐covered lakes influences sub‐ice currents, biological activity and shoreline morphology. Perennially ice‐covered lakes in the McMurdo Dry Valleys, Antarctica, provide an excellent natural laboratory to study hydrologic–limnologic interactions under ice cover. For a 2 h period on 17 December 2012, we injected a lithium chloride tracer into Andersen Creek, a pro‐glacial stream flowing into Lake Hoare. Over 4 h, we collected 182 water samples from five stream sites and 15 ice boreholes. Geochemical data showed that interflow travelled West of the stream mouth along the shoreline and did not flow towards the lake interior. The chemistry of water from Andersen Creek was similar to the chemistry of water below shoreline ice. Additional evidence for Westward flow included the morphology of channels on the ice surface, the orientation of ripple marks in lake sediments at the stream mouth and equivalent temperatures between Andersen Creek and water below shoreline ice. Streamwater deflected to the right of the mouth of the stream, in the opposite direction predicted by the Coriolis force. Deflection of interflow was probably caused by the diurnal addition of glacial runoff and stream discharge to the Eastern edge of the lake, which created a strong pressure gradient sloping to the West. This flow directed stream momentum away from the lake interior, minimizing the impact of stream momentum on sub‐ice currents. It also transported dissolved nutrients and suspended sediments to the shoreline region instead of the lake interior, potentially affecting biological productivity and bedform development. Copyright © 2014 John Wiley & Sons, Ltd.     

Nested‐Scale Nutrient Flux in a Mixed‐Land‐Use Urbanizing Watershed

To improve quantitative understanding of mixed‐land‐use impacts on nutrient yields, a nested‐scale experimental watershed study design (n = 5) was applied in a 303(d), clean water act impaired urbanizing watershed of the lower Missouri River Basin, USA. From 2010 to 2013, water samples (n = 858 sample days per site) were analysed for total inorganic nitrogen (TIN‐N), nitrite (NO₂–N) nitrate (NO₃–N), ammonia (NH₃–N), and total phosphorus (TP‐P). Annual, seasonal, and monthly flow‐weighted concentrations (FWCs) and nutrient yields were estimated. Mean nutrient concentrations were highest where agricultural land use comprised 58% of the drainage area (NH₃ = 0.111 mg/l; NO₂ = 0.045 mg/l; NO₃ = 0.684 mg/l, TIN = 0.840 mg/l; TP = 0.127 mg/l). Average TP‐P increased by 15% with 20% increased urban land use area. Highly variable annual precipitation was observed during the study with highest nutrient yields during 2010 (record setting wet year) and lowest nutrient yields during 2012 (extreme drought year). Annual TIN‐N and TP‐P yields exceeded 10.3 and 2.04 kg ha^?1 yr^?1 from the agricultural dominated headwaters. Mean annual NH₃–N, NO₂–N, NO₃–N, TIN‐N, and TP‐P yields were 0.742, 0.400, 4.24, 5.38, and 0.979 kg ha^?1 yr^?1, respectively near the watershed outlet. Precipitation accounted for the majority of the explained variance in nutrient yields (R² values from 0.68 to 0.85). Nutrient yields were also dependent on annual precipitation of the preceding year (R² values from 0.87 to 0.91) thus enforcing the great complexity of variable mixed‐land‐use mediated source‐sink nutrient yield relationships. Study results better inform land managers and best management practices designed to mitigate nutrient pollution issues in mixed‐land‐use freshwater ecosystems. Copyright © 2015 John Wiley & Sons, Ltd.     

Predicting the displacement of triple pendulum™ bearings in a full‐scale shaking experiment using a three‐dimensional element

The accuracy of a series spring model to predict the peak displacement and displacement history of Triple Pendulum? (TP) bearings in a strongly shaken, full‐scale building is evaluated in this paper. The series spring model was implemented as a self‐contained three‐dimensional TP bearing element in OpenSees and is now available for general use. The TP bearing element contains the option for constant friction or a generalized friction model that accounts for the effect of instantaneous velocity and compression load on the friction coefficient. Comparison between numerical simulation and experimental data of a five‐story steel moment frame building shows that the peak displacement of isolation system can generally be predicted with confidence using a constant friction coefficient model. The friction coefficient model accounting for the effect of axial load and velocity leads to minor improvement over the constant friction coefficient models in some cases. Copyright © 2013 John Wiley & Sons, Ltd.     

Implementation of a second‐order paraxial boundary condition for a water‐saturated layered half‐space in plane strain

A half‐space finite element and a transmitting boundary are developed for a water‐saturated layered half‐space using a paraxial boundary condition. The exact dynamic stiffness of a half‐space in plane strain is derived and a second‐order paraxial approximation of the stiffness is obtained. A half‐space finite element and a transmitting boundary are then formulated. The development is verified by comparison of the dynamic stiffness of impermeable and permeable rigid strip foundations with other published results. The advantage of using the paraxial boundary condition in comparison with the rigid boundary condition is examined. It is shown that the paraxial boundary condition offers significant gain and the resulting half‐space finite element and transmitting boundary can represent the effects of a water‐saturated layered half‐space with good accuracy and efficiency. In addition, the numerical method described herein maintains the strengths and advantages of the finite element method and can be easily applied to demanding problems of soil–structure interaction in a water‐saturated layered half‐space. Copyright © 2010 John Wiley & Sons, Ltd.     

Life of a bifurcation in a gravel‐bed braided river

Channel bifurcation is a key element in braided rivers, determining the water and sediment distribution and hence controlling the morphological evolution. Recent theoretical and experimental findings, as well as field observations, showed that bifurcations in gravel‐bed braided rivers are often asymmetrical and highly unstable. In this paper field data are presented on a bifurcation in the Tagliamento River, northeast Italy. The planform configuration of the bifurcation and its temporal evolution was monitored by an automatic digital camera during a series of seven floods with different magnitudes. This remote sensing technique allowed a high temporal resolution (pictures were acquired every hour) that was proved to be essential in a highly dynamic system as the one considered here. Digitized maps of the channels provided information on the location of the bifurcation, the width of the anabranches, the angle between them, along with the occurrence and migration of sediment bars. Data were acquired at two different water levels, giving the possibility to compare low and high flow conditions. The monitored bifurcation is largely unstable and shows sudden changes in the water distribution, mainly driven by the bar migrating in the upstream channel and entering the distributaries. A relationship between width asymmetry and flood magnitude was observed, confirming previous analyses. Moreover, recent theoretical findings were applied, in order to test the possibility to estimate general trends in bifurcation evolution. The analysis pointed out the relevance of a correct assessment of the characteristic temporal scales, as the bifurcation evolves on a timescale similar to that of bar migration and flood duration. Understanding the interactions between these processes is therefore crucial in order to increase the ability to model and predict the morphological evolution of a braided network. Copyright © 2012 John Wiley & Sons, Ltd.     

Thomsen's parameters from p‐wave measurements in a spherical sample‡

The aim of this paper is to understand the seismic anisotropy of the overburden shale in an oilfield in the North West Shelf of Western Australia. To this end, we first find the orientation of the symmetry axis of a spherical shale sample from measurements of ultrasonic P‐wave velocities in 132 directions at the reservoir pressure. After transforming the data to the symmetry axis coordinates, we find Thomsen's anisotropy parameters δ and ? using these measurements and measurements of the shear‐wave velocity along the symmetry axis from a well log. To find these anisotropy parameters, we use a very fast simulated re‐annealing algorithm with an objective function that contains only the measured ray velocities, their numerical derivatives and the unknown elasticity parameters. The results show strong elliptical anisotropy in the overburden shale. This approach produces smaller uncertainty of Thomsen parameter δ than more direct approaches.     

Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm

Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.     

Changes in the saltation flux following a step‐change in macro‐roughness

The effect of a step change in macro‐roughness on the saltation process under sediment supply limited conditions was examined in the atmospheric boundary layer. For an array of roughness elements of roughness density λ = 0.045 (λ = total element frontal area/total surface area of the array) the horizontal saltation flux was reduced by 90% (±7%) at a distance of ≈150 roughness element heights into the array. This matches the value predicted using an empirical design model and provides confidence that it can be effectively used to engineer roughness arrays to meet sand flux reduction targets. Measurements of the saltation flux characteristics in the vertical dimension, including: saltation layer decay (e‐folding) height and particle size, revealed that with increasing distance into the array, the rate of mass flux change with increasing height decreased notably, and (geometric) mean particle diameter decreased. The distribution of the saltation mass flux in the vertical remains exponential in form with increasing distance into the roughness array, and the e‐folding height increases as well as increasing at a greater rate as particle diameter diminishes. The increase in e‐folding height suggests the height of saltating particles is increasing along with their mean speed. This apparent increase in mean speed is likely due to the preferential removal, or sequestration, of the slower moving particles across the size spectrum, as they travel through the roughness array. Copyright © 2018 John Wiley & Sons, Ltd.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Impact–entrainment relationship in a saltating cloud

The problem of impact–entrainment relationship is one of the central issues in understanding saltation, a primary aeolian transport mode. By using particle dynamic analyser measurement technology the movement of saltating particles at the very near‐surface level (1 mm above the bed) was detected. The impacting and entrained particles in the same impact–entrainment process were identified and the speeds, angle with respect to the horizontal, and energy of the impacting and entrained sand cloud were analysed. It was revealed that both the speed and angle of impacting and entrained particles vary widely. The probability distribution of the speed of impacting and entrained particles in the saltating cloud is best described by a Weibull distribution function. The mean impact speed is generally greater than the mean lift‐off speed except for the 0·1–0·2 mm sand whose entrainment is significantly influenced by air drag. Both the impact and lift‐off angles range from 0° to 180°. The mean lift‐off angles range from 39° to 94° while the mean impact angles range from 40° to 78°, much greater than those previously reported. The greater mean lift‐off and especially the mean impact angles are attributed to mid‐air collisions at the very low height, which are difficult to detect by conventional high‐speed photography and are generally ignored in the existing theoretical simulation models. The proportion of backward‐impacting particles also evidences the mid‐air collisions. The impact energy is generally greater than the entrainment energy except for the 0·1–0·2 mm sand. There exists a reasonably good correlation of the mean speed, angle and energy between the impacting and entrained cloud in the impact–entrainment process. The results presented in this paper deserve to be considered in modelling saltation. Copyright © 2002 John Wiley & Sons, Ltd.     

Evaluation of a dynamic multi‐class sediment transport model in a catchment under soil‐conservation agriculture

Soil erosion models are essential tools for the successful implementation of effective and adapted soil conservation measures on agricultural land. Therefore, models are needed that predict sediment delivery and quality, give a good spatial representation of erosion and deposition and allow us to account for various soil conservation measures. Here, we evaluate how well a modified version of the spatially distributed multi‐class sediment transport model (MCST) simulates the effectiveness of control measures for different event sizes. We use 8 year runoff and sediment delivery data from two small agricultural watersheds (0·7 and 3·7 ha) under optimized soil conservation. The modified MCST model successfully simulates surface runoff and sediment delivery from both watersheds; one of which was dominated by sheet and the other was partly affected by rill erosion. Moreover, first results of modelling enrichment of clay in sediment delivery are promising, showing the potential of MCST to model sediment enrichment and nutrient transport. In general, our results and those of an earlier modelling exercise in the Belgian Loess Belt indicate the potential of the MCST model to evaluate soil erosion and deposition under different agricultural land uses. As the model explicitly takes into account the dominant effects of soil‐conservation agriculture, it should be successfully applicable for soil‐conservation planning/evaluation in other environments. Copyright © 2008 John Wiley & Sons, Ltd.