For soil profile results, data trends were visually evaluated for each toxin (AEG, BMAA, and DAB) and the depth of the soil horizon. Each of the outcomes was log-transformed to statistically prove that the data followed a normal distribution. A series of linear bivariate and multivariate regressions were conducted to evaluate the crude and the adjusted associations between the different outcomes and soil horizon depth, controlling for the presence of the remaining soil toxins in the latter. P-values were reported and tested for each association at the level of significance of 0.05. All statistical analyses were performed using IBM SPSS Statistics 22.