It depends very much on the symmetry of the material. In general, there is always a sort of mixture contributions either from surface dipoles, bulk dipoles (non centrosymmetric), bulk quadrupoles, spatial dispersion, and even magnetic contributions for magnetic materials. The separation of these effects e.g. determining how much comes from the surface and bulk can be challenging and requires careful analysis of the symmetry e.g. group theory by looking at the tensor components and nonlinear modelling especially using the so called SBHM. I recommend the paper:
where it is shown that SHG in Si111 is dominated by interface dipole and also bulkquadrupole with relatively the same strength if you measure SHG in reflection. Meanwhile, for SHG from zincblende which are non centrosymmetric such as GaAs the main contribution is clearly bulk dipoles as is evident from SBHM measurements which must include Fresnel coefficients. This is along with group theory because the second order nonlinear susceptibility from dipolar contributions does not vanish in such system contrary to Si. I recommend you to read more about this in: