In this paper, following the idea of Glavan and Lin [Phys. Rev. Lett. 124, 081301 (2020)] of rescaling the Gauss-Bonnet (GB) coupling constant to α/(D − 4), we obtain the Bardeen black hole solution of 4D Einstein-Gauss-Bonnet gravity. We analyze the horizon structure to determine the effect of GB parameter α on the minimum and maximum cutoff values of mass, M 0 , and magnetic monopole charge, g 0 , for the existence of a black hole horizon. The analytical expressions for thermodynamical quantities, namely Hawking temperature T + , entropy S + , Gibbs free energy F + and specific heat C + associated with the black hole horizon are determined. In the analysis of heat capacity, we find that there exists a critical value of horizon radius, r c + , corresponding to the local maximum of Hawking temperature, at which heat capacity diverges, confirming the second-order phase transition.