We consider a homogeneous fluid of viscosity ν confined within an oblate spheroidal cavity, x 2 /a 2 + y 2 /a 2 + z 2 /(a 2 (1 − E 2 )) = 1, with eccentricity 0 < E < 1. The spheroidal container rotates rapidly with an angular velocity Ω 0 , which is fixed in an inertial frame and defines a small Ekman number E = ν/(a 2 |Ω 0 |), and undergoes weak latitudinal libration with frequencyω|Ω 0 | and amplitude Po|Ω 0 |, where Po is the Poincaré number quantifying the strength of Poincaré force resulting from latitudinal libration. We investigate, via both asymptotic and numerical analysis, fluid motion in the spheroidal cavity driven by latitudinal libration. When |ω − 2/(2 − E 2 )| O(E 1/2 ), an asymptotic solution for E 1 and Po 1 in oblate spheroidal coordinates satisfying the no-slip boundary condition is derived for a spheroidal cavity of arbitrary eccentricity without making any prior assumptions about the spatial-temporal structure of the librating flow. In this case, the librationally driven flow is nonaxisymmetric with amplitude O(Po), and the role of the viscous boundary layer is primarily passive such that the flow satisfies the no-slip boundary condition. When |ω − 2/(2 − E 2 )| O(E 1/2 ), the librationally driven flow is also non-axisymmetric but latitudinal libration resonates with a spheroidal inertial mode that is in the form of an azimuthally travelling wave in the retrograde direction. The amplitude of the flow becomes O(Po/E 1/2 ) at E 1 and the role of the viscous boundary layer becomes active in determining the key property of the flow. An asymptotic solution for E 1 describing the librationally resonant flow is also derived for an oblate spheroidal cavity of arbitrary eccentricity. Three-dimensional direct numerical simulation in an oblate spheroidal cavity is performed to demonstrate that, in both the non-resonant and resonant cases, a satisfactory agreement is achieved between the asymptotic solution and numerical simulation at E 1.