$\DeclareMathOperator\GL{GL}\DeclareMathOperator\Mat{Mat}$Consider the general linear group
$$
\GL(n) = \left\lbrace
\left(\begin{array}{cc}
A & C \\ 
M & B
\end{array}\right) \text{ with } A\in \Mat(k,k),\: B\in \Mat(n-k,n-k),\: M\in \Mat(n-k,k),\: C\in \Mat(k,n-k)
\right\rbrace
$$
of $n\times n$ invertible matrices, and the following $n\times n$ matrix 
$$
J = \left(\begin{array}{cc}
I_{k} & 0 \\ 
0 & 0
\end{array}\right)
$$
where $I_k$ is the $k\times k$ identity matrix. Consider the action of $\GL(n)$ on the projective space $\mathbb{P}^N$ of $n\times n$ symmetric matrices modulo scalar given by $(P,S)\mapsto PSP^T$. Then 
$$
PJP^T = 
\left(\begin{array}{cc}
AA^T & AM^T \\ 
MA^T & MM^t
\end{array}\right)
$$
Now, consider the subgroup $F\subset G$ defined by imposing $M = 0$ and $AA^T = cI_k$ for some $c\neq 0$ in the base field (which we can assume algebraically closed and of characteristic zero). Let $X_k = G/F$ be the orbit of $I_k$. Then $X_k$ has dimension $\frac{2nk-k^2+k-2}{2}$. 

I would like to ask if anyone knows a method to compute the Picard group of $X_k$. This should be $\mathbb{Z}$ for $k = 1$ and $\mathbb{Z}/n\mathbb{Z}$ for $k = n$. Thank you.