Transient‐state groundwater flow in various geometries with wells at arbitrary positions: analytical methods

We have developed a method for analytically solving the porous medium flow equation in many different geometries for horizontal (two‐dimensional), homogeneous and isotropic aquifers containing impermeable boundaries and any number of pumping or injection wells located at arbitrary positions within the system. Solutions and results are presented for rectangular and circular aquifers but the method presented here is easily extendible to many geometries. Results are also presented for systems where constant head boundary conditions can be emulated internal to the aquifer boundary. Recommendations for extensions of the present work are briefly discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

基于土地适宜性评价与元胞自动机模型的城市开发边界划定方法

城市开发边界的划定可以在一定程度上引导城市空间的良性扩张。为了研究边界划定的方法,以我国某港口城市作为研究区域,在空间增长模拟的方向进行探索。研究中借助土地适宜性评价和元胞自动机边界工具,构建既符合生态资源环境要求又符合城市发展需要的边界划定,同时,为了实现城市用地健康和动态管控举措,提出了建立健全土地用地管理机制的要求。     

中国大陆及其邻区强震活动与活动地块关系研究

从活动地块假说出发 ,在活动地块研究的基础上 ,探讨了中国大陆及邻区活动地块与强震活动的关系。研究指出 ,主要构造变形和强烈地震大都发生在活动地块边界。在占总面积 17%的活动地块边界上 ,集中了全部的 8级以上巨大地震和 86 %的 7级以上大地震 ,其释放能量占全部总能量的 95 %以上 ,表明中国大陆及其邻区活动地块边界带控制了绝大部分的强地震。从活动地块的整体来看 ,强震活动不仅显示出显著的韵律性特征 ,而且其高、低起伏基本上与中国大陆地区一致 ,只是强震活跃时段有时稍长于中国大陆。各轮回强震活动都有各自活动的主体地区 ,反映了不同活跃期内地块的不同活动方式。文中还从现今地壳运动角度 ,讨论了活动地块运动速率与强地震活动水平之间的可能联系。     

Two contrasted P–T–time paths of coronitic metanorites of the French Massif Central: are reaction textures reliable guides to metamorphic histories?

Metanorites from two eclogitized metagabbros of the Hercynian French Massif Central preserve coronitic textures of hornblende, garnet, quartz and/or kyanite produced at the expense of the primary magmatic assemblage orthopyroxene and plagioclase. Using a petrogenetic grid in the CFMASH system, two possible P–T evolutions for the origin of the coronas are evaluated. The sequence of reactions involving the formation of Hbl (–Ky) ± Grt and Qtz coronitic assemblages is consistent with an isobaric cooling at high pressure (c. 1–2 GPa) under hydrated conditions. However, this P–T path, inferred by using only petrographical observations, is inconsistent with the geochronological constraints: emplacement of the gabbro at 490 Ma and high‐pressure metamorphism at 410 Ma. In order to reconcile petrographical observations with geochronological constraints, we propose a discontinuous two‐stage evolution involving a change in water activity with time. (1) Emplacement and cooling of the norite at low pressure under anhydrous conditions, at 490 Ma. (2) During the Hercynian orogeny, the norite experienced an increase in pressure and temperature under fluid‐present conditions. Adding water to the system implies a dramatic change in the petrogenetic grid topology, restricting the orthopyroxene–plagioclase assemblage only to high temperatures. Therefore, the breakdown of the unstable magmatic assemblage, through apparent retrograde reactions, occurred along the prograde P–T path which never crossed the equilibrium boundaries of these reactions.     

The migration of fluid-filled grain boundaries in recrystallizing synthetic bischofite: first results of in-situ high-pressure, high-temperature deformation experiments in transmitted light

The effect of fluids on recrystallization behaviour is well known; however, the detailed microscale distribution of fluid in grain boundaries and the influence of fluid on grain boundary migration are still unresolved. In this study, in‐situ deformation experiments in transmitted light microscopy were undertaken, as this allows continuous and direct observation of the whole range of processes involved in fluid‐assisted grain boundary migration. A new see‐through deformation apparatus was developed to enable the control of fluid pressure. Bischofite containing small amounts of aqueous fluid was deformed at temperatures between 50 and 90 °C, over a range of fluid pressure from 0.5 to 1 MPa, and strain rates of 5 × 10^?6 to 1 × 10^?4 s^?1. The rates of grain boundary migration were measured at different temperatures and strain rates. Detailed observations during and after the deformation illustrate the evolution of migrating fluid‐filled grain boundaries and show that the incorporation of fluids from inclusions as well as their pinch‐off is dependent on the grain boundary velocity, the thickness of the grain boundary and the size and shape of the inclusions. Direct evidence is presented for the contraction of the grain boundary fluids into isolated inclusions after equilibrium conditions are attained.     

Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Study of the crust and mantle using magnetic surveys by Magsat and other satellites

Among the first measurements made from near-Earth orbiting satellites were measurements of the magnetic field. The sources of that field lie both within the Earth, in its core and crust, and in the surrounding ionosphere and magnetosphere. This article summarizes some of the methodology and results for studies of the Earth’s mantle and crust. Mantle conductivity studies can be made either by studying signals impressed on the Earth from outside, e.g., the ionosphere or magnetosphere, or by studying signals originating in the core and transmitted through the mantle. Crustal field studies begin with a careful selection of the data and subsequent removal of core and external fields by some sort of filtering. Average maps from different local times sometimes differ, presumably due to the remaining presence of fields of external origin. Several techniques for further filtering are discussed. Where large-area aeromagnetic maps are available, crustal maps derived from satellite data can be compared with upward continued data. In general, the comparisons show agreement, with some differences, particularly in and near the auroral belts. The satellite data are further reduced by various methods of inverse and forward modelling, sometimes including reduction to the pole (RTP). These techniques are generally unstable at the equator. Common methods of stabilizing the inversions include principle components analysis and ridge regression. Because of the presence of the core field, the entire crustal contribution from the field is not known. Also, there is a basic nonuniqueness to the inverse solutions. Nevertheless, magnetizations that are interpretable can be derived.     

Nonequilibrium thermodynamics of pressure solution

This paper is concerned with the thermodynamic theory of solution and precipitation processes in wet crustal rocks and with the mechanism of steady pressure-solution slip in contact zones, such as grain-to-grain contacts, fracture surfaces, and permeable gouge layers, that are infiltrated by a mobile aqueous solution phase. A local dissipation jump condition at the phase boundary is fundamental to identifying the thermodynamic force driving the solution and precipitation process and is used here in setting up linear phenomenological relations to model near-equilibrium phase transformation kinetics. The local thermodynamic equilibrium of a stressed pure solid in contact with its melt or solution phase is governed by Gibbs's relation, which is rederived here, in a manner emphasizing its independence of constitutive assumptions for the solid while neglecting surface tension and diffusion in the solid. Fluid-infiltrated contact zones, such as those formed by rough surfaces, cannot generally be in thermodynamic equilibrium, especially during an ongoing process of pressure-solution slip, and the existing equilibrium formulations are incorrect in overlooking dissipative processes tending to eliminate fluctuations in superficial free energies due to stress concentrations near asperities, defects, or impurities. Steady pressure-solution slip is likely to exhibit a nonlinear dependence of slip rate on shear stress and effective normal stress, due to a dependence of the contact-zone state on the latter. Given that this dependence is negligible within some range, linear relations for pressure-solution slip can be derived for the limiting cases of diffusion-controlled and interface-reaction-controlled rates. A criterion for rate control by one of these mechanisms is set by the magnitude of the dimensionless quantityk/2C _pD, wherek is the interfacial transfer coefficient, is the mean diffusion path length,C _p is the solubility at pressurep, andD is the mass diffusivity.