Magnetic resonance sounding applied to aquifer characterization

Magnetic resonance sounding (MRS) is distinguished from other geophysical tools used for ground water investigation by the fact that it measures a magnetic resonance signal generated directly from subsurface water molecules. An alternating current pulse energizes a wire loop on the ground surface and the MRS signal is generated; subsurface water is indicated, with a high degree of reliability, by nonzero amplitude readings. Measurements with varied pulse magnitudes then reveal the depth and thickness of water saturated layers. The hydraulic conductivity of aquifers can also be estimated using boreholes for calibration. MRS can be used for both predicting the yield of water supply wells and for interpolation between boreholes, thereby reducing the number of holes required for hydrogeological modeling. An example of the practical application of MRS combined with two-dimensional electrical imaging, in the Kerbernez and Kerien catchments area of France, demonstrates the efficiency of the technique.     

Fitting a stochastic partial differential equation to aquifer data

The steady state two dimensional groundwater flow equation with constant transmissivities was studied by Whittle in 1954 as a stochastic Laplace equation. He showed that the correlation function consisted of a modified Bessel function of the second kind, order 1, multiplied by its argument. This paper uses this pioneering work of Whittle to fit an aquifer head field to unequally spaced observations by maximum likelihood. Observational error is also included in the model. Both the isotropic and anisotropic cases are considered. The fitted field is then calculated on a two dimensional grid together with its standard deviation. The method is closely related to the use of two-dimensional splines for fitting surfaces to irregularly spaced observations.     

The response of idealized aquifer/river systems to climate change

Abstract

A study of the effect of changes in climate on aquifer storage and river recharge using a simple model of an idealized aquifer/river system shows the combined influence of aquifer properties and climate change scenario on the system response. The study shows that changes in the seasonal distribution of recharge may have a critical effect on low flows in rivers supported by baseflow. However, rivers supported by slowly responding aquifers may show a considerable delay in response to climate change allowing an opportunity for water resources planning over an extended period.     

Spatial measurements of stream baseflow,a relevant method for aquifer characterization and permeability evaluation. Application to a hard‐rock aquifer,the Oman ophiolite

The structure, functioning and hydrodynamic properties of aquifers can be determined from an analysis of the spatial variability of baseflow in the streams with which they are associated. Such analyses are based on simple low‐cost measurements. Through interpreting the hydrological profiles (Q = f(A)) it is possible to locate the aquifer(s) linked to the stream network and to determine the type of interrelated flow, i.e. whether the stream drains or feeds the aquifer. Using an analytical solution developed for situations with a positive linear relationship, i.e. where the baseflow increases linearly with increasing catchment size, it is also possible to estimate the permeability of the aquifer(s) concerned at catchment scale. Applied to the hard‐rock aquifers of the Oman ophiolite, this method shows that the ‘gabbro’ aquifer is more permeable than the ‘peridotite’ aquifer. As a consequence the streams drain the peridotites and ‘leak’ into the gabbro. The hydrological profiles within the peridotite are linear and positive, and indicate homogeneity in the hydrodynamic properties of these formations at the kilometre scale. The permeability of the peridotite is estimated at 5 · 10^?7 to 5 · 10^?8 m/s. Copyright © 2004 John Wiley & Sons, Ltd.     

Delineating alluvial aquifer heterogeneity using resistivity and GPR data

Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with >2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.     

Contribution of magnetic resonance sounding to aquifer characterization and recharge estimate in semiarid Niger

To improve the knowledge of the regionally important Continental Terminal 3 (CT3) aquifer in south-western Niger, fifteen magnetic resonance soundings (MRS) were carried out in December 2005 in the vicinity of wells and boreholes. The output MRS geophysical parameters, i.e. water content and decay constants versus depth, were compared to hydrogeological characteristics, i.e. water table depth, total porosity, specific yield and transmissivity estimated from direct measurements, pumping tests and transient groundwater modelling. The MRS-determined parameters were then used to estimate the rates of groundwater recharge.Contained in poorly consolidated Tertiary sandstones, the CT3 aquifer's water table has continuously risen by 4 m in total over the past four decades. Additionally, a significant portion of this increase has occurred in the past decade alone, with an annual rise now ranging between 0.1 and 0.3 m depending on the monitored well. Increase in groundwater recharge due to land clearance and deforestation explains this situation. According to previous estimations, the pre-clearing recharge ranged from 1 to 5 mm per year in 1950–60 s, while more recent recharge rates (1990s–2000s) range from 20 to 50 mm per year. These recharge values are directly affected by estimated aquifer specific yield value, while the spatial variation of rates of water table rise can be attributed to large scale hydrodynamic heterogeneities in the aquifer. However, few field measurements were available to confirm these assumptions.The main results of this study are: (1) The water table depth and aquifer transmissivity are estimated from MRS output parameters with an average accuracy of ± 10% and ± 9% respectively. (2) The MRS-determined water content is linked to both the total porosity and the specific yield of the aquifer, but no quantitative formulation can be proposed as yet. (3) Using the average MRS-determined water content over the investigated area, i.e. 13%, the groundwater recharge rates can be estimated to be ~ 2 mm per year in the 1950–1960s (pre-clearing period), and ~ 23 mm per year for the last decade. (4) The variations in specific yield and transmissivity cannot explain by themselves the spatial variability of the rise of the water table. (5) The ranges in transmissivity and water content obtained from MRS are more realistic than the groundwater modelling outputs. Therefore, MRS could be used to better constrain the aquifer parameters in groundwater modelling with a dense site network.Finally, this work illustrates how MRS can successfully improve characterisation and transient multi-year groundwater balance of commonly found sedimentary aquifers, particularly when integrated with well observations and pumping tests.     

The incentive program for fishermen to collect marine debris in Korea

A significant amount of sea-based marine debris is generated and deposited in Korean coastal waters. The Korean government has removed 66,000 tons of deposited marine debris since 1999. However, to fully address the problem, the generation of marine debris needs to be prevented. As an alternative, Korean government has initiated an incentive program for fishermen, who collect fishing gear or other marine debris while fishing. Although a considerable amount of marine debris is collected and removed from major fishing grounds, the program may not incentivize long-term sustainable behavior, resulting in further discarding of marine debris when the program is completed.     

提高地磁绝对观测数据精度的有效途径

经过长期观察发现用CZM-2型核旋仪(及其同类仪器)进行定点观测读数时,运用分量仪线圈“单向”旋转的方法进行分量观测,再加上做好一些其它辅助措施,能够充分消除人为因素和仪器本身误差所带来的影响,明显提高仪器的观测精度,达到使观测资料更加准确的目的。     

Using existing data to estimate aquifer properties, Great Lakes Region, USA

To determine specific storage and porosity, areally limited and time-consuming aquifer tests are frequently done. Hydrogeologic studies often do not have the resources to collect such data and rely on existing data sources for aquifer properties. An alternative tool for determining these aquifer properties is the analysis of earth tides. The objective of this study was to determine whether existing water-level and barometric-pressure data could be used to determine aquifer properties, such as porosity and specific storage, on a regional scale. In this study, national databases from the Great Lakes Region were queried for continuous records of groundwater-level and barometric-pressure data. Records from 37 selected wells were then analyzed for barometric efficiency and earth-tide responses. Specific-storage (S(s) ) and porosity values were determined, and the quality of the results were assessed with a measure of the "goodness of fit" (percent variance) of reconstruction of the response. Records from wells completed in several aquifer systems were analyzed with varying degrees of success. Aquifer S(s) values ranging from 5.9 × 10(-8) to 3.8 × 10(-6) /m were derived, with percent variance of reconstruction ranging from 1% to 78%. Comparisons with aquifer and laboratory testing of S(s) and porosity are favorable if the percent variance of reconstruction is above about 30%. Although the earth-tide-analysis method is not suitable for every situation, the S(s) and porosity of aquifers can, in many places, be estimated with existing water-level and barometric-pressure data or with data that are relatively inexpensive to collect.     

Numerical simulation of subsurface characterization methods: application to a natural aquifer analogue

Information from an outcrop is used as an analogue of a natural heterogeneous aquifer in order to provide an exhaustive data set of hydraulic properties. Based on this data, two commonly used borehole based investigation methods are simulated numerically. For a scenario of sparse sampling of the aquifer, the process of regionalization of the borehole hydraulic conductivity values is simulated by application of a deterministic interpolation approach and conditioned stochastic simulations. Comparison of the cumulative distributions of particle arrival times illustrates the effects of the sparse sampling, the properties of the individual investigation methods and the regionalization methods on the ability to predict flow and transport behaviour in the real system (i.e. the exhaustive data set).     

Statistic inversion of multi-zone transition probability models for aquifer characterization in alluvial fans

Understanding the heterogeneity arising from the complex architecture of sedimentary sequences in alluvial fans is challenging. This paper develops a statistical inverse framework in a multi-zone transition probability approach for characterizing the heterogeneity in alluvial fans. An analytical solution of the transition probability matrix is used to define the statistical relationships among different hydrofacies and their mean lengths, integral scales, and volumetric proportions. A statistical inversion is conducted to identify the multi-zone transition probability models and estimate the optimal statistical parameters using the modified Gauss–Newton–Levenberg–Marquardt method. The Jacobian matrix is computed by the sensitivity equation method, which results in an accurate inverse solution with quantification of parameter uncertainty. We use the Chaobai River alluvial fan in the Beijing Plain, China, as an example for elucidating the methodology of alluvial fan characterization. The alluvial fan is divided into three sediment zones. In each zone, the explicit mathematical formulations of the transition probability models are constructed with optimized different integral scales and volumetric proportions. The hydrofacies distributions in the three zones are simulated sequentially by the multi-zone transition probability-based indicator simulations. The result of this study provides the heterogeneous structure of the alluvial fan for further study of flow and transport simulations.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Influence of aquifer heterogeneity and return flow on pumping test data interpretation

Analytical solutions of drawdown in unconfined aquifers are widely applied for determining the specific yield, S_y, and the horizontal and the vertical hydraulic conductivity K_r and K_z, respectively. In many previous studies, estimates of S_y and K_z were observed to be highly variable and physically unrealistic. This has been attributed to the conceptualization of flow above the declining water table and aquifer heterogeneity in the applied models. We present the analysis of time-drawdown data from a pumping test instrumented with depth-differentiated observation piezometers arranged in clusters. Applying homogeneous anisotropic aquifer models in combination with nonlinear least squares parameter identification techniques, the data were analyzed in different groups: analysis of data from individual piezometer clusters and simultaneous analysis of the entire data set from all piezometer clusters (global analysis). From the cluster analyses, estimates of S_y and K_z exhibit large variances and depart from a priori estimates inferred from the hydrostratigraphy. Parameter estimates from the global analysis do not fall within the parameter bounds (minimum and maximum values) defined by the cluster analyses. While heterogeneity appears to be the important reason for large parameter variances, we discuss the influence of rarely considered aquifer return flow on drawdown and the inconsistent results from the cluster and global analyses. We corroborate our findings with data on hydraulic gradients, slug test data, and results from the application of a more realistic numerical flow model.     

Trench infiltration for managed aquifer recharge to permeable bedrock

Managed aquifer recharge to permeable bedrock is increasingly being utilized to enhance resources and maintain sustainable groundwater development practices. One such target is the Navajo Sandstone, an extensive regional aquifer located throughout the Colorado Plateau of the western United States. Spreading‐basin and bank‐filtration projects along the sandstone outcrop's western edge in southwestern Utah have recently been implemented to meet growth‐related water demands. This paper reports on a new cost‐effective surface‐infiltration technique utilizing trenches for enhancing managed aquifer recharge to permeable bedrock. A 48‐day infiltration trench experiment on outcropping Navajo Sandstone was conducted to evaluate this alternative surface‐spreading artificial recharge method. Final infiltration rates through the bottom of the trench were about 0·5 m/day. These infiltration rates were an order of magnitude higher than rates from a previous surface‐spreading experiment at the same site. The higher rates were likely caused by a combination of factors including the removal of lower permeability soil and surficial caliche deposits, access to open vertical sandstone fractures, a reduction in physical clogging associated with silt and biofilm layers, minimizing viscosity effects by maintaining isothermal conditions, minimizing chemical clogging caused by carbonate mineral precipitation associated with algal photosynthesis, and diminished gas clogging associated with trapped air and biogenic gases. This pilot study illustrates the viability of trench infiltration for enhancing surface spreading of managed aquifer recharge to permeable bedrock. Published in 2010 by John Wiley & Sons, Ltd.     

On the way to airborne gravelometry based on 3D spatial data derived from images

This paper presents basic tests to develop an airborne photogrammetric methodology that derives grain size characteristics of gravel bed rivers. The data acquisition was done using a lightweight action cam and a hand-held digital single lens reflex camera. Image processing comprised the structure from motion technique and multiview-stereo algorithms to obtain digital elevation models of non-cohesive gravel beds. Laboratory results indicate that the method accuracy is about four to six times lower than laser-scan data when based on action cam data. The accuracy of digital elevation models computed via photos taken by the reflex camera is almost of the same range as the laser data. Field experiments were done to test the performance of image based gravelometry against manual surface sampling. For this application the action cam was mounted to a low-cost quadrocopter, while the reflex camera was operated by hand. Results indicate that this combination has a high potential to generate data from which characteristic grain size parameters can be estimated.     

Prediction of groundwater‐induced flooding in a chalk aquifer for future climate change scenarios

Current climate change models for the southeast UK predict changing rainfall patterns, with increased incidence of extreme events. The chalk aquifer in the UK and northern France is susceptible to groundwater‐induced flooding under such conditions. In this methodological study we apply a frequency domain analysis approach to the chalk aquifer to derive a transfer function between effective rainfall and groundwater level from 7 years of monitoring data from the North Heath Barn site, near Brighton. The derived transfer function was calibrated and validated against monitoring data and then used to predict groundwater level for rainfall models for high, medium and low emission scenarios from the UKCP09 database. The derived transfer function is most closely comparable to the linear aquifer model, despite evidence for both matrix and fracture or karst water flow in the chalk, with transmissivity and unconfined storativity at the catchment scale of 1548 m² day^?1 and 1.6 × 10^?2. The application of the transfer function to UKCP09 rainfall data suggests that groundwater‐induced flooding may be about four times more frequent by 2040–2069 compared with 1961–1990 and seven times more frequent by 2070–2099. The model data also suggest an increase in the duration of groundwater minima relative to the reference period. Compared to deterministic modelling which requires detailed knowledge of aquifer heterogeneity and processes, the transfer function approach, although with limitations, is simpler, incorporating these factors into the analysis through frequency and phase coefficients, and thus may have the potential for groundwater risk assessment in other areas. Copyright © 2015 John Wiley & Sons, Ltd.     

基于VB/MATLAB的前兆资料小波分析程序的实现

在研究近年来国内前兆资料小波处理的基础上,开发出一种基于VB／MATLAB的地震前兆资料处理软件,可在VB应用程序里直接调用MATLAB小波、去噪等函数,使分析处理更方便、快捷,同时以实例说明其实现过程。