Using continuum extrapolated lattice data we trace a family of running couplings in three-flavour QCD over a large range of scales from about 4 to 128 GeV. The scale is set by the finite space time volume so that recursive finite size techniques can be applied, and Schrödinger functional (SF) boundary conditions enable direct simulations in the chiral limit. Compared to earlier studies we have improved on both statistical and systematic errors. Using the SF coupling to implicitly define a reference scale 1 MS , a reliable estimate of the truncation error requires non-perturbative data for a sufficiently large range of values of α s =ḡ 2 /(4π). In the present work we reach this precision by studying scales that vary by a factor 2 5 = 32, reaching down to α s ≈ 0.1. We here provide the details of our analysis and an extended discussion. a