Lateral dynamic compliance of pile embedded in poroelastic half space

The time-harmonic response of a pile in a poroelastic half space and under lateral loadings is studied. By treating the pile as a one-dimensional structure and the half-space as a three-dimensional poroelastic continuum, the dynamic interaction between a pile and a poroelastic medium is formulated as a Fredholm integral equation of the second kind. Green's functions for a distributed lateral force field acting inside a poroelastic half space is an important ingredient of this paper. Numerical results for lateral dynamic compliance functions are presented to illustrate the dynamic characteristics of a pile in a poroelastic half space.     

Radiating Vortices in Geophysical Fluid Dynamics

The dynamics of a single vortex on a beta-plane is discussed in this paper. A barotropic, an equivalent barotropic, one-and-a half and two-layer models are considered. The momentum and energy balances are used to describe the evolution of a vortex. A quasi-stationary balance of the Rossby, Zhukovsky-Kutta forces and the force induced by Rossby-wave radiation, describes the dynamics of the barotropic vortex. A net Coriolis force occurs if the fluid is stratified. The difference between the dynamics of cyclones and anticyclones results directly from the Coriolis force acting on a single vortex in a stratified fluid.All vortices radiate Rossby waves in the quasigeostrophic approximation but intense anticyclones propagate steadily in a one-and-a half layer model. A critical amplitude that bounds radiating and steadily propagating anticyclones is found. Steady propagation of anticyclones in general is impossible in a two-layer fluid due to the radiation of a barotropic Rossby-wave. Some solutions of solitary wave type which are known for a two-layer model, survive owing to wave interference.A single vortex can extract energy from a Rossby wave if synchronism conditions are satisfied. The wave interference again plays a crucial role in this case. The wave interference also determines the energy exchange of vortices located at larger distances. If the distance between the vortices is shorter than the length of the radiated waves, modon may be formed due to a small energy loss.The unbounded monotonic variation of the planetary vorticity is a characteristic feature of a beta-plane approximation. As a result, a single vortex propagates up to a 'rest latitude' where it disappears. The evolution of a single barotropic vortex over bottom topography provides another example of a background vorticity distribution with a local extremum above hills (valleys) or ridges (troughs). Physics of its movement differs from a beta-plane case, but if a vortex lies over broad topography, equations are similar and the evolution of a vortex manifests the same typical features. Particularly, a cyclonic vortex tends to drift to the top of a hill or a ridge. An anticyclonic vortex, on the contrary, slides to the bottom of a valley or a trough.An interaction of a barotropic vortex with a broad mean flow is tractable qualitatively on the basis of previous results. Numerical examples illustrating absorption of a small vortex by a larger one and a vortex movement across the flow, are direct analogies of the vortex evolution over a hill and a ridge, respectively. At the same time, strong influence of strain drastically changes the vortex structure.     

Impact of sensor measurement error on sensor positioning in water quality monitoring networks

This paper studies the impact of sensor measurement error on designing a water quality monitoring network for a river system, and shows that robust sensor locations can be obtained when an optimization algorithm is combined with a statistical process control (SPC) method. Specifically, we develop a possible probabilistic model of sensor measurement error and the measurement error model is embedded into a simulation model of a river system. An optimization algorithm is used to find the optimal sensor locations that minimize the expected time until a spill detection in the presence of a constraint on the probability of detecting a spill. The experimental results show that the optimal sensor locations are highly sensitive to the variability of measurement error and false alarm rates are often unacceptably high. An SPC method is useful in finding thresholds that guarantee a false alarm rate no more than a pre-specified target level, and an optimization algorithm combined with the thresholds finds a robust sensor network.     

Late Quaternary Activity and Paleoearthquakes along the Nanyukou Segment of the Northern Piedmont Fault of Wutai Mountain

The northern piedmont fault of Wutai Mountain is located at the north of the Shanxi Graben system,which is the dominating fault of the south boundary of the Fanshi-Daixian depression. This paper discusses the fault activity and paleoearthquakes around the Nanyukou segment of the northern piedmont fault of Wutai Mountain during the late Quaternary through field investigation along the fault,measuring geomorphic deformation and excavating trenches at some important sites. From Nanyukou to the southwest of Shanhui,we find obviously dislocated alluvial fans,with strong neotectonic movement at these sites. Since nearly 20ka,the vertical average slip rate is 1. 55mm /a to 2. 0mm /a. However,since nearly 6ka,it has reached as high as 2. 3mm /a,which is twice that on other segments. 2 trenches were excavated around Nanyukou with 6 events discovered. The referenced ages of the events are before 7600a,6700a ～ 7600a,5321a ～ 5575a,4400a ～ 5400a,4200a ～ 4400a and after 1600a B. P. with approximate recurrence interval 1400a. The latest event is likely to be the earthquake occurring at 512 A. D. ,so it is necessary to do further work to verify this in the future.     

INTERPRETATION AND DEPTH OF INVESTIGATION OF GRADIENT MEASUREMENTS IN DIRECT CURRENT GEOELECTRICS*

Gradient measurements in a homogeneous electrical primary field can easily be interpreted for simple models. The simplified solution (conducting or resistant body in a homogeneous space in a homogeneous electrical field) is often sufficiently accurate, as comparisons with the exact solution (body of finite resistivity in a homogeneous half-space in a quasihomogeneous electrical field) show. The exact geometry of the body cannot be determined by gradient measurements; the same anomaly of apparent resistivity can be caused by different bodies. In particular, the similarity between a sphere and a cube of the same volume is very high. There is a distinct influence of the resistivity of the overburden: the higher this resistivity is, the stronger is the effect caused by a buried body. If a deviation of 10% of the apparent resistivity is assumed as the lower boundary at which a buried body can be detected by gradient measurements, the depth of investigation for a three-dimensional body is approximately equal to its width; in the two-dimensional case the thickness of the overburden can be twice the width. If the overburden has a resistivity which is higher than the resistivity of the substratum, these depths are greater. The greatest possible depth is approximately three times the width of the body.     

Boundary shock waves in a conducting gas

Summary In this paper, we have postulated, in a conducting gas, the occurrence of a boundary shock wave, considered as a quick transition region in which the viscous effects are compressive and are confined to a thin layer adjacent to a surface. The analogues of the well known Rankine-Hugoniot jump relations, Prandtl relation, and other properties of a shock wave are derived and a discussion about the existence of a boundary shock wave is made.     

Some consequences of the differences in seismic properties in old and young continental plates for isostatic compensation

We present a short review of seismological data and show that a difference of upper mantle structure is clearly resolved between old and young continental plates. An Archean plate is approximately 1 second faster than a Paleozoic plate for P waves. Such a change in plate structure implies that the density of an Archean lithosphere should be greater than the density of a Paleozoic lithosphere, if they have the same chemical composition. Applying isostasy, we would expect a Precambrian plate to sink and its surface to stand at a level 8 km beneath that of a Paleozoic plate. We first check that the weight of a crustal column is not smaller for a Precambrian crust than for a Paleozoic crust. Explosion seismology data show that isostatic equilibrium is nearly achieved in the crust. Therefore a compensating mechanism has to be acting in the subcrustal lithosphere. Two hypotheses are presented. Provided that the chemical composition of the lower lithosphere does not vary with time, a simple condition for a continental plate not to sink when aging, would be that the density of the lower lithosphere be ?0.1 g/cm³ smaller than the density of the asthenosphere. Another condition that would prevent plates to sink would be a lowering of the density of the lower lithosphere with age due for example to a variation of its chemical composition and to a reduction in the mean atomic weight. Both hypotheses imply that isostatic equilibrium of old continental plates is maintained by a geochemical mechanism. Although thermal evolution is a key factor for explaining the evolution of continents, it has to be complemented by a hypothesis on the change of the chemistry of the lower continental lithosphere with respect to the asthenosphere and/or with respect to age.     

Contamination of a well in a uniform background flow

In this paper we describe the transport of pollution in groundwater in the neighbourhood of a well in a uniform background flow. We compute the rate at which contaminated particles reach the well as a function of the place of the source of pollution. The motion of a particle in a dispersive flow is seen as a random walk process. The Fokker-Planck equation for the random motion of a particle is transformed using the complex potential for the advective flow field. The resulting equation is solved asymptotically after a stretching transformation. Finally, the analytical solution is compared with results from Monte Carlo simulations with the random walk model. The method can be extended to arbitrary flow fields. Then by a numerical coordinate transformation the analytical results can still be employed.     

西宁地区第四纪晚期孢粉组合及古环境

西宁地区位于黄土高原与青藏高原的接壤处,是环境和气候变迁的敏感地带。本文对西宁市北川河河谷地区黄土沉积的钻孔剖面进行了孢粉分析,获得了比较丰富的孢粉数据,结果呈现出以松、桦、蒿及藜等木本及草本植物花粉为主的孢粉组合。据各类植物花粉量在该剖面垂向上的明显变化,由下而上可划分为5个孢粉带。研究表明,该剖面的孢粉分析数据及资料不仅记录了75kaB.P以来大量的古地理及古环境的变迁,也揭示了西宁地区晚更新世以来古植被及古气候的演替及变化,其演变分别经历了5个阶段:第一阶段为森林型针叶林,反应了温爽轻湿的气候;第二阶段为针叶林,为温凉轻润的气候;第三阶段为针叶林夹稀树草原,属温和湿润夹温凉轻湿的气候;第四阶段为森林草原,反应了温凉轻爽的气候;第五阶段为针叶林,反应了温和湿润的气候。     

SEQUENTIAL STIFFNESS DESIGN FOR SEISMIC DRIFT RANGES OF A SHEAR BUILDING–PILE–SOIL SYSTEM

A shear building supported by a prescribed pile–soil system is subjected to bedrock earthquake input. A new design procedure is presented for generating a sequence of stiffness designs satisfying the constraints on interstorey drifts. The mean peak interstorey drifts of the shear building subjected to a set of spectrum-compatible ground motions at the bedrock are evaluated by a modal combination rule. Tuning of the fundamental natural period of a shear building with a fixed base with that of a shear beam ground results in a non-monotonic sequence of stiffness designs with respect to a ground stiffness parameter and previous approaches cannot be applied to such a problem. This difficulty in finding such a non-monotonic sequence is overcome by utilizing the ground stiffness parameter and the superstructure stiffness parameter alternately in multiple design phases and by developing a new multi-phase perturbation technique. Fundamental characteristics of this sequence of stiffness designs and the effect of ground stiffnesses on the design of the shear building are disclosed. It is further shown that the stiffness contour method is also useful for the design procedure such that a scattering effect in the estimates of ground stiffnesses is taken into account. The usefulness of the proposed procedure of sequential stiffness design and contour line method is demonstrated through several sequential design examples.     

A stochastic methodology for risk assessment of a large earthquake when a long time has elapsed

We propose a stochastic methodology for risk assessment of a large earthquake when a long time has elapsed from the last large seismic event. We state an approximate probability distribution for the occurrence time of the next large earthquake, by knowing that the last large seismic event occurred a long time ago. We prove that, under reasonable conditions, such a distribution is exponential with a rate depending on the asymptotic slope of the cumulative intensity function corresponding to a nonhomogeneous Poisson process. As it is not possible to obtain an empirical cumulative distribution function of the waiting time for the next large earthquake, an estimator of its cumulative distribution function based on existing data is derived. We conduct a simulation study for detecting scenario in which the proposed methodology would perform well. Finally, a real-world data analysis is carried out to illustrate its potential applications, including a homogeneity test for the times between earthquakes.     

Analytical configurations of a force-free magnetic cylinder in the solar wind

The problem of a smooth field configuration, which should be an initial configuration in modeling (using the method of coarse particles) the problem of a stationary solar wind flow around a magnetic cloud in the case of a spatially two-dimensional statement (when a magnetic cloud is considered as a force-free magnetic cylinder with a finite radius) is considered. It has been indicated that such a statement is possible only when the magnetic field in the solar wind is parallel to the cylinder axis. The method for finding the magnetic field of a force-free cylinder with a finite radius, when some field component is specified and another component is determined based on this one (which makes it possible to construct fields with preassigned properties), has been proposed. The variant for constructing the initial field configuration in the transition region around a cylinder has been proposed. This variant makes it possible to gradually pass from homogeneous crossed fields in the solar wind to a force-free magnetic and zero electric fields within a cylinder, an electric field being potential and orthogonal to a magnetic field (in the reference system related to a magnetic cloud).     

Frequency-independent impedances of soil-structure systems in horizontal and rocking modes

The frequency-independent foundation impedances, commonly used in soil-structure dynamic interaction problems, are developed for a circular footing resting on a homogeneous halfspace. As they ignore the structure attached to the foundation, the error introduced in the structural response may be 50 per cent or more in the neighbourhood of the fundamental frequency of the soil-structure system. The present study proposes a new method developed for most dynamic soil-structure interaction problems. The key idea is to retain for the frequency-independent impedances values computed for the fundamental frequency of the soil-structure system; thus these values include the dynamic characteristics of the whole soil-structure system and lead to a satisfactory approximation of the exact solution over a wide frequency range. The method is developed here for the horizontal and rocking modes of a structure with a circular base resting on a homogeneous halfspace. Numerical applications are given for a simple linear oscillator in order to make possible a thorough parametric study. The response of some idealized building-foundation systems to harmonic excitation or to a seismic input is next examined in order to illustrate the efficiency of the proposed model.     

On a Physically Realistic Fast Dynamo

As a step towards a physically realistic model of a fast dynamo, we study numerically a kinematic dynamo driven by convection in a rapidly rotating cylindrical annulus. Convection maintains the quasi-geostrophic balance whilst developing more complicated time-dependence as the Rayleigh number is increased. We incorporate the effects of Ekman suction and investigate dynamo action resulting from a chaotic flow obtained in this manner. We examine the growth rate as a function of magnetic Prandtl number Pm, which is proportional to the magnetic Reynolds number. Even for the largest value of Pm considered, a clearly identifiable asymptotic behaviour is not established. Nevertheless the available evidence strongly suggests a fast dynamo process.     

Obituary – Emerita Professor Valerie Anne Hall BSc PhD FSA FHEA (1946–2016)

Professor Emerita Valerie Hall died on 28th July, 2016, in Belfast City Hospital after a long battle with cancer. Valerie was a renowned and highly esteemed inter-disciplinarian, who began her career in botany and palynology but made a truly international mark in the field of cryptotephrochronology. Between us, we knew her as a student, a supervisor, an examiner, a mentor, a colleague, a role model, and, above all, a wonderful, kind-hearted and generous friend. She touched the hearts and earned the respect of Quaternary scientists around the world. Here, to honour her memory, we reflect briefly on the great woman and scientist that she was.     

Holocene grassland vegetation, climate and human impact in central eastern Inner Mongolia

~~Holocene grassland vegetation, climate and human impact in central eastern Inner Mongolia1.Ye,D.Z.,Chou,J.E,Liu,J.Y.et al.,Causes of sand-stormy weather in northern China and control measures,Acta Geographica Sinica(in Chinese),2000,55(5):513-520. 2.Qiu,X.E,Zeng,Y.,Miao,Q.L.,Temporal spatial distribution as well as tracks and source areas of sand dust storms in China,Acta Geographica Sinica(in Chinese),2001,56(3):316-322. 3.Cui,H.T,Kong,Z.C,Preliminary analysis on the cli…     

Radioactive high level waste insight modelling for geological disposal facilities

Within this paper we present a simplified analytical model to provide insight into the key performance measures of a generic disposal system for high level waste within a geological disposal facility. The model assumes a low solubility waste matrix within a corrosion resistant disposal container surrounded by a low permeability buffer. Radionuclides migrate from the disposal area through a porous geosphere to the biosphere and give a radiological dose to a receptor. The system of equations describing the migration is transformed into Laplace space and an approximation used to determine peak values for the radionuclide mass transfer rate entering the biosphere. Results from the model are compared with those from more detailed numerical models for key radionuclides in the UK high level waste inventory. Such an insight model can provide a valuable second line of argument to assist in confirming the results of more detailed models and build confidence in the safety case for a geological disposal facility.     

Shear partitioning in the Philippines: Constraints from Philippine Fault and global positioning system data

Abstract The Philippine Fault is a major left-lateral structure formed in an island arc setting. It accommodates a component of the oblique convergence between the Philippine Sea Plate and the Philippine archipelago. This observation is quantified through a series of global positioning satellite experiments between 1991 and 1996. The formation of the Fault marks the onset of a new geodynamic regime in the Philippine region. In the central Philippines, this event corresponds to the creation of a new tectonic boundary separating the Philippine Mobile Belt and the Philippine Sea Plate, following the latter's kinematic reorganization that occurred around 4 Ma ago. During this event, the Philippine Sea Plate changed its relative movement with respect to Eurasia from a northward to a north-westward motion, favoring the formation of a Philippine Fault–Philippine Trench system under a shear partitioning mechanism.     

Study on Paleoearthquakes of the Eastern Liupanshan Piedmont Fault Zone

Based on trench excavation and detailed field investigation, the authors have ascertained that 6 paleoearthquake events have occurred in the Haizixia area of the Eastern Liupanshan Piedmont fault zone since 46Ka B.P. The recurrence time of 6 events were about 15000a, 20250a, 14750a, 12150a, 8550a and 4000a B.P., respectively, the recurrence intervals were about 15000a, 5500a, 2600a, 3600a and 4550a, respectively. According to M-D empirical relations, the vertical dislocation of events 1 and 2 corresponded to the dislocation of the earthquake with M8.0; the vertical dislocation of events 3-6, was corresponded to the dislocation of earthquakes with M6.5~7.5.     

Improved water table dynamics in MODFLOW

The standard formulation of a block-centered finite-difference model, such as MODFLOW, uses the center of the cell as the location of a cell node. Simulations of a dynamic water table can be improved if the node of a cell containing the water table is located at the water table rather than at the center of the cell. The LPF package of MOD-FLOW-2000 was changed to position a cell's node at the water table in convertible cells with a water table. Improved accuracy in the upper regions of an unconfined aquifer is demonstrated for pumping from a partially penetrating well. The change introduces a nonlinearity into the solution of the flow equations that results in slightly slower convergence of the flow solution, 7% slower in the presented demonstration. Accuracy of simulations is improved where vertical flow is dominated by a moving water table, but not when a large water table gradient dominates over the water table movement.