Lincoln H. Pitcher

Efficient meltwater drainage through supraglacial streams and rivers on the southwest Greenland ice sheet

Significance Meltwater runoff from the Greenland ice sheet is a key contributor to global sea level rise and is expected to increase in the future, but it has received little observational study. We used satellite and in situ technologies to assess surface drainage conditions on the southwestern ablation surface after an extreme 2012 melting event. We conclude that the ice sheet surface is efficiently drained under optimal conditions, that digital elevation models alone cannot fully describe supraglacial drainage and its connection to subglacial systems, and that predicting outflow from climate models alone, without recognition of subglacial processes, may overestimate true meltwater release from the ice sheet. Thermally incised meltwater channels that flow each summer across melt-prone surfaces of the Greenland ice sheet have received little direct study. We use high-resolution WorldView-1/2 satellite mapping and in situ measurements to characterize supraglacial water storage, drainage pattern, and discharge across 6,812 km2 of southwest Greenland in July 2012, after a record melt event. Efficient surface drainage was routed through 523 high-order stream/river channel networks, all of which terminated in moulins before reaching the ice edge. Low surface water storage (3.6 ± 0.9 cm), negligible impoundment by supraglacial lakes or topographic depressions, and high discharge to moulins (2.54–2.81 cm⋅d−1) indicate that the surface drainage system conveyed its own storage volume every <2 d to the bed. Moulin discharges mapped inside ∼52% of the source ice watershed for Isortoq, a major proglacial river, totaled ∼41–98% of observed proglacial discharge, highlighting the importance of supraglacial river drainage to true outflow from the ice edge. However, Isortoq discharges tended lower than runoff simulations from the Modèle Atmosphérique Régional (MAR) regional climate model (0.056–0.112 km3⋅d−1 vs. ∼0.103 km3⋅d−1), and when integrated over the melt season, totaled just 37–75% of MAR, suggesting nontrivial subglacial water storage even in this melt-prone region of the ice sheet. We conclude that (i) the interior surface of the ice sheet can be efficiently drained under optimal conditions, (ii) that digital elevation models alone cannot fully describe supraglacial drainage and its connection to subglacial systems, and (iii) that predicting outflow from climate models alone, without recognition of subglacial processes, may overestimate true meltwater export from the ice sheet to the ocean.

Derivation of High Spatial Resolution Albedo from UAV Digital Imagery: Application over the Greenland Ice Sheet

Measurements of albedo are a prerequisite for modelling surface melt across the Earths cryosphere, yet available satellite products are limited in spatial and/or temporal resolution. Here, we present a practical methodology to obtain centimetre resolution albedo products with accuracies of 5% using consumer-grade digital camera and unmanned aerial vehicle (UAV) technologies. Our method comprises a workflow for processing, correcting and calibrating raw digital images using a white reference target, and upward and downward shortwave radiation measurements from broadband silicon pyranometers. We demonstrate the method with a set of UAV sorties over the western, K-sector of the Greenland Ice Sheet. The resulting albedo product, UAV10A1, covers 280 km2, at a resolution of 20 cm per pixel and has a root-mean-square difference of 3.7% compared to MOD10A1 and 4.9% compared to ground-based broadband pyranometer measurements. By continuously measuring downward solar irradiance, the technique overcomes previous limitations due to variable illumination conditions during and between surveys over glaciated terrain. The current miniaturization of multispectral sensors and incorporation of upward facing radiation sensors on UAV packages means that this technique will likely become increasingly attractive in field studies and used in a wide range of applications for high temporal and spatial resolution surface mapping of debris, dust, cryoconite and bioalbedo and for directly constraining surface energy balance models.

Fluvial morphometry of supraglacial river networks on the southwest Greenland Ice Sheet

Extensive, complex supraglacial river networks form on the southwest Greenland ice sheet (GrIS) surface each melt season. These networks are the dominant pathways for surface meltwater transport on this part of the ice sheet, but their fluvial morphometry has received little study. This paper utilizes high-resolution (2 m) WorldView-1/2 images, digital elevation models, and GIS tools to present a detailed morphometric characterization (river number, river length, Strahler stream order, width, depth, bifurcation ratio, braiding index, drainage density, slope, and relief ratio) for 523 GrIS supraglacial river networks. A new algorithm is presented to determine Strahler stream order in supraglacial environments. Results show that (1) Supraglacial river networks are broadly similar to terrestrial landscapes in that they follow Horton’s laws (river number, mean river length, and slope versus stream order), widen downstream, and have comparable mean bifurcation ratios (3.7 ± 1.9) and braiding indices; (2) unlike terrestrial systems, supraglacial drainage densities (0.90–4.75 km/km2) have no correlation with elevation relief, but instead display a weakly inverse correlation with ice surface elevation; (3) both well-developed (e.g., fifth-order) and discrete (e.g., first-order) supraglacial river networks form on the ice sheet, with the latter associated with short flow distances upstream of a terminal moulin; (4) mean river flow widths increase substantially, but flow depths only modestly, with increasing stream order. Viewed collectively, the 523 supraglacial river networks studied here display fluvial morphometries both similar and dissimilar to terrestrial systems, with moulin capture an important physical process driving the latter.

Direct measurements of meltwater runoff on the Greenland ice sheet surface

Significance Meltwater runoff is an important hydrological process operating on the Greenland ice sheet surface that is rarely studied directly. By combining satellite and drone remote sensing with continuous field measurements of discharge in a large supraglacial river, we obtained 72 h of runoff observations suitable for comparison with climate model predictions. The field observations quantify how a large, fluvial supraglacial catchment attenuates the magnitude and timing of runoff delivered to its terminal moulin and hence the bed. The data are used to calibrate classical fluvial hydrology equations to improve meltwater runoff models and to demonstrate that broad-scale surface water drainage patterns that form on the ice surface powerfully alter the timing, magnitude, and locations of meltwater penetrating into the ice sheet. Meltwater runoff from the Greenland ice sheet surface influences surface mass balance (SMB), ice dynamics, and global sea level rise, but is estimated with climate models and thus difficult to validate. We present a way to measure ice surface runoff directly, from hourly in situ supraglacial river discharge measurements and simultaneous high-resolution satellite/drone remote sensing of upstream fluvial catchment area. A first 72-h trial for a 63.1-km2 moulin-terminating internally drained catchment (IDC) on Greenland’s midelevation (1,207–1,381 m above sea level) ablation zone is compared with melt and runoff simulations from HIRHAM5, MAR3.6, RACMO2.3, MERRA-2, and SEB climate/SMB models. Current models cannot reproduce peak discharges or timing of runoff entering moulins but are improved using synthetic unit hydrograph (SUH) theory. Retroactive SUH applications to two older field studies reproduce their findings, signifying that remotely sensed IDC area, shape, and supraglacial river length are useful for predicting delays in peak runoff delivery to moulins. Applying SUH to HIRHAM5, MAR3.6, and RACMO2.3 gridded melt products for 799 surrounding IDCs suggests their terminal moulins receive lower peak discharges, less diurnal variability, and asynchronous runoff timing relative to climate/SMB model output alone. Conversely, large IDCs produce high moulin discharges, even at high elevations where melt rates are low. During this particular field experiment, models overestimated runoff by +21 to +58%, linked to overestimated surface ablation and possible meltwater retention in bare, porous, low-density ice. Direct measurements of ice surface runoff will improve climate/SMB models, and incorporating remotely sensed IDCs will aid coupling of SMB with ice dynamics and subglacial systems.

AirSWOT measurements of river water surface elevation and slope: Tanana River, AK

Fluctuations in water surface elevation (WSE) along rivers have important implications for water resources, flood hazards, and biogeochemical cycling. However, current in situ and remote sensing methods exhibit key limitations in characterizing spatiotemporal hydraulics of many of the worlds river systems. Here we analyze new measurements of river WSE and slope from AirSWOT, an airborne analogue to the Surface Water and Ocean Topography (SWOT) mission aimed at addressing limitations in current remotely sensed observations of surface water. To evaluate its capabilities, we compare AirSWOT WSEs and slopes to in situ measurements along the Tanana River, Alaska. Root-mean-square error is 9.0 cm for WSEs averaged over 1 km2 areas and 1.0 cm/km for slopes along 10 km reaches. Results indicate that AirSWOT can accurately reproduce the spatial variations in slope critical for characterizing reach-scale hydraulics. AirSWOTs high-precision measurements are valuable for hydrologic analysis, flood modeling studies, and for validating future SWOT measurements.

CryoSheds: a GIS modeling framework for delineating land-ice watersheds for the Greenland Ice Sheet

Choice of watershed delineation technique is an important source of uncertainty for cryo-hydrologic studies of the Greenland Ice Sheet (GrIS), with different methods yielding different watersheds for a common pour point. First, this paper explores this uncertainty for the Akuliarusiarsuup Kuua River Northern Tributary, Western Greenland. Next, a standardized, semi-automated modeling framework for generating land-ice watersheds for GrIS land-terminating ice (henceforth referred to as CryoSheds) using geographic information systems (GIS) hydrologic modeling tools is presented. The framework uses ArcGIS and the ArcPy geoprocessing library to delineate two types of land-ice watersheds, namely those defined by: (1) a hydraulic pressure potential with varying water to ice overburden pressure ratios (k-value), which determines theoretical flow paths from the hydrostatic equation, using surface and bedrock digital elevation models (DEMs) and (2) a surface topography DEM alone. Lastly, a demonstration of the CryoSheds method is presented for seven remotely sensed proglacial pour points along the Aussivigssuit River (AR), Western Greenland, and its largest tributaries. GrIS meltwater runoff from these seven nested land-ice watersheds is estimated using Modele Atmospherique Regional (MAR) v.3.2 and runoff uncertainties due to watershed delineation parameter selection is estimated.

Threats of future climate change and land use to vulnerable tree species native to Southern California

Author(s): Riordan, EC; Gillespie, TW; Pitcher, LH; Pincetl, SS; Jenerette, GD; Pataki, DE | Abstract: Climate and land-use changes are expected to drive high rates of environmental change and biodiversity loss in Mediterranean ecosystems this century. This paper compares the relative future impacts of land use and climate change on two vulnerable tree species native to Southern California (Juglans californica and Quercus engelmannii) using species distribution models. Under the Intergovernmental Panel for Climate Changes A1B future scenario, high levels of both projected land use and climate change could drive considerable habitat losses on these two already heavily-impacted tree species. Under scenarios of no dispersal, projected climate change poses a greater habitat loss threat relative to projected land use for both species. Assuming unlimited dispersal, climate-driven habitat gains could offset some of the losses due to both drivers, especially in J. californica which could experience net habitat gains under combined impacts of both climate change and land use. Quercus engelmannii, in contrast, could experience net habitat losses under combined impacts, even under best-case unlimited dispersal scenarios. Similarly, projected losses and gains in protected habitat are highly sensitive to dispersal scenario, with anywhere from g 60% loss in protected habitat (no dispersal) to g 170% gain in protected habitat (unlimited dispersal). The findings underscore the importance of dispersal in moderating future habitat loss for vulnerable species. Copyright

Correction to Supporting Information for Smith et al., Direct measurements of meltwater runoff on the Greenland ice sheet surface

EARTH, ATMOSPHERIC, AND PLANETARY SCIENCES Correction to Supporting Information for “Direct measurements of meltwater runoff on the Greenland ice sheet surface,” by Laurence C. Smith, Kang Yang, Lincoln H Pitcher, Brandon T. Overstreet, Vena W. Chu, Åsa K. Rennermalm, Jonathan C. Ryan, Matthew G. Cooper, Colin J. Gleason, Marco Tedesco, Jeyavinoth Jeyaratnam, Dirk van As, Michiel R. van den Broeke, Willem Jan van de Berg, Brice Noël, Peter L. Langen, Richard I. Cullather, Bin Zhao, Michael J. Willis, Alun Hubbard, Jason E. Box, Brittany A. Jenner, and Alberto E. Behar, which was first published December 5, 2017; 10.1073/pnas.1707743114 (Proc Natl Acad Sci USA 114:E10622–E10631). The authors note that in the SI Appendix, page 23, line 570, “tp = Ct(Lc) 0.3 ” should instead appear as “tp = Ct(LLc) .” Additionally, on page 24 of the SI Appendix, line 613, Eq. 1 should instead appear as: