Dependence of upper atmosphere photochemistry on the shape of the diurnal cycle of the photolysis rates

A photochemical model of the atmosphere constitutes a non–linear, non–autonomous dynamical system, enforced by the Earth’s rotation. Some studies have shown that the region of the mesopause tends towards non–linear responses such as period-doubling cascades and chaos. In these studies, simple approximations for the diurnal variations of the photolysis rates are assumed. The goal of this article is to investigate what happens if the more realistic, calculated photolysis rates are introduced. It is found that, if the usual approximations—sinusoidal and step fiunctions—are assumed, the responses of the system are similar: it converges to a 2–day periodic solution. If the more realistic, calculated diurnal cycle is introduced, a new 4–day subharmonic appear.     

Interpretation of tracer tests performed in fractured rock of the Lange Bramke basin, Germany

 Two multitracer tests performed in one of the major cross-fault zones of the Lange Bramke basin (Harz Mountains, Germany) confirm the dominant role of the fault zone in groundwater flow and solute transport. Tracers having different coefficients of molecular diffusion (deuterium, bromide, uranine, and eosine) yielded breakthrough curves that can only be explained by a model that couples the advective–dispersive transport in the fractures with the molecular diffusion exchange in the matrix. For the scale of the tests (maximum distance of 225 m), an approximation was used in which the influence of adjacent fractures is neglected. That model yielded nearly the same rock and transport parameters for each tracer, which means that the single-fracture approximation is acceptable and that matrix diffusion plays an important role. The hydraulic conductivity of the fault zone obtained from the tracer tests is about 1.5×10^–2 m/s, whereas the regional hydraulic conductivity of the fractured rock mass is about 3×10^–7 m/s, as estimated from the tritium age and the matrix porosity of about 2%. These values show that the hydraulic conductivity along the fault is several orders of magnitude larger than that of the remaining fractured part of the aquifer, which confirms the dominant role of the fault zones as collectors of water and conductors of fast flow. Received, April 1997 Revised, January 1998, August 1998 Accepted, August 1998     

A detailed analysis of the February 1996 aftershock sequence in the eastern Pyrenees, France

Following the 1996 February 18 M _L = 5.2 earthquake in the Agly massif in the eastern French Pyrenees, we installed a temporary network of seismometers around the epicentre. In this paper, we analyse 336 well-located aftershocks recorded from February 19 to February 23 by 18 temporary stations and two permanent stations located less than 35  km from the epicentre. Most aftershocks have been located with an accuracy better than 1.5  km in both horizontal and vertical positions. Their spatial distribution suggests the reactivation of a known fault system. We determined 39 fault-plane solutions using P -wave first motions. Despite their diversity, the focal mechanisms yield an E–W subhorizontal T-axis. We also determined fault-plane solutions and principal stress axes using the method developed by Rivera & Cisternas (1990 ) for the 15 best-recorded events. We obtain a pure-shear-rupture tectonic regime under N–S subhorizontal compression and E–W subhorizontal extension. These principal stress axes, which explain the focal mechanisms for at least 75 per cent of the 39 aftershocks, are different from the axes deduced from the main shock. The post-earthquake stress field caused by the main-shock rupture, modelled as sinistral strike slip on three vertical fault segments, is computed for various orientations and magnitudes of the regional stress field, assumed to be horizontal. The aftershock distribution is best explained for a compressive stress field oriented N30°E. Most aftershocks concentrate where the Coulomb failure stress change increases by more than 0.2  MPa. The diversity of aftershock focal mechanisms, poorly explained by this model, may reflect the great diversity in the orientations of pre-existing fractures in the Agly massif.     

Finite element analysis of flow patterns near geological lenses in hydrodynamic and hydrothermal systems

We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method.
  Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately.
  Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.     

用长周期面波波形拟合及P波初动方向估计中等地震的震源机制

提出利用周期范围４５-１００ｓ间的面波波形拟合与Ｐ波初动符号相结合,估计中小地震震源机制的方法．给出两个倾滑为主和三个走滑为主的震例．其中最小的是１９９６年１２月１６日北京顺义ＭＬ４．５级地震.比较了１９９６年２月５日云南丽江Ｍｓ６．０级余震的机制解与美国哈佛大学ＣＭＴ结果的相似程度,各震例记录图与理论图的相似程度,以及初动分布与主压力轴Ｐ和主动力轴Ｔ分布区的一致程度．从比较结果看,得到的震源机制解有参考价值,开发的软件有可用性．     

Elastic waves in nonhomogeneous media under 2D conditions: II. Numerical implementation

The purpose of this work is to present three methods of analysis for elastic waves propagating in two dimensional, elastic nonhomogeneous media. The first step, common to all methods, is a transformation of the governing equations of motion so that derivatives with respect to the material parameters no longer appear in the differential operator. This procedure, however, restricts analysis to a very specific class of nonhomogeneous media, namely those for which Poisson's ratio is equal to 0.25 and the elastic parameters are quadratic functions of position. Subsequently, fundamental solutions are evaluated by: (i) conformal mapping in conjunction with wave decomposition, which in principle allows for both vertical and lateral heterogeneities; (ii) wave decomposition into pseudo-dilatational and pseudo-rotational components, which results in an Euler-type equation for the transformed solution if medium heterogeneity is a function of one coordinate only; and (iii) Fourier transformation followed by a first order differential equation system solution, where the final step involving inverse transformation from the wavenumber domain is accomplished numerically. Finally, in the companion paper numerical examples serve to illustrate the above methodologies and to delineate their range of applicability.     

横观各向同性圆柱土样轴对称Biot固结的解析解

给出了横观各向同性圆柱土样轴对称Biot固结的封闭级数解,该解由满足非齐次边界条件的弹性力学解与满足齐次边界条件的渗流拟动态解的叠加构成。通过算例分析,探讨了横观各向土体的Biot固结规律和Mandel-Cryer效应。     

Second-Order Solutions for Random Interfacial Waves Under Steady Uniform Currents

1.Introduction Owingtothesimilaritiesbetweenthesurfacewaveandtheinterfacialwave,itisnaturaltoapply themethodsdevelopedforsurfacewavestothestudyofinterfacialwavesasreviewedbyUmeyama (1998;2000).Recently,Song(2004)derivedsecond ordersolutionsforrandominterfacialwavesat aconstantdepthinatwo layerfluidsystemwitharigidlidusinganexpansiontechnique,analogousto thatusedbyLonguet Higgins(1963)andSharmaandDean(1979),tostudyrandomsurfacewaves. Inthispaper,Song’sresultsareextendedtoamoregeneralcaseoft…     

A generalization of extendedQ-mode factor analysis to data matrices with variable row sums

The extended method of Q-mode factor analysis developed by Miesch for data matrices with constant row sums is generalized to data matrices with variable row sums. With the algorithm provided it is possible to compute factor scores in the metric of the original data and compute goodness-of-fit statistics and model geological systems unconstrained by constancy of row sums of data points.