# The SEDIGISM survey: Molecular cloud morphology

K. Neralwar, D. Colombo, A. Duarte-Cabral, J. Urquhart, M. Mattern, F. Wyrowski, K. Menten, P. Barnes, Á. Sánchez-Monge, A. Rigby, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

K. Neralwar, D. Colombo, A. Duarte-Cabral, J. Urquhart, M. Mattern, et al. The SEDIGISM survey: Molecular cloud morphology. Astronomy and Astrophysics - A&A, 2022, 664, pp.A84. 10.1051/0004-6361/202142513. [hal-03799782](https://hal.science/hal-03799782)

HAL Id: hal-03799782

https://hal.science/hal-03799782

Submitted on 6 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
# The SEDIGISM survey: Molecular cloud morphology

## II. Integrated source properties

K. R. Neralwar<fnref n="1" />,*, D. Colombo<fnref n="1" />, A. Duarte-Cabral<fnref n="2" />, J. S. Urquhart<fnref n="3" />, M. Mattern<fnref n="4" />, F. Wyrowski<fnref n="1" />, K. M. Menten<fnref n="1" />, P. Barnes<fnref n="5" />,<fnref n="6" />, Á. Sánchez-Monge<fnref n="7" />, A. J. Rigby<fnref n="2" />, P. Mazumdar<fnref n="1" />, D. Eden<fnref n="8" />, T. Csengeri<fnref n="9" />, C. L. Dobbs<fnref n="10" />, V. S. Veena<fnref n="7" />, S. Neupane<fnref n="1" />, T. Henning<fnref n="11" />, F. Schuller<fnref n="1" />,<fnref n="12" />, S. Leurini<fnref n="13" />, M. Wienen<fnref n="1" />, A. Y. Yang<fnref n="1" />, S. E. Ragan<fnref n="2" />, S. Medina<fnref n="1" />, and Q. Nguyen-Luong<fnref n="14" />

Received 22 October 2021 / Accepted 4 May 2022

## ABSTRACT

The Structure, Excitation, and Dynamics of the Inner Galactic Interstellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution ¹³CO (2–1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds (MCs). In a previous paper, we classified the SEDIGISM clouds into four morphologies. However, how the properties of the clouds vary for these four morphologies is not well understood. Here, we use the morphological classification of SEDIGISM clouds to find connections between the cloud morphologies, their integrated properties, and their location on scaling relation diagrams. We observe that ring-like clouds show the most peculiar properties, having, on average, higher masses, sizes, aspect ratios, and velocity dispersions, compared to other morphologies. We speculate that this is related to the physical mechanisms that regulate their formation and evolution; for example, turbulence from stellar feedback can often result in the creation of bubble-like structures. We also see a trend of morphology with the virial parameter, whereby ring-like, elongated, clumpy, and concentrated clouds have virial parameters in decreasing order. Our findings provide a foundation for a better understanding of MC behaviour, based on their measurable properties.

**Key words.** ISM: clouds – local interstellar matter – ISM: bubbles – submillimeter: ISM

## 1. Introduction

Molecular clouds (MCs) are some of the densest and coldest regions in the interstellar medium (ISM) and the exclusive sites of star formation in galaxies (Dobbs et al. 2014; Heyer & Dame 2015; Miville-Deschênes et al. 2017). They are turbulent structures with masses of the order of $10^2–10^7 M_\odot$ and most of this mass is concentrated in high-mass giant molecular cloud complexes (>10⁵ $M_\odot$ for the Milky Way) (Roman-Duval et al. 2010). Generally speaking, most MCs are magnetically super-critical, though not by a large margin (Crutcher et al. 2010; Crutcher 2012; Dobbs et al. 2014). Moreover, magnetic fields vary from

cloud to cloud, and some clouds may be magnetically supported (Chapman et al. 2011; Barnes et al. 2015). MCs span a range of sizes between ~1 and ~200 pc, and present a hierarchical structure with dense clumps and denser (proto-stellar or pre-stellar) cores (Blitz & Stark 1986; Rosolowsky et al. 2008; Ballesteros-Paredes et al. 2020). The surface densities of MCs have a wide range, in other words, ~1 to more than 1000 $M_\odot$ pc⁻<fnref n="2" /> (Barnes et al. 2018). A typical value for the surface densities of clouds is ~100 $M_\odot$ pc⁻<fnref n="2" /> (Roman-Duval et al. 2010; Rebolledo et al. 2012) and it may change with the Galactic environment (e.g. M51; Colombo et al. 2013; Sun et al. 2018).

The correlations between velocity dispersion and other cloud properties, such as surface density and the virial parameter, have led to various interpretations explaining their structure. The historically widely used Larson (1981) paper concluded that clouds

<footnote n="1">¹ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: kneralwar@mpi.fr-bonn.mpg.de</footnote>

<footnote n="2">² School of Physics & Astronomy, Cardiff University, Queen's Building, The Parade, Cardiff CF24 3AA, UK</footnote>

<footnote n="3">³ Centre for Astrophysics and Planetary Science, University of Kent, Canterbury CT2 7NH, UK</footnote>

<footnote n="4">⁴ Laboratoire d'Astrophysique (AIM), CEA, CNRS, Université Paris-Saclay, Université Paris Diderot, Sorbonne Paris Cité, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="5">⁵ Space Science Institute, 4765 Walnut St Suite B, Boulder, CO 80301, USA</footnote>

<footnote n="6">⁶ School of Science and Technology, University of New England, Armidale NSW 2351, Australia</footnote>

<footnote n="7">⁷ I. Physikalisches Institut, Universität zu Köln, Zülpicher Strasse 77, 50937 Cologne, Germany</footnote>

<footnote n="8">⁸ Astrophysics Research Institute, Liverpool John Moores University, IC2, Liverpool Science Park, 146 Brownlow Hill, Liverpool, L3 5RF, UK</footnote>

<footnote n="9">⁹ Laboratoire d'astrophysique de Bordeaux, CNRS, Univ. Bordeaux, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="10">¹⁰ School of Physics and Astronomy, University of Exeter, Stocker Road, Exeter EX4 4QL, UK</footnote>

<footnote n="11">¹¹ Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany</footnote>

<footnote n="12">¹² Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany</footnote>

<footnote n="13">¹³ INAF – Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047 Selargius (CA), Italy</footnote>

<footnote n="14">¹⁴ The American University of Paris, 2bis, Passage Landrieu 75007 Paris, France</footnote>

<footnote>* Member of the International Max Planck Research School (IMPRS) for Astronomy and Astrophysics at the Universities of Bonn and Cologne.</footnote>
are self-gravitating objects in virial equilibrium and the gravity acting on clouds is counterbalanced by internal forces, giving an average virial parameter around unity (Solomon et al. 1987; Fukui et al. 2008; Roman-Duval et al. 2010). Another work, Ballesteros-Paredes et al. (2011), argued that clouds are in free-fall collapse, but that the uncertainties and scatter in the size-linewidth relation allow for virial equilibrium to be a consistent mechanism. Field et al. (2011) present clouds as marginally self-gravitating pressure-confined objects. A more recent work, Vázquez-Semadeni et al. (2019), presents the global hierarchical collapse (GHC) scenario as a description of the processes involved throughout the life cycle of MCs, from their formation in the diffuse medium to their destruction by massive stars. It proposes a 'moderated' free-fall to maintain the cloud structure, which takes into account the interplay at different scales, and includes processes such as stellar feedback, which affect the overall cloud evolution.

The filamentary nature of MCs is ubiquitous in the Galaxy and most of the star-forming clumps are situated in large-scale filamentary structures (Schneider & Elmegreen 1979; Arzoumanian et al. 2011; Polychroni et al. 2013; Hacar et al. 2013; Li et al. 2016; Olmi et al. 2016; Kainulainen et al. 2017; Bresnahan et al. 2018; Mattern et al. 2018; Ballesteros-Paredes et al. 2020; Arzoumanian et al. 2022). The evolution of filaments is influenced by various mechanisms, such as ISM turbulence, shocked flows, supernovae feedback, or Galactic shear (Chen et al. 2020; Colombo et al. 2021). These filaments can be further classified into various types – for instance, bones or giant molecular filaments (GMFs) – based on their properties (e.g. aspect ratio; Zucker et al. 2018). The variations in filament properties might be a result of the Galactic environment and the scales at which they evolve. For example, recent studies attribute the formation of GMFs to galactic shear and large-scale gas motions, considering these filaments as a subset of giant MCs (Koda et al. 2006; Duarte-Cabral & Dobbs 2016).

The formation and evolution of MCs is still not completely understood (Ballesteros-Paredes et al. 2020). Similarly, there is no unique global property of MCs that can define their ability to form high-mass stars. Modern approaches towards studying MCs include an analysis of their three-dimensional structure (Zucker et al. 2021) and their morphological classification (Yuan et al. 2021) into filamentary and non-filamentary structures. Relating the cloud morphology to its properties and Galactic environment may provide us with a better understanding of MCs and star formation. It could also help us to understand the physical processes behind the structure of clouds (Arzoumanian et al. 2022).

Stellar feedback plays an important role in regulating the star formation in a galaxy and shaping the neighbouring MCs. Stars can impart mechanical energy in the form of turbulence, stellar winds and shocks, and radiative energy, which heats up the gas, leading to the creation and destruction of MCs and their substructures. Feedback drives the mass, energy, momentum, and metal enrichment into the surrounding ISM, inducing new star formation events (Ballesteros-Paredes et al. 2020; Schneider et al. 2020; Beuther et al. 2022). On the other hand, it can also shred the natal molecular gas, limiting further star formation (Geen et al. 2016; Urquhart et al. 2021). Feedback often results in the formation of interstellar bubbles (Deharveng et al. 2010; Jayasinghe et al. 2019), which are some of the most morphologically complex structures in the Galaxy. For instance, the bubble RCW 49 (Schneider et al. 2020; Rodgers et al. 1960) shows an expanding shell that is decoupled from the ambient ISM (Tiwari et al. 2021). Similarly, a fragmented ring of molecular gas is detected around the ionised region surrounding the

bubble RCW 120 (Zavagno et al. 2007; Rodgers et al. 1960), the presence of which is attributed to the collect and collapse mode (Elmegreen & Lada 1977; Zavagno et al. 2006; Zhou et al. 2020; Luisi et al. 2021). These bubbles and shells are prime examples of the ring-like clouds discussed in this work.

In this paper, we try to find a connection between the morphology of MCs and their properties. The paper is organised as follows: Sect. 2 provides a brief overview of the survey data and the cloud catalogue used for the analysis. Section 3 gives a brief description of the methods used for the classification of clouds into different morphologies, and the resulting morphological classes and cloud samples. We have discussed these in detail in Neralwar et al. (2022, hereafter referred to as Paper I). In Sect. 4, we analyse the integrated properties of MCs for the four morphological classes. We also check if the results from the different cloud samples agree with each other for the different morphologies. In Sect. 5, we study Larson's and Heyer's scaling relations, and in Sect. 6 we discuss the results and the physical processes leading to the observed morphologies. We summarise our main findings in Sect. 7.

## 2. Data

Structure, Excitation, and Dynamics of the Inner Galactic Inter-Stellar Medium (SEDIGISM) is a southern hemisphere inner Galaxy survey that probes the moderately dense (≈10³ cm⁻³) ISM. It covers −60° ≤ l ≤ +18° and |b| ≤ 0.5° in several molecular lines, mainly the *J* = 2–1 transitions of ¹³CO and C¹⁸O. The observations for the survey were carried out using the 12m Atacama Pathfinder EXperiment (APEX, Güsten et al. 2006) telescope. A detailed description of the observations, data reduction, and data quality checks is provided in Schuller et al. (2017) and Schuller et al. (2021).

Duarte-Cabral et al. (2021, hereafter called DC21), used the first data release (DR1; Schuller et al. 2021) of the SEDIGISM survey to obtain a catalogue of 10663 clouds using the *scimes*<fnref n="1" /> algorithm (detailed description in Colombo et al. 2015, 2019). The data have a full width at half maximum (FWHM) beam size of 28'' with a 1 σ sensitivity of 0.8–1.0 K at 0.25 km s⁻<fnref n="1" /> spectral resolution. In Paper I, we classify these clouds based on their morphology (see Sect. 3) and thus add more information to the original catalogue<fnref n="2" />.

This cloud catalogue (hereafter called the SEDIGISM cloud catalogue) presents the directly measured properties and the derived properties for the clouds identified with *scimes*. Each cloud is assigned an ID from *scimes* and a name using its Galactic coordinates. The directly measured properties include the cloud's position, velocity, velocity dispersion, and size, and the intensity associated with each cloud. The catalogue also lists near, far, and adopted kinematic distances to the clouds, their reliability and uncertainties, and the presence of star formation tracers. These properties were used to derive other properties such as the aspect ratio, mass, surface density, and the virial parameter. The catalogue also provides the values of the properties obtained after beam deconvolution, which are used in the current work.

## 3. Morphological cloud classification

We study the integrated cloud properties for different morphologies and this demands that the clouds be highly resolved and have

<footnote n="1">¹ https://github.com/Astroua/SCIMES</footnote>

<footnote n="2">² https://sedigism.mpifr-bonn.mpg.de/cgi-bin-seg/SEDIGISM_DATABASE.cgi</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Examples of different cloud morphologies. The images are integrated intensity (moment 0) maps of the <sup>13</sup>CO (2–1) transition. The numbers in the right bottom corner of the images represent the *J*₁ and the *J*₂ moments (*J* plots algorithm), respectively. The colour bars represent the <sup>13</sup>CO integrated intensity in K km s<sup>-1</sup>. The white contours represent the cloud edge. The white circles at the bottom left of the figures represent the telescope beam size.</figcaption></figure>

reliable distance measurements. Thus, we use the clouds belonging to the science sample from DC21 to study their properties. The clouds in the science sample have reliable distance measurements, are well resolved (cloud area >3Ω<sub>beam</sub>), and do not lie on the latitude edges of the survey (edge = 0; in the SEDIGISM catalogue).

Paper I presents the classification of MCs from the SEDIGISM cloud catalogue based on their morphology. It uses two methods for this: *J* plots (Jaffa et al. 2018) and by-eye (visual) classification. It is important to note that *J* plots is an automated algorithm that uses the principal moment of inertia of a structure to determine its morphology. It distinguishes the structures based on their degree of elongation and degree of central concentration, thus classifying them into three types: cores, filaments, and bubbles (rings). In Paper I, we also conducted a visual analysis of the integrated intensity maps of the clouds, which resulted in four morphological classes: (i) ring-like clouds (e.g. Fig. 1a), (ii) elongated clouds (e.g. Fig. 1b), (iii) concentrated clouds (e.g. Fig. 1c), and (iv) clumpy clouds (e.g. Fig. 1d). There are 298 clouds that could not be classified into any of the above classes.

The two classification methods were used to identify two samples from the SEDIGISM cloud catalogue. These are the visually classified (VC) and morphologically reliable (MR) samples. Both of these samples were obtained from the SEDIGISM full sample – a sample of 10663 MCs (detailed description in DC21). The VC sample was obtained by categorising the clouds into the four morphologies using by-eye classification (i.e. the SEDIGISM full sample excluding unclassified clouds), whereas the MR sample contains clouds that have the same morphology under both *J* plot and by-eye classifications. This means that the MR sample includes ring-like, elongated and concentrated clouds that are identified by *J* plots as bubbles, filaments,

and cores, respectively. The clumpy clouds that are recognised as either filaments or cores by the *J* plots algorithm also belong to the MR sample. A detailed description of the classifications and the resulting samples is provided in Sects. 3 and 4.1 of Paper I.

We restrict our MR and VC samples in this paper (Table 1) to only the clouds from the science sample (DC21), and use this sub-sample of clouds for further analysis. These samples are used to understand the connection between cloud properties, their morphologies, and the effects of the Galactic environment. It is worth noting that the quality of classification for the MR sample is better than for the VC sample, as it is classified using both visual and automated classifications. However, this leads to the rejection of a large amount of data and results in a small sample size for some morphologies (e.g. ring-like clouds). In light of this, we use both samples<fnref n="3" /> for our analysis (Table 1).

## 4. Results

In this section, we analyse the connection between cloud morphologies and their other properties, that is to say, cloud mass, surface density, radius, velocity dispersion, length, aspect ratio, and virial parameter. The properties used throughout this paper are the deconvolved properties from DC21. The cloud radius (effective radius and deconvolved radius) is calculated (in DC21) from the exact area of the cloud, assuming a spherical geometry. Therefore, the concept of radius may be inefficient in correctly describing the size of clouds with complex structures, which have a non-spherical geometry. An alternative measure of the major and minor axes of a cloud could be the length and width

<footnote n="3">³ We remind the readers that the VC and MR samples used in this paper contain only the science sample clouds and therefore differ from the VC and MR samples presented in Paper I.</footnote>
<table><caption>**Table 1.** Quantitative description of the clouds.</caption><thead><tr><th>Morphological class</th><th>MR sample</th><th>VC sample</th></tr></thead><tbody><tr><td>Ring-like cloud</td><td>156</td><td>961</td></tr><tr><td>Elongated cloud</td><td>3607</td><td>3676</td></tr><tr><td>Concentrated cloud</td><td>528</td><td>1094</td></tr><tr><td>Clumpy cloud</td><td>729</td><td>733</td></tr><tr><td>Total</td><td>5020</td><td>6464</td></tr></tbody></table>

**Notes.** The table presents the clouds in the visually classified (VC) and morphologically reliable (MR) samples for the different morphologies. We only present the clouds contained in the science sample. For the full SEDIGISM sample, refer to Paper I.

obtained using the geometric medial axis. The longest running central spine through a 2D projected cloud mask is considered to determine the cloud length, whereas the width is twice the average distance from the spine to the cloud edge. The cloud length (length<sub>MA</sub>) and the aspect ratio (AR<sub>MA</sub>) used in this work were obtained using the medial axis technique (DC21). The medial axis length is free of the assumption that clouds have a particular geometry and hence provides a different estimate of the cloud size (especially for elongated structures) compared to the cloud radius. Our implementation of the medial axis analysis does not consider the intensity distribution of the cloud, and therefore we do not trace the density-weighted spine.

The average integrated intensity of a cloud is obtained from the observed flux and is converted into the cloud mass M using the ¹³CO (2–1) to H₂ conversion factor; X<sub>CO</sub> = 1<sup>+1</sup><sub>-0.5</sub> × 10<sup>21</sup> cm<sup>-2</sup> (K km s<sup>-1</sup>)<sup>-1</sup> (Schuller et al. 2017). This mass is then used to obtain the surface density and the virial parameter for the clouds (DC21). The virial parameter is calculated assuming that clouds are spheres (Rosolowsky & Leroy 2006), but this is not consistent with the different cloud morphologies that we see. Moreover, the virial parameter is a crude property calculated using the balance between the kinetic and gravitational energies in the clouds, and might not be an accurate representation of the gravitational state of all the clouds. Nevertheless, we discuss it for the different morphologies as it is a widely used parameter in star formation literature (e.g. DC21).

The violin plots in Figs. 2–5 present distributions of the properties for the different cloud morphologies using the VC and MR samples. The main statistics (i.e. the median, interquartile range, skewness and kurtosis) of cloud properties for the different morphologies for both samples are reported in Table 2. The similar values (especially the medians) between the two samples, in most cases, indicate that they follow each other closely. The two samples show the most similar distributions for elongated and clumpy clouds (p-value = 1, Table A.1). The p-values presented throughout this work were obtained using the two sample Kolmogorov–Smirnov (KS) test. A detailed description of the KS test and p-values is presented in Appendix A. We compare the p-values between the VC and MR samples using the distributions of cloud properties in Table A.1. We find a p-value > 0.01 in ~85% of the cases, which suggests that the VC and MR samples usually show a similar distribution.

The differences between the distributions of the cloud properties for the various morphological classes can be identified visually from the violin plots (Figs. 2–5) and confirmed using the p-values from the KS test (Table A.2). The distributions for the different morphologies differ from each other for each cloud property (typically p-value < 0.01, Table A.2); for instance,

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Distribution of radius (top) and velocity dispersion (bottom) for the VC (blue) and MR (orange) samples. The violin plots present the density of the data at different values, which is smoothed through kernel density estimation. The horizontal lines at the ends and middle of the plots represent the extreme and the median values of the distributions, respectively.</figcaption></figure>

the mass distribution of concentrated and ring-like clouds differ from each other.

Ring-like clouds show the most different properties compared to other morphologies. These clouds show higher average values and different distributions of radius, velocity dispersion, aspect ratio, and length, compared to other morphologies (Figs. 2 and 3; see ‘Median’ values in Table 2). They also have a higher average cloud mass than elongated and clumpy clouds, but the average surface density is comparable to elongated clouds and lower than the other morphologies (Fig. 4). Ring-like clouds also show slightly higher average values for the virial parameter compared to the other morphologies (Fig. 5). However, we again caution readers that the virial parameter might not correctly represent a cloud’s gravitational state for all the morphologies.

Concentrated clouds present a morphology quite different to ring-like clouds. These clouds have a high density at the centre, leading to a near-circular structure, and therefore exhibit smaller (average) lengths and aspect ratios (Fig. 3). They also have a lower average mass than other morphologies, but their small size results in a higher average surface density than other structures (Fig. 4).
<table><caption>**Table 2.** Statistics of integrated properties for various cloud morphologies.</caption><thead><tr><th rowspan="2">Subset</th><th colspan="4">Visually classified sample</th><th colspan="4">Morphologically reliable sample</th></tr><tr><th>Median</th><th>IQR</th><th>Skewness</th><th>Kurtosis</th><th>Median</th><th>IQR</th><th>Skewness</th><th>Kurtosis</th></tr></thead><tbody><tr><td>M [×10<sup>3</sup> M<sub>☉</sub>]</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>1.27</td><td>3.25</td><td>52.3</td><td>3451.1</td><td>1.17</td><td>2.83</td><td>13.2</td><td>294.6</td></tr><tr><td>Ring-like</td><td>2.51</td><td>7.96</td><td>14.1</td><td>283.3</td><td>1.40</td><td>5.89</td><td>2.8</td><td>8.3</td></tr><tr><td>Elongated</td><td>1.09</td><td>2.67</td><td>16.6</td><td>422.3</td><td>1.09</td><td>2.65</td><td>16.9</td><td>437.0</td></tr><tr><td>Concentrated</td><td>1.06</td><td>2.01</td><td>32.1</td><td>1046.9</td><td>1.05</td><td>1.93</td><td>7.2</td><td>71.4</td></tr><tr><td>Clumped</td><td>2.04</td><td>4.71</td><td>5.4</td><td>39.8</td><td>2.04</td><td>4.71</td><td>5.5</td><td>41.2</td></tr><tr><td>Σ [M<sub>☉</sub> pc<sup>-2</sup>]</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>73.42</td><td>42.40</td><td>5.0</td><td>67.1</td><td>72.15</td><td>40.29</td><td>3.7</td><td>28.1</td></tr><tr><td>Ring-like</td><td>71.94</td><td>52.35</td><td>2.7</td><td>10.7</td><td>59.49</td><td>34.30</td><td>2.2</td><td>4.9</td></tr><tr><td>Elongated</td><td>67.04</td><td>32.09</td><td>5.1</td><td>52.6</td><td>67.22</td><td>32.03</td><td>5.1</td><td>53.6</td></tr><tr><td>Concentrated</td><td>92.29</td><td>49.68</td><td>8.2</td><td>144.7</td><td>101.73</td><td>55.33</td><td>1.8</td><td>4.9</td></tr><tr><td>Clumped</td><td>88.99</td><td>47.33</td><td>2.8</td><td>11.6</td><td>88.99</td><td>47.33</td><td>2.8</td><td>11.5</td></tr><tr><td>R [pc]</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>2.32</td><td>2.26</td><td>3.1</td><td>30.7</td><td>2.23</td><td>2.15</td><td>3.7</td><td>43.6</td></tr><tr><td>Ring-like</td><td>3.24</td><td>3.44</td><td>1.4</td><td>3.3</td><td>2.71</td><td>3.25</td><td>1.1</td><td>0.8</td></tr><tr><td>Elongated</td><td>2.24</td><td>2.17</td><td>4.3</td><td>56.2</td><td>2.23</td><td>2.17</td><td>4.3</td><td>57.6</td></tr><tr><td>Concentrated</td><td>1.90</td><td>1.51</td><td>3.2</td><td>24.3</td><td>1.86</td><td>1.40</td><td>1.7</td><td>4.7</td></tr><tr><td>Clumped</td><td>2.67</td><td>2.37</td><td>1.7</td><td>4.1</td><td>2.67</td><td>2.34</td><td>1.7</td><td>4.1</td></tr><tr><td>σ<sub>v</sub>[km s<sup>-1</sup>]</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>0.76</td><td>0.53</td><td>6.9</td><td>158.5</td><td>0.74</td><td>0.49</td><td>2.6</td><td>12.0</td></tr><tr><td>Ring-like</td><td>0.99</td><td>0.89</td><td>2.2</td><td>10.5</td><td>0.84</td><td>0.86</td><td>1.9</td><td>4.3</td></tr><tr><td>Elongated</td><td>0.73</td><td>0.48</td><td>2.6</td><td>11.7</td><td>0.73</td><td>0.48</td><td>2.5</td><td>11.3</td></tr><tr><td>Concentrated</td><td>0.70</td><td>0.39</td><td>19.4</td><td>522.1</td><td>0.72</td><td>0.38</td><td>1.9</td><td>5.5</td></tr><tr><td>Clumped</td><td>0.83</td><td>0.60</td><td>2.5</td><td>11.5</td><td>0.83</td><td>0.60</td><td>2.6</td><td>11.6</td></tr><tr><td>AR<sub>MA</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>4.95</td><td>3.66</td><td>1.6</td><td>4.6</td><td>4.85</td><td>3.46</td><td>1.7</td><td>5.3</td></tr><tr><td>Ring-like</td><td>7.11</td><td>4.72</td><td>1.5</td><td>3.9</td><td>8.18</td><td>5.71</td><td>1.6</td><td>4.4</td></tr><tr><td>Elongated</td><td>5.25</td><td>3.45</td><td>1.5</td><td>3.5</td><td>5.19</td><td>3.41</td><td>1.5</td><td>3.6</td></tr><tr><td>Concentrated</td><td>3.03</td><td>2.04</td><td>2.9</td><td>18.7</td><td>2.75</td><td>1.89</td><td>1.5</td><td>4.6</td></tr><tr><td>Clumped</td><td>4.60</td><td>3.25</td><td>1.3</td><td>2.7</td><td>4.59</td><td>3.21</td><td>1.3</td><td>2.8</td></tr><tr><td>α<sub>vir</sub></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>1.25</td><td>1.31</td><td>8.4</td><td>125.7</td><td>1.25</td><td>1.29</td><td>9.1</td><td>144.1</td></tr><tr><td>Ring-like</td><td>1.50</td><td>1.68</td><td>4.3</td><td>27.1</td><td>1.74</td><td>1.80</td><td>2.4</td><td>6.1</td></tr><tr><td>Elongated</td><td>1.31</td><td>1.34</td><td>8.9</td><td>140.4</td><td>1.31</td><td>1.33</td><td>8.9</td><td>139.5</td></tr><tr><td>Concentrated</td><td>0.98</td><td>0.90</td><td>8.5</td><td>105.7</td><td>0.98</td><td>0.93</td><td>7.8</td><td>94.1</td></tr><tr><td>Clumped</td><td>1.08</td><td>1.18</td><td>9.6</td><td>135.7</td><td>1.08</td><td>1.18</td><td>9.6</td><td>135.0</td></tr><tr><td>Length [pc]</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>All</td><td>7.60</td><td>9.43</td><td>3.6</td><td>37.2</td><td>7.32</td><td>8.77</td><td>4.4</td><td>53.9</td></tr><tr><td>Ring-like</td><td>12.96</td><td>15.21</td><td>1.5</td><td>2.9</td><td>10.84</td><td>15.20</td><td>1.3</td><td>1.5</td></tr><tr><td>Elongated</td><td>7.71</td><td>8.87</td><td>5.0</td><td>66.2</td><td>7.69</td><td>8.82</td><td>5.0</td><td>67.5</td></tr><tr><td>Concentrated</td><td>4.73</td><td>4.87</td><td>3.8</td><td>25.1</td><td>4.38</td><td>4.13</td><td>2.4</td><td>9.9</td></tr><tr><td>Clumped</td><td>8.39</td><td>10.00</td><td>2.1</td><td>6.2</td><td>8.34</td><td>9.95</td><td>2.1</td><td>6.2</td></tr></tbody></table>

**Notes.** The properties are mass (*M*), surface density (Σ), radius (*R*), velocity dispersion (σ<sub>v</sub>), aspect ratio (*AR*<sub>MA</sub>), virial parameter (α<sub>vir</sub>), and length (*length*<sub>MA</sub>). The median is an estimate of the average of a distribution. The interquartile range (IQR) is the difference between the 75th percentile and the 25th percentile of the distribution and is a measure of the variability in the data. Skewness is a measure of symmetry of a distribution, and the positive and negative values signify the data being skewed left and right, respectively. Kurtosis describes the outliers or extremities in the distribution with a high value representing a highly tailed distribution.

Elongated clouds and clumpy clouds have comparable radius, velocity dispersion, length, and aspect ratio<sup>4</sup> distributions

(Figs. 2 and 3), with elongated clouds having more tailed distributions towards both the lower and higher values (see kurtosis

<footnote n="4">⁴ There are 25 elongated and 12 clumpy clouds with an aspect ratio <1.5 in the VC sample. These clouds are low resolution structures that might have been classified incorrectly or have unreliable medial axes</footnote>

<footnote>measurements. We do not exclude these structures as they constitute only a very small percentage of the total sample and do not affect our analysis.</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 3.** Distribution of length (top) and aspect ratio (bottom). The symbols and conventions follow Fig. 2.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Distribution of mass (top) and surface density (bottom). The symbols and conventions follow Fig. 2.</figcaption></figure>

in Table 2). The lower kurtosis of clumpy clouds may be a result of their smaller sample size. The length distributions (Fig. 3) of these two morphologies are similar to those of all clouds, which might be because most of the clouds have an elongated structure. Clumpy clouds have more mass but similar lengths to elongated clouds, resulting in higher surface densities (Fig. 4).

We also checked if the distribution of the clouds in the Galaxy is influenced by environmental factors, namely the Galactocentric radius/distance ($R_{gal}$) and the Galactic height ($z_{gal}$) (see Fig. 6). We did not find any specific environmental trend that separates a particular cloud morphology from others. All cloud morphologies show a Gaussian-like distribution in the $R_{gal}$ histograms (Fig. 6, top). The small sample size at a smaller $R_{gal}$ (Fig. 6; top) is due to the rejection of clouds near the Galactic centre (the science sample in DC21 excludes clouds with unreliable distances), whereas the small sample size at a higher $R_{gal}$ (which for SEDIGISM corresponds to large heliocentric distances) is mostly due to the sensitivity and resolution limitations of the telescope. Unsurprisingly, a large number of MCs are situated at low Galactic heights (Fig. 6, middle), which may be a consequence of the high density of molecular gas in the Galactic

plane. Moreover, all the morphologies seem to be distributed similarly across the Galactic height.

In order to identify any potential observational biases, we also investigated how the different morphologies are distributed as a function of heliocentric distance ($R_d$). The $R_d$ histograms show that ring-like clouds have a higher concentration closer to us (Fig. 6; bottom). This could be due to a projection effect, leading to distortion in the observed structure of the clouds at large distances. Also, the most distant ring-like clouds could be unresolved and thus get classified into other morphologies. This is possibly reflected in the distribution of concentrated clouds, which shows slightly higher number densities at a larger $R_d$. Hence, this is most likely a consequence of observational bias and not a physical phenomenon. The distribution of other morphologies is similar along the heliocentric distance, with a dip in the number of clouds at $R_d \approx 6$ kpc. The clouds in this region lie very close to the tangent velocity and are assigned a higher distance value (i.e. the tangent distance) to avoid distance ambiguity (see Sect. 4.1 in DC21).

We observed that clouds of all morphologies are abundant at all $R_d$ (Fig. 6; bottom). However, the slight trends in the distributions indicate the presence of observational biases. We discuss
<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Distribution of virial parameter for the VC and MR samples. The symbols and conventions follow Fig. 2.</figcaption></figure>

the observational limitations that potentially affect our analysis in Appendix B. The influence of biases is most visible in the directly measured properties, namely mass, radius, and length. For example, smaller and less massive clouds can be detected in our vicinity, but the technical limitation of telescopes minimises their detection at large distances. We also discuss the influence of the environmental factors, that is R<sub>gal</sub> and z<sub>gal</sub> , in Appendix C. We find that smaller clouds are usually seen at low Galactic heights. However, clouds at all Galactic heights show similar aspect ratios. A general observation is that the clouds’ properties follow the global trends (as seen in Figs. 2–5) irrespective of the Galactic environment and the observational biases.

**5. Scaling relations**

Larson's size-linewidth relation (Larson 1981), also known as
Larson's first law, is one of the most studied scaling relations,
which has provided observational constraints on the dynamics of
MCs for the past 40 yr. The relation σ<sub>v</sub> = 1.10 L<sup>0.38</sup> was obtained
by analysing a large number of MCs of different sizes, from sub-
parsec-sized MCs to those measuring a few hundred parsecs.
It has a similar form to Kolmogorov's law for turbulence and
implies that larger clouds have a higher velocity dispersion. It
was further updated for MCs using analysis from Solomon et al.
(1987), as discussed in Colombo et al. (2019).

We present the size-linewidth relation for the various cloud
morphologies in Fig. 7. All cloud morphologies follow the gen-
eral trend of high velocity dispersions for large clouds. However,
it should be noted that the ‘size’ in terms of an equivalent
radius is a better approximation for some clouds than others.
The high values of velocity dispersion for ring-like clouds with
large radii (log(R) > 0.5) might be attributed to the momen-
tum imparted due to stellar feedback. A larger radius means that
more of the gas (shell) has travelled farther away from the cen-
tre, resulting in a higher velocity difference between the centre
and the edge. Although we expect the ring-like clouds to have
a vertical offset given their higher velocity dispersions, this is
not the case and they still lie along Larson’s relation (Fig. 7).
We performed principal component analysis (PCA) for different
cloud morphologies on the Larson’s plot (Fig. 7) to obtain the

<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Histogram (normalised) for distribution of clouds with respect to Galactocentric radius (top) Galactic height (z<sub>gal</sub>) (middle), and heliocentric distance (bottom). The histograms represent the ratio of clouds in each bin to the total number of clouds in all the bins. The bin width is 1 kpc for the heliocentric distance and Galactocentric radius plots, and 1 pc for the z<sub>gal</sub> histogram. Blue represents the VC sample and orange represents the MR sample. The cloud morphologies *from top to bottom* represent ring-like clouds (Ring), elongated clouds (Elon), concentrated clouds (Conc), and clumpy clouds (Clumpy).</figcaption></figure>

confidence ellipses (1-sigma) and compared the cloud distribu-
tions with the Larson's relation. The slopes of the confidence
ellipses for the clouds (Table 3) may indicate that elongated
clouds and clumpy clouds (slope = 0.48) follow Larson's first law
(slope = 0.50) more closely than other morphologies. However,
the large scatter in the data (Fig. 7) prevents us from confirming
this.

Larson's second relation correlates the velocity dispersion
to cloud mass, implying that clouds are approximately in virial
equilibrium, while the third relation shows the anti-correlation
between cloud mean density and size. These relations, together
with the analysis of Solomon et al. (1987), assert that clouds have
a constant surface mass density. Heyer et al. (2009) reanalysed
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** Size-linewidth relation ($\sigma_v$ versus R) using the MR sample. The dashed line represents Larson's first relation (Larson 1981; Solomon et al. 1987). The different colours represent four cloud classes (Sect. 3); ring-like clouds are represented in blue, elongated clouds in red, concentrated clouds in green, and clumpy clouds in magenta. The contours include 68% of the data (1-sigma level) and the small circles represent centroids (means) of the distributions. The expected slope for the relation is 0.5 (dashed black line) (Larson 1981).</figcaption></figure>

<table><caption>**Table 3.** Slopes and their standard deviations (SD) recovered from PCA on Larson's scaling relation ($\sigma_v \propto R^{\text{slope}}$) for the different cloud morphologies.</caption><thead><tr><th>Cloud</th><th>Ring</th><th>Elongated</th><th>Concentrated</th><th>Clumpy</th></tr></thead><tbody><tr><td>Slope</td><td>0.63</td><td>0.48</td><td>0.41</td><td>0.48</td></tr><tr><td>SD</td><td>0.01</td><td>&lt;0.01</td><td>0.01</td><td>0.01</td></tr></tbody></table>

**Notes.** The values were obtained using the bootstrap technique, where different samples were generated by varying the values of x and y parameters of the scaling relation within their errors. We calculated the slope for each sample and report the mean and standard deviation of these slopes in the table. The results were obtained using 10 000 random samples generated using random seed 0. The expected value of the slope is 0.5 (Larson 1981).

the clouds in the first quadrant of the Galaxy and found that these clouds do not have a constant surface mass density and they are not necessarily gravitationally bound. We show Heyer's relation ($\sigma_v^2/R \propto \Sigma$) for different cloud morphologies in Fig. 8 and report the slopes of the corresponding 1-sigma PCA ellipses in Table 4. The different cloud morphologies show a trend with respect to the virial parameter, with ring-like clouds present towards a region of higher virial parameter. This is understood to mean that ring-like clouds are effectively expanding and thus not gravitationally bound. Towards a higher surface density and a lower velocity dispersion (i.e. a lower virial parameter), the regions are correspondingly occupied by elongated, clumpy, and

<table><caption>**Table 4.** Slopes recovered from PCA on the Heyer scaling relation ($\sigma_v^2/R \propto \Sigma^{\text{slope}}$).</caption><thead><tr><th>Cloud</th><th>Ring</th><th>Elongated</th><th>Concentrated</th><th>Clumpy</th></tr></thead><tbody><tr><td>Slope</td><td>3.17</td><td>4.50</td><td>3.46</td><td>3.87</td></tr><tr><td>SD</td><td>0.21</td><td>0.12</td><td>0.20</td><td>0.21</td></tr></tbody></table>

**Notes.** The terminology follows Table 3. The expected value of the slope is 1 (Heyer et al. 2009).

then concentrated clouds. This is in agreement with our previous analysis of concentrated clouds being dense, gravitationally bound compact objects. We also observe that none of the cloud distributions follow $\alpha_{vir} = 1$ (Heyer et al. 2009), but the bulk of clouds could be virialised if they were under a relatively mild external pressure or were contracting at some fraction of free fall (Fig. 8).

## 6. Discussion

The differences in properties of different morphologies can be associated with various physical processes that play a role in the formation and evolution of clouds. The large size (as measured by the length and radius in Figs. 2 and 3) and aspect ratio (Fig. 3) for ring-like clouds can be attributed to their formation process. These clouds could be associated with infrared bubbles, which form due to stellar winds and radiation feedback. We find in Paper I that 605 ring-like clouds show association with infrared
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Scaling relation between $\sigma_p^2/R$ and gas surface density ($\Sigma$) using the MR sample. The black solid lines and dashed lines (from bottom to top) represent $\alpha_{vir} = 1$ in the absence of external pressure and at an external pressure $P_{ext} = 1, 10$, and $100 M_{\odot} \text{pc}^{-3} \text{km}^2 \text{s}^{-2}$, respectively. The symbols and conventions follow Fig. 7.</figcaption></figure>

bubbles from the Milky Way Project (MWP). Stellar feedback could impart high velocity and turbulence to the shells, which then travel large distances to form large rings with high velocity dispersion (Snell et al. 1980; Liu et al. 2019). Ring-like clouds have the highest mass (on average). However, the large sizes of these clouds result in them having lower surface densities than other structures. Overall, ring-like clouds show the most extreme properties compared to the global distribution of clouds, as well as other morphologies.

The high mass and surface density (Fig. 4) of the clumpy clouds can be explained in two possible ways. The first hypothesis is that these are multiple overlapping clouds, considered to be a single cloud by *scimes*. The other hypothesis is that these clouds evolve similar to globular filaments (Schneider & Elmegreen 1979), having multiple dense regions in the filamentary structure. The second hypothesis is supported by higher star formation efficiencies (refer Paper I) for these clouds when compared to elongated clouds. Clumpy clouds have length distributions similar to elongated clouds, which could be due to a common formation or evolution mechanism, leading to filamentary structures of comparable sizes. It could also be an artefact of the cloud selection criteria used by the *scimes* algorithm, which is oriented towards selecting clouds with similar sizes and could thus present a single MC as smaller individual clouds, or combine small MCs and present them as a single cloud.

A major observation across all the data samples and throughout the different Galactic environments is that most of the MCs are elongated (Table 1; histograms in Figs. C.1, C.2, C.9, C.10, B.1 and B.2). The average properties of elongated clouds lie

between those of ring-like clouds and concentrated clouds. The large sample size for elongated clouds might be the reason for their values being non-extreme. We also find that the Galactic environment does not significantly influence the trends in the integrated properties for most of the morphologies (see the ridge plots in Appendix B).

## 7. Summary

In this paper, we study the connection between the morphology of clouds and their integrated properties. We also analyse how the different cloud morphologies are distributed in the Galaxy. Our analysis is based on the catalogue of MCs from the SEDIGISM survey (presented in DC21; along with the integrated properties). Clouds were classified into four types based on their morphology (in Paper I). Our main findings are as follows:

- Elongated clouds are the most abundant type of clouds in the Galaxy, with their average properties lying between ring-like and concentrated clouds;

- Concentrated clouds have the smallest median sizes and highest median surface densities of all cloud categories;

- Ring-like clouds show the most discrepant properties compared to the global population of clouds, with higher average mass, radius, length, aspect ratio, and velocity dispersion than other morphologies;

- The different morphology classes are not distributed differently in terms of Galactocentric distance ($R_{gal}$) and Galactic height ($z_{gal}$);
- The distributions of clouds do not show significant differences on Larson's size-linewidth plot for the different morphologies;

- Heyer's scaling relation shows the evolution of morphologies with the virial parameter, with ring-like, elongated, clumpy, and concentrated clouds lying in order of decreasing virial parameter.

In conclusion, we find some connections between the integrated properties of the MCs and their morphology. The ring-like clouds show an extreme behaviour compared to other cloud morphologies, which may be a consequence of their nature, in other words, their formation due to stellar feedback. Moreover, the global trends in the properties for different morphologies are followed even in separate smaller regions of the Galaxy (i.e. throughout the distance bins (Appendix C)), irrespective of the observational biases (Appendix B). Even though we did not find any particular trend in the distribution of different types of clouds as a function of Galactocentric distances or Galactic height, we have not looked into other environments within the Galaxy, possibly more critical to the evolution of clouds, such as spiral arms or the Galactic centre. The categorisation of clouds into different morphologies could thus give way to various potential follow-up studies to help understand the role of the Galaxy in shaping the formation and evolution of MCs.

**Acknowledgement.** The authors thank the referee for their valuable comments on the draft, which helped improve the quality of the paper. This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX) under programmes 092.F-9315 and 193.C-0584. APEX is a collaboration among the Max-Planck-Institut für Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory. This publication uses data generated via the Zooniverse.org platform, development of which is funded by generous support, including a Global Impact Award from Google, and by a grant from the Alfred P. Sloan Foundation. A part of this work is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. K.R.N. thanks Henrik Beuther for the valuable suggestions and feedback on the paper. K.R.N. would also like to thank Antonio Hernandez for his careful reading of the manuscript. DC acknowledges support from the German Deutsche Forschungsgemeinschaft, DFG project number SFB956A. A.D.C. acknowledges the support from the Royal Society University Research Fellowship (URF/R1/191609). H.B. acknowledges support from the European Research Council under the Horizon 2020 Framework Programme via the ERC Consolidator Grant CSF-648505. H.B. also acknowledges support from the Deutsche Forschungsgemeinschaft in the Collaborative Research Center (SFB 881) "The Milky Way System" (sub-project B1). C.L.D. acknowledges funding from the European Research Council for the Horizon 2020 ERC consolidator grant project ICYBOB, grant number 818940.

## References

Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6

Arzoumanian, D., Russeil, D., Zavagno, A., et al. 2022, A&A, 660, A56

Ballesteros-Paredes, J., Hartmann, L. W., Vázquez-Semadeni, E., Heitsch, F., & Zamora-Avilés, M. A. 2011, MNRAS, 411, 65

Ballesteros-Paredes, J., André, P., Hennebelle, P., et al. 2020, Space Sci. Rev., 216, 76

Barnes, P., Li, D., Telesco, C., et al. 2015, MNRAS, 453, 2622

Barnes, P. J., Hernandez, A. K., Muller, E., & Pitts, R. L. 2018, ApJ, 866, 19

Beuther, H., Schneider, N., Simon, R., et al. 2022, A&A, 659, A77

Blitz, L., & Stark, A. A. 1986, ApJ, 300, L89

Bresnahan, D., Ward-Thompson, D., Kirk, J. M., et al. 2018, A&A, 615, A125

Chapman, N. L., Goldsmith, P. F., Pineda, J. L., et al. 2011, ApJ, 741, 21

Chen, C.-Y., Mundy, L. G., Ostriker, E. C., Storm, S., & Dhabal, A. 2020, MNRAS, 494, 3675

Colombo, D., Schinnerer, E., Hughes, A., et al. 2013, in American Astronomical Society Meeting Abstracts, 221, 349.16

Colombo, D., Rosolowsky, E., Ginsburg, A., Duarte-Cabral, A., & Hughes, A. 2015, MNRAS, 454, 2067

Colombo, D., Rosolowsky, E., Duarte-Cabral, A., et al. 2019, MNRAS, 483, 4291

Colombo, D., König, C., Urquhart, J. S., et al. 2021, A&A, 655, L2

Crutcher, R. M. 2012, ARA&A, 50, 29

Crutcher, R. M., Wandelt, B., Heiles, C., Falgarone, E., & Troland, T. H. 2010, ApJ, 725, 466

Deharveng, L., Schuller, F., Anderson, L. D., et al. 2010, A&A, 523, A6

Dobbs, C. L., Krumholz, M. R., Ballesteros-Paredes, J., et al. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 3

Duarte-Cabral, A., & Dobbs, C. L. 2016, MNRAS, 458, 3667

Duarte-Cabral, A., Colombo, D., Urquhart, J. S., et al. 2021, MNRAS, 500, 3027

Elmegreen, B. G. & Lada, C. J. 1977, ApJ, 214, 725

Field, G. B., Blackman, E. G., & Keto, E. R. 2011, MNRAS, 416, 710

Fukui, Y., Kawamura, A., Minamidani, T., et al. 2008, ApJS, 178, 56

Geen, S., Hennebelle, P., Tremblin, P., & Rosdahl, J. 2016, MNRAS, 463, 3129

Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13

Hacar, A., Tafalla, M., Kauffmann, J., & Kovács, A. 2013, A&A, 554, A55

Heyer, M., & Dame, T. M. 2015, ARA&A, 53, 583

Heyer, M., Krawczyk, C., Duval, J., & Jackson, J. M. 2009, ApJ, 699, 1092

Jaffa, S. E., Whitworth, A. P., Clarke, S. D., & Howard, A. D. P. 2018, MNRAS, 477, 1940

Jayasinghe, T., Dixon, D., Povich, M. S., et al. 2019, MNRAS, 488, 1141

Kainulainen, J., Stutz, A. M., Stanke, T., et al. 2017, A&A, 600, A141

Koda, J., Sawada, T., Hasegawa, T., & Scoville, N. Z. 2006, ApJ, 638, 191

Larson, R. B. 1981, MNRAS, 194, 809

Li, G.-X., Urquhart, J. S., Leurini, S., et al. 2016, A&A, 591, A5

Liu, M., Li, D., Krčo, M., et al. 2019, ApJ, 885, 124

Luisi, M., Anderson, L. D., Schneider, N., et al. 2021, Sci. Adv., 7, eabe9511

Mattern, M., Kauffmann, J., Csengeri, T., et al. 2018, A&A, 619, A166

Miville-Deschênes, M.-A., Murray, N., & Lee, E. J. 2017, ApJ, 834, 57

Neralwar, K. R., Colombo, D., Duarte-Cabral, A., et al. 2022, A&A, 663, A56 (Paper I)

Olmi, L., Cunningham, M., Elia, D., & Jones, P. 2016, A&A, 594, A58

Polychroni, D., Schisano, E., Elia, D., et al. 2013, ApJ, 777, L33

Rebolledo, D., Wong, T., Leroy, A., Koda, J., & Donovan Meyer, J. 2012, ApJ, 757, 155

Rodgers, A. W., Campbell, C. T., & Whiteoak, J. B. 1960, MNRAS, 121, 103

Roman-Duval, J., Jackson, J. M., Heyer, M., Rathborne, J., & Simon, R. 2010, ApJ, 723, 492

Rosolowsky, E., & Leroy, A. 2006, Publ. ASP, 118, 590

Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, ApJ, 679, 1338

Schneider, S., & Elmegreen, B. G. 1979, ApJS, 41, 87

Schneider, N., Simon, R., Guevara, C., et al. 2020, Publ. ASP, 132, 104301

Schuller, F., Csengeri, T., Urquhart, J. S., et al. 2017, A&A, 601, A124

Schuller, F., Urquhart, J. S., Csengeri, T., et al. 2021, MNRAS, 500, 3064

Snell, R. L., Loren, R. B., & Plambeck, R. L. 1980, ApJ, 239, L17

Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730

Sun, J., Leroy, A. K., Schruba, A., et al. 2018, ApJ, 860, 172

Tiwari, M., Karim, R., Pound, M. W., et al. 2021, ApJ, 914, 117

Urquhart, J. S., Figura, C., Cross, J. R., et al. 2021, MNRAS, 500, 3050

Vázquez-Semadeni, E., Palau, A., Ballesteros-Paredes, J., Gómez, G. C., & Zamora-Avilés, M. 2019, MNRAS, 490, 3061

Yuan, L., Yang, J., Du, F., et al. 2021, ApJS, 257, 51

Zavagno, A., Deharveng, L., Comerón, F., et al. 2006, A&A, 446, 171

Zavagno, A., Pomarès, M., Deharveng, L., et al. 2007, A&A, 472, 835

Zhou, J., Zhou, D., Esimbek, J., et al. 2020, ApJ, 897, 74

Zucker, C., Battersby, C., & Goodman, A. 2018, ApJ, 864, 153

Zucker, C., Goodman, A., Alves, J., et al. 2021, ApJ, 919, 35
## Appendix A: Kolmogorov?Smirnov test and p-value

The two-sample KS test is a non-parametric goodness-of-fit test that compares the cumulative distribution of two datasets by calculating the p-value and the maximum difference between them. We use the null hypothesis that the two distributions are obtained from the same distribution. Hence, a low p-value (typically below 0.01) leads to the rejection of the null hypothesis, and it is an indication that the two samples belong to different distributions. All the p-values throughout this paper were obtained from the KS test. We used the ks_2samp module from the python scipy package to implement the KS test on our distributions.

We used the KS test for two major tasks. The first was to confirm whether the distribution of the properties for a particular morphology is the same for the VC and MR samples. We confirm that the VC and MR samples are consistent with each other in Table A.1. The second task was to check if different morphologies show different distribution for a given property. The differences in the properties for the different morphologies can be visually distinguished in the violin plots (Fig. 3–5) and ridge plots (appendix B and C), and they are statistically confirmed in Table A.2.

<table><caption>Table A.1: P-values for the two-sample KS test conducted on the VC and MR samples for each property and morphology.</caption><thead><tr><th>Cloud</th><th>M</th><th>Σ</th><th>R</th><th>σ<sub>v</sub></th><th>AR</th><th>α<sub>vir</sub></th><th>Length</th></tr></thead><tbody><tr><td>Total</td><td>7.7e-02</td><td>1.1e-01</td><td>1.1e-01</td><td>7.5e-02</td><td>2.7e-01</td><td>1.0e+00</td><td>8.8e-02</td></tr><tr><td>Ring-like</td><td>4.5e-03</td><td>3.8e-05</td><td>3.5e-02</td><td>4.4e-02</td><td>6.8e-03</td><td>1.5e-01</td><td>3.9e-02</td></tr><tr><td>Elongated</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td></tr><tr><td>Concentrated</td><td>9.0e-01</td><td>1.8e-03</td><td>2.9e-01</td><td>6.3e-01</td><td>3.5e-03</td><td>1.0e+00</td><td>4.2e-02</td></tr><tr><td>Clumpy</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td><td>1.0e+00</td></tr></tbody></table>

**Notes.** The rows represent the cloud morphologies whereas the columns represent the cloud properties. The symbols follow Table 2.

<table><caption>Table A.2: P-values for the two-sample KS test conducted on the different morphologies for each cloud property (VC sample).</caption><thead><tr><th>Cloud</th><th>M</th><th>Σ</th><th>R</th><th>σ<sub>v</sub></th><th>AR</th><th>α<sub>vir</sub></th><th>Length</th></tr></thead><tbody><tr><td>Tot-Ring</td><td>5.9e-15</td><td>3.3e-02</td><td>5.9e-15</td><td>5.9e-15</td><td>5.9e-15</td><td>5.3e-10</td><td>5.9e-15</td></tr><tr><td>Tot-Elon</td><td>2.4e-04</td><td>2.9e-15</td><td>1.1e-01</td><td>1.2e-03</td><td>8.5e-09</td><td>1.2e-03</td><td>4.6e-01</td></tr><tr><td>Tot-Conc</td><td>1.6e-07</td><td>1.8e-15</td><td>1.8e-15</td><td>4.1e-10</td><td>1.8e-15</td><td>1.8e-15</td><td>1.8e-15</td></tr><tr><td>Tot-Clumpy</td><td>5.4e-10</td><td>8.8e-29</td><td>5.7e-06</td><td>3.3e-05</td><td>2.1e-02</td><td>3.3e-05</td><td>1.7e-02</td></tr><tr><td>Ring-Elon</td><td>2.1e-15</td><td>1.5e-13</td><td>2.1e-15</td><td>2.1e-15</td><td>2.1e-15</td><td>8.6e-06</td><td>2.1e-15</td></tr><tr><td>Ring-Conc</td><td>6.9e-32</td><td>1.3e-30</td><td>2.7e-54</td><td>1.1e-43</td><td>1.8e-191</td><td>2.7e-28</td><td>5.0e-103</td></tr><tr><td>Ring-Clumpy</td><td>1.4e-03</td><td>1.7e-15</td><td>2.1e-07</td><td>8.2e-07</td><td>1.7e-15</td><td>1.1e-14</td><td>1.7e-15</td></tr><tr><td>Elon-Conc</td><td>5.6e-03</td><td>9.6e-90</td><td>3.5e-12</td><td>1.3e-04</td><td>1.7e-141</td><td>1.6e-24</td><td>6.8e-50</td></tr><tr><td>Elon-Clumpy</td><td>3.0e-16</td><td>4.1e-59</td><td>3.7e-08</td><td>3.4e-09</td><td>5.8e-07</td><td>6.5e-09</td><td>1.9e-02</td></tr><tr><td>Conc-Clumpy</td><td>1.2e-15</td><td>1.4e-01</td><td>1.2e-15</td><td>5.3e-15</td><td>1.2e-15</td><td>9.8e-03</td><td>1.2e-15</td></tr></tbody></table>

**Notes.** The rows represent the two morphologies for which the KS test was conducted. Tot, Ring, Elon, Conc, and Clumpy refer to the total, ring-like, elongated, concentrated, and clumpy clouds, respectively. The columns represent the property for which the p-value was obtained. The symbols follow Table 2.
# Appendix B: Testing the observational biases affecting cloud properties and morphological trends

The sensitivity and resolution limitations of the telescope may result in the incorrect classification of cloud morphology at large distances ($R_d$). We discuss the effects of these technical limitations in detail in this section using the ridge plots for heliocentric distance ($R_d$). In general, ridge plots are a visualisation tool that show the distribution of a numerical value for various groups. In our case, they can be considered violin plots for cloud properties in separate distance bins. The clouds are separated into nine distance bins with a 2 kpc width each, covering 0–18 kpc, based on their proximity to the Sun. Fig. B.1 shows the number of structures in each distance bin from 0–24 kpc.

The cloud properties that are most affected by the technical limitations are mass and size (radius and length). They show visible changes in the median as well as the distribution range across the distance bins. Fig. B.2 shows that the lower limit on mass increases with $R_d$. At a smaller $R_d$, low-mass clouds can be detected due to the high resolution and sensitivity of the telescope, but as $R_d$ increases, the low-mass clouds are perceived as background noise. Similarly, we also find that the lower limit on radius (Fig. B.4) and length (Fig. B.8) increases with distance, following the same principle as mass distributions.

The surface density (Fig. B.3), aspect ratio (Fig. B.6), and velocity dispersion are not strongly affected by the technical limitations. The distributions for the virial parameter (Fig. B.7) show the most overlap compared to those of other properties. The decrease in the virial parameter with an increasing $R_d$ can be attributed to observational biases at play (discussed in appendix C of DC21).

We study the cloud properties in a two-dimensional position-position space, and it could lead to the projection effects affecting our analysis. For example, our cloud length and width do not represent the original dimensions of the MCs, but rather a projected length and width. Similarly, the cloud morphology is also influenced by the positioning of clouds with respect to the telescope's line of sight. However, given the large size of our data sample, the general trends in the results are unlikely to be affected, and thus we ignore the projection effects in this work.

It is interesting to note that, regardless of the changes in the properties with distance, the properties show similar trends across all the distance bins for the different morphologies. This validates our conclusions on the trends found for the full samples (see Sec. 4).
<figure><img src="image_11.png" /><figcaption>Fig. B.1: Mass ridge plot for the VC sample with $R_d$ bins.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>Fig. B.2: Mass ridge plot for the MR sample with $R_d$ bins.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Fig. B.3: Surface density ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Fig. B.4: Radius ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Fig. B.5: Velocity dispersion ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>Fig. B.6: Aspect ratio ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>Fig. B.7: Virial parameter ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Fig. B.8: Length ridge plots with $R_d$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
# Appendix C: Dependence of cloud properties on the Galactic environment

The properties of MCs might show changes within different Galactic environments and based on their position in the Galaxy. We studied changes in the integrated properties of MCs based on their location in the Galaxy for the VC and MR samples using ridge plots. The data in each bin were plotted separately for each morphology in order to understand the environmental influence on each of the cloud properties for the different morphologies. The distributions of elongated clouds and clumpy clouds for the two samples are nearly identical as only a few of these clouds from the VC sample are not a part of the MR sample. The clouds are separated into distance bins based on two distance<fnref n="5" /> types. The first is the Galactocentric radius ($R_{gal}$), that is to say, the distance of a cloud from the centre of the Milky Way. The second is the Galactic height ($z_{gal}$), in other words, the vertical height of a cloud from the Galactic plane. The cloud properties, namely mass, surface density, radius, velocity dispersion, aspect ratio, virial parameter, and length, were studied using ridge plots. The low sample size of morphologies in some bins led to sporadic peaks in the plots and also prohibited an in-depth study.

## Appendix C.1: Galactocentric radius ($R_{gal}$) bins

This section presents the ridge plots with Galactocentric radius ($R_{gal}$) bins. The data are divided into eight bins of width 1 kpc from 2–10 kpc<fnref n="6" />. The central part of the Galaxy (< 1 kpc) is excluded because the local gas motions influence the overall properties, causing a change in the dynamics of the region. The decrease in sample size at a large distance from the Galactic centre sets an upper limit of 10 kpc for $R_{gal}$ bins (Fig. C.1 and C.2).

The mass distribution for the two samples (Fig. C.1 and C.2) show the highest values for ring-like clouds and clumpy clouds. The surface density distributions (Fig. C.3) follow the global trends (Sec. 4), showing high values for concentrated and clumpy clouds. The properties pertaining to the size of the clouds show similar features across most of the $R_{gal}$ bins. These are the radius (Fig. C.4) and length (Fig. C.8), which show the highest values for ring-like clouds. The aspect ratio distribution (Fig. C.6) also follows similar trends. These high values can be explained as a consequence of stellar feedback expanding the bubbles (see Sec. 6).

## Appendix C.2: Galactic height ($z_{gal}$) bins

In this section, we study the clouds by distributing them into $z_{gal}$ bins. The highest $z_{gal}$ value for the clouds is ≈ 230 pc, however, < 1 % clouds lie above 160 pc. Thus, we binned the clouds into 8 bins of width 20 pc covering 0–160 pc. Clouds near the Galactic plane have a larger span of sizes (both radius and length), as seen in Figs. C.12 and C.16. Smaller clouds at higher Galactic heights might be obscured due to the foreground gas in the Galactic plane. The sensitivity and resolution of the telescope could also prevent the detection of smaller clouds at large distances. Clouds at all Galactic heights exhibit almost similar aspect ratio distributions (Fig. C.14), demonstrating that they show similar average shapes for a given morphology. We also see that less massive clouds are mostly located near the Galactic plane (Figs. C.9 and C.10). The virial parameter distribution has a mild trend, suggesting a lower average value with increasing Galactic height (Fig. C.15). To summarise, we detect a greater number of smaller and less massive clouds in the Galactic plane than at large Galactic heights, but this could be due to the foreground molecular gas, telescope limitations, or both.

<footnote n="5">⁵ The distances are provided as $R_{gal}$ and $z_{gal_kpc}$ in the cloud catalogue of DC21.</footnote>

<footnote n="6">⁶ The MR sample has ring-like clouds up to $R_{gal}$ = 8 kpc only.</footnote>
<figure><img src="image_19.png" /><figcaption>Fig. C.1: Mass ridge plot for the VC sample with $R_{\text{gal}}$ bins. The vertical dashed lines represent the medians of the distributions.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>Fig. C.2: Mass ridge plot for the MR sample with $R_{\text{gal}}$ bins. The vertical dashed lines represent the medians of the distributions.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>Fig. C.3: Surface density ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>Fig. C.4: Radius ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_23.png" /><figcaption>Fig. C.5: Velocity Dispersion ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>Fig. C.6: Aspect ratio ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_25.png" /><figcaption>Fig. C.7: Virial Parameter ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>Fig. C.8: Length ridge plots with $R_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_27.png" /><figcaption>Fig. C.9: Mass ridge plot for the VC sample with $z_{\text{gal}}$ bins.</figcaption></figure>
<figure><img src="image_28.png" /><figcaption>Fig. C.10: Mass ridge plot for the MR sample with $z_{\text{gal}}$ bins.</figcaption></figure>

<figure><img src="image_29.png" /><figcaption>Fig. C.11: Surface density ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_30.png" /><figcaption>Fig. C.12: Radius ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>Fig. C.13: Velocity Dispersion ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_32.png" /><figcaption>Fig. C.14: Aspect ratio ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>

<figure><img src="image_33.png" /><figcaption>Fig. C.15: Virial parameter ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>
<figure><img src="image_34.png" /><figcaption>Fig. C.16: Length ridge plots with $z_{\text{gal}}$ bins. Left (solid): VC sample. Right (dashed): MR sample.</figcaption></figure>