(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We investigate the dependence of galaxy structure on a variety of galactic and environmental parameters for ~500,000 galaxies at z<0.2, taken from the Sloan Digital Sky Survey data release 7 (SDSS-DR7). We utilise bulge-to-total stellar mass ratio, (B/T)_*, as the primary indicator of galactic structure, which circumvents issues of morphological dependence on waveband. We rank galaxy and environmental parameters in terms of how predictive they are of galaxy structure, using an artificial neural network approach. We find that distance from the star forming main sequence (Delta_SFR), followed by stellar mass (M_*), are the most closely connected parameters to (B/T)_*, and are significantly more predictive of galaxy structure than global star formation rate (SFR), or any environmental metric considered (for both central and satellite galaxies). Additionally, we make a detailed comparison to the Illustris hydrodynamical simulation and the LGalaxies semi-analytic model. In both simulations, we find a significant lack of bulge-dominated galaxies at a fixed stellar mass, compared to the SDSS. This result highlights a potentially serious problem in contemporary models of galaxy evolution.

Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.

We present Gemini/GNIRS spectroscopy of the Seyfert 2 galaxy NGC 4388, with simultaneous coverage from 0.85 - 2.5 $\mu$m. Several spatially-extended emission lines are detected for the first time, both in the obscured and unobscured portion of the optical narrow line region (NLR), allowing us to assess the combined effects of the central continuum source, outflowing gas and shocks generated by the radio jet on the central 280 pc gas. The HI and [FeII] lines allow us to map the extinction affecting the NLR. We found that the nuclear region is heavily obscured, with E(B-V) ~1.9 mag. To the NE of the nucleus and up to ~150 pc, the extinction remains large, ~1 mag or larger, consistent with the system of dust lanes seen in optical imaging. We derived position-velocity diagrams for the most prominent lines as well as for the stellar component. Only the molecular gas and the stellar component display a well-organized pattern consistent with disk rotation. Other emission lines are kinematically perturbed or show little evidence of rotation. Extended high-ionization emission of sulfur, silicon and calcium is observed to distances of at least 200 pc both NE and SW of the nucleus. We compared flux ratios between these lines with photoionization models and conclude that radiation from the central source alone cannot explain the observed high-ionization spectrum. Shocks between the radio-jet and the ambient gas are very likely an additional source of excitation. We conclude that NGC 4388 is a prime laboratory to study the interplay between all these mechanisms.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

We present the stellar mass functions (SMF) and mass densities of galaxies, and their spheroid and disk components in the local (z~0.1) universe over the range 8.9 <= log(M/M_solar) <= 12 from spheroid+disk decompositions and corresponding stellar masses of a sample of over 600,000 galaxies in the SDSS-DR7 spectroscopic sample. The galaxy SMF is well represented by a single Schechter function (M* = 11.116+/-0.011, alpha = -1.145+/-0.008), though with a hint of a steeper faint end slope. The corresponding stellar mass densities are (2.670+/-0.110), (1.687+/-0.063) and (0.910+/-0.029)x10^8 M_solar Mpc^-3 for galaxies, spheroids and disks respectively. We identify a crossover stellar mass of log(M/M_solar) = 10.3+/-0.030 at which the spheroid and disk SMFs are equal. Relative contributions of four distinct spheroid/disk dominated sub-populations to the overall galaxy SMF are also presented. The mean disk-to-spheroid stellar mass ratio shows a five fold disk dominance at the low mass end, decreasing monotonically with a corresponding increase in the spheroidal fraction till the two are equal at a galaxy stellar mass, log(M/M_solar)=10.479+/-0.013, the dominance of spheroids then grows with increasing stellar mass. The relative numbers of composite disk and spheroid dominated galaxies show peaks in their distributions, perhaps indicative of a preferred galaxy mass. Our characterization of the low redshift galaxy population provides stringent constraints for numerical simulations to reproduce.

We analyze the stellar absorption features in high signal-to-noise ratio near-infrared (NIR) spectra of the nuclear region of 12 nearby galaxies, mostly spirals. The features detected in some or all of the galaxies in this sample are the TiO (0.843 $\mu$m\ and 0.886 $\mu$m), VO (1.048 $\mu$m), CN (1.1 $\mu$m\ and 1.4 $\mu$m), H$\rm _2$O (1.4 $\mu$m\ and 1.9 $\mu$m) and CO (1.6 $\mu$m\ and 2.3 $\mu$m) bands. The C$\rm _2$ (1.17 $\mu$m\ and 1.76 $\mu$m) bands are generally weak or absent, although C$\rm _2$ (1.76 $\mu$m) may be weakly present in the mean galaxy spectrum. A deep feature near 0.93 $\mu$m, likely caused by CN, TiO and/or ZrO, is also detected in all objects. Fitting a combination of stellar spectra to the mean spectrum shows that the absorption features are produced by evolved stars: cool giants and supergiant stars in the early- or thermally-pulsing asymptotic giant branch (E-AGB or TP-AGB) phases. The high luminosity of TP-AGB stars, and the appearance of VO and ZrO features in the data, suggest that TP-AGB stars dominate these spectral features. However, a contribution from other evolved stars is also likely. Comparison with evolutionary population synthesis models shows that models based on empirical libraries that predict relatively strong NIR features provide a more accurate description of the data. However, none of the models tested accurately reproduces all of the features observed in the spectra. To do so, the models will need to not only improve the treatment of TP-AGB stars, but also include good quality spectra of red giant and E-AGB stars. The uninterrupted wavelength coverage, high S/N, and quantity of features we present here will provide a benchmark for the next generation of models aiming to explain and predict the NIR properties of galaxies.

We examine the stellar velocity dispersions (sigma) of a sample of 48 galaxies, 35 of which are spirals, from the Palomar nearby galaxy survey. It is known that for ultra-luminous infrared galaxies (ULIRGs) and merger remnants thesigma derived from the near-infrared CO band-heads is smaller than that measured from optical lines, while no discrepancy between these measurements is found for early-type galaxies. No such studies are available for spiral galaxies - the subject of this paper. We used cross-dispersed spectroscopic data obtained with the Gemini Near-Infrared Spectrograph (GNIRS), with spectral coverage from 0.85 to 2.5um, to obtain sigma measurements from the 2.29 $\mu$m CO band-heads (sigma_{CO}), and the 0.85 um calcium triplet (sigma_{CaT}). For the spiral galaxies in the sample, we found that sigma_{CO} is smaller than sigma_{CaT}, with a mean fractional difference of 14.3%. The best fit to the data is given by sigma_{opt} = (46.0+/-18.1) + (0.85+/-0.12)sigma_{CO}. This "sigma discrepancy" may be related to the presence of warm dust, as suggested by a slight correlation between the discrepancy and the infrared luminosity. This is consistent with studies that have found no sigma-discrepancy in dust-poor early-type galaxies, and a much larger discrepancy in dusty merger remnants and ULIRGs. That sigma_{CO}$ is lower than sigma_{opt} may also indicate the presence of a dynamically cold stellar population component. This would agree with the spatial correspondence between low sigma_{CO} and young/intermediate-age stellar populations that has been observed in spatially-resolved spectroscopy of a handful of galaxies.

We present spectra of the nuclear regions of 50 nearby (D = 1 - 92 Mpc, median = 20 Mpc) galaxies of morphological types E to Sm. The spectra, obtained with the Gemini Near-IR Spectrograph on the Gemini North telescope, cover a wavelength range of approximately 0.85-2.5 microns at R~1300--1800. There is evidence that most of the galaxies host an active galactic nucleus (AGN), but the range of AGN luminosities (log (L2-10 keV [erg/s]) = 37.0-43.2) in the sample means that the spectra display a wide variety of features. Some nuclei, especially the Seyferts, exhibit a rich emission-line spectrum. Other objects, in particular the type 2 Low Ionisation Nuclear Emission Region galaxies, show just a few, weak emission lines, allowing a detailed view of the underlying stellar population. These spectra display numerous absorption features sensitive to the stellar initial mass function, as well as molecular bands arising in cool stars, and many other atomic absorption lines. We compare the spectra of subsets of galaxies known to be characterised by intermediate-age and old stellar populations, and find clear differences in their absorption lines and continuum shapes. We also examine the effect of atmospheric water vapor on the signal-to-noise ratio achieved in regions between the conventional NIR atmospheric windows, of potential interest to those planning observations of redshifted emission lines or other features affected by telluric H2O. Further exploitation of this data set is in progress, and the reduced spectra and data reduction tools are made available to the community.

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

(Abridged) Distant galaxy clusters provide important tests of the growth of large scale structure in addition to highlighting the process of galaxy evolution in a consistently defined environment at large look back time. We present a sample of 22 distant (z>0.8) galaxy clusters and cluster candidates selected from the 9 deg2 footprint of the overlapping X-ray Multi Mirror (XMM) Large Scale Structure (LSS), CFHTLS Wide and Spitzer SWIRE surveys. Clusters are selected as extended X-ray sources with an accompanying overdensity of galaxies displaying optical to mid-infrared photometry consistent with z>0.8. Nine clusters have confirmed spectroscopic redshifts in the interval 0.8<z<1.2, four of which are presented here for the first time. A further 11 candidate clusters have between 8 and 10 band photometric redshifts in the interval 0.8<z<2.2, while the remaining two candidates do not have information in sufficient wavebands to generate a reliable photometric redshift. All of the candidate clusters reported in this paper are presented for the first time. Those confirmed and candidate clusters with available near infrared photometry display evidence for a red sequence galaxy population, determined either individually or via a stacking analysis, whose colour is consistent with the expectation of an old, coeval stellar population observed at the cluster redshift. We further note that the sample displays a large range of red fraction values indicating that the clusters may be at different stages of red sequence assembly. We compare the observed X-ray emission to the flux expected from a suite of model clusters and find that the sample displays an effective mass limit M200 ~ 1e14 Msolar with all clusters displaying masses consistent with M200 < 5e14 Msolar. This XMM distant cluster study represents a complete sample of X-ray selected z>0.8 clusters.

We present secure [NII[ detections in two mm-bright, strongly lensed objects at high redshift, APM08279+5255 (z=3.911) and MM18423+5938 (z=3.930), using the IRAM Plateau de Bure Interferometer. Due to its ionization energy [NII] is a good tracer of the ionized gas phase in the interstellar medium. The measured fluxes are S([NII])=(4.8+/-0.8) Jy km/s and (7.4+/-0.5) Jy km/s respectively, yielding line luminosities of L([NII]) =(1.8+/-0.3) x 10^9 \mu^{-1} Lsun for APM08279+5255 and L([NII]) =(2.8+/-0.2) x 10^9 \mu^{-1} Lsun for MM18423+5938. Our high-resolution map of the [NII] and 1 mm continuum emission in MM18423+5938 clearly resolves an Einstein ring in this source, and reveals a velocity gradient in the dynamics of the ionized gas. A comparison of these maps with high-resolution EVLA CO observations enables us to perform the first spatially-resolved study of the dust continuum-to-molecular gas surface brightness (Sigma_{FIR} Sigma_{CO}^N, which can be interpreted as the star formation law) in a high-redshift object. We find a steep relation (N=1.4+/-0.2), consistent with a starbursting environment. We measure a [NII]/FIR luminosity ratio in APM0828+5255 and MM18423+5938 of 9.0 x 10^{-6} and 5.8 x 10^{-6}, respectively. This is in agreement with the decrease of the [NII]/FIR ratio at high FIR luminosities observed in local galaxies.

We have used Gemini North together with the NIRI-ALTAIR adaptive optics imager in the H and K bands to explore the core of the Galactic globular cluster M71 (NGC 6838). We obtained proper motions for 217 stars and have resolved its internal proper motion dispersion. Using a 3.8 year baseline, the proper motion dispersion in the core is found to be 179 +/- 17 microarcsec/yr. We find no evidence of anisotropy in the motions and no radial variation in the proper motions with respect to distance from the cluster center. We also set an upper limit on any central black hole to be ~150 Msun at 90% confidence level.

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We study the 64 square degrees W1 field, the largest of the CFHTLS-Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet made. 2.66 million galaxy shapes are measured, using a KSB pipeline verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio {\nu}>3.5, consistent with predictions of a {\Lambda}CDM model. Of these peaks, 126 lie within 3.0' of a BCG identified from multicolor optical imaging in an independent, red sequence survey. We also identify 7 counterparts for massive clusters previously seen in X-ray emission within 6 square degrees XMM-LSS survey. With photometric redshift estimates for the source galaxies, we use a tomographic lensing method to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we confirm 85 groups/clusters with \chi^2 reduced < 3.0, at a mean redshift <z_c> = 0.36 and velocity dispersion <\sigma_c> = 658.8 km/s. Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to apply these same techniques to map clusters in much larger volumes and thus tightly constrain cosmological models.

We present high resolution imaging of the low order (J=1 and 2) CO line emission from the z = 3.93 submillimeter galaxy (SMG) MM18423+5938 using the Expanded Very Large Array, and optical and near-IR imaging using the Canada-France-Hawaii Telescope. This SMG with a spectroscopic redshift was thought to be gravitationally lensed given its enormous apparent brightness. We find that the CO emission is consistent with a complete Einstein ring with a major axis diameter of ~ 1.4", indicative of lensing. We have also identified the lens galaxy as a very red elliptical coincident with the geometric center of the ring and estimated its photometric redshift z~1.1. A first estimate of the lens magnification factor is m~12. The luminosity L'_CO(1-0) of the CO(1-0) emission is 2.71+/-0.38 x10^11/m K km s^-1 pc^2, and, adopting the commonly used conversion factor for ULIRGs, the molecular gas mass is M(H_2) = 2.2 x10^11/m M_sol, comparable to unlensed SMGs if corrected by m ~ 12. Our revised estimate of the far-IR luminosity of MM18423+5938 is 2 x10^13/m < L_FIR < 3 x10^14/m L_sol, comparable to that of ULIRGs. Further observations are required to quantify the star formation rate in MM18423+5938 and to constrain the mass model of the lens in more detail.

Using a combined analysis of strong lensing and galaxy dynamics, we characterize the mass distributions and M/L ratios of galaxy groups, which form an important transition regime in Lambda-CDM cosmology. By mapping the underlying mass distribution, we test whether groups are dark matter dominated as hypothesized by the standard cosmogony, or isothermal as observed in baryon rich field galaxies. We present our lensing + galaxy dynamics formalism built around the dark matter dominant NFW and Hernquist distributions, compared against the Isothermal Sphere observed in galaxy scale objects. We show that mass measurement in the core of the group (r ~ 0.2 r_{vir}), determined jointly from a lens model and from differential velocity dispersion estimates, may effectively distinguish between these density distributions. We apply our method to MOS observations of two groups, SL2SJ1430+5546 and SL2SJ1431+5533, drawn from our CFHTLS lens catalog. With the measured lensing and dynamical masses, combined with a maximum likelihood estimator built around our model, we estimate the concentration index characterizing each density distribution and the corresponding virial mass of each group. Our results indicate that both groups are dark matter dominant, and reject the Isothermal distribution at >>3 sigma level. For both groups, the estimated i-band M/L ratios of ~260 Msun/Lsun, are similar to other published values for groups. The Gaussian distributions of the velocities of their member galaxies support a high degree of virialization. The differences in their virial masses, 2.8 and 1.6 x 10^14 Msun, and velocity dispersions, 720 and 560 km/s respectively, may indicate however that each group is at a different stage of transition to a cluster. We aim to populate this important transition regime with additional results from ongoing observations of the remaining lensing groups in our catalog.

We have developed a new method, K2, optimized for the detection of galaxy clusters in multicolor images. Based on the Red Sequence approach, K2 detects clusters using simultaneous enhancements in both colors and position. The detection significance is robustly determined through extensive Monte-Carlo simulations and through comparison with available cluster catalogs based on two different optical methods, and also on X-ray data. K2 also provides quantitative estimates of the candidate clusters' richness and photometric redshifts. Initially K2 was applied to 161 sq deg of two color gri images of the CFHTLS-Wide data. Our simulations show that the false detection rate, at our selected threshold, is only ~1%, and that the cluster catalogs are ~80% complete up to a redshift of 0.6 for Fornax-like and richer clusters and to z ~0.3 for poorer clusters. Based on Terapix T05 release gri photometric catalogs, 35 clusters/sq deg are detected, with 1-2 Fornax-like or richer clusters every two square degrees. Catalogs containing data for 6144 galaxy clusters have been prepared, of which 239 are rich clusters. These clusters, especially the latter, are being searched for gravitational lenses -- one of our chief motivations for cluster detection in CFHTLS. The K2 method can be easily extended to use additional color information and thus improve overall cluster detection to higher redshifts. The complete set of K2 cluster catalogs, along with the supplementary catalogs for the member galaxies, are available on request from the authors.

The existence of strong lensing systems with Einstein radii (Re) covering the full mass spectrum, from ~1-2" (produced by galaxy scale dark matter haloes) to >10" (produced by galaxy cluster scale haloes) have long been predicted. Many lenses with Re around 1-2" and above 10" have been reported but very few in between. In this article, we present a sample of 13 strong lensing systems with Re in the range 3"- 8", i.e. systems produced by galaxy group scale dark matter haloes, spanning a redshift range from 0.3 to 0.8. This opens a new window of exploration in the mass spectrum, around 10^{13}- 10^{14} M_{sun}, which is a crucial range for understanding the transition between galaxies and galaxy clusters. Our analysis is based on multi-colour CFHTLS images complemented with HST imaging and ground based spectroscopy. Large scale properties are derived from both the light distribution of the elliptical galaxies group members and weak lensing of the faint background galaxy population. On small scales, the strong lensing analysis yields Einstein radii between 2.5" and 8". On larger scales, the strong lenses coincide with the peak of the light distribution, suggesting that mass is traced by light. Most of the luminosity maps have complicated shapes, indicating that these intermediate mass structures are dynamically young. Fitting the reduced shear with a Singular Isothermal Sphere, we find sigma ~ 500 km/s and an upper limit of ~900 km/s for the whole sample. The mass to light ratio for the sample is found to be M/L_i ~ 250 (solar units, corrected for evolution), with an upper limit of 500. This can be compared to mass to light ratios of small groups (with sigma ~ 300 km/s and galaxy clusters with sigma > 1000 km/s, thus bridging the gap between these mass scales.