Multi-jet production has been investigated at HERA using the ZEUS detector. Three- and four-jet final states have been measured1 in photoproduction (PHP). The three- and four-jet cross sections have been compared with two leading-logarithm parton-shower Monte Carlo models, with and without multi-parton interactions. Two- and three-jet production in deep inelastic (DIS) ep scattering has been measured2 for 10 < Q2 < 100 GeV2 in the region of small Bjorken x (10-4 < xBj < 10-2) and compared to the predictions based on next-to-leading-order QCD calculations. In the third analysis3 described here, angular correlations between the jets and the proton beam direction were investigated in three-jet production in DIS to separate the contributions from the different color configurations. The measurements are compared to fixed-order perturbative calculations assuming the values of the color factors CF, CA and TF as derived from a variety of gauge groups. The photoproduction analysis used a data sample with a luminosity of 121 pb-1. Both DIS analyses used 82 pb-1. Jets were reconstructed with the kT cluster algorithm in the longitudinally invariant inclusive mode.