In April 2012, cores were retrieved from the deepest part of Lake Muzelle using a UWITECH piston corer with hammer action. The master core sequence (MUZ12 / IGSN:IEFRA00A4) is composed of two coring holes. MUZ1201 – (IGSN: IEFRA094Y / IEFRA094Z) was retrieved from the lake bottom with a 90-mm-diameter tube, and MUZ1202 was retrieved with a 63-mm tube. A short gravity core, MUZ12P5 (IGSN : IEFRA00A2), was also taken to provide a well-preserved sample of the water–sediment interface.

In June 2014, a set of nine short gravity cores (MUZ14P1 to MUZ14P9) distributed across the deep portion of the lake basin were retrieved using a UWITECH piston corer, allowing us to better understand sediment partitioning throughout the entire lake system.

All cores were split in half lengthwise and archived for re-use in other projects.

A detailed macroscopic description allowed to determine the different sedimentary facies and to establish the stratigraphic correlations between cores. A composite core was built using clearly identifiable layers in the overlapping sections of cores MUZ1201 and MUZ1202 to obtain a 283.5-cm-long sediment sequence. Geochemical analysis was carried out using an Avaatech X-ray Fluorescence Core Scanner at a resolution of 1 mm.

The sediment sequence chronology of lake Muzelle is based on ten 14C measures on terrestrial macro-remains, hort-lasting radionuclides measures and secular paleomagnetic fluctuations. the 14C measures were carried out using Accelerator Mass Spectrometry (AMS) technology.

AVAATECH XRF Core Scanner 2 / source Rhodium / Analytical settings were adjusted at 10Kv 30Kv / A62 Winaxil version 4.5.2, Winaxilbatch version 1.8.1 (Canberra) / Standards (Yes, Certified standards Names) : MESS-3, SARM-4

Measurements were performed every "1mm" resolution", using a XRF Avaatech core scanner

(EDYTEM laboratory, CNRS-Université de Savoie, France). X-Ray beam was generated with a rhodium anode and a 125 μm beryllium window which allow voltage range of 7 to 50 kV and a current range 0 to 2 mA. Each individual power spectrum was transformed by deconvolution process into relative contents (intensities) expressed in counts (cts).