The adsorption and desorption of Kr on graphite at temperatures in the range 60 K to 88 K, was systematically investigated using a combination of several simulation techniques including: grand canonical Monte Carlo (GCMC), canonical kinetic-Monte Carlo (C-kMC) and the Mid-Density scheme (MDS). Particular emphasis was placed on the gas-solid, gas-liquid and liquid-solid 2D phase transitions. For temperatures below the bulk triple point, the transition from a 2D-liquid-like monolayer to a 2D-solid-like state is manifested as a sub-step in the isotherm. A further increase in the chemical potential leads to another rearrangement of the 2D-solidlike state from a disordered structure to an ordered structure that is signalled by (1) another sub-step in the monolayer region and (2) a spike in the plot of the isosteric heat versus density at loadings close to the dense monolayer coverage concentration. Whenever a 2D transition occurs in a grand canonical isotherm it is always associated with a hysteresis, a feature that is not widely recognized in the literature. We studied in details this hysteresis with the analysis of the canonical isotherm, obtained with C-kMC, which exhibits a van der Waals (vdW) type loop with a vertical segment in the middle. We complemented the hysteresis loop and the vdW curve with the analysis of the equilibrium transition obtained with the MDS, and found that the equilibrium transition coincides exactly with the vertical segment of the C-kMC isotherm, indicating the coexistence of two phases at equilibrium. We also analysed adsorption at higher layers and found that the 2D-coexistence is also observed, provided that the temperature is well below the triple point. Finally the 2D-critical temperatures were obtained for the first three layers and they are in good agreement with the experimental data in the literature.