Mining activities have contaminated many riverine floodplains with arsenic (As). When floodplain soils become anoxic under water-saturated conditions, As can be released from the solid phase. Several microbially-driven As solubilization processes and numerous influential factors were recognized in the past. However, the interplay and relative importance of soil properties and the influence of environmental factors such as temperature remain poorly understood, especially considering the (co)variation of soil properties in a floodplain. We conducted anoxic microcosm experiments at 10, 17.5, and 25 °C using 65 representative soils from the mining-impacted Ogosta River floodplain in Bulgaria. To investigate the processes of As solubilization and its quantitative variation we followed the As and Fe redox dynamics in the solid and the dissolved phase and monitored a range of other solution parameters including pH, Eh, dissolved organic C, and dissolved Mn. We related soil properties to dissolved As observed after 20 days of microcosm incubation to identify key soil properties for As solubilization. Our results evidenced reductive dissolution of As-bearing Fe(III)-oxyhydroxides as the main cause for high solubilization. The availability of nutrients, most likely organic C as the source of energy for microorganisms, was found to limit this process. Following the vertical nutrient gradient common in vegetated soil, we observed several hundred μM dissolved As after 1–2 weeks for some topsoils (0–20 cm), while for subsoils (20–40 cm) with comparable total As levels only minor solubilization was observed. While high Mn contents were found to inhibit As solubilization, the opposite applied for higher temperature (Q10 2.3–6.1 for range 10–25 °C). Our results suggest that flooding of nutrient-rich surface layers might be more problematic than water-saturation of nutrient-poor subsoil layers, especially in summer floodings when soil temperature is higher than in winter or spring.

In 2005/6, nearly 3000 moss samples from (semi-)natural location across 16 European countries were collected for nitrogen analysis. The lowest total nitrogen concentrations in mosses (<0.8%) were observed in northern Finland and northern UK. The highest concentrations (≥1.6%) were found in parts of Belgium, France, Germany, Slovakia, Slovenia and Bulgaria. The asymptotic relationship between the nitrogen concentrations in mosses and EMEP modelled nitrogen deposition (averaged per 50 km × 50 km grid) across Europe showed less scatter when there were at least five moss sampling sites per grid. Factors potentially contributing to the scatter are discussed. In Switzerland, a strong (r² = 0.91) linear relationship was found between the total nitrogen concentration in mosses and measured site-specific bulk nitrogen deposition rates. The total nitrogen concentrations in mosses complement deposition measurements, helping to identify areas in Europe at risk from high nitrogen deposition at a high spatial resolution.

This study estimated current benzo(a)pyrene (BaP) concentration levels, population exposure and potential health impacts of exposure to ambient air BaP in Europe. These estimates were done by combining the best available information from observations and chemical transport models through the use of spatial interpolation methods. Results show large exceedances of the European target value for BaP in 2012 over large areas, particularly in central-eastern Europe. Results also show large uncertainties in the concentration estimates in regions with a few or no measurement stations. The estimation of the population exposure to BaP concentrations and its health impacts was limited to 60% of the European population, covering only the modelled areas which met the data quality requirement for modelling of BaP concentrations set by the European directive 2004/107/EC. The population exposure estimate shows that 20% of the European population is exposed to BaP background ambient concentrations above the EU target value and only 7% live in areas with concentrations under the estimated acceptable risk level of 0.12 ng m−3. This exposure leads to an estimated 370 lung cancer incidences per year, for the 60% of the European population included in the estimation. Emissions of BaP have increased in the last decade with the increase in emissions from household combustion of biomass. At the same time, climate mitigation policies are promoting the use of biomass burning for domestic heating. The current study shows that there is a need for more BaP measurements in areas of low measurement density, particularly where high concentrations are expected, e.g. in Romania, Bulgaria, and other Balkan states. Furthermore, this study shows that the health risk posed by PAH exposure calls for better coordination between air quality and climate mitigation policies in Europe.

Magnetic properties of road dusts from 26 urban sites in Bulgaria are studied. Temporal variations of magnetic susceptibility (χ) during eighteen months monitoring account for approximately 1/3rd of the mean annual values. Analysis of heavy metal contents and magnetic parameters for the fraction d < 63 μm reveal significant correlations (p < 0.05) between χ and Fe, Mn and PLI index. The highest negative correlation (R2 = −0.84) is observed between the ratio ARM/χ and Pb content. It suggests that Pb is related to brake/tyre wear emissions, releasing larger particles and higher Pb during slow driving – braking. Bulk χ values of road dusts per city show significant correlation with population size and mean annual NO2 concentration on a log-normal scale. The results demonstrate the applicability of magnetic measurements of road dusts for estimation of mean NO2 levels at high spatial density, which is important for pollution modelling and health risk assessment.

In this study Pb isotope signatures were used to identify the provenance of contaminant metals and establish patterns of downstream sediment dispersal within the River Maritsa catchment, which is impacted by the mining of polymetallic ores. A two-fold modelling approach was undertaken to quantify sediment-associated metal delivery to the Maritsa catchment; employing binary mixing models in tributary systems and a composite fingerprinting and mixing model approach in the wider Maritsa catchment. Composite fingerprints were determined using Pb isotopic and multi-element geochemical data to characterize sediments delivered from tributary catchments. Application of a mixing model allowed a quantification of the percentage contribution of tributary catchments to the sediment load of the River Maritsa. Sediment delivery from tributaries directly affected by mining activity contributes 42-63% to the sediment load of the River Maritsa, with best-fit regression relationships indicating that sediments originating from mining-affected tributaries are being dispersed over 200 km downstream.

In recent times many benthic indices have been proposed to assess the ecological quality of marine waters worldwide. In this study we compared single metrics and multi-metric methods to assess coastal and transitional benthic status along human pressure gradients in five distinct environments across Europe: Varna bay and lake (Bulgaria), Lesina lagoon (Italy), Mondego estuary (Portugal), Basque coast (Spain) and Oslofjord (Norway). Hence, 13 single metrics (abundance, number of taxa, and several diversity and sensitivity indices) and eight of the most common indices used within the European Water Framework Directive (WFD) for benthic assessment were selected: index of size spectra (ISS), Benthic assessment tool (BAT), Norwegian quality index (NQI), Multivariate AMBI (M-AMBI), Benthic quality index (BQI), (Benthic ecosystem quality index (BEQI), Benthic index based on taxonomic sufficiency (BITS), and infaunal quality index (IQI). Within each system, sampling sites were ordered in an increasing pressure gradient according to a preliminary classification based on professional judgement. The different indices are largely consistent in their response to pressure gradient, except in some particular cases (i.e. BITS, in all cases, or ISS when a low number of individuals is present). Inconsistencies between indicator responses were most pronounced in transitional waters (i.e. IQI, BEQI), highlighting the difficulties of the generic application of indicators to all marine, estuarine and lagoonal environments. However, some of the single (i.e. ecological groups approach, diversity, richness) and multi-metric methods (i.e. BAT, M-AMBI, NQI) were able to detect such gradients both in transitional and coastal environments, being these multi-metric methods more consistent in the detection than single indices. This study highlights the importance of survey design and good reference conditions for some indicators. The agreement observed between different methodologies and their ability to detect quality trends across distinct environments constitutes a promising result for the implementation of the WFD’s monitoring plans. Moreover, these results have management implications, regarding the dangers of misclassification, uncertainty in the assessment, use of conflicting indices, and testing and validation of indices.

The sources of marine litter (ML) pollution along the Bulgarian Black Sea coast in 2019 were identified and categorized to: Public litter, Fishing, Shipping, Sanitary and sewage, Fly tipped, Medical and Non-sourced, following Veiga et al., 2016 approach. Largest was the contribution of Public litter - 48.3% and Non-sourced - 38.5% of the total. Most Public litter items were attributed to recreational activities - 82.9% Public litter Recreation and less to Public litter Smoking-related - 17.1%. Regarding to Non-sourced, the following sub-categories seemed to contribute clearly more: Land (run off) - 8.79% and Tourism beach users - 8.58%. Considerably lower was the contribution of Fishing, Offshore and Shipping. None of the items were attributed to Sewage related. The ML pollution on the individual beaches (10 sites) exhibited largest contribution of Public litter with most significant input of Public litter Recreation, varying from 65.9 to 86.7% between beaches.

Cities located on Eastern Mediterranean is exposed to both local and distant anthropogenic and natural sources. In this study, to determine the effect of these sources on Particulate Matter (PM) concentrations in a coastal city, particulate matter with diameters less than 2.5 μm (PM2.5) and with diameters between 2.5 and 10 μm (PM2.5−10) were collected once in a two-day period for 24 h between July 2014 and July 2015 in downtown Antalya which is located on the Mediterranean coast of Turkey. Antalya is one of the fast growing city of Turkey with a population of 2.3 million. Samples were analyzed using an energy-dispersive X-ray fluorescence for 15 elements (Na, Mg, Al, Si, S, K, Ca, Ti, V, Mn, Fe, Cu, Zn, As, Pb). Statistical parameters were calculated for all measured elements in fine (PM2.5) and coarse (PM2.5−10) fraction. Crustal and marine elements, such as Ti, Ca, Al and Na were abundant in the course fraction. Only S was found in higher concentration in the fine fraction. Monthly variation of Crustal Enrichment Factor (EFC) results of Si showed that the area was under influence of non-local crustal dust especially during spring and late summer. EFC also indicated that during winter season, fine fraction K was due to local wood combustion. Source regions of S was determined using Potential Source Contribution Function (PSCF) and compared with previous studies conducted at a rural site of Antalya approximately twenty years ago. Most of the source regions affecting S concentrations at the Eastern Mediterranean region were found out to be same: western Anatolia, Marmara region, the Aegean Sea coasts of Greece and some parts of Bulgaria and Romania. However, due to decrease in SO2 emissions over the northeast coast of Black Sea and between Caspian Sea and Ukraine, the region was not turned up to be a source region.

This study aims to investigate the changes in atmospheric deposition trends in Bulgaria, studied using the moss biomonitoring technique since 1995. For the first time, a paired (site-wise) comparison was performed after a critical review of the sampling networks and adjusting for location, the distance between the sampling points, and moss species. Data from the 2005/2006 and 2015/2016 moss surveys were chosen as instrumental neutron activation analysis was employed in both to determine the content of 34 elements (Al, As, Ba, Br, Са, Ce, Cl, Со, Cr, Cs, Fe, Hf, I, K, La, Mn, Na, Nd, Ni, Rb, Sb, Sc, Se, Sr, Ta, Tb, Th, Ti, Tm, U, V, W, Yb, Zn). In addition, Cd, Cu, and Pb were determined using complementary analytical methods: inductively coupled plasma atomic emission spectroscopy in 2015/2016 and atomic absorption spectroscopy in 2005/2006. For the subset of 57 routinely sampled locations in Bulgaria, hierarchical clustering on principal components and multiple factor analysis (MFA) were applied to assess the spatial trends in the 10 years elapsed between the surveys, as well as to characterise the origin of the determined elements. Elevation and distance between the sampling points were used as additional variables in the multiple factor analysis plane to ascertain their effect on the overall variance in the datasets. Distribution maps were constructed to illustrate the deposition patterns for the pollutant Pb. The results were consistent with decreased industrial output in the country, increased coal combustion and transport pollution, and construction of roads.

The recent environmental debates in the European Green Deal to reduce pollution determine the reconsideration of the relationship between economic development and environmental quality. In this context, this paper explains the GHG emissions in few Central and Eastern European countries (CEEs) (Bulgaria, Poland, Hungary, Slovak Republic, Czech Republic and Romania) in the period 1990–2019 and proposes suitable economic policies to reduce pollution. Environmental Kuznets Curve (EKC) and renewable energy Kuznets Curve (RKC) are considered in a methodological framework based on panel threshold and dynamic panel models. The analysis based on GHG emissions from all economic sectors and GHG emissions from agriculture indicates similarities and differences related to the impact of economic indicators on pollution. The results of estimations indicated an inverted N-shaped relationship between GDP and GHG emissions and an N-shaped pattern between value added in agriculture and pollution. The U-shaped RKC was confirmed for overall economy and for agriculture. The policy proposals based on empirical findings promote sustainable development in the analysed states by reducing pollution.