We systematically investigate the impact of the different terms of the Skyrme energy density functional of the effective nucleon-nucleon interaction, and of its associated nuclear matter (NM) properties, on the density distributions of spherical nuclei. Twenty five Skyrme force parameterizations are examined simultaneously, covering a broad range of each characteristic parameter and NM property. The diffuseness and the neutron-skin thickness are found to be the most sensitive density quantities to the force parameterization. The diffuseness is indicated to decrease with increasing the central zero-range and the effective mass terms of the effective force, and the power σ of its density dependent term, as well as with the coefficient of the NM symmetry energy (asym) and its density slope (L) at saturation density, and the incompressibility (Ko). In contrast, the proton and neutron diffuseness tend to increase with increasing the spin–orbit force and the isoscalar effective nucleon-mass (m*), and to increases slightly with the density dependence parameters other than the power σ. Opposite impacts are pointed out for the different parts of the finite-range, and J2 tensor terms on the proton and neutron density. While the neutron-skin thickness tends to increase significantly upon increasing the central zero-range and spin–orbit force terms, asym, L, and Ko, and to increase slightly with the finite-range and J2 tensor terms, and σ, it decreases with the effective-mass term, the density-dependence exchange parameter, and with the indicated isoscalar effective mass. The proton and neutron radii exhibit decreasing behavior with the central zero-range and the spin–orbit terms, and with Ko, and m*. Increasing asym and L indicate slightly less (larger) proton (neutron) radius.