Do transport taxes reduce air pollution in the top 10 countries with the highest transport tax revenues? A country-specific panel data analysis

Environmental Science and Pollution Research - Recently, policymakers have been developing some policies to improve environmental quality. Environmental taxes are the most important of these...     

Impacts of regional development on emissions in China’s transport sector

Environmental Science and Pollution Research - The CO2 emissions in China’s transport sector increased from 349.00 Mt in 2005 to 723.87 Mt in 2017. Thus, a number of climate change policies...     

Long-range transport episodes of fine particles in southern Finland during 1999–2007

The frequency, strength and sources of long-range transport (LRT) episodes of fine particles (PM_2.5) were studied in southern Finland using air quality monitoring results, backward air mass trajectories, remote sensing of fire hot spots, transport and dispersion modelling of smoke and chemical analysis of particle samples (black carbon, monosaccharide anhydrides, oxalate, succinate, malonate, SO₄^2?, NO₃^?, K⁺ and NH₄⁺). At an urban background site in Helsinki, the daily WHO guideline value (24-h PM_2.5 mean 25 μg m^?3) was exceeded during 1–7 LRT episodes per year in 1999–2007. The 24-h mean maximum concentrations varied between 25 and 49 μg m^?3 during the episodes, which was 3–6 times higher than the local mean concentration (8.7 μg m^?3) in 1999–2007. The highest particle concentrations (max. 1-h mean 163 μg m^?3) and the longest episodes (max. 9 days) were mainly caused by the emissions from open biomass burning, especially during springs and late-summers in 2002 and 2006. During the period 2001–2007, the satellite remote sensing of active fire hot spots and transport and dispersion modelling of smoke indicated that approximately half of the episodes were caused partly by the emissions from wildfires and/or agricultural waste burning in fields in Eastern Europe, especially in Russia, Belarus and Ukraine. Other episodes were mainly caused by the LRT of ordinary anthropogenic pollutants, e.g. from energy production, traffic, industry and wood combustion. During those ‘other episodes’, air masses also arrived from Eastern Europe, including Poland. The highest concentrations of biomass-burning tracers, such as monosaccharide anhydrides (levoglucosan + mannosan + galactosan) and K⁺, were observed during open biomass-burning episodes, but quite high values were also measured during some winter episodes due to wood combustion emissions. Our results indicate that open biomass burning in Eastern Europe causes high fine particle concentration peaks in large areas of Europe almost every year.     

Modelling base cations in Europe—sources,transport and deposition of calcium

The deposition of the base cations calcium, magnesium and potassium from the atmosphere needs to be quantified in the calculation of the total deposited acidity in the critical loads approach. Of these base cations, calcium has been found to be the most important in terms of mass deposited. However, the sources of calcium to the atmosphere are not well understood. Recently, the first spatially disaggregated inventory of industrial calcium emissions for Europe was presented by Lee and Pacyna (1998) who estimated a total European emission of 0.7–0.8 Mt yr^-1. However, it is thought that wind blown dust from soils contributes a substantial fraction to the deposition of calcium. In this work, the source strength of calcium from arid regions within the EMEP modelling domain was estimated using the global mineral dust emission data base of Tegen and Fung (1994) and an estimation of the calcium content of soils. This results in a “natural” calcium emission of 6 Mt yr^-1. A long-range transport model, TRACK, was used to calculate the wet and dry deposition of calcium arising from these industrial and natural sources to the UK which resulted in a total deposition of 29–30 kt yr^-1. Of this annual deposition, 0.6–0.7 kt arises from cement manufacturing, 0.02–0.03 kt from iron and steel manufacturing, 0.8–0.83 kt from a large point source power generation, and 28 kt from power generation from a small boiler plant. The natural emissions of calcium from arid regions result in a deposition of calcium to the UK of 0.5 kt yr^-1. The measured wet deposition of calcium to the UK is 89 kt yr^-1 and the estimated dry deposition 14 kt yr^-1. The short-fall in the modelled deposition of calcium is thus of the order of 70 kt yr^-1, which is suggested to arise from wind-blown dust from agricultural land in the UK and mainland Europe. The estimated emissions, and thus modelled deposition are rather uncertain, such that estimating deposition of calcium attributable to agricultural soil emissions by differencing has a large uncertainty. However, this is the first such study of its kind for Europe and represents a first step towards understanding the sources of calcium and their contribution to mitigating deposited acidity from acidifying pollutants such as sulphur dioxide, nitrogen oxides and ammonia.     

Persistence and transport of bacteria and viruses in groundwater — a conceptual evaluation

The results of an interdisciplinary research program on the subsurface persistence and transport of pathogenic bacteria and viruses are evaluated on the basis of an expanded advection-dispersion model which considers the persistence of these microorganisms under the conditions found within an aquifer, the retardation by adsorption-desorption processes and the role of filtration processes. The model indicates that the principal controls are filtration processes, the microorganisms being fixed on the filter media are ultimately eliminated or inactivated by biological, chemical and physical processes.     

Solute transport in crystalline rocks at Aspö--I: geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted Asp? Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours-days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption K(d)s are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.     

The long-range transport of aerosols from northern China to Hong Kong – a multi-technique study

The results of the inorganic and organic analyses of aerosol samples collected on the east and west sides of Hong Kong during a dust episode (9–10 May 1996) are reported. The origin of the dust was traced to Northern China. The dust reached Hong Kong by way of the East China Sea. The characteristics of the inorganic elements and organic compounds were quite different from the non-episodic samples collected on 1–2 April 1996, EPD (Environmental Protection Department, Special Administrative Region, Hong Kong, China) results for April–May 1994, and our early studies (Zheng et al., 1997. Atmospheric Environment 31(2), 227–237.). Results from X-ray spectrometry showed pronounced increase in the relative abundance of Al, Fe, Ca, S and Cl in the dust samples compared to the non-episodic samples. The high abundance of Cl in the dust samples suggested the aerosols experienced long-range transport by way of the sea. ICP-MS analysis revealed higher concentrations of Fe, Ca, S and Pb in the episodic samples relative to the values measured during April–May 1994 by EPD. The high Ca content in the soil samples is a characteristic of northern Chinese crustal material (Liu et al., 1985). Hong Kong aerosols are characterized by high octadecenoic acid concentration due to heavy urbanization and Chinese-style stir-fry cooking. A much lower C_18:1/C_18:0 ratio was found in the episodic samples, however, suggesting the aerosols were transported from a long distance. The high ratio of ⩾C₂₀/<C₂₀ in n-alkanoic acids, the high input of n-alkanes and n-alkanols from plant waxes, and the unusually similar distribution of the organic compounds in the east and west samples suggested the existence of non-local sources on 9–10 May 1996. The compositions and distributions of lipids in the aerosol samples collected during the episode corresponded well with those of the eolian dust samples over the Atlantic and Pacific Oceans. Back trajectories and low altitude (<3 km) mesoscale flow modeling were performed, suggesting the existence of a mesoscale atmospheric structure off the east China coast, which could be responsible for the dust episode. Surface charts indicated the presence of suspended dust near Shanghai on 8 May 1996. Taiwan also experienced a similar episode on 8–9 May 1996. This integrated, multi-technique approach provides clear evidence that the 9–10 May episode was attributed to a consequence of Asian Dust. This is the first scientific report of Asian Dust in Hong Kong.     

Stochastic simulation of meteorological variability for long-range atmospheric transport—II. Long-term statistical models

Methods are developed for simulating meteorological inputs for statistical atmospheric transport models. Stochastic models are presented for the annual mean wind vector, the horizontal diffusion coefficient, the fraction of time it is raining, and the quantity of precipitation. The models are derived by aggregating the statistical models for short-term meteorological processes given in the companion paper of Small and Samson (Atmospheric Environment, 1989, 23, 2813–2824). Normal distributions are used for each variable. The normal distribution is derived directly for the mean wind vector based on the normality of short-term winds, and serves as an acceptable approximation for the other variables given their low interannual coefficient of variation. Correlation between variables is considered, including correlation between winds and precipitation, between the fraction of time it is raining and the total quantity of precipitation at a site, and between the annual precipitation quantity at multiple sites. Development and confirmation of the interannual distributions is supported by long-term meteorological data, including a nine-year backward trajectory record computed for a site in western Pennsylvania. Application of the method is demonstrated with a statistical atmospheric LRT model used to predict distributions of annual precipitation sulfates in eastern North America. Predicted interannual coefficients of variation of wet SO₄²⁻ concentrations range from about 10 to 15 per cent, similar to the level of interannual variation inferred from long-term data sets.     

Monitoring OH-PCBs in PCB transport worker’s urine as a non-invasive exposure assessment tool

In this study, we analyzed hydroxylated polychlorinated biphenyls (OH-PCBs) in urine of both PCB transport workers and PCB researchers. A method to monitor OH-PCB in urine was developed. Urine was solid-phase extracted with 0.1% ammonia/ methanol (v/v) and glucuronic acid/sulfate conjugates and then decomposed using β-glucuronidase/arylsulfatase. After alkaline digestion/derivatization, the concentration of OH-PCBs was determined by HRGC/HRMS-SIM. In the first sampling campaign, the worker’s OH-PCB levels increased several fold after the PCB waste transportation work, indicating exposure to PCBs. The concentration of OH-PCBs in PCB transport workers’ urine (0.55~11 μg/g creatinine (Cre)) was higher than in PCB researchers’ urine (< 0.20 μg/g Cre). However, also a slight increase of OH-PCBs was observed in the researchers doing the air sampling at PCB storage area. In the second sampling, after recommended PCB exposure reduction measures had been enacted, the worker’s PCB levels did not increase during handling of PCB equipment. This suggests that applied safety measures improved the situation. Hydroxylated trichlorobiphenyls (OH-TrCBs) were identified as a major homolog of OH-PCBs in urine. Also, hydroxylated tetrachlorobiphenyls (OH-TeCBs) to hydroxylated hexachlorobiphenyls (OH-HxCBs) were detected. For the sum of ten selected major indicators, a strong correlation to total OH-PCBs were found and these can possibly be used as non-invasive biomarkers of PCB exposure in workers managing PCB capacitors and transformer oils. We suggest that monitoring of OH-PCBs in PCB management projects could be considered a non-invasive way to detect exposure. It could also be used as a tool to assess and improve PCB management. This is highly relevant considering the fact that in the next 10 years, approx. 14 million tons of PCB waste need to be managed. Also, the selected populations could be screened to assess whether exposure at work, school, or home has taken place.     

A reactive transport model for mercury fate in soil—application to different anthropogenic pollution sources

Soil systems are a common receptor of anthropogenic mercury (Hg) contamination. Soils play an important role in the containment or dispersion of pollution to surface water, groundwater or the atmosphere. A one-dimensional model for simulating Hg fate and transport for variably saturated and transient flow conditions is presented. The model is developed using the HP1 code, which couples HYDRUS-1D for the water flow and solute transport to PHREEQC for geochemical reactions. The main processes included are Hg aqueous speciation and complexation, sorption to soil organic matter, dissolution of cinnabar and liquid Hg, and Hg reduction and volatilization. Processes such as atmospheric wet and dry deposition, vegetation litter fall and uptake are neglected because they are less relevant in the case of high Hg concentrations resulting from anthropogenic activities. A test case is presented, assuming a hypothetical sandy soil profile and a simulation time frame of 50 years of daily atmospheric inputs. Mercury fate and transport are simulated for three different sources of Hg (cinnabar, residual liquid mercury or aqueous mercuric chloride), as well as for combinations of these sources. Results are presented and discussed with focus on Hg volatilization to the atmosphere, Hg leaching at the bottom of the soil profile and the remaining Hg in or below the initially contaminated soil layer. In the test case, Hg volatilization was negligible because the reduction of Hg²⁺ to Hg⁰ was inhibited by the low concentration of dissolved Hg. Hg leaching was mainly caused by complexation of Hg²⁺ with thiol groups of dissolved organic matter, because in the geochemical model used, this reaction only had a higher equilibrium constant than the sorption reactions. Immobilization of Hg in the initially polluted horizon was enhanced by Hg²⁺ sorption onto humic and fulvic acids (which are more abundant than thiols). Potential benefits of the model for risk management and remediation of contaminated sites are discussed.     

Inventory and substance flow analysis of polybrominated diphenyl ethers in the Nigerian transport sector—end-of-life vehicles policy and management

Recently, certain polybrominated diphenyl ethers (PBDEs) have been listed as persistent organic pollutants (POPs) in the Stockholm Convention. In this study, a preliminary material and substance flow analysis of commercial pentabromodiphenyl ether (c-PentaBDE) was conducted for motor vehicles—a major use sector for POP-PBDE in polyurethane (PUR) foam—for Nigeria. The methodology of the Stockholm Convention PBDE inventory guidance was used for the calculation of c-PentaBDE. Material/substance flow analysis was conducted applying the STAN software. The time frame for this analysis was 1980–2010, considering that this was the period when POP-PBDE-containing vehicles were largely imported into Nigeria.It is estimated that the approximately 19 million passenger cars imported from 1980 to 2010 contained ca. 270 t of POP-PBDEs in ca. 401,000 t of PUR foam. A major share of cars from the USA and only a small share of cars from Europe and Asia were impacted. This simplified material and substance flow of PUR foams and POP-PBDEs in motor vehicles demonstrated the potential for environmental/human contamination and pollution of recycling/reuse for Nigeria and other developing countries.The study developed the first preliminary inventory of end-of-life vehicles for Nigeria, following which the environment ministry has taken up the important issue of end-of-life vehicles management. Considering that a range of other pollutants are contained in vehicles (e.g. heavy metals, flame retardants, PCBs, chlorofluorocarbons, hydrofluorocarbons and waste oil), this initiated activity should finally lead to an integrated management of pollutants and resources from the transport sector.     

Groundwater chemistry of the Okélobondo uraninite deposit area (Oklo, Gabon): two-dimensional reactive transport modelling

The stability of uranium-bearing minerals in natural environments is of interest to evaluate the feasibility of radioactive waste repositories. The uraninite bodies, UO2(s), in the Oklo district (Gabon) are the result of a natural fission process, which took place 1970 Ma ago. These deposits can be regarded as natural analogues for spent fuel. One of the uraninite bodies, the Okélobondo deposit, is located at a depth of 300 m. Groundwater samples from boreholes located at shallow depths (100-200 m) show neutral to basic pH, anoxic conditions (Eh = 0.10 to -0.05 V) and are saturated with respect to uraninite. In contrast, deeper samples collected in the vicinity of the ore body are oxidising (Eh = 0.32-0.47 V), slightly basic (pH = 7.0-8.5) and undersaturated with respect to uraninite. These oxidising conditions at depth, if present under repository conditions, may affect the stability of uranium oxide. In order to improve our understanding of the observed site geochemistry, the available information on the lithology and groundwater flow was integrated in a reactive transport model. The chemical composition and the pH-Eh values of the water sampled above and in the western side of the Okélobondo deposit can be explained by the interaction of meteoric recharge with pelites, dolomites and sandstones. The dissolution of Fe(II)-silicates and the oxidation of the Fe(II)-aqueous species maintained the pH-Eh distribution along the Fe(2+)-Fe(OH)3(am) equilibrium, with the result that uraninite does not dissolve. This may explain the lower uranium content in the water samples from pelites and dolomites above the Okélobondo deposit. The high Mn/Fe ratio and the high pH-Eh values of the water sampled at depth, close to the Okélobondo deposit, suggest a control by the Mn(2+)-MnOOH(s) equilibrium. This control is attributed to the dissolution of a large rhodochrosite, MnCO3(s), and manganite, MnOOH(s) deposit in the recharge area on the eastern side.     

The assessment of the soil–plant-animal transport of the risk elements at the locations affected by brown coal mining

Environmental Science and Pollution Research - The North Bohemian Brown Coal Basin (Czech Republic) is suggested as a source of significant pollution in the surrounding environment with various...     

Applied tracer tests in fractured rock: Can we predict natural gradient solute transport more accurately than fracture and matrix parameters?

Applied tracer tests provide a means to estimate aquifer parameters in fractured rock. The traditional approach to analysing these tests has been using a single fracture model to find the parameter values that generate the best fit to the measured breakthrough curve. In many cases, the ultimate aim is to predict solute transport under the natural gradient. Usually, no confidence limits are placed on parameter values and the impact of parameter errors on predictions of solute transport is not discussed. The assumption inherent in this approach is that the parameters determined under forced conditions will enable prediction of solute transport under the natural gradient. This paper considers the parameter and prediction uncertainty that might arise from analysis of breakthrough curves obtained from forced gradient applied tracer tests. By adding noise to an exact solution for transport in a single fracture in a porous matrix we create multiple realisations of an initial breakthrough curve. A least squares fitting routine is used to obtain a fit to each realisation, yielding a range of parameter values rather than a single set of absolute values. The suite of parameters is then used to make predictions of solute transport under lower hydraulic gradients and the uncertainty of estimated parameters and subsequent predictions of solute transport is compared. The results of this study show that predictions of breakthrough curve characteristics (first inflection point time, peak arrival time and peak concentration) for groundwater flow speeds with orders of magnitude smaller than that at which a test is conducted can sometimes be determined even more accurately than the fracture and matrix parameters.     

Mesoscale modelling of vertical atmospheric transport in the Alps associated with the advection of a tropopause fold – a winter ozone episode

The role of vertical atmospheric transport (VAT) for a winter ozone episode (January 1988) in the Alps, observed both on the mountain crest (3580 m asl) and in the lee of the Alps (209 m asl), is investigated. Numerical simulations were carried out with a doubly nested mesoscale model (14 km horizontal resolution), characterised by a refined orography scheme. A comparison between the modelled vertical profiles with features observed in the aerological soundings over Payerne indicated that the model is able to reproduce the fold in the Alps down to 650 hPa. Frictional processes are found to be responsible for the fragmentation of the lower parts of the fold when it is advected over the mountains, thus enhancing ozone concentrations at the mountain crest with likely cross-tropopause mixing. Vertically propagating gravity waves and strong mesoscale vertical winds (nordfoehn) are responsible for the further downward transport of ozone-rich air to the lower troposphere. This mechanism leads to ozone concentrations up to 70 ppb(v) in the leeside. The study is important because tropospheric ozone trends in the lee of the Alps have been inferred, in the past, from ozone peaks associated with such nordfoehn conditions. It is now shown that nordfoehn may entrain air which has been recently exchanged across the tropopause.     

Solute transport in crystalline rocks at Aspö--II: blind predictions,inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish Asp? Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na(+), 85Sr(2+), 47Ca(2+)and more strongly sorbing 86Rb(+), 133Ba(2+), 137Cs(+). Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2-15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted K(d) values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25-60, but agree within a factor of 3-5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning K(d)s is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.     

Influence of usage and chemical–physical properties on the atmospheric transport and deposition of pesticides to agricultural regions of Manitoba,Canada

This study quantified the masses of 14 pesticides deposited as wet (precipitation) versus dry (gaseous and particle) atmospheric deposition at a research farm in southwestern Manitoba, Canada. The concentration in air of these pesticides was also measured. Total bulk deposition amounts (wet + dry) ranged from 0.009 to 2.3 μg m^?2 for the 12 pesticides detected, and for the six pesticides with both wet and dry deposition detections, dry deposition contributed 12–51% of the total deposition over the crop growing season. Although not applied at the site, eight herbicides registered for use in Canada, as well as lindane (γ-HCH), were all frequently detected (92–100%) in the 12 air samples analyzed during the crop growing season, with by-product isomer α-HCH (75%), clopyralid (50%) and atrazine (8%) detected to a lesser extent. The chemical’s physicochemical properties and the relative mean mass of each agricultural pesticide applied in the province of Manitoba and in a 13 km radius were significant parameters in explaining the trends in the concentrations of pesticides detected in our samples. The important contribution of dry deposition to total pesticide deposition warrants greater attention in arid and semi-arid areas such as the Prairie Region of Canada, also because under a changing climate this region is estimated to experience more severe droughts while the more favorable conditions predicted for pest infestations could lead to increased pesticide applications in agricultural and urban areas.     

Sorption–desorption and transport of trimethoprim and sulfonamide antibiotics in agricultural soil: effect of soil type,dissolved organic matter,and pH

Use of animal manure is a main source of veterinary pharmaceuticals (VPs) in soil and groundwater through a series of migration processes. The sorption–desorption and transport of four commonly used VPs including trimethoprim (TMP), sulfapyridine, sulfameter, and sulfadimethoxine were investigated in three soil layers taken from an agricultural field in Chongming Island China and two types of aqueous solution (0.01 M CaCl₂ solution and wastewater treatment plant effluent). Results from sorption–desorption experiments showed that the sorption behavior of selected VPs conformed to the Freundlich isotherm equation. TMP exhibited higher distribution coefficients (K _d?=?6.73–9.21) than other sulfonamides (K _d?=?0.03–0.47), indicating a much stronger adsorption capacity of TMP. The percentage of desorption for TMP in a range of 8–12 % is not so high to be considered significant. Low pH (<pK _a of tested VPs) and rich soil organic matter (e.g., 0–20 cm soil sample) had a positive impact on sorption of VPs. Slightly lower distribution coefficients were obtained for VPs in wastewater treatment plant (WWTP) effluent, which suggested that dissolved organic matter might affect their sorption behavior. Column studies indicated that the transport of VPs in the soil column was mainly influenced by sorption capacity. The weakly adsorbed sulfonamides had a high recovery rate (63.6–98.0 %) in the leachate, while the recovery rate of TMP was only 4.2–10.4 %. The sulfonamides and TMP exhibited stronger retaining capacity in 20–80 cm and 0–20 cm soil samples, respectively. The transport of VPs was slightly higher in the columns leached by WWTP effluent than by CaCl₂ solution (0.01 M) due to their sorption interactions.