Hydrothermal Quartz Vein Formation, Revealed by Coupled SEM-CL Imaging and Fluid Inclusion Microthermometry: Shuteen Complex, South Gobi, Mongolia

Abstract. Scanning electron microscopy-cathodoluminescence (SEM-CL) imaging of vein quartz in the Cu-mineralised, Shuteen Complex (South Gobi, Mongolia) has revealed a complex history of crystal growth, dissolution and microfracture healing, associated with several hydrothermal events that could not be detected using other observational techniques (e.g. transmitted/reflected light microscopy, back-scattered electron imaging, or secondary electron imaging).
The quartz initially grew as CL-bright/grey crystals in a 345±30C liquid reservoir, as inferred by the analysis of primary liquid fluid inclusions (average Th of 343C; 6.6∼7.7 wt% NaCl_eq). Quartz precipitation occurred at the edge of the crystals as reservoir fluids cooled to 260±25C, as indicated by micron-scale CL-dark/CL-bright quartz growth bands containing abundant fluid inclusions (with an average T_h values of 261C). Pressure fluctuations were the likely cause of dissolution, as SEM-CL imaging reveals the quartz have corroded or rounded crystal edges, and precipitation of later quartz into open space. SEM-CL imaging shows the quartz contains healed microfractures that trapped low salinity fluids (3.9 wt% NaC1_eq) with Th values of 173±15C.
SEM-CL imaging provides a means of deciphering the thermal and chemical evolution of the fossil Shuteen hydrothermal system, and the nature of hydrothermal quartz vein-forming processes, by facilitating the correlation of distinct fluid inclusion populations and their relative chronology, with specific hydrothermal events.     

Evolution of Hydrothermal System at the Dizon Porphyry Cu-Au Deposit, Zambales, Philippines

Abstract. Evolution of hydrothermal system from initial porphyry Cu mineralization to overlapping epithermal system at the Dizon porphyry Cu‐Au deposit in western central Luzon, Zambales, Philippines, is documented in terms of mineral paragen‐esis, fluid inclusion petrography and microthermometry, and sulfur isotope systematics. The paragenetic stages throughout the deposit are summarized as follows; 1) stockwork amethystic quartz veinlets associated with chalcopyrite, bornite, magnetite and Au enveloped by chlorite alteration overprinting biotite alteration, 2) stockwork quartz veinlets with chalcopyrite and pyrite associated with Au and chalcopyrite and pyrite stringers in sericite alteration, 3) stringer quartz veinlets associated with molybdenite in sericite alteration, and 4) WNW‐trending quartz veins associated with sphalerite and galena at deeper part, while enargite and stibnite at shallower levels associated with advanced argillic alteration. Chalcopyrite and bornite associated with magnetite in quartz veinlet stockwork (stage 1) have precipitated initially as intermediate solid solution (iss) and bornite solid solution (bnss), respectively. Fluid inclusions in the stockwork veinlet quartz consist of gas‐rich inclusions and polyphase inclusions. Halite in polyphase inclusions dissolves at temperatures ranging from 360d?C to >500d?C but liquid (brine) and gas (vapor) do not homogenize at <500d?C. The maximum pressure and minimum temperature during the deposition of iss and bnss with stockwork quartz veinlets are estimated to be 460 bars and 500d?C. Fluid inclusions in veinlet stockwork quartz enveloped in sericite alteration (stage 2) consist mainly of gas‐rich inclusions and polyphase inclusions. In addition to the possible presence of saturated NaCl crystals at the time of entrapment of fluid inclusions that exhibit the liquid‐vapor homogenization temperatures lower than the halite dissolution temperatures in some samples, wide range of temperatures of halite dissolution and liquid‐vapor homogenization of polyphase inclusions from 230d?C to >500d?C and from 270d?C to >500d?C, respectively, suggests heterogeneous entrapment of gaseous vapor and hypersaline brine. The minimum pressure and temperature are estimated to be about 25 bars and 245d?C. Fluid inclusions in veinlet quartz associated with molybdenite (stage 3) are dominated by gas‐rich inclusions accompanied with minor liquid‐rich inclusions that homogenize at temperatures between 350d?C and 490d?C. Fluid inclusions in vuggy veinlet quartz associated with stibnite (stage 4) consist mainly of gas‐rich inclusions with subordinate polyphase inclusions that do not homogenize below 500d?C. Fluid inclusions in veinlet quartz associated with galena and sphalerite (stage 4) are composed of liquid‐rich two‐phase inclusions, and they homogenize into liquid phase at temperatures ranging widely from 190d?C to 300d?C (suggesting boiling) and the salinity ranges from 1.0 wt% to 3.4 wt% NaCl equivalent. A pressure of about 15 bars is estimated for the dilute aqueous solution of 190d?C from which veinlet quartz associated with galena and sphalerite precipitated. In addition to a change in temperature‐pressure regime from lithostatic pressure during the deposition of iss and bnss with stockwork quartz veinlets to hydrostatic pressure during fracture‐controlled quartz veinlet associated with galena and sphalerite, a decrease in pressure is supposed to have occurred due to unroofing or removal of the overlying piles during the temperature decrease in the evolution of hydrothermal system. The majority of the sulfur isotopic composition of sulfides ranges from ±0 % to +5 %. Sulfur originated from an iso‐topically uniform and homogeneous source, and the mineralization occurred in a single hydrothermal system.     

Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.     

Influence of the Neotethys rifting on the development of the Dampier Sub-basin (North West Shelf of Australia), highlighted by subsidence modelling

During the Late Palaeozoic and the Mesozoic, the development and evolution of the North West Shelf of Australia have been mostly driven by rifting phases associated with the break-up of Gondwana. These extensional episodes, which culminated in the opening of the Neotethys Ocean during the Permo-Carboniferous and a series of abyssal plains during the Jurassic-Cretaceous, are characterised by different stress regimes and modes of extension, and therefore had distinctive effects on the margin, and particularly on the Northern Carnarvon Basin.Interpretation of 3D and 2D seismic data enables a structural and stratigraphic analysis of the Late Palaeozoic sediments deposited in the proximal part of the Dampier Sub-basin (Mermaid Nose). Based on their seismic characters, stratigraphic relationship, internal patterns, lateral continuity, and architecture, these units are associated here with the Pennsylvanian?–Early Sakmarian glaciogenic Lyons Group and the Sakmarian–Artinskian Callytharra Formation. The former were deposited in a half-graben whose development is associated with the onset of the Neotethys rifting, and the latter is characterised by restricted deposition, inversion of prograding patterns, and uplift.The integration of seismo-stratigraphic characterisation of the Late Palaeozoic sequences and Mesozoic data from one exploration well (Roebuck-1) enables the construction of subsidence curves for the Mermaid Nose and the interpretation of its geohistory.The tectonic subsidence curves show a striking Permo-Carboniferous rifting phase related to the Neotethys rifting and a discrete Late Jurassic–Early Cretaceous event coeval with the opening and the spreading of the Argo Abyssal Plain.This result points out the predominance of the effects of the Permo-Carboniferous Neotethys episode, whereas the extension related to the Argo Abyssal Plain rifting that occurred later and closer to the studied area, had only limited effects on the subsidence of the proximal Dampier Sub-basin. Therefore, it supports a tectonic model with two distinct modes of extension for the Late Palaeozoic (widespread) and the Mesozoic (localised) rifting phases.     

Some considerations concerning seismic geomorphology and paleoseismology

Seismic geomorphology studies landforms which developed in connection with earthquakes. Among them, two different end members may be distinguished: 1) seismo-tectonic landforms, including surface faults and fractures, land uplift and subsidence at different scales, surface bulges, elongate ridges, and any other permanent ground deformations directly related to tectonic stress, and 2) seismo-gravitational landforms, such as landslides, deep-seated gravitational slope deformations, sinkholes, and fissures due to sediment compaction or liquefaction and sand blows, connected with both seismic shaking and gravitational stress.A clear-cut distinction between the two categories of landforms is not always easy to make (and in many instances not really useful), while there are, in many cases, ground effects that might be (and should be) considered as simultaneous combinations of seismo-tectonic and seismo-gravitational processes. This applies especially to surface fracturing and faulting which could be the combined result of tectonic stress, stress produced by seismic shaking, and gravitational stress.The objective of this paper is to review selected case histories mainly from Italy and the Mediterranean region, in order to show the importance of a comprehensive study of earthquake-generated landforms for understanding the seismicity level of the area under investigation. We argue that in earthquake prone areas, seismic landforms often constitute typical patterns (seismic landscapes) whose recognition, mapping and paleoseismological analysis may help in the evaluation of seismic hazards.     

Anomalous lithospheric structure of Northern Juan de Fuca plate — a consequence of oceanic rift propagation?

Anomalous crustal and upper mantle structure of northern Juan de Fuca plate is revealed from wide-angle seismic and gravity modelling. A 2-D velocity model is produced for refraction line II of the 1980 Vancouver Island Seismic Project (VISP80). The refraction data were recorded on three ocean bottom seismometers (OBSs) deployed at the ends and middle of a 110 km line oriented parallel to the North American continental margin. The velocity model is constructed via ray tracing and conforms to first-arrival amplitude observations and travel time picks of direct, converted and reflected phases. Between sub-sediment depths of 3 to 11 km, depths normally associated with the lower crust and upper oceanic mantle, the final model shows that compressional-wave velocities decrease significantly from southeast to northwest along the profile. At sub-sediment depths of 11 km at the northwestern end of the profile, P-wave velocities are as low as 7.2 km/s. A complementary 2-D gravity model using the geometry of the velocity model and velocity–density relationships characteristic of oceanic crust is produced. The high densities required to match the gravity field indicate the presence of peridotites containing 25–30% serpentine by volume, rather than excess gabbroic crust, within the deep low velocity zone. Anomalous travel time delays and unusual reflection characteristics observed from proximal seismic refraction and reflection experiments suggest a broader zone of partially serpentinized peridotites coincident with the trace of a pseudofault. We propose that partial serpentinization of the upper mantle is a consequence of slow spreading at the tip of a propagating rift.     

Directional dependence of P- and S-wave propagation and polarization in foliated rocks from the Kola superdeep well: Evidence from laboratory measurements and calculations based on TOF neutron diffraction

We have measured P- and S-wave velocities on two amphibolite and two gneiss samples from the Kola superdeep borehole as a function of pressure (up to 600 MPa) and temperature (up to 600 °C). The velocity measurements include compressional (Vp) and shear wave velocities (Vs1, Vs2) propagating in three orthogonal directions which were in general not parallel to inherent rock symmetry axes or planes. The measurements are accompanied by 3D-velocities calculations based on lattice preferred orientation (LPO) obtained by TOF (Time Of Flight) neutron diffraction analysis which allows the investigation of bulk volumes up to several cubic centimetres due to the high penetration depth of neutrons. The LPO-based numerical velocity calculations give important information on the different contribution of the various rock-forming minerals to bulk elastic anisotropy and on the relations of seismic anisotropy, shear wave splitting, and shear wave polarization to the structural reference frame (foliation and lineation). Comparison with measured velocities obtained for the three propagation directions that were not in accordance with the structural frame of the rocks (foliation and lineation) demonstrate that for shear waves propagating through anisotropic rocks the vibration directions are as important as the propagation directions. The study demonstrates that proper measurement of shear wave splitting by means of two orthogonal polarized sending and receiving shear wave transducers is only possible when their propagation and polarization directions are parallel and normal to foliation and lineation, respectively.     

Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

Determination of groundwater availability in shallow arid region aquifers utilizing GIS technology: a case study in Hada Al-Sham, Western Saudi Arabia

A simple, physically based method is developed in this paper to assist in the allocation of areas with high groundwater potential and for the determination of maximum allowed pumping rate to ensure proper groundwater management. This method utilizes the aquifer physical properties as well as GIS technology to accomplish this purpose. The design of this method was considered to be applicable in areas with little data, such as in most arid regions. This technique was applied to a catchment in an arid environment where qualitative as well as quantitative analyses of the results were undertaken. Locations of available groundwater and rates of maximum allowable pumping were compared with observations and experiments in the field and a good agreement was found. It was concluded that the best groundwater location was in the alluvial area, which represents only 16% of the total aquifer, which is a typical case in arid region catchments. The rate of maximum pumping was estimated to be 65 m³/h. However, to benefit 55% of the area, the maximum pumping rate should only be 40 m³/h with an average rate throughout the area (55%) of about 24 m³/h.This revised version was published online in December 2004 with corrections to the category.