To track metamorphic devolatilization along the top of subducted slabs, we adopted the geotherms 12 for the subduction zones of northwestern and southeastern Japan (Fig. 1).