The transition from subduction to continental collision: crustal structure in the North Canterbury region, New Zealand

The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.     

Spectral analysis using orthonormal functions with a case study on the sea surface topography

b
Spherical harmonics are orthonormalized using the Gram-Schmidt process in a function space. The problem of linear dependence of spherical harmonics over the oceans is studied using the Gram matrices and consequently three sets of orthonormal (ON) functions have been constructed. For the process an efficient formula for computing inner products of spherical harmonics has been developed. Important spectral properties of the ON functions are addressed. The ON functions may be used for representing the sea surface topography (SST) in the analysis of satellite altimeter data. The geoid error can be transformed to a representation by the ON functions and hence the comparison of powers of the geoid error and the SST signal only over the oceans is possible, leading to a better way of determining the cut-off frequency of the SST in the simultaneous solution using satellite altimeter data. As a case study, the modified Levitus SST is expanded into the ON functions. The results show that 99.90 per cent of that signal's energy is contained within degree 24 of the orthonormal functions. Such expansions also render better spectral behaviour of oceanic signals as compared to that from spherical harmonic expansions. The study shows that these generalized Fourier functions are suitable for spectral analyses of oceanic signals and they can be applied to future altimetric mission where the geoid and the SST are to be recovered.     

Asymptotic eigenfrequencies of the Earth's normal modes

We derive asymptotic formulae for the toroidal and spheroidal eigenfrequencies of a SNREI earth model with two discontinuities, by considering the constructive interference of propagating SH and P-SV body waves. For a model with a smooth solid inner core, fluid outer core and mantle, there are four SH and 10 P-SV ray parameters regimes, each of which must be examined separately. The asymptotic eigenfrequency equations in each of these regimes depend only on the intercept times of the propagating wave types and the reflection and transmission coefficients of the waves at the free surface and the two discontinuities. If the classical geometrical plane-wave reflection and transmission coefficients are used, the final eigenfrequency equations are all real. In general, the asymptotic eigenfrequencies agree extremely well with the exact numerical eigenfrequencies; to illustrate this, we present comparisons for a crustless version of earth model 1066A.     

On the frequency contents of local events: source or path effect?

Local events which occurred in the Timanfaya volcanic field are analysed and their spectral contents are interpreted as being due to the effect of a finely layered medium with large impedance contrast, the presence of which is well known from well-logging data. By considering the layered medium as a filter, its spectral response (amplitude and phase) is computed from different randomly generated values of its parameters (total thickness, mean thickness of the individual layers and impedance contrast), showing that this stratigraphic filter acts as a low-pass filter. This effect can be correctly modelled by means of an exponential correlation function. However, in some cases resonant peaks appear in the spectrum superimposed on the low-pass filter effect, which can only be attributed to resonant scattering. An analysis of the phase and group delay associated with the amplitudes is made and strongly suggests that it is necessary to incorporate the theory of resonant scattering for the retrieval of medium properties. From the point of view of the amplitudes, the superposition of both effects, low-pass filter and resonant peaks, allows us to reproduce the gross features that are usually attributed to source effects. The need for a careful analysis of the local site response previous to any interpretation of the spectral features in terms of source characteristics is thus emphasized.     

Non-linear optimization for seismic traveltime tomography

This paper presents a non-linear algorithmic approach for seismic traveltime. It is based on large-scale optimization using non-linear least-squares and trust-region methods. These methods provide a natural way to stabilize algorithms based on Newton's iteration for non-linear minimization. They also correspond to an alternative (and often more efficient) view of the Levenberg-Marquardt method. Numerical experience on synthetic data and on real borehole-to-borehole problems are presented. In particular, results produced by the new algorithm are compared with those of Ivansson (1985) for the Kråkemåla experiment.     

乡村聚落地理研究的国外动态与国内趋势

本文根据国内外乡村聚落地理研究的内容变化和理论进展，初步划分了乡村聚落地理研究的主要发展阶段。在此基础上，结合国外乡村聚落地理研究的动态及我国乡村地区发展的现实，对我国乡村聚落地理研究进行了回顾与评价，并提出若干重点研究内容与趋势。     

反演地壳磁化强度等效源法的模型验证及实际资料试算

等效源法曾被广泛应用于卫星磁异常解释,而本文是将等效源法引用到航空△T磁场的解释中,用其研究区域地质构造。根据等效源法的理论,推导出球坐标系下球棱柱体的△T磁场公式,采用线性反演的方法进行模型验证,同时将该方法用于秦巴地区航空△T磁异常的试算,得到了一些有意义的结论。利用该方法可由所求得的磁场强度来得到居里深度分布。     

我国首例接触极化曲线野外实测结果

使用自行研制的野外测试设备,在云南省武定县迤纳厂铜铁矿床上进行了接触极化曲线法试验性的测试。初步证实,接触极化曲线法确有反映导电矿体成分和储量的能力,它的应用范围也不是非常局狭的。     

城乡一体化“现代地籍”建设关键问题研究

论文分析了现行地籍管理信息系统中存在的问题和发展趋势,对现代地籍的概念进行界定和阐述,进而研究探讨了城乡一体化“现代地籍”建设过程中的一些关键技术问题,并总结出可行的解决方案。