Splash-saltation trajectories of soil particles under wind-driven rain

It is usually recognized that relatively large amounts of soil particles cannot be transported by raindrop splashes under windless rain. However, the splash-saltation process can cause net transportation in the prevailing wind direction since variations in splash-saltation trajectory due to the wind are expected in wind-driven rain. Therefore, determining the combined effect of rain and wind on the process should enable improvement of the estimation of erosion for any given prediction technique. This paper presents experimental data on the effects of slope aspect, slope gradient, and horizontal wind velocity on the splash-saltation trajectories of soil particles under wind-driven rain. In a wind tunnel facility equipped with a rainfall simulator, the rains driven by horizontal wind velocities of 6, 10, and 14 m s⁻¹ were allowed to impact three agricultural soils packed into 20×55 cm soil pans placed at both windward and leeward slopes of 7%, 15%, and 20%. Splash-saltation trajectories were measured by trapping the splashed particles at distances downwind on a 7-m uniform slope segment in the upslope and downslope directions, respectively, for windward and leeward slopes. Exponential decay curves were fitted for the mass distribution of splash-saltation sediment as a function of travel distance, and the average splash-saltation trajectory was derived from the average value of the fitted functions. The results demonstrated that the average trajectory of a raindrop-induced and wind-driven soil particle was substantially affected by the wind shear velocity, and it had the greatest correlation (r=0.96 for all data) with the shear velocity; however, neither slope aspect nor slope gradient significantly predicted the splash-saltation trajectory. More significantly, a statistical analysis conducted with nonlinear regression model of C₁(u_*²/g) showed that average trajectory of splash saltation was approximately three times greater than that of typical saltating sand grain.     

The reflection of plane waves in a poroelastic half-space saturated with inviscid fluid

This paper discusses surface displacements, surface strain, rocking, and energy partitioning during reflection-of-plane waves in a fluid-saturated poroelastic half-space. The medium is modeled by Biot's theory, and is assumed to be saturated with inviscid fluid. A linear porosity-modulus relation based on experimental data on sandstones is used to determine the material parameters for Biot's model. Numerical results in terms of angle of incident waves and Poisson's ratio are illustrated for various porosities and degrees of solid frame stiffness. The results show that the amount of solid frame stiffness controls the response of a fluid-saturated porous system. A poroelastic medium with essentially dry-frame stiffness behaves like an elastic medium, and the influence of pore fluid increases as dry-frame stiffness is reduced. The effects of a second P-wave become noticeable in poroelastic media with low dry-frame stiffness.     

A 3D BEM-FEM methodology for simulation of high speed train induced vibrations

This work presents an efficient methodology for the analysis of vibrations in a railroad track system, induced by the passage of conventional and high-speed trains. To this end the Boundary Element Method is used to model the soil-tie system within the framework of impulse response techniques. Conventional Finite Element Methods along with Newmark's integration is used for the modeling of the rail system. The two methods are coupled at the tie-rail interface and the solution is obtained following a staggered, time marching scheme in an efficient manner. The methodology accounts for Soil-Structure Interaction and traveling wave effects. In addition, this work identifies the parameters that affect the efficient modeling of railroad track systems as they pertain to the proposed solution methodology.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Short‐term intertidal bar mobility on a ridge‐and‐runnel beach,Merlimont, northern France

Digital elevation models and topographic pro?les of a beach with intertidal bar and trough (ridge‐and‐runnel) morphology in Merlimont, northern France, were analysed in order to assess patterns of cross‐shore and longshore intertidal bar mobility. The beach exhibited a pronounced dual bar–trough system that showed cross‐shore stationarity. The bars and troughs were, however, characterized by signi?cant longshore advection of sand under the in?uence of suspension by waves and transport by strong tide‐ and wind‐driven longshore currents. Pro?le changes were due in part to the longshore migration of medium‐sized bedforms. The potential for cross‐shore bar migration appears to be mitigated by the large size of the two bars relative to incident wave energy, which is modulated by high vertical tidal excursion rates on this beach due to the large tidal range (mean spring tidal range = 8·3 m). Cross‐shore bar migration is also probably hindered by the well‐entrenched troughs which are maintained by channelled high‐energy intertidal ?ows generated by swash bores and by tidal discharge and drainage. The longshore migration of intertidal bars affecting Merlimont beach is embedded in a regional coastal sand transport pathway involving tidal and wind‐forced northward residual ?ows affecting the rectilinear northern French coast in the eastern English Channel. Copyright © 2004 John Wiley & Sons, Ltd.     

Knickpoint recession rate and catchment area: the case of uplifted rivers in Eastern Scotland

Knickpoint behaviour is a key to understanding both the landscape responses to a base‐level fall and the corresponding sediment fluxes from rejuvenated catchments, and must be accommodated in numerical models of large‐scale landscape evolution. Knickpoint recession in streams draining to glacio‐isostatically uplifted shorelines in eastern Scotland is used to assess whether knickpoint recession is a function of discharge (here represented by its surrogate, catchment area). Knickpoints are identified using DS plots (log slope versus log downstream distance). A statistically significant power relationship is found between distance of headward recession and catchment area. Such knickpoint recession data may be used to determine the values of m and n in the stream power law, E = KA^mSⁿ. The data have too many uncertainties, however, to judge definitively whether they are consistent with m = n = 1 (bedrock erosion is proportional to stream power and KPs should be maintained and propagate headwards) or m = 0·3, n = 0·7 (bedrock incision is proportional to shear stress and KPs do not propagate but degrade in place by rotation or replacement). Nonetheless, the E Scotland m and n values point to the dominance of catchment area (discharge) in determining knickpoint retreat rates and are therefore more consistent with the stream power law formulation in which bedrock erosion is proportional to stream power. Copyright © 2005 John Wiley & Sons, Ltd.     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.