Silica depleted melting of pelites. Petrogenetic grid and application to the Susqueda Aureole, Spain

Metapelitic rocks in the low pressure contact metamorphic aureole around the Susqueda igneous complex, Spain show a number of features that make them an ideal testing ground for the modelling of silica‐undersaturated melting. Rocks in the aureole experienced localized depletion in silica by the segregation of quartz veins during a pre‐anatectic, regional cordierite‐andalusite grade metamorphic event. These rocks were then intruded by gabbroic to dioritic rocks of the Susqueda igneous complex that formed a migmatitic contact metamorphic aureole in the country rocks. This migmatisation event caused quartz‐saturated hornfels and restite formation in rocks that had experienced no quartz vein segregation in the previous regional metamorphic event, but silica‐undersaturated melting in those rocks that were previously depleted in silica. Silica‐undersaturated melting is investigated using a new petrogenetic P–T projection and equilibrium pseudosections calculated in the KFMASH and NCKFMASH systems, respectively. The grid considers quartz absent equilibria and a range of phases that form typically in silica‐undersaturated bulk compositions, for example corundum. It is shown that the quartz‐rich precursors in the Susqueda contact aureole produced about 10% melt during contact metamorphism. However, most of this melt was extracted leaving behind rocks with restitic bulk compositions and minor leucosome segregation. It is suggested that the melt mixed with the host igneous rocks causing an apparent magmatic zoning from diorite in the centre of the complex to tonalite at the margins. In contrast, the quartz‐poor precursors (from which the quartz veins segregated) melted in the silica‐undersaturated field producing a range of assemblages including peritectic corundum and spinel. Melting of the silica‐undersaturated rocks produced only negligible melt and no subsequent melt loss.     

Pillow talk

Three distinct types of pillows and pillow lava sequences with different modes of origin have been recognized in the extrusive sequences comprising the upper parts of ophiolite complexes that represent the mafic portion of the floor of an Early Cretaceous back-arc basin in southern Chile. One type of pillow formed by non-explosive submarine effusion. A second type formed by magmatic intrusion into pre-existing aquagene tuff formed by explosive eruption. The third type of pillow occurs within dikes, forming pillowed dikes, possibly as a result of vapor streaming within a cooling dike. Where studied in southern Chile, aquagene tuffs and intrusive pillows decrease and water-lain pillows increase in relative abundance from north to south. This variation corresponds with a north-to-south decrease in both the relative volume of extrusives to extensional dikes and the range and volume of differentiated rocks, suggesting a southward increase in rate of extension relative to rate of magma supply within the spreading ridges at which the ophiolites formed. In the northern part of the original basin where the rate of extension was small relative to the rate of magma supply, magma remained in magma chambers longer, resulting in a greater range and volume of differentiated rocks. The larger volume of more differentiated, cooler and more viscous magmas, in conjunction with the likelihood of more violent eruption of volatile-rich differentiates, may have been responsible for the large volume of aquagene tuff in the northern part of the original basin. These observations in southern Chile suggest that ophiolites which contain a great abundance of aquagene tuffs and intrusive pillow lavas formed in tectonic environments in which the rate of extension was small relative to the rate of magma supply (island arcs, embryonic marginal basins). Ophiolites with predominantly water-lain pillowed and massive lavas formed in tectonic environments in which the rate of extension was large relative to the rate of magma supply (mid-ocean ridges, mature back-arc basins). Thus geologic field data may supplement geochemical data as a tool in distinguishing the original igneous-tectonic environments in which ophiolites originate.     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

Using amphibole phenocrysts to track vapor transfer during magma crystallization and transport: An example from Mount St. Helens,Washington

In order to evaluate and further constrain models for volatile movement and vapor enrichment of magma stored at shallow levels, amphibole phenocrysts from 2004–2005 Mount St. Helens dacite were analyzed for major and selected trace elements (Li, Cu, Zn, Mn, and REE) and Li isotopes. Several recent studies have examined fluid-mobile trace element abundances in phencryst phases and melt inclusions as a means of tracking volatile movement within subvolcanic magmatic systems, and high Li contents in plagioclase phenocrysts from 1980 and 2004 Mount St. Helens dacites have been interpreted as evidence that shallow magma was fluxed by a Li-bearing vapor phase prior to eruption.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

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.     

Horizontal and vertical components of earthquake ground motions at liquefiable sites

Field observations on ground motions from recent earthquakes imply that current knowledge is limited with regard to relating vertical and horizontal motions at liquefiable sites. This paper describes a study with the purpose of clarifying this emerging issue to some extent. A series of numerical analyses is carried out on a liquefiable soil deposit with a verified, fully coupled, nonlinear procedure. It is shown that the transformation of vertical motions in the deposit differs considerably from the transformation of horizontal motions. Both the amplitude and frequency content of the horizontal motions are strongly dependent on the shaking level or the associated nonlinear soil behavior. The transfer function for vertical motions is however likely to be independent of the intensity of input motions; no reduction in the amplitude occurs even in the case of strong shaking. The results are shown to be in consistence with the laboratory observations on shaking table tests and recent field observations that less nonlinearity exists for vertical motions. It is also shown that the possibility exists for using information on spectral ratios between the horizontal and vertical surface motions to quickly identify in situ soil behavior and liquefaction that are not readily covered by conventional field or laboratory experimentation procedures.     

Crystal features of supergene gold at Hannan South,Western Australia

Deep weathering of basalt-hosted sulphidic gold mineralisation has resulted in remobilisation of the gold and the formation of a supergene deposit. The deposit occurs beneath an ephemeral lake system and is at least partially saturated by highly saline groundwater. A general downward movement of an iron redox front associated with the weathering has controlled the distribution of the gold and its morphology. It is unusual in that the high fineness gold crystals formed are coarse and well-preserved. Small octahedra, up to 50 μm, are the most abundant form of gold crystals but large, up to 3 mm, octahedral plates are common. Crystal morphology variations described here suggest that the gold is initially precipitated as euhedral octahedra and plates. Subsequent dissolution and recrystallisation associated with fluctuations in the chemical environment above the redox front has led to the development of the more common irregular dendritic or wire gold forms observed in other supergene deposits.     

Characterization of groundwater chemistry under the influence of lithologic and anthropogenic factors along a climatic gradient in Upper Cauvery basin, South India

Hydrogeological and climatic effect on chemical behavior of groundwater along a climatic gradient is studied along a river basin. ‘Semi-arid’ (500–800 mm of mean annual rainfall), ‘sub-humid’ (800–1,200 mm/year) and ‘humid’ (1,200–1,500 mm/year) are the climatic zones chosen along the granito-gneissic plains of Kabini basin in South India for the present analysis. Data on groundwater chemistry is initially checked for its quality using NICB ratio (<±5 %), EC versus TZ+ (~0.85 correlation), EC versus TDS and EC versus TH analysis. Groundwater in the three climatic zones is ‘hard’ to ‘very hard’ in terms of Ca–Mg hardness. Polluted wells are identified (>40 % of pollution) and eliminated for the characterization. Piper’s diagram with mean concentrations indicates the evolution of CaNaHCO₃ (semi-arid) from CaHCO₃ (humid zone) along the climatic gradient. Carbonates dominate other anions and strong acids exceeded weak acids in the region. Mule Hole SEW, an experimental watershed in sub-humid zone, is characterized initially using hydrogeochemistry and is observed to be a replica of entire sub-humid zone (with 25 wells). Extension of the studies for the entire basin (120 wells) showed a chemical gradient along the climatic gradient with sub-humid zone bridging semi-arid and humid zones. Ca/Na molar ratio varies by more than 100 times from semi-arid to humid zones. Semi-arid zone is more silicaceous than sub-humid while humid zone is more carbonaceous (Ca/Cl ~14). Along the climatic gradient, groundwater is undersaturated (humid), saturated (sub-humid) and slightly supersaturated (semi-arid) with calcite and dolomite. Concentration–depth profiles are in support of the geological stratification i.e., ~18 m of saprolite and ~25 m of fracture rock with parent gneiss beneath. All the wells are classified into four groups based on groundwater fluctuations and further into ‘deep’ and ‘shallow’ based on the depth to groundwater. Higher the fluctuations, larger is its impact on groundwater chemistry. Actual seasonal patterns are identified using ‘recharge–discharge’ concept based on rainfall intensity instead of traditional monsoon–non-monsoon concept. Non-pumped wells have low Na/Cl and Ca/Cl ratios in recharge period than in discharge period (Dilution). Few other wells, which are subjected to pumping, still exhibit dilution chemistry though water level fluctuations are high due to annual recharge. Other wells which do not receive sufficient rainfall and are constantly pumped showed high concentrations in recharge period rather than in discharge period (Anti-dilution). In summary, recharge–discharge concept demarcates the pumped wells from natural deep wells thus, characterizing the basin.