Abstract--A new approach to the interpretation of M6ssbauer spectra of Fe3+-phyllosilicates having vacant trans-octahedra is based on (1) crystal structure simulation methods that allow for the size and the shape ofa Fe3+-octahedron as a function of the nearest surrounding cations; and (2) calculations of electric field gradients (EFG) on Fe 3+ in terms of the ionic point-charge model Calculations were performed by direct summation within the region of radius -<50/~. Coordinates for the anions in the coordination octahedra have been assigned to take into account the nearest cationic environment. Atomic coordinates for the rest of the summation volume are those for the average unit cell. EFG calculations for cation combinations responsible for the visible quadrupole splitting Avi~ in the spectra of nontronite, "red" muscovite, and celadonite have led to good agreement between Avis and A~ c. Computer fitting of the nontronite and celadonite spectra based on EFG calculations for the rest of the possible cation combinations suggests that the distribution of tetrahedral cations in nontronite obeys the Loewenstein rule, and in celadonite, the distribution of R 3+ and R 2+ over cis-octahedra is predominantly ordered, in agreement with electron diffraction and infrared spectroscopy data. The M6ssbauer spectrum of one of the glauconites suggested the presence of celadonite-like and muscovite-like domains in its 2:1 layers.