A kinematic statistical method is proposed to determine the position for Chang'E-3 (CE-3) lunar lander. This method uses both ranging and VLBI measurements to the lander for a continuous arc, combing with precise knowledge about the motion of the moon as provided by planetary ephemeris, to estimate the lander's position on the lunar surface with high accuracy. Accuracy analyses are carried out with simulation data using the software developed at Shanghai Astronomical Observatory in this study to show that measurement errors will dominate the position accuracy. Application of lunar digital elevation model (DEM) as constraints in the lander positioning is also analyzed. Simulations show that combing range/doppler and VLBI data, single epoch positioning accuracy is at several hundred meters level, but with ten minutes data accumulation positioning accuracy is able to be achieved with several meters. Analysis also shows that the information given by DEM can provide constraints in positioning, when DEM data reduce a 3-dimensional positioning problem to 2-dimensional. Considering the Sinus Iridum, CE-3 lander's planned landing area, has been observed with dedicated details during the CE-1 and CE-2 missions, and its regional DEM model accuracy may be higher than global models, which will certainly support CE-3's lander positioning. The first Chinese lunar explorer CE-1 was launched on October 24, 2007, and following CE-2 was launched on October 1, 2010 [1][2][3]. According to the long-term schedule, Chinese unmanned lunar exploration plan is divided into three stages: orbiting, landing and returning stages. In the orbiting stage, the CE-1 and CE-2 satellites fly around the moon and take photos about the landing area. In the landing stage, the main task is to softland on the lunar surface and to explore automatically [4]. A lunar lander and a rover are needed in this stage. How to determine the positions of the lander and rover on the lunar surface is very important to the following Chinese lunar exploration project. The rover moves very slowly and not far away from the lander, so the position determination of the lander is more critical. Dealing with position determination of a satellite, the dynamic statistical orbit determination (OD) method and single point positioning are widely used in lunar exploration [1][2][3]5,6]. OD can be viewed as a special case of general parameter estimation problem, where the parameters characterizing the orbit of a satellite, have to be determined from observations using the least squares method. The main error sources are from observations as well as the force model, and in order to get high OD accuracy, forces exerting on the satellite must be modeled precisely. Comparing with dynamic OD method, single point positioning does not need to model forces, and the position of the satellite can be