\ We calculated statistical average of thermal data to speculate regional thermal structure of the forearc area of the Japanese Islands. The three thermal statistical averages show a difference of a high thermal regime in the western part of forearc inner zone and a low in the Kanto forearc outer zone. The Kanto zone marks 18 K km(-1) for mean geothermal gradient, 44 MW m(-2) for mean heat flow, while the western inner zone shows 27 K km-1 for mean geothermal gradient, 63 MW m(-2) for mean heat flow. The geothermal gradients of the Nobi Plain and the Osaka Plain in the western inner zone are 29 and 36 K km(-1), respectively, while the value of the Kanto Plain in the Kanto zone is 21 K km(-1). Taking account of the effect of accumulation of sediments, we see the difference in the thermal regime between the plains and conclude that the difference is significant. Heat flux in the crust depends on the volume of granite rich in radioactive elements. There are few granitic rocks in the Kanto zone, while granitic rocks are dominant in the western inner zone. The heat flow of 20 mW m(-2) is attributed to the granitic rocks of about 8 km in thickness. There are two oceanic plate subductions of the Pacific plate and the Philippine Sea plate under the Kanto zone, while only the Philippine Sea plate has been subducting under the western inner zone. The model simulation based on thermal and subduction model shows a heat flow ranging 50-60 mW m(-2) in the southwest Japan forarc area and a low value of about 20 mW m-2 in the northeast Japan forearc area. The heat flux from the cooling oceanic lithosphere depends on the age of plate. The Shikoku Basin, a part of the Philippine Sea plate, off the western inner zone is 15-30 Ma, while the Pacific plate off the Kanto zone is 122-132 Ma. Theoretically, heat flux values of 15 and 50 Ma oceanic plates range 60-120 mW m(-2) and those of 122-132 Ma could be about IOMWM-2. If the heat flux contribution from the Philippine Sea plate under the Kanto zone is smaller than the plate under the western inner zone, there could be a thermal regime difference in order of several tens of mW m(-2). Conclusively, the cause of the difference of heat flux could be the uneven granitic rocks distribution and/or the difference of heat flux between the two subducting plate. (c) 2005 Published by Elsevier B.V.