Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Modelling of deformation response and chemo‐mechanical coupling in chalk

A constitutive relation is derived for describing the mechanical response of chalk. The approach is based on a phenomenological framework which employs chemo‐plasticity. The properties of the material are assumed to be affected by the physico‐chemical processes that occur through the interaction between the skeleton and the pore fluid. The underlying mechanism is discussed by invoking a micromechanical analysis. The performance of the framework is illustrated by examining the evolution of mechanical characteristics in the presence of different pore fluids. Copyright © 2006 John Wiley & Sons, Ltd.     

CCD联合平差方法及其初步实例检验

在常规照相天体测量工作中 ,对暗天体的定位是采用逐级定标星过渡的方法实现的。由于星表系统差和局部差的存在 ,最终定位结果可能包含难以把握的误差积累。在CCD小视场观测情况下时常难以找到足够数量的定标星。鉴于此 ,推导了CCD观测联合平差方法的严格矢量表达式 ,以期通过共同星的过渡 ,实现大天区统一平差。原理上此方法可实现相对于河外天体的定位 ,以削弱星表系统差和局部差对定位结果的影响。讨论了CCD联合平差方法实现中的问题 ,包括底片常数模型选择、切点位置改正、共同星的较差改正和法方程解算方法等 ,并分析了具体的处理措施。用实例初步检验了CCD联合平差方法的效果 ,表明通过对多张CCD观测的联合平差 ,可以削弱单张CCD观测参考星数目过少、单张CCD观测参考星分布不均、个别参考星存在位置偏差和局部区域参考星存在系统性位置偏差等不利因素对归算结果的影响 ,进而达到扩大视场和提高归算精度的目的     

Metamorphic evolution of exhumed middle to lower crustal rocks in the Mojave Extensional Belt, southern California, USA

The Waterman Metamorphic Complex of the central Mojave Desert was exposed as a consequence of early Miocene detachment-dominated extension. However, it has evidence consistent with a more extensive geological history that involves collision of a crustal fragment(s), tectonic thickening by overthrusting and two periods of extension. The metamorphic complex contains granitoid intrusives and felsic mylonitic gneisses as well as polymetamorphic rocks that include marble, calc-silicate, quartzite. mafic granulite, pyribolite, amphibolite, migmatite and biotite schist. The latter group of rocks was affected by an initial series of high-grade metamorphic events (M1 and M2) and a localized lower grade overprint (M3). The initial metamorphism (M1) can be separated into two stages along its high-grade P–T path: M1a, a granulite facies metamorphism at 800–850° C and 7.5–9 kbar and Mlb, an upper amphibolite facies overprint at 750–800° C and 10–12 kbar. M1a developed mineral assemblages and textures consistent with granulite facies conditions at a reduced activity of H₂O and is associated with intense ductile deformation (D1) and minor local partial melting. M1b overprinted the granulite assemblages with a series of hydrous phases under conditions of increasing pressure and H₂O activity and is accompanied by little or no deformation. M2 developed at lower pressures and temperatures (650–750° C, 4.5–5.5 kbar) and is distinguished by a second local overprint of hydrous phases that reflects an input of aqueous fluids probably associated with the intrusion of a series of granitic dykes and veins. Effects of M3 are confined to the Mitchel detachment zone, an anastomosing early Miocene detachment fault, and are characterized by local ductile/brittle deformation (D2) of the pre-existing high-grade rocks and granitoid intrusives and by the production of mylonites and mylonitic gneisses under greenschist facies conditions (300–350° C, 3–5 kbar). The initial overprint (M1a) represents metamorphism, devolatilization and minor partial melting of supracrustal rocks under granulite facies conditions as a consequence of tectonic and, possibly, magmatic thickening. The increasing pressure transition of M1a to M1b reflects a period of continued compressional tectonism, thrusting and influx of H₂O, in part, locally related to crystallization of partial melts. The near isothermal decompression between M1b and M2 probably represents a pre-112-Ma extensional episode that may have been the result of a decompressional readjustment of a thickened crust. Following the initial extensional event, the metamorphic complex remained at depths of 10–17 km for at least 90 Ma until it was uplifted following Miocene extension. M3 develops locally in response to this second extensional period resulting from the early Miocene detachment faulting.     

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.     

Timing and tectonic setting of Grenvillian metamorphism—constraints from a transect along Georgian Bay,Ontario*

Systematic mapping of a transect along the well-exposed shores of Georgian Bay, Ontario, combined with the preliminary results of structural analysis, geochronology and metamorphic petrology, places some constraints on the geological setting of high-grade metamorphism in this part of the Central Gneiss Belt. Correlations within and between map units (gneiss associations) have allowed us to recognize five tectonic units that differ in various aspects of their lithology, metamorphic and plutonic history, and structural style. The lowest unit, which forms the footwall to a regional decollement, locally preserves relic pre-Grenvillian granulite facies assemblages reworked under amphibolite facies conditions during the Grenvillian orogeny. Tectonic units above the decollement apparently lack the early granulite facies metamorphism; out-of-sequence thrusting in the south produced a duplex-like structure. Two distinct stages of Grenvillian metamorphism are apparent. The earlier stage (c. 1160–1120 Ma) produced granulite facies assemblages in the Parry Sound domain and upper amphibolite facies assemblages in the Parry Island thrust sheet. The later stage (c. 1040–1020 Ma) involved widespread, dominantly upper amphibolite facies metamorphism within and beneath the duplex. Deformation and metamorphism recently reported from south and east of the Parry Sound domain at c. 1100–1040 Ma have not yet been documented along the Georgian Bay transect. The data suggest that early convergence was followed by a period of crustal thickening in the orogenic core south-east of the transect area, with further advance to the north-west during and after the waning stages of this deformation.     

山东省济宁强磁异常区深部铁矿初步验证及其意义

山东省济宁磁异常是一个重、磁同源体,面积大于100 km2,磁异常峰值为3800nT。钻探验证在孔深1041.57～1796.54m位置发现铁矿体,矿体总厚度74.04～220m,磁性铁平均品位15.89～25.19%。矿石类型有条带状方解磁铁石英岩和条带状磁铁石英大理岩,矿石的主要组成矿物为石英、方解石、磁铁矿、磁赤铁矿、菱铁矿。矿体产于济宁岩群浅变质岩系中,矿床特征与条带状铁建造（BIF）铁矿或鞍山式铁矿有明显区别,铁矿成因类型属与千枚岩、变质中酸性火山岩、大理岩有关的沉积变质型铁矿床。该区铁矿资源潜力巨大。     

中国大别山东南缘首次发现大坝陨坑构造

简介首次在大别山东南缘安徽省境内发现的大坝陨坑构造。经初步评价,知陨坑呈椭圆形,长轴呈北北东向,长约19km,短轴长约12km,最大坑深约2km,是一个有中央隆起区的复杂型陨坑。在卫片上陨坑显示环形影像,地貌形态为一洼地。陨坑基座保存尚好,可对它直观和追索陨坑构造边界。形成于230Ma左右的各类撞击变质岩石系列齐全,其中含有柯石英等典型的撞击变质矿物及撞击碎理等超微构造,特别是在陨坑基座内壁普遍发现有鉴别陨坑构造最可靠的标志——干裂自角砾岩,都证明大坝环形影像是一个典型的陨星撞击坑构造。它的发现,具有很大的科学及经济意义,对今后褶皱山区寻找和研究陨坑构造具有示范和指导作用,同时大大丰富了建设大别山世界地质公园的地质依据。无庸置疑此发现将促进大别山旅游业的发展及陨星撞击科普知识的传播。     

基于Poulos弹性理论的被动桩改进算法

传统的Poulos弹性理论仅适合于均质土中土体侧向位移时桩的性状分析,无法考虑土的层状特性。通过引入层状地基中作用一水平集中力的广义Mindlin解和地面作用有竖向荷载时的应力和位移通解,对Poulos方法进行了改进,使之扩展到多层土中,还用于研究堆载条件下的被动桩变形和受力响应。算例分析表明,改进弹性理论要比Poulos方法更为严密、合理,提高了计算精度,应用范围也更广。     

Modelling of periodic wave transformation in the inner surf zone

In this paper, we analyse the ability of the nonlinear shallow-water (NSW) equations to predict wave distortion and energy dissipation of periodic broken waves in the inner surf zone. This analysis is based on the weak-solution theory for conservative equations. We derive a new one-way model, which applies to the transformation of non-reflective periodic broken waves on gently sloping beaches. This model can be useful to develop breaking-wave parameterizations (in particular broken-wave celerity expression) in both time-averaged wave models and time-dependent Boussinesq-type models. We also derive a new wave set-up equation which provides a simple and explicit relation between wave set-up and energy dissipation. Finally, we compare numerical simulations of both, the NSW model and the simplified one-way model, with spilling wave breaking experiments and we find a good agreement.