An experimental study of grain scale melt segregation mechanisms in two common crustal rock types

Creation of pathways for melt to migrate from its source is the necessary first step for transport of magma to the upper crust. To test the role of different dehydration‐melting reactions in the development of permeability during partial melting and deformation in the crust, we experimentally deformed two common crustal rock types. A muscovite‐biotite metapelite and a biotite gneiss were deformed at conditions below, at and above their fluid‐absent solidus. For the metapelite, temperatures ranged between 650 and 800 °C at P_c=700 MPa to investigate the muscovite‐dehydration melting reaction. For the biotite gneiss, temperatures ranged between 850 and 950 °C at P_c=1000 MPa to explore biotite dehydration‐melting under lower crustal conditions. Deformation for both sets of experiments was performed at the same strain rate (ε.) 1.37×10^?5 s^?1. In the presence of deformation, the positive ΔV and associated high dilational strain of the muscovite dehydration‐melting reaction produces an increase in melt pore pressure with partial melting of the metapelite. In contrast, the biotite dehydration‐melting reaction is not associated with a large dilational strain and during deformation and partial melting of the biotite gneiss melt pore pressure builds more gradually. Due to the different rates in pore pressure increase, melt‐enhanced deformation microstructures reflect the different dehydration melting reactions themselves. Permeability development in the two rocks differs because grain boundaries control melt distribution to a greater extent in the gneiss. Muscovite‐dehydration melting may develop melt pathways at low melt fractions due to a larger volume of melt, in comparison with biotite‐dehydration melting, generated at the solidus. This may be a viable physical mechanism in which rapid melt segregation from a metapelitic source rock can occur. Alternatively, the results from the gneiss experiments suggest continual draining of biotite‐derived magma from the lower crust with melt migration paths controlled by structural anisotropies in the protolith.     

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.     

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.     

均匀土-箱基-结构相互作用体系的计算分析

采用通用有限元程序ANSYS，针对捱动台试验中的均匀土-箱基-结构试验进行了三维有限元分析，计算中土体的本构模型采用等效线性模型，利用面-面接触单元考虑土体与基础交界面的状态非线性。计算表明，基础底面和土体发生滑移，基础侧面和土体之间发生了滑移和脱离，上部结构柱顶加速度反应主要由基础转动引起的摆动分量组成，通过与试验结果的对照研究，二者得出的规律基本一致，验证了采用的计算模型与分析方法的合理性，为进一步计算研究和实际工程应用奠定了基础。     

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

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

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

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

Development of site-specific design ground motions in Western and Eastern North America

This paper presents results recently obtained for generating site-specific ground motions needed for design of critical facilities. The general approach followed in developing these ground motions using either deterministic or probabilistic criteria is specification of motions for rock outcrop or very firm soil conditions followed by adjustments for site-specific conditions. Central issues in this process include development of appropriate attenuation relations and their uncertainties, differences in expected motions between Western and Eastern North America, and incorporation of site-specific adjustments that maintain the same hazard level as the control motions, while incorporating uncertainties in local dynamic material properties. For tectonically active regions, such as the Western United States (WUS), sufficient strong motion data exist to constrain empirical attenuation relations for M up to about 7 and for distances greater than about 10–15 km. Motions for larger magnitudes and closer distances are largely driven by extrapolations of empirical relations and uncertainties need to be substantially increased for these cases.

For the Eastern United States (CEUS), due to the paucity of strong motion data for cratonic regions worldwide, estimation of strong ground motions for engineering design is based entirely on calibrated models. The models are usually calibrated and validated in the WUS where sufficient strong motion data are available and then recalibrated for applications to the CEUS. Recalibration generally entails revising parameters based on available CEUS ground motion data as well as indirect inferences through intensity observations. Known differences in model parameters such as crustal structure between WUS and CEUS are generally accommodated as well. These procedures are examined and discussed.     

Seismic analysis of base-isolated liquid storage tanks using the BE–FE–BE coupling technique

A three-dimensional soil–structure–liquid interaction problem is numerically simulated in order to analyze the dynamic behavior of a base-isolated liquid storage tank subjected to seismic ground motion. A dynamic analysis of a liquid storage tank is carried out using a hybrid formulation, which combines the finite shell elements for structures and the boundary elements for liquid and soil. The system is composed of three parts: the liquid–structure interaction part, the soil–foundation interaction part, and the base-isolation part. In the liquid–structure interaction part, the tank structure is modeled using the finite elements and the liquid is modeled using the internal boundary elements, which satisfy the free surface boundary condition. In the soil–foundation interaction part, the foundation is modeled using the finite elements and the half-space soil media are modeled using the external boundary elements, which satisfy the radiation condition in the infinite domain. Finally, above two parts are connected with the base-isolation system to solve the system's behavior. Numerical examples are presented to demonstrate the accuracy of the developed method, and an earthquake response analysis is carried out to demonstrate the applicability of the developed technique. The properties of a real LNG tank located in the west coast of Korea are used. The effects of the ground and the base-isolation system on the behavior of the tank are analyzed.