Spectral and integrated solar reflectance of nonhomogeneous snowpacks are derived assuming surface reflection of direct radiation and subsurface multiple scattering. For surface reflection, a bidirectional reflectance distribution function derived for an isotropic Gaussian faceted surface is considered, and for subsurface multiple scattering an approximate solution of the radiative transfer equation is studied. Solar radiation incident on the snowpack is decomposed into direct and atmospherically scattered radiation. Spectral attenuation coefficients of ozone, carbon dioxide, water vapor, aerosol, and molecular scattering are included in the calculation of incident solar radiation. Illustrative numerical results are given for a case of North American winter atmospheric condition. The calculated dependence of spectrally integrated directional reflectance (or albedo) on solar elevation is in qualitative agreement with available observations.