  We introduce the concept of Random Sequential Renormalization (RSR) for
arbitrary networks. RSR is a graph renormalization procedure that locally
aggregates nodes to produce a coarse grained network. It is analogous to the
(quasi-)parallel renormalization schemes introduced by C. Song {\it et al.}
(Nature {\bf 433}, 392 (2005)) and studied more recently by F. Radicchi {\it et
al.} (Phys. Rev. Lett. {\bf 101}, 148701 (2008)), but much simpler and easier
to implement. In this first paper we apply RSR to critical trees and derive
analytical results consistent with numerical simulations. Critical trees
exhibit three regimes in their evolution under RSR: (i) An initial regime
$N_0^{\nu}\lesssim N<N_0$, where $N$ is the number of nodes at some step in the
renormalization and $N_0$ is the initial size. RSR in this regime is described
by a mean field theory and fluctuations from one realization to another are
small. The exponent $\nu=1/2$ is derived using random walk arguments. The
degree distribution becomes broader under successive renormalization --
reaching a power law, $p_k\sim 1/k^{\gamma}$ with $\gamma=2$ and a variance
that diverges as $N_0^{1/2}$ at the end of this regime. Both of these results
are derived based on a scaling theory. (ii) An intermediate regime for
$N_0^{1/4}\lesssim N \lesssim N_0^{1/2}$, in which hubs develop, and
fluctuations between different realizations of the RSR are large. Crossover
functions exhibiting finite size scaling, in the critical region $N\sim
N_0^{1/2} \to \infty$, connect the behaviors in the first two regimes. (iii)
The last regime, for $1 \ll N\lesssim N_0^{1/4}$, is characterized by the
appearance of star configurations with a central hub surrounded by many leaves.
The distribution of sizes where stars first form is found numerically to be a
power law up to a cutoff that scales as $N_0^{\nu_{star}}$ with
$\nu_{star}\approx 1/4$.
