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.     

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.     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

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.     

POLARIX: a pathfinder mission of X-ray polarimetry

Since the birth of X-ray astronomy, spectral, spatial and timing observation improved dramatically, procuring a wealth of information on the majority of the classes of the celestial sources. Polarimetry, instead, remained basically unprobed. X-ray polarimetry promises to provide additional information procuring two new observable quantities, the degree and the angle of polarization. Polarization from celestial X-ray sources may derive from emission mechanisms themselves such as cyclotron, synchrotron and non-thermal bremsstrahlung, from scattering in aspheric accreting plasmas, such as disks, blobs and columns and from the presence of extreme magnetic field by means of vacuum polarization and birefringence. Matter in strong gravity fields and Quantum Gravity effects can be studied by X-ray polarimetry, too. POLARIX is a mission dedicated to X-ray polarimetry. It exploits the polarimetric response of a Gas Pixel Detector, combined with position sensitivity, that, at the focus of a telescope, results in a huge increase of sensitivity. The heart of the detector is an Application-Specific Integrated Circuit (ASIC) chip with 105,600 pixels each one containing a full complete electronic chain to image the track produced by the photoelectron. Three Gas Pixel Detectors are coupled with three X-ray optics which are the heritage of JET-X mission. A filter wheel hosting calibration sources unpolarized and polarized is dedicated to each detector for periodic on-ground and in-flight calibration. POLARIX will measure time resolved X-ray polarization with an angular resolution of about 20 arcsec in a field of view of 15 × 15 arcmin and with an energy resolution of 20% at 6 keV. The Minimum Detectable Polarization is 12% for a source having a flux of 1 mCrab and 10⁵ s of observing time. The satellite will be placed in an equatorial orbit of 505 km of altitude by a Vega launcher. The telemetry down-link station will be Malindi. The pointing of POLARIX satellite will be gyroless and it will perform a double pointing during the earth occultation of one source, so maximizing the scientific return. POLARIX data are for 75% open to the community while 25% + SVP (Science Verification Phase, 1 month of operation) is dedicated to a core program activity open to the contribution of associated scientists. The planned duration of the mission is one year plus three months of commissioning and SVP, suitable to perform most of the basic science within the reach of this instrument. A nice to have idea is to use the same existing mandrels to build two additional telescopes of iridium with carbon coating plus two more detectors. The effective area in this case would be almost doubled.     

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.     

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...