We present the anisotropic clustering of emission line galaxies ( ELGs ) from the Sloan Digital Sky Survey IV ( SDSS-IV ) extended Baryon Oscillation Spectroscopic Survey ( eBOSS ) Data Release 16 ( DR16 ) . Our sample is composed of 173,736 ELGs covering an area of 1170 deg ^ { 2 } over the redshift range 0.6 \leq z \leq 1.1 . We use the Convolution Lagrangian Perturbation Theory in addition to the Gaussian Streaming Redshift-Space Distortions to model the Legendre multipoles of the anisotropic correlation function . We show that the eBOSS ELG correlation function measurement is affected by the contribution of a radial integral constraint that needs to be modelled to avoid biased results . To mitigate the effect from unknown angular systematics , we adopt a modified correlation function estimator that cancels out the angular modes from the clustering . At the effective redshift , z _ { eff } = 0.85 , including statistical and systematical uncertainties , we measure the linear growth rate of structure f \sigma _ { 8 } ( z _ { eff } ) = 0.35 \pm 0.10 , the Hubble distance D _ { H } ( z _ { eff } ) / r _ { drag } = 19.1 ^ { +1.9 } _ { -2.1 } and the comoving angular diameter distance D _ { M } ( z _ { eff } ) / r _ { drag } = 19.9 \pm 1.0 . These results are in agreement with the Fourier space analysis , leading to consensus values of : f \sigma _ { 8 } ( z _ { eff } ) = 0.315 \pm 0.095 , D _ { H } ( z _ { eff } ) / r _ { drag } = 19.6 ^ { +2.2 } _ { -2.1 } and D _ { M } ( z _ { eff } ) / r _ { drag } = 19.5 \pm 1.0 , consistent with \Lambda CDM model predictions with Planck parameters .