Corner plot showing the distribution of ISM properties (Galactocentric radius R, midplane density n_{H, mid}, velocity dispersion σ, vertical scale height z_scale and vertical velocity gradient ∂_zv_z) at the SNR locations averaged over various length scales. We average over quadratic apertures with side lengths L = 0.2,  0.5,  1.0,  1.5 and 2.0 kpc. As expected the mean of n_{H, mid} decreases with R, with large scatter at R > 2 kpc. Compared to the ICs, the mean density profile has steepened. The velocity dispersion is roughly constant σ ∼ 10 km/s throughout the disk with considerable scatter. The disk scale height follows the scaling behavior predicted by vertical hydrostatic equilibrium of the gas $z_{\text{scale}}\propto \sigma /n_{\text{H}}^{0.5}$, but with a slightly lower normalization, likely due to the presence of stars. The vertical velocity gradient on average is zero with increasing scatter towards larger R. The maximum values reveal a preferred scale of $\left|{\partial }_z{v}_z\right|\sim \sigma /z_{\text{scale}}\sim t_{\text{ff}}^{-1}$ indicating that these expansion and contraction motions correspond to gravitational breathing modes of the disk. The velocity gradients peak around n_{H, mid} ∼ 1 cm⁻³, namely, the breathing of the disk coincides with the presence of multi-phase gas (see Fig. A.1).