Modelling and inversion of single-ended refraction data from the shot gathers of multifold deep seismic reflection profiling—an approach for deriving the shallow velocity structure

A multifold crustal-scale deep seismic near-vertical reflection profile generates a large number of single-ended shot gathers, which provide redundant data sets because of overlapping coverage of the shallow refractors. We present an approach for deriving the shallow velocity structure by modelling and inversion of single-ended seismic refraction first arrival traveltime data. We apply this method to a data set acquired with a 12-km long spread with 100 m spacing of shots and receivers, of the Neoproterozoic Marwar basin in the NW Indian shield. The approach is shown to be quite successful for delineating the shallow refractor depths, steep dips and velocities, even in the absence of regular reverse refraction profiles. The study reveals two-layered sedimentary formations, Malani volcanics and a complicated basement configuration of the Marwar basin, and provides a measure of resolution and uncertainty of the estimated model parameters. A seismic section of the near-trace gather is found to be qualitatively consistent with the derived structural features of the basin. The relative highs and lows, observed in the Bouguer gravity profile, further corroborate the derived velocity model. The present approach can be especially useful in offshore areas and elsewhere, where the single-ended multifold seismic profiles are the only available data sets.     

New insights into the P- and S-wave velocity structure of the D" discontinuity beneath the Cocos plate

Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary.     

Focal mechanisms of recent earthquakes in the Southern Korean Peninsula

We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( M_L = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ₁) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σ_Hmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ₂) and minimum (σ₃) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ₂–σ₃ plane, suggesting a permutation of σ₂ and σ₃. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.     

The radio frequency dielectric properties of the stratified UG1–UG2 geological unit in the Bushveld Complex

A new measurement technique enables the complex dielectric properties of the geological strata comprising the UG1–UG2 (Upper Group 1–Upper Group 2) unit of the Bushveld Complex in South Africa to be determined with unprecedented detail at radio frequencies (RF). Results of non-destructive laboratory measurements of representative diamond drill core samples from the UG1–UG2 unit are presented at 25 MHz. These data establish that the UG1 and UG2 chromitite layers are embedded in rock strata (norite, pyroxenite and anorthosite) which are translucent in the HF spectral band, whereas the chromitite layers themselves exhibit significant velocity contrast, making them good radar reflectors. The data presented here is useful for calibration of the radar system, and for predicting the range and resolution performance of borehole radars operating in both the hanging and footwalls of the economically important platiniferous UG2 reef.     

基于线弹性位错模型反演1997年西藏玛尼Mw7.5级地震的干涉雷达同震形变场——Ⅱ滑动分布反演

1997年11月8日西藏M_w7.5级玛尼地震是干涉雷达技术应用于地震观测以来的一次重要事件.在第一部分中，我们应用广泛使用的Okada线弹性位错模型，假设断层的各个分段滑动量均匀，反演得到断层各个分段的几何参数和均匀滑动量.本部分的反演进一步去除滑动均匀假设，并利用更能反映断层真实状态的角形元位错模型（线弹性），在第一部分反演得到断层几何的基础上，反演断层面的静态位错分布.反演结果表明，线弹性滑动分布模型能够更好地解释观测数据，进一步提高反演的数据拟合程度.最终得到了断层面上的走滑和倾滑位错分布.首次得到的断层面滑动分布显示断层面滑动在浅部（0～12 km）比较集中，地震破裂长度约170 km，最大左旋走滑位移达4.8 m；反演结果还表明局部段落存在较大倾滑位移，量值达到1.9 m，这在断层模型中是不能忽略的，它可能是断层两侧形变不对称的原因之一；反演得到的标量地震矩为2.18×10²⁰ N·m，相当于矩震级M_w7.5，与Velasco等利用地震波形反演得到的结果一致.     

NUMERICAL STUDY ON THE MIXING OF UNSORTED SEDIMENT PARTICLES DISCHARGED INTO A CROSS-FLOW BY MULTIPHASE PARTICLE-IN-CELL （MP-PIC） METHOD

The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.     

Automatic Interpretation of Magnetic Data Using Euler Deconvolution with Nonlinear Background

The voluminous gravity and magnetic data sets demand automatic interpretation techniques like Naudy, Euler and Werner deconvolution. Of these techniques, the Euler deconvolution has become a popular choice because the method assumes no particular geological model. However, the conventional approach to solving Euler equation requires tentative values of the structural index preventing it from being fully automatic and assumes a constant background that can be easily violated if the singular points are close to each other. We propose a possible solution to these problems by simultaneously estimating the source location, depth and structural index assuming nonlinear background. The Euler equation is solved in a nonlinear fashion using the optimization technique like conjugate gradient. This technique is applied to a published synthetic data set where the magnetic anomalies were modeled for a complex assemblage of simple magnetic bodies. The results for close by singular points are superior to those obtained by assuming linear background. We also applied the technique to a magnetic data set collected along the western continental margin of India. The results are in agreement with the regional magnetic interpretation and the bathymetric expressions.     

Transient deformation following the 30 January 1997 dike intrusion at Kīlauea volcano,Hawai’i

On 30 January 1997 an intrusion on Kīlauea volcano opened a new fissure within the East Rift Zone (ERZ) at Nāpau Crater, 3 km uprift from the ongoing eruptions at Pu’u ’Ō’ō. The fissure eruption lasted 22 h and opened a 5.1 km long, nearly vertical dike 1.9 m, extending from the surface to a depth of 2.4 km (Owen et al. 2000b). During the eruption, the lava pond at Pu’u ’Ō’ō drained, and eruptions ceased there. Pu’u ’Ō’ō eventually refilled in late February and eruptions resumed there on 28 March 1997. Continuous GPS data show a large transient following the 30 January 1997 dike intrusion. After lengthening 40 cm during the initial eruption, the baseline between two stations spanning the ERZ lengthened an additional 10 cm over the following 6 months. A coastal station KAEP also exhibited transient deformation, as it continued to move southward (5 cm) over the same 6-month period. The baseline between two stations spanning Kīlauea’s summit caldera contracted sharply during the eruption, but gradually recovered to slightly longer than its previous length 2 months after the intrusion. We use the extended network inversion filter (McGuire and Segall 2003) to invert continuous GPS data for volume change of a spherical pressure source under Kīlauea’s summit, opening distribution on a nearly vertical dike in the ERZ and potential slip on a decollement 9 km beneath the south flank. Following the 30 January intrusion, rift extension continued below the initial dike intrusion for the duration of the transient. Decollement slip, regardless of its assumed depth, is not required to fit the data. The modeled transient summit reinflation and rift opening patterns under Nāpau crater coincide with changes in observed behavior of Pu’u ’Ō’ō’s lava pond. Rift opening accelerated while Pu’u ’Ō’ō eruptions paused and began to decelerate after the lava pond reappeared nearly a month after the Nāpau eruption. The transient deformation is interpreted as resulting from shallow accommodation of the new dike volume.     

Source mechanism and source parameters of May 28, 1998 earthquake,Egypt

On May 28, 1998, a moderate size earthquake of mb 5.5 occurred offshore the northwestern part of Egypt (latitude 31.45°N and longitude 27.64°E). It was widely felt in the northern part of Egypt. Being the largest well-recorded event in the area for which seismic data from the global digital network are available, it provides an excellent opportunity to study the tectonic process and present day stress field occurring along the offshore Egyptian coast. The source parameters of this event are determined using three different techniques: modeling of surface wave spectral amplitudes, regional waveform inversion, and teleseismic body waveform inversion. The results show a high-angle reverse fault mechanism generally trending NNW–SSE. The P-axis trends ENE–WSW consistently with the prevailed compression stress along the southeastern Hellenic arc and southwestern part of the Cyprean arc. This unexpected mechanism is most probably related to a positive inversion of the NW trending offshore normal faults and confirms an extension of the back thrusting effects towards the African margin. The estimated focal depth ranges from 22 to 25 km, indicating a lower crustal origin earthquake owing to deep-seated tectonics. The source time function indicates a single source with rise time and total rupture duration of 2 and 5 s, respectively. The seismic moment (M _o) and the moment magnitude (M _w) determined by the three techniques are 1.03 × 10¹⁷ Nm, 5.28; 1.24 × 10¹⁷ Nm, 5.33; and 1.68 × 10¹⁷ Nm, 5.42; respectively. The calculated fault radius, stress drop, and the average dislocation assuming a circular fault model are 7.2 km, 0.63 Mpa, and 0.11 m, respectively.