We present 2000 mock galaxy catalogs for the analysis of baryon acoustic oscillations in the Emission Line Galaxy ( ELG ) sample of the Extended Baryon Oscillation Spectroscopic Survey Data Release 16 ( eBOSS DR16 ) . Each mock catalog has a number density of 6.7 \times 10 ^ { -4 } h ^ { 3 } Mpc ^ { -3 } , covering a redshift range from 0.6 to 1.1 . The mocks are calibrated to small-scale eBOSS ELG clustering measurements at scales of around 10 h ^ { -1 } Mpc . The mock catalogs are generated using a combination of GaLAxy Mocks ( GLAM ) simulations and the Quick Particle-Mesh ( QPM ) method . GLAM simulations are used to generate the density field , which is then assigned dark matter halos using the QPM method . Halos are populated with galaxies using a halo occupation distribution ( HOD ) . The resulting mocks match the survey geometry and selection function of the data , and have slightly higher number density which allows room for systematic analysis . The large-scale clustering of mocks at the baryon acoustic oscillation ( BAO ) scale is consistent with data and we present the correlation matrix of the mocks .