  We compute the two-loop $\beta$-function of scalar and spinorial quantum electrodynamics as well as pure Yang-Mills and quantum chromodynamics using the background field method in a fully quadridimensional setup using Implicit Regularization (IREG). Moreover, a thorough comparison with dimensional approaches such as conventional dimensional regularization (CDR) and dimensional reduction (DRED) is presented. Subtleties related to Lorentz algebra contractions/symmetric integrations inside divergent integrals as well as renormalisation schemes are carefully discussed within IREG where the renormalisation constants are fully defined as basic divergent integrals to arbitrary loop order. Moreover, we confirm the hypothesis that momentum routing invariance in the loops of Feynman diagrams implemented via setting well-defined surface terms to zero deliver non-abelian gauge invariant amplitudes within IREG just as it has been proven for abelian theories. 
Two-loop renormalisation of gauge theories in $4D$ Implicit   Regularisation and connections to dimensional methods