Point source identification using a simple permutation test: a case study of elevated PCDD/F levels in ambient air and soil and their relation to the distance to a local municipal solid waste incinerator

A point source for environmental pollution may be identified by the spatial association of samples collected in the vicinity of the emission source. In this study, we used a simple permutation test to explore the spatial correlations between the ambient air and soil polychlorinated dibenzo-p-dioxin/dibenzofuran (PCDD/F) congener profile concentrations and their distances to a municipal solid waste incinerator (MSWI) in northern Taiwan. Kendall’s τ was applied for the correlation between the PCDD/F concentrations of a total of 9 airborne and 12 soil sampling sites and their distances to the incinerator. The results showed that the correlation was highly negatively associated for soil PCDD/Fs (τ = −0.515, p-value = 0.011) and was marginal for airborne PCDD/Fs (τ = −0.667, p-value = 0.055), conditional on the seasonal wind attributions. For the 17 individual congener profiles, most of the air and soil samples had a negative association, with 7 and 10 reaching statistical significance (p-value < 0.05) or borderline significance (p-value < 0.10), respectively. The permutation test provides a simple and robust statistical method for the identification of a potential pollution source from a limited sample size. Our analytical results show that dioxin emissions from the investigated MSWI had a definite environmental impact on the surrounding area.     

Point source identification using a simple permutation test: a case study of elevated PCDD/F levels in ambient air and soil and their relation to the distance to a local municipal solid waste incinerator

A point source for environmental pollution may be identified by the spatial association of samples collected in the vicinity of the emission source. In this study, we used a simple permutation test to explore the spatial correlations between the ambient air and soil polychlorinated dibenzo-p-dioxin/dibenzofuran (PCDD/F) congener profile concentrations and their distances to a municipal solid waste incinerator (MSWI) in northern Taiwan. Kendall’s τ was applied for the correlation between the PCDD/F concentrations of a total of 9 airborne and 12 soil sampling sites and their distances to the incinerator. The results showed that the correlation was highly negatively associated for soil PCDD/Fs (τ = −0.515, p-value = 0.011) and was marginal for airborne PCDD/Fs (τ = −0.667, p-value = 0.055), conditional on the seasonal wind attributions. For the 17 individual congener profiles, most of the air and soil samples had a negative association, with 7 and 10 reaching statistical significance (p-value < 0.05) or borderline significance (0.05 < p-value < 0.10), respectively. The permutation test provides a simple and robust statistical method for the identification of a potential pollution source from a limited sample size. Our analytical results show that dioxin emissions from the investigated MSWI had a definite environmental impact on the surrounding area.     

Polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans in sediments from Hong Kong

Concentrations of 2,3,7,8-substituted polychlorinated dibenzo-p-dioxins (PCDDs) and polychlorinated dibenzofurans (PCDFs) were determined in 14 sediment samples collected from four sites in the Mai Po Marshes Nature Reserve (within a RAMSAR Site) and from another six sites in Victoria Harbour and along the Hong Kong coastline. Elevated levels of PCDDs, and particularly OCDD, were detectable in all samples collected from the Mai Po Marshes and five of the six sites. In contrast to PCDDs, PCDFs were mainly found in sediment samples collected from industrial areas (Kwun Tong and To Kwa Wan) in Victoria Harbour. PCDD/F levels and congener profiles in the samples from the Mai Po Marshes Nature Reserve in particular show strong similarities to those reported in studies which have attributed similar elevated PCDD concentrations to nonanthropogenic PCDD sources.     

Early contamination of European flounder (Platichthys flesus) by PCDD/Fs and dioxin-like PCBs in European waters

Contamination levels and profiles of 7 polychlorinated-p-dioxins, 10 polychlorinated furans (PCDD/Fs) and 12 dioxin-like polychlorinated biphenyls (dl-PCBs) were investigated in juvenile European flounder (Platichthys flesus) captured in different nursery areas in the northeastern Atlantic coast across its geographical distribution range. The toxic equivalent concentrations (WHO-TEQ_fish) were also determined in order to evaluate which P. flesus population was more exposed to dioxin-like toxicity. Juveniles caught in the Sørfjord (Norway) showed the lowest WHO-TEQ_fish concentration (0.052 pg WHO-TEQ_fish g⁻¹ wet weight) whereas the highest value was observed in fish from the Wadden Sea (The Netherlands; 0.291 pg WHO-TEQ_fish g⁻¹ ww), mainly due to the greater contribution of 2,3,7,8-tetrachlorodibenzo-p-dioxin, the most toxic congener. Nonetheless, when comparing the results with existent tissue residue-based toxicity benchmarks, no adverse effects resulting from PCDD/Fs and dl-PCBs are expected to occur in flounder from the studied systems.     

Polychlorierte Dibenzo-p-dioxine, -furane (PCDD/F) und Analoga in Oberflächensedimenten und datierten aquatischen Sedimentkernen

Polychlorinated Dibenzo-p-dioxins, -furans (PCDD/F) and Analogues in Surface Sediments and Dated Aquatic Sediment Cores Aquatic sediments are the final sinks of unintentionally formed PCDD/Fs. Still not known are the conditions of formation of PCDD/Fs in the 8120 years old sediments in Japan and in 350 to 400 years old sediments in the Baltic. Sediments are means of source identification if transformation processes are taken into consideration in view of hydrological conditions. PCDD/Fs have been analysed 30 km from the source in sediments and over a distance of 50 km in benthos indicating a remarkable bioaccumulation potential. Elaborated data of the PCDD/F sediment contamination are missing as well as information on internal diagenetic processes, interactions at the sediment/water interface, and on the aerobic and anaerobic biodegradation and bioavailability of PCDD/Fs in sediments. Quality criteria for sediments and aquatic ecosystems remain to be elaborated.     

Dioxin-like compounds in sediments from the Daliao River Estuary of Bohai Sea: distribution and their influencing factors

The concentrations, compositional profiles, and potential ecological risk of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs), polychlorinated biphenyls (PCBs) and polychlorinated naphthalenes (PCNs) in sediments of the Daliao River Estuary were investigated. Total concentrations of PCDD/Fs, PCBs, and PCNs were in the range of 11.3-133.2 ng/kg dry weight (dw), 1 971-37 632 ng/kg dw and 33.1-284.4 ng/kg dw, respectively. The total TEQ values varied from 0.37 to 4.08 ng/kg dw, with the dominant contributions by PCDD/Fs, then by PCBs and PCNs. The spatial distributions of PCDD/Fs, PCBs and PCNs in the river estuary were much related to hydrodynamic conditions. The risk of contamination to the deeper sea was increased in the flood seasons. Moreover, our data confirmed that both organic matter in sediments and molecular properties of dioxin-like compounds were the factors which strongly influenced the partition behavior of these dioxin-like compounds between sediments and water phase in the estuarine zone.     

Temporal trends in dioxins (polychlorinated dibenzo-p-dioxin and dibenzofurans) and dioxin-like polychlorinated biphenyls in Baltic herring (Clupea harengus)

Polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) and dioxin-like polychlorinated biphenyl (dl-PCBs) concentrations in Baltic herring (Clupea harengus) have been relatively stable since the mid to late 1990s. It is unclear why concentrations in Baltic herring are not following the observed decreases in other environmental matrices. Here, changes in long-term temporal trends in Baltic herring were examined. A number of biological variables were examined alongside the temporal trends to investigate whether fish biology e.g., growth (length, weight, age), lipid content, reproductive phase or fishing date may provide an explanation for the temporal trends observed.     

Assessment of risk of PCDD/Fs and dioxin-like PCBs in marine and freshwater fish in Pearl River Delta,China

Fish consumption is known to be beneficial to human health. However since the age of industrialization, the released/disposed chemical pollutants into water systems make fish a source of various environmental toxicants to humans. In oceanic cities with heavy industrial activities, fish products contribute the greatest proportion of exposure to pollutants. In this study, risks and potential effects of dioxins to health of coastal populations in the Pearl River Delta were assessed. Concentrations of polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), and dioxin-like polychlorinated biphenyls (PCBs) were measured in common fish species purchased at local markets. Concentrations of total dioxins in fish ranged from 0.481 to 9.05 pg TEQ/g wet weight were similar to the lesser concentrations reported for fish from other countries. The greatest concentrations of dioxins were measured in mandarin fish, a carnivorous freshwater fish. Exposure of murine primary leydig and ovarian cells to 2,3,7,8-tetrachlorinated dibenzo-p-dioxins (2,3,7,8-TCDD) reduced the synthesis of progesterone, testosterone and/or estrogen. The reductions were probably via inhibitory effects on the expression of the steroidogenic enzymes, steroidogenic acute regulatory protein (StAR) and cytochrome P450 side-chain cleavage enzyme (P450scc). Based on these reproductive parameters, the concentrations of dioxins and dioxin-like residues represent a moderate health risk due to consumption of fish.     

Preliminary health risk assessment for polybrominated diphenyl ethers and polybrominated dibenzo-p-dioxins/furans in seafood from Guangzhou and Zhoushan, China

Dietary intake is one of the important routes of human exposure to brominated flame retardants (BFRs) such as polybrominated diphenyl ethers (PBDEs). The use of PBDEs may also result in exposure to polybrominated dibenzo-p-dioxins and dibenzofurans (PBDDs/DFs), as these compounds are impurities in technical mixtures of BFRs and can also be formed unintentionally by the same processes that generate chlorinated dioxins. This study determined the concentrations of polybrominated compounds in common seafood in Guangzhou and Zhoushan, and assessed the health risks of these chemicals via consumption of contaminated seafood. Seafood samples (fish, bivalves, shrimp, crab, and cephalopods) purchased from local markets in 2003 and 2004 were analyzed for PBDEs and PBDDs/DFs. The highest concentration of total PBDEs (46.3 ng g⁻¹ lipid wt.) was detected in fish from Guangzhou, in which BDEs 47 and 209 were the two predominant congeners. The total daily intakes of PBDEs, PBDDs, and PBDFs were, 946, 6.39, and 6.54 pg kg⁻¹ body weight (bw) in Guangzhou, and 489, 4.99, and 7.65 pg kg⁻¹ bw in Zhoushan, respectively. The hazard ratios for PBDDs and PBDFs were both greater than unity, indicating that these compounds may pose some health risks to the local population.     

Seasonal rainfall–runoff relationships in a lowland forested watershed in the southeastern USA

Hydrological processes of lowland watersheds of the southern USA are not well understood compared to a hilly landscape due to their unique topography, soil compositions, and climate. This study describes the seasonal relationships between rainfall patterns and runoff (sum of storm flow and base flow) using 13 years (1964–1976) of rainfall and stream flow data for a low‐gradient, third‐order forested watershed. It was hypothesized that runoff–rainfall ratios (R/P) are smaller during the dry periods (summer and fall) and greater during the wet periods (winter and spring). We found a large seasonal variability in event R/P potentially due to differences in forest evapotranspiration that affected seasonal soil moisture conditions. Linear regression analysis results revealed a significant relationship between rainfall and runoff for wet (r² = 0·68; p < 0·01) and dry (r² = 0·19; p = 0·02) periods. Rainfall‐runoff relationships based on a 5‐day antecedent precipitation index (API) showed significant (r² = 0·39; p < 0·01) correspondence for wet but not (r² = 0·02; p = 0·56) for dry conditions. The same was true for rainfall‐runoff relationships based on 30‐day API (r² = 0·39; p < 0·01 for wet and r² = 0·00; p = 0·79 for dry). Stepwise regression analyses suggested that runoff was controlled mainly by rainfall amount and initial soil moisture conditions as represented by the initial flow rate of a storm event. Mean event R/P were higher for the wet period (R/P = 0·33), and the wet antecedent soil moisture condition based on 5‐day (R/P = 0·25) and 30‐day (R/P = 0·26) prior API than those for the dry period conditions. This study suggests that soil water status, i.e. antecedent soil moisture and groundwater table level, is important besides the rainfall to seasonal runoff generation in the coastal plain region with shallow soil argillic horizons. Copyright © 2011 John Wiley & Sons, Ltd.     

Turbidity as a control on phytoplankton biomass and productivity in estuaries

In many coastal plain estuaries light attenuation by suspended sediments confines the photic zone to a small fraction of the water column, such that light limitation is a major control on phytoplankon production and turnover rate. For a variety of estuarine systems (e.g. San Francisco Bay, Puget Sound, Delaware Bay, Hudson River plume), photic-zone productivity can be estimated as a function of phytoplankton biomass times mean irradiance of the photic zone. Net water column productivity also varies with light availability, and in San Francisco Bay net productivity is zero (estimated respiratory loss of phytoplankton balances photosynthesis) when the ratio of photic depth (Z_p) to mixed depth (Z_m) is less than about 0.2. Thus whenever Z_p: Z_m < 0.2, the water column is a sink for phytoplankton production.Much of the spatial and temporal variability of phytoplankton biomass or productivity in estuaries is explained by variations in the ratio of photic depth to mixed depth. For example, phytoplankton blooms often coincide with stratification events that reduce the depth of the surface mixed layer (increase Z_p: Z_m). Shallow estuarine embayments (high Z_p: Z_m) are often characterized by high phytoplankton biomass relative to adjacent channels (low Z_p: Z_m). Many estuaries have longitudinal gradients in productivity that mirror the distribution of suspended sediments: productivity is low near the riverine source of sediments (low Z_p: Z_m) and increases toward the estuary mouth where turbidity decreases. Some of these generalizations are qualitative in nature, and detailed understanding of the interaction between turbidity and estuarine phytoplankton dynamics requires improved understanding of vertical mixing rates and phytoplankton respiration.     

Seasonal variations in pCO2 and its controlling factors in surface seawater of the northern East China Sea

Surface partial pressure of CO₂ (pCO₂), temperature, salinity, nutrients, and chlorophyll a were measured in the East China Sea (ECS; 31°30′–34°00′N to 124°00′–127°30′E) in August 2003 (summer), May 2004 (spring), October 2004 (early fall), and November 2005 (fall). The warm and saline Tsushima Warm Current was observed in the eastern part of the survey area during four cruises, and relatively low salinity waters due to outflow from the Changjiang (Yangtze River) were observed over the western part of the survey area. Surface pCO₂ ranged from 236 to 445 μatm in spring and summer, and from 326 to 517 μatm in fall. Large pCO₂ (values >400 μatm) occurred in the western part of the study area in spring and fall, and in the eastern part in summer. A positive linear correlation existed between surface pCO₂ and temperature in the eastern part of the study area, where the Tsushima Warm Current dominates; this correlation suggests that temperature is the major factor controlling surface pCO₂ distribution in that area. In the western part of the study area, however, the main controlling factor is different and seasonally complex. There is large transport in this region of Changjiang Diluted Water in summer, causing low salinity and low pCO₂ values. The relationship between surface pCO₂ and water stability suggests that the amount of mixing and/or upwelling of CO₂-rich water might be the important process controlling surface pCO₂ levels during spring and fall in this shallow region. Sea–air CO₂ flux, based on the application of a Wanninkhof [1992. Relationship between wind speed and gas exchange over the ocean. Journal of Geophysical Research 97, 7373–7382] formula for gas transfer velocity and a set of monthly averaged satellite wind data, were −5.04±1.59, −2.52±1.81, 1.71±2.87, and 0.39±0.18 mmol m⁻² d⁻¹ in spring, summer, early fall, and fall, respectively, in the northern ECS. The ocean in this study area is therefore a carbon sink in spring and summer, but a weak source or in equilibrium with the atmosphere in fall. If the winter flux value is assumed to have been the mean of autumnal and vernal values, then the northern ECS absorbs about 0.013 Pg C annually. That result suggests that the northern ECS is a net sink for atmospheric CO₂, a result consistent with previous studies.     

Evaluation of temporal consistency of snow depth drivers of a Rocky Mountain watershed in southern Alberta

Collecting spatially representative data over large areas is a challenge within snow monitoring frameworks. Identifying consistent trends in snow accumulation properties enables increased sampling efficiency by minimizing field collection time and/or remote sensing costs. Seasonal snowpack depth estimations during mid-winter and melt onset conditions were derived from airborne Lidar over the West Castle Watershed in the southern Canadian Rockies on three dates. Each dataset was divided into five sets of snow depth driver classes: elevation, aspect, topographic position index, canopy cover and slope. Datasets were quality controlled by eliminating snow depth values above the 99th percentile value, which had a negligible effect on average snow depths. Consistent trends were observed among driver classes with peak snow accumulation occurring within the treeline ecotone, north-facing aspects, open canopies, topographic depressions and areas with low slope angle. Although mid-winter class trends for each driver were similar and watershed-scale snow depth distributions were significantly correlated (0.76, p < .01), depth distributions within the same driver class of the three datasets were not correlated due to recent snowfall events, redistribution and settling processes. Trends in driver classes during late season melt onset were similar to mid-winter conditions but watershed scale distribution correlation results varied with seasonality (0.68 mid-winter 2014 and melt onset 2016; 0.65 mid-winter 2017 and melt onset 2016, p < .1). This is due to the differing stages of accumulation or ablation and the upward migration in the 0°C isotherm during spring, when snow depth can be declining in valley bottoms while still increasing at higher elevations. The observed consistency in depth driver controls can be used to guide future integrated snow monitoring frameworks.     

A high-resolution analyser for the measurement of ammonium in oligotrophic seawater

In this work, we describe a high-resolution fluorometric shipboard analyser and an improved method to determine NH₄⁺ in oligotrophic seawater. The limit of detection is <5 nM, calculated with 95% confidence level using the weighted regression line applied to the standard addition method using real samples prepared with low nutrient seawater from the Atlantic. The results are summarised and cross-compared with spiked artificial seawater (ASW) and spiked Milli-Q water samples. The analyser has a precision of ±1–4% with a high performance over a wide range from 5 nM to 25 μM. The methodology of NH₄⁺ analysis is based on the fluorescent product formed between o-pthaldialdehyde and NH₄⁺ in the presence of sulfite. Due to the high resolution of the developed system, we were able to study in depth the sensitivity of the method to salinity, amines, amino acids and potential interferences from particles/algae. The method was found to be sensitive to salinity variations, reducing the signal by up to 85% at 5 nM; this effect decreased at higher concentrations of ammonium. It was noted that the interference from amines at low concentrations was negligible; however, at either high amino acid or high amine concentrations, the signal was depressed. To test for the effect of particles on the system, the system was tested with samples containing phytoplankton (Dunaliella primolecta) cells at different concentrations prepared with ASW to simulate the effect of a phytoplankton bloom. This experiment assessed the potential impact of both particles and other potential fluorescence interferences from cells and/or ammonium leaching from cells. This experiment showed that a phytoplankton bloom could potentially have an impact of up to 12% on the signal of interest. Thus, we propose that this method is suitable for oligotrophic environments rather than coastal and eutrophic environments. The reagent was found to be stable for 17 days and standards of 1 μM were stable for 6 days under laboratory conditions. The developed analyser was successfully demonstrated in the North Atlantic Ocean, in an area of oligotrophic, low NH₄⁺ oceanic waters.     

Scaling model experiments on propagation of electromagnetic waves using salt solution to simulate the Earth's crust: Loma Prieta and Kobe earthquakes

Propagation of electromagnetic (EM) waves from an earthquake focus in the conductive Earth has been investigated using 1/1,000,000 scaling models taking earth-ionosphere and ocean-Moho plane parallel-plate waveguides into account. Microwaves at a frequency, ω_m, a million times higher than that of seismic EM signal (SEMS), ω, were generated at the model focus. They are propagated in a salt solution modeling the earth's crust and reflected by ocean, fault planes, ionosphere and Moho plane all made by aluminum. Distribution of EM power was mapped by scanning a detector antenna over the model Earth's surface. The skin depth, δ, calculated by the exact skin depth equation, 1/δ=ω(μ/2)^1/2 [(1+(1/ωρ)²)^1/2 −1]^1/2 where dielectric constant, and permeability, μ are the same but resistivity, ρ, 10⁻⁶ times smaller than that of Earth, gave 10⁻⁶ times small skin depth validating the model scaling index. Images for evanescent and wave-ripple standing waves disturbed by normal, strike-slip and dip-slip conductive fault planes have been obtained using an aluminum plate. The co-circular contour map above the epicenter due to evanescence was pushed to the north east direction from the epicenter by the presence of ocean for the Loma Prieta earthquake, while to north direction for the Kobe earthquake. The intensity of EM ULF emissions for the Loma Prieta earthquake is discussed quantitatively.     

Comparison of methods for quantifying active layer dynamics and bedload discharge in armoured gravel‐bed rivers

Several methods were employed in the Ardennian rivers (Belgium) to determine the depth of the active layer mobilized during floods and to evaluate the bedload discharge associated with these events. The use of scour chains has shown that the depth of the active layer is systematically less than the b‐axis of the average particle size (D₅₀) of the elements which compose the surface layer of the riffles. This indicates that only a partial transport exists during low magnitude floods. The bedload discharge has been evaluated by combining data obtained using the scour chains technique and the distance covered by tracers. Quantities of sediment transported during frequent floods are relatively low (0·02 t km^–2) due to the armour layer which protects the subsurface material. These low values are also related to the fact that the distance calculated for mobilized bedload only applies to tracers fitted with PIT (passive integrated transponder)‐tags (diameter > 20 mm), whereas part of the bedload discharge is composed of sand and fine gravel transported over greater distances than the pebbles. The break‐up of the armour layer was observed only once, for a decennial discharge. During this event, the bedload discharge increased considerably (2 t km^–2). The use of sediment traps, data from dredging and a Helley–Smith sampler confirm the low bedload transport in Ardennian rivers in comparison to the bedload transport in other geomorphological contexts. This difference is explained by the presence of an armoured layer but also by the imbricated structures of flat bed elements which increase the resistance to the flow. Finally, the use of the old iron industry wastes allowed to quantify the thickness of the bed reworked over the past centuries. In the Lembrée River, the river‐bed contains slag elements up to a depth of about 50 cm, indicating that exceptional floods may rework the bed to a considerable depth. Copyright © 2012 John Wiley & Sons, Ltd.     

Isotopic and geochemical characterization of salinization in the shallow aquifers of a reclaimed subsiding zone: The southern Venice Lagoon coastland

The coastal plain bordering the southern Venice Lagoon is a reclaimed lowland characterized by high subsidence rate, and ground level and water-table depth below sea level. In this agricultural region, where the surface hydrologic network is entirely artificially controlled by irrigation/drainage canals, salinization problems have long been encountered in soils and groundwaters. Here we use isotopic and geochemical tracers to improve our understanding of the origin of salinization and mineralization of the semi-confined aquifer (0–40 m), and the freshwater inputs to this hydrological system. Water samples have been collected at different seasons in the coastal Adriatic Sea, lagoon, rivers and irrigation canals, as well as in the semi-confined aquifer at depths between 12 and 35 m (14 boreholes), and in the first confined aquifer (three boreholes drilled between 40 and 80 m depth). Stable isotopes (δ¹⁸O and δD) and conductivity profiles show that direct saline intrusion from the sea or the lagoon is observed only in a restricted coastal strip, while brackish groundwaters are found over the entire topographic and piezometric depression in the centre of the study area. Fresh groundwaters are found only in the most western zone. The sharp isotopic contrast between the western and central regions suggests disconnected hydrological circulations between these two parts of the shallow aquifer. The border between these two regions also corresponds to the limits of the most strongly subsiding zone.Our results can be interpreted in terms of a four end-member mixing scheme, involving (1) marine water from the lagoon or the open sea, (2) alpine and pre-alpine regional recharge waters carried either by the main rivers Adige, Bacchiglione and Brenta (irrigation waters) or by the regional groundwater circulation, (3) local precipitation, and (4) evaporated waters infiltrated from the surface. Infiltration from the surface is also revealed by the stratification of the electrical conductivity profiles, showing that the brackish groundwaters are overlain by a shallow layer of less saline water all over the central depression. In the first confined aquifer, the groundwaters have isotopic compositions similar to the deep groundwaters of the Venetian confined aquifers (40–400 m depth). The isotopic data and the Br/Cl ratio show that the origin of the salinization of the phreatic aquifer can be ascribed to seawater intrusion alone, with no indication of the involvement of deep brines (identified at 450 m depth) in the process.The chemical composition of the saline and brackish groundwaters is characterized by an excess of sodium and a deficit of calcium compared to conservative mixing between fresh groundwaters and seawater. This suggests that the phreatic aquifer is progressively freshening, as a consequence of the beneficial influence of the extensive irrigation/drainage network, including raised canals acting as a hydraulic barrier along the coast. This freshening tendency may have been lasting since the reclamation in the mid-twentieth century, and has probably been accelerated by the ban on groundwater abstraction since the 1970s.     

Spin-up of a two-layer fluid in a rotating cylinder

Abstract

we report the results of experiments on the spin-up of two layers of immiscible fluid with a free upper surface in a rotating cylinder over a wide range of internal Froude numbers. Observations of the evolution of the velocity field by particle tracking indicates that spin-up of the azimuthal velocity in the upper layer take much longer than in a homogeneous fluid. Initially, spin-up occurs at a rate comparable to that of homogeneous fluid but, at high internal Froude number, a second phase follows in which the remaining lative motion decays much more slowly. Quantitative comparison of these measurements to the theory of Pedlosky (1967) shows good agreement.

Visualization of the interface displacement during spin-up detected the presence of transient azimuthal variations in the interface elevation over a wide range of Froude (F), Ekman (E), and Rossby (ε) number. nalysis of the occurrence of the asymmetric variations using the parameter space (Q, F), where Q = E ^1/2/ε, suggested by the baroclinic instability theory and experiments of Hart (1972), showed that the flow was stable for Q > 0.06 with no discernable dependence on F. This result, together with the prediction of Pedlosky's theory that radial gradient of potential vorticity in the two layers have opposite signs, suggests at the baroclinic instability mechanism was responsible for the asymmetries. The location and timing of these instabilities may account for the discrepancies between the observations and the Pedlosky (1967) theory.     

Downstream composition changes of acidic volcanic waters discharged into the Banyupahit stream, Ijen caldera, Indonesia

The crater lake of Kawah Ijen volcano contains extremely low pH (<0.4) waters with high SO₄ (70000 mg/kg), Cl (21000 mg/kg), F (1500 mg/kg), Al (5000 mg/kg), Fe (2000 mg/kg) and trace metal (Cu 0.5, Zn 4, Pb 3 mg/kg) contents. These brines seep outward through the western crater rim and reappear on the other side as streamlets, which form the headwaters of the Banyupahit stream. The Banyupahit first mixes with fresh rivers and thermal springs in the Ijen caldera and then irrigates a coastal agricultural plain which is 30 km from the summit crater.We discuss the downstream composition changes affecting the Banyupahit waters by using stable isotope, chemical and mineralogical data collected from sites along the stream length. The saturation of the stream waters with respect to minerals was evaluated with SOLVEQ and WATEQ4F and compared with the geochemical observations. An aluminous mineralogy (alunogen, pickeringite, tamarugite and kalinite) develops in the upper part of the Banyupahit due to concentration of the headwaters by evaporation. Downstream attenuation of dissolved element concentrations results principally from dilution and from mineral precipitation. The stream pH changes from 0 at the source to >4 close to the mouth. The δD and δ¹⁸O values and the relative SO₄–Cl–F contents of the Banyupahit waters indicate that the tributaries are mostly meteoric. Dissolved SO₄ in the acidic stream come only from the crater lake seepages and are not involved later in microbially mediated reactions, as shown by their δ³⁴S and δ¹⁸O values. Re-equilibration of the stream SO₄ oxygen-isotope composition with H₂O from tributaries does not occur.Calcium, SiO₂, Al, Fe, K and SO₄ behave non-conservatively in the stream waters. Gypsum, silica (amorphous or poorly ordered), a basic aluminum hydroxysulfate (basaluminite?), K-jarosite and amorphous ferric hydroxide may exert a solubility control on these elements along the entire stream length, or in certain stream sections, consistent with the thermochemical model results. Downstream concentration trends and mineral saturation levels suggest that precipitation of Sr-, Pb-rich barite and celestite consume Ba, Sr and Pb, whereas dissolved Cu, Pb and Zn may adsorb onto solid particles, especially after the junctions of the acidic stream with non-acidic rivers. We calculated that significant fluxes of SO₄, F, Cl, Al, SiO₂, Ti, Mn and Cu may reach the irrigation system, possibly causing serious environmental impacts such as soil acidification and induration.     

Numerical Simulations of Seasonally Oscillated Groundwater Dynamics in Coastal Confined Aquifers

Studies investigating the effects of inland recharge on coastal groundwater dynamics were carried out typically in unconfined aquifers, with few in confined aquifers. This study focused on the groundwater dynamics in confined aquifers with seasonally sinusoidally fluctuated inland groundwater head and constant sea level by numerical simulations. It is known that the mixing zone (MZ) of saltwater wedge in response to the seasonal oscillations of inland groundwater head swings around the steady-state MZ. However, our simulation results indicate that even the most landward freshwater-saltwater interface over a year is seaward from the steady-state location when the hydraulic conductivity K is ≤10⁻⁴ m/s under certain boundary conditions with given parameter values. That is, seasonal oscillations of inland groundwater head may reduce seawater intrusion in confined coastal aquifers when K ≤ 10⁻⁴ m/s. Sensitivity analysis indicates that for aquifers of K ≤ 10⁻⁴ m/s, the larger the inland head fluctuation amplitude is, the less the seawater intrudes. This is probably due to the reason that the seawater intrusion time decreases with the increase of fluctuation amplitude when K ≤ 10⁻⁴ m/s. Numerical simulations demonstrate that seasonal inland groundwater head oscillations promote the annual averaged recirculated seawater discharge across the seaward boundary.