Open Access
Issue
A&A
Volume 615, July 2018
Article Number A114
Number of page(s) 16
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201732410
Published online 24 July 2018

© ESO 2018

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0;), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

The environment plays a fundamental role in shaping galaxy evolution. Since the work of Dressler (1980), it is clear that clusters of galaxies are dominated by early-type systems (ellipticals and lenticulars) which are relatively rare in the field (see also Whitmore et al. 1993). It is also well established that late-type galaxies in dense environments have a lower atomic gas content than their isolated counterparts (Haynes & Giovanelli 1984; Cayatte et al. 1990; Solanes et al. 2001; Gavazzi et al. 2005) with several recent indications that this occurs also for the molecular gas phase (Fumagalli et al. 2009; Boselli et al. 2014a) and the dust content (Cortese et al. 2010, 2012a). The gas removal, which begins in the outer regions and progresses inwards, leading to truncated discs (Cayatte et al. 1994; Boselli et al. 2006), induces a decrease of the star formation activity, as is systematically observed in nearby clusters (e.g. Kennicutt 1983; Gavazzi et al. 1998; Gómez et al. 2003; Boselli et al. 2014b).

Since the discovery of these systematic differences between galaxies in rich environments and those in the field, observers, modelers and simulators have made great strides in identifying the dominant physical mechanism responsible for the gas removal and for the subsequent quenching of the star formation activity, as reviewed in Boselli & Gavazzi (2006, 2014). The detailed observations of local galaxies with strong signs of an ongoing perturbation (Gavazzi et al. 2001; Kenney et al. 2004, 2014; Boselli et al. 2006, 2016a; Vollmer et al. 2006, 2008a,b, 2009, 2012; Sun et al. 2007; Scott et al. 2012), the analysis of local samples of cluster galaxies (Gavazzi et al. 1998, 2010, 2013; Boselli et al. 2008a,b, 2014b, 2016b), as well as finely tuned models and simulations of representative objects (Roediger & Hensler 2005; Roediger & Brüggen 2007, 2008; Tonnesen & Bryan 2009) suggest that the dominant process in local clusters is the ram pressure (Gunn & Gott 1972) exerted by the dense intracluster medium on the interstellar medium of galaxies moving at high velocity within the cluster. On the other hand, the statistical analysis of large samples of galaxies extracted from the SDSS spanning a wide range of density, combined with hydrodynamic cosmological simulations, favour a less violent interaction able to remove only the hot gas halo and thus stop the the infall of fresh material, thus reducing the activity of star formation only once the cold gas on the galaxy discshas been fully transformed into stars (starvation – Larson et al. 1980). Other processes, such as galaxy harassment (Moore et al. 1998), might have significantly contributed during the accretion of galaxies in massive clusters via small groups (pre-processing, Dressler 2004) and may still be dominant under specific conditions.

So far, however, there have been no systematic efforts to study the fate of the stripped gas. There are two major reasons for this situation: 1) the difficulty in observing the gas once removed from the galactic disc, needed to provide stringent constraints on models and simulations; and 2) the difficulty in simulating the dynamical evolution of the gas properties within a complex and inhomogeneous cluster environment (density, turbulence, temperature) in its different phases (atomic, molecular, ionised, hot) and at different scales (from molecular clouds on pc scales to tails 50–100 kpc long). One of the first detections of stripped material in cluster late-type galaxies was the observations of ~50 kpc-long radio continuum tails (synchrotron emission) in three galaxies at the periphery of A1367 (Gavazzi & Jaffe 1985, 1987; Gavazzi et al. 1995). Their cometary shape has been interpreted as clear evidence of an ongoing ram pressure stripping event. Deep narrow-band Hα observations of the same galaxies have highlighted long tails of ionised gas (Gavazzi et al. 2001), suggesting that the cold atomic hydrogen might change phase once in contact with the hot intracluster medium via, for example, thermal evaporation (Cowie & Songaila 1977). Indeed, the observations of tails of neutral atomic gas (HI) are still quite uncommon: only a few cases are known in nearby clusters (Virgo – Chung et al. 2007; A1367 – Scott et al. 2012). The situation changed recently thanks to the advent of large panoramic detectors mounted on 4-m-class telescopes with narrow-band filters that allowed the detection of tails of ionised gas after very deep exposures in several galaxies in Virgo (Yoshida et al. 2002; Kenney et al. 2008; Boselli et al. 2016a), Coma (Yagi et al. 2010; Fossati et al. 2012), Norma (Sun et al. 2007; Zhang et al. 2013), and A1367 (Gavazzi et al. 2001; Boselli & Gavazzi 2014; Yagi et al. 2017), and a further observational constraint came from X-ray observations, which allowed the detection of hot gas within these tails (Sun et al. 2006). Information is also provided by long-slit (Yoshida et al. 2012; Yagi et al. 2013) or IFU spectroscopic observations (Fumagalli et al. 2014; Fossati et al. 2016; Poggianti et al. 2017; Bellhouse et al. 2017; Fritz et al. 2017; Consolandi et al. 2017), which are fundamental for understanding the kinematics of the stripped gas and its chemical composition and physical state. Recently, also CO observations have been made possible for the detection of the molecular gas phase (Jáchym et al. 2013, 2014, 2017; Verdugo et al. 2015).

A surprising result of these recent studies is that only in a few cases does the stripped gas collapse to form stars outside the disc of the perturbed galaxy. This generally happens within compact H II regions dominated by young stellar populations (Hester et al. 2010; Fumagalli et al. 2011, 2014; Fossati et al. 2016; Consolandi et al. 2017). There are indeed several instances within the Virgo cluster where the stripped gas remains diffuse in its atomic neutral (Boissier et al. 2012) or ionised (Boselli et al. 2016a) phase and does not collapse to form new stars, while in other objects a star formation event occurs (Hester et al. 2010; Fumagalli et al. 2011; Arrigoni Battaia et al. 2012; Kenney et al. 2014). This observational evidence is in contrast with the results of models and simulations which systematically predict the formation of new stars outside the galaxy disc (Kapferer et al. 2009; Tonnesen & Bryan 2010, 2012), an indication that the physical prescriptions used in the models still need to be refined.

VESTIGE (A Virgo Environmental Survey Tracing Ionised Gas Emission; Boselli et al. 2018, Paper I) is a deep narrow-band Hα imaging survey of the Virgo cluster carried out with MegaCam at the CFHT. This survey is providing us with a unique opportunity to study the fate of the stripped gas in cluster galaxies, being the first complete narrow-band imaging survey of a nearby cluster up to its virial radius (covering an area of 104o2). Thanks to its unique sensitivity (f(H α) ~ 4 × 10−17 erg s−1 cm−2 − 5σ detection limit for point sources and Σ (H α) ~ 2 × 10−18 erg s−1 cm−2 arcsec−2 − 1σ detection limit at 3 arcsec resolution for extended sources) and subarcsecond image quality, VESTIGE will detect the diffuse ionised gas stripped from the perturbed galaxies, as well as compact H II regions of luminosity L(H α) ≥ 1036 erg s−1. VESTIGE will also benefit from the large number of available multifrequency data for the Virgo cluster, from the X-ray to the radio, which are necessary for a coherent and complete analysis of the multiphase gas (Boselli et al. 2018). The use of Hα data with respect to other star formation tracers is crucial if we want to study the process of star formation on timescales of ≲10 Myr (Boselli et al. 2009; Boquien et al. 2014). This is necessary whenever the perturbing process is rapid (≲500 Myr), as indeed expected within a rich cluster such as Virgo (Boselli et al. 2016b).

This paper is dedicated to the study of the spiral (SA(s)c) galaxy NGC 4254 (M99) located at the northern periphery of the cluster at ≃1.0 Mpc from M87 (≃ 0.7 × Rvir; see Fig. 1). The presence of a tail of HI gas ≃250 kpc in length indicates that the galaxy has been gravitationally perturbed after a rapid encounter with another Virgo cluster member (Haynes et al. 2007). The interaction is thought to be responsible for the creation of a free-floating HI cloud of gas not associated to any stellar component: the so called dark galaxy Virgo HI21 (Minchin et al. 2005). Dedicated simulations suggest that this interaction occurred ~280–750 Myr ago (Vollmer et al. 2005; Bekki et al. 2005; Duc & Bournaud 2008). The same simulations indicate that, during a rapid galaxy–galaxy encounter (galaxy harassment – Moore et al. 1998), only the gaseous component is perturbed leaving the stellar disc intact (Duc & Bournaud 2008). What happens to the stripped gas? Why is it still in its cold atomic phase? Being still cold, can it collapse and form new stars outside the galaxy disc and, if so, under what conditions? NGC 4254 is thus an ideal laboratory for studying the star formation process in extreme and unusual conditions, and another candidate to extend previous studies to different environments (Boquien et al. 2007, 2009, 2010, 2011; Fumagalli et al. 2011; Arrigoni Battaia et al. 2012; Jáchym et al. 2013, 2014, 2017; Kenney et al. 2014; Lisenfeld et al. 2016).

In this paper, we search for and study the properties of extraplanar H II regions formed after the interaction of the galaxy within the hostile cluster environment. We do this using the unique set of multifrequency data sensitive to the youngest stellar populations available for this representative galaxy combined with SED fitting models. The narrow-band Hα observationsand data reduction are described in Sect. 2. The multifrequency dataset used in the analysis is described in Sect. 3, while Sects. 4 and 5 describe the identification of the extraplanar H II regions and the derivation of their physical parameters, respectively. A discussion is given in Sect. 6. Given its position within the cluster, we assume the galaxy to be at a distance of 16.5 Mpc (Gavazzi et al. 1999; Mei et al. 2007; Blakeslee et al. 2009). All magnitudes are given in the AB system.

thumbnail Fig. 1

Colour ugi RGB image of the galaxy NGC 4254 obtained using the NGVS data (Ferrarese et al. 2012). At the distance of the galaxy (16.5 Mpc), 1 arcmin = 4.8 kpc. The star-forming regions analysed in this work are the blue blobs in the south-west outside the disc of the galaxy.

2 Observations and data reduction

The galaxy NGC 4254 was observed during the second run of pilot observations of the VESTIGE survey in spring 2016 (see Boselli et al. 2018 for details). Observations were carried out with MegaCam at the CFHT using the new narrow-band filter MP9603 (λc = 6590Å; Δ λ = 104Å, T = 93%) which includes observations at the redshift of the galaxy (vel = 2404 km s−1), the Hα line (λ = 6563Å), and the two [NII] lines (λ = 6548, 6583Å)1. The camera is composed of 40 CCDs with a pixel scale of 0.187 arcsec pixel−1. The galaxy was observed using the pointing-macro QSO LDP-CCD7 especially designed for the Elixir-LSB data reduction pipeline. This macro is composed of seven different pointings overlapping over a region of 40 × 30 arcmin2. To secure the determination of the stellar continuum, the galaxy was also observed in the r-band filter. The integration time for each single pointing within the macro was of 66 s in the r-band and 660 s in the narrow-band filter. The macro was run two times, thus the total integration time in the central 40 × 30 arcmin2 of the combined image is of 924 s in r-band and 9240 s in narrow-band.

Data reduction was completed using the standard procedures adopted for the VESTIGE data (extensively described in Paper I). The MegaCam images were reduced using Elixir-LSB, a data reduction pipeline optimised to detect the diffuse emission of extended low-surface-brightness features associated with the perturbed galaxies by removing any contribution of scattered light in the different science frames. This procedure, tuned to detect low-surface-brightness features, is perfectly adapted for the narrow-band frames whenever the images are background dominated, as is the case for those obtained in this work (see Boselli et al. 2016a). The photometric calibration of the r-band filter was derived using the standard MegaCam calibration procedures (Gwyn 2008). Those in the narrow-band filter were derived as described in Fossati et al. (in prep.). The typical uncertainty is 0.01 mag in both bands. The single images were stacked and an astrometric correction was applied using the MegaPipe pipeline (Gwyn 2008). As for the VESTIGE survey, the sensitivity is of f(H α) ~ 4 × 10−17 erg s−1 cm−2 for point sources (5σ) and Σ (H α) ~ 2 × 10−18 erg s−1 cm−2 arcsec−2 for extendedsources (1σ detection limit at ~3 arcsec resolution), while the depth in the r-band is 24.5 AB mag (5σ) for point sources and 25.8 AB mag arcsec−2 (1σ) for extended sources. The mean image quality is of 0.60 arcsec in the r-band and 0.56 arcsec in the narrow-band filter.

3 Multifrequency data

The Virgo cluster region has been the target of several multifrequency blind surveys. These data are crucial for the identification of the stellar population dominating any extraplanar star-forming region.

3.1 GALEX

Two very deep GALEX exposures are available for the galaxy NGC 4254 in the FUV (λc = 1539Å, Δ λ = 442Å; integration time 18 131 s) and NUV (λc = 2316Å. Δ λ = 1060Å; integration time 27 726 s) bands (Boselli et al. 2011). For these exposure times, the limiting magnitude for point sources is FUV ~ 25.3 AB mag and NUV ~ 24.9 AB mag (Voyer et al. 2014), while for extended sources is ~29.7 mag arcsec−2 in both bands. The GALEX images have a pixel size of 1.5 arcsec, and an angular resolution of 4.0 arcsec FWHM in the FUV and 5.6 arcsec in the NUV. The photometric accuracy of the instrument is of 0.05 AB mag in FUV and 0.03 in NUV, while the astrometric accuracy is of 0.5 arcsec in both bands (Morrissey et al. 2007).

3.2 Optical

Very deep optical images in the u, g, i, and z bands are available from the NGVS survey taken with MegaCam at the CFHT (Ferrarese et al. 2012). Their depth is ~1 mag deeper than the VESTIGE r-band images. The photometric and astrometric accuracy, and angular resolution of NGVS in these photometric bands is comparable to that reached in the r- and in the narrow-band filters.

3.3 Near- and far-IR

NGC 4254 has been observed in the near- and far-infrared by Spitzer with IRAC at 3.6, 4.5, 5.8, and 8.0 μm (Ciesla et al. 2014) and MIPS at 24 and 70 μm (Bendo et al. 2012) during the SINGS survey (Kennicutt et al. 2003), and by Herschel with PACS at 100 and 160 μm (Cortese et al. 2014) and SPIRE at 250, 350, and 500 μm (Ciesla et al. 2012) as part of the HeViCS (Davies et al. 2010) and HRS (Boselli et al. 2010) surveys. Near-IR data are also available from the WISE survey (Wright et al. 2010). Given the compact nature of the extraplanar H II regions analysed in this work, and their weak emission, the following analysis will use only the four IRAC and MIPS 24 μm bands, where the angular resolution and the sensitivity of the instruments allow the detection or the determination of stringent upper limits for the H II regionsobserved in the UV and optical bands. None of the sources has been detected at longer wavelengths.

3.4 HI and CO

The discovery of an ongoing perturbation of the disc of NGC 4254 due to a flyby encounter with another Virgo cluster member hasbeen made possible by HI observations. VLA observations first revealed the presence of HI blobs located along a diffuse feature starting at the south-west edge of the stellar disc and extending ~30 kpc to the north-west (Phookun et al. 1993). The galaxy has been observed in HI also at Effelsberg (Vollmer et al. 2005), Parkes (Wong et al. 2006), and with the VLA (Chung et al. 2009). More recent low-angular-resolution HI observations taken with the Arecibo radio telescope during the ALFALFA (Giovanelli et al. 2005) survey of the Virgo cluster indicate that this tail of HI gas is much more extended than previously thought (~250 kpc, to be compared to the isophotal radius of the galaxy which is 13 kpc) and links NGC 4254 to the tidal debris Virgo HI21 (Davies et al. 2004; Minchin et al. 2005; Bekki et al. 2005). The total mass of HI gas in this tail is estimated at (4.8 ± 0.6) × 108 M and its column density Σ (HI) ~ 3 × 1018 cm−2 (Haynes et al. 2007), while the HI blobs in the tail detected at the VLA have Σ (HI) ~ 1019 cm −2 at a resolution of ≃30 arcsec. The asymmetric distribution of the gas over the disc of the galaxy, with the presence of a low surface brightness tail in the northern part (Phookun et al. 1993), and the presence of a radio continuum extended tail (Kantharia et al. 2008), suggest that the galaxy is also suffering ram pressure stripping.

NGC 4254 has also been observed in 12CO(2–1) during the HERACLES survey of nearby galaxies using the IRAM 30m radio telescope, with an angular resolution of 13 arcsec, a velocity resolution of 2.6 km s−1, and a typical RMS of 20–25 mK (TMB scale; Leroy et al. 2009). Unfortunately, the final map covers only approximately one-third of the extraplanar H II regions of NGC 4254, and does not detect any of them.

4 Extraplanar H II regions

Figure 2 shows the multifrequency images of NGC 4254. The FUV and NUV images indicate the presence of several blue compact objects in the south-west periphery of the galaxy extending up to ~20 kpc outside the optical disc. These star forming complexes have already been noticed by Thilker et al. (2007) who identified them as part of a standard extended UV disc. Among these we selected, in the FUV GALEX image, those located at the edge of the prominent south-west spiral arm or outside the stellar disc of NGC 4254, with a compact morphology, and avoiding extended diffuse structures. We then rejected those with typical morphologies of background galaxies as seen in the higher-quality optical images. We identify 60 of these compact regions as depicted in Fig. 3. These objects are all detected in the NUV image, 30 of them also in the VESTIGE Hα continuum-subtracted image, while 34, 26,24, 21, and 16 in the ugriz images, respectively. The VESTIGE and NGVS images indicate that these sources are very compact since they are only barely resolved in these subarcsecond resolution frames. Although this sample is not complete in any sense, we are confident that it is statistically representative of the extraplanar young stellar regions associated with NGC 4254.

The FUV GALEX image alone does not allow us to make robust conclusions on whether these blue compact regions are all associated to the galaxy or are foreground or background objects. However, their relative distribution with respect to NGC 4254 and to the HI tail of stripped gas detected either in the VLA (Chung et al. 2009) or in the ALFALFA HI data (Haynes et al. 2007), as shown in Fig. 4, strongly suggests that they are star forming regions formed within the stripped gas of the galaxy. It is surprising that these star forming regions are only present within ≲20 kpc of the edge of the stellar disc of the galaxy, while they are totally absent in the rest of the HI tail. We can estimate the contamination of background objects using the FUV GALEX number counts derived by Xu et al. (2005). At the limiting magnitude of our detections (FUV ≃ 24 mag) the number of background sources is ~1000 deg−2, thus ~6–7 over the ~ 6′ × 4′ region outside NGC 4254 analysed in this work (~10%). We also measured the mean background source density down to a limiting magnitude of FUV ≤ 24 mag and detected with a signal-to-noise ratio (SN) > 5 within the whole GALEX FUV frame (0.5° radius) once the bright Virgo galaxies (NGC 4254, NGC 4262, NGC 4298) were masked. This density is ~2300 deg−2, thus the expected contamination is ~15 objects in the studied region. Since the background distribution slightly changes within the FUV frame, we also measured the density of FUV-emitting sources with properties similar to those analysed in this work in a corona centred on NGC 4254 of inner radius 210 arcsec and outer radius 400 arcsec, where all the studied regions are located. Considering as reference the density within this corona, the number of expected contaminants is approximately ten objects.

The coordinates of these regions are given in Table A.1. To study the stellar populations of these compact regions we extract their fluxes in the different photometric bands using the apertures shown in Fig. 3. All these apertures have a radius rap ≳ 4 arcsec to encompass the total flux of the emitting region in the different photometric bands, including those that are diffraction limited (FUV, NUV, 24 μm Spitzer/MIPS; see Table A.12). At the same time they are sufficiently small (rap ≤ 10 arcsec, where rap is the radius of theregions plotted in Figs. 2 and 3 encompassing all the FUV emission in the GALEX frame) to limit the contamination of background or local point sources. We notice, however, that despite the small size of the aperture, in some cases the high-resolution MegaCam images in Hα or in the broad bands reveal the presence of more than a single source unresolved by GALEX. None of these regions is detected at wavelengths longer than 24 μm. To obtain the most accurate estimate of the flux of each emitting region, or measure a stringent upper limit, the sky background is derived in the same apertures randomly displaced around the region itself once the other sources, and the diffuse emission of the galaxy is masked. This procedure is repeated 1000 times, minimising any uncertainty on the sky background determination. As in Fossati et al. (2018), the uncertainties on the fluxes are obtained as the quadratic sum of the uncertainties on the flux counts and the uncertainties on the background (rms of the bootstrap iterations). The uncertainties on the flux counts are derived assuming a Poissonian distribution for the source photons. We assume those sources with SN > 3 as being detected. For the undetected sources, we estimate an upper limit to the flux. Flux densities, corrected for Galactic extinction following Schlafly & Finkbeiner (2011) (E(BV) = 0.034) combined with the reddening law of Fitzpatrick (1999), with their uncertainty and 1σ upper limits, are listed in Table A.1. For consistency with the other photometric bands, the Hα flux is given as LyC (in units of μJy), the flux in the Lyman continuum pseudo filter PSEUDOLyC derived as described in Boselli et al. (2016b): LyC [mJy]=1.07×1037×L(Hα) [erg s1]D2 [Mpc].\begin{equation*} {\textrm{LyC}\,\rm{[mJy]} = \frac{1.07 \times 10^{-37} \times {L(\textrm{H}\alpha)} ~\rm{[erg ~s^{-1}]}}{ \it{D}^2 ~\rm{[Mpc]}}}.\end{equation*}(1)

thumbnail Fig. 2

Multifrequency images of the galaxy NGC 4254 (north is up, east is left). Upper panels, from left to right: continuum-subtracted Hα (VESTIGE), FUV (GALEX), NUV (GALEX); lower panels: u-band (NGVS), IRAC 3.6 μm (Spitzer), MIPS 24 μm (Spitzer). The extraplanar star-forming regions marked with green circles are evident in FUV image at the south-west of the galaxy. The linear size of each single image corresponds 40 × 35 kpc2 at the distance of the galaxy.

thumbnail Fig. 3

FUV/NUV GALEX colour (left) and continuum-subtracted Hα (VESTIGE; right) magnified images of the extraplanar star-forming regions (marked with green circles) in the south-west quadrant of NGC 4254. The UV colour image shows that these regions have very blue colours, and are thus dominated by young stellar populations (≲100 Myr). Half of them, however, are undetected in Hα, indicating that their typical age is 10 Myr ≲ age ≲ 100 Myr.

thumbnail Fig. 4

VLA HI contours at a column density Σ(HI) = 1019 cm−2 (from Chung et al. 2009; red) and the Arecibo HI contours at a column density Σ (HI) = 1−2 × 1019 cm−2 (from Haynes et al. 2007; white) are superposed on the GALEX FUV∕NUV colour image. The right panel is a magnified view of the boxed region marked in the left panel with a dashed line. The typical velocity of the different HI structures detected at the VLA are given. The yellow and cyan ellipses at the lower left corner of the two images show the beam sizes at Arecibo (3.3′ × 3.8′) and at the VLA (30.86′′ × 28.07′′).

5 Physical parameters

5.1 Hα luminosities

The narrow-band filter encompasses the Hα line at λ = 6563Å and the two [NII] lines at λ = 6548 and λ = 6583Å. To measure the Hα luminosity of these extraplanar H II regions shown in Fig. 3 we have first to remove the [NII] contribution. Since no spectroscopic data are available, we apply a mean standard correction. The typical [NII]λ6583/Hα ratio in extragalactic H II region is [NII]λ6583/Hα ≲ 0.3 (e.g. McCall et al. 1985; Kewley et al. 2001; Sánchez et al. 2015), with a ratio decreasing radially from the nucleus to the outer disc (Kennicutt et al. 1989; Ho et al. 1997). H II regions formed within the stripped gas of some perturbed galaxies have been observed in nearby clusters. Spectroscopic observations of these peculiar regions give quite different estimates of the [NII]λ6583/Hα line ratio, ranging from 0.1 to 1 in the tails of two galaxies in the Coma cluster (Yoshida et al. 2012), to ≃0.4 in ESO 137-001 (Fossati et al. 2016), and ≃0.2 in the perturbed galaxy JO206 recently observed with MUSE by Poggianti et al. (2017). Given that the extraplanar H II regions of NGC 4254 have been formed after the collapse of the stripped gas which has been removed from the outer disc of the galaxy after a rapid encounter with another cluster member, we expect a low-metallicity environment. Indeed, the mean metallicity of NGC 4254 is 12 + log O/H = 8.73 (Hughes et al. 2013), but the galaxy has a steep metallicity gradient suggesting that the metallicity in the outer disc is 12 + log (O/H) ~ 8.5 (Skillman et al. 1996). This metallicity also corresponds to that observed in the outer discs of spiral galaxies (12 + log (O/H) ~ 8.45; Sánchez-Menguiano et al. 2016; Bresolin 2017). We thus assume for the correction [NII]λ6583/Hα = 0.2 (Pettini & Pagel 2004). We also assume a mean dust attenuation of A(H α) = 0.7, a value consistent with those observed in the extraplanar H II regions of the same perturbed galaxies (0 ≲ A(Hα) ≲ 1 mag – Fossati et al. 2016; Poggianti et al. 2017), in the outer discs of late-type galaxies (0 ≲ A(Hα) ≲ 2 mag – Kennicutt et al. 1989; Sánchez et al. 2015), or in UV extended discs (A(Hα) ≃ 0.7 mag – Bresolin et al. 2009, 2012).

Figure 5 shows the distribution of the Hα luminosity of the detected and undetected extraplanar H II regions of NGC 4254. Figure 5 also shows how the distribution would shift assuming different corrections for dust attenuation (A(Hα) = 0.0–1.5) or [NII] contamination ([NII]λ6583/Hα = 0.0–1.5). The luminosity of the H II regions is in the range 1037L(H α) ≲ 1039 (erg s−1), and corresponds thus to the typical luminosity of giant or super-giant Galactic or extragalactic H II regions (Lee et al. 2011).

5.2 Colour analysis

The blue colours of the H II regions seen in Figs. 2 and 3 indicate that they are dominated by very young stellar populations. The fact that they are all detected in the two UV bands, while only a few are in the Hα and u-band, and most undetected at longer wavelengths, suggests that they have been formed by a single coeval, short-lived burst of star formation. To estimate their typical age we first compare their observed (H α–FUV) and (FUV–NUV) age-sensitive colours to the typical synthetic colours of H II regions derived assuming a star formation history defined as: SFR(t)=t×etτMyr1,\begin{equation*} \textrm{SFR}(t) = t \times {\textrm{e}}^{-\frac{t}{\tau}}\,M_{\odot}\,\textrm{yr}^{-1} ,\end{equation*}(2)

the form of which is shown in Fig. 6 (upper panel). The τ e-folding time is set to 3 Myr, consistent with the typical age of giant (L(Hα) ≃ 1037 erg s−1) H II regions in the Milky Way or in other nearby galaxies (Copetti et al. 1985; Tremblin et al. 2014). This approach is similar to the one adopted by Fumagalli et al. (2011) for dating the H II regions of IC3418 and Boquien et al. (2007) for the intergalactic star-formingregions around NGC 5291. We then estimate the typical colour of the synthetic H II regions assuming that the burst of star formation occurred at different epochs (Fig. 6, lower panel). The synthetic colours are derived using the CIGALE SED fitting code (Noll et al. 2009), assuming the Bruzual & Charlot (2003) population synthesis models for the stellar emission with a Salpeter IMF and an updated version of the Draine & Li (2007) physical models of dust emission (see Boselli et al. 2016b for details). To use the Hα emission as a constraint in the fit, we convert the extinction-corrected Hα fluxes into numbers of ionising photons as in Boselli et al. (2016b). We assume an escape fraction of the ionising radiation equal to zero, and that the ionising radiation is not absorbed by dust before ionising the gas (e.g. Boselli et al. 2009). The models are determined assuming different values for the dust attenuation of the stellar continuum (E(BV) = 0.0, 0.1, 0.2), as expected if their metallicity and dust content change from region to region. Geometrical effects on the stellar continuum related to the age of the stellar population are taken into account using the Calzetti et al. (2000) attenuation law. To combine the Hα line flux emission with the broad-band monochromatic flux FUV (in units of mag) in a colour index we define the H α-FUV colour index as: HαFUV=2.5×log(LyC)+20FUV,\begin{equation*} {\textrm{H}\alpha-\textrm{FUV} = -2.5 \times {\textrm{log}}(\textrm{LyC}) + 20 - \textrm{FUV}}, \end{equation*}(3)

where LyC is the flux in the Lyman continuum pseudo filter PSEUDOLyC (in units of μJy) as defined in Eq. (1). Figure 7 shows how the H α-FUV and FUV–NUV colour indices vary as a function of the age of the burst. Clearly the two colour indices are sensitive to the age of the H II regions, the Hα-FUV to ages ≲150 Myr, and the FUV–NUV to ages in the range 500 ≲age ≲ 1000 Myr. Only the former can be used on the present set of data: the FUV–NUV colour index, indeed, gradually increases from FUV–NUV ≃ 0 mag at age = 0 Myr to FUV–NUV ≃ 0.6 mag at age = 1000 Myr, an overly limited variation compared to the observed distribution of the FUV–NUV colour index and its large uncertainty (≃0.2 mag). Figure 7 suggests that the age of the extraplanar H II regions of NGC 4254 is ≲100 Myr. This is also the case for the H II regions undetected in Hα (age ≳ 10 Myr) and in the u-band (age ≲ 200 Myr).

The observed colours of these H II regions can also be compared to the evolutionary tracks derived by CIGALE in different photometric bands sensitive to young stellar populations, Hα, FUV, NUV, and u (Fig. 8). Again, despite the large photometric uncertainties on the data, or the large number of undetected sources in some of the photometric bands, the colours of the H II regions are all consistent with very young ages (0–50 Myr). Figure 8 also shows the presence of several outliers with colours significantly different from those predicted by the evolutionary tracks. A clear example are those with FUV–NUV > 0.4 in the upper panel. The visual inspection of those detected in Hα on the high-resolution MegaCam images revealed the presence of multiple sources unresolved in the GALEX frames. If these single H II regions have different star formation histories, it is plausible that their UV colour in the unresolved GALEX images does not follow the proposed evolutionary tracks.

thumbnail Fig. 5

Upper panel: Hα luminosity distribution of the extraplanar H II regions of NGC 4254. The black shaded histogram indicates the Hα detected regions, the red histogram the distribution of the upper limits. The vertical dashed line shows the detection limit of VESTIGE, the vertical dotted line the median value for the Hα luminosity distribution of the detected sources. The error bars show how the median value of the distribution would shift assuming different corrections for dust attenuation (A(H α = 0.0−1.5) or [NII] contamination ([NII]λ6583/Hα = 0.0−0.5). Lower panel: distribution of the inferred stellar mass. The black open histogram shows the distribution of all the H II regions.

thumbnail Fig. 6

Upperpanel: delayed star formation history with an e-folding time of 3 Myr assumed to fit the observed SED of the extraplanar H II regions of NGC 4254. Lower panel: typical age of their dominant stellar population is the age of the burst (80 Myr – red, 50 Myr – green, 20 Myr – blue).

thumbnail Fig. 7

Left panels: variations of the Hα–FUV (upper) and FUV–NUV (lower) age-sensitive colour indices as a function of time derived for the star formation history given in Eq. (1). Different colours are used for three different dust attenuations: E(BV) = 0.0 (blue), 0.1 (green), and 0.2 (red). The black solid line shows the median of the observed colour distribution, the grey lines the 16% and the 84% quartiles of the distribution. Right panels: distribution of the two colour indices for the detected (black) and undetected (red – lower limit) H II regions.

5.3 SED fitting

We estimate the age of each single H II region using the CIGALE SED fitting code assuming the same configuration as given above, that is, the star formation history given in Eq. (1), the Bruzual & Charlot (2003) population synthesis models derived for a Salpeter IMF, coupled with the Draine & Li (2007) dust models. We fit all (11) photometric bands from the FUV to the Spitzer/MIPS 24 μm with the exception of the r-band where the contamination of the emittingHα line can be dominant combined with the Hα emission using the PSEUDOLyC pseudo filter as in Boselli et al. (2016b), and treating upper limits as in Sawicki (2012). Upper limits in the far-infrared bands are used to constrain the dust attenuation in the star-forming regions. The grid of models is created using the parameters listed in Table 1 (244 620 models).

The output of the fit (star formation rates, stellar masses, ages, and their uncertainties, self-consistently determined assuming the star formation history given in Eq. (1)) for each single H II region is given in Table 2. To quantify how the ages of the regions depend on the assumed dust attenuation and [NII] contamination corrections on the Hα flux we compare in Fig. 9 two sets of values derived assuming in the first one [NII]λ6583/Hα = 0.1 and A(H α) = 0.0, and in the second [NII]λ6583/Hα = 0.2 and A(H α) = 0.7. We also test the reliability of the output parameters and their uncertainties given in Table 2 by creating a photometric mock catalogue introducing into the observed data some extra noise, randomly distributed according to a Gaussian curve of standard deviation equal to the median error for each band, and re-fit the data using the same SED models (see Giovannoli et al. 2011; Boquien et al. 2012; Ciesla et al. 2014 for details). We follow this step twice, the first time keeping the relatively uncertain [NII] contamination and Hα attenuation given in Sect. 5.1 ([NII]λ6583/Hα = 0.2, A(H α) = 0.7), and the second time assuming more extreme values ([NII]λ6583/Hα = 0.1, A(H α) = 0.0); we compare the output of the fit in Fig. 10. Figures 9 and 10 clearly show that, despite the number of parameters used to create the model SEDs, the large uncertainty on the photometric data and that on the corrections on the [NII] contamination and on the Hα attenuation, the stellar masses and the ages of the H II regions derived by the SED fitting code do not vary significantly. The mean age of the stellar population is systematically underestimated whenever A(H α) is overestimated; the opposite is true when [NII]/Hα is overestimated. Star formation rates are instead very uncertain once they are below 10−4 M yr−1. The SED fitting analysis thus confirms the results obtained from the colour analysis done in the previous section, that is, that these H II regions have been formed ≲100 Myr ago, with a few exceptions, some of which might be background galaxies.

thumbnail Fig. 8

Age-sensitive colour–colour diagrams. Upper left panel: Hα–FUV vs. FUV–NUV, lower left panel: NUV–u vs. FUV–NUV, lower right panel: 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, and green open squares regions undetected in Hα and in the u-band. 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(BV) = 0.0 (blue), 0.1 (green), and 0.2 (red), respectively.

Table 1

Input parameters for CIGALE.

5.4 Dynamical age

It has been shown that the size of the H II regions is tightly related to their dynamical age, defined as tdyn ~ rHII∕ΔV, where rHII is the radius of the H II region and ΔV its expansion velocity (Ambrocio-Cruz et al. 2016). We do not have any dynamical data for these H II regions, but thanks to the excellent quality of the Hα image we can measure their angular size and compare their stellar age-size relation to the dynamical age-size relation determined in the LMC by Ambrocio-Cruz et al. (2016). Upon doing this we make the assumption that the dynamical evolution of the extraplanar H II regions of NGC 4254 is similar to that of typical disc H II regions despite the conditions of the surrounding medium (gas density and temperature) probably being significantly different from those generally encountered in the interstellar medium of late-type galaxies (Tonnesen & Bryan 2012). For a fair comparison we need to measure the size of the H II regions as consistently as possible with Ambrocio-Cruz et al. (2016).

We measure the size of the few point-like, symmetric H II regions of the sample (14/60 objects) with a single counterpart within the selected UV region. This last condition is required to avoid H II regions formed by distinct associations of newly formed stars, for which a measure of a diameter might be meaningless. This measurement is performed using the GALFIT code (Peng et al. 2002, 2010) run on the narrow-band Hα image by assuminga Gaussian profile for the H II region and using an empirical PSF image derived from 138 non-saturated stars in the narrow-band image of NGC 4254. We adopt as diameter of the H II region the diameter including 80% of the total flux. When σ is the widthparameter of the GALFIT Gaussian model, the diameter including 80% of the total flux is Diameter = 1.3 ×σ. The diameters of these H II regions are given in Table 2. We then compare the age-size relation of the extraplanar H II of NGC 4524 (where here the age is the one derived from SED fitting) to the dynamical age vs. size relation of H II regions in the LMC in Fig. 113. We also notice that none of these regions are outliers in Fig. 8.

Despite a possible large uncertainty on the measure of the diameter of the H II regions due to the size close to the seeing limit and to a possible blending of multiple star clusters, Figure 11 shows that the age derived from the SED fitting analysis is larger than the dynamical age expected for H II regions of similar size. This suggests that the ages derived from CIGALE are probably upper limits, as is indeed expected given the limits of the stellar population synthesis codes to derive ages younger than a few million years.

Table 2

Parameters derived from the SED fitting.

thumbnail Fig. 9

Relationship between the age determined assuming [NII]λ6583/Hα = 0.2 and A(Hα) = 0.7 (X-axis) and that assuming [NII]λ6583/Hα = 0.1 and A(Hα) = 0.0 (Y -axis). Filled dots indicates the H II regions where the SED fitting gives reduced χr2<6$\chi^2_r < 6$, empty dots χr26$\chi^2_r \geq 6$. The black solid lines show the 1:1 relations.

6 Discussion

6.1 Detection of the stripped material

The first result of the present analysis is a further confirmation that the detection of the material stripped from a perturbed galaxy in rich environments requires multi-frequency observations to trace the different gas, dust, and stellar components. There are, indeed, several examples where the stripped material is mainly cold HI gas not associated with any process of star formation (Chung et al. 2007; Boissier et al. 2012; Scott et al. 2012). This gas can also be in the cold molecular phase whenever giant molecular clouds associated with star-forming regions are present (Verdugo et al. 2015; Jáchym et al. 2017), although this gas phase is not always detected (Jáchym et al. 2013), or in warm H2 (Sivanandam et al. 2014). Some of the tails are also detected as diffuse structures in Hα (ionised gas; Gavazzi et al. 2001; Yoshida et al. 2002; Cortese et al. 2006; Yagi et al. 2007, 2010, 2017; Sun et al. 2007; Kenney et al. 2008; Fossati et al. 2012, 2016; Zhang et al. 2013; Boselli & Gavazzi 2014; Fumagalli et al. 2014), in X-rays (Sun et al. 2006, 2007, 2010), or in radio continuum (Gavazzi & Jaffe 1985, 1987; Gavazzi et al. 1995). Others are detected mainly in UV and optical bands because of the presence of compact stellar regions formed after the interaction, as observed in the wake of IC 3418 in Virgo (Hester et al. 2010; Fumagalli et al. 2011; Kenney et al. 2014) or the jellyfish galaxies of Poggianti et al. (2017).

In the case of NGC 4254, the stripped material has been detected in HI (Phookun et al. 1993; Haynes et al. 2007), while the compact H II regions are clearly visible only in the GALEX UV bands, or in the deep Hα image, but only within ≲20 kpc of the edge of the galaxy disc, while they are totally absent in the outer HI tail. Due to their very young age, they are almost undetected in the deep NGVS optical bands. NGC 4254 would thus be absent in a selection of perturbed galaxies such as the one done in the optical B-band by Poggianti et al. (2017) (GASP survey).

thumbnail Fig. 10

Relationship between the stellar mass (left panel), the star formation rate (centre panel), and the mean age of the stellar population (right panel) derived with CIGALE on the observed data (abscissa) and that derived on the mock catalogue as described in Sect. 3 (ordinate). Black filled dots are for Hα fluxes corrected for [NII] contamination and dust attenuation assuming [NII]λ6583/Hα = 0.2 and A(Hα) = 0.7 (see Sect. 5.1), blue open circles assuming [NII]λ6583/Hα = 0.1 and A(Hα) = 0.0. The black solid lines show the 1:1 relations.

thumbnail Fig. 11

Dynamical age vs. size relation derived for H II regions in the LMC (red symbols, Ambrocio-Cruz et al. 2016) is compared to the age vs. size relation derived from SED fitting for the extraplanar H II regions inNGC 4254 (black filled dots). The vertical dashed line indicates the limiting seeing of the r- and narrow-band images. The black solid line shows the best fit for H II regions in the LMC: log(age) = 0.86(±0.03) × log(Diameter) −1.58(±0.06), R2 = 0.88.

6.2 Impact of a stochastic IMF on the dating of the H II regions

The analysis indicates that the observed regions have Hα luminosities of 1037L(H α) ≲ 1039 erg s−1, which is typical of giant H II regions in the Milky Way or in nearby galaxies, and ages ≲ 100 Myr. Given their luminosity and stellar mass (Mstar ≳ 104 M), these H II regions straddle the boundary of mass at which stochastic sampling of the IMF starts to play a role (Fumagalli et al. 2011; da Silva et al. 2012; Cerviño et al. 2013). At this luminosity, the H II regions should include ~1000 stars (Lee et al. 2011) and quite a few ionising objects (Elmegreen 2000; Koda et al. 2012). Nevertheless, Mstar ≃ 104 M is close to the limit where the undersampling of the IMF starts to be important for the determination of the ionising stellar component using population synthesis models (Cerviño et al. 2003, 2013; Cerviño & Luridiana 2004; da Silva et al. 2012). If the IMF is undersampled, however, the observed Hα luminosity at constant age is on average less than that predicted by our fully sampled population synthesis models. Thus, the ages of the H II regions may be on average overestimated for the smallest H II regions. This effect would not alter our main conclusion that the extraplanar H II regions of NGC 4254 are young (≲100 Myr).

6.3 Star formation process within extraplanar H II regions

IFU observations of the stripped gas of cluster galaxies indicate that star formation only occurs when the velocity dispersion in the tail is sufficiently low (σ ~ 25−50 km s−1; Fumagalli et al. 2014; Fossati et al. 2016; Consolandi et al. 2017) to allow the formation of giant molecular clouds (GMCs). These velocities dispersions are comparable to those observed over the disc of the galaxy by Wilson et al. (2011). ALMA observations are necessary to detect the molecular gas in these compact regions and derive its typical column density and velocity dispersion. Indeed, none of the few H II regions mapped with HERA are detected in the 12CO(2–1) line (Leroy et al. 2009). This is not surprising if the extraplanar star forming regions of NGC 4254 have been formed within GMCs having properties similar to those observed within the Milky Way or in nearby spiral galaxies, for example, masses in the range 5 × 104MGMC ≲ 5 × 106 M, sizes 5 ≲ RGMC ≲ 30 pc, and thus gas column densities Σ (H2) ≃ 170 M pc−2 and volume densities 50 ≲ n(H2) ≲ 500 cm−3 (Solomon et al. 1987; Engargiola et al. 2003). At a distance of 16.5 Mpc the filling factor of these typical GMCs is ≃10%, thus the total CO emission is sufficiently diluted not to be detected within a 13 arcsec beam. Furthermore, the CO emission should also be low because of the expected low metallicity of the gas.

At these column densities, the typical collapse time (free fall time) for a cloud with a spherically symmetric distribution of mass is: tff=3π32Gρ5 Myr,\begin{equation*} {t_{ff} = \sqrt{\frac{3\pi}{32G\rho}} \simeq \rm{5\,Myr}} ,\vspace*{-4pt}\end{equation*}(4)

where ρ is the mean density of the gas. This timescale can be compared with the time since the first interaction of the galaxy with a nearby companion, which dynamical simulations place between 280 Myr (Vollmer et al. 2005) and 750 Myr ago (Duc & Bournaud 2008), and with the typical age of the stellar populations of the H II regions (≲100 Myr). The stripped gas must have reached the typical column density of GMC in ≲40 Myr (see next section) to produce the observed extraplanar H II regions. The densest regions detected at the VLA in the HI tail have column densities of Σ(HI) ≃ 1019 cm−2 (Fig. 4; Chung et al. 2009) and are unresolved at the ~30 arcsec beam resolution. Assuming a filling factor of ~10%, and a clumpy distribution, it is possible that the HI gas reaches column densities similar to those encountered in the outer discs of late-type galaxies producing molecular clouds where star formation takes place. Indeed, star formation seems to follow the typical Schmidt–Kennicutt relation in tidal debries (Boquien et al. 2011). However, the physical conditions of the stripped gas within a hot intergalactic medium (gravitational potential, pressure, temperature, turbulence) change significantly with respect to the typical ISM in galactic discs, and it is still unclear why in these outer regions the gas collapses to form giant H II regions while at further distances from the galaxy along the tail or in other similar objects within the same cluster it does not (e.g. NGC 4569; Boselli et al. 2016a). In particular, in these extraplanar H II regions the lack of a dominant galactic potential well should reduce the hydrostatic pressure of the gas, and thus prevent the formation of molecular gas (Blitz & Rosolowsky 2006). However, turbulence and shock fronts produced by the gravitational perturbation or by the ram pressure that the galaxy is suffering while entering into the cluster can compress the gas and lead to H2 formation (Tonnesen & Bryan 2010, 2012). In NGC 4254 the stripping of the gas seems to occur mostly edge-on. Simulations suggest thatin this geometrical configuration the gas, located in a thin disc before the interaction, is dispersed to a lesser extent than in a face-on stripping process, thus favouring the formation of molecular clouds (Roediger & Brüggen 2006).

Figure 12 shows the relation between the age of the stellar population of the different H II regions and their projected distance from the nucleus of the galaxy. There is no obvious trend between the two quantities, implying that there is no strong differential evolution of these regions after the perturbation. We recall, however, that all these regions are located within ≃35 kpc in projected distance from the galaxy nucleus, corresponding to ≃2 optical radii of the galaxy. This small distance range probably corresponds to a small dynamic range in stellar age which is difficult to resolve with a simple SED fitting analysis, particularly with a sample characterised by a low detection rate at several wavelengths.

6.4 Fate of the stripped material

Considering the line-of-sight motion of NGC 4254 within the cluster we can estimate the time necessary to travel ~30 kpc from the galaxy, the mean projected distance of the H II regions. NGC 4254 is a member of the northern cloud of the Virgo cluster (Gavazzi et al. 1999) and has a relative velocity of ~1400 km s−1 with respect to the cluster centre. Assuming that the galaxy is moving on the plane of the sky with the typical velocity dispersion of the cluster (≃800 km s−1; Boselli et al. 2014b), it would have traveled ~30 kpc in ~40 Myr. In the same time the galaxy would have crossed ~55 kpc along the line of sight. This rough estimate suggests that the furthest H II regions are at RHII ~60 kpc from the nucleus of the galaxy. At this distance these H II regions are still gravitationally bound to NGC 4254. Indeed, the truncation radius RT of the galaxy, which lies at a projected distance of RM87 ~ 1 Mpc from the core of the cluster (M87) is (Binney & Tremaine 2008): RT=RM87×(MNGC4254,DM3×MVirgo,DM)1/3,\begin{equation*} {R_T = R_{\textrm{M87}} \times \left(\frac{M_{\textrm{NGC4254,DM}}}{3 \times M_{\textrm{Virgo,DM}}}\right)^{1/3}}, \end{equation*}(5)

where MNGC4254,DM is the total (dark matter) mass of NGC 4254 and MVirgo,DM that of the Virgo cluster. Assuming MVirgo,DM = 4 × 1014 M (McLaughlin 1999) and MNGC 4254,DM = 1012 M as derived from Behroozi et al. (2010) for a stellar mass of NGC 4254 of MNGC4254,star = 2.4 × 1010 M, the truncation radius of NGC 4254 is RT ≃ 95 kpc. These regions are within the truncation radius of the galaxy and might thus fall back on the disc. However, the HI observations of Phookun et al. (1993) suggest that the HI clouds associated with these regions have a recessional velocity along the line of sight that is ~200 km s−1 lower than that of the galaxy (see Fig. 4), corresponding to ~340 km s−1 in three-dimensional space. This value is close to the escape velocity from a spherically symmetric solid body given by the relation: ve=2GMNGC4254,DMRHII,\begin{equation*} {v_e = \sqrt{\frac{2GM_{\textrm{NGC4254,DM}}}{R_{\textrm{HII}}}}}, \end{equation*}(6)

that for the furthest H II regions (RHII ~ 60 kpc) is ~390 km s−1. Furthermore, the galaxy moving towards the cluster centre is undoubtedly undergoing a strong ram pressure event as clearly seen in the asymmetric distribution of the radio continuum emission (Kantharia et al. 2008) and of the HI gas over its disc, with an extended tail in the northern direction (Phookun et al. 1993; Kantharia et al. 2008). This external pressure might have favoured the formation of GMCs (Blitz & Rosolowsky 2006), although some clouds might have been stripped from the wester spiral arm during the interaction. Here, ram pressure is indeed made more efficient than elsewhere in the disc because of the combined effect of the edge-on motion of the galaxy within the cluster and its clockwise rotation. The HI gas blobs associated with the extraplanar H II regions are suffering the same ram pressure event affecting the galaxy, as indeed suggested by their lower recessional velocity with respect to NGC 4254 (≃200 km s−1, see Fig. 4), preventing them from falling back onto the stellar disc. Given these considerations, while keeping in mind the large uncertainties on all the rough calculations made above, it is possible that some of the extraplanar H II regions will become free-floating objects within the cluster, as indeed indicated by the dynamical simulations of Duc & Bournaud (2008). In that case, these compact stellar objects could contribute to the diffuse intracluster light of Virgo (Mihos et al. 2005, 2009, 2017). Because of their compact nature and their single-age stellar populations, these regions could be the progenitors of globular clusters, diffuse star clusters (DSC), and ultra compact dwarf galaxies (UCD), mainly observed around the dominant massive elliptical galaxies in nearby clusters (Peng et al. 2006; Liu et al. 2015, 2016). Their stellar masses (see Table 2 and Fig. 5) are less than those of typical UCDs and globular clusters, but these latter objects have been formed several Gyr ago, when the typical gas content of galaxies was significantly larger than in similar objects in the local universe. It is thus plausible that high-speed fly-by encounters of gas-rich galaxies in forming clusters at high redshift might have produced more massive compact objects than those observed around NGC 4254, giving birth to globular clusters, DSCs, and UCDs such as those observed in the local universe. The properties of the extraplanar H II regions of NGC 4254 in terms of stellar mass are comparable to those observed in other Virgo cluster compact sources formed after a gravitational perturbation between galaxies (M 49 – VCC 1249; Arrigoni Battaia et al. 2012) or by ram pressure stripping (IC 3418 – Hester et al. 2010; Fumagalli et al. 2011; Kenney et al. 2014). They are remarkably similars to those of the ultra-compact high-velocity cloud AGC 226067, a possible stripping remnant in the Virgo cluster, characterised by a single stellar population with an age of ≃7–50 Myr (Sand et al. 2017).

thumbnail Fig. 12

Left panel: relationship between the age of the different H II regions as derived from CIGALE (in Myr) and their projected distance from the nucleus of NGC 4254 (in kpc). Filled or empty dots indicate the H II regions where the SED fitting gives reduced χr2<6$\chi^2_r < 6$ or χr26$\chi^2_r\,{\geq}\,6$, respectively. Blue big filled dots indicate the median values within different distance bins. The red long-dashed vertical line indicates the23.5 mag arcsec−2 i-band isophotal radius of the galaxy given in Cortese et al. (2012b). Right panel: age distribution of all the H II regions.

7 Conclusion

Very deep Hα images obtained during the VESTIGE survey of the Virgo cluster (f(Hα) ~ 4 × 10−17 erg s−1 cm−2 for point sources at 5σ and Σ (H α) ~ 2 × 10−18 erg s−1 cm−2 arcsec−2 for extended sources at 1σ), combined with UV GALEX and optical NGVS data, reveal the presence of 60 compact (70–500 pc) star-forming regions outside the optical disc of the massive spiral galaxy NGC 4254 (M99). These regions are located along the tail of HI gas harassed from the galaxy after a gravitational interaction with another Virgo cluster member that simulations indicate occurred between 280 and 750 Myr ago (Vollmer et al. 2005; Duc & Bournaud 2008). The analysis of their stellar populations, as well as their physical sizes, consistently indicate that these regions are young (≲100 Myr). These observations offer new observational constraints for models and simulations of star formation in the stripped gas of cluster galaxies. Consistent with other observations of perturbed galaxies, these data indicate that the process of star formation outside the disc of galaxies is very episodic, lasts for a few million years, and occurs only in compact regions. The velocity dispersion of the gas in these compact regions is expected to be sufficiently low to allow matter to collapse and form GMCs able to shield the gas from the external heating due to the hot intracluster medium. A complete multiphase modelling of the star formation process within these extreme regions requires more observational constraints on the different gas phases (HI, H2, X-rays), as well as high-resolution spectroscopy for the determination of their kinematic and physical properties.

Table A.1

Flux densities of the H II regions.

Acknowledgements

We thank the anonymous referee for useful comments on the manuscript. We are grateful to M. Haynes for providing uswith the HI data of NGC 4254, to C. Spengler for the use of GALFIT, and to the whole CFHT team who assisted us in the preparation and in the execution of the observations and in the calibration and data reduction: Todd Burdullis, Daniel Devost, Bill Mahoney, Nadine Manset, Andreea Petric, Simon Prunet, Kanoa Withington. We acknowledge financial support from “ProgrammeNational 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). MB was supported by MINEDUC-UA projects, code ANT 1655 and ANT 1656. MS acknowledges support from the NSF grant 1714764 and the Chandra Award GO6-17111X. MF acknowledges support by the science and technology facilities council [grant number ST/P000541/1].

References

  1. Ambrocio-Cruz, P., Le Coarer, E., Rosado, M., et al. 2016, MNRAS, 457, 2048 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arrigoni Battaia, F., Gavazzi, G., Fumagalli, M., et al. 2012, A&A, 543, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Behroozi, P. S., Conroy, C., & Wechsler, R. H. 2010, ApJ, 717, 379 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bekki, K., Koribalski, B. S., & Kilborn, V. A. 2005, MNRAS, 363, L21 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bellhouse, C., Jaffé, Y. L., Hau, G. K. T., et al. 2017, ApJ, 844, 49 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bendo, G. J., Galliano, F., & Madden, S. C. 2012, MNRAS, 423, 197 [NASA ADS] [CrossRef] [Google Scholar]
  7. Binney, J., & Tremaine, S. 2008, in Galactic Dynamics, 2nd edn. (Princeton: Princeton University Press) [Google Scholar]
  8. Blakeslee, J. P., Jordán, A., Mei, S., et al. 2009, ApJ, 694, 556 [NASA ADS] [CrossRef] [Google Scholar]
  9. Blitz, L., & Rosolowsky, E. 2006, ApJ, 650, 933 [NASA ADS] [CrossRef] [Google Scholar]
  10. Boissier, S., Boselli, A., Duc, P.-A., et al. 2012, A&A, 545, A142 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Boquien, M., Duc, P.-A., Braine, J., et al. 2007, A&A, 467, 93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Boquien, M., Duc, P.-A., Wu, Y., et al. 2009, AJ, 137, 4561 [NASA ADS] [CrossRef] [Google Scholar]
  13. Boquien, M., Duc, P.-A., Galliano, F., et al. 2010, AJ, 140, 2124 [NASA ADS] [CrossRef] [Google Scholar]
  14. Boquien, M., Lisenfeld, U., Duc, P.-A., et al. 2011, A&A, 533, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Boquien, M., Buat, V., Boselli, A., et al. 2012, A&A, 539, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Boquien, M., Buat, V., & Perret, V. 2014, A&A, 571, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517 [NASA ADS] [CrossRef] [Google Scholar]
  18. Boselli, A., & Gavazzi, G. 2014, A&ARv, 22, 74 [NASA ADS] [CrossRef] [Google Scholar]
  19. Boselli, A., Boissier, S., Cortese, L., et al. 2006, ApJ, 651, 811 [NASA ADS] [CrossRef] [Google Scholar]
  20. Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008a, ApJ, 674, 742 [NASA ADS] [CrossRef] [Google Scholar]
  21. Boselli, A., Boissier, S., Cortese, L., & Gavazzi, G. 2008b, A&A, 489, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Boselli, A., Boissier, S., Cortese, L., et al. 2009, ApJ, 706, 1527 [NASA ADS] [CrossRef] [Google Scholar]
  23. Boselli, A., Eales, S., Cortese, L., et al. 2010, PASP, 122, 261 [NASA ADS] [CrossRef] [Google Scholar]
  24. Boselli, A., Boissier, S., Heinis, S., et al. 2011, A&A, 528, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Boselli, A., Cortese, L., Boquien, M., et al. 2014a, A&A, 564, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Boselli, A., Voyer, E., Boissier, S., et al. 2014b, A&A, 570, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Boselli, A., Cuillandre, J. C., Fossati, M., et al. 2016a, A&A, 587, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Boselli, A., Roehlly, Y., Fossati, M., et al. 2016b, A&A, 596, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Boselli, A., Fossati, M., Ferrarese, L., et al. 2018, A&A, 614, A56 (Paper I) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Bresolin, F. 2017, Outskirts of Galaxies, 434, 145 [NASA ADS] [CrossRef] [Google Scholar]
  31. Bresolin, F., Ryan-Weber, E., Kennicutt, R. C., & Goddard, Q. 2009, ApJ, 695, 580 [NASA ADS] [CrossRef] [Google Scholar]
  32. Bresolin, F., Kennicutt, R. C., & Ryan-Weber, E. 2012, ApJ, 750, 122 [NASA ADS] [CrossRef] [Google Scholar]
  33. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  34. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  35. Cayatte, V., van Gorkom, J. H., Balkowski, C., & Kotanyi, C. 1990, AJ, 100, 604 [NASA ADS] [CrossRef] [Google Scholar]
  36. Cayatte, V., Kotanyi, C., Balkowski, C., & van Gorkom, J. H. 1994, AJ, 107, 1003 [NASA ADS] [CrossRef] [Google Scholar]
  37. Cerviño, M., & Luridiana, V. 2004, A&A, 413, 145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Cerviño, M., Luridiana, V., Pérez, E., Vílchez, J. M., & Valls-Gabaud, D. 2003, A&A, 407, 177 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Cerviño, M., Román-Zúñiga, C., Luridiana, V., et al. 2013, A&A, 553, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Chung, A., van Gorkom, J. H., Kenney, J. D. P., & Vollmer, B. 2007, ApJ, 659, L115 [NASA ADS] [CrossRef] [Google Scholar]
  41. Chung, A., van Gorkom, J. H., Kenney, J. D. P., Crowl, H., & Vollmer, B. 2009, AJ, 138, 1741 [NASA ADS] [CrossRef] [Google Scholar]
  42. Ciesla, L., Boselli, A., Smith, M. W. L., et al. 2012, A&A, 543, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Ciesla, L., Boquien, M., Boselli, A., et al. 2014, A&A, 565, A128 [Google Scholar]
  44. Consolandi G., Gavazzi G., & Fossati M., et al. 2017, A&A, 606, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Copetti, M. V. F., Pastoriza, M. G., & Dottori, H. A. 1985, A&A, 152, 427 [NASA ADS] [Google Scholar]
  46. Cortese, L., Gavazzi, G., Boselli, A., et al. 2006, A&A, 453, 847 [Google Scholar]
  47. Cortese, L., Davies, J. I., Pohlen, M., et al. 2010, A&A, 518, L49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Cortese, L., Ciesla, L., Boselli, A., et al. 2012a, A&A, 540, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Cortese, L., Boissier, S., Boselli, A., et al. 2012b, A&A, 544, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Cortese, L., Fritz, J., Bianchi, S., et al. 2014, MNRAS, 440, 942 [NASA ADS] [CrossRef] [Google Scholar]
  51. Cowie, L. L., & Songaila, A. 1977, Nature, 266, 501 [Google Scholar]
  52. da Silva, R. L., Fumagalli, M., & Krumholz, M. 2012, ApJ, 745, 145 [NASA ADS] [CrossRef] [Google Scholar]
  53. Davies, J., Minchin, R., Sabatini, S., et al. 2004, MNRAS, 349, 922 [NASA ADS] [CrossRef] [Google Scholar]
  54. Davies, J. I., Baes, M., Bendo, G. J., et al. 2010, A&A, 518, L48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [NASA ADS] [CrossRef] [Google Scholar]
  56. Dressler, A. 1980, ApJ, 236, 351 [NASA ADS] [CrossRef] [Google Scholar]
  57. Dressler, A. 2004, Clusters of Galaxies: Probes of Cosmological Structure and Galaxy Evolution, 206 [Google Scholar]
  58. Duc, P.-A., & Bournaud, F. 2008, ApJ, 673, 787 [NASA ADS] [CrossRef] [Google Scholar]
  59. Elmegreen, B. G. 2000, ApJ, 539, 342 [NASA ADS] [CrossRef] [Google Scholar]
  60. Engargiola, G., Plambeck, R. L., Rosolowsky, E., & Blitz, L. 2003, ApJS, 149, 343 [NASA ADS] [CrossRef] [Google Scholar]
  61. Ferrarese, L., Côté, P., Cuillandre, J.-C., et al. 2012, ApJS, 200, 4 [NASA ADS] [CrossRef] [Google Scholar]
  62. Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
  63. Fossati, M., Gavazzi, G., Boselli, A., & Fumagalli, M. 2012, A&A, 544, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Fossati, M., Fumagalli, M., Boselli, A., et al. 2016, MNRAS, 455, 2028 [Google Scholar]
  65. Fossati, M., Mendel, J.T., Boselli, A., et al. 2018, A&A, 614, A57 (Paper II) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Fritz, J., Moretti, A., Gullieuszik, M., et al. 2017, ApJ, 848, 132 [NASA ADS] [CrossRef] [Google Scholar]
  67. Fumagalli, M., Krumholz, M. R., Prochaska, J. X., Gavazzi, G., & Boselli, A. 2009, ApJ, 697, 1811 [NASA ADS] [CrossRef] [Google Scholar]
  68. Fumagalli, M., Gavazzi, G., Scaramella, R., & Franzetti, P. 2011a, A&A, 528, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Fumagalli, M., da Silva, R. L., & Krumholz, M. R. 2011b, ApJ, 741, L26 [NASA ADS] [CrossRef] [Google Scholar]
  70. Fumagalli, M., Fossati, M., Hau, G. K. T., et al. 2014, MNRAS, 445, 4335 [Google Scholar]
  71. Gavazzi, G., & Jaffe, W. 1985, ApJ, 294, L89 [NASA ADS] [CrossRef] [Google Scholar]
  72. Gavazzi, G., & Jaffe, W. 1987, A&A, 186, L1 [NASA ADS] [Google Scholar]
  73. Gavazzi, G., Contursi, A., Carrasco, L., et al. 1995, A&A, 304, 325 [NASA ADS] [Google Scholar]
  74. Gavazzi, G., Catinella, B., Carrasco, L., Boselli, A., & Contursi, A. 1998, AJ, 115, 1745 [NASA ADS] [CrossRef] [Google Scholar]
  75. Gavazzi, G., Boselli, A., Scodeggio, M., Pierini, D., & Belsole, E. 1999, MNRAS, 304, 595 [NASA ADS] [CrossRef] [Google Scholar]
  76. Gavazzi, G., Boselli, A., Mayer, L., et al. 2001, ApJ, 563, L23 [NASA ADS] [CrossRef] [Google Scholar]
  77. Gavazzi, G., Boselli, A., Donati, A., Franzetti, P., & Scodeggio, M. 2003, A&A, 400, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Gavazzi, G., Boselli, A., van Driel, W., & O’Neil, K. 2005, A&A, 429, 439 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Gavazzi, G., Fumagalli, M., Cucciati, O., & Boselli, A. 2010, A&A, 517, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Gavazzi, G., Fumagalli, M., Fossati, M., et al. 2013, A&A, 553, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Giovanelli, R., Haynes, M. P., Kent, B. R., et al. 2005, AJ, 130, 2598 [NASA ADS] [CrossRef] [Google Scholar]
  82. Giovannoli, E., Buat, V., Noll, S., Burgarella, D., & Magnelli, B. 2011, A&A, 525, A150 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Gómez, P. L., Nichol, R. C., Miller, C. J., et al. 2003, ApJ, 584, 210 [NASA ADS] [CrossRef] [Google Scholar]
  84. Gunn, J. E., & Gott, J. R., III 1972, ApJ, 176, 1 [Google Scholar]
  85. Gwyn, S. D. J. 2008, PASP, 120, 212 [NASA ADS] [CrossRef] [Google Scholar]
  86. Haynes, M. P., & Giovanelli, R. 1984, AJ, 89, 758 [NASA ADS] [CrossRef] [Google Scholar]
  87. Haynes, M. P., Giovanelli, R., & Kent, B. R. 2007, ApJ, 665, L19 [NASA ADS] [CrossRef] [Google Scholar]
  88. Hester, J. A., Seibert, M., Neill, J. D., et al. 2010, ApJ, 716, L14 [Google Scholar]
  89. Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJ, 487, 579 [NASA ADS] [CrossRef] [Google Scholar]
  90. Hughes, T. M., Cortese, L., Boselli, A., Gavazzi, G., & Davies, J. I. 2013, A&A, 550, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Jáchym, P., Kenney, J. D. P., Ržuička, A., et al. 2013, A&A, 556, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Jáchym, P., Combes, F., Cortese, L., Sun, M., & Kenney, J. D. P. 2014, ApJ, 792, 11 [NASA ADS] [CrossRef] [Google Scholar]
  93. Jáchym, P., Sun, M., Kenney, J. D. P., et al. 2017, ApJ, 839, 114 [NASA ADS] [CrossRef] [Google Scholar]
  94. Kantharia, N. G., Rao, A. P., & Sirothia, S. K. 2008, MNRAS, 383, 173 [NASA ADS] [CrossRef] [Google Scholar]
  95. Kapferer, W., Sluka, C., Schindler, S., Ferrari, C., & Ziegler, B. 2009, A&A, 499, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Kenney, J. D. P., van Gorkom, J. H., & Vollmer, B. 2004, AJ, 127, 3361 [NASA ADS] [CrossRef] [Google Scholar]
  97. Kenney, J. D. P., Tal, T., Crowl, H. H., Feldmeier, J., & Jacoby, G. H. 2008, ApJ, 687, L69 [Google Scholar]
  98. Kenney, J. D. P., Geha, M., Jáchym, P., et al. 2014, ApJ, 780, 119 [NASA ADS] [CrossRef] [Google Scholar]
  99. Kennicutt, R. C., Jr. 1983, AJ, 88, 483 [NASA ADS] [CrossRef] [Google Scholar]
  100. Kennicutt, R. C., Jr., Keel, W. C., & Blaha, C. A. 1989, AJ, 97, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  101. Kennicutt, R. C., Jr., Armus, L., Bendo, G., et al. 2003, PASP, 115, 928 [NASA ADS] [CrossRef] [Google Scholar]
  102. Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 2001, ApJ, 556, 121 [Google Scholar]
  103. Koda, J., Yagi, M., Boissier, S., et al. 2012, ApJ, 749, 20 [NASA ADS] [CrossRef] [Google Scholar]
  104. Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 [NASA ADS] [CrossRef] [Google Scholar]
  105. Lee, J. H., Hwang, N., & Lee, M. G. 2011, ApJ, 735, 75 [NASA ADS] [CrossRef] [Google Scholar]
  106. Leroy, A. K., Walter, F., Bigiel, F., et al. 2009, AJ, 137, 4670 [NASA ADS] [CrossRef] [Google Scholar]
  107. Lisenfeld, U., Braine, J., Duc, P. A., et al. 2016, A&A, 590, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  108. Liu, C., Peng, E. W., Côté, P., et al. 2015, ApJ, 812, 34 [NASA ADS] [CrossRef] [Google Scholar]
  109. Liu, Y., Peng, E. W., Lim, S., et al. 2016, ApJ, 830, 99 [NASA ADS] [CrossRef] [Google Scholar]
  110. McCall, M. L., Rybski, P. M., & Shields, G. A. 1985, ApJS, 57, 1 [NASA ADS] [CrossRef] [Google Scholar]
  111. McLaughlin, D. E. 1999, ApJ, 512, L9 [NASA ADS] [CrossRef] [Google Scholar]
  112. Mei, S., Blakeslee, J. P., Côté, P., et al. 2007, ApJ, 655, 144 [NASA ADS] [CrossRef] [Google Scholar]
  113. Mihos, J. C., Harding, P., Feldmeier, J., & Morrison, H. 2005, ApJ, 631, L41 [NASA ADS] [CrossRef] [Google Scholar]
  114. Mihos, J. C., Janowiecki, S., Feldmeier, J. J., Harding, P., & Morrison, H. 2009, ApJ, 698, 1879 [NASA ADS] [CrossRef] [Google Scholar]
  115. Mihos, J. C., Harding, P., Feldmeier, J. J., et al. 2017, ApJ, 834, 16 [NASA ADS] [CrossRef] [Google Scholar]
  116. Minchin, R., Davies, J., Disney, M., et al. 2005, ApJ, 622, L21 [NASA ADS] [CrossRef] [Google Scholar]
  117. Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 [NASA ADS] [CrossRef] [Google Scholar]
  118. Morrissey, P., Conrow, T., Barlow, T. A., et al. 2007, ApJS, 173, 682 [Google Scholar]
  119. Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 507, 1793 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2002, AJ, 124, 266 [Google Scholar]
  121. Peng, E. W., Côté, P., Jordán, A., et al. 2006, ApJ, 639, 838 [NASA ADS] [CrossRef] [Google Scholar]
  122. Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097 [NASA ADS] [CrossRef] [Google Scholar]
  123. Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  124. Phookun, B., Vogel, S. N., & Mundy, L. G. 1993, ApJ, 418, 113 [NASA ADS] [CrossRef] [Google Scholar]
  125. Poggianti, B. M., Moretti, A., Gullieuszik, M., et al. 2017, ApJ, 844, 48 [Google Scholar]
  126. Roediger, E., & Brüggen, M. 2006, MNRAS, 369, 567 [NASA ADS] [CrossRef] [Google Scholar]
  127. Roediger, E., & Brüggen, M. 2007, MNRAS, 380, 1399 [NASA ADS] [CrossRef] [Google Scholar]
  128. Roediger, E., & Brüggen, M. 2008, MNRAS, 388, L89 [NASA ADS] [CrossRef] [Google Scholar]
  129. Roediger, E., & Hensler, G. 2005, A&A, 433, 875 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  130. Sánchez, S. F., Pérez, E., Rosales-Ortega, F. F., et al. 2015, A&A, 574, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  131. Sánchez-Menguiano, L., Sánchez, S. F., Pérez, I., et al. 2016, A&A, 587, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  132. Sand, D. J., Seth, A. C., Crnojević, D., et al. 2017, ApJ, 843, 134 [NASA ADS] [CrossRef] [Google Scholar]
  133. Sawicki, M. 2012, PASP, 124, 1208 [Google Scholar]
  134. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  135. Scott, T. C., Cortese, L., Brinks, E., et al. 2012, MNRAS, 419, L19 [NASA ADS] [CrossRef] [Google Scholar]
  136. Sivanandam, S., Rieke, M. J., & Rieke, G. H. 2014, ApJ, 796, 89 [NASA ADS] [CrossRef] [Google Scholar]
  137. Skillman, E. D., Kennicutt, R. C., Jr., Shields, G. A., & Zaritsky, D. 1996, ApJ, 462, 147 [NASA ADS] [CrossRef] [Google Scholar]
  138. Solanes, J. M., Manrique, A., García-Gómez, C., et al. 2001, ApJ, 548, 97 [NASA ADS] [CrossRef] [Google Scholar]
  139. Solomon, P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 319, 730 [NASA ADS] [CrossRef] [Google Scholar]
  140. Sun, M., Jones, C., Forman, W., et al. 2006, ApJ, 637, L81 [NASA ADS] [CrossRef] [Google Scholar]
  141. Sun, M., Donahue, M., & Voit, G. M. 2007, ApJ, 671, 190 [Google Scholar]
  142. Sun, M., Donahue, M., Roediger, E., et al. 2010, ApJ, 708, 946 [NASA ADS] [CrossRef] [Google Scholar]
  143. Thilker, D. A., Bianchi, L., Meurer, G., et al. 2007, ApJS, 173, 538 [NASA ADS] [CrossRef] [Google Scholar]
  144. Tonnesen, S., & Bryan, G. L. 2009, ApJ, 694, 789 [NASA ADS] [CrossRef] [Google Scholar]
  145. Tonnesen, S., & Bryan, G. L. 2010, ApJ, 709, 1203 [NASA ADS] [CrossRef] [Google Scholar]
  146. Tonnesen, S., & Bryan, G. L. 2012, MNRAS, 422, 1609 [NASA ADS] [CrossRef] [Google Scholar]
  147. Tremblin, P., Anderson, L. D., Didelon, P., et al. 2014, A&A, 568, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  148. Verdugo, C., Combes, F., Dasyra, K., Salomé, P., & Braine, J. 2015, A&A, 582, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  149. Vollmer, B., Huchtmeier, W., & van Driel, W. 2005, A&A, 439, 921 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  150. Vollmer, B., Soida, M., Otmianowska-Mazur, K., et al. 2006, A&A, 453, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  151. Vollmer, B., Soida, M., Chung, A., et al. 2008a, A&A, 483, 89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  152. Vollmer, B., Braine, J., Pappalardo, C., & Hily-Blant, P. 2008b, A&A, 491, 455 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  153. Vollmer, B., Soida, M., Chung, A., et al. 2009, A&A, 496, 669 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  154. Vollmer, B., Soida, M., Braine, J., et al. 2012, A&A, 537, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  155. Voyer, E. N., Boselli, A., Boissier, S., et al. 2014, A&A, 569, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  156. Xu, C. K., Donas, J., Arnouts, S., et al. 2005, ApJ, 619, L11 [NASA ADS] [CrossRef] [Google Scholar]
  157. Yagi, M., Komiyama, Y., Yoshida, M., et al. 2007, ApJ, 660, 1209 [NASA ADS] [CrossRef] [Google Scholar]
  158. Yagi, M., Yoshida, M., Komiyama, Y., et al. 2010, AJ, 140, 1814 [Google Scholar]
  159. Yagi, M., Gu, L., Fujita, Y., et al. 2013, ApJ, 778, 91 [NASA ADS] [CrossRef] [Google Scholar]
  160. Yagi, M., Yoshida, M., Gavazzi, G., et al. 2017, ApJ, 839, 65 [NASA ADS] [CrossRef] [Google Scholar]
  161. Yoshida, M., Yagi, M., Okamura, S., et al. 2002, ApJ, 567, 118 [NASA ADS] [CrossRef] [Google Scholar]
  162. Yoshida, M., Yagi, M., Komiyama, Y., et al. 2012, ApJ, 749, 43 [NASA ADS] [CrossRef] [Google Scholar]
  163. Whitmore, B. C., Gilmore, D. M., & Jones, C. 1993, ApJ, 407, 489 [NASA ADS] [CrossRef] [Google Scholar]
  164. Wilson, C. D., Warren, B. E., Irwin, J., et al. 2011, MNRAS, 410, 1409 [NASA ADS] [Google Scholar]
  165. Wong, O. I., Ryan-Weber, E. V., Garcia-Appadoo, D. A., et al. 2006, MNRAS, 371, 1855 [NASA ADS] [CrossRef] [Google Scholar]
  166. Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [Google Scholar]
  167. Zhang, B., Sun, M., Ji, L., et al. 2013, ApJ, 777, 122 [NASA ADS] [CrossRef] [Google Scholar]

1

Hereafter we refer to the Hα + [NII] band simply as Hα, unless otherwise stated.

2

An aperture correction has been applied to the only region detected at 24 μm.

3

The diameters of the H II regions in the LMC are isophotal. Given the small size of the H II region on the image, and the unknown surface-brightness limit used in Ambrocio-Cruz et al. (2016), we decided to use the Diameter definition given above as a representative measure of the size of these extraplanar H II regions. We caution that part of the systematic misalignment observed in Fig. 11 can be due to the different definition of the size used for the two sets of data.

All Tables

Table 1

Input parameters for CIGALE.

Table 2

Parameters derived from the SED fitting.

Table A.1

Flux densities of the H II regions.

All Figures

thumbnail Fig. 1

Colour ugi RGB image of the galaxy NGC 4254 obtained using the NGVS data (Ferrarese et al. 2012). At the distance of the galaxy (16.5 Mpc), 1 arcmin = 4.8 kpc. The star-forming regions analysed in this work are the blue blobs in the south-west outside the disc of the galaxy.

In the text
thumbnail Fig. 2

Multifrequency images of the galaxy NGC 4254 (north is up, east is left). Upper panels, from left to right: continuum-subtracted Hα (VESTIGE), FUV (GALEX), NUV (GALEX); lower panels: u-band (NGVS), IRAC 3.6 μm (Spitzer), MIPS 24 μm (Spitzer). The extraplanar star-forming regions marked with green circles are evident in FUV image at the south-west of the galaxy. The linear size of each single image corresponds 40 × 35 kpc2 at the distance of the galaxy.

In the text
thumbnail Fig. 3

FUV/NUV GALEX colour (left) and continuum-subtracted Hα (VESTIGE; right) magnified images of the extraplanar star-forming regions (marked with green circles) in the south-west quadrant of NGC 4254. The UV colour image shows that these regions have very blue colours, and are thus dominated by young stellar populations (≲100 Myr). Half of them, however, are undetected in Hα, indicating that their typical age is 10 Myr ≲ age ≲ 100 Myr.

In the text
thumbnail Fig. 4

VLA HI contours at a column density Σ(HI) = 1019 cm−2 (from Chung et al. 2009; red) and the Arecibo HI contours at a column density Σ (HI) = 1−2 × 1019 cm−2 (from Haynes et al. 2007; white) are superposed on the GALEX FUV∕NUV colour image. The right panel is a magnified view of the boxed region marked in the left panel with a dashed line. The typical velocity of the different HI structures detected at the VLA are given. The yellow and cyan ellipses at the lower left corner of the two images show the beam sizes at Arecibo (3.3′ × 3.8′) and at the VLA (30.86′′ × 28.07′′).

In the text
thumbnail Fig. 5

Upper panel: Hα luminosity distribution of the extraplanar H II regions of NGC 4254. The black shaded histogram indicates the Hα detected regions, the red histogram the distribution of the upper limits. The vertical dashed line shows the detection limit of VESTIGE, the vertical dotted line the median value for the Hα luminosity distribution of the detected sources. The error bars show how the median value of the distribution would shift assuming different corrections for dust attenuation (A(H α = 0.0−1.5) or [NII] contamination ([NII]λ6583/Hα = 0.0−0.5). Lower panel: distribution of the inferred stellar mass. The black open histogram shows the distribution of all the H II regions.

In the text
thumbnail Fig. 6

Upperpanel: delayed star formation history with an e-folding time of 3 Myr assumed to fit the observed SED of the extraplanar H II regions of NGC 4254. Lower panel: typical age of their dominant stellar population is the age of the burst (80 Myr – red, 50 Myr – green, 20 Myr – blue).

In the text
thumbnail Fig. 7

Left panels: variations of the Hα–FUV (upper) and FUV–NUV (lower) age-sensitive colour indices as a function of time derived for the star formation history given in Eq. (1). Different colours are used for three different dust attenuations: E(BV) = 0.0 (blue), 0.1 (green), and 0.2 (red). The black solid line shows the median of the observed colour distribution, the grey lines the 16% and the 84% quartiles of the distribution. Right panels: distribution of the two colour indices for the detected (black) and undetected (red – lower limit) H II regions.

In the text
thumbnail Fig. 8

Age-sensitive colour–colour diagrams. Upper left panel: Hα–FUV vs. FUV–NUV, lower left panel: NUV–u vs. FUV–NUV, lower right panel: 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, and green open squares regions undetected in Hα and in the u-band. 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(BV) = 0.0 (blue), 0.1 (green), and 0.2 (red), respectively.

In the text
thumbnail Fig. 9

Relationship between the age determined assuming [NII]λ6583/Hα = 0.2 and A(Hα) = 0.7 (X-axis) and that assuming [NII]λ6583/Hα = 0.1 and A(Hα) = 0.0 (Y -axis). Filled dots indicates the H II regions where the SED fitting gives reduced χr2<6$\chi^2_r < 6$, empty dots χr26$\chi^2_r \geq 6$. The black solid lines show the 1:1 relations.

In the text
thumbnail Fig. 10

Relationship between the stellar mass (left panel), the star formation rate (centre panel), and the mean age of the stellar population (right panel) derived with CIGALE on the observed data (abscissa) and that derived on the mock catalogue as described in Sect. 3 (ordinate). Black filled dots are for Hα fluxes corrected for [NII] contamination and dust attenuation assuming [NII]λ6583/Hα = 0.2 and A(Hα) = 0.7 (see Sect. 5.1), blue open circles assuming [NII]λ6583/Hα = 0.1 and A(Hα) = 0.0. The black solid lines show the 1:1 relations.

In the text
thumbnail Fig. 11

Dynamical age vs. size relation derived for H II regions in the LMC (red symbols, Ambrocio-Cruz et al. 2016) is compared to the age vs. size relation derived from SED fitting for the extraplanar H II regions inNGC 4254 (black filled dots). The vertical dashed line indicates the limiting seeing of the r- and narrow-band images. The black solid line shows the best fit for H II regions in the LMC: log(age) = 0.86(±0.03) × log(Diameter) −1.58(±0.06), R2 = 0.88.

In the text
thumbnail Fig. 12

Left panel: relationship between the age of the different H II regions as derived from CIGALE (in Myr) and their projected distance from the nucleus of NGC 4254 (in kpc). Filled or empty dots indicate the H II regions where the SED fitting gives reduced χr2<6$\chi^2_r < 6$ or χr26$\chi^2_r\,{\geq}\,6$, respectively. Blue big filled dots indicate the median values within different distance bins. The red long-dashed vertical line indicates the23.5 mag arcsec−2 i-band isophotal radius of the galaxy given in Cortese et al. (2012b). Right panel: age distribution of all the H II regions.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.