On imperfection of elasticity in the Earth's interior

Summary. Attempts are made to estimate the distribution of damping with depth, based on the observed damping of P waves, of free vibrations, and of surface waves and partly on a discussion by D. L. Anderson and R. S. Hart. Results are consistent in order of magnitude, but could easily be wrong by a factor of 2. Comments are added on the need for further work. Stress is laid on the damping of the free nutation, which indicates that α in the modified Lomnitz law is about 0.2 on an average.     

Expansion of a high-frequency time-harmonic wavefield given on an initial surface into Gaussian beams

Summary. Realistic anelastic attenuation laws are usually formulated as convolution operators, but this representation is intractable for time-domain synthetic seismogram methods such as the finite difference method. An approach based on Padé approximants provides a convenient, accurate reformulation of general anelastic laws in differential form. The resulting differential operators form a uniformly convergent sequence of increasing order in the time derivative, and all are shown to be causal, stable and dissi- pative. In the special case of frequency-independent Q , all required coefficients for the operators are obtained in closed form in terms of Legendre polynomials.
Low-order approximants are surprisingly accurate. Finite-difference impulse responses for a plane wave in a constant- Q medium, calculated with the fifth-order convergent, are virtually indistinguishable from the exact solution. The formulation is easily generalized to non-scalar waves. Moreover, this method provides a framework for incorporating amplitude-dependent attenuation into numerical simulations.     

Normal modes of the viscoelastic earth

Summary A uniformly valid linear viscoelastic rheology is described which takes the form of a 'generalized' Burgers' body and which appears capable of reconciling the behaviour of the Earth's mantle across the complete spectrum of geodynamic time-scales. This spectrum is bracketed by the short time-scales of body wave and free oscillation seismology on which anelastic effects are dominant, and the long time-scale of mantle convection on which the Earth behaves viscously. The parameters of the model which control the viscous response are fixed by post-glacial rebound data whereas those which govern the anelasticity are to be determined by fitting the model to observations of seismic Q. The paper is concerned primarily with a discussion of the normal mode spectrum of the Earth as a generalized Burgers' body. Focusing upon the homogeneous model, it includes an initial analysis of the accuracy of first-order perturbation theory as a method of calculating the respective Q s of the elastic gravitational free oscillations. Also considered are the quasi-static modes of relaxation which only exact eigenanalysis can reveal. The importance of these modes is assessed within the context of a discussion of the effect of viscoelasticity upon the efficiency of Chandler wobble excitation.     

Waveform inversion of mantle Love waves:the Born seismogram approach

Summary. Normal mode theory, extended to the slightly laterally heterogeneous earth by the first-order Born approximation, is applied to the waveform inversion of mantle Love wave (200–500 s) for the Earth's lateral heterogeneity at l = 2 and a spherically symmétric anelasticity ( Q _μ) structure. The data are from the Global Digital Seismograph Network (GDSN). The l =2 pattern is very similar to the results of other studies that used either different méthods, such as phase velocity measurements and multiplet location measurements, or a different data set, such as mantle Rayleigh waves from different instruments. The results are carefully analysed for variance reduction and are most naturally explained by heterogeneity in the upper 420 km. Because of the poor resolution of the data set for the deep interior, however, a fairly large heterogeneity in the transition zones, of the order of up to 3.5 per cent in shear wave velocity, is allowed. It is noteworthy that Love waves of this period range cannot constrain the structure below 420 km and thus any model presented by similar studies below this depth are likely to be constrained by Rayleigh waves (spheroidal modes) only.
The calculated modal Q values for the obtained Q _μ model fall within the error bars of the observations. The result demonstrates the discrepancy of Rayleigh wave Q and Love wave Q and indicates that care must be taken when both Rayleigh and Love wave data, including amplitude information, are inverted simultaneously.
Anomalous amplitude inversions of G2 and G3, for example, are observed for some source-receiver pairs. This is due to multipathing effects. One example near the epicentral region, which is modelled by the obtained l = 2 heterogeneity, is shown.     

近岸二类水体生物光学模型参数优化

为了优化生物光学模型中的悬浮物后向散射模型系数（b^_bx）和双向反射因子（f/Q）,基于珠江口、韩江河口、徐闻珊瑚礁保护区实测数据（遥感反射率、叶绿素a质量浓度、悬浮物质量浓度、CDOM吸收系数）,采用模拟退火算法同时优化了b_bx以及f/Q（N＝43）,并开发了生物光学正演模型（N＝22）。文章优化的f/Q为0.104 9,悬浮物在参考波长531 nm处后向散射系数b_bx为0.268×[TSS]^0.295。通过输入实测的水质数据模拟出531 nm处的遥感反射率,并外推到其他波段（412、443、490、551、667 nm）。模拟得到的531 nm处遥感反射率与实测值的均方根误差RMSE＝0.000 8,N＝22,除了模拟的667 nm处遥感反射率误差较大外（RMSE＝0.003 6,N＝22）,其他波段都具有较好的精度（RMSE＜0.002 3）。结果表明,该研究区与其他研究区（如黄、东海）的bbx和f/Q不一致,这2个参数需要根据研究区的数据进行优化。验证结果说明迭代优化的b_bx以及f/Q能够适用于广东近岸河口生物光学参数模拟。     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

The coupling and attenuation of nearly resonant multiplets in the Earth's free oscillation spectrum

Summary. The effect of attenuation on the coupling of nearly resonant multiplets in the Earth's free oscillation spectrum is investigated and numerical results are presented for several of the most strongly coupled low frequency multiplet pairs. The coupling influences considered are those of the Coriolis forces due to rotation and of the Earth's hydrostatic ellipticity of figure. It is found that the effects of attenuation (in particular the difference in Q ⁻¹ for the two multiplets) can significantly change the splitting diagrams and the degree to which coupling takes place. The Q values for the coupled singlets are, in general, all different and lie between the two Q values of the uncoupled multiplet pair. In addition it is shown that the diagonal sum rule may be readily extended to sets of coupled multiplets.     

城市节庆活动的整合与可持续发展——以昆明市为例

通过分析＇99世博会前后节庆活动的发展对旅游业带来的影响,研究了昆明城市节庆活动的整合和可持续发展问题.分析表明：通过举办不同类型、不同规模、不同时间长短的系列节庆活动,对节庆活动进行整合,可以实现城市节庆活动的可持续发展,从而促进城市作为旅游目的地的品牌化进程.在此基础上归纳了一个城市节庆开发的CSD理论框架.     

土地利用对街头抢劫影响的空间分异模式：以H市为例

选取H市中心城区2009年和2010年发生的街头抢劫案件为研究对象,构建地理加权泊松回归GWPR模型,探讨不同性质土地利用对案件影响的空间差异性。研究表明,各类用地与街头抢劫的关系存在显著的空间异质性。其中,位于中心区外围的传统居住区和部分新建社区等居住用地对街头抢劫存在显著的负向影响;商业服务业设施用地对街头抢劫的影响总体上显著性不强,但与其中的商业设施点类型密切相关,独栋高层的商业综合体会使犯罪数量下降,而临街商铺形式的商业设施点则是犯罪高发区域;工业用地仅对老城区的街头抢劫产生正向影响;公园绿地总体上与街头抢劫案件呈负相关性,但城市开敞绿地多的地方案件数量也较多;交通设施用地与街头抢劫的关系则受到交通设施、道路等级、路网形态的影响。最后通过聚类分析发现,土地利用对案件影响的主导因素的空间分异模式与不同的城市功能分区有一定的对应关系,可以对每个功能分区根据其街头抢劫发生的主导因素制定相应的土地利用规划导则。     

Equal arc ratio projection and a new spherical triangle quadtree model

A projection family, named Equal Arc Ratio Projection (EARP), is proposed in this paper, which has two kinds of geometric cases. Coordinates in the projections are determined by the ratios between four special arcs. The projection for spherical octahedron is derived, which is available for Dutton's Quaternary Triangular Mesh (QTM), and both projections of Goodchild's and Otoo's for QTM are found to be two special implementations of EARP. The icosahedral projection and thus the corresponding spherical triangle quadtree partition model are also developed. The partition model is compared with several other global grid models.     

Modelling the compliance of crustal rock—I. Response of shear-wave splitting to differential stress

We show that seismic shear waves may be used to monitor the in situ stress state of deep inaccessible rocks in the crust. The most widespread manifestation of the stress-related behaviour of seismic waves is the shear-wave splitting (shear-wave birefringence) observed in almost all rocks, where the polarizations of the leading split shear waves are usually subparallel to the direction of the local maximum horizontal stress. It has been recognized that such shear-wave splitting is typically the result of propagation through distributions of stress-aligned fluid-filled microcracks and pores, known as extensive-dilatancy anisotropy or EDA. This paper provides a quantitative basis for the EDA hypothesis. We model the evolution of anisotropic distributions of microcracks in triaxial differential stress, where the driving mechanism is fluid migration along pressure gradients between neighbouring microcracks and pores at different orientations to the stress field. This leads to a non-linear anisotropic poroelasticity (APE) model for the stress-sensitive behaviour of fluid-saturated microcracked rocks. A companion paper shows that APE modelling matches a range of observed phenomena and is a good approximation to the equation of state of a stressed fluid-saturated rock mass.     

An experiment on the elevation accuracy improvement of photogrammetrically derived DEM

This paper focuses on a topic barely considered in the literature: how to improve the accuracy of a given Digital Elevation Model (DEM) irrespective of its lineage by identifying its most suspicious values (also denoted here as outliers). Methods tailored to a specific procedure and source (contour maps, remote sensing image, etc.) exist but they are not valid for other cases. This is a problem for both the producer and end user. The results of a comparison of two methods using six DEMs intended to be representative of different landscapes are reported here. Both methods have been applied to each DEM, producing a number of height candidates to be analysed. Assuming that all candidates are wrong, their elevations have been blindly replaced by interpolated heights, simulating the behaviour of the inexperienced user. The improved (or degraded) DEMs are then compared against the ground truth, and updated accuracy figures are calculated. The RMSE can diminish by 2 to 8% of the original value by changing less than 1% of the elevations in the dataset.     

SELECTED GRAPHS OF WORLDWIDE,CLOUDLESS SOLAR RADIATION INCIDENT ON ALL SLOPES

The amounts of global solar radiation (Q + q) incident on all slope inclinations and latitudes have been computed for both summer and winter. These calculations include the hour-by-hour amounts, as well as the daily sums of clear-day shortwave radiation on southern, western, and northern slopes. Some of the results are presented in graphical, three-dimensional form in order to gain an improved worldwide perspective of the varied responses of impinging Q + q as a complex function of different latitudes, slope angles and directions, and solar azimuths and altitudes under changing seasonal regimes.     

Determination of Lg-wave attenuation using single-station seismograms: a case study for Zimbabwe

In this paper we show how the quality factor Q may be calculated using a single seismograph station and a number of events recorded on analogue seismograms. We followed Nuttli's (1973) method and extended it to one seismograph station. Using the single station Bulawayo (BUL), we determined a mean Q value of 650 for Zimbabwe. Furthermore, we considered different propagation paths over Zimbabwe as the seismic waves travelled to BUL and found a low Q value of 350 for the Deka fault zone. the Q value of 650 obtained in this study agrees well within error with that of 603 reported by Chow et al. (1980) using the multistation, multi-event method, and lies within the range of Q values (454–759) reported by Xie & Mitchell (1990) using a 'back-project'method to image large-scale lateral variations of Lg coda Q . the results obtained are important as the Q value constitutes part of the input data in seismic hazard calculations. the method may be used to determine Q in regions where there may be only one station with reliable analogue seismogram data.     

南极长城站地区气象要素的中期振荡特征

本文用多维谱分析的方法研究了1985年1月-1987年12月南极长城站地面气象要素的谱特征。该地区的气象要素与北半球一样,也普遍存在3-5天、准一周、准二周和3-4周振荡周期。在长城站所在的乔治王岛地区,从地面到平流层中下部,气压系统都存在着显著的准二周振荡,这种振荡是由上往下传播的。无论冬、夏季,准二周和准一周振荡都是由西往东传播的,而扰动的南北传播方向则与季节及振荡周期长度有关。南极地区的中期振荡特征与500hPa极地涡旋的活动和强度变化有密切关系。     

An individual cognitive evacuation behaviour model for agent-based simulation: a case study of a large outdoor event

Route selection is a complex issue in simulating the pedestrian’s walking behaviour during the evacuation. In many current studies, the agent-based model has been widely used to simulate the individualized behaviour, and growing concern is given to the cognitive aspects in route selection. However, the influence of incomplete spatial knowledge and the preferences for different routing strategies has not been fully considered. To represent more reliable route choice by different pedestrians, this study presents an individual cognitive evacuation behaviour model which is mainly focusing on the individual differences in routing strategy and spatial knowledge. The model consists of two major components, namely the multiple routing strategies and personalized spatial cognitive road network. Based on a review of the previous studies, the potential factors that may affect the evacuation behaviour were discussed. The quantitative relationships between these potential factors and the pedestrian’s routing preference, as well as the possibilities of recognizing different roads, were explored using the regression method. As a case study, the proposed model was implemented using the data collected from a field survey of a real outdoor event. When applied in an agent-based simulation, the implemented model will support more reliable simulation of route choices corresponding to the pedestrians’ behaviour in the real event.     

A new class of stratified viscoelastic models by analytical techniques

Multilayer, spherically stratified, self-gravitating relaxation models with a large number of layers (more than 100) can be dealt with analytically. Relaxation processes are studied for both Heaviside surface loads and tidal forcings. Simulations of the relaxation process of a realistic earth model with an incompressible Maxwell rheology show that models containing about 30 to 40 layers have reached continuum limits on all timescales and for all harmonic degrees up to at least 150 whenever an elastic lithosphere is present, irrespective of the viscosity profile in the mantle. In particular, fine-graded stratification of the shallow layers proves to be important for high harmonic degrees in these models. The models produce correct long-time (fluid) limits. It is shown that differences in the transient behaviour of the various models are due to the applied volume-averaging procedure of the rheological parameters. Our earlier proposed hypothesis that purported shortcomings in the fundamental physics of (discrete) normal-mode theory are artificial consequences of numerical inaccuracies, theoretical misinterpretations and the use of incomplete sets of normal modes is reinforced by the results presented. We show explicitly that the models produce both continuous behaviour resulting from continuous rheological stratifications and discrete behaviour resulting from sharp density contrasts, as at the outer surface and the core-mantle boundary. The differences between volume-averaged models and fixed-boundary contrast models are outlined. Reducing many-layer models with a volume-averaging procedure before employing a normal mode analysis is both economical and highly accurate on all timescales and for all spherical harmonic degrees. The procedure minimizes the chances of missing contributing modes, while using models with more layers will not result in any substantial increase of accuracy.     

The British who are not British and the immigration policies that are not: the case of Hong Kong

"While 'closed-door' immigration policies are adopted by most countries, 'exceptionalist' legislation is often made to permit entry of special immigrant groups. An example is the British Nationality (Hong Kong) Act 1990, which was passed in the run-up to the change in sovereignty of Hong Kong in 1997. Britain's increasingly restrictive immigration policies prior to 1990 [have] resulted in the majority of Hong Kong citizens having British nationality (as British Dependent Territories citizens) but without the right of abode in the U.K. The 1990 Act conferred full British citizenship status on 50,000 heads of households in Hong Kong." The authors conclude that "in a world of marked global inequalities, immigration pressure will become even more extreme and is likely to produce an increasing number of cases of exceptionalist immigration legislation in countries with both ?open' and 'closed'-door policies.'     

黄土高原多沙粗沙区侵蚀模型探讨

国内外曾研究出不同类型的坡面侵蚀计算模型,但还没有一个成熟的流域侵蚀计算模型,至今流域侵蚀量仍以输沙量替代。本文充分利用黄河中游不同区域具有大量淤地坝的条件,再通过侵蚀影响因素机理的分析、在侵蚀形态类型区的划分等基础上探讨侵蚀变权模型的建设。首先选择了155个“闷葫芦”淤地坝,采集每个坝的年均淤积量,以及相应流域的侵蚀影响因素：植被盖度、降雨量、沟谷密度、切割深度、地表组成物质、>15°的坡耕地面积比等资料,采用变权形式,建立侵蚀强度宏观估算模型。最后用年降水量与年输沙量的关系对模型进行改进。