Two-region non-Darcian flow toward a well in a confined aquifer

We have derived an analytical solution for two-region flow toward a well in a confined aquifer based on a linearization method. The two-region flow includes Izbash non-Darcian flow near the well and Darcian flow in the rest of the aquifer. The wellbore storage is also considered. The type curves in the non-Darcian and Darcian flow domains are obtained by a numerical Laplace inversion method incorporated in MATLAB programs. We have compared our results with the one-region Darcian flow model (Theis). Our solutions agree with those of Sen [Sen Z. Type curves for two-region well flow. J Hydr Eng 1988;114(12):1461–84] which were obtained using the Boltzmann transform at late times for fully turbulent flow, while some difference has been found at early and moderate times. We have defined a dimensionless non-Darcian hydraulic conductivity term which is shown to be a key parameter for analyzing the two-region flow. A smaller dimensionless non-Darcian hydraulic conductivity results in a larger drawdown in the non-Darcian flow region at late times. However, the dimensionless non-Darcian hydraulic conductivity does not affect the slope of the dimensionless drawdown versus the logarithmic dimensionless time in the non-Darcian flow region at late times. The dimensionless non-Darcian hydraulic conductivity does not affect the late time drawdown in the Darcian flow region.     

Non-Darcian flow to a well in an aquifer–aquitard system

In this study, we use a linearization procedure and a finite difference method to solve non-Darcian flow to a well in an aquifer–aquitard system. The leakage effect is considered. Flow in the aquifer is assumed to be non-Darcian and horizontal, whereas flow in the aquitard is assumed to be Darcian and vertical. The Izbash equation [Izbash SV. O filtracii V Kropnozernstom Materiale. USSR: Leningrad; 1931 [in Russian]] is employed to describe the non-Darcian flow. The wellbore storage is also considered in this study. An approximate semi-analytical solution has been obtained by the linearization procedure, and a numerical solution has been obtained by using a finite difference method. The previous solutions for Darcian flow case and non-Darcian flow case without leakage can be described as special cases of the new solutions. The error caused by the linearization procedure has also been analyzed. The relative error caused by the linearization procedure is nearly 100% at early times, and decreases to zero at late times. We have also compared the results in this study with Wen et al. [Wen Z, Huang G, Zhan H. A numerical solution for non-Darcian flow to a well in a confined aquifer using the power law function. J Hydrol, 2008d [in revision]] in which the leakage effect is not considered, and Hantush and Jacob [Hantush MS, Jacob CE. Non-steady radial flow in an infinite leaky aquifer. Trans Am Geophys Union 1955;36(1):95–100] who investigated a similar problem in Darcian flow case. The comparison of this study and Wen et al. (2008d) indicates the dimensionless drawdown in the aquifer with leakage is less than that without leakage, and the leakage has little effect at early times. The comparison between the results of this study and that of Hantush and Jacob (1955) indicates that the dimensionless drawdown in the aquifer for non-Darcian flow is larger at early times and smaller at late times, than their counterparts for Darcian flow. A larger dimensionless non-Darcian conductivity k_D results in a smaller dimensionless drawdown in the aquifer at late times, and leads to a larger dimensionless drawdown in the aquifer at early times. A smaller dimensionless leakage parameter B_D results in a smaller drawdown at late times, and the leakage does not affect the early-time drawdown. The analysis of the dimensionless drawdown inside the well has also been included in this study when the wellbore storage is considered.     

On the integral equation of the continuous-state stochastic reservoir with seasonally varying autocorrelated inflows

The paper deals with continuous-state reservoirs in discrete time, in terms of the (continuous) Moran transform U_t = Z_t + X_t of the (discontinuous) storage levels Z_t. In the non-seasonal case, if the inflows X_t have a first-order Markov structure, the modified storage levels U_t have a second-order Markov structure. It is shown that, for each value of u, the limit p(u,v) as t → ∞ of the joint density p_{t + 1} (u, v) of U_{t+ 1} at u and U_t at v satisfies a certain non-standard first-order linear integral equation of the form

or, equivalently, a standard second-order equation of the form

In the seasonal case, with two seasons, the corresponding bivariate density p(u,v) for consecutive seasons satisfies a second-order integral equation of the form

the kernel being a function of the season-to-season inflow transition densities.     

Imposex in Morula Granulata as Bioindicator of Tributyltin (TBT) Contamination in the Dampier Archipelago, Western Australia

This study was undertaken to examine the use of imposex in the tropical neogastropod whelk Morula granulata as bioindicator of tributyltin (TBT) contamination. Samples were collected from sites throughout the Dampier Archipelago; both impact sites close to intensive shipping activity and remote control sites. TBT in this area originates mainly from commercial vessels using the Port of Dampier that are >25 m in length and exempt from legislation controlling the use of TBT in antifouling paints. A field survey was undertaken in July 1997; 100 M. granulata were taken at each of 18 sites with varying vessel activity, and examined for imposex. The level of TBT contamination in the Archipelago was low; ranging from <0.3 to 25 ng Sn l⁻¹ in water,<0.3 to 80 ng Sn g−¹ in oyster tissue (Saccostrea cuccullata – prey of M. granulata), and <0.3 to 33 ng Sn g⁻¹ in whelks. Percentage of imposex in M. granulata ranged from 0% to 57%, and was correlated with distance to vessel activity (r=0.646, p<0.05). The only significant relationship of imposex with butyltin concentrations was with TBT in oysters (r=0.826, p<0.05), and dibutyltin (DBT) in whelks (r=0.783, p<0.05). It is postulated that this is because of interspecies variation in ability to metabolize TBT to DBT and MBT (monobutyltin), as well as induction of imposex in earlier life stages, so that imposex in adults could be a reflection of prior contamination. M. granulata is a suitable bioindicator of TBT contamination in the Dampier Archipelago, which could be used throughout the Indo-Pacific given its abundance in the region.     

An analytical solution for non-Darcian flow in a confined aquifer using the power law function

We have developed a new method to analyze the power law based non-Darcian flow toward a well in a confined aquifer with and without wellbore storage. This method is based on a combination of the linearization approximation of the non-Darcian flow equation and the Laplace transform. Analytical solutions of steady-state and late time drawdowns are obtained. Semi-analytical solutions of the drawdowns at any distance and time are computed by using the Stehfest numerical inverse Laplace transform. The results of this study agree perfectly with previous Theis solution for an infinitesimal well and with the Papadopulos and Cooper’s solution for a finite-diameter well under the special case of Darcian flow. The Boltzmann transform, which is commonly employed for solving non-Darcian flow problems before, is problematic for studying radial non-Darcian flow. Comparison of drawdowns obtained by our proposed method and the Boltzmann transform method suggests that the Boltzmann transform method differs from the linearization method at early and moderate times, and it yields similar results as the linearization method at late times. If the power index n and the quasi hydraulic conductivity k get larger, drawdowns at late times will become less, regardless of the wellbore storage. When n is larger, flow approaches steady state earlier. The drawdown at steady state is approximately proportional to r¹⁻ⁿ, where r is the radial distance from the pumping well. The late time drawdown is a superposition of the steady-state solution and a negative time-dependent term that is proportional to t^{(1−n)/(3−n)}, where t is the time.     

Alpha-effect in flux ropes and sheets

The alpha-effect is studied at large magnetic Reynolds numbers (R) for certain simple steady flows. The flux sheets created by a two-dimensional cellular flow sustain an of order R^−1/2 at large values of R. For an analogous helical axisymmetric motion the mean electromotive force gives a finite nonzero limit of as R → ∞, this being a consequence of the existence of an intense axial flux rope. The latter estimate may be typical of steady three-dimensional motion.     

PCBs and chlorinated pesticides in the atmosphere and aquatic organisms of Ross Island, Antarctica

PCBs, p,p′-DDT, p,p′-DDE and lindane (γ-hexachlorocyclohexane) were monitored in the lower atmosphere of Ross Island, in Antarctica for 2 yr. Geometrical means were 15.2 pg m⁻³ for PCBs, 2.0 pg m⁻³ for p,p′-DDT, 1.0 pg m⁻³ for p,p′-DDE and 25.8 pg m⁻³ for lindane. Atmospheric levels of lindane were positively correlated with temperature, and a significant difference was found between spring-summer and summer-winter concentrations. No season related differences were found for the other chlorinated hydrocarbons, possibly owing to their lower vapour pressure and the cold climate. Periods with increased atmospheric levels of PCBs and DDT compounds were recorded. Lindane, p,p′-DDE and PCBs were present in fish and zooplankton sampled close to Ross Island. Pollutant levels in the zooplankton (on an extractable fat basis) were highest during the Antarctic spring and autumn and were inversely correlated to their fat content.     

Towards a more ecologically relevant assessment of the impact of heavy metals on the photosynthesis of the seagrass,Zostera capricorni

This in situ study used photosynthetic activity (measured as chlorophyll a fluorescence) and photosynthetic pigment concentrations to assess the effect of copper, cadmium, lead and zinc on the seagrass Zostera capricorni. Custom-made portable in situ exposure (PIE) chambers were developed so seagrasses could be dosed within the meadow. Z. capricorni was exposed to 0.1 and 1 mg l⁻¹ of metal solutions for 10 h. During this time and for the subsequent four-day recovery period, the effective quantum yield of photosystem II (PS II) (ΔF/Fm^′) was measured. While the results were variable, copper and zinc exposed samples had a depressed ΔF/Fm^′ during the exposure period. Samples exposed to zinc recovered to pre-exposure levels but those exposed to copper did not. Cadmium and lead did not impact on the chlorophyll a fluorescence and the chlorophyll pigment data supported these findings. This study presents an innovative new application of chlorophyll a fluorescence stress assessment.     

A spatial analysis of structural controls on Karst groundwater geochemistry at a regional scale

The coupled spatial investigation of the geometrical and geochemical properties of a chalk karstic aquifer provides information on the degree to which geologic structure controls aquifer functioning and groundwater quality. Major ion concentrations in the chalk aquifer of the Haute-Normandie region (France) were measured at a high spatial resolution (more than 100 sampling sites over a 6000 km² area) and mapped. The first observation is a continuity of the geochemical properties, in spite of the karstic properties of the aquifer principal components analysis of geochemical maps revealed two types of spatial distributions: ions with an autochthonous origin (Ca²⁺, HCO₃), and ions with a principally allochthonous origin (Cl⁻, Na⁺, , ). Mg²⁺ was categorised as both autochthonous (chalk dissolution) and allochthonous (brought in by infiltration of Tertiary deposits). To better understand the spatial distribution of the geochemistry, the aquifer geochemistry was compared to the physical properties of the aquifer, in particular aquifer thickness (representing aquifer geometry) and piezometric level (representing aquifer flow). Use of spatial correlation between the geochemical and the geometrical properties provided insight regarding the directional structure of the data and give evidence of directional relations between geochemical and geometrical properties. The degree of mineralisation (principally composed of Ca²⁺ and ions) increased along the direction of flow, corresponding to an increase in chalk dissolution rate along the flowpath. The steepest mineralisation gradients were related to an increase in the Mg/Ca ratio, evidence of longer residence times and corresponding to zones where aquifer flow capacity is limited because of a decrease of the thickness of the flow section (anticlines or faults). These results highlight the dominant role played by the geometry and the structural context in controlling aquifer geochemistry.     

The expansion and decline of charophyte communities in lakes within the Sejny Lake District (north-eastern Poland) and changes in water chemistry

Over the last 25 years, considerable changes have been observed in the littoral vegetation of 11 lakes within the Sejny Lake District. In eight lakes (Dowcień, Jurkowo, Kunis, Miałkie, Pilwie, Płaskie, Wiłkokuk, Zelwa) where Chara species dominated, their communities declined or disappeared completely. In four of these lakes (Dowcień, Jurkowo, Kunis and Miałkie), charophytes were replaced by communities of the class Potametea (mainly Nupharo-Nymphaeetum albae and Ceratophylletum demersi). In three other lakes (Pilwie, Płaskie and Wiłkokuk), the area covered by charophyte communities decreased. In Lake Zelwa, Charetum rudis and Charetum jubatae disappeared while Charetum tomentosae and Charetum fragilis expanded. A total area of charophyte communities remained approximately the same. A considerable increase in the water total hardness and concentration was noted in each of the eight lakes. In most of them increased levels of dissolved organic matter (measured as COD-KMnO₄) and pH were detected as well. However, the increase in the concentration was smaller in Lakes Wiłkokuk and Zelwa in which only slight changes in the vegetation were observed. In the above two lakes, the concentration of dissolved organic matter decreased and low water colour values are now noted. In three lakes (Długie, Dmitrowo and Gajlik), the charophyte communities tended to expand under the condition of moderate or high (Lake Długie) concentration, while dissolved organic matter concentrations were low. These lakes are also characterized by low values of colour. Charophyte communities may persist in lakes for a long period of time even when there is a moderate increase in concentration until the colour of water distinctly increases (>20 mg Pt L⁻¹).     

Analysis of the time-scaling behaviour in the sequence of the aftershocks of the Bovec (Slovenia) April 12, 1998 earthquake

Analysis of the temporal fluctuations, described by a power spectral density S(f) decaying as 1/f with exponent >0, in the sequence of aftershocks of the Bovec (Slovenia) earthquake of April 12, 1998 has been performed, using the Allan Factor (AF) method. The sequence of the occurrence times of the events with threshold magnitude M_th≥2.0 is characterised by a scale-invariant behaviour from the time scale τ2·10⁴ s with a scaling coefficient 0.9, evaluated by a least-square method. By gradually increasing the threshold magnitude up to M_th=2.9, the AF curves, associated, respectively with the processes of selected events with magnitude M≥M_th, indicate a monotonic decrease of the value of the scaling exponent , in the range of time scales stretching from 10² to 5·10⁵ s. For all the thresholds considered, the AF curves increase monotonically with the counting time and are well fitted by an increasing power-law function for counting times greater than approximately 2·10⁴ s. This monotonic power-law increase indicates the presence of fluctuations on many time scales and therefore of fractal clustering.     

Antibiotic Resistance and Plasmid Profiles of Vibrio Isolates from Cultured Silver Sea Bream, Sparus sarba

A total of 51 potential pathogenic vibrios were isolated from moribund silver sea bream (Sparus sarba) collected from fish farms in Hong Kong. Using the API 20E system and the scheme of Alsina and Blanch (1994), 7 species were identified from all isolates. These species were Vibrio alginolyticus (24 strains), Vibrio vulnificus (12 strains), Vibrio parahaemolyticus (7 strains), Vibrio logei (4 strains), Vibrio pelagius II (2 strains), Vibrio fluvialis (1 strain) and Vibrio meditterranei (1 strain). The three dominant species (V. alginolyticus, V. vulnificus and V. parahaemolyticus) were confirmed to be virulent to sea bream by experimental challenge.

All isolates were screened for plasmid DNA by agarose gel electrophoresis and tested for susceptibility to 16 antibiotics by the agar dilution method. Of the 51 isolates examined, all strains were sensitive to ceftriaxone, streptomycin, nalidixic acid and rifampicin, and almost all were sensitive to ceftazidime, netilimicin, chloramphenicol and sulfamethoxazole. Most isolates were resistant to ampicillin (60.8%), cefuroxime (66.7%), amikacin (55%), kanamycin (58.8%) and trimethoprim (76.5%). Fifteen of the 51 isolates harboured 1–4 plasmids, with sizes ranging from 9 to 123 kb. Both the plasmids and the associated antibiotic resistance (ampicillin, cefuroxime and trimethoprim) of 9 isolates could be transferred to a recipient by single-step conjugation. However, the frequencies were very low, ranging from 10⁻¹¹ to 10⁻⁹. The present results indicate that resistance to these antibiotics is chromosomal.     

A case study on fuel oil contamination in a mangrove swamp in Hong Kong

Mangroves commonly found along tropical and subtropical coastlines are susceptible to oil pollution. In December 2000, around 500 1 m tall Kandelia candel saplings at the age of 3–5 years old located at the foreshore region of Sheung Pak Nai swamp, Hong Kong SAR, were found to be damaged by oil pollution. More than 80% of the saplings were either dead or washed away and leaving less than 5% healthy saplings with dense green leaves. Elevated concentrations of light n-alkanes (ranging from n-C₁₄ to n-C₂₀), pristane and phytane were recorded in surface sediments collected in December 2000. The ratio between light and total n-alkanes was 0.4. The total petroleum hydrocarbons (60–80 μg g⁻¹ TPH) and unresolved complex mixtures (60–70 μg g⁻¹ UCM) were higher than the background values of other mangrove sediments in Hong Kong, which were 40 and 20 μg g⁻¹, respectively. In certain root zone sediments, TPH concentrations were above 1000 μg g⁻¹. These results suggest that surface sediments in Sheung Pak Nai were contaminated by petroleum oil, most likely by illegal discharge of fuel oil which occurred between 1998 and 2002. One year later, in December 2001, unhealthy saplings had recovered and re-grown. The concentrations of TPH and UCM in sediments declined to around 40 μg g⁻¹, pristane and phytane dropped by 80%, and the ratio of light to total n-alkanes was 0.15, suggesting that residual oil in sediments was weathered leading to a remarkable recovery of the unhealthy saplings.     

Organochlorine pesticide residues in sediments and biota from the coastal area of Dar es Salaam city,Tanzania

Sediment and biota samples were collected from Msimbazi and Kizinga rivers and from the coastal marine environment of Dar es Salaam during both dry and wet seasons. The samples were analyzed for various organochlorine pesticide residues using GC-ECD and GC-MS. Dieldrin, p,p^′-DDT, p,p^′-DDE, p,p^′-DDD, o,p^′-DDT and γ-HCH were detected at significantly greater concentrations above the method detection limits. Recoveries of pesticide residues ranged 86.5–120% in sediments and 62–102% in biota. The average concentrations of total DDT in sediments for the two seasons were almost the same. Biota samples showed significant difference in levels of residues depending on mode of feeding and age of analyzed biota. p,p^′-DDT to total DDT ratios in all matrices indicated recent use of DDT. The levels of residues in sediments suggest possible adverse effects to humans consuming biota that are directly exposed to the sediments. This effect, which is associated with bio-concentration of residues in the tissues of edible aquatic biota, should not affect wildlife. In spite of concern on the adverse effect to humans posed by pesticide residues in sediments, edible biota from the examined water bodies were found to be safe for human consumption.     

The Impact of the Herbicide Diuron on Photosynthesis in Three Species of Tropical Seagrass

The impact and recovery from exposure to the herbicide diuron [DCMU; 3-(3′,4′-dichlorophenyl)-1,1-dimethylurea] was assessed for three tropical seagrasses, maintained in outdoor aquaria over a 10-day period. Photosynthetic stress was detected using chlorophyll a fluorescence, measured with a Diving-PAM (pulse amplitude modulated fluorometer). Exposure to 10 and 100 μg l⁻¹ diuron resulted in a decline in effective quantum yield (ΔF/F_m^′) within 2 h of herbicide exposure in Cymodocea serrulata, Halophila ovalis and Zostera capricorni. Effective quantum yield also declined over the first 24 h of exposure in H. ovalis at even lower diuron concentrations (0.1 and 1.0 μg l⁻¹). Effective quantum yield in H. ovalis and Z. capricorni was significantly depressed at all diuron concentrations (0.1–100 μg l⁻¹) after 5 days exposure, whereas effective quantum yield in C. serrulata was only significantly lower in plants exposed to highest diuron concentrations (10 and 100 μg l⁻¹). Effective quantum yield depression was present 5 days after plants exposed to 10 and 100 μg l⁻¹ diuron were returned to fresh seawater. These results indicate that exposure to herbicide concentrations present in nearshore Queensland sediments present a potential risk to seagrass functioning.     

Steady state analysis of unsaturated flow above a shallow water-table aquifer drained by ditches

The one-dimensional pressure head profile above a fixed water-table was studied for different steady infiltration rates. As shown in previous studies, when the infiltration rate (q_in) is less than the soil saturated hydraulic conductivity (Ks), this profile can be divided into two parts: (1) from the water-table surface (z₀) to an elevation z_γ, the pressure head varies from 0 to a value h_γ such as K(h_γ)=−q_in; (2) above the elevation z_γ the pressure head is constant and equal to h_γ. Above the water-table the zone where the pressure head is variable has been called ‘transition zone’. Its height is shown to be highly variable and to depend on soil properties as well as on the infiltration rate. This transition zone is not the ‘capillary fringe’ as defined by Gillham (Gillham R.W., 1984. The capillary fringe and its effect on water-table response. J. Hydrol. 67, 307–324). Numerical experiments performed with HYDRUS-2D^® for the case of a water-table drained by parallel ditches have shown that the height of the transition zone is similar in the one-dimensional profile and in the two-dimensional system as long as the local slope of the water-table is small. This result is important since in a two-dimensional system, the transition zone is the place where all the horizontal component of the unsaturated flow occurs. The ratio of the horizontal component of the unsaturated flow vs the total horizontal flow in both the unsaturated and saturated zones has been computed. For a given soil, this ratio decreases as the infiltration rate increases; for a given infiltration rate, the soil with the thinnest transition zone transfers the largest amount of water above the water-table.     

川06井水位固体潮效应变化初探

通过分析承压井孔储水效应对井水位固体潮效应(井潮)的影响, 重点研究了含水层储水系数S和导水系数T对潮汐排水响应的影响, 以探讨井潮变化形成的初步原因. 以川06井为例, 通过对实测井孔水位数据进行潮汐分析, 得到井孔储水效应引起的相对振幅A和相位移η, 进而分析A和η随时间变化的规律以及两者之间的关系. 结果表明, 井潮变化主要受井孔储水效应变化影响; 在潮汐排水响应阶段, 井潮主要受含水层导水系数T(或渗透系数K)影响, 而受储水系数S影响较小.      

Clam burrowing behaviour: Inhibition by copper-enriched sediment

Burrowing behaviour is adaptive and allows clams to escape predation; yet the effects of potentially toxic metals on such behaviour have not been adequately investigated. In natural marine sediment contaminated with copper the time for littleneck clams (Protothaca staminea) to achieve complete burial was recorded. Above a threshold of 5.8 μg g⁻¹ Cu added to dry sediment, the time for 50% of the clams to burrow (ET₅₀) increased logarithmically with increasing sediment copper concentration according to:

logET₅₀ = 0.15 (CU) - 1.37 (n = 4, r² = 0.98)

where ET₅₀=time in hours for 50% of clams to burrow and Cu=μg g⁻¹ Cu in dry sediment. Previously exposed clams had both a lower threshold to Cu and a longer reburrowing time (ET₅₀). Clams exposed to sediment mixed with Chelex-100®-sorbed copper showed no significant change in burrowing time. Bioassays based on claim burrowing behaviour can measure both bioeffectiveness of sediment-sorbed metals and a sublethal effect with ecological meaning.     

Organochlorines in Hong Kong Fish

Muscle samples from 15 species of fish (n=1) purchased from markets in Hong Kong and 10 liver samples of tilapia (Tilapia mossambica) collected from the Shing Mun River were analysed for organochlorines (polychlorinated biphenyls, ΣDDTs, hexachlorobenzene, hexachlorocyclohexanes, chlordanes, mirex and dieldrin). Polychlorinated biphenyls (PCBs) were detected in 15 market fish samples but the levels were very low (around 1 ng g⁻¹). PCB levels in tilapia livers collected from Tai Wai (29.3–65.1 ng g⁻¹) were higher than those from Fo Tan (3.5–23.2 ng g⁻¹) suggesting that there may be some local point soucres. ΣDDTs were detected in all samples, ranging from 3.3 to 75.6 ng g⁻¹ in the market fish and from 7.1 to 88.8 ng g⁻¹ in tilapia. The DDE/DDT ratios in the market fish (0.12–0.75) showed higher variability than those of tilapia (0.30–0.46), suggesting that some of the market fish may have been collected from areas where DDT was recently used. Results of this preliminary study show that organochlorine levels in Hong Kong market fish are low and do not cause any concern for human consumption. An on-going monitoring program, however, is recommended.