Toxicities of ten metals to Crassostrea gigas and Mytilus edulis embryos and Cancer magister larvae

Pacific oyster (Crassostrea gigas Thunberg) embryos, bay mussel (Mytilus edulis Linnaeus) embryos, and Dungeness crab (Cancer magister Dana) zoea I stage larvae were exposed to arsenic, cadmium, chromium, copper, lead, mercury, nickel, selenium, silver, and zinc salts. The two bivalve species were exposed through the completion of embryogenesis (48 h), and the effects of the metals upon abnormal development were monitored. The crab larval stages were exposed to the metal solutions for 96 h. Concentrations (μg l_.^?1 as metal) which caused 50% abnormal development in the two species of bivalve larvae and which were fatal to 50% of the crab larvae ranged from 5.3 (oyster) to 49 (crab) for copper; 5.8 (mussel) to 8.2 (crab) for mercury; 14 (mussel) to 55 (crab) for silver; 119 (oyster) to 456 (crab) for zinc; 232 (crab) to 3000 (mussel) for arsenic; 247 (crab) to 1200 (mussel) for cadmium; 349 (oyster) to 4360 (crab) for nickel; 476 (mussel) to 758 (oyster) for lead; 3440 (crab) to 4538 (oyster) for hexavalent chromium; and selenium greater than 10 000 for the three species. The acute toxicity testing of C. gigas embryo together with two California native species confirmed its usefulness for predicting toxicity to the latter species.     

A note on topological matrices and characteristic polynomials of fourth-magnitude drainage net maps

I suggest that a connection matrix and a characteristic polynomial, which can be formed for any drainage net, using Shreve link magnitudes, may be used in studying drainage nets. I show how the matrix and the polynomial can be formed, using selected drainage nets as examples.     

Dynamical geochemistry of the Hawaiian plume

This paper presents a simple dynamical model for melting and trace element distribution in the Hawaiian mantle plume. I model the plume as a partially molten stagnation point flow against the oceanic lithosphere, and obtain solutions for the temperature, melt migration rate, and trace element concentration within it. Trace element concentrations in the melt exceed simple batch melting predictions by up to 70%. The magnitude of this effect depends strongly on the solid-melt partition coefficientK. Trace elements with differentK therefore experience a “dynamical fractionation” within the plume, and incompatible trace element ratios such asLa/Ce always exceed the batch melting predictions. I suggest a simple model for plume-lithosphere interaction in which melts from these two sources mix in proportions determined by thermodynamic constraints. The model can explain the composition of basalts from Haleakala if the degree of melting of the lithosphereF₁ decreases with time from roughly 10% for tholeiites to 2% for alkalic basalts. These values are considerably higher than previous estimates ofF₁ < 1%, and imply correspondingly smaller and more realistic values ( 10 km) for the thickness of the melted part of the lithosphere. Partial melting of additional depleted sources such as the asthenosphere is therefore not required by the Haleakala data. Estimates ofF₁ are highly sensitive to the values chosen for the partition coefficients, however, and should therefore be interpreted with caution.     

In Situ Measurement and Numerical Simulation of Oxygen Limited Biotransformation

Enhanced subsurface biorestoration is rapidly becoming recognized as a valuable tool for the restoration of hydrocarbon-contaminated aquifers and sediments. Previous field and laboratory studies at a former wood creosoting facility near Conroe, Texas, have indicated that insufficient oxygen is the primary factor limiting the biotransformation of polynuclear aromatics (PNAs) in sediments and ground water at this site. A series of laboratory experiments and field push-pull injection tests were performed as part of this project to: (1) study the effect of low oxygen concentrations on the biotransformation of PNAs; (2) identify the minimum concentration of PNAs that could be achieved through the addition of oxygen alone; (3) confirm that enhanced subsurface biorestoration is feasible at this site; and (4) test an existing numerical model of the biotransformation process (BIOPLUME). The laboratory studies demonstrated that biotransformation of the PNAs was not inhibited at dissolved oxygen concentrations as low as 0.7 mg/L although this work did suggest that there may be a minimum PNA concentration of 30 to 70 μg/L total PNAs below which biotransformation was inhibited. The field push-pull tests confirmed that addition of oxygen was effective in enhancing the subsurface biodegradation of the PNAs. The minimum concentration achieved using oxygen alone was approximately 60 μg/L total PNAs. Minimal biotransformation of these compounds was observed without oxygen addition. The numerical model BIOPLUME was tested against monitoring data from the field experiments and appears to provide a good approximation of the biodegradation process.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 2. coupled model of water and salt balances

A long-term salt balance model is coupled with the small catchment water balance model presented in Part 1 of this series of papers. The salt balance model was designed as a simple robust, conceptually based model of the fundamental salt fluxes and stores in forested and cleared catchments. The model has four interdependent stores representing salt storage in the unsaturated zone, the deep permanent saturated groundwater system, the near-stream perched groundwater system and in a ‘salt bulge’ just above the permanent water-table. The model has performed well in simulations carried out on Salmon and Wights, two small experimental catchments in south-west Western Australia. When applied to Wights catchment the salt balance model was able to predict the stream salinities prior to clearing of native forests, and the increased salinities after the clearing.     

Upscaling unsaturated hydraulic parameters for flow through heterogeneous anisotropic sediments

We compare two methods for determining the upscaled water characteristics and saturation-dependent anisotropy in unsaturated hydraulic conductivity from a field-scale injection test. In both approaches an effective medium approximation is used to reduce a porous medium of M textures to an equivalent homogenous medium. The first approach is a phenomenological approach based on homogenization and assumes that moisture-based Richards’ equation can be treated like the convective–dispersive equation (CDE). The gravity term, dK_z(θ)/d(θ), analogous to the vertical convective velocity in the CDE, is determined from the temporal evolution of the plume centroid along the vertical coordinate allowing calculation of an upscaled K_z(θ). As with the dispersion tensor in the CDE, the rate of change of the second spatial moment in 3D space is used to calculate the water diffusivity tensor, D(θ), from which an upscaled K(θ) is calculated. The second approach uses the combined parameter scale inverse technique (CPSIT). Parameter scaling is used first to reduce the number of parameters to be estimated by a factor M. Upscaled parameters are then optimized by inverse modeling to produce an upscaled K(θ) characterized by a pore tortuosity–connectivity tensor, L. Parameters for individual textures are finally determined from the optimized parameters by inverse scaling using scale factors determined a priori. Both methods produced upscaled K(θ) that showed evidence of saturation dependent anisotropy. Flow predictions with the STOMP simulator, parameterized with upscaled parameters, were compared with field observations. Predictions based on the homogenization method were able to capture the mean plume behavior but could not reproduce the asymmetry caused by heterogeneity and lateral spreading. The CPSIT method captured the effects of heterogeneity and anisotropy and reduced the mean squared residual by nearly 90% compared to local-scale and upscaled parameters from the homogenization method. The Pacific Northwest National Laboratory is operated for the US Department of Energy by Battelle under Contract DE-AC05-76RL01830.     

Assessing the performance of a physically based hydrological model using a proxy‐catchment approach in an agricultural environment

Physically based models are useful frameworks for testing intervention strategies designed to reduce elevated sediment loads in agricultural catchments. Evaluating the success of these strategies depends on model accuracy, generally established by a calibration and evaluation process. In this contribution, the physically based SHETRAN model was assessed in two similar U.K. agricultural catchments. The model was calibrated on the Blackwater catchment (18 km²) and evaluated in the adjacent Kit Brook catchment (22 km²) using 4 years of 15 min discharge and suspended sediment flux data. Model sensitivity to changes in single and multiple combinations of parameters and sensitivity to changes in digital elevation model resolution were assessed. Model flow performance was reasonably accurate with a Nash–Sutcliffe efficiency coefficient of 0.78 in Blackwater and 0.60 in Kit Brook. In terms of event prediction, the mean of the absolute percentage of difference (μAbs_diff) between measured and simulated flow volume (Qv), peak discharge (Qp), sediment yield (Sy), and peak sediment flux (Sp) showed larger values in Kit Brook (48% [Qv], 66% [Qp], 298% [Sy], and 438% [Sp]) compared with the Blackwater catchment (30% [Qv], 41% [Qp], 106% [Sy], and 86% [Sp]). Results indicate that SHETRAN can produce reasonable flow prediction but performs less well in estimation of sediment flux, despite reasonably similar hydrosedimentary behaviour between catchments. The sensitivity index showed flow volume sensitive to saturated hydraulic conductivity and peak discharge to the Strickler coefficient; sediment yield was sensitive to the overland flow erodibility coefficient and peak sediment flux to raindrop/leaf soil erodibility coefficient. The multiparameter sensitivity analysis showed that different combinations of parameters produced similar model responses. Model sensitivity to grid resolution presented similar flow volumes for different digital elevation model resolutions, whereas event peak and duration (for both flow and sediment flux) were highly sensitive to changes in grid size.     

Exchange of Nitrogen and Phosphorus Between a Shallow Lagoon and Coastal Waters

West Falmouth Harbor, a shallow lagoon on Cape Cod, has experienced a threefold increase in nitrogen load since the mid- to late 1990s due to input from a groundwater plume contaminated by a municipal wastewater treatment plant. We measured the exchange of nitrogen and phosphorus between the harbor and the coastal waters of Buzzards Bay over several years when the harbor was experiencing this elevated nitrogen load. During summer months, the harbor not only retained the entire watershed nitrogen load but also had a net import of nitrogen from Buzzards Bay. During the spring and fall, the harbor had a net export of nitrogen to Buzzards Bay. We did not measure the export in winter, but assuming the winter net export was less than 112 % of the load, the harbor exported less than half of the watershed nitrogen load on an annual basis. For phosphorus, the harbor had a net import from coastal waters in the spring and summer months and a net export in the fall. Despite the large increase in nitrogen load to the harbor, the summertime import of phosphorus from Buzzards Bay was sufficient to maintain nitrogen limitation of primary productivity during the summer. Our findings illustrate that shallow systems dominated by benthic producers have the potential to retain large terrestrial nitrogen loads when there is sufficient supply of phosphorus from exchange with coastal waters.     

Mapping subvertical discontinuities in rock cuts using a 400-MHz ground penetrating radar antenna

Hidden subvertical discontinuities oriented parallel to subparallel to the exposed faces of outcropping sandstone were effectively mapped at three different study sites in central Missouri using a ground-penetrating radar system (GPR) equipped with a 400-MHz monostatic antenna and a survey wheel. At each site, a suite of 2-D ground-penetrating radar profiles were acquired along multiple closely spaced traverses on relatively smooth exposed rock surfaces. Time-zero correction was applied to the raw GPR data which were then processed using band-pass filtering, range and display gain, color transformation, and deconvolution techniques. Pseudo 3D images of each identified discontinuity at each site were constructed based on the interpretation of the nonmigrated ground-penetrating radar profiles. These pseudo 3D images were hand-migrated and transformed into true 3D images which depict variable depths at “perpendicular horizontal distance” to each discontinuity relative to the exposed rock face. The results demonstrate that GPR can be used to detect and map hidden discontinuities. This information can then be used for rock slope stability analysis and rock engineering purposes.