Проведен анализ негидростатических напряжений в недрах Марса для двух типов моделей неоднородной упругости: модели с литосферой и модели с литосферой и возможными областями подплавления в ней. Численное моделирование системы уравнений упругого равновесия гравитирующей планеты проведено с шагом 1×1° по широте и долготе до глубины 1000 км. Граничными условиями задачи служат данные топографии и гравитационного поля планеты, определяемые по отношению к референсной поверхности, за которую принимается равновесный референсный сфероид. В качестве критерия выбора возможных эпицентров марсотрясений приняты большие значения максимальных касательных напряжений на фоне существенных растягивающих напряжений. Независимо от типа модели неоднородной упругости, зоны высоких сдвиговых и растягивающих напряжений в коре и мантии выявлены под ударными бассейнами Эллада и Аргир, равнинами Ацидалийское море, Аркадия и долиной Маринера.