Radon anomalies and volcanic eruptions

A well-documented case of ²²²Rn anomaly preceding the eruptive activity of Karymsky volcano (Kamchatka) was recently reported in the literature. Stimulated by this example, we have attempted to utilize the available data on radon emanation from rocks, its solubility, and its circulation in waters to discuss how a ²²²Rn anomaly can be produced by magma approaching the surface. It is shown that the most likely process of radon release is the flushing of gases through pore fluids. Heating of extensive fracture surfaces by high-temperature gases may also be important. In order to survive in detectable amounts after moving distances greater than a few meters, radon must be transported by fast-moving fluids, such as those rising toward the surface through a fracture or in the ascending limbs of fluids convecting in very porous rocks under high temperature gradients. The pattern observed at Karymsky volcano is interpreted on the basis of these relations.     

A numerical study of damping

Diagonal damping matrices were computed for three systems which have non-proportional damping matrices. These diagonal damping matrices were computed on three bases, as follows: 1. After normalizing the equations of motion by the modal matrix, the diagonal terms are retained ignoring the non-diagonal terms. 2. Diagonal damping matrix is established by the optimization algorithm which minimizes the mean square error of the frequency response. 3. Diagonal damping is determined from the normalized differential equation by matching the peaks of the coupled and uncoupled system. The frequency responses for the three cases of one of the three systems are presented together with a comparison of the energy dissipation.     

The evolution of the mantle source beneath Mt. Etna (Sicily,Italy): from the 600 ka tholeiites to the recent trachybasaltic magmas

ABSTRACT

The spatial/temporal proximity of Mt. Etna to the Hyblean Plateau and the Aeolian slab makes the discussion on the nature of its mantle source/s extremely controversial. In this study, a detailed geochemical overview of the entire Mt. Etna evolutionary sequence and a comparison with the magmatism of the Hyblean Plateau was proposed to: (i) simulate the composition of Mt. Etna tholeiitic to alkaline primitive magmas in equilibrium with a fertile mantle source; (ii) model the nature, composition and evolution of the mantle source from the tholeiitic stage (600 ka) to present magmatism. According to our simulations, two amphibole + phlogopite-bearing spinel lherzolite sources are able to explain the wide range of Etnean primary magmas. The enrichment in LILE, ⁸⁷Sr/⁸⁶Sr, Rb and H₂O of the magmas emitted after 1971 (but also discontinuously generated in both historic and prehistoric times) are caused by different melting proportions of amphibole and phlogopite in a modally and compositionally homogeneous mantle domain, with melting degrees analogous to those required to produce magmas erupted prior to 1971. The behaviour of the hydrous phases during melting could be ascribed to a variable H₂O/CO₂ activity in the mantle source, in turn related to the heat/fluxes supply from the asthenospheric upwelling beneath Mt. Etna. All these considerations, strengthened by numerical models, are then merged to review the complex Pliocene/Lower Pleistocene to present day’s geodynamic evolution of eastern Sicily.     

Comparison of Integrated Geodetic Data Models and Satellite Radar Interferograms to Infer Magma Storage During the 1991–1993 Mt. Etna Eruption

We compare the results obtained from the modelling of EDM, GPS, levelling and tilt data measured in the first part of the 1991–1993 eruption at Etna to the InSAR data acquired during the second part. The geodetic changes are very marked in the first half of the eruption and constrain a deflation source located at a few kilometers of depth ( 3 km b.s.l.), in agreement with other independent geophysical evidence. SAR data, available during the second part of the eruption, were analysed for different time intervals in the second part of the eruption. The interpretation of SAR interferograms reveals a large-scale but less marked deflation of the volcano that could be caused by a deeper source. This second source is in accord with a second deeper anomaly revealed by recent seismic investigations. The combination of geodetic data modelling and SAR images suggests a complex plumbing system composed at least of two possible storage regions located at different depths.     

Interbasin underflow between closed Altiplano basins in Chile

Interbasin ground water movement of 200 to 240 L/sec occurs as underflow beneath a mountainous surface water divide separating the topographically higher Salar de Michincha from the topographically lower Salar de Coposa internally drained basins in the Altiplano of northern Chile. Salt-encrusted flats (salars) and saline lakes occur on the lowest parts of the basin floors and comprise the principal evaporative discharge areas for the basins. Because a surface water divide separates the basins, surface water drainage boundaries do not coincide with ground water drainage boundaries. In the region, interbasin ground water movement is usually not recognized, but occurs for selected basins, and at places is an important component of ground water budgets. With increasing development of water for mining industry and potential exportation of ground water from the Altiplano for use at coastal cities, demonstration and quantification of interbasin movement is important for assessment of sustainable ground water development in a region of extreme aridity. Recognition and quantification of interbasin ground water underflow will assist in management of ground water resources in the arid Chilean Altiplano environment.     

On the application of SAR interferometry to geomorphological studies: estimation of landform attributes and mass movements

This paper presents two examples of application of Synthetic Aperture Radar (SAR) interferometry (InSAR) to typical geomorphological problems. The principles of InSAR are introduced, taking care to clarify the limits and the potential of this technique for geomorphological studies. The application of InSAR to the quantification of landform attributes such as the slope and to the estimation of landform variations is investigated. Two case studies are presented. A first case study focuses on the problem of measuring landform attributes by interferometric SAR data. The interferometric result is compared with the corresponding one obtained by a Digital Elevation Model (DEM). In the second case study, the use of InSAR for the estimation of landform variations caused by a landslide is detailed.     

Morphometric Analysis on the Size, Shape and Areal Distribution of Glacial Cirques in the Maritime Alps (Western French-Italian Alps)

The morphometry of 432 glacial cirques in the Maritime Alps (Western French‐Italian Alps), studied over several years of fieldwork, was analysed with the use of a geographical information system. Some of the parameters automatically evaluated from digital elevation models required an objective and relatively new definition. In particular, cirque length was measured along a line that, from the threshold midpoint, splits the cirque into two equivalent surfaces; cirque width was automatically drawn as the longest line inscribed in the cirque and perpendicular to the length line. Significant correlations were found among the different factors and parameters analysed. In particular, cirque shape analysis showed that cirques develop allometrically in the three dimensions, i.e. more in length and width than in altitudinal range. Nevertheless cirques of the Maritime Alps have a regular, almost circular shape (mean L/W value = 1.07). The correlations among length, width and area are all very high (r²= 0.8–0.9). In terms of size, cirques show a wide range in area from 0.06 to 5.2 km² with a mean value of 0.4 km². The largest cirques are found on SSW‐facing slopes and at high elevations. Small cirques can be found at all altitudes but all those at high elevation are part of compound cirques at the main head valleys. Most cirques (37%) are characterized by a northern aspect; NE and SW are also frequent directions.     

Comparing the adequacy of alternative ground motion intensity measures for the estimation of structural responses

Practical methods for the probability‐based seismic assessment of structures make use of estimates of demands produced by earthquakes of different intensities. The uncertainties associated with these estimates are highly dependent on the variable adopted as the intensity measure (IM, e.g., PGA, spectral acceleration, etc.). This generates the need to compare the efficiency of an originally adopted IMwith that of a new candidate. This implies comparing the dispersion of the demand measure (DM, e.g., maximum interstorey drift ratio, ductility demand, etc.) conditional to each of the two IMs. In order to obtain the demand estimates in a conventional way, a full set of dynamic response analyses should be performed for each IM under scrutiny, i.e., multiple records scaled at several fixed values of each IM. The procedure developed here serves to accelerate this comparison avoiding the effort required to evaluate the dynamic responses of the structure for all the ground motion time histories considered every time that a new IM is adopted. For this purpose, use is made of available results of analyses performed for a different (i.e., the original) IM. Two methods are proposed: the direct method involves performing a regression of the results obtained from the original analyses, taking the candidate IM as the independent variable. The indirect method involves rebuilding the probability density function of the DM given a defined value of the candidate IM by means of the total probability theorem, using the results of the original analyses and certain data relating the two IMs. The proposed methods have been tested by application to several SDOF systems with different periods and different cyclic‐response backbone curves. The conditions affecting their approximation are explored, and some criteria to improve them are identified. The procedure can also be used to determine the optimum value of a parameter to be used in a parameter‐based IM. Copyright © 2004 John Wiley & Sons, Ltd.