Let $C_n=\frac1{n+1}\binom{2n}n$ denote the Catalan numbers.
This question is motivated by the (unanswered) MO post by Alexander Burstein and my own (answered by Fedor Petrov) MO post. Specifically, Shapiro's convolution formula states that $$\sum_{k=0}^nC_{2k}C_{2n-2k}=4^nC_n. \tag1$$
QUESTION. Is there a $q$-analogue to this "elegant" identity in equation (1)?