We explore synergies between the space-based Wide-Field Infrared Survey Telescope (WFIRST) and the ground-based Rubin Observatory Legacy Survey of Space and Time (LSST). In particular, we consider a scenario where the currently envisioned survey strategy for WFIRST's High Latitude Survey (HLS), i.e., 2000 square degrees in four narrow photometric bands is altered in favor of a strategy that combines rapid coverage of the LSST area (to full LSST depth) in one band. We find that a 5-month WFIRST survey in the W-band can cover the full LSST survey area providing high-resolution imaging for >95% of the LSST Year 10 gold galaxy sample. We explore a second, more ambitious scenario where WFIRST spends 1.5 years covering the LSST area. For this second scenario we quantify the constraining power on dark energy equation of state parameters from a joint weak lensing and galaxy clustering analysis, and compare it to an LSST-only survey and to the Reference WFIRST HLS survey. Our survey simulations are based on the WFIRST exposure time calculator and redshift distributions from the CANDELS catalog. Our statistical uncertainties account for higher-order correlations of the density field, and we include a wide range of systematic effects, such as uncertainties in shape and redshift measurements, and modeling uncertainties of astrophysical systematics, such as galaxy bias, intrinsic galaxy alignment, and baryonic physics. Assuming the 5-month WFIRST wide scenario, we find a significant increase in constraining power for the joint LSST+WFIRST wide survey compared to LSST Y10 (FoM(Wwide)= 2.4 FoM(LSST)) and compared to LSST+WFIRST HLS (FoM(Wwide)= 5.5 FoM(HLS)).

We simulate the scientific performance of the Wide-Field Infrared Survey Telescope (WFIRST) High Latitude Survey (HLS) on dark energy and modified gravity. The 1.6 year HLS Reference survey is currently envisioned to image 2000 deg$^2$ in multiple bands to a depth of $\sim$26.5 in Y, J, H and to cover the same area with slit-less spectroscopy beyond z=3. The combination of deep, multi-band photometry and deep spectroscopy will allow scientists to measure the growth and geometry of the Universe through a variety of cosmological probes (e.g., weak lensing, galaxy clusters, galaxy clustering, BAO, Type Ia supernova) and, equally, it will allow an exquisite control of observational and astrophysical systematic effects. In this paper we explore multi-probe strategies that can be implemented given WFIRST's instrument capabilities. We model cosmological probes individually and jointly and account for correlated systematics and statistical uncertainties due to the higher order moments of the density field. We explore different levels of observational systematics for the WFIRST survey (photo-z and shear calibration) and ultimately run a joint likelihood analysis in N-dim parameter space. We find that the WFIRST reference survey alone (no external data sets) can achieve a standard dark energy FoM of >300 when including all probes. This assumes no information from external data sets and realistic assumptions for systematics. Our study of the HLS reference survey should be seen as part of a future community driven effort to simulate and optimize the science return of WFIRST.

We present a weak-lensing analysis of X-ray galaxy groups and clusters selected from the XMM-XXL survey using the first-year data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. Our joint weak-lensing and X-ray analysis focuses on 136 spectroscopically confirmed X-ray-selected systems at 0.031 < z < 1.033 detected in the 25sqdeg XXL-N region. We characterize the mass distributions of individual clusters and establish the concentration-mass (c-M) relation for the XXL sample, by accounting for selection bias and statistical effects, and marginalizing over the remaining mass calibration uncertainty. We find the mass-trend parameter of the c-M relation to be \beta = -0.07 \pm 0.28 and the normalization to be c200 = 4.8 \pm 1.0 (stat) \pm 0.8 (syst) at M200=10^{14}Msun/h and z = 0.3. We find no statistical evidence for redshift evolution. Our weak-lensing results are in excellent agreement with dark-matter-only c-M relations calibrated for recent LCDM cosmologies. The level of intrinsic scatter in c200 is constrained as \sigma(\ln[c200]) < 24% (99.7% CL), which is smaller than predicted for the full population of LCDM halos. This is likely caused in part by the X-ray selection bias in terms of the relaxation state. We determine the temperature-mass (Tx-M500) relation for a subset of 105 XXL clusters that have both measured HSC lensing masses and X-ray temperatures. The resulting Tx-M500 relation is consistent with the self-similar prediction. Our Tx-M500 relation agrees with the XXL DR1 results at group scales, but has a slightly steeper mass trend, implying a smaller mass scale in the cluster regime. The overall offset in the Tx-M500 relation is at the $1.5\sigma$ level, corresponding to a mean mass offset of (34\pm 20)%. We also provide bias-corrected, weak-lensing-calibrated M200 and M500 mass estimates of individual XXL clusters based on their measured X-ray temperatures.

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $\xi_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<\theta<56'$ for $\xi_+$ and $28'<\theta<178'$ for $\xi_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega_m=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints.

We analyze the clustering of galaxies in the first public data release of the HSC Subaru Strategic Program. Despite the relatively small footprints of the observed fields, the data are an excellent proxy for the deep photometric datasets that will be acquired by LSST, and are therefore an ideal test bed for the analysis methods being implemented by the LSST DESC. We select a magnitude limited sample with $i<24.5$ and analyze it in four redshift bins covering $0.15\lesssim z \lesssim1.5$. We carry out a Fourier-space analysis of the two-point clustering of this sample, including all auto- and cross-correlations. We demonstrate the use of map-level deprojection methods to account for fluctuations in the galaxy number density caused by observational systematics. Through an HOD analysis, we place constraints on the characteristic halo masses of this sample, finding a good fit up to scales $k_{\rm max}=1\,{\rm Mpc}^{-1}$, including both auto- and cross-correlations. Our results show monotonically decreasing average halo masses, which can be interpreted in terms of the drop-out of red galaxies at high redshifts for a flux-limited sample. In terms of photometric redshift systematics, we show that additional care is needed in order to marginalize over uncertainties in the redshift distribution in galaxy clustering, and that these uncertainties can be constrained by including cross-correlations. We are able to make a $\sim3\sigma$ detection of lensing magnification in the HSC data. Our results are stable to variations in $\sigma_8$ and $\Omega_c$ and we find constraints that agree well with measurements from Planck and low-redshift probes. Finally, we use our pipeline to study the clustering of galaxies as a function of limiting flux, and provide a simple fitting function for the linear galaxy bias for magnitude limited samples as a function of limiting magnitude and redshift. [abridged]

We apply the Fourier Power Function Shapelets (FPFS) shear estimator to the first year data of the Hyper Suprime-Cam survey to construct a shape catalog. The FPFS shear estimator has been demonstrated to have multiplicative bias less than one percent in the absence of blending, regardless of complexities of galaxy shapes, smears of point spread functions (PSFs) and contamination from noise. The blending bias is calibrated with realistic image simulations, which include the impact of neighboring objects, using the COSMOS Hubble Space Telescope images. Here we carefully test the influence of PSF model residual on the FPFS shear estimation and the uncertainties in the shear calibration. Internal null tests are conducted to characterize potential systematics in the FPFS shape catalog and the results are compared with those measured using a catalog where the shapes were estimated using the re-Gaussianization algorithms. Furthermore, we compare various weak lensing measurements between the FPFS shape catalog and the re-Gaussianization shape catalog and conclude that the weak lensing measurements between these two shape catalogs are consistent with each other within the statistical uncertainty.

We present the first measurement of cross-correlation between the lensing potential, reconstructed from cosmic microwave background (CMB) {\it polarization} data, and the cosmic shear field from galaxy shapes. This measurement is made using data from the POLARBEAR CMB experiment and the Subaru Hyper Suprime-Cam (HSC) survey. By analyzing an 11~deg$^2$ overlapping region, we reject the null hypothesis at 3.5$\sigma$\ and constrain the amplitude of the {\bf cross-spectrum} to $\widehat{A}_{\rm lens}=1.70\pm 0.48$, where $\widehat{A}_{\rm lens}$ is the amplitude normalized with respect to the Planck~2018{} prediction, based on the flat $\Lambda$ cold dark matter cosmology. The first measurement of this {\bf cross-spectrum} without relying on CMB temperature measurements is possible due to the deep POLARBEAR map with a noise level of ${\sim}$6\,$\mu$K-arcmin, as well as the deep HSC data with a high galaxy number density of $n_g=23\,{\rm arcmin^{-2}}$. We present a detailed study of the systematics budget to show that residual systematics in our results are negligibly small, which demonstrates the future potential of this cross-correlation technique.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

We report the weak lensing measurement of filaments between Sloan Digital Sky Survey (SDSS) III/Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy pairs at $z\sim0.55$, using the Subaru Hyper Suprime-Cam (HSC) first-year galaxy shape catalogue. Despite of the small overlap of $140$ deg$^2$ between these surveys we detect the filament lensing signal at 3.9$\sigma$ significance, which is the highest signal-to-noise lensing measurement of filaments between galaxy-scale halos at this redshift range. We derive a theoretical prediction and covariance using mock catalogues based on full-sky ray-tracing simulations. We find that the intrinsic scatter of filament properties and the fluctuations in large scale structure along the line-of-sight are the primary component of the covariance and the intrinsic shape noise from source galaxies no longer limits our lensing measurement. This fact demonstrates the statistical power of the HSC survey due to its deep observations and high number density of source galaxies. Our result is consistent with the theoretical prediction and supports the "thick" filament model. As the HSC survey area increases, we will be able to study detailed filament properties such as the dark matter distributions and redshift evolution of filaments.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

We use the full-sky ray-tracing weak lensing simulations to generate 2268 mock catalogues for the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalogue. Our mock catalogues take into account various effects as in the real data: the survey footprints, inhomogeneous angular distribution of source galaxies, statistical uncertainties in photometric redshift (photo-$z$) estimate, variations in the lensing weight, and the statistical noise in galaxy shape measurements including both intrinsic shapes and the measurement errors. We then utilize our mock catalogues to evaluate statistical uncertainties expected in measurements of cosmic shear two-point correlations $\xi_{\pm}$ with tomographic redshift information for the HSC survey. We develop a quasi-analytical formula for the Gaussian sample variance properly taking into account the number of source pairs in the survey footprints. The standard Gaussian formula significantly overestimates or underestimates the mock results by $50\%$ level. We also show that different photo-$z$ catalogues or the six disconnected fields, rather than a consecutive geometry, cause variations in the covariance by $\sim 5\%$. The mock catalogues enable us to study the chi-square distribution for $\xi_{\pm}$. We find the wider distribution than that naively expected for the distribution with the degrees-of-freedom of data vector used. Finally, we propose a method to include non-zero multiplicative bias in mock shape catalogue and show the non-zero multiplicative bias can change the effective shape noise in cosmic shear analyses. Our results suggest an importance of estimating an accurate form of the likelihood function (and therefore the covariance) for robust cosmological parameter inference from the precise measurements.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

The Large Synoptic Survey Telescope (LSST) Dark Energy Science Collaboration (DESC) will use five cosmological probes: galaxy clusters, large scale structure, supernovae, strong lensing, and weak lensing. This Science Requirements Document (SRD) quantifies the expected dark energy constraining power of these probes individually and together, with conservative assumptions about analysis methodology and follow-up observational resources based on our current understanding and the expected evolution within the field in the coming years. We then define requirements on analysis pipelines that will enable us to achieve our goal of carrying out a dark energy analysis consistent with the Dark Energy Task Force definition of a Stage IV dark energy experiment. This is achieved through a forecasting process that incorporates the flowdown to detailed requirements on multiple sources of systematic uncertainty. Future versions of this document will include evolution in our software capabilities and analysis plans along with updates to the LSST survey strategy.

We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

We constrain the scaling relation between optical richness ($\lambda$) and halo mass ($M$) for a sample of SDSS redMaPPer galaxy clusters within the context of the {\it Planck} cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, $P(\ln \lambda| M)$. To model the abundance and the stacked lensing profiles, we use an emulator specifically built to interpolate the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, which is calibrated based on a suite of high-resolution $N$-body simulations. We apply our method to 8,312 SDSS redMaPPer clusters with $20\le \lambda \le 100$ and $0.10\le z_{\lambda}\le0.33$, and show that the log-normal distribution model for $P(\lambda|M)$, with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, $\left\langle \ln{\lambda}\right\rangle(M)=A+B\ln(M/M_{\rm pivot})$, with $A=3.207^{+0.044}_{-0.046}$ and $B=0.993^{+0.041}_{-0.055}$ (68\%~CL), where the pivot mass scale $M_{\rm pivot}=3\times 10^{14} h^{-1}M_\odot$, and the scatter $\sigma_{\mathrm{\ln\lambda}|M}=\sigma_0+q\ln(M/M_{\rm pivot})$ with $\sigma_0=0.456^{+0.047}_{-0.039}$ and $q=-0.169^{+0.035}_{-0.026}$. We find that a large scatter in halo masses is required at the lowest richness bins ($20\le \lambda \lesssim 30$) in order to reproduce the measurements. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of intrinsically low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of $30\le \lambda \le 100$.

Using $\sim$140 deg$^2$ Subaru Hyper Suprime-Cam (HSC) survey data, we stack the weak lensing (WL) signal around five Planck clusters found within the footprint. This yields a 15$\sigma$ detection of the mean Planck cluster mass density profile. The five Planck clusters span a relatively wide mass range, $M_{\rm WL,500c} = (2-30)\times10^{14}\,M_\odot/h$ with a mean mass of $M_{\rm WL,500c} = (4.15\pm0.61)\times10^{14}\,M_\odot/h$. The ratio of the stacked Planck Sunyaev-Zel'dovich (SZ) mass to the stacked WL mass is $ \langle M_{\rm SZ}\rangle/\langle M_{\rm WL}\rangle = 1-b = 0.80\pm0.14$. This mass bias is consistent with previous WL mass calibrations of Planck clusters within the errors. We discuss the implications of our findings for the calibration of SZ cluster counts and the much discussed tension between Planck SZ cluster counts and Planck $\Lambda$CDM cosmology.

We present optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.

We present a catalog of 182 galaxy clusters detected through the Sunyaev-Zel'dovich effect by the Atacama Cosmology Telescope in a contiguous 987.5 deg$^{2}$ field. The clusters were detected as SZ decrements by applying a matched filter to 148 GHz maps that combine the original ACT equatorial survey with data from the first two observing seasons using the ACTPol receiver. Optical/IR confirmation and redshift measurements come from a combination of large public surveys and our own follow-up observations. Where necessary, we measured photometric redshifts for clusters using a pipeline that achieves accuracy $\Delta z/(1 + z)=0.015$ when tested on SDSS data. Under the assumption that clusters can be described by the so-called Universal Pressure Profile and its associated mass-scaling law, the full signal-to-noise > 4 sample spans the mass range $1.6 < M^{\rm UPP}_{\rm 500c}/10^{14}{\rm M}_{\odot}<9.1$, with median $M^{\rm UPP}_{\rm 500c}=3.1 \times 10^{14}$ M$_{\odot}$. The sample covers the redshift range $0.1 < z < 1.4$ (median $z = 0.49$) and 28 clusters are new discoveries (median $z = 0.80$). We compare our catalog with other overlapping cluster samples selected using the SZ, optical,and X-ray wavelengths. We find the ratio of the UPP-based SZ mass to richness-based weak-lensing mass is $\langle M^{\rm UPP}_{\rm 500c} \rangle / \langle M^{\rm \lambda WL}_{\rm 500c} \rangle = 0.68 \pm 0.11$. After applying this calibration, the mass distribution for clusters with $M_{\rm 500c} > 4 \times 10^{14}$ M$_{\odot}$ is consistent with the number of such clusters found in the South Pole Telescope SZ survey.

Future high-resolution measurements of the cosmic microwave background (CMB) will produce catalogs of tens of thousands of galaxy clusters through the thermal Sunyaev-Zel'dovich (tSZ) effect. We forecast how well different configurations of a CMB Stage-4 experiment can constrain cosmological parameters, in particular the amplitude of structure as a function of redshift $\sigma_8(z)$, the sum of neutrino masses $\Sigma m_{\nu}$, and the dark energy equation of state $w(z)$. A key element of this effort is calibrating the tSZ scaling relation by measuring the lensing signal around clusters. We examine how the mass calibration from future optical surveys like the Large Synoptic Survey (LSST) compares with a purely internal calibration using lensing of the CMB itself. We find that, due to its high-redshift leverage, internal calibration gives constraints on cosmological parameters comparable to the optical calibration, and can be used as a cross-check of systematics in the optical measurement. We also show that in contrast to the constraints using the CMB lensing power spectrum, lensing-calibrated tSZ cluster counts can detect a minimal $\Sigma m_{\nu}$ at the 3-5$\sigma$ level even when the dark energy equation of state is freed up.

We present and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

We present wide-field (167 deg$^2$) weak lensing mass maps from the Hyper Supreme-Cam Subaru Strategic Program (HSC-SSP). We compare these weak lensing based dark matter maps with maps of the distribution of the stellar mass associated with luminous red galaxies. We find a strong correlation between these two maps with a correlation coefficient of $\rho=0.54\pm0.03$ (for a smoothing size of $8'$). This correlation is detected even with a smaller smoothing scale of $2'$ ($\rho=0.34\pm 0.01$). This detection is made uniquely possible because of the high source density of the HSC-SSP weak lensing survey ($\bar{n}\sim 25$ arcmin$^{-2}$). We also present a variety of tests to demonstrate that our maps are not significantly affected by systematic effects. By using the photometric redshift information associated with source galaxies, we reconstruct a three-dimensional mass map. This three-dimensional mass map is also found to correlate with the three-dimensional galaxy mass map. Cross-correlation tests presented in this paper demonstrate that the HSC-SSP weak lensing mass maps are ready for further science analyses.

We develop a method to simulate galaxy-galaxy weak lensing by utilizing all-sky, light-cone simulations and their inherent halo catalogs. Using the mock catalog to study the error covariance matrix of galaxy-galaxy weak lensing, we compare the full covariance with the "jackknife" (JK) covariance, the method often used in the literature that estimates the covariance from the resamples of the data itself. We show that there exists the variation of JK covariance over realizations of mock lensing measurements, while the average JK covariance over mocks can give a reasonably accurate estimation of the true covariance up to separations comparable with the size of JK subregion. The scatter in JK covariances is found to be $\sim10\%$ after we subtract the lensing measurement around random points. However, the JK method tends to underestimate the covariance at the larger separations, more increasingly for a survey with a higher number density of source galaxies. We apply our method to the the Sloan Digital Sky Survey (SDSS) data, and show that the 48 mock SDSS catalogs nicely reproduce the signals and the JK covariance measured from the real data. We then argue that the use of the accurate covariance, compared to the JK covariance, allows us to use the lensing signals at large scales beyond a size of the JK subregion, which contains cleaner cosmological information in the linear regime.

We present an optically-selected cluster catalog from the Hyper Suprime-Cam (HSC) Subaru Strategic Program. The HSC images are sufficiently deep to detect cluster member galaxies down to $M_*\sim 10^{10.2}M_\odot$ even at $z\sim 1$, allowing a reliable cluster detection at such high redshifts. We apply the CAMIRA algorithm to the HSC Wide S16A dataset covering $\sim 232$ deg$^2$ to construct a catalog of 1921 clusters at redshift $0.1<z<1.1$ and richness $\hat{N}_{\rm mem}>15$ that roughly corresponds to $M_{\rm 200m}\gtrsim 10^{14}h^{-1}M_\odot$. We confirm good cluster photometric redshift performance, with the bias and scatter in $\Delta z/(1+z)$ being better than 0.005 and 0.01 over most of the redshift range, respectively. We compare our cluster catalog with large X-ray cluster catalogs from XXL and XMM-LSS surveys and find good correlation between richness and X-ray properties. We also study the miscentering effect from the distribution of offsets between optical and X-ray cluster centers. We confirm the high ($>0.9$) completeness and purity for high mass clusters by analyzing mock galaxy catalogs.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

We measure the gravitational lensing shear signal around dark matter halos hosting CMASS galaxies using light sources at $z\sim 1$ (background galaxies) and at the surface of last scattering at $z\sim 1100$ (the cosmic microwave background). The galaxy shear measurement uses data from the CFHTLenS survey, and the microwave background shear measurement uses data from the {\it Planck} satellite. The ratio of shears from these cross-correlations provides a purely geometric distance measurement across the longest possible cosmological lever arm. This is because the matter distribution around the halos, including uncertainties in galaxy bias and systematic errors such as miscentering, cancels in the ratio for halos in thin redshift slices. We measure this distance ratio in three different redshift slices of the CMASS sample, and combine them to obtain a $17\%$ measurement of the distance ratio, $r=0.390^{+0.070}_{-0.062}$ at an effective redshift of $z=0.53$. This is consistent with the predicted ratio from the {\it Planck} best-fit $\Lambda$CDM cosmology of $r=0.419$.

We measure a combination of gravitational lensing, galaxy clustering, and redshift-space distortions called $E_G$. The quantity $E_G$ probes both parts of metric potential and is insensitive to galaxy bias and $\sigma_8$. These properties make it an attractive statistic to test $\Lambda$CDM, General Relativity and its alternate theories. We have combined CMASS DR11 with CFHTLenS and recent measurements of $\beta$ from RSD analysis, and find $E_G(z = 0.57) = 0.42 \pm 0.056$, an 13\% measurement in agreement with the prediction of general relativity $E_G(z = 0.57) = 0.396 \pm 0.011$ using the Planck 2015 cosmological parameters. We have corrected our measurement for various observational and theoretical systematics. Our measurement is consistent with the first measurement of $E_G$ using CMB lensing in place of galaxy lensing (Pullen et. al. 2015a) at small scales, but shows 2.8$\sigma$ tension when compared with their final results including large scales. This analysis with future surveys will provide improved statistical error and better control over systematics to test General Relativity and its alternate theories.

Mass calibration uncertainty is the largest systematic effect for using clusters of galaxies to constrain cosmological parameters. We present weak lensing mass measurements from the Canada-France-Hawaii Telescope Stripe 82 Survey for galaxy clusters selected through their high signal-to-noise thermal Sunyaev-Zeldovich (tSZ) signal measured with the Atacama Cosmology Telescope (ACT). For a sample of 9 ACT clusters with a tSZ signal-to-noise greater than five the average weak lensing mass is $\left(4.8\pm0.8\right)\,\times10^{14}\,\mathrm{M}_\odot$, consistent with the tSZ mass estimate of $\left(4.70\pm1.0\right)\,\times10^{14}\,\mathrm{M}_\odot$ which assumes a universal pressure profile for the cluster gas. Our results are consistent with previous weak-lensing measurements of tSZ-detected clusters from the Planck satellite. When comparing our results, we estimate the Eddington bias correction for the sample intersection of Planck and weak-lensing clusters which was previously excluded.

The next generation weak lensing surveys (i.e., LSST, Euclid and WFIRST) will require exquisite control over systematic effects. In this paper, we address shear calibration and present the most realistic forecast to date for LSST/Euclid/WFIRST and CMB lensing from a stage 4 CMB experiment (CMB S4). We use the CosmoLike code to simulate a joint analysis of all the two-point functions of galaxy density, galaxy shear and CMB lensing convergence. We include the full Gaussian and non-Gaussian covariances and explore the resulting joint likelihood with Monte Carlo Markov Chains. We constrain shear calibration biases while simultaneously varying cosmological parameters, galaxy biases and photometric redshift uncertainties. We find that CMB lensing from CMB S4 enables the calibration of the shear biases down to 0.2% - 3% in 10 tomographic bins for LSST (below the ~0.5% requirements in most tomographic bins), down to 0.4% - 2.4% in 10 bins for Euclid and 0.6% - 3.2% in 10 bins for WFIRST. For a given lensing survey, the method works best at high redshift where shear calibration is otherwise most challenging. This self-calibration is robust to Gaussian photometric redshift uncertainties and to a reasonable level of intrinsic alignment. It is also robust to changes in the beam and the effectiveness of the component separation of the CMB experiment, and slowly dependent on its depth, making it possible with third generation CMB experiments such as AdvACT and SPT-3G, as well as the Simons Observatory.

We search for the lensing signal of massive filaments between 135,000 pairs of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey. We develop a new estimator that cleanly removes the much larger shear signal of the neighboring LRG halos, relying only on the assumption of spherical symmetry. We consider two models: a thick filament model constructed from ray-tracing simulations for $\Lambda$CDM model, and a thin filament model which models the filament by a string of halos along the line connecting the two LRGs. We show that the filament lensing signal is in nice agreement with the thick simulation filament, while strongly disfavoring the thin model. The magnitude of the lensing shear due to the filament is below $10^{-4}$. Employing the likelihood ratio test, we find a 4.5$\sigma$ significance for the detection of the filament lensing signal, corresponding to a null hypothesis fluctuation probability of $3 \times 10^{-6}$. We also carried out several null tests to verify that the residual shear signal from neighboring LRGs and other shear systematics are minimized.

We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$\sigma$ in the weak lensing and 4.4$\sigma$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

We perform a joint analysis of the abundance, the clustering and the galaxy-galaxy lensing signal of galaxies measured from Data Release 11 of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (SDSS III-BOSS) in our companion paper, Miyatake et al. (2014). The lensing signal was obtained by using the shape catalog of background galaxies from the Canada France Hawaii Telescope Legacy Survey, which was made publicly available by the CFHTLenS collaboration, with an area overlap of about 105 deg$^2$. We analyse the data in the framework of the halo model in order to fit halo occupation parameters and cosmological parameters ($\Omega_{\rm m}$ and $\sigma_8$) to these observables simultaneously, and thus break the degeneracy between galaxy bias and cosmology. Adopting a flat $\Lambda$CDM cosmology with priors on $\Omega_b h^2$, $n_{\rm s}$ and $h$ from the analysis of WMAP 9-year data, we obtain constraints on the stellar mass-halo mass relation of galaxies in our sample. Marginalizing over the halo occupation distribution parameters and a number of other nuisance parameters in our model, we obtain $\Omega_{\rm m}=0.310^{+0.019}_{-0.020}$ and $\sigma_8=0.785^{+0.044}_{-0.044}$ (68% confidence). We demonstrate the robustness of our results with respect to sample selection and a variety of systematics such as the halo off-centering effect and possible incompleteness in our sample. Our constraints are consistent, complementary and competitive with those obtained using other independent probes of these cosmological parameters. The cosmological analysis is the first of its kind to be performed at a redshift as high as $0.53$.

A joint analysis of the clustering of galaxies and their weak gravitational lensing signal is well-suited to simultaneously constrain the galaxy-halo connection as well as the cosmological parameters by breaking the degeneracy between galaxy bias and the amplitude of clustering signal. In a series of two papers, we perform such an analysis at the highest redshift ($z\sim0.53$) in the literature using CMASS galaxies in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Eleventh Data Release (SDSS-III/BOSS DR11) catalog spanning 8300~deg$^2$. In this paper, we present details of the clustering and weak lensing measurements of these galaxies. We define a subsample of 400,916 CMASS galaxies based on their redshifts and stellar mass estimates so that the galaxies constitute an approximately volume-limited and similar population over the redshift range $0.47\le z\le 0.59$. We obtain a signal-to-noise ratio $S/N\simeq 56$ for the galaxy clustering measurement. We also explore the redshift and stellar mass dependence of the clustering signal. For the weak lensing measurement, we use existing deeper imaging data from the CFHTLS with publicly available shape and photometric redshift catalogs from CFHTLenS, but only in a 105~deg$^2$ area which overlaps with BOSS. This restricts the lensing measurement to only 5,084 CMASS galaxies. After careful systematic tests, we find a highly significant detection of the CMASS weak lensing signal, with total $S/N\simeq 26$. These measurements form the basis of the halo occupation distribution and cosmology analysis presented in More et al. (Paper II).

GALSIM is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GALSIM software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GALSIM meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA's Euclid mission, and NASA's WFIRST-AFTA mission. The GALSIM project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GALSIM repository can be found at https://github.com/GalSim-developers/GalSim .

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is an image analysis competition that aims to test algorithms to measure weak gravitational lensing from astronomical images. The challenge started in October 2013 and ends 30 April 2014. The challenge focuses on testing the impact on weak lensing measurements of realistically complex galaxy morphologies, realistic point spread function, and combination of multiple different exposures. It includes simulated ground- and space-based data. The details of the challenge are described in [15], and the challenge website and its leader board can be found at http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/, respectively.

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present a Subaru weak lensing measurement of ACT-CL J0022.2-0036, one of the most luminous, high-redshift (z=0.81) Sunyaev-Zel'dovich (SZ) clusters discovered in the 268 deg^2 equatorial region survey of the Atacama Cosmology Telescope. For the weak lensing analysis using i'-band images, we use a model-fitting (Gauss-Laguerre shapelet) method to measure shapes of galaxy images, where we fit galaxy images in different exposures simultaneously to obtain best-fit ellipticities taking into account the different PSFs in each exposure. We also take into account the astrometric distortion effect on galaxy images by performing the model fitting in the world coordinate system. To select background galaxies behind the cluster at z=0.81, we use photometric redshift (photo-z) estimates for every galaxy derived from the co-added images of multi-passband Br'i'z'Y, with PSF matching/homogenization. After a photo-z cut for background galaxy selection, we detect the tangential weak lensing distortion signal with a total signal-to-noise ratio of about 3.7. By fitting a Navarro-Frenk-White model to the measured shear profile, we find the cluster mass to be M_200\bar{\rho}_m = [7.5^+3.2_-2.8(stat.)^+1.3_-0.6(sys.)] x 10^14 M_\odot/h. The weak lensing-derived mass is consistent with previous mass estimates based on the SZ observation, with assumptions of hydrostatic equilibrium and virial theorem, as well as with scaling relations between SZ signal and mass derived from weak lensing, X-ray, and velocity dispersion, within the measurement errors. We also show that the existence of ACT-CL J0022.2-0036 at z=0.81 is consistent with the cluster abundance prediction of the \Lambda-dominated cold dark matter structure formation model. We thus demonstrate the capability of Subaru-type ground-based images for studying weak lensing of high-redshift clusters.