Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations

The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.     

A New Method for the Estimation of Avalanche Distance Exceeded Probabilities

A crucial point in any methodology for avalanche hazard assessment is the evaluation of avalanche distance exceeded probability, i.e., the annual probability that any assigned location along a given path is reached or exceeded by an avalanche. Typically this problem is faced by estimating the snow volume in the starting zone that is likely to accumulate an average every T years by statistical analysis of snowfall record, and then using this volume as input to an appropriately calibrated avalanche dynamics model to determine the runout distancesfor this design event. This methodology identifies the areas that canbe affected by an avalanche for the considered value of the return period (i.e. the average interval of time for a certain event to repeat itself), ¯T. However, it does not allow us to evaluate the actual avalanche encounter probability for any given point in the runout zone. In the present work this probability is computed by numerical integration of the expression P(x) = ∫₀ ^∞ P^*(V)f(V) dV, where f is the probabilitydensity function (PDF) of the avalanche release volume V, and P^* is the probability of the point x being reached or passed by an avalanche if the release volume is V; this latter probability is calculated by avalanche dynamics simulations. The procedure is implemented using a one-dimensional hydraulic-continuum avalanche dynamic model, calibrated on data from different Italian Alpine ranges, and is applied to a real world hazard mapping problem.     

The Holocene Secche di Lazzaro phreatomagmatic succession (Stromboli,Italy): evidence of pyroclastic density current origin deduced by facies analysis and AMS flow directions

The edifice of Stromboli volcano gravitationally collapsed several times during its volcanic history (>100 ka–present). The largest Holocene event occurred during the final stage of the Neostromboli activity (∼13–5 ka), and was accompanied by the emplacement of phreatomagmatic and lahar deposits, known as the Secche di Lazzaro succession. A stratigraphic and paleomagnetic study of the Secche di Lazzaro deposits allows the interpretation of the emplacement and the eruptive processes. We identify three main units within the succession that correspond to changing eruption conditions. The lower unit (UA) consists of accretionary lapilli-rich, thinly bedded, parallel- to cross-stratified ash deposits, interpreted to indicate the early stages of the eruption and emplacement of dilute pyroclastic density currents. Upward, the second unit (UB) of the deposit is more massive and the beds thicker, indicating an increase in the sedimentation rate from pyroclastic density currents. The upper unit (UC) caps the succession with thick, immediately post-eruptive lahars, which reworked ash deposited on the volcano’s slope. Flow directions obtained by Anisotropy of Magnetic Susceptibility (AMS) analysis of the basal bed of UA at the type locality suggest a provenance of pyroclastic currents from the sea. This is interpreted to be related to the initial base-surges associated with water–magma interaction that occurred immediately after the lateral collapse, which wrapped around the shoulder of the sector collapse scar. Upward in the stratigraphy (upper beds of UA and UB) paleoflow directions change and show a provenance from the summit vent, probably related to the multiple collapses of a vertical, pulsatory eruptive column.     

LIDAR Data Filtering and DTM Interpolation Within GRASS

LIDAR (Light Detection and Ranging) is one of the most recent technologies in surveying and mapping. LIDAR is based on the combination of three different data collection tools: a laser scanner mounted on an aircraft, a Global Positioning System (GPS) used in phase differential kinematic modality to provide the sensor position and an Inertial Navigation System (INS) to provide the orientation. The laser sends towards the ground an infrared signal, which is reflected back to the sensor. The time employed by the signal, given the aircraft position and attitude, allows computation of the earth point elevation. In standard conditions, taking into account the flight (speed 200–250 km/hour, altitude 500–2,000 m) and sensor characteristics (scan angle ± 10–20 degrees, emission rate 2,000–50,000 pulses per second), earth elevations are collected within a density of one point every 0.5–3 m. The technology allows us therefore to obtain very accurate (5–20 cm) and high resolution Digital Surface Models (DSM). For many applications, the Digital Terrain Model (DTM) is needed: we have to automatically detect and discard from the previous DSM all the features (buildings, trees, etc.) present on the terrain. This paper describes a procedure that has been implemented within GRASS to construct DTMs from LIDAR source data.     

IFKIS-Hydro: an early warning and information system for floods and debris flows

IFKIS-Hydro is an information and warning system for hydrological hazards in small- and medium-scale catchments. The system collects data such as weather forecasts, precipitation measurements, water level gauges, discharge simulations and local observations of event-specific phenomena. In addition, IFKIS-Hydro incorporates a web-based information platform, which serves as a central hub for the submission and overview of data. Special emphasis is given to local information. This is accomplished particularly by human observers. In medium-scale catchments, discharge forecast models have an increasing importance in providing valuable information. IFKIS-Hydro was developed in several test regions in Switzerland and the first results of its application are available now. The system is constantly extended to additional regions and may become the standard for warning systems in smaller catchments in Switzerland.     

Local And Global Diffusion Along Resonant Lines in Discrete Quasi-integrable Dynamical Systems

We detect and measure diffusion along resonances in a quasi-integrable symplectic map for different values of the perturbation parameter. As in a previously studied Hamiltonian case (Lega et al., 2003) results agree with the prediction of the Nekhoroshev theorem. Moreover, for values of the perturbation parameter slightly below the critical value of the transition between Nekhoroshev and Chirikov regime we have also found a diffusion of some orbits along macroscopic portions of the phase space. Such a diffusion follows in a spectacular way the peculiar structure of resonant lines.     

A novel methodological approach for land subsidence prediction through data assimilation techniques

Computational Geosciences - Anthropogenic land subsidence can be evaluated and predicted by numerical models, which are often built over deterministic analyses. However, uncertainties and...     

Macroseismic effects highlight site response in Rome and its geological signature

A detailed analysis of the earthquake effects on the urban area of Rome has been conducted for the L??Aquila sequence, which occurred in April 2009, by using an online macroseismic questionnaire. Intensity residuals calculated using the mainshock and four aftershocks are analyzed in light of a very accurate and original geological reconstruction of the subsoil of Rome based on a large amount of wells. The aim of this work is to highlight ground motion amplification areas and to find a correlation with the geological settings at a subregional scale, putting in evidence the extreme complexity of the phenomenon and the difficulty of making a simplified model. Correlations between amplification areas and both near-surface and deep geology were found. Moreover, the detailed scale of investigation has permitted us to find a correlation between seismic amplification in recent alluvial settings and subsiding zones, and between heard seismic sound and Tiber alluvial sediments.     

Adiabatic chaos in the spin–orbit problem

We provide evidences that the angular momentum of a symmetric rigid body in a spin–orbit resonance can perform large scale chaotic motions on time scales which increase polynomially with the inverse of the oblateness of the body. This kind of irregular precession appears as soon as the orbit of the center of mass is non-circular and the angular momentum of the body is far from the principal directions with minimum (maximum) moment of inertia. We also provide a quantitative explanation of these facts by using the theory of adiabatic invariants, and we provide numerical applications to the cases of the 1:1 and 1:2 spin–orbit resonances.