Velocity and density field reconstructions of the volume of the universe within 0.05c derived from the Cosmicflows-3 catalog of galaxy distances has revealed the presence of a filamentary structure extending across ~ 0.11c. The structure, at a characteristic redshift of 12,000 km/s, has a density peak coincident with the celestial South Pole. This structure, the largest contiguous feature in the local volume and comparable to the Sloan Great Wall at half the distance, is given the name the South Pole Wall.

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

This study is a part of the Cosmicflows-4 project with the aim of measuring the distances of more than ~10,000 spiral galaxies in the local universe up to ~15,000 km/s. New HI linewidth information has come primarily from the Arecibo Legacy Fast ALFA Survey. Photometry of our sample galaxies has been carried out in optical (SDSS u, g, r, i, z) and infrared (WISE W1 and W2) bands. Inclinations have been determined using an online graphical interface accessible to a collaboration of citizen scientists. Galaxy distances are measured based on the correlation between the rotation rate of spirals and their absolute luminosity, known as the Tully-Fisher Relation (TFR). In this study, we present the calibration of the TFR using a subsample of ~600 spirals located in 20 galaxy clusters. Correlations among such observables as color, surface brightness, and relative HI content are explored in an attempt to reduce the scatter about the TFR with the goal of obtaining more accurate distances. A preliminary determination of the Hubble constant from the distances and velocities of the calibrator clusters is H0=76.0+-1.1(stat.)+-2.3(sys.) km/s/Mpc.

This article explores in depth a watershed concept to partition the universe, introduced in \cite{2019MNRAS.489L...1D} and applied to the {\it Cosmicflows-3} observational dataset. We present a series of tests conducted with cosmological dark matter simulations. In particular we are interested in quantifying the evolution with redshift of large scale structures when defined as segmented basins of attraction. This new dynamical definition in the field of measuring standard rulers demonstrates robustness since all basins show a density contrast $\delta$ above one (mean universe density) independently of the simulation spatial resolution or the redshift. Another major finding is that density profiles of the basins show universality in slope. Consequently, there is a unique definition of what is a gravitational watershed at large scale, that can be further used as a probe for cosmology studies.

Type Ia Supernovae have yet again the opportunity to revolutionize the field of cosmology as the new generation of surveys are acquiring thousands of nearby SNeIa opening a new era in cosmology: the direct measurement of the growth of structure parametrized by $fD$. This method is based on the SNeIa peculiar velocities derived from the residual to the Hubble law as direct tracers of the full gravitational potential caused by large scale structure. With this technique, we could probe not only the properties of dark energy, but also the laws of gravity. In this paper we present the analytical framework and forecasts. We show that ZTF and LSST will be able to reach 5\% precision on $fD$ by 2027. Our analysis is not significantly sensitive to photo-typing, but known selection functions and spectroscopic redshifts are mandatory. We finally introduce an idea of a dedicated spectrograph that would get all the required information in addition to boost the efficiency to each SNeIa so that we could reach the 5\% precision within the first two years of LSST operation and the few percent level by the end of the survey.

Secular extragalactic parallax caused by the solar system's velocity relative to the cosmic microwave background rest frame may be observable as a dipole proper motion field with amplitude $78~\mu$as yr$^{-1}$ Mpc. Nearby galaxies also exhibit proper motions caused by their transverse peculiar velocities that prevent detection of secular parallax for any single galaxy, although a statistical detection may be made instead. Such a detection could constrain the local Hubble parameter. We present methods to measure secular parallax using correlated extragalactic proper motions and find a first limit on the secular parallax amplitude using proper motions of 232 nearby galaxies from Gaia Data Release 2. The recovered dipole has insignificant upper limit of 3500 $\mu$as yr$^{-1}$ Mpc. This measurement will be improved by larger sample size and reduced proper motion uncertainties in future data releases. Using the local peculiar velocity field derived from Cosmicflows-3, we simulate galaxy proper motions and predict that a significant detection ($5-10\sigma$) of the secular parallax amplitude will be possible by Gaia's end of mission. The detection is contingent on proper motions of nearby ($<5$ Mpc), bright ($G<15$ mag) galaxies, and corresponds to an insignificant upper limit on the Hubble parameter. We further investigate the implications of our simulations for the study of transverse peculiar velocities, which we find to be consistent with large scale structure theory. The peculiar velocity field additionally results in low-multipole correlated proper motions on the order of $0.3~\mu$as yr$^{-1}$ that may be confounded with other cosmological proper motion measurements, such as limits on the gravitational wave background and the anisotropy of the Hubble expansion.

Tools are provided at the Extragalactic Distance Database website that provide relationships between the distances and velocities of galaxies based on smoothed versions of the velocity fields derived by the Cosmicflows program.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

The universe region obscured by the Milky Way is very large and only future blind large HI redshift, and targeted peculiar surveys on the outer borders will determine how much mass is hidden there. Meanwhile, we apply for the first time two independent techniques to the galaxy peculiar velocity catalog $CosmicFlows-3$ in order to explore for the kinematic signature of a specific large-scale structure hidden behind this zone : the Vela supercluster at cz $\sim 18,000$,km s$^{-1}$ . Using the gravitational velocity and density contrast fields, we find excellent agreement when comparing our results to the Vela object as traced in redshift space. The article provides the first kinematic evidence of a major mass concentration (knot of the Cosmic Web) located in the direction behind Vela constellation, pin-pointing that the Zone of Avoidance should be surveyed in detail in the future .

To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shear tensor corresponding to the direction of greatest collapse obtained by reconstructing the full velocity field in the linear regime. Here we extend that work by both considering an additional 5 planes beyond the 5 examined previously and by examining the alignment with respect to the quasi-linear field, a more sophisticated reconstruction, which is a better approximation on smaller (quasi-linear) scales. Our analysis recovers the previous result while not finding a significantly tight alignment with the additional 5 planes. However, the additional 5 planes normals also do not appear to be randomly oriented. We conclude that this could either be due to the normals of the new planes being poorly defined and described; the quasi-linear field at those locations being poorly constrained; or different formation mechanisms for the orientation of planes of dwarf galaxies.

The emerging light from a galaxy is under the influence of its own interstellar medium, as well as its spatial orientation. Considering a sample of 2,239 local spiral galaxies in optical (SDSS u, g, r, i, and z) and infrared bands (WISE W1, W1), we study the dependency of the global intrinsic attenuation in spiral galaxies on their morphologies, sizes, and spatial inclinations. Reddening is minimal at the extremes of low mass and gas depletion and maximal in galaxies that are relatively massive and metal-rich and still retain substantial gas reserves. A principal component constructed from observables that monitor galaxy mass, relative HI content to old stars, and infrared surface brightness is strongly correlated with the amplitude of obscuration. We determine both a parametric model for dust obscuration and a non-parametric model based on the Gaussian process formalism. An average dust attenuation curve is derived for wavelengths between 0.36 and 4.5 microns.

This letter presents a new approach using the cosmic peculiar velocity field to characterize the morphology and size of large scale structures in the local Universe. The algorithm developed uses the three-dimensional peculiar velocity field to compute flow lines, or streamlines. The local Universe is then partitioned into volumes corresponding to gravitational basins, also called watersheds, among the different end-points of the velocity flow lines. This new methodology is first tested on numerical cosmological simulations, used as benchmark for the method, and then applied to the {\it Cosmic-Flows} project observational data in order to to pay particular attention to the nearby superclusters including ours. More extensive tests on both simulated and observational data will be discussed in an accompanying paper.

We present HI 21 cm spectroscopy from the GBT for the host galaxies of 31 nearby AGNs with direct M$_{\textrm{BH}}$ measurements from reverberation mapping. These are the first published HI detections for 12 galaxies, and the spectral quality is generally an improvement over archival data for the remainder of the sample. We present measurements of emission-line fluxes, velocity widths, and recessional velocities from which we derive HI mass, total gas mass, and redshifts. Combining M$_{\textrm{GAS}}$ with constraints on M$_{\textrm{STARS}}$ allows exploration of the baryonic content of these galaxies. We find a typical M$_{\textrm{GAS}}$/M$_{\textrm{STARS}}$ fraction of 10%, with a few reaching $\sim$30-50%. We also examined several relationships between M$_{\textrm{STARS}}$, M$_{\textrm{GAS}}$, M$_{\textrm{BH}}$, baryonic mass, and morphological type. We find a weak preference for galaxies with larger M$_{\textrm{GAS}}$ to host more massive black holes. We also find gas-to-stellar fractions to weakly correlate with later types in unbarred spirals, with an approximately constant fraction for barred spirals. Consistent with previous studies, we find declining M$_{\textrm{GAS}}$/M$_{\textrm{STARS}}$ with increasing M$_{\textrm{STARS}}$, with a slope suggesting the gas reservoirs have been replenished. Finally, we find a clear relationship for M$_{\textrm{BH}}$-M$_{\textrm{BARY}}$ with a similar slope as M$_{\textrm{BH}}$-M$_{\textrm{STARS}}$ reported by Bentz & Manne-Nicholas (2018). The dwarf Seyfert NGC 4395 appears to follow this relationship as well, even though it has a significantly higher gas fraction and smaller M$_{\textrm{BH}}$ than the remainder of our sample.

Cosmicflows-3 distances and inferred peculiar velocities of galaxies have permitted the reconstruction of the structure of over and under densities within the volume extending to 0.05c. This study focuses on the under dense regions, particularly the Local Void that lies largely in the zone of obscuration and consequently has received limited attention. Major over dense structures that bound the Local Void are the Perseus-Pisces and Norma-Pavo-Indus filaments separated by 8,500 km/s. The void network of the universe is interconnected and void passages are found from the Local Void to the adjacent very large Hercules and Sculptor voids. Minor filaments course through voids. A particularly interesting example connects the Virgo and Perseus clusters, with several substantial galaxies found along the chain in the depths of the Local Void. The Local Void has a substantial dynamical effect, causing a deviant motion of the Local Group of 200-250 km/s. The combined perturbations due to repulsion from the Local Void and attraction toward the Virgo Cluster account for about 50% of the motion of the Local Group in the rest frame given by the cosmic microwave background.

(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.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

This article explores three usual estimators, noted as $v_{12}$ of the pairwise velocity, $\psi_1$ and $\psi_2$ of the observed two-point galaxy peculiar velocity correlation functions. These estimators are tested on mock samples of {\it Cosmicflows-3} dataset \citep{Tully:2016aa} , derived from a numerical cosmological simulation, and also on a number of constrained realizations of this dataset. Observational measurements errors and cosmic variance are taken into consideration in the study. The result is a local measurement of $f \sigma_8 = 0.43 \left( \pm 0.03 \right)_\mathrm{obs} \left( \pm 0.11 \right)_\mathrm{cosmic}$ out to $z=0.05$, in support of a $\Lambda$CDM cosmology.

In the upcoming decade cadenced wide-field imaging surveys will increase the number of identified $z<0.3$ Type~Ia supernovae (SNe~Ia) from the hundreds to the hundreds of thousands. The increase in the number density and solid-angle coverage of SNe~Ia, in parallel with improvements in the standardization of their absolute magnitudes, now make them competitive probes of the growth of structure and hence of gravity. The peculiar velocity power spectrum is sensitive to the growth index $\gamma$, which captures the effect of gravity on the linear growth of structure through the relation $f=\Omega_M^\gamma$. We present the first projections for the precision in $\gamma$ for a range of realistic SN peculiar-velocity survey scenarios. In the next decade the peculiar velocities of SNe~Ia in the local $z<0.3$ Universe will provide a measure of $\gamma$ to $\pm 0.01$ precision that can definitively distinguish between General Relativity and leading models of alternative gravity.

A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $\Lambda$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate biases inherent to peculiar velocity analyses, such as the Homogeneous Malmquist bias or the log-normal distribution of peculiar velocities. The statistical uncertainty on the recovered velocity field is about 150 km/s depending on the location, and we study systematics coming from the selection function and calibration of distance indicators. The resulting velocity field and related density fields recover the cosmography of the Local Universe which is presented in an unprecedented volume of universe 10 times larger than previously reached. This methodology open the doors to reconstruction of initial conditions for larger and more accurate constrained cosmological simulations. This work is also preparatory to larger peculiar velocity datasets coming from Wallaby, TAIPAN or LSST.

The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constrained by the standard cosmological model and the observational data. The quasi-linear density field is the geometric mean and variance of the fully non-linear density fields of the simulations. The main nearby clusters (Virgo, Centaurus, Coma), superclusters (Shapley, Perseus-Pisces) and voids (Dipole Repeller) are robustly reconstructed. Galaxies are born biased with respect to the underlying dark matter distribution. Using our quasi-linear framework we demonstrate that the luminosity-weighted density field derived from the 2M++ redshift compilations is non-linearly biased with respect to the matter density field. The bias diminishes in the linear regime.

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration.

The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion events are defined as the moment a halo becomes a substructure, namely when it crosses its host's virial radius. We quantify the cosmic web at each point by the shear tensor and examine where, with respect to its eigenvectors, such accretion events occur in cold ($\Lambda$CDM) and warm (1keV sterile neutrino WDM) dark matter cosmological models. We find that the CDM and WDM subhalos are preferentially accreted along the principal axis of the shear tensor corresponding to the direction of weakest collapse. The beaming strength is modulated by the host and subhalo masses and by the redshift at which the accretion event occurs. Although strongest for the most massive hosts and subhalos at high redshift, the preferential infall is found to be always aligned with the axis of weakest collapse, thus we say that it has universal nature. We compare the strength of beaming in the WDM cosmology with the one found in the $\Lambda$CDM scenario. While the main findings remain the same, the accretion in the WDM model for the most massive host halos appears more beamed than in $\Lambda$CDM cosmology across all the redshifts.

The three-dimensional gravitational velocity field within z~0.1 has been modeled with the Wiener filter methodology applied to the Cosmicflows-3 compilation of galaxy distances. The dominant features are a basin of attraction and two basins of repulsion. The major basin of attraction is an extension of the Shapley concentration of galaxies. One basin of repulsion, the Dipole Repeller, is located near the anti-apex of the cosmic microwave background dipole. The other basin of repulsion is in the proximate direction toward the 'Cold Spot' irregularity in the cosmic microwave background. It has been speculated that a vast void might contribute to the amplitude of the Cold Spot from the integrated Sachs-Wolfe effect.

The network of filaments with embedded clusters surrounding voids seen in maps derived from redshift surveys and reproduced in simulations has been referred to as the cosmic web. A complementary description is provided by considering the shear in the velocity field of galaxies. The eigenvalues of the shear provide information on whether a region is collapsing in three dimensions, the condition for a knot, expanding in three-dimensions, the condition for a void, or in the intermediate condition of a filament or sheet. The structures that are quantitatively defined by the eigenvalues can be approximated by iso-contours that provide a visual representation of the cosmic velocity (V) web. The current application is based on radial peculiar velocities from the Cosmicflows-2 collection of distances. The three-dimensional velocity field is constructed using the Wiener filter methodology in the linear approximation. Eigenvalues of the velocity shear are calculated at each point on a grid. Here, knots and filaments are visualized across a local domain of diameter ~0.1c.

In the standard (LCDM) model of cosmology the universe has emerged out of an early homogeneous and isotropic phase. Structure formation is associated with the growth of density irregularities and peculiar velocities. Our Local Group is moving with respect to the cosmic microwave background (CMB) with a velocity 631+/-20 km s-1 and participates in a bulk flow that extends out to distances of at least 20,000 km s-1. Since the discovery of the CMB dipole, the implicit assumption was that excesses in the abundance of galaxies induce the Local Group motion. Yet, underdense regions push as much as overdensities attract but they are deficient of light and consequently difficult to chart. It was suggested a decade ago that an underdensity in the northern hemisphere roughly 15,000 km s-1 away is a significant actor in the local flow. Here we report on kinematic evidence for such an underdensity. We map the large scale 3D velocity field using a Wiener filter reconstruction from the Cosmicflows-2 dataset of peculiar velocities, and identify the attractors and repellers that dominate the local dynamics. We show here that the local flow is dominated by a single attractor -associated with the Shapley Concentration- and a single previously unidentified repeller. Multipole expansion of the local flow provides further support for the existence and role played by the attractor and repeller. The bulk flow (i.e. dipole moment) is closely (anti)aligned with the repeller at a distance of 16,000+/-4,500 km s-1. The expansion eigenvector of the shear tensor (quadrupole moment) is closely aligned with the Shapley Attractor out to 7,000 km s-1. The close alignment of the local bulk flow with the repeller provides further support for its dominant role in shaping the local flow. This Dipole Repeller is predicted to be associated with a void in the distribution of galaxies.

Cosmography, the study and making of maps of the universe or cosmos, is a field where visual representation benefits from modern three-dimensional visualization techniques and media. At the extragalactic distance scales, visualization is contributing in understanding the complex structure of the local universe, in terms of spatial distribution and flows of galaxies and dark matter. In this paper, we report advances in the field of extragalactic cosmography obtained using the SDvision visualization software in the context of the Cosmicflows Project. Here, multiple visualization techniques are applied to a variety of data products: catalogs of galaxy positions and galaxy peculiar velocities, reconstructed velocity field, density field, gravitational potential field, velocity shear tensor viewed in terms of its eigenvalues and eigenvectors, envelope surfaces enclosing basins of attraction. These visualizations, implemented as high-resolution images, videos, and interactive viewers, have contributed to a number of studies: the cosmography of the local part of the universe, the nature of the Great Attractor, the discovery of the boundaries of our home supercluster of galaxies Laniakea, the mapping of the cosmic web, the study of attractors and repellers.

We use a sample of Milky Way (MW) analogs for which we have stellar and disk gas mass measurements, published measurements of halo gas masses of the MW and of similar galaxies, and the well-established value of the cosmological baryon fraction to place a lower bound on the mass of the Galaxy of $7.7\times10^{11} M_\odot$ and estimate that the mass is likely to be $\ge 1.2\times10^{12}$ M$_\odot$. Although most dynamical analyses yield measurements consistent with these results, several recent studies have advocated for a total mass well below $10^{12} M_\odot$. We reject such low mass estimates because they imply a Galactic baryon matter fraction significantly above the universal value. Convergence between dynamical mass estimates and those based on the baryonic mass is an important milestone in our understanding of galaxies.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The goodness-of-fit (GoF) of the Cosmicflows-2 (CF2) database of peculiar velocities with the LCDM standard model of cosmology is presented. Standard application of the Chi^2 statistics of the full database, of its 4,838 data points, is hampered by the small scale nonlinear dynamics which is not accounted for by the (linear regime) velocity power spectrum. The bulk velocity constitutes a highly compressed representation of the data which filters out the small scales non-linear modes. Hence the statistics of the bulk flow provides an efficient tool for assessing the GoF of the data given a model. The particular approach introduced here is to use the (spherical top-hat window) bulk velocity extracted from the Wiener filter reconstruction of the 3D velocity field as a linear low pass filtered highly compressed representation of the CF2 data. An ensemble 2250 random linear realizations of the WMAP/LCDM model has been used to calculate the bulk velocity auto-covariance matrix. We find that the CF2 data is consistent with the WMAP/LCDM model to better than the 2 sigma confidence limits. This provides a further validation that the CF2 database is consistent with the standard model of cosmology.

The Cosmicflows database of galaxy distances that in the 2nd edition contained 8,188 entries is now expanded to 17,669 entries. The major additions are 2,257 distances that we have derived from the correlation between galaxy rotation and luminosity with photometry at 3.6 microns obtained with Spitzer Space Telescope and 8,885 distances based on the Fundamental Plane methodology from the 6dFGS collaboration. There are minor augmentations to the Tip of the Red Giant Branch and Type Ia supernova compilations. A zero point calibration of the supernova luminosities give a value for the Hubble Constant of 76.2 +-3.4(r) +-2.7(s) km/s/Mpc. Alternatively, a restriction on the peculiar velocity monopole term representing global infall/outflow implies H_0 = 75 +-2 km/s/Mpc.

Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in order to measure the compatibility of those estimates with LambdaCDM. This is achieved using an ensemble of local group (LG) look-alikes collected from a set of Con- strained Simulations (CSs) of the local Universe, and a standard unconstrained LambdaCDM. The latter allows us to build a control sample of LG-like pairs and to single out the influence of the environment in our results. We find that neither estimate is at odds with LambdaCDM; how- ever, whereas CSs favour higher values of vtan , the reverse is true for estimates based on LG samples gathered from unconstrained simulations, overlooking the environmental element

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

Using 3.6 and 4.5$\mu$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

Superclusters of galaxies can be defined kinematically from local evaluations of the velocity shear tensor. The location where the smallest eigenvalue of the shear is positive and maximal defines the center of a basin of attraction. Velocity and density fields are reconstructed with Wiener Filter techniques. Local velocities due to the density field in a restricted region can be separated from external tidal flows, permitting the identification of boundaries separating inward flows toward a basin of attraction and outward flows. This methodology was used to define the Laniakea Supercluster that includes the Milky Way. Large adjacent structures include Perseus-Pisces, Coma, Hercules, and Shapley but current kinematic data are insufficient to capture their full domains. However there is a small region trapped between Laniakea, Perseus-Pisces, and Coma that is close enough to be reliably characterized and that satisfies the kinematic definition of a supercluster. Because of its shape, it is given the name the Arrowhead Supercluster. This entity does not contain any major clusters. A characteristic dimension is ~25 Mpc and the contained mass is only ~10^15 Msun.

Recent observational studies have demonstrated that the majority of satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard ({\Lambda}CDM) model of cosmology if not even in contradiction to its generic predictions. Guided by {\Lambda}CDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in general agreement with predictions of the {\Lambda}CDM model.

The cosmic web that characterizes the large-scale structure of the Universe can be quantified by a variety of methods. For example, large redshift surveys can be used in combination with point process algorithms to extract long curvilinear filaments in the galaxy distribution. Alternatively, given a full 3D reconstruction of the velocity field, kinematic techniques can be used to decompose the web into voids, sheets, filaments and knots. In this paper we look at how two such algorithms - the Bisous model and the velocity shear web - compare with each other in the local Universe (within 100 Mpc), finding good agreement. This is both remarkable and comforting, given that the two methods are radically different in ideology and applied to completely independent and different data sets. Unsurprisingly, the methods are in better agreement when applied to unbiased and complete data sets, like cosmological simulations, than when applied to observational samples. We conclude that more observational data is needed to improve on these methods, but that both methods are most likely properly tracing the underlying distribution of matter in the Universe.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/).

Full sky surveys of peculiar velocity are arguably the best way to map the large scale structure out to distances of a few times 100 Mpc/h. Using the largest and most accurate ever catalog of galaxy peculiar velocities "Cosmicflows-2", the large scale structure has been reconstructed by means of the Wiener filter and constrained realizations assuming as a Bayesian prior model the LCDM model with the WMAP inferred cosmological parameters. The present paper focuses on studying the bulk flow of the local flow field, defined as the mean velocity of top-hat spheres with radii ranging out to R=500 Mpc/h. The estimated large scale structures, in general, and the bulk flow, in particular, are determined by the tension between the observational data and the assumed prior model. A prerequisite for such an analysis is the requirement that the estimated bulk flow is consistent with the prior model. Such a consistency is found here. At R=50(150) Mpc/h the estimated bulk velocity is 250+/-21 (239+/-38) km/s. The corresponding cosmic variance at these radii is 126(60)km/s, which implies that these estimated bulk flows are dominated by the data and not by the assumed prior model. The estimated bulk velocity is dominated by the data out to R~200 Mpc/h, where the cosmic variance on the individual Supergalactic Cartesian components (of the r.m.s. values) exceeds the variance of the Constrained Realizations by at least a factor of 2. The supergalactic SGX and SGY components of the CMB dipole velocity are recovered by the Wiener filter velocity field down to a very few km/s. The SGZ component of the estimated velocity, the one that is most affected by the Zone of Avoidance, is off by 126 km/s (an almost 2 sigma discrepancy).

We identify some of the most HI massive and fastest rotating disk galaxies in the local universe with the aim of probing the processes that drive the formation of these extreme disk galaxies. By combining data from the Cosmic Flows project, which has consistently reanalyzed archival galaxy HI profiles, and 3.6$\mu$m photometry obtained with the Spitzer Space Telescope, with which we can measure stellar mass, we use the baryonic Tully-Fisher (BTF) relationship to explore whether these massive galaxies are distinct. We discuss several results, but the most striking is the systematic offset of the HI-massive sample above the BTF. These galaxies have both more gas and more stars in their disks than the typical disk galaxy of similar rotational velocity. The "condensed" baryon fraction, $f_C$, the fraction of the baryons in a dark matter halo that settle either as cold gas or stars into the disk, is twice as high in the HI-massive sample than typical, and almost reaches the universal baryon fraction in some cases, suggesting that the most extreme of these galaxies have little in the way of a hot baryonic component or cold baryons distributed well outside the disk. In contrast, the star formation efficiency, measured as the ratio of the mass in stars to that in both stars and gas, shows no difference between the HI-massive sample and the typical disk galaxies. We conclude that the star formation efficiency is driven by an internal, self-regulating process, while $f_C$ is affected by external factors. We also found that the most massive HI detected galaxies are located preferentially in filaments. We present the first evidence of an environmental effect on galaxy evolution using a dynamical definition of a filament.

The Spitzer Survey of Stellar Structure in Galaxies (S4G) is the largest available database of deep, homogeneous middle-infrared (mid-IR) images of galaxies of all types. The survey, which includes 2352 nearby galaxies, reveals galaxy morphology only minimally affected by interstellar extinction. This paper presents an atlas and classifications of S4G galaxies in the Comprehensive de Vaucouleurs revised Hubble-Sandage (CVRHS) system. The CVRHS system follows the precepts of classical de Vaucouleurs (1959) morphology, modified to include recognition of other features such as inner, outer, and nuclear lenses, nuclear rings, bars, and disks, spheroidal galaxies, X patterns and box/peanut structures, OLR subclass outer rings and pseudorings, bar ansae and barlenses, parallel sequence late-types, thick disks, and embedded disks in 3D early-type systems. We show that our CVRHS classifications are internally consistent, and that nearly half of the S4G sample consists of extreme late-type systems (mostly bulgeless, pure disk galaxies) in the range Scd-Im. The most common family classification for mid-IR types S0/a to Sc is SA while that for types Scd to Sm is SB. The bars in these two type domains are very different in mid-IR structure and morphology. This paper examines the bar, ring, and type classification fractions in the sample, and also includes several montages of images highlighting the various kinds of "stellar structures" seen in mid-IR galaxy morphology.

We introduce SoFiA, a flexible software application for the detection and parameterization of sources in 3D spectral-line datasets. SoFiA combines for the first time in a single piece of software a set of new source-finding and parameterization algorithms developed on the way to future HI surveys with ASKAP (WALLABY, DINGO) and APERTIF. It is designed to enable the general use of these new algorithms by the community on a broad range of datasets. The key advantages of SoFiA are the ability to: search for line emission on multiple scales to detect 3D sources in a complete and reliable way, taking into account noise level variations and the presence of artefacts in a data cube; estimate the reliability of individual detections; look for signal in arbitrarily large data cubes using a catalogue of 3D coordinates as a prior; provide a wide range of source parameters and output products which facilitate further analysis by the user. We highlight the modularity of SoFiA, which makes it a flexible package allowing users to select and apply only the algorithms useful for their data and science questions. This modularity makes it also possible to easily expand SoFiA in order to include additional methods as they become available. The full SoFiA distribution, including a dedicated graphical user interface, is publicly available for download.

Cosmic Flows is an international multi-element project with the goal to map motions of galaxies in the Local Universe. Kinematic information from observations in the radio HI line and photometry at optical or near-infrared bands are acquired to derive the large majority of distances that are obtained through the luminosity-linewidth or Tully-Fisher relation. This paper gathers additional observational radio data, frequently unpublished, retrieved from the archives of Green Bank, Parkes and Arecibo telescopes. Extracted HI profiles are consistently processed to produce linewidth measurements. Our current "All-Digital HI Catalog" contains a total of 20,343 HI spectra for 17,738 galaxies with 14,802 galaxies with accurate linewidth measurement useful for Tully-Fisher galaxy distances. This addition of 4,117 new measurements represents an augmentation of 34\% compared to our last release.

Galaxies congregate in clusters and along filaments, and are missing from large regions referred to as voids. These structures are seen in maps derived from spectroscopic surveys that reveal networks of structure that are interconnected with no clear boundaries. Extended regions with a high concentration of galaxies are called 'superclusters', although this term is not precise. There is, however, another way to analyse the structure. If the distance to each galaxy from Earth is directly measured, then the peculiar velocity can be derived from the subtraction of the mean cosmic expansion, the product of distance times the Hubble constant, from observed velocity. The peculiar velocity is the line-of-sight departure from the cosmic expansion and arises from gravitational perturbations; a map of peculiar velocities can be translated into a map of the distribution of matter. Here we report a map of structure made using a catalogue of peculiar velocities. We find locations where peculiar velocity flows diverge, as water does at watershed divides, and we trace the surface of divergent points that surrounds us. Within the volume enclosed by this surface, the motions of galaxies are inward after removal of the mean cosmic expansion and long range flows. We define a supercluster to be the volume within such a surface, and so we are defining the extent of our home supercluster, which we call Laniakea. A key component of this paper is an accompanying movie that can be viewed (also in 3D) and downloaded at http://irfu.cea.fr/laniakea or http://vimeo.com/pomarede/laniakea

We present the compilation catalogue of redshift-independent distances included in the HyperLEDA database. It is actively maintained to be up-to-date, and the current version counts 6640 distance measurements for 2335 galaxies compiled from 430 published articles. Each individual series is recalibrated onto a common distance scale based on a carefully selected set of high-quality measurements. This information together with data on HI line-width, central velocity dispersion, magnitudes, diameters, and redshift is used to derive a homogeneous distance estimate and physical properties of galaxies, such as their absolute magnitudes and intrinsic size.

This paper involves a data release of the observational campaign: Cosmicflows with Spitzer (CFS). Surface photometry of the 1270 galaxies constituting the survey is presented. An additional ~ 400 galaxies from various other Spitzer surveys are also analyzed. CFS complements the Spitzer Survey of Stellar Structure in Galaxies, that provides photometry for an additional 2352 galaxies, by extending observations to low galactic latitudes (|b|<30 degrees). Among these galaxies are calibrators, selected in K band, of the Tully-Fisher relation. The addition of new calibrators demonstrate the robustness of the previously released calibration. Our estimate of the Hubble constant using supernova host galaxies is unchanged, H0 = 75.2 +/- 3.3 km/s/Mpc. Distance-derived radial peculiar velocities, for the 1935 galaxies with all the available parameters, will be incorporated into a new data release of the Cosmicflows project. The size of the previous catalog will be increased by 20%, including spatial regions close to the Zone of Avoidance.

In order to explore local large-scale structures and velocity fields, accurate galaxy distance measures are needed. We now extend the well-tested recipe for calibrating the correlation between galaxy rotation rates and luminosities -- capable of providing such distance measures -- to the all-sky, space-based imaging data from the Wide-field Infrared Survey Explorer (WISE) W1 ($3.4\mu$m) and W2 ($4.6\mu$m) filters. We find a linewidth to absolute magnitude correlation (known as the Tully-Fisher Relation, TFR) of $\mathcal{M}^{b,i,k,a}_{W1} = -20.35 - 9.56 (\log W^i_{mx} - 2.5)$ (0.54 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.76 - 9.74 (\log W^i_{mx} - 2.5)$ (0.56 magnitudes rms) from 310 galaxies in 13 clusters. We update the I-band TFR using a sample 9% larger than in Tully & Courtois (2012). We derive $\mathcal{M}^{b,i,k}_I = -21.34 - 8.95 (\log W^i_{mx} - 2.5)$ (0.46 magnitudes rms). The WISE TFRs show evidence of curvature. Quadratic fits give $\mathcal{M}^{b,i,k,a}_{W1} = -20.48 - 8.36 (\log W^i_{mx} - 2.5) + 3.60 (\log W^i_{mx} - 2.5)^2$ (0.52 magnitudes rms) and $\mathcal{M}^{b,i,k,a}_{W2} = -19.91 - 8.40 (\log W^i_{mx} - 2.5) + 4.32 (\log W^i_{mx} - 2.5)^2$ (0.55 magnitudes rms). We apply an I-band -- WISE color correction to lower the scatter and derive $\mathcal{M}_{C_{W1}} = -20.22 - 9.12 (\log W^i_{mx} - 2.5)$ and $\mathcal{M}_{C_{W2}} = -19.63 - 9.11 (\log W^i_{mx} - 2.5)$ (both 0.46 magnitudes rms). Using our three independent TFRs (W1 curved, W2 curved and I-band), we calibrate the UNION2 supernova Type Ia sample distance scale and derive $H_0 = 74.4 \pm 1.4$(stat) $\pm\ 2.4$(sys) kms$^{-1}$ Mpc$^{-1}$ with 4% total error.

We combine data from the Spitzer Survey for Stellar Structure in Galaxies (S$^4$G), a recently calibrated empirical stellar mass estimator from Eskew et al., and an extensive database of HI spectral line profiles to examine the baryonic Tully-Fisher (BTF) relation. We find 1) that the BTF has lower scatter than the classic Tully-Fisher (TF) relation and is better described as a linear relationship, confirming similar previous results, 2) that the inclusion of a radial scale in the BTF decreases the scatter but only modestly, as seen previously for the TF relation, and 3) that the slope of the BTF, which we find to be $3.5\pm 0.2$ ($\Delta$ log $M_{baryon}/\Delta$ log $v_c$), implies that on average a nearly constant fraction ($\sim 0.4$) of all baryons expected to be in a halo are "condensed" onto the central region of rotationally supported galaxies. The condensed baryon fraction, $M_{baryon}/M_{total}$, is, to our measurement precision, nearly independent of galaxy circular velocity (our sample spans circular velocities, $v_c$, between 60 and 250 km s$^{-1}$, but is extended to $v_c\sim 10$ km s$^{-1}$ using data from the literature). The observed galaxy-to-galaxy scatter in this fraction is generally $\le$ a factor of 2 despite fairly liberal selection criteria. These results imply that cooling and heating processes, such as cold vs. hot accretion, mass loss due to stellar winds, and AGN driven feedback, to the degree that they affect the global galactic properties involved in the BTF, are independent of halo mass for galaxies with $10 < v_c < 250$ km/s and typically introduce no more than a factor of two range in the resulting $M_{baryon}/M_{total}$. Recent simulations by Aumer et al. of a small sample of disk galaxies are in excellent agreement with our data, suggesting that current simulations are capable of reproducing the global properties of individual disk galaxies.

Peculiar velocities, obtained from direct distance measurements, are data of choice to achieve constrained simulations of the Local Universe reliable down to a scale of a few Megaparsecs. Unlike redshift surveys, peculiar velocities are direct tracers of the underlying gravitational field as they trace both baryonic and dark matter. This paper presents the first attempt to use solely observational peculiar velocities to constrain cosmological simulations of the nearby universe. In order to set up Initial Conditions, a Reverse Zel'dovich Approximation (RZA) is used to displace constraints from their positions at z=0 to their precursors' locations at higher redshifts. An additional new feature replaces original observed radial peculiar velocity vectors by their full 3D reconstructions provided by the Wiener-Filter (WF) estimator. Subsequently, the Constrained Realization of Gaussian fields technique (CR) is applied to build various realizations of the Initial Conditions. The WF/RZA/CR method is first tested on realistic mock catalogs built from a reference simulation similar to the Local Universe. These mocks include errors on peculiar velocities, on data-point positions and a large continuous zone devoid of data in order to mimic galactic extinction. Large scale structures are recovered with a typical accuracy of 5 Megaparsecs/h in position, the best realizations reaching a 2-3 Mpc/h precision, the limit imposed by the RZA linear theory. Then, the method is applied to the first observational radial peculiar velocity catalog of the project Cosmicflows. This paper is a proof of concept that the WF/RZA/CR method can be applied to observational peculiar velocities to successfully build constrained Initial Conditions.

Cosmicflows-2 is a compilation of distances and peculiar velocities for over 8000 galaxies. Numerically the largest contributions come from the luminosity-linewidth correlation for spirals, the TFR, and the related Fundamental Plane relation for E/S0 systems, but over 1000 distances are contributed by methods that provide more accurate individual distances: Cepheid, Tip of the Red Giant Branch, Surface Brightness Fluctuation, SNIa, and several miscellaneous but accurate procedures. Our collaboration is making important contributions to two of these inputs: Tip of the Red Giant Branch and TFR. A large body of new distance material is presented. In addition, an effort is made to assure that all the contributions, our own and those from the literature, are on the same scale. Overall, the distances are found to be compatible with a Hubble Constant H_0 = 74.4 +-3.0 km/s/Mpc. The great interest going forward with this data set will be with velocity field studies. Cosmicflows-2 is characterized by a great density and high accuracy of distance measures locally, falling to sparse and coarse sampling extending to z=0.1.

We use the FUV fluxes measured with the GALEX to study the star formation properties of galaxies collected in the "Local Orphan Galaxies" catalog (LOG). Among 517 LOG galaxies having radial velocities V(LG) < 3500 km/s and Galactic latitudes |b|> 15 degr, 428 objects have been detected in FUV. We briefly discuss some scaling relations between the specific star formation rate (SSFR) and stellar mass, HI-mass, morphology, and surface brightness of galaxies situated in extremely low density regions of the Local Supercluster. Our sample is populated with predominantly late-type, gas-rich objects with the median morphological type of Sdm. Only 5% of LOG galaxies are classified as early types: E, S0, S0/a, however, they systematically differ from normal E and S0 galaxies by lower luminosity and presence of gas and dust. We find that almost all galaxies in our sample have their SSFR below 0.4 [Gyr^{-1}]. This limit is also true even for a sample of 260 active star-burst Markarian galaxies situated in the same volume. The existence of such a quasi-Eddington limit for galaxies seems to be a key factor which characterizes the transformation of gas into stars at the current epoch.

The large scale structure of the universe is a complex web of clusters, filaments, and voids. Its properties are informed by galaxy redshift surveys and measurements of peculiar velocities. Wiener Filter reconstructions recover three-dimensional velocity and total density fields. The richness of the elements of our neighborhood are revealed with sophisticated visualization tools. A key component of this paper is an accompanying movie that can be viewed and downloaded at : http://irfu.cea.fr/cosmography or http://vimeo.com/pomarede/cosmography .The ability to translate and zoom helps the viewer follow structures in three dimensions and grasp the relationships between features on different scales while retaining a sense of orientation. The ability to dissolve between scenes provides a technique for comparing different information, for example, the observed distribution of galaxies, smoothed representations of the distribution accounting for selection effects, observed peculiar velocities, smoothed and modeled representations of those velocities, and inferred underlying density fields. The agreement between the large scale structure seen in redshift surveys and that inferred from reconstructions based on the radial peculiar velocities of galaxies strongly supports the standard model of cosmology where structure forms from gravitational instabilities and galaxies form at the bottom of potential wells.

We report on measurements of the disk central surface brightnesses (mu0) at 3.6 microns for 438 galaxies selected by distance and absolute magnitude cutoffs from the 2350+ galaxies in the Spitzer Survey of Stellar Structure in Galaxies (S4G), one of the largest and deepest homogeneous mid-infrared datasets of nearby galaxies. Our sample contains nearly 3 times more galaxies than the most recent study of the mu0 distribution. We demonstrate that there is a bimodality in the distribution of mu0. Between the low and high surface brightness galaxy regimes there is a lack of intermediate surface brightness galaxies. Caveats invoked in the literature from small number statistics to the knowledge of the environmental influences, and possible biases from low signal to noise data or corrections for galaxy inclination are investigated. Analyses show that the bimodal distribution of mu0 cannot be due to any of these biases or statistical fluctuations. It is highly probable that galaxies settle in two stable modes: a dark matter dominated mode where the dark matter dominates at all radii - this gives birth to low surface brightness galaxies - and a baryonic matter dominated mode where the baryons dominate the dark matter in the central parts - this gives rise to the high surface brightness disks. The lack of intermediate surface brightness objects suggests that galaxies avoid (staying in) a mode where dark matter and baryons are co-dominant in the central parts of galaxies.

Distance measures on a coherent scale around the sky are required to address the outstanding cosmological problems of the Hubble Constant and of departures from the mean cosmic flow. The correlation between galaxy luminosities and rotation rates can be used to determine distances to many thousands of galaxies in a wide range of environments potentially out to 200 Mpc. Mid-infrared (3.6 microns) photometry with the Spitzer Space Telescope is particularly valuable as the source of the luminosities because it provides products of uniform quality across the sky. From a perch above the atmosphere, essentially the total magnitude of targets can be registered in exposures of a few minutes. Extinction is minimal and the flux is dominated by the light from old stars which is expected to correlate with the mass of the targets. In spite of the superior photometry, the correlation between mid-infrared luminosities and rotation rates extracted from neutral hydrogen profiles is slightly degraded from the correlation found with I band luminosities. A color correction recovers a correlation that provides comparable accuracy to that available at I band (~20% 1sigma in an individual distance) while retaining the advantages identified above. Without the color correction the relation between linewidth and [3.6] magnitudes is M^{b,i,k,a}_{[3.6]} = -20.34 - 9.74 (log W_{mx}^{i} -2.5). This description is found with a sample of 213 galaxies in 13 clusters that define the slope and 26 galaxies with Cepheid or tip of the red giant branch distances that define the zero point. A color corrected parameter M_{C_{[3.6]}} is constructed that has reduced scatter: M_{C_{[3.6]}} = -20.34 - 9.13 (log W_{mx}^{i} -2.5). Consideration of the 7 calibration clusters beyond 50 Mpc, outside the domain of obvious peculiar velocities, provides a preliminary Hubble Constant estimate of H_0=74+/-5 km/s/Mpc.

The Reverse Zeldovich Approximation (RZA) is a reconstruction method which allows to estimate the cosmic displacement field from galaxy peculiar velocity data and to constrain initial conditions for cosmological simulations of the Local Universe. In this paper, we investigate the effect of different observational errors on the reconstruction quality of this method. For this, we build a set of mock catalogues from a cosmological simulation, varying different error sources like the galaxy distance measurement error (0 - 20%), the sparseness of the data points, and the maximum catalogue radius (3000 - 6000 km/s). We perform the RZA reconstruction of the initial conditions on these mock catalogues and compare with the actual initial conditions of the simulation. We also investigate the impact of the fact that only the radial part of the peculiar velocity is observationally accessible. We find that the sparseness of a dataset has the highest detrimental effect on RZA reconstruction quality. Observational distance errors also have a significant influence, but it is possible to compensate this relatively well with Wiener Filter reconstruction. We also investigate the effect of different object selection criteria and find that distance catalogues distributed randomly and homogeneously across the sky (such as spiral galaxies selected for the Tully-Fisher method) allow for a higher reconstruction quality than if when data is preferentially drawn from massive objects or dense environments (such as elliptical galaxies). We find that the error of estimating the initial conditions with RZA is always dominated by the inherent non-linearity of data observed at z=0 rather than by the combined effect of the observational errors. Even an extremely sparse dataset with high observational errors still leads to a good reconstruction of the initial conditions on a scale of about 5 Mpc/h.

We propose a new method to recover the cosmological initial conditions of the presently observed galaxy distribution, which can serve to run constrained simulations of the Local Universe. Our method, the Reverse Zeldovich Approximation (RZA), can be applied to radial galaxy peculiar velocity data and extends the previously used Constrained Realizations (CR) method by adding a Lagrangian reconstruction step. The RZA method consists of applying the Zeldovich approximation in reverse to galaxy peculiar velocities to estimate the cosmic displacement field and the initial linear matter distribution from which the present-day Local Universe evolved.We test our method with a mock survey taken from a cosmological simulation. We show that the halo peculiar velocities at z = 0 are close to the linear prediction of the Zeldovich approximation, if a grouping is applied to the data to remove virial motions. We find that the addition of RZA to the CR method significantly improves the reconstruction of the initial conditions. The RZA is able to recover the correct initial positions of the velocity tracers with a median error of only 1.36 Mpc/h in our test simulation. For realistic sparse and noisy data, this median increases to 5 Mpc/h. This is a significant improvement over the previous approach of neglecting the displacement field, which introduces errors on a scale of 10 Mpc/h or even higher. Applying the RZA method to the upcoming high-quality observational peculiar velocity catalogues will generate much more precise constrained simulations of the Local Universe.

In previous works we proposed the Reverse Zeldovich Approximation (RZA) method, which can be used to estimate the cosmological initial conditions underlying the galaxy distribution in the Local Universe using peculiar velocity data. In this paper, we apply the technique to run constrained cosmological simulations from the RZA-reconstructed initial conditions, designed to reproduce the large-scale structure of the Local Universe. We test the method with mock peculiar velocity catalogues extracted from a reference simulation. We first reconstruct the initial conditions of this reference simulation using the mock data, and then run the reconstructed initial conditions forward in time until z=0. We compare the resulting constrained simulations with the original simulation at z=0 to test the accuracy of this method. We also compare them with constrained simulations run from the mock data without the addition of RZA, i.e. using only the currently established constrained realizations (CR) method. Our re-simulations are able to correctly recover the evolution of the large-scale structure underlying the data. The results show that the addition of RZA to the CR method significantly improves both the reconstruction of the initial conditions and the accuracy of the obtained constrained resimulations. Haloes from the original simulation are recovered in the re-simulations with an average accuracy of about 2 Mpc/h on their position and a factor of 2 in mass, down to haloes with a mass of approx 10^14 M_odot/h. In comparison, without RZA the re-simulations recover only the most massive haloes with masses of about 5x10^14 M_odot/h and higher, and with a systematic shift on their position of about approx 10 Mpc/h due to the cosmic displacement field. We show that with the additional Lagrangian reconstruction step introduced by the RZA, this shift can be removed.

A nearby friable cloud in Ursa Majoris contains 270 galaxies with radial velocities 500 < VLG < 1500 km s^-1 inside the area of RA= [11h; 13h] and DEC= [+40deg; +60deg]. At present, 97 galaxies of them have individual distance estimates. We use these data to clarify the structure and kinematics of the UMa complex. According to Makarov & Karachentsev (2011), most of the UMa galaxies belong to seven bound groups, which have the following median parameters: velocity dispersion of 58 km s^-1, harmonic projected radius of 300 kpc, virial mass of 2.10^12 Msol, and virial- mass-to-K-band-luminosity of 27Msol/Lsol. Almost a half of the UMa cloud population are gas-rich dwarfs (Ir, Im, BCD) with active star formation seen in the GALEX UV-survey. The UMa groups reside within 15-19 Mpc from us, being just at the same distance as Virgo cluster. The total virial mass of the UMa groups is 4.10^13 Msol, yielding the average density of dark matter in the UMa cloud to be Omega_m = 0.08, i.e. a factor three lower than the cosmic average. This is despite the fact that the UMa cloud resides in a region of the Universe that is an apparent overdensity. A possible explanation for this is that most mass in the Universe lies in the empty space between clusters. Herewith, the mean distances and velocities of the UMa groups follow nearly undisturbed Hubble flow without a sign of the 'Z-wave" effect caused by infall toward a massive attractor. This constrains the total amount of dark matter between the UMa groups within the cloud volume.

A visualization of three-dimensional structures and cosmic flows is presented using information from the Extragalactic Distance Database V8k redshift catalog and peculiar velocities from the Cosmicflows-1 survey. Structures within a volume bounded at 8000 km/s on the cardinal Supergalactic axes are explored in terms of both the display of the positions of the 30124 galaxies of the catalog and its reconstructed luminosity density field, corrected to account for growing incompleteness with distance. Cosmography of the Local Universe is discussed with the intent to identify the most prominent structures, including voids, galaxy clusters, filaments and walls. The mapping also benefits from precise distance measures provided through the Cosmicflows-1 observational program. Three-dimensional visualizations of the coherent flows of galaxies in the nearby universe are presented using recent results obtained on the reconstruction of cosmic flows with the Wiener Filter approach. The three major components of the Milky Way motion, namely the expulsion from the Local Void, the infall toward the Virgo Cluster, and the bulk flow of the historic Local Supercluster toward the Great Attractor are illustrated using different visualization techniques and analyzed in the light of the cosmography derived from the V8k redshift and Cosmicflows-1 distance catalogs.

The availability of photometric imaging of several thousand galaxies with the Spitzer Space Telescope enables a mid-infrared calibration of the correlation between luminosity and rotation in spiral galaxies. The most important advantage of the new calibration in the 3.6 micron band, IRAC ch.1, is photometric consistency across the entire sky. Additional advantages are minimal obscuration, observations of flux dominated by old stars, and sensitivity to low surface brightness levels due to favorable backgrounds. Through Spitzer cycle 7 roughly 3000 galaxies had been observed and images of these are available at the Spitzer archive. In cycle 8 a program called Cosmic Flows with Spitzer has been initiated that will increase by 1274 the available sample of spiral galaxies with inclinations greater than 45 degrees from face-on suitable for distance measurements. This paper describes procedures based on the photometry package Archangel that are being employed to analyze both the archival and the new data in a uniform way. We give results for 235 galaxies, our calibrator sample for the Tully-Fisher relation. Galaxy magnitudes are determined with uncertainties held below 0.05 mag for normal spiral systems. A subsequent paper will describe the calibration of the [3.6] luminosity-rotation relation.

This paper builds on a calibration of the SNIa absolute distance scale begun with a core of distances based on the correlation between galaxy rotation rates and optical Ic band photometry. This new work extends the calibration through the use of mid-infrared photometry acquired at 3.6 microns with Spitzer Space Telescope. The great virtue of the satellite observations is constancy of the photometry at a level better than 1% across the sky. The new calibration is based on 39 individual galaxies and 8 clusters that have been the sites of well observed SNIa. The new 3.6 micron calibration is not yet as extensively based as the Ic band calibration but is already sufficient to justify a preliminary report. Distances based on the mid-infrared photometry are 2% greater in the mean than reported at Ic band. This difference is only marginally significant. The Ic band result is confirmed with only a small adjustment. Incorporating a 1% decrease in the LMC distance, the present study indicates Ho = 75.2 +/- 3.0 km/s/Mpc.

Cosmic Flows is a program to determine galaxy distances for 30,000 galaxies with systematic errors below 2%, almost ten times the number currently known and a five-fold improvement in systematics. The resultant velocity field will provide input for constrained local universe simulations: CLUES (www.clues-project.org). The observed and the simulated universe are then comparatively studied. This synergy of observations and theory distinguishes the program, and should lead to fundamental discoveries regarding the sources of deviations from the expansion of the universe. Specifically, the program should give a definitive answer to one of the most outstanding unsolved problem in cosmology: the cause of the motion of 630 km/s of our Galaxy manifested in the microwave background dipole. This paper presents current results with particular emphasis on the "great attractor" reconstruction.

In order to measure distances with minimal systematics using the correlation between galaxy luminosities and rotation rates it is necessary to adhere to a strict and tested recipe. We now derive a measure of rotation from a new characterization of the width of a neutral Hydrogen line profile. Additionally, new photometry and zero point calibration data are available. Particularly the introduction of a new linewidth parameter necessitates the reconstruction and absolute calibration of the luminosity-linewidth template. The slope of the new template is set by 267 galaxies in 13 clusters. The zero point is set by 36 galaxies with Cepheid or Tip of the Red Giant Branch distances. Tentatively, we determine H0 = 75 km s-1 Mpc-1. Distances determined using the luminosity-linewidth calibration will contribute to the distance compendium Cosmicflows-2.

This paper presents an analysis of the local peculiar velocity field based on the Wiener Filter reconstruction method. We used our currently available catalog of distance measurements containing 1,797 galaxies within 3000 km/s: Cosmicflows-1. The Wiener Filter method is used to recover the full 3D peculiar velocity field from the observed map of radial velocities and to recover the underlying linear density field. The velocity field within a data zone of 3000 km/s is decomposed into a local component that is generated within the data zone and a tidal one that is generated by the mass distribution outside that zone. The tidal component is characterized by a coherent flow toward the Norma-Hydra-Centaurus (Great Attractor) region while the local component is dominated by a flow toward the Virgo Cluster and away from the Local Void. A detailed analysis shows that the local flow is predominantly governed by the Local Void and the Virgo Cluster plays a lesser role. The analysis procedure was tested against a mock catalog. It is demonstrated that the Wiener Filter accurately recovers the input velocity field of the mock catalog on the scale of the extraction of distances and reasonably recovers the velocity field on significantly larger scales. The Bayesian Wiener Filter reconstruction is carried out within the ?CDM WMAP5 framework. The Wiener Filter reconstruction draws particular attention to the importance of voids in proximity to our neighborhood. The prominent structure of the Local Supercluster is wrapped in a horseshoe collar of under density with the Local Void as a major component.

We notice that nearby galaxies having high negative peculiar velocities are distributed over the sky very inhomogeneously. A part of this anisotropy is caused by the "Local Velocity Anomaly", i.e. by the bulk motion of nearby galaxies away from the Local Void. But a half of the fast-flying objects reside within a small region RA = [11.5h, 13.0h], Dec. = [+20\circ, +40\circ], known as the Coma I cloud. According to Makarov & Karachentsev (2011), this complex contains 8 groups, 5 triplets, 10 pairs and 83 single galaxies with the total mass of 4.7\star10^13M\odot. We use 122 galaxies in the Coma I region with known distances and radial velocities VLG < 3000 km/s to draw the Hubble relation for them. The Hubble diagram shows a Z-shape effect of infall with an amplitude of +200 km/s on the nearby side and -700 km/s on the back side. This phenomena can be understood as the galaxy infall towards a dark attractor with the mass of \sim 2\star10^14M\odot situated at a distance of 15 Mpc from us. The existence of large void between the Coma and Virgo clusters affects probably the Hubble flow around the Coma I also.

Within the "Cosmic Flows" project, I-band photometry of 524 relatively nearby galaxies has been carried out over the course of several years with the University of Hawaii 2.2m telescope and a camera with a 7.5 arcmin field of view. The primary scientific purpose was to provide global magnitudes and inclinations for galaxies for the purpose of measuring distances through the correlation between galaxy luminosities and rotation rates. The accuracy on a total magnitude is 0.08 mag. The observations typically extend to 7-8 exponential disk scalelengths so the data are useful for studies of the structural properties of galaxies.

The neutral hydrogen properties of 1,822 galaxies are being studied with the Green Bank 100m and the Parkes 64m telescopes as part of the 'Cosmic Flows' program. Observed parameters include systemic velocities, profile line widths, and fluxes. The line width information can be combined with optical and infrared photometry to obtain distances. The 1,822 HI observations complement an inventory of archives. All told, HI line width information is available for almost all of five samples: (i) luminosity-line width correlation calibrators, (ii) zero-point calibrators for the supernova Ia scale, (iii) a dense local sample of spiral galaxies with M_{Ks} < -21 within 3,000 km/s, (iv) a sparser sample of 60 $\mu$m selected galaxies within 6,000 km/s that provides all-sky coverage of our extended supercluster complex, and (v) an even sparser sample of flat galaxies, extreme edge-on spirals, extending in a volume out to 12,000 km/s. The HI information for 13,941 galaxies, whether from the archives or acquired as part of the Cosmic Flows observational program, is uniformly re-measured and made available through the Extragalactic Distance Database web site.

The work presented in this paper aims at restricting the input parameter values of the semi-analytical model used in GALICS and MOMAF, so as to derive which parameters influence the most the results, e.g., star formation, feedback and halo recycling efficiencies, etc. Our approach is to proceed empirically: we run lots of simulations and derive the correct ranges of values. The computation time needed is so large, that we need to run on a grid of computers. Hence, we model GALICS and MOMAF execution time and output files size, and run the simulation using a grid middleware: DIET. All the complexity of accessing resources, scheduling simulations and managing data is harnessed by DIET and hidden behind a web portal accessible to the users.

An important component of the Extragalactic Distance Database (EDD) at http://edd.ifa.hawaii.edu is a group of catalogs related to the measurement of HI line profile parameters. One of these is the All Digital HI catalog which contains an amalgam of information from new data and old. The new data results from observations with Arecibo and Parkes telescopes and with the Green Bank Telescope (GBT), including continuing input since the award of the NRAO Cosmic Flows Large Program. The old data has been collected from archives, wherever available, particularly the Cornell University Digital HI Archive, the Nancay Telescope extragalactic HI archive, and the Australia Telescope archive. The catalog currently contains information on ~15,000 profiles relating to ~13,000 galaxies. The channel - flux per channel files, from whatever source, are carried through a common pipeline. The derived parameter of greatest interest is W_m50, the profile width at 50% of the mean flux. After appropriate adjustment, the parameter W_mx is derived, the linewidth which statistically approximates the peak to peak maximum rotation velocity before correction for inclination, 2 V_max sin(i).

A database can be accessed on the web at http://edd.ifa.hawaii.edu that was developed to promote access to information related to galaxy distances. The database has three functional components. First, tables from many literature sources have been gathered and enhanced with links through a distinct galaxy naming convention. Second, comparisons of results both at the levels of parameters and of techniques have begun and are continuing, leading to increasing homogeneity and consistency of distance measurements. Third, new material are presented arising from ongoing observational programs at the University of Hawaii 2.2m telescope, radio telescopes at Green Bank, Arecibo, and Parkes and with Hubble Space Telescope. This new observational material is made available in tandem with related material drawn from archives and passed through common analysis pipelines.

The Tully-Fisher relation is a correlation between the luminosity and the HI 21cm line width in spiral galaxies (LLW relation). It is used to derive galaxy distances in the interval 7 to 100 Mpc. Closer, the Cepheids, TRGB and Surface Brightness Fluctuation methods give a better accuracy. Further, the SNIa are luminous objects still available for distance measurement purposes, though with a dramatically lower density grid of measurements on the sky. Galaxies in clusters are all at the same distance from the observer. Thus the distance of the cluster derived from a large number of galaxies (N) has an error reduced according to the square root of N. However, not all galaxies in a cluster are suitable for the LLW measurement. The selection criteria we use are explained hereafter; the important point being to avoid Malmquist bias and to not introduce any systematics in the distance measurement.

The peculiar velocity of the Local Group of galaxies manifested in the Cosmic Microwave Background dipole is found to decompose into three dominant components. The three components are clearly separated because they arise on distinct spatial scales and are fortuitously almost orthogonal in their influences. The nearest, which is distinguished by a velocity discontinuity at ~7 Mpc, arises from the evacuation of the Local Void. We lie in the Local Sheet that bounds the void. Random motions within the Local Sheet are small. Our Galaxy participates in the bulk motion of the Local Sheet away from the Local Void. The component of our motion on an intermediate scale is attributed to the Virgo Cluster and its surroundings, 17 Mpc away. The third and largest component is an attraction on scales larger than 3000 km/s and centered near the direction of the Centaurus Cluster. The amplitudes of the three components are 259, 185, and 455 km/s, respectively, adding collectively to 631 km/s in the reference frame of the Local Sheet. Taking the nearby influences into account causes the residual attributed to large scales to align with observed concentrations of distant galaxies and reduces somewhat the amplitude of motion attributed to their pull. On small scales, in addition to the motion of our Local Sheet away from the Local Void, the nearest adjacent filament, the Leo Spur, is seen to be moving in a direction that will lead to convergence with our filament. Finally, a good distance to an isolated galaxy within the Local Void reveals that this dwarf system has a motion of at least 230 km/s away from the void center. Given the velocities expected from gravitational instability theory in the standard cosmological paradigm, the distance to the center of the Local Void must be at least 23 Mpc from our position. The Local Void is large!

This paper presents MoLUSC, a new method for generating mock galaxy catalogs from a large scale ($\approx 1000^3$ Mpc$^3$) dark matter simulation, that requires only modest CPU time and memory allocation. The method uses a small-scale ($\approx 256^3$ Mpc$^3$) dark matter simulation on which the \galics semi-analytic code has been run in order to define the transformation from dark matter density to galaxy density transformation using a probabilistic treatment. MoLUSC is then applied to a large-scale dark matter simulation in order to produce a realistic distribution of galaxies and their associated spectra. This permits the fast generation of large-scale mock surveys using relatively low-resolution simulations. We describe various tests which have been conducted to validate the method, and demonstrate a first application to generate a mock Sloan Digital Sky Survey redshift survey.

One way to access the aggregated power of a collection of heterogeneous machines is to use a grid middleware, such as DIET, GridSolve or NINF. It addresses the problem of monitoring the resources, of handling the submissions of jobs and as an example the inherent transfer of input and output data, in place of the user. In this paper we present how to run cosmological simulations using the RAMSES application along with the DIET middleware. We will describe how to write the corresponding DIET client and server. The remainder of the paper is organized as follows: Section 2 presents the DIET middleware. Section 3 describes the RAMSES cosmological software and simulations, and how to interface it with DIET. We show how to write a client and a server in Section 4. Finally, Section 5 presents the experiments realized on Grid'5000, the French Research Grid, and we conclude in Section 6.

The length of the three-dimensional filaments observed in the fourth public data-release of the SDSS is measured using the {\em local skeleton} method. It consists in defining the set of points where the gradient of the smoothed density field is extremal along its isocontours, with some additional constraints on local curvature to probe actual ridges in the galaxy distribution. A good fit to the mean filament length per unit volume, $\cal{L}$, in the SDSS survey is found to be ${\cal{L}}=(52500\pm6500) (L/{\rm Mpc})^{-1.75\pm0.06}\rm{Mpc}/(100 \rm{Mpc})^{3}$ for $8.2 \leq L \leq 16.4$ Mpc, where $L$ is the smoothing length in Mpc. This result, which deviates only slightly, as expected, from the trivial behavior ${\cal{L}} \propto L^{-2}$, is in excellent agreement with a $\Lambda$CDM cosmology, as long as the matter density parameter remains in the range $0.25 < \Omega_{\rm matter} < 0.4$ at one sigma confidence level, considering the universe is flat. These measurements, which are in fact dominated by linear dynamics, are not significantly sensitive to observational biases such as redshift distortion, edge effects, incompleteness, and biasing between the galaxy distribution and the dark matter distribution. Hence it is argued that the local skeleton is a rather promising and discriminating tool for the analysis of filamentary structures in three-dimensional galaxy surveys.

We present distance measurements to 71 high redshift type Ia supernovae discovered during the first year of the 5-year Supernova Legacy Survey (SNLS). These events were detected and their multi-color light-curves measured using the MegaPrime/MegaCam instrument at the Canada-France-Hawaii Telescope (CFHT), by repeatedly imaging four one-square degree fields in four bands. Follow-up spectroscopy was performed at the VLT, Gemini and Keck telescopes to confirm the nature of the supernovae and to measure their redshift. With this data set, we have built a Hubble diagram extending to z=1, with all distance measurements involving at least two bands. Systematic uncertainties are evaluated making use of the multi-band photometry obtained at CFHT. Cosmological fits to this first year SNLS Hubble diagram give the following results : Omega_M = 0.263 +/- 0.042(stat) +/- 0.032(sys) for a flat LambdaCDM model; and w = -1.023 +/- 0.090(stat) +/- 0.054(sys) for a flat cosmology with constant equation of state w when combined with the constraint from the recent Sloan Digital Sky Survey measurement of baryon acoustic oscillations.

We present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

In order to investigate the properties of large-scale structures of galaxies in the universe, we present an analysis of their spatial distribution at z$<$0.033. We used the LEDA extragalactic database containing over 1 million of galaxies covering the all-sky and the SDSS data included in the public release DR1, yielding to a sample of around 134,000 galaxies having a measured redshift in two survey areas representing 690 sq. degrees. The results of the study are 2D, 3D maps and magnitude number counts of galaxies, drawn from B-band samples. Movies and high resolution figures are available on http://www-obs.univ-lyon1.fr/cosmologie/helene.html

This paper presents the first I-band photometric catalog of the brightest galaxies extracted from the Deep Near Infrared Survey of the Southern Sky (DENIS) An automatic galaxy recognition program has been developed to build this provisional catalog. The method is based on a discriminating analysis. The most discriminant parameter to separate galaxies from stars is proved to be the peak intensity of an object divided by its array. Its efficiency is better than 99%. The nominal accuracy for galaxy coordinates calculated with the Guide Star Catalog is about 6 arcseconds. The cross-identification with galaxies available in the Lyon-Meudon Extragalactic DAtabase (LEDA) allows a calibraton of the I-band photometry with the sample of Mathewson et Al. Thus, the catalog contains total I-band magnitude, isophotal diameter, axis ratio, position angle and a rough estimate of the morphological type code for 20260 galaxies. The internal completeness of this catalog reaches magnitude $I_{lim}=14.5$, with a photometric accuracy of $\sim 0.18m$. 25% of the Southern sky has been processed in this study. This quick look analysis allows us to start a radio and spectrographic follow-up long before the end of the survey.