The East Asian continent was formed by multiple Precambrian cratons and orogenic belts through its geological history (Figure 1a). In the Mesozoic, the eastern part of East Asia deformed in response to the subduction of the Paleo-Pacific Plate (Li et al., 2019). Some stable blocks, such as the North China Craton (NCC), were thermally rejuvenated (Zhu et al., 2011). The Cenozoic tectonics of East Asia has been dominated by the India-Eurasian collision, which led to the formation of the Tibetan Plateau (Yin & Harrison, 2000). The diffuse deformation and volcanism make East Asia one of the best examples of complex intracontinental tectonics. Temperature, which controls rock rheology and magmatism, is a key for understanding such continental tectonics (Artemieva, 2006). However, the three-dimensional thermal structure of East Asian lithosphere remains uncertain.Direct measurements of lithospheric temperature at depth are not feasible. However, it is possible to indirectly estimate the temperature of continental lithosphere (Goes et al., 2000). In one approach, the lithospheric temperature is calculated by solving the steady-state heat conduction equation based on the near-surface heat flow, thermal conductivity, and heat production data (Artemieva & Mooney, 2001;Artemieva, 2019). This method is reasonably accurate for stable cratons where the assumption of steady-state conductive heat transfer holds. But errors may occur for Cenozoic orogenic belts, where steady-state conduction within the lithosphere may have not been reached. Compared with other methods, this method is more accurate for near-surface temperature. For the