An evaluation of contaminant migration patterns at two waste disposal sites on fractured porous media in terms of the equivalent porous medium (EPM) model

Contamination has occurred many non-indurated and bedrock systems wherein the groundwater flows almost exclusively through a network of connected, open fractures. The matrix surrounding the fractures often possesses porosity which allows contaminant diffusion into the matrix. If the diffusion rates are fast relative to the fracture groundwater velocity, transport effects may be predicted by considering the system to be an equivalent porous medium (EPM). The rapidity with which fracture/immobile-matrix equilibrium is established will be determined in part by the: fracture aperture (2b); interfracture spacing (2B); porosity in the immobile matrix (θ_im); and the matrix diffusion coefficient (D′). Two systems which are characterized by very different values of the above parameters have been studied by our laboratories. At Alkali Lake, Oregon, the EPM approach describes contaminant transport well. At Bayview Park, Ontario, the EPM approach is not appropriate. Several features of the two sites are compared to illustrate the different nature of these two sites. These features include: (1) natural characteristics of the groundwater systems; (2) contaminant distributions; (3) observed transport; and (4) computed fracture/immobile-matrix diffusion times.     

Modelling of silica diffusion experiments with Si in Boom Clay

A mathematical model describing the dissolution of nuclear glass directly disposed in clay combines a first-order dissolution rate law with the diffusion of dissolved silica in clay. According to this model, the main parameters describing the long-term dissolution of the glass are ηR, the product of the diffusion accessible porosity η and the retardation factor R, and the apparent diffusion coefficient D_app of dissolved silica in clay.For determining the migration parameters needed for long-term predictions, four Through-Diffusion (T-D) experiments and one percolation test have been performed on undisturbed clay cores. In the Through-Diffusion experiments, the concentration decrease after injection of ³²Si (radioactive labelled silica) was measured in the inlet compartment. At the end of the T-D experiments, the clay cores were cut in thin slices and the activity of labelled silica in each slice was determined. The measured activity profiles for these four clay cores are well reproducible.Since no labelled silica could be detected in the outlet compartments, the Through-Diffusion experiments are fitted by two In-Diffusion models: one model assuming linear and reversible sorption equilibrium and a second model taking into account sorption kinetics. Although the kinetic model provides better fits, due to the sufficiently long duration of the experiments, both models give approximately similar values for the fit parameters. The single percolation test leads to an apparent diffusion coefficient value about two to three times lower than those of the Through-Diffusion tests.Therefore, dissolved silica appears to be strongly retarded in Boom Clay. A retardation factor R between 100 and 300 was determined. The corresponding in situ distribution coefficient K_d is in the range 25–75 cm³ g⁻¹. The apparent diffusion coefficient of dissolved silica in Boom Clay is estimated between 2×10⁻¹³ and 7×10⁻¹³ m² s⁻¹. The pore diffusion coefficient is in the range from 6×10⁻¹¹ to 1×10⁻¹⁰ m² s⁻¹.     

Phosphorus characterization in sediments impacted by septic effluent at four sites in central Canada

Characterization of phosphorus (P) enriched solids was undertaken in the sediments below four mature septic system infiltration beds, where previous monitoring of phosphate (PO₄) concentrations in the groundwater had indicated that substantial retention of P was occurring in the vadose zone. At each site, zones of sediment P enrichment were identified by an acid extraction procedure. Acid extractable sediment P concentrations were found to be 2–5 times higher than background values, within narrow discrete zones generally 10–30 cm in thickness, located within one meter of the infiltration pipes. Back scattered electron images of the P enriched zones indicated that the P solids occurred as distinct authigenic grains (up to 300 μm diameter) and as grain coatings. Microprobe analyses indicated predominantly Fe–P in calcareous sediments (Cambridge and Langton) and Al–Fe–P in non-calcareous sediments (Muskoka and Harp Lake). Porewater analyses indicated that the zones of P accumulation were closely associated with zones of redox change characterized by the conversion of effluent NH₄⁺ to NO₃⁻. The data suggests that a substantial amount of the septic derived P is being attenuated by mineral precipitation reactions that occur rapidly after the effluent encounters subsurface sediments. Reductive dissolution of ferric (oxy)hydroxide minerals as a consequence of reducing environments near the infiltrations pipes, the release of Fe²⁺ in solution and subsequent conversion of Fe²⁺ to Fe³⁺ may promote the precipitation of ferric or ferrosoferric PO₄ minerals. In sediments with limited buffering capacity (calcite deficient), the decrease in pH resulting from effluent oxidation may cause Al (oxy)hydroxide dissolution and subsequent precipitation of Al–P rich phases. These precipitation reactions appear to have the capacity to immobilize a substantial amount of septic derived P (25–99% at these sites) for a considerable period of time.     

Actinide geochemistry: From the molecular level to the real system

Geochemical processes leading to either mobilization or retention of radionuclides in an aquifer system are significantly influenced by their interaction with rock, sediment and colloid surfaces. Therefore, a sound safety assessment of nuclear waste disposal requires the elucidation and quantification of those processes. State-of-the-art analytical techniques as e.g. laser- and X-ray spectroscopy are increasingly applied to study solid–liquid interface reactions to obtain molecular level speciation insight.We have studied the sorption of trivalent lanthanides and actinides onto aluminium oxides, hydroxides and purified clay minerals by the time-resolved laser fluorescence spectroscopy and X-ray-absorption spectroscopy. Chemical constitution and structure of surface bound actinides are proposed based on spectroscopic information. Open questions still remain with regard to the exact nature of mineral surface ligands and the mineral/water interface. Similarities of spectroscopic data obtained for M(III) sorbed onto γ-alumina, and clay minerals suggest the formation of very comparable inner-sphere surface complexes such as S–O–An(III)(OH)_x^(2 − x)(H₂O)_5 − x at pH > 5. Those speciation data are found consistent with those predicted by surface complexation modelling. The applicability of data obtained for pure mineral phases to actinide sorption onto heterogeneously composed natural clay rock is examined by experiments and by geochemical modelling. Good agreement of experiment and model calculations is found for U(VI) and trivalent actinide/lanthanide sorption to natural clay rock. The agreement of spectroscopy, geochemical modelling and batch experiments with natural rock samples and purified minerals increases the reliability in model predictions.The assessment of colloid borne actinide migration observed in various laboratory and field studies calls for detailed information on actinide–colloid interaction. Kinetic stabilization of colloid bound actinides can be due to inclusion into inorganic colloid matrix or by macromolecular rearrangement in case of organic, humic/fulvic like colloids. Only a combination of spectroscopy, microscopy and classical batch sorption experiments can help to elucidate the actinide–colloid interaction mechanisms and thus contribute to the assessment of colloids for radionuclide migration.     

Adsorption of chlorpyrifos,penconazole and metalaxyl from aqueous solution by modified clays

Sorption of three pesticides (chlorpyrifos, metalaxyl and penconazole) has been measured on a commercial clay montmorillonite and on the same mineral modified with either of two cationic-surfactant micelles. Both micelle–clay complexes, commercial names Cloisite 20A and Cloisite 30B, showed a good capacity to sorb all three pesticides from water, whereas their sorption on the natural montmorillonite was not described by an isotherm. Modelling sorption on both micelle–clay complexes showed that the Freundlich sorption constant (K _F) was higher for chlorpyrifos on Cloisite 20A (K _F = 7.76) than on Cloisite 30B (K _F = 5.91), whereas the sorption of metalaxyl was stronger on Cloisite 30B (K _F = 1.07) than on Cloisite 20A (K _F = 0.57). Moreover the micelle–clay complex Cloisite 20A also showed a good affinity for penconazole, the maximum quantity adsorbed (q _m) of 6.33 mg g^?1 being 45% more than that on Cloisite 30B. Single-batch adsorption of each pesticide onto both micelle–clay complexes was studied using the Freundlich isotherm for chlorpyrifos and metalaxyl and the Langmuir isotherm for penconazole. The Cloisite 20A micelle–clay complex was predicted to require 23% less adsorbent to treat certain volumes of wastewater containing 30 mg L^?1 chlorpyrifos, 43% more to treat metalaxyl similarly and 57% less to treat penconazole compared with Cloisite 30B.     

不同粒径泥沙理化特性对磷吸附过程的影响

以北京大兴南海子湖表层沉积物为研究对象,测试分析了0.147～0.246 mm(细砂)、0.074～0.147 mm(极细砂)、0.0385～0.0740 mm(粉粒)和<0.0385 mm(粉粒粘粒混合物)4种粒径泥沙对磷的吸附行为,并采用相关分析及逐步回归分析探讨不同粒径沉积物中有机质(OM)、Fe、Al、Ca、Mn和TP含量对磷吸附过程的影响。结果表明,二级动力学方程和Langmuir模型能较好地描述南海子不同粒径泥沙的吸附动力学及等温吸附过程(R2>0.90)。粒径对单位质量泥沙吸附磷量具有明显影响,粉粒粘粒混合物>粉粒>细砂>极细砂。总体上,泥沙有机质(OM)、TP、Fe、Al、Ca和Mn含量随粒径的减小而增大,且粘粒对其影响较大。不同粒径泥沙(OM)、Fe、Al、Ca和Mn含量之间存在极显著正相关关系(P<0.01),且均对单位质量泥沙最大吸附量(Xm)和饱和吸附量(Cse)具有正效应,其中Al含量对该参数的影响更为显著。这说明泥沙对磷的吸附行为可能受到粒径和化学成分的共同影响。     

Sorption–desorption of imidacloprid onto a lacustrine Egyptian soil and its clay and humic acid fractions

Sorption–desorption of the insecticide imidacloprid 1-[(6-chloro-3-pyridinyl)-methyl]-N-nitro-2-imidazolidinimine onto a lacustrine sandy clay loam Egyptian soil and its clay and humic acid (HA) fractions was investigated in 24-h batch equilibrium experiments. Imidacloprid (IMDA) sorption–desorption isotherms onto the three sorbents were found to belong to a non-linear L-type and were best described by the Freundlich model. The value of the IMDA adsorption distribution coefficient, Kd_ads, varied according to its initial concentration and was ranged 40–84 for HA, 14–58 for clay and 1.85–4.15 for bulk soil. Freundlich sorption coefficient, Kf_ads, values were 63.0, 39.7 and 4.0 for HA, clay and bulk soil, respectively. The normalized soil Koc value for imidacloprid sorption was ～800 indicating its slight mobility in soils. Nonlinear sorption isotherms were indicated by 1/n_ads values <1 for all sorbents. Values of the hysteresis index (H) were <1, indicating the irreversibility of imidacloprid sorption process with all tested sorbents. Gibbs free energy (ΔG) values indicated a spontaneous and physicosorption process for IMDA and a more favorable sorption to HA than clay and soil. In conclusion, although the humic acid fraction showed the highest capacity and affinity for imidacloprid sorption, the clay fraction contributed to approximately 95% of soil-sorbed insecticide. Clay and humic acid fractions were found to be the major two factors controlling IMDA sorption in soils. The slight mobility of IMDA in soils and the hysteresis phenomenon associated with the irreversibility of its sorption onto, mainly, clay and organic matter of soils make its leachability unlikely to occur.     

Tritium and Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and ³⁶Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this “fast flow” in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891–989] and suggest that fast flow in fractures with minimal fracture–matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and ³⁶Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Sorption of alpha and beta hydrophobic endosulfan in a Vertisol from southeast region of Turkey

Endosulfan has been applied to control numerous insects in a variety of food and non-food crops. Limited information is available on dynamics of this pesticide in the soil. The objective of this research was to determine the adsorption–desorption behavior of the alpha (α) and beta (β) endosulfan in a Vertisol from the southeast region of Turkey, where cotton is the main crop in the large irrigated lowlands. The α and β endosulfan were adsorbed considerably and Freundlich adsorption–desorption isotherms fitted the α and β endosulfan data (R² > 0.98). Freundlich adsorption coefficients (K_f) for the α endosulfan ranged between 21.63 and 16.33 while for the β endosulfan they were between 14.01 and 17.98 for the Ap and Bw2 horizons. The difference of K_f values of α and β endosulfan for two horizons were explained with the slight difference in the amount of organic matter and clay, but considerable difference in Fe contents of the two horizons. Alpha and β endosulfan K_fd values were 118.03 and 45.81 for the Ap and 48.08 and 68.71 for the Bw2 horizons. Higher adsorption and desorption behavior of the endosulfan isomers for the same horizon was attributed to poor physical bonding between the endosulfan molecule and the surfaces of fundamental soil particles. This fact is thought to increase the effective use of endosulfan in agriculture with a possibility of its movement to the surface and groundwater in the Vertisol studied.     

Effect of sediment properties on the sorption of C12-2-LAS in marine and estuarine sediments

Linear alkylbenzene sulfonates (LAS) are anionic high production volume surfactants used in the manufacture of cleaning products. Here, we have studied the effect of the characteristics of marine and estuarine sediments on the sorption of LAS. Sorption experiments were performed with single sediment materials (pure clays and sea sand), with sediments treated to reduce their organic carbon content, and with field marine and estuarine sediments. C₁₂-2-LAS was used as a model compound. Sorption to the clays montmorillonite and kaolinite resulted in non-linear isotherms very similar for both clays. When reducing the organic content, sorption coefficients decreased proportionally to the fraction removed in fine grain sediments but this was not the case for the sandy sediment. The correlation of the sediment characteristics with the sorption coefficients at different surfactant concentrations showed that at concentrations below 10 μg C₁₂-2-LAS/L, the clay content correlated better with sorption, while the organic fraction became more significant at higher concentrations.     

An integrated assessment of seawater intrusion in a small tropical island using geophysical,geochemical, and geostatistical techniques

In this study, geophysics, geochemistry, and geostatistical techniques were integrated to assess seawater intrusion in Kapas Island due to its geological complexity and multiple contamination sources. Five resistivity profiles were measured using an electric resistivity technique. The results reveal very low resistivity <1 Ωm, suggesting either marine clay deposit or seawater intrusion or both along the majority of the resistivity images. As a result, geochemistry was further employed to verify the resistivity evidence. The Chadha and Stiff diagrams classify the island groundwater into Ca-HCO₃, Ca-Na-HCO₃, Na-HCO₃, and Na-Cl water types, with Ca-HCO₃ as the dominant. The Mg²⁺/Mg²⁺+Ca²⁺, HCO₃ ^?/anion, Cl^?/HCO₃ ^?, Na⁺/Cl^?, and SO₄ ^2?/Cl^? ratios show that some sampling sites are affected by seawater intrusion; these sampling sites fall within the same areas that show low-resistivity values. The resulting ratios and resistivity values were then used in the geographical information system (GIS) environment to create the geostatistical map of individual indicators. These maps were then overlaid to create the final map showing seawater-affected areas. The final map successfully delineates the area that is actually undergoing seawater intrusion. The proposed technique is not area specific, and hence, it can work in any place with similar completed characteristics or under the influence of multiple contaminants so as to distinguish the area that is truly affected by any targeted pollutants from the rest. This information would provide managers and policy makers with the knowledge of the current situation and will serve as a guide and standard in water research for sustainable management plan.     

Mass removal of chlorinated ethenes from rough-walled fractures using permanganate

In situ chemical oxidation (ISCO) employing permanganate is an emerging technology that has been successful at enhancing mass removal from DNAPL source zones in unconsolidated media at the pilot-scale. The focus of this study was to evaluate the applicability of flushing a permanganate solution across two single vertical fractures in a laboratory environment to remove free phase DNAPL. The fracture experiments were designed to represent a portion of a larger fractured aquifer system impacted by a near-surface DNAPL spill over a shallow fractured rock aquifer. Each fracture was characterized by hydraulic and tracer tests, and the aperture field for one of the fractures was mapped using a co-ordinate measurement machine. Following DNAPL emplacement, a series of water and permanganate flushes were performed. To support observations from the fracture experiments, a set of batch experiments was conducted. The data from both fracture experiments showed that the post-oxidation effluent concentration was not impacted by the oxidant flush; however, changes in the aperture distribution, flow field, and flow rate were observed. These changes resulted in a significant decrease to the mass loading from the fractures, and were attributed to the build-up of oxidation by-products (manganese oxides and carbon dioxide) within the fracture which was corroborated by the batch experiment data and visual examination of the walls of one fracture. These results provide insight into the potential impact that a permanganate solution and oxidation by-products can have on the aperture distribution within a fracture and on DNAPL mass transfer rates. A permanganate flush or injection completed within a fractured rock aquifer may lead to the development of an insoluble product adjacent to the DNAPL which results in the reduction or complete elimination of advective regions near the DNAPL and reduces mass transfer rates. This outcome would have significant implications on the plume generating potential of the remaining DNAPL.     

Dissecting the spatial scales of mercury accumulation in Ontario lake sediment

Total mercury concentration was analyzed in 171 lakes from pre-industrial (>30 cm depth; Hg_{pre-industrial}) and present-day sediments (0.5–1 cm; Hg_present-day). Numerous hot or cold spots of sediment mercury enrichment (Hg EF; Hg_{pre-industrial}/Hg_present-day) were evident as determined by local tests of autocorrelation, although in most cases, the maximum correlation among sites was not the nearest neighbor, indicating a strong influence of watershed characteristics. Hg EF was correlated with the area of open water (ha) (r = 0.91, p = 0.035), mine tailings (r = 0.94, p = 0.019), and organic deposits in surficial geology of the watershed (r = −0.91, p = 0.034). Through use of local rather than global regression coefficients, R² increased from 0.20 (p = 0.005) to 0.60 (p = 0.013). A broad spatial pattern (>500 km) observed only in Hg_{pre-industrial} was best explained by mean annual precipitation (shared variance = 3.5%), while finer spatial patterns only observed in Hg_present-day and Hg EF were best explained by pH (average shared variance = 10.8%).     

Sorptive interactions and binding of δ‐endotoxin protein from bacillus thuringiensis subsp. kurstaki inforest soils

Abstract

The adsorption, desorption and binding of the insecticidal protein from Bacillus thuringiensis subsp. kurstaki (Btk toxin) onto autoclaved sandy and clay loam forest soils were studied at 23°C in a buffer medium (pH 10.2) using the precipitated protein mixture (active + inactive) obtained from a commercial Btk formulation. The active protein in the buffer solution was quantified by ELISA technique. Maximum adsorption of the toxin onto the sandy (301 μg/g) and clay (474 μg/g) loam soils was found to occur after 3 and 4 hours of agitation, respectively. Adsorption of the toxin was higher in the clay loam soil than in sandy loam. Adsorption parameters were calculated using the Freundlich and linear isotherm equations. The K_F and 1/n values for the soils were 1.12 and 1.48 (sandy), and 20.42 and 0.874 (clay), respectively, indicating stronger affinity of the toxin for the clay compared to the sandy loam soil. The linear model showed deviations at higher concentrations, nevertheless using the best fit, K_D and K_OC values were computed for the two soils. For sandy loam, the K_D and K_OC values were 9.38 and 391, respectively; the corresponding values for clay loam were 13.19 and 425, confirming the higher sorption affinity of the toxin for clay loam. The adsorption data did not fit the Langmuir equation because of heterogeneity of the soil surface. Desorption studies showed that more than half of the adsorbed toxic protein remained firmly attached to sandy (162.6 μg/g or 54.5%) and clay (314.0 μg/g or 67.4%) loam soils after six 0.5‐h washes (total 3.0 h wash time). Although the toxin appears to be a non‐leacher, its lateral mobility, soil persistence and biological consequences, including bioavailability of the bound residues, are poorly understood and require further investigation.     

Kinetics of metribuzin release from bentonite-polymer composites in water

A series of bentonite polymer-composites (BPCs) loaded with metribuzin were studied for their controlled release in aqueous medium. The release of active ingredient from BPCs was significantly lower as compared to commercial metribuzin formulation. The results revealed that the cumulative metribuzin release was highest (81%) from the BPCs containing 8% clay (commercial bentonite) and 2% metribuzin which correspond to the lowest (14 days) half-life values i.e., time required for 50% release of active ingredient (t_1/2). The metribuzin release from the BPCs decreased with increased concentration of clays in polymer matrix and the release was further decreased with BPCs prepared with pure nano-bentonite. BPCs containing 12% clay and 2% metribuzin showed maximum t_1/2 values i.e., 25 and 51 days for commercial bentonite and pure nano-bentonite as clay sources, respectively. The differential behaviour in the metribuzin release rates from BPCs was ascribed due to variations in crosslinking of metribuzin in the composites. As metribuzin release was found to be slower in BPCs compared to commercial formulation, it could be used for control of weeds tailored to different crops.     

A simplified approach to modelling underground migration of radionuclides from contaminated river sediments

Accidental releases of waste water from the first Czechoslovak nuclear power plant, A1, caused contamination of sediments of the Dudváh river, flowing into the Vah river, in Slovakia. Rather high concentrations of ¹³⁷Cs and ⁹⁰Sr (2150Bq dm⁻³ and 215Bq dm⁻³, respectively) were found in bottom sediments of a former channel of the re-engineered river body at a distance of about 250 m from a village, Siladice, with water-supply wells. In order to assess the possibility of contamination of the wells, underground migration of both radionuclides from the contaminated area was simulated using an original layered convection-diffusion model. K_d values determined in laboratory experiments were used. The analysis of the hydrological situation in the area reveals that the critical condition is a dominant horizontal groundwater flow near the water table in the direction from the Váh bank to Siladice, in the periods when the contaminated body lies under the water table. The simulation calculated under conservative conditions showed that the contamination of water-supply wells would not exceed permissible concentration limits.     

Bioaccumulative characteristics of hexabromocyclododecanes in freshwater species from an electronic waste recycling area in China

Hexabromocyclododecanes (HBCDs) are now emerging ubiquitous contaminants due to their wide usage, persistence and toxicities. To investigate the bioaccumulative characteristics of HBCDs, sediments, Winkle (Littorina littorea), crucian carp (Carassius carassius) and loach (Misgurnus anguillicaudatus) were collected from two streams near an E-waste dismantling site in China, and HBCD exposure test was then conducted on Chinese rare minnow. The concentration of HBCDs was 14 ng g⁻¹ dry weight in sediments, 186, 377 and 1791 ng g⁻¹ lipid weight in winkle, crucian carp and loach, respectively. γ-HBCD was found to be the dominant diastereoisomer in the sediments (63% of total HBCDs). However, α-HBCD was selectively accumulated in the biotic samples and contributed to 77%, 63% and 63% of total HBCDs in winkle, crucian carp and loach, respectively. Moreover, an enrichment of (−)-enantiomers of α- and γ-HBCD were found in the winkle. The reverse results were observed in the crucian carp and loach. Similar observations of diastereoisomeric and enantiomeric composition were obtained in Chinese rare minnow with those found in the crucian carp and loach. These results indicate that the freshwater species from the streams are contaminated by HBCDs. α-HBCD can be selectively accumulated in organisms and the accumulative characteristics are enantioselective among species.     

Depleted N carryover, leaching and uptake for three years of irrigated corn

Deep percolation of nitrate can contribute to the deterioration of groundwater resources. Leaching of nitrate is a complex process affected by fertilizer and irrigation practices, efficiency of N use by the crop, and how the soil's water holding capacity and water transmission properties are affected by soil texture. Depleted (¹⁵NH₄)₂SO₄ fertilizer at N rates of 0, 125, 250 and 375 kg ha⁻¹ was applied annually for 3 years to continuous corn grown within three different water regimes. This time period and the labeled N permitted an evaluation of N use efficiency by the crop and NO₃ leaching and carryover on a Weld silty clay loam, a fine-textured soil, typical of the “hardland” soils of the semi-arid Great Plains. Three water regimes, W₁ ( 1.5 ET), W₂ ( ET) and W₃ ( 0.8 ET), were used. Beneath each plot within each water regime, Duke-Haise vacuum trough extractors were installed under undisturbed soil profiles at 1.22-m depth to measure weekly percolate and the NO₃ concentration in the percolate. The corn was harvested in the fall in the dent stage to measure the total above-ground biomass N uptake. Soil profiles (1.8 m) were sampled annually in the fall after crop harvest to determine NO₃---N in the soil or carryover.Great variability was encountered in measuring the amount of extractor water and its NO₃ content under each water regime, which made estimates of N03 leaching losses unreliable. Also, the variability demonstrates formidable problems in quantifying percolation losses with vacuum trough extractors under undisturbed fine-textured soil profiles. With the highest N rate of 376 kg ha⁻¹ yr⁻¹ and within the water regime W₁, where leaching was expected to be greatest, only 1% of the cumulative labeled N applied was found in extractor waters and most movement of the labeled N into extractors occurred the third year. The 125-kg-ha⁻¹ yr⁻¹ fertilizer N rate significantly increased the crop yield over the unfertilized plots without increasing residual NO₃---N accumulation; whereas fertilizer N rates of > 125 kg ha⁻¹ yr⁻¹ did not appreciably increase plant yields over the 125-kg-ha⁻¹-N rate, but did appreciably increase residual NO₃.     

Mathematical prediction of imidacloprid persistence in two Croatian soils with different texture,organic matter content and acidity under laboratory conditions

In the present laboratory study, persistence of imidacloprid (IMI) as a function of initial insecticide concentration and soil properties in two Croatian soils (Krk sandy clay and Istria clay soils) was studied and described mathematically. Upon fitting the obtained experimental data for the higher concentration level (5 mg/kg) to mathematical models, statistical parameters (R ², scaled root mean squared error and χ ² error) indicated that the single first-order kinetics model provided the best prediction of IMI degradation in the Krk sandy clay soil, while in the Istria clay soil biphasic degradation was observed. At the lower concentration level (0.5 mg/kg), the biphasic models Gustafson and Holden models as well as the first-order double exponential model fitted the best experimental data in both soils. The disappearance time (DT₅₀) values estimated by the single first-order double exponential model (from 50 to 132 days) proved that IMI can be categorized as a moderately persistent pesticide. In the Krk sandy clay soil, resulting DT₅₀ values tended to increase with an increase of initial IMI concentration, while in the Istria clay soil, IMI persistence did not depend on the concentration. Organic matter of both experimental soils provided an accelerating effect on the degradation rate. The logistic model demonstrated that the effect of microbial activity was not the most important parameter for the biodegradation of IMI in the Istria clay soil, where IMI degradation could be dominated by chemical processes, such as chemical hydrolysis. The results pointed that mathematical modeling could be considered as the most convenient tool for predicting IMI persistence and contributes to the establishment of adequate monitoring of IMI residues in contaminated soil. Furthermore, IMI usage should be strictly controlled, especially in soils with low organic matter content where the risk of soil and groundwater contamination is much higher due to its longer persistence and consequent leaching and/or moving from soil surface prior to its degradation.