This paper describes the development of a method to characterize the temperature profile of silicon wafers in high-temperature environments. Monocrystalline wafers are implanted on one surface with B and P ions, which diffuse into the wafer at different rates based on the temperature-dependent diffusivity of the ions during a 30 min soak in the high-temperature environment between 1000 and 1400 °C. The use of two different dopant species, instead of one, yields higher sensitivity of the measured resistance to changes in temperature in this high-temperature range. The ρ-T relation is simulated using TSUPREM-4 and calibrated using a furnace of known temperature. The final sheet resistance varies spatially between 50 and 250 Ω/sq, and can be related to the temperature of each part of the wafer during the soak step, with a sensitivity of ∼0.5 Ω/°C, and a two-dimensional temperature map can be extracted. The method is demonstrated on wafers by characterizing the hot zone of a high-temperature furnace.