Compositional zoning of fluid inclusions in the Archaean Junction gold deposit,Western Australia: A process of fluid ‐ wall‐rock interaction?

The Junction gold deposit, in Western Australia, is an orogenic gold deposit hosted by a differentiated, iron‐rich, tholeiitic dolerite sill. Petrographic, microthermometric and laser Raman microprobe analyses of fluid inclusions from the Junction deposit indicate that three different vein systems formed at three distinct periods of geological time, and host four fluid‐inclusion populations with a wide range of compositions in the H₂O–CO₂–CH₄–NaCl ± CaCl₂ system. Pre‐shearing, pre‐gold, molybdenite‐bearing quartz veins host fluid inclusions that are characterised by relatively consistent phase ratios comprising H₂O–CO₂–CH₄ ± halite. Microthermometry suggests that these veins precipitated when a highly saline, >340°C fluid mixed with a less saline ≥150°C fluid. The syn‐gold mineralisation event is hosted within the Junction shear zone and is associated with extensive quartz‐calcite ± albite ± chlorite ± pyrrhotite veining. Fluid‐inclusion analyses indicate that gold deposition occurred during the unmixing of a 400°C, moderately saline, H₂O–CO₂ ± CH₄ fluid at pressures between 70 MPa and 440 MPa. Post‐gold quartz‐calcite‐biotite‐pyrrhotite veins occupy normal fault sets that slightly offset the Junction shear zone. Fluid inclusions in these veins are predominantly vapour rich, with CO₂?CH₄. Homogenisation temperatures indicate that the post‐gold quartz veins precipitated from a 310 ± 30°C fluid. Finally, late secondary fluid inclusions show that a <200°C, highly saline, H₂O–CaCl₂–NaCl–bearing fluid percolated along microfractures late in the deposit's history, but did not form any notable vein type. Raman spectroscopy supports the microthermometric data and reveals that CH₄–bearing fluid inclusions occur in syn‐gold quartz grains found almost exclusively at the vein margin, whereas CO₂–bearing fluid inclusions occur in quartz grains that are found toward the centre of the veins. The zonation of CO₂:CH₄ ratios, with respect to the location of fluid inclusions within the syn‐gold quartz veins, suggest that the CH₄ did not travel as part of the auriferous fluid. Fluid unmixing and post‐entrapment alteration of the syn‐gold fluid inclusions are known to have occurred, but cannot adequately account for the relatively ordered zonation of CO₂:CH₄ ratios. Instead, the late introduction of a CH₄–rich fluid into the Junction shear zone appears more likely. Alternatively, the process of CO₂ reduction to CH₄ is a viable and plausible explanation that fits the available data. The CH₄–bearing fluid inclusions occur almost exclusively at the margin of the syn‐gold quartz veins within the zone of high‐grade gold mineralisation because this is where all the criteria needed to reduce CO₂ to CH₄ were satisfied in the Junction deposit.     

An iterative pressure‐stabilized fractional step algorithm in saturated soil dynamics

This paper presents an iterative, incremental pressure‐stabilized fractional step algorithm for coupled hydro‐mechanical problems with mixed formulations of the displacement–pressure ( u –p) model in saturated soil dynamics that allows the use of finite elements with equal low order of interpolation approximation of u and p. In comparison with the original fractional step algorithm, the distinct features of the proposed algorithm lie in its enhanced stability owing to the introduction of both an iteration procedure and a finite increment calculus (FIC) process into the algorithm. The introduction of the iterative procedure makes the velocity term satisfy the momentum conservation equation in an implicit sense and allows much larger time step sizes to be used than those limited in existing explicit and semi‐implicit versions of the algorithm. The introduction of the FIC process removes the dependence of the stability of the proposed algorithm on the time step size, as a result it allows to using the incremental version of the algorithm and evades the minimum time step size requirement presented in the existing versions of the fractional step algorithm that restricts the application of the algorithm to saturated soil dynamics problems with high frequencies. Numerical experiments demonstrate the effectiveness and improved performance of the proposed iterative pressure‐stabilized fractional step algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

Primary carbonate/CO2 inclusions in sapphirine-bearing granulites from central Sri Lanka

High-density CO₂-rich fluid inclusions from a sapphirine-bearing granulite (Hakurutale, Sri Lanka) have been studied by microthermometry, Raman spectrometry and SEM analysis. Based on textural evidence, two groups of inclusions can be identified: primary, negative crystal shaped inclusions (group I) and pseudo-secondary inclusions, which experienced a local, limited post-trapping modification (group II). Both groups contain magnesite as a daughter mineral, occurring in a relatively constant fluid/solid inclusion volume ratio (vol_solid =0.15 total volume). CO₂ densities for group I and II differ only slightly. Both groups contain a fluid, which was initially trapped at peak metamorphic conditions as a homogeneous (CO₂+MgCO₃) mixture. Thermodynamic calculations suggest that such a fluid (CO₂+15 vol% MgCO₃) is stable under granulite facies conditions. After trapping, magnesite separated upon cooling, while the remaining CO₂ density suffered minor re-adjustments. A model isochore based on the integration of the magnesite molar volume in the CO₂ fluid passes about 1.5–2 kbar below peak metamorphic conditions. This remaining discrepancy can be explained by the possible role of a small quantity of additional water.     

Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.     

The effects of vigorous mixing in a convective model of zonal flow on the ice giants

Previous studies have used models of three-dimensional (3D) Boussinesq convection in a rotating spherical shell to explain the zonal flows on the gas giants, Jupiter and Saturn. In this paper we demonstrate that this approach can also generate flow patterns similar to those observed on the ice giants, Uranus and Neptune. The equatorial jets of Uranus and Neptune are often assumed to result from baroclinic cloud layer processes and have been simulated with shallow layer models. Here we show that vigorous, 3D convection in a spherical shell can produce the retrograde (westward) equatorial flows that occur on the ice giants as well as the prograde (eastward) equatorial flows of the gas giants. In our models, the direction of the equatorial jet depends on the ratio of buoyancy to Coriolis forces in the system. In cases where Coriolis forces dominate buoyancy, cylindrical Reynolds stresses drive prograde equatorial jets. However, as buoyancy forces approach and exceed Coriolis forces, the cylindrical nature of the flow is lost and 3D mixing homogenizes the fluid's angular momentum; the equatorial jet reverses direction, while strong prograde jets form in the polar regions. Although the results suggest that conditions involving strong atmospheric mixing are responsible for generating the zonal flows on the ice giants, our present models require roughly 100 and 10 times the internal heat fluxes observed on Uranus and Neptune, respectively.     

Statistical and geostatistical investigations into the effects of the Gabcikovo hydropower plant on the groundwater resources of northwest Hungary / Analyses statistiques et géostatistiques des effets de la centrale hydroélectrique de Gabcikovo sur les ressources en eaux souterraines du Nord-Ouest de la Hongrie

Abstract

Abstract The construction of the Gabcikovo hydropower plant and the diversion of the Danube River over 25 km into an artificial channel in 1992 influenced the groundwater regime of the region considerably. Statistical and geostatistical methods are used to quantify changes of different groundwater characteristics on the Hungarian side of the river based on observations in the time period 1960–2000. External drift kriging was used to interpolate groundwater levels and the other related variables. While mean groundwater levels did not change appreciably, there are significant changes in the variability. Standard deviations of the groundwater levels and the amplitude of the annual cycle decreased near the old river bed of the Danube. The water-level fluctuations of the Danube have a decreased influence on the groundwater dynamics. Interrelationships between water levels in wells have also changed.     

Efficient methods of discharge measurements in rivers and streams based on the probability concept

Regularities exist in fluid flows and can be represented by a set of constants. These constants are functions of the parameter of a probability distribution that exhibits resilience and stability under various flow conditions. Together, these regularities form a network and interact with each other, such that if one is known then the others can be determined from it. The regularities and their network explain the various fluid‐flow phenomena and can be used in analysis of rivers and streams. For example, they can be used as the basis to develop simple and efficient methods of discharge measurements as presented herein, which only require velocity sampling at a single point on a water surface or a few points on a single vertical. Because of their simplicity and the short time requirement, these methods can be easily automated for collecting discharge data in unsteady, high flows that are badly needed for real‐time flow forecasting and design of flood control structures, and for advancing the fundamental, scientific knowledge in hydrology. Copyright © 2005 John Wiley & Sons, Ltd.     

Thermal history of the Hodgkinson Province and Laura Basin,Queensland: Multiple cooling episodes identified from apatite fission track analysis and vitrinite reflectance data

Apatite fission track analysis and vitrinite reflectance data from outcrop and well samples in the Hodgkinson Province and Laura Basin reveal regional Cretaceous cooling. Apatite fission track analysis appears to define two discrete cooling episodes, in the mid‐Cretaceous (110–100 Ma) and Late Cretaceous (80–70 Ma), although in most samples data allow only definition of a single episode. Rocks now at outcrop cooled from Cretaceous palaeotemperatures generally between 50 and 130°C in the south of the region, and from >100°C in the north. Some samples from the Hodgkinson Province also show evidence for an Early Jurassic cooling episode, characterised by maximum palaeotemperatures varying from at least 95°C (from apatite fission track analysis) to ～200–220°C (from vitrinite reflectance), with cooling beginning at around 200 Ma. Apatite fission track analysis data do not reveal the earlier event in the Laura Basin, but on the basis of vitrinite reflectance data from Permian? units this event is also inferred to have affected the pre‐Jurassic basin units in this region. The regional extent of the Cretaceous cooling episode in the Hodgkinson Province suggests that the elevated palaeotemperatures in this region were most likely due to greater depth of burial, with subsequent cooling due to kilometre‐scale denudation. For a palaeogeothermal gradient of 30°C/km and a palaeosurface temperature of 25°C the total degree of Cretaceous cooling experienced by the samples corresponds to removal of between ～0.8 and >3.0 km of Triassic and younger section removed by denudation, beginning some time between ca 110 and 80 Ma. Higher palaeogradients would require correspondingly lower amounts of removed section. The geology of the Laura Basin suggests that an explanation of the observed Cretaceous palaeotemperatures in this region in terms of deeper burial may be untenable. Heating due to hot fluid flow may be a more realistic mechanism for producing the observed Cretaceous palaeothermal effects in the Laura Basin.     

Ion winds in Saturn's southern auroral/polar region

We present profiles of the line-of-sight (l.o.s.) ionospheric wind velocities in the southern auroral/polar region of Saturn. Our velocities are derived from the measurement of Doppler shifting of the H₃⁺ν₂Q(1,0⁻) line at 3.953 microns. The data for this study were obtained using the facility high-resolution spectrometer CSHELL on the NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii, during the night of February 6, 2003 (UT). The l.o.s. velocity profiles finally derived are consistent with an extended region of the upper atmosphere sub-corotating with the planet: the ion velocities in the inertial reference are only 1/3 of those expected for full planetary corotation. We discuss the results in the light of recent proposals for the kronian magnetosphere, and suggest that, in this region, Saturn's ion winds may be under solar wind control.     

Regular and chaotic motion in softened gravitational systems

The stability of the dynamical trajectories of softened spherical gravitational systems is examined, both in the case of the full N -body problem and that of trajectories moving in the gravitational field of non-interacting background particles. In the latter case, for   N 10 000  , some trajectories, even if unstable, had exceedingly long diffusion times, which correlated with the characteristic e-folding time-scale of the instability. For trajectories of   N ≈100 000  systems this time-scale could be arbitrarily large – and thus appear to correspond to regular orbits. For centrally concentrated systems, low angular momentum trajectories were found to be systematically more unstable. This phenomenon is analogous to the well-known case of trajectories in generic centrally concentrated non-spherical smooth systems, where eccentric trajectories are found to be chaotic. The exponentiation times also correlate with the conservation of the angular momenta along the trajectories. For N up to a few hundred, the instability time-scales of N -body systems and their variation with particle number are similar to those of the most chaotic trajectories in inhomogeneous non-interacting systems. For larger N (up to a few thousand) the values of the these time-scales were found to saturate, increasing significantly more slowly with N . We attribute this to collective effects in the fully self-gravitating problem, which are apparent in the time variations of the time-dependent Liapunov exponents. The results presented here go some way towards resolving the long-standing apparent paradoxes concerning the local instability of trajectories. This now appears to be a manifestation of mechanisms driving evolution in gravitational systems and their interactions – and may thus be a useful diagnostic of such processes.