
LAB thinning dominates; Moho relief partially cancels; topography modulates near axis.

Mature rift signature: broad uplift + axial graben + thinned crust (Moho 22 km) + mantle dome.

Active broad low (B&R) → mature rift with magmatic high (EAR) → failed rift with central gabbroic high (Keweenawan).

Western lows (active extension) — Midcontinent Rift high (failed) — Appalachian low (orogenic root). After USGS FS 78-95 (public domain).
import xarray as xr
import numpy as np
from scipy.signal import find_peaks
url = "https://www.ngdc.noaa.gov/geomag/EMAG2/EMAG2_V3_20170530.nc"
emag = xr.open_dataset(url)
profile = emag["z"].sel(lat=47.0, method="nearest").sel(lon=slice(-132, -126))
distances_km = (profile.lon - (-129.0)) * 111 * np.cos(np.radians(47))
peaks, _ = find_peaks(profile.values, prominence=80, distance=20)
# Match to GPTS (Cande & Kent 1995); fit slope through origin.
The fourth canonical data-access block of the course.

Find peaks → match to GPTS → linear fit → spreading rate. One line, in principle.
Six stages: stable craton → incipient → mature → break-up → slow MOR → fast MOR.
For each, predict sign + magnitude of: Bouguer · heat flow · topography · · magma · seismicity. 5 minutes.

Same six observables, six qualitatively different signatures. EAR spatially preserves stages 2–4.

Trend after Bell et al. 2022 (CC-BY 4.0): fast → shallow continuous AMC; slow → deep, episodic.
Cascadia = ridge offshore + trench onshore + arc through the city. Drive distance from Seattle.
scipy.signal.find_peaks succeeds in the data-rich middle: clean profiles, strong signal, smoothly varying noise.You must know enough geophysics to recognise when the tool will lie.
Open question: Why do some continental rifts break to spreading and others fail?
See you next time for L28: convergent margins.