Long-period seismic source inversions using global tomographic models

We investigate the effect of laterally varying earth structure on centroid moment tensor inversions using fundamental mode mantle waves. Theoretical seismograms are calculated using a full formulation of surface wave ray theory. Calculations are made using a variety of global tomographic earth models. Results are compared with those obtained using the so-called great-circle approximation, which assumes that phase corrections are given in terms of mean phase slowness along the great circle, and which neglects amplitude effects of heterogeneity. Synthetic tests suggest that even source parameters which fit the data very well may have large errors due to incomplete knowledge of lateral heterogeneity. The method is applied to 31 shallow, large earthquakes. For a given earthquake, the focal mechanisms calculated using different earth models and different forward modelling techniques can significantly vary. We provide a range of selected solutions based on the fit to the data, rather than one single solution. Difficulties in constraining the dip-slip components of the seismic moment tensor often produce overestimates of seismic moment, leading to near vertical dip-slip mechanisms. This happens more commonly for earth models not fitting the data well, confirming that more accurate modelling of lateral heterogeneity can help to constrain the dip-slip components of the seismic moment tensor.     

Seismic modeling of Tertiary sandstone clinothems, Spitsbergen

Two-dimensional seismic modelling has been undertaken on an overall progradational succession of sloping mudstone and sandstone units from the Palaeogene of Spitsbergen. The modelling shows that the main geometric features of the section would be resolved at 1500 m depth (with frequencies below 60 Hz, which is common in seismic data at these depths). However, interference between the base and top of lithological units gives lateral amplitude variations and discrepancies between the seismic image and the geometrical model. This is particularly prominent in low-frequency models. Terminations of reflectors, resembling toplap and onlap, may be interpreted, but are artefacts of the general convergence of lithological units present in the geometrical model. The geological section causes a seismic pattern resembling sigmoid progradational seismic facies. Two-dimensional seismic modelling is an efficient tool in bridging the gap between outcrop observations and subsurface data. Hence, modelled outcrop sections are important as reference points' for improved seismic stratigraphic interpretation.     

A method for calculating synthetic seismograms in laterally varying media

Summary An effective algorithm for computing synthetic seismograms in laterally inhomogeneous media has been developed. The method, based on zero-order asymptotic ray theory, is primarily intended for use in refraction and reflection studies and provides an economical means of seismic modelling.
A given smoothed velocity-depth-distance model is divided into small squares with constant seismic parameters and first-order interfaces are represented by an arbitrary number of dipping linear segments. The computation of ray propagation and amplitudes through such a model does not involve complicated analytic expressions and therefore minimizes computer time.
Amplitudes are determined by geometrical spreading of spherical wave-fronts and energy partitioning at interfaces. Synthetic seismograms calculated for laterally homogeneous models are in good agreement with those obtained by the Reflectivity Method.     

Teleseismic waveform modelling with a one-way wave equation

It is now widely accepted that elastic properties of the continental lithosphere and the underlying sublithospheric mantle are both anisotropic and laterally heterogeneous at a range of scales. To fully exploit modern three-component broad-band array data sets requires the use of comprehensive modelling tools. In this work, we investigate the use of a wide-angle, one-way wave equation to model variations in teleseismic 3-D waveforms due to 2-D elastic heterogeneity and anisotropy. The one-way operators are derived based on a high-frequency approximation of the square-root operator and include the effects of wave propagation as well as multiple scattering. Computational cost is reduced through a number of physically motivated approximations. We present synthetic results from simple 1-D (layer over a half-space) and 2-D (subduction zone) models that are compared with reference solutions. The algorithm is then used to model data from an array of broad-band seismograph stations deployed in northwestern Canada as part of the IRIS-PASSCAL/LITHOPROBE CANOE experiment. In this region radial-component receiver functions show a clear continental Moho and the presence of crustal material dipping into the mantle at the suture of two Palaeo-Proterozoic terranes. The geometry of the suture is better defined on the transverse component where subduction is associated with a ∼10 km thick layer exhibiting strong elastic anisotropy. The modelling reproduces the main features of the receiver functions, including the effects of anisotropy, heterogeneity and finite-frequency scattering.     

Wide-angle seismic velocities in heterogeneous crust

Seismic velocities measured by wide-angle surveys are commonly used to constrain material composition in the deep crust. Therefore, it is important to understand how these velocities are affected by the presence of multiscale heterogeneities. The effects may be characterised by the scale of the heterogeneity relative to the dominant seismic wavelength (λ); what is clear is that heterogeneities of all scales and strengths bias wide-angle velocities to some degree. Waveform modelling was used to investigate the apparent wide-angle P -wave velocities of different heterogeneous lower crusts. A constant composition (50 per cent felsic and 50 per cent ultramafic) was formed into a variety of 1- and 2-D heterogeneous arrangements and the resulting wide-angle seismic velocity was estimated. Elastic, 1-D models produced the largest velocity shift relative to the true average velocity of the medium (which is the velocity of an isotropic mixture of the two components). Thick (width > λ) horizontal layers, as a result of Fermat's Principle, provided the largest increase in velocity; thin (width ≪λ) vertical layers produced the largest decrease in velocity. Acoustic 2-D algorithms were shown to be inadequate for modelling the kinematics of waves in bodies with multiscale heterogeneities. Elastic, 2-D modelling found velocity shifts (both positive and negative) that were of a smaller magnitude than those produced by 1-D models. The key to the magnitude of the velocity shift appears to be the connectivity of the fast (and/or slow) components. Thus, the models with the highest apparent levels of connectivity between the fast phases, the 1-D layers, produced the highest-magnitude velocity shifts. To understand the relationship between measured seismic velocities and petrology in the deep crust it is clear that high-resolution structural information (which describes such connectivity) must be included in any modelling.     

Synthetic tsunami mareograms for realistic oceanic models

We study how the tsunami mode is generated by a scaled double-couple seismic source, and how it propagates in realistic oceanic models. The method developed and used is the direct extension to tsunami waves propagating in multilayered oceanic media of the well-known Haskell method. The most intensive tsunamis may be expected from sources located within the crust in the deep-water parts of the ocean. The extension to laterally heterogeneous structures shows that, if the thickness of the ocean liquid layer diminishes, the maximum amplitude of the tsunami wave train increases.     

Scattering characteristics in heterogeneously fractured reservoirs from waveform estimation

Offset-dependent characteristics of seismic scattering are useful for characterizing fractured reservoirs. We use two models that have different background medium properties and different azimuthal AVO responses to study elastic wave propagation and scattering in gas-saturated, heterogeneously fractured reservoirs. Heterogeneous fracture density distributions are built through stochastic modelling. Synthetic seismograms are generated by 3-D finite difference modelling, and waveforms along crack-normal and strike directions are considered in this paper. The multiple signal classification (MUSIC) frequency estimator is used in waveform estimation to provide frequency-domain attributes related to seismic wave scattering by fracture heterogeneity. Our results indicate that the strength of the scattering field is a function of the background medium. The strength also increases with increasing fracture scatterer density and with decreasing correlation length of spatial variations of fracture density. The scattering field is weak at the top of the fractured reservoir. The first-order results are dominated by velocity anisotropy of the mean fracture density field. However, the base of the fractured reservoir corresponds to a strong scattering field on which fracture heterogeneity has a larger effect and is characterized by the loss of coherence.     

Seismic wave properties in time-dependent porosity homogeneous media

It is quantified the properties of seismic waves in fully saturated homogeneous porous media within the framework of Sahay's modified and reformulated poroelastic theory. The computational results comprise amplitude attenuation, velocity dispersion and seismic waveforms. They show that the behaviour of all four waves modelled as a function of offset, frequency, porosity, fluid viscosity and source bandwidth depicts realistic dissipation within the sonic–ultrasonic band. Therefore, it appears that there is no need to include material heterogeneity to model attenuation. By inference it is concluded that the fluid viscosity effects may be enhanced by dynamic porosity.     

利用场模型拓扑表达GIS中的地理目标

该文阐述为了实现在GIS中描述自然界里带有模糊不确定性的地理目标 ,扩充GIS的模糊查询功能 ,首先基于集合论思想探讨了确定性地理目标的代数解析表达与其拓扑表达 ,指出了确定性点、线、面之间拓扑表达的构成机理 ;然后 ,基于场模型建立了模糊地理目标的空间表达 ,分析了模糊地理目标的位置不确定性 ,进而建立了模糊地理目标的拓扑表达模型 ,即点集拓扑内部、边界和外部。分析表明 ,公认的Egenhofer模型是该文模型在地理目标不带有误差或不确定性情况下的特例。最后 ,与Clementini(1996 )提出的模型做了比较分析 ,表明了该文模型的合理性。     

Numerical modelling of wavefields in three-dimensional inhomogeneous media

Summary. Numerical modelling is one of the most efficient methods for an investigation of the relationship between structural features and peculiarities of observed wavefields. It is practically the only method for 2-D and 3-D inhomogeneous media.
An algorithm based on ray theory has been developed for calculations of travel times and amplitudes of seismic waves in 3-D inhomogeneous media with curved interfaces. It was applied for numerical modelling of kinematic and dynamic characteristics of seismic waves propagating in laterally inhomogeneous media.
Travel-time and amplitude patterns were studied in the 2-D and 3-D models of a geosyncline, in which velocity distribution was given by an analytical function of the coordinates. For a more complicated model representing a subducting high-velocity lithospheric plate in a transition zone between oceanic and continental upper mantle, the velocity distribution was given by discrete values on a 2-D non-rectangular grid. It was shown that when a source was placed above the lithospheric plate, a shadow zone appeared along a strike of the structure, i.e. in the direction which is perpendicular to a strong lateral velocity gradient. Travel-time residuals were calculated along the seismological profile for a 3-D velocity distribution in the upper mantle beneath Central Asia, obtained as a result of inversion of travel times by the Backus-Gilbert method. They were found to be in a good agreement with the observed data.     

基于单元复形的数字图像多级别层次表达

基于复形理论定义了数字图像空间的拓扑元素及其性质,在此基础上提出一套完备的保持拓扑等价性的层次表达数字图像的数学模型体系框架,并验证了层次表达结构中的Jordan曲线定理.同时,基于单元复形扩展模型,对SPOT影像实施了渐进式离散分割,有效地利用影像蕴含的空间信息,获得了比最大似然分类法更优的分割结果.     

The first-order statistical moment of the seismic moment tensor

Summary. If a complex earthquake is assumed to be a set of individual, randomly oriented elementary pure double couple sources, the solution for the seismic moment of the complex event projected on the mean trend of the fault will perforce be comprised of sources of both double couple and compensated linear vector dipole (CLVD) types. We investigate the statistical properties of these two components of seismic sources in terms of the invariants of the seismic moment tensor of a realistic set of synthetic earthquakes. It is very likely that the size of the CLVD component is two to three orders of magnitude smaller than that of the double couple component.     

Chaotic ray behaviour in regional seismology

Rays propagating through strongly laterally varying media exhibit chaotic behaviour. This means that initially close rays diverge exponentially, rather than according to a power law. This chaotic behaviour is especially pronounced if the medium contains laterally varying interfaces. By studying simple 2-D and 3-D versions of models with laterally varying interfaces, the importance of chaotic ray behaviour is determined. A model of the Moho below Germany produces sharp variations with epicentral distance of the number of arrivals. In addition, the number of caustics grows dramatically: up to 1200 caustics are present between a distance of 0 and 800 km. Using the theory of Hamiltonian systems, a more in-depth study of the chaotic character of the ray equations is obtained. It is found that for realistic heterogeneous models most of the relevant rays will exhibit chaotic behaviour. The degree of chaos is quantified in terms of predictability horizons. Beyond the predictability horizons ray tracing cannot be carried out accurately. For the models under consideration, the length from the source to the predictability horizon has an order of magnitude of 1000 km. The chaotic behaviour of the rays makes it necessary to use extensions of asymptotic ray theory, such as Maslov theory, to compute seismic waveforms. It is shown that pseudo-caustics, an important obstacle in computing Maslov synthetics, are a generic feature of the 2-D laterally varying models that are studied. Eventually, the use of asymptotic methods is restricted because of the inaccuracy in the computation of the ray paths.     

Integration of linear and areal hierarchies for continuous multi-scale representation of road networks

Spatial data can be represented at different scales, and this leads to the issue of multi-scale spatial representation. Multi-scale spatial representation has been widely applied to online mapping products (e.g., Google Maps and Yahoo Maps). However, in most current products, multi-scale representation can only be achieved through a series of maps at fixed scales, resulting in a discontinuity (i.e., with jumps) in the transformation between scales and a mismatch between the available scales and users' desired scales. Therefore, it is very desirable to achieve smoothly continuous multi-scale spatial representations. This article describes an integrated approach to build a hierarchical structure of a road network for continuous multi-scale representation purposes, especially continuous selective omission of roads in a network. In this hierarchical structure, the linear and areal hierarchies are constructed, respectively, using two existing approaches for the linear and areal patterns in a road network. Continuous multi-scale representation of a road network can be achieved by searching in these hierarchies. This approach is validated by applying it to two study areas, and the results are evaluated by both quantitative analysis with two measures (i.e., similarity and average connectivity) and visual inspection. Experimental results show that this integrated approach performs better than existing approaches, especially in terms of preservation of connectivity and patterns of a road network. With this approach, efficient and continuous multi-scale selective omission of road networks becomes feasible.     

Urban expansion in Auckland,New Zealand: a GIS simulation via an intelligent self-adapting multiscale agent-based model

ABSTRACT

Abstract: When modelling urban expansion dynamics, cellular automata models focus mostly on the physical environments and cell neighbours, but ignore the ‘human’ aspect of the allocation of urban expansion cells. This limitation is overcome here using an intelligent self-adapting multiscale agent-based model. To simulate the urban expansion of Auckland, New Zealand, a total of 15 urban expansion drivers/constraints were considered over two periods (2000–2005, 2005–2010). The modelling takes into consideration both a macro-scale agent (government) and micro-scale agents (residents of three income levels), and their multi-level interactions. In order to achieve reliable simulation results, ABM was coupled with an artificial neural network to reveal the learning process and heterogeneity of the multi-sub-residential agents. The ANN-ABM accurately simulated the urban expansion of Auckland at both the global and local scales, with kappa simulation value at 0.48 and 0.55, respectively. The validated simulation result shows that the intelligent and self-adapting ANN-ABM approach is more accurate than an ABM with a general type of agent model (kappa simulation = 0.42) at the global scale, and more accurate than an ANN-based CA model (kappa simulation = 0.47) at the local scale. Simulation inaccuracy stems mostly from the outdated master land use plan.     

Paraxial ray approximations in the computation of seismic wavefields in inhomogeneous media

Summary. Several important applications of the paraxial ray approximation (PRA) to numerical modelling of high-frequency seismic body wavefields are discussed. The PRA can be used to evaluate the displacement vector not only directly on the ray, as in the standard ray method. but also approximately in the vicinity of this ray. The PRA also offers simple ways of approximate evaluation of paraxial rays, situated in the vicinity of the central ray, and of two-point ray tracing. A very important application of the PRA consists in a simple, fast and effective Computation of body-wave synthetic seismograms in general, 3-D, laterally inhomogeneous, layered structures. Examples of synthetic seismograms for 3-D structures, computed using the PRA, are presented.     

Evaluating representation and scale error in the maximal covering location problem using GIS and intelligent areal interpolation

A common problem in location-allocation modeling is the error associated with the representation and scale of demand. Numerous researchers have investigated aggregation errors associated with using different scaled data, and more recently, error associated with the geographic representation of model objects has also been studied. For covering problems, the validity of using polygon centroid representations of demand has been questioned by researchers, but the alternative has been to assume that demand is uniformly distributed within areal units. The spatial heterogeneity of demand within areal units thus has been modeled using one of two extremes – demand is completely concentrated at one location or demand is uniformly distributed. This article proposes using intelligent areal interpolation and geographic information systems to model the spatial heterogeneity of demand within spatial units when solving the maximal covering location problem. The results are compared against representations that assume demand is either concentrated at centroids or uniformly distributed. Using measures of scale and representation error, preliminary results from the test study indicate that for smaller scale data, representation has a substantial impact on model error whereas at larger scales, model error is not that different for the alternative representations of the distribution of demand within areal units.     

景观格局与生态过程相互关系研究进展

景观格局的形成反映了不同的景观生态过程,与此同时景观格局又在一定程度上影响着景观的演变过程。 空间格局、生态过程与尺度之间的相互关系是景观生态学研究的核心内容。本文对景观格局与生态过程的研究方 法、景观格局与生态过程相互关系等方面的国内外研究进展进行了综述。最后分析了目前景观格局与生态过程相 互关系研究中存在的问题与不足、未来的研究热点及发展趋势。