In this paper we study the property of asymptotic completeness in (massive) Euclidean lattice quantum field theories. We use the methods of Spencer and Zirilli [2] to prove, under suitable hypothesis, two-body asymptotic completeness, i.e., for the energy range just above the two-particle threshold.This paper is dedicated to the problem of asymptotic completeness (AC) in Euclidean lattice quantum field theories. In our analysis we followed closely the methods of Spencer and Zirilli [2], who proved, probably for the first time, asymptotic completeness for a true relativistic quantum field theoretical model for a limited range of energies. In spite of the higher technical difficulties posed by the lack of Lorentz invariance and other problems we are able to reproduce the same results of [2] in our lattice context. We have namely proven two-body asymptotic completeness, i.e., asymptotic completeness up to energies just above the two-particle threshold. These methods rely basically on the exponential decay of the Bethe-Salpeter kernel in space variables, and can be applied to the analysis of bound states and resonances for weakly coupled models, as in [15,16] in the continuum, or [44,47,45,46] in lattice models. Also, and with more work, the methods can be extended to the analysis of three-particle bound states and AC following, for instance, [17] and [13], respectively. Limitations on the energy range are found, unfortunately, in all the proofs of asymptotic completeness performed along these lines on the continuum. This deplorable state of affairs urges the QFT community to develop new ideas and techniques to deal with the spectral problems of QFT, but this is not our subject here. See, however, [8] (and also [9, 10]) for a proposal involving the nuclearity criterion. For a discussion of the problem of asymptotic completeness in QED in the context of perturbation theory, see [11].We believe our results are interesting not only due to their connection to QFT problems, as mentioned above. Some of the models we consider are, in fact, models of classical statistical mechanical spin systems (for instance, the Ising model) and the spectral properties of the transfer matrix reflect on corrections to the exponential decay of correlations as, for instance, the Ornstein-Zernike corrections (see [22,33,34,35]).In our work many adaptations to the lattice context were necessary, which increased considerably the technical complications involved. Let us briefly discuss some of them. In their work, for instance, Spencer and Zirilli [2] restricted a good part of their analysis to the zero-momentum sector of the energy-momentum spectrum and then invoked Lorentz covariance to extend their results to nonzero momenta. This strategy of argumentation simplifies many computations but cannot be applied to a situation where Lorentz covariance is lacking. Another major source of complications involved a series of space-time changes of variables intended to express some four point functions and the Bethe-Salpeter kernel in terms of "centre of...