Tunneling of seismic body waves through thin high-velocity layers in complex structures

Summary The hybrid ray-reflectivity method is applied to the problem of the transmission of the reflected wave field through a thin high-velocity layer (or through a thin stack of high velocity layers), situated in the overburden of the reflector. In the hybrid ray-reflectivity method, the standard ray method is applied in the smooth parts of the model, and the reflectivity method is used locally at the thin high-velocity layer. With the exception of small epicentral distances, the standard ray method itself fails in such computations. The reason is that a considerable part of the energy for overcritical angles of incidence may be tunneled through the thin high-velocity layer along complex ray-paths, corresponding to inhomogeneous waves. The reflectivity method, applied locally at the thin high-velocity layer, automatically includes all inhomogeneous wave contributions. Thus, the hybrid ray-reflectivity method removes fully the limitations of the standard ray method, but still retains its main advantages, such as its applicability to 2-D and 3-D complex layered structures, flexibility, and low-cost computations. In the numerical examples, the hybrid ray-reflectivity synthetic seismograms are compared with standard ray synthetic seismograms and with full reflectivity computations. The numerical examples show that the hybrid ray-reflectivity method describes the tunneling of seismic energy through a thin high-velocity layer with sufficient accuracy.     

Simple method for combined studies of macroseismic intensities and focal mechanisms

A simple method is presented for the computation of theoretical models of the macroseismic field, approximately valid close to the epicentre of a weak crustal earthquake. It is assumed that the intensity is logarithmically proportional to the energy flux of a complete directS wave. A circular source is used, whose energy-flux directivity is weak and thus simply predictable. The focal mechanism influences the solution through standard far-field double-couple radiation patterns. For the wave propagation in the layered crust the ray method is used, and a simple absorption correction is applied. Conversion coefficients at the earth's surface are included. To speed up repeated computations of the theoretical macroseismic fields for varying focal mechanisms, the ray quantities are computed (and stored) separately. This makes the program fast and simple enough even for routine applications on small microcomputers, whenever observed macroseismic fields, focal mechanisms, and hypocentre locations need joint interpretation.     

The contribution of magnetospheric currents toSq

Since the discovery of the magnetosphere, it has been known that the currents flowing in the magnetosphere contribute toSq, the regular daily variation in the earth's surface magnetic field. The early models, however, were not very accurate in the vicinity of the earth. The magnetospheric contribution toSq has therefore been recalculated by direct integration over the three major magnetospheric current systems; magnetopause, tail and ring. The finite electrical conductivity of the earth, which increases the horizontal and decreases the vertical components of the magnetospheric field at the earth's surface, has been taken into account. The magnetospheric currents are found to contribute 12 nanotesla to the day to night difference in the mid-latitudeSq pattern for steady, quiet magnetospheric conditions. They also contribute to the annual variation in the surface field and must be considered an important source of the observed day to day variation in theSq pattern.     

Identification and quantification of mixed sources of oil spills based on distributions and isotope profiles of long-chain n-alkanes

Combined with quantitative determination of concentration and isotopic composition of petroleum hydrocarbons, weathering simulation experiments on artificially mixed oils and their two end-member oils are performed for identification and quantification of mixed sources. The >C₁₈n-alkanes show no appreciable losses during a short-term weathering process. An approach based on distribution of long-chain n-alkanes (>C₁₈) is suggested for estimating the contribution proportion of each source in mixed oils. Stable carbon isotope profile of individual n-alkanes is a powerful tool to differentiate sources of oil spills, but unavailable to accurately allocate each contribution due to a relatively large analytical error.     

Point source radiation in inhomogeneous anisotropic structures

The ray formulae for the radiation from point sources in unbounded inhomogeneous isotropic as well as anisotropic media consist of two factors. The first one depends fully on the type and orientation of the source and on the parameters of the medium at the source. We call this factor the directivity function. The second factor depends on the parameters of the medium surrounding the source and this factor is the well-known geometrical spreading. The displacement vector and the radiation pattern defined as a modulus of the amplitude of the displacement vector measured on a unit sphere around the source are both proportional to the ratio of the directivity function and the geometrical spreading.For several reasons it is desirable to separate the two mentioned factors. For example, there are methods in exploration seismics, which separate the effects of the geometrical spreading from the observed wave field (so-called true amplitude concept) and thus require the proposed separation. The separation also has an important impact on computer time savings in modeling seismic wave fields generated by point sources by the ray method. For a given position in a given model, it is sufficient to calculate the geometrical spreading only once. A multitude of various types of point sources with a different orientation can then be calculated at negligible additional cost.In numerical examples we show the effects of anisotropy on the geometrical spreading, the directivity and the radiation pattern. Ray synthetic seismograms due to a point source positioned in an anisotropic medium are also presented and compared with seismograms for an isotropic medium.     

Macroalgae blooms and δN in subtropical coastal lagoons from the Southeastern Gulf of California: Discrimination among agricultural, shrimp farm and sewage effluents

Macroalgae blooms of Gracilaria vermiculophylla, Hypnea spinella and Spyridia filamentosa have been found in coastal lagoons in the SE Gulf of California. Agriculture, livestock, shrimp and poultry farms and sewage contribute anthropogenic nitrogen to the systems. The δ¹⁵N of these sources, water column and macroalgae were studied in order to identify the N supply for macroalgae blooms. δ¹⁵N of three species of macroalgae (4.3-13.6‰) were enriched compared to the water column ( 3.7-6.8‰), probably because of fractioning from the macroalgae. δ¹⁵N of POM (1.4-10.3‰) was similar to the water column but the relationship was unclear. Depending on the site, macroalgae showed different δ¹⁵N values since some sites receive more or less influence from one given source of the associated watershed, which is reflected in the different δ¹⁵N values of the macroalgae of the same system and in the relative contributions of the sources.     

Grundlagen zur numerischen Behandlung des allgemeinen Zweischichtenproblems in der Geoelektrik

Zusammenfassung Es wird das allgemeine Zweischichtenproblem in der Geoelektrik neu bearbeitet, und im besonderen werden zwei umfassende Folgen von Spezialfällen einer numerischen Behandlung zugänglich gemacht.

---

Summary A contribution to a convenient numerical calculation of the field of a point source of current situated on the earth's surface above an inclined plane is given.

     

Derivatives of reflection point coordinates with respect to model parameters

The motivation for this paper is to provide expressions for first-order partial derivatives of reflection point coordinates, taken with respect to model parameters. Such derivatives are expected to be useful for processes dealing with the problem of estimating velocities for depth migration of seismic data.The subject of the paper is a particular aspect of ray perturbation theory, where observed parameters—two-way reflection time and horizontal components of slowness, are constraining the ray path when parameters of the reference velocity model are perturbed. The methodology described here is applicable to general rays in a 3D isotropic, heterogeneous medium. Each ray is divided into a shot ray and a receiver ray, i.e., the ray portions between the shot/receiver and the reflection point, respectively. Furthermore, by freezing the initial horizontal slowness of these subrays as the model is perturbed,elementary perturbation quantities may be obtained, comprising derivatives of ray hit positions within theisochrone tangent plane, as well as corresponding time derivatives. The elementary quantities may be estimated numerically, by use of ray perturbation theory, or in some cases, analytically. In particular, when the layer above the reflection point is homogeneous, explicit formulas can be derived. When the elementary quantities are known,reflection point derivatives can be obtained efficiently from a set of linear expressions.The method is applicable for a common shot, receiver or offset data sorting. For these gather types, reflection point perturbationlaterally with respect to the isochrone is essentially different. However, in theperpendicular direction, a first-order perturbation is shown to beindependent of gather type.To evaluate the theory, reflection point derivatives were estimated analytically and numerically. I also compared first-order approximations to true reflection point curves, obtained by retracing rays for a number of model perturbations. The results are promising, especially with respect to applications in sensitivity analysis for prestack depth migration and in velocity model updating.     

The information content of high-frequency seismograms and the near-surface geologic structure of “hard rock” recording sites

Due to hardware developments in the last decade, the high-frequency end of the frequency band of seismic waves analyzed for source mechanisms has been extended into the audio-frequency range (>20 Hz). In principle, the short wavelengths corresponding to these frequencies can provide information about the details of seismic sources, but in fact, much of the signal is the site response of the nearsurface. Several examples of waveform data recorded at hard rock sites, which are generally assumed to have a flat transfer function, are presented to demonstrate the severe signal distortions, includingf _max, produced by near-surface structures. Analysis of the geology of a number of sites indicates that the overall attenuation of high-frequency (>1 Hz) seismic waves is controlled by the whole-path-Q between source and receiver but the presence of distinctf _max site resonance peaks is controlled by the nature of the surface layer and the underlying near-surface structure. Models of vertical decoupling of the surface and nearsurface and horizontal decoupling of adjacent sites on hard rock outcrops are proposed and their behaviour is compared to the observations of hard rock site response. The upper bound to the frequency band of the seismic waves that contain significant source information which can be deconvolved from a site response or an array response is discussed in terms off _max and the correlation of waveform distortion with the outcrop-scale geologic structure of hard rock sites. It is concluded that although the velocity structures of hard rock sites, unlike those of alluvium sites, allow some audio-frequency seismic energy to propagate to the surface, the resulting signals are a highly distorted, limited subset of the source spectra.     

基于气枪震源信号和背景噪声的2016年呼图壁6.2级地震前后波速变化特征初步研究

本文选取呼图壁气枪发射台周边流动台记录的2016年6月至2017年12月气枪震源信号和背景噪声信号,分别采用互相关时延检测法和移动窗互谱法研究2016年12月呼图壁6.2级地震前后波速变化特征,得到以下初步结果:呼图壁6.2级地震前后,距6.2级地震60 km的流动台记录的气枪震源信号未发现明显的波速变化,距6.2级地...     

Probabilistic Seismic Hazard Assessment for Iraq Using Complete Earthquake Catalogue Files

Probabilistic seismic hazard analysis (PSHA) has been carried out for Iraq. The earthquake catalogue used in the present study covers an area between latitude 29°–38.5° N and longitude 39°–50° E containing more than a thousand events for the period 1905–2000. The entire Iraq region has been divided into thirteen seismogenic sources based on their seismic characteristics, geological setting and tectonic framework. The completeness of the seismicity catalogue has been checked using the method proposed by Stepp (1972). The analysis of completeness shows that the earthquake catalogue is not complete below M_s=4.8 for all of Iraq and seismic source zones S1, S4, S5, and S8, while it varies for the other seismic zones. A statistical treatment of completeness of the data file was carried out in each of the magnitude classes. The Frequency Magnitude Distributions (FMD) for the study area including all seismic source zones were established and the minimum magnitude of complete reporting (M_c) were then estimated. For the entire Iraq the M_c was estimated to be about M_s=4.0 while S11 shows the lowest M_c to be about M_s=3.5 and the highest M_c of about M_s=4.2 was observed for S4. The earthquake activity parameters (activity rate , b value, maximum regional magnitude m_max) as well as the mean return period (R) with a certain lower magnitude m_min m along with their probability of occurrence have been determined for all thirteen seismic source zones of Iraq. The maximum regional magnitude m_max was estimated as 7.87 ± 0.86 for entire Iraq. The return period for magnitude 6.0 is largest for source zone S3 which is estimated to be 705 years while the smallest value is estimated as 9.9 years for all of Iraq.The large variation of the b parameter and the hazard level from zone to zone reflects crustal heterogeneity and the high seismotectonic complexity. The seismic hazard near the source boundaries is directly and strongly affected by the change in the delineation of these boundaries. The forces, through which the geological structure along the plate boundary in Eastern and Northeastern Iraq are evolved, are still active causing stress-strain accumulation, deformation and in turn producing higher probabilities of earthquake activity. Thus, relatively large destructive earthquakes are expected in this region. The study is intended to serve as a reference for more advanced approaches and to pave the path for the probabilistic assessment of seismic hazard in this region.     

鄱阳湖湿地典型植被群落地下水—土壤—植被—大气系统界面水分通量及水源组成

地下水-土壤-植被-大气系统(GSPAC)界面水分传输是湿地生态水文过程研究的关键.本文选取鄱阳湖湿地高位滩地的2种典型植被群落:茵陈蒿(Artemisia capillaris)和芦苇(Phragmites australis)群落为研究对象,运用HYDRUS-1D垂向一维数值模拟,量化了湿地GSPAC系统界面水分通量,阐明了典型丰水年(2012年)和枯水年(2013年)鄱阳湖湿地植被群落的蒸腾用水规律和水源组成.结果表明:(1)茵陈蒿和芦苇群落土壤-大气界面的年降水入渗量为1570~1600 mm,主要集中在雨季4-6月,占年总量的60%;植物-大气界面的年蒸腾总量分别为346~470 mm和926~1057 mm,其中7-8月植被生长旺季最大,占年总量的40%~46%;地下水-根区土壤界面的向上补给水量受不同水文年水位变化的影响显著,地下水年补给量分别为15~513 mm和277~616 mm,主要发生在蒸散发作用强烈和地下水埋深较浅的时段.(2)植被蒸腾用水分为生长初期(4-6月)和生长旺季(7-10月)2个阶段,丰水年植被的整个生长期蒸腾用水充足,枯水年植被生长旺季的蒸腾用水受到严重水分胁迫,实际蒸腾量仅为潜在蒸腾量的一半左右.(3)不同水文年湿地植被生长旺季的水源贡献不同:丰水年茵陈蒿群落以地下水补给为主,芦苇群落以湖水和地下水补给为主;枯水年茵陈蒿群落以降水和前期土壤水储量为主,芦苇群落以地下水补给为主.本研究结果有助于揭示湿地植被的水分利用策略,为阐明湖泊水情变化与植被演替的作用机理提供参考依据.     

1973-2013年红碱淖水域水质变化及驱动力分析

对1973-2013年8期Landsat MSS、TM、ETM⁺、OLI影像进行了辐射定标、大气校正、辐射归一化和波段运算等处理;利用归一化差值植被指数（NDVI）,分期提取了红碱淖水域面积,分析了湖水水质及红碱淖周围植被变化.研究表明：40 a间红碱淖水域面积呈阶段性萎缩趋势,1990s后萎缩速率加剧,水域面积总体缩小45.7%;湖区周边NDVI波动性增加显示了水退草进的变化趋势;湖面NDVI值的骤增,暗示叶绿素a或悬浮物浓度增加,间接表明湖区水质变差.根据40 a来水域面积变化,红碱淖的演变进程可依据湖泊面积动态度划分为稳定期（1973-1994年）和萎缩期（1994-2013年）两个阶段,气候暖干化是影响稳定期（1973-1994年）湖泊变化的主要因素,在萎缩期（1994-2013年）气候暖干化叠加人类活动是湖泊水量减少和水质变差的诱因,高强度的人为干扰如人工筑坝、灌溉耗水和煤炭开采是红碱淖水域面积锐减的主要原因.     

三峡水库香溪河库湾拟多甲藻(Peridiniopsis)的昼夜垂直迁移行为对碳磷分布的响应

于2015年2月26-27日在香溪河库湾高岚河段拟多甲藻(Peridiniopsis)水华发生地昼夜监测其在水体中的垂直迁移,分别在6个时间段(8:00、12:00、16:00、20:00、0:00、4:00)分层取样分析拟多甲藻在水中的叶绿素a浓度、细胞密度、碳酸酐酶(CA)与碱性磷酸酶(AP)活性;分析总磷(TP)、溶解性总磷(DTP)、溶解性反应磷(SRP)与溶解性有机磷(DOP)浓度在水中的垂直分布.结果显示,水华发生期拟多甲藻为绝对优势种,监测的6个时段中,拟多甲藻呈显著的昼夜垂直迁移规律,正午12:00,表层细胞密度最大,可达5.2×10~6cells/L,凌晨0:00表层细胞密度最低,仅为0.4×10~6cells/L;TP与DTP浓度无垂直分布规律,但SRP与DOP浓度呈现显著的垂直分布规律,SRP浓度随水深增加而减小,DOP浓度随水深增加而增加;CA活性日间高于夜间,日间CA活性呈垂直分布规律,随水深增加而减小,但夜间无垂直分布规律;AP活性在水中昼夜呈极其显著的垂直分布规律,随着水深增加逐渐增加.相关性分析结果表明,拟多甲藻细胞密度与水深、DOP浓度呈极显著负相关,水深与CA和AP活性分别呈极其显著负相关与正相关,拟多甲藻细胞密度与CA和AP浓度分别呈极其显著正相关与负相关.因此,拟多甲藻的垂直迁移行为可能是对日间迁移至表层获取光能和无机碳来源、夜间迁移至水下获取磷源的一种适应.     

多层倾斜层状介质P-SV转换波转换点计算的最佳角度搜索法

转换波转换点的位置对转换波道集的抽取和叠加都非常重要.目前对于单一倾斜反射界面或水平层状介质模型,已经能较好地计算转换波转换点的位置.本文针对多层倾斜层状介质,提出最佳角度搜索法求取P-SV转换波转换点的位置,此方法是通过搜索最符合Snell定律的入射角和反射角来确定最佳的P-SV转换波转换点的位置.在搜索最佳转换点位...     

Spatial Scaling of Effective Modulus and Correlation of Deformation Near the Critical Point of Fracturing

Many observations point to the lithosphere being metastable and close to a critical mechanical point. Exercises in modelling deformation, past or present, across subsurface reservoirs need to take account of this criticality in an efficient way. Using a renormalization technique, the spatial scaling of effective elastic modulus is derived for 2-D and 3-D bodies close to the critical point of through-going fracturing. The resulting exponent, d_μ, of spatial scaling of effective modulus with size, , takes the values ~ −2.5 and −4.2 in two- and three-dimensional space, respectively. The exponents are compatible with those for scaling of effective modulus with fracture density near the percolation threshold determined by other workers from numerical experiments; the high absolute values are also approximately consistent with empirical data from a) fluctuations in depth of a seismic surface; b) `1/k' scaling of heterogeneities observed in one-dimensional well-log samples; c) spatial correlation of slip displacements induced by water injection. The effective modulus scaling modifies the spatial correlation of components of displacement or strain for a domain close to the critical point of fracturing. This correlation function has been used to geostatistically interpolate components of the strain tensor across subsurface reservoirs with the prime purpose of predicting fracture densities between drilled wells. Simulations of strain distributions appear realistic and can be conditioned to surface depths and observations at wells of fracture-related information such as densities and orientations, welltest permeabilities, changes in well-test permeabilities, etc.     

Interpretation of magnetic and gravity gradient tensor data using normalized source strength – A case study from McFaulds Lake,Northern Ontario,Canada

In this paper, we present a case study on the use of the normalized source strength (NSS) for interpretation of magnetic and gravity gradient tensors data. This application arises in exploration of nickel, copper and platinum group element (Ni‐Cu‐PGE) deposits in the McFaulds Lake area, Northern Ontario, Canada. In this study, we have used the normalized source strength function derived from recent high resolution aeromagnetic and gravity gradiometry data for locating geological bodies. In our algorithm, we use maxima of the normalized source strength for estimating the horizontal location of the causative body. Then we estimate depth to the source and structural index at that point using the ratio between the normalized source strength and its vertical derivative calculated at two levels; the measurement level and a height h above the measurement level. To discriminate more reliable solutions from spurious ones, we reject solutions with unreasonable estimated structural indices. This method uses an upward continuation filter which reduces the effect of high frequency noise. In the magnetic case, the advantage is that, in general, the normalized magnetic source strength is relatively insensitive to magnetization direction, thus it provides more reliable information than standard techniques when geologic bodies carry remanent magnetization. For dipping gravity sources, the calculated normalized source strength yields a reliable estimate of the source location by peaking right above the top surface. Application of the method on aeromagnetic and gravity gradient tensor data sets from McFaulds Lake area indicates that most of the gravity and magnetic sources are located just beneath a 20 m thick (on average) overburden and delineated magnetic and gravity sources which can be probably approximated by geological contacts and thin dikes, come up to the overburden.