Oceanic heat flux (Fw) is the vertical heat flux that is transmitted to the base of sea ice. It is the main source of sea ice bottom melting. The residual method was adopted to study oceanic heat flux under sea ice. The data acquired by 28 ice mass balance buoys (IMBs) deployed over the period of 2004 to 2013 in the Arctic Ocean were used. Fw values presented striking seasonal and spatial variations. The average summer Fw values for the Canada Basin, Transpolar Drift, and Multiyear Ice area were 16.8, 7.7, and 5.9 Wm−2, respectively. The mean summer Fw for the whole Arctic was 10.1 Wm−2, which was equivalent to a bottom melt of 0.4m. Fw showed an autumn peak in November in the presence of the near-surface temperature maximum (NSTM). The average Fw for October to December was 3.7Wm−2. And the average Fw for January to March was 1.0Wm−2, which was approximately one third of the average Fw in the presence of NSTM. The summer Fw was almost wholly attributed to the incident solar radiation that enters the upper ocean through leads and the open water. Fw calculated through the residual method using IMB data was compared with that calculated through the parameterization method using Autonomous Ocean Flux Buoy data. The results revealed that the Fw provided by the two methods were consistent when the sea ice concentration exceeded 70% and mixing layer temperature departure from freezing point was less than 0.15℃. Otherwise, the Fw yielded by the residual method was approximately one third smaller than that provided by the parameterization method.