Transport of moderately sorted gravel at low bed shear stresses: The role of fine sediment infiltration

A reliable estimation of sediment transport in gravel‐bed streams is important for various practical engineering and biological studies (e.g., channel stability design, bed degradation/aggradation, restoration of spawning habitat). In the present work, we report original laboratory experiments investigating the transport of gravel particles at low bed shear stresses. The laboratory tests were conducted under unsteady flow conditions inducing low bed shear stresses, with detailed monitoring of the bed topography using a laser scanner. Effects of bed surface arrangements were documented by testing loose and packed bed configurations. Effects of fine sediments were examined by testing beds with sand, artificial fine sand or cohesive silt infiltrated in the gravel matrix. Analysis of the experimental data revealed that the transport of gravel particles depends upon the bed arrangement, the bed material properties (e.g., size and shape, consolidation index, permeability) and the concentration of fine sediments within the surface layer of moving grains. This concentration is directly related to the distribution of fine particles within the gravel matrix (i.e., bottom‐up infiltration or bridging) and their transport mode (i.e., bedload or suspended load). Compared to loose beds, the mobility of gravel is reduced for packed beds and for beds clogged from the bottom up with cohesive fine sediments; in both cases, the bed shear stress for gravel entrainment increases by about 12%. On the other hand, the mobility of gravel increases significantly (bed shear stress for particle motion decreasing up to 40%) for beds clogged at the surface by non‐cohesive sand particles. Copyright © 2017 John Wiley & Sons, Ltd.     

Depth of closure: New calculation method based on sediment data

Obtaining depth of closure (DoC) in an accurate manner is a fundamental issue for coastal engineering, since good results for coastal structures and beach nourishment depend mainly on DoC. Currently, there are two methods for obtaining the DoC, mathematical formulations and profile surveys. However, these methods can incur important errors if one does not take into account the characteristics and morphology of the area, or if one does not have a sufficiently long time series. In this work the DoC is obtained from the break in the trend of the sediment with the depth, that is, in general with the increase of the depth a decrease in the size of the sediment takes place. However, at one point this tendency changes and the size increases, and then decreases again. When comparing the point where the minimum sediment size occurs before the increase, it is observed that the error incurred is small compared to other methods. If the Standard Deviation of Depth Change (SDDC) method is considered as the most accurate method, the error incurred by the proposed method is less than 7%. In addition, it can be seen that the dispersion of the sediment method always occurs outside the zone of bar movement. Whereas in the methods of profiles survey (using 2 cm precision profiles), sometimes the DoC is obtained within the active zone of bar movement. In addition, where the relative minimum of the median sediment size is found, and the sizes of 0.063 and 0.125 mm predominate in the composition of the sample. Therefore, this new method allows the precise location of the DoC to be obtained in a fast and simple way. Furthermore, this method has the advantage that it is not affected by the modifications that may be experienced by both the study area and the cross-shore beach profile.     

Simulating rainfall time-series: how to account for statistical variability at multiple scales?

Daily rainfall is a complex signal exhibiting alternation of dry and wet states, seasonal fluctuations and an irregular behavior at multiple scales that cannot be preserved by stationary stochastic simulation models. In this paper, we try to investigate some of the strategies devoted to preserve these features by comparing two recent algorithms for stochastic rainfall simulation: the first one is the modified Markov model, belonging to the family of Markov-chain based techniques, which introduces non-stationarity in the chain parameters to preserve the long-term behavior of rainfall. The second technique is direct sampling, based on multiple-point statistics, which aims at simulating a complex statistical structure by reproducing the same data patterns found in a training data set. The two techniques are compared by first simulating a synthetic daily rainfall time-series showing a highly irregular alternation of two regimes and then a real rainfall data set. This comparison allows analyzing the efficiency of different elements characterizing the two techniques, such as the application of a variable time dependence, the adaptive kernel smoothing or the use of low-frequency rainfall covariates. The results suggest, under different data availability scenarios, which of these elements are more appropriate to represent the rainfall amount probability distribution at different scales, the annual seasonality, the dry-wet temporal pattern, and the persistence of the rainfall events.     

On the significance of the climate-dataset time resolution in characterising wind-driven rain and simultaneous wind pressure. Part I: scalar approach

Stochastic Environmental Research and Risk Assessment - The joint action of wind-driven rain and wind pressure is the main cause of water penetration in building facades, which causes various...     

Equilibrium modeling of the beach profile on a macrotidal embayed low tide terrace beach

Eleven-year long time series of monthly beach profile surveys and hourly incident wave conditions are analyzed for a macrotidal Low Tide Terrace beach. The lower intertidal zone of the beach has a pluriannual cycle, whereas the upper beach profile has a predominantly seasonal cycle. An equilibrium model is applied to study the variation of the contour elevation positions in the intertidal zone as a function of the wave energy, wave power, and water level. When forcing the model with wave energy, the predictive ability of the equilibrium model is around 60% in the upper intertidal zone but decreases to 40% in the lower intertidal zone. Using wave power increases the predictive ability up to 70% in both the upper and lower intertidal zones. However, changes around the inflection point are not well predicted. The equilibrium model is then extended to take into account the effects of the tide level. The initial results do not show an increase in the predictive capacity of the model, but do allow the model free parameters to represent more accurately the values expected in a macrotidal environment. This allows comparing the empirical model calibration in different tidal environment. The interpretation of the model free parameter variation across the intertidal zone highlights the behavior of the different zones along the intertidal beach profile. This contributes to a global interpretation of the four model parameters for beaches with different tidal ranges, and therefore to a global model applicable at a wide variety sites.     

Formation events of shoreline sand waves on a gravel beach

Kilometric-scale shoreline sand waves (KSSW) have been observed in the north-east flank of the Dungeness Cuspate Foreland (southeastern coast of the UK). They consist of two bumps separated by embayments with a 350–450-m spacing. We have analysed 36 shoreline surveys of 2-km length using the Discrete Fourier Transformation (DFT), from 2005 to 2016, and seven topographic surveys encompassing the intertidal zone, from 2010 to 2016. The data set shows two clear formation events. In order to test the role of high-angle waves on the KSSW formation, the 10-year wave series is propagated from the wave buoy located at 43 m depth up to a location in front of the undulations at 4 m depth using the SWAN wave model. The dominating SW waves arrive with a very high incidence angle (～ 80°) while the NE waves arrive almost shore normal. The ratio R, which measures the degree of dominance of high-angle waves with respect to low-angle waves, correlates well with the shoreline DFT magnitude values of the observed wavelength undulations. In particular, the highest R values coincide with the formation events. Finally, a linear stability model based on the one-line approximation is applied to the Dungeness profile and the 10-year propagated wave series. It predicts accurately the formation moments, with positive growth rates in the correct order of magnitude for wavelengths similar to the observed ones. All these results confirm that the shoreline undulations in Dungeness are self-organized and that the underlying formation mechanism is the high-angle wave instability. The two detected formation events provide a unique opportunity to validate the existing morphodynamic models that include such instability.     

A first look at the effects of ionospheric signal bending on a globally processed GPS network

This study provides a first attempt at quantifying potential signal bending effects on the GPS reference frame, coordinates and zenith tropospheric delays (ZTDs). To do this, we homogeneously reanalysed data from a global network of GPS sites spanning 14 years (1995.0–2009.0). Satellite, Earth orientation, tropospheric and ground station coordinate parameters were all estimated. We tested the effect of geometric bending and dTEC bending corrections, which were modelled at the observation level based, in part, on parameters from the International Reference Ionosphere 2007 model. Combined, the two bending corrections appear to have a minimal effect on site coordinates and ZTDs except for low latitude sites. Considering five days (DOY 301–305, 28 October–1 November 2001) near ionospheric maximum in detail, they affect mean ZTDs by up to ~1.7 mm at low latitudes, reducing to negligible levels at high latitudes. Examining the effect on coordinates in terms of power-spectra revealed the difference to be almost entirely white noise, with noise amplitude ranging from 0.3 mm (high latitudes) to 2.4 mm (low latitudes). The limited effect on station coordinates is probably due to the similarity in the elevation dependence of the bending term with that of tropospheric mapping functions. The smoothed z-translation from the GPS reference frame to ITRF2005 changes by less than 2 mm, though the effect combines positively with that from the second order ionospheric refractive index term. We conclude that, at the present time, and for most practical purposes, the geometric and dTEC bending corrections are probably negligible at current GPS/reference frame precisions.     

Analysis of soil drying incorporating a constitutive model for curling

This paper focuses on the shrinkage behavior of soil specimens involving sand, kaolinite, and kaolinite/sand mixtures subjected to desiccation under controlled conditions. Both, free and restrained shrinkage conditions are studied. The experiments show that pure soils do not curl upon unrestrained shrinkage; however, (under the same conditions) kaolinite/sand mixtures exhibited a marked curling. Furthermore, the mixture with the higher sand content broke through the middle of the sample after displaying a significant curling. Soils subjected to restricted shrinkage developed cracks with slight curling. To simulate the observed behavior, a mechanical model able to reproduce the detachment of the soil sample from the mold is proposed in this work and implemented in a fully coupled hydro-mechanical finite-element code. It is concluded that suction and differential shrinkage are key factors influencing the curling behavior of soils. The proposed framework was able to satisfactorily explain and reproduce the different stages and features of soil behavior observed in the experiments.

     

The Radial Basis Functions Method for Improved Numerical Approximations of Geological Processes in Heterogeneous Systems

Mathematical Geosciences - A robust, high order modeling approach is introduced, based on the finite difference-based radial basis functions method, for solving the groundwater flow equation in the...     

Subaquatic slope instabilities: The aftermath of river correction and artificial dumps in Lake Biel (Switzerland)

River engineering projects are developing rapidly across the globe, drastically modifying water courses and sediment transfer. Investigation of the impact of engineering works focuses usually on short-term impacts, thus a longer-term perspective is still missing on the effects that such projects have. The ‘Jura Water Corrections’ – the largest river engineering project ever undertaken in Switzerland – radically modified the hydrological system of Lake Biel in the 19th and 20th Century. The deviation of the Aare River into Lake Biel more than 140 years ago, in 1878, thus represents an ideal case study to investigate the long-term sedimentological impacts of such large-scale river rerouting. Sediment cores, along with new high-resolution bathymetric and seismic reflection datasets were acquired in Lake Biel to document the consequences of the Jura Water Corrections on the sedimentation history of Lake Biel. Numerous subaquatic mass transport structures were detected on all of the slopes of the lake. Notably, a relatively large mass transport complex (0·86 km²) was observed on the eastern shore, along the path of the Aare River intrusion. The large amount of sediment delivered by the Aare River since its deviation into the lake likely caused sediment overloading resulting in subaquatic mass transport. Alternatively, the dumping since 1963 in a subaquatic landfill of material excavated during the second phase of river engineering, when the channels flowing into and out of Lake Biel were widened and deepened, might have triggered the largest mass transport, dated to 1964 or 1965. Additional potential triggers include two nearby small earthquakes in 1964 and 1965 (M_W 3·9 and 3·2, respectively). The data for this study indicate that relatively large mass transports have become recurrent in Lake Biel following the deviation of the Aare River, thus modifying hazard frequency for the neighbouring communities and infrastructure.