The interaction of Cu with dissolved organic matter (DOM, extracted from an organic forest floor) was investigated and the resulting data was evaluated in terms of their uncertainty. The speciation of Cu over 'free' Cu (as analysed by diffusive gradients in thin films (DGT)), dissolved Cu-DOM complexes and precipitated Cu-DOM was determined as a function of pH (3.5, 4.0 and 4.5) and Cu/C ratio. The dissolved organically bound fraction was highest at pH 4.5, but this fraction decreased with increasing Cu/C ratio, which was observed for all pH levels. In the range of Cu/C = 7 x 10(-5) -2.3 x 10(-2) (mol/mol) the precipitated fraction was very small. The speciation of both Al and Fe was not affected by increasing Cu concentrations. From a continuous distribution model using the Scatchard approach, we calculated the optimal fit and corresponding upper and lower 95% uncertainty bounds of the overall stability constants (K(o)) with the shuffled complex evolution Metropolis (SCEM) algorithm. Although the optimal equation fitted the data very well, the uncertainty of the, according to literature, most reliable approach to establish stability constants, was still large. Accordingly, the usually reported intrinsic stability constants exhibited large uncertainty ranging from log K(i) = 6.0-7.1 (optimal 6.7) for pH 3.5, log K(i) = 6.5-7.1 (optimal 6.8) for pH 4.0, and log K(i) = 6.4-7.2 (optimal 6.8) for pH 4.5 and showed only little effect of pH.