Vent conditions for expected eruptions at Vesuvius

Determining consistent sets of vent conditions for next expected eruptions at Vesuvius is crucial for the simulation of the sub-aerial processes originating the volcanic hazard and the eruption impact. Here we refer to the expected eruptive scales and conditions defined in the frame of the EC Exploris project, and simulate the dynamics of magma ascent along the volcanic conduit for sub-steady phases of next eruptions characterized by intensities of the Violent Strombolian (VS), Sub-Plinian 2 (SP2), and Sub-Plinian 1 (SP1) scale. Sets of conditions for the simulations are determined on the basis of the bulk of knowledge on the past history of Vesuvius [Cioni, R., Bertagnini, A., Santacroce, R., Andronico, D., Explosive activity and eruption scenarios at Somma–Vesuvius (Italy): towards a new classification scheme. Journal of Volcanology and Geothermal Research, this issue.]. Volatile contents (H₂O and CO₂) are parameterized in order to account for the uncertainty in their expected amounts for a next eruption. In all cases the flow in the conduit is found to be choked, with velocities at the conduit exit or vent corresponding to the sonic velocity in the two-phase non-equilibrium magmatic mixture. Conduit diameters and vent mixture densities are found to display minimum overlapping between the different eruptive scales, while exit gas and particle velocities, as well as vent pressures, largely overlap. Vent diameters vary from as low as about 5 m for VS eruptions, to 35–55 m for the most violent SP1 eruption scale. Vent pressures can be as low as less than 1 MPa for the lowest volatile content employed of 2 wt.% H₂O and no CO₂, to 7–8 MPa for highest volatile contents of 5 wt.% H₂O and 2 wt.% CO₂ and large eruptive scales. Gas and particle velocities at the vent range from 100–250 m/s, with a tendency to decrease, and to increase the mechanical decoupling between the phases, with increasing eruptive scale. Except for velocities, all relevant vent quantities are more sensitive to the volatile content of the discharged magma for the highest eruptive scales considered.     

Lava flow hazard at Nyiragongo Volcano,DRC

Mt. Nyiragongo is one of the most dangerous volcanoes in the world for the risk associated with the propagation of lava flows. In 2002 several vents opened along a huge system of fractures, pouring out lava which reached and destroyed a considerable part of Goma, a town of about 500,000 inhabitants on the shore of Lake Kivu. In a companion paper (Favalli et al. in Bull Volcanol, this issue, 2008) we employed numerical simulations of probable lava flow paths to evaluate the lava flow hazard on the flanks of the volcano, including the neighbouring towns of Goma (DRC) and Gisenyi (Rwanda). In this paper we use numerical simulations to investigate the possibility of significantly reducing the lava flow hazard in the city through the construction of protective barriers. These barriers are added to the DEM of the area as additional morphological elements, and their effect is evaluated by repeating numerical simulations with and without the presence of barriers. A parametric study on barrier location, size, shape and orientation led to the identification of barriers which maximize protection while minimizing their impact. This study shows that the highest hazard area corresponding to eastern Goma, which was largely destroyed by lava flows in 2002, cannot be effectively protected from future lava flows towards Lake Kivu and should be abandoned. On the contrary, the rest of the town can be sheltered from lava flows by means of two barriers that deviate or contain the lava within the East Goma sector. A proposal for the future development of the town is formulated, whereby “new” Goma is completely safe from the arrival of lava flows originating from vents outside its boundaries. The proposal minimizes the risk of further destruction in town due to future lava flows.     

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.     

Optimal weakening and damping using polynomial control for seismically excited nonlinear structures

This paper presents an approach for the optimal design of a new retrofi t technique called weakening and damping that is valid for civil engineering inelastic structures.An alternative design methodology is developed with respect to the existing ones that is able to determine the locations and the magnitude of weakening and/or softening of structural elements and adding damping while insuring structural stability.An optimal polynomial controller that is a summation of polynomials in nonlinear states is used...     

Parameterisation of surface radiation flux at an Antarctic site

During the Antarctic summer 1994/95 the values of downward and upward flux densities of both solar and terrestrial radiation were recorded at 1200 m for a period of 1 month on the Reeves Nevè Glacier (lat 74°39′S, long 161°35′E), near the Italian base of Terra Nova Bay. The relations proposed by Swinbank [Swinbank, W.C., 1963. Long-wave radiation from clear skies. Q. J. R. Meteorol. Soc. 89, 339–348], Idso and Jackson [Idso, S.B., Jackson, R.D., 1969. Thermal radiation from the atmosphere. J. Geophys. Res. 74, 5397–5403] and Deacon [Deacon, E.L., 1970. The derivation of Swinbank's long-wave radiation formula. Q. J. R. Meteorol. Soc. 96, 313–319] associating the long-wave atmospheric radiation flux only to the air temperature at screen level were tested in extreme Antarctic climatological condition. A relation between the long-wave radiation flux and both screen air temperature and cloud cover fraction in accordance to the height of the cloud base was defined using the Kasten and Czeplak relationship that relates the solar radiation flux and the cloud cover index. The study of the incoming short-wave radiation flux from the clear sky and that reflected by the surrounding snow cover allowed for highlighting the role of surface geometry on the albedo measurements.     

The mass distribution in the innermost regions of spiral galaxies

We use high-spatial resolution (100 pc) rotation curves of 83 spiral galaxies to investigate the mass distribution of their innermost kpc. We show that, in this region, the luminous matter completely accounts for the gravitational potential and no dark component is required. The derived I-band disk mass-to-light ratios agree well with those obtained from population synthesis models and correlate with color in a similar way. We find strict upper limits of 10⁷ M for the masses of compact bodies at the center of spirals, ruling out that these systems host the remnants of the quasar activity.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

The role of constraints on joint hypocentre determination in seismotectonic studies

In this paper we show how the performance of the joint hypocentre determination (JHD) method can be improved, leading to reduced instability in cases close to singularity. The method has been extended by imposing a number of constraints introduced by other authors, and adding a new one. We tested the stability of the method and the relative advantages of the various constraints by simulating a geometrical space distribution of hypocentres recorded by a regional seismic network. We then applied this method to deep earthquakes that occurred in the Southern Tyrrhenian Sea subduction zone and to the seismicity of the Northern Apennines, which is moderately deeper than the typical shallow seismicity of the Apennines. The results obtained from the analysis of synthetic data and actual earthquakes confirm that the JHD method produces less scatter in the hypocentral determinations with respect to the standard locations. The main conclusion of our study is that we can significantly reduce the systematic mislocations that result from applying JHD to very clustered seismicity if we introduce the appropriate set of constraints.     

Redox-controlled iron isotope fractionation during magmatic differentiation: an example from the Red Hill intrusion, S. Tasmania

This study presents accurate and precise iron isotopic data for 16 co-magmatic rocks and 6 pyroxene–magnetite pairs from the classic, tholeiitic Red Hill sill in southern Tasmania. The intrusion exhibits a vertical continuum of compositions created by in situ fractional crystallisation of a single injection of magma in a closed igneous system and, as such, constitutes a natural laboratory amenable to determining the causes of Fe isotope fractionation in magmatic rocks. Early fractionation of pyroxenes and plagioclase, under conditions closed to oxygen exchange, gives rise to an iron enrichment trend and an increase in $ f_{{{\text{O}}_{2} }} $ of the melt relative to the Fayalite–Magnetite–Quartz (FMQ) buffer. Enrichment in Fe³⁺/ΣFe_melt is mirrored by δ⁵⁷Fe, where ^VIFe²⁺-bearing pyroxenes partition ⁵⁷Fe-depleted iron, defining an equilibrium pyroxene-melt fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{px}} - {\text{melt}}}} \le - 0.25\,\permille \times 10^{6} /T^{2} $ . Upon magnetite saturation, the $ f_{{{\text{O}}_{2} }} $ and δ⁵⁷Fe of the melt fall, commensurate with the sequestration of the oxidised, ⁵⁷Fe-enriched iron into magnetite, quantified as $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{melt}}}} = + 0.20\,\permille \times 10^{6} /T^{2} $ . Pyroxene–magnetite pairs reveal an equilibrium fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{px}}}} \approx + 0.30\,\permille $ at 900–1,000?°C. Iron isotopes in differentiated magmas suggest that they may act as an indicator of their oxidation state and tectonic setting.     

Effects of trace metal and exposure to air on serotonin and dopamine levels in tissues of the mussel Perna perna

We evaluated levels of serotonin (5HT) and dopamine (DOPA) in muscle and digestive glands of the mussel, Perna perna, collected at different times of day; exposed to air for 24 h, followed by re-submersion; and after exposure to different metals. Mussels collected at different periods of day showed little oscillation in 5HT and DOPA levels. Mussels exposed to metals showed significant changes in 5HT and DOPA levels in digestive gland and muscle, as did mussels exposed to air. Our data suggest that analyses of 5HT and DOPA in tissues of mussels could serve as a tool to evaluate the presence and effects of heavy metal contamination in mussels. Care in data interpretation is required, however, since other environmental factors such as exposure of mussels to air (i.e. at low tides) can also cause changes in DOPA and 5HT levels. Additional research is necessary to separate such natural environmental effects from effects of contaminants.