We study the problem of constructing ε-coresets for the (k, z)-clustering problem in a doubling metric M (X, d). An ε-coreset is a weighted subset S ⊆ X with weight function w : S → R ≥0 , such that for any k-subset C ∈ [X] k , it holds thatWe present an efficient algorithm that constructs an ε-coreset for the (k, z)-clustering problem in M (X, d), where the size of the coreset only depends on the parameters k, z, and the doubling dimension ddim(M ). To the best of our knowledge, this is the first efficient -coreset construction of size independent of |X| for general clustering problems in doubling metrics.To this end, we establish the first relation between the doubling dimension of M (X, d) and the shattering dimension (or VC-dimension) of the range space induced by the distance d. Such a relation is not known before, since one can easily construct instances in which neither one can be bounded by (some function of) the other. Surprisingly, we show that if we allow a small (1 ± )-distortion of the distance function d (the distorted distance is called the smoothed distance function), the shattering dimension can be upper bounded by O( −O(ddim(M )) ). For the purpose of coreset construction, the above bound does not suffice as it only works for unweighted spaces. Therefore, we introduce the notion of τ -error probabilistic shattering dimension, and prove a (drastically better) upper bound of O(ddim(M )·log(1/ε)+log log 1 τ ) for the probabilistic shattering dimension for weighted doubling metrics. As it turns out, an upper bound for the probabilistic shattering dimension is enough for constructing a small coreset. We believe the new relation between doubling and shattering dimensions is of independent interest and may find other applications.Furthermore, we study robust coresets for (k, z)-clustering with outliers in a doubling metric. We show an improved connection between α-approximation and robust coresets. This also leads to improvement upon the previous best known bound of the size of robust coreset for Euclidean space [Feldman and Langberg, STOC 11]. The new bound entails a few new results in clustering and property testing.As another application, we show constant-sized (ε, k, z)-centroid sets in doubling metrics can be constructed by extending our coreset construction. Prior to our result, constant-sized centroid sets for general clustering problems were only known for Euclidean spaces. We can apply our centroid set to accelerate the local search algorithm (studied in [Friggstad et al., FOCS 2016]) for the (k, z)-clustering problem in doubling metrics.