利用现有重力场模型求定CHAMP卫星加速度计修正参数

CHAMP卫星加速度计数据的标定是通过确定其尺度因子和偏差参数来完成的.本文基于能量守恒方程,给出利用现有重力场模型标定CHAMP卫星加速度数据的基本原理和数学模型;提出相邻历元间差分算法,大大简化了观测方程,同时避免积分常量的计算.该算法既能同时解算尺度因子和偏差参数,也可任意求解其中之一.基于实测的CHAMP卫星加速度数据,利用EGM96模型和最新公布的EIGEN-2模型进行计算与比较,验证该方法的有效性.     

中国大陆地区中强地震前加卸载响应比异常时间尺度的统计研究

加卸载响应比理论是近年来提出的地震预测方法。本文进行了加卸载响应比异常时间尺度的统计研究，其中包括中国大陆地区5．0～8．1级的部分中强地震共30个震例，并得出加卸载响应比异常的时间尺度与未来地震震级之间的拟合函数。结果表明，地震前加卸载响应比异常的时间尺度与未来地震的震级具有正变关系，即震级越高，地震前加卸载响应比异常的时间尺度越长。根据加卸载响应比异常的时间尺度与震级之间的关系可以估计未来地震的发震时间，同时，可以确定加卸载响应比时空扫描过程中时问长度的大小。     

2-D crustal Poisson's ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.     

再论各向异性介质转换波地震学,第一部分:理论

我们业已研发了计算各向异性、非均质介质中P- SV转换波(C-波)的转换点和旅行时的新理论。据此 可以利用诸如相似性分析、迪克斯模型建模、克契 霍夫求和等常规方法来完成各向异性的处理和各向 异性处理,并使各向异性的处理成为可能。这里将 我们的新发展分作两部分来介绍。第一部分为理 论,第二部分为对速度分析和参数计算的应用。第 一部分理论包括转换点的计算和动校正的分析。     

Converted-wave seismology in anisotropic media revisited,Part I: Basic theory

We have developed new basic theories for calculating the conversion point and the travel time of the P-SV converted wave (C-wave) in anisotropic, inhomogeneous media. This enables the use of conventional procedures such as semblance analysis, Dix-type model building and Kirchhoff summation, to implement anisotropic processing, and makes anisotropic processing affordable. Here we present these new developments in two parts: basic theory and application to velocity analysis and parameter estimation. This part deals with the basic theory, including both conversion-point calculation and moveout analysis. Existing equations for calculating the PS-wave (C-wave) conversion point in layered media with vertical transverse isotropy (VTI) are strictly limited to offsets about half the reflector depth (an offset-depth ratio, xlz, of 0.5), and those for calculating the C-wave traveltimes are limited to offsets equal to the reflector depth (x/z=l.0). In contrast, the new equations for calculating the conversion-point extend into offsets about three-times the reflector depth (x/z=3.0), those for calculating the C-wave traveltimes extend into offsets twice the reflector depth (x/z=2.0). With the improved accuracy, the equations can help in C-wave data processing and parameter estimation in anisotropic, inhomogeneous media. This work is funded by the Edinburgh Anisotropy Project (EAP) of the British Geological Survey. First author: Xiangyang Li, Mr. Li is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project in the British Geological Survey. He also holds a honorary professorship in multicomponent seismology at the School of Geosciences, University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984) in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992) in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute. Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent seismology.     

2-D crustal Poisson’s ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L ₁ are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson’s ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson’s ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson’s ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson’s ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson’s ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson’s ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region. Foundation item: Key Project from China Earthquake Administration and the Project (95-11-02-01) from Ministry of Science and Technology (2001DIA10003). Contribution No. RCEG200401, Geophysical Exploration Center, China Earthquake Administration.     

Photosynthesis-transpiration coupling model at canopy scale in terrestrial ecosystem

At the hypothesis of big leaf, an ecosystem photosynthesis-transpiration coupling cycle model was established by the scaled SMPT-SB model from single leaf to canopy, and model parameterization methods were discussed. Through simulating the canopy light distribution, canopy internal conductance to CO2 can be scaled from single leaf to canopy by integrating to canopy using the relationship between single internal conductance and photosynthetic photon flux density. Using the data observed by eddy covariance method from the Changbai Mountains site of ChinaFLUX, the application of the model at the canopy scale was examined. Under no water stress, the simulated net ecosystem photosynthesis rate fitted with the observed data very well, the slope and R2 of the line regression equation of the observed and simulated values were 0.7977 and 0.8892, respectively (n = 752), and average absolute error was 3.78 μmol CO2 m-2s-1; the slope, R2 and average absolute error of transpiration rate were 0.7314, 0.4355 and 1.60mmol H2O m-2 s-1, respectively (n = 752). The relationship between canopy photosynthesis,transpiration and external environmental conditions was discussed by treating the canopy as a whole and neglecting the comprehensive feedback mechanism within canopy, and it was noted that the precipitation course affected the transpiration rate simulation badly. Compared to the models based on eco-physiological processes, the SMPT-SB model was simple and easy to be used. And it can be used as a basic carbon and water coupling model of soil-plant-atmosphere continuum.     

Changes in the spatial scale of Beijing UHI and urban development

The seasonal and interannual variations of Beijing urban heat island (UHI) are investigated in this paper using the temperature data from 1960 to 2000 at 20 meteorological stations in the Beijing region, and then the relationship between the intensity and spatial scale of UHI and Beijing urbanization indices is analyzed and discussed. Main conclusions are the followings. First, Beijing UHI shows obvious seasonal variations, and it is strongest in winter, next in spring and autumn, and least in summer. The seasonal variation of the UHI mainly occurs in the urban area. The UHI intensity at the center of Beijing is more than 0.8℃ in winter, and only 0.5℃ in summer. Second, the intensity of Beijing HUI exhibits a clear interannual warming trend with its mean growth rate (MGR) being 0.3088℃/10 a. The MGR of HUI is largest in winter, next in spring and autumn, and least in summer, and the urban temperature increase makes a major contribution to the growth of HUI intensity. Third, since the Reform and Opening, the urbanization indices have grown several ten times or even one hundred times, the intensity of HUI has increased dramatically, and its spatial scale also expanded distinctively along with the expansion of urban architectural complexes. Fourth, the interannual variation of urbanization indices is very similar with that of HUI intensity, and their linear correlation coefficients are significant at a more than 0.001 confidence level.     

Stability and accuracy analysis of the central difference method for real‐time substructure testing

The central difference method (CDM) that is explicit for pseudo‐dynamic testing is also believed to be explicit for real‐time substructure testing (RST). However, to obtain the correct velocity dependent restoring force of the physical substructure being tested, the target velocity is required to be calculated as well as the displacement. The standard CDM provides only explicit target displacement but not explicit target velocity. This paper investigates the required modification of the standard central difference method when applied to RST and analyzes the stability and accuracy of the modified CDM for RST. Copyright © 2005 John Wiley & Sons, Ltd.     

A practical method for estimating dynamic soil stiffness on surface of multi‐layered soil

It is important to estimate the influence of layered soil in soil–structure interaction analyses. Although a great number of investigations have been carried out on this subject, there are very few practical methods that do not require complex calculations. In this paper, a simple and practical method for estimating the horizontal dynamic stiffness of a rigid foundation on the surface of multi‐layered soil is proposed. In this method, waves propagating in the soil are traced using the conception of the cone model, and the impulse response function can be calculated directly and easily in the time domain with a good degree of accuracy. The characteristics of the impedance, that is the transformed value to the frequency domain of the obtained impulse response, are studied using two‐ to four‐layered soil models. The cause of the fluctuation of impedance is expressed clearly from its relation to reflected waves from the lower layer boundary in the model. Copyright © 2005 John Wiley & Sons, Ltd.