Su comentario parece ir muy bien encaminado. El $\nabla f \ne 0$ nos dice que las curvas de nivel son submanifolds incrustados. Podemos excluir las curvas cerradas utilizando la condición de gradiente junto con las propiedades topológicas de $\mathbb R^2$ (por ejemplo, el teorema de la curva de Jordan + el teorema del valor extremo), pero es posible que los conjuntos de niveles sean una unión de múltiples líneas; y el número de líneas/componentes conectados es no necesariamente un invariante de la función. Consideremos
$$f(x,y) = x^2 + \eta(x) y$$
donde $\eta$ es la función de bacheo
$$ \eta(x) = \begin{cases}\exp( x^2/(x^2 - 1)) & |x| < 1 \\ 0 & \text{otherwise.}\end{cases}$$
Se puede comprobar fácilmente que $f$ es suave y sin puntos críticos. Los conjuntos de niveles $f^{-1}(c)$ para $c \ge 1$ tienen tres componentes conectados mientras que para $c<1$ tienen uno:
En cierto sentido, el cambio de topología se produce en el infinito, por lo que los argumentos de la teoría de Morse que funcionan en el caso compacto no se aplican aquí: como $c \nearrow 1$ la región $f < c$ "le crecen los cuernos". Aquí hay un primer plano de lo que estoy hablando: