This study introduces HARPA, a high-rate seismic phase association framework based on an optimal-transport metric. It jointly estimates earthquake parameters and wave-speed structure from dense, overlapping arrivals and opens a path to passive-source seismic tomography.