Simultaneous iterative time-domain sparse deconvolution to teleseismic receiver functions

In this paper, we present a new approach to estimate high-resolution teleseismic receiver functions using a simultaneous iterative time-domain sparse deconvolution. This technique improves the deconvolution by using reweighting strategies based on a Cauchy criterion. The resulting sparse receiver functions enhance the primary converted phases and its multiples. To test its functionality and reliability, we applied this approach to synthetic experiments and to seismic data recorded at station ABU, in Japan. Our results show Ps conversions at approximately 4.0 s after the primary P onset, which are consistent with other seismological studies in this area. We demonstrate that the sparse deconvolution is a simple, efficient technique in computing receiver functions with significantly greater resolution than conventional approaches.     

SEAWAT 2000: modelling unstable flow and sensitivity to discretization levels and numerical schemes

A systematic analysis shows how results from the finite difference code SEAWAT are sensitive to choice of grid dimension, time step, and numerical scheme for unstable flow problems. Guidelines to assist in selecting appropriate combinations of these factors are suggested. While the SEAWAT code has been tested for a wide range of problems, the sensitivity of results to spatial and temporal discretization levels and numerical schemes has not been studied in detail for unstable flow problems. Here, the Elder-Voss-Souza benchmark problem has been used to systematically explore the sensitivity of SEAWAT output to spatio-temporal resolution and numerical solver choice. A grid size of 0.38 and 0.60% of the total domain length and depth respectively is found to be fine enough to deliver results with acceptable accuracy for most of the numerical schemes when Courant number (Cr) is 0.1. All numerical solvers produced similar results for extremely fine meshes; however, some schemes converged faster than others. For instance, the 3rd-order total variation-diminishing method (TVD3) scheme converged at a much coarser mesh than the standard finite difference methods (SFDM) upstream weighting (UW) scheme. The sensitivity of the results to Cr number depends on the numerical scheme as expected.     

Modeling doline populations with logistic growth functions

Initiation and evolution of a large karst depression population consisting of parent and daughter dolines were spatially and temporally modeled using logistic growth functions. Logistic growth models are well suited for analyzing doline population initiation and evolution because they reflect the density‐dependent growth mechanisms present in the evolution of karst depressions. Seven assumptions based upon previous studies were refined into mathematical statements and tested using more than 2000 dolines from a subpopulation of both parent and daughter karst depressions on the Western Highland Rim and Pennyroyal Plain of Tennessee and Kentucky. Logistic growth models quantify the initiation and evolution of doline populations and interface well with recent models describing the evolution of three‐dimensional conduit systems. Logistic growth models should apply to modeling other doline populations in karst terranes with hydraulically efficient, three‐dimensional conduit systems. Copyright © 2006 John Wiley & Sons, Ltd.     

变形监测网的质量测度

对变形监测网络的质量测度－精度、灵敏度和可区分度进行了讨论，重点论述了它们之间的关系及应用的条件，并以算例加以论证。     

水稻水分生产函数时空变异规律研究

基于对水稻Jensen模型中敏感指数在全生育期变化规律的认识,以生长曲线函数建立了水稻敏感指数累积函数,分析了其特征,从而解决了不同长时段敏感指数转换计算问题。较全面地揭示了水稻水分生产函数及其敏感指数累积函数中主要参数随气象条件及土壤因子变化的规律。通过参照作物需水量及其频率以及不同地区土壤有效含水量为媒介,建立了对水分敏感指标在不同水文年份 (时间)和不同地区 (空间)进行预报的数学模型,据此提出了水稻水分生产函数在时、空两方面插补、延长、移用与扩展的理论与方法。借助于参照作物需水量等值线图及土壤分布图,探讨了绘制水分生产函数及其主要参数等值线图的原理和方法。     

Short-term precursor of the geomagnetic horizontal field transfer functions

IntroductionThe history of studying the changes in deep geoelectric conductivity of an earthquake preparation using the parameters of geomagnetic short-period variations has been about 30 years. At first, Japanese calculated the real transfer functions as well as the a and b coefficients of Parkinson vector by taking the measurement of graph and using empirical formula. There are some examples, such as the 1923 Kanto M=7.9 earthquake (Yanagihara, 1972), the 1966 Tashkent M=5.5 earthquake, t…     

地下水有机污染的水文地球化学标志物探讨——以河南油田为例

经过野外现场调查和取样分析及室内研究得知南阳油田地下水已遭受不同程度的有机物污染，且污染范围可能进一步向油田南部扩散。根据近似地下水流线方向上地下水中总油质量浓度和Fe,Mn等无机组分的变化势态，结合含水介质化学分析结果，发现硫酸盐、Fe和Mn可作为地下水有机污染和地球化学标志物。在含水介质中Fe和Mn质量浓度较高的地方，地下水有机污染物降解速度快，含水介质中Fe和Mn的氧化物和氢氧化物的还原作用导致了含水层介质中Fe和Mn的缺乏和地下水中溶解Fe和Mn的积聚；在含水介质中Fe和Mn质量浓度低的地方，地下水中的有机物质量浓度并没有降低，相应地地下水中溶解Fe和Mn的质量浓度也很低。同时，由于有机污染物的存在使地下水中硫酸盐被还原，导致城下水中硫酸盐质量浓度偏低，且地下水中Fe对有机物污染的敏感性比Mn强。     

Evolutionary properties of stochastic models of earthquake accelerograms: Their dependence on magnitude and distance

An approach to generate artificial earthquakeaccelerograms on hard soil sites is presented. Eachtime-history of accelerations is considered as arealization of a non-stationary gaussian stochasticprocess, with statistical parameters depending onmagnitude and source-to-site distance. In order tolink the values of these parameters for each groundmotion record with the corresponding magnitude andsource-to-site distance, semi-empirical functionalrelations called generalized attenuationfunctions are determined. The set of realground-motion time histories used to obtain thesefunctions correspond to shocks generated at differentsources and recorded at different sites in thevicinity of the southern coast of Mexico. The resultsshow significant dispersion in the parameters of themodel adopted, which reflect that associated with thereal earthquakes included in the sample employed.The problem of conditional simulation of artificialacceleration time histories for prescribed intensitiesis briefly presented, but its detailed study is leftfor a companion paper. The criteria and modelsproposed are applied to generate two families ofartificial acceleration records for recurrenceintervals of 100 and 200 years at a specific sitelocated in the region under study. The results shownin this article correspond to acceleration timehistories recorded on firm ground for earthquakesgenerated at the subduction zone that runs along thesouthern coast of Mexico, and cannot be generalized tocases of earthquakes generated at other sources orrecorded at other types of local conditions. Thismeans that the methods and functional forms presentedhere are applicable to these other cases, but thevalues of the parameters that characterize thosefunctions may differ from those presented here.     

Neuro‐control of seismically excited steel structure through sensitivity evaluation scheme

The neuro‐controller training algorithm based on cost function is applied to a multi‐degree‐of‐freedom system; and a sensitivity evaluation algorithm replacing the emulator neural network is proposed. In conventional methods, the emulator neural network is used to evaluate the sensitivity of structural response to the control signal. To use the emulator, it should be trained to predict the dynamic response of the structure. Much of the time is usually spent on training of the emulator. In the proposed algorithm, however, it takes only one sampling time to obtain the sensitivity. Therefore, training time for the emulator is eliminated. As a result, only one neural network is used for the neuro‐control system. In the numerical example, the three‐storey building structure with linear and non‐linear stiffness is controlled by the trained neural network. The actuator dynamics and control time delay are considered in the simulation. Numerical examples show that the proposed control algorithm is valid in structural control. Copyright © 2001 John Wiley & Sons, Ltd.