  Using a data sample of $2.93~fb^{-1}$ of $e^+e^-$ collisions collected at
$\sqrt{s}=3.773 GeV$ in the BESIII experiment, we perform an analysis of the
decay $D^0\rightarrow K_{S}^{0} K^{+} K^{-}$. The Dalitz plot is analyzed using
$1856\pm 45$ flavor-tagged signal decays. We find that the Dalitz plot is well
described by a set of six resonances: $a_0(980)^0$, $a_0(980)^+$, $\phi(1020)$,
$a_2(1320)^+$, $a_2(1320)^-$ and $a_0(1450)^-$. Their magnitudes, phases and
fit fractions are determined as well as the coupling of $a_0(980)$ to
$K\bar{K}$, $g_{K\bar{K}}=3.77\pm 0.24\text{(stat.)}\pm0.35\text{(sys.)} GeV$.
  The branching fraction of the decay $D^0\rightarrow K_{S}^{0} K^{+} K^{-}$ is
measured using $11660\pm 118$ untagged signal decays to be $(4.51\pm
0.05\text{(stat.)}\pm 0.16\text{(sys.)})10^{-3}$. Both measurements are limited
by their systematic uncertainties.
