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^0rightarrow K_{S}^{0} K^{+} K^{-}$. The Dalitz plot is analyzed using $1856pm 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 $Kbar{K}$, $g_{Kbar{K}}=3.77pm 0.24text{(stat.)}pm0.35text{(sys.)} GeV$. The branching fraction of the decay $D^0rightarrow K_{S}^{0} K^{+} K^{-}$ is measured using $11660pm 118$ untagged signal decays to be $(4.51pm 0.05text{(stat.)}pm 0.16text{(sys.)})10^{-3}$. Both measurements are limited by their systematic uncertainties.