  Zonal flows in rapidly-rotating celestial objects such as the Sun, gas or ice
giants form in a variety of surface patterns and amplitudes. Whereas the
differential rotation on the Sun, Jupiter and Saturn features a super-rotating
equatorial region, the ice giants, Neptune and Uranus harbour an equatorial jet
slower than the planetary rotation. Global numerical models covering the
optically thick, deep-reaching and rapidly rotating convective envelopes of gas
giants reproduce successfully the prograde jet at the equator. In such models,
convective columns shaped by the dominant Coriolis force typically exhibit a
consistent prograde tilt. Hence angular momentum is pumped away from the
rotation axis via Reynolds stresses. Those models are found to be strongly
geostrophic, hence a modulation of the zonal flow structure along the axis of
rotation, e.g. introduced by persistent latitudinal temperature gradients,
seems of minor importance. Within our study we stimulate these thermal
gradients and the resulting ageostrophic flows by applying an axisymmetric and
equatorially symmetric outer boundary heat flux anomaly ($Y_{20}$) with
variable amplitude and sign. Such a forcing pattern mimics the thermal effect
of intense solar or stellar irradiation. Our results suggest that the
ageostrophic flows are linearly amplified with the forcing amplitude $q^\star$
leading to a more pronounced dimple of the equatorial jet (alike Jupiter). The
geostrophic flow contributions, however, are suppressed for weak $q^\star$, but
inverted and re-amplified once $q^\star$ exceeds a critical value. The inverse
geostrophic differential rotation is consistently maintained by now also
inversely tilted columns and reminiscent of zonal flow profiles observed for
the ice giants. Analysis of the main force balance and parameter studies
further foster these results.
