Dissolution and depletion of ferromagnesian minerals from Holocene tephra layers in an acid bog,New Zealand,and implications for tephra correlation

This study examines the depletion of ferromagnesian silicate minerals from a sequence of thin, distal, mainly rhyolitic tephra layers of Holocene age preserved in an acid peat bog (Kopouatai), North Island, New Zealand. The rate of such depletion has been fast, as indicated by the complete loss of biotite from one tephra layer (Kaharoa Tephra), in which it is normally dominant, in only ca. 770 yr. Chemical dissolution is advocated as the likely cause for the depletion, with amphiboles and other mineral grains commonly showing etch pits, microcaves, and other characteristic surface solution features. Theoretical thermodynamic and kinetic models show a marked increase in the rate of dissolution of all ferromagnesian minerals under conditions of low pH (< 4), but that where silica concentrations in solution are high the relative proportions of minerals remaining are unaffected. However, where concentrations of dissolved silica are low, as in most bog environments, the relative proportions of ferromagnesian minerals are affected as well as absolute amounts being decreased. Amphiboles are depleted relative to pyroxenes, consistent with kinetic studies. The results show that the identification and correlation of tephras on the basis of relative abundances of ferromagnesian minerals alone may be unreliable, and emphasise the need to use multiple criteria in such studies.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Generation of random waves in time-dependent extended mild-slope equations using a source function method

We develop techniques of numerical wave generation in the time-dependent extended mild-slope equations of Suh et al. [1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Engineering 32, 91–117] and Lee et al. [2003. Extended mild-slope equation for random waves. Coastal Engineering 48, 277–287] for random waves using a source function method. Numerical results for both regular and irregular waves in one and two horizontal dimensions show that the wave heights and the frequency spectra are properly reproduced. The waves that pass through the wave generation region do not cause any numerical disturbances, showing usefulness of the source function method in avoiding re-reflection problems at the offshore boundary.     

南极布兰斯菲尔德海峡表层沉积物的酸解烃

本文研究了南极布兰斯菲尔德海峡表层沉积物的酸解烃分布特征，并佐以热释汞和蚀变碳酸盐的分析。从酸解烃中甲烷含量的低值和Ｃ１／（Ｃ２＋Ｃ３）比值，以及乙烯的存在等判断，布兰斯菲尔德海峡沉积物的烃类气体来自于生物成因。     

Parabolic Wave Propagation Modelling in Orthogonal Coordinate Systems

This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild-slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the confonnal system. The computational results show that the confonnal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical confonnal grid system.     

低硼氢氧化镁的制备及热分解动力学参数

本文以水合氯化镁、工业氨水为原料制备出了低硼含量的纯度 >99%的氢氧化镁和氧化镁。利用差热 -热重 (DTA- TG)、X-射线衍射 (XRD)、透射电子显微镜 (TEM)等手段 ,对氢氧化镁的热稳定性、相结构、表观状态以及脱水动力学参数进行了分析。结果为 :氢氧化镁以束状纤维和片状粉末共存 ,在 6 2 3K基本完成脱水反应生成氧化镁 ,其脱水反应的表观活化能 E为135 .4 k J· mol-1,指前因子 A为 5 .15× 10 10 ,反应级数 n为 1.10     

黄海海域陆相中生界地震反射特征及靶区优选

利用前人对黄海海域油气勘探成果，作者对北黄海盆地，南黄海盆地北部坳陷陆相中生界进行了研究。结果认为：黄海海域陆相中生界分布广、厚度大、具有良好的生、储、盖组合，是该区下步油气战略性调查的首选目标。并提出北黄海陆相中生界油气勘探应以东部凹陷为主，同时兼顾中部凹陷。南黄海北部坳陷近期油气勘探重点应选择北部凹陷和加大东北凹陷的研究力度。南黄海南部坳陷陆相中生界油气勘探重点放在南五和南七凹陷。     

一种带自由面的三维船体贴体网络的数值生成

本文提出了一个生成自由面的三维船体贴体网格的数值方法。文中以静水面以上的船体型线作为描述非定常自由面的计算网格域，进一步完善了我们以往开发的无法考虑兴波与静水面以上的船体型线相互作用对船舶粘性流动影响的网格生成方法。静水面以上的初始自由面网格分别沿船面法线η及切线ζ方向的曲线长的指数进行分布，即可极方便地调节自由面在船面及静水面附近的网格间距，也保证了在数值离散三维Ｐｏｉｓｓｏｎ方程时具有相当好的收敛性。这种方法的特点是计算量相当小，尤其能极方便地控制网格在船面簿粘性层中及静水面附近的分布。以Ｓｅｒｉｅｓ６０船模作为计算算例，带自由面的网格分布是相当满意的