Investigating galactic fountains in M101

Aashiya Anitha Shaji, Françoise Combes, Anne-Laure Melchior, Anaëlle Hallé

► To cite this version:

Aashiya Anitha Shaji, Françoise Combes, Anne-Laure Melchior, Anaëlle Hallé. Investigating galactic fountains in M101. Astronomy & Astrophysics - A&A, 2025, 703, pp.A252. 10.1051/0004-6361/202556351. hal-05396258

HAL Id: hal-05396258

https://hal.science/hal-05396258v1

Submitted on 3 Dec 2025

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
# Investigating galactic fountains in M101
## Insights from UV emission, ionised gas, and neutral gas

Aashiya Anitha Shaji<fnref n="1" />,*, Françoise Combes<fnref n="1" />,<fnref n="2" />, Anne-Laure Melchior<fnref n="1" />, and Anaëlle Hallé<fnref n="1" />

Received 10 July 2025 / Accepted 30 September 2025

## ABSTRACT

Spiral galaxy disks are thought to exist in a quasi-stationary state, between fresh gas accretion from cosmic filaments and disk star formation, self-regulated through massive-star feedback. Our goal here is to quantify these feedback processes and probe their efficiency. While star formation can be traced at 10 Myr time-scales through Hα emission, the signature of OB stars, and at 100 Myr scale with UV emission, the gas surface density is traced by HI emission for the atomic phase. We chose to investigate feedback processes using fountain effects in M101, a nearby well-observed face-on galaxy. Face-on studies are complementary to the more frequent edge-on observations of these fountains in the literature. We use high-resolution data from THINGS for the HI emission, GALEX for UV, and SITELLE-SIGNALS for the Hα tracer. We identified 20 new HI holes, in addition to the 52 holes found in 1993. We study in more detail the nine holes satisfying strong criteria to be true fountain effects, compute their physical properties, and derive their energy balance. Only one small HI hole still contains Hα and young stars inside, while the largest hole of 2.4 kpc and the oldest age (94 Myr) is deprived of Hα and UV. For face-on disks, the possibility to study simultaneously the HI shell morphology, the stellar association, and kinematic evidence is of primordial importance. In M101, we quantified how stellar feedback is responsible for carving the observed cavities in the atomic gas disk, and how it can expel neutral gas above the disk, which is then unavailable for star formation during up to 100 Myr.

**Key words.** galaxies: evolution – galaxies: individual: M101 – galaxies: kinematics and dynamics – galaxies: spiral – galaxies: star formation

## 1. Introduction

Feedback processes are pivotal regulators of galaxy evolution, governing the baryon cycle between stars, the interstellar medium (ISM), and the circumgalactic medium. The interplay of gas accretion, outflow, and recycling controls star formation rates (SFRs), disk stability, and chemical enrichment (e.g. Somerville & Davé 2015). Stellar feedback, primarily from supernovae (SNe), stellar winds, and intense radiation fields, can expel metal-rich gas from the disk. When this gas cools and condenses after ejection, it can fall back and fuel subsequent star formation, forming a closed-loop mechanism known as the galactic fountain (Shapiro & Field 1976).

In the galactic fountain framework, clustered SN explosions and massive OB stellar winds heat the ISM to temperatures of order ~$10^6$ K, launching gas vertically (upto 1 kpc, Shapiro & Field 1976) at velocities that may or may not exceed the galaxy's escape speed. If the velocity is insufficient for escape, the ejected gas cools radiatively or mixes with cooler ISM phases, recombines into atomic hydrogen (HI), and can later form molecular clouds (H$_2$) under favourable conditions (Bregman 1980). These clouds, often seen as high-velocity complexes such as those in the Milky Way (van Woerden et al. 1999) and M31 (Davies 1975; Westmeier et al. 2007), eventually rain back onto the disk with velocities ~100 km s$^{-1}$ (Shapiro & Field 1976) over timescales of ~80–100 Myr (33 Myr for cooling and another 47 Myr to return to the disk, Houck & Bregman

1990). This cycle contributes to disk turbulence, redistributes metals, and regulates star formation over galactic scales (Armillotta et al. 2016).

Galactic fountains have been observed most directly in edge-on galaxies such as NGC 891 ($i \ge 89^\circ$; Oosterloo et al. 2007) and NGC 2403 ($i \sim 60^\circ$; Fraternali et al. 2002), where extra-planar HI and ionised gas trace vertical fountain flows above the plane. However, nearly face-on systems such as NGC 628 offer a complementary view. In such galaxies, vertical motions are largely hidden but the in-plane manifestations – HI holes, expanding shells, and luminous clusters in ultraviolet (UV) – can be identified in the disk as the launch sites of fountain activity (Kamphuis & Briggs 1992).

M101 (NGC 5457) is a nearby face-on spiral galaxy at a distance of 6.4 ± 0.7 Mpc (Shappee & Stanek 2011), with a modest inclination of 18° (Walter et al. 2008, see Table 1 for other properties). It hosts vigorous star formation (~2.4 $M_\odot$ yr$^{-1}$; Leroy et al. 2013), a flocculent barred structure, and no prominent bulge. At its distance, 1" corresponds to 31 pc see Table 1. Early HI studies by (Kamphuis 1993, henceforth K93) on this galaxy, using ~15" resolution Westerbork data, identified 52 HI holes and two high-velocity complexes. One of these holes, located on the outer disk, was interpreted as a superbubble requiring ~1000 SNe to account for its energy (Kamphuis et al. 1991). More recently, Chakraborti & Ray (2011) discovered a distinct supershell in the inner disk of M101, traced by UV emission, likely powered by clustered SNe. These studies established M101 as a promising case for investigating stellar feedback and the resulting ISM structures.

<footnote n="1">¹ Observatoire de Paris, LUX, CNRS, PSL University, Sorbonne University, 75014 Paris, France</footnote>

<footnote n="2">² Collège de France, 11 place Marcelin Berthelot, 75231 Paris, France</footnote>

<footnote>* Corresponding author: aashiya.shaji@obspm.fr</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Multi-wavelength views of NGC 5457 or M101. Left: False-colour RGB composite of M101 using H1 (THINGS, red), Hα (SITELLE, green), and FUV (GALEX, blue). The colour mapping is subtractive rather than additive, highlighting variations in neutral gas phase (H1), recently star-forming regions (Hα), and older star formation (FUV). Right: False-colour RGB composite of M101 using SITELLE's SN1, SN2, and SN3 filters. The mosaic combines four pointings (each covering 11' × 11') to capture the full extent of the galaxy.</figcaption></figure>

<table><caption>**Table 1.** Properties of galaxy M101 (NGC 5457) from the literature.</caption><thead><tr><th>Galaxy Properties</th><th>Values</th><th>References</th></tr></thead><tbody><tr><td>Distance (Mpc)</td><td>6.4</td><td>1</td></tr><tr><td>Inclination (deg)</td><td>18</td><td>2</td></tr><tr><td>Position Angle (deg)</td><td>39</td><td>3</td></tr><tr><td>r<sub>25</sub> (arcmin)</td><td>11.99</td><td>4</td></tr><tr><td>Scale length, r<sub>e</sub><sup>gas</sup> (kpc)</td><td>3.57</td><td>5</td></tr><tr><td>Scale length, r<sub>e</sub><sup>star</sup> (kpc)</td><td>3.89</td><td>6</td></tr><tr><td>H1 Mass (M<sub>☉</sub>)</td><td>14.17 × 10<sup>9</sup></td><td>7</td></tr><tr><td>SFR (M<sub>☉</sub> yr<sup>-1</sup>)</td><td>2.4</td><td>8</td></tr></tbody></table>

**Notes.** (1) Shappee & Stanek (2011); (2) Walter et al. (2008); (3) Bottema (1993); (4) Schruba et al. (2011); (5) Schruba et al. (2011); (6) Casasola et al. (2017); (7) Walter et al. (2008); (8) Leroy et al. (2013).

While UV, H1 and ionised gas observations of M101 have each been studied independently, no systematic multi-wavelength attempt has been made to identify galactic fountain activity and form a complete view of the feedback loop. This study presents the first such effort: using high-resolution data from H1 21 cm, UV, and optical integral field unit (IFU) spectroscopy, we aim to identify observational signatures of galactic fountains in M101 from a face-on perspective.

The features we look for as signatures of galactic fountain activity span multiple wavelengths and gas phases. In the ultraviolet regime, we expect to detect compact star-forming complexes whose massive stars are capable of driving fountains via their stellar winds and SN explosions. These appear as clumps in the GALEX FUV intensity maps and trace recent massive star formation on ≲100 Myr timescales, compared to ≲10 Myr for Hα emission (Keel et al. 2004; Rampazzo et al. 2022). The UV thus

provides a longer-baseline locator for potential fountain launch sites.

The SNe from these associations inject energy into the surrounding interstellar medium, driving outflows of hot, ionised gas. These can be traced with SITELLE kinematics via velocity-broadened or multi-peaked Hα profiles and with coherent morphologies such as shells and bubble-like structures in narrow-band maps (Camps-Fariña et al. 2017; Rupke 2018, and references therein). Recent PHANGS studies provide high-resolution, multi-phase benchmarks for feedback-driven bubbles: from molecular-gas measurements of sizes, expansion speeds, energetics and lifetimes (Watkins et al. 2023a), to energy-balance analyses linking the mechanical output of young massive stars with the kinetic energy of turbulent, ionised gas (Egorov et al. 2023), and a panchromatic, spatially resolved census of bubbles and their associated stellar populations in NGC 628 (Watkins et al. 2023b).

These energetic events also displace the atomic hydrogen that previously occupied the region. As a result, observations in the 21 cm line of H1 often reveal these sites as gaps, depressions, or cavities in the atomic gas distribution. Such H1 holes are commonly seen in nearby spiral galaxies (Brinks & Bajaja 1986; Deul & den Hartog 1990; Puche et al. 1992; Kamphuis 1993; Wilcots & Miller 1998; Fraternali et al. 2002; Boomsma et al. 2008; Bagetakos et al. 2011; Pokhrel et al. 2020), including the Milky Way (Heiles 1979, 1984; Hu 1981). While some of these holes appear to be genuine voids, others are known to be filled with either hot ionised gas or colder molecular material. Their physical origin can vary, but many are linked to disk processes such as clustered SN explosions, stellar winds, or, in some cases, gas infall from the halo or circumgalactic medium (e.g. Putman et al. 2012). A recent work by Lara-López et al. (2023) concluded that supernovae remnants (SNRs) concentrate on the rims of H1 holes and can trigger star formation there.
Through this study, we aim to determine which disk cavities in M101 are consistent with feedback-driven fountains and to assess their sizes, possible ages, and implications for the cycling of gas between the disk and halo. We describe our data sources in Section 2 and methodology in Section 3. Our results are presented in Section 4, with their implications discussed in Section 5.

## 2. Data description

This study makes use of three archival, multi-wavelength datasets to trace young stellar populations, atomic hydrogen structures, and ionised gas morphology and kinematics, respectively. Figure 1 is a hybrid map combining atomic gas, ionised gas, and ultraviolet emission.

### 2.1. GALEX FUV and NUV

The Galaxy Evolution Explorer (GALEX, Martin et al. 2005) data used in this study were obtained from the Guest Investigator Cycle 3 with ID: G13_050008_NGC5457 in the far-UV (FUV: 1500–1800 Å) and near-UV (NUV: 1750–2800 Å) bands, with exposure lengths of 13293.4 and 13294.4 seconds, respectively. These data were retrieved from the MAST Archive and have spatial resolutions of 4.3" (FUV) and 5.3" (NUV). For this analysis, we use the intensity maps (fd-int.fits.gz and nd-int.fits.gz), which are provided in units of counts s⁻<fnref n="1" /> pixel⁻<fnref n="1" /> with a pixel scale of 1.5". These images are used primarily for morphological identification of UV-bright stellar associations near H I holes and SITELLE shells.

### 2.2. The H I nearby galaxy survey

The H I Nearby Galaxy Survey (THINGS, Walter et al. 2008), conducted using the NRAO Very Large Array, provides high-resolution 21 cm H I observations of nearby galaxies. We retrieved the data for M101 (NGC 5457) from the THINGS public archive<fnref n="1" />, along with its moment maps for integrated flux, velocity field, and velocity dispersion.

For our analysis, we used the data cube with natural weighting, which provides low noise and has a spectral resolution of 5.2 km s⁻<fnref n="1" /> with a total of 69 velocity channels. We also used the maps obtained from the cube of robust weighting to cross-check hole boundaries and structure in ambiguous cases, as it has higher resolution but also higher noise levels. The image dimensions are 2048 × 2048 pixels, with a pixel size of 1.0". The root-mean-squared (RMS) noise in one channel map is 0.46 mJy beam⁻<fnref n="1" />, which corresponds to an H I column density sensitivity of ~3.2 × 10²⁰ cm⁻². The synthesised beam size of the natural (respectively, robust) cube is 10.82" × 10.17" (7.49" × 6.07") with a position angle of -67.0° (-56.6° measured north to east; Walter et al. 2008).

Table 2 compares the observational properties of the original K93 dataset with the THINGS data used in this study.

### 2.3. SIGNALS-SITELLE

M101 is one of the galaxies observed in the Star-formation, Ionised Gas, and Nebular Abundances Legacy Survey (SINGS; Rousseau-Nepton et al. 2019), conducted at the Canada-France-Hawaii Telescope (CFHT) using the Spectromètre Imageur à Transformée de Fourier pour l'Etude en Long et en

<table><caption>**Table 2.** Comparison of H I observational properties.</caption><thead><tr><th>Properties</th><th>K93</th><th>THINGS</th></tr></thead><tbody><tr><td>Channel spacing (km s<sup>-1</sup>)</td><td>4.1</td><td>5.2</td></tr><tr><td>No. of Channels</td><td>127</td><td>69</td></tr><tr><td>Bandwidth (MHz)</td><td>2.5</td><td>1.56</td></tr><tr><td>Spatial Resolution</td><td>13" × 16"</td><td>7.49" × 6.07"</td></tr><tr><td>RMS Noise (mJy/beam)</td><td>0.5</td><td>0.46</td></tr><tr><td>Exposure time</td><td>16 × 12 h</td><td>10 h</td></tr></tbody></table>

Large de raies d'Émission (SITELLE) imaging spectrograph (Drissen et al. 2019). Because of its proximity and extent, the galaxy was observed across four fields with three spectral filters: SN1 (363–386 nm), SN2 (482–513 nm), and SN3 (647–685 nm), at spectral resolutions of 1000 for SN1 and SN2, and 5000 for SN3, where 'SN' stands for SITELLE Narrow-band. The reduced data, retrieved from the Canadian Astronomy Data Centre, is available as HDF5 files for each filter and field.

The emission lines covered comprise [O II] λλ3726, 3729 (SN1), Hβ at 4861 Å and [O III] λλ4959, 5007 (SN2), and Hα at 6563 Å, [N II] λλ6548, 6584, and [S II] λλ6717, 6731 (SN3). The deep frame, which comes with the data cube and is the sum of all individual interferometric frames, provides a high signal-to-noise image of the field per filter. Amosaic using these narrowband filters in the optical regime is shown in the right panel of Fig. 1.

The instrumental line shape of SITELLE is a cardinal sine (sinc) function, and the physical information of an emitting region is embedded in a sinc-Gauss (the result of the convolution of a sinc with a Gaussian) profile (Martin et al. 2016). We therefore analyse the SITELLE data cubes using *Outils de Réduction de Cubes Spectraux* (ORCS, Martin et al. 2015), a software designed for SITELLE data reduction and analysis. In this work, we used the ORCS-generated maps of flux, velocity and dispersion, along with their uncertainties, for the Hα and [S II] lines, which were produced after background subtraction and simultaneous fitting of a single sinc-Gauss profile to all emission lines in the SN3 data cube without any pixel binning. Accordingly, when referring to 'Hα velocity' or 'Hα dispersion' throughout this paper, we mean the velocity and dispersion maps derived from the global SN3 fit.

## 3. Methodology

In this section, we describe how regions hosting fountain-driven feedback features are identified using a multi-wavelength analysis of H I, ultraviolet, and ionised gas data. Following Bagetakos et al. (2011) and Pokhrel et al. (2020), we forgo automated H I-hole detection in favour of visual inspection owing to the irregular morphologies of H I structures in galaxy disks. We identified candidate galactic fountain sites in M101 via two approaches mentioned in Section 3.1, then subjected them to a scoring scheme in Section 3.2 to identify true holes driven by feedback, and finally estimated certain observed and derived properties as mentioned in Section 3.3 for each of these high-scoring sites hosting a galactic fountain.

### 3.1. Preliminary candidates

As a preliminary step, we constructed an initial catalogue of H I cavities in M101 using two complementary routes, before any

<footnote n="1">¹ https://www2.mpcia-hd.mpg.de/THINGS/Data.html</footnote>
assessment of whether these are bona fide feedback-driven holes (Section 3.2). Our goal at this stage was to identify all plausible cavity candidates visible in the H I moment-zero map, regardless of their eventual classification.

First, we revisited the 52 H I holes catalogued by K93 whose diameters range from ~0.8 to 6.6 kpc. These were originally detected using data with ~15" resolution, limiting sensitivity to structures larger than ~600 pc. With the improved resolution of the THINGS dataset (6", corresponding to ~200 pc at M101's distance of 6.4±0.7 Mpc; Shappee & Stanek 2011), we revisited these holes using higher-quality H I, SITELLE Hα, and GALEX UV data. To do so, we overlaid the positions and extents of the K93 holes as ellipses or circles on the THINGS H I moment maps, SITELLE Hα emission maps, and GALEX FUV intensity maps using CARTA (Comrie et al. 2021). Out of these 52 holes, one hole lies outside the THINGS field-of-view, so we drop it early on. The rest of these 51 structures are then subjected to a scoring scheme described in Section 3.2.

In the second approach, we started with a visual inspection of the SITELLE colour image given in the right panel of Fig. 1, which was produced from the individual SN1, SN2 and SN3 deep frames using the IRAF rgbsun routine (Tody 1986), with final adjustments made in GIMP<fnref n="2" /> by Laurent Drissen (private communication). Owing to the stratified, onion-like ionisation structure of photoionised nebulae (see e.g. Morisset 2018), and to the fact that each SN filter samples different ionisation lines, spatially resolved SITELLE-IFU data reveal distinct ionisation zones in such an RGB composite: low-ionisation lines dominate the outer rims, while higher-ionisation lines appear towards the interior. We exploit this to visually identify circular or elliptical rim- or bubble-like H II regions in the RGB image. Each candidate is then cross-matched with the GALEX FUV map to confirm the presence of a young stellar cluster inside or near the centre. Finally, we inspected the THINGS H I natural and robust moment-zero map at that location for a corresponding local minimum. This process yielded 20 new, generally smaller, cavity candidates. This approach also found back some of the previously-catalogued holes from K93.

Combining the two approaches, we obtained a preliminary set of 71 H I cavities. These are all subjected, in the next stage, to the uniform scoring procedure described in Section 3.2 to assess whether they are genuine holes and whether their properties are consistent with galactic fountain-driven origins. We also generated [S II]/Hα maps, as shock-ionised shells from SNR often exhibit elevated [S II]/Hα > 0.5.

## 3.2. Quality assurance scoring scheme

The scoring schemes of Bagetakos et al. (2011) and Pokhrel et al. (2020) were designed to identify H I holes, without explicitly distinguishing between those created by stellar feedback and those of dynamical origin. In this work, we adapt and extend their approach to ensure that only feedback-driven holes are selected, with additional criteria tailored to identifying sites of galactic fountain activity. Each hole is assigned points based on its morphological, kinematic, and multi-wavelength signatures, as depicted in Figure 2.

One point is assigned if the hole shows a well-defined flux minimum in the H I integrated flux map – quantified as a ≥50% drop in surface brightness in the line profile extracted along the position angle of the ellipse – suggesting coherent depletion (Panels a and b of Fig. 2). Another point is added if this depres-

sion remains spatially consistent across at least three consecutive channels in the THINGS data cube, helping to mitigate false positives caused by turbulence or noise. Additional points are given depending on the number of contiguous velocity-channels where the depression is seen: 1–2 (5–10 km s⁻¹), 3–4 (11–20 km s⁻¹), and >4 (>20 km s⁻¹) channels receive +1, +2, and +3 points respectively. An extra point is added if the hole shows an elevated velocity dispersion in the H I moment-2 map compared to the local surroundings (Panel c in Figure 2). If the hole lies completely within a spiral arm, it receives +1, as spiral-arm holes are more likely to be star-formation-driven.

Next, we extracted position–velocity (PV) diagrams along the spiral arm direction for each hole using pvextractor (Ginsburg et al. 2016). In ambiguous cases, we rotated the cut across multiple angles (from 0° to 180°) to optimise visibility. Based on the classification scheme of Brinks & Bajaja (1986, henceforth referred to as Brinks type), the holes are classified into three types as illustrated in row h Figure 2. A type I is known as a fully blown-out hole with neither receding nor approaching components visible on either side of the PV diagram (top panel). In type II holes, gas deviation is observable on one side of the PV diagram (middle panel). Such cavities are called partially blown-out holes. The third type is termed intact holes. Here, both receding and approaching sides are visible, forming an elliptical signature in PV space (bottom panel). Type III holes generally evolve into larger type I holes over time. We assigned +1 point for sharp PV edges in type I holes. Type II holes receive +1 if a well-defined receding or approaching ‘bump’ is visible. If this bump spans more than two channels, an additional +2 is given. This extended bump criterion, although rarely observed, was included for completeness.

To evaluate associated star formation, we assessed UV (Panel e of Fig. 2) and Hα (Panel f of Fig. 2) emission in and around each hole. If a UV-bright clump lies inside the H I hole, we assigned +2 points; if it lies near the edge, +1 point. We assigned +2 if an Hα clump lies inside, +1 if near the edge. Additionally, if the Hα emission within the hole showed higher velocity dispersion (Panel g of Fig. 2) than its surroundings, we gave it an extra point. Two extra points were given if the hole’s edge contains bright rims or arcs (bubble-like structure) in the false-colour RGB image created using SITELLE deep frames (right of Fig. 1), when scoring the holes found through approach one.

Each hole is thus scored across multiple observational dimensions: H I morphology, kinematics, ionised gas features, and embedded stellar populations, out of a total of 20 points. Any structures with a quality score, Q, less than 10 points have been discarded, whereas those with more than 15 points are considered highly likely to host galactic fountains. We applied this score to all candidate holes; all 20 newly discovered holes have crossed the threshold of Q ≥ 10, hence we compare them with the earlier K93 catalogue in Section 5.1. Those scoring Q ≥ 15 are then characterised in detail through the following observed and derived parameters.

## 3.3. Estimation of properties

### 3.3.1. Observable properties

For each shell with points above 15, we record: central coordinates (RA, Dec; ±3"), heliocentric velocity $v_{hel}$ (channel uncertainty ±5.2 km s⁻¹), semi-major or minor axes $b_{maj}$, $b_{min}$, axial ratio $b_{min}/b_{maj}$, position angle (±30°), Brinks type, expansion velocity $v_{exp}$ (see below), and mean H I flux density ($S_{avg}$).

<footnote n="2">² GNU Image Manipulation Program.</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Example of our multi-wavelength quality assurance scoring for candidate H I holes. Panels a–g: Each panel corresponds to the same hole (Hole 1) with the red ellipse marking its boundary and annotations in light blue giving the score for each feature. (a) H I moment-0 map showing a cavity. (b) H I line profile along the blue dashed line in panel a, used to quantify a 50% surface brightness depression. (c) H I velocity-dispersion (moment-2) map highlighting locally elevated $\sigma_{HI}$. (d1–d3) Individual H I channel map at representative velocities; the hole is visible in >4 consecutive channels. (e) GALEX FUV map with bright emission coincident with the hole. (f) H$\alpha$ emission map with strong flux within the hole. (g) H$\alpha$ velocity-dispersion map with elevated $\sigma_{H\alpha}$ inside the hole. Bottom row: Position–velocity (PV) diagrams illustrating the three morphological classes of holes following Brinks & Bajaja (1986). The PV panels here (left to right) show Hole 3 (Type I), Hole 2 (Type II), and Hole 1 (Type III). Hole 2 also suffers from a tidal interaction, see Section 4.</figcaption></figure>

We take $v_{\text{hel}}$ to be the velocity channel at which the hole is clearly visible (Panel d2 of Fig. 2). The calculation of expansion velocity depends on the type of hole. For type I holes, $v_{\text{exp}}$ can-

not be determined, but its upper limit can be deduced from the average velocity dispersion within the hole. For partially blown-out holes, $v_{\text{exp}} = v_{\text{hel}} - v_{\text{bump}}$ where $v_{\text{bump}}$ is the velocity at the
<table><caption>**Table 3.** Physical properties observed in the nine H I holes that are sites hosting galactic fountains.</caption><thead><tr><th>Name</th><th>RA<br>(hh:mm:ss)</th><th>Dec<br>(dd:mm:ss)</th><th>b<sub>maj</sub><br>(arcsec)</th><th>b<sub>min</sub><br>(arcsec)</th><th>PA<br>(degrees)</th><th>b<sub>min</sub>/b<sub>min</sub></th><th>S<sub>avg</sub>Δv<br>(Jy/beam * m/s)</th><th>Brinks<br>(Type)</th><th>v<sub>exp</sub><br>(km s<sup>-1</sup>)</th><th>v<sub>hel</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>Hole 1</td><td>14:03:54.84</td><td>+54:21:41.78</td><td>34</td><td>30</td><td>37</td><td>0.9</td><td>86.7</td><td>1</td><td>54</td><td>275</td></tr><tr><td>Hole 2</td><td>14:03:41.69</td><td>+54:29:38.29</td><td>33</td><td>25</td><td>65</td><td>0.8</td><td>54.9</td><td>2</td><td>41</td><td>296</td></tr><tr><td>Hole 3</td><td>14:03:27.13</td><td>+54:17:35.74</td><td>45</td><td>33</td><td>0</td><td>0.7</td><td>86.2</td><td>3</td><td>13</td><td>224</td></tr><tr><td>Hole 4</td><td>14:02:25.19</td><td>+54:27:44.45</td><td>30</td><td>23</td><td>84</td><td>0.8</td><td>132.0</td><td>3</td><td>28</td><td>255</td></tr><tr><td>Hole 5</td><td>14:03:33.88</td><td>+54:18:29.59</td><td>9</td><td>8</td><td>90</td><td>0.9</td><td>91.3</td><td>0</td><td>9</td><td>248</td></tr><tr><td>Hole 6</td><td>14:02:47.53</td><td>+54:14:54.31</td><td>17</td><td>6</td><td>90</td><td>0.4</td><td>132.4</td><td>0</td><td>10</td><td>183</td></tr><tr><td>Hole 7</td><td>14:02:17.63</td><td>+54:18:31.60</td><td>13</td><td>8</td><td>90</td><td>0.6</td><td>167.8</td><td>0</td><td>8</td><td>183</td></tr><tr><td>Hole 8</td><td>14:02:37.50</td><td>+54:28:57.82</td><td>20</td><td>7</td><td>259</td><td>0.4</td><td>96.1</td><td>0</td><td>7</td><td>265</td></tr><tr><td>Hole 9</td><td>14:02:14.72</td><td>+54:25:47.95</td><td>12</td><td>6</td><td>90</td><td>0.5</td><td>57.2</td><td>0</td><td>9</td><td>235</td></tr></tbody></table>

Notes. Brinks type of the smaller holes has been set to zero, considering their position-velocity diagrams are too ambiguous.

<table><caption>**Table 4.** Derived properties of the nine H I holes that are sites hosting galactic fountains.</caption><thead><tr><th>Name</th><th>Diameter, d<br>(kpc)</th><th>Galactocentric<br>Distance, R<br>(kpc)</th><th>Kinetic<br>Age, t<sub>kin</sub><br>(Myr)</th><th>N(HI)<br>(10<sup>20</sup> cm<sup>-2</sup>)</th><th>Effective<br>thickness, l(r)<br>(pc)</th><th>E<sub>Mc</sub><br>(erg)</th><th>E<sub>Ch</sub><br>(erg)</th><th>Mass<br>Displaced, M<sub>HI</sub><br>(M<sub>☉</sub>)</th></tr></thead><tbody><tr><td>Hole 1</td><td>2.03</td><td>12.76</td><td>18.8</td><td>8.6</td><td>31.1</td><td>8.37e+56</td><td>3.96e+56</td><td>9.71e+08</td></tr><tr><td>Hole 2</td><td>1.80</td><td>19.45</td><td>21.9</td><td>5.5</td><td>13.2</td><td>4.31e+56</td><td>2.91e+56</td><td>1.01e+09</td></tr><tr><td>Hole 3</td><td>2.43</td><td>8.33</td><td>93.6</td><td>8.6</td><td>55.0</td><td>9.47e+54</td><td>5.00e+55</td><td>9.44e+08</td></tr><tr><td>Hole 4</td><td>1.67</td><td>20.43</td><td>29.9</td><td>13.2</td><td>11.6</td><td>3.24e+56</td><td>4.19e+56</td><td>2.22e+09</td></tr><tr><td>Hole 5</td><td>0.55</td><td>8.35</td><td>30.6</td><td>9.1</td><td>54.9</td><td>1.71e+53</td><td>3.10e+53</td><td>1.16e+07</td></tr><tr><td>Hole 6</td><td>0.68</td><td>14.19</td><td>33.8</td><td>13.2</td><td>25.9</td><td>1.09e+54</td><td>2.40e+54</td><td>6.62e+07</td></tr><tr><td>Hole 7</td><td>0.69</td><td>17.07</td><td>42.8</td><td>16.7</td><td>17.9</td><td>1.05e+54</td><td>3.61e+54</td><td>1.26e+08</td></tr><tr><td>Hole 8</td><td>0.78</td><td>19.97</td><td>55.8</td><td>9.6</td><td>12.3</td><td>7.57e+53</td><td>3.66e+54</td><td>1.55e+08</td></tr><tr><td>Hole 9</td><td>0.56</td><td>20.51</td><td>31.2</td><td>5.7</td><td>11.5</td><td>5.31e+53</td><td>1.13e+54</td><td>3.69e+07</td></tr></tbody></table>

bump (where the gas is deviating from the surroundings, rightmost panel in row h of Fig. 2). For the intact holes, we use the average difference between $v_{hel}$ and the velocities of the gas on the approaching and receding sides of the hole. The uncertainty of the calculation is again the velocity resolution (5.2 km s⁻¹). We determined the $S_{avg}$ by averaging the flux values on the integrated H I map over a region of size twice the $b_{maj}$ and $b_{min}$ of that of the hole itself. The errors here are estimated to be of the order 10%, arising from the uncertainty about the position of the hole centre.

### 3.3.2. Derived properties

The derived properties are as follows. Most of the following equations, if not specifically cited, are adopted from the past H I hole studies conducted by Bagetakos et al. (2011) and Pokhrel et al. (2020):

- Diameter, $d = 2\sqrt{b_{maj}b_{min}}$.

- The galactocentric distance in pc of a hole with centre $(\alpha, \delta)$ to the galaxy centre $(\alpha_0, \delta_0)$ in radians, can be obtained using the distance to the galaxy (D) in pc, position angle ($\theta$), and inclination ($i$):

$$R = D \sqrt{(x'')^2 + (y'')^2}, \qquad (1)$$

where $x'' = x \sin \theta + y \cos \theta$, $y'' = \frac{y \cos \theta - x \cos \theta}{\cos i}$,

and $x = (\alpha - \alpha_0) \cos \delta_0$, $y = \delta - \delta_0$.

- The kinetic age in Myr can be calculated from the diameter in parsec and the expansion velocity in km s⁻¹,

$$t_{\text{kin}} = 0.978 \frac{d/2}{v_{\text{exp}}}. \qquad (2)$$

- H I column density,

$$N_{\text{HI}} = 1.823 \times 10^{18} \sum_{i} \left[ \frac{S_{\text{avg}}^{i}}{1.66 \times 10^{-3} B_{\text{maj}} B_{\text{min}}} \cdot \Delta v \right]. \qquad (3)$$

Here, $S_{avg}^i$ is the mean H I flux density around the hole in mJy/beam, in each velocity channel $i$. $B_{maj}$ and $B_{min}$ are beam sizes in arcsec, $\Delta v$ channel width in km s⁻¹.

- The effective thickness, $l(r)$, of the atomic hydrogen disk is calculated from the inclination of the galaxy, $i$, as:

$$l(r)[\text{pc}] = \frac{Z_0 \sqrt{2\pi}}{\cos i}. \qquad (4)$$

Here, the scale height $Z_0$ is derived from the H I velocity dispersion (≈12 km s⁻¹), assuming hydrostatic equilibrium of the atomic gas in the gravitational potential of the stellar disk, with negligible dark matter (Combes & Becquaert 1997):

$$Z_0[\text{pc}] = \frac{v_{\text{disp}}}{\sqrt{4\pi G\rho(r)}}, \qquad (5)$$

$$\text{where } \rho(r) = \frac{\Sigma^*(r)}{h_z^*} = \frac{\Sigma_0^* e^{-r/r_e^*}}{h_z^*}, \qquad (6)$$
is the volumic stellar density, $\Sigma^*(r)$ the stellar surface density, $r_e^*$ the stellar exponential radial scale and $h_z^*$ is the constant scale height of stellar disks (e.g. Bottema 1993). We adopted the stellar exponential radial scale from Casasola et al. (2017), which utilised the IRAC data to calculate the stellar scale length of M101 to be 2.09 arcmin $\approx 3.89$ kpc. We deduced $\Sigma_0^*$ to be 381 M$_\odot$ pc$^{-2}$ from their study. We also considered $h_z^* = r_e/7.3 \approx 0.53$ kpc following Leroy et al. (2008), who adopted the average flattening ratio of 7.3 from Kregel et al. (2002). This is based on the assumption that the height of the stellar disk is independent of its radius (van der Kruit & Searle 1982).

- The energy required to power the hole can be calculated using two different formulae. First, using Chevalier (1974), which is more dependent on the diameter of the hole, and the second using McCray & Kafatos (1987), which has high dependence on the expansion velocity.

$$E_{\text{Ch}}[\text{erg}] = 5.3 \times 10^{43} n_0^{1.12} [\text{cm}^{-3}] \left(\frac{d[\text{pc}]}{2}\right)^{3.12} v_{\text{exp}}^{1.4} [\text{kms}^{-1}], \quad (7)$$

$$E_{\text{Mc}}[\text{erg}] = n_0 \left(\frac{d[\text{pc}]}{194}\right)^2 \left(\frac{v_{\text{exp}}}{5.7}\right)^3 \times 10^{51}. \quad (8)$$

Here we took the density of the ambient medium, $n_0$, as the particle density of the atomic hydrogen, $n_{\text{HI}}$, by ignoring contributions from He and H$_2$. The midplane H I volume density ($n_{\text{HI}}$) is calculated using the column density as follows:

$$n_{\text{HI}} = \frac{N_{\text{HI}}}{3.08 \times 10^{18} \ell(r)}. \quad (9)$$

- H I mass of the hole:

$$M_{\text{HI}}[M_{\odot}] = 0.0245 n_{\text{HI}} V, \quad (10)$$

where the volume of the hole is given by

$$V[\text{pc}^3] = \frac{4}{3}\pi\left(\frac{d}{2}\right)^3. \quad (11)$$

## 4. Results

Out of the 72 H I holes surveyed (52 from K93 plus 20 newly identified cavities), nine satisfied our quality threshold ($Q \ge 15$) and exhibit the most compelling multi-wavelength signatures of galactic fountains. Figure 3 maps these nine high-confidence sites onto the H I moment-0 image of M101 as red ellipses. Four of the 52 ‘classic’ K93 holes (dashed, red ellipses) and five of the new bubbles (solid, red ellipses) pass our scoring cut, while holes with intermediate scores ($10 < Q < 15$) are shown in purple for completeness.

Tables 3 and 4 list the observed and derived properties of these nine holes, respectively. Briefly, their diameters span 0.5–2.4 kpc (median 0.7 kpc), expansion velocities run from 7 to 54 km s$^{-1}$ (median 10 km s$^{-1}$), and kinetic ages range from 18 to 93 Myr (median ~31 Myr). The implied H I masses displaced were $1 \times 10^7 - 2 \times 10^9 M_\odot$ (median: $1.5 \times 10^8 M_\odot$), and the required energies (via Chevalier’s formula) lie between $10^{53}$ and $10^{56}$ erg (median: $3.5 \times 10^{54}$ erg). These best-score holes tend to coincide with prominent spiral-arm H II complexes (e.g. NGC 5461, 5462, 5455), confirming that fountain-driving feedback is concentrated in the disk’s active star-forming regions.

Our two discovery routes yield a complementary population. The four largest cavities (diameters 1.6–2.4 kpc; Figure 4)

<figure><img src="image_4.png" /><figcaption>Fig. 3. H I moment-0 map of NGC 5457 (M101), overlaid with the locations of highly likely sites of galactic fountains. Regions are categorized based on quality flags (Q): red dashed ellipses mark holes from K93 with $Q \ge 15$, red solid ellipses show newly identified holes with $Q \ge 15$, purple dashed ellipses correspond to K93 holes with $10 < Q < 15$, and purple solid ellipses mark new detections with $10 < Q < 15$. The background blue-scale image shows the integrated H I emission, with brightness scaled between the 10th and 98th percentiles. The nine regions in red are the identified sites of galactic fountains.</figcaption></figure>

were all catalogued by K93. In contrast, the five smaller holes (0.5–0.7 kpc; Figure A.1) emerged only from our SITELLE-first approach. This hierarchy of scales demonstrates that high-resolution multi-wavelength data can reveal more modest but still fountain-driven structures that previous H I surveys missed. For these smaller bubbles, we could not reliably assign Brinks types via PV diagrams; instead, we adopted the local H I velocity dispersion as an upper limit on the expansion speed.

Hole 3 is the largest cavity (2.4 kpc diameter) and also the oldest (~93 Myr), located on the south-eastern arm adjacent to NGC 5461. Its PV diagram shows a fully “blown-out” type I signature, with no coherent shell walls ($v_{\text{exp}} \le 13$ km s$^{-1}$). The lack of interior H$\alpha$ or FUV emission is consistent with its advanced age (the progenitor OB stars appear to have died). The filamentary clumps at its rim and a two-component fit to the H$\alpha$ line (whose centroid separation matches $v_{\text{exp}}$), point to triggered star formation along the shell edge. This is also adjacent to where the 2023 SN went off.

Hole 1 is the second largest cavity (2.03 kpc diameter) and the largest expansion velocity of 54 km s$^{-1}$. It was already discovered by Kamphuis et al. (1991) as the first clear case of an expanding H I shell associated to an H I hole, corresponding to the energy of a thousand SNe. Situated on the north-eastern arm, Hole 2 exhibits the second largest expansion velocities, ~40 km s$^{-1}$, and was one of the holes inferred from both approaches. Although bubble morphology is evident on the rim of this hole with SITELLE RGB, the H$\alpha$ line does not decompose cleanly into two components, suggesting geometric projection effects or rapid shell fragmentation at this evolutionary stage. Considering the exceptionally high velocity dispersion in this hole (cf second row and second column in Figure 4, and
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Multi-wavelength views of the four biggest sites of galactic fountains in M101 marked by red ellipses. From top: H I velocity map overlayed with its 5th and 95th flux contours in black; H I velocity dispersion; FUV flux map; Hα flux map; Hα velocity dispersion map; the ratio of [S II] to Hα flux overlayed with the positions of optically identified SNRs by Matonick & Fesen (1997) in black symbols.</figcaption></figure>

the middle panel in bottom row of Fig. 2), it is possible that this hole is suffering in addition from a tidal interaction from the M101 dwarf companion NGC 5477, as proposed by Combes (1991), or a collision with a large extra-galactic gas cloud as proposed by van der Hulst & Sancisi (1988), since it corresponds exactly to the perturbed region.

Among the small bubbles, Hole 6 stands out as the sole H I cavity that coincides precisely with an optical rim (Figure A.1). With a diameter of just 0.6 kpc, it shows elevated velocity dispersion in both H I moment-2 and Hα σ maps at its core, and hosts a compact FUV cluster at its centre.

# 5. Discussion

## 5.1. *20 New H* <u>*i*</u>*o*le*<u>*h*</u>*o*ls

Approach 2 yields 20 new and smaller holes that were missed by the 15" resolution of the K93 studies. Here, we discuss their size distribution, age, and shape using Figure 5. The new holes peaked at smaller radii (150–400 pc), whereas the original sample is biased towards radii ≥ 400 pc. This size-distribution shift directly reflects the factor-of-2 improvement in spatial resolution from THINGS (6") versus the earlier 15" WSRT data and demonstrates how higher resolution reveals a hierarchy of feedback structures.
<figure><img src="image_6.png" /><figcaption>Fig. 5. Comparison of properties between the 20 new H I holes in blue and the 52 K93 holes in grey. The medians of each sample are plotted as dashed lines. As expected, the new holes span smaller diameters and are younger than those of the K93 sample. However, these younger holes do not seem to have a more circular shape than the older ones.</figcaption></figure>

These new holes expand at only 5–17 km s⁻¹, far slower than the 25–45 km s⁻¹ typical of the K93 sample. Because the PV diagrams of sub-kpc cavities rarely show distinct approaching or receding bumps, we estimate each small hole’s expansion velocity simply as the H I moment-2 dispersion averaged over its area. This procedure naturally yields a lower $v_{exp}$ than for the K93 holes, due to the measurement method: tiny shells cannot be kinematically resolved as well as 1 kpc-scale bubbles. Consequently, when we computed the $t_{kin}$ (as we did for the high Q sites), they appeared systematically similar in age (median ~30 Myr) to their larger, faster K93 holes (~35 Myr). This is clear from Figure 6, which confirms that we have an observational bias of lower velocities for the newly identified holes. We note that the error bars on $v_{exp}$ for these small holes are much larger, since the estimated expansion velocity is very close to the H I typical velocity dispersion of 10 km s⁻¹. These large relative errors of up to ~100% prevent us from estimating their corresponding energy with precision.

Both samples show a wide range of axial ratios (0.4–1.0). Ideally, we expected that the differential rotation would shear holes into an elliptical shape with time. This would also mean that the shape of the bigger holes of the K93 sample should be more elliptical than the younger ones found in this work. Consequently, we compared the axial ratio and kinetic age of holes (Figure 7) and found no statistically significant correlation. This lack of correlation might be due to the low shear in the rotation curve of this late-type galaxy, with a very light bulge. The rotation is almost a solid-body one until a radius of 10 kpc (Bosma et al. 1981).

We also examined of these newly found holes as a function of the radius from the centre of the galaxy. There is a weak trend suggesting that the age of the H I hole increases with galactocentric distance. This indicates that inner-disk shells are, on average, younger, consistent with more rapid dissolution in regions of denser gas. However, the scatter is large, and the trend is not yet robust, as shown in Figure 8.

If we assume that all of these holes were indeed powered by type II SNe, we can plot their “expansion” velocities (which in our case is an upper limit) against their diameters as done in Figure 9 and reversely use Eq. (7) to determine the energy required to power them. We take the $n_0$ value to be 0.35 cm⁻³, which is the midplane density discussed by Kuntz et al. (2003).

<figure><img src="image_7.png" /><figcaption>Fig. 6. Expansion velocities of the new and old holes plotted against their kinetic age. The size of any symbol is proportional to the diameter of the H I hole it represents. Even if their sizes are bigger, the relatively higher expansion velocities of the K93 sample holes push them to lie within a younger age area of the plot. The slower expansion velocities, close to the local H I velocity dispersion (~10 km s⁻¹), render the kinetic age of the new sample very similar to that of the old.</figcaption></figure>

The following equation is based on the hydrodynamical description of the evolution of SNe modelled by Chevalier (1974), where they related the time with the radius (R) of the holes and expansion velocity $v_{exp}$ to be:

$$t[\text{Myr}] = 0.304R/v_{\text{exp}}. \qquad (12)$$

We use this to draw constant time and energy lines on Figure 9 as was done by Puche et al. (1992) for the nearby dwarf galaxy Holmberg II. As in their case, this equation gives ages 3 times lower than the kinetic ages we calculate. Furthermore, it is important to note here that these estimates of energy are not necessarily related to the energy estimated in Table 4 to power the four largest fountain-sites.

As a final step, we study the radial trends of the distribution of both the new and old H I holes in the galaxy as opposed to the SFR. Figure 10 compares the radial surface density of H I holes
<figure><img src="image_8.png" /></figure>

**Fig. 7.** Kinetic age vs axial ratio of the new (blue) and old (grey) holes, with the size of the symbols being representative of the hole diameter. This denotes a very weak negative correlation between axial ratio and kinetic age in our sample. This statistically insignificant relationship confirms that one cannot correlate the shape of the holes to its age.

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

**Fig. 8.** Variation of age with galactocentric distance. The Pearson *r* value = 0.2 with a *p*-value of 0.08. This indicates a weak positive linear correlation between age and galactocentric distance. In other words, as galactocentric distance increases, age tends to increase slightly, but the relationship is not strong.

to the dust-corrected FUV SFR profile (Refer to Appendix A for FUV-derived SFR calculation). We found a moderate anti-correlation (*r* = -0.57, *p* = 0.017) between SFR surface density, ΣSFR(R) and surface density of the H I holes. Regions with the highest ongoing star formation (inner 3–5 kpc) surprisingly host fewer large H I cavities, perhaps because vigorous feedback has already disrupted old shells or because high SFR goes hand in hand with high gas density, which can either refill the holes more rapidly or stall hole expansion before large cavities develop. This apparent deficit could also reflect a detection bias, since feedback-driven holes near galaxy centres are more difficult to distinguish from the underlying lack of atomic gas in regions dominated by molecular gas. Beyond 5 kpc, the number of holes rises as the SFR surface density tails off, suggesting that past episodes of clustered SNe have carved long-lasting cavities in the outer disk.

We noted that several new H I holes rarely coincided with the bubble feature visible in SITELLE RGB imagery. This might be

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

**Fig. 9.** Upper limits of the expansion velocity (given by the average velocity dispersion inside the hole) vs their diameters. The coloured curves are traced assuming the holes are powered by SNe, with total energy labelled on each curve. The dashed grey lines correspond to different ages (see text).

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

**Fig. 10.** Radial distribution of H I-hole counts compared to FUV-derived, dust-corrected star-formation-rate surface density. The two profiles exhibit a statistically significant moderate anti-correlation (Pearson *r* = -0.57, *p* = 0.017), indicating that regions of elevated recent star formation tend to host fewer H I holes.

due to the different time scales involved, and also to the previous star formation site and consequent SNe having triggered new star formation in the adjacent region. Another simpler explanation would be the fact that the part of the disk that pushed away the neutral gas (thereby creating cavities) has now moved due to the rotational motion.

## 5.2. Fountains in face-on galaxies

The combination of shell morphology, stellar associations, and kinematic evidence builds a compelling case that stellar feedback is responsible for carving these cavities in M101's neutral gaseous disk. Comparing this to the edge-on galaxies where one usually detects extraplanar gas or line-of-sight lagging halos as evidence of fountain flows (e.g. Rand 1998), here we see the imprints of those flows on the disk and the launch sites of the
<table><caption>**Table 5.** Properties of the nine galactic fountains.</caption><thead><tr><td>Hole</td><td>Max. Height [pc]</td><td>Time of flight [Myr]</td></tr></thead><tbody><tr><td>1</td><td>985.03</td><td>36.48</td></tr><tr><td>2</td><td>1767.22</td><td>86.21</td></tr><tr><td>3</td><td>134.19</td><td>20.64</td></tr><tr><td>4</td><td>1368.89</td><td>97.78</td></tr><tr><td>5</td><td>93.14</td><td>20.70</td></tr><tr><td>6</td><td>219.22</td><td>43.84</td></tr><tr><td>7</td><td>253.94</td><td>63.49</td></tr><tr><td>8</td><td>322.58</td><td>92.16</td></tr><tr><td>9</td><td>444.55</td><td>98.79</td></tr></tbody></table>

Notes. Col. 2: Maximum vertical height $Z_{max}$ from the $Z_0$ Eq. (5) by replacing $v_{disp}$ with $v_{exp}$. Col. 3: Total time of fountain flight as $2Z_{max}/v_{exp}$.

fountains. This proves that studying face-on galaxies offers a complementary perspective.

All nine high-quality fountain candidates coincide with regions of recent star formation. In each case, we detected $H\alpha$ emission within the cavity or along its rim. The former case indicates that the hole is formed in the last few Myrs and the responsible stars are still alive. In the latter case, where $H\alpha$ is rim-dominant, they are almost always coincident with FUV emissions in the centre. From this, we infer that the hole was created hundreds of Myr ago: the ionising population has aged, yet secondary ‘contagious’ star formation along the compressed shell continues to power localised $H\alpha$ along the rims. In every fountain site, the previous generation of stars has successfully displaced the neutral gas, revealing clear depressions in H I moment-zero map, accompanied by elevated velocity dispersion in both H I and $H\alpha$ (moment-2).

Re-examining Figure 4, the four largest cavities exhibit classical super-shell or super-bubble morphologies. Two of these (leftmost and rightmost panels) were previously classified as a super-bubble and a super-shell by Kamphuis et al. (1991) and Chakraborti & Ray (2011), respectively. Their implied energy requirements ~$10^{56-57}$ erg exceed by two orders of magnitude the budget of our smaller fountain sites.

Matonick & Fesen (1997) carried out an optical survey of M101 and catalogued 93 SNRs, classifying their morphologies as stellar, filled, diffuse, arc or shell. When we overlaid their positions on the [S II]/$H\alpha$ map (bottom panel of Figs. 4 and A.1), five of our nine high-quality fountain candidates lie within a few hundred parsecs of these SNRs. In particular, Holes 1 and 3 are delineated by multiple remnants, and additionally, several of our intermediate Q-score detections (purple ellipses) also coincide with SNRs along their rims.

The multi-wavelength dataset of these nine sites does not show any strong evidence of inflowing gas. However, we can deduce that from other observations of edge-on galaxies, where we see gas above the disk, and also from the amplitude of the expansion velocity, which is far below the escape velocity estimated to be $240 \text{ km s}^{-1}$ (Kuntz et al. 2003), the gas expelled by the SNe (at median velocities of $34.5 \text{ km s}^{-1}$) must return to the disk. To more precise the maximum height at which the H I gas can be expelled, and the time-scale of its travel outside the plane, we have computed these quantities, from $v_{exp}$ and the restoring force from the stellar disk, in Table 5. This gives an order of the magnitude of the feedback efficiency when the gas is unavailable for star formation.

### 5.3. Mechanical energy

All of the mechanical energy imparted to the ISM comes from OB stars via ionising photons, stellar winds or SNe explosions (Abbott 1982; McCray & Kafatos 1987; Agertz et al. 2013; Kim & Ostriker 2015).

Of the two energy estimates in this study, we have shown in Table 4 that both equations from Chevalier (1974) and McCray & Kafatos (1987) give comparable results, despite their different sensitivity to expansion velocities and diameters. It is possible to use these estimates to calculate the number of SNe that produced these features. According to the formula by McCray & Kafatos (1987), the radius of an H I shell is given by:

$$R_s[\text{pc}] = 97 \left( \frac{N_* E_{51}}{n_0} \right)^{1/5} t_7^{3/5}, \quad (13)$$

$$N_* = \frac{R_*^5 n_0}{97^5 t_7^3 E_{51}}, \quad (14)$$

where $N_*$ is the number of SNe, $t_7$ is the age of the shell in $10^7$ years, $E_{51}$ is the energy deposited by the SNe in units of $10^{51}$ erg and $n_0$ is the density of the ambient interstellar medium in $\text{cm}^{-3}$ which can be considered to be equal to $n_H$ as we did in Section 3.3.2. We found $N_*$ between 100 and $10^4$, for the nine holes considered here in detail.

## 6. Conclusions

We have investigated galactic fountains in the almost face-on, spiral galaxy M101, comparing ionised, neutral gas and star clusters mapped through $H\alpha$, H I-21 cm and FUV. We use high-resolution data from THINGS for the H I emission, NGS-GALEX for UV, and SITELLE-SIGNALS IFU for the $H\alpha$ tracer. The $H\alpha$ gas is ionised by OB stars of a lifetime shorter than 10 Myr, and the FUV emission remains after 100 Myr, dating the last star-forming event. Supernovae explosions accompanying the star formation blow away the neutral gas, digging H I-holes in the galaxy plane, and sending material above the plane, which then infalls back to the plane like a fountain (Fraternali & Binney 2006, 2008). This delays star formation and plays a feedback role. We revisited the 52 H I holes catalogued by K93, and found 20 more holes, essentially in the outer parts of the galaxy. We have characterised the feedback processes and quantified their efficiency through measurements of sizes, ages, mechanical energies, and displaced gas masses. These 20 new holes are smaller in size and expansion velocity, mainly due to the increased spatial and velocity resolution. There is a deficit of small holes in the inner regions. This could be a selection bias, since small holes are easier to detect in the outer parts, where the gas density is lower, while in the inner parts, the small holes are more rapidly refilled or their expansion gets stalled. All H I hole properties, such as their size, axial ratio, their atomic gas column density and their kinetic energy, based on the expansion velocity $v_{exp}$ of the bubble, are studied as a function of their distance to the galaxy centre. There is no correlation between the axis ratio and the distance from the centre, maybe because this late-type galaxy has almost solid-body rotation until 10 kpc, meaning little shearing of the bubbles. Their age is estimated from their size and $v_{exp}$, assumed constant during the bubble expansion. The energy of the bubble is computed from $v_{exp}$, its diameter, and the volumic density of the gas $n_0$ in which it expands, from two methods (Chevalier 1974; McCray & Kafatos 1987). We derive from these energies that the holes were created by several SNe between 100 and $10^4$. We mapped in more
detail the nine holes satisfying strong criteria to be a true foun-
tain effect. Most of these holes contain FUV emission inside, and
some Hα at their borders, or somewhat offset, indicating conta-
gious star formation. Young stars and Hα are found inside only
one small hole, showing that these are relatively rare. The largest
hole of 2.4 kpc and the oldest age (94 Myr) has no Hα nor UV
emission, and is consistent with an old fountain effect, when the
star formation has terminated. Globally, the holes are relatively
old, and it is difficult to find in their centre the very star formation
at the origin of the hole. Therefore, no correlation can be found
between the present star formation and the mechanical energy
of the hole. The kinetic age of H I holes is remarkably constant
with radius, around 50 Myr, and does not vary with size and $v_{exp}$.
The number of holes is smaller in the galaxy centre, where star
formation is higher, which may look like a paradox, but is in fact
a selection bias. It is more difficult to see H I holes where the gas
is denser, and the rotation period is shorter, since holes are then
smaller, and destroyed by star formation feedback. Finally, we
note that some regions might also be perturbed by tidal interac-
tions, as is likely for Hole 2, since M101 is highly perturbed, as
is obvious from its remarkable lopsidedness.

Acknowledgements. AS would like to thank Laurent Drissen and Thomas Martin for their valuable help with analysing the SITELLE data. This work was supported by the Action Thématique Cosmology-Galaxies (ATCG) of the CNRS/INSU PN Astro. This research was based on observations obtained at the CFHT, which is operated from the summit of Mauna Kea by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the Uni- versity of Hawaii. The authors wish to recognise and acknowledge the very significant cultural role that the summit of Mauna Kea has always had within the indigenous Hawaiian community. The authors are most grateful to have the opportunity to conduct observations from this mountain. The observations were obtained with SITELLE, a joint project between Université Laval, ABB- Bomem, Université de Montréal, and the CFHT, with funding support from the Canada Foundation for Innovation (CFI), the National Sciences and Engineer- ing Research Council of Canada (NSERC), Fonds de Recherche du Québec – Nature et Technologies (FRQNT), and CFHT. The collaboration is grateful to the FRQNT, CFHT, the Canada Research Chair program, the National Science foundation NSF – 2109124, Natural Sciences and Engineering Research Coun- cil of Canada NSERC – RGPIN-2023-03487, the Swedish Research Council, the Swedish National Space Board, the Royal Society, and the Newton Fund via the award of a Royal Society-Newton Advanced Fellowship (NAF\R1\180403), FAPESC, CNPq, FAPESP (2014/11156-4), FAPESB (7916/2015), and CONA- CyT (CB2015-254132).

References

Abbott, D. C. 1982, ApJ, 263, 723

Agertz, O., Kravtsov, A. V., Leitner, S. N., & Gnedin, N. Y. 2013, ApJ, 770, 25

Armillotta, L., Fraternali, F., & Marinacci, F. 2016, MNRAS, 462, 4157

Bagetakos, I., Brinks, E., Walter, F., et al. 2011, AJ, 141, 23

Boomsma, R., Oosterloo, T. A., Fraternali, F., van der Hulst, J. M., & Sancisi, R. 2008, A&A, 490, 555

Bosma, A., Goss, W. M., & Allen, R. J. 1981, A&A, 93, 106

Bottema, R. 1993, A&A, 275, 16

Bregman, J. N. 1980, ApJ, 236, 577

Brinks, E., & Bajaja, E. 1986, A&A, 169, 14

Camps-Fariña, A., Zaragoza-Cardiel, J., Beckman, J. E., et al. 2017, MNRAS, 468, 4134

Casasola, V., Cassarà, L. P., Bianchi, S., et al. 2017, A&A, 605, A18

Chakraborti, S., & Ray, A. 2011, ApJ, 728, 24

Chevalier, R. A. 1974, ApJ, 188, 501

Combes, F. 1991, A&A, 243, 109

Combes, F., & Becquaert, J. F. 1997, A&A, 326, 554

Comrie, A., Wang, K. S., Hsu, S. C., et al. 2021, Astrophysics Source Code Library [record ascl:2103.031]

Davies, R. D. 1975, MNRAS, 170, 45

Deul, E. R., & den Hartog, R. H. 1990, A&A, 229, 362

Drissen, L., Martin, T., Rousseau-Nepton, L., et al. 2019, MNRAS, 485, 3930

Egorov, O. V., Kreckel, K., Glover, S. C. O., et al. 2023, A&A, 678, A153

Fraternali, F., & Binney, J. J. 2006, MNRAS, 366, 449

Fraternali, F., & Binney, J. J. 2008, MNRAS, 386, 935

Fraternali, F., van Moorsel, G., Sancisi, R., & Oosterloo, T. 2002, AJ, 123, 3124

Ginsburg, A., Robitaille, T., & Beaumont, C. 2016, Astrophysics Source Code Library [record ascl:1608.010]

Hao, C.-N., Kennicutt, R. C., Johnson, B. D., et al. 2011, ApJ, 741, 124

Heiles, C. 1979, ApJ, 229, 533

Heiles, C. 1984, ApJS, 55, 585

Houck, J. C., & Bregman, J. N. 1990, ApJ, 352, 506

Hu, E. M. 1981, ApJ, 248, 119

Kamphuis, J. J. 1993, Ph.D. Thesis, University of Groningen, Netherlands

Kamphuis, J., & Briggs, F. 1992, A&A, 253, 335

Kamphuis, J., Sancisi, R., & van der Hulst, T. 1991, A&A, 244, L29

Keel, W. C., Holberg, J. B., & Treuthardt, P. M. 2004, AJ, 128, 211

Kim, C.-G., & Ostriker, E. C. 2015, ApJ, 802, 99

Kregel, M., van der Kruit, P. C., & de Grijs, R. 2002, MNRAS, 334, 646

Kuntz, K. D., Snowden, S. L., Pence, W. D., & Mukai, K. 2003, ApJ, 588, 264

Lara-López, M. A., Pilyugin, L. S., Zaragoza-Cardiel, J., et al. 2023, A&A, 669, A25

Leroy, A. K., Walter, F., Brinks, E., et al. 2008, AJ, 136, 2782

Leroy, A. K., Walter, F., Sandstrom, K., et al. 2013, AJ, 146, 19

Martin, D. C., Fanson, J., Schiminovich, D., et al. 2005, ApJ, 619, L1

Martin, T., Drissen, L., & Joncas, G. 2015, ASP Conf. Ser., 495, 327

Martin, T. B., Prunet, S., & Drissen, L. 2016, MNRAS, 463, 4223

Matonick, D. M., & Fesen, R. A. 1997, ApJS, 112, 49

McCray, R., & Kafatos, M. 1987, ApJ, 317, 190

Morisset, C. 2018, https://doi.org/10.5281/zenodo.1206115

Oosterloo, T., Fraternali, F., & Sancisi, R. 2007, AJ, 134, 1019

Pokhrel, N. R., Simpson, C. E., & Bagetakos, I. 2020, AJ, 160, 66

Puche, D., Westpfahl, D., Brinks, E., & Roy, J.-R. 1992, AJ, 103, 1841

Putman, M. E., Peek, J. E. G., & Joung, M. R. 2012, ARA&A, 50, 491

Rampazzo, R., Mazzei, P., Marino, A., et al. 2022, A&A, 664, A192

Rand, R. J. 1998, PASA, 15, 106

Rousseau-Nepton, L., Martin, R. P., Robert, C., et al. 2019, MNRAS, 489, 5530

Rupke, D. S. N. 2018, Galaxies, 6, 138

Schruba, A., Leroy, A. K., Walter, F., et al. 2011, AJ, 142, 37

Shapiro, P. R., & Field, G. B. 1976, ApJ, 205, 762

Shappee, B. J., & Stanek, K. Z. 2011, ApJ, 733, 124

Somerville, R. S., & Davé, R. 2015, ARA&A, 53, 51

Tody, D. 1986, SPIE Conf. Ser., 627, 733

van der Hulst, T., & Sancisi, R. 1988, AJ, 95, 1354

van der Kruit, P. C., & Searle, L. 1982, A&A, 110, 61

van Woerden, H., Peletier, R. F., Schwarz, U. J., Wakker, B. P., & Kalberla, P. M. W. 1999, ASP Conf. Ser., 166, 1

Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, AJ, 136, 2563

Watkins, E. J., Kreckel, K., Groves, B., et al. 2023a, A&A, 676, A67

Watkins, E. J., Barnes, A. T., Henny, K., et al. 2023b, ApJ, 944, L24

Westmeier, T., Braun, R., Brüns, C., Kerp, J., & Thilker, D. A. 2007, New Astron. Rev., 51, 108

Wilcots, E. M., & Miller, B. W. 1998, AJ, 116, 2363
## Appendix A: SFR calculation using GALEX data

We use the intensity (fd(nd)-int), sky-background (fd(nd)-skybg) and weight (fd(nd)-wt) from the GALEX data products to calculate the radial distribution of SFR of M101. The intensity and the sky-background maps are in units of counts per pixel, while the weight maps are dimensionless. The higher the weight of a pixel, the more confident we are in its intensity value.

The flux measurements have been made from the background-subtracted versions of the intensity images, whose units are in counts per second (CPS) by multiplying with a factor of 1.40 × 10¹⁵ for FUV. This gives the flux, $F_ν$, in erg s⁻¹ cm⁻² Å⁻¹. The FUV luminosity is obtained by:

$$L_ν[\text{erg s}^{-1}] = 4πD_L^2 F_ν, \qquad (A.1)$$

where $D_L$ is the luminosity distance to galaxy M101 in cm, $ν = 442$ is the FUV bandwidth. This luminosity needs to be corrected for dust attenuation, for which we must first calculate the FUV and NUV magnitudes using the background-subtracted intensity in counts per second, $I$:

$$m_{\text{FUV}} = -2.5 \log_{10}(I_{\text{FUV}}) + 18.82, \qquad (A.2)$$

$$m_{\text{NUV}} = -2.5 \log_{10}(I_{\text{NUV}}) + 20.08. \qquad (A.3)$$

The attenuation in FUV is thus obtained by using the following formula from Hao et al. (2011):

$$A_{\text{FUV}} = 3.83 [(m_{\text{FUV}} - m_{\text{NUV}}) - 0.022]. \qquad (A.4)$$

The luminosity is then corrected using:

$$L_{\text{FUV}}^{\text{corr}} = 10^{0.4A_{\text{FUV}}} L_{\text{FUV}}. \qquad (A.5)$$

To obtain the dust-corrected SFR from the FUV luminosity, we calculate:

$$\text{SFR}_{\text{FUV}}^{\text{corr}}[\text{M}_\odot \text{yr}^{-1}] = \frac{L_{\text{FUV}}^{\text{corr}}}{10^{43.35}}. \qquad (A.6)$$
<figure><img src="image_12.png" /><figcaption>Fig. A.1. Same as Figure 4 for the smaller five sites of galactic fountains in M101.</figcaption></figure>