Complex phenomena characterize the intercalation of ions inside stratified crystals. Their comprehension is crucial in view of exploiting the intercalation mechanism to change the transport properties of the crystal or obtaining a fine control of crystal delamination. In particular, the relationship between the concentration and nature of intercalated ions and surface structural modifications of the host stratified crystal is still under debate. Here, we discuss a theoretical effort to provide a rationale for some structural changes observed on the highly oriented pyrolytic graphite (HOPG) surface after electrochemical treatment in perchloric and sulphuric acid solutions. The formation of the so-called nano-protrusions on the basal plane of intercalated graphite was previously observed with scanning tunneling microscopy (STM). In this work, we employed both STM and density functional theory (DFT) simulations to elucidate the physical and chemical mechanisms driving the emergence of these nano-protrusions. The DFT results show that, in a bilayer graphene system, the presence of a single ion can generate a nano-protrusion with 2.49 Å height and 21.27 Å width. In the deformed area, the C-C bond length is stretched by about 2.5% more than the normal graphene bond. These values are of the same dimensional scale as those reported in previous STM experimental results.25 However, the simulated STM images obtained by increasing the amount of intercalated ions per area suggest that the presence of more than one ion is needed for the deformation of the uppermost graphite layer during the early stages of intercalation. In contrast, in a multilayer graphene system, no significant surface deformation is detected when ions are intercalated between the third and fourth layers. Charge analysis indicates an altered distribution of the charges as a consequence of the intercalation. The charge transfer from graphene layers to the intercalated ions results in a surface layer more prone to oxidation.