Geometrically each of the iso-lines of the dimple is an ellipse -they are congruent-, and, each colour should represent a particular height, probably steps in the gradient with equal steps in the height.
Then, in this case, there is a cut at the side, so it would only be elliptical arces instead of full ellipses.
So the problem is, articulating the scale parameter of the basic object, making each of the width of the arches at equal heights follow the Gaussian shape as an envelope.
Like those in that
old topic.
The interpolate option can only give linear interpolation as an lpe, and, as an extension the only option is to add an exponent.
Each of those can be used to represent isolines of solids of revolution, which with the linear option would produce cones, with the exponent some more curved surface, if the resulted new isolines are repositioned in equal heights.
If calculated right, this may give a close match to the Gaussian shape but I assume it would be too complicated to make two interpolation matching eachother with their exponents -if possible at all, as the imputs are limited by decimals.
And, once those iso-lines are being made, somehow each of two should be paired, to form a fill-able area,
Maybe would need a bit more tweaking to avoid the zero-width gap rendering issue.
Radial gradients can represent regular cone's isolines.
If they are viewed from the top.
Like in the attached example, with a diffuse lighting filter it may get clearer.