There are four parts mainly involved: 1)development and estimation on the International Terrestrial Reference Frame (ITRF), the main characteristics of the most new ITRF2005; 2)definition, construction and maintenance of a local geodetic reference coordinate frame, present situation and development of the local geodetic reference coordinate frame in Europe, America and north America; 3)the modernization of the vertical datum, the progress in the establishment of a modern vertical datum in United States of America and Europe, and the way how to establish a global vertical datum; 4)some special problems should be concerned in the construction of a geodetic reference coordinate frame, for example, the contribution of VLBI and SLR, the way to reduce the influence of the periodical variation of the frame points on the precise positioning, and the improvement suggestion for the computation of zero tide system.

The Hori and YabukiMatu’ura inverse methods are used to estimate the distribution of interseismic slip deficit
on the border land between Philippine Sea plate and the southwest Japan with GPS data from Jan.,1997 to Dec.,1999.
The estimated results show that maximum slip rate is about 7 cm/a, which is consistent with the Philippine Sea plate convergence rate.
The 50% contribution of strong coupling is confined to the depths of 10-30 km.
By compring the two results it can be seen that there are some different effects caused by the collision between eastern and central
Japan, because of different analytic process and constrain condition. It can be concluded that Hori inverse method is more applicable and impersonal to estimate slip distribution in the area with dense observation networks,
such as the southwest Japan.

　 Ｔｈｅ ｐｒｅｓｅｎｔ-ｄａｙ ｃｒｕｓｔａｌ ｄｅｆｏｒｍａｔｉｏｎ ｉｎ ｔｈｅ ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ ｉｓ ｓｔｕｄｉｅｄ ｗｉｔｈ ＧＰＳ ｄａｔａ ｏｆ ３７ｓｔａｔｉｏｎｓ ｄｕｒｉｎｇ １９９７－２００５． Ｔｈｅ ｖｅｌｏｃｉｔｉｅｓ ｒｅ-ｌａｔｅｄ ｔｏ Ｅｕｒａｓｉａ ｉｓ ｇｅｎｅｒａｌｌｙ ｓｍａｌｌ， １．５８ｍｍ／ａｏｎ ａｖｅｒａｇｅ， ｉｔ ｉｓ ｓｈｏｗｉｎｇ ｔｈａｔ ｔｈｅ ｒｅｇｉｏｎ ｉｓ ｇｅｎｅｒａｌｌｙ ｓｔａｂｌｅ． Ｈｏｒｉｚｏｎｔａｌ ｖｅｌｏｃｉｔｉｅｓ ｏｆ ５ｓｔａｔｉｏｎｓ ａｐｐｅａｒ ｔｏ ｂｅ ｓｉｇｎｉｆｉｃａｎｔ ａｔ ９５％ ｃｏｎ-ｆｉｄｅｎｃｅｌｅｖｅｌ， ｈｏｗｅｖｅｒ ｏｔｈｅｒ ｓｔａｔｉｏｎｓ ｓｈｏｗ ｉｎ-ｓｉｇｎｉｆｉｃａｎｔ ｍｏｔｉｏｎ． Ｓｅｖｅｒａｌ ｓｔａｔｉｏｎｓ ｉｎ ｔｈｅ ｓｏｕｔｈｗｅｓｔｅｒｎ ｃｏｒｎｅｒ ｏｆ ｔｈｅ ｓｔｕｄｉｅｄ ｒｅｇｉｏｎ ａｒｅｃｏｎｓｉｓｔｅｎｔｌｙ ｍｏｖｉｎｇ ｔｏｗａｒｄｓ ｔｈｅ ｓｏｕｔｈｅａｓｔ， ｔｈｉｓ ｉｓ ｍｏｓｔ ｐｒｏｂａｂｌｙ ａ ｄｉｓｔａｎｔ ｅｆｆｅｃｔ ｏｆ ｔｈｅ ｅａｓｔｗａｒｄ ｍｏｔｉｏｎ ｏｆ Ｑｉｎｇｈａｉ-Ｔｉｂｅｔ ｐｌａｔｅａｕ． Ｓｍａｌｌ ｖｅｌｏｃｉｔｉｅｓ ｏｆ ｔｈｅ ｓｔａｔｉｏｎｓ ｎｅａｒ ｔｈｅ ｅｐｉ-ｃｅｎｔｅｒｓ ｏｆ ｔｈｅ １９９９ Ｍｗ７．１ ａｎｄ ２００２ Ｍｗ７．３ ｄｅｅｐ-ｆｏｃｕｓ ｅａｒｔｈｑｕａｋｅｓ ｉｎ Ｗａｎｇｑｉｎｇ Cｏｕｎｔｙ，ｓｈｏｗ ｔｈａｔ ｔｈｅ ２ ｅａｒｔｈｑｕａｋｅｓ ｄｉｄ ｎｏｔ ｃａｕｓｅ ｂｒｏａｄ ｄｅ-ｆｏｒｍａｔｉｏｎ ｏｆ ｓｈａｌｌｏｗ ｃｒｕｓｔ．Ｔｈｅ ｈｏｒｉｚｏｎｔａｌ ｖｅｌｏｃｉｔｉｅｓ ｏｆ ｔｈｅ ｓｔａｔｉｏｎｓ ｉｎ ｔｈｅ Ｎｏｒｔｈｅａｓｔ- ｅｒｎ Ｈｅｉｌｏｎｇｊｉａｎｇ ｐｒｏｖｉｎｃｅ ｍａｙ ｒｅｆｌｅｃｔ ｔｈｅ ｅｆ-ｆｅｃｔ ｏｆ ｔｈｅ Ｎｏｒｔｈ Ａｍｅｒｉｃａ ｐｌａｔｅ’ｓ ｉｎ ｖａｓｉｏｎ ｏｎｔｈｅ Ｅｕｒａｓｉａ ｐｌａｔｅ．

The tidal factors from the EW and NS components of ground tilt in Fujian from 1990 to 2004

are
calculated. The fulltime plots for the characteristic quantities (dilatancy, shearing strain, media anisotropy) are drawn and their relations with midstrong earthquake are analyzed. The results show that the ocean tide is a main factor to determine the characteristics of tilt tide in the shore area of Fujian and the dynamic combination characteristic quantities of tilt tide correspond to Fujian earthquake (M＞4.5) or Taiwan earthquake (M＞7.0) well.

The statistic properties of parametric components of least

square kernel estimator of nonlinear semiparametric models under certain regular conditions is analyzed, the results show that the estimation of parametric component is biased and variance is asymptotic normal,finally the asymptotic variance formulea are deduced.

On the basis of the framework of L ^P norm, the mathematic model of classic phase unwrapping methods is studied. Phase unwrapping methods are classed into three types. The topographic SAR interferogram in Bam region of Iran is analyzed with these phase unwrapping methods.The results show that the Goldstein filtering method can improve effectively the distribution of residuals and branch cuts and raise the visual effects of the interferogram as well as the ratio of signal to noise.
The L ¹norm methods based on network currents can achieve an optimal resolving, but the efficiency is rather low, the

L ²norm method can obtain better optimal resoloring and has better efficiency,

L ⁰norm algorithm performed best in efficiency, but with discontinuity and uncertainty.

Several expressions of horizontal displacements have been analyzed. A new idea in which the horizontal displacements are represented based on the new type of geodetic coordinate system is put forward. According to the transformation formulas between the new geodetic coordinate and the geodetic longtude and latitude, the horizontal displacement components along the northern(southern) and eastern(western) direction on the ellipsoidal surface can be expressed with the new system.Taking the length as the coordinate parameter,the geodetic computation can be realized simply and conveniently and the first fundamental form on the ellipsoidal surface is brief in the new geodetic coordinate system, so the planimetric position and horizontal displacements can be expressed in the same type of coordinate system. It would be beneficial to the strain field analysis.

An analysis of influence from the uncorrected satellite orbit errors on the GPS/VRS network based positioning is carried out. A precise realtime method for using the IGS ultra rapid products (IGU) and the GPS broadcast ephemeris to calculate the VRS orbit corrections is presented which is suited for the realtime correction of precise orbit of VRS in
GPS/VRS reference station network based positioning. Test data acquired from both the SGRSN (Sichuan GPS Reference Station Network) and SCIGN (Southern California integrated GPS network) are used to evaluate the performance of the modeling techniques. The new method is proved to be more precise and reliable as compared with the existing conventional networkbased orbit error interpolation method. It is shown that 0.004×10 ^-6

relative accuracy is reached, i.e. the influence from the orbit bias for the RTK positioning within 100km area can be of submillimeter level.