The claim is that if you apply gram schmidt to the columns of a matrix whose entries are i.i.d normal the resulting distribution is Haar measure on orthogonal matrices. Gram-Schmidt gives you an orthogonal matrix so invariance of the distribution under orthogonals is the issue. It is true because the orthogonals preserve i.i.d normals and also all elements of the gram-schmidt procedure, i.e. inner products. Suppose you have a $2 \times 2$ matrix $X$ , the general case being the same, with columns $x_1, x_2$. Gram schmidt gives an orthogonal matrix with columns $\frac {x_1}{\Vert x_1 \Vert},\frac {x_2 - \langle x_1,x_2 \rangle \frac {x_1}{\Vert x_1 \Vert}}{\Vert \text {the numerator }\Vert }$. Call that $X_{GS}$. If O is an orthogonal matrix $OX$ has columns $Ox_1, Ox_2$, and when you apply gram schmidt to it you get $\frac {Ox_1}{\Vert Ox_1 \Vert},\frac {Ox_2 - \langle Ox_1,Ox_2 \rangle \frac {Ox_1}{\Vert x_1 \Vert}}{\Vert O\text {the numerator }\Vert }$. Using orthogonality of $O$ once this has the same distribution as $X_{GS}$ since $Ox_i$ are i.i.d. with the same distribution as $x_1,x_2$,and that determines the distribution, and using it again the matrix is equal to $\frac {Ox_1}{\Vert x_1 \Vert},\frac {Ox_2 - \langle x_1,x_2 \rangle \frac {Ox_1}{\Vert x_1 \Vert}}{\Vert \text {the numerator }\Vert }$ since norms and inner products are the same. The last expression is $OX_{GS}$, and this shows that the distribution of $X_{GS}$ is invariant under orthogonal matrices.