Coastal floodplains are commonly underlain by sulfidic sediments and coastal acid sulfate soils (CASS). Oxidation of sulfidic sediments leads to increases in acidity and mobilisation of trace metals, resulting in an increase in the concentrations of conducting ions in sediment and pore water. The distribution of these sediments on floodplains is highly heterogeneous. Accurately identifying the distribution of CASS is essential for developing targeted management strategies. One approach is the use of digital soil mapping (DSM) using ancillary information. Proximal sensing instruments such as an EM38 can provide data on the spatial distribution of soil salinity, which is associated with CASS, and can be complemented by digital elevation models (DEM). We used EM38 measurements of the apparent soil electrical conductivity (ECa) in the horizontal and vertical modes in combination with a high resolution DEM to delineate the spatial distribution of CASS. We used a fuzzy k-means algorithm to cluster the data. The fuzziness exponent, number of classes (k) and distance metric (i.e. Euclidean, Mahalanobis and diagonal) were varied to determine a set of parameters to identify CASS. The mean-squared prediction error variance of the class mean of various soil properties (e.g. EC1:5 and pH) was used to identify which of these metrics was suitable for further analysis (i.e. Mahalanobis) and also determine the optimal number of classes (i.e. k = 4). The final map is consistent with previously defined soil–landscape units generated using traditional soil profile description, classification and mapping. The DSM approach is amenable for evaluation on a larger scale and in order to refine CASS boundaries previously mapped using the traditional approach or to identify CASS areas that remain unmapped.