The simplest type of dune is the transverse one, which propagates with invariant profile orthogonally to a fixed wind direction. Here we show numerically and with a linear stability analysis that transverse dunes are unstable with respect to along-axis perturbations in their profile and decay on the bedrock into barchan dunes. Any forcing modulation amplifies exponentially with growth rate determined by the dune turnover time. We estimate the distance covered by a transverse dune before fully decaying into barchans and identify the patterns produced by different types of perturbation.PACS numbers: 45.70.Qj, 05.65.+b, 45.70.Mg Wind directionality and sand availability are the main factors dictating dune morphology. Bimodal and multidirectional wind systems form longitudinal and star dunes, respectively [1]. Under unimodal winds, two types of dune may occur, depending on the amount of sand: crescent-shaped barchans, evolving on the bedrock, and transverse dunes, which appear when the ground is covered with sand [1][2][3]. Studies of dune genesis have focused on the growth of longitudinal sand-wave instabilities of a plane leading to transverse dunes -which migrate downwind with invariant profile orthogonal to the transport direction [4]. The stability of the transverse dune shape, however, has remained a long-standing open issue, of relevance for several areas of aeolian research and planetary sciences. As shown in water tank experiments, a transverse sand ridge of finite length evolving on the bedrock destabilizes and decays into barchans when subjected to a stream of nearly constant direction [5].The complete quantitative study of transverse dune evolution requires a mathematical modeling that combines the description of the average turbulent wind field with a model for sand transport in three dimensions [6][7][8][9]. Here we adapt this model in order to investigate systematically, for the first time, the stability of a transverse dune under unidirectional wind. The dune model consists of iteratively performing the calculations listed in the steps which follow. (i) Wind -the average wind shear stress field (τ ) over the terrain is calculated from the equation,where τ 0 is the wind shear stress over the flat ground, and the shear stress perturbation due to the local topography,τ , is computed by solving the three-dimensional analytical equations of Weng et al. [10]. Since this wind model is only valid for smooth surfaces, the calculation must be adapted in order to account for flow separation at the dune brink. For each longitudinal slice of the dune, a separation streamline, s(x, y), is introduced at the dune lee, where x and y are the directions longitudinal and perpendicular to the wind, respectively. The wind model is solved, then, for the envelope h s (x, y) = max{h(x, y), s(x, y)} comprising the dune surface, h(x, y), and the separation streamlines at the dune lee; these define the so-called separation bubble, inside which the wind shear is set to zero [7]. The shape of s(x, y) is approximated by a thir...