Three-dimensional elastic wave scattering and diffraction due to a rigid cylinder embedded in an elastic medium by a point source

The formal solutions of displacement field to the problem of elastic wave scattering and diffraction due to an infinitely long rigid cylinder embedded in an infinite elastic medium by an impulsive point source have been obtained in the integral form. The integrals for the reflected and the diffracted waves both in the shadow zone and in the illuminated zone are evaluated asymptotically for the early time motion by the Reisdue-Cagniard method and the Saddle-point-Cagniard method.Numerical results of the diffractedP, S andPS waves at a fixed circum-distance from the surface of the rigid cylinder show noticeably that (1) the energy partition for the diffractedS wave is small in comparison with that for the diffractedP wave, (2) the wave form of the diffractedS wave is broader and more diffused than that of the diffractedP wave, (3) the direction of the radial motions of the diffractedP andS waves varies as a function of the observational point, and (4) the energy partition for the diffractedP wave is much smaller than that for the direct or the reflectedP waves.This paper has been presented at the 46th Annual International Meeting of Society of Exploration Geophysicists in Houston, Texas, Oct. 28, 1976.     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

Determining magma flow in sills,dykes and laccoliths and their implications for sill emplacement mechanisms

Three-dimensional seismic data from the Faeroe-Shetland Basin provides detailed information on the relationships between sills, dykes, laccoliths and contemporaneous volcanic activity. The data shows that sills are predominantly concave upwards, being complete or partial versions of radially or bilaterally symmetrical forms that possess flat inner saucers connected to a flat outer rim by a steeply inclined sheet. Such morphologies are only partially modified by pre-existing faults. Sills can be sourced from dykes or the steep climbing portions of deeper sills. Both sills and dykes can provide magma to overlying volcanic fissures and sills can be shown to feed shallow laccoliths. Magma flow patterns, as revealed by opacity rendering, suggest that sills propagate upwards and outwards away from the magma feeder. As an individual sill can consist of several leaves emplaced at different stratigraphic levels, and as a sill or dyke can provide magma to volcanic fissures, other sills and laccoliths, the data suggests that neutral buoyancy concepts may not provide a complete explanation for the mechanism and level of sill emplacement. Instead, the data suggests that the presence of lithological contrasts, particularly ductile horizons such as overpressured shales may permit sill formation at any level below the neutrally buoyant level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Ken Thomson–deceased, April 2007     

Comments on “Steady- and transient-state inversion in hydrogeology by successive flux estimation” by P. Pasquier and D. Marcotte

Pasquier and Marcotte [Pasquier P, Marcotte D. Steady- and transient-state inversion in hydrogeology by successive flux estimation. Adv Wat Res 2006;29:1934–52] propose some modifications to the Comparison Model Method (CMM), in order to apply it to transient 3D ground water flow data for conductivity identification. We present some remarks on that paper to improve the comprehension of the basic features of the CMM and of the real value of the novelties introduced by Pasquier and Marcotte.     

激光扫描测量系统的应用研究

介绍了激光扫描测量技术的原理，并就激光扫描测量在工业中堆体的实时变化监测的技术和方法进行了探讨和研究，提出了一种基于激光扫描测量堆体变化的快速监测方法。实验表明：这种测量技术与传统的测量方法相比具有较强的实时性和较高的测量精度。     

用航摄像对制作全息地形立体图的新技术

全息地形立体图是将传统的光学机械立体投影术与全息术相结合，用航摄像对合成为一张具有一定焦深的近像面全息图。它既具有像面全息图衍射效率高，可以用扩展的白光光源再现的优点，也具有体全息图一定的波长和角度选择性。本文介绍了拍摄这种立体图的装置和方法，并从理论上对实验结果进行了一些分析，最后指出了它的应用前景。     

鲜水河断裂带古构造残余应力场对大地震的控制

以岩体正交异性弹性理论为基础,用X射线法,在鲜水河断裂带测区沿7条测线,测量了岩体三维残余应力,绘成区域残余应力等值线图、主应力线图及其应变能密度等值线图,探讨了其对带内大地震的控制作用     

Groebner-basis solution of the three-dimensional resection problem (P4P)

The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points P_i{X_i,Y_i,Z_i}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan.     

Steady- and transient-state inversion in hydrogeology by successive flux estimation

A calibration method to solve the groundwater inverse problem under steady- and transient-state conditions is presented. The method compares kriged and numerical head field gradients to modify hydraulic conductivity without the use of non-linear optimization techniques. The process is repeated iteratively until a close match with piezometric data is reached. The approach includes a damping factor to avoid divergence and oscillation of the solution in areas of low hydraulic gradient and a weighting factor to account for temporal head variation in transient simulations. The efficiency of the method in terms of computing time and calibration results is demonstrated with a synthetic field. It is shown that the proposed method provides parameter fields that reproduce both hydraulic conductivity and piezometric data in few forward model solutions. Stochastic numerical experiments are conducted to evaluate the sensitivity of the method to the damping function and to the head field estimation errors.     

Three-dimensional mapping of salt stores in the Murray-Darling Basin,Australia: 2. Calculating landscape salt loads from airborne electromagnetic and laboratory data

Salt is widespread in the Australian landscape – in soil, regolith and groundwater – leading to concerns that land management practices may be putting much agricultural land and important water resources at risk of salinisation. Defining the location and nature of salt stores is an important first step in understanding the processes leading to salinity of soils, streams and groundwater resources, and predicting areas that may be at risk. Airborne geophysics can define subsurface salt stores and mobilisation pathways. Airborne electromagnetics (AEM) can map the three-dimensional conductivity structure of the landscape but does not, intrinsically, quantify the amount of salt. Salinity, moisture content, porosity and mineralogy all contribute to the electromagnetic signal, and each can vary significantly throughout the landscape.In the Lower Balonne catchment, Queensland, the relationship between AEM and the amount of salt in the landscape was quantified using laboratory analyses of pore fluids from core samples. A general statistical relationship was established between AEM conductivity and salt load (defined as the product of pore fluid salinity, porosity and moisture content)—with a significant positive correlation although data were generally widely dispersed. Comparison of calculated salt load with borehole electromagnetic logs gives insight into the factors contributing to dispersion in the AEM data.The relationship transforms bulk conductivity to salt load in 5 m layers, allowing the generation of a three-dimensional map of the salt load. This is a powerful tool for identifying areas that may require monitoring and management interventions to reduce salinity risk. An example is given of salt loads beneath an established irrigation area in the Lower Balonne catchment, Queensland.