Source parameters of some recent earthquakes in the Gulf of Aqaba, Egypt

Gulf of Aqaba is recognized as an active seismic zone where many destructive earthquakes have occurred. The estimation of source parameters and coda Q attenuation are the main target of this work. Fifty digital seismic events in eight short-period seismic stations with magnitude 2.5–5.2 are used. Most of these events occurred at hypocentral depths in the range of 7–20 km, indicating that the activity was restricted in the upper crust. Seismic moment, M _o, source radius, r, and stress drop, Δσ, are estimated from P- and S-wave spectra using the Brune’s seismic source model. The average seismic moment generated by the whole sequence of events was estimated to be 4.6E?+?22 dyne/cm. The earthquakes with higher stress drop occur at 10-km depth. The scaling relation between the seismic moment and the stress drop indicates a tendency of increasing seismic moment with stress drop. The seismic moment increases with increasing the source radius. Coda waves are sensitive to changes in the subsurface due to the wide scattering effects generating these waves. Single scattering model of local earthquakes is used to the coda Q calculation. The coda with lapse times 10, 20, and 30 s at six central frequencies 1.5, 3, 6, 12, 18, 24 Hz are calculated. The Q _c values are frequency dependent in the range 1–25 Hz, and are approximated by a least squares fit to the power law [ $ {Q_c}(f) = {Q_o}{(f/{f_o})^\eta } $ ]. The average of Q _c values increases from 53?±?10 at 1.5 Hz to 700?±?120 at 24 Hz. The average of Q _o values ranges from 13?±?1 at 1.5 Hz to 39?±?4 at 24 Hz. The frequency exponent parameter η ranges between 1.3?±?0.008 and 0.9?±?0.001.     

Strain amplitude-dependent attenuation of P and S waves in dry and water-saturated sandstone under confining pressure

The dependence of P- and S-wave attenuation on strain amplitude and frequency has been studied experimentally in dry and water-saturated sandstone samples under a confining pressure of 20 MPa. Attenuation of P and S reflections was measured at a frequency of 1 MHz in a strain range of ε ～(0.3?2.0) × 10⁶. The measured P-wave attenuation (Q_P^?1) in dry sandstone and S-wave attenuation (Q_S^?1) in dry and saturated samples turned out to be inversely proportional to strain amplitude while Q_P^?1 in saturated sandstone showed no strain dependence. The frequency-dependent attenuation spectra in dry and saturated sandstone differed considerably in S waves but were generally similar for P reflections. Strain amplitude variations were found out to influence the frequency dependence of attenuation and to shift the relaxation spectra of S waves toward high frequencies relative to those of P waves. As strain increased, the S-wave attenuation peak in saturated sandstone became notably (to 40%) narrower. The unusual strain amplitude-dependent behavior of attenuation may be due to joint action of viscoelastic and microplastic mechanisms. The reported results may be useful for improving geological interpretation of acoustic and seismic data.     

Source parameters of earthquakes in the reservoir-triggered seismic (RTS) zone of Koyna–Warna, Western India

New empirical relations are derived for source parameters of the Koyna–Warna reservoir-triggered seismic zone in Western India using spectral analysis of 38 local earthquakes in the magnitude range M _L 3.5–5.2. The data come from a seismic network operated by the CSIR-National Geophysical Research Institute, India, during March 2005 to April 2012 in this region. The source parameters viz. seismic moment, source radius, corner frequency and stress drop for the various events lie in the range of 10¹³–10¹⁶ Nm, 0.1–0.4 km, 2.9–9.4 Hz and 3–26 MPa, respectively. Linear relationships are obtained among the seismic moment (M ₀), local magnitude (M _L), moment magnitude (M _w), corner frequency (fc) and stress drop (?σ). The stress drops in the Koyna–Warna region are found to increase with magnitude as well as focal depths of earthquakes. Interestingly, accurate depths derived from moment tensor inversion of earthquake waveforms show a strong correlation with the stress drops, seemingly characteristic of the Koyna–Warna region.     

Raman spectroscopy of natural silica in Chicxulub impactite,MexicoSpectroscopie Raman de silices naturelles dans l'impactite de Chicxulub,Mexique

A series of natural silica impactite samples from Chicxulub (Mexico) was investigated by Raman microprobe (RMP) analysis. The data yield evidence for high-pressure shock metamorphism in the rock. The impactite contains three polymorphs of silica: the original α-quartz, and two high-pressure varieties – coesite and disordered quartz representing various degrees of crystallinity. We found systematic changes in frequencies and half-widths of the Raman bands, caused by increasing irregularities of bond-lengths and bond-angles and a general breaking-up of the structure as a result of impact events. Therefore, RMP is an adequate tool for measuring the crystallinity of disordered quartz. The half-width Γ and the frequency ω of the symmetric SiOSi stretching vibrational band (A₁ mode) of the SiO₄ tetrahedra are the most amenable parameters for estimating the degree of crystallinity. In well-crystallized quartz, Γ=5 cm^?1 and ω=464 cm^?1, while in highly disordered quartz this line shifts up to ω=455 cm^?1 and broadens up to Γ=30 cm^?1. The Raman lineshapes appear to depend strongly on the degree of lattice disorder subsequent to impact events. To cite this article: M. Ostroumov et al., C. R. Geoscience 334 (2002) 21–26     

Source parameters and scaling relations for small earthquakes in the Kachchh region of Gujarat,India

The scaling relationships for stress drop and corner frequency with respect to magnitude have been worked out using 159 accelerograms from 34 small earthquakes (M _w 3.3–4.9) in the Kachchh region of Gujarat. The 318 spectra of P and S waves have been analyzed for this purpose. The average ratio of P- to S-wave corner frequency is found to be 1.19 suggestive of higher corner frequency for P wave as compared to that for S wave. The seismic moments estimated from P waves, M ₀(P), range from 1.98 × 10¹⁴ N m to 1.60 × 10¹⁶ N m and those from S waves, M ₀(S), range from 1.02 × 10¹⁴ N m to 3.4 × 10¹⁶ N m with an average ratio, M ₀(P)/M ₀(S), of 1.11. The total seismic energy varies from 1.83 × 10¹⁰ J to 2.84 × 10¹³ J. The estimated stress drop values do not depend on earthquake size significantly and lie in the range 30–120 bars for most of the events. A linear regression analysis between the estimated seismic moment (M ₀) and corner frequency (f _c) gives the scaling relation M ₀ f _c ³  = 7.6 × 10¹⁶ N m/s³. The proposed scaling laws are found to be consistent with similar scaling relations obtained in other seismically active regions of the world. Such an investigation should prove useful in seismic hazard and risk-related studies of the region. The relations developed in this study may be useful for the seismic hazard studies in the region.     

Use of site amplification and anelastic attenuation for the determination of source parameters of the Sikkim earthquake of September 18, 2011, using far-field strong-motion data

An earthquake of magnitude 6.9 (M _w) occurred in the Sikkim region of India on September 18, 2011. This earthquake is recorded on strong-motion network in Uttarakhand Himalaya located about 900 km away from the epicenter of this earthquake. In this paper acceleration record from six far-field stations has been used to compute the source parameters of this earthquake. The acceleration spectra of ground motion at these far-field stations are strongly affected by both local site effects and near-site anelastic attenuation. In the present work the spectrum of S-phase recorded at these far-field stations has been corrected for anelastic attenuation at both source and site and the site amplification terms. Site amplifications at different stations and near-site shear wave attenuation factor have been computed by the technique of inversion of acceleration spectra given by Joshi et al. (Pure Appl Geophys 169:1821–1845, 2012a). For estimation of site amplification and shear wave quality factor [Q _β (f)] at the recording sites, ten local events recorded at various stations between July 2011 and December 2011 have been used. The obtained source spectrum from acceleration records is compared with the theoretical source spectrum defined by Brune (J Geophys Res 76:5002, 1970) at each station for both horizontal components of the records. Iterative forward modeling of theoretical source spectrum gives the average estimate of seismic moment (M _o), source radius (r _o) and stress drop (Δσ) as (3.2 ± 0.8) × 10²⁶ dyne cm, 13.3 ± 0.8 km and 59.2 ± 8.8 bars, respectively, for the Sikkim earthquake of September 18, 2011.     

CO<Subscript>2</Subscript> content of andesitic melts at graphite-saturated upper mantle conditions with implications for redox state of oceanic basalt source regions and remobilization of reduced carbon from subducted eclogite

We have performed experiments to determine the effects of pressure, temperature and oxygen fugacity on the CO₂ contents in nominally anhydrous andesitic melts at graphite saturation. The andesite composition was specifically chosen to match a low-degree partial melt composition that is generated from MORB-like eclogite in the convective, oceanic upper mantle. Experiments were performed at 1–3 GPa, 1375–1550?°C, and fO₂ of FMQ ?3.2 to FMQ ?2.3 and the resulting experimental glasses were analyzed for CO₂ and H₂O contents using FTIR and SIMS. Experimental results were used to develop a thermodynamic model to predict CO₂ content of nominally anhydrous andesitic melts at graphite saturation. Fitting of experimental data returned thermodynamic parameters for dissolution of CO₂ as molecular CO₂: ln(K ⁰) = ?21.79?±?0.04, ΔV ⁰?=?32.91?±?0.65 cm³mol^?1, ΔH ⁰?=?107?±?21 kJ mol^?1, and dissolution of CO₂ as CO₃ ^2?: ln(K ⁰ ) = ?21.38?±?0.08, ΔV ⁰?=?30.66?±?1.33 cm³ mol^?1, ΔH ⁰?=?42?±?37 kJ mol^?1, where K ⁰ is the equilibrium constant at some reference pressure and temperature, ΔV ⁰ is the volume change of reaction, and ΔH ⁰ is the enthalpy change of reaction. The thermodynamic model was used along with trace element partition coefficients to calculate the CO₂ contents and CO₂/Nb ratios resulting from the mixing of a depleted MORB and the partial melt of a graphite-saturated eclogite. Comparison with natural MORB and OIB data suggests that the CO₂ contents and CO₂/Nb ratios of CO₂-enriched oceanic basalts cannot be produced by mixing with partial melts of graphite-saturated eclogite. Instead, they must be produced by melting of a source containing carbonate. This result places a lower bound on the oxygen fugacity for the source region of these CO₂-enriched basalts, and suggests that fO₂ measurements made on cratonic xenoliths may not be applicable to the convecting upper mantle. CO₂-depleted basalts, on the other hand, are consistent with mixing between depleted MORB and partial melts of a graphite-saturated eclogite. Furthermore, calculations suggest that eclogite can remain saturated in graphite in the convecting upper mantle, acting as a reservoir for C.     

Coda <Emphasis Type="Italic">Q</Emphasis> in Qeshm Island,south of Iran,using aftershocks of the Qeshm earthquake of November 27, 2005

Coda wave attenuation is estimated for Qeshm Island which is located in the southeastern part of Zagros. For this purpose, the aftershocks of Qeshm earthquake in November 27, 2005, recorded within an epicentral distance less than 100 km, have been used. More than 829 earthquakes were recorded by a local temporary network consisting of 16 short period stations installed after a week after the main shock for ～10 weeks. The coda quality factor, Q _c, was estimated using the single-backscattering model in frequency bands of 0.5–24 Hz. In this research, lateral and vertical variations of coda Q in Qeshm Island are explored. In Qeshm Island, absence of significant lateral variation of coda Q is observed. To investigate the attenuation variation with depth, the coda Q value was calculated for coda time windows with different lengths (5, 10, 15, 20, 25, and 30 s). It is observed that coda Q increases with depth. However, in our study area, the rate of increase of coda Q with depth is not uniform. Beneath Qeshm Island, the rate of increase of coda Q is greater at depths less than ～40 km compared with those of larger depths. This is indicating the existence of a low attenuation anomalous structure under the ～40-km depth which may be correlated with the Moho depth in this region. The average frequency relation for this region is Q _c = 36 ± 1.2f ^{0.94 ± 0.039} at a 5 s-lapse time window length and Q _c = 110 ± 1.8f ^{0.88 ± 0.09} at a 30-s lapse time window length.     

Study of the 26 December 2011 Aswan earthquake,Aswan area,South of Egypt

The source process and parameters for a moderate earthquake of magnitude Ml 4.1 that occurred on the Kalabsha fault at the Aswan area are analyzed. The derived focal mechanisms of this event and other two aftershocks using polarities of P, SV, and SH waves show strike-slip fault with minor vertical movement of normal type. The solutions give two nodal planes trending ENE–WSW and NNW–SSE in close agreement with the surface traces of the faults crossing the area. The movement is right lateral along the first plane while left lateral along the second one. The rupture process characterization of this event has been investigated by using the empirical Green’s function deconvolution method. By inversion only for the P wave part of the records of these three events (main and other two aftershocks), the source time function for the master events and the azimuthally variations in the (RSTF) pulse amplitude are retrieved for estimating the rupture directivities. The estimated rupture direction is combined with the P-wave focal mechanisms for the three events to identify the fault plane solution for these earthquakes. Based on the width, amplitudes, and numbers of the isolated source time functions, a complex bi-lateral rupture of the studied earthquake is delineated. The source parameters of the master event is calculated and the derived corner frequencies f _o for P-wave spectra show a value of 6.6 Hz; the seismic moment (M _o ) is 4.2?×?10²² Nm; the average displacement (U) is 0.5 m; fault radius (r) 40 m; the average value of the stress drops (Δσ) is 0.6 Mpa, and the moment magnitude (M _w ) is 4.4.     

Spectral Irradiance, Phytoplankton Community Composition and Primary Productivity in a Salt Marsh Estuary, North Inlet, South Carolina, USA

We investigated spatial and temporal changes in spectral irradiance, phytoplankton community composition, and primary productivity in North Inlet Estuary, South Carolina, USA. High concentrations of colored dissolved organic matter (CDOM) were responsible for up to 84 % of the attenuation of photosynthetically available radiation (PAR). Green-yellow wavelengths were the predominant colors of light available at the two sampling sites: Clam Bank Creek and Oyster Landing. Vertical attenuation coefficients of PAR were 0.7–2.1 m^?1 with corresponding euphotic zone depths of 1.5–6.7 m. Phytoplankton biomass (as chlorophyll a [chl a]) varied seasonally with a summer maximum of 16 μg chl a l^?1 and a winter minimum of 1.4 μg chl a l^?1. The phytoplankton community consisted mainly of diatoms, prasinophytes, cryptophytes and haptophytes, with diatoms and prasinophytes accounting for up to 67 % of total chl a. Changes in phytoplankton community composition showed strongest correlations with temperature. Light-saturated chl a-specific rates of photosynthesis and daily primary productivity varied with season and ranged from 1.6 to 14 mg C (mg chl a) ^?1?h^?1 (32–803 mg C m^?3?day^?1). Calculated daily rates added up to an annual carbon fixation rate of 84 g C m^?3?year^?1. Overall, changes in phytoplankton community composition and primary productivity in North Inlet showed a strong dependence on temperature, with PAR and spectral irradiance playing a relatively minor role due to short residence times, strong tidal forcing and vertical mixing.     

Structural variation of babingtonite depending on cation distribution at the octahedral sites

Babingtonite, Ca₂Fe²⁺Fe³⁺[Si₅O₁₄(OH)] (Z?=?2, space group $ P\overline{1} $ ) from Yakuki mine (Japan), Grönsjöberget (Sweden), Kandivali Quarry (India), Baveno Quarry (Italy), Bråstad Mine (Norway), and Kouragahana (Japan), and manganbabingtonite, Ca₂(Mn²⁺, Fe²⁺)Fe³⁺[Si₅O₁₄(OH)], from Iron Cap mine (USA) were studied using electron-microprobe analysis (EMPA), ⁵⁷Fe Mössbauer analysis and single-crystal X-ray diffraction methods to determine the cation distribution at M1 and M2 and to analyze its effect on the crystal structure of babingtonite. Although all studied babingtonite crystals are relatively homogeneous, chemical zonation due to mainly Fe ? Mn substitution is observed in manganbabingtonite. Mössbauer spectra consist of two doublets with isomer shift (I.S.)?=?1.16–1.22 mm/s and quadrupole splitting (Q.S.)?=?2.33–2.50 mm/s and with I.S.?=?0.38–0.42 mm/s and Q.S.?=?0.82–0.90 mm/s, assigned to Fe²⁺ and Fe³⁺ at the M1 and M2 octahedral sites, respectively. The determined ratio of Fe²⁺/total Fe in manganbabingtonite (0.26) was smaller than that in the others (0.35–0.44) because of high Mn²⁺ content instead of Fe²⁺. The unit-cell parameters of babingtonite are a?=?7.466–7.478, b?=?11.624–11.642, c?=?6.681–6.690 Å, α?=?91.53–91.59, β?=?93.86–93.94, γ?=?104.20–104.34º, and V?=?560.2–562.3 Å³, and those of manganbabingtonite are a?=?7.4967(3), b?=?11.6632(4), c?=?6.7014(2) Å, α?=?91.602(2), β?=?93.989(2), γ?=?104.574(3)º, and V =565.09(5) ų. Structural refinements converged to R ₁ values of 1.64–3.16 %. The <M1-O> distance was lengthened due to the substitution of large octahedral cations such as Mn²⁺ for Fe²⁺. The increase of the M1-O8, M1-O8’ and M1-O13 lengths with mean ionic radii is slightly more pronounced than of the other M1-Oi lengths. The lengthened M1-O13 distance leads the positive correlation between Si5-O15-Si1 angle and M1-O13 distance. The increase of Si2-O3-Si1 and Si5-O12-Si4 angles due to the increase of mean ionic radius of M2 is also observed.     

Assessment of soil losses and siltation of the K’sob hydrological system (semiarid area—East Algeria)

Soil losses and siltation of the hydrological system (watershed–dam) of K’sob were obtained using direct and indirect methods. The Wadi K’sob watershed of 1,484 km², average slope of 0.14, and average elevation of 1,060 m is located in a semiarid climate. The average annual rainfall is 341 mm and the mean annual water discharge is 0.89 m³/s. Data from the Medjez gauging station located 6 km upstream of the dam, are the daily liquid flow and instantaneous concentrations of suspended sediments. Over a time period from 1973 to 2010, the relationship between water and sediment discharges is quantified by the equation: Q _s?=?5.6 Q ^1.31. Thus, in view of the availability data on a daily scale, the assessment of soil erodibility of the K’sob watershed was used to estimate specific soil losses of 203 t?km^?2?year^?1or 301,000 t eroded annually from the K’sob basin. The bathymetric measurements of the sediment volumes deposited in the K’sob dam, has quantified the annual siltation of 0.8 hm³, corresponding to an average erodibility of the K’sob watershed of 809 t?km^?2?year^?1. However, when adding the volume of sediment removed by the dredging operation and de-silting by the valves during heavy floods, the value of soil losses is 2,780 t?km^?2?year^?1. The indirect assessment of soil erodibility of the basin was obtained by applying two models: the quantitative geomorphological analysis (QGA) and PISA model (prediction of silting in the artificial reservoirs, in Italian: Previsioni dell’Interimento nei Serbatoi Artificiali) using physical and climatic factors in the watershed. The obtained results by QGA method underestimate specific soil losses of 524 t?km^?2?year^?1. The PISA model gives a value of 2,915 t?km^?2?year^?1, which is close to the value obtained by bathymetric measurements. This study concludes that PISA model is most suitable to estimate soil loss and siltation of the K’sob hydrological system.     

Seismicity and source parameters of local earthquakes in Bilaspur region of Himachal Lesser Himalaya

The pattern of local seismicity (110 events) and the source parameters of 26 local events (1.0?≤?Mw?≤?2.5) that occurred during May 2008 to April 2009 in Bilaspur region of Himachal Lesser Himalaya were determined. The digital records available from one station have been used to compute the source parameters and f _max based on the Brune source model (1970) and a high-frequency diminution factor (Boore 1983) above f _max. The epicentral distribution of events within 30 km of local network is broadly divided into three clusters of seismic activity: (1) a cluster located to the south of the Jamthal (JAMT) station and falls to the north of the Main Boundary Thrust (MBT) which seems to reflect the contemporary local seismicity of the segment of the MBT, (2) an elongated zone of local seismicity NE–SW trending, delineated NE of JAMT station that falls in the Lesser Himalaya between the MBT and the Main Central Thrust, and (3) NE–SW trending zone of local seismic activity located at about 10 km east of NHRI station and about 15 km northeast of NERI station and extending over a distance of about 20 km. Majority of events occur at shallow depths up to 20 km, and the maximum number of events occurs in the focal depth range between 10 and 15 km. The entire seismic activity is confined to the crust between 5 and 45 km. The average values of these source parameters range from 3.29?×?10¹⁷ to 3.73?×?10¹⁹?dyne-cm for seismic moment, 0.1 to 9.7 bars for stress drops, and 111.78 to 558.92 m for source radii. The average value of f _max for these events varies from 7 to 18 Hz and seems to be source dependent.     

P–V–T equation of state of Mn3Al2Si3O12 spessartine garnet

The thermoelastic parameters of synthetic Mn₃Al₂Si₃O₁₂ spessartine garnet were examined in situ at high pressure up to 13 GPa and high temperature up to 1,100 K, by synchrotron radiation energy dispersive X-ray diffraction within a DIA-type multi-anvil press apparatus. The analysis of room temperature data yielded K ₀ = 172 ± 4 GPa and K ₀ ^′  = 5.0 ± 0.9 when V _0,300 is fixed to 1,564.96 Å³. Fitting of P–V–T data by means of the high-temperature third-order Birch–Murnaghan EoS gives the thermoelastic parameters: K ₀ = 171 ± 4 GPa, K ₀ ^′  = 5.3 ± 0.8, (?K _0,T /?T) _P  = ?0.049 ± 0.007 GPa K^?1, a ₀ = 1.59 ± 0.33 × 10^?5 K^?1 and b ₀ = 2.91 ± 0.69 × 10^?8 K^?2 (e.g., α _0,300 = 2.46 ± 0.54 × 10^?5 K^?1). Comparison with thermoelastic properties of other garnet end-members indicated that the compression mechanism of spessartine might be the same as almandine and pyrope but differs from that of grossular. On the other hand, at high temperature, spessartine softens substantially faster than pyrope and grossular. Such softening, which is also reported for almandine, emphasize the importance of the cation in the dodecahedral site on the thermoelastic properties of aluminosilicate garnet.     

Thermo-elastic behaviour of Be2BO3OH (hambergite) up to 7 GPa and 1,100 K

The thermo-elastic behaviour of Be₂BO₃(OH)_0.96F_0.04 (i.e. natural hambergite, Z = 8, a = 9.7564(1), b = 12.1980(2), c = 4.4300(1) Å, V = 527.21(1) ų, space group Pbca) has been investigated up to 7 GPa (at 298 K) and up to 1,100 K (at 0.0001 GPa) by means of in situ single-crystal X-ray diffraction and synchrotron powder diffraction, respectively. No phase transition or anomalous elastic behaviour has been observed within the pressure range investigated. P?V data fitted to a third-order Birch–Murnaghan equation of state give: V ₀ = 528.89(4) ų, K _T0 = 67.0(4) GPa and K′ = 5.4(1). The evolution of the lattice parameters with pressure is significantly anisotropic, being: K _T0(a):K _T0(b):K _T0(c) = 1:1.13:3.67. The high-temperature experiment shows evidence of structure breakdown at T > 973 K, with a significant increase in the full-width-at-half-maximum of all the Bragg peaks and an anomalous increase in the background of the diffraction pattern. The diffraction pattern was indexable up to 1,098 K. No new crystalline phase was observed up to 1,270 K. The diffraction data collected at room-T after the high-temperature experiment showed that the crystallinity was irreversibly compromised. The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α ₀ + α ₁ T ^?1/2. The refined parameters for Be₂BO₃(OH)_0.96F_0.04 are: α ₀ = 7.1(1) × 10^?5 K^?1 and α ₁ = ?8.9(2) × 10^?4 K ^?1/2 for the unit-cell volume, α ₀(a) = 1.52(9) × 10^?5 K^?1 and α ₁(a) = ?1.4(2) × 10^?4 K ^?1/2 for the a-axis, α ₀(b) = 4.4(1) × 10^?5 K^?1 and α ₁(b) = ?5.9(3) × 10^?4 K ^?1/2 for the b-axis, α ₀(c) = 1.07(8) × 10^?5 K^?1 and α ₁(c) = ?1.5(2) × 10^?4 K ^?1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α ₀(a):α ₀(b):α ₀(c) = 1.42:4.11:1. The main deformation mechanisms in response to the applied temperature, based on Rietveld structure refinement, are discussed.     

Textural and micro-analytical insights into mafic–felsic interactions during the Oruanui eruption,Taupo

Quenched juvenile mafic inclusions (enclaves) are an occasional but informative component in the deposits of large felsic eruptions. Typically, the groundmasses of these inclusions rapidly crystallize as the mafic magma is chilled against a more voluminous, cooler felsic host, providing a physical and chemical record of the nature and timing of mafic–felsic interactions. We examine mafic inclusions of two compositional lineages (tholeiitic and calc-alkaline) from deposits of the 25.4 ka Oruanui eruption (Taupo, New Zealand). 2-D quantitative textural data from analysis of back-scattered electron images reveal a marked diversity in the groundmass textures of the inclusions, including median crystal sizes (amphibole: 14–45 µm; plagioclase: 21–75 µm) and aspect ratios (amphibole: 1.7–4.2; plagioclase: 2.1–4.0), area number densities (amphibole: 122–2660 mm^?2; plagioclase: 117–2990 mm^?2), area fractions (?) of minerals (?_plag?=?23–45%, ?_amph?=?0–28%, ?_cpx?=?0–6%, ?_oxides?=?0.6–5.5%), and the relative abundance of plagioclase and amphibole (?_plag/?_amph?=?1.0–4.6). Textural parameters vary more significantly within, rather than between, the two compositional lineages, and in some cases show marked variations across individual clasts, implying that each inclusion’s cooling history, rather than bulk composition, was the dominant control on textural development. Groundmass mineral compositions are also diverse both within and between inclusions (e.g. plagioclase from An_34–92, with typical intra-clast variability of ~?20 mol%), and do not correlate with bulk chemistry. Diverse groundmass textures and mineral and glass chemistries are inferred to reflect complex interplay of a range of factors including the degree and rate of undercooling, bulk composition, water content and, possibly, intensive variables. Our data are inconsistent with breakup of a crystallizing ponded mafic layer at the base of the Oruanui melt-dominant body, instead implying that each inclusion partially crystallized as a discrete body with a unique cooling history. Extensive ingestion of mush-derived macro-crystals suggests that mechanical breakup of mafic feeder dikes occurred within a transition zone between the mush and melt-dominant magma body. In this zone, the mush lacked yield strength, as has been inferred from field studies of narrow (meters to few tens of meters) mush-melt transition zones preserved in composite intrusions. Evidence for plastic deformation of inclusions during eruption and the abundance of fresh residual glass in inclusions from all eruptive phases suggest that the inclusions formed syn-eruptively, and must have been formed recurrently at multiple stages throughout the eruption.     

Some source and medium properties of the Vrancea seismic region, Romania

Simple spectral theory of seismic sources was used to determine source parameters directly related to medium properties (stress drop, seismic efficiency and fracture energy) and quality factors of the Vrancea (Romania) seismic region. The results show an increase in maximum static stress drop, maximum seismic efficiency and fracture energy with depth. The seismic efficiency is magnitude independent, but the stress drop is magnitude independent only for events with M_L > 3.8; below this value, the logarithm of stress drop increases quasi-linearly with magnitude. In the depth interval 50–160 km the stress drop increases with a slope of about 2–3 bar/km. The fracture energy per unit area of the fault has values of the order of 10⁵–10⁸ erg/cm².The frequency independent quality factors indicate that the attenuation of P waves is generally higher than that of S waves and that Q_p values are in agreement with recent tectonic models for the Vrancea region: total decoupling of the slab now sinking gravitationally is present only in the southwestern part of the Vrancea region, as suggested by the spatial position of intermediate depth hypocenters.     

An 8-year record of gas geochemistry and isotopic composition of methane during baseline sampling at a groundwater observation well in Alberta (Canada)

Variability in baseline groundwater methane concentrations and isotopic compositions was assessed while comparing free and dissolved gas sampling approaches for a groundwater monitoring well in Alberta (Canada) over an 8-year period. Methane concentrations in dissolved gas samples (n?=?12) were on average 4,380?±?2,452 μg/L, yielding a coefficient of variation (CV) >50 %. Methane concentrations in free gas samples (n?=?12) were on average 228,756?±?62,498 ppm by volume, yielding a CV of 27 %. Quantification of combined sampling, sample handling and analytical uncertainties was assessed via triplicate sampling (CV of 19 % and 12 % for free gas and dissolved gas methane concentrations, respectively). Free and dissolved gas samples yielded comparable methane concentration patterns and there was evidence that sampling operations and pumping rates had a marked influence on the obtained methane concentrations in free gas. δ¹³C_CH4 and δ²H_CH4 values of methane were essentially constant (?78.6?±?1.3 and ?300?±?3?‰, respectively) throughout the observation period, suggesting that methane was derived from the same biogenic source irrespective of methane concentration variations. The isotopic composition of methane constitutes a robust and highly valuable baseline parameter and increasing δ¹³C_CH4 and δ²H_CH4 values during repeat sampling may indicate influx of thermogenic methane. Careful sampling and analytical procedures with identical and repeatable approaches are required in baseline-monitoring programs to generate methane concentration and isotope data for groundwater that can be reliably compared to repeat measurements once potential impact from oil and gas development, for example, may occur.     

Polarized electronic absorption spectra of colourless chalcocyanite,CuSO4, with a survey on crystal fields in Cu2+ minerals

Polarized electronic absorption spectra of colourless chalcocyanite, CuSO₄, have been measured using microscope-spectrometric techniques. The spectra are characterized by a structured and clearly polarized band system in the near-infrared spectral range with components centred at 11,720, 10,545, 9,100, and 7,320 cm^?1, which have been assigned to crystal field d–d transitions of Cu²⁺ cations in pseudo-tetragonally elongated CuO₆ polyhedra with point symmetry C _i ( \(\bar{1}\) ). The polarization behaviour is interpreted based on a D ₂(C ₂″) pseudo-symmetry. Crystal field calculations were performed for the actual triclinic point symmetry by applying the Superposition Model of crystal fields, as well as in terms of a ‘classic’ pseudo-tetragonal crystal field approach yielding the parameters Dq _(eq) = 910, Dt = 395, and Ds = 1,336 cm^?1, corresponding to a cubically averaged Dq _cub = 679 cm^?1. A comparative survey on crystal fields in Cu²⁺ minerals shows that the low overall crystal field strength in chalcocyanite, combined with a comparatively weak pseudo-tetragonal splitting of energy levels, is responsible for its unique colourless appearance among oxygen-based Cu²⁺ minerals. The weak crystal field in CuSO₄ can be related to the lower position of the SO₄ ^2? anion compared to, e.g. the H₂O molecule in the spectrochemical series of ligands.     

Equations of state of magnesium silicates anhydrous B and superhydrous B

High-pressure single crystal X-ray diffraction experiments of phase anhydrous B and superhydrous B have been carried out to 7.3 and 7.7?GPa, respectively, at room temperature. Fitting a third-order Birch-Murnaghan equation of state to the P-V data yields values of V ₀?=?838.86?±?0.04?ų, K_T,0?=?151.5?±?0.9?GPa and K′?=?5.5?±?0.3 for Anhy-B and V ₀?=?624.71?± 0.03?ų, K_T,0?=?142.6?±?0.8?GPa and K′?=?5.8?±?0.2 for Shy-B. A similar analysis of the axial compressibilities in Anhy-B reveals that the c-axis is most compressible (K_c?=?137?±?3?GPa), the b-axis is least compressible (K_b?=?175?±?4?GPa), and the a-axis is intermediate (K_a?=?148?±?1?GPa). In Shy-B, the a-axis is most compressible (K_a?=?135?±?1?GPa), followed by the b- and c-axes which have similar compressibilities (K_b?=?146?±?3?GPa; K_c?=?148?±?3?GPa). The fact that the b-axis of Shy-B is approximately 16% more compressible than Anhy-B is primarily due to differences in the O-T layer in which the H atoms are located and the linkages with the adjacent O layers. The rigid edge-sharing chains of MgO₆ and SiO₆ octahedra in the O layer control compressibility along the a- and c-axes in both structures. The net result is a reduction in the overall anisotropic compression from ～22% in Anhy-B to ～9% in Shy-B.