With unique atomic stacking structures and outstanding performances, amorphous alloys have been formulated during the past several decades via various solidification techniques. [1,2] The cooling rates required to produce amorphous alloys are generally as high as 10 4 -10 5 °C s À1 . [3,4] This imposes restrictions on manufacturing amorphous alloys into giant-sized and geometrically complex parts. Nevertheless, amorphous alloys exhibit viscous-like behavior in the supercooled liquid region (SLR). Micron-sized amorphous powders can be fabricated from most alloy systems, even those having low glass-forming ability. [4] These provide a practical approach to densifying amorphous powders via powder metallurgy (P/M) processes. [5,6] Among the various P/M techniques available, spark plasma sintering (SPS) is promising for consolidating amorphous powders due to its fast heating rate and short cycle time. [5] In the past decade, consolidation of amorphous powders using the SPS process has been reported for Zr-based, [7,8] Ti-based, [9][10][11] Fe-based, [12] and Cu-based [13] alloys.The amorphous powders are heated to the SLR and densified by viscous flow. [10] However, the crystallization tendency is more intense with the increased temperature. [14] Therefore, the research focuses on fast densification while maintaining the amorphous structure. Amorphous powders' densification and crystallization dynamics race each other. Temperature is the key factor affecting this competitive process. Local overheating can lead to partial crystallization, resulting in severe embrittlement. For example, Nowak et al. reported that excessive overheating is generated in the vicinity of particle necks during the initial stages and is responsible for partial devitrification. [15] Similar phenomena have been widely observed in Fe-based [16] and Al-based amorphous alloys. [17] Despite the significant advantages of SPS, its further popularization is hindered by severe temperature and stress gradients, [18] especially for complex-shaped and large-sized parts. Therefore, it is of great concern to figure out the temperature field evolution in the amorphous alloys prepared by SPS. However, to the authors' knowledge, the current understanding of the temperature distribution in amorphous alloys prepared by SPS is limited.The temperature gradients in the sintering molds and samples can be divided into inhomogeneity at the macro and particle scales. In terms of the macroscopic temperature field, Song et al. developed a coupled thermal-electric-mechanical model to simulate the SPS process of Fe powders, [19] which concerns