Possible open-system (hydraulic) pingos in and around the Argyre impact region of Mars

We report the observation of possible (hydraulic) open-system pingos (OSPs ) at the mid latitudes (∼37°S) in and around the Argyre impact-basin. OSPs are perennial (water)–ice cored mounds; they originate and evolve in periglacial and pro-glacial landscapes on Earth where intra- or sub-permafrost water under hydraulic/artesian pressure uplifts localised sections of surface or near-surface permafrost that then freezes in-situ. We invoke three lines of evidence in support of our analogue-based interpretation: (1) similarities of shape, size and summit traits between terrestrial OSPs and the Martian mounds; (2) clustered distribution and the slope-side location of the mounds, consistent with terrestrial permafrost-environments where OSPs are found; and, (3) spatially-associated landforms putatively indicative of periglacial and glacial processes on Mars that characterise OSP landscapes on Earth. This article presents five OSP candidate-locations and nests these mound locations within a new geological map of the Argyre impact-basin and margins. It also presents three periglacial hypotheses about the possible origin of the water required to develop the mounds. Alternative (non-periglacial) formation-hypotheses also are considered; however, we show that their robustness is not equal to that of the periglacial ones.

Predictive modelling of mineral scaling, corrosion and the performance of solute geothermometers in a granitoid-hosted, enhanced geothermal system

Experience is lacking with mineral scaling and corrosion in enhanced geothermal systems (EGS) in which surface water is circulated through hydraulically stimulated crystalline rocks. As an aid in designing EGS projects we have conducted multicomponent reactive-transport simulations to predict the likely characteristics of scales and corrosion that may form when exploiting heat from granitoid reservoir rocks at ∼200 °C and 5 km depth. The specifications of an EGS project at Basel, Switzerland, are used to constrain the model. The main water–rock reactions in the reservoir during hydraulic stimulation and the subsequent doublet operation were identified in a separate paper (Alt-Epping et al., 2013b). Here we use the computed composition of the reservoir fluid to (1) predict mineral scaling in the injection and production wells, (2) evaluate methods of chemical geothermometry and (3) identify geochemical indicators of incipient corrosion. The envisaged heat extraction scheme ensures that even if the reservoir fluid is in equilibrium with quartz, cooling of the fluid will not induce saturation with respect to amorphous silica, thus eliminating the risk of silica scaling. However, the ascending fluid attains saturation with respect to crystalline aluminosilicates such as albite, microcline and chlorite, and possibly with respect to amorphous aluminosilicates. If no silica-bearing minerals precipitate upon ascent, reservoir temperatures can be predicted by classical formulations of silica geothermometry. In contrast, Na/K concentration ratios in the production fluid reflect steady-state conditions in the reservoir rather than albite–microcline equilibrium. Thus, even though igneous orthoclase is abundant in the reservoir and albite precipitates as a secondary phase, Na/K geothermometers fail to yield accurate temperatures. Anhydrite, which is present in fractures in the Basel reservoir, is predicted to dissolve during operation. This may lead to precipitation of pyrite and, at high exposure of anhydrite to the circulating fluid, of hematite scaling in the geothermal installation. In general, incipient corrosion of the casing can be detected at the production wellhead through an increase in H2(aq) and the enhanced precipitation of Fe-bearing aluminosilicates. The appearance of magnetite in scales indicates high corrosion rates.

A thermodynamic model for di-trioctahedral chlorite from experimental and natural data in the system MgO-FeO-Al2O3-SiO2-H2O. Applications to P-T sections and geothermometry

We present a new thermodynamic activity-composition model for di-trioctahedral chlorite in the system FeO–MgO–Al2O3–SiO2–H2O that is based on the Holland–Powell internally consistent thermodynamic data set. The model is formulated in terms of four linearly independent end-members, which are amesite, clinochlore, daphnite and sudoite. These account for the most important crystal-chemical substitutions in chlorite, the Fe–Mg, Tschermak and di-trioctahedral substitution. The ideal part of end-member activities is modeled with a mixing-on-site formalism, and non-ideality is described by a macroscopic symmetric (regular) formalism. The symmetric interaction parameters were calibrated using a set of 271 published chlorite analyses for which robust independent temperature estimates are available. In addition, adjustment of the standard state thermodynamic properties of sudoite was required to accurately reproduce experimental brackets involving sudoite. This new model was tested by calculating representative P–T sections for metasediments at low temperatures (<400 °C), in particular sudoite and chlorite bearing metapelites from Crete. Comparison between the calculated mineral assemblages and field data shows that the new model is able to predict the coexistence of chlorite and sudoite at low metamorphic temperatures. The predicted lower limit of the chloritoid stability field is also in better agreement with petrological observations. For practical applications to metamorphic and hydrothermal environments, two new semi-empirical chlorite geothermometers named Chl(1) and Chl(2) were calibrated based on the chlorite + quartz + water equilibrium (2 clinochlore + 3 sudoite = 4 amesite + 4 H2O + 7 quartz). The Chl(1) thermometer requires knowledge of the (Fe3+/ΣFe) ratio in chlorite and predicts correct temperatures for a range of redox conditions. The Chl(2) geothermometer which assumes that all iron in chlorite is ferrous has been applied to partially recrystallized detrital chlorite from the Zone houillère in the French Western Alps.

Komatiites constrain molybdenum isotope composition of the Earth's mantle

In order to estimate the Mo isotope composition and Mo abundance in the Bulk Silicate Earth (BSE), a total of thirty komatiite samples from five localities on three continents were analyzed using an isotope dilution double spike technique. Calculated Mo concentrations of the emplaced komatiite lavas range from 25±325±3 to 66±22 ng/g66±22 ng/g, and the inferred Mo concentrations in the deep mantle sources of the komatiites range between 17±417±4 and 30±12 ng/g30±12 ng/g, with an average value of 23±7 ng/g23±7 ng/g (2SE). This average value represents our best estimate for the Mo concentration in the BSE; it is identical, within the uncertainty, to published previous estimates of 39±16 ng/g39±16 ng/g, but is at least a factor of 2 more precise. The Mo isotope compositions of the komatiite mantle sources overlap within uncertainty and range from View the MathML sourceδMo98=−0.04±0.28 to 0.11±0.10‰0.11±0.10‰, with an average of 0.04±0.06‰0.04±0.06‰ (2SE). This value is analytically indistinguishable from published Mo isotope compositions of ordinary and enstatite chondrites and represents the best estimate for the Mo isotope composition of the BSE. The inferred δ98Mo for the BSE is therefore lighter than the suggested average of the upper continental crust (0.3 to 0.4‰). Thus, from the mass balance standpoint, a reservoir with lighter Mo isotope composition should exist in the Earth's mantle; this reservoir can potentially be found in subducted oceanic crust. The similarity of δ98Mo between chondritic meteorites and estimates for the BSE from this study indicates that during the last major equilibration between Earth's core and mantle, i.e., the one that occurred during the giant impact that produced the Moon, chemical and isotopic equilibrium of Mo between Fe metal of the core and the silicate mantle was largely achieved.

Geometry and kinematics of the Roisan-Cignana Shear Zone, and the orogenic evolution of the Dent Blanche Tectonic System (Western Alps)

The Dent Blanche Tectonic System (DBTS) is a composite thrust sheet derived from the previously thinned passive Adriatic continental margin. A kilometric high-strain zone, the Roisan-Cignana Shear Zone (RCSZ) defines the major tectonic boundary within the DBTS and separates it into two subunits, the Dent Blanche s.s. nappe to the northwest and the Mont Mary nappe to the southeast. Within this shear zone, tectonic slices of Mesozoic and pre-Alpine meta-sediments became amalgamated with continental basement rocks of the Adriatic margin. The occurrence of high pressure assemblages along the contact between these tectonic slices indicates that the amalgamation occurred prior to or during the subduction process, at an early stage of the Alpine orogenic cycle. Detailed mapping, petrographic and structural analysis show that the Roisan-Cignana Shear Zone results from several superimposed Alpine structural and metamorphic stages. Subduction of the continental fragments is recorded by blueschist-facies deformation, whereas the Alpine collision is reflected by a greenschist facies overprint associated with the development of large-scale open folds. The postnappe evolution comprises the development of low-angle brittle faults, followed by large-scale folding (Vanzone phase) and finally brittle extensional faults. The RCSZ shows that fragments of continental crust had been torn off the passive continental margin prior to continental collision, thus recording the entire history of the orogenic cycle. The role of preceding Permo-Triassic lithospheric thinning, Jurassic rifting, and ablative subduction processes in controlling the removal of crustal fragments from the reactivated passive continental margin is discussed. Results of this study constrain the temporal sequence of the tectono-metamorphic processes involved in the assembly of the DBTS, but they also show limits on the interpretation. In particular it remains difficult to judge to what extent precollisional rifting at the Adriatic continental margin preconditioned the efficiency of convergent processes, i.e. accretion, subduction, and orogenic exhumation.

Modelling the system behaviour of wet snow avalanches using an expert system approach for risk management on high alpine traffic roads

The presented approach describes a model for a rule-based expert system calculating the temporal variability of the release of wet snow avalanches, using the assumption of avalanche triggering without the loading of new snow. The knowledge base of the model is created by using investigations on the system behaviour of wet snow avalanches in the Italian Ortles Alps, and is represented by a fuzzy logic rule-base. Input parameters of the expert system are numerical and linguistic variables, measurable meteorological and topographical factors and observable characteristics of the snow cover. Output of the inference method is the quantified release disposition for wet snow avalanches. Combining topographical parameters and the spatial interpolation of the calculated release disposition a hazard index map is dynamically generated. Furthermore, the spatial and temporal variability of damage potential on roads exposed to wet snow avalanches can be quantified, expressed by the number of persons at risk. The application of the rule base to the available data in the study area generated plausible results. The study demonstrates the potential for the application of expert systems and fuzzy logic in the field of natural hazard monitoring and risk management.

Origin and history of ureilitic material in the solar system: The view from asteroid 2008 TC3 and the Almahata Sitta meteorite

Asteroid 2008 TC3 (approximately 4m diameter) was tracked and studied in space for approximately 19h before it impacted Earth's atmosphere, shattering at 44-36km altitude. The recovered samples (>680 individual rocks) comprise the meteorite Almahata Sitta (AhS). Approximately 50-70% of these are ureilites (ultramafic achondrites). The rest are chondrites, mainly enstatite, ordinary, and Rumuruti types. The goal of this work is to understand how fragments of so many different types of parent bodies became mixed in the same asteroid. Almahata Sitta has been classified as a polymict ureilite with an anomalously high component of foreign clasts. However, we calculate that the mass of fallen material was 0.1% of the pre-atmospheric mass of the asteroid. Based on published data for the reflectance spectrum of the asteroid and laboratory spectra of the samples, we infer that the lost material was mostly ureilitic. Therefore, 2008 TC3 probably contained only a few percent nonureilitic materials, similar to other polymict ureilites except less well consolidated. From available data for the AhS meteorite fragments, we conclude that 2008 TC3 samples essentially the same range of types of ureilitic and nonureilitic materials as other polymict ureilites. We therefore suggest that the immediate parent of 2008 TC3 was the immediate parent of all ureilitic material sampled on Earth. We trace critical stages in the evolution of that material through solar system history. Based on various types of new modeling and re-evaluation of published data, we propose the following scenario. (1) The ureilite parent body (UPB) accreted 0.5-0.6Ma after formation of calcium-aluminum-rich inclusions (CAI), beyond the ice line (outer asteroid belt). Differentiation began approximately 1Ma after CAI. (2) The UPB was catastrophically disrupted by a major impact approximately 5Ma after CAI, with selective subsets of the fragments reassembling into daughter bodies. (3) Either the UPB (before breakup), or one of its daughters (after breakup), migrated to the inner belt due to scattering by massive embryos. (4) One daughter (after forming in or migrating to the inner belt) became the parent of 2008 TC3. It developed a regolith, mostly 3.8Ga ago. Clasts of enstatite, ordinary, and Rumuruti-type chondrites were implanted by low-velocity collisions. (5) Recently, the daughter was disrupted. Fragments were injected or drifted into Earth-crossing orbits. 2008 TC3 comes from outer layers of regolith, other polymict ureilites from deeper regolith, and main group ureilites from the interior of this body. In contrast to other models that have been proposed, this model invokes a stochastic history to explain the unique diversity of foreign materials in 2008 TC3 and other polymict ureilites.

Ecosystem regimes and responses in a coupled ancient lake system from MIS 5b to present: the diatom record of lakes Ohrid and Prespa

We reconstruct the aquatic ecosystem interactions since the last interglacial period in the oldest, most diverse, hydrologically connected European lake system, by using palaeolimnological diatom and selected geochemistry data from Lake Ohrid “DEEP site” core and equivalent data from Lake Prespa core, Co1215. Driven by climate forcing, the lakes experienced two adaptive cycles during the last 92 ka: "interglacial and interstadial" and "glacial" cycle. The short-term ecosystems reorganizations, e.g. regime shifts within these cycles substantially differ between the lakes, as evident from the inferred amplitudes of variation. The deeper Lake Ohrid shifted between ultra oligo- and oligotrophic regimes in contrast to the much shallower Lake Prespa, which shifted from a deeper, (oligo-) mesotrophic to a shallower, eutrophic lake and vice versa. Due to the high level of ecosystem stability (e.g. trophic state, lake level), Lake Ohrid appears relatively resistant to external forcing, such as climate and environmental change. Recovering in a relatively short time from major climate change, Lake Prespa is a resilient ecosystem. At the DEEP site, the decoupling between the lakes' response to climate change is marked in the prolonged and gradual changes during the MIS 5/4 and 2/1 transitions. These response differences and the lakes' different physical and chemical properties may limit the influence of Lake Prespa on Lake Ohrid. Regime shifts of Lake Ohrid due to potential hydrological change in Lake Prespa are not evident in the data presented here. Moreover, a complete collapse of the ecosystems functionality and loss of their diatom communities did not happen in either lake for the period presented in the study.

Characterising the three-dimensional ozone distribution of a tidally locked Earth-like planet

We simulate the 3D ozone distribution of a tidally locked Earth-like exoplanet using the high-resolution, 3D chemistry climate model CESM1(WACCM) and study how the ozone layer of a tidally locked Earth (TLE) (ΩTLE = 1/365 days) differs from that of our present-day Earth (PDE) (ΩPDE = 1/1 day). The middle atmosphere reaches a steady state a symptotically within the first 80 days of the simulation. An upwelling, centred on the subsolar point, is present on the day side while a downwelling, centred on the antisolar point, is present on the night side. In the mesosphere, we find similar global ozone distributions for the TLE and the PDE, with decreased ozone on the day side and enhanced ozone on the night side. In the lower mesosphere, a jet stream transitions into a large-scale vortex around a low-pressure system, located at low latitudes of the TLE night side. In the middle stratosphere, the concentration of odd oxygen is approximately equal to that of the ozone [(Ox) ≈ (O3)]. At these altitudes, the lifetime of odd oxygen is ~16 h and the transport processes significantly contribute to the global distribution of stratospheric ozone. Compared to the PDE, where the strong Coriolis force acts as a mixing barrier between low and high latitudes, the transport processes of the TLE are governed by jet streams variable in the zonal and meridional directions. In the middle stratosphere of the TLE, we find high ozone values on the day side, due to the increased production of atomic oxygen on the day side, where it immediately recombines with molecular oxygen to form ozone. In contrast, the ozone is depleted on the night side, due to changes in the solar radiation distribution and the presence of a downwelling. As a result of the reduced Coriolis force, the tropical and extratropical air masses are well mixed and the global temperature distribution of the TLE stratosphere has smaller horizontal gradients than the PDE. Compared to the PDE, the total ozone column global mean is reduced by ~19.3 %. The day side and the night side total ozone column means are reduced by 23.21 and 15.52 %, respectively. Finally, we present the total ozone column (TOC) maps as viewed by a remote observer for four phases of the TLE during its revolution around the star. The mean TOC values of the four phases of the TLE vary by up to 23 %.

The isotope composition of selenium in chondrites constrains the depletion mechanism of volatile elements in solar system materials

Solar nebula processes led to a depletion of volatile elements in different chondrite groups when compared to the bulk chemical composition of the solar system deduced from the Sun's photosphere. For moderately-volatile elements, this depletion primarily correlates with the element condensation temperature and is possibly caused by incomplete condensation from a hot solar nebula, evaporative loss from the precursor dust, and/or inherited from the interstellar medium. Element concentrations and interelement ratios of volatile elements do not provide a clear picture about responsible mechanisms. Here, the abundance and stable isotope composition of the moderately- to highly-volatile element Se are investigated in carbonaceous, ordinary, and enstatite chondrites to constrain the mechanism responsible for the depletion of volatile elements in planetary bodies of the inner solar system and to define a δ(82/78)Se value for the bulk solar system. The δ(82/78)Se of the studied chondrite falls are identical within their measurement uncertainties with a mean of −0.20±0.26‰ (2 s.d., n=14n=14, relative to NIST SRM 3149) despite Se abundance depletions of up to a factor of 2.5 with respect to the CI group. The absence of resolvable Se isotope fractionation rules out a kinetic Rayleigh-type incomplete condensation of Se from the hot solar nebula or partial kinetic evaporative loss on the precursor material and/or the parent bodies. The Se depletion, if acquired during partial condensation or evaporative loss, therefore must have occurred under near equilibrium conditions to prevent measurable isotope fractionation. Alternatively, the depletion and cooling of the nebula could have occurred simultaneously due to the continuous removal of gas and fine particles by the solar wind accompanied by the quantitative condensation of elements from the pre-depleted gas. In this scenario the condensation of elements does not require equilibrium conditions to avoid isotope fractionation. The results further suggest that the processes causing the high variability of Se concentrations and depletions in ordinary and enstatite chondrites did not involve any measurable isotope fractionation. Different degrees of element depletions and isotope fractionations of the moderately-volatile elements Zn, S, and Se in ordinary and enstatite chondrites indicate that their volatility is controlled by the thermal stabilities of their host phases and not by the condensation temperature under canonical nebular conditions.