Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

区域地下水信息管理系统研究

本文通过《河南省地下水动态信息管理系统》的建立,讨论了区域地下水信息管理系统的设计目标,数据结构和软件实现,着重介绍了系统的软件设计特点,最后指出了这种管理系统在水文地质部门中具有广泛的应用前景。     

A realistic model for point-sources imaged on array detectors: The model and initial results

We have constructed a computer model for simulation of point-sources imaged on two-dimensional detectors. An attempt has been made to ensure that the model produces data that mimic real data taken with 2-D detectors. To be realistic, such simulations must include randomly generated noise of the appropriate type from all sources (e.g. source, background, and detector). The model is generic and accepts input values for parameters such as pixel size, read noise, source magnitude, and sky brightness. Point-source profiles are then generated with noise and detector characteristics added via our model. The synthetic data are output as simple integrations (onedimensional), as radial slices (two-dimensional), and as intensity-contour plots (three-dimensional). Each noise source can be turned on or off so that they can be studied separately as well as in combination to yield a realistic view of an image. This paper presents the basic properties of the model and some examples of how it can be used to simulate the effects of changing image position, image scale, signal strength, noise characteristics, and data reduction procedures.Use of the model has allowed us to confirm and quantify three points: 1) The use of traditionalsize apertures for photometry of faint point-sources adds substantial noise to the measurement which can significantly degrade the quality of the observation; 2) The number of pixels used to estimate the background is important and must be considered when estimating errors; and 3) The CCD equation normally used by the astronomical community consistently overestimates the signal-to-noise obtainable by a measurement while a revised equation, discussed here, provides a better estimator.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

用于湿地气候效应模拟的三江平原下垫面数据获取

三江平原是我国沼泽湿地的重要分布区。在20世纪的后50 a中,三江平原在经过几次大规模的开荒后,其下垫面发生了明显的阶段性变化,湿地大面积消失,农田面积不断增加。为了认识20世纪下半叶不同时期三江平原湿地的气候效应,拟采用第三代区域气候模式(RegCM3)来开展研究。而目前被广泛用于RegCM3的GL-CC数据,所反映的研究区下垫面几乎不存在沼泽湿地,这与事实严重不符,所以必须获取三江平原有沼泽湿地存在的真实的下垫面数据,才能得到可靠的湿地气候效应的模拟结果。分别选择了4期(1954年、1985年、1996年和2000年)有阶段代表性的航片、TM影像数据为主要数据源,利用ERDAS8.7和Arcview3.2等软件,首先得到4期研究区的土地利用数据,然后通过土地利用数据类型向模式数据类型的转换,从4期研究区土地利用数据中提取出了三江平原下垫面的模式数据。经过误差分析,认为获取的4期三江平原下垫面模式数据准确的反映出研究区几个时期的湿地下垫面实况,可以作为三江平原湿地气候效应模拟的基础数据。     

Bayesian Spatial Modelling of Gamma Ray Count Data

Gamma ray logging is a method routinely employed by geophysicists and environmental engineers in site geology evaluations. Modelling of gamma ray data from individual boreholes assists in the local identification of major lithological changes; modelling these data from a network of boreholes assists with lithological mapping and spatial stratigraphic correlation. In this paper we employ Bayesian spatial partition models to analyse gamma ray data spatially. In particular, a spatial partition is defined via a Voronoi tessellation and the mean intensity is assumed constant in each cell of the partition. The number of vertices generating the tessellation as well as the locations of vertices are assumed unknown, and uncertainty about these quantities is described via a hierarchical prior distribution. We describe the advantages of the spatial partition modelling approach in the context of smoothing gamma ray count data and describe an implementation that may be extended to the fitting of a more general model than a constant mean within each cell of the partition. As an illustration of the methodology we consider a data set collected from a network of eight boreholes, which is part of a geophysical study to assist in mapping the lithology of a site. Gamma ray logs are linked with geological information from cores and the spatial analysis of log data assists with predicting the lithology at unsampled locations.     

Gold mining in Ghana's forest reserves: a report on the current debate

This paper examines the debate surrounding a recent decision made by the Ghanaian government to permit gold exploration – and potentially, mining – in 'protected' forest reserves. In 2001, four mining companies were awarded mineral exploration concessions in forested regions of the country, and have since put forward applications to mine for gold. Notwithstanding the sharp divide in opinion on the issue, the continued uncertainty surrounding the implications of the proposed activities makes further research on the ground imperative in the short term. Work aiming to elicit indigenous perspectives on the projects, as well as research that facilitates dialogue between and/or among stakeholder parties, should be prioritized.     

采矿引起的边坡倾倒滑移变形机理与变形安全性分析研究—以抚顺西露天矿边坡为例

在矿山开发、水电建设过程中倾倒变形问题日渐突出。倾倒变形边坡一般具有反倾边坡结构,变形的过程和机理比较复杂、涉及的变形岩体范围较大、危害严重。论文以抚顺西露天矿边坡为例,通过岩体位移监测资料的跟踪分析、底面摩擦模拟试验,以岩体力学、松散介质力学原理为理论基础,探讨了采矿影响下边坡倾倒滑移体的变形机理,分析认为采矿引起的倾倒滑移变形问题是在特定的地质构造和采矿工程条件下产生的,并提出了岩体变形安全性分析方法。通过工程实例,验证了方法的实用性,并针对实际问题提出了变形整治对策。     

吾祠煤矿区构造特征及矿井生产布局探讨

通过对漳平市吾祠煤矿区的煤层地质构造特征分析，探讨了复杂构造条件下的矿井生产布局，主要是煤矿区段水平划分、煤炭采区的巷道布置、采煤方法及有关地质硐探等问题，对矿山开发有一定的指导意义。