Chirp sub-bottom profiler 2026-10-25
A shallow marine acoustic profiling system transmitting a swept-frequency pulse and using matched processing to resolve interfaces below the seabed. Fine vertical seismic resolution must be balanced against penetration and seismic attenuation.
The asthenosphere exists because the mantle's ability to deform is strongly dependent on temperature, pressure, composition, water and deformation mechanism. Cooling from the surface creates a strong lithosphere. Beneath it, hotter mantle deforms by material creep at much lower stress on geological timescales. The asthenosphere is predominantly solid weak rock, not a global ocean of magma. On the much shorter timescale of a seismic wave, it still supports S waves.
The hot shallow mantle can be close to its solidus. Water weakening of rocks and, in some settings, a small connected melt fraction can further lower effective viscosity. At greater depths, changes in pressure, homologous temperature, mineral structure and water storage modify the deformation laws; increasing depth does not imply indefinitely decreasing viscosity. Explaining a weak upper-mantle interval does not require every location to contain the same amount of melt.
There are several physically distinct boundaries to test. The thermal lithosphere-asthenosphere boundary describes a thermal transition, the mechanical boundary a contrast in long-term strength or viscosity, and the seismic lithosphere-asthenosphere boundary a change in seismic velocity or other wave properties. They need not have identical depths. A smooth geothermal gradient can explain broad thermal weakening, but a very sharp decrease of seismic velocity may require changes in melt fraction, hydration, composition or seismic anisotropy. A sharp shallow reflector within a craton might instead be an internal lithospheric boundary.[ref]
Begin with a thermal null model. A half-space cooling model gives
If a thermal boundary corresponds to a fixed temperature , its depth is
so it deepens as the square root of plate age. Test that prediction against independently determined seafloor ages, bathymetry and heat flux. An oceanic plate cooling model adds finite thickness and basal heat supply, allowing the old-age depth to approach a finite value. For continents, include radiogenic heating, crustal thickness, the history of continental rifting or tectonic thickening, and advective heating by magma or mantle flow. A cold, dry craton and a recently extended volcanic region should not share the same thermal profile.
Then calculate the observable consequences of competing mechanisms rather than merely assigning a boundary to a map:
Compare models against the observed boundary depth and its sharpness, including uncertainty and resolution. An age–depth trend with matching heat flux favours thermal control; a sharp discontinuity that fails the thermal prediction but coincides with an electrical anomaly can support a melt- or hydration-related transition. Neither observation alone proves a unique cause. The spatially varying boundary is best understood as the result of thermal evolution interacting with composition, fluids and deformation history.
Probably yes in connected fractures or transiently opened pore space, especially where a transform fault admits circulating seawater; a uniform water-filled porous layer several kilometres thick is not guaranteed. Porosity and permeability of a porous medium generally decrease as cracks close or seal, while ongoing faulting can reopen pathways. Fluid can occur at high pressure and temperature, including a supercritical state; it need not resemble cold surface water.
Several complementary tests address presence and origin:
No one anomaly uniquely identifies water. Seismic velocity can also reflect cracks, composition and temperature; electrical conductivity can reflect melt or conductive minerals. Serpentinization demonstrates incorporation of water into minerals, which must be distinguished from fluid currently occupying pore space. At an oceanic transform fault, an external seawater supply is the natural first hypothesis, but isotope ratios and chemical budgets can test contributions from deep fluids rather than presuming the source. A research example of coupled deformation and hydrothermal alteration is Mass transfer and fluid flow during detachment faulting and development of an oceanic core complex, Atlantis Massif by Chiara Boschi and others; a detachment system illustrates the methods without proving identical flow at every transform fault.
At a convergent margin, the incoming oceanic lithosphere is a chemically and mechanically layered package: sediments, oceanic crust and lithospheric mantle. Subduction does not deliver this whole package unchanged to one destination. Solids can be accreted, underplated, exhumed or recycled deeply, while released fluids and melts carry selected components to the overriding plate. The outcome depends on incoming composition and hydration, slab age, convergence rate, thermal structure, interface coupling and subsequent tectonic history.
Figure 1. . Original process diagram. Solid-material pathways are separated from fluid and melt transfer. The arrows represent possible routes, not universal proportions or fixed reaction depths; temporary storage and later remobilization are possible. Details.
The incoming material already contains a history. Sediments include terrigenous debris, biogenic carbonates, siliceous deposits, clay minerals, pore water and organic matter. Oceanic crust consists mainly of basalt, sheeted dikes and gabbro; hydrothermal alteration can add water, carbonates and mobile elements before subduction. Beneath it, residual peridotite forms the lithospheric mantle. Faulting at a spreading environment, a transform boundary or an outer-rise bending region can let seawater enter and cause serpentinization. The lithospheric mantle is not uniformly hydrated to its full thickness. Consequently, sediment thickness alone does not measure the incoming volatile inventory.
Sediments: separation starts near the trench. Compaction and increasing effective pressure expel pore fluids, which can escape through the forearc, transport dissolved constituents and influence pore pressure on the plate interface. Some sediments are scraped off into an accretionary prism, shortened by thrusts and retained in the upper plate. Others pass beneath the leading edge, then detach and undergo tectonic underplating beneath the forearc. Deeper underplated material may return to the surface during tectonic exhumation, preserving a record of burial and deformation. At a sediment-starved erosive margin, little prism may form, and tectonic erosion can entrain material from the overriding plate as well. Thus net continental gain and net continental loss are both possible; they can even alternate at one margin.
Sediments that remain attached to the slab undergo metamorphism, dehydration, and potentially decarbonation or carbonate dissolution. The released fluids can transport boron, barium, strontium, lead and other soluble constituents. Carbon can leave as dissolved species, remain in carbonate, or occur as reduced carbon depending on redox conditions and reaction history. Along sufficiently hot paths, sediment-derived melts or mixtures of sediment and altered rock can transfer material into the mantle wedge; buoyant mixed material may rise as a mélange diapir. Cold paths can instead retain more solid material and some volatile-bearing phases to greater depths. These are alternative reaction and transport pathways, not a single compulsory sequence. A carbon-budget study explicitly includes fluid dissolution, sediment diapirism and storage in the upper plate, showing why volcanic carbon dioxide alone cannot close the whole subduction carbon budget.[ref]
Oceanic crust: metamorphic transformation and chemical filtration. Altered basalt and gabbro are progressively buried. Depending on their pressure-temperature paths and reaction kinetics, they can pass through blueschist facies or warmer assemblages and eventually form eclogite, dominated by garnet and omphacite. Eclogitization increases density and can strengthen negative slab buoyancy, but a kinetic delay can leave metastable rock. It is incorrect to assign every slab the same sequence or the same fixed depths for these changes.
During burial, breakdown of hydrous phases releases water-rich fluids; different mineral assemblages retain different quantities of water. The fluids migrate through or along the interface and into the mantle wedge, where they cause metasomatism and lower the solidus of peridotite. Flux melting of the mantle wedge generates much of the primary magma feeding a volcanic arc. This is different from saying that the entire descending oceanic crust melts beneath the arc. Most of its refractory mass can remain solid. Partial melting of oceanic crust is more plausible for unusually hot or young slabs, or other favorable thermal conditions, and produces a distinct crustal contribution.
The depth of dehydration is controlled by thermal structure and composition. Global thermal-petrological calculations, for example, predict more efficient deep water transport in old, fast, cold slabs than in hot ones, with substantial differences between volcanic crust, gabbro and hydrated peridotite. Such a calculation supports conditional deep retention; its modeled proportions should not be imposed on every margin.[ref]
Mechanical detachment adds more branches. Upper-crustal slices, seamounts or oceanic plateaus may be accreted or underplated instead of continuing down the slab. A complete or partial ophiolite can be emplaced onto another plate by obduction in suitable convergence histories. Most intact oceanic crust is not saved this way. High-pressure fragments may be entrained in an interface channel or mélange and exhumed; their small preserved volume should not be mistaken for the dominant global mass flux. Magma generated after slab-derived transfer may also pond and crystallize beneath the upper crust by magmatic underplating, which is a different process from mechanical tectonic underplating of solid slices.
Lithospheric mantle: residual rock and a variable water carrier. Most incoming lithospheric mantle is depleted peridotite left after extraction of ridge melts. Its cold, dense, mechanically coherent portion contributes to slab pull and generally sinks with the plate. Hydrated portions contain serpentine and possibly other hydrous minerals. Breakdown of these phases releases water at conditions that depend on the thermal path, feeding further metasomatism of the mantle wedge and sometimes dehydration embrittlement. The latter is a proposed contributor to intermediate-depth earthquakes, not an exclusive explanation for all such events. Very cold paths may transfer some water onward in high-pressure hydrous phases or as hydrogen incorporated in nominally anhydrous minerals.
Do not confuse the incoming hydrated lithospheric mantle with newly serpentinized forearc mantle in the overriding plate. The latter receives slab-derived fluids and can temporarily store their water. Hydrated, relatively buoyant interface material can facilitate return flow and tectonic exhumation, while the colder deeper bulk of the lithospheric mantle continues down. Some peridotite may also be accreted or obducted with an ophiolite. These exchanges mean that even an exposed serpentinite requires structural and chemical context before it can be assigned to its original plate.
Deep recycling and eventual return. The residual eclogite, remaining sediments and lithospheric mantle can enter the deeper mantle. A slab can temporarily stagnate around the mantle transition zone or penetrate into the lower mantle, depending on buoyancy, phase transformations, geometry and mantle viscosity. Crustal material can remain as compositionally distinct domains, mix, or react with surrounding peridotite; it need not retain the near-surface eclogite assemblage indefinitely as pressure rises. Subsequent mantle convection can return recycled components toward spreading regions or mantle plumes, where partial melting incorporates them into new basalt. A primary study of Hawaiian basalt illustrates a specific model in which melts derived from recycled crust react with peridotite to form a pyroxenite source. It is an interpretation of combined chemical evidence, not a unique label attached to every enriched basalt.[ref]
How the pathways can be observed. Seismic reflection profiles, drilling, wedge geometry and balanced structural sections constrain accretion, fluid escape and tectonic erosion. High-pressure metamorphism, mineral zoning, radiometric dating and pressure-temperature-time paths in exhumed rocks reveal burial, reaction and return. Their reconstructable histories are especially useful when deformation and fluid inclusions are linked to the dated mineral growth. Earthquake distributions, seismic tomography, seismic attenuation and electrical conductivity constrain the geometry and physical state of the slab and mantle wedge, although fluids, melts and temperature can produce nonunique signals.
Arc magmas carry complementary chemical evidence: relatively mobile barium, boron, lead or strontium can trace fluid-mediated transfer, whereas sediment-sensitive thorium and isotope ratios can help identify other contributions. Element mobility changes with fluid composition and pressure-temperature conditions, so no one element is an infallible tracer. Boron isotope ratios, radiogenic strontium isotope ratios and lead isotope ratios should be interpreted with measured incoming sediments and altered oceanic crust. Arc gases and forearc fluid discharge constrain volatile return; ocean-island basalts, mantle xenoliths and diamond inclusions provide possible evidence of deeper recycling. Each observation samples only part of the system.
For any element , a useful control-volume statement is
Here permanent accretion is counted outside the chosen reactive control volume, return includes fluid, melt and exhumed-solid exports, and transient underplating or reservoir growth contributes to storage until exported. For a different boundary choice, the terms must be regrouped consistently. The three incoming layers can have very different transfer efficiencies for bulk rock, water, carbon and a particular trace element. A small mass of escaping fluid can dominate an element's return while most of the solid slab continues into the deep mantle.
The star is near the southwestern end of the main fault slip patch, so a much longer part of the finite-source rupture propagates northeastward. For comparable distances, ground conditions and buildings, the northeastern site receives the stronger, shorter direct pulse; the southwestern site receives the weaker, longer direct shaking. This is rupture directivity.[ref]
For a segment of length , rupture velocity and seismic wave speed , successive arrival times span approximately
Forward arrivals overlap and concentrate energy into a shorter pulse. Backward arrivals spread out. The supplied finite-source rupture is not perfectly unilateral, so the prediction concerns the dominant contribution. Seismic radiation patterns, seismic attenuation, sediment amplification and building response can alter the ranking of actual damage or total coda duration.
Let a positive incident stress wavelet be and measure return times relative to a receiver just above the reservoir. The top interface gives a negative P return of amplitude and no basal P return. For a separately generated and measured S wave, the top return is also negative, , while the basal return is positive, at the lower interface.
If is reservoir thickness, the basal return arrives or after the top return; no numerical is supplied. Including downward and upward transmission through the top boundary, but ignoring seismic attenuation, spreading and seismic multiples, its observed S-wave amplitude is
Thus the synthetic S trace has an earlier negative wavelet and a later, smaller positive wavelet; the synthetic P trace has only the earlier negative wavelet. The figure chooses solely to separate the pulses visually. Receiver electronics or a displacement convention can reverse the displayed trace polarity, so the convention must accompany the sketch.
A marine hydrophone streamer with a purely normal-incidence P source does not directly record the hypothetical S trace. A seabed multicomponent receiver and suitable S excitation or oblique-incidence seismic mode conversion are required. This distinguishes the calculated solid-layer response from what a particular acquisition system can measure.
Figure 1.
Idealized normal-incidence P and S reflection wavelets, including transmission through the reservoir top
. Details.
Quality is degraded by vessel and flow noise, poor receiver coupling or position uncertainty, insufficient acquisition aperture and fold, coarse sampling and aliasing, and inaccurate source timing. Finite source bandwidth limits vertical seismic resolution; shallow-water reverberations, bubble pulses, source/receiver ghosts and seismic multiples overlap primary reflections. In of water with sound speed about , the water-layer round-trip time is only about .
Seismic attenuation removes high-frequency energy, geometric spreading weakens distant returns, and heterogeneous lithology, scattering and seismic anisotropy complicate travel paths and amplitudes. Incorrect seismic velocity models misplace reflectors during seismic migration; dipping structure and out-of-plane arrivals particularly challenge two-dimensional processing. Thin layers can cause wavelet interference and tuning. Finally, even noiseless data cannot recover a normal-incidence P reflection from the stipulated zero-contrast sandstone–salt boundary: that is absent signal rather than poor recording.
Use a high-resolution multichannel seismic reflection survey, with closely spaced crossing lines or a suitable small three-dimensional survey over the turbine locations. It should map the seabed, the sediment–sandstone boundary, local faults, buried channels, and lateral changes relevant to foundations. Seismic reflection responds to acoustic impedance contrasts and can image an irregular boundary without requiring it to be a first-arrival refractor. Here even the seabed has a useful normal-incidence reflection coefficient of despite its zero P-velocity contrast.
A seismic refraction survey can provide complementary seismic velocities and geotechnical information, particularly using seabed multicomponent receivers for S waves. Its limitations as the principal imaging method here include the missing water–sediment P refractor, required offsets of hundreds of metres to obtain useful basement first arrivals, and weaker constraints on thin layers or velocity inversions. Variable bathymetry and sediment thickness further complicate simple refraction interpretations.
A shallow-towed boomer seismic source or small sparker seismic source is a sensible initial seismic source for penetrating about 20 m of sediment beneath a 50 m water column. A chirp sub-bottom profiler can add finer imaging of the upper cover where penetration permits; penetration depends on grain size, consolidation, gas, and seismic attenuation. A pilot line should establish the minimum useful source energy and bandwidth. Large low-frequency air-gun arrays are unnecessary for this shallow target. Small coastal vessels, sea state, tow stability and acoustic background constrain usable frequencies and operations. Marine mammals and other local receptors motivate a low-energy effective source, monitoring and appropriate operational mitigation.[ref]
The expected normal-incidence two-way travel times are for the seabed and for basement, separated by . Source duration, recording sample interval and bandwidth must distinguish them. As an initial vertical seismic resolution estimate, : gives about in the sediment, although usable resolution depends on the whole seismic wavelet and signal-to-noise ratio.
Choose shot and receiver spacing from the shortest useful wavelength, target size, maximum reflector dip, seismic migration needs and achievable navigation accuracy. For a plane event of horizontal seismic ray parameter , avoiding spatial aliasing requires approximately . Source spacing and receiver spacing also set midpoint coverage and seismic acquisition fold; adequate offsets enable moveout correction and velocity estimation. For nearly horizontal shallow reflectors a few-metre shot interval can be practical, whereas steep events at high frequency require denser spatial sampling. These values must be designed together rather than imposing one spacing for every target. Cross-line spacing must resolve lateral features between turbine sites, not merely along each track.
At a vessel speed of , one shot per second gives 2 m shot spacing and ample separation relative to a roughly target return, subject to the longer record needed for seismic multiples and source recovery. Shallow towing limits seismic source ghosts and receiver ghosts, but requires stable tow depth. Accurate positioning, tides, measured water sound speed, bathymetry and checks for seabed seismic multiples are essential for correct target depths.
The shorter, stronger northern pulses and broader southern pulses suggest propagation from south to north, towards the northern station. The useful evidence is the difference in pulse widths, not just the earlier northern arrivals: the latter chiefly reflect station distance.
As rupture reaches successively more northerly patches, each patch radiates a new P wave and S wave pulse. Towards the north, the later emission time is partly cancelled by a shorter propagation path; those contributions arrive closer together. Towards the south, later emission and a longer path add, stretching the signal. This is rupture directivity. It refers to the advance of the slipping region along the fault, not necessarily the direction of the local fault slip vector. A thrust can slip up-dip while its rupture front advances in a direction with a northward component.
More generally, emission from a patch at distance along the rupture arrives at
The forward direction has a smaller arrival-time spread. This is the same travel-time compression captured by isochrone-based rupture directivity calculations.[ref] Larger northern amplitude supports the interpretation but is not decisive on its own because geometrical spreading, seismic attenuation, seismic radiation patterns and site response also affect amplitude. With only two traces, propagation, filtering and site effects remain alternatives to test; the idealized exercise is most naturally explained by northward directivity.