In this paper, a novel massively space-time parallel finite element (FE) transformer model is proposed for power system electromagnetic transient simulations. The method utilizes the reduced magnetic vector potential (RMVP) formulation to solve the electromagnetic field in the transformer, and the 3-D tetrahedral edge element is applied to discretize the RMVP formulation. The discretized nonlinear system is solved by adaptive transmission line modeling (TLM) and matrix-free conjugate gradient iterator, which allows parallelism at the elemental-level and eliminates the need for Newton-Raphson (NR) iterations over the large-scale global matrix. The paper also integrates the parallel-in-time method with the finite element solver to further parallelize the model at space-time level, and the GPU implementation of the model realized a speedup of over 55 and high accuracy, compared with a commercial FE software.