# Asymmetries in random motions of neutral hydrogen gas in spiral galaxies

P. Adamczyk, P. Amram, L. Chemin, B. Epinat, J. Braine, F. Combes, W. J.

G. de Blok

► To cite this version:

P. Adamczyk, P. Amram, L. Chemin, B. Epinat, J. Braine, et al. Asymmetries in random motions of neutral hydrogen gas in spiral galaxies. Astronomy and Astrophysics - A&A, 2023, 678, 10.1051/0004-6361/202346750. insu-04505559

HAL Id: insu-04505559

https://insu.hal.science/insu-04505559

Submitted on 15 Mar 2024

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Asymmetries in random motions of neutral hydrogen gas in spiral galaxies

P. Adamczyk, P. Amram, L. Chemin, B. Epinat, J. Braine, F. Combes, W. J.

G. de Blok

► To cite this version:

P. Adamczyk, P. Amram, L. Chemin, B. Epinat, J. Braine, et al. Asymmetries in random motions of neutral hydrogen gas in spiral galaxies. Astronomy & Astrophysics - A&A, 2023, 678, 10.1051/0004-6361/202346750. insu-04505559

HAL Id: insu-04505559

https://insu.hal.science/insu-04505559v1

Submitted on 15 Mar 2024

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Asymmetries in random motions of neutral hydrogen gas in spiral galaxies

P. Adamczyk<fnref n="1" />,<fnref n="2" />, P. Amram<fnref n="2" />, L. Chemin<fnref n="3" />, B. Epinat<fnref n="2" />,<fnref n="4" />, J. Braine<fnref n="5" />, F. Combes<fnref n="6" />, and W. J. G. de Blok<fnref n="7" />,<fnref n="8" />,<fnref n="9" />

Received 26 April 2023 / Accepted 8 June 2023

## ABSTRACT

**Context.** The velocity dispersion ellipsoid of gas in galactic discs is usually assumed to be isotropic. Under this approximation, no projection effect occurs in the random motions of gas, as traced by the line-of-sight velocity dispersion. However, it has been recently shown that random motions of the neutral hydrogen gas of the Triangulum galaxy (M 33) exhibit a bisymmetric perturbation which is aligned with the minor axis of the galaxy, suggesting a projection effect.

**Aims.** To investigate if perturbations in the velocity dispersion of nearby discs are comparable to those of M 33, the sample is extended to 32 galaxies from The H I Nearby Galaxy Survey (THINGS) and the Westerbork H I Survey of Spiral and Irregular Galaxies (WHISP).

**Methods.** We studied velocity asymmetries in the disc planes by performing Fourier transforms of high-resolution H I velocity dispersion maps corrected for beam-smearing effects, and we measured the amplitudes and phase angles of the Fourier harmonics.

**Results.** In all velocity dispersion maps, we find strong perturbations of first, second, and fourth orders. The strongest asymmetry is the bisymmetry, which is predominantly associated with the presence of spiral arms. The first order asymmetry is generally orientated close to the disc major axis, and the second and fourth order asymmetries are preferentially orientated along intermediate directions between the major and minor axes of the discs. These results are evidence that strong projection effects shape the H I velocity dispersion maps. The most likely source of systematic orientations is the anisotropy of velocities, through the projection of streaming motions that are stronger along one of the planar directions in the discs. Moreover, systematic phase angles of asymmetries in the H I velocity dispersion could arise from tilted velocity ellipsoids, that is when the velocities are correlated. We expect a larger incidence of correlation between the radial and tangential velocities of H I gas with $|ρ_{Rθ}| \sim 0.6$, which could be tested against the kinematics of the youngest stellar populations of the Milky Way.

**Conclusions.** H I velocity dispersions cannot be considered devoid of projection effects. The systematic orientations of asymmetries can be explained by the projection of unresolved streaming motions mainly arising from spiral arms. Our methodology is a powerful tool to constrain the dominant direction of streaming motions and thus the shape of the velocity ellipsoid of H I gas, which is de facto anisotropic at the angular scales probed by the observations. The next step is to study the shape of the velocity ellipsoids of molecular and ionised gas and their link with galaxy mass and/or morphology, in addition to extending the sample size.

**Key words.** galaxies: fundamental parameters – galaxies: kinematics and dynamics – galaxies: spiral – galaxies: structure – galaxies: ISM

## 1. Introduction

In discs of galaxies, large-scale perturbations such as bars, spiral arms, warps, or lopsidedness make the orbits non-circular. They generate streaming motions, which are asymmetric radial and tangential velocities in the stellar and gaseous media in galactic discs (e.g. Visser 1980). This is beautifully illustrated from an observational viewpoint via the kinematics of millions of individual stars in the Large Magellanic Cloud (LMC), for which the Gaia Collaboration used astrometry to deduce the velocity fields of both the radial and tangential components of

a disc galaxy (Gaia Collaboration 2016, 2018a,b, 2021a), and from which significant asymmetries are observed in the velocity maps, as caused by the bar and spiral arms of the LMC (Gaia Collaboration 2021b).

The study of kinematic asymmetries has become relevant in galactic dynamics, as it impacts our knowledge of the structure of dark matter halos and the comparison with simulations made in a cosmological context. Interstellar gas is thought to follow circular orbits as it undergoes dissipative collisions. If asymmetries are important, then there is no guarantee that rotation curves perfectly trace the mass distribution of galaxies.

<footnote n="1">¹ Centro de Astronomia, Universidad de Antofagasta, Avda. U. de Antofagasta, 02800 Antofagasta, Chile
e-mail: p.adamczyk@free.fr</footnote>

<footnote n="2">² Aix-Marseille Univ., CNRS, CNES, LAM, 38 Rue Frédéric Joliot Curie, 13338 Marseille, France
e-mail: philippe.amram@lam.fr</footnote>

<footnote n="3">³ Instituto de Astrofísica, Universidad Andes Bello, Fernandez Concha 700, Las Condes, Santiago RM, Chile
e-mail: astro.chemin@gmail.com</footnote>

<footnote n="4">⁴ Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Highway, Kamuela, HI 96743, USA</footnote>

<footnote n="5">⁵ Laboratoire d'Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="6">⁶ Observatoire de Paris, LERMA, Collège de France, CNRS, PSL University, Sorbonne University, 75014 Paris, France</footnote>

<footnote n="7">⁷ ASTRON, Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands</footnote>

<footnote n="8">⁸ Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands</footnote>

<footnote n="9">⁹ Department of Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa</footnote>
# Asymmetries in random motions of neutral hydrogen gas in spiral galaxies

P. Adamczyk<fnref n="1" />,<fnref n="2" />, P. Amram<fnref n="2" />, L. Chemin<fnref n="3" />, B. Epinat<fnref n="2" />,<fnref n="4" />, J. Braine<fnref n="5" />, F. Combes<fnref n="6" />, and W. J. G. de Blok<fnref n="7" />,<fnref n="8" />,<fnref n="9" />

Received 26 April 2023 / Accepted 8 June 2023

## ABSTRACT

**Context.** The velocity dispersion ellipsoid of gas in galactic discs is usually assumed to be isotropic. Under this approximation, no projection effect occurs in the random motions of gas, as traced by the line-of-sight velocity dispersion. However, it has been recently shown that random motions of the neutral hydrogen gas of the Triangulum galaxy (M 33) exhibit a bisymmetric perturbation which is aligned with the minor axis of the galaxy, suggesting a projection effect.

**Aims.** To investigate if perturbations in the velocity dispersion of nearby discs are comparable to those of M 33, the sample is extended to 32 galaxies from The H I Nearby Galaxy Survey (THINGS) and the Westerbork H I Survey of Spiral and Irregular Galaxies (WHISP).

**Methods.** We studied velocity asymmetries in the disc planes by performing Fourier transforms of high-resolution H I velocity dispersion maps corrected for beam-smearing effects, and we measured the amplitudes and phase angles of the Fourier harmonics.

**Results.** In all velocity dispersion maps, we find strong perturbations of first, second, and fourth orders. The strongest asymmetry is the bisymmetry, which is predominantly associated with the presence of spiral arms. The first order asymmetry is generally orientated close to the disc major axis, and the second and fourth order asymmetries are preferentially orientated along intermediate directions between the major and minor axes of the discs. These results are evidence that strong projection effects shape the H I velocity dispersion maps. The most likely source of systematic orientations is the anisotropy of velocities, through the projection of streaming motions that are stronger along one of the planar directions in the discs. Moreover, systematic phase angles of asymmetries in the H I velocity dispersion could arise from tilted velocity ellipsoids, that is when the velocities are correlated. We expect a larger incidence of correlation between the radial and tangential velocities of H I gas with $|ρ_{Rθ}| \sim 0.6$, which could be tested against the kinematics of the youngest stellar populations of the Milky Way.

**Conclusions.** H I velocity dispersions cannot be considered devoid of projection effects. The systematic orientations of asymmetries can be explained by the projection of unresolved streaming motions mainly arising from spiral arms. Our methodology is a powerful tool to constrain the dominant direction of streaming motions and thus the shape of the velocity ellipsoid of H I gas, which is de facto anisotropic at the angular scales probed by the observations. The next step is to study the shape of the velocity ellipsoids of molecular and ionised gas and their link with galaxy mass and/or morphology, in addition to extending the sample size.

**Key words.** galaxies: fundamental parameters – galaxies: kinematics and dynamics – galaxies: spiral – galaxies: structure – galaxies: ISM

## 1. Introduction

In discs of galaxies, large-scale perturbations such as bars, spiral arms, warps, or lopsidedness make the orbits non-circular. They generate streaming motions, which are asymmetric radial and tangential velocities in the stellar and gaseous media in galactic discs (e.g. Visser 1980). This is beautifully illustrated from an observational viewpoint via the kinematics of millions of individual stars in the Large Magellanic Cloud (LMC), for which the Gaia Collaboration used astrometry to deduce the velocity fields of both the radial and tangential components of

a disc galaxy (Gaia Collaboration 2016, 2018a,b, 2021a), and from which significant asymmetries are observed in the velocity maps, as caused by the bar and spiral arms of the LMC (Gaia Collaboration 2021b).

The study of kinematic asymmetries has become relevant in galactic dynamics, as it impacts our knowledge of the structure of dark matter halos and the comparison with simulations made in a cosmological context. Interstellar gas is thought to follow circular orbits as it undergoes dissipative collisions. If asymmetries are important, then there is no guarantee that rotation curves perfectly trace the mass distribution of galaxies.

<footnote n="1">¹ Centro de Astronomia, Universidad de Antofagasta, Avda. U. de Antofagasta, 02800 Antofagasta, Chile
e-mail: p.adamczyk@free.fr</footnote>

<footnote n="2">² Aix-Marseille Univ., CNRS, CNES, LAM, 38 Rue Frédéric Joliot Curie, 13338 Marseille, France
e-mail: philippe.amram@lam.fr</footnote>

<footnote n="3">³ Instituto de Astrofísica, Universidad Andes Bello, Fernandez Concha 700, Las Condes, Santiago RM, Chile
e-mail: astro.chemin@gmail.com</footnote>

<footnote n="4">⁴ Canada-France-Hawaii Telescope, 65-1238 Mamalahoa Highway, Kamuela, HI 96743, USA</footnote>

<footnote n="5">⁵ Laboratoire d'Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote n="6">⁶ Observatoire de Paris, LERMA, Collège de France, CNRS, PSL University, Sorbonne University, 75014 Paris, France</footnote>

<footnote n="7">⁷ ASTRON, Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PD Dwingeloo, The Netherlands</footnote>

<footnote n="8">⁸ Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands</footnote>

<footnote n="9">⁹ Department of Astronomy, University of Cape Town, Private Bag X3, Rondebosch 7701, South Africa</footnote>
Asymmetries could affect the inner slope of gaseous rotation curves and require a model including asymmetries to determine the scale parameters of dark matter halos (Hayashi & Navarro 2006; Oman et al. 2019). Strategies have been developed to quantify asymmetries in gaseous velocity fields in the last 25 years. One approach is to measure harmonic velocity components that overlap with the axisymmetric motions (Schoenmakers et al. 1997; Krajnović et al. 2006). Early works focussed on the large-scale kinematic lopsidedness in H I velocity fields, a perturbation of first order that causes the kinematic centre to drift with radius. In this case, rotation curves from the approaching and receding halves are not consistent with each other. This could be due to a misalignment between the axes of the disc and the host dark matter halo (Swaters et al. 1999). Trachternach et al. (2008) estimated asymmetries for a large sample of H I velocity fields to study the elongation of the gravitational potential. They found that the asymmetries do not alter the mass distribution inferred from the axisymmetric rotational motions. Spekkens & Sellwood (2007) measured bisymmetric perturbations in CO and Hα velocity fields to highlight the importance of bisymmetric flows caused by a stellar bar on the inner shape of the rotation curve. Chemin et al. (2016) modelled the harmonics observed in the Hα velocity field of a grand design spiral through direct derivation of gravitational potentials and assessed the impact of asymmetries on the structure of the host dark matter halo.

While studying streaming motions in ordered velocity fields has become routine, little is known about asymmetries in maps of apparently random motions of interstellar gas, as traced by the line-of-sight velocity dispersion ($\sigma_{los}$). The velocity dispersion can be used to assess the dynamical heating of discs, and to estimate the support due to pressure and how it compares with rotational motions (e.g. Combes & Becquart 1997; Koyama & Ostriker 2009; Bershady et al. 2010; Oh et al. 2015). The main reason dispersion has not been studied extensively is that $\sigma_{los}$ mixes turbulence, asymmetric drift, shear, superposition of multiple elliptical orbits, and unresolved motions on relatively large scales. The local velocity dispersion ellipsoid of gas is also usually considered isotropic because the gaseous component dissipates energy through collisions between clouds in all directions. Within this hypothesis, no projection effect occurs in the random motions of gas, as traced by $\sigma_{los}$, unlike that of collisionless stars. Furthermore, the projection of the velocity ellipsoid can be degenerate if the ellipsoid is anisotropic. The finite velocity and spatial resolution further complicate the interpretation of maps of $\sigma_{los}$.

In a recent study of the Local Group spiral M33, Chemin et al. (2020) found perturbations in the random motions through Fourier transforms of $\sigma_{los}$ maps of 21-cm H I data at various angular resolutions. Bisymmetry dominated in the outer half of the disc with a phase aligned with the minor axis of M33. The question of whether this is a typical characteristic of spiral galaxies in general remains open. One possibility could be that these observations are due to the influence of an anisotropic velocity ellipsoid, which harbours a dominant radial component, as suggested by these authors. The goal of this study is to determine whether what was observed in M33 is a general feature of spiral discs.

Our goal is to systematically search for asymmetries in velocity dispersion maps, identify their origin, and evaluate the consequences for galactic dynamics. In this paper, we present the first census of asymmetries in the velocity dispersion of H I gas in local, massive discs. The samples used for this analysis are described in Sect. 2. The asymmetries are measured through

fast Fourier transforms (FFTs) of velocity dispersion maps from interferometric data with a careful treatment of the effect of beam smearing (BS hereafter), as detailed in Sect. 3. The properties of the asymmetries are presented in Sect. 4, and the discussion of their origin is given in Sect. 5. Finally, we provide a synthesis and conclusions in Sect. 6.

## 2. Selection of a working sample

Several factors limit the modelling of spatially resolved velocity dispersions. The smearing induced by the finite telescope beam can produce systematic asymmetries in velocity dispersion maps (see Sect. 3). For a given instrumental configuration of spatial and spectral resolution and sampling, the BS effect increases with the distance and inclination of the galaxies.

For these reasons, our sample is from The high-resolution H I Nearby Galaxy Survey (THINGS, Walter et al. 2008). With galaxy distances between 2 and 15 Mpc, THINGS yields linear resolutions up to ~900 pc. The velocity resolution is 2.6 or 5.2 km s⁻<fnref n="1" />, and the pixel size 1.5", except for the galaxy NGC 2403 (1"). We have used the three moment maps of the integrated H I emission (0th moment), line-of-sight velocity and velocity dispersion (1st and 2nd moments, respectively), which are made available by the THINGS collaboration<fnref n="1" />. In particular, we used the data obtained with natural weighting. We however verified that the conclusions of this analysis are not changed by using higher resolution kinematics, from maps obtained with robust weighting. A discussion of the impact of the angular resolution can be found in Sect. 5.2.3.

To perform the Fourier analysis of velocity dispersion fields, robust constraints on the geometry of the discs are necessary, and particularly the variation of the inclination and the position angle of the discs as a function of galactocentric radius. The disc warping parameters for 19 of the 34 THINGS galaxies have been measured with tilted-ring models of the velocity fields by de Blok et al. (2008). Among several methods applied to obtain velocity fields from the H I data cubes, these authors adopted results from Hermite h3 polynomials fittings, due to their stability at low signal-to-noise ratio (S/N). They explained that two masks were successively applied to filter low-quality regions from the velocity fields. The first one consisted in rejecting H I profiles (i) for which the fitted maximum intensity was lower than $3\sigma_{ch}$, with $\sigma_{ch}$ being the average noise in the profile outside of the emission line, and (ii) for which the dispersion of the fitted function was lower than the channel separation. The second mask consisted in a sigma-clipping on the H I column density maps to suppress noise pixels and to exclude regions inside $r_{min}$ and outside $r_{max}$, which are respectively the innermost and outermost radii of the tilted-ring models fitted of the final velocity fields. They finally measured the warp parameters and rotation curves, adopting a sampling of two points per synthesised beam size.

We masked the moment maps following de Blok et al. (2008) and adopted the parameters of their tilted-ring models. This makes our gaseous and kinematic distributions, and the galaxy cylindrical frames fully consistent with their analysis. A minor difference with their study is that we considered sampling the radial profiles with one data-point per beam only, to minimise the correlation between adjacent rings. Finally, among the 19 galaxies, we selected the 15 massive, regular galaxies, excluding dwarfs or strongly disturbed galaxies (NGC 2366, DDO 154, IC 2574 and NGC 4826).

<footnote n="1">¹ https://www2.mpc-hd.mpg.de/THINGS/Data.html</footnote>
Asymmetries could affect the inner slope of gaseous rotation curves and require a model including asymmetries to determine the scale parameters of dark matter halos (Hayashi & Navarro 2006; Oman et al. 2019). Strategies have been developed to quantify asymmetries in gaseous velocity fields in the last 25 years. One approach is to measure harmonic velocity components that overlap with the axisymmetric motions (Schoenmakers et al. 1997; Krajnović et al. 2006). Early works focussed on the large-scale kinematic lopsidedness in H I velocity fields, a perturbation of first order that causes the kinematic centre to drift with radius. In this case, rotation curves from the approaching and receding halves are not consistent with each other. This could be due to a misalignment between the axes of the disc and the host dark matter halo (Swaters et al. 1999). Trachternach et al. (2008) estimated asymmetries for a large sample of H I velocity fields to study the elongation of the gravitational potential. They found that the asymmetries do not alter the mass distribution inferred from the axisymmetric rotational motions. Spekkens & Sellwood (2007) measured bisymmetric perturbations in CO and Hα velocity fields to highlight the importance of bisymmetric flows caused by a stellar bar on the inner shape of the rotation curve. Chemin et al. (2016) modelled the harmonics observed in the Hα velocity field of a grand design spiral through direct derivation of gravitational potentials and assessed the impact of asymmetries on the structure of the host dark matter halo.

While studying streaming motions in ordered velocity fields has become routine, little is known about asymmetries in maps of apparently random motions of interstellar gas, as traced by the line-of-sight velocity dispersion ($\sigma_{los}$). The velocity dispersion can be used to assess the dynamical heating of discs, and to estimate the support due to pressure and how it compares with rotational motions (e.g. Combes & Becquart 1997; Koyama & Ostriker 2009; Bershady et al. 2010; Oh et al. 2015). The main reason dispersion has not been studied extensively is that $\sigma_{los}$ mixes turbulence, asymmetric drift, shear, superposition of multiple elliptical orbits, and unresolved motions on relatively large scales. The local velocity dispersion ellipsoid of gas is also usually considered isotropic because the gaseous component dissipates energy through collisions between clouds in all directions. Within this hypothesis, no projection effect occurs in the random motions of gas, as traced by $\sigma_{los}$, unlike that of collisionless stars. Furthermore, the projection of the velocity ellipsoid can be degenerate if the ellipsoid is anisotropic. The finite velocity and spatial resolution further complicate the interpretation of maps of $\sigma_{los}$.

In a recent study of the Local Group spiral M33, Chemin et al. (2020) found perturbations in the random motions through Fourier transforms of $\sigma_{los}$ maps of 21-cm H I data at various angular resolutions. Bisymmetry dominated in the outer half of the disc with a phase aligned with the minor axis of M33. The question of whether this is a typical characteristic of spiral galaxies in general remains open. One possibility could be that these observations are due to the influence of an anisotropic velocity ellipsoid, which harbours a dominant radial component, as suggested by these authors. The goal of this study is to determine whether what was observed in M33 is a general feature of spiral discs.

Our goal is to systematically search for asymmetries in velocity dispersion maps, identify their origin, and evaluate the consequences for galactic dynamics. In this paper, we present the first census of asymmetries in the velocity dispersion of H I gas in local, massive discs. The samples used for this analysis are described in Sect. 2. The asymmetries are measured through

fast Fourier transforms (FFTs) of velocity dispersion maps from interferometric data with a careful treatment of the effect of beam smearing (BS hereafter), as detailed in Sect. 3. The properties of the asymmetries are presented in Sect. 4, and the discussion of their origin is given in Sect. 5. Finally, we provide a synthesis and conclusions in Sect. 6.

## 2. Selection of a working sample

Several factors limit the modelling of spatially resolved velocity dispersions. The smearing induced by the finite telescope beam can produce systematic asymmetries in velocity dispersion maps (see Sect. 3). For a given instrumental configuration of spatial and spectral resolution and sampling, the BS effect increases with the distance and inclination of the galaxies.

For these reasons, our sample is from The high-resolution H I Nearby Galaxy Survey (THINGS, Walter et al. 2008). With galaxy distances between 2 and 15 Mpc, THINGS yields linear resolutions up to ~900 pc. The velocity resolution is 2.6 or 5.2 km s⁻<fnref n="1" />, and the pixel size 1.5", except for the galaxy NGC 2403 (1"). We have used the three moment maps of the integrated H I emission (0th moment), line-of-sight velocity and velocity dispersion (1st and 2nd moments, respectively), which are made available by the THINGS collaboration<fnref n="1" />. In particular, we used the data obtained with natural weighting. We however verified that the conclusions of this analysis are not changed by using higher resolution kinematics, from maps obtained with robust weighting. A discussion of the impact of the angular resolution can be found in Sect. 5.2.3.

To perform the Fourier analysis of velocity dispersion fields, robust constraints on the geometry of the discs are necessary, and particularly the variation of the inclination and the position angle of the discs as a function of galactocentric radius. The disc warping parameters for 19 of the 34 THINGS galaxies have been measured with tilted-ring models of the velocity fields by de Blok et al. (2008). Among several methods applied to obtain velocity fields from the H I data cubes, these authors adopted results from Hermite h3 polynomials fittings, due to their stability at low signal-to-noise ratio (S/N). They explained that two masks were successively applied to filter low-quality regions from the velocity fields. The first one consisted in rejecting H I profiles (i) for which the fitted maximum intensity was lower than $3\sigma_{ch}$, with $\sigma_{ch}$ being the average noise in the profile outside of the emission line, and (ii) for which the dispersion of the fitted function was lower than the channel separation. The second mask consisted in a sigma-clipping on the H I column density maps to suppress noise pixels and to exclude regions inside $r_{min}$ and outside $r_{max}$, which are respectively the innermost and outermost radii of the tilted-ring models fitted of the final velocity fields. They finally measured the warp parameters and rotation curves, adopting a sampling of two points per synthesised beam size.

We masked the moment maps following de Blok et al. (2008) and adopted the parameters of their tilted-ring models. This makes our gaseous and kinematic distributions, and the galaxy cylindrical frames fully consistent with their analysis. A minor difference with their study is that we considered sampling the radial profiles with one data-point per beam only, to minimise the correlation between adjacent rings. Finally, among the 19 galaxies, we selected the 15 massive, regular galaxies, excluding dwarfs or strongly disturbed galaxies (NGC 2366, DDO 154, IC 2574 and NGC 4826).

<footnote n="1">¹ https://www2.mpc-hd.mpg.de/THINGS/Data.html</footnote>
<table><caption>**Table 1.** Parameters of the galaxy sample.</caption><thead><tr><th>Galaxy name (1)</th><th>α (J2000) (hh mm ss) (2)</th><th>δ (J2000) (dd mm ss) (3)</th><th>D (Mpc) (4)</th><th>V<sub>sys</sub> (km s<sup>-1</sup>) (5)</th><th>Incl. (deg) (6)</th><th>PA (deg) (7)</th><th>B<sub>maj</sub> (arcsec pc<sup>-1</sup>) (8)</th><th>B<sub>min</sub> (arcsec pc<sup>-1</sup>) (9)</th><th>B<sub>PA</sub> (deg) (10)</th><th>Δ<sub>V</sub> (km s<sup>-1</sup>) (11)</th></tr></thead><tbody><tr><td>NGC 925</td><td>02 27 16.5</td><td>+33 34 44</td><td>9.2</td><td>546</td><td>66</td><td>287</td><td>5.9/263</td><td>5.7/254</td><td>31</td><td>2.6</td></tr><tr><td>NGC 2403</td><td>07 36 51.1</td><td>+65 36 03</td><td>3.2</td><td>133</td><td>63</td><td>124</td><td>8.8/136</td><td>7.7/119</td><td>25</td><td>5.2</td></tr><tr><td>NGC 2841</td><td>09 22 02.6</td><td>+50 58 35</td><td>14.1</td><td>634</td><td>74</td><td>153</td><td>11.1/756</td><td>9.4/641</td><td>-12</td><td>5.2</td></tr><tr><td>NGC 2903</td><td>09 32 10.1</td><td>+21 30 04</td><td>8.9</td><td>556</td><td>65</td><td>204</td><td>15.3/659</td><td>13.3/573</td><td>-51</td><td>5.2</td></tr><tr><td>NGC 2976</td><td>09 47 15.3</td><td>+67 55 00</td><td>3.6</td><td>1</td><td>65</td><td>335</td><td>7.4/129</td><td>6.4/112</td><td>72</td><td>5.2</td></tr><tr><td>NGC 3031</td><td>09 55 33.1</td><td>+69 03 55</td><td>3.6</td><td>-40</td><td>59</td><td>330</td><td>12.9/225</td><td>12.4/216</td><td>80</td><td>2.6</td></tr><tr><td>NGC 3198</td><td>10 19 55.0</td><td>+45 32 59</td><td>13.8</td><td>661</td><td>72</td><td>215</td><td>13.0/867</td><td>11.6/774</td><td>-59</td><td>5.2</td></tr><tr><td>NGC 3521</td><td>11 05 48.6</td><td>-00 02 09</td><td>10.7</td><td>804</td><td>73</td><td>340</td><td>14.1/730</td><td>11.2/580</td><td>-62</td><td>5.2</td></tr><tr><td>NGC 3621</td><td>11 18 16.5</td><td>-32 48 51</td><td>6.6</td><td>729</td><td>65</td><td>345</td><td>15.9/508</td><td>10.2/326</td><td>4</td><td>5.2</td></tr><tr><td>NGC 3627</td><td>11 20 15.0</td><td>+12 59 30</td><td>9.3</td><td>708</td><td>62</td><td>173</td><td>10.6/477</td><td>8.9/400</td><td>-48</td><td>5.2</td></tr><tr><td>NGC 4736</td><td>12 50 53.0</td><td>+41 07 13</td><td>4.7</td><td>307</td><td>41</td><td>296</td><td>10.2/232</td><td>9.1/207</td><td>-23</td><td>5.2</td></tr><tr><td>NGC 5055</td><td>13 15 49.2</td><td>+42 01 45</td><td>10.1</td><td>497</td><td>59</td><td>102</td><td>10.1/493</td><td>8.7/425</td><td>-40</td><td>5.2</td></tr><tr><td>NGC 6946</td><td>20 34 52.2</td><td>+60 09 14</td><td>5.9</td><td>44</td><td>33</td><td>243</td><td>6.0/171</td><td>5.6/160</td><td>7</td><td>1.3</td></tr><tr><td>NGC 7331</td><td>22 37 04.1</td><td>+34 24 57</td><td>14.7</td><td>818</td><td>76</td><td>168</td><td>6.1/433</td><td>5.6/398</td><td>34</td><td>5.2</td></tr><tr><td>NGC 7793</td><td>23 57 49.7</td><td>-32 35 28</td><td>3.9</td><td>226</td><td>50</td><td>290</td><td>15.6/295</td><td>10.9/206</td><td>11</td><td>2.6</td></tr><tr><td>UGC 01256</td><td>01 47 53.9</td><td>+27 25 55</td><td>8.8</td><td>426</td><td>70</td><td>69</td><td>25.6/1090</td><td>11.6/494</td><td>-1.0</td><td>16.5</td></tr><tr><td>UGC 01913</td><td>02 27 16.9</td><td>+33 34 44</td><td>10.8</td><td>553</td><td>59</td><td>287</td><td>16.7/872</td><td>8.6/449</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 02455</td><td>02 59 42.5</td><td>+25 14 19</td><td>7.5</td><td>373</td><td>42</td><td>208</td><td>27.4/1417</td><td>11.2/580</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 04284</td><td>08 14 40.1</td><td>+49 03 42</td><td>10.7</td><td>547</td><td>60</td><td>170</td><td>13.1/855</td><td>10.2/666</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 04305</td><td>08 19 04.3</td><td>+70 43 18</td><td>5.3</td><td>158</td><td>51</td><td>195</td><td>12.3/607</td><td>11.6/572</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 04325</td><td>08 19 20.5</td><td>+50 00 35</td><td>10.3</td><td>506</td><td>68</td><td>60</td><td>15.1/490</td><td>11.8/383</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 04499</td><td>08 37 41.5</td><td>+51 39 09</td><td>13.5</td><td>687</td><td>81</td><td>151</td><td>14.8/937</td><td>11.7/741</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 05414</td><td>10 03 57.2</td><td>+40 45 27</td><td>10.2</td><td>604</td><td>54</td><td>216</td><td>17.7/583</td><td>11.4/375</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 05721</td><td>10 32 17.2</td><td>+27 40 08</td><td>6.6</td><td>532</td><td>62</td><td>273</td><td>29.7/1236</td><td>13.0/540</td><td>1.2</td><td>4.1</td></tr><tr><td>UGC 05789</td><td>10 39 09.5</td><td>+41 41 13</td><td>13.1</td><td>738</td><td>63</td><td>38</td><td>18.4/863</td><td>11.9/558</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07323</td><td>12 17 30.2</td><td>+45 37 09</td><td>8.6</td><td>516</td><td>52</td><td>23</td><td>17.5/728</td><td>12.7/528</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07766</td><td>12 35 57.7</td><td>+27 57 35</td><td>9.7</td><td>814</td><td>65</td><td>328</td><td>22.1/589</td><td>9.3/283</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07831</td><td>12 39 59.3</td><td>+61 36 33</td><td>4.4</td><td>146</td><td>70</td><td>302</td><td>13.0/773</td><td>8.6/511</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07853</td><td>12 41 32.9</td><td>+41 09 04</td><td>8.6</td><td>538</td><td>58</td><td>212</td><td>19.8/1006</td><td>12.9/655</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 08490</td><td>13 29 36.6</td><td>+58 25 14</td><td>5.5</td><td>202</td><td>59</td><td>185</td><td>13.5/490</td><td>11.3/410</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 11891</td><td>22 03 33.7</td><td>+43 44 56</td><td>12.3</td><td>461</td><td>43</td><td>130</td><td>16.5/854</td><td>11.4/590</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 12632</td><td>23 29 58.7</td><td>+40 59 25</td><td>10.5</td><td>422</td><td>37</td><td>20</td><td>18.4/472</td><td>12.2/313</td><td>0.0</td><td>4.1</td></tr></tbody></table>

**Notes.** Sources from THINGS are the NGC names, those from WHISP the UGC names. (1) Name of the object; (2) and (3) Right ascension and Declination (J2000) of the galactic centre; (4) Galaxy distance from Walter et al. (2008) for THINGS, and deduced from the systemic velocity taken in NED corrected from Virgo infall, assuming $H_0 = 67.8$ km s$^{-1}$ Mpc$^{-1}$ for WHISP; (5)–(7) Systemic velocity, mean inclination and position angle from de Blok et al. (2008); (8) and (9) Major and minor axes of the synthesised beam; (10) Beam position angle; (11) Spectral sampling.

Chemin et al. (2020) showed that the angular resolution did not affect the detection of the bisymmetry and/or anisotropy observed in H I velocity dispersion maps of M 33, as it was seen at 70, 100, and 490 pc resolution. To increase the size of our sample, we included sources from the Westerbork H I Survey of Spiral and Irregular Galaxies (WHISP; van der Hulst et al. 2001). The reduction pipeline of the WHISP observations provides data at resolutions of $14'' \times 14''/\sin\delta$, where $\delta$ is the galaxy declination, $30'' \times 30''$, and $60'' \times 60''$, for a spectral resolution of $2.06$ to $16.5$ km s$^{-1}$. At distances out to more than 100 Mpc, the implied linear resolution for many WHISP galaxies is beyond the kiloparsec scale. Considering the balance between spatial resolution and sensitivity, and, due to the necessity to work with an appropriate spatial resolution, we only selected the sources with a physical resolution better than 1.5 kpc (Adamczyk 2021). This corresponds to lower S/N measurements than for THINGS targets, and consequently to a smaller extent of the gas distribution and kinematics. We further required the neutral gas to cover seven times the size of the beam (as suggested in Bosma 1978, to derive rotation curves), in the same inclination range as THINGS, and with no obvious signs of tidal interactions. In the end, 17 WHISP galaxies were added to the THINGS sample. We produced moment maps for them by means of CAMEL, a Python

tool<fnref n="2" /> described in Epinat et al. (2012). We then performed tilted-ring models of WHISP data cubes using the package 3D-Based Analysis of Rotating Object via Line Observations (3D-Barolo, Di Teodoro & Fraternali 2015). This software creates a model data cube from input centre of the rings, inclination and position angle, systemic velocity, scale height of the disc, and velocity dispersion, which is fitted to the observations. To determine the kinematic centre, we first kept the inclination and position angle fixed, with the centre of radial rings and the systemic velocity as free parameters. Then, we fixed the centre and let the inclination and position angles vary to get the warp parameters.

The final sample of galaxies is presented in Table 1 which lists the coordinates, distance, systemic velocity, the mean inclination and position angle of the discs, and the properties of the synthesised beam. With mean source distances of 8 and 9 Mpc for THINGS and WHISP respectively, the mean beam $B_{maj}$ of 11" and 18" correspond to average linear resolutions of 470 pc and 80 pc, respectively. The THINGS sample is the main reference and the WHISP sample is used to check the observed trends are not linked to a particular telescope. Section 4 shows that

<footnote n="2">² CAMEL stands for Cube Analysis: Moment maps of Emission Lines, see https://gitlab.lam.fr/bepinat/CAMEL</footnote>
<table><caption>**Table 1.** Parameters of the galaxy sample.</caption><thead><tr><th>Galaxy name<br>(1)</th><th>α (J2000)<br>(hh mm ss)<br>(2)</th><th>δ (J2000)<br>(dd mm ss)<br>(3)</th><th>D<br>(Mpc)<br>(4)</th><th>V<sub>sys</sub><br>(km s<sup>-1</sup>)<br>(5)</th><th>Incl.<br>(deg)<br>(6)</th><th>PA<br>(deg)<br>(7)</th><th>B<sub>maj</sub><br>(arcsec pc<sup>-1</sup>)<br>(8)</th><th>B<sub>min</sub><br>(arcsec pc<sup>-1</sup>)<br>(9)</th><th>B<sub>PA</sub><br>(deg)<br>(10)</th><th>Δ<sub>V</sub><br>(km s<sup>-1</sup>)<br>(11)</th></tr></thead><tbody><tr><td>NGC 925</td><td>02 27 16.5</td><td>+33 34 44</td><td>9.2</td><td>546</td><td>66</td><td>287</td><td>5.9/263</td><td>5.7/254</td><td>31</td><td>2.6</td></tr><tr><td>NGC 2403</td><td>07 36 51.1</td><td>+65 36 03</td><td>3.2</td><td>133</td><td>63</td><td>124</td><td>8.8/136</td><td>7.7/119</td><td>25</td><td>5.2</td></tr><tr><td>NGC 2841</td><td>09 22 02.6</td><td>+50 58 35</td><td>14.1</td><td>634</td><td>74</td><td>153</td><td>11.1/756</td><td>9.4/641</td><td>-12</td><td>5.2</td></tr><tr><td>NGC 2903</td><td>09 32 10.1</td><td>+21 30 04</td><td>8.9</td><td>556</td><td>65</td><td>204</td><td>15.3/659</td><td>13.3/573</td><td>-51</td><td>5.2</td></tr><tr><td>NGC 2976</td><td>09 47 15.3</td><td>+67 55 00</td><td>3.6</td><td>1</td><td>65</td><td>335</td><td>7.4/129</td><td>6.4/112</td><td>72</td><td>5.2</td></tr><tr><td>NGC 3031</td><td>09 55 33.1</td><td>+69 03 55</td><td>3.6</td><td>-40</td><td>59</td><td>330</td><td>12.9/225</td><td>12.4/216</td><td>80</td><td>2.6</td></tr><tr><td>NGC 3198</td><td>10 19 55.0</td><td>+45 32 59</td><td>13.8</td><td>661</td><td>72</td><td>215</td><td>13.0/867</td><td>11.6/774</td><td>-59</td><td>5.2</td></tr><tr><td>NGC 3521</td><td>11 05 48.6</td><td>-00 02 09</td><td>10.7</td><td>804</td><td>73</td><td>340</td><td>14.1/730</td><td>11.2/580</td><td>-62</td><td>5.2</td></tr><tr><td>NGC 3621</td><td>11 18 16.5</td><td>-32 48 51</td><td>6.6</td><td>729</td><td>65</td><td>345</td><td>15.9/508</td><td>10.2/326</td><td>4</td><td>5.2</td></tr><tr><td>NGC 3627</td><td>11 20 15.0</td><td>+12 59 30</td><td>9.3</td><td>708</td><td>62</td><td>173</td><td>10.6/477</td><td>8.9/400</td><td>-48</td><td>5.2</td></tr><tr><td>NGC 4736</td><td>12 50 53.0</td><td>+41 07 13</td><td>4.7</td><td>307</td><td>41</td><td>296</td><td>10.2/232</td><td>9.1/207</td><td>-23</td><td>5.2</td></tr><tr><td>NGC 5055</td><td>13 15 49.2</td><td>+42 01 45</td><td>10.1</td><td>497</td><td>59</td><td>102</td><td>10.1/493</td><td>8.7/425</td><td>-40</td><td>5.2</td></tr><tr><td>NGC 6946</td><td>20 34 52.2</td><td>+60 09 14</td><td>5.9</td><td>44</td><td>33</td><td>243</td><td>6.0/171</td><td>5.6/160</td><td>7</td><td>1.3</td></tr><tr><td>NGC 7331</td><td>22 37 04.1</td><td>+34 24 57</td><td>14.7</td><td>818</td><td>76</td><td>168</td><td>6.1/433</td><td>5.6/398</td><td>34</td><td>5.2</td></tr><tr><td>NGC 7793</td><td>23 57 49.7</td><td>-32 35 28</td><td>3.9</td><td>226</td><td>50</td><td>290</td><td>15.6/295</td><td>10.9/206</td><td>11</td><td>2.6</td></tr><tr><td>UGC 01256</td><td>01 47 53.9</td><td>+27 25 55</td><td>8.8</td><td>426</td><td>70</td><td>69</td><td>25.6/1090</td><td>11.6/494</td><td>-1.0</td><td>16.5</td></tr><tr><td>UGC 01913</td><td>02 27 16.9</td><td>+33 34 44</td><td>10.8</td><td>553</td><td>59</td><td>287</td><td>16.7/872</td><td>8.6/449</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 02455</td><td>02 59 42.5</td><td>+25 14 19</td><td>7.5</td><td>373</td><td>42</td><td>208</td><td>27.4/1417</td><td>11.2/580</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 04284</td><td>08 14 40.1</td><td>+49 03 42</td><td>10.7</td><td>547</td><td>60</td><td>170</td><td>13.1/855</td><td>10.2/666</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 04305</td><td>08 19 04.3</td><td>+70 43 18</td><td>5.3</td><td>158</td><td>51</td><td>195</td><td>12.3/607</td><td>11.6/572</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 04325</td><td>08 19 20.5</td><td>+50 00 35</td><td>10.3</td><td>506</td><td>68</td><td>60</td><td>15.1/490</td><td>11.8/383</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 04499</td><td>08 37 41.5</td><td>+51 39 09</td><td>13.5</td><td>687</td><td>81</td><td>151</td><td>14.8/937</td><td>11.7/741</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 05414</td><td>10 03 57.2</td><td>+40 45 27</td><td>10.2</td><td>604</td><td>54</td><td>216</td><td>17.7/583</td><td>11.4/375</td><td>-1.0</td><td>2.1</td></tr><tr><td>UGC 05721</td><td>10 32 17.2</td><td>+27 40 08</td><td>6.6</td><td>532</td><td>62</td><td>273</td><td>29.7/1236</td><td>13.0/540</td><td>1.2</td><td>4.1</td></tr><tr><td>UGC 05789</td><td>10 39 09.5</td><td>+41 41 13</td><td>13.1</td><td>738</td><td>63</td><td>38</td><td>18.4/863</td><td>11.9/558</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07323</td><td>12 17 30.2</td><td>+45 37 09</td><td>8.6</td><td>516</td><td>52</td><td>23</td><td>17.5/728</td><td>12.7/528</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07766</td><td>12 35 57.7</td><td>+27 57 35</td><td>9.7</td><td>814</td><td>65</td><td>328</td><td>22.1/589</td><td>9.3/283</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07831</td><td>12 39 59.3</td><td>+61 36 33</td><td>4.4</td><td>146</td><td>70</td><td>302</td><td>13.0/773</td><td>8.6/511</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 07853</td><td>12 41 32.9</td><td>+41 09 04</td><td>8.6</td><td>538</td><td>58</td><td>212</td><td>19.8/1006</td><td>12.9/655</td><td>0.0</td><td>4.1</td></tr><tr><td>UGC 08490</td><td>13 29 36.6</td><td>+58 25 14</td><td>5.5</td><td>202</td><td>59</td><td>185</td><td>13.5/490</td><td>11.3/410</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 11891</td><td>22 03 33.7</td><td>+43 44 56</td><td>12.3</td><td>461</td><td>43</td><td>130</td><td>16.5/854</td><td>11.4/590</td><td>-1.0</td><td>4.1</td></tr><tr><td>UGC 12632</td><td>23 29 58.7</td><td>+40 59 25</td><td>10.5</td><td>422</td><td>37</td><td>20</td><td>18.4/472</td><td>12.2/313</td><td>0.0</td><td>4.1</td></tr></tbody></table>

**Notes.** Sources from THINGS are the NGC names, those from WHISP the UGC names. (1) Name of the object; (2) and (3) Right ascension and Declination (J2000) of the galactic centre; (4) Galaxy distance from Walter et al. (2008) for THINGS, and deduced from the systemic velocity taken in NED corrected from Virgo infall, assuming $H_0 = 67.8$ km s⁻¹ Mpc⁻¹ for WHISP; (5)–(7) Systemic velocity, mean inclination and position angle from de Blok et al. (2008); (8) and (9) Major and minor axes of the synthesised beam; (10) Beam position angle; (11) Spectral sampling.

Chemin et al. (2020) showed that the angular resolution did not affect the detection of the bisymmetry and/or anisotropy observed in H I velocity dispersion maps of M 33, as it was seen at 70, 100, and 490 pc resolution. To increase the size of our sample, we included sources from the Westerbork H I Survey of Spiral and Irregular Galaxies (WHISP; van der Hulst et al. 2001). The reduction pipeline of the WHISP observations provides data at resolutions of 14"×14"/sin δ, where δ is the galaxy declination, 30"×30", and 60"×60", for a spectral resolution of 2.06 to 16.5 km s⁻¹. At distances out to more than 100 Mpc, the implied linear resolution for many WHISP galaxies is beyond the kiloparsec scale. Considering the balance between spatial resolution and sensitivity, and, due to the necessity to work with an appropriate spatial resolution, we only selected the sources with a physical resolution better than 1.5 kpc (Adamczyk 2021). This corresponds to lower S/N measurements than for THINGS targets, and consequently to a smaller extent of the gas distribution and kinematics. We further required the neutral gas to cover seven times the size of the beam (as suggested in Bosma 1978, to derive rotation curves), in the same inclination range as THINGS, and with no obvious signs of tidal interactions. In the end, 17 WHISP galaxies were added to the THINGS sample. We produced moment maps for them by means of CAMEL, a Python

tool<fnref n="2" /> described in Epinat et al. (2012). We then performed tilted-ring models of WHISP data cubes using the package 3D-Based Analysis of Rotating Object via Line Observations (3D-Barolo, Di Teodoro & Fraternali 2015). This software creates a model data cube from input centre of the rings, inclination and position angle, systemic velocity, scale height of the disc, and velocity dispersion, which is fitted to the observations. To determine the kinematic centre, we first kept the inclination and position angle fixed, with the centre of radial rings and the systemic velocity as free parameters. Then, we fixed the centre and let the inclination and position angles vary to get the warp parameters.

The final sample of galaxies is presented in Table 1 which lists the coordinates, distance, systemic velocity, the mean inclination and position angle of the discs, and the properties of the synthesised beam. With mean source distances of 8 and 9 Mpc for THINGS and WHISP respectively, the mean beam $B_{maj}$ of 11" and 18" correspond to average linear resolutions of 470 pc and 80 pc, respectively. The THINGS sample is the main reference and the WHISP sample is used to check the observed trends are not linked to a particular telescope. Section 4 shows that

<footnote n="2">² CAMEL stands for Cube Analysis: Moment maps of Emission Lines, see https://gitlab.lam.fr/bepinat/CAMEL</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. Density and velocity maps of NGC 2841. From left to right we show: the observed flux density map, the observed velocity field, the observed velocity dispersion, the BS model, and the velocity dispersion map corrected from the BS effect.</figcaption></figure>

the WHISP data fully support the trends found with THINGS with data only.

## 3. Methodology

Observations yield the integrated flux of gas emission lines, the line-of-sight velocity and velocity dispersion. This observed dispersion cannot be modelled directly as it results from the convolution of the random motion in the plane of galaxies ($\sigma$) with the Line Spread Function of the instrumental device ($\sigma_{LSF}$) and the broadening from thermal processes occurring in the gas ($\sigma_{T}$). The observed dispersion is then expressed by:

$$ \sigma_{los} = (\sigma^2 + \sigma_{LSF}^2 + \sigma_{T}^2)^{1/2}. \qquad (1) $$

Additionally, due to the finite angular resolution of observations, $\sigma$ is affected by a smearing effect. This comes from unresolved velocity gradients in the line-of-sight velocity, which broadens emission lines.

Below we describe how we modelled the BS effect, and how a corrected velocity dispersion has been inferred. Then, we describe how Fourier transforms constrain the amplitudes and phases of asymmetries in the velocity dispersion maps.

### 3.1. Modelling the effect of BS in a velocity dispersion map

Beam smearing impacts the flux, velocity field, and velocity dispersion maps extracted from datacubes in a different way. Its main effect on velocity fields is to weaken the inner gradient of the rotation curve by mixing data from adjacent regions in the disc. de Blok et al. (2008) showed that BS does not have a significant impact on THINGS galaxies rotation velocities. They built mock datacubes, smoothed the cubes to the THINGS resolution, and derived the velocity fields and rotation curves from the mock datacubes. Compared to the input used to make the datacube, deviations were less than ~1 km s⁻¹.

For velocity dispersion maps, one expects BS to be more prominent where large variations of radial velocities are observed locally, for instance due to rotation curve inner gradients and to variations of $\cos(\theta)$, $\theta$ being the azimuth angle (see Eq. (C.1)). We need to model the impact of BS on velocity dispersion maps by taking into account the variations of velocities in the two spatial dimensions. We use the tool MocKinG<fnref n="3" /> based on an analytical formula presented in Epinat et al. (2010) who derived the observed velocity dispersion in a galaxy following:

$$ \sigma^2 = \sigma_{corr}^2 + \sigma_{bs}^2, \qquad (2) $$

$$ \sigma_{corr}^2 = \frac{\iint_{pix} [\sigma_{loc}^2 M] \otimes \text{PSF} \, ds}{\iint_{pix} M \otimes \text{PSF} \, ds}, \qquad (3) $$

and

$$ \sigma_{bs}^2 = \frac{\iint_{pix} [V_{los}^2 M] \otimes \text{PSF} \, ds}{\iint_{pix} M \otimes \text{PSF} \, ds} - \left( \frac{\iint_{pix} [V_{los} M] \otimes \text{PSF} \, ds}{\iint_{pix} M \otimes \text{PSF} \, ds} \right)^2, \qquad (4) $$

where $\sigma_{loc}$ is the local velocity dispersion map, $V_{los}$ is the line-of-sight velocity field, $M$ is the line flux map, $\otimes$PSF represents the two-dimension convolution by the PSF, and $\iint_{pix} ds$ integrates over the surface of the pixel. Equation (3) shows that $\sigma_{corr}$ is impacted by BS on the local velocity dispersion $\sigma_{loc}$ whereas Eq. (4) shows that unresolved velocity shears in the first moment map create an artificial line broadening. Our goal is to study $\sigma_{corr}$ by correcting the observed velocity dispersion $\sigma$ for the term $\sigma_{bs}$ of Eq. (4). The maps of first two moments ($M$ and $V_{los}$) in this equation should be at high-resolution and free from BS to properly account for velocity variations inside the beam and pixels. In practice, such high-resolution maps do not exist so the observed flux is usually used for $M$, and either a high-resolution model (see e.g. Epinat et al. 2010) or an observed velocity field are used for $V_{los}$. We used MocKinG with the observed H I flux maps and velocity fields to compute the BS correction. It is thus asymmetric by construction and accounts for both circular and non-circular motions from the kinematics. This leads to a fair correction, although not perfect, and is discussed in Sect. 5.2.3, 5.2.4, 5.3.1, and Appendix B. These BS dispersion maps were modelled assuming a bi-dimensional Gaussian synthesised beam (see Table 1). The effect of the observed dirty beams, which are not entirely elliptical, is addressed in Sect. 5.2.2.

An example of BS modelling with the galaxy NGC 2841 is shown in Fig. 1. This fast-rotating galaxy (second panel) has an inclination of ~74°. The fourth panel shows the corresponding $\sigma_{bs}$, which, once subtracted to the observed dispersion (middle panel), yields the corrected velocity dispersion map (right panel) from which asymmetries are measured. The general behaviour of the smearing effect is thus that $\sigma_{bs}$ is larger at low radius and exhibits a typical X-shape pattern. On average, $\sigma_{bs}$ decreases with radius, but non-negligible values are also observed near non-axisymmetric density and velocity features (see Sect. 5.2.1), underlining the benefits of modelling the smearing from 2D data.

<footnote n="3">³ MocKinG stands for modelling Kinematics of Galaxies, see https://gitlab.lam.fr/bepinat/MocKinG</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1. Density and velocity maps of NGC 2841. From left to right we show: the observed flux density map, the observed velocity field, the observed velocity dispersion, the BS model, and the velocity dispersion map corrected from the BS effect.</figcaption></figure>

the WHISP data fully support the trends found with THINGS with data only.

## 3. Methodology

Observations yield the integrated flux of gas emission lines, the line-of-sight velocity and velocity dispersion. This observed dispersion cannot be modelled directly as it results from the convolution of the random motion in the plane of galaxies ($\sigma$) with the Line Spread Function of the instrumental device ($\sigma_{LSF}$) and the broadening from thermal processes occurring in the gas ($\sigma_{T}$). The observed dispersion is then expressed by:

$$ \sigma_{\text{los}} = (\sigma^2 + \sigma_{\text{LSF}}^2 + \sigma_{\text{T}}^2)^{1/2}. \qquad (1) $$

Additionally, due to the finite angular resolution of observations, $\sigma$ is affected by a smearing effect. This comes from unresolved velocity gradients in the line-of-sight velocity, which broadens emission lines.

Below we describe how we modelled the BS effect, and how a corrected velocity dispersion has been inferred. Then, we describe how Fourier transforms constrain the amplitudes and phases of asymmetries in the velocity dispersion maps.

### 3.1. Modelling the effect of BS in a velocity dispersion map

Beam smearing impacts the flux, velocity field, and velocity dispersion maps extracted from datacubes in a different way. Its main effect on velocity fields is to weaken the inner gradient of the rotation curve by mixing data from adjacent regions in the disc. de Blok et al. (2008) showed that BS does not have a significant impact on THINGS galaxies rotation velocities. They built mock datacubes, smoothed the cubes to the THINGS resolution, and derived the velocity fields and rotation curves from the mock datacubes. Compared to the input used to make the datacube, deviations were less than ~1 km s⁻¹.

For velocity dispersion maps, one expects BS to be more prominent where large variations of radial velocities are observed locally, for instance due to rotation curve inner gradients and to variations of $\cos(\theta)$, $\theta$ being the azimuth angle (see Eq. (C.1)). We need to model the impact of BS on velocity dispersion maps by taking into account the variations of velocities in the two spatial dimensions. We use the tool MocKinG<fnref n="3" /> based on an analytical formula presented in Epinat et al. (2010) who derived the observed velocity dispersion in a galaxy following:

$$ \sigma^2 = \sigma_{\text{corr}}^2 + \sigma_{\text{bs}}^2, \qquad (2) $$

$$ \sigma_{\text{corr}}^2 = \frac{\iint_{\text{pix}} [\sigma_{\text{loc}}^2 M] \otimes \text{PSF} \, \text{d}s}{\iint_{\text{pix}} M \otimes \text{PSF} \, \text{d}s}, \qquad (3) $$

and

$$ \sigma_{\text{bs}}^2 = \frac{\iint_{\text{pix}} [V_{\text{los}}^2 M] \otimes \text{PSF} \, \text{d}s}{\iint_{\text{pix}} M \otimes \text{PSF} \, \text{d}s} - \left( \frac{\iint_{\text{pix}} [V_{\text{los}} M] \otimes \text{PSF} \, \text{d}s}{\iint_{\text{pix}} M \otimes \text{PSF} \, \text{d}s} \right)^2, \qquad (4) $$

where $\sigma_{\text{loc}}$ is the local velocity dispersion map, $V_{\text{los}}$ is the line-of-sight velocity field, $M$ is the line flux map, $\otimes$PSF represents the two-dimension convolution by the PSF, and $\iint_{\text{pix}}$ ds integrates over the surface of the pixel. Equation (3) shows that $\sigma_{\text{corr}}$ is impacted by BS on the local velocity dispersion $\sigma_{\text{loc}}$ whereas Eq. (4) shows that unresolved velocity shears in the first moment map create an artificial line broadening. Our goal is to study $\sigma_{\text{corr}}$ by correcting the observed velocity dispersion $\sigma$ for the term $\sigma_{\text{bs}}$ of Eq. (4). The maps of first two moments ($M$ and $V_{\text{los}}$) in this equation should be at high-resolution and free from BS to properly account for velocity variations inside the beam and pixels. In practice, such high-resolution maps do not exist so the observed flux is usually used for $M$, and either a high-resolution model (see e.g. Epinat et al. 2010) or an observed velocity field are used for $V_{\text{los}}$. We used MocKinG with the observed H I flux maps and velocity fields to compute the BS correction. It is thus asymmetric by construction and accounts for both circular and non-circular motions from the kinematics. This leads to a fair correction, although not perfect, and is discussed in Sect. 5.2.3, 5.2.4, 5.3.1, and Appendix B. These BS dispersion maps were modelled assuming a bi-dimensional Gaussian synthesised beam (see Table 1). The effect of the observed dirty beams, which are not entirely elliptical, is addressed in Sect. 5.2.2.

An example of BS modelling with the galaxy NGC 2841 is shown in Fig. 1. This fast-rotating galaxy (second panel) has an inclination of ~74°. The fourth panel shows the corresponding $\sigma_{\text{bs}}$, which, once subtracted to the observed dispersion (middle panel), yields the corrected velocity dispersion map (right panel) from which asymmetries are measured. The general behaviour of the smearing effect is thus that $\sigma_{\text{bs}}$ is larger at low radius and exhibits a typical X-shape pattern. On average, $\sigma_{\text{bs}}$ decreases with radius, but non-negligible values are also observed near non-axisymmetric density and velocity features (see Sect. 5.2.1), underlining the benefits of modelling the smearing from 2D data.

<footnote n="3">³ MocKinG stands for modelling Kinematics of Galaxies, see https://gitlab.lam.fr/bepinat/MocKinG</footnote>
<figure><img src="image_3.png" /><figcaption>Fig. 2. Distribution of the fraction of missing values for the radial bins in the THINGS galaxies.</figcaption></figure>

## 3.2. Fourier series modelling of velocity dispersion maps

In order to quantify asymmetries, we performed FFTs of the maps of H I random motions. Discrete FFTs of the corrected $\sigma_{\text{los}}^2$ maps were performed using the FFTPACK package of the SciPy libraries in Python (see Appendix A). In cylindrical coordinates, the discrete FFT of the velocity dispersion is:

$$ \sigma_{\text{asym}}^2 = \sum_{k=0}^{N-1} \sigma_k^2 \cos(k(\theta - \phi_k)), \qquad (5) $$

where $\sigma_k$ and $\phi_k$ are the amplitude and phase angle of an asymmetry of order k, N is the number of orders in the decomposition, equivalent to the number of elements along the considered ring, and $\theta$ is the azimuthal angle in the galaxy plane measured from the semi-major axis of the receding disc half. With this FFT formalism, we make the choice of studying the azimuthal asymmetry, that is we cannot study the asymmetry with respect to the disc plane in the vertical direction. Nevertheless, the present study includes radial, azimuthal and vertical components projected on the sky plane, assuming either that contributions along the direction orthogonal to the disc plane are averaged across the line of sight, which is in practice more valid for face-on than for edge-on galaxies, or that discs are infinitely thin. We use the squared velocity dispersion because $\sigma_{\text{los}}^2$ is a simple sum of quadratic terms (Eqs. (1) and (2), see also Appendix C, Eq. (C.12)). Consistent results are obtained when FFTs of the linear dispersion are calculated.

In practice, we decomposed a galactic disc into a series of concentric rings whose geometry was defined by the tilted-ring models (Sect. 2). Given the different angular sizes of the galaxies, the adopted ring width of one full width at half beam power leads to differing numbers of rings per galaxy, from 7 for NGC 3627 to more than 140 for NGC 2403. Such a non-uniformity in the number of radial bins among galaxies has no visible consequence on the results (see Sect. 4).

For each dispersion map, we first subtracted $\sigma_{\text{LSF}}^2$ and $\sigma_{\text{T}}^2$ from the squared observed dispersion, as well as the corresponding 2D map of $\sigma_{\text{bs}}^2$ described in Sect. 3.1. All pixels with resulting negative quadratic velocity dispersion after these subtractions were discarded from the maps at this stage of the process,

because such values are unphysical. The instrumental dispersion is $\sigma_{\text{LSF}} = \Delta_V/2.35$ and $\sigma_{\text{LSF}} = 1.2\Delta_V/2.35$ for the THINGS and WHISP galaxies, respectively, where $\Delta_V$ is the velocity channel width listed in Table 1. H I is a mixture of cool (~100 K) and warm (~5000–8000 K) gas. In M 33, the H I velocity dispersion was sometimes narrower than $\sigma_T$ if the warm gas was assumed to dominate (Chemin et al. 2020). Hence, we do not consider $\sigma_T$ here. We note that this has no consequence hereafter, because considering gas as a warm neutral medium is equivalent to subtracting quadratically $\sigma_T \sim 6 \text{ km s}^{-1}$ from the axisymmetric term $\sigma_0$. This latter term, measured as the mean dispersion of a given ring, is then subtracted quadratically from all pixel values inside the considered ring. Therefore, the observed $\sigma_{\text{asym}}^2 - \sigma_0^2$ are centred on 0, and can be negative (see Eq. (5)). We also point out that $\sigma_T$ could be asymmetric and vary over the disc. In such a case, fluctuations of $\sigma_T$ as a function of the position have been accounted in the measurements of $\sigma_{\text{asym}}$, though it is not possible to disentangle these effects from those arising from other local motions of, for example, gravitational origins, without being able to measure locally the gas temperature.

We then sorted the squared velocity dispersions with increasing values of azimuth and apply the FFT, leading to a harmonic decomposition with $N/2$ terms, where N is the number of pixels in the considered radial bin. Incomplete coverage of azimuths as caused by missing pixels (not-a-number values) could be seen by the FFT as artificial perturbations. In the THINGS sample, about 77% of the 883 available tilted rings show less than 1% of missing values of all available pixels, and 95% less than 4% of missing azimuthal angles (Fig. 2). The missing factor is thus low, and we verified that it has no impact in the analysis. Within WHISP data, we rejected 30% of the initial 203 rings because they had more than 40% of missing values. Within the 144 remaining rings, 15% of the pixels have missing values. Even though no impact was detected, we replaced the missing values by the azimuthally averaged dispersion before the modelling.

## 3.3. Generating toy models to study the impact of unresolved ordered velocity variations

As discussed in Sect. 3.1, BS is present in our data and is corrected in this study using observed line flux maps and velocity fields, which limits the accuracy of the correction. In order to study the impact of residual BS effects, we perform toy models for cases with barely resolved motions due to large-scale axisymmetric rotation of various strengths, with and without additional asymmetric velocity perturbations on local scales, under both isotropic and anisotropic hypotheses. We build toy models to produce data cubes, velocity, and $\sigma_{\text{los}}$ fields, and calculate the FFT of $\sigma_{\text{los}}$. This enables us to generate different configurations and identify the conditions for which systematic values in the distributions of $\sigma_k$ and/or $\phi_k$ could occur. Our toy models have no dynamical basis, yet they are very useful to assess the effects of BS, anisotropy, and streamings in $\sigma_{\text{los}}$. The full description of the toy models is presented in Appendix B.

We produced mock datacubes of 400 × 400 pixels (scale of 1", or ~50 pc at a typical distance of 10 Mpc), with 200 spectral elements with a $3 \text{ km s}^{-1}$ velocity sampling using $5 \times 10^6$ uniformly distributed points. To first order, the galaxy is assumed to be an axisymmetric rotating disc to which velocity perturbations can be added. Two rotation curves, one with a weak velocity gradient and a moderate velocity plateau, the other with a steep inner gradient quickly reaching a high velocity plateau, were used. Sharp planar velocity perturbations are produced by a bisymmetric spiral pattern, with five possible inner angles. No
<figure><img src="image_3.png" /><figcaption>Fig. 2. Distribution of the fraction of missing values for the radial bins in the THINGS galaxies.</figcaption></figure>

## 3.2. Fourier series modelling of velocity dispersion maps

In order to quantify asymmetries, we performed FFTs of the maps of H I random motions. Discrete FFTs of the corrected $\sigma_{\text{los}}^2$ maps were performed using the FFTPACK package of the SciPy libraries in Python (see Appendix A). In cylindrical coordinates, the discrete FFT of the velocity dispersion is:

$$ \sigma_{\text{asym}}^2 = \sum_{k=0}^{N-1} \sigma_k^2 \cos(k(\theta - \phi_k)), \qquad (5) $$

where $\sigma_k$ and $\phi_k$ are the amplitude and phase angle of an asymmetry of order $k$, $N$ is the number of orders in the decomposition, equivalent to the number of elements along the considered ring, and $\theta$ is the azimuthal angle in the galaxy plane measured from the semi-major axis of the receding disc half. With this FFT formalism, we make the choice of studying the azimuthal asymmetry, that is we cannot study the asymmetry with respect to the disc plane in the vertical direction. Nevertheless, the present study includes radial, azimuthal and vertical components projected on the sky plane, assuming either that contributions along the direction orthogonal to the disc plane are averaged across the line of sight, which is in practice more valid for face-on than for edge-on galaxies, or that discs are infinitely thin. We use the squared velocity dispersion because $\sigma_{\text{los}}^2$ is a simple sum of quadratic terms (Eqs. (1) and (2), see also Appendix C, Eq. (C.12)). Consistent results are obtained when FFTs of the linear dispersion are calculated.

In practice, we decomposed a galactic disc into a series of concentric rings whose geometry was defined by the tilted-ring models (Sect. 2). Given the different angular sizes of the galaxies, the adopted ring width of one full width at half beam power leads to differing numbers of rings per galaxy, from 7 for NGC 3627 to more than 140 for NGC 2403. Such a non-uniformity in the number of radial bins among galaxies has no visible consequence on the results (see Sect. 4).

For each dispersion map, we first subtracted $\sigma_{\text{LSF}}^2$ and $\sigma_{\text{T}}^2$ from the squared observed dispersion, as well as the corresponding 2D map of $\sigma_{\text{bs}}^2$ described in Sect. 3.1. All pixels with resulting negative quadratic velocity dispersion after these subtractions were discarded from the maps at this stage of the process,

because such values are unphysical. The instrumental dispersion is $\sigma_{\text{LSF}} = \Delta_V/2.35$ and $\sigma_{\text{LSF}} = 1.2\Delta_V/2.35$ for the THINGS and WHISP galaxies, respectively, where $\Delta_V$ is the velocity channel width listed in Table 1. H I is a mixture of cool (~100 K) and warm (~5000–8000 K) gas. In M 33, the H I velocity dispersion was sometimes narrower than $\sigma_T$ if the warm gas was assumed to dominate (Chemin et al. 2020). Hence, we do not consider $\sigma_T$ here. We note that this has no consequence hereafter, because considering gas as a warm neutral medium is equivalent to subtracting quadratically $\sigma_T \sim 6 \text{ km s}^{-1}$ from the axisymmetric term $\sigma_0$. This latter term, measured as the mean dispersion of a given ring, is then subtracted quadratically from all pixel values inside the considered ring. Therefore, the observed $\sigma_{\text{asym}}^2 - \sigma_0^2$ are centred on 0, and can be negative (see Eq. (5)). We also point out that $\sigma_T$ could be asymmetric and vary over the disc. In such a case, fluctuations of $\sigma_T$ as a function of the position have been accounted in the measurements of $\sigma_{\text{asym}}$, though it is not possible to disentangle these effects from those arising from other local motions of, for example, gravitational origins, without being able to measure locally the gas temperature.

We then sorted the squared velocity dispersions with increasing values of azimuth and apply the FFT, leading to a harmonic decomposition with $N/2$ terms, where $N$ is the number of pixels in the considered radial bin. Incomplete coverage of azimuths as caused by missing pixels (not-a-number values) could be seen by the FFT as artificial perturbations. In the THINGS sample, about 77% of the 883 available tilted rings show less than 1% of missing values of all available pixels, and 95% less than 4% of missing azimuthal angles (Fig. 2). The missing factor is thus low, and we verified that it has no impact in the analysis. Within WHISP data, we rejected 30% of the initial 203 rings because they had more than 40% of missing values. Within the 144 remaining rings, 15% of the pixels have missing values. Even though no impact was detected, we replaced the missing values by the azimuthally averaged dispersion before the modelling.

## 3.3. Generating toy models to study the impact of unresolved ordered velocity variations

As discussed in Sect. 3.1, BS is present in our data and is corrected in this study using observed line flux maps and velocity fields, which limits the accuracy of the correction. In order to study the impact of residual BS effects, we perform toy models for cases with barely resolved motions due to large-scale axisymmetric rotation of various strengths, with and without additional asymmetric velocity perturbations on local scales, under both isotropic and anisotropic hypotheses. We build toy models to produce data cubes, velocity, and $\sigma_{\text{los}}$ fields, and calculate the FFT of $\sigma_{\text{los}}$. This enables us to generate different configurations and identify the conditions for which systematic values in the distributions of $\sigma_k$ and/or $\phi_k$ could occur. Our toy models have no dynamical basis, yet they are very useful to assess the effects of BS, anisotropy, and streamings in $\sigma_{\text{los}}$. The full description of the toy models is presented in Appendix B.

We produced mock datacubes of 400 × 400 pixels (scale of 1", or ~50 pc at a typical distance of 10 Mpc), with 200 spectral elements with a 3 km s$^{-1}$ velocity sampling using 5 × 10$^6$ uniformly distributed points. To first order, the galaxy is assumed to be an axisymmetric rotating disc to which velocity perturbations can be added. Two rotation curves, one with a weak velocity gradient and a moderate velocity plateau, the other with a steep inner gradient quickly reaching a high velocity plateau, were used. Sharp planar velocity perturbations are produced by a bisymmetric spiral pattern, with five possible inner angles. No
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Maps and histograms of toy models with a steep rotation curve ($v_0 = 250$ km s$^{-1}$, $r_s = 20''$ and $\gamma = 2$) with a spiral perturbation. We show (top), from left to right, the density map, the velocity field, and the BS corrected dispersion maps in the isotropic and anisotropic cases, obtained with an inclination of 60° and with $\phi_{sp} = 135°$. The bottom panels show amplitude and phase histograms of orders $k = 2$ and $k = 4$. Blue and green histograms correspond to the cases with local uniform isotropy with large-scale asymmetries, and local uniform anisotropy with large-scale asymmetries, respectively.</figcaption></figure>

lopsidedness ($k = 1$ mode) is introduced for simplicity. Therefore, the perturbed models have intrinsically anisotropic velocities on large scales, as illustrated by the ellipsoid elongated in the radial dimension in Appendix B. We also synthesise velocity anisotropy in the toy models by modifying the shape of the velocity ellipsoid locally, that is by choosing a uniform radial bias $\sigma_{\theta} = 0.7\sigma_{R}$ with null covariance, with $\sigma_{R} = 8$ km s$^{-1}$ and $\sigma_{z} = 5$ km s$^{-1}$, in opposition to isotropic velocity distributions ($\sigma_{R} = \sigma_{\theta} = \sigma_{z} = 8$ km s$^{-1}$). The choice of $\sigma_{\theta} = 0.7\sigma_{R}$ comes from the radial bias seen of young stellar populations in the disc of the Milky Way (Gaia Collaboration 2023). This is an additional effect to the streaming-driven velocity anisotropy, and switching it off and on in the toy models enables to assess other specific sources of anisotropy not accounted for by the streaming perturbations in the planar velocity components.

For each particle, velocities are drawn randomly with a mean velocity and a velocity dispersion for each component. The mean in radial, azimuthal and vertical velocities are computed from the rotation curve and the perturbation. The velocities are then projected along the line-of-sight for three inclinations (45°, 60°, and 75°). For each inclination, cubes were created for 12 orientations of the spiral perturbation. Gaussian smoothing was then applied to mimic BS (8 pixels FWHM, corresponding to ~400 pc), before extracting moment maps. Beam-smearing corrections were applied as for observed data. Because of projection effects, a perturbation along the radial (azimuthal) direction is mainly seen along the minor (major) axis in the line-of-sight velocity fields and velocity dispersion maps (see Fig. B.2). FFTs of $\sigma_{\text{los}}^2$ were calculated for 35 independent rings, yielding profiles and distributions of $\sigma_k$ and $\phi_k$, as in Sect. 4 using the 36 models (3 inclinations, 12 orientations of the perturbation). Figure 3 shows a velocity perturbation (case where the velocity perturbation vector is orientated towards the closest position of the spiral) for the steep rotation curve. From the maps, we can see residual BS effects due to both large-scale rotation and to the unresolved perturbation. The velocity perturbation is more orientated azimuthally in the centre than in the outer parts. The effect of the radial bias is also clearly seen along the minor axis in both the velocity dispersion

map and in the histograms of the second order phases. Other cases are shown in Appendix B. An analysis of BS residuals induced by large-scale rotation is provided in Sect. 5.2.4 and the analysis of projection effects from asymmetric perturbations in Sect. 5.3.1.

# 4. Asymmetries in velocity dispersion maps

Coefficients of the FFTs were calculated for the 883 and 144 radial rings of THINGS and WHISP sources, respectively. Of the first 20 orders, the dominant asymmetries are those up to the fourth harmonic, in agreement with Chemin et al. (2020), so we limit the analysis to $k \le 4$ hereafter.

## 4.1. A case study: NGC 2841

The results of the Fourier analysis for all galaxies from the THINGS sample are presented in Appendix D. To illustrate examples of results, Fig. 4 shows the process and results for NGC 2841, from $\sigma_{\text{corr}}^2$ to the FFT and the individual FFT components. The normalised phase angles $\phi_k/T_k$, where $T_k = 2\pi/k$, are shown, with values of 0.5 and 1 corresponding to half and a full period. The normalised phases $\phi_1/T_1 = 0, 0.5, 1$ and $\phi_k/T_k = 0$ for other orders are aligned along the major axis of the galaxy, and $\phi_2/T_2 = 0.5$ along the minor axis. In NGC 2841, $\sigma_{\text{corr}}^2$ is high along a cross shaped pattern, with strong $k = 2$ and $k = 4$ perturbations and the $k = 2$ perturbation grows stronger with radius. The residual $\sigma_{\text{res}}^2$ is weak, showing that the first 4 orders reproduce the structure in $\sigma_{\text{corr}}^2$. The galaxy is lopsided at small radii (e.g. Baldwin et al. 1980) and this is detected by the $k = 1$ mode of the FFT. The phase angle of the bisymmetry does not vary much beyond $R \sim 170''$, at a value of ~0.7 times the period of the $k = 2$ asymmetry.

## 4.2. Census of asymmetries in H¹ random motions

### 4.2.1. General trends

Figure 5 presents the distributions of the amplitudes and normalised phase angles of the Fourier modes for our reference
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Maps and histograms of toy models with a steep rotation curve ($v_0 = 250$ km s$^{-1}$, $r_s = 20''$ and $\gamma = 2$) with a spiral perturbation. We show (top), from left to right, the density map, the velocity field, and the BS corrected dispersion maps in the isotropic and anisotropic cases, obtained with an inclination of 60° and with $\phi_{sp} = 135°$. The bottom panels show amplitude and phase histograms of orders $k = 2$ and $k = 4$. Blue and green histograms correspond to the cases with local uniform isotropy with large-scale asymmetries, and local uniform anisotropy with large-scale asymmetries, respectively.</figcaption></figure>

lopsidedness ($k = 1$ mode) is introduced for simplicity. Therefore, the perturbed models have intrinsically anisotropic velocities on large scales, as illustrated by the ellipsoid elongated in the radial dimension in Appendix B. We also synthesise velocity anisotropy in the toy models by modifying the shape of the velocity ellipsoid locally, that is by choosing a uniform radial bias $\sigma_{\theta} = 0.7\sigma_{R}$ with null covariance, with $\sigma_{R} = 8$ km s$^{-1}$ and $\sigma_{z} = 5$ km s$^{-1}$, in opposition to isotropic velocity distributions ($\sigma_{R} = \sigma_{\theta} = \sigma_{z} = 8$ km s$^{-1}$). The choice of $\sigma_{\theta} = 0.7\sigma_{R}$ comes from the radial bias seen of young stellar populations in the disc of the Milky Way (Gaia Collaboration 2023). This is an additional effect to the streaming-driven velocity anisotropy, and switching it off and on in the toy models enables to assess other specific sources of anisotropy not accounted for by the streaming perturbations in the planar velocity components.

For each particle, velocities are drawn randomly with a mean velocity and a velocity dispersion for each component. The mean in radial, azimuthal and vertical velocities are computed from the rotation curve and the perturbation. The velocities are then projected along the line-of-sight for three inclinations (45°, 60°, and 75°). For each inclination, cubes were created for 12 orientations of the spiral perturbation. Gaussian smoothing was then applied to mimic BS (8 pixels FWHM, corresponding to ~400 pc), before extracting moment maps. Beam-smearing corrections were applied as for observed data. Because of projection effects, a perturbation along the radial (azimuthal) direction is mainly seen along the minor (major) axis in the line-of-sight velocity fields and velocity dispersion maps (see Fig. B.2). FFTs of $\sigma_{\text{los}}^2$ were calculated for 35 independent rings, yielding profiles and distributions of $\sigma_k$ and $\phi_k$, as in Sect. 4 using the 36 models (3 inclinations, 12 orientations of the perturbation). Figure 3 shows a velocity perturbation (case where the velocity perturbation vector is orientated towards the closest position of the spiral) for the steep rotation curve. From the maps, we can see residual BS effects due to both large-scale rotation and to the unresolved perturbation. The velocity perturbation is more orientated azimuthally in the centre than in the outer parts. The effect of the radial bias is also clearly seen along the minor axis in both the velocity dispersion

map and in the histograms of the second order phases. Other cases are shown in Appendix B. An analysis of BS residuals induced by large-scale rotation is provided in Sect. 5.2.4 and the analysis of projection effects from asymmetric perturbations in Sect. 5.3.1.

# 4. Asymmetries in velocity dispersion maps

Coefficients of the FFTs were calculated for the 883 and 144 radial rings of THINGS and WHISP sources, respectively. Of the first 20 orders, the dominant asymmetries are those up to the fourth harmonic, in agreement with Chemin et al. (2020), so we limit the analysis to $k \le 4$ hereafter.

## 4.1. A case study: NGC 2841

The results of the Fourier analysis for all galaxies from the THINGS sample are presented in Appendix D. To illustrate examples of results, Fig. 4 shows the process and results for NGC 2841, from $\sigma_{\text{corr}}^2$ to the FFT and the individual FFT components. The normalised phase angles $\phi_k/T_k$, where $T_k = 2\pi/k$, are shown, with values of 0.5 and 1 corresponding to half and a full period. The normalised phases $\phi_1/T_1 = 0, 0.5, 1$ and $\phi_k/T_k = 0$ for other orders are aligned along the major axis of the galaxy, and $\phi_2/T_2 = 0.5$ along the minor axis. In NGC 2841, $\sigma_{\text{corr}}^2$ is high along a cross shaped pattern, with strong $k = 2$ and $k = 4$ perturbations and the $k = 2$ perturbation grows stronger with radius. The residual $\sigma_{\text{res}}^2$ is weak, showing that the first 4 orders reproduce the structure in $\sigma_{\text{corr}}^2$. The galaxy is lopsided at small radii (e.g. Baldwin et al. 1980) and this is detected by the $k = 1$ mode of the FFT. The phase angle of the bisymmetry does not vary much beyond $R \sim 170''$, at a value of ~0.7 times the period of the $k = 2$ asymmetry.

## 4.2. Census of asymmetries in H¹ random motions

### 4.2.1. General trends

Figure 5 presents the distributions of the amplitudes and normalised phase angles of the Fourier modes for our reference
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Examples of FFT results with the galaxy NGC 2841. *Top*: observed squared velocity dispersion and its modelling through FFT up to the order 4, and their residuals. *Middle panel*: individual orders of the FFT projected in the plane of the galaxy. *Bottom*: amplitudes (**left**) and phase angles (**right**) of FFT coefficients as a function of the radius.</figcaption></figure>

THINGS sample (blue histograms), and the WHISP sample appended to the THINGS sources (green histograms). Table 2 lists the mean, median and standard deviation of the amplitudes for the 883 rings of the THINGS sample. Globally, the distributions of amplitudes have comparable shapes. With an average of 11 km s⁻¹, the strongest Fourier mode is of second order, and the weakest is of third order (8 km s⁻¹). Most *k* = 2 amplitudes are between 2.5 and 7.5 km s⁻¹ with a few rings below 2.5 km s⁻¹, possibly due to noise.

For *k* = 1, we see in Fig. 5 (bottom) that the phase angle peaks at $\phi_1/T_1$ = 0.45 and 0.95 with dips at 0.15 and 0.75. The bisymmetric mode is maximal around $\phi_2/T_2$ = 0.65, and is minimum at 0. The *k* = 3 distribution is relatively flat with a maximum located at $\phi_3/T_3$ = 0.15. The *k* = 4 histogram is strongly peaked at $\phi_4/T_4$ = 0.5. The uncertainties are approximately the bin size of the distributions (0.1 $T_k$).

<table><caption>**Table 2.** Properties of Fourier amplitudes for the THINGS galaxies, as measured from 883 tilted rings.</caption><thead><tr><th>Order amplitude</th><th>Mean (km s<sup>-1</sup>)</th><th>Median (km s<sup>-1</sup>)</th><th>Standard deviation (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>σ<sub>1</sub></td><td>9.3</td><td>7.1</td><td>7.1</td></tr><tr><td>σ<sub>2</sub></td><td>10.9</td><td>8.0</td><td>7.7</td></tr><tr><td>σ<sub>3</sub></td><td>7.9</td><td>6.6</td><td>5.1</td></tr><tr><td>σ<sub>4</sub></td><td>8.8</td><td>7.0</td><td>6.1</td></tr></tbody></table>

The distributions of Fig. 5 may be biased by a few galaxies with more radial rings than others. To assess the impact of different numbers of rings, we measured new profiles of amplitudes
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Examples of FFT results with the galaxy NGC 2841. *Top*: observed squared velocity dispersion and its modelling through FFT up to the order 4, and their residuals. *Middle panel*: individual orders of the FFT projected in the plane of the galaxy. *Bottom*: amplitudes (**left**) and phase angles (**right**) of FFT coefficients as a function of the radius.</figcaption></figure>

THINGS sample (blue histograms), and the WHISP sample appended to the THINGS sources (green histograms). Table 2 lists the mean, median and standard deviation of the amplitudes for the 883 rings of the THINGS sample. Globally, the distributions of amplitudes have comparable shapes. With an average of 11 km s⁻¹, the strongest Fourier mode is of second order, and the weakest is of third order (8 km s⁻¹). Most *k* = 2 amplitudes are between 2.5 and 7.5 km s⁻¹ with a few rings below 2.5 km s⁻¹, possibly due to noise.

For *k* = 1, we see in Fig. 5 (bottom) that the phase angle peaks at $\phi_1/T_1 = 0.45$ and 0.95 with dips at 0.15 and 0.75. The bisymmetric mode is maximal around $\phi_2/T_2 = 0.65$, and is minimum at 0. The *k* = 3 distribution is relatively flat with a maximum located at $\phi_3/T_3 = 0.15$. The *k* = 4 histogram is strongly peaked at $\phi_4/T_4 = 0.5$. The uncertainties are approximately the bin size of the distributions ($0.1 T_k$).

<table><caption>**Table 2.** Properties of Fourier amplitudes for the THINGS galaxies, as measured from 883 tilted rings.</caption><thead><tr><th>Order amplitude</th><th>Mean (km s<sup>-1</sup>)</th><th>Median (km s<sup>-1</sup>)</th><th>Standard deviation (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>σ<sub>1</sub></td><td>9.3</td><td>7.1</td><td>7.1</td></tr><tr><td>σ<sub>2</sub></td><td>10.9</td><td>8.0</td><td>7.7</td></tr><tr><td>σ<sub>3</sub></td><td>7.9</td><td>6.6</td><td>5.1</td></tr><tr><td>σ<sub>4</sub></td><td>8.8</td><td>7.0</td><td>6.1</td></tr></tbody></table>

The distributions of Fig. 5 may be biased by a few galaxies with more radial rings than others. To assess the impact of different numbers of rings, we measured new profiles of amplitudes
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Results of the FFT measurements of the H I velocity dispersion maps. The histograms show the number of rings as a function of amplitude (top row) and normalised phase angle (bottom row) of the FFT harmonics of the THINGS sample (883 rings, in blue), the interpolated distribution with 20 rings per galaxy for THINGS (280 rings, in red, normalised to the maximum values of the blue histograms), and the WHISP sample (144 rings, in green). Phase angles are normalised to the period $T_k$ of each order.</figcaption></figure>

and phase angles using 20 equally spaced rings per galaxy as all galaxies but NGC 3627 have at least 20 independent rings. The initial profiles were interpolated at the 20 new radii. This yields new distributions built on 280 radii for 14 galaxies (NGC 3627 was excluded for this assessment). These histograms are shown as red lines in Fig. 5. The agreement between the blue and red histograms shows that the unequal number of radial bins has little impact on the results.

We can estimate the significance of the dips and peaks of the phase distributions. Assuming that the gas velocity ellipsoids are isotropic, the velocity dispersion maps should exhibit asymmetries randomly distributed over the plane of the sky. The 883 rings observed for the THINGS sample should thus yield, on average, $N = 88$ counts in each of the 10 bins of the histograms. If we quantify a confidence level as $\zeta = \sqrt{N}$, then the uncertainty is $\zeta \sim 9$. Keeping only the bins showing an excess of at least 3$\zeta$, the peaks observed at 0.45 and 0.95 for $k = 1$ are detected at a level of 6.3$\zeta$ and 3.3$\zeta$ (147 and 119 counts, respectively), 5.5$\zeta$ at 0.65 for $k = 2$ (140 counts), 4.8$\zeta$ at 0.15 for $k = 3$ (133 counts), and up to 11$\zeta$ for the bin at 0.45 of $k = 4$ (194 counts). Now keeping only the bins showing a deficit of at least 3$\zeta$, the minima observed at 0.15 and 0.75 for $k = 1$ correspond to a confidence level of 4.2$\zeta$ (49 counts), 3.7$\zeta$ at 0.95 for $k = 2$ (53 counts), and 4.5 to 5.3$\zeta$ at 0.05, 0.15, 0.25 and 0.95 for $k = 4$ (45, 42, 46 and 38 counts, respectively). Therefore, the systematic phase angles are not consistent with a random fluctuation of the orientations of asymmetries in the galaxies. We note that in the case of $k = 3$ no trend is seen around the unique bin exceeding 3$\zeta$. It may be that this peak is caused by chance due to the limited statistics, and not by a systematic effect, unlike the peaks or dips seen in the other asymmetries. With the interpolated profiles at 20 bins for each galaxy (red histograms of Fig. 5), we find less significant peaks and dips than with the 883 rings due to lower statistics, but the trends are preserved.

### 4.2.2. Correlations between the Fourier modes

The upper row of Fig. 6 compares the Fourier amplitudes as a function of the phase angle $\phi_2/T_2$ for the strongest orders

($k = 1, 2, 4$) for THINGS galaxies. The colour code is the number of tilted rings within bin widths of 0.05 for $\phi_2/T_2$, and 1.75 for the amplitude difference. The density of tilted rings is highest around a difference of 0 km s$^{-1}$ (in agreement with values given in Table 2), irrespective of the value of $\phi_2/T_2$, which implies correlated amplitudes. We measure a Pearson correlation coefficient of 0.8 ± 0.1 between the amplitudes of orders 2 and 4.

The bottom row of Fig. 6 compares the phase angle differences between orders $\Delta\phi_{m,n} = \phi_m - \phi_n$, and shows the differences within one period of the highest order, that is within $\pm T_m/2$, with $m > n$. The distributions of $\Delta\phi_{2,1}$ and $\Delta\phi_{4,2}$ are peaked and symmetric around zero, implying a correlation between the $k = 1$ and $k = 2$ modes and between $k = 2$ and $k = 4$. These correlations suggest that $k = 4$ is a harmonic of the second order perturbation.

### 4.3. Correlations with the H I flux density

Several processes may induce correlations or anti-correlations between gas velocity dispersion and gas density. A large velocity dispersion can arise from unresolved motions, such as insufficient spatial resolution or the presence of unresolved multiple peak profiles. These unresolved velocity gradients may result from various physical processes, including steep density gradients, gas compression in density waves, instabilities associated with spiral-like features, or starburst outflows. The faintest regions may also display a high velocity dispersion because the S/N affects the profile widths, noisier profiles appearing broader. Leaving aside the aforementioned observational caveats (low resolution, low S/N) and focussing on the resolved regions, correlations between dispersion and density can have different origins. Large-scale star formation is clearly linked to H I content but this is not true for small scales (e.g. Zhou et al. 2018). As long as the H I density is below the density necessary to gravitationally collapse and form molecular hydrogen, the H I gas clouds have no reason to cool and the velocity dispersion to decrease. In some galaxies (e.g. NGC 4214, Wilcots & Thurow 2001), the largest H$\alpha$ velocity dispersion is observed in the diffuse ionised gas regions which often has a
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Results of the FFT measurements of the H I velocity dispersion maps. The histograms show the number of rings as a function of amplitude (top row) and normalised phase angle (bottom row) of the FFT harmonics of the THINGS sample (883 rings, in blue), the interpolated distribution with 20 rings per galaxy for THINGS (280 rings, in red, normalised to the maximum values of the blue histograms), and the WHISP sample (144 rings, in green). Phase angles are normalised to the period $T_k$ of each order.</figcaption></figure>

and phase angles using 20 equally spaced rings per galaxy as all galaxies but NGC 3627 have at least 20 independent rings. The initial profiles were interpolated at the 20 new radii. This yields new distributions built on 280 radii for 14 galaxies (NGC 3627 was excluded for this assessment). These histograms are shown as red lines in Fig. 5. The agreement between the blue and red histograms shows that the unequal number of radial bins has little impact on the results.

We can estimate the significance of the dips and peaks of the phase distributions. Assuming that the gas velocity ellipsoids are isotropic, the velocity dispersion maps should exhibit asymmetries randomly distributed over the plane of the sky. The 883 rings observed for the THINGS sample should thus yield, on average, $N = 88$ counts in each of the 10 bins of the histograms. If we quantify a confidence level as $\zeta = \sqrt{N}$, then the uncertainty is $\zeta \sim 9$. Keeping only the bins showing an excess of at least 3$\zeta$, the peaks observed at 0.45 and 0.95 for $k = 1$ are detected at a level of 6.3$\zeta$ and 3.3$\zeta$ (147 and 119 counts, respectively), 5.5$\zeta$ at 0.65 for $k = 2$ (140 counts), 4.8$\zeta$ at 0.15 for $k = 3$ (133 counts), and up to 11$\zeta$ for the bin at 0.45 of $k = 4$ (194 counts). Now keeping only the bins showing a deficit of at least 3$\zeta$, the minima observed at 0.15 and 0.75 for $k = 1$ correspond to a confidence level of 4.2$\zeta$ (49 counts), 3.7$\zeta$ at 0.95 for $k = 2$ (53 counts), and 4.5 to 5.3$\zeta$ at 0.05, 0.15, 0.25 and 0.95 for $k = 4$ (45, 42, 46 and 38 counts, respectively). Therefore, the systematic phase angles are not consistent with a random fluctuation of the orientations of asymmetries in the galaxies. We note that in the case of $k = 3$ no trend is seen around the unique bin exceeding 3$\zeta$. It may be that this peak is caused by chance due to the limited statistics, and not by a systematic effect, unlike the peaks or dips seen in the other asymmetries. With the interpolated profiles at 20 bins for each galaxy (red histograms of Fig. 5), we find less significant peaks and dips than with the 883 rings due to lower statistics, but the trends are preserved.

### 4.2.2. Correlations between the Fourier modes

The upper row of Fig. 6 compares the Fourier amplitudes as a function of the phase angle $\phi_2/T_2$ for the strongest orders

($k = 1, 2, 4$) for THINGS galaxies. The colour code is the number of tilted rings within bin widths of 0.05 for $\phi_2/T_2$, and 1.75 for the amplitude difference. The density of tilted rings is highest around a difference of $0 \text{ km s}^{-1}$ (in agreement with values given in Table 2), irrespective of the value of $\phi_2/T_2$, which implies correlated amplitudes. We measure a Pearson correlation coefficient of 0.8 ± 0.1 between the amplitudes of orders 2 and 4.

The bottom row of Fig. 6 compares the phase angle differences between orders $\Delta\phi_{m,n} = \phi_m - \phi_n$, and shows the differences within one period of the highest order, that is within $\pm T_m/2$, with $m > n$. The distributions of $\Delta\phi_{2,1}$ and $\Delta\phi_{4,2}$ are peaked and symmetric around zero, implying a correlation between the $k = 1$ and $k = 2$ modes and between $k = 2$ and $k = 4$. These correlations suggest that $k = 4$ is a harmonic of the second order perturbation.

### 4.3. Correlations with the H I flux density

Several processes may induce correlations or anti-correlations between gas velocity dispersion and gas density. A large velocity dispersion can arise from unresolved motions, such as insufficient spatial resolution or the presence of unresolved multiple peak profiles. These unresolved velocity gradients may result from various physical processes, including steep density gradients, gas compression in density waves, instabilities associated with spiral-like features, or starburst outflows. The faintest regions may also display a high velocity dispersion because the S/N affects the profile widths, noisier profiles appearing broader. Leaving aside the aforementioned observational caveats (low resolution, low S/N) and focussing on the resolved regions, correlations between dispersion and density can have different origins. Large-scale star formation is clearly linked to H I content but this is not true for small scales (e.g. Zhou et al. 2018). As long as the H I density is below the density necessary to gravitationally collapse and form molecular hydrogen, the H I gas clouds have no reason to cool and the velocity dispersion to decrease. In some galaxies (e.g. NGC 4214, Wilcots & Thurow 2001), the largest H$\alpha$ velocity dispersion is observed in the diffuse ionised gas regions which often has a
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Comparisons of amplitudes and phase angles between the Fourier $k = 1, 2$ and 4 orders, shown as differences of amplitudes versus $\phi_2/T_2$ (upper row) and histograms of phase angles difference normalised to $\pi$ (bottom row). The comparisons between the orders $k = 1$ and $k = 2$, $k = 1$ and $k = 4$, and $k = 2$ and $k = 4$ are shown in the left, middle, and right columns, respectively. In the upper row, the amplitude difference is colour coded with the logarithmic number of the tilted rings.</figcaption></figure>

low H I column density. Around massive star forming regions, the large velocity gradients might be explained by outflows. In low H I column density regions far from massive stars (e.g. NGC 2366, Hunter et al. 2001), in the absence of a local source, the velocity width might be induced by long-range turbulent pressure.

Using NGC 2841 as an example, Fig. 1 shows that the gas density ($\Sigma_{\text{los}}$) is principally distributed in bright rings and outer spiral arms, whereas larger velocity dispersion ($\sigma_{\text{los}}$) is along a centred cross-shaped structure. In other words, the gas density is either anti-correlated or correlated with velocity dispersion, sometimes fainter in large velocity dispersion regions (e.g. in the galaxy centre), brighter at low velocity dispersion (e.g. in the NE and SW quadrants), but also brighter at larger velocity dispersion (e.g. NW and SE quadrants). Across the sample, similar features are observed, and velocity dispersion seems correlated to spiral arms. Therefore, due to the wide diversity of processes, a galaxy-by-galaxy or a pixel-by-pixel comparison between the flux density and the velocity dispersion is beyond the scope of this paper. We rather aim at identifying global trends with respect to our analysis of the velocity dispersion.

The correlation between gas density and velocity dispersion has therefore been investigated for the whole sample using FFTs of the density maps, in the same way as for the velocity dispersion. The phase angle differences $\Delta\phi_k \equiv \phi_k(\Sigma) - \phi_k(\sigma)$ are shown in Fig. 7. As in NGC 2841, we observe that all orders except $k = 4$ exhibit correlated (peaks) and uncorrelated (troughs) phases. An anti-correlation is observed for $k = 1$ (dip at $\Delta\phi_k \sim 0$), and a correlation for $k = 2$ (peak at $\Delta\phi_k \sim 0$). These findings suggest that the asymmetries observed in the velocity dispersion are

related to those of the gas distribution that are mainly induced by spiral arms, bars, warps and lopsidedness.

### 4.4. FFT results versus galaxy properties

Our sample is somewhat biased by the fact that the smallest galaxies, which are usually the latest type and faintest ones, are also the closest<fnref n="4" />. To look for an effect, we divided the sample into two sub-samples of seven galaxies each around the median value of each parameter: morphological type, optical radius, absolute magnitude, and metallicity. Class I (II) is the sub-sample with the more distant (nearest), or the brightest (faintest), or the earliest (latest) type, or the largest (smallest) galaxies.

In Fig. 8, we show histograms of $\sigma_2$, $\phi_2$, and $\phi_4$ for classes I (red) and II (blue) of absolute magnitude, which is the parameter providing the largest difference between classes I and II. As expected, $\sigma_2$ is lower for small galaxies. For $k = 2$, the incidence of $\phi_2/T_2 \sim 0.5$ is low (high, respectively) for class I galaxies (class II). $\phi_2/T_2$ in class II galaxies agrees with the observations of M 33 (Chemin et al. 2020), which is a class II galaxy. For orders $k = 1$ and 3, the phase does not depend significantly on the class. $\phi_4/T_4$ has a bigger peak for class I than for class II but the difference is smaller for other parameters than the absolute magnitudes. The $\phi_4/T_4$ distribution is almost flat for faint galaxies and peaked for bright ones. This could be due to residual BS (see Sect. 5.2.4) or because higher gas density contrasts

<footnote n="4">⁴ It is a well-known observational merit factor bias which consists in ‘filling the field of view’ of the instrument as much as possible.</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Comparisons of amplitudes and phase angles between the Fourier $k = 1, 2$ and 4 orders, shown as differences of amplitudes versus $\phi_2/T_2$ (upper row) and histograms of phase angles difference normalised to $\pi$ (bottom row). The comparisons between the orders $k = 1$ and $k = 2$, $k = 1$ and $k = 4$, and $k = 2$ and $k = 4$ are shown in the left, middle, and right columns, respectively. In the upper row, the amplitude difference is colour coded with the logarithmic number of the tilted rings.</figcaption></figure>

low H I column density. Around massive star forming regions, the large velocity gradients might be explained by outflows. In low H I column density regions far from massive stars (e.g. NGC 2366, Hunter et al. 2001), in the absence of a local source, the velocity width might be induced by long-range turbulent pressure.

Using NGC 2841 as an example, Fig. 1 shows that the gas density ($\Sigma_{\text{los}}$) is principally distributed in bright rings and outer spiral arms, whereas larger velocity dispersion ($\sigma_{\text{los}}$) is along a centred cross-shaped structure. In other words, the gas density is either anti-correlated or correlated with velocity dispersion, sometimes fainter in large velocity dispersion regions (e.g. in the galaxy centre), brighter at low velocity dispersion (e.g. in the NE and SW quadrants), but also brighter at larger velocity dispersion (e.g. NW and SE quadrants). Across the sample, similar features are observed, and velocity dispersion seems correlated to spiral arms. Therefore, due to the wide diversity of processes, a galaxy-by-galaxy or a pixel-by-pixel comparison between the flux density and the velocity dispersion is beyond the scope of this paper. We rather aim at identifying global trends with respect to our analysis of the velocity dispersion.

The correlation between gas density and velocity dispersion has therefore been investigated for the whole sample using FFTs of the density maps, in the same way as for the velocity dispersion. The phase angle differences $\Delta\phi_k \equiv \phi_k(\Sigma) - \phi_k(\sigma)$ are shown in Fig. 7. As in NGC 2841, we observe that all orders except $k = 4$ exhibit correlated (peaks) and uncorrelated (troughs) phases. An anti-correlation is observed for $k = 1$ (dip at $\Delta\phi_k \sim 0$), and a correlation for $k = 2$ (peak at $\Delta\phi_k \sim 0$). These findings suggest that the asymmetries observed in the velocity dispersion are

related to those of the gas distribution that are mainly induced by spiral arms, bars, warps and lopsidedness.

### 4.4. FFT results versus galaxy properties

Our sample is somewhat biased by the fact that the smallest galaxies, which are usually the latest type and faintest ones, are also the closest<fnref n="4" />. To look for an effect, we divided the sample into two sub-samples of seven galaxies each around the median value of each parameter: morphological type, optical radius, absolute magnitude, and metallicity. Class I (II) is the sub-sample with the more distant (nearest), or the brightest (faintest), or the earliest (latest) type, or the largest (smallest) galaxies.

In Fig. 8, we show histograms of $\sigma_2$, $\phi_2$, and $\phi_4$ for classes I (red) and II (blue) of absolute magnitude, which is the parameter providing the largest difference between classes I and II. As expected, $\sigma_2$ is lower for small galaxies. For $k = 2$, the incidence of $\phi_2/T_2 \sim 0.5$ is low (high, respectively) for class I galaxies (class II). $\phi_2/T_2$ in class II galaxies agrees with the observations of M 33 (Chemin et al. 2020), which is a class II galaxy. For orders $k = 1$ and 3, the phase does not depend significantly on the class. $\phi_4/T_4$ has a bigger peak for class I than for class II but the difference is smaller for other parameters than the absolute magnitudes. The $\phi_4/T_4$ distribution is almost flat for faint galaxies and peaked for bright ones. This could be due to residual BS (see Sect. 5.2.4) or because higher gas density contrasts

<footnote n="4">⁴ It is a well-known observational merit factor bias which consists in ‘filling the field of view’ of the instrument as much as possible.</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Histograms of phase angles (normalised by $\pi$) difference between gas density and velocity dispersion $\phi_k(\Sigma) - \phi_k(\sigma)$ for $k = 1, 2$ and 4 (left, middle, and right panels, respectively).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Histograms of the velocity dispersion amplitudes $\sigma_k$ ($k = 2$, left panel) and of the normalised phase angles $\phi_k/T_k$ ($k = 2$, middle panel and 4, right panel) for the 7 galaxies brighter and fainter than $M = -20.7$ mag (red and blue colours, respectively) for the THINGS sample. The sub-samples are the same as when splitting between high (red) and low (blue) projected maximum velocity (see Sect. 5.2.4). The grey histograms are the total distributions from Fig. 5.</figcaption></figure>

are expected in earlier type spirals, inducing larger perturbations
(see Sect. 5.3.1).

**5. Discussion**

Our FFT analysis of the velocity dispersion of the H I gas shows that common values for $\phi_k/T_k$ are (1) $\phi_1/T_1 \approx 0.45$ or 0.95, (2) $\phi_2/T_2 \approx 0.35$ and 0.65 and (3) $\phi_4/T_4 \approx 0.45$, with no major structure in the $k=3$ mode. These peaks correspond to major axis of the 1st order asymmetry and half period of the 4th order asymmetry. Interestingly, finding $k=2$ asymmetries projected preferentially near the major axis is significantly less likely.

In the absence of instrumental effects such as BS at any scale, and assuming gas is isotropic, we expect the principal axes of the galaxies to be randomly orientated with respect to the spiral pattern or bar, so the distributions of $\phi_k/T_k$ should be uniform. In this section, we discuss the origin of the alignment of phases around particular angles.

**5.1. Systematic effects in the FFT measurements and the tilted-ring models**

**5.1.1. Significance of FFT coefficients**

In Appendix A, we describe the method used to infer the accu-
racy $d(\sigma_k^2)$ of a measured amplitude $\sigma_k^2$ (Eq. (A.5)), and its sig-
nificance $s_k = \sigma_k^2/d(\sigma_k^2)$. For this exercise, we assume a rea-
sonable uncertainty $d\sigma = 1 \text{ km s}^{-1}$ on the standard deviation of
the line spread function. By keeping only the rings for which

the significance is greater than 3, we find that amongst the 883
rings of the THINGS sample, 482 exceed this threshold for the
order $k = 1$, 587 for $k = 2$, 431 for $k = 3$, and 536 for
$k = 4$.

We further divided the rings into high and low significance
samples using the median significance of each order. The median
values of the significance $s_k$ are 3.3, 5.2, 2.9, and 3.9 for $k =$
1, 2, 3 and 4, respectively. The resulting histograms are given
in Fig. 9, with the low(high)-significance samples shown in red
(blue, respectively). A first result is that apart larger tails in the
distributions of high-significance samples, the amplitudes are
not much affected by the division into two sub-samples. Since
some high amplitudes rings are located in the innermost part
of galaxies where the number of pixels is much lower than
in the outskirt, we observed that some high amplitudes rings
(red tail) are less significant than some of the low amplitude
ones (blue distribution at low amplitudes). On the other hand,
the low-significance samples show essentially random phases,
except for order $k = 4$, while the high-significance show the
trends identified before ($\phi_1/T_1 \approx 0.45$, $\phi_2/T_2 \approx 0.65$ and
$\phi_4/T_4 \approx 0.45$) but more clearly. In other words, the rings of
lower significance are consistent with having asymmetries ran-
domly distributed in the sky plane (except for $k = 4$), while
those of higher significance seem to be more representative of
the systematically orientated asymmetries in the velocity dis-
persion maps. Moreover, the fact that both sub-samples have
$\phi_4$ centrally peaked around the half period is a hint that low
amplitude residuals due to BS still affect the data, despite the
correction applied before deriving the FFTs. We address this in
Sect. 5.2.
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Histograms of phase angles (normalised by $\pi$) difference between gas density and velocity dispersion $\phi_k(\Sigma) - \phi_k(\sigma)$ for $k = 1, 2$ and 4 (**left**, *middle*, and *right panels*, respectively).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Histograms of the velocity dispersion amplitudes $\sigma_k$ ($k = 2$, **left panel**) and of the normalised phase angles $\phi_k/T_k$ ($k = 2$, **middle panel** and 4, **right panel**) for the 7 galaxies brighter and fainter than $M = -20.7$ mag (red and blue colours, respectively) for the THINGS sample. The sub-samples are the same as when splitting between high (red) and low (blue) projected maximum velocity (see Sect. 5.2.4). The grey histograms are the total distributions from Fig. 5.</figcaption></figure>

are expected in earlier type spirals, inducing larger perturbations
(see Sect. 5.3.1).

**5. Discussion**

Our FFT analysis of the velocity dispersion of the H I gas shows that common values for $\phi_k/T_k$ are (1) $\phi_1/T_1 \approx 0.45$ or 0.95, (2) $\phi_2/T_2 \approx 0.35$ and 0.65 and (3) $\phi_4/T_4 \approx 0.45$, with no major structure in the $k=3$ mode. These peaks correspond to major axis of the 1st order asymmetry and half period of the 4th order asymmetry. Interestingly, finding $k=2$ asymmetries projected preferentially near the major axis is significantly less likely.

In the absence of instrumental effects such as BS at any scale, and assuming gas is isotropic, we expect the principal axes of the galaxies to be randomly orientated with respect to the spiral pattern or bar, so the distributions of $\phi_k/T_k$ should be uniform. In this section, we discuss the origin of the alignment of phases around particular angles.

**5.1. Systematic effects in the FFT measurements and the tilted-ring models**

**5.1.1. Significance of FFT coefficients**

In Appendix A, we describe the method used to infer the accu-
racy $d(\sigma_k^2)$ of a measured amplitude $\sigma_k^2$ (Eq. (A.5)), and its sig-
nificance $s_k = \sigma_k^2/d(\sigma_k^2)$. For this exercise, we assume a rea-
sonable uncertainty $d\sigma = 1 \text{ km s}^{-1}$ on the standard deviation of
the line spread function. By keeping only the rings for which

the significance is greater than 3, we find that amongst the 883
rings of the THINGS sample, 482 exceed this threshold for the
order $k = 1$, 587 for $k = 2$, 431 for $k = 3$, and 536 for
$k = 4$.

We further divided the rings into high and low significance
samples using the median significance of each order. The median
values of the significance $s_k$ are 3.3, 5.2, 2.9, and 3.9 for $k =$
1, 2, 3 and 4, respectively. The resulting histograms are given
in Fig. 9, with the low(high)-significance samples shown in red
(blue, respectively). A first result is that apart larger tails in the
distributions of high-significance samples, the amplitudes are
not much affected by the division into two sub-samples. Since
some high amplitudes rings are located in the innermost part
of galaxies where the number of pixels is much lower than
in the outskirt, we observed that some high amplitudes rings
(red tail) are less significant than some of the low amplitude
ones (blue distribution at low amplitudes). On the other hand,
the low-significance samples show essentially random phases,
except for order $k = 4$, while the high-significance show the
trends identified before ($\phi_1/T_1 \approx 0.45$, $\phi_2/T_2 \approx 0.65$ and
$\phi_4/T_4 \approx 0.45$) but more clearly. In other words, the rings of
lower significance are consistent with having asymmetries ran-
domly distributed in the sky plane (except for $k = 4$), while
those of higher significance seem to be more representative of
the systematically orientated asymmetries in the velocity dis-
persion maps. Moreover, the fact that both sub-samples have
$\phi_4$ centrally peaked around the half period is a hint that low
amplitude residuals due to BS still affect the data, despite the
correction applied before deriving the FFTs. We address this in
Sect. 5.2.
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Histograms of amplitudes (upper row) and normalised phase angles (bottom row) resulting from the FFT analysis for the THINGS sample when splitting the sample in two, using the median significance of each order (see Sect. 5.1.1). The blue and red distributions are for the rings with a significance respectively greater and smaller than the median value.</figcaption></figure>

### 5.1.2. Deviations in the tilted-ring model projection parameters

We now test the robustness of our results with respect to variations of kinematic centre, inclination and position angle. On the one hand, we changed the position of the kinematic centre of the rings by adding a random value selected from a Gaussian distribution of standard deviation of 1" and 3", respectively for the right ascension and declination, as expected from their average uncertainties (Trachternach et al. 2008). On the other hand, we ignored warps by using the mean inclination and position angle for each galaxy, as given in de Blok et al. (2008).

The new distributions of phases are then compared to those of Fig. 5, by averaging the positive differences of histogram heights of all bins. Overall, the trends described in Sect. 4 are still present. Varying the kinematic centre leads to differences of 10%, 7.2%, 5.7% and 7.7%, for *k* = 1, 2, 3 and 4, respectively. Not taking into account the disc warps has little impact as well, although the variations are more important, of 8%, 15.6%, 10% and 22.8% for *k* = 1, 2, 3 and 4, respectively. We can thus exclude the possibility that the systematic phase angles of the asymmetries in the HI velocity dispersion are caused by incorrect kinematic parameters.

## 5.2. Systematic signatures associated to BS

Beam smearing might provide residuals in the velocity dispersion beyond our correction process if locally the velocity gradient is not resolved as shown from Eq. (4), or if the radio beam is not well represented by the 2D Gaussian model described in Sect. 3.1. An incomplete BS correction could contribute to or create velocity dispersion asymmetries.

### 5.2.1. Properties of the BS pattern

To numerically explore the effect of BS, we first created a mock galaxy having both exponential flux and rotation curve models. In order to work with the same amplitude of the projected velocity field and thus with a similar BS induced velocity dispersion at zeroth order, we use a maximum projected velocity

of 100 km s⁻¹ independent from the inclination. The Gaussian beam width is set to 10" × 15", and in order to sample correctly the exponential disc, we adopt a pixel size of 1" and a scale-length of 50". Figure 10 shows (left) the BS induced velocity dispersion maps at 30° (top) and 70° (bottom) inclinations. The central and right panels show the amplitudes and phase angles of the FFT decompositions. In our mock galaxy and PSF models, the BS pattern can be described using only orders *k* = 2 and *k* = 4. The ratio between the amplitude of these two orders vary with the inclination, with the *k* = 4 strength being predominant at large inclination. This latter is maximum around the radius where the maximum rotation curve is reached. The phases φ₂ and φ₄ are rather constant with radius, with their value depending on inclination. At fixed inclination, this constant value is related to the position angle of the galactic disc relative to that of the beam. In case both position angles match, or if the beam is circular, we have exactly φ₂ = 0.5 *T*₂ and φ₄ = 0.5 *T*₄. For an inclination of 70°, φ₂ and φ₄ are always in the range of 0.4 *T*₂−0.6 *T*₂ and 0.4 *T*₄−0.6 *T*₄ respectively, with φ₄ being close to its half period for most position angles.

Figure 11 shows that asymmetries in σbs are centrally peaked and lower than the observed asymmetries. The phase angles of σbs are similar to those seen in the galaxies. Subtle differences are observed, like the peak of probability at φ₂/T₂ ~ 0.5 not observed in galaxies, or the stronger and narrower peak at φ₄/T₄ = 0.5. Finally, we also compared the results of the FFTs obtained with and without BS correction, and found negligible differences. These findings indicate either that our data are sparsely affected by BS, or that the BS correction is underestimated. In Sects. 5.2.3 and 5.2.4, we further investigate the impact of using higher resolution data to infer BS correction and the amplitude of BS residuals respectively.

### 5.2.2. Uncertainty on the beam shape

The observed dirty beam of interferometric observations is complex. In this section, we measure σbs using the observed dirty beam for 12 of the galaxies from THINGS. Maps of the dirty beam are mainly composed by a combination of two Gaussian functions, a narrow one to describe the PSF core, similar to the
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Histograms of amplitudes (upper row) and normalised phase angles (bottom row) resulting from the FFT analysis for the THINGS sample when splitting the sample in two, using the median significance of each order (see Sect. 5.1.1). The blue and red distributions are for the rings with a significance respectively greater and smaller than the median value.</figcaption></figure>

### 5.1.2. Deviations in the tilted-ring model projection parameters

We now test the robustness of our results with respect to variations of kinematic centre, inclination and position angle. On the one hand, we changed the position of the kinematic centre of the rings by adding a random value selected from a Gaussian distribution of standard deviation of 1" and 3", respectively for the right ascension and declination, as expected from their average uncertainties (Trachternach et al. 2008). On the other hand, we ignored warps by using the mean inclination and position angle for each galaxy, as given in de Blok et al. (2008).

The new distributions of phases are then compared to those of Fig. 5, by averaging the positive differences of histogram heights of all bins. Overall, the trends described in Sect. 4 are still present. Varying the kinematic centre leads to differences of 10%, 7.2%, 5.7% and 7.7%, for *k* = 1, 2, 3 and 4, respectively. Not taking into account the disc warps has little impact as well, although the variations are more important, of 8%, 15.6%, 10% and 22.8% for *k* = 1, 2, 3 and 4, respectively. We can thus exclude the possibility that the systematic phase angles of the asymmetries in the HI velocity dispersion are caused by incorrect kinematic parameters.

## 5.2. Systematic signatures associated to BS

Beam smearing might provide residuals in the velocity dispersion beyond our correction process if locally the velocity gradient is not resolved as shown from Eq. (4), or if the radio beam is not well represented by the 2D Gaussian model described in Sect. 3.1. An incomplete BS correction could contribute to or create velocity dispersion asymmetries.

### 5.2.1. Properties of the BS pattern

To numerically explore the effect of BS, we first created a mock galaxy having both exponential flux and rotation curve models. In order to work with the same amplitude of the projected velocity field and thus with a similar BS induced velocity dispersion at zeroth order, we use a maximum projected velocity

of 100 km s⁻¹ independent from the inclination. The Gaussian beam width is set to 10" × 15", and in order to sample correctly the exponential disc, we adopt a pixel size of 1" and a scale-length of 50". Figure 10 shows (left) the BS induced velocity dispersion maps at 30° (top) and 70° (bottom) inclinations. The central and right panels show the amplitudes and phase angles of the FFT decompositions. In our mock galaxy and PSF models, the BS pattern can be described using only orders *k* = 2 and *k* = 4. The ratio between the amplitude of these two orders vary with the inclination, with the *k* = 4 strength being predominant at large inclination. This latter is maximum around the radius where the maximum rotation curve is reached. The phases φ₂ and φ₄ are rather constant with radius, with their value depending on inclination. At fixed inclination, this constant value is related to the position angle of the galactic disc relative to that of the beam. In case both position angles match, or if the beam is circular, we have exactly φ₂ = 0.5*T₂ and φ₄ = 0.5*T₄. For an inclination of 70°, φ₂ and φ₄ are always in the range of 0.4*T₂−0.6*T₂ and 0.4*T₄−0.6*T₄ respectively, with φ₄ being close to its half period for most position angles.

Figure 11 shows that asymmetries in σbs are centrally peaked and lower than the observed asymmetries. The phase angles of σbs are similar to those seen in the galaxies. Subtle differences are observed, like the peak of probability at φ₂/T₂ ~ 0.5 not observed in galaxies, or the stronger and narrower peak at φ₄/T₄ = 0.5. Finally, we also compared the results of the FFTs obtained with and without BS correction, and found negligible differences. These findings indicate either that our data are sparsely affected by BS, or that the BS correction is underestimated. In Sects. 5.2.3 and 5.2.4, we further investigate the impact of using higher resolution data to infer BS correction and the amplitude of BS residuals respectively.

### 5.2.2. Uncertainty on the beam shape

The observed dirty beam of interferometric observations is complex. In this section, we measure σbs using the observed dirty beam for 12 of the galaxies from THINGS. Maps of the dirty beam are mainly composed by a combination of two Gaussian functions, a narrow one to describe the PSF core, similar to the
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Results of the FFTs applied to the BS dispersion maps $\sigma_{bs}$, as built from disc models projected at 30° (top row) and 70° (bottom row). Left panel: map of $\sigma_{bs}$. The synthesised beam adopted in the toy model is shown on the bottom left corner. Middle panel: amplitudes of the Fourier modes for orders $k = 0$ to $k = 4$. Right panel: phase angles of the Fourier modes. A horizontal dashed line shows the half period of the Fourier modes. The disc scale-length is shown as a white contour (left) and as a vertical dashed line (middle and right).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Histograms of amplitudes (upper row) and normalised phase angles (bottom row) resulting from the FFT analysis for the THINGS sample (in blue) and for the BS dispersion maps $\sigma_{bs}$ (in red).</figcaption></figure>

one of our model, and a second one with larger full width at half maximum to describe the PSF wings. The extent of dirty beam maps is often larger than those of the galaxies but the beam power at large angular distance from its centre is very low. We apply a circular mask to the dirty beam map to deal mostly with

positive PSF values, and get rid of small variations at large radius from the beam centre that could induce spurious effects in our analysis. We vary the radius of the masks from 5 to 50 pixels, in order to be consistent with our 2D Gaussian model which has an extent of 5 to 10 pixels within the THINGS sample. Using a
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Results of the FFTs applied to the BS dispersion maps $\sigma_{bs}$, as built from disc models projected at 30° (top row) and 70° (bottom row). Left panel: map of $\sigma_{bs}$. The synthesised beam adopted in the toy model is shown on the bottom left corner. Middle panel: amplitudes of the Fourier modes for orders $k = 0$ to $k = 4$. Right panel: phase angles of the Fourier modes. A horizontal dashed line shows the half period of the Fourier modes. The disc scale-length is shown as a white contour (left) and as a vertical dashed line (middle and right).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Histograms of amplitudes (upper row) and normalised phase angles (bottom row) resulting from the FFT analysis for the THINGS sample (in blue) and for the BS dispersion maps $\sigma_{bs}$ (in red).</figcaption></figure>

one of our model, and a second one with larger full width at half maximum to describe the PSF wings. The extent of dirty beam maps is often larger than those of the galaxies but the beam power at large angular distance from its centre is very low. We apply a circular mask to the dirty beam map to deal mostly with

positive PSF values, and get rid of small variations at large radius from the beam centre that could induce spurious effects in our analysis. We vary the radius of the masks from 5 to 50 pixels, in order to be consistent with our 2D Gaussian model which has an extent of 5 to 10 pixels within the THINGS sample. Using a
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Comparison of the BS models using optical and cm data for the galaxy NGC 3521. **Left panel:** Hα data based BS model. **Middle panel:** H1 data based BS model, trimmed at the same extent as the optical model. **Right panel:** ratio of the BS models, defined as $BS(Hα)/BS(H1)$.</figcaption></figure>

dirty beam map truncated at the same size and scale as our 2D Gaussian model, the resulting $\sigma_{bs}$ maps are very comparable to those obtained with the Gaussian beam. For example, within a cut of five pixels around the beam centroid, the amplitude of $\sigma_{bs}$ is ~10% lower on average than for the Gaussian shape throughout the THINGS sample. However, when the spatial extent of the dirty beam is larger, the patterns in $\sigma_{bs}$ are more prominent, leading to larger amplitudes. For instance, with a mask cut at a radius of 20 pixels, the amplitudes are twice larger than our model, on average. As for the phase angles, the minor difference between the two beam shapes is that the peaks are narrower and slightly larger for $k=2$ and 4 in the case dirty beam.

### 5.2.3. Kinematics at higher angular resolution

Daigle et al. (2006) and Dicaire et al. (2008) published Hα data cubes for respectively 28 and 37 galaxies from the SINGS sample (Kennicutt et al. 2003). Those observations were led with a Fabry-Perot interferometer at the 1.6 m telescope of the Observatoire du Mont Megantic (OMM, Canada), and the 3.6 m telescopes at the Canada-France-Hawaii Telescope and the European Southern Observatory (ESO, La Silla, Chile). In our sample, all galaxies were observed at OMM except NGC 3621 and NGC 7793, observed at ESO. In this section, we benefit from the higher angular resolution of the Hα velocity fields to estimate, and correct from, the BS effects on H1 THINGS velocity fields. The angular resolution of the Hα data is seeing limited, ~1" at ESO and $\le 3''$ at OMM, which represents a gain of ~2 to 13 with respect to THINGS data. The underlying hypothesis is that the Hα kinematics have a similar behaviour as the H1 one. Nevertheless, the flux distribution of the ionised gas is more peaked due to the star forming regions. To model the Hα BS map, we used the Hα velocity field weighted by the flux distribution of the neutral gas. Since the kinematical maps available from Daigle et al. (2006) and Dicaire et al. (2008) were obtained using adaptative smoothing, we compute new maps from the sky-removed data cubes. We first apply a 2-pixels Gaussian smoothing to increase the Hα S/N and then extract kinematical maps using the same barycentre method used in the two original papers. We also reproject the THINGS flux density maps on the same grid as Hα maps in order to match the positions and pixel size of both datasets. To compute the BS maps, we apply the same methodology as defined in Sect. 3.1, using the reprojected H1 flux distributions, the 2-pixels Gaussian smoothed Hα velocity field and the single 2D Gaussian PSF to describe the new H1 angular resolution.

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Pixel-by-pixel comparison of the BS contribution using Hα SINGS and H1 THINGS data. The y-axis shows the ratio of the BS amplitudes $Hα/H1$ and in the x-axis is the amplitude of the H1 model.</figcaption></figure>

In Fig. 12, we show the resulting $Hα \sigma_{bs}$ map (left panel), compared to the H1 $\sigma_{bs}$ map (middle panel), as well as the ratio of the two maps, for an example galaxy (NGC 3521). The BS pattern is similar in both cases, with the typical cross-shape in the centre. While the BS correction inferred from H1 data is less peaked in the centre than for $Hα$, the latter also decreases more slowly with radius. This explains the large values in the normalised ratio found far from the centre. Beam-smearing amplitudes are higher when modelling with the ionised gas. As the amplitude may be an important concern in our modelling, we study galaxy by galaxy the amplitude ratio for the two different cases.

In Fig. 13, a comparison of the BS contributions is made for the THINGS sample. For each galaxy, a 2D histogram with the amplitude of the H1 BS correction on the x-axis and the ratio of the BS corrections in $Hα$ to H1 on the y-axis, is made. Then, the histograms were normalised by the number of pixels in the maps, to avoid being biased by objects with a larger number of pixels. Figure 13 is then obtained by summing all the 2D histograms. It shows that the ratio is >1, which is not surprising given that the Hα BS benefits from a higher angular resolution. On average, the H1 BS is 2.3 times smaller than the Hα BS. We also observe a trend at low H1 BS with high ratio induced either by noise or by discontinuities between regions in the $Hα$ data. Indeed, a continuously increasing ratio with decreasing H1 BS can be modelled by assuming that the Hα BS is equal to the quadratic sum of
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Comparison of the BS models using optical and cm data for the galaxy NGC 3521. **Left panel:** Hα data based BS model. **Middle panel:** H1 data based BS model, trimmed at the same extent as the optical model. **Right panel:** ratio of the BS models, defined as $BS(Hα)/BS(H1)$.</figcaption></figure>

dirty beam map truncated at the same size and scale as our 2D Gaussian model, the resulting $\sigma_{bs}$ maps are very comparable to those obtained with the Gaussian beam. For example, within a cut of five pixels around the beam centroid, the amplitude of $\sigma_{bs}$ is ~10% lower on average than for the Gaussian shape throughout the THINGS sample. However, when the spatial extent of the dirty beam is larger, the patterns in $\sigma_{bs}$ are more prominent, leading to larger amplitudes. For instance, with a mask cut at a radius of 20 pixels, the amplitudes are twice larger than our model, on average. As for the phase angles, the minor difference between the two beam shapes is that the peaks are narrower and slightly larger for $k=2$ and 4 in the case dirty beam.

### 5.2.3. Kinematics at higher angular resolution

Daigle et al. (2006) and Dicaire et al. (2008) published Hα data cubes for respectively 28 and 37 galaxies from the SINGS sample (Kennicutt et al. 2003). Those observations were led with a Fabry-Perot interferometer at the 1.6 m telescope of the Observatoire du Mont Megantic (OMM, Canada), and the 3.6 m telescopes at the Canada-France-Hawaii Telescope and the European Southern Observatory (ESO, La Silla, Chile). In our sample, all galaxies were observed at OMM except NGC 3621 and NGC 7793, observed at ESO. In this section, we benefit from the higher angular resolution of the Hα velocity fields to estimate, and correct from, the BS effects on H1 THINGS velocity fields. The angular resolution of the Hα data is seeing limited, ~1" at ESO and $\le 3''$ at OMM, which represents a gain of ~2 to 13 with respect to THINGS data. The underlying hypothesis is that the Hα kinematics have a similar behaviour as the H1 one. Nevertheless, the flux distribution of the ionised gas is more peaked due to the star forming regions. To model the Hα BS map, we used the Hα velocity field weighted by the flux distribution of the neutral gas. Since the kinematical maps available from Daigle et al. (2006) and Dicaire et al. (2008) were obtained using adaptative smoothing, we compute new maps from the sky-removed data cubes. We first apply a 2-pixels Gaussian smoothing to increase the Hα S/N and then extract kinematical maps using the same barycentre method used in the two original papers. We also reproject the THINGS flux density maps on the same grid as Hα maps in order to match the positions and pixel size of both datasets. To compute the BS maps, we apply the same methodology as defined in Sect. 3.1, using the reprojected H1 flux distributions, the 2-pixels Gaussian smoothed Hα velocity field and the single 2D Gaussian PSF to describe the new H1 angular resolution.

<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Pixel-by-pixel comparison of the BS contribution using Hα SINGS and H1 THINGS data. The y-axis shows the ratio of the BS amplitudes $Hα/H1$ and in the x-axis is the amplitude of the H1 model.</figcaption></figure>

In Fig. 12, we show the resulting $Hα \sigma_{bs}$ map (left panel), compared to the H1 $\sigma_{bs}$ map (middle panel), as well as the ratio of the two maps, for an example galaxy (NGC 3521). The BS pattern is similar in both cases, with the typical cross-shape in the centre. While the BS correction inferred from H1 data is less peaked in the centre than for $Hα$, the latter also decreases more slowly with radius. This explains the large values in the normalised ratio found far from the centre. Beam-smearing amplitudes are higher when modelling with the ionised gas. As the amplitude may be an important concern in our modelling, we study galaxy by galaxy the amplitude ratio for the two different cases.

In Fig. 13, a comparison of the BS contributions is made for the THINGS sample. For each galaxy, a 2D histogram with the amplitude of the H1 BS correction on the x-axis and the ratio of the BS corrections in $Hα$ to H1 on the y-axis, is made. Then, the histograms were normalised by the number of pixels in the maps, to avoid being biased by objects with a larger number of pixels. Figure 13 is then obtained by summing all the 2D histograms. It shows that the ratio is >1, which is not surprising given that the Hα BS benefits from a higher angular resolution. On average, the H1 BS is 2.3 times smaller than the Hα BS. We also observe a trend at low H1 BS with high ratio induced either by noise or by discontinuities between regions in the $Hα$ data. Indeed, a continuously increasing ratio with decreasing H1 BS can be modelled by assuming that the Hα BS is equal to the quadratic sum of
the H I BS and a constant term of ~3 to 10 km s⁻¹. This constant term can be related to pixel-to-pixel variations either induced by noise in the Hα data, or by true local velocity variations or discontinuities unseen in H I. Amplitudes vary within the sample, due to various intrinsic velocity gradients in each galaxy, as well as its distance and geometry.

Despite these findings, we decided to continue working with the H I σ<sub>bs</sub> for several reasons. First, the spatial extent of H I gas is larger, making it possible to evaluate the local BS contribution of spiral arms or other features that are not located in the inner disc regions. A second reason is that the Hα gas does not trace the same interstellar component, therefore the geometry of the discs might differ, and the flux distribution are noticeably different. Furthermore, it has been shown that the imprint of the systematic phase along the minor axis of the galaxy M 33 did not depend on the angular resolution (Chemin et al. 2020), while a range of 7 in the size of the synthesised beam was probed. Using H I maps to derive σ<sub>bs</sub> should thus not be an issue in the analysis (see also Sect. 5.2.4).

To assess the impact of BS under more conservative condi-
tions, we then subtracted respectively twice and thrice the H I BS
from the THINGS velocity dispersions. The result of this process
is that a large number of rings are removed in this case. Indeed,
in order to avoid sampling issues in the FFT, we removed rings
having more than 5% of pixels with undefined value after the BS
correction, due to this quadratic correction being larger than the
actual measured dispersion. The median fraction of such pixels
per ring is 0.7%, 4.0% and 16.6% when subtracting once, twice
and thrice the BS contribution, respectively. Only 621 and 258
tilted rings remain in the two latter cases, respectively, with less
than 5% of pixels with an undefined value, decreasing drastically
the statistics. With these particularly conservative BS subtrac-
tions, some galaxies totally disappear from the remaining rings,
which indicates that the velocity dispersion would be totally
dominated by BS, which should not be the case owing to the
(high) H I resolution.

In Fig. 14, we show the phase angles for *k* = 2 and *k* = 4 resulting from the subtraction of once or twice the contribution σ<sub>bs</sub> (red and blue histograms, respectively). We note that the amplitude distributions are not strongly affected. As mentioned before, more than 200 rings were removed, which explains why the average number is lower for the red histograms. For *k* = 2, the peak seen in the blue distribution is no more present in the red one. This behaviour is also observed for *k* = 1 (not shown here). For *k* = 4, in both cases we observe the strong peak at *φ*<sub>4</sub>/T<sub>4</sub> ~ 0.5. The interpretation is not straightforward, as the observed phase angle can be explained either by a residual BS signature or by the correlation of the asymmetry with the BS. Nevertheless, the fact that a strong peak remains in *φ*<sub>4</sub> tends to indicate that BS alone cannot explain the observed trends.

We also smoothed the Hα data at the average resolution of
the H I, 11" for all the sample with the goal to verify that the
same level of smoothing of the gradients implies a similar BS
effect. The results are as expected, being that the normalised
residuals for each galaxy show a distribution peaked on 0, indi-
cating that both modelling are consistent with each other.

Similarly, to verify the consistency of the findings described
above, we also performed the Fourier analysis of an example
galaxy (NGC 2841) using the 6" angular resolution THINGS
velocity and density maps, as obtained from robust weightings.
The amplitude of the BS is larger than for natural weightings,
which is due to the combined effects from the higher angular
resolution with the noisier velocity field for robust weightings
than for natural weightings. This finding matches perfectly the

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Comparison of the normalised phase angles histograms of orders *k* = 2 and *k* = 4 (left and right panels, respectively) in the velocity dispersion by subtracting once and twice the contribution from the BS effect (blue and red colours, respectively).</figcaption></figure>

result found with the Hα kinematics. Nevertheless, no signifi-
cant difference in the strength and phase angles of the Fourier
modes is observed with respect to the natural weighting for this
galaxy.

5.2.4. Amplitude of BS residuals

The correlation between BS and observed signatures found in
Sect. 5.2.1 may indicate that BS residuals remain after BS cor-
rection, due to the use of observed flux maps and velocity
fields rather than high resolution data to infer the correction (see
Sect. 3.1). To quantify residuals, we use the toy models described
in Sect. 3.3 and Appendix B. We focus on the isotropic case with
no perturbation and a steep inner gradient of the rotation curve,
for which the impact of BS is the largest, and added realistic
noise to dispersion maps (cf. Appendix A). Correcting for BS
as described in Sect. 3.1 reduces the strengths of orders k = 2
and k = 4 by a factor between 2 and 4 and is more efficient for
less inclined galaxies and for order k = 2. We also show that the
maximum strength of orders k = 2 and k = 4 of BS residuals
is reached at a lower radius than without BS correction, and that
it decreases for decreasing inclinations. For inclinations below
60°, residual BS cannot account for the k = 2 signature at radii
greater than a few beam FWHM and at radii >10 beam FWHM
for k = 4. At 75°; however, the BS signature remains well above
the uncertainties because (i) the projected velocity gets higher
than at lower inclination, and (ii) the spatial resolution along the
minor axis gets poor, inducing more hidden velocity gradients
within the beam and pixel size. We also performed an FFT anal-
ysis separately for the 50% largest and smallest radii and found,
as expected, that the amplitude of orders k = 2 and k = 4 is
reduced (by a factor larger than two) and that the significance of
the peak in phase histograms of the fourth order is much lower
(by a factor around 4) at large radii. For the model with a weaker
velocity gradient, when noise is added, all the systematic phase
angles disappear and the distributions become flat. The impact of
BS induced by large-scale rotation depends on the physical reso-
lution, the rotation curve shape, the projected maximum velocity,
galaxy inclination and high resolution flux distribution. A com-
plete study of BS on velocity dispersion maps depending on all
these parameters is beyond the scope of the present study.

About half galaxies in our sample have projected maximum
rotation velocities above 150 km s⁻¹, as inferred from the profile
width at 20% of the peak intensity given in Walter et al. (2008),
and about 60% have an inclination larger than 60°, which may
therefore present BS residuals signatures in the fourth order.
We first split the sample between inner and outer rings. The
the H I BS and a constant term of ~3 to 10 km s⁻¹. This constant term can be related to pixel-to-pixel variations either induced by noise in the Hα data, or by true local velocity variations or discontinuities unseen in H I. Amplitudes vary within the sample, due to various intrinsic velocity gradients in each galaxy, as well as its distance and geometry.

Despite these findings, we decided to continue working with the H I σ<sub>bs</sub> for several reasons. First, the spatial extent of H I gas is larger, making it possible to evaluate the local BS contribution of spiral arms or other features that are not located in the inner disc regions. A second reason is that the Hα gas does not trace the same interstellar component, therefore the geometry of the discs might differ, and the flux distribution are noticeably different. Furthermore, it has been shown that the imprint of the systematic phase along the minor axis of the galaxy M 33 did not depend on the angular resolution (Chemin et al. 2020), while a range of 7 in the size of the synthesised beam was probed. Using H I maps to derive σ<sub>bs</sub> should thus not be an issue in the analysis (see also Sect. 5.2.4).

To assess the impact of BS under more conservative condi-
tions, we then subtracted respectively twice and thrice the H I BS
from the THINGS velocity dispersions. The result of this process
is that a large number of rings are removed in this case. Indeed,
in order to avoid sampling issues in the FFT, we removed rings
having more than 5% of pixels with undefined value after the BS
correction, due to this quadratic correction being larger than the
actual measured dispersion. The median fraction of such pixels
per ring is 0.7%, 4.0% and 16.6% when subtracting once, twice
and thrice the BS contribution, respectively. Only 621 and 258
tilted rings remain in the two latter cases, respectively, with less
than 5% of pixels with an undefined value, decreasing drastically
the statistics. With these particularly conservative BS subtrac-
tions, some galaxies totally disappear from the remaining rings,
which indicates that the velocity dispersion would be totally
dominated by BS, which should not be the case owing to the
(high) H I resolution.

In Fig. 14, we show the phase angles for *k* = 2 and *k* = 4 resulting from the subtraction of once or twice the contribution σ<sub>bs</sub> (red and blue histograms, respectively). We note that the amplitude distributions are not strongly affected. As mentioned before, more than 200 rings were removed, which explains why the average number is lower for the red histograms. For *k* = 2, the peak seen in the blue distribution is no more present in the red one. This behaviour is also observed for *k* = 1 (not shown here). For *k* = 4, in both cases we observe the strong peak at *φ*<sub>4</sub>/T<sub>4</sub> ~ 0.5. The interpretation is not straightforward, as the observed phase angle can be explained either by a residual BS signature or by the correlation of the asymmetry with the BS. Nevertheless, the fact that a strong peak remains in *φ*<sub>4</sub> tends to indicate that BS alone cannot explain the observed trends.

We also smoothed the Hα data at the average resolution of
the H I, 11" for all the sample with the goal to verify that the
same level of smoothing of the gradients implies a similar BS
effect. The results are as expected, being that the normalised
residuals for each galaxy show a distribution peaked on 0, indi-
cating that both modelling are consistent with each other.

Similarly, to verify the consistency of the findings described
above, we also performed the Fourier analysis of an example
galaxy (NGC 2841) using the 6" angular resolution THINGS
velocity and density maps, as obtained from robust weightings.
The amplitude of the BS is larger than for natural weightings,
which is due to the combined effects from the higher angular
resolution with the noisier velocity field for robust weightings
than for natural weightings. This finding matches perfectly the

<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Comparison of the normalised phase angles histograms of orders *k* = 2 and *k* = 4 (left and right panels, respectively) in the velocity dispersion by subtracting once and twice the contribution from the BS effect (blue and red colours, respectively).</figcaption></figure>

result found with the Hα kinematics. Nevertheless, no signifi-
cant difference in the strength and phase angles of the Fourier
modes is observed with respect to the natural weighting for this
galaxy.

5.2.4. Amplitude of BS residuals

The correlation between BS and observed signatures found in
Sect. 5.2.1 may indicate that BS residuals remain after BS cor-
rection, due to the use of observed flux maps and velocity
fields rather than high resolution data to infer the correction (see
Sect. 3.1). To quantify residuals, we use the toy models described
in Sect. 3.3 and Appendix B. We focus on the isotropic case with
no perturbation and a steep inner gradient of the rotation curve,
for which the impact of BS is the largest, and added realistic
noise to dispersion maps (cf. Appendix A). Correcting for BS
as described in Sect. 3.1 reduces the strengths of orders k = 2
and k = 4 by a factor between 2 and 4 and is more efficient for
less inclined galaxies and for order k = 2. We also show that the
maximum strength of orders k = 2 and k = 4 of BS residuals
is reached at a lower radius than without BS correction, and that
it decreases for decreasing inclinations. For inclinations below
60°, residual BS cannot account for the k = 2 signature at radii
greater than a few beam FWHM and at radii >10 beam FWHM
for k = 4. At 75°; however, the BS signature remains well above
the uncertainties because (i) the projected velocity gets higher
than at lower inclination, and (ii) the spatial resolution along the
minor axis gets poor, inducing more hidden velocity gradients
within the beam and pixel size. We also performed an FFT anal-
ysis separately for the 50% largest and smallest radii and found,
as expected, that the amplitude of orders k = 2 and k = 4 is
reduced (by a factor larger than two) and that the significance of
the peak in phase histograms of the fourth order is much lower
(by a factor around 4) at large radii. For the model with a weaker
velocity gradient, when noise is added, all the systematic phase
angles disappear and the distributions become flat. The impact of
BS induced by large-scale rotation depends on the physical reso-
lution, the rotation curve shape, the projected maximum velocity,
galaxy inclination and high resolution flux distribution. A com-
plete study of BS on velocity dispersion maps depending on all
these parameters is beyond the scope of the present study.

About half galaxies in our sample have projected maximum
rotation velocities above 150 km s⁻¹, as inferred from the profile
width at 20% of the peak intensity given in Walter et al. (2008),
and about 60% have an inclination larger than 60°, which may
therefore present BS residuals signatures in the fourth order.
We first split the sample between inner and outer rings. The
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Histograms of normalised phase angles resulting from the FFT analysis for the THINGS sample splitting the rings as a function of radius (inner part in blue, outer part in red).</figcaption></figure>

FFT amplitudes of the innermost rings are 27% to 43% larger than for the outermost rings, depending on the order. The phase angle distributions are shown in Fig. 15, with the innermost and outermost rings shown in blue and red, respectively. These distributions are different. For *k* = 2, the distribution for the outermost rings is similar to that of the global distribution whereas the distribution does not peak on 0.5 φ₂/T₂ for the inner rings, as would be expected if the dispersion maps were dominated by BS residuals. The inner and outer regions both show a prominent peak *k* = 4 close to the half period, the peak being slightly offset towards lower phase angles in the outer parts. We also split the sample in low versus high projected maximum velocities, leading to the same sub-samples as for the absolute magnitude (see Sect. 4.4 and Fig. 8). As already discussed, the peak observed in the phase angle histograms of the fourth order is much more pronounced when projected velocities are high. This may be a proof that BS residuals are affecting our results at least for the fourth order. Nevertheless, given the large values of order four strengths σ₄ compared to what is expected from residuals in toy models, it may also be that such a signature is intrinsic to massive galaxies with high gas density contrasts. On the other hand, the peak in the second order around half period of galaxies with low projected velocities cannot be attributed to BS.

### 5.2.5. A possible link between velocity gradients and elevated velocity dispersion

We have shown that BS may explain the signatures we observe if not properly taken into account, as there are strong correlations between the modelled BS pattern and the observed velocity dispersion. Nevertheless, we have also shown that the amplitude of BS induced by large-scale motions is expected to be too low to have a significant impact on the data. The fact that BS seems to be coupled to observed signatures in the phase angles of asymmetries may indicate that the regions with the strongest velocity gradients or discontinuities, especially in the outer parts of galaxies, are also regions with an intrinsic large velocity dispersion at a higher level than what would be induced by BS from the observed gradients. These regions with large velocity gradients are often related to perturbations like bars, spiral arms, inter-arms and warps, especially at large radius, which means that there might be locally induced turbulence responsible for perturbations in both the velocity field and the velocity dispersion. This may also be that BS occurs on much smaller scale than that reached within our observations and that what we observe is not related to large-scale motions. A complementary way to investigate the effects of unresolved gradients of asymmetric velocities on the velocity dispersion is addressed in Sect. 5.3.1.

### 5.3. Projection effects in random motions

Finding systematic phase angles of asymmetries in velocity dispersion maps is a surprising result. The projection of asymmetric patterns along the line-of-sight should be randomly distributed. Flat distributions of phase angles were then expected in Fig. 5 from a theoretical perspective (see discussion in Appendix C). Privileged phases thus imply a dependency of the patterns on the orientation of the disc principal axes with respect to the observer, because the major axis of discs has been chosen at the origin of the azimuthal angles (θ = 0 along the major axis). In other words, this suggests the presence of projection effects in the velocity dispersion maps. In this section, we thus want to study the possible origin(s) of projection effects in the velocity dispersion. This is achieved by investigating the effect of asymmetries in the ordered motions on the random motions (Sect. 5.3.1). We also address the effect of correlated velocity components (Sect. 5.3.2).

#### 5.3.1. Effects from the asymmetric ordered motions

It is usual to assume that the velocity ellipsoid of gas is isotropic in galactic discs. Under this assumption, no deprojection of data is required, so that the observed velocity dispersion is a direct proxy of the gaseous random motions. However, this assumption is very idealised, and can only apply to axisymmetric kinematics, which rarely occurs in reality. Indeed, the velocity fields (Vᵢ₀s) always exhibit disturbances, as caused by the projection of asymmetric ordered motions (Vᵣ, Vθ, Vz) due to, for example, bars, rings, spiral arms, or warps. The effect of any of such disturbances in Vᵢ₀s should also propagate to σᵢ₀s because of the implied gradients in one or several directions in the plane. We test this hypothesis by studying the toy models with asymmetric kinematics perturbations described in Sect. 3.3 and Appendix B in both isotropic and anisotropic cases. We mainly focus hereafter on cases with the weakest velocity gradient (shown in Fig. B.2) to reduce possible residual BS signatures related to unresolved rotation (see Sect. 5.2.4) and emphasise on local perturbations.

In the axisymmetric and isotropic model (used as reference), the strength of the Fourier coefficients of the various orders are comparable (~0.5 km s⁻¹ on average). In all other models, σ₁ and σ₃ barely vary from this value, while the *k* = 2 and *k* = 4 components are significantly larger, being mostly >2 km s⁻¹, except in the axisymmetric case with uniform local anisotropy, where σ₄ is very comparable to the odd amplitudes. The distributions of the phase angles for the negligible modes *k* = 1 and *k* = 3 do not show systematic peaks. The even orders thus always dominate the dispersion asymmetries, and the bisymmetry is the strongest perturbation, on average.

Any model with a velocity ellipsoid showing a radial bias σθ < σR creates a bisymmetry theoretically aligned with the minor axis. The strength of *k* = 2 asymmetries in σᵢ₀s within mock data obtained from models combining effects from the uniform radial bias and Vᵣ and/or Vθ streamings is between 20 and 40% larger than those without the radial bias, on average, while the amplitude of the mode *k* = 4 is not affected. The distribution of φ₂ is sharply peaked at π/2 for the axisymmetric case of anisotropy with the radial bias, and the peak is slightly enlarged when there are asymmetries in Vᵣ and/or Vθ, regardless of their strength. Moreover, in this latter case with the velocity perturbation exclusively set along the azimuthal direction (ΔVᵣ = 0, ΔVθ ≠ 0), the incidence of φ₂ = 0 or 1 remains extremely low. These findings can be explained with the help of Eq. (3)
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Histograms of normalised phase angles resulting from the FFT analysis for the THINGS sample splitting the rings as a function of radius (inner part in blue, outer part in red).</figcaption></figure>

FFT amplitudes of the innermost rings are 27% to 43% larger than for the outermost rings, depending on the order. The phase angle distributions are shown in Fig. 15, with the innermost and outermost rings shown in blue and red, respectively. These distributions are different. For *k* = 2, the distribution for the outermost rings is similar to that of the global distribution whereas the distribution does not peak on 0.5 φ₂/T₂ for the inner rings, as would be expected if the dispersion maps were dominated by BS residuals. The inner and outer regions both show a prominent peak *k* = 4 close to the half period, the peak being slightly offset towards lower phase angles in the outer parts. We also split the sample in low versus high projected maximum velocities, leading to the same sub-samples as for the absolute magnitude (see Sect. 4.4 and Fig. 8). As already discussed, the peak observed in the phase angle histograms of the fourth order is much more pronounced when projected velocities are high. This may be a proof that BS residuals are affecting our results at least for the fourth order. Nevertheless, given the large values of order four strengths σ₄ compared to what is expected from residuals in toy models, it may also be that such a signature is intrinsic to massive galaxies with high gas density contrasts. On the other hand, the peak in the second order around half period of galaxies with low projected velocities cannot be attributed to BS.

### 5.2.5. A possible link between velocity gradients and elevated velocity dispersion

We have shown that BS may explain the signatures we observe if not properly taken into account, as there are strong correlations between the modelled BS pattern and the observed velocity dispersion. Nevertheless, we have also shown that the amplitude of BS induced by large-scale motions is expected to be too low to have a significant impact on the data. The fact that BS seems to be coupled to observed signatures in the phase angles of asymmetries may indicate that the regions with the strongest velocity gradients or discontinuities, especially in the outer parts of galaxies, are also regions with an intrinsic large velocity dispersion at a higher level than what would be induced by BS from the observed gradients. These regions with large velocity gradients are often related to perturbations like bars, spiral arms, inter-arms and warps, especially at large radius, which means that there might be locally induced turbulence responsible for perturbations in both the velocity field and the velocity dispersion. This may also be that BS occurs on much smaller scale than that reached within our observations and that what we observe is not related to large-scale motions. A complementary way to investigate the effects of unresolved gradients of asymmetric velocities on the velocity dispersion is addressed in Sect. 5.3.1.

### 5.3. Projection effects in random motions

Finding systematic phase angles of asymmetries in velocity dispersion maps is a surprising result. The projection of asymmetric patterns along the line-of-sight should be randomly distributed. Flat distributions of phase angles were then expected in Fig. 5 from a theoretical perspective (see discussion in Appendix C). Privileged phases thus imply a dependency of the patterns on the orientation of the disc principal axes with respect to the observer, because the major axis of discs has been chosen at the origin of the azimuthal angles (θ = 0 along the major axis). In other words, this suggests the presence of projection effects in the velocity dispersion maps. In this section, we thus want to study the possible origin(s) of projection effects in the velocity dispersion. This is achieved by investigating the effect of asymmetries in the ordered motions on the random motions (Sect. 5.3.1). We also address the effect of correlated velocity components (Sect. 5.3.2).

#### 5.3.1. Effects from the asymmetric ordered motions

It is usual to assume that the velocity ellipsoid of gas is isotropic in galactic discs. Under this assumption, no deprojection of data is required, so that the observed velocity dispersion is a direct proxy of the gaseous random motions. However, this assumption is very idealised, and can only apply to axisymmetric kinematics, which rarely occurs in reality. Indeed, the velocity fields (Vᵢ₀s) always exhibit disturbances, as caused by the projection of asymmetric ordered motions (Vᵣ, Vθ, Vz) due to, for example, bars, rings, spiral arms, or warps. The effect of any of such disturbances in Vᵢ₀s should also propagate to σᵢ₀s because of the implied gradients in one or several directions in the plane. We test this hypothesis by studying the toy models with asymmetric kinematics perturbations described in Sect. 3.3 and Appendix B in both isotropic and anisotropic cases. We mainly focus hereafter on cases with the weakest velocity gradient (shown in Fig. B.2) to reduce possible residual BS signatures related to unresolved rotation (see Sect. 5.2.4) and emphasise on local perturbations.

In the axisymmetric and isotropic model (used as reference), the strength of the Fourier coefficients of the various orders are comparable (~0.5 km s⁻¹ on average). In all other models, σ₁ and σ₃ barely vary from this value, while the *k* = 2 and *k* = 4 components are significantly larger, being mostly >2 km s⁻¹, except in the axisymmetric case with uniform local anisotropy, where σ₄ is very comparable to the odd amplitudes. The distributions of the phase angles for the negligible modes *k* = 1 and *k* = 3 do not show systematic peaks. The even orders thus always dominate the dispersion asymmetries, and the bisymmetry is the strongest perturbation, on average.

Any model with a velocity ellipsoid showing a radial bias σθ < σR creates a bisymmetry theoretically aligned with the minor axis. The strength of *k* = 2 asymmetries in σᵢ₀s within mock data obtained from models combining effects from the uniform radial bias and Vᵣ and/or Vθ streamings is between 20 and 40% larger than those without the radial bias, on average, while the amplitude of the mode *k* = 4 is not affected. The distribution of φ₂ is sharply peaked at π/2 for the axisymmetric case of anisotropy with the radial bias, and the peak is slightly enlarged when there are asymmetries in Vᵣ and/or Vθ, regardless of their strength. Moreover, in this latter case with the velocity perturbation exclusively set along the azimuthal direction (ΔVᵣ = 0, ΔVθ ≠ 0), the incidence of φ₂ = 0 or 1 remains extremely low. These findings can be explained with the help of Eq. (3)
from Chemin et al. (2020), where having $k = 2$ asymmetries in $\sigma_{los}$ aligned with the minor axis is a genuine imprint of a radially biased velocity ellipsoid (see also Eq. (C.11), but with null cross terms). Thus, a local and uniform anisotropy of the velocity ellipsoid with a radial bias drives the phases of asymmetries in $\sigma_{los}$ more strongly than the effect of the velocity anisotropy induced at larger scale by the projection of $V_R$ and $V_\theta$ gradients. This is because the radial bias is applied to all velocities in the mock data, while the $V_R$ and $V_\theta$ streamings are only restricted to the locations of the perturbations. If a tangential bias $\sigma_\theta > \sigma_R$ had been assumed, the bisymmetry would have been exclusively aligned along the disc major axis, to the detriment of the minor axis. Because THINGS and WHISP galaxies do not show a $\phi_2$ distribution strongly concentrated around $\pi/2$, anisotropic gaseous velocity ellipsoids consistent with a radial bias uniformly distributed everywhere in discs, like the one presented here, seem very unlikely. A solution to avoid the $\phi_2/T_2$ concentration at 0.5 is to have either an axis ratio $\sigma_\theta/\sigma_R$ which varies with radius, or tilted velocity ellipsoids. Section 5.3.2 provides a hint of how this goal could be achieved.

The impact of velocity anisotropy arising from asymmetries in $V_R$ (or $V_\theta$, respectively) is to greatly decrease the likelihoods of finding $\phi_2$ close to zero ($\pi/2$), in the case the velocity perturbation is along the radial axis, that is $\Delta V_\theta = 0$ and $\Delta V_R \neq 0$ (azimuthal axis, $\Delta V_R = 0$, $\Delta V_\theta \neq 0$). This is because there is little projection of $V_\theta$ ($V_R$, respectively) onto the minor (major) axis. These results remain valid regardless of the uniform isotropy/radial bias assumptions. If there were an equal number of tilted rings dominated by radial and azimuthal velocity streamings, this would lead to distributions of $\phi_2$ with minima around 0 and $\pi/2$, and conversely with maxima $\phi_2 = 0.25\pi$ and $0.75\pi$, which values are not far from those observed in the sample within the quoted uncertainties.

Then, still for models with asymmetries in $V_R$ and/or $V_\theta$, the distributions of $\sigma_4$ and $\phi_4$ are very similar in both cases of uniform isotropy or radial bias. Thus, the uniform velocity anisotropy has negligible effect on the fourth order perturbation of $\sigma_{los}$, and this latter is an harmonics of the second order perturbation in our toy models. The information about the orientation of the bisymmetry (i.e. along the major or minor axis) is nevertheless lost modulo $\pi/2$ in the fourth order<fnref n="5" />. Moreover, it is obvious that asymmetries have destroyed the systematic peak at $\phi_4 = 0.5 T_4$ seen in the reference axisymmetric models. This is because this peak could be caused by small residuals from BS correction, as mentioned previously. Unless our velocity asymmetries in the toy models are significantly overestimated with respect to the reality, the occurrence of a systematic peak at $\phi_4 = 0.5 T_4$ in the observations indicates that the fourth order asymmetry is very likely an harmonic of a second order which would have preferential phase angles $\phi_2$ close to either $0.25\pi$ and $0.75\pi$, which values remains close to the observed $\phi_2 \sim 0.35\pi$ and $0.65\pi$ within the quoted uncertainty of $0.1 T_2$.

We note that combining a steep velocity gradient top the perturbations (see Fig. B.3) does not change qualitatively the previous results for the order $k = 2$ though the contrast between high and low probability is a bit reduced, especially in the isotropic case with the velocity perturbation along $V_R$, due to residual BS inducing a larger probability for null phase angles whatever the direction of the perturbation. Anisotropy always dominates order $k = 2$ when present. For order $k = 4$, a peak related to residual

BS clearly remains around $\phi_4/T_4 = 0.5$, which increases the trend already observed with weak rotation-induced BS.

Finally, although not shown here, we observe that when the BS correction is made with flux density and velocity maps at the highest angular resolution of the mock data (native resolution of ~50 pc, before spatial smoothing), no significant residuals due to BS remain in the distribution of phase angles. We note here that this value of 50 pc is only applicable to our toy models, and not to real observations, because we do not know precisely the angular scale of velocity streamings in galaxies. This indicates that the velocity anisotropy from the streamings can only affect the phase angles when the streamings are not properly resolved. Therefore, our method consisting in searching for systematic effects in the phase angles of $\sigma_{los}$ from observations of 'limited' resolutions like THINGS is a powerful tool to predict the directions in the plane towards which the perturbations and anisotropy are more significant. In particular, it implies that the lopsidedness perturbation of H I kinematics in THINGS is more likely to be stronger along the tangential dimension than along the radial direction, hence a preferentially large-scale tangential bias in the velocity anisotropy, in order to explain the higher incidence of $\phi_1$ along the major axis of the galaxies. Similarly, it is less probable to find a bisymmetry perturbation stronger in the tangential direction, to explain the lower incidence of $\phi_2$ towards the major axis of the galaxies. The observations of $\phi_2$ aligned near the disc minor axis in lower mass galaxies (Sect. 4.4), similarly to the outer part of M 33, would thus likely be explained by dominant radial streamings, hence a stronger large-scale radial bias in the velocity anisotropy (Chemin et al. 2020). Another consequence of the analysis is that velocity asymmetries in the first moment maps of galaxies must be observed on a scale of about the PSF width, though of reduced amplitude. We thus expect that systematic effects are present in the phase angles of $V_{los}$ asymmetries, by analogy with those evidenced in $\sigma_{los}$, tracing the direction of dominant bias of the velocity anisotropy from the perturbations. Comparing the properties of asymmetries of $V_{los}$ and $\sigma_{los}$ is beyond the scope of this article, however.

### 5.3.2. Effects from anisotropic velocity ellipsoids with tilted axes

Until now, the tests on the shape of the velocity ellipsoid has only considered components of the velocity tensor which are independent from each other. However, the asymmetries in $V_R$, $V_\theta$ and $V_z$ all result from the same perturbed gravitational potential. By consequence, the components must show a certain degree of correlation, at least at the vicinity of perturbations. Correlated velocities are evidenced in the kinematics of stellar populations of the Milky Way, through tilted axes of the stellar velocity ellipsoid. The observed tilt angle of the axes of the radial and vertical velocity distribution of Galactic halo stars can inform on the 3D shape of the mass distribution and the dark matter halo (Siebert et al. 2008; Smith et al. 2009; Wegg et al. 2019), while that of the radial and tangential velocity distribution (the vertex deviation) is caused by a bisymmetric perturbation in the disc, such as a bar, a spiral structure, a warp, or more simply an elliptical disc (Kuijken & Tremaine 1991; Debattista et al. 2019).

We can test the effect of correlations between H I velocities by means of the uniform velocity anisotropy presented above, by considering the possibility of tilted ellipsoids, as we do for stars. This has a consequence on the projection of the ellipsoid on the line-of-sight, and following the complete mathematical description from Appendix C, the squared velocity dispersion is

<footnote n="5">⁵ A 4th-order mode could be generated as an harmonic of a second order mode with a phase of either $\phi_2 = \phi_4$ or $\phi_2 = \phi_4 + \pi/2$.</footnote>
from Chemin et al. (2020), where having $k = 2$ asymmetries in $\sigma_{los}$ aligned with the minor axis is a genuine imprint of a radially biased velocity ellipsoid (see also Eq. (C.11), but with null cross terms). Thus, a local and uniform anisotropy of the velocity ellipsoid with a radial bias drives the phases of asymmetries in $\sigma_{los}$ more strongly than the effect of the velocity anisotropy induced at larger scale by the projection of $V_R$ and $V_\theta$ gradients. This is because the radial bias is applied to all velocities in the mock data, while the $V_R$ and $V_\theta$ streamings are only restricted to the locations of the perturbations. If a tangential bias $\sigma_\theta > \sigma_R$ had been assumed, the bisymmetry would have been exclusively aligned along the disc major axis, to the detriment of the minor axis. Because THINGS and WHISP galaxies do not show a $\phi_2$ distribution strongly concentrated around $\pi/2$, anisotropic gaseous velocity ellipsoids consistent with a radial bias uniformly distributed everywhere in discs, like the one presented here, seem very unlikely. A solution to avoid the $\phi_2/T_2$ concentration at 0.5 is to have either an axis ratio $\sigma_\theta/\sigma_R$ which varies with radius, or tilted velocity ellipsoids. Section 5.3.2 provides a hint of how this goal could be achieved.

The impact of velocity anisotropy arising from asymmetries in $V_R$ (or $V_\theta$, respectively) is to greatly decrease the likelihoods of finding $\phi_2$ close to zero ($\pi/2$), in the case the velocity perturbation is along the radial axis, that is $\Delta V_\theta = 0$ and $\Delta V_R \neq 0$ (azimuthal axis, $\Delta V_R = 0$, $\Delta V_\theta \neq 0$). This is because there is little projection of $V_\theta$ ($V_R$, respectively) onto the minor (major) axis. These results remain valid regardless of the uniform isotropy/radial bias assumptions. If there were an equal number of tilted rings dominated by radial and azimuthal velocity streamings, this would lead to distributions of $\phi_2$ with minima around 0 and $\pi/2$, and conversely with maxima $\phi_2 = 0.25\pi$ and $0.75\pi$, which values are not far from those observed in the sample within the quoted uncertainties.

Then, still for models with asymmetries in $V_R$ and/or $V_\theta$, the distributions of $\sigma_4$ and $\phi_4$ are very similar in both cases of uniform isotropy or radial bias. Thus, the uniform velocity anisotropy has negligible effect on the fourth order perturbation of $\sigma_{los}$, and this latter is an harmonics of the second order perturbation in our toy models. The information about the orientation of the bisymmetry (i.e. along the major or minor axis) is nevertheless lost modulo $\pi/2$ in the fourth order<fnref n="5" />. Moreover, it is obvious that asymmetries have destroyed the systematic peak at $\phi_4 = 0.5 T_4$ seen in the reference axisymmetric models. This is because this peak could be caused by small residuals from BS correction, as mentioned previously. Unless our velocity asymmetries in the toy models are significantly overestimated with respect to the reality, the occurrence of a systematic peak at $\phi_4 = 0.5 T_4$ in the observations indicates that the fourth order asymmetry is very likely an harmonic of a second order which would have preferential phase angles $\phi_2$ close to either $0.25\pi$ and $0.75\pi$, which values remains close to the observed $\phi_2 \sim 0.35\pi$ and $0.65\pi$ within the quoted uncertainty of $0.1 T_2$.

We note that combining a steep velocity gradient top the perturbations (see Fig. B.3) does not change qualitatively the previous results for the order $k = 2$ though the contrast between high and low probability is a bit reduced, especially in the isotropic case with the velocity perturbation along $V_R$, due to residual BS inducing a larger probability for null phase angles whatever the direction of the perturbation. Anisotropy always dominates order $k = 2$ when present. For order $k = 4$, a peak related to residual

BS clearly remains around $\phi_4/T_4 = 0.5$, which increases the trend already observed with weak rotation-induced BS.

Finally, although not shown here, we observe that when the BS correction is made with flux density and velocity maps at the highest angular resolution of the mock data (native resolution of ~50 pc, before spatial smoothing), no significant residuals due to BS remain in the distribution of phase angles. We note here that this value of 50 pc is only applicable to our toy models, and not to real observations, because we do not know precisely the angular scale of velocity streamings in galaxies. This indicates that the velocity anisotropy from the streamings can only affect the phase angles when the streamings are not properly resolved. Therefore, our method consisting in searching for systematic effects in the phase angles of $\sigma_{los}$ from observations of 'limited' resolutions like THINGS is a powerful tool to predict the directions in the plane towards which the perturbations and anisotropy are more significant. In particular, it implies that the lopsidedness perturbation of H I kinematics in THINGS is more likely to be stronger along the tangential dimension than along the radial direction, hence a preferentially large-scale tangential bias in the velocity anisotropy, in order to explain the higher incidence of $\phi_1$ along the major axis of the galaxies. Similarly, it is less probable to find a bisymmetry perturbation stronger in the tangential direction, to explain the lower incidence of $\phi_2$ towards the major axis of the galaxies. The observations of $\phi_2$ aligned near the disc minor axis in lower mass galaxies (Sect. 4.4), similarly to the outer part of M 33, would thus likely be explained by dominant radial streamings, hence a stronger large-scale radial bias in the velocity anisotropy (Chemin et al. 2020). Another consequence of the analysis is that velocity asymmetries in the first moment maps of galaxies must be observed on a scale of about the PSF width, though of reduced amplitude. We thus expect that systematic effects are present in the phase angles of $V_{los}$ asymmetries, by analogy with those evidenced in $\sigma_{los}$, tracing the direction of dominant bias of the velocity anisotropy from the perturbations. Comparing the properties of asymmetries of $V_{los}$ and $\sigma_{los}$ is beyond the scope of this article, however.

### 5.3.2. Effects from anisotropic velocity ellipsoids with tilted axes

Until now, the tests on the shape of the velocity ellipsoid has only considered components of the velocity tensor which are independent from each other. However, the asymmetries in $V_R$, $V_\theta$ and $V_z$ all result from the same perturbed gravitational potential. By consequence, the components must show a certain degree of correlation, at least at the vicinity of perturbations. Correlated velocities are evidenced in the kinematics of stellar populations of the Milky Way, through tilted axes of the stellar velocity ellipsoid. The observed tilt angle of the axes of the radial and vertical velocity distribution of Galactic halo stars can inform on the 3D shape of the mass distribution and the dark matter halo (Siebert et al. 2008; Smith et al. 2009; Wegg et al. 2019), while that of the radial and tangential velocity distribution (the vertex deviation) is caused by a bisymmetric perturbation in the disc, such as a bar, a spiral structure, a warp, or more simply an elliptical disc (Kuijken & Tremaine 1991; Debattista et al. 2019).

We can test the effect of correlations between H I velocities by means of the uniform velocity anisotropy presented above, by considering the possibility of tilted ellipsoids, as we do for stars. This has a consequence on the projection of the ellipsoid on the line-of-sight, and following the complete mathematical description from Appendix C, the squared velocity dispersion is

<footnote n="5">⁵ A 4th-order mode could be generated as an harmonic of a second order mode with a phase of either $\phi_2 = \phi_4$ or $\phi_2 = \phi_4 + \pi/2$.</footnote>
given by:

$$ \sigma_{\text{los}}^2 = \left( \frac{\sigma_{\theta}^2 - \sigma_{\text{R}}^2}{2} \cos 2\theta + \sigma_{\text{R}\theta}^2 \sin 2\theta \right) \sin^2 i \\ + \left( \sigma_{\text{R}z}^2 \cos \theta + \sigma_{\theta z}^2 \sin \theta \right) \sin 2i \\ + \sigma_z^2 \cos^2 i + \frac{\sigma_{\theta}^2 + \sigma_{\text{R}}^2}{2} \sin^2 i, \qquad (6) $$

where $\sigma_{\text{R}\theta}^2$, $\sigma_{\theta z}^2$ and $\sigma_{\text{R}z}^2$ are the covariance terms marking the correlation between the components of the velocity tensor in the cylindrical frame of a disc. This equation can then be expressed as a sum of three terms of orders 0, 1 and 2, with an analogue formalism as the one we used to study asymmetries:

$$ \sigma_{\text{los}}^2 = k_0 + k_1 \cos(\theta - \alpha_1) + k_2 \cos(2(\theta - \alpha_2)), \qquad (7) $$

where $k_i$ are the amplitudes of these terms and $\alpha_i$ are their phases, that all depend on the various velocity dispersion terms. More specifically, we have

$$ \tan 2\alpha_2 = 2\sigma_{\text{R}\theta}^2/(\sigma_{\theta}^2 - \sigma_{\text{R}}^2). \qquad (8) $$

Assuming no tilted axes for the uniform velocity anisotropy in the previous section left no choice than having a dispersion bisymmetry systematically aligned with the minor or major axis ($\alpha_2 = \pi/2$ or 0, respectively), depending on the sign of $\sigma_{\theta}^2 - \sigma_{\text{R}}^2$. This is at odds with the observation that the second order phase angles $\phi_2 = 0.35T_2$ and $\phi_2 = 0.65T_2$ are particularly more likely. Equations (7) and (8) thus imply that it is possible to get various $\phi_2$ orientations if the velocity anisotropy is not as uniform through the disc, and more importantly if the covariance varies. In the Milky Way disc, the measurement of the axis ratio of the tangential-to-radial velocity dispersions of young stars as measured with Gaia spectro-astrometric kinematics, is $\sigma_{\theta}/\sigma_{\text{R}} = 0.66$ (Gaia Collaboration 2023). Since young stars kinematics are expected to be comparable to that of the gas, we further assume for this exercise that the velocity anisotropy of the planar components of the interstellar gas for any galaxies in our sample is comparable to the value found for young stars in the MW disc. Reminding now that $\sigma_{\text{R}\theta}^2 = \rho_{\text{R}\theta}\sigma_{\text{R}}\sigma_{\theta}$, with $\rho_{\text{R}\theta}$ the correlation coefficient, an appropriate choice of covariance and correlation coefficient which explain the higher incidence of $\phi_2 \equiv \alpha_2 = 0.35\pi$ and $0.65\pi$ for gas is $\rho_{\text{R}\theta} = \pm 0.6$. Furthermore, Eq. (7) implies that correlations between the vertical and the two planar components should also be considered to get systematic phases of first order phases. In order to explain those seen in Fig. 5, we would thus need to have $\alpha_1 \sim 0$ or $\pi$ to match the observations, that is, $\tan \alpha_1 = \sigma_{\theta z}^2/\sigma_{\text{R}z}^2 \sim 0$, following Eq. (C.13). This is possible if $\rho_{\theta z}\sigma_{\theta} \ll \rho_{\text{R}z}\sigma_{\text{R}}$. Here again, studies of the kinematics of young stars in the Galaxy should help to test these predictions, which is relevant only if the shape of gaseous velocity ellipsoids are comparable to those of young stars in nearby discs. It would then be interesting to investigate how these predictions on $\rho_{\text{R}\theta}$ compare with the observed planar correlations of young stars of the Galactic disc, and how $\rho_{\theta z}\sigma_{\theta}$ compare with $\rho_{\text{R}z}\sigma_{\text{R}}$ for the same stellar populations. To our knowledge, none of such distributions have ever been measured yet for the Milky Way.

# 6. Summary and concluding remarks

For this work, we have performed the first systematic search and characterisation of asymmetries in velocity dispersion maps of H I gas in nearby galaxies. We used the best H I data available

in the archives, namely 32 galaxies from THINGS and WHISP samples with a spatial resolution better than 1.5 kpc and a neutral gas distribution extended by at least seven times the resolution, to carry out this pioneering work. We performed a Fourier analysis of the observed second moment maps previously corrected from both the line spread function and the two-dimensional BS contribution tracing the velocity gradients occurring at a large angular scale. This allowed us to measure the strength and phase angle for each of the first four Fourier harmonics, which are the most important perturbations in the random motions of H I gas.

We find a wide range of strength of asymmetries, up to ~35 km s⁻¹. The shape of the likelihood distributions of the amplitude of the Fourier modes are similar, with a peak around 5–10 km s⁻¹, and an extended tail towards larger values. Overall, the strongest Fourier mode is the bisymmetry ($k = 2$), with a mean amplitude of ~11 km s⁻¹, followed by the lopsidedness $k = 1$ asymmetry. On the other hand, the shape of the likelihood distributions of phase angles strongly depends on the Fourier mode, and the distributions present significant variations. For the $k = 1$ asymmetry, the probability to find a phase angle near the major axis of the discs is larger, while finding an asymmetry oriented near 0.15 and 0.75 × 2π is less probable. For $k = 2$, the likelihood is larger at phase angles of ~0.35 and 0.65 × π, and lower near the disc major axis, whereas that of the $k = 4$ mode is concentrated at a value of 0.5 × π/2, implying an orientation mainly lying at an equivalent separation from the disc major and minor axes. These systematic phase angles are robust against the number of titled rings and against the significance of the detected signal. Using a uniform number of rings for all galaxies lowers the significance of the result, but the trends for the shape of the probability distributions are preserved. More uniform likelihood distributions were expected for randomly distributed perturbations in the galaxies and with isotropic velocity ellipsoids for the H I gas. This is evidence that strong projection effects imprint on the velocity dispersion maps of H I gas in nearby galaxies.

Finding the origin(s) of the asymmetries and its systematics is a difficult task because of the nature of the velocity dispersion itself, which stems from instrumental effects and various processes inherent to the galactic random motions (dynamical, hydrodynamical, and thermal). In that aspect, we have shown that large-scale-induced BS might produce non-negligible residuals, especially in the fourth order for highly inclined galaxies with high projected rotation velocities. Typical signatures of BS on phase angles are orientated near the major axis ($k = 1$), the minor axis ($k = 2$), and at $0.5 \times \pi/2$ ($k = 4$). Nevertheless, these orientations are not fully consistent with the observations, especially for the second order, and the strength of Fourier modes in the BS contribution is much smaller than the observed asymmetries. This suggests that unresolved velocity streamings could be responsible for some variations in the likelihood distributions of the phase angles. The correlations found between Fourier decomposition of H I gas density and velocity dispersion maps suggest that the observed bisymmetry, and thus streaming motions, are induced by spiral arms, bars, or warps.

We thus further addressed the impact on the random motions from such velocity streamings present in the ordered motions $V_R$ and $V_{\theta}$. To do so, toy models mocking discs with a barely resolved spiral perturbation in both the density and kinematics were made, from which mock moment maps of discs seen under various projection angles, strengths, and orientations, of the perturbations were derived, and Fourier analyses were performed, by analogy to the observations. When streamings in the ordered motions are stronger in the radial dimension, hence in the presence of a radially biased velocity anisotropy on large scales, they
given by:

$$ \sigma_{\text{los}}^2 = \left( \frac{\sigma_{\theta}^2 - \sigma_R^2}{2} \cos 2\theta + \sigma_{R\theta}^2 \sin 2\theta \right) \sin^2 i \\ + \left( \sigma_{Rz}^2 \cos \theta + \sigma_{\theta z}^2 \sin \theta \right) \sin 2i \\ + \sigma_z^2 \cos^2 i + \frac{\sigma_{\theta}^2 + \sigma_R^2}{2} \sin^2 i, \qquad (6) $$

where $\sigma_{R\theta}^2$, $\sigma_{\theta z}^2$ and $\sigma_{Rz}^2$ are the covariance terms marking the correlation between the components of the velocity tensor in the cylindrical frame of a disc. This equation can then be expressed as a sum of three terms of orders 0, 1 and 2, with an analogue formalism as the one we used to study asymmetries:

$$ \sigma_{\text{los}}^2 = k_0 + k_1 \cos(\theta - \alpha_1) + k_2 \cos(2(\theta - \alpha_2)), \qquad (7) $$

where $k_i$ are the amplitudes of these terms and $\alpha_i$ are their phases, that all depend on the various velocity dispersion terms. More specifically, we have

$$ \tan 2\alpha_2 = 2\sigma_{R\theta}^2/(\sigma_{\theta}^2 - \sigma_R^2). \qquad (8) $$

Assuming no tilted axes for the uniform velocity anisotropy in the previous section left no choice than having a dispersion bisymmetry systematically aligned with the minor or major axis ($\alpha_2 = \pi/2$ or 0, respectively), depending on the sign of $\sigma_{\theta}^2 - \sigma_R^2$. This is at odds with the observation that the second order phase angles $\phi_2 = 0.35T_2$ and $\phi_2 = 0.65T_2$ are particularly more likely. Equations (7) and (8) thus imply that it is possible to get various $\phi_2$ orientations if the velocity anisotropy is not as uniform through the disc, and more importantly if the covariance varies. In the Milky Way disc, the measurement of the axis ratio of the tangential-to-radial velocity dispersions of young stars as measured with Gaia spectro-astrometric kinematics, is $\sigma_{\theta}/\sigma_R = 0.66$ (Gaia Collaboration 2023). Since young stars kinematics are expected to be comparable to that of the gas, we further assume for this exercise that the velocity anisotropy of the planar components of the interstellar gas for any galaxies in our sample is comparable to the value found for young stars in the MW disc. Reminding now that $\sigma_{R\theta}^2 = \rho_{R\theta}\sigma_R\sigma_{\theta}$, with $\rho_{R\theta}$ the correlation coefficient, an appropriate choice of covariance and correlation coefficient which explain the higher incidence of $\phi_2 \equiv \alpha_2 = 0.35\pi$ and $0.65\pi$ for gas is $\rho_{R\theta} = \pm 0.6$. Furthermore, Eq. (7) implies that correlations between the vertical and the two planar components should also be considered to get systematic phases of first order phases. In order to explain those seen in Fig. 5, we would thus need to have $\alpha_1 \sim 0$ or $\pi$ to match the observations, that is, $\tan \alpha_1 = \sigma_{\theta z}^2/\sigma_{Rz}^2 \sim 0$, following Eq. (C.13). This is possible if $\rho_{\theta z}\sigma_{\theta} \ll \rho_{Rz}\sigma_R$. Here again, studies of the kinematics of young stars in the Galaxy should help to test these predictions, which is relevant only if the shape of gaseous velocity ellipsoids are comparable to those of young stars in nearby discs. It would then be interesting to investigate how these predictions on $\rho_{R\theta}$ compare with the observed planar correlations of young stars of the Galactic disc, and how $\rho_{\theta z}\sigma_{\theta}$ compare with $\rho_{Rz}\sigma_R$ for the same stellar populations. To our knowledge, none of such distributions have ever been measured yet for the Milky Way.

# 6. Summary and concluding remarks

For this work, we have performed the first systematic search and characterisation of asymmetries in velocity dispersion maps of H I gas in nearby galaxies. We used the best H I data available

in the archives, namely 32 galaxies from THINGS and WHISP samples with a spatial resolution better than 1.5 kpc and a neutral gas distribution extended by at least seven times the resolution, to carry out this pioneering work. We performed a Fourier analysis of the observed second moment maps previously corrected from both the line spread function and the two-dimensional BS contribution tracing the velocity gradients occurring at a large angular scale. This allowed us to measure the strength and phase angle for each of the first four Fourier harmonics, which are the most important perturbations in the random motions of H I gas.

We find a wide range of strength of asymmetries, up to ~35 km s⁻¹. The shape of the likelihood distributions of the amplitude of the Fourier modes are similar, with a peak around 5–10 km s⁻¹, and an extended tail towards larger values. Overall, the strongest Fourier mode is the bisymmetry ($k = 2$), with a mean amplitude of ~11 km s⁻¹, followed by the lopsidedness $k = 1$ asymmetry. On the other hand, the shape of the likelihood distributions of phase angles strongly depends on the Fourier mode, and the distributions present significant variations. For the $k = 1$ asymmetry, the probability to find a phase angle near the major axis of the discs is larger, while finding an asymmetry oriented near 0.15 and 0.75 × 2π is less probable. For $k = 2$, the likelihood is larger at phase angles of ~0.35 and 0.65 × π, and lower near the disc major axis, whereas that of the $k = 4$ mode is concentrated at a value of 0.5 × π/2, implying an orientation mainly lying at an equivalent separation from the disc major and minor axes. These systematic phase angles are robust against the number of titled rings and against the significance of the detected signal. Using a uniform number of rings for all galaxies lowers the significance of the result, but the trends for the shape of the probability distributions are preserved. More uniform likelihood distributions were expected for randomly distributed perturbations in the galaxies and with isotropic velocity ellipsoids for the H I gas. This is evidence that strong projection effects imprint on the velocity dispersion maps of H I gas in nearby galaxies.

Finding the origin(s) of the asymmetries and its systematics is a difficult task because of the nature of the velocity dispersion itself, which stems from instrumental effects and various processes inherent to the galactic random motions (dynamical, hydrodynamical, and thermal). In that aspect, we have shown that large-scale-induced BS might produce non-negligible residuals, especially in the fourth order for highly inclined galaxies with high projected rotation velocities. Typical signatures of BS on phase angles are orientated near the major axis ($k = 1$), the minor axis ($k = 2$), and at 0.5 × π/2 ($k = 4$). Nevertheless, these orientations are not fully consistent with the observations, especially for the second order, and the strength of Fourier modes in the BS contribution is much smaller than the observed asymmetries. This suggests that unresolved velocity streamings could be responsible for some variations in the likelihood distributions of the phase angles. The correlations found between Fourier decomposition of H I gas density and velocity dispersion maps suggest that the observed bisymmetry, and thus streaming motions, are induced by spiral arms, bars, or warps.

We thus further addressed the impact on the random motions from such velocity streamings present in the ordered motions $V_R$ and $V_{\theta}$. To do so, toy models mocking discs with a barely resolved spiral perturbation in both the density and kinematics were made, from which mock moment maps of discs seen under various projection angles, strengths, and orientations, of the perturbations were derived, and Fourier analyses were performed, by analogy to the observations. When streamings in the ordered motions are stronger in the radial dimension, hence in the presence of a radially biased velocity anisotropy on large scales, they