Skyrmions are topological spin textures of interest for fundamental science and applications. Previous theoretical studies have focused on systems with broken bulk inversion symmetry, where skyrmions are stabilized by easy-axis anisotropy. We investigate here systems that break surfaceinversion symmetry, in addition to possible broken bulk inversion. This leads to two distinct Dzyaloshinskii-Moriya (DM) terms with strengths D ⊥ , arising from Rashba spin-orbit coupling (SOC), and D from Dresselhaus SOC. We show that skyrmions become progressively more stable with increasing D ⊥ /D , extending into the regime of easy-plane anisotropy. We find that the spin texture and topological charge density of skyrmions develops nontrivial spatial structure, with quantized topological charge in a unit cell given by a Chern number. Our results give a design principle for tuning Rashba SOC and magnetic anisotropy to stabilize skyrmions in thin films, surfaces, interfaces and bulk magnetic materials that break mirror symmetry.Recently there has been a surge of interest in skyrmions in chiral magnetic materials 1-3 , ranging from fundamental science to potential device applications. A skyrmion is a spin texture characterized by a topological invariant that, in metallic magnets, gives rise to the topological Hall effect 4,5 and may also have implications for nonFermi liquid behavior 6 . The ability to write and erase individual skyrmions 7 , along with their topological stability, small size, and low depinning current density 8 , paves the way for potential information storage and processing applications.Experiments have focussed primarily on skyrmions in non-centrosymmetric crystals with broken bulk inversion symmetry, e.g., metals like MnSi, FeGe and insulators like Cu 2 OSeO 3 . In these bulk materials, the skyrmion crystal (SkX) phase is stable only in a very limited region of the magnetic field (H), temperature (T ) phase diagram 9-13 . On the other hand, the skyrmion phase is found to be stable over a much wider region of (T, H) in thin films of the same materials 10,14-16 , even extending down to T = 0 in some cases 14,16 . (A class of two dimensional (2D) systems shows atomic-scale skyrmions 17 arising from competing local interactions, distinct from the spin-orbit induced chiral interactions that we focus on here.)A key question that we address is this paper is: How can we enhance the domain of stability of skyrmion spin textures? We are motivated in part by the thin film experiments, and also by the possibility of chiral magnetism in new 2D systems like oxide interfaces [18][19][20] . We show how the SkX become progressively more stable over ever larger regions in parameter space of field H and magnetic anisotropy A, as the effects of broken surface inversion dominate over those of broken bulk inversion. The key parameter responsible for this behavior is the ratio D ⊥ /D of the strength of the chiral magnetic interaction arising from broken bulk inversion (D ) to that arising from broken surface inversion (D ⊥ ); see Fig...