While reading a paper by Owens (arXiv:1906.05913) about embeddings of rational homology balls in the complex projective plane, I found out the following somewhat unexpected number theory corollary (here $F(n)$ denotes the $n$-th Fibonacci number, with $F(1)=1$ and $F(2)=1$).
If $n>1$, the square of the odd Fibonacci number $F(2n+1)$ can be written as the sum of exactly $F(2n+1)+1$ nonzero squares.
I will put an answer with the proof. However, I am curious about whether a proof of this result can be obtained by purely numer-theoretic methods.