Using electrical resistivity tomography to assess the effectiveness of managed aquifer recharge in a salinized coastal aquifer

Over 40 years, the detrital aquifer of the Plana de Castellón (Spanish Mediterranean coast) has been subjected to seawater intrusion because of long dry periods combined with intensive groundwater exploitation. Against this backdrop, a managed artificial recharge (MAR) scheme was implemented to improve the groundwater quality. The large difference between the electrical conductivity (EC) of the ambient groundwater (brackish water due to marine intrusion) and the recharge water (freshwater) meant that there was a strong contrast between the resistivities of the brackish water saturated zone and the freshwater saturated zone. Electrical resistivity tomography (ERT) can be used for surveying similar settings to evaluate the effectiveness of artificial recharge schemes. By integrating geophysical data with lithological information, EC logs from boreholes, and hydrochemical data, we can interpret electrical resistivity (ER) with groundwater EC values and so identify freshwater saturated zones. Using this approach, ERT images provided a high-resolution spatial characterization and an accurate picture of the shape and extent of the recharge plume of the MAR site. After 5 months of injection, a freshwater plume with an EC of 400–600 μS/cm had formed that extended 400 m in the W-E direction, 250 m in the N-S direction, and to a depth of 40 m below piezometric level. This study also provides correlations between ER values with different lithologies and groundwater EC values that can be used to support other studies.     

Geostatistical Evaluation of Spatial Variation Related to Groundwater Quality Database: Case Study for Arak Plain Aquifer,Iran

Geostatistical methods are one of the advanced techniques to interpolate groundwater quality data. Geostatistical interpolation techniques employ both the mathematical and the statistical properties of the measured points. Compiling the data distribution on spatial and temporal domain is of crucial importance in order to evaluate its quality and safety. The main purpose of this paper is to assess groundwater quality of Arak plain, Iran, by an unbiased interpolated method so called Kriging. Therefore, seven quality variables of Arak plain aquifer including TDS, SAR, EC, Na⁺, TH, Cl^?, and SO₄ ^2? have been analyzed, studied, and interpreted statistically and geostatistically. Utilized data in this study were collected from 97 water well samples in Arak plain, in 2012. After normalizing data, variogram as a geostatistical tool for defining spatial regression was calculated and experimental variograms have been plotted by GS⁺ software, then the best theoretical model was fitted to each variogram based on minimum RSS error. Cross validation was used to determine the accuracy of the estimated data. The uncertainty of the method could be well assessed via this method since the method not only gave the average error (around 0 in this study) but also gave the standard deviation of the estimations. Therefore, more than 3800 points were estimated by ordinary Kriging algorithm in places which have not been sampled. Finally, estimation maps of groundwater quality were prepared and map of estimation variance, EV, has been presented to assess the quality of estimation in each estimated point. Results showed that the Kriging method is more accurate than the traditional interpolation algorithms not honoring the spatial properties of the database.     

Nutrient dynamics between sediment and overlying water in the inner part of Izmir Bay, Eastern Aegean

In this article a concept is described in order to predict and map the occurrence of benthic communities within and near the German Exclusive Economic Zone (EEZ) of the North Sea. The approach consists of two work steps: (1) geostatistical analysis of abiotic measurement data and (2) calculation of benthic provinces by means of Classification and Regression Trees (CART) and GIS-techniques. From bottom water measurements on salinity, temperature, silicate and nutrients as well as from punctual data on grain size ranges (0–20, 20–63, 63–2,000 μ) raster maps were calculated by use of geostatistical methods. At first the autocorrelation structure was examined and modelled with help of variogram analysis. The resulting variogram models were then used to calculate raster maps by applying ordinary kriging procedures. After intersecting these raster maps with punctual data on eight benthic communities a decision tree was derived to predict the occurrence of these communities within the study area. Since such a CART tree corresponds to a hierarchically ordered set of decision rules it was applied to the geostatistically estimated raster data to predict benthic habitats within and near the EEZ.     

Historical variation of elements with respect to different geochemical fractions in recent sediments from Pigeon Lake, Alberta, Canada

Geochemical analysis of elements and organic matter were conducted on vertical profiles of the recent sediments from Pigeon Lake, Alberta, Canada, to determine historical variations in elemental content of the sediments as related to their geochemical fractions. The elements are grouped according to their affinity with different geochemical fractions, by using cluster analysis and sequential extraction experiments. As a result, four elemental fractions were identified: clastic mineral detritus; carbonate; organic; and elements that show less similarity to the previous groups perhaps due to anthropogenic input or the influence of other fractions, such as oxyhydroxides. Following the identification of geochemical fractions in the sediments, a three-step normalizing method was applied using parameters that represent each geochemical fraction. These normalizing techniques appear to be important in verifying whether the variation of elements is indeed the result of anthropogenic and/or natural activities. The normalized data are correlated with the recent history of human activity and natural events near Pigeon Lake. Given the age of the lake sediments, this correlation indicates that the depth profiles of elements after being normalized to the organic and carbonate fractions reflect the variation of detrital input into the lake. However, the former mainly corresponds to the coarse-grained clastic minerals originating from high-energy erosion as the result of natural events (e.g., flooding), whereas the latter corresponds to the low-energy erosion of the fine particles (enriched in lithophile elements) due to deforestation in the drainage basin. Normalizing to the clastic mineral detritus fraction results in the increase of heavy metals in the uppermost part of the sediment profiles, which coincides with industrial activities during the past two decades in central Alberta. However, the concentration of these elements is negligible, as compared to the quantities released by geogenic processes (erosion).     

Methylated arsenic, antimony and tin species in soils

Methylated species of antimony, arsenic and tin were examined in urban soils of the Ruhr basin, near the cities of Duisburg and Essen, Germany. The main aim of this study was to investigate the occurrence of mono-, di- and trimethylated species of these elements in urban soils. The influence of historical and present land use upon the species content was examined. The distribution of inorganic As, Sb and Sn and their methylated species along the profile depth was investigated. As, Sb and Sn speciation was performed by pH-gradient hydride generation purge and trap gas chromatography, followed by inductively-coupled plasma mass spectrometry (HG-PT-GC/ICP-MS). Species' structures were confirmed by GC-EI/MS-ICP-MS. Monomethylated Sb and As were the dominant species detected: the concentration of these metal(loid) species varied between <0.07-56 microg kg(-1) per dry mass. All dimethylated species and monomethyltin concentrations were between <0.01-7.6 microg kg(-1) per dry mass, and for the trimethylated species of all examined elements, concentrations between <0.001-0.63 microg kg(-1) per dry mass were detected. The highest organometal(loid) concentrations were observed in agricultural soils and garden soils; lower concentrations were found in the soils of abandoned industrial sites (wasteland, primary forest and grassland) and a flood plain soil of the Rhine. This result can be ascribed to both the cultivation and the increased biological activity of the agricultural soils, and the generally higher contamination, the disturbed structure and the artificial substrates (deposits from industrial sources) of the abandoned industrial soils. Due to periodical sedimentation, the flood plain profile was the only one where no depth dependence of organometal(loid) species concentration was detected. The other soil profiles showed a decrease of species content with increasing depth; this was particularly noticeable in soils with a clear change from a horizon with an organic character towards a mineral horizon, i.e. decreasing vitality from profile top to bottom. It is not as yet clear whether the organometal(loid) species are formed in the mineral horizons of the profiles or whether they are displaced from the organic, biologically-active horizons towards the mineral horizons. Field studies revealed that soil parameters like pH, water content or temperature did not correlate significantly with the degree of biomethylation observed. In contrast to the lower in vitro biomethylation efficiency of Sb vs. As in microbial incubations, we consistently detected higher proportions of transformed Sb compounds in situ in soil samples. These data may indicate a need to re-examine the currently accepted model of Sb biogeochemical cycling in the real environment.     

Assessment of polycyclic aromatic hydrocarbon influx and sediment contamination in an urbanized estuary

Sediments from the Patos Lagoon Estuary in Southern Brazil and sludge from incoming effluents were assessed for the distribution of polycyclic aromatic hydrocarbons (PAHs). Shallow sediments adjoining the City of Rio Grande were found to be contaminated by PAHs mainly near urban effluent discharge, as well as in the port area. Effluents clustered into four groups according to their sources (sewage, industrial, runoff, and mixed), with each demonstrating different contributions of PAHs to the estuary. There was a predominance of runoff and mixed sources. Navigation activity was the second most important source of PAHs to sediments. The PAHs ratio identified the origin of these contaminants as essentially pyrolytic. The impact of PAHs as a result of uncontrolled disposal or accidental discharge of PAH-rich residues was suggested for several points. These points were primarily near gas stations and motor workshops. In about 30% of sampled sediments, the concentration of benzo[a]pyrene surpassed the Threshold Effects Level adopted for marine environments.     

Surface geophysical exploration: developing noninvasive tools to monitor past leaks around Hanford’s tank farms

A characterization program has been developed at Hanford to image past leaks in and around the underground storage tank facilities. The program is based on electrical resistivity, a geophysical technique that maps the distribution of electrical properties of the subsurface. The method was shown to be immediately successful in open areas devoid of underground metallic infrastructure, due to the large contrast in material properties between the highly saline waste and the dry sandy host environment. The results in these areas, confirmed by a limited number of boreholes, demonstrate a tendency for the lateral extent of the underground waste plume to remain within the approximate footprint of the disposal facility. In infrastructure-rich areas, such as tank farms, the conventional application of electrical resistivity using small point-source surface electrodes initially presented a challenge for the resistivity method. The method was then adapted to directly use the buried infrastructure, specifically the steel-cased wells that surround the tanks, as “long” electrodes for both transmission of electrical current and measurements of voltage. Overcoming the drawbacks of the long electrode method has been the focus of our work over the past 7 years. The drawbacks include low vertical resolution and limited lateral coverage. The lateral coverage issue has been improved by supplementing the long electrodes with surface electrodes in areas devoid of infrastructure. The vertical resolution has been increased by developing borehole electrode arrays that can fit within the small-diameter drive casing of a direct push rig. The evolution of the program has led to some exceptional advances in the application of geophysical methods, including logistical deployment of the technology in hazardous areas, development of parallel processing resistivity inversion algorithms, and adapting the processing tools to accommodate electrodes of all shapes and locations. The program is accompanied by a full set of quality assurance procedures that cover the layout of sensors, measurement strategies, and software enhancements while insuring the integrity of stored data. The data have been shown to be useful in identifying previously unknown contaminant sources and defining the footprint of precipitation recharge barriers to retard the movement of existing contamination.     

Comparison of stochastic and deterministic methods for mapping groundwater level spatial variability in sparsely monitored basins

In sparsely monitored basins, accurate mapping of the spatial variability of groundwater level requires the interpolation of scattered data. This paper presents a comparison of deterministic interpolation methods, i.e. inverse distance weight (IDW) and minimum curvature (MC), with stochastic methods, i.e. ordinary kriging (OK), universal kriging (UK) and kriging with Delaunay triangulation (DK). The study area is the Mires Basin of Mesara Valley in Crete (Greece). This sparsely sampled basin has limited groundwater resources which are vital for the island’s economy; spatial variations of the groundwater level are important for developing management and monitoring strategies. We evaluate the performance of the interpolation methods with respect to different statistical measures. The Spartan variogram family is applied for the first time to hydrological data and is shown to be optimal with respect to stochastic interpolation of this dataset. The three stochastic methods (OK, DK and UK) perform overall better than the deterministic counterparts (IDW and MC). DK, which is herein for the first time applied to hydrological data, yields the most accurate cross-validation estimate for the lowest value in the dataset. OK and UK lead to smooth isolevel contours, whilst DK and IDW generate more edges. The stochastic methods deliver estimates of prediction uncertainty which becomes highest near the southeastern border of the basin.     

Status of the Riparian Ecosystem in the Upper San Pedro River, Arizona: Application of an Assessment Model

A portion of Arizona’s San Pedro River is managed as a National Riparian Conservation Area but is potentially affected by ground-water withdrawals beyond the conservation area borders. We applied an assessment model to the Conservation Area as a basis for monitoring long-term changes in riparian ecosystem condition resulting from changes in river water availability, and collected multi-year data on a subset of the most sensitive bioindicators. The assessment model is based on nine vegetation bioindicators that are sensitive to changes in surface water or ground water. Site index scores allow for placement into one of three condition classes, each reflecting particular ranges for site hydrology and vegetation structure. We collected the bioindicator data at 26 sites distributed among 14 reaches that had similar stream flow hydrology (spatial flow intermittency) and geomorphology (channel sinuosity, flood-plain width). Overall, 39% of the riparian corridor fell within condition class 3 (the wettest condition), 55% in condition class 2, and 6% in the driest condition class. Condition class 3 reaches have high cover of herbaceous wetland plants (e.g., Juncus and Schoenoplectus spp.) along the perennial stream channel and dense, multi-aged Populus-Salix woodlands in the flood plain, sustained by shallow ground water in the stream alluvium. In condition class 2, intermittent stream flows result in low cover of streamside wetland herbs, but Populus-Salix remain abundant in the flood plain. Perennial wetland plants are absent from condition class 1, reflecting highly intermittent stream flows; the flood plain is vegetated by Tamarixa small tree that tolerates the deep and fluctuating ground water levels that typify this reach type. Abundance of herbaceous wetland plants and growth rate of Salix gooddingii varied between years with different stream flow rates, indicating utility of these measures for tracking short-term responses to hydrologic change. Repeat measurement of all bioindicators will indicate long-term trends in hydro-vegetational condition.     

Risk assessment of contaminated sediments in river basins--theoretical considerations and pragmatic approach

In order to ensure chemical quality of river basins, measures as required in the Water Framework Directive will also have to address contaminated sediments with subsequent monitoring of their successful application. Financial resources need to be allocated as part of a river basin plan to those contaminated sites that pose the biggest risk to the river basin. In order to differentiate between areas with elevated contaminant levels ("areas of concern") and those sites from which contaminated sediments can become resuspended and transported with the river, affecting the water phase and downstream sites in the catchment ("areas of risk"), the dynamics of sediment and suspended matter need to come into focus. Hydrological data have to be combined with concentration of suspended matter and its contaminant concentration to allow assessment of particle bound contaminant load. Each of these kinds of data, however, are subject to uncertainties--e.g. due to natural variability, heterogeneity of the matrix, challenges during sampling and chemical analyses of the suspended matter. Considering these uncertainties throughout the traceability chains of data collection, use of different lines of evidence and tools like fuzzy logic will increase the confidence of potentially costly management decisions.     

Mercury in bottom sediments of the Amur River, its flood-plain lakes and estuary, Eastern Siberia

Mercury (Hg) is an element of a special concern in the Amur River basin, where numerous cinnabar deposits and manifestations have been prospected. Moreover, the territory is under heavy anthropogenic pressure due to intensive economic development that includes activities accompanied by noticeable emissions of Hg to the environment through poor waste management practices and accidental emergency discharges. Yet, information on Hg distribution and behavior in this region is scarce and inadequate. In order to evaluate Hg levels and fate in this vast territory, surveys of river, lake, and estuarine bottom sediments, as integral indicators of environmental status, were carried out in 1990, 1991, 1997, and 2004. The results showed the following: (1) stagnation of the Russian economy in the 1990s has resulted in a noticeable decrease of the Hg content in the Amur River sediments to the basin pristine level of about 0.05 mg kg^???1; (2) Hg distribution in the sediment depth proves the element redox-dependent behavior; (3) in some cases, Hg enrichment may be related to the long-term anthropogenic emission; (4) Hg concentration in bottom sediments was found to increase in the following order—the Amur River mouth, the estuary, and the Sea of Okhotsk, showing the weakly non-conservative Hg behavior during estuarine water mixing.     

Levels, sources and spatiotemporal variation of nutrients and micropollutants in small streams of a Mediterranean River basin

In this study, nutrients, trace metals and priority pesticide compounds were investigated for the first time in water and sediment samples in streams of the Evrotas River basin (S.E. Greece) from 2006 to 2008. The most important sources of contamination were from the entry of pesticides and nutrients into surface waters and sediments as a result of the intensive agricultural activity as well as from the uncontrolled disposal of olive mill and citrus processing wastewaters. Aquatic risk assessment revealed that all insecticides detected showed high risk, suggesting adverse effects on the stream biota. Among the metals analyzed, Cr, Ni and Ba presented the highest concentrations in sediments, however, due to natural geological processes. Multivariate statistical techniques applied for data compression, exploration and interpretation proved to be useful tools for identifying the most critical pollutants affecting the surface water quality. The findings of this study suggest that the inclusion of streams with small catchment areas into WFD monitoring and assessment programs is essential, especially those of the Mediterranean region.     

Fractionation of heavy metals and phosphorus in suspended sediments of the Yamuna river,India

The Yamuna river is the largest tributary of the Ganges river system. It originates in the Himalayas and flows through a varied geological terrain encompassing a large basin area. Metals Fe, Mn, Pb, Zn, Cu in different chemical fractions of suspended sediments such as exchangeable, carbonates, Fe–Mn oxides, organics and residual fractions were studied. Phosphorus associated with different chemical forms are discussed. The metals are mostly associated with residual fractions in the sediments followed by organics, Fe–Mn oxides, exhangeable and carbonates. Intensive use of chemical fertilizers and pesticides in agriculture in the basin affects the high inorganic phosphorus content in sediments.     

Longitudinal heterogeneity of sediment characteristics during southwest monsoon season in hyper-eutrophic Krishnagiri reservoir,India

Krishnagiri reservoir is a hyper-eutrophicated reservoir located in Krishnagiri district which is one of the drought-prone districts in Tamil Nadu, India. The reservoir water is being used for various purposes such as irrigation, drinking, fish rearing, livestock rearing, and recreation. Since there is no an evidence of investigation on bottom sediments in Krishnagiri reservoir, the present study was carried out during southwest monsoon season in 2008. This study examined the physical and chemical characteristics of the bottom sediments such as composition, redox potential, moisture content, organic carbon, organic matter, total phosphorus, and total iron at 15 locations in the reservoir. Phosphorus fractionation study was carried out to find out different fractions such as loosely adsorbed phosphorous, iron and aluminium-bound phosphorus, calcium-bound phosphorous, and organic phosphorous. Results indicated that there was spatial variation in the composition of sediments and low values of redox potential. The significant positive correlation exists between the organic carbon and organic phosphorus concentration. The lacustrine zone of the reservoir showed high accumulation of total phosphorus and total iron when compared to riverine and transition zones. This study concludes an allogenic origin of majority of inorganic phosphorus in the reservoir during the study period and this might have been derived from the catchment during the erosion process. The high concentration of surface sediment phosphorus clearly indicates a greater threat of eutrophication in Krishnagiri reservoir.     

Post-remediation evaluation of a LNAPL site using electrical resistivity imaging

Present understanding of the earth's subsurface is most often derived from samples at discrete points (wells) and interpolations or models that interpret the space between these points. Electrical resistivity imaging techniques have produced an improved capability to map contaminants (especially NAPLs--NonAqueous Phase Liquids) away from traditional wells using actual field data. Electrical resistivity image data, confirmed by drilling, have demonstrated that LNAPL (Light NAPL--less dense than water, such as gasoline) contaminants exist outside of a delineated and remediated area in Golden, Oklahoma. The data also demonstrate that LNAPL exists between monitoring and remediation wells which indicate low contaminant levels when sampled. Additionally, the electrical images provided the drilling location with the highest concentration of hydrocarbon ever found on the site, even after two phases of remediation work had been performed, although the sampling protocols varied. The results indicate that current methods of post-remediation site characterization are inadequate for complete site characterization.     

Levels and distribution of PCDD/Fs, dl-PCBs, and organochlorine pesticides in sediments from the lower reaches of the Haihe River basin, China

The concentrations of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs), dioxin-like polychlorinated biphenyls (dl-PCBs), and organochlorine pesticides (OCPs) in surface sediments from the lower reaches of the Haihe River basin, northern China, were determined by high-resolution gas chromatograph–high-resolution mass spectrometer. The concentrations of 2,3,7,8-substituted PCDD/Fs, dl-PCBs, and total OCPs [sum of hexachlorobenzene, hexachlorocyclohexanes, and dichlorodiphenyltrichloroethane (DDTs)] in 17 sediment samples were in the ranges of 11.6–1,180,924 pg/g dry weight (dw), 18.7–50,017 pg/g dw, and 1.7–35,280 ng/g dw, respectively. The contamination levels in the samples varied significantly between the different sites. Abnormally high concentrations of PCDD/Fs, dl-PCBs, and some OCPs were found in sediments from the lower reaches of the main channel of the Haihe River and the Dagu Drainage River, which were attributed to the historical production of pentachlorophenol and other pesticides near these locations. High levels of DDTs were detected in the Yongding New River sediment, which were likely to have originated from the discharge of wastewater from a dicofol factory upstream. In samples taken from other sites, the concentrations of these pollutants were at levels comparable to those documented in other areas of China. This preliminary investigation suggests that historical pesticide production in the Haihe River basin has contributed significantly to the contamination of this aquatic ecosystem and that further attention to this issue is warranted.     

Nature of distribution of mercury in the sediments of the river Yamuna (tributary of the Ganges), India

Suspended Particulate Matter (SPM), surface (bed sediments) and short length cores of sediments collected from the largest tributary of the river Ganges, namely the river Yamuna, were analysed for total mercury as well as its fractionation in various size and chemical sites in the sediments following standard procedures. Also, attempts were made to determine the vertical distribution in sediments in relation to the recent timescale of a few decades. Our observations indicate that the SPM in general showed higher levels of total mercury compared to the surface sediments while at places the enhancement could be by a factor of 10, say around 25 microg g(-1) in the downstream region that integrates the industrial midstream and agricultural downstream terrain near its confluence with the Ganges. Surface sediments in the upstream direction near the Himalayan foothills and SPM in the lower reaches showed significant high Index of Geoaccumulation (Igeo) as defined by Müller. Size fractionation studies indicate that the finer fraction preferentially showed higher levels of mercury while in the lower reaches of the river, the total mercury is equitably distributed among all size fractions. The proportion of the residual fraction of mercury in relation to mobile fractions, in general decreases downstream towards its confluence with the Ganges river. In sediment cores, the vertical distribution show systematic peaks of mercury indicating that addition of this toxic metal to the aquatic system is in direct proportion to the increase in various types of human activities such as thermal power plants, land use changes (urbanisation) in the midstream region and intensive fertiliser application in lower reaches of this vast river basin.     

Landscape changes as a factor affecting the course and consequences of extreme floods in the Otava river basin, Czech Republic

The paper presents the analysis of anthropogenical modifications of the landscape in relation to the course and consequences of floods. The research was conducted in the Otava river basin which represents the core zone of the extreme flood in August 2002 in Central Europe. The analysis was focused on the key indicators of landscape modification potentially affecting the runoff process - the long-term changes of land-use, changes of land cover structure, land drainage, historical shortening of the river network and the modifications of streams and floodplains. The information on intensity and spatial distribution of modifications was derived from different data sources - historical maps, available GIS data, remote sensing and field mapping. The results revealed a high level of spatial diversity of anthropogenical modifications in different parts of the river basin. The intensive modifications in most of indicators were concentrated in the lowland region of the river basin due to its agricultural use; however important changes were also recorded in the headwater region of the basin. The high spatial diversity of the modifications may result in their varying effect on the course and consequences of floods in different parts of the river basin. This effect is demonstrated by the cluster analysis based on the matrix of indicators of stream and floodplain modification, physiogeographical characteristics and geomorphological evidences of the flood in August 2002, derived from the individual thematic layers using GIS.     

Background Levels of Metals in Soils, McMurdo Station, Antarctica

McMurdo Station is the largest research station in Antarctica, with a population that ranges each year from 250 to 1200 people. Because of its size and 40-year history of use, a number of locations around the station have become contaminated with wastes. Soils and sediments in these areas have been shown to contain elevated levels of petroleum-related products, PCBs, other organics, and metals. While some remedial investigations have been conducted, background levels of metals in soils have not been determined. This paper reports on background levels of metals in a natural basalt-derived soil (gray soil) and scoria (soft porous rock used as fill, red soil) near McMurdo Station using two fundamentally different analytical procedures, concentrated acid extraction/analysis and total metals. These data facilitate determining the extent and levels of metal contamination near McMurdo Station and provide reference levels of metals for comparison with existing and future remediation data. There were statistically significant differences between metals concentrations in both gray and red soils, and no correlations between the level of extracted versus total metal. Generally, only a small fraction of a metal was extractable.     

Distribution, accumulation and mobility of mercury in superficial sediment samples from Tianjin, northern China

Seventeen sediment samples at three representative sites (the Yuqiao Reservoir, the Haihe River and the Haihe River Estuary) in Tianjin, northern China, were analyzed to investigate the pollution status, accumulation and mobility of mercury (Hg). The results show that the Haihe River has to be considered as moderate to strongly contaminated with Hg (2 < mean I(geo) = 2.35 < 3), where total Hg contents were ca. 3 to 24 orders of magnitude greater than the regional background value. The sediments collected near a coal-fired power plant in an urban area were found to have very high Hg concentrations, which were significantly related to Hg emissions from coal-fired utility boilers. In the Yuqiao Reservoir, the surface sediments have to be considered as unpolluted with Hg (mean I(geo) = -0.05 < 0) and the Haihe River Estuary sediments have to be considered as unpolluted to moderately polluted with Hg (0 < mean I(geo) = 0.18 < 1). Sediment-bound Hg in the Yuqiao Reservoir and the Haihe River Estuary was found to be predominantly associated with the organo-chelated phase of the sediment (38.3% and 50.5%, respectively). However, unlike the Yuqiao Reservoir and the Haihe River Estuary, Hg in the Haihe River sediments existed mainly as sulfide Hg and elemental Hg, which accounted for 54.2% and 30.7% of total Hg, respectively. The availability of this element seemed restricted. The majority of Hg contamination in the Haihe River sediments had been attributed to historic and modern atmospheric deposition and Hg released from the Haihe River sediments didn't seem to be an important pollutant pathway into the Haihe River Estuary. The results provide new insights into Hg contamination in this region.