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.     

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.     

Lower crustal granulite xenoliths from the Pannonian Basin,Hungary, Part 2: Sr–Nd–Pb–Hf and O isotope evidence for formation of continental lower crust by tectonic emplacement of oceanic crust

Mafic granulite xenoliths from the lower crust of the Pannonian Basin are dominated by LREE-depleted bulk-rock compositions. Many of these have MORB-like ¹⁴³Nd/¹⁴⁴Nd but ⁸⁷Sr/⁸⁶Sr is elevated relative to most MORBs. Their '¹⁸O values cover a wide range from +3.8 to +9.5‰. A group of LREE-enriched mafic granulites have higher ⁸⁷Sr/⁸⁶Sr (0.704-0.708) and lower ¹⁴³Nd/¹⁴⁴Nd (0.5128-0.5124), with higher '¹⁸O values on average (+7.8 to +10.6‰) than the LREE-depleted granulites. The LREE-enriched granulites are, however, isotopically similar to newly discovered metasedimentary granulite xenoliths. A sublinear correlation in )Hf-)Nd isotope space has a shallower slope than the crust-mantle array, with the metasedimentary rocks forming the low )Hf end member; the radiogenic end is restricted to the LREE-depleted granulites and these overlap the field of MORB. Pb isotopes for the LREE-depleted samples are less radiogenic on average than those of the LREE-enriched and metasedimentary xenoliths, and metasedimentary granulites have consistently higher ²⁰⁸Pb/²⁰⁴Pb. The wide range in '¹⁸O over a restricted range in Nd and Sr isotope values, in combination with the predominance of LREE-depleted trace-element compositions, is consistent with an origin as a package of hydrothermally altered oceanic crust. The existence of '¹⁸O values lower than average MORB and/or mantle peridotite requires that at least some of these rocks were hydrothermally altered at high temperature, presumably in the oceanic lower crust. The low ¹⁴³Nd/¹⁴⁴Nd of the LREE-enriched mafic granulites cannot be explained by simple mixing between a LREE-depleted melt and an enriched component, represented by the recovered metasediments. Instead, we interpret these rocks as the metamorphic equivalent of the shallowest levels of the ocean crust where pillow basalts are intimately intercalated with oceanic sediments. A possible model is accretion of oceanic crustal slices during subduction and convergence, followed by high-grade metamorphism during the Alpine orogeny.     

A coupled,dynamical and chemical model for the prestellar core L1544: Comparison of modeled and observed C<Superscript>18</Superscript>O,HCO<Superscript>+</Superscript>, and CS emission spectra

We consider radiative transfer in C¹⁸O, HCO⁺, and CS molecular lines in a spherically symmetrical, coupled, dynamical and chemical model of a prestellar core whose evolution is determined by ambipolar diffusion. Theoretical and observed line profiles are compared for the well-studied core L1544, which may be a collapsing protostellar cloud. We study the relationship between the line shapes and model parameters. The structure of the envelope and kinematic parameters of the cloud are the most important factors determining the shape of the lines. Varying the input model parameters for the radiative transfer—the kinetic temperature and microturbulent velocity—within the limits imposed by observations does not result in any substantial variations of the line profiles. The comparison between the model and observed spectra indicates that L1544 displays a flattened structure, and is viewed at an oblique angle. A two-dimensional model is needed to reproduce this structure.     

Possible refugia for reefs in times of environmental stress

This paper investigates the refuge potential of (1) upwelling areas, (2) coral areas at medium depth, and (3) offshore bank and island reefs in a scenario of increased global warming, and thus increased sea surface temperature (SST) and increased solar UV radiation. (1) Observations on coral health and water temperature in the subtropical Atlantic (Eleuthera and Cat Island, Bahamas) and Indian Ocean (Sodwana Bay, South Africa) suggest a link between cool water delivered by upwelling and coral health. After the 1998 bleaching event, caused by strong SST anomalies, coral health and recovery from the previous year's bleaching was significantly better on the narrow southern Cat Island shelf (70% of corals healthy) where the presence of cold water was observed, which was attributed to small-scale upwelling, than on the wide northern Eleuthera shelf (44% of corals healthy), where downwelling of hot bank waters was believed to have damaged corals. In South Africa, regular, short-term upwelling events in five summers reduced SST to well below bleaching level. (2) In the northern Red Sea (Safaga Bay) and in South Africa (Sodwana Bay), wide areas with either coral frameworks or non-framework communities exist. Calculations show that if the top 10 m (20 m) of the ocean became inhospitable to corals, still 50.4% (17.5%) of the coral area would remain intact in the Red Sea and 99% (40%) in South Africa. (3) Offshore bank and island reefs investigated in the Turks, Caicos, and Mouchoir Banks and Grand and Little Cayman showed high rates of mortality and coral diseases. The most remote sites (Mouchoir Bank) were not the healthiest. Refuge areas appear to exist in (1) and (2), but in (3) only if vigorous water-circulation is encountered.