In this article, we introduce, for the first time, multivariate symmetrized and perturbed hyperbolic tangent-activated convolution-type operators in three forms. We present their approximation properties, that is, their quantitative convergence to the unit operator via the multivariate modulus of continuity. We continue with the multivariate global smoothness preservation of these operators. We present, in detail, the related multivariate iterative approximation, as well as, multivariate simultaneous approximation, and their combinations. Using differentiability in our research, we produce higher rates of approximation, and multivariate simultaneous global smoothness preservation is also achieved.