I think if you try this a bit, you'll find that to get the points where the degree is lowest, you can just apply the same permutation to the points for every value of $i$. The procedure of making a Cremona transformation at points $p_7$, $p_8$, and $p_9$, and then reordering the points so those come first, will achieve the sequence of transformations that you want.
Once you do that, to get the degree and multiplicity of $C_i$ is just a matter of applying some explicit $10 \times 10$ matrix a bunch of times. This matrix will have a $3 \times 3$ Jordan block associated to the eigenvalue $1$ (probably, I didn't check), and so you get the growth of degrees that you want. If you use 10 points instead, it will have an eigenvalue greater than $1$.
To illustrate, here are the first few $C_i$. I write "$d;m_1,m_2,\dots,m_9$" for a curve of degree $d$ with multiplicities $m_i$. (Sorry, don't know how to make a table.)
Each row I get by making a Cremona transformation at the last three points, computing the new degree, and then moving those three points to the beginning. This has the convenient property that I always use the points of least multiplicity, like you want. Going from one row to the next is a linear map (the one saying what the new degree/mults are, composed with the one for permuting). The degrees are growing quadratically, and figuring out an explicit formula is simple linear algebra.
You can also do this using facts about the Weyl group of a certain root system. The Asterisque volume by Dolgachev and Ortland is a good source. Another good reference is "Automorphisms of Blow-ups of $\mathbf P^2$" by McMullen, where you will see a similar procedure. (The application of the "Coxeter element" there corresponds to doing the same thing, but with a different permutation rule, so it won't always hit the three points of lowest multiplicity. Nevertheless the degrees of the $C_i$ grow unboundedly.)