We develop an algorithm for automatic identification of fault zone trapped waves in data recorded by seismic fault zone arrays. Automatic S picks are used to identify time windows in the seismograms for subsequent search for trapped waves. The algorithm calculates five features in each seismogram recorded by each station: predominant period, 1 s duration energy (representative of trapped waves)…