We present the first continental-scale seismic model of the lithosphere and underlying mantle beneath Southeast Asia obtained from adjoint waveform tomography (often referred to as full-waveform inversion or FWI), using seismic data filtered at periods from 20 to 150 s. Based on >3,000 hr of analyzed waveform data gathered from ∼13,000 unique source-receiver pairs, we image isotropic P-wave velocity, radially anisotropic S-wave velocity and density via an iterative non-linear inversion that begins from a 1-D reference model. At each iteration, the full 3-D wavefield is determined through an anelastic Earth, accommodating effects of topography, bathymetry and ocean load. Our data selection aims to maximize sensitivity to deep structure by accounting for body wave arrivals separately. SASSY21, our final model after 87 iterations across seven period bands, is able to explain true-amplitude data from events and receivers not included in the inversion. The trade-off between inversion parameters is estimated through an analysis of the Hessian-vector product. SASSY21 reveals detailed anomalies down to the mantle transition zone, including multiple subduction zones. The most prominent feature is the (Indo-)Australian plate descending beneath Indonesia, which is imaged as one continuous slab along the 180° curvature of the Banda Arc. The tomography confirms the existence of a hole in the slab beneath Mount Tambora and locates a high S-wave velocity zone beneath northern Borneo that may be associated with subduction termination in the mid-late Miocene. A previously undiscovered feature beneath the east coast of Borneo is also revealed, which may be a signature of post-subduction processes, delamination or underthrusting from the formation of Sulawesi.Plain Language Summary Southeast Asia is one of the world's most tectonically active regions, as evidenced by frequent large earthquakes and volcanic eruptions. We present a large-scale 3-D seismic structural model of this region down to a depth of 800 km that reveals a variety of primary features, including beneath the poorly understood islands of Borneo and Sulawesi. This is possible thanks to the use of a sizable data set of earthquakes recorded by a large number of permanent and temporary stations located in Southeast Asia, and advanced imaging methodology that is better able to capture the true physics of seismic wave propagation compared to more traditional methods. Our new model is capable of resolving variations in seismic properties associated with ongoing subduction (when one tectonic plate descends into the mantle below another plate), particularly along the northern margin of the Australian plate beneath the Sunda Arc. More subtle anomalies associated with remnant subduction, which correspond to plate fragments that remain once subduction stops, can also be imaged. These results are important for achieving a better understanding of the subduction cycle, which plays a central role in plate tectonics, and has important implications for, among other things, the evolut...