This paper presents an experimental and theoretical study of the distribution of carbon atoms in the octahedral interstitial sites of the face-centered cubic (fcc) phase of the iron-carbon system. The experimental part of the work consists of Mossbauer measurements in Fe-C alloys with up to about 12 atomic percent C, which are interpreted in terms of two alternative models for the distribution of C atoms in the interstitial sites. The theoretical part combines an analysis of the chemical potential of C based on the quasichemical approximation to the statistical mechanics of interstitial solutions, with three-dimensional Monte Carlo simulations. The latter were performed by assuming a gas like mixture of C atoms and vacancies (Va) in the octahedral interstitial sites. The number of C-C, C-Va and Va-Va pairs calculated using Monte Carlo simulations are compared with those given by the quasichemical model. Furthermore, the relative fraction of the various Fe environments were calculated and compared with those extracted from the Mossbauer spectra. The simulations reproduce remarkably well the relative fractions obtained assuming the Fe(8)C(1-y) model for Mossbauer spectra, which includes some blocking of the nearest neighbour interstitial sites by a C atom. With the new experimental and theoretical information obtained in the present study, a critical discussion is reported of the extent to which such blocking effect is accounted for in current thermodynamic models of the Fe-C fcc phase. Abstract PACS Codes: 2.70.Uu, 76.