Screening 31 endocrine-disrupting pesticides in water and surface sediment samples from Beijing Guanting reservoir

For screening 31 potential or suspected endocrine-disrupting pesticides in water and surface sediments, a multiresidue analysis method based on gas chromatography with electron capture detection (GC/ECD) was developed. Solid phase extraction (SPE) technology with Oasis^® HLB cartridge was also applied in sample extraction. The relevant mean recoveries were 70–103% and 71–103% for water and sediment, respectively. Relative standard deviations (RSD) are 2.0–7.0%, 4.0–8.0% for water and sediment, respectively. Thirty one pesticides (-HCH, β-HCH, γ-HCH, δ-HCH, hexachlorobenzene (HCB), aldrin, heptachlor, endosulfan I & II, p,p′-DDD, o,p′-DDT, p,p′-DDT, p,p′-DDE, endrin aldehyde, endosulfan sulphate, methoxychlor, hepachlor epoxide, -chlordane, γ-chlordane, dieldrin, endrin, dicofol, acetochlor, alachlor, metolachlor, chlorpyriphos, nitrofen, trifluralin, cypermethrin, fenvalerate, deltamethrin) in water and surface sediment samples from Beijing Guanting reservoir were analyzed. Concentrations of pesticides ranged from 7.59 to 36.0 ng g⁻¹ on a dry wt. basis for sediment samples, from 279.3 to 2740 ng l⁻¹ for pore waters and from 48.8 to 890 ng l⁻¹ for water samples, respectively, with a mean concentration of 10.7 ng g⁻¹ in sediment, 735 ng l⁻¹ in pore water and 295 ng l⁻¹ in water, respectively. The data obtained provides information on the levels and sources of endocrine-disrupting pesticides in Guanting reservoir. These results underscore the need to improved environmental protection measures in order to reduce the exposure of the population and aquatic biota to these endocrine-disrupting compounds.     

Organochlorine pesticides accumulation and degradation products in vegetation samples of a contaminated area in Galicia (NW Spain)

The content of 21 organochlorine pesticides were studied in vegetation samples of a highly contaminated area by isomers of hexachlorocyclohexane (HCH) located close to a former industrial area in Galicia (NW Spain). Five species of plants were collected at different points of the contaminated area and the different parts of the plants were separated in order to study differences in accumulation capabilities. Samples were extracted employing microwave energy followed by a clean-up step using solid phase extraction and finally determined by GC–ECD. The results obtained show that the most abundant pesticides are HCHs isomers, being the main isomers β-HCH and -HCH in all samples whereas δ-HCH and γ-HCH were at lower levels. Some other pesticides such as p,p′-DDT, p,p′-DDD and p,p′-DDE were also present in much lower amount in some of the samples. Several degradation products of HCH were also identified in some samples by GC–MS.     

Effect of some heavy metals and soil humic substances on the phytochelatin production in wild plants from silver mine areas of Guanajuato, Mexico

Phytochelatins (PCs) were determined in the wild plants, focusing on their relationship with the levels of heavy metals and humic substances (HS) in soil. Ricinus communis and Tithonia diversifolia were collected from several sites in Guanajuato city (Mexico), which had long been the silver and gold mining center. The analysis of PCs in root extracts was carried out by liquid chromatography (derivatization with monobromobimane). Total Ag, Cd, Cu and Pb in plant roots and in soil samples, as well as soil HS were determined. The association of metals with HS in soils was evaluated by size exclusion chromatography (SEC) with UV and ICP-MS detection. The results obtained revealed the induction of PCs in R. communis but not in T. diversifolia. The levels of Cd and Pb in plant roots presented strong positive correlation with PC-2 (r = 0.9395, p = 0.005; r = 0.9573, p = 0.003, respectively), indicating that these two metals promote PCs induction in R. communis. On the other hand, the inverse correlation was found between soil HS and metal levels in roots of R. communis (Cu > Pb > Cd > Ag), in agreement with the decreasing affinity of these metals to HS. Importantly, the inverse correlation between soil HS and plant PC-2 was observed (r = −0.7825, p = 0.066). These results suggest that metals strongly bound to HS could be less bioavailable to plants, which in turn would limit their role in the induction of PCs. Indeed, the SEC elution profiles showed Pb but not Cd association with HS and the correlation between metal in soil and PC-2 in plant was statistically significant only for Cd (r = 0.7857, p = 0.064). Based on these results it is proposed that the role of heavy metals in PCs induction would depend on their uptake by R. communis, which apparently is controlled by the association of metals with soil HS. This work provides further evidence on the role of environmental conditions in the accumulation of heavy metals and phytochelatin production in plants.     

Levels of organochlorine pesticides in human milk from central Taiwan

The present study determined the residues of organochlorine pesticides (OCPs) in human milk collected in central Taiwan between December 2000 and November 2001. The OCPs were analyzed by GC/MS for 36 human milk samples from healthy women ranging between 20 and 36 years of age. The predominant OCPs were p,p′-DDE, p,p′-DDT, -CHL, heptachlor epoxide, heptachlor, β-HCH, and γ-HCH, with median levels of 228, 19, 7.4, 4.0, 2.3, 1.2, and 0.8 ng/g lipid, respectively. The residues of OCPs in human milk from central Taiwan were comparable to those described in results from Sweden, the United Kingdom, and Japan, and were significantly lower than those from investigations in Asian countries, including China, Thailand, Indonesia, and Vietnam. Low DDE/DDT ratio (mean = 13.6, SD = 6.54) indicated that residual OCPs in human milk mainly originate from past exposure. A notable decrease in DDT levels (∑ DDT = 333 ng/g lipid) in human milk was found in this study compared to results from the previous two decades (∑ DDT = 3595 ng/g lipid). Hypothetically, the level of -CHL was significantly associated with total TEQ levels in Taiwanese human milk because of the sources of food contaminant, i.e. animal fat. Based on low OCP levels in Taiwanese human milk and low estimated median daily intake of total DDTs for a breastfed infant (1358 ng/kg/day) with the assumption of an infant weighting 4 kg and consuming 699 g milk per day in the first month after birth, the Taiwanese policy of breast-feeding promotion was supported.     

Cold season CH4, CO2 and N2O fluxes from freshwater marshes in northeast China

Cold season (winter and thaw) CH₄, CO₂ and N₂O fluxes from freshwater marshes (47°35′N, 133°31′E, Northeast China) were measured, using the static chamber method. The mean CH₄ and CO₂ fluxes from Carex lasiocarpa (Cl) were 0.5 ± 0.19 and 6.23 ± 1.36 mg C m⁻² h⁻¹, respectively, and those from Deyeuxia angustifoli (Da) were 0.18 ± 0.15 and 5.22 ± 2.48 mg C m⁻² h⁻¹, respectively in winter. There was no significant difference between Cl and Da (p > 0.05). The contributions of winter CH₄ fluxes were about 5.5% and 3% in the Cl and Da, respectively. Marshes are an important potential N₂O sink in winter season in northeast China. During thaw, the CH₄ and CO₂ emissions rapidly increased, 4.5–6 times of winter emissions. Wetland became a source of N₂O. Cold season gases flux from northern wetlands play an important role in the seasonal gas exchange.     

Bioavailability of zinc in runoff water from roofing materials

Corrosion and runoff from zinc-coated materials and outdoor structures is an important source for the dispersion of zinc in the environment. Being part of a large inter-disciplinary research project, this study presents the bioavailability of zinc in runoff water immediately after release from the surface of 15 different commercially available zinc-based materials exposed to the urban environment of Stockholm, Sweden. Runoff water was analysed chemically and evaluated for its possible environmental impact, using both a biosensor test with the bacteria Alcaligenes eutrophus (Biomet^®) and the conventional 72 h growth inhibition test with the green alga Raphidocelis subcapitata. Chemical speciation modelling revealed that most zinc (94.3–99.9%) was present as the free Zn ion, the most bioavailable speciation form. These findings were confirmed by the results of the biosensor test (Biomet^®) which indicated that all zinc was indeed bioavailable. Analysis of the ecotoxicity data also suggested that the observed toxic effects were due to the presence of Zn²⁺ ions. Finally, regression analysis showed that, for this type of runoff samples, the rapid screening biosensor was capable of predicting (a) the total amount of zinc present in the runoff samples (R² of 0.93–0.98; p<0.05) and (b) the observed 72 h-EbC₅₀s (R² of 0.69–0.97; p<0.05).     

Multi-generation studies with Chironomus riparius--effects of low tributyltin concentrations on life history parameters and genetic diversity

While toxicological data are available for numerous chemicals from standard tests, little is known on effects of pollutants over several generations or regarding chronic effects of chemicals on genetic diversity. Within the experiments, effects of the model pollutant tributyltin (TBT) were investigated over eleven generations at a sublethal TBT concentration of 4.46 μg as Sn kg⁻¹ sediment dw on life-cycle parameters and genetic variability of Chironomus riparius. Moreover, the adaptation potential towards TBT was determined. This experimental design enables the identification of TBT effects on life-cycle parameters and the determination of a potential extinction risk caused by chronic exposure. Furthermore, effects on the genetic structure can be determined, which are not predictable based solely on knowledge of the toxic mode of action of the chemical.

Genetic variety was determined via microsatellite analysis, measuring individual length differences of highly variable satellite DNA fragments. For the identification of changes in tolerances towards the stressor, acute and chronic toxicity experiments were conducted.

During the multi-generation study, significant effects on development and reproduction were determined. For some generations, the emergence was significantly (p < 0.05) delayed under TBT exposure. Reproduction seems to be a sensitive parameter as well, whereby females laid significantly larger egg masses (p < 0.05) in the latter generations. TBT did not affect the population growth rate nor the genetic variability, while clear deviations from the Hardy–Weinberg equilibrium appeared. The study also provides strong evidence for the acquirement of a higher tolerance towards the stressor in the TBT-exposed group.     

An investigation of adsorption at the air-water and soil-water interfaces for non-micellizing ethylene oxide-propylene oxide surfactants

Adsorption at the air–water interface and soil sorption from aqueous solution have been investigated for a group of ethylene oxide (EO)–propylene oxide (PO) block copolymeric surfactants. The group which have a common structural formula of EO_m PO_n EO_m is distinguished by the fact that they have large critical micelle concentration (CMC) values and therefore do not readily form micelles at common environmental concentrations and temperatures. Adsorption at the air–water interface is readily shown to be driven by the size of the hydrophobic PO block. The size of the reduction in surface tension produced by a common concentration of 10⁻⁵ mol dm⁻³ linearly increases with the size of the PO block as does the efficiency of adsorption at the air–water interface as measured by pC₂₀ – the negative logarithm of the surfactant concentration that produces a reduction in surface tension of 20 mN m⁻¹. Soil sorption data have also been captured for these compounds and the data are readily fitted to the Freundlich adsorption isotherm. However soil sorption is shown to be inversely related to the molecular mass of the molecules and appears to be related to the size of the hydrophilic EO blocks in the molecule.     

Anaerobic biodegradation of 4-alkylphenols in a paddy soil microcosm supplemented with nitrate

Anaerobic degradation of phenol, p-cresol, 4-n-propylphenol (n-PP), 4-i-propylphenol (i-PP), 4-n-butylphenol (n-BP) and 4-sec-butylphenol (sec-BP) was observed in a paddy soil supplemented with nitrate. We detected the metabolites 4′-hydroxypropiophenone (HPP) from n-PP, 4-i-propenylphenol from i-PP, and 4-(1-butenyl)phenol and 4′-hydroxybutyrophenone (HBP) from n-BP. Compared with the original soils, Betaproteobacteria became predominant in the microcosm during the degradation of phenol and p-cresol whereas no remarkable change was observed in the community degrading propylphenols and butylphenols. The microcosm, however, did not degrade 4-t-butylphenol (t-BP), 4-t-octylphenol (t-OP) and 4-n-octylphenol (n-OP). Paddy soil supplemented with sulfate or iron (III) as electron acceptors did not degrade phenol and 4-alkylphenols with the exception of the degradation of p-cresol in sulfate-reducing conditions. It was demonstrated for the first time that anaerobic microbial degradation of alkylphenols, in a paddy soil supplemented with nitrate as an electron acceptor, occurred via oxidation of the alpha carbon in the alkyl chain.     

Structure relationship for catalytic dechlorination rate of dichlorobenzenes in water

Three isomers of dichlorobenzene (o-, m- and p-DCB) were dechlorinated by Pd/Fe catalyst in aqueous solutions through catalytic reduction. The dechlorination reaction took place on the surface site of the catalyst via a pseudo-first-order kinetics, and resulted in benzene as the final reduction product. The rate constants of the reductive dechlorination for the three dichlorobenzenes (DCBs) in the presence of Pd/Fe as a catalyst were measured experimentally. In all cases, the reaction rate constants were found to increase with the decrease in the Gibbs free energy of the formation of DCBs. The reaction rate constant for o-, m- and p-DCBs in the presence of 0.020% (w/w) Pd/Fe at 25 °C was determined to be 0.0213, 0.0223, and 0.0254 min⁻¹, respectively. While the activation energy of each dechlorination reaction was measured to be 102.5, 96.6 and 80.0 kJ mol⁻¹ for o-, m- and p-DCBs, respectively. The results demonstrated that p-DCBs were reduced more easily than o- or m-DCBs, and the order of the tendency of the dechlorination was p-DCB > m-DCB > o-DCB. The presented data show the catalytic reduction using Pd/Fe as a catalyst is a fast and easy approach for the dechlorination of DCBs.     

Comparison of PDFs, closure schemes and turbulence parameterisations in Lagrangian stochastic models

Six one-dimensional models, based on the Ito-type stochastic equation, are presented and compared. Four of these take into account up to the fourth order moment of vertical velocity fluctuations, and two up to the third order moment. Four models make use of a bi-Gaussian probability density function (PDF) and the other two are based on a Gram-Charlier series expansion truncated to the third or fourth order. All the models were run with a parameterisation of input turbulence (i.e. w², w³, and τ profiles). Concerning the fourth order moment w⁴, two different parameterisations were considered. Comparisons are made between ground-level concentrations, plume height and plume width observed in the Willis and Deardorff water tank experiments and those predicted by the different models here considered. The goal of this study was to find the models that give greater confidence in their applicability in dispersion studies and to verify the importance of considering the fourth order moment. The main conclusions are: simulation results largely depend on the turbulence parameterisation chosen; the Gram-Charlier PDF gives the best agreement with observations; some combinations of models and turbulence parameterisations perform well in simulating the shape of the ground-level concentration (g.1.c.) trend but fail in correctly simulating the form of the plume (plume height and vertical width); in the case of the Gram-Charlier PDF, the fourth order model reproduced the vertical plume width better than the third order one, whereas the two schemes yielded similar g.1.c. distributions.     

Tissue distribution of polychlorinated naphthalenes (PCNs) and non-ortho chlorinated biphenyls (non-ortho CBs) in harbour porpoises (Phocoena phocoena) from Swedish waters

Polychlorinated naphthalenes (PCNs) and non-ortho chlorinated biphenyls (non-ortho CBs) were analysed in blubber, nuchal fat, liver, muscle, kidney and brain of three male harbour porpoises (Phocoena phocoena) from the west coast of Sweden. To estimate spatial variation, PCNs and non-ortho CBs were analysed in six blubber samples collected at different anatomical sites of each animal. Highest wet weight concentrations of ΣPCNs were detected in the lipid rich tissues (blubber and nuchal fat) and liver (520–730 and 520 pg/g, respectively) and lowest in brain (22 pg/g). TetraCNs were most abundant in muscle, kidney and brain, while the hexaCNs were most abundant in the lipid rich tissues and liver. The highest lipid weight concentration recorded (11 ng/g) was for the hexaCN congeners no. 66/67 in liver. These coeluting hexaCN congeners accounted for 80–100% of total hexaCNs in all tissues examined.

Concentrations of Σnon-ortho CBs were highest in lipid rich tissues (220–280 pg/g wet weight). Non-ortho CB no. 77 and 169 constituted between 62–86% and 4.9–9.3%, respectively, of total Σnon-ortho CBs. No major variation of Σnon-ortho CB concentrations was found between the six different blubber sites but higher ΣPCN concentrations (wet weight) were found dorsally at the peduncle. Toxic equivalent concentrations (TEQs) showed that non-ortho CB no. 126 was the main contributor to total TEQs in all tissues, except liver in which hexaCN congener nos. 66/67 contributed to about 50% of total TEQs.     

Bioconcentration of metals in the moss Scleropodium purum in the area surrounding a power plant. A geotopographical predictive model for mercury

Samples of the moss Scleropodium purum collected in 1995 and 1997 were used to biomonitor the deposition of metals in the area surrounding a thermal power plant. Significantly higher levels of Cu (p<0.05), Fe (p<0.01), As and Hg (p<0.001) were found in the 1997 samples than in the 1995 samples, due to changes in atmospheric conditions. The influence on bioconcentration of the orientation of the sampling sites relative to the source of emission was studied. It was found that the increase recorded in 1997 generally occurred in the sampling sites in the south east of the study area. Analysis of the effect of distance from the source of emission revealed that the increase in metal levels in 1997 took place close to the power station (10–30 km). Finally, multiple regression analysis was used to construct a model that related different topographical variables to the concentrations of Hg in moss. The model, constructed using the data collected in both sampling periods, included the orientation of the sampling sites relative to the source of emission as well as the height of sampling sites in 1995 and the distance from the emission point in 1997. The model allowed us to determine the extent of the area affected by deposition and to establish the magnitude of deposition.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Modelling the environmental degradation of water contaminants. Kinetics and mechanism of the riboflavin-sensitised-photooxidation of phenolic compounds

The aerobic visible-light-photosensitised irradiation of methanolic solutions of either of the phenolic-type contaminants model compounds (ArOH) p-phenylphenol (PP), p-nitrophenol (NP) and phenol (Ph), and for two additional phenolic derivatives, namely p-chlorophenol (ClP) and p-methoxyphenol (MeOP), used in some experiments, was carried out. Employing the natural pigment riboflavin (Rf) as a sensitiser, the degradation of both the ArOH and the very sensitiser was observed. A complex mechanism, common for all the ArOH studied, operates. It involves superoxide radical anion (O₂^√−) and singlet molecular oxygen (O₂(¹Δ_g)) reactions. Maintaining Rf in sensitising concentrations levels (≈0.02 mM), the mechanism is highly dependent on the concentration of the ArOH. Kinetic experiments of oxygen and substrate consumption, static fluorescence, laser flash photolysis and time-resolved phosophorescence detection of O₂(¹Δ_g) demonstrate that at ArOH concentrations in the order of 10 mM, no chemical transformation occurs due to the complete quenching of Rf singlet excited state. When ArOH is present in concentrations in the order of mM or lower, O₂^√− is generated from the corresponding Rf radical anion, which is produced by electron transfer reaction from the ArOH to triplet excited Rf. The determined reaction rate constants for this step show a fairly good correlation with the electron-donor capabilities for Ph, PP, NP, ClP and MeOP. In this context, the main oxidative species is O₂^√−, since O₂(¹Δ_g) is quenched in an exclusive physical fashion by the ArOH. The production of O₂^√− regenerates Rf impeding the total degradation of the sensitiser. This kinetic scheme could partially model the fate of ArOH in aquatic media containing natural photosensitisers, under environmental conditions.     

Simulation of dispersion of a power plant plume using an adaptive grid algorithm

A new dynamic adaptive grid algorithm has been developed for use in air quality modeling. This algorithm uses a higher order numerical scheme—the piecewise parabolic method (PPM)—for computing advective solution fields; a weight function capable of promoting grid node clustering by moving grid nodes; and a conservative interpolation equation using PPM for redistributing the solution field after movement of grid nodes. Applications of the algorithm to a model problem, in which emissions from a point source disperse through the atmosphere in time, reflect that the algorithm is able to capture not only the regional ozone plume distribution, but also the small-scale plume structure near the source. In contrast, the small-scale plume structure was not captured in the corresponding static grid solution. Performance achieved in model problem simulations indicates that the algorithm has the potential to provide accurate air quality modeling solutions at costs that may be significantly less than those incurred in obtaining equivalent static grid solutions.     

Indirect identification of an unstable intermediate in γ-HCH degradation by Pseudomonas paucimobilis UT26

While a postulated intermediate, 1,3,4,5-tetrachlorocyclohexa-1,4-diene(1,4-TCDN), is one of the key-compounds in γ-hexachlorocyclohexane(γ-HCH) degradation by Pseudomonas paucimobilis UT26(UT26), it has neither been isolated nor identified yet. Thus, an indirect identification was attempted to demonstrate its existence by a chemical oxidation of the microbial degradation products in situ.

1,2,4,5-Tetrachlorobenzene(1,2,4,5-TeCB) was newly produced by the chemical oxidation with hydrogen peroxide under iron(II) of W medium containing degradation products in situ of γ-1,3,4,5,6-pentachlorocyclohexene(γ-PCCH) by Escherichia coli MV1184(E. coli) which possesses the dehydrochlorinating activity of UT26. The chemical oxidation required the preceding enzymatic reactions to γ-PCCH by E. coli for the production of 1,2,4,5-TeCB. In facts, 1,2,4,5-TeCB was not produced from γ-PCCH by either of the chemical or enzymatic reaction alone.

Considering these results and the fact that hydrogen peroxide is known to dehydrogenate chemical compounds, it was concluded that 1,2,4,5-TeCB was produced directly from the unstable intermediate, 1,4-TCDN, on yielded from γ-PCCH by an enzymatic dehydrochlorination.     

Thyroid status in juvenile alligators (Alligator mississippiensis) from contaminated and reference sites on Lake Okeechobee,Florida, USA

Exposure to environmental contaminants has been shown to alter normal thyroid function in various wildlife species, including the American alligator (Alligator mississippiensis). Abnormalities in circulating levels of the thyroid hormone thyroxine (T₄) have been reported in juvenile alligators from several contaminated lakes in Florida. To further elucidate these functional thyroid abnormalities, this study examines the structure of thyroids and circulating T₄ concentrations from juvenile alligators collected from three sites of varying contamination on Lake Okeechobee, Florida. The following variables were used to characterize thyroid morphology: epithelial cell height, width and area, percent colloid, and follicle area. These variables were compared among study sites and between genders. No difference was detected in epithelial cell height, epithelial cell area, or follicle area among the sites, whereas significant differences in epithelial cell width (p=0.02) and percent colloid (p=0.008) were found. Animals from the most contaminated site (Belle Glade) had significantly greater epithelial cell widths and significantly less colloid present in their follicles compared to animals from the reference site (West). Gender did not have a significant interaction with site for any variable measured. Thyroxine (T₄) concentrations were elevated in the intermediately contaminated site (Conservation Area 3A) compared to the other sites (p<0.0001). It is proposed that the disruptions seen in Lake Okeechobee alligators are due to disruptions at both the thyroid and extra-thyroidal tissues.     

Ambient Air Hydroxyl Radical Concentrations: Measurements and Model Predictions

This review critically examines the limited number of recent near surface measurements of hydroxyl radicals in urban and rural areas. Predictions of average tropospheric, hemispheric and/or zonal concentrations of hydroxyl radicals derived from methylchloroform and carbon monoxide concentrations and models are discussed. The limitations of one- and two-dimensional model predictions of zonal and seasonal hydroxyl concentrations receive detailed consideration. Of especial concern on a regional scale are cloud cover effects on the penetration of ultraviolet radiation below 320 μm which is critical to the subsequent reactions forming hydroxyl radicals.     

One-dimensional simulation of solute transfer in saturated-unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous-heterogeneous, saturated-unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection-dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term. The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.