A three-dimensional density model of the crust and uppermost mantle is determined by the inversion of a set of GOCE gravity and gradients residual anomalies beneath the eastern margin of the Tibetan Plateau and its surrounding areas.In our work,we choose five independent gravity gradients(Txx,Tzz,Txy,Txz,Tyz) to perform density inversion.Objective function is given based on Tikhonov regularization theory.Seismic S-wave velocities play the role of initial constraint for the inversion based on a relationship between density and S-wave velocity.Damped Least Square method is used during the inversion.The final density results offer some insights into understanding the underlying geodynamic processes:(1)Low densities in the margin of the Tibet,along with low wave velocity and resistivity results,yield conversions from soft and weak Tibet to the hard and rigid cratons.(2)The lowest densities are found in the boundary of the plateau,instead of the whole Tibet indicates that the effects of extrusion stress environment in the margin affect the changes of the substance there.The substances and environments conditioning for the earthquake preparations and strong deformation in this transitional zone.(3)Evident low-D anomaly in the upper and middle crust in the Lasha terrane and Songpan-Ganzi terrane illustrated the eastward sub-ducted of southeastern Tibet,which could be accounts for the frequent volcano and earthquakes there. A three-dimensional density model of the crust and uppermost mantle is determined by the inversion of a set of GOCE gravity and gradients residual anomalies beneath the eastern margin of the Tibetan Plateau and its surrounding areas. In our work, we choose five independent gravity gradients (Txx, Tzz, Txy, Txz, Tyz) to perform density inversion. Objective function is given based on Tikhonov regularization theory. Seismic S-wave velocities play the role of initial constraint for the inversion based on a relationship between density and S-wave velocity.Damped Least Square method is used during the inversion. The final density results offer some insights into understanding the underlying geodynamic processes: (1) Low densities in the margin of the Tibet, along with low wave velocity and resistivity results, yield conversions from soft and weak Tibet to the hard and rigid cratons. (2) The lowest densities are found in the boundary of the plateau, instead of the whole Tibet said that the effects of ext rusion stress environment in the margin impact the changes of the substance there.The substances and environments conditioning for the earthquake preparations and strong deformation in this transitional zone. (3) Evident low-D anomaly in the upper and middle crust in the Lasha terrane and Songpan-Ganzi terrane illustrated the eastward sub-ducted of southeastern Tibet, which could be accounts for the frequent volcano and earthquakes there.