Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods

Abstract

The resolution of the freshwater and saline water aquifers in a coastal terrain (Mahanadi Basin, India) is updated. We analysed electrical borehole log data at four sites and compared the water resistivity regime of the freshwater and saline water zones obtained from electrical borehole logging, with the resistivity regime obtained by interpreting vertical electrical sounding (VES) data. The multilayer VES data interpretation is modified to a simple model, containing only the freshwater zone and the saline water zone. The composite geophysical parameters of the freshwater and saline water zones, in particular the resistivity and longitudinal unit conductance regime, are identified. The resolution obtained from the composite geophysical data analyses is very clear and convincing. The composite longitudinal unit conductance regime of the saline water zones is very high compared to that of the freshwater zones. This makes the identification of the two aquifers easy and increases its reliability. A technique which enables analysis of composite geophysical data of freshwater and saline water zones at VES sites in the vicinity of the borehole log sites is proposed. The significance of longitudinal unit conductance in resolving the freshwater and saline water aquifers is illustrated graphically. The proposed technique is validated by correlating the longitudinal unit conductance and resistivity with the total dissolved solids. The efficiency of the technique is validated by carrying out discriminant function analysis.

Citation Hodlur, G. K., Dhakate, R., Sirisha, T. & Panaskar, D. B. (2010) Resolution of freshwater and saline water aquifers by composite geophysical data analysis methods. Hydrol. Sci. J. 55(3), 414–434.     

Impact of the water salinity on the hydraulic conductivity of fen peat

Coastal peatlands represent an interface between marine and terrestrial ecosystems; their hydrology is affected by salt and fresh water inflow alike. Previous studies on bog peat have shown that pore water salinity can have an impact on the saturated hydraulic conductivity (K_s) of peat because of chemical pore dilation effects. In this study, we aimed at quantifying the impact of higher salinities (up to 3.5% NaCl) on K_s of fen peat. Two experiments employing a constant‐head upward‐flow permeameter and differing in measurement and salinity change duration were conducted. Additionally, a third experiment to determine the impact of water salinity on the release of dissolved organic carbon (DOC) of the studied peat type was carried out. The results show a decrease of K_s with time, which does not depend on the water salinity but is differently shaped for different peat types. We assume pore clogging due to a conglomerate of physical, chemical, and biological processes, which rather depend on water movement rate and time than on water salinity. However, an increased water salinity did increase the DOC release. We conclude that salinity‐dependent behaviour of K_s is a function of peat chemistry and that for some peat types, salinity may only affect the DOC release without having a pronounced impact on water flow.     

内蒙古达里诺尔湖流域地表水和地下水环境同位素特征及补给关系

为了探明达里诺尔湖流域地表水与地下水的氢(H)、氧(O)同位素的变化特征及相互补给关系,于2013年对达里诺尔湖及其周围的河水、井水、泉水中H、O同位素进行了取样分析,并结合总溶解性固体悬浮物(TDS)和区域水文地质对达里诺尔湖流域的补给关系进行讨论分析.结果表明:1)河水、泉水、井水中H、O同位素的值基本落在全球雨水线上,湖水H、O同位素落在全球雨水线的右下方,说明河水、井水、泉水没有发生蒸发分馏,而湖水则发生较大程度的蒸发分馏;对达里诺尔湖流域地表水与地下水的H、O同位素进行回归模拟,得出该区域的蒸发趋势线方程:δD=4.8753δ18O-20.139(n=32,R2=0.9968).蒸发线表明,这些水样具有相同水源的特征.2)从实地考察发现,泉水补给河水,泉水和河水补给湖水,同时井水、泉水和河水有相似的δD、δ18O和TDS值,且不随季节变化而变化,推断达里诺尔湖附近地下水补给湖水;区域水文地质条件亦证明达里诺尔湖周边地下水补给湖水.     

Pathlines of pollutant migration in groundwater of the Visakhapatnam urban area,India

The groundwater regime in the north‐western part of the Visakhapatnam urban area was polluted as early as 1981 by discharge of untreated industrial effluent from a Hindustan Polymers Limited (HPL) plant. A total dissolved solids (TDS) concentration of surficial effluent up to 6500 mg/l and of groundwater in the range 3000–4200 mg/l has been reported in the environs of the HPL plant during May 1992. Groundwater occurs under water table conditions. The groundwater flow model was simulated as a single layer aquifer, using MODFLOW and FLOWPATH computer models. Aquifer parameters were estimated and the value of effective porosity assumed. The water table configuration was mapped in May 1981 using 33 observation wells. The calibrated hydraulic head distribution was used to compute the velocity field using the effective porosity values. The flow paths of groundwater migration from the source locations up to 2002 AD were predicted. Pathlines of particles in the groundwater indicated a predominant north‐east and south‐west migration of groundwater pollution in the area. Copyright © 1999 John Wiley & Sons, Ltd.     

西南地区前兆观测异常特征和地震短期预测方法研究

通过对西南地区1970～2001年39组52次6级以上地震及四川部分5级以上地震震例资料的剖析,进一步研究西南地区各类前兆中、短期异常判据和指标;研究西南地区强震孕育前兆异常由中期向短期(临震)过渡、演化的标志和判据并作物理解释;并进行了西南地区前兆短期预测方法及强震短期预测时、空、强指标和判据的研究和西南地区地下流体新观测方法及灵敏测项的研究。     

Investigation of ionospheric TEC changes related to the 2008 Wenchuan earthquake based on statistic analysis and signal detection

Ionospheric TEC (total electron content) time series are derived from GPS measurements at 13 stations around the epicenter of the 2008 Wenchuan earthquake. Defining anomaly bounds for a sliding window by quartile and 2-standard deviation of TEC values, this paper analyzed the characteristics of ionospheric changes before and after the destructive event. The Neyman-Pearson signal detection method is employed to compute the probabilities of TEC abnormalities. Result shows that one week before the Wenchuan earthquake, ionospheric TEC over the epicenter and its vicinities displays obvious abnormal disturbances, most of which are positive anomalies. The largest TEC abnormal changes appeared on May 9, three days prior to the seismic event. Signal detection shows that the largest possibility of TEC abnormity on May 9 is 50.74%, indicating that ionospheric abnormities three days before the main shock are likely related to the preparation process of the M_S8.0 Wenchuan earthquake.     

Total waterborne carbon export and DOC composition from ten nested subarctic peatland catchments—importance of peatland cover,groundwater influence,and inter‐annual variability of precipitation patterns

Waterborne carbon (C) export from terrestrial ecosystems is a potentially important flux for the net catchment C balance and links the biogeochemical C cycling of terrestrial ecosystems to their downstream aquatic ecosystems. We have monitored hydrology and stream chemistry over 3 years in ten nested catchments (0.6–15.1 km²) with variable peatland cover (0%–22%) and groundwater influence in subarctic Sweden. Total waterborne C export, including dissolved and particulate organic carbon (DOC and POC) and dissolved inorganic carbon (DIC), ranged between 2.8 and 7.3 g m^–2 year^–1, representing ~10%–30% of catchment net ecosystem exchange of CO₂. Several characteristics of catchment waterborne C export were affected by interacting effects of peatland cover and groundwater influence, including magnitude and timing, partitioning into DOC, POC, and DIC and chemical composition of the exported DOC. Waterborne C export was greater during the wetter years, equivalent to an average change in export of ~2 g m^–2 year^–1 per 100 mm of precipitation. Wetter years led to a greater relative increase in DIC export than DOC export due to an inferred relative shift in dominance from shallow organic flow pathways to groundwater sources. Indices of DOC composition (SUVA₂₅₄ and a₂₅₀/a₃₆₅) indicated that DOC aromaticity and average molecular weight increased with catchment peatland cover and decreased with increased groundwater influence. Our results provide examples on how waterborne C export and DOC composition might be affected by climate change. Copyright © 2012 John Wiley & Sons, Ltd.     

High-resolution OSL chronology of a well-preserved inland dune in the Lys valley (Sint-Martens-Latem,NW Belgium)

The inland dune “Molenberg” (Lys valley, NW Belgium) is an exceptionally well-preserved dune of more than 10 m high. To determine the phases of dune formation and identify possible levels for human occupation, this dune was sampled at closely spaced vertical intervals for Optically Stimulated Luminescence (OSL) dating (every 50 cm) and grain-size analysis (intervals of 10–25 cm) over its entire depth, including the top of the underlying sediments. The resulting ages show dune formation during the late glacial period (GI-1 and GS-1), with a probable extensive deposition during the Allerød (GI-1a-c) and Younger Dryas (GS-1). Interestingly, the OSL dataset for the main body of the dune shows significant cyclical age inversions that may be caused by reworking of sediment causing incomplete bleaching of sand grains. This may further indicate the potential of high-sampling resolution OSL to contribute to the reconstruction of complex dune dynamics.     

Biogeochemical and ecological functioning of the low-salinity water lenses in the region of the Rhone River freshwater influence,NW Mediterranean Sea

A quasi-synoptic survey of a eulerian small grid was carried out twice during mid-spring 2002 in the Gulf of Lions, NW Mediterranean Sea. Analysis of hydrological core parameters reveal for the first time, in the region of freshwater influence (ROFI) of the Rhone River, the presence of low-salinity water (LSW) lenses. The present work details the biogeochemical and ecological functioning of the two LSW lenses from a combined analysis of nutrients and organic matter content, taxonomic assemblages of phytoplankton, primary productivity measurements and nitrogen regeneration fluxes. During the first survey, the lens observed is only detached in part from the Rhone River plume and is considered as a juvenile lens. In contrast, the second lens is totally detached from the plume forming a confined 3D structure with a large vertical development and is considered as having a more advanced maturity. A second survey, 4 days later, provided the opportunity to propose a complete sequence of ecological functioning of the LSW lenses, from their formation to a late stage of maturity just before dispersion. Nitrate contents and dissolved organic matter remain in high concentrations during the juvenile stages while the little available phosphate is rapidly exhausted. In such, an unbalanced-nutrient environment picoplankton appear to out-compete bacterioplankton for phosphate and other resources such as ammonium. In turn, the dominance of such prokaryotic assemblages could have involved the rapid development of microzooplankton. The sudden increase in phosphate observed in a more advanced stage of lens maturity is attributed to intense P-regeneration driven by the large abundance of microzooplankton. This top-down control does not seem to enable the prokaryotic assemblages to bloom again but the high concentrations of phosphate and nitrate favour the development of larger phytoplankton. These autotrophic communities rapidly exhaust nutrients and then decline in the confined ecosystem of the LSW lens that subsequently evolves towards oligotrophic conditions. In the observed last stage, small-sized phytoplankton again dominates the microbial assemblages and the trophic regime is mainly based on regenerated production. The long lifetime of the LSW lens compared to that of the Rhone River plume enables a deep transformation of the initial characteristics of the Rhone River freshwaters both in terms of nutrients balance and microbial food webs, especially through trophic cascading effects. Hence, we show that the LSW lenses have the potential to transfer nutrients and organic matter offshore but the nature (inorganic vs. organic) of the dominant compound exported will depend on the maturity of the lens at the time of its dispersion.     

建设科学、有效的地震前兆监测台网提高捕捉地震前兆信息的能力--试论当前地震监测台网的三大弊病与三点建议

地震预报本质上是一门观测科学。本文分析了当前我国地震前兆监测台网存在的主要弊病,给出了建设具有中国特色、科学有效的地震前兆监测台网的建议。     

Temperature profiles in the earth of importance to deep electrical conductivity models

Deep in the Earth, the electrical conductivity of geological material is extremely dependent on temperature. The knowledge of temperature is thus essential for any interpretation of magnetotelluric data in projecting lithospheric structural models. The measured values of the terrestrial heat flow, radiogenic heat production and thermal conductivity of rocks allow the extrapolation of surface observations to a greater depth and the calculation of the temperature field within the lithosphere. Various methods of deep temperature calculations are presented and discussed. Characteristic geotherms are proposed for major tectonic provinces of Europe and it is shown that the existing temperatures on the crust-upper mantle boundary may vary in a broad interval of 350–1,000°C. The present work is completed with a survey of the temperature dependence of electrical conductivity for selected crustal and upper mantle rocks within the interval 200–1,000°C. It is shown how the knowledge of the temperature field can be used in the evaluation of the deep electrical conductivity pattern by converting the conductivity-versustemperature data into the conductivity-versus-depth data.     

Probabilistic assessment of earthquake recurrence in the January 26, 2001 earthquake region of Gujrat, India

Kutch region of Gujrat is one of the most seismic prone regions of India. Recently, it has been rocked by a large earthquake (M _w = 7.7) on January 26, 2001. The probabilities of occurrence of large earthquake (M≥6.0 and M≥5.0) in a specified interval of time for different elapsed times have been estimated on the basis of observed time-intervals between the large earthquakes (M≥6.0 and M≥5.0) using three probabilistic models, namely, Weibull, Gamma and Lognormal. The earthquakes of magnitude ≥5.0 covering about 180 years have been used for this analysis. However, the method of maximum likelihood estimation (MLE) has been applied for computation of earthquake hazard parameters. The mean interval of occurrence of earthquakes and standard deviation are estimated as 20.18 and 8.40 years for M≥5.0 and 36.32 and 12.49 years, for M≥6.0, respectively, for this region. For the earthquakes M≥5.0, the estimated cumulative probability reaches 0.8 after about 27 years for Lognormal and Gamma models and about 28 years for Weibull model while it reaches 0.9 after about 32 years for all the models. However, for the earthquakes M≥6.0, the estimated cumulative probability reaches 0.8 after about 47 years for all the models while it reaches 0.9 after about 53, 54 and 55 years for Weibull, Gamma and Lognormal model, respectively. The conditional probability also reaches about 0.8 to 0.9 for the time period of 28 to 40 years and 50 to 60 years for M≥5.0 and M≥6.0, respectively, for all the models. The probability of occurrence of an earthquake is very high between 28 to 42 years for the magnitudes ≥5.0 and between 47 to 55 years for the magnitudes ≥6.0, respectively, past from the last earthquake (2001).     

Suspended sediment and total dissolved solid yield patterns at the headwaters of Urumqi River,northwestern China: a comparison between glacial and non‐glacial catchments

In order to understand the differences in the suspended sediment and total dissolved solid (TDS) yield patterns between the glacial and non‐glacial catchments at the headwaters of Urumqi River, northwestern China, water samples were collected from a glacier catchment and an empty cirque catchment within the region, during three melting seasons from 2006 to 2008. These samples were analyzed to estimate suspended sediment and TDS concentrations, fluxes and erosion rates in the two adjoining catchments. There were remarked differences in suspended sediment and TDS yield patterns between the two catchments. Suspended sediment concentrations were controlled mainly by the sediment source, whereas TDS concentrations were primarily related to the hydrologic interaction with soil minerals. Generally, the glacial catchment had much higher suspended sediment and TDS yields, together with higher denudation rates, than the non‐glacial catchment. Overall, glacial catchment was mainly dominated by physical denudation process, whereas the non‐glacial catchment was jointly influenced by physical and chemical denudation processes. The observed differences in material delivery patterns were mainly controlled by the runoff source and the glacial processes. The melting periods of glacier and snow were typically the most important time for the suspended sediment and TDS yields. Meanwhile, episodic precipitation events could generate disproportionately large yields. Subglacial hydrology dynamics, glaciers pluck and grind processes could affect erodibility, and the large quantities of dust stored on the glacier surface provided additional sources for suspended sediment transport in the glacial catchment. These mechanisms imply that, in response to climate change, the catchment behaviour will be modified significantly in this region, in terms of material flux. Copyright © 2013 John Wiley & Sons, Ltd.     

The relationships between particulate and soluble alkaline phosphatase activities and the concentration of phosphorus dissolved in the seawater of Toulon Bay (NW Mediterranean)

The activities of particulate and soluble phosphatase were analyzed monthly for 1 year in the coastal ecosystems of the NW Mediterranean Sea. The mean contribution of the particulate activity increased from 56% at an MUF-P concentration of 30 μM to 77% at 0.04 μM. This particulate activity was negatively correlated with the DIP, DOP and TDP concentrations when the activities were related to the seawater volume, chlorophyll a or the protein concentration. The TDP correlations were highly significant (p: 0.001). The DOP correlations were significant (p: 0.04) and became highly significant (p: 0.009) at low DIP concentrations (<0.13 μM). The DIP correlations were significant (p: 0.04) only at low DOP concentrations (<0.18 μM). Thus, the effects of seawater DIP and DOP were found to be linked. The soluble activity exhibited distinct phosphatase fractions with high (0.5–29.5 μM) and low (0.02–2 μM) Km values, but none exhibited significant correlations with phosphorus compounds.     

Regionalization and dependence of coda Q on frequency and lapse time in the seismically active Peloritani region (northeastern Sicily,Italy)

The Peloritani region is one of the most seismically active regions in Italy and, consequently, the quantification of attenuation of the medium plays an important role for seismic risk evaluation. Moreover, it is necessary for the prediction of earth ground motion and future seismic source studies. An in depth analysis has been made here to understand the frequency and lapse time dependence of attenuation characteristics of the region by using the coda of local earthquakes. A regionalization is likewise performed in order to investigate the spatial variation of coda Q across the whole region. Finally, our results are jointly interpreted with those obtained from recently published 3D velocity tomographies for further insights.     

Conductivity structure of crust in the Tangshan seismic area and the possibility of exploring potential seismic sources by magnetotelluric method

It is shown by the result of digital magnetotelluric soundings in the Tangshan seismic area and its surrounding regions that the crust under the surface conductive sediments is divided into two layers,i. e., the resistive upper crust and the conductive lower crust. The upper crust wherein the Tangshan main shock and most of the aftershocks occurred is a convex lens-like body which is cut by faults at the east, south and west sides. The focus of the mainshock was located at the position of maximum thickness of the resistive upper crust while the spatial variation of Curie point isothermal surface and the deepest limit of the depths of aftershocks coincide with the downward depression of the bottom of the resistive upper crust. Thus, the Tangshan main shock and most of its aftershocks were related closely to the resistive upper crust from the view points of either vertical layering or lateral variations. And there were only a very few aftershocks in the conductive lower crust. The mechanical property of the rocks transforms from being brittle in the upper crust into ductile in the lower crust mainly due to the combination of different factors, e.g., increase of confining pressure, change in minerals, rise in temperature as well as stabilization of slips by high pore—pressure. A small amount of water and a rise in temperature may lead to a decrease of the electric resistivity within the rock while a change in the static pressure and mineral content within the rock causes very little change in the electric resistivity. Thus it is deduced that a resistive upper crust and a more conductive lower crust from the view points of either vertical or lateral variations are related to the brittle and ductile properties respectively. Hence it is possible that there is a relationship between the electric structure of the crust and zones of potential seismic hazard. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 354–363, 1991. The project is supported by the Chinese Joint Seismological Science Foundation, and the Swedish International Development Authority.     

Scale‐dependency of hydraulic conductivity: an example from Thorne Moor,a raised mire in South Yorkshire,UK

The hydraulic conductivity (K) of many low permeability materials is strongly scale‐dependent. In raised mires and other types of peat deposit the effects of features such as abandoned infilled ditches, root holes and localized woody material, cause K to be heterogeneous and scale‐dependent. Despite this, field measurements are routinely made using auger hole (slug) tests at the scale of only a few tens of centimetres. Such measurements are locally valid, but where the regional subsurface movement of water through peat bogs is simulated using groundwater models, typically at the scale of hundreds of metres, they give rise to a systematic underestimate of flows and an overprediction of water table elevations. Until now, techniques to obtain values at a scale sufficiently large to include the effects of localized features of higher permeability have not been applied routinely. Research at Thorne Moor, a large raised mire, demonstrates that the K of peat varies over several orders of magnitude when measured at different scales, using a variety of techniques. Laboratory and auger hole tests cannot be relied upon to provide results that represent the hydraulic conductivity of large expanses of peatland. This has significant implications for the management and long‐term restoration of peatlands where both regional and local control of water levels is crucial. For groundwater models to be used successfully to plan such schemes, it is essential to apply the K values relevant to the scale of the simulation. This paper describes and tests novel techniques, using ditches, for the derivation of K at large scales which overcome many of the problems that have been identified with conventional techniques and are capable of producing estimates that are appropriate to the application of physically based regional flow models. Copyright © 2004 John Wiley & Sons, Ltd.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

2005年云南会泽、文山5.3级地震短临预报及再认识

2005年8月5日、8月13日云南会泽、文山分别发生5．3级地震。震前,前兆各学科出现不同程度的异常变化,异常台项达30多项,并呈现非线性上升趋势。据此,震前对这两次地震做出了短临预报。震后对两次地震的震前异常做了进一步分析研究,认为异常主要分布在距离震中20～180km范围内,集中出现在震前2个月。水温异常对发震地点有一定的指示意义。