Gravitational time delays provide a powerful one-step measurement of H0, independent of all other probes. One key ingredient in time-delay cosmography are high-accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (105–106 CPU hours and ∼0.5–1 yr, respectively). Major improvements in modelling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, we develop an automated modelling pipeline and apply it to a sample of 31 lens systems, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in the Fermat potential (proportional to time delay) with respect to modelling choices. We find that for 10/30 lenses, our models are cosmography or nearly cosmography grade (<3 per cent and 3–5 per cent variations). For 6/30 lenses, the models are close to cosmography grade (5–10 per cent). These results utilize informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modelling sequence and options. In conclusion, we show that uniform cosmography grade modelling of large strong lens samples is within reach.

ABSTRACT We report the results of the STRong lensing Insights into the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates with ESO Faint Object Spectrograph and Camera 2 on the NTT and Echellette Spectrograph and Imager on Keck, confirming 10 new lensed quasars and 10 quasar pairs. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged (DESJ0345−2545, z = 1.68), and one is quadruply imaged (quad: DESJ0053−2012, z = 3.8). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on Southern Astrophysical Research Telescope or Near InfraRed Camera 2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcsec. After spectroscopic follow-up, we extract multi-epoch grizY photometry of confirmed lensed quasars and contaminant quasar + star pairs from DES data using parametric multiband modelling, and compare variability in each system’s components. By measuring the reduced χ2 associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems. Based on our spectroscopic follow-up, this variability information improves selection of lensed quasars and quasar pairs from 34-45 per cent to 51–70 per cent, with most remaining contaminants being star-forming galaxies. Using mock lensed quasar light curves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample, explained by the magnification bias and hence lower luminosity/more variable sources in quads.