A boundary layer analysis was carried out to investigate the coupled phenomena of heat and mass transfer by natural convection from concentrated heat and mass sources embedded in saturated porous media. Both line and point source problems were treated. The boundary layer equations based on Darcy's law and Boussinesq approximation were solved by means of similarity transformation to obtain the details of velocity, temperature and concentration distributions above a concentrated heat source. Two important parameters, namely the Lewis number Le and the buoyancy ratio N were identified to conduct a series of numerical integrations. For the case of small Le, a substance diffuses further away from the plume centerline, such that the mass transfer influences both velocity and temperature profiles over a wide range. For large Le, on the other hand, the substance diffuses within a narrow range along the centerline. Naturally, the influence of mass transfer is limited to the level of the centerline velocity, so that a peaky velocity profile appears for positive N whereas a velocity defect emerges along the centerline for negative N. For such cases of large Let the temperature profiles are found to be fairly insensitive to Le.