SPECULOOS-South, an observatory composed of four independent 1m robotic telescopes, located at ESO Paranal, Chile, started scientific operation in January 2019. This Southern Hemisphere facility operates as part of SPECULOOS, an international network of 1m-class telescopes surveying for transiting terrestrial planets around the nearest and brightest ultra-cool dwarfs. To automatically and efficiently process the observations of SPECULOOS-South, and to deal with the specialised photometric requirements of ultra-cool dwarf targets, we present our automatic pipeline. This pipeline includes an algorithm for automated differential photometry and an extensive correction technique for the effects of telluric water vapour, using ground measurements of the precipitable water vapour. Observing very red targets in the near-infrared can result in photometric systematics in the differential lightcurves, related to the temporally-varying, wavelength-dependent opacity of the Earth's atmosphere. These systematics are sufficient to affect the daily quality of the lightcurves, the longer time-scale variability study of our targets and even mimic transit-like signals. Here we present the implementation and impact of our water vapour correction method. Using the 179 nights and 98 targets observed in the I+z' filter by SPECULOOS-South since January 2019, we show the impressive photometric performance of the facility (with a median precision of ~1.5 mmag for 30-min binning of the raw, non-detrended lightcurves) and assess its detection potential. We compare simultaneous observations with SPECULOOS-South and TESS, to show that we readily achieve high-precision, space-level photometry for bright, ultra-cool dwarfs, highlighting SPECULOOS-South as the first facility of its kind.

We report the discovery of NGTS-4b, a sub-Neptune-sized planet transiting a 13th magnitude K-dwarf in a 1.34d orbit. NGTS-4b has a mass M=$20.6\pm3.0$M_E and radius R=$3.18\pm0.26$R_E, which places it well within the so-called "Neptunian Desert". The mean density of the planet ($3.45\pm0.95$g/cm^3) is consistent with a composition of 100% H$_2$O or a rocky core with a volatile envelope. NGTS-4b is likely to suffer significant mass loss due to relatively strong EUV/X-ray irradiation. Its survival in the Neptunian desert may be due to an unusually high core mass, or it may have avoided the most intense X-ray irradiation by migrating after the initial activity of its host star had subsided. With a transit depth of $0.13\pm0.02$%, NGTS-4b represents the shallowest transiting system ever discovered from the ground, and is the smallest planet discovered in a wide-field ground-based photometric survey.

We report the discovery and tomographic detection of WASP-174b, a planet with a near-grazing transit on a 4.23-d orbit around a $V$ = 11.9, F6V star with [Fe/H] = 0.09 $\pm$ 0.09. The planet is in a moderately misaligned orbit with a sky-projected spin-orbit angle of $\lambda$ = 31$^{\circ}$ $\pm$ 1$^{\circ}$. This is in agreement with the known tendency for orbits around hotter stars to be misaligned. Owing to the grazing transit the planet's radius is uncertain, with a possible range of 0.8-1.8 R$_{\rm Jup}$. The planet's mass has an upper limit of 1.3 M$_{\rm Jup}$. WASP-174 is the faintest hot-Jupiter system so far confirmed by tomographic means.

We present here SPECULOOS, a new exoplanet transit search based on a network of 1m-class robotic telescopes targeting the $\sim$1200 ultracool (spectral type M7 and later) dwarfs bright enough in the infrared ($K$-mag $\leq 12.5$) to possibly enable the atmospheric characterization of temperate terrestrial planets with next-generation facilities like the $\textit{James Webb Space Telescope}$. The ultimate goals of the project are to reveal the frequency of temperate terrestrial planets around the lowest-mass stars and brown dwarfs, to probe the diversity of their bulk compositions, atmospheres and surface conditions, and to assess their potential habitability.

We present the discovery of NGTS-3Ab, a hot Jupiter found transiting the primary star of an unresolved binary system. We develop a joint analysis of multi-colour photometry, centroids, radial velocity (RV) cross-correlation function (CCF) profiles and their bisector inverse slopes (BIS) to disentangle this three-body system. Data from the Next Generation Transit Survey (NGTS), SPECULOOS and HARPS are analysed and modelled with our new blendfitter software. We find that the binary consists of NGTS-3A (G6V-dwarf) and NGTS-3B (K1V-dwarf) at <1 arcsec separation. NGTS-3Ab orbits every P = 1.675 days. The planet radius and mass are R_planet = 1.48+-0.37 R_J and M_planet = 2.38+-0.26 M_J, suggesting it is potentially inflated. We emphasise that only combining all the information from multi-colour photometry, centroids and RV CCF profiles can resolve systems like NGTS-3. Such systems cannot be disentangled from single-colour photometry and RV measurements alone. Importantly, the presence of a BIS correlation indicates a blend scenario, but is not sufficient to determine which star is orbited by the third body. Moreover, even if no BIS correlation is detected, a blend scenario cannot be ruled out without further information. The choice of methodology for calculating the BIS can influence the measured significance of its correlation. The presented findings are crucial to consider for wide-field transit surveys, which require wide CCD pixels (>5 arcsec) and are prone to contamination by blended objects. With TESS on the horizon, it is pivotal for the candidate vetting to incorporate all available follow-up information from multi-colour photometry and RV CCF profiles.

A high accuracy photometric reduction method is needed to take full advantage of the potential of the transit method for the detection and characterization of exoplanets, especially in deep crowded fields. In this context, we present DECPHOT, a new deconvolution-based photometry algorithm able to deal with a very high level of crowding and large variations of seeing. It also increases the resolution of astronomical images, an important advantage for the discrimination of false positives in transit photometry.

A new method is presented for determining the Point Spread Function (PSF) of images that lack bright and isolated stars. It is based on the same principles as the MCS (Magain, Courbin, Sohy, 1998) image deconvolution algorithm. It uses the information contained in all stellar images to achieve the double task of reconstructing the PSFs for single or multiple exposures of the same field and to extract the photometry of all point sources in the field of view. The use of the full information available allows to construct an accurate PSF. The possibility to simultaneously consider several exposures makes it very well suited to the measurement of the light curves of blended point sources from data that would be very difficult or even impossible to analyse with traditional PSF fitting techniques. The potential of the method for the analysis of ground-based and space-based data is tested on artificial images and illustrated by several examples, including HST/NICMOS images of a lensed quasar and VLT/ISAAC images of a faint blended Mira star in the halo of the giant elliptical galaxy NGC5128 (Cen A).

A high accuracy photometry algorithm is needed to take full advantage of the potential of the transit method for the characterization of exoplanets, especially in deep crowded fields. It has to reduce to the lowest possible level the negative influence of systematic effects on the photometric accuracy. It should also be able to cope with a high level of crowding and with large scale variations of the spatial resolution from one image to another. A recent deconvolution-based photometry algorithm fulfills all these requirements, and it also increases the resolution of astronomical images, which is an important advantage for the detection of blends and the discrimination of false positives in transit photometry. We made some changes to this algorithm in order to optimize it for transit photometry and used it to reduce NTT/SUSI2 observations of two transits of OGLE-TR-113b. This reduction has led to two very high precision transit light curves with a low level of systematic residuals, used together with former photometric and spectroscopic measurements to derive new stellar and planetary parameters in excellent agreement with previous ones, but significantly more precise.

We present a new method to estimate time delays from light curves of lensed quasars. The method is based on chi^2 minimization between the data and a numerical model light curve. A linear variation can be included in order to correct for slow long-term microlensing effects in one of the lensed images. An iterative version of the method can be applied in order to correct for higher order microlensing effects. The method is tested on simulated light curves. When higher order microlensing effects are present the time delay is best constrained with the iterative method. Analysis of a published data set for the lensed double Q0957+561 yields results in agreement with other published estimates.

A method for spatial deconvolution of spectra is presented. It follows the same fundamental principles as the ``MCS image deconvolution algorithm'' (Magain, Courbin, Sohy, 1998) and uses information contained in the spectrum of a reference Point Spread Function (PSF) to spatially deconvolve spectra of very blended sources. An improved resolution rather than an infinite one is aimed at, overcoming the well known problem of ``deconvolution artefacts''. As in the MCS algorithm, the data are decomposed into a sum of analytical point sources and a numerically deconvolved background, so that the spectrum of extended sources in the immediate vicinity of bright point sources may be accurately extracted and sharpened. The algorithm has been tested on simulated data including seeing variation as a function of wavelength and atmospheric refraction. It is shown that the spectra of severely blended point sources can be resolved while fully preserving the spectrophotometric properties of the data. Extended objects ``hidden'' by bright point sources (up to 4-5 magnitudes brighter) can be accurately recovered as well, provided the data have a sufficiently high total signal-to-noise ratio (200-300 per spectral resolution element). Such spectra are relatively easy to obtain, even down to faint magnitudes, within a few hours of integration time with 10m class telescopes.

A new method for improving the resolution of astronomical images is presented. It is based on the principle that sampled data cannot be fully deconvolved without violating the sampling theorem. Thus, the sampled image should not be deconvolved by the total Point Spread Function, but by a narrower function chosen so that the resolution of the deconvolved image is compatible with the adopted sampling. Our deconvolution method gives results which are, in at least some cases, superior to those of other commonly used techniques: in particular, it does not produce ringing around point sources superimposed on a smooth background. Moreover, it allows to perform accurate astrometry and photometry of crowded fields. These improvements are a consequence of both the correct treatment of sampling and the recognition that the most probable astronomical image is not a flat one. The method is also well adapted to the optimal combination of different images of the same object, as can be obtained, e.g., from infrared observations or via adaptive optics techniques.