Generalized ‘dependencies’ and parameter variance

The autodependence, (a special case of the — now quite obsolete — dependences, which had been introduced for very specialized astrometric purposes) is proportional to the parameter variance which is the expectation of the variance of the systematic error of a function evaluated with estimated parameters.     

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 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.     

Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

Drainage spacing regularity on a fault-block: A case study from the eastern Ruahine Range

Abstract:  Recent research has indicated river basin outlets draining linear sections of large, uplifting mountain belts often show a regularity of spacing, transverse to the main structural trend. A morphometric analysis of part of the Ruahine Range, on the North Island was undertaken to test whether drainage regularity may exist in smaller, younger mountain ranges. The ratio, R , of the half-width of the mountain belt, W , and the outlet spacing, S , was used to characterize drainage networks on the eastern side of the range. The spacing ratio for the range of 1.31 is lower than R results from studies of larger mountain belts ( R  = 1.91–2.23). We suggest the cause of this lower ratio is related to eastward migration of the Ruahine drainage divide.     

Some Observations on the Use of Data from Historical Documents in Debris-Flow Studies

The knowledge of past events is important for the assessment of debris-flow hazard. Amongst the sources of information, documents from historical archives are particularly important in sites where the debris flows cause damage to urban areas and transportation routes. The paper analyses the availability of historical documents on debris flows in Northeastern Italy and discusses factors that can influence the building of time series from archive data both at regional and single basin scales. An increased number of debris flows was observed in the studied region for the last decades. This could be due both to an increased frequency of the events and to a larger availability of information: the analysis carried out indicates that the latter factor is probably the most influencing. The importance of factors, which affect the collection of data, including the conservation of documents and the presence and fragility of the elements at risk, is stressed in view of a wise use of historical data on debris flows.     

深层地震勘探激发井深探讨

深层地震勘探野外资料采集的关键是在保证原始资料一定信噪比的基础上,同时兼顾分辨率,激发井深是决定资料品质的重要因素之一。通过对虚反射滤波作用分析,提出了充分利用虚反射能量设计最佳激发井深,增强下传地震波能量的采集方法。理论分析和试验资料表明,在高速层下4～7m的最佳激发岩性深度是深层勘探的理想井深。