The paper is devoted to a behavior investigation of the upper bound of deviation of functions from Zygmund classes from their biharmonic Poisson integrals. Systematic research in this direction was conducted by a number of Ukrainian as well as foreign scientists. But most of the known results relate to an estimation of deviations of functions from different classes from operators that were constructed based on triangular l-methods of the Fourier series summation (Fejer, Valle Poussin, Riesz, Rogozinsky, Steklov, Favard, etc.). Concerning the results relating to linear methods of the Fourier series summation, given by a set of functions of natural argument (Abel-Poisson, Gauss-Weierstrass, biharmonic and threeharmonic Poisson integrals), in this direction the progress was less notable. This may be due to the fact that the above-mentioned linear methods the Fourier series summation are solutions of corresponding integral and differential equations of elliptic type. And, therefore, they require more time-consuming calculations in order to obtain some estimates, that are suitable for a direct use for applied purposes. At the same time, in the present paper we investigate approximative characteristics of linear positive Poisson-type operators on Zygmund classes of functions. According to the well-known results by P.P. Korovkin, these positive linear operators realize the best asymptotic approximation of functions from Zygmund classes. Thus, the estimate obtained in this paper for the deviation of functions from Zygmund classes from their biharmonic Poisson integrals (the least studied and most valuable among all linear positive operators) is relevant from the viewpoint of applied mathematics.