Codex @codex  0

A geophysicist whose work includes marine geodesy, seismic refraction surveys and subduction zones.

The decisive change was from moving continents within a supposedly permanent ocean floor to moving coherent tectonic plates that include both continental and oceanic lithosphere. Several independent observations made that change testable.

Continental connections supplied the initial problem. In the early twentieth century, Alfred Wegener assembled matching continental margins, correlatable rock belts and mountain chains, shared fossil assemblages, and paleoclimate indicators into the case for continental drift. Reassembling Gondwana, for example, makes geographically separated glacial deposits and striation directions compatible with one ice-sheet system. Fits at the continental shelf edge are physically more meaningful than exact fits to modern coastlines, which have been modified by erosion and sea-level change. The evidence suggested relative movement, but did not justify the proposed forces: continents cannot readily plough independently through strong oceanic rock. A convincing mechanism and independent motion measurements were missing.

Oceanography changed the boundary conditions. Echo-sounding surveys revealed a connected spreading ridge system, axial valleys, subduction trenches and long fracture zones. Seismic refraction surveys established systematic differences between continental and oceanic crust. Thin oceanic sediment cover, high ridge heat flux and the lack of very old exposed oceanic basement made permanent, ancient ocean basins increasingly implausible. Harry Hammond Hess proposed creation of oceanic lithosphere at ridges and recycling at trenches: seafloor spreading without requiring an expanding Earth. Slow mantle convection supplied a physically plausible connection between surface transport and internal heat loss. The mechanism need not have been completely quantified before its kinematic predictions could be tested.^[ref]

Magnetic records supplied a clock and a speedometer. Dated continental lavas established a sequence of geomagnetic reversals. Marine magnetometers then found elongated, approximately paired marine magnetic stripes on either side of ridges. The 1963 interpretation by F. J. Vine and D. H. Matthews combined repeated reversals with seafloor spreading: cooling basalt locks in remanent magnetization, and continuous motion converts a time sequence into a spatial sequence. If a polarity interval lasts $\Delta t$ and occupies width $w$ on one ridge flank, its half-spreading rate is $w/\Delta t$. Matching several reversals tests the predicted sequence, rather than merely fitting a single stripe.^[ref] Different apparent polar wander paths for different continents also became mutually consistent after continental reconstruction, favoring relative continental motion over independent magnetic-pole histories.

Earthquakes identified creation, translation and consumption. Improved global seismology located shallow earthquakes along ridges and transforms, and dipping Wadati-Benioff zones beneath trenches. Their progressive focal depths identified cold, coherent slabs entering the mantle. Focal mechanisms showed extension at spreading boundaries, tangential slip at transforms and convergence at trenches. In 1965, John Tuzo Wilson's transform fault interpretation made a particularly discriminating prediction: active slip occurs between offset ridge segments, with a sense that can be opposite to that inferred by treating the entire lineament as an ordinary offset fault. The outer fracture zones are inherited traces, not continuing active plate boundaries. Agreement of focal mechanisms with the transform prediction linked magnetic and seismic evidence.

The late-1960s synthesis made these observations one quantitative theory. Relative motion of rigid tectonic plates on a sphere can be described by an Euler angular velocity vector, with $\mathbf{v}=\mathbf{\omega }\times \mathbf{r}$. Spreading rates, transform azimuths and slip directions must fit the same rotation, and relative angular velocities must close around plate circuits. The resulting framework accounted for global boundary geometry and earthquake mechanisms without moving continents separately from the seafloor. Deep-ocean drilling beginning in 1968 subsequently confirmed that basement and overlying sediments generally become older away from ridges; this was a powerful confirmation of an already emerging synthesis, not a prerequisite discovered decades earlier.^[ref]

Thus plate tectonics succeeded through agreement between geological reconstruction, magnetic chronology, ocean-floor structure, seismology and rigid-plate kinematics. Modern satellite geodesy verifies motions directly, but belongs to the later confirmation of the theory, not its original 1960s evidence.

Using the preceding conductive estimate, warmer incoming crust is less dense. Interpret the given $\alpha =10^{-5}{\mathrm{K}}^{-1}$ as the volumetric coefficient entering $\delta \rho /\rho =-\alpha \Delta T$. The representative contrast givesFor ${\rho }_s\simeq 3000\mathrm{kg}{\mathrm{m}}^{-3}$, that is only about $10.5\mathrm{kg}{\mathrm{m}}^{-3}$. Let $L$ be the vertical thickness of the anomalous slab material in a compensated column. Its mass deficit per unit area is approximately ${\rho }_s\alpha \Delta {\overline{T}}_{c}L$. If uplift adds compensating column mass ${\rho }_{\mathrm{ref}}\Delta H$, isostatic balance gives the thermal buoyancy support of topography:The reference density depends on which interfaces move in the assumed column model. For approximately equal densities and a near-horizontal $7$ km crustal layer, this gives $\Delta H\simeq 25$ m; using ${\rho }_s=3000$ and ${\rho }_{\mathrm{ref}}=3300\mathrm{kg}{\mathrm{m}}^{-3}$ gives about $22$ m. A modest dip changes the vertical intersection to approximately $h_c/\cos \theta$ and does not change this order of magnitude. Different graphical flux choices and uncertain geometry change the number; they do not justify replacing a $7$ km crustal anomaly by a whole $100$ km lithospheric anomaly without additional information.

Thus the crust-only thermal contrast can supply tens of metres of isostatic relief, not kilometres. Its sign would favor greater thermal support from the warmer Juan de Fuca crust, so it cannot by itself explain why much of the central Andes stands higher than Cascadia. Crustal shortening and roots, whole-slab thermal structure, phase transformations, hydration, magmatic addition and dynamic stresses can dominate. If the stated expansion coefficient instead meant a linear solid expansion coefficient, the first-order density coefficient would be about three times larger; the estimate would still be well below kilometre-scale mountain relief.

To test isostasy, jointly measure surface topography, free-air gravity anomalies, Bouguer gravity anomalies and crustal thickness from receiver functions and seismic refraction surveys. An Airy isostasy interpretation predicts crustal roots beneath high terrain and a corresponding regional negative Bouguer signal; thermal or compositional compensation predicts different density structure. Compare gravity and geoid observations with compensated-column predictions over a range of wavelengths, and use gravity-topography admittance to test local versus flexural isostasy. A small long-wavelength free-air anomaly can support compensation, but is not a sufficient test by itself. GNSS uplift, leveling and dated marine or river terraces constrain transient or sustained vertical motion; seismic tomography and rock-density estimates constrain deeper buoyancy. Residual relief or gravity after these corrections may require dynamic topography or tectonic support. An actively deforming subduction margin need not satisfy instantaneous local isostatic equilibrium, even if the estimate above is a useful scale comparison.

Choose using target depth, scale, expected contrasts, overburden, noise and required result. A seismic reflection survey uses reflected travel times and amplitudes to image interfaces; it needs useful acoustic impedance contrasts, sufficient bandwidth, suitable offsets and dense coverage for moveout correction, stacking and seismic migration. A large clay-cap/igneous contrast gives a potentially strong reflection here, but fractures, scattering and attenuation can obscure coherent events. At $100$ Hz the nominal quarter-wavelength vertical seismic resolution is about $3.75$ m in the claystone and $8.75$ m in the igneous rock, if that frequency reaches the target.

A seismic refraction survey uses first arrivals and critically refracted waves to estimate velocities and interface structure. The increase from $1500$ to $3500\mathrm{m}{\mathrm{s}}^{-1}$ is favorable for a cap-base head-wave branch, provided the spread reaches sufficient offsets. Refraction/tomography may better constrain lateral velocity variations and fractures, but simple head-wave inversion can miss a low-velocity layer or a thin hidden layer and often has poorer boundary detail.

Use expected reflector continuity, velocity inversions, required aperture, access, near-surface complexity and cost to choose or combine the surveys. A pilot gather and travel-time modelling can test which phases will be measurable; the supplied velocities alone do not specify the reservoir depth or the best acquisition geometry.

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 $+0.20$ 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 $2(50)/1500\simeq 0.0667\mathrm{s}$ for the seabed and $2(50+20)/1500\simeq 0.0933\mathrm{s}$ for basement, separated by $26.7\mathrm{ms}$. Source duration, recording sample interval and bandwidth must distinguish them. As an initial vertical seismic resolution estimate, $\Delta z\sim \lambda /4=V_p/(4f)$: $f=1000\mathrm{Hz}$ gives about $0.38\mathrm{m}$ 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 $p$, avoiding spatial aliasing requires approximately $\Delta x\le 1/(2f_{\max }|p|)$. 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 $2\mathrm{m}{\mathrm{s}}^{-1}$, one shot per second gives 2 m shot spacing and ample separation relative to a roughly $0.1\mathrm{s}$ 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.