The discovery of multiply-imaged gravitationally lensed QSOs is fundamental to many astronomical and cosmological studies. However, these objects are rare and challenging to discover due to requirements of high-angular resolution astrometric, multiwavelength photometric and spectroscopic data. This has limited the number of known systems to a few hundred objects. We aim to reduce the constraints on angular resolution and discover multiply-imaged QSO candidates by using new candidate selection principles based on unresolved photometric time-series and ground-based images from public surveys. We selected candidates for multiply-imaged QSOs based on low levels of entropy computed from Catalina unresolved photometric time-series or Euclidean similarity to known lenses in a space defined by the wavelet power spectra of Pan-STARSS DR2 or DECaLS DR7 images, combined with multiple {\it Gaia} DR2 sources or large astrometric errors and supervised and unsupervised learning methods. We then confirmed spectroscopically some candidates with the Palomar Hale, Keck-I, and ESO/NTT telescopes. Here we report the discovery and confirmation of seven doubly-imaged QSOs and one likely double quasar. This demonstrates the potential of combining space-astrometry, even if unresolved, with low spatial-resolution photometric time-series and/or low-spatial resolution multi-band imaging to discover multiply-imaged lensed QSOs.

We have used spectra of 181 projected quasar pairs at separations $\le 1.5$ arcmins from the Sloan Digital Sky-Survey Data Release 12 in the redshift range of 2.5 to 3.5 to probe the proximity regions of the foreground quasars. We study the proximity effect both in the longitudinal as well as in the transverse directions, by carrying out a comparison of the \lya absorption lines originating from the vicinity of quasars to those originating from the general inter-galactic medium at the same redshift. We found an enhancement in the transmitted flux within 4 Mpc to the quasar in the longitudinal direction. However, the trend is found to be reversed in the transverse direction. In the longitudinal direction, we derived an excess overdensity profile showing an excess up to $r\leq 5$ Mpc after correcting for the quasar's ionization, taking into account the effect of low spectral resolution. This excess overdensity profile matches with the average overdensity profile in the transverse direction without applying any correction for the effect of the quasar's ionization. Among various possible interpretations, we found that the anisotropic obscuration of the quasar's ionization seems to be the most probable explanation. This is also supported by the fact that all of our foreground quasars happen to be Type-I AGNs. Finally, we constrain the average quasar's illumination along the transverse direction as compared to that along the longitudinal direction to be $\leq$27\% (3$\sigma$ confidence level).