A theoretical study on terrestrial gravimetric data refinement by Earth gravity models

The idea of this paper is to present estimators for combining terrestrial gravity data with Earth gravity models and produce a high‐quality source of the Earth's gravity field data through all wavelengths. To do so, integral and point‐wise estimators are mathematically developed, based on the spectral combination theory, in such a way that they combine terrestrial data with one and/or two Earth gravity models. The integral estimators are developed so that they become biased or unbiased to a priori information. For testing the quality of the estimators, their global mean square errors are generated using an Earth gravity model08 model and one of the recent products of the gravity field and steady‐state ocean circulation explorer mission. Numerical results show that the integral estimators have smaller global root mean square errors than the point‐wise ones but they are not efficient practically. The integral estimator of the biased type is the most suited due to its smallest global root mean square error comparing to the rest of the estimators. Due largely to the omission errors of Earth gravity models the point‐wise estimators are not sensitive to the Earth gravity model commission error; therefore, the use of high‐degree Earth gravity models is very influential for reduction of their root mean square errors. Also it is shown that the use of the ocean circulation explorer Earth gravity model does not significantly reduce the root mean square errors of the presented estimators in the presence of Earth gravity model08. All estimators are applied in the region of Fennoscandia and a cap size of 2° for numerical integration and a maximum degree of 2500 for generation of band‐limited kernels are found suitable for the integral estimators.     

地震与核爆识别的小波包分量比方法

频谱分析法在核爆与地震识别中具有广泛的应用．但是频谱分析方法是稳态方法,即使采用Gabor变换,也因时-频窗口形状不变而分辨串较低．为提高时-频分辨率,本文将小波变换理论用于乌鲁木齐台记录的地震与核爆事件的分析,并提出了识别核爆和天然地震的小波包分量比判据．通过对加拿大黄刀地震台记录的印度地下核爆的分析,进一步验证了小波包分量比判据对核爆和地震的识别具有较高的识别效率．结果表明：对于地震信号,其小波包分量比U03/U1一般都大于1.0,而对于核爆信号,比值U03／U13一般都小于1.0.     

安庆4.8级地震序列的地震波谱变化研究

采用安徽区域数字地震台网记录的安庆4.8级地震序列波形资料,开展地震波谱研究,求取了零频位移谱值、应力降和地震矩等震源参数。分析了纵、横波零频谱值比和应力降随时间的变化趋势,并讨论了震源参数之间的关系。结果表明,安庆4.8级主震后,纵、横波零频位移谱值比和应力降呈现一定的变化规律,可为主震—余震序列的最大余震的预测提供一定的参考依据,同时发现震源参数与震级有一定相关。     

Evaluation of codified elastic design spectrum models for regions of low-to-moderate seismicity

Design spectrum (DS) model is typically specified in a seismic code of practice for structural design. In a region of low-to-moderate seismicity where seismic code does not exist, a DS model in a well established code of practice is usually adopted, while the suitability of such model has seldom been evaluated. In this article, the elastic DS models for reference (rock) site stipulated in six major codes of practice (AS1170.4–2007, EN1998-1:2004, GB50011–2010, IBC–2012, NBCC–2010 and NZS1170.5:2004) have been compared and scrutinized. Three cities of low-to-moderate seismicity, namely, Melbourne (Australia), Hong Kong (China) and Karlsruhe (Germany), have been selected for illustrative purposes. Particular emphasis has been put on the parameterization scheme for DS model. It is found that huge discrepancies (over 100%) exist among the models, especially at the long period range, due to differences in spectral shapes and the recommended corner periods, which would lead to undesirable effects on the use of the displacement-based seismic design approach. It is urged that the values of corner periods should be determined specifically and cautiously based on the regional seismicity pattern and local geological conditions.     

Air flow calculation of reverse circulation drilling technique with air DTH hammer

Air DTH hammer has been successfully applied in minor-caliber solid mineral exploration,water-well drilling and other drilling areas.In order to expand the applications of the technology,the authors further studied the principle and analyzed the mechanism of reverse circulation drilling technique with air DTH hammer to get the perfect assembles of equipments by optimizing working parameters.No parameter seemed more important than the air volume because it could maintain the working performance stability.The minimum air volume is related to the parameters such as depth and pressure,which was calculated under the actual conditions.It was solved for the air injection flow tables of the air DTH Hammer working at the different pressures.According to the data tables,operators could adjust the air volume to meet the demand on this technique,which had a realistic guiding significance.So it could build up a set of systematic and complete hi-technique.     

Coseismic displacement estimate of the 2013 <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>7.0 Lushan,China earthquake based on the simulation of near-fault displacement field

Usually, GPS observation provides direct evidence to estimate coseismic displacement. However, GPS stations are scattered, sparse and cannot provide a detailed distribution of coseismic displacement. Strong ground motion records share the same disadvantages as GPS in estimating coseismic displacement. Estimations from InSAR data can provide displacement distributions; however, the resolution of such methods is limited by the analysis techniques. The paper focuses on estimating the coseismic displacement of the M _S7.0 Lushan earthquake on April 20, 2013 using a simulation of the wave field based on the elastic wave equation instead of a quasi-static equation. First, the media and source models were constructed by comparing the simulated velocity and the record velocity of the ground motion. Then simulated static displacements were compared with GPS records. Their agreement validates our results. Careful analysis of the distribution of simulated coseismic displacements near the fault reveals more details of the ground motion. For example, an uplift appears on the hanging wall of the fault, rotation is associated with the horizontal displacement, the fault strike and earthquake epicenter provide the main control on motion near the faults, and the motion on the hanging wall is stronger than that on the footwall. These results reveal additional characteristics of the ground motion of the Lushan earthquake.     

The assimilation of spectral sensing and the WOFOST model for the dynamic simulation of cadmium accumulation in rice tissues

The accurate detection of heavy metal-induced stress on crop growth is important for food security and agricultural, ecological and environmental protection. Spectral sensing offers an efficient and undamaged observation tool to monitor soil and vegetation contamination. This study proposed a methodology for dynamically estimating the total cadmium (Cd) accumulation in rice tissues by assimilating spectral information into WOFOST (World Food Study) model. Based on the differences among ground hyperspectral data of rice in three experiments fields under different Cd concentration levels, the spectral indices MCARI1, NREP and RH were selected to reflect the rice stress condition and dry matter production of rice. With assimilating these sensitive spectral indices into the WOFOST + PROSPECT + SAIL model to optimize the Cd pollution stress factor f_wi, the dynamic dry matter production processes of rice were adjusted. Based on the relation between dry matter production and Cd accumulation, we dynamically simulating the Cd accumulation in rice tissues. The results showed that the method performed well in dynamically estimating the total amount of Cd accumulation in rice tissues with R² over 85%. This study suggests that the proposed method of integrating the spectral information and the crop growth model could successfully dynamically simulate the Cd accumulation in rice tissues.     

Introducing sensor spectral response into the classification process

Many sensors have their bands overlapped and therefore do not set a normal space. If a spectral distance is measured, as in first-order statistical classifiers, the direct consequence is that the result will not be the most accurate. Image classification processes are independent of the spectral response function of the sensor, so this overlap is usually ignored during image processing. This paper presents a methodology that introduces the spectral response function of sensors into the classification process to increase its accuracy. This process takes place in two steps: first, incident energy values of the sensors are reconstructed; second, the energy of the bands is set in an orthonormal space using a matrix singular value decomposition. Sensors with and without overlapping spectral bands were simulated to evaluate the reconstruction of energy values. The whole process was implemented on three types of images with medium, high and very high spatial resolution obtained with the sensors ASTER, IKONOS and DMC camera, respectively. These images were classified by ISODATA and minimum distance algorithms. The ISODATA classifier showed well-defined features in the processed images, while the results were less clear in the original images. At the same time, the minimum distance classifier showed that overall accuracy of the processed images increased as the maximum tolerance distance decreased compared to the original images.     

A spectral gradient difference based approach for land cover change detection

Change detection with remotely sensed imagery plays an important role in land cover mapping, process analysis and dynamic information services. Euclidean distance, correlation and other mathematic metrics between spectral curves have been used to calculate change magnitude in most change detection methods. However, many pseudo changes would also be detected because of inter-class spectral variance, which remains a significant challenge for operational remote sensing applications. In general, different land cover types have their own spectral curves characterized by typical spectral values and shapes. These spectral values are widely used for designing change detection algorithms. However, the shape of spectral curves has not yet been fully considered. This paper proposes to use spectral gradient difference (SGD) to quantitatively describe the spectral shapes and the differences in shape between two spectra. Change magnitude calculated in the new spectral gradient space is used to detect the change/no-change areas. Then, a chain model is employed to represent the SGD pattern both qualitatively and quantitatively. Finally, the land cover change types are determined by pattern matching with the knowledgebase of reference SGD patterns. The effectiveness of this SGD-based change detection approach was verified by a simulation experiment and a case study of Landsat data. The results indicated that the SGD-based approach was superior to the traditional methods.