It is the most appropriate time to characterize the Earth-like exoplanets in order to detect biosignature beyond the Earth because such exoplanets will be the prime targets of big-budget missions like JWST, Roman Space Telescope, HabEx, LUVOIR, Thirty Meter Telescope, Extremely Large Telescope, etc. We provide models for the transmission spectra of Earth-like exoplanets by incorporating the effects of multiple scattering. For this purpose we numerically solve the full multiple-scattering radiative transfer equations instead of using Beer-Bouguer-Lambert’s law, which does not include the diffuse radiation due to scattering. Our models demonstrate that the effect of this diffuse transmission radiation can be observationally significant, especially in the presence of clouds. We also calculate the reflection spectra and polarization phase curves of Earth-like exoplanets by considering both cloud-free and cloudy atmospheres. We solve the 3D vector radiative transfer equations numerically and calculate the phase curves of albedo and disk-integrated polarization by using appropriate scattering phase matrices and integrating the local Stokes vectors over the illuminated part of the disks along the line of sight. We present the effects of the globally averaged surface albedo on the reflection spectra and phase curves as the surface features of such planets are known to significantly dictate the nature of these observational quantities. Synergic observations of the spectra and phase curves will certainly prove to be useful in extracting more information and reducing the degeneracy among the estimated parameters of terrestrial exoplanets. Thus, our models will play a pivotal role in driving future observations.