The Indian National Gas Hydrate Program (NGHP) Expedition 01 presented the unique opportunity to constrain regional heat flow derived from seismic observations by using drilling data in three regions on the continental margin of India. The seismic bottom simulating reflection (BSR) is a well-documented feature in hydrate bearing sediments, and can serve as a proxy for apparent heat flow if data are available to estimate acoustic velocity and density in water and sediments, thermal conductivity, and seafloor temperature. Direct observations of temperature at depth and physical properties of the sediment obtained from drilling can be used to calibrate the seismic observations, decreasing the uncertainty of the seismically-derived estimates. Anomalies in apparent heat flow can result from a variety of sources, including sedimentation, erosion, topographic refraction and fluid flow. We constructed apparent heat flow maps for portions of the Krishna-Godavari (K-G) basin, the Mahanadi basin, and the Andaman basin and modeled anomalies using 1-D conductive thermal models.Apparent heat flow values in the Krishna-Godavari (K-G) basin and Mahanadi basin are generally 0.035 to 0.055 watts per square meter (W/m 2 ). The borehole data show an increase in apparent heat flow as water depth increases from 900 to 1500 m. In the SW part of the seismic grid, 1D modeling of the effect of sedimentation on heat flow shows that ~50% of the observed increase in apparent heat flow with increasing water depth can be attributed to trapping of sediments behind a "toe-thrust" ridge that is forming along the seaward edge of a thick, rapidly accumulating deltaic sediment pile. The remainder of the anomaly can be explained either by a decrease in thermal conductivity of the sediments filling the slope basin or by lateral advection of heat through fluid flow along stratigraphic horizons within the basin and through flexural faults in the crest of the anticline. Such flow probably plays a role in bringing methane into the ridge formed by the toe-thrust. Because of the small anomaly due to this process and the uncertainty in thermal conductivity, we did not model this process explicitly. In the NE part of the K-G basin seismic grid, a number of local heat flow lows and highs are observed, which can be attributed to topographic refraction and to local fluid flow along faults, respectively. No regional anomaly can be resolved. Because of lack of continuity between the K-G basin sites within the seismic grid and those ~70 km to the NE in water depths of 1200 to 1500 m, we do not speculate on the reason for higher heat flow at these depths. The Mahanadi basin results, while limited in geographic extent, are similar to those for the K-G basin.The Andaman basin exhibits much lower apparent heat flow values, ranging from 0.015 to 0.025 W/m 2 . Heat flow here also appears to increase with increasing water depth. The very low heat flow here is among the lowest heat flow observed anywhere and gives rise to a very thick hydrate stability zone in the...