I believe that it is "well known" that the following two statistics on Dyck paths have symmetric joint distribution:
- number of returns to the axis $RET(D)$
- height of the first peak (or length of the last descent) $HFP(D)$
That is: $\sum_{D} x^{RET(D)}y^{HFP(D)} = \sum_{D} x^{HFP(D)}y^{RET(D)}$
However, I could not find a reference for that. Might it be due to Kreweras?