In the past decade, numerous studies have successfully mapped thousands of former charcoal production sites (also called relict charcoal hearths) manually using digital elevation model (DEM) data from various forested areas in Europe and the north‐eastern USA. The presence of these sites causes significant changes in the soil physical and chemical properties, referred to as legacy effects, due to high amounts of charcoal that remain in the soils. The overwhelming amount of charcoal hearths found in landscapes necessitates the use of automated methods to map and analyse these landforms. We present a novel approach based on open source data and software, to automatically detect relict charcoal hearths in large‐scale LiDAR datasets (visualized with Simple Local Relief Model). In addition, the approach simultaneously provides both general as well as domain‐specific information, which can be used to further study legacy effects. Different versions of the methodology were fine‐tuned on data from north‐western Connecticut and subsequently tested on two different areas in Connecticut. The results show that these perform adequate, with F1‐scores ranging between 0.21 and 0.76, although additional post‐processing was needed to deal with variations in LiDAR quality. After testing, the best performing version of the prediction model (with an average F1‐score of 0.56) was applied on the entire state of Connecticut. The results show a clear overlap with the known distribution of charcoal hearths in the state, while new concentrations were found as well. This shows the usability of the approach on large‐scale datasets, even when the terrain and LiDAR quality varies.