A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(6) m(3) s(-1)) for the overturning and 1.8 and 1.4 (x 10(15) W) for heat transport near 10 degrees S and 10 degrees N, respectively. A dynamical decomposition of the overturning and heat transport shows that the rime-varying Ekman flow plus its barotropic compensation can explain a large parr of the seasonal variat...