This question asks us to redraw the region of integration for $I$ and read off the new order $dy\,dx$. Let's identify the exact region traced by the two given pieces and check each option by matching both the region description and the numeric value of $I$.
The first piece $\int_0^1\int_0^{\sqrt y} dx\,dy$ covers all points with $0\le x\le\sqrt y$ and $0\le y\le1$; squaring the bound $x\le\sqrt y$ gives $x^2\le y$, so this piece is the region between the parabola $y=x^2$ and the line $y=1$, for $0\le x\le1$.
The second piece $\int_1^2\int_{\sqrt{y-1}}^1 dx\,dy$ covers points with $\sqrt{y-1}\le x\le1$ and $1\le y\le2$; squaring gives $y-1\le x^2$, i.e. $y\le x^2+1$, so this piece is the region between $y=1$ and $y=x^2+1$, again for $0\le x\le1$.
The two pieces meet exactly at $y=1$, so together they form one continuous region: for every $x$ from $0$ to $1$, $y$ runs from $x^2$ up to $x^2+1$. No split at $x=\tfrac12$ is needed.
Computing the original integral directly confirms the check: $\int_0^1\sqrt y\,dy=\tfrac23$ and $\int_1^2\left(1-\sqrt{y-1}\right)dy=\tfrac13$, so $I=\tfrac23+\tfrac13=1$. Evaluating option B the same way, $\int_0^1\left[(x^2+1)-x^2\right]dx=\int_0^1 1\,dx=1$, the same value.
Let's summarize:
So the correct expression after changing the order of integration is $\int_0^1\int_{x^2}^{x^2+1}dy\,dx$, option B.
\[ \boxed{I = \int_0^1 \int_{x^2}^{x^2+1} dy\, dx} \]