Altimetry data and the elastic stress tensor of subduction zones

The stress field in the lithosphere caused by the distribution of density anomalies associated to the geoidal undulations observed by the GEOS-3 and SEASAT Earth satellites in the Tonga region was studied. Different models of the lithosphere were generated with different assumptions on the density distribution and geometry, all generating a geoid profile almost identical to the observed one. The first model is the Airy isostatic hypothesis which consists of a crust of density 2.85 laying on a lithosphere of density 3.35. The models obtained with different compensation depths give residual shortwavelength anomalies of the order of several tens of mgal and several tens of meters geoidal undulations. It indicates that there is no isostasy of the Airy type in the Tonga region because the observed geoid has very smooth undulation of about 25 m over a distance of 2000 km. The Pratt isostatic hypothesis is used in a model consisting of a crust of variable density laying on a lithosphere of higher density. This model gives smaller residual anomalies but still shows that there is no isostasy of the Pratt type in the Tonga region because the observed geoidal undulation are much smaller and smoother than the residual undulations associated to the Pratt model of isostasy.