Celestial surfaces and quadratic forms

We start with a classical problem of Kummer and Darboux about describing surfaces that contain many circles and explain an answer that uses polynomial solutions of quadratic forms over arbitrary fields.