Chenxiao Yu, Min Guo, Guo Qing, Dongwei Zhao, Lechi Zhang, Zhenyu Xu, Anjie Cao, Jie Yang, Wensheng Lin, Wenchi Cheng, Qinghe Du, Li Li
As wireless networks evolve from 5G toward 6G, the complexity of the electromagnetic environment increases sharply. Spectrum usage expands significantly into millimetre-wave (mmWave) and terahertz (THz) high-frequency bands. Network node density and mobility increase markedly. Moreover, communication-sensing-computation functions are deeply integrated. Accurate, real-time, full-band Electromagnetic Spectrum Maps (ESMs) have become a core infrastructure for 6G spectrum situational awareness, Dynamic Spectrum Access (DSA), interference coordination, and Integrated Sensing and Communication (ISAC). However, while a growing body of recent work extends radio mapping to multi-band and temporal domains, the predominant focus of existing Radio Map research remains the two-dimensional spatial power distribution at a single fixed frequency—essentially a degenerate special case of ESM after the frequency and time dimensions are collapsed—and no existing survey unifies 3D spatial construction, time-varying prediction, and full 6G system integration under a shared 4D formalism. This paper focuses on the three core research dimensions of ESMs, i.e., 3D spatial ESM construction, dynamic time-varying ESM modelling and prediction, and ESM integration with 6G systems. Under a unified four-dimensional ESM framework (space × frequency × time × power), we clarify the hierarchical relationships among ESM/SEM/REM/Radio Map/Channel Knowledge Maps (CKMs). Then, we systematically review 3D ESM construction, dynamic ESM modelling and prediction, and the integration of ESM with CKM/Digital Twin Networks (DTNs)/ISAC. Finally, we identify five, core open problems that constrain the development of the field to provide a systematic reference for 6G intelligent spectrum management research.