Gadolinium (Gd) has been proven to be such a complex magnetic system that its static universality class varies as a function of the reduced temperature (|ɛ| = |T−TC|/TC). Though many works have been done, physical mechanisms of this phenomenon have not been fully elucidated. In this work, we analyzed the static magnetization data, M(H), around the ferromagnetic-paramagnetic phase transition temperature (TC) based on Banerjee's criteria, the modified Arrott plot and Kouvel–Fisher techniques, the critical isotherm analysis, and the scaling hypothesis to reconsider the phase transition, magnetic order, and isotropic and uniaxial behaviors of a polycrystalline Gd bulk. The M(H) data analyses with H variations up to 140 kOe indicated that Gd still keeps the nature of a second-order phase transition and short-range magnetic order. In the critical region of |ɛ| = 3.4 × 10−3–6.9 × 10−2, we observed the critical-behavior transition from a mixed isotropic and uniaxial ferromagnet for H = 3–10 kOe through an isotropic 3D Heisenberg ferromagnet for H = 10–100 kOe to a uniaxial 3D Ising ferromagnet for H = 100–140 kOe. The results demonstrate the possibility of using high fields to identify the universality class of Gd. The fluctuations and preferred reorientation of spins in the hexagonal close-packed lattice at temperatures T > 230 K are thought to cause these interesting phenomena.