In AD 1362, a major storm surge drowned wide areas of cultivated medieval marshland along the north-western coast of Germany and turned them into tidal flats. This study presents a new methodological approach for the reconstruction of changing coastal landscapes developed from a study site in the Wadden Sea of North Frisia.Initially, we deciphered long-term as well as event-related short-term geomorphological changes, using a geoscientific standard approach of vibracoring, analyses of sedimentary, geochemical and microfaunal palaeoenvironmental parameters and radiocarbon dating. In a next step, Direct Push (DP)-based Cone Penetration Testing (CPT) and the Hydraulic Profiling Tool (HPT) were applied at vibracore locations to obtain in situ high-resolution stratigraphic data. In a last step, multivariate linear discriminant analysis (LDA) was successfully applied to efficiently identify different sedimentary facies (e.g., fossil marsh or tidal flat deposits) from the CPT and HPT test dataset, to map the facies' lateral distribution, also in comparison to reflection seismic measurements and test their potential to interpolate the borehole and CPT/HPT data. The training dataset acquired for the key site from coring and DP sensing finally allows an automated facies classification of CPT/HPT data obtained elsewhere within the study area. The new methodological approach allowed a detailed reconstruction of the local coastal landscape development in the interplay of natural marsh formation, medieval land reclamation and storm surge-related land losses.