Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

藏南卡鲁雄峰枪勇冰川新冰期冰川发育探讨

山地冰川对气候变化响应敏感,全新世以来,冰川呈加速退缩趋势。冰川地貌是冰川变化最直接的证据。枪勇冰川发源于卡鲁雄峰,位于喜马拉雅山中段北坡雨影区,雅鲁藏布江以南,冰川末端在海拔5000 m附近,冰川面积约为6.63 km2。枪勇沟从外至内共保存有4期冰碛垄,其中以"大枪勇错"外侧的新冰期冰碛垄最为高大,终碛垄高约30 m,且形态清晰。文章基于野外观察及谷歌影像（Google Earth）上冰碛垄的分布范围,恢复并计算了新冰期时枪勇冰川面积约为10.0 km2,较现代冰川扩张了1.5倍。采用VOLTA模型计算了现代冰川的体积为0.21±0.04 km3,平均厚度约为31.7 m。基于纵剖线模型模拟了新冰期冰川表面高程,新冰期冰川平均厚度约为58.5 m,冰川体积约0.59±0.12 km3,是现代冰川体积的2.8倍,体积减小较面积退缩更为剧烈。利用积累区面积比率法（Accumulation Area Ratio Method,简称AAR）,采用AAR值为0.7±0.05,基于现代冰川表面高程和模拟的新冰期冰川表面高程,计算现代冰川和新冰期冰川的物质平衡线高度（Equilibrium Line Altitude,简称ELA）分别为5780±70 m和5660±100 m,新冰期时平衡线高度降低了约120 m。而不利用模型恢复新冰期表面高程时,平衡线降低值为240 m,可能高估约120 m。新冰期时,枪勇冰川扩张受区域降温控制,冰川前进可能是对4.2 ka冷事件的响应。冰碛垄形成时代和大枪勇错孢粉数据表明,冰川可能在2.5 ka左右已经退缩,2.2 ka时退缩加剧。     

Heterogeneous characteristics of streambed saturated hydraulic conductivity of the Touchet River,south eastern Washington,USA

Traditionally a streambed is treated as a layer of uniform thickness and low saturated hydraulic conductivity (K) in surface‐ and ground‐water studies. Recent findings have shown a high level of spatial heterogeneity within a streambed and such heterogeneity directly affects surface‐ and ground‐water exchange and can have ecological implications for biogeochemical transformations, nutrient cycling, organic matter decomposition, and reproduction of gravel spawning fish. In this study a detailed field investigation of K was conducted in two selected sites in Touchet River, a typical salmon spawning stream in arid south eastern Washington, USA. In‐stream slug tests were conducted to determine K following the Bouwer and Rice method. For the upper and lower sites, each 50 m long and 9 m wide and roughly 20 m apart, a sampling grid of 5 m longitudinally and 3 m transversely was used. The slug tests were performed for each horizontal coordinate at 0·3–0·45, 0·6–0·75, 0·9–1·05 and 1·2–1·35 m depth intervals unless a shallower impenetrable obstruction was encountered. Additionally, water levels were measured to obtain vertical hydraulic gradient (VHG) between each two adjacent depth intervals. Results indicated that K ranged over three orders of magnitude at both the upper and lower sites and differed between the two sites. At the upper site, K did not differ significantly among different depth intervals based on nonparametric statistical tests for mean, median, and empirical cumulative distribution, but the spatial pattern of K varied among different depth intervals. At the lower site, K for the 0·3–0·45 m depth interval differed statistically from those at other depth intervals, and no similar spatial pattern was found among different depth intervals. Zones of upward and downward water flow based on VHG also varied among different depth intervals, reflecting the complexities of the water flow regime. Detailed characterization of the streambed as attempted in this study should be helpful in providing information on spatial variations of streambed hydraulic properties as well as surface‐ and ground‐water interaction. Copyright © 2009 John Wiley & Sons, Ltd.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

桶形基础沉贯室内模型试验研究

桶形基础采用负压沉贯法施工,其贯入阻力与压桩,打桩等施工方法显著不同。本文介绍了在粉土和粉质粘土地层中的一组模型试验情况,认为负压法施工可显著降低粉土的抗贯阻力,但在粉质粘土中减阻效果不明显。     

A multi-equation spatial econometric model,with application to EU manufacturing productivity growth

A multi-equation spatial econometric model is used to explain variations across EU regions in manufacturing productivity growth based on recent theoretical developments in urban economics and economic geography. The paper shows that temporal and spatial parameter homogeneity is an unrealistic assumption, contrary to what is typically assumed in the literature. Constraints are imposed on parameters across time periods and between core and peripheral regions of the EU, with the significant loss of fit providing overwhelming evidence of parameter heterogeneity, although the final model does highlight increasing returns to scale, which is a central feature of contemporary theory.      

Assessing landslide movements in volcanic islands using near-shore marine geophysical data: south Pico island,Azores

Marine geophysical data from around the submarine flanks of volcanic islands can potentially help to resolve whether large-scale instability of an edifice has been geologically recently active. We use geophysical data to investigate part of the coast of Pico Island of the Azores where, above sea-level, a major slump of Topo volcano has been interpreted previously from arcuate escarpments and a rugged irregular topography seaward of them. Multibeam echo-sounder data collected offshore of this feature show remarkably little evidence for slump fault movements in the island's submarine slope. Mid-slope benches, like those associated with the Hilina Slump of Kilauea, are absent. The high-resolution data extends onto the island's shelf, allowing us to evaluate evidence for continuing activity there. In particular, as the shelf's rock platform will have been last modified by surf erosion during the postglacial period of sea-level transgression, it provides a reference surface of intermediate age (7–19 ka) that can potentially reveal whether any movements occurred in pre-historic times. Where the arcuate escarpments are continued offshore, the modern seabed shows no bathymetric evidence for active faults where the shelf rock platform crops out in the multibeam data. Elsewhere, mobile shelf sediments could be disguising evidence for active faulting so we examined boomer profiles able to image the rock platform beneath them. The data reveal a platform that is steep (6.6°) compared with the dips of platforms that we have studied previously around the coast of adjacent Faial Island and steeper than the platform outside the proposed slump. This suggests that it was created by coastal erosion over a shorter period and hence is consistent with a younger age of the coastline. As with the multibeam data, where escarpments are continued offshore onto the shelf, the rock surface imaged with these boomer data also shows no clear evidence of major slump-related fault displacements. This study therefore illustrates how high-resolution boomer seismic and multibeam data could usefully contribute to hazard assessment of volcanic islands, by helping to evaluate areas with no historical movements. Explanations to reconcile the onshore and offshore data here are also put forward.     

APPLICATION OF THE CLUSTERING METHOD IN ANALYSING SHALLOW WATER MASSES AND MODIFIED WATER MASSES IN THE HUANGHAI SEA AND EAST CHINA SEA

The existing various T-S diagram methods of analysing water masses have a good effect on the ocean water mass analysis, but they have many limitations in analysing water masses in the shallow sea. The main point is that it is difficult to determine the core of original water type in the shallow sea area, so a large error may be carried into the determination of the boundaries of water mass. For example, when water masses are analyzed with the conventional analysis method of concentration mixture, we usually obtain a very     

Finite element analysis for a finite conductivity fracture in an infinite poroelastic medium

In the technology of oil recovery, oil production rate can be increased by generation of a vertical conductive fracture adjacent to the well-bore. In this paper the seepage flow and isothermal deformation in both the oil formation and the fracture are studied by modelling the formation as a two-dimensional infinite poroelastic medium and the conductive fracture as a one-dimensional poroelastic material, saturated by a one-phase compressible fluid. The plane strain condition is employed. Solutions for a growing conductive fracture and a stationary conductive fracture in the infinite medium are obtained by means of the finite element method based on a variational principle for the formation which can impose the governing equations of the fracture. Infinite elements are used outside the finite element domain. Numerical results indicate that the injection rate, the applied pressure and the crack mouth opening displacement at the well-bore oscillate during the propagation of the conductive fracture. The production rate of a well with the conductive fracture is compared with that of a well without the conductive fracture. Finally, a new definition of the conductivity coefficient for the conductive fracture is presented. Copyright © 1999 John Wiley & Sons, Ltd.