Layer-stripping reverse-time migration1

We present a layer-stripping method of migration for irregularly layered media in which first-order velocity discontinuities separate regions of constant or smoothly varying velocity. We use the reverse-time method to migrate seismic data layer by layer, from the surface downwards. As part of the migration of a given layer, the bottom boundary of the layer is defined based on power in the migrated signal, and a seismic section is collected along it. This new section serves as the boundary condition for migration in the next layer. This procedure is repeated for each layer, with the final image formed from the individual layer images. Layer-stripping migration consists of three steps: (1) layer definition, (2) wavefield extrapolation and imaging, and (3) boundary determination. The migration scheme when used with reverse-time extrapolation is similar to datuming with an imaging condition. The reverse-time method uses an explicit fourth-order time, tenth-order space, finite-difference approximation to the scalar wave equation. The advantages of layer-stripping reverse-time migration are: (1) it preserves the benefits of the reverse-time method by handling strong velocity contrasts between layers and steeply dipping structures; (2) it reduces computer memory and saves computation time in high-velocity layers, and (3) it allows interpretational control of the image. Post-stack layer-stripping reverse-time migration is illustrated with a synthetic CMP data example. Prestack migration is illustrated with a synthetic data set and with a marine seismic reflection profile across the Santa Maria Basin and the Hosgri Fault in central California.     

Mantle transition zone beneath the Caribbean-South American plate boundary and its tectonic implications

We analyzed receiver-function data recorded by a temporary broadband array deployed as part of the BOLIVAR project and the permanent seismic network of Venezuela to study the mantle transition zone structure beneath the Caribbean-South American plate boundary and Venezuela. Significant topography on both the 410-km and the 660-km discontinuities was clearly imaged in the CCP (common-conversion-point) stacked images. Beneath the southeastern Caribbean, the 410-km is featured by a narrow (～ 200 km EW) ～ 25-km uplift extending in the NS direction around 63° west, while the 660-km is depressed by ～ 20 km in a narrow region slightly west to the uplift, a scenario that is more consistent with westward descent of the oceanic South American plate rather than a break-off of NNW dipping proto-Caribbean oceanic lithosphere along the El Pilar Fault. We also found a thick transition zone beneath the Falcon region in northwestern Venezuela, possibly associated with the subducted Nazca plate. A flat 410-km was observed beneath the Guayana shield, suggesting that the shield has a stable and moderately deep keel, which has little effect on the underlying transition zone structure.     

Receiver function imaging in strongly laterally heterogeneous crust: Synthetic modeling of BOLIVAR data

We resolve a large (~20 km) discrepancy in Moho depth determined from PdS receiver functions (RFs) and from active source seismic profiling in the complex Caribbean-South American plate boundary zone in eastern Venezuela. As part of the BOLIVAR experiment 20 broadband stations were deployed along an active source profile to record teleseisms. Using the extremely heterogeneous crustal model obtained from active source data, we generated 2D finite-difference elastic wave synthetics and from them calculated receiver functions and CCP stacks. We compare the observations with synthetic sections that have been spatially sampled at 0.25 km to 40 km. The densely sampled synthetics show that several events in the field data that were originally interpreted as the Moho are multiple reflections within sedimentary basins. Where the Moho has the steepest dip under the plate boundary the CCP stacks fail to image the Moho well, regardless of the density of spatial sampling. A suitable spatial sampling criterion for clearly imaging the lower crust and Moho is to overlap Fresnel zones by 50% at Moho depth, which for the 1 Hz receiver functions examined here, requires an instrument spacing of 15–20 km, with the actual field data density ranging from 20 km to 100 km.     

Composition of widespread volcanic glass in deep-sea sediments of the Southern Pacific Ocean: an Antarctic source inferred

Widespread Plio-Pleistocene (2.43-0.06 Ma) tephra zones recognised in deep-sea cores from high latitudes (>60°) in the Southern Pacific Ocean were thought to have originated from calc-alkaline rhyolitic eruptions in New Zealand, some 5000 km distant. Electron microprobe analyses of the glasses reveal a wide diversity of alkalic felsic compositions, as well as minor components of basic and intermediate glasses, incompatible with a New Zealand Neogene source but similar to contemporaneous eruptives from the Antarctic region. Most tephra zones are trachytic; seven zones are peralkaline rhyolite. The rhyolitic zones represent a deep-sea record of widespread silicic eruptions from continental Antarctica, possibly Marie Byrd Land. The extent of these rhyolitic zones suggest a greater frequency of large explosive eruptions in Antarctica than previously documented. The coarse grain size of some of the shards (up to 3 mm), their great distance from the closest sources (>1600 km for some cores), and the presence of nonvolcanic ice-rafted debris indicate some of the glasses, especially the more basic compositions, may have been ice-rafted, contrary to previous suggestions of a fallout origin.     

Imag(in)ing the continental lithosphere

This paper is primarily concerned with seismically imaging details in the mantle at an intermediate scale length between the large scales of regional and global tomography and the small scales of reflection profiles and outcrops. This range is roughly 0.1–1 km < a < 10–10² km, where a is the scale. We consider the implications of several models for mantle evolution in a convecting mantle, and possible scales present in the non-convecting tectosphere. Reflection seismic evidence shows that the structures preserved from continental accretion within and at the margins of the Archean cratons are subduction related, and we use subduction as an analog for scales left by past events. In modern orogenic belts we expect to find subduction structures, small scale upper mantle convection structures, and basalt extraction structures. We examine some of the scales that are likely formed by orogenic processes.We also examine the seismic velocity and density contrasts expected between various upper mantle constituents, including fertile upper mantle, depleted upper mantle, normal and eclogitized oceanic crust, and fertile mantle with and without partial melt. This leads directly to predicting the size of seismic signals that can be produced by specular conversion, and scattering from layers and objects with these contrasts.We introduce an imaging scheme that makes use of scattered waves in teleseismic receiver functions to make a depth migrated image of a pseudo-scattering coefficient. Image resolution is theoretically at least an order of magnitude better than traveltime tomography. We apply the imaging scheme to three data sets from 1) the Kaapvaal craton, 2) the Cheyenne Belt, a Paleoproterozoic suture between a protocontinent and an island arc, and 3) the Jemez Lineament, a series of aligned modern volcanic structures at the site of a Proterozoic suture zone. The Kaapvaal image, although not defining a unique base of the tectosphere, shows complicated “layered” events in the region defined as the base of the tectosphere in tomography images. The image of the transition zone discontinuities beneath the Kaapvaal craton is remarkable for clarity. The migrated receiver function image of the upper mantle beneath the Cheyenne belt is complicated, more so than the tomography image, and may indicate limitations in the receiver function imaging system. In contrast the Jemez Lineament image shows large-amplitude negative-polarity layered events beneath the Moho to depths of 120 km, that we interpret as melt-containing sills in the upper mantle. These sills presumably feed the Quaternary–Neogene regional basaltic volcanic field.     

Tag loss and movement of stocked yearling rainbow trout (Oncorhynchus mykiss) in southwestern Australia

Rainbow trout, Oncorhynchus mykiss, have been stocked in southwestern Australia since the 1930s. Trials at a research station maintained 868 tagged yearling O. mykiss (age: 14 months; standard length: 208 mm) for up to 375 days. The mean tag loss rate was 0.12% ± 0.040% per day for both single- and double-tagged fish. The tag loss rate was used to interpret tag return data from 7030 single-tagged yearlings stocked into three river systems in southwestern Australia. Recaptures indicated that O. mykiss survived into their second season in all rivers, with the maximum time at liberty of 21 months (longevity of at least 35 months). Upstream and downstream movements of O. mykiss increased with time, although 29% of recaptures were reported from the stocking site. Significantly more O. mykiss were recaptured from the Warren River, and differences remained even when corrected for stocking levels; angler effort had a marginal effect on return rates. Higher return rates from the Warren River may be due to availability of suitable habitat (stream cover, cooler water) for O. mykiss. Results allow the stocking regime to be reviewed to improve fishery performance.     

Significant genetic differentiation among meroplanktonic barrel jellyfish Rhizostoma pulmo (Cnidaria: Scyphozoa) in the Mediterranean Sea

Molecular data have shown that jellyfishes are more geographically restricted and evolutionarily divergent than previously thought. We examined genetic variation and divergence within the meroplanktonic barrel jellyfish Rhizostoma pulmo in the Mediterranean Sea; specific sampling areas were the northern Adriatic, western Mediterranean and Tunisian coast. A total of 19 sampling sites and 68 sequences of the mtDNA cytochrome c oxidase subunit I (COI) gene were used. Of the 68 COI sequences, 45 were newly collected specimens which originated from nine sampling sites along the Tunisian coast. A total of 24 haplotypes were obtained and the specimens sampled were characterised by relatively high levels of haplotype diversity (h = 0.866) and low levels of nucleotide diversity (π = 0.004). Haplotype network analysis showed the presence of three distinct phylogenetic lineages (populations), with separate geographic ranges in the northern Adriatic, western Mediterranean and Tunisian coast. The observed genetic differentiation between these three lineages was supported by the presence of significant genetic differentiation between the 19 populations (F_ST = 0.757, p < 0.001). The high level of genetic differentiation detected in the barrel jellyfish investigated could be attributed to either intrinsic and/or extrinsic barriers to genetic exchange between different populations that may have adapted to different environmental conditions.