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.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 max net
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) | DELTA_T=15 # years gmt max net GPS_cascadia
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