## Letter to the Editor

## North Polar Spur : gaseous plume(s) from star-forming regions at ∼ 3-5 kpc from Galactic Center?

Eugene Churazov 1 2 , Ildar I. Khabibullin , 3 , 1 , 2 , Andrei M. Bykov 4 , Nikolai N. Chugai , Rashid A. Sunyaev 5 2 , 1 , Victor P. Utrobin 5 , 6 , and Igor I. Zinchenko 7

- 1 Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany
- 2 Space Research Institute (IKI), Profsoyuznaya 84/32, Moscow 117997, Russia
- 3 Universitäts-Sternwarte, Fakultät für Physik, Ludwig-Maximilians-Universität München, Scheinerstr.1, 81679 München, Germany
- 4 Ioffe Institute, Politekhnicheskaya st. 26, Saint Petersburg 194021, Russia
- 5 Institute of Astronomy, Russian Academy of Sciences, 48 Pyatnitskaya str., Moscow 119017, Russia
- 6 NRC 'Kurchatov Institute', acad. Kurchatov Square 1, Moscow 123182, Russia
- 7 Institute of Applied Physics of the Russian Academy of Sciences, 46 Ul'yanov str., Nizhny Novgorod 603950, Russia

November 4, 2024

## ABSTRACT

We argue that the North Polar Spur (NPS) and many less prominent structures are formed by gaseous metal-rich plumes associated with star-forming regions (SFRs). The SFRs located at the tangent to the 3-5 kpc rings might be particularly relevant to NPS. A multitemperature mixture of gaseous components and cosmic rays rises above the Galactic disk under the action of their initial momentum and buoyancy. Eventually, the plume velocity becomes equal to that of the ambient gas, which rotates with different angular speed than the stars in the disk. As a result, the plumes acquire characteristic bent shapes. An ad hoc model of plumes' trajectories shows an interesting resemblance to the morphology of structures seen in the radio continuum and X-rays.

Key words. ISM: jets and outflows - X-rays: binaries - plasmas - acceleration of particles

## 1. Introduction

North Polar Spur (NPS, Tunmer 1958) is an elongated bright structure first discovered in all-sky radio images along with other similar structures dubbed Loops (e.g. Dickinson 2018, for a recent review). Unlike other Loops, NPS (the northern part of Loop I) is visible in soft X-rays across its full extent (Bunner et al. 1972; Egger &amp; Aschenbach 1995; Kataoka et al. 2018; Predehl et al. 2020; LaRocca et al. 2020; Kataoka et al. 2021).

reached from, e.g., the Faraday Rotation and X-ray absorption measurements (see a comprehensive discussion in Lallement 2023). Recent discoveries of the giant structures in gamma-rays (Fermi bubbles, Su et al. 2010) and soft X-rays in the southern Galactic sky (eROSITA, Predehl et al. 2020), the "bubbles" appear to be more symmetric with respect to the Galactic center, strengthening the case for the central energy release as a viable solution (e.g. Sarkar 2024).

Given the NPS's large angular size, ≳ 90 deg, the first proposed explanation (e.g. Hanbury Brown et al. 1960; Berkhuĳsen 1971) was a shell of a very nearby, 𝑑 ≲ 100 pc, supernova remnant. Alternatively, a scenario in which NPS represents a distant, 𝑑 ≳ 10 kpc, Galaxy-scale structure has been put forward, with the possible energetic outburst in the Galactic center being considered as the primary source behind it (Sofue 1977, 1994). The latter might be caused by an intense star-formation episode in the vicinity of the Galactic Center or by the activity of the supermassive black hole Sgr A* (e.g. Sofue 1977, 1994; Sarkar et al. 2015; Sofue &amp; Kataoka 2021; Yang et al. 2022). Both scenarios suggest that NPS is powered by a shock wave propagating in the hot medium (either interstellar or circumgalactic). This shock is responsible for the acceleration of radio-emitting relativistic electrons, as well as the heating and compression of the X-ray-emitting hot gas.

Although these two scenarios involve very different physical sizes, time scales and required energetics, it turned out to be very difficult to unequivocally distinguish between them based on the available data, with contradicting conclusions having been

In this Letter, we consider another scenario for NPS formation motivated by the morphological and spectral properties of the soft X-ray emission measured in the course of the eROSITA all-sky survey (Predehl et al. 2021; Sunyaev et al. 2021). We propose that NPS is produced by a break-out of the massive star formation regions associated with the end of the Galactic bar, whichlocates its base at 𝑑 ∼ 5kpcfromus(e.g., Bland-Hawthorn &amp; Gerhard 2016). Qualitatively similar arguments have recently been discussed in the work of Shimoda &amp; Asano (2024). Flows of enriched hot gas move upward from the Galactic disk and get entrained in the relative rotation of the circumgalactic medium above the disk, resulting in the spiral-like structure, appearing like a loop in the sky. In this model, it is the advection of the enriched gas and relativistic particles, rather than a shock, that is responsible for the appearance of NPS, predicting the high metal abundance of X-ray emitting gas and lack of evolutionary signatures in the direction perpendicular to the NPS edge.

Fig. 1. X-ray and radio images of the NPS region. The left and middle panels: the 0.52-0.61 keV and 0.7-0.96 keV eROSITA X-ray images, respectively. For comparison, the right panel shows the 408 MHz radio image from the all-sky continuum surveys (Haslam et al. 1981; Remazeilles et al. 2015). The images in galactic coordinates are shown in stereographic projection for a better view of the regions near the Galactic plane and the Galactic North Pole. The X-ray images are particle background subtracted, and exposure corrected. The brightest compact sources and galaxy clusters have been masked and the resulting image convolved with 𝜎 = 20 ′ Gaussian. Still prominently visible in the image are the RS Oph (bright dot near the bottom of the left image) and the stray light halo around Sco X-1 (at the right edge of the middle image). The white dashed lines show the wedge used for the extraction of radial profiles. The numbers indicate the distance (in degrees) from the wedge center, which is at ( 𝑙, 𝑏 ) = ( 338 ◦ , 32 ◦ ) in Galactic coordinates.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-ad541458-9eff-4d1c-911c-eea1d1294585.jpg)

## 2. X-ray picture

For X-ray analysis, we use the data accumulated by eROSITA (Predehl et al. 2021) telescope onboard the SRG X-ray observatory (Sunyaev et al. 2021) in the course of its all-sky survey in the Eastern Galactic hemisphere (i.e. the Galactic longitude range 0 ◦ &lt; 𝑙 &lt; 180 ). Details on the data preparation and analy-◦ sis are given in Appendix A. Images extracted from the eROSITA all-sky survey in the 0.52-0.61 keV and 0.7-0.96 keV bands are shown in Fig. 1 (left and middle panels, respectively). The former band includes the triplet of He-like oxygen O VII near 574 eV, while the bright lines of helium-like neon Ne IX (at 905 and 921 eV) and neon-like iron Fe XVII (between 725 and 826 eV) fall into the latter. For plasma in collisional ionization equilibrium (CIE), the peak line emissivities in these bands are at temperatures ∼ 0 2 keV and . ∼ 0 6 keV, respectively. .

brightened shell. Similarly, the 0.7-0.96 keV surface brightness is composed of a fainter diffuse emission surrounding a bright central "filament".

Rather complex X-ray morphology of NPS is obvious from these images and the comparison between them. One can readily see that the correlation length of X-ray structures is larger in the direction along NPS for both bands but a simple combination of concentric spherical shells is not a good description of NPS. For example, the surface brightness in the 0.52-0.61 keV band dominated by the O VII line, has a sharp outer edge but a rather flat surface brightness across the inner region, unlike an edge-

More quantitatively this is illustrated via profiles of the X-ray emission in several narrow bands. The white dashed lines in Fig. 1 outline a wedge that is aligned with the NPS "outer" boundary. Radial profiles extracted from this wedge are shown in Fig. 2. The steep rise, corresponding to the "edge" of NPS, is visible in all bands. However, further inside NPS, the profiles do not follow a typical edge-brightened behavior characteristic for a spherical shock associated with an instantaneous point-like energy release in a uniform medium (see, e.g., profiles at energies below 0.61 keV). The observed profiles could be affected by the upstream gas density gradients, Non-Equlibrium Ionization (NEI) effects downstream of the shock, and, especially, by the energy release mode. For example, in a 'steady wind' scenario (see, e.g. Sarkar et al. 2015), the surface brightness profile is rather flat.

In the model without the shock that we pursue here, this behavior suggests that the volume occupied by the plasma emitting at these energies resembles a flattened sheet rather than a spherical shell. In harder bands, the profiles feature a steep jump, followed by a gradual increase of intensity before reaching a peak at 𝑥 ≈ 35 degrees in this plot. The center of the wedge is at ( 𝑙, 𝑏 ) = ( 338 ◦ , 32 ◦ ) in Galactic coordinates.

An alternative to the local SNR or GC-related scenarios, which associate the steep rise with a shock, is the assumption that NPS is a multi-temperature gas, where a hotter ( &gt; 0.3 keV) component, responsible for the emission in iron lines above 0.7 keV, coexists with a warm (e.g., 0.15-0.2 keV) plasma, producing the bulk of the emission below 0.7 keV and most prominently in oxygen lines.

The sharpness of the NPS edge is considered as a strong argument in favor of the shock scenario (e.g. Predehl et al. 2020). However, there are examples of astrophysical objects when very sharp edges are seen in the hot X-ray-emitting plasma. The most relevant are 'cold fronts' in galaxy clusters (see Markevitch &amp; Vikhlinin 2007; Zuhone &amp; Roediger 2016, for reviews). These cold fronts are contact discontinuities rather than shocks, which might retain their sharpness due to magnetic fields and/or accelerated gas motions along the front (see, e.g. Churazov &amp; Inogamov 2004). Similar effects might be relevant to the rising plumes discussed here.

Accurate temperature and abundance measurements for gas with temperature below 1 keV (when pure bremsstrahlung continuum is very weak and subdominant to other continuum components) are highly model-dependent. As is clear from Fig.2, the ratio of fluxes in different bands varies strongly across NPS. This means that the best-fitting temperature strongly depends on the region used for spectrum extraction. It further implies that a mixture of components with different temperatures (and different ionization states, especially in the shock-based scenario) will inevitably be present. The abundance determination for the gas in the relevant temperature range (0.1-0.6 keV) is highly degenerate with the emission measure even when the gas is in the Collisional Ionization Equilibrium (CIE). In practice, it is difficult to measure the abundance if it is larger than 0.1 Solar. The effects of nonequilibrium ionization that are expected in the shock scenario, further affect the 'apparent' temperature and abundance. These issues are likely the reason for different temperature/abundance measurements available in the literature (e.g. Kataoka et al. 2013; Yamamoto et al. 2022) for various patches of the NPS. The exceptionally uniform and sensitive eROSITA all-sky data give a clear view of the amplitude of spectral variations across such large regions for the first time. We defer the extended spectral analysis of the NPS to future work.

In what follows, we consider a scenario in which such multitemperature gas is venting from the active star-forming regions in the central parts of the Galactic disk.

## 3. Morphological model

Here we consider a simple morphological model that is based on a broad-brush picture of the evolution of a massive star-forming region in the Galactic disk (e.g. Norman &amp; Ikeuchi 1989; Kim et al. 2020). The first shock waves of supernova explosions originating from the most massive newborn stars in such a region produce hot and dense gas strongly enriched with metals of the ejecta. Initially confined by the surrounding dense and cold gas, the hot gas manages to break through and form an (adiabatically cooled) plume that moves away from the disk and interacts with the interstellar (ISM) and circumgalactic (CGM) media forming some sort of chimneys and plumes (see Appendix D). The trajectories of rising plumes might be complicated, especially when the gas is multi-phase and magnetized. They are also affected by the pressure, magnetic fields, and velocities in the ambient ISM and/or CGM (e.g. Faucher-Giguère &amp; Oh 2023).

For illustration purposes, we consider an ad hoc model that encapsulates all the complicated physics in three parameters.

Namely, the initial (vertical) velocity of the plume 𝑣 𝑧, 0 and two effective scale heights 𝑧 𝜙 and 𝑧 𝑧 . Here 𝑧 𝜙 describes how the plume that initially moves together with stars and gas close to the disk plane, decelerates and joins the pressure-supported hot halo gas, which is either non-rotating or rotating slower than the stellar disk. In turn, the 𝑧 𝑧 parameter controls the decline of the vertical (perpendicular to the disk) velocity component. As a result, the two velocity components have the following dependence on 𝑧 :

$$\text{$\text{$manvxy}$} \\ \text{$\text{the most} \quad v_{z}(z)=v_{z,0} e^{-z/z_{z}$ and $v_{\phi}(z)=v_{rot}$} \left [ 1 - e^{-z/z_{\phi}$} \right ]. \quad \quad ( 1 )$$

A flat rotation curve with 𝑣 rot = 220 km s -1 is adopted for the sake of simplicity, the effective scale heights are set to 𝑧 𝜙 = 2 kpc and 𝑧 𝑧 = 5 kpc, and 𝑣 𝑧, 0 = 200 km s -1 . The choice of these parameters is rather arbitrary, but it suffices for illustration purposes and the final answer is not dramatically sensitive to it.

For the source of the plume, e.g. a star-forming region, we assume a circular trajectory in the disk plane with the velocity 𝑣 𝜙, src . In our illustrative model, the following two versions are considered:

$$\underset { \text{and} \, \equivfty } { \text{ss NPS.} } \quad ( a ) \ v _ { \phi, \text{src} } = v _ { \text{rot} } \text{ and } ( b ) \ v _ { \phi, \text{src} } = \Omega R,$$

where 𝑅 is the distance from the Galactic Center. Version (a) corresponds to the source moving with the same velocity as stars, while version (b) mimics a pattern motion. The former case is relevant when the active phase of a given star-forming region associated with the dense gas is very long and the active region always moves together with the gas (with the velocity 𝑣 rot ). The latter case corresponds to a situation when the spatial distribution of the plume sources is linked to the Galactic bar and/or spiral arms. For illustration of this case, we adopted the value Ω = 33 km s -1 kpc -1 derived for the bar pattern speed (e.g. Clarke &amp; Gerhard 2022). As a further simplification, we assume that none of these parameters depend on the position of the plume within the galaxy.

The above relations (1-2) can be trivially integrated to derive a position of the gas lump released some time 𝑡 ago relative to the current position of the "source". By design, plume trajectories appear in 3D as spirals rising from the disk from the current location of the plume source (see Appendix B and Fig. B.1). We note here that for a strongly underdense lump, the buoyancy force is solely set by pressure gradients in the ambient gas. These gradients themselves depend on the rotation velocity of the halo gas and can have both radial and vertical components. In the above model, we explicitly neglect the former component, although in real systems it can force the plume to move radially. In particular, for slow (fast) rotation of the halo gas, underdense lumps can move to larger (smaller) radii (see Fig. B.2).

As a final step, we project the trajectory onto the sky as seen from the position of the Sun. This modifies the appearance of the plume due to different distances from the Sun to various segments of the plume and the Sun's motion relative to the source. Anexample of plume trajectories for a set of "sources" chosen by hand is shown in Fig 3. A table with the positions of sources is given in the Appendix B. The trajectories have been integrated for ∼ 140 Myr, which is approximately one rotation period around the Galactic Center near the Sun for the adopted value 𝑣 rot .

Figure 3 illustrates the key morphological properties of the plumes. Naturally, the trajectories remain confined to the regions defined by the radial distance to the source from the Galactic Center, while the notable left-right asymmetry is imposed by the rotation direction of the Galaxy. All trajectories have rising parts

Fig. 2. Left: Radial profiles of the NPS X-ray emission in several energy bands. The wedge used for the flux extraction is shown in Fig. 1. The bright edge of the NPS is at 𝑥 ≈ 48 degrees. The contributions of the instrumental background and distant (unresolved) sources (CXB) have been subtracted, leaving only the flux from the Galaxy. Right: The same profiles as in the left plot, but after subtracting the mean level at 𝑥 = 60 -70 degrees, renormalizing by the flux at 𝑥 = 35 degrees, and lightly smoothing with a Gaussian filter ( 𝜎 = 0 4 . ◦ ). In addition, the radio surface brightness profile (Haslam, 408 MHz) is shown with the red dashed curve after applying the same procedure (subtraction, renormalization, and smoothing) as used for the X-ray profiles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-8665b9ca-0a11-4928-89d7-2bebc346c7d9.jpg)

close to the Galactic Plane, with the direction of the curvature set by the mutual position of the source and the Sun and the direction of the Galaxy rotation. This is clearly illustrated by the family of trajectories for sources located at the 5 kpc ring around GC. The morphological difference is even more pronounced for sources located at larger distances from GC (cases marked at Cygnus X and Vela in Table B.1). For the Cygnus X region, the plume's apparent trajectory rises to the Galactic Poles before forming a spiral. For the Vela region, the plume is bent much earlier and follows an almost horizontal line at | 𝑏 | ∼ 20 -25 degrees.

If NPS is indeed a gaseous plume, its base should be associated with the star formation regions that are currently at Galactic longitudes 𝑙 ∼ 20 ◦ -40 ◦ (see Fig. 3). Incidentally, this range of longitudes corresponds to the tangential direction to the 3-5 kpc rings around the Galactic Center and indeed hosts the most active star-forming regions in the Galaxy, including a prototypical mini-starburst complex W43 (see Appendix E). Below we discuss this scenario in more detail.

## 4. Discussion

Motivated by the morphological model presented above, we assume that the leading edge of NPS is located at a tangent to the 5 kpc ring. In this case, the distance of the NPS base from the Sun is 𝐷 Sun NPS , ≈ 6 6 kpc. Of course, some parts of the plumes . can be closer to the Sun, while some further away.

Given the apparent size of NPS on the sky, one can assume that its physical size (the "depth" along the line of sight) is of the same order as the distance from the Sun, i.e. 𝑙 ∼ 𝐷 Sun NPS. Based , on that, we can estimate the hot gas density from the observed peak surface brightness 𝐼 𝑋 in the 0.7-1.05 keV band. This yields a proton number density

$$\text{mem} \quad \text{mem} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad  \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{)} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \vec { \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(-} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad 0 \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \hat \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{. } \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{-} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad ) \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \end{} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \\ \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{/ \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \tag \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \tt \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{ (} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad. \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{) \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \coprod \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \text{(} \quad \$$

where 𝑍 is the metal abundance with respect to the Solar value (the abundances relative to hydrogen from Asplund et al. 2009, are used here). In this derivation, we assumed that the gas is in collisional ionization equilibrium with temperature 0 6 keV . and used APEC model (Foster et al. 2012) to predict the emissivity. The metallicity dependence in the above expression is valid for 𝑍 / 𝑍 ⊙ ≳ 0 1. . The corresponding cooling time of the gas with temperature 𝑘𝑇 ≈ 0 6 keV can be estimated us-. ing the cooling function of Sutherland &amp; Dopita (1993) as GLYPH&lt;16&gt; 𝑛 GLYPH&lt;17&gt; -1 GLYPH&lt;16&gt; 𝑍 GLYPH&lt;17&gt; - / 1 2

$$\underset { \mathfrak { d } \nu n = } { \text{most} } \ t _ { \text{cool} } \sim 2 \times 1 0 ^ { 9 } \left ( \frac { n } { 1 0 ^ { - 3 } \text{ cm} ^ { - 3 } } \right ) ^ { - 1 } \left ( \frac { Z } { Z _ { \odot } } \right ) ^ { - 1 / 2 } _ { \cdots \cdots } \, y r.$$

While the 3D geometry of NPS is uncertain, we assume that its volume is ≈ 𝑙 3 (i.e. its line-of-sight extension is comparable to its transverse size) to get an estimate of the total energy (enthalpy) 𝐸 ∼ 𝑙 3 ( 𝑛 𝑝 + 𝑛 𝑒 ) 5 2 𝑘𝑇 ∼ 2 × 10 55 erg for 𝑙 = 6 kpc. Assuming that the integrated energy release associated with the formation of 1 solar mass of stars is ∼ 10 49 erg, it needs ∼ 2 × 10 6 𝑀 ⊙ of gas to be converted into stars. Assuming the star-formation rate of 0 1 . 𝑀 ⊙ yr -1 , it takes ∼ 20 Myr to generate enough energy for powering the entire NPS. The long cooling time and preferential accumulation of trajectories in the NPS region imply that it could be a result of cumulative contributions from many sources.

In the simplest version of the above scenario, the gas in NPS is in pressure equilibrium with the ambient medium. Assuming that the temperature of the ambient medium is ∼ 0 15 keV, the . regions filled with ∼ 0 6 keV plasma should have . ∼ 4 times lower density. Yet, NPS appears much brighter than the diffuse emission of the Galaxy. There are two plausible reasons for that. One is that for lines of Fe XVII, Ne IX, and Ne X, the gas in the halo is

Fig. 3. Top. Sample trajectories of the gaseous plumes superposed on the radio continuum map at 408 MHz map (Haslam et al. 1981; Remazeilles et al. 2015). These trajectories correspond to an arbitrarily chosen set of plume "sources" (see Table B.1). The red and blue colors correspond to cases when sources have the same azimuthal velocity in the disk as stars. They correspond to the plume sources currently observed at a given Galactic longitude but different distances from the Sun: the red and blue colors mark sources that are closer and further away, respectively. The white color corresponds to sources that move with the pattern speed. Three properties of the trajectories are worth mentioning: (i) the form "spirals" winded in the central region of the Galaxy well above the plane, (ii) East-West asymmetry set by the Galaxy rotation direction, and (iii) trajectories tend to overlap in certain regions of the sky even if they come from well-separated sources. Bottom. Polarized synchrotron emission map (e.g. Planck Collaboration et al. 2020) with a subset of trajectories superposed. The circles correspond to time tags (every 30 Myr) in this model. In the bottom panel, these tags are labeled.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-e9932ac8-62ef-4333-ab8f-92cf3ab0b777.jpg)

simply too cool to produce a bright emission. The second reason is that the plumes are plausibly much more metal-rich than the halo gas. This is particularly important for O VII lines, which are present in the spectrum of the halo. The gas in the plumes could easily be a factor of 10 or more metal-rich and even in the absence of shock, it can shine prominently (see Fig. 4). In Appendix C we argue that the plume metallicity is not much larger than the

Solar one. This implicitly suggests that the abundance in the halo is ≲ 0 1. .

An important implication of the model is the global asymmetry of the Milky Way diffuse X-ray and radio emission. The asymmetry is present even if the distribution of plumes' sources is itself symmetric. In 3D, the plumes can form a symmetric pattern, but the direction of the Galaxy rotation produces the appar-

Fig. 4. The spectrum of a large region inside NPS (red points) in comparison with the typical "Milky Way" spectrum well outside NPS (black points). Contributions of the detector background and CXB have been subtracted. The blue and green lines illustrate a few characteristic models. The solid blue line shows the APEC spectrum with the temperature 𝑇 w = 0 16 keV and abundance . 𝑍 / 𝑍 ⊙ = 0 05, using for the abundance . ratios from Asplund et al. (2009). The blue dashed line shows the same model with 3 times larger abundance, i.e. 𝑍 / 𝑍 ⊙ = 0 15. . The green solid line shows the APEC model with the temperature 𝑇 h = 0 5 keV, . 𝑍 / 𝑍 ⊙ = 0 05, and the emission measure multiplied by factor . ( 𝑇 w / 𝑇 h ) 2 . The dashed green line shows the same model with 𝑍 / 𝑍 ⊙ = 0 7. . A comparison of the histograms and the dashed lines shows that overabundances in the range of 3-10 are needed to reproduce enhanced brightness of the NPS compared to the Galaxy (assuming comparable pressures and linear sizes of emitting regions).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-5df5a453-9be4-4307-bf79-c8b279cbfe18.jpg)

ent asymmetry when viewed from the Sun's position. Some extra asymmetry could come from global gas motions in the Galactic halo (e.g. Mou et al. 2023), e.g. caused by interaction with the satellite galaxies, but it comes on top of the "natural" East-West asymmetry. The North-South asymmetry, which is clearly visible in both radio and X-ray sky, in this model is attributed to the asymmetry of the primary sources of the plume. When the hot gas finds its way through the dense cold gas in the disk, it may preferentially go to one side rather than forming a symmetric structure on both sides.

The morphological model considered above, treated trajectories of individual gas plumes as independent, neglecting the hydrodynamic nature of the flows and, in particular, possible gas mixing instabilities that would naturally arise in such a situation. The sharpness of the NPS's outer edge as well as its overall filamentary appearance might be a manifestation of the suppression of these instabilities by the well-ordered and sufficiently strong magnetic field.

Indeed, the synchrotron radio emission from the NPS region is known to be polarized (e.g. Sun et al. 2015, for a recent analysis). The polarization implies that the magnetic field is ordered. In the model discussed above, these structures can be especially prominent plumes located at different distances from the Galactic Center. The ordering is caused by stretching the field lines during certain phases of plume evolution and/or by the halo gas motions rather than by shocks (as illustrated in Appendix D).

On a more speculative side, we note that in this model the gas vented from the disk tends to accumulate high above the plane in the general direction of the Galactic Center. From this point of view, a question arises whether structures like the eROSITA bubbles (Predehl et al. 2020) could be produced by the same mechanism. In this case, these structures can evolve on a longer time scale than implied by the shock-driven scenario. In fact, different regions across NPS could come from gas lumps having different ages (and different star-forming regions) and, therefore, can have different properties.

In the proposed scenario, the NPS might be a Galactic analog of the magnetic structures directly observed to reach significant distances above the disks in many nearby galaxies, e.g. NGC 4217 (Stein et al. 2020). On the other hand, NPS might also be a case similar to the off-disk parts of the so-called anomalous spiral arms in NGC 4258, which are believed to be powered by the interaction of the relativistic jet with the galactic disk (Zeng et al. 2023, and references therein). In all these cases, we might witness signatures of the disk-halo interaction via intense feedback episodes resulting in complex multi-phase and magnetized interface region (e.g. a review by Beck 2015).

One observational test that can falsify the plume scenario is the measurements of the non-equilibrium ionization (NEI) signatures in the gas, which are pertinent to the shock scenario (e.g., Yamamoto et al. 2022). Indeed, for the gas electron density ∼ 10 -3 cm -3 and the downstream plasma velocity ∼ 200 km s -1 , the ionization parameter 𝜏 = 𝑛 𝑒 × ∼ 𝑡 1 5 . × 10 11 cm -3 s is expected at a distance of 1 kpc downstream of the shock. Such scales are easily resolved for any distance to NPS (for the "nearby" shock scenario the ionization parameter is even lower at the same angular distance from the shock) and can be derived from the spectra. The complication here is that NPS does not have a simple "layered" structure that allows for clean separation of shells at different 3D distances from the edge. As a result, the lines of O VII, O VIII, Fe XVII, Ne IX, and Ne X are always present in a proportion that is difficult to predict unless the geometry of the gas distribution is specified. We discuss the NEI scenario in a forthcoming publication. There we adopt the shock-driven scenario for NPS, consider the implications for the X-ray spectra, and compare the two scenarios.

If the model outlined above is correct, NPS and similar structures offer a possibility to probe the rotation pattern of the hot gas in the Milky Way. Future large grasp microcalorimetric missions like Line Emission Mapper (LEM, Kraft et al. 2022) capable of mapping the entire NPS region in the course of its all-sky survey (Khabibullin et al. 2023a) will be instrumental in measuring abundances and velocities using emission lines of elements from carbon to iron across the 0.2-2 keV energy band to narrow the range of plausible models.

## 5. Conclusions

While the North Polar Spur is usually attributed to a shock front associated with the activity of our Galactic Center or a local SNR, we discuss an alternative scenario. It posits that metalenriched plumes rise above the disk from active star-forming regions. Interactions with the hot halo gas give these plumes the appearance of bent spirals. Their shapes sensitively depend on the rotation pattern of the hot gas above the disk. These plumes are mostly hotter than the ambient gas and can be metal-rich compared to the hot gas in the Milky Way halo. In this model, the NPS itself is associated with the star formation within a 3-5

kpc distance from the Galactic Center. Fainter plumes might be associated with other star-forming regions in the Galaxy.

## Acknowledgments

This work is partly based on observations with the eROSITA telescope onboard SRG space observatory. The SRG observatory was built by Roskosmos in the interests of the Russian Academy of Sciences represented by its Space Research Institute (IKI) in the framework of the Russian Federal Space Program, with the participation of the Deutsches Zentrum für Luft- und Raumfahrt (DLR). The eROSITA X-ray telescope was built by a consortium of German Institutes led by MPE, and supported by DLR. The SRG spacecraft was designed, built, launched, and is operated by the Lavochkin Association and its subcontractors. The science data are downlinked via the Deep Space Network Antennae in Bear Lakes, Ussurĳsk, and Baikonur, funded by Roskosmos.

The development and construction of the eROSITA X-ray instrument was led by MPE, with contributions from the Dr. Karl Remeis Observatory Bamberg &amp; ECAP (FAU ErlangenNuernberg), the University of Hamburg Observatory, the Leibniz Institute for Astrophysics Potsdam (AIP), and the Institute for Astronomy and Astrophysics of the University of Tübingen, with the support of DLR and the Max Planck Society. The Argelander Institute for Astronomy of the University of Bonn and the Ludwig Maximilians Universität Munich also participated in the science preparation for eROSITA. The eROSITA data were processed using the eSASS/NRTA software system developed by the German eROSITA consortium and analyzed using proprietary data reduction software developed by the Russian eROSITA Consortium.

IK acknowledges support by the COMPLEX project from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program grant agreement ERC-2019-AdG 882679.

## References

Asplund, M., Grevesse, N., Sauval, A. J., &amp; Scott, P. 2009, ARA&amp;A, 47, 481 Barros, D. A., Lépine, J. R. D., &amp; Dias, W. S. 2016, A&amp;A, 593, A108 Baumgartner, V. &amp; Breitschwerdt, D. 2013, A&amp;A, 557, A140

Beck, R. 2015, A&amp;A Rev., 24, 4

Berkhuĳsen, E. M. 1971, A&amp;A, 14, 359

Bland-Hawthorn, J. &amp; Gerhard, O. 2016, ARA&amp;A, 54, 529

- Bunner, A. N., Coleman, P. L., Kraushaar, W. L., &amp; McCammon, D. 1972, ApJ, 172, L67
- Bykov, A. M. 2014, A&amp;A Rev., 22, 77
- Churazov, E., Forman, W., Jones, C., &amp; Böhringer, H. 2000, A&amp;A, 356, 788
- Churazov, E. &amp; Inogamov, N. 2004, MNRAS, 350, L52

Churazov, E. M., Khabibullin, I. I., Bykov, A. M., et al. 2021, MNRAS, 507, 971

Clarke, J. P. &amp; Gerhard, O. 2022, MNRAS, 512, 2171

- de Avillez, M. A. 2000, MNRAS, 315, 479

De Buizer, J. M., Lim, W., Karnath, N., &amp; Radomski, J. T. 2024, ApJ, 963, 55

Dickinson, C. 2018, Galaxies, 6, 56

- Drew, J. E., Monguió, M., &amp; Wright, N. J. 2019, MNRAS, 486, 1034
- Egger, R. J. &amp; Aschenbach, B. 1995, A&amp;A, 294, L25

Faucher-Giguère, C.-A. &amp; Oh, S. P. 2023, ARA&amp;A, 61, 131

- Foster, A. R., Ji, L., Smith, R. K., &amp; Brickhouse, N. S. 2012, ApJ, 756, 128

Fukui, Y., Inoue, T., Hayakawa, T., &amp; Torii, K. 2021, PASJ, 73, S405

- Fukui, Y., Ohama, A., Hanaoka, N., et al. 2014, ApJ, 780, 36
- Gull, S. F. &amp; Northover, K. J. E. 1973, Nature, 244, 80
- H. E. S. S. Collaboration, Abdalla, H., Abramowski, A., et al. 2018, A&amp;A, 612, A1
- Hanbury Brown, R., Davies, R. D., &amp; Hazard, C. 1960, The Observatory, 80, 191
- Haslam, C. G. T., Klein, U., Salter, C. J., et al. 1981, A&amp;A, 100, 209
- Heiles, C. 1990, ApJ, 354, 483
- Hodges-Kluck, E. J., Miller, M. J., &amp; Bregman, J. N. 2016, ApJ, 822, 21
- Kataoka, J., Sofue, Y., Inoue, Y., et al. 2018, Galaxies, 6, 27
- Kataoka, J., Tahara, M., Totani, T., et al. 2013, ApJ, 779, 57
- Kataoka, J., Yamamoto, M., Nakamura, Y., et al. 2021, ApJ, 908, 14
- Khabibullin, I., Churazov, E., &amp; Sunyaev, R. 2022, MNRAS, 509, 6068

Khabibullin,

I.,

Galeazzi,

M.,

Bogdan,

A., et

al.

2023a, arXiv

e-prints, arXiv:2310.16038

- Khabibullin, I. I., Churazov, E. M., Bykov, A. M., Chugai, N. N., &amp; Sunyaev, R. A. 2023b, MNRAS, 521, 5536
- Khabibullin, I. I., Churazov, E. M., Chugai, N. N., et al. 2024, arXiv e-prints, arXiv:2401.17261
- Kim, C.-G., Ostriker, E. C., Somerville, R. S., et al. 2020, ApJ, 900, 61
- Kobayashi, C., Karakas, A. I., &amp; Lugaro, M. 2020, ApJ, 900, 179
- Kobayashi, C., Umeda, H., Nomoto, K., Tominaga, N., &amp; Ohkubo, T. 2006, ApJ, 653, 1145
- Kraft, R., Markevitch, M., Kilbourne, C., et al. 2022, arXiv e-prints, arXiv:2211.09827

Lallement, R. 2023, Comptes Rendus Physique, 23, 1 LaRocca, D. M., Kaaret, P., Kuntz, K. D., et al. 2020, ApJ, 904, 54 Li, J. J., Immer, K., Reid, M. J., et al. 2022, ApJS, 262, 42 Mac Low, M.-M. &amp; McCray, R. 1988, ApJ, 324, 776 Markevitch, M. &amp; Vikhlinin, A. 2007, Phys. Rep., 443, 1 Motte, F., Nony, T., Louvet, F., et al. 2018, Nature Astronomy, 2, 478 Mou, G., Sun, D., Fang, T., et al. 2023, Nature Communications, 14, 781 Nguyen Luong, Q., Motte, F., Schuller, F., et al. 2011, A&amp;A, 529, A41 Norman, C. A. &amp; Ikeuchi, S. 1989, ApJ, 345, 372 Planck Collaboration, Akrami, Y., Ashdown, M., et al. 2020, A&amp;A, 641, A4 Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&amp;A, 647, A1 Predehl, P., Sunyaev, R. A., Becker, W., et al. 2020, Nature, 588, 227 Remazeilles, M., Dickinson, C., Banday, A. J., Bigot-Sazy, M. A., &amp; Ghosh, T. 2015, MNRAS, 451, 4311

- Sarkar, K. C. 2024, A&amp;A Rev., 32, 1
- Sarkar, K. C., Nath, B. B., &amp; Sharma, P. 2015, MNRAS, 453, 3827
- Schulreich, M. M. &amp; Breitschwerdt, D. 2022, MNRAS, 509, 716
- Shimoda, J. &amp; Asano, K. 2024, ApJ, 973, 78
- Sofue, Y. 1977, A&amp;A, 60, 327
- Sofue, Y. 1994, ApJ, 431, L91
- Sofue, Y. &amp; Kataoka, J. 2021, MNRAS, 506, 2170

Stein, Y., Dettmar, R. J., Beck, R., et al. 2020, A&amp;A, 639, A111

Su, M., Slatyer, T. R., &amp; Finkbeiner, D. P. 2010, ApJ, 724, 1044

- Sun, X. H., Landecker, T. L., Gaensler, B. M., et al. 2015, ApJ, 811, 40

Sunyaev, R., Arefiev, V., Babyshkin, V., et al. 2021, A&amp;A, 656, A132

Sutherland, R. S. &amp; Dopita, M. A. 1993, ApJS, 88, 253

Tomisaka, K. 1998, MNRAS, 298, 797

Tunmer, H. 1958, Philosophical Magazine, 3, 370

Woosley, S. E. &amp; Weaver, T. A. 1995, ApJS, 101, 181

Yamamoto, M., Kataoka, J., &amp; Sofue, Y. 2022, MNRAS, 512, 2034

- Yang, H. Y. K., Ruszkowski, M., &amp; Zweibel, E. G. 2022, Nature Astronomy, 6, 584
- Yang, R.-Z. &amp; Wang, Y. 2020, A&amp;A, 640, A60
- Zeng, Y., Wang, Q. D., &amp; Fraternali, F. 2023, MNRAS, 526, 483
- Zhang, B., Moscadelli, L., Sato, M., et al. 2014, ApJ, 781, 89
- Zhang, C., Churazov, E., &amp; Schekochihin, A. A. 2018, MNRAS, 478, 4785
- Zuhone, J. A. &amp; Roediger, E. 2016, Journal of Plasma Physics, 82, 535820301

Table B.1. A set of "sources" used to illustrate typical trajectories of gaseous plumes in Fig. 3. In the first two columns, the apparent Galactic Coordinates of the current source position are given. The next two columns give the distance (in kpc) from the Earth and the Galactic center. These positions coupled with equations 1-2 were used to generate trajectories shown in Fig. 3.

| 𝑙     | 𝑏   |   𝐷 Earth |   𝐷 GC | Comments            |
|-------|-----|-----------|--------|---------------------|
| 23 ◦  | 0 ◦ |       6   |    3.6 |                     |
| 23 ◦  | 0 ◦ |      12   |    5.4 | ∼ 5-kpc-ring, far   |
| 0 ◦   | 0 ◦ |       3   |    5.3 | ∼ 5-kpc-ring, near  |
| 0 ◦   | 0 ◦ |      13   |    4.7 | ∼ 5-kpc-ring, far   |
| 332 ◦ | 0 ◦ |      10   |    4.7 | ∼ 5-kpc-ring, far   |
| 332 ◦ | 0 ◦ |       3   |    5.8 | ∼ 5-kpc-ring, near  |
| 82 ◦  | 0 ◦ |       1.4 |    8.2 | Cygnus region, near |
| 82 ◦  | 0 ◦ |       3.4 |    8.5 | Cygnus region, far  |
| 262 ◦ | 0 ◦ |       1   |    8.5 | Vela region, near   |
| 262 ◦ | 0 ◦ |       5   |   10   | Vela region, far    |
| 45 ◦  | 0 ◦ |       3   |    6.5 |                     |
| 45 ◦  | 0 ◦ |       9   |    6.7 |                     |

## Appendix A: X-ray data analysis

Data from all four consecutive scans are combined together after filtering for the periods of enhanced solar activity, which causes a strongly elevated level of the instrumental background. For the imaging analysis, the data taken with all seven Telescope Modules (TMs) are combined, while for the spectral analysis, only five TMs protected by the on-chip filter (i.e. TMs 1-4, and 6) are used. Data reduction, filtering, vignetting-correction, and background subtraction are performed in the same way as was done in the previous studies exploring Galactic diffuse X-ray sources (Churazov et al. 2021; Khabibullin et al. 2022, 2023b, 2024). In particular, the energy-dependent contribution of the instrumental background is modeled and subtracted based on the calibration data accumulated via observations with the "Closed Filter Wheel" configuration, while corrections for exposure time and vignetting are conducted so that the data are characterized by Field-of-View averaged response matrices.

## Appendix B: Plumes in 3D

In this section, we outline two simple scenarios for plume trajectories. In the first scenario, a plume is produced by a continuous source. The plume initially moves together with stars/gas in the disk and eventually joins the motion of the halo gas. In the second scenario, a trajectory of a massless buoyant bubble is considered.

## Appendix B.1: Plumes from a moving source

The first case corresponds to the model described by equations 1 and 2. Fig. B.1 illustrates the shapes of the "plumes" in a stationary 3D frame. The blue curve shows the case when the source of the plume moves with the same velocity as the gas in the disk. As the plume rises, it gradually joints the slowly rotating gas in the disk. The red curve shows modest modifications that appear when the source of the plume moves slower than the gas and stars in the disk. By design, the rising plumes retain the same distance from the Galactic Center in projection to the disk.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-e2a68163-cca7-44e4-b833-32a3e106548f.jpg)

-4

Fig. B.1. An example of a 3D trajectory of a plume rising above the disk according to equations 1 and 2 and integrated over 140 Myr. Only one side of the plume is shown. The black circle depicts a circle in the disk plane with a radius of ∼ 5 kpc. The blue line corresponds to case (a), namely the source of the plume moves together with the gas/stars in the disk plane. The red curve is the case (b) when the source of the plume, i.e. an area of active star formation, moves relative to the gas. At this distance from the GC, the difference between these two cases is not large.

## Appendix B.2: Trajectory of a massless buoyant bubble

Here we consider the trajectory of a single (massless) bubble moving in a stratified atmosphere of the Galaxy. Similarly to galaxy clusters (Gull &amp; Northover 1973; Churazov et al. 2000), the velocity of the bubble is set by the balance of pressure gradients and the drag force acting on the bubble, i.e.

$$u ^ { 2 } \frac { \vec { u } } { u } \frac { 1 } { 2 } \rho A C _ { d } = - \nabla P V,$$

where fi 𝑢 is the 3D velocity vector of the bubble velocity with respect to the halo gas ( 𝑢 = | fi| 𝑢 ), 𝜌 and 𝑃 are the halo gas density and pressure, respectively, 𝐴 and 𝑉 are the area and volume of the bubble, and 𝐶 𝑑 is the appropriate drag coefficient (e.g. Zhang et al. 2018). The bubble size is assumed to be smaller than the pressure scale height and the velocity remains subsonic. For an atmosphere in hydrostatic equilibrium, 1 / 𝜌 ∇ 𝑃 = -∇ 𝜙 and, therefore, fi 𝑢 can be directly calculated from the gravitational potential 𝜙 , which in our case can include a contribution from halo gas rotation. To this end, we use the approximation of the MW potential from Barros et al. (2016) to which we add a centrifugal term ( -𝑉 2 ℎ ln 𝑅 ) for a cylindrical rotation, where 𝑅 is the distance from the MW center in the disk plane and 𝑉 ℎ is the gas rotation velocity. The rotation velocity of the gas is added to fi 𝑢 to get the bubble velocity in the non-rotating frame. The effects of rotation on the bubble trajectory are illustrated in Fig. B.2. This figure shows that different configurations are possible depending on the gas rotation pattern. The three solutions shown in the figure, cover the plausible range of halo rotations. For a non-rotating halo (green curve) the bubble moves along the potential gradient created by the Milky Way mass distribution. The blue and

11

9:

10

Fig. B.2. 3Dtrajectory of a (single) small and massless bubble rising under the action of buoyancy in a rotating but otherwise static atmosphere. The initial bubble is released at a distance of 5 kpc from the center (just above the disk plane). The green line shows the case of a non-rotating atmosphere. In this case, the bubble moves away from the Galaxy center in the disk and then switches to a more radial trajectory. The blue curve shows the case of a slowly rotating halo ( 𝑉 ℎ ∼ 100 km s -1 ). The bubble is now involved in rotation and motion towards larger radii. Finally, the red curve illustrates the case of fast rotation of the halo gas ( 𝑉 ℎ ∼ 200 km s -1 ). In this case, the centrifugal force is strong, and an inverted pressure gradient pushes the bubble closer to the rotation axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-7e6100dd-59b3-4a9d-8e5d-81925d47bbed.jpg)

red lines cover the cases of moderate (100 km s -1 ) and extreme (200 km s -1 ) halo rotation speed. In particular, the case of the halo rotation velocity with ∼ 180 km s -1 (Hodges-Kluck et al. 2016) should be more close to the red curve. Of course, the rotation pattern can be much more complicated than the cylindrical pattern considered here.

We also emphasize that Fig. B.2 (unlike Fig. B.1) shows trajectories of individual "bubbles" rather than "plumes" that are formed by many bubbles released at different times. Assuming that the source of the bubbles is moving together with stars in the disk, the plume trajectories (a combination of bubbles released at different times) are shown in Fig. B.3. Qualitatively, they resemble the plumes shown in Fig. B.1, except for the presence of radial migration.

## Appendix C: Enhanced metallicity of plumes

The enhanced metallicity of the NPS hot material enriched by supernovae is important for the proposed scenario. The thermal energy of the NPS 𝐸 th ∼ 2 × 10 55 erg (see §4) suggests the explosion of 𝑁 sn ∼ 2 × 10 4 core-collapse supernovae (CCSN) with a typical kinetic energy of 10 51 erg (lifetime of a star-forming region is short for SNe Ia to explode). We rely on two sets of oxygen nucleosynthesis calculations for CCSN progenitors in the range of 11 - 40 𝑀 ⊙ by Woosley &amp; Weaver (1995) (WW95) and Kobayashi et al. (2006) (K06). In K06, the data on 11 𝑀 ⊙ progenitors are lacking. We, therefore, adopt the ejected mass of oxygen from WW95 for this particular progenitor. The abundance of other elements relative to oxygen is assumed to follow (approximately) the solar composition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-55327b0e-71b7-4d0e-840a-8b9e8e1ffd6b.jpg)

10

Fig. B.3. 3D trajectory of a sequence of bubbles (i.e., a plume). Individual bubbles follow the trajectories shown in Fig. B.2. The source of the bubbles moves with the same velocity as stars in the disk. Overall, plume trajectories are qualitatively similar to those shown in Fig. B.3.

The average mass of ejected oxygen per CCSN ( 𝑚 O) is inferred assuming Salpeter initial mass function 𝑑𝑁 𝑑𝑚 / ∝ 𝑚 -2 35 . in the mass range of 0.1 - 100 𝑀 ⊙ . We find comparable values of 𝑚 O, 2.2 𝑀 ⊙ and 2.6 𝑀 ⊙ , for WW95 and K06 data, respectively. The average 𝑚 O = 2 4 . 𝑀 ⊙ multiplied by the CCSN number 𝑁 𝑠𝑛 provides us with the total amount of oxygen synthesized by supernovae of the star-forming region 𝑀 O ≈ 5 × 10 4 𝑀 ⊙ . The synthesized oxygen mixed with the gas of the solar composition (solar abundance 𝑋 ( O ) ⊙ = 0 01) produces enhanced overall metallic-. ity. The total mass of the mixture ( 𝑀 hot ) is fixed by the thermal energy of the X-ray-emitting gas 𝐸 th = ( 5 2 / ) 𝑘𝑇 𝑀 x hot /( 𝜇𝑚 𝑝 ) . For the average molecular weight 𝜇 = 0 61 . and temperature 𝑘𝑇 x = 0 6 keV one obtains . 𝑀 hot = 4 × 10 6 𝑀 ⊙ . The expected oxygen abundance of the NPS hot gas is then 𝑋 ( O ) = 𝑋 ( O ) ⊙ + 𝑀 O / 𝑀 hot ≈ 0 02, twice the solar value. .

The effect of the enhanced oxygen abundance of the NPS can be expressed formally as 𝑋 ( O ) = 𝑋 ( O ) ⊙ ( 1 + 𝑦 ) with the factor 𝑦 responsible for the synthesized oxygen and estimated for NPS to be about unity.

Let the abundance of a certain element (El) relative to oxygen produced by CCSNe with respect to the solar abundance is 𝜙 =(El/O)/(El/O) ⊙ . The overabundance factor of this element with respect to solar in the NPS material, assuming the initial abundance in the star-forming region to be solar, is then 𝑓 = 1 + 𝑦𝜙 .

The relative abundances 𝜙 of key elements, viz. C, N, Ne, Mg, Si, S, and Fe, generated by the CCSNe, can be derived from spectroscopic stellar data of [El/Fe] = log(El/Fe) - log(El/Fe) ⊙ for metal-poor stars ([Fe/H] &lt; -2), in which metallicity domain the nucleosynthesis by CCSNe dominates. We rely on [El/Fe] vs. [Fe/H] data compiled by Kobayashi et al. (2020). For C and N the scatter of [El/Fe] values relative to the average for different stars is very large, of ∼ 0.5 dex; in other cases, the scatter is of 0.2-0.3 dex. Anyway, for each element we obtain an average value of [El/Fe] with a typical error of 0.1 dex, which is converted to 𝜙 with a relative error of ≈ 20%. In the case of Ne, lacking stellar data, we rely on the theoretical prediction that CCSNe nucleosynthesis results in the solar ratio of Ne/O (Kobayashi et al. 2020). For the

adopted value of 𝑦 = 1, the estimated overabundance factors in the NPS material turn out to be 𝑓 = 1 5 for C, N, and S, . 𝑓 = 1 6 . for Mg and Si, 𝑓 = 2 for O and Ne, and 𝑓 = 1 25 in the case of . Fe - all factors with a relative error of approximately 20%.

## Appendix D: The chimney model

The clustering of core-collapse supernovae produced at the end of the evolution of massive stars in OB associations and compact clusters has a profound effect on the interstellar medium and the galactic ecology (Mac Low &amp; McCray 1988; Heiles 1990; de Avillez 2000; Kim et al. 2020). Multiple correlated supernovae and the powerful winds of OB stars create superbubbles and eventually (for powerful enough systems) produce the chimney-type structures (Norman &amp; Ikeuchi 1989). We consider such superbubbles as potential sources of the plumes. To study the effect of the galactic magnetic field on the superbubble breakthrough from the disk to halo Tomisaka (1998) performed 3D MHD simulations for different assumptions on the large-scale magnetic field structure. The magnetic field with the strength of ∼ 5 𝜇 G with a broad ∼ kpc scale height can confine a superbubble of a modest kinetic luminosity 3 × 10 37 erg s -1 for about 20 Myr within the scale height of | 𝑧 |∼ 300 pc from the disk, while in a model with the field that follows the scaling 𝐵 ∝ 𝜌 1 2 / (with the midplane field magnitude of 5 𝜇 G) the superbubble blows out to the halo. The evolution of superbubbles in density stratified disks that blow out into galactic halo may be a subject of Rayleigh-Taylor instabilities (Baumgartner &amp; Breitschwerdt 2013; Schulreich &amp; Breitschwerdt 2022). A parsec-resolution multi-phase simulation of the local star-forming galactic disk with the account for SNe feedback effects revealed that the hot galactic outflows (with gas temperatures above 10 6 K) may carry about 10%-20% of the energy and 30%-60% of the metal mass injected by supernovae (Kimetal. 2020). These values are broadly consistent with the energy requirements needed to produce the X-ray plumes discussed in the text.

Non-thermal particles accelerated by shocks from SNe and stellar winds at the active phase of the evolution of superbubble which is about 10 Myrs (Bykov 2014). Then relativistic particles will be blown out to the low halo with the frozen-in magnetic fields of the chimney-type plasma outflow. Relativistic electrons of energy ∼ 100 GeV would have in the halo a lifetime of about 10 Myrs due to synchrotron - Compton losses and will radiate in magnetic fields of a few 𝜇 G at frequencies of &lt; ∼ 50 GHz. The synchrotron-Compton cooling time for the relativistic electron of energy E GeV (measured in GeV) in the galactic halo can be estimated as

$$t _ { \text{syn} } & \sim 1 0 ^ { 9 } \frac { \mathcal { E } _ { \text{GeV} } ^ { - 1 } } { 1 + 0. 1 B _ { \mu G } ^ { 2 } } \, y r. & ( D. 1 ) & \text{dobe $D$} \\ & \text{the core} \\ & \text{high-ma}.$$

## Appendix E: The role of W43

Below, we discuss in brief the mini-starburst complex W 43 as a prototype powerhouse located in the Galactic molecular ring assuming that an ensemble of such complexes may supply the plumes over a long period of ∼ hundred million years. The molecular complex W 43 with the estimated mass ∼ 7 × 10 6 𝑀 ⊙ is presumably located at the connecting point of the ScutumCentaurus Galactic arm and the Galactic bar (Nguyen Luong et al. 2011). W 43 excited a giant H II region dubbed NGC 3603. The equivalent diameter of the coherent complex of molecular clouds is ∼ 140 pc and it is surrounded by an atomic gas envelope

Fig. D.1. Possible formation scenario of a NPS-like structure. During the initial phase after the breakout (left), a shock driven by the hot and metal-rich gas of a superbubble propagates through the ambient gas. At later stages (right) the gas is stretched and bent into a plume but possibly maintains the layered structure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00752v2-figure-ce387973-ba38-4fdc-8d03-a278b88ea6a7.jpg)

of about two times larger diameter. Nguyen Luong et al. (2011) noted that the transition from circular to elliptical orbits in the spiral arm and bar potentials might cause high-velocity streams and thus efficient star formation episodes of current rate at least ∼ 0 01 . 𝑀 ⊙ yr -1 as estimated from 8 𝜇 mluminosity measurements but it is likely increasing to ∼ 0 1 . 𝑀 ⊙ yr -1 . On longer time scales ∼ 100 Myrs of interest for this paper multiple episodes of high star formation rate can be expected from colliding molecular clouds entrained by the large scale streams.

The giant HII region NGC 3603 which is associated with the complex produces the Lyman continuum photon rate of about 4 × 10 51 photons s -1 is the most powerful HII region in the Galaxy (De Buizer et al. 2024). A compact dense cluster of young massive stars HD 97950 (also known as W 43 cluster) powering NGC 3603 has an estimated dynamic stellar mass of 18,000 𝑀 ⊙ including more than 70 O and 3 WR stars. The distance to the powerful starburst region NGC 3603 of 7.2 ± 0.1 kpc was estimated by Drew et al. (2019) from GAIA DR2 stellar parallax measurements (see also De Buizer et al. (2024) for a recent discussion). From the VLBA trigonometric parallax measurements of masers toward W43 Zhang et al. (2014) derived a distance of 5.49 + 0 39 . -0 34 . kpc to the complex (see also Li et al. 2022, for another distance estimate of ∼ 4 8 kpc). The origin and ages of molecular . clumps and stellar population in NGC 3603 are still the subject of debate De Buizer et al. (2024). ALMA observations revealed that the core mass function in W43-MM1 is much shallower in the high-mass range than the standard initial mass function (Motte et al. 2018). Fukui et al. (2014, 2021) discussed the formation of cluster W 43 due to a collision of two molecular clouds about 1 Myr ago. This age is consistent with that of the highest mass stellar population in the cluster. The suggestion by Nguyen Luong et al. (2011) about the origin of the star-forming complex W 43/NGC 3603 due to the collision of molecular clouds following their circular orbits in the Scutum-Centaurus arm and the clouds at the elliptical orbits of the Galactic bar makes the complex W 43 as well as other complexes along the Molecular Ring where the Galactic arms meet the bar to be favorable locations for the intense star formation episodes in the past.

Very high energy emission source HESS J1848-018 (H. E. S. S. Collaboration et al. 2018) can possibly be associ-

ated with W 43. From Fermi-LAT detection of extended high energy gamma-ray source source in the direction of W 43 Yang &amp;Wang (2020) estimated the total cosmic ray energy within the source to explain the detected flux to be ( 2 3 . ± 0 3 . ) × 10 48 erg.
E. Churazov I. Khabibullin A. M. Bykov N. N. Chugai R. A. Sunyaev V. P. Utrobin I. I. Zinchenko
North Polar Spur: gaseous plume(s) from star-forming regions at $\sim$3-5 kpc from Galactic Center?
We argue that the North Polar Spur (NPS) and many less prominent structures
are formed by gaseous metal-rich plumes associated with star-forming regions
(SFRs). The SFRs located at the tangent to the 3-5~kpc rings might be
particularly relevant to NPS. A multi-temperature mixture of gaseous components
and cosmic rays rises above the Galactic disk under the action of their initial
momentum and buoyancy. Eventually, the plume velocity becomes equal to that of
the ambient gas, which rotates with different angular speed than the stars in
the disk. As a result, the plumes acquire characteristic bent shapes. An ad hoc
model of plumes' trajectories shows an interesting resemblance to the
morphology of structures seen in the radio continuum and X-rays.