Let $h_d$ be the number of $SL_{2}(\mathbb{Z})$ classes of primitive binary quadratic forms of discriminant $d$. It's natural to impose the hypothesis that $d$ is not at square, as we do below.
In Carl Ludwig Siegel's paper titled The Average Measure of Quadratic Forms With Given Discriminant and Signature Siegel cites two formulae given by Gauss in Disquisitiones Arithmeticae:
(a) $\displaystyle\sum\limits_{d= -N }^1 h_d \sim \frac{\pi}{18 \zeta(3)}N^{3/2}$
(b) $\displaystyle\sum\limits_{d = 1}^N h_d \log{\epsilon}_d \sim \frac{{\pi}^2}{18 \zeta(3)}N^{3/2}$
Where $N > 0$ and $\epsilon_{d} = \frac{1}{2}(t + u \sqrt{d})$ where $(t,u)$ is the smallest positive solution to $t^2 - ud^2 = 4$.
(Actually, Gauss restricts to consideration to binary quadratic forms with even middle coefficient correspondingly arrives at different formulae, but they're essentially the same as those above).
Siegel gives two proofs of these formulae: one proceeding from Dirichlet's class number formula together with character sum estimates due to Polya and Landau, and one via a direct lattice point counting argument.
In light of the facts that (i) I haven't heard anyone say that Gauss's was the one to discover the class number formula and (ii) the character sum estimates seem outside of the scope of Gauss's work, I imagine that his argument was via lattice point counting. Do we have any evidence otherwise? (I checked Gauss's book and he doesn't describe his methods there.)