一种改进的基于非高斯性最大化的预测反褶积算法

The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.     

Source depletion and extent of melting in the Tongan sub-arc mantle

The fluid immobile High Field Strength Elements (HFSE) Nb and Ta can be used to distinguish between the effects of variable extents of melting and prior source depletion of the Tongan sub-arc mantle. Melting of spinel lherzolite beneath the Lau Basin back-arc spreading centres has the ability to fractionate Nb from Ta due to the greater compatibility of the latter in clinopyroxene. The identified spatial variation in plate velocities and separation of melt extraction zones, combined with extremely depleted lavas make Tonga an ideal setting in which to test models for arc melt generation and the role of back-arc magmatism.We present new data acquired by laser ablation-ICPMS of fused sample glasses produced without the use of a melt fluxing agent. The results show an arc trend towards strongly sub-chondritic Nb/Ta (< 17) with values as low as 7.2. Melting models show that large degree melts of depleted MORB mantle fail to reproduce the observed Nb/Ta. Alternatively, incorporation of residual back-arc mantle that has undergone less than 1% melting into the sub-arc melting regime reproduces arc values. However, the extent of partial melting required to produce the composition of the Lau Basin back-arc basalts averages 7%. This apparent discrepancy can be explained if only the lowermost 4 km of the residua from the mantle melt column beneath the back-arc is added to the source of arc magmas. We have identified that the degree of arc/back-arc coupling displayed in the rock record provides an index of the depth of hydrous melting beneath the arc. In this case, this would imply a depth of ~ 75 km for generation of arc magmas, indicating that hydrous melting in the mantle wedge is triggered by the breakdown of hydrous phases in the subducting slab.     

Modeling fluvial response to in‐stream woody vegetation: implications for stream corridor restoration

River restoration and bank stabilization programs often use vegetation for improving stream corridor habitat, aesthetic and function. Yet no study has examined the use of managed vegetation plantings to transform a straight, degraded stream corridor into an ecologically functional meandering channel. Experimental data collected using a distorted Froude‐scaled flume analysis show that channel expansion and widening, thalweg meandering and riffle and pool development are possible using discrete plantings of rigid, emergent vegetation, and the magnitudes of these adjustments depend on the shape of the vegetation zone and the density of the vegetation. These experimental results were verified and validated using a recently developed numerical model, and model output was then used to discuss mechanistically how rivers respond to the introduction of in‐stream woody vegetation. Finally, a hybrid method of meander design is proposed herein where managed vegetation plantings are used to trigger or force the desired morphologic response, transforming a straight, degraded reach into a more functional meandering corridor. It is envisioned that such numerical models could become the primary tool for designing future stream restoration programs involving vegetation and assessing the long‐term stability of such activities. Copyright © 2007 John Wiley & Sons, Ltd.     

Validity of the Bouchet’s complementary relationship at 102 observatories across China

Bouchet in 1963 hypothesized that for large homogeneous land surface with minimum advection of heat and moisture, there exists a 1:1 complementary relationship of potential and actual evaporation coupled through land-atmosphere feedbacks. The complementary relationship has been widely used to estimate regional actual evaporation and explain the pan evaporation paradox. We examine the standardized potential evaporation (potential evaporation divided by wet environment evaporation) at 102 observatories at different elevations across China. Generally, the relationship is appropriate at the low elevations (<1000 m). With the increase of elevation, vapor transfer power becomes much less than radiation energy budget because of lower vapor pressure deficit and stronger global solar radiation. As a result, at the high elevations (over 1000 m), the excess energy resulted by limited moisture availability is not enough to be converted into drying power of the air. This result suggests that the complementary relationship is asymmetric at the high elevations. Supported by the Presidential Special Award Foundation, the Chinese Academy of Sciences (Grant No. O7R70020SD) and the National Key Technology R & D Program (Grant No. 2006BAC08B0408)     

COMPUTER MODELING OF VERTICAL SEISMIC PROFILING*

There are numerous modeling techniques commonly employed for the computer simulation of seismic wave propagation. The capabilities of these techniques vary according to the theoretical foundations and subsequent approximations upon which the algorithms are based. This paper constitutes a comparative review of seven modeling methods: geometric ray theory, asymptotic ray theory, generalized ray theory, Kirchhoff wave theory, Fourier synthesis, finite differences, and finite elements. These methods can be categorized as ray or wave, acoustic or elastic, and can be contrasted according to their relative abilities to simulate such behavior as wave interference effects, diffractions, and mode conversions. As is implied by their names, geometric ray theory and asymptotic ray theory are both ray methods. The other five methods provide wave theory simulations. Geometric ray theory and Kirchhoff wave theory are normally implemented in acoustic form, while the other methods are readily adapted for computing elastic theory solutions. Generalized ray theory and Fourier synthesis are more limited in the complexity of geological model they can accommodate than are the other techniques. The methods which typically demand the greatest computer resources are the finite-difference and the finite-element techniques. All methods can incorporate at least some multiple events. Diffractions, however, are only inherent in the solutions computed by Kirchhoff wave theory, finite differences and finite elements. Attenuation is readily incorporated in both the Fourier synthesis and the finite-element methods. As an example of the application of seismic modeling, a geological model representative of a typical petroleum exploration target is used to compare vertical seismic profiles calculated by different modeling methods.     

A statistical study of pressure changes in the troposphere and lower stratosphere after strong solar flares

Summary The average change in height of several constant pressure levels in the troposphere and lower stratosphere after a strong solar flare is described. The analysis covers the northern hemisphere north of 10°N and is based on a sample of 81 carefully selected flares from the period July 1957 through December 1959.The statistical significance of the results is tested by drawing a comparison with the results obtained when one applies exactly the same analysis to three samples of 81 random key-dates selected from the same period.Properties of the reaction pattern are (1) it is pronounced at high as well as low latitudes, (2) it consists of alternating cells of positive and negative height change, (3) it is established within six hours after the flare, and (4) the maximum response occurs near the tropopause.It appears more likely that the atmospheric reaction pattern can be attributed to very energetic solar particles rather than to enhanced ultraviolet radiation.     

豫01井水位异常与地震的关系分析

针对2018年青海德令哈台尕海井水位破年变异常,系统梳理和分析该井相关资料,发现该井两套水位仪运行稳定、数据一致性非常好,而这段时期异常与降雨和气压不相关,以及井孔附近地区的视应力值较高,表明该地区地壳应力水平进一步增强。利用水位变幅、去趋势和一阶差分分析,结果显示,该项异常信度较高,反映该区域构造活动信息,是周边两次中等强度地震的共同前兆,属于中短期异常。     

高分辨率地震成像研究--21世纪地震学发展的一个重要趋势

随着宽频带地震观测技术的不断发展,全球及区域数字地震观测台网的密度正在不断加大,用于流动地震观测的宽频带地震仪的数量将继续迅速增加。与此同时,宽频带数字地震观测资料的积累速率按指数增长,这意味着震源波裂过程及地球内部结构的高分辨率地震成像研究将成为２１世纪地震学发展的重要趋势。地震学与地球动力学研究的关系将日趋紧密,宽频带流动地震观测的作用和重要性必将与日俱增。     

Statistical assessment of seismicity patterns in Italy: Are they precursors of subsequent events?

The variations of seismicity rate in Central Apenninesprior to the sequence started in September, 1997 (at00:33 UTC, M _L5.6) has been analysedby statistical methods, with the purpose of pointingup eventual periods of quiescence. The analysis wascarried out on the instrumental catalogue of theIstituto Nazionale di Geofisica (ING), covering theperiod from January 1975 to March 1998. In apreliminary phase, the catalogue was declustered usingthe Reasenberg algorithm. After that, eventualmagnitude shifts due to variations in the modalitiesof observation have been individuated and corrected.The subsequent analysis, carried out making use of theZmap software package, has put in evidence thatthe sequence of September 1997 was preceded bya 2.5 year period characterised by absence of eventsof magnitude larger than 3.2, in an area approximately20 × 40 km wide, including the epicentre of themain shock. The statistical methodology shows thatonly 1/10³ of the space-time volumes analysed inthis study, exhibited quiescence of the same level.The study of seismicity rate change correlated toprevious main shocks in a larger area of CentralApennines shows that none of them were preceded by aseismic quiescence, specially close to the epicentreof the main shock, and lasting until the time ofoccurrence of the main shock as in the 1997 case.Actually, we found other patterns of precursoryquiescence with different time or space distribution.We conclude that precursory quiescence is a realfeature of Central Apennines seismicity, but it isdifficult to define a simple hypothesis, which appliesto the generality of cases and can be tested beforeimplementation in a system of earthquake riskmitigation.