Ground deformations in urban areas can be the result of a combination of multiple factors and pose several hazards to infrastructures and human lives. In order to monitor these phenomena, ...Interferometric Synthetic Aperture Radar (InSAR) techniques are applied. The obtained signals record the overlapping of the phenomena, and their separation is a relevant issue. In this framework, we explored a new multi-method approach based on the combination of Principal Component Analysis (PCA), Independent Component Analysis (ICA) and Hierarchal Clustering (HC) on the standardized results to distinguish the main trends and seasonal signals embedded in the time series of ground displacements, to understand spatial-temporal patterns, to correlate ground deformation phenomena with geological and anthropogenic factors, and to recognize the specific footprints of different ground deformation phenomena. This method allows us to classify the ground deformations at the site scale in the metropolitan area of Naples, which is affected by uplift cycles, subsidence, cavity instabilities and sinkholes. At the local scale, the results allow a kinematic classification using the extracted components and considering the effect of the radius of influence generated by each cavity, as it is performed from a theoretical point of view when the draw angle is considered. According to the results, among the classified cavities, 2% were assigned to subsidence and 11% to uplift kinematics, while the remaining were found to be stable. Furthermore, our results show that the centering of the Spatial-PCA (S-PCA) is representative of the region’s main trend, whereas Temporal-PCA (T-PCA) gives information about the displacement rates identified by each component.
This study investigates the dynamics of a spherical projectile impact onto a granular bed via numerical simulations by discrete element method (DEM). The granular bed is modeled as an assembly of ...polydisperse spherical particles and the projectile is represented by a rigid sphere. The DEM model is used to investigate the cratering process, including the dynamics of the projectile and energy transformation and dissipation. The cratering process is illustrated by tracking the motion of the projectile and granular particles in the bed. The numerical results show that the dynamics of the projectile follows the generalized Poncelet law that the final penetration depth is a power-law function of the falling height. The numerical results can match well the experimental data reported in the literature, demonstrating the reliability of the DEM model in analyzing the impact of a spherical projectile on a granular bed. Further analyses illustrate that the impact process consists of three main stages, namely the impact, penetration and collapse, as characterized by the evolution of projective velocity, strong force chains and crater shape. The initial kinetic and potential energy of the projectile is dissipated mainly by inter-particle friction which governs the projectile dynamics. The stopping time of projectile decreases as the initial impact velocity increases. The final penetration depth scales as one-third the power of total falling height and is inversely proportional to the macroscopic granular friction coefficient.
A block impact model based on the elasto-viscoplastic macro element approach is developed for regular base prismatic blocks. This model upgrades a previously conceived model for spherical boulders ...introducing (i) a rotational degree of freedom; (ii) a moment-rotation relationship; (iii) a toppling mechanism. The model can improve rockfall simulations by considering block geometry and the exchange between translational and rotational energies.
Model parameters were calibrated by using laboratory tests on vertical impacts. Parametric analyses were carried out to investigate for both vertical and inclined impacts, the effects of block shape and orientation. The influence of these factors on the impact force, the maximum penetration depth as well as the exchange between translational and rotational energies is discussed. A comparison with the available results for small and large scale laboratory tests shows model capabilities and put in evidence the “nonlinear” relationship between maximum acceleration (or equivalently the maximum contact force) and impact translational velocity.
For vertical impacts the trend of the maximum penetration depth is a function of the block shape. Prismatic blocks can experience larger values of maximum penetration than spherical blocks characterized by coincident masses and kinetic energies. In case of bouncing of a prismatic block the increment of normal maximum displacement with respect to spherical blocks ranges from about 66% for triangular base prisms to 132% for hexagonal base blocks. In case of no bouncing, the increments range from about 82% for triangular blocks to −32% for hexagonal blocks.
Maximum normal forces also depend on block shape and orientation. In case of a vertex impact with no bouncing, triangular blocks show a decrement in the maximum force of about 43% with respect to the spherical block. The increment of initial block angular velocity generates a reduction in both maximum penetration depths and impact forces.
We present a new computationally efficient methodology to estimate the probability of rainfall-induced slope failure based on mechanical probabilistic slope stability analyses coupled with a ...hydrogeological model of the upslope area. The model accounts for: (1) uncertainty of geotechnical and hydrogeological parameters; (2) rainfall precipitation recorded over a period of time; and (3) the effect of upslope topography. The methodology provides two key outputs: (1) time-varying conditional probability of slope failure; and (2) an estimate of the absolute frequency of slope failure over any time period of interest. The methodology consists of the following steps: first, characterising the uncertainty of the slope geomaterial strength parameters; second, performing limit equilibrium method stability analyses for the realisations of the geomaterial strength parameters required to calculate the slope probability of failure by a Monte Carlo Simulation. The stability analyses are performed for various phreatic surface heights. These phreatic surfaces are then matched to a phreatic surface time series obtained from the 1D Hillslope-Storage Boussinesq model run for the upslope area to generate Factor of Safety (FoS) time series. A time-varying conditional probability of failure and an absolute frequency of slope failure can then be estimated from these FoS time series. We demonstrate this methodology on a road slope cutting in Nepal where geotechnical tests are not readily conducted. We believe this methodology improves the reliability of slope safety estimates where site investigation is not possible. Also, the methodology enables practitioners to avoid making unrealistic assumptions on the hydrological input. Finally, we find that the time-varying failure probability shows marked variations over time as a result of the monsoon wet–dry weather.
Highlights
Probabilistic slope stability analyses are coupled with a hydrogeological hillslope model to estimate the probability and frequency of rainfall-induced slope failure.
The model accounts for the uncertainty about rainfall using a time-dependent method, and for uncertainty relative to the geomaterial properties.
The model is tested on a road cut slope in Nepal (mountainous area subject to a monsoon season) finding that the cut slope will fail every other year.
Time-varying failure probability shows marked variations over time as a result of the monsoon wet–dry weather.
The findings indicate that it is important to use a time-dependent system to represent rainfall variability for slope failure probability analysis.
This paper discusses a strategy to identify failure conditions in geomaterials simulated by elastoplastic constitutive laws. The main objective is to express different forms of failure through the ...same formalism. For this purpose, we use a set of material instability indices combining the concepts of loss of controllability and critical hardening modulus with a simple, but versatile, elastoplastic model for soils and soft rocks. This choice has allowed us to (i) compute the instability indices in analytical form, (ii) capture the implications of non-normality and prior deposition/lithification history and (iii) inspect a broad range of failure modes (e.g., brittle and ductile failure, static liquefaction and compaction banding). It is shown that, although each mode of failure has its own specific features, they can all be encapsulated in a unified mathematical representation. To obtain these results, the instability moduli must reflect the static/kinematic constraints that generate the failure process at stake. Thus, the instability indices are expressed as functions of both the hardening modulus and additional terms of kinematic origin, with the latter terms reflecting a control-dependence of the plastic response. Such results describe a procedure for achieving a unified definition of failure in elastoplastic geomaterials, which is closely linked to the theory of controllability and encompasses the intuitive notions of ‘hardening’ and ‘softening’ as particular cases.
In this paper the onset of instabilities in elastoplastic materials is theoretically studied and a conceptual basis for understanding the physical implications of a loss of uniqueness and/or ...existence of the incremental response is provided. For this purpose, the concept of test controllability is reinterpreted and mixed stress–strain loading programmes are accounted for. A set of scalar indices, the moduli of instability, related with the inception of an unstable response is introduced and their dependency on the loading programme is explicitly illustrated. The paper shows that the use of these newly defined scalar measures provides support for an alternative definition for mechanical stability, which is closely related with the mathematical notions of existence and uniqueness of the predicted incremental response. In the final section, some mathematical properties of the moduli of instability are discussed, suggesting a novel reinterpretation of other well established theories and providing additional tools for the future application of the proposed framework.
The presence of trees along the slope and block fragmentation at impact strongly affect rockfall dynamics and hazard as a consequence. However, these phenomena are rarely simulated explicitly in ...rockfall studies. We performed rockfall simulations by using the 3D rockfall simulator Hy-Stone, modeling both the presence of trees and fragmentation through specific algorithms implemented in the code. By comparing these simulations with a more classical approach that attempts to account implicitly for such phenomena in the model parameters and by using a new probabilistic rockfall hazard analysis (PRHA) method, we were able to quantify the impact of these phenomena on the design of countermeasures and on hazard.
The dynamic fragmentation of jointed rock blocks during rockslide avalanches has been investigated by discrete element method simulations for a multiple arrangement of a rock block sliding over a ...simple slope geometry. The rock blocks are released along an inclined sliding plane and subsequently collide onto a flat horizontal plane at a sharp kink point. The contact force chains generated by the impact appear initially at the bottom frontal corner of the rock block and then propagate radially upward to the top rear part of the block. The jointed rock blocks exhibit evident contact force concentration and discontinuity of force wave propagation near the joint, associating with high energy dissipation of granular dynamics. The corresponding force wave propagation velocity can be less than 200 m/s, which is much smaller than that of an intact rock (1,316 m/s). The concentration of contact forces at the bottom leads to high rock fragmentation intensity and momentum boosts, facilitating the spreading of many fine fragments to the distal ends. However, the upper rock block exhibits very low rock fragmentation intensity but high energy dissipation due to intensive friction and damping, resulting in the deposition of large fragments near the slope toe. The size and shape of large fragments are closely related to the orientation and distribution of the block joints. The cumulative fragment size distribution can be well fitted by the Weibull's distribution function, with very gentle and steep curvatures at the fine and coarse size ranges, respectively. The numerical results of fragment size distribution can match well some experimental and field observations.
Key Points
The dynamic fragmentation characteristics of jointed rock blocks during rockslide‐avalanches have been analyzed by discrete element method
The force concentration at the bottom leads to high fragmentation intensity and momentum boosts, increasing the fragment runout distance
The final size, shape, and spatial distribution of large fragments are closely related to the initial joint orientation and distribution
Extremely Energetic Rockfalls De Blasio, Fabio Vittorio; Dattola, Giuseppe; Crosta, Giovanni Battista
Journal of geophysical research. Earth surface,
October 2018, 2018-10-00, 20181001, Letnik:
123, Številka:
10
Journal Article
Recenzirano
Odprti dostop
Extremely energetic rockfalls (EERs) are defined here as rockfalls for which a combination of both large volume and free fall height of hundreds of meters results in energy larger than about 80 GJ ...released in a short time. Examples include several events worldwide. In contrast to low energy rockfalls where block disintegration is limited, in EERs the impact after free fall causes immediate release of energy much like an explosion. The resulting air blast can snap trees hundreds of meters ahead of the fall area. Pulverized rock at high speed can abrade vegetation in a process of sandblasting, and particles suspended by the blast and the subsequent debris cloud may travel farther than the impact zone, blanketing vast areas. Using published accounts and new data, we introduce physically based models formulated on analogies with explosions and explosive fragmentation to describe EERs. Results indicate that a portion of the initial potential energy of the block is spent in rock disintegration at impact (typically 0.2%–18%), while other sources of energy loss (air drag, seismic, sound, and ground deformation) are negligible; consequently, more than 80% of the potential energy is converted to kinetic energy of the fragmented block (ballistic projection, shock wave, sand blast, and dust cloud). We also propose simple estimates for the flow of the dust cloud associated with an EER and its long settling time. The areal extent of the affected zone is estimated from the energy balance and an empirical power law relationship.
Key Points
EER (extremely energetic rockfalls) are characterized by explosive fragmentation. Information on EERs from different sources is provided
Physical models are suggested for the description of different sequences of EERs, including air blast and rock powder cloud propagation
Comparison with field data is provided, and a new model to compute fragmentation energy from particle size distribution is proposed
The process and dynamics of rock fragmentation during the collapse of rockfalls and rock avalanches is a poorly developed topic. The most severe fragmentation often leads to the formation of a rock ...dust that rises to form a cloud suspended in the air. The understanding of fragmentation processes is hampered by the environmental disturbances that alter the dust cloud deposit shortly after deposition. Here, we study the fragmentation of the October 2017 Pousset rockfall, detached from a NNE facing steep bedrock wall in the permafrost zone, that involved 8,300m
3
of metamorphic rock and fell about 800 m. The collapse generated large boulders which rolled downslope and a thick and large dust cloud. The source and deposit were investigated, and dust cloud material was sampled at different locations to reconstruct an exponential thickness distribution and perform grain size characterization. The fragmentation energy was estimated by integrating the spectrum of the grains assuming that the fragmentation energy is proportional to the generated area. The fragmentation energy was found to be about 0.4% of the initial potential energy. Most probable fragmentation points and block deposition areas were evaluated and positioned by means of the HyStone 3D rockfall simulator. Furthermore, we calculated the flow rate of the suspended powder generated by the fragmentation process and compared the results with observations available for the evolution of the phenomenon and the collected samples. The Pousset event, in its relatively simple dynamics, may be a good testing ground to address the current theories of rockfall and rock avalanche fragmentation and dust cloud behavior.