A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

The Effect of Undetected Barriers on Groundwater Drawdown and Recovery

In large-scale pumping projects, such as mine dewatering, predictions are often made about the rate of groundwater level recovery after pumping has ceased. However, these predictions may be impacted by geological uncertainty—including the presence of undetected impermeable barriers. During pumping, an impermeable barrier may be undetected if it is located beyond the maximum extent of the cone of depression; yet it may still control drawdown during the recovery phase. This has implications for regional-scale modeling and monitoring of groundwater level recovery. In this article, non-dimensional solutions are developed to show the conditions under which a barrier may be undetected during pumping but still significantly impact groundwater level recovery. The magnitude of the impact from an undetected barrier will increase as the ratio of pumping rate to aquifer transmissivity increases. The results are exemplified for a hypothetical aquifer with an unknown barrier 3 km from a pumping well. The difference in drawdown between a model with and without a barrier may be <1 m in the 10 years while pumping is occurring, but up to 50 m after pumping has ceased.     

Sensitivity of reflection seismic data to oil-column height in high-porosity sandstones

Tuning is the effect of interference between the reflections from the top and bottom of a thin layer on the amplitude of the composite reflection. For a homogeneous sandstone reservoir containing an oil column overlying brine, interference between the reflection from the top reservoir and the oil/water contact is a function of the height of the oil column. If the properties of the sandstone do not vary across the oil/water contact, the SS, PS and SP reflection coefficients from the oil/water contact are small in comparison to the PP reflection coefficient. This allows analytic expressions for the effective PP and PS reflection coefficients from the reservoir to be derived that include all P‐wave multiples within the oil column. For a given source/receiver offset, the component of the wavevector inside the oil column normal to the interface is larger for the PPPP reflection than for the PPPS reflection, due to the asymmetry in the raypath for the PPPS reflection. The PPPS reflection is therefore useful for determining oil‐column heights larger than that discriminated by the PPPP reflection, especially when used at wider offsets. A convenient classification of the AVO response of hydrocarbon‐bearing sandstone reservoirs overlain by shale is the scheme of Rutherford and Williams. Class 1 sands have higher acoustic impedance for normal incidence than the overlying shale, Class 2 sands have nearly the same acoustic impedance as the shale and Class 3 sands have lower acoustic impedance. Synthetic shot gathers calculated for these three classes as a function of oil‐column height show that a combination of the PPPP and the PPPS amplitudes can be plotted as a tuning trajectory, which can be used to determine the oil‐column height. This method is most sensitive for reservoirs that belong to AVO classes 1 and 2, and therefore may be useful in AVO analysis of Class 1 and 2 reservoirs where the traditional AVO indicators (developed for Class 3 reservoirs) do not work very well. These results demonstrate the usefulness of shear waves recorded in the marine environment at wide offsets.     

Sulfur Hexafluoride and Potassium Bromide as Groundwater Tracers for Managed Aquifer Recharge

Sulfur hexafluoride (SF₆) is an established tracer for use in managed aquifer recharge projects. SF₆ exsolves from groundwater when it encounters trapped air according to Henry's law. This results in its retardation relative to groundwater flow, which can help determine porous media saturation and flow dynamics. SF₆ and the conservative, nonpartitioning tracer, bromide (Br⁻ added as KBr), were introduced to recharge water infiltrated into stacked glacial aquifers in Thurston County, Washington, providing the opportunity to observe SF₆ partitioning. Br⁻, which is assumed to travel at the same velocity as the groundwater, precedes SF₆ at most monitoring wells (MWs). Average groundwater velocity in the unconfined aquifer in the study area ranges from 3.9 to 40 m/d, except in the southwestern corner where it is slower. SF₆ in the shallow aquifer exhibits an average retardation factor of 2.5 ± 3.8, suggesting an air-to-water ratio on the order of 10⁻³ to 10⁻² in the pore space. Notable differences in tracer arrival times at adjacent wells indicate very heterogeneous conductivity. One MW exhibits double peaks in concentrations of both tracers with different degrees of retardation for the first and second peaks. This suggests multiple flowpaths to the well with variable saturation. The confining layer between the upper two aquifers appears to allow intermittent connection between aquifers but serves as an aquitard in most areas. This study demonstrates the utility of SF₆ partitioning for evaluating hydrologic conditions at prospective recharge sites.     

A Simple Method for Simulating Groundwater Interactions with Fens to Forecast Development Effects

Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.     

Bulk density of ordinary chondrite meteorites and implications for asteroidal internal structure

Abstract— The bulk densities of 82 samples of 72 ordinary chondrites (OCs) were measured to an accuracy of ～1% using a modified Archimedian method. We found the average bulk density to be 3.44 ± 0.19 g/cm³ for the H group, 3.40 ± 0.15 g/cm³ for the L group, and 3.29 ± 0.17 g/cm³ for the LL group (± represent 1σ). Bulk density measurements of 11 pieces of one fall, Pultusk (H group), were also found to vary considerably (3.31 to 3.63 g/cm³). To investigate controls on bulk density within the OCs, we compared density with bulk chemical composition (the ratio of Fe metal to total Fe; the ratio of total Fe to SiO₂; the ratio of FeO to total Fe and MgO). Within each OC group, bulk chemical composition is nearly invariant whereas bulk density varies from about 3.0 to 3.8 g/cm³. Slight but systematic differences in average density between the H, L, and LL groups presumably relate to differences in metallic Fe abundance. However, considerable overlap between OC groups and the wide range of bulk densities within each group suggest differences in porosity dominantly control variations of density within the OC subgroups.     

Artifact size and spatial process: Macro- and microartifacts in a mississippian house

Despite their evident utility in archaeological analysis, microartifacts (those artifacts smaller than 2 mm) have, as a class, been used only sparingly by archaeologists and then only as if they were larger artifacts. This article explores the variable information microartifacts contain, using a case study from the Loy Site, a Mississippian village in East Tennessee. We show that microartifacts provide different information from that obtained from macroartifacts. We also demonstrate that comparisons among different size classes provide new information on dynamic aspects of artifact deposition. Especially in spatial analysis, comparing distributions of different size artifacts enhances our interpretative capabilities over any single artifact size class.     

Magmatic 87Sr/86Sr relicts in hydrothermally altered quartz diorites (Brabant Massif,Belgium) and the role of epidote as a Sr filter

The porphyritic quartz diorites of the Caledonian Brabant Massif have been totally altered. Ca, Rb, Sr, Zr, Ce, Y measurements and Sr-Nd isotopic analyses were performed on the Quenast plug and the Lessines sill, in an attempt to study the relative mobility of Sr and evaluate the extent, direction and magnitude of the ⁸⁷Sr/⁸⁶Sr alterations. Sr electron microprobe analyses of epidote were also carried out to assess its role in the Sr distribution.The initial ⁸⁷Sr/⁸⁶Sr ratio is shown to have had an unsteady behaviour during the studied water/rock interactions since it has been sometimes enhanced, sometimes depressed and occasionally not modified. The possibility and magnitude of the ⁸⁷Sr contamination turn out to be strictly related to the degree of Sr accommodation in the secondary minerals. Epidote in particular has proved to be the main trap for the hydrothermal Sr and this mineral is thus regarded as the major controlling factor of ⁸⁷Sr hydrothermal contamination. The epidote-poor rocks (albite+chlorite-rich rocks) seem to have been unaffected by any Sr interchange with the aqueous solutions. Therefore, as alteration quickly follows the crystallization of the magma, their initial ⁸⁷Sr/ ⁸⁶Sr ratio, which is deduced from an isochron, might be a primary petrogenetic feature enabling interpretation of the genesis of their parental magmas. On the other hand, in the epidote-rich rocks, this ratio has been readily altered; it could thus generally be used only to trace the origin of the hydrothermal solutions. As a consequence, these rocks should not be selected for dating an alteration event by the Rb-Sr method.     

Dialogue, metaphors of dialogue and understandings of geography

Increasingly, human and physical geographers alike describe their research practices as involving dialogue of some kind. However, the widespread popularity of the term belies some very different understandings of its meaning and methodological implications. In this paper we reflect on these different conceptions of 'dialogical' geography as a way, first, of illustrating the broad range of understandings of research methodology now current in the discipline and, second, of identifying their implications for three long–standing controversies in the discipline: relativism and the truth of geographical knowledge, the associated dualisms between subject/object, nature/society, and realism/constructionism, and the potential for unity between human and physical geography. We argue that, while dialogue is a potentially fruitful way of understanding and practising geography, that the defence of it, like the attacks on it, is often misconceived.