Discussion on the feature of strong earthquake: Orderly distribution in time, space and intensity before the Western Kunlun Mountain Pass M=8.1 earthquake

Introduction The Western Kunlun Mountain Pass M=8.1 earthquake occurred on November 14, 2001 is the other M=8 earthquake occurred 50 years after Dangxiong, Tibet M=8.0 earthquake in Chinese mainland. The earthquake has caused the attention of the seismologists in the following aspects: 1) The fracture length is more than 400 km, which is far away from the estimated length by the statistic empirical function between the magnitude and the fracture length (WANG, et al, 2002); 2) The aftersh…     

盆地岩石圈结构与油气成藏及分布

本文综述了大陆岩石圈研究现状和克拉通盆地、裂谷盆地和前陆盆地的岩石圈结构特征,指出在古裂谷、褶皱带或区域性深断裂等陆壳构造薄弱带上发育起来的多期叠合盆地,具有很好的含油气前景。大型含油气盆地往往存在地幔上隆、地壳减薄和地壳内低速层,盆地基底沉降与盖层沉积厚度较大。适度的后期构造活动改造和岩浆活动有利于沉积盆地内油气生成与保存。     

Melt/biotite B/B isotopic fractionation and the boron local environment in the structure of volcanic glasses

This paper is focused on the role of boron coordination in determining the ¹¹B/¹⁰B isotopic fractionation between melt/glass and biotite at magmatic temperatures. For this purpose, three evolved volcanic rocks from Roccastrada, Mt. Amiata, and Mt. Cimini belonging to the Neogene-Quaternary magmatism of central Italy were studied. In these samples, the measured boron biotite-glass partition coefficient ranges between 0.004 and 0.011, indicating that boron behaves as an incompatible element during biotite crystallization. The ¹¹B magic-angle spinning nuclear magnetic resonance (NMR) spectra reveal the presence of trigonal BO_3/2 units, tetrahedral BO_4/2⁻ sites, and three-coordinated BO_2/2O⁻ species containing one nonbridging oxygen. The relative contributions of these different boron sites were estimated by spectral deconvolution, and it was observed that the fraction of trigonally coordinated boron decreases with increasing K₂O concentration in the glass. The ¹¹B/¹⁰B isotopic fractionation between biotite and melt/glass was observed to be large even at magmatic temperatures and was found to be 1.0066 (Roccastrada sample), 1.00535 (Mt. Amiata sample), and 1.00279 (Mt. Cimini sample). Fractionation is mostly related to the relative amount of trigonal and tetrahedral boron sites in the glass network rather than to other processes, including the speciation of hydrous species in the glass structure. The measured α values are significantly higher than the calculated ones obtained using the reduced partition function ratios (RPFRs) for B(OH)₃ and B(OH)₄⁻ as reported by Kakihana et al. (1977) and the abundance of trigonal and tetrahedral boron obtained by ¹¹B NMR spectra. Furthermore, a nonlinear relationship is observed between the percentage of BO₄ in the glass structure and the measured 1000lnα, suggesting that the approximation of monomeric B(OH)₃ and B(OH)₄⁻ species contributions through ideal mixing in calculating the RPFRs in polyanions (Oi et al., 1989) probably does not apply to silicate glasses.The large B isotopic fractionation measured between glass and biotite and its dependence on the boron coordination in the glass are a limitation to the use of δ¹¹B in the mineral to characterize magmas. Nonetheless, the high incompatible behavior of boron in the most common magmatic minerals rules out that fractional crystallization significantly modified the B isotopic composition of the melt.     

Evidence for fullerenes in solid bitumen from pillow lavas of Proterozoic age from Mítov (Bohemian Massif, Czech Republic)

Andesitic pillow lavas containing biogenic, solid bitumen (SB) are a constituent of a Neoproterozoic volcanosedimentary sequence (Teplá-Barrandian unit, Bohemian Massif) in the Mítov area of the Czech Republic. A black shale formation that is crosscut by these andesitic basalts is 565 Ma old. Carbon disulfide extracts of two powdered samples of SB contain 0.2 and 0.3 ppm of C₆₀, respectively, as determined by high-pressure liquid chromatography. The peak assignment based on retention time is fully supported by high-resolution electron ionization mass spectrometry (EI-MS). No C₇₀ was detected, nor was C₆₀ found in two other SB samples from this locality. Other investigated carbonaceous samples from Bohemia (coals and anthracites of Upper Paleozoic age and anthraxolite, graphitoids, and graphite of Upper Proterozoic age) did not contain fullerenes at concentrations above the detection limit of 0.01 ppm. The absence of C₆₀ in these samples was confirmed by EI-MS. The proposed mechanism of fullerene formation involves a primary algal phase, generation of a hydrocarbonaceous mixture in the course of thermal evolution of the sedimentary series, and their high-temperature transformation related to the extrusion of basalt. An important feature for fullerene conservation was the enclosure of fullerenes in SB with a structure similar to glasslike carbon, where the fullerene was protected against oxidation.     

Pressure-Dependence of Rare Earth Element Distribution in Amphibolite- and Granulite- Grade Garnets. A LA-ICP-MS Study

Using an excimer (KrF) laser ablation ICP-MS system, we studied the distribution of REE in garnets from metapelites and metabasites from Ivrea-Verbano (Western Alps, Italy) and from the Peña Negra Anatectic Complex (Central Iberia), finding systematic variations that correlate well with the metamorphic grade. Chondrite-normalized REE patterns of garnets from amphibolite-grade metapelites have lower-than-chondrite levels from La to Sm, a very small or no Eu anomaly, and a steep rise in the abundance of heavy REE as the atomic number increases. Metapelitic garnets from the amphibolite-granulite transition have a marked Eu negative anomaly and are enriched in MREE such that Sm is 10-15 times chondrite and the pattern is almost flat from Dy to Yb-Lu. In garnets from granulite-grade metapelites, the intensity of the Eu anomaly and the relative concentration of Nd, Sm, Gd and Tb increase, with almost flat chondrite-normalized patterns from Sm to Lu. Garnets from mafic granulites are remarkably similar to those of metapelitic garnets equilibrated at the same pressure, except for the Eu anomaly. The apparent paradox of enhanced uptake of larger REE ions with increasing pressure is attributed to the 3M²⁺ 2REE³⁺+ vacancy substitution, which produces a net decrease in the dimensions of the unit-cell of garnet. Variations in REE patterns depend essentially on the pressure and have little dependence on either temperature, bulk-composition of garnet, or REE whole-rock composition, so they could represent a new approach for geobarometric studies. The best numerical parameter to express pressure-related variations of REE distribution in garnets is the Gd/Dy ratio which does not seem perceptibly affected by disequilibrium partitioning. The regression equation between GASP pressure and the average Gd/Dy^garnet is P = 3.6 + 5.6 Gd/Dy. This equation seems to be reliable for garnets: (1)equilibrated within a pressure range of 4-9 kbar, (2) coexisting with modal monazite; and (3) with unit-cell dimensions under 11.46 Å.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.