The influence of groundwater abstraction on interpreting climate controls and extreme recharge events from well hydrographs in semi-arid South Africa

There is a scarcity of long-term groundwater hydrographs from sub-Saharan Africa to investigate groundwater sustainability, processes and controls. This paper presents an analysis of 21 hydrographs from semi-arid South Africa. Hydrographs from 1980 to 2000 were converted to standardised groundwater level indices and rationalised into four types (C1–C4) using hierarchical cluster analysis. Mean hydrographs for each type were cross-correlated with standardised precipitation and streamflow indices. Relationships with the El Niño–Southern Oscillation (ENSO) were also investigated. The four hydrograph types show a transition of autocorrelation over increasing timescales and increasingly subdued responses to rainfall. Type C1 strongly relates to rainfall, responding in most years, whereas C4 notably responds to only a single extreme event in 2000 and has limited relationship with rainfall. Types C2, C3 and C4 have stronger statistical relationships with standardised streamflow than standardised rainfall. C3 and C4 changes are significantly (p <?0.05) correlated to the mean wet season ENSO anomaly, indicating a tendency for substantial or minimal recharge to occur during extreme negative and positive ENSO years, respectively. The range of different hydrograph types, sometimes within only a few kilometres of each other, appears to be a result of abstraction interference and cannot be confidently attributed to variations in climate or hydrogeological setting. It is possible that high groundwater abstraction near C3/C4 sites masks frequent small-scale recharge events observed at C1/C2 sites, resulting in extreme events associated with negative ENSO years being more visible in the time series.

     

WIP: A Web‐based program for indexing planar features in quartz grains and its usage

Planar deformation features (PDFs) in quartz are the most important diagnostic features that allow the unambiguous identification of impact structures on Earth. In order to confirm that these features (that are characterized by planar character and form along specific crystallographic planes) are indeed PDFs, they need to be properly investigated and indexed. Following universal‐stage measurements, the process of indexing is usually performed manually, using a Wulff stereonet and following a strict procedure, which is time consuming and error prone. In this article, we present WIP, a new Web‐based program for indexing planar deformation features in quartz. The correctness of our program is shown by its application to measurements that had previously been indexed manually. The observed minor differences, especially in the absolute frequency percentage of PDFs, are negligible and not significant enough to influence the estimation of shock pressure that could be calculated from the indexed results. Usability of this program is shown using the spatial relationships between a statistically significant number of 278 quartz grains with 409 sets of PDFs analyzed within the area (~35 mm²) of a single thin section of a meta‐greywacke from the Bosumtwi impact structure. Our program is not only more accurate and faster than the manual (graphical) method but also removes the human error from the plotting process and allows control of several parameters, such as the value of estimated measurement error used in the indexing calculation or method of aggregated error handling. The program also provides information about the angles between the planes of the measured PDF sets present in a grain, which allows determination of the angles between (for example) indexed {} and {} sets.     

Heterogeneity preserving upscaling for heat transport in fractured geothermal reservoirs

In simulation of fluid injection in fractured geothermal reservoirs, the characteristics of the physical processes are severely affected by the local occurence of connected fractures. To resolve these structurally dominated processes, there is a need to develop discretization strategies that also limit computational effort. In this paper, we present an upscaling methodology for geothermal heat transport with fractures represented explicitly in the computational grid. The heat transport is modeled by an advection-conduction equation for the temperature, and solved on a highly irregular coarse grid that preserves the fracture heterogeneity. The upscaling is based on different strategies for the advective term and the conductive term. The coarse scale advective term is constructed from sums of fine scale fluxes, whereas the coarse scale conductive term is constructed based on numerically computed basis functions. The method naturally incorporates the coupling between solution variables in the matrix and in the fractures, respectively, via the discretization. In this way, explicit transfer terms that couple fracture and matrix solution variables are avoided. Numerical results show that the upscaling methodology performs well, in particular for large upscaling ratios, and that it is applicable also to highly complex fracture networks.     

Differential Effects of Bivalves on Sediment Nitrogen Cycling in a Shallow Coastal Bay

In coastal ecosystems, suspension-feeding bivalves can remove nitrogen though uptake and assimilation or enhanced denitrification. Bivalves may also retain nitrogen through increased mineralization and dissimilatory nitrate reduction to ammonium (DNRA). This study investigated the effects of oyster reefs and clam aquaculture on denitrification, DNRA, and nutrient fluxes (NO _x , NH₄ ⁺, O₂). Core incubations were conducted seasonally on sediments adjacent to restored oyster reefs (Crassostrea virginica), clam aquaculture beds (Mercenaria mercenaria) which contained live clams, and bare sediments from Smith Island Bay, Virginia, USA. Denitrification was significantly higher at oyster reef sediments and clam aquaculture site than bare sediment in the summer; however, DNRA was not enhanced. The clam aquaculture site had the highest ammonium production due to clam excretion. While oyster reef and bare sediments exhibited seasonal differences in rate processes, there was no effect of season on denitrification, or dissimilatory nitrate reduction to ammonium (DNRA) or ammonium flux at the clam aquaculture site. This suggests that farm management practices or bivalve metabolism and excretion may override seasonal differences. When water column nitrate concentration was elevated, denitrification increased in clam aquaculture site and oyster reef sediments but not in bare sediment; DNRA was only stimulated at the clam aquaculture site. This, along with a significant and positive relationship between denitrification and sediment organic matter, suggests that labile carbon limited nitrate reduction at the bare sediment site. Bivalve systems can serve as either net sinks or sources of nitrogen to coastal ecosystems, depending mainly on the type of bivalve, location, and management practices.     

Analysis of thermal parameters of Triassic,Permian and Carboniferous sedimentary rocks in central Poland

The paper shows thermal tests results (thermal conductivity, specific heat) carried out during the research project to recognize the potential of hot dry rocks for heat and electricity production in Poland. Analysis of such parameters of the rock medium has fundamental significance in search of structures for location of enhanced geothermal systems (EGS). Analyses of thermal conductivity as well as specific heat determine the possibility of effective extraction of heat from hot rock formations and its efficient transfer to the land surface. Reservoir rock should be characterized by the lowest possible porosity and permeability, and the highest possible thermal conductivity. In order to recognize reservoir parameters of rocks that form potential reservoirs for EGS, 300 samples of sedimentary rocks were taken from 11 wells located in central Poland. For samples, sizes of which enabled measurement of thermal parameters, tests of thermal conductivity and specific heat were carried out. Independently, measurements of porosity were made. The porosity magnitude determines values of measured thermal parameters. Measurements of thermal conductivity of rocks were made for 24 samples collected from Triassic, Permian and Carboniferous deposits in six wells. Measurements of specific heat on rock samples were made for 20 rock samples. Among sedimentary rocks, principally sandstones and limestones with reservoir parameters favourable for this type of systems are considered to be petrogeothermal reservoirs and are often characterized by favourable thermal parameters for EGS. Copyright © 2014 John Wiley & Sons, Ltd.     

Late Devensian deglaciation of the Tyne Gap Palaeo‐Ice Stream,northern England

The deglacial history of the central sector of the last British–Irish Ice Sheet is poorly constrained, particularly along major ice‐stream flow paths. The Tyne Gap Palaeo‐Ice Stream (TGIS) was a major fast‐flow conduit of the British–Irish Ice Sheet during the last glaciation. We reconstruct the pattern and constrain the timing of retreat of this ice stream using cosmogenic radionuclide (¹⁰Be) dating of exposed bedrock surfaces, radiocarbon dating of lake cores and geomorphological mapping of deglacial features. Four of the five ¹⁰Be samples produced minimum ages between 17.8 and 16.5 ka. These were supplemented by a basal radiocarbon date of 15.7 ± 0.1 cal ka BP, in a core recovered from Talkin Tarn in the Brampton Kame Belt. Our new geochronology indicates progressive retreat of the TGIS from 18.7 to 17.1 ka, and becoming ice free before 16.4–15.7 ka. Initial retreat and decoupling of the TGIS from the North Sea Lobe is recorded by a prominent moraine 10–15 km inland of the present‐day coast. This constrains the damming of Glacial Lake Wear to a period before ∼18.7–17.1 ka in the area deglaciated by the contraction of the TGIS. We suggest that retreat of the TGIS was part of a regional collapse of ice‐dispersal centres between 18 and 16 ka.

     

Expert engagement in participatory processes: translating stakeholder discourses into policy options

This paper demonstrates an innovative role for experts in supporting participatory policy processes with an application to landslide risk management in the Italian town of Nocera Inferiore. Experts co-produce risk mitigation options based on their specialized knowledge taking account of local knowledge and values by directly coupling stakeholder discourses with option design. Drawing on the theory of plural rationality and based on a literature review, interviews and a public questionnaire, stakeholder discourses are elicited on the landslide risk problem and its solution. Armed with the discourses and in close interaction with stakeholders, experts provide a range of technical mitigation options, each within a given budget constraint. These options are subsequently deliberated in the participatory process with the intent of reaching compromise recommendations for landslide risk mitigation. As we show in an accompanying paper, “Compromise not consensus. Designing a participatory process for landslide risk mitigation” (this issue), the provision of multiple co-produced policy options enhances stakeholder deliberation by respecting legitimate differences in values and worldviews.

     

Optimal design of an equipment isolation system with nonlinear hysteretic behaviour

Vibration isolation is well recognized as an effective mitigation strategy for acceleration‐sensitive equipment subjected to earthquake. In the present paper, an equipment isolation system with nonlinear hysteretic behaviour is proposed and a methodology for the optimal design is developed. An integrable constitutive model, derived from the mathematical Duhem hysteresis operator, is adopted for the isolation system. The optimization procedure is defined through a dual‐criteria approach that involves a transmissibility criterion combined with an energy performance criterion: the former consists in limiting the absolute acceleration of the isolated equipment below an allowable threshold value; the latter, in maximizing the ratio between the energy dissipation due to hysteresis and the input energy to reduce the isolator displacements. The seismic effectiveness of the nonlinear hysteretic isolation system is numerically investigated under natural accelerograms with different frequency content and increasing levels of excitation. Both ground‐mounted and floor‐mounted equipment items are considered in the analyses; in the second case, the dynamic interaction between the equipment and its supporting structure is taken into account in the design of the isolation system, and its effects on the isolation performance and the structural response are discussed. Comparisons in terms of effectiveness and robustness with a linear isolation system with viscoelastic behaviour are eventually provided. Copyright © 2013 John Wiley & Sons, Ltd.