Deterioration of water quality is a major problem world widely according to many international non-governmental organizations (NGO). As one of the European Union (EU) countries, Bulgaria is also obliged by EU legislation to maintain best practices in assessing surface water quality and the efficiency of wastewater treatment processes. For these reasons studies were undertaken to utilize ecotoxicological (Microtox^®, Phytotoxkit F^TM, Daphtoxkit F^TM), instrumental (to determine pH, electrical conductivity (EC), chemical oxygen demand, total suspended solids (TSS), total nitrogen (N) and phosphorus (P), chlorides, sulphates, Cr, Co, Cu, Cd, Ba, V, Mn, Fe, Ni, Zn, Se, Pb), as well as advanced chemometric methods (partial least squares–discriminant analysis (PLS-DA)) in data evaluation to comprehensively assess wastewater treatment plants’ (WWTPs) effluents and surface waters quality around 21 major Bulgarian cities. The PLS-DA classification model for the physicochemical parameters gave excellent discrimination between WWTP effluents and surface waters with 93.65% correct predictions (with significant contribution of EC, TSS, P, N, Cl, Fe, Zn, and Se). The classification model based on ecotoxicological data identifies the plant test endpoints as having a greater impact on the classification model efficiency than bacterial, or crustaceans’ endpoints studied.

The goal of the present study is to assess the soil quality in Bulgaria using (i) an appropriate set of soil quality indicators, namely primary nutrients (C, N, P), acidity (pH), physical clay content and potentially toxic elements (PTEs: Cu, Zn, Cd, Pb, Ni, Cr, As, Hg) and (ii) respective data mining and modeling using chemometrical and geostatistical methods. It has been shown that five latent factors are responsible for the explanation of nearly 70% of the total variance of the data set available (principal components analysis) and each factor is identified in terms of its contribution to the formation of the overall soil quality—the mountain soil factor, the geogenic factor, the ore deposit factor, the low nutrition factor, and the mercury-specific factor. The obtained soil quality patterns were additionally confirmed via hierarchical cluster analysis. The spatial distribution of the patterns throughout the whole Bulgarian territory was visualized via the mapping of the factor scores for all identified latent factors. The mapping of identified soil quality patterns was used to outline regions where additional measures for the monitoring of the phytoavailability of PTEs were required. The suggested regions are located near to thermoelectric power plants and mining and metal production facilities and are characterized by intensive agricultural activity.

The temporal dynamics of anthropogenic impacts on the Pchelina Reservoir is assessed based on chemical element analysis of three sediment cores at a depth of about 100–130 cm below the surface water. The ¹³⁷Cs activity is measured to identify the layers corresponding to the 1986 Chernobyl accident. The obtained dating of sediment cores gives an average sedimentation rate of 0.44 cm/year in the Pchelina Reservoir. The elements’ depth profiles (Ti, Mn, Fe, Zn, Cr, Ni, Cu, Mo, Sn, Sb, Pb, Co, Cd, Ce, Tl, Bi, Gd, La, Th and U_nat) outline the Struma River as the main anthropogenic source for Pchelina Reservoir sediments. The principal component analysis reveals two groups of chemical elements connected with the anthropogenic impacts. The first group of chemical elements (Mn, Fe, Cr, Ni, Cu, Mo, Sn, Sb and Co) has increasing time trends in the Struma sediment core and no trend or decreasing ones at the Pchelina sampling core. The behavior of these elements is determined by the change of the profile of the industry in the Pernik town during the 1990s. The second group of elements (Zn, Pb, Cd, Bi and U_nat) has increasing time trends in Struma and Pchelina sediment cores. The increased concentrations of these elements during the whole investigated period have led to moderate enrichments for Pb and U_nat, and significant enrichments for Zn and Cd at the Pchelina sampling site. The moderately contaminated, according to the geoaccumulation indexes, Pchelina Reservoir surface sediment samples have low ecotoxicity.

Surface water quality strongly depends on anthropogenic activity. Among the main anthropogenic sources of this activity are the wastewater treatment plant (WWTP) effluents. The discharged loads of nutrients and suspended solids could provoke serious problems for receiving water bodies and significantly alter the surface water quality. This study presents inventory analysis and chemometric assessment of WWTP effluents based on the mandatory monitoring data. The comparison between the Bulgarian WWTPs and previously reported data from other countries reveals that discharged loads from investigated WWTPs are lower. This is particularly valid for total suspended solids (TSS). The low TSS loads are the reason for the deviations of the typical calculated WWTP effluent ratios of Bulgarian WWTPs compared to the WWTPs worldwide. The performed multivariate analysis reveals the hidden factors that determine the content of WWTP effluents. The source apportioning based on multivariate curve resolution analysis provides detailed information for source contribution profiles of the investigated WWTP effluent loads and elucidate the difference between WWTPs included in this study.

Wastewater treatment plant (WWTP) installations are designed and operated to reduce the quantity of pollutants emitted to surface waters receiving treated wastewaters. In this work, we used classical instrumental studies (to determine chemicals and parameters under obligations put with Directive 91/271/EEC), ecotoxicological tools (Sinapis alba root growth inhibition (SA-RG) and Heterocypris incongruens mortality (MORT) and growth inhibition (GRINH)) and multivariate statistical analysis to gain information on feature profiles of WWTPs’ effluents and the possible burden of surface water bodies receiving treated wastewaters in eleven locations of Bulgaria. Initial screening of results has shown that only phosphorus content exceeds the admissible level in 5 out 11 WWTP effluents, while As, Cr, Cu, Mn, and Zn show exceedance at several locations. The multivariate statistical analysis reveals the discriminating water quality parameters and outlines the ability of Heterocypris incongruens to evaluate the ecotoxicological potential of different groups of waters.