Maximal Seshadri constants on arbitrary polarized surfaces

We prove that, for every smooth integral complex projective surface S and every ample line bundle L, the multipoint Seshadri constant at r very general points equals $\sqrt{L^2/r}$ for every sufficiently large integer r. This resolves positively the qualitative Nagata–Biran conjecture for surfaces.