The spatial and temporal variability of gravity waves (GWs) potential energy (Ep) over South America (SA) was examined by analyzing temperature profiles obtained through the utilization of Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) from January 2002 to December 2021. We used the empirical orthogonal function (EOF) analysis to decompose GWs parameters and to analyze the GW variations over SA. We considered the first three eigenmodes (EOF1, EOF2, and EOF3) and their principal components (PC1, PC2, and PC3) of the EOF decomposition, which accounts for ∼80–90% of the total GWs variation over SA. Further, we analyzed the coupled variation of Ep and zonal mean wind (U) to verify their inter-dependencies using the singular value decomposition (SVD). The spatial variation showed that different localized mechanisms generate GWs at different sectors of the continent. The EOF1 of Ep comprised more than 50%, the EOF2 ∼20–25%, and the third ∼10–15% of the total GWs variation. The positive variation of GWs energy in the EOF1 is localized in the tropical region from the lower stratosphere to the lower mesosphere and southward below 1.5° S in the upper mesosphere. The spectral analysis of GWs energy showed biannual, annual, semiannual, and 11-year variations at different eigenvectors. Relative Ep (REp) showed an asymmetric hemispheric response to solar flux over South America. The REp response to QBO showed a modulating effect below 70 km and a positive response above 70 km. There is a good positive correlation between the temporal component of EOF2 of Ep and the quasi-biennial oscillation (QBO) at 30 mb and 50 mb in the PC2 temporal variation.Graphical