Active tectonic morphology and submarine deformation of the northern Gulf of Eilat/Aqaba from analyses of multibeam data

A high-resolution marine geophysical study was conducted during October-November 2006 in the northern Gulf of Aqaba/Eilat, providing the first multibeam imaging of the seafloor across the entire gulf head spanning both Israeli and Jordanian territorial waters. Analyses of the seafloor morphology show that the gulf head can be subdivided into the Eilat and Aqaba subbasins separated by the north-south-trending Ayla high. The Aqaba submarine basin appears starved of sediment supply, apparently causing erosion and a landward retreat of the shelf edge. Along the eastern border of this subbasin, the shelf is largely absent and its margin is influenced by the Aqaba Fault zone that forms a steep slope partially covered by sedimentary fan deltas from the adjacent ephemeral drainages. The Eilat subbasin, west of the Ayla high, receives a large amount of sediment derived from the extensive drainage basins of the Arava Valley (Wadi ’Arabah) and Yutim River to the north–northeast. These sediments and those entering from canyons on the south-western border of this subbasin are transported to the deep basin by turbidity currents and gravity slides, forming the Arava submarine fan. Large detached blocks and collapsed walls of submarine canyons and the western gulf margin indicate that mass wasting may be triggered by seismic activity. Seafloor lineaments defined by slope gradient analyses suggest that the Eilat Canyon and the boundaries of the Ayla high align along north- to northwest-striking fault systems—the Evrona Fault zone to the west and the Ayla Fault zone to the east. The shelf–slope break that lies along the 100 m isobath in the Eilat subbasin, and shallower (70–80 m isobaths) in the Aqaba subbasin, is offset by approx. 150 m along the eastern edge of the Ayla high. This offset might be the result of horizontal and vertical movements along what we call the Ayla Fault on the east side of the structure. Remnants of two marine terraces at 100 m and approx. 150 m water depths line the southwest margin of the gulf. These terraces are truncated by faulting along their northern end. Fossil coral reefs, which have a similar morphological appearance to the present-day, basin margin reefs, crop out along these deeper submarine terraces and along the shelf–slope break. One fossil reef is exposed on the shelf across the Ayla high at about 60–63 m water depth but is either covered or eroded in the adjacent subbasins. The offshore extension of the Evrona Fault offsets a fossil reef along the shelf and extends south of the canyon to linear fractures on the deep basin floor.     

Origin and role of major strike-slip transfers during plate collision in the central Mediterranean

Variations in the crustal structure along the northern African plate margin have caused different modes of collision with Eurasia. Lateral density variations along the central Mediterranean collision zone are expressed in a change of the angle of the downbending African Plate and lead to the formation of strike-slip transfers in these transition zones that are roughly perpendicular to the trend of the collisional zone. In some cases these transfer zones are developed into hinge faults, while in others they can be developed into transform faults. This process governs the segmentation of the collision zone in the central Mediterranean region south of the Maghrebian thrust belt in Tunisia and Sicily through the Calabrian Arc to the northeastern Hellenic Arc, extending further to the Cyprian Arc and to the Taurus-Zagros chain.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part II: Pre-critical slowness match

Part II of this paper is a direct continuation of Part I, where we consider the same types of orthorhombic layered media and the same types of pure-mode and converted waves. Like in Part I, the approximations for the slowness-domain kinematical characteristics are obtained by combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. In Part I, the wide-angle ray was set to be the critical ray (‘critical slowness match’), whereas in Part II we consider a finite long offset associated with a given pre-critical ray (‘pre-critical slowness match’). Unlike the critical slowness match, the approximations in the pre-critical slowness match are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. Moreover, for the pre-critical slowness match, there is no need to distinguish between the high-velocity layer and the other, low-velocity layers. The form of the approximations in both critical and pre-critical slowness matches is the same, where only the wide-angle power series coefficients are different. Comparing the approximated kinematical characteristics with those obtained by exact numerical ray tracing, we demonstrate high accuracy. Furthermore, we show that for all wave types, the accuracy of the pre-critical slowness match is essentially higher than that of the critical slowness match, even for matching slowness values close to the critical slowness. Both approaches can be valuable for implementation, depending on the target offset range and the nature of the subsurface model. The pre-critical slowness match is more accurate for simulating reflection data with conventional offsets. The critical slowness match can be attractive for models with a dominant high-velocity layer, for simulating, for example, refraction events with very long offsets.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part I: 1D ray propagation‡

Anisotropy in subsurface geological models is primarily caused by two factors: sedimentation in shale/sand layers and fractures. The sedimentation factor is mainly modelled by vertical transverse isotropy (VTI), whereas the fractures are modelled by a horizontal transversely isotropic medium (HTI). In this paper we study hyperbolic and non‐hyperbolic normal reflection moveout for a package of HTI/VTI layers, considering arbitrary azimuthal orientation of the symmetry axis at each HTI layer. We consider a local 1D medium, whose properties change vertically, with flat interfaces between the layers. In this case, the horizontal slowness is preserved; thus, the azimuth of the phase velocity is the same for all layers of the package. In general, however, the azimuth of the ray velocity differs from the azimuth of the phase velocity. The ray azimuth depends on the layer properties and may be different for each layer. In this case, the use of the Dix equation requires projection of the moveout velocity of each layer on the phase plane. We derive an accurate equation for hyperbolic and high‐order terms of the normal moveout, relating the traveltime to the surface offset, or alternatively, to the subsurface reflection angle. We relate the azimuth of the surface offset to its magnitude (or to the reflection angle), considering short and long offsets. We compare the derived approximations with analytical ray tracing.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part I: Critical slowness match

Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy.     

Moveout approximation for horizontal transversely isotropic and vertical transversely isotropic layered medium. Part II: effective model‡

We use residual moveouts measured along continuous full azimuth reflection angle gathers, in order to obtain effective horizontal transversely isotropic model parameters. The angle gathers are generated through a special angle domain imaging system, for a wide range of reflection angles and full range of phase velocity azimuths. The estimation of the effective model parameters is performed in two stages. First, the background horizontal transversely isotropic (HTI)/vertical transversely isotropic (VTI) layered model is used, along with the values of reflection angles, for converting the measured residual moveouts (or traveltime errors) into azimuthally dependent normal moveout (NMO) velocities. Then we apply a digital Fourier transform to convert the NMO velocities into azimuthal wavenumber domain, in order to obtain the effective HTI model parameters: vertical time, vertical compression velocity, Thomsen parameter delta and the azimuth of the medium axis of symmetry. The method also provides a reliability criterion of the HTI assumption. The criterion shows whether the medium possesses the HTI type of symmetry, or whether the azimuthal dependence of the residual traveltime indicates to a more complex azimuthal anisotropy. The effective model used in this approach is defined for a 1D structure with a set of HTI, VTI and isotropic layers (with at least one HTI layer). We describe and analyse the reduction of a multi‐layer structure into an equivalent effective HTI model. The equivalent model yields the same NMO velocity and the same offset azimuth on the Earth's surface as the original layered structure, for any azimuth of the phase velocity. The effective model approximates the kinematics of an HTI/VTI layered structure using only a few parameters. Under the hyperbolic approximation, the proposed effective model is exact.     

Traveltime approximation in vertical transversely isotropic layered media

The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging.     

Volcanic gaps due to oblique consumption of aseismic ridges

The present-day consumption of oceanic ridges and other buoyant rises and fragments at circum-Pacific subduction zones, and presumably elsewhere, are closely related to existing gaps in volcanism. Examples are the gaps associated with the Nazca, Juan Fernandez, Cocos, Marcus-Necker and Louisville ridges. The buoyancy of these ridges breaks the continuity of the subducted plate, which may lead to reduced water supply required for melting of magma, and therefore create temporary volcanic gaps. The oblique consumption of these ridges causes the gap to migrate with time. This mechanism may be useful in interpreting time-space patterns of past volcanic chains associated with subduction in terms of the consumption of the disruptive oceanic plateaus and ridges.     

A case study of sub-basalt imaging in land region covered with basalt flows

In this study, a set of 20 2D seismic lines, acquired over the Golan Heights basaltic plateau, was processed and analysed. Although the data were acquired and processed by standard techniques, in some cases good-quality seismic images were obtained under several hundred metres of basalts. We describe how the seismic characteristics of the top basalt layer were defined and show the effect of the numerous widespread volcanic sources on the quality of the final images. The new data reveal the first images of the sedimentary sequence under the basalt flows, and indicate that strands of the Dead Sea Transform extend into this area. The entire region was found to be very deformed. Several attractive traps for hydrocarbon exploration were also identified on the output sections.