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.