Geometric corrections on sidescan sonar images based on bathymetry. Application with SeaMARC II and Sea Beam data

Acoustic backscatter images of the seafloor obtained with sidescan sonar systems are displayed most often using a flat bottom assumption. Whenever this assumption is not valid, pixels are mapped incorrectly in the image frame, yielding distorted representations of the seafloor. Here, such distortions are corrected by using an appropriate representation of the relief, as measured by the sonar that collected the acoustic backscatter information. In addition, all spatial filtering operations required in the pixel relocation process take the sonar geometry into account. Examples of the process are provided by data collected in the Northeastern Pacific over Fieberling Guyot with the SeaMARC II bathymetric sidescan sonar system and the Sea Beam multibeam echo-sounder. The nearly complete (90%) Sea Beam bathymetry coverage of the Guyot serves as a reference to quantify the distortions found in the backscatter images and to evaluate the accuracy of the corrections performed with SeaMARC II bathymetry. As a byproduct, the processed SeaMARC II bathymetry and the Sea Beam bathymetry adapted to the SeaMARC II sonar geometry exhibit a 35m mean-square difference over the entire area surveyed.On leave at the Naval Research Laboratory, Code 7420, Washington D.C. 20375-5350.     

Characteristics of SeaMARC II phase data

The dispersion of SeaMARC II phase-difference samples is discussed. They appear to be a function of signal direction, range, noise level, and backscatter strength of the bottom. Field data from a lava flow area and from sedimented areas at different depths are compared. The temporal distribution of the phase-difference samples was skewed and asymmetrical about the model. The angular distribution was symmetrical about the mode, with some phase wrap-around. The field data show the presence of a complicated noise interference field. The amount of phase-difference dispersion was larger than that calculated by using a simple Gaussian isotropic noise model, possibly suggesting an additional phase-dispersion process caused by bottom roughness. The method used to produce bathymetry data from the phase-difference samples was evaluated in light of the phase-difference sample distribution     

Using inverse theory and seabeam bathymetry to improve seismic refraction data corrections

An array consisting of ocean bottom seismometer and on-bottom hydrophones, was used to conduct a seismic experiment on 0.4 Ma crust east of the Juan de Fuca Ridge. Seismic sources were large (>50 kg) explosive charges detonated by SUS devices set to explode at 1829 or 2438 m nominal depth. The objectives of the experiment were to determine the compressional wave velocity and attenuation structures of the uppermost 500 m depth. The relative positions of shots and receivers were originally determined by treating each shot-receiver pair independently, via raytracing of various water waves. Due to the reflection of some of these water waves by the rough bottom, significant scatter resulted, preventing a determination of a physically realizable velocity-depth function. A new method is described that co-locates shot and receiver positions, including receiver depths consistent withseabeam bathymetry, using only the water waves that do not interact with the bottom. Several potential pitfalls are outlined using this method. A stable solution could only be achieved by discarding shots located well outside the array. The water path corrections were applied to the refracted arrivals, again using theseabeam bathymetry. The joint inversion location procedure, along with the use of precise gridded bathymetry, reduced the travel time scatter to a level whereby a velocity-depth function could be determined. The results, using only the hydrophone data, indicate an initial velocity at the seafloor of 2.7 km s^-1 with gradients from 4.6 s^-1 slowly decreasing to 4.1 s^-1 at 679 m depth. This velocity is similar to others conducted over very young oceanic crust, and can be interpreted as being due to a high porosity at the surface, due to cracks, fissures, and open pores, which rapidly diminish with depth.     

A SeaMARC II survey of recent submarine volcanism near Easter Island

SeaMARC II side-scan sonar data reveal that a large area of seafloor north and west of Easter Island has been disrupted by recent submarine volcanism. A large volcanic area begins approximately 60 km WNW of the island and extends for over 130 km to the west. The volcanic field is characterized by high backscatter intensity in the side-scan sonar records and is elevated 400–1000 m above the N-S seafloor fabric that surrounds it. This field, the Abu Volcanic Field, covers at least 2500 km² and appears to consist of recent lava flows and small volcanoes. Backscatter intensity of the Abu Volcanic Field is similar to that of the adjacent ridge flank which is less than 0.4 Ma, suggesting a similar age for its formation. Two additional areas of high backscatter immediately north of Easter Island cover a combined area of over 300 km². The sidescan sonar records show that these features are clearly of volcanic origin and are not debris flows from the nearby island. The flows are nearly 300 m thick and are morphologically similar to subaerial pahoehoe lava shields. Their high backscatter indicates that they are also the products of relatively recent submarine volcanic activity. The presence of these large areas of recent volcanism in the vicinity of Easter Island has important implications for the various models that have been proposed to explain the origin of the Easter Seamount Chain. In addition, the similar ages of Easter Island and the Easter Microplate suggest that the presence of a hotspot near or beneath this fast-spreading portion of the East Pacific Rise about 4.5 m.y. ago may have initiated the large-scale rift propagation that created the microplate.     

Submarine canyon development in the Izu-Bonin forearc: A SeaMARC II and seismic survey of Aoga Shima Canyon

SeaMARC II sidescan (imagery and bathymetry) and seismic data reveal the morphology, sedimentary processes, and structural controls on submarine canyon development in the central Izu-Bonin forearc, south of Japan. Canyons extend up to 150 km across the forearc from the trench-slope break to the active volcanic arc. The canyons are most deeply incised (1200–1700 m) into the gentle gradients (1–2°) upslope on the outer arc high (OAH) and lose bathymetric expression on the steep (6–18°) inner trench-slope. The drainage patterns indicate that canyons are formed by both headward erosion and downcutting. Headward erosion proceeds on two scales. Initially, pervasive small-scale mass wasting creates curvilinear channels and pinnate drainage patterns. Large-scale slumping, evidenced by abundant crescent-shaped scarps along the walls and tributaries of Aoga Shima Canyon, occurs only after a channel is present, and provides a mechanism for canyon branching. The largest slump has removed >16 km³ of sediment from an 85 km² area of seafloor bounded by scarps more than 200 m high and may be in the initial stages of forming a new canyon branch. The northern branch of Aoga Shima Canyon has eroded upslope to the flanks of the arc volcanoes allowing direct tapping of this volcaniclastic sediment source. Headward erosion of the southern branch is not as advanced but the canyon may capture sediments supplied by unconfined (non-channelized) mass flows.Oligocene forearc sedimentary processes were dominated by unconfined mass flows that created sub-parallel and continuous sedimentary sequences. Pervasive channel cut-and-fill is limited to the Neogene forearc sedimentary sequences which are characterized by migrating and unconformable seismic sequences. Extensive canyon formation permitting sediment bypassing of the forearc by canyon-confined mass flows began in the early Miocene after the basin was filled to the spill points of the OAH. Structural lows in the OAH determined the initial locus of canyon formation, and outcropping basement rocks have prevented canyon incision on the lower slope. A major jog in the canyon axis, linear tributaries, and a prominent sidescan lineament all trend NW-NNW, reflecting OAH basement influence on canyon morphology. This erosional fabric may reflect joint/fracture patterns in the sedimentary strata that follow the basement trends. Once the canyons have eroded down to more erosion-resistant levels, channel downcutting slows relative to lateral erosion of the canyon walls. This accounts for the change from a narrow canyon axis in the thickly sedimented forearc basin to a wider, more rugged canyon morphology near the OAH. About 9500 km³ of sediment has been eroded from the central, 200 km long, segment of the Izu-Bonin forearc by the formation of Aoga Shima, Myojin Sho and Sumisu Jima canyons. The volume of sediment presently residing in the adjacent trench, accretionary wedge, and lower slope terrace basin accounts for <25% of that eroded from the canyons alone. This implies that a large volume (>3500 km³ per 100 km of trench, ignoring sediments input via forearc bypassing) has been subducted beneath the toe of the trench slope and the small accretionary prism. Unless this sediment has been underplated beneath the forearc, it has recycled arc material into the mantle, possibly influencing the composition of arc volcanism.     

Morphology of the Ebro fan valleys from SeaMARC and sea beam profiles

The northern continental slope off the Ebro Delta has a badland topography indicating major slope erosion and mass movement of material that deposits sediment into a ponded lobe. The southern slope has a low degree of mass movement activity and slope valleys feed channel levee-complexes on a steep continental rise. The last active fan valley is V-shaped with little meandering and its thalweg merges downstream with the Valencia Valley. The older and larger inactive channel-levee complex is smoother, U-shaped, and meanders more than the active fan valley.     

SeaMARC II study of a giant submarine slump on the northern Chile continental slope

Abstract

A giant submarine slump, encompassing a 91‐km by 26‐km block, occurring on the continental slope offshore Iquique, Chile, was identified during a SeaMARC II survey. Utilizing SeaMARC II side‐scan imagery, bathymetry, and seismic reflection data, five morphostructural zones of the slump were identified: the fissured zone, scar zone, tensional depression, central block, and front zone. The fissured zone was developed on the crown of the slump; the scar zone is characterized by scars with the crescent‐shaped slip surfaces and throws ranging from 200 m to 50 m. The tensional depression zone is marked by an area voided by mass slumping, while the central block morphology was formed by uplift. The front zone is comprised of both compressional and tensional subzones. The compressional subzone is characterized by a relative topographic low, on the middle slope, whereas the extensional subzone is characterized by a convex pattern of alternated ridges and hollows, which may represent the debris of the slump on the lower slope. The formation of the slump was strongly influenced by the subduction of the Nazca plate beneath the Chile continental margin, which resulted in the subsidence of the continental slope with a resultant increase in the slope gradient and pore‐water pressure in the sedimentary layers. Slump formation was further facilitated by the development of a complex fault system associated with the subduction and by the triggering effect of earthquakes in the area.     

Simultaneous operation of the Sea Beam multibeam echo-sounder andthe SeaMARC II bathymetric sidescan sonar system

An experiment aboard the Scripps Institution of Oceanography's RV Thomas Washington has demonstrated the seafloor mapping advantages to be derived from combining the high-resolution bathymetry of a multibeam echo-sounder with the sidescan acoustic imaging plus wide-swath bathymetry of a shallow-towed bathymetric sidescan sonar. To a void acoustic interference between the ship's 12-kHz Sea Beam multibeam echo-sounder and the 11-12-kHz SeaMARC II bathymetric sidescan sonar system during simultaneous operations, Sea Beam transmit cycles were scheduled around SeaMARC II timing events with a sound source synchronization unit originally developed for concurrent single-channel seismic, Sea Beam, and 3.5-kHz profile operations. The scheduling algorithm implemented for Sea Beam plus SeaMARC II operations is discussed, and the initial results showing their combined seafloor mapping capabilities are presented     

Progressive deformation of an evaporite-bearing accretionary complex: SeaMARC I,SeaBeam and piston-core observations from the Mediterranean Ridge

The Mediterranean Ridge is an arcuate ridge of deformed sediment caught up in the convergent plate margin between the African plate and the Aegean. An intensive campaign of SeaMARC I and SeaBeam surveys followed by piston coring has been conducted along the contact between undeformed turbidites of the Sirte Abyssal Plain and folded and faulted sediments of the Mediterranean Ridge. Along the outer edge of the Ridge, surficial sediments have been deformed into sinusoidal ridges and troughs (wavelengths 0.5–2 km, amplitude 20–150 m), which we interpret as folds. In plan view, the ridge and the trough fabric parallels the NW-SE trending regional contours, suggesting that the folds formed in response to compression orthogonal to the Mediterranean Ridge. The outermost ridge is shedding a debris apron out onto the abyssal plain, implying that uplift and deformation are ongoing. We show that the geometry of the outermost folds can be produced by elastic bending of a packet of 5–10 relatively strong layers, each 10–20 m thick, interbedded between weaker layers; we equate the strong layers with gypsum beds in the Messinian upper evaporites. Folding the seafloor from a flat layer into the observed ridge and trough topography would shorten the layer by less than 2%. Two percent shortening (equals two percent thickening) is insufficient to create the observed relief of the Mediterranean Ridge even if the entire sediment column down to basement were involved; we infer that additional shortening/thickening is accommodated by thrust faulting above a decollement at the top of the Messinian salt layer. At distances > 15 km from the deformation front and more than 500 m from the abyssal plain, sharp-edged, fine-grained side-scan lineations with very little vertical relief cut across the kilometer-scale ridge and trough topography. These fine-grained lineations fall in two groups trending N/S to NNE/SSW and ~ENE. We interpret these lineaments as traces of conjugate strike-slip faults formed in the same compressional regime which formed the NW/SE trending folds. The onset of strike-slip faulting may coincide with the cessation of imbricate thrust fan development above the initial salt-controlled decollement surface. The following characteristics of the Mediterranean Ridge are attributed to the presence of evaporites in the incoming sedimentary section: (1) initial deformation by folding rather than thrust faulting; (2) narrow taper; (3) rapid rate of outward growth; (4) karstification.     

Small-scale sediment transport patterns and bedform morphodynamics: new insights from high-resolution multibeam bathymetry

New multibeam echosounder and processing technologies yield sub-meter-scale bathymetric resolution, revealing striking details of bedform morphology that are shaped by complex boundary-layer flow dynamics at a range of spatial and temporal scales. An inertially aided post processed kinematic (IAPPK) technique generates a smoothed best estimate trajectory (SBET) solution to tie the vessel motion-related effects of each sounding directly to the ellipsoid, significantly reducing artifacts commonly found in multibeam data, increasing point density, and sharpening seafloor features. The new technique was applied to a large bedform field in 20–30 m water depths in central San Francisco Bay, California (USA), revealing bedforms that suggest boundary-layer flow deflection by the crests where 12-m-wavelength, 0.2-m-amplitude bedforms are superimposed on 60-m-wavelength, 1-m-amplitude bedforms, with crests that often were strongly oblique (approaching 90°) to the larger features on the lee side, and near-parallel on the stoss side. During one survey in April 2008, superimposed bedform crests were continuous between the crests of the larger features, indicating that flow detachment in the lee of the larger bedforms is not always a dominant process. Assessment of bedform crest peakedness, asymmetry, and small-scale bedform evolution between surveys indicates the impact of different flow regimes on the entire bedform field. This paper presents unique fine-scale imagery of compound and superimposed bedforms, which is used to (1) assess the physical forcing and evolution of a bedform field in San Francisco Bay, and (2) in conjunction with numerical modeling, gain a better fundamental understanding of boundary-layer flow dynamics that result in the observed superimposed bedform orientation.     

Three-dimensional SeaMARC II,gravity, and magnetics study of large-offset rift propagation at the Pito Rift,Easter microplate

We present results from a SeaMARC II bathymetry, gravity, and magnetics survey of the northern end of the large-offset propagating East Rift of the Easter microplate. The East Rift is offset by more than 300 km from the East Pacific Rise and its northern end has rifted into approximately 3 Ma lithosphere of the Nazca Plate forming a broad (70–100 km) zone of high (up to 4 km) relief referred to as the Pito Rift. This region appears to have undergone distributed and asymmetric extension that has been primarily accommodated tectonically, by block faulting and tilting, and to a lesser degree by seafloor spreading on a more recently developed magmatic accretionary axis. The larger fault blocks have dimensions of 10–15 km and have up to several km of throw between adjacent blocks suggesting that isostatic adjustments occur on the scale of the individual blocks. Three-dimensional terrain corrected Bouguer anomalies, a three-dimensional magnetic inversion, and SeaMARC II backscatter data locate the recently developed magmatic axis in an asymmetric position in the western part of the rift. The zone of magmatic accretion is characterized by an axis of negative Bouguer gravity anomalies, a band of positive magnetizations, and a high amplitude magnetization zone locating its tip approximately 10 km south of the Pito Deep, the deepest point in the rift area. Positive Bouguer gravity anomalies and negative magnetizations characterize the faulted area to the east of the spreading axis supporting the interpretation that this area consists primarily of pre-existing Nazca plate that has been block faulted and stretched, and that no substantial new accretion has occurred there. The wide zone of deformation in the Pito Rift area and the changing trend of the fault blocks from nearly N-S in the east to NW-SE in the west may be a result of the rapidly changing kinematics of the Easter microplate and/or may result from ridge-transform like shear stresses developed at the termination of the East Rift against the Nazca plate. The broad zone of deformation developed at the Pito Rift and its apparent continuation some distance south along the East Rift has important implications for microplate mechanics and kinematic reconstructions since it suggests that initial microplate boundaries may consist in part of broad zones of deformation characterized by the formation of lithospheric scale fault blocks, and that what appear to be pseudofaults may actually be the outer boundaries of tectonized zones enclosing significant amounts of stretched pre-existing lithosphere.     

Morphology of Shatsky Rise oceanic plateau from high resolution bathymetry

Newly collected, high resolution multi-beam sonar data are combined with previous bathymetry data to produce an improved bathymetric map of Shatsky Rise oceanic plateau. Bathymetry data show that two massifs within Shatsky Rise are immense central volcanoes with gentle flank slopes declining from a central summit. Tamu Massif is a slightly elongated, dome-like volcanic edifice; Ori Massif is square shaped and smaller in area. Several down-to-basin normal faults are observed on the western flank of the massifs but they do not parallel the magnetic lineations, indicating that these faults are probably not related to spreading ridge faulting. Moreover, the faults are observed only on one side of the massifs, which is contrary to expectations from a mechanism of differential subsidence around the massif center. Multi-beam data show many small secondary cones with different shapes and sizes that are widely-distributed on Shatsky Rise massifs, which imply small late-stage magma sources scattered across the surface of the volcanoes in the form of lava flows or explosive volcanism. Erosional channels occur on the flanks of Shatsky Rise volcanoes due to mass wasting and display evidence of down-slope sediment movement. These channels are likely formed by sediments spalling off the edges of summit sediment cap.     

A coupled-mode system with application to nonlinear water waves propagating in finite water depth and in variable bathymetry regions

A non-linear coupled-mode system of horizontal equations is presented, modelling the evolution of nonlinear water waves in finite depth over a general bottom topography. The vertical structure of the wave field is represented by means of a local-mode series expansion of the wave potential. This series contains the usual propagating and evanescent modes, plus two additional terms, the free-surface mode and the sloping-bottom mode, enabling to consistently treat the non-vertical end-conditions at the free-surface and the bottom boundaries. The present coupled-mode system fully accounts for the effects of non-linearity and dispersion, and the local-mode series exhibits fast convergence. Thus, a small number of modes (up to 5–6) are usually enough for precise numerical solution. In the present work, the coupled-mode system is applied to the numerical investigation of families of steady travelling wave solutions in constant depth, corresponding to a wide range of water depths, ranging from intermediate depth to shallow-water wave conditions, and its results are compared vs. Stokes and cnoidal wave theories, as well as with fully nonlinear Fourier methods. Furthermore, numerical results are presented for waves propagating over variable bathymetry regions and compared with nonlinear methods based on boundary integral formulation and experimental data, showing good agreement.     

Swath bathymetry: principles of operation and an analysis of errors

The principles of swath bathymetry are described, and the main cause of depth error is identified as acoustic interference, particularly from the sea surface. An error analysis is presented which gives the relationship among depth errors, the signal-to-interference ratio, the grazing angle, receiver spacing, and area resolution. It permits a prediction of when its measurement of depth can meet the accuracies required for nautical charting. Ways of reducing multipath interference and of minimizing its effect when it does occur are discussed. Particularly important are area averaging, the use of widely spaced receivers with ambiguities resolved by the vernier technique, and phase tracking for avoiding bias problems     

The East Pacific Rise and its flanks 8–18° N: History of segmentation,propagation and spreading direction based on SeaMARC II and Sea Beam studies

SeaMARC II and Sea Beam bathymetric data are combined to create a chart of the East Pacific Rise (EPR) from 8°N to 18°N reaching at least 1 Ma onto the rise flanks in most places. Based on these data as well as SeaMARC II side scan sonar mosaics we offer the following observations and conclusions. The EPR is segmented by ridge axis discontinuities such that the average segment lengths in the area are 360 km for first-order segments, 140 km for second-order segments, 52 km for third-order segments, and 13 km for fourth-order segments. All three first-order discontinuities are transform faults. Where the rise axis is a bathymetric high, second-order discontinuities are overlapping spreading centers (OSCs), usually with a distinctive 3:1 overlap to offset ratio. The off-axis discordant zones created by the OSCs are V-shaped in plan view indicating along axis migration at rates of 40–100 mm yr^–1. The discordant zones consist of discrete abandoned ridge tips and overlap basins within a broad wake of anomalously deep bathymetry and high crustal magnetization. The discordant zones indicate that OSCs have commenced at different times and have migrated in different directions. This rules out any linkage between OSCs and a hot spot reference frame. The spacing of abandoned ridges indicates a recurrence interval for ridge abandonment of 20,000–200,000 yrs for OSCs with an average interval of approximately 100,000 yrs. Where the rise axis is a bathymetric low, the only second-order discontinuity mapped is a right-stepping jog in the axial rift valley. The discordant zone consists of a V-shaped wake of elongated deeps and interlocking ridges, similar to the wakes of second-order discontinuities on slow-spreading ridges. At the second-order segment level, long segments tend to lengthen at the expense of neighboring shorter segments. This can be understood if segments can be approximated by cracks, because the propagation force at a crack tip is directly proportional to crack length.There has been a counter-clockwise change in the direction of spreading on the EPR between 8 and 18° N during the last 1 Ma. The cumulative change has been 3°–6°, producing opening across the Orozco and Siqueiros transform faults and closing across the Clipperton transform. The instantaneous present-day Cocos-Pacific pole is located at approximately 38.4° N, 109.5° W with an angular rotation rate of 2.10° m.y.^–1 This change in spreading direction explains the predominance of right-stepping discontinuities of orders 2–4 along the Siqueiros-Clipperton and Orozco-Rivera segments, but does not explain other aspects of segmentation which are thought to be linked to patterns of melt supply to the ridge axis.There are 23 significant seamount chains in the mapped area and most are created very near the spreading axis. Nearly all of the seamount chains have trends which fall between the absolute and relative plate motion vectors.     

Controls on crustal accretion along the back-arc East Scotia Ridge: constraints from bathymetry and gravity data

This study investigates crustal accretion processes along the East Scotia Ridge (ESR), an intermediate-rate back-arc spreading center with ten segments (E1–E10) that strike north–south. Mantle Bouguer anomaly (MBA) was calculated for the ESR region using satellite-derived and shipboard data sources. De-trended MBA (MBA_det) was determined by removing a residual plane from the MBA map, and ΔMBA_det was defined as the along-segment change in MBA_det. ΔMBA_det, as well as segment-averaged values of Na₈, Fe₈, and ⁸⁷Sr/⁸⁶Sr obtained from the published literature, generally appear to be better correlated with dist_sst (the distance from each segment center to the nearest point on the South Sandwich Trench) than with spreading rate. For each of the northern segments E2 through E6, MBA_det has a central low. MBA_det values also form a broad, longer-wavelength low from segments E2 through E6. Generally speaking, these findings are consistent with earlier studies such as Livermore et al. (Earth Planet Sci Lett 150:261–275, 1997) in suggesting that the region around segment E2 is a center for focused accretion along the ESR. On the other hand, southern segments E7 and E8 have central MBA_det highs, and MBA_det decreases somewhat linearly from segment E7 to E9, notwithstanding intrasegment variations. The quasi-linear MBA_det trend along these ESR segments is similar to that observed for the southernmost Lau spreading centers (e.g., Martinez and Taylor in Nature 416:417–420, 2002). Overall, plate boundary geometry and three-dimensional mantle flow may play a significant role in melting processes along the ESR, especially if the spreading center is processing geochemically heterogeneous South Atlantic mantle.     

Morphology and structure of the Camarinal Sill from high-resolution bathymetry: evidence of fault zones in the Gibraltar Strait

The Gibraltar Strait is the very narrow neck which connects the Atlantic Ocean and the Mediterranean Sea. The causes and mode of its opening at the end of the Messinian Salinity Crisis are still a matter of debate, and models based on eustatic rise and/or topographic lowering due to either erosion or faulting are generally evoked. We investigated the presence of faults based on a morphological and structural analysis of the Camarinal Sill, the shallowest passage in the Gibraltar Strait (<100 m water depth in places). This sill connects the Spanish and Moroccan shelves, and probably represents a structural high inherited from the Miocene compressive tectonics which took place in the external zones of the Betic-Rif orogenic arc. Our high-resolution bathymetric data enabled us to identify and interpret the origin of major morphological features in the area, including canyons, channels and a landslide, which we name the Tarifa landslide. Topographic arguments suggest that the Camarinal Sill is crossed by two main E-W- to ENE-WSW-directed fault zones which bound areas with different distribution, orientation and slopes of both scarps and crests. We name these the Hercules and Tarik fault zones, north and south of the sill respectively. The Hercules fault zone probably incorporates a normal movement component, whereas kinematic indicators are poor along the Tarik fault zone. The age of faulting is poorly constrained in both cases. Together with existing evidence of faults onland, the presence of these fault zones implies that they could be responsible for, or have contributed to, the opening of the Gibraltar Strait.     

Passive systems theory with narrow-band and linear constraints: Part II -- Temporal diversity

This paper is part of a series of three papers studying passive tracking problems arising in navigation and positioning applications. The basic question here lies with the determination of the position and dynamics of a point source being tracked by an omnidirectional observer, through demodulation of the Doppler effect induced on the radiated signals by the relative motions. A simple model, fitting a finite parameter nonlinear estimation context, is developed, the receiver designed, and its mean-square error performance studied. It is shown that, besides the speed and angle estimation, simultaneous global range passive tracking is possible. The signal model precludes range acquisition from synchronous measurement of the absolute phase reference: the global range estimation is attained by processing the higher order temporal modulations (varying Doppler). Quantifying the statistical and geometric performance tradeoffs, the work presents simple expressions and graphical displays that can be used as design tools in practical passive tracking problems. A subsequent paper considers the space/ time coupling issues, generalizing the study to the context where a moving source is tracked by a directional array.     

Tectonic shortening and gravitational spreading in the Gulf of Cadiz accretionary wedge: Observations from multi-beam bathymetry and seismic profiling

The Gulf of Cadiz lies astride the complex plate boundary between Africa and Eurasia west of the Betic-Rif mountain belt. We report on the results of recent bathymetric swathmapping and multi-channel seismic surveys carried out here. The seafloor is marked by contrasting morphological provinces, spanning the SW Iberian and NW Moroccan continental margins, abyssal plains and an elongate, arcuate, accretionary wedge. A wide variety of tectonic and gravitational processes appear to have shaped these structures. Active compressional deformation of the wedge is suggested by folding and thrusting of the frontal sedimentary layers as well as basal duplexing in deeper internal units. There is evidence for simultaneous gravitational spreading occurring upslope. The very shallow mean surface and basal slopes of the accretionary wedge (1° each) indicate a very weak decollement layer, geometrically similar to the Mediterranean Ridge accretionary complex. Locally steep slopes (up to 10°) indicate strongly focused, active deformation and potential gravitational instabilities. The unusual surface morphology of the upper accretionary wedge includes “raft-tectonics” type fissures and abundant sub-circular depressions. Dissolution and/or diapiric processes are proposed to be involved in the formation of these depressions.