A hybrid machine learning model which combines a shallow convolutional neural network and a long short-term memory network (CNN--LSTM), has been developed to automate the detection of kink oscillations in coronal plasma loops within large volumes of high-cadence sequences of imaging data. The network was trained on a set of 10,000 synthetic data cubes designed to mimic sequences of coronal images, achieving an accuracy greater than 98\% on this synthetic dataset. The model was then applied to detect kink oscillations in real data cubes of coronal active regions observed with SDO/AIA in the 171~Å channel. This dataset consisted of 50 samples with visually detected kink oscillations and 128 samples without. Each sample covered an area of 260$\times$260~pixels in the spatial domain and a duration of 30~min with a 12~s cadence in the time domain. Both off-limb and on-disk regions of interest were used. The data were pre-processed by median filtering in the time domain, and Gaussian smoothing and Contrast Limited Adaptive Histogram Equalization in the spatial domain. In the real dataset, the performance of the model was 83.7\%.The model is fully available in open access. We regard the CNN--LSTM model developed as a first step toward creating robust tools for routine solar coronal data mining in the context of coronal oscillation study.

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis-Chandrasekhar-Fermi (DCF) technique. However, we identify the major problem of the DCF to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF, the new technique uses not the dispersion of the polarization angle and line of sight velocities, but increments of these quantities given by the structure functions. To address the variety of the astrophysical conditions for which our technique can be applied, we consider the turbulence in both media with magnetic pressure larger than the gas pressure corresponding e.g. to molecular and the gas pressure larger than the magnetic pressure corresponding to the warm neutral medium. We provide general expressions for arbitrary admixture of Alfv\'en, slow and fast modes in these media and consider in detail the particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our Differential Measure Approach (DMA), unlike the DCF, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data and is much more stable in the presence of large scale variations induces of non-turbulent nature. In parallel, our study uncover the deficiencies of the earlier DCF research.

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals \Delta[X/H] from this two-parameter fit. The rms residuals range from ~0.01-0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow comparison of populations controlled for differences in metallicity and [\alpha/Fe]. Relative to the main disk (R=3-13 kpc, |Z|<2 kpc), we find nearly identical abundance patterns in the outer disk (R=15-17 kpc), 0.05-0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4-1 dex) of multiple elements in \omega Cen. Residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

There exist isolated elliptical galaxies, whose dynamics can be modelled without resorting to dark matter or MOND, e.g. NGC 7507. The isolated elliptical NGC 5812 is another object to investigate a possible role of isolation. We use globular clusters (GCs) and the galaxy light itself as dynamical tracers to constrain its mass profile. We employ Gemini/GMOS mask spectroscopy, apply the GMOS reduction procedures provided within IRAF, measure GC velocities by cross correlation methods and extract the line-of-sight kinematics of galaxy spectra using the tool pPXF. We identify 28 GCs with an outermost galactocentric distance of 20 kpc, for which velocities could be obtained. Furthermore, 16 spectra of the integrated galaxy light out to 6 kpc have been used to model the central kinematics. These spectra provide evidence for a disturbed velocity field. We construct spherical Jeans models for the galaxy light and apply tracer mass estimators for the globular clusters. With the assumptions inherent to the mass estimators, MOND is compatible with the mass out to 20 kpc. However, a dark matter free galaxy is not excluded. We find one globular cluster with an estimated mass of 1.6x10**7 solar masses, the first Ultra Compact Dwarf in an isolated elliptical. We put NGC 5812 into the general context of dark matter or alternative ideas in elliptical galaxies. The case for a MONDian phenomenology also among early-type galaxies has become so strong that deviating cases appear astrophysically more interesting than agreements. The baryonic Tully Fisher relation (BTFR) as predicted by MOND is observed in some samples of early-type galaxies, in others not. However, in cases of galaxies that deviate from the MONDian prediction, data quality and data completeness are often problematic. (abridged)

We present the first detailed chemical analysis from APOGEE-2S observations of stars in six regions of recently discovered substructures in the outskirts of the Magellanic Clouds extending to 20 degrees from the LMC center. We also present, for the first time, the metallicity and alpha-abundance radial gradients of the LMC and SMC out to 11 degrees and 6 degrees, respectively. Our chemical tagging includes 13 species including light, alpha, and Fe-peak elements. We find that the abundances of all of these chemical elements in stars populating two regions in the northern periphery - along the northern "stream"-like feature - show good agreement with the chemical patterns of the LMC, and thus likely have an LMC origin. For substructures located in the southern periphery of the LMC, we find more complex chemical and kinematical signatures, indicative of a mix of LMC-like and SMC-like populations. However, the southern region closest to the LMC shows better agreement with the LMC, whereas that closest to the SMC shows a much better agreement with the SMC chemical pattern. When combining this information with 3-D kinematical information for these stars, we conclude that the southern region closest to the LMC has likely an LMC origin, whereas that closest to the SMC has an SMC origin, and the other two southern regions have a mix of LMC and SMC origins. Our results add to the evidence that the southern substructures of the LMC periphery are the product of close interactions between the LMC and SMC, and thus likely hold important clues that can constrain models of their detailed dynamical histories.

We study the bar pattern speeds and corotation radii of 225 barred galaxies, using IFU data from MaNGA and the Tremaine-Weinberg method. Our sample, which is divided between strongly and weakly barred galaxies identified via Galaxy Zoo, is the largest that this method has been applied to. We find lower pattern speeds for strongly barred galaxies than for weakly barred galaxies. As simulations show that the pattern speed decreases as the bar exchanges angular momentum with its host, these results suggest that strong bars are more evolved than weak bars. Interestingly, the corotation radius is not different between weakly and strongly barred galaxies, despite being proportional to bar length. We also find that the corotation radius is significantly different between quenching and star forming galaxies. Additionally, we find that strongly barred galaxies have significantly lower values for R, the ratio between the corotation radius and the bar radius, than weakly barred galaxies, despite a big overlap in both distributions. This ratio classifies bars into ultrafast bars (R < 1.0; 11% of our sample), fast bars (1.0 < R < 1.4; 27%) and slow bars (R > 1.4; 62%). Simulations show that R is correlated with the bar formation mechanism, so our results suggest that strong bars are more likely to be formed by different mechanisms than weak bars. Finally, we find a lower fraction of ultrafast bars than most other studies, which decreases the recently claimed tension with {\Lambda}CDM. However, the median value of R is still lower than what is predicted by simulations.

Propagation of cosmic rays (CRs) in turbulent and magnetized astrophysical media is a long-standing problem that requires both understanding of the properties of turbulent magnetic fields and their interaction with energetic particles. This review focuses on selected recent theoretical findings made based on the progress in understanding and simulating magnetohydrodynamic (MHD) turbulence. In particular, we address the problem of perpendicular and parallel propagation of CRs and identify the conditions when the perpendicular propagation is superdiffusive and diffusive. For the parallel diffusion, we discuss the problems of the traditionally used diffusion mechanism arising from pitch angle scattering and the possible solutions provided by the recently identified ``mirror diffusion" in the presence of turbulent magnetic mirrors.

Several models have been proposed to describe the dynamics of epithelial tissues undergoing morphogenetic changes driven by apical constriction pulses, which differ in where the constriction is applied, either at the perimeter or medial regions. To help discriminate between these models, using the vertex model for epithelial dynamics, we analysed the impact of where the constriction is applied on the final geometry of the active cell that is reducing its apical size. We find that medial activity, characterised by a reduction in the reference area in the vertex model, induces symmetry breaking and generates anisotropic cell shapes, while isotropic cell shapes and larger contractions occur when the reference perimeter in the model is reduced. When plasticity is included, sufficiently slow processes of medial contractile activity, compared with typical apical constriction pulses, can also achieve significant cell contraction. Finally, we apply the model to describe the active apical contractile pulses observed during cellular mitotic events within the epithelial enveloping cell layer in the developing annual killifish Austrolebias nigripinnis, being able to quantitatively describe the temporal evolution of cell shape changes when perimeter activity and area plasticity are included. A global fit of all parameters of the vertex model is provided.

The challenge of consistent identification of internal structure in galaxies - in particular disc galaxy components like spiral arms, bars, and bulges - has hindered our ability to study the physical impact of such structure across large samples. In this paper we present Galaxy Zoo: 3D (GZ: 3D) a crowdsourcing project built on the Zooniverse platform which we used to create spatial pixel (spaxel) maps that identify galaxy centres, foreground stars, galactic bars and spiral arms for 29831 galaxies which were potential targets of the MaNGA survey (Mapping Nearby Galaxies at Apache Point Observatory, part of the fourth phase of the Sloan Digital Sky Surveys or SDSS-IV), including nearly all of the 10,010 galaxies ultimately observed. Our crowd-sourced visual identification of asymmetric, internal structures provides valuable insight on the evolutionary role of non-axisymmetric processes that is otherwise lost when MaNGA data cubes are azimuthally averaged. We present the publicly available GZ:3D catalog alongside validation tests and example use cases. These data may in the future provide a useful training set for automated identification of spiral arm features. As an illustration, we use the spiral masks in a sample of 825 galaxies to measure the enhancement of star formation spatially linked to spiral arms, which we measure to be a factor of three over the background disc, and how this enhancement increases with radius.

A recent publication (Li et al. 2021) discovered one of the largest filamentary neutral hydrogen features dubbed Cattail from high resolution FAST observations that might be a new galactic arm of our own Milky Way. However in the analysis, it was suggested that this neutral hydrogen feature is cold despite having 12km/s total linewidth. We evaluate the probability whether the Cattail is actually cold neutral media via the newly developed Velocity Decomposition Algorithm (Yuen et al. 2021a) and Force Balancing Model (Ho et al. 2021a). We discovered that even with the inclusion of the galactic shear term, the feature is still at the unstable neutral media regime. Moreover, we also discover that the Cattail is two disjoint features in caustics space, suggesting that the Cattail might have two different turbulent systems. We check the spectra of the individual system separated via VDA to confirm this argument. We do not exclude the existence of smaller scale cold media being embedded within this structure.

We investigate how the stellar rotational support changes as a function of spatially resolved stellar population age ($\rm D_n4000$) and relative central stellar surface density ($\Delta \Sigma_1$) for MaNGA isolated/central disk galaxies. We find that the galaxy rotational support $\lambda_{R_\mathrm{e}}$ varies smoothly as a function of $\Delta \Sigma_1$ and $\rm D_n4000$. $\rm D_n4000$vs.$\Delta \Sigma_1$ follows a "J-shape", with $\lambda_{R_\mathrm{e}}$ contributing to the scatters. In this "J-shaped" pattern rotational support increases with central $\rm D_n4000$ when $\Delta \Sigma_1$is low but decreases with$\Delta \Sigma_1$when$\Delta \Sigma_1$ is high. Restricting attention to low-$\Delta \Sigma_1$ (i.e, large-radius) galaxies, we suggest that the trend of increasing rotational support with $\rm D_n4000$ for these objects is produced by a mix of two different processes, a primary trend characterized by growth in $\lambda_{R_\mathrm{e}}$ along with mass through gas accretion, on top of which disturbance episodes are overlaid, which reduce rotational support and trigger increased star formation. An additional finding is that star forming galaxies with low $\Delta \Sigma_1$ have relatively larger radii than galaxies with higher $\Delta \Sigma_1$ at fixed stellar mass. Assuming that these relative radii rankings are preserved while galaxies are star forming then implies clear evolutionary paths in central $\rm D_n4000$ vs. $\Delta \Sigma_1$. The paper closes with comments on the implications that these paths have for the evolution of pseudo-bulges vs. classical-bulges. The utility of using $\rm D_n4000$-$\Delta \Sigma_1$ to study $\lambda_{R_\mathrm{e}}$ reinforces the notion that galaxy kinematics correlate both with structure and with stellar-population state, and indicates the importance of a multi-dimensional description for understanding bulge and galaxy evolution.

The APOGEE Open Cluster Chemical Abundances and Mapping (OCCAM) survey is used to probe the chemical evolution of the s-process element cerium in the Galactic disk. Cerium abundances were derived from measurements of Ce II lines in the APOGEE spectra using the Brussels Automatic Code for Characterizing High Accuracy Spectra (BACCHUS) in 218 stars belonging to 42 open clusters. Our results indicate that, in general, for Ages &lt; 4 Gyr, younger open clusters have higher [Ce/Fe] and [Ce/$\alpha$-element] ratios than older clusters. In addition, metallicity segregates open clusters in the [Ce/X]-Age plane (where X can be H, Fe, and the $\alpha$-elements O, Mg, Si, or Ca). These metallicity-dependant relations result in [Ce/Fe] and [Ce/$\alpha$] ratios with age that are not universal clocks. Radial gradients of [Ce/H] and [Ce/Fe] ratios in open clusters, binned by age, were derived for the first time, with d[Ce/H]dR$_{GC}$ being negative, while d[Ce/Fe]/dR$_{GC}$ is positive. [Ce/H] and [Ce/Fe] gradients are approximately constant over time, with the [Ce/Fe] gradient becoming slightly steeper, changing by $\sim$+0.009 dex-kpc$^{-1}$-Gyr$^{-1}$. Both the [Ce/H] and [Ce/Fe] gradients are shifted to lower values of [Ce/H] and [Ce/Fe] for older open clusters. The chemical pattern of Ce in open clusters across the Galactic disk is discussed within the context of s-process yields from AGB stars, $\sim$Gyr time delays in Ce enrichment of the interstellar medium, and the strong dependence of Ce nucleosynthesis on the metallicity of its AGB stellar sources.

Magnetic fields play a crucial role in star formation, yet tracing them becomes particularly challenging, especially in the presence of outflow feedback in protostellar systems. We targeted the star-forming region L1551, notable for its apparent outflows, to investigate the magnetic fields. These fields were probed using polarimetry observations from the Planck satellite at 353 GHZ/849 $\mu$m, the SOFIA/HAWC+ measurement at 214 $\mu$m, and the JCMT/SCUPOL 850 $\mu$m survey. Consistently, all three measurements show that the magnetic fields twist towards the protostar IRS 5. Additionally, we utilized the Velocity Gradients Technique (VGT) on the $^{12}$CO (J = 1-0) emission data to distinguish the magnetic fields directly associated with the protostellar outflows. These were then compared with the polarization results. Notably, in the outskirts of the region, these measurements generally align. However, as one approaches the center of IRS 5, the measurements tend to yield mostly perpendicular relative orientations. This suggests that the outflows might be dynamically significant from a scale of approximately $\sim0.2$~pc, causing the velocity gradient to change direction by 90 degrees. Furthermore, we discovered that the polarization fraction $p$ and the total intensity $I$ correlate as $p \propto I^{-\alpha}$. Specifically, $\alpha$ is approximately $1.044\pm0.06$ for SCUPOL and around $0.858\pm0.15$ for HAWC+. This indicates that the outflows could significantly impact the alignment of dust grains and magnetic fields in the L1551 region.

The recent exquisite Gaia astrometric, photometric, and radial velocity (RV) measurements resulted in a substantial advancement for the determination of the orbits for old star clusters, including the oldest Milky Way globular clusters (MW GCs). The main goal of this paper is to use the Gaia DR3 and the VVVX measurements to obtain the orbits for nearly a dozen new Galactic GC candidates that have been poorly studied or previously unexplored. We use the Gaia DR3 and VVVX databases to identify bonafide members of the Galactic GC candidates: VVV-CL160, Patchick122, Patchick125, Patchick126, Kronberger99, Kronberger119, Kronberger143, ESO92-18, ESO93-08, Gaia2, and Ferrero54. The relevant mean cluster physical parameters are derived (distances, Galactic coordinates, proper motions, RVs). We measure accurate mean RVs for the GCs VVV-CL160 and Patchick126, using observations acquired at the Gemini-South telescope with the IGRINS high-resolution spectrograph. Orbits for each cluster are then computed using the GravPot16 model, assuming typical Galactic bar pattern speeds. We reconstruct the orbits for these clusters for the first time. These include star clusters with retrograde and prograde orbital motions, both in the Galactic bulge and disk. Orbital properties, such as the mean time-variations of perigalactic and apogalactic distances, eccentricities, vertical excursions from the Galactic plane, and Z-components of the angular momentum are obtained for our sample. Our main conclusion is that, based on the orbital parameters, Patchick125 and Patchick126 are genuine MW bulge/halo GCs; Ferrero54, Gaia2 and Patchick122 are MW disk GCs. The orbits of Kronberger99, Kronberger119, Kronberger143, ESO92-18, and ESO93-08 are more consistent with old MW disk open clusters. VVV-CL160 falls very close to the Galactic centre, but reaches larger distances beyond the Sun, thus its origin is still unclear.

The goal of the Open Cluster Chemical Abundances and Mapping (OCCAM) survey is to constrain key Galactic dynamic and chemical evolution parameters by the construction and analysis of a large, comprehensive, uniform data set of infrared spectra for stars in hundreds of open clusters. This sixth contribution from the OCCAM survey presents analysis of SDSS/APOGEE Data Release 17 (DR17) results for a sample of stars in 150 open clusters, 94 of which we designate to be "high quality'' based on the appearance of their color-magnitude diagram. We find the APOGEE DR17-derived [Fe/H] values to be in good agreement with those from previous high resolution spectroscopic open cluster abundance studies. Using a subset of the high quality sample, the Galactic abundance gradients were measured for 16 chemical elements, including [Fe/H], for both Galactocentric radius ($R_{GC}$) and guiding center radius ($R_{Guide}$). We find an overall Galactic [Fe/H] vs $R_{GC}$ gradient of $-0.073 \pm 0.002$ dex/kpc over the range of 6 &lt; R_{GC} &lt; 11.5 kpc, and a similar gradient is found for [Fe/H] versus $R_{Guide}$. Significant Galactic abundance gradients are also noted for O, Mg, S, Ca, Mn, Na, Al, K and Ce. Our large sample additionally allows us to explore the evolution of the gradients in four age bins for the remaining 15 elements.

Chemical abundances are an essential tool in untangling the Milky Way's enrichment history. However, the evolution of the interstellar medium abundance gradient with cosmic time is lost as a result of radial mixing processes. For the first time, we quantify the evolution of many observational abundances across the Galactic disk as a function of lookback time and birth radius, $R_\text{birth}$. Using an empirical approach, we derive $R_\text{birth}$ estimates for 145,447 APOGEE DR17 red giant disk stars, based solely on their ages and [Fe/H]. We explore the detailed evolution of 6 abundances (Mg, Ca ($\alpha$), Mn (iron-peak), Al, C (light), Ce (s-process)) across the Milky Way disk using 87,426 APOGEE DR17 red giant stars. We discover that the interstellar medium had three fluctuations in the metallicity gradient $\sim 9$,$\sim 6$, and$\sim4$ Gyr ago. The first coincides with the end of high-$\alpha$ sequence formation around the time of the Gaia-Sausage-Enceladus disruption, while the others are likely related to passages of the Sagittarius dwarf galaxy. A clear distinction is found between present-day observed radial gradients with age and the evolution with lookback time for both [X/Fe] and [X/H], resulting from the significant flattening and inversion in old populations due to radial migration. We find the [Fe/H]--[$\alpha$/Fe] bimodality is also seen as a separation in the $R_\text{birth}$--[X/Fe] plane for the light and $\alpha$-elements. Our results recover the chemical enrichment of the Galactic disk over the past 12 Gyr, providing tight constraints on Galactic disk chemical evolution models.

A coronal loop of 290~Mm length, observed at 171~Å with SDO/AIA on February 6th 2024 near AR 13571, is found to oscillate with two significantly different oscillation periods, $48.8 \pm 6.1$~min and $4.8\pm 0.3$~min. The oscillations occur in the time intervals without detected flares or eruptions. Simultaneously, near the Northern footpoint of the oscillating loop, we detect a $49.6 \pm 5.0$-min periodic variation of the average projected photospheric magnetic field observed with SDO/HMI. The shorter-period decayless oscillation is attributed to the eigen-mode, standing kink oscillation of the loop, while the longer-period oscillation may be the oscillatory motion caused by the periodic footpoint driver. The photospheric long-period process can also drive the short-period, eigen oscillation of the loop via the self-oscillatory, \lq\lq violin\rq\rq\, mechanism, in which a transverse oscillation is excited by an external quasi-steady flow. This finding indicates that the most powerful, lower-frequency spectral components of photospheric motions, which are well below the Alfvénic/kink cutoff, can reach the corona.

The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R$\,$136 has been shown to be the main source of energy creating multiple pc-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214$\,\mu$m. The B-field morphology is complex, showing bending structures around R$\,$136. In addition, we use high spectral and angular resolution [\textsc{CII}] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B-field morphology and shows evidence of multiple expanding shells. Our B-field strength maps, estimated using the Davis-Chandrasekhar-Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350$\,\mu$G at 89, 154, and 214$\,\mu$m, respectively. We estimated that the majority of the 30 Doradus clouds are sub-critical and sub-Alfvénic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B-field is sufficient to hold the cloud structure integrity under feedback from R$\,$136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique (VGT) using [\textsc{CII}] and CO(2-1) is likely to confirm these results.

Polarization of interstellar dust emission is a powerful probe of dust properties and magnetic field structure. Yet studies of external galaxies are hampered by foreground dust contribution. The aim of this study is to separate the polarised signal from the Large Magellanic Cloud (LMC) from that of the Milky Way (MW) in order to construct a wide-field, spatially complete map of dust polarization using the Planck 353 GHz data. To estimate the foreground polarization direction, we used velocity gradients in HI spectral line data and assessed the performance of the output by comparing to starlight extinction polarization. We estimate the foreground intensity using dust properties derived previously from the Planck data and we assume the foreground polarization to be uniform and equal to the average of the MW around the galaxy. After foreground removal, the geometry of the plane-of-the-sky magnetic field in the LMC tends to follow the structure of the atomic gas in the LMC. This is notably the case along the molecular ridges extending south and south-east of the 30 Doradus star-forming complex, and along more diffuse southern arm extending towards the Small Magellanic Cloud. There is also an alignment between the magnetic field and the outer arm in the western part. The median polarization fraction in the LMC is slightly lower than that observed for the MW while the anti-correlation between the polarization angle dispersion function and the polarization fraction is slightly larger. Overall, polarization fraction distribution is similar to that observed in the MW.

We report the evidence for the existence of the universal, continuous turbulent cascade of velocity fluctuations with Kolmogorov -5/3 slope spanning 6 orders of length scales, from $10^4$ pc down to $10^{-2}$ pc. This was achieved by applying our innovative technique of separating density and velocity fluctuations to a set of spectroscopic surveys featuring various galactic spectral lines. This unified velocity cascade involves different interstellar phases from diffuse galactic media to dense self-gravitating clouds and persists despite interstellar phase transitions. However, the turbulent density fluctuations do not show this universality as the value of the spectral slope changes in different interstellar phases. This agrees with the expectation of compressible turbulence theory and demonstrates that the density is only an indirect tracer of interstellar turbulence. We report that the density fluctuations for clouds and filaments that are preferentially parallel to magnetic fields exhibit the spectral slope of -2. The universal grand turbulent velocity cascade that is established in our paper has significant implications for key galactic physical processes, including star formation, cosmic ray transport etc. We anticipate our result to be a starting point for in vitro models of multiphase interstellar turbulence studies with a significant impact for modeling of spiral galaxies.