Panels a–c: ¹³CO integrated intensity maps (blue contours) superimposed on the Spitzer-IRAC 8 μm emission (gray). The integrated-velocity ranges are shown on top of each image. The pink arrows mark the pillars, while the red dashed arrows indicate two gas flows. The blue contour levels start from 0.2 K km s⁻¹(5σ) in a step of 1.6 K km s⁻¹. The pink dashed circles represent M 16, N19, and RCW165. The green star marks the NGC 6611 cluster ionizing M 16. Panels d–f: C¹⁸O integrated intensity images overlaid with the ATLASGAL 870 μm (red contours) and 8 μm emission (gray). The blue contour levels start from 0.15 K km s⁻¹(5σ) in a step of 0.75 K km s⁻¹, while the redcontour levels start from 0.06 Jy beam⁻¹ (3σ) in a step of 0.3 Jy beam⁻¹. The green dashed lines mark the directions of position–velocity diagrams in Fig. 10.