We investigate melt generation in a slowly extending lithosphere with the aim of understanding the spatial and temporal relationships between magmatism and preexisting rift systems. We present numerical models that consider feedback between melt generation and lithospheric deformation, and we incorporate three different damage mechanisms: brittle damage, creep damage, and melt damage. Melt conditions are calculated with a Helmholtz free energy minimization method, and the energy equation is solved self-consistently for latent heat and shear heating effects. Using a case of a slowly extending (1-1.5 mm/yr) continental lithosphere with a relatively low surface heat flow (~50 mW/m 2 ), we show that melt-rich shear bands are nucleated at the bottom of the lithosphere as a result of shear heating and damage mechanisms. Upon further deformation, melt zones intersect creep damage zones, thus forming channels that may be used for the melt to migrate upward. If a preexisting structure resides only in the brittle crust, it does not control the path of melt migration to the surface, and melt-filled channels propagate from the bottom upwards, independently of upper crustal structures. In contrast, a preexisting weak structure that reaches a critical depth of 20 km allows fast (~2 Ma) propagation of melt-filled channels that link melt damage from the bottom of the lithosphere to near-surface structures. Our model results may explain the short time scale, volume, and magma extraction from the asthenosphere through a low surface heat flow lithosphere, such as observed, for example, in the Harrat Ash-Shaam volcanic field (northwestern Arabia), which developed in the Arabian Plate and is spatially linked to the Azraq-Sirhan Graben.