A non-linear stress-stiffness model for geomaterials at small to intermediate strains

A double exponential fitting model (DEFM) capable of expressing the non-linear stress-stiffness relationship of geomaterials has been proposed by Shibuya et al. (1997). The model comprises two material constants; the elastic stiffness at very small strains and the strength, together with other free parameters to determine the complete stress-stiffness relationship. In this paper, the capability of the original function used for DEFM in simulating the tangent stiffness-stress relationship of geomaterials is first discussed. Second, the methods for determining the free model parameters, as well as its conversion to obtain a stress-strain relationship are proposed. The applicability of DEFM to predicting non-linear stress-stiffness relationship is examined in detail in a total of forty-nine fitting cases of compression test data on sedimentary rock, artificial soft rock and soft clay. It is found that the DEFM is effective in expressing the non-linear stress-stiffness relationship of various kinds of geomaterials at small to intermediate strains, say less than 0.5%. The superiority of this model compared to other fitting models currently in use is also demonstrated in some of the fitting cases.     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Origin of Gold—Bearing Fluid and Its Initiative Localization Mechanism in Xiadian Gold Deposit,Shandong Province

The composition of quartz inclusions and trace elements in ore indicate that gold-bearing fluid in the Xiadian gold deposit,Shandong Province,stemmed from both mantle and magma,belonging to a composite origin.Based on theoretical analysis and high temperature and high pressure experimental studies,gold-bearing fluid initiative localization mechanism and the forming environment of ore-host rocks are discussed in the present paper.The composite fluid extracted gold from rocks because of its expanding and injecting forces and injecting forces and flew through ore-conducive structures,leading to the breakup of rocks.The generation of ore-host faults and the precipitation of gold-bearing fluid occurred almost simultaneously.This study provides fur-ther information about the relationships between gold ore veins and basic-ultrabasic vein rocks and intermediate vein rocks,the spatial distribution of gold ore veins and the rules governing the migration of ore fluids.     

县市级政府决策信息系统的开发应用:专业GIS 软件SK set V4 .0功能介绍———以贵州沿河土家族自治县资源遥感信息系统为例

SK setV4.0是中国地质科学院岩溶地质研究所GIS中心2001年4月开发成功的、具有自主板权的专业GIS平台软件。它集成了现今较为先进的组件、控件重用，面向对象编程等技术，并以美国Inprise公司的BorladnC^ Builder编译软件为开发平台，提供了全中文菜单环境，不仅具有基本信息和多图层查询、相关分析动态监测、多因子分析、图斑编辑、自动演示、系统维护等功能，而且还具有快速响应的工具条、友好的对话框和较强的纠错机制，很适合缺乏专业GIS人才的县市级政府管理部门使用。文中还以贵州省沿河县资源遥感信息系统为例，对SK setV4.0的功能和使用方法进行了具体的介绍。     

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.     

The Quaternary evolution of the Gulf of Corinth, central Greece: coupling between surface processes and flow in the lower continental crust

The Gulf of Corinth in central Greece is an active normal fault zone with particularly clear evidence of isostatic footwall uplift, constrained by Quaternary marine terraces, and hanging-wall subsidence and sedimentation. It is bounded to the south by a Pliocene to Early Pleistocene sedimentary basin, which is now eroding into the Gulf. Previous work has suggested that the relief across this region has increased dramatically since the Early Pleistocene, due to the isostatic response to increased rates of footwall erosion and hanging-wall sedimentation. It is indeed assumed here that incision accompanying the draw-down of global sea-level at 0.9 Ma, during the first major Pleistocene glaciation, initiated the erosion of the basin south of the Gulf of Corinth and so abruptly increased the sedimentation rate in the Gulf. The resulting transient thermal and isostatic response to these changes is modelled, with the subsiding depocentre and eroding sediment source coupled by flow in the lower continental crust. The subsequent enhancement of relief, involving an increase in bathymetry from near zero to 900 m and 500 m of uplift of the eroding land surface in the sediment source, is shown to be a direct consequence of this change. The model is sensitive to the effective viscosity of the lower crust, and can thus resolve this parameter by matching observations. A value of 6×10¹⁹ Pa s is indicated, suggesting a viscosity at the Moho no greater than 10¹⁸ Pa s. Similar transient topographic effects caused by increased rates of sedimentation and erosion are likely to be widespread within the geological record, suggesting that this coupling process involving flow in the weak lower crust may be of major geological and geomorphological importance.     

Physical properties of ultrahigh-pressure metamorphic rocks from the Sulu terrain, eastern central China: implications for the seismic structure at the Donghai (CCSD) drilling site

About 30 samples representing major lithologies of Sulu ultrahigh-pressure (UHP) metamorphic rocks were collected from surface exposures and exploration wells, and compressional (Vp) and shear wave (Vs) velocities and their directional dependence (anisotropy) were determined over a range of constant confining pressures up to 600 MPa and temperatures ranging from 20 to 600 °C. Samples range in composition from acidic to ultramafic. P- and S-wave velocities measured at 600 MPa vary from 5.08 to 8.64 km/s and 2.34 to 4.93 km/s, respectively. Densities are in the range from 2.60 to 3.68 g/cm³. To make a direct tie between seismic measurements (refraction and reflection) and subsurface lithologies, the experimental velocity data (corresponding to shallow depths) were used to calculate velocity profiles for the different lithologies and profiles of reflection coefficients at possible lithologic interfaces across the projected 5000-m Chinese Continental Scientific Drilling Program (CCSD) crustal segment. Comparison of calculated in situ velocities with respective intrinsic velocities suggests that the in situ velocities at shallow depths are lowered by an increased abundance of open microcracks. The strongly reflective zone beneath the Donghai drill site can be explained by the impedance contrasts between the different lithologies. Contacts between eclogite/peridotite and felsic rocks (gt-gneiss, granitic gneiss), in particular, may give rise to strong seismic reflections. In addition, shear-induced (lattice preferred orientation (LPO)-related) seismic anisotropy can increase reflectivity. For the explanation of the high velocity bodies (>6.4 km/s) around 1000 m and below 3200-m depth, large proportions of eclogite/peridotite (about 40 and 30 vol.%, respectively) are needed.     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.     

Al hematites prepared by homogeneous precipitation of oxinates: material characterization and determination of the Morin transition

A novel method for synthesis of aluminium hematites, based upon the homogeneous precipitation of Fe and Al oxinates in various proportions, is presented. The precursor precipitates are heated in air at 700?°C. X-ray diffraction, thermal analyses, BET, FTIR, optical reflection analysis, TEM and Mössbauer spectroscopy at room temperature and 80?K of the resulting products indicate that single-phase hematites are formed with structural Al substitution of up to 10 at%. Interestingly, the particle size (>100?nm) is not substantially reduced by the Al content. Although it remains difficult to obtain a homogeneously distributed Al substitution in the final hematite, this processing line offers a unique opportunity to separate the effects of grain size and Al substitution on the Morin transition temperature (T _M) of Al hematite. From the comparison between the present hematites and a series of Al-substituted hematites with lepidocrocite as precursor, it could be shown that the effect on T _M, associated with a change of a factor 10 in grain size, is about 1/3 of the effect caused by a change of 10 in the degree of substitution. Finally, it is suggested that proper thermal treatments under different conditions of the same precursors are likely to produce spinel phases.