Analytical Green's functions in the frequency domain are presented for the three-dimensional diffusion equation in an unbounded, halfspace, slab and layered media. These proposed expressions take into account the conduction and convection phenomena, assuming that the system is subjected to spatially sinusoidal harmonic heat line sources and do not require any type of discretization of the space domain. The application of time and spatial Fourier transforms along the two horizontal directions allows the solution of the three-dimensional time convection-diffusion equation for a heat point source to be obtained as a summation of one-dimensional responses. The problem is recast in the time domain by means of inverse Fourier transforms using complex frequencies in order to avoid aliasing phenomenon. Further, no restriction is placed on the source time dependence, since the static response is obtained by limiting the frequency to zero and the high frequency contribution to the response is small.The proposed functions have been verified against analytical time domain solutions, known for the case of an unbounded medium, and the Boundary Element Method solutions for the case of the half-space, slab and layered media. q