Records of past environmental conditions in shell carbonate are usually derived from compositional analysis (i.e. trace elements, stable oxygen, carbon, and nitrogen isotopes) performed along the direction of the shell's growth and thus through time. However, compositional variations within isochronous parts of the shell can distort the environmental record and are difficult to assess without extensively mapping the whole shell. Here we apply Laser Induced Breakdown Spectroscopy (LIBS) to efficiently map the elemental change throughout the growth increments of three mollusc shells (Conomurex fasciatus, Ostrea edulis, Anomalocardia flexuosa). We employ an automated LIBS setup to map the Mg/Ca composition of whole shell sections with over 2000 data points per hour. By assessing the spatial variability of Mg/Ca intensity ratios this method has the potential to mitigate distorted results while increasing the resolution of derived palaeoenvironmental information.

Elemental analysis of biogeochemical archives is an established technique used to study climate in a range of applications, including ocean circulation, glacial/interglacial climates, and anthropogenic climate change. Data from mollusc archives are especially important because of their global abundance and sub-annual resolution. Despite this potential, they are underrepresented among palaeoclimate studies, due to enigmatic physiological influences skewing the elemental record. Understanding the patterns behind these influences will improve data interpretation and lead to the development of new climate proxies. Here, we show for the first time that extensive spatial mapping of multiple mollusc specimens using Laser Induced Breakdown Spectroscopy (LIBS) across a wider region can resolve enigmatic patterns within the elemental record caused by physiological influences. 2D elemental (Mg/Ca) maps of whole limpet shells (Patella caerulea) from across the Mediterranean revealed patterns of variability within individual mollusc records as well as within isochronous parts of specimens. By registering and quantifying these patterns, we established previously uninterpretable correlations with temperature (R² > 0.8, p < 0.01). This outcome redefines the possibilities of accessing sub-annual climate proxies and presents the means to assess annual temperature ranges using oxygen isotope analysis requiring only 2 samples per shell.