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.
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
DELTA_T=15 # years gmt max net GPS_cascadia.xyz -T$DELTA_T -R-125/-122/42/49 -I0.2deg > cascadia_max.txt gmt max net gps_data
Note: Values are highly dependent on Δt and regional tectonic setting. Data: 500 GPS stations, 2008–2023 Command: Always validate outputs with synthetic data (e
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)
