Systems and methods are provided for generating a model for detection of seismic events. In this regard, one or more processors may receive from one or more stations located along an underwater optical route, one or more time series of polarization states of a detected light signal during a time period. The one or more processors may transform the one or more time series of polarization states into one or more spectrums in a frequency domain. Seismic activity data for the time period may be received by the one or more processors, where the seismic activity data include one or more seismic events detected in a region at least partially overlapping the underwater optical route. The one or more processors then generate a model for detecting seismic events based on the one or more spectrums and the seismic activity data.