During the breakup of continents extension is commonly accommodated by connected networks of fluid filled fractures (dykes) and by faults. Despite the importance of these two extension mechanisms their spatial relationship in three dimensions is poorly understood primarily because it is difficult to quantify the subsurface distribution of faulting and intrusion. In order to address this problem, we conduct a quantitative analysis of the spatial distribution and clustering of earthquakes and volcanic vents in the Main Ethiopian Rift in East Africa in order to understand how extension by faulting and intrusion is distributed throughout the volcanic rift. We use fractal analysis of earthquake epicentres in order to infer the 2D characteristics of the subsurface fault network, and directly test our model results against the 3D distribution of earthquake hypocentres. Our results show that fractal analysis of these features is a reliable means to characterise the 3D properties of the fault network. In addition, the strong similarity between the properties of the fault network derived from earthquakes and properties of the magma-filled fracture network derived from fractal analysis of volcanic vents strongly suggests that these are genetically linked. We then explore their spatial link using computation of earthquake and vent density, which shows that the zone of seismicity is generally around 20-30-km-wide, while the zone of vents is narrower and centred within the zone of seismicity. This spatial relationship suggests that the faults, which form rift axial grabens, are induced above a narrower and central zone of diking. We also demonstrate significant along-rift variation in degree of magmatism and faulting with regions of increased degree of diking inferred from a higher cone density characterised by reduced degree of faulting.