Freshwater geographies: Experimenting with knowing and doing geography differently

The fair and effective governance of freshwater is an increasingly prominent issue in New Zealand. Emerging from a complex of cultural, economic and biophysical narratives, freshwater geographies are multiple, varied and increasingly acknowledged as worthy of interdisciplinary scrutiny. In this commentary, we reflect on a series of generative spaces that we – as group of postgraduate geographers (plus supporting staff) – created to engage with the multiplicity of freshwater meanings both within and beyond the academy. Through this evolving epistemic‐political project, we significantly reframed our own understandings about what freshwater ‘is’ and how it ought to be governed. By pursuing a deeper understanding of how the world gets made, we expand our ability to know and make it differently.     

Relationship among titanium,rare earth elements,U–Pb ages and deformation microstructures in zircon: Implications for Ti-in-zircon thermometry

A zircon grain in an orthopyroxene–garnet–phlogopite–zircon–rutile-bearing xenolith from Udachnaya, Siberia, preserves a pattern of crystallographic misorientation and subgrain microstructure associated with crystal–plastic deformation. The zircon grain records significant variations in titanium (Ti) from 2.6 to 30 ppm that corresponds to a difference in calculated Ti-in-zircon temperatures of over several hundred degrees Celsius. The highest Ti concentration is measured at subgrain centres (30 ppm), and Ti is variably depleted at low-angle boundaries (down to 2.6 ppm). Variations in cathodoluminescence coincide with the deformation microstructure and indicate localised, differential enrichment of rare earth elements (REE) at low-angle boundaries. Variable enrichment of U and Th and systematic increase of Th/U from 1.61 to 3.52 occurs at low-angle boundaries. Individual SHRIMP-derived U–Pb ages from more deformed zones (mean age of 1799 ± 40, n = 22) are systematically younger than subgrain cores (mean age of 1851 ± 65 Ma, n = 7), and indicate that open system behaviour of Ti–Th–U occurred shortly after zircon growth, prior to the accumulation of significant radiogenic Pb. Modelling of trace-element diffusion distances for geologically reasonable thermal histories indicates that the observed variations are ~ 5 orders of magnitude greater than can be accounted for by volume diffusion. The data are best explained by enhanced diffusion of U, Th and Ti along deformation-related fast-diffusion pathways, such as dislocations and low-angle (< 5°) boundaries. These results indicate chemical exchange between zircon and the surrounding matrix and show that Ti-in-zircon thermometry and U–Pb geochronology from deformed zircon may not yield information relating to the conditions and timing of primary crystallisation.     

Archean to Recent aeolian sand systems and their sedimentary record: Current understanding and future prospects

The sedimentary record of aeolian sand systems extends from the Archean to the Quaternary, yet current understanding of aeolian sedimentary processes and product remains limited. Most preserved aeolian successions represent inland sand‐sea or dunefield (erg) deposits, whereas coastal systems are primarily known from the Cenozoic. The complexity of aeolian sedimentary processes and facies variability are under‐represented and excessively simplified in current facies models, which are not sufficiently refined to reliably account for the complexity inherent in bedform morphology and migratory behaviour, and therefore cannot be used to consistently account for and predict the nature of the preserved sedimentary record in terms of formative processes. Archean and Neoproterozoic aeolian successions remain poorly constrained. Palaeozoic ergs developed and accumulated in relation to the palaeogeographical location of land masses and desert belts. During the Triassic, widespread desert conditions prevailed across much of Europe. During the Jurassic, extensive ergs developed in North America and gave rise to anomalously thick aeolian successions. Cretaceous aeolian successions are widespread in South America, Africa, Asia, and locally in Europe (Spain) and the USA. Several Eocene to Pliocene successions represent the direct precursors to the present‐day systems. Quaternary systems include major sand seas (ergs) in low‐lattitude and mid‐latitude arid regions, Pleistocene carbonate and Holocene–Modern siliciclastic coastal systems. The sedimentary record of most modern aeolian systems remains largely unknown. The majority of palaeoenvironmental reconstructions of aeolian systems envisage transverse dunes, whereas successions representing linear and star dunes remain under‐recognized. Research questions that remain to be answered include: (i) what factors control the preservation potential of different types of aeolian bedforms and what are the characteristics of the deposits of different bedform types that can be used for effective reconstruction of original bedform morphology; (ii) what specific set of controlling conditions allow for sustained bedform climb versus episodic sequence accumulation and preservation; (iii) can sophisticated four‐dimensional models be developed for complex patterns of spatial and temporal transition between different mechanisms of accumulation and preservation; and (iv) is it reasonable to assume that the deposits of preserved aeolian successions necessarily represent an unbiased record of the conditions that prevailed during episodes of Earth history when large‐scale aeolian systems were active, or has the evidence to support the existence of other major desert basins been lost for many periods throughout Earth history?     

Neoglacial (<3000 years) till and flutes at Saskatchewan Glacier,Canadian Rocky Mountains,formed by subglacial deformation of a soft bed

Understanding the processes that deposit till below modern glaciers provides fundamental information for interpreting ancient subglacial deposits. A process‐deposit‐landform model is developed for the till bed of Saskatchewan Glacier in the Canadian Rocky Mountains. The glacier is predominantly hard bedded in its upper reaches and flows through a deep valley carved into resistant Palaeozoic carbonates but the ice margin rests on a thick (<6 m) soft bed of silt‐rich deformation till that has been exposed as the glacier retreats from its Little Ice Age limit reached in 1854. In situ tree stumps rooted in a palaeosol under the till are dated between ca 2900 and 2700 yr bp and record initial glacier expansion during the Neoglacial. Sedimentological and stratigraphic observations underscore the importance of subglacial deformation of glaciofluvial outwash deposited in front of the advancing glacier and mixing with glaciolacustrine carbonate‐rich silt to form a soft bed. The exposed till plain has a rolling drumlinoid topography inherited from overridden end moraines and is corrugated by more than 400 longitudinal flute ridges which record deformation of the soft bed and fall into three genetically related types: those developed in propagating incipient cavities in the lee of large subglacial boulders embedded in deformation till, and those lacking any originating boulder and formed by pressing of wet till up into radial crevasses under stagnant ice. A third type consists of U‐shaped flutes akin to barchan dunes; these wrap around large boulders at the downglacier ends of longitudinal scours formed by the bulldozing of boulders by the ice front during brief winter readvances across soft till. Pervasive subglacial deformation during glacier expansion was probably facilitated by large boulders rotating within the soft bed (‘glacioturbation’).     

Application of the Ordered Weighted Averaging (OWA) method to the Caspian Sea conflict

This study proposes a promising allocation mechanism of the Caspian Sea natural resources, which are presently shared among five countries. To date, these nations have been unable to reach an allocation agreement. We apply a methodology to propose the most appropriate solution under different risk attitudes of the states. This research is different from other studies regarding the Caspian Sea negotiations in that it employs risk-based fuzzy multi attribute decision making methods for simulating the risk attitudes or optimism/pessimism degrees of the decision makers. The ordered weighted averaging (OWA) approach, which considers the optimism/pessimism degree quantitatively, is used to take into account the effects of different risk attitudes of the negotiators on the final outcome. We demonstrate how one could obtain a range of alternatives under different multi attribute and risk attitudes. The induced OWA (IOWA) method is also used to determine the relative power of these states bordering the Caspian Sea by considering several attributes, including different risk attitudes of agents. Results indicate that taking into account the risk attitude (prone, neutral, averse) of the states can affect the overall ranking of the proposed solutions. The findings from this study may facilitate negotiation regarding the most preferred allocation mechanism for the Caspian Sea.     

Brief note on effects of thermal pulses on elastic moduli and Q of rock

Preliminary results of a current experimental study on relationships between elastic moduli, internal friction (Q^?1) and rock structure are compatible with lunar in situ observations which indicate that mechanisms which lower the elastic moduli of rock need not contribute significantly to the internal friction. In terrestrial rock samples which were thermally pulsed in vacuum, elastic moduli decreased andQ's increased. These results imply that the thermal pulsing induces damage such as microcracks (which lower the elastic moduli), eliminates crack surface contaminants and may destroy other possible sites of internal friction such as grain boundary mechanisms and hydrous minerals (increasingQ).     

Using micro-catchment experiments for multi-local scale modelling of nature-based solutions

The Q-natural flood management project has co-developed with the Environment Agency 18 monitored micro-catchments (~1 km²) in Cumbria, UK installing calibrated flumes aimed at quantifying the potential shift in observed flows resulting from a range of nature-based-solutions installed by local organizations. The small-scale reduces the influence of variability characterizing larger catchments that would otherwise mask any such shifts, which we attempt to relate to a shift in model parameters. This paper demonstrates an approach to applying donor-parameter-shifts obtained from modelling two of the paired micro-catchments to a much larger scale, in order to understand the potential for improved distributed modelling of nature-based solutions in the form of additional tree-planting. The models include a rainfall-runoff model, Dynamic Topmodel, and a 2D hydrodynamic model, JFlow, permitting analysis of changes in hillslope processes and channel hydrodynamics resulting from a range of distributed measures designed to emulate natural hydrological processes that evaporate, store or infiltrate flows. We report on attempts to detect shift in hydrological response using one of the paired-micro-catchment moorland versus forestry sites in Lorton using Dynamic Topmodel. A donor-parameter-shift approach is used in a hypothetical experiment to represent new woodland in a much larger catchment, although testing all combinations of spatial planting strategies, responses to multiple-extremes, failure-modes and changes to synchronization becomes intractable to support good decision making. We argue that the problem can be re-framed to use donor-parameter-shifts at multi-local-scale catchments above communities known to be at risk, commensurate with most of the evidence of NbS impacts being effective at the small scale (ca. 10 km²). This might lead to more effective modelling to help catchment managers prioritize those communities-at-risk where there is more evidence that NbS might be effective.     

Zircon U-Pb-Hf constraints from Gongga Shan granites on young crustal melting in eastern Tibet

The Gongga Shan batholith is a complex granitoid batholith on the eastern margin of the Tibetan Plateau with a long history of magmatism spanning from the Triassic to the Pliocene. Late Miocene-Pliocene units are the youngest exposed crustal melts within the entire Asian plate of the Tibetan Plateau.Here, we present in-situ zircon Hf isotope constraints on their magmatic source, to aid the understanding of how these young melts were formed and how they were exhumed to the surface. Hf isotope signatures of Eocene to Pliocene zircon rims(ε_Hf(t)=-4 to +4), interpreted to have grown during localised crustal melting, are indicative of melting of a Neoproterozoic source region, equivalent to the nearby exposed Kangding Complex. Therefore, we suggest that Neoproterozoic crust underlies this region of the Songpan-Ganze terrane, and sourced the intrusive granites that form the Gongga Shan batholith. Localised young melting of Neoproterozoic lower or middle crust requires localised melt-fertile lithologies. We suggest that such melts may be equivalent to seismic and magnetotelluric low-velocity and high-conductivity zones or "bright spots" imaged across much of the Tibetan Plateau. The lack of widespread exposed melts this age is due either to the lack of melt-fertile rocks in the middle crust, the very low erosion level of the Tibetan plateau, or to a lack of mechanism for exhuming such melts. For Gongga Shan, where some melting is younger than nearby thermochronological ages of low temperature cooling, the exact process and timing of exhumation remains enigmatic, but their location away from the Xianshuihe fault precludes the fault acting as a conduit for the young melts. We suggest that underthrusting of dry granulites of the lower Indian crust(Archaean shield) this far northeast is a plausible mechanism to explain the uplift and exhumation of the eastern Tibetan Plateau.     

Alternative water management options to reduce vulnerability for climate change in the Netherlands

Urbanization, land subsidence and sea level rise will increase vulnerability of the urbanized low-lying areas in the western part of the Netherlands. In this article possibilities are explored to reduce vulnerability of these areas by implementing alternative water management options. Two main water management fields are distinguished, water supply and flood control. A four-component vulnerability framework is presented that includes threshold capacity, coping capacity, recovery capacity, and adaptive capacity. By using the vulnerability framework it is shown that current water supply and flood control strategies in the Netherlands focus on increasing threshold capacity by constructing higher and stronger dikes, improved water storage and delivery infrastructure. A complete vulnerability decreasing strategy requires measures that include all four capacities. Flood damage reduction, backup water supply systems and emergency plans are measures that can contribute to increasing coping capacity. Recovery capacity can be increased by multi-source water supply, insurance, or establishing disaster funds. Adaptive capacity can be developed by starting experiments with new modes of water supply and urbanization. Including all four components of the vulnerability framework enables better understanding of water and climate related vulnerability of urban areas and enables developing more complete water management strategies to reduce vulnerability.     

Finite Element Solvers for Biot’s Poroelasticity Equations in Porous Media

We study and compare five different combinations of finite element spaces for approximating the coupled flow and solid deformation system, so-called Biot’s equations. The permeability and porosity fields are heterogeneous and depend on solid displacement and fluid pressure. We provide detailed comparisons among the continuous Galerkin, discontinuous Galerkin, enriched Galerkin, and two types of mixed finite element methods. Several advantages and disadvantages for each of the above techniques are investigated by comparing local mass conservation properties, the accuracy of the flux approximation, number of degrees of freedom (DOF), and wall and CPU times. Three-field formulation methods with fluid velocity as an additional primary variable generally require a larger number of DOF, longer wall and CPU times, and a greater number of iterations in the linear solver in order to converge. The two-field formulation, a combination of continuous and enriched Galerkin function space, requires the fewest DOF among the methods that conserve local mass. Moreover, our results illustrate that three out of the five methods conserve local mass and produce similar flux approximations when conductivity alteration is included. These comparisons of the key performance indicators of different combinations of finite element methods can be utilized to choose the preferred method based on the required accuracy and the available computational resources.