Рассматривается новый алгоритм сейсмической томографии, основанный на совместном использовании данных объемных и поверхностных волн. Алгоритм является модификацией существующего кода LOTOS [Koulakov, 2009], разработанного для томографической инверсии по данным объемных волн от локальных землетрясений. Входные данные предлагаемой программы - времена пробега Р- и S-волн и дисперсионные кривые групповых скоростей волн Рэлея и Лява. Преимущества одновременного использования волн разных типов связаны с тем, что их свойства являются взаимодополняющими. В частности, объемные волны имеют хорошее разрешение на значительных глубинах, где расположены очаги землетрясений, а разрешение для поверхностных волн лучше на малых глубинах. Данные по поверхностным волнам выделяются из сейсмического шума и не зависят от распределения землетрясений, что позволяет более надежно контролировать геометрию системы наблюдения. При инверсии вклад данных разных типов определяется весами при соответствующих уравнениях. Представляемый алгоритм совместной инверсии был использован при анализе данных, собранных для региона кальдеры Тоба (Северная Суматра, Индонезия) на базе двух временных сетей сейсмических станций - IRIS и PASSCAL (1995 г.); GFZ и LAKE TOBA (2008 г.).
