顺序注射平台-分光光度法测定水中总磷

参照希思迪在线总磷分析仪和荷兰Skalar在线过程分析仪测定水中总磷的方法,对顺序注射平台分析技术进行改进,从而建立了钼酸铵分光光度法测定水中正磷酸盐的快速分析方法,该方法集微功耗、微量试剂消耗于一体。注射平台具有模块化、投入式探头等结构特征,仪器包含3个检测量程,适用于环境监测等领域的实时在线分析。实验选择磷酸蓝的测定吸收波长为880 nm,当总磷质量浓度为0.10~2.00 mg/L时,浓度与吸光度呈线性关系,线性回归方程y=0.970 3x-0.025 3,相关系数为0.999 5,方法的检出限为1.36×10-2mg/L。     

The impact of an industrial complex on freshly deposited sediments, Chener Rahdar river case study, Shiraz, Iran

Concentrations of elements (As, Co, Cu, Ni, Mo, Pb, V, and Zn) are studied in the sediments of two adjacent stretches of Chenar Rahdar river. The first stretch (S1) is influenced by urban and arable land wastewater, and the second (S2) is mainly loaded with industrial effluents. The average abundance order of heavy metals content in S1 sediments is Ni > V > Zn > Cu > Co > As > Pb > Mo and in S2 sediments is Ni > Zn > V > Cu > Mo > Pb > Co > As. The maximum average concentration for these heavy metals (except for As) occurs in the S2 sediments. The contamination factor (CF) base of background in S1 for eight analyzed elements is moderate. The CF for Cu, Zn, and Pb in S2 sediments is considerable. The highest CF in S1 and S2 sediments is observed for Mo (CF = 10.95 and 12.41) and indicates very high contamination. The application of modified degree of contamination values (mC_d) indicates low and high degree of contamination (1.89–4.15) in S1 and S2, respectively. Calculated enrichment factors (EF) reveal enrichment of Mo and As in S1 and Zn, Cu, Mo, and Pb in S2 compared to the average abundances of background level. The maximum EF for Mo is 7.61 (significant enrichment), while Pb, Zn, and Cu with maximum EF between 2 and 5 indicate moderate contamination. Principal component analysis (PCA) shows distinctly different elemental associations in S1 and S2 sediments. The strong association of Zn, Co, Ni, Sc, Cu, Al and Fe in S1 suggests a similar source. The results of PCA for Zn, Pb, Mo and Cu in S2 (componente2) indicate that these metals are influenced by anthropogenic activity. Also, high loading heavy metals with OC (0.97) indicate that organic carbon plays a significant role in the distribution and sorption of these heavy metals in the sediments. Factor analysis indicates that As and Mo behave differently in sediment samples.     

Modelling Dissolved Organic Carbon Turnover in Humic Lake Örträsket, Sweden

The organic carbon balance of a lake with high input of allochthonous organic carbon is modelled integrating physical, chemical and biological processes. The physical model captures the behaviour of real thermal stratification in the lake for different flow situations during the period 1993–1997. The dissolved organic carbon model is based on simulated trajectories of water parcels. By tracking parcels, account is kept of environmental factors such as temperature and radiation as well as DOC quality for each parcel. The DOC concentration shows seasonal variations primarily dependent on inflow. The organic matter degradation (bacterial- and photodegradation) in the lake amounts to 1.5–2.5 mg C l^–1 yr^–1, where photooxidation is responsible for approximately 10%. The estimated DIC production in the lake is large compared to sediment mineralisation and primary production. The main conclusion is that the model with the selected parameterisations of the degradation processes reasonably well describes the DOC dynamics in a forest lake.     

离子色谱法测定丙烯酸丁酯生产废水中丙烯酸和对甲基苯磺酸

建立了用离子色谱法测定丙烯酸丁酯生产废水中丙烯酸和对甲基苯磺酸的新方法。采用IonPac AS11-HC阴离子色谱柱,淋洗液为30mmol/L的NaOH淋洗液,流速1ml/min,抑制电流75mA。丙烯酸根的检测下限为0.44mg/L(峰面积)、0.46mg/L(峰高),保留时间、峰高、峰面积的变异系数分别为0.03%、0.50%、0.50%,峰高、峰面积的线性回归系数R分别为0.9929、0.9954。对甲基苯磺酸根的检测下限为0.49mg/L(峰面积)、0.48mg/L(峰高),保留时间、峰高、峰面积的变异系数分别为0.02%、0.48%、0.47%,峰高、峰面积的线性回归系数R分别为0.9990、0.9998。常规无机阴离子对丙烯酸、对甲基苯磺酸的测定没有干扰。对实际废水样品进行分析,丙烯酸、对甲基苯磺酸的回收率分别为101.6%～104.4%、95.4%～96.8%。方法简便、灵敏、可靠。     

Inverse Modeling for Aquatic Source and Transport Parameters Identification and its Application to Fukushima Nuclear Accident

Inverse modeling technique based on nonlinear least square regression method (LSRM) is developed for the identification of aquatic source and transport parameters. Instantaneous line source release model in two-dimensional domain and continuous point source release model in three-dimensional domain are used for the purpose. Case studies have been carried out for both types of releases to illustrate their application. Error analysis has been carried out to identify the maximum error that can be tolerated in the input concentration data used in the inverse model and to specify the minimum number of sampling points to generate such input data. The LSRM is compared with the well-established correlation coefficient optimization method for instantaneous line source release model, and good comparison is observed between them. The LSRM is used to quantitatively estimate the releases of different radionuclides into the Pacific Ocean which has resulted due to the discharge of highly radioactive liquid effluent from the affected Daiichi Nuclear Power Station at Fukushima in Japan. The measured concentrations of these radionuclides in seawater samples collected from two sampling points near Fukushima are used for the estimation. The average release works out to be 1.09?×?10¹⁶ for ¹³¹I, 3.4?×?10¹⁵ Bq for ¹³⁴Cs, and 3.57?×?10¹⁵ Bq for ¹³⁷Cs. Very good agreement is observed between the releases estimated in this study and those estimated by other different agencies.     

Modelling Temperature Trends in New Zealand

If global warming is accelerating, then one might expect temperatures for most stations to be accelerating and perhaps variability to be increasing. In this study, we examine 57 New Zealand temperature time series for evidence of non-linearity and changing variability. These correspond to time series for annual minima, annual means and annual maxima for 19 stations. Estimation is by an extended least-squares method. We find a surprising diversity of behaviour of these series – presumably reflecting their different geographic factors as well as series length. We give evidence of regions where temperatures are decreasing. For series where a linear trend is significant, it is downwards in about one third of the cases. This proportion was higher in the South Island, especially for series of minima. Where a non-linear trend is significant, temperatures are decelerating in about one half of the cases. The ratio of downward to upward trends is highest among annual maxima and South Island minima and smallest in annual means. Where a linear trend in the variability is significant, it is decreasing in 13 cases and increasing in 5 cases, although possibly this is partly due to poorer quality data last century. Where a non-linear trend in the variability is significant, variability is decelerating in about two thirds of the cases. The results are used to project upper and lower return levels of minima, means and maxima for each of the series to the year 2010.     

Comparative rice seed toxicity tests using filter paper,growth pouch-tm,and seed tray methods

Paper substrate, especially circular filter paper placed inside a Petri dish, has long been used for the plant seed toxicity test (PSTT). Although this method is simple and inexpensive, recent evidence indicates that it gives results that are significantly different from those obtained using a method that does not involve paper, especially when testing metal cations. The study compared PSTT using three methods: filter paper, Growth Pouch-TM, and seed tray. The Growth Pouch-TM is a commercially available device. The seed tray is a newly designed plastic receptacle placed inside a Petri dish. The results of the Growth Pouch-TM method showed no toxic effects on rice for Ag up to 40 mg L^–1 and Cd up to 20 mg L^–1. Using the seed tray method, IC₅₀ (50% inhibitory effect concentration) values were 0.55 and 1.4 mg L^–1 for Ag and Cd, respectively. Although results of filter paper and seed tray methods were nearly identical for NaF, Cr(VI), and phenol, the toxicities of cations Ag and Cd were reduced by using the filter paper method; IC₅₀ values were 22 and 18 mg L^–1, respectively. The results clearly indicate that paper substrate is not advisable for PSTT.     

Desert Water Harvesting to Benefit Wildlife: A Simple, Cheap, And Durable Sub-Surface Water Harvester for Remote Locations

A sub-surface desert water harvester was constructed in the sagebrush steppe habitat of south-central Idaho, U.S.A. The desert water harvester utilizes a buried micro-catchment and three buried storage tanks to augment water for wildlife during the dry season. In this region, mean annual precipitation (MAP) ranges between about 150–250 mm (6″–10″), 70% of which falls during the cold season, November to May. Mid-summer through early autumn, June through October, is the dry portion of the year. During this period, the sub-surface water harvester provides supplemental water for wildlife for 30–90 days, depending upon the precipitation that year. The desert water harvester is constructed with commonly available, “over the counter” materials. The micro-catchment is made of a square-shaped, 20 mL. “PERMALON” polyethylene pond liner (approximately 22.9 m × 22.9 m = 523 m²) buried at a depth of about 60 cm. A PVC pipe connects the harvester with two storage tanks and a drinking trough. The total capacity of the water harvester is about 4777 L (1262 U.S. gallons) which includes three underground storage tanks, a trough and pipes. The drinking trough is refined with an access ramp for birds and small animals. The technology is simple, cheap, and durable and can be adapted to other uses, e.g. drip irrigation, short-term water for small livestock, poultry farming etc. The desert water harvester can be used to concentrate and collect water from precipitation and run-off in semi-arid and arid regions. Water harvested in such a relatively small area will not impact the ground water table but it should help to grow small areas of crops or vegetables to aid villagers in self-sufficiency.     

Some selected heavy metal concentrations in water,sediment, and oysters in the Er-Ren estuary,Taiwan: chemical fractions and the implications for biomonitoring

Studies of heavy metal contamination and ecological risk in estuaries are an important emerging area of environmental science. However, there have been few detailed studies of heavy metal contamination that concern the spatial variation of heavy metal levels in water, sediment, and oyster tissue. Because of the effective uptake of heavy metals, cultured oysters are a cheap and effective subject for study. This study, conducts an experiment in the Er-Ren river to examine the biological uptake of heavy metals in farmed, cultured oysters. The distribution of copper, zinc, lead, cadmium, and arsenic concentrations in water, sediment, and oysters from the Er-Ren river is also evaluated. By sequential extraction of the sediments, the following order of mobilities is found for heavy metals Pb?>?Cd?>?As?>?Zn?>?Cu. The highest percentages of heavy metals are found in the residual phase. The mean uptake rates for young oysters are 7.24 mg kg^?1 day^?1 for Cu and 94.52 mg kg^?1 day^?1 for Zn, but that for adult oyster is 10.79 mg kg^?1 day^?1 for Cu and 137.24 mg kg^?1 day^?1 for Zn. With good policies and management, the establishment of cultured oyster frames in these contaminated tributaries and near shore environments is a potential method for removing Cu and Zn and protecting the coast.     

Spectrophotometric determination of platinum(IV) in alloys, complexes, environmental, and pharmaceutical samples using 4-[N,N-(diethyl)amino] benzaldehyde thiosemicarbazone

4-[N,N-(Diethyl)amino] benzaldehyde thiosemicarbazone (DEABT) is proposed as an analytical reagent for the spectrophotometric determination of platinum(IV). The DEABT forms 1:2 yellow complex with Pt(IV), which is sparingly soluble in water and completely soluble in water–ethanol–DMF medium. The Pt(IV)–DEABT complex shows maximum absorbance at 405 nm. Beer’s law is valid up to 7.80 μg cm^???3, and optimum concentration range for the determination of platinum(IV) is 0.48–7.02 μg cm^???3. The molar absorptivity and Sandell’s sensitivity of the method are found to be 1.755 × 10⁴ dm³ mol^???1 cm^???1 and 0.0012 μg cm^???2, respectively. The relative error and coefficient of variation (n?=?6) for the method does not exceed ±0.43% and 0.35%, respectively. Since the method tolerates a number of metal ions commonly associated with platinum, it can be employed for the determination of platinum in environmental samples, pharmaceutical samples, alloys, catalysts, and complexes. The method is rapid as the Pt(IV)–DEABT complex is soluble in water–ethanol–DMF medium and not requiring any time consuming extraction method for the complex.     

京津冀地区经济发展与PM10浓度变化的灰色关联分析

近几年,京津冀地区因灰霾天气频发而备受关注。通过分析京津冀地区PM10浓度变化与主要经济指标的关系,查找京津冀地区空气污染原因,并为今后污染控制措施提供参考。采用灰色关联分析法,对2007—2012年京津冀地区PM10与GDP、钢铁产量、煤炭消费总量、原油消费总量、天然气消费总量5项经济指标进行了分析,结果表明,在与本地经济指标的关联度方面,京津冀地区、北京市、天津市的PM10与煤炭消费总量相关性高;河北省略有不同,钢铁产量为关联度最高的经济指标。在与京津冀地区经济指标的关联度方面,相比本地,北京市、河北省均呈现较大差异性,天津市差异性较小,说明北京市、河北省空气质量受自身影响较大,天津市空气质量受区域传输影响明显。建议整个京津冀地区均应重点加大"控煤"力度,天津市还需加强区域联防联控力度,河北省应加快钢铁企业转型升级,压缩过剩产能。     

Estimation of the Thermal Sensitivity of a Dwelling to Variations in Architectural and Behavioural Parameters

The provision of energy for households is a major contributor to greenhouse gas emissions in the domestic sector. There is significant scope for energy savings and thus emission reduction in this sector. This paper constructs a bridge between thermal analysis and environmental assessment in the occupancy stage of the dwelling life cycle, approaching both methods as being on a common system Environment–Dwelling–Inhabitant. The importance of the local level in the thermal analysis and assessment of either a dwelling or an inhabitant's behaviour is demonstrated. It is shown that the researcher can choose between natural experiment, computer experiment and mathematical modelling to analyze the system. Such a choice is discussed for two particular methods of thermal analysis. Thermal Nomograms are the result of mathematical modelling, and the Energy Score Sheet is the result of a computer experiment in energy rating. Despite being developed for the purpose of thermal analysis, these methods can also be utilized for environmental assessment at the local level. The discussion centres on Australian conditions where the majority of the population resides in a relatively benign climate wherein behavioural alterations have significant potential for energy savings and environmental impact. The purpose of the paper is to emphasise benefits apart from energy minimisation to promote the use of energy efficient housing strategies.     

Modeling the Effect of Liming on Calcium Concentration in Swedish Lakes

To abate the severe acidification of Swedish lakes are about 7% directly limed (AL-lakes). These are important to include when assessing acidification. Consequently, we need to estimate their chemical status as if they were not lime-affected. We hypothesize that there is a spatial dependence for calcium (Ca) and magnesium (Mg) concentrations. Spatial variation is determined using variograms. We study lakes within two 150 ∗ 150 km quadrates (EMEP150-grid) in the middle and southern parts of Sweden. We model Ca in AL-lakes using the Ca and Mg concentrations in nearby, unaffected lakes (UL-lakes). The mean Ca/Mg for the three or seven closest UL-lakes is used. We assume that Ca/Mg is constant for UL-lakes and that Mg unaltered by liming. For UL-lakes, Ca and Mg are spatially dependent. For AL-lakes, there is a spatial dependence for Mg, but, as anticipated, not for Ca. The modeled Ca for UL-lakes show a good correspondence with measured Ca; r ² > 0.6; slopes close to 1 and intercepts almost 0. Using three or seven nearby lakes when modeling Ca give similar results. Modeled Ca was higher than measured Ca for 36% (three closest lakes) and 38% (seven closest) of the modeled AL-lakes in the middle part of Sweden. The corresponding results for the southern part were 10% and 9%, respectively. This model is an acceptable estimator of Ca in limed lakes and for estimating critical loads and acidification status on an EMEP-scale.     

Measures of flow variability for great lakes tributaries

Design of monitoring programs for load estimation is often hampered by the lack of existing chemical data from which to determine patterns of flux variance, which determine the sampling program requirements when loads are to be calculated using flux-dependent models like the Beale Ratio Estimator. In contrast, detailed flow data are generally available for the important tributaries. For pollutants from non-point sources there is often a correlation between flow and pollutant flux. Thus, measures of flow variability might be calibrated to flux variability for well-known watersheds, after which flow variability could be used as a proxy for flux variability to estimate sampling needs for tributaries for which adequate chemical observations are lacking.Three types of measures of flow variability were explored: ratio measures, which are of the form q _x/q_y, where q _xis the flow corresponding to the percentile x, and y=100–x; spread measures, of the form (q _x–q_y)/q_m, where q _mis the median flow; and the coefficient of variation of the logs of flows. In the latter, flows are log transformed because flow distributions are often approximately log-normal. Three ratio measures were evaluated, based on the percentiles (10,90), (20,80), and (25,75). The analogous spread measures were also evaluated; the spread measure based on percentiles (25,75) is derived from the commonly used fourth spread of non-parametric statistics. The ratio measures and the spread measures are scale independent, and thus are measures only of the shape of the distribution. The coefficient of variation is also scale independent, but in log space.Values of these measures of flow variability for 120 Great Lakes tributaries are highly intercorrelated, although the relationship is often non-linear. The coefficient of variation of the log of the flows is also well correlated with the coefficient of variation of fluxes of suspended solids, total phosphorus, and chloride, for a smaller set of rivers where the existence of abundant chemical data allows comparison.Tributaries with abnormal distributions often show up as outliers when one measure of flow variability is plotted against another. Several examples are discussed.     

Biofilms coating Schoenoplectus californicus as indicators of water quality in the Río de la Plata Estuary (Argentina)

N-ethyl-3-carbazolecarboxaldehyde-3-thio- semicarbazone (ECCT) as an new analytical reagent used for the development of a highly sensitive extractive spectrophotometric method for the determination of copper(II). The ECCT forms a greenish-yellow colored 1:1 (M:L) complex with copper(II) at pH 3.0, which is well extracted into n-butanol and shows maximum absorbance at 380nm. The color of the complex is stable for more than forty eight hours. The system obey Beer's law in the range 0.4–3.6 with 2.243 × 10⁴lmol⁻¹cm⁻¹, 2.83 × 10^-3μgcm⁻² molar absorptivity and Sandell's sensitivity respectively. The regression coefficient is 0.412 with 0.99 correlation coefficient. The precision and accuracy of the method was checked by finding the relative standard deviation (0.422%). This developed method has been successfully employed for the determination of copper(II) in environmental and pharmaceutical samples. The method is evaluated by analyzing samples from the bureau of analyzed samples (BCS 233, 266, 216/1, 207 and 179) and by inter comparison of experimental values using AAS.     

Expres: An expert system for assessing the fate of pesticides in the subsurface

The expert system described here is designed to aid regulatory personnel in their assessment of the potential for pesticides to contaminate the soil and groundwater environment. The expert system, known as EXPRES (EXpert system for Pesticide Regulatory Evaluation Simulations), consists of existing numerical models which are used to simulate the transport and transformation of pesticides in the unsaturated zone, coupled with a knowledge-based system that guides the user through the choice of all the necessary information for characterizing the physical, meteorological, hydrogeological, pedological and agricultural settings of typical agricultural regions across Canada, as required by the pesticide model. The expert system is designed to be used as a management tool to aid in policy decisions. Thus, its purpose is to provide an assessment of the potential hazards and to identify if further study is warranted.     

Calibration and verification of a dissolved oxygen management model for a highly polluted river with extreme flow variations in Pakistan

A dissolved oxygen (DO) model is calibrated and verified for a highly polluted River Ravi with large flow variations. The model calibration is done under medium flow conditions (431.5 m³/s), whereas the model verification is done using the data collected during low flow conditions (52.6 m³/s). Biokinetic rate coefficients for carbonaceous biochemical oxygen demand (CBOD) and nitrogenous biochemical oxygen demand (NBOD) (i.e, K _cr and K _n ) are determined through the measured CBOD and ammonia river profiles. The calculated values of K _cr and K _n are 0.36 day^?1 and 0.34 day^?1, respectively. The close agreement between the DO model results and the field values shows that the verified model can be used to develop DO management strategies for the River Ravi. The biokinetic coefficients are known to vary with degree of treatment (DOT) and therefore need to be adjusted for a rational water quality management model. The effect of this variation on level of treatment has been evaluated by using the verified model to attain a DO standard of 4 mg/L in the river using the biokinetic rate coefficients as determined during the model calibration and verification process. The required DOT in this case is found to be 96 %, whereas the DOT is 86 % if adjusted biokinetic rate coefficients are used to reflect the effect of wastewater treatment. The cost of wastewater treatment is known to increase exponentially as the removal efficiency increases; therefore, the use of appropriate biokinetic coefficients to manage the water quality in rivers is important.     

Chemical characterization of precipitation at a coastal site in Trombay, Mumbai, India

Precipitation is the best scavenger for the particulates and dissolved gaseous pollutants present in the atmosphere. The chemical composition of precipitation is dominated by a number of in-cloud and below-cloud scavenging processes. The present study is aimed at analyzing the chemical composition of rainwater in the relatively less industrialized part of Mumbai. The pH of rainwater in this region ranges from 4.8 to 6.4. The percentage contributions of ions were calculated and the major contributing ions were calcium (28%), chloride (23%), sodium (18%), sulfate (14%), magnesium (11%), ammonium (4%), potassium (1%), and nitrate (1%). The correlation coefficient is highest for Na and Cl (r ²?=?0.99), giving a clear indication of contribution from sea salt. Sulfate and nitrate ions also show a very good correlation (r ²?=?0.90), which may be due to their coemission from fossil fuel combustion. Acidification caused by these ions is neutralized by Ca, Mg, and NH₄ ions. The neutralization effect due to these ions is validated by calculating the neutralization factor (NF). The NF values are in the order Ca > Mg > NH₄. The major source contributors for the ions in precipitation are sea salt (Na, Cl, and K) and fossil fuel combustion (SO₄ and NO₃). These assumptions are supported by the values of wet-only ratio, enrichment factor, and percent sea salt fraction.     

Validation and Applications of a Chemical Equilibrium Model for Struvite Precipitation

Removal of nutrients like nitrogen and phosphorus from wastewater can be accomplished by precipitating these as the mineral struvite (NH₄MgPO₄. 6H₂O). Predicting struvite precipitation potential, yield, and purity is important for designers and operators of reactors for struvite precipitation. In this paper, a mathematical model of this precipitation process is developed using physicochemical equilibrium expressions, mass balance equations for nitrogen, phosphorous and magnesium, and charge balance. The model was simulated to explicitly solve for equilibrium concentrations of eighteen species that included dissolved (three), ionic (ten), and solid (five) species for a given set of initial concentrations of ammonium–nitrogen, magnesium and phosphate–phosphorus, and pH. The model simulations were validated against literature experimental data, which used synthetic as well as actual wastewater, and data from our experiments. The model satisfactorily predicted most data. Struvite fraction in the precipitate ranged from 27% to 100%. The purity of struvite in the precipitate and the pH that maximizes struvite fraction was dependent on the initial concentrations of ammonium, magnesium, and phosphate. Optimum pH and struvite fraction was, respectively, 8.5 and 29.3% for an equimolar mixture of ammonia, magnesium, and phosphate and 9.8 and 98.3% for 10:1.7:3.4 mM ratio. Struvite fraction in the precipitate increased as magnesium became limiting or as ammonia to phosphate ratio increased and magnesium to phosphate ratio decreased. Since the struvite component is only a fraction of the total solids, it is erroneous to report the total precipitate produced as being struvite as is conventionally done.     

Factors affecting metribuzin retention in Algerian soils and assessment of the risks of contamination

Metribuzin is a widely used herbicide around the world but it could lead to soil and water contamination. Metribuzin retention on a silty–clay agricultural soil of Algeria was studied in laboratory batch experiments to assess the contamination risk of the groundwater. Factors conditioning the fate of metribuzin were investigated: soil nature, metribuzin formulation, NPK fertilizer, and soil pH. Freundlich sorption isotherms gave the coefficients K _F between 1.2 and 4.9 and 1/n _a between 0.52 and 0.93. The adsorption is directly dependent on organic and clay soil contents. Formulated metribuzin (Metriphar) reduces the adsorption (K _F?=?1.25) compared to pure metribuzin (K _F?=?2.81). The addition of an NPK fertilizer decreases the soil pH (6.67 for the soil without fertilizer and 5.86 for 2 % of fertilizer) and increases metribuzin adsorption (K _F is 4.83 for 2 % of fertilizer). The pH effect on the adsorption is corroborated in experiments changing the soil pH between 5 (K _F is 4.17) and 8 (K _F is 1.57) under controlled conditions. Desorption isotherms show a hysteresis and only 30 to 40 % of the initially adsorbed metribuzin is released. The estimated GUS index is ≥2.8 for a DT50?≥?30 days. K _F values and the hysteresis show that metribuzin is little but strongly retained on the soil. Formulated metribuzin and addition of fertilizer affect the retention. However, the GUS index indicates a high mobility and a significant risk of leaching. The most appropriate risk management measure would be an important increase in organic matter content of the soil by addition of organic amendments.