We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

The Beijing-Arizona Sky Survey (BASS) is a wide and deep imaging survey to cover a 5400 deg$^2$ area in the Northern Galactic Cap with the 2.3m Bok telescope using two filters ($g$ and $r$ bands). The Mosaic $z$-band Legacy Survey (MzLS) covers the same area in $z$ band with the 4m Mayall telescope. These two surveys will be used for spectroscopic targeting of the Dark Energy Spectroscopic Instrument (DESI). The BASS survey observations were completed in 2019 March. This paper describes the third data release (DR3) of BASS, which contains the photometric data from all BASS and MzLS observations between 2015 January and 2019 March. The median astrometric precision relative to {\it Gaia} positions is about 17 mas and the median photometric offset relative to the PanSTARRS1 photometry is within 5 mmag. The median $5\sigma$ AB magnitude depths for point sources are 24.2, 23.6, and 23.0 mag for $g$, $r$, and $z$ bands, respectively. The photometric depth within the survey area is highly homogeneous, with the difference between the 20\% and 80\% depth less than 0.3 mag. The DR3 data, including raw data, calibrated single-epoch images, single-epoch photometric catalogs, stacked images, and co-added photometric catalogs, are publicly accessible at \url{http://batc.bao.ac.cn/BASS/doku.php?id=datarelease:home}.

We investigate the distribution of companion galaxies around quasars using {\em Hubble Space Telescope} ({\em HST}) Advanced Camera for Surveys Wide Field Camera (ACS/WFC) archival images. Our master sample contains 532 quasars which have been observed by {\em HST} ACS/WFC, spanning a wide range of luminosity $(-31<M_i(z=2)<-23)$ and redshift ($0.3<z<3$). We search for companions around the quasars with projected distance of $10\text{ kpc}<d<100\text{ kpc}$. PSF subtraction is performed to enhance the completeness for close companions. The completeness is estimated to be high $(>90\%)$ even for the faintest companions of interest. The number of physical companions is estimated by subtracting a background density from the number density of projected companions. We divide all the companions into three groups (faint, intermediate and bright) according to their fluxes. A control sample of galaxies is constructed to have similar redshift distribution and stellar mass range as the quasar sample using the data from {\em HST} deep fields. We find that quasars and control sample galaxies have similar numbers of faint and bright companions, while quasars show a $3.7\sigma$ deficit of intermediate companions compared to galaxies. The numbers of companions in all three groups do not show strong evolution with redshift, and the number of intermediate companions around quasars decreases with quasar luminosity. Assuming that merger-triggered quasars have entered the final coalescence stage during which individual companions are no longer detectable at large separations, our result is consistent with a picture in which a significant fraction of quasars is triggered by mergers.

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17.

Among more than 200 quasars known at $z\gtrsim6$, only one object, J0100+2802 (z=6.327), was found hosting a $>10^{10}M_\odot$ super-massive black hole (SMBH). In order to investigate the host galaxy properties of J0100+2802, we performed multi-band ALMA observations, aiming at mapping the dust continuum, [CII] and CO(6-5) emission lines with sub-kiloparsec scale resolution, as well as detecting high-J CO lines in CO(11-10), CO(10-9), and CO(7-6). The galaxy size is measured to be $R_{\rm major}=3.6\pm0.2$ kpc from the high resolution continuum observations. No ordered motion on kilo-parsec scales was found in both [CII] and CO(6-5) emissions. The velocity dispersion is measured to be 161$\pm$7 km s$^{-1}$, $\sim$3 times smaller than that estimated from the local M-$\sigma$ relation. In addition, we found that the CO emission is more concentrate (a factor of 1.8$\pm$0.4) than the [CII] emission. Together with CO(2-1) detected by VLA, we measured the CO Spectral Line Energy Distribution (SLED), which is best fit by a two-components model, including a cool component at $\sim24$ K with a density of $n_{\rm (H_2)}$=10$^{4.5}$ cm$^{-3}$, and a warm component at $\sim224$ K with a density of $n_{\rm (H_2)}$=10$^{3.6}$ cm$^{-3}$. We also fitted the dust continuum with a graybody model, which indicates that it has either a high dust emissivity $\beta\gtrsim2$ or a hot dust temperature $T_{\rm dust}\gtrsim60$ K, or a combination of both factors. The highly excited CO emission and hot dust temperature suggest that the powerful AGN in J0100+2802 could contribute to the gas and dust heating although future observations are needed to confirm this.

We present the results of the Extremely Luminous Quasar Survey in the $3\pi$ survey of the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS; PS1). This effort applies the successful quasar selection strategy of the Extremely Luminous Survey in the Sloan Digital Sky Survey footprint ($\sim12,000\,\rm{deg}^2$) to a much larger area ($\sim\rm{21486}\,\rm{deg}^2$). This spectroscopic survey targets the most luminous quasars ($M_{1450}\le-26.5$; $m_{i}\le18.5$) at intermediate redshifts ($z\ge2.8$). Candidates are selected based on a near-infrared JKW2 color cut using WISE AllWISE and 2MASS photometry to mainly reject stellar contaminants. Photometric redshifts ($z_{\rm{reg}}$) and star-quasar classifications for each candidate are calculated from near-infrared and optical photometry using the supervised machine learning technique random forests. We select 806 quasar candidates at $z_{\rm{reg}}\ge2.8$ from a parent sample of 74318 sources. After exclusion of known sources and rejection of candidates with unreliable photometry, we have taken optical identification spectra for 290 of our 334 good PS-ELQS candidates. We report the discovery of 190 new $z\ge2.8$ quasars and an additional 28 quasars at lower redshifts. A total of 44 good PS-ELQS candidates remain unobserved. Including all known quasars at $z\ge2.8$, our quasar selection method has a selection efficiency of at least $77\%$. At lower declinations $-30\le\rm{Decl.}\le0$ we approximately triple the known population of extremely luminous quasars. We provide the PS-ELQS quasar catalog with a total of 592 luminous quasars ($m_{i}\le18.5$, $z\ge2.8$). This unique sample will not only be able to provide constraints on the volume density and quasar clustering of extremely luminous quasars, but also offers valuable targets for studies of the intergalactic medium.

We have designed the Extremely Luminous Quasar Survey (ELQS) to provide a highly complete census of unobscured UV-bright quasars during the cosmic noon, $z=2.8-5.0$. Here we report the discovery of 70 new quasars in the ELQS South Galactic Cap (ELQS-S) quasar sample, doubling the number of known extremely luminous quasars in $4,237.3\,\rm{deg}^2$ of the SDSS footprint. These observations conclude the ELQS and we present the properties of the full ELQS quasar catalog, containing 407 quasars over $11,838.5\,\rm{deg}^2$. Our novel ELQS quasar selection strategy resulted in unprecedented completeness at the bright end and allowed us to discover 109 new quasars in total. This marks an increase of $\sim36\%$ (109/298) to the known population at these redshifts and magnitudes, while we further are able to retain a selection efficiency of $\sim80\%$. On the basis of 166 quasars from the full ELQS quasar catalog, who adhere to the uniform criteria of the 2MASS point source catalog, we measure the bright-end quasar luminosity function (QLF) and extend it one magnitude brighter than previous studies. Assuming a single power law with exponential density evolution for the functional form of the QLF, we retrieve the best fit parameters from a maximum likelihood analysis. We find a steep bright-end slope of $\beta\approx-4.1$ and we can constrain the bright-end slope to $\beta\leq-3.4$ with $99\%$ confidence. The density is well modeled by the exponential redshift evolution, resulting in a moderate decrease with redshift ($\gamma\approx-0.4$).

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

This is the fourth paper in a series of publications aiming at discovering quasars at the epoch of reionization. In this paper, we expand our search for $z\sim 7$ quasars to the footprint of the Dark Energy Survey (DES) Data Release One (DR1), covering $\sim 5000$ deg$^2$ of new area. We select $z\sim 7$ quasar candidates using deep optical, near-infrared (near-IR) and mid-IR photometric data from the DES DR1, the VISTA Hemisphere Survey (VHS), the VISTA Kilo-degree Infrared Galaxy (VIKING) survey, the UKIRT InfraRed Deep Sky Surveys -- Large Area Survey (ULAS) and the unblurred coadds from the Wide-field Infrared Survey Explore ($WISE$) images (unWISE). The inclusion of DES and unWISE photometry allows the search to reach $\sim$ 1 magnitude fainter, comparing to our $z \gtrsim 6.5$ quasar survey in the northern sky (Wang et al. 2018). We report the initial discovery and spectroscopic confirmation of six new luminous quasars at $z>6.4$, including an object at $z=7.02$, the fourth quasar yet known at $z>7$, from a small fraction of candidates observed thus far. Based on the recent measurement of $z \sim 6.7 $ quasar luminosity function using the quasar sample from our survey in the northern sky, we estimate that there will be $\gtrsim$ 55 quasars at $z > 6.5$ at $M_{1450} < -24.5$ in the full DES footprint.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

We present the final results from our survey of luminous $z \sim $ 5.5 quasars. This is the first systematic quasar survey focusing on quasars at $z \sim$ 5.5, during the post-reionization epoch. It has been challenging to select quasars at $5.3 < z < 5.7$ using conventional color selections, due to their similar optical colors to those of late-type stars, especially M dwarfs. We developed a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR, and mid-IR photometry, using data from the Sloan Digital Sky Survey (SDSS), PanSTARR1 (PS1), the UKIRT Infrared Deep Sky Surveys - Large Area Survey, the UKIRT Hemisphere Survey, the VISTA Hemisphere Survey, and the Wide Field Infrared Survey Explorer ($WISE$), covering $\sim$ 11000 deg$^2$ of high galactic latitude sky. In this paper, we present the discovery of 15 new quasars at $z\sim 5.5$. Together with results from Yang et al. (2017), our survey provides a complete, flux-limited sample of 31 quasars at $5.3 \le z \le 5.7$. We measure the quasar spatial density at $z \sim 5.5$ and $M_{1450} < -26.2$. Our result is consistent with the rapid decline of quasar spatial density from $z = 5$ to 6, with $k=-0.66 \pm 0.05$ ($\rho(z) \propto 10^{kz}$). In addition, we present a new survey using optical colors only from the full PS1 area for luminous quasars at $z = 5.0 - 5.5$, beyond the SDSS footprint, and report the preliminary results from this survey, including 51 new quasars discovered at $4.61\le z \le5.71$.

The Beijing-Arizona Sky Survey (BASS) is a new wide-field legacy imaging survey in the northern Galactic cap using the 2.3m Bok telescope. The survey will cover about 5400 deg$^2$ in the $g$ and $r$ bands, and the expected 5$\sigma$ depths (corrected for the Galactic extinction) in the two bands are 24.0 and 23.4 mag, respectively. BASS started observations in January 2015, and has completed about 41% of the whole area as of July 2016. The first data release contains both calibrated images and photometric catalogs obtained in 2015 and 2016. The depths of single-epoch images in the two bands are 23.4 and 22.9 mag, and the full depths of three epochs are about 24.1 and 23.5 mag, respectively.

We present the rest-frame UV and optical photometry and morphology of low-redshift broad-line quasar host galaxies from the Sloan Digital Sky Survey Reverberation Mapping project. Our sample consists of 103 quasars at $z<0.8$, spanning a luminosity range of $-25\le M_g\le -17$ mag. We stack the multi-epoch images in the $g$ and $i$ bands taken by the Canada-France-Hawaii Telescope. The combined $g$-band ($i$-band) images reach a $5\sigma$ depth of 26.2 (25.2) mag, with a typical PSF size of $0.7"$ ($0.6"$). Each quasar is decomposed into a PSF and a S\'ersic profile, representing the central AGN and the host galaxy components, respectively. The systematic errors of the measured host galaxy flux in the two bands are 0.23 and 0.18 mag. The relative errors of the measured galaxy half-light radii ($R_e$) are about 13%. We estimate the rest-frame $u$- and $g$-band flux of the host galaxies, and find that the AGN-to-galaxy flux ratios in the $g$ band are between 0.9 to 4.4 (68.3% confidence). These galaxies have high stellar masses $M_\ast = 10^{10}\sim10^{11}\, M_\odot$. They have similar color with star-forming galaxies at similar redshifts, in consistent with AGN positive feedback in these quasars. We find that the $M_*-M_\text{BH}$ relation in our sample is shallower than the local $M_\text{Bulge}-M_\text{BH}$ relation. The S\'ersic indices and the $M_*-R_e$ relation indicate that the majority of the host galaxies are disk-like.

We present initial results from the first systematic survey of luminous $z\sim 5.5$ quasars. Quasars at $z \sim$ 5.5, the post-reionization epoch, are crucial tools to explore the evolution of intergalactic medium, quasar evolution and the early super-massive black hole growth. However, it has been very challenging to select quasars at redshifts 5.3 $\le z \le$ 5.7 using conventional color selections, due to their similar optical colors to late-type stars, especially M dwarfs, resulting in a glaring redshift gap in quasar redshift distributions. We develop a new selection technique for $z \sim$ 5.5 quasars based on optical, near-IR and mid-IR photometric data from Sloan Digital Sky Survey (SDSS), UKIRT InfraRed Deep Sky Surveys - Large Area Survey (ULAS), VISTA Hemisphere Survey (VHS) and Wide field Infrared Survey Explorer (WISE). From our pilot observations in SDSS-ULAS/VHS area, we have discovered 15 new quasars at 5.3 $\le z \le$ 5.7 and 6 new lower redshift quasars, with SDSS z band magnitude brighter than 20.5. Including other two $z \sim$ 5.5 quasars already published in our previous work, we now construct an uniform quasar sample at 5.3 $\le z \le$ 5.7 with 17 quasars in a $\sim$ 4800 square degree survey area. For further application in a larger survey area, we apply our selection pipeline to do a test selection by using the new wide field J band photometric data from a preliminary version of the UKIRT Hemisphere Survey (UHS). We successfully discover the first UHS selected $z \sim$ 5.5 quasar.

This is the second paper in a series on a new luminous z ~ 5 quasar survey using optical and near-infrared colors. Here we present a new determination of the bright end of the quasar luminosity function (QLF) at z ~ 5. Combined our 45 new quasars with previously known quasars that satisfy our selections, we construct the largest uniform luminous z ~ 5 quasar sample to date, with 99 quasars in the range 4.7 <= z < 5.4 and -29 < M1450 <= -26.8, within the Sloan Digital Sky Survey (SDSS) footprint. We use a modified 1/Va method including flux limit correction to derive a binned QLF, and we model the parametric QLF using maximum likelihood estimation. With the faint-end slope of the QLF fixed as alpha = -2.03 from previous deeper samples, the best fit of our QLF gives a flatter bright end slope beta = -3.58+/-0.24 and a fainter break magnitude M*1450 = -26.98+/-0.23 than previous studies at similar redshift. Combined with previous work at lower and higher redshifts, our result is consistent with a luminosity evolution and density evolution (LEDE) model. Using the best fit QLF, the contribution of quasars to the ionizing background at z ~ 5 is found to be 18% - 45% with a clumping factor C of 2 - 5. Our sample suggests an evolution of radio loud fraction with optical luminosity but no obvious evolution with redshift.

High-redshift quasars are important tracers of structure and evolution in the early universe. However, they are very rare and difficult to find when using color selection because of contamination from late-type dwarfs. High-redshift quasar surveys based on only optical colors suffer from incompleteness and low identification efficiency, especially at $z\gtrsim4.5$. We have developed a new method to select $4.7\lesssim z \lesssim 5.4$ quasars with both high efficiency and completeness by combining optical and mid-IR Wide-field Infrared Survey Explorer (WISE) photometric data, and are conducting a luminous $z\sim5$ quasar survey in the whole Sloan Digital Sky Survey (SDSS) footprint. We have spectroscopically observed 99 out of 110 candidates with $z$-band magnitudes brighter than 19.5 and 64 (64.6\%) of them are quasars with redshifts of $4.4\lesssim z \lesssim 5.5$ and absolute magnitudes of $-29\lesssim M_{1450} \lesssim -26.4$. In addition, we also observed 14 fainter candidates selected with the same criteria and identified 8 (57.1\%) of them as quasars with $4.7<z<5.4$ . Among 72 newly identified quasars, 12 of them are at $5.2 < z < 5.7$, which leads to an increase of $\sim$36\% of the number of known quasars at this redshift range. More importantly, our identifications doubled the number of quasars with $M_{1450}<-27.5$ at $z>4.5$, which will set strong constraints on the bright end of the quasar luminosity function. We also expand our method to select quasars at $z\gtrsim5.7$. In this paper we report the discovery of four new luminous $z\gtrsim5.7$ quasars based on SDSS-WISE selection.