At page 763 of Trace theorems for vector fields - R. Monti, D. Morbidelli the Carnot-Carathéodory metric $d$ induced on $\mathbb{R}^2$ by the vector fields $$X_1 = \partial_x \ \text{ and } \ X_2 = |x|^\alpha \partial_y, \ \alpha >0$$ is computed rather explicitely. In particular they say:
Lemma. Let $\lambda >0$. For all $(x,y), (\xi, \eta) \in \mathbb{R}^2$ with $|x|\ge|\xi|$ $$d((x,y),(\xi,\eta)) \simeq |x-\xi| + \frac{|y-\eta|}{|x|^\alpha} \quad \text{ if } \quad |x|^{\alpha +1} \ge \lambda |y-\eta|$$and $$d((x,y),(\xi,\eta)) \simeq |x-\xi| + |y-\eta|^{1/(\alpha+1)} \quad \text{ if } \quad |x|^{\alpha +1} < \lambda |y-\eta|$$where the equivalence constants depend on $\lambda$.
What I would like to know is if there exists a similar way to estimate the CC metric $\hat{d}$ induced on $\mathbb{R}^3$ by the vector fields (called Martinet field - see e.g. Control of Nonholonomic Systems: from Sub-Riemannian Geometry to Motion Planning - F. Jean, page 18) $$X_1 = \partial_x \ \text{ and } \ X_2 = \partial_y + x^2 \partial_z.$$