Karst spring responses examined by process-based modeling

Ground water in karst terrains is highly vulnerable to contamination due to the rapid transport of contaminants through the highly conductive conduit system. For contamination risk assessment purposes, information about hydraulic and geometric characteristics of the conduits and their hydraulic interaction with the fissured porous rock is an important prerequisite. The relationship between aquifer characteristics and short-term responses to recharge events of both spring discharge and physicochemical parameters of the discharged water was examined using a process-based flow and transport model. In the respective software, a pipe-network model, representing fast conduit flow, is coupled to MODFLOW, which simulates flow in the fissured porous rock. This hybrid flow model was extended to include modules simulating heat and reactive solute transport in conduits. The application of this modeling tool demonstrates that variations of physicochemical parameters, such as solute concentration and water temperature, depend to a large extent on the intensity and duration of recharge events and provide information about the structure and geometry of the conduit system as well as about the interaction between conduits and fissured porous rock. Moreover, the responses of solute concentration and temperature of spring discharge appear to reflect different processes, thus complementing each other in the aquifer characterization.     

Transition probability/Markov chain analyses of DNAPL source zones and plumes

At sites where a dense nonaqueous phase liquid (DNAPL) was spilled or released into the subsurface, estimates of the mass of DNAPL contained in the subsurface from core or monitoring well data, either in the nonaqueous or aqueous phase, can be highly uncertain because of the erratic distribution of the DNAPL due to geologic heterogeneity. In this paper, a multiphase compositional model is applied to simulate, in detail, the DNAPL saturations and aqueous-phase plume migration in a highly characterized, heterogeneous glaciofluvial aquifer, the permeability and porosity data of which were collected by researchers at the University of Tübingen, Germany. The DNAPL saturation distribution and the aqueous-phase contaminant mole fractions are then reconstructed by sampling the data from the forward simulation results using two alternate approaches, each with different degrees of sampling conditioning. To reconstruct the DNAPL source zone architecture, the aqueous-phase plume configuration, and the contaminant mass in each phase, one method employs the novel transition probability/Markov chain approach (TP/MC), while the other involves a traditional variogram analysis of the sampled data followed by ordinary kriging. The TP/MC method is typically used for facies and/or hydraulic conductivity reconstruction, but here we explore the applicability of the TP/MC method for the reconstruction of DNAPL source zones and aqueous-phase plumes. The reconstructed geometry of the DNAPL source zone, the dissolved contaminant plume, and the estimated mass in each phase are compared using the two different geostatistical modeling approaches and for various degrees of data sampling from the results of the forward simulation. It is demonstrated that the TP/MC modeling technique is robust and accurate and is a preferable alternative compared to ordinary kriging for the reconstruction of DNAPL saturation patterns and dissolved-phase contaminant plumes.     

Magnetic cleansing of weathered/tarry oiled feathers--the role of pre-conditioners

Iron powder has previously been demonstrated to be effective in the removal, via magnetic harvesting, of a wide variety of oil contaminants from feathers and plumage. This study investigates the efficacy of magnetic cleansing for the removal from feathers of tarry contamination that has been allowed to weather. Clusters of feathers from Mallard duck (Anas platyrhnchos) and Little Penguin (Eudyptula minor) were completely immersed in a tarry contaminant and allowed to weather from one to fourteen days. The contaminant was removed using a magnetic cleansing protocol and the removal efficacy assessed gravimetrically. For one, seven and fourteen days of weathering, a final removal (after fourteen treatments) of more than 99% and 97% was achieved for duck feathers and penguin feathers, respectively. Repeating the experiments (for a seven-day weathering period) for both duck and penguin feathers, with the judicious application of a pre-conditioner (olive oil), further improved removal efficacy. A convenient method to screen for improved pre-conditioning agents is suggested.     

Boron isotopes as an artificial tracer

A field study was conducted using a combination of intrinsic and artificial tracers to estimate travel times and dilution during transport of infiltrate from a reclaimed water infiltration basin to nearby monitoring wells. A major study objective was to validate boric acid enriched in (10)B as an artificial tracer. Basin 10E at the Rio Hondo Spreading Grounds in Whittier, California, was the site of the test. The basin normally receives a mixture of treated municipal waste water, purchased State Project water, and local runoff from the San Gabriel River. Approximately 3.5 kg of (10)B-enriched boric acid was dispersed among 2.05 x 10(5) m(3) of basin water to initiate the experiment. The resultant median delta(11)B in the infiltration basin was -71 per thousand. Prior to tracer addition, the basin water had an intrinsic delta(11)B of +2 per thousand. Local monitoring wells that were used to assess travel times had delta(11)B values of +5 per thousand and +8 per thousand at the time of tracer addition. Analytic results supported an assumption that boron is conserved during ground water transport and that boron enriched in (10)B is a useful artificial tracer. Several intrinsic tracers were used to reinforce the boric acid tracer findings. These included stable isotopes of oxygen (delta(18)O) and hydrogen (deltaD), sulfate concentration, and the boron to chloride ratio. Xenon isotopes, (136)Xe and (124)Xe, also supported boron isotope results. Xenon isotopes were added to the recharge basin as dissolved gases by investigators from the Lawrence Livermore National Laboratory.     

Influence of small-scale heterogeneities on contaminant transport in fractured crystalline rock

We present a sequence of purely advective transport models that demonstrate the influence of small-scale geometric inhomogeneities on contaminant transport in fractured crystalline rock. Special weight is placed on the role of statistically generated variable fracture apertures. The fracture network geometry and the aperture distribution are based on information from an in situ radionuclide retardation experiment performed at Grimsel test site (Swiss Alps). The obtained breakthrough curves are fitted with the advection dispersion equation and continuous-time random walks (CTRW). CTRW is found to provide superior fits to the late-arrival tailing and is also found to show a good correlation with the velocity distributions obtained from the hydraulic models. The impact of small-scale heterogeneities, both in fracture geometry and aperture, on transport is shown to be considerable.     

A sediment quality triad assessment of the impact of copper mine tailings disposal on the littoral sedimentary environment in the Atacama region of northern Chile

A sediment quality triad (SQT) assessment was made of the impact of copper mine tailings disposal on littoral meiofaunal assemblages in the Atacama region of northern Chile. This situation is unusual in that the disposal is direct into the high-energy coastal system and not via a river estuary or other low-energy environment. This situation also allows for the examination of the impact of copper mine tailings in the absence of confounding effects from other pollutants. The three components of the SQT were: 1. an analysis of the bioavailable metals in both the sedimentary porewater and the adjacent seawater, 2. a microcosm bioassay of both sediments and seawaters using meiofaunal assemblages, and 3. quantitative field samples of the meiofaunal assemblages. Twelve study sites with varying degrees of impact were used, including three reference sites. The study identified that both the meiofaunal assemblage densities and taxa diversities decrease with increasing levels of bioavailable copper, that the Foraminifera and Harpacticoida are sensitive to copper, and that otoplanid Turbellaria are often characteristic of impacted sites; tailings also have both chemical and physical impacts on the environment. In some cases the physical impact of tailings is more important in excluding some organisms e.g. the interstitial polychaete, Saccocirrus sonomacus, from a site than is their chemical impact.     

A 3D hydrodynamic fate and transport model for herbicides in Sacca di Goro coastal lagoon (Northern Adriatic)

Sacca di Goro is a shallow coastal microtidal lagoon with a surface area of 26 km2, and an average depth of about 1.5m. Fresh water pollutant loads from Po River branches and several drainage canals lead to anthropogenic eutrophication, frequent summer anoxia crises and chemical contamination. Such events not only affect the lagoon ecosystem but also cause serious economic losses, the lagoon being the second largest producer of clams in Italy. The present work aims at using a fate model coupled with COHERENS 3D hydrodynamic model to simulate and to explain the spatial distribution and temporal variations of s-triazines herbicides in the Sacca di Goro lagoon. The simulation results of spatial and temporal dynamic behaviour of atrazine, simazine and terbuthylazine have been compared with experimental data obtained during an annual monitoring programme.     

Survival under chronic stress from sediment load: spatial patterns of hard coral communities in the southern islands of Singapore

Six reef sites were chosen along the west coasts of Singapore's southern islands, to: (1) quantitatively assess and compare coral community composition and structure, and recruitment rates, (2) assess the relationship between the aforementioned patterns and the environmental conditions, and (3) provide insights on potential processes that incorporate history at these study sites. Chronic exposure to high sediment load was the most obvious form of anthropogenic stress. Recruitment rates on ceramic tiles were low (1.4+/-1.0-20+/-14.7 recruits m(-2)year(-1)) but decreased towards the main island of Singapore as did hard coral cover and coral density. Coral fauna consisted of genera generally found in deeper waters (e.g., fungiids, foliose Oxypora, Leptoseris, and Echinopora) or those well adapted to turbid waters (e.g., Porites, Pectinia, Leptastrea, Montipora). Light extinction coefficient (K) and % live coral cover (%LCC) showed a strong and inverse curvilinear relationship (%LCC=13.60 *K(-3.40)). Similarly, the rate of sediment deposition (DFSPM) (Recruitment rate, RR=1.51-0.17 *DFSPM) and water clarity (RR=3.56-2.92 *K) exhibited strong and inverse relationships with recruitment rates. Although measured levels of the down-ward flux of suspended particulate matter and suspended solids were well within "normal" levels recorded in the literature, it was the proportion of benthic space, generic coral composition, and site history that offered compelling evidence of chronic exposure to increased sediment load. Clearly a reduction in both water clarity and live-coral cover has taken place since monitoring efforts began in the early 1970s, in fact coral cover has more than halved at all sites examined since the 1980s and benthic space was predominantly occupied by dead corals covered with sediment and filamentous algae.