Complexity and non-linearity in earth surface processes - concepts, methods and applications

Complexity has long been recognized and is increasingly becoming mainstream in geomorphology. However, the relative novelty of various concepts and techniques associated to it means that ambiguity continues to surround complexity. In this commentary, we present and discuss a variety of recent contributions that have the potential to help clarify issues and advance the use of complexity in geomorphology.

Microscale mapping of alteration conditions and potential biosignatures in basaltic-ultramafic rocks on early Earth and beyond

Subseafloor environments preserved in Archean greenstone belts provide an analogue for investigating potential subsurface habitats on Mars. The c. 3.5-3.4 Ga pillow lava metabasalts of the mid-Archean Barberton greenstone belt, South Africa, have been argued to contain the earliest evidence for microbial subseafloor life. This includes candidate trace fossils in the form of titanite microtextures, and sulfur isotopic signatures of pyrite preserved in metabasaltic glass of the c. 3.472 Ga Hooggenoeg Formation. It has been contended that similar microtextures in altered martian basalts may represent potential extraterrestrial biosignatures of microbe-fluid-rock interaction. But despite numerous studies describing these putative early traces of life, a detailed metamorphic characterization of the microtextures and their host alteration conditions in the ancient pillow lava metabasites is lacking. Here, we present a new nondestructive technique with which to study the in situ metamorphic alteration conditions associated with potential biosignatures in mafic-ultramafic rocks of the Hooggenoeg Formation. Our approach combines quantitative microscale compositional mapping by electron microprobe with inverse thermodynamic modeling to derive low-temperature chlorite crystallization conditions. We found that the titanite microtextures formed under subgreenschist to greenschist facies conditions. Two chlorite temperature groups were identified in the maps surrounding the titanite microtextures and record peak metamorphic conditions at 315 ± 40°C (XFe3+(chlorite) = 25-34%) and lower-temperature chlorite veins/microdomains at T = 210 ± 40°C (lower XFe3+(chlorite) = 40-45%). These results provide the first metamorphic constraints in textural context on the Barberton titanite microtextures and thereby improve our understanding of the local preservation conditions of these potential biosignatures. We suggest that this approach may prove to be an important tool in future studies to assess the biogenicity of these earliest candidate traces of life on Earth. Furthermore, we propose that this mapping approach could also be used to investigate altered mafic-ultramafic extraterrestrial samples containing candidate biosignatures.

Lithium isotope constraints on crust–mantle interactions and surface processes on Mars

Lithium abundances and isotope compositions are reported for a suite of martian meteorites that span the range of petrological and geochemical types recognized to date for Mars. Samples include twenty-one bulk-rock enriched, intermediate and depleted shergottites, six nakhlites, two chassignites, the orthopyroxenite Allan Hills (ALH) 84001 and the polymict breccia Northwest Africa (NWA) 7034. Shergottites unaffected by terrestrial weathering exhibit a range in δ7Li from 2.1 to 6.2‰, similar to that reported for pristine terrestrial peridotites and unaltered mid-ocean ridge and ocean island basalts. Two chassignites have δ7Li values (4.0‰) intermediate to the shergottite range, and combined, these meteorites provide the most robust current constraints on δ7Li of the martian mantle. The polymict breccia NWA 7034 has the lowest δ7Li (−0.2‰) of all terrestrially unaltered martian meteorites measured to date and may represent an isotopically light surface end-member. The new data for NWA 7034 imply that martian crustal surface materials had both a lighter Li isotope composition and elevated Li abundance compared with their associated mantle. These findings are supported by Li data for olivine-phyric shergotitte NWA 1068, a black glass phase isolated from the Tissint meteorite fall, and some nakhlites, which all show evidence for assimilation of a low-δ7Li crustal component. The range in δ7Li for nakhlites (1.8 to 5.2‰), and co-variations with chlorine abundance, suggests crustal contamination by Cl-rich brines. The differences in Li isotope composition and abundance between the martian mantle and estimated crust are not as large as the fractionations observed for terrestrial continental crust and mantle, suggesting a difference in the styles of alteration and weathering between water-dominated processes on Earth versus possibly Cl–S-rich brines on Mars. Using high-MgO shergottites (>15 wt.% MgO) it is possible to estimate the δ7Li of Bulk Silicate Mars (BSM) to be 4.2 ± 0.9‰ (2σ). This value is at the higher end of estimates for the Bulk Silicate Earth (BSE; 3.5 ± 1.0‰, 2σ), but overlaps within uncertainty.

Transient Earth system responses to cumulative carbon dioxide emissions: linearities, uncertainties, and probabilities in an observation-constrained model ensemble

Information on the relationship between cumulative fossil CO2 emissions and multiple climate targets is essential to design emission mitigation and climate adaptation strategies. In this study, the transient response of a climate or environmental variable per trillion tonnes of CO2 emissions, termed TRE, is quantified for a set of impact-relevant climate variables and from a large set of multi-forcing scenarios extended to year 2300 towards stabilization. An  ∼ 1000-member ensemble of the Bern3D-LPJ carbon–climate model is applied and model outcomes are constrained by 26 physical and biogeochemical observational data sets in a Bayesian, Monte Carlo-type framework. Uncertainties in TRE estimates include both scenario uncertainty and model response uncertainty. Cumulative fossil emissions of 1000 Gt C result in a global mean surface air temperature change of 1.9 °C (68 % confidence interval (c.i.): 1.3 to 2.7 °C), a decrease in surface ocean pH of 0.19 (0.18 to 0.22), and a steric sea level rise of 20 cm (13 to 27 cm until 2300). Linearity between cumulative emissions and transient response is high for pH and reasonably high for surface air and sea surface temperatures, but less pronounced for changes in Atlantic meridional overturning, Southern Ocean and tropical surface water saturation with respect to biogenic structures of calcium carbonate, and carbon stocks in soils. The constrained model ensemble is also applied to determine the response to a pulse-like emission and in idealized CO2-only simulations. The transient climate response is constrained, primarily by long-term ocean heat observations, to 1.7 °C (68 % c.i.: 1.3 to 2.2 °C) and the equilibrium climate sensitivity to 2.9 °C (2.0 to 4.2 °C). This is consistent with results by CMIP5 models but inconsistent with recent studies that relied on short-term air temperature data affected by natural climate variability.

The tectonometamorphic evolution of the Sesia–Dent Blanche nappes (internal Western Alps): review and synthesis

This study reviews and synthesizes the present knowledge on the Sesia–Dent Blanche nappes, the highest tectonic elements in the Western Alps (Switzerland and Italy), which comprise pieces of pre-Alpine basement and Mesozoic cover. All of the available data are integrated in a crustal-scale kinematic model with the aim to reconstruct the Alpine tectono-metamorphic evolution of the Sesia–Dent Blanche nappes. Although major uncertainties remain in the pre-Alpine geometry, the basement and cover sequences of the Sesia–Dent Blanche nappes are seen as part of a thinned continental crust derived from the Adriatic margin. The earliest stages of the Alpine evolution are interpreted as recording late Cretaceous subduction of the Adria-derived Sesia–Dent Blanche nappes below the South-Alpine domain. During this subduction, several sheets of crustal material were stacked and separated by shear zones that rework remnants of their Mesozoic cover. The recently described Roisan-Cignana Shear Zone of the Dent Blanche Tectonic System represents such a shear zone, indicating that the Sesia–Dent Blanche nappes represent a stack of several individual nappes. During the subsequent subduction of the Piemonte–Liguria Ocean large-scale folding of the nappe stack (including the Roisan-Cignana Shear Zone) took place under greenschist facies conditions, which indicates partial exhumation of the Dent Blanche Tectonic System. The entrance of the Briançonnais micro-continent within the subduction zone led to a drastic change in the deformation pattern of the Alpine belt, with rapid exhumation of the eclogite-facies ophiolite bearing units and thrust propagation towards the foreland. Slab breakoff probably was responsible for allowing partial melting in the mantle and Oligocene intrusions into the most internal parts of the Sesia–Dent Blanche nappes. Finally, indentation of the Adriatic plate into the orogenic wedge resulted in the formation of the Vanzone back-fold, which marks the end of the pervasive ductile deformation within the Sesia–Dent Blanche nappes during the earliest Miocene.

Land system science and sustainable development of the earth system: A global land project perspective

Land systems are the result of human interactions with the natural environment. Understanding the drivers, state, trends and impacts of different land systems on social and natural processes helps to reveal how changes in the land system affect the functioning of the socio-ecological system as a whole and the tradeoff these changes may represent. The Global Land Project has led advances by synthesizing land systems research across different scales and providing concepts to further understand the feedbacks between social-and environmental systems, between urban and rural environments and between distant world regions. Land system science has moved from a focus on observation of change and understanding the drivers of these changes to a focus on using this understanding to design sustainable transformations through stakeholder engagement and through the concept of land governance. As land use can be seen as the largest geo-engineering project in which mankind has engaged, land system science can act as a platform for integration of insights from different disciplines and for translation of knowledge into action.

The behaviour of incompatible elements during hydrous melting of metasomatized peridotite at 4–6 GPa and 1000 °C–1200 °C

Trace element behavior during hydrous melting of a metasomatized garnet–peridotite was examined at pressures of 4–6 GPa and temperatures of 1000 °C–1200 °C, conditions appropriate for fluid penetrating the mantle wedge atop the subducting slab. Experiments were performed in a rocking multi-anvil apparatus using a diamond-trap setup. The compositions of the fluid and melt phases were measured using the cryogenic LA-ICP-MS technique. The water-saturated solidus of the K-lherzolite composition is located between 900 °C and 1000 °C at 4 GPa and between 1000 °C and 1100 °C at 5 and 6 GPa. The partition coefficients between fluid or melt and clinopyroxene reveal an asymmetric MREE trough with a minimum at Dy. The clinopyroxene in equilibrium with aqueous fluids is characterized by DUfluid–cpx > DThfluid–cpx while DUmelt–cpx tends to be similar to DThmelt–cpx. The partition coefficients between fluid or melt and garnet reveal very strong light to heavy REE fractionation, DLa/DLu from 95 (hydrous melt) to 1600 (aqueous fluid). The LILE are highly incompatible with partition coefficients > 50. The behavior of HFSE are decoupled, with DZr,Hf close to 1 while DNb,Ta > 10. Garnet is characterized by DUmelt/fluid–garnet < DThmelt/fluid–garnet. A comparison of our experimental partitioning results for trivalent cations as well as the results from the literature and the calculations carried out using the lattice strain model adapted to the presence of water in the bulk system indicates that H2O in the fluid or melt phase has a prominent effect on trace element partitioning. Garnet in mantle rocks in equilibrium with an aqueous fluid is characterized by significantly higher Do(3 +) for REE in the X site of the garnet compared with the partitioning values of the optimal cation in garnet in equilibrium with hydrous melts. Our data show for the first time that the change in the nature of the mobile phase (fluid vs. melt) does affect the affinities of trace elements into the garnet crystal at conditions below the second critical endpoint of the system. The same also applies for clinopyroxene, although this is less clear. Consequently, our new data allow for refinements in predictive modeling of element transfer from the slab to the mantle wedge and of possible compositions of metasomatized mantle that sources OIB magmatism.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

MRD's approach: Ideally suited to support Future Earth research agenda in mountains

Both Future Earth and Mountain Research and Development (MRD) aim to support production and dissemination of knowledge for sustainable development. As shown in Future Earth’s Strategic Research Agenda 2014, the global research community has begun to acknowledge its societal role and the need for a new type of research in which scientists link disciplines and coproduce transformation knowledge with stakeholders. Future Earth has defined three research themes that conceptualize the issues to be dealt with at the same time as the way in which this should be done. In many ways, MRD’s policy has made the journal a forerunner of Future Earth’s stipulated “step-change in research”. Indeed, MRD’s section policies aim to support similar contents and ways of producing these forms of knowledge. MRD publishes “systems knowledge” in its MountainResearch section, “target knowledge” in its MountainAgenda section, and “transformation knowledge” in its MountainDevelopment section. Each of these sections has dedicated review criteria to assess and enhance the quality of the knowledge presented in the papers. In this poster, we provide examples from each of the three sections of what the knowledge types look like, how they are assessed, and how they contribute to the three Future Earth themes.

Complex geomorphologic assemblage of terrains in association with the banded terrain in Hellas basin, Mars

Hellas basin acts as a major sink for the southern highlands of Mars and is likely to have recorded several episodes of sedimentation and erosion. The north-western part of the basin displays a potentially unique Amazonian landscape domain in the deepest part of Hellas, called “banded terrain”, which is a deposit characterized by an alternation of narrow band shapes and inter-bands displaying a sinuous and relatively smooth surface texture suggesting a viscous flow origin. Here we use high-resolution (HiRISE and CTX) images to assess the geomorphological interaction of the banded terrain with the surrounding geomorphologic domains in the NW interior of Hellas to gain a better understanding of the geological evolution of the region as a whole. Our analysis reveals that the banded terrain is associated with six geomorphologic domains: a central plateau named Alpheus Colles, plain deposits (P1 and P2), reticulate (RT1 and RT2) and honeycomb terrains. Based on the analysis of the geomorphology of these domains and their cross-cutting relationships, we show that no widespread deposition post-dates the formation of the banded terrain, which implies that this domain is the youngest and latest deposit of the interior of Hellas. Therefore, the level of geologic activity in the NW Hellas during the Amazonian appears to have been relatively low and restricted to modification of the landscape through mechanical weathering, aeolian and periglacial processes. Thermophysical data and cross-cutting relationships support hypotheses of modification of the honeycomb terrain via vertical rise of diapirs such as ice diapirism, and the formation of the plain deposits through deposition and remobilization of an ice-rich mantle deposit. Finally, the observed gradual transition between honeycomb and banded terrain suggests that the banded terrain may have covered a larger area of the NW interior of Hellas in the past than previously thought. This has implications on the understanding of the evolution of the deepest part of Hellas.

Eclogitic diamonds from variable crustal protoliths in the northeastern Siberian craton: trace elements and coupled δ13C – δ18O signatures in diamonds and garnet inclusions

Diamonds of eclogitic assemblages are dominant in the placer diamond deposits of the northeastern Siberian platform. In this study we present new trace elements and stable isotopes (δ13C and δ18O) data for alluvial diamonds and their garnet inclusions from this locality. Cr-rich garnets of peridotitic affinity in the studied diamonds have a narrow range of δ18O values from 5.7‰ to 6.2‰, which is largely overlapping with the accepted mantle range. This narrow range suggests that the garnet inclusions showing different REE patterns and little variations in oxygen isotopes may have formed by different processes involving fluid/melts that, however, were in oxygen isotopic equilibrium with the mantle. The trace element composition of the eclogitic garnet inclusions supports a crustal origin for at least the high-Ca garnets, which show flat HREE patterns and in some cases a positive Eu-anomaly. High-Ca eclogitic garnets generally show heavier oxygen isotope compositions (δ18O 6.5–9.6‰) than what is observed in low-Ca garnets (δ18O 5.7–7.4‰). The variability in oxygen isotopes and trace elements is suggested to be inherited from contrasting crustal protoliths. The relationship between the high δ18O values of inclusions and the low δ13C values of the host diamonds implies that the high-Ca garnet inclusions were derived from intensely hydrated (e.g., δ18O > 7‰) and typically oxidised basaltic rock close to the seawater interface, and that the carbon for diamonds was closely associated with this protolith.

The supergene thorium and rare-earth element deposit at Morro do Ferro, Poços de Caldas, Minas Gerais, Brazil

The thorium and rare-earth element (Th-REE) deposit at Morro do Ferro formed under supergene lateritic weathering conditions. The ore body consists of shallow NW-SE elongated argillaceous lenses that extend from the top of the hill downwards along its south-eastern slope. The deposit is capped by a network of magnetite layers which protected the underlying highly weathered, argillaceous host rock from excessive erosion. The surrounding country rocks comprise a sequence of subvolcanic phonolite intrusions that have been strongly altered by hydrothermal and supergene processes. From petrological, mineralogical and geochemical studies, and mass balance calculations, it is inferred that the highly weathered host rock was originally carbonatitic in composition, initially enriched in Th and REEs compared to the surrounding silicate rocks. The intrusion of the carbonatite caused fenitic alteration in the surrounding phonolites, consisting of early potassic alteration followed by a vein-type Th-REE mineralization with associated fluorite, carbonate, pyrite and zircon. Subsequent weathering has completely decomposed the carbonatite forming a residual supergene enrichment of Th and REEs. Initial weathering of the carbonatite has created a chemical environment that might have been conductive to carbonate and phosphate complexing of the REEs in groundwaters. This may have appreciably restricted the dissolution of primary REE phases. Strongly oxidic weathering has resulted in a fractionation between Ce and the other light rare earth elements (LREEs). Ce3+ is oxidized to Ce4+ and retained together with Th by secondary mineral formation (cerianite, thorianite), and by adsorption on poorly crystalline iron- and aluminium-hydroxides. In contrast, the trivalent LREEs are retained to a lesser degree and are thus more available for secondary mineral formation (Nd-lanthanite) and adsorption at greater depths down the weathering column. Seasonally controlled fluctuations of recharge waters into the weathering column may help to explain the observed repetition of Th-Ce enriched zones underlain by trivalent LREE enriched zones.

Radiochemical Determination of Rare Earth Elements in Proton-Irradiated Lead–Bismuth Eutectic

Various types of proton-irradiated lead–bismuth eutectic (LBE) samples from the MEGAPIE prototype spallation target were analyzed concerning their content of 148Gd, 173Lu, and 146Pm by use of α- and γ-spectrometry. A radiochemical separation procedure was developed to isolate the lanthanide fraction and to prepare thin samples for α-ray measurement. The results prove a substantial depletion of these three elements in bulk samples, whereas accumulation on the LBE/steel-interfaces was observed. The amount of material accumulated on surfaces was roughly estimated by relating the values measured on the sample surfaces to the total surface of the inner target walls. The amount of 148Gd, 173Lu, and 146Pm was then quantified by summing up the contributions from every sample type. The results show a reasonable agreement with theoretical predictions. The obtained results are of utmost importance for the evaluation of the performance of high-power spallation targets, especially concerning the residual nuclide production, the physicochemical behavior of the produced radionuclides during operation, and in terms of an intermediate or final disposal.

Earth system commitments due to delayed mitigation

As long as global CO₂ emissions continue to increase annually, long-term committed Earth system changes grow much faster than current observations. A novel metric linking this future growth to policy decisions today is the mitigation delay sensitivity (MDS), but MDS estimates for Earth system variables other than peak temperature (ΔT max) are missing. Using an Earth System Model of Intermediate Complexity, we show that the current emission increase rate causes a ΔT max increase roughly 3–7.5 times as fast as observed warming, and a millenial steric sea level rise (SSLR) 7–25 times as fast as observed SSLR, depending on the achievable rate of emission reductions after the peak of emissions. These ranges are only slightly affected by the uncertainty range in equilibrium climate sensitivity, which is included in the above values. The extent of ocean acidification at the end of the century is also strongly dependent on the starting time and rate of emission reductions. The preservable surface ocean area with sufficient aragonite supersaturation for coral reef growth is diminished globally at an MDS of roughly 25%–80% per decade. A near-complete loss of this area becomes unavoidable if mitigation is delayed for a few years to decades. Also with respect to aragonite, 12%–18% of the Southern Ocean surface become undersaturated per decade, if emission reductions are delayed beyond 2015–2040. We conclude that the consequences of delaying global emission reductions are much better captured if the MDS of relevant Earth system variables is communicated in addition to current trends and total projected future changes.