Abstract. Efficient and robust landslide mapping and volume estimation is essential to rapidly infer landslide spatial distribution, to quantify the role of triggering events on landscape changes and to assess direct and secondary landslide-related geomorphic hazards. Many efforts have been made during the last decades to develop landslide areal mapping methods, based on 2D satellite or aerial images, and to constrain empirical volume-area (V-A) allowing in turn to offer indirect estimates of landslide volume. Despite these efforts, some major issues remain including the uncertainty of the V-A scaling, landslide amalgamation and the under-detection of reactivated landslides. To address these issues, we propose a new semi-automatic 3D point cloud differencing method to detect geomorphic changes, obtain robust landslide inventories and directly measure the volume and geometric properties of landslides. This method is based on the M3C2 algorithm and was applied to a multi-temporal airborne LiDAR dataset of the Kaikoura region, New Zealand, following the Mw 7.8 earthquake of 14 November 2016. We demonstrate that 3D point cloud differencing offers a greater sensitivity to detect small changes than a classical difference of DEMs (digital elevation models). In a small 5 km2 area, prone to landslide reactivation and amalgamation, where a previous study identified 27 landslides, our method is able to detect 1431 landslide sources and 853 deposits with a total volume of 908,055 ± 215,640 m3 and 1,008,626 ± 172,745 m3, respectively. This high number of landslides is set by the ability of our method to detect subtle changes and therefore small landslides with a carefully constrained lower limit of 20 m2 (90 % with A