Recent oceanographic field measurements and high-resolution numerical modelling studies have revealed intense, transient, submesoscale motions characterised by a horizontal length scale of 100-10,000 m. This submesoscale activity increases in the fall and winter when the mixed layer (ML) depth is at its maximum. In this study, the submesoscale motions associated with a large-scale anticyclonic gyre in the central Gulf of Taranto were examined using realistic submesoscale-permitting simulations. We used realistic flow field initial conditions and multiple nesting techniques to perform realistic simulations, with very-high horizontal resolutions (> 200 m) in areas with submesoscale variability. Multiple downscaling was used to increase resolution in areas where instability was active enough to develop multi-scale interactions and produce 5-km-diameter eddies. To generate a submesoscale eddy, a 200-m resolution was required. The submesoscale eddy was formed through small-scale baroclinic instability in the rim of a large-scale anticyclonic gyre leading to large vertical velocities and rapid restratification of the ML in a time-scale of days. The submesoscale eddy was confirmed by observational data from the area and we can say that for the first time we have a proof that the model reproduces a realistic