We present the characteristics of the Damped Lyman-$alpha$ (DLA) systems found in the data release DR16 of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) of the Sloan Digital Sky Survey (SDSS). DLAs were identified using the convolutional neural network (CNN) of~cite{Parks2018}. A total of 117,458 absorber candidates were found with $2 leq zdla leq 5.5$ and $19.7 leq lognhi leq 22$, including 57,136 DLA candidates with $lognhi geq 20.3$. Mock quasar spectra were used to estimate DLA detection efficiency and the purity of the resulting catalog. Restricting the quasar sample to bright forests, i.e. those with mean forest fluxes $meanflux>2timesfluxunit$, the completeness and purity are greater than 90% for DLAs with column densities in the range $20.1leq lognhi leq 22$.