Mare basalts collected at the Apollo 15 landing site can be classified into two groups. Based on differing whole-rock major element chemistry, these groups are the quartz-normative basalt suite and the olivine-normative basalt suite. In this study we use modelling of Fe-Mg interdiffusion in zoned olivine crystals to investigate the magmatic environments in which the zonation was formed, be that within the lunar crust or during cooling within a surficial lava flow, helping to understand the thermal histories of the two basalt suites. Interdiffusion of Fe-Mg in olivine was modelled in 29 crystals in total, from six olivine-normative basalt thin sections and from three quartz-normative basalt thin sections. We used a dynamic diffusion model that includes terms for both crystal growth and intracrystalline diffusion during magma cooling. Calculated diffusion timescales range from 5 to 24 days for quartz-normative samples, and 6 to 91 days for olivine-normative samples. Similarities in diffusion timescales point to both suites experiencing similar thermal histories and eruptive processes. The diffusion timescales are short (between 5 and 91 days), and compositional zonation is dominated by crystal growth, which indicates that the diffusion most likely took place during cooling and solidification within lava flows at the lunar surface. We used a simple conductive cooling model to link our calculated diffusion timescales with possible lava flow thicknesses, and from this we estimate that Apollo 15 lava flows are a minimum of 3–6 m thick. This calculation is consistent with flow thickness estimates from photographs of lava flows exposed in the walls of Hadley Rille at the Apollo 15 landing site. Our study demonstrates that diffusion modelling is a valuable method of obtaining information about lunar magmatic environments recorded by individual crystals within mare basalt samples.