Active mud volcanoes on- and offshore eastern Makran, Pakistan

To study the activity, source and driving force of the venting of fluidized mud in the coastal Makran, we carried out reconnaissance surveys of two active onshore mud volcano fields (Chandragup and Jebel-u-Ghurab) and of a newly born (March 1999) offshore mud volcano (Malan Island). All studied on- and offshore mud volcanoes line up along one southwest/northeast-trending structural lineament, the axis of the Dhak Anticline. Isotopic data point to a bacterial origin of the gas (mainly methane). Mixed benthic foraminiferal faunas and calcareous nannofloras suggest that the source level of the extruded mud is at a sub-surface depth of 2-3 km. Observed mud discharge varied between 0 and 1.4 m³ h^-1 and gas discharge rates between negligible amounts to 1 m³ s^-1. Mud temperatures in the crater lake of Chandragup I are near-ambient temperatures. They rise slightly (Б.5°C) during episodes of modest mud outflow. An area of 160,000 m³ of soft mud was vigorously extruded from the sea floor at a water depth of 10 m within days after 15 March 1999, forming Malan Island. The island was destroyed within a few months after its birth by deep-reaching wave action during the SW monsoon. This was possibly aided by local subsidence of the mud volcano structure due to the volume loss following rapid degassing and mud extrusion.     

Intergrowth and twinning in opal-CT lepispheres

Under conditions of unobstructed diagenetic growth, opal-CT usually forms spherical aggregates, so-called lepispheres. From a scanning electron microscope study of deep-sea porcelanites, the development of these structures is described, from isolated opal-CT platelets via complete lepispheres to large coalescent aggregates. The conspicuous symmetric ultrastructure of lepispheres—reflected by the constancy of angles between interpenetrating crystal blades—is crystallographically determined: The lepispheres consist of groups of opal-CT blades intergrown according to the (30–34) and (10–16) twinning laws of tridymite. The same rules of intergrowth can be applied to the morphology of hydrothermally synthesized opal-CT lepispheres.     

Geology and geochemistry of D–O–C isotope systematics of the Qolqoleh gold deposit, Northwestern Iran: Implications for ore genesis

The Qolqoleh gold deposit is located in northwestern part of the Sanandaj–Sirjan metamorphic belt, northwestern Iran. Igneous and sedimentary units exposed in the area have undergone greenschist metamorphism. The area was affected by a NE–SW trending shear zone and subsequent deformation. Two different types of mineralization are distinguished in the Qolqoleh gold deposit based on geological–structural conditions indicated by microtextural analysis: ductile and then brittle. Ore-forming processes are divided into three stages: Early (I), Middle (II) and Late (III), which include quartz–pyrite (I), sulfides and gold (II) and carbonate veinlets (III), respectively. The stage I fluids are characterized by δ¹⁸O = 15.5‰ at 440 ºC, and are thought to be deep-sourced metamorphic waters; the stage III fluids, with δ¹⁸O = 1.6‰, are shallow-sourced meteoric waters; whereas, the stage II fluids, with δ¹⁸O = 13.1‰, are a mixture of deep-sourced metamorphic and shallow-sourced meteoric fluids. Based on comparisons of the D–O–C isotopic systematics, the ore-forming fluids with characteristic high δ¹⁸O and δ¹³C and low δD originated from metamorphic devolatilization of Cretaceous volcano-sedimentary (felsic to mafic metavolcanic rocks–shale–carbonate–carbonaceous chert) sequences, locally rich in organic matter. During late Cretaceous continental collision of the Afro-Arabian continent and the Iranian microcontinent, a crustal slab consisting of felsic to mafic metavolcanic rocks, carbonate, shale and carbonaceous chert was underthrust northwards beneath the central Iranian microcontinent along the Zagros fault. During further contraction, deformation was localized in reverse oblique-slip structures with vergence toward south; shear zones generally follow contacts between more competent and less competent rock units. Metamorphic devolatilization of this underthrust slab is the source of the ore-forming fluids that generated the Au ore belt, which includes the Qolqoleh gold deposit.     

Comparison of plant species diversity with different plant communities in deciduous forests

Species diversity is one of the most important indices used for evaluating the sustainability of forest communities. This study aims to characterize the forest communities and to identify and compare the plant species diversity in the study area. For this purpose, 152 relevés were sampled by a randomized-systematic method, using the Braun-Blanquet scale. Classification of the vegetation was conducted by the twinspan algorithm. Four communities, including Querco-Carpinetum betulii, Carpineto-Fagetum Oriental, Rusco-Fagetum Oriental and Fagetum Oriental were recognized. Species richness, Shannon, and Simpson indices were applied to quantify diversity of the different communities. Turkey test was used to investigate the differences in the species richness, diversity and evenness indices among the different communities. The results illustrate that Querco-Carpinetum betulii and Carpineto-Fagetum Oriental communities are significantly more diverse than Rusco-Fagetum Oriental and Fagetum Oriental communities. The spatial structure of the releves becomes more ‘homogenous’ and the dominance structure changes: the proportion of beech-forest species is gradually increasing. At the same time, the number of species per unit area decreases constantly, reaching eventually the value comparable to that recorded for hornbeam forest. Generally, species diversity is inversely correlated with the dominance of shade tolerant climax species.     

Magnetic stratigraphy of the Permian-Triassic traps in the Kotui River valley (Siberian Platform): New paleomagnetic data

The Permian-Triassic effusive traps in the Kotui River valley (Siberian Platform) were studied in detail by a paleomagnetic method, resulting in a summarized magnetic stratigraphic section of the studied sequence. The presence of the reversed polarity zone corresponding to the Khardakh Formation was argued for the section basement. Inside the Kogotok Group, the boundary between the direct and reversed polarity zones occurs in the lower part of the Onkuchak Formation and mismatches the boundary between the formations of the Kogotok Group, as was accepted before. These results contradict the presence of the transition record between the direct and reversed polarity zones in basalts of the lower part of the Onkuhcak Formation. The strong transitional interval between the Ivakinsk Formation and the upper parts of the Nadezhda Formation of the Norilsk section has no analogs in the Kotui section. This means the absence of any significant volcanic activity in the Maimecha-Kotui province during the intense eruptions in the Norilsk region resulting in origination of the greater part of the lower formations of the trap sequence. These data and possible correlations of the traps of the Kotui River valley and Norilsk region indicate that the Norilsk and Maimecha-Kotui sections were temporally overlapped.     

Numerical simulation of crack propagation in layered formations

In the present study, a boundary element method based on the higher order displacement discontinuity formulation is presented to solve the general problem of hydraulic fracture propagation in layered formations. Displacement collocation technique is employed to model the higher order displacement variation along the crack and the special crack tip element near its ends. The hydraulic fracture propagation and its interaction with the layer interface in non-homogenous rock materials are studied by the proposed semi-analytical (hybridized boundary element-boundary collocation) method. The maximum tangential stress criterion (or σ-criterion) of fracture mechanics considering different elastic constants (Young modulus and Poisson’s ratio) is used to obtain the fracture path. The fracture propagation from stiff to soft and soft to stiff media for cracks having different inclination angles is modeled, and the effects of elastic constants on the hydraulic fracture propagation is studied. The results show that if the hydraulic fracture originates in the stiffer layer, its capability to cross the layer increases and is vice versa for the softer material. The comparison of the results gained from the numerical method with those in the literature show a good performance of the method in the case of propagation of hydraulic fracture in layered formations.     

Relationship between porosity and permeability with stress using pore volume compressibility characteristic of reservoir rocks

Effective stress is one of the most important parameters which strongly affects pore volume compressibility curve, which is used in assessing the reservoir rock properties. It causes change in porosity as well as permeability of the reservoir rocks. In the present study, pore volume compressibility characteristics of the reservoir rocks at different effective stresses were used to derive the relationship of porosity and permeability with the effective stress. To this end, analytical processes for deriving the porosity–stress and permeability–stress relationships are discussed in relation to the reservoir rocks. As a result, a porosity–stress formulation is proposed which is in good agreement with the experimental data. Also, two formulas are proposed for permeability–stress relationship; one on the base of Kozeny–Carman permeability–porosity model and the other one is based on a differential form of permeability–porosity relationship. After calibrating the required coefficients for one sandstone and three limestones, it was concluded that the first permeability–stress model is the upper bound correlation while the latter is the lower bound. Furthermore, it is shown that the latter has better agreement to the real experimental data of the sansdstone samples, while the first one is close to the experimental observations from limestone samples. Also, it is concluded that structure of pores is a key factor on permeability–stress relationship, so that there is a significant difference between the experimental data and the proposed relationship for a limestone sample with vuggy pore spaces.     

Least-Squares Wavelet Analysis of Unequally Spaced and Non-stationary Time Series and Its Applications

Least-squares spectral analysis, an alternative to the classical Fourier transform, is a method of analyzing unequally spaced and non-stationary time series in their first and second statistical moments. However, when a time series has components with low or high amplitude and frequency variability over time, it is not appropriate to use either the least-squares spectral analysis or Fourier transform. On the other hand, the classical short-time Fourier transform and the continuous wavelet transform do not consider the covariance matrix associated with a time series nor do they consider trends or datum shifts. Moreover, they are not defined for unequally spaced time series. A new method of analyzing time series, namely, the least-squares wavelet analysis is introduced, which is a natural extension of the least-squares spectral analysis. This method decomposes a time series to the time–frequency domain and obtains its spectrogram. In addition, the probability distribution function of the spectrogram is derived that identifies statistically significant peaks. The least-squares wavelet analysis can analyze any non-stationary and unequally spaced time series with components of low or high amplitude and frequency variability, including datum shifts, trends, and constituents of known forms, by taking into account the covariance matrix associated with the time series. The outstanding performance of the proposed method on synthetic time series and a very long baseline interferometry series is demonstrated, and the results are compared with the weighted wavelet Z-transform.