I've thought of a way to turn it into a purely topological problem.
We make a definition that, given a differentiable manifold ##M##, two differentiable functions ##f,g:M\to\mathbb R## are topographically equivalent if there exists a map ##\phi:M\to M## such that ##\phi## is a bijective diffeomorphism and ##f\circ \phi = g##.
The OP asks what branches of mathematics deals with the question of whether two scalar fields ##f,g## on a sphere are topographically equivalent.
For a differentiable map ##f:S^2\to R## define the ##f##-contour topology ##T_f## on ##S^2## to be the topology generated by the sub-base ##B_f## consisting of all connected components of pre-images ##f^{-1}((a,b))## of open intervals ##(a,b)## in ##\mathbb R##.
That sub-base may seem intimidatingly big, given ##a,b## can be any real numbers. Fortunately, we can show that the topology it generates is the same as that generated by a much smaller sub-base ##B'_f##, in which the intervals ##(a,b)## are required to be of the form ##(k\cdot 2^{-n}, (k+1)\cdot 2^{-n})## for some integers ##k,n##. Then for any interval ##(a',b')## with arbitrary real ##a',b'##, we define ##C## to be the set of pairs of integers ##(k,n)## such that ##a'\le k\cdot 2^{-n} < (k+1)\cdot 2^{-n}\le b'## and we have:
$$(a',b') = \bigcup_{(k,n)\in C} (k\cdot 2^{-n}, (k+1)\cdot 2^{-n}) $$
Hence:
$$ f^{-1}((a',b')) =
f^{-1}\left( \bigcup_{(k,n)\in C} (k\cdot 2^{-n}, (k+1)\cdot 2^{-n}) \mathbb R \right)
= \bigcup_{(k,n)\in C} f^{-1}\left((k\cdot 2^{-n}, (k+1)\cdot 2^{-n}) \mathbb R \right)
$$
so any element of ##B_f## is a union of elements of ##B'_f##, so the topology generated by ##B_f## must be contained in that generated by ##B'_f##. Since ##B'_f\subset B_f## (because ##f## is continuous), the reverse containment also holds, so the two topologies must be identical.
The elements of the sub-base ##B'## are like connected components of isolines (aka contours, aka level sets), eg a closed contour ring of height 250m around a peak of height 280m, except that they are \textit{bands} rather than lines. For instance, a band ringing the peak that covers all ground that has height between 240m and 260m. These bands can be made as narrow as we like, and thereby approximate a contour line with arbitrarily-close accuracy.
I hypothesise that functions ##f,g## are topographically equivalent iff there exists a map ##\phi:S^2\to S^2## that is a homeomorphism when its domain and range have the topologies ##T_f## and ##T_g## respectively, and that ##\phi## will also be a homeomorphism when its domain and range have the metric topologies, so that it is the homeomorphism required under the definition of `topographically equivalent'.
Some constraints may be needed to rule out cases where ##\phi## is a homeomorphism in the contour topology but not a diffeomorphism in the metric topology of ##S^2##. My instinct is that that would only fail in pathological cases. If so, the constraints would be minimal and intuitive, such as requiring that all contour lines have finite length.
If the hypothesis is correct, the question can be translated from one of topography to one of topology, with perhaps a little bit of calculus on manifolds thrown in (for the diffeomorphism).