
Hyper SuprimeCam (HSC) is a widefield 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 300night multiband imaging survey of the highlatitude sky. The survey
includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad
bands ($grizy$), with a $5\,\sigma$ pointsource 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 report the discovery of a new gravitationally lensed quasar, SDSS
J131339.98+515128.3, at a redshift of 1.875 with an image separation of 1.24".
The lensing galaxy is clearly detected in visiblelight followup observations.
We also identify three absorptionline doublets in the spectra of the lensed
quasar images, from which we measure the lens redshift to be 0.194. Like
several other known lenses, the lensed quasar images have different continuum
slopes. This difference is probably the result of reddening and microlensing in
the lensing galaxy. The lensed quasar was selected by correlating Sloan Digital
Sky Survey (SDSS) spectroscopic quasars with Two Micron All Sky Survey (2MASS)
sources and choosing quasars that show nearinfrared (IR) excess. The nearIR
excess can originate, for example, from the contribution of the lensing galaxy
at nearIR wavelengths. We show that the nearIR excess technique is indeed an
efficient method to identify lensed systems from a large sample of quasars.

We present measurements of the normalised redshiftspace threepoint
correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS)
main galaxy sample. We have applied our "npt" algorithm to both a
volumelimited (36738 galaxies) and magnitudelimited sample (134741 galaxies)
of SDSS galaxies, and find consistent results between the two samples, thus
confirming the weak luminosity dependence of Q_z recently seen by other
authors. We compare our results to other Q_z measurements in the literature and
find it to be consistent within the full jackknife error estimates. However,
we find these errors are significantly increased by the presence of the ``Sloan
Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the
3point correlation function (3PCF) by 70% on large scales (s>=10h^1 Mpc). If
we exclude this supercluster, our observed Q_z is in better agreement with that
obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity
of these higherorder correlation functions to largescale structures in the
Universe. This analysis highlights that the SDSS datasets used here are not
``fair samples'' of the Universe for the estimation of higherorder clustering
statistics and larger volumes are required. We study the shapedependence of
Q_z(s,q,theta) as one expects this measurement to depend on scale if the large
scale structure in the Universe has grown via gravitational instability from
Gaussian initial conditions. On small scales (s <= 6h^1 Mpc), we see some
evidence for shapedependence in Q_z, but at present our measurements are
consistent with a constant within the errors (Q_z ~ 0.75 +/ 0.05). On scales
>10h^1 Mpc, we see considerable shapedependence in Q_z.

We study the luminosity and color dependence of the galaxy 2point
correlation function in the Sloan Digital Sky Survey, starting from a sample of
200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation
function w(r_p), which is directly related to the real space \xi(r). The
amplitude of w(r_p) grows continuously with luminosity, rising more steeply
above the characteristic luminosity L_*. Redder galaxies exhibit a higher
amplitude and steeper correlation function at all luminosities. The correlation
amplitude of blue galaxies increases continuously with luminosity, but the
luminosity dependence for red galaxies is less regular, with bright red
galaxies more strongly clustered at large scales and faint red galaxies more
strongly clustered at small scales. We interpret these results using halo
occupation distribution (HOD) models assuming concordance cosmological
parameters. For most samples, an HOD model with two adjustable parameters fits
the w(r_p) data better than a powerlaw, explaining inflections at r_p ~ 13
Mpc/h as the transition between the 1halo and 2halo regimes of \xi(r). The
implied minimum mass for a halo hosting a central galaxy above a luminosity
threshold L grows as M_min ~ L at low luminosities and more steeply above L_*.
The mass at which an average halo has one satellite galaxy brighter than L is
M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at
fixed halo mass) in which central galaxies lie far above a Schechter function
extrapolation of the satellite population. HOD models nicely explain the joint
luminositycolor dependence of w(r_p) in terms of the color fractions of
central and satellite populations as a function of halo mass. The inferred HOD
properties are in good qualitative agreement with theoretical predictions.

The statistical property of the weak lensing fields is studied quantitatively
using the raytracing simulations. Motivated by the empirical lognormal model
that characterizes the probability distribution function(PDF) of the
threedimensional mass distribution excellently, we critically investigate the
validity of lognormal model in the weak lensing statistics. Assuming that the
convergence field, $\kappa$, is approximately described by the lognormal
distribution, we present analytic formulae of convergence for the onepoint PDF
and the Minkowski functionals. Comparing those predictions with raytracing
simulations in various cold dark matter models, we find that the onepoint
lognormal PDF can describe the nonGaussian tails of convergence fields
accurately up to $\nu\sim10$, where $\nu$ is the level threshold given by
$\nu\equiv\kappa/\var^{1/2}$, although the systematic deviation from lognormal
prediction becomes manifest at higher source redshift and larger smoothing
scales. The lognormal formulae for Minkowski functionals also fit to the
simulation results when the source redshift is low. Accuracy of the
lognormalfit remains good even at the small angular scales, where the
perturbation formulae by Edgeworth expansion break down. On the other hand,
lognormal models does not provide an accurate prediction for the statistics
sensitive to the rare events such as the skewness and the kurtosis of
convergence. We therefore conclude that the empirical lognormal model of the
convergence field is safely applicable as a useful cosmological tool, as long
as we are concerned with the nonGaussianity of $\nu\simlt5$ for low source
redshift samples.