Agroforestry land-use systems in the Andean region have great socioeconomical and biophysical relevance due to the abundance of products and services they provide. Biomass estimation in these systems constitutes a priority concern as it facilitates assessment of carbon sink potential and functionality for biomass production. In this paper, a set of equations were fitted to enable easy and reliable estimation of the total aboveground biomass of four frequently used species in Andean agroforestry systems: Acacia melanoxylon L., Alnus acuminata Kunth., Buddleja coriacea Remy. and Polylepis racemosa Ruiz&Pav. The best models for each biomass component (stem, thick branches, thin branches and leaves) per species were fitted simultaneously according to SUR methodology (seemingly unrelated regressions). All models showed high goodness of fit statistics and more than 70% of the observed variation in biomass components was explained by the independent variables. The inclusion of height as a predictive variable in the models improved their predictive reliability and expanded the application range. The models developed here are useful for assessing the sustainability of agroforestry systems and could support governmental or non-governmental forest conservation incentive programs and initiatives.