Post-caldera activity in the Alban Hills volcanic district (Italy): 40Ar/39Ar geochronology and insights into magma evolution

We present 24 ⁴⁰Ar/³⁹Ar ages for the youngest volcanic products from the Alban Hills volcanic district (Rome). Combined with petrological data on these products, we have attempted to define the chronology of the most recent phase of activity and to investigate the magma evolution of this volcanic district. The early, mainly explosive activity of the Alban Hills spanned the interval from 561±1 to 351±3 ka. After approximately 50-kyr of dormancy, a mainly effusive phase of activity took place, accompanied by the strombolian activity of a small central edifice (Monte delle Faete). This second phase of activity spanned the interval 308±2 to 250±1 ka. After another dormancy period of approximately 50-kyr, a new, hydromagmatic phase of activity started at 200 ka at several centers located to the southwest of the Monte delle Faete edifice. After an initial recurrence period of approximately 50-kyr, which also characterized this new phase of activity, the longest dormancy period (approximately 80-kyr) in the history of the volcanic district preceded the start of the activity of the Albano and Giuturna centers at 70±1 ka. Results of our study suggest a quasi-continuous magmatic activity feeding hydromagmatic centers with a new acme of volcanism since around 70 ka. Based on data presented in this paper, we argue that the Alban Hills should not be considered an extinct volcanic district and a detailed re-assessment of the volcanic hazard for the area of Rome is in order. Published online: 4 April 2003 Editorial responsibility: J. Donnelly-Nolan     

Numerical viscoelastic modelling by the spectral Laguerre method

Seismic wave propagation in a viscoelastic media can be described by a system of integro-differential equations. The solution of such equations requires special methods when using finite-difference techniques in the time domain. In the frequency domain, the integral terms are represented by complex elastic parameters. This paper presents an efficient algorithm for viscoelastic modelling based on the integral Laguerre transform for the approximation of temporal derivatives and for the calculation of convolution integrals. For the calculation of spatial derivatives, it is possible to use various methods: finite-difference and finite-element techniques, spectral and pseudo-spectral methods. We then obtain a system of algebraic equations with a matrix independent of the parameter m , i.e. the degree of the Laguerre polynomials. In this case, only the right-hand side of the system has recurrent dependence on the parameter m , which is an analogue of the temporal frequency in the frequency domain. The obtained system with a large number of right-hand sides can be solved using fast methods, where the matrix is transformed only once, as opposed to the frequency-domain approach, when the matrix is transformed for each temporal frequency.     

Frequency-dependent anisotropy due to meso-scale fractures in the presence of equant porosity

Fractured rock is often modelled under the assumption of perfect fluid pressure equalization between the fractures and equant porosity. This is consistent with laboratory estimates of the characteristic squirt-flow frequency. However, these laboratory measurements are carried out on rock samples which do not contain large fractures. We consider coupled fluid motion on two scales: the grain scale which controls behaviour in laboratory experiments and the fracture scale. Our approach reproduces generally accepted results in the low- and high-frequency limits. Even under the assumption of a high squirt-flow frequency, we find that frequency-dependent anisotropy can occur in the seismic frequency band when larger fractures are present. Shear-wave splitting becomes dependent on frequency, with the size of the fractures playing a controlling role in the relationship. Strong anisotropic attenuation can occur in the seismic frequency band. The magnitude of the frequency dependence is influenced strongly by the extent of equant porosity. With these results, it becomes possible in principle to distinguish between fracture- and microcrack-induced anisotropy, or more ambitiously to measure a characteristic fracture length from seismic data.     

Dublin days: the 2003 NAM

Mark E Bailey and Tom Ray give an overview of the National Astronomy Meeting in Dublin, where record numbers of astronomers went to talk, listen and network.     

盆地岩石圈结构与油气成藏及分布

本文综述了大陆岩石圈研究现状和克拉通盆地、裂谷盆地和前陆盆地的岩石圈结构特征,指出在古裂谷、褶皱带或区域性深断裂等陆壳构造薄弱带上发育起来的多期叠合盆地,具有很好的含油气前景。大型含油气盆地往往存在地幔上隆、地壳减薄和地壳内低速层,盆地基底沉降与盖层沉积厚度较大。适度的后期构造活动改造和岩浆活动有利于沉积盆地内油气生成与保存。     

He and Ne isotopes in oceanic crust: implications for noble gas recycling in the mantle

In an attempt to determine the helium and neon isotopic composition of the lower oceanic crust, we report new noble gas measurements on 11 million year old gabbros from Ocean Drilling Program site 735B in the Indian Ocean. The nine whole rock samples analyzed came from 20 to 500 m depth below the seafloor. Helium contents vary from 3.3×10⁻¹⁰ to 2.5×10⁻⁷ ccSTP/g by crushing and from 5.4×10⁻⁸ to 2.4×10⁻⁷ ccSTP/g by melting. ³He/⁴He ratios vary between 2.2 and 8.6 Ra by crushing and between 2.9 and 8.2 by melting. The highest R/Ra ratios are similar to the mean mid-ocean ridge basalt (MORB) ratio of 8±1. The lower values are attributed to radiogenic helium from in situ α-particle production during uranium and thorium decay. Neon isotopic ratios are similar to atmospheric ratios, reflecting a significant seawater circulation in the upper 500 m of exposed crust at this site. MORB-like neon, with elevated ²⁰Ne/²²Ne and ²¹Ne/²²Ne ratios, was found in some high temperature steps of heating experiments, but with very small anomalies compared to air. These first results from the lower oceanic crust indicate that subducted lower oceanic crust has an atmospheric ²⁰Ne/²²Ne ratio. Most of this neon must be removed during the subduction process, if the ocean crust is to be recirculated in the upper mantle, otherwise this atmospheric neon will overwhelm the upper mantle neon budget. Similarly, the high (U+Th)/³He ratio of these crustal gabbros will generate very radiogenic ⁴He/³He ratios on a 100 Ma time scale, so lower oceanic crust cannot be recycled into either MORB or oceanic island basalt without some form of processing.     

Computer Simulation of Shaped Charge Detonation Using Meshless Particle Method

Shaped charges are widely applied in the defense industry and civil engineering for penetrating hardened targets. Simulations of shaped charges are generally based on traditional grid-based numerical methods. Different from the existing literature, the SPH method is applied in this paper to simulate the detonation process of shaped charge explosive. The combination of the meshless, Lagrangian and particle nature inherent in the SPH method is very attractive in dealing with large deformations and large inhomogeneities in extremely transient situations such as detonation process. The SPH method and some related numerical aspects are briefly discussed first, followed by applications in simulating shaped charge detonations. The effects of different detonation cavity angles and different charge head lengths are investigated. The results demonstrate that major shaped charge detonation phenomena are captured with revealing observations.