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

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.

Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constraint on the cosmological parameters. For this purpose, the accurate shape measurement of background galaxies is absolutely important in which any systematic error in the measurement should be carefully corrected. One of the main systematic error comes from photon noise which is Poisson noise of flux from the atmosphere(noise bias). We investigate how the photon noise makes a systematic error in shear measurement within the framework of ERA method we developed in earlier papers and gives a practical correction method. The method is tested by simulations with real galaxy images with various conditions and it is confirmed that it can correct $80 \sim 90\%$ of the noise bias except for galaxies with very low signal to noise ratio.

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.

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

Highly precise weak lensing shear measurement is required for statistical weak gravitational lensing analysis such as cosmic shear measurement to achieve severe constrain on the cosmological parameters. For this purpose any systematic error in the measurement should be corrected. One of the main systematic error comes from Pixel noise which is Poisson noise of flux from atmosphere. We investigate how the pixel noise makes systematic error in shear measurement based on ERA method and develop the correction method. This method is tested by simulations with various conditions and it is confirmed that the correction method can correct $80 \sim 90\%$ of the systematic error except very low signal to noise ratio galaxies.

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 improve the ERA(Ellipticity of Re-smeared Artificial image) method of PSF(Point Spread Function) correction in weak lensing shear analysis in order to treat realistic shape of galaxies and PSF. This is done by re-smearing PSF and the observed galaxy image smeared by a RSF(Re-Smearing Function), and allows us to use a new PSF with a simple shape and to correct PSF effect without any approximations and assumptions. We perform numerical test to show that the method applied for galaxies and PSF with some complicated shapes can correct PSF effect with systematic error less than 0.1%. We also apply ERA method for real data of Abell 1689 cluster to confirm that it is able to detect the systematic weak lensing shear pattern. The ERA method requires less than 0.1 or 1 second to correct PSF for each object in numerical test and real data analysis, respectively.

Weak gravitational lensing causes subtle changes in the apparent shapes of galaxies due to the bending of light by the gravity of foreground masses. By measuring the shapes of large numbers of galaxies (millions in recent surveys, up to tens of billions in future surveys) we can infer the parameters that determine cosmology. Imperfections in the detectors used to record images of the sky can introduce changes in the apparent shape of galaxies, which in turn can bias the inferred cosmological parameters. In this paper we consider the effect of two widely discussed sensor imperfections: tree-rings, due to impurity gradients which cause transverse electric fields in the Charge-Coupled Devices (CCD), and pixel-size variation, due to periodic CCD fabrication errors. These imperfections can be observed when the detectors are subject to uniform illumination (flat field images). We develop methods to determine the spurious shear and convergence (due to the imperfections) from the flat-field images. We calculate how the spurious shear when added to the lensing shear will bias the determination of cosmological parameters. We apply our methods to candidate sensors of the Large Synoptic Survey Telescope (LSST) as a timely and important example, analyzing flat field images recorded with LSST prototype CCDs in the laboratory. We find that tree-rings and periodic pixel-size variation present in the LSST CCDs will introduce negligible bias to cosmological parameters determined from the lensing power spectrum, specifically w, {\Omega}m and {\sigma}8.

We study the extent to which very bright (-23.0 < MUV < -21.75) Lyman-break selected galaxies at redshifts z~7 display detectable Lya emission. To explore this issue, we have obtained follow-up optical spectroscopy of 9 z~7 galaxies from a parent sample of 24 z~7 galaxy candidates selected from the 1.65 sq.deg COSMOS-UltraVISTA and SXDS-UDS survey fields using the latest near-infrared public survey data, and new ultra-deep Subaru z'-band imaging (which we also present and describe in this paper). Our spectroscopy has yielded only one possible detection of Lya at z=7.168 with a rest-frame equivalent width EW_0 = 3.7 (+1.7/-1.1) Angstrom. The relative weakness of this line, combined with our failure to detect Lya emission from the other spectroscopic targets allows us to place a new upper limit on the prevalence of strong Lya emission at these redshifts. For conservative calculation and to facilitate comparison with previous studies at lower redshifts, we derive a 1-sigma upper limit on the fraction of UV bright galaxies at z~7 that display EW_0 > 50 Angstrom, which we estimate to be < 0.23. This result may indicate a weak trend where the fraction of strong Lya emitters ceases to rise, and possibly falls between z~6 and z~7. Our results also leave open the possibility that strong Lya may still be more prevalent in the brightest galaxies in the reionization era than their fainter counterparts. A larger spectroscopic sample of galaxies is required to derive a more reliable constraint on the neutral hydrogen fraction at z~7 based on the Lya fraction in the bright galaxies.

We generalize ERA method of PSF correction for more realistic situations. The method re-smears the observed galaxy image(galaxy image smeared by PSF) and PSF image by an appropriate function called Re-Smearing Function(RSF) to make new images which have the same ellipticity with the lensed (before smeared by PSF) galaxy image. It has been shown that the method avoids a systematic error arising from an approximation in the usual PSF correction in moment method such as KSB for simple PSF shape. By adopting an idealized PSF we generalize ERA method applicable for arbitrary PSF. This is confirmed with simulated complex PSF shapes. We also consider the effect of pixel noise and found that the effect causes systematic overestimation.

We present an analysis of the impact of the tree rings seen in the candidate sensors of the Large Synoptic Survey Telescope (LSST) on galaxy-shape measurements. The tree rings are a consequence of transverse electric fields caused by circularly symmetric impurity gradients in the silicon of the sensors. They effectively modify the pixel area and shift the photogenerated charge around, displacing the observed photon positions. The displacement distribution generates distortions that cause spurious shears correlated with the tree-rings patterns, potentially biasing cosmic shear measurements. In this paper we quantify the amplitude of the spurious shear caused by the tree rings on the LSST candidate sensors, and calculate its 2-point correlation function. We find that 2-point correlation function of the spurious shear on an area equivalent to the LSST field of view is order of about $10^{-13}$, providing a negligible contribution to the 2-point correlation of the cosmic shear signal. Additional work is underway, and the final results and analysis will be published elsewhere (Okura et al. (2015), in prep.)

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 developed a new method that uses ellipticity defined by 0th order moments (0th-ellipticity) for weak gravitational lensing shear analysis. Although there is a strong correlation between the ellipticity calculated using this approach and the usual ellipticity defined by the 2nd order moment, the ellipticity calculated here has a higher signal-to-noise ratio because it is weighted to the central region of the image. These results were confirmed using data for Abell 1689 from the Subaru telescope. For shear analysis, we adopted the ellipticity of re-smeared artificial image (ERA) method for point spread function (PSF) correction, and we tested the precision of this 0th-ellipticity with simple simulation, then we obtained the same level of precision with the results of ellipticity defined by quadrupole moments. Thus, we can expect that weak lensing analysis using 0 shear will be improved in proportion to the statistical error.

We propose a new method for Point Spread Function (PSF) correction in weak gravitational lensing shear analysis using an artificial image with the same ellipticity as the lensed image. This avoids the systematic error associated with the approximation in PSF correction used in previous approaches. We test the new method with simulated objects which have Gaussian or Cersic profiles smeared by a Gaussian PSF, and confirm that there is no systematic error.

We examined the anisotropic point spread function (PSF) of Suprime-Cam data utilizing dense star field data. We decomposed the PSF ellipticities into three components, the optical aberration, atmospheric turbulence, and chip-misalignment in an empirical manner, and evaluated the amplitude of each component. We found that, for long-exposure data, the optical aberration has the largest contribution to the PSF ellipticities, which could be modeled well by a simple analytic function based on the lowest-order aberration theory. The statistical properties of PSF ellipticities resulting from the atmospheric turbulence were investigated by using the numerical simulations. The simulation results are in a reasonable agreement with the observed data. It is also found that the optical PSF can be well corrected by the standard correction method with a polynomial fitting function. However, for the atmospheric PSF, its correction is affected by the common limitation caused by sparse sampling of PSFs due to a limited number of stars. We also examined the effects of the residual PSF anisotropies on Suprime-Cam cosmic shear data. We found that the shape and amplitude of the B-mode shear variance are broadly consistent with those of the residual PSF ellipticities measured from the dense star field data. This indicates that most of the sources of residual systematic are understood, which is an important step for cosmic shear statistics to be a practical tool of the precision cosmology.

This is the third paper on the improvements of systematic errors in our weak lensing analysis using an elliptical weight function, called E-HOLICs. In the previous papers we have succeeded in avoiding error which depends on ellipticity of background image. In this paper, we investigate the systematic error which depends on signal to noise ratio of background image. We find that the origin of the error is the random count noise which comes from Poisson noise of sky counts. Random count noise makes additional moments and centroid shift error, and those 1st orders are canceled in averaging, but 2nd orders are not canceled. We derived the equations which corrects these effects in measuring moments and ellipticity of the image and test their validity using simulation image. We find that the systematic error becomes less than 1% in the measured ellipticity for objects with $S/N>3$.

We have developed a new method(E-HOLICs) of estimating gravitational shear by adopting an elliptical weight function to measure background galaxy images in our previous paper. Following to the previous paper where isotropic Point Spread Function(PSF) correction is calculated, we consider an anisotropic PSF correction in this paper in order to apply E-HOLICs for real data. A an example, E-HOLICs is applied to Subaru data of massive and compact galaxy clusters A370, and is able to detect double peaks in the central region of the cluster consistent with the analysis of strong lensing. We also study the systematic error in E-HOLICs using STEP II simulation. In particular we consider the dependences of signal to noise ratio "S/N" of background galaxies in the shear estimation. Although E-HOLICs does improve systematic error due to the ellipticity dependence as shown in part 1, a systematic error due to the S/N dependence remains, namely E-HOLICs underestimates shear when background galaxies with low S/N objects are used. We discuss a possible improvement of the S/N dependence.

We develop a new method to estimate gravitational shear by adopting an elliptical weight function to measure background galaxy images. In doing so, we introduce a new concept of "zero plane" which is an imaginal source plane where shapes of all sources are perfect circles, and regard the intrinsic shear as the result of an imaginal lensing distortion. This makes the relation between the observed shear, the intrinsic shear and lensing distortion more simple and thus higher-order calculation more easy. The elliptical weight function allows us to measure the mutiplemoment of shape of background galaxies more precisely by weighting highly to brighter parts of image and moreover to reduce systematic error due to insufficient expansion of the weight function in the original approach of KSB. Point Spread Function(PSF) correction in E-HOLICs methods becomes more complicated than those in KSB methods. In this paper we studied isotropic PSF correction in detail. By adopting the lensing distortion as the ellipticity of the weight function, we are able to show that the shear estimation in E-HOLICs method reduces to solve a polynomial in the absolute magnitude of the distortion. We compare the systematic errors between our approach and KSB using STEP2 simulation. It is confirmed that KSB method overestimate the input shear for images with large ellipticities, and E-HOLICs correctly estimate the input shear even for such images. Anisotropic PSF correction and analysis of real data will be presented in forthcoming paper.

We obtain the projected mass distributions for two Subaru/Suprime-Cam fields in the southwest region (r\simlt 60') of the Coma cluster (z=0.0236) by weak lensing analysis and detect eight subclump candidates. We quantify the contribution of background large-scale structure (LSS) on the projected mass distributions using SDSS multi-bands and photometric data, under the assumption of mass-to-light ratio for field galaxies. We find that one of eight subclump candidates, which is not associated with any member galaxies, is significantly affected by LSS lensing. The mean projected mass for seven subclumps extracted from the main cluster potential is <M_2D^(corr)> = (5.06\pm1.30)10^12h^-1 M_sun after a LSS correction. A tangential distortion profile over an ensemble of subclumps is well described by a truncated singular-isothermal sphere model and a truncated NFW model. A typical truncated radius of subclumps, r_t\simeq 35 h^-1 kpc, is derived without assuming any relations between mass and light for member galaxies. The radius coincides well with the tidal radius, \sim42 h^-1 kpc, of the gravitational force of the main cluster. Taking into account the incompleteness of data area, a projection effect and spurious lensing peaks, it is expected that mass of cluster substructures account for 19 percent of the virial mass, with 13 percent statistical error. The mass fraction of cluster substructures is in rough agreement with numerical simulations.

We investigate the possibility to use higher order moments of gravitational lensed images in the weak lensing analysis. For this purpose we employ spin-2 component of HOLICs(Higher Order Lensing Characteristics) developed by us. We test the weak lensing analysis with spin-2 HOLICs using actual, ground based Subaru observations of the massive galaxy cluster A1689(z=0.183). It turns out that spin-2 HOLICs of order up to 8 are sufficiently applicable for weak lensing analysis after correcting PSF anisotropy as well as isotropic PSF smearing.

We have developed a method for measuring higher-order weak lensing distortions of faint background galaxies, namely the weak gravitational flexion, by fully extending the Kaiser, Squires & Broadhurst method to include higher-order lensing image characteristics (HOLICs) introduced by Okura, Umetsu, & Futamase. We take into account explicitly the weight function in calculations of noisy shape moments and the effect of higher-order PSF anisotropy, as well as isotropic PSF smearing. Our HOLICs formalism allows accurate measurements of flexion from practical observational data in the presence of non-circular, anisotropic PSF. We test our method using mock observations of simulated galaxy images and actual, ground-based Subaru observations of the massive galaxy cluster A1689 ($z=0.183$). From the high-precision measurements of spin-1 first flexion, we obtain a high-resolution mass map in the central region of A1689. The reconstructed mass map shows a bimodal feature in the central $4'\times 4'$ region of the cluster. The major, pronounced peak is associated with the brightest cluster galaxy and central cluster members, while the secondary mass peak is associated with a local concentration of bright galaxies. The refined, high-resolution mass map of A1689 demonstrates the power of the generalized weak lensing analysis techniques for quantitative and accurate measurements of the weak gravitational lensing signal.

We study a possibility to use the octopole moment of gravitationally lensed images as a direct measure of the third-order weak gravitational lensing effect, or the gravitational flexion. It turns out that there is a natural relation between flexion and certain combinations of octopole/higher-multipole moments which we call the Higher Order Lensing Image's Characteristics (HOLICs). This will allow one to measure directly flexion from observable octopole and higher-multipole moments of background images. We show based on simulated observations how the use of HOLICs can improve the accuracy and resolution of a reconstructed mass map, in which we assume Gaussian uncertainties in the shape measurements estimated using deep i'-band data of blank fields observed with Suprime-Cam on the Subaru telescope.