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.     

Architecture of a distal pre-vegetation braidplain: Cambrian middle member of the Wood Canyon Formation,southern Marble Mountains,California, USA

Architectural element analysis and detailed mapping of a 300 m along-strike exposure of the middle member Wood Canyon Formation, southern Marble Mountains, California, USA, provides new evidence for extensive braided–fluvial channel-belt deposits with adjacent overbank environments. Three-dimensional models constructed using ‘Structure from Motion’ techniques, combined with field-based observations, allowed interpretation of outcrop-scale trends, barforms, channel fills and fine-scale features. The ca 80 m thick member is divisible into five distinct units, including units M1 to M3 that form the bulk of the stratigraphy. Units are defined by stacking patterns of three facies associations (Facies Association 1 to Facies Association 3), each representing the product of a subenvironment within the fluvial system. In Facies Association 1, stacked cosets, interpreted as low-relief fluvial bars and channel fills, preserve vertical-accretion and downstream-accretion elements under unimodal north-north-west palaeoflow, with minor lateral accretion near bar edges. Deposits of Facies Association 2 to Facies Association 3, linked to overbank environments, are found only in unit M2, in the middle 27 m of the middle member. Floodplains, represented by Facies Association 2, include crumbly red-orange intervals of fine to medium-grained sandstone and thinner sets of cross-bedding than Facies Association 1, interbedded with thicker cross-stratification indicative of overbank splay or overland flow aggradation from adjacent channel belts during flood stage. Possible aeolian beds of Facies Association 3 preserve broad festooned trough cross-strata that average 23 cm in thickness; their small size, medium-grained sandstone and iron oxide cement suggest a high water table. The diverse assemblage of interpreted subenvironments, paired with bedform and facies patterns, implies a perennial fluvial system that gradually built large sand bars as the channel belt migrated and avulsed across an unconfined braided–fluvial reach, leaving the overbank area on its flanks subject to weathering and aeolian transport. Despite the occurrence of strata deposited in low-energy and ponded settings, and a marine influence proposed for nearby sections of middle member, no ichnofossils were encountered.     

Unification Of Non-Dimensional Solutions To Asymptotic Equations For Plumes Of Different Shape

A model of buoyant fluid rise is developed in aLagrangian framework. Results from the model arevalidated by comparison with laboratory and fieldexperiments. The model is sufficiently general toinclude geometries of any cross-sectional shape. Exact solutions to the asymptotic equations forcontinuous discharges from line and point sources andfor an instantaneous discharge from a point source areconsidered. Prismatic, cylindrical andspherical shapes approximate these three geometries,respectively. Accommodation of these shapes withinthe same general model allows for direct comparison. It is shown that, for discharges into a cross-flowthat may be stratified or unstratified, thenon-dimensional trajectory, thickness and dilution canbe uniquely specified using three parameters. Theseare the non-dimensional size of the source, therelative importance of the initial fluxes of momentumand buoyancy and the number of orthogonal axes throughwhich entrainment can occur. Such non-dimensionalresults are particularly useful for examining thosescenarios for which there are limited experimentaldata.     

The rise of forced plumes in a stably stratified crossflow

Solutions are found describing the rise of turbulent buoyant plumes in a neutral or stably stratified crossflow from a source having finite fluxes of buoyancy, momentum and volume. Plumes from such sources are known as forced plumes. The solutions arise from a unified theory having one disposable parameter, the entrainment constant, and they describe a wide range of behaviour of plumes and jets in a crossflow. Solutions for buoyant plumes indicate that, for any given source, heights of rise are reduced with increased crossflow velocity and increased stratification while for constant environmental parameters, heights of rise are increased with increased fluxes of source buoyancy and momentum and are decreased with source radius. Maximum dilutions occur in plumes from small sources with relatively large buoyancy fluxes in light crossflows with small stratification.     

The anomalous tides near Broad Sound

Observations of tidal current and height, in conjunction with theoretical mathematical models are used to investigate the propagation of the tide near Broad Sound, a narrowing estuary situated on a wide section of continental shelf toward the southern end of the Great Barrier Reef. The observations indicate that the dense offshore reefs severely inhibit tidal flow, with the result that tides flood toward Broad Sound from the north and from the south, along the main lagoon. There is a local magnification of the semi-diurnal tides within Broad Sound itself.Models of flow across reefs confirm the effectiveness of dense, shallow, and broad reefs in acting as a barrier to the tide. The diffraction of tides through large gaps in the reef is modelled using conformal mapping techniques and with the inclusion of energy leakage, the diffraction model predicts magnification of the semi-diurnal tidal heights by a factor of about 4 and a phase lag of 3 h on the shelf near Broad Sound, these values being consistent with observation. The observed convergence of the tide close to, and within Broad Sound itself is consistent with the proximity of the semi-diurnal tidal period to the natural period for flow in Broad Sound, considered as a narrowing estuary. This results in further amplification, by an additional factor of about 1.5, so that the tides in Broad Sound are increased by a factor of between 5 and 6, altogether, compared with those elsewhere on the east Australian coast.     

Sea-floor spreading rate changes in the South Atlantic

A computer correlation technique was used to deduce the spreading history of the Mid-Atlantic Ridge from 5 magnetic profiles between 28°S and 43°S. In general, several possible histories are indicated for each profile involving changes of spreading rate and faulting, some of which are easily overlooked by the visual method. The only spreading history that was consistent will all the profiles required spreading at approximately 2.2 cm yr^-1 from 11 m.y.b.p. to approximately 5.5 m.y.b.p., followed by a decrease in rate to 1.7 cm yr^-1 relative to the Vine (1966) magnetic reversal model based on the South Pacific. Comparison of the data with other reported spreading rate discontinuities suggests that the South Pacific may be reponsible for the reported spreading rate changes.