Global quantum discord (GQD) in two‐dimensional (2D) transverse‐field Ising models on N × N square lattices at zero temperature and finite temperatures are studied. First, the GQD Gtrue(ρ2×2true) defined on the reduced density matrix ρ2×2 of a 2 × 2 sub‐square in the center of the N × N lattices, with N=4, 5, 6, 7, ∞ is studied. It is found that Gtrue(ρ2×2true) shows little size dependence in the strong‐field and weak‐field limit, and presents dramatic size dependence in the middle‐field region. The results are explained by short‐range correlations in non‐critical regions and long‐range correlations in critical regions of the models, respectively. Then the scaling behavior of the GQD Gtrue(true|ψg〉true) of the ground states for the entire N × N squares is studied. It is found that Gtrue(true|ψg〉true) contains a 1D contribution (corresponding to the side length of the squares) and a 2D contribution (corresponding to the area of the squares). Furthermore, the effect of thermodynamic fluctuations on GQD by investigating the thermal‐state GQD Gtrue(e−HkBTtrue) is considered. In high‐temperature regions, GQD is always destroyed by thermodynamic fluctuations for any strength of the magnetic field. However, in low‐temperature regions, some thermal enhancement of GQD under weak fields and some thermal robustness of GQD under strong fields is observed. The results are explained by the low‐lying energy states and the energy gap of the models.