(9) Time-series along tracks

A common application in many scientific disciplines involves plotting one or several time-series as as “wiggles” along tracks. Marine geophysicists often display magnetic anomalies in this manner, and seismologists use the technique when plotting individual seismic traces. In our example we will show how a set of Geosat sea surface slope profiles from the south Pacific can be plotted as “wiggles” using the pswiggle program. We will embellish the plot with track numbers, the location of the Pacific-Antarctic Ridge, recognized fracture zones in the area, and a “wiggle” scale. The Geosat tracks are stored in the file tracks.txt, the ridge in ridge.xy, and all the fracture zones are stored in the multiple segment file fz.xy. The profile id is contained in the segment headers and we wish to use the last data point in each of the track segments to construct an input file for pstext that will label each profile with the track number. We know the profiles trend approximately N40E so we want the labels to have that same orientation (i.e., the angle with the baseline must be 50). We do this by extracting the last record from each track and select segment header as textstring when running the output through pstext. Note we offset the positions by -0.05 inch with -D in order to have a small gap between the profile and the label:

#!/bin/bash
#		GMT EXAMPLE 09
#
# Purpose:	Make wiggle plot along track from geoid deflections
# GMT modules:	gmtconvert, pswiggle, pstext, psxy
# Unix progs:	
#
ps=example_09.ps
gmt pswiggle tracks.txt -R185/250/-68/-42 -K -Jm0.13i -Ba10f5 -BWSne+g240/255/240 -G+red \
	-G-blue -Z2000 -Wthinnest -S240/-67/500/@~m@~rad --FORMAT_GEO_MAP=dddF > $ps
gmt psxy -R -J -O -K ridge.xy -Wthicker >> $ps
gmt psxy -R -J -O -K fz.xy -Wthinner,- >> $ps
# Take label from segment header and plot near coordinates of last record of each track
gmt convert -El tracks.txt | gmt pstext -R -J -F+f10p,Helvetica-Bold+a50+jRM+h \
	-D-0.05i/-0.05i -O >> $ps

The output shows the sea-surface slopes along 42 descending Geosat tracks in the Eltanin and Udintsev fracture zone region in a Mercator projection.

../_images/example_09.png

Time-series as wiggles along a track.