Predictive potential of Perzyna viscoplastic modelling for granular geomaterials

This paper reappraises Perzyna-type viscoplasticity for the constitutive modelling of granular geomaterials, with emphasis on the simulation of rate/time effects of different magnitude. An existing elasto-plastic model for sands is first recast into a Perzyna viscoplastic formulation and then calibrated/validated against laboratory test results on Hostun sand from the literature. Notable model features include (1) enhanced definition of the viscous nucleus function and (2) void ratio dependence of stiffness and viscous parameters, to model the pycnotropic behaviour of granular materials with a single set of parameters, uniquely identified against standard creep and triaxial test results. The comparison between experimental data and numerical simulations points out the predicative capability of the developed model and the complexity of defining a unique viscous nucleus function to capture sand behaviour under different loading/initial/boundary and drainage conditions. It is concluded that the unified viscoplastic simulation of both drained and undrained response is particularly challenging within Perzyna's framework and opens to future research in the area. The discussion presented is relevant, for instance, to the simulation of multiphase strain localisation phenomena, such as those associated to slope stability problems in variably saturated soils.     

On the rheological characterisation of liquefied sands through the dam-breaking test

This paper concerns the rheological characterisation of liquefied sands as non-Newtonian Bingham fluids. For this purpose, dam-breaking laboratory tests are often executed and interpreted, offering a viable option to identify the properties of fluidised water-soil mixtures. However, limited attention has been devoted so far to clarify what variables and measurements would allow unambiguous calibration of Bingham parameters, namely, the viscosity η and the yield stress τ_y. The numerical results of parametric studies based on the particle finite element method (PFEM) are critically inspected to gain deeper insight into the problem. First, it is confirmed that multiple η − τ_y pairs may reproduce the same experimental evidence when formed by only one measurement—usually, the post–dam-breaking displacement of the bottom toe (tip) of the liquefied mass. Then, two alternative procedures are proposed for unambiguous identification of both η and τ_y: one is based on monitoring the evolving aspect ratio of the fluid mass during free, gravity-driven flow; the other relies on a slightly different dam-breaking test, also including impact against a rigid obstacle. In particular, the latter approach reduces the relevant duration of the test, reducing the possible influence of reconsolidation effects on the calibration of rheological parameters.     

Analysis of lateral loading of pile groups using embedded beam elements with interaction surface

The numerical simulation of soil-pile interaction problems, by means of full 3D finite element models, involves a large number of degrees of freedom (DOF) and difficulties during the mesh generation process. In order to reduce the unknowns and simplify and properly analyze such class of geotechnical problems, the so-called embedded beam elements (EBE) have recently been developed. In a preceding contribution of the authors, an improved EBE formulation, which brings into play the soil-pile interaction surface, was proposed with the aim to localize material plasticity in the soil surrounding the pile. This embedded beam model couples two different finite elements, each described by distinct kinematics (ie, solid and beam). The coupling is incorporated in the formulation by means of kinematical constrains established over the solid and beam displacement fields on the interaction surface. One of the main advantages of the embedded elements is that the addition of beams structural members immersed within the 3D soil model does not represent a constraint for the solid mesh, which can be adopted independently from the beam mesh. In this paper, the lateral loading of pile groups is studied by means of the proposed EBE approach with elasto-plastic interfaces. In order to represent a rigid cap, a master node and a special set of kinematical restrictions are incorporated into the formulation. The paper presents results obtained by means of the present formulation compared against other well-established analysis methods and test results published in the literature, for both elastic and elasto-plastic cases.     

Aysén seismic swarm (January 2007) in southern Chile: analysis using Joint Hypocenter Determination

A seismic swarm at the Aysén fjord started in January, 2007, reaching its highest activity between January and April, 2007. It was punctuated by the Mw5.3 23rd January event and the Mw6.2 21st April event which triggered a water-wave due to massive landslides; both with dextral strike slip focal mechanisms. More than 100 events were relocated with the Joint Hypocenter Determination method to study the depth distribution of the events of the swarm recorded by the Chile Ridge Subduction Project local network. The events are associated to a small area of approximately 7?×?7 km located in the Aysén fjord, on one of the main branches of the Liquiñe–Ofqui fault. Distribution in depth is located between 0 and 8 km and presents a west high dip with almost NS strike, consistent with the main trending of the Liquiñe–Ofqui fault and the published focal mechanisms. These characteristics, among others, suggest that the Liquiñe–Ofqui fault is active in this region and thus seismic risk has to be re-determined.     

Trends of tropical deforestation in Southeast Mexico

We assessed deforestation in Southeast Mexico (a surface area of 29 000 km² in seven states) through the comparison of land use/land cover maps at a scale of 1:250 000. This facilitated mapping of the land use/land cover change (LULCC) processes and calculation of the rates of change and the change matrix for the period 1978–2000. An original method was used to assess the accuracy of the LULCC map. The verification sites were selected through a stratified random sampling and were corroborated with aerial photographs for 1978 and 2000. Error matrices were elaborated using both hard and fuzzy set approaches in order to take into account the errors related to generalization of the map in fragmented landscapes. The results showed an average annual deforestation rate of 1.1 per cent which represents an average annual loss of 190 000 ha of forest, or an estimated total reduction of 4.2 million ha over 22 years. Furthermore, deforestation processes are concentrated in some areas such as Yucatan and Chiapas states, which registered major forest conversions to grassland and slash‐burning. The overall accuracy of the LULCC map, assessed with hard and fuzzy set approaches, was 72 per cent and 88 per cent respectively.     

Mapping land cover changes in Mexico, 1976–2000 and applications for guiding environmental management policy

Land use conversion typically implicates deforestation and fragmentation of primary land cover types, which invariably translates into impoverishment of both natural and cultural capital. Understanding where conversion is taking place crucially underpins sound environmental policy instruments to prevent these enormous social and economic costs. This paper examines 30 years of semi‐detailed (1:250 000) land cover mapping in Mexico. Pre‐existing analogue databases describing land cover patterns in the 1970s and 1990s were reviewed, corrected, reorganized and transformed into a digital format. Current land cover patterns were depicted by conducting updated reinterpretation on Landsat ETM+ imagery. Digital cartographic overlaying was performed and the results were used to construct a spatially explicit land use/land cover change (LULCC) database with an additional accuracy assessment procedure. The value of the results of this analysis is also seen in the light of their direct applications for identifying critical watershed trends, for guiding the allocation of financial funds for sound land use planning and for assessing the effectiveness of established protected areas. This effort highlights the importance of new and more effective geographical approaches to depict, understand and contribute to informed measures to mitigate ongoing negative trends in land cover and climatic changes.     

Multifractality in Seismicity Spatial Distributions: Significance and Possible Precursory Applications as Found for Two Cases in Different Tectonic Environments

We explore fractal properties of two observed seismicity distributions prior to the 2003 M _w 7.4 Colima, Mexico and 1992 M _w 7.3 Landers, USA earthquakes, together with several mathematical fractal distributions and two non-fractal ones, in order to estimate minimum reliable sample sizes, determine whether fractality for observed seismicity is essentially different from random uniform distributions, and explore the possibility of extracting premonitory information from fractal characteristics of seismicity before large earthquakes. Sample sizes above 800 events for whole catalogs appear to be sufficient to maintain ordered multifractality and to yield dimension estimates that vary smoothly and reliably. Fractal estimates appear to be best for whole catalogs that include aftershocks. The fractal characteristics of spatial distributions of seismicity are essentially different from those of the uniform random distribution, which is the null hypothesis of a non-fractal distribution with minimum information. The fractal dimensions and afractality measures of seismicity distributions change with time and show distinctive behaviors associated with foreshocks and main events, although these behaviors are different for each example. Results suggest the possibility of a priori identification of foreshocks to large earthquakes. A combination of fractal dimension and afractality measures over time may be helpful in large earthquake premonitory studies.     

Simulation and properties of a non-homogeneous spring-block earthquake model with asperities

The spring-block model proposed by Olami, Feder and Christensen (OFC) has several properties that are similar to those observed in real seismicity. In this paper we propose a modification of the original model in order to take into account that in a real fault there are several regions with different properties (non-homogeneity). We define regions in the network that is reminiscent of the real seismic fault, with different sizes and elastic parameter values. We obtain the Gutenberg-Richter law for the synthetic earthquake distributions of magnitude and the stair-shaped plots for the cumulative seismicity. Again, as in the OFC-homogeneous case, we obtain the stability for the cumulative seismicity stair-shaped graphs in the long-term situation; this means that the straight line slopes that are superior bounds of the staircases have a behavior akin to the homogeneous case. We show that with this non-homogeneous OFC model it is possible to include the asperity concept to describe high-stress zones in the fault.     

Influence of substrate tectonic heritage on the evolution of composite volcanoes: Predicting sites of flank eruption, lateral collapse, and erosion

This paper aims to aid understanding of the complicated interplay between construction and destruction of volcanoes, with an emphasis on the role of substrate tectonic heritage in controlling magma conduit geometry, lateral collapse, landslides, and preferential erosion pathways. The influence of basement structure on the development of six composite volcanoes located in different geodynamic/geological environments is described: Stromboli (Italy), in an island arc extensional tectonic setting, Ollagüe (Bolivia–Chile) in a cordilleran extensional setting, Kizimen (Russia) in a transtensional setting, Pinatubo (Philippines) in a transcurrent setting, Planchon (Chile) in a compressional cordilleran setting, and Mt. Etna (Italy) in a complex tectonic boundary setting. Analogue and numerical modelling results are used to enhance understanding of processes exemplified by these volcanic centres. We provide a comprehensive overview of this topic by considering a great deal of relevant, recently published studies and combine these with the presentation of new results, in order to contribute to the discussion on substrate tectonics and its control on volcano evolution. The results show that magma conduits in volcanic rift zones can be geometrically controlled by the regional tectonic stress field. Rift zones produce a lateral magma push that controls the direction of lateral collapse and can also trigger collapse. Once lateral collapse occurs, the resulting debuttressing produces a reorganization of the shallow-level magma migration pathways towards the collapse depression. Subsequent landslides and erosion tend to localize along rift zones. If a zone of weakness underlies a volcano, long-term creep can occur, deforming a large sector of the cone. This deformation can trigger landslides that propagate along the destabilized flank axis. In the absence of a rift zone, normal and transcurrent faults propagating from the substrate through the volcano can induce flank instability in directions respectively perpendicular and oblique to fault strike. This destabilization can evolve to lateral collapse with triggering mechanisms such as seismic activity or magmatic intrusion.