Ecological footprint method is an effective approach for quantitative evaluatation of the regional sustainable development. Dynamic changes of ecological footprint(EF) and ecological capacity(EC) were analyzed based on ecological footprint method and land use sustainability was assessed using EF, EC and other indicators of Xuzhou City, Jiangsu Province from 1999 to 2008. The results showed that the per capita ecological footprint of Xuzhou City increased from 1.13 hm² in 1989 to 2.57 hm² in 2008, and the per capita ecological capacity decreased from 0.58 hm² in 1989 to 0.45 hm² in 2008. The per capita ecological deficit increased from 0.55 hm² in 1989 to 2.10 hm² in 2008. Intensification of land-use is low in Xuzhou City. The author used a linear regression model to predict the population, and predicted the ecological footprint of 2010-2015 by gray prediction model. The multi-objective linear programming model was built to optimize the land use structure in Xuzhou City, which minimizes ecological deficit and maximizes economic benefits as the objective function, and sets up seven variables and nine constraints. The land use structure was optimized in Xuzhou City according to optimal scheme.