Lithosphere density structure beneath the eastern margin of the Tibetan Plateau and its surrounding areas derived from GOCE gradients data

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 earthouakes there.

Improved Implicit SPH Method for simulating free surface flows of power law fluids

Implicit smoothed particle hydrodynamics method has been proposed to overcome the problem that only very small time steps can be used for high viscosity fluids(such as power law fluids)in order to obtain a stable simulation.However,the pressure field is difficult to simulate correctly with this method because the numerical high-frequency noise on the pressure field cannot be removed.In this study,several improvements,which are the diffusive term in the continuity equation,the artificial viscosity and a simplified physical viscosity term in the momentum equation,are introduced,and a new boundary treatment is also proposed.The linear equations derived from the momentum equation are large-scale,sparse and positive definite but unsymmetrical,therefore,Conjugate Gradient Squared(CGS)is used to solve them.For the purpose of verifying the validity of the proposed method,Poiseuille flows with Newtonian and power law fluids are solved and compared with exact solution and traditional SPH.Drops of different fluid properties impacting a rigid wall are also simulated and compared with VOF solution.All the numerical results obtained by the proposed method agree well with available data.The proposed method shows the higher efficiency than traditional SPH and the less numerical noise on the pressure field and better stability than implicit SPH.