Empirical data from a 6-month mesocosms experiment were used to assess the ability and performance of two machine learning (ML) models, including artificial neural network (NN) and random forest (RF), to predict temporal bioavailability changes of complex chemical mixtures in contaminated soils amended with compost or biochar. From the predicted bioavailability data, toxicity response for relevant ecological receptors was then forecasted to establish environmental risk implications and determine acceptable end-point remediation. The dataset corresponds to replicate samples collected over 180 days and analysed for total and bioavailable petroleum hydrocarbons and heavy metals/metalloids content. Further to this, a range of biological indicators including bacteria count, soil respiration, microbial community fingerprint, seeds germination, earthworm's lethality, and bioluminescent bacteria were evaluated to inform the environmental risk assessment. Parameters such as soil type, amendment (biochar and compost), initial concentration of individual compounds, and incubation time were used as inputs of the ML models. The relative importance of the input variables was also analysed to better understand the drivers of temporal changes in bioavailability and toxicity. It showed that toxicity changes can be driven by multiple factors (combined effects), which may not be accounted for in classical linear regression analysis (correlation). The use of ML models could improve our understanding of rate-limiting processes affecting the freely available fraction (bioavailable) of contaminants in soil, therefore contributing to mitigate potential risks and to inform appropriate response and recovery methods.