The ultra‐high pressure coronitic and pseudomorphous reactions in a metagranodiorite from the Brossasco‐Isasca Unit,Dora‐Maira Massif,western Italian Alps: a petrographic study and equilibrium thermodynamic modelling

Metagranodiorite samples from the Brossasco‐Isasca Unit, Dora‐Maira Massif, western Alps, show pseudomorphous and coronitic textures where igneous minerals were partially replaced by ultra‐high pressure (UHP) metamorphic assemblages. The original magmatic paragenesis consisted of quartz, plagioclase, K‐feldspar, biotite and minor phases. During UHP metamorphism, the plagioclase (site P) was replaced by zoisite, jadeite, quartz, K‐feldspar and kyanite, and coronitic reactions developed between biotite and adjacent minerals. At the original igneous biotite–quartz contact (site A), a single corona of poorly zoned garnet is developed, whereas at the biotite–K‐feldspar (site B) and biotite–plagioclase (site C) contacts, composite coronas are formed. Integration of results from petrographic observations, calculations of mineral stoichiometry and thermodynamic calculations of mineral stability has allowed the determination of the metamorphic reactions involved and the estimation of the metamorphic conditions, which reached as high as 24 kbar and 650 °C. Accurate microanalysis by energy‐dispersive spectroscopy (EDS) and statistical analysis of the data allowed us to identify, for the first time in a natural Na‐pyroxene of metagranitoid rocks, the end‐member Ca‐Eskola.     

Continuous soil CO<Superscript>2</Superscript> and discrete plume SO<Subscript>2</Subscript> measurements at Mt. Etna (Italy) during 1997–2000: a contribution to volcano monitoring

Continuous monitoring of soil CO₂ dynamic concentration (which is proportional to the CO₂ flux through the soil) was carried out at a peripheral site of Mt. Etna during the period November 1997–September 2000 using an automated station. The acquired data were compared with SO₂ flux from the summit craters measured two to three times a week during the same period. The high frequency of data acquisition with both methods allowed us to analyze in detail the time variations of both parameters. Anomalous high values of soil CO₂ dynamic concentration always preceded periods of increased flux of plume SO₂, and these in turn were followed by periods of summit eruptions. The variations were modeled in terms of gas efflux increase due to magma ascent to shallow depth and its consequent depressurization and degassing. This model is supported by data from other geophysical and volcanological parameters. The rates of increase both of soil CO₂ dynamic concentration and of plume SO₂ flux are interpreted to be positively correlated both to the velocity of magma ascent within the volcano and to lava effusion rate once magma is erupted at the surface. Low rates of the increase were recorded before the nine-month-long 1999 subterminal eruption. Higher rates of increase were observed before the violent summit eruption of September-November 1999, and the highest rates were observed during shorter and very frequent spike-like anomalies that preceded the sequence of short-lived but very violent summit eruptions that started in late January 2000 and continued until late June of the same year. Furthermore, the time interval between the peaks of CO₂ and SO₂ in a single sequence of gas anomalies is likely to be controlled by magma ascent velocity.Editorial responsibility: H. Shinohara     

Defining a model of 3D seismogenic sources for Seismic Hazard Assessment applications: The case of central Apennines (Italy)

Geology-based methods for Probabilistic Seismic Hazard Assessment (PSHA) have been developing in Italy. These methods require information on the geometric, kinematic and energetic parameters of the major seismogenic faults. In this paper, we define a model of 3D seismogenic sources in the central Apennines of Italy. Our approach is mainly structural-seismotectonic: we integrate surface geology data (trace of active faults, i.e. 2D features) with seismicity and subsurface geological–geophysical data (3D approach). A fundamental step is to fix constraints on the thickness of the seismogenic layer and deep geometry of faults: we use constraints from the depth distribution of aftershock zones and background seismicity; we also use information on the structural style of the extensional deformation at crustal scale (mainly from seismic reflection data), as well as on the strength and behaviour (brittle versus plastic) of the crust by rheological profiling. Geological observations allow us to define a segmentation model consisting of major fault structures separated by first-order (kilometric scale) structural-geometric complexities considered as likely barriers to the propagation of major earthquake ruptures. Once defined the 3D fault features and the segmentation model, the step onward is the computation of the maximum magnitude of the expected earthquake (M _max). We compare three different estimates of M _max: (1) from association of past earthquakes to faults; (2) from 3D fault geometry and (3) from geometrical estimate corrected by earthquake scaling laws. By integrating all the data, we define a model of seismogenic sources (seismogenic boxes), which can be directly used for regional-scale PSHA. Preliminary applications of PSHA indicate that the 3D approach may allow to hazard scenarios more realistic than those previously proposed.     

Investigation of insured earthquake damage

The assessment of the loss potential caused by natural perils is a very important task for all insurance companies working in hazard-prone markets. It has to be based on two crucial items: the frequency of events and the investigation of their effects on the insured portfolio.This article deals with the second aspect, i.e. an evaluation of the insured damage caused by two earthquakes, namely those occurring near Albstadt, Germany, on 3 September 1978, and in central Chile on 3 March 1985. The results of the analysis of the earthquake in central Chile enable the mean damage ratio (damage in relation to the value) to be related to the height and the type of construction of the buildings affected. The Albstadt earthquake data permit an illustration of the effects of the type of subsoil on the mean damage ratio. The damage to individual buildings can be described by a lognormal distribution. Possible applications of these results are mentioned.     

In-situ observations of the latitudinal gradients of the solar wind parameters during 1976 and 1977

Interplanetary observations from Helios 1, Helios 2, and IMP-8 spacecraft during 1976 and 1977, namely the early portion of solar cycle 21, have been used to investigate the latitudinal gradients of the solar wind parameters with respect to the angular displacement from the current sheet inferred from synoptic HAO white-light maps of the solar corona at 1.75 solar radii. A latitudinal belt of ±25 deg around the current sheet has been investigated. Large gradients for solar wind flow speed, proton density and temperature have been found. Smoother gradients were also found for particle flux, kinetic, gravitational and thermal energy density flux. All these gradients revealed to become smoother going towards the solar cycle's maximum. Neither latitudinal nor temporal variations were identified for magnetic and thermal energy density. A remarkable result of this study is that the momentum flux density and the total energy flux density which other authors found to be independent of any longitudinal stream structure were also found to be independent of any latitudinal structure. Moreover, these two parameters did not show any temporal variation during the period of interest.     

A seasonal model of the Saturnian upper troposphere: Comparison with Voyager infrared measurements

An extension of the seasonal climate model of R. D. Cess and J. Caldwell (1979, Icarus, 38, 349–357) to Saturn's upper troposphere is presented. The ring-modulated latitudinal dependence of the insolation, the ring thermal emission, the oblateness of the planet, the orbit eccentricity, and the latitudinal variation of the internal heat flux are taken into account. Calculations agree closely with the temperature—latitude profiles retrieved from Voyager IRIS measurements at atmospheric levels located above the 0.2-bar pressure level; they reproduce the observed large-scale hemispheric asymmetry which is then shown to result from the seasonally variable insolation. Aerosol absorption is found to be the dominant source of atmospheric solar heating in the troposphere and the model suggests an aerosol mean unit optical depth around the 0.25-bar level in the equatorial region and around the 0.35-bar level at other latitudes. The model fails to predict the retrieved temperature—latitude profiles below the 0.3-bar level. This discrepancy is attributed to the existence of clouds at these levels which are responsible for an additional far-infrared opacity not taken into account in the temperature retrieval. The cloud-top altitude would be about 0.3 bar except in the 20 to 40°N region where these clouds would be confined below the 0.6-bar level. The poor correlation between infrared measurements and visible images is discussed and a possible model of Saturn's cloud structure is proposed.     

The major features of the crustal structure in north-eastern Venezuela from deep wide-angle seismic observations and gravity modelling

Venezuela is located on the plate boundary zone between the South American continent and the Caribbean plate. A relative movement of 2 cm/year is accommodated by a system of strike–slip faults running from the Andes to the Gulf of Paria. The Interior Range, a moderate-height mountain range, separates the Oriental Basin from the Caribbean. To the south, predominantly Precambrian rocks are outcropping in the Guayana Shield south of the Orinoco River. Results of deep wide-angle seismic measurements for the region were obtained during field campaigns in 1998 (ECOGUAY) for the Guayana Shield and in 2001 (ECCO) for the Oriental Basin. The total crustal thickness decreases from 45 km beneath the Guayana Shield, to 39 km at the Orinoco River, and 36 km close to El Tigre, in the center of the Oriental Basin. The average crustal velocity decreases in the same sense from 6.5 to 5.95 km/s. Detailed information was obtained on the velocity distribution within the Oriental Basin. Velocities are as low as 2.2 km/s for the uppermost 2 km, 4.5 km/s down to 4 km in depth, and a maximum depth of 13 km was derived for material with seismic velocities up to 5.9 km/s, interpreted as the base of the sedimentary basin. A gravimetric model confirms the structures derived from the seismic data. Discrete increases in sedimentary thickness along the basin may be associated to extension processes during the passive margin phase in the Cretaceous, or during earlier extension phases.     

Effects of buoyancy and mechanical layering on collisional deformation of continental lithosphere: Results from physical modeling

We use two suites of lithospheric-scale physical experiments to investigate the manner in which deformation of the continental lithosphere is affected by both (1) variations of lithospheric density (quantified by the net buoyant mass per area in the lithospheric mantle layer, M_B), and (2) the degree of coupling between the crust and lithospheric mantle (characterized by a modified Ampferer ratio, Am). The dynamics of the experiments can be characterized with a Rayleigh–Taylor type ratio, _CLM. Models with a positively buoyant lithospheric mantle layer (M_B > 0 and _CLM > 0) result in distributed root formation and a wide deformation belt. In contrast, models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust (M_B < 0, 0 > _CLM > ≈ − 0.2, and Am > ≈ 10^{− 3}) exhibit localized roots and narrow deformation belts. Syncollisional delamination of the model lithospheric mantle layer and a wide deformation belt is exhibited in models with negatively buoyant lithospheric mantle layers weakly coupled to the crust (M_B < 0, _CLM < 0, and Am < ≈ 10^{− 3}). Syncollisional delamination of the continental lithosphere may initiate due to buoyancy contrasts within the continental plate, instead of resulting from wedging by the opposing plate. Rayleigh–Taylor instabilities dominate the style of deformation in models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust and a slow convergence rate (M_B < 0 and _CLM > ≈ − 0.2). The degree of coupling (Am) between the model crust and lithospheric mantle plays a lesser role in both the style of lower-lithospheric deformation and the width of the crustal deformed zone with increasing density of the lithospheric mantle layer.     

Boron isotopic fractionation between minerals and fluids: New insights from in situ high pressure-high temperature vibrational spectroscopic data

Equilibrium boron isotopic fractionations between trigonal B(OH)₃ and tetragonal B(OH)₄⁻ aqueous species have been calculated at high P-T conditions using measured vibrational spectra (Raman and IR) and force-field modeling to compute reduced partition function ratios for B-isotopic exchange following Urey’s theory. The calculated isotopic fractionation factor at 300 K, α_3/4 = 1.0176(2), is slightly lower than the formerly calculated value of α_3/4 = 1.0193 (Kakihana and Kotaka, 1977), due to differences in the determined vibrational frequencies. The effect of pressure on α_3/4 up to 10 GPa and 723 K is shown to be negligible relative to temperature or speciation (pH) effects. Implications for the interpretation of boron fractionation in experimental and natural systems are discussed. We also show that the relationship between seawater-mineral B isotope fractionation and pH can be expressed using two variables, α_3/4 on one hand, and the pK_a of the boric acid-borate equilibrium on the other hand. This latter value is given by the equilibrium of boron species in water for the carbonate-water exchange, but could be governed by mineral surface properties in the case of clays. This may allow defining intrinsic paleo-pHmeters from B isotope fractionation between carbonate and authigenic minerals. Finally, it is shown that fractionation of boron isotopes can be rationalized in terms of the changes in 1) coordination of B from trigonal to tetrahedral in both fluids and minerals; and 2) the ligand nature around B from OH⁻ in the fluid and some hydrous minerals to non-hydrogenated O in many minerals. Relationships are established that allow predicting the isotopic fractionation factor of B between minerals and fluid.