Based on direct numerical simulations and symmetry arguments, we show that the large-scale Fourier modes are useful tools to describe the flow structures and dynamics of flow reversals in Rayleigh-Bénard convection (RBC). We observe that during the reversals, the amplitude of one of the large-scale modes vanishes, while another mode rises sharply, very similar to the "cessation-led" reversals observed earlier in experiments and numerical simulations. We find anomalous fluctuations in the Nusselt number during the reversals. Using the structures of the RBC equations in the Fourier space, we deduce two symmetry transformations that leave the equations invariant. These symmetry transformations help us in identifying the reversing and non-reversing Fourier modes.PACS numbers: 47.55. 47.27.De, Many experiments [1][2][3][4][5][6][7] and numerical simulations [5,[8][9][10] on turbulent convection reveal that the velocity field of the system reverses randomly in time (also see review articles [11]). This phenomenon, known as "flow reversal", remains ill understood. This process gains practical importance due to its similarities with the magnetic field reversals in geodynamo and solar dynamo [12]. In this letter, we study the dynamics and symmetries of flow reversals in turbulent convection using the large-scale Fourier modes of the velocity and temperature fields.The experiments and simulations performed to explore the nature of flow reversals are typically for an idealized convective system called Rayleigh-Bénard convection (RBC) in which a fluid confined between two plates is heated from below and cooled at the top. Detailed measurements show that the first Fourier mode vanishes abruptly during some reversals [3,4]. These reversals are referred to as "cessation-led". Recently Sugiyama et al. [5] performed RBC experiments on water in a quasi two-dimensional box, and observed flow reversals with the flow profile dominated by a diagonal large-scale roll and two smaller secondary rolls at the corners. They attribute the flow reversals to the growth of the two smaller corner rolls as a result of plume detachments from the boundary layers.Several theoretical studies performed to understand reversals in RBC provide important clues. Broadly, these works involve either stochasticity (e.g., "stochastic resonance" [8, 13]), or low-dimensional models with noise [14,15]. Mishra et al.[10] studied the large-scale modes of RBC in a cylindrical geometry and showed that the dipolar mode decreases in amplitude and the quadrupolar mode increases during the cessation-led reversals. Regarding dynamo, low-dimensional models constructed using the large-scale modes and symmetry arguments reproduced dynamo reversals successfully [7,16].The theoretical models described above only focus on the large-scale modes. Here too, they provide limited information about these modes due to small number of measuring probes. In this letter we compute large-scale and intermediate-scale Fourier modes accurately using the complete flow profile. This helps us ...