Transverse momentum dependent (TMD) distributions depend on the pair of scaling parameters and their evolution is given by a pair of coupled equations. I present the analysis of the TMD evolution equations and their solution with the emphasis on their two-dimensional structure. It results in a new viewpoint on TMD evolution, both from the technical and interpretation sides. I formulate the non-perturbative definition of zeta-prescription and introduce the notion of optimal TMD distribution. I demonstrate that the updated form of TMD evolution produces lesser theoretical uncertainty and improves agreement with the data.