Biomass additivity is a desirable characteristic of a system of equations for predicting components and total biomass, since equations independently adjusted generate biologically inconsistent results. The aim of this study was to fit and compare three methods for modelling biomass: (i) total biomass individual regression, (ii) total biomass regression function calculated as the sum of separate biomass components, and (iii) simultaneous equations of biomass components based on Nonlinear Seemingly Unrelated Regression. A total of 208 trees of Eucalyptus dunnii and Eucalyptus grandis were harvested and destructively sampled to record above-ground biomass. Results indicate that a system of equations adjusted by simultaneous equations provides accurate biomass estimations, guaranteeing additivity. This model system showed good fit and good prediction performance, given that the correlation coefficient was higher than 97 % for total above-ground biomass, for both species; whereas root mean square error was 23.9 kg and 30.2 kg for E. grandis and E. dunnii, respectively. A system of biomass equations was developed for each eucalyptus species, such that the sum of the estimations of the biomass components equaled the estimate of above-ground biomass. Results showed that the systems of equations have high potential for improving the accuracy of individual tree above-ground biomass estimates for both species.La aditividad de las ecuaciones empeladas para predecir los componentes y la biomasa total es una característica deseable de un sistema de ecuaciones. Las ecuaciones ajustadas independientemente generan resultados biológicamente inconsistentes. El objetivo de este estudio fue ajustar y comparar tres métodos para estimar la biomasa de rodales de Eucalyptus: (i) regresión clásica individual, (ii) estimación de biomasa total calculada como la suma de las ecuaciones individuales de las fracciones, y (iii) ecuaciones simultáneas aparentemente no relacionadas. Un total de 208 árboles de Eucalyptus dunnii y Eucalyptus grandis fueron cosechados y muestreados destructivamente para registrar la biomasa aérea. Los resultados indican que el sistema de ecuaciones simultáneas proporciona estimaciones precisas de biomasa, garantizando la aditividad. Este sistema mostró un alto ajuste y desempeño de predicción, dado que registró un coeficiente de determinación mayor al 97 % para la biomasa total, en ambas especies, mientras que la raíz del error cuadrático medio fue 23,9 kg y 30,2 kg para E. grandis y E. dunnii, respectivamente. Se desarrolló un sistema de ecuaciones de biomasa para cada especie de eucalipto, de modo que la suma de las estimaciones de las fracciones equivalía a la estimación de biomasa total. Los resultados muestran que los sistemas de ecuaciones tienen un gran potencial para mejorar la precisión de las estimaciones individuales de biomasa aérea de árboles para ambas especies.