The global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale.

We study the intrinsic shapes of a sample of over 400 quiescent galaxies in the cores of the Virgo and Fornax clusters with luminosities $10^{6} \leq L_{g}/L_{\odot} \leq 10^{8}$. Similar to satellites of the Local Group and Centaurus A, these faint, low surface brightness cluster galaxies are best described as a family of thick ($C/A > 0.5$), oblate-triaxial spheroids. However, the large sample size allows us to show that the flattening of their stellar distributions depends both on luminosity and on the presence of a nuclear star cluster. Nucleated satellites are thicker at all luminosities compared to their non-nucleated counterparts, and fainter galaxies are systematically thicker as well, regardless of nucleation. Once nucleation is accounted for, we find no evidence that the environment the satellites live in plays a relevant role in setting their three-dimensional structure. We interpret both the presence of stellar nuclei and the associated thicker shapes as the result of preferential early and rapid formation, effectively making these faint nucleated galaxies the first generation of cluster satellites.

We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised.

Using deep, high resolution optical imaging from the Next Generation Virgo Cluster Survey we study the properties of nuclear star clusters (NSCs) in a sample of nearly 400 quiescent galaxies in the core of Virgo with stellar masses $10^{5}\lesssim M_{*}/M_{\odot} \lesssim10^{12}$. The nucleation fraction reaches a peak value $f_{n}\approx90\%$ for $M_{*} \approx 10^{9} M_{\odot}$ galaxies and declines for both higher and lower masses, but nuclei populate galaxies as small as $M_{*} \approx5\times10^{5} M_{\odot}$. Comparison with literature data for nearby groups and clusters shows that at the low-mass end nucleation is more frequent in denser environments. The NSC mass function peaks at $M_{NSC}\approx7\times10^{5} M_{\odot}$, a factor 3-4 times larger than the turnover mass for globular clusters (GCs). We find a nonlinear relation between the stellar masses of NSCs and of their host galaxies, with a mean nucleus-to-galaxy mass ratio that drops to $M_{NSC}/M_{*}\approx3.6\times10^{-3}$ for $M_{*} \approx 5\times10^{9} M_{\odot}$ galaxies. Nuclei in both more and less massive galaxies are much more prominent: $M_{NSC}\propto M_{*}^{0.46}$ at the low-mass end, where nuclei are nearly 50% as massive as their hosts. We measure an intrinsic scatter in NSC masses at fixed galaxy stellar mass of 0.4 dex, which we interpret as evidence that the process of NSC growth is significantly stochastic. At low galaxy masses we find a close connection between NSCs and GC systems, including a very similar occupation distribution and comparable total masses. We discuss these results in the context of current dissipative and dissipationless models of NSC formation.

We have observed the late-type peculiar galaxy NGC 4424 during VESTIGE, a blind narrow-band Halpha[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT. The presence of a 110 kpc long HI tail in the S direction indicates that this galaxy is undergoing a ram pressure stripping event. The deep narrow-band image revealed the presence of a low surface brightness ionised gas tail ~10 kpc long extending from the centre of the galaxy to the NW direction, in the direction opposite to the HI tail. Chandra and XMM X-rays data do not show any compact source in the nucleus nor the presence of an extended tail of hot gas, while MUSE spectroscopy indicates that the gas is photo-ionised in the inner regions and shock-ionised in the outer parts. IFU spectroscopy confirms that the ionised gas is kinematically decoupled from the stellar component and indicates the presence of two kinematically distinct structures in the stellar disc. The analysis of the SED of the galaxy indicates that the activity of star formation has been totally quenched in the outer disc ~ 250-280 Myr ago, while only reduced by ~80% in the central regions. All this observational evidence suggests that NGC 4424 is the remnant of an unequal-mass merger occurred <= 500 Myr ago, when the galaxy was already a member of the Virgo cluster, now undergoing a ram pressure stripping event which has removed the gas and quenched the activity of star formation in the outer disc. The tail of ionised gas probably results from the outflow produced by a central starburst fed by the collapse of gas induced by the merging episode. This outflow is sufficiently powerful to overcome the ram pressure induced by the intracluster medium on the disc of the galaxy crossing the cluster. This analysis thus suggests that feedback can participate in the quenching process of galaxies in high-density regions.

We derive stellar population parameters for a representative sample of ultracompact dwarfs (UCDs) and a large sample of massive globular clusters (GCs) with stellar masses $\gtrsim$ 10$^{6}$ $M_{\odot}$ in the central galaxy M87 of the Virgo galaxy cluster, based on model fitting to the Lick-index measurements from both the literature and new observations. After necessary spectral stacking of the relatively faint objects in our initial sample of 40 UCDs and 118 GCs, we obtain 30 sets of Lick-index measurements for UCDs and 80 for GCs. The M87 UCDs have ages $\gtrsim$ 8 Gyr and [$\alpha$/Fe] $\simeq$ 0.4 dex, in agreement with previous studies based on smaller samples. The literature UCDs, located in lower-density environments than M87, extend to younger ages and smaller [$\alpha$/Fe] (at given metallicities) than M87 UCDs, resembling the environmental dependence of the Virgo dE nuclei. The UCDs exhibit a positive mass-metallicity relation (MZR), which flattens and connects compact ellipticals at stellar masses $\gtrsim$ 10$^{8}$ $M_{\odot}$. The Virgo dE nuclei largely follow the average MZR of UCDs, whereas most of the M87 GCs are offset towards higher metallicities for given stellar masses. The difference between the mass-metallicity distributions of UCDs and GCs may be qualitatively understood as a result of their different physical sizes at birth in a self-enrichment scenario or of galactic nuclear cluster star formation efficiency being relatively low in a tidal stripping scenario for UCD formation. The existing observations provide the necessary but not sufficient evidence for tidally stripped dE nuclei being the dominant contributors to the M87 UCDs.

During pilot observations of the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the CFHT, we have observed the spiral galaxy NGC 4254 (M99). Deep Halpha+[NII] narrow-band and GALEX UV images revealed the presence of 60 compact (70-500 pc radius) star forming regions up to ~ 20 kpc outside the optical disc of the galaxy. These regions are located along a tail of HI gas stripped from the disc of the galaxy after a rapid gravitational encounter with another Virgo cluster member that simulations indicate occurred 280-750 Myr ago. We have combined the VESTIGE data with multifrequency data from the UV to the far-infrared to characterise the stellar populations of these regions and study the star formation process in an extreme environment such as the tails of stripped gas embedded in the hot intracluster medium. The colour, spectral energy distribution (SED), and linear size consistently indicate that these regions are coeval and have been formed after a single burst of star formation that occurred ~< 100 Myr ago. These regions might become free floating objects within the cluster potential well, and be the local analogues of compact sources produced after the interaction of gas-rich systems that occurred during the early formation of clusters.

Substructure in globular cluster (GC) populations around large galaxies is expected in galaxy formation scenarios that involve accretion or merger events, and it has been searched for using direct associations between GCs and structure in the diffuse galaxy light, or with GC kinematics. Here, we present a search for candidate substructures in the GC population around the Virgo cD galaxy M87 through the analysis of the spatial distribution of the GC colors.~The study is based on a sample of $\sim\!1800$ bright GCs with high-quality $u,g,r,i,z,K_s$ photometry, selected to ensure a low contamination by foreground stars or background galaxies.~The spectral energy distributions of the GCs are associated with formal estimates of age and metallicity, which are representative of its position in a 4-D color-space relative to standard single stellar population models.~Dividing the sample into broad bins based on the relative formal ages, we observe inhomogeneities which reveal signatures of GC substructures.~The most significant of these is a spatial overdensity of GCs with relatively young age labels, of diameter $\sim\!0.1$\,deg ($\sim\!30\,$kpc), located to the south of M87.~The significance of this detection is larger than about 5$\sigma$ after accounting for estimates of random and systematic errors.~Surprisingly, no large Virgo galaxy is present in this area, that could potentially host these GCs.~But candidate substructures in the M87 halo with equally elusive hosts have been described based on kinematic studies in the past.~The number of GC spectra available around M87 is currently insufficient to clarify the nature of the new candidate substructure.

We describe a program to measure surface brightness fluctuation (SBF) distances to galaxies observed in the Next Generation Virgo Cluster Survey (NGVS), a photometric imaging survey covering $104~deg^2$ of the Virgo cluster in the ${u}^*,g,i,z$ bandpasses with the Canada-France Hawaii Telescope. We describe the selection of the sample galaxies, the procedures for measuring the apparent $i$-band SBF magnitude $\bar{i}$, and the calibration of the absolute $\bar{M}_i$ as a function of observed stellar population properties. The multi-band NGVS data set provides multiple options for calibrating the SBF distances, and we explore various calibrations involving individual color indices as well as combinations of two different colors. Within the color range of the present sample, the two-color calibrations do not significantly improve the scatter with respect to wide-baseline, single-color calibrations involving $u^{*}$. We adopt the ${u}^*{-}z$ calibration as reference for the present galaxy sample, with an observed scatter of 0.11 mag. For a few cases that lack good ${u}^*$ photometry, we use an alternative relation based on a combination of $g{-}i$ and $g{-}z$ colors, with only a slightly larger observed scatter of 0.12 mag. The agreement of our measurements with the best existing distance estimates provides confidence that our measurements are accurate. We present a preliminary catalog of distances for 89 galaxies brighter than $B_T\approx13.0$ mag within the survey footprint, including members of the background M and W Clouds at roughly twice the distance of the main body of the Virgo cluster. The extension of the present work to fainter and bluer galaxies is in progress.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band H$\alpha$+[NII] imaging survey carried out with MegaCam at the CFHT. During pilot observations we have observed NGC 4330, an intermediate mass, edge-on star forming spiral currently falling into the core of the Virgo cluster. New deep observations revealed a low surface brightness 10 kpc tail exhibiting a peculiar filamentary structure. The filaments are remarkably parallel one another and clearly indicate the direction of motion of the galaxy in the Virgo potential. Motivated by the detection of these features, indicating ongoing gas stripping, we collected literature photometry in 15 bands from the far-UV to the far-IR and deep optical long slit spectroscopy using the FORS2 instrument at the ESO Very Large Telescope. Using a newly developed Monte Carlo code that jointly fits spectroscopy and photometry, we reconstructed the star formation histories in apertures along the major axis of the galaxy. Our results have been validated against the output of CIGALE, a fitting code which has been previously used for similar studies. We found a clear outside-in gradient with radius of the time when the quenching event started: the outermost radii have been stripped 500 Myr ago, while the stripping has reached the inner 5 kpc from the center in the last 100 Myr. Regions at even smaller radii are currently still forming stars fueled by the presence of HI and H2 gas. When compared to statistical studies of the quenching timescales in the local Universe we find that ram pressure stripping of the cold gas is an effective mechanism to reduce the transformation times for galaxies falling into massive clusters. Future systematic studies of all the active galaxies observed by VESTIGE in the Virgo cluster will extend these results to a robust statistical framework.

The Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) is a blind narrow-band Halpha+[NII] imaging survey carried out with MegaCam at the Canada-France-Hawaii Telescope. The survey covers the whole Virgo cluster region from its core to one virial radius (104 deg^2). The sensitivity of the survey is of f(Halpha) ~ 4 x 10^-17 erg sec-1 cm^-2 (5 sigma detection limit) for point sources and Sigma (Halpha) ~ 2 x 10^-18 erg sec^-1 cm^-2 arcsec^-2 (1 sigma detection limit at 3 arcsec resolution) for extended sources, making VESTIGE the deepest and largest blind narrow-band survey of a nearby cluster. This paper presents the survey in all its technical aspects, including the survey design, the observing strategy, the achieved sensitivity in both the narrow-band Halpha+[NII] and in the broad-band r filter used for the stellar continuum subtraction, the data reduction, calibration, and products, as well as its status after the first observing semester. We briefly describe the Halpha properties of galaxies located in a 4x1 deg^2 strip in the core of the cluster north of M87, where several extended tails of ionised gas are detected. This paper also lists the main scientific motivations of VESTIGE, which include the study of the effects of the environment on galaxy evolution, the fate of the stripped gas in cluster objects, the star formation process in nearby galaxies of different type and stellar mass, the determination of the Halpha luminosity function and of the Halpha scaling relations down to ~ 10^6 Mo stellar mass objects, and the reconstruction of the dynamical structure of the Virgo cluster. This unique set of data will also be used to study the HII luminosity function in hundreds of galaxies, the diffuse Halpha+[NII] emission of the Milky Way at high Galactic latitude, and the properties of emission line galaxies at high redshift.

We present a photometric study of the dwarf galaxy population in the core region ($< r_{\rm vir}/4$) of the Fornax galaxy cluster based on deep $u'g'i'$ photometry from the Next Generation Fornax Cluster Survey. All imaging data were obtained with the Dark Energy Camera mounted on the 4-meter Blanco telescope at the Cerro-Tololo Interamerican Observatory. We identify 258 dwarf galaxy candidates with luminosities $-17 < M_{g'} < -8$ mag, corresponding to typical stellar masses of $9.5\gtrsim \log{\cal M}_{\star}/M_\odot \gtrsim 5.5$, reaching $\sim\!3$ mag deeper in point-source luminosity and $\sim\!4$ mag deeper in surface-brightness sensitivity compared to the classic Fornax Cluster Catalog. Morphological analysis shows that surface-brightness profiles are well represented by single-component S\'ersic models with average S\'ersic indices of $\langle n\rangle_{u',g',i'}=(0.78-0.83) \pm 0.02$, and average effective radii of $\langle r_e\rangle_{u',g',i'}\!=(0.67-0.70) \pm 0.02$ kpc. Color-magnitude relations indicate a flattening of the galaxy red sequence at faint galaxy luminosities, similar to the one recently discovered in the Virgo cluster. A comparison with population synthesis models and the galaxy mass-metallicity relation reveals that the average faint dwarf galaxy is likely older than ~5 Gyr. We study galaxy scaling relations between stellar mass, effective radius, and stellar mass surface density over a stellar mass range covering six orders of magnitude. We find that over the sampled stellar mass range several distinct mechanisms of galaxy mass assembly can be identified: i) dwarf galaxies assemble mass inside the half-mass radius up to $\log{\cal M}_{\star}$ ~8.0, ii) isometric mass assembly in the range $8.0 < \log{\cal M}_{\star}/M_\odot < 10.5$, and iii) massive galaxies assemble stellar mass predominantly in their halos at $\log{\cal M}_{\star}$ ~10.5 and above.

We present an analysis of 39 nuclei and their early-type hosts in the Virgo Cluster using ten broadband filters: F300W, F475W, F850LP, F160W, $u^*griz$, and $K_s$. We describe the Virgo Redux program, which provides high-resolution UV and NIR imaging. Combining this data with optical and NIR imaging from the ACS Virgo Cluster Survey and the Next Generation Virgo Cluster Survey, we estimate masses, metallicities and ages using simple stellar population (SSP) models. For 19 nuclei, we compare to SSP parameters derived from Keck and Gemini spectra and find reasonable agreement between the photometric and spectroscopic metallicity: the RMS scatter is 0.3 dex. We reproduce the nucleus-galaxy mass fraction of $0.33^{+0.09}_{-0.07}$ percent for galaxy stellar masses $10^{8.4}-10^{10.3} M_\odot$ with a typical precision of $\sim$35% for the nuclei masses. Based on available model predictions, there is no single preferred formation scenario for nuclei, suggesting that nuclei are formed stochastically through a mix of processes. Nuclei metallicities are statistically identical to those of their hosts, appearing $0.07 \pm 0.3$ dex more metal-rich on average -- although, omitting galaxies with unusual origins, nuclei are $0.20\pm0.28$ dex more metal-rich. Nuclei appear to be $0.56 \pm 0.12$ dex more metal rich than ultra-compact dwarf galaxies (UCDs) at fixed mass. We find no clear age difference between nuclei and their galaxies, with nuclei displaying a broad range of ages. Interestingly, we find that the most massive nuclei may be flatter and more closely aligned with the semi-major axes of their hosts, suggesting that they formed through predominantly dissipative processes.

Large samples of globular clusters (GC) with precise multi-wavelength photometry are becoming increasingly available and can be used to constrain the formation history of galaxies. We present the results of an analysis of Milky Way (MW) and Virgo core GCs based on five optical-near-infrared colors and ten synthetic stellar population models. For the MW GCs, the models tend to agree on photometric ages and metallicities, with values similar to those obtained with previous studies. When used with Virgo core GCs, for which photometry is provided by the Next Generation Virgo cluster Survey (NGVS), the same models generically return younger ages. This is a consequence of the systematic differences observed between the locus occupied by Virgo core GCs and models in panchromatic color space. Only extreme fine-tuning of the adjustable parameters available to us can make the majority of the best-fit ages old. Although we cannot exclude that the formation history of the Virgo core may lead to more conspicuous populations of relatively young GCs than in other environments, we emphasize that the intrinsic properties of the Virgo GCs are likely to differ systematically from those assumed in the models. Thus, the large wavelength coverage and photometric quality of modern GC samples, such as used here, is not by itself sufficient to better constrain the GC formation histories. Models matching the environment-dependent characteristics of GCs in multi-dimensional color space are needed to improve the situation.

We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening.

Low surface brightness galaxies (LSBGs) represent a significant percentage of local galaxies but their formation and evolution remain elusive. They may hold crucial information for our understanding of many key issues (i.e., census of baryonic and dark matter, star formation in the low density regime, mass function). The most massive examples - the so called giant LSBGs - can be as massive as the Milky Way, but with this mass being distributed in a much larger disk. Malin 1 is an iconic giant LSBG, perhaps the largest disk galaxy known. We attempt to bring new insights on its structure and evolution on the basis of new images covering a wide range in wavelength. We have computed surface brightness profiles (and average surface brightnesses in 16 regions of interest), in six photometric bands (FUV, NUV, u, g, i, z). We compared these data to various models, testing a variety of assumptions concerning the formation and evolution of Malin 1. We find that the surface brightness and color profiles can be reproduced by a long and quiet star-formation history due to the low surface density; no significant event, such as a collision, is necessary. Such quiet star formation across the giant disk is obtained in a disk model calibrated for the Milky Way, but with an angular momentum approximately 20 times larger. Signs of small variations of the star-formation history are indicated by the diversity of ages found when different regions within the galaxy are intercompared.For the first time, panchromatic images of Malin 1 are used to constrain the stellar populations and the history of this iconic example among giant LSBGs. Based on our model, the extreme disk of Malin 1 is found to have a long history of relatively low star formation (about 2 Msun/yr). Our model allows us to make predictions on its stellar mass and metallicity.

We present an analysis of high-quality photometry for globular clusters (GCs) in the Virgo cluster core region, based on data from the Next Generation Virgo Cluster Survey (NGVS) pilot field, and in the Milky Way (MW) based on VLT/X-Shooter spectrophotometry. We find significant discrepancies in color-color diagrams between sub-samples from different environments, confirming that the environment has a strong influence on the integrated colors of GCs. GC color distributions along a single color are not sufficient to capture the differences we observe in color-color space. While the average photometric colors become bluer with increasing radial distance to the cD galaxy M87, we also find a relation between the environment and the slope and intercept of the color-color relations. A denser environment seems to produce a larger dynamic range in certain color indices. We argue that these results are not due solely to differential extinction, IMF variations, calibration uncertainties, or overall age/metallicity variations. We therefore suggest that the relation between the environment and GC colors is, at least in part, due to chemical abundance variations, which affect stellar spectra and stellar evolution tracks. Our results demonstrate that stellar population diagnostics derived from model predictions which are calibrated on one particular sample of GCs may not be appropriate for all extragalactic GCs. These results advocate a more complex model of the assembly history of GC systems in massive galaxies that goes beyond the simple bimodality found in previous decades.

The central region of the Virgo cluster of galaxies contains thousands of globular clusters (GCs), an order of magnitude more than the numbers found in the Local Group. Relics of early star formation epochs in the universe, these GCs also provide ideal targets to test our understanding of the Spectral Energy Distributions (SEDs) of old stellar populations. Based on photometric data from the Next Generation Virgo cluster Survey (NGVS) and its near-infrared counterpart NGVS-IR, we select a robust sample of 1846 GCs with excellent photometry and spanning the full range of colors present in the Virgo core. The selection exploits the well defined locus of GCs in the uiK diagram and the fact that the globular clusters are marginally resolved in the images. We show that the GCs define a narrow sequence in 5-dimensional color space, with limited but real dispersion around the mean sequence. The comparison of these SEDs with the predictions of eleven widely used population synthesis models highlights differences between models, and also shows that no single model adequately matches the data in all colors. We discuss possible causes for some of these discrepancies. Forthcoming papers of this series will examine how best to estimate photometric metallicities in this context, and compare the Virgo globular cluster colors with those in other environments.

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $\alpha=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $\alpha = -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $\alpha=-1.21\pm0.05$.

We present an investigation of the globular cluster (GC) systems of NGC 3607 and NGC 3608 as part of the ongoing SLUGGS survey. We use wide-field imaging data from the Subaru telescope in the g, r, and i filters to analyse the radial density, colour and azimuthal distributions of both GC systems. With the complementary kinematic data obtained from the Keck II telescope, we measure the radial velocities of a total of 81 GCs. Our results show that the GC systems of NGC 3607 and NGC 3608 have a detectable spatial extent of ~ 15, and 13 galaxy effective radii, respectively. Both GC systems show a clear bimodal colour distribution. We detect a significant radial colour gradient for the GC subpopulations in both galaxies. NGC 3607 exhibits an overabundance of red GCs on the galaxy minor axis and NGC 3608 shows a misalignment in the GC subpopulation position angles with respect to the galaxy stellar component. With the aid of literature data, we discuss several relationships between the properties of GC systems and their host galaxies. A one-to-one relation between the ellipticities of red GCs and the galaxy stellar light emphasises the evolutionary similarities between them. In our sample of four slowly rotating galaxies with kinematically decoupled cores, we observe a higher ellipticity for the blue GC subpopulation than their red counterparts. Also, we notice the flattening of negative colour gradients for the blue GC subpopulations with increasing galaxy stellar mass. Finally, we discuss the formation scenarios associated with the blue GC subpopulation.

(Abridged) We investigate the intrinsic shapes of low-luminosity galaxies in the central 300 kpc of the Virgo cluster using deep imaging obtained as part of the NGVS. We build a sample of nearly 300 red-sequence cluster members in the yet unexplored $-14 < M_{g} < -8$ magnitude range. The observed distribution of apparent axis ratios is then fit by families of triaxial models with normally-distributed intrinsic ellipticities and triaxialities. We develop a Bayesian framework to explore the posterior distribution of the model parameters, which allows us to work directly on discrete data, and to account for individual, surface brightness-dependent axis ratio uncertainties. For this population we infer a mean intrinsic ellipticity E=0.43, and a mean triaxiality T=0.16. This implies that faint Virgo galaxies are best described as a family of thick, nearly oblate spheroids with mean intrinsic axis ratios 1:0.94:0.57. We additionally attempt a study of the intrinsic shapes of Local Group satellites of similar luminosities. For the LG population we infer a slightly larger mean intrinsic ellipticity E=0.51, and the paucity of objects with round apparent shapes translates into more triaxial mean shapes, 1:0.76:0.49. We finally compare the intrinsic shapes of NGVS low-mass galaxies with samples of more massive quiescent systems, and with field, star-forming galaxies of similar luminosities. We find that the intrinsic flattening in this low-luminosity regime is almost independent of the environment in which the galaxy resides--but there is a hint that objects may be slightly rounder in denser environments. The comparable flattening distributions of low-luminosity galaxies that have experienced very different degrees of environmental effects suggests that internal processes are the main drivers of galaxy structure at low masses--with external mechanisms playing a secondary role.

Due to longer dynamical timescales, the outskirts of early-type galaxies retain the footprint of their formation and assembly. Under the popular two-phase galaxy formation scenario, an initial in-situ phase of star formation is followed by minor merging and accretion of ex-situ stars leading to the expectation of observable transitions in the kinematics and stellar populations on large scales. However, observing the faint galactic outskirts is challenging, often leaving the transition unexplored. The large scale, spatially-resolved stellar kinematic data from the SAGES Legacy Unifying Galaxies and GlobularS (SLUGGS) survey are ideal for detecting kinematic transitions. We present kinematic maps out to 2.6 effective radii on average, kinemetry profiles, measurement of kinematic twists and misalignments, and the average outer intrinsic shape of 25 SLUGGS galaxies. We find good overall agreement in the kinematic maps and kinemetry radial profiles with literature. We are able to confirm significant radial modulations in rotational versus pressure support of galaxies with radius so that the central and outer rotational properties may be quite different. We also test the suggestion that galaxies may be more triaxial in their outskirts and find that while fast rotating galaxies were already shown to be axisymmetric in their inner regions, we are unable to rule out triaxiality in their outskirts. We compare our derived outer kinematic information to model predictions from a two-phase galaxy formation scenario. We find that the theoretical range of local outer angular momentum agrees well with our observations, but that radial modulations are much smaller than predicted.

We report on the properties of the most massive ultra-compact dwarf galaxy (UCD) in the nearby Virgo Cluster of galaxies using imaging from the Next Generation Virgo Cluster Survey (NGVS) and spectroscopy from Keck/DEIMOS. This object (M59-UCD3) appears to be associated with the massive Virgo galaxy M59 (NGC 4621), has an integrated velocity dispersion of 78 km/s, a dynamical mass of $3.7\times10^8 M_\odot$, and an effective radius ($R_e$) of 25 pc. With an effective surface mass density of $9.4\times10^{10} M_\odot/kpc^2$, it is the densest galaxy in the local Universe discovered to date, surpassing the density of the luminous Virgo UCD, M60-UCD1. M59-UCD3 has a total luminosity of $M_{g'}=-14.2$ mag, and a spectral energy distribution consistent with an old (14 Gyr) stellar population with [Fe/H]=0.0 and [$\alpha$/Fe]=+0.2. We also examine deep imaging around M59 and find a broad low surface brightness stream pointing towards M59-UCD3, which may represent a tidal remnant of the UCD progenitor. This UCD, along with similar objects like M60-UCD1 and M59cO, likely represents an extreme population of tidally stripped galaxies more akin to larger and more massive compact early-type galaxies than to nuclear star clusters in present-day dwarf galaxies.

Mass-to-light versus colour relations (MLCRs), derived from stellar population synthesis models, are widely used to estimate galaxy stellar masses (M$_*$) yet a detailed investigation of their inherent biases and limitations is still lacking. We quantify several potential sources of uncertainty, using optical and near-infrared (NIR) photometry for a representative sample of nearby galaxies from the Virgo cluster. Our method for combining multi-band photometry with MLCRs yields robust stellar masses, while errors in M$_*$ decrease as more bands are simultaneously considered. The prior assumptions in one's stellar population modelling dominate the error budget, creating a colour-dependent bias of up to 0.6 dex if NIR fluxes are used (0.3 dex otherwise). This matches the systematic errors associated with the method of spectral energy distribution (SED) fitting, indicating that MLCRs do not suffer from much additional bias. Moreover, MLCRs and SED fitting yield similar degrees of random error ($\sim$0.1-0.14 dex) when applied to mock galaxies and, on average, equivalent masses for real galaxies with M$_* \sim$ 10$^{8-11}$ M$_{\odot}$. The use of integrated photometry introduces additional uncertainty in M$_*$ measurements, at the level of 0.05-0.07 dex. We argue that using MLCRs, instead of time-consuming SED fits, is justified in cases with complex model parameter spaces (involving, for instance, multi-parameter star formation histories) and/or for large datasets. Spatially-resolved methods for measuring M$_*$ should be applied for small sample sizes and/or when accuracies less than 0.1 dex are required. An Appendix provides our MLCR transformations for ten colour permutations of the $grizH$ filter set.

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$.

We perform in-depth dynamical modelling of the luminous and dark matter (DM) content of the elliptical galaxy NGC 1407. Our strategy consists of solving the spherical Jeans equations for three independent dynamical tracers: stars, blue GCs and red GCs in a self-consistent manner. We adopt a maximum-likelihood Markov-Chain Monte Carlo fitting technique in the attempt to constrain the inner slope of the DM density profile (the cusp/core problem), and the stellar initial mass function (IMF) of the galaxy. We find the inner logarithmic slope of the DM density profiles to be $\gamma = 0.6\pm0.4$, which is consistent with either a DM cusp ($\gamma = 1$) or with a DM core $(\gamma = 0)$. Our findings are consistent with a Salpeter IMF, and marginally consistent with a Kroupa IMF. We infer tangential orbits for the blue GCs, and radial anisotropy for red GCs and stars. The modelling results are consistent with the virial mass--concentration relation predicted by $\Lambda$CDM simulations. The virial mass of NGC 1407 is $\log$ $M_{\rm vir} = 13.3 \pm 0.2 M_{\odot}$, whereas the stellar mass is $\log M_* = 11.8 \pm 0.1 M_{\odot}$. The overall uncertainties on the mass of NGC 1407 are only 5 per cent at the projected stellar effective radius. We attribute the disagreement between our results and previous X-ray results to the gas not being in hydrostatic equilibrium in the central regions of the galaxy. The halo of NGC 1407 is found be DM dominated, with a dynamical mass-to-light ratio of $M/L=260_{-100} ^{+174} M_{\odot}/L_{\odot, B}$. However, this value can be larger up to a factor of 3 depending on the assumed prior on the DM scale radius.

We introduce and provide the scientific motivation for a wide-field photometric and spectroscopic chemodynamical survey of nearby early-type galaxies (ETGs) and their globular cluster (GC) systems. The SLUGGS (SAGES Legacy Unifying Globulars and GalaxieS) survey is being carried out primarily with Subaru/Suprime-Cam and Keck/DEIMOS. The former provides deep gri imaging over a 900 arcmin^2 field-of-view to characterize GC and host galaxy colors and spatial distributions, and to identify spectroscopic targets. The NIR CaII triplet provides GC line-of-sight velocities and metallicities out to typically ~8 R_e, and to ~15 R_e in some cases. New techniques to extract integrated stellar kinematics and metallicities to large radii (~2-3 R_e) are used in concert with GC data to create two-dimensional velocity and metallicity maps for comparison with simulations of galaxy formation. The advantages of SLUGGS compared with other, complementary, 2D-chemodynamical surveys are its superior velocity resolution, radial extent, and multiple halo tracers. We describe the sample of 25 nearby ETGs, the selection criteria for galaxies and GCs, the observing strategies, the data reduction techniques, and modeling methods. The survey observations are nearly complete and more than 30 papers have so far been published using SLUGGS data. Here we summarize some initial results, including signatures of two-phase galaxy assembly, evidence for GC metallicity bimodality, and a novel framework for the formation of extended star clusters and ultracompact dwarfs. An integrated overview of current chemodynamical constraints on GC systems points to separate, in-situ formation modes at high redshifts for metal-poor and metal-rich GCs.

We present far-UV spectroscopy from the Cosmic Origins Spectrograph on the Hubble Space Telescope of a cool, star-forming filament in the core of Abell 1795. These data, which span 1025A - 1700A, allow for the simultaneous modeling of the young stellar populations and the intermediate-temperature (10^5.5 K) gas in this filament, which is far removed (~30 kpc) from the direct influence of the central AGN. Using a combination of UV absorption line indices and stellar population synthesis modeling, we find evidence for ongoing star formation, with the youngest stars having ages of 7.5 +/- 2.0 Myr and metallicities of 0.4 +/- 0.2 Zsun. The latter is consistent with the local metallicity of the intracluster medium. We detect the O VI (1038) line, measuring a flux of 4.0 +/- 0.9 x 10^-17 erg s^-1 cm^-2. The O VI (1032) line is redshifted such that it is coincident with a strong Galactic H2 absorption feature, and is not detected. The measured O VI (1038) flux corresponds to a cooling rate of 0.85 +/- 0.2 (stat) +/- 0.15 (sys) Msun/yr at ~10^5.5 K, assuming that the cooling proceeds isochorically, which is consistent with the classical X-ray luminosity-derived cooling rate in the same region. We measure a star formation rate of 0.11 +/- 0.02 Msun/yr from the UV continuum, suggesting that star formation is proceeding at 13 +/- 3% efficiency in this filament. We propose that this inefficient star formation represents a significant contribution to the larger-scale cooling flow problem.

Optimised population synthesis provides an empirical method to extract the relative mix of stellar evolutionary stages and the distribution of atmospheric parameters within unresolved stellar systems, yet a robust validation of this method is still lacking. We here provide a calibration of population synthesis via non-linear bound-constrained optimisation of stellar populations based upon optical spectra of mock stellar systems and observed Galactic Globular Clusters (GGCs). The MILES stellar library is used as a basis for mock spectra as well as templates for the synthesis of deep GGC spectra from Schiavon et al. (2005). Optimised population synthesis applied to mock spectra recovers mean light-weighted stellar atmospheric parameters to within a mean uncertainty of 240 K, 0.04 dex, and 0.03 dex for T_eff, log(g), and [Fe/H], respectively. Decompositions of both mock and GGC spectra confirm the method's ability to recover the expected mean light-weighted metallicity in dust-free conditions (E[B-V] < 0.15) with uncertainties comparable to evolutionary population synthesis methods. Dustier conditions require either appropriate dust-modelling when fitting to the full spectrum, or fitting only to select spectral features. We derive light-weighted fractions of stellar evolutionary stages from our population synthesis fits to GGCs, yielding on average a combined 25+/-6 per cent from main sequence and turnoff dwarfs, 64+/-7 per cent from subgiant, red giant and asymptotic giant branch stars, and 15+/-7 per cent from horizontal branch stars and blue stragglers. Excellent agreement is found between these fractions and those estimated from deep HST/ACS CMDs. Overall, optimised population synthesis remains a powerful tool for understanding the stellar populations within the integrated light of galaxies and globular clusters.

We present an extenstive literature compilation of age, metallicity, and chemical abundance pattern information for the 41 Galactic globular clusters (GGCs) studied by Schiavon et al. (2005). Our compilation constitutes a notable improvement over previous similar work, particularly in terms of chemical abundances. Its primary purpose is to enable detailed evaluations of and refinements to stellar population synthesis models designed to recover the above information for unresolved stellar systems based on their integrated spectra. However, since the Schiavon sample spans a wide range of the known GGC parameter space, our compilation may also benefit investigations related to a variety of astrophysical endeavours, such as the early formation of the Milky Way, the chemical evolution of GGCs, and stellar evolution and nucleosynthesis. For instance, we confirm with our compiled data that the GGC system has a bimodal metallicity distribution and is uniformly enhanced in the alpha-elements. When paired with the ages of our clusters, we find evidence that supports a scenario whereby the Milky Way obtained its globular clusters through two channels, in situ formation and accretion of satellite galaxies. The distributions of C, N, O, and Na abundances and the dispersions thereof per cluster corroborate the known fact that all GGCs studied so far with respect to multiple stellar populations have been found to harbour them. Finally, using data on individual stars, we also confirm that the atmospheres of stars become progressively polluted by CN(O)-processed material after they leave the main sequence and uncover evidence which suggests the alpha-elements Mg and Ca may originate from more than one nucleosynthetic production site. [abridged]

We present stellar age profiles for 64 Virgo cluster disk galaxies whose analysis poses a challenge for current galaxy formation models. Our results can be summarized as follows: first, and contrary to observations of field galaxies, these cluster galaxies are distributed almost equally amongst the three main types of disk galaxy luminosity profiles (I/II/III), indicating that the formation and/or survival of Type II breaks is suppressed within the cluster environment. Second, we find examples of statistically-significant inversions ("U-shapes") in the age profiles of all three disk galaxy types, reminescent of predictions from high-resolution simulations of classically-truncated Type II disks in the field. These features characterize the age profiles for only about a third (<36%) of each disk galaxy type in our sample. An even smaller fraction of cluster disks (~11% of the total sample), exhibit age profiles which decrease outwards (i.e., negative age gradients). Instead, flat and/or positive age gradients prevail (>50%) within our Type I, II and III sub-samples. These observations thus suggest that while stellar migrations and inside-out growth can play a significant role in the evolution of all disk galaxy types, other factors contributing to the evolution of galaxies can overwhelm the predicted signatures of these processes. We interpret our observations through a scenario whereby Virgo cluster disk galaxies formed initially like their bretheren in the field but which, upon falling into the cluster, were transformed into their present state through external processes linked to the environment. Current disk galaxy formation models fail to reproduce these results, thus calling for adequate hydrodynamical simulations of dense galaxy environments, for which the current paper provides many constraints. [Abridged]

We present griz and H-band surface brightness profiles and bulge-disk decompositions for a morphologically-broad sample of 286 Virgo cluster galaxies. The H-band data come from a variety of sources including our survey of 171 galaxies at the UH 2.2-m, CFHT and UKIRT telescopes, and another 115 galaxies from the Two-Micron All-Sky Survey and GOLDMine archives. The optical data for all 286 Virgo galaxies were extracted from Sloan Digital Sky Survey (SDSS) images. The H-band and SDSS griz data were analyzed in a homogeneous manner using our own software, yielding a consistent set of deep, multi-band surface brightness profiles for each galaxy. Average surface brightness profiles per morphological bin were created in order to characterize the variety of galaxy light profiles across the Hubble sequence. The 1D bulge-disk decomposition parameters, as well as non-parametric galaxy measures, such as effective radius, effective surface brightness and light concentration, are presented for all 286 galaxies in each of the five optical/near-infrared wavebands. The profile decompositions account for bulge and disk components, spiral arms, nucleus and atmospheric blurring. The Virgo spiral galaxy bulges typically have a Sersic index n ~ 1, while elliptical galaxies prefer n ~ 2. No galaxy spheroid requires n > 3. The light profiles for 70% of the Virgo elliptical galaxies reveal the presence of both a spheroid and disk component. A more in-depth discussion of the structural parameter trends can be found in McDonald et al (2009b). The data provided here should serve as a base for studies of galaxy structure and stellar populations in the cluster environment. The galaxy light profiles and bulge-disk decomposition results are available at the Centre de Donnees astronomiques de Strasbourg (CDS; http://cds.u-strasbg.fr) and the author's own website (http://www.astro.queensu.ca/virgo).

We use a combination of deep optical (gri) and near-infrared (H) photometry to study the radially-resolved colours of a broad sample of 300 Virgo cluster galaxies. For most galaxy types, we find that the median g-H colour gradient is either flat (gas-poor giants and gas-rich dwarfs) or negative (i.e., colours become bluer with increasing radius; gas-poor dwarfs, spirals, and gas-poor peculiars). Later-type galaxies typically exhibit more negative gradients than early-types. Given the lack of a correlation between the central colours and axis ratios of Virgo spiral galaxies, we argue that dust likely plays a small role, if at all, in setting those colour gradients. We search for possible correlations between galaxy colour and photometric structure or environment and find that the Virgo galaxy colours become redder with increasing concentration, luminosity and surface brightness, while no dependence with cluster-centric radius or local galaxy density is detected (over a range of ~2 Mpc and ~3-16 Mpc^-2, respectively). However, the colours of gas-rich Virgo galaxies do correlate with neutral gas deficiency, such that these galaxies become redder with higher deficiencies. Comparisons with stellar population models suggest that these colour gradients arise principally from variations in stellar metallicity within these galaxies, while age variations only make a significant contribution to the colour gradients of Virgo irregulars. A detailed stellar population analysis based on this material is presented in Roediger et al (2011b; arXiv:1011.3511).

We use a combination of deep optical and near-infrared light profiles for a morphologically diverse sample of Virgo cluster galaxies to study the radially-resolved stellar populations of cluster galaxies over a wide range of galaxy structure. We find that, in the median, the age gradients of Virgo galaxies are either flat (lenticulars and Sa-Sb spirals) or positive (ellipticals, Sbc+Sc spirals, gas-rich dwarfs, and irregulars), while all galaxy types have a negative median metallicity gradient. Comparison of the galaxy stellar population diagnostics (age, metallicity, and gradients thereof) against structural and environmental parameters also reveals that the ages of gas-rich systems depend mainly on their atomic gas deficiencies. Conversely, the metallicities of Virgo gas-poor galaxies depend on their concentrations, luminosities, and surface brightnesses. The stellar population gradients of all Virgo galaxies exhibit no dependence on either their structure or environment. We interpret these stellar population data for Virgo galaxies in the context of popular formation and evolution scenarios, and suggest that gas-poor giants grew hierarchically (through dissipative starbursts), gas-poor dwarfs have descended from at least two different production channels (e.g., environmental transformation and merging), while spirals formed inside-out, but with star formation in the outskirts of a significant fraction of the population having been quenched due to ram pressure stripping. (Abridged)