Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Planar deformation features and impact glass in inclusions from the Vredefort Granophyre,South Africa

Abstract— The Vredefort Granophyre represents impact melt that was injected downward into fractures in the floor of the Vredefort impact structure, South Africa. This unit contains inclusions of country rock that were derived from different locations within the impact structure and are predominantly composed of quartzite, feldspathic quartzite, arkose, and granitic material with minor proportions of shale and epidiorite. Two of the least recrystallized inclusions contain quartz with single or multiple sets of planar deformation features. Quartz grains in other inclusions display a vermicular texture, which is reminiscent of checkerboard feldspar. Feldspars range from large, twinned crystals in some inclusions to fine‐grained aggregates that apparently are the product of decomposition of larger primary crystals. In rare inclusions, a mafic mineral, probably biotite or amphibole, has been transformed to very fine‐grained aggregates of secondary phases that include small euhedral crystals of Fe‐rich spinel. These data indicate that inclusions within the Vredefort Granophyre were exposed to shock pressures ranging from <5 to 8–30 GPa. Many of these inclusions contain small, rounded melt pockets composed of a groundmass of devitrified or metamorphosed glass containing microlites of a variety of minerals, including K‐feldspar, quartz, augite, low‐Ca pyroxene, and magnetite. The composition of this devitrified glass varies from inclusion to inclusion, but is generally consistent with a mixture of quartz and feldspar with minor proportions of mafic minerals. In the case of granitoid inclusions, melt pockets commonly occur at the boundaries between feldspar and quartz grains. In metasedimentary inclusions, some of these melt pockets contain remnants of partially melted feldspar grains. These melt pockets may have formed by eutectic melting caused by inclusion of these fragments in the hot (650 to 1610 °C) impact melt that crystallized to form the Vredefort Granophyre.     

Late Weichselian glacial maximum on Andøya, North Norway

Morphological, seismic and lithostratigraphic investigations of a moraine deposit at Bleik (the Bleik moraine), northern Andøya, show short-distance transported till overlying long-distance transported predominantly glaciofluvial ice-marginal deposits. Consolidated glaciomarine sediments from a core at present sea-level, c . 400 m distally to the moraine complex, contain 31 species of foraminifera, among which Cassidulina reniforme, Islandiella helenae and Trifarina fluens dominate, and fragments of the molluscs Mya truncata and Astarte sp. and the echinoid Strongylocentrotus sp. Amino acid analyses of the foraminifera Cibicides lobatulus and the mollusc Mya truncata indicate ages between 22,000 and 16,000 BP. Radiocarbon dating of fragments of Mya truncata from the upper part of the core gave an age of 17,940 ± 245 BP, while a dating of unidentified shell fragments from the lower part gave an infinite age (>40,000 BP). The sediment was probably disturbed by icebergs beyond the end moraine zone, and the radiocarbon and amino acid dating of Mya truncata therefore represent a maximum age for this process. This new evidence indicates two phases with a higher relative sea-level than at present at Bleik, c . 18,000 and >40,000 BP. The Bleik moraine probably represents the early Late Weichselian glacial maximum ( c . 22,000 BP), while the underlying deglaciation deposit and associated beach-ridge (Bruvollen) is of pre-Late Weichselian age. Moraine ridges 3–4 km to the south of Bleik probably indicate advances of local glaciers between 22,000 and 18,000 BP.     

Isolation of resonance in acoustic backscatter from elastic targetsusing adaptive estimation schemes

It has been shown that at certain frequencies the acoustic backscatter from elastic targets exhibits certain resonance behavior which closely relates to the physical properties of the target such as dimension, thickness, and composition. The purpose of this paper is to develop an automated approach for identifying the presence of resonance in the acoustic backscatter from an unknown underwater target by isolating the resonance part from the specular contribution. An adaptive transversal filter structure is used to estimate the specular part of the backscatter and consequently the error signal would provide an estimate of the resonance part. An important aspect of this scheme lies in the fact that it does not require an underlying model for the elastic return. The adaptation rule is based upon fast recursive least squares (RLS) learning. The approach taken in this paper is general in the sense that it can be applied to targets of unknown geometry and thickness and, further, does not require any a priori information about the target and/or the environment. Test results on acoustic data are presented which indicate the effectiveness of the proposed approach