We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10 8 -10 12 M ), dust mass (10 6 -10 9 M ), and star-formation rates (0.01-100M yr −1 ), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 ± 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065-0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.

No abstract

We present the dust mass function (DMF) of 15 750 galaxies with redshift $z$ &lt; 0.1, drawn from the overlapping area of the GAMA and H-ATLAS surveys. The DMF is derived using the density corrected Vmax method, where we estimate Vmax using: (i) the normal photometric selection limit (pVmax) and (ii) a bivariate brightness distribution (BBD) technique, which accounts for two selection effects. We fit the data with a Schechter function, and find $M^{*}=(4.65 \pm 0.18)\times 10^{7}\,h^2_{70}\, \mathrm{ M}_{\odot }$, α = (−1.22 ± 0.01), $\phi ^{*}=(6.26 \pm 0.28)\times 10^{-3}\,h^3_{70}\,\rm Mpc^{-3}\,dex^{-1}$. The resulting dust mass density parameter integrated down to 104 M⊙ is Ωd = (1.11 ± 0.02) × 10−6 which implies the mass fraction of baryons in dust is $f_{m_\mathrm{ b}}=(2.40\pm 0.04)\times 10^{-5}$; cosmic variance adds an extra  7–17 per cent uncertainty to the quoted statistical errors. Our measurements have fewer galaxies with high dust mass than predicted by semi-analytic models. This is because the models include too much dust in high stellar mass galaxies. Conversely, our measurements find more galaxies with high dust mass than predicted by hydrodynamical cosmological simulations. This is likely to be from the long time-scales for grain growth assumed in the models. We calculate DMFs split by galaxy type and find dust mass densities of Ωd = (0.88 ± 0.03) × 10−6 and Ωd = (0.060 ± 0.005) × 10−6 for late types and early types, respectively. Comparing to the equivalent galaxy stellar mass functions (GSMF) we find that the DMF for late types is well matched by the GSMF scaled by (8.07 ± 0.35) × 10−4.