Background: Estimating forest productivity is critical for effective management and site assessment. The dominant height is used to calculate the Site Index (SI), which is commonly used to assess forest productivity. In this study, an algebraic difference approach was used to develop a dominant height model incorporating the rainfall effect for Eucalyptus grandis x Eucalyptus urophylla (E. grandis x E. urophylla). The dataset consists of 75 permanent sample plots ranging in age from 0.5 to 11 years, as well as 7 rainfall stations spread across plantations in Coastal Zululand, South Africa. Using fixed and mixedeffects in the predictor function, twelve candidate models were derived from the Bertalanffy-Richards, Lundqvist-Korf, McDill-Amateis, and Hossfeld growth functions. A continuous-time autoregressive error structure was used to account for serial autocorrelation in the longitudinal unbalanced data. Model fit statistics and graphical methods were used to evaluate the candidate models. Results:The addition of the rainfall effect increased model precision by 37%. The mixed-effects formulation produced 18% more precision when compared to similar models with all parameters fixed. Due to their compatibility with expected biological behaviour and good performance on validation data, mixed-effects models based on Lundqvist-Korf and McDill-Amateis functions were chosen as the final models. Conclusion:Unlike similar models that do not take rainfall into account, these models can capture the effects of severe rainfall conditions such as drought and can thus be used in short-rotation pulp forests with fluctuating rainfall.