DELTA_T=15 # years gmt max net GPS_cascadia.xyz -T$DELTA_T -R-125/-122/42/49 -I0.2deg > cascadia_max.txt
gmt blockmean net_mag.xyz -Rregion -I0.5deg -C > blocked.xyz gmt gmtinfo blocked.xyz -C -Ine | Magnitude Range | Interpretation | |----------------|----------------| | < 5 mm/yr | Stable craton, far-field | | 5–20 mm/yr | Diffuse plate boundary | | 20–80 mm/yr | Active fault creep or interseismic loading | | > 80 mm/yr | Coseismic displacement, glacial isostatic adjustment (extreme) | gmt max net
Maximum Net Displacement: 312.4 ± 3.2 mm Location: -124.31°E, 44.23°N (near Cape Blanco, OR) Direction: N68°E (subduction convergence direction) DELTA_T=15 # years gmt max net GPS_cascadia
gmt max net gps_data.xyz -T5 -R-120/-110/30/40 -I0.1d -Egrid.nc -Sgaussian -V gmt max net —whether as a scripted workflow or a future built-in module—provides a robust, reference-frame-aware method to locate extreme net displacements in geodetic data. By combining GMT’s vector math, filtering, and statistical tools, it delivers actionable insights for earthquake geology, crustal dynamics, and hazard assessment. Note to users: The exact implementation depends on your GMT version (6.x recommended). Always validate outputs with synthetic data (e.g., using gmt grdmath to create a dipole field) before processing real observations. Always validate outputs with synthetic data (e