I am self-studying some econometrics/linear regression model. In the text (Greene 2018), it is assumed that $\text{E}[\mathbf{\varepsilon}\mathbf{\varepsilon}'\mid \mathbf{X}]=\sigma^2\mathbf{I}$. Then the book states that, by using the variance decomposition formula, we find $$ \text{Var}[\mathbf{\varepsilon}] = \text{E}[\text{Var}[\mathbf{\varepsilon}\mid\mathbf{X}]] + \text{Var}[\text{E}[\mathbf{\varepsilon}\mid\mathbf{X}]] = \sigma^2\mathbf{I}. $$
Here, $\mathbf{\varepsilon}$ is an $n\times1$ column vector of disturbances and $\mathbf{X}$ is the $n\times K$ data matrix.
I cannot understand how equality is derived. Could someone please explain it for me? Thanks a lot in advance.