Deep MIR surveys have revealed numerous strongly star-forming galaxies at redshift z~<2. Their MIR fluxes are produced by a combination of continuum and PAH emission features. The PAH features can dominate the total MIR flux, but are difficult to measure without spectroscopy. We aim to study star-forming galaxies by using a blind spectroscopic survey at MIR wavelengths to understand evolution of their star formation rate (SFR) and specific SFR up to z~=0.5, by paying particular attention to their PAH properties. We conducted a low-resolution (R~=50) slitless spectroscopic survey at 5-13um of 9um flux-selected sources (>0.3 mJy) around the North Ecliptic Pole with the Infrared Camera (IRC) onboard AKARI. After removing 11 AGN candidates by using the IRC photometry, we identified 48 PAH galaxies with PAH 6.2, 7.7, and 8.6um features at z<0.5. The rest-frame optical-MIR spectral energy distributions (SEDs) based on CFHT and AKARI/IRC imaging covering 0.37-18um were produced, and analysed in conjunction with the PAH spectroscopy. We defined the PAH enhancement by using the luminosity ratio of the 7.7um PAH feature over the 3.5um stellar component of the SEDs. The rest-frame SEDs of all PAH galaxies have a universal shape with stellar and 7.7um bumps, except that the PAH enhancement significantly varies as a function of the PAH luminosities. We identified a PAH-enhanced population at z~>0.35, whose SEDs and luminosities are typical of luminous infrared galaxies. They show particularly larger PAH enhancement at high luminosity, implying that they are vigorous star-forming galaxies with elevated specific SFR. Our composite starburst model that combines a very young and optically very thick starburst with a very old population can successfully reproduce most of their SED characteristics, although we could not confirm this optically thick component from our spectral analysis.

$\mathit{Herschel}$ extragalactic surveys offer a unique opportunity to efficiently select a significant number of rare and massive dusty objects, and thus gain insight into the prodigious star-forming activity that takes place in the very distant Universe. To search for $z\geq4$ dusty star-forming galaxies, in this work we consider red SPIRE objects with fluxes rising from 250 $\mu$m to $500\:\mu$m (so-called "500 $\mu$m-risers"). We aim to implement a novel method to obtain a statistical sample of "500 $\mu$m-risers" and fully evaluate our selection inspecting different models of galaxy evolution. We consider one of the largest and deepest ${\it Herschel}$ surveys, the Herschel Virgo Cluster Survey. We develop a novel selection algorithm which links the source extraction and spectral energy distribution fitting. We select 133 "500 $\mu$m-risers" over 55 deg$^{2}$, imposing the criteria: $S_{500}>S_{350}>S_{250}$, $S_{250}>13.2$ mJy and $S_{500}>$30 mJy. Differential number counts are in a fairly good agreement with models, displaying better match than other existing samples. In order to interpret the statistical properties of selected sources, which has been proven as a very challenging task due the complexity of observed artefacts, we make end-to-end simulations including physical clustering and lensing. The estimated fraction of strongly lensed sources is $24^{+6}_{-5}\%$ based on models. We present the faintest known statistical sample of "500 $\mu$m-risers" and show that noise and strong lensing have crucial impact on measured counts and redshift distribution of selected sources. We estimate the flux-corrected star formation rate density at $4<z<5$ with the "500 $\mu$m-risers" and found it close to the total value measured in far-infrared. It indicates that colour selection is not a limiting effect to search for the most massive, dusty $z>4$ sources.