The effect of fluid salinity on infiltration-driven contact metamorphism of carbonate rocks

Calculated phase equilibria involving minerals and H₂O–CO₂–NaCl fluid lead to predictions of how infiltration of rock by H₂O–NaCl fluids with X _NaCl in the range 0–0.3 (0–58 wt% NaCl) drives the reactions calcite + quartz = wollastonite + CO₂ and dolomite = periclase + calcite + CO₂. Calculations focus on metamorphism in four aureoles that together are representative of the normal P–T conditions and processes of infiltration-driven contact metamorphic reactions. The effect of salinity on the spatial extent of oxygen isotope alteration was also computed. The time-integrated input fluid flux (q°) that displaces the mineral reaction front an increment of distance along the flow path always increases with increasing X _NaCl. For input fluids with salinity up to approximately five times that of seawater (X _NaCl ≤ 0.05), values of q° required to explain the spatial extent of decarbonation reaction are no more than 1.1–1.5 times that computed assuming the input fluid was pure H₂O. For more saline fluids, values of q° may be up to 1.4–7.9 times that for pure H₂O. Except for reaction in the presence of halite and vapor (V), infiltration of H₂O–NaCl fluids expands the region of oxygen isotope alteration relative to the size of the region of mineral reaction. The expansion is significant only for saline fluids with X _NaCl ≥ ~0.1. Immiscible fluid phase separation and differential loss of the liquid (L) or V phase from the mineral reaction site increases the amount of reactive fluid required to advance the mineral reaction front compared to conditions under which equilibration of minerals and fluid is attained with no loss of L or V. Decarbonation reactions driven by infiltration of fluids with even modest seawater-like salinity can explain the occurrence of salt-saturated fluid and solid halide inclusions in contact metamorphosed carbonate rocks.     

Geochemical syntheses among the cratonic,off-cratonic and orogenic garnet peridotites and their tectonic implications

Garnet-bearing mantle peridotites, occurring as either xenoliths in volcanic rocks or lenses/massifs in high-pressure and ultrahigh-pressure terrenes within orogens, preserve a record of deep lithospheric mantle processes. The garnet peridotite xenoliths record chemical equilibrium conditions of garnet-bearing mineral assemblage at temperatures (T) ranging from ~700 to 1,400°C and pressures (P) > 1.6–8.9 GPa, corresponding to depths of ~52–270 km. A characteristic mineral paragenesis includes Cr-bearing pyropic garnet (64–86 mol% pyrope; 0–10 wt% Cr₂O₃), Cr-rich diopside (0.5–3.5 wt% Cr₂O₃), Al-poor orthopyroxene (0–5 wt% Al₂O₃), high-Cr spinel (Cr/(Cr + Al) × 100 atomic ratio = 2–86) and olivine (88–94 mol% forsterite). In some cases, partial melting, re-equilibration involving garnet-breakdown, deformation, and mantle metasomatism by kimberlitic and/or carbonatitic melt percolations are documented. Isotope model ages of Archean and Proterozoic are ubiquitous, but Phanerozoic model ages are less common. In contrast, the orogenic peridotites were subjected to ultrahigh-pressure (UHP) metamorphism at temperature ranging from ~700 to 950°C and pressure >3.5–5.0 GPa, corresponding to depths of >110–150 km. The petrologic comparisons between 231 garnet peridotite xenoliths and 198 orogenic garnet peridotites revealed that (1) bulk-rock REE (rare earth element) concentrations in xenoliths are relatively high, (2) clinopyroxene and garnet in orogenic garnet peridotites show a highly fractionated REE pattern and Ce-negative anomaly, respectively, (3) Fo contents of olivines for off-cratonic xenolith are in turn lower than those of orogenic garnet and cratonic xenolith but mg-number of garnet for orogenic is less than that of off-cratonic and on-cratonic xenolith, (4) Al₂O₃, Cr₂O₃, CaO and Cr# of pyroxenes and chemical compositions of whole rocks are very different between these garnet peridotites, (5) orogenic garnet peridotites are characterized by low T and high P, off-cratonic by high T and low P, and cratonic by medium T and high P and (6) garnet peridotite xenoliths are of Archean or Proterozoic origin, whereas most of orogenic garnet peridotites are of Phanerozoic origin. Taking account of tectonic settings, a new orogenic garnet peridotite exhumation model, crust-mantle material mixing process, is proposed. The composition of lithospheric mantle is additionally constrained by comparisons and compiling of the off-cratonic, on-cratonic and orogenic garnet peridotite.     

Impact of internal waves on the acoustic field at a coastal station off Paradeep,east coast of India

Internal Wave (IW) characteristics and the impact of IW on acoustic field have been studied utilizing the hourly time series of temperature and salinity data collected at a coastal site off Paradeep (north Bay of Bengal) during 24–25 October 2008. The IW characteristics, viz. period (t _per ), velocity (C _vel ), wavelength (L), and wave numbers (k), are found to be 2.133–34.72 h, 0.135 km h⁻¹, 0.37–6.2 km and 2.70–0.16 cycles km⁻¹, respectively. The semi-diurnal tidal forces are predominant than diurnal as well as at other frequencies and its contribution is about 64% towards the total potential energy (E ₀ = 3.34 J m⁻²). Sound velocity perturbations with space and time in the presence of IW field are examined from Garrettt-Munk (GM) model. Transmission loss anomaly for optimized source-receiver configuration at the depth of 53 m and range of 9 km has been computed from acoustic modelling. The loss in the acoustic transmission is found to be 38.4 dB in the presence of low-frequency IW field.     

Synthetic ANaB(Na x Li1 − x Mg1)CMg5Si8O22(OH)2 (with x = 0.6, 0.2 and 0) P21/m amphiboles at high pressure: a synchrotron infrared study

The high-pressure behavior of three synthetic amphiboles crystallized with space group P2₁/m at room conditions in the system Li₂O–Na₂O–MgO–SiO₂–H₂O has been studied by in situ synchrotron infrared absorption spectroscopy. The amphiboles have compositions ^ANa ^B(Na _x Li_1 − x Mg₁) ^CMg₅ Si₈ O₂₂(OH)₂ with x = 0.6, 0.2 and 0.0, respectively. The high-P experiments up to 32 GPa were carried out on the U2A beamline at Brookhaven National Laboratory (NY, USA) using a diamond anvil cell under non-hydrostatic or quasi-hydrostatic conditions. The two most intense absorption bands in the OH-stretching infrared spectra can be assigned to two non-equivalent O–H dipoles in the P2₁/m structure, bonded to the same local environment ^M1M3Mg₃–OH–^ANa, and pointing toward two differently kinked tetrahedral rings. In all samples these bands progressively merge to give a unique symmetrical absorption with increasing pressure, suggesting a change in symmetry from P2₁/m to C2/m. The pressure at which the transition occurs appears to be linearly correlated to the aggregate B-site dimension. The infrared spectra collected for amphibole ^B(Na_0.2Li_0.8Mg₁) in the frequency range 50 to 1,400 cm⁻¹ also show a series of changes with increasing pressure. The data reported here support the inference of Iezzi et al. (Am Miner 91:479–482, 2006a) regarding a new high-pressure amphibole polymorph.     

Alternative water resources in granitic rock: a case study from Kinmen Island,Taiwan

Kinmen Island is a small, tectonically stable, granitic island that has been suffering from a scarcity of fresh water resources due to excessive annual evapotranspiration over annual precipitation. Recent studies further indicate that shallow (0–70 m) sedimentary aquifers, the major sources of groundwater supply, have already been over-exploited. Therefore, this preliminary study is to investigate the existence of exploitable water resources that can balance the shortage of fresh water on this island. Site characterization data are obtained from island-wide geophysical surveys as well as small-scale tests performed in a study area formed by three deep (maximum depth to 560 m) vertical boreholes installed in mid-east Kinmen northeast to Taiwu Mountain. Vertical fracture frequency data indicate that the rock body is fractured with a spatially correlated pattern, from which three major fracture zones (depths 0–70, 330–360, and below 450 m) can be identified. Geologic investigations indicate that the deepest fracture zone is caused by the large-scale, steeply dipping Taiwushan fault. This fault may have caused a laterally extensive low-resistivity zone, a potential fractured aquifer, near Taiwu Mountain. The middle fracture zone is induced by the Taiwushan fault and intersects the fault approximately 21 m southeast of the study area below a depth of 350 m. Slug testing results yield fracture transmissivity varying from 4.8 × 10⁻⁷ to 2.2 × 10⁻⁴ m²/s. Cross-hole tests have confirmed that hydraulic connectivity of the deeper rock body is controlled by the Taiwushan fault and the middle fracture zone. This connectivity may extend vertically to the sedimentary aquifers through high-angle joint sets. Despite the presence of a flow barrier formed by doleritic dike at about 300 m depth, the existence of fresh as well as meteoric water in the deeper rock body manifests that certain flow paths must exist through which the deeper fractured aquifers can be connected to the upper rock body. Therefore, groundwater stored within the Taiwushan fault and the associated low-resistivity zone can be considered as additional fresh water resources for future exploitation.     

Stability and chemical equilibrium of amphibole in calc-alkaline magmas: an overview, new thermobarometric formulations and application to subduction-related volcanoes

This work focuses on a rigorous analysis of the physical–chemical, compositional and textural relationships of amphibole stability and the development of new thermobarometric formulations for amphibole-bearing calc-alkaline products of subduction-related systems. Literature experimental results (550–1,120°C, <1,200 MPa, −1 ≤ ΔNNO ≤ +5), H₂O–CO₂ solubility models, a multitude of amphibole-bearing calc-alkaline products (whole-rocks and glasses, representing 38 volcanoes worldwide), crustal and high-P (1–3 GPa) mantle amphibole compositions have been used. Calcic amphiboles of basalt-rhyolite volcanic products display tschermakitic pargasite (37%), magnesiohastingsite (32%) and magnesiohornblende (31%) compositions with aluminium number (i.e. Al# = ^[6]Al/Al_T) ≤ 0.21. A few volcanic amphiboles (~1%) show high Al# (>0.21) and are inferred to represent xenocrysts of crustal or mantle materials. Most experimental results on calc-alkaline suites have been found to be unsuitable for using in thermobarometric calibrations due to the high Al# (>0.21) of amphiboles and high Al₂O₃/SiO₂ ratios of the coexisting melts. The pre-eruptive crystallization of consistent amphiboles is confined to relatively narrow physical–chemical ranges, next to their dehydration curves. The widespread occurrence of amphiboles with dehydration (breakdown) rims made of anhydrous phases and/or glass, related to sub-volcanic processes such as magma mixing and/or slow ascent during extrusion, confirms that crystal destabilization occurs with relatively low T–P shifts. At the stability curves, the variance of the system decreases so that amphibole composition and physical–chemical conditions are strictly linked to each other. This allowed us to retrieve some empirical thermobarometric formulations which work independently with different compositional components (i.e. Si^*, Al_T, Mg^*, ^[6]Al^*) of a single phase (amphibole), and are therefore easily applicable to all types of calc-alkaline volcanic products (including hybrid andesites). The Si^*-sensitive thermometer and the fO₂–Mg^* equation account for accuracies of ±22°C (σ_est) and 0.4 log units (maximum error), respectively. The uncertainties of the Al_T-sensitive barometer increase with pressure and decrease with temperature. Near the P–T stability curve, the error is <11% whereas for crystal-rich (porphyritic index i.e. PI > 35%) and lower-T magmas, the uncertainty increases up to 24%, consistent with depth uncertainties of 0.4 km, at 90 MPa (~3.4 km), and 7.9 km, at 800 MPa (~30 km), respectively. For magnesiohornblendes, the ^[6]Al^*-sensitive hygrometer has an accuracy of 0.4 wt% (σ_est) whereas for magnesiohastingsite and tschermakitic pargasite species, H₂O_melt uncertainties can be as high as 15% relative. The thermobarometric results obtained with the application of these equations to calc-alkaline amphibole-bearing products were finally, and successfully, crosschecked on several subduction-related volcanoes, through complementary methodologies such as pre-eruptive seismicity (volcano-tectonic earthquake locations and frequency), seismic tomography, Fe–Ti oxides, amphibole–plagioclase, plagioclase–liquid equilibria thermobarometry and melt inclusion studies. A user-friendly spreadsheet (i.e. AMP-TB.xls) to calculate the physical–chemical conditions of amphibole crystallization is also provided.     

A Connectivity Index for Discrete Fracture Networks

Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model are analysed. The average number of intersections per fracture, X_f, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that X_f is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a practical example is used to demonstrate the method of computing the connectivity index.     

P–V–T equation of state of Ca3Al2Si3O12 grossular garnet

The thermoelastic parameters of synthetic Ca₃Al₂Si₃O₁₂ grossular garnet were examined in situ at high-pressure and high-temperature by energy dispersive X-ray diffraction, using a Kawai-type multi-anvil press apparatus coupled with synchrotron radiation. Measurements have been conducted at pressures up to 20 GPa and temperatures up to 1,650 K: this P, T range covered the entire high-P, T stability field of grossular garnet. The analysis of room temperature data yielded V _0,300 = 1,664 ± 2 ?³ and K ₀ = 166 ± 3 GPa for K^￠₀ K^{\prime}_{0} fixed to 4.0. Fitting of our P–V–T data by means of the high-temperature third order Birch–Murnaghan or the Mie–Grüneisen–Debye thermal equations of state, gives the thermoelastic parameters: (∂K _0,T /∂T) _P  = −0.019 ± 0.001 GPa K⁻¹ and α _0,T  = 2.62 ± 0.23 × 10⁻⁵ K⁻¹, or γ ₀ = 1.21 for fixed values q ₀ = 1.0 and θ ₀ = 823 (Isaak et al. Phys Chem Min19:106–120, 1992). From the comparison of fits from two different approaches, we propose to constrain the bulk modulus of grossular garnet and its pressure derivative to K _T0 = 166 GPa and K^￠_T0 K^{\prime}_{T0}  = 4.03–4.35. Present results are compared with previously determined thermoelastic properties of grossular-rich garnets.     

Titanium in phengite: a geobarometer for high temperature eclogites

Phengite chemistry has been investigated in experiments on a natural SiO₂–TiO₂-saturated greywacke and a natural SiO₂–TiO₂–Al₂SiO₅-saturated pelite, at 1.5–8.0 GPa and 800–1,050°C. High Ti-contents (0.3–3.7 wt %), Ti-enrichment with temperature, and a strong inverse correlation of Ti-content with pressure are the important features of both experimental series. The changes in composition with pressure result from the Tschermak substitution (Si + R²⁺ = Al^IV + Al^VI) coupled with the substitution: Al^VI + Si = Ti + Al^IV. The latter exchange is best described using the end-member Ti-phengite (KMgTi[Si₃Al]O₁₀(OH)₂, TiP). In the rutile-quartz/coesite saturated experiments, the aluminoceladonite component increases with pressure while the muscovite, paragonite and Ti-phengite components decrease. A thermodynamic model combining data obtained in this and previous experimental studies are derived to use the equilibrium MgCel + Rt = TiP + Cs/Qz as a thermobarometer in felsic and basic rocks. Phengite, rutile and quartz/coesite are common phases in HT-(U)HP metamorphic rocks, and are often preserved from regression by entrapment in zircon or garnet, thus providing an opportunity to determine the T–P conditions of crystallization of these rocks. Two applications on natural examples (Sulu belt and Kokchetav massif) are presented and discussed. This study demonstrates that Ti is a significant constituent of phengites that could have significant effects on phase relationships and melting rates with decreasing P or increasing T in the continental crust.     

Geophysical investigation for shallow subsurface geotechnical problems of Mokattam area,Cairo, Egypt

Nine vertical electrical soundings of Schlumberger configuration were measured with AB/2 = 1–500 m. Manual and computerized interpretation were done to detect the subsurface stratigraphy of the study area. The results show that the subsurface section consists of alternated units of limestone, clay, marly limestone and dolomitic limestone and the thickness of clay unit ranged from 10 to 40 m. Nine dipole–dipole sections have also been constructed to give a clearer picture of the subsurface at the study area. The length of each dipole–dipole section is 235 m, with a electrode spacing ranging between 5 and 25 m. The Res2Dinv software was used for processing and interpretation of field data. The dipole–dipole sections at the upper plateau display high resistivity values at most parts of the plateau. Twelve shallow seismic refraction profiles are measured at selected locations for the dipole sections to define the interface between the fractured limestone and the upper surface of the clay layer. Each profile consists of 24 geophones with a geophone spacing of 2–3 m. Interpretation of seismic data indicates that the surface layer of the upper plateau consists of fractured limestone with a velocity range of 1.16–1.56 km/s and another layer of compacted clay with a velocity range of 1.38–1.88 km/s. Furthermore, the surface layer of the middle plateau consists of marl and marly limestone with a velocity about 2.1 km/s and its underlying layer consists of massive limestone with a velocity of 4.94 km/s.     

Relation Between Normalized Axial Stress and Failure Strain in Heterogeneous Carbonate Rocks Exhibiting Large Axial Strains

A semi-analytical equation for the modeling of stress–strain relationship for heterogeneous carbonate rocks exhibiting large axial strains (ε_af > 1%) is formulated. The equation is derived by modifying the stress–strain model based on Haldane’s distribution proposed by Palchik (2006) for carbonate rocks exhibiting ε _af ≤ 1%. The developed exponential model is used to relate normalized axial stress (σ _a/σ _c) over the whole pre-failure strain range to current axial strain (ε _a) and failure strain (ε _af). For carbonate rocks exhibiting ε _af > 1%, the value of pre-calculated parameter δ involved in the stress–strain model is not constant, but dependent on the failure strain value (ε _af). The normalized stress–strain model can be used to calculate the failure strain in terms of uniaxial compressive strength and stress–strain measurement at one point only. The advantages of the failure strain model and ways of its use in engineering practice are discussed.     

Ab initio quantum mechanical study of γ-AlOOH boehmite: structure and vibrational spectrum

The structure and vibrational spectrum of boehmite have been investigated at the quantum-mechanical level with the CRYSTAL code, using a Gaussian-type basis set and the B3LYP Hamiltonian. Three space groups are considered in this study: Cmcm, Cmc2₁, P2₁/c. Cmcm turns out to correspond to a transition state, whereas Cmc2₁ and P2₁/c are minimum energy structures. The difference among them is the position of H atoms only, the Al-O frame being essentially the same. Harmonic frequencies at the Γ point have been computed. The comparison between calculated and experimental frequencies shows a good agreement for the Al-O part of the spectrum (under 790 cm⁻¹). For the Al-OH bending modes (800–1,300 cm⁻¹) an absolute differences of 50–100 cm⁻¹ is observed; for the OH stretching modes (3,200–3,500 cm⁻¹) it increases to 120–200 cm⁻¹: anharmonicity is large because OH groups are involved in strong hydrogen bonds.     

Vulnerability of buildings to debris flow impact

Quantitative risk assessments (QRAs) for landslide hazards are increasingly being executed to determine an unmitigated level of risk and compare it with risk tolerance criteria set by the local or federal jurisdiction. This approach allows urban planning with a scientific underpinning and provides the tools for emergency preparedness. Debris-flow QRAs require estimates of the hazard probability, spatial and temporal probability of impact (hazard assessment) and vulnerability of the elements at risk. The vulnerability term is perhaps the most difficult to estimate confidently because (a) human death in debris flows is most commonly associated with building damage or collapse and is thus an indirect consequence and (b) the type and scale of building damage is very difficult to predict. To determine building damage, an intensity index (I _DF) was created as the product of maximum expected flow depth d and the square of the maximum flow velocity v (I _DF = dv ²). The I _DF surrogates impact force and thus correlates with building damage. Four classes of building damage were considered ranging from nuisance flood/sedimentation damage to complete destruction. Sixty-six well-documented case studies in which damage, flow depth and flow velocity were recorded or could be estimated were selected through a search of the global literature, and I _DF was plotted on a log scale against the associated damage. As expected, the individual damage classes overlap but are distinctly different in their respective distributions and group centroids. To apply this vulnerability model, flow velocity and flow depth need to be estimated for a given building location and I _DF calculated. Using the existing database, a damage probability (P _DF) can then be computed. P _DF can be applied directly to estimate the likely insurance loss or associated loss of life. The model presented here should be updated with more case studies and is therefore made openly available to international researchers who can access it at .     

P–V–T equations of state of MgO and thermodynamics

A simplest equation within the framework of the Mie-Grüneisen–Einstein approach is considered. Pressure estimation values are presented that are derived by conventional arithmetic and algebraic calculations as a function of temperature and volume. The equation under consideration complies with the Mie-Grüneisen–Debye model at high temperature. Different versions of an equation of state (EoS) of MgO proposed by Speziale et al. (J Geophys Res 106B:515–528, 2001) as a pressure standard at high temperatures are subject to analyses. In the literature, at least four versions of Speziale et al. EoS of MgO are discussed; the discrepancy between them reaching a few GPa at T > 2,000 K and P > 100 GPa. Our analyses of these equations suggest that the volume dependence of the Debye temperature is accepted arbitrarily and does not agree with the definition of the Grüneisen parameter, γ = −(∂lnΘ/∂lnV) _T . Pressure as a function of temperature and volume in the Mie-Grüneisen–Einstein approach or the Gao pressure calculator can be used to estimate true pressure at compression x = V/V ₀ < 1 with the Speziale et al. EoS of MgO.     

The influence of degradation of the swamp and alpine meadows on CH<Subscript>4</Subscript> and CO<Subscript>2</Subscript> fluxes on the Qinghai-Tibetan Plateau

CH₄ and CO₂ fluxes from a high-cold swamp meadow and an alpine meadow on the Qinghai-Tibetan Plateau, subject to different degrees of degradation, were measured over a 12-month period. Air temperature, soil temperature and moisture, and the depths of the water table and thawing-freezing layer were determined. For swamp meadows, the greater the degradation, the lesser the carbon efflux. CH₄ emissions at the nondegraded swamp meadow site were 1.09–3.5 and 2.5–11.27 times greater, and CO₂ emissions 1.08–1.69 and 1.41–4.43 times greater, respectively, than those from moderately and severely degraded sites. For alpine meadows, the greater the degradation, the greater the CH₄ consumption and CO₂ emissions. CH₄ consumption at the severely degraded alpine meadow site was 6.6–21 and 1.1–5.25 times greater, and CO₂ emissions 1.05–78.5 and 1.04–6.28 times greater, respectively, than those from the nondegraded and moderately degraded sites. The CH₄ and CO₂ fluxes at both sites were significantly correlated (R ² > 0.59, P < 0.05) with air temperature, soil temperature, and topsoil (0–5 cm depth) moisture, indicating these to be the main environmental factors affecting such fluxes.     

Calcic amphiboles in calc-alkaline and alkaline magmas: thermobarometric and chemometric empirical equations valid up to 1,130°C and 2.2 GPa

The following article presents constraints of the stability of Mg-rich (Mg/(Mg + Fe²⁺) > 0.5) calcic amphibole in both calc-alkaline and alkaline magmas, testing of previous thermobarometers, and formulation of new empirical equations that take into consideration a large amount of literature data (e.g. more than one thousand amphibole compositions among experimental and natural crystals). Particular care has been taken in choosing a large number of natural amphiboles and selecting quality experimental data from literature. The final database of experimental data, composed of 61 amphiboles synthesized in the ranges of 800–1,130°C and 130–2,200 MPa, indicates that amphibole crystallization occurs in a horn-like P–T stability field limited by two increasing curves (i.e. the thermal stability and an upper limit), which should start to bend back to higher pressures. Among calcic amphiboles, magnesiohornblendes and tschermakitic pargasites are only found in equilibrium with calc-alkaline melts and crystallize at relatively shallow conditions (P up to ~1 GPa). Kaersutite and pargasite are species almost exclusively found in alkaline igneous products, while magnesiohastingsite is equally distributed in calc-alkaline and alkaline rocks. The reliability of previous amphibole applications was checked using the selected experimental database. The results of this testing indicate that none of the previous thermobarometers can be successfully used to estimate the P, T and fO₂ in a wide range of amphibole crystallization conditions. Multivariate least-square analyses of experimental amphibole compositions and physico-chemical parameters allowed us to achieve a new thermobarometric model that gives reasonably low uncertainties (T ± 23.5°C, P ± 11.5%, H₂O_melt ± 0.78wt%) for calc-alkaline and alkaline magmas in a wide range of P–T conditions (up to 1,130°C and 2,200 MPa) and ∆NNO values (±0.37 log units) up to 500 MPa. The ^AK-^[4]Al relation in amphibole can be readily used to distinguish crystals of calc-alkaline liquids from those of alkaline magmas. In addition, several chemometric equations allowing to estimate the anhydrous composition of the melts in equilibrium with amphiboles of calc-alkaline magmas were derived.     

Magnetic fabric development in a highly anisotropic magnetite-bearing ductile shear zone (Seve Nappe Complex,Scandinavian Caledonides)

Magnetite-bearing mylonitic garnet–micaschists close to the major suture between the Baltica and Iapetus terranes (Seve Nappe Complex, Scandinavian Caledonides) show very high anisotropy of magnetic susceptibility (AMS) with corrected degree of anisotropy (P′) up to 4.8. Three different magnetic fabric types can be distinguished. They correspond to protomylonite (type I, P′ < 2), mylonite (type II, 2 < P′ < 3), and ultramylonite (type III, P′ > 3), respectively. The orientation of the ellipsoid axes from all applied magnetic fabric methods in this study is similar with shallow dips of the metamorphic foliation toward WSW and subhorizontal, mostly NW–SE trending mineral lineation. Differences between subfabrics were minimized under high shear strain as all markers tend to align parallel with the shear plane. The very high anisotropies and mostly oblate ellipsoid shapes of type III correlate with high magnetic susceptibility (k _mean up to 55 × 10⁻³ SI units) and are related to the concentration of magnetite aggregates with shape-preferred orientation. They show a distinct field dependence of magnetic susceptibility of up to 10% in the k _max-direction. We attribute this field dependence to a “memory” of high strains in the domain walls of the crystals acquired during synkinematic magnetite growth during shear zone fabric development at temperatures of 550–570°C.     

Single-crystal EPR and DFT study of a <Superscript>VI</Superscript>Al–O<Superscript>−</Superscript>–<Superscript>VI</Superscript>Al center in jeremejevite: electronic structure and <Superscript>27</Superscript>Al hyperfine constants

Single-crystal electron paramagnetic resonance (EPR) spectra of a gem-quality jeremejevite, Al₆B₅O₁₅(F, OH)₃, from Cape Cross, Namibia, reveal an S = 1/2 hole center characterized by an ²⁷Al hyperfine structure arising from interaction with two equivalent Al nuclei. Spin-Hamiltonian parameters obtained from single-crystal EPR spectra at 295 K are as follows: g ₁ = 2.02899(1), g ₂ = 2.02011(2), g ₃ = 2.00595(1); A ₁/g _e β _e  = −0.881(1) mT, A ₂/g _e β _e  = −0.951(1) mT, and A ₃/g _e β _e  = −0.972(2) mT, with the orientations of the g ₃- and A ₃-axes almost coaxial and perpendicular to the Al–O–Al plane; and those of the g ₁- and A ₁-axes approximately along the Al–Al and Al–OH directions, respectively. These results suggest that this aluminum-associated hole center represents hole trapping on a hydroxyl oxygen atom linked to two equivalent octahedral Al³⁺ ions, after the removal of the proton (i.e., a ^VIAl–O⁻–^VIAl center). Periodic ab initio UHF and DFT calculations confirmed the experimental ²⁷Al hyperfine coupling constants and directions, supporting the proposed structural model. The ^VIAl–O⁻–^VIAl center in jeremejevite undergoes the onset of thermal decay at 300 °C and is completely bleached at 525 °C. These data obtained from the ^VIAl–O⁻–^VIAl center in jeremejevite provide new insights into analogous centers that have been documented in several other minerals.     

Early Cretaceous highly positive <Emphasis Type="Italic">ε</Emphasis><Subscript>Nd</Subscript> felsic volcanic rocks from the Hinggan Mountains,NE China: origin and implications for Phanerozoic crustal growth

An early Cretaceous (135 ± 2 Ma) felsic volcanic suite of dacite and rhyolite from Huolinhe, NE China is characterized by large ion lithophile element and light REE enrichment and high field strength element (HFSE, e.g., Nb and Ta) and Ti–P depletion, and bulk silicate earth-like Sr [⁸⁷Sr/⁸⁶Sr(i) = 0.70409–0.70481], quite radiogenic Nd [ε _Nd(t) = +3.98 to +5.88], Pb [e.g., ²⁰⁶Pb/²⁰⁴Pb(i) = 18.46–18.55] and Hf [ε _Hf(t) ~+9.2] isotope compositions. Compared with contemporaneous mafic rocks in the region, these felsic rocks have even higher Nd and Hf isotopic ratios, precluding an origin through differentiation of coeval mantle-derived magmas. Isotope calculation results suggest that these magmas were derived from a preexistent mixture composed of mainly juvenile crust (70–80%), and a subordinate recycled crustal component (20–30%) having highly radiogenic Sr and Pb and unradiogenic Nd and Hf. About 25–30% melting of such a mixed source produced the primary dacitic magma. The rhyolites, which have relatively low MgO, FeO*, Al₂O₃, CaO, TiO₂, P₂O₅, Na₂O, Ba, Sr, REE, HFSE and Y, were differentiates of the dacites after removal of a fractional assemblage of hornblende + plagioclase + K-feldspar + apatite + zircon. Considering the prolonged events (from 262 to 130 Ma) that produced such highly positive ε _Nd felsic igneous rocks in the region, we prefer a pre-Mesozoic crustal growth model related to arc accretion associated with the Paleo-Asian Ocean subduction.