Let me decompose your domain into two region, by splitting your domain using the line $y=1$:
$D_1$ is the region that bounded above by $y=x^2+1$, bounded below by $y=1$, bounded from the left by $x=-1$ and bounded on the right by $x=1$.
$$\iint_{D_1} x y \ dy dx=0$$
because for every $(x,y) \in D_1$,$(-x,y) \in D_1$.
$D_2 = D \setminus D_1$.
$D_2$ can also be described as the region that is bounded above by $y=1$, bounded below by $y=-1$, bounded on the left by $x=-1$ and bounded on the right by $y^2=x$.
For this region, for every $(x,y) \in D_2$, $(x,-y) \in D_2$. Hence
$$\iint_{D_2} x y \ dy dx=0$$
Hence
$$\iint_{D} x y \ dy dx=\iint_{D_1} x y \ dy dx+ \iint_{D_2} x y \ dy dx= 0$$