A Dislocation Model of Seismic Wave Attenuation and Micro-creep in the Earth: Harold Jeffreys and the Rheology of the Solid Earth

—A microphysical model of seismic wave attenuation is developed to provide a physical basis to interpret temperature and frequency dependence of seismic wave attenuation. The model is based on the dynamics of dislocation motion in minerals with a high Peierls stress. It is proposed that most of seismic wave attenuation occurs through the migration of geometrical kinks (micro-glide) and/or nucleation/migration of an isolated pair of kinks (Bordoni peak), whereas the long-term plastic deformation involves the continuing nucleation and migration of kinks (macro-glide). Kink migration is much easier than kink nucleation, and this provides a natural explanation for the vast difference in dislocation mobility between seismic and geological time scales. The frequency and temperature dependences of attenuation depend on the geometry and dynamics of dislocation motion both of which affect the distribution of relaxation times. The distribution of relaxation times is largely controlled by the distribution in distance between pinning points of dislocations, L, and the observed frequency dependence of Q, Q, Q∝ω^α is shown to require a distribution function of P(L)∝L ^-m with m=4-2α The activation energy of Q ^?1 in minerals with a high Peierls stress corresponds to that for kink nucleation and is similar to that of long-term creep. The observed large lateral variation in Q ^?1 strongly suggests that the Q ^?1 in the mantle is frequency dependent. Micro-deformation with high dislocation mobility will (temporarily) cease when all the geometrical kinks are exhausted. For a typical dislocation density of ～ 10⁸ m^?2, transient creep with small viscosity related to seismic wave attenuation will persist up to the strain of ～ 10^?6, thus even a small strain (～ 10^?6?10^?4) process such as post-glacial rebound is only marginally affected by this type of anelastic relaxation. At longer time scales continuing nucleation of kinks becomes important and enables indefinitely large strain, steady-state creep, causing viscous behavior.     

Melting of forsterite Mg2SiO4 up to 15 GPa

The melting curve of forsterite has been studied by static experiment up to a pressure of 15 GPa. Forsterite melts congruently at least up to 12.7 GPa. The congruent melting temperature is expressed by the Kraut-Kennedy equation in the following form: T_m(K)=2163 (1+3.0(V₀ ? V)/V₀), where the volume change with pressure was calculated by the Birch-Managhan equation of state with the isothermal bulk modulus K₀ = 125.4 GPa and its pressure derivative K′ = 5.33. The triple point of forsterite-β-Mg₂SiO₄-liquid will be located at about 2600°C and 20 GPa, assuming that congruent melting persists up to the limit of the stability field of forsterite. The extrapolation of the previous melting data on enstatite and periclase indicates that the eutectic composition of the forsterite-enstatite system should shift toward the forsterite component with increasing pressure, and there is a possibility of incongruent melting of forsterite into periclase and liquid at higher pressure, although no evidence on incongruent melting has been obtained in the present experiment.     

The peculiarities of magma rise in presence of water

The viscosity of basalts (quartz and olivine tholeiite) was studied under pressure in dry conditions and in the presence of water. In dry conditions at 1400°C when pressure increases to 20 kbars the viscosity reduces by a factor of 2. In conditions of water saturation of basalt melts at 5 kbar the viscosity is smaller by a factor of ～ 50 than that in dry conditions. In the water undersaturated conditions when water content is fixed (3.3% H₂O) in melt the viscosity considerably decreases with pressure and takes intermediate value between those in dry and water saturated conditions. Experimental data recently obtained permit us to consider the peculiarities of physical properties of magma in the presence of water on a new base. Ascending magma can reach critical velocities of transition to the turbulent regime under negligible pressure drop, as a result of low viscosity. It is known at present that water influences on the viscosity of acidic melt under pressure of 1–8 kbars and at temperatures between 800–1200°C. Various authors gave physico-chemical evaluation of the dynamics of granite melts on the basis of these data. The viscosity of basalt melts and their dynamics under normal pressure is also well-known. The known new experimental data of basaltic melt viscosity under pressure in dry conditions (Kushiro et al., 1976;Khitarov et al., 1978) and in the presence of water (Khitarov et al., 1976) embrace broader intervals of physico-chemical conditions as on the pressure (up to 20–30 kbar) as well on the content of water (from 3% up to 12 %). These data permitted to evaluate on a new base the dynamics of magmatic melts under pressure.     

Experimental studies of melt-peridotite reactions at 1–2 GPa and 1250–1400°C and their implications for transforming the nature of lithospheric mantle and for high-Mg signatures in adakitic rocks

Experiments of the melt-peridotite reaction at pressures of 1 and 2 GPa and temperatures of 1250–1400°C have been carried out to understand the nature of the peridotite xenoliths in the Mesozoic high-Mg diorites and basalts of the North China Craton,and further to elucidate the processes in which the Mesozoic lithospheric mantle in this region was transformed.We used Fuxin alkali basalt,Feixian alkali basalt,and Xu-Huai hornblende-garnet pyroxenite as starting materials for the reacting melts,and lherzolite xenoliths and synthesized harzburgite as starting materials for the lithospheric mantle.The experimental results indicate that:(1)the reactions between basaltic melts and lherzolite and harzburgite at 1–2 GPa and 1300–1400°C tended to dissolve pyroxene and precipitate low-Mg#olivine(Mg#=83.6–89.3),forming sequences of dunite-lherzolite(D-L)and duniteharzburgite(D-H),respectively;(2)reactions between hornblende-garnet pyroxenite and lherzolite at 1 GPa and 1250°C formed a D-H sequence,whereas reactions at 2 GPa and 1350°C formed orthopyroxenite layers and lherzolite;and(3)the reaction between a partial melt of hornblende-garnet pyroxenite and harzburgite resulted in a layer of orthopyroxenite at the boundary of the pyroxenite and harzburgite.The reacted melts have higher MgO abundances than the starting melts,demonstrating that the melt-peridotite reactions are responsible for the high-Mg#signatures of andesites or adakitic rocks.Our experimental results support the proposition that the abundant peridotite and pyroxenite xenoliths in western Shandong and the southern Taihang Mountains might have experienced multiple modifications in reaction to a variety of melts.We suggest that melt-peridotite reactions played important roles in transforming the nature of the Mesozoic lithospheric mantle in the region of the North China Craton.     

Great-circle Rayleigh wave attenuation and group velocity,Part IV: Regionalization and pure-path models for shear velocity and attenuation

Pure-path averages for group velocities and specific attenuation have been calculated from individual observations and from path averages for two regionalizations; one original to this study and the other previously devised by Wu. Both are based on four upper-mantle provinces: ocean basin, continent, island arc and mid-ocean ridge. Pure-path group velocities and specific attenuation have also been calculated for combinations of regions and provide well separated regional measurements for such composite regions.Shear-velocity models for pure and combined regions have been derived by a controlled Monte Carlo inversion procedure and indicates that a low-velocity zone is required beneath the oceans, but is not required beneath continents. Models have been produced for pure and combined ocean, ocean-ridge, continent and continent-arc provinces.Q^?1_R determined from pure-path average group velocities and attenuation coefficients has been regionalized successfully for 2- and 3-region combinations. The resulting pure-path Q^?1_R for continents is much lower than that for ocean basins and ocean-ridge provinces. Inversion of Q^?1_R for ocean-ridge provinces shows that the average Q_β for the upper 200 km of these regions is between 85 and 100.     

The geochemistry and petrogenesis of the volcanic rocks of Carriacou,Grenadine Islands,West Indies

Carriacou is small volcanic island located near the southern end of the Lesser Antillean chain. Field relationships, petrography and geochemistry of the Tertiary lavas, outcropping in the southern half of the island, are used to identify the rocks present and to determine their petrogenesis and assess their significance within the island arc.Six main volcanic units are present. From oldest to youngest, these are the clinopyroxene-phyric basalt (CPB) sequence, the amphibole-phyric andesite (APA) sequence, the clinopyroxene-megaphyric basalt (CMB) sequence, the olivine-microphyric basalt (OMB) sequence, the clinopyroxene phyric andesite (CPA) sequence, and the amphibole-megaphyric andesite (AMA) sequence. Volcaniclastic deposits are associated with the APA, CMB, and AMA sequences. The APA sequence is calcalkaline, whereas the other five sequences are tholeiitic.Sr isotope and rare earth element (REE) data suggest that these volcanic rocks were derived from partial melts of garnet-peridotite generated deep within the mantle. The OMB lavas have the highest temperature assemblages of intratelluric minerals and the least evolved chemical characteristics, and are considered to be closest in composition to a parental melt. Phenocryst assemblages and chemical variation suggest that the andesite sequences have been derived from the mafic melts by low pressure fractional crystallization of approximately 20% clinopyroxene and 20% olivine, plus smaller amounts of plagioclase and amphibole. The high concentrations of incompatible and compatible elements and the high⁸⁷Sr/⁸⁶Sr ratios may indicate that subduction is slower in the southern part of the arc, and fluids released during slab dehydration rich in Incompatible trace elements, in Radiogenic strontium, and in Silica (IRS), have modified the parental melts.     

Viscosity and structural changes of albite (NaAlSi3O8) melt at high pressures

Viscosity of anhydrous albite melt, determined by the falling-sphere method in the solid-media, piston-cylinder apparatus, decreases with increasing pressure from 1.13 × 10⁵ P at 1 atm to 1.8 × 10⁴ P at 20 kbar at 1400°C. The rate of decrease in viscosity is larger between 12 and 15 kbar than in other pressure ranges examined. The density of the quenched albite melt increases with increasing pressure of quenching from 2.38 g/cm³ at 1 atm to 2.53 g/cm³ at 20 kbar. The rate of increase in density is largest at pressures between 15 and 20 kbar. The melting curve of albite shows an inflexion at about 16 kbar. These observations strongly suggest that structural changes of albite melt would take place effectively at pressures near 15 kbar. Melt of jadeite (NaAlSi₂O₆) composition shows very similar changes in viscosity and density and a melting curve inflexion at pressures near 10 kbar. Difference in pressure for the suggested effective structural changes of albite and jadeite melts is 5–6 kbar, which is nearly the same as that between the subsolidus reaction curves nepheline + albite= 2jadeite and albite=jadeite + quartz. The structural changes of the melts are, however, continuous and begin to take place at pressures lower than those of the crystalline phases.     

High-temperature elasticity of the fluorite-structure compounds CaF2, SrF2 and BaF2

The elastic moduli of single-crystal CaF₂, SrF₂ and BaF₂ have been determined by the ultrasonic pulse superposition technique as a function of temperature from T = 298 to T = 650°K. These new data are consistent with other data obtained by ultrasonic pulse techniques in the region of room temperature and are superior to previous high-temperature data from resonance experiments. The elastic moduli (c) are represented by quadratic functions in T over the experimental temperature range with the curvature in the same sense for all the moduli. Evaluation of the temperature derivatives of the elastic moduli at constant volume indicates that the dominant temperature effect is extrinsic for (?K_S/?T)_P and intrinsic for (?μ/?T)_P, where K_S and μ are the isotropic bulk and shear moduli, respectively. For the series CaF₂SrF₂BaF₂, |(?c/?T)_p| decreases with increasing molar volume for all moduli; however there are no theoretical or empirical grounds on which to derive a simple relationship between (?c/?T)_P and crystallographic parameters.     

Shear modulus and Q of forsterite and dunite near partial melting from forced-oscillation experiments

An apparatus designed to determine the complex shear modulus of rock samples by forced torsion oscillations at high temperature and in the seismic frequency band 0.003–30 Hz is briefly described. Measurements were performed on natural dunite from Åheim (Norway) up to 1400°C and on polycrystalline forsterite up to 1500°C at 1 atm pressure. The two materials were chosen to study, by comparison, the effect of melt on the elasticity and anelasticity of mantle rocks.Between 1000 and 1200°C the absolute values of the shear modulus G are almost equal for both materials. Above 1200°C G for natural dunite decreases progressively with temperature and at 1400°C and 1 Hz reaches $\text{～}\text{1}\text{3}$ of its value at 1100°C. In contrast, G of pure forsterite depends little on temperature. For petrological reasons, supported by simultaneous measurements of the electric resistivity, there is strong evidence that the decrease of G in dunite above 1200°C is due to melt from the lower melting components of the dunite. Based on different models estimates of the melt fraction are made.At high temperature, in both materials Q^?1 is characterized by a monotonic decrease with frequency according to ω^?α, with α ≈ 0.25. An apparent activation energy of 38±5 kcal mol^?1 for forsterite and 48±8 kcal mol^?1 for dunite was found with no significant change in the regime of partial melting. From this it is concluded that Q^?1, even at partial melting, is dominated by solid state high temperature background absorption. There is no indication from these experiments for a constant-Q-band at low seismic frequencies or an increase of Q proportional to frequency as suggested by some seismologists. The present results are in good qualitative agreement with those for Young's modulus obtained previously by strain retardation experiments.     

Volume relaxation of compacted SiO2 glass: a model for the conservation of natural diaplectic glasses

Relaxation experiments on vitreous SiO₂, densified at 900°C, 20 kbar to a permanent density of 2.36 g cm^?3 yield an activation energy of69 ± 5kcal mol^?1 in the temperature range 700–800°C. The data can be used to estimate minimum cooling rates, maximum depositional temperatures or maximum thicknesses of shock-metamorphosed rocks containing diaplectic glasses, if relaxation of a glass densified under dynamic conditions is assumed to be similar to that of statically densified glass.     

Thermal equation of state of majorite with MORB composition

In situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press at SPring-8 on majoritic garnet synthesized from natural mid-ocean ridge basalt (MORB), whose chemical composition is close to the average of oceanic crust, at 19 GPa and 2200 K. Pressure-volume-temperature data were collected using a newly developed high-pressure cell assembly to 21 GPa and 1273 K. Data were fit to the high-temperature Birch-Murnaghan equation of state, with fixed values for the ambient cell volume (V₀ = 1574.14(4) Å³) and the pressure derivative of the isothermal bulk modulus (K^′T = 4). This yielded an isothermal bulk modulus of K_T0 = 173(1) GPa, a temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.022(5) GPa K⁻¹, and a volumetric coefficient of thermal expansivity α = a + bT with values of a = 2.0(3) × 10⁻⁵ K⁻¹ and b = 1.0(5) × 10⁻⁸ K⁻². The derived thermoelastic parameters are very similar to those of pyrope. The density of subducted oceanic crust compared to pyrolitic mantle at the conditions in Earth's transition zone (410-660 km depth) was calculated using these results and previously reported thermoelastic parameters for MORB and pyrolite mineral assembledges. These calculations show that oceanic crust is denser than pyrolitic mantle throughout the mantle transition zone along a normal geotherm, and the density difference is insensitive to temperature at the pressures in lower part of the transition zone.     

Apatite/liquid partition coefficients for the rare earth elements and strontium

Sixteen sets of apatite/liquid partition coefficients (D^ap/liq) for the rare earth elements (REE; La, Sm, Dy, Lu) and six values for Sr were experimentally determined in natural systems ranging from basanite to granite. The apatite + melt (glass) assemblages were obtained from starting glasses artificially enriched in REE, Sr and fluorapatite components; these were run under dry and hydrous conditions of 7.5–20 kbar and 950–1120°C in a solid-media, piston-cylinder apparatus. An SEM-equipped electron microprobe was used for subsequent measurement of REE and Sr concentrations in coexisting apatites and quenched glasses. The resulting partition coefficient patterns resemble previously determined apatite phenocryst/groundmass concentration ratios in the following respects: (1) the rare earth patterns are uniformly concave downward (i.e., the middle REE are more compatible in apatite than the light and heavy REE); (2) D_REE^ap/liq is much higher for silicic melts than for basic ones; and (3) strontium (and therefore Eu²⁺) is less concentrated by apatite than are the trivalent REE. The effects of both temperature and melt composition on D_REE^ap/liq are systematic and pronounced. At 950°C, for example, a change in melt SiO₂ content from 50 to 68 wt.% causes the average REE partition coefficient to increase from ～7 to ～30. A 130°C increase in temperature, on the other hand, results in a two-fold decrease in D_REE^ap/liq. Partitioning of Sr is insenstitive to changes in melt composition and temperature, and neither the Sr nor the REE partition coefficients appear to be affected by variations in pressure or H₂O content of the melt.The experimentally determined partition coefficients can be used not only in trace element modelling, but also to distinguish apatite phenocrysts from xenocrysts in rocks. Reported apatite megacryst/host basalt REE concentration ratios [12], for example, are considerably higher than the equilibrium partition coefficients, which suggest that in this particular case the apatite is actually xenocrystic.A reversal experiment incorporated in our study yielded diffusion profiles of REE in apatite, from which we extracted a REEαCa interdiffusion coefficient of 2–4×10^?14 cm²/s at 1120°C. Extrapolated downward to crustal temperatures, this low value suggests that complete REE equilibrium between felsic partial melts and residual apatite is rarely established.     

High-pressure disproportionation of Co2SiO4 spinel and implications for Mg2SiO4 spinel

Co₂SiO₄ spinel has been found to disproportionate into its isochemically mixed oxides with rocksalt and rutile structures at pressures between 170 and 190 kbar and temperatures between 1400 and 1800°C in a diamond-anvil press. The exact disproportionation pressure is not certain due to transient increases in pressure during the local and rapid heating by a continuous YAG laser. The slope of the phase boundary between the spinel phase and the mixed oxides is calculated to be?33 ± 20bar/deg. This negative slope is consistent with the observed anomalously large entropy of CoO (relative to its isostructural oxides) in entropy vs.(MV)^?1/2 systematics, whereM is the formula weight andV the molar volume. The sign of the slope for a phase boundary in the disproportionation of spinel depends on the values of entropy of the rocksalt oxides as well as the inverse character exhibited in the spinel phases. The normal entropy of MgO suggests that the phase boundary for the disproportionation of Mg₂SiO₄ spinel has positive slope.     

Frequency-Dependent Attenuation of Coda Waves in the Crust in Southwest Anatolia (Turkey)

The attenuation of coda waves in the earth’s crust in southwest (SW) Anatolia is estimated by using the coda wave method, which is based on the decrease of coda wave amplitude in time and distance. A total of 159 earthquakes were recorded between 1997 and 2010 by 11 stations belonging to the KOERI array. The coda quality factor Q _c is determined from the properties of scattered coda waves in a heterogeneous medium. Firstly, the quality factor Q ₀ (the value of Q _c at 1 Hz.) and its frequency dependency η are determined from this method depending on the attenuation properties of scattered coda waves for frequencies of 1.5, 3.0, 6.0, 8.0, 12 and 20 Hz. Secondly, the attenuation coefficients (δ) are estimated. The shape of the curve is controlled by the scattering and attenuation in the crustal volume sampled by the coda waves. The average Q _c values vary from 110 ± 15 to 1,436 ± 202 for the frequencies above. The Q ₀ and η values vary from 63 ± 7 to 95 ± 10 and from 0.87 ± 0.03 to 1.04 ± 0.09, respectively, for SW Anatolia. In this region, the average coda Q–f relation is described by Q _c = (78 ± 9)f ^0.98±0.07 and δ = 0.012 km^?1. The low Q ₀ and high η are consistent with a region characterized by high tectonic activity. The Q _c values were correlated with the tectonic pattern in SW Anatolia.     

Static compression of hydrous silicate melt and the effect of water on planetary differentiation

High pressure experiments using the sink/float method have bracketed the density of hydrous iron-rich ultrabasic silicate melt from 1.35 to 10.0 GPa at temperatures from 1400 to 1860 °C. The silicate melt composition was a 50–50 mixture of natural komatiite and synthetic fayalite. Water was added in the form of brucite Mg(OH)₂ and was present in the experimental run products at 2 wt.% and 5 wt.% levels as confirmed by microprobe analyses of total oxygen. Buoyancy marker spheres were olivines and garnets of known composition and density. The density of the silicate melt with 5 wt.% water at 2 GPa and 1500 °C is 0.192 g cm^? 3 less than the anhydrous form of this melt at the same P and T. This density difference gives a partial molar volume of water in silicate melt of ～ 7 cm³ mol^? 1, which is similar to previous studies at high pressure. The komatiite–fayalite liquids with 0 and 2 wt.% H₂O, have extrapolated density crossovers with equilibrium liquidus olivine at 8 and 9 GPa respectively, but there is no crossover for the liquid with 5 wt.% H₂O. These results are consistent with the hypothesis that dense hydrous melts could be gravitationally stable atop the 410 km discontinuity in the Earth. The results also support the notion that equilibrium liquidus olivine could float in an FeO-rich hydrous martian magma ocean. Extrapolation of the data suggests that FeO-rich hydrous melt could be negatively buoyant in the Earth's D″-region or atop the core–mantle-boundary (CMB), although experiments at higher pressure are needed to confirm this prediction.     

Short-period PKP phases and the anelastic mechanism of the inner core

Short-period seismograms are synthesized for PKP phases in anelastic Earth models. The synthetics were constructed using a synthetic technique valid at grazing incidence, a source-time function appropriate for deep-focus earthquakes, and an instrument response for either a short-period WWSSN or SRO seismograph. The agreement between predicted and observed amplitudes and spectral ratios requires neither a low-Q_α zone at 0.2–2 Hz nor a low or negative P-velocity gradient at the bottom of the outer core. Thin low-Q_α zones beneath the inner core boundary fit spectral ratio data that sample the upper 200 km of the inner core but fail to fit data that sample the lower inner core. Only a model having Q_α^?1?[0.003, 0.004] at 0.2–2 Hz, nearly constant with depth in the inner core, satisfies all of the spectral ratio and amplitude data. The assumption of a bulk viscosity of 10-10³ Pa s for the liquid phase of a partially molten inner core combined with the observation of low shear attenuation in the inner core at frequencies less than 0.005 Hz limit the physical parameters associated with two possible attenuation mechanisms: (1) fluid flow and viscous relaxation due to ellipsoidally shaped inclusions of melt, and (2) the solid-liquid phase transformation induced by the stress change during the passage of a seismic wave. Both mechanisms require an order of 0.1% partial melt to reproduce the observed Q_α^?1. In the outer core, the time constant of the mechanism of phase transformation is predicted to be 10⁴–10⁶ s. Confirmation of small shear attenuation in the inner core in the frequency band of seismic body waves would favor the mechanism of phase transformation.     

Experimental crystallization of chrome spinel in FAMOUS basalt 527-1-1

FAMOUS basalt 527-1-1 (a high-Mg oceanic pillow basalt) has three generations of spinel which can be distinguished petrographically and chemically. The first generation (Group I) have reaction coronas and are high in Al₂O₃. The second generation (Group II) have no reaction coronas and are high in Cr₂O₃ and the third generation (Group III) are small, late-stage spinels with intermediate Al₂O₃ and Cr₂O₃. Experimental synthesis of spinels from fused rock powder of this basalt was carried out at temperatures of 1175–1270°C and oxygen fugacities of 10^?5.5 to 10^?10 atm at 1 atm pressure. Spinel is the liquidus phase at oxygen fugacities of 10^?8.5 atm and higher but it does not crystallize at any temperature at oxygen fugacities less than 10^?9.5. The composition of our spinels synthesized at 1230–1250°C and 10^?9 atmf_O2 are most similar to the high-Cr spinels (Group II) found in the rock. Spinels synthesized at 1200°C and 10^?8.5 atm_O2 are chemically similar to the Group III spinels in 527-1-1. We did not synthesize spinel at any temperature or oxygen fugacity that are similar to the high-Al (Group I) spinel found in 527-1-1. These results indicate that the high-Cr (Group II) spinel is the liquidus phase in 527-1-1 at low pressure and Group III spinel crystallize below the liquidus (～1200°C) after eruption of the basalt on the sea floor. The high-Al spinel (Group I) could have crystallized at high pressure or from a magma enriched in Al and perhaps Mg compared to 527-1-1.     

Thermal diffusivity of rhyolitic glasses and melts: effects of temperature, crystals and dissolved water

Thermal diffusivity (D) was measured using laser-flash analysis on pristine and remelted obsidian samples from Mono Craters, California. These high-silica rhyolites contain between 0.013 and 1.10?wt% H₂O and 0 to 2?vol% crystallites. At room temperature, D _glass varies from 0.63 to 0.68?mm²?s^?1, with more crystalline samples having higher D. As T increases, D _glass decreases, approaching a constant value of ??0.55?mm²?s^?1 near 700?K. The glass data are fit with a simple model as an exponential function of temperature and a linear function of crystallinity. Dissolved water contents up to 1.1?wt% have no statistically significant effect on the thermal diffusivity of the glass. Upon crossing the glass transition, D decreases rapidly near ??1,000?K for the hydrous melts and ??1,200?K for anhydrous melts. Rhyolitic melts have a D _melt of ??0.51?mm²?s^?1. Thermal conductivity (k?=?D·??·C _P) of rhyolitic glass and melt increases slightly with T because heat capacity (C _P) increases with T more strongly than density (??) and D decrease. The thermal conductivity of rhyolitic melts is ??1.5?W?m^?1?K^?1, and should vary little over the likely range of magmatic temperatures and water contents. These values of D and k are similar to those of major crustal rock types and granitic protoliths at magmatic temperatures, suggesting that changes in thermal properties accompanying partial melting of the crust should be relatively minor. Numerical models of shallow rhyolite intrusions indicate that the key difference in thermal history between bodies that quench to obsidian, and those that crystallize, results from the release of latent heat of crystallization. Latent heat release enables bodies that crystallize to remain at high temperatures for much longer times and cool more slowly than glassy bodies. The time to solidification is similar in both cases, however, because solidification requires cooling through the glass transition in the first case, and cooling only to the solidus in the second.     

Geobarometry of ultramafic xenoliths from Loihi Seamount,Hawaii, on the basis of CO2 inclusions in olivine

Abundant fluid inclusions in olivine of dunite xenoliths (～1–3 cm) in basalt dredged from the young Loihi Seamount, 30 km southeast of Hawaii, are evidence for three coexisting immiscible fluid phases—silicate melt (now glass), sulfide melt (now solid), and dense supercritical CO₂ (now liquid + gas)—during growth and later fracturing of some of these olivine crystals. Some olivine xenocrysts, probably from disaggregation of xenoliths, contain similar inclusions.Most of the inclusions (2–10 μm) are on secondary planes, trapped during healing of fractures after the original crystal growth. Some such planes end abruptly within single crystals and are termed pseudosecondary, because they formed during the growth of the host olivine crystals. The “vapor” bubble in a few large (20–60 μm), isolated, and hence primary, silicate melt inclusions is too large to be the result of simple differential shrinkage. Under correct viewing conditions, these bubbles are seen to consist of CO₂ liquid and gas, with an aggregate ? = ～ 0.5–0.75 g cm^?3, and represent trapped globules of dense supercritical CO₂ (i.e., incipient “vesiculation” at depth). Some spinel crystals enclosed within olivine have attached CO₂ blebs. Spherical sulfide blebs having widely variable volume ratios to CO₂ and silicate glass are found in both primary and pseudosecondary inclusions, demonstrating that an immiscible sulfide melt was also present.Assuming olivine growth at ～ 1200°C and hydrostatic pressure from a liquid lava column, extrapolation of CO₂P-V-T data indicates that the primary inclusions were trapped at ～ 220–470 MPa (2200–4700 bars), or ～ 8–17 km depth in basalt magma of ? = 2.7 g cm^?3. Because the temperature cannot change much during the rise to eruption, the range of CO₂ densities reveals the change in pressure from that during original olivine growth to later deformation and rise to eruption on the sea floor. The presence of numerous decrepitated inclusions indicates that the inclusion sample studied is biased by the loss of higher-density inclusions and suggests that some part of these olivine xenoliths formed at greater depths.     

Borehole study of compressional and shear attenuation of basalt flows penetrated by the Brugdan and William wells on the Faroes shelf

We investigated the seismic attenuation of compressional (P‐) and converted shear (S‐) waves through stacked basalt flows using short‐offset vertical seismic profile (VSP) recordings from the Brugdan (6104/21–1) and William (6005/13–1A) wells in the Faroe‐Shetland Trough. The seismic quality factors (Q) were evaluated with the classical spectral ratio method and a root‐mean‐square time‐domain amplitude technique. We found the latter method showed more robust results when analysing signals within the basalt sequence. For the Brugdan well we calculated effective Q estimates of 22–26 and 13–17 for P‐ and S‐waves, respectively, and 25–33 for P‐waves in the William well. An effective Q_S/Q_P ratio of 0.50–0.77 was found from a depth interval in the basalt flow sequence where we expect fully saturated rocks. P‐wave quality factor estimates are consistent with results from other VSP experiments in the North Atlantic Margin, while the S‐wave quality factor is one of the first estimates from a stacked basalt formation using VSP data. Synthetic modelling demonstrates that seismic attenuation for P‐ and S‐waves in the stacked basalt flow sequence is mainly caused by one‐dimensional scattering, while intrinsic absorption is small.