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.     

Pan evaporation modelling and changing attribution analysis on the Tibetan Plateau (1970–2012)

Pan evaporation (E_p) is an important indicator of water and energy and the decline of E_p has been reported in many regions over the last decades. The climate and E_p are dependent on each other. In this study, the temporal trends of E_p and main E_p drivers, namely mean air temperature (T_a), wind speed (u), global solar radiation (R_s), net long‐wave radiation(R_nl) and vapour pressure deficit (D) from 1970 to 2012, were calculated on the basis of 26 meteorological stations on the Tibetan Plateau. The arithmetic average of E_p from 26 stations decreased with the rate of ?11.91 mm a^?2; the trends of R_s, R_nl, T_a, u and D were ?1.434 w m^?2 decade^?1, 0.2511 w m^?2 decade^?1, 0.3590°C decade^?1, ?0.2376 m s^?1 decade^?1 and 9.523 Pa decade^?1, respectively. The diffuse irradiance is an essential parameter to model E_p and quantify the contribution of climatic factors to changing E_p. 60 724 observations of R_s and diffuse solar irradiance (R_d) from seven of the 26 stations were used to develop the correlation between the diffuse fraction (R_d/R_s), and the clearness index (R_s/R_o). On the basis of the estimation of the diffuse component of R_s and climatic data, we modified the PenPan model to estimate Chinese micro‐pan evaporation (E_p) and assess the attribution of E_p dynamics using partial derivatives. The results showed that there was a good agreement between the observed and calculated daily E_p values. The observed decrease in E_p was mostly due to declining wind speed (?13.7 mm a^?2) with some contributions from decreasing solar irradiance (?3.1 mm a^?2); and the increase of temperature had a large positive effect (4.55 mm a^?2) in total whilst the increase of R_nl had insignificant effect (0.35 mm a^?2) on E_p rates. The change of E_p is the net result of all the climatic variables. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling Rupture Dynamics of a Planar Fault in 3-D Half Space by Boundary Integral Equation Method: An Overview

In this article, we first reviewed the method of boundary integral equation (BIEM) for modelling rupture dynamics of a planar fault embedded in a 3-D elastic half space developed recently (ZHANG and CHEN, 2005a,b). By incorporating the half-space Green's function, we successfully extended the BIEM, which is a powerful tool to study earthquake rupture dynamics on complicated fault systems but limited to full-space model to date, to half-space model. In order to effectively compute the singular integrals in the kernels of the fundamental boundary integral equation, we proposed a regularization procedure consisting of the generalized Apsel-Luco correction and the Karami-Derakhshan algorithm to remove all the singularities, and developed an adaptive integration scheme to efficiently deal with those nonsingular while slowly convergent integrals. The new BIEM provides a powerful tool for investigating the physics of earthquake dynamics. We then applied the new BIEM to investigate the influences of geometrical and physical parameters, such as the dip angle (δ) and depth (h) of the fault, radius of the nucleation region (R_asp), slip-weakening distance (D_c), and stress inside (T_i) and outside (T_e) the nucleation region, on the dynamic rupture processes on the fault embedded in a 3-D half space, and found that (1) overall pattern of the rupture depends on whether the fault runs up to the free surface or not, especially for strike-slip, (2) although final slip distribution is influenced by the dip angle of the fault, the dip angle plays a less important role in the major feature of the rupture progress, (3) different value of h, δ, R_asp, T_e, T_i and D_c may influence the balance of energy and thus the acceleration time of the rupture, but the final rupture speed is not controlled by these parameters.     

Hydrogen isotope composition of natural gases from the Tarim Basin and its indication of depositional environments of the source rocks

By measuring carbon and hydrogen isotope compositions for C1, C2 and C3 of 74 gas samples, natural gases from the Tarim Basin can be divided into six groups on the basis of their origins: (1) coal-type gas derived from coal measures; (2) coal-type gas generated from the T-J lacustrine mudstones; (3) oil-type gas derived from the Cambrian and low Ordovician marine source rocks; (4) oil-type gas from the source rocks deposited in the marine-transitional facies; (5) mixing gas between gas derived from the Carboniferous transitional source rocks and the Mesozoic humic gas, and (6) mixing gases of thermal genetic gas and little deep gas in the Southwest depression of the Tarim Basin. The δ D values of methane in natural gases originating from different type kerogens are affected by both palaeo-environments of the source rock formation (kerogen types) and thermal maturity, with sedimentary environment (kerogen type) as the main controlling factor. Under the similar thermal maturity, the hydrogen isotope composition of methane is more enriched in deuterium in marine environments than lacustrine one. With the increase of thermal maturity and the increase of carbon atomic numbers of gaseous alkanes, the hydrogen isotopes become enriched in deuterium. The δ D values of ethane and propane (δ D2, δ D3) are controlled mainly by thermal maturity and to a lesser degree by sedimentary environment of the source rock formation. The partial reversal of hydrogen isotopes for gaseous alkanes would be related to the microbial oxidation, mixing of sapropelic and humic gases and / or mixing of gases from similar kerogen sources with various thermal maturities. In the oil-type gas, the sulfate reduction reaction would result in the reversed order of δ D1 and δ D2 (e.g. δ D1>δ D2).     

Fractal approach of the temporal earthquake distribution in the Hellenic arc-trench system

The time clustering of earthquakes occurring in the Hellenic arc-trench system is quantitatively analyzed by means of the fractal dimension,D, of their time distribution in the time intervals of 1950–1985 (M _s >-4.5) and 1964–1985 (M _s 4.0). The results obtained imply that scale-invariant clustering holds over very large scale lengths of time,T, with 2²–2⁸T (in min) 2²⁰–2²², depending on the seismotectonic segment considered. In all segments a common feature is the relation between theD ₁,D ₂ andD ₃-values found for shallow, intermediate-depth and all-depth shocks, respectively:D ₃>D₁>D₂. TheD-values found for shallow shocks range between 0.137 and 0.191 with the exception of the Ionian Islands and Cretan segments where anomalously high values (D=0.221–0.251) have been determined. We discuss possible seismotectonic interpretations of the results.     

An extended theory of equi-risk line for a general release rule

The expression of equi-risk line derived by the authors represents the relationship between discharge capacityy ₀ ^u and storage capacityz ₀ ^u to keep flood frequency under a certain risk level represented by the return periodT, i.e.,z ₀/z ₀ ^u ={(y ₀ ^u –y ₀ ^u )/y ₀ ^u } ^S , wherey ₀ ^u andz ₀ ^u areT-year probability peak discharge and total volume of a hydrograph. The shape parametersS is evaluated in this paper for various release rules of the storage facilities and correlations of durations and peaks of hydrographs. The expression forS is: , whereS ₀ andS are the values ofS forp=0 and , andp is the exponent of a general storage-release relation,q=az' ^p, wherea is the storage constant, andz' andq are the volume of stored water and the corresponding release. The values ofS ₀ andS are expressed in terms of the correlation coefficient of durations and peaks of inflow hydrographs.     

Carbon and hydrogen isotopic characteristics of natural gases from the Luliang and Baoshan basins in Yunnan Province, China

The Luliang and Baoshan basins are two small ba- sins in Yunnan Province. In the recent ten years or so, there have been found a number of natural gas pools of commercial importance in the two basins. Although the gas pools are small in size, the natural …     

Properties of electric turbulence in the polar cap ionosphere

Small-scale (scales of ∼0.5–256 km) electric fields in the polar cap ionosphere are studied on the basis of measurements of the Dynamics Explorer 2 (DE-2) low-altitude satellite with a polar orbit. Nineteen DE-2 passes through the high-latitude ionosphere from the morning side to the evening side are considered when the IMF z component was southward. A rather extensive polar cap, which could be identified using the ɛ-t spectrograms of precipitating particles with auroral energies, was formed during the analyzed events. It is shown that the logarithmic diagrams (LDs), constructed using the discrete wavelet transform of electric fields in the polar cap, are power law (μ ∼ s ^α). Here, μ is the variance of the detail coefficients of the signal discrete wavelet transform, s is the wavelet scale, and index α characterizes the LD slope. The probability density functions P(δE, s) of the electric field fluctuations δE observed on different scales s are non-Gaussian and have intensified wings. When the probability density functions are renormalized, that is constructed of δE/s ^γ, where γ is the scaling exponent, they lie near a single curve, which indicates that the studied fields are statistically self-similar. In spite of the fact that the amplitude of electric fluctuations in the polar cap is much smaller than in the auroral zone, the quantitative characteristics of field scaling in the two regions are similar. Two possible causes of the observed turbulent structure of the electric field in the polar cap are considered: (1) the structure is transferred from the solar wind, which is known to have turbulent properties, and (2) the structure is generated by convection velocity shears in the region of open magnetic field lines. The detected dependence of the characteristic distribution of turbulent electric fields over the polar cap region on IMF B _y and the correlation of the rms amplitudes of δE fluctuations with IMF B _z and the solar wind transfer function (B _y ² + B _z ²)^1/2sin(θ/2), where θ is the angle between the geomagnetic field and IMF reconnecting on the dayside magnetopause when IMF B _z < 0, together with the absence of dependence on the IMF variability are arguments for the second mechanism.     

Characteristics of low-frequency oscillation intensity of air-sea turbulent heat fluxes over the northwest Pacific

Based on the daily turbulent heat fluxes and related meteorological variables datasets (1985–2006) from Objectively Analyzed air-sea Fluxes (OAFlux) Project of Woods Hole Oceanographic Institution (WHOI), characteristics of low-frequency oscillation intensity of air-sea turbulent heat fluxes over the northwest Pacific are analyzed by linear perturbation method and correlation analysis. It can be concluded that: 1) the distribution of low-frequency oscillation intensity of latent heat flux (LHF) over the northwest Pacific is mainly affected by that of low-frequency oscillation intensity of anomalous air-sea humidity gradient (Δq′) as well as mean air-sea humidity gradient (), while the distribution of low-frequency oscillation intensity of sensible heat flux (SHF) is mainly affected by that of low-frequency oscillation intensity of anomalous air-sea temperature gradient (ΔT′). 2) The low-frequency oscillation of turbulent heat fluxes over the northwest Pacific is the strongest in winter and the weakest in summer. And the seasonal transition of low-frequency oscillation intensity of LHF is jointly influenced by those of low-frequency oscillation intensity of Δq′, low-frequency oscillation intensity of anomalous wind speed (U′), and mean wind speed (Ū), while the seasonal transition of low-frequency oscillation intensity of SHF is mainly influenced by those of low-frequency oscillation intensity of ΔT′ and Ū. 3) Over the tropical west Pacific and sea areas north of 20°N, the low-frequency oscillation of LHF (SHF) is mainly influenced by atmospheric variables q _a ′ (T _a ′) and U′, indicating an oceanic response to overlying atmospheric forcing. In contrast, over the tropical eastern and central Pacific south of 20°N, q _s ′ (T _s ′) also greatly influences the low-frequency oscillation of LHF (SHF). Supported by National Natural Science Foundation of China (Grant No. 40675028) and National Basic Research Program of China (Grant No. 2006CB403600)     

Regional characteristics of coda attenuation after the lancang-Gengma earthquake

In the light of the single scattering model of coda originating from local earthquakes, and based on the aftershock coda registered respectively at the 4 short period stations installed near the foci shortly after theM7.6 Lancang andM7.2 Gengma earthquakes, this paper has tentatively calculated the rate of amplitude attenuation and theQ _c-value of the coda in the Lancang and Gengma areas using a newly-founded synthetic determination method. Result of the study shows the rate of coda amplitude attenuation demonstrates remarkable regional differences respectively in the southern and northern areas. The southern area presents a faster attenuation (Q _c=114), whereas the northern area shows a slower attenuation (Q _c=231). The paper also discusses the reasons causing such differences. Result of the study also suggests a fairly good linear relation between the coda source factorA _o(f) and the seismic moment and the magnitude. Using the earthquake scaling law, the following formulas can be derived: lgM ₀=lgA ₀(f)+17.6,M _D=0.67lgA ₀(f)+1.21 and logM ₀=1.5M _D+15.79. In addition, the rates of amplitude attenuationβ _s andβ _m are respectively calculated using the single scattering and multiple scattering models, and the ratioβ _s/β_m=1.20−1.50 is found for the results respectively from the two models. Finally, the mean free pathL of the S-wave scattering in the southern and northern areas are determined to be 54 km and 122 km respectively by the relations which can distinguish between the inherentQ _i and scatteringQ _s, testify to this areas having lowQ-values correspond to stronger scatterings. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 71–82, 1992. This study is partly supported by the Seismological Science Foundation of the State Seismological Bureau of China, and the present English version of the paper is translated from its Chinese original by Wenyi Xia, Seismological Bureau of Yunnan Province.     

Variations of the geomagnetic field and the geomagnetic effects of three earthquakes ( M s>5.0) in Jiangsu Province, China

Since 1979 the repeated observations and experiments of geomagnetic total intensity and vertical component have been carried out for ten years in the geomagetic network which is located in Jiangsu Province, China. Three earthquakes aboveM _s 5.0 occurred during the decade, and some seismomagnetic effects were observed. The observation results show that the anomalies of the vertical geomagnetic component can’t be observed untill some months before the earthquake (M _s>5.0) in this area. In this paper it is suggested that a densely distributed network for continuous observation of geomagnetic vertical component may catch seismomagnetic anomalies and thus improve earthquake prediction in the light of the geomagnetic measurements of the mid — or — low latitude locations. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 80–87, 1991. This study is sponsored by the Chinese Joint Seismological Science Foundation.     

Fluid mixing as the mechanism of formation of the Dajing Cu-Sn-Ag-Pb-Zn ore deposit,Inner Mongolia

Dajing Cu-Sn-Ag-Pb-Zn ore deposit, in the Inner Mongolia Autonomous Region of China, is a fissure-filling hydrothermal ore deposit. The δD values of quartz-hosted inclusion water are centered at −100%.– −130%.. The δ³⁴S values of sulfide ore minerals and δ¹³ C values of carbonate gangue minerals vary from −0.3%. to 2.6%. and from −2.9%. to −7.0%., respectively. Integrated isotopic data point to two major contributions to the mineralizing fluid that include a dominant meteoric-derived groundwater, and sulfur and carbon species from hypogene magma. Linear trends are exhibited on the gaseous H₂O versus CO₂ plot, and plots of CO, N₂, CH₄, and C₂H₆. It is shown by quantitative simulation that magma degassing cannot explain the linear trends. Hence, these linear trends are interpreted in terms of mixing of CO₂-rich magmatic fluid with meteoric-derived groundwater. The groundwater circulated in Paleozoic sedimentary rocks and absorbed CO, N₂, CH₄, C₂H₆ and radiogenic Ar from organic matter. Cooling effects resulting from mixing have caused the precipitation of ore minerals.     

Modelling snowpack surface temperature in the Canadian Prairies using simplified heat flow models

Three practical schemes for computing the snow surface temperature T_s, i.e. the force–restore method (FRM), the surface conductance method (SCM), and the Kondo and Yamazaki method (KYM), were assessed with respect to T_s retrieved from cloud‐free, NOAA‐AVHRR satellite data for three land‐cover types of the Paddle River basin of central Alberta. In terms of R², the mean T_s, the t‐test and F‐test, the FRM generally simulated more accurate T_s than the SCM and KYM. The bias in simulated T_s is usually within several degrees Celsius of the NOAA‐AVHRR T_s for both the calibration and validation periods, but larger errors are encountered occasionally, especially when T_s is substantially above 0 °C. Results show that the simulated T_s of the FRM is more consistent than that of the SCM, which in turn was more consistent than that of the KYM. This is partly because the FRM considers two aspects of heat conduction into snow, a stationary‐mean diurnal (sinusoidal) temperature variation at the surface coupled to a near steady‐state ground heat flux, whereas the SCM assumes a near steady‐state, simple heat conduction, and other simplifying assumptions, and the KYM does not balance the snowpack heat fluxes by assuming the snowpack having a vertical temperature profile that is linear. Copyright © 2005 John Wiley & Sons, Ltd.     

Stable isotopic studies of dinosaur eggshells from the Nanxiong Basin, South China

Stable carbon and oxygen isotopes have been observed in dinosaur eggshell samples,identified as Macroolithus yaotunensis,collected in two sections of the Nanxiong Basin,Guangdong Province,South China.Multiple positiveδ18O perturbations that occurred during the K/T transition of about 150 ka suggest that there may have been at least three periods of the extreme dry climate with a mean annual air temperature over 27℃.δ13C ranges from -8.37‰ to -11.60‰ of the dinosaur eggshells,added to the metabolic enrichment of 16‰,may indicate that the dinosaurs represented by this type of eggshells probably consumed the plant species of C3 type with the δ13C values of -24.3‰- -27.6‰.Moreover,the trend toward more negative δ13C values of the eggshells during the K/T transition could be an indirect consequence of a lasting increase in atmospheric CO2 concentration.     

The influence of stiffness of the true triaxial testing machine on the rupturing and acoustic emission of rocks and its correlation with seismic activities

This paper introduces how the ratioR of the characteristic stiffness of rock samplevs. the stiffness of testing machines would influence the rupturing process and the acoustic emission (AE) on the part of the tested rock samples. Result of the experiment shows: WhenR>0.20, the rock sample would rupture abruptly; whereas whenR<0.20, the rock samples would rupture slowly. When the samples rupture abruptly, the time-dependent variation of the AE rate takes such a pattern:peak value—stable low values—rises to the maximum value (concentration)—drops back to the minimum value (quiescence)—(rises again)—ruptures. Moreover, smallerR-value tends to be associated with longer quiescence and vice versa. WhenR>1.50, no pre-failure quiescence is detected. When the rock samples rupture slowly, the variation pattern of the AE rate (after the stress has increased to more than 50% of the rupturing stress) is as the following:stable low (or high) values—rises (or drops) to its maximum (or minimum) values and then continues for some time—ruptures. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 223–233, 1991. This study is supported by the Chinese Joint Seismological Science Foundation.     

Reference evapotranspiration trends and their sensitivity to climatic change on the Tibetan Plateau (1970–2009)

Evapotranspiration is an important component of the water and energy balance. It is dependent on climate. Precipitation, solar radiation, temperature, humidity, and wind all contribute to the rate of evapotranspiration. In this study, the temporal trends of reference evapotranspiration (ET_ref) and four main ET_ref drivers, namely, mean air temperature (T_a), wind speed (u₂), net radiation (R_n) and actual vapour pressure (e_a) from 1970 to 2009, were calculated based on 75 meteorological stations on the Tibetan Plateau. The results showed that the ET_ref on the Tibetan Plateau decreased on average by 0.6909 mm a^‐1a^‐1 from 1970 to 2009. T_a and e_a showed an increasing trend, whereas u₂ and R_n exhibited a decreasing trend. To explore the underlying causes of the ET_ref variation, an attribution analysis was performed to quantify the contribution of T_a, u₂, R_n and e_a, which showed that the changes in u₂, R_n and e_a produced the negative effect, whereas T_a produced the positive effect on ET_ref rates. The changes in u₂ were found to produce the largest decrease (?0.7 mm) in ET_ref, followed by e_a (?0.4 mm) and R_n (?0.1 mm). Although the significant increase in T_a had a large positive effect (0.51 mm) on ET_ref rates, changes in the other three variables each reduced ET_ref rates, resulting in an overall negative trend in ET_ref. Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of maximal daily average values of relativistic electron fluxes in geostationary orbit during the magnetic storm recovery phase

The relation of the maximal daily average values of the relativistic electron fluxes with an energy higher than 2 MeV, obtained from the measurements on GOES geostationary satellites, during the recovery phase of magnetic storms to the solar wind parameters and magnetospheric activity indices has been considered. The parameters of Pc5 and Pi1 geomagnetic pulsations and the relativistic electron fluxes during the prestorm period and the main phase of magnetic storms have been used together with the traditional indices of geomagnetic activity (A _E, K _p, D _st). A simple model for predicting relativistic electron fluxes has been proposed for the first three days of the magnetic storm recovery phase. The predicted fluxes of the outer radiation belt relativistic electrons well correlate with the observed values (R ∼ 0.8–0.9).     

Seasonal climate patterns and their influence on calibration of the Hargreaves-Samani equation

Abstract

Annual patterns in climate parameters were studied to evaluate how these influence the quality of reference evapotranspiration (ETo) estimates obtained from the Hargreaves-Samani (HS) equation, since the method only uses the measured temperature directly. The work evaluates how these patterns can be used to improve the HS ETo estimates. Ten-year moving averages from a set of California Irrigation Management Information System (CIMIS) stations were used to evaluate the relationships between solar radiation (R_s), temperature (T) and ETo. The results indicate that T treads behind solar radiation and its value peaks some 25 days later. Thus, the main irrigation season in the Mediterranean climate (1 May–30 September) can be divided into three phases: increasing R_s and T; decreasing R_s with increasing T; and decreasing R_s and T. Non-univocal annual cycles were observed between R_s and T, ETo and R_s, and ETo and T. These annual patterns result in important seasonal changes in the ratio between the HS and Penman-Monteith (FAO PM) ETo estimates. The changes are particularly important during the irrigation season, where the FAO PM initially calculates greater ETo values than the HS methodology, and from the end of May to early September, where the HS equation overestimates the ETo values (by 17 mm, or 3%). These patterns obtained from 2000–2009 data were used to calibrate and improve HS ETo estimates at new sites for the 2010–2011 period. Calibration based on the proposed seasonal region-wide FAO PM/HS ETo ratios improved both the bias (decreased from 0.40 to 0.36 mm d^-1) and r² (increased from 0.67 to 0.87) of the ETo estimates for the irrigation season. The proposed methodology can be easily applied to other regions, even when the existing weather stations are sparse.

Editor Z.W. Kundzewicz     

Variation ofS q-ranges andL 2 (H) with the dip latitude in Indian region

Summary A modified method of computing theS _q -range in terms ofR _x inH, D andZ has been suggested. Mean quiet-day daily rangesR _x , have been computed for five Indian stations usingH, D andZ data for selected years. The seasonal and latitudinal variation ofR _x (H), R _x (Z) andR _x (D) are discusses. Thee-season maximum inR _x (H) andR _x (Z) is attributed to the decrease in the distance between the foci during equinoxes; the electrojet and theS _q -foci movement with seasons have little influence onS _q (D). It is inferred that the electrojet current is independent of the worldwideS _q current system and stands with its own return currents.The variation of lunar semi-diurnal tide inH[L ₂ (H)] with the dip latitude in Indian region shows a secondary peak at 9° N dip latitude. This secondary peak in theL ₂ (H) is termed as Lunar secondary electrojet, and it is suggested that this is possibly produced by magneto-acoustic oscillations due to the drift motion of the charged particles that produce the primary jet in a direct transverse to itself.     

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.