Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.     

Limiting the soil degradational impacts of wildfire in pine and eucalyptus forests in Portugal: A comparison of alternative post-fire management practices

In newly burnt and unburnt pine and eucalyptus forest in Portugal, overland flow and soil losses were monitored to assess the impacts of the following post-fire treatments: application of different quantities of logging litter; rip-ploughing compared with minimum tillage prior to planting eucalyptus seedlings; and clearance of pine needles and vegetation. Eucalyptus logging litter reduced soil losses by up to 95 per cent. The impact of pine logging litter was equivocal, but removal of pine needles increased soil losses elevenfold. Implications for soil longevity, soil quality and land management strategy are discussed.     

H2O activity in concentrated NaCl solutions at high pressures and temperatures measured by the brucite-periclase equilibrium

 H₂O activities in concentrated NaCl solutions were measured in the ranges 600°–900° C and 2–15 kbar and at NaCl concentrations up to halite saturation by depression of the brucite (Mg(OH)₂) – periclase (MgO) dehydration equilibrium. Experiments were made in internally heated Ar pressure apparatus at 2 and 4.2 kbar and in 1.91-cm-diameter piston-cylinder apparatus with NaCl pressure medium at 4.2, 7, 10 and 15 kbar. Fluid compositions in equilibrium with brucite and periclase were reversed to closures of less than 2 mol% by measuring weight changes after drying of punctured Pt capsules. Brucite-periclase equilibrium in the binary system was redetermined using coarsely crystalline synthetic brucite and periclase to inhibit back-reaction in quenching. These data lead to a linear expression for the standard Gibbs free energy of the brucite dehydration reaction in the experimental temperature range: ΔG° (±120J)=73418–134.95T(K). Using this function as a baseline, the experimental dehydration points in the system MgO−H₂O−NaCl lead to a simple systematic relationship of high-temperature H₂O activity in NaCl solution. At low pressure and low fluid densities near 2 kbar the H₂O activity is closely approximated by its mole fraction. At pressures of 10 kbar and greater, with fluid densities approaching those of condensed H₂O, the H₂O activity becomes nearly equal to the square of its mole fraction. Isobaric halite saturation points terminating the univariant brucite-periclase curves were determined at each experimental pressure. The five temperature-composition points in the system NaCl−H₂O are in close agreement with the halite saturation curves (liquidus curves) given by existing data from differential thermal analysis to 6 kbar. Solubility of MgO in the vapor phase near halite saturation is much less than one mole percent and could not have influenced our determinations. Activity concentration relations in the experimental P-T range may be retrieved for the binary system H₂O-NaCl from our brucite-periclase data and from halite liquidus data with minor extrapolation. At two kbar, solutions closely approach an ideal gas mixture, whereas at 10 kbar and above the solutions closely approximate an ideal fused salt mixture, where the activities of H₂O and NaCl correspond to an ideal activity formulation. This profound pressure-induced change of state may be characterized by the activity (a) – concentration (X) expression: a _{H 2}O=X _{H 2}O/(1+αX _NaCl), and a _NaCl=(1+α)^(1+α)[X _NaCl/(1+αX _NaCl)]^(1+α). The parameter α is determined by regression of the brucite-periclase H₂O activity data: α=exp[A–B/ϱ_{H 2}O ]-CP/T, where A=4.226, B=2.9605, C=164.984, and P is in kbar, T is in Kelvins, and ϱ_{H 2}O is the density of H₂O at given P and T in g/cm³. These formulas reproduce both the H₂O activity data and the NaCl activity data with a standard deviation of ±0.010. The thermodynamic behavior of concentrated NaCl solutions at high temperature and pressure is thus much simpler than portrayed by extended Debye-Hückel theory. The low H₂O activity at high pressures in concentrated supercritical NaCl solutions (or hydrosaline melts) indicates that such solutions should be feasible as chemically active fluids capable of coexisting with solid rocks and silicate liquids (and a CO₂-rich vapor) in many processes of deep crustal and upper mantle metamorphism and metasomatism. Received: 1 September 1995 / Accepted: 24 March 1996     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

Radiometric characteristics of heavy mineral deposits along the west coast of South Africa

During the last decade, exploration and mining of modern-Tertiary heavy mineral beach and raised beach sands along the west coast of South Africa has developed into a major industry. High resolution radiometric techniques have demonstrated their use as a quantitative indicator of total heavy mineral concentration (THM) and also have the ability to discriminate between sediments derived from different provenance terrains. Results indicate that it is possible to calculate the total heavy mineral concentration from high resolution radiometric measurements, but the ability to quantify the concentration of individual mineral fractions, such as ilmenite, requires further refinement of the method. Radiometric characteristics of the light mineral fraction made it possible to distinguish between sediments from mixed marine-aeolian palaeoplacers and active present beach placers. The effects of in situ alteration of the ilmenite fraction in the palaeoplacers were also reflected by subtle contrasts in the radiometric character of the deposits. Furthermore, the uniform radiometric character of the sample population indicates a common provenance for the heavy mineral suite and supports previous results which indicate the metamorphites of the mid-Proterozoic Namaqua Province as the primary source terrain. The high resolution radiometric techniques discussed offer a broad spectrum of applications in sedimentology and have the potential to aid heavy mineral exploration. With further refinement it can possibly be used quantitatively in grade control during mining and beneficiation of heavy mineral deposits. Received: 19 July 1996 / Accepted: 7 January 1997     

Early Archean processes: evidence from the South African Kaapvaal craton and its greenstone belts (extended abstract)

Geologie en Mijnbouw -     

H2O activity in concentrated KCl and KCl-NaCl solutions at high temperatures and pressures measured by the brucite-periclase equilibrium

H₂O activities in supercritical fluids in the system KCl-H₂O-(MgO) were measured at pressures of 1, 2, 4, 7, 10 and 15  kbar by numerous reversals of vapor compositions in equilibrium with brucite and periclase. Measurements spanned the range 550–900 °C. A change of state of solute KCl occurs as pressures increase above 2 kbar, by which H₂O activity becomes very low and, at pressures of 4 kbar and above, nearly coincident with the square of the mole fraction (x _H2O). The effect undoubtedly results primarily from ionic dissociation as H₂O density (ρ_H2O) approaches 1 gm/cm³, and is more pronounced than in the NaCl-H₂O system at the same P-T-X conditions. Six values of solute KCl activity were yielded by terminal points of the isobaric brucite-periclase T-x _H2O curves where sylvite saturation occurs. The H₂O mole fraction of the isobaric invariant assemblage brucite-periclase-sylvite-fluid is near 0.52 at all pressures, and the corresponding temperatures span only 100 °C between 1 and 15 kbar. This remarkable convergence of the isobaric equilibrium curves reflects the great influence of pressure on lowering of both KCl and H₂O activities. The H₂O and KCl activities can be expressed by the formulas: a _H2O = γ_H2O[x _H2O+(1 + (1 + α)x _KCl)], and a _KCL = γ_KCl[(1 + α)x _KCl/(x _H2O +(1 + α)x _KCl)]^(1 + α), where α is a degree of dissociation parameter which increases from zero at the lowest pressures to near one at high pressures and the γ's are activity coefficients based on an empirical regular solution parameter W: ln γ_i = (1 − x_i)²W. Least squares fitting of our H₂O and KCl activity data evaluates the parameters: α = exp(4.166 −2.709/ρ_H2O) − 212.1P/T, and W = (−589.6 − 23.10P) /T, with ρ_H2O in gm/cm³, P in kbar and T in K. The standard deviation from the measured activities is only ± 0.014. The equations define isobaric liquidus curves, which are in perfect agreement with previous DTA liquidus measurements at 0.5–2 kbar, but which depart progressively from their extrapolation to higher pressures because of the pressure-induced dissociation effect. The great similarity of the NaCl-H₂O and KCl-H₂O systems suggests that H₂O activities in the ternary NaCl-KCl-H₂O system can be described with reasonable accuracy by assuming proportionality between the binary systems. This assumption was verified by a few reconnaissance measurements at 10 kbar of the brucite-periclase equilibrium with a Na/(Na + K) ratio of 0.5 and of the saturation temperature for Na/(Na + K) of 0.35 and 0.50. At that pressure the brucite-periclase curves reach a lowest x _H2O of 0.45 and a temperature of 587 °C before salt saturation occurs, values considerably lower than in either binary. This double-salt eutectic effect may have a significant application to natural polyionic hypersaline solutions in the deep crust and upper mantle in that higher solute concentrations and very low H₂O activities may be realized in complex solutions before salt saturation occurs. Concentrated salt solutions seem, from this standpoint, and also because of high mechanical mobility and alkali-exchanging potential, feasible as metasomatic fluids for a variety of deep-crust and upper mantle processes. Received: 9 August 1996 / Accepted: 15 November 1996     

Geochemical exploration for platinum-group elements in the Bushveld Complex, South Africa

Analyses of stream sediment and soil samples from the Bushveld Complex, South Africa have revealed enhanced precious metal concentrations, which can be related both to mining activities and the presence of hidden concentrations of platinum-group elements (PGEs) and gold. The economically important PGE deposits hosted by the Upper Critical Zone of the Rustenburg Layered Suite are revealed by a high PGE and Au content in the overlying soils. A second zone of elevated precious metal concentrations straddles the boundary between the Main and Upper Zones and has to date been traced for more than 100 km. This zone follows the igneous layering of the Rustenburg Layered Suite and is offset by the Brits Graben. It is therefore thought to be the reflection of a magmatic PGE-Au mineralisation. Received: 31 May 1996 / Accepted: 7 January 1997