Richard Forster | Search USGS

USGS2022

Hydrologic modeling of a perennial firn aquifer in southeast Greenland

A conceptual model, based on field observations and assumed physics of a perennial firn aquifer near Helheim Glacier (southeast Greenland), is evaluated via steady-state 2-D simulation of liquid water flow and energy transport with phase change. The simulation approach allows natural representation of flow and energy advection and conduction that occur in vertical meltwater recharge through the unsaturated zone and in lateral flow within the saturated aquifer. Agreement between measured and simulated aquifer geometry, temperature, and recharge and discharge rates confirms that the conceptual field-data-based description of the aquifer is consistent with the primary physical processes of groundwater flow, energy transport and phase change. Factors that are found to control simulated aquifer configuration include surface temperature, meltwater recharge rate, residual total-water saturation and capillary fringe thickness. Simulation analyses indicate that the size of perennial firn aquifers depends primarily on recharge rates from surface snowmelt. Results also imply that the recent aquifer expansion, likely due to a warming climate, may eventually produce lakes on the ice-sheet surface that would affect the surface energy balance.

Journal of Glaciology↗

USGS2020

Hydrologic properties of a highly permeable firn aquifer in the Wilkins Ice Shelf, Antarctica

We present measurements of the density, hydraulic conductivity, and specific discharge of a widespread firn aquifer in Antarctica, within the Wilkins Ice Shelf. At the field site, the aquifer is 16.2 m thick, starting at 13.4 m from the snow surface and transitioning from water‐saturated firn to ice at 29.6 m. Hydraulic conductivity derived from slug tests show a geometric mean value of 1.4 ± 1.2 × 10 −4 m s −1 , equivalent to permeability of 2.6 ± 2.2 × 10 −11 m 2 . A borehole dilution test indicates an average specific discharge value of 1.9 ± 2.8 × 10 −6 m s −1 . Ground‐penetrating radar profiles and a groundwater flow model show the aquifer is draining laterally into a large nearby rift. Our findings indicate that the firn aquifer in the vicinity of the field site is likely not in a steady state and its presence likely contributed to past ice shelf instability.

Geophysical Research Letters↗

USGS2014

Prevalence of pure versus mixed snow cover pixels across spatial resolutions in alpine environments: implications for binary and fractional remote sensing approaches

Remote sensing of snow-covered area (SCA) can be binary (indicating the presence/absence of snow cover at each pixel) or fractional (indicating the fraction of each pixel covered by snow). Fractional SCA mapping provides more information than binary SCA, but is more difficult to implement and may not be feasible with all types of remote sensing data. The utility of fractional SCA mapping relative to binary SCA mapping varies with the intended application as well as by spatial resolution, temporal resolution and period of interest, and climate. We quantified the frequency of occurrence of partially snow-covered (mixed) pixels at spatial resolutions between 1 m and 500 m over five dates at two study areas in the western U.S., using 0.5 m binary SCA maps derived from high spatial resolution imagery aggregated to fractional SCA at coarser spatial resolutions. In addition, we used in situ monitoring to estimate the frequency of partially snow-covered conditions for the period September 2013–August 2014 at 10 60-m grid cell footprints at two study areas with continental snow climates. Results from the image analysis indicate that at 40 m, slightly above the nominal spatial resolution of Landsat, mixed pixels accounted for 25%–93% of total pixels, while at 500 m, the nominal spatial resolution of MODIS bands used for snow cover mapping, mixed pixels accounted for 67%–100% of total pixels. Mixed pixels occurred more commonly at the continental snow climate site than at the maritime snow climate site. The in situ data indicate that some snow cover was present between 186 and 303 days, and partial snow cover conditions occurred on 10%–98% of days with snow cover. Four sites remained partially snow-free throughout most of the winter and spring, while six sites were entirely snow covered throughout most or all of the winter and spring. Within 60 m grid cells, the late spring/summer transition from snow-covered to snow-free conditions lasted 17–56 days and averaged 37 days. Our results suggest that mixed snow-covered snow-free pixels are common at the spatial resolutions imaged by both the Landsat and MODIS sensors. This highlights the additional information available from fractional SCA products and suggests fractional SCA can provide a major advantage for hydrological and climatological monitoring and modeling, particularly when accurate representation of the spatial distribution of snow cover is critical.

Remote Sensing↗