This paper presents a new heroic computing method for unstructured, low-order, finite-element, implicit nonlinear wave simulation: 1.97 PFLOPS (18.6% of peak) was attained on the full K computer when solving a 1.08T degrees-of-freedom (DOF) and 0.270T-element problem. This is 40.1 times more DOF and elements, a 2.68-fold improvement in peak performance, and 3.67 times faster in time-to-solution compared to the SC14 Gordon Bell finalist's state-ofthe-art simulation. The method scales up to the full K computer with 663,552 CPU cores with 96.6% sizeup efficiency, enabling solving of a 1.08T DOF problem in 29.7 s per time step. Using such heroic computing, we solved a practical problem involving an area 23.7 times larger than the state-of-the-art, and conducted a comprehensive earthquake simulation by combining earthquake wave propagation analysis and evacuation analysis. Application at such scale is a groundbreaking accomplishment and is expected to change the quality of earthquake disaster estimation and contribute to society. Categories: Time-to-solution, Scalability, Peak performance I. CONTRIBUTIONS OF SUPERCOMPUTERS TO REDUCING EARTHQUAKE DISASTERS A. Overview and importance of the problem An earthquake can affect many people. The 2011 Tohoku Earthquake in Japan killed 20,000 and more than 200,000 people are still in temporary housing. The loss of lives, damage to the economy, and catastrophic damage are fresh in our memory. This damage occurred in Japan, a country that leads the world in earthquake disaster mitigation, and there are concerns over similar disasters in earthquake-prone mega-cities such as Los Angeles, San Francisco, and Tokyo. Reliable earthquake disaster estimation plays an important role in mitigating such disasters. Physics-based comprehensive earthquake simulation is the only way to make reliable estimations of such infrequent and untestable events, and is