(本题满分 12 分) 计算曲面积分 $$ I = \iint_ {\Sigma} 2 x ^ {3} \mathrm {d} y \mathrm {d} z + 2 y ^ {3} \mathrm {d} z \mathrm {d} x + 3 (z ^ {2} - 1) \mathrm {d} x \mathrm {d} y, $$ 其中 $\Sigma$ 是曲面 $z = 1 - x^{2} - y^{2}(z\geqslant 0)$ 的上侧