Geology topics

Daniel Mcgrath

Publications and source records attributed to Daniel Mcgrath.

12 recordsLinked to original sources

USGS2018

A snow density dataset for improving surface boundary conditions in Greenland ice sheet firn modeling

The surface snow density of glaciers and ice sheets is of fundamental importance in converting volume to mass in both altimetry and surface mass balance studies, yet it is often poorly constrained. Site-specific surface snow densities are typically derived from empirical relations based on temperature and wind speed. These parameterizations commonly calculate the average density of the top meter of snow, thereby systematically overestimating snow density at the actual surface. Therefore, constraining surface snow density to the top 0.1 m can improve boundary conditions in high-resolution firn-evolution modeling. We have compiled an extensive dataset of 200 point measurements of surface snow density from firn cores and snow pits on the Greenland ice sheet. We find that surface snow density within 0.1 m of the surface has an average value of 315 kg m −3 with a standard deviation of 44 kg m −3 , and has an insignificant annual air temperature dependency. We demonstrate that two widely-used surface snow density parameterizations dependent on temperature systematically overestimate surface snow density over the Greenland ice sheet by 17–19%, and that using a constant density of 315 kg m −3 may give superior results when applied in surface mass budget modeling.

Frontiers in Earth Science↗

USGS2017

Observationally constrained surface mass balance of Larsen C ice shelf, Antarctica

The surface mass balance (SMB) of the Larsen C ice shelf (LCIS), Antarctica, is poorly constrained due to a dearth of in situ observations. Combining several geophysical techniques, we reconstruct spatial and temporal patterns of SMB over the LCIS. Continuous time series of snow height (2.5–6 years) at five locations allow for multi-year estimates of seasonal and annual SMB over the LCIS. There is high interannual variability in SMB as well as spatial variability: in the north, SMB is 0.40 ± 0.06 to 0.41 ± 0.04 m w.e. year −1 , while farther south, SMB is up to 0.50 ± 0.05 m w.e. year −1 . This difference between north and south is corroborated by winter snow accumulation derived from an airborne radar survey from 2009, which showed an average snow thickness of 0.34 m w.e. north of 66° S, and 0.40 m w.e. south of 68° S. Analysis of ground-penetrating radar from several field campaigns allows for a longer-term perspective of spatial variations in SMB: a particularly strong and coherent reflection horizon below 25–44 m of water-equivalent ice and firn is observed in radargrams collected across the shelf. We propose that this horizon was formed synchronously across the ice shelf. Combining snow height observations, ground and airborne radar, and SMB output from a regional climate model yields a gridded estimate of SMB over the LCIS. It confirms that SMB increases from north to south, overprinted by a gradient of increasing SMB to the west, modulated in the west by föhn-induced sublimation. Previous observations show a strong decrease in firn air content toward the west, which we attribute to spatial patterns of melt, refreezing, and densification rather than SMB.

Larsen C Ice Shelf↗

USGS2017

Fracture propagation and stability of ice shelves governed by ice shelf heterogeneity

Tabular iceberg calving and ice shelf retreat occurs after full‐thickness fractures, known as rifts, propagate across an ice shelf. A quickly evolving rift signals a threat to the stability of Larsen C, the Antarctic Peninsula's largest ice shelf. Here we reveal the influence of ice shelf heterogeneity on the growth of this rift, with implications that challenge existing notions of ice shelf stability. Most of the rift extension has occurred in bursts after overcoming the resistance of suture zones that bind together neighboring glacier inflows. We model the stresses in the ice shelf to determine potential rift trajectories. Calving perturbations to ice flow will likely reach the grounding line. The stability of Larsen C may hinge on a single suture zone that stabilizes numerous upstream rifts. Elevated fracture toughness of suture zones may be the most important property that allows ice shelves to modulate Antarctica's contribution to sea level rise.

Larsen C Ice Shelf↗

USGS2017

Hypsometric control on glacier mass balance sensitivity in Alaska and northwest Canada

Glacier hypsometry provides a first‐order approach for assessing a glacier's response to climate forcings. We couple the Randolph Glacier Inventory to a suite of in situ observations and climate model output to examine potential change for the ∼27,000 glaciers in Alaska and northwest Canada through the end of the 21st century. By 2100, based on Representative Concentration Pathways (RCPs) 4.5–8.5 forcings, summer temperatures are predicted to increase between +2.1 and +4.6°C, while solid precipitation (snow) is predicted to decrease by −6 to −11%, despite a +9 to +21% increase in total precipitation. Snow is predicted to undergo a pronounced decrease in the fall, shifting the start of the accumulation season back by ∼1 month. In response to these forcings, the regional equilibrium line altitude (ELA) may increase by +105 to +225 m by 2100. The mass balance sensitivity to this increase is highly variable, with the most substantive impact for glaciers with either limited elevation ranges (often small (<1 km 2 ) glaciers, which account for 80% of glaciers in the region) or those with top‐heavy geometries, like icefields. For more than 20% of glaciers, future ELAs, given RCP 6.0 forcings, will exceed the maximum elevation of the glacier, resulting in their eventual demise, while for others, accumulation area ratios will decrease by >60%. Our results highlight the first‐order control of hypsometry on individual glacier response to climate change, and the variability that hypsometry introduces to a regional response to a coherent climate perturbation.

Alaska, British Columbia, Yukon↗

USGS2017

Geometry, mass balance and thinning at Eklutna Glacier, Alaska: an altitude-mass-balance feedback with implications for water resources

We analyzed glacier surface elevations (1957, 2010 and 2015) and surface mass-balance measurements (2008–2015) on the 30 km 2 Eklutna Glacier, in the Chugach Mountains of southcentral Alaska. The geodetic mass balances from 1957 to 2010 and 2010 to 2015 are −0.52 ± 0.46 and −0.74 ± 0.10 m w.e. a −1 , respectively. The glaciological mass balance of −0.73 m w.e. a −1 from 2010 to 2015 is indistinguishable from the geodetic value. Even after accounting for loss of firn in the accumulation zone, we found most of the mass loss over both time periods was from a broad, low-slope basin that includes much of the accumulation zone of the main branch. Ice-equivalent surface elevation changes in the basin were −1.0 ± 0.8 m a −1 from 1957 to 2010, and −0.6 ± 0.1 m a −1 from 2010 to 2015, shifting the glacier hypsometry downward and resulting in more negative mass balances: an altitude-mass-balance feedback. Net mass loss from Eklutna Glacier accounts for 7 ± 1% of the average inflow to Eklutna Reservoir, which is entirely used for water and power by Anchorage, Alaska's largest city. If the altitude-mass-balance feedback continues, this ‘deglaciation discharge dividend’ is likely to increase over the short-term before it eventually decreases due to diminishing glacier area.

Alaska↗

USGS2017

Hydrologic and geomorphic changes resulting from episodic glacial lake outburst floods: Rio Colonia, Patagonia, Chile

Glacial lake outburst floods (GLOFs) are a prominent but poorly understood cryospheric hazard in a warming climate. We quantify the hydrologic and geomorphic response to 21 episodic GLOFs that began in April 2008 using multitemporal satellite imagery and field observations. Peak discharge exiting the source lake became progressively muted downstream. At ~40–60 km downstream, where the floods entered and traveled down the main stem Rio Baker, peak discharges were generally < 2000 m 3 s −1 , although these flows were still >1–2 times the peak annual discharge of this system, Chile's largest river by volume. As such, caution must be applied to empirical relationships relating lake volume to peak discharge, as the latter is dependent on where this observation is made along the flood path. The GLOFs and subsequent periods of free drainage resulted in > 40 m of incision, the net removal of ~25 × 10 6 m 3 of sediment from the source lake basin, and a nonsteady channel configuration downstream. These results demonstrate that GLOFs sourced from low‐order tributaries can produce significant floods on major main stem rivers, in addition to significantly altering sediment dynamics.

Cachet‐Colonia‐Baker Valley, Patagonia, Rio Coloni↗

USGS2017

Links between atmosphere, ocean, and cryosphere from two decades of microseism observations on the Antarctic Peninsula

The lack of landmasses, climatological low pressure, and strong circumpolar westerly winds between the latitudes of 50°S to 65°S produce exceptional storm‐driven wave conditions in the Southern Ocean. This combination makes the Antarctic Peninsula one of Earth's most notable regions of high‐amplitude wave activity and thus, ocean‐swell‐driven microseism noise in both the primary (direct wave‐coastal region interactions) and secondary (direct ocean floor forcing due to interacting wave trains) period bands. Microseism observations are examined across 23 years (1993–2015) from Palmer Station (PMSA), on the west coast of the Antarctic Peninsula, and from East Falkland Island (EFI). These records provide a spatially integrative measure of both Southern Ocean wave amplitudes and the interactions between ocean waves and the solid Earth in the presence of sea ice, which can reduce wave coupling with the continental shelf. We utilize a spatiotemporal correlation‐based approach to illuminate how the distribution of sea ice influences seasonal microseism power. We characterize primary and secondary microseism power due to variations in sea ice and find that primary microseism energy is both more sensitive to sea ice and more capable of propagating across ocean basins than secondary microseism energy. During positive phases of the Southern Annular Mode, sea ice is reduced in the Bellingshausen Sea and overall storm activity in the Drake Passage increases, thus strongly increasing microseism power levels.

Antarctic Penninsula↗

USGS2017

Hydrologic impacts of changes in climate and glacier extent in the Gulf of Alaska watershed

High‐resolution regional‐scale hydrologic models were used to quantify the response of late 21st century runoff from the Gulf of Alaska (GOA) watershed to changes in regional climate and glacier extent. NCEP Climate Forecast System Reanalysis data were combined with five Coupled Model Intercomparison Project Phase 5 general circulation models (GCMs) for two representative concentration pathway (RCP) scenarios (4.5 and 8.5) to develop meteorological forcing for the period 2070–2099. A hypsographic model was used to estimate future glacier extent given assumed equilibrium line altitude (ELA) increases of 200 and 400 m. GCM predictions show an increase in annual precipitation of 12% for RCP 4.5 and 21% for RCP 8.5, and an increase in annual temperature of 2.5°C for RCP 4.5 and 4.3°C for RCP 8.5, averaged across the GOA. Scenarios with perturbed climate and glaciers predict annual GOA‐wide runoff to increase by 9% for RCP4.5/ELA200 case and 14% for the RCP8.5/ELA400 case. The glacier runoff decreased by 14% for RCP4.5/ELA200 and by 34% for the RCP8.5/ELA400 case. Intermodel variability in annual runoff was found to be approximately twice the variability in precipitation input. Additionally, there are significant changes in runoff partitioning and increases in snowpack runoff are dominated by increases in rain‐on‐snow events. We present results aggregated across the entire GOA and also for individual watersheds to illustrate the range in hydrologic regime changes and explore the sensitivities of these results by independently perturbing only climate forcings and only glacier cover.

Alaska, British Columbia, Yukon↗

USGS2016

Latest Pleistocene and Holocene glacial events in the Colonia valley, Northern Patagonia Icefield, southern Chile

The Northern Patagonia Icefield (NPI) is the primary glaciated terrain worldwide at its latitude (46.5–47.5°S), and constraining its glacial history provides unique information for reconstructing Southern Hemisphere paleoclimate. The Colonia Glacier is the largest outlet glacier draining the eastern NPI. Ages were determined using dendrochronology, lichenometry, radiocarbon, cosmogenic 10 Be and optically stimulated luminescence. Dated moraines in the Colonia valley defined advances at 13.2 ± 0.95, 11.0 ± 0.47 and 4.96 ± 0.21 ka, with the last being the first constraint on the onset of Neoglaciation for the eastern NPI from a directly dated landform. Dating in the tributary Cachet valley, which contains an ice-dammed lake during periods of Colonia Glacier expansion, defined an advance at ca. 2.95 ± 0.21 ka, periods of advancement at 810 ± 49 cal a BP and 245 ± 13 cal a BP, and retreat during the intervening periods. Recent Colonia Glacier thinning, which began in the late 1800s, opened a lower-elevation outlet channel for Lago Cachet Dos in ca. 1960. Our data provide the most comprehensive set of Latest Pleistocene and Holocene ages for a single NPI outlet glacier and expand previously developed NPI glacial chronologies.

Colonia valley, Northern Patagonia Iceﬁeld↗

USGS2015

End-of-winter snow depth variability on glaciers in Alaska

A quantitative understanding of snow thickness and snow water equivalent (SWE) on glaciers is essential to a wide range of scientific and resource management topics. However, robust SWE estimates are observationally challenging, in part because SWE can vary abruptly over short distances in complex terrain due to interactions between topography and meteorological processes. In spring 2013, we measured snow accumulation on several glaciers around the Gulf of Alaska using both ground- and helicopter-based ground-penetrating radar surveys, complemented by extensive ground truth observations. We found that SWE can be highly variable (40% difference) over short spatial scales (tens to hundreds of meters), especially in the ablation zone where the underlying ice surfaces are typically rough. Elevation provides the dominant basin-scale influence on SWE, with gradients ranging from 115 to 400 mm/100 m. Regionally, total accumulation and the accumulation gradient are strongly controlled by a glacier's distance from the coastal moisture source. Multiple linear regressions, used to calculate distributed SWE fields, show that robust results require adequate sampling of the true distribution of multiple terrain parameters. Final SWE estimates (comparable to winter balances) show reasonable agreement with both the Parameter-elevation Relationships on Independent Slopes Model climate data set (9–36% difference) and the U.S. Geological Survey Alaska Benchmark Glaciers (6–36% difference). All the glaciers in our study exhibit substantial sensitivity to changing snow-rain fractions, regardless of their location in a coastal or continental climate. While process-based SWE projections remain elusive, the collection of ground-penetrating radar (GPR)-derived data sets provides a greatly enhanced perspective on the spatial distribution of SWE and will pave the way for future work that may eventually allow such projections.

Alaska↗

USGS2015

Oceanic and atmospheric forcing of Larsen C Ice-Shelf thinning

The catastrophic collapses of Larsen A and B ice shelves on the eastern Antarctic Peninsula have caused their tributary glaciers to accelerate, contributing to sea-level rise and freshening the Antarctic Bottom Water formed nearby. The surface of Larsen C Ice Shelf (LCIS), the largest ice shelf on the peninsula, is lowering. This could be caused by unbalanced ocean melting (ice loss) or enhanced firn melting and compaction (englacial air loss). Using a novel method to analyse eight radar surveys, this study derives separate estimates of ice and air thickness changes during a 15-year period. The uncertainties are considerable, but the primary estimate is that the surveyed lowering (0.066 ± 0.017 m yr −1 ) is caused by both ice loss (0.28 ± 0.18 m yr −1 ) and firn-air loss (0.037 ± 0.026 m yr −1 ). The ice loss is much larger than the air loss, but both contribute approximately equally to the lowering because the ice is floating. The ice loss could be explained by high basal melting and/or ice divergence, and the air loss by low surface accumulation or high surface melting and/or compaction. The primary estimate therefore requires that at least two forcings caused the surveyed lowering. Mechanisms are discussed by which LCIS stability could be compromised in the future. The most rapid pathways to collapse are offered by the ungrounding of LCIS from Bawden Ice Rise or ice-front retreat past a "compressive arch" in strain rates. Recent evidence suggests that either mechanism could pose an imminent risk.

The Cryosphere↗

USGS2015

Documenting 35 years of land cover change: Lago Cachet Dos drainage, Chile

The U.S. Geological Survey (USGS) Special Applications Science Center is monitoring temporal changes at the Colonia Glacier and Lago Cachet Dos, Northern Patagonia Icefield of southern Chile. This location is one of the newest international sites in the USGS Global Fiducial Program (GFP)—a program which provides systematic monitoring of dynamic and environmentally critical areas with high-resolution imagery (http://gfp.usgs.gov/). In 2008, Lago Cachet Dos began experiencing glacial lake outburst floods (GLOFs) during which the entire pool of water (about 200 million cubic meters) rapidly drains from the lake and flows south-southeast through the Colonia Glacier. These catastrophic events cause massive erosion of valley-fill deposits and consequent upstream expansion of Lago Cachet Dos towards Lago Cachet Uno. Panchromatic and multispectral images for 1979, 2007, and 2014 highlight the dramatic changes that have occurred at this site over a 35-year period. The lake was smallest in 1979, when the Colonia Glacier was at its maximum extent during the study period. Between 1979 and 2007, the glacier shrank causing an increase in the surface area of the lake. The size of the lake increased substantially, from 2.98 square kilometers (km 2 ) in 1979 to 4.41 km 2 in 2014, primarily due to erosion of valley-fill deposits upstream of its northern edge by the 15 GLOFs that occurred between April 2008 and February 2014. Ongoing studies of the Colonia Glacier and Lago Cachet Dos are focused on providing real-time monitoring of Lago Cachet Dos lake levels, understanding the history of advances and retreats of the Colonia Glacier, and determining the physical mechanisms and hazards associated with the GLOFs that come from Lago Cachet Dos.

Lago Cachet Dos↗