Biological and Geochemical Record of Anthropogenic Impacts in Recent Sediments from Lake Pupuke, Auckland City, New Zealand

Frozen sediment cores from Lake Pupuke in Auckland City, New Zealand, contain a high resolution decadal to annual scale record of changing lake paleoenvironments and geochemistry that reflects changing landuse and hydrology in the catchment over the past c. 190 years. A reliable chronology is available from AMS ¹⁴C and ²¹⁰Pb dating of the sediments, with the timing of the older part of the record supported by the first appearance of pollen of introduced Pinus and Plantago lanceolata associated with European settlement of Auckland in the early 1840s. Diatom stratigraphy, sediment elemental and carbon isotope geochemistry reflect changes in sediment sources and lake conditions commensurate with European development of the Lake Pupuke catchment, in particular enhanced algal productivity controlled by the influx of nutrients after c. 1920 AD. Attempts to prevent nuisance algal blooms in 1933, 1934 and 1939 using CuSO₄ addition produced Cu spikes in the sediment that allowed confirmation of the accuracy of the ²¹⁰Pb chronology. Hence, the elemental and isotopic composition of the Lake Pupuke sediments reflect the timing of significant anthropogenic activities, rather than climatic variations, that have occurred within the watershed over the past c. 190 years. The comparison of records of land use change in the catchment with the multi-proxy record of changes in the sediments usually allowed unambiguous identification of the signatures of change and their causes.     

Computation of Large Anisotropic Seismic Heterogeneities (CLASH)

A general tomographic technique is designed in order (i) to operate in anisotropic media; (ii) to account for the uneven seismic sampling and (iii) to handle massive data sets in a reasonable computing time. One modus operandi to compute a 3-D body wave velocity model relies on surface wave phase velocity measurements. An intermediate step, shared by other approaches, consists in translating, for each period of a given mode branch, the phase velocities integrated along ray paths into local velocity perturbations. To this end, we develop a method, which accounts for the azimuthal anisotropy in its comprehensive form. The weakly non-linear forward problem allows to use a conjugate gradient optimization. The Earth's surface is regularly discretized and the partial derivatives are assigned to the individual grid points. Possible lack of lateral resolution, due to the inescapable uneven ray path coverage, is taken into account through the a priori covariances on parameters with laterally variable correlation lengths. This method allows to efficiently separate the 2ψ and the 4ψ anisotropic effects from the isotropic perturbations. Fundamental mode and overtone phase velocity maps, derived with real Rayleigh wave data sets, are presented and compared with previous maps. The isotropic models concur well with the results of Trampert & Woodhouse. Large 4ψ heterogeneities are located in the tectonically active regions and over the continental lithospheres such as North America, Antarctica or Australia. At various periods, a significant 4ψ signature is correlated with the Hawaii hotspot track. Finally, concurring with the conclusions of Trampert & Woodhouse, our phase velocity maps show that Rayleigh wave data sets do need both 2ψ and 4ψ anisotropic terms.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

Seismic Hazards Assessment for Radioactive Waste Disposal Sites in Regions of Low Seismic Activity

A comparative tectonic quiescence and lack of earthquakes make the stable centres of continents attractive for siting long-term radioactive waste storage facilities. The low rates of deformation in such regions, however, make it difficult to characterize their long-term seismotectonic behaviour, leading to uncertain estimates for the very low probability hazard estimates required by society. In an attempt to overcome the deficiency of both contemporary seismicity and paleoseismic data in central Canada, we have used earthquake histories from regions with similar seismotectonic characteristics from around the world. Substituting space for time, we estimate a long-term rate per 10⁶ km² of 0.004 magnitude ≥6 earthquakes per annum, of which 33–100% might rupture to the surface.     

我国西部1∶5万地形图空白区地物解译研究——以西藏那曲其香错地区为例

我国西部部分地区1∶5万地形图数据尚未完成,给区域和国家经济发展带来极大不便。空间技术、计算机技术和信息技术的高速发展为遥感技术应用提供了技术支持,高分辨率影像获取能力的提高为大比例尺地物特征的提取提供了数据可能。本文采用SPOT5遥感影像,以1∶5万地形图图式标准为依据,进行青藏高原其香错地区地物判读,结果表明SPOT5能基本满足1∶5万地形图制图要求。为进一步提高判读结果,可采用多时相、多分辨率影像进行1∶5万地形图地物判读。     

地学自由曲面径向基函数网络重建的详细算法

径向基函数(RadialBasisFunction,简称RBF)神经网络是一种理想的地学离散数据网格化工具,能够适应各种不同分布形式和边界条件的数据,收敛速度较快,可以逼近任何复杂曲面。这里详细介绍了RBF神经网络的算法。适当径向基函数的形式和偏差系数是使用RBF神经网络作地学曲面重建的关键。大量的实际数据验证结果表明,当选用Gauss型径向基函数时,一般可获得比较理想的网格化效果,同时具备"曲面平滑"和"拟合度高"的特点。