The dynamics of vortex lattices in stirred Bose-Einstein condensates have been studied at finite temperatures. The decay of the vortex lattice was observed non-destructively by monitoring the centrifugal distortions of the rotating condensate. The formation of the vortex lattice could be deduced from the increasing contrast of the vortex cores observed in ballistic expansion. In contrast to the decay, the formation of the vortex lattice is insensitive to temperature change.PACS 03.75.Fi, 67.40.Vs, 32.80.Pj Gaseous Bose-Einstein condensates (BEC) have become a testbed for many-body theory. The properties of a condensate at zero temperature are accurately described by a nonlinear Schrödinger equation. More recently, theoretical work on the ground state properties of condensates [1] has been extended to rotating condensates containing one or several vortices [2] and their dynamics. These include vortex nucleation [3,4], crystallization of the vortex lattice [5], and decay [6,7]. Experimental study has focused mainly on the nucleation of vortices, either by stirring condensates directly with a rotating anisotropy [8][9][10] or creating condensates out of a rotating thermal cloud [11]. Here we report the first quantitative investigation of vortex dynamics at finite temperature. The crystallization and decay of a vortex lattice have been studied and a striking difference is found between the two processes: while the crystallization is essentially temperature independent, the decay rate increases dramatically with temperature.The method used to generate vortices has been outlined in previous work [9,12]. Condensates of up to 75 million sodium atoms (> 80% condensate fraction) were prepared in a cigar-shaped Ioffe-Pritchard magnetic trap using evaporative cooling. The radial and axial trap frequencies of ω x = 2π × (88.8 ± 1.4) Hz, ω y = 2π × (83.3 ± 0.8) Hz, and ω z = 2π × (21.1 ± 0.5) Hz, respectively, corresponded to a radial trap asymmetry of ǫ r = (ω2)% and an aspect ratio of A = ω x /ω z = (4.20 ± 0.04). The relatively large value of the radial trap asymmetry is due to gravitational sag and the use of highly elongated pinch coils (both estimated to contribute equally). The radio frequency used for evaporation was held at its final value to keep the temperature of the condensate roughly constant throughout the experiment. The condensate's ThomasFermi radii, chemical potential, and peak density were R r = 28 µm, R z = 115 µm, 300 nK, and 4 × 10 14 cm −3 , respectively, corresponding to a healing length ξ ≃ 0.2 µm.Vortices were produced by spinning the condensate for 200 ms along its long axis with a scanning, blue-detuned laser beam (532 nm) [8,13]. For this experiment two symmetric stirring beams were used (Gaussian waist w = 5.3 µm, stirring radius 24 µm). The laser power of 0.16 mW per beam corresponded to an optical dipole potential of 310 nK. After the stirring beams were switched off, the rotating condensate was left to equilibrate in the static magnetic trap for various hold times. As in our previous work,...