Evaluation of a probabilistic approach to simulation of alkalinity and electrical conductivity at a river bank filtration site

In river bank filtration, impurities present in the river water travel with the bank filtrate towards the pumping well. During this passage, certain types of impurities, such as turbidity, total coliform, and so forth, may get attenuated; however, it is interesting to note that some of the instant raw river water quality parameters, such as alkalinity and electrical conductivity, increase after the passage of water through the porous medium. This occurs because water, when passing through the soil pores, absorbs many of the solutes that cause an increase in alkalinity and electrical conductivity. Measurements at a river bank filtration site for a year showed that alkalinity of 116–32 mg l^?1 in river water increased to 222.4–159.9 mg l^?1 in the river bank filtered water. Likewise, the electrical conductivity increased from 280–131 μS cm^?1 to 462–409.6 μS cm^?1. This study uses a probabilistic approach for investigating the variation of alkalinity and electrical conductivity of source water that varies with the natural logarithm of the concentration of influent water. The probabilistic approach has the potential of being used in simulating the variation of alkalinity and electrical conductivity in river bank filtrate. Copyright © 2011 John Wiley & Sons, Ltd.     

Garnet crystal plasticity in the continental crust,new example from south Madagascar

Garnet (10 vol.%; pyrope contents 34–44 mol.%) hosted in quartzofeldspathic rocks within a large vertical shear zone of south Madagascar shows a strong grain‐size reduction (from a few cm to ～300 μm). Electron back‐scattered diffraction, transmission electron microscopy and scanning electron microscope imaging coupled with quantitative analysis of digitized images (PolyLX software) have been used in order to understand the deformation mechanisms associated with this grain‐size evolution. The garnet grain‐size reduction trend has been summarized in a typological evolution (from Type I to Type IV). Type I, the original porphyroblasts, form cm‐sized elongated grains that crystallized upon multiple nucleation and coalescence following biotite breakdown: biotite + sillimanite + quartz = garnet + alkali feldspar + rutile + melt. These large garnet grains contain quartz ribbons and sillimanite inclusions. Type I garnet is sheared along preferential planes (sillimanite layers, quartz ribbons and/or suitably oriented garnet crystallographic planes) producing highly elongated Type II garnet grains marked by a single crystallographic orientation. Further deformation leads to the development of a crystallographic misorientation, subgrains and new grains resulting in Type III garnet. Associated grain‐size reduction occurs via subgrain rotation recrystallization accompanied by fast diffusion‐assisted dislocation glide. This plastic deformation of garnet is associated with efficient recovery as shown by the very low dislocation densities (10¹⁰ m^?3 or lower). The rounded Type III garnet experiences rigid body rotation in fine‐grained matrix. In the highly deformed samples, the deformation mechanisms in garnet are grain‐size‐ and shape‐dependent: dislocation creep is dominant for the few large grains left (>1 mm; Type II garnet), rigid body rotation is typical for the smaller rounded grains (300 μm or less; Type III garnet) whereas diffusion creep may affect more elliptic garnet (Type IV garnet). The P–T conditions of garnet plasticity in the continental crust (≥950 °C; 11 kbar) have been identified using two‐feldspar thermometry and GASP conventional barometry. The garnet microstructural and deformation mechanisms evolution, coupled with grain‐size decrease in a fine‐grained steady‐state microstructure of quartz, alkali feldspar and plagioclase, suggests a separate mechanical evolution of garnet with respect to felsic minerals within the shear zone.     

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.     

The role of eddies in solute transport and recovery in rock fractures: Implication for groundwater remediation

A better understanding of solute transport and retention mechanism in rock fractures has been challenging due to difficulty in their direct observations in microscale rough‐walled fractures. Six representative troughs in a rough‐walled fracture were selected for microscale observations of eddy formation with increasing flow velocity and its effect on spatiotemporal changes of solute concentration. This experimental study was enabled by a microscale visualization technique of micro particle image velocimetry. With increasing flow velocity (Re ≤ 2.86), no eddies were generated, and solutes along the main streamlines transported rapidly, whereas those near the wall moved slowly. A larger amount of solutes remained trapped at all troughs at Re = 2.86 than Re < 1. For Re = 8.57, weak eddies started to be developed at the troughs on the lee side, which little contributed to overall solute flushing in the fracture. Accordingly, a large of amount of water was needed for solute flushing. The flow condition of 1 < Re < 10, before a full development of eddies, was least favourable in terms of time and amount of remediation fluid required to reach a target concentration. After large eddies were fully developed at troughs on the lee side for Re = 17.13, solutes were substantially reduced by eddies with less amount of water. Fully developed eddies were found to enhance solute transport and recovery, as opposed to a general consensus that eddies trap and delay solutes. Direct inflow into troughs on the stoss side also made a great contribution to solute flushing out of the troughs. This study indicates that fully developed eddies or strong inflows at troughs are highly possible to form for Re > 10 and this flow range could be favourable for efficient remediation.     

Late-Hercynian shearing during crystallization of granitoid magmas (Sila massif,southern Italy): regional implications

Shearing of regional extent, involving granitoids and underlying mid-crustal rocks of the Sila massif (Calabria, Italy), is analysed in this paper. The deformed granitoids are affected by a wide NNW-SSE oriented deformation zone, stretching for about 60 km, from the neighbourhood of Cecita Lake to Cropani village. Meso- and micro-structures in granitoids, close to the boundary with underlying migmatitic paragneiss, indicate that deformation developed from melt-present to solid-state conditions. Simultaneous tectonics and magmatism activated a plutonic accretionary process at mid-crustal levels. This took place at about 300 Ma and involved hybrid magmas with a dominat contribution from a mantle source. The deformation regime remained steady for a long time during magma crystallization and cooling in subsolidus conditions. The regional top-to-the-W sense of shear in the present geographic coordinates, recorded in the deformed granitoids, seems geometrically consistent with the coeval direction of maximum extension found in another sector of the southern Hercynian belt, suggesting the original position of the Sila basement in this context. Magmatic ativity ended with the intrusion of mafic and felsic magams affected by a very weak deformation, ongoing during the final strain increments of the late-Hercynian stage.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

河岸带蜿蜒性与植被密度对潜流驻留时间的复合效应

潜流驻留时间是反映潜流交换、溶质迁移转化、生物组成、生态过程的重要特征变量。应用双循环可控式河岸带模型，以NaCl为示踪剂，深入研究了河岸带蜿蜒性与植被密度对潜流驻留时间的复合效应；应用量纲分析和主成分分析法，探析了潜流驻留时间的关键影响因子及其作用程度和敏感性。结果表明：①河岸带植被密度对潜流驻留时间具有促进作用，但存在一定的阈值。②潜流驻留时间分布与河岸蜿蜒形态呈类镜像效应；潜流驻留时间在河岸带蜿蜒波不同位置差异明显，迎水面平均驻留时间约为背水面的1/2。③河岸带蜿蜒性与植被密度是影响潜流驻留时间的关键因子，总作用程度达91.07%。潜流驻留时间对蜿蜒性最敏感，对植被密度敏感性较弱，对复合因子的敏感性介于两者之间；迎水面内组合因子对潜流驻留时间具有共同促进效应，背水面内组合因子对潜流驻留时间具有抑制和促进双重效应。     

An unusual zone of seismic coupling in the Bonin arc: The 1972 Hachijo-Oki earthquakes and related seismicity

The 1972 February and December Hachijo-Oki earthquakes (M _s=7.3 and 7.4), in the northernmost part of the Izu-Bonin subduction zone, are the only major events (M _s>7.0) in the Bonin arc for the past 80 years. Relocation of the hypocenters, using one smaller event having a wellconstrained focal depth as a master event, shows that the depth of the February event is 10 km shallower than that of the December event. We have determined the rupture process for both events by minimizing the error in waveform between observed and synthetic seismograms. Although the number of available stations are limited, the depth range of the major energy release for the December event extends deeper than for the February one. The rupture propagated up-dip for both events. It is likely that the rupture zone of the two events overlapped, and that the December event ruptured the deeper part. This suggestion is consistent with the observation that the aftershock zones of both events overlap with that of the December event shifted landward. The waveforms of the December event have a smaller high frequency component than those of the February event, suggesting that the stress at the thrust zone became more uniform or reduced after the February event.No thrust type smaller event occurred near the rupture zone. Instead, theP-axes of smaller events are parallel to the dip of the slab and theirT-axes dip to the southwest. Focal depths of these events estimated byP-wave forward modeling are generally between 40–50 km and located beneath the thrust zone. We thus interpret them as the events within the Pacific slab near the zone ruptured by the two major events. The stress concentration around the rupture zone of the major events is suggested to have triggered these slab events. After the occurrence of the large events, the slab events are concentrated near the deeper portion of the rupture zone. These events may have been caused by the loading of the down-dip compressional stress near the down-dip end of the rupture zone due to the rupture. The occurrence of the doublet of large earthquakes and a number of down-dip compressional events beneath their rupture zones in a shallow portion of the subducting slab indicates an unusual zone of seismic coupling in the Bonin arc, most of which is seismically quiescent.     

Fractal reconstruction of sea-floor topography

Sea-floor bathymetric profiles exhibit features at many different scales of length; this suggests that they could be described as fractals. An algorithm interpolating a fractal line between points has been used to reconstruct bathymetric profiles from a few data points. In general, this fractal line has the same Fourier amplitude spectrum as real bathymetry, and, if the parameters of the interpolation are suitably chosen, it has a very similar appearance. The success of this fractal reconstruction algorithm for the sea-floor raises the possibility that it could be used to extrapolate, from data collected at one scale, the properties of the sea-floor at finer scales, and that similar techniques could be used to interpolate a surface between bathymetric profiles. The fractal character is a sign that the processes that shape the sea-floor are scale invariant and suggests that the renormalization group technique could be used to model these processes.