Computation of solitary waves during propagation and runup on a slope

A numerical time-simulation algorithm for analysing highly nonlinear solitary waves interacting with plane gentle and steep slopes is described by employing a mixed Eulerian–Lagrangian method. The full nonlinear free surface conditions are considered here in a Lagrangian frame of reference without any analytical approximations, and thus the method is valid for very steep waves including overturning. It is found that the runup height is crucially dependent on the wave steepness and the slope of the plane. Pressures and forces exerted on impermeable walls of different inclinations (slopes) by progressive shallow water solitary waves are studied. Strong nonlinear features in the form of pronounced double peaks are visible in the time history of pressure and force signals with increasing heights of the oncoming solitary waves. The effect of nonlinearity is less pronounced as the inclination of the wall decreases with respect to the bottom surface.     

Refraction Effect in Geothermal Heat Flow Due to a 3-D Prism in a Two-Layered Earth

A mathematical model of the disturbance of the stationary geothermal field due to a three-dimensional perturbing body embedded near the surface of a two layered earth. The theoretical analysis is based on the generalized theory of the double-layer potential, similar to the boundary integral method likely to direct current geoelectricity problems. Special attention is paid to cases in which a pan of the body's surface is tangential to the planar boundaries of the substratum or of the surface of the earth The numerical calculations are performed for a three-dimensional prism using methods similar to the boundary element method (BEM). Numerous graphs are shown for the disturbance of the temperature field and of the heat flow on the surface of the earth and along vertical profiles (boreholes).     

Multiscale Gravitational Field Recovery from GPS-Satellite-to-Satellite Tracking

The purpose of GPS-satellite-to-satellite tracking (GPS-SST) is to determine the gravitational potential at the earth's surface from measured ranges (geometrical distances) between a low-flying satellite and the high-flying satellites of the Global Positioning System (GPS). In this paper, GPS-satellite-to-satellite tracking is reformulated as the problem of determining the gravitational potential of the earth from given gradients at satellite altitude. The uniqueness and stability of the solution are investigated. The essential tool is to split the gradient field into a normal part (i.e. the first-order radial derivative) and a tangential part (i.e. the surface gradient). Uniqueness is proved for polar, circular orbits corresponding to both types of data (first radial derivative and/or surface gradient). In both cases gravity recovery based on satellite-to-satellite tracking turns out to be an exponentially ill-posed problem. Regularization in terms of spherical wavelets is proposed as an appropriate solution method, based on the knowledge of the singular system. Finally, the extension of this method is generalized to a nonspherical earth and a non-spherical orbital surface, based on combined terrestrial and satellite data.     

Vertical dynamic response of a disk on a saturated poroelastic half-space

This paper considers the vertical dynamic response of a disk on a saturated poroelastic half-space. Firstly the pressure-solid displacement form of the harmonic equations of motion for a poroelastic solid are developed from the form of the equations originally presented by Biot. These equations are solved by a new method. Then the mixed boundary value problem for the vertical harmonic vibration of a disk on a poroelastic half-space is studied. The two types of drainage conditions at the surface of the poroelastic half-space are considered: (a) the surface of the poroelastic half-space is assumed to be completely pervious both within and exterior to the plate; (b) The interface between the plate and the poroelastic half-space is assumed to be impervious and the exterior region is assumed to be pervious. By using the Hankel transform techniques, the paper develops the governing dual integral equations. These governing integral equations are further reduced to systems of standard Fredholm integral equations of the second kind by Abel transform.     

E3D,the Euro3D visualization tool I: Description of the program and its capabilities

We present the first version of E3D, the Euro3D visualization tool for data from integral field spectroscopy. We describe its major characteristics, based on the proposed requirements, the current state of the project, and some planned future upgrades. We show examples of its use and capabilities. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Euro3D data format: A common FITS data format for integral field spectrographs

As integral field spectrographs have become more common around the world and in Europe in particular, the need for a common data format has been recognized. Here, we present the Euro3D format that is adapted as a post‐instrumentalsignature‐removal format for all instruments within the Euro3D network. It follows the FITS standard, and includes several extensions, all of them being binary FITS tables. This article is intended to give a comprehensive overview, but does not attempt to serve as a full definition document. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

埋设地下光缆对武都地电影响的分析

着重探讨武都地电台测区埋设地下光缆后北西道地电阻率测值急剧上升的原因.通过实际观测、理论分析和实验,总结出了这种变化的机制,定量地解释了埋设地下光缆对地电观测的影响,从而为以后判断地电测值变化是干扰还是地震异常提供了依据.     

A BRIEF INTRODUCTION TO THE NEW METHOD FOR RIVER PROFILE ANALYSIS: Integral Approach

The topography and geomorphology of active orogens result from the interaction of tectonics and climate. In most orogens, a fluvial channel is most sensitive to the coupling between tectonics, lithology, and climate. Meanwhile, the related signals have been recorded by both the drainage geometry and channel longitudinal profile. Thus, how to extract tectonic information from fluvial channels has been a focused issue in geologic and geomorphologic studies. The well known stream-power river incision model bridges the gap between tectonic uplift, river incision and channel profile change, making it possible to retrieve rock uplift pattern from river profiles. In this model, the river incision rate depends on the rock erodibility, contributing drainage area and river gradient. The steady-state form of the river incision model predicts a power-law scaling between the drainage area and channel gradient. Via a linear regression to the log-transformed slope-area data, the slope and intercept are channel concavity and steepness indices, respectively. The concavity relates to lithology, climatic setting and incision process while the channel steepness can be used to map the spatial pattern of rock uplift. For its simple calculation process, the slope-area analysis has been widely used in the study of tectonic geomorphology during past decades. However, to calculate river slope, the coarse channel elevation data must be smoothed, re-sampled, and differentiated without any reasonable smooth window or rigid mathematical fundamentals. One may lose important information and derive stream-power parameters with high uncertainties. In this paper, we introduce the integral approach, a procedure that has been widely used in the latest four years and demonstrated to be a better method for river profile analysis than the traditional slope-area analysis. Via the integration to the steady-state form of the stream-power river incision equation, the river longitudinal profile can be converted into a straight line of which the independent variable is the integral quantity χ with the unit of distance and the dependent variable is the relative channel elevation. We can calculate the linear correlation coefficient between elevation and χ based on a series of concavity values and find the best linear fit to be the reasonable channel concavity index. The slope of the linear fit to the χ value and elevation is simply related to the ratio of the uplift rate to the erodibility. Without calculating channel slope, the integral approach makes up for the drawback of the slope-area analysis. Meanwhile, via the integral approach, a steady-state river profile can be expressed as a continuous function, which can provide theoretical principle for some geomorphic parameters (e.g., slope-length index, hypsometric integral). In addition, we can determine the drainage network migration direction using this method. Therefore, the integral approach can be used as a better method for tectonogeomorphic research.     

2017年8月9日精河6.6级地震前LURR异常时空演化特征

基于加卸载响应比理论,对2017年8月9日新疆精河6.6级地震前加卸载响应比时空演化特征和异常区的孕震积分时程曲线进行分析。结果表明,在精河6.6级地震前,研究区出现了明显的加卸载响应比高值异常区域,异常区呈现环状特征;而且,异常区域随时间逐渐向震源附近迁移,异常演化过程中异常区域内部和周围发生了6次5级以上地震,其中,2次6级地震发生在异常高值减弱过程中,异常峰值到发震的时间分别为11、22个月,大于以往研究的理论发震时间,这可能是受构造应力区域性特征影响的缘故。     

Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.