The proper orthogonal decomposition (POD) method is used to construct a set of basis functions for spanning the ensemble of data in a certain least squares optimal sense. Compared with the singular value decomposition (SVD), the POD basis functions can capture more energy in the forecast ensemble space and can represent its spatial structure and temporal evolution more effectively. After the analysis variables are expressed by a truncated expansion of the POD basis vectors in the ensemble space, the control variables appear explicitly in the cost function, so that the adjoint model, which is used to de-rive the gradient of the cost function with respect to the control variables, is no longer needed. The application of this new technique significantly simplifies the data assimilation process. Several as-similation experiments show that this POD-based explicit four-dimensional variational data assimila-tion method performs much better than the usual ensemble Kalman filter method on both enhancing the assimilation precision and reducing the computation cost. It is also better than the SVD-based ex-plicit four-dimensional assimilation method, especially when the forecast model is not perfect and the forecast error comes from both the noise of the initial filed and the uncertainty of the forecast model.