Regions of maximum shear and tension-compression stresses in the Martian interior have been revealed using the three-level compensation model. Nonequilibrium relief, density anomalies at the crust−mantle boundary, and density anomalies at the base of the lithosphere are the sources of the anomalous gravitational field. The thickness of elastic lithosphere positioned on a weak layer that has partially lost its elastic properties varies from 150 to 500 km. The weakening of the layer under the lithosphere is simulated
by a tenfold lower value of the shear modulus down to the core boundary. In general, the stresses for the threelevel compensation model differ from the values obtained for the two-level model (nonequilibrium relief and density anomalies at the crust−mantle boundary are the sources of the anomalous gravitational field) by 5−10%. Considerable differences between the models of two-level and three-level compensation are revealed beneath Hellas and Argyre regions.