Glaciological Research Articles

Applying the engineering geological model approach to quaternary glacial geology in western Tasmania

Applying the engineering geological model approach to quaternary glacial geology in western Tasmania

Recovering groundwater for wetlands from an anthropogenic aquifer

Freshwater wetlands are groundwater-dependent ecosystems that require groundwater for saturation, for wetland plants and creatures, for maintenance of wetland soils, and thermal buffering. With worldwide wetland area in decline for decades if not centuries, finding and restoring wetlands provides enormous ecosystem and public benefits, yet so often these projects fail to yield self-sustaining wetland ecosystems. One reason is that restored wetlands are often built in places that are neither wet enough nor possess the underlying geology to sustain them, and they dry out or require continual (expensive!) water inputs. Massachusetts is making the best of a challenging situation for the declining cranberry farming industry: while competition from less expensive land and more productive varietals shifts cranberry production to other locations, everything under historic cranberry farms is ripe for resilient wetland restoration projects. These low-lying water-rich areas are underlain by glacial geology (peats and clays) that are ideal for holding water, they possess historic seed banks of wetland plants and large accumulations of organic and hydric soils, and are currently sought-after by a statewide restoration program, for which these results provide critical information for restoration design, enabling practitioners to maximize the capture and residence time of groundwater inputs to sustain the future wetland. In this paper, we investigate the human legacy of cranberry farming on the surface of a wetland as it has created a unique hydrogeologic unit: the anthropogenic aquifer. Water moves through an anthropogenically constructed aquifer in specific and predictable ways that were engineered to favor a monoculture of cranberry plants on the surface of what once was a peatland. In order to restore this landscape to a functioning freshwater wetland, every property of the anthropogenic aquifer must be reversed. We detail observational, thermal, hydrologic, geologic and isotopic evidence for the location of groundwater inflows to Foothills Preserve in southeastern Massachusetts. The specific properties of the Anthropogenic aquifer, and the location and magnitude of groundwater discharge at this location provide crucial information for practitioners when designing plans for a self-sustaining, resilient restored freshwater wetland on this and future sites.

Imaging mass‐wasting sliding surfaces within complex glacial deposits along coastal cliffs using geophysics

AbstractThis study presents a multidisciplinary survey combining geological fieldwork and geophysical data to better constrain the parameters influencing the morphology and behaviour of a retreating coastal cliff. Erosion rates are spatially highly variable and hard to predict because of the manifold parameters acting on them. Among these parameters, rock resistance exerts a paramount influence on cliff retreat. Characterizing the rock resistance distribution along a coastal region requires the mapping of several key subsurface properties including the bulk lithology, faulting, fracturing, or weathering. This is a difficult and expensive task because of the high spatial variability of these factors linked to the spatial complexity of the geology. Geophysical methods can be used to tackle this challenge by quickly providing the 3D visualization and distribution of these parameters within the subsurface. A fast‐eroding portion of the Norfolk coast (UK) at West Runton is investigated using a multidisciplinary approach, combining ground‐penetrating radar, electrical resistivity tomography (ERT), cone penetration tests, and outcrop studies. The results allowed us to build a 3D geological and geophysical model of a highly complex area of glacial geology. It forms part of a relict glaciotectonic thrust‐tip moraine and sand basin sequence. The surfaces interpreted on radar data are associated with strong resistivity contrasts on the ERT data. These contrasts have been attributed to petrophysical variations between the lithological units. The base of the sand basin is marked by a low‐permeability clay bed. Its low shear strength is likely to be more susceptible to failure, hereby accelerating the erosion rate of an already fast‐eroding sand basin. The resulting model can be used as input for locally constraining the ground parameters in coastal recession and erosion models.

An Assessment of Geophysical Survey Techniques for Characterising the Subsurface Around Glacier Margins, and Recommendations for Future Applications

Geophysical surveys provide an efficient and non-invasive means of studying subsurface conditions in numerous sedimentary settings. In this study, we explore the application of three geophysical methods to a proglacial environment, namely ground penetrating radar (GPR), seismic refraction and multi-channel analysis of surface waves (MASW). We apply these geophysical methods to three glacial landforms with contrasting morphologies and sedimentary characteristics, and we use the various responses to assess the applicability and limitations of each method for these proglacial targets. Our analysis shows that GPR and seismic (refraction and MASW) techniques can provide spatially extensive information on the internal architecture and composition of moraines, but careful survey designs are required to optimise data quality in these geologically complex environments. Based on our findings, we define a number of recommendations and a potential workflow to guide future geophysical investigations in analogous settings. We recommend the initial use of GPR in future studies of proglacial environments to inform (a) seismic survey design and (b) the selection of seismic interpretation techniques. We show the benefits of using multiple GPR antenna frequencies (e.g., 25 and 100 MHz) to provide decimetre scale imaging in the near surface (e.g., < 15 m) while also enabling signal penetration to targets at up to ∼40 m depth (e.g., bedrock). This strategy helps to circumvent changes in radar signal penetration resulting from variations in substrate conductivity or abundant scatterers. Our study also demonstrates the importance of combining multiple geophysical methods together with ground-truthing through sedimentological observations to reduce ambiguity in interpretations. Implementing our recommendations will improve geophysical survey practice in the field of glacial geology and allow geophysical methods to play an increasing role in the interpretation of glacial landforms and sediments.

Effect of Composition on Engineering Behavior of Clay Tills

Glacial geology, marine environment, and arid climate govern the composition of clay tills. The main purpose of this work is to develop a clear understanding of the engineering behavior of compacted clay till under soil suction and applied stress. The results indicate moderate water adsorption due to the presence of clay minerals (26% corrensite, 10% illite, and 8% kaolinite) with Ca2+ as the dominant cation and a flocculated fabric in a slightly basic (pH = 7.5) pore water. The water retention curve comprised four transition points that are associated with capillary water drainage from large pores (air entry value of 2 kPa and residual suction value of 20 kPa) and small pores (air entry value of 700 kPa and residual suction value of 5 × 104 kPa). Beyond the last value, vapor flow is dominant and removes the adsorbed water by evaporation. The ratio of soil volume change to water volume change best described the s-shaped shrinkage path that also comprised four stages, namely: from most large pores with low volume change; from remaining large and most small pores along with almost equal volume change; from some small pores with low volume change; and from the rest of the small pores with no volume change. Likewise, the s-shaped swelling potential curve comprised three stages and correlated well with bimodal hydraulic conductivity curve, that is, slow initial swelling (unsaturated hydraulic conductivity around 10−14 m/s) due to expansion of peripheral clay in lumps; high primary swelling (unsaturated hydraulic conductivity of up to 10−11 m/s) due to thickened double layer of most colloids; and slow secondary swelling (albeit unsaturated hydraulic conductivity around 10−10 m/s) due to expansion of remaining particles. Soil compression (compression index of 0.164) was due to a gradual reduction in number of large pores, whereas rebound (swelling index of 0.047) was due to water adsorption on clay with part of the deformations recovered. Finally, the consolidation rate was related to saturated hydraulic conductivity, which varied by three orders of magnitude.

Lyell, the Geikies and Croll's observations on terrestrial glacial sediments and landforms

ABSTRACTWithin the context of the work and achievements of James Croll, this paper reviews the records of direct observations of glacial landforms and sediments made by Charles Lyell, Archibald and James Geikie and James Croll himself, in order to evaluate their contributions to the sciences of glacial geology and Quaternary environmental change. The paper outlines the social and physical environment of Croll's youth and contrasts this with the status and experiences of Lyell and the Geikies. It also outlines the character and role of the ‘Glasgow School’ of geologists, who stimulated Croll's interest into the causes of climate change and directed his focus to the glacial and ‘interglacial’ deposits of central Scotland. Contributions are outlined in chronological order, drawing attention to: (i) Lyell's high-quality observations and interpretations of glacial features in Glen Clova and Strathmore and his subsequent rejection of the glacial theory in favour of processes attributed to floating icebergs; (ii) the significant impact of Archibald Geikie's 1863 paper on the ‘glacial drift of Scotland’, which firmly established the land-ice theory; (iii) the fact that, despite James Croll's inherent dislike of geology and fieldwork, he provided high-quality descriptions and interpretations of the landforms and sediments of central Scotland in order to test his theory of climate change; and (iv) the great communication skills of James Geikie, enhanced by contacts and evidence from around the world. It is concluded that whilst direct observations of glacial landforms and sediments were critical to the long-term development of the study of glaciation, the acceptance of this theory was dependent also upon the skills, personality and status of the Geikies and Croll, who developed and promoted the concepts. Sadly, the subsequent rejection of the land-ice concept by Lyell resulted in the same factors challenging the acceptance of the glacial theory.

Glacial erosion by the Trift glacier (Switzerland): Deciphering the development of riegels, rock basins and gorges

A long-lasting question in glacial geology is how and how fast glaciers were able to shape the distinctive landscapes of the Alps. This study contributes to the understanding on the formation of overdeepened basins, especially the processes and the amount of time involved. We examine the remarkably high (150 m) cross-valley bedrock riegel and the associated overdeepening located in front of the Trift glacier in the central Swiss Alps. A combined approach of field survey with measurements of two cosmogenic nuclides, 10Be and in-situ 14C, and a numerical model was used to determine the spatial glacial erosion patterns on the bedrock riegel. Ten samples were taken along two transects; one perpendicular to the glacier flow direction, from outside of the Little Ice Age (LIA) extent down to the centre of the riegel, and the other following the former ice-flow direction across the riegel. Analysis of measured nuclide concentrations shows that the sample outside of the LIA was constantly exposed since the retreat of the Egesen stadial Trift glacier (~11.5 ka). The samples inside the LIA extent indicate a distinct trend of increasing glacial erosion rates from 0 mm/a near the LIA ice margin to high erosion all across the top of the riegel. The resulting minimum glacial erosion rates from samples on the riegel are 0.5–1.1 mm/a (10Be) and 0.6–>1.8 mm/a (in-situ 14C) which correspond to minimum erosion depths of 1.6–>3 m (10Be) and 1–>5 m (14C). The extremely low nuclide concentrations measured at the riegel highlight the substantial erosion (predominantly abrasion) of the bedrock surface during late Holocene glacier coverage. Field observations suggest that the formation of the overdeepening and, as a consequence, the riegel is due to a combination of valley shape, bedrock structures, glacier confluence and hydrology. We further hypothesise that the gorge is a key factor responsible for this impressive overdeepening, by lowering the threshold for the subglacial meltwater, effectively decoupling the height of the riegel from the depth of the overdeepening.

Zastosowanie ziemskich izotopów kosmogenicznych w datowaniu glin lodowcowych – potencjał, ograniczenia i perspektywy badawcze

Application of in situ cosmogenic nuclides in dating glacial till – research potential, limitations and prospects. A b s t r a c t. Determination of numerical age for glacial till is one of the most intriguing issues in glacial geology and geomorphology. Till layers are unequivocal evidence of ice sheets advances and retreats, and therefore they are very important for Quaternary palaeoclimatic reconstructions. Herein, we present a review dedicated to the methodology of glacial till dating with in situ cosmogenic nuclides. This relatively novel approach enables direct dating of deposition of till and its exposure after deglaciation. Getting reliable numerical age (in years) of till opens a new prospect for reconstructions of past glaciations and palaeogeographic studies. Despite some limitations of applying in situ cosmogenic exposure and burial dating, we argue that this technique may be successfully employed in studies of glacial deposits in Poland.

Developing 3D geological and hydrogeological models for the Low Level Waste Repository site, west Cumbria, UK

This case study presents the recent development of new 3D geological models for the Low Level Waste Repository in Cumbria through the review of historical and new datasets. The site is located in an area of complex Late Pleistocene glacial geology and geomorphology, which control near-surface groundwater movement, and benefits from extensive historical site characterization data. This work emphasizes the value of a quality review of existing borehole records, drilling techniques and supplementary information to build a predictive geological model. The integration of different datasets has proved crucial to our understanding of the ground conditions, which are characterized by complex Quaternary sequences, and the depositional processes responsible for their formation. The updated geological interpretation suggests that the Low Level Waste Repository is located on a former ice-dammed lake plain. This is contrary to the previous interpretation of multiple and extensive till units in the area. Although the permeabilities may be similar in both interpretations, the extents may be different. Understanding the depositional processes and extent and the hydrogeological characteristics of the deposits is key to developing a hydrogeological model that is appropriate for the assessment of the environmental safety of the site, both now and into the future. Thematic collection: This article is part of the Ground models in engineering geology and hydrogeology collection available at: https://www.lyellcollection.org/cc/Ground-models-in-engineering-geology-and-hydrogeology

Major Ice Sheet Change in the Weddell Sea Sector of West Antarctica Over the Last 5,000 Years

AbstractUntil recently, little was known about the Weddell Sea sector of the West Antarctic Ice Sheet. In the last 10 years, a variety of expeditions and numerical modelling experiments have improved knowledge of its glaciology, glacial geology and tectonic setting. Two of the sector's largest ice streams rest on a steep reverse‐sloping bed yet, despite being vulnerable to change, satellite observations show contemporary stability. There is clear evidence for major ice sheet reconfiguration in the last few thousand years, however. Knowing precisely how and when the ice sheet has changed in the past would allow us to better understand whether it is now at risk. Two competing hypotheses have been established for this glacial history. In one, the ice sheet retreated and thinned progressively from its Last Glacial Maximum position. Retreat stopped at, or very near, the present position in the Late Holocene. Alternatively, in the Late Holocene, the ice sheet retreated significantly upstream of its present grounding line and then advanced to today's location due to glacial isostatic adjustment and ice shelf and ice rise buttressing. Both hypotheses point to data and theory in their support yet neither has been unequivocally tested or falsified. Here we review geophysical evidence to determine how each hypothesis has been formed, where there are inconsistencies in the respective glacial histories, how they may be tested or reconciled, and what new evidence is required to reach a common model for the Late Holocene ice sheet history of the Weddell Sea sector of West Antarctica.

Determining the ice volume distribution of glaciers using helicopter-borneGPR and glaciological modeling

Determining the ice volume distribution of glaciers using helicopter-borneGPR and glaciological modeling

Glacial geology and cosmogenic-nuclide exposure ages from the Tucker Glacier - Whitehall Glacier confluence, northern Victoria Land, Antarctica

We describe glacial-geological observations and cosmogenic-nuclide exposure ages from the vicinity of the present grounding line of Tucker Glacier, a large alpine glacier flowing from the mountains of northern Victoria Land, East Antarctica, into the outer Ross Sea. These data are relevant for constraining the extent of ice sheet expansion and retreat in the Ross Sea, and associated eustatic sea level impact, between the Last Glacial Maximum (LGM) and the present. In addition, a terrestrial geological record of ice thickness change from this region could provide evidence for or against the hypothesis that rapid eustatic sea-level rise during meltwater pulse 1A (“MWP-1A”) at 14.6 ka was in part the result of rapid, large-scale thinning or breakup of a marine-based portion of the LGM ice sheet in the outer Ross Sea. Glacial-geological observations, exposure ages on glacial deposits, and a novel application of in-situ-produced cosmogenic 14C in quartz-bearing bedrock to identify the limits of LGM ice cover in the absence of direct geomorphic evidence, show that Tucker Glacier near its present grounding line was 300 to 350 m thicker than present during the LGM and thinned steadily between 17 to 5 ka. The largest possible rapid thickness change in the time period 14 to 15 ka that could be accommodated by the exposure-age data is ∼50 m, which is a small fraction of that predicted for the western Ross Sea by model simulations of the Antarctic contribution to MWP-1A. There do exist possible scenarios in which hypothesized marine ice sheet collapse in the outer Ross Sea during MWP-1A might not be recorded by ice thickness changes at Tucker Glacier. However, our record of ice thickness changes spanning this time period is the closest such record to the outer Ross Sea that is likely to exist, and it agrees with all exposure-age deglaciation chronologies from other regions of the Ross embayment in providing no evidence for such an event.

Polar Conference and Glaciological Symposium in Sochi, September 2018

.

Glacial geology of Saalian relief around Midwolda, eastern Groningen, the Netherlands

AbstractDue to canal-digging activities in 2011 and 2014, two small and one large temporary exposure, all ranging from 4 to 5 m in depth, were studied with respect to the sedimentology and structural geology, in the glacial ridge of Midwolda, Groningen, the Netherlands. The lowermost unit consists of clay of Elsterian age and is composed of glaciolacustrine and turbiditic deposits (Peelo Formation). These show synsedimentary deformations due to loading, as well as post-sedimentary Saalian glaciotectonic deformations, consisting of folding, and faulting structures. The overlying Saalian till sequence consists of two main units. The lower unit, with clear features of a subglacial deformation zone (e.g. lateral heterogeneity), has a local origin and strongly resembles the underlying Elsterian clay. Glacial tectonic and morphological observations indicate a primary NE–SW ice-flow direction. The second till layer has a sandy texture and high crystalline gravel content, while glacial-tectonic indicators point to a NW–SE ice-flow direction. The deformation of the till layers has caused a repetition and mixing of till layers, due to the last ice movement. The NW–SE ice movement is supported by the morphology as well as data from erratic gravel counts. Correlation with geological cross-sections strongly suggests regional subsurface control on ice-sheet behaviour.

Applicability of an ultra-long-range terrestrial laser scanner to monitor the mass balance of Muz Taw Glacier, Sawir Mountains, China

Glacier mass balance is a key component of glacier monitoring programs. Information on the mass balance of Sawir Mountains is poor due to a dearth of in-situ measurements. This paper introduces the applicability of an ultra-long-range terrestrial laser scanner (TLS) to monitor the mass balance of Muz Taw Glacier, Sawir Mountains, China. The Riegl VZ ® -6000 TLS is exceptionally well-suited for measuring snowy and icy terrain. Here, we use TLS to create repeated high spatiotemporal resolution DEMs, focusing on the annual mass balance (June 2, 2015 to July 25, 2016). According to TLS-derived high spatial resolution point clouds, the front variation (glacier retreat) of Muz Taw Glacier was 9.3 m. The mean geodetic elevation change was 4.55 m at the ablation area. By comparing with glaciological measurements, the glaciological elevation change of individual stakes and the TLS-derived geodetic elevation change of corresponding points matched closely, and the calculated balance was −3.864±0.378 m w.e.. This data indicates that TLS provides accurate results and is therefore suitable to monitor mass balance evolution of Muz Taw Glacier.

Statistical Models to Predict and Assess Spatial and Temporal Low‐Flow Variability in New England Rivers and Streams

In the northern hemisphere, summer low flows are a key attribute defining both quantity and quality of aquatic habitat. I developed one set of models for New England streams/rivers predicting July/August median flows averaged across 1985 to 2015 as a function of weather, slope, % imperviousness, watershed storage, glacial geology and soils. These models performed better than most USGS models for summer flows developed at a statewide scale. I developed a second set of models predicting interannual differences in summer flows as a function of differences in air temperature, precipitation, the North Atlantic Oscillation Index (NAO), and lagged NAO. Use of difference equations eliminated the need for transformations and accounted for serial autocorrelations at lag 1. The models were used in sequence to estimate time series for monthly low flows and for two derived flow metrics (tenth percentile (Q10) and minimum 3-in-5 year average flows). The first metric is commonly used in assessing risk to low flow conditions over time while the second has been correlated with increased probability of localized extinctions for brook trout. The flow metrics showed increasing trends across most of New England for 1985-2015. However, application of summer flow models with average and extreme climate projections to the Taunton River, MA, a sensitive watershed undergoing rapid development, projected that low flow metrics will decrease over the next 50 years.

Mountain glaciers under a changing climate

Mountain glaciers are found around the world in ranges such as the Himalaya, the Andes and the European Alps. The majority of mountain glaciers world‐wide are shrinking. However, the rugged alpine topography through which these glaciers flow governs their dynamics and impacts on the regional climate systems that modify glacier mass balance. As a result, the response of mountain glaciers to climate change is difficult to predict, and highly spatially variable even across one mountain range, particularly where orography controls precipitation distributions. To understand how mountain glaciers behave and change, geologists combine many different techniques based on direct observations and dating of glacial geology, measurements of present‐day glaciers, and predictive numerical (computer) models. Recent advances in these techniques and their applications to glacial environments have demonstrated that the glacial geological record is a rich archive of information about how climate has changed in the past, and gives greater confidence in predictions of glacier change in the future, which is required if populations living in glacerised catchments are able to adapt to the rapid response of glaciers to a changing climate.

Global Research Patterns on Ground Penetrating Radar (GPR)

The article deals with the analysis of worldwide research patterns concerning ground penetrating radar (GPR) during 1995–2014. To do this, the Thomson Reuters’ Science Citation Index Expanded (SCI-EXPANDED) and the Social Sciences Citation Index accessed via the Web of Science Core Collection were the two bibliographic databases taken as a reference. We pay attention to the document typology and language, the publication trend and citations, the subject categories and journals, the collaborations between authors, the productivity of the authors, the most cited articles, the countries and the institutions involved, and other hot issues. Concerning the main research subfields involving GPR use, there were five, physical–mathematical, sedimentological–stratigraphical, civil engineering/engineering geology/cultural heritage, hydrological (HD), and glaciological (GL), subfields.

Tool wear analysis of pressurized face TBM drives in the glacial geology of the Pacific Northwest

AbstractFor underground construction projects in the United States and Canada it is standard procedure to use a Geotechnical Baseline Report (GBR) to contractually define subsoil conditions. The GBR sets baselines based on which tunneling contractors develop bids and plan the works. Baseline values for soil abrasiveness are a focus especially where drives with pressurized‐face Tunnel Boring Machines (TBM) beneath the groundwater table and in unstable face conditions require changing the cutterhead tools under hyperbaric conditions or in pre‐constructed safe havens. Several laboratory procedures exist that can be used for providing soil abrasiveness baselines in the context of the GBR. However, none of them cover all the soil characteristics that are relevant in causing tool wear. Also, other factors need to be considered for wear rate prediction. Analyzing the performance of previous TBM drives is a proven way to gain insight into the wear system behavior. This paper presents correlation analyses of geotechnical conditions, TBM operational data, and tool wear measurements from several TBM drives in the metropolitan areas of Seattle and Vancouver, B.C. These drives with earth pressure balance and slurry TBMs include various tool types and were conducted in glacial and interglacial deposits that are considered highly abrasive.

Background Concentrations of Arsenic in Ohio Soils: Sources and Influencing Factors

Arsenic concentrations—representative of natural background conditions in Ohio soils—were compiled from the USGS National Geochemical Survey database, Ohio EPA studies, and private environmental site investigations. The resulting dataset of 2,783 samples collected from 1,116 locations was then correlated with soil data from USDA-NRCS SSURGO, glacial and bedrock geology, and analyses of other metals at these specific locations to identify factors most likely to affect arsenic concentrations in soil. Bedrock geology, particularly Devonian-aged materials and black shales, was found to significantly correlate to arsenic concentrations in soil. However, this correlation is complicated by both glaciation and post-glaciation erosion. Approximately 70% of Ohio’s bedrock is covered in glacial materials, such as till and outwash. As glaciers advanced across Ohio, bedrock materials were eroded and deposited farther south. Arsenic concentrations in Ohio soils tend to be highest where Devonian bedrock materials were deposited by glaciation. Following glacial deposition of materials, arsenic can be eroded or leached from the deposits. Deposits of wind-blown loess and materials deposited by water (as opposed to ice) had lower arsenic concentrations than glacial tills and outwash plains that underwent minimal erosion following glacial deposition.