Reaction microstructures in corundum‐ and kyanite‐bearing mafic mylonites from the Takahama metamorphic rocks,western Kyushu,Southwest Japan

High‐grade mylonites occur in the Takahama metamorphic rocks, a member of the high‐pressure low‐temperature type Nagasaki Metamorphic Rocks, western Kyushu, Japan. Mafic layers within the mylonites retain reaction microstructures consisting of margarite aggregates armoring both corundum and kyanite. The following retrograde reaction well accounts for the microstructures in the CaO–Al₂O₃–SiO₂–H₂O system: 3Al₂O₃ + 2Al₂SiO₅ + 2Ca₂Al₃Si₃O₁₂(OH) + 3H₂O = 2Ca₂Al₈Si₄O₂₀(OH)₄ (corundum + kyanite + clinozoisite + fluid = margarite). Mass balance analyses and chemical potential modeling reveal that the chemical potential gradients present between kyanite and corundum have likely driven the transport of the CaO and SiO₂ components. The mylonitization is considered to take place chronologically after peak metamorphism and before the above reaction, based on the following features: approximately constant thickness of the margarite aggregates, random orientation of margarite, and local modification of garnet composition at a boudin neck that formed during mylonitization. The estimated peak temperature of 640°C and the pressure–temperature conditions of the above reaction indicate that the mylonitization took place at temperature between 530 and 640°C at pressures higher than 1.2 GPa, approximately equivalent to the depth of the lower crust of island arcs.     

Microstructure characteristics of illite from Chuanlingou Formation of Changcheng System in Jixian County,Tianjin City

The microstructure charateristics of illite from the Chuanlinggou Formation of Changcheng System (Chch) in Jixian County, Tianjin City has been studied by means of high-resolution transmission electron microscopy (HRTEM), selected-area electron diffraction (SAED), X-ray energy dispersive spectrum (EDS) and X-ray powder diffraction (XRD). The Kübler index of "illite crystallinity" is 0.505°△2θ, which indicates that the host rock is in a middle diagenetic stage. The chemical analyses of EDS for illite studied indicate a heavily absent in interlayer cation and an average chemical formula of K0.57(Al1.80Mg0.42Fe2 0.12)∑=2.34 (Si2.92Al1.08)∑=4O10(OH)2. It is found, from one-dimensional lattice images, that the layers of illite not only stack in a flat way but also in a curving way. A "matting fabric" illite structure results from stacking faults. Combined with SAED analysis the illite studied can be affirmed as 1M illite. The two-dimensional lattice images are obtained from [100] and [110] incidences, whose lattice images have the same d-values but different intersecting angles. The two-dimensional lattice image with [100] incidence is orthogonal to each other, whilst that with [110] incidence is oblique. This paper provides some important structure information of authigenic clay minerals for the well-known mesoproterozoic section of Jixian County.     

Model of vibrational level populations of Herzberg states of oxygen molecules at heights of the lower thermosphere and mesosphere

The paper presents a model of the kinetics of electronically excited O₂(c¹Σ_u⁻,v), O₂(A′³Δ_u,v), O₂(A³Σ_u⁺,v) molecules at heights of the lower thermosphere and mesosphere with allowance for electronic excitation transfer processes during molecular collisions. The model is used to calculate the relative O₂(A³Σ_u⁺,v) and O₂(A′³Δ_u,v) populations at heights of 80–110 km. The calculated populations are compared with the available literature results on experimental estimates, and good agreement is obtained. It is shown how the increase in the quenching rates of the considered states by oxygen atoms affects the calculation results.     

Determination of thermodynamic properties of (Fe,Mg)-pyroxenes at 1000 K by the emf method

Emf measurements were made on the cell Pt|Fe,(Fe,Mg)_xSi₂O₆,SiO₂|(ZrO₂)_0.85(CaO)_0.15|Fe,FeO|Pt at 1000 K. Using the present data, the standard free energy of formation of ferrosilite (compound FeSiO₃), from the component oxides FeO and SiO₂, is calculated to be −6.35 ± 0.80 kJ/mol. The activity-composition relation for pyroxene solid solution shows that it has a positive deviation from ideality at 1000 K. The present results are compared with the results of other workers.ΔG_mix andΔG^ex are calculated and plotted againstN_FeSiO3.     

Phase relations of titan-phlogopite,K2Mg4TiAl2Si6O20(OH)4: a refractory phase in the upper mantle?

Experimental data on the stability of titan-phlogopite [K₂Mg₄TiAl₂Si₆O₂₀(OH)₄] are presented which show it to be stable to substantially higher temperatures than normal phlogopite [K₂Mg₆Al₂Si₆O₂₀(OH)₄]. A qualitative model to explain the role of titan-phlogopite during magma generation is put forward. Breakdown of titan-phlogopite during melting at depth (> 150km) on subducted lithospheric slabs is believed responsible for the concomitant increase of K and Ti observed in magmas erupted during orogenic volcanism. At lower pressures (up to about 10 kbar) beneath mid-oceanic ridges, titan-phlogopite is predicted to behave as a refractory phase during partial melting in the mantle, especially if H₂O-excess conditions pertain, although at higher pressures in this environment it would almost certainly behave as a low-melting component.     

Source parameters of bohai earthquake, July 18, 1969 and yongshan earthquake, May 11, 1974 determined by synthetic seismograms of teleseismic P waves

The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M _o=0.9×10¹⁹Nm,L=59.9km.V _R=3.5km/s, ∧ _R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M _o=1.3×10¹⁸Nm,L=48.8km,V _R=3km/s, ∧ _R =−10°, where φ is strike, δ dip angle, λ slip angle,M _o seismic moment,L rupture length,V _R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ _R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.     

Crystal structure and compositional variation of Angra dos Reis fassaite

The crystal structure of fassaite from the Angra dos Reis meteorite has been determined by least-squares refinement of three-dimensional X-ray data to anR value of 3.3%. The pyroxene is monoclinic, space groupC2/c, with unit-cell dimensionsa = 9.738(1),b = 8.874(2),c = 5.2827(5)Å, β = 105.89(1)°, andV = 439.1(1)ų. Average bond lengths are (Si,Al)-O = 1.651, M1-O = 2.061, and M2-O = 2.489Å. The distribution of iron and magnesium between M1 and M2 suggests a temperature of equilibration greater than 1000°C.Electron microprobe analysis of several fassaite grains reveals small but statistically significant variations of (Mg + Si) versus (Al - Ti). The range of fassaite composition may be represented byEn₃Hd₂₂TiCpx₆(Di_53±2CaTs_16?2) whereEn=Mg₂Si₂O₆,Hd=CaFeSi₂O₆,TiCpx=CaTiAl₂O₆,Di=CaMgSi₂O₆,CaTs=CaAl₂SiO₆. Most fassaite analyses calculated on the basis of four cations yielded greater than six anions, suggesting that part of the titanium or chromium might be reduced to Ti³⁺ or Cr²⁺.     

SSQ-1 Digital Tape Recording Horizontal Pendulum Tiltmeter

The SSQ-1 Digital Tape Recording Horizontal Pendulum Tiltmeter is an instrument of high sensitivity to detect the changes of ground tilt. It uses a quartz horizontal pendulum held by Z?llner bifilar suspension to sense vertical displacement, the eddy-current transducer on the pendulum converts its displacements into electric signals. Then a microcomputer is used for data acquisition and printing as well as digit tape recording. And at the same time, a filtering pen recorder is used for visible recording. The scale value of the instrument is calibrated by computer using the known tilt angle of the bulging plate — mercury cup. The scale value is 0.3–0.5 m(″)/mV. The following are the testing results for the earth tide observation obtained from the east-west component at Baijiatan Seismic Station, Beijing: r(O₁): 0.6490±0.0179 α(O₁): −3.83°±1.58° r(K₁): 0.8049±0.0128 α(K₁): 1.40°±0.90° r(M₂): 0.6699±0.0040 α(M₂): −0.27°±0.34° r(S₂): 0.7316±0.0075 α(S₂): 2.83°±0.58° r(M₃): 0.8497±0.0964 α(M₃): 1.61°±6.51° The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 353–359, 1993. In this work also participated Mr. Huai-Wen PAN and Wei-Jin ZHANG. And we here give special thanks to Yun-Zao XI, Qin-Wen XI, Bo-Xong TANG, Yi-Hui CHEN and An-Xu WU for their kind help in the development and test of the instrument. This paper is translated by Mr. Zhong ZHENG.     

An anomalous subauroral red arc on 4 August, 1972: comparison of ISIS-2 satellite data with numerical calculations

This study compares the Isis II satellite measurements of the electron density and temperature, the integral airglow intensity and volume emission rate at 630 nm in the SAR arc region, observed at dusk on 4 August, 1972, in the Southern Hemisphere, during the main phase of the geomagnetic storm. The model results were obtained using the time dependent one-dimensional mathematical model of the Earth’s ionosphere and plasmasphere (the IZMIRAN model). The major enhancement to the IZMIRAN model developed in this study to explain the two component 630 nm emission observed is the analytical yield spectrum approach to calculate the fluxes of precipitating electrons and the additional production rates of N⁺₂, O⁺₂, O⁺(⁴S), O⁺(²D), O⁻(²P), and O⁺(²P) ions, and O(¹D) in the SAR arc regions in the Northern and Southern Hemispheres. In order to bring the measured and modelled electron temperatures into agreement, the additional heating electron rate of 1.05 eV cm⁻³ s⁻¹ was added in the energy balance equation of electrons at altitudes above 5000 km during the main phase of the geomagnetic storm. This additional heating electron rate determines the thermally excited 630 nm emission observed. The IZMIRAN model calculates a 630 nm integral intensity above 350 km of 4.1 kR and a total 630 nm integral intensity of 8.1 kR, values which are slightly lower compared to the observed 4.7 kR and 10.6 kR. We conclude that the 630 nm emission observed can be explained considering both the soft energy electron excited component and the thermally excited component. It is found that the inclusion of N₂(v > 0) and O₂(v > 0) in the calculations of the O⁺(⁴S) loss rate improves the agreement between the calculated N_e and the data on 4 August, 1972. The N₂(v > 0) and O₂(v > 0) effects are enough to explain the electron density depression in the SAR arc F-region and above F2 peak altitude. Our calculations show that the increase in the O⁺ + N₂ rate factor due to the vibrationally excited nitrogen produces the 5–19% reductions in the calculated quiet daytime peak density and the 16–24% decrease in NmF2 in the SAR arc region. The increase in the O⁺ + N₂ loss rate due to vibrationally excited O₂ produces the 7–26% decrease in the calculated quiet daytime peak density and the 12–26% decrease in NmF2 in the SAR arc region. We evaluated the role of the electron cooling rates by low-lying electronic excitation of O₂(a¹δ_g) and O₂(b¹σ_g⁺), and rotational excitation of O₂, and found that the effect of these cooling rates on T_e can be considered negligible during the quiet and geomagnetic storm period 3–4 August, 1972. The energy exchange between electron and ion gases, the cooling rate in collisions of O(³P) with thermal electrons with excitation of O(¹D), and the electron cooling rates by vibrational excitation of O₂ and N₂ are the largest cooling rates above 200 km in the SAR arc region on 4 August, 1972. The enhanced IZMIRAN model calculates also number densities of N₂(B³π_g⁺), N₂(C³π_u), and N₂(A³σ_u⁺) at several vibrational levels, O(¹S), and the volume emission rate and integral intensity at 557.7 nm in the region between 120 and 1000 km. We found from the model that the integral integral intensity at 557.7 nm is much less than the integral intensity at 630 nm.     

Echo intensity data as a directional antenna for observing processes above sloping ocean bottoms

Relative ‘echo intensity’ data (dI) from a bottom-mounted four-beam 300 kHz acoustic Doppler current profiler (ADCP) are used to infer propagation of vigorous processes above a continental slope. The 3- to 60-m horizontal beam spread and the 2-Hz sampling allow the distinction of different arrival times t _i , i = 1,..., 4, at different distances in the acoustic beams from sharp changes in dI-content associated with frontal non-linear and turbulent bores or ‘waves’. The changes in dI are partially due to variations in amounts of resuspended material carried by the near-bottom turbulence and partially due to the fast variations in density stratification (‘stratified turbulence’), as inferred from 1-Hz sampled thermistor string data above the ADCP. Such bores are observed to pass the mooring up to 80 m above the bottom, having typical propagation speeds c = 0.15–0.5 m s⁻¹, as determined from dI(t _i ). Particle speeds in the immediate environment of a bore amount to |u|_env=c ± 0.05 m s⁻¹, the equality being a necessary condition for kinematic instability, whilst the maximum particle speeds amount |u|_max = 1.2–2c. The dI-determined directions of up-, down- and alongslope processes are all to within ±10° of the ADCP’s beam-spread averaged current (particle velocity) data.     

Pyroxene-garnet solid-solution equilibria in the systems Mg4Si4O12Mg3Al2Si3O12 and Fe4Si4O12Fe3Al2Si3O12 at high pressures and temperatures

Pyroxene-garnet solid-solution equilibria have been studied in the pressure range 41–200 kbar and over the temperature range 850–1,450°C for the system Mg₄Si₄O₁₂Mg₃Al₂Si₃O₁₂, and in the pressure range 30–105 kbar and over the temperature range 1,000–1,300°C for the system Fe₄Si₄O₁₂Fe₃Al₂Si₃O₁₂. At 1,000°C, the solid solubility of enstatite (MgSiO₃) in pyrope (Mg₃Al₂Si₃O₁₂) increases gradually to 140 kbar and then increases suddenly in the pressure range 140–175 kbar, resulting in the formation of a homogeneous garnet with composition Mg₃(Al_0.8Mg_0.6Si_0.6)Si₃O₁₂. In the MgSiO₃-rich field, the three-phase assemblage of β- or γ-Mg₂SiO₄, stishovite and a garnet solid solution is stable at pressures above 175 kbar at 1,000°C. The system Fe₄Si₄O₁₂Fe₃Al₂Si₃O₁₂ shows a similar trend of high-pressure transformations: the maximum solubility of ferrosilite (FeSiO₃) in almandine (Fe₃Al₂Si₃O₁₂) forming a homogeneous garnet solid solution is 40 mol% at 93 kbar and 1,000°C.If a pyrolite mantle is assumed, from the present results, the following transformation scheme is suggested for the pyroxene-garnet assemblage in the mantle. Pyroxenes begin to react with the already present pyrope-rich garnet at depths around 150 km. Although the pyroxene-garnet transformation is spread over more than 400 km in depth, the most effective transition to a complex garnet solid solution takes place at depths between 450 and 540 km. The complex garnet solid solution is expected to be stable at depths between 540 and 590 km. At greater depths, it will decompose to a mixture of modified spinel or spinel, stishovite and garnet solid solutions with smaller amounts of a pyroxene component in solution.     

Oxysulfates of copper, sodium, and potassium in the lava flows of the 2012–2013 Tolbachik Fissure Eruption

Samples from the surface of lava flows discharged by the 2012–2013 Tolbachik Fissure Eruption were found to contain oxysulfates of copper, sodium, and potassium: K₂Cu₃O(SO₄)₂ (fedotovite), NaKCu₂O(SO₄)₂, and Na₃K₅Cu₈O₄(SO₄)₈. The last two phases have no naturally occurring or synthetic analogues that we are aware of. They form flattened crystals of prismatic to long-prismatic habits. The crystals of Na₃K₅Cu₈O₄(SO₄)₈ have a chemical composition corresponding to the empirical formula Na_2.22K_5.47Cu_8.02S_8.05O₃₆. An X-ray analysis of this compound showed that it has a monoclinic symmetry, P2/c, a = 13.909(4), b = 4.977(1), c = 23.525(6) Å, β = 90.021(5)°, V = 1628.3(7) ų. The crystal structure was determined by direct techniques and refined to yield R ₁ for 3955 reflexes//web// with F ² > 4σF. The compound NaKCu₂O(SO₄)₂ also belongs to the monoclinic system, P2/c, a = 14.111(4), b = 4.946(1), c = 23.673(6) Å, β = 92.052(6)°, V = 1651.1(8) ų. The structure was determined by direct techniques to yield a tentative structural model that has been refined up to R ₁ = 0.135 for 4088 reflexes with F ² > 4σF. The crystal structure of Na₃K₅Cu₈O₄(SO₄)₈ is based on chains of [O₂Cu₄]⁴⁺ consisting of rib-coupled oxy-centered tetrahedrons of (OCu₄)⁶⁺. The chains are surrounded by sulfate radicals, resulting in columns of {[O₂Cu₄](SO₄)₄}^4? aligned along the b axis. The interchain space contains completely ordered positions of Na⁺ and K⁺ cations. The principle underlying the connection of NaKCu₂O(SO₄)₂ columns in the crystal structure of {[O₂Cu₄](SO₄)₄}^4? is different, in view of the relation Na:K = 1 as contrasted with 3:5 for the compound Na₃K₅Cu₈O₄(SO₄)₈. The presence of oxy-centered tetrahedrons in the structure of these new compounds furnishes an indirect hint at the importance of polynuclear copper-oxygen radicals with centering oxygen atoms as forms of transport of copper by volcanic gases.     

Internal waves in a randomly stratified fluid

Abstract

We discuss the propagation of internal waves in a rotating stratified unbounded fluid with randomly varying stability frequency, N. The first order smoothing approximation is used to derive the dispersion relation for the mean wave field when N is of the form N ² = N _o ²(1 + ?μ), where μ is a centered stationary random function of either depth (z) or time (t), N _o = constant and O < ?² ≦ 1. Expressions are then derived for the change in phase speed and growth rate due to the random fluctuations μ; in particular, attention is focused on the behaviour of these expressions for short and long correlation lengths (case μ = μ(z)) and times (case μ = μ(t)). For the case μ = μ(z), which represents a model for the temperature and salinity fine-structure in the ocean, the appropriate statistics of the fluctuations observed at station P (50°N, 145°W) have been incorporated into the theory to estimate the actual importance of the effects due to these random fluctuations. It is found that the phase speed of the mean wave decreases significantly if (i) the wavelength is short compared to g/N_o ² or (ii) the wave number vector is essentially horizontal and the wave frequency is very close to N _o. Also, the random fluctuations cause a significant growth (decay) in the amplitude of a wave propagating upwards (downwards) through a depth of a few kilometers. However, in the direction of energy propagation, the kinetic energy is conserved. Finally, it is shown that the average effect of the depth dependent fluctuations at station P is to slightly decrease the stability frequency and the magnitude of the group velocity.     

High-pressure phase relations in the system CaAl4Si2O11–NaAl3Si3O11 with implication for Na-rich CAS phase in shocked Martian meteorites

High-pressure phase relations in the system NaAl₃Si₃O₁₁–CaAl₄Si₂O₁₁ were examined at 13–23 GPa and 1600–1900 °C, using a multianvil apparatus. A Ca-aluminosilicate with CaAl₄Si₂O₁₁ composition, designated CAS phase, is stable above about 13 GPa at 1600 °C. In the system NaAl₃Si₃O₁₁–CaAl₄Si₂O₁₁, the CAS phase dissolving NaAl₃Si₃O₁₁ component coexists with jadeite, corundum and stishovite below 22 GPa, above which the CAS phase coexists with Na-rich calcium ferrite, corundum and stishovite. At 1600 °C, the solubility of NaAl₃Si₃O₁₁ component in the CAS solid solution increases with increasing pressure up to about 50 mol% at about 22 GPa, above which the solubility decreases with pressure. The maximum solubility of NaAl₃Si₃O₁₁ component in the CAS phase increases with temperature up to around 70 mol% at 1900 °C at 22 GPa. The dissociation of NaAlSi₂O₆ jadeite to NaAlSiO₄ calcium ferrite plus stishovite occurs at about 22 GPa. Lattice parameters of the CAS phase with the hexagonal Ba-ferrite structure change with increase of the NaAl₃Si₃O₁₁ component: a-axis decreases and c-axis slightly increases, resulting in decrease of molar volume. Enthalpies of the CAS solid solutions were measured by high-temperature drop-solution calorimetry techniques. The results show that enthalpy of hypothetical NaAl₃Si₃O₁₁ CAS phase is much higher than the mixture of NaAlSi₂O₆ jadeite, corundum and stishovite and is close to that of the mixture of NaAlSiO₄ calcium ferrite, corundum and stishovite. When we adopt the Na:Ca ratio of 75:25 of the natural Na-rich CAS phase in a shocked Martian meteorite, Zagami, the phase relations determined above suggest that the natural CAS phase crystallized from melt at pressure around 22 GPa and temperature close to or higher than 2000–2200 °C. The inferred P, T conditions are consistent with those estimated using other high-pressure minerals in the shocked meteorite.     

Fundamental geodetic parameters of the earth's figure and the structure of the earth's gravity field derived from satellite data

Summary Using the geocentric constant GM=398 601.3 × 10 ⁹ m ³s ^–2 , the known value of the angular velocity of the Earth's rotation , Stokes' constants J _n ^(k) and S _n ^(k) upto n=21 (zonal), n=16 (tesseral and sectorial) [2], the geocentric co-ordinates and heights above sea-level of SAO satellite stations [2], the following will be derived: the potential on the geoid W_o, the scale factor for lengths R_o=GM/W_o, the radius-vector of the surface W=W_o, the parameters of the best-fitting Earth tri-axial ellipsoid, and the components of the deflections of the vertical with respect to the geocentric rotational IAG ellipsoid (Lucerne 1967), as well as to the best-fitting geocentric tri-axial ellipsoid. Some of the differences in the structure of the gravity field over the Northern and Southern Hemispheres will be given, and the mean values of gravity over the equatorial zone, determined from the dynamics of satellite orbits, on the one hand, and from terrestrial gravity data, on the other, will be compared.Presented at the Fifteenth IUGG General Assembly, Moscow, July 30 — August 14, 1971.     

Distribution of ionospheric O<Superscript>+</Superscript> ion in synchronous altitude region

Based on satellite observation data, using dynamics equation, the ionospheric O⁺ ion’s distribution in the synchronous altitude region for different geomagnetic activity indexK _p is studied by theoretical modeling and numerical analyzing, and semi-empirical models for the O⁺ ion’s density and flux versus longitude in the synchronous altitude region for differentK _p are given. The main results show that in the synchronous altitude region: (i) The O⁺ ion’s density and flux in day-side are larger than those in nightside. (ii) With longitude changing, the higher the geomagnetic activity indexK _p is, the higher the O⁺ ion’s density and flux, and their variation amplitude will be. The O⁺ ion’s density and flux whenK _{p ≥} 6 will be about ten times as great as that whenK _p = 0. (iii) WhenK _p = 0 orK _{p ≥} 6, the O⁺ ion’s density reaches maximum at longitudes 120° and 240° respectively, and minimum in the magnetotail. WhenK _p = 3−5, the O⁺ ion’s density gets to maximum at longitude 0°, and minimum in the magnetotail. However, the O⁺ ion’s flux reaches maximum at longitude 120° and 240° respectively, and minimum in the magnetotail for anyK _p value.     

The solubility mechanisms of volatiles in silicate melts and their relations to crystal-andesite liquid equilibria

The anionic structure of magmatic liquids has been estimated at 1 atm and at pressures corresponding to those of the upper mantle. These estimates are based predominantly on spectroscopic data on binary metal oxide-silica and ternary metal oxide-silica-alumina melts. Structural information on melt compositions in aluminate-silica joins has been used to provide detailed information on the role of Al³⁺ in natural magma at atmospheric and high pressure.Regardless of pressure, andesitic melts may be described as combinations of chain, sheet, and three-dimensional network units. Nearly all Al³⁺ in the magmatic liquid resides in the three-dimensional network units. This Al³⁺ is locally charge-balanced with Na⁺, K⁺, Ca²⁺, and Mg²⁺. In the latter two cases, Al³⁺ and Si⁴⁺ are ordered, whereas for Na⁺ and K⁺, Si⁴⁺ and Al³⁺ are randomly mixed. Solution of water in natural magma results in the formation of new nonbridging oxygens in addition to OH groups attached to Si⁴⁺ and metal cations.On the basis of determined solution mechanisms of CO₂ and H₂O in silicate melts, thermodynamic properties of HO+CO₂, fluids and hydrous silicate melts and melting phase relations in peridotite-H₂O-CO₂, systems, it is found that natural andesitic magma in equilibrium with spinel Iherzolite in the upper mantle (10–20 kbar) must contain at least 5–7 wt.% H₂O. Andesitic magma with 5–7 wt.% H₂O in solution may be described as a mixture of Al-free three-dimensional units, sheets, and chains with a small proportion (less than 10%) of monomers.     

Mineralogy and crystal structure determination of Mg-fillowite

Accurate crystal structure determination of complex phosphate Mg-fillowite.—Na₂Ca(Mn₄Mg₂Fe)₇[PO₄]₆ has been finished. Mg-fillowite occurs in muscovite-pegmatite in the Altay area, Xinjiang, China. It is of the trigonal lattice, with space group R-3, unit cell a=1.5143(3)nm, c=4.3191(2) nm, V=8.5736 nm³, Z = 18. The R-factor of the determination accuracy is R (I<2 σ(I)) = 0.0776. The cationic polyhedra consist of two kinds of structural units: compound column and screw column. In three-dimensional space, they are linked with each other through corner sharing or edge sharing to build the entire structure. There are 45 atoms in an asymmetric unit and 720 atoms in a unit cell. In terms of types and connection manners of the coordination polyhedra, the degree of complicity is rarely seen in all mineral crystal structures. Study of such a structure is significant for the classification of phosphate and the building rule of coordination polyhedra in a complex crystal structure.     

Study of calibration function for surface wave magnitude of DK1 seismographs

ＳｔｕｄｙｏｆｃａｌｉｂｒａｔｉｏｎｆｕｎｃｔｉｏｎｆｏｒｓｕｒｆａｃｅｗａｖｅｍａｇｎｉｔｕｄｅｏｆＤＫ１ｓｅｉｓｍｏｇｒａｐｈｓＦＥＮＧＸＵＥ（薛峰）ＹＯＮＧＺＨＡＯ（赵永）ＣｅｎｔｅｒｆｏｒＡｎａｌｙｓｉｓａｎｄＰｒｅｄｉｃｔｉｏｎ，Ｓｔａｔ...     

Prediction of rockburst probability given seismic energy and factors defined by the expert method of hazard evaluation (MRG)

In this paper we suggest that conditional estimator/predictor of rockburst probability (and rockburst hazard, P ^T (t)) can be approximated with the formula P ^T (t) = P ₁(θ ₁)…P _N (θ _N )·P _dyn ^T (t), where P _dyn ^T (t) is a time-dependent probability of rockburst given only the predicted seismic energy parameters, while P _i (θ _i ) are amplifying coefficients due to local geologic and mining conditions, as defined by the Expert Method of (rockburst) Hazard Evaluation (MRG) known in the Polish mining industry. All the elements of the formula are (approximately) calculable (on-line) and the resulting P ^T value satisfies inequalities 0 ≤ P ^T (t) ≤ 1. As a result, the hazard space (0–1) can be always divided into smaller subspaces (e.g., 0–10⁻⁵, 10⁻⁵–10⁻⁴, 10⁻⁴–10⁻³, 10⁻³–1), possibly named with symbols (e.g., A, B, C, D, …) called “hazard states” — which saves the prediction users from worrying of probabilities. The estimator P ^T can be interpreted as a formal statement of (reformulated) Comprehensive Method of Rockburst State of Hazard Evaluation, well known in Polish mining industry. The estimator P ^T is natural, logically consistent and physically interpretable. Due to full formalization, it can be easily generalized, incorporating relevant information from other sources/methods.