Terrestrial heat flux was measured in fjords, in boreholes, and in offshore wells at sites across the convergent margin of southwestern British Columbia from the continental shelf landward to the Garibaldi Volcanic Belt. Temperatures in the offshore wells were corrected for drilling disturbances, and formation thermal conductivities were modeled using measurements on cuttings and downhole geophysical logs. Marine measurements in the fjords were corrected for the large effects of refraction as well as aperiodic temperature variations in the bottom waters. There were excellent agreement between marine measurements and those from nearby onshore boreholes. The heat flux above the subducting Juan de Fuca plate steadily decreases landward from over 50 mW m-2 on the shelf to 25 mW m-2 seaward of the Garibaldi Volcanic Belt. An abrupt increase to 80 mW m-2 over a distance of 20 km is centered 30 km seaward of the volcanic zone. Very large variations in heat flux occur locally within the Pleistocene volcanic area, the result of advective cooling of intrusive magmas. The measured heat generation of crustal samples along the entire profile is low, 0.6--0.8 μW m-3. A landward dipping, seismically reflective zone above the subducting oceanic plate beneath Vancouver Island appears to be nearly profile is low, 0.6--0.8 μW m-3. A landward dipping, seismically reflective zone above the subducting oceanic plate beneath Vancouver Island appears to be nearly isothermal. It is postulated that dehydration of the subducting oceanic crust at above approximately 450¿C absorbs heat and produces water which flows updip along this zone in the overlying subduction complex, effectively redistributing the heat seaward to where the water is reabsorbed in hydration processes. A relatively cool crustal wedge lies above the deeper subducting oceanic crust, and at its thick, landward side the abrupt increase in surficial heat flux must caused by a shallow (10 km depth) heat source produced by ascending magma. ¿ American Geophysical Union 1988