The R{\overline{3}} c \to R{\overline{3}} m transition in nitratine, NaNO3, and implications for calcite, CaCO3

The temperature dependences of the crystal structure and superstructure intensities in sodium nitrate, mineral name nitratine, NaNO₃, were studied using Rietveld structure refinements based on synchrotron powder X-ray diffraction. Nitratine transforms from $R{\overline{3}} c\;\hbox{to}\;R{\overline{3}} m$ at T _c = 552(1) K. A NO₃ group occupies, statistically, two positions with equal frequency in the disordered $R{\overline{3}} m$ phase, but with unequal frequency in the partially ordered $R{\overline{3}} c$ phase. One position for the NO₃ group is rotated by 60° or 180° with respect to the other. The occupancy of the two orientations in the $R{\overline{3}} c$ phase is obtained from the occupancy factor, x, for the O1 site and gives rise to the order parameter, S = 2x ? 1, where S is 0 at T _c and 1 at 0 K. The NO₃ groups rotate in a rapid process from about 541 to T _c, where the a axis contracts. Using a modified Bragg–Williams model, a good fit was obtained for the normalized intensities (that is, normalized, NI^1/2) for the (113) and (211) reflections in $R{\overline{3}} c\hbox {\,NaNO}_{3},$ and indicates a second-order transition. Using the same model, a reasonable fit was obtained for the order parameter, S, and also supports a second-order transition.     

The calcium oxotellurate(IV) nitrates Ca5Te4O12(NO3)2(H2O)2 and Ca6Te5O15(NO3)2

Single crystals of two novel calcium oxotellurate(IV) nitrates were grown under hydrothermal conditions and were structurally characterized by X-ray diffraction. Ca $_5$ Te $_4\text {O}_{12}$ (NO $_3$ ) $_2$ (H $_2$ O) $_2$ [ $Cc$ , $Z=4$ , $a=25.258(3)$ Å, $b=5.7289(7)$ Å, $c=17.0066(19)$ Å, $\beta =124.377(2)^{\circ}$ , $R[F^2 > 2\sigma (F^2)]=0.043$ , 4083 $F^2$ data, 281 parameters] can be described as a non-classic order/disorder (OD) structure, which fulfills the basic principle of OD theory, viz. local equivalence of polytypes, but does not strictly follow the vicinity condition (VC) of OD theory. The structure is made up from an alternating stacking of non-polar layers composed of isolated [TeO $_3$ ] units and Ca $^{2+}$ ions and polar layers containing NO $_3^-$ ions and water molecules. The electron lone-pairs of the [TeO $_3$ ] units protrude into the free space of the anion/water layers. The crystal under investigation was a non-classic OD-twin of domains of a maximum degree of order (MDO). At the twin plane a fragment of the second MDO polytype is located. The main building blocks of Ca $_6$ Te $_5\text {O}_{15}$ (NO $_3$ ) $_2$ [ $P2_1/c$ , $Z=4$ , $a=15.494(2)$ Å, $b=5.6145(7)$ Å, $c=39.338(4)$ Å, $\beta =142.480(5)^{\circ}$ , $R[F^2 > 2\sigma (F^2)]=0.043$ , 3026 $F^2$ data, 307 parameters] are isolated [TeO $_3$ ] units and Ca $^{2+}$ ions which are connected to a three-dimensional framework perforated by channels in which the N atoms of the nitrate anions are located and the electron lone-pairs of the [TeO $_3$ ] units protrude. The structure can topologically be derived from the structure of Ca $_5$ Te $_4\text {O}_{12}$ (NO $_3$ ) $_2$ (H $_2$ O) $_2$ by removing the water molecules and connecting the CaTeO $_3$ layers with additional [TeO $_3$ ] units and Ca $^{2+}$ ions.     

The P–V–T equation of state of CaPtO3 post-perovskite

Orthorhombic post-perovskite CaPtO₃ is isostructural with post-perovskite MgSiO₃, a deep-Earth phase stable only above 100 GPa. Energy-dispersive X-ray diffraction data (to 9.4 GPa and 1,024 K) for CaPtO₃ have been combined with published isothermal and isobaric measurements to determine its P–V–T equation of state (EoS). A third-order Birch–Murnaghan EoS was used, with the volumetric thermal expansion coefficient (at atmospheric pressure) represented by α(T) = α₀ + α₁(T). The fitted parameters had values: isothermal incompressibility, $ K_{{T_{0} }} $  = 168.4(3) GPa; $ K_{{T_{0} }}^{\prime } $  = 4.48(3) (both at 298 K); $ \partial K_{{T_{0} }} /\partial T $  = ?0.032(3) GPa K^?1; α₀ = 2.32(2) × 10^?5 K^?1; α₁ = 5.7(4) × 10^?9 K^?2. The volumetric isothermal Anderson–Grüneisen parameter, δ _T , is 7.6(7) at 298 K. $ \partial K_{{T_{0} }} /\partial T $ for CaPtO₃ is similar to that recently reported for CaIrO₃, differing significantly from values found at high pressure for MgSiO₃ post-perovskite (?0.0085(11) to ?0.024 GPa K^?1). We also report axial P–V–T EoS of similar form, the first for any post-perovskite. Fitted to the cubes of the axes, these gave $ \partial K_{{aT_{0} }} /\partial T $  = ?0.038(4) GPa K^?1; $ \partial K_{{bT_{0} }} /\partial T $  = ?0.021(2) GPa K^?1; $ \partial K_{{cT_{0} }} /\partial T $  = ?0.026(5) GPa K^?1, with δ _T  = 8.9(9), 7.4(7) and 4.6(9) for a, b and c, respectively. Although $ K_{{T_{0} }} $ is lowest for the b-axis, its incompressibility is the least temperature dependent.     

Petrological cannibalism: the chemical and textural consequences of incremental magma body growth

The textures of minerals in volcanic and plutonic rocks testify to a complexity of processes in their formation that is at odds with simple geochemical models of igneous differentiation. Zoning in plagioclase feldspar is a case in point. Very slow diffusion of the major components in plagioclase means that textural evidence for complex magmatic evolution is preserved, almost without modification. Consequently, plagioclase affords considerable insight into the processes by which magmas accumulate in the crust prior to their eventual eruption or solidification. Here, we use the example of the 1980–1986 eruptions of Mount St. Helens to explore the causes of textural complexity in plagioclase and associated trapped melt inclusions. Textures of individual crystals are consistent with multiple heating and cooling events; changes in total pressure (P) or volatile pressure ( $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O ) are less easy to assess from textures alone. We show that by allying textural and chemical analyses of plagioclase and melt inclusions, including volatiles (H₂O, CO₂) and slow-diffusing trace elements (Sr, Ba), to published experimental studies of Mount St. Helens magmas, it is possible to disambiguate the roles of pressure and temperature to reconstruct magmatic evolutionary pathways through temperature–pressure–melt fraction (T– $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O –F) space. Our modeled crystals indicate that (1) crystallization starts at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  > 300 MPa, consistent with prior estimates from melt inclusion volatile contents, (2) crystal cores grow at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  = 200–280 MPa at F = 0.65–0.7, (3) crystals are transferred to $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  = 100–130 MPa (often accompanied by 10–20 °C of heating), where they grow albitic rims of varying thicknesses, and (4) the last stage of crystallization occurs after minor heating at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  ~ 100 MPa to produce characteristic rim compositions of An₅₀. We hypothesize that modeled $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O decreases in excess of ~50 MPa most likely represent upward transport through the magmatic system. Small variations in modeled $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O , in contrast, can be effected by fluxing the reservoir with CO₂-rich vapors that are either released from deeper in the system or transported with the recharge magma. Temperature fluctuations of 20–40 °C, on the other hand, are an inevitable consequence of incremental, or pulsed, assembly of crustal magma bodies wherein each pulse interacts with ancestral, stored magmas. We venture that this “petrological cannibalism” accounts for much of the plagioclase zoning and textural complexity seen not only at Mount St. Helens but also at arc magmas generally. More broadly we suggest that the magma reservoir below Mount St. Helens is dominated by crystal mush and fed by frequent inputs of hotter, but compositionally similar, magma, coupled with episodes of magma ascent from one storage region to another. This view both accords with other independent constraints on the subvolcanic system at Mount St. Helens and supports an emerging view of many active magmatic systems as dominantly super-solidus, rather than subliquidus, bodies.     

High-pressure elastic behavior of Ca4La6(SiO4)6(OH)2 a synthetic rare-earth silicate apatite: a powder X-ray diffraction study up to 9.33 GPa

The compression behavior of a synthetic Ca₄La₆(SiO₄)₆(OH)₂ has been investigated to about 9.33 GPa at 300 K using in situ angle-dispersive X-ray diffraction and a diamond anvil cell. No phase transition has been observed within the pressure range investigated. The values of zero-pressure volume V ₀, K ₀, and $K_{0}^{'}$ refined with a third-order Birch–Murnaghan equation of state are V ₀ = 579.2 ± 0.1 Å³, K ₀ = 89 ± 2 GPa, and $K_{0}^{'} = 10.9 \pm 0.8$ . If $K_{0}^{'}$ is fixed at 4, K ₀ is obtained as 110 ± 2 GPa. Analysis of axial compressible modulus shows that the a-axis (K _a0 = 79 ± 2 GPa) is more compressible than the c-axis (K _c0 = 121 ± 7 GPa). A comparison between the high-pressure elastic response of Ca₄La₆(SiO₄)₆(OH)₂ and the iso-structural calcium apatites is made. The possible reasons of the different elastic behavior between Ca₄La₆(SiO₄)₆(OH)₂ and calcium apatites are discussed.     

Flow rule in a high-cycle accumulation model backed by cyclic test data of 22 sands

The flow rule used in the high-cycle accumulation (HCA) model proposed by Niemunis et al. (Comput Geotech 32: 245, 2005) is examined on the basis of the data from approximately 350 drained long-term cyclic triaxial tests (N = 10⁵ cycles) performed on 22 different grain-size distribution curves of a clean quartz sand. In accordance with (Wichtmann et al. in Acta Geotechnica 1: 59, 2006), for all tested materials, the “high-cyclic flow rule (HCFR)”, i.e., the ratio of the volumetric and deviatoric strain accumulation rates \(\dot{\varepsilon}_{\rm{v}}^{{\rm acc}}/\dot{\varepsilon}_{\rm{q}}^{{\rm acc}}\) , was found dependent primarily on the average stress ratio η ^av = q ^av/p ^av and independent of amplitude, soil density and average mean pressure. The experimental HCFR can be fairly well approximated by the flow rule of the modified Cam-clay (MCC) model. Instead of the critical friction angle \(\varphi_{\rm{c}}\) which enters the flow rule for monotonic loading, the HCA model uses the MCC flow rule expression with a slightly different parameter \(\varphi_{\rm{cc}}\) . It should be determined from cyclic tests. \(\varphi_{\rm{cc}}\) and \(\varphi_{\rm{c}}\) are of similar magnitude but not always identical, because they are calibrated from different types of tests. For a simplified calibration in the absence of cyclic test data, \(\varphi_{\rm{cc}}\) may be estimated from the angle of repose \(\varphi_{\rm{r}}\) determined from a pluviated cone of sand (Wichtmann et al. in Acta Geotechnica 1: 59, 2006). However, the paper demonstrates that the MCC flow rule with \(\varphi_{\rm{r}}\) does not fit well the experimentally observed HCFR in the case of coarse or well-graded sands. For an improved simplified calibration procedure, correlations between \(\varphi_{\rm{cc}}\) and parameters of the grain-size distribution curve (d ₅₀,   C _u) have been developed on the basis of the present data set. The approximation of the experimental HCFR by the generalized flow rule equations proposed in (Wichtmann et al. in J Geotech Geoenviron Eng ASCE 136: 728, 2010), considering anisotropy, is also discussed in the paper.     

Crystallographic preferred orientation in wüstite (FeO) through the cubic-to-rhombohedral phase transition

Magnesiowüstite, (Mg_0.08Fe_0.88)O, and wüstite, Fe_0.94O, were compressed to ~36?GPa at ambient temperature in the diamond anvil cell (DAC) at the Advanced Light Source. X-ray diffraction patterns were taken in situ in radial geometry in order to study the evolution of crystallographic preferred orientation through the cubic-to-rhombohedral phase transition. Under uniaxial stress in the DAC, {100}_c planes aligned perpendicular to the compression direction. The {100}_c in cubic became { $\left\{ {10\bar 14} \right\}$ }_r in rhombohedral and remained aligned perpendicular to the compression direction. However, the {101}_c and {111}_c planes in the cubic phase split into { ${10{\bar{1}}4}$ }_r and { ${11{\bar{2}}0}$ }_r, and (0001)_r and { ${10{\bar{1}}1}$ }_r, respectively, in the rhombohedral phase. The { ${11{\bar{2}}0}$ }_r planes preferentially aligned perpendicular to the compression direction while { ${10{\bar{1}}4}$ }_r oriented at a low angle to the compression direction. Similarly, { ${10{\bar{1}}1}$ }_r showed a slight preference to align more closely perpendicular to the compression direction than (0001)_r. This variant selection may occur because the 〈 ${10{\bar{1}}4}$ 〉_r and [0001]_r directions are the softer of the two sets of directions. The rhombohedral texture distortion may also be due to subsequent deformation. Indeed, polycrystal plasticity simulations indicate that for preferred { ${10{\bar{1}}4}$ }〈 ${1{\bar{2}}10}$ 〉_r and { ${11{\bar{2}}0}$ }〈 ${{\bar{1}}101}$ 〉_r slip and slightly less active { ${10{\bar{1}}1}$ }〈 ${{\bar{1}}2{\bar{1}}0}$ 〉_r slip, the observed texture pattern can be obtained.     

The dilatant behaviour of sand–pile interface subjected to loading and stress relief

Property and behaviour of sand–pile interface are crucial to shaft resistance of piles. Dilation or contraction of the interface soil induces change in normal stress, which in turn influences the shear stress mobilised at the interface. Although previous studies have demonstrated this mechanism by laboratory tests and numerical simulations, the interface responses are not analysed systematically in terms of soil state (i.e. density and stress level). The objective of this study is to understand and quantify any increase in normal stress of different pile–soil interfaces when they are subjected to loading and stress relief. Distinct element modelling was carried out. Input parameters and modelling procedure were verified by experimental data from laboratory element tests. Parametric simulations of shearbox tests were conducted under the constant normal stiffness, constant normal load and constant volume boundary conditions. Key parameters including initial normal stress ( $ \sigma_{{{\text{n}}0}}^{\prime } $ ), initial void ratio (e ₀), normal stiffness constraining the interface and loading–unloading stress history were investigated. It is shown that mobilised stress ratio ( $ \tau /\sigma_{\text{n}}^{\prime } $ ) and normal stress increment ( $ \Updelta \sigma_{\text{n}}^{\prime } $ ) on a given interface are governed by $ \sigma_{{{\text{n}}0}}^{\prime } $ and e ₀. An increase in $ \sigma_{{{\text{n}}0}}^{\prime } $ from 100 to 400 kPa leads to a 30 % reduction in $ \Updelta \sigma_{\text{n}}^{\prime } $ . An increase in e ₀ from 0.18 to 0.30 reduces $ \Updelta \sigma_{\text{n}}^{\prime } $ by more than 90 %, and therefore, shaft resistance is much lower for piles in loose sands. A unique relationship between $ \Updelta \sigma_{\text{n}}^{\prime } $ and normal stiffness is established for different soil states. It can be applied to assess the shaft resistance of piles in soils with different densities and subjected to loading and stress relief. Fairly good agreement is obtained between the calculated shaft resistance based on the proposed relationship and the measured results in centrifuge model tests.     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

The thermodynamics of the I \overline{1} –P \overline{1} phase transition in Ca-rich plagioclase from an assessment of the spontaneous strain

In-situ powder diffraction measurements between 90 and 935?K on four anorthite-rich plagioclase samples (An₁₀₀, An₉₆Ab₄, An₈₉Ab₁₁ and An₇₈Ab₂₂) were used to determine the detailed evolution of these samples through the $I \overline{1} $ – $P \overline{1} $ phase transition. The c-type reflections indicative of $P \overline{1} $ symmetry were detected only in An₁₀₀, An₉₆Ab₄, whereas deviations in the evolution of the unit-cell parameters with temperature were observed in all samples, most prominently in the β unit-cell angle. The c-type reflections disappear at ~510 and ~425?K in An₁₀₀ and An₉₆Ab₄ respectively, and their intensity decreases according to a tricritical trend $ I^{2} \propto \left( {T - T_{\text{c}} } \right) $ . The cell parameter changes were used to determine the spontaneous strains arising from the transition which were modelled with Landau theory, allowing for low-temperature quantum saturation, in order to determine the thermodynamic behaviour. In An₁₀₀ tricritical behaviour was observed [T _c?=?512.7(4)?K; θ_s?=?394(4)] in good agreement with previous studies, and the c-type superlattice reflections indicative of $P \overline{1} $ symmetry persist up to the T _c determined from the spontaneous strain, and then disappear. The evolution of the spontaneous strain in An₉₆Ab₄ is tricritical at low temperatures [T _c?=?459(1) K, θ_s?=?396(5)] up to the temperature of disappearance of c-type reflections, but becomes second order beyond ~440?K. In An₈₉Ab₁₁ the strain displays second-order behaviour throughout [T _c?=?500(1) and θ_s?=?212(5)], and the c-type reflections are not detected in the powder diffraction patterns at any temperature. The apparent discrepancy between the absence of c-type reflections in temperature ranges where the cell parameters display significant spontaneous strain is resolved through consideration of the sizes of the anti-phase domains within the crystals. It is deduced that the tricritical phase transition occurs in well-ordered crystals with large domains in which the behavior of individual domains is dominant (i.e. in pure anorthite) or where the $P \overline{1} $ distortions within the domains are large enough to dominate the structural coherency strains between the domains. When both the magnitude of the $P \overline{1} $ pattern of displacements of the tetrahedral framework become smaller and the influence of the structural coherency between anti-phase domains becomes significant, the thermodynamic behavior becomes 2nd-order in character, the c-type reflections disappear, and the orientation of the spontaneous strain changes.     

Thermal equation of state of synthetic orthoferrosilite at lunar pressures and temperatures

Iron-rich orthopyroxene plays an important role in models of the thermal and magmatic evolution of the Moon, but its density at high pressure and high temperature is not well-constrained. We present in situ measurements of the unit-cell volume of a synthetic polycrystalline end-member orthoferrosilite (FeSiO₃, fs) at simultaneous high pressures (3.4–4.8 GPa) and high temperatures (1,148–1,448 K), to improve constraints on the density of orthopyroxene in the lunar interior. Unit-cell volumes were determined through in situ energy-dispersive synchrotron X-ray diffraction in a multi-anvil press, using MgO as a pressure marker. Our volume data were fitted to a high-temperature Birch–Murnaghan equation of state (EoS). Experimental data are reproduced accurately, with a  $\varDelta P$ Δ P  standard deviation of 0.20 GPa. The resulting thermoelastic parameters of fs are: V ₀ = 875.8 ± 1.4 Å³, K ₀ = 74.4 ± 5.3 GPa, and $\frac{{\text d}K}{{\text d}T} = -0.032 \pm 0.005\,\hbox{GPa K}^{-1}$ d K d T = - 0.032 ± 0.005 GPa K - 1 , assuming ${K}^{\prime}_{0} = 10 $ K 0 ′ = 10 . We also determined the thermal equation of state of a natural Fe-rich orthopyroxene from Hidra (Norway) to assess the effect of magnesium on the EoS of iron-rich orthopyroxene. Comparison between our two data sets and literature studies shows good agreement for room-temperature, room-pressure unit-cell volumes. Preliminary thermodynamic analyses of orthoferrosilite, FeSiO₃, and orthopyroxene solid solutions, (Mg_1?x Fe _x ) SiO₃, using vibrational models show that our volume measurements in pressure–temperature space are consistent with previous heat capacity and one-bar volume–temperature measurements. The isothermal bulk modulus at ambient conditions derived from our measurements is smaller than values presented in the literature. This new simultaneous high-pressure, high-temperature data are specifically useful for calculations of the orthopyroxene density in the Moon.     

P wave attenuation structure beneath the Kanto district and surrounding area, Japan

Spectral ratios of teleseismic P waves for 15 deep (>200 km) earthquakes recorded at 146 High-Sensitivity Seismographic Network stations in the Kanto district and its surrounding area, eastern Japan, were inverted for attenuation parameter $ t_P^{ * } $ . The dataset consisted of good-quality vertical-component seismograms, whose P phases were handpicked. The P wave spectral ratios with high signal-to-noise ratios were calculated up to 1 Hz for all the station pairs, linear regressed, and then inverted for $ t_P^{ * } $ using the technique of least squares . The result showed that the active volcanic areas were clearly characterized by high $ t_P^{ * } $ . In contrast, $ t_P^{ * } $ varied in the nonvolcanic areas. The present result on the $ t_P^{ * } $ distribution was roughly consistent with the shallow part (<30 km) of an attenuation structure, which has been previously obtained based on 3-D tomography by using records of high-frequency (around 5 Hz) P waves from local earthquakes. This suggested that the present method of $ t_P^{ * } $ estimation is valid. The advantage and possible application to other areas were also discussed.     

Crystal chemistry of Fe3+-bearing (Mg,Fe)SiO3 perovskite: a single-crystal X-ray diffraction study

Magnesium silicate perovskite is the predominant phase in the Earth’s lower mantle, and it is well known that incorporation of iron has a strong effect on its crystal structure and physical properties. To constrain the crystal chemistry of (Mg, Fe)SiO₃ perovskite more accurately, we synthesized single crystals of Mg_0.946(17)Fe_0.056(12)Si_0.997(16)O₃ perovskite at 26 GPa and 2,073 K using a multianvil press and investigated its crystal structure, oxidation state and iron-site occupancy using single-crystal X-ray diffraction and energy-domain Synchrotron Mössbauer Source spectroscopy. Single-crystal refinements indicate that all iron (Fe²⁺ and Fe³⁺) substitutes on the A-site only, where \( {\text{Fe}}^{ 3+ } /\Upsigma {\text{Fe}}\sim 20\,\% \) based on Mössbauer spectroscopy. Charge balance likely occurs through a small number of cation vacancies on either the A- or the B-site. The octahedral tilt angle (Φ) calculated for our sample from the refined atomic coordinates is 20.3°, which is 2° higher than the value calculated from the unit-cell parameters (a = 4.7877 Å, b = 4.9480 Å, c = 6.915 Å) which assumes undistorted octahedra. A compilation of all available single-crystal data (atomic coordinates) for (Mg, Fe)(Si, Al)O₃ perovskite from the literature shows a smooth increase of Φ with composition that is independent of the nature of cation substitution (e.g., \( {\text{Mg}}^{ 2+ } - {\text{Fe}}^{ 2+ } \) or \( {\text{Mg}}^{ 2+ } {\text{Si}}^{ 4+ } - {\text{Fe}}^{ 3+ } {\text{Al}}^{ 3+ } \) substitution mechanism), contrary to previous observations based on unit-cell parameter calculations.     

Impact of Data Density and Geostatistical Simulation Technique on the Estimation of Residence Times in a Synthetic Two-dimensional Aquifer

Connectivity patterns of heterogeneous porous media are important in the estimation of groundwater residence time distributions (RTDs). Understanding the connectivity patterns of a hydraulic conductivity ( \(K\) ) field often requires knowledge of the entire aquifer, which is not practical. As such, the method used to estimate unknown \(K\) values using known \(K\) values is important. This study investigates how varying levels of conditioning data and four simulation techniques, one multi-Gaussian and three multi-point, are able to recreate key \(K\) field features and connectivity patterns of a synthetic two-dimensional bimodal distributed ln( \(K\) ) field with highly connected high \(K\) features. These techniques are then assessed in the context of RTD estimation. It was found that the multi-Gaussian technique presented a bias towards earlier travel times with increased conditioning data. This was due to the inability of the method to recreate multiple scales of connecting features. Of the multi-point methods investigated, the facies method was unable to predict early arrival times. The use of a continuous variable training image produced good fits to the observed residence time distribution with a high number of conditioning points. The ability of the methods to predict the shape of residence time distributions appears to be related to their ability to reproduce the connection patterns of higher \(K\) features.     

Biradical states of oxygen-vacancy defects in ��-quartz: centers E_{2}^{\prime \prime } and E_{4}^{\prime \prime }

Several new radiation defects with total electron spin S?=?1 occurring in electron-irradiated, synthetic ??-quartz have been observed by using electron paramagnetic resonance spectroscopy. These defects are considered to be biradicals, i.e., pairs of S?=?1/2 species. The concentration of these centers depends on the condition of the fast-electron irradiation. They have different decay behaviors that allow measurements of any individual species especially when it predominates over the others. The primary spin Hamiltonian parameter matrices g ₁, g ₂, D have now been determined for two similar defects, which herein are labeled $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ . Inter-electron distances estimated by using the magnetic dipole model, suggest that the structures of centers $ E_{2}^{\prime \prime } $ and $ E_{4}^{\prime \prime } $ both involve the unpaired electrons each located in orbitals of two silicon atoms next to a common oxygen vacancy but which have slightly different Si?CSi distances at 0.90 and 0.79?nm, respectively. This model is consistent with previous DFT calculations of the triplet configurations with local energetic minima. Observed decay behaviors suggest a transformation of centers $ E_{2,4}^{\prime \prime } $ to the analogous $ E_{1}^{\prime \prime } $ center. These triplet centers in quartz provide new insights into the structures of analogous defects in amorphous silica.     

Experimental Investigation of Seepage Properties of Fractured Rocks Under Different Confining Pressures

The effectiveness of transmitting underground water in rock fractures is strongly influenced by the widths of the fractures and their interconnections. However, the geometries needed for water flow in fractured rock are also heavily controlled by the confining pressure conditions. This paper is intended to study the seepage properties of fractured rocks under different confining pressures. In order to do this, we designed and manufactured a water flow apparatus that can be connected to the electro-hydraulic servo-controlled test system MTS815.02, which provides loading and exhibits external pressures in the test. Using this apparatus, we tested fractured mudstone, limestone and sandstone specimens and obtained the relationship between seepage properties and variations in confining pressure. The calculation of the seepage properties based on the collection of water flow and confining pressure differences is specifically influenced by non-Darcy flow. The results show that: (1) The seepage properties of fractured rocks are related to confining pressure, i.e. with the increase of confining pressure, the permeability $ k $ decreases and the absolute value of non-Darcy flow coefficient $ \beta $ increases. (2) The sandstone coefficients $ k $ and $ \beta $ range from $ 1.03 \times 10^{ - 18} $ to $ 1.53 \times 10^{ - 17} $  m² and $ - 1.13 \times 10^{17} $ to $ - 2.35 \times 10^{18} $  m^?1, respectively, and exhibit a greater change compared to coefficients of mudstone and limestone. (3) From the regression analysis of experimental data, it is concluded that the polynomial function is a better fit than the power and logarithmic functions. The results obtained can provide an important reference for understanding the stability of rock surrounding roadways toward prevention of underground water gushing-out, and for developing underground resources (e.g. coal).     

Vigrishinite, Zn2Ti4 − x Si4O14(OH,H2O,□)8, a new mineral from the Lovozero alkaline complex, Kola Peninsula, Russia

A new mineral vigrishinite, epistolite-group member and first layer titanosilicate with species-defining Zn, was found at Mt. Malyi Punkaruaiv, in the Lovozero alkaline complex, Kola Peninsula, Russia. It occurs in a hydrothermally altered peralkaline pegmatite and is associated with microcline, ussingite, aegirine, analcime, gmelinite-Na, and chabazite-Ca. Vigrishinite forms rectangular or irregularly shaped lamellae up to 0.05 × 2 × 3 cm flattened on [001]. They are typically slightly split and show blocky character. The mineral is translucent to transparent and pale pink, yellowish-pinkish or colorless. The luster is vitreous. The Mohs’ hardness is 2.5–3. Vigrishinite is brittle. Cleavage is {001} perfect. D _meas = 3.03(2), D _calc = 2.97 g/cm³. The mineral is optically biaxial (?), α = 1.755(5), β = 1.82(1), γ = 1.835(8), 2V _meas = 45(10)°, 2V _calc = 50°. IR spectrum is given. The chemical composition (wt %; average of 9 point analyses, H₂O is determined by modified Penfield method) is as follows: 0.98 Na₂O, 0.30 K₂O, 0.56 CaO, 0.05 SrO, 0.44 BaO, 0.36 MgO, 2.09 MnO, 14.39 ZnO, 2.00 Fe₂O₃, 0.36 Al₂O₃, 32.29 SiO₂, 29.14 TiO₂, 2.08 ZrO₂, 7.34 Nb₂O₅, 0.46 F, 9.1 H₂O, ?0.19 O=F₂, total is 101.75. The empirical formula calculated on the basis of Si + Al = 4 is: H_7.42(Zn_1.30Na_0.23Mn_0.22Ca_0.07Mg_0.07K_0.05Ba_0.02)_Σ1.96(Ti_2.68Nb_0.41Fe _0.18 ³⁺ Zr_0.12)_Σ3.39(Si_3.95Al_0.05)_{Σ4 20.31}F_0.18. The simplified formula is: Zn₂Ti_4?x Si₄O₁₄(OH,H₂O,□)₈ (x < 1). Vigrishinite is triclinic, space group P $\bar 1$ , a = 8.743(9), b = 8.698(9), c = 11.581(11)Å, α = 91.54(8)°, β = 98.29(8)°, γ = 105.65(8)°, V = 837.2(1.5) ų, Z = 2. The strongest reflections in the X-ray powder pattern (d, Å, ?I[hkl]) are: 11.7-67[001], 8.27-50[100], 6.94-43[0 $\bar 1$ 1, $\bar 1$ 10], 5.73–54[1 $\bar 1$ 1, 002], 4.17-65[020, $\bar 1$ $\bar 1$ 2, 200], and 2.861-100[3 $\bar 1$ 0, 2 $\bar 2$ 2, 004, 1 $\bar 3$ 1]. The crystal structure model was obtained on a single crystal, R = 0.171. Vigrishinite and murmanite are close in the structure of the TiSiO motif, but strongly differ from each other in part of large cations and H-bearing groups. Vigrishinite is named in honor of Viktor G. Grishin (b. 1953), a Russian amateur mineralogist and mineral collector, to pay tribute to his contribution to the mineralogy of the Lovozero Complex. The type specimen is deposited in the Fersman Mineralogical Museum of Russian Academy of Sciences, Moscow.     

Concentration dependence of water diffusion in obsidian and dacitic melts at high-pressures

The diffusion of water in natural obsidian and model dacitic melts (Ab90Di8Wo2, mol %) has been studied at water vapor pressure up to 170 MPa, temperatures of 1200°C, H₂O contents in melts up to ～6 wt % using a high gas pressure apparatus equipped with a unique internal device. The experiments were carried out by a new low-gradient technique with application of diffusion hydration of a melt from fluid phase. The water solubility in the melts and its concentration along $C_{H_2 O} $ diffusion profiles were determined using IR microspectrometry by application of the modified Bouguer-Beer-Lambert equation. A structural-chemical model was proposed to calculate and predict the concentration dependence of molar absorption coefficients of the hydroxyl groups (OH^?) and water molecules (H₂O) in acid polymerized glasses (quenched melts) in the obsidian-dacite series. The water diffusion coefficients $D_{H_2 O} $ were obtained by the mathematical analysis of concentration profiles and the analytical solution of the second Fick diffusion law using the Boltzman-Matano method. It was shown experimentally that $D_{H_2 O} $ exponentially increases with a growth of water concentration in the obsidian and dacitic melts within the entire range of diffusion profiles. Based on the established quantitative correlation between $D_{H_2 O} $ and viscosity of such melts, a new method was developed to predict and calculate the concentration, temperature, and pressure dependences of $D_{H_2 O} $ in acid magmatic melts (obsidian, rhyolite, albite, granite, dacite) at crustal T, P parameters and water concentrations up to 6 wt %.     

Depth dependence and exponential models of permeability in alluvial-fan gravel deposits

To determine depth dependence of permeability in various geologic deposits, exponential models have often been proposed. However, spatial variability in hydraulic conductivity, K, rarely fits this trend in coarse alluvial aquifers, where complex stratigraphic sequences follow unique trends due to depositional and post-depositional processes. This paper analyzes K of alluvial-fan gravel deposits in several boreholes, and finds exponential decay in K with depth. Relatively undisturbed gravel cores obtained in the Toyohira River alluvial fan, Sapporo, Japan, are categorized by four levels of fine-sediment packing between gravel grains. Grain size is also analyzed in cores from two boreholes in the mid-fan and one in the fan-toe. Profiles of estimated conductivity, $ \overline{K} $ , are constructed from profiles of core properties through a well-defined relation between slug-test results and core properties. Errors in $ \overline{K} $ are eliminated by a moving-average method, and regression analysis provides the decay exponents of $ \overline{K} $ with depth. Moving-average results show a similar decreasing trend in only the mid-fan above ～30-m depth, and the decay exponent is estimated as ≈0.11 m^?1, which is 10- to 1,000-fold that in consolidated rocks. A longitudinal cross section is also generated by using the profiles to establish hydrogeologic boundaries in the fan.     

Experimental delineation of the “e” ⇌ I\bar 1 and “e” ⇌ C\bar 1 transformations in intermediate plagioclase feldspars

Approximately 125 hydrothermal annealing experiments have been carried out in an attempt to bracket the stability fields of different ordered structures within the plagioclase feldspar solid solution. Natural crystals were used for the experiments and were subjected to temperatures of ～650°C to ～1,000°C for times of up to 370 days at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =600 bars, or \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =1,200 bars. The structural states of both parent and product materials were characterised by electron diffraction, with special attention being paid to the nature of type e and type b reflections (at h+k=(2n+1), l=(2n+1) positions). Structural changes of the type C \(\bar 1\) → I \(\bar 1\) , C \(\bar 1\) → “e” structure, I \(\bar 1\) → “e” and “e” structure → I \(\bar 1\) have been followed. There are marked differences between the ordering behaviour of crystals with compositions on either side of the C \(\bar 1\) ? I \(\bar 1\) transition line. In the composition range ～ An₅₀ to ～ An₇₀ the e structure appears to have a true field of stability relative to I \(\bar 1\) ordering, and a transformation of the type I \(\bar 1\) ? e has been reversed. It is suggested that the e structure is the more stable ordered state at temperatures of ～ 800°C and below. For compositions more albite-rich than ～ An₅₀ the upper temperature limit for long range e ordering is lower than ～ 750°C, and there is no evidence for any I \(\bar 1\) ordering. The evidence for a true stability field for “e” plagioclase, which is also consistent with calorimetric data, necessitates reanalysis both of the ordering behaviour of plagioclase crystals in nature and of the equilibrium phase diagram for the albite-anorthite system. Igneous crystals with compositions of ～ An₆₅, for example, probably follow a sequence of structural states C \(\bar 1\) → I \(\bar 1\) → e during cooling. The peristerite, Bøggild and Huttenlocher miscibility gaps are clearly associated with breaks in the albite, e and I \(\bar 1\) ordering behaviour but their exact topologies will depend on the thermodynamic character of the order/disorder transformations.