Time and frequency domain channel models are proposed for nanonetworks utilizing the terahertz band (0.1–10 THz) for wireless communication. Nanonetworks are formed by tiny nanodevices which consist of nanoscale (molecular scale) components. Channel models capturing the unique peculiarities of the THz band are needed for designing proper physical layer techniques and for accurate performance analysis. Existing channel models have included the free space path loss and the molecular absorption loss, which is significant in the THz band. This paper theoretically analyzes scattering including multiple scattering referring to a sequence of scattering events from small particles, such as aerosols. Both the frequency and the impulse responses are derived. It is shown that the small particle scattering can result into significant additional loss that needs to be taken into account with the loss depending on the density and size distribution of the particles. It is shown that multiple scattering leads to a long tail in the impulse response. As most of the physical layer proposals for nanonetworks are based on the on-off keying, the channel response to pulse waveforms is specifically considered.