A molal-based model for strong acid chemistry at low temperatures (<200 to 298 K)

Geochemical processes occurring in cold environments on Earth, Mars, and Europa have elicited considerable interest in the application of geochemical models to subzero temperatures. Few existing geochemical models explicitly include acid chemistry and those that do are largely restricted to temperatures ≥0°C or rely on the mole-fraction scale rather than the more common molal scale. This paper describes (1) use of the Clegg mole-fraction acid models to develop a molal-based model for hydrochloric, nitric, and sulfuric acids at low temperatures; (2) incorporation of acid chemistry and nitrate minerals into the FREZCHEM model; (3) validation and limitations of the derived acid model; and (4) simulation of hypothetical acidic brines for Europa.The Clegg mole-fraction acid models were used to estimate activities of water and mean ionic activity coefficients that serve as the database for estimating molal Pitzer-equation parameters for HCl (188 to 298 K), HNO₃ (228 to 298 K), and H₂SO₄ (208 to 298 K). Model eutectics for HNO₃ and H₂SO₄ agree with experimental measurements to within ± 0.2°C. In agreement with previous work, the experimental freezing point depression (fpd) data for pure HCl at subzero temperatures were judged to be flawed and unreliable. Three alternatives are discussed for handling HCl chemistry at subzero temperatures. In addition to defining the solubility of solid-phase acids, this work also adds three new nitrate minerals and six new acid salts to the FREZCHEM model and refines equilibria among water ice, liquid water, and water vapor over the temperature range from 180 to 298 K. The final system is parameterized for Na-K-Mg-Ca-H-Cl-SO₄-NO₃-OH-HCO₃-CO₃-CO₂-H₂O.Simulations of hypothetical MgSO₄-H₂SO₄-H₂O and Na₂SO₄-MgSO₄-H₂SO₄-H₂O brines for Europa demonstrate how freezing can convert a predominantly salt solution into a predominantly acid solution at subzero temperatures. This result has consequences for the effects of salinity, acidity, and temperature as limiting factors for potential life on Europa. Strong acidity would limit life-forms to highly acidophilic organisms.     

Genetic significance of fluid inclusions in the CSA Cu–Pb–Zn deposit, Cobar, Australia

CSA mine exploits a ‘Cobar-type’ Cu–Pb–Zn±Au±Ag deposit within a cleaved and metamorphosed portion of the Cobar Supergroup, central New South Wales. The deposit comprises systems of ‘lenses’ that encompass veins, disseminations and semi-massive to massive Cu–Pb–Zn ores. The systems and contained lenses truncate bedding, are approximately coplanar with regional cleavage and similarly oriented shear zones and plunge parallel to the elongation lineation. Systems have extreme vertical continuity (>1000 m), short strike length (400 m) and narrow width (100 m), exhibit vertical and lateral ore-type variation and have alteration haloes. Models of ore formation include classical hydrothermalism, structurally controlled remobilisation and polymodal concepts; syntectonic emplacement now holds sway.Fluid inclusions were examined from quartz±sulphide veins adjacent to now-extracted ore, from coexisting quartz–sulphide within ore, and from vughs in barren quartz veins. Lack of early primary inclusions precluded direct determination of fluids associated with D₂–D₃ ore and vein emplacement. Similarly, decrepitation (by near-isobaric heating) of the two oldest secondary populations precluded direct determination of fluid phases immediately following D₂–D₃ ore and vein emplacement. Post-decrepitation outflow (late D₃ to early post-D₃) is recorded by monophase CH₄ inclusions. Entrained outflow of deeply circulated meteoric fluid modified the CH₄ system; modification is recorded by H₂O+CH₄ and H₂O+(trace CH₄) secondary populations and by an H₂O+(trace CH₄) primary population. The contractional tectonics (D₂–D₃) of ore emplacement was superseded by relaxational tectonics (D_4P) that facilitated meteoric water penetration and return flow.Under D₂ prograde metamorphism, entrapment temperatures (Tt) and pressures (Pt) for pre-decrepitation secondary inclusions are estimated as Tt300–330 °C and Pt1.5–2 kbar≈Plith (the lithostatic pressure). Decrepitation accompanied peak metamorphism (T350–380 °C) in mid- to late-D₃, while in late-D₃ to early post-D₃, essentially monophase CH₄ secondary inclusions were entrapped at Tt350 °C and Pt=1.5–2 kbar≈Plith. Subsequently, abundant CH₄ and entrained meteoric water were entrapped as H₂O+CH₄ secondaries under slowly decreasing temperature (Tt330–350 °C) and constant pressure (Pt1.5–2 kbar). Finally, with increasingly dominant meteoric outflow, H₂O+(trace CH₄) populations record decreasing temperatures (Tt>300 to <350 down to 275–300 °C) at pressures of Phydrostatic<Pt (1 kbar) <Plith (1.5 kbar).The populations of inclusions provide insight into fluid types, flow regimes and P–T conditions during parts of the deposit's evolution. They indirectly support the role of basin-derived CH₄ fluids in ore formation, but provide no insight into a basement-sourced ore-forming fluid. They fully support post-ore involvement of meteoric water. The poorly constrained entrapment history is believed to span 10 Ma from 395 to 385 Ma.     

Effects of pressure on aqueous chemical equilibria at subzero temperatures with applications to Europa

Pressure plays a critical role in controlling aqueous geochemical processes in deep oceans and deep ice. The putative ocean of Europa could have pressures of 1200 bars or higher on the seafloor, a pressure not dissimilar to the deepest ocean basin on Earth (the Mariana Trench at 1100 bars of pressure). At such high pressures, chemical thermodynamic relations need to explicitly consider pressure. A number of papers have addressed the role of pressure on equilibrium constants, activity coefficients, and the activity of water. None of these models deal, however, with processes at subzero temperatures, which may be important in cold environments on Earth and other planetary bodies. The objectives of this work were to (1) incorporate a pressure dependence into an existing geochemical model parameterized for subzero temperatures (FREZCHEM), (2) validate the model, and (3) simulate pressure-dependent processes on Europa. As part of objective 1, we examined two models for quantifying the volumetric properties of liquid water at subzero temperatures: one model is based on the measured properties of supercooled water, and the other model is based on the properties of liquid water in equilibrium with ice.The relative effect of pressure on solution properties falls in the order: equilibrium constants(K) > activity coefficients (γ) > activity of water (a_w). The errors (%) in our model associated with these properties, however, fall in the order: γ > K > a_w. The transposition between K and γ is due to a more accurate model for estimating K than for estimating γ. Only activity coefficients are likely to be significantly in error. However, even in this case, the errors are likely to be only in the range of 2 to 5% up to 1000 bars of pressure. Evidence based on the pressure/temperature melting of ice and salt solution densities argue in favor of the equilibrium water model, which depends on extrapolations, for characterizing the properties of liquid water in electrolyte solutions at subzero temperatures, rather than the supercooled water model. Model-derived estimates of mixed salt solution densities and chemical equilibria as a function of pressure are in reasonably good agreement with experimental measurements.To demonstrate the usefulness of this low-temperature, high-pressure model, we examined two hypothetical cases for Europa. Case 1 dealt with the ice cover of Europa, where we asked the question: How far above the putative ocean in the ice layer could we expect to find thermodynamically stable brine pockets that could serve as habitats for life? For a hypothetical nonconvecting 20 km icy shell, this potential life zone only extends 2.8 km into the icy shell before the eutectic is reached. For the case of a nonconvecting icy shell, the cold surface of Europa precludes stable aqueous phases (habitats for life) anywhere near the surface. Case 2 compared chemical equilibria at 1 bar (based on previous work) with a more realistic 1460 bars of pressure at the base of a 100 km Europan ocean. A pressure of 1460 bars, compared to 1 bar, caused a 12 K decrease in the temperature at which ice first formed and a 11 K increase in the temperature at which MgSO₄·12H₂O first formed. Remarkably, there was only a 1.2 K decrease in the eutectic temperatures between 1 and 1460 bars of pressure. Chemical systems and their response to pressure depend, ultimately, on the volumetric properties of individual constituents, which makes every system response highly individualistic.     

Energetic Charged Particles Above Thunderclouds

The French government has committed to launch the satellite TARANIS to study transient coupling processes between the Earth’s atmosphere and near-Earth space. The prime objective of TARANIS is to detect energetic charged particles and hard radiation emanating from thunderclouds. The British Nobel prize winner C.T.R. Wilson predicted lightning discharges from the top of thunderclouds into space almost a century ago. However, new experiments have only recently confirmed energetic discharge processes which transfer energy from the top of thunderclouds into the upper atmosphere and near-Earth space; they are now denoted as transient luminous events, terrestrial gamma-ray flashes and relativistic electron beams. This meeting report builds on the current state of scientific knowledge on the physics of plasmas in the laboratory and naturally occurring plasmas in the Earth’s atmosphere to propose areas of future research. The report specifically reflects presentations delivered by the members of a novel Franco-British collaboration during a meeting at the French Embassy in London held in November 2011. The scientific subjects of the report tackle ionization processes leading to electrical discharge processes, observations of transient luminous events, electromagnetic emissions, energetic charged particles and their impact on the Earth’s atmosphere. The importance of future research in this area for science and society, and towards spacecraft protection, is emphasized.     

Solar modulation in surface atmospheric electricity

The solar wind modulates the flux of galactic cosmic rays impinging on Earth inversely with solar activity. Cosmic ray ionisation is the major source of air's electrical conductivity over the oceans and well above the continents. Differential solar modulation of the cosmic ray energy spectrum modifies the cosmic ray ionisation at different latitudes, varying the total atmospheric columnar conductance. This redistributes current flow in the global atmospheric electrical circuit, including the local vertical current density and the related surface potential gradient. Surface vertical current density and potential gradient measurements made independently at Lerwick Observatory, Shetland, from 1978 to 1985 are compared with modelled changes in cosmic ray ionisation arising from solar activity changes. Both the lower troposphere atmospheric electricity quantities are significantly increased at cosmic ray maximum (solar minimum), with a proportional change greater than that of the cosmic ray change.     

Geology and tectonics of Neoproterozoic salt diapirs and salt sheets in the eastern Willouran Ranges,South Australia

Allochthonous salt structures and associated primary and secondary minibasins are exposed in Neoproterozoic strata of the eastern Willouran Ranges, South Australia. Detailed geologic mapping using high‐quality airborne hyperspectral remote‐sensing data and satellite imagery, combined with a qualitative structural restoration, are used to elucidate the evolution of this complex, long‐lived (>250 Myr) salt system. Field observations and interpretations at a resolution unobtainable from seismic or well data provide a means to test published models of allochthonous salt emplacement and associated salt‐sediment interaction derived from subsurface data in the northern Gulf of Mexico. Salt diapirs and sheets are represented by megabreccias of nonevaporite lithologies that were originally interbedded with evaporites that have been dissolved and/or altered. Passive diapirism began shortly after deposition of the Callanna Group layered evaporite sequence. A primary basin containing an expulsion‐rollover structure and megaflap is flanked by two vertical diapirs. Salt flowed laterally from the diapirs to form a complex, multi‐level canopy, now partly welded, containing an encapsulated minibasin and capped by suprasalt basins. Salt and minibasin geometries were modified during the Late Cambrian–Ordovician Delamerian Orogeny (ca. 500 Ma). Small‐scale structures such as subsalt shear zones, fractured or mixed ‘rubble zones’ and thrust imbricates are absent beneath allochthonous salt and welds in the eastern Willouran Ranges. Instead, either undeformed strata or halokinetic drape folds that include preserved diapir roof strata are found directly below the transition from steep diapirs to salt sheets. Allochthonous salt first broke through the diapir roofs and then flowed laterally, resulting in variable preservation of the subsalt drape folds. Lateral salt emplacement was presumably on roof‐edge thrusts or, because of the shallow depositional environment, via open‐toed advance or extrusive advance, but without associated subsalt deformation.     

Land cover classification by an artificial neural network with ancillary information

Abstract

Remote sensing is an important source of land cover data required by many GIS users. Land cover data are typically derived from remotely–sensed data through the application of a conventional statistical classification. Such classification techniques are not, however, always appropriate, particularly as they may make untenable assumptions about the data and their output is hard, comprising only the code of the most likely class of membership. Whilst some deviation from the assumptions may be tolerated and a fuzzy output may be derived, making more information on class membership properties available, alternative classification procedures are sometimes required. Artificial neural networks are an attractive alternative to the statistical classifiers and here one is used to derive a fuzzy classification output from a remotely–sensed data set that may be post–processed with ancillary data available in a GIS to increase the accuracy with which land cover may be mapped. With the aid ancillary information on soil type and prior knowledge of class occurrence the accuracy of an artificial neural network classification was increased by 29–93 to 77–37 per cent. An artificial neural network can therefore be used generate a fuzzy classification output that may be used with other data sets in a GIS, which may not have been available to the producer of the classification, to increase the accuracy with which land cover may be classified.