This paper examines the convergence in $L_{\lp}$, for $1\leq \lp \leq 2$, and almost surely of the numbers in a varying environment Galton-Watson process with a countable set of types. The harmonic functions of the mean matrices play a central role in the analysis. Many previously studied models provide particular cases.