A new data assimilation method called the explicit four-dimensional variational (4DVAR) method is proposed. In this method, the singular value decomposition (SVD) is used to construct the orthogonal basis vectors from a forecast ensemble in a 4D space. The basis vectors represent not only the spatial structure of the analysis variables but also the temporal evolution. After the analysis variables are expressed by a truncated expansion of the basis vectors in the 4D space, the control variables in the cost function appear explicitly, so that the adjoint model, which is used to derive the gradient of cost function with respect to the control variables, is no longer needed. The new technique significantly simplifies the data assimilation process. The advantage of the proposed method is demonstrated by several experiments using a shallow water numerical model and the results are compared with those of the conventional 4DVAR. It is shown that when the observation points are very dense, the conventional 4DVAR is better than the proposed method. However, when the observation points are sparse, the proposed method performs better. The sensitivity of the proposed method with respect to errors in the observations and the numerical model is lower than that of the conventional method. data assimilation, four-dimensional variation, explicit method, singular value decomposition, shallow water equationThe four-dimensional variational data assimilation (4DVAR) has been a very successful technique and used in operational numerical weather prediction (NWP) of some weather forecast centers [1,2] . In this method the optimal estimate of initial condition of a forecast model is obtained by fitting the forecasts to observations within a time window. The attractive features of 4DVAR include: (1) the full-model is set as a strong dynamical constraint, and (2) it has the ability to assimilate the data at multiple time. However, the control variables (initial state) are expressed implicitly in the cost function. In order to compute gradient of the cost function with respect to the control variables, one has to integrate the adjoint model of the forecast model. But coding the adjoint for the 4DVAR and maintaining the adjoint, updated with the model upgrading, are extremely labor-intensive, especially when the forecast model is nonlinear and the model physics contain parameterized discontinuities [3,4] . Some researchers try to avoid integrating the adjoint model or reducing the expensive computation [5][6][7] . But the linear or adjoint model is still required in the methods mentioned above. So the three-dimensional variational data assimilation (3DVAR) becomes the common practice in many numerical weather forecast centers. The 3DVAR can be considered as a simplification of the 4DVAR, but it has lost the two advantages of the 4DVAR mentioned above. In general, the analysis results rely heavily on the information from the background field. However, as we know, in 3DVAR the background error covariance matrix is usually simplified and not flo...