Shear wave anisotropy in the crust around the San Andreas fault near Parkfield: spatial and temporal analysis

We systematically analysed shear wave splitting (SWS) for seismic data observed at a temporary array and two permanent networks around the San Andreas Fault (SAF) Observatory at Depth. The purpose was to investigate the spatial distribution of crustal shear wave anisotropy around the SAF in this segment and its temporal behaviour in relation to the occurrence of the 2004 Parkfield M 6.0 earthquake. The dense coverage of the networks, the accurate locations of earthquakes and the high-resolution velocity model provide a unique opportunity to investigate anisotropy in detail around the SAF zone. The results show that the primary fast polarization directions (PDs) in the region including the SAF zone and the northeast side of the fault are NW–SE, nearly parallel or subparallel to the SAF strike. Some measurements on the southwest side of the fault are oriented to the NNE–SSW direction, approximately parallel to the direction of local maximum horizontal compressive stress. There are also a few areas in which the observed fast PDs do not fit into this general pattern. The strong spatial variations in both the measured fast PDs and time delays reveal the extreme complexity of shear wave anisotropy in the area. The top 2–3 km of the crust appears to contribute the most to the observed time delays; however substantial anisotropy could extend to as deep as 7–8 km in the region. The average time delay in the region is about 0.06 s. We also analysed temporal patterns of SWS parameters in a nearly 4-yr period around the 2004 Parkfield main shock based on similar events. The results show that there are no appreciable precursory, coseismic, or post-seismic temporal changes of SWS in a region near the rupture of an M 6.0 earthquake, about 15 km away from its epicentre.     

Simulating three-dimensional seismograms in 2.5-dimensional structures by combining two-dimensional finite difference modelling and ray tracing

Finite difference (FD) simulation of elastic wave propagation is an important tool in geophysical research. As large-scale 3-D simulations are only feasible on supercomputers or clusters, and even then the simulations are limited to long periods compared to the model size, 2-D FD simulations are widespread. Whereas in generally 3-D heterogeneous structures it is not possible to infer the correct amplitude and waveform from 2-D simulations, in 2.5-D heterogeneous structures some inferences are possible. In particular, Vidale & Helmberger developed an approach that simulates 3-D waveforms using 2-D FD experiments only. However, their method requires a special FD source implementation technique that is based on a source definition which is not any longer used in nowadays FD codes. In this paper, we derive a conversion between 2-D and 3-D Green tensors that allows us to simulate 3-D displacement seismograms using 2-D FD simulations and the actual ray path determined in the geometrical optic limit. We give the conversion for a source of a certain seismic moment that is implemented by incrementing the components of the stress tensor.
Therefore, we present a hybrid modelling procedure involving 2-D FD and kinematic ray-tracing techniques. The applicability is demonstrated by numerical experiments of elastic wave propagation for models of different complexity.     

Broad-band calibration of marine seismic sources used by R/V Polarstern for academic research in polar regions

Air guns and air-gun arrays of different volumes are used for scientific seismic surveys with R/V Polarstern in polar regions. To assess the potential risk of these research activities on marine mammal populations, knowledge of the sound pressure field of the seismic sources is essential. Therefore, a broad-band (0–80 kHz) calibration study was conducted at the Heggernes Acoustic Range, Norway. A GI (2.4 l), a G (8.5 l) and a Bolt gun (32.8 l) were deployed as single sources, 3 GI (7.4 l), 3 G (25.6 l) and 8 VLF™ Prakla-Seismos air guns (24.0 l) as arrays. Each configuration was fired along a line of 3–4 km length running between two hydrophone chains with receivers in 35, 100, 198 and 263 m depth. Peak-to-peak, zero-to-peak, rms and sound exposure levels (SEL) were analysed as functions of range. They show the typical dipole-like directivity of marine seismic sources with amplitude cancellation close to the sea surface, higher amplitudes in greater depths, and sound pressure levels which continuously decrease with range. Levels recorded during the approach are lower than during the departure indicating a shadowing effect of Polarsterns's hull. Backcalculated zero-to-peak source levels range from 224–240 dB re 1 μPa @ 1 m. Spectral source levels are highest below 100 Hz and amount to 182–194 dB re 1 μPa Hz^–1. They drop off continuously with range and frequency. At 1 kHz they are ∼30 dB, at 80 kHz ∼60 dB lower than the peak level. Above 1 kHz amplitude spectra are dominated by Polarstern's self-noise. From the rms and sound exposure levels of the deepest hydrophone radii for different thresholds are derived. For a 180 dB rms-level threshold radii maximally vary between 200 and 600 m, for a 186 dB SEL threshold between 50 and 300 m.     

基于波动方程预测和双曲Radon变换联合压制表面多次波

基于波动方程预测的表面多次波压制方法可处理复杂地下介质的地震资料,但计算成本较高.基于滤波的多次波压制方法计算效率较高,但其成功应用仅局限于一次波和多次波有明显时差差别的地震数据,对来自速度逆转等复杂介质数据则较难获得满意的压制效果.本文将波动方程预测的反馈迭代法和滤波法有效结合,采用GPU(图形处理器)和CPU协同并行加速计算粗略预测表面多次波,随后在双曲Radon域比较分析原始数据和预测的多次波,设计合理有效的Butterworth型自适应滤波器,滤出原始数据Radon域中的多次波能量,进行Radon反变换后,在时空域将多次波从原始数据中减去,得多次波压制结果.文中对理论模拟的单炮数据、复杂的SMAART模型以及实际地震数据进行了计算,结果表明,结合基于波动方程预测和双曲Radon变换的方法有效突破了两种方法各自的局限性,可高效高精度地压制复杂地下介质的表面多次波.     

A central difference method with low numerical dispersion for solving the scalar wave equation

In this paper, we propose a nearly‐analytic central difference method, which is an improved version of the central difference method. The new method is fourth‐order accurate with respect to both space and time but uses only three grid points in spatial directions. The stability criteria and numerical dispersion for the new scheme are analysed in detail. We also apply the nearly‐analytic central difference method to 1D and 2D cases to compute synthetic seismograms. For comparison, the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method are used to model acoustic wavefields. Numerical results indicate that the nearly‐analytic central difference method can be used to solve large‐scale problems because it effectively suppresses numerical dispersion caused by discretizing the scalar wave equation when too coarse grids are used. Meanwhile, numerical results show that the minimum sampling rate of the nearly‐analytic central difference method is about 2.5 points per minimal wavelength for eliminating numerical dispersion, resulting that the nearly‐analytic central difference method can save greatly both computational costs and storage space as contrasted to other high‐order finite‐difference methods such as the fourth‐order Lax‐Wendroff correction scheme and the fourth‐order staggered‐grid finite‐difference method.     

Research Note: Internal/external seismic source wavefield separation and cancellation

A compact representation is obtained for the separation of a scalar wavefield on a closed surface into parts due to internal and external sources. The formula assumes that the total field and its gradients are known on the surface, as is the exact Green function of the medium. The derivation involves four rather straightforward applications of Green’s theorem or the representation theorem, though it is a remarkable result in that waves from either source that traverse the boundary many times are appropriately separated. The intermediate results at the four steps of the derivation also shed light on the possibility of acoustic shielding from unwanted sources without knowledge of the Green function for the medium.     

南海障碍层的季节变化及其与海面通量的关系

对气候平均态的温度和盐度分析证实,夏、秋季南海南部上层海洋存在障碍层现象.决定障碍层出现的关键因素是净淡水通量,而障碍层的深度和厚度分布则受风场的显著影响,障碍层受混合层和等温层发展的调制.夏季,南海南部上层海洋在东南向Ek-man水平输运以及东侧下降运动双重因素作用下,较淡的水体在南海东南侧堆积,混合层底部高温水脱离混合层保留在等温层中,造成障碍层在南海东南侧最为深厚,达到30 m.最厚障碍层出现的位置和最厚"南海暖水"出现的位置几乎重合,障碍层的"热障"作用促进了"南海暖水"的发展.     

污水人工湿地处理对滨海生态系统修复研究进展

我国滨海生态系统受到水污染和生态破坏的双重压力。滨海区污水人工湿地处理既能有效去除BOD、氮磷等污染物,又能通过形成淡水帷幕控制海水入侵危害,实现了中水资源的综合利用和滨海生态系统的生态修复。     

河北及邻区平均波速比变化特征分析

用和达法计算河北及邻区地震波速比并按照不同震级不同窗口作滑动平均。通过对唐山地区波速比变化特征与中强地震对应关系的定性分析,发现M_L≥1.5地震波速比曲线与M_L≥4.0的地震对应较好,提出异常指标进行定量统计得到该区波速比异常与地震的对应关系和预报效能;按照同样的思路对河北另外3个地区进行分析,发现各分区波速比异常既有共性特征也有差异。