We present an empirical chemical evolution model that explains the distribution of metals in the interstellar medium (ISM) and the circumgalactic medium (CGM) of galaxies based on the UniverseMachine and NeutralUniverseMachine models in the framework of $\Lambda$CDM structure formation. We parameterize the fractions of outflowing metals returned and mixed into the multi-phase ISM of the star-forming regions ($f_{\rm H2}$) and into the neutral gas regions ($f_{\rm HI}$); metal production, transfer, and dilution are caused by star formation, galaxy mergers, and gas inflow from the inter-galactic medium, respectively, with rates determined by the (Neutral)UniverseMachine models. Using a Markov Chain Monte Carlo algorithm, we explore the posterior distributions of metal return and mixing consistent with observed mass-metallicity relations in HII regions (at $0Comment: 27 pages, 17 figures, accepted for publication in ApJ