Decompression and H<Subscript>2</Subscript>O exsolution driven crystallization and fractionation: development of a new model for low-pressure fractional crystallization in calc-alkaline magmatic systems

Magma ascent, decompression-induced H₂O exsolution and crystallization is now recognized as an important process in hydrous subduction zone magmas. During the course of such a process calculations suggest that the ascent rate of a degassing and crystallizing mafic magma will be greater than crystal settling velocities. Thus, any crystals formed as a consequence of volatile exsolution will remain suspended in the magma. If the magma erupts before the percentage of suspended crystals reaches the critical crystallinity value for mafic magma (~55 vol.%) it will produce the commonly observed crystal rich island arc basalt lava. If the magma reaches its critical crystallinity before it erupts then it will stall within the crust. Extension of compaction experiments on a 55 vol.% sand-Karo syrup suspension at different temperatures (and liquid viscosities) to the likely viscosities of interstitial andesitic to dacitic liquid within such a stalled magma suggest that small amounts (up to ~10%) can be expelled on a time scale of 1–10 years. The expelled liquid can create a new intermediate to silicic body of magma that is related to the original mafic magma via fractional crystallization. The short time scale for liquid expulsion indicate that decompression-induced H₂O exsolution and crystallization can be an important mechanism for fractional crystallization. Based on this assumption a general model of decompression-induced crystallization and fractionation is proposed that explains many of the compositional, mineralogical and textural features of Aleutian (and other andesites).     

Eutrophication,health risk assessment and spatial analysis of water quality in Gucheng Lake,China

The water characteristics of the Gucheng Lake, such as eutrophication, health and spatial distribution, were investigated. On the basis of the trophic state index (TSI) and entropy weight, a synthesized trophic state index (STSI) model was established to assess lake eutrophication condition through calculating STSI, choosing TP, TN, COD, BOD and NH₃-N as trophic variables. The STSI ranged from 50.58 to 62.44, which showed that the water has been between eutrophic and supereutrophic. A histogram was applied to health risk assessment which was analyzed from carcinogenic substances (Cr⁺⁶, As and Cd) and non-carcinogenic substances (hydroxybenzene, Pb, Hg, CN⁻ and NH₃), and the results showed that the former was much greater than the latter for effect. The total risk for each resident caused by all pollutants ranged from 5.18E-05 to 8.34E-05, which is far higher than the standard, recommended by Sweden Bureau of Environment Protection and Holland Ministry of Building and Environment Protection (1.0E-05). Cluster analysis was used to detect similarities and dissimilarities among the seven sampling sites and explain the observed clustering in terms of affected conditions. Twenty-one variables were used to divide seven sampling sites into three groups, namely, north lake, south lake and lake center.     

Microstructure characteristics and forming environment of late Quaternary Period loess in the Loess Plateau of China

Since the late Quaternary Period, the rapid rising of the northeastern Qinghai-Tibetan Plateau has greatly changed the East Asian climate and formed different regional monsoon climates. Various microstructures of the loess in northern China have formed under a variety of climatic conditions. The scanning electron microscope (SEM) and energy spectrum analysis are used in this paper to study the relationship between the indexes of loess microstructures and the forming climates, and the indexes include the diameter of sand grains, Ca/Fe value and the characteristic shapes of the loess from SEM microstructures pictures. In terms of their respective climatic circumstances, the major microstructures of regional loess may be classified into nine categories according to their forms and particle sizes, such as loess formed in near-sand source, loess of granule sand dropping, and loess formed in hot and rainy places, i.e. with the loess SEM pictures, climates of about 15 ka b.p. in different places may be evaluated. In conclusion, the clear and distinguishable loess microstructures and their indexes in the Loess Plateau of China can also deduce the evolution of climate, such as warm or cold and wet or dry, and even sandstorm activities when the loess is deposited.     

Efficient Modeling of Gravity Fields Caused by Sources with Arbitrary Geometry and Arbitrary Density Distribution

We present a brief review of gravity forward algorithms in Cartesian coordinate system, including both space-domain and Fourier-domain approaches, after which we introduce a truly general and efficient algorithm, namely the convolution-type Gauss fast Fourier transform (Conv-Gauss-FFT) algorithm, for 2D and 3D modeling of gravity potential and its derivatives due to sources with arbitrary geometry and arbitrary density distribution which are defined either by discrete or by continuous functions. The Conv-Gauss-FFT algorithm is based on the combined use of a hybrid rectangle-Gaussian grid and the fast Fourier transform (FFT) algorithm. Since the gravity forward problem in Cartesian coordinate system can be expressed as continuous convolution-type integrals, we first approximate the continuous convolution by a weighted sum of a series of shifted discrete convolutions, and then each shifted discrete convolution, which is essentially a Toeplitz system, is calculated efficiently and accurately by combining circulant embedding with the FFT algorithm. Synthetic and real model tests show that the Conv-Gauss-FFT algorithm can obtain high-precision forward results very efficiently for almost any practical model, and it works especially well for complex 3D models when gravity fields on large 3D regular grids are needed.     

Stochastic evaluation of simple pairing approaches to reconstruct incomplete rainfall time series

Two-station pairing approaches are routinely used to infill missing information in incomplete rainfall databases. We evaluated the performance of three simple methodologies to reconstruct incomplete time series in presence of variable nonlinear correlation between data pairs. Nonlinearity stems from the statistics describing the marginal peak-over-threshold (POT) values of rainfall events. A Monte Carlo analysis was developed to quantitatively assess expected errors from the use of chronological pairing (CP) with linear and nonlinear regression and frequency pairing (FP). CP is based on a priori selection of regression functions, while FP is based on matching the probability of non-exceedance of an event from one time series with the probability of non-exceedance of a similar event from another time series. We adopted a generalized Pareto (GP) model to describe POT events, and a t-copula algorithm to generate reference nonlinearly correlated pairs of random temporal distributions distributed according with the GP model. The results suggest that the optimal methodology strongly depends on GP statistics. In general, CP seems to provide the lowest errors when GP statistics were similar and correlation became linear; we found that a power-2 function performs well for the selected statistics when the number of missing points is limited. FP outperforms the other methods when POT statistics are different and variables are markedly nonlinearly correlated. Ensemble-based results seem to be supported by the analysis of observed precipitation at two real-world gauge stations.     

Long-term Geomechanical Observation in the Jeroným Mine

Geomechanical measurements have been carried out in the medieval Jeroným Mine since 2001. At first quarterly observation and later also continual monitoring by distributed measurement network were used. The network includes a number of different sensors that have been installed in selected places of underground spaces. There have not been any manifestations documenting dangerous changes of stability during the period of monitoring (2001–2015). The results of accomplished measurements from all convergence profiles confirm the total stability of all underground spaces. However, a number of places have been documented as potentially unstable. Measured values of selected parameters are presented in the paper, namely, movements along fractures, changes of stress tensor and fluctuations in the level of mine water table.