In the present work, we study the oscillation dynamics of 3 mm-diameter bubbles generated through an orifice submerged in viscous liquids. The viscosity of those liquids is varied to change the behavior of the rising bubble. The details of the rising motion and shape oscillation of the bubbles are measured using a combination of high speed, high-resolution imaging, and an accurate digital image processing technique. Direct Numerical Simulations that mimic the experimental conditions are also performed using a front-tracking technique, called the Local Front Reconstruction Method. The predictions of the bubble shape and rising velocity obtained by the numerical simulations show good agreement with the experimental results. Our experimental and numerical results show that the oscillation frequency and the damping rate at lower modes can be predicted using available theoretical models found in the literature. However, discrepancies arise between our results with the theoretical predictions at higher order oscillation modes. We conclude that the discrepancies are due to the influence of rising motion and the vortex wave, which is not considered in the theoretical models.