A four-dimensional variational method has been developed that assimilates expandable bathythermograph and Tropical Atmosphere-Ocean subsurface temperature data into the Hamburg Ocean Primitive Equation model. The method decreases the misfit between model and observed ocean temperatures by adjusting the surface forcing. The main goal of the assimilation scheme is to improve ocean analyses in the tropical Pacific. As a first study, only wind stress is adjusted in the assimilation. In two identical twin experiments it is demonstrated that the scheme works well in the equatorial Pacific. The scheme is capable of reducing errors in the ocean analysis that originate either from the wind stress forcing or the initial state. The impact of model errors on the data assimilation is investigated in a experiment with real observations. In this experiment the temperature innovations near the equator are comparable to those of an optimal interpolation data assimilation scheme. ¿ 2001 American Geophysical Union