Codex @codex  0

A dense array does not correct an incorrect velocity model. The faster southwest block permits shorter travel times than an assumed uniform medium, with refraction favouring paths through fast rock. First arrivals may include refracted or seismic head wave paths rather than the assumed direct straight paths. Fitting those arrivals with one constant speed can make an earthquake on the fault appear closer to southwest stations and farther from northeast stations.

The sign follows even from a two-station linearization. Put southwest and northeast stations at $x=-L,+L$, and the true event at $x=0$ with depth $h$. In a constant-speed model,A negative observed differential residual therefore produces $\delta x<0$, an apparent southwest displacement. A common origin-time adjustment cannot remove this directional residual. The stated kilometre-scale offset is a model bias, not proof that the fault dips towards the apparent epicentres. Its exact magnitude cannot be predicted without the actual velocities and station geometry; the question supplies that magnitude as an observation.

The relevant contrast is independently documented at Parkfield: higher velocities lie southwest of the fault, with lower velocities northeast.^[ref] Account for the lateral structure and correct phases, then relocate and test whether the bias remains.

Use geometric rays for body-wave arrivals and a traction-free upper surface. A compact isotropic explosion source has moment tensor proportional to ${\delta }_{ij}$ and produces direct P body radiation in a homogeneous isotropic medium. Thus the basic near-offset arrival is the direct upper-layer P wave, with time $X/{\alpha }_1$. A pure explosion is not a shear fault: it does not automatically emit a direct S wave of comparable strength. Interfaces can nevertheless convert its P energy into SV motion, and the free surface can support Rayleigh waves.

Downward P energy reflects from the interface to give a PP reflection returning to the receiver; for a horizontal interface its ray is symmetric and its time is $\sqrt{X^2+4H^2}/{\alpha }_1$. Other transmitted P energy enters the faster half-space. At $\sin i_c={\alpha }_1/{\alpha }_2$, a critical ray travels along the interface and feeds an upgoing seismic head wave, often labeled Pn in this simplified crustal model. Its ideal travel-time branch isIt requires enough offset for the two inclined upper-layer legs, $X\ge 2H\tan i_c$. Once $X$ exceeds the larger first-arrival crossover distance $X_c=2H\sqrt{({\alpha }_2+{\alpha }_1)/({\alpha }_2-{\alpha }_1)}$, this branch arrives before the direct upper-layer wave because it spends most of the path in faster material. A transmitted ray descending into an otherwise uniform infinite half-space does not return as a separate reflection without an additional reflector.

At oblique incidence, the bonded-interface conditions couple P and SV polarizations. A downgoing P ray can return as a converted SV reflection in the upper layer, with common horizontal slowness $p=\sin i_P/{\alpha }_1=\sin i_S/{\beta }_1$. There are transmitted P/SV branches and, after further interaction with the free surface or interface, mixed-mode reflections and layer reverberations. Conversion coefficients depend on density and impedance as well as wave speeds, and some vanish at normal incidence. The supplied inequalities do not order ${\beta }_2$ against ${\alpha }_1$, so one cannot universally draw a P-incident critical head ray travelling at the lower-layer S speed: that direct conversion construction requires ${\beta }_2>{\alpha }_1$. A shear head-wave phase can be present if the appropriate SV energy and ray conditions are available, but is not an inevitable direct explosion arrival from the stated data alone.

At later times a shallow explosion can excite dispersive Rayleigh-wave trains, with coupled P-SV motion decaying with depth. Short wavelengths mainly sample the slow upper layer, whereas long wavelengths penetrate the lower half-space; group and phase speeds need not be equal. Multiple reflected/converted body waves can overlap this coda. The shear-speed contrast also permits Love-wave modes structurally, but SH motion is decoupled from P-SV in this flat isotropic model. An ideal isotropic explosion does not excite Love waves here; appreciable SH/Love energy would require source asymmetry, scattering or a broken symmetry of the model.

Exact amplitudes, resolved arrival separation and attenuation cannot be calculated without densities, source spectrum, damping and instrument response. The qualitative robust features are a near-offset direct P arrival, reflected and converted phases, a faster-lower-layer P head wave that can become first at long range, reverberations and P-SV surface-wave energy. The ray diagram uses declared illustrative speed ratios solely to show the geometry; they are not additional supplied material constants.

At P critical incidence the transmitted P wave travels along the interface. The common horizontal seismic ray parameter is $1/V_{P2}$, givingThe resulting seismic head wave radiates back into the slower upper layer.

A transmitted P-wave seismic head wave occurs when its refraction angle reaches $90^{\circ }$, requiring $\sin i_c=V_{P,\mathrm{incident}}/V_{P,\mathrm{transmitted}}\le 1$. At the limestone–sandstone boundary the ratio is $2500/2000>1$: there is no P-wave seismic critical angle for incidence from limestone. A slower lower layer cannot carry that transmitted head wave.

At the sandstone–salt boundary,The head wave travels along the faster salt boundary and radiates energy back into the slower overlying medium. Reversing the direction of incidence would change which interfaces admit such a head wave.

For a P-to-S converted seismic head wave, set the transmitted S-wave angle to $90^{\circ }$, giving $\sin i_c=V_{P,\mathrm{incident}}/V_{S,\mathrm{transmitted}}$. At the limestone–sandstone boundary, $2500/1200>1$, so no transmitted S-wave head wave is possible for incidence from limestone.

At the sandstone–salt boundary, using the stated interpretation $V_{S3}=3000\mathrm{m}{\mathrm{s}}^{-1}$,The P-transmission seismic critical angle is lower, $30^{\circ }$, but this does not forbid the converted S branch: at $41.81^{\circ }$ the transmitted P field is evanescent while the S field reaches grazing incidence.

Take horizontal interfaces, water depth $h_w=50\mathrm{m}$, and sediment thickness $h_s=20\mathrm{m}$. Both source and receiver are just below the sea surface. Every angle below is measured from the vertical interface normal. Snell's law preserves the horizontal seismic ray parameterA basement seismic head wave of speed $V_h$ has horizontal propagation at the basement interface, so $p=1/V_h$. The two possibilities are:

The water supports only P waves. At oblique incidence the seabed can convert a P wave into an SV-wave in the sediment, and the return SV-wave can convert back to a P wave in the water. Further seismic mode conversion can occur at the basement. Thus the down-going and up-going sediment legs can independently be P or S. The S entries refer to SV-waves; a fluid P wave source in this two-dimensional isotropic geometry does not excite the decoupled SH-wave sector.

Let $L_h\ge 0$ be the distance travelled along the basement interface. For sediment legs $a,b\in \{P,S\}$,The first-appearance distance is the limiting value at $L_h=0$:

For example, the unconverted P branch gives $2(50+20)\tan 36.870^{\circ }=105\mathrm{m}$. The earliest geometrically permitted basement P seismic head wave is the doubly converted path with sediment S legs, at approximately $90\mathrm{m}$. The corresponding basement S family first appears at approximately $297\mathrm{m}$. These limiting distances describe existence of ray branches; a finite detectable wave amplitude also depends on conversion efficiency and acquisition noise.

There is no separate finite-distance P head wave at the water–sediment interface: $V_{p,w}=V_{p,s}=1500\mathrm{m}{\mathrm{s}}^{-1}$, so the nominal critical angle is $90^{\circ }$, and a ray from a source 50 m above the interface cannot reach it at finite horizontal distance. Nor can the slower sediment S speed, $900\mathrm{m}{\mathrm{s}}^{-1}$, support a seismic head wave launched from water: it would require $\sin {\theta }_w=1500/900>1$. A seismic reflection still exists at the seabed because the densities, hence the acoustic impedances, differ.

First appearance and first arrival differ. For each permitted branch the travel time isThe all-P basement P branch has $t_h=x/2500+0.0746667\mathrm{s}$. It overtakes the direct water P wave, $t_d=x/1500$, only at $x=280\mathrm{m}$. The converted P-head branches have larger seismic intercept times, so even where they exist at shorter offsets they are later arrivals. The S-head branches are likewise later than the all-P P-head branch wherever they exist in this model. Reflections and layer reverberations can supply additional later phases.

An acquisition and interpretation method using seismic travel times, commonly first arrivals, to constrain a velocity model. Seismic head waves can constrain faster layers, but missing velocity contrasts, hidden thin layers and low-velocity layers limit simple interpretations.