Measurement of multijet production at low Bjorken x with ZEUS at HERA Inclusive dijet and trijet production in deep inelastic ep scattering was measured with ZEUS at HERA using an integrated luminossity of 82pb-1. The data cover the kinematic range of 10 < Q2 < 100 GeV2 and Bjorken x of 10-4 < xBj < 10-2. Jets were identified in the hadronic centre-of-mass (HCM) frame using the kT cluster algorithm in the longitudinally invariant inclusive mode. Measurements of dijet and trijet differential cross sections are presented as functions of Q2, xBj, jet transverse energy and jet pseudorapidity. To further examine low-xBj dynamics, multi-differential cross sections as functions of the jet correlations in transverse momenta, azimuthal angles and pseudorapidity are also presented. Calculations at O(alpha_s^3) generally describe the trijet data well and improve the description of the dijet data compared to the calculation at O(alpha_s^2).