AbstractWe present results of Herschel PACS imaging spectroscopy data toward ten massive young stellar objects taken as part of the WISH project. Our sample consists of four high mass protostellar objects (HMPOs), two hot molecular cores (HMCs), and four ultracompact HII regions (UCHIIs), and the spectra cover a broad range of wavelengths (55 to 210 μm) imaged over an ~50” field with 5×5 spaxels. By fitting the continua utilizing a modified black-body formula we estimate mass-weighted dust temperature and column density distributions of warm dust and find that UCHII regions are hottest and HMCs are most deeply embedded. We also estimate rotational temperature and column density distributions of warm CO gas using the rotational diagram analysis, which are comparable over targets in contrast to continuum results. By comparing high J CO line fluxes to the RATRAN estimates of centrally heated envelope models, we find that majority of warm CO originates from bipolar outflow shocks.