An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

Spectroscopic evidence for pressure-induced phase transitions in diopside

 Raman spectra of diopside were collected from atmospheric pressure to 71 GPa. The pressure dependences of 22 modes were determined. Changes occurred in the spectra at three different pressures. First, at approximately 10 GPa, the two Raman modes at 356 and 875 cm⁻¹ disappeared, while the mode at 324 cm⁻¹ split into two modes, diverging at this pressure with significantly different pressure shifts; second, at approximately 15 GPa, a small (1 to 2 cm⁻¹) drop in several of the frequencies was observed accompanied by changes in the pressure dependency of some of the modes; and third, above 55 GPa, the modes characteristic of chains of tetrahedrally coordinated silicon disappeared, while those for octahedrally coordinated silicon appeared. The first change at 10 GPa appears to be a C2/c to C2/c transition involving a change in the Ca coordination. The third change above 55 GPa appears to be a change in the silicon coordination. At 15 GPa, it is suggested that a change in compressional mechanism takes place. Received: 14 November 2000 / Accepted: 9 January 2002     

Fluid permeability of sedimentary rocks in a complete stress–strain process

The permeability of sedimentary rocks during triaxial compression tests was investigated to relate it to the complete strain–stress process. It was found that the permeability was not constant, but varied with the stress and strain states in the rocks. Prior to the peak strength, the permeability decreases with increasing load. A dramatic increase in permeability occurs during the strain softening period. In the present study, in situ measurements of fluid flow and pressure in floor strata was carried out in a double longwall mining face in the Yangzhuang colliery. These measurements show that both the strata pressure and the position with respect to the mining face influence the hydrogeologic properties. The permeability increased in the floor strata behind the mining face because those mining induced fractures opened as the strata pressure decreased. To better understand this change in hydraulic behavior around the mining faces, 3-D numerical modeling was carried out. The model provides the general picture of the stress distribution and failure zone both in the floor and roof strata. The field and model results demonstrate the importance of changes in the stress and strain states on the hydrogeology of a site.     

Three Alkali-Rich Intrusive Rock Belts Newly Discovered in the Mt. Kunlun-Mt. Altun Region, Xinjiang, China

The newly discovered three alkali-rich intrusive rock belts in the Mt. Kunlun Mt.Altun region of southern Xinjiang are the Lapeiquan-Yitunbulak alkali-rich intrusive rock belt,the Gez-Taxkorgan alkali-rich intrusive rock belt and the Beilisai-Abulash alkali-rich intrusive rock belt. The former two belts were formed during the Yanshanian period, and the third one was formed during the Himalayan period, which is the youngest alkali-rich intrusive rock belt in China. The discovery of the alkali-rich intrusive rock belts is of great significance in shedding light on the history of tectono-magmatic activities in this region.     

Aeromagnetic signatures over western Marie Byrd Land provide insight into magmatic arc basement, mafic magmatism and structure of the Eastern Ross Sea Rift flank

Aeromagnetic signatures over the Edward VII Peninsula (E7) provide new insight into the largely ice-covered and unexplored eastern flank of the Ross Sea Rift (RSR). Positive anomalies, 10–40 km in wavelength and with amplitudes ranging from 50 to 500 nT could reveal buried Late Devonian(?)–Early Carboniferous Ford Granodiorite plutons. This is suggested by similar magnetic signature over exposed, coeval Admiralty Intrusives of the Transantarctic Mountains (TAM). Geochemical data from mid-Cretaceous Byrd Coast Granite, contact metamorphic effects on Swanson Formation and hornblende-bearing granitoid dredge samples strengthen this magnetic interpretation, making alternative explanations less probable. These magnetic anomalies over formerly adjacent TAM and western Marie Byrd Land (wMBL) terranes resemble signatures typically observed over magnetite-rich magmatic arc plutons. Shorter wavelength (5 km) 150 nT anomalies could speculatively mark mid-Cretaceous mafic dikes of the E7, similar to those exposed over the adjacent Ford Ranges. Anomalies with amplitudes of 100–360 nT over the Sulzberger Bay and at the margin of the Sulzberger Ice Shelf likely reveal mafic Late Cenozoic(?) volcanic rocks emplaced along linear rift fabric trends. Buried volcanic rock at the margin of the interpreted half-graben-like “Sulzberger Ice Shelf Block” is modelled in the Kizer Island area. The volcanic rock is marked by a coincident positive Bouguer gravity anomaly. Late Cenozoic volcanic rocks over the TAM, in the RSR, and beneath the West Antarctic Ice Sheet exhibit comparable magnetic anomaly signature reflecting regional West Antarctic Rift fabric. Interpreted mafic magmatism of the E7 is likely related to mid-Cretaceous and Late Cenozoic regional crustal extension and possible mantle plume activity over wMBL. Magnetic lineaments of the E7 are enhanced in maximum horizontal gradient of pseudo-gravity, vertical derivative and 3D Euler Deconvolution maps. Apparent vertical offsets in magnetic basement at the location of the lineaments and spatially associated mafic dikes and volcanic rocks result from 2.5D magnetic modelling. A rift-related fault origin for the magnetic lineaments, segmenting the E7 region into horst and graben blocks, is proposed by comparison with offshore seismic reflection, marine gravity, on-land gravity, radio-echo sounding, apatite fission track data and structural geology. The NNW magnetic lineament, which we interpret to mark the eastern RSR shoulder, forms the western margin of the “Alexandra Mountains horst”. This fundamental aeromagnetic feature lies on strike with the Colbeck Trough, a prominent NNW half-graben linked to Late Cretaceous(?) and Cenozoic(?) faulting in the eastern RSR. East–west and north–north–east to NE magnetic trends are also imaged. Magnetic trends, if interpreted as reflecting the signature of rift-related normal faults, would imply N–S to NE crustal extension followed by later northwest–southeast directed extension. NW–SE extension would be compatible with Cenozoic(?) oblique RSR rifting. Previous structural data from the Ford Ranges have, however, been interpreted to indicate that both Cretaceous and Cenozoic extensions were N–S to NE–SW directed.     

淤积泥砂对水平地运动作用时刚性坝面动压力的影响

文章基于固体介质，流体介质，液固两相饱和介质的时域显示有限元波动分析方法，研究了可压缩库水条件下淤积泥砂对水平地运动作用时刚性坝面上动压力的影响，研究中分析了作为两相介质处理的淤积泥砂的饱和度、厚度、渗透系数和孔隙率对坝面动压力的影响，还比较了作为流程两相介质，单相固体介质、单相重流体介质处理的不同淤积泥砂层模型间的计算结果差异。     

Spatial and temporal variability of pore‐water pressures in residual soil slopes in a tropical climate

It is critical to understand and quantify the temporal and spatial variability in hillslope hydrological data in order to advance hillslope hydrological studies, evaluate distributed parameter hydrological models, analyse variability in hydrological response of slopes and design efficient field data sampling networks. The spatial and temporal variability of field‐measured pore‐water pressures in three residual soil slopes in Singapore was investigated using geostatistical methods. Parameters of the semivariograms, namely the range, sill and nugget effect, revealed interesting insights into the spatial structure of the temporal situation of pore‐water pressures in the slopes. While informative, mean estimates have been shown to be inadequate for modelling purposes, indicator semivariograms together with mean prediction by kriging provide a better form of model input. Results also indicate that significant temporal and spatial variability in pore‐water pressures exists in the slope profile and thereby induces variability in hydrological response of the slope. Spatial and temporal variability in pore‐water pressure decreases with increasing soil depth. The variability decreases during wet conditions as the slope approaches near saturation and the variability increases with high matric suction development following rainfall periods. Variability in pore‐water pressures is greatest at shallow depths and near the slope crest and is strongly influenced by the combined action of microclimate, vegetation and soil properties. Copyright © 2002 John Wiley & Sons, Ltd.     

含预制软弱带的岩石破裂过程的数值模拟

考虑到岩石脆性破坏过程中介质的不均匀性之特点，对岩石样品中含预制弱介质条带的岩石样品破坏过程进行二维有限元数值模拟，并对弱介质带的破坏贯通过程、新的断层的产生和有关的地震活动进行了研究．数值模型展示了岩石从变形、微观破坏到整体破坏的全过程以及微震活动的时空分布特征．应力、应变和微震活动的时空分布形象地描绘了岩石变形的局部化和时空迁移等现象，这与实际地壳中所观测到的现象是一致的．此外，模拟结果与实验观测结果也是一致或相似的．      

地震前兆时空非均匀性指标Cv值的实验检验

在双向加载条件下，对两种构造和介质不同的岩石标本进行变形实验，测量了标本表面不同部位的应变．在此基础上，利用王晓青和陈学忠等提出的描述地震前兆非均匀分布的参量————Cv值，分析了标本变形失稳过程中应变异常的分布，以期从实验角度检验Cv值方法并探讨其物理意义．研究表明，Cv值的变化与岩石变形特征的变化有关，是描述前兆分布非均匀性的一种有效指标；Cv值在失稳发生前先异常上升，然后下降，并在Cv值恢复过程中或恢复后伴随着失稳事件的发生．因此，Cv值是一种有效的预报指标．