We introduce the Phoenix Simulations, a suite of highly resolved cosmological simulations featuring hydrodynamics, primordial gas chemistry, Population III and II star formation and feedback, UV radiative transfer, and saved outputs with ∆t=200 kyr. The suite samples 73,523 distinct primordial star formation events within 3,313 distinct regions, forming 4,546 second-generation enriched star clusters by z ≥ 12 within a cumulative 156.25 Mpc 3 volume. The regions that lead to enriched star formation contain up to 167 primordial stars, with 78.7 % of regions having experienced multiple types of primordial supernovae. The extent of a primordial region, measured by its metal-rich surrounding cloud, is highly variable: the average region has radius ∼ 3 kpc, with 95 % confidence limit on the distribution of measured radii is ∼ 5−7 kpc. For continuing star formation, we find that the metallicity distribution of second generation stars is similar to that of subsequent Population II star formation, with both distributions spanning hyper metal-deficient ([Z/H]∼ −7) to super-solar ([Z/H]∼ 0.8). We find that the metallicity of second generation stars has no strong dependence on the configuration of progenitor supernovae, with the mean metallicity of second-generation stars having −1.73 <[Z/H]< −2.15. Finally, we create an interpretable regression model to predict the radius of metal-rich influence of Population III star systems within the first 7-18 Myr after the first light. The model predicts the radius with R 2 0.4 and mean squared error ≤ 0.06. The probability distribution function of predicted radii compares well to that of observed radii with Jensen-Shannon distance 0.2 for all modelled times.