Freely decaying 2D turbulence in magnetized electron columns has been measured using a phosphor screen͞CCD camera diagnostic. The relaxation towards axisymmetric equilibrium is hindered by medium-sized coherent vorticity holes which persist for hundreds of column rotations despite the strong negative rotational shear. These holes tend to settle into symmetric configurations in the column, give a negative skew to the measured vorticity fluctuations, and are eventually destroyed due to a slow outward radial creep. PACS numbers: 47.27.Jv, 47.32.Cc, 52.25.Wz, 52.35.Ra Coherent structures are important in the dynamics and transport of many flows. In two-dimensional (2D) freely decaying Navier-Stokes turbulence, computer simulations have demonstrated that self-coherent vortices emerge from random initial conditions and dominate the late-time relaxation [1]. Experiments using thin soap films and electrolyte have observed similar vortices, despite non-2D effects associated with the boundary layers and free surfaces [2,3]. Theoretical studies have analyzed coherent structures, such as drift-wave vortices, in various plasma systems [4,5]. Experimentally, the existence of such plasma structures is generally inferred from point probe measurements [6][7][8][9], and it is often presumed that these structures are destroyed by shear in the background flow. 2D quasigeostrophic models of sheared planetary atmospheres exhibit vortices [10,11], and computer simulations [11] and driven fluid experiments [12] have demonstrated the persistence of prograde vortices.In this paper, we characterize freely decaying 2D turbulence in magnetized electron columns, and find that coherent density holes (i.e., prograde vortices) survive for hundreds of column rotations. An initially hollow column exhibits instabilities, transport, and fine-scale turbulence [13], then relaxes to an axisymmetric minimum enstrophy metaequilibrium with essentially no energy loss, the "selective decay" consistent with hypothesis [14].Here, we show that this relaxation is slow because it is hindered by the longevity of elliptical density holes. The orientation and shape of these holes agree with an inviscid equilibrium model for vorticity patches in a background shear flow [15]. Hole interactions lead to a preferred configuration of two diametrically opposed holes on a distorted core (a "tripole"). On a time scale of hundreds of column bulk rotations, these holes creep radially outward and are eventually destroyed. Because of these coherent holes, the observed fluctuation relaxation rate is about 50 times slower than expected from simple passive tracer mixing, and the measured fluctuations are negatively skewed from Gaussian. Figure 1 shows the experimental device with the imaging diagnostic. Electrons from a tungsten filament are trapped in a series of conducting cylinders (radius R w 3.05 cm) enclosed in a vacuum chamber (ϳ10 29 torr). The electrons are contained axially by negative voltages (2150 V) on the two end cylinders, and confined radially by a uniform axia...