It appears that the zeta map which I used to answer [Vince Vatter's initial question](http://mathoverflow.net/questions/131585) and which I describe in [my answer there](http://mathoverflow.net/questions/131585#131767), see also page 50 of [Jim Haglund's book](http://math.upenn.edu/~jhaglund/books/qtcat.pdf), indeed solves also this problem:

As discussed in a comment above, the zeta map has the property that it sends

- the number of $0$'s in an "area sequence" $a = (a_1,\ldots,a_n)$ (item $u$ in [Stanley's list](http://www-math.mit.edu/~rstan/ec/catalan.pdf), see also Property $D$ in [my answer to the other question](http://mathoverflow.net/questions/131585#131767)) to the number of initial north (or up) steps, and
- the number of $i$'s for which all $i$'s appear before all $i+1$'s within $a$ to the number of returns to $0$ (excluding the very last return).

On the other hand, David's bijection between rooted planar trees and Dyck paths sends

- the number of children of the root to the number of $0$'s in the corresponding area sequence, and 
- the number of crucial vertices in a tree (excluding the root which is always, except for $n=1$, crucial) to the number of $i$'s for which all $i$'s appear before all $i+1$'s within the corresponding area sequence.

Combining his bijection with the zeta map yields therefore a bijection between rooted planar trees and Dyck paths that sends the bistatistic

(number of children of the root, number of crucial vertices)

to the bistatistic

(number of initial north steps, number of returns to $0$).

We thus proved that $\sum c(p,q,n)x^py^q$ is indeed equal to the Tutte polynomial $T_n(x,y)$.