Capture-zone design in an aquifer influenced by cyclic fluctuations in hydraulic gradients

Design of a groundwater pumping and treatment system for a wood-treatment facility adjacent to the tidally influenced Fraser River estuary required the development of methodologies to account for cyclic variations in hydraulic gradients. Design of such systems must consider the effects of these cyclic fluctuations on the capture of dissolved-phase contaminants. When the period of the cyclic fluctuation is much less than the travel time of the dissolved contaminant from the source to the discharge point, the hydraulic-gradient variations resulting from these cycles can be ignored. Capture zones are then designed based on the average hydraulic gradient determined using filter techniques on continuous groundwater-level measurements. When the period of cyclic fluctuation in hydraulic gradient is near to or greater than the contaminant travel time, the resulting hydraulic-gradient variations cannot be ignored. In these instances, procedures are developed to account for these fluctuations in the capture-zone design. These include proper characterization of the groundwater regime, assessment of the average travel time and period of the cyclic fluctuations, and numerical techniques which allow accounting for the cyclic fluctuations in the design of the capture zone. Electronic Publication     

Identifying and quantifying urban recharge: a review

The sources of and pathways for groundwater recharge in urban areas are more numerous and complex than in rural environments. Buildings, roads, and other surface infrastructure combine with man-made drainage networks to change the pathways for precipitation. Some direct recharge is lost, but additional recharge can occur from storm drainage systems. Large amounts of water are imported into most cities for supply, distributed through underground pipes, and collected again in sewers or septic tanks. The leaks from these pipe networks often provide substantial recharge. Sources of recharge in urban areas are identified through piezometry, chemical signatures, and water balances. All three approaches have problems. Recharge is quantified either by individual components (direct recharge, water-mains leakage, septic tanks, etc.) or holistically. Working with individual components requires large amounts of data, much of which is uncertain and is likely to lead to large uncertainties in the final result. Recommended holistic approaches include the use of groundwater modelling and solute balances, where various types of data are integrated. Urban recharge remains an under-researched topic, with few high-quality case studies reported in the literature. Electronic Publication     

An experimental study of grain scale melt segregation mechanisms in two common crustal rock types

Creation of pathways for melt to migrate from its source is the necessary first step for transport of magma to the upper crust. To test the role of different dehydration‐melting reactions in the development of permeability during partial melting and deformation in the crust, we experimentally deformed two common crustal rock types. A muscovite‐biotite metapelite and a biotite gneiss were deformed at conditions below, at and above their fluid‐absent solidus. For the metapelite, temperatures ranged between 650 and 800 °C at P_c=700 MPa to investigate the muscovite‐dehydration melting reaction. For the biotite gneiss, temperatures ranged between 850 and 950 °C at P_c=1000 MPa to explore biotite dehydration‐melting under lower crustal conditions. Deformation for both sets of experiments was performed at the same strain rate (ε.) 1.37×10^?5 s^?1. In the presence of deformation, the positive ΔV and associated high dilational strain of the muscovite dehydration‐melting reaction produces an increase in melt pore pressure with partial melting of the metapelite. In contrast, the biotite dehydration‐melting reaction is not associated with a large dilational strain and during deformation and partial melting of the biotite gneiss melt pore pressure builds more gradually. Due to the different rates in pore pressure increase, melt‐enhanced deformation microstructures reflect the different dehydration melting reactions themselves. Permeability development in the two rocks differs because grain boundaries control melt distribution to a greater extent in the gneiss. Muscovite‐dehydration melting may develop melt pathways at low melt fractions due to a larger volume of melt, in comparison with biotite‐dehydration melting, generated at the solidus. This may be a viable physical mechanism in which rapid melt segregation from a metapelitic source rock can occur. Alternatively, the results from the gneiss experiments suggest continual draining of biotite‐derived magma from the lower crust with melt migration paths controlled by structural anisotropies in the protolith.     

An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

Rate of strike-slip motion on the Amanos Fault (Karasu Valley, southern Turkey) constrained by K–Ar dating and geochemical analysis of Quaternary basalts

The left-lateral Amanos Fault follows a 200-km-long and up to 2-km-high escarpment that bounds the eastern margin of the Amanos mountain range and the western margin of the Karasu Valley in southern Turkey, just east of the northeastern corner of the Mediterranean Sea. Regional kinematic models have reached diverse conclusions as to the role of this fault in accommodating relative motion between either the African and Arabian, Turkish and African, or Turkish and Arabian plates. Local studies have tried to estimate its slip rate by K–Ar dating Quaternary basalts that erupted within the Amanos Mountains, flowed across it into the Karasu Valley, and have since become offset. However, these studies have yielded a wide range of results, ranging from 0.3 to 15 mm a⁻¹, which do not allow the overall role and significance of this fault in accommodating crustal deformation to be determined. We have used the Cassignol K–Ar method to date nine Quaternary basalt samples from the vicinity of the southern part of the Amanos Fault. These basalts exhibit a diverse chemistry, which we interpret as a consequence varying degrees of partial melting of their source combined with variable crustal contamination. This dating allows us to constrain the Quaternary slip rate on the Amanos fault to 1.0 to 1.6 mm a⁻¹. The dramatic discrepancies between past estimates of this slip rate are partly due to technical difficulties in K–Ar dating of young basalts by isotope dilution. In addition, previous studies at the key locality of Hacılar have unwittingly dated different, chemically distinct, flow units of different ages that are juxtaposed. This low slip rate indicates that, at present, the Amanos Fault takes up a small proportion of the relative motion between the African and Arabian plates, which is transferred southward to the Dead Sea Fault Zone. It also provides strong evidence against the long-standing view that its slip continues offshore to the southwest along a hypothetical left-lateral fault zone located south of Cyprus.     

Mantle structure and rifting processes in the Baikal–Mongolia region: geophysical data and evidence from xenoliths in volcanic rocks

The origin of the Baikal rift zone (BRZ) has been debated between the advocates of passive and active rifting since the 1970s. A re-assessment of the relevant geological and geophysical data from Russian and international literature questions the concept of broad asthenospheric upwelling beneath the rift zone that has been the cornerstone of many “active rifting” models. Results of a large number of early and recent studies favour the role of far-field forces in the opening and development of the BRZ. This study emphasises the data obtained through studies of peridotite and pyroxenite xenoliths brought to the surface by alkali basaltic magmas in southern Siberia and central Mongolia. These xenoliths are direct samples of the upper mantle in the vicinity of the BRZ. Of particular importance are suites of garnet-bearing xenoliths that have been used to construct P–T- composition lithospheric cross-sections in the region for the depth range of 35–80 km.Xenolith studies have shown fundamental differences in the composition and thermal regime between the lithospheric mantle beneath the ancient Siberian platform (sampled by kimberlites) and beneath younger mobile belts south of the platform. The uppermost mantle in southern Siberia and central Mongolia is much hotter at similar levels than the mantle in the Siberian craton and also has significantly higher contents of ‘basaltic’ major elements (Ca, Al, Na) and iron, higher Fe/Si and Fe/Mg. The combination of the moderately high geothermal gradient and the fertile compositions in the off-cratonic mantle appears to be a determining factor controlling differences in sub-Moho seismic velocities relative to the Siberian craton. Chemical and isotopic compositions of the off-cratonic xenoliths indicate small-scale and regional mantle heterogeneities attributed to various partial melting and enrichment events, consistent with long-term evolution in the lithospheric mantle. Age estimates of mantle events based on Os–Sr–Nd isotopic data can be correlated with major regional stages of crustal formation and may indicate long-term crust–mantle coupling. The ratios of ^143/144Nd in many LREE-depleted xenoliths are higher than those in MORB or OIB source regions and are not consistent with a recent origin from asthenospheric mantle.Mantle xenoliths nearest to the rift basins (30–50 km south of southern Lake Baikal) show no unequivocal evidence for strong heating, unusual stress and deformation, solid state flow, magmatic activity or partial melting that could be indicative of an asthenospheric intrusion right below the Moho. Comparisons between xenoliths from older and younger volcanic rocks east of Lake Baikal, together with observations on phase transformations and mineral zoning in individual xenoliths, have indicated recent heating in portions of the lithospheric mantle that may be related to localised magmatic activity or small-scale ascent of deep mantle material. Overall, the petrographic, P–T, chemical and isotopic constraints from mantle xenoliths appear to be consistent with recent geophysical studies, which found no evidence for a large-scale asthenospheric upwarp beneath the rift, and lend support to passive rifting mechanism for the BRZ.     

Pseudotachylytes in the southern border fault of the Cenozoic intracontinental Teletsk basin (Altai, Russia)

The Cenozoic intracontinental Teletsk basin in the Central Asian Altai Mountains is composed of a complexly structured northern and a more simple southern sub-basin. These sub-basins formed in two distinct kinematic stages when first the NNW-striking Teletsk- and then the NE-striking West-Sayan shear zones became reactivated in the Cenozoic under dominant NS-oriented horizontal compression. Although the entire Teletsk basin strikes roughly NS, the southern sub-basin is parallel to the NNW-trending, amphibolite facies Teletsk ductile shear zone, while the northern sub-basin is NS-striking and flanked by differently structured, greenschist facies basement. Basement reactivation closely controlled the southern sub-basin formation, but this is less clear for the northern sub-basin. Contrasts between northern and southern basement and the exclusive occurrence of pseudotachylytes along the margins of the southern basin are explored for their contribution to the formation of the Teletsk basin with two distinct sub-basins.In the ductile shear fabric of the basement flanking the southern sub-basin, concordantly interleaved pseudotachylytes and isolated breccia lenses reflect local brittle deformation along the ductile fabric. The genetic link between breccia lenses and pseudotachylyte occurrences was defined by microstructural investigation. It allows to explore their possible development in a dextral strike–slip zone. These rocks occur in a large fault-bounded segment of the basement. The geometry of the structures in the segment is comparable with a dextral strike–slip sidewall-ripout structure along the Teletsk shear zone. Seismic slip related to pseudotachylytes is attributed to the sudden stress release on the NNW-striking Teletsk shear zone, when the latter became unconstrained by reactivation of the NE-trending West-Sayan fault zone at its northern boundary. The boundary of the sidewall-ripout structure was reactivated as a large listric fault in a later stage. The northern sub-basins roughly develop along an NS strike and are assumed to reflect reactivation of the ductile shear zone underneath the variably structured greenschist facies basement outcropping along the flanks of the sub-basin.     

Ground deformation of the southern sector of the Aeolian islands volcanic arc from geodetic data

The deformation pattern and the dynamics of the southern sector of the Aeolian archipelago are investigated. A study on the ground deformation, measured over the last 20 years in the trilateration geodetic network between the islands of Vulcano and Lipari, has been conducted. Analysis of the relative displacements and the uniform strain tensor parameters, as well as the comparison between areal dilatation and the vertical variations deduced by precise levelling, allow distinguishing different phases associated both with the regional dynamics and the local volcanic context of the area. These phases, however, appear to be closely interrelated. The analysis of the deformation pattern allows to constrain the predominance of a roughly E–W trending extension and a N–S contraction at a regional scale. This regime is consistent with right-lateral movements along a NW–SE striking fault system.     

The presence, characteristics and earthquake implications of the St. Lawrence fault zone within and near Lake Ontario (Canada–USA)

Questions persist concerning the earthquake potential of the populous and industrial Lake Ontario (Canada–USA) area. Pertinent to those questions is whether the major fault zone that extends along the St. Lawrence River valley, herein named the St. Lawrence fault zone, continues upstream along the St. Lawrence River valley at least as far as Lake Ontario or terminates near Cornwall (Ontario, Canada)–Massena (NY, USA). New geological studies uncovered paleotectonic bedrock faults that are parallel to, and lie within, the projection of that northeast-oriented fault zone between Cornwall and northeastern Lake Ontario, suggesting that the fault zone continues into Lake Ontario. The aforementioned bedrock faults range from meters to tens of kilometers in length and display kinematically incompatible displacements, implying that the fault zone was periodically reactivated in the study area. Beneath Lake Ontario the Hamilton–Presqu'ile fault lines up with the St. Lawrence fault zone and projects to the southwest where it coincides with the Dundas Valley (Ontario, Canada). The Dundas Valley extends landward from beneath the western end of the lake and is marked by a vertical stratigraphic displacement across its width. The alignment of the Hamilton–Presqu'ile fault with the St. Lawrence fault zone strongly suggests that the latter crosses the entire length of Lake Ontario and continues along the Dundas Valley.The Rochester Basin, an east–northeast-trending linear trough in the southeastern corner of Lake Ontario, lies along the southern part of the St. Lawrence fault zone. Submarine dives in May 1997 revealed inclined layers of glaciolacustrine clay along two different scarps within the basin. The inclined layers strike parallel to the long dimension of the basin, and dip about 20° to the north–northwest suggesting that they are the result of rigid-body rotation consequent upon post-glacial faulting. Those post-glacial faults are growth faults as demonstrated by the consistently greater thickness, unit-by-unit, of unconsolidated sediments on the downthrown (northwest) side of the faults relative to their counterparts on the upthrown (southeast) side. Underneath the western part of Lake Ontario is a monoclinal warp that displaces the glacial and post-glacial sediments, and the underlying bedrock–sediment interface. Because of the post-glacial growth faults and the monoclinal warp the St. Lawrence fault zone is inferred to be tectonically active beneath Lake Ontario. Furthermore, within the lake it crosses at least five major faults and fault zones and coexists with other neotectonic structures. Those attributes, combined with the large earthquakes associated with the St. Lawrence fault zone well to the northeast of Lake Ontario, suggest that the seismic risk in the area surrounding and including Lake Ontario is likely much greater than previously believed.