Geochemical modeling of magma mixing and magma reservoir volumes during early episodes of Kīlauea Volcano’s Pu‘u ‘Ō‘ō eruption

Geochemical modeling of magma mixing allows for evaluation of volumes of magma storage reservoirs and magma plumbing configurations. A new analytical expression is derived for a simple two-component box-mixing model describing the proportions of mixing components in erupted lavas as a function of time. Four versions of this model are applied to a mixing trend spanning episodes 3–31 of Kilauea Volcano’s Puu Oo eruption, each testing different constraints on magma reservoir input and output fluxes. Unknown parameters (e.g., magma reservoir influx rate, initial reservoir volume) are optimized for each model using a non-linear least squares technique to fit model trends to geochemical time-series data. The modeled mixing trend closely reproduces the observed compositional trend. The two models that match measured lava effusion rates have constant magma input and output fluxes and suggest a large pre-mixing magma reservoir (46±2 and 49±1 million m³), with little or no volume change over time. This volume is much larger than a previous estimate for the shallow, dike-shaped magma reservoir under the Puu Oo vent, which grew from ∼3 to ∼10–12 million m³. These volumetric differences are interpreted as indicating that mixing occurred first in a larger, deeper reservoir before the magma was injected into the overlying smaller reservoir. Electronic Supplementary Material Supplementary material is available at and is accessible for authorized users.     

Finding regions of influence on SO2 and SO= 4 daily concentration measurements at four sites in Spain

Contamination by the pollutants SO₂ and SO⁼₄ was analyzed for the 1989–1992 period at four regional stations in Spain under the auspices of the EMEP-BAPMON program. The evolution of the time series of the daily pollution has also been assessed, and high mean concentrations at La Cartuja and Logroño observed, with values of 3.8 and 4.5 g m⁻³ for SO₂, respectively. Maximum annual concentrations were recorded in 1989, when SO₂ reached values of 6.24, 5.39, 5.71, and 9.30 g m⁻³ for the stations of La Cartuja, San Pablo de los Montes, Roquetas, and Logroño, respectively. This work attempts to establish a relationship between the concentrations of the pollutants - both SO₂ gas and SO⁼₄ aerosol - and the zones of emission or persistence of these long-range transported pollutants. In this way, those regions showing a greater impact on the air quality in each season have been determined. To achieve this, the trajectories of the air masses carrying away the pollution to each of the receiving stations were considered and followed by a sectorial analysis. Nonparametric statistical methods were implemented to contrast the chemical homogeneity among the different sectors. The criterion that several homogeneous sectors form a chemically homogeneous region was used. To improve this sectorial analysis, we have proposed a new technique based on the Potential-Source-Contribution Function (PSCF). Starting out from a set of specified regions, considered to be chemically homogeneous domains, it is possible to determine the likelihood that an air mass with particular characteristics (e.g., that a value of the daily concentration higher than the mean recorded at the station has been obtained) will arrive at a given station after having crossed one of the previously defined regions. Using this technique, it is possible to determine the source regions through which the air masses circulate and bring high pollution concentrations to the studied stations. Thanks to the PSCF, these statistical methods offer, through a sectorial analysis, the possibility to pass from a qualitative to a more quantitative view.     

The Role of the Magnetic Field in the Ascension of the 22nd and 23rd Cycles of Activity

The first small active regions of the last two new cycles seem to be formed, in its very first, preparatory stage, close to the equator with uncertain magnetic orientation, and as secondary products at the periphery of equatorial `magnetic bubbles' of the new magnetic flux still belonging to the old cycle. At the real beginning of the new cycle, at the periphery of its higher-latitude `magnetic bubbles', the opposite applies to the secondary regions of the ending cycle. It seems that two modes of magnetic flux supply the visible layers of the photosphere during the cycle transition phase: the emergence of the magnetic flux from the depth of the convective zone, and the mutual interaction of the earlier emerging fluxes, or the action of the `local dynamo', giving rise to the new local magnetic field concentrations.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.     

GIOVE-A initial signal analysis

In late December 2005 the GIOVE-A test satellite was launched by the European Space Agency (ESA) to secure the frequencies for the Galileo system and to provide a platform for testing the new navigation signals. We performed an initial assessment of these signals using the 30 m deep space antenna of the DLR ground station in Weilheim (Germany). The antenna gain raised the signals above the noise level, thus allowing a detailed analysis even without knowledge of the ranging codes. The present paper covers the analysis of the L1/E1 signals, which includes a discussion of the spectrum, the time domain signal and a decoding of the spreading codes for the Open Service.     

Analysis of the precipitation and cloudiness associated with COLs occurrence in the Iberian Peninsula

Summary The Iberian Peninsula is one of the regions in the world with higher occurrence of cut-off low systems (COL). The aim of this paper is to analyse the weather events (rainfall and cloudiness layer) associated to COLs in the Iberian Peninsula with tools not previously used: (a) the use of the new multidecadal COLs database developed by Nieto et al (2005) that permit us to study a 41 years period (1958–1998), (b) the checking of the expected weather events (rainfall and cloudiness layer) associated with COLs in a conceptual model (Winkler et al, 2005) and (c) the extensive use of radiosoundings to analyse convective instability in areas inside and close to the COL. Two points of view are used to make the analysis: (1) a source oriented method, when a particular COL is followed and its associated precipitation and cloudiness is analysed over four quadrants in which Iberia was divided and (2) a receptor oriented method, when the precipitation associated to COLs is analysed in given areas, defined by patterns of precipitation. Results reveal that the precipitation and cloudiness patterns associated to COLs in the conceptual model reproduce quite well the main characteristics found over the Iberian Peninsula. The generalized idea that most of the COLs produce intense convective rainfall is show to be misleading. Convective phenomena are important usually when the centre of the COL is located on the Mediterranean region. Most of the rainfall associated with COLs comes from the baroclinic shield; specially in cases located over the west half of the Iberian Peninsula. It is shown that nearly 30% of COLs do not induce any rainfall; most of them located in the southern half of the Peninsula, and mainly during autumn. Only 30% of COLs produce generalized rainfall over the whole analysed territory, being most of them (about 90%) located over the western half of the Iberian Peninsula.     

Petrology and geochronology of basalt breccia from the 1996 earthquake swarm of Loihi seamount, Hawaii: magmatic history of its 1996 eruption

 Samples of basalt were collected during the Rapid Response cruise to Loihi seamount from a breccia that was probably created by the July to August 1996 Loihi earthquake swarm, the largest swarm ever recorded from a Hawaiian volcano. ²¹⁰Po–²¹⁰Pb dating of two fresh lava blocks from this breccia indicates that they were erupted during the first half of 1996, making this the first documented historical eruption of Loihi. Sonobuoys deployed during the August 1996 cruise recorded popping noises north of the breccia site, indicating that the eruption may have been continuing during the swarm. All of the breccia lava fragments are tholeiitic, like the vast majority of Loihi's most recent lavas. Reverse zoning at the rim of clinopyroxene phenocrysts, and the presence of two chemically distinct olivine phenocryst populations, indicate that the magma for the lavas was mixed just prior to eruption. The trace element geochemistry of these lavas indicates there has been a reversal in Loihi's temporal geochemical trend. Although the new Loihi lavas are similar isotopically and geochemically to recent Kilauea lavas and the mantle conduits for these two volcanoes appear to converge at depth, distinct trace element ratios for their recent lavas preclude common parental magmas for these two active volcanoes. The mineralogy of Loihi's recent tholeiitic lavas signify that they crystallized at moderate depths (∼8–9 km) within the volcano, which is approximately 1 km below the hypocenters for earthquakes from the 1996 swarm. Taken together, the petrological and seismic evidence indicates that Loihi's current magma chamber is considerably deeper than the shallow magma chamber (∼3–4 km) in the adjoining active shield volcanoes. Received: 21 August 1997 / Accepted: 15 February 1998     

Tableau de répartition stratigraphique des Grands Foraminifères caractéristiques du Crétacé moyen de la région Méditerranéenne

Le Groupe de Travail Européen des Grands Foraminifères présente un tableau de répartition stratigraphique de 42 espèces bien définies du Crétacé moyen de la région méditerranéenne. La répartition stratigraphique proposée pour chaque espèce est fondée soit sur des propres observations soit sur des données de la littérature supprimer et est contrôlée par la présence d'Ammonites ou de Foraminifères planctoniques.A summarizing account is presented of the deliberations of the research group for Large Foraminifers of the IGCP Project “Mid-Cretaceous Events”. Large Foraminifers are of incontestable value for dating carbonate platform sequences owing to the absence of many other diagnostic groups of organísms. A table of stratigraphical distributions for 42 species is presented.     

Crystal clots,amphibole fractionation and the evolution of calc-alkaline magmas

Petrographic and microprobe investigations of calc-alkaline (CA) rocks from the High Cascade Range (i.e., Mt. St. Helens, Mt. Jefferson, Crater Lake and Mt. Shasta) of western North America show that crystal clots represent primary igneous phase assemblages and are not products of amphibole reactions with melt. For each eruptive complex, crystal clots display diverse modal proportions even within a single eruptive unit. Nevertheless, in all cases the crystal-clot minerals are also represented in the rock as phenocrysts or microphenocrysts. Basalts contain clots of ol+plag+mgt, ol+mgt, cpx+ plag+mgt, cpx+mgt and plag+mgt; andesites, clots of cpx+mgt, opx+mgt, cpx+opx+plag+mgt, cpx+plag+mgt, opx+plag+mgt and plag±mgt; and dacites, clots of opx+mgt, cpx+opx+plag+ mgt, opx+plag+mgt, amph+plag+mgt±ilm, amph+mgt±ilm and plag±mgt. The bulk compositions of most of these clot assemblages could not have been derived from amphibole percursors. Although some amphiboles in dacitic rocks display a breakdown reaction of amph=plag+cpx+opx +mag, these mineral clusters, unlike those of clots, typically have a relict amphibole crystal outline and a fine-grained metamorphic texture. Plagioclase grains in the mineral clusters lack oscillatory zoning which is typical of crystal clot plagioclase grains. The euhedral to subhedral shapes of most clot minerals and the oscillatory zoning present in most clot plagioclase grains are not likely to have formed from the breakdown of amphibole. Crystal clots are also observed in Hawaiian and ocean floor basalts, although amphibole fractionation has not been proposed for those lavas. Magnetite fractionation may be the controlling process limiting iron enrichment in CA magmas rather than amphibole fractionation. Textural evidence indicates that magnetite is an early-forming phase in CA magmas. V, which is concentrated in magnetite, shows a strong decrease with increasing silica in many CA rocks, supporting a magnetite fractionation model.Hawaii Institute of Geophysics Contrib. No. 969     

Infrared spectrum and structural role of titanium in synthetic Ti-garnets

Ti-spessartite, Ti-andradite and its indium homolog, Ca₃In₂(Si, Ti)₃O₁₂, have been synthesized and investigated by infrared (ir) spectroscopy. The use of isotopic species (⁴⁶Ti-⁵⁰Ti) gives the unequivocal proof that some of the additional bands observed in the 800–600 cm^?1 region are due to TiO₄ tetrahedra. For Si-Ti replacements up to 10 mol%, the localization of Ti on the available tetrahedral sites depends on the nature of the cations. For the In garnet all (or nearly all) Ti is located on tetrahedral sites; in Ti-andradite Ti is distributed over tetrahedral and octahedral sites, tetrahedral sites being thus occupied, in part by Ti (ir band near 700 cm^?1) and in part by Fe (ir band near 650 cm^?1); and Ti-spessartite the presence of Ti on tetrahedral sites is doubtful, these sites being essentially occupied by Al.