The surface heat budget over the Weddell Sea ice cover in 1996 was studied on the basis of data from Argos buoys equipped with meteorological sensors. In addition, a thermodynamic sea ice model, satellite-based data on the sea ice concentration, sonar results on ice thickness distribution, and output from large-scale meteorological models were all utilized. Applying the buoy data, the sensible heat flux over sea ice was calculated by Monin-Obukhov theory using the gradient method, and the latent heat flux was obtained by the bulk method. A second estimate for the surface fluxes was obtained from the thermodynamic sea ice model, which was forced by the buoy observations. The results showed a reasonable agreement. The dominating component in the heat budget over ice was the net longwave radiation, which had a mean annual cooling effect of -28 W m-2. This was balanced by the net shortwave radiation (annual mean 13 W m-2), the sensible (13 W m-2) and latent (-3 W m-2) heat fluxes, and the conductive heat flux through the ice (5 W m-2). The regional surface fluxes over the fractured ice cover were estimated using the buoy data and Special Sensor Microwave Imager (SSMI)-derived ice concentrations. In winter the regional surface sensible heat flux was sensitive to the ice concentration and thickness distribution. The estimate for the area-averaged formation rate of new ice in leads in winter varies from 0.05 to 0.21 m per month depending on the SSMI processing algorithm applied. Countergradient fluxes occurred 8--10% of the time. The buoy observations were compared with the operational analyses of the European Centre for Medium-Range Weather Forecasts (ECMWF) and the reanalyses of the National Centers for Environmental Prediction (NCEP)/National Center for Atmospheric Research (NCAR). The 2 m air temperature and surface temperature were 3.5¿ and 4.4¿C too high, respectively, in the ECMWF and 3.2¿ and 3.0¿C too low in the NCEP/NCAR fields, but the models reproduced the synoptic-scale temperature variations well. The errors seem to be related to the cloud cover and the surface boundary conditions. Neither of the models recognizes leads in the ice pack, and the ice and snow thicknesses are often far from reality. The distribution of the cloud cover in the both models differed a lot from observations.