Conditional returns distributions generated by a GARCH process, which are important for many applications in market risk assessment and portfolio optimization, are typically generated via simulation. This paper extends previous research on analytic moments of GARCH returns distributions in several ways: we consider a general GARCH model -- the GJR specification with a generic innovation distribution; we derive analytic expressions for the first four conditional moments of the forward return, of the forward variance, of the aggregated return and of the aggregated variance -- corresponding moments for some specific GARCH models largely used in practice are recovered as special cases; we derive the limits of these moments as the time horizon increases, establishing regularity conditions for the moments of aggregated returns to converge to normal moments; and we demonstrate empirically that some excellent approximate predictive distributions can be obtained from these analytic moments, thus precluding the need for time-consuming simulations.