Even though this question already has a great answer, I think there's more to say.
Let's rephrase the question informally as "what do the constants in the ball-box theorem (and its generalizations) depend on?" After all, the original question wants constants which are uniform over a sequence, so if we can make explicit the dependence of our constants in the right way, perhaps we can show that it is uniform over an appropriate sequence. I'll answer this more vague question to the best of my ability.
As the question points out, the original proofs of ball-box type theorems rested on some kind of implicit function theorem--and this step greatly influenced the estimates. As I will discuss below, there are applications where this is not good enough. Fortunately, Tao and Wright gave an alternative approach to these results which deftly side-steps this and other problems and paves the way to stronger and more general results. From a technical standpoint, their approach opens the door to improving the dependence of constants in results which generalize the ball-box theorem. Results which are rooted in Tao and Wright's approach are not usually explicitly tied back to the original ball-box theorem (for reasons I'll discuss), and so it can be hard for someone starting out to connect all these results. I hope this post may help.
In what follows, I'll often refer to two important papers on the subject:
- [NSW] Nagel, Stein, and Wainger, Balls and metrics defined by vector fields. I. Basic
properties, Acta Math. 155 (1985), no. 1-2, 103–147. MR 793239
- [TW] Tao and Wright, $L^p$ improving bounds for averages along curves, J. Amer. Math. Soc. 16 (2003), no. 3, 605–638.
To begin with, let's discuss some of [NSW] and in what ways it is not sufficient for some applications. Then we'll turn to results rooted in [TW], which are in many ways optimal.
Let $V_1,\ldots, V_r$ be smooth vector fields on a manifold $M$. Define the sub-Riemannian ball of radius $1$, centered at $x\in M$, in terms of $V=V_1,\ldots, V_r$ by:
$$B_V(x,1):=\left\{ y\in M : \exists \gamma:[0,1]\rightarrow M, \gamma(0)=x, \gamma(1)=y, \gamma'(t)=\sum_{j=1}^r a_j(t) V_j(\gamma(t)), \left\| \sum_{j} |a_j|^2\right\|_{L^\infty}<1\right\}.$$
Set
$$B_V(x,\delta):=B_{\delta V_1,\ldots, \delta V_r}(x,1).$$
Let $W_1,\ldots, W_r$ be smooth vector fields on a manifold $M$ of dimension $n$, satisfying H\"ormander's condition of order $m$. To each $W_j$ assign the formal degree $1$. If $V$ has formal degree $e$, we assign to $[W_j, V]$ the formal degree $e+1$. Let $(X_1,d_1),\ldots, (X_q,d_q)$ be an enumeration of all such vector fields of formal degree $\leq m$. By the assumption of H\"ormander's condition, $X_1,\ldots, X_q$ span the tangent space to $M$ at every point. We let $(X,d)$ denote the list $(X_1,d_1),\ldots, (X_q,d_q)$. A key property of these vector fields, which we will come back to, is:
$$[X_i,X_j]=\sum_{d_k\leq d_i+d_j} c_{i,j}^k X_k, \quad c_{i,j}^k\in C^\infty(M),\quad (1)$$
indeed, if $d_i+d_j\leq m$, this follows from the Jacobi identity, while if $d_i+d_j>m$ this follows from the fact that $X_1,\ldots, X_q$ span the tangent space.
We set
$$B_{(X,d)}(x,\delta):=B_{\delta^{d_1}X_1,\ldots, \delta^{d_q}X_q}(x,1).$$
[NSW] showed the balls $B_{(X,d)}(x,\delta)$ and $B_W(x,\delta)$ were comparable, so let's use the former instead.
Given $x\in M$, $\delta>0$, pick $j_1=j_1(x,\delta),\ldots, j_n=j_n(x,\delta)$ so that
$$\left| \det \left( \delta^{d_{j_1}} X_{j_1}(x) |\cdots| \delta^{d_{j_n}} X_{j_n}(x)\right)\right|$$
is maximal among all such choices of $j_1,\ldots, j_n$. Set
$$\Psi_{x,\delta}(t_1,\ldots, t_n):= \exp\left( t_1 X_{j_1}+\cdots+t_n X_{j_n}\right)x.$$
The Ball-Box Theorem [NSW] (informally): $B_{(X,d)}(x,\delta)$ is ``comparable'' to the box $$\Psi_{x,\delta}( [-\delta^{d_{j_1}}, \delta^{d_{j_1}}] \times \cdots\times [-\delta^{d_{j_n}},\delta^{d_{j_n}}])$$
Instead of trying to make $B_{(X,d)}(x,\delta)$ look like a small rectangle, it is often more convienient to make it look like a ball of essentially unit radius. So let's set
$$\Phi_{x,\delta}(t_1,\ldots, t_n):= \Psi_{x,\delta}(\delta^{d_{j_1}}t_1,\ldots, \delta^{d_{j_n}} t_n).$$
Another Ball-Box-type Theorem [NSW]: There exists $\eta,\xi\approx 1$ such that
$$B_{(X,d)}(x,\xi\delta) \subseteq \Phi_{x,\delta}(B^n(\eta))\subseteq B_{(X,d)}(x,\delta). $$
Futhermore, the vector fields $\Phi_{x,\delta}^{*} \delta_{d_j} X_j$ are $C^\infty$ and span the tangent space uniformly in $\delta$. I.e., the vector fields $\Phi_{x,\delta}^{*} \delta W_1,\ldots, \Phi_{x,\delta}^{*} \delta W_r$ satisfy H\"ormander's condition uniformly in $\delta$. In short, the map $\Phi_{x,\delta}$ ``rescales'' the case of $\delta$ small to the case $\delta=1$. The implicit constants are independent of $x\in M$ and $\delta\in (0,\delta_0]$ where $\delta_0>0$ is some small number.
It will be more convienient in a moment to work on the Euclidean unit ball $B^n(1)$ instead of $B^n(\eta)$, and this can be achieved by replacing $\Phi_{x,\delta}(t)$ with $\Phi_{x,\delta}(\eta t)$. In this incarnation, it says that the ball-box theorem means the sub-Riemannian ball $B_{(X,d)}(x,\delta)$ ``looks like'' the Euclidean ball $B^n(1)$ in an appropriate coordinate system in a nice way. This is why we don't usually refer to generalizations of this result as ball-box theorems: they're more like ball-ball theorems, though that isn't very descriptive.
There are two ways that the above theorems are not quite good enough for some applications in harmonic analysis:
- Above, we discussed balls of the form $B_{\delta W_1,\ldots, \delta W_r}(x,1)$. But what if we instead wanted to study $B_{\delta_1 W_1,\ldots, \delta_r W_r}(x,1)$ and obtain results which are uniform for $\delta_1,\ldots, \delta_r$ small? What assumptions on $W_1,\ldots, W_r$ allow us to obtain uniform estimates? This so-called multi-parameter situaiton arose in some questions and is easily seen to beyond the methods of [NSW]. The methods of [NSW] are good enough if $\delta_1,\ldots, \delta_r$ are assumed to be weakly-comparable (i.e., $\delta_j^N\lesssim \kappa \delta_k$ for all $j,k$), but not if the $\delta$s are unrestricted. One could also create other situations by making the vector fields depend on $\delta$ in a complicated way, and it is desirable to study these as well.
- The conclusions of the above ball-box-type theorems are clearly invariant under arbitrary $C^2$ diffeomorphisms. But the methods of [NSW] are not: they rely on estimates of the $C^m$ norms of the coeffients of the vector filds in some fixed coordinate system.
In Section 4 of [TW], Tao and Wright sketched another approach to the main results of [NSW]; despite the fact that [TW] was working in the weakly comparable setting and they acknoledged that the methods of [NSW] were sufficient to prove what they needed. Perhaps for this reason, these methods were described without much explantion for why they were better. It took more than 10 years before it was explicitly pointed out that their methods avoided the two above issues.
To describe what is possible using their methods, a change in perspective is useful. Set
$$X^\delta_j := \delta^{d_j} X_j,$$
$$c_{i,j}^{k,\delta}:=
\begin{cases}\delta^{d_i+d_j-d_k} c_{i,j}^k & d_k\leq d_i+d_j\\
0&\text{otherwise.}
\end{cases}$$
Multiplying both sides of (1) by $\delta^{d_i+d_j}$ we have
$$[X_i^{\delta}, X_j^{\delta}] = \sum_k c_{i,j}^{k,\delta} X_k^{\delta},$$
where $c_{i,j}^{k,\delta}$ is $C^\infty$ uniformly in $\delta$.
So here is our new perspective. Let $Z_1,\ldots, Z_q$ be $C^1$ vector fields on a $C^2$ manifold $M$ of dimension $n$, which span the tangent space at every point. Suppose
$$[Z_i,Z_j]=\sum_{k} c_{i,j}^k Z_k.$$
We will make sure our results only depend on quantities like:
$$\sup_{x} \sum_{|\alpha|\leq l} |Z^{\alpha} c_{i,j}^k(x)|,\quad (2)$$
for some $l$, where $\alpha$ ranges over ordered multi-indices (since $Z_1,\ldots, Z_q$ may not commute). Thus, whatever results we prove for $Z_1,\ldots, Z_q$ will apply to $Z_1= X_1^{\delta},\ldots, Z_q=X_q^{\delta}$ uniformly in $\delta$. So we can ignore the parameter $\delta$ entirely now. Also, (2) is a quantity which is invariant under arbitrary $C^2$ diffeomorphisms, so if we only use quantities like this and other diffeomorphic invariant quantities, our results will be invariant under arbitrary $C^2$ diffeomorphisms.
Theorem: Fix a point $x_0\in M$ and take $Z_1,\ldots, Z_q$ as above. All implicit constants which follow will be as desrcibed above. There exists a map $\Phi:B^n(1)\rightarrow B_Z(x_0,1)$ which is a $C^2$ diffeomorphism onto its image which is an open subset of $B_Z(x_0,1)$. There exists $\xi\approx 1$ such that
$$B_Z(x_0,\xi)\subseteq \Phi(B^n(1))\subseteq B_{Z}(x_0,1)$$
Furthermore,
$$\| \Phi^{*} Z_j\|_{C^m}\lesssim 1, \forall m,$$
i.e., $\Phi^{*}Z_1,\ldots, \Phi^{*} Z_q$ are smooth ``uniformly''. These vector fields also span the tangent space ''uniformly''
$$\sup_{j_1,\ldots, j_n} \left| \det\left( \Phi^{*} Z_{j_1} | \cdots | \Phi^{*} Z_{j_n} \right)\right|\approx 1.$$
When applied to the vector fields $\delta^{d_1}X_1,\ldots, \delta^{d_q}X_q$, this yeilds the ball-box type theorem discussed above. But it can also be applied to the multi-parameter setting. Or in the setting of a sequence of vector fields as in the original question. If the hypotheses of the Theorem hold uniformly, then the conclusions do as well. It is also invariant under arbitary $C^2$ diffeomorphisms.
I wasn't explicit above exactly what all the constants depend on. This can be found in this paper. This might be the kind of paper the original question asked for--its main theorem can be applied to a sequence of distributions to sometimes give uniform results.
In this theorem, one can try to use the map $\Phi_{x,\delta}$ for $\Phi$. This works to a certain extent (see this paper, joint with Stovall). However, for sharper results, a different $\Phi$ is used.
It can be quite useful to have the flexibility to use a $\Phi$ other than $\Phi_{x,\delta}$. As just mentioned, it sometimes allows for sharper results. Also, one can ask the same sort of generalized "ball box theorems" in other cateogries. For example, one could want a result on complex manifolds using holomorphic maps. This is also possible--see this paper. It is better to think of these as "scaling" theorems than "ball box" theorems, though, I think.