By Givens' rotation $R(1,2;\theta)$ I mean a matrix which has the $$\begin{pmatrix} \hphantom{-}\cos \theta &\sin \theta \cr -\sin \theta & \cos \theta \end{pmatrix}$$ $2 \times 2$ block at the first and second rows and columns and a diagonal matrix for the remaining dimension. Similarly define $R(i,j;\theta)$ to be the matrix which is a rotation along the $i \wedge j$ 2-plane, $i < j$. Define the following map $$ \phi_{I,J}: T^t \to SO(n); \phi(\theta_1,\ldots, \theta_t) = R(i_1,j_1;\theta_1) \ldots R(i_t,j_t;\theta_t)$$.

My question is given a sequence of pairs $(I,J)=(i_1,j_1),(i_2,j_2),\ldots, (i_t,j_t)$, is there a way to determine whether every element in $SO(n)$ can be expressed as a product of the form $$ \phi_{I,J}(\theta_1,\ldots,\theta_t)$$ for some $\theta_1, \ldots, \theta_t$?

In fact one sufficient condition for a product of such sequence to contain an open set is obtained simply by looking at a neighborhood of the identity and taking derivatives. But does the local surjectivity imply global surjectivity? If not, then how big of a neighborhood of the identity does the image of $\phi$ cover?

For people who are comfortable with probability: imagine the $\theta_j$'s to be iid uniform $[0,2\pi)$ random variables. If the image of $\phi_{I,J}$ does contain an open set of $SO(n)$, what does the singular set of the density of the resulting measure on $SO(n)$ look like?

Answer to any of the above question will be considered a correct answer to this entry. Helps are deeply appreciated!

edit note: thanks to Robert Bryant for clarifying about the local condition.