Journal of the Korean Astronomical Society https://doi.org/10.5303/JKAS.2018.00.0.1 00: 1 ∼ 99, 2018 April pISSN: 1225-4614 · eISSN: 2288-890X Published under Creative Commons license CC BY-SA 4.0 http://jkas.kas.org

STAR FORMATION ACTIVITYOF GALAXIES UNDERGOING RAM PRESSURE STRIPPINGINTHE VIRGO CLUSTER Jae Yeon Mun1, Ho Seong Hwang2, Myung Gyoon Lee1, Aeree Chung3, Hyein Yoon3,4,5, and Jong Chul Lee2 1Department of Physics and Astronomy, Seoul National University, Gwanak-gu, Seoul 151-742, Korea; [email protected] 2Korea Astronomy and Space Science Institute, Daejeon 305-348, Korea; [email protected] 3Department of Astronomy, Yonsei University, 50 Yonsei-ro, Seodaemun-gu, Seoul 03722, Korea 4Sydney Institute for Astronomy, School of Physics, A28, The University of Sydney, Sydney, NSW 2006, Australia 5ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D) Received XXX; accepted YYY

Abstract: We study galaxies undergoing ram pressure stripping in the Virgo cluster to examine whether we can identify any discernible trend in their star formation activity. We first use 48 galaxies undergoing different stages of stripping based on HI morphology, HI deficiency, and relative extent to the stellar disk, from the VIVA survey. We then employ a new scheme for galaxy classification which combines HI mass fractions and locations in projected phase space, resulting in a new sample of 365 galaxies. We utilize a variety of star formation tracers, which include g − r, WISE [3.4] − [12] colors, and starburstiness that are defined by stellar mass and star formation rates to compare the star formation activity of galaxies at different stripping stages. We find no clear evidence for enhancement in the integrated star formation activity of galaxies undergoing early to active stripping. We are instead able to capture the overall quenching of star formation activity with increasing degree of ram pressure stripping, in agreement with previous studies. Our results suggest that if there is any ram pressure stripping induced enhancement, it is at best locally modest, and galaxies undergoing enhancement make up a small fraction of the total sample. Our results also indicate that it is possible to trace galaxies at different stages of stripping with the combination of HI gas content and location in projected phase space, which can be extended to other galaxy clusters that lack high-resolution HI imaging. Key words: galaxies: clusters: general − galaxies: clusters: individual (Virgo) − galaxies: evolution − galaxies: star formation

1.I NTRODUCTION & Kenney 2008; Jaff´eet al. 2015). However, there have been observational and theoretical studies insin- Star formation and nuclear activity of galaxies are uating the possibility of ram pressure induced star for- strongly affected by the cluster environment (von der mation enhancement prior to quenching (e.g., Bothun Linden et al. 2010; Hwang et al. 2012a) through various & Dressler 1986; Fujita & Nagashima 1999; Kronberger physical mechanisms: ram pressure stripping (Gunn & et al. 2008; Kapferer et al. 2009). One-sided tails of Hα Gott 1972), thermal evaporation (Cowie & Songaila emission or young stellar knots extending beyond disks 1977), strangulation (Larson et al. 1980; Bekki et al. of galaxies in distant clusters have shown to be linked 2002), galaxy harassment (Moore et al. 1996), and cu- to interactions with the ICM (e.g., Cortese et al. 2007; mulative galaxy-galaxy hydrodynamic/gravitational in- Hester et al. 2010; Smith et al. 2010), giving further teractions (Park & Hwang 2009). Among them, ram support to the scenario of ram pressure induced star pressure stripping is one of the most effective ways formation enhancement. While not all galaxies subject to remove gas from galaxies. Ram pressure arises to ram pressure stripping show such features, those that arXiv:2101.07472v2 [astro-ph.GA] 24 Mar 2021 due to the hot X-ray gas embedded in the intraclus- exhibit extended tails of optically bright gas are known ter medium (ICM), and pushes against the interstel- as “jellyfish” galaxies. As of recent, many observational lar medium (ISM) of a galaxy as the galaxy travels studies of jellyfish galaxies have captured the enhance- through the galaxy cluster. As a galaxy approaches ment of star formation activity in their disks and tails the cluster core, the strength of ram pressure increases, (Vulcani et al. 2018; Ramatsoku et al. 2019). In par- and a greater amount of its gas - the fuel necessary for ticular, Vulcani et al. (2018) categorized their sample star formation - is stripped (Gunn & Gott 1972). In of jellyfish galaxies from the GAs Stripping Phenom- general, ram pressure stripping (RPS) is expected to ena survey (GASP; Poggianti et al. 2017) into different quench star formation in cluster galaxies (e.g., Koop- stripping stages - initial, moderate, and extreme strip- mann & Kenney 2004b; Boselli & Gavazzi 2006; Crowl ping, along with truncated disks - based on their Hα Corresponding author: H. S. Hwang morphology. Upon comparing the star formation ac- 1 2 Mun et al. tivity of their sample to field galaxies of similar stel- of HI-detections within the “stripped” and “virialized” lar mass, they show that jellyfish galaxies undergoing regions. Moreover, they investigated the stellar pop- moderate to extreme stripping show the strongest en- ulations of their galaxy sample using NUV − r color hancement of star formation. As such, it is of great to find that the bluest galaxies were generally found interest to probe the change in star formation activ- outside the “stripped” and “virialized” zones, while the ity of a galaxy as it undergoes different stages of ram redder galaxies seemed to dominate those two regions. pressure stripping. While they were not able to take HI morphology into Given its close proximity, the Virgo cluster has of- account for their analysis, their results confirm the gen- ten been the target of observational studies of interstel- eral trend that as blue, star-forming galaxies fall into lar medium (ISM) - intracluster medium (ICM) inter- a cluster, its gas content will be gradually stripped by actions (e.g., Boselli et al. 2018). Koopmann & Kenney ram pressure as it settles, with its star formation being (2004b) studied the Hα morphology of Virgo cluster and quenched in the process. isolated spiral galaxies to not only divide cluster galax- Given that observational studies of jellyfish galax- ies into different categories - normal, anemic, enhanced, ies have shown that a stage of star formation enhance- and truncated - but also to identify the history of envi- ment could occur prior to its gas content being com- ronmental processes that galaxies may have experienced pletely quenched, this then brings up the question of in the cluster environment. While Hα is a useful probe whether a similar analysis can be done for the Virgo of the current star formation activity, the distribution of cluster. Yoon et al. (2017) provides an ideal sample to neutral hydrogen (HI) of a galaxy offers a clearer picture answer such a question, as they have used the VIVA of the environmental processes a galaxy has undergone. (Chung et al. 2009) sample to categorize galaxies in the Given that HI gas is usually extended out towards the Virgo cluster into different stages of stripping based on outer disk region of a galaxy, it is often the first com- their HI morphology, HI deficiency (defHI, a measure of ponent to show perturbations due to external processes how deficient in HI a galaxy is relative to its field coun- such as tidal interactions and ram pressure stripping. terpart of the same size, Haynes & Giovanelli 1984), Chung et al. (2009) obtained high-resolution HI data and relative extent to the optical disk. The availabil- of a sample of 53 late type galaxies in the Virgo clus- ity of high-resolution HI imaging plays a crucial role ter, among which 46 of them were selected based on the in their study, as the peculiarities in the HI morphol- classification in Koopmann & Kenney (2004b) to probe ogy allows a direct probe of the environmental processes galaxies with a wide range of star formation properties. that the galaxy may have been subject to. Yoon et al. Known as the VLA Imaging survey of Virgo galaxies (2017) also briefly examined the star formation proper- in Atomic gas (VIVA), Chung et al. (2009) identified a ties of the categorized galaxies using the Hα equivalent trend in HI morphology with respect to the projected lo- width and Hα concentration from Koopmann & Ken- cation of galaxies within the cluster. Galaxies near or in ney (2004a). Their results (see Yoon et al. 2017, Fig. the core were often found with symmetric HI disks trun- 8) showed reduced star formation and increased con- cated with respect to their stellar disks, while galaxies centration of star formation in the galaxy center, along in the outskirts were generally HI-rich in comparison with decreasing HI gas content from early through post (but still HI poor compared to field galaxies). On the stripping stages. other hand, galaxies located at intermediate distances In this work, we aim to build upon the work of from the cluster core (i.e., M87; 0.6 . dM87 . 1 Mpc) Yoon et al. (2017) by examining the star formation were observed to have one-sided HI tails pointing away properties of their sample of galaxies in more detail, from the core, which led Chung et al. (2007) to suggest using various types of star formation tracers. To make that such galaxies have only recently began falling into up for the small sample size in their study, we expand the cluster, and that the impact of ram pressure on the the sample by using the Extended Virgo Cluster Cat- gas content of galaxies extends out towards intermedi- alog (EVCC; Kim et al. 2014). Since high-resolution ate distances. HI data are not available for the entirety of the EVCC, It is then possible to infer from the results of pre- we make use of HI mass fraction (relative HI mass to vious studies that the combination of HI morphology stellar mass, MHI/M*) and location in projected phase and projected location from the cluster center can offer space. While MHI/M* can work as a first-order indi- a closer look into the orbital history of a galaxy within cator of gas stripping, the location of galaxies in phase a galaxy cluster. As of recent, many studies have incor- space assures that the origin of stripping is likely to be porated the use of projected clustercentric distance and of ram pressure (Rhee et al. 2017; Jaff´eet al. 2015; Yoon line-of-sight velocity with respect to the cluster mean et al. 2017). in the form of a phase-space diagram to identify galax- This paper is organized as follows. In Section 2, ies at different stages in their orbits (Oman et al. 2013; we describe the sample and explain the derivation of Oman & Hudson 2016; Rhee et al. 2017). For example, physical properties used to examine the star formation Jaff´eet al. (2015) studied Abell 963, a galaxy cluster at activity of galaxies in the Virgo cluster. In Section 3, z = 0.2, by combining the availability of HI detections we explain how the different stripping classes were de- and location in projected phase space. They found that fined in both Yoon et al. (2017) and this study, follow- HI-detected galaxies were generally located in the “re- ing up with the analysis of the star formation activity cent infalls” region, whereas there was a noticeable lack per class. In Section 4, we discuss the limitations of SF Activity of Galaxies Undergoing RPS in Virgo 3 integrated properties in capturing local star formation (2003) IMF. The synthetic SEDs are computed at dif- enhancements in Virgo galaxies, and how the simplifi- ferent redshifts with δz = 0.001 as the furthest galaxy in cation of the classification scheme using HI gas fraction our sample is located at z ∼ 0.01. The BC03 templates and location in phase space allowed the trend of quench- provided by Le Phare have three metallicities and seven ing of star formation activity to be observed with a exponentially decreasing star formation models (SFR ∝ larger sample. We summarize our results in Section 5. e−t/τ ) with τ = 0.1, 0.3, 1, 2, 3, 5, 10, 15, and 30 Gyr. Throughout this paper, we adopt a Salpeter (Salpeter We allow the age of the stellar population to vary be- 1955) initial mass function (IMF) and assume a Λ cold tween 0 and 13 Gyr. Dust extinction was applied to dark matter cosmology with ΩM = 0.3, ΩΛ = 0.7, and the SPS templates using the extinction law of Calzetti -1 -1 H0 = 70 km s Mpc . We also adopt a distance of 16.5 et al. (2000), with values of E(B −V ) varying from 0 to Mpc, or (m − M) = 31.1, to the Virgo cluster (Mei et 0.6. Adopting different values of the stellar population al. 2007). All photometric magnitudes used are given parameters gives rise to a distribution of different esti- in the AB magnitude system. mates for the stellar mass of each galaxy. We adopt the median of the distribution as the stellar mass of each 2.D ATA galaxy. 2.1. EVCC We compare our stellar mass estimates with those The EVCC is comprised of fundamental information given by the SDSS MPA/JHU DR7 VAGC (Kauffmann such as cluster membership, morphology, radial veloci- et al. 2003), where stellar masses were measured based ties, and ugriz photometry of a total of 1589 galaxies on a SED fit of SDSS ugriz photometry with the SPS spread throughout the main body to the outskirts of models of Bruzual & Charlot (2003). The stellar mass the Virgo cluster. More precisely, the spatial coverage estimates from the MPA/JHU DR7 VAGC are based on of the EVCC extends out to ∼4 times the virial radius the Kroupa IMF. For consistency throughout this study, (r200 = 1.55 Mpc, McLaughlin 1999; Ferrarese et al. stellar mass estimates computed via Le Phare and those 2012) of the cluster, for which Kim et al. (2014) ex- given by the MPA/JHU DR7 VAGC are converted to cluded background galaxies with radial velocities larger a Salpeter (1955) IMF by multiplying the former by than 3000 km s-1. a factor of 100.25 (Chabrier 2003) and the latter by a With regards to the ugriz photometry, the Sloan factor of 100.2 (Kroupa 2001). Both conversion factors Digital Sky Survey (SDSS) photometric pipeline is were taken from Bernardi et al. (2010). We show the known to be unreliable for large and bright galaxies due results of the comparison in the left panel of Figure 1, to issues with deblending and sky subtraction (Abaza- which shows good agreement between the two stellar jian et al. 2004, 2009; West et al. 2010). To account for mass estimates with a median offset of about 0.08 dex. such shortcomings, Kim et al. (2014) performed their While there is in general a good agreement between the own Source Extractor (SExtractor, Bertin & Arnouts two methods, the rms of the difference between our esti- 1996) photometry on ugriz images from SDSS Data mates and those given by the MPA/JHU DR7 VAGC is Release 7 (DR7). As shown by Figure 13 in Kim et about 0.23 dex when including all the galaxies in com- al. (2014), the growing discrepancy between the pho- mon. The absolute uncertainty arising from panchro- tometric magnitudes with increasing angular size of a matic SEDs is approximately 0.3 dex (Conroy 2013), galaxy further confirms that the total flux derived from which suggests that our stellar mass estimates are over- the SDSS pipeline does not give an accurate represen- all consistent with the measurements given by the SDSS tation of the total light from nearby galaxies such as MPA/JHU DR7 VAGC. those in the Virgo cluster. We thus attest that physi- cal parameters derived using EVCC ugriz photometry, 2.2.2. Star Formation Rates & Mid-Infrared Colors such as stellar mass, are reliable measurements of the We use mid-infrared photometric data from the Wide- stellar population of the galaxies in the Virgo cluster. field Infrared Survey Explorer (WISE, Wright et al. 2010) to derive SFRs and examine mid-infrared col- 2.2. Physical Parameters of Galaxies ors. The WISE survey covers the entire sky in the mid- To examine the star formation properties of Virgo clus- infrared, offering four-band photometric data (3.4, 4.6, ter galaxies, we proceed by building a master catalog 12, and 22µm) for over 747 million objects. We used based on the EVCC. We add stellar masses, star forma- a matching tolerance of 600 (∼FWHM of the PSF at tion rates (SFRs), mid-infrared photometric data, along 22µm) between the EVCC and AllWISE Source Cata- with neutral hydrogen masses to the catalog. Here, we log1. This matching tolerance is small enough to ensure briefly describe each of the physical parameters. that we avoid matching our sample with the wrong mid- 2.2.1. Stellar Masses infrared counterparts. In the cross-matching process, we found some issues with the coordinates given by the The stellar masses of EVCC galaxies are computed us- EVCC. EVCC provides two different coordinates, one ing Le Phare (Arnouts et al. 1999; Ilbert et al. 2006). derived from SExtractor, another given by the SDSS We used EVCC ugriz photometry and redshifts to build fiber locations. The coordinates given by SExtractor synthetic spectral energy distributions (SEDs) from the were offset from the center of nearby galaxies, result- stellar population synthesis (SPS) models of Bruzual & Charlot (BC03; 2003), which are based on the Chabrier 1http://wise2.ipac.caltech.edu/docs/release/allwise/ 4 Mun et al.

M*, LePhare vs. M*, SDSS SFRIR vs. SFRSDSS

101 1011

e 0 r 10 10

a 10 R h I P R e L F , S * 9

M 10 10 1

108 10 2 S S S

D 2 1 0 1 S

S 10 10 10 10 D , 2 2 * 10 S 10 R M

F /

0 S 0

e 10 10 r /

a R h I P 2 2 e 10 R 10 L F , * S

M 108 109 1010 1011 10 2 10 1 100 101 M*, SDSS SFRSDSS

Figure 1. Comparison diagram of stellar masses measured via an SED fit of EVCC ugriz photometry with Le Phare (left), and star formation rates derived using total IR luminosities with those given in the SDSS MPA/JHU DR7 VAGC. The dashed lines in the top panels show the one-to-one relations. ing in an offset as large as ∼13.6800 in one particular a signal-to-noise ratio cut of S/N ≥ 3 in each band. case. We found that coordinate offsets occur for large To compute SFRs, we use WISE 22µm flux den- spiral galaxies given the presence of many clumps in sities, which were shown to be reliable SFR indicators their disks, despite the fact that the photometry given in Hwang et al. (2012c) and Lee et al. (2013). We de- by EVCC was carefully conducted. Given that the rive total infrared (IR) luminosities from WISE 22µm galaxies with the largest offsets were those with high- luminosities using the SED templates of Chary & El- resolution HI imaging used in Yoon et al. (2017), it was baz (2001). We then convert the total IR luminosities necessary to ensure that these galaxies could be found 00 into SFRs with the Kennicutt (1998) relation, assuming within a matching tolerance of 6 . We thus replaced a Salpeter IMF: SFR (M yr-1) = 1.72 × 10-10 L the coordinates of Yoon et al.’s 48 galaxies with those IR IR (L ). We remove active galactic nuclei (AGN) from given by the NASA Extragalactic Database (NED). For the sample to avoid contamination in the mid-infrared. 2 galaxies with no match, the NED-given coordinates We identify 21 AGN via the Baldwin-Phillips-Terlevich are replaced with WISE-given coordinates of the clos- (BPT) classification provided by the MPA/JHU DR7 est WISE detection in the AllWISE Source Catalog. VAGC. We also use the WISE mid-infrared color-color There were 6 additional galaxies in the EVCC that were selection criteria defined by Jarrett et al. (2011) and not found within the matching tolerance, which were Mateos et al. (2012) to identify one additional AGN. manually searched for and added to the master cata- A total of 24 AGN, including 2 composite galaxies log. The AllWISE Source Catalog provides both point with large offsets from the one-to-one relation in Figure source profile-fitting and elliptical aperture magnitudes. 1 (via 2-sigma clipping), were identified and removed Because we are dealing with nearby galaxies, some of from the sample. them are well-resolved to the point that we generally require larger apertures to obtain an accurate measure- We then compare our SFR estimates with those in ment of their total brightness. We thus prioritize the the MPA/JHU DR7 VAGC (Brinchmann et al. 2004) in use of elliptical aperture magnitudes for galaxies in the the right panel of Figure 1. The SFRs in the MPA/JHU EVCC. For those with no elliptical magnitudes given DR7 VAGC were derived from extinction and aperture- by the AllWISE Source Catalog (i.e., relatively small corrected Hα luminosities. For AGN and galaxies with galaxies), we use point source profile-fitting magnitudes weak emission lines, they measure the SFRs using the instead. We also select reliable flux density values with 4000 A( Dn4000) break. There seems to be a correlation between the two measurements despite a large scatter. SF Activity of Galaxies Undergoing RPS in Virgo 5

Table 1 The VIVA sample defined by different stages of gas stripping

Class defHI defHI HI morphology and deficiency compared to field galaxies (Range) (Median) (1) (2) (3) (4) Class 0 -0.81 ∼ 0.38 0.13 No definite signs of gas stripping due to the ICM (N = 13)

Class I -0.43 ∼ 0.41 0.02 One-sided HI feature; (N = 7) relatively comparable to field galaxies in HI gas content

Class II 0.12 ∼ 1.16 0.76 Highly asymmetric HI disk; quite deficient in HI, (N = 10) with an average of ∼17% of HI mass compared to those of field galaxies.

Class III 0.82 ∼ 2.25 1.42 Symmetric and severely truncated HI disk; (N = 10) extremely deficient in HI, with <4% of the HI mass of a field counterpart.

Class IV 0.51 ∼ 1.17 0.79 Symmetric HI disk with marginal truncation (N = 8) within the radius of the stellar disk; lower HI surface density than other classes; quite deficient in HI with on average ∼15% of the HI mass of a field counterpart.

(1) RPS class, (2) HI deficiency range, (3) median HI deficiency, (4) HI morphology and relative HI deficiency compared to those of field galaxies. Table adapted from Yoon et al. (2017).

We caution the reader that the two measurements are 2005; Haynes et al. 2011, 2018). The latest release, the not guaranteed to show a tight correlation, as it is not ALFALFA Extragalactic HI Source catalog (Haynes et ideal to derive physical parameters for galaxies in our al. 2018), provides HI data for a total of ∼31,500 ex- sample’s redshift range, due to the small size of the tragalactic sources detected out to z < 0.06. We cross- SDSS fibers. (Kewley et al. 2005). As such, the plot match the catalog with the EVCC within a matching should be referenced simply for comparison purposes. tolerance of 400, giving a total of 579 matched objects. The median offset between SFRIR and SFRSDSS is ap- Among the 579 matched objects, there are a total of 39 proximately 0.03 dex. The data with the largest offsets galaxies with RPS classes defined by Yoon et al. (2017). can primarily be explained by two factors. Upon vi- For the rest of the RPS class sample without HI masses sual inspection of SDSS images and examination of op- from ALFALFA, we adopt HI masses from the VIVA tical size estimates, we were able to conclude that most survey. The HI detection limit for the ALFALFA sur- 6 galaxies with measured SFRIR larger than SFRSDSS vey is about 10 M , and the stellar mass of our sam- 7 11 were either optically large or contained inner structures ple ranges from about 2×10 M to 6×10 M . This (e.g., bar, bulge). We attribute the likely cause of this indicates that the lowest valid value of the relative HI- discrepancy to the fact that SDSS fiber observations to-stellar mass ratio for our sample is approximately were limited to the central part of the galaxies with 2×10-6, which is equivalent to a logarithmic value of low star formation activity, thus leading to an underes- -5.8. Our lowest measured value gives a logarithmic timation of the SFR despite aperture corrections. On value of -2.97, suggesting that all our measurements of the other hand, we found that most galaxies with SFRIR the HI mass fractions are valid. smaller than SFRSDSS were blue irregular galaxies. This makes sense as UV bright galaxies that consist entirely 3.R ESULTS of young stellar populations (and thus do not contain 3.1. Star Formation Properties of Galaxies Undergoing much dust) would result in a very weak detection in Different Stages of RPS WISE 22µm. This would in turn result in a higher 3.1.1. Definition of RPS Classes measurement of the SFR in Hα and a lower measure- ment of the SFR in the infrared. Yoon et al. (2017) categorize a total of 35 galaxies into different stages of stripping and identify an additional class of 13 control galaxies that are not subject to strip- 2.2.3. Atomic Gas Masses ping. Characterization of each class is briefly described We adopt the atomic gas masses, MHI, from the Arecibo in Table 1, adapted from Table 1 and Section 2 in Yoon Legacy Fast ALFA Survey (ALFALFA, Giovanelli et al. et al. 2017. Among the stripping classes, Classes I, II, 6 Mun et al.

3 3 E Field 2 B 2

Recent 1 A 1 D l C l c c Ancient Intermediate First

/ 0 / 0

v v D C 1 A 1

B EVCC 2 2 Class 0 Class I E Class II 3 3 Class III Class IV

0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 Rproj / r200 Rproj / r200

Figure 2. (Left) Projected phase-space diagram with regions A ∼ E defined according to time since infall into the cluster center, closely following the classification in Rhee et al. (2017). Clustercentric distance and velocity are normalized by the cluster virial radius r200 and cluster velocity dispersion σcl, respectively. First infallers are located in Region A, recent infallers in Region B, intermediate infallers in Region C, and ancient infallers in Region D. Region E is defined by field galaxies. (Right) Projected phase-space diagram of RPS class galaxies and the rest of the EVCC sample shown as reference (grey open symbols). and III represent galaxies undergoing early, active, and diagram. post stripping stages, respectively. On the other hand, galaxies that show no strong evidence of stripping due In the left panel of Figure 2, we label different re- to the ICM are utilized as a control sample, labeled as gions in projected phase space according to time since Class 0 (N = 13), whereas galaxies that appear to have infall, using a modification of the classification adopted lost HI gas at all radii due to a mechanism other than by Rhee et al. (2017). According to the expected or- ram pressure stripping, such as starvation (Larson et bit of a radially infalling galaxy (e.g., Yoon et al. 2017; al. 1980), are labeled as Class IV (N = 8). Rhee et al. 2017; Jaff´eet al. 2015), a galaxy under the influence of the cluster potential will slowly accelerate Yoon et al. utilize these five classes in combination from the outskirts to eventually virialize in the cluster with projected phase space to better trace the orbital core after multiple core crossings. Classes I ∼ III over- histories of HI stripped galaxies within the Virgo clus- all follow the expected trend of radially infalling galax- ter. Given the large scatter in projected phase space ies within projected phase space, as shown in the right due to projection effects (e.g., Oman et al. 2013; Rhee et panel of Figure 2. Class I galaxies are located at high ve- al. 2017), the availability of high-resolution HI imaging locity and low to intermediate clustercentric distances, plays an important role in identifying the orbital stage which in combination with their HI morphology, show of the selected galaxies. We plot the locations of RPS that they are likely to be on their first infall into the class galaxies, along with other EVCC galaxies, in pro- cluster core. Class II galaxies are at similar locations in jected phase space in the right panel of Figure 2. In the projected phase space as Class I galaxies, with a select diagram, we can clearly see that our reference sample few found to be traveling at relatively low velocities. is not uniformly distributed in projected phase space, For the latter, Yoon et al. suggest that tidal interac- such that they preferentially inhabit regions with posi- tions may cause a shift in projected phase space, and tive clustercentric velocities. We attribute the cause of thus caution is advised when interpreting their locations this non-uniformity to the limited sky coverage of the with respect to their orbital histories. Although some EVCC. Yoon et al. find reference galaxies traveling at Class III galaxies are found in the same region in pro- positive velocities to be consistent with the NGC 5353/4 jected phase space as Class II/I objects, their heavily filament and the W-M sheet (Kim et al. 2016), which truncated HI disks suggest that they have already expe- are found to be located behind the Virgo cluster. On rienced their first core crossing, and are now falling out the other hand, Kim et al. (2016) also find filamentary of the core. HI morphology alone suggests that Class structures located in front of and/or in the vicinity of III galaxies are at a later stripping stage than Class II Virgo, to be elongated towards the cluster, indicating galaxies are. However, their relative locations in pro- cluster infall. However, these structures fall outside of jected phase space imply that a subsample of Class II the sky coverage of the EVCC, which explains the ab- galaxies are at a later stage in their orbital histories rel- sence of galaxies in the bottom half of the phase-space ative to some Class III galaxies. Yoon et al. therefore SF Activity of Galaxies Undergoing RPS in Virgo 7

102 EVCC Class 0 Class I 101

) Class II 1 Class III r

y Class IV

100 M (

R I R

F 10 1 S

10 2 0.2 0.4

g - r 0.6 0.8 2 1 0 1 2 [3.4] - [12] 108 109 1010 1011 1012 M* (M ) Figure 3. Total IR-derived star formation rates, optical, and mid-infrared colors (from top to bottom) as a function of stellar mass for Virgo galaxies with defined RPS classes from Yoon et al. (2017) and the remaining EVCC galaxies with no defined classes plotted as reference. The solid black lines define a reference sequence for each star formation activity tracer, to normalize for the effects of different masses. The reference sequence for SFRIR is the star-forming main sequence defined by Elbaz et al. (2007), and the dashed lines are the lower and upper bounds to the fit. The solid line in the middle panel is 10 the red sequence, which was derived based on a best linear fit of non-RPS Class galaxies with M* > 1.5×10 M and g − r > 0.58. The solid line in the bottom panel is adopted from the peak of the mid-infrared color distribution of star-forming galaxies in Hwang et al. (2012b). suggest that Class III galaxies do not always occur after ogy and location in projected phase space enriches the Class II. Among the Class III galaxies, there are a sub- information one can obtain with regards to the orbital sample of them found at radii beyond the cluster virial history of galaxies within a galaxy cluster. radius (i.e., Rproj/r200 = 1). Yoon et al. state that their low velocities with respect to the cluster, on top of them 3.1.2. Star Formation Properties of Galaxies with Different being located outside the cluster virial radius, suggest RPS Classes that they have fallen out of the cluster core after first We examine the star formation properties of galaxies infall, hence labeling them as a “backsplashing” popu- undergoing different stages of ram pressure stripping, lation. While Class IV galaxies are generally found at using star formation tracers sensitive to a broad range low to intermediate clustercentric distances, they are of star-forming timescales. We use SFRIR, g − r, and rather spread out in projected phase space. Class 0 WISE [3.4] − [12] colors. We plot these star formation galaxies also exhibit a large scatter in projected phase tracers against stellar mass in Figure 3. For SFRIR, we space, although more than half of them are found in the adopt the main sequence defined in Elbaz et al. (2007), first infalls (i.e., Region A) region. With the exception which was derived from a best linear fit of star-forming of galaxies subject to other stripping mechanisms such SDSS galaxies at 0.015 ≤ z ≤ 0.1. The solid line rep- as tidal interactions and starvation, or those that do resents the main sequence, with the dashed lines rep- not seem to show any signs of gas depletion, Yoon et resenting the lower and upper bounds to the main se- al. demonstrate that the combination of HI morphol- quence fit. The optical color g − r can also work as a 8 Mun et al.

Cluster members 200 20 50 Class 0 10 100 20

0 0 0 0 0 0 Class I 200 20 50 10 100 20

0 0 0 0 0 0 Class II 200 20 50 10 100 20 N N N

0 0 0 0 0 0 Class III 200 20 50 10 100 20

0 0 0 0 0 0 Class IV 200 20 50 10 100 20

0 0 0 0 0 0 10 4 10 3 10 2 10 1 100 101 0.4 0.2 0.0 0.2 0.4 4 2 0 2 RSB = sSFR/sSFRMS (g - r) ([3.4] - [12])

Figure 4. Starburstiness (RSB), ∆(g − r), ∆([3.4] − [12]) histograms (from left to right). Vertical lines denote the locus of the reference sequence (i.e., star-forming main sequence for RSB, red sequence for g − r, star-forming sequence for [3.4] − [12]). Certain members, for which cluster membership is given by EVCC, with no defined RPS classifications are plotted together for reference. The arrows represent the median values for each sample. Bin counts were increased by a factor of 5 (10 for the ∆(g − r) distribution) for RPS galaxies for a fair comparison with the reference sample. The scale of the original bin counts are labeled on the right hand side of each panel. The representative (either the mean or the median) errors for each parameter are shown in the top right corner of each panel. tracer of star formation activity (Strateva et al. 2001; RPS class galaxies are overall blue relative to the red Blanton et al. 2003). We define a red sequence using sequence, with Class IV galaxies noticeably redder in the best linear fit of non-RPS class galaxies with M* comparison. As for mid-infrared colors, Classes I and 10 > 1.5×10 M and g − r > 0.58. The mid-infrared II galaxies are found to be distributed around the star- [3.4] − [12] color is a tracer of recent star formation forming sequence, despite undergoing RPS. This seems activity in a galaxy, up to ∼2 Gyr timescales (e.g., Ko to contradict the notion given by their total IR-based et al. 2013). It is also known that 12µm and 3.4µm SFRs, as the two classes overall exhibit low SFRs with luminosities show strong correlations with SFRs and respect to the star-forming main sequence. This dispar- stellar masses, respectively. As such, [3.4] − [12] color ity can be explained by the fact that mid-infrared col- works as a proxy of SFRs per stellar mass (i.e. specific ors trace longer star-forming timescales than those of SFRs) (Donoso et al. 2012; Li et al. 2007). We adopt total IR-derived star formation rates (up to ∼100 Myr; the peak of the mid-infrared color distribution of star- Kennicutt & Evans 2012). On the other hand, Class forming galaxies (i.e., [3.4] − [12] = 1) from Hwang et III and IV galaxies overall exhibit lower values of mid- al. (2012b) as the star-forming sequence. infrared colors (primarily corresponding to lower 12µm emission), which is not surprising considering their HI Largely due to small number statistics of the RPS morphology and range of HI deficiencies. In contrast class sample, it is difficult to identify a trend in the star to the clustering of galaxies around the star-forming formation activity of RPS class galaxies from Figure sequence, we also identify a passive sequence located 3 alone. Overall, galaxies in the Virgo cluster appear around [3.4] − [12] ∼ -2.0. According to Figure 4c from to have lower current star formation rates relative to Ko et al. (2013), where they compare their sample of the star-forming main sequence. There are a select few quiescent red sequence galaxies with single stellar popu- Class 0 and I galaxies found above the main sequence. lation (SSP) models that include mid-infrared emission However, there are also some others that appear to be from AGB stars, the mid-infrared colors of the passive less star-forming compared to other RPS class galax- galaxies are roughly consistent with stellar populations ies that are far more deficient in HI. We note that the with a mean age of &5 Gyr. Those with [3.4] − [12] main sequence adopted here was based on a best linear . -2.0 can be described by SSP models without AGB fit of galaxies located at 0.015 ≤ z ≤ 0.1 (see Elbaz dust. Such galaxies with relatively old stellar popula- et al. 2007, Figure 18), and thus may not be an accu- tions are found to lie on the optical red sequence and rate representation of the main sequence for galaxies at exhibit much lower star formation rates. z . 0.01, which is the case for galaxies in the EVCC. As such, we do not use the main sequence as an ab- We note that RPS class galaxies exhibit different solute definition of the star-forming main sequence for ranges in stellar mass, which shows that there is a need the EVCC sample, but rather as a reference sequence to minimize mass effects. We thus choose to exam- for computing starburstiness. In terms of optical color, ine the distributions of starburstiness, ∆(g − r), and SF Activity of Galaxies Undergoing RPS in Virgo 9

∆([3.4] − [12]). Starburstiness (RSB = sSFR/sSFRMS; formation tracers systematically change with stripping Elbaz et al. 2011) is a measure of the relative excess classes, such differences are not strongly supported by of the specific star formation rate (sSFR = SFR/M*) statistical tests. We suspect that this is largely due to of a galaxy compared to that of a main sequence star- small number statistics, and the fact that we are using forming galaxy of the same mass. ∆(g − r) is a measure integrated photometry to trace the star formation ac- of how blue (i.e., star-forming) a galaxy is relative to tivity in these galaxies. Moreover, we also consider the a red sequence galaxy of the same mass, with negative possibility that if there is any enhancement in star for- values indicating bluer galaxies and positive values red- mation induced by ram pressure stripping, it is unlikely der galaxies. ∆([3.4] − [12]) measures the deficiency to last for a long enough period of time to be traceable in IR color (which is a proxy for deficiency in star for- by the star formation tracers used in this study. We mation) of a galaxy relative to a galaxy of the same discuss the implications of our results in further detail mass on the star-forming sequence, with positive values in Section 4. We observe the most statistically signifi- signifying galaxies rich in IR color. Using these param- cant differences for RPS classes with large disparities in eters thus allows a fair comparison of the star formation HI deficiencies, relative HI extent, and HI morphology, tracer distributions between different RPS classes. such as Class 0 and III, and Class 0 and IV. Figure 4 shows the resulting histograms of the star formation tracer distributions - starburstiness, ∆(g − 3.2. Star Formation Properties of Galaxies Classified by r), and ∆([3.4] − [12]) - from left to right. We denote HI Mass Fractions and Positions in Phase Space the median of each distribution with downward arrows. While we are able to capture the overall quenching of We overlay the distribution of each RPS class on top of star formation activity with increasing degree of gas that of other cluster members with no RPS classifica- stripping with RPS classes, the differences among the tion, as reference. In general, the resulting distributions stripping classes are not statistically significant. As of each star formation activity indicator seem to be con- high-resolution HI imaging is only available for a limited sistent with one another. Classes 0 and I seem to make sample of the EVCC, small number statistics do come up the actively star-forming end (i.e., higher starbursti- into play. Moreover, because Yoon et al. (2017)’s clas- ness, bluer colors, and higher mid-IR colors) of the clus- sification is primarily based on HI morphology, we do ter population, whereas Classes II and III tend to lie require a more quantitative classification scheme on top in the range of moderately to passively star-forming. of a larger sample. We thus aim to expand the sample On the other hand, Class IV galaxies appear to be the by re-defining new stripping groups using a combination most quenched in star formation activity for all star for- of HI mass fraction (HI-to-stellar mass ratios; MHI/M*) mation tracers. We also note the presence of galaxies and location in projected phase space. While Yoon et with lower starburstiness and ∆([3.4] − [12]) than RPS al. categorized galaxies into different stripping classes class galaxies, which can be attributed to the sequence based on HI properties, and then traced the orbital his- of quiescent galaxies observed in the bottom panel of tories of such galaxies using location in projected phase Figure 3. However, the differences between Classes I space, we utilize both HI mass fractions and expected through III do not show to be very clear, hinting at position in projected phase space according to radial in- the limitations of small number statistics. We perform fall to define a new classification scheme. We choose to Kolmogorov-Smirnov (K-S) tests to examine the star use HI-to-stellar mass ratios as they are not dependent formation activity distributions in more detail. on a separate sample of isolated galaxies, like HI de- We are not able to observe any noticeable enhance- ficiencies are. Moreover, given that HI-to-stellar mass ment in star formation activity for galaxies undergoing ratios are more commonly used to measure HI gas con- early to active stripping, i.e., Classes I and II. Despite a tent, we believe that HI-to-stellar mass ratios would be shift in the median of the ∆(g − r) distribution towards a better choice if one were to compare the results of bluer colors for Class I relative to that of Class 0, we this study to those of other studies. We visually in- are not able to reject the hypothesis that the two distri- spect all EVCC galaxies with HI detections to identify butions are drawn from the same population (i.e., the and remove a total of 11 spheroidal galaxies that may null hypothesis), even at a >2σ level, for the two classes have been subject to morphological quenching (Martig for any of the star formation activity indicators. Taking et al. 2009) from the sample. This leads to a total of into account that the HI deficiency ranges of Class 0 and 519 spiral and irregular galaxies with HI detections in I galaxies are overall comparable to one another (see Ta- the EVCC, which includes all 48 of Yoon et al.’s RPS ble 1), it is not surprising that the distinction between class galaxies. the two classes is statistically insignificant. Classes I and II also show no statistically significant differences 3.2.1. Definition of HI Stripping Groups from one another. On the other hand, we are able to re- To define new HI stripping groups, we first plot the HI ject the null hypothesis for Class 0 and III, and Class 0 mass fraction as a function of stellar mass, and perform and IV at a >2σ level for all star formation tracers. The a best linear fit of Class 0 and I galaxies to define a results of the K-S test on the ∆([3.4] − [12]) distribution reference sequence (Eq. 1), as shown in Figure 5. De- gives pKS ≈ 0.001 for Classes 0 and III, and for Classes spite differences in HI morphology, Class 0 and I are 0 and IV, indicating a significance of >3σ. Although similar in HI deficiencies (see Table 1) and thus the two there seems to be a trend that the distributions of star classes can be thought to represent the HI-rich end of 10 Mun et al.

HI-normal to rich Class 0 HI-poorish Class I 1 HI-poor 1 Class II Class III Class IV 0 0 * * M M

/ /

1 1 I I H H M M

g g o o

l 2 l 2

3 3

4 4

8 9 10 11 12 8 9 10 11 12 log M* log M*

Figure 5. (Left) HI mass fraction as a function of stellar mass, for which each of the shaded regions mark the range of HI mass fractions of each group, as defined by Eqs. 2 through 4. The blue line is the reference sequence as defined by a best linear fit of Class 0 and I galaxies, defined by Eq. 1. (Right) HI mass fraction as a function of stellar mass for the RPS class sample, to verify whether the selected ranges of ∆(log MHI/M*) for each group are consistent with those of the RPS class sample.

galaxies in the Virgo cluster. As similarly done in Fig- 0.7 ure 4, we measure ∆(log MHI/M*), the deficiency of the HI mass fraction of a galaxy with respect to that of a galaxy on the reference sequence with the same mass. 0.6 We then use Gaussian Mixture Modeling (GMM) to identify three groups defined by richness in their HI gas content: HI-normal to rich, HI-poorish, and HI-poor. 0.5 The GMM results are shown in Figure 6. We adopt the intersections of the Gaussian components, i.e., the 0.4 dashed grey lines in Figure 6, as the boundaries between the different groups. The equations defining the ranges of ∆(log MHI/M*) for each group are given by Eqs. 2, 0.3 3, and 4, denoting HI-normal to rich, HI-poorish, and Probability HI-poor, respectively.

0.2 MHI log = −0.915 log M∗ + 8.478 (1) M∗ 0.1 HI-normal to rich:

0.0 MHI (2) 2 1 0 log + 0.915 log M∗ ≥ (8.478 − 0.42) M∗ (log MHI / M*) HI-poorish: Figure 6. ∆(log MHI/M*) histogram with the resulting Gaus- sian Mixture Model (GMM) overlaid on top. Each Gaussian ( MHI log + 0.915 log M∗ < (8.478 − 0.42) (3) component specifies the mean and range of a potential HI M∗ MHI log + 0.915 log M∗ ≥ (8.478 − 1.179) group. The blue component represents the HI-normal to M∗ rich galaxies, the green represents the HI-poorish galaxies, and the red the HI-poor. The grey dashed lines mark the HI-poor: intersections between the different Gaussian components. MHI (4) log + 0.915 log M∗ < (8.478 − 1.179) M∗ SF Activity of Galaxies Undergoing RPS in Virgo 11

HI-normal-EAB 3 HI-poorish-ABC E HI-poor-BCD HI-normal to rich 2 HI-poorish B HI-poor A 1

l D c C

/

0 v D C 1 B A 2 E 3

0.0 0.5 1.0 1.5 2.0 2.5 3.0 Rproj / r200

Figure 7. Projected phase-space diagram of all newly defined HI stripping groups. The black solid lines indicate the boundaries between different regions of time since infall as shown in Figure 2. The final group candidates are shown as filled color symbols, whereas those that did not satisfy the condition of expected location in phase space are shown as open symbols. Red open symbols (HI-poor galaxies located outside of Regions B through D) residing in the region defined by grey solid lines are selected as potential backsplashing candidates, which are discussed briefly in Section 4.2.2. Table 2 Categorization of the EVCC into different HI stripping groups

HI Stripping Group ∆(log MHI/M*) Location in phase space (Range) (1) (2) (3) HI-normal-EAB Eqn. 2 Regions E, A, or B (N = 132) HI-poorish-ABC Eqn. 3 Regions A, B, or C (N = 168) HI-poor-BCD Eqn. 4 Regions B, C, or D (N = 65)

(1) HI stripping group, (2) HI mass fraction range, and (3) location in phase space.

With these definitions, we classify galaxies in the thus affirm that the selected ranges of ∆(log MHI/M*) EVCC with HI detections into a total of three groups, for each group, as given by the results of GMM, are well which are visualized in the left panel of Figure 5. On founded. the right, we verify whether the selected ranges of HI To distinguish galaxies at different stages in their mass fractions are consistent with those of the RPS class orbital histories, we employ the use of projected phase sample. Despite the presence of some scatter, most of space as a second means of categorizing galaxies into HI Class 0 and I galaxies are well represented as HI-normal stripping groups, as shown in Figure 7. We utilize the to rich galaxies. The selected range of HI mass frac- same definitions of regions in phase space as shown in tions for the HI-poorish group is consistent with those the left panel of Figure 2. Galaxies that are likely to be of Class II and III galaxies, and we confirm that more on their first infall are expected to be found in Region than half of the Class III objects are categorized as HI- A. For those that are approaching their first core cross- poor galaxies. With the exception of one galaxy, Class ing, or have already passed and are falling out of the IV galaxies are generally well described as HI-poorish, core, they are likely to be found in Region B. Galax- which is expected given that their HI deficiency range ies that have gone through multiple core crossings will overlaps with that of Class II galaxies (see Table 1). We eventually settle and virialize within the cluster poten- 12 Mun et al. tial, finding themselves at low velocities and cluster- subgroup, 1088 km s-1 (Cluster A; Mei et al. 2007). centric distances, such as Regions C or D. However, as We adopt values of r200 = 1.55 Mpc (McLaughlin 1999; -1 confirmed by Yoon et al. (2017), it is possible that some Ferrarese et al. 2012) and σcl = 593 km s (Mei et al. galaxies may reach further out after outfall, otherwise 2007). known as backsplashing galaxies. We briefly examine the possibility of finding backsplashing candidates with 3.2.2. Star Formation Properties of New HI Stripping our group categorization in Section 4.2.2. Results from Groups Jaff´eet al. (2015) suggest that star-forming galaxies We examine the star formation properties of the newly relatively rich in HI gas content are found to be recent defined HI stripping groups, using a larger sample. Us- infallers, whereas passive galaxies with low gas content ing the same star formation tracers as before, we show are found in the virialized zone. This suggests that our results in Figure 8. As seen similarly with Figure galaxies are expected to lose their gas and become in- 3, we are not able to observe a specific trend in the HI creasingly quenched in star formation activity as they stripping groups, especially given that a disparity in the fall into the cluster core. As for galaxies that are located range of stellar mass is apparent between them. For at relatively high velocities and large clustercentric dis- example, a large fraction of HI-poorish-ABC and HI- tances, they are less subject to the cluster potential poor-BCD galaxies appear to dominate the low-mass and thus are expected to retain a comparable amount end, whereas the high-mass end seems to be increas- of HI mass as that of a field galaxy. We may find such ingly dominated by HI-normal-EAB galaxies. In fact, galaxies in Region E. However, we stress that the HI- there are some HI-normal-EAB galaxies found below normal to rich galaxies located in Region E may not be the star-forming main sequence in the low-mass end. an ideal control sample of galaxies outside of the clus- Upon investigating them, we find that they belong to ter potential well, as those presumed to be background galaxy groups, which suggests that they are undergo- galaxies were excluded from the EVCC with a radial ing pre-processing prior to cluster infall (Kourkchi & velocity cut, as mentioned in Section 2. Keeping the Tully 2017). To minimize mass effects, we perform the typical infall trajectory along with the expected trend same in-depth analysis as done in Section 3.1.2, where in star formation activity and gas content in mind, we we minimize biases due to mass effects by plotting the finalize our HI group sample based on their expected histograms of starburstiness, ∆(g − r), and ∆([3.4] − position in phase space, in combination with their HI [12]). Moreover, for a fair comparison between the dis- mass fractions, resulting in a total of 365 galaxies. The tributions of different star formation tracers, we per- 8 10 final HI group sample is presented in Table 2, labeled as form a stellar mass cut of 2×10 M ≤ M* ≤ 6×10 HI-normal-EAB, HI-poorish-ABC, and HI-poor-BCD. M , which is roughly the stellar mass range of galaxies While there were many other possible combinations of with measured SFRIR (as plotted in the top panel of the two parameters we could have used, our classifi- Figure 8), excluding a select few outliers on the high- cation scheme was chosen such that galaxies of each mass end. group fell in the expected region of phase space as those The resulting histograms are shown in Figure 9. of corresponding RPS classes (see Figure 2). For ex- The starburstiness distribution shows that the median ample, HI-normal-EAB galaxies were chosen such that values are similar for HI-normal-EAB and HI-poorish- they roughly lie in the same region as Class 0 and I ABC groups, while the median of the HI-poor-BCD galaxies. HI-poorish-AB were meant to represent Class group is shifted towards lower values of RSB. The K-S II galaxies, whereas HI-poor-BCD represent Class III. tests reveal that all 3 groups show differences from one Amongst the final sample, there is a total of 43 galax- another on the level of >2σ. We believe the lack of ies (out of 48) that comprise the RPS class sample in a highly significant discrepancy between the 3 groups Yoon et al. (2017). There are 19 galaxies categorized may be due to the fact that the WISE 22µm detec- as HI-normal-EAB, which consist of RPS Classes 0 - II, tion limit is not deep enough to result in a large sam- and IV. The HI-poorish-ABC group includes 17 RPS ple of galaxies compared to those of other star forma- class galaxies, which consist of a mixture of Classes I - tion tracers. On the other hand, we see a much clearer IV. The HI-poor-BCD group consists of 6 RPS Class III trend in the change in star formation activity for ∆(g galaxies. The remaining galaxy is one of the confirmed − r); the median ∆(g − r) increases towards redder backsplashing candidates in Yoon et al. (2017), which colors with decreasing HI gas content. Our K-S test is included in our analysis of potential backsplashing results suggest that we can reject the null hypothesis candidates in Section 4.2.2. at a >3σ level for the ∆(g − r) distribution for all 3 When plotting the projected phase-space diagram, groups. In contrast, we observe slightly different re- we take on the same assumptions as done in Yoon et al. sults with the ∆([3.4] − [12]) distribution. There is (2017). Here, we assume that the center of the Virgo a shift in the median ∆([3.4] − [12]) towards larger cluster is at M87 (B¨ohringeret al. 1994). We normal- values, which can suggest an enhancement of star for- ize the projected clustercentric distance by r200 and the mation activity in HI-poorish-ABC galaxies, relative to line-of-sight clustercentric velocity by the cluster disper- that of HI-normal-EAB galaxies. We then observe a sion, σcl. The clustercentric velocity is calculated us- quenching of the star formation activity of HI-poor- ing the difference between a galaxy’s velocity along the BCD galaxies with the median shifted towards nega- line-of-sight and the mean radial velocity of the M87 tive values of ∆([3.4] − [12]). However, the supposed SF Activity of Galaxies Undergoing RPS in Virgo 13

102 HI-normal-EAB HI-poorish-ABC HI-poor-BCD 101 ) 1 r y

100 M (

R I

R 1 F 10 S

10 2 0.0

0.5 g - r

2 1 0 1 [3.4] - [12] 108 109 1010 1011 1012 M* (M ) Figure 8. Star formation rates derived from total IR luminosities, optical, and mid-infrared colors (from top to bottom) as a function of stellar mass for the EVCC sample with newly defined HI stripping groups. The reference sequences shown here are identical to the ones used in Figure 3.

40 30 HI-normal-EAB 60 HI-poorish-ABC 20 40 20 10 20

0 0 0 30 HI-poor-BCD 40 8 10 2 × 10 M M* 6 × 10 M 60 2 × 108M M 6 × 1010M 8 10 * 2 × 10 M M* 6 × 10 M 20 40 N N 20 N 10 20

0 0 0 40 30 Backsplashing 60

20 40 20 10 20

0 0 0 10 2 10 1 100 101 0.2 0.0 0.2 0.4 2 1 0 1 2 RSB = sSFR/sSFRMS (g - r) ([3.4] - [12])

Figure 9. Starburstiness (RSB), ∆(g − r), ∆([3.4] − [12]) histograms (from left to right) for new HI stripping groups. The downward arrows mark the median of each distribution. The vertical lines here are the same lines defined in Figure 4. Stellar mass range of the subsamples used are displayed in the middle panel of each histogram. Additionally, we show the distribution of potential backsplashing candidates in the bottom panel. 14 Mun et al. enhancement of star formation activity observed by the galaxy NGC 4294 (Class I). Furthermore, pixel-by-pixel shift in the median is not expected to be meaningful analyses of FUV images by Boissier et al. (2012) show as the K-S test does not yield significant results and little to no star formation activity in the stripped HI the difference in the median values are only slightly out gas of a select few Class I and II galaxies, with NGC of the range of error. On the other hand, we observe 4522 (Class II) being an exception. However, the au- a statistically significant difference between HI-normal- thors come to the conclusion that despite the presence EAB and HI-poor-BCD galaxies, and HI-poorish-ABC of a mild enhancement in star formation, it is unlikely and HI-poor-BCD galaxies with a p-value < 0.004 for to affect the global star formation activity and thus will both. In general, our results suggest that there is no not be discernible when considering integrated quanti- statistically significant difference between HI-normal- ties or large samples of galaxies as a whole. As RPS EAB and HI-poorish-ABC galaxies. However, we are class galaxies show to be undergoing quenching on a able to reject the null hypothesis at a >3σ level for the global scale as a result of ongoing gas stripping, any optical color distribution. We also observe a significant local enhancement in star formation is unlikely to be difference between HI-poorish-ABC and HI-poor-BCD discernible with integrated photometry. This explains galaxies for all star formation tracers. Overall, we are why, on top of small number statistics, the differences able to capture the quenching of star formation activ- in star formation activity between different RPS classes ity with increasing deficiency in HI gas content with a are not statistically significant. larger sample compared to our analysis in Section 3.1.2. On the other hand, we observe a significant dif- ference on a >3σ level between Classes 0 and III, and 4.D ISCUSSION Classes 0 and IV. Use of integrated photometry works 4.1. Locally Enhanced Star Formation in Galaxies particularly well in cases where galaxies showing large Undergoing Gas Stripping disparities in HI gas content and star formation activity We examine the star formation activity of galaxies un- are compared with one another. Moreover, Yoon et al. dergoing ram pressure stripping in the Virgo cluster us- (2017) suggest that Class IV galaxies are likely to have ing various tracers. We initially perform a follow-up undergone starvation (Larson et al. 1980), given that study of Yoon et al. (2017), given the availability of they show gas depletion at all radii in the disk instead high-resolution HI-imaging thanks to the VIVA survey of the expected asymmetry arising from ram pressure (Chung et al. 2009). While our results confirm the pre- stripping in their HI morphology. Such galaxies would viously known trend of star formation quenching with have their star formation quenched throughout the disk, increasing degree of ram pressure stripping (Jaff´eet which would be apparent when traced with integrated al. 2015; Yoon et al. 2017), the difference between the properties. distributions of galaxies undergoing early through post While enhancements in star formation may be stripping is not statistically significant. However, we prominent for galaxies undergoing gas stripping in other do observe a statistically significant distinction between galaxy clusters (e.g., Vulcani et al. 2018; Ramatsoku et Class 0 and III, and Class 0 and IV galaxies. al. 2019), it is evident that the data used in this study do not show such a global enhancement for galaxies in 4.1.1. Implications from Observational Studies Based on the Virgo cluster. If one intends to capture the local en- the Virgo Cluster hancement of star formation in Virgo galaxies, spatially The lack of a conspicuous enhancement in star forma- resolved data would be necessary, as done in Koopmann tion activity is not too surprising given the results of & Kenney (2004b). Results of Lee et al. (2017) further previous studies on the Virgo cluster. Koopmann & prove the advantage of using spatially resolved data - Kenney (2004b) utilize the spatial distributions of Hα they confirm strong enhancement in both far-ultraviolet and R-band emission to investigate the environmental (FUV) and Hα along CO compression within the stel- effects on the star formation activity on 52 galaxies in lar disk for three Class II galaxies, which support our the Virgo cluster, among which 39 of them are RPS hypothesis that galaxies undergoing active stripping are class galaxies. Most of the Hα-truncated spirals show likely to experience some sort of star formation enhance- signs of having been subject to ram pressure stripping, ment. where many of them exhibit normal to slightly enhanced In addition, Kenney et al. (2014) suggests that star formation rates in the inner galaxy disks. While whereas galaxies with one-sided tails of young stellar they find asymmetric enhancements in star formation knots are often observed in more distant and massive at the outer edge of the Hα disk for a few galaxies, (i.e., denser ICM and higher peak ram and thermal including NGC 4654 (Class II) and NGC 4405 (Class pressures) clusters such as the Coma cluster, the Virgo III), which is indicative of ongoing ICM pressure, these cluster itself is not massive enough to completely strip enhancements are at most globally modest and local. spiral galaxies. However, such features can be observed They instead find that galaxies with the largest global in lower mass galaxies such as dwarf irregulars, as Ken- enhancement in star formation are those that have un- ney et al. (2014) show with their analysis of IC 3418. dergone low-velocity tidal interactions and/or gas ac- IC 3418 was not included in the EVCC, and moreover cretion. An example of such a galaxy in their sample has been shown to have weak HI (Chung et al. 2009; is NGC 4299 (Class I), which shows evidence of under- Kenney et al. 2014) and mid-infrared emission (checked going both RPS and tidal interactions due to its pair with WISE 12µm and 22µm images). Thus, our overall SF Activity of Galaxies Undergoing RPS in Virgo 15 results are unlikely to have changed significantly even sensitive to timescales of ∼100 Myr (Kennicutt & Evans if we had included dwarf irregular galaxies such as IC 2012), which implies that our choice of star formation 3418, given our choice of star formation tracers. tracers was appropriate. However, it is entirely possible that any enhancement due to RPS may occur for even 4.1.2. Impact of RPS on the Global Star Formation Activity shorter timescales, such that it would be more advan- of Galaxies as Predicted by Simulations tageous to use star formation tracers sensitive to even Results of previous observational studies, along with shorter timescales (∼10 Myr), such as Hα. Barring that ours, suggest that a global enhancement in star forma- possibility, our choice of star formation tracers do not tion triggered by RPS is rather mild and challenging to seem to be of an issue in terms of sensitivity to star detect when considering integrated quantities or large formation timescales. samples of galaxies as a whole. Studies such as Kenney It is also worth noting that previous studies have et al. (2014) suggest that the degree of enhancement found that while RPS can enhance the global SFR for depends on a variety of factors, such as the stellar mass albeit a short period of time, such changes were not of the galaxy, total mass of the host galaxy cluster, discernible in galaxy colors. For example, Fujita & Na- density of the ICM, etc. Moreover, recent simulation gashima (1999) study the change in the SFR along with studies (e.g., Oman & Hudson 2016; Rhee et al. 2017) a galaxy’s optical B − V color to find that while the tracing the trajectories of infalling cluster galaxies in SFR increases by at most a factor of 2, the color stays phase space have shown that the relative velocity of overall the same. The authors suggest that ram pres- galaxies, along with the ambient ICM density, changes sure compression does not induce observable star forma- with the galaxy’s location with respect to the cluster tion activity as the intrinsic scatters of color and lumi- center. As such, the impact of RPS on the star forma- nosity among galaxies are larger. In fact, after stripping tion activity is not a trivial matter, given that it is also occurs, the color of their model galaxy rapidly turns red. closely coupled with the orbital history of the galaxy. Kapferer et al. (2009) also report similar results, in that Disagreements with regards to the degree or even despite significant enhancement of the overall SFR due the possibility of global enhancement are also found in to ram pressure, a blue star-forming disk galaxy is ex- simulation studies, given the strong dependency on ini- pected to transform into a red, passively evolving sys- tial conditions. Recent RPS simulations such as Bekki tem under <250 Myr. Steinhauser et al. (2016) studied (2014) and Steinhauser et al. (2016) suggest that global the SDSS u − i color evolution of their model galaxies enhancement of star formation induced by ram pres- to find similar results as well, with the rate of evolution sure is modest at best (e.g., not on the order of a factor towards redder colors dependent on the strength of ram of ∼10). Furthermore, despite disagreements with one pressure. As such, while it is feasible to trace the overall another in the specific conditions for enhancement to quenching of star formation activity with galaxy colors, occur, both Bekki (2014) and Steinhauser et al. (2016) results of previous studies seem to suggest that galaxy show that global enhancement only takes place under colors are not sensitive to short-term enhancements of rigorous conditions. As such, the number of galaxies un- the SFR. dergoing enhancement likely make up a small fraction of the cluster population, making it less distinguishable 4.2. Depletion of HI Gas Content According to Time when considering the cluster population as a whole. Since Infall in Phase Space Just as galaxies undergoing even a modest degree As the RPS classification given in Yoon et al. (2017) is of global enhancement can be rare in galaxy clusters, as limited to galaxies with high-resolution HI imaging, we insinuated by recent simulation studies, it is also worth- build a larger sample of galaxies undergoing ram pres- while exploring the possibility that the enhancement sure stripping in the cluster by using a combination timescale is rather short, which may bring the choice of of HI mass fractions and positions in projected phase star formation tracers used in this study into question. space. We examine a total of three different groups of However, there is no clear consensus with regards to galaxies in the EVCC, with the final sample totaling the timescale over which enhancement lasts as reported to 365: 1) HI-normal-EAB: HI-normal to rich galaxies by simulation studies. Despite this, one could never- located at large clustercentric distances or high veloc- theless try to constrain the time period for which the ities, 2) HI-poorish-ABC: HI-poorish galaxies likely to enhancement stage can be best captured. Given that be on their first infall into the cluster core, and 3) HI- the compression of interstellar gas due to ram pressure poor-BCD: HI-poor galaxies that are in the process of is expected to be strongest at orbital pericenter (e.g., virializing within the cluster core. Vollmer 2009; Bekki 2014), it would be reasonable to as- sume that a modest starburst could occur around this 4.2.1. Overall Quenching of Star Formation with stage of a galaxy’s orbital history. Using models based Decreasing HI Gas Content and Increasing Time on galaxies in the Virgo cluster, Vollmer (2009) claims Since Infall that the time window during which one could identify Given the disparity in HI gas content and location in perturbations due to RPS around peak ram pressure is projected phase space, one may expect to see relatively ∼300 Myr. With this in mind, RPS itself is expected to quenched star formation in HI-poorish-ABC galaxies occur on a relatively short timescale of 100 Myr (e.g., compared to that of the HI-normal-EAB galaxies. How- Roediger 2009). SFRs based on TIR luminosities are ever, as mentioned in Section 3.2.2, the two groups are 16 Mun et al. not necessarily distinct from one another. One possible ters has long been discussed since Balogh et al. (2000) explanation for the lack of a significant discrepancy is and Mamon et al. (2004). After the first core crossing, that while it is expected for HI-poorish-ABC galaxies to galaxies may travel out as far as 2.5 times the virial ra- become quenched as they infall into the cluster core, the dius of the cluster (Gill et al. 2005). Such galaxies are number of galaxies that have actually begun quenching then expected to make another subsequent infall into may not comprise of a large fraction of the group. In the cluster core, eventually joining the galaxies that fact, while HI-poorish-ABC galaxies are likely to have are virializing in the cluster core. We select a poten- begun losing their gas due to the ICM, the effects of tial group of backsplashers by choosing HI-poor galax- gas loss could have yet to be reflected in the tracers of ies residing in the region |∆v|/σcl ≤ 1 and Rproj/r200 ≥ star formation activity used in this study. If the latter 1.5, resulting in a total of 49 galaxies (refer to Fig- is true, then it is not surprising that the star forma- ure 7). Amongst the selected candidates, there is one tion activity of HI-poorish-ABC and HI-normal-EAB Class III galaxy present, NGC 4064, which was con- galaxies are hardly distinguishable from one another. firmed as a backsplasher in Yoon et al. (2017). Yoon et Simulation results reported in Oman & Hudson (2016) al. identified a total of four backsplashing galaxies in suggest that star formation quenching occurs after a de- projected phase space, of which NGC 4457 is classified lay time of typically 3.5 - 5 Gyr, measured from the first as HI-poorish-ABC; NGC 4580 and NGC 4293 are HI- crossing of 2.5 rvir (∼3.4 r200), which can be thought poor-BCD and do not fall within the selected region of of as the first time the galaxy enters the cluster. This backsplashing candidates. The confirmed backsplash- onset of quenching usually corresponds to times near ing candidates in Yoon et al. (2017) are Class III galax- or shortly after first pericentric approach. As a signifi- ies located near ∆v/σcl ∼ 0 and at roughly 1.2 - 1.7 cant fraction of HI-poorish-ABC galaxies still have yet r200 in Figure 2. Compared to these galaxies, we have to make their first pericentric passage, it is likely that chosen a sample of backsplashing candidates that are they have yet to have their star formation quenched, overall travelling at relatively larger velocities with re- which can explain the lack of discrepancy in their star spect to the cluster, and some are found even further formation activity with that of HI-normal-EAB galax- out to the outskirts around 2.5 - 3 r200. We examine ies. On the other hand, we are able to confirm with a the star formation activity of these candidates in com- >2σ confidence that HI-poor-BCD galaxies are likely to parison to the three HI stripping groups. We show the be at a different stage in their orbital histories and thus resulting distributions of each star formation tracer in are a distinct population of galaxies within the Virgo the bottom panels of Figure 9. We perform the same cluster, compared to HI-normal-EAB and HI-poorish- stellar mass cut for each star formation tracer as de- ABC galaxies. scribed in Section 3.2.2 for consistency. While we defined sharp boundaries between differ- Overall, the selected backsplashing group seems ent regions in projected phase space, categorized ac- to represent an intermediate population between HI- cording to time since infall as shown by the simulation poorish-ABC and HI-poor-BCD galaxies. Taking into results of Rhee et al. (2017), we stress that the locus account that backsplashing may occur past the clus- of these boundary lines are not to be taken as absolute ter’s virial radius after first core passage, the trend seen dividing lines between distinct populations of galaxies in Figure 9 seems to make sense. Although the trend in phase space. They are merely to be used as a means seen in the starburstiness distribution is unclear due of identifying galaxies likely to be undergoing ram pres- to the limited number of galaxies with SFR estimates sure stripping given their HI gas fractions in projected (i.e., WISE 22µm detections), we are able to observe phase space. To confirm that our results are not sig- increased star formation quenching as a galaxy makes nificantly affected by the locus of the boundary lines multiple pericentric approaches within the Virgo clus- in phase space, we perform tests by making different ter. If they are true backsplashers, we surmise that they modifications to the boundaries. For example, we have will eventually join the HI-poor-BCD galaxies near the tried lowering the upper limit of the boundary line of cluster core, and become passive as they continue to lose Region E from ∆v/σcl = 3 to 2.5, which led to a de- their gas. We also consider the possibility that these crease in the selection of HI-poorish-ABC galaxies. We galaxies may not be a true backsplashing population, also tightened the selection of HI-poor-BCD galaxies by and may very well have lost their gas through means limiting Regions B and C to lower velocities and clus- other than ram pressure stripping. When galaxies are tercentric distances. In spite of these alterations, there accreted into a nearby galaxy cluster, some may fall in were no significant changes observable in the star for- as an isolated galaxy, but it is also likely that a signifi- mation distributions. We thus confirm that our results cant fraction of them may fall in as a group. Simulation 14.5 -1 are not significantly changed even if we make different results of McGee et al. (2009) show that a 10 h M modifications to the boundaries of different regions in cluster at z = 0 has accreted about 40% of its galaxies 9 -1 phase space. (M* > 10 h M ) from groups with masses greater 13 -1 than 10 h M . Prior to infall, galaxies in groups 4.2.2. Potential Backsplashing Population are likely to have their star formation quenched due We also examine a potential backsplashing population to mechanisms such as galaxy-galaxy interactions and of galaxies in the Virgo cluster. The existence of a galaxy-group tidal interactions, a phenomenon known backsplashing population in the outskirts of galaxy clus- as “pre-processing” (Fujita 2004). These galaxies thus SF Activity of Galaxies Undergoing RPS in Virgo 17 may appear to have lost even more gas with respect to do not rule out the possibility that these galaxies their location in projected phase space. Taking into ac- might have lost their gas through pre-processing count the presence of filamentary structures such as the rather than ram pressure stripping. Leo II A and B filaments (Kim et al. 2016) around the vicinity of the Virgo cluster, pre-processing is definitely We utilized the availability of multi-wavelength a possibility that we cannot rule out when it comes data for galaxies in the Virgo cluster to study the star to interpreting the gas content of galaxies in projected formation properties tracing different timescales. How- phase space. ever, the parameters we extracted were based on inte- grated photometry, which limited a detailed study of 5.S UMMARYAND CONCLUSIONS the effects of ram pressure stripping on the ISM of the We investigate the star formation properties of galaxies galaxies. To study the comprehensive picture of ram undergoing ram pressure stripping in the Virgo clus- pressure stripping on the star formation activity of clus- ter, building upon the analysis of the orbital histories ter galaxies, it is needed to make use of high spatial res- of VIVA galaxies with high-resolution HI imaging as olution data obtainable via integral field spectroscopy reported in Yoon et al. (2017). We use starburstiness, (IFS) observations. IFS data allows spatially resolved optical, and mid-infrared colors to examine the star for- mapping of emission lines from which one can obtain mation activity of galaxies undergoing different stages diagnostic line ratios to constrain the ionization mech- of ram pressure stripping. Due to the limited number of anism of the ISM, i.e., distinguish whether the observed galaxies with high-resolution HI imaging, we extend our emission originates from AGN feedback or star forma- study to the EVCC sample using HI mass fraction and tion (e.g., Merluzzi et al. 2013). In particular, simulta- position in projected phase space to categorize galaxies neous observations of Hα and Hβ emission lines enable into different stages of stripping. We summarize our the derivation of dust extinction maps that can be stud- main results below. ied along with Hα-derived SFR maps to study the star formation activity of galaxies on local scales (Nelson et 1. We identify a trend of decreasing star formation ac- al. 2016; Jafariyazani et al. 2019). tivity with increasing degree of stripping, as con- firmed by previous studies. However, we do not ACKNOWLEDGMENTS find strong evidence for star formation enhance- ment in galaxies undergoing early through active We thank the referees for their insightful comments stripping. We attribute the lack of a strong dis- that helped improve the paper. This study was sup- cernment towards the insensitivity of integrated ported by the National Research Foundation of Korea photometry towards local star formation enhance- (NRF) grant funded by the Korean Government (NRF- ments and the possibility of a relatively short en- 2019R1A2C2084019). Funding for the SDSS and SDSS- hancement stage. II has been provided by the Alfred P. Sloan Founda- 2. To make up for the shortcomings of small num- tion, the Participating Institutions, the National Sci- ber statistics, we extend our study towards the ence Foundation, the U.S. Department of Energy, the EVCC sample, using HI mass fractions and posi- National Aeronautics and Space Administration, the tion in projected phase space. We identify three HI Japanese Monbukagakusho, the Max Planck Society, stripping groups - HI-normal to rich, HI-poorish, and the Higher Education Funding Council for Eng- and HI-poor - represented by different ranges of HI land. The SDSS Web Site is http://www.sdss.org/. mass fractions. We then perform a second means The SDSS is managed by the Astrophysical Research of categorization by expected location in projected Consortium for the Participating Institutions. The Par- phase space according to radial infall into the clus- ticipating Institutions are the American Museum of ter core. Natural History, Astrophysical Institute Potsdam, Uni- versity of Basel, University of Cambridge, Case West- 3. While we are not able to observe a strong discern- ern Reserve University, University of Chicago, Drexel ment in the star formation properties between the University, Fermilab, the Institute for Advanced Study, HI-normal-EAB and HI-poorish-ABC groups, the the Japan Participation Group, Johns Hopkins Univer- HI-poor-BCD group shows to be clearly quenched sity, the Joint Institute for Nuclear Astrophysics, the in star formation with respect to the other two Kavli Institute for Particle Astrophysics and Cosmol- groups. We attribute the lack of an obvious dis- ogy, the Korean Scientist Group, the Chinese Academy tinction between HI-normal-EAB and HI-poorish- of Sciences (LAMOST), Los Alamos National Labora- ABC galaxies to the late onset of star formation tory, the Max-Planck-Institute for Astronomy (MPIA), quenching. the Max-Planck-Institute for Astrophysics (MPA), New 4. We investigate a potential backsplashing popula- Mexico State University, Ohio State University, Univer- tion, selected from HI-poor galaxies located at low sity of Pittsburgh, University of Portsmouth, Princeton velocities and large clustercentric distances beyond University, the United States Naval Observatory, and the virial radius of the cluster. Distributions of the University of Washington. This publication makes their star formation activity show that they may use of data products from the Wide-field Infrared Sur- represent an intermediate population between HI- vey Explorer, which is a joint project of the University poorish-ABC and HI-poor-BCD galaxies. We also of California, Los Angeles, and the Jet Propulsion Lab- 18 Mun et al. oratory/California Institute of Technology, and NEO- Chung, A., van Gorkom, J. H., Kenney, J. D. P., et al. 2009, WISE, which is a project of the Jet Propulsion Lab- VLA Imaging of Virgo Spirals in Atomic Gas (VIVA). I. oratory/California Institute of Technology. WISE and The Atlas and the H I Properties, AJ, 138, 1741 NEOWISE are funded by the National Aeronautics and Conroy, C. 2013, Modeling the Panchromatic Spectral En- Space Administration. ergy Distributions of Galaxies, ARAA, 51, 393 Cortese, L., Marcillac, D., Richard, J., et al. 2007, The REFERENCES strong transformation of spiral galaxies infalling into mas- sive clusters at z ∼ 0.2, MNRAS, 376, 157 Abazajian, K., Adelman-McCarthy, J. K., Ag¨ueros,M. A., Cowie, L. L. & Songaila, A. 1977, Thermal evaporation of et al. 2004, The Second Data Release of the Sloan Digital gas within galaxies by a hot intergalactic medium, Na- Sky Survey, AJ, 128, 502 ture, 266, 501 Abazajian, K. N., Adelman-McCarthy, J. K., Ag¨ueros, M. A., et al. 2009, The Seventh Data Release of the Sloan Crowl, H. H. & Kenney, J. D. P. 2008, The Stellar Popu- Digital Sky Survey, ApJS, 182, 543 lations of Stripped Spiral Galaxies in the Virgo Cluster, Arnouts, S., Cristiani, S., Moscardini, L., et al. 1999, Mea- AJ, 136, 1623 suring and modelling the redshift evolution of clustering: Donoso, E., Yan, L., Tsai, C., et al. 2012, Origin of 12 the Hubble Deep Field North, MNRAS, 310, 540 µm Emission across Galaxy Populations from WISE and Balogh, M. L., Navarro, J. F., & Morris, S. L. 2000, The Ori- SDSS Surveys, ApJ, 748, 80 gin of Star Formation Gradients in Rich Galaxy Clusters, Elbaz, D., Daddi, E., Le Borgne, D., et al. 2007, The rever- ApJ, 540, 113 sal of the star formation-density relation in the distant Bekki, K., Couch, W. J., & Shioya, Y. 2002, Passive Spiral universe, A&A, 468, 33 Formation from Halo Gas Starvation: Gradual Transfor- Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, mation into S0s, ApJ, 577, 651 GOODS-Herschel: an infrared main sequence for star- Bekki, K. 2014, Galactic star formation enhanced and forming galaxies, A&A, 533, A119 quenched by ram pressure in groups and clusters, MN- Ferrarese, L., Cˆot´e,P., Cuillandre, J.-C., et al. 2012, The RAS, 438, 444 Next Generation Virgo Cluster Survey (NGVS). I. Intro- Bernardi, M., Shankar, F., Hyde, J. B., et al. 2010, Galaxy duction to the Survey, ApJS, 200, 4 luminosities, stellar masses, sizes, velocity dispersions as Fujita, Y. & Nagashima, M. 1999, Effects of Ram Pressure a function of morphological type, MNRAS, 404, 2087 from the Intracluster Medium on the Star Formation Rate Bertin, E., & Arnouts, S. 1996, SExtractor: Software for of Disk Galaxies in Clusters of Galaxies, ApJ, 516, 619 source extraction., A&AS, 117, 393 Fujita, Y. 2004, Pre-Processing of Galaxies before Entering Blanton, M. R., Hogg, D. W., Bahcall, N. A., et al. 2003, a Cluster, PASJ, 56, 29 The Broadband Optical Properties of Galaxies with Red- Gill, S. P. D., Knebe, A., & Gibson, B. K. 2005, The evolu- shifts 0.02

Jafariyazani, M., Mobasher, B., Hemmati, S., et al. 2019, the Virgo cluster. I. How far can galaxies bounce out of Spatially Resolved Properties of Galaxies from CAN- clusters?, A&A, 414, 445 DELS+MUSE: Radial Extinction Profile and Insights on Martig, M., Bournaud, F., Teyssier, R., et al. 2009, Morpho- Quenching, ApJ, 887, 204 logical Quenching of Star Formation: Making Early-Type Jaff´e,Y. L., Smith, R., Candlish, G. N., et al. 2015, BUD- Galaxies Red, ApJ, 707, 250 HIES II: a phase-space view of H I gas stripping and star Mateos, S., Alonso-Herrero, A., Carrera, F. J., et al. 2012, formation quenching in cluster galaxies, MNRAS, 448, Using the Bright Ultrahard XMM-Newton survey to de- 1715 fine an IR selection of luminous AGN based on WISE Jarrett, T. H., Cohen, M., Masci, F., et al. 2011, The colours, MNRAS, 426, 3271 Spitzer-WISE Survey of the Ecliptic Poles, ApJ, 735, 112 McGee, S. L., Balogh, M. L., Bower, R. G., et al. 2009, The Kapferer, W., Sluka, C., Schindler, S., et al. 2009, The effect accretion of galaxies into groups and clusters, MNRAS, of ram pressure on the star formation, mass distribution 400, 937 and morphology of galaxies, A&A, 499, 87 McLaughlin, D. E. 1999, Evidence in Virgo for the Universal Kauffmann, G., Heckman, T. M., White, S. D. M., et al. Dark Matter Halo, ApJL, 512, L9 2003, Stellar masses and star formation histories for 105 Mei, S., Blakeslee, J. P., Cˆot´e,P., et al. 2007, The ACS galaxies from the Sloan Digital Sky Survey, MNRAS, 341, Virgo Cluster Survey. XIII. SBF Distance Catalog and the 33 Three-dimensional Structure of the Virgo Cluster, ApJ, Kennicutt, R. C. 1998, The Global Schmidt Law in Star- 655, 144 forming Galaxies, ApJ, 498, 541 Merluzzi, P., Busarello, G., Dopita, M. A., et al. 2013, Kennicutt, R. C., & Evans, N. J. 2012, Star Formation in ACCESS - V. Dissecting ram-pressure stripping through the Milky Way and Nearby Galaxies, ARAA, 50, 531 integral-field spectroscopy and multiband imaging, MN- Kenney, J. D. P., Geha, M., J´achym, P., et al. 2014, Trans- RAS, 429, 1747 formation of a Virgo Cluster Dwarf Irregular Galaxy by Moore, B., Katz, N., Lake, G., et al. 1996, Galaxy harass- Ram Pressure Stripping: IC 3418 and Its Fireballs, ApJ, ment and the evolution of clusters of galaxies, Nature, 780, 119 379, 613 Kewley, L. J., Jansen, R. A., & Geller, M. J. 2005, Aperture Nelson, E. J., van Dokkum, P. G., Momcheva, I. G., et al. Effects on Star Formation Rate, Metallicity, and Redden- 2016, Spatially Resolved Dust Maps from Balmer Decre- ing, PASP, 117, 227 ments in Galaxies at z ∼ 1.4, ApJL, 817, L9 Kim, S., Rey, S.-C., Jerjen, H., et al. 2014, The Extended Oman, K. A., Hudson, M. J., & Behroozi, P. S. 2013, Disen- Virgo Cluster Catalog, ApJS, 215, 22 tangling satellite galaxy populations using orbit tracking Kim, S., Rey, S.-C., Bureau, M., et al. 2016, Large-scale Fil- in simulations, MNRAS, 431, 2307 amentary Structures around the Virgo Cluster Revisited, Oman, K. A. & Hudson, M. J. 2016, Satellite quenching ApJ, 833, 207 time-scales in clusters from projected phase space mea- Ko, J., Hwang, H. S., Lee, J. C., et al. 2013, The Mid- surements matched to simulated orbits, MNRAS, 463, infrared and Near-ultraviolet Excess Emissions of Quies- 3083 cent Galaxies on the Red Sequence, ApJ, 767, 90 Park, C., & Hwang, H. S. 2009, Interactions of Galaxies in Koopmann, R. A., & Kenney, J. D. P. 2004, Massive Star the Galaxy Cluster Environment, ApJ, 699, 1595 Formation Rates and Radial Distributions from Hα Imag- Poggianti, B. M., Moretti, A., Gullieuszik, M., et al. 2017, ing of 84 Virgo Cluster and Isolated Spiral Galaxies, ApJ, GASP. I. Gas Stripping Phenomena in Galaxies with 613, 851 MUSE, ApJ, 844, 48 Koopmann, R. A., & Kenney, J. D. P. 2004, Hα Morpholo- Ramatsoku, M., Serra, P., Poggianti, B. M., et al. 2019, gies and Environmental Effects in Virgo Cluster Spiral GASP - XVII. H I imaging of the jellyfish galaxy JO206: Galaxies, ApJ, 613, 866 gas stripping and enhanced star formation, MNRAS, 487, Kourkchi, E. & Tully, R. B. 2017, Galaxy Groups Within 4580 3500 km s-1, ApJ, 843, 16 Rhee, J., Smith, R., Choi, H., et al. 2017, Phase-space Anal- Kronberger, T., Kapferer, W., Ferrari, C., et al. 2008, On ysis in the Group and Cluster Environment: Time Since the influence of ram-pressure stripping on the star forma- Infall and Tidal Mass Loss, ApJ, 843, 128 tion of simulated spiral galaxies, A&A, 481, 337 Roediger, E. 2009, Ram pressure stripping of disk galaxies Kroupa, P. 2001, On the variation of the initial mass func- in galaxy clusters, Astronomische Nachrichten, 330, 888 tion, MNRAS, 322, 231 Salpeter, E. E. 1955, The Luminosity Function and Stellar Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, The Evolution., ApJ, 121, 161 evolution of disk galaxies and the origin of S0 galaxies, Smith, R. J., Lucey, J. R., Hammer, D., et al. 2010, Ultra- ApJ, 237, 692 violet tails and trails in cluster galaxies: a sample of can- Lee, B., Chung, A., Tonnesen, S., et al. 2017, The effect of didate gaseous stripping events in Coma, MNRAS, 408, ram pressure on the molecular gas of galaxies: three case 1417 studies in the Virgo cluster, MNRAS, 466, 1382 Steinhauser, D., Schindler, S., & Springel, V. 2016, Simula- Lee, J. C., Hwang, H. S., & Ko, J. 2013, The Calibration of tions of ram-pressure stripping in galaxy-cluster interac- Star Formation Rate Indicators for WISE 22 µm-Selected tions, A&A, 591, A51 Galaxies in the Sloan Digital Sky Survey, ApJ, 774, 62 Strateva, I., Ivezi´c, Z.,ˇ Knapp, G. R., et al. 2001, Color Sep- Li, H.-N., Wu, H., Cao, C., et al. 2007, Morphological De- aration of Galaxy Types in the Sloan Digital Sky Survey pendence of Mid-Infrared Properties of SDSS Galaxies in Imaging Data, AJ, 122, 1861 the Spitzer SWIRE Survey, AJ, 134, 1315 Vollmer, B. 2009, A holistic view on ram pressure stripping Mamon, G. A., Sanchis, T., Salvador-Sol´e,E., et al. 2004, in the Virgo cluster. The first complete model-based time The origin of HI-deficiency in galaxies on the outskirts of sequence, A&A, 502, 427 20 Mun et al. von der Linden, A., Wild, V., Kauffmann, G., et al. 2010, Star formation and AGN activity in SDSS cluster galax- ies, MNRAS, 404, 1231 Vulcani, B., Poggianti, B. M., Gullieuszik, M., et al. 2018, Enhanced Star Formation in Both Disks and Ram- pressure-stripped Tails of GASP Jellyfish Galaxies, ApJL, 866, L25 West, A. A., Garcia-Appadoo, D. A., Dalcanton, J. J., et al. 2010, H I-Selected Galaxies in the Sloan Digital Sky Survey. I. Optical Data, AJ, 139, 315 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, The Wide-field Infrared Survey Explorer (WISE): Mission Description and Initial On-orbit Performance, AJ, 140, 1868 Yoon, H., Chung, A., Smith, R., et al. 2017, A History of HI Stripping in Virgo: A Phase-space View of VIVA Galax- ies, ApJ, 838, 81