ViCTORIA project: MeerKAT H I observations of the ram pressure stripped galaxy NGC 4523

pressure


Introduction
The environment plays a major role in shaping galaxy evolution.Galaxies inhabiting rich clusters are predominantly quiescent and relaxed systems, while those located in the field are star-forming, late-type objects (morphology segregation, Dressler 1980;Dressler et al. 1997;Whitmore et al. 1993).At the same time, late-type systems inhabiting dense regions are generally gas-poor and with a reduced star formation activity (e.g., Boselli & Gavazzi 2006).Several mechanisms able to affect the evolution of galaxies in high-density regions have been proposed in the literature to explain these differences.They include gravitational perturbations induced by the interaction with nearby companions (Merritt 1983), with the gravitational potential well of the cluster dark matter halo (Byrd & Valtonen 1990), or their combined effect after several high-speed fly-by encounters (galaxy harassment, Moore et al. 1998).Hydrodynamic interactions can also occur between the cold interstellar medium (ISM) of galaxies and the surrounding hot (T ∼ 10 7 -10 8 K) intracluster medium (ICM) trapped within the gravitational potential well of clusters and groups and emitting in X-rays (e.g., Sarazin 1986;Sun 2012).They include heat conduction between the two media at different temperatures (thermal evaporation, Cowie & Songaila 1977), mixing in the instabilities formed at the interface between the galaxy ISM and the surrounding ICM (e.g., viscous stripping, Nulsen 1982), the end of infall of fresh gas once a galaxy becomes a satellite of a larger system (starvation, Larson et al. 1980), and the pressure exerted by the ICM on galaxies moving at high velocity (500-1000 km s −1 ) within it (ram pressure stripping, Gunn & Gott 1972).All these mechanisms are able to perturb the gas content and distribution of galaxies and as a consequence alter their production of new stars (e.g., Boselli & Gavazzi 2006, 2014;Cortese et al. 2021;Boselli et al. 2022).Over long timescales, once the gas reservoir is removed from the disc, the activity of star formation is reduced, and galaxies become quiescent (e.g., Peng et al. 2010;Boselli et al. 2023).
While the general picture of galaxy evolution in rich environments is now well understood, several open questions still need to be answered.First of all, it is still unclear which of these mechanisms is dominant in different structures of the cosmic web, from filaments, loose and compact groups, to rich clusters of galaxies.Since the efficiency of all these processes is tightly connected to the properties of the high-density regions (gas density and temperature, galaxy number density, velocity dispersion, infall rate, etc.; Boselli & Gavazzi 2006), their relative weight can also significantly change with cosmic time.Furthermore, the effects of these mechanisms are still not fully understood.The major difficulties are related to the fact that they concern a multi-phase medium, from cold atomic and molecular gas to ionised and hot plasma, and it acts on a wide range of scales, from individual H ii regions and giant molecular clouds (GMC, 50 pc scales) to entire clusters of galaxies (1 Mpc scales).For the same reasons, hydrodynamic simulations of these perturbing mechanisms are also challenging (e.g., Roediger 2009).One particular aspect that deserves further investigation is the fate of the stripped gas.The most recent observations seem to indicate that the gas is principally stripped from the galaxy discs as cold atomic hydrogen and, once mixed with the surrounding ICM, it first becomes ionised and then becomes hot (e.g., Boselli et al. 2016aBoselli et al. , 2021Boselli et al. , 2022;;Poggianti et al. 2019;Campitiello et al. 2021;Sun et al. 2021; Bartolini et al. 2022), thus contributing to the pollution of the diffuse ICM (e.g., Longobardi et al. 2020a,b).Part of the molecular gas can also be removed during the interaction (Fumagalli et al. 2009;Boselli et al. 2014b;Zabel et al. 2019Zabel et al. , 2022)).The stripped gas can collapse into GMCs to form new stars (Jáchym et al. 2014(Jáchym et al. , 2017(Jáchym et al. , 2019;;Moretti et al. 2020a,b), although this process is not ubiquitous (e.g., Gavazzi et al. 2001;Boselli et al. 2016a;Pedrini et al. 2022).
The local Universe is a perfect target for environmental studies.First of all, it gives us the end point of the evolutionary path leading to local galaxies, and different samples selected according to different criteria can be used to drive strong statistical results.It also allows us to study at high physical resolution different samples at various wavelengths, giving us a unique view of the effects of the perturbations on the different galaxy components (gas in its different phases, dust, stars, cosmic rays, etc.).Finally, dwarf galaxies, the most fragile objects to external perturbations and to internal energy release because of their shallow gravitational potential well, are only here accessible at all frequencies.For these reasons the Virgo cluster, the closest cluster of galaxies, has been the target of several untargeted multifrequency surveys.These include GALEX UV observations (GALEX Ultraviolet Virgo Cluster Survey, GUViCS, Boselli et al. 2011), deep optical imaging (Next Generation Virgo cluster Survey, NGVS, Ferrarese et al. 2012), mid-(Wide-field Infrared Survey Explorer, WISE, Wright et al. 2010) and far-IR imaging (Herschel Virgo Cluster Survey, HeViCS, Davies et al. 2010).A deep untargeted Hα narrowband imaging of the cluster, especially designed to study the effects of the perturbations on the star formation activity of galaxies and to detect extended low surface brightness ionised gas features formed during the interactions of galaxies with their surrounding environment, has been almost completed (Virgo Environmental Survey Tracing Ionised Gas Emission, VESTIGE, Boselli et al. 2018a).Given their untargeted nature, these surveys are ideally designed to identify galaxies undergoing a perturbation and the very nature of the perturbing mechanism on a statistically complete sample, thus minimising any possible bias related to selection effects present in targeted observations.
The Virgo cluster has been covered also in the radio domain by the NRAO VLA Sky Survey (NVSS, Condon et al. 1998) at 1.4 GHz, and in the H i line by the H i Parkes All Sky Survey (HIPASS, Barnes et al. 2001) and Arecibo Legacy Fast ALFA survey (ALFALFA, Giovanelli et al. 2005) surveys.The data collected by these projects have been crucial for identifying spiral galaxies with an enhanced radio continuum emission in the cluster environment with respect to the field (Gavazzi & Boselli 1999) and for deriving the radial decrease of the H i gas content of galaxies as a function of the distance from the cluster centre (Gavazzi et al. 2013).This last measurement is critical for comparison with the predictions of hydrodynamic cosmological simulations.The untargeted nature of the ALFALFA survey also allowed us to discover H i in extended tails of gas associated with perturbed galaxies (Haynes et al. 2007) or free floating within the cluster (Haynes 2007;Kent et al. 2007Kent et al. , 2009;;Kent 2010).While being sensitive to diffuse gas emission, the ALFALFA data are still relatively shallow and have a low angular resolution (3.2 ).They are thus not optimal for resolved studies.Highquality interferometric data are available only for ∼50 galaxies in the H i 21 cm (VIVA, Chung et al. 2009) and 12 CO(2-1), 13 CO(2-1), and C 18 O(2-1) lines (VERTICO, Brown et al. 2021).For this reason we are undertaking an untargeted survey of the Virgo cluster in the radio domain using different groundbased facilities: MeerKAT in the continuum at 1.4 GHz and in the H i line at 21 cm, and Low Frequency Array (LOFAR) A92, page 2 of 24 1.6 × 10 9 M 2 PA (b)  19 • ± 2 • TW incl. (b)  32 TW log f (Hα + [NII]) (d)  −12.22 ± 0.02 erg s −1 cm −2 2 SFR (a)  0.08 M yr −1 2 Notes.All quantities are scaled to a distance of 16.5 Mpc, the mean distance of the cluster, for a fair comparison with other VESTIGE works.The uncertainty on the distance can be assumed to be ∼1 Mpc, the virial radius of the Virgo cluster.at 144 MHz (Edler et al. 2023).These surveys will allow us to study with unprecedented sensitivity and angular resolution the radio emission properties of cluster galaxies, adding thus a new major milestone in the study of galaxy evolution in rich environments.
Our radio survey of the Virgo cluster (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN: ViCTORIA) will be presented in a future communication (de Gasperin et al., in prep.).Here we present the first results obtained from a pilot observation of ∼10 deg 2 of Virgo with MeerKAT.With only 42 min of integration per pointing, the excellent quality of the MeerKAT data allowed us to discover an extended tail of H i gas associated with the star-forming galaxy NGC 4523 (see Table 1), which is undergoing ram-pressure stripping as discussed in this publication.The paper is structured as follows: in Sect. 2 we present the MeerKAT observations and the data reduction, in Sect. 3 the multifrequency data gathered during other surveys, including unpublished narrowband Hα (VESTIGE) and far-ultraviolet (FUV; ASTROSAT UltraViolet Imaging Telescope (UVIT)) imaging.In Sect. 4 we compare the H i data to those available in the literature, and in Sect. 5 we combine them with those available at other frequencies for a multifrequency analysis.Discussion and conclusions are given in Sects.6 and 7, respectively.Throughout this work we assume the galaxy to be at the distance of the cluster (16.5 Mpc, Gavazzi et al. 1999;Mei et al. 2007;Blakeslee et al. 2009;Cantiello et al. 2018).

MeerKAT observations and data reduction
ViCTORIA pilot observations of the Virgo cluster were carried out using MeerKAT (Camilo et al. 2018;Mauch et al. 2020) on September 10 and 12, 2022.These observations were programmed to test the feasibility of a large untargeted radio survey at 20 cm covering the a large fraction of the Virgo cluster, planned to be carried out during 2023 (∼60 deg 2 , with a total of 125 h allocated; PI F. de Gasperin).The observations are part of the ViCTORIA project.The project includes data from MeerKAT and from the low-and high-band systems of LOFAR (Edler et al. 2023), and it will greatly increase the sensitivity and resolution of the wide-field coverage of Virgo in radio continuum between 42-1700 MHz and in the 21 cm line.Polarisation data from MeerKAT observations will also be obtained.
The present pilot observations concern ten partially overlapping fields covering an L-shaped pattern centred on M87 and stretching towards both the north and east.The data were taken with the MeerKAT 32 k mode, covering the frequency range 856-1712 MHz with 32 768 channels of width 26.123 kHz.We splitted the ten fields in two groups of five.Each group was observed in a single ∼4.5 hour-long MeerKAT observing session, during which we cycled nine times through the five fields integrating for 280 s on each field at each visit.We observed a phase calibrator for 2 min in between fivefield cycles (i.e.every ∼23 min) and observed the bandpass and polarisation calibrators for 10 min each at the end of the observing session.The total integration time was 42 min for each of our Virgo fields, which thanks to the excellent MeerKAT uv coverage and to our observing strategy ensured that all the relevant angular scales were fully covered for each field.
For the purpose of this paper, we reduced the MeerKAT pilot data using the Containerized Automated Radio Astronomy Calibration (CARACal) pipeline (Józsa et al. 2020).We adopted an identical data reduction strategy as in Serra et al. (2023), and refer to that paper for a detailed description.Here we summarise the main data reduction steps and highlight any differences compared to that paper.Given that our focus is on H i science, we only processed a subset of the MeerKAT data (HH and VV correlations only), focusing on the frequency range 1395.48-1437.22MHz with 2× binning in frequency.This resulted in 800 channels with a width of 52.246 kHz, which corresponds to ∼11.0 km s −1 for H i at redshift Following automated flagging of radio frequency interference (RFI) we performed standard cross-calibration.This included solving for antenna-based, time-independent delays and bandpass, and for antenna-based, frequency-independent, time-dependent complex gains (one solution every 23 min).We applied these calibration terms to the target fields, and flagged RFI in the resulting visibilities.We then performed standard frequency-independent, phase-only self-calibration on each field independently, adopting a 2 min solution interval.The imagingcalibration iterations are automated, with a continuum sourcefinding strategy fine-tuned to obtain a model of M87 of sufficient quality for our goal of obtaining good H i cubes.Science-quality radio continuum images were obtained with a separate data reduction using a wider frequency band, and will be described in a separate paper (Edler et al., in prep.).For each field, we subtracted the radio continuum emission in two steps: i) we Fourier transformed and subtracted the continuum clean components from the self-calibrated visibilities; and ii) we fitted 3 rd -order polynomials to the visibility spectra, and subtracted the result from the data.The latter step was performed together with Doppler correction to a common barycentric frequency grid.We imaged the resulting continuum-subtracted visibilities with i) Brigg's robust = −0.5 and 6 uv tapering, and ii) Brigg's robust = 0 and 15 uv tapering.This resulted in H i cubes with a resolution of 12 × 21 (PA = 170 deg) and 27 × 39 (PA = 178 deg), respectively, and A92, page 3 of 24 typical noise level ∼0.65 mJy beam −1 per 11-km s −1 -wide channel for each field and at both resolutions -within ∼10% of the expected value.The 3 fields closest to M87 have higher noise.
Finally, we mosaiced the H i cubes.Due to the limited overlap between the fields included in this pilot observations, the noise in the mosaic did not improve significantly over that of the single fields.We used the Source Finding Application (SoFiA, Serra et al. 2015;Westmeier et al. 2021) to detect and parameterise H i sources in the mosaic.The uncertainty on the derived parameters is a combination of statistical uncertainties and a ∼10% uncertainty on the flux scale.
One of our H i detections was NGC 4523, the target of this work, located at a projected angular distance of 2.87 • from M87 in the north direction.The noise of the HI mosaic cubes at the position of this galaxy was 0.75 mJy beam −1 at both resolutions of 12 × 21 and 27 × 39 .This was slightly worse than the noise level mentioned above for the single fields because none of the observed fields was centred on this galaxy, combined with the limited overlap between adjacent pointings in our MeerKAT pilot data.The corresponding 3σ column density sensitivity was 1.6 × 10 20 cm −2 at high resolution and 3.9 × 10 19 cm −2 at low resolution assuming a line width of 25 km s −1 .During the following analysis we subsequently used both the low-and the highresolution data cubes, the former to study the extended low column density features such as the extended H i gas tail, the latter to study the gas kinematics in the inner regions minimising beam smearing effects.

Multifrequency data
The H i MeerKAT data presented in this work were combined with those available at other frequencies to study the perturbing mechanism and its effect on NGC 4523.These include narrowband Hα imaging, sensitive to the distribution of the young stellar population (age 10 Myr, Kennicutt 1998;Boselli et al. 2009) and to the presence of diffuse ionised gas possibly removed during the interaction of the galaxy with the surrounding environment (e.g., Gavazzi et al. 2001Gavazzi et al. , 2018;;Yoshida et al. 2002Yoshida et al. , 2012;;Yagi et al. 2007Yagi et al. , 2010Yagi et al. , 2017;;Sun et al. 2007Sun et al. , 2010;;Fossati et al. 2012Fossati et al. , 2018;;Boselli et al. 2016aBoselli et al. , 2018b;;Poggianti et al. 2017), UV data, sensitive to the distribution of young stars (age ∼ 100 Myr, Kennicutt 1998;Boselli et al. 2009), and optical imaging tracing the distribution of the bulk of the stellar population.

Optical imaging
High-quality optical images of NGC 4523 in the u, g, i, z filters are available thanks to NGVS, a broad-band imaging survey of the entire Virgo cluster carried out at the Canada France Hawaii Telescope (CFHT) with MegaCam (Ferrarese et al. 2012).These data have been gathered using a specific observing strategy and reduced with the Elixir LSB pipeline (Ferrarese et al. 2012) especially tailored to detect possible extended, low surface brightness structures (shells, tidal tails) possibly produced during the interaction of galaxies with their surrounding environment (Duc et al. 2015).Thanks to this optimised data analysis, NGVS is sensitive to extended features of surface brightness µ(g) 29 mag arcsec −2 .The optical data will be analysed in Sect.5.1.

VESTIGE narrow-band Hα imaging
High-quality narrow-band Hα imaging data are available thanks to VESTIGE, a deep untargeted survey of the Virgo cluster car-ried out with MegaCam at the CFHT (Boselli et al. 2018a).NGC 4523 was observed in the narrow-band Hα filter MP9603 (λ c = 6591 Å; ∆λ = 106 Å) with an integration time of 2 h, and 12 min in the broad-band r-filter to secure the subtraction of the stellar continuum (Boselli et al. 2019).The camera has a pixel scale of 0.187 arcsec pixel −1 .The data were reduced using Elixir LSB (Ferrarese et al. 2012), a specific procedure especially tailored to detect low surface brightness and extended features such as those produced during the interaction of galaxies with the surrounding environment.The data were reduced as described in Boselli et al. (2018a), securing an accurate astrometric and photometric calibration of the data (few % uncertainty).The survey has a sensitivity of f (Hα) 4 × 10 −17 erg s −1 cm −2 (5σ) for point sources and Σ(Hα) 2 × 10 −18 erg s −1 cm −2 arcsec −2 (1σ after smoothing the data to ∼3 resolution) for extended sources, and has been carried out under excellent seeing conditions (FWHM = 0.76 ± 0.07 ).It is thus perfectly suited to detect star-forming regions of luminosity L(Hα) 10 36.5 erg s −1 .These extremely low Hα luminosities 1 are comparable to those emitted by a single early-B star, among the young stars the one with the lowest mass and temperature able to ionise the gas (Sternberg et al. 2003).The Hα data will be analysed in Sect.5.3.

ASTROSAT/UVIT and GALEX UV imaging
We also compared the H i data to UV data obtained during the ASTROSAT/UVIT legacy survey of the Virgo cluster (PI. A. Boselli) and the GUViCS (A GALEX Ultraviolet Virgo Cluster Survey) survey (Boselli et al. 2011).ASTROSAT/UVIT (Agrawal 2006;Tandon et al. 2020) imaging data of NGC 4523 were obtained thanks to a pointed observation of the galaxy in the FUV filter BaF2 (λ c = 1541 Å; ∆λ = 380 Å), within the field of view of the instrument ( 28 ) with an angular resolution of 1.5 .The observations were carried out with an integration time of 3635 s, reaching a typical surface brightness of µ(FUV) 25.9 AB mag arcsec −2 .The data were reduced as described in Tandon et al. (2020) using a zero point of z p = 17.771 mag.The astrometry of the field was checked against the accurate NGVS imaging data (Ferrarese et al. 2012).The galaxy has been also observed during the GUViCS survey with GALEX in the near-ultraviolet (NUV) band (λ c = 2316 Å, ∆λ = 1060 Å), with an integration time of 1629 s, able to reach a surface brightness limit of 27 AB mag arcsec −2 .The angular resolution of the NUV band image is of 5 .The FUV UVIT and NUV GALEX images of the galaxy will be analysed in Sect.5.3.

Derived parameters
The MeerKAT data were used to extract the H i parameters of the galaxy: the integrated H i flux S HI = 23.4 ± 0.2 Jy km s −1 , the flux-weighted recessional velocity V HI = 271 ± 6 km s −1 (see Sect. 5.2), and the line width of the integrated spectrum measured at 50% of the peak WHI 50 = 121 ± 6 km s −1 .These numbers are in line with those derived in the literature (see Table 2).
The SHI flux derived in this work can be converted in a H i gas mass of M HI = 1.50 × 10 9 M using Eq. ( 9) of Haynes & Giovanelli (1984), and used to estimate a (5) Bottinelli et al. (1990).
H i-deficiency parameter HI-de f = 0.39 using the recent calibration of Cattorini et al. (2023).This calibration, which is based on the optical diameter, is optimal for the present purpose since it is the only one which extends to the dwarf systems regime.This H i-deficiency is slightly larger than the typical dispersion of the scaling relation used to define this parameter in galaxies living in low-density environments (σ ∼ 0.3, Cattorini et al. 2023).Of the measured H i content, 67% (9.8 × 10 8 M ) is located on the galaxy, 33% (4.8 × 10 8 M ) in the tail outside the stellar disc at the 24.5 mag arcsec −2 i-band isophote, as shown in Fig. 1.

Analysis
5.1.Gas distribution and origin of the extended H i tail Figure 1 shows the H i gas distribution (moment-0 map) of NGC 4523.The same gas distribution is compared to that of the stars in Fig. 2. At the angular resolution of 27 × 39 the data reach a 3σ column density of N(HI) = 3.9 × 10 19 cm −2 for a line width of 25 km s −1 .The H i is distributed asymmetrically, with an extended tail of low-column-density gas (N(HI) 2.5 × 10 20 cm −2 ) stretching out to 10 kpc (projected distance) from the edge of the stellar disc measured at the i-band 24.5 mag arcsec −2 in the north direction.We stress that this 10 kpc length is a lower limit to the extension of the tail because the galaxy has an important velocity component along the line-of-sight (LoS) (∼800 km s −1 ) with respect to the cluster (∼1040 km s −1 , Kashibadze et al. 2020).
The possible presence of similar features in deep optical images is crucial for the identification of the dominant perturbing mechanism.As shown in Figs. 1 and 2, the deep NGVS frames indicate that stars are also distributed asymmetrically, as often observed in galaxies of this morphological type (SAB(s)m).The i-band 24.5 mag arcsec −2 isophote has a circular shape, while it gets elliptical (b/a = 0.73, where a and b are the major and minor diameters) and elongated in the north-south direction (PA ∼ 15 • ) at the surface brightness limit of 26.5 mag arcsec −2 .The northern side is more elongated (147.6 ) than the southern one (120.8 ) when measured from the galaxy centre.This asymmetric stellar distribution at a low surface brightness limit in the same direction as the HI gas tail suggests that the galaxy might have been perturbed by a gravitational interaction.
There are several arguments suggesting that tidal interactions cannot be the main or only cause of the observed asymmetric distribution of the H i out to large radius, and that ram pressure must play an important role.First of all, the deep optical images do not show any evident structure such as streams, tidal tails, or shells -which might have formed during a gravitational perturbation (Duc et al. 2015) -down to a surface bright-ness limit of µ(g) 29 mag arcsec −2 across the full extent of the H i tail.Such features cannot be totally ruled out if their surface brightness is below the detection limit of the NGVS and VESTIGE data, as some simulations suggest (∼33 mag arcsec −2 Mancillas et al. 2019), but observations of a large sample of galaxies gathered with MegaCam at the CFHT rather indicate that tidal features are rare at a surface brightness limit below 27.5 mag arcsec −2 (Sola et al. 2022).Furthermore, the presence of the H i tail combined with the relatively large content of atomic gas on the galaxy disc and the lack of ionised gas in the tail, which requires 100 Myr to get mixed with the surrounding hot IGM and change of phase (e.g., Boselli et al. 2022), suggest that the perturbing process is still ongoing, as also suggested by our simulations.It is thus unlikely that the tidal streams, if formed, have had enough time to disappear (optical features are relatively long lived: 0.7-4 Gyr according to Mancillas et al. 2019).Finally, the galaxy is located at the northern periphery of the cluster and is falling into it for the first time.Here the number of massive galaxies which might have recently perturbed NGC 4523 is very limited: IC 800 (VCC 1532), a late-type system of comparable stellar mass, is located at a projected distance d proj ∼ 55 kpc to the north, but has a relative LoS velocity with respect to NGC 4523 of >2000 km s −1 .Assuming as relative distance between the two objects their projected distance (lower limit) we can measure the gravitational acceleration exerted by this perturber on the stellar disc of NGC 4523 and compare it with that keeping the matter linked to its gravitational potential well as indicated in Henriksen & Byrd (1996) and conclude that an efficient gravitational perturbation with this object is unlikely.This is also the case for the galaxies IC 3522 (VCC 1585) east of NGC 4523 (d proj ∼ 75 kpc, with a difference in the LoS velocity of only ∆(v) ∼ 400 km s −1 but with a low stellar mass of M star 8.5 × 10 7 M ), IC 797 (VCC 1393) west of NGC 4523 (d proj ∼ 135 kpc, ∆(v) ∼ 1800 km s −1 ), NGC 4540 to the north-east (d proj ∼ 135 kpc, ∆(v) ∼ 1000 km s −1 ), or the two massive galaxies to the south (NGC 4501, d proj ∼ 250 kpc, ∆(v) ∼ 2000 km s −1 , NGC 4548, d proj ∼ 230 kpc, ∆(v) ∼ 200 km s −1 ).We recall, however, that in a cluster environment where the velocity dispersion of galaxies is of the order of 1000 km s −1 the relative position of the perturbed object and of its perturber can drastically change (∼1 kpc Myr −1 ), making the identification of any possible perturber very challenging.A further possibility is that the tail is a remnant of gas accreted by NGC 4523 after a minor merging event with a gas-rich system.This picture, however, seems ruled out by the regular velocity field of the gaseous component (see Sect. 5.2) which would keep a perturbed shape if the merging event is recent with respect to the galaxy revolution time ( 500 Myr).We also remind that merging events in clusters such as Virgo are very unlikely given the high velocity dispersion of galaxies (Makino & Hut 1997;Boselli & Gavazzi 2006).
All these arguments suggest that the H i gas tail has formed during a hydrodynamic interaction with the surrounding hot ICM (namely, ram pressure stripping).The diffuse stellar emission associated to the tail (see Fig. 2) has a blue colour indicating that it is composed of stars of age 300 Myr and might thus have been formed in the stripped gas (see Sect. 7.3).
The projected orientation of the tail, approximately in the opposite direction of the cluster centre, suggests that the galaxy is at its first infall into the cluster.More in detail, the fact that the H i gas column density has a steep gradient at the southeast edge of the stellar disc (Figs. 1 and 2), opposite to the north-west maximum extension of the tail, suggests that the A92, page 5 of 24  infall vector has a south-east component on the plane of the sky.
Combined with the fact that the galaxy is blue-shifted relative to the cluster (its recessional velocity is v sys = 271 km s −1 while the mean velocity of the main body of Virgo, cluster A, is A92, page 6 of 24 to estimate the velocity of the galaxy with respect to the ICM necessary to strip atomic hydrogen with a typical column density of Σ HI 5 M pc −2 .This gives V ⊥ 1000 km s −1 , and suggests that the velocity of the galaxy on the plane of the sky is V sky ∼ 600 km s −1 .The shape and the orientation of the spiral arms indicate that the galaxy is rotating anticlockwise on the plane of the sky by assuming trailing arms.Given that the northern side is blueshifted and the southern side red-shifted, this indicates that the eastern side of the disc is the one closer to the observer (see Fig. 3).Therefore, as H i is pushed away from the disc towards north-west by ram pressure, it continues rotating anticlockwise until, on the eastern side of the galaxy, its motion is directed against the ram-pressure wind.This explains the orientation of the low column density H i tail to the north as well as the compressed H i contours on the north-eastern side of the tail.It should also lead to a piling up of H i on the northern side of the galaxy.In Sect.5.2 we see what consequences this has on the H i kinematics.

Gas kinematics
The integrated H i profile of NGC 4523 is fairly asymmetric (Fig. 4).It has significantly more flux below than above systemic, and it is steeper at low velocity near the peak of the emission (∼200 km s −1 ) and shallower at higher velocities (∼340 km s −1 ).The H i velocity field (Fig. 5, left panels) is also asymmetric due to the presence of the low column density tail to the north.We note that the northern side of the galaxy is blueshifted compared to systemic.Therefore, the piling-up of gas on the northern side of the galaxy discussed in the previous section explains the steep rise of the integrated spectrum at low velocities (Fig. 4).
The position-velocity diagram drawn along the kinematical major axis (PA = 19 • , Fig. 6) shows that the gas stripped from the disc and now located in the northern tail decelerates to velocities closer to the systemic velocity of the galaxy; that is to say, the H i LoS velocity increases towards systemic with increasing distance from the centre along the tail.This is consistent with a red-shifting ram-pressure wind caused by the galaxy's blue-shift relative to the ICM.Because of the orientation of the galaxy relative to the wind, the stripping process occurs mainly edge-on.The gas removed at the south-eastern edge of the disc remains gravitationally bound to the galaxy, while that at the south-western edge gets stripped, producing the extended tail along the LoS in the integrated profile.
To derive the main kinematical parameters of the galaxy and quantify the impact of the external perturbation on the H i kinematics, we fitted the high-resolution H i velocity field of Fig. 5 following the 2D fitting method described in Epinat et al. (2008).This method is based on the Levenberg-Marquardt nonlinear least-square algorithm and has been upgraded using the MocKinG software 2 .We fitted the H i velocity field, assuming an analytic axisymmetric velocity field purely based on circular motion, with a solid-body rise in the inner regions and a flat part in the outer regions.The parameters of this modelled velocity field (i.e., the radius where the velocity field changes from solidbody to flat, and the circular velocity of the flat part) are optimised during the fit in minimising the residual velocity field (the difference between the observed and the model velocity fields).
All other model parameters included in the fit (central position and velocity, inclination, position angle) are assumed to be independent of radius.Given that we fit the velocity field by weighting pixels based on their H i flux intensity, the best-fitting parameter values are dominated by the inner galaxy regions, where the effects of the environmental interaction under investigation are minimal.We assessed the robustness of the fit by weighting the pixels using different methods, and found that the output parameters are remarkably stable.
Having the observed velocity field and the parameter determined from the fit of the velocity field, we can derive an observed rotation curve for both sides of the galaxy.For this purpose, we defined an angular sector of inclusion of 67.5 • around the major axis (on the galaxy plane) in order to minimise possible contamination from radial motions3 .Within this sector, we calculated the mean velocity along circular arcs (in the galaxy plane) whose geometry is defined by our best-fitting model, correcting the velocity values with the cosine of the angular distance from the major axis on the galaxy plane.The left panel of Fig. 7 shows that the resulting rotation curve exhibits a systematic difference between the approaching and receding side in the inner regions (R 2 kpc).This discrepancy can be removed by shifting the kinematic centre of the rotating disc by 4.4 (350 pc) (see Fig. 7, right panel).Such a small shift of the dynamical centre is consistent with the idea that the external perturbation only marginally affects the innermost regions of the disc, where the gravitational potential well is deep.In order to reach the best matching between both sides of the solid body shape of the rotation curve, we shifted the coordinates of the centre to the position RA = 12:33:47.47and Dec = 15:09:53.8,at ∼14 arcsec to the north-west of the photometric centre (RA(J2000) = 12:33:47.95and Dec = 15:10:05.7;NGVS catalogue, Ferrarese, priv. comm.).We repeated the fit of the velocity field and obtained v sys = 271 km s −1 , incl.= 32.2• , PA = 19.4• , and v flat = 103 km s −1 .This is our final model.Its parameters are given in Table 1 and   An accurate inspection of the rotation curve derived using the above best-fitting model parameters (Fig. 7, right panel) clearly shows that the receding side of the rotating disc begins to be perturbed at R 2.5 kpc from the galaxy centre.Here the rotational velocity begins to be lower than expected for a solid body rotation, it is fairly flat in between 3.5 R/kpc 6, and then steeply increases up to 9.5 kpc, dropping again at larger radius.The approaching side follows a solid body rotation with a steep gradient for R 3 kpc, then has a flatter but fairly constant slope up to 13 kpc.A mild increase of the rotation curve up to these radii is not surprising in a low-mass galaxy such as NGC 4523.What is, however, unusual is the abrupt change of slope observed at R 3 kpc and the presence of two well defined components inside and outside this radius.These kind of asymmetries are expected in galaxies undergoing a ram pressure stripping event (Kronberger et al. 2008a).Finally, the approaching side rotational velocity decreases more rapidly than Keplerian for R 13 kpc.All these peculiarities in the rotation curve might be related to the gas compression on the receding side and gas removal from the galaxy disc on the approaching side.
The difference between our data and the best model can be associated to the external perturbation.We thus subtracted the best-fitting model velocity field from the observed velocity field (see Fig. 8).The figure shows that within a large circular inner A92, page 8 of 24 region approximately delimited by the i-band 23.5 mag arcsec −2 isophote, both sides of the disc show residual velocities close to zero km s −1 .This suggests that the inner velocity field is not strongly perturbed by the interaction with the Virgo environment.In a region delimited between the 25 and 23.5 mag arcsec −2 isophotes, the residual velocities are abruptly positive on both sides of the galaxy along the minor axis.This difference with respect to the model is stronger on the west side, where the tail is more extended.Outside the i-band 25 mag arcsec −2 isophote we observe the following.First, within a wide double cone of ∼45 • around the minor axis in the west direction, the residual velocities are again close to zero km s −1 .Because we are around the minor axis of the galaxy, this shows that radial motions (if any) are not significant.Second, within the complementary double cone of ∼45 • around the major axis, the residual velocities are negative.This is evident on the approaching side (north), where the H i tail is located, but it is observable in a narrow region on the receding side (south), too.Finally, in the region in between these two double cones the residual velocities are typically positive.
Figure 8 gives some interesting information on the H i kinematics on small physical scales, too.It shows that the peaks in the residual do not correlate with those observed in the velocity dispersion map shown in Fig. 5 (right panel).This implies that the H i kinematics is not dominated by the turbulence within the star forming regions, but by local motions due to the external pressure that the galaxy is suffering (large scale shearing).
A92, page 9 of 24 To summarise, the H i gas distribution, the velocity field, the integrated H i line profile, and the position-velocity diagrams can be explained if the galaxy is entering into the cluster from the back, north-west side and interacting with the surrounding ICM through a ram-pressure wind opposite to the infall and non perpendicular to the stellar disc (Fig. 3).The gas on the southernwest side of the disc, which has a LoS velocity of ∼340 km s −1 , is removed from the disc mainly along the LoS (and is thus red-shifted by ram pressure), producing the red-shifted shallow shoulder of the integrated H i spectrum and the compressed H i contours visible in Fig. 2. The gas in the northern side, which is rotating against the ICM, is more efficiently removed, producing the northern diffuse tail.The LoS velocity of H i in the tail is pushed back towards systemic, likely producing the steep shoulder of the integrated H i spectrum.In the outer regions this gas is still rotating but on a different plane than the one of the stars.

Gas vs star formation
Figure 9 compares the distribution of the ionised gas traced by the Balmer Hα line with that of the cold atomic hydrogen.Despite its sensitivity, the Hα narrow-band image does not show any diffuse ionised gas emission associated with the H i gas tail in the northern side of the disc.A few compact H ii regions, however, are detected at the north-east edge of the H i gas tail outside the stellar disc defined by the i-band surface brightness isophote of µ(i) = 24.5 mag arcsec −2 .They have been identified as those objects with an excess in emission in the narrow-band Hα image not associated to background galaxies in the broadband r-image.These regions are unresolved in the low resolution H i gas map, and only a few of them are detected in the shallow FUV ASTROSAT/UVIT image (Fig. 10).The GALEX NUV band image, which is sensitive to the emission of a some-

Simulations
We simulated the galaxy using a 3D N-body code with two components: a non-collisional component that simulates the stellar bulge/disc and the dark halo, and a collisional component that simulates the ISM.The non-collisional component consists of 81920 particles that simulate the galactic halo, bulge, and disc.The characteristics of the different galactic components are approximately adapted to the observed properties.We adopted a model where the ISM is simulated as a collisional component, that is, as discrete particles that possess a mass and a radius and can have inelastic collisions (sticky particles).The 20 000 particles of the collisional component represent gas cloud complexes that evolve in the gravitational potential of the galaxy.During the disc evolution, the particles can have inelastic collisions, the outcome of which (coalescence, mass exchange, or fragmentation) is simplified following Wiegel (1994).This results in an effective gas viscosity in the disc.The model galaxy is somewhat smaller and has a somewhat steeper rotation curve than NGC 4523.
As the galaxy moves through the ICM, its gas clouds are accelerated by ram pressure.Within the galaxy's inertial system, its clouds are exposed to a wind caused by the galaxy's motion through the ICM.The effect of ram pressure on the clouds is simulated by an additional force on the clouds in the wind direction.Only clouds that are not protected by other clouds against the wind are affected.Since the gas cannot develop instabilities, the influence of turbulence on the stripped gas is not included in the model.The mixing of the intracluster medium into the ISM is very crudely approximated by a finite penetration length A92, page 10 of 24 of the intracluster medium into the ISM, that is, up to this penetration length the clouds undergo an additional acceleration caused by ram pressure.A scheme for star formation was implemented where stars are formed during cloud collisions and then evolve as non-collisional particles (see Vollmer et al. 2012).These newly formed star particles carry their time of formation.Model Hα maps were produced with all star particles whose ages are smaller than 10 Myr.The UV emission of a star particle in the ASTROSAT/UVIT FUV band is modelled by the FUV flux from single stellar population models from STARBURST99 (Leitherer et al. 1999).The total FUV distribution is then the extinction-free distribution of the FUV emission of the newly created star particles.Since NGC 4523 is located far away from the cluster centre (∼ √ 2 × 0.83 Mpc ∼ 1.2 Mpc), we assumed that the galaxy's velocity and the encountered ICM density do not significantly vary during the last 360 Myr.With a velocity of ∼1000 km s −1 this corresponds to a distance of 380 kpc.We also tested slowly increasing temporal ram-pressure profiles.Since these simulations gave similar results, we only discuss the constant ram-pressure simulations.All gas clouds beyond the galaxy's effective radius are assumed to have a constant gas surface density of 10 M pc −2 .The surface density of the gas clouds increases towards the galaxy centre (e.g., Vollmer et al. 2008b).The surface density of atomic hydrogen was calculated via the gas density (Sect.5.1 of Vollmer et al. 2008b).For the different simulations we assumed a galaxy velocity of 1000 km s −1 and varied the ICM density n ICM and the angle between the ram pressure wind and the galaxy's disc plane α.The parameters of the different simulations are presented in Table 3.
The timesteps of interest for the α = 45 • simulations are 40-50 Myr later than those of the α = 10 • and α = 30 • simulations: in the case of a more face-on stripping (α = 45 • ) more time is needed to push the ISM towards the north and west.For the projection angles we used an inclination of i = 40 • and a position angle of PA = 20 • .The eastern half of the galactic disc is the near side.Moreover, the azimuthal viewing angle az (e.g., Vollmer et al. 2012) was chosen to best reproduce the observations.The galaxy's normalised velocity vectors are (−0.06,−0.09, −0.99) for the α = 45 All six models show the same characteristics, which are similar to those of NGC 4523 (see Fig. 11 for the low resolution view of model A; the high resolution view of the same model as well both low-and high resolution images of all other models are shown in Appendix A): (i) the east-west and north-west H i asymmetry, (ii) a relatively dense outer gas arm at the south-east edge of the H i disc pointing towards the north-east direction, (iii) an overdense gas arm to the north-east, (iv) an asymmetric velocity field in the north-south direction with increasing velocities to the south and an almost constant velocity to the north with a slight increase to large radii, and (v) an asymmetric velocity field in the east-west direction with higher velocity than expected by symmetry in the north-western quadrant.The observed characteristics, which are not reproduced by the models, are the A92, page 11 of 24 more extended to the east compared to the observed H i distribution and that of models A, E, and F. The western extent of the observed H i distribution is best reproduced by models A and F.
The velocity field along the galaxy's major axis is also somewhat different between the model: whereas the northern sides of models A and E better reproduce the observations, the southern side is slightly better reproduced by model D. The southern lowest radial velocities of models E and F are significantly lower than the corresponding observed velocities.In addition, in models E and F the eastern half of the zero-velocity contour is horizontal or slightly bent to the north, contrary to observations.Only models A and F show H ii regions to the north-east of the stellar disc as it is observed (lower left panel of Fig. 11).The south-eastern star-forming arm observed in Hα and UV is reproduced by models A, D, and F (lower right panels of Fig. 11, and the figures in the Appendix).
We conclude that models A, D, E, and F reproduce the observed characteristics of NGC 4523 in a satisfactory way.Model A (α = 30 • ) is somewhat preferred because it simultaneously reproduces the main characteristics of the observed H i distribution and kinematics and the observed distribution of H ii regions outside the optical disc.
If the galaxy's velocity within the ICM is 1000 km s −1 , the encountered ICM density is 1-2 × 10 −4 cm −3 .The model ICM density is consistent with the value derived from X-ray observations (Schindler et al. 1999).The galaxy moves to the south-east with a velocity of 400 ± 300 km s −1 .For the position of NGC 4523 with respect to M 87 we assumed (−188, 810, 834) kpc.Based on these numbers and the orbit modelling of Vollmer et al. (2001), we estimated NGC 4523's closest approach to the cluster centre (M 87) to be D ∼ 500-600 kpc in about 1 Gyr.The galaxy's velocity at closest approach will be ∼1300 km s −1 .The galaxy will thus undergo a maximum ram pressure that will be five to ten times higher than the current ram pressure.Since the current stripping radius is about 8 kpc, it is probable that NGC 4523 will lose most, if not all, of its ISM within the next ∼1.5 Gyr.

The MeerKAT survey of the Virgo cluster
Pilot observations of the Virgo cluster undertaken with the MeerKAT radio telescope are confirming, once again, the power of deep H i 21 cm observations in identifying galaxies undergoing a perturbation with their surrounding environment.With only 42 min of exposure per field, MeerKAT is able to detect H i tails of stripped gas at column densities as low as 3.9 × 10 19 cm −2 (0.3 M pc −2 ; 3σ assuming a 25 km s −1 line width) at a spectral resolution of 11 km s −1 and an angular resolution of ∼30 , sufficient to resolve most of the galaxies at the distance of the cluster (16.5 Mpc, Gavazzi et al. 1999;Mei et al. 2007).We thus expect that, once completed, the full survey will detect several star-forming systems and resolve most of them.Given the tight relation between the atomic gas content and the star formation activity of galaxies (e.g., Boselli et al. 2001), valid also within the Virgo cluster (e.g., Gavazzi et al. 2013;Boselli & Gavazzi 2014), we expect to detect ∼200 of the 384 galaxies detected in Hα by the VESTIGE survey (Boselli et al. 2023).Not all Hα emitting sources will be detected because of the limited coverage of the cluster (60 deg 2 for the proposed MeerKAT observations vs. 104 deg 2 for the VESTIGE survey) and sensitivity, able to detect only objects with H i gas masses M(HI) 2 × 10 6 M (3σ over a 100 km s −1 linewidth).Despite this limitation, the MeerKAT observations will increase by at least a factor of four the number of H i resolved galaxies observed by the VIVA survey in Virgo (Chung et al. 2009), extending the parameter space to dwarf irregular and early-type gas-rich systems.

Ram pressure at the periphery of the cluster
The extraordinary data gathered for NGC 4523 are a further confirmation that ram pressure is an efficient mechanism responsible for the stripping of the atomic gas content of galaxies of stellar mass M star 10 9 M entering a cluster of mass M 200 10 14 M (Boselli et al. 2022).The fact that the galaxy is located at 2.87 • from M87, the centre of the cluster, thus at R ∼ 0.85 × r 200 (see Table 1), suggests once again that in similar environments ram pressure stripping is efficient at least up to the virial radius of the cluster (Boselli et al. 2022).As mentioned in Sect.5, the orientation of the H i tail in the direction opposite to the cluster centre suggests that the galaxy is at its first infall into the cluster.Since the amount of H i already removed by the interaction with the ICM, if any, is still minor as deduced from the H i-deficiency parameter, this implies that ram pressure stripping is a rapid process, as already deduced by targeted observations of other representative objects (Vollmer et al. 2004;Boselli & Gavazzi 2006;Pappalardo et al. 2010;Abramson et al. 2011;Fossati et al. 2018), of large statistical samples (Crowl & Kenney 2008;Boselli et al. 2008aBoselli et al. ,b, 2016b;;Boselli & Gavazzi 2014;Vollmer 2009;Gavazzi et al. 2013;Vulcani et al. 2020), or of tuned hydrodynamic simulations (Vollmer et al. 2006(Vollmer et al. , 2008a(Vollmer et al. ,b, 2012(Vollmer et al. , 2018(Vollmer et al. , 2021;;Boselli et al. 2021).A recent stripping process during its first infall at the periphery of the cluster is also consistent with the lack of a ionised gas tail.Indeed, the stripped cold gas requires a few hundreds million years to mix with the hot ICM and change of phase, as indicated by the hydrodynamic simulations of IC 3476 (Boselli et al. 2021).
Simulations suggest that under specific conditions ram pressure stripping can trigger the star formation activity of the perturbed galaxies (Fujita & Nagashima 1999;Bekki & Couch 2003;Steinhauser et al. 2012;Bekki 2014;Henderson & Bekki 2016;Steyrleithner et al. 2020;Lee et al. 2020;Troncoso-Iribarren et al. 2020).Figure 12 shows the position of NGC 4523 on the main sequence relation defined by all the star-forming objects detected by the VESTIGE survey (Boselli et al. 2023).The total star formation rate of NGC 4523 and those of the other VESTIGE galaxies are derived assuming a Chabrier IMF after correcting the data for [NII] contamination and dust attenuation as described in Boselli et al. (2023).The galaxy is located on the main sequence relation of gas-rich systems (HI-def ≤ 0.4), here considered as unperturbed objects.In NGC 4523, as in most of the ram pressure stripped galaxies identified within the cluster, the overall star formation activity is not increased (Boselli et al. 2022(Boselli et al. , 2023)).Given that the ram pressure exerted on the gas disc is close to edge-on, and that this geometrical configuration is the most favourable for enhancing the activity of star formation (Kronberger et al. 2008b;Boselli et al. 2022), it is conceivable that the increased activity observed in the GASP sample (Vulcani et al. 2018) or in other ram pressure stripped galaxies in Coma and A1367 (Molnár et al. 2022;Boselli et al. 2022) might be due to selection effects (optically selected tails might favour galaxies undergoing a more violent or more advanced gas stripping phase than H i selected objects), or to the more hostile environment typical of these more massive clusters (M 200 10 15 M in Coma vs. M 200 10 14 M in Virgo).However, episodes of increased activity have been observed also A92, page 13 of 24 Fig. 12. Location of NGC 4523 on the main sequence relation (filled cyan square) defined by star-forming galaxies in the Virgo cluster (adapted from Boselli et al. 2023).Late-type galaxies (≥Sa) are indicated by black symbols, early-type galaxies by red symbols.Filled dots are for H i gas rich objects (HI-def ≤ 0.4), while H i deficient objects (HI-def > 0.4) by empty circles.Star formation rates have been derived assuming stationary conditions.The Y-axis on the right side gives the corresponding Hα luminosities corrected for dust attenuation.The black solid line shows the best fit obtained for gas-rich star-forming systems.The horizontal dotted and dashed lines indicate the corresponding SFR derived using the number of ionising photons produced by an O3 class III and a B0 class III star as derived using the model atmospheres of Sternberg et al. (2003).The typical uncertainty in the data is shown at the lower left corner.
in Virgo (e.g., IC 3476, Boselli et al. 2021) and in other local low-mass clusters (JO206, Ramatsoku et al. 2019).Their short duration and the moderate increase in the star formation activity makes a negligible effect on a complete star-forming population (Boselli et al. 2022(Boselli et al. , 2023)).If compared to the predictions of the 2D models of gas stripping in cluster galaxies presented in Boselli & Gavazzi (2006, 2014), Boselli et al. (2023), the position of NGC 4523 on the main sequence relation and its H i gas deficiency (HI-def = 0.39) suggest that the galaxy has recently started its interaction with the surrounding ICM ( 100 Myr, see Fig. 9 of Boselli et al. 2023; notice that the model plotted in this figure is for a galaxy with a rotational velocity of 70 km s −1 , similar to the one of NGC 4523).The star formation activity is marginally reduced probably because most of the stripped gas is located in the outer disc where star formation does not necessarily take place.This evolutionary picture is consistent with the simulations presented in Sect.6.

Star formation in the tail
Twelve star-forming regions have been detected in the continuum-subtracted Hα narrow-band image outside the i-band surface brightness isophote of µ(i) = 24.5 mag arcsec −2 .They are all within the H i gas tail of NGC 4523.These regions are located only on the north-east side of the tail (approaching side) where ram pressure stripping is expected to be the most efficient given the combined effect of rotation and the galaxy journey within the cluster, infalling from the north-east and back side.It is possible that these H ii regions have been formed within the turbulent structures at the interface of the two media (ISM and ICM), as indeed suggested by our simulations.We can age-date these H ii regions following the prescription given in Boselli et al. (2018b) by comparing their observed (Hα-FUV) and (FUV-NUV) age-sensitive colours to the typical synthetic colours of H ii regions as described in Boselli et al. (2018b).
We calculated the colours of these regions after extracting their fluxes as described in Fossati et al. (2018).As in Boselli et al. (2018b) we run the CIGALE code (Boquien et al. 2019) to derive star formation rates, stellar masses and ages of each individual region using the flux densities measured in the VESTIGE Hα, UVIT BaF2 FUV, GALEX NUV, CFHT ugriz, and Spitzer IRAC1 and IRAC2 bands, but derived these physical parameters only for a few of them given the very low detection rate in the optical and near-IR bands (see Table 4).
As expected, the typical colour of these regions suggests that they are recent since they are composed of stars with a typical age of 30 Myr, Fig. 13.These ages are comparable to those derived for the H ii regions formed within the tails of other perturbed galaxies in Virgo, such as the harassed gas of NGC 4254 (Boselli et al. 2018b) or in the ram pressure stripped gas of IC 3418 (Fumagalli et al. 2011).For comparison, if the galaxy is moving on the plane of the sky at ∼1000 km s −1 , a velocity comparable to the velocity dispersion of cluster A (the main body of the Virgo cluster, σ = 799 km s −1 , Boselli et al. 2014a) or to that of the late-type systems within the whole Virgo cluster (σ = 1150 km s −1 , Boselli & Gavazzi 2006), it would take 10 Myr to travel 10 kpc, the mean distance of these H ii regions from the outer stellar disc.Similarly, in 30 Myr the galaxy would make ∼10% of its revolution.
Worth noticing is also the fact that star formation within the tails of ram pressure stripped gas is not ubiquitous.The VESTIGE survey, which has covered the whole cluster with exquisite sensitivity and angular resolution, has identified tails without any star-forming region (e.g., NGC 4569, Boselli et al. 2016a), others with a few (NGC 4388, Yoshida et al. 2002;NGC 4330, Fossati et al. 2018;IC 3476, Boselli et al. 2021).There are also examples of galaxies with several H ii regions outside their stellar discs but without any apparent associated diffuse ionised gas tail (e.g., IC 3418, Hester et al. 2010;Fumagalli et al. 2011; A92, page 14 of 24 Fig. 13.Age-sensitive colour-colour diagrams.Upper left: Hα-FUV vs. FUV−NUV, lower left: NUV−u vs. FUV−NUV, lower right: NUV−u vs. Hα-FUV.Black filled dots are for regions detected in all bands, red filled triangles are lower limits in Hα-FUV, blue filled triangles upper limits in NUV − u, while the cyan filled square the mean value for the diffuse extended tail.The blue, green, and red filled dots and solid lines indicate the expected colours at different ages (in Myr) for the assumed star formation history for an E(B − V) = 0.0, 0.1, and 0.2, respectively.The green numbers are the ages in Myr for E(B − V) = 0.1.Jáchym et al. 2013;Kenney et al. 2014;Hota et al. 2021), at least at the sensitivity of the available observations.The process of star formation and the conditions under which this occurs in the tails of stripped material is still poorly known.Star formation requires gas cooling in the molecular phase to take place.The molecular gas is hardly removed during a ram pressure stripping event since it is principally located within giant molecular clouds with a limited cross section with respect to the external pressure (Boselli et al. 2022).Molecular gas in clumpy regions has been observed in the tails of several ram pressure stripped galaxies (e.g., Vollmer et al. 2005Vollmer et al. , 2008b;;Jáchym et al. 2014Jáchym et al. , 2017Jáchym et al. , 2019;;Moretti et al. 2018Moretti et al. , 2020b)).It is thus conceivable that the cold atomic gas located in the tails, under some still unclear conditions, collapses to form new stars.This might happen in the turbulent regions formed by the instabilities at the edges of the stripped cold ISM and the surrounding hot ICM.The mixing of the different gas phases can induce cooling of the ICM into dense clouds (Tonnesen & Bryan 2021), while magnetic fields can keep the cold gas confined within filaments, thus affecting the cooling process (Tonnesen & Stone 2014;Ruszkowski et al. 2014;Müller et al. 2021).Furthermore, it is also unknown which is the efficiency with which the molecular gas component is converted into stars within ram pressure stripped tails (e.g., Boissier et al. 2012;Villanueva et al. 2022).
The completion of the MeerKAT survey of the Virgo cluster will certainly boost our understanding of this important physical process.Despite the presence of the bright radio galaxy M87 (Virgo A) in the field, the quality of the data is excellent in terms of sensitivity and angular resolution, and will be increased by a factor of ∼2 once the survey will be completed thanks to an increased overlap of pointings.ViCTORIA will provide us with critical information on the atomic gas content, the main gas component stripped during the interaction since loosely bound to the gravitational potential well of the galaxy.Combined with that of the other gas phases (molecular, VERTICO survey, Brown et al. 2021;ionised, VESTIGE survey, Boselli et al. 2018a;hot, eROSITA, Predehl et al. 2021) H i data will be critical to understand under which conditions the stripped gas can cool and condense to form stars or heat up changing of phase and becoming first ionised, then hot gas.This complete set of multifrequency data on the different gas phases will thus be crucial for constraining tuned hydrodynamic simulations of interacting galaxies now able to resolve giant molecular clouds and star-forming regions within the tails.ViCTORIA will also provide us with a unique set of radio continuum data at different frequencies, where the emission is sensitive to the energy loss of relativistic electrons in weak magnetic fields, two other major ingredients in the study of galaxy evolution in rich environments.

Conclusion
During pilot observations of the Virgo cluster carried out with the MeerKAT radio telescope we discovered an extended (10 kpc projected distance) low column density atomic gas tail (N(HI) 2.5 × 10 20 cm −2 ) in a dwarf galaxy (NGC 4523, M star = 1.6 × 10 9 M ) at the northern periphery of the cluster (R 0.85 × r 200 ).The multifrequency analysis of the data consistently suggests that NGC 4523 is suffering a ram pressure stripping event which has already removed part of the atomic gas.The overall star formation activity of the galaxy, estimated using new Hα narrow-band imaging gathered during the VESTIGE survey, is not significantly reduced if compared to that of similar unperturbed objects.A few compact H ii regions, whose age has been estimated comparing the prediction of spectral energy distribution (SED) models to their observed Hα, ASTROSAT/UVIT FUV, GALEX NUV, and NGVS u colours ( 30 Myr), have been formed within the tail of stripped gas.NGC 4523 is a further example indicating that ram pressure stripping is efficient in intermediate mass clusters such as Virgo (M 200 10 14 M ) up to the virial radius, confirming once again the main role of this hydrodynamic process in shaping galaxy evolution in rich environments.
The results presented in this work are a further demonstration of the power of untargeted surveys of nearby clusters in identifying galaxies undergoing a perturbation.Once completed, the MeerKAT survey of the Virgo cluster will provide us with a unique view of the effects of the environment on the cold atomic gas of infalling galaxies.Combined with similar information gathered thanks to the VESTIGE survey for the ionised gas component, the H i data will be crucial to understand the effects induced by the perturbations on the process of star formation.They will be necessary also to understand the fate of the stripped gas, the formation of possible compact sources outside the stellar discs of perturbed galaxies, and ultimately the pollution of the ICM.

Fig. 1 .
Fig. 1.Low-resolution (27 × 39 ; left) and high-resolution (12 × 21 ; right) H i gas distribution of NGC 4523.Green contours are at column densities of N(HI) = 2 n × 3.9 × 10 19 cm −2 with n = 0, 1, . . . 5. The black contour shows the i-band isophote at 24.5 mag arcsec −2 .The open dot shows the kinematic centre of the galaxy, the ellipse indicates the size of the beam, the arrow the direction of the cluster centre (M87), located at 0.83 Mpc projected distance.

Fig. 2 .
Fig. 2. Pseudo-colour image of the galaxy NGC 4523 obtained combining the NGVS (Ferrarese et al. 2012) and VESTIGE (Boselli et al. 2018a) optical u and g in the blue channel, the r and the narrow-band Hα in the green, and the i and the continuum-subtracted Hα in the red.Red contours are shown at column densities of N(HI) = 2 n × 3.9 × 10 19 cm −2 with n = 0, 1, . . . 5. The cyan contours show the i-band isophote at 24.5 mag arcsec −2 .The yellow arrow indicates the position of the cluster centre (M87), located at 0.83 Mpc projected distance.

Fig. 3 .
Fig.3.Cartoon of the orbit of the galaxy and its orientation within the cluster.The yellow ellipse represents the stellar disc, the red contour the H i distribution, the black shaded region the trajectory of the galaxy during its first infall into the cluster.The position of the cluster centre (M87) and the direction of the wind are indicated by the black arrows.α is the angle between the wind direction and the plane of the stellar disc.

Fig. 4 .
Fig. 4. Integrated H i line profile of NGC 4523.The red arrow indicates the systemic heliocentric velocity (271 km s −1 ).

Fig. 5 .
Fig. 5. Low-resolution (27 × 39 ; upper panels) and high-resolution (12 × 21 ; lower panels) H i gas velocity field (moment-1 , left) and velocity dispersion (moment 2, right) of NGC 4523.The black contour shows the stellar distribution as traced by the i-band 24.5 mag arcsec −2 isophote.The open dot show the kinematic centre of the galaxy, the black ellipse indicates the beam size and the arrow shows the direction to the cluster centre (M87).The thick white contour represents the systemic velocity of 271 km s −1 , and the thinner white contours are drawn at intervals of 10 km s −1 .

Fig. 8 .
Fig. 8. High-resolution residual velocity field considering a twoslope rotation curve.The black contours show the i-band 23.5 and 25.0 mag arcsec −2 isophotes, the black open dot the kinematic centre.

Fig. 9 .
Fig. 9. H i gas distribution (green contours) overlayed on the continuum-subtracted Hα image of NGC 4523 derived using the VESTIGE narrow-band imaging data.Green contours are shown at column densities of N(HI) = 2 n × 3.9 × 10 19 cm −2 with n = 0, 1, . . . 5. The colour scale on the right gives the Hα surface brightness in units of erg s −1 cm −2 arcsec −2 .The H ii regions located outside the stellar disc and analysed in Sect.6 are identified with yellow circles.The yellow arrow indicates the direction of the cluster centre (M87), located at a projected distance of 0.83 Mpc.
a) Galaxy velocity within the ICM.(b)ICM density.(c)Angle between the ram pressure wind and the disc plane.(d)Time since the beginning of the external pressure.(e)Azimuthal viewing angle.stripping, (−0.31, −0.11, −0.95) for the α = 30 • stripping, and (−0.44, −0.52, −0.73) for the α = 10 • stripping (the parameters of the orbit of model A are given in the Appendix).Therefore, the galaxy's dominant velocity component is the radial (z) component in all simulations.For α = 45 • and α = 10 • the galaxy moves approximately to the south-east, whereas it moves more to the east for α = 30 • .The associated velocities on the plane of the sky are ∼100, ∼330, and ∼680 km s −1 for α = 45 • , 30 • , and 10 • , respectively.

Fig. 11 .
Fig. 11.Output of the simulations for model A. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n = 0, 1, 2, . ..) Mpc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9 × 10 19 × 10 0.5×n (n = 0, 1, 2 . ..) cm −2 .
acknowledge financial support from 'Programme National de Cosmologie and Galaxies' (PNCG) funded by CNRS/INSU-IN2P3-INP, CEA and CNES, France, and from "Projet International de Coopération Scientifique" (PICS) with Canada funded by the CNRS, France.This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration and of the GOLDMine database (http: //goldmine.mib.infn.it/)(Gavazzi et al. 2003).Part of the data published here have been reduced using the CARACal pipeline, partially supported by ERC Starting grant number 679627 FORNAX, MAECI Grant Number ZA18GR02, DST-NRF Grant Number 113121 as part of the ISARP Joint Research Scheme, and BMBF project 05A17PC2 for D-MeerKAT.Information about CARACal can be obtained online under the URL: https://caracal.readthedocs.io.M.B. gratefully acknowledges support by the ANID BASAL project FB210003 and from the FONDECYT regular grant 1211000.N.Z. is supported through the South African Research Chairs Initiative of the Department of Science and Technology and National Research Foundation.

Fig. A. 2 .
Fig. A.2. Output of the simulations for model B. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n=0,1,2,..) M pc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9×10 19 × 10 0.5×n (n=0,1,2...) cm −2 .

Fig. A. 3 .
Fig. A.3.Output of the simulations for model C. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n=0,1,2,..) M pc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9×10 19 × 10 0.5×n (n=0,1,2...) cm −2 .

Fig. A. 4 .
Fig. A.4.Output of the simulations for model D. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n=0,1,2,..) M pc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9×10 19 × 10 0.5×n (n=0,1,2...) cm −2 .

Fig. A. 5 .
Fig. A.5.Output of the simulations for model E. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n=0,1,2,..) M pc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9×10 19 × 10 0.5×n (n=0,1,2...) cm −2 .

Fig. A. 6 .
Fig. A.6.Output of the simulations for model F. Upper panels, left: Model H i moment 0 map, the contour levels show the model stellar distribution at 10 × 10 0.5×n (n=0,1,2,..) M pc −2 .Right: Model H i moment 1 map, the white contour levels are from −70 to 70 km s −1 in steps of 10 km s −1 , the black contour the model stellar distribution at 10 M pc −2 .The empty dot shows the kinematic centre.Lower left panel: Model H i contours on model Hα distribution.Lower right panel: Model H i contours on model FUV distribution.The H i contour levels on the Hα and the FUV images are N(HI) = 3.9×10 19 × 10 0.5×n (n=0,1,2...) cm −2 .

Table 2 .
HI measurements of NGC 4523 in the literature.

Table 4 .
Boselli et al. 202310) regions outside the stellar disc.=0.02 mag) as inBoselli et al. (2018b).Star formation rates are derived from the Hα luminosity using the calibration ofCalzetti et al. (2010).Given the low star formation regime, these numbers should be taken as indicative (seeBoselli et al. 2023for details).Both SFR and M star have been derived assuming a Chabrier IMF.