Strontium, oxygen, and hydrogen isotope geochemistry of hydrothermally altered and weathered rocks in DSDP Hole 504B, Costa Rica Rift

DSDP Hole 504B was drilled into 6 Ma crust, about 200 km south of the Costa Rica Rift, Galapagos Spreading Center, penetrating 1.35 km into a section that can be divided into four zones—Zone I: oxic submarine weathering; Zone II: anoxic alteration; Zones III and IV: hydrothermal alteration to greenschist facies. In Zone III there is intense veining of pillow basalts. Zone IV consists of altered sheeted dikes. Isotopic geochemical signatures in relation to the alteration zones are recorded in Hole 504B, as follows:

Zone	Depth(m)	Average87Sr/86Sr	Average δ18O (%o)	Average δD (%o)
I	275–550	0.7032	7.3	−63
II	550–890	0.7029	6.5	−45
III	890–1050	0.7035	5.6	−31
IV	1050–1350	0.7032	5.5	−36

Full-size table

Alteration temperatures are as low as 10°C in Zones I and II based on oxygen isotope fractionation. Strontium isotopic data indicate that a circulation of seawater is much more restricted in Zone II than in Zone I. Fluid inclusion measurements of vein quartz indicate the alteration temperature was mainly 300 ± 20°C in Zones III and IV, which is consistent with secondary mineral assemblages.The strontium, oxygen, and hydrogen isotopic compositions of hydrothermal fluids which were responsible for the greenschist facies alteration in Zones III and IV are estimated to be 0.7037, 2‰, and 3‰, respectively. Strontium and oxygen isotope data indicate that completely altered portions of greenstones and vein minerals were in equilibrium with modified seawater under low water/rock ratios (in weight) of about 1.6. This value is close to that of the end-member hydrothermal fluids issuing at 21°N EPR.Basement rocks are not completely hydrothermally altered. About 32% of the greenstones in Zones III and IV have escaped alteration. Thus 1 g of fresh basalt including the 32% unaltered portion are required in order to make 1 g of end-member solution from fresh seawater in water-rock reactions.     

Transformation of phosphorus species in settling seston and during early sediment diagenesis

Sequential P extraction was combined with electron microscop and X-ray spectroscopy to characterise various P species and to study their transformation in settling seston and in recent sediment. During early diagenesis most of the particulate P formed in the water was redissolved. No net transformation into species that would resist dissolution was observed.It was shown that

Regional distribution and temporal trends of some contaminants in the water of the Baltic Sea

During a cruise on board RV Gauss in May/June 1988, joint investigations into organochlorine compounds, dissolved trace metals, petroleum hydrocarbons and basic hydrography were carried out at representative stations of the Baltic Monitoring Programme (BMP). The aim of the cruise was to study distribution patterns and — using previous data — to establish temporal trends if at all discernible.Each group of contaminants investigated showed specific characteristics, with differences even between compounds within the same group. The differences are due to:

Re-visiting large historical earthquakes in the Colombian Eastern Cordillera

A re-assessment of the historic seismicity of the central sector of the Colombian Eastern Cordillera (EC) is made by revision of bibliographic sources, by calibration with modern instrumental earthquakes, and by interpretations in terms of current knowledge of the tectonics and seismicity of the region. Throughout the process we have derived an equation to estimate M_w for shallow crustal earthquakes in Colombia using the length of isoseismal VIII, L_VIII:

A multidisciplinary effort to assign realistic source parameters to models of volcanic ash-cloud transport and dispersion during eruptions

During volcanic eruptions, volcanic ash transport and dispersion models (VATDs) are used to forecast the location and movement of ash clouds over hours to days in order to define hazards to aircraft and to communities downwind. Those models use input parameters, called “eruption source parameters”, such as plume height H, mass eruption rate Ṁ, duration D, and the mass fraction m₆₃ of erupted debris finer than about 4 or 63 μm, which can remain in the cloud for many hours or days. Observational constraints on the value of such parameters are frequently unavailable in the first minutes or hours after an eruption is detected. Moreover, observed plume height may change during an eruption, requiring rapid assignment of new parameters. This paper reports on a group effort to improve the accuracy of source parameters used by VATDs in the early hours of an eruption. We do so by first compiling a list of eruptions for which these parameters are well constrained, and then using these data to review and update previously studied parameter relationships. We find that the existing scatter in plots of H versus Ṁ yields an uncertainty within the 50% confidence interval of plus or minus a factor of four in eruption rate for a given plume height. This scatter is not clearly attributable to biases in measurement techniques or to well-recognized processes such as elutriation from pyroclastic flows. Sparse data on total grain-size distribution suggest that the mass fraction of fine debris m₆₃ could vary by nearly two orders of magnitude between small basaltic eruptions ( 0.01) and large silicic ones (> 0.5). We classify eleven eruption types; four types each for different sizes of silicic and mafic eruptions; submarine eruptions; “brief” or Vulcanian eruptions; and eruptions that generate co-ignimbrite or co-pyroclastic flow plumes. For each eruption type we assign source parameters. We then assign a characteristic eruption type to each of the world's  1500 Holocene volcanoes. These eruption types and associated parameters can be used for ash-cloud modeling in the event of an eruption, when no observational constraints on these parameters are available.     

Sur l'étude de l'intensité solaire en Égypte et en particulier de la région spectrale de courtes longueurs d'onde

Résumé Des mesures actinométriques de l'intensité solaire dans les domaines spectraux: <425, <525, >630 m ainsi que celle dans l'ensemble du spectre solaire (totale) ont été effectuées à Giza (lat. 30°N) à l'aide d'un actinomètre Linke-Feussner, pendant une année complète pour tous les jours le ciel étant parfaitement clair.D'après ces mesures on peut conclure que: l'intensité solaire totale en été est inférieure à sa valeur en hiver en raison de l'accroissement de la teneur en vapeur d'eau atmosphérique en été; la distribution spectrale de l'énergie solaire, au cours de l'année peut être considérée comme constante, pour une masse d'air déterminée; l'intensité solaire, arrivant au sol, dans la région spectrale en-dessous de =425 m peut donner une indication des fluctuations de l'épaisseur réduite de l'ozone.

---

Summary Actinometric measurements of the solar intensity in the spectral regions: <425, <525, >630 m as well as the total solar intensity for the whole spectrum of the sun have been undertaken at Giza (lat. 30°N) with a Linke-Feussner actinometer, for a whole year and for every day the sky was perfectly clear.From these measurements, one can conclude that: the total solar intensity in summer is less than that in winter because of the increase of atmospheric water vapor content; the spectral distribution of the solar energy, during the year could be considered constant for a certain air mass; the solar intensity, reaching the earth's surface, in the spectral region below =425 m can give an indication to the fluctuations in the reduced thickness of the ozone layer.

     

Ueber Herde elastischer Wellen in isotropen,homogenen Medien

Zusammenfassung 1) Es werden Multipollösungen der skalaren Wellengleichung ² f/t ² – c² div gradf=0 betrachtet. Einerseits kann man solche Lösungen direkt durch Kugelfunktionenn-ter Ordnung ausdrücken, anderseits aus der Einpollösungf=1/p F(t–p/c) durch Differentiation nachn Richtungen erhalten. Es wird der Zusammenhang zwischen den Ergebnissen der beiden Verfahren gezeigt. — 2) Für die Energiedichte und den Energiefluss durch Kugelflächen bei kleinen elastischen Verschiebungen werden Ausdrücke in Kugelkoordinaten angegeben. — 3) Für die Wellengleichung grad div –b ² rot rot werden rotationsfreie Multipollösungen angegeben und Ausdrücke für Energiedichte und Energiefluss hergeleitet. — 4) Das gleiche wird für divergenzfreie Multipollösungen durchgeführt. — 5) Es werden Multipole betrachtet, die weder rotationsfrei noch divergenzfrei sind. Als Spezialfälle werden Multipole mit zeitlich begrenzter und solche mit periodischer Erregung gezeigt, ferner Lösungen der Wellengleichung, die sowohl rotationsfrei wie divergenzfrei sind. — 6) Es wird gezeigt, wie man die elastischen Wellen, die im Sinne vonStokes von einem Herdgebiet endlicher Ausdehnung ausgehen, näherungsweise durch elastische Multipole darstellen kann. — 7) Es wird angedeutet, wie man durch Messung von Komponenten von oder u.s.w. in Punkten im Innern des Mediums die Erregung und Energie von elastischen Multipolen bestimmen kann. Ferner wird auf den Fall hingewiesen, wo ein rotationsfreier Einpol sich im Innern eines Halbraumes befindet und die Messungen an seiner Oberfläche ausgeführt werden.

---

Summary (On foci of elastic waves in isotropic homogeneous media) — 1) Multiplets as solutions of the scalar wave equation ² f/t ² – c² div gradf=0 are considered. Such solutions can be obtained either directly by aid of spherical harmonics of ordern, or by differentiating the single polef=1/p F(t–p/c) with respect ton directions. The relations between the results of those two procedures are shown. — 2) In the case of small elastic displacements , the density of energy and the flow of energy through spherical surfaces are expressed by spherical coordinates. — 3) Multiplets which satisfy the equation of motion =a ² grad div b ² curl curl and the equation curl = 0 are given, and expressions for the density and flow of energy are found. — 4) The same is done with multiplets satisfying the equation of motion and the equation div = 0. — 5) General multiplets which satisfy the equation of motion are treated. As special cases, multiplets with excitation of finite length and multiplets with periodic excitation are considered, furthermore solutions of the equation of motion and of the equations curl = 0 and div = 0 are given. — 6) It is shown how elastic waves whose origin is a region of finite extension in the sense given byStokes, can be approximated by elastic multiplets. — 7) Some indications are given on the problem of how to find the functions of excitation and the energy of an elastic multiplet by measuring components of or etc., at points in the interior of the medium. The same problem is considered in the case of the single elastic pole. = grad 1/p F (t–p/a), if the measurements are made at the surface of an elastic half space.

     

Organic thermometry for chondritic parent bodies

A unique spectroscopic feature has been identified in a study of twenty-five different samples of meteoritic insoluble organic matter (IOM) spanning multiple chemical classes, groups, and petrologic types, using carbon X-ray Absorption Near Edge Structure (XANES) spectroscopy. The intensity of this feature, a 1s − σ exciton, appears to provide a precise measure of parent body metamorphism. The intensity of this exciton is also shown to correlate well with a large negative paramagnetic shift observed through solid state ¹³C NMR. Experiments reveal that upon heating primitive IOM is transformed into material that is indistinguishable from that in thermally processed chondrites, including the development of the 1s − σ exciton. A thermo-kinetic expression is derived from the experimental data that allows the intensity of the 1s − σ exciton to be used to estimated the effective temperature integrated over time. A good correlation is observed between the intensity of the 1s − σ exciton and previously published microRaman spectral data. These data provide a self-consistent organic derived temperature scale for the purpose of calibrating Raman based thermometric expressions.     

Temperature variability in central Mexico and its possible association to solar activity

Minimum extreme temperature variability from five meteorological stations in the central part of Mexico covering a period from 1920 to 1990 is examined. We found a correlation coefficient (r=0.65) between these temperature records and geomagnetic activity. Furthermore, by performing spectral analysis peaks were obtained with similar periodicities to those found in the sunspot number, the magnetic solar cycle, cosmic ray fluxes and geomagnetic activity; all of these phenomena are modulated by solar activity. Signals with periodicities comparable to those observed in El Niño and the Quasi-Biennial Oscillation were also identified. We conclude that the solar signal is probably present in the minimum extreme temperature record of the central part of Mexico.     

Elastic properties of garnet solid-solution series

The elastic constants of sixteen garnet specimens of wide variety in chemical composition are accurately determined by means of the rectangular parallelpiped resonance method. The dependence of the elastic properties on chemical composition is analyzed using the present data and those for seven garnets investigated by other authors. The property X_i of a garnet solid solution i is given by a linear addition law in terms of the mole fraction n_ij of component j; X_i = Σn_ijX_j where the X_j's are the properties of the end-members j (j = pyrope, almandine, spessartine, grossular and andradite). The X_j's are determined for density ρ, bulk modulus K, and shear moduli C_s = (C₁₁ ? C₁₂)/2 and C₄₄. No systematic deviation is observed from the linear addition law for the elastic moduli nor for other quantities such as the elastic wave velocities. The extrapolated elastic moduli (Mbar) of the end-members are:

     

15.

L'ordre du coefficient de récombinaison dans la coucheE     

Kyryl Borissov Serafimov 《Pure and Applied Geophysics》1966,64(1):133-144

Résumé On propose une méthode de détermination de l'ordre du coefficient équivalent de récombinaison en n'utilisant que les données ionosphériques. A cette fin on évalue l'apport de la vitesse des changements passagersdN _m E/dt dans la densité électroniqueN _m E. L'évaluation se fait en présence d'une série de valeurs choisies de , on observe laquelle des courbes définis _m (t) dans cette étude coïncidera de plus près avec _m (t) mesurée. Les changements de montrés sur la figure 4 sont obtenus de cette manière. La valeur diurne de varie de (0.5÷1) 10^–7 cm³ sec^–1; peu après le lever et peu avant le coucher du soleil devient >10^–7 cm³ sec^–1. Lorsque cos>0, d'après la règle décrcit et devient 10^–8 cm³ sec^–1. Les explications des changements de obtenus de cette manière sont données en rendant compte des changements par rapport aux ions atomiques et moléculaires, à la disparition rapide des ions positifs pendant le coucher du soleil en présence d'un coefficient de récombinaison plus grand et des changements de température eventuels. La symétrie ou l'asymétrie des fréquences critiquesf ₀ E quand cos sont égaux permettent dans la marche diurne de juger de l'ordre de . Toutes ces évaluations indiquent également des valeurs de environ 10^–7 cm³ sec^–1. En précisant des mesures def ₀ E il est possible de définir non seulement l'ordre, mais aussi la valeur de elle-même.

---

Summary A method for the determination of the order of the equivalent recombination coefficient is suggested, by using ionospheric data only. The increase in the speed of temporary changes indN _m /dt in the electronic densityN _m E is estimated for this purpose. The estimation is done with series of selected for the sake of expediency values of , following at the same time which curves _m (t), determined in the course of work, will coincide most closely with the _m (t) measured. The changes in , shown in figure 4, have been obtained in this way. The diurnal value of is in the range of (0.5 to 1) 10^–7 cm³ sec, being >10^–7 cm³ sec a little after sunrise and a little prior to sunset. At cos<0, by rule decreases and becomes 10^–8 cm³ sec. Explanations of the changes thus obtained in are indicated, taking into consideration the relation of atomic and molecular ions, the rapid disappearance of positive ions at sunset with a higher recombination coefficient and eventual temperature changes. From the symmetry or asymmetry of the critical frequencies off ₀ E at equal cos in the course of the day it can also be judged for the order of . All those estimates show values of in the range of 10^–7 cm³ sec. In the case of precise measurements off ₀ E, not only the determination of the order but also the real value of is possible.

     

16.

The system Mg₂SiO₄MgOH₂O at high pressures and temperatures — Possible hydrous magnesian silicates in the mantle transition zone     

S. Akimoto  M. Akaogi 《Physics of the Earth and Planetary Interiors》1980,23(4):268-275

The system Mg₂SiO₄MgOH₂O was investigated at pressures between 85 and 160 kbar and at temperatures between 750 and 1200°C. In runs for a gel with Mg/Si ratio of 3 and with 4.0 wt.-percent H₂O, a dense hydrous magnesian silicate, denoted phase B by Ringwood and Major, was found at pressures from about 100 kbar to at least 160 kbar in the whole temperature range studied. In the following table the crystallographic parameters and chemical formula of phase B, determined in this study, are compared with those of the other dense hydrous silicates in Mg₂SiO₄MgOH₂O.

	Almandine	Pyrope	Spessartine	Grossular	Andradite
$\text{K}$	1.779 ± 0.008	1.730 ± 0.009	1.742 ± 0.009	1.691 ± 0.008	1.379 ± 0.017
$\text{C}{}_{\mathrm{<mtext>s</mtext>}}$	0.981 ± 0.004	0.925 ± 0.004	0.964 ± 0.004	1.106 ± 0.004	0.979 ± 0.007
$\text{C}{}_{\mathrm{<mtext>44</mtext>}}$	0.958 ± 0.005	0.919 ± 0.005	0.937 ± 0.005	1.017 ± 0.006	0.827 ± 0.010

     

17.

Differences between Mean Sea Levels for the Pacific,Atlantic and Indian Oceans From Topex/Poseidon Altimetry     

Burša  Milan  Kouba  Jan  Müller  Achim  Raděj  Karel  True  Scott A.  Vatrt  Viliam  Vojtíšková  Marie 《Studia Geophysica et Geodaetica》1999,43(1):1-6

Geopotential values W of the mean equipotential surfaces representing the mean ocean topography were computed on the basis of four years (1993 - 1996) TOPEX/POSEIDON altimeter data: W = 62 636 854.10m ² s ^–2 for the Pacific (P), W = 62 636 858.20m ² s ^–2 for the Atlantic (A), W = 62 636 856.28m ²s^–2 for the Indian (I) Oceans. The corresponding mean separations between the ocean levels were obtained as follows: A – P = – 42 cm, I– P = – 22 cm, I – A = 20 cm, the rms errors came out at about 0.3 cm. No sea surface topography model was used in the solution.     

18.

Outlier detection with partial information: application to emergency mapping     

Davide D’Alimonte  Dan Cornford 《Stochastic Environmental Research and Risk Assessment (SERRA)》2008,22(5):613-620

This paper, addresses the problem of novelty detection in the case that the observed data is a mixture of a known ‘background’ process contaminated with an unknown other process, which generates the outliers, or novel observations. The framework we describe here is quite general, employing univariate classification with incomplete information, based on knowledge of the distribution (the probability density function, pdf) of the data generated by the ‘background’ process. The relative proportion of this ‘background’ component (the prior ‘background’ probability), the pdf and the prior probabilities of all other components are all assumed unknown. The main contribution is a new classification scheme that identifies the maximum proportion of observed data following the known ‘background’ distribution. The method exploits the Kolmogorov–Smirnov test to estimate the proportions, and afterwards data are Bayes optimally separated. Results, demonstrated with synthetic data, show that this approach can produce more reliable results than a standard novelty detection scheme. The classification algorithm is then applied to the problem of identifying outliers in the SIC2004 data set, in order to detect the radioactive release simulated in the ‘joker’ data set. We propose this method as a reliable means of novelty detection in the emergency situation which can also be used to identify outliers prior to the application of a more general automatic mapping algorithm.

---

Phase	Composition	Space group	Cell parameters				Density (g cm?3)
			$\text{a}$	$\text{b}$	$\text{c}$	β
			(Å)	(Å)	(Å)	(°)
Chondrodite	Mg5Si2O10H2	$\text{P2}{}_1\text{/c}$	7.914	4.752	10.350	108.71	3.06
Clinohumite	Mg9Si4O18H2	$\text{P2}{}_1\text{/c}$	13.695	4.747	10.284	100.64	3.14
Phase A	Mg7Si2O14H6	$\text{P6}{}_3$	7.860		9.573		2.96
Phase B	Mg23Si8O42H6	$\text{P2}{}_1\text{/c}$	10.600	14.098	10.092	104.05	3.32

Davide D’AlimonteEmail:

Dan Cornford (Corresponding author)Email:

     

L'ordre du coefficient de récombinaison dans la coucheE

Résumé On propose une méthode de détermination de l'ordre du coefficient équivalent de récombinaison en n'utilisant que les données ionosphériques. A cette fin on évalue l'apport de la vitesse des changements passagersdN _m E/dt dans la densité électroniqueN _m E. L'évaluation se fait en présence d'une série de valeurs choisies de , on observe laquelle des courbes définis _m (t) dans cette étude coïncidera de plus près avec _m (t) mesurée. Les changements de montrés sur la figure 4 sont obtenus de cette manière. La valeur diurne de varie de (0.5÷1) 10^–7 cm³ sec^–1; peu après le lever et peu avant le coucher du soleil devient >10^–7 cm³ sec^–1. Lorsque cos>0, d'après la règle décrcit et devient 10^–8 cm³ sec^–1. Les explications des changements de obtenus de cette manière sont données en rendant compte des changements par rapport aux ions atomiques et moléculaires, à la disparition rapide des ions positifs pendant le coucher du soleil en présence d'un coefficient de récombinaison plus grand et des changements de température eventuels. La symétrie ou l'asymétrie des fréquences critiquesf ₀ E quand cos sont égaux permettent dans la marche diurne de juger de l'ordre de . Toutes ces évaluations indiquent également des valeurs de environ 10^–7 cm³ sec^–1. En précisant des mesures def ₀ E il est possible de définir non seulement l'ordre, mais aussi la valeur de elle-même.

---

Summary A method for the determination of the order of the equivalent recombination coefficient is suggested, by using ionospheric data only. The increase in the speed of temporary changes indN _m /dt in the electronic densityN _m E is estimated for this purpose. The estimation is done with series of selected for the sake of expediency values of , following at the same time which curves _m (t), determined in the course of work, will coincide most closely with the _m (t) measured. The changes in , shown in figure 4, have been obtained in this way. The diurnal value of is in the range of (0.5 to 1) 10^–7 cm³ sec, being >10^–7 cm³ sec a little after sunrise and a little prior to sunset. At cos<0, by rule decreases and becomes 10^–8 cm³ sec. Explanations of the changes thus obtained in are indicated, taking into consideration the relation of atomic and molecular ions, the rapid disappearance of positive ions at sunset with a higher recombination coefficient and eventual temperature changes. From the symmetry or asymmetry of the critical frequencies off ₀ E at equal cos in the course of the day it can also be judged for the order of . All those estimates show values of in the range of 10^–7 cm³ sec. In the case of precise measurements off ₀ E, not only the determination of the order but also the real value of is possible.

     

The system Mg2SiO4MgOH2O at high pressures and temperatures — Possible hydrous magnesian silicates in the mantle transition zone

The system Mg₂SiO₄MgOH₂O was investigated at pressures between 85 and 160 kbar and at temperatures between 750 and 1200°C. In runs for a gel with Mg/Si ratio of 3 and with 4.0 wt.-percent H₂O, a dense hydrous magnesian silicate, denoted phase B by Ringwood and Major, was found at pressures from about 100 kbar to at least 160 kbar in the whole temperature range studied. In the following table the crystallographic parameters and chemical formula of phase B, determined in this study, are compared with those of the other dense hydrous silicates in Mg₂SiO₄MgOH₂O.

Differences between Mean Sea Levels for the Pacific,Atlantic and Indian Oceans From Topex/Poseidon Altimetry

Geopotential values W of the mean equipotential surfaces representing the mean ocean topography were computed on the basis of four years (1993 - 1996) TOPEX/POSEIDON altimeter data: W = 62 636 854.10m ² s ^–2 for the Pacific (P), W = 62 636 858.20m ² s ^–2 for the Atlantic (A), W = 62 636 856.28m ²s^–2 for the Indian (I) Oceans. The corresponding mean separations between the ocean levels were obtained as follows: A – P = – 42 cm, I– P = – 22 cm, I – A = 20 cm, the rms errors came out at about 0.3 cm. No sea surface topography model was used in the solution.     

Outlier detection with partial information: application to emergency mapping

This paper, addresses the problem of novelty detection in the case that the observed data is a mixture of a known ‘background’ process contaminated with an unknown other process, which generates the outliers, or novel observations. The framework we describe here is quite general, employing univariate classification with incomplete information, based on knowledge of the distribution (the probability density function, pdf) of the data generated by the ‘background’ process. The relative proportion of this ‘background’ component (the prior ‘background’ probability), the pdf and the prior probabilities of all other components are all assumed unknown. The main contribution is a new classification scheme that identifies the maximum proportion of observed data following the known ‘background’ distribution. The method exploits the Kolmogorov–Smirnov test to estimate the proportions, and afterwards data are Bayes optimally separated. Results, demonstrated with synthetic data, show that this approach can produce more reliable results than a standard novelty detection scheme. The classification algorithm is then applied to the problem of identifying outliers in the SIC2004 data set, in order to detect the radioactive release simulated in the ‘joker’ data set. We propose this method as a reliable means of novelty detection in the emergency situation which can also be used to identify outliers prior to the application of a more general automatic mapping algorithm.

---

Viscosity of a Teide phonolite in the welding interval

The viscosity of a natural phonolitic composition with variable amounts of H₂O has been experimentally determined. The starting materials were crystal-free phonolitic glasses from Montaña Blanca, situated within the Las Cañadas caldera of Teide. Dry phonolitic melt viscosities were determined using concentric cylinder viscometry in the low viscosity range. The glassy quench products of these runs were then hydrated by high pressure synthesis in a piston–cylinder apparatus to generate a suite of samples with water contents ranging from 0.02 to 3.75 wt%. Samples thus hydrated were quenched rapidly and prepared (cut and polished) for the determination of water contents by infrared spectroscopy before and after experimental viscometry. The viscosities of the melts (dry and hydrated) were determined at 1 bar using a micropenetration technique. Samples were stable under the measurement conditions up to 3.75 wt% H₂O. Homogeneity of water content was confirmed by infrared spectroscopy and total water contents were calculated using absorptivity coefficients for compositions extremely close to that investigated here. The variation of viscosity as a function of water content and temperature can be described in the high viscosity interval of relevance to many welding processes by the non-Arrhenian expression:

(1)

log₁₀ η=−5.900−0.286 ln (H₂O)+(10775.4−394.8(H₂O))/(T−148.7+21.65 ln (H₂O))

whereas the high viscosity range alone is adequately described by the Arrhenian expression

(2)

log₁₀ η=−10.622−0.738 ln (H₂O)+(17114.3−590.4(H₂O))×1/T

where η is the viscosity in Pa s, H₂O is the water content in wt% and T is the temperature in K.These results are particularly useful for the scaling of conditions extant during the welding of phonolitic products of Montaña Blanca. The welding of glassy phonolitic rocks is enhanced by the lower viscosity of these melts with respect to calcalkaline rhyolites. The ratio of viscosities of phonolitic to calcalkaline rhyolitic melts is a complex function of temperature and water content and reaches up to 10^4.5 at 0.1 wt% H₂O and 500°C. Abundant evidence of welding and remobilisation of pyroclastic and spatter products of Teide system volcanism are consistent with these experimental observations.     

Contrasting compositions of Mariana Trough fallout tephra and Mariana Island arc volcanics: a fractional crystallization link

Discrete Quaternary (<400 ka) tephra fallout layers (mostly <1 cm thick) within the siliceous oozes of the central Mariana Trough at 18°N are characterized by medium-K to high-K subalkalic volcanic glasses (K₂O=0.8–3.2 wt.%) with high large-ion lithophile elements (LILE)/high-field-strength elements (HFSE) ratios and Nb depletion (Ba/La35; Ba/Zr3.5; La/Nb4) typical for convergent margin volcanic rocks. Compositional zoning within layers ranges from basaltic to dacitic (SiO₂=48–71 wt.%; MgO=0.7–6.5 wt.%); all layers contain basaltic andesites. The tephra layers are interpreted as single explosive eruptive events tapping chemically zoned reservoirs, the sources being the Mariana arc volcanoes (MAV) due to their proximity (100–400 km) and similar element ratios (MAV: Ba/La=36±7; Ba/Zr=3.5±0.9). The glasses investigated, however, contrast with the contemporaneous basaltic to dacitic lavas of the MAV by being more enriched in TiO₂ (1.2 wt.%; MAV0.8 wt.%), FeO^* (10 wt.%, MAV8–9 wt.%), K₂O (1.1 wt.%; MAV0.8 wt.%) and P₂O₅ (0.4 wt.%; MAV0.2 wt.%). (Semi-)Incompatible trace elements (including Rare Earth Elements (REE)) of the basaltic-andesitic and dacitic glasses match those of the dacitic MAV lavas, which became enriched by fractional crystallization. Moreover, the glasses follow a tholeiitic trend of fractionation in contrast to MAV transitional trends and have a characteristic P₂O₅ trend that reaches a maximum of 0.6 wt.% P₂O₅ at 57 wt.% SiO₂, whereas MAV lavas increase linearly in P₂O₅ from 0.1 to 0.3 wt.% with increasing silica. Both explosive and effusive series are interpreted to have evolved in common magma reservoirs by convective fractionation. Similar parental magmas are suggested to have separated into coexisting Si-andesitic to dacitic and basaltic melts by in situ crystallization. The differentiated melt is interstitial in an apatite-saturated crystalline mush of plag+px±ox±ol at the cooler chamber margins in contrast to the less differentiated basaltic to basaltic-andesitic magmas, which are not yet saturated in apatite and occupy the chamber interior. Reinjection of interstitial melt into the chamber interior and mixing with larger melt fractions of the interior liquid (mixing ratios about 1: 8–9) can explain the paradoxical behavior of apatite-controlled P and MREE variation in the basaltic andesite glasses and their MAV dacite-like fractionation patterns. The process may also account for the exclusively tholeiitic trend of fractionation of the glass shard series, but in situ crystallization alone cannot cause their absolute enrichment in (semi-)incompatible elements. The newly mixed melt is suggested to form the basaltic end member of the glass shard series. However, it must have become physically separated from the main MAV magma body (possibly by density-driven convective fractionation) in order to allow for further evolution of the contrasting geochemical paths as well as differentiation.