We perform the stress field analysis of different zones on the surface and in the lithosphere of Mars in order to reveal the areas of maximum shear and extensional stresses as potential marsquakes’ sources. Joint analysis of recent gravity and topography data truncated to degree 90 spherical harmonic in a frame of the static method is applied. The source of gravity anomalies is assumed to be the topographic loading and density anomalies at the crust-mantle boundary. Only non-equilibrium components of gravity and topography fields are considered: an outer surface of the hydrostatic equilibrium form of the planet is taken as a reference surface for topography and gravity field of Mars. Numerical calculations of extension-compression stresses and maximum shear stresses are carried out for an interior structure model of Mars having 150-300 km thick lithosphere, with a 1х1 arc-deg spherical grid and up to 1000 km depth.