This thesis studies the dynamics of a rotating compressible gas sphere, driven by internal convection, as a model for the dynamics on the giant planets. We develop a new general circulation model for the Jovian atmosphere, based on the MITgcm dynamical core augmenting the nonhydrostatic model. The grid extends deep into the planet's interior allowing the model to compute the dynamics of a whole sphere of gas rather than a spherical shell (including the strong variations in gravity and the equation of state). Dierent from most previous 3D convection models, this model is anelastic rather than Boussinesq and thereby incorporates the full density variation of the planet.We show that the density gradients caused by convection drive the system away from an isentropic and therefore barotropic state as previously assumed, leading to signicant baroclinic shear. This shear is concentrated mainly in the upper levels and associated with baroclinic compressibility eects. The interior ow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy uxes, generating the observed equatorial superrotation. Heat uxes align with the axis of rotation, contributing to the observed at meridional emission. We show the transition from weak convection cases with symmetric spiraling columnar modes similar to those found in previous analytic linear theory, to more turbulent cases which exhibit similar, though less regular and solely cyclonic, convection columns which manifest on the surface in the form of waves embedded within the superrotation. We develop a mechanical understanding of this system and scaling laws by studying simpler congurations and the dependence on physical properties such as the rotation period, bottom boundary location and forcing structure.These columnar cyclonic structures propagate eastward, driven by dynamics similar to that of a Rossby wave except that the restoring planetary vorticity gradient 3 is in the opposite direction, due to the spherical geometry in the interior. We further study these interior dynamics using a simplied barotropic annulus model, which shows that the planetary vorticity radial variation causes the eddy angular momentum ux divergence, which drives the superrotating equatorial ow. In addition we study the interaction of the interior dynamics with a stable exterior weather layer, using a quasigeostrophic two layer channel model on a beta plane, where the columnar interior is therefore represented by a negative beta eect. We nd that baroclinic instability of even a weak shear can drive strong, stable multiple zonal jets. For this model we nd an analytic nonlinear solution, truncated to one growing mode, that exhibits a multiple jet meridional structure, driven by the nonlinear interaction between the eddies. Finally, given the density eld from our 3D convection model we derive the high order gravitational spectra of Jupiter, which is a measurable quantity for the upcoming JUNO mission to Jupiter. Glenn Flierl, my...