If you look for a probability interpretation of QM, you get a continuity equation of the form $$\partial_{t}\rho+ \vec{\nabla}\cdot\vec{j}=0,$$ where $$\rho\left(\vec{x},t\right)\equiv\psi^*\left(\vec{x},t\right)\psi\left(\vec{x},t\right)$$ and $$\vec{j}\equiv\frac{\hbar}{2im}\left(\psi^*\vec{\nabla}\psi -
\psi\vec{\nabla}\psi^*\right).$$ Integrating the continuity equation over a finite volume $V$ and using the divergence theorem, you can find $$\frac{\mathrm{d}}{\mathrm{d}t}\int_V \mathrm{d}^3x\psi^*\psi=-\frac{\hbar}{2im}\int_\sigma \mathrm{d}\sigma\left(\psi^*\partial_n\psi-\psi\partial_n \psi^*\right),$$ where $\sigma\equiv\partial V$ is the smooth boundary surface of $V$ and $\partial_n\left(\cdot\right)$ indicates the normal derivative (directional derivative taken in the direction normal to the surface). If you let the volume $V$ tend to the whole space, $\mathbb{R}^3$, since your wavefunction is in the Hilbert space $\mathscr{L}^2\left(\mathbb{R}^3\right)$, last equation becomes $$\frac{\mathrm{d}}{\mathrm{d}t}\int_V \mathrm{d}^3x\psi^*\psi\equiv\frac{\mathrm{d}}{\mathrm{d}t}\left|\psi\left(\vec{x},t\right)\right|^2=0,$$ so you can see that the preservation in time of the norm of the wavefunction is required in order to get a probability interpretation of the theory.
Note that, at this level, you're not allowing decays, for example.