Two-dimensional magnetism in CrI3 links to its stacking structural order. In bulk CrI3, coexistence of monoclinic and rhombohedral phases has been observed across wide temperature range. Direct characterization of the magnetic textures on the exposed surface of bulk CrI3 has been challenging but meaningful. Here, we have employed a custom glovebox-assisted magnetic force microscope to observe the evolution of magnetic domains on the bare CrI3 surface under variable temperatures and magnetic fields for the first time. Below the Curie temperature, robust parallel stripe magnetic domains spontaneously emerge. The root mean square values of the localized weak pinning domain images at various temperatures indicate a critical transition behavior consistent with 3D Ising model with Tc ∼ 59 K. Images obtained by varying the magnetic field at low temperatures, showing magnetization saturation-reentry cycles, reveal larger hysteresis and stronger disorder compared to the VSM measurement, indicating local structural phase inhomogeneity. Two in-plane stripe magnetic domains oriented locally at an angle of 51-degrees compete under an out-of-plane magnetic field, presumably due to different sliding orientations induced by the magnetic field during the transition from a monoclinic to a rhombohedral phase. This study provides an experimental reference for future spin density wave manipulation in CrI3.