Singularities in $(3+1)$-dimensional Horava-Lifshitz (HL) theory of gravity are studied. These singularities can be divided into scalar, non-scalar curvature, and coordinate singularities. Because of the foliation-preserving diffeomorphisms of the theory, the number of scalars that can be constructed from the extrinsic curvature tensor $K_{ij}$, the 3-dimensional Riemann tensor and their derivatives is much large than that constructed from the 4-dimesnional Riemann tensor and its derivatives in general relativity (GR). As a result, even for the same spacetime, it may be singular in the HL theory but not in GR. Two representative families of solutions with projectability condition are studied, one is the (anti-) de Sitter Schwarzschild solutions, and the other is the Lu-Mei-Pope (LMP) solutions written in a form satisfying the projectability condition - the generalized LMP solutions. The (anti-) de Sitter Schwarzschild solutions are vacuum solutions of both HL theory and GR, while the LMP solutions with projectability condition satisfy the HL equations coupled with an anisotropic fluid with heat flow. It is found that the scalars $K$ and $K_{ij}K^{ij}$ are singular only at the center for the de Sitter Schwarzschild solution, but singular at both the center and $ r = (3M/|Λ|)^{1/3}$ for the anti-de Sitter Schwarzschild solution. The singularity at $ r = (3M/|Λ|)^{1/3}$ is absent in GR. In addition, all the generalized LMP solutions have two scalar curvature singularities, located at either $r = 0$ and $r=r_{s} > 0$, or $r=r_{1}$ and $r= r_{2}$ with $r_{2} > r_{1} > 0$, or $r=r_{s} > 0$ and $r = \infty$, depending on the choice of the free parameter $λ$.

Revtex4, six figures. Version to appear in Phys. Lett. B