EDIT___(UPDATE ON QUESTION 2)
Ira Gessel has now carried out computer calculations for all l<1500; here's a summary of his remarkable results. Consider the triples (r1,r2,r3) where the squares of r1,r1,r2 and r3 sum to l. To avoid duplicates, normalize each such triple so that r1,r2 and r3 are positive with r2>r3. Ira finds that for each such triple there is a unique second (normalized) triple (s1,s2,s3) such that the power series C(r1,r1,r2,r3)+C(2s1,s2+s3,s2-s3)+C is divisible by x^(l^2). Furthermore (r1,r2,r3)-->(s1,s2,s3) is an involution on the set of normalized triples with at most 1 fixed point.
When l=7 mod 16, the argument I gave above when l=23 then shows that for each normalized triple (r1,r2,r3) and the corresponding (s1,s2,s3) we have the identity C=C(r1,r1,r2,r3)+C(2s1,s2+s3,s2-s3); note that the squares of 2s1, s2+s3, and s2-s3 sum to 2l. So, for example, when l=1447 we get 17 distinct formulae for C.
When l=15 mod 16, it seems certain that once again C(r1,r1,r2,r3) and C(s1,s2,s3) sum to C. This could be proved by extending Ira's calculations to prove divisibility by x^(d+1) where d=l(l+1)(l+1), as in my treatment of l=31 and l=47. But I think this extension is unnecessary, and that one may deduce the identities simply from the divisibility by x^(l^2) established by Ira. If my idea for demonstrating this works out I'll post it as a comment.
But great mysteries remain. Why should this all be true? And can one describe the mysterious involution (r1,r2,r3)-->(s1,s2,s3) explicitly? By the way, it's known that the number of normalized triples is odd or even according as n is 7 or 15 mod 16. The proof of this goes back to Hasse; one shows that the number of triples is h/4 where h is the class number of Q(Root(-2l)) and uses results of Gauss on representations by sums of 3 squares, together with some genus theory for binary quadratic forms.