Messier 8 (M8) is one of the brightest HII regions in the sky. We collected an extensive dataset comprising multiple submillimeter spectral lines from neutral and ionized carbon and from CO. Based on this dataset, we aim to understand the morphology of M8 and that of its associated photodissociation region (PDR) and to carry out a quantitative analysis of the physical conditions of these regions such as kinetic temperatures and volume densities. We used the Stratospheric Observatory For Infrared Astronomy (SOFIA), the Atacama Pathfinder Experiment (APEX) 12m, and the Institut de Radioastronomie Millimetrique (IRAM) 30m telescopes to perform a comprehensive imaging survey of the emission from the fine structure lines of CII and CI and multiple rotational transitions of carbon monoxide (CO) isotopologs within 1.3x1.3pc around the dominant Herschel 36 (Her 36) system, which is composed of at least three massive stars. To further explore the morphology of the region, we compared archival infrared, optical, and radio images of the nebula with our newly obtained fine structure line and CO data, and in particular with the velocity information these data provide. We performed a quantitative analysis, using both LTE and non-LTE methods to determine the abundances of some of the observed species, kinetic temperatures, and volume densities. Bright CO, CII and CI emission have been found toward the HII region and the PDR in M8. Our analysis places the bulk of the molecular material in the background of the nebulosity illuminated by the bright stellar systems Her 36 and 9 Sagitarii. Since the emission from all observed atomic and molecular tracers peaks at or close to the position of Her 36, we conclude that the star is still physically close to its natal dense cloud core and heats it. A veil of warm gas moves away from Her 36 toward the Sun and its associated dust contributes to the foreground extinction in the region. One of the most prominent star forming regions in M8, the Hourglass Nebula, is particularly bright due to cracks in this veil close to Her 36. We obtain H_2_ densities ranging from ~10^4^-10^6^cm^-3^ and kinetic temperatures of 100-150K in the bright PDR caused by Her 36 using radiative transfer modeling of various transitions of CO isotopologs.

Cone search capability for table J/A+A/615/A158/list (List of fits data cubes)