We investigate the effect of varying the ion temperature gradient (ITG) and toroidal equilibrium scale sheared flow on ion-scale turbulence in the outer core of MAST by means of local gyrokinetic simulations. We show that nonlinear simulations reproduce the experimental ion heat flux and that the experimentally measured values of the ITG and the flow shear lie close to the turbulence threshold. We demonstrate that the system is subcritical in the presence of flow shear, i.e., the system is formally stable to small perturbations, but transitions to a turbulent state given a large enough initial perturbation. We propose that the transition to subcritical turbulence occurs via an intermediate state dominated by low number of coherent long-lived structures, close to threshold, which increase in number as the system is taken away from the threshold into the more strongly turbulent regime, until they fill the domain and a more conventional turbulence emerges. We show that the properties of turbulence are effectively functions of the distance to threshold, as quantified by the ion heat flux. We make quantitative comparisons of correlation lengths, times, and amplitudes between our simulations and experimental measurements using the MAST BES diagnostic. We find reasonable agreement of the correlation properties, most notably of the correlation time, for which significant discrepancies were found in previous numerical studies of MAST turbulence. * ferdinand.vanwyk@physics.ox.ac.uk † highcock@chalmers.se ‡ alex.schekochihin@physics.ox.ac.uk arXiv:1704.02830v2 [physics.plasm-ph] 1 Aug 2017 2 V 0 dV B · ∇φ is the toroidal magnetic flux, V is the volume enclosed by the flux surface, B is the magnetic field, φ is the toroidal angle, and ψ tor,LCFS is the toroidal flux enclosed by the last closed flux surface [see figure 1(b)], ψ pol = (1/2π) 2 V 0 dV B · ∇θ is the poloidal magnetic flux, θ is the poloidal angle, and ψ pol,LCFS is the poloidal flux enclosed by the LCFS. 4 http://w3.pppl.gov/transp/ c , including two cases that match the experimental level of heat flux. This is a considerable improvement over previous nonlinear gyrokinetic simulations of this MAST discharge [51], which overpredicted τ SYNTH c by two orders of magnitude. Examining figure 24(d), we see that (δn i /n i ) SYNTH rms increases with increasing κ T or de-